CN-121993186-A - Method and device for determining layering yield of gas well, storage medium and processor

CN121993186ACN 121993186 ACN121993186 ACN 121993186ACN-121993186-A

Abstract

The application relates to the field of natural gas exploitation and discloses a method, a device, a storage medium and a processor for determining layering output of a gas well, wherein the method comprises the steps of establishing a first corresponding relation between the minor axis dimension of each gas-bearing layer and the minor axis dimension of a first gas-bearing layer in each gas-bearing layer, and establishing a second corresponding relation between the major axis dimension of each gas-bearing layer and the minor axis dimension of the first gas-bearing layer; the method comprises the steps of determining dynamic reserves of each gas-bearing layer according to a first corresponding relation and a second corresponding relation, determining bottom hole flow pressures corresponding to each gas-bearing layer according to wellhead casing pressure of a gas well and related parameters of a gas well pipe column, and determining yield of each gas-bearing layer based on the dynamic reserves of each gas-bearing layer and the bottom hole flow pressures corresponding to each gas-bearing layer. According to the scheme provided by the application, the accuracy of calculating the yield of each longitudinal gas-bearing layer can be greatly improved by constructing the gas well layered yield determination scheme based on the oval seepage boundary.

Inventors

CHENG GANG
ZHU HANQING
CHEN KEFEI
LIU RUOHAN
CAO ZHENGLIN
GUO ZHI
JI GUANG
Han Jiangchen
WANG QIFENG
MENG DEWEI
WANG GUOTING
Cheng Minhua

Assignees

中国石油天然气股份有限公司

Dates

Publication Date: 20260508
Application Date: 20241106

Claims (15)

1. A method for determining the production of a gas well in layers, the method comprising: Regarding the cross section of each gas-containing layer corresponding to the ground as an ellipse, and determining the major axis dimension, the minor axis dimension and the thickness dimension of each gas-containing layer; Establishing a first corresponding relation between the minor axis dimension of each gas-containing layer and the minor axis dimension of a first gas-containing layer in each gas-containing layer based on the major axis dimension, the minor axis dimension and the thickness dimension of each gas-containing layer, and establishing a second corresponding relation between the major axis dimension of each gas-containing layer and the minor axis dimension of the first gas-containing layer; Determining the dynamic reserves of each gas-bearing layer according to the first corresponding relation and the second corresponding relation; Determining the corresponding bottom hole flow pressure of each gas-bearing layer according to the wellhead casing pressure of the gas well and the related parameters of the gas well pipe column; and determining the yield of each gas bearing layer based on the dynamic reserve of each gas bearing layer and the corresponding bottom hole flow pressure of each gas bearing layer.
2. The method of determining the production of gas well formations according to claim 1, wherein after the determining the production of each gas bearing layer, the method further comprises: Determining the simulated production of the gas well according to the production of each gas-bearing layer; and under the condition that the error between the simulated production of the gas well and the measured production of the gas well meets the preset condition, adjusting the short axis size of the first gas-bearing layer.
3. The method of determining the production of gas well formations according to claim 1, wherein after the determining the production of each gas bearing layer, the method further comprises: And determining the elliptical gas leakage range of each gas-containing layer based on the short axis dimension of the first gas-containing layer, the first corresponding relation and the second corresponding relation.
4. The method of determining the production of gas well formations according to claim 1, wherein the first gas bearing formation is the gas bearing formation closest to the surface of each gas bearing formation.
5. The method of determining the production of gas well formations according to claim 1, wherein establishing a first correspondence between the minor axis dimension of each gas bearing formation and the minor axis dimension of the first gas bearing formation, and establishing a second correspondence between the major axis dimension of each gas bearing formation and the minor axis dimension of the first gas bearing formation, based on the major axis dimension, the minor axis dimension, and the thickness dimension of each gas bearing formation, comprises: Fitting based on the short axis dimension of each gas-containing layer and the thickness dimension of each gas-containing layer to obtain a first fitting result, and establishing a first corresponding relation between the short axis dimension of each gas-containing layer and the short axis dimension of the first gas-containing layer according to the first fitting result; fitting is carried out based on the long axis size of each gas-containing layer and the short axis size of each gas-containing layer, a second fitting result is obtained, and a second corresponding relation between the long axis size of each gas-containing layer and the short axis size of the first gas-containing layer is established according to the second fitting result and the first corresponding relation.
6. The method of determining gas well stratified production as claimed in claim 1 wherein the dynamic reserves of each gas bearing layer are determined based on the following equation: The method comprises the steps of obtaining a dynamic reserve of a kth gas-containing layer, wherein G k is the dynamic reserve of the kth gas-containing layer, the unit is 10 4 m 3 ;L k is the long axis size of the kth gas-containing layer, the unit is m based on the second corresponding relation, W k is the short axis size of the kth gas-containing layer, the unit is m based on the first corresponding relation, h k is the thickness size of the kth gas-containing layer, the unit is m, phi k is the porosity of the kth gas-containing layer, S gik is the original gas-containing saturation of the kth gas-containing layer, and B gik is the natural gas volume coefficient corresponding to the original formation pressure of the kth gas-containing layer.
7. The method of determining gas well production in layers of claim 1, wherein the gas well string related parameters include a natural gas deviation factor at an average in-annulus temperature and an average in-annulus pressure, the method further comprising: Determining the bottom hole flow pressure corresponding to each gas bearing layer based on the following formula: Wherein, P wjk is the bottom hole flow pressure corresponding to the kth layer of gas-containing layer when the gas well is put into production, the unit is MPa, P cj is the well head casing pressure of the gas well when the gas well is put into production, the unit is MPa, r g is the relative density of natural gas, and H k is the depth of the target part of the kth layer of gas-containing layer, the unit is m; The unit is K, which is the average temperature in the annular space; natural gas deviation factor at average pressure in annulus at day j of production for gas well.
8. The method for determining the stratified production of a gas well as claimed in claim 7, wherein the natural gas deviation factor under the average pressure in the annulus at the j th day of gas well production is obtained based on the average pressure in the annulus at the j th day of gas well production The average annular pressure at the j th day of gas well production is obtained based on the following process: determining the bottom hole flow pressure corresponding to the second gas-containing layer at the initial time of gas well production; determining the association relation between the average pressure in the annular space at the j th day of gas well production and the bottom hole stream pressure corresponding to the second gas-containing layer at the j-1 th day of gas well production; and determining the average pressure in the annular space at the j th day of gas well production based on the bottom hole stream pressure corresponding to the second gas-containing layer at the initial time of gas well production and the association relation.
9. The method of determining the production of gas well formations according to claim 1, wherein the production of each gas bearing formation is determined based on the following equation before the production pressure drop propagates to the percolation boundary: In the above-mentioned formula(s), Wherein q jk is the yield of the gas-containing layer of the kth layer when the gas well is put into production, the unit is M 3 /d;τ j is the open time length of the gas well when the gas well is put into production, and the unit is h; P jk is the stratum pressure of the kth layer of gas-containing layer when the gas well is put into production on the jth day, and is obtained based on the dynamic reserve of the kth layer of gas-containing layer, wherein the unit is Mpa; P wjk is the bottom hole flow pressure corresponding to the kth layer gas-containing layer, mu is the natural gas viscosity, mu is mPas, Z jk is the natural gas deviation factor of the kth layer gas-containing layer, T k is the stratum temperature of the kth layer gas-containing layer, K is the standard state pressure, P sc is the Mpa, K is the effective permeability of the gas-containing layer, mD is the thickness dimension of the kth layer gas-containing layer, H k is M, Z sc is the standard state natural gas deviation factor, T sc is the standard state temperature, K is the short axis dimension of the kth layer gas-containing layer, obtained based on the first corresponding relation, W k is M, R w is the radius of a shaft, M is the molecular weight of air, R g is the relative density of natural gas, beta is the high-speed turbulence coefficient, M2 is the general purpose gas constant, and R is the standard state natural gas deviation factor, T3896K/(kK.m).
10. The method of determining the production of gas well formations according to claim 1, wherein after the production pressure drop has passed to the percolation boundary, the production of each gas bearing formation is determined based on the following equation: In the above-mentioned formula(s), Wherein q jk is the yield of the gas-containing layer of the kth layer when the gas well is put into production, the unit is M 3 /d;τ j is the open time length of the gas well when the gas well is put into production, and the unit is h; P jk is the stratum pressure of the kth layer of gas-containing layer when the gas well is put into production on the jth day, and is obtained based on the dynamic reserve of the kth layer of gas-containing layer, wherein the unit is Mpa; P wjk is the bottom hole flow pressure corresponding to the kth layer of gas-containing layer when the gas well is put into production, mu is the viscosity of natural gas, mu is the natural gas deviation factor of the kth layer of gas-containing layer when Z jk is the jth layer of gas well is put into production, T k is the stratum temperature of the kth layer of gas-containing layer, K is the standard state pressure, P sc is the Mpa, K is the effective permeability of the gas-containing layer, mu is the thickness dimension of the kth layer of gas-containing layer, M is the standard state natural gas deviation factor, Z sc is the standard state temperature, T sc is the standard state temperature, K is the short axis dimension of the kth layer of gas-containing layer, and is obtained based on the first corresponding relation, R w is the radius of a shaft, L k is the long axis dimension of the kth layer of gas-containing layer, K is obtained based on the second corresponding relation, M is the unit is the air molecular weight, K is the effective permeability of gas-containing layer is the unit of kg/mol, Z sc is the standard state natural gas deviation factor, T sc is the standard state temperature is the short axis dimension of the kth layer of gas, K is obtained based on the first corresponding relation, W k is the standard state gas is the K, R is the K, R is the general constant, K is the relative density of K is 3 , K is the K.
11. A gas well layered production determination apparatus, comprising: the acquisition module is used for regarding the section of each gas-containing layer corresponding to the ground as an ellipse and determining the major axis size, the minor axis size and the thickness size of each gas-containing layer; The corresponding relation establishing module is used for establishing a first corresponding relation between the minor axis dimension of each gas-containing layer and the minor axis dimension of a first gas-containing layer in each gas-containing layer based on the major axis dimension, the minor axis dimension and the thickness dimension of each gas-containing layer, and establishing a second corresponding relation between the major axis dimension of each gas-containing layer and the minor axis dimension of the first gas-containing layer; The dynamic reserve determining module is used for determining the dynamic reserve of each gas-bearing layer according to the first corresponding relation and the second corresponding relation; the bottom hole flow pressure determining module is used for determining the bottom hole flow pressure corresponding to each gas-containing layer according to the wellhead casing pressure of the gas well and the related parameters of the gas well pipe column; And the yield calculation module is used for determining the yield of each gas-bearing layer based on the dynamic reserves of each gas-bearing layer and the corresponding bottom hole flow pressure of each gas-bearing layer.
12. The gas well stratified production determination device of claim 11, further comprising: and the adjustment module is used for determining the simulated production of the gas well according to the production of each gas-containing layer, and adjusting the short axis size of the first gas-containing layer under the condition that the error between the simulated production of the gas well and the measured production of the gas well meets the preset condition.
13. The gas well stratified production determination device of claim 11, further comprising: And the air release range calculation module is used for determining the elliptical air release range of each air-bearing layer based on the short axis dimension of the first air-bearing layer, the first corresponding relation and the second corresponding relation.
14. A processor configured to perform the method of determining gas well production according to any one of claims 1 to 10.
15. A machine-readable storage medium having instructions stored thereon, which when executed by a processor, cause the processor to be configured to perform the method of gas well layered production determination of any of claims 1 to 10.

Description

Method and device for determining layering yield of gas well, storage medium and processor Technical Field The invention relates to the technical field of natural gas development, in particular to a method for determining layered production of a gas well, a device for determining layered production of the gas well, a machine-readable storage medium and a processor. Background Often gas reservoirs develop multiple gas bearing layers longitudinally, such as a lens lithologic trap gas reservoir of the Erdos basin. In the development of such reservoirs, to increase single well production, multiple gas bearing formations are often shot at the same time in the longitudinal direction. However, in the production process, in order to save production cost, a pressure gauge and a flow meter are not generally arranged for each longitudinal gas-containing layer, so that only the yield of a single well mouth can be measured, and the yields of each longitudinal gas-containing layer cannot be directly measured respectively. The acquisition of the yield of each longitudinal gas-bearing layer is a basis for evaluating the distribution of residual gas, and has important significance for adjusting the exploitation scheme. Therefore, in order to obtain the yield of each longitudinal gas-bearing layer, methods for indirectly determining the yield of each longitudinal gas-bearing layer, such as a stratum coefficient method, a mutation theory method, a numerical simulation method, a gas production profile test method and the like, are currently proposed. However, these methods are less reliable and it is difficult to accurately determine the production of each of the longitudinal gas bearing layers. Disclosure of Invention The invention aims to solve the problem that the prior art is difficult to accurately determine the production of each longitudinal gas bearing layer, and provides a method for determining the layered production of a gas well, a device for determining the layered production of the gas well, a machine-readable storage medium and a processor. To achieve the above object, a first aspect of the present invention provides a method for determining the production of a gas well in layers, the method comprising: Regarding the cross section of each gas-containing layer corresponding to the ground as an ellipse, and determining the major axis dimension, the minor axis dimension and the thickness dimension of each gas-containing layer; Establishing a first corresponding relation between the minor axis dimension of each gas-containing layer and the minor axis dimension of a first gas-containing layer in each gas-containing layer based on the major axis dimension, the minor axis dimension and the thickness dimension of each gas-containing layer, and establishing a second corresponding relation between the major axis dimension of each gas-containing layer and the minor axis dimension of the first gas-containing layer; Determining the dynamic reserves of each gas-bearing layer according to the first corresponding relation and the second corresponding relation; Determining the corresponding bottom hole flow pressure of each gas-bearing layer according to the wellhead casing pressure of the gas well and the related parameters of the gas well pipe column; and determining the yield of each gas bearing layer based on the dynamic reserve of each gas bearing layer and the corresponding bottom hole flow pressure of each gas bearing layer. In an embodiment of the present application, after the determining the production of each gas-bearing layer, the method further includes: Determining the simulated production of the gas well according to the production of each gas-bearing layer; and under the condition that the error between the simulated production of the gas well and the measured production of the gas well meets the preset condition, adjusting the short axis size of the first gas-bearing layer. In an embodiment of the present application, after the determining the production of each gas-bearing layer, the method further includes: And determining the elliptical gas leakage range of each gas-containing layer based on the short axis dimension of the first gas-containing layer, the first corresponding relation and the second corresponding relation. In the embodiment of the application, the first gas-bearing layer is the gas-bearing layer closest to the ground in each gas-bearing layer. In an embodiment of the present application, establishing a first correspondence between the minor axis dimension of each gas-containing layer and the minor axis dimension of the first gas-containing layer, and establishing a second correspondence between the major axis dimension of each gas-containing layer and the minor axis dimension of the first gas-containing layer based on the major axis dimension, the minor axis dimension, and the thickness dimension of each gas-containing layer includes: Fitting based on the short axis dimension of each gas-contai