On a generalized diffusion equation arising in petroleum engineering

Transcription

1 Al-Homidan et al. Advances in Difference Equations 213, 213:349 R E S E A R C H Open Access On a generalized diffusion equation arising in petroleum engineering Suliman Al-Homidan 1, Ryad A Ghanam 2 and Nasser-eddine Tatar 1* * Correspondence: tatarn@kfupm.edu.sa 1 Department of Mathematics and Statistics, King Fahd University of Petroleum and Minerals, Dhahran, 31261, Saudi Arabia Full list of author information is available at the end of the article Abstract In this work we investigate a model which describes diffusion in petroleum engineering. The original model, which already generalizes the standard one usual diffusion equation) to a non-local model taking into account the memory effects, is here further extended to cope with many other different possible situations. Namely, we consider the Hilfer fractional derivative which by its nature interpolates the Riemann-iouville fractional derivative and the Caputo fractional derivative the one that has been studied previously). At the same time, this kind of derivative provides us with a whole range of other types of fractional derivatives. We treat both the Neumann boundary conditions case and the Dirichlet boundary conditions case and find explicit solutions. In addition to that, we also discuss the case of an infinite reservoir. 1 Introduction In order to modernize the public water service of the town of Dijon, France, Henry Darcy made several experiments and wrote an interesting document which soon formed the basis of the theory of fluid conduction. His work was about the download flow of water through filter sands. He established a diffusion) equation, well known nowadays as Darcy s law, which plays the same role as Fourier s law in heat conduction theory and Ohm s law in the electricity conduction theory see [1, 2). Of course, this law has its limitations. It is restricted to the situations where the flow through the pores can be modeled as Stokes flows. In particular, for large values of the Reynolds number or at high values of flow rates) Darcy s law is not valid anymore. It losses its accuracy. Apart from that, Darcy s law has proved its usefulness in many applications such as the recovery of fuel from underground oil reservoirs. It defines the rock permeability which controls the directional movement of the flow rate). Many researchers have shown interest in this law and have worked extensively on extending and generalizing it to more complicated situations and in different contexts. For instance, they observed that this flow is linear : Darcy himself assumed that the flow is weak and as such the pressure drop is linearly related to the flow discharge rate. This motivated them to extend this law to the nonlinear case. They also considered the case when the inertial forces and/or the deformation of the solid) are not negligible as compared with those arising from viscosity. In particular, they generalized the law from Stokes to non-stokes flows in porous media. In fact, it has been found that more general physical formulations may be obtained by including new parameters and forces in the equation 213 Al-Homidan et al.; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons Attribution icense which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

2 Al-Homidan et al. Advances in Difference Equations 213, 213:349 Page 2 of 14 itself. Moreover, we may mention here the transition from a single phase to multiphase and from fixed to growing porous media. In [3 Caputo modified the standard diffusion equation p t = k 2 p x 2 by introducing a fractional derivative in the sense of Caputo) in Darcy s law [1 to take into account the memory effects that may prevail in the fluid. The equation became p t = k 2 x 2 c D α +p), <α <1, 1) where c D α + denotes the Caputo fractional derivative defined by C D α +f ) t)= 1 Ɣ1 α) t s) α f s) ds, t >,<α <1. Similar equations have been derived in other situations like the one in [4which describes anomalous subdiffusion particles and p is the probability density function) and in [5 when studying the Stokes first problem see also [6). A natural question which arises is: what about the commonly used Riemann-iouville fractional derivative, R D α +f ) 1 d t)= Ɣ1 α) dt t s) α f s) ds, t >,<α <1, or other types of fractional derivatives? It is well known by now that the Caputo fractional derivative is often preferred over the Riemann-iouville fractional derivative at least for two main reasons. The Caputo derivative allowsustouse the usualinitial data like px,)=f x)), whereas the Riemann- iouville derivative requires the prior knowledge of D α 1 px,)=gx) when the order of the equation is α say < α < 1). In addition to that, the Caputo derivative of a constant is zero which is not the case for the Riemann-iouville derivative. In the present note, we consider a generalized fractional derivative which encompasses both fractional derivatives as special cases and provides a whole range of other types of fractional derivatives in between. This derivative was introduced by Hilfer in [7 9 α,β D + f ) t)= I β1 α) d ) dt I1 α)1 β) f t), < α <1, β 1, where I α +f ) t)= 1 Ɣα) t s) α 1 f s) ds, t >,α >, and therefore we name it after him. This kind of fractional derivative has already proved its usefulness see [7 9), and we are witnessing a growing interest in it. The parameter β when equal to zero gives the Riemann-iouville derivative and when it takes the value one gives the Caputo derivative. For < β < 1, we obtain many fractional derivatives interpolating these two types of well-known derivatives.

3 Al-Homidan et al. Advances in Difference Equations 213, 213:349 Page 3 of 14 Using the aplace-fourier transform, we find the explicit solution for our equation under appropriate boundary conditions. We prove here that, in our case, all these derivatives when considered in diffusion equation 1) lead to similar solutions. Consequently, there is no major benefit in treating separately these derivatives unless there is a need to unify these treatments. For more on fractional calculus and more on interesting fractional differential equations and treatments, we refer the reader to [1 15. In the next section, we present some material needed in our arguments. In Section 3 we consider the modified diffusion equation with Neumann boundary conditions and indicatehow tofind anexplicitsolution. Section 4 contains the Dirichlet boundary conditions case.finally,in Section 5 we treat the infinite reservoir case. Our results are illustrated by some graphs. 2 Preliminaries In this section we present some definitions and results which will be needed later in our arguments see [12 15formore). Definition 1 The Riemann-iouville fractional integral of order α of f is defined by I α +f ) t)= 1 Ɣα) when the right-hand side exists. t s) α 1 f s) ds, t >,α > Definition 2 The Riemann-iouville fractional derivative of order α of f is defined by R D α +f ) 1 d t)= Ɣ1 α) dt when the right-hand side exists. Note that R D α +f ) t)= d I 1 α f ) t). dt t s) α f s) ds, t >,<α <1 Definition 3 The Caputo fractional derivative of order α of f is defined by C D α +f ) t)= 1 Ɣ1 α) t s) α f s) ds, t >,<α <1 the prime here is for the derivative) when the right-hand side exists. Note that C D α +f ) t)= I 1 α ddt ) f t). The relationship between these two types of derivatives is given by the following theorem. Theorem 1 [12 We have R D α +f ) t)= C D α +f ) t)+ t α Ɣ1 α) f +), t >,<α <1.

4 Al-Homidan et al. Advances in Difference Equations 213, 213:349 Page 4 of 14 Definition 4 The Hilfer fractional derivative of f of order α and type β is defined by α,β D + f ) t)= I β1 α) d ) dt I1 α)1 β) f t), < α <1, β 1 whenever the right-hand side exists. Note that when β =, D α, + f ) t)= d dt I 1 α f ) t), which is the Riemann-iouville fractional derivative see Definition 2), and when β =1, D α,1 + f ) t)= I 1 α ddt ) f t), which is the Caputo fractional derivative see Definition 3). For < α < 1, the aplace transforms of these derivatives are given by [ R D α +f ) t) s)=s α [ f t) s) I 1 α + f ) +), 2) [C D α +f ) t) s)=s α [ f t) s) s α 1 f +), 3) [ D α,β + f ) t) s)=s α [ f t) s) s βα 1) I 1 α)1 β) + f ) +), β 1. 4) It is clear that the differences in these aplace transforms are in the initial data f + ), I 1 α + f )+ )andi 1 α)1 β) + f )+ ) this last one is natural initial data for the Hilfer derivative). The natural space for the Hilfer fractional derivative is where C α,β 1 γ [a, b={ f C 1 γ [a, b, D α,β a + f C 1 γ [a, b }, C 1 γ a, b= { g :a, b R :t a) 1 γ gt) C[a, b }. However, here, since our problem is of order one, solutions must be much smoother. 3 Theproblem We shall investigate the following linear generalized fractional diffusion problem: p x, t)= k t + px, t)), x, t >, px,)=p i, x, q = [ ka μ x Dα,β + p), t >, x= q = [ ka μ x Dα,β + p), t >. x= 2 μφc t x 2 D α,β 5) The different parameters μ, φ, C t, k and A are positive constants which account for the viscosity, porosity, total compressibility, permeability and the area, respectively. This model describes the flow of oil in a finite reservoir. A closely related and interesting work is in [16 but for the normal form of the diffusion.

5 Al-Homidan et al. Advances in Difference Equations 213, 213:349 Page 5 of 14 This problem, for the Caputo fractional derivative, was derived and studied in [17. Here we consider the more general Hilfer) fractional derivative D α,β + which covers both the Riemann-iouville derivative and the Caputo derivative in addition to many others for β in the interval, 1)). Following the argument in [17, we introduce the dimensionless variables: P D = p i p p i, t D = ηt and X 2 D = x,whereη = k μφc t. The next two lemmas will, in particular, show that we obtain the same explicit solution using the aplace transform as in the case of Caputo derivative see [17). emma 1 If t = Ct D, then D μ,ν t px, t)=c μ D μ,ν t D px, t), < μ <1, ν 1, x, t >. Proof From Definition 4,wehave,for<μ <1, ν 1, D μ,ν t px, t) = I ν1 μ) DI 1 μ)1 ν) px, t) t τ) ν1 μ) 1 [ 1 = Ɣ[ν1 μ) τ Ɣ[1 μ)1 ν) et t = Ct D,weseethat D μ,ν t px, t) =A CtD Put τ = Cξ, then clearly τ τ s) 1 μ)1 ν) 1 px, s) ds dτ. Ct D t) ν1 μ) 1 [ τ τ s) 1 μ)1 ν) 1 px, s) ds dτ. τ D μ,ν t px, t) D = AC ν1 μ) 1 t D ξ) ν1 μ) 1 Cξ) = AC ν1 μ) 1 D Put also s = Cu,itappearsthat [ Cξ Cξ s) 1 μ)1 ν) 1 px, s) ds Cdξ t D ξ) ν1 μ) 1 [ Cξ Cξ s) 1 μ)1 ν) 1 px, s) ds ξ dξ. D μ,ν t px, t)=ac ν1 μ) 1 C 1 μ)1 ν) 1 D t D ξ) ν1 μ) 1 [ ξ ξ D = C μ A t D ξ) ν1 μ) 1 [ ξ ξ = C μ I ν1 μ) DI 1 μ)1 ν) px, t D )=C μ D μ,ν t D px, t). ξ u) 1 μ)1 ν) 1 px, Cu)Cdu dξ ξ u) 1 μ)1 ν) 1 px, u) du dξ The proof is complete. emma 2 Assume that f t) is continuous on [,A for some A >,then lim t + Iα +f t)=, α >.

6 Al-Homidan et al. Advances in Difference Equations 213, 213:349 Page 6 of 14 Proof et t A,then t s) α 1 f s) ds t s) α 1 f s) ds M t s) α 1 ds M α t s) t = M α α tα, where M is a bound for f t)on[,a. In view of our problem 5), the solution px, t) must be at least absolutely continuous. Its continuity nearby zero implies by emma 2 that I+ 1 αf + )=I 1 α)1 β) + f + ) =. Therefore all three derivatives have equal aplace transforms see 2)-4)) if in addition f + )=). This fact together with emma 1 shows that the argument in [17 is valid without any changes. Observe that the change of variables has lead to px D,)=,whichbyTheorem1 implies that both the Riemann-iouville derivative and the Caputo derivative are equal. Another way to see this fact is from the next two lemmas. Definition 5 Afunctionf t) 1 a, b)thespaceofsummablefunctions)issaidtohavea summable fractional derivative D α +f if I1 α + f AC[a, b) the space of absolutely continuous functions). emma 3 [18, p.48 If f t) has a summable fractional derivative R D β +fwith<β <1, then we have I+ α R D β +f t)=iα β + f f + ) t) I1 β t α 1. Ɣα) emma 4 [12,p.95 If f t) C[a, b, < α <1,then C D α +Iα f )t)=f t), t [a, b. Now, Definition 3, Definition 4 and emma 3 imply that D α,β f t)=i β1 α) DI 1 α)1 β) f t)= C D 1 β1 α) I 1 α)1 β) f t) = C D 1 β1 α) I 1 β1 α) α f t) = C D 1 β1 α) [I 1 β1 α)r D α f t)+ I 1 α f + ) Ɣ1 β1 α)) t β1 α). 6) This holds if f t) has a summable fractional derivative D α +f t) thatisi1 α + f AC[, T), which is the case when f itself is absolutely continuous on [, T forsomet >.Inour case, px, t) is absolutely continuous. emma 2 ensures that I 1 α + f ) =. Therefore, D α,β + f t)=c D 1 β1 α) I 1 β1 α)r D α + f t)=r D α + f t) by emma 4 if R D α + f t) is continuous. This is the case if f + )=andf is continuously differentiable because I 1 α + maps C[, T intoc[, T italsomapsc γ [, T intoc[, T, see [19).

7 Al-Homidan et al. Advances in Difference Equations 213, 213:349 Page 7 of 14 4 Dirichlet boundary conditions Consider the problem p t = k 2 D α,β p)+f x, t), x, t >, x 2 px,)=p i, x, p, t)=gt), t >, p, t)=ht), t >, 7) where gt), ht) are continuously differentiable and f x, t) is a forcing term. Putting px, t) =Px, t) +ux, t), with ux, t)=gt)+ x [ ht) gt), we see that implies Clearly, [ 2 Px, t)+ux, t) = k D α,β Px, t)+d α,β ux, t) ) + f x, t) t x 2 P 2 x, t)=k D α,β Px, t) ) + kd α,β 2 u t x 2 x u + f x, t). 8) 2 t 2 u x = and u 2 t = g t)+ x [ h t) g t). 9) Moreover, p i = px,)=px,)+ux,)=px,)+g) + x [ h) g). So Px,)=p i ux,)=p i g) x [ h) g). 1) In addition to that, we have P,t) =P, t) =. 11) Therefore from 8)-11), Px, t) satisfies P t = k 2 D α,β P)+f x, t), x, t >, x 2 Px,)=P x), x, P, t)=p, t)=, t >, 12) where f x, t)=f x, t) u t = f x, t) g t) x [ h t) g t)

8 Al-Homidan et al. Advances in Difference Equations 213, 213:349 Page 8 of 14 and P x)=p i g) x [ h) g). Definition 6 The finite Fourier sine transform of f x), <x <,is defined by F s n)= and its inverse is f x) sin nπx dx, f x)= 2 n=1 F s n) sin nπx. Theorem 2 The following holds: ) v F s = nπ x F cv), where F c is the finite Fourier cosine transform and 2 ) v F s = nπ ) v x 2 F c = n2 π 2 F x 2 s v)+ nπ [ v, t) v, t) cos nπ. Applying first the aplace transform to the equation in 12), we find spx, s) Px,)=k 2 x 2 [ s α Px, s) s βα 1) I 1 α)1 β) P ) +) + f x, s), 13) where Pand f denote the aplace transforms of P and f, respectively. In view of emma 2,equation13)reducesto spx, s)=ks α 2 P x 2 x, s)+ f x, s)+p x). 14) Next, we apply the Fourier sine transform to 14)toget s Pn, s)= n2 π 2 k 2 s α Pn, s)+f s [ f x, s)+p x), where the hat is for the Fourier sine transform. Therefore, Pn, s)= s α F s [ f x, s)+p x) s 1 α + n2 π 2 k 2. 15) From the relation F s a + bx)= [ a 1) n a + b) nπ

9 Al-Homidan et al. Advances in Difference Equations 213, 213:349 Page 9 of 14 it follows that F s P x) ) = [ pi g) 1) n p i h) ). 16) nπ Furthermore, it is easy to see that F s [ f x, s) [ = F s f x, s) ḡ s) x h s) ḡ s) ) = ˆ f [ n, s) F s ḡ s)+ x h s) ḡ s) ) = ˆ f n, s) [ḡ s) 1) n h s). 17) nπ Using the inverse of the Fourier sine transform, we find Px, s)= 2 n=1 Pn, s)sin nπx, and then, inverting the aplace transform, we obtain Px, t)= 2 n=1 ) Pn, 1 nπx s)sin. 18) In view of 15)we have 1[ Pn, [ s) = 1 s α F s [ f [ x, s) s + 1 α F s [P x) s 1 α + n2 π 2 k s 2 1 α + n2 π 2 k 2 Now, appealing to the formula [ s = 1 α [ ˆ f n, s) s 1 α s 1 α + n2 π 2 k nπ 2 ḡ s) 1) n h s)) s 1 α + n2 π 2 K 2 [ + 1 pi g) 1) n p i h) )) s α nπ s 1 α + n2 π 2 k 2. 19) [ t β 1 E α,β ±at α ) = sα β s α a, we see that 1 [ Further s α s 1 α + n2 π 2 k 2 = E 1 α,1 n2 π 2 ) k t 1 α. 2) 2 [ s 1 α ˆ f n, s) t = ˆf n, t s)e s 1 α + n2 π 2 k 1 α,1 n2 π 2 ) k s 1 α ds 21) 2 2

10 Al-Homidan et al. Advances in Difference Equations 213, 213:349 Page 1 of 14 and [ s 1 α ḡ s) 1) n h s)) s 1 α + n2 π 2 k 2 = [ g t s) 1) n h t s) E 1 α,1 n2 π 2 k 2 s 1 α ) ds. 22) Taking into account 2)-22) in19), we obtain 1[ Pn, s) = Finally, relation 18)leadsto {ˆf n, t s) [ g t s) 1) n h t s) } E 1 α,1 n2 π 2 ) k s 1 α ds 2 + [ pi g) 1) n p i h) ) E 1 α,1 n2 π 2 k t 1 α nπ 2 ). 23) Px, t)= 2 {ˆf n, t s) [ g t s) 1) n h t s) } n=1 E 1 α,1 n2 π 2 ) k s 1 α ds sin nπx [ pi g) 1) n p i h) ) nπ n=1 E 1 α,1 n2 π 2 ) k t 1 α sin nπx 2, 24) and px, t) =Px, t) +ux, t) gives the explicit solution of problem 7). 5 Explicit solution in an infinite reservoir In this section we consider a similar problem but in an infinite reservoir, which is the case in many practical situations. p 2 x, t)=k D α,β px, t)) + f x, t), < x <, t >, t x 2 px,)=gx), < x <, p±, t)=, t >. 25) Figure 1 The graph of Px, t)whenα =.2.

11 Al-Homidan et al. Advances in Difference Equations 213, 213:349 Page 11 of 14 Figure 2 The graph of Px, t)whenα =.5. Figure 3 The graph of Px, t)whenα =.8. Figure 4 The graph of Px, t)whenα =.95. The aplace transform with respect to the time variable applied to the equation in 25) gives s px, s) px,)=k 2 x 2 [ sα px, s) s βα 1) I 1 α)1 β) p ) x, +) + f x, s). As px,) = gx) andi 1 α)1 β) p)x,) =, by the continuity of px, t) nearbyt =see emma 2), we get s px, s)=ks α 2 p x 2 x, s)+ f x, s)+gx). 26)

12 Al-Homidan et al. Advances in Difference Equations 213, 213:349 Page 12 of 14 Figure 5 The pressure curve at various distances from the well. Figure 6 The pressure curve at various times. Next, we apply the Fourier transform ˆf y)= 1 2π + f x)eixy dx)tothislastequation26) to arrive at sˆ py, s)= ks α y 2 ˆ py, s)+ˆ f y, s)+ĝy). Therefore [ s + ks α y 2ˆ py, s)=ˆ f y, s)+ĝy),

13 Al-Homidan et al. Advances in Difference Equations 213, 213:349 Page 13 of 14 and hence ˆ py, s) = ˆ f y, s) s + ks α y + ĝy) 2 s + ks α y. 27) 2 Now, the inverse aplace transform of 27)yields ˆpy, t) = E 1 α,1 ky 2 t s) 1 α)ˆf y, s) ds + ĝy)e 1 α,1 ky 2 t 1 α), and the inverse Fourier transform gives px, t)= + E 1 α,1 ky 2 t s) 1 α)ˆf ) y, s) ds e iyx dy ĝy)e 1 α,1 ky 2 t 1 α) e iyx dy. Finally we mention that this type of fractional derivative was also studied in [2 and refer the reader to [21, 22 for other interesting problems and methods. The graphs in Figures 1-6 correspond to different orders and the following data: f x, t)=, Px,) = 3,, P, t)=3,2,p1,6, t) = 2,, k =7.5937e+5. Competing interests The authors declare that they have no competing interests. Authors contributions SA and NT solved the problems. RAG participated in establishing the explicit solution and carried out the numerical treatment. All authors read and approved the final manuscript. Author details 1 Department of Mathematics and Statistics, King Fahd University of Petroleum and Minerals, Dhahran, 31261, Saudi Arabia. 2 Department of Mathematics, University of Pittsburgh at Greensburg, Greensburg, PA 1561, USA. Acknowledgements The authors are grateful for the financial support and the facilities provided by King Fahd University of Petroleum and Minerals and King Abdulaziz City for Science and Technology KACST) NSTIP research grant No. 11-OI Received: 5 June 213 Accepted: 7 November 213 Published: 3 Dec 213 References 1. Darcy, H: es Fontaines Publiques de la Ville de Dijon. Dalmont, Paris 1856) 2. King, M: Darcy s law and the field equations of the flow of underground fluids. Bull. Int. Assoc. Sci. Hydrol. 21), ) 3. Caputo, M: Diffusion of fluids in porous media with memory. Geothermics 28, ) 4. Metzler, R, Klafter, J: The random walk s guide to anomalous diffusion: a fractional dynamics approach. Phys. Rep. 339, ) 5. Salim, TO, El-Kahlout, A: Solution of fractional order Rayleigh-Stokes equations. Adv. Theor. Appl. Mech. 15), ) 6. anglands, TAM: Solution of a modified fractional diffusion equation. Physica A 367, ) 7. Hilfer, R: Fractional time evolution. In: Hilfer, R ed.) Applications of Fractional Calculus in Physics, p. 87. World Scientific, Singapore 2) 8. Hilfer, R: Fractional calculus and regular variation in thermodynamics. In: Hilfer, R ed.) Applications of Fractional Calculus in Physics, p World Scientific, Singapore 2) 9. Hilfer, R: Experimental evidence for fractional time evolution in glass materials. Chem. Phys. 284, ) 1. Gorenflo, R, Mainardi, F: Fractional calculus: integral and differential equations of fractional orders. In: Fractals and Fractional Calculus in Continuum Mechanics. Springer, New York 1997) 11. Herzallah, MAE, El-Sayed, AMA, Baleanu, D: On the fractional-order diffusion-wave process. Rom. J. Phys ), ) 12. Kilbas, AA, Srivastava, HM, Trujillo, TJ: Theory and Applications of Fractional Differential Equations. North-Holland Mathematics Studies, vol. 24, van Mill, J ed.). Elsevier, Amsterdam 26) 13. Miller, KS, Ross, B: An Introduction to the Fractional Calculus and Fractional Differential Equations. Wiley, New York 1993)

14 Al-Homidan et al. Advances in Difference Equations 213, 213:349 Page 14 of Oldham, KB, Spanier, J: The Fractional Calculus. Academic Press, New York 1974) 15. Podlubny, I: Fractional Differential Equations. Mathematics in Sciences and Engineering, vol Academic Press, San Diego 1999) 16. Sandev, T, Metzler, R, Tomovski, Z: Fractional diffusion equation with a generalized Riemann-iouville time fractional derivative. J. Phys. A, Math. Theor. 44, ) 17. Malek, NA, Ghanam, RA, Al-Homidan, S: Sensitivity of the pressure distribution to the fractional order α in the fractional diffusion equation submitted) 18. Samko, SG, Kilbas, AA, Marichev, OI: Fractional Integrals and Derivatives. Theory and Applications. Gordon & Breach, Yverdon 1993) 19. Kou, C, iu, J, Ye, Y: Existence and uniqueness of solutions for the Cauchy-type problems of fractional differential equations. Discrete Dyn. Nat. Soc. 21,Article ID ) 2. Furati, KM, Kassim, MD, Tatar, N-e: Existence and uniqueness for a problem involving Hilfer fractional derivative. Comput. Math. Appl. 646), ) 21. Baleanu, D, Diethelm, K, Scalas, E, Trujillo, JJ: Fractional Calculus Models and Numerical Methods. Series on Complexity, Nonlinearity and Chaos. World Scientific, Singapore 212) 22. Baleanu, D, Machado, JAT, uo, AC: Fractional Dynamics and Control. Springer, Berlin 212) / Cite this article as: Al-Homidan et al.: On a generalized diffusion equation arising in petroleum engineering. Advances in Difference Equations 213, 213:349