Method for determining recovery factor of deep coal bed methane wells

Abstract

A method for determining a recovery factor of a deep coal bed methane (CBM) well includes steps of: (1) collecting an original formation pressure and pressure-volume-temperature (PVT) experimental data of the CBM well, and establishing a relationship table between pressure and deviation factor; (2) obtaining isothermal adsorption experimental data, an original adsorbed gas content and an original free gas content for the deep coal rock; (3) according to the isothermal adsorption experimental data, determining Langmuir's pressure and Langmuir's volume; (4) determining an abandoned formation pressure in CBM well exploitation; (5) obtaining a deviation factor at the original formation pressure and a deviation factor at the abandoned formation pressure; and (6) according to the original formation pressure, the original adsorbed gas content, the original free gas content, the Langmuir's pressure, the abandoned formation pressure, the above deviation factors, the recovery factor of the CBM well is calculated.

Description

CROSS REFERENCE OF RELATED APPLICATION

The present invention claims priority under 35 U.S.C. 119 (a-d) to CN 202311456103.4, filed Nov. 3, 2023.

BACKGROUND OF THE PRESENT INVENTION

Field of Invention

The present invention relates to the oil and gas field development research, and more particularly to a method for determining a recovery factor of a deep coal bed methane well.

Description of Related Arts

In recent years, the development technology of deep coal bed methane (CBM) has made a breakthrough, and deep CBM has become an important energy resource like conventional natural gas. Practice shows that there are both free gas and adsorbed gas in deep CBM. At present, the prediction methods of the recovery factor of CBM mainly include analogy method, desorption method, isothermal adsorption curve method and gas reservoir numerical simulation method. Among them, the analogy method is more arbitrary and less reliable, which is mainly used in the early development stage of lack of basic data. The desorption method and the isothermal adsorption curve method are widely used in the development of shallow CBM, but these two methods are only able to be used to determine the recovery factor of adsorbed gas, but are unable to consider the influence of free gas on the recovery factor of CBM, so they are not suitable for determining the recovery factor of the deep CBM in which free gas and adsorbed gas coexist. The gas reservoir numerical simulation method is able to comprehensively consider the influence of heterogeneous geological characteristics and seepage characteristics of coal bed. However, this method needs to collect a large number of data. The research process is complex and the application is difficult. Moreover, the reliability of the simulation results is greatly affected by the understanding of the researchers on geological conditions of the target area, and the experience and technical level of the researchers. Therefore, it is urgent to provide a new method for determining a recovery factor of a deep coal bed methane well.

SUMMARY OF THE PRESENT INVENTION

An object of the present invention is to provide a method for determining a recovery factor of a deep coal bed methane (CBM) well. The method is able to solve the problem that the existing technology does not consider the influence of free gas on the recovery factor of CBM, and is not suitable for the determination of the recovery factor of the deep CBM well in which free gas and adsorbed gas coexist.

Accordingly, a method for determining a recovery factor of a deep coal bed methane (CBM) well comprises steps of:

• • (step 1) collecting an original formation pressure P_i and pressure-volume-temperature (PVT) experimental data of the CBM well, obtaining deviation factors corresponding to different pressures at a formation temperature of CBM, and establishing a relationship table between the pressures and the deviation factors;
  • (step 2) obtaining isothermal adsorption experimental data, an adsorbed gas content G_0ad and a free gas content G_0f at the original formation pressure and an original temperature for the deep coal rock;
  • (step 3) according to the isothermal adsorption experimental data, determining Langmuir's pressure P_L and Langmuir's volume V_L;
  • (step 4) according to a depth of burial of the CBM well, determining an abandoned formation pressure P_ab in CBM well exploitation;
  • (step 5) according to the relationship table between the pressures and the deviation factors in the (step 1), obtaining a deviation factor Z_i at the original formation pressure P_i and a deviation factor Z_ab at the abandoned formation pressure P_ab; and (step 6) according to the original formation pressure P_i in the (step 1), the adsorbed gas content G_0ad and the free gas content G_0f at the original formation pressure and the original temperature in the (step 2), the Langmuir's pressure P_L in the (step 3), the abandoned formation pressure P_ab in the (step 4), the deviation factor Z_i at the original formation pressure P_i and the deviation factor Z_ab at the abandoned formation pressure P_ab in the (step 5), by a model of

$R=\frac{G_{\mathrm{of}}\times \left(1-\frac{P_{\mathrm{ab}}{Z}_i}{P_i{Z}_{\mathrm{ab}}}\right)+G_{\mathrm{oad}}\times \left(1-\frac{P_{\mathrm{ab}}}{P_i}\times \frac{P_i\times P_L}{P_{\mathrm{ab}}+P_L}\right)}{G_{\mathrm{oad}}+G_{\mathrm{of}}},$

• • calculating the recovery factor of the CBM well.

Preferably, the (step 3) specifically comprises:

$\mathrm{taking}{y}_{\left(i\right)}=\frac{1}{V_{g\left(i\right)}}\mathrm{and}{x}_{\left(i\right)}=\frac{1}{P_{\left(i\right)}},\mathrm{wherein}i=1,2,\dots ,n;$

• • according to adsorbed gas volumes V_g(1), V_g(2), . . . , V_g(n) corresponding to different pressures P₍₁₎, P₍₂₎, . . . , P_(n) at the formation temperature of the CBM, obtaining a series of observation points (y_(i), x_(i)); and
  • obtaining a line equation by performing linear fitting on the observation points, wherein the Langmuir's volume V_L is equal to a reciprocal of an intercept of the line equation, and the Langmuir's pressure P_L is equal to a slope of the line equation divided by the intercept of the line equation.

Preferably, in the (step 4), the abandoned formation pressure P_ab in CBM well exploitation is determined by analog method or empirical formula method.

Preferably, in the (step 5), according to the relationship table between the pressures and the deviation factors in the (step 1), the deviation factor Z_i at the original formation pressure P_i and the deviation factor Z_ab at the abandoned formation pressure P_ab are obtained by interpolation method;

or according to the relationship table between the pressures and the deviation factors in the (step 1), a fitted function relationship Z=f(P) between the deviation factors and the pressures is established by taking the deviation factors as a dependent variable and the pressures as an independent variable, and then by the fitted function relationship, obtaining the deviation factor Z_i=f(P_i) at the original formation pressure P_i, and the deviation factor Z_ab=f(P_ab) at the abandoned formation pressure P_ab.

Beneficial effects as follows.

(1) Based on the PVT experimental data and the isothermal adsorption experimental data of the CBM, the present invention provides a new method for determining a recovery factor of a deep coal bed methane (CBM) well, which is able to effectively solving the problems of the analog method, desorption method, isothermal adsorption curve method and gas reservoir numerical simulation method. The method provided by the present invention not only considers the influence of free gas and adsorbed gas exploitation on the recovery factor of CBM, but also avoids the application difficulty of gas reservoir numerical simulation method. It is simple and convenient in application, easy to understand and realize, has strong operability, is effective and practical, and has good popularization and application value.

(2) The recovery factor is an essential index for development plan preparation, development benefit evaluation and development feasibility demonstration of CBM, so the determination of recovery factor has very important practical value in the mine field.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a flow chart of a method for determining the recovery factor of deep coal bed methane wells according to a preferred embodiment of the present invention.

FIG. 2 is a linear fitting diagram of observation point data in the process of determining Langmuir's pressure and Langmuir's volume.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT

Referring to FIG. 1, a flow chart of a method for determining a recovery factor of a deep coal bed methane (CBM) well according to a preferred embodiment of the present invention is illustrated. The method comprises steps of:

(step 1) collecting an original formation pressure P_i and pressure-volume-temperature (PVT) experimental data of the CBM well, obtaining deviation factors corresponding to different pressures at a formation temperature of CBM, and establishing a relationship table between the different pressures and the deviation factors;

(step 2) obtaining isothermal adsorption experimental data, an adsorbed gas content G_0ad and a free gas content G_0f at the original formation pressure and an original temperature for the deep coal rock;

(step 3) according to the isothermal adsorption experimental data, determining Langmuir's pressure P_L and Langmuir's volume V_L, wherein:

• • an isothermal adsorption equation of CBM is

$\begin{array}{cc}{V}_g=\frac{V_L\times P}{P+P_L},& \left(1\right)\end{array}$

• • by taking the reciprocal at both sides of the equation (1), an equation (2) is obtained as follows

$\begin{array}{cc}\frac{1}{V_g}=\frac{P+P_L}{V_L\times P},& \left(2\right)\end{array}$

• • and then based on the equation (2), an equation (3) is obtained as follows

$\begin{array}{cc}\frac{1}{V_g}=\frac{1}{V_L}+\frac{P_L}{V_L}\times \frac{1}{P};& \left(3\right)\\ \mathrm{taking}{y}_{\left(i\right)}=\frac{1}{V_{g\left(i\right)}},\mathrm{and}& \left(4\right)\\ x_{\left(i\right)}=\frac{1}{P_{\left(i\right)}},& \left(5\right)\end{array}$

• • here, i is isothermal adsorption experimental data point number, i=1, 2, . . . , n, V_g(i) (m³/t) is adsorbed gas volume at the i^th experimental data point, P_(i) (MPa) is pressure at the i^th experimental data point, V_L (m³/t) is Langmuir's volume, P_L (MPa) is Langmuir's pressure,
  • according to adsorbed gas volumes V_g(1), V_g(2), . . . , V_g(n) corresponding to different pressures P₍₁₎, P₍₂₎, . . . , P_(n) at the formation temperature of CBM in the isothermal adsorption experimental data, a series of observation points (y_(i), x_(i)) are obtained by the equations (4) and (5), linear fitting is performed on the observation points, wherein according to the equation (3), it is able to be known that the Langmuir's volume V_L is reciprocal of an intercept of a line equation obtained by linear fitting, and the Langmuir's pressure P_L is equal to a slope of the fitted line equation divided by the intercept of the fitted line equation;
  • (step 4) according to a buried depth of the CBM well, determining an abandoned formation pressure P_ab of CBM well exploitation, which comprises determining the abandoned formation pressure by analogy method or empirical formula method, wherein preferably, the abandoned formation pressure is calculated by a Meck empirical formula of P_ab=2.149×10⁻³ D, here, P_ab (MPa) is the abandoned formation pressure, D (m) is the buried depth of the CBM well;
  • (step 5) according to the relationship table between the different pressures and the deviation factors in the (step 1), obtaining deviation factors Z_i and Z_ab by interpolation method, or according to the relationship table between the different pressures and the deviation factors in the (step 1), by taking the deviation factors as a dependent variable and the pressures as an independent variable, establishing fitted function relationship Z=f(P) between the different pressures and the deviation factors, and then by the fitted function relationship, obtaining the deviation factor Z_i=f(P_i) at the original formation pressure P_i, and the deviation factor Z_ab=f(P_ab) at the abandoned formation pressure P_ab; and
  • (step 6) according to the original formation pressure P_i in the (step 1), the adsorbed gas content G_0ad and free gas content G_0f at the original formation pressure and original temperature in the (step 2), the Langmuir's pressure P_L in the (step 3), the abandoned formation pressure P_ab in the (step 4), and the deviation factor Z at the original formation pressure P_i and the deviation factor Z_ab at the abandoned formation pressure P_ab in the (step 5), using a model of

$R=\frac{G_{\mathrm{of}}\times \left(1-\frac{P_{\mathrm{ab}}{Z}_i}{P_i{Z}_{\mathrm{ab}}}\right)+G_{\mathrm{oad}}\times \left(1-\frac{P_{\mathrm{ab}}}{P_i}\times \frac{P_i\times P_L}{P_{\mathrm{ab}}+P_L}\right)}{G_{\mathrm{oad}}+G_{\mathrm{of}}}$

• •  to obtain the recovery factor of the CBM well, wherein:

according to the isothermal adsorption equation (1) of CBM, an adsorbed gas volume V_gi at the original formation pressure P_i is obtained and expressed by an equation of

$\begin{array}{cc}{V}_{\mathrm{gi}}=\frac{V_L\times P_i}{P_i+P_L},& \left(6\right)\end{array}$

• • an adsorbed gas volume V_gab at the abandoned formation pressure P_ab is expressed by an equation of

$V_{\mathrm{gab}}=\frac{V_L\times P_{\mathrm{ab}}}{P_{\mathrm{ab}}+P_L},$

• •  so an absorbed gas cumulative production G_pa at the abandoned formation pressure P_ab is expressed by an equation of

$\begin{array}{cc}{G}_{\mathrm{pa}}=V_{\mathrm{gi}}-V_{\mathrm{gab}}=\frac{V_L\times P_i}{P_i+P_L}-\frac{V_L\times P_{\mathrm{ab}}}{P_{\mathrm{ab}}+P_L},& \left(7\right)\end{array}$

• • an adsorbed gas recovery factor R_a at the abandoned formation pressure P_ab is expressed by an equation of

$\begin{array}{cc}{R}_a=G_{\mathrm{pa}}/V_{\mathrm{gi}},& \left(8\right)\end{array}$

• • it is able to be obtained from the equations (6), (7) and (8) that

$\begin{array}{cc}{R}_a=\frac{\frac{V_L\times P_i}{P_i+P_L}\times \frac{V_L\times P_{\mathrm{ab}}}{P_{\mathrm{ab}}+P_L}}{\frac{V_L\times P_i}{P_i+P_L}}=1-\frac{P_{\mathrm{ab}}}{P_i}\frac{P_i+P_L}{P_{\mathrm{ab}}+P_L},& \left(9\right)\end{array}$

• • here, the equation (9) is a model for calculating the adsorbed gas recovery factor at the abandoned formation pressure P_ab,
  • a material balance equation of free gas is expressed as

$\begin{array}{cc}\frac{P}{Z}=\left(1-\frac{G_P}{G_0}\right)\frac{P_i}{Z_i},& \left(10\right)\end{array}$

• • a free gas recovery factor R_f at the abandoned formation pressure P_ab is expressed as an equation of

$\begin{array}{cc}{R}_f=\frac{G_P}{G_0}=1-\frac{P_{\mathrm{ab}}{Z}_i}{P_i{Z}_{\mathrm{ab}}},& \left(11\right)\end{array}$

• • when the free gas content at the original formation pressure and original temperature is G_0f, and the free gas recovery factor at the abandoned formation pressure P_ab is R_f, then a free gas cumulative production G_pf at the abandoned formation pressure is expressed as an equation of

$\begin{array}{cc}{G}_{\mathrm{pf}}=G_{0f}\times R_f,& \left(12\right)\end{array}$

• • it is able to be obtained from the equations (11) and (12) that

$\begin{array}{cc}{G}_{\mathrm{pf}}=G_{0f}\times \left(1-\frac{P_{\mathrm{ab}}{Z}_i}{P_i{Z}_{\mathrm{ab}}}\right),& \left(13\right)\end{array}$

• • when the adsorbed gas content at the original formation pressure and original temperature is G_0ad, and the adsorbed gas recovery factor at the abandoned formation pressure P_ab is R_a, then the adsorbed gas cumulative production G_pa at the abandoned formation pressure is expressed as an equation of

$\begin{array}{cc}{G}_{\mathrm{pa}}=G_{0\mathrm{ad}}\times R_a,& \left(14\right)\end{array}$

• • it is able to be obtained from the equations (9) and (14) that

$\begin{array}{cc}{G}_{\mathrm{pa}}=G_{0\mathrm{ad}}\times \left(1-\frac{P_{\mathrm{ab}}}{P_i}\frac{P_i+P_L}{P_{\mathrm{ab}}+P_L}\right),& \left(15\right)\end{array}$

• • it is able to be obtained from the equations (13) and (15) that a total cumulative production G_p of the free gas and the adsorbed gas at the abandoned formation pressure P_ab is expressed as

$\begin{array}{cc}{G}_p=G_{\mathrm{pf}}+G_{\mathrm{pa}}=G_{0f}\times \left(1-\frac{P_{\mathrm{ab}}{Z}_i}{P_i{Z}_{\mathrm{ab}}}\right)+G_{0\mathrm{ad}}\times \left(1-\frac{P_{\mathrm{ab}}}{P_i}\frac{P_i+P_L}{P_{\mathrm{ab}}+P_L}\right),& \left(16\right)\end{array}$

• • then, the recovery factor R of CBM at the abandoned formation pressure P_ab is expressed as an equation of

$\begin{array}{cc}R=G_p/\left(G_{0a}+G_{0f}\right)=\frac{G_{0f}\times \left(1\frac{P_{\mathrm{ab}}{Z}_i}{P_i{Z}_{\mathrm{ab}}}\right)+G_{0\mathrm{ad}}\times \left(1-\frac{P_{\mathrm{ab}}{P}_i+P_L}{P_i{P}_{\mathrm{ab}}+P_L}\right)}{G_{0\mathrm{ad}}+G_{0f}},& \left(17\right)\end{array}$

• • wherein the equation (17) is a computation expression of the recovery factor of CBM at the abandoned formation pressure P_ab,
  • here, P_i (MPa) is the original formation pressure, V_g(t) (m³/t) is the adsorbed gas volume at the original formation pressure P_i, P_ab (MPa) is the abandoned formation pressure, P (MPa) is the formation pressure, V_g (m³/t) is the adsorbed gas volume at the formation pressure P, G_pa (m³/t) is the adsorbed gas cumulative production at the abandoned formation pressure, R_a is the adsorbed gas recovery factor at the abandoned formation pressure and is a decimal fraction, G_pf (m³/t) is the free gas cumulative production at the abandoned formation pressure, R_f is the free gas recovery factor at the abandoned formation pressure and is a decimal fraction, z is the gas deviation factor at the formation pressure P and is a decimal fraction, Z_i is the gas deviation factor at the original formation pressure and is a decimal fraction, Z_ab is the gas deviation factor at the abandoned formation pressure P_ab and is a decimal fraction, G_P (m³) is the natural gas cumulative production, G₀ (m³) is the natural gas original reserve, G_0f (m³/t) is the free gas content at the original formation pressure and original temperature, G_0a (m³/t) is the adsorbed gas content at the original formation pressure and original temperature, R is the recovery factor of CBM at the abandoned formation pressure P_ab and is a decimal fraction, P_L (MPa) is Langmuir's pressure, and V_L (m³/t) is Langmuir's volume.

According to the preferred embodiment of the present invention, the original formation pressure and PVT experimental data are collected to establish the relationship table between the pressures and the deviation factors; the isothermal adsorption experimental data, the original adsorbed gas content and the original free gas content at the original formation pressure and original temperature for the deep coal rock are collected; the Langmuir's pressure and the Langmuir's volume are determined by the isothermal adsorption experimental data for the CBM well; the abandoned formation pressure of the CBM well is determined; the derivation factor at the original formation pressure and the derivation factor at the abandoned formation pressure are determined; based on the above original formation pressure, the original adsorbed gas content and the original free gas content at the original formation pressure and original temperature, the Langmuir's pressure, the abandoned formation pressure, the derivation factor at the original formation pressure and the derivation factor at the abandoned formation pressure, the recovery factor of the CBM well is determined. The method provided by the present invention is simple, easy to understand, operable, effective and practical, and has good popularization and application value.

In addition, the present invention is explained with the embodiment as follows, but the protective scope of the present invention is not limited to this embodiment.

First Embodiment

(1) The collected basic data of the CBM well are as follows. The original formation pressure P_i=28 MPa. The relationship between the pressures and the deviation factors of the PVT parameters is shown in Table 1.

TABLE 1
Relationship table of the pressures and the deviation factors
Pressure P	Derivation factor Z	Apparent pressure
(MPa)	(Decimal fraction)	P/Z (MPa)
31	0.9133	33.943
28	0.8805	31.800
25	0.8716	28.683
22	0.8609	25.555
19	0.8631	22.014
16	0.8619	18.564
13	0.8715	14.917
10	0.8770	11.403
8.5	0.8843	9.612
6	0.9047	6.632
3	0.9405	3.190

(2) The collected isothermal adsorption experimental data for deep CBM are shown in Table 2. At the original formation pressure and original temperature, the adsorbed gas content G_0ad is 18.34 m³/t, and the free gas content G_0f is 7.66 m³/t.

TABLE 2
Isothermal adsorption experimental data
and observation data for deep CBM
	P (MPa)	Vg (m3/t)	y(i) = 1/Vg(i)	x(i) = 1/P(i)
	28	18.34	0.0545	0.0357
	26	18.20	0.0549	0.0385
	24	18.04	0.0554	0.0417
	22	17.85	0.0560	0.0455
	20	17.63	0.0567	0.0500
	18	17.37	0.0576	0.0556
	16	17.05	0.0587	0.0625
	14	16.66	0.0600	0.0714
	12	16.17	0.0618	0.0833
	10	15.53	0.0644	0.1000
	8	14.65	0.0683	0.1250
	6	13.38	0.0747	0.1667
	4	11.43	0.0875	0.2500
	2	7.94	0.1259	0.5000

(3) According to adsorbed gas volumes V_g(1), V_g(2), . . . , V_g(n) corresponding to different pressures P₍₁₎, P₍₂₎, . . . , P_(n) at the formation temperature of CBM, when

$y_{\left(i\right)}=\frac{1}{V_{g\left(i\right)}},$

x_(i)=1/P_(i), i=1, 2, . . . , n, a series of observation points (y_(i), x_(i)) (as shown in Table 2) are obtained. Referring to FIG. 2, the linear fitting is performed on the observation points, the Langmuir's volume V_L is reciprocal of the intercept of the line equation obtained by the linear fitting, that is,

$V_L=\frac{1}{0.049}=20.4082$

m³/t, and the Langmuir's pressure P_L is equal to the slope of the fitted line equation divided by the intercept of the fitted line equation, that is,

$P_L=\frac{0.1536}{0.049}=3.1347\mathrm{MPa}.$

(4) The abandoned formation pressure in the CBM well exploitation is determined by analogy method or empirical formula method. In this embodiment, the abandoned formation pressure is calculated by a Meck empirical formula of P_ab=2.149×10⁻³ D, here, P_ab (MPa) is the abandoned formation pressure, D (m) is the buried depth of the CBM well and is about 2850 m, so the abandoned formation pressure is

$P_{\mathrm{ab}}=2.149\times 10^{-3}D=2.149\times {10}^{-3}\times 2850=6.1246\mathrm{MPa}.$

(5) According to the relationship table between the pressures and the deviation factors in the step (1), by linear interpolation method, it is obtained that the gas derivation factor Z_i=0.8805 at the original formation pressure P_i=28 MPa, and the gas derivation factor Z_ab=0.9035 at the abandoned formation pressure P_ab=6.1246 MPa.

(6) According to the original formation pressure P_i=28 MPa in the step (1), the adsorbed gas content G_0ad=18.34 m³/t and the free gas content G_0f=7.66 m³/t at the original formation pressure and original temperature in the step (2), the Langmuir's pressure P_L=3.1347 MPa in the step (3), the abandoned formation pressure P_ab=6.1246 MPa in the step (4), and the deviation factor Z_i=0.8805 at the original formation pressure P_i and the deviation factor Z_ab=0.9035 at the abandoned formation pressure P_ab in the step (5), by a model of

$R=\frac{G_{0f}\times \left(1-\frac{P_{\mathrm{ab}}{Z}_i}{P_i{Z}_{\mathrm{ab}}}\right)+G_{0\mathrm{ad}}\times \left(1-\frac{P_{\mathrm{ab}}}{P_i}\times \frac{P_i+P_L}{P_{\mathrm{ab}}+P_L}\right)}{G_{0\mathrm{ad}}+G_{0f}},$

the recovery factor of the CBM well is obtained, namely,

$\begin{array}{c}R=\frac{G_{0f}\times \left(1-\frac{P_{\mathrm{ab}}{Z}_i}{P_i{Z}_{\mathrm{ab}}}\right)+G_{0\mathrm{ad}}\times \left(1-\frac{P_{\mathrm{ab}}}{P_i}\times \frac{P_i+P_L}{P_{\mathrm{ab}}+P_L}\right)}{G_{0\mathrm{ad}}+G_{0f}}\\ =\frac{\begin{array}{c}7.66\times \left(1-\frac{6.1246\times 0.8805}{28\times 0.9035}\right)+\\ 18.34\times \left(1-\frac{6.1246}{28}\times \frac{28+3.1347}{6.1246+3.1347}\right)\end{array}}{18.34+7.66}\\ =0.4183.\end{array}$

Therefore, according to the first embodiment of the present invention, the recovery factor of the CBM well is 41.83%.

The above embodiment is the preferred embodiment of the present invention, but not the limitation to other embodiments of the present invention. Any other changes that do not deviate from the present invention shall be equivalent replacement modes and shall be included in the protective scope of the present invention.

Claims

1. A method for determining a recovery factor of a deep coal bed methane (CBM) well, the method comprising steps of: (step 1) collecting an original formation pressure Pi and pressure-volume-temperature (PVT) experimental data of the CBM well, obtaining deviation factors corresponding to different pressures at a formation temperature of CBM, and establishing a relationship table between the pressures and the deviation factors;(step 2) obtaining isothermal adsorption experimental data, an adsorbed gas content G0ad and a free gas content G0f at the original formation pressure and an original temperature for deep coal rock;(step 3) according to the isothermal adsorption experimental data, determining Langmuir's pressure PL and Langmuir's volume VL;(step 4) according to a depth of burial of the CBM well, determining an abandoned formation pressure Pab in CBM well exploitation;(step 5) according to the relationship table between the pressures and the deviation factors in the (step 1), obtaining a deviation factor Zi at the original formation pressure Pi and a deviation factor Zab at the abandoned formation pressure Pab; and(step 6) according to the original formation pressure Pi in the (step 1), the adsorbed gas content G0ad and the free gas content G0f at the original formation pressure and the original temperature in the (step 2), the Langmuir's pressure PL in the (step 3), the abandoned formation pressure Pab in the (step 4), and the deviation factor Zi at the original formation pressure Pi and the deviation factor Zab at the abandoned formation pressure Pab in the (step 5), by a model of

2. The method according to claim 1, wherein the (step 3) specifically comprises:

3. The method according to claim 1, wherein in the (step 4), the abandoned formation pressure Pab in CBM well exploitation is determined by analog method or empirical formula method.

4. The method according to claim 1, wherein in the (step 5), according to the relationship table between the pressures and the deviation factors in the (step 1), the deviation factor Zi at the original formation pressure Pi and the deviation factor Zab at the abandoned formation pressure Pab are obtained by interpolation method; or according to the relationship table between the pressures and the deviation factors in the (step 1), a fitted function relationship Z=f(P) between the deviation factors and the pressures is established by taking the deviation factors as a dependent variable and the pressures as an independent variable, and then by the fitted function relationship, obtaining the deviation factor Zi=f(Pi) at the original formation pressure Pi, and the deviation factor Zab=f(Pab) at the abandoned formation pressure Pab.