Next Article in Journal
Numerical Simulation of Diffusion Regularity and Parameter Optimization of Shaft Grouting Slurry
Next Article in Special Issue
An Unsupervised Clustering Method for Selection of the Fracturing Stage Design Based on the Gaussian Mixture Model
Previous Article in Journal
Nitric Oxide-Releasing NO–Curcumin Hybrid Inhibits Colon Cancer Cell Proliferation and Induces Cell Death In Vitro
Previous Article in Special Issue
Investigation into Hydraulic Fracture Propagation Behavior during Temporary Plugging and Diverting Fracturing in Coal Seam
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Shale Gas Productivity Prediction Model Considering Time-Dependent Fracture Conductivity

1
School of Petroleum Engineering, Yangtezu University, Wuhan 430100, China
2
The Multiphase Flow Laboratory of Gas Lift Innovation Center, China National Petroleum Corporation, Wuhan 430100, China
3
Key Laboratory of Drilling and Production Engineering for Oil and Gas, Wuhan 430100, China
4
Engineering Institute of Engineering Technology, Xinjiang Oilfield Company, Karamay 834000, China
5
Research Institute of Exploration and Development, Petrochina Changqing Oilfield Company, Xi’an 710016, China
*
Author to whom correspondence should be addressed.
Processes 2022, 10(5), 801; https://doi.org/10.3390/pr10050801
Submission received: 23 March 2022 / Revised: 11 April 2022 / Accepted: 17 April 2022 / Published: 19 April 2022

Abstract

:
Conventional shale gas productivity prediction techniques consider fracture conductivity to be a fixed value, but in actual production processes, conductivity changes with time. Therefore, this paper proposed a capacity prediction method that considers time-dependent conductivity and validates its accuracy using commercial simulators. First, relevant parameters were obtained by fitting the improved long-term conductivity test, and then the shale gas seepage model was established using the EDFM method. The laboratory test results showed that the order of significance affecting the conductivity retention rate was fracturing fluid viscosity > sand concentration > fracturing fluid retention time; the calculation results of the production prediction model show that the flow and the pressure curves that corresponded to constant conductivity and variable conductivity were to some extent different. In the presence of complex fractures and natural fractures, the increase in the variable conductivity production curve was smaller than that of the constant conductivity production curve. This study provides some guidance for field production.

1. Introduction

Shale gas is an important unconventional natural gas resource [1], and the multi-stage fracturing of horizontal wells is a key technology used to increase the level of production in the development of shale gas fields [2]. Due to the difference between shale reservoirs and conventional reservoirs [3], the question of how to accurately evaluate and predict the productivity of multi-stage fracturing wells remains a research hotspot. Traditional gas well dynamic analyses and productivity evaluations often focus on the influence of the physical properties of reservoirs, and they assign a certain value to fracture conductivity [4,5]. However, in the actual production process, the diversion capacity constantly changes with time, and its impact on the production capacity cannot be ignored.
There are many factors that affect the conductivity of hydraulic fractures, including the properties of the fracturing fluid [6], the proppant [7] itself, and its interaction with the reservoir rock [8]. As the production progresses, the proppant is broken and embedded [9], the fractures are gradually closed, and the conductivity is gradually reduced [10]. Studies have shown that proppant embedment is more severe when the rock is exposed to a water-softened state, potentially losing up to 88% of conductivity [11]. Obadare et al. proposed a fracture conductivity prediction model that considered dynamic data as they argued that changes in conductivity are related to the degree of breakage within fracturing fluid gel [12]. Yu et al. studied the change law concerning the conductivity of propped fractures alongside the proppant loading concentration; the stress sensitivity of conductivity was weaker when the proppant loading concentration was higher [13]. Wang et al. believed that the placement form of a proppant had a great influence on conductivity [14].
The mainstream yield forecasting methods mainly include analytical methods (including semi-analytical methods), numerical simulation methods, and some machine learning methods [15,16,17]. The models of analytical methods are too simplified, and they do not have high calculation accuracy. The machine learning methods are similar to the black box model, and their interpretability is poor. Numerical simulation methods can simulate complex flow states, and they are widely used in commercial simulators, but at present time-dependent flow conductivity has not been paid much attention.
In response to the above problems, here, we mainly carried out the following two aspects of work: (1) a conductivity test was conducted to study the influence mechanism of different factors on conductivity, and we proposed methods to improve the formula concerning time-dependent conductivity; and (2) based on the EDFM, a shale gas productivity prediction model considering time-dependent conductivity was established, and the model was extended and applied.

2. Time-Dependent Effect of Fracture Conductivity

With gas production, the net stress around the fracture increases during the production process; phenomena such as proppant fragmentation, cuttings plugging, and proppant embedding cause the fracture to gradually close; and the fracture conductivity decreases with time. This effect is called time-dependent fracture conductivity [18].

2.1. Common Formulas and Problems

There are many factors that affect fracture conductivity, and it is difficult to establish a mathematical model to characterize it through theoretical research. At present, the conductivity function relationship is obtained by fitting the experimental data mainly through laboratory tests. At present, the exponential relationship is commonly used to characterize the change in fracture conductivity with time [19]:
F c t = F c 0 ( 1 + η e t / c η )
where F c t is the conductivity of the fracture at time t; η and c is the fitting coefficient; F c 0 is the conductivity of the fracture at the initial time. However, Equation (1) is inconvenient to use for the following reasons:
(1)
An inconsistent experiment process:
Figure 1 shows the long-term conductivity test results of shale reservoirs [20,21]. Figure 1a shows the results of a conventional long-term conductivity test in which F c 0 was equal to the initial conductivity at the start of the test. The gas was injected throughout the test, and the damage to the conductivity caused by the immersion of the fracturing fluid was not considered. In fact, due to the influence of factors such as fracturing and shut-in, the fracture rate of the proppant changes after being soaked in fracturing fluid, and this change affects the change law of conductivity with time.
Figure 1b shows the results of the long-term conductivity test considering the influence of the liquid. The test process was divided into three stages: I. Gas was injected during the test to measure the conductivity at the initial moment; II. Fracturing fluid was injected to simulate the well shut-in, and the process of soaking the proppant in fracturing fluid was carried out; and III. Gas was injected again, and the initial moment of conductivity after the fracturing fluid was soaked was measured and used as F c 0 in Equation (1). Additionally, we used the process of Figure 1b to conduct experiments in this paper;
(2)
Variable experiment parameters:
The practicability of Equation (1) is poor. When the test conditions (i.e., field construction parameters) are changed, the test needs to be carried out again to determine the values of η and c in Equation (1) under the new test conditions, which undoubtedly increases the time it takes to obtain results.

2.2. Long-Term Conductivity Test Experiment

In view of the above problems, we obtained an improved method for the equation through a long-term conductivity test experiment. The experimental design employed a response surface method. The basic idea behind this methodology is to vary multiple parameters at the same so that maximum inference can be attained with minimum cost, and the method can quantify the relationship between the experimental results and the experimental parameters.
Notably, closure pressure, proppant type, and proppant particle size can be considered constant for a given reservoir condition. Combined with the previous experimental results, we chose the sand concentration, fracturing fluid viscosity, and fracturing fluid residence time as variables. The experimental scheme is shown in Table 1.
During the experiment, the temperature was kept at a constant of 25°, the closing pressure was 50 MPa, the liquid flow rate was 5 mL/min, and the gas flow rate was 300 mL/min.
In addition, the conductivity retention rate was defined as the ratio of the conductivity before and after fracturing fluid soaking:
R f g = F c 2 F c 1
The experimental results are shown in Table 2.
The experimental data were processed using Design Expert software, and the regression equation for the prediction of the retention rate of conductivity was obtained:
R f g = 41.83 + 7.09 A 10.96 B 2 C + 0.2316 A B + 0.0238 A C + 0.2610 B C + 0.5755 A 2 + 1.34 B 2 0.6028 C 2
At the same time, after obtaining A and B in Formula (1) based on the response surface test results, nonlinear polynomial fitting was performed on these values using the test parameters. Finally, the relationship between A and B and sand concentration, fracturing fluid viscosity, and fracturing fluid residence time was obtained:
η = 5.868 0.244 A 0.141 B 0.545 C + 0.061 A B + 0.028 A C 0.022 B C 0.027 A 2 0.011 B 2 + 0.014 C 2
c = 0.036 + 0.021 A 0.003 B 0.014 C + 0.0004 A B + 0.0002 A C 0.0002 B C 0.002 A 2 + 0.0005 B 2 + 0.0005 C 2
where A is the sand concentration, B is the fracturing fluid viscosity, and C is the fracturing fluid retention time.
In addition, according to the analysis results of the response surface model, the influences of the three factors on the response value were ranked as fracturing fluid viscosity > sand concentration > fracturing fluid residence time. The coefficient of determination, R2 = 0.9471, of the regression model indicates that the model fit well and the model could be used for prediction.

2.3. Amendment to the Formula of Conductivity

As mentioned above, there are two main aspects to the correction of Equation (1); firstly, the initial conductivity in Equation (1) was corrected using Equation (3) and then using Equations (4) and (5) to modify η   c , and the final formula for the change in conductivity with time can be written as:
F c t = F c 0 R f g ( 1 + η e t / c η )
where R f g is the retention rate of conductivity; η c is the coefficient corrected using Equations (4) and (5).

3. Shale Gas Physicomathematical Statements

The high cost of developing unconventional oil and gas resources also means that higher-quality predictive models of production capacity are required to design and manage production processes. The complex migration mechanism of shale gas reservoirs and the development of natural fractures make productivity prediction a challenge.
This section provides a theoretical description of the established shale gas seepage model. Referring to the research results of Wang B et al., the single-phase gas seepage equation in the matrix and fractures can be written as [22]:
t ( b ϕ m ) + t [ ( 1 ϕ m ) ρ 0 q g ] = ( b k m μ g p m ) + b q b q m f
t ( b ϕ f ) = ( b k f μ g p f ) + b q + b q m f
where b is the inverse formation volume factor, equal to ρ g / ρ 0 ; ρ g is the density of the shale gas in the reservoir, g/cm3; ϕ m ϕ f are the matrix and the fracture porosity, respectively; q g is the adsorption quality for shale gas, cm3/s; k m k f are the corrected matrix and the fracture permeability, respectively, k = f 1 f 2 f 3 k 0 , D; μ g is the viscosity of the shale gas in the reservoir, mPa·s; q is the source and sink, cm3/s; q m f is the leakage term between matrix and fractures, cm3/s.
In addition, to build the model, the following equations are required:

3.1. Gas Adsorption

In terms of organic matter, kerogen, clay minerals in shale, and a large amount of adsorbed gas on the pore surface, there are various isotherm adsorption curves with different forms to describe the isotherm adsorption process of different adsorption systems under different assumptions [23]. The Langmuir isotherm adsorption equation is usually used to describe the adsorption process of shale gas:
V = V L p p L + p
The mass of gas adsorbed per unit rock volume is:
q g = V L p p L + p ρ s ρ s c
where ρ s is the density of shale under formation pressure, kg/m3; V L is the Langmuir volume, m3/kg; p L is the Langmuir pressure, MPa.

3.2. Gas Slip and Spread

Due to the influence of slippage and diffusion, the flow of shale gas in matrix pores no longer follows Darcy’s law. Florence et al. proposed a correction factor for shale gas permeability that can be used to obtain shale gas permeability when slip and Kundsen diffusion are considered [24]:
f 1 = ( 1 + α k K n ) ( 1 + 4 K n 1 + K n )
where α k is the rarefaction coefficient parameter:
α k = 128 15 π 2 tan 1 ( 4 K n 0.4 )
Additionally, K n is the Kundsen constant:
K n = λ d

3.3. Gas Non-Darcy Flow

When the flow velocity in a fracture is too high, there is no longer a linear relationship between the pressure difference and the flow velocity, that is, the high-speed non-Darcy effect occurs. The permeability is corrected using the related concept proposed by Barree et al. [25]:
f 2 = 2 1 + 1 + 4 ρ β ( k 0 μ ) 2 | p |
where β is the empirical Forchheimer coefficient:
β = 3.2808 1.485 × 10 9 ( k 0 × 10 15 ) 1.021

3.4. Time-Dependent Conductivity

In the first section of the paper, the laboratory test of the conductivity test that was carried out is described, and the variation law of the conductivity with time that was obtained is described. Assuming that the fracture width is constant, Equation (6) was used to characterize the change in fracture conductivity:
f 3 = R f g ( 1 + η e t / c η )

4. The EDFM Numerical Model

Lee and Moinfar et al. proposed the embedded discrete fracture model (EDFM), which mathematically treats fractures as well-like sources and “embeds” them into the matrix, avoiding local mesh refinement. Additionally, it can better adapt to the situation of complex fracture distribution [26].

4.1. Numerical Discrete

Based on the mathematical model of seepage established above, it is discretized using the finite volume method. The matrix Equation (7) can be written as:
ϕ m V m Δ t ( b n + 1 b n ) + ( 1 ϕ m ) V m q g ρ 0 Δ t ( p m n + 1 p m n ) = ( V m k m μ g n + 1 b n + 1 p m n + 1 ) + V m b n + 1 q n + 1 V m b n + 1 q m f n + 1
Equation (8) can be written as:
ϕ f V f Δ t ( b n + 1 b n ) = ( V f k f μ g n + 1 b n + 1 p f n + 1 ) + V f b n + 1 q n + 1 + V f b n + 1 q m f n + 1
For the meaning of the parameter of Equations (17) and (18), refer to Equations (7) and (8).

4.2. Non-Neighbor Connections

For the two-dimensional reservoir model, the fluid in the fracture flowed along the tangential direction of the fracture wall, that is, a one-dimensional flow, and the fracture grid was a line element. Flow exchange occurred between the matrix grid block and the fracture unit in the form of:

4.2.1. Fracture–Matrix NNC

For the fracture–matrix intersection shown in Figure 2 below, the wellhead treatment method in reservoir numerical simulation was used to calculate the channel flow between the fracture and the matrix grid block intersecting with the fracture. Peaceman established a wellhead model for vertical wells, and Lee et al. introduced a well source-like flow coefficient to calculate leakage flow [27]:
For two-dimensional fluid seepage, the channeling flow between the fracture grid and the matrix grid can be expressed by the following formula:
q m f = T I k ( p m n + 1 p f n + 1 )
where T I k is flow coefficient:
T I k = 2 k m l k h f μ d ¯
where k m is the matrix permeability, D; l k is the step size of the crack mesh, m; h f is the height of fracture, m; d ¯ is the average distance from the matrix grid block to the fracture, m.
d ¯ = 1 A l x y d A
where A is the area of the matrix grid, m2; l x y is the distance from area element to fracture in matrix mesh, m.

4.2.2. Fracture–Fracture NNC

As shown in Figure 3, multiple fractures intersected for flow exchange, and the volume at the node was much smaller than the grid size of other fractures. If the node also participates in the calculation of the governing equation, the calculation speed will be greatly reduced. Therefore, the conductivity factor between fracture elements was used to determine the flow exchange at the fracture node.
When the fracture height was constant, the flow exchange of intersecting fractures at nodes satisfied the material balance equation, namely:
i = 1 n w f i h f V i o = 0
where w f i is the width of fracture, m; V i o is the velocity of fracture, m/s.
According to Darcy’s law, the flow of fluid flowing into node O through the i-th fracture can be obtained:
Q i o = T i o ( p f o p f i )
where T i o is the conductivity factor:
T i o = 2 k f i w f i h f μ l i
where k f i is the permeability of fracture, D; l i is the length of fracture grid, m.
Finally, the conductivity factor between the intersecting fractures was determined:
T i j = T j o T i o n = 1 m T n o
where T i j is the conductivity factor between fractures; m is the total number of intersecting fractures.

4.2.3. Fracture–Well NNC

According to the method proposed by Peaceman et al., the flow conductivity factor between the fracture unit and the well can be expressed as:
W I = 2 π k f w f ln ( r e / r w )
r e = 0.14 ( Δ x ) 2 + ( Δ y ) 2
where r e is the equivalent well radius, m; r w is the well radius, m; Δ x Δ y are the lengths of the grid in the x and y directions, m, respectively.

4.3. Analysis of Model Calculation Results

4.3.1. Validation

In order to verify the accuracy of the model, the calculation results of the model were compared with the Uncong (Unconventional Oil and Gas Simulator, Ennosoft) commercial simulator. Uncong is a reservoir numerical simulation software. It also uses embedded discrete fractures to simulate fracturing fractures and natural fractures. Physical processes such as Kundsen diffusion and Langmuir adsorption can be considered. The calculation parameters are shown in Table 3.
Figure 4 shows the cumulative production calculation for one year of production with the model set to constant conductivity. The error between the model production calculation results and the Uncong calculation results was 3.11%. In the early stage, the model calculation results showed good consistency with the Uncong calculation. As production progressed, the model calculations were slightly lower than the Uncong calculation. Errors may have occurred due to differences in solution methods and in differences in some mechanism models.
In addition, in order to facilitate the comparison of formation pressure changes (Figure 5), the authors set formation pressure contours in the Uncong calculation results. It can be seen from the settlement results that the two maintained good consistency.

4.3.2. Analysis of Flow and Pressure Changes in Fractures

The formation pressure was set to 30 MPa, the bottom hole pressure was 15 MPa, the initial fracture permeability was 10D, η = 0.9, c = 0.008, the model size was 300 × 300, the number of matrix grids was 100 × 100, and the production time was 1 year. The variation law of fracture flow and the pressure in different production times were studied under constant conductivity and variable conductivity conditions. Figure 6 shows the formation pressure changes corresponding to different production times.
Figure 7 shows the pressure distribution at the near-wellbore location corresponding to different times. In both cases, the pressure near the wellbore exhibited a distinct “funnel” shape. In the early stage of production (50d), the “funnel” curves of the two cases almost overlapped, but the “funnel” that corresponded to the constant conductivity pressure curve was deeper, indicating that the flow into the wellbore was larger under the constant conductivity, so the drop in pressure was larger [28].
With the progress of production, the “funnel” of the variable conductivity pressure curve gradually separated from the “funnel” of the constant conductivity pressure curve. This is because the fracture conductivity gradually decreased under the condition of variable conductivity, and less fluid flowed into the fracture. The pressure around the fracture was higher, so the “funnel” of the variable conductivity pressure curve was thinner than the “funnel” of the constant conductivity pressure curve.
Figure 8 displays the flow distribution curve of the well point location. In the early stage of production, the flow rates at each well point in the two cases were quite different, and both appeared as a straight line. With the increase in time, the inter-fracture interference became stronger, and the fractures at both ends were less disturbed and had a larger oil drainage area. This meant that the production rate was higher, displaying an upward trend at both ends.
In addition, the connecting line of the middle well point was approximately a straight line in both cases, but the upturn of both ends of the flow curve of variable conductance was not obvious. This also means that, in actual production, if the change in conductivity over time is considered, the production of the fractures at both ends may be lower than the value expected by reservoir engineers.

4.3.3. Effects of Complex Fractures

Several branch fractures were added around the artificial fractures to simulate production under complex fracture morphology. The parameter settings were consistent with the previous section, the production time was one year, and the calculation results regarding pressure and cumulative output are shown in the figure below.
Figure 9 shows the cloud map of formation pressure under the conditions of constant conductivity and variable conductivity for one year. It can be seen from the figure that the range and the degree of the formation pressure drop was smaller when the production with variable conductivity was used. This is because the initial production of shale reservoirs mainly comes from artificial fractures. As the production progressed, the proppant was embedded and broken, resulting in the gradual closure of artificial fractures, and the decrease in conductivity led to an increase in the flow resistance of fluids, so the formation pressure dropped more slowly [29].
Figure 10 shows a comparison of one year of cumulative production under constant and variable conductivity conditions. It can be seen from the figure that the cumulative output that corresponded to the variable conductivity was smaller than the cumulative output that corresponded to the constant conductivity, and the final output differed by 19%, indicating that the time-varying effect of conductivity on the actual output cannot be ignored.

4.3.4. Effects of Natural Fractures

In order to consider the production situation under the condition of natural fractures, 50 natural fractures were randomly distributed in the reservoir, and the rest of the calculation parameters and grid division were consistent with the previous section.
Figure 11 displays the cloud map of formation pressure under the condition of natural fractures. Compared with the situation that did not consider natural fractures, natural fractures were shown to strengthen the connection between artificial fractures and the reservoir, which was equivalent to increasing the matrix permeability to a certain extent, so the pressure wave propagated faster in the reservoir and the drop in pressure was more obvious.
On the other hand, although the presence of natural fractures resulted in a faster drop in pressure, there was still a significant difference in the formation pressure distribution between the two cases (constant versus variable conductivity). This is consistent with the conclusion obtained above: as production progresses, the conductivity gradually decreases. The greater the decline rate or the faster the decline rate, the shorter the stable production period of the gas well and the earlier the decline in production.
Figure 12 shows the comparison of one year of cumulative production under the conditions of constant conductivity and variable conductivity. Compared with the previous section, the existence of natural fractures increased the production in both cases (constant conductivity vs. variable conductivity) by 14% and 13%, respectively. The cumulative yield that corresponded to constant conductivity and the cumulative yield that corresponded to variable conductivity differed by 21%.

5. Application Case

Using data from Well W1 in the China Shale Gas Field to verify the productivity prediction model, the depth of the well was 2531 m, the formation pressure was 35 MPa, the permeability was 2.08 × 10−4 mD, the porosity was 5%, the compressibility was 1.05 × 10−4 MPa−1, the fracture length was 120 m, the height was 40 m, and the width was 0.003 m. The fracturing parameters of this well were used in Equations (4)–(6) to calculate R f g = 72.01%, η = 0.4, c = 0.003. The production performance of the well was evaluated using the established productivity prediction model, considering the effect of variable conductivity. The comparison results are shown in Figure 13.
In the early stage of production, the decline in conductivity was the largest, so compared with the calculation results of the constant conductivity model, the calculation results of the variable conductivity model were most similar to the actual production data in the early stage. Although the conductivity decreased in the later stage of production, the magnitude was small, so the difference between the calculation results of the constant conductivity model and the variable conductivity model in the later stage of production was small.
Generally speaking, the error between the calculation result of the variable conductivity model and the actual situation was 18%, which means the variable conductivity model was more consistent with the actual situation, and the production capacity forecast needs to consider the change in conductivity over time.

6. Conclusions

(1) The order of significance of the test factors on the conductivity retention rate is: fracturing fluid viscosity > sand concentration > fracturing fluid retention time.
(2) According to the experimental data, the relationship between R f g , η , and c and the fracturing fluid viscosity, sand concentration, and fracturing fluid residence time was fitted. An improved method for the formula of the dependent conductivity was proposed, and this formula is more practical.
(3) A prediction model of shale gas production was established, and its accuracy was verified using the commercial simulator Uncong (the error of the calculation result was 3.11%). The analysis results of the flow change law in the fracture show that, compared with the constant conductivity flow curve, the upturn at both ends of the variable conductivity flow curve was not obvious; the analysis results of the pressure change law in the fracture show that the “funnel” of the pressure curve of variable conductivity under the same time was thinner and shorter.
(4) We investigated the impact on production performance in the presence of complex fractures and in the presence of natural fractures. In both cases, the cumulative production of the constant conductivity model and the variable conductivity model could be increased, but the cumulative production increase corresponding to the variable conductivity model was smaller than that of the constant conductivity model.
(5) The practical application results show that the production prediction model considering the time-dependent conductivity was closer to the production data, with an error of 18%. This shows that time-dependent conductivity cannot be ignored in the production forecast.

Author Contributions

Conceptualization, Y.P.; methodology, Y.P.; software, Y.P.; validation, Y.P.; formal analysis, Y.P.; investigation, Y.X.; resources, Y.X.; data curation, Z.Y.; writing—original draft preparation, Y.P.; writing—review and editing, C.W.; visualization, Z.Y.; supervision, R.L.; project administration, R.L.; funding acquisition, R.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by [National Natural Science Foundation of China] grant number [62173049].

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kuuskraa, V.A. Unconventional natural gas, industry savior or bridge? In Proceedings of the EIA Energy Outlook and Modeling Conference, Washington, DC, USA, 27 March 2006. [Google Scholar]
  2. Li, X.; Hu, S.; Cheng, K. Suggestions from the development of fractured shale gas in North America. Pet. Explor. Addit. Dev. 2007, 34, 392–400. [Google Scholar]
  3. Zhang, J.; Jin, Z.; Yuan, M. Reservoiring Mechanism of Shale Gas Additionally, Its Distribution. Nat. Gas Ind. 2004, 24, 15–18. [Google Scholar]
  4. Zhang, L.; Xu, B.; Xin, C.; Qiao, X.; Mu, J.; Xu, Y.; Han, C. Production forecasting model of shale gas considering the main fractures. Nat. Gas Geosci. 2019, 30, 247–256. [Google Scholar] [CrossRef]
  5. Li, Y.; Liu, X.; Hu, Z.; Duan, X.; Zhang, J.; Zhan, H. Summary of numerical models for predicting productivity of shale gas horizontal wells. Adv. Earth Sci. 2020, 35, 350–362. [Google Scholar]
  6. Kelvin, A.; John, W.; Derek, E. Laboratory Investigation of Impact of Slickwater Composition on Multiphase Permeability Evolution in Tight Sandstones. SPE Prod. Oper. 2021. [Google Scholar] [CrossRef]
  7. Antonia, S.; Victor, T.; Ulrich, B. Assessment of Polymer-Induced Clogging Using Microfluidics. SPE J. 2021, 26, 3793–3804. [Google Scholar]
  8. Jesse, G.; Zhu, D.; Hill, A.D. Impairment of Fracture Conductivity in the Eagle Ford Shale Formation. SPE Prod. Oper. 2018, 33, 637–653. [Google Scholar]
  9. Jansen, T.; Zhu, A.D. The Effect of Rock Mechanical Properties on Fracture Conductivity for Shale Formations. In Proceedings of the SPE Hydraulic Fracturing Technology Conference, The Woodlands, TX, USA, 3–5 February 2015. [Google Scholar]
  10. Abhinav, M.; Chandra, S.; Carl, H. A Study of Propped-Fracture Conductivity: Impairment Mechanisms Under Laboratory Conditions. In Proceedings of the SPWLA 58th Annual Logging Symposium, Oklahoma City, OK, USA, 17–21 June 2017. [Google Scholar]
  11. Zhang, J.J.; Ouyang, L.; Hill, A.D. Experimental and Numerical Studies of Reduced Fracture Conductivity due to Proppant Embedment in Shale Reservoirs. In Proceedings of the SPE Annual Technical Conference and Exhibition, Amsterdam, The Netherlands, 27–29 October 2014. [Google Scholar]
  12. Awoleke, O.O.; Zhu, D. New propped-fracture-conductivity models for tight gas sands. SPE J. 2016, 21, 1508–1517. [Google Scholar] [CrossRef]
  13. Yu, J.; Wang, J.; Wang, S. Propped Fracture conductivity Evolution during Reservoir Drawdown. In Proceedings of the 55th U.S. Rock Mechanics/Geomechanics Symposium, Houston, TX, USA, 18–25 June 2021. [Google Scholar]
  14. Wang, J.; Elsworth, D.; Denison, M.K. Proppant transport in a propagating hydraulic fracture and the evolution of transport properties. In Proceedings of the 53rd U.S. Rock Mechanics/Geomechanics Symposium, New York, NY, USA, 23–26 June 2019. [Google Scholar]
  15. Guo, X.; Zhou, C. Diffusion Seepage Model for Fractured Horizontal Well in Shale Gas Reservoir. J. Southwest Pet. Univ. 2015, 37, 38–44. [Google Scholar]
  16. Wei, Y.; Wang, J.; Yu, W.; Qi, Y.; Miao, J.; Yuan, H.; Liu, C. A smart productivity evaluation method for shale gas wells based on 3D fracture network model. Pet. Explor. Dev. 2021, 48, 787–796. [Google Scholar] [CrossRef]
  17. Pan, Y.; Wang, Y.; Che, M.; Liao, R.; Zheng, H. Post-fracturing production prediction and fracturing parameter optimization of horizontal wells based on grey relation projection random forest algorithm. J. Xi’Shiyou Univ. 2021, 36, 71–76. [Google Scholar] [CrossRef]
  18. Wang, X.; Feng, D.; Li, X.; Yan, Y.; Guo, X.; Qiao, X.; Xue, B.; Ma, C.; Wang, Q.; Lei, K.; et al. Productivity evalution of tight gas well with time-dependent mechanism: A case study of Yan’an gasfield. Acta Pet. Sin. 2019, 40, 1358–1367. [Google Scholar]
  19. Wen, Q.; Zhang, S.; Wang, L.; Liu, Y.; Li, X. The effect of proppant embedment upon the long-term conductivity of fractures. J. Pet. Sci. Eng. 2007, 55, 221–227. [Google Scholar] [CrossRef]
  20. Zhao, Y.; Zhang, S.; Xiao, F.; He, J.; Dong, Y.; Zhao, W. Experimental study on long-term flow conductivity of fractured fractures in different types of reservoirs. Sci. Technol. Eng. 2017, 17, 192–197. [Google Scholar]
  21. Yuan, X.; Xu, D.; Chen, S.; He, X.; Cheng, N.; Chen, J. The damage of fracturing fluid intrusion to the conductivity of shale reservoir. Sci. Technol. Eng. 2020, 20, 3591–3597. [Google Scholar]
  22. Wang, B.; Fidelibus, C. An Open-Source Code for Fluid Flow Simulations in Unconventional Fractured Reservoirs. Geosciences 2021, 11, 106. [Google Scholar] [CrossRef]
  23. Gregg, S.J.; Sing, K.S. Adpsorption, Surface Area and Porosity; Academic Press: London, UK, 1982. [Google Scholar]
  24. Beskok, A.; Karniadakis, G.E. Report: A model for flows in channels, pipes, and ducts at micro and nanoscales. Microscale Thermophys. Eng. 1999, 3, 43–77. [Google Scholar]
  25. Barree, R.D.; Conway, M.W. Beyond beta factors: A complete model for Darcy, Forchheimer, and trans-Forchheimer flow in porous media. In Proceedings of the SPE Annual Techical Conference and Exhibition, Houston, TX, USA, 26–29 September 2004. [Google Scholar]
  26. Lee, S.H.; Lough, M.F.; Jensen, C.L. Hierarchical modeling of flow in naturally fractured formations with multiple length scales. Water Resour. Res. 2001, 37, 443–455. [Google Scholar] [CrossRef]
  27. Peaceman, D.W. Interpretation of well-block pressures in numerical reservoir simulation. SPE J. 1978, 19, 183–194. [Google Scholar]
  28. Li, Z.; Wu, X.; Han, G.; Wang, R.; Zheng, L. Semi-analytical model of multi-stage fractured horizontal well productivity considering time-dependent fracture conductivity. Oil Drill. Prod. Technol. 2019, 41, 354–362. [Google Scholar]
  29. Chen, X. Fracturing Fluid Flowback Simulation of Fractured Horizontal Well in Shale Gas Reservoirs. Doctoral Dissertation, Southwest Petroleum University, Chengdu, China, 2019. [Google Scholar]
Figure 1. Long-term conductivity test. (a) Long-term conductivity test. (b) Long-term conductivity test (considering liquids).
Figure 1. Long-term conductivity test. (a) Long-term conductivity test. (b) Long-term conductivity test (considering liquids).
Processes 10 00801 g001
Figure 2. Flow exchange between matrix mesh and fracture mesh.
Figure 2. Flow exchange between matrix mesh and fracture mesh.
Processes 10 00801 g002
Figure 3. Fracture–fracture flow exchange.
Figure 3. Fracture–fracture flow exchange.
Processes 10 00801 g003
Figure 4. Cumulative production (model vs. Uncong).
Figure 4. Cumulative production (model vs. Uncong).
Processes 10 00801 g004
Figure 5. Formation pressure (model vs. Uncong). (a) Result from Uncong. (b) Result from model.
Figure 5. Formation pressure (model vs. Uncong). (a) Result from Uncong. (b) Result from model.
Processes 10 00801 g005
Figure 6. Formation pressure at different production times. (a) Pressure of production, 50d (constant conductivity). (b) Pressure of production, 50d (variable conductivity). (c) Pressure of production, 180d (constant conductivity). (d) Pressure of production, 180d (variable conductivity). (e) Pressure of production, 1y (constant conductivity). (f) Pressure of production, 1y (variable conductivity).
Figure 6. Formation pressure at different production times. (a) Pressure of production, 50d (constant conductivity). (b) Pressure of production, 50d (variable conductivity). (c) Pressure of production, 180d (constant conductivity). (d) Pressure of production, 180d (variable conductivity). (e) Pressure of production, 1y (constant conductivity). (f) Pressure of production, 1y (variable conductivity).
Processes 10 00801 g006
Figure 7. Pressure distribution near wellbore.
Figure 7. Pressure distribution near wellbore.
Processes 10 00801 g007
Figure 8. Near well flow distribution.
Figure 8. Near well flow distribution.
Processes 10 00801 g008
Figure 9. Formation pressure (constant conductivity vs. variable conductivity).
Figure 9. Formation pressure (constant conductivity vs. variable conductivity).
Processes 10 00801 g009
Figure 10. Cumulative production (constant conductivity vs. variable conductivity).
Figure 10. Cumulative production (constant conductivity vs. variable conductivity).
Processes 10 00801 g010
Figure 11. Formation pressure (constant conductivity vs. variable conductivity).
Figure 11. Formation pressure (constant conductivity vs. variable conductivity).
Processes 10 00801 g011
Figure 12. Cumulative production (constant conductivity vs. variable conductivity).
Figure 12. Cumulative production (constant conductivity vs. variable conductivity).
Processes 10 00801 g012
Figure 13. Comparison between actual production data and model prediction data of Well W1.
Figure 13. Comparison between actual production data and model prediction data of Well W1.
Processes 10 00801 g013
Table 1. Conductivity Damage Response Surface Factor Level.
Table 1. Conductivity Damage Response Surface Factor Level.
LevelSand
Concentration
(kg/m2)
Fluid Viscosity
(mPa·s)
Fluid Retention
Time
(h)
−15.5110
07312
18.5514
Table 2. Results of the response surface test.
Table 2. Results of the response surface test.
NumberSand Concentration (kg/m2)Fluid Viscosity (mPa·s)Fluid Retention Time (h)Conductivity Retention Rate (%)
1731243.03
2731239.96
35.511247.61
4731242.38
58.531447.43
6711056.53
7751032.89
8711451.71
98.511260.59
108.531051.08
115.551226.43
12731241.85
135.531432.46
14731241.88
158.551240.34
16751429.11
175.531036.21
Table 3. Calculated parameters.
Table 3. Calculated parameters.
ParametersValueParametersValue
Model size200 × 200mFracture permeability10,000 mD
Number of grids100 × 100Fracture porosity1
Matrix porosity0.03Fracture length100 m
Matrix permeability0.0001 mDFracture width0.003 m
Langmuir pressure9 MPaInitial formation pressure30 MPa
Langmuir volume0.004 m3/kgBottom hole pressure15 MPa
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pan, Y.; Xu, Y.; Yang, Z.; Wang, C.; Liao, R. Shale Gas Productivity Prediction Model Considering Time-Dependent Fracture Conductivity. Processes 2022, 10, 801. https://doi.org/10.3390/pr10050801

AMA Style

Pan Y, Xu Y, Yang Z, Wang C, Liao R. Shale Gas Productivity Prediction Model Considering Time-Dependent Fracture Conductivity. Processes. 2022; 10(5):801. https://doi.org/10.3390/pr10050801

Chicago/Turabian Style

Pan, Yuan, Yiwen Xu, Ze Yang, Chunli Wang, and Ruiquan Liao. 2022. "Shale Gas Productivity Prediction Model Considering Time-Dependent Fracture Conductivity" Processes 10, no. 5: 801. https://doi.org/10.3390/pr10050801

APA Style

Pan, Y., Xu, Y., Yang, Z., Wang, C., & Liao, R. (2022). Shale Gas Productivity Prediction Model Considering Time-Dependent Fracture Conductivity. Processes, 10(5), 801. https://doi.org/10.3390/pr10050801

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop