Integration of Multi-Region Material Balance Equation with Binomial Productivity Equation for Performance Prediction

Liu, Yan; Zhu, Langtao

doi:https://doi.org/10.1155/2022/5169602

Geofluids

On this page

Abstract Introduction Conclusions Appendix Data Availability Conflicts of Interest References Copyright Related Articles

Special Issue

Flow and Transport Dynamics in Fractured Porous Media 2021

View this Special Issue

Research Article | Open Access

Volume 2022 | Article ID 5169602 | https://doi.org/10.1155/2022/5169602

Integration of Multi-Region Material Balance Equation with Binomial Productivity Equation for Performance Prediction

Yan Liu¹and Langtao Zhu¹

Academic Editor: Mohammed Fattah

Received14 Jan 2022

Revised28 Feb 2022

Accepted10 Mar 2022

Published30 Mar 2022

Abstract

In order to overcome the limitations of conventional material balance methods, a new multi-region material balance method, considering the strong heterogeneity and low-velocity non-Darcy flow characteristics of bio-reef gas reservoir, is proposed to estimate the development performance of reef reservoirs. Firstly, the reef reservoir is divided into two regions, the reef cover high-permeability development area and the reef edge low-permeability supplying area. Then, based on one-dimensional non-linear seepage mechanism and mass conservation principle, the multi-region supplying material balance equation is established. Finally, combined with the binomial productivity equation, the development performance can be easily estimated at different curves. The application results indicate that the low-permeability replenishment area plays an important role, especially after entering the decline stage, on the development performance of high-permeability region. Compared with the traditional material balance method, the new method can help for better understanding of unbalanced development in different regions and the replenishment effect of low-permeability region.

1. Introduction

A typical gas well needs to go through three periods, including constant gas rate production, production decline, and pressurization production [1]. At the stage of constant rate production, traditional material balance, production dynamic fitting, and numerical simulation are three main methods, which are usually adopted to predict the dynamic of gas reservoirs [2]. However, these technologies are not applicable to bio-reef gas reservoirs for the heterogeneity and the low-velocity non-Darcy flow characteristics. (1)Traditional material balance is a simple and convenient method to evaluate the dynamic reserves and only can carry out simple development performance prediction [3, 4]. However, the application indicates that the traditional material balance cannot accomplish to accurately describe the production characteristics of different types of reservoirs(2)Numerical simulations are common methods used to analyze the dynamic characteristics [5, 6]. However, long well spacing and sparse well pattern make it difficult to build accurate geological models, and it is inconvenient to be applied to evaluate the production characteristics of low-permeability reservoirs(3)Production performance fitting methods are another method to reveal the low-permeability reservoir seepage law [7, 8]. It requires more detailed geological information and reservoir fluid parameters, which cannot be adopted to quickly evaluate the production characteristics of heterogeneous reservoirs

Recently, the multi-region material balance equation (MMBE) is proposed. The evaluation result of this method matches well with the numerical simulation, which verifies its reliability. Chengtai Gao et al. [9–11] carried out indoor seepage experiment to analyze the performance of multi-region reservoirs, and initially realized dynamic prediction based on constant rate production period data, but lacks a complete theoretical basis. Based on the concept of multi-region replenishment, Shiqing Cheng et al. [12–16] conducted an in-depth dynamic evaluation of the low-permeability gas reservoirs, but early predicting the dynamic development characteristics of a single well cannot realize. In addition, the MMBE, considering the low-velocity non-Darcy flow of low-permeability reservoirs, was used to evaluate the elastic flooding and dissolved gas flooding recovery efficiency [17–20]. However, the reservoir is simplified to a number of annular bands which made the theoretical foundation imperfect.

In this work, a new MMBE method is developed to carry out gas well performance prediction. Compared with the previous work [9–22], this new method initially introduced the effects of low-velocity non-Darcy seepage and combined the MMBE with binomial productivity equation for development performance prediction. The specific steps of this method are as follows: firstly, based on the architectural model of reef body and low-velocity non-Darcy seepage law in low-permeability region, the heterogeneous bio-reef reservoirs could be simplified to 1D reservoir model with two regions. Then, considering the low-velocity non-Darcy seepage in low-permeability region, a new mathematic method derived from heterogeneous multi-region model is proposed to estimate the reservoir parameter. Subsequently, integrating the multi-region material balance equation with binomial productivity equation, an integrated approach assessing the gas well development performance is derived. Finally, practical application was conducted, and some important conclusions are obtained. The flow chart of using this method to carry out development performance prediction can be seen in Figure 1.

2. Dynamic Prediction Method with Multi-Region Replenishment

2.1. Reservoir Characteristics of Typical Bio-Reef Gas Reservoirs

The bio-reef gas reservoir presents the characteristics of multi-stage sedimentation; the reef core and reef front are impervious regions. The reef cover with high-permeability is a development area, and the reef edge with low-permeability is a replenishment area, shown in Figure 2. The low-permeability region exhibits strong non-Darcy flow characteristics.

As shown in Figure 3, the fluid seepage of low-permeability reservoir is divided into two stages: (1)Affected by gas-solid interaction, the first stage shows low-velocity non-Darcy seepage characteristics, and the gas seepage velocity closes to 0. For the convenience of research, we assume that the seepage velocity is equal to 0, when the pressure gradient is lower than , non-Darcy seepage factor, shown in Figure 3 and Eq. (1)(2)When the pressure gradient is higher than , the second stage behaves linear Darcy flow. As shown in Figure 3 and Eq. (1), there is a positive correlation between seepage velocity and pressure gradient

2.2. Multi-Region Material Balance Equation

2.2.1. Physical Model

The reef configuration shows that the gas reservoir is divided into two parts: reef cover high-permeability development area and reef edge low-permeability supply area, as shown in Figure 4. For the convenience of research, a typical one-dimensional two-zone seepage model is used to approximately describe the configuration characteristics of the typical reef. The assumptions of the model are as follows: (1)The reservoirs in high-permeability area and low-permeability area are horizontally homogeneous with equal thickness, ignoring the influence of gravity and water(2)The reservoir in the low-permeability area is dense and has strong gas-solid interaction, showing the characteristics of low-velocity non-Darcy seepage, while the reservoir in the reef cover has high-permeability and the seepage meets the characteristics of Darcy seepage(3)The supply of the two sections is controlled by the pressure difference between the two regions, the supply gas volume is directly proportional to the pressure difference between the two areas, and the supply flow distance is equal to , shown in Figure 4(4)The length of high-permeability area and low-permeability area is and , respectively, the width is , and the formation thickness is (5)In the high-permeability area and the low-permeability area, the geological reserves are and , respectively, and the permeability is and , respectively(6)The gas production rate in the high-permeability area is

2.2.2. Mathematical Model and Solution

Reservoir heterogeneity has a greater impact on the development of bio-reef gas reservoirs and fluid migration between different regions. The low-permeability area shows obvious low-velocity non-Darcy seepage characteristics, which makes the pressure drop in the reef edge low-permeability area slower than that in the reef cover high-permeability area, and the formation pressure difference between the two areas is obvious.

Affected by the physical property difference between the two areas, the gas well production, in the high-permeability area, shows two obvious development stages. In the initial stage, only the fluid in the reef cover high-permeability development area flows to the near well region. When the pressure gradient between the two areas is greater than the low-velocity non-Darcy flow factor G, the fluid in the reef edge low-permeability area starts to replenish to the reef cover high-permeability area. (1)Early seepage stage in high-permeability region

Affected by the low-velocity non-Darcy seepage characteristics in the low-permeability area, and when the pressure gradient between the two areas is less than the low-velocity non-Darcy flow factor G. Therefore, only the fluid in the reef covers high-permeability area seepages, and the fluid in the low-permeability area does not flow.

Based on the principle of mass balance, the initial material balance equation in the early seepage stage can be obtained. For the specific formula derivation process, see Appendix A. (2)Late replenishment stage in low-permeability region

When the cumulative gas production volume is greater than , see Figure 3, the pressure gradient between the two areas is greater than the low-velocity non-Darcy flow factor G, and the fluid in the low-permeability area begins to supply to the high-permeability area.

Based on the principle of zonal replenishment, the material balance equations of different regions can be obtained. For the specific formula derivation process, see Appendix B.

2.3. Dynamic Prediction Method

Assuming that the development is divided into two development stages, constant rate production and decreasing rate production, and there is only one well in the high-permeability area. After productivity testing, we can get the binomial productivity equation of the gas well: where is bottom hole flowing pressure, MPa; and are the binomial capacity equation coefficients.

By synthesizing binomial production equation and multi-region material balance equation, it can be able to predict the development performance of gas wells. (1)Calculation of constant rate production period (1)Calculation of high-permeability development zone

Suppose a single well is initially produced at a constant rate until the bottom hole pressure drops to the lowest value . According to Eq. (5), the reservoir pressure in the high-permeability development zone at the end of the constant rate production period is determined by the following equation:

The reservoir pressure interval [, ] is divided into uniform intervals with node ( =0, 1, 2, ..., ). The production rate and pressure , in the constant rate production period, are brought into Eq. (2) or Eq. (3) to obtain the cumulative production time , the cumulative gas production volume , and the apparent formation pressure () of high-permeability area.

And the remaining reserves in the high-permeability development zone can be obtained through the following formula: (2)Calculation of low-permeability recharge area

Furthermore, the cumulative gas production , the cumulative production time , and the gas production rate , in the constant rate production period, are brought into Eq. (4). And we can get the apparent formation pressure (/) of low-permeability area.

According to the following Eq. (8), the cumulative gas volume supply from the low-permeability area to the high-permeability development area is:

In addition, the remaining reserves in the low-permeability development zone can be obtained by the following formula: (2)Calculation of production rate decline period (1)Calculation of high-permeability development zone

Assuming that the economic limit production is , and the abandoned reservoir pressure can be obtained by the following formula:

Furthermore, the reservoir pressure and abandoned formation pressure interval [, ] are divided into any equal cells, and the node pressure is (). Similar to the calculation of constant rate production period, the average reservoir pressure and average production in interval can be obtained, and the relationship between those parameters is determined by Eq. (5), that is.

Solving Eq. (11) can be obtained:

Incorporating the average production rate and node reservoir pressure into Eq. (2), we can get the cumulative production time , the cumulative gas production volume , and the apparent formation pressure () of high-permeability area. And the remaining reserves in the high-permeability development zone can be obtained through the following formula: (2)Calculation of low-permeability recharge area

Furthermore, cumulative gas production volume , cumulative production time , and average production rate are taken into Eq. (3), and the apparent formation pressure () of the low-permeability replenishment area can be obtained. The cumulative gas replenishment volume from the low-permeability replenishment area to the high-permeability development area can be obtained by the following equation:

In addition, the remaining reserves in the low-permeability development zone can be obtained by the following formula:

2.4. The Limitations of This Method

Take into account that the physical and mathematical models have many presuppositions, the proposed method has the following three limitations: (1)Reservoirs in the high- and low-permeability zones are simplified into two independent homogeneous whole, which cannot reflect the internal seepage characteristics of the two areas. And the evaluation parameters can only reflect the overall characteristics of each seepage unit(2)This method can only evaluate the reserves in the low-permeability area and the average distance from the low-permeability area to the high-permeability area. And the distribution characteristics of low-permeability reserves can only be determined by combining geological analysis(3)As the binomial productivity equation will change with the formation pressure drop, introducing the initial binomial productivity equation to predict the development performance has relatively low accuracy

3. Curve Characteristics and Practical Application

3.1. Curve Characteristics

Based on the above proposed new multi-region material balance, the relationship curve between the apparent reservoir pressure and the cumulative gas production volume in the different areas can be obtained, shown in Figures 5 and 6. Those curves illustrate obvious two-stage segmentation characteristics. (1)In the first stage, only the fluid, in the high-permeability zone, seepage to the bottom of the gas well. In the high-permeability development zone, the first decline characteristic of the curve is mainly affected by the reserves in the high-permeability development zone, and the decline rate is equal to . In the low-permeability development zone, the curve value is constant as the initial value(2)At the second stage, the reserves in the low-permeability area began to be recharged to the high-permeability area. Affected by the reserves in the low-permeability recharge zone, the decline rate closes to in the high-permeability zone. Compared with the traditional material balance method, the recoverable reserves evaluation result of this model is significantly lower, and the difference is , as shown in Figure 5. And the decline rate of low-permeability zone is equal to

3.2. Practical Application

Taking well W of a bio-reef gas reservoir as an example, the initial reservoir pressure is 68.69 MPa, and the average daily production rate is 470,000 m³/d. After productivity test, the open flow rate is 1,990,000m³/d, and the binomial productivity equation is: (1) curve fitting

Figure 7 illustrates the fitting results of this well and shows the relationship curve between and cumulative gas production volume in the high-permeability development area of this well. Applying the new multi-region material balance method to fit the values of different times, we can get relevant reservoir parameters of different regions. The evaluation results show that the reserves in the high-permeability development zone are equal to , the reserves in the low-permeability recharge zone are , the low-velocity non-Darcy seepage factor G in the low-permeability recharge zone is 0.008Mpa/m, and the seepage resistance coefficient A is equal to 0.2.

As shown in Figure 8, the cumulative gas volume supplied from low-permeability region behaves a linear increase trend. When the cumulative gas production reaches , the fluid in the low-permeability zone begins to supply to high-permeability. When the gas well produced natural gas, the low-permeability area could supply natural gas to the high-permeability area. (2)Gas well performance prediction

The lowest bottom hole flow pressure of this well is 16 MPa, the economic limit production is , and the gas rate of constant rate production period is . Based on the new multi-region material balance and traditional material balance methods, we can get the dynamic forecast curve of this well, as shown in Figure 9. Using the new multi-region material balance equation to predict the development performance, the constant rate production time is 10.3 years, which is 0.65 years lower than prediction results of the traditional material balance method.

According to the new multi-region material balance method proposed in this work, the change law of average formation pressure and supplying gas rate can be obtained, as shown in Figure 10. The curve characteristics can be divided into three stages: (1)First is the low-velocity non-Darcy seepage effect stage. The reservoir pressure and bottom hole flowing pressure , in the high-permeability development area, show rapid decline. However, affected by the low-velocity non-Darcy seepage law, the reservoir pressure () in the low-permeability recharge area remains constant, and the gas rate supplied from low-permeability region is equal to 0(2)Second is the later stage of constant rate production period. In this period, the pressure gradient between the low-permeability and the high-permeability area is greater than the low-velocity non-Darcy seepage factor G. With the cumulative gas volume increasing, the reservoir pressure () in the high-permeability area, the reservoir pressure () in the low-permeability area, the bottom hole pressure, and the cumulative gas volume (Q_C), supplying from the low-permeability zone to the high-permeability development zone, illustrates linear rapidly changing characteristics(3)Third is the stage of decline in production. During the production decline stage, the gas production rate decreases rapidly. The reservoir pressure in high-permeability area, the reservoir pressure in low-permeability area, and the bottom hole flowing pressure show a slow decline. At the same time, the cumulative gas volume supplying from the low-permeability zone to the high-permeability development zone shows a slow increase

In summary, as production progresses, the low-permeability replenishment area will gradually play an important role. Especially after entering the decline stage, the proportion of the supplied gas in the cumulative production gas volume is even greater. In addition, the total cumulative gas production volume is estimated to be , and the final recovery degree is 58.8%.

4. Summary and Conclusions

In this work, a new multi-region material balance method for dynamic prediction of bio-reef gas reservoir is presented, which comprehensively considers the heterogeneity of reefs, low-permeability reservoir supply, and initially introduced the effects of low-velocity non-Darcy seepage. And this method further combined the MMBE with binomial productivity equation for development performance prediction. In addition, practical application was also conducted. Then, some important conclusions are drawn as following: (1)Through fitting the curve of typical bio-reef gas well, dynamic reserves in high/low-permeability zones, the resistance coefficient, and the non-Darcy seepage factor in the low-permeability zone can be obtained. Compared with the traditional material balance method, this method can more precisely describe the development dynamics of the low-permeability zone(2)Combined the new multi-region material balance with productivity equations, we can achieve early dynamic prediction. The high/low-permeability zone reservoir pressure, bottom hole flowing pressure, and the recharge gas volume can be achieved by curves features in the constant rate production period(3)Compared with the traditional material balance method, the evaluation of the constant rate production time is relatively low. It can more accurately predict the reserves of the low-permeability zone, providing a theoretical basis for the determination of the production well working system, and the deployment/adjustment of the development well pattern

Appendix

A. Early seepage stage in high-permeability region

At the early seepage stage in high-permeability region, the material balance equation of the high-permeability development zone is:

Now with the definition of gas compressibility:

The chain rule for derivatives is:

For the spatial derivative of across the physical model of average thickness , Eq. (2) can be rewritten as:

Initial conditions:

The results for the seepage in the high-permeability are:

Then, with the real gas equation of state:

And the definition of cumulative gas produced:

Eq. (A.6) can be written in the form:

Where: where is molecular weight; is universal gas constant, MPa.m³/mol.K; is reservoir pressure, MPa; is initial reservoir pressure, MPa; is gas deviation factor, decimal; is pressure under standard conditions, MPa; is the reservoir temperature, K; is the temperature under standard conditions, K; is the fluid density in the high-permeability zone, kg/m³; is the cumulative gas production volume starting to supply from the low-permeability area to the high-permeability area, 10⁸m³.

B. Late replenishment stage in low-permeability region

Based on the one-dimensional non-linear equivalent analytical model of multi-region material balance, the material balance equation, at the late replenishment stage in low-permeability region, can be described approximately by the pair of equations:

sssHere Eqs. (A.11) and (A.12) describe the material balance in the high-permeability zone and the low-permeability zone separately.where is the fluid density in the high-permeability area, kg/m³; is the fluid density in the low-permeability area, kg/m³; is the equivalent permeability between low-permeability area and high-permeability area, mD; is the low-velocity non-Darcy seepage factor, MPa/m; is time, day; is the fluid density under standard conditions, kg/m³; is the gas production rate under standard conditions, 10⁴m³/day.

Initial conditions: where is the gas density under the original conditions, kg/m³.

Eq. (A.11)–Eq. (A.14) can be transformed into a second-order homogeneous linear differential equation with constant coefficients:

Among them, the parameters , , and are, respectively, expressed as: where is the natural gas compression factor, MPa^-1; is the seepage resistance coefficient, decimal.

From the above formula, we can get:

Nomenclature

:	The gas seepage velocity (m/s)
:	Reservoir permeability (mD)
:	Pressure difference across the core (MPa)
:	Core length (m)
:	Fluid viscosity (mPa·s;)
:	Molecular weight
:	The universal gas constant (MPa.m³/mol.K)
:	Reservoir pressure (MPa)
:	Gas deviation factor, decimal
:	The reservoir temperature (K)
:	The fluid density (kg/m³)
:	The cumulative gas production volume, (10⁸m³)
:	The low-velocity non-Darcy seepage factor (MPa/m)
:	Time (day)
:	The gas production rate (10⁴m³/day)
:	The natural gas compression factor, (MPa^-1)
:	The seepage resistance coefficient (decimal)
:	The first-order coefficient of the binomial capacity equation coefficients
:	The second-order coefficient of the binomial capacity equation coefficients.

Subscripts

i:	Gas reservoir initial conditions
:	The standard conditions
1:	The high-permeability zone
2:	The low-permeability zone
:	Natural gas
:	Bottom hole
:	Starting to supply gas from the low-permeability area to the high-permeability area, 10⁸m³.

Superscripts

:	the equivalent parameter between low-permeability area and high-permeability area.

Data Availability

The [Bio-reef sedimentary architecture map, pictures, pressure data] data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

S. Zhuangzhi, L. J. Du Zhimin, L. Yingchuan, and Z. Xu, “Prediction method of production performance of gas wells,” Natural Gas Industry, vol. 24, pp. 78–81, 2004.
View at: Google Scholar
Y. X. Wang, “Decline lows and performance prediction of gas well production,” Natural Gas Industry, vol. 15, no. 2, pp. 52–55, 1995.
View at: Google Scholar
R. A. Galecio and V. A. Huerta, “Traditional material-balance analysis from a stochastic point of view: predicting reservoir performance,” in Paper presented at the SPE Latin American and Caribbean Petroleum Engineering Conference, Lima, Peru, 2010.
View at: Publisher Site | Google Scholar
G. Penuela, A. Ordonez, and A. Bejarano, “A generalized material balance equation for coal seam gas reservoirs,” in Paper presented at the SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, 1998.
View at: Publisher Site | Google Scholar
D. W. Peaceman, Fundamentals of Numerical Reservoir Simulation, Elsevier, 2000.
T. Ertekin, J. H. Abou-Kassem, and G. R. King, Basic Applied Reservoir Simulation, Society of Petroleum Engineers, Richardson, TX, 2001.
A. H. El-Banbi, Analysis of Tight Gas Well Performance, Texas A&M University, 1998.
D. M. Anderson, M. Nobakht, S. Moghadam, and L. Mattar, “Analysis of Production Data from Fractured Shale Gas Wells,” in SPE unconventional gas conference, Pittsburgh, Pennsylvania, USA, 2010.
View at: Publisher Site | Google Scholar
G. Chengtai, “The material balance equation with a supplying region and its use in gas reservoir,” Petroleum Expoloration and Development, vol. 20, no. 5, pp. 53–61, 1993.
View at: Google Scholar
C. Gao and T. Lu, “Multi-region material balance method and its use in the performance prediction and optimal well pattern design of edge water reservoir,” Petroleum Exploration and Development, vol. 33, no. 1, pp. 103–106, 2006.
View at: Google Scholar
H. Yong, Z. Huayin, and C. Jianjun, “Physical simulated study on the gas-supplying mechanism of high-low permeable serial connection-formations,” Natural Gas Geoscience, vol. 18, no. 3, pp. 469–472, 2007.
View at: Google Scholar
C. Shiqing, Y. Xiuxiang, and X. I. E. Linfeng, “Multi-region calculation of dynamic reserve and pressure of constant volume reservoirs using material balance,” Petroleum Drilling Techniques, vol. 35, no. 3, pp. 66–68, 2007.
View at: Google Scholar
M. E. Lord and R. E. Collins, “Detecting compartmented gas reservoirs through production performance,” in SPE annual technical conference and exhibition, Dallas, Texas, 1991.
View at: Publisher Site | Google Scholar
J. Hagoort, Fundamentals of Gas Reservoir Engineering, Elsevier, 1988.
L. Mattar and D. Anderson, “Dynamic material balance (oil or gas-in-place without shut-ins),” in Canadian International Petroleum Conference, Calgary, Alberta, 2005.
View at: Google Scholar
Y. Ling, G. Chengtai, and G. Weixin, “Application of unbalance recovery to the gas field in central Shanganning basin,” Journal of Xi'an Petroleum Institute, vol. 14, pp. 13–15, 1999.
View at: Google Scholar
Y. Y. Zhou and X. W. Che, “Regional material balance method for low-permeability reservoir,” Journal of Xi'an Shiyou University, vol. 25, no. 6, pp. 42–46, 2010.
View at: Google Scholar
C. Guang-mei and C. H. E. Xiong-wei, “Domain material balance method for low-permeability gas-oil reservoir by depletion drive process,” Xinjiang Petroleum Geology, vol. 32, no. 2, pp. 157–159, 2011.
View at: Google Scholar
W. E. Yi, R. A. Qiquan, L. I. Ran, Y. U. Jiangru, and D. O. Jiaxin, “Determination of dynamic reserves of fractured horizontal wells in tight oil reservoirs by multi-region material balance method,” Petroleum Exploration and Development, vol. 43, no. 3, pp. 490–498, 2016.
View at: Publisher Site | Google Scholar
H. D. Sun, “A balance equation of gas materials with recharge capacity and the performance prediction,” Acta Petrolei Sinica, vol. 32, no. 4, pp. 683–686, 2011.
View at: Google Scholar
M. Dejam, H. Hassanzadeh, and Z. Chen, “Pre-Darcy flow in porous media,” Water Resources Research, vol. 53, no. 10, pp. 8187–8210, 2017.
View at: Publisher Site | Google Scholar
L. He, H. Mei, X. Hu, M. Dejam, Z. Kou, and M. Zhang, “Advanced flowing material balance to determine original gas in place of shale gas considering adsorption hysteresis,” SPE Reservoir Evaluation & Engineering, vol. 22, no. 4, pp. 1282–1292, 2019.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Yan Liu and Langtao Zhu. 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.

PDF Download Citation

Download other formats

Order printed copies

Views

180

Downloads

488

Citations