# Reserve Estimation from Early Time Production Data in Geopressured Gas Reservoir: Gas Production of Cumulative Unit Pressure Drop Method

Authors: Hongfeng Wang; Xiaoping Li; Hedong Sun; Guangren Meng; Wen Cao; Songbai Zhu
Journal: Geofluids
 (2021)

---

## Abstract

There is high uncertainty in reserve estimation during the early development of deep ultrahigh pressure gas reservoirs, largely because it remains challenging in accurately determining the formation compressibility. To overcome this, starting from the definition of compressibility, a novel gas production of cumulative unit pressure drop analysis method was established, of which the effectiveness was proven by applications in calculating the reserves of three gas reservoirs. It has been found that, in the limiting case, i.e., when the formation pressure dropped to the normal atmospheric pressure, the dimensionless gas production of the cumulative unit pressure drop was the reciprocal of the initial formation pressure. Besides, the relationship curve of the dimensionless gas production of the cumulative unit pressure drop and pressure drop was a straight line in the medium term, extending the straight line and intersecting the vertical line passing through the original formation pressure point, and the reserves can be determined according to the intersection point and the initial formation pressure. However, due to the influence of natural gas properties, the value needs further correction, and the correction coefficient depends on the pseudocritical temperature of natural gas. Specifically, when the pseudocritical temperature is given, the correction coefficient would be close to the minimum value of the natural gas deviation factor. When the pseudocritical temperature is more than 1.9 and less than 3.0, the minimum deviation factor would be between 0.90 and 1.0, and the higher the pseudocritical temperature, the closer the ratio is to 1.0.


---

## Body

## 1. Introduction

The proportion of natural gas in primary energy consumption has been continuously increasing since the beginning of the 21st century, demonstrating its importance as a global strategic resource and a livelihood material. According to existing reports [1], 61% of natural gas reserves which were discovered worldwide in the past decade are distributed in deep zones and are characterized by strong ground stress, high heterogeneity, temperature, and pressure. High pressure gas reservoirs have downward curving p/Z~Gp behavior. Precise estimation of the reserves of such gas reservoirs is crucial for reservoir engineers. However, reserve estimation’s basic method is the material balance method [2]. The formation compressibility is one of the key parameters in the material balance equation of a high pressure gas reservoir but challenging to determine; the error in estimates could even exceed 100% [3]. For many years, the petroleum industry has relied on Hall’s correlation [4] for estimating pore volume compressibility. This correlation was developed from measurements on seven consolidated limestone and five consolidated sandstone samples. Compressibilities, however, are highly affected by reservoir type and overburden conditions.According to the classical theory, the material balancep/Z~Gp curve of a high pressure gas reservoir has a two-segment characteristic [5], with the second straight line segment occurring late (apparent formation pressure drops by 6%-38%). However, for high pressure and ultrahigh pressure gas reservoirs at the early stage of development, even if the production test time is up to 1 year and the drawdown scope reaches 3%-38% of the initial formation pressure or even higher, the starting point deviating from the early straight line segment would still not occur, and thus, the starting condition of using the material balance method to calculate reserves cannot be met. If conditions are not fully met, arbitrary applications of this method will result in serious deviation in the calculation [6]. In recent years, tremendous efforts have been spent in designing methods independent of formation compressibility, including the type curve match analysis and nonlinear regression methods. Table 1 indicates the scope of various types of information on this subject, listing the relevant papers that have been published in the petroleum literature.Table 1
Published studies of reserve estimate problem of high pressure gas reservoir.
No.Ref.YearAuthorSolution methodNeed pore volume compressibilityAdvantages and disadvantages1[3]1971HammerlindlPolyline analysis methodHammerlindl methodAverage formation compressibility✓Gradient>1.13MPa/100m2Material balance✓3[7]1983ChenChen Yuanqian method✓Volumetric gas reservoirs, pressure coefficient at the transition point is 1.2~1.34[8]2001GanGan-BlasingameХNo pore volume compressibility is required, and it can be calculated whether the curve inflection point appears or not5[9]1981RamagostLinear regression methodRamagost-Farshad✓Volumetric gas reservoirs6[10]1981RoachRoachХPore volume compressibility can be calculated, but it is sensitive to the original pressure data7[11]1994PostonImproved Roach methodХIt is suitable for water drive gas reservoirs and can calculate the reserves before water invasion, the size of water invasion, and the effective compressibility coefficient8[12]1993Becerra-ArteagaBecerra-ArteagaХNeed gradient9[2]1963HavlenaHavlena-OdehХSuitable for water intrusion gas reservoirs, sensitive to original pressure and early data; if it is a closed gas reservoir, the rock compressibility coefficient can be calculated10[13]1993YuanqianNonlinear regression methodBinary regressionХKnowing the variation law of natural gas deviation coefficient; this method can calculate dynamic reserves and compressibility11[14]2008GonzalezNonlinear regression methodQuadratic production modelХDimensionless linear coefficient is less than 0.412[15]2011QinTrinomial approximationХThe shape of the curve is not easy to control and may cause large errors13[16]2017JiaoLimit formХVolumetric gas reservoirs14[17]2019SunPower function formХVolumetric gas reservoirs15[18]1991AmbasthaType curve matchingAmbastha methodХVolumetric gas reservoirs, original pressure, and temperature are known in type curve matching, strong multisolution16[19]1998FetkovichFetkovich methodХConsidering the influence of nonreservoir dissolved gas, inversely calculate the compressibility and the size of the water body17[14]2008GonzalezGonzalez methodQuadratic production modelХMainly used to verify the results of other methods18[6]2020SunSingle log match analysisХVolumetric gas reservoirs19[20]2001MarhaendrajanaMultiwell modern production decline analysis method✓Without considering the impact of water invasion, only reserves of connected well groups can be calculated20[21]1998WalshTrial analysis method✓It can be used to calculate the volume of water bodies and calculate the water-soluble gas reservesThis has provided new solutions to the reserve estimate problem of a high pressure gas reservoir when the starting condition is met in the middle and late stages of development. However, a simple but practical method is still missing for the reserve estimate in the early development stage.To bridge the gap, starting with the definition of formation compressibility and the material balance equation of volumetric gas reservoir, we established an analysis method for the reserve estimate based on gas production of cumulative unit pressure drop and evaluated its performance in three gas reservoirs.


## 2. Gas Production of Cumulative Unit Pressure Drop

### 2.1. Isothermal Compressibility of Gas

In general, the isothermal coefficient of compressibility of a material is defined as(1)Cf=−1V∂V∂pT,where Cf denotes isothermal compressibility of gas in 1/MPa, with p denoting pressure in MPa and V denoting volume in 108 m3.


### 2.2. Gas Production of Cumulative Pressure Drop

Gas production per unit pressure drop is defined as cumulative gas production per unit drop of average formation pressure. Gas production of cumulative unit pressure drop is defined as the ratio of cumulative gas production to cumulative drawdown of formation pressure from the initial formation pressure condition to the current formation pressure condition.When equation (1) is applied to a gas reservoir,
(2)Cf=−1GdGdpT.When separation variables are integrated,(3)Gp=G∫ppiCfdp=GC_fpi−p,where GC_f is the gas production of cumulative unit pressure drop, 108 m3/MPa. When the formation pressure drops to 0.101325 MPa, there is
(4)Gp=G=GC_fpi,i.e., when the formation pressure drops to normal atmospheric pressure, the dimensionless gas production of the cumulative unit pressure drop is the reciprocal of the initial formation pressure.


### 2.3. Variation Law of Gas Production of Cumulative Unit Pressure Drop

The material balance equation of a volumetric gas reservoir can be expressed as(5)pZ=pZi1−GpG,where p is the pressure, MPa; Gp is the cumulative gas production, 108 m3; G is the reserves, 108 m3; and the subscript i represents the initial state. Equation (5) is deformed, and the dimensionless cumulative gas production is
(6)GpG=1−p/Zp/Zi.When the pressure drops topj, the dimensionless cumulative gas production is
(7)GpGj=1−p/Zjp/Zi.When the pressure drops topj+1, the dimensionless gas production is
(8)GpGj+1=1−p/Zj+1p/Zi.The dimensionless gas production per unit pressure drop (assumingpj−pj+1=1) is
(9)GpGj+1−GpGj=p/Zjp/Zi−p/Zj+1p/Zi.The dimensionless gas production of cumulative unit pressure drop is(10)Gp/Gjpi−pj=1−p/Zj/p/Zipi−pj.The theoretical calculation results show that with the depletion of formation pressure, affected by the deviation factor, the gas production per unit pressure drop first increases and then decreases, while the dimensionless gas production of cumulative unit pressure drop increases gradually. When the formation pressure drops to atmospheric pressure, its value is the reciprocal of the initial formation pressure, as shown in Figure1.Figure 1
Gas production per unit pressure drop and dimensionless gas production of cumulative unit pressure drop (formation pressure is 40 MPa, formation temperature is 373.15 K, and gravity of gas is 0.6).Figure2 shows the variation of gas production of cumulative unit pressure drop with the drop of formation pressure at different temperatures, when the formation pressure is 40 MPa and the gravity of gas is 0.6. The relation curve of gas production of cumulative unit pressure drop and cumulative pressure drop is linear at early and middle stages; when the linear relation is prolonged to intersect with the vertical axis, the ratio a of gas production of cumulative unit pressure drop to intersection value is approximately equal to the minimum deviation factor. In Figure 2(a), the temperature is 333.15 K, a =0.025/0.0284 = 0.88 ≈ 0.85, and the relative error is 3.6%; in Figure 2(b), the temperature is 373.15 K, a =0.025/0.0264 = 0.95 ≈ 0.94, and the relative error is 1.1%; in Figure 2(c), the temperature is 413.15 K, a =0.025/0.0267 = 0.94 ≈ 0.96, and the relative error is 2.1%.Figure 2
Effect of temperature on gas production of cumulative pressure drop.
(a)
333.15 K(b)
373.15 K(c)
413.15 K


### 2.4. Minimum Deviation Factor of Natural Gas

The minimum deviation factor of natural gas depends on the pseudocritical property [22], as shown in Figure 3. When the pseudoreduced temperature is more than 1.9, the minimum deviation factor of natural gas is above 0.90.Figure 3
Minimum deviation factor at different pseudoreduced pressures.When the pseudoreduced temperature is unknown, it can be calculated based on the gravity of natural gas, as given by(11)Tpc=168+325γg−12.5γg21.8.


### 2.5. Calculation Method of Reserves

To sum up, the Cartesian plot of gas production of the cumulative unit pressure drop and the cumulative pressure drop is drawn, the linear regression is carried out, the straight line is prolonged to intersect with the vertical axis (as shown in Figure2), and the product of the intersection ordinate value and the minimum deviation factor gives the reserve.


## 2.1. Isothermal Compressibility of Gas

In general, the isothermal coefficient of compressibility of a material is defined as(1)Cf=−1V∂V∂pT,where Cf denotes isothermal compressibility of gas in 1/MPa, with p denoting pressure in MPa and V denoting volume in 108 m3.


## 2.2. Gas Production of Cumulative Pressure Drop

Gas production per unit pressure drop is defined as cumulative gas production per unit drop of average formation pressure. Gas production of cumulative unit pressure drop is defined as the ratio of cumulative gas production to cumulative drawdown of formation pressure from the initial formation pressure condition to the current formation pressure condition.When equation (1) is applied to a gas reservoir,
(2)Cf=−1GdGdpT.When separation variables are integrated,(3)Gp=G∫ppiCfdp=GC_fpi−p,where GC_f is the gas production of cumulative unit pressure drop, 108 m3/MPa. When the formation pressure drops to 0.101325 MPa, there is
(4)Gp=G=GC_fpi,i.e., when the formation pressure drops to normal atmospheric pressure, the dimensionless gas production of the cumulative unit pressure drop is the reciprocal of the initial formation pressure.


## 2.3. Variation Law of Gas Production of Cumulative Unit Pressure Drop

The material balance equation of a volumetric gas reservoir can be expressed as(5)pZ=pZi1−GpG,where p is the pressure, MPa; Gp is the cumulative gas production, 108 m3; G is the reserves, 108 m3; and the subscript i represents the initial state. Equation (5) is deformed, and the dimensionless cumulative gas production is
(6)GpG=1−p/Zp/Zi.When the pressure drops topj, the dimensionless cumulative gas production is
(7)GpGj=1−p/Zjp/Zi.When the pressure drops topj+1, the dimensionless gas production is
(8)GpGj+1=1−p/Zj+1p/Zi.The dimensionless gas production per unit pressure drop (assumingpj−pj+1=1) is
(9)GpGj+1−GpGj=p/Zjp/Zi−p/Zj+1p/Zi.The dimensionless gas production of cumulative unit pressure drop is(10)Gp/Gjpi−pj=1−p/Zj/p/Zipi−pj.The theoretical calculation results show that with the depletion of formation pressure, affected by the deviation factor, the gas production per unit pressure drop first increases and then decreases, while the dimensionless gas production of cumulative unit pressure drop increases gradually. When the formation pressure drops to atmospheric pressure, its value is the reciprocal of the initial formation pressure, as shown in Figure1.Figure 1
Gas production per unit pressure drop and dimensionless gas production of cumulative unit pressure drop (formation pressure is 40 MPa, formation temperature is 373.15 K, and gravity of gas is 0.6).Figure2 shows the variation of gas production of cumulative unit pressure drop with the drop of formation pressure at different temperatures, when the formation pressure is 40 MPa and the gravity of gas is 0.6. The relation curve of gas production of cumulative unit pressure drop and cumulative pressure drop is linear at early and middle stages; when the linear relation is prolonged to intersect with the vertical axis, the ratio a of gas production of cumulative unit pressure drop to intersection value is approximately equal to the minimum deviation factor. In Figure 2(a), the temperature is 333.15 K, a =0.025/0.0284 = 0.88 ≈ 0.85, and the relative error is 3.6%; in Figure 2(b), the temperature is 373.15 K, a =0.025/0.0264 = 0.95 ≈ 0.94, and the relative error is 1.1%; in Figure 2(c), the temperature is 413.15 K, a =0.025/0.0267 = 0.94 ≈ 0.96, and the relative error is 2.1%.Figure 2
Effect of temperature on gas production of cumulative pressure drop.
(a)
333.15 K(b)
373.15 K(c)
413.15 K


## 2.4. Minimum Deviation Factor of Natural Gas

The minimum deviation factor of natural gas depends on the pseudocritical property [22], as shown in Figure 3. When the pseudoreduced temperature is more than 1.9, the minimum deviation factor of natural gas is above 0.90.Figure 3
Minimum deviation factor at different pseudoreduced pressures.When the pseudoreduced temperature is unknown, it can be calculated based on the gravity of natural gas, as given by(11)Tpc=168+325γg−12.5γg21.8.


## 2.5. Calculation Method of Reserves

To sum up, the Cartesian plot of gas production of the cumulative unit pressure drop and the cumulative pressure drop is drawn, the linear regression is carried out, the straight line is prolonged to intersect with the vertical axis (as shown in Figure2), and the product of the intersection ordinate value and the minimum deviation factor gives the reserve.


## 3. Case Analysis

### 3.1. Small Ultrahigh Pressure Gas Reservoir: Anderson “L” Gas Reservoir

The basic parameters of the American Anderson “L” gas reservoir [23] are as follows: the buried depth is 3404.5 m, the pressure coefficient is 1.907 MPa/100 m, the initial pressure is 65.55 MPa, the temperature is 403 K, the gravity of natural gas is 0.656, and the volumetric reserves are 19.68×108m3, with the producing history listed as in Table 2.Table 2
Producing history of Anderson “L” gas reservoir (Duggan, 1972).
pGpp/ZΔpGp/ΔpMPa108 m3MPaMPa108 m3/MPa65.5480.00045.5200.000064.0660.11845.181.4820.079661.8460.49244.593.7020.132959.2600.96644.096.2880.153657.4471.27643.658.1010.157555.2201.64743.0710.3280.159552.4212.25742.3113.1270.171951.0622.62041.9214.4860.180948.2773.14641.0517.2710.182246.3403.51940.4019.2080.183245.0573.82739.9820.4910.186839.7415.16337.9225.8070.200132.8606.83633.6332.6880.209129.6138.38931.9135.9350.233425.8559.68929.0239.6930.244122.38710.93126.2143.1610.2533The relation curve of the cumulative pressure drop and the gas production of the cumulative unit pressure drop is given in Figure4; linear regression is carried out for the data points, the straight line is prolonged to intersect with the vertical axis, and the intersection ordinate value is 0.31×108m3/MPa, which is multiplied by the initial formation pressure 65.548 MPa, and the reserve is obtained as 20.29×108m3. The minimum deviation factor is 0.935 based on Figure 3. Therefore, the reserves of the gas reservoir is estimated as 20.29×0.935, i.e., 19.00×108m3. The reserves calculated by this method are 3.5% smaller than that by the volumetric method. The p/Z curve method, which gave an estimation of 26.6×108m3, is shown in Figure 5.Figure 4
Reserves of Anderson “L” gas reservoir calculated by the method of gas production of cumulative unit pressure drop.Figure 5
Reserves of Anderson “L” gas reservoir calculated by thep/Z curve method.


### 3.2. Medium Ultrahigh Pressure Gas Reservoir: Louisiana Offshore Gas Reservoir

The basic parameters of the American Louisiana offshore ultrahigh pressure gas reservoir [23] are as follows: the buried depth is 4055 m, the pressure coefficient is 1.95 MPa/100 m, the initial pressure is 78.903 MPa, the temperature is 401.55 K, the gravity of the natural gas is 0.60, the pseudoreduced temperature is 1.96, the initial water saturation is 0.22, and the reserves are 131.6×108m3, with the producing history listed as in Table 3.Table 3
Producing history of Louisiana offshore gas reservoir (Ramagost, 1981).
pGpp/ZΔpGp/ΔpMPa108 m3MPaMPa108 m3/MPa78.9030.00052.7430.0000.000073.5952.80951.1785.3090.529169.8518.10450.0009.0530.895263.79715.17847.96815.1061.004759.11621.99446.18419.7881.111554.51028.71944.31724.3941.177350.88334.08242.68728.0201.216347.20841.06240.90831.6951.295544.04445.48539.25534.8601.304840.17651.63337.06238.7281.333237.29455.99135.28341.6101.345634.47461.06833.37244.4301.374531.02666.75430.87247.8771.394328.75169.63129.10050.1521.3884The relation curve of the cumulative pressure drop and gas production of the cumulative unit pressure drop is given in Figure6. The intersection ordinate value is 1.788×108m3/MPa, which is multiplied by the initial formation pressure 78.903 MPa leading to the reserves of 141.07×108m3. The pseudoreduced temperature of the gas reservoir is 1.96, and the minimum deviation factor is 0.93. Therefore, the reserve of the gas reservoir is 141.07×0.93, i.e., 131.2×108m3. The reserves calculated by this method are highly consistent with that calculated by the dynamic method in the references [9]. If the reserves are calculated directly by the p/Z curve method, the calculated results are 161.2×108m3, as shown in Figure 7, which are much larger than the values calculated by the new method.Figure 6
Reserves of Louisiana offshore gas reservoir calculated by the method of gas production of cumulative unit pressure drop.Figure 7
Reserves of Louisiana offshore gas reservoir calculated by thep/Z curve method.


### 3.3. Large High Pressure Gas Reservoir: X Gas Reservoir

The basic parameters of the X gas reservoir [16] are as follows: the middle buried depth is 3750 m, the pressure coefficient is 2.0 MPa/100 m, the initial formation pressure is 74.48 MPa, the minimum deviation factor is 0.93, and the reserves calculated by the volumetric method are 2091.5×108m3, with the producing history listed as in Table 4.Table 4
Producing history of X gas reservoir.
pGpp/ZΔpGp/ΔpMPa108 m3MPaMPa108 m3/MPa74.480.0053.060074.332.5853.010.1517.2072.6734.9652.421.8119.3168.55119.7450.955.9320.1963.53226.4849.1310.9520.6858.45338.6047.1816.0321.1253.84451.0445.2520.6421.8549.53555.2543.2624.9522.2545.70659.4641.2928.7822.9141.71763.9639.0132.7723.3138.06868.1736.6936.4223.8434.59972.3834.2439.8924.3831.011076.5931.4743.4724.7727.831181.0828.7946.6525.3224.731285.2925.9849.7525.8321.521389.5022.8952.9626.2418.461493.1619.7856.0226.6515.551593.8916.7258.9327.05The relation curve of the cumulative pressure drop and gas production of the cumulative unit pressure drop is given in Figure8. The intersection ordinate value is 29.344×108m3/MPa, which is multiplied by the initial formation pressure 74.48 MPa leading to the reserves of 2185.54×108m3. The minimum deviation factor is 0.93. Therefore, the reserve of the gas reservoir is 2185.54×0.93, i.e., 2032.55×108m3, which is 2.8% lower than the reserves calculated by the volumetric method. If the reserves are calculated directly by the p/Z method, the calculated result is 2460.4×108m3, as shown in Figure 9, which is 17.6% higher than the reserves calculated by the static method.Figure 8
Reserves of X gas reservoir calculated by the method of gas production of cumulative unit pressure drop.Figure 9
Reserves of X gas reservoir calculated by thep/Z method.


### 3.4. Comparison with Calculated Results of Other Methods without Considering Formation Compressibility

As seen in Table5, the calculated results of the method introduced in the paper are close to those of the nonlinear regression method and the semilogarithmic type curve match method, and the three cases further verify the correctness of the method. When the gas reservoir is in the early stage of development and there is no obvious inflection point in the p/Z curve, this method can be used to approximately estimate the reserves of the gas reservoir.Table 5
Comparison of reserves calculated by different methods.
Gas reservoirVolumetric method reserves (108 m3)Reserves (108m3)p/Z curve methodNonlinear regression methodSemilog type curve match methodMethod introduced in the paperAnderson “L”19.6826.619.920.019.0Louisiana offshore131.6161.2125.3130.0131.2X2091.52460.42032.02000.02032.55


## 3.1. Small Ultrahigh Pressure Gas Reservoir: Anderson “L” Gas Reservoir

The basic parameters of the American Anderson “L” gas reservoir [23] are as follows: the buried depth is 3404.5 m, the pressure coefficient is 1.907 MPa/100 m, the initial pressure is 65.55 MPa, the temperature is 403 K, the gravity of natural gas is 0.656, and the volumetric reserves are 19.68×108m3, with the producing history listed as in Table 2.Table 2
Producing history of Anderson “L” gas reservoir (Duggan, 1972).
pGpp/ZΔpGp/ΔpMPa108 m3MPaMPa108 m3/MPa65.5480.00045.5200.000064.0660.11845.181.4820.079661.8460.49244.593.7020.132959.2600.96644.096.2880.153657.4471.27643.658.1010.157555.2201.64743.0710.3280.159552.4212.25742.3113.1270.171951.0622.62041.9214.4860.180948.2773.14641.0517.2710.182246.3403.51940.4019.2080.183245.0573.82739.9820.4910.186839.7415.16337.9225.8070.200132.8606.83633.6332.6880.209129.6138.38931.9135.9350.233425.8559.68929.0239.6930.244122.38710.93126.2143.1610.2533The relation curve of the cumulative pressure drop and the gas production of the cumulative unit pressure drop is given in Figure4; linear regression is carried out for the data points, the straight line is prolonged to intersect with the vertical axis, and the intersection ordinate value is 0.31×108m3/MPa, which is multiplied by the initial formation pressure 65.548 MPa, and the reserve is obtained as 20.29×108m3. The minimum deviation factor is 0.935 based on Figure 3. Therefore, the reserves of the gas reservoir is estimated as 20.29×0.935, i.e., 19.00×108m3. The reserves calculated by this method are 3.5% smaller than that by the volumetric method. The p/Z curve method, which gave an estimation of 26.6×108m3, is shown in Figure 5.Figure 4
Reserves of Anderson “L” gas reservoir calculated by the method of gas production of cumulative unit pressure drop.Figure 5
Reserves of Anderson “L” gas reservoir calculated by thep/Z curve method.


## 3.2. Medium Ultrahigh Pressure Gas Reservoir: Louisiana Offshore Gas Reservoir

The basic parameters of the American Louisiana offshore ultrahigh pressure gas reservoir [23] are as follows: the buried depth is 4055 m, the pressure coefficient is 1.95 MPa/100 m, the initial pressure is 78.903 MPa, the temperature is 401.55 K, the gravity of the natural gas is 0.60, the pseudoreduced temperature is 1.96, the initial water saturation is 0.22, and the reserves are 131.6×108m3, with the producing history listed as in Table 3.Table 3
Producing history of Louisiana offshore gas reservoir (Ramagost, 1981).
pGpp/ZΔpGp/ΔpMPa108 m3MPaMPa108 m3/MPa78.9030.00052.7430.0000.000073.5952.80951.1785.3090.529169.8518.10450.0009.0530.895263.79715.17847.96815.1061.004759.11621.99446.18419.7881.111554.51028.71944.31724.3941.177350.88334.08242.68728.0201.216347.20841.06240.90831.6951.295544.04445.48539.25534.8601.304840.17651.63337.06238.7281.333237.29455.99135.28341.6101.345634.47461.06833.37244.4301.374531.02666.75430.87247.8771.394328.75169.63129.10050.1521.3884The relation curve of the cumulative pressure drop and gas production of the cumulative unit pressure drop is given in Figure6. The intersection ordinate value is 1.788×108m3/MPa, which is multiplied by the initial formation pressure 78.903 MPa leading to the reserves of 141.07×108m3. The pseudoreduced temperature of the gas reservoir is 1.96, and the minimum deviation factor is 0.93. Therefore, the reserve of the gas reservoir is 141.07×0.93, i.e., 131.2×108m3. The reserves calculated by this method are highly consistent with that calculated by the dynamic method in the references [9]. If the reserves are calculated directly by the p/Z curve method, the calculated results are 161.2×108m3, as shown in Figure 7, which are much larger than the values calculated by the new method.Figure 6
Reserves of Louisiana offshore gas reservoir calculated by the method of gas production of cumulative unit pressure drop.Figure 7
Reserves of Louisiana offshore gas reservoir calculated by thep/Z curve method.


## 3.3. Large High Pressure Gas Reservoir: X Gas Reservoir

The basic parameters of the X gas reservoir [16] are as follows: the middle buried depth is 3750 m, the pressure coefficient is 2.0 MPa/100 m, the initial formation pressure is 74.48 MPa, the minimum deviation factor is 0.93, and the reserves calculated by the volumetric method are 2091.5×108m3, with the producing history listed as in Table 4.Table 4
Producing history of X gas reservoir.
pGpp/ZΔpGp/ΔpMPa108 m3MPaMPa108 m3/MPa74.480.0053.060074.332.5853.010.1517.2072.6734.9652.421.8119.3168.55119.7450.955.9320.1963.53226.4849.1310.9520.6858.45338.6047.1816.0321.1253.84451.0445.2520.6421.8549.53555.2543.2624.9522.2545.70659.4641.2928.7822.9141.71763.9639.0132.7723.3138.06868.1736.6936.4223.8434.59972.3834.2439.8924.3831.011076.5931.4743.4724.7727.831181.0828.7946.6525.3224.731285.2925.9849.7525.8321.521389.5022.8952.9626.2418.461493.1619.7856.0226.6515.551593.8916.7258.9327.05The relation curve of the cumulative pressure drop and gas production of the cumulative unit pressure drop is given in Figure8. The intersection ordinate value is 29.344×108m3/MPa, which is multiplied by the initial formation pressure 74.48 MPa leading to the reserves of 2185.54×108m3. The minimum deviation factor is 0.93. Therefore, the reserve of the gas reservoir is 2185.54×0.93, i.e., 2032.55×108m3, which is 2.8% lower than the reserves calculated by the volumetric method. If the reserves are calculated directly by the p/Z method, the calculated result is 2460.4×108m3, as shown in Figure 9, which is 17.6% higher than the reserves calculated by the static method.Figure 8
Reserves of X gas reservoir calculated by the method of gas production of cumulative unit pressure drop.Figure 9
Reserves of X gas reservoir calculated by thep/Z method.


## 3.4. Comparison with Calculated Results of Other Methods without Considering Formation Compressibility

As seen in Table5, the calculated results of the method introduced in the paper are close to those of the nonlinear regression method and the semilogarithmic type curve match method, and the three cases further verify the correctness of the method. When the gas reservoir is in the early stage of development and there is no obvious inflection point in the p/Z curve, this method can be used to approximately estimate the reserves of the gas reservoir.Table 5
Comparison of reserves calculated by different methods.
Gas reservoirVolumetric method reserves (108 m3)Reserves (108m3)p/Z curve methodNonlinear regression methodSemilog type curve match methodMethod introduced in the paperAnderson “L”19.6826.619.920.019.0Louisiana offshore131.6161.2125.3130.0131.2X2091.52460.42032.02000.02032.55


## 4. Conclusions and Suggestions

(1)
Since it is difficult to accurately determine the formation compressibility in the material balance equation of high pressure gas reservoirs, the method without considering the formation compressibility is recommended for reserve estimation(2)
When the formation pressure drops to the normal atmospheric pressure, the dimensionless gas production of the cumulative unit pressure drop is the reciprocal of the initial formation pressure. The relation curve of the dimensionless gas production of the cumulative unit pressure drop and pressure drop is a straight line in the medium term. The reserve correction coefficient is related to temperature(3)
The method of gas production of cumulative unit pressure drop, the semilog curve match analysis technique, and the nonlinear regression method can be used in turn to calculate the reserves based on the length of production time. If the production time is short, the reserves can be roughly estimated by the method of gas production of the cumulative unit pressure drop. Based on the reserve estimate results, combined with material balance equation, the effective compressibility coefficient can be calculated, and then, the water influx can be calculated(4)
For the water drive gas reservoir, the reserves calculated by the method of gas production of the cumulative unit pressure drop are relatively optimistic; in the future, the evaluation method of water drive gas reservoirs should be further studied to reduce the ambiguity of analysis results



---

*Source: 9926983-2021-06-01.xml*