Transient Pressure Behavior of a Horizontal Well in a Naturally Fractured Gas Reservoir with Dual-Permeability Flow and Stress Sensitivity Effect

Wang, Guiqin; Zhang, Rui; Cui, Liangliang

doi:https://doi.org/10.1155/2022/5993701

Geofluids

On this page

Abstract Introduction Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Special Issue

Nanoconfined Fluid Occurrence and Flow Behavior in Unconventional Gas/Oil Reservoirs

View this Special Issue

Research Article | Open Access

Volume 2022 | Article ID 5993701 | https://doi.org/10.1155/2022/5993701

Transient Pressure Behavior of a Horizontal Well in a Naturally Fractured Gas Reservoir with Dual-Permeability Flow and Stress Sensitivity Effect

Guiqin Wang,¹Rui Zhang,¹and Liangliang Cui¹

Academic Editor: Zheng Sun

Received28 Feb 2022

Accepted11 Apr 2022

Published30 Apr 2022

Abstract

Dual-permeability flow and stress sensitivity effect are two fundamental issues that have been widely investigated in transient pressure analysis for horizontal wells. However, few attempts have been made to simulate the combined effects of dual-permeability flow and stress-dependent fracture permeability on the pressure transient dynamics of a horizontal well in a naturally fractured gas reservoir. In this approach, an analytical model is proposed to integrate the complexities of pressure-dependent PVT properties, dual-permeability flow behavior, and stress-dependent fracture permeability characteristics. The nonlinearity of the mathematical model is weakened by using Pedrosa’s transform formulation. Then, the Laplace integral transformation and separation of variables are applied to solve the model. Based on the solution of the mathematical model, a series of new-type curves are drawn to make a precise observation of different flow regimes. The main differences between the proposed model and the traditional models are discussed, and the effects of the permeability modulus of fractures, storability ratio, interporosity flow factor, and skin factor on transient pressure response are also examined. The results show that there are obvious differences in transient pressure dynamic curves between the proposed model and traditional models. The stress sensitivity effect plays a significant role in the intermediate flow period and the late-time pseudoradial flow period. The dual-permeability flow behavior mainly affects the early transient and interporosity flow stages. The proposed model can accurately simulate the transient pressure behaviors of a horizontal well in a naturally fractured gas reservoir with a dual-permeability flow and stress sensitivity effect. The novel model can be used to interpret pressure signals with accurate matching results and more reasonable interpreted parameters.

1. Introduction

Transient pressure analysis of a horizontal well in a naturally fractured gas reservoir is greatly affected by fracture seepage parameters and stress-dependent fracture permeability. For naturally fractured reservoirs, the fractures are always with heterogeneities [1, 2]. An experimental study on fracture stress sensitivities proves that stress-dependent fracture permeability significantly affects the transient pressure response [3]. So, it is pretty essential to propose a comprehensive model to capture the transient behavior of horizontal wells in naturally fractured gas reservoirs.

Because it is often impossible to describe the complex fractures precisely, continuum models are proposed to capture the flow behavior of this kind of reservoir. Much research has been done on theoretical models of vertical wells in naturally fractured reservoirs. Barenblatt et al. [4] proposed a classic double-porosity and single-permeability model to study vertical well production in porous media reservoirs. This model assumes that a naturally fractured reservoir is composed of two completely overlapping continua, porous matrix, and fractures. Warren and Root [5] expanded Barenblatt et al.’s approach to cover the independent physical properties of fracture and matrix. In their model, naturally fractured formation is formed by matrix blocks, which is separated by uniform and orthogonal fractures. Besides, the pseudosteady-state interporosity flow is firstly adopted to simulate mass transfer between fracture and matrix systems. After that, Kazemi et al. [6], de Swaan [7], Raghavan and Ohaeri [8], Serra et al. [9], Jalali and Ershaghi [10], Wu and Pruess [11], Bui et al. [12], Wu et al. [13], Kuchuk et al. [14], Jia et al. [2], and Wang et al. [15] proposed their own dual-media model for naturally fractured reservoirs with the consideration of transient interporosity flow behavior. These models assume that the fracture system is the only flow pathway directly connected with wellbore by ignoring the flow from the matrix system to the wellbore. Because the dual-porosity and single-permeability model is no longer applicable in naturally fractured reservoirs, a great deal of work have instead been directed at using dual-porosity and dual-permeability models [14, 16–18] to describe the flow behavior, which assumes that both the fracture and matrix systems are the flow pathway directly connected with the wellbore, also considered the pseudosteady-state and transient-state interporosity flow between matrix and fracture systems. However, most of the models ignored the effect of the stress-dependent fracture permeability on transient pressure response. Although a great deal of work has been done on theoretical models of naturally fractured gas wells considering stress sensitivity of fracture permeability, most of them are restricted to the dual-porosity and single-permeability flow problem [3, 15, 19–21].

In recent years, horizontal wells have been increasingly applied to some naturally fractured gas reservoirs. Research on the transient pressure behavior of this kind of well has become increasingly popular among engineers [22–30]. However, transient pressure analysis for horizontal wells is commonly performed assuming that permeability for natural fractures remains constant, which might not be physically applicable for stress-sensitive reservoirs. Besides, dual-permeability flow is seldom considered in their models. In general, the big challenge of analyzing the transient pressure response of a horizontal well is that the dual-permeability flow behavior and stress-dependent fracture permeability should be all incorporated in the mathematical models.

This paper presented a novel semianalytical model to examine the combined effects of dual-permeability flow behavior and stress-dependent fracture permeability on the transient pressure response of a horizontal well in a naturally fractured gas reservoir. The nonlinearity of the governing equations caused by the stress sensitivity of fracture permeability is eliminated using Pedrosa’s [31] transform formulation. With Laplace transform and separation of variables, we got the analytical solution of the mathematical model. A series of new transient pressure dynamic curves are drawn to observe different flow regimes based on the solution. Then, differences between the proposed model and traditional models are discussed and the effects of some critical parameters on transient pressure response are also analyzed with the proposed model.

2. Methodology

2.1. Model Assumption

As shown in Figure 1, the naturally fractured gas reservoir is composed of fracture and matrix systems and the physical properties of the two systems are independent. A radial cylindrical dual-porosity and dual-permeability medium reservoir is considered in which a single horizontal well is located at the center, completely penetrating the formation. The matrix/fracture flow is schematically described in Figure 2. In this study, the proposed model assumes that both the fracture and matrix systems are the flow pathway directly connected with the wellbore and fluids in the fracture and matrix systems first flow into the horizontal wellbore, followed by the matrix-fracture interporosity flow. Some simplifying physical model assumptions for the derivation of the governing equation are listed as the horizontal well produced with the constant production rate in a naturally fractured gas reservoir. The external boundaries of the top and bottom are assumed to be closed, and the lateral boundary is assumed to be infinite. The matrix-fracture interporosity flow in the reservoir is described by the pseudosteady-state model [5, 25, 26]. Fluid flow follows the law of Darcy seepage, and stress-dependent fracture permeability is considered. Also, capillary and gravity forces are neglected to simplify the model.

2.2. Mathematical Model

The PVT properties, such as fluid viscosity and volume factor of the gas phase, are quite sensitive to formation pressure. In this section, the pseudopressure transformation is used to capture pressure-dependent PVT properties and reduce the nonlinearity of governing differential equations. The definitions of pseudopressure and pseudotime are given by where is the pressure, MPa; is the pseudopressure of the fracture, MPa²/(mPa·s); is the pseudo pressure of the matrix, MPa²/(mPa·s); is the time, h; is the pseudotime; is the gas viscosity, mPa·s; is the total compressibility coefficient, MPa⁻¹; and is the gas compressibility factor.

To describe the degree of stress sensitivity and its influence to fracture permeability, the concept of pseudopermeability modulus is defined as [19, 26]

Equation (2) can be further written as where is the permeability of fracture, mD; is the initial permeability of fracture, mD; is the initial pseudopressure, MPa²/(mPa·s); and is the pseudopermeability modulus of the fracture, (mPa·s)/MPa².

To establish and solve the model, a radial cylindrical system () is used to describe the flow of the fracture and matrix system. With consideration of the dual-porosity and dual-permeability flow behavior and stress-dependent fracture permeability, the governing differential equation of the complex system can be described as follows:

For the fracture system,

For the matrix system,

For the initial condition,

In the inner boundary condition,

The top and bottom boundaries are assumed to be closed and given by

The lateral boundary condition is assumed to be infinite and expressed as where is the total compressibility of the fracture, MPa⁻¹; is the total compressibility of the fracture, MPa⁻¹; is the reservoir thickness, m; is the initial horizontal permeability of the fracture, mD; is the initial vertical permeability of the fracture, mD; is the horizontal permeability of the matrix, mD; is the vertical permeability of the matrix, mD; is pressure, MPa; is the pressure at standard condition, MPa; is the surface gas production rate, 10⁴ m³/d; is the radial distance, m; is temperature, K; is the temperature at standard condition, K; is the vertical distance from the bottom, m; is the vertical distance of the horizontal well from the bottom, m; is a variable in the direction, m; is the porosity of fracture; is the porosity of the matrix; and is the geometric shape factor of matrix block, m⁻².

To make the equations homogeneous, some dimensionless variables are defined and tabulated in Table 1. Taking the dimensionless variables into equations (4)–(9), one can obtain the dimensionless differential equations.

For the fracture system,

For the matrix system,

In the initial condition,

In the inner boundary condition,

In the outer boundary conditions,

2.3. Solution to the Mathematical Model

It should be noted that equations (10) and (13) are strongly nonlinear with the consideration of the stress sensitivity of fracture permeability. This is because the stress-dependent fracture permeability is a function of the pseudopressure of the fracture. However, the pressure of fracture is an unknown parameter. Therefore, the mathematical model cannot be solved analytically. In this work, the Pedrosa [31] variable substitution and regular perturbation method are firstly deployed to alleviate the nonlinearity. Then, the Laplace transformation and separation of variables are adopted to address the linearized model. Thus, the model can be solved in the Laplace space and the Stehfest and Harald [32] numerical inversion is used to calculate the pressure in real space.

2.3.1. Linearization of the Flow Equation

To linearize the mathematical model, the Pedrosa transformation is employed in this section and given by where is an intermediate variable called the perturbation deformation function.

After the Pedrosa transformation, equations (10)–(16) can be rewritten as

According to the regular perturbation theory, equations (18)–(24) can be simplified as

2.3.2. Solution of the Proposed Model

To derive the analytical solution of the model, the mathematical model is translated into the Laplace domain with respect to : where is the Laplace transform variable.

With equations (32) and (33) and taking the Laplace transform of equations (25)–(31), one can obtain the dimensionless mathematical model in the Laplace space:

For the fracture system,

For the matrix system, where , , and .

In the inner boundary condition,

In the outer boundary conditions,

This section uses the separation of variables to solve the dual-porosity and dual-permeability modeling of a horizontal well in a naturally fractured reservoir. With the separation of variables, the dimensionless pseudopressure in the Laplace space can be separated by [25]

Substituting equation (40) into equations (34) and (35), one can obtain

According to equation (41), the fluid flow in the horizontal direction can be written as

And the fluid flow in the vertical direction is

Without consideration of the flow in the -direction and taking equation (42) into equations (34) and (35), we have

Under the infinite external boundary of side, the solutions of equations (45) and (46) can be expressed by

Substituting equations (47) and (48) into (45) and (46), we have

Because the modeling must have solutions, the coefficients and cannot be zero, so the term in equation (50) can be given by

With equation (49)–(51), the general solutions of equations (45) and (46) can be expressed by

Considering the fluid flow in the -direction, the general solutions of equations (52) and (53) can be given as

Combined with the boundary conditions, the terms and in equation (56) are

The general solution of equation (44) can be expressed by

Substituting equation (60) into equations (37) and (38), we have so the solution in the vertical direction is

Combined with equations (56), (57), and (62), the pressure solution for a three-dimensional volumetric source can be obtained and expressed by

According to the superposition principle, the dimensionless pressure in the Laplace space for constant-rate production can be obtained by integrating equation (63) along with the direction of horizontal wellbore.

Equation (64) is the dimensionless pseudopressure solution expression considering the skin effect. The wellbore storage effect can be incorporated in the abovementioned solution using Duhamel’s principle [25], and the bottom-hole pressure solution is

in equation (65) is the dimensionless pressure solution in the Laplace space, and in real space, the dimensionless pressure can be obtained using Stehfest numerical inversion [32]. After that, utilizing the inverse transformation of equation (66), the bottom-hole pressure response for a horizontal well incorporating the stress-dependent permeability of fracture system can be obtained:

3. New-Type Curves

In this work, a novel semi-analytical model is presented to examine the combined effects of dual-permeability flow behavior and stress-dependent fracture permeability on the transient pressure response of a horizontal well in a naturally fractured gas reservoir. A series of new transient pressure dynamic curves are drawn to observe different flow regimes based on the solution of the proposed model. The differences between the proposed model and traditional models are discussed, and the effects of some critical parameters on transient pressure response are also analyzed with the proposed model. Furthermore, it provides an efficient method for field engineers and related research and further to interpret pressure signals with accurate matching results and more reasonable interpreted parameters. After that, the effects of some critical parameters, including the dimensionless permeability modulus, the storability ratio of fracture, interporosity flow factor of matrix system into fracture system, and skin factor, on the characteristics of the type curves are examined and analyzed.

3.1. Flow Regime Identification

Figure 3 depicts the standard-type curves of wellbore pressure responses for a horizontal well located at a naturally fractured gas reservoir considering dual-permeability flow and stress-dependent fracture permeability. Basic data used to generate the type curves are listed in Table 2. An entire transient flow process is clearly shown, and the following six main flow stages can be recognized: (1)Pure wellbore storage stage: this period is characterized by a slope of 1 on the pressure and pressure derivative curves and governed by the wellbore storage coefficient, (2)Skin effect transition stage: the pressure derivative curve exhibits like a “hump.” The peak of the “hump” is dominated by the skin factor, (3)Early radial flow stage: this period is present whenever the wellbore storage coefficient, , and the horizontal wellbore length, , are suitable. During this period, the pressure derivative curve is a horizontal line with a value of “”(4)Early linear flow stage: it can be identified by a half-slope trend on the pressure derivative curve. During this period, gas flows linearly from the formation to natural fractures(5)Interporosity flow stage from the matrix system to the fracture system: this process is characterized by a concave in the pressure derivative curve. The concave shape is controlled by the storability ratio of fracture, , and the interporosity flow coefficient from the matrix system to the fracture system, (6)External boundary response stage: during this period, the pressures in the matrix and fracture systems reach a dynamic balance state. This period is marked by a slope of 0.5 on the pressure derivative curve without considering stress-dependent fracture permeability. However, the pressure derivative curve is no longer a horizontal line with a value of “0.5” but exhibits an upward tendency due to the effect of stress-dependent fracture permeability

3.2. Comparisons with the Traditional Models

So far, few attempts have been made to quantify the combined effects of dual-permeability flow behavior and stress-dependent fracture permeability on the transient pressure behavior of a horizontal well in fractured gas reservoirs. The main difference between the proposed and traditional models [25, 26] is that the dual-permeability flow behavior and stress-dependent fracture permeability are all incorporated in the new model. In this section, we simultaneously simulated the pressure response of both the proposed and conventional models using the same group of formation and well parameters in Table 2.

Figure 4 shows the comparison results of transient responses for the new model with the solution presented by Nie et al. [25]. The main difference between the two models is that stress-dependent fracture permeability is not considered in the model of Nie et al. [25]. As shown in Figure 4, fracture permeability stress sensitivity is found to significantly affect the middle stream flow period and the late-time pseudoradial flow period. There exist obvious differences during late-time pseudoradial flow period. The dimensionless pressure derivative curve exhibits a horizontal line with a value of 0.5 in the model of Nie et al. [25]; however, the derivative curve is no longer horizontal but bends upward in the new model. In addition, the location of the dimensionless pressure and pressure derivative curves during the middle stream flow period is higher in our model. This is because an additional pressure drop will be required to maintain a constant flow rate when the stress-sensitivity effect is taken into account.

We also compared the proposed model with the solution of a single-permeability model [26] based on the same reservoir properties and fracture parameters. The stress-dependent fracture permeability is considered both in two models. As can be seen in Figure 5, the combined effects of dual-permeability flow behavior and stress-dependent fracture permeability play a significant role in the early radial flow stage, the early linear flow stage, and the interporosity flow stage of the matrix system to the fracture system. Dimensionless pressure of the proposed model is lower than that of the single-permeability model [26] during the three flow periods. Besides, the concave in the pressure derivative curve of the dual-permeability model is shallower than that of single-permeability model. This is because the single-permeability model assumes that the fracture system is the only flow pathway directly connected with wellbore; however, the fluid supply from the matrix system to the wellbore is not considered. The dual-permeability flow behavior will accelerate energy supplement in the matrix during production compared with the single-permeability model.

3.3. Sensitivity Analysis

Based on the proposed model, the influences of stress sensitivity of fracture permeability, storativity ratio of the fracture, interporosity flow factor, and skin factor on pressure response are discussed. Except for the parameters analyzed, other parameters are the same and are shown in Table 2.

3.3.1. Effect of the Permeability Modulus of the Fracture

The stress-sensitivity effect can be determined with the dimensionless permeability modulus, . Figure 6 shows the pressure and pressure derivative curves for , 0.06, and 0.08. As stated, with increasing the value of , the slope of the derivative curves increases during the intermediate and late time period. This is because the pressure drop increases as the dimensionless permeability modulus increases and fluid flow will be difficult and more gas is left in the reservoir. The stress sensitivity reflects the damage of permeability, and a larger dimensionless permeability modulus will increase the damage of permeability. Consequently, the permeability stress sensitivity of the fracture decreases the cumulative production. So, it is believed that reasonable producing pressure differential is excellent for reducing the negative effect of stress sensitivity on gas productivity in the development of fractured gas reservoirs.

3.3.2. Effect of the Storability Ratio of the Fracture

The effect of the storability ratio of the fracture, , on transient behavior is shown in Figure 7. As shown in Figure 7, the storability ratio of the fracture not only determines the duration and the depth of the concave but also has a significant effect on the early flow regimes (early radial and early linear flow stage). It can be clearly observed that the larger the is, the deeper and wider the concave in dimensionless pressure derivative curve. In addition, the dimensionless pressure curve becomes higher with the increase of the storability ratio of the fracture. This is because the storativity ratio of fracture reflects the relative capacity of fluid stored in the fracture system; a smaller storability ratio of the fracture is the response of relative abundant reserves in the matrix system. The pressure drop should increase to maintain the constant production rate when increasing the storability ratio of the fracture.

3.3.3. Effect of the Interporosity Flow Factor

The effect of the interporosity flow factor of the matrix system to the fracture system, , on pressure response is shown in Figure 8. According to the definition of the interporosity flow factor, the represents the starting time of the flow exchange from the matrix system to the fracture system. The larger the is, the earlier the time of the interporosity flow period is. Besides, the late-time pseudoradial flow period would be masked if the is large enough.

3.3.4. Effect of the Skin Factor

The effect of the skin factor, , on pressure response is shown in Figure 9. As shown in Figure 9, skin factor plays a significant role in the early transient flow period. A larger skin factor leads to a higher location of dimensionless pressure and pressure derivative curves. This is because the incremental value of the skin factor results in the increasing additional filtration resistance and the skin effect transition period will last longer. The larger the skin factor is, the slower the pressure wave propagates to the external boundary and the larger the pressure drop is.

4. Conclusions

This paper provided a semianalytical model to investigate the combined effects of dual-permeability flow behavior and stress-dependent fracture permeability on the transient pressure response for a horizontal well in a naturally fractured gas reservoir. The main conclusions of this work are as follows: (i)A horizontal production well in a naturally fractured gas reservoir with consideration of stress sensitivity effect may exhibit six flow stages: pure wellbore storage stage, skin effect transition stage, early radial flow stage, early linear flow stage, interporosity flow stage, and external boundary response stage(ii)The stress-dependent fracture permeability imposes effects on the intermediate flow period and the late-time pseudoradial flow period; the existence of dual-permeability flow behavior can make the stress sensitivity effect more significant(iii)The storability ratio of the fracture, , mainly affects the duration and the depth of the concave, and a larger leads to a deeper and wider concave in dimensionless pressure derivative curve. In addition, the storability ratio plays a significant role in the early radial flow stage and early linear flow stage. The interporosity flow factor, , mainly affects the starting time of the flow exchange from the matrix to the fracture. The larger the value of , the earlier the occurrence of the interporosity flow period(iv)The proposed model is suitable for various naturally fractured gas reservoirs and can interpret pressure signals with accurate matching results and more reasonable interpreted parameters

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors acknowledge that this study was partially funded by the 2021 Innovation and Entrepreneurship Training Projects for College Students in Shanxi Province (no. S202110719144).

References

Y. Wu, L. Cheng, S. Fang, S. Huang, and P. Jia, “A green element method-based discrete fracture model for simulation of the transient flow in heterogeneous fractured porous media,” Advances in Water Resources, vol. 136, p. 103489, 2020.
View at: Publisher Site | Google Scholar
P. Jia, L. Cheng, S. Huang, and Y. Wu, “A semi-analytical model for the flow behavior of naturally fractured formations with multi-scale fracture networks,” Journal of Hydrology, vol. 537, pp. 208–220, 2016.
View at: Publisher Site | Google Scholar
X. Tian, L. Cheng, R. Cao et al., “A new approach to calculate permeability stress sensitivity in tight sandstone oil reservoirs considering micro-pore-throat structure,” Journal of Petroleum Science and Engineering, vol. 133, pp. 576–588, 2015.
View at: Publisher Site | Google Scholar
G. E. Barenblatt, I. P. Zheltov, and I. N. Kochina, “Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks [strata],” Journal of Applied Mathematics and Mechanics, vol. 24, no. 5, pp. 1286–1303, 1960.
View at: Publisher Site | Google Scholar
J. E. Warren and P. J. Root, “The behavior of naturally fractured reservoirs,” SPE Journal, vol. 3, no. 3, pp. 245–255, 1963.
View at: Publisher Site | Google Scholar
H. Kazemi, M. S. Seth, and G. W. Thomas, “The interpretation of interference tests in naturally fractured reservoirs with uniform fracture distribution,” SPE Journal, vol. 9, no. 4, pp. 463–472, 1969.
View at: Publisher Site | Google Scholar
A. O. De Swaan, “Analytic solutions for determining naturally fractured reservoir properties by well testing,” SPE Journal, vol. 16, no. 3, pp. 117–122, 1976.
View at: Publisher Site | Google Scholar
R. Raghavan and C. U. Ohaeri, “Unsteady flow to a well produced at constant pressure in a fractured reservoir,” in SPE California regional meeting, Bakersfield, California, 1981, SPE paper no. SPE-9902-MS.
View at: Publisher Site | Google Scholar
K. V. Serra, A. C. Reynolds, and R. Raghavan, “New pressure transient analysis methods for naturally fractured reservoirs(includes associated papers 12940 and 13014),” Journal of Petroleum Technology, vol. 35, no. 12, pp. 2271–2283, 1983.
View at: Publisher Site | Google Scholar
Y. Jalali and I. Ershaghi, “Pressure transient analysis of heterogeneous naturally fractured reservoirs,” in SPE California regional meeting, Ventura, California, 1987, SPE paper no. SPE 16341-MS.
View at: Publisher Site | Google Scholar
Y. Wu and K. Pruess, “A multiple-porosity method for simulation of naturally fractured petroleum reservoirs,” SPE Reservoir Engineering, vol. 3, no. 1, pp. 327–336, 1988.
View at: Publisher Site | Google Scholar
T. D. Bui, D. D. Mamora, and W. J. Lee, “Transient pressure analysis for partially penetrating wells in naturally fractured reservoirs,” in SPE rocky mountain regional/low permeability reservoirs symposium and exhibition, Denver, Colorado, 2000, SPE paper no. SPE 60289-MS.
View at: Publisher Site | Google Scholar
Y. Wu, H. Liu, and G. S. Bodvarsson, “A triple-continuum approach for modeling flow and transport processes in fractured rock,” Journal of Contaminant Hydrology, vol. 73, no. 1-4, pp. 145–179, 2004.
View at: Publisher Site | Google Scholar
F. Kuchuk, D. Biryukov, and T. Fitzpatrick, “Fractured-reservoir modeling and interpretation,” SPE Journal, vol. 20, no. 5, pp. 983–1004, 2015.
View at: Publisher Site | Google Scholar
S. Wang, L. Cheng, S. Huang et al., “A semi-analytical method for modeling two-phase flow behavior in fractured carbonate oil reservoirs,” Journal of Energy Resources Technology, vol. 141, no. 7, pp. 072902.1–072902.11, 2019.
View at: Publisher Site | Google Scholar
B. X. Hu and H. Huang, “Stochastic analysis of reactive solute transport in heterogeneous, fractured porous media: a dual-permeability approach,” Transport in Porous Media, vol. 48, no. 1, pp. 1–39, 2002.
View at: Publisher Site | Google Scholar
T. M. Al-Shaalan, L. S. K. Fung, and A. H. Dogru, “A scalable massively parallel dual-porosity dual-permeability simulator for fractured reservoirs with super-k permeability,” in SPE annual technical conference and exhibition, Denver, Colorado, 2003, SPE paper no. SPE 84371-MS.
View at: Publisher Site | Google Scholar
A. P. G. Van Heel, P. M. Boerrigter, and J. J. Van Dorp, “Thermal and hydraulic matrix-fracture interaction in dual-permeability simulation,” SPE Reservoir Evaluation and Engineering, vol. 11, no. 4, pp. 735–749, 2008.
View at: Publisher Site | Google Scholar
L. Jiang, T. Liu, and D. Yang, “A semianalytical model for predicting transient pressure behavior of a hydraulically fractured horizontal well in a naturally fractured reservoir with non-darcy flow and stress-sensitive permeability effects,” SPE Journal, vol. 24, no. 3, pp. 1322–1341, 2019.
View at: Publisher Site | Google Scholar
P. Yu, D. Dempsey, and R. Archer, “Techno-economic feasibility of enhanced geothermal systems (EGS) with partially bridging multi-stage fractures for district heating applications,” Energy Conversion and Management, vol. 257, p. 115405, 2022.
View at: Publisher Site | Google Scholar
V. Celis, R. Silva, M. Ramones, J. Guerra, and G. da Prat, “A new model for pressure transient analysis in stress sensitive naturally fractured reservoirs,” SPE Advanced Technology, vol. 2, no. 1, pp. 126–135, 1994.
View at: Publisher Site | Google Scholar
E. Ozkan, R. Raghavan, and S. D. Joshi, “Horizontal-well pressure analysis,” SPE Formation Evaluation, vol. 4, no. 4, pp. 567–575, 1989.
View at: Publisher Site | Google Scholar
E. J. Kuchuk, P. A. Goode, D. J. Wilkinson, and R. K. M. Thambynayagam, “Pressure-transient behavior of horizontal wells with and without gas cap or aquifer,” SPE Formation Evaluation, vol. 6, no. 1, pp. 86–94, 1991.
View at: Publisher Site | Google Scholar
A. Hashemi, L. M. Nicolas, and A. C. Gringarten, “Well-test analysis of horizontal wells in gas/condensate reservoirs,” SPE Reservoir Evaluation and Engineering, vol. 9, no. 1, pp. 86–99, 2006.
View at: Publisher Site | Google Scholar
R. Nie, Y. Meng, Y. Jia, F. X. Zhang, X. T. Yang, and X. N. Niu, “Dual porosity and dual permeability modeling of horizontal well in naturally fractured reservoir,” Transport in Porous Media, vol. 92, no. 1, pp. 213–235, 2012.
View at: Publisher Site | Google Scholar
X. Li, L. Cao, C. Luo, B. Zhang, J. Q. Zhang, and X. H. Tan, “Characteristics of transient production rate performance of horizontal well in fractured tight gas reservoirs with stress-sensitivity effect,” Journal of Petroleum Science and Engineering, vol. 158, pp. 92–106, 2017.
View at: Publisher Site | Google Scholar
S. Wang, L. Cheng, Y. Xue et al., “A semi-analytical method for simulating two-phase flow performance of horizontal volatile oil wells in fractured carbonate reservoirs,” Energies, vol. 11, p. 10, 2018.
View at: Publisher Site | Google Scholar
Z. Sun, B. Huang, Y. Li, H. Lin, S. Shi, and W. Yu, “Nanoconfined methane flow behavior through realistic organic shale matrix under displacement pressure: a molecular simulation investigation,” Journal of Petroleum Exploration and Production Technologies, vol. 12, no. 4, 2022.
View at: Publisher Site | Google Scholar
P. Yu, D. Dempsey, and R. Archer, “A three-dimensional coupled thermo-hydro-mechanical numerical model with partially bridging multi-stage contact fractures in horizontal-well enhanced geothermal system,” International Journal of Rock Mechanics and Mining Sciences, vol. 143, article 104787, 2021.
View at: Publisher Site | Google Scholar
Z. Sun, S. Wang, H. Xiong, K. Wu, and J. Shi, “Optimal nanocone geometry for water flow,” AIChE Journal, vol. 68, no. 3, article e17543, 2022.
View at: Publisher Site | Google Scholar
O. A. Pedrosa, “Pressure transient response in stress-sensitive formations,” in SPE California regional meeting, Oakland, California, 1986, SPE paper no. SPE-15115-MS.
View at: Publisher Site | Google Scholar
Stehfest and Harald, “Algorithm 368: numerical inversion of Laplace transforms [D5],” Communications of the ACM, vol. 13, no. 1, pp. 47–49, 1970.
View at: Publisher Site | Google Scholar
F. Kuchuk and Biryukov, “Pressure-transient behavior of continuously and discretely fractured reservoirs,” SPE Reservoir Evaluation and Engineering, vol. 17, no. 1, pp. 82–97, 2014.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Guiqin Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies