The existence of a TPG can generate a relatively high pressure gradient in the process of fluid flow in porous media in low-permeable reservoirs, and neglecting the QPGTs in the governing equations, by assuming a small pressure gradient for such a problem, can cause a significant error in predicting the formation pressure. Based on these concerns, in consideration of the QPGT, a moving boundary model of radial flow in low-permeable reservoirs with the TPG for the case of a constant flow rate at the inner boundary is constructed. Due to strong nonlinearity of the mathematical model, a numerical method is presented: the system of partial differential equations for the moving boundary problem is first transformed equivalently into a closed system of partial differential equations with fixed boundary conditions by a spatial coordinate transformation method; and then a stable, fully implicit finite difference method is used to obtain its numerical solution. Numerical result analysis shows that the mathematical models of radial flow in low-permeable reservoirs with TPG must take the QPGT into account in their governing equations, which is more important than those of Darcy’s flow; the sensitive effects of the QPGT for the radial flow model do not change with an increase of the dimensionless TPG.
1. Introduction
Due to a continuously decreasing crude oil output from conventional reservoirs and a high international gasoline price in recent years, unconventional reservoirs such as low-permeable reservoirs and shale oil and gas reservoirs have become urgent development resources in the petroleum industry. Consequently, considerable attention has been paid to the relevant research on the kinematic principles for the fluid flow behavior in these unconventional reservoirs [1–5] at present. Abundant experimental and theoretical analyses [6–13] have demonstrated that the fluid flow, in low-permeable porous media, does not obey the classical Darcy’s law: the seepage velocity is not proportional to the formation pressure gradient, and there exists a threshold pressure gradient (TPG) λ; the fluid flow happens only if the formation pressure gradient is larger than TPG.
Much research on these relevant moving boundary models has been conducted [14–22]. The computed formation pressure distributions corresponding to these moving boundary problems of the fluid flow in the porous media with TPG show big difference from the ones based on Darcy’s law (see Figure 1): the formation pressure gradient is much steeper, it decreases until up to zero at a certain value of a dimensionless distance from a well, that is, the position of a moving boundary, and the pressure distribution curve shows a property of compact support [21]; whereas, for Darcy’s flow problem, the formation pressure drop can propagate to any infinite distance transiently according to the exact analytical solution [20], and the formation pressure gradient is much more smooth. Actually, the pressure distribution difference between Darcy’s flow and fluid flow in the porous media with TPG can be explained through the angle α between the dimensionless formation pressure curve and the dimensionless distance at the place of moving boundary, as shown in Figure 1: for the case of the existence of TPG, the tangent of α is equal to the derivative of pressure with respect to distance at the place of moving boundary, that is, the TPG; and, as α decreases gradually and is equal to zero, it becomes a limit case, which corresponds to Darcy’s flow; that is, TPG is equal to zero. And if a constant value of TPG is given, the tangent of α will not change transiently. Therefore, it can be concluded that the formation pressure gradient is relatively much higher for the fluid flow in porous media with large TPG. The main physical reason is that the presence of TPG can make the formation pressure drop propagate more slowly, which causes a larger pressure gradient in a relatively shorter pressure disturbed distance from a production well.
Comparison of computed dimensionless formation pressure distribution curves.
It is well known that the rock porosity and fluid density are dependent on the formation pressure, and their formula are commonly expressed in the nonlinear forms of exponential functions [23]. Therefore, in the mathematical modeling, the deduced governing equation, by substituting these equations of state and a kinematic equation into a mass balance equation, always contains a nonlinear quadratic pressure gradient term (QPGT), where the coefficient of this term is proportional to the fluid compressibility. In general, the nonlinear QPGT is usually neglected by assuming small fluid compressibility or a small formation pressure gradient [24]. The error by neglecting the nonlinear term may be acceptable for most routine engineering applications in the development of conventional reservoirs. However, it has been realized that the linearization is not applicable for large values of time [25, 26], and certain operations such as hydraulic fracturing, high injection or production rates of wellbore, and large pressure pulse testing can generate a high formation pressure gradient.
For conventional models of fluid flow in porous media based on Darcy’s law, Odeh and Babu [27] studied the analytical solutions of the nonlinear model of slightly compressible fluid flow through porous media with consideration of QPGT by a Laplace transform and pointed out that the nonlinear solutions showed the difference in the pressure change for injection and pumping conditions. Wang and Dusseault [28] presented an analytical solution for pore pressure coupled with deformation in a porous medium by taking the quadratic gradient term into account, and the deviations from the existing solutions were identified in cases of high pressure gradients. Chakrabarty et al. [29, 30] gave the analytical dimensionless pressure solutions for radial flow systems by a Laplace transform and concluded that the linearized analysis by neglecting QPGT could lead to serious errors in some cases such as high injection rates. Braeuning et al. [31] studied the effect of QPGT on the variable-rate well-tests, and concluded that the linearization error depended on the magnitude of wellbore damage, pseudoskin, and a nonlinear flow parameter. Tong et al. [25, 26] presented the exact analytical solutions of nonlinear transient flow models and dual-porosity models including QPGTs by the generalized Weber transform and Hankel transform; it was concluded that the effect of the QPGT should be considered in large time well testing in reservoir engineering. Li et al. [32] took the effects of both QPGT and wellbore storage into account to build a mathematical flow model in a fractal multilayer reservoir and presented its analytical solution by a Laplace transform. Dewei et al. [33] obtained the analytical solution of a mathematical model of transient seepage flow including QPGT by a Laplace transform, and they demonstrated that QPGT could not be ignored for the unsteady flow model with consideration of a wellbore skin effect. Bai et al. [24] constructed a nonlinear dual-porosity model incorporating QPGT in fracture space and obtained its analytical solution using a Hankel transform, and their analysis showed that the presented nonlinear model appeared to be suitable for simulating naturally fractured reservoirs, subjected to a high injection or production rate or significant fracture compressibility. Nie et al. [34, 35] studied nonlinear flow models with a QPGT for both a double-porosity reservoir and a triple-porosity reservoir and obtained a solution through a variable substitution for linearization; it was found that the influence of QPGT was very distinct, and the parameter values of nonlinear model explanation law were much accurate than those of linear model explanation. Yao et al. [36] established the mathematical model of transient flow in a double-porosity fractal reservoir by considering QPGT and solved the model by a Laplace transform; it was demonstrated that the relative errors, caused by ignoring QPGT in the constructed model, may amount to 10%. Guo and Nie [37] discussed the origins of nonlinear flow issues in underground formations by using quadratic pressure gradients to deduce diffusion equations and concluded that the nonlinear models could describe fluid flow in underground formations realistically and accurately.
As far as we know, for the moving boundary problems of fluid flow in porous media with TPG, which are based on modified Darcy’s law [8] and have been widely involved in engineering applications for low-permeable reservoirs, the effect of QPGT has not yet been taken into account in their governing equations [14–22]. However, this nonlinear QPGT may play an important role in the temporal and spatial pressure profiles due to a relatively high formation pressure gradient, generated in fluid flow in porous media with TPG. What is more, in modern times, with the development of advanced analysis methods and improved resolutions of pressure measurement machines [29], it is necessary to study quantitatively the effect of QPGT on these moving boundary problems.
Hence, based on these concerns, in this paper, QPGT is incorporated in the governing equation for a basic moving boundary problem of radial flow in an infinite reservoir with TPG for the case of a constant flow rate at the inner boundary. Owing to the existence of the moving boundary and the nonlinearity of partial differential equations, it is really difficult to obtain its exact analytical solution. Here, we adopted a spatial coordinate transform based finite difference method to solve the nonlinear moving boundary model for the radial flow case. The numerical method has been strictly verified for solving such moving boundary problems with good accuracy by a one-dimensional flow case in the recently published literature [21]. Furthermore, the effect of this QPGT on the numerical solutions of the moving boundary problem can be discussed and analyzed quantitatively, by using the numerical results with respect to different values of a dimensionless TPG.
2. Mathematical Modeling by Considering QPGT
The problem considered involves radial flow in an infinite reservoir with TPG for the case of a constant flow rate at the inner boundary; the porous medium is homogeneous, isotropic, and isothermal; the single-phase horizontal flow does not have any gravity effect; for a basic problem, skin effect and wellbore storage are not considered here; finally, the fluid and porous medium are slightly compressible.
The fluid density is as follows [20]:(1)ρ=ρiexp-Cfpi-p,where ρ is the fluid density, kg·m^{−3}; ρi is the initial fluid density, kg·m^{−3}; pi is the initial pressure, MPa; p is the pressure, MPa; and Cf is the compression coefficient of the fluid, MPa^{−1}.
The porosity of the porous medium is as follows [20]:(2)ϕ=ϕiexp-Cϕpi-p,where ϕ is the porosity of the porous medium, fraction; ϕi is the initial porosity, fraction; and Cϕ is the compression coefficient of the porosity, MPa^{−1}.
The modified Darcy’s law for the fluid flow in the porous medium with TPG is as follows [8]:(3)υ=0,0≤dpdr≤λ,-kμ·dpdr-λ,dpdr>λ,where k is the permeability of the porous medium, 10^{−3}·μ m^{2}; μ is the fluid viscosity, mPa·s; r is the radial distance, m; υ is the seepage velocity, m·s^{−1}; and λ is the TPG, MPa·m^{−1}.
The continuous equation for the radial flow in the porous medium is as follows [15]:(4)-1r∂∂rrρυ=∂ρϕ∂t,0≤r≤st,where t is time, h and s is the moving boundary, m.
Substituting (1)–(3) into (4), the governing (mass balance) equation for the radial flow in the low-permeable reservoir, considering the nonlinear QPGT, can be deduced as follows (see Appendix):(5)∂2p∂r2+1r·∂p∂r-λr+Cf·∂p∂r2=μϕiCtk·∂p∂t,where Ct is the total compression coefficient, MPa^{−1}.
The initial condition is as follows: (6)s0=0,pt=0=pi.
The inner boundary conditions with constant flow rate are (7)2πkhrμ∂p∂r-λr=rw=q·B,where q is the constant flow rate, m^{3}·d^{−1}; h is the reservoir thickness, m; rw is the wellbore radius, m; and Bis the volume factor, dimensionless.
The moving boundary conditions are: (8)pr=st=pi,∂p∂rr=st=λ.
The moving boundary conditions (8) physically mean that the seepage flow only happens at the area near the wellbore inside the moving boundary, where the formation pressure gradient is larger than TPG; however, outside the moving boundary, the formation pressure gradient is smaller than TPG (equal to zero), so there is no seepage flow behavior, and the formation pressure also keeps the initial pressure; in addition, on the moving boundary, the pressure gradient is just equal to TPG; with the time increasing, the moving boundary also moves outside gradually. The existence of moving boundary is the main difference between the models of fluid flow in porous media with TPG and the classical Darcy’s flow models.
Equations (5)–(8) together form the moving boundary problem of radial flow in the infinite reservoir with TPG for the case of a constant flow rate at the inner boundary, in consideration of QPGT.
Define the following dimensionless variables: (9)rD=rrw,tD=3.6×10-3ktμϕiCtrw2,δ=srw,PD=khpi-p1.842qμB,λD=khrwλ1.842qμB,αD=1.842qμBCfkh,where rD is the dimensionless radial distance, tD is the dimensionless time, PD is the dimensionless pressure, αD is the dimensionless compressibility, λD is the dimensionless TPG, and δ is the dimensionless moving boundary. Consider (10)1rD∂PD∂rD+∂2PD∂rD2+λDrD-αD∂PD∂rD2=∂PDrD,tD∂tD,1≤rD≤δtD,(11)PDtD=0=0,(12)δ0=0,(13)rD∂PD∂rDrD=1=-1+λD,(14)∂PD∂rDrD=δtD=-λD,(15)PDrD=δtD=0.
Equations (10)–(15) form the dimensionless mathematical model, considering QTPG, for the radial flow in the infinite reservoir with TPG for the case of a constant flow rate at the inner boundary.
From (15), we have(16)PDδtD,tD=0.
Differentiating two sides of (16), with respect to tD, we have (17)∂PD∂tDrD=δtD+∂PD∂rDrD=δtD·∂δ∂tD=0.
Substituting (14) into (17) yields(18)∂PD∂tDrD=δtD=λD·∂δ∂tD.
Letting xD=δtD on both sides of (10) and then substituting (14) yield(19)∂PD∂tDrD=δtD=∂2PD∂rD2rD=δtD-αD·λD2.
By (18) and (19), the velocity of the moving boundary can be written as follows:(20)∂δ∂tD=1λD∂2PD∂rD2rD=δtD-αD·λD.
From (20), it can be seen that the existence of the QPGT (the dimensionless compressibility αD is not equal to zero) can slow down the velocity of the moving boundary of the radial flow in low-permeable reservoir with TPG.
3. Numerical Method
Due to the existence of the moving boundary in the dimensionless mathematical model, the flow region is not fixed and expands outward continuously with time increasing. In order to overcome this difficulty in the space discretization for the transient flow region with the moving boundary, a spatial coordinate transformation method is used as follows [21, 38]:(21)yrD,tD=rD-1δtD-1,1≤rD≤δtD.
Through (21), the dynamic flow region for the moving boundary problem 0,δtD can be transformed to a fixed region [0,1], and the dimensionless pressure PDrD,tD can be transformed as a new unknown function η of two variables y and tD, that is, ηy,tD equivalently; correspondingly, the differential variables can be transformed, respectively, as follows: (22)∂PD∂rD=∂η∂y·1δ-1,(23)∂2PD∂rD2=∂2η∂y2·1δ-12,(24)∂PD∂tD=-∂η∂y·yδ-1·∂δ∂tD+∂η∂tD.
Substituting (22)–(24) into (10) yields(25)∂η∂y·1δ-1·1yδ-1+1+∂2η∂y2·1δ-12+λDyδ-1+1=∂η∂tD-∂η∂y·yδ-1·∂δ∂tD.
Substituting (22)–(24) into (13)–(15) and (20), respectively, yields (26)∂η∂yy=0=-1+λD·δ-1,(27)∂η∂yy=1=-λDδ-1,(28)ηy=1=0,(29)∂δ∂tD=1λD∂2η∂y2y=1·1δ-12-αD·λD.
From (26), we have (30)δ-1=-11+λD·∂η∂yy=0.
Substituting (29) and (30) into (25) to cancel the variables ∂δ/∂tD and δ yields(31)∂η∂y·-11+λD·∂η∂yy=02-y1+λDsssss·-y1+λD·∂η∂yy=0+1ssssssss·∂η∂yy=021λD∂2η∂y2y=1·1+λD2-αD·λDssssssss·∂η∂yy=02+∂2η∂y2·-y1+λD·∂η∂yy=0+1∂η∂yy=0+λD11+λD2∂η∂yy=03-αD∂η∂y2·-y1+λD·∂η∂yy=0+1∂η∂yy=0=∂η∂tD·11+λD2·∂η∂yy=03·-y1+λD·∂η∂yy=0+1.
From (21), (11) can be transformed into(32)ηtD=0=0.
From (26) and (27), the following equation can be obtained [21]:(33)∂η∂yy=0-1+λDλD∂η∂yy=1=0.
Equations (31)–(33) and (28) together form a closed system of partial differential equations with the fixed boundary conditions with respect to PD(y,tD), which is equivalent to the dimensionless mathematical model, that is, (10)–(15). The new partial differential equation (31) shows strong nonlinearity, indirectly indicating the strong nonlinearity of the original, untransformed moving boundary problem incorporating the QPGT. Here, a stable, fully implicit finite difference method [21, 39] is used to obtain its numerical solutions. The first derivative is replaced by a first-order forward difference; the second derivative is replaced by the second-order central difference [21], and then the difference equation for (31) can be written as follows:(34)ηi+1j+1-ηij+1Δy-11+λD·η1j+1-η0j+1Δy2sssssssssssss-iΔy1+λD·-iΔy1+λD·η1j+1-η0j+1Δy+1sssssssssssss·η1j+1-η0j+1Δy21λDηN-2j+1-2ηN-1j+1Δy2·1+λD2ssssssssssssssssss-αD·λD·η1j+1-η0j+1Δy2+ηi-1j+1-2ηij+1+ηi+1j+1Δy2·-iΔy1+λD·η1j+1-η0j+1Δy+1η1j+1-η0j+1Δy+λD11+λD2η1j+1-η0j+1Δy3-αDηi+1j+1-ηij+1Δy2·-iΔy1+λD·η1j+1-η0j+1Δy+1η1j+1-η0j+1Δy-ηij+1-ηijΔtD·11+λD2·η1j+1-η0j+1Δy3·-iΔy1+λD·η1j+1-η0j+1Δy+1=0,where N denotes the total number of spatial grid subintervals with the same length; Δy is the length of a grid subinterval, which is equal to 1/N; i denotes the index of the spatial grid from the well; j denotes the index of a time step; and ΔtD denotes the time step size [21].
From (28), we have(35)ηNj+1=0.
Then, from (34) and (35), the difference equation corresponding to the (N-1)th spatial grid can be expressed as follows: (36)-ηN-1j+1Δy·-11+λD·η1j+1-η0j+1Δy2-N-1Δy1+λDsssssssss·-N-1Δy1+λD·η1j+1-η0j+1Δy+1sssssssss·η1j+1-η0j+1Δy21λDηN-2j+1-2ηN-1j+1Δy2·1+λD2sssssssssssss-αD·λD·η1j+1-η0j+1Δy2+ηN-2j+1-2ηN-1j+1Δy2·-N-1Δy1+λD·η1j+1-η0j+1Δy+1×η1j+1-η0j+1Δy+λD11+λD2η1j+1-η0j+1Δy3-αDηN-1j+1Δy2·-N-1Δy1+λD·η1j+1-η0j+1Δy+1×η1j+1-η0j+1Δy-ηN-1j+1-ηN-1jΔtD·11+λD2·η1j+1-η0j+1Δy3·-N-1Δy1+λD·η1j+1-η0j+1Δy+1=0.
The difference equation of (33) is as follows:(37)η1j+1-η0j+1Δy+1+λDλD·ηN-1j+1Δy=0.
From (32), the initial conditions are obtained as follows: (38)ηi0=0,i=0,1,2,…,N-1.
Equations (34), (36), and (37) together form N difference equations at the (j+1)th time step and also contain N unknown variables PDij+1i=0,1,2,…,N-1. The Newton-Raphson iterative method [21, 39] is used to numerically solve these nonlinear difference equations; when PDij+1i=0,1,2,…,N-1 are numerically solved for, let j be equal to j+1, and, in the same manner, the numerical solutions of N difference equations with respect to N unknown variables PDij+2i=0,1,2,…,N-1 at the j+2th time step can also be numerically solved for; the rest can be deduced by analogy [21, 39]. Eventually, numerical solutions for the transformed nonlinear partial differential equation system with respect to ηy,tD can be obtained.
The difference equation of (21) is (39)rDij+1=iΔy·δj+1-1+1.
The difference equation of (30) is(40)δj+1-1=-11+λD·η1j+1-η0j+1Δy.
Substituting (40) into (39) yields(41)rDij+1=-i·η1j+1-η0j+11+λD+1.
By (41), numerical solutions of ηy,tD can be transformed as the ones of PD(rD,tD) in the process of numerical solutions at every time step [21]. From (41), it can be further concluded that the presented spatial coordinate transformation lets the time-dependent space discretization in the spatial coordinate of rD be transformed as a time-independent space discretization in the new spatial coordinate of y and then makes the numerical solutions by the finite difference method more applicable and simpler [21].
4. Results and Discussions4.1. Effect of QPGT under Different Values of TPG
Figures 2–4 show the effect of the QPGT on numerical solutions of the model, with respect to the dimensionless formation pressure distribution, dimensionless transient wellbore pressure, and dimensionless transient distance of moving boundary, respectively, under different values of the dimensionless TPG. From Figures 2–4, it can be clearly seen that with an increase of the dimensionless TPG, the effect of the QPGT on the mathematical model solutions become more and more obvious; the reason can be explained as follows: the bigger the TPG, the steeper the formation pressure gradient, which generates the solution deviations, resulting from neglecting the QPGT in the governing equation, more seriously. Moreover, from Figures 2–4, it can also be indicated that the dimensionless transient wellbore pressure for the case of neglecting the QPGT is larger than that for the case of considering the QPGT; and the dimensionless transient distance of moving boundary for the case of neglecting the QPGT is also larger than that for the case of considering the QPGT. The QPGT can slow down the velocity of the moving boundary in low-permeable reservoir with TPG.
The effect of QPGT on dimensionless formation pressure distribution under different values of dimensionless TPG.
The effect of QPGT on dimensionless transient wellbore pressure under different values of dimensionless TPG.
The effect of QPGT on dimensionless transient distance of moving boundary under different values of dimensionless TPG.
Tables 1 and 2 show the selected data of calculated relative errors εr from the already computed dimensionless formation pressure distribution and dimensionless transient wellbore pressure between the two cases considering and not considering the QPGTs, under different values of dimensionless TPG, which correspond to Figures 2 and 3, respectively. The relative error εr is equal to the absolute error divided by the magnitude of the exact value, that is, the solutions of the models considering the QPGT. In order to clearly figure out the change behavior of relative errors, the curves regarding the relative errors by the data from Tables 1 and 2 are also plotted in Figures 5 and 6, respectively.
Comparison data for the computed dimensionless formation pressure distribution and relative errors.
Computed dimensionless formation pressure distribution PD (rD, 1×105)
rD
λD=0.034
λD=0.067
λD=0.1
PD (αD=0.005)
PD (αD=0)
εr/%
PD (αD=0.005)
PD (αD=0)
εr/%
PD (αD=0.005)
PD (αD=0)
εr/%
1
23.48
25.17
7.20
25.73
28.88
12.24
27.32
32.65
19.51
50
12.55
13.73
9.40
14.42
16.94
17.48
15.17
19.72
29.99
100
9.09
10.05
10.56
9.68
11.99
23.86
9.06
13.28
46.58
150
6.45
7.34
13.80
5.72
7.91
38.29
3.61
7.63
111.36
200
4.22
5.03
19.19
2.11
4.15
96.68
0
2.4
*
250
2.24
2.98
33.04
0
0.65
*
0
0
*
300
0.45
1.14
153.33
0
0
*
0
0
*
350
0
0
*
0
0
*
0
0
*
^{*}denotes the null data.
Comparison data for the computed dimensionless transient wellbore pressure and relative errors.
Relative error curves for computing dimensionless formation pressure distribution by neglecting quadratic pressure gradient term.
Relative error curves for computing dimensionless transient wellbore pressure by neglecting quadratic pressure gradient term.
From Figure 5, it can be clearly indicated that, for any case with the same value of the dimensionless distance, the larger the dimensionless TPG, the larger the relative error of the dimensionless formation pressure; furthermore, the larger dimensionless TPG can make the relative errors grow more quickly with an increase of dimensionless distance. For the cases with the values of the dimensionless TPG λD=0.034, 0.067, and 0.1, the relative errors have largely exceeded 5% in the whole region.
From Figure 6, it can be clearly indicated that, for any case with the same value of dimensionless time, the larger the dimensionless TPG, the larger the relative error of the dimensionless transient wellbore pressure; and the larger dimensionless TPG can make the relative errors grow more quickly with an increase of dimensionless time. For the cases with the values of the dimensionless TPG λD=0.067 and 0.1, the relative errors have largely exceeded 5% in the whole time.
From Figures 5 and 6, it can also be concluded that the relative errors corresponding to the case with the value of the dimensionless TPG λD=0.034, which is relatively close to Darcy’s flow, stay a lower level for the change of relative errors of the dimensionless transient wellbore pressure with the increase of dimensionless time, whereas the change of relative errors of the dimensionless formation pressure distribution with the increase of dimensionless distance is still very considerable.
In conclusion, the QPGT has a great effect on the mathematical model solutions of radial flow in low-permeable reservoirs with TPG, especially for a large value of the dimensionless TPG; the larger the dimensionless distance and dimensionless time, the bigger the effect of the QPGT on the dimensionless formation pressure and the dimensionless transient wellbore pressure. In summary, the mathematical models of radial flow in low-permeable reservoirs with TPG should take the QPGT into account in their governing equations, which is more important than those of Darcy’s flow.
4.2. Sensitive Effect of QPGT: The Parameter of Dimensionless Compressibility
Figures 7, 8, and 9 show the sensitive effect of the QPGT on the formation pressure distribution, transient wellbore pressure, and transient distance of moving boundary, respectively. From these figures, it can be seen that the larger the value of the dimensionless compressibility, the smaller the values of the dimensionless formation pressure, the dimensionless transient wellbore pressure, and the dimensionless transient distance of the moving boundary. Moreover, the sensitive effects of the QPGT on the formation pressure distribution, transient wellbore pressure, and transient distance of the moving boundary for the radial flow model do not change with an increase of the dimensionless TPG.
The sensitive effect of QPGT on dimensionless formation pressure distribution.
The sensitive effect of QPGT on dimensionless transient wellbore pressure.
The sensitive effect of QPGT on dimensionless transient distance of moving boundary.
5. Conclusions
The existence of a TPG can lead to relatively steep formation pressure gradients, and, thus, it is not appropriate to neglect the QPGT for the fluid flow in porous media with the TPG. In consideration of the QPGT in the governing equation, a numerical method is presented for investigation of the effect of QPGT on a moving boundary problem of radial flow in low-permeable reservoirs with TPG for the case of a constant flow rate at the inner boundary. Numerical result analysis shows that the QPGT plays an important role in the model solutions with respect to the formation pressure distribution, transient wellbore pressure, and transient distance of the moving boundary. The mathematical models of radial flow in low-permeable reservoirs with TPG should take the QPGT into account in their governing equations, which is more important than those of Darcy’s flow. Besides, the sensitive effects of the QPGT on the numerical solutions for the radial flow model do not change with an increase of the dimensionless TPG. The presented research supports theoretical foundations for further improving technologies of well testing and numerical simulation in developing low-permeable reservoirs.
Appendix
Equation (1) can be rewritten as follows: (A.1)p=1Cflnρ-1Cflnρi+pi.
Differentiating two sides of (A.1) with respect to r, we have (A.2)∂p∂r=1Cf·ρ∂ρ∂r.
Equation (A.2) can be rewritten as follows:(A.3)∂ρ∂r=Cf·ρ∂p∂r.
In the same manner as above, from (1) and (2), the following equations can also be deduced as(A.4)∂ρ∂t=Cf·ρ∂p∂t,∂ϕ∂t=Cϕ·ϕ∂p∂t.
The left-hand side of (4) can be expanded as follows:(A.5)-1r·∂∂rr·ρ·υ=kμ·1r·∂∂rr·ρ·∂p∂r-λ=kμ·1r·ρ·∂p∂r-λ+kμ·1r·r·∂ρ∂r·∂p∂r-λ+kμ·1r·r·ρ·∂2p∂r2=kμ·1r·ρ·∂p∂r-λ+kμ·1r·r·Cf·ρ·∂p∂r·∂p∂r-λ+kμ·1r·r·ρ·∂2p∂r2=kρμ∂2p∂r2+1r·∂p∂r-λr+Cf·∂p∂r·∂p∂r-λ=kρμ×∂2p∂r2+1r·∂p∂r-λr︸MainTerm+Cf·∂p∂r2︸QuadraticGradientTerm-Cf·λ·∂p∂r︸SmallTerm.
Because λ≪1 and Cf≪1, the small term in the right-hand side of (A.5) can be neglected, and then (A.5) can be rewritten as follows:(A.6)-1r·∂∂rr·ρ·υ=kρμ∂2p∂r2+1r·∂p∂r-λr︸MainTerm+Cf·∂p∂r2︸QuadraticGradientTerm.
In (A.6), the QPGT is retained for the deduction of the governing equation.
Expanding the right-hand side of (4) yields(A.7)∂ρϕ∂t=ρ∂ϕ∂t+ϕ∂ρ∂t.
Substituting (A.4) into the right-hand side of (A.7) yields(A.8)∂ρϕ∂t=ρ·Cϕ·ϕ·∂p∂t+ϕ·Cf·ρ·∂p∂t=ρ·ϕ·∂p∂t·Cϕ+Cf=ρ·ϕ·∂p∂t·Ct.
Substituting (A.6) and (A.8) into (4), the governing equation in consideration of the QPGT can be obtained as follows: (A.9)kρμ∂2p∂r2+1r·∂p∂r-λr︸MainTerm+Cf·∂p∂r2︸QuadraticGradientTerm=ρ·ϕ·∂p∂t·Ct.
Equation (A.9) can be equivalently simplified, by canceling the variable ρ in its both sides, as follows:(A.10)∂2p∂r2+1r·∂p∂r-λr+Cf·∂p∂r2=μϕCtk∂p∂t.
Conflict of Interests
Wenchao Liu and Jun Yao declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors would like to acknowledge the funding by the Projects Grants nos. 11102237 and 51404232, both sponsored by the Natural Science Foundation of China (NSFC), and the Program for Changjiang Scholars and Innovative Research Team inUniversity (Grant no. IRT1294).
HuangY. Z.YangZ. M.HeY.WangX. W.An overview on nonlinear porous flow in low permeability porous MediaMonteiroaP. J. M.RycroftcC. H.BarenblattG. I.A mathematical model of fluid and gas flow in nanoporous mediaYuR. Z.BianY. N.LiY.ZhangX. W.YanJ.WangH. C.WangK. J.Non-Darcy flow numerical simulation of XPJ low permeability reservoirGuoJ.ZhangS.ZhangL.QingH.LiuQ.Well testing analysis for horizontal well with consideration of threshold pressure gradient in tight gas reservoirsZengB. Q.ChengL. S.LiC. L.Low velocity non-linear flow in ultra-low permeability reservoirCaiJ. C.YuB. M.A discussion of the effect of tortuosity on the capillary imbibition in porous mediaYaoY. D.GeJ. L.Characteristics of non-Darcy flow in low-permeability reservoirsCivanF.Xiang'anY.HaoguangW.LijuanZ.RenbaoZ.YongpanZ.Low pressure gas percolation characteristic in ultra-low permeability porous mediaHaoF.ChengL. S.HassanO.HouJ.LiuC. Z.FengJ. D.Threshold pressure gradient in ultra-low permeability reservoirsWeiX.QunL.ShushengG.ZhimingH.HuiX.Pseudo threshold pressure gradient to flow for low permeability reservoirsCaiJ.HuX.StandnesD. C.YouL.An analytical model for spontaneous imbibition in fractal porous media including gravityCaiJ. C.A fractal approach to low velocity non-Darcy flow in a low permeability porous mediumPascalH.Nonsteady flow through porous media in the presence of a threshold gradientWuY. S.PruessK.WitherspoonP. A.Flow and displacement of Bingham non-Newtonian fluids in porous mediaFuquanS.CiqunL.FanhuaL.Transient pressure of percolation through one dimension porous media with threshold pressure gradientG. Q. FengQ. G. LiuG. Z. ShiZ. H. LinAn unsteady seepage flow model considering kickoff pressure gradient for low-permeability gas reservoirsLuJ.GhedanS.Pressure behavior of vertical wells in low-permeability reservoirs with threshold pressure gradientLuJ.Pressure behavior of uniform-flux hydraulically fractured wells in low-permeability reservoirs with threshold pressure gradientLiuW.YaoJ.WangY.Exact analytical solutions of moving boundary problems of one-dimensional flow in semi-infinite long porous media with threshold pressure gradientYaoJ.LiuW. C.ChenZ. X.Numerical solution of a moving boundary problem of one-dimensional flow in semi-infinite long porous media with threshold pressure gradientLiuW.YaoJ.ChenZ.Analytical solution of a double moving boundary problem for nonlinear flows in one-dimensional semi-infinite long porous media with low permeabilityMarshallS. L.Nonlinear pressure diffusion in flow of compressible liquids through porous mediaBaiM.MaQ.RoegiersJ.A nonlinear dual-porosity modelCaoX. L.TongD. K.WangR. H.Exact solutions for nonlinear transient flow model including a quadratic gradient termTongD. K.ZhangH. Q.WangR. H.Exact solution and its behavior characteristic of nonlinear dual-porosity modelOdehA. S.BabuD. K.Comparison of solutions of the nonlinear and linearized diffusion equationsWangY.DusseaultM. B.The effect of quadratic gradient terms on the borehole solution in poroelastic MediaChakrabartyC.Farouq AliS. M.TortikeW. S.Effects of the nonlinear gradient term on the transient pressure solution for a radial flow systemChakrabartyC.AliS. M. F.TortikeW. S.Analytical solutions for radial pressure distribution including the effects of the quadratic-gradient termBraeuningS.JelmertT. A.VikS. A.The effect of the quadratic gradient term on variable-rate well-testsLiW.LiX.LiS.LiQ.The similar structure of solutions in fractal multilayer reservoir including a quadratic gradient termDeweiM.AilinJ.ChengyeJ.QianZ.DongboH.Research on transient flow regulation with the effect of quadratic pressure gradientNieR. S.GeF.LiuY. L.The researches on the nonlinear flow model with quadratic pressure gradient and its application for double porosity reservoirProceedings of the International Forum on Porous Flow and Applications: Flow in Porous Media: From Phenomena to Engineering and Beyond2009Wuhan, ChinaNieR. S.JiaY. L.YuJ.LiuY. L.The transient well test analysis of fractured- vuggy triple-porosity reservoir with the quadratic pressure gradient termProceedings of the Latin American and Caribbean Petroleum Engineering Conference2009Cartagena, ColombiaCartagena de IndiasYaoY. D.WuY. S.ZhangR. L.The transient flow analysis of fluid in a fractal, double-porosity reservoirGuoJ. C.NieR. S.Nonlinear flow model for well production in an underground formationCrankJ.BurdenR. L.FairesJ. D.