A New Method to Calculate the Recoverable Reserves and Recovery Ratio of Shale Gas Reservoir

Maolin Zhang, Yizhong Zhang^*, Long Yang, Haiyan Mei, Pengxiang Shen, Jingnan Ge

Electronics & Information College of Yangtze University, Wuhan, China.

DOI: 10.4236/ojogas.2017.24016                

Abstract

At present, the most commonly used method to calculate recoverable reserves of shale gas reservoirs is the production decline curve method, which is only applicable to the production wells with decreasing production. For the production wells with constant productions, the material balance method could be used to calculate the recoverable reserves, but it does not consider the particularity of shale gas reservoirs. Under this circumstance, this paper deduces a new calibration method for recoverable reserves of shale gas reservoirs with double porosity medium. The method is based on the material balance equation, which has a wide application range. The method cannot only calculate the recoverable reserves and recoveries of free gas in fractures, but also calculate recoverable reserves and recoveries of free gas and adsorbed gas in the matrix. The calculated results show that the free gas in the shale gas reservoir is mainly in the fractures. The recovery of the free gas in the fractures is relatively high and the recovery of the free gas and adsorbent gas in the matrix is rather low. Gas reservoir recovery and recoverable reserves are more sensitive to abandonment pressure, therefore, in the production process, the reduction of the abandonment pressure of gas reservoir should be attempted. The matrix porosity, the adsorption phase density and the crack compression coefficient have a greater effect on the recovery rate. And the study of these parameters needs to be strengthened in order to obtain more accurate parameters.

Keywords

Recoverable Reserves, Fracture (Matrix), Recovery Ratio, Controlled Reserves, Material Balance

1. Introduction

Shale gas reservoir is a kind of unconventional reservoir where the natural gas is stored in the dark shale or high-carbon shale and it is a typical in-situ accumulated reservoir [1] . The natural gas is stored in the reservoir formation in the form of free gas and adsorbed gas. Free gas is mainly distributed in the fracture pores and the matrix pores, while the adsorbed gas is mainly distributed on the inner surface of micro-porous matrix system [2] . Chen proposed pressure decline method, material balance method, and production decline method to calculate the conventional gas reservoir recoverable reserves [3] . Chen used the principle of material balance to obtain the recoverable reserves of coalbed gas reservoir with constant volume [4] . However, the methods are not applicable for shale gas reservoirs as status of gas in the coalbed gas reservoirs is different from that of shale gas reservoirs. Ilk proposed an exponential law reduction method that can be used to calculate the recoverable reserves of shale gas reservoirs [5] . Valko proposed an extended exponential reduction method for shale gas reservoirs, which is an empirical method [6] . Duong proposed empirical analysis of the decline in gas production in fractured shale gas reservoirs and calculated recoverable reserves using this method [7] . Duan found that the use of power-law exponential decay methods in the early production of gas wells can predict gas production, recoverable reserves and declining rates [8] . Bai found that shale gas production decline curve can be used to learn most of the hyperbolic decline [9] . Chen introduced the hyperbolic yield reduction model to predict recoverable reserves of shale gas reservoirs [10] . Cao predicted the recoverable reserves of shale gas reservoirs by using a power law exponential diminution model [11] .

The widely used recoverable reserves calculation method now is the pressure drop method and the production decline curve method. But the pressure drop method does not consider the particularity of the shale gas reservoir. Thus, the production decline curve method is only applicable to the well with gradually decrease output and it is only an empirical method. Therefore, from the view of mass conservation, this paper establishes the material balance equation of shale gas reservoir, using the equation the controlled reserves, the recoverable reserves and recovery rate could be obtained when the formation pressure reaches the abandonment pressure. Moreover, this paper can calculate the recovery of free gas in the matrix, free gas in the cracks and adsorbed gas in the matrix which does not require the production data and pressure data. It is very convenient for low porosity and low permeability shale gas reservoirs.

2. Methods

The development of the shale gas reservoir is an isothermal process, so the shale gas adsorption isotherm [12] is Equation (1).

$V_{\text{E}}=\frac{V_{\text{L}}p}{p_{\text{L}}+p}$ (1)

Adsorption gas is mainly adsorbed on the pore surface of the shale matrix [13] , which will cause that the pore size occupied by the gas is smaller than the actual pore size. Ambrose [14] and Cao [15] define the apparent porosity of the adsorption phase, which is shown in Equation (2).

${\varphi }_{\text{a}}=4.6545\times {10}^{-5}M\frac{{\rho }_{\text{h}}}{{\rho }_{\text{s}}}\left(V_{\text{L}}\frac{p}{p+p_{\text{L}}}\right)$ (2)

The reserves of free gas in the matrix at the initial formation pressure can be expressed as Equation (3).

$G_{\text{mi}}=\frac{V_{\text{m}}\left[{\varphi }_{\text{mi}}\left(1-s_{\text{mwi}}\right)-{\varphi }_{\text{ai}}\right]}{B_{\text{gi}}}$ (3)

At initial formation pressure, the reserves of free gas in the fracture can be expressed as Equation (4).

$G_{\text{fi}}=\frac{V_{\text{fp}}\left(1-s_{\text{fwi}}\right)}{B_{\text{gi}}}$ (4)

At the initial formation pressure, the reserves of adsorbed gas in the matrix are Equation (5).

$G_{\text{ai}}=V_{\text{m}}{\rho }_{\text{h}}{V}_{\text{Ei}}$ (5)

When the formation pressure decreases and desorption of the adsorbed gas takes place, the shale matrix shrinks and deforms which cause the increase of the matrix porosity. When the fluid flows out, the effective stress increases and the rock skeleton volume is compressed, which cause the matrix porosity to decrease. Therefore, the matrix porosity is affected by the effects of both shrinkage of matrix and increase of stress [16] .

The relationship between rock deformation and reservoir formation pressure could be explained by the Bangham [17] solid deformation theory and the Langmuir isothermal adsorption model, which is shown in Equation (6).

$\Delta \epsilon =\frac{{\rho }_{\text{h}}RT{V}_{\text{L}}}{E{V}_0}\times \left[\ln \left(p_{\text{L}}+p_{\text{i}}\right)-\ln \left(p_{\text{L}}+p\right)\right]$ (6)

The total deformation of the shale matrix is the sum of the matrix shrinkage and the elastic expansion of the rock, as shown in Equation (7).

$\Delta {\epsilon }_{\text{z}}=\frac{{\rho }_{\text{h}}RT{V}_{\text{L}}}{E{V}_0}\times \left[\ln \left(p_{\text{L}}+p_{\text{i}}\right)-\ln \left(p_{\text{L}}+p\right)\right]-c_{\text{m}}\left(p_{\text{i}}-p\right)$ (7)

Seidle (1985) [18] proposed an equation to explain the relationship between the matrix porosity and the deformation of the matrix, which is shown in Equation (8).

$\frac{{\varphi }_{\text{m}}}{{\varphi }_{\text{mi}}}=1+\left(1+\frac{2}{{\varphi }_{\text{mi}}}\right)\left\{\frac{{\rho }_{\text{h}}RT{V}_{\text{L}}}{E{V}_0}\times \left[\ln \left(p_{\text{L}}+p_{\text{i}}\right)-\ln \left(p_{\text{L}}+p\right)\right]-c_{\text{m}}\left(p_{\text{i}}-p\right)\right\}$ (8)

So the matrix porosity after deformation can be obtained from Equation (9).

${\varphi }_{\text{m}}=\left\{1+\left(1+\frac{2}{{\varphi }_{\text{mi}}}\right)\left\{\frac{{\rho }_{\text{h}}RT{V}_{\text{L}}}{E{V}_0}\times \left[\ln \left(p_{\text{L}}+p_{\text{i}}\right)-\ln \left(p_{\text{L}}+p\right)\right]-c_{\text{m}}\left(p_{\text{i}}-p\right)\right\}\right\}{\varphi }_{\text{mi}}$ (9)

As the pressure decreases, the rock particles deform and the initial water expands, and the change in the pore volume [19] in the fracture is given by Equation (10).

$\Delta V=V_{\text{fp}}\left[\left(c_{\text{f}}+c_{\text{w}}{s}_{\text{fwi}}\right)\left(p_{\text{i}}-p\right)\right]$ (10)

The definition of the coefficient of compressibility can be shown in Equation (11).

$\begin{array}{c}{c}_{\text{mx}}=-\frac{1}{V_{\text{mp}}}\frac{\Delta V_{\text{mp}}}{\Delta p}=\frac{{\varphi }_{\text{mi}}-{\varphi }_{\text{m}}}{{\varphi }_{\text{mi}}\left(p_{\text{i}}-p\right)}\\ =\left(1+\frac{2}{{\varphi }_i}\right)\left\{\frac{{\rho }_{\text{h}}RT{V}_L}{E{V}_0\left(p_{\text{i}}-p\right)}\times \left[\ln \left(p_{\text{L}}+p\right)-\ln \left(p_{\text{L}}+p_{\text{i}}\right)\right]+c_{\text{m}}\right\}\end{array}$ (11)

Based on the elastic deformation of matrix rock particles and initial water [20] , the material equilibrium equation of formation water is shown in Equation (12).

$V_{\text{mp}}\left[1-c_{\text{mx}}\left(p_{\text{i}}-p\right)\right]s_{\text{mw}}=V_{\text{mp}}{s}_{\text{mwi}}+V_{\text{mp}}{c}_{\text{w}}{s}_{\text{mwi}}\left(p_{\text{i}}-p\right)$ (12)

The irreducible water saturation in the matrix at the present formation pressure is shown in Equation (13).

$s_{\text{mw}}=\frac{s_{\text{mwi}}\left[1+c_{\text{w}}\left(p_{\text{i}}-p\right)\right]}{1-c_{\text{mx}}\left(p_{\text{i}}-p\right)}$ (13)

At the present formation pressure, the residual reserves of adsorbed gas in the matrix is Equation (14).

$G_{\text{ar}}=V_{\text{m}}{\rho }_{\text{h}}{V}_{\mathrm{E}}$ (14)

At the present formation pressure, the residual free gas reserves in the matrix can be expressed as Equation (15).

$G_{\text{mr}}=\frac{V_{\text{m}}\left[{\varphi }_{\text{m}}\left(1-s_{mw}\right)-{\varphi }_{\text{a}}\right]}{B_{\text{g}}}$ (15)

At the present formation pressure, the residual free gas reserves in the fracture type can be shown in Equation (16).

$G_{\text{fr}}=\frac{V_{\text{fp}}\left(1-s_{\text{fwi}}\right)-\Delta V}{B_{\text{g}}}$ (16)

From the conservation of mass, The sum of the original matrix free gas reserves, the original crack free gas reserves, the original matrix adsorbed gas reserves is equal to the sum of the remaining free reserves of the matrix, the remaining free reserves of the fissures, the residual gas reserves of the matrix and the cumulative production.

So the Equation (17) can be obtained when the Equations (3), (4), (5), (10), (14), (15) and (16) are obtained.

$\begin{array}{c}{G}_{\text{p}}=\frac{V_{\text{fp}}\left(1-s_{fwi}\right)}{B_{\text{gi}}}+\frac{V_{\text{m}}\left[{\varphi }_{\text{mi}}\left(1-s_{\text{mwi}}\right)-{\varphi }_{\text{ai}}\right]}{B_{\text{gi}}}+V_{\text{m}}{\rho }_{\text{h}}{V}_{\text{Ei}}-V_{\text{m}}{\rho }_{\text{h}}{V}_{\text{E}}\\ -\frac{V_{\text{m}}\left[{\varphi }_{\text{m}}\left(1-s_{\text{mw}}\right)-{\varphi }_{\text{a}}\right]}{B_{\mathrm{g}}}-\frac{V_{\text{fp}}\left(1-s_{\text{fwi}}\right)-V_{\text{fp}}\left[\left(c_{\text{f}}+c_{\text{w}}{s}_{\text{fwi}}\right)\left(p_{\text{i}}-p\right)\right]}{B_{\text{g}}}\end{array}$ (17)

Then the Equation (18) can be deduced.

$\begin{array}{c}{G}_{\text{p}}=V_{\text{fp}}\left(\frac{\left(1-s_{\text{fwi}}\right)}{B_{gi}}-\frac{\left(1-s_{\text{fwi}}\right)}{B_{\text{g}}}+\frac{\left(c_{\text{f}}+c_{\text{w}}{s}_{fwi}\right)\left(p_{\text{i}}-p\right)}{B_{\text{g}}}\right)\\ +V_{\text{m}}\left(\frac{\left[{\varphi }_{\text{mi}}\left(1-s_{\text{mwi}}\right)-{\varphi }_{\text{ai}}\right]}{B_{\text{gi}}}+{\rho }_h{V}_{Ei}-{\rho }_{\text{h}}{V}_{\text{E}}-\frac{\left[{\varphi }_{\text{m}}\left(1-s_{\text{mw}}\right)-\varphi \right]}{B_g}\right)\end{array}$ (18)

Substituting Equations (3) and (4) into Equation (18) will deduce Equation (19).

$\begin{array}{c}{G}_{\text{p}}=\frac{G_{\text{fi}}{B}_{\text{gi}}}{\left(1-s_{\text{fwi}}\right)}\left(\frac{\left(1-s_{\text{fwi}}\right)}{B_{\text{gi}}}-\frac{\left(1-s_{\text{fwi}}\right)}{B_{\text{g}}}+\frac{\left(c_{\text{f}}+c_{\text{w}}{s}_{\text{fwi}}\right)\left(p_i-p\right)}{B_{\text{g}}}\right)\\ +\frac{G_{\text{mi}}{B}_{\text{gi}}}{\left[{\varphi }_{\text{mi}}\left(1-s_{\text{mwi}}\right)-{\varphi }_{ai}\right]}\left(\frac{\left[{\varphi }_{\text{mi}}\left(1-s_{\text{mwi}}\right)-{\varphi }_{\text{ai}}\right]}{B_{\text{gi}}}+{\rho }_{\mathrm{h}}{V}_{\text{Ei}}-{\rho }_{\text{h}}{V}_{\text{E}}-\frac{\left[{\varphi }_{\text{m}}\left(1-s_{\text{mw}}\right)-{\varphi }_{\text{a}}\right]}{B_{\text{g}}}\right)\end{array}$ (19)

Now some symbols can be defined to simplify the Equation (19).

$F=\frac{B_{\text{gi}}}{\left[{\varphi }_{\text{mi}}\left(1-s_{\text{mwi}}\right)-{\varphi }_{\text{ai}}\right]}\left(\frac{\left[{\varphi }_{\text{mi}}\left(1-s_{\text{mwi}}\right)-{\varphi }_{\text{ai}}\right]}{B_{gi}}+{\rho }_h{V}_{Ei}-{\rho }_{\text{h}}{V}_{\text{E}}-\frac{\left[{\varphi }_{\text{m}}\left(1-s_{\text{mw}}\right)-{\varphi }_{\text{a}}\right]}{B_{\text{g}}}\right)$

$Y=G_{\text{p}}$

$M=\frac{B_{\text{gi}}}{\left(1-s_{\text{fwi}}\right)}\left(\frac{\left(1-s_{fwi}\right)}{B_{\text{gi}}}-\frac{\left(1-s_{\text{fwi}}\right)}{B_{\text{g}}}+\frac{\left(c_{\text{f}}+c_{\text{w}}{s}_{\text{fwi}}\right)\left(p_i-p\right)}{B_{\mathrm{g}}}\right)$ (20)

Substituting Equation (20) into Equation (19), the Equation (21) can be deduced.

$\frac{Y}{M}=G_{\text{fi}}+G_{\text{mi}}\frac{F}{M}$ (21)

From the above equations, a straight curve could be obtained if the value of Y/M is plotted against F/M in the Cartesian coordinate system. The intercept of this curve is the free gas reserves in the fracture and the slope of this curve is the free gas reserves in the matrix. The total reserves controlled could be expressed by Equation (22).

$G=G_{\text{fi}}+G_{\text{mi}}+\frac{G_{\text{mi}}{B}_{\text{gi}}}{\left[{\varphi }_{\text{mi}}\left(1-s_{mwi}\right)-{\varphi }_{\text{ai}}\right]}{\rho }_{\text{h}}{V}_{\text{Ei}}$ (22)

When the formation pressure reaches the abandonment pressure, the well is shut down and the production is stopped. The cumulative production amount at this time is the recoverable reserves of the reservoir. Substituting Equations (3), (4) and (5) into Equation (17) will yield Equation (23).

$\begin{array}{c}{G}_{\text{p}}=\frac{G_{\text{mi}}{B}_{\text{gi}}}{\left[{\varphi }_{\text{mi}}\left(1-s_{\text{mwi}}\right)-{\varphi }_{\text{ai}}\right]}\left(\frac{\left[{\varphi }_{\text{mi}}\left(1-s_{\text{mwi}}\right)-{\varphi }_{\text{ai}}\right]}{B_{\text{gi}}}-\frac{\left[{\varphi }_{\text{m}}\left(1-s_{\text{mw}}\right)-{\varphi }_{\text{a}}\right]}{B_{\text{g}}}\right)\\ +\frac{G_{\text{fi}}{B}_{\text{gi}}}{\left(1-s_{\text{fwi}}\right)}\left(\frac{\left(1-s_{fwi}\right)}{B_{gi}}-\frac{\left(1-s_{\text{fwi}}\right)}{B_g}+\frac{\left(c_f+c_{\text{w}}{s}_{\text{fwi}}\right)\left(p_{\text{i}}-p\right)}{B_{\text{g}}}\right)+\frac{G_{\text{ai}}}{{\rho }_{\text{h}}{V}_{\text{Ei}}}\left({\rho }_{\text{h}}{V}_{\text{Ei}}-{\rho }_{\text{h}}{V}_{\text{E}}\right)\end{array}$ (23)

And because the volume factor of the gas at present formation pressure can be expressed by Equation (24).

$B_g=\frac{zT{p}_{\text{sc}}}{z_{\text{sc}}{T}_{\text{sc}}p}$ (24)

When the formation pressure reaches the abandonment pressure, substituting Equations (1) and (24) into Equation (23) will yield Equation (25).

$\begin{array}{c}{G}_{\text{R}}=G_{\text{mi}}\left(1-\frac{\left[{\varphi }_{\text{m}}\left(1-s_{\text{mw}}\right)-{\varphi }_{\mathrm{a}}\right]}{\left[{\varphi }_{mi}\left(1-s_{\text{mwi}}\right)-{\varphi }_{\text{ai}}\right]}\frac{p_{\text{a}}/z_{\text{a}}}{p_{\text{i}}/z_{\text{i}}}\right)\\ +G_{\text{fi}}\left(1-\frac{p_{\text{a}}/z_{\text{a}}}{p_{\text{i}}/z_{\text{i}}}\left(1-\frac{\left(c_{\text{f}}+c_{\text{w}}{s}_{\text{fwi}}\right)\left(p_{\text{i}}-p_{\text{a}}\right)}{\left(1-s_{\text{fwi}}\right)}\right)\right)+G_{\text{ai}}\left(\frac{p_{\text{L}}}{p_{\text{i}}}\frac{p_{\text{i}}-p_{\text{a}}}{p_{\text{L}}+p_{\text{a}}}\right)\end{array}$ (25)

The recovery of free gas in the matrix can be expressed by Equation (26).

$E_{Rm}=1-\frac{\left[{\varphi }_{\text{m}}\left(1-s_{\text{mw}}\right)-{\varphi }_{\text{a}}\right]}{\left[{\varphi }_{\text{mi}}\left(1-s_{\text{mwi}}\right)-{\varphi }_{\text{ai}}\right]}\frac{p_{\text{a}}/z_{\text{a}}}{p_{\text{i}}/z_{\text{i}}}$ (26)

The recoverable reserves of free gas in the matrix can be expressed by Equation (27).

$G_{\text{Rm}}=G_{\text{mi}}{E}_{\text{Rm}}$ (27)

The recovery of free gas in the fracture can be expressed by Equation (28).

$E_{\text{Rf}}=1-\frac{p_{\text{a}}/z_{\text{a}}}{p_{\text{i}}/z_{\text{i}}}\left(1-\frac{\left(c_f+c_{\text{w}}{s}_{fwi}\right)\left(p_{\text{i}}-p_{\text{a}}\right)}{\left(1-s_{fwi}\right)}\right)$ (28)

The recoverable reserves of free gas in fractures can be expressed by Equation (29).

$G_{\text{Rf}}=G_{\text{fi}}{E}_{\text{Rf}}$ (29)

The recovery of adsorbed gas in the matrix can be expressed by Equation (30).

$E_{\text{Ra}}=\frac{p_{\mathrm{L}}}{p_{\text{i}}}\frac{p_{\text{i}}-p_{\text{a}}}{p_L+p_{\text{a}}}$ (30)

The recoverable reserves of adsorbed gas in the matrix can be expressed by Equation (31).

$G_{\text{Ra}}=G_{\text{ai}}{E}_{Ra}$ (31)

3. Results

The basic formation parameters of a certain well of a shale gas reservoir are listed as follows.

$p_{\text{i}}=24.13\text{MPa}$ , $B_{\text{gi}}=4.82\times {10}^{-3}$ , $s_{\text{mwi}}=0.25$ , $c_{\text{m}}=4.35\times {10}^{-4}{\text{MPa}}^{-\text{1}}$ , $c_{\text{f}}=2.04\times {10}^{-2}{\text{MPa}}^{-\text{1}}$ , $c_{\text{w}}=4.35\times {10}^{-4}{\text{MPa}}^{-\text{1}}$ , $S_{\text{fwi}}=0$ , ${\phi }_{\text{i}}=0.112$ , $V_{\text{L}}=11.32{\text{m}}^3/\text{t}$ , $p_L=2.41\text{MPa}$ , ${\rho }_s=0.34\text{g}/{\text{cm}}^{\text{3}}$ , ${\rho }_{\text{h}}=2.65\text{g}/{\text{cm}}^{\text{3}}$ , $E=26800\text{MPa}$ , $T=366.49\text{K}$ , $R=8.134\times {10}^{-6}\text{MPa}/\left(\text{K}\cdot \text{mol}\right)$ , $V_0=0.0245{\text{m}}^3/\text{mol}$ , $M=20\text{g}/\text{mol}$ , $p_{\text{a}}=3\text{MPa}$

The production data are shown in Table 1.

Table 1. The production data for shale gas reservoir.

Based on the production data in Table 1, the single well reserve, recoverable reserve and recovery could be calculated by using the material balance equation (see Figure 1). The results are shown in Table 2.

Figure 1. The results for this article using material balance equation.

Table 2 shows that the free gas mainly exists in the fractures, as the calculated recovery of free gas in the fracture is 93.37%, and the calculated recovery of free gas in the matrix is 56.06%. The calculated recovery of adsorbed gas in the matrix is the lowest, which is 39.04% and the total recovery of the gas is 62.13%.

Now assuming different abandonment pressures and calculating the recovery and recoverable reserves under different waste pressures, the results are shown in Table 3.

Table 3. The results of recovery for different abandonment pressure calculation.

Table 3 shows the calculation results by applying different abandonment pressure and it could be observed that as the abandonment pressure decreases, both of the recovery of the gas and recoverable reserves increase, and they are sensitive to the abandonment pressure.

Assuming different matrix porosity, then the recovery, controlled reserves and recoverable reserves under different matrix porosity can be calculated. And the results are shown in Table 4.

Table 4. The calculation results of recovery for different matrix porosity.

Table 4 shows that with the increase of matrix porosity, the recovery of free gas in the matrix and the total recovery are increasing, while the recovery of free gas in the cracks and the recovery of adsorbed gas in the matrix are unchanged.

In addition, with the increase of matrix porosity, the controlled reserves of free gas in the matrix are increasing, the recoverable reserves of free gas in the matrix are increasing, the controlled reserves of free gas in the cracks are decreasing, the recoverable reserves of free gas in the cracks are also decreasing, the controlled reserves of the adsorbed gas in the matrix are decreasing, the recoverable reserves of the adsorbed gas in the matrix are decreasing, the total controlled reserves are reduced and the total recoverable reserves are also reduced.

Therefore, the porosity of the matrix not only affects free gas recovery in the matrix and total recovery, but also has a great influence on the calculated reserves and recoverable reserves of each part. But there are no effects on the recovery of free gas in the fracture and the recovery of adsorbed gas in the matrix.

Now assuming different density of the adsorption phase and calculating the recovery, controlled reserves and re-coverable reserves under different density of the adsorption phase, the results are shown in Table 5.

Table 5. The calculation results of recovery for different density of the adsorption phase.

Table 5 shows that with the increase of adsorption phase density, the recovery of free gas in the matrix and the total recovery are increasing, while the recovery of free gas in the cracks and the recovery of adsorbed gas in the matrix are unchanged.

In addition, with the increase of adsorption phase density, the controlled reserves of free gas in the matrix are in-creasing, the recoverable reserves of free gas in the matrix are increasing, the controlled reserves of free gas in the cracks are decreasing, the recoverable reserves of free gas in the cracks are also decreasing, the controlled reserves of the adsorbed gas in the matrix are decreasing, the recoverable reserves of the adsorbed gas in the matrix are decreasing, the total controlled reserves are reduced and the total recoverable reserves are also reduced.

Therefore, the density of the adsorption phase not only affects the free gas recovery in the matrix and total recovery, but also has a great influence on the calculated reserves and recoverable reserves of each part. But there are no effects on the recovery of free gas in the fracture and the recovery of adsorbed gas in the matrix.

Now assuming different coefficient of compressibility of the fracture and calculating the recovery, controlled re-serves and recoverable reserves under different coefficient of compressibility of the fracture, the results are shown in Table 6.

Table 6. The calculation results of recovery for different coefficient of compressibility of the fracture.

Table 6 shows that with the increase of coefficient of compressibility in the fracture, the recovery of free gas in the cracks is increasing, while the recovery of free gas in the cracks and the recovery of adsorbed gas in the matrix are un-changed, the total recovery is decreasing.

In addition, with the increase of coefficient of compressibility in the fracture, the controlled reserves of free gas in the matrix are increasing, the recoverable reserves of free gas in the matrix are increasing, the controlled reserves of free gas in the cracks are decreasing, the recoverable reserves of free gas in the cracks are also decreasing, the controlled reserves of the adsorbed gas in the matrix are increasing, the recoverable reserves of the adsorbed gas in the matrix are increasing, the total controlled reserves are increasing and the total recoverable reserves are also increasing.

Therefore, the coefficient of compressibility in the fracture not only affects the free gas recovery in the cracks and total recovery, but also has a great influence on the calculated reserves and recoverable reserves of each part. But there are no effects on the recovery of free gas in the matrix and the recovery of adsorbed gas in the matrix.

4. Conclusions

In this paper, a calibration method of controlled reserves, recoverable reserves and recovery of shale gas reservoirs is established, which is widely applicable. It is shown from a realistic example that this method is simple and practically applicable. The calculated results show that the free gas in the shale gas reservoir mainly exists in the fracture and the proportion of the free gas in the matrix is small. But the adsorbed gas occupies most of the reserves. The recovery of adsorbed gas in the matrix is low, so in the later development, the matrix adsorption gas has great potential.

It also reveals that the recoverable reserves and recoveries are very sensitive to the abandonment pressure. Therefore, during the development process, the gas flow resistance should be minimized and the abandonment pressure should be reduced to obtain higher recovery and recoverable reserves. Besides, the porosity in matrix and the adsorption phase density affect the recovery of free gas in the matrix, and the fracture compression coefficient affects the recovery of free gas in the fracture. These three parameters have a significant effect on the calculation of reserves and recoverable reserves. Therefore, it is necessary to strengthen the study of these three parameters in future work in order to obtain more accurate values for the calculation of reserves.

Nomenclature

$p$ is the reservoir formation pressure, MPa;

$p_{\text{i}}$ is the reservoir initial formation pressure, MPa;

$V_{\text{E}}$ is the adsorption amount under the reservoir formation pressure, m³/t;

$V_{\text{L}}$ is the langmuir volume, m³/t;

$p_{\text{L}}$ is the langmuir pressure, MPa;

${\varphi }_{\text{a}}$ is the porosity of adsorption phase, f;

$M$ is the molar weight, g/mol;

${\rho }_{\text{h}}$ is the density of the rock, g/cm³;

${\rho }_{\text{s}}$ is the density of the adsorption phase, g/cm³;

$B_{\text{gi}}$ is the volume factor of the gas at the initial formation pressure, m³/m³;

$s_{\text{mwi}}$ is the initial water saturation of the matrix at the initial formation pressure, f;

${\varphi }_{\text{mi}}$ is the matrix porosity at the initial formation pressure, f;

${\varphi }_{\text{ai}}$ is the porosity of adsorption phase at the initial formation pressure, f;

$G_{\text{mi}}$ is the reserves of the free gas in the matrix at the initial formation pressure, 10⁸ m³;

$V_{\text{m}}$ is the total volume of the matrix, 10⁸ m³;

$s_{\text{fwi}}$ is the initial water saturation in the fracture, f;

$G_{\text{fi}}$ is the reserves of free gas in the fracture at initial pressure,10⁸ m³;

$V_{\text{fp}}$ is the fracture pore volume, 10⁸ m³;

$G_{\text{ai}}$ is the reserves of adsorbed gas at the initial formation pressure, 10⁸ m³;

$V_{\text{Ei}}$ is the adsorption amount at the initial formation pressure, m³/t;

$\Delta \epsilon$ is the deformation of the rock caused by the desorption of gas, dimensionless;

$R$ is the gas constant, MPa/(K・mol);

$T$ is the temperature of the reservoir formation, K;

$V_0$ is the molar based specific volume, m³/mol;

$E$ is the Young’s modulus, MPa;

$\Delta {\epsilon }_z$ is the elastic expansion of the rock under effective stress and matrix shrinkage, dimensionless;

$c_{\text{m}}$ is the elastic coefficient of compressibility for the rock, MPa⁻¹;

${\varphi }_{\text{m}}$ is the matrix porosity after deformation;

$\Delta V$ is the change of the pore volume in the fracture, 10⁸ m³;

$c_{\text{f}}$ is the coefficient of compressibility of the fracture, MPa⁻¹;

$c_{\text{w}}$ is the coefficient of compressibility of the irreducible water, MPa⁻¹;

$V_{\text{mp}}$ is the Matrix pore volume, 10⁸ m³;

$s_{\text{mw}}$ is the saturation of bound water under pressure, f;

$G_{\text{ar}}$ is the remaining reserves of adsorbed gas in the matrix at present formation pressure, 10⁸ m³;

$G_{fr}$ is the remaining reserves of free gas in the fracture under present formation pressure, 10⁸ m³;

$G_{\text{R}}$ is the recoverable reserves at abandoned pressure, 10⁸ m³;

$B_{\text{g}}$ is the volume factor of the gas at present formation pressure, m³/m³;

$p_{\text{a}}$ is the abandoned pressure, MPa;

$E_{\text{Rm}}$ is the. recovery of free gas in matrix, f;

$G_{Rm}$ is the recoverable reserves of free gas in the matrix, 10⁸ m³;

$E_{\text{Rf}}$ is the recovery of free gas in fractures, f;

$G_{\text{Rf}}$ is the recoverable reserves of free gas in fractures, 10⁸ m³;

$G_{\text{Ra}}$ is the recoverable gas in the matrix, 108 m³.

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	Zhang, J., Jin, Z. and Yuan, M. (2004) Reservoiring Mechanism of Shale Gas and Its Distribution. Natural Gas Industry, 24, 15-18.
[2]	Hay, D. (2012) Unconventional Hydrocarbon Resources in China and the Prospect of Exploration and Development. Petroleum Exploration and Development, 39, 139-146. https://doi.org/10.1016/S1876-3804(12)60026-3
[3]	Chen, Y. (2000) Calculation Methods of Recoverable Reserves of Oilfields. Xinjiang Petroleum Geology, 5, 15-24.
[4]	Chen, Y. and Hu, J. (2008) Derivation of Methods for Estimating OGIP and Recoverable Reserves and Recovery Ratio of Saturated Coal-Seam Gas Reservoirs. Oil & Gas Geology, 29, 151-156.
[5]	Ilk, D., Rushing, J. A., Perego, A.D., et al. (2008) Exponential vs. Hyperbolic Decline in Tight Gas Sands: Understanding the Origin and Implications for Reserve Estimates Using Arps’ Decline Curves. Society of Petroleum Engineers, Dallas.
[6]	Valkó, P.P. (2009) Assigning Value to Stimulation in the Barnett Shale: A Simultaneous Analysis of 7000 Plus Production Hystories and Well Completion Records. Society of Petroleum Engineers, Dallas.
[7]	Duong, A.N. (2010) An Unconventional Rate Decline Approach for Tight and Fracture-Dominated Gas Wells. Society of Petroleum Engineers, Dallas. https://doi.org/10.2118/137748-MS
[8]	Duan, Y.G., Cao, T.K., Wang, R. and Wei, M.Q. (2013) Analysis of Power Factor Index of Shale Gas Production. Journal of Southwest Petroleum University (Natural Science Edition), 35, 172-176.
[9]	Bai, Y.H., Yang, H. and Chen, G.H. (2013) Method for Determining Key Parameters in Typical Curve of Shale Gas Production. Special Oil and Gas Reservoirs, 20, 65-68.
[10]	Chen, Y.Q., Li, J., Qi, Y.D. and Zhang, X. (2014) Determination Methods of Geological Resources, Recoverable Resources and Well Controlled Recoverable Reserves in Shale Gas Reservoir. Xin Jiang Petroleum Geology, 35, 547-551.
[11]	Cao, A., Cao, G.J. and Zhang, X. (2015) Prediction of Recoverable Reserves in Unconventional Gas Reservoirs by the Power-Law Loss Ratio Rate Decline Relation. Liaoning Chemical Industry, 44, 315-317.
[12]	Fu, X.C., Shen, W.X. and Yao, T.Y. (2006) Physical Chemistry. Higher Education Press, Beijing, 360-375.
[13]	Xiong, Y., Xiong, W.L. and Liu, Q.G. (2015) Calculation of Shale Gas Reserves Based on the Volume of Adsorbed Gas. Geological Science and Technology Information, 34, 139-143.
[14]	Ambrose, R.J., Hartman, R.C., Diaz Campos, M., Akkutlu, I.Y. and Sondergeld, C. (2010) New Pore-Scale Considerations for Shale Gas in Place Calculations. In SPE Unconventional Gas Conference. Society of Petroleum Engineers, Dallas. https://doi.org/10.2118/131772-MS
[15]	Cao, T., Duan, Y., Wang, R., Fang, Q. and Wei, M. (2011) Evaluation of Shale Gas Reserves Considering Multi-Component Adsorption. Reservoir Evaluation & Development, 1, 68-70.
[16]	Liu, B., Yin, H., Wang, X. and Li, Q. (2013) Material Balance Equation with Revised Rock Compressibility for Shale Gas Reserve Calculation. Oil & Gas Geology, 34, 471-474.
[17]	Tan, M.H. and Huang, Y.Y. (1985) Superficial Physical Chemistry. China Building Industry Press, Beijing, 50-54.
[18]	Seidle, J.R. and Huitt, L.G. (1995) Experimental Measurement of Coal Matrix Shrinkage Due to Gas Desorption and Implications for Cleat Permeability Increases. In International Meeting on Petroleum Engineering. Society of Petroleum Engineers, Dallas. https://doi.org/10.2118/30010-MS
[19]	Zhang, L., Chen, G., Zhao, Y., Liu, Q. and Zhang, H. (2013) A Modified Material Balance Equation for Shale Gas Reservoir and a Calculation Method of Shale Gas Reserves. Gas Industry, 33, 66-70.
[20]	King, G.R. (1993) Material-Balance Techniques for Coal-Seam and Devonian Shale Gas Reservoirs with Limited Water Influx. SPE Reservoir Engineering, 8, 67-72. https://doi.org/10.2118/20730-PA