Method for determining recovery factor of deep coal bed methane wells
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.
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 InventionThe 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 ArtsIn 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 INVENTIONAn 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 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 the 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), 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
-
- calculating the recovery factor of the CBM well.
Preferably, the (step 3) specifically comprises:
-
- according to adsorbed gas volumes Vg(1), Vg(2), . . . , Vg(n) corresponding to different pressures P(1), P(2), . . . , 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 VL is equal to a reciprocal of an intercept of the line equation, and the Langmuir's pressure PL 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 Pab 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 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.
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.
Referring to
(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 different 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 the deep coal rock;
(step 3) according to the isothermal adsorption experimental data, determining Langmuir's pressure PL and Langmuir's volume VL, wherein:
-
- an isothermal adsorption equation of CBM is
-
- by taking the reciprocal at both sides of the equation (1), an equation (2) is obtained as follows
-
- and then based on the equation (2), an equation (3) is obtained as follows
-
- here, i is isothermal adsorption experimental data point number, i=1, 2, . . . , n, Vg(i) (m3/t) is adsorbed gas volume at the ith experimental data point, P(i) (MPa) is pressure at the ith experimental data point, VL (m3/t) is Langmuir's volume, PL (MPa) is Langmuir's pressure,
- according to adsorbed gas volumes Vg(1), Vg(2), . . . , Vg(n) corresponding to different pressures P(1), P(2), . . . , 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 VL is reciprocal of an intercept of a line equation obtained by linear fitting, and the Langmuir's pressure PL 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 Pab 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 Pab=2.149×10−3 D, here, Pab (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 Zi and Zab 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 Zi=f(Pi) at the original formation pressure Pi, and the deviation factor Zab=f(Pab) 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 free gas content G0f at the original formation pressure and 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 Z at the original formation pressure Pi and the deviation factor Zab at the abandoned formation pressure Pab in the (step 5), using a model of
-
- to obtain the recovery factor of the CBM well, wherein:
according to the isothermal adsorption equation (1) of CBM, an adsorbed gas volume Vgi at the original formation pressure Pi is obtained and expressed by an equation of
-
- an adsorbed gas volume Vgab at the abandoned formation pressure Pab is expressed by an equation of
-
- so an absorbed gas cumulative production Gpa at the abandoned formation pressure Pab is expressed by an equation of
-
- an adsorbed gas recovery factor Ra at the abandoned formation pressure Pab is expressed by an equation of
-
- it is able to be obtained from the equations (6), (7) and (8) that
-
- here, the equation (9) is a model for calculating the adsorbed gas recovery factor at the abandoned formation pressure Pab,
- a material balance equation of free gas is expressed as
-
- a free gas recovery factor Rf at the abandoned formation pressure Pab is expressed as an equation of
-
- when the free gas content at the original formation pressure and original temperature is G0f, and the free gas recovery factor at the abandoned formation pressure Pab is Rf, then a free gas cumulative production Gpf at the abandoned formation pressure is expressed as an equation of
-
- it is able to be obtained from the equations (11) and (12) that
-
- when the adsorbed gas content at the original formation pressure and original temperature is G0ad, and the adsorbed gas recovery factor at the abandoned formation pressure Pab is Ra, then the adsorbed gas cumulative production Gpa at the abandoned formation pressure is expressed as an equation of
-
- it is able to be obtained from the equations (9) and (14) that
-
- it is able to be obtained from the equations (13) and (15) that a total cumulative production Gp of the free gas and the adsorbed gas at the abandoned formation pressure Pab is expressed as
-
- then, the recovery factor R of CBM at the abandoned formation pressure Pab is expressed as an equation of
-
- wherein the equation (17) is a computation expression of the recovery factor of CBM at the abandoned formation pressure Pab,
- here, Pi (MPa) is the original formation pressure, Vg(t) (m3/t) is the adsorbed gas volume at the original formation pressure Pi, Pab (MPa) is the abandoned formation pressure, P (MPa) is the formation pressure, Vg (m3/t) is the adsorbed gas volume at the formation pressure P, Gpa (m3/t) is the adsorbed gas cumulative production at the abandoned formation pressure, Ra is the adsorbed gas recovery factor at the abandoned formation pressure and is a decimal fraction, Gpf (m3/t) is the free gas cumulative production at the abandoned formation pressure, Rf 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, Zi is the gas deviation factor at the original formation pressure and is a decimal fraction, Zab is the gas deviation factor at the abandoned formation pressure Pab and is a decimal fraction, GP (m3) is the natural gas cumulative production, G0 (m3) is the natural gas original reserve, G0f (m3/t) is the free gas content at the original formation pressure and original temperature, G0a (m3/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 Pab and is a decimal fraction, PL (MPa) is Langmuir's pressure, and VL (m3/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 Pi=28 MPa. The relationship between the pressures and the deviation factors of the PVT parameters is shown in Table 1.
(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 G0ad is 18.34 m3/t, and the free gas content G0f is 7.66 m3/t.
(3) According to adsorbed gas volumes Vg(1), Vg(2), . . . , Vg(n) corresponding to different pressures P(1), P(2), . . . , P(n) at the formation temperature of CBM, when
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
m3/t, and the Langmuir's pressure PL is equal to the slope of the fitted line equation divided by the intercept of the fitted line equation, that is,
(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 Pab=2.149×10−3 D, here, Pab (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
(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 Zi=0.8805 at the original formation pressure Pi=28 MPa, and the gas derivation factor Zab=0.9035 at the abandoned formation pressure Pab=6.1246 MPa.
(6) According to the original formation pressure Pi=28 MPa in the step (1), the adsorbed gas content G0ad=18.34 m3/t and the free gas content G0f=7.66 m3/t at the original formation pressure and original temperature in the step (2), the Langmuir's pressure PL=3.1347 MPa in the step (3), the abandoned formation pressure Pab=6.1246 MPa in the step (4), and the deviation factor Zi=0.8805 at the original formation pressure Pi and the deviation factor Zab=0.9035 at the abandoned formation pressure Pab in the step (5), by a model of
the recovery factor of the CBM well is obtained, namely,
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: R = G 0 f × ( 1 - P ab Z i P i Z ab ) + G 0 ad × ( 1 - P ab P i × P i + P L P ab + P L ) G 0 ad + G 0 f,
- (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
- calculating the recovery factor of the CBM well.
2. The method according to claim 1, wherein the (step 3) specifically comprises: taking y ( i ) = 1 V g ( i ) and x ( i ) = 1 / P ( i ),, wherein i = 1, 2, …, n;
- according to adsorbed gas volumes Vg(1), Vg(2),..., Vg(n) corresponding to different pressures P(1), P(2),..., 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 VL is equal to a reciprocal of an intercept of the line equation, and the Langmuir's pressure PL is equal to a slope of the line equation divided by the intercept of the line equation.
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.
Type: Application
Filed: Sep 3, 2024
Publication Date: Dec 26, 2024
Inventors: Bo Hu (Zhengzhou city), Yongyi Zhou (Zhengzhou city), Kui Chen (Zhengzhou city), Linsong Liu (Zhengzhou city), Yaonan Yu (Zhengzhou city), Yanling Luo (Zhengzhou city), Lu Cui (Zhengzhou city)
Application Number: 18/823,654