Next Article in Journal
Optimizing Clinical Workflow Using Precision Medicine and Advanced Data Analytics
Previous Article in Journal
Development and Optimization of Chromatographic Conditions for the Determination of Selected B Vitamins in Pharmaceutical Products
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Well-Controlled Dynamic Reserves Calculation of Abnormally High-Pressured Gas Reservoirs Developed with Edge and Bottom Aquifers

1
School of Earth Science and Engineering, Xi’an Shiyou University, Xi’an 710065, China
2
School of Petroleum Engineering, Xi’an Shiyou University, Xi’an 710065, China
3
China National Oil and Gas Exploration and Development Company Ltd., Beijing 100034, China
4
Research Institute of Petroleum Exploration and Development, PetroChina, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Processes 2023, 11(3), 938; https://doi.org/10.3390/pr11030938
Submission received: 14 February 2023 / Revised: 9 March 2023 / Accepted: 16 March 2023 / Published: 19 March 2023
(This article belongs to the Section Energy Systems)

Abstract

:
A material balance equation (MBE) is usually employed to estimate well-controlled dynamic reserves, but a more accurate calculation of these reserves is crucial. For abnormally pressured reservoirs with aquifers, the abnormally pressure and aquifer are difficult to determine, which leads to inaccurate calculations of well-controlled dynamic reserves. In this work, the MBE was established. Then, the correction function and abnormally pressure effect were deduced to determine well-controlled dynamic reserves calculations for weak aquifers and strong aquifers, respectively. For a weak aquifer, a pressure correction was determined to exist at a pressure coefficient of 1.2, and then the correction function was established to calculate the dynamic reserves. For a strong aquifer, the MBE derived by Walsh was employed to calculate the dynamic reserves and water influx. The relationship between the abnormally pressure effect and the dimensionless formation pressure was proposed. Classical, abnormally high-pressured gas reservoirs, Anderson L and Amu Darya, were selected to validate the accuracy and applicability of the newly developed methods. The results show that an abnormally high-pressure effect has a great influence on the dynamic reserves calculation in the early stage and the water influx has a stronger influence on the dynamic reserves calculation in the late stage.

1. Introduction

Due to the effect of abnormally high pressure with bottom water influx, the relationship between p/Z and cumulative gas production volume is no longer a straight line as shown in normally pressured reservoirs, which renders the calculation of dynamic reserves difficult [1,2,3,4]. The MBE (material balance equation) has been indicated to be one of the most accurate calculation methods for well-controlled dynamic reserves of different types of gas reservoirs, whose basic principle is the conservation of matter [5,6,7]. Many material balance equations (MBEs) with different forms have been proposed by researchers to calculate the dynamic reserves involved with abnormally high-pressured or water influx issues [8,9,10,11,12].
Regarding abnormally high-pressure effects, Ramagost and Farshad (1981) considered the relationship between p/Z and cumulative gas production volume in the abnormally high-pressured phase to be a straight line. However, with a decrease in the extraction of gas and the formation pressure, the elastic expansion of gas and rock particles and the recompaction of reservoir rock are inevitable, which can supplement gas reservoir energy and reduce the formation pressure decline rate. Thus, the first straight line with a lower pressure decline is formed. As production continues downward to the normally pressured, a second straight line with a higher slope is formed. To solve the complex problem of calculation, certain scholars have proposed the form of a binomial, such as Chen and Hu (1993) [13] who proposed a dynamic reserves calculation equation in binomial form for constant volume, abnormally high-pressured reservoirs. Gonzales et al. (2008) [14] developed a binomial MBE for abnormally high-pressured reservoirs without aquifers in phases of cumulative gas production. An MBE for water-driven gas reservoirs considering the overall compressibility of rocks was established by Li (2008) [15]. The compressibility factor is expressed as a function of production characteristics, such as declining rate and gas productivity index, so only known production data are required for reserves calculation [16]. Based on the MBE of abnormally high-pressured reservoirs, a new MBE containing the squared phase of cumulative production is proposed, the parameters of which can be determined by applying production history data using binary regression analysis [17]. A new binomial MBE for predicting well-controlled dynamic reserves of gas reservoirs is proposed, exploiting the similarity between the curve shape of Ph-GpD and the IPR in abnormally high-pressured reservoirs [18]. Jiao et al. (2017) [19] proposed a new concept to describe the energy of elastic rock and water influx using the linear relationship of cumulative production, which simplifies the MBE and converts it to a multivariate nonlinear equation.
The effective compressibility of abnormally high-pressured reservoirs is related to the effective pore compressibility of rock and the compressibility and saturation of formation irreducible water, and it is difficult to accurately obtain these parameters. To solve these problems, Ramagost and Farshad (1981) [20] considered the influence of fluid and pore space to estimate the compressibility of formation rock (Cf) and immobile water (Cw); then, the effective compressibility factor of formation was determined. Pletcher (1981) [21] reorganized the MBE of a geo-pressured gas reservoir proposed by Ramagost and Farshad (1981) and derived a different linear expression, by which the original gas in place (OGIP) and the effective elastic driving force of rock or water can be simultaneously estimated. Fetkovich et al. (1998) [22] established a more complete MBE formula that can be applied to any type of reservoir. Yale et al. (1993) [23] established compressibility factor plates for unconsolidated cores, brittle cores, and strongly consolidated cores by core experiments. The empirical formula for calculating rock compressibility was established in abnormally high-pressured gas reservoirs [24]. Rahman et al. (2006) [25] proposed a rigorous MBE for geo-pressured hydrocarbon reservoirs using the origin expressions of compressibility of pore, water, and oil. The compressibility factor is no longer considered a constant but a function of pressure. A high-pressure experimental system was designed by the authors of [26], and a series of high-pressure compressibility factor tests of pure water, nitrogen, and rocks under different water saturations were carried out.
With the production of natural gas, the formation pressure declines and diffuses to the external natural waters via elastic expansion, which causes the elastic expansion of formation water and renders the calculation of dynamic reserves difficult. Therefore, to accurately calculate the water influx, certain scholars have proposed different forms of calculation methods [27]. For example, Liao (1990) [28] employed the unary, Lagrangian, three-point interpolation method and least squares curve to regress the dimensionless, natural water influx calculation formula of hemispherical flow, planar radial flow, and linear flow, thus greatly improving the calculation speed and accuracy of water influx. Fetkovich et al. (1998) [22] introduced a cumulative compressibility including pore compressibility, water compressibility, gas solubility, and any limited aquifer influx into the BME. Yang et al. (2008) [29] discussed the dynamic reserves of well control based on numerical simulation and verified the effect of the aquifer ratio on the dynamic reserves of water-driven gas reservoirs. Based on the principle of the material balance of water-driven gas reservoirs, a mathematical model of material balance was established, a new aquifer ratio parameter was introduced, and the relationship between the cumulative production of gas and that of water was obtained under well-controlled dynamic reserves and aquifer ratios that are unknown [30,31].
In this work, to solve the problems of abnormally high-pressure effects and water influx, the MBEs of abnormally high-pressured gas reservoirs with different aquifer size ratios are derived. For a weak aquifer ratio, by determining the inflection point of the curve of an abnormally high-pressure effect, a new calculation of well-controlled dynamic reserves has been established. With different aquifer size ratios, the total elastic expansion factor of the reservoir is accurately obtained by regression fitting the relationship between the dimensionless formation pressure (ppD) and the abnormally pressure effect (ĈeΔp). Next, the binomial equation for calculating well-controlled dynamic reserves and water influx is obtained. The accuracy of these methods is verified using the Anderson L gas reservoir, and the applicability and convenience of these methods are verified using the Amu Darya gas reservoir as an example.

2. Materials and Methods

2.1. MBE of Abnormally High-Pressured Reservoirs

The energy during the exploitation of reservoirs is mainly derived from the elastic expansion of natural gas, compaction of the reservoir and elastic expansion of rock particles, elastic expansion of formation immobile water, and the influx of an external aquifer. The MBE of the reservoir can be obtained by describing the above process in mathematical language:
G p B g = G ( B g B gi ) + V gi C f Δ p + V gi S wi C w Δ p + W e W p B w
where G is the initial dynamic reserves of the reservoir under standard conditions, 108 m3; Gp is cumulative gas production under standard conditions, 108 m3; Bgi is the volume factor of natural gas in the initial state, m3/m3; Vgi is the total pore volume of a gas field in the initial state, 108 m3; Swi is irreducible water saturation, dimensionless; Cf and Cw are the elastic factor of formation rock and water, MPa−1; Δp is formation pressure decline, MPa; We is cumulative natural water influx, m3; Wp is the cumulative water yield, m3; and Bw is the water volume factor of formation, m3/m3.
Considering Δp = pi − p and Vgi = GBgi/(1 − Swi), rearranging Equation (1) yields:
G = G p B g ( B g B gi ) + ( C w S wi + C f 1 S wi ) B gi Δ p + ( W e W p B w V gi ) B gi
Define the formation effective compressibility factor Ce and water influx volume factor ω:
C e = C w S wi + C f 1 S wi
ω = W e W p B w V gi B gi
According to the expression of the natural gas volume factor:
B g = p sc Z T p T sc
B gi = p sc Z i T p i T sc
where Z and Zi represent the natural gas deviation factor under the current reservoir conditions and original reservoir conditions, respectively, dimensionless.
Introducing Equation (2) to Equation (6) into Equation (1) and organizing it yields:
p Z ( 1 C e Δ p ω ) = p i Z i ( 1 G p G )
As shown in Equation (7), the difference between the pressure decline formula of abnormally high-pressured reservoirs and that of constant-volume gas reservoirs is that the former takes into account the impact of the abnormally high-pressure effect CeΔp and water influx of reservoir ω. Therefore, the determination of the above two parameters is critical to the well-controlled dynamic reserves estimation of reservoirs by the pressure decline method.
In this work, the aquifer size ratio, which is defined as the ratio of the volume of aquifer to the total gas in place at the reservoir condition, was introduced to depict the strength of aquifer energy. Herein, an aquifer ratio Mwg below 5 is classified as a weak aquifer, while a strong aquifer is classified as an aquifer ratio above 5 (Mwg ≥ 5).

2.2. Calculation of Dynamic Reserves with a Weak Aquifer Ratio

When the reservoirs have a weak aquifer, the calculation method of the Chatas water influx is approximately valid:
j Q ( Δ t D j ) Δ p e j Q D 0 · ( p i p )
Q D 0 = 1 2 [ ( r aq r ws ) 2 1 ]
where raq is the outer edge radius of natural water, m.
Herein:
W e = π ( r aq 2 r ws 2 ) h w ϕ C e ( p i p ) = V aq C e ( p i p )
where Vaq is the pore volume in natural waters, m3.
The abnormally high-pressure stage is generally in the middle and early stages of development. For weak aquifers, most reservoirs still do not produce formation water (Wp = 0). Herein, introducing Equation (9) into Equation (7) yields:
p Z ( 1 C ^ e Δ p ) = p i Z i ( 1 G p G )
where C ^ e = ( 1 + M wg ) C e is deemed the effective compressibility factor of an abnormally high-pressured reservoir, MPa−1.
Therefore, the parameters of the abnormally high-pressure effect and aquifer size ratio can be combined, and Equation (10) has the same form as the pressure decline equation of the reservoir.

2.2.1. Establishment of the Reserves Estimation Equation

The pressure decline curve of an abnormally high–pressured, confined gas reservoir can be simplified into two straight sections, as shown in Figure 1. Section I and section II represent the two stages of abnormally high-pressured and normal pressure, respectively, in the development process, and the inflection point of the two straight lines is recorded as C.
According to the MBE, describing straight-line Section I yields:
p C Z C = p i Z i ( 1 G pC G pes )
where pC is the formation pressure corresponding to the inflection point, MPa; Zc is the gas deviation factor under the condition of the inflection point, dimensionless; Gpc is the cumulative gas recovery corresponding to the formation pressure at the inflection point, 108 m3; and Gpes is the virtual dynamic reserves of the reservoir, 108 m3.
The pressure decline equation of the reservoir can be written as:
p C Z C [ 1 C ^ e ( p i p C ) ] = p i Z i ( 1 G pC G real )
where Greal is the real dynamic reserves of the reservoir, 108 m3.
The reserves correction equation of abnormally high-pressured gas reservoirs can be obtained by combining Equations (11) and (12). When calculating well-controlled dynamic reserves of the reservoir, it is only necessary to determine the virtual dynamic reserves using the data of linear section Ⅰ and then combine it with Equation (13) to obtain the real reserves:
G real = G pes 1 + C ^ e ( p i p C ) p i / Z i p C / Z C 1

2.2.2. Determination of Inflection Point Position

The relationship between the abnormally pressure effect and the pressure coefficient is obtained by rearranging the pressure decline formula of the reservoir:
C ^ e Δ p = 1 ( 1 G p G ) p / Z p i / Z i
In order to analyze the variation rule of the abnormal high-pressure effect Ĉp, five typical abnormally high-pressured gas reservoirs N. Ossun, Louis. Offshore, Stafford, GOM, and Miocene S. Louis were taken as examples (Table 1). ĈeΔp values at different pressure stages were calculated according to Equation (14). The relationship between Ĉp and pressure coefficient (α) was analyzed to determine the variation rule of abnormally high-pressure effect (Figure 2). The pressure coefficient (α) is defined in Equation (15). The historical production data of the five typical abnormally high-pressured gas reservoirs is shown in Appendix A.
α = p i / ( C · D )
where C is the hydrostatic pressure gradient, MPa/m; D is the central depth of reservoirs, m.
The relationship between ĈeΔp and α shows two different regular patterns in Figure 2: when the pressure coefficient was greater than 1.2, ĈeΔp significantly changed with α and showed a good linear negative correlation (refer to Table 2 for the regression results); when the pressure coefficient was less than 1.2, the change in ĈeΔp with α was slow. It can be seen that α = 1.2 should be taken as the cutoff point (position of inflection point C) between the abnormally high-pressured stage and the normal pressure stage in the pressure decline diagram. The gas reservoirs mentioned in this paper all had pressure coefficients greater than 1.2, so they are all exceptionally high-pressured reservoirs.
After determining the position of the inflection point α = 1.2 by replacing the pressure coefficient and deviation coefficient of inflection point C in Equation (13) (recorded as P1.2 and Z1.2, respectively), a new reserve estimation formula of the pressure decline method for abnormally high-pressured reservoirs can be obtained as Equation (16):
G real = G pes 1 + C ^ e ( p i p 1.2 ) p i / Z i p 1.2 / Z 1.2 1
where Greal is the real dynamic reserves corresponding to the inflection point α = 1.2, 108 m3.

2.3. Calculation of Dynamic Reserves with a Strong Aquifer Ratio

In this paper, the early stage of gas reservoir development is defined as the point when the recovery degree of reservoirs is less than 10%, and the late stage is defined as the point when the recovery degree of reservoirs is more than 50%. When reservoirs have a strong aquifer (Mwg > 5), the water influx process of reservoirs cannot be simplified according to the previously described method. Herein, the MBE will have the phases of the abnormally high-pressure effect and water influx. To accurately calculate well-controlled dynamic reserves, first, one of the variables needs to be determined.
Walsh (1996) gave a new form of the MBE of the reservoir, considering water influx (Equation (17)), which has the advantage of being able to regress both dynamic reserves and water influx parameters.
F = G E t + W e
Inside,
F = G p B g + W p B w
E t = E g + E w [ B gi ( S wi + M wg ) B wi ( 1 S wi ) ] + E f [ B gi ( 1 + M wg ) 1 S wi ]
E g = B g B gi
E w = B w B wi = C w ( p i p )
E f = V gi V g V gi = C f ( p i p )
where F is the volume of accumulated natural gas and formation water in the state of the gas reservoir, m3; Et is the total elastic expansion factor of the reservoir, m3/m3; Eg is the elastic expansion factor of natural gas, m3/m3; Ew is the elastic expansion factor of formation water, m3/m3; Ef is the elastic expansion factor of rock, m3/m3; and Bw is the volume factor of formation water, m3/m3. The subscript i indicates the initial state.
The total elastic expansion factor Et of the reservoir is accurately calculated as follows:
E t = E g + E w [ B gi ( S wi + M wg ) B wi ( 1 S wi ) ] + E f [ B gi ( 1 + M wg ) 1 S wi ] = B g B gi + B gi ( p i p ) 1 S wi [ C w ( S wi + M wg ) + C f ( 1 + M wg ) ] = p sc Z T p T sc [ 1 p / Z p i / Z i ( 1 C ^ e Δ p ) ]
According to Equation (23), the key to obtain Et is still the determination of abnormally high-pressured term (ĈeΔp), and the recompaction of rocks and the elastic expansion of rock particles will occur in reservoirs, and the elastic expansion of bound water in reservoirs will also occur, leading to the change of effective elastic compression coefficient (Ce). In order to solve this problem, the dimensionless apparent formation pressure (ppD) is defined in Equation (24).
p pD = p / Z p i / Z i
Subsequently, the numerical simulation model of abnormally high-pressured gas reservoirs is established; the compressibility factor of rock and formation water is input through Equations (21) and (22), respectively; and the water influx is defined as unsteady by the Carter Tracy aquifer method [32]. After calculation, the data under different aquifer size ratios were obtained (Figure 3). The relationship between ĈeΔp and ppD can be obtained by fitting the regression of Cases 1 to 4 (Equation (25)). After ĈeΔp is determined, Et can be further calculated, and then the well-controlled dynamic reserves of the reservoir can be calculated in Equation (17).
C e ^ Δ p = 3.8475 p pD 2 7.7496 p pD + 3.9019

3. Results and Discussion

3.1. Verification against Field Case

The classic Anderson L reservoir was selected to verify the reliability of the method. The reservoir space is mainly composed of secondary pores with a few fractures, the structural types of it include complex faults and anticlinal structures. It is a typical weakly elastic water-drive abnormally high-pressured condensate gas reservoir. The initial formation pressure of the reservoir is 65.548 MPa, the burial depth is 3403.7 m, the irreducible water saturation of the formation is 0.35, and the compressibility of formation rock and irreducible water are 2.176 × 10−3 MPa−1 and 4.410 × 10−3 MPa−1, respectively. The gas production data are shown in Table 3.
To compare the characteristics and reliability of different methods, we took the calculation results of the volumetric method as the true value of calculation, and compared the reserves calculation methods in this paper and the literature with the volumetric method to calculate the relative errors. Table 4 lists part of the calculation results.
A comparison of the calculation results in the Table 4 reveals the following findings:
  • The errors of the methods that can be applied to the early dynamic reserves calculation in the literature were too large, which proves the necessity of correction. The error of the calculation results of Equations (16) and (25) was 5.59% and 4.88%, respectively; therefore, Equations (16) and (25) were feasible;
  • Among all the methods that can be used in the early dynamic reserves calculations, Equation (16) had the highest accuracy with a weak aquifer, and Equation (25) had a higher accuracy with a strong aquifer. More importantly, these methods are simple in form and convenient in calculation, which has stronger practicality.

3.2. Validation against Case

In this section, the proposed new dynamic reserves evaluation method is validated by B-15 well in the Amu Darya gas reservoir in Turkmenistan as an example. The validation case was taken from the Callovian-Oxfordian Stage of the Middle and Upper Jurassic carbonate gas field on the right bank of the Amu Darya, Turkmenistan. The reservoir space is mainly composed of secondary pores and fractures which are extensively developed. The initial formation pressure of well B-15 is 54.34 MPa; the burial depth of the reservoir is 3100 m; the initial pressure coefficient is 1.80; and the aquifer ratio is 13. Presently, the formation pressure is 46.557 MPa, and the pressure coefficient is 1.53. The production history data for well B-15 since production are shown in Table 5 and the current degree of recovery is close to 70%.
Bourgoyne (1990) [33] proposed a method for calculating binomial reserves. Through fitting results, it was found that the method had multiple solutions, and the error between two different formulas could not be ignored (Figure 4). The regression factors were a1 = 0.6528, b1 = 0.0371, a2 = 0.6144, and b2 = 0.0466; the geological reserves were G1 = 26.857 × 108 m3 and G2 = 24.935 × 108 m3. This finding shows that the method does suffer from a high degree of ambiguity when applied to early reservoir development.
The pseudo reserves correction method was utilized for the calculation, and the virtual dynamic reserves were calculated as Gpes = 55.007 × 108 m3 according to the fitting of Equation (11) (Figure 5), and the final real dynamic reserves were Greal-a = 25.766 × 108 m3 after correction according to Equation (16). We used abnormally high-pressured reservoirs for the whole period production data to verify the accuracy of Equation (16) and the real dynamic reserves were Greal-b = 25.473 × 108 m3. Then, we found that the method used linear extrapolation, it was less likely to cause error accumulation, and it was suitable for the calculation of early dynamic reserves.
The relationship curve between F and Et was drawn (Figure 6) according to Equation (17), fitting the gas field dynamic reserves G = 24.254 × 108 m3; the current cumulative water influx was We = 1.2 × 105 m3.

3.3. Discussion

In this paper, we showed the different calculation results by Equation (11), (16), (17), and Bourgoyne’s Model (Figure 7). Among them, Bourgoyne’s Model showed ambiguity in the early dynamic reserves calculation, and the error of each solution could not be ignored. Therefore, Bourgoyne’s Model is not suitable for the dynamic reserves calculation from mid to late development. However, both Equation (16) and Equation (17) had preferable accuracy. When calculating the dynamic reserves of gas reservoirs, Equation (16) mainly considered the influence of the abnormally high-pressure effect (ĈeΔp), and Equation (17) mainly considered the influence of water influx. Therefore, Equation (16) is recommended when the abnormal high-pressure effect dominates in the early stage, and Equation (17) is recommended when the edge and bottom water invade in the late stage. When calculating dynamic gas reservoir reserves with different aquifer ratios, we recommend Equation (16) for a weak aquifer ratio and recommend Equation (17) for a strong aquifer ratio.
When the abnormally high-pressure effect is not considered, the calculation results in Figure 7 showed that Equation (11) was 52.9% and 53.2% above the Bourgoyne’s Model and Equation (16), respectively, indicating that the abnormally high-pressure effect cannot be disregarded. With the different aquifer size ratios, Equation (17) was 6.8% and 6.2% below Bourgoyne’s Model, and Equation (16), respectively; therefore, the aquifer size ratio cannot be underestimated in well-controlled dynamic reserves calculations.
In conclusion, the influence of the abnormally high-pressure effect should be considered in the early calculation of gas reservoir dynamic reserves in the development process of abnormally high-pressured edge and bottom water reservoirs. The influence of water influx should be considered when calculating the dynamic reserves of gas reservoirs.

4. Conclusions

(1)
A new reserves estimation method was proposed for dynamic reserves estimation of abnormally high-pressured reservoirs with edge and bottom water influx. The problem of an abnormally pressure effect was addressed by solving the inflection point problem of the traditional virtual reserves correction method. The problem of water influx was resolved by fitting the empirical formulas between ĈeΔp and ppD, and the formulas for the dynamic reserves and water influx of the rearranged gas reservoir were obtained. These methods are simpler in form and more applicable;
(2)
With a weak aquifer ratio, the boundary point between the abnormally high-pressured stage and normally pressure stage of edge and bottom water reservoir development was the pressure coefficient, which was equal to 1.2. In the abnormally high-pressured stage, the pressure correction term ĈeΔp was approximately linearly negatively related to the corresponding pressure coefficient α;
(3)
Field cases showed that the calculated reserves by the newly proposed reserves estimation method were in accordance with actual reserves estimation methods. In the field case in Anderson L, the relative errors of the two new methods of dynamic reserves calculation were 5.59% and 4.88%, which demonstrated that the dynamic reserves estimation method was effective in abnormally high-pressured reservoirs;
(4)
In the field case of the Amu Darya B reservoir, the influence of the abnormally pressure effect on the dynamic reserves calculation was 52.9% and 53.2%, which mainly affected the dynamic reserves calculation of the early stage. The effect of the water influx on the dynamic reserves calculation was 6.8% and 6.2%, which mainly affected the dynamic reserves calculation of the late stage.

Author Contributions

Conceptualization, Y.C. and X.L.; Formal analysis, Y.C. and X.L.; Funding acquisition, Y.C. and C.T.; Investigation, C.G. and T.L.; Methodology, Y.C. and C.G.; Project administration, Y.C. and T.L.; Software development, X.L.; Visualization, Y.C. and C.G.; Writing-original draft Y.C.; Writing—review and editing Y.C. All authors have read and agreed to the published version of the manuscript.

Funding

The authors would like to appreciate the PetroChina Technological Research Project (2021DJ3301); Scientific Research Project of the Shaanxi Provincial Department of Education, China (20JK0848).

Institutional Review Board Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Appendix A

Table A1. Historical production data of the N. Ossun reservoir [33].
Table A1. Historical production data of the N. Ossun reservoir [33].
Numberp/MPaZGp/108 m3α
161.521.4730.001.65
261.001.4650.191.63
357.391.4000.931.54
451.151.2882.941.37
547.441.2194.411.27
641.821.1306.791.12
737.861.0757.891.01
832.970.9679.540.88
928.30.88711.360.76
Table A2. Historical production data of the Louis. Offshore reservoir [20].
Table A2. Historical production data of the Louis. Offshore reservoir [20].
Numberp/MPaZGp/108 m3α
178.921.4960.001.99
273.611.4382.811.85
369.861.3978.101.76
463.811.33015.181.61
559.131.28021.991.49
654.521.23028.721.37
750.891.19234.081.28
847.221.15441.061.19
944.051.12245.491.11
1040.181.08451.631.01
1137.301.05755.990.94
Table A3. Historical production data of the Stafford reservoir [20].
Table A3. Historical production data of the Stafford reservoir [20].
Numberp/MPaZGp/108 m3α
149.651.1840.001.83
248.101.1670.201.78
346.351.1480.431.71
444.341.1260.701.64
542.981.1120.941.59
643.061.1131.041.59
740.941.0911.211.51
839.441.0761.471.46
936.871.0501.721.36
1031.180.9992.351.15
1125.320.9563.010.93
1221.490.9353.470.79
1319.550.9283.730.72
Table A4. Historical production data of the GOM reservoir [33].
Table A4. Historical production data of the GOM reservoir [33].
Numberp/MPaZGp/108 m3α
184.481.68690.002.18
281.031.63940.112.09
376.131.57190.241.96
471.721.51080.381.85
566.891.4440.521.72
665.101.41920.651.68
761.581.37050.771.59
859.651.34390.931.54
956.131.32681.081.45
1052.551.27841.231.35
1149.101.23291.381.26
1246.691.20171.511.20
1341.721.1391.701.07
Table A5. Historical production data of the Miocene S. Louis reservoir [33].
Table A5. Historical production data of the Miocene S. Louis reservoir [33].
Numberp/MPaZGp/108 m3α
175.751.4890.002.27
257.331.2720.821.72
348.711.1741.551.46
443.101.1701.771.29
533.961.0782.121.02

References

  1. Oscar, M.O.; Fernando, S.V.; Jorge, A.V.; Arévalo, V.J.A. Advances in the production mechanism diagnosis of gas reservoirs through material balance studies. In Proceedings of the SPE Annual Technical Conference and Exhibition, Houston, TX, USA, 26–29 September 2004. [Google Scholar] [CrossRef]
  2. Zhang, L.X.; He, Y.X.; Guo, C.Q.; Yang, Y. Dynamic material balance method for estimating gas in place of abnormally high-pressured gas reservoirs. Lithosphere 2021, 2021, 6669012. [Google Scholar] [CrossRef]
  3. Lyu, X.; Zhu, G.; Liu, Z. Well-controlled dynamic hydrocarbon reserves calculation of fracture–cavity karst carbonate reservoirs based on production data analysis. J. Petrol. Explor. Prod. Technol. 2020, 10, 2401–2410. [Google Scholar] [CrossRef] [Green Version]
  4. Zhang, Y.; Ge, H.; Zhao, K. Simulation of pressure response resulted from non-uniform fracture network communication and its application to interwell-fracturing interference in shale oil reservoirs. Geomech. Geophys. Geo-Energy Geo-Resour. 2022, 8, 114. [Google Scholar] [CrossRef]
  5. Johnson, N.L.; Currie, S.M.; Ilk, D.; Blasingame, T.A. A simple methodology for direct estimation of gas-in-place and reserves using rate-time data. In Proceedings of the SPE Rocky Mountain Petroleum Technology Conference, Denver, CO, USA, 14–16 April 2009. [Google Scholar] [CrossRef] [Green Version]
  6. Kabir, C.S.; Parekh, B.; Mustafa, M.A. Material-balance analysis of gas reservoirs with diverse drive mechanisms. In Proceedings of the SPE Annual Technical Conference and Exhibition, Society of Petroleum Engineers, Houston, TX, USA, 28–30 September 2015. [Google Scholar] [CrossRef]
  7. Miao, Y.; Zhao, C.; Zhou, G. New rate-decline forecast approach for lowpermeability gas reservoirs with hydraulic fracturing treatments. J. Petrol. Sci. Eng. 2020, 190, 107112. [Google Scholar] [CrossRef]
  8. King, G.R. Material-balance techniques for coal-seam and devonian shale gas reservoirs with limited water Influx. SPE Reserv. Eng. 1993, 8, 67–72. [Google Scholar] [CrossRef]
  9. Walsh, M.P. Discussion of Application of material balance for high pressure gas reservoirs. In Proceedings of the SPE Gas Technology Symposium, Calgary, AB, Canada, 28 April–1 May 1996. [Google Scholar] [CrossRef]
  10. Pletcher, J.L. Improvements to reservoir material-balance methods. SPE Reserv. Eval. Eng. 2002, 5, 49–59. [Google Scholar] [CrossRef] [Green Version]
  11. Sun, Z.; Shi, J.; Zhang, T.; Wu, K.; Miao, Y.; Feng, D.; Sun, F.; Han, S.; Wang, S.; Hou, C.; et al. The modified gas-water two phase version flowing material balance equation for low permeability CBM reservoirs. J. Petrol. Sci. Eng. 2018, 165, 726–735. [Google Scholar] [CrossRef]
  12. Zhang, Y.J.; Ge, H.K.; Shen, Y.H. Evaluating the potential for oil recovery by imbibition and time-delay effect in tight reservoirs during shut-in. J. Pet. Sci. Eng. 2019, 184, 106557. [Google Scholar] [CrossRef]
  13. Chen, Y.Q.; Hu, J.G. A new method for determining geological reserves and effective compression factors in abnormally high-pressured gas reservoirs. Gas Ind. 1993, 13, 53–58. [Google Scholar]
  14. Gonzales, F.E.; Ilk, D.; Blasingame, T.A. A quadratic cumulative production model for the material balance of an abnormally pressured gas reservoir. In Proceedings of the SPE Western Regional and Pacific Section AAPG Joint Meeting, Bakersfield, CA, USA, 31 March–2 April 2008. [Google Scholar] [CrossRef] [Green Version]
  15. Li, C.L. Study on the calculation method of water intrusion in gas reservoirs. Xinjiang Pet. Geol. 2003, 430–431. [Google Scholar]
  16. Feng, G.Q.; Liu, C.D.; Ran, L.; Xiao, M.H. Comparison of reserves calculation methods for abnormally high-pressured gas reservoirs. J. Chongqing Univ. Sci. Technol. (Nat. Sci. Ed.) 2010, 12, 13–16. [Google Scholar] [CrossRef]
  17. Ding, X.F.; Liu, Z.B.; Pan, D.Z. New method for calculating geological reserves and cumulative effective compression factor of abnormally high-pressured gas reservoirs. Acta Pet. Sin. 2010, 31, 626–628, 632. [Google Scholar]
  18. Bian, X.Q.; Du, Z.M.; Chen, J. Binomial material balance equation predicts reserves in abnormally high-pressured gas reservoirs. J. Southwest Pet. Univ. (Sci. Technol. Ed.) 2010, 32, 75–78. [Google Scholar] [CrossRef]
  19. Jiao, Y.W.; Xia, J.; Liu, P.C.; Zhang, J.; Li, B.Z.; Tian, Q.; Wu, Y.P. New material balance analysis method for abnormally high-pressured gas-hydrocarbon reservoir with water influx. Int. J. Hydrogen Energy 2017, 42, 18718–18727. [Google Scholar] [CrossRef]
  20. Ramagost, B.P.; Farshad, F.F. P/Z abnormally pressured gas reservoirs. In Proceedings of the SPE Annual Technical Conference and Exhibition, San Antonio, TX, USA, 4–7 October 1981. [Google Scholar] [CrossRef]
  21. Roach, R.H. Analyzing Geopressured Reservoirs—A Material Balance Technique; Unsolicited paper; SPE: Richardson, TX, USA, December 1981. [Google Scholar]
  22. Fetkovich, M.J.; Reese, D.E.; Whitson, C.H. Application of a general material balances for high-pressure gas reservoirs. SPE J. 1991, 3, 3–13. [Google Scholar] [CrossRef] [Green Version]
  23. Yale, D.P.; Nabor, G.W.; Russell, J.A. Application of variable formation compressibility for improved reservoir analysis. In Proceedings of the SPE Annual Technical Conference and Exhibition, Houston, TX, USA, 3–6 October 1993. [Google Scholar] [CrossRef]
  24. Wang, S.W.; Stevenson, V.M.; Ohaeri, D.H. Analysis of overpressured reservoirs with a new material balance method. In Proceedings of the SPE Annual Technical Conference and Exhibition, Houston, TX, USA, 3–6 October 1999. [Google Scholar] [CrossRef]
  25. Rahman, N.M.A.; Anderson, D.M.; Mattar, L. New, rigorous material balance equation for gas flow in a compressible formation with residual fluid saturation. In Proceedings of the SPE Gas Technology Symposium, Calgary, AB, Canada, 15–17 May 2006. [Google Scholar] [CrossRef]
  26. Zhang, J.Z.; Gao, T.R. Compressibility of abnormally-pressured gas reservoirs and its effect on reserves. ACS Omega 2021, 6, 26221–26230. [Google Scholar] [CrossRef] [PubMed]
  27. Zhang, Y.; Zou, Y.; Zhang, Y. Experimental Study on Characteristics and Mechanisms of Matrix Pressure Transmission Near the Fracture Surface During Post-Fracturing Shut-In in Tight Oil Reservoirs. J. Pet. Sci. Eng. 2022, 219, 111133. [Google Scholar] [CrossRef]
  28. Liao, Y.T. Regression equation for calculating natural water influx. Oil Explor. Dev. 1990, 1, 71–74. [Google Scholar]
  29. Yang, E.L.; Zhang, J.G.; Song, K.P.; Huang, B. Numerical simulation study of low permeability gas reservoir at Daqing five stations. Gas Ind. 2008, 3, 102–104, 147–148. [Google Scholar] [CrossRef]
  30. Wang, Z.H.; Wang, Z.D.; Zeng, F.H.; Guo, P.; Xu, X.Y.; Deng, D. The material balance equation for fractured vuggy gas reservoirs with bottom water-drive combining stress and gravity effects. J. Nat. Gas Sci. Eng. 2017, 44, 96–108. [Google Scholar] [CrossRef]
  31. Xue, T.; Ji, Y.; Cheng, L.B.; Lv, C.S.; Yang, Y.X.; Lan, Z.K. A new method for determining dynamic storage capacity and water energy in water-driven gas reservoirs. In Petroleum Geology & Oilfield Development in Daqing; CNKI, Tongfang Knowledge Network Technology Co., Ltd.: Beijing, China, 2022; pp. 1–6. [Google Scholar] [CrossRef]
  32. Carter, R.D.; Tracy, G.W. An improved method for calculating water influx. Trans. AIME 1960, 219, 415–417. [Google Scholar] [CrossRef]
  33. Bourgoyne, A.T., Jr. Shale Water as a Pressure Support Mechanism in Gas Reservoirs Having Abnormal Formation Pressure. J. Pet. Sci. 1990, 3, 305. [Google Scholar] [CrossRef]
Figure 1. The relationship between p/Z and cumulative gas production volume of the abnormally high-pressured reservoirs.
Figure 1. The relationship between p/Z and cumulative gas production volume of the abnormally high-pressured reservoirs.
Processes 11 00938 g001
Figure 2. Relationship between the abnormally high-pressure effect and the pressure coefficient.
Figure 2. Relationship between the abnormally high-pressure effect and the pressure coefficient.
Processes 11 00938 g002
Figure 3. Fitting diagram of the relationship between gas fields ĈeΔp and ppD.
Figure 3. Fitting diagram of the relationship between gas fields ĈeΔp and ppD.
Processes 11 00938 g003
Figure 4. Pressure decline diagram of the methods of Bourgoyne (1990).
Figure 4. Pressure decline diagram of the methods of Bourgoyne (1990).
Processes 11 00938 g004
Figure 5. Pressure decline diagram of the pseudo reserves correction method in well B-15. (a). Early historical production data; (b) Historical production data for the whole period.
Figure 5. Pressure decline diagram of the pseudo reserves correction method in well B-15. (a). Early historical production data; (b) Historical production data for the whole period.
Processes 11 00938 g005
Figure 6. Well-controlled dynamic reserves estimation fitting diagram of well B-15.
Figure 6. Well-controlled dynamic reserves estimation fitting diagram of well B-15.
Processes 11 00938 g006
Figure 7. Comparison diagram of dynamic reserves obtained by three kinds of methods.
Figure 7. Comparison diagram of dynamic reserves obtained by three kinds of methods.
Processes 11 00938 g007
Table 1. Initial production data of the five typical abnormally high-pressured gas reservoirs.
Table 1. Initial production data of the five typical abnormally high-pressured gas reservoirs.
Gas FieldDepth/mInitial Pressure/MPaInitial Pressure Coefficient (αi)
N. Ossun3810.061.521.65
Louis. Offshore4053.878.921.99
Stafford2764.549.651.83
GOM3962.484.482.18
Miocene S. Louis.3962.475.742.27
Table 2. Gas reservoir relational regression equation of ĈeΔp and α (α ≥ 1.2).
Table 2. Gas reservoir relational regression equation of ĈeΔp and α (α ≥ 1.2).
Gas FieldRegression EquationCorrelation Coefficient (R2)
N. OssunĈeΔp = −0.1600α + 0.26090.9577
Louis. OffshoreĈeΔp = −0.1706α + 0.30350.9789
StaffordĈeΔp = −0.2871α + 0.52730.9565
GOMĈeΔp = −0.1900α + 0.40330.9777
Miocene S. Louis.ĈeΔp = −0.2754α + 0.61760.9175
Table 3. Historical production data of the Anderson L gas reservoir.
Table 3. Historical production data of the Anderson L gas reservoir.
Numberp/MPaZGp/108 m3α
165.5481.4400.0001.964
264.0661.4180.1111.920
361.8461.3870.4651.854
459.2601.3440.9131.776
557.4471.3161.2061.722
655.2201.2821.5581.655
752.4211.2392.1351.571
851.0631.2182.4781.530
948.2771.1762.9761.447
1046.3401.1473.3301.389
1145.0571.1273.6221.350
1239.7411.0484.8881.191
1332.8600.9776.4820.985
1429.6130.9287.9700.887
1525.8550.8919.2220.775
1622.3870.85410.4260.671
Table 4. Comparison of dynamic reserves estimation results of different methods.
Table 4. Comparison of dynamic reserves estimation results of different methods.
MethodsReserves/108 m3Relative Error/%Apply to Early Stage
This articleEquation (16)21.0355.59Yes
Equation (25)20.8944.88Yes
LiteratureBourgoyne (1990)21.0725.87No
Yale (1993)21.5218.03Yes
Chen (1993)21.1246.03No
Gan (2001)21.3797.31No
Oscar (2004)21.6898.87No
Elsharkawy (1996)21.8049.45Yes
Volumetric methodBian19.922//
Table 5. Historical production data of well B-15.
Table 5. Historical production data of well B-15.
Numberp/MPaGp/108 m3Zα
154.3490.0001.2761.80
254.3030.01061.2571.79
351.1521.4071.2191.68
448.4053.0151.1761.59
546.5574.3811.1471.53
643.7256.2051.1351.45
740.5018.7521.1211.34
837.12710.5581.1081.32
935.14613.1061.0911.18
1033.17916.7281.0801.06
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cheng, Y.; Luo, X.; Lu, T.; Guo, C.; Tan, C. Well-Controlled Dynamic Reserves Calculation of Abnormally High-Pressured Gas Reservoirs Developed with Edge and Bottom Aquifers. Processes 2023, 11, 938. https://doi.org/10.3390/pr11030938

AMA Style

Cheng Y, Luo X, Lu T, Guo C, Tan C. Well-Controlled Dynamic Reserves Calculation of Abnormally High-Pressured Gas Reservoirs Developed with Edge and Bottom Aquifers. Processes. 2023; 11(3):938. https://doi.org/10.3390/pr11030938

Chicago/Turabian Style

Cheng, Youyou, Xiang Luo, Tianyuan Lu, Chunqiu Guo, and Chengqian Tan. 2023. "Well-Controlled Dynamic Reserves Calculation of Abnormally High-Pressured Gas Reservoirs Developed with Edge and Bottom Aquifers" Processes 11, no. 3: 938. https://doi.org/10.3390/pr11030938

APA Style

Cheng, Y., Luo, X., Lu, T., Guo, C., & Tan, C. (2023). Well-Controlled Dynamic Reserves Calculation of Abnormally High-Pressured Gas Reservoirs Developed with Edge and Bottom Aquifers. Processes, 11(3), 938. https://doi.org/10.3390/pr11030938

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