Geofluids
Volume 2021 (2021), Article ID 6637322, 7 pages
https://doi.org/10.1155/2021/6637322

Research Article

A Novel Methodology for Production Decline Analysis in Tight Gas Reservoirs: A Case Study from Sulige Tight Gas Field

Xiaofeng Li,1,2 Lili Liu,1,2 Jiaojiao Liu,1,2 Zhichao Li,1,2 and Xiaoliang Zhao³

¹PetroChina Changqing Oilfield Company Research Institute of Exploration and Development, Xi’an, China
²National Engineering Laboratory for Exploration and Development of Low Permeability Oil/Gas Field, Xi’an, China
³College of Petroleum Engineering, China University of Petroleum–Beijing, China

Correspondence should be addressed to Jiaojiao Liu; liujj69_cq@petrochina.com.cn

Received 6 December 2020; Revised 11 March 2021; Accepted 31 March 2021; Published 26 April 2021

Academic Editor: Bao Jia

Copyright © 2021 Xiaofeng Li 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.

Abstract.

Study shows that decline exponents of gas wells in Sulige Tight Gas Field vary during their life cycle. However, decline analysis methods with variable decline exponent have not been developed so far. To address this problem, a simple-yet-effective method is proposed. In this method, the fracture linear flow regime and the channel linear flow regime, which appear two straight lines on the coordinate axis ( vs. ), respectively, are used to carry out the decline analysis. Each regime has a different slope and intercept (in the form of vs. ), which leads to a useful graphical technique for predicting gas rate. The graphical technique is verified by matching actual gas rate and predicting future gas rate trend. Theoretically, the influence of the two slopes and intercepts in the graphical technique on the rate decline is also studied. Aiming at decline analysis for tight gas wells with variable decline exponents, this paper proposes a novel method using linear flow characteristics, which ingeniously avoids the establishment of an empirical method to deal with variable decline exponent. The method in this study can help for better understanding of decline analysis of tight gas wells in a theoretical manner.

1. Introduction

As we all know, Arps decline-curve method [1] is not applicable for tight gas reservoir. The main reasons are that its assumptions are violated in tight gas wells. For example, tight gas wells rarely reach boundary-dominated flow even after several years of production, while Arps decline-curve method demands boundary-dominated flow. Moreover, Kupchenko [2] proved that decline exponent is variable for fractured tight gas wells. In linear flow regimes, decline exponent will be bigger than 1, until flow regime enters into boundary-dominated flow. For gas wells in some extreme low-permeability gas reservoirs, the decline exponent will be bigger than 1 in their life cycle.

To deal with the new problem, new methods are proposed. Duong [3] introduced an empirically decline model based on long-term linear flow for tight gas reservoirs. ILK et al. [4] proposed a method named Power law Exponential Decline to predict the reserves of unconventional reservoirs. For A. N. Duong method and ILK method, there are many parameters to be determined in the fitting process, so multiple solutions are inevitable. Matter et al. [5] proposed modified power law exponential decline. In 2009, Valko [6] proposed stretched exponential production decline method for shale gas. Joshi et al. [7] proposed a method which combined Duong model with hyperbolic decline model. Yu [8] presented a new improved methodology to determine a more accurate to be used in the Arps decline curve analysis for tight gas reservoirs by developing relationship between Qcum, Qcum, -, qt, and . In general, the above methodologies are mainly empirical. Neal and Mian [9] presented a predictive technique by introducing linear flow equation. Dought and Moridis [10] developed a simple, Excel-based tool for the analysis of the complex problem of gas production from a hydraulically fractured tight/shale gas reservoir, based on curve fitting a semianalytical solution to production decline data. In Alem et al.’s [11] opinions, the traditional hyperbolic decline equation can be used to predict recovery from tight gas plays by selecting the right decline exponent in correspondence with flow regime. Mienzan and Asumadu [12] present a new set of rate-decline type curves to analyze and predict gas well performance. The rate-decline type curves have been developed based on a semianalytical model. However, the method they proposed whether is applicable to tight gas is not mentioned.

Unlike the approaches mentioned above, this paper developed a new approach to handle rate decline with variable decline exponent by utilizing fracture linear flow equation and channel flow equation, rather than establishing empirical equations.

In this paper, firstly, the field characteristic is presented. Then, we employ a numerical model to investigate gas rate and decline exponents in Sulige gas field. Third, the theory basis to draw the novel approach in this paper is elaborated. Forth, practical production data is used to validate the approach. Fifth, sensitivity analysis of slope() and intercept() is performed. Finally, a discussion about the novel approach is presented.

2. Field Characteristic

Sulige gas field is a large lithologic gas reservoir, characterized by braided river development. The effective reservoir formations are mainly isolated and stripped, with strong heterogeneity. Average formation thickness is 6~7 m, porosity is 7.4~8.3%, and average permeability is 0.036 md. Reservoir depth ranges from 3450 to 3730 m; average formation pressure is 30 MPa, and formation temperature is 110°C.

3. Decline Exponent Investigation

For the purpose of investigating decline exponents for tight gas wells in Sulige gas field, we commence our study by using numerical simulation method. According to reservoir characterization and fracturing treatment of Sulige tight gas field, reservoir model can be simplified into rectangle reservoir models centered a fracture for the convenience of analysis (shown in Figure 1). Physical properties used in the reservoir model present in field characteristic. Gas wells produce at constant bottom flowing pressure.

Figure 1: The rectangle reservoir model.

Three cases are simulated at different formation permeability, with certain fracture half-length, fracture conductivity, and drainage area. Simulated gas rate data are shown in Figure 2. By using Equation (1) [13], decline exponent can be obtained by calculating simulated gas rate data (shown in Figure 3).

Figure 2: Rate vs. time with varied permeability.

Figure 3: Decline exponent vs. time with varied permeability.

As Figure 3 shows, for permeability less than 0.01 md, decline exponents will be bigger than 1 for a long time; while for permeability between 0.01 md and 0.1 md, decline exponents will be bigger than 1 in initial phase of production, smaller than 1 in late phase of production. However, decline exponents for traditional Arps decline method are between 0 and 1. Therefore, traditional Arps decline method will not be applicable.

4. Theory Basis

As the invalidation of Arps model, we turned to find new method. According to relevant reservoir knowledge, for narrow and striped tight formations, flow regime can be divided into three regimes, which are linear flow regime, elliptical flow regime, and pseudoradial flow, boundary-dominated flow. Here, linear flow is of importance for the new method.

Theoretically, there may be three kinds of linear flow, fracture-related linear flow [14], channel-related linear flow [15], Watterbargen linear flow [16], equations show in Table 1. Owing to the fracture length always smaller than channel width, fracture linear flow and channel linear flow will appear in most tight gas wells in Sulige tight gas field. Therefore, in this paper, we mainly highlight fracture linear flow and channel linear flow. From Table 1, we can see, any kind of linear flow, the inverse dimensionless rate vs. dimensionless square time is linear in coordinate system. Hence, it is possible to utilize the relationship to perform decline analysis. The following is detailed information.

Table 1: Working equation for different linear flow.

Linear flow
Working equation

Fracture linear flow

Channel linear flow

Wattenbarger linear flow

4.1. Fracture Linear Flow

For fracture linear flow, dimensionless for constant pressure will be given:

According to the definition of dimensionless variables, Equation (1) can be simplified to

where

4.2. Channel Linear Flow

For channel linear flow, dimensionless for constant pressure will be given:

Through dimensioning, Equation (5) can be simplified to

where

Equation (3) and Equation (6) show that and are linearly related during the fracture and channel linear flow.

The total volume produced from to can be found by integrating Equation (6):

Looking at Equation (7), for channel flow, the slope can be rewritten as

where and are, respectively, defined as:

can be treated as a constant independent of a particular well and independent of a particular reservoir, and is independent of particular well but dependent on its pressure difference, which impose minor effect on the linear relationship between and . Therefore, the variation in among different wells in the same field would be attributable to changes in permeability, thickness, and channel width. For a given tight gas reservoir, porosity varies little compared with permeability and thickness and can be treated as a constant. As a result, the differences in observed slopes for wells in a particular area can often be primarily attributed to variations in the term .

Equations (3) and (6) both embodies linear relationship in the form of vs. theoretically. Simulated production performance for long stripped formation verified that relationship. Figure 4 clearly displays these two linear relationships.

Figure 4: Cases with varied slope same intercept.

5. Practical Application

In order to validate the theory established, practical application is essential. Equations (6) and (8) were applied to predict and match gas well production performance from Sulige tight gas field. Well information is provided in Table 2. By plotting rate/time production data vs. , determining the and , then utilizes and to match rate/time in a coordinate system. Plots ( vs. ) for well A, well B, well C, well D, are presented in Figures 5–8, which exhibit straight-line behavior. Aided by straight-line, matching curves ( vs. ) for well A, well B, well C, and well D are presented in Figures 9–12. Figures 5–8 suggest that linear flow was dominant for well A, well B, and well C. In addition, taking D well for example, launching prediction for well performance in the future during linear flow is also plausible. For D well, linear relationship was generated only by year 1 and year 2, and year 3 production data is projected in Figures 8 and 12.

Table 2: Reservoir and rock data.

Well name
Pi (MPa)
(°C)
(mD)
(m)
(%)
Sw (%)

Well A
31.2
110
0.021
5.4
6.2
55

Well B
30.9
112.3
0.058
12.5
7.2
57

Well C
31.2
115.1
0.079
15.8
6.8
54

Well D
31.8
112.4
0.046
16.8
6.4
56.5

Figure 5: vs. matching plot for A well.

Figure 6: vs. matching plot for B well.

Figure 7: vs. matching plot for C well.

Figure 8: vs. matching plot for D well.

Figure 9: vs. matching plot for A well.

Figure 10: vs. matching plot for B well.

Figure 11: vs. matching plot for C well.

Figure 12: vs. matching plot for D well.

6. Sensitivity Analysis of and

Equation (8) can be employed to describe each gas well’s production performance during channel linear flow. Slope and intercept dominate gas rate decline tendency. To examine the effect of these two variables on production decline, a sensitivity analysis was performed. Figures 13~16 show the results. Figure 13 shows the vs. curves where is varied but is constant, Figure 14 corresponds vs. curves. Figure 15 shows the vs. curves where is varied but is constant; Figure 16 is corresponding vs. curves. Figure 13 shows that is indicative of the strength of well; a small predicts a better production performance. Figure 15 shows that intercept is only an indication of near-wellbore conditions.

Figure 13: Cases with varied slope same intercept.

Figure 14: Rate/time curve with varied slope same intercept.

Figure 15: Cases with varied intercept and same slope.

Figure 16: Rate/time curve with varied intercept and same slope.

7. Discussion

In this paper, we focus on taking advantage of linear behaviors in tight gas wells to predict well performance, and the results show well-suited. It is smart to bypass establishing a model considering variable decline exponent. Also, utilizing linear flow model to promote decline analysis is more convincing than empirical methods. However, the idea does not involve the circumstance when some gas wells enter into boundary-dominated flow. Combining the approach developed in this paper with Arps model will be a good idea. Besides, given the point of this paper is decline analysis, so distinguishing channel linear flow from Watterbarger linear flow is not discussed in detail. The reason we highlight channel linear flow in this paper is that the fracture length for most gas wells in Sulige tight gas field is smaller than the drainage boundary. Hence, Watterbarger linear flow can be excluded in most cases. The only possibility for long-term linear flow would be channel linear flow.

8. Summary and Conclusions

(1)An effective novel approach for predicting future gas rate and EUR has been developed for tight gas wells by using linear flow(2)The inverse of rate vs. the square of time is observed to fit a straight line in Sulige tight gas field and same gas field, which can forecast well performance accurately before pseudoradial flow(3)The slope obtained from fitting can be used to indicate the quality of gas formation(4)It is possible to match the entire production history for tight gas wells by combining the approach developed in this paper with Arps model

Nomenclature

: | Decline exponent, dimensionless
: | Decline rate, %
Di: | Initial decline rate, %
: | Defined in Equation (10)
: | Defined in Equation (11)
: | Total system compressibility, Pa⁻¹
: | Cumulative gas produced
: | Formation thickness, m
: | Formation permeability, md
: | Slope of inverse rate vs square root of time
: | Real gas pseud pressure change, Pa/s
: | Standard pressure, Pa
: | Flow rate, measured at standard conditions
: | Dimensionless rate, (when fracture linear flow), (when channel linear flow)
: | Inverse rate change caused by skin effects
: | Time, s or day
: | Dimensionless time, (when fracture linear flow), (when channel linear flow)
: | Reservoir temperature, K
: | Standard temperature, K
: | Fracture half-length, m
: | Channel width, m
: | Gas viscosity, pa·s
: | Initial Pressure, MPa
: | Effective Porosity, %
: | Initial Gas Saturation, %.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was financially supported by the National Science Technology Major Project (No. 2016ZX05047004003).

References

J. J. Arps, “Analysis of decline curves,” Transactions of the AIME, vol. 160, no. 1, pp. 228–247, 1945.

C. L. Kupchenko, B. W. Gault, and L. Mattar, “Tight gas production performance using decline curves,” in CIPC/SPE Gas Technology Symposium 2008 Joint Conference, Calgary Alberta, Canada, June 2008.

A. N. Duong, “An unconventional rate decline approach for tight and fracture-dominated gas wells,” in Canadian Unconventional Resources and International Petroleum Conference, Calgary, Alberta, Canada, October 2010.

D. Ilk, J. A. Rushing, A. D. Perego, and T. A. Blasingame, “Exponential vs. hyperbolic decline in tight gas sands -understanding the origin and implications for reserve estimates using Arps’ decline curves,” in SPE Annual Technical Conference and Exhibition, Denver, CO, USA, September 2008.

L. Matter and S. Moghadam, “Modified power law exponential decline for tight gas,” in Canadian International Petroleum Conference (CIPC), Calgary, Alberta, Canada, June 2009.

P. Valko, “Assigning value to stimulation in the Barnett shale-a simultaneous analysis of 7000 plus production histories and well completion records,” in SPE Hydraulic Fracturing Technology Conference, The Woodlands, TX, USA, January 2009.

K. Joshi and W. J. Lee, “Comparison of various deterministic forecasting technique in shale gas reserviors,” in SPE Hydraulic Fracturing Technology Conference, The Woodlands, TX, USA, February 2013.

S. Yu, “Best practice of using empirical methods for production forecast and EUR estimation in tight/shale gas reservoirs,” in SPE Unconventional Resources Conference Canada, Calgary, Alberta, Canada, November 2013.

D. B. Neal and M. A. Mian, “Early-time tight gas production forecasting technique improves reserves and reservoir description,” SPE Formation Evaluation, vol. 4, no. 1, pp. 25–32, 1989.

C. Dought and G. J. Moridis, “The use of the bimodal production decline curve for the analysis of hydraulically fractured shale/tight gas reservoirs,” in Proceedings of the 6th Unconventional Resources Technology Conference, Houston, TX, USA, July 2018.

M. Alem, T. Baig, and A. Muggeridge, “Predicting the performance of tight gas reservoirs,” in SPE Europec featured at 81st EAGE Conference and Exhibition, London, England, UK, June 2019.

F. K. Mienzan and D. Y. Asumadu, “Gas rate-decline type curve analysis using a semi-analytical model,” in SPE Nigeria Annual International Conference and Exhibition, Virtual, August 2020.

J. P. Spivey, J. H. Frantz Jr., J. R. Williamson, and W. K. Sawyer, “Applications of the transient hyperbolic exponent,” in SPE Rocky Mountain Petroleum Technology Conference, Keystone, CO, USA, May 2001.

T. Chen and R. Kry, “Pressure buildup from a vertically fractured low-permeability formation(includes associated paper 12938),” Journal of Petroleum Technology, vol. 34, no. 4, pp. 917–924, 1982.

L. Larsen and M. Hovdan, “Analyzing well test data from linear reservoirs by conventional methods,” in SPE Annual Technical Conference and Exhibition, Dallas, TX, USA, September 1987.

J. A. Arévalo-Villagrán, R. A. Wattenbarger, F. Samaniego-Verduzco, and T. T. Pham, “Production analysis of long-term linear flow in tight gas reservoirs: case histories,” in Paper SPE 71516 presented at the 2001 SPE Annual Technical Conference and Exhibition, New Orleans, LA, USA, September-October 2001.