Rate Decline Analysis for Limited-Entry Well in Abnormally High-Pressured Composite Naturally Fractured Gas Reservoirs

Mingtao Wu 1,2, Xiaodong Wang 1,*, Wenqi Zhao 3, Lun Zhao 3, Meng Sun 3 and Hai Zhou 2 1 School of Energy Resources, China University of Geosciences (Beijing), Beijing 100083, China; wumt@cugb.edu.cn 2 Beijing Key Laboratory of Unconventional Natural Gas Geological Evaluation and Development Engineering, China University of Geosciences (Beijing), Beijing 100083, China; zhouhai@cugb.edu.cn 3 Research Institute of Petroleum Exploration and Development, PetroChina, Beijing 100083, China; zhaowenqi@petrochina.com.cn (W.Z.); zhaolun@petrochina.com.cn (L.Z.); mickfoleypia@163.com (M.S.) * Correspondence: wxd_cug@cugb.edu.cn; Tel.: +86-10-8232-2754


Introduction
Whether CO 2 geological sequestration, disposal of high-level radioactive waste, or exploitation of geothermal and hydrocarbon resources, all of them are associated with fractured reservoirs [1][2][3].Therefore, it is increasingly significant to understand fluid flow behavior through naturally fractured reservoir in hydrology, environment, and petroleum engineering.Especially, with the development of natural gas exploration and exploitation in the world, the proven reserves and production of naturally fractured gas reservoirs are increasing year by year [4].
The majority of naturally fractured gas reservoirs exhibit abnormally high-pressured as a result of disequilibrium compaction [5].In China, Kela-2 gas field in Tarim Basin and Moxi gas field in Sichuan Basin which have been developed for several years are both abnormally high-pressured fractured gas reservoirs.Due to the characteristics of abnormal pressure, the rock and fluid properties of naturally fractured gas reservoirs are commonly pressure-dependent.
With microseismic monitoring and transient pressure analysis, it is demonstrated that some naturally fractured gas reservoirs may exhibit composite properties in the radial direction.That is, a large number of natural fractures exist near the wellbore area while the outer region far away from well has fewer ones [6].Therefore, a composite model is appropriate to analyze the production performance of the well.Composite reservoir models were originally targeted to model special reservoir configuration consisting of two or more concentric zones with different diffusivities or fluid properties.As early as 1960, Hurst [7] developed the analytical solution for radially composite reservoirs in which the storativity ratio between two regions are identical.Two decades later, Olarewaju and Lee [8,9] proposed an analytical solution to analyze pressure transient behavior for wells in finite composite reservoir system.Besides the composite formation with single porosity media in inner and outer regions, many researchers [10][11][12][13][14] also developed the mathematical model for composite dual-porosity reservoirs.Subsequently, Kuchuk and Habashy [15] solved fluid flow problems in composite reservoirs by reflection and transmission theory of electromagnetics.In recent years, composite model is applied to describe stimulated region volume around well in unconventional reservoirs [16][17][18].
Unfortunately, most literature aimed at studying fully penetrated wells in composite reservoirs.Nevertheless, some wells in actual fields were partially penetrated especially in thick formations.For the research of partially penetrated wells, extensive works were carried out in hydrology and petroleum engineering.Muskat [19] first investigated steady flow through porous media for partially penetrated vertical well.Nisle [20] and Brons and Marting [21] studied the effect of partial penetration on the pressure drawdown and buildup of oil wells in an isotropic layer with no-flow boundary conditions at the top and bottom of the layer.Odeh [22] established the steady-state flow model of partially penetrated well by Fourier transforms.Kazemi and Seth [23] studied the combined effect of anisotropy and stratification with crossflow on pressure response for partially penetrated well.With Laplace and Hankel transforms, Streltsova [24][25][26] solved the unsteady-state flow problems for partially penetrated well.The Laplace transform was also used by Dougherty and Babu [27] to analyze hydraulics problem for a partially penetrated well in a dual-porosity reservoir.Bui et al. [28] developed an analytical solution to analyze transient pressure behavior of partially penetrated wells in naturally fractured reservoirs with infinite radial extent.Based on the solution of Bui et al. [28], a set of pressure derivative type curves were generated by Slimani and Tiab [29].Recently, Mishra et al. [30] studied the radial flow to a partially penetrated well with storage in an anisotropic confined aquifer.Biryukov and Kuchuk [31] presented an analytical pressure transient solution for a limited-thickness open cylindrical interval on a nonpermeable cylindrical wellbore.In their model, the mixed boundary value problems were investigated.Javandel [32] proposed a semianalytical solution for partial penetration in two layer aquifers.Dejam et al. [33] studied the effect of a constant top pressure on the pressure transient analysis of a partially penetrated well in an infinite-acting fractured reservoir with wellbore storage and skin factor effects.However, most of them merely investigated transient pressure behavior for partially penetrated well in oil reservoirs or aquifers.Despite the great efforts presented in the aforementioned literature, the rate decline analysis model for partially penetrated wells in composite gas reservoir is lacking, particularly in stress-sensitive composite naturally fractured gas reservoir.
The objective of this paper is to develop a semianalytical model for rate decline analysis of partially penetrated well in abnormally high-pressured composite naturally fractured gas reservoirs.By the definition of pseudopressure function and pseudotime factor, fluid and rock pressure-dependent properties were taken into consideration.Laplace and finite Fourier cosine transforms were used to solve the 2D diffusivity equation.Based on the proposed model, the effects of prevailing factors on production performance were investigated.Finally, some remarkable conclusions are drawn.

Preliminaries
Fracture tends to close when pressure drops during the process of fluid extraction.Thus, the permeability of fracture in naturally fractured gas reservoirs may notably decrease [34][35][36][37].The explicit relationship can be determined by a stress-sensitive experiment [35], and the corresponding experimental results are shown in Figure 1.Based on the curve matching of experimental results, a common approach to account for permeability variations can be written as: where, p i and p denote initial and current reservoir pressure, respectively; k i and k denote permeability at initial and current reservoir pressure, respectively; γ k denotes permeability stress-sensitivity coefficient and is assumed to be constant.

Preliminaries
Fracture tends to close when pressure drops during the process of fluid extraction.Thus, the permeability of fracture in naturally fractured gas reservoirs may notably decrease [34][35][36][37].The explicit relationship can be determined by a stress-sensitive experiment [35], and the corresponding experimental results are shown in Figure 1.Based on the curve matching of experimental results, a common approach to account for permeability variations can be written as: where, pi and p denote initial and current reservoir pressure, respectively; ki and k denote permeability at initial and current reservoir pressure, respectively; γk denotes permeability stress-sensitivity coefficient and is assumed to be constant.

Figure 1.
Experimental data [35] matching for permeability of naturally fractured gas reservoir.
Furthermore, gas flow through porous media is essentially nonlinear because the properties of gas are strongly dependent on pressure.Conventionally, pseudopressure and pseudotime methods are used to linearize the diffusivity equation of gas [38][39][40][41][42].In this paper, novel definitions of pseudopressure function and pseudotime factor are proposed incorporating pressure-dependent porosity and permeability as follows respectively: where, pp denotes gas pseudopressure function; ki and k denote permeability at initial and current reservoir pressure, respectively, following Equation (1); Zi and Z denote gas deviation factor at initial and current reservoir pressure, respectively; μgi and μg denote gas viscosity at initial and current reservoir pressure, respectively; β denotes gas pseudotime factor; φi and φ denote porosity at initial and current reservoir pressure, respectively; cti and ct denote total compressibility at initial and current reservoir pressure, respectively.Experimental data [35] matching for permeability of naturally fractured gas reservoir.
Furthermore, gas flow through porous media is essentially nonlinear because the properties of gas are strongly dependent on pressure.Conventionally, pseudopressure and pseudotime methods are used to linearize the diffusivity equation of gas [38][39][40][41][42].In this paper, novel definitions of pseudopressure function and pseudotime factor are proposed incorporating pressure-dependent porosity and permeability as follows respectively: where, p p denotes gas pseudopressure function; k i and k denote permeability at initial and current reservoir pressure, respectively, following Equation (1); Z i and Z denote gas deviation factor at initial and current reservoir pressure, respectively; µ gi and µ g denote gas viscosity at initial and current reservoir pressure, respectively; β denotes gas pseudotime factor; ϕ i and ϕ denote porosity at initial and current reservoir pressure, respectively; c ti and c t denote total compressibility at initial and current reservoir pressure, respectively.The calculations of pseudopressure and pseudotime factor need to solve material balance equation for abnormally high-pressured naturally fractured gas reservoir [43]: where, G p denotes cumulative production; G f and G t denote original gas in place that is stored in fracture system and total original gas in place, respectively; c m and c f denote matrix and fracture system compressibility, respectively.

Physical Model and Basic Assumptions
As shown in Figure 2, a partially penetrated vertical well is located in the center of a circular naturally fractured gas reservoir.The reservoir is assumed to have two regions: Inner region and outer region, in which reservoir properties are drastically different.The radius of the inner region is R 1 , and the radius of the whole system is R e .The inner region consists of fractures and matrix while the outer region consists of matrix only, provided that the transfer flow between matrix and fractures in inner region is pseudosteady state.Besides, both of the two regions possess permeability anisotropies.Apart from the aforementioned description, some other assumptions are elaborated as: (1) The reservoir is horizontal with uniform thickness of h; (2) The upper and bottom boundaries of reservoir are impermeable, and outer boundary is also closed; (3) The water phase is assumed to be immobile, which means the gas production process is single phase flow; (4) The well is a line source and only part of it produces, and its perforation interval is h w ; (5) The effect of gravity and capillary pressure are ignored.The calculations of pseudopressure and pseudotime factor need to solve material balance equation for abnormally high-pressured naturally fractured gas reservoir [43]: where, Gp denotes cumulative production; Gf and Gt denote original gas in place that is stored in fracture system and total original gas in place, respectively; cm and cf denote matrix and fracture system compressibility, respectively.

Physical Model and Basic Assumptions
As shown in Figure 2, a partially penetrated vertical well is located in the center of a circular naturally fractured gas reservoir.The reservoir is assumed to have two regions: Inner region and outer region, in which reservoir properties are drastically different.The radius of the inner region is R1, and the radius of the whole system is Re.The inner region consists of fractures and matrix while the outer region consists of matrix only, provided that the transfer flow between matrix and fractures in inner region is pseudosteady state.Besides, both of the two regions possess permeability anisotropies.Apart from the aforementioned description, some other assumptions are elaborated as: (1) The reservoir is horizontal with uniform thickness of h; (2) The upper and bottom boundaries of reservoir are impermeable, and outer boundary is also closed; (3) The water phase is assumed to be immobile, which means the gas production process is single phase flow; (4) The well is a line source and only part of it produces, and its perforation interval is hw; (5) The effect of gravity and capillary pressure are ignored.

Model Solution
Based on the assumptions described in Section 2.2, the dimensionless governing equations can be established and the detailed derivation is presented in Appendix A. As shown in Appendix B, by

Model Solution
Based on the assumptions described in Section 2.2, the dimensionless governing equations can be established and the detailed derivation is presented in Appendix A. As shown in Appendix B, by the Laplace transforms with respect to t D and the finite Fourier cosine transforms with respect to z D , the line source solution for dimensionless bottom-hole pressure is obtained as following: where, s is the time variable in Laplace domain, and the expressions of E n , F n , G n , H n , q n , T n , f (s) are as following: As shown above, f (s) for the pseudosteady-state matrix-fracture mass exchange is slightly different from the expression of Warren-Root's model [44] due to the consideration of anisotropy.It is noticed that if κ 1 equals to 1, the expression of f (s) simplifies to the same form as isotropic dual-porosity mediums: On the basis of superposition principle [45], the relationship between the dimensionless bottom-hole pseudopressure solution at a fixed constant rate inner boundary and the dimensionless flow rate solution at a fixed bottom-hole pressure can be written as: With a numerical Laplace inversion algorithm proposed by Stehfest [46], the solution of dimensionless production rate in real time domain can be obtained, and the dimensional production rate can be calculated based on the following equation: where, q g denotes production rate; k fhi denotes horizontal fracture permeability at initial reservoir pressure; α p is the unit conversion coefficient.

Model Verification
Since there is no relevant literature on partially penetrated well in composite naturally fractured gas reservoir, to validate the accuracy of the proposed model, the comparison of this model with an analytical solution of completely penetrated well in composite reservoir is implemented.
Although Prado and Da [10] mainly aimed at analyzing transient pressure, the dimensionless production rate solution can also be obtained based on the relationship of Equation ( 14).Therefore, we consider the circumstance of h w = h in our model, which means the well is penetrated completely.We compare the results obtained from our model with Prado and Da's.The dimensionless parameters used for validation are presented in Table 1.All of the dimensionless parameters are defined in Appendix A. Figure 3 shows the comparison of production rate curves.The lines represent the results obtained from the proposed model, and the dots represent the results obtained from Prado and Da's method.As we can see, there is a good agreement between these two models, which indicates that our proposed model for partially penetrated well in composite naturally fractured reservoir is accurate.Although Prado and Da [10] mainly aimed at analyzing transient pressure, the dimensionless production rate solution can also be obtained based on the relationship of Equation (14).Therefore, we consider the circumstance of hw = h in our model, which means the well is penetrated completely.We compare the results obtained from our model with Prado and Da's.The dimensionless parameters used for validation are presented in Table 1.All of the dimensionless parameters are defined in Appendix A. Figure 3 shows the comparison of production rate curves.The lines represent the results obtained from the proposed model, and the dots represent the results obtained from Prado and Da's method.As we can see, there is a good agreement between these two models, which indicates that our proposed model for partially penetrated well in composite naturally fractured reservoir is accurate.
q wD , this model −dq wD /dlnt D , this model q wD , by Prado and Da's method −dq wD /dlnt D , by Prado and Da' method

Parameters Sensitivity Analysis
Relevant parameters used to conduct sensitivity analysis are presented in Table 2.

Inner Region Radius
Figure 4 shows the rate decline curves and cumulative production curves for four different values of inner region radius under constant bottom-hole pressure.It can be seen from Figure 4 that the inner region radius has an obvious effect on the production performance.Due to the high permeability of the naturally fractured inner region, the production rate is higher at the larger inner region size in the initial stage.Nevertheless, the rate decreases rapidly after 700 days.This is mainly because the gas reserves in the inner region have been produced to a large extent, for cases with larger inner region radius, the outer region radius is smaller, leading to reduced gas reserves.Figure 5 illustrates the evolution of pseudotime with respect to production time, which embodies the degree of deviation from liquid flow behavior.

Mobility Ratio Between Outer and Inner Regions
Figures 6 and 7 present the impact of mobility ratio between outer and inner region regions on the production performance curves and pseudotime factor respectively.In these scenarios, we keep the fracture permeability k fh the same while changing the value of k m2h in order to create different mobility ratio.It can be seen that with all other parameters kept constant, the production curves and pseudotime curves exhibit similar characteristics as in Figures 4 and 5.The larger the value of M, the higher production rate, as M reflects the flow capacity of region 2 compared to region 1.Thus, there will be more supplement into region 1 to hold a high production rate when M is bigger.With the increase of mobility ratio, the diffusivity ratio between the outer and inner regions will also increase, and the pressure wave propagates faster, which leads to a notable decrease of pseudotime factor.

Mobility Ratio Between Outer and Inner Regions
Figures 6 and 7 present the impact of mobility ratio between outer and inner region regions on the production performance curves and pseudotime factor respectively.In these scenarios, we keep the fracture permeability kfh the same while changing the value of km2h in order to create different mobility ratio.It can be seen that with all other parameters kept constant, the production curves and pseudotime curves exhibit similar characteristics as in Figures 4 and 5.The larger the value of M, the higher production rate, as M reflects the flow capacity of region 2 compared to region 1.Thus, there Pseudo-time factor t, day

Mobility Ratio Between Outer and Inner Regions
Figures 6 and 7 present the impact of mobility ratio between outer and inner region regions on the production performance curves and pseudotime factor respectively.In these scenarios, we keep the fracture permeability kfh the same while changing the value of km2h in order to create different mobility ratio.It can be seen that with all other parameters kept constant, the production curves and pseudotime curves exhibit similar characteristics as in Figures 4 and 5.The larger the value of M, the higher production rate, as M reflects the flow capacity of region 2 compared to region 1.Thus, there Pseudo-time factor t, day The effect of inner region radius on pseudotime factor.
Appl.Sci.2019, 9, x; doi: FOR PEER REVIEW www.mdpi.com/journal/applsciwill be more supplement into region 1 to hold a high production rate when M is bigger.With the increase of mobility ratio, the diffusivity ratio between the outer and inner regions will also increase, and the pressure wave propagates faster, which leads to a notable decrease of pseudotime factor.

Fracture Permeability Anisotropy Factor
Figures 8 and 9 display the effect of fracture permeability anisotropy factor on production performance curves and pseudotime factor respectively.It is worth noting that fracture permeability will be more supplement into region 1 to hold a high production rate when M is bigger.With the increase of mobility ratio, the diffusivity ratio between the outer and inner regions will also increase, and the pressure wave propagates faster, which leads to a notable decrease of pseudotime factor.

Fracture Permeability Anisotropy Factor
Figures 8 and 9 display the effect of fracture permeability anisotropy factor on production performance curves and pseudotime factor respectively.It is worth noting that fracture permeability anisotropy factor mainly affects production rate in the early stage.The increase of anisotropy factor causes the slightly decrease of production rate.In late stage, the effect of κ 1 is weak.
Appl.Sci.2019, 9, x; doi: FOR PEER REVIEW www.mdpi.com/journal/applscianisotropy factor mainly affects production rate in the early stage.The increase of anisotropy factor causes the slightly decrease of production rate.In late stage, the effect of κ1 is weak.

Permeability Stress-Sensitivity Coefficient
Figures 10 and 11 depict the effect of permeability stress-sensitivity coefficient on production performance curves and pseudotime factor respectively.For larger stress-sensitivity coefficient, there will be a rapid rate decline and smaller cumulative production.This is because that stress sensitivity decreases the permeability of natural fracture and increases flowing resistance.With the increase of anisotropy factor mainly affects production rate in the early stage.The increase of anisotropy factor causes the slightly decrease of production rate.In late stage, the effect of κ1 is weak.

Permeability Stress-Sensitivity Coefficient
Figures 10 and 11 depict the effect of permeability stress-sensitivity coefficient on production performance curves and pseudotime factor respectively.For larger stress-sensitivity coefficient, there will be a rapid rate decline and smaller cumulative production.This is because that stress sensitivity decreases the permeability of natural fracture and increases flowing resistance.With the increase of The effect of fracture permeability anisotropy factor on pseudotime factor.

Permeability Stress-Sensitivity Coefficient
Figures 10 and 11 depict the effect of permeability stress-sensitivity coefficient on production performance curves and pseudotime factor respectively.For larger stress-sensitivity coefficient, there will be a rapid rate decline and smaller cumulative production.This is because that stress sensitivity decreases the permeability of natural fracture and increases flowing resistance.With the increase of stress-sensitivity coefficient, the pseudotime factor will decrease, which means the formation energy depletes faster.

Well Perforation Height
Figures 12 and 13 show the effect of well perforation height on production performance curves and pseudotime factor respectively.As we know, the perforation height determines the area of wellbore into which fluid flows.The smaller the hw, the bigger pressure drop and flow resistance.Appl.Sci.2019, 9, x; doi: FOR PEER REVIEW www.mdpi.com/journal/applscistress-sensitivity coefficient, the pseudotime factor will decrease, which means the formation energy depletes faster.

Well Perforation Height
Figures 12 and 13 show the effect of well perforation height on production performance curves and pseudotime factor respectively.As we know, the perforation height determines the area of wellbore into which fluid flows.The smaller the hw, the bigger pressure drop and flow resistance.

Well Perforation Height
Figures 12 and 13 show the effect of well perforation height on production performance curves and pseudotime factor respectively.As we know, the perforation height determines the area of wellbore into which fluid flows.The smaller the h w , the bigger pressure drop and flow resistance.Thus, the small perforation height leads to low production rate, and slow formation energy depletion as depicted in Figure 13.It is desirable to improve perforation height in premise of avoiding water encroachment.
Appl.Sci.2019, 9, x; doi: FOR PEER REVIEW www.mdpi.com/journal/applsciThus, the small perforation height leads to low production rate, and slow formation energy depletion as depicted in Figure 13.It is desirable to improve perforation height in premise of avoiding water encroachment.

Field Data Matching
In order to prove the practical application of the proposed model, the production data of an actual case from Moxi gas field in Sichuan Basin was used.The well was partially penetrated and a schedule of constant well bottom-hole pressure was adopted at the initial production stage.Thus, the Appl.Sci.2019, 9, x; doi: FOR PEER REVIEW www.mdpi.com/journal/applsciThus, the small perforation height leads to low production rate, and slow formation energy depletion as depicted in Figure 13.It is desirable to improve perforation height in premise of avoiding water encroachment.

Field Data Matching
In order to prove the practical application of the proposed model, the production data of an actual case from Moxi gas field in Sichuan Basin was used.The well was partially penetrated and a schedule of constant well bottom-hole pressure was adopted at the initial production stage.Thus, the

Field Data Matching
In order to prove the practical application of the proposed model, the production data of an actual case from Moxi gas field in Sichuan Basin was used.The well was partially penetrated and a schedule of constant well bottom-hole pressure was adopted at the initial production stage.Thus, the gas flow rate was characterized by a decreasing trend for a long time.The following equation was used to describe the relationship between porosity and reservoir pressure: where, ϕ i and ϕ denote porosity at initial and current reservoir pressure, respectively, following Equation (3); c ϕ denotes formation compressibility, and it can be calculated by the following equation: where, c m and c f denote matrix and fracture system compressibility, respectively, following Equation (4); ϕ m and ϕ f denote matrix and fracture system compressibility, respectively.All input parameters are listed in Table 3. Figure 14 illustrates a comparison between the field production rate data and the matching results obtained from the proposed model.In Figure 14, the solid line represents the matching result considering formation compressibility.The dash line represents the matching result not considering formation compressibility.That is, the formation compressibility equals to zero and the porosity is constant.As we can see from Figure 14, the real field data represented by the red circles is in good agreement with the result considering formation compressibility.The matching data is not good for the model without considering formation compressibility at late time.Therefore, it is necessary to consider formation compressibility and porosity stress-sensitivity in abnormally high-pressure gas reservoir production performance analysis.Based on Equation (A5), the control equations can be further transformed as following: where, κ 1 is the permeability anisotropy factor of fracture, which is defined in Table A1.
In the outer region, the single porosity model is employed.Thus, the control equation for outer region is as following: where, The analogous assumptions are made as: Thus, the control equation for outer region can be written as: By using the dimensionless variables listed in Table A1, the governing equations can be transformed into dimensionless form as follows.
Governing equations of inner region: Governing equations of outer region: Correspondingly, the initial conditions are: Inner boundary conditions are: Outer boundary conditions are: Interface boundary conditions are:

Appendix B. Derivation of the Solution
Firstly, the Laplace transforms with respect to t D for dimensionless fracture and matrix pseudopressure are defined as following, respectively: To eliminate the variable z D in the control equations, the finite Fourier cosine transforms are defined as follows: The general solutions of Equations (A28) and (A29) are: Combining the conditions of Equations (A30) to (A32), the solution in Laplace and Fourier domain is: E n I 0 (r D q n ) + F n K 0 (r D q n ) F n K 1 (q n ) − E n I 1 (q n ) cos(u n z wD ) (A35) where, E n = q n G n K 1 (R 1D q n ) + MT n H n K 0 (R 1D q n ) (A36) By the utility of inverse Fourier transforms, we can obtain: q 0 E 0 I 0 (r D q 0 )+F 0 K 0 (r D q 0 ) F 0 K 1 (q 0 )−E 0 I 1 (q 0 ) +2 ∞ n=1 sin( 12 u n h wD ) 1 2 u n h wD 1 q n E n I 0 (r D q n )+F n K 0 (r D q n ) F n K 1 (q n )−E n I 1 (q n ) cos(u n z wD ) cos(u n z D ) (A40) For limited-entry well, the equivalent pressure point can be written as follows by integral-averaging method.s p wD = 1 q 0 E 0 I 0 (r D q 0 )+F 0 K 0 (r D q 0 ) F 0 K 1 (q 0 )−E 0 I 1 (q 0 ) + 2 q n E n I 0 (r D q n )+F n K 0 (r D q n ) F n K 1 (q n )−E n I 1 (q n ) cos 2 (u n z wD ) (A41)

Figure 2 .
Figure 2. A schematic of partially penetrated well in composite naturally fractured gas reservoir.

Figure 2 .
Figure 2. A schematic of partially penetrated well in composite naturally fractured gas reservoir.

Figure 3 .
Figure 3.A comparison of the dimensionless production rate curves between this model and literature's method.

Figure 3 .
Figure 3.A comparison of the dimensionless production rate curves between this model and literature's method.

Figure 4 .
Figure 4.The effect of inner region radius on production performance curves.

Figure 5 .
Figure 5.The effect of inner region radius on pseudotime factor.

Figure 4 .Figure 4 .
Figure 4.The effect of inner region radius on production performance curves.

Figure 5 .
Figure 5.The effect of inner region radius on pseudotime factor.

Figure 6 .
Figure 6.The effect of mobility ratio on production performance curves.

Figure 7 .
Figure 7.The effect of mobility ratio on pseudotime factor.

Figure 6 .
Figure 6.The effect of mobility ratio on production performance curves.

Figure 7 .
Figure 7.The effect of mobility ratio on pseudotime factor.

Figures 8 and 9 1 M =0. 2 M =0. 3 M =0. 4 Figure 7 .
Figures8 and 9display the effect of fracture permeability anisotropy factor on production performance curves and pseudotime factor respectively.It is worth noting that fracture permeability

Figure 8 .
Figure 8.The effect of fracture permeability anisotropy factor on production performance curves.

Figure 9 .
Figure 9.The effect of fracture permeability anisotropy factor on pseudotime factor.

Figure 8 .
Figure 8.The effect of fracture permeability anisotropy factor on production performance curves.

Figure 8 .
Figure 8.The effect of fracture permeability anisotropy factor on production performance curves.

Figure 9 .
Figure 9.The effect of fracture permeability anisotropy factor on pseudotime factor.

Figure 10 .
Figure 10.The effect of stress-sensitivity coefficient on production performance curves.

Figure 11 .
Figure 11.The effect of stress-sensitivity coefficient on pseudotime factor.

1 γ 1 Figure 10 .
Figure 10.The effect of stress-sensitivity coefficient on production performance curves.

Figure 10 .
Figure 10.The effect of stress-sensitivity coefficient on production performance curves.

Figure 11 .
Figure 11.The effect of stress-sensitivity coefficient on pseudotime factor.

Figure 12 .
Figure 12.The effect of well perforation height on production performance curves.

Figure 13 .
Figure 13.The effect of well perforation height on pseudotime factor.

Figure 12 .
Figure 12.The effect of well perforation height on production performance curves.

Figure 12 .
Figure 12.The effect of well perforation height on production performance curves.

Figure 13 .
Figure 13.The effect of well perforation height on pseudotime factor.

Figure 13 .
Figure 13.The effect of well perforation height on pseudotime factor.

1 0s sin 1 2 u n h wD 1 2
cos(u n z D )dz D , p m2D (u n ) = p m2D cos(u n z D )dz D (A25)Thus, the inverse Fourier transforms are:p f D (z D ) = ∞ n=0 cos(u n z D ) N(n) p f D (u n ), p m2D (z D ) = transform presented by Equation (A25), we can obtain: n z D )dz D = − 1 u n h wD cos(u n z wD ) (A31) p f D = p m2D , ∂ p f D ∂r D = M ∂ p m2D ∂r D , r D = R 1D .(A32)

Table 1 .
The dimensionless parameters for model validation.

Table 1 .
The dimensionless parameters for model validation.

Table 2 .
The reservoir, fluid, and well parameters for sensitivity analysis.

Table A1 .
The definitions of dimensionless parameters.