Application of Homotopy Analysis Method in One-Dimensional Instability Phenomenon Arising in Inclined Porous Media
1Department of Applied Mathematics & Humanities, S. V. National Institute of Technology, Surat, India
During secondary oil recovery process when water is injected in inclined oil formatted area then phenomenon of instability occurs due to viscosity difference of water and oil. The non-linear partial differential equation for this instability phenomenon have been obtained. The Homotopy analysis method has been applied to this governing equation by using appropriate initial and boundary conditions. The guess value of saturation of injected water has been satisfying its initial and boundary conditions. The numerical value and graphical presentation are given by using Maple software and it is concluded that the saturation of injected water is increasing during instability phenomenon in inclined porous media when length of fingers and time increases.
At a glance: Figures
Keywords: instability phenomenon, homotopy analysis method, porous matrix, secondary oil recovery process
American Journal of Applied Mathematics and Statistics, 2014 2 (3),
Received Jaunary 21, 2014; Revised March 06, 2014; Accepted April 24, 2014Copyright © 2013 Science and Education Publishing. All Rights Reserved.
Cite this article:
- Patel, Kajal K., M. N. Mehta, and Twinkle R. Singh. "Application of Homotopy Analysis Method in One-Dimensional Instability Phenomenon Arising in Inclined Porous Media." American Journal of Applied Mathematics and Statistics 2.3 (2014): 106-114.
- Patel, K. K. , Mehta, M. N. , & Singh, T. R. (2014). Application of Homotopy Analysis Method in One-Dimensional Instability Phenomenon Arising in Inclined Porous Media. American Journal of Applied Mathematics and Statistics, 2(3), 106-114.
- Patel, Kajal K., M. N. Mehta, and Twinkle R. Singh. "Application of Homotopy Analysis Method in One-Dimensional Instability Phenomenon Arising in Inclined Porous Media." American Journal of Applied Mathematics and Statistics 2, no. 3 (2014): 106-114.
|Import into BibTeX||Import into EndNote||Import into RefMan||Import into RefWorks|
Since last few decades, the study of multiphase flow through porous media has great important in petroleum technology for oil recovery process. In oil recovery process oil can recover through three different stages like primary, secondary and tertiary or enhanced oil recovery. Primary oil recovery is the process in which oil comes out from oil formatted area by natural pressure and it produces only 10% to 15% of oil by natural pressure of the reservoir. Remaining oil can recover by injecting water, gas or any chemical. Hence this process is called secondary oil recovery process. In this paper, we consider the phenomenon of instability occurs during water injection in secondary oil recovery process as shown in Figure 1 when water is injected in oil formatted area. Due to the viscosity difference of water and oil, instead of regular displacement of whole front (common interface) the protuberances will occurs with irregular fingers in size and shape. Hence it is called fingering phenomenon. Injection of water is the principal form of the secondary oil recovery because supply of water is often plentiful, inexpensive and it is usually more stable frontal displacement than other form of secondary oil recovery. Due to water injection, oil will displace towards the production well in this way remaining oil can recover in secondary recovery process. In 1972, Bear considered two types of frontal displacement when two or more fluids in motion occupy the porous media. Firstly: Stable displacement and secondly: instable displacement . The stability of water flood depends on the mobility ratio between oil and water, heterogeneity of the porous medium and the interaction of several forces.
Instabilities may occur in both miscible and immiscible processes, and originate on the interface between two fluids (e.g. oil and water).For displacements according to the Muskat model , in which oil and water flow in separate macroscopic regions and gravitational forces are absent. Chuoke et al.  determined a condition for instability is that the mobility of the displacing water be higher than that of the displaced oil. He also accounted for capillary effects by defining an effective interfacial tension between the fluids in a porous medium analogous to the interfacial tension in the bulk fluids. The first application of the Hele-Shaw cell in viscous fingering problems were made by Saffman and Taylor . As both the velocity of the fluid in a porous medium and the velocity of a fluid in a Hele-Shaw cell obey Darcy’s law, Saffman and Taylor suggested that the behavior of a Hele-Shaw cell could represent flow in a porous medium . This phenomenon is analytically discussed first time by Schidegger and Johnson . They have defined a different approach to this phenomenon in which they have called the statistical of fingers as shown in Figure 2. In the statistical treatment of fingers only the average cross section area occupied by the schematic fingers had been considered, while size and shape of the individual fingers were neglected and the capillary mean pressure of the finger were also ignored. Verma  has obtained the statistical behavior of fingering in a heterogeneous porous media with mean capillary pressure. Patel and Mehta  have studied this phenomenon by converting the governing equation of fingering phenomenon in the form of Burger's equation and solution obtained in analytical form. Kinjal and Mehta  have given a series solution of nonlinear partial differential equation for fingering phenomenon in heterogeneous as well as in homogeneous porous media. Joshi and Mehta (2009) have discussed solution of governing equation by the Group Invariant method in homogenous porous media . Vyas and Mehta have discussed instability phenomenon with magnetic field effect . Meher and Mehta  have been obtained an approximate solution of governing equation of instability phenomenon by Adomain Decomposition method. Since governing equation of instability phenomenon is nonlinear partial differential equation of diffusion type, hence it is appropriate to apply Homotopy analysis method to find an approximate solution of the problem.
2. Statement of the Problem
For the sake of mathematical study, we consider one-dimensional instability phenomenon for which we choose a piece of cylindrical porous matrix from large natural field area and take vertical cross-sectional area of this small finite incline cylindrical porous matrix as a rectangle. Which is incline at small angle with x-axis and its three sides are impermeable expect for one end of cylinder which is designated as common interface at x = 0 as shown in Figure 2. Since the water is injected in oil formatted area and water and oil are flowing through porous matrix then it will satisfy Darcy's law  for low Reynolds number. Hence cylindrical pieces of incline porous matrix is filled with the oil and water is injected at x = 0 hence both will satisfy the equation of continuity.
When water is injected at common interface x = 0 then instead of regular displacement of the whole front, small protuberance will occurs due to the water following through interconnected capillaries in the porous matrix due to capillary pressure the irregular fingers occurred with irregular sizes and shapes. To find the area occupy by water in form of saturation of water in different irregular fingers as shown in Figure 2 by white dotted fingers. Schidegger and Johnson  suggested schematic presentation of fingers by replacing irregular fingers by rectangular fingers for mathematical study as shown in Figure 3. But still it is difficult to find the saturation occupied in schematic fingers for given length x and time t > 0. Hence for one-dimensional study of instability phenomenon by taking average cross-sectional area occupied by schematic finger is considered as saturation of injected water for the length of fingers x for given time t > 0 which is rectangular shape as shown by Figure 4. Here the injection of water is taking place at x= 0 to native formation of oil formatted (dark area) porous media as shown in Figure 4. We consider here that saturation of water at common interface occupied cross-sectional area is very small saturation. It is also considered for mathematical study that the initial saturation at common interface as average cross-sectional area is very small where. The comparative study have been carried out for instability phenomenon in homogeneous porous matrix with and without inclination.
This paper presents the instability phenomenon which occurs during the displacement process of two immiscible fluids through homogeneous porous media in secondary oil recovery process.
3. Mathematical Formulation
Since water is injected at common interface in inclined homogeneous porous matrix contenting oil which will displaced by injecting water. Hence water and oil both will satisfy Darcy's law given by Bear  which gives velocities of water and oil respectively as
Since in inclined porous matrix gravitational effect g and angle of inclination play important role in velocities of injected water and native oil so the second terms are added in above equations.
Where and represents velocity of injected water and native oil respectively, K is the permeability of the homogenous porous medium, is the angle of inclination of porous media measured with the horizontal line in positive direction, and are relative permeability of the injected water and native oil, which are the function of their saturations and ,andare pressures of injecting water and displaced native oil, and are constant densities of injected water and native oil, g is the acceleration due to gravity, andare the constant kinematic viscosities of injected water and native oil.
Since water and oil are following in a porous matrix through interconnected capillaries during the phenomenon of instability, due to capillary pressure of water and oil. Injected water and displaced native oil which will satisfy equation of continuity as given by Scheidegger  as
Where is the porosity of the homogeneous porous medium which is considered here as a constant.
It is also given that the sum of the saturation of the injected water and native oil is equal to unity  (i.e. fully saturated) written as
Further for the instability phenomenon it is necessary to understand the role of capillary pressure. Hence when less viscous water is injected at x = 0 in oil formatted inclined porous matrix of length x = L, the water can flow through interconnected capillaries made by effective pores and it is due to the difference of pressure of native oil and injected water. Which is also a function of saturation of injected water defined by Scheidegger  as
But Mehta  consider that the injected water linearly flowing through small interconnected capillary in homogeneous porous media. Hence he consider capillary pressure as linear function of saturation given by
Here minus sign shows the direction of capillary pressure opposite to the saturation of the injected water. This relation has been considered by many author [6, 7, 12, 16, 20] in their research work. Schidegger  and Bear  had defined this relationship between permeability and saturation of water and oil respectively as
where which is constant.
Now substituting the value of seepage velocity of injected water and velocity of native oil from equation (1) and (2) into the continuity equations of (3) and (4) the following coupled partial differential equations as
Eliminating from equations (9) by using relation (6) we get
Adding equation (10) and (11) and using relation (5) and integrating with respect to x we get,
Where V(t) is a constant of integration.
After simplification of (12), we get
Using above value of in equation (11) we get
Now Oroveanu  gave relation of mean pressure by
where is the constant mean pressure.
To determine the value of V(t) we use relation (16) in equation (13) and after simplification we get,
Hence using above value of V(t) in (14) and after simplification we get,
For more simplification of above equation of motion we use relation of and from equation (6) and (8) respectively in equation (18) we get,
Which is non-linear partial differential equation of motion for the saturation of injected water during instability phenomenon occurring in secondary oil recovery process in incline homogeneous porous matrix.
To solve nonlinear equation (19), it is necessary to choose appropriate boundary and initial condition at common interface x = 0. As par given statement, we consider here that when water is injected at common interface x = 0 the saturation of injected water is Sw0 which is very small for given time t > 0.
It is also considered in statement that the initial saturation of injected water is .
i.e. at t = 0 for x > 0, L is length of porous matrix where
4. The Solution with HAM
To find the solution of non-linear partial differential equation (19) for phenomenon of instability together with the condition (20)-(21). Choosing dimensionless variable
The equation (19) together with condition (20-21) will be converted into dimensionless form as
where , and units which is constant.
To solve the equation (22) together with condition (23) and (24) we use Homotopy analysis method. For this method first necessary step is to guess value of saturation which satisfies above conditions (23) and (24) as
Let the saturation of injected water at common interface X = 0 for embedded parameter q = 0 is then
Let us consider non-linear partial differential equation (22) can be expressed as,
Where is a nonlinear operator, is considered as unknown function which represent saturation of injected water for any distance X for given time for. We use auxiliary linear operator and initial first approximation of the saturation of injected water for q = 0 will be to construct the corresponding zeroth-order deformation equation. As the auxiliary linear operator which satisfies, where C1 is arbitrary constant. We construct a homotopy as ,
Enforcing the Homotopy (28) to be zero ,
Establish the zero-order deformation equation of instability phenomenon as ,
Where denote initial first approximation of the saturation of injected water of the exact solution of the saturation of injected water which we want to find. For which is an auxiliary parameter, is an auxiliary function , is an embedding parameter and is an auxiliary linear operator with the property
Since auxiliary parameter and auxiliary function is nonzero.
For andwe will get respectively,
According to (28) as the embedding parameter q increases from 0 to 1, solution varies continuously from the initial guess of the saturation of the injected water to the required solution (solution gradually improved from X = 0 to the end of rectangle as shown in Figure 4). As  considered saturation of injected water and its solution is assumed as
i.e. the saturation of injected water is a function of distance X, and time T for any parametric value q is expressed as, the saturation of injected fluid at common interface and sum of schematic presentation of fingers for saturation of injected fluid as given by Scheidegger for different value of parameter q is expressed as saturation of injected fluid at time T=0 is and sum of saturation of injected fluid fingers are , ,… at different time T for different value of parameter q. Here, the series (35) is called homotopy-series; the series (31) is called homotopy series solution of and is called the mth-order derivative of S*. Auxiliary parameter in homotopy-series (31) can be regard as iteration factor and is widely used in numerical computations. It is well known that the properly chosen iteration factor can ensure the convergence of homotopy series (31) is depending upon the value of , one can ensure that convergent of homotopy series, solution simply by means of choosing the proper value of as shown by Liao [8, 9, 10, 11]. If the auxiliary linear operator, the initial guess, the auxiliary parameter, the auxiliary function are so properly chosen, the series (31) converges at q=1.
Hence the saturation of injected fluid can be expressed as,
andwill be calculated by equation (38).
This must be one of the solution of original non-linear partial differential equation (22) for saturation of injected fluid. According to the definition (32), the governing equation can be deduced from the zero-order deformation equation (26). Define the vector
Differentiating equation (30) m times with respect to the embedding parameter q and then setting q=1 and finally dividing them by m!, we have the so-called mth order deformation equation for saturation of injected fluid water S will be as 
It should be emphasized that for is governed by the linear equation (32) with the linear boundary condition that come from the original problem, which can solve by symbolic computation software Maple as bellow. The rule of solution expression as given by equation (22) and equation (26), the auxiliary function independent of can be chosen as.
According to (30) and taking inverse of equation (34) we get,
In this way, we get for successively by using Maple software as,
Using initial guess value for the saturation of injected water from equation (22) and successive saturation of injected water into the native oil form (39) and (32) etc. and using in equation (34), we get
Whereare given by equation (28) and (29) respectively represents saturation of the injected water into the native oil at any time T for the horizontal distance X for T>0. The solution is an infinite series solution, which represents the approximate value of saturation of injected water for time T>0. Which is convergent at for auxiliary parameter.
Thus the saturation of injected water into the native oil is expressed in terms of exponentials of function of X and time T>0, which depends on first guess value of the solution (25) only. Equation (41) represents the saturation of injected water in oil formatted region of instability phenomenon in homogeneous porous media with small inclination when injected water is injected with external injecting force at common interface during secondary oil recovery process.
5. Numerical and Graphical Presentation with Inclination
Maple coding has been used to obtain numerical and graphical presentations of equation (41). Figure 5 represents the graph of vs. distance X for time T = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0 and fixed value, and . Table 1 indicates the numerical value for Figure 5. All tabular values used for the graphical representations of the instability phenomenon, which shows the behavior of the saturation of injected water. The convergence of the Homotopy series (41) is dependent upon the value of convergence-parameter [8, 9, 10, 11]. Therefore we choose proper value of the convergence-parameter to obtain convergent Homotopy-series solution .
Table 1. Numerical Values of the Saturation of Injected Water in Homogeneous Porous Matrix with Small Inclination θ=15
6. Deduction of instability phenomenon in homogeneous porous media without inclination
For discussing this case, we put in equation 1 and 2 so that the equation (22) becomes
As consider the same condition (23) and (24) and same guess value (25).
Above equation (42) is nonlinear partial differential equation of motion for the saturation of injected water during instability phenomenon occurring in secondary oil recovery process in homogeneous porous matrix. The solution of equation (42) can be obtained as .
Equation (43) represents the saturation of injected water during instability phenomenon in homogeneous porous matrix without inclination.
7. Numerical and Graphical Presentation without inclination
Numerical and graphical presentations of equation (43). Figure 6 represents the graph of vs. distance X for time T = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0 and fixed value, and Table II indicates the numerical value for Figure 6. The convergence of the Homotopy series (41) is dependent upon the value of convergence-parameter [8, 9, 10, 11]. Therefore we choose proper value of the convergence-parameter to obtain convergent Homotopy-series solution .
Table 2. Numerical Values of the Saturation of Injected Water in Homogeneous Porous Matrix Without Inclination
8. Comparative Study of Saturation of Injected Water in Instability Phenomenon with and without Inclined Homogeneous Porous Matrix
The saturation of injected water in horizontal homogeneous porous matrix will be faster than saturation of injected water in inclined homogeneous porous matrix due to absents of gravitational force and angle of inclination. The comparative study of instability phenomenon also shows by the Table 3 that numerical values of saturation of injected water in horizontal homogeneous porous matrix is more than numerical value in inclined homogeneous porous matrix.
Table 3. Comparative Numerical Value of Saturation of Injected Water With and Without Inclined Porous Matrix
The equation (41) represents solution of equation (22) together with condition (23) and (24) by using Homotopy analysis method. Which gives saturation of injected water in inclined porous matrix during instability phenomenon occurred in secondary oil recovery process which satisfy both condition (23) and (24). The solution is in form of exponential function (X) and function of time T for given perturb parameter. Hence solution shows that saturation of injected water increases exponentially as distance X (length of average cross-sectional area of schematic fingers) increases for given T > 0. From given finger 5 we can concluded that the saturation of injected water is not steadily increasing for small distance X due to the angle of inclination and gravitational effect but after some small distance the saturation of water is steadily exponentially increases as distance X increases for given time T > 0. When it is flowing through interconnected capillaries and external injecting force.
In deduction, for horizontal porous matrix for inclination angle θ=0 gives governing equation (42) together with same condition (23) and (24), the solution (43) represents the saturation of injected water in horizontal porous matrix with same guess value of injected water saturation. The graph of solution (43) given by Figure 6 is exponentially increases from common interface X = 0 because here gravitational force and angle of inclination does not play any role.
The comparative study of saturation of injected water for inclined and without incline porous matrix have been given by Figure 5 and 6. Which shows that saturation of injected water is increasing as distance X increase for given time T > 0. But due inclination the saturation of injecting water is slightly less than the saturation of injected water for case of without inclination for some distance X increases for given time T > 0.
Hence overall we can conclude that when water is injected in oil formatted area at common interface then instability phenomenon occurred and saturation of injected water is increase as distance X increase for given time T > 0 in case of incline porous matrix and horizontal porous matrix. But due to inclination the saturation of injected water is slightly less comparative to without inclination. Which is quite obvious for physical as well as experimental results.
|||Bear, J., Dynamics of Fluids in Porous Media, American Elsevier Publishing Co., New York, pp. 687-701, 1972.|
|||Booth, R.J., Miscible flow through porous media .PhD thesis, University of Oxford page-5-10, 2008.|
|||Chuoke, R.L., van Meurs, P. and van der Poel, C., The Instability of Slow, Immiscible Viscous Liquid-Liquid Displacements in Permeable Media, Trans., AIME 216, 188-94, 1959.|
|||Dake. L. P., The Practice of Reservoir Engineering, Elsevier Science B. V., Amsterdam, 1994.|
|||Darcy H., Les Fontaines publiques de lo ville de Dijon, Dalmont, Paris, 1856.|
|||Joshi M.S. and Mehta M.N., Solution by Group Invariant method of Instability phenomenon arising in fluid flow through porous media. International Journal of Engineering Research and Industrial Applications, Vol. 2, No. 1, 35-48, 2009.|
|||Kinjal R. P., Mehta M. N. and Patel T. R., An approximate solution of Instability phenomenon in heterogeneous porous media with mean pressure., Gen. Math. Notes., Vol. 10 (2), 9-21, 2012.|
|||Liao S. J., A kind of approximate solution technique which does not depend upon small parameters (II): an application in fluid mechanics, Int. J. Nonlinear. Mech, 32-815, 1997.|
|||Liao, S. J., An explicit, totally analytic approximation of Blasius viscous ﬂow problems. Intl J. Non-Linear Mech. 34, 759-778, 1999 a.|
|||Liao S. J., On the Homotopy analysis method for nonlinear problems. Appl. Math. Comput. Vol. 147, pp. 499-513, 2004.|
|||Liao, S. J., The proposed homotopy analysis technique for the solution of nonlinear problems. PhD thesis, Shanghai Jiao Tong University, 1992.|
|||Meher R. K. and Mehta M. N. Meher S. K., Instability phenomenon arising in double phase flow through porous media with capillary pressure., Int. J. of Appl. Math and Mech. 7 (15): 97-112, 2011.|
|||Mehta M. N., Asymptotic expansion of fluid flow through porous media, Ph.D. thesis South Gujarat University, Surat, India, 1977.|
|||Muskat M., Physical principles of oil production, Mc-graw-hill Book co., Inc New York, 1949.|
|||Oroveanu T., Scurgerea Fluiidelor Prin Medii Poroase Neomogene, Editura Academiei Republicii Populare Romine 92, 328, 1963.|
|||Patel T. R. and Mehta M. N., A classical solution of the Burger's equation arising into the instability phenomenon in double phase flow through homogenous porous media, International Journal of Applied Mathematics, Bulgaria, Vol. 25 (1), 2008.|
|||Saffman, P.G. and Taylor, G., The Penetration of a Fluid into a Porous Medium or Hele-Shaw Cell Containing a more Viscous Liquid, Proc. Roy. Soc. London A 245, 312-329, 1958.|
|||Scheidegger A. E., The physics of flow through Porous Media, University of Toronto Press, 1960.|
|||Scheidegger A. E. and Johnson E. F., The statistical behaviour of instabilities in displacement process in porous media., Canadian J. Physics, 39-326, 1961.|
|||Vyas N. B., Patel K. R., Mehta M. N. and Patel T. R., Power series solution of instability phenomena in double phase flow through porous media under magnetic field effect, Int. J. of Appl. Math and Mech. 7 (16): 1-10, 2011.|
|||Verma A. P., Statistical behavior of fingering in a displacement process in heterogeneous porous medium with capillary pressure. Can. J. Phys., 47 (3), 319-324, 1969.|