Next Article in Journal
Wind Tunnel Tests of Wake Characteristics for a Scaled Wind Turbine Model Based on Dynamic Similarity
Next Article in Special Issue
Inversion Study on Parameters of Cascade Coexisting Gas-Bearing Reservoirs in Huainan Coal Measures
Previous Article in Journal
Strategic Planning of Oil and Gas Companies: The Decarbonization Transition
Previous Article in Special Issue
Applicability Analysis of Pre-Stack Inversion in Carbonate Karst Reservoir
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Differentiation and Prediction of Shale Gas Production in Horizontal Wells: A Case Study of the Weiyuan Shale Gas Field, China

1
PetroChina Research Institute of Petroleum Exploration and Development, Beijing 100083, China
2
China National Shale Gas Research and Development (Experiment) Center, Langfang 065007, China
*
Authors to whom correspondence should be addressed.
Energies 2022, 15(17), 6161; https://doi.org/10.3390/en15176161
Submission received: 6 July 2022 / Revised: 1 August 2022 / Accepted: 1 August 2022 / Published: 25 August 2022
(This article belongs to the Special Issue Advances in Methane Production from Coal, Shale and Other Tight Rocks)

Abstract

:
The estimated ultimate recovery (EUR) of shale gas is an important index for evaluating the production capacity of horizontal wells. The Weiyuan shale gas field has wells with considerable EUR differentiation, which hinders the prediction of the production capacity of new wells. Accordingly, 121 wells with highly differentiated production are used for analysis. First, the main control factors of well production are identified via single-factor and multi-factor analyses, with the EUR set as the production capacity index. Subsequently, the key factors are selected to perform the multiple linear regression of EUR, accompanied by the developed method for well production prediction. The thickness and drilled length of Long 111 (Substratum 1 of Long 1 submember, Lower Silurian Longmaxi Formation) are demonstrated to have the uttermost effects on the well production, while several other factors also play important roles, including the fractured horizontal wellbore length, gas saturation, brittle mineral content, fracturing stage quantity, and proppant injection intensity. The multiple linear regression method can help accurately predict EUR, with errors of no more than 10%, in wells that have smooth production curves and are free of artificial interference, such as casing deformation, frac hit, and sudden change in production schemes. The results of this study are expected to provide certain guiding significances for shale gas development.

1. Introduction

The recovery technology for shallow (<3500 m) marine shale gas in China has been mature in terms of supporting techniques and fit-for-purpose equipment [1]. Starting from scratch, the shale gas industry of China reached an annual production of 1.0 × 1010 m3 within six years, which was doubled to approach a historical breakthrough in eight years [2]. Shale gas reservoirs are found in various sedimentary environments in China [3], including marine (50% of the whole country’s shale gas resources), marine–continental transitional (20%), and continental (30%) environments [4,5,6,7].
Located in the southern Sichuan Basin [8,9], the Weiyuan shale gas field has been in exploration and development since 2009, with a cumulative proven shale gas resource in place of 4.277 × 1011 m3 [10]. In 2021, the shale gas production climbs up to 4.12 × 109 m3, which accounts for 32% of the total shale gas production of China National Petroleum Corporation (CNPC). In this context, it is of critical importance to calculate the estimated ultimate recovery (EUR) of gas wells for accurately estimating the potential of shale gas recovery and realizing large-scale economic development.
EUR is an important index for evaluating the production capacity of a shale gas well, and it is affected by numerous factors [11,12], which can be analyzed by various methods [13,14]. Most popular methods include the grey correlation approach [15,16] and numerical simulations [17]. Efforts have been made to quantitatively investigate the correlation between the engineering parameters and the fracture network effectiveness via comparison of reservoir attributes and analysis of well performances [18,19]. Moreover, numerous methods have been developed for calculating the EUR of shale gas wells, including the empirical method, the modern production decline analysis, the analytical method, the linear flow analysis, and the big data approach [20]. Specifically, the pressure transient and Blasingame production decline analysis were used to calculate the EUR of a well with short-term production [21]. Moreover, innovative models for the multi-stage fractured horizontal well [22] and those based on the deep neural network [23] are useful ways of predicting the production of shale gas wells. Niu et al. [24] forecast the production capacity of shale gas wells of the Weiyuan gas field via multiple regression. Furthermore, Yin et al. [25] developed a production capacity assessment model for shale gas wells by building and solving the dual-porosity mathematical model of shale gas reservoirs. Accurately calculating EUR of wells and analyzing the main control factors of the production capacity of shale gas wells are of great significance for improving well production and guiding the deployment of shale gas recovery.
In the Weiyuan shale gas field, the production of horizontal wells is highly differentiated from well to well, and thus systematic efforts should be made to clarify the main control factors of gas well production capacity. In this research, single-factor analysis and grey correlation analysis are used to identify key factors affecting the production capacity of wells, while multiple linear regression is used to predict the production capacity of a single well. The findings of this research are expected to guide shale gas well operation in the Weiyuan gas field and enable high-efficiency recovery of shale gas.

2. Geological Setting

The Weiyuan gas field is located in the Weiyuan and Zizhong Counties in the southwestern Sichuan Basin. It lies in the Southwest Sichuan gentle fold zone of the Central Sichuan uplift and presents itself as a large dome anticline (Figure 1). The main layers for shale gas recovery are within the Wufeng Formation–Longmaxi Formation. The Wufeng–Longmaxi shale is buried at 2000–4000 m, with a thickness of 180–600 m. The interval from the Wufeng Formation to the Long 1 Sub-Member of the Longmaxi Formation is the target interval, with a shale thickness of 43.9–54.8 m, averaging 46.4 m. The Long 1 Sub-Member is further divided into the Long 111–114 small layers from the bottom to top, among which the Long 111 is the main layer for recovery, with a thickness of 1.7–7.0 m (4.6 m on average). The shale in the interval from the Wufeng Formation to the Long 1 Sub-Member has a TOC of 2.7–3.6%, averaging 3.2%; a porosity of 5.2–6.7%, averaging 5.9%; a gas saturation of 32.7–84.6%, averaging 64.7%; a gas content of 3.3–8.5 m3/t, averaging 5.5 m3/t; a brittle mineral content of 60–82%, averaging 74%; and a pressure coefficient of 1.2–2.0 in the favorable zone for production capacity building.

3. Data and Analysis Methods

By the end of 2021, a total of 396 wells have been brought into production in the Weiyuan block and these wells present an average testing daily production of 2.16 × 105 m3/d and an average single-well EUR of 8.9 × 107 m3. A total of 121 wells that were brought into production within the recent three years and featured considerably differentiated production are selected as the samples for analysis, which present the average testing daily production of 2.488 × 105 m3/d, the average first-year daily production of 9.4 × 104 m3/d, and the EUR of (0.31–2.22) × 108 m3 (averaging 1.06 × 108 m3; Figure 2; Table 1). Among these wells, low-production wells with the EUR lower than 0.6 × 108 m3 account for 8%; medium-production wells with the EUR of (0.6–1.2) × 108 m3, 59%; and high-production wells with the EUR above 1.2 × 108 m3, 33%. A total of 62 wells are found with EUR below 1.0 × 108 m3, which account for about 50% of the wells brought into production. The overall distribution of EUR is even and the production of wells is greatly differentiated.

3.1. Factors Affecting Gas Well Production

Various factors affect the production capacity of shale gas wells, which include geological and engineering factors. The geological factors include the thickness of the Long 111 interval, TOC, brittle mineral content, gas content, gas saturation, and pressure coefficient, while the engineering factors include the drilled length of the Long 111 interval, the fractured horizontal wellbore length, fracturing stage quantity, liquid injection intensity, and proppant injection intensity.

3.2. Grey Correlation Analysis

The production capacity of shale gas wells is subjected to the joint effects of multiple factors, which may interact with each other. The basic regression method fails to deliver a satisfactory analysis of the main control factors, as it cannot quantitatively capture the root causes of the differentiated production capacity of wells [26]. Hence, the grey correlation analysis is performed to quantify the intensity of impacts of each factor on the well production capacity and correspondingly identify the main control factors of the production capacity of shale gas horizontal wells in the Weiyuan gas field.
The grey system theory was developed by the Chinese scholar Julong Deng in 1982. This theory deals with the “limited-sample” system with poor data (lean in information) with partially known and partially unknown information [27]. The grey correlation analysis is an important analytical method for the grey system theory and is highly valuable in investigating the correlations among variables [28]. When no strict mathematical relationships exist between each influential factor and the total result, the grey correlation method can deliver an effective analysis of data and describe the strength, magnitude, and ranking of relationships between factors [29]. The grey correlation analysis determines the closeness of the relationship, according to the geometric similarity between the reference sequence and analysis sequence, and can be used to determine how intensive the effects of each factor are on the result [30].
The grey correlation analysis is performed after normalizing the above data of the factors. The temporal sequence of the EUR is set as the reference sequence X0, while the temporal sequences of each factor, namely the Long 111 thickness, TOC, brittle mineral content, gas content, gas saturation, pressure coefficient, drilled Long 111 length, fractured horizontal wellbore length, fracturing stage quantity, liquid injection intensity, and proppant injection intensity, are used as the analysis sequences. The following matrix is built:
X 0 , X 1 , , Xn = x 0 1 x 1 1 x n 1 x 0 2 x 1 2 x n 2 x 0 m x 1 m x n m
Then, the absolute difference between the elements of each analysis sequence and the reference sequence is calculated one after another:
Δ 0 i k = x 0 k x i k k = 1 , , m , i = 1 , , n
Accordingly, the absolute difference matrix is formed:
Δ 01 1 Δ 02 1 Δ 0 n 1 Δ 01 2 Δ 02 2 Δ 0 n 2 Δ 01 m Δ 02 m Δ 0 n m
The maximum value in the absolute difference matrix is defined as the maximum difference, while the minimum one is identified as the minimum difference:
Δ m a x = max Δ 0 i k
Δ m i n = min Δ 0 i k
The correlation coefficient is computed as below:
ξ 0 i k = Δ m i n + ρ Δ m a x Δ 0 i k + ρ Δ m a x
where ρ is the resolution ratio determining the effects of Δ m a x on the data conversion; ξ 0 i k is the correlation coefficient with a positive value of no more than 1. The function of ρ is to control the influence of the ∆(max) on the overall correlation degree, so that the resolution of the correlation degree between different factors at the longitudinal level is as large as possible. A lower ρ results in promoted differentiation of correlation coefficients. In this research, ρ is empirically set as 0.5.
The grey correlation degrees between factors are computed:
γ 0 i = 1 N k = 1 N ξ 0 i k   i = 1 , 2 , 3 , , n
where γ 0 i is the correlation degree, a measure of the correlation between the reference sequence X 0 and the analysis sequence X i i = 1 , 2 , 3 , , n .
The weight coefficients of each factor are determined using the following equation, according to the correlation degrees:
W i = γ 0 i i = 1 n γ 0 i

3.3. Principles of the Multiple Regression Algorithm

The multiple linear regression (MLR) analysis is a mathematical analysis method based on the correlation between the dependent and independent variables [31,32]. With a dependent variable related to multiple factors, MLR can be performed for analysis [33], in which an MLR model is developed via the weighted summation of these factors [34]. It is assumed that the dependent variable y can be expressed as a linear combination of m independent variables:
y = β 0 + β 1 x 1 + β 2 x 2 + + β i x m + ε
where y is the dependent variable; x 1 ,   x 2 x i are the independent variables; β 0 ,   β 1 ,   β 2 β m are the regression coefficients; and ε is the random error with the normal distribution N 0 , δ 2 .
Providing that there are N sets of observation data (samples), the theoretical equation shall be:
y = X β + ϵ
where:
y = y 1 y 2 y n
X = 1 x 11 1 x 21 x 1 m x 2 m 1 x n 1 x n m
ϵ = ϵ 1 ϵ 2 ϵ n
β = β 0 β 1 β m
In accordance with the least-square principle, the difference between the regression estimate and the true value (also known as the residual error) is calculated and the regression coefficients are computed by searching the least sum of squares of residual errors. Hence, the regression equation is obtained.
The test of significance can be performed for the regression coefficients of the MLR equation to validate whether or not a linear relationship exists between variables. On one hand, the regression coefficients that are not all zero stand for significant linear relationships between variables. On the other hand, regression coefficients that are all zero represent no significant linear relationship.

4. Results and Discussion

4.1. Single-Factor Analysis

As is shown in Figure 3, EUR is found with relatively large positive correlations with the Long 111 thickness, gas saturation, drilled length of the Long 111, fractured horizontal wellbore length, and fracturing stage quantity, whereas it presents smaller correlations with TOC, brittle mineral content, gas content, pressure coefficient, liquid injection intensity, and proppant injection intensity.
The Pearson correlation analysis is implemented on EUR and 11 influential factors after the Z-score normalization to reveal the linear relationships among these factors and those between each factor and the EUR (Table 2 and Figure 4). Significant correlations are seen between EUR and Long 111 thickness, brittle mineral content, gas saturation, drilled Long 111 length, fractured horizontal wellbore length, fracturing stage quantity, and proppant injection intensity, whereas weak relationships are observed between EUR and TOC, gas content, pressure coefficient, and liquid injection intensity.

4.2. Grey Correlation Analysis

The correlation degrees between each factor and EUR are shown in Table 3, while the corresponding weight coefficients are presented in Figure 5. Clearly, the thickness and drilled length of Long 111 have the highest effects on the EUR, followed by the second tier composed of the fractured horizontal wellbore length, gas saturation, brittle mineral content, and fracturing stage quantity, and the third tier consisting of the proppant injection intensity, liquid injection intensity, TOC, pressure coefficient, and gas content. According to the results, for the first seven factors, the results are consistent with single-factor analysis and grey correlation analysis, which indicates that the results are reliable.

4.3. MLR Analysis

According to the results of the grey correlation analysis of the factors affecting well production, seven key factors that are highly correlated with the EUR are identified and selected for the MLR analysis. These factors are the Long 111 thickness, brittle mineral content, gas saturation, drilled Long 111 length, fractured horizontal wellbore length, fracturing stage quantity, and proppant injection intensity.
As shown in Table 4, the calculated significance is 1.1 × 10−15, far lower than the 0.05 significance level, and the constant for the MLR is non-zero, which means that the regression performance of the MLR equation is satisfactory and the results of the linear regression analysis between the EUR and these factors are reliable. The MLR fitting formula of the EUR of a shale gas horizontal well in relation to the above factors can be written as(Table 5): the EUR of a shale gas well = −1.008 + 0.143 × Long 111 thickness + 0.006 × brittle mineral content + 0.001 × gas saturation + 1.6 × 10−4 × drilled Long 111 length − 4.9 × 10−6 × fractured horizontal wellbore length + 0.018 × fracturing stage quantity + 0.022 × proppant injection intensity.

4.4. Application

The EUR of 20 shale gas wells of the Weiyuan gas field is calculated using the MLR production capacity model presented above and compared with that computed using the analytical method. The analytical method considers the shale reservoir’s physical properties, fracture parameters, fracture conductivity, and other model parameters. Moreover, it establishes the shale gas horizontal well model, with the model parameters historically fitted to correct the model parameters. After the historical fitting is completed, the production is predicted. Generally, the analytical method is used to calculate EUR for shale gas and the results are reliable in most cases. The results show (Table 6) that production predictions of 4 wells are overestimated, those of 5 are underestimated, and those of 11 are found with relatively accurate estimates, with errors no more than 10%. The matching rate of the model is around 50%.
Furthermore, the production performance of the wells presenting accurate estimates (Figure 6) are typically found with continuous production and pressure curves, with no artificial interference at the time being brought into production and during the production (e.g., engineering-induced casing deformation, frac hit, and sudden change in the production scheme). In contrast, the typical production performances of the overestimated wells (Figure 7) reveal that the wells are subjected to frac hit from the adjacent well and present discontinuous production curves; accordingly, the actual EUR far deviates from the ideal EUR predicted by the MLR model. Moreover, the typical production performances of the underestimated wells (Figure 8) show a change in the production scheme around the 350th day of production, and the production rate and pressure remain stable for a rather long period; as a result, the predicted EUR is lower than the actual EUR. The above analysis suggests that effective screening of sample data is required to select wells with continuous production decline for the MLR analysis in an attempt to avoid artificial interference and improve the accuracy of the MLR model.

5. Conclusions

(1)
First, the main geological and engineering factors that control the production capacity of horizontal wells in the Weiyuan shale gas field are clarified via the single-factor analysis and grey correlation analysis, with 121 horizontal wells used as the samples. The primary control factors are the thickness and drilled length of Long 111, while the secondary control factors are the fractured horizontal wellbore length, gas saturation, brittle mineral content, fracturing stage quantity, and proppant injection intensity.
(2)
For the production forecast of wells based on the multiple linear regression, one needs to select wells with continuous production curves and that are free of artificial interference when being brought into production or during production, such as casing deformation, frac hit, and sudden change in production schemes. For such wells, the model presented in this research can accurately predict the EUR of wells. For the wells newly brought into production, the model presented in this research can rapidly predict the ideal production capacity and the EUR of wells, with errors of no more than 10% compared to the analytical results.

Author Contributions

Conceptualization, L.K., W.G. and X.Z.; Data curation, L.K., W.G., Y.L. and Z.S.; Formal analysis, X.Z., Y.L. and Z.S.; Investigation, L.K., W.G., Y.L. and Z.S.; Methodology, L.K., X.Z. and Y.L.; Project administration, W.G. and X.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by [Study on productivity evaluation and development technology policy optimization of Marine and continental transitional shale gas] grant number [2021DJ2005] and [Study on shale gas optimization and scale efficiency development technology in southern Sichuan] grant number [2022KT1204].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xu, C.H.; Zhong, Y.H.; Wang, C.S.; Ke, L.; Chou, W.D.; Zhao, Q. A brief analysis of key points and early production characteristics of shale gas exploration and development in China. J. Geol. 2021, 45, 197–212. [Google Scholar]
  2. Zou, C.N.; Zhao, Q.; Cong, L.Z.; Wang, H.Y.; Shi, Z.S.; Wu, J.; Pan, S.Q. Development progress, potential and prospect of shale gas in China. Nat. Gas Ind. 2021, 41, 1–14. [Google Scholar]
  3. Feng, Y.; Ning, Z.F.; Liu, H.Q. Fractal characteristics of shales from a shale gas reservoir in the Sichuan Basin, China. Fuel 2014, 115, 378–384. [Google Scholar]
  4. Zou, C.N.; Yang, Z.; Sun, S.S.; Zhao, Q.; Bai, W.H.; Liu, H.L.; Pan, S.Q.; Wu, S.T.; Yuan, Y.L. Exploring petroleum inside source kitchen: Shale oil and gas in Sichuan Basin. Sci. China Earth Sci. 2020, 50, 903–920. [Google Scholar] [CrossRef]
  5. Zhao, W.Z.; Jia, A.L.; Wei, Y.S.; Wang, J.L.; Zhu, H.Q. Progress in shale gas exploration in China and prospects for future development. China. Pet. Explor. 2020, 25, 31–44. [Google Scholar]
  6. Li, Y.; Chen, J.-Q.; Yang, J.-H.; Liu, J.-S.; Tong, W.-S. Determination of shale macroscale modulus based on microscale measurement: A case study concerning multiscale mechanical characteristics. Pet. Sci. 2022, 19, 1262–1275. [Google Scholar] [CrossRef]
  7. Li, Y.; Wang, Z.; Pan, Z.; Niu, X.; Meng, S. Pore structure and its fractal dimensions of transitional shale: A cross-section from east margin of the ordos basin, china. Fuel 2019, 241, 417–431. [Google Scholar] [CrossRef]
  8. Liang, F.; Zhang, Q.; Lu, B.; Jiang, W.; Xiong, X.L.; Chen, P.; Jiang, R.; Ma, C. Lithofacies and distribution of Wufeng-Longmaxi organic-rich shale and its impact on shale gas production in Weiyuan shale gas play, southern Sichuan Basin, China. Acta Sedimentol. Sin. 2022; in press. [Google Scholar]
  9. Ma, X.H.; Li, X.Z.; Liang, F.; Wan, Y.J.; Shi, Q.; Wang, Y.H.; Zhang, X.W.; Che, M.G.; Guo, W.; Guo, W. Dominating factors on well productivity and development strategies optimization in Weiyuan shale gas play, Sichuan Basin, SW China. Pet. Explor. Dev. 2020, 47, 555–563. [Google Scholar] [CrossRef]
  10. He, X.; Wu, J.F.; Yong, R.; Zhao, S.X.; Zhou, X.J.; Zhang, J.D.; Zhang, D.L.; Zhong, C.X. Accumulation conditions and key exploration and development technologies of marine shale gas field in Changning-Weiyuan block, Sichuan Basin. Acta Pet. Sin. 2021, 42, 259–272. [Google Scholar]
  11. Ma, X. "Extreme utilization" development theory of unconventional natural gas. Pet. Explor. Dev. 2021, 48, 326–336. [Google Scholar] [CrossRef]
  12. Guo, X.S. Controlling factors on shale gas accumulations of Wufeng-Longmaxi Formations in Pingqiao shale gas field in Fuling area, Sichuan Basin. Nat. Gas Geosci. 2019, 30, 1–10. [Google Scholar] [CrossRef]
  13. Wang, P.W.; Zou, C.; Li, X.J.; Jiang, L.W.; Li, J.J.; Mei, J.; Zhang, Z.; Li, Q.F. Main geological controlling factors of shale gas enrichment and high yield in Zhaotong demonstration area. Acta Pet. Sin. 2018, 39, 744–753. [Google Scholar]
  14. Nie, H.K.; Chen, Q.; Zhang, G.R.; Sun, C.X.; Wang, P.W.; Lu, Z.Y. An overview of the characteristic of typical Wufeng–Longmaxi shale gas fields in the Sichuan Basin, China. Nat. Gas Ind. B 2021, 8, 217–230. [Google Scholar] [CrossRef]
  15. He, C.; Wan, Y.J.; Su, Y.H.; Luo, W.J.; Geng, X.Y. Comprehensive Analysis of Influencing Factors for Shale Gas Well Production in Weiyuan Gas Field. In Proceedings of the International Field Exploration and Development Conference 2019, Xi’an, China, 16 October 2019. [Google Scholar]
  16. Chai, Y. Study on the Main Controlling Factors of Shale Gas Production Based on Grey Correlation Method. J. Chon. Uni. Sci. Tech. Nat. Sci. Ed. 2018, 20, 32–34. [Google Scholar]
  17. Jayakumar, R.; Rai, R. Impact of Uncertainty in Estimation of Shale-Gas-Reservoir and Completion Properties on EUR Forecast and Optimal Development Planning: A Marcellus Case Study. SPE Res. Eval. Eng. 2014, 17, 60–73. [Google Scholar] [CrossRef]
  18. Chen, X.; Xu, J.L.; Li, J.; Xiao, J.F.; Zhong, S.C. An Analysis of Main Controlling Factors of Production for Horizontal Shale Gas Well in Weiyuan Block. J. Southwest. Pet. Uni. Sci. Tech. Ed. 2020, 41, 63–74. [Google Scholar]
  19. Jia, C.Y.; Jia, A.L.; He, D.B.; Wei, Y.S.; Qi, D.Y.; Wang, J.L. Key factors influencing shale gas horizontal well production. Nat. Gas Ind. 2017, 37, 80–88. [Google Scholar]
  20. Zhao, X.; Rui, Z.; Liao, X.; Zhang, R. A simulation method for modified isochronal well testing to determine shale gas well productivity. J. Nat. Gas Sci. Eng. 2015, 27, 479–485. [Google Scholar] [CrossRef]
  21. Xu, Y.; Liu, Q.; Li, X.; Meng, Z.; Yang, S.; Tan, X. Pressure transient and Blasingame production decline analysis of hydraulic fractured well with induced fractures in composite shale gas reservoirs. J. Nat. Gas Sci. Eng. 2021, 94, 104058. [Google Scholar] [CrossRef]
  22. Cui, Y.; Jiang, R.; Wang, Q.; Liu, X. Production performance analysis of multi-fractured horizontal well in shale gas reservoir considering space variable and stress-sensitive fractures. J. Pet. Sci. Eng. 2021, 207, 109071. [Google Scholar] [CrossRef]
  23. Han, D.; Kwon, S. Application of Machine Learning Method of Data-Driven Deep Learning Model to Predict Well Production Rate in the Shale Gas Reservoirs. Energies 2021, 14, 3629. [Google Scholar] [CrossRef]
  24. Niu, W.; Lu, J.; Sun, Y. A Production Prediction Method for Shale Gas Wells Based on Multiple Regression. Energies 2021, 14, 1461. [Google Scholar] [CrossRef]
  25. Yin, D.; Wang, D.; Zhang, C.; Duan, Y. Shale Gas Productivity Predicting Model and Analysis of Influence Factor. Open Pet. Eng. J. 2015, 8, 203–207. [Google Scholar]
  26. Guo, W. Differentiation Analysis on Shale Gas Production of the Horizontal Wells in Sichuan Weiyuan Block. Sci. Tech. Eng. 2018, 18, 228–233. [Google Scholar]
  27. Deng, J.L. Introduction to the grey system theory. Inner Mongolia Elec. Power. 1993, 3, 51–52. [Google Scholar]
  28. Fang, B.; Hu, J.; Xu, J.; Zhang, Y. A semi-analytical model for horizontal-well productivity in shale gas reservoirs: Coupling of multi-scale seepage and matrix shrinkage. J. Pet. Sci. Eng. 2020, 195, 107869. [Google Scholar] [CrossRef]
  29. Tan, X.R.; Deng, J.L. A new method for the multi-variable statistical analysis. Stat. Res. 1995, 3, 46–48. [Google Scholar]
  30. Liu, Y.; Zhang, X.; Shi, J.; Guo, W.; Kang, L.; Yu, R.; Sun, Y.; Wang, Z.; Pan, M. A reservoir quality evaluation approach for tight sandstone reservoirs based on the gray correlation algorithm: A case study of the Chang 6 layer in the W area of the as oilfield, Ordos Basin. Energy Explor. Exploit. 2021, 39, 1027–1056. [Google Scholar] [CrossRef]
  31. Zhou, Y.X.; Zhao, A.K.; Yu, Q.; Zhang, D.; Zhang, Q.; Lei, Z.H. A new method for evaluating favorable shale gas exploration areas based on multi-linear regression analysis: A case study of marine shales of Wufeng-Longmaxi Formations, Upper Yangtze Region. Sediment. Geol. Teth. Geol. 2021, 41, 387–397. [Google Scholar]
  32. Liu, Y.; Lai, F.; Zhang, H.; Tan, Z.; Wang, Y.; Zhao, X.; Tan, X. A novel mineral composition inversion method of deep shale gas reservoir in Western Chongqing. J. Pet. Sci. Eng. 2021, 202, 108528. [Google Scholar]
  33. Tang, L.; Song, Y.; Li, Q.; Pang, X.; Jiang, Z.; Li, Z.; Tang, X.; Yu, H.; Sun, Y.; Fan, S.; et al. A Quantitative Evaluation of Shale Gas Content in Different Occurrence States of the Longmaxi Formation: A New Insight from Well JY-A in the Fuling Shale Gas Field, Sichuan Basin. Acta Geol. Sin. Engl. Ed. 2019, 93, 400–419. [Google Scholar] [CrossRef]
  34. Li, C.Y. Analysis of Affecting Factors and Productivity Prediction of Shale Gas Well in ZT Block. Master’s Thesis, Xi’an Shiyou University, Xi’an, China, 2020. [Google Scholar]
Figure 1. Structural map of the Weiyuan shale gas field.
Figure 1. Structural map of the Weiyuan shale gas field.
Energies 15 06161 g001
Figure 2. Probability distribution of the EUR of the 121 sample wells.
Figure 2. Probability distribution of the EUR of the 121 sample wells.
Energies 15 06161 g002
Figure 3. Single-factor regression of factors affecting EUR.
Figure 3. Single-factor regression of factors affecting EUR.
Energies 15 06161 g003
Figure 4. Heat map of the single-factor analysis results.
Figure 4. Heat map of the single-factor analysis results.
Energies 15 06161 g004
Figure 5. Weight coefficients characterizing the effects of each factor on the EUR of wells.
Figure 5. Weight coefficients characterizing the effects of each factor on the EUR of wells.
Energies 15 06161 g005
Figure 6. The production history of Well W5 (with the accurate prediction).
Figure 6. The production history of Well W5 (with the accurate prediction).
Energies 15 06161 g006
Figure 7. The production history of Well W1 (overestimated).
Figure 7. The production history of Well W1 (overestimated).
Energies 15 06161 g007
Figure 8. The production history of Well W17 (underestimated).
Figure 8. The production history of Well W17 (underestimated).
Energies 15 06161 g008
Table 1. Parameter statistics of the sample wells of the Weiyuan shale gas field.
Table 1. Parameter statistics of the sample wells of the Weiyuan shale gas field.
ParameterUnitSampling QuantityMinMaxAverageStandard Deviation
EUR108 m31210.312.221.060.38
Long 111 interval thicknessm1212.407.605.371.37
TOC%1213.306.755.220.53
Brittle mineral content%1216396835.43
Gas contentm3/t1215.4310.307.341.03
Gas saturation%1217184763.60
Pressure coefficientdimensionless1211.402.051.740.18
Drilled length of the Long 111 intervalm12148225151593344.75
Fractured horizontal wellbore lengthm12189925771664328.16
Fracturing stage quantityStages1211636234.36
Liquid injection intensitym3/m1211648274.13
Proppant injection intensityt/m1211.083.502.080.52
Table 2. Results of the Pearson correlation analysis.
Table 2. Results of the Pearson correlation analysis.
EURLong 111 Interval ThicknessTOCBrittle Mineral ContentGas ContentGas SaturationPressure CoefficientDrilled Length of the Long 111 IntervalFractured Horizontal Wellbore LengthFracturing Stage QuantityLiquid Injection IntensityProppant Injection Intensity
EURCorrelation1.000.663 **(0.08)0.325 **(0.14)0.482 **0.140.526 **0.499 **0.298 **0.080.292 **
Significance 0.000.370.000.120.000.140.000.000.000.370.00
Long 111 Interval ThicknessCorrelation0.663 **1.000.100.338 **0.030.709 **0.262 **0.524 **0.494 **0.08−0.188 *0.481 **
Significance0.00 0.300.000.730.000.000.000.000.390.040.00
TOCCorrelation(0.08)0.101.00(0.15)0.584 **0.286 **0.435 **(0.17)(0.16)(0.10)0.100.246 **
Significance0.370.30 0.100.000.000.000.060.090.270.290.01
Brittle Mineral ContentCorrelation0.325 **0.338 **(0.15)1.00−0.190 *0.250 **0.224 *0.282 **0.120.09(0.06)(0.06)
Significance0.000.000.10 0.040.010.010.000.180.310.540.51
Gas ContentCorrelation(0.14)0.030.584 **−0.190*1.000.030.595 **(0.12)(0.06)(0.12)(0.04)0.232*
Significance0.120.730.000.04 0.770.000.190.530.190.640.01
Gas SaturationCorrelation0.482 **0.709 **0.286 **0.250 **0.031.000.180.316 **0.363 **0.11(0.09)0.421 **
Significance0.000.000.000.010.77 0.050.000.000.230.330.00
Pressure CoefficientCorrelation0.140.262 **0.435 **0.224 *0.595 **0.181.000.01(0.11)(0.04)(0.04)(0.13)
Significance0.140.000.000.010.000.05 0.900.210.690.680.17
Drilled Length of the Long 111 IntervalCorrelation0.526 **0.524 **(0.17)0.282 **(0.12)0.316 **0.011.000.753 **0.348 **(0.17)0.336 **
Significance0.000.000.060.000.190.000.90 0.000.000.060.00
Fractured Horizontal Wellbore LengthCorrelation0.499 **0.494 **(0.16)0.12(0.06)0.363 **(0.11)0.753 **1.000.534 **−0.24 **0.403 **
Significance0.000.000.090.180.530.000.210.00 0.000.010.00
Fracturing Stage QuantityCorrelation0.298 **0.08(0.10)0.09(0.12)0.11(0.04)0.348 **0.534 **1.000.35 **(0.17)
Significance0.000.390.270.310.190.230.690.000.00 0.000.06
Liquid Injection IntensityCorrelation0.08−0.188 *0.10(0.06)(0.04)(0.09)(0.04)(0.17)−0.24 **0.35 **1.00(0.11)
Significance0.370.040.290.540.640.330.680.060.010.00 0.22
Proppant Injection IntensityCorrelation0.292 **0.481 **0.246 **(0.06)0.232 *0.421 **(0.13)0.336 **0.403 **(0.17)(0.11)1.00
Significance0.000.000.010.510.010.000.170.000.000.060.22
Notes: * represents a significant correlation; ** represents an extremely significant correlatio.
Table 3. Calculated correlation degrees between each factor and the EUR.
Table 3. Calculated correlation degrees between each factor and the EUR.
FactorLong 111 Interval ThicknessTOCBrittle Mineral ContentGas ContentGas SaturationPressure CoefficientDrilled Length of the Long 111 IntervalFractured Horizontal Wellbore LengthFracturing Stage QuantityLiquid Injection IntensityProppant Injection Intensity
Correlation Degree0.8100.7370.7710.7290.7740.7300.7900.7860.7660.7440.744
Ranking1951141023687
Table 4. Variance analysis of the MLR model.
Table 4. Variance analysis of the MLR model.
Variance AnalysisSum of SquaresDegree of FreedomMean SquareFSignificance
Regression8.93371.27617.6971.10 × 10−15
Residual8.1481130.072
Total17.081120
Table 5. Results of the MLR calculation.
Table 5. Results of the MLR calculation.
ConstantLong 111 Interval ThicknessBrittle Mineral ContentGas SaturationDrilled Long 111 LengthFractured Horizontal Wellbore LengthFracturing Stage QuantityProppant Injection Intensity
−1.0080.1430.0060.0011.6 × 10−4−4.9 × 10−60.0180.022
Table 6. Comparison of the results between MLR and the analytical method.
Table 6. Comparison of the results between MLR and the analytical method.
Prediction PerformanceWell No.Testing Daily Production (104 m3/d)First-Year Daily Production (104 m3/d)EUR Based on the Analytical Method (108 m3)EUR Based on MLR (108 m3)Relative Error
OverestimationW114.375.140.640.8430.3%
W224.9311.540.981.2527.5%
W328.5710.431.081.2010.8%
W428.5610.420.901.0010.5%
Accurate EstimationW533.1213.151.411.549.4%
W613.986.420.890.968.4%
W717.586.840.870.893.0%
W843.5514.541.581.601.4%
W916.417.790.900.900.7%
W1031.4313.881.361.360.0%
W1117.458.400.970.97−0.3%
W12218.280.900.88−2.2%
W1332.6511.831.351.31−3.1%
W1421.0413.501.301.20−8.0%
W1522.3810.721.321.19−9.7%
UnderestimationW1628.1910.321.080.88−18.6%
W1731.0113.021.311.05−20.1%
W1835.5519.421.711.28−25.2%
W1931.0310.721.521.10−27.9%
W2040.4516.801.971.35−31.3%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kang, L.; Guo, W.; Zhang, X.; Liu, Y.; Shao, Z. Differentiation and Prediction of Shale Gas Production in Horizontal Wells: A Case Study of the Weiyuan Shale Gas Field, China. Energies 2022, 15, 6161. https://doi.org/10.3390/en15176161

AMA Style

Kang L, Guo W, Zhang X, Liu Y, Shao Z. Differentiation and Prediction of Shale Gas Production in Horizontal Wells: A Case Study of the Weiyuan Shale Gas Field, China. Energies. 2022; 15(17):6161. https://doi.org/10.3390/en15176161

Chicago/Turabian Style

Kang, Lixia, Wei Guo, Xiaowei Zhang, Yuyang Liu, and Zhaoyuan Shao. 2022. "Differentiation and Prediction of Shale Gas Production in Horizontal Wells: A Case Study of the Weiyuan Shale Gas Field, China" Energies 15, no. 17: 6161. https://doi.org/10.3390/en15176161

APA Style

Kang, L., Guo, W., Zhang, X., Liu, Y., & Shao, Z. (2022). Differentiation and Prediction of Shale Gas Production in Horizontal Wells: A Case Study of the Weiyuan Shale Gas Field, China. Energies, 15(17), 6161. https://doi.org/10.3390/en15176161

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