A method for determining recovery factor of a layered gas reservoir under different well interference conditions

Abstract

The application discloses a kind of prediction methods for determining recovery ratio of stratified gas reservoir under different well interference conditions, and belongs to the technical field of oil and gas field development.It includes collecting and arranging gas reservoir basic data;Establish well control area under different well interference, and, based on the existing production pressure difference condition gas reservoir waste gas pressure determination formula;Finally, consider the determination of stratified gas reservoir recovery ratio of well control area.The present application determines the accurate gas reservoir abandonment condition by definite integral principle under the maximum production pressure difference, overcomes the problem that gas reservoir abandonment pressure is difficult to obtain;And based on the existing well pattern condition, establish the determination method of gas reservoir recovery ratio, achieve the purpose of quickly and accurately obtaining gas reservoir recovery ratio.

Description

A method for predicting the recovery rate of layered gas reservoirs under different inter-well disturbance conditions. Technical Field

[0001] This invention relates to a method for predicting recovery rate, and more particularly to a method for predicting recovery rate of layered gas reservoirs under different inter-well interference conditions, belonging to the field of oil and gas field development technology. Background Technology

[0002] Carbonate layered gas reservoirs have complex reservoir spaces and diverse types. They are classified into fractured gas reservoirs, fracture-pore gas reservoirs, and porous gas reservoirs, with fracture-pore and porous gas reservoirs accounting for over 80% and contributing significantly to production. However, the control range of different gas wells varies depending on the mutual interference between wells, resulting in different combinations. This inevitably presents technical challenges in determining the actual well-controlled area and formation pressure.

[0003] Recovery rate (RRR) refers to the ratio of recovered natural gas to the original reserves of a gas reservoir, and it is an important indicator for evaluating the effectiveness of gas reservoir development. RRR evaluation research is also a crucial aspect of natural gas reservoir development. Currently, the main methods for predicting RRR are the empirical method and the material balance method (abandoned pressure method), used to evaluate the rationality of current methods and the effectiveness of gas reservoir development. The empirical method selects a corresponding RRR value range based on the "Methods for Estimating Recoverable Natural Gas Reserves" (SY / T 6098-2010), as shown in Figure 1; the material balance method uses formula E... R =1 - (Pa / Za) / (Pi / Zi), where: ER is the recovery rate, in MPa; Pa / Za is the apparent abandonment pressure; Pi / Zi is the original apparent formation pressure, in MPa; Pa / Za is the abandonment pressure condition of the gas reservoir, which can only be measured by downhole tools using conventional material balance methods, thus allowing the calculation of the recovery rate E. R However, in this prediction process, the recovery rate of the gas reservoir is mainly related to the abandoned pressure of the gas reservoir and the reserves within the well network control range, which presents the following technical problems:

[0004] I. The material balance method requires high accuracy of pressure data, with at least two static pressure test points at different times;

[0005] Second, the material balance method cannot determine the uncontrolled reserves of a gas reservoir, which is not conducive to expanding the controlled reserves range with subsequent supplementary wells.

[0006] Third, the gas reservoir abandonment pressure in the material balance method is determined by empirical methods, which has certain data errors and therefore cannot guarantee the accuracy of the recovery rate.

[0007] Furthermore, the existing technology CN1464429, a method for calculating gas reservoir recovery rate and recoverable reserves, aims to consider the dynamic reserve ratio and the difference in gas saturation between the original and abandoned stages of water-driven gas reservoirs when calculating gas reservoir recovery rate and recoverable reserves. It proposes a new method for calculating gas reservoir recovery rate and recoverable reserves—the "corrected volumetric method"—deriving a new calculation formula and improving the theory of the original method. CN114021821A, a gas reservoir recovery rate prediction method based on multiple regression, measures gas reservoir recovery rate by simulating actual reservoir gas extraction methods and reconstructing the formation seepage process. It is suitable for measuring the recovery rate of various types of gas reservoirs and gas reservoirs with different development methods and stages. In other words, the aforementioned existing technologies do not solve the problem of determining the recovery rate of layered gas reservoirs under different inter-well interference conditions.

[0008] Therefore, it is necessary to improve or optimize the current methods for predicting the recovery rate of layered gas reservoirs, accurately evaluate the production and development effects of gas reservoirs, improve the economic benefits of gas field development, and at the same time play an important guiding role in improving the utilization of gas reservoir reserves and increasing the recovery rate. Summary of the Invention

[0009] This invention aims to overcome the shortcomings of existing technologies by proposing a predictive method for determining the recovery rate of layered gas reservoirs under different inter-well interference conditions. This method solves the problem of existing material balance methods requiring a large number of pressure measurement points at different times. In this technical solution, the accurate gas reservoir abandonment conditions are solved using the definite integral principle under the maximum production pressure differential, overcoming the difficulty in obtaining the gas reservoir abandonment pressure. Furthermore, a method for determining the gas reservoir recovery rate is established based on existing well network conditions, achieving the goal of obtaining the gas reservoir recovery rate quickly and accurately.

[0010] To achieve the above technical objectives, the following technical solution is proposed:

[0011] The primary objective of this technical solution is to propose a method for predicting the recovery rate of layered gas reservoirs under different inter-well interference conditions, comprising the following steps:

[0012] S1: Based on the basic data of the gas reservoir, determine the well-controlled area under the existing well network conditions, and establish the first relationship about the well-controlled area based on the natural gas seepage theory;

[0013] S2: Based on the basic data of the gas reservoir, determine the abandoned pressure of the gas reservoir under the existing production pressure differential conditions, and establish a second relationship regarding the abandoned pressure of the gas reservoir;

[0014] S3: Based on the first and second relations, establish the third relation on the recovery rate to determine the recovery rate under the existing well network conditions and the existing production pressure differential conditions.

[0015] In step S1, the first relationship regarding the well-controlled area is established for three different cases, including:

[0016] S11: The control ranges of the two gas wells do not intersect, i.e., r1 + r2 <d;

[0017] That is, the control range of the gas reservoir is the sum of the control areas of the two gas wells, and the resulting well control area is shown in the following formula (1);

[0018] (1)

[0019] S12: The control ranges of the two gas wells completely overlap, i.e., r2-r1≥d;

[0020] That is, the control range of the gas reservoir is the control area S2 of the gas well with a larger radius, and the resulting well control area is shown in the following formula (2);

[0021] (2)

[0022] S13: The control ranges of the two gas wells intersect, i.e., r2-r1 <d<r1+r2;

[0023] The control area of ​​a gas reservoir cannot be simply defined as S1 + S2, as this would fall into the trap of the single-well aggregation method and lead to an overestimation of the calculated reserves. To address this, this invention considers the superposition of pressure fields, i.e., solving for the area of ​​the boundary figure. The control area is equal to the sum of the areas of the two gas wells on either side. The area of ​​a single gas well can be considered as two symmetrical parts along the X-axis. According to the principles of calculus, the control areas of the two gas wells are respectively:

[0024] ;

[0025] The resulting well control area is shown in formula (3);

[0026] (3)

[0027] In equations (1), (2), and (3), r1 and r2 are the well control radii of the two gas wells, in meters; d is the distance between the two gas wells, in meters; S1 and S2 are the control areas of the two gas wells, and S is the control area of ​​the gas reservoir, both in meters. 2 ;x m Let be any point on the control boundary of the gas reservoir, in meters (m).

[0028] In step S2, a second relationship regarding the gas reservoir abandonment pressure is established for three different cases, including:

[0029] S21: The control ranges of the two gas wells do not intersect, i.e., r1 + r2 <d;

[0030] Based on the principle of calculus, the control area of ​​the two gas wells is processed into a circular slice, and the gas reservoir abandonment pressure of the two gas wells is shown in the following formulas (4) and (5), respectively.

[0031] (4)

[0032] (5)

[0033] S22: The control ranges of the two gas wells completely overlap, i.e., r2-r1≥d;

[0034] The abandoned pressure of the gas reservoir in the controlled area is shown in equation (6);

[0035] (6)

[0036] S23: The control ranges of the two gas wells intersect, i.e., r2-r1 <d<r1+r2;

[0037] The abandoned pressure of the gas reservoir in the controlled area is shown in equation (7);

[0038] (7)

[0039] In equations (4), (5), (6), and (7), P1 and P2 are the reservoir abandonment pressures of the two gas wells, respectively, and P is the reservoir abandonment pressure of the controlled area. These are the limiting bottomhole flowing pressures of the two gas wells, P. e This refers to the original formation pressure of the gas reservoir, and the unit is MPa. R is the wellbore radius, and R is the gas reservoir radius; both are in meters (m).

[0040] In step S3, based on the principle of mass balance, the expression for the gas reservoir recovery rate is as follows:

[0041] ;

[0042] in, The recovery rate of the integrated gas reservoir is expressed as % (%). These represent the abandonment pressure and corresponding deviation factor of the complete gas reservoir, respectively. These are the original formation pressure of the gas reservoir and the original natural gas deviation factor, respectively, both dimensionless.

[0043] Furthermore, when considering the uncontrolled reserves of gas reservoirs under existing well network conditions, the expression for gas reservoir recovery rate is improved as follows:

[0044] (8)

[0045] In equation (8), S R The gas reservoir area is expressed in km². 2 E represents the recovery rate of the controlled portion of the gas well, in units of %.

[0046] Then, establish a third relationship regarding the recovery rate for three different scenarios, including:

[0047] S31: The control ranges of the two gas wells do not intersect, i.e., r1 + r2 <d;

[0048] The gas reservoir has two control areas, and the recovery rates of the two gas wells are obtained respectively. Then, the area weighting is used to obtain the recovery rate of the entire gas reservoir, as shown in the following formula (9);

[0049] (9)

[0050] S32: The control ranges of the two gas wells completely overlap, i.e., r2-r1≥d;

[0051] The gas reservoir has only one control range. The gas reservoir recovery rate under the current conditions is obtained as shown in the following formula (10);

[0052] (10)

[0053] Z3: The control ranges of the two gas wells intersect, i.e., r2-r1 <d<r1+r2;

[0054] The gas reservoir is a control range, and the recovery rate of the gas reservoir under the current conditions is obtained as shown in the following formula (11);

[0055] (11)

[0056] In equations (9), (10) and (11), Z is the deviation factor of the complete gas reservoir, and Z1 and Z2 are the natural gas deviation factors of the two gas wells, respectively, which are dimensionless.

[0057] The second objective of this technical solution is to provide: a computer-readable storage medium storing a computer program, wherein when the computer program is executed by a processor, it implements the steps of the above-mentioned method for predicting the recovery rate of layered gas reservoirs under different inter-well interference conditions.

[0058] The third objective of this technical solution is to provide an information data processing terminal for a method of predicting the recovery rate of layered gas reservoirs under different inter-well interference conditions.

[0059] The beneficial technical effects of adopting this technical solution are as follows:

[0060] 1. This invention solves the technical problem of the material balance method requiring pressure test data points at different times, and the formation pressure under different inter-well interference conditions can be determined by the determination method provided in this invention.

[0061] 2. This invention provides three methods for determining the controlled area of ​​gas reservoirs under different inter-well interference conditions, accurately determining the range of uncontrolled reserves, which can be used to guide subsequent supplementary wells to expand the control range of reserves. This is a major improvement on the material balance method.

[0062] 3. By using the definite integral principle under the maximum production pressure difference, the present invention solves the accurate gas reservoir abandonment condition, overcoming the problem of difficultly obtaining the gas reservoir abandonment pressure.

[0063] 4. By calculating the gas reservoir control area and abandonment pressure under different interference conditions, the present invention establishes a set of prediction methods for the recovery factor of layered gas reservoirs, which can be used for the recovery factor prediction of different complex well patterns. It is more systematic in method and more scientific and accurate in data solution than the traditional gas reservoir material balance method, providing a theoretical basis for quickly and accurately obtaining the gas reservoir recovery factor. BRIEF DESCRIPTION OF THE DRAWINGS

[0064] FIG. 1 shows the recovery factor ranges of different types of gas reservoirs in the prior art.

[0065] FIG. 2 is a flow block diagram related to the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0066] The technical solutions in the embodiments of the present invention will be clearly and completely described below. Apparently, the described embodiments are only a part of the embodiments of the present invention, rather than all of the embodiments. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0067] Embodiment 1

[0068] When the control ranges of two gas wells involved in this embodiment do not intersect, that is, r1 + r2 < d; a prediction method for determining the recovery factor under different well interference conditions in a layered gas reservoir is provided. As shown in FIG. 2, it specifically includes:

[0069] [[ID=??]]

[0070] [[ID=2??]] (1)

[0071] In the formula, r1 and r2 are the well control radii of the two gas wells, with the unit of m; d is the distance between the two gas wells, with the unit of m; S1 and S2 are the control ranges (i.e., control areas) of the two gas wells, and S is the control area of the gas reservoir, all with the unit of m 2 .

[0072] B. Determine the abandonment pressure of the gas reservoir under the existing production pressure difference, and establish a second relational expression for the abandonment pressure of the gas reservoir. Among them, according to the calculus principle, the control areas of the two gas wells are processed by circular ring slicing, and the abandonment pressures of the two gas wells are respectively as follows: It should be noted that there seems to be some unclear or incomplete content in the original text, such as the specific expressions in steps A and B that are not fully presented. The translation is based on the existing text as accurately as possible.

[0073]

[0074]

[0075] In the formula, P1 and P2 are the reservoir abandonment pressures of the two gas wells, respectively. These are the limiting bottomhole flowing pressures of the two gas wells, P. e This refers to the original formation pressure of the gas reservoir, and the unit is MPa. R is the wellbore radius, and R is the gas reservoir radius; both are in meters (m).

[0076] C. Based on the first and second relations, establish the third relation regarding the recovery rate to determine the recovery rate under the existing well network conditions and the existing production pressure differential conditions; wherein, based on the principle of material balance, the expression for the gas reservoir's recovery rate is as follows:

[0077] ;

[0078] In the formula, The recovery rate of the integrated gas reservoir is expressed as % (%). These represent the abandonment pressure and corresponding deviation factor of the complete gas reservoir, respectively. These are the original formation pressure of the gas reservoir and the original natural gas deviation factor, respectively, both dimensionless.

[0079] Since the gas reservoir has two control areas, the recovery rates of the two gas wells are obtained separately. Then, the area weighting is used to obtain the recovery rate of the entire gas reservoir, as shown in the following formula.

[0080] surface

[0081] In the formula, These represent the recovery rates for the two gas wells, respectively.

[0082] Example 2

[0083] This embodiment involves two gas wells whose control ranges completely overlap, i.e., r2-r1≥d; it provides a prediction method for determining the recovery rate of a layered gas reservoir under different inter-well interference conditions, as shown in Figure 2, specifically including:

[0084] A. Determine the well control area under the existing well network conditions. Based on the natural gas seepage theory, establish the first relationship regarding the well control area. Among them, the control range of the gas reservoir is the control area S2 of the gas well with the larger radius. The obtained well control area is shown in the following formula.

[0085]

[0086] Let \(r_1\) and \(r_2\) be the well control radii of two gas wells, with the unit of m; \(d\) be the distance between the two gas wells, with the unit of m; \(S_2\) be the control area of the gas well with a larger radius, and \(S\) be the control range (i.e., control area) of the gas reservoir, both with the unit of \(m^2\). 2 ;

[0087] B. Determine the abandonment pressure of the gas reservoir under the existing production pressure difference condition, and establish a second relationship regarding the abandonment pressure of the gas reservoir; among them, the abandonment pressure of the gas reservoir for the control area is shown in the following formula;

[0088]

[0089] In the formula, \(P_1\) and \(P_2\) are the abandonment pressures of the gas reservoirs of the two gas wells respectively, \(P\) is the abandonment pressure of the gas reservoir for the control area, and the unit of all is MPa; \(S_1\) is the control area of the gas well with a smaller radius, with the unit of \(m^2\). 2 ;

[0090] C. Establish a third relationship regarding the recovery factor based on the first relationship and the second relationship, and determine the recovery factor under the existing well pattern condition and the existing production pressure difference condition; among them, according to the principle of material balance, the expression of the recovery factor of the gas reservoir is as follows:

[0091] ;

[0092] In the formula, is the recovery factor of the integral gas reservoir, with the unit of %; are the abandonment pressure and the corresponding deviation factor of the integral gas reservoir respectively, are the original formation pressure of the gas reservoir and the original gas deviation factor respectively, dimensionless.

[0093] Since there is only one control range for the gas reservoir, the recovery factor of the gas reservoir under the current condition is obtained as shown in the following formula;

[0094]

[0095] In the formula, \(R\) refers to the radius of the gas reservoir, with the unit of m; \(Z\) is the deviation factor of the integral gas reservoir, dimensionless.

[0096] Example 3

[0097] For the two gas wells involved in this example, their control ranges intersect, that is, \(r_2 - r_1 < d < r_1 + r_2\); provide a prediction method for determining the recovery factor under different well interference conditions in a layered gas reservoir, as shown in Figure 2, specifically including:

[0098] A. Determine the well control area under the existing well network conditions. Based on the natural gas seepage theory, establish the first relationship regarding the well control area. However, the control range of the gas reservoir cannot be simply S1 + S2, as this would fall into the trap of the single-well aggregation method and lead to an overestimation of the calculated reserves. For this situation, this embodiment considers the superposition of pressure fields, i.e., solving for the area of ​​the boundary figure. The control area is equal to the sum of the areas of the two wells on the left and right sides. The area of ​​a single gas well can be considered as two symmetrical parts along the X-axis. According to the principles of calculus, the control areas of the two gas wells are respectively:

[0099] ;

[0100] The resulting well control area is shown in the following formula;

[0101]

[0102] In the formula, r1 and r2 are the well control radii of the two gas wells, respectively, in meters; d is the distance between the two gas wells, in meters; S1 and S2 are the control areas of the two gas wells, respectively, and S is the control range (i.e., control area) of the gas reservoir, both in meters. 2 ;x m Let be any point on the control boundary of the gas reservoir, in meters (m).

[0103] B. Determine the reservoir abandonment pressure under the existing production pressure differential conditions, and establish a second relationship regarding the reservoir abandonment pressure; wherein, the reservoir abandonment pressure of the two gas wells is shown in the following formula;

[0104]

[0105] In the formula, P1 and P2 are the reservoir abandonment pressures of the two gas wells, respectively, and P is the reservoir abandonment pressure of the controlled area, with the unit being MPa.

[0106] Z. Based on the first and second relations, establish the third relation regarding the recovery rate to determine the recovery rate under the existing well network conditions and the existing production pressure differential conditions; wherein, based on the principle of material balance, the expression for the gas reservoir's recovery rate is as follows:

[0107] ;

[0108] In the formula, The recovery rate of the integrated gas reservoir is expressed as % (%). These represent the abandonment pressure and corresponding deviation factor of the complete gas reservoir, respectively. These are the original formation pressure of the gas reservoir and the original natural gas deviation factor, respectively, both dimensionless.

[0109] Since the gas reservoir is a controlled area, the recovery rate of the gas reservoir under the current conditions is obtained as shown in the following formula;

[0110]

[0111] In the formula, Z is the deviation factor of the complete gas reservoir, which is dimensionless; R refers to the gas reservoir radius, and the unit is m.

[0112] Example 4

[0113] Based on Examples 1-3, this example improves the gas recovery rate expression to consider the uncontrolled reserve range of gas reservoirs under existing well network conditions:

[0114] ;

[0115] In the formula, S R The gas reservoir area is expressed in km². 2 E represents the recovery rate of the controlled portion of the gas well, in units of %.

[0116] Example 5

[0117] Based on Example 1, this example uses two pure gas wells (Well 1 and Well 2) in a gas reservoir in the Sichuan Basin to further illustrate the invention.

[0118] The gas wells are all 5400m deep, with ultimate bottomhole flowing pressures of 6.7MPa and 8MPa, respectively. The corresponding control radius is 670m and 820m, respectively, and the controlled reserves are 569 million cubic meters and 1.709 billion cubic meters, respectively. The geological reserves of the block are 4 billion cubic meters, and the area is 9.8km². 2 The original formation pressure was 59 MPa.

[0119] 1. Determine the well control area under the existing well network conditions: The total control area of ​​the two gas wells is the sum of the areas of the two gas wells, i.e., 3.14 (670). 2 +820 2 = 3.52km 2 ;

[0120] 2. Determine the reservoir abandonment pressure under the existing production pressure differential: Based on the principle of calculus, the control area of ​​the two gas wells is processed into a circular slice, and the reservoir abandonment pressures of the two gas wells are obtained as 10MPa and 12MPa, respectively.

[0121] 3. The recovery rate under the current well network and production pressure differential conditions is obtained: Based on the recovery rates of 85.7% and 81.2% within the control range of the two gas wells, respectively, the area-weighted recovery rate of the entire gas reservoir is then obtained as 28.8%.

[0122] Example 6

[0123] Based on Example 2, this example uses two pure gas wells (Well 3 and Well 4) in a gas reservoir in the Sichuan Basin to further illustrate the present invention.

[0124] The gas wells are all 5400m deep, with ultimate bottomhole flowing pressures of 7MPa and 8.4MPa, respectively. The corresponding control radius is 1020m and 2050m, respectively, and the controlled reserves are 1.713 billion cubic meters and 4.309 billion cubic meters, respectively. The geological reserves of the block are 8 billion cubic meters, and the area is 23.2 km². 2 The original formation pressure was 59 MPa.

[0125] 1. Determine the well control area under the existing well network conditions: The combined control area of ​​the two gas wells is the area of ​​Well 4, i.e., 3.14 × 2050. 2 =13.19km 2 ;

[0126] 2. Determine the reservoir abandonment pressure under the existing production pressure differential: Considering the superposition effect of inter-well production pressure differential, the reservoir abandonment pressure of both gas wells is 11 MPa.

[0127] 3. The recovery rate under the current well network and production pressure differential conditions is obtained: Based on the fact that only one gas well has a recovery rate of 86.2% within its control range, the area-weighted recovery rate of the entire gas reservoir is then obtained as 48.3%.

[0128] Example 7

[0129] Based on Example 3, this example uses two pure gas wells (Well 5 and Well 6) in a gas reservoir in the Sichuan Basin to further illustrate the present invention.

[0130] The gas wells are all 5400m deep, with ultimate bottomhole flowing pressures of 6.6MPa and 7.4MPa, respectively. The corresponding control radius is 1100m and 1168m, respectively, and the controlled reserves are 1.832 billion cubic meters and 2.814 billion cubic meters, respectively. The geological reserves of the block are 75 billion cubic meters, and the area is 28.3 km². 2 The original formation pressure was 59 MPa.

[0131] 1. Determine the well control area under the existing well network conditions: The control range is equal to the sum of the areas of the two gas wells on the left and right sides. The area of ​​a single gas well can be considered as two symmetrical parts along the x-axis. According to the principle of calculus, the control area of ​​the two gas wells is 16.08 km². 2 ;

[0132] 2. Determine the reservoir abandonment pressure under the existing production pressure differential: Considering the superposition effect and the area weighting principle, the reservoir abandonment pressure of both gas wells is 12.8 MPa;

[0133] 3. The recovery rate under the current well network and production pressure differential conditions is obtained: Based on the fact that the recovery rate of one of the two gas wells within the control range is 84.8%, the area weighted average yields a recovery rate of 47.5% for the entire gas reservoir.

Claims

1. A method for predicting the recovery rate of a layered gas reservoir under different inter-well disturbance conditions, characterized in that, It includes the following steps: S1: According to the basic data of the gas reservoir, determine the well control area under the existing well pattern conditions, and establish a first relationship about the well control area; among them, the first relationship about the well control area is established in three cases respectively, including: S11: The control ranges of two gas wells do not intersect, that is, r1 + r2 < d; the obtained well control area is shown in the following formula (1). S12: The control ranges of the two gas wells completely overlap, i.e., r2-r1≥d; the resulting well control area is shown in the following formula (2); S13: The control ranges of two gas wells intersect, i.e., r2 - r1 < d < r1 + r2; the control areas of the two gas wells are respectively: The resulting well control area is shown in formula (3) below. In Formulas (1), (2), and (3), r1 and r2 are the well control radii of two gas wells, with the unit of m; d is the distance between the two gas wells, with the unit of m; S1 and S2 are the controlled areas of the two gas wells, and S is the controlled area of the gas reservoir, all with the unit of m 2 ; x m is an arbitrary point on the control boundary of the gas reservoir, with the unit of m; S2: Based on the basic data of the gas reservoir, determine the abandonment pressure of the gas reservoir under the existing production pressure difference conditions, and establish a second relationship about the abandonment pressure of the gas reservoir; among them, the second relationship about the abandonment pressure of the gas reservoir is established in three cases, including: S21: The control ranges of the two gas wells do not intersect, that is, r1 + r2 < d; the abandonment pressures of the two gas wells are shown in the following Formulas (4) and (5) respectively; S22: The control ranges of the two gas wells completely overlap, i.e., r2-r1≥d; the gas reservoir abandonment pressure of the control area is shown in the following formula (6); S23: The controlled areas of two gas wells intersect, i.e., r2 - r1 < d < r1 + r2; the abandonment pressure of the gas reservoir in the controlled area is shown in the following formula (7); In equations (4), (5), (6), and (7), P1 and P2 are the reservoir abandonment pressures of the two gas wells, respectively, and P is the reservoir abandonment pressure of the controlled area. These are the limiting bottomhole flowing pressures of the two gas wells, P. e This refers to the original formation pressure of the gas reservoir, and the unit is MPa. Where R is the wellbore radius and the gas reservoir radius is the gas reservoir radius; both are in meters. S3: Based on the first and second relationships, establish the third relationship regarding the recovery rate to determine the recovery rate under existing well network conditions and existing production pressure differential conditions. The expression for the gas reservoir recovery rate, based on the principle of material balance, is as follows: ;in, The recovery rate of the integrated gas reservoir is expressed as % (%). These represent the abandonment pressure and corresponding deviation factor of the complete gas reservoir, P. e Z e These are the original formation pressure of the gas reservoir and the original natural gas deviation factor, respectively, both dimensionless.

2. The method for predicting the recovery rate of a layered gas reservoir under different inter-well interference conditions according to claim 1, characterized in that, In step S1, when establishing the first relationship about the well control area, it includes based on the natural gas seepage theory.

3. The method for predicting the recovery rate of a layered gas reservoir under different inter-well interference conditions according to claim 1, characterized in that, In step S3, when there is an uncontrolled reserve range in the gas reservoir under the existing well pattern conditions, the recovery factor expression of the gas reservoir is as follows: In equation (8), S R The gas reservoir area is expressed in km². 2 E represents the recovery rate of the controlled portion of the gas well, in units of %.

4. The method for determining the recovery rate of an integrated gas reservoir considering well-controlled area according to claim 1, characterized in that, The third relationship about the recovery factor is established in three cases respectively, including: S31: The control ranges of two gas wells do not intersect, that is, r1 + r2 < d; there are two control ranges in the gas reservoir, and the recovery factors of the two gas wells are obtained respectively, and then the recovery factor of the integrated gas reservoir is obtained, as shown in the following formula (9). S32: The control ranges of the two gas wells completely overlap, i.e., r2-r1≥d; the gas reservoir has only one control range, and the recovery rate of the gas reservoir under the current conditions is obtained as shown in the following formula (10); S33: The controlled areas of two gas wells intersect, i.e., r2 - r1 < d < r1 + r2; the gas reservoir is one controlled area, and the recovery factor of the gas reservoir under the current conditions is obtained as shown in the following formula (11); In equations (9), (10) and (11), Z is the deviation factor of the complete gas reservoir, and Z1 and Z2 are the natural gas deviation factors of the two gas wells, respectively, which are dimensionless.

5. A computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, it implements the prediction method for determining the recovery factor under different well interference conditions in a layered gas reservoir according to any one of claims 1-4.

6. An information data processing terminal for the prediction method of determining the recovery factor under different well interference conditions in a layered gas reservoir according to any one of claims 1-4.

Images