Method, system, equipment and medium for judging water invasion and water seal strength of water-containing gas reservoir
By constructing a calculation model of water seal gas volume and converting it into a matter equilibrium equation of water-fighting reservoir, combined with the dynamic characteristic parameters of the energy region, the problem of water invasion and water seal strength in the existing technology is solved, and the problem of water invasion and water seal strength in the water invasion area is realized. Accurate determination and recovery rate of water-driven gas reservoirs are improved.
Patent Information
- Application Number
- CN202310966964.0
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2023-08-02
- Publication Date
- 2025-02-14
AI Technical Summary
The prior art fails to fully consider the determination and quantitative research on the rules of water invasion and water sealing of water flood reservoirs, which leads to the invasion of water invasion areas where water invasion is not accurately determined.
By obtaining the permeability parameters of the water invasion area, a calculation model of the water seal gas volume is constructed, and it is converted into the equilibrium equation of the water-fighting reservoir, a theoretical diagram is drawn and the energy area is divided, and the water invasion and water in the water invasion area is determined based on dynamic characteristic parameters. Sealing strength.
Accurate judgment of water invasion and water seal strength of water displaced reservoirs is achieved, and the influence of reservoir heterogeneity is taken into account, and the accuracy and credibility of the calculation model of water seal volume is improved. Countermeasures are taken in a timely manner to delay water invasion and reduce Small water sealing effect, thereby improving the recovery rate of water-domain reservoirs.
Smart Images

Figure CN119442557A_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to the technical field of oil and gas exploration and development, and in particular to a method, system, equipment and medium for determining water invasion and water seal strength of a water-containing gas reservoir. Background Art
[0002] This section is intended to provide a background or context to the embodiments recited in the claims. No admission is made that the description herein is prior art by inclusion in this section.
[0003] The efficient development of water-driven gas reservoirs has always been a difficult problem in the oil and gas industry. Over the years, domestic and foreign scholars have done a lot of research on water invasion laws, dynamic reserves, and water prevention and control in water-driven gas reservoirs. Among them, the material balance method has always been the basic method for gas reservoir engineering research. This method has been widely used in gas reservoir engineering. With the deepening of research on water-driven gas reservoirs in recent years, the material balance method has been continuously improved and optimized. It is mainly used for physical simulation and mathematical model solution of water invasion laws in water-driven gas reservoirs. For example, physical simulation methods have been used to carry out water invasion of porous, fractured, pore-cave and fracture-cave gas reservoirs. The law visualization experiment witnessed the seepage law of water in the reservoir, and further demonstrated the water control measures and methods to improve the recovery rate of heterogeneous reservoirs; based on the material balance principle of water-driven gas reservoirs, considering the influence of water-sealed gas, an automatic fitting mathematical model for the solution method of dynamic reserves, water-sealed gas volume, and water invasion volume of water-driven gas reservoirs was established; under different conditions, based on the material balance principle, different methods were used to successively establish dynamic reserve calculation models for water-driven gas reservoirs, all of which were effectively demonstrated, but the above methods did not consider the determination and quantitative research of water invasion and water sealing laws of water-driven gas reservoirs under heterogeneous reservoir conditions. In summary, the existing technology does not fully consider the determination and quantitative research of water invasion and water sealing laws of water-driven gas reservoirs, which leads to the problem that the strength of water invasion and water sealing in water-invaded gas reservoirs cannot be determined. Summary of the invention
[0004] In view of the above problems, embodiments of the present invention provide a method, system, device and medium for determining water invasion and water seal strength in a water-containing gas reservoir.
[0005] In a first aspect, an embodiment of the present invention provides a method for determining water invasion and water seal strength in a water-containing gas reservoir, comprising:
[0006] Obtaining permeability parameters of the water intrusion zone, and using the permeability parameters to construct a water-sealed air volume calculation model;
[0007] Performing data conversion on the water-seal gas volume calculation model to obtain a material balance equation for a water-drive gas reservoir;
[0008] Draw a theoretical map according to the material balance equation of the water-drive gas reservoir, divide the theoretical map into energy regions;
[0009] The water invasion and water sealing strength of the water-containing gas reservoir in the water invasion zone are determined according to the dynamic characteristic parameters in the energy region.
[0010] According to an embodiment of the present invention, the step of constructing a water-sealed air volume calculation model using the permeability parameter includes:
[0011] Selecting target permeability parameters of a high permeability layer and a low permeability layer from the permeability parameters, calculating the target permeability parameters, and obtaining the water body seepage rate;
[0012] Obtaining the water-sealed gas volume of the low permeability layer and the gas volume coefficient under the current formation conditions, performing ratio calculation on the target permeability parameter, the gas volume coefficient under the current formation conditions and the water-sealed gas volume to obtain the pore volume ratio;
[0013] Convert the expression of the pore volume ratio to obtain the expression of the water-sealed gas volume;
[0014] Obtaining the produced water volume of the gas reservoir, adding and calculating the target permeability parameter and the produced water volume of the gas reservoir to obtain the water invasion volume of the gas reservoir;
[0015] Setting the permeability ratio of the high permeability layer and the low permeability layer, combining the water body seepage, the pore volume ratio, the expression of the water-sealed gas volume, the gas reservoir water intrusion and the permeability ratio to obtain a target equation for the water-sealed gas volume;
[0016] Obtaining the gas volume coefficient, the dynamic reserves of the gas reservoir, and the water storage volume coefficient of the gas reservoir under the original formation conditions, and calculating the gas reservoir equation corresponding to the water storage volume coefficient of the gas reservoir according to the gas volume coefficient, the dynamic reserves of the gas reservoir, and the water storage volume coefficient of the gas reservoir under the original formation conditions;
[0017] A water-sealed gas volume calculation model is generated according to the target equation for the water-sealed gas volume and the gas reservoir equation.
[0018] According to an embodiment of the present invention, generating a water-sealed gas volume calculation model according to the target equation of the water-sealed gas volume and the gas reservoir equation includes:
[0019] According to the target equation of the water-sealed gas volume, a reservoir characteristic constant is set, and the reservoir characteristic constant, the target equation and the gas reservoir equation are combined to obtain a water-sealed gas volume calculation model;
[0020] The reservoir characteristic constant is expressed as:
[0021]
[0022] Wherein, A represents the reservoir characteristic constant, V represents the pore volume ratio, and K represents the permeability ratio;
[0023] The water seal gas volume calculation model is expressed as:
[0024]
[0025] Among them, G F represents the water-sealed gas volume, G represents the dynamic reserves of the gas reservoir, and B g0 represents the gas volume coefficient under the original formation conditions, B g represents the gas volume coefficient under the current formation conditions, A represents the reservoir characteristic constant, and t represents the water volume coefficient of the gas reservoir.
[0026] According to an embodiment of the present invention, the data conversion of the water-seal gas volume calculation model to obtain the material balance equation of the water-drive gas reservoir includes:
[0027] Obtaining a gas volume coefficient in the water-seal gas volume calculation model, and updating the gas volume coefficient according to a preset formation pressure to obtain an updated gas volume coefficient;
[0028] A dimensionless pseudo-pressure is generated according to the updated gas volume coefficient, and a material balance equation of a water-drive gas reservoir is generated according to the dimensionless pseudo-pressure and the water-seal gas volume calculation model.
[0029] According to an embodiment of the present invention, generating a material balance equation for a water-drive gas reservoir according to the dimensionless pseudo-pressure and the water-seal gas volume calculation model includes:
[0030] Obtaining water invasion characteristic constants, current cumulative gas production and gas reservoir production degree, extracting gas reservoir water storage volume coefficient and gas reservoir dynamic reserves in the water-sealed gas volume calculation model, and establishing an equation based on the water invasion characteristic constants, the current cumulative gas production, the gas reservoir production degree, the gas reservoir water storage volume coefficient and the gas reservoir dynamic reserves;
[0031] A material balance equation for a water-driven gas reservoir is established based on the equation, the dimensionless pseudo-pressure and the water-seal gas volume calculation model.
[0032] According to an embodiment of the present invention, the step of drawing a theoretical diagram according to the material balance equation of the water drive gas reservoir includes:
[0033] Performing parameter analysis on the material balance equation of the water drive gas reservoir to obtain relevant parameters;
[0034] The dimensionless pseudo-pressure and gas reservoir recovery degree in the relevant parameters are obtained, and based on the material balance equation of the water-driven gas reservoir, a theoretical chart is generated by using a preset production decline analysis method with the dimensionless pseudo-pressure as a column vector and the gas reservoir recovery degree as a row vector.
[0035] According to an embodiment of the present invention, dividing the theoretical map to obtain energy regions includes:
[0036] Extracting a reservoir characteristic constant and a water invasion characteristic constant from the material balance equation of the water drive gas reservoir, and performing a numerical judgment on the reservoir characteristic constant when the water invasion characteristic constant is not 0;
[0037] When the reservoir characteristic constant is 0, the theoretical map is divided into an energy replenishment zone;
[0038] When the reservoir characteristic constant is greater than 0 and less than 1, the theoretical map is divided into an energy replenishment dominant area;
[0039] When the reservoir characteristic constant is 1, the theoretical map is divided on the equilibrium line;
[0040] When the reservoir characteristic constant is greater than 1, the theoretical map is divided into an energy isolation dominant area;
[0041] The energy replenishment area, the energy replenishment leading area, the balance line and the energy isolation leading area together constitute an energy area.
[0042] In a second aspect, an embodiment of the present invention provides a water invasion and water seal strength determination system for a water-containing gas reservoir, characterized in that it includes:
[0043] A calculation model building module, used to obtain the permeability parameters of the water intrusion area and build a water-sealed air volume calculation model using the permeability parameters;
[0044] A data conversion module, used for performing data conversion on the water-seal gas volume calculation model to obtain a material balance equation for a water-drive gas reservoir;
[0045] An energy region generation module is used to draw a theoretical map according to the material balance equation of the water-drive gas reservoir, divide the theoretical map, and obtain energy regions;
[0046] The strength determination module is used to determine the water invasion and water sealing strength of the water-containing gas reservoir in the water invasion area according to the dynamic characteristic parameters in the energy area.
[0047] In a third aspect, an embodiment of the present invention provides an electronic device, comprising:
[0048] processor;
[0049] a memory for storing instructions executable by the processor;
[0050] The processor is configured to execute the instructions to implement a method for determining water invasion and water seal strength in a water-containing gas reservoir as described in the first aspect.
[0051] In a fourth aspect, an embodiment of the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements a method for determining water invasion and water seal strength of a water-containing gas reservoir as described in the first aspect.
[0052] Compared with the prior art, the above technical solution of the present invention has the following beneficial effects:
[0053] The embodiment of the present invention obtains the permeability parameters of the water invasion zone and constructs a water-sealed gas volume calculation model using the permeability parameters, taking into account the influence of reservoir heterogeneity during the water invasion process, and ensuring that the water-sealed gas volume calculation model is more accurate; by performing data conversion on the water-sealed gas volume calculation model, the water-sealed gas volume is introduced into the material balance equation, and a water-drive gas reservoir material balance equation considering water-sealed gas conditions is established, so that the water-drive gas reservoir material balance equation is more accurate and credible; a theoretical map is drawn through the water-drive gas reservoir material balance equation, and the theoretical map is divided to obtain energy regions, which combines the actual conditions of water invasion and water sealing, making the theoretical map more realistic; the water invasion and water sealing strength of the water-containing gas reservoir in the water invasion zone are determined through the dynamic characteristic parameters in the energy region, the determination and quantitative research of the water invasion and water sealing laws of the water-drive gas reservoir are fully considered, and corresponding countermeasures are taken in time to delay water invasion, reduce the water sealing effect to a large extent, and thus improve the recovery rate of the water-containing gas reservoir. BRIEF DESCRIPTION OF THE DRAWINGS
[0054] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.
[0055] Figure 1 The flowchart of the method for determining water invasion and water seal strength of a water-containing gas reservoir according to the first embodiment of the present invention is shown;
[0056] Figure 2 It shows a schematic structural diagram of a water-driven gas reservoir physical model under water-sealed gas conditions according to the first embodiment of the present invention;
[0057] Figure 3 A schematic diagram showing the structure of a theoretical diagram of the first embodiment of the present invention;
[0058] Figure 4 The functional module diagram of the water invasion and water seal strength determination system of the water-containing gas reservoir according to the third embodiment of the present invention is shown;
[0059] Figure 5 A schematic diagram of the composition structure of an electronic device for implementing the method for determining water invasion and water seal strength in a water-containing gas reservoir according to the fourth embodiment of the present invention is shown. DETAILED DESCRIPTION
[0060] The present invention will be further described below in conjunction with the embodiments shown in the accompanying drawings.
[0061] It should be noted that, in the absence of conflict, the embodiments and features in the embodiments of the present application can be combined with each other. The present invention will be described in detail below with reference to the accompanying drawings and in combination with the embodiments.
[0062] The present invention proposes a method for determining water invasion and water seal strength of water-containing gas reservoirs. Based on the permeability parameters of the water invasion zone, a water seal gas volume calculation model is constructed, and the material balance equation of the water-driven gas reservoir is converted according to the water seal gas volume calculation model. A theoretical plate is drawn according to the material balance equation of the water-driven gas reservoir, and the energy region is divided using the theoretical plate. The water invasion and water seal strength of the water-containing gas reservoir in the water invasion zone is determined according to the dynamic characteristic parameters in the energy region. Compared with the traditional method, the method for determining water invasion and water seal strength of water-containing gas reservoirs fully considers the determination and quantitative research of water invasion and water seal laws of water-driven gas reservoirs, and can further determine the strength of water invasion and water seal of water-containing gas reservoirs in the water invasion zone.
[0063] Embodiment 1
[0064] like Figure 1 As shown, the present invention proposes a method for determining water invasion and water seal strength of a water-containing gas reservoir, comprising the following steps:
[0065] S1. Obtaining permeability parameters of the water intrusion area, and using the permeability parameters to construct a water-sealed air volume calculation model;
[0066] In the embodiment of the present invention, the permeability parameters of the water invasion zone in the preset water drive gas reservoir physical model are obtained, such as Figure 2 The physical model of the water-driven gas reservoir under water-sealed gas conditions shown in the figure includes a high permeability layer and a low permeability layer on the plane of the gas reservoir reservoir, and water invading the high permeability layer to partially isolate the gas in the low permeability layer, namely, water-sealed gas. After the formation pressure drops to a certain level, the water-sealed gas expands and drives the formation water, and the volume of the water-sealed gas changes continuously with the formation pressure. The permeability parameters include water seepage rate, permeability, seepage area and seepage distance, etc.
[0067] In the embodiment of the present invention, the method of constructing a water-sealed air volume calculation model using the permeability parameter includes:
[0068] Selecting target permeability parameters of a high permeability layer and a low permeability layer from the permeability parameters, calculating the target permeability parameters, and obtaining the water body seepage rate;
[0069] Obtaining the water-sealed gas volume of the low permeability layer and the gas volume coefficient under the current formation conditions, performing ratio calculation on the target permeability parameter, the gas volume coefficient under the current formation conditions and the water-sealed gas volume to obtain the pore volume ratio;
[0070] Convert the expression of the pore volume ratio to obtain the expression of the water-sealed gas volume;
[0071] Obtaining the produced water volume of the gas reservoir, adding and calculating the target permeability parameter and the produced water volume of the gas reservoir to obtain the water invasion volume of the gas reservoir;
[0072] Setting the permeability ratio of the high permeability layer and the low permeability layer, combining the water body seepage, the pore volume ratio, the expression of the water-sealed gas volume, the gas reservoir water intrusion and the permeability ratio to obtain a target equation for the water-sealed gas volume;
[0073] Obtaining the gas volume coefficient, the dynamic reserves of the gas reservoir, and the water storage volume coefficient of the gas reservoir under the original formation conditions, and calculating the gas reservoir equation corresponding to the water storage volume coefficient of the gas reservoir according to the gas volume coefficient, the dynamic reserves of the gas reservoir, and the water storage volume coefficient of the gas reservoir under the original formation conditions;
[0074] A water-sealed gas volume calculation model is generated according to the target equation for the water-sealed gas volume and the gas reservoir equation.
[0075] In the embodiment of the present invention, the target infiltration parameter is calculated to obtain the water seepage rate, including:
[0076] The water seepage volume is calculated using the following formula:
[0077]
[0078]
[0079] Among them, ω H represents the water seepage rate of the high permeability layer (H), ω I represents the water seepage of the low permeability layer (I), K H represents the permeability of the target permeability parameter of the high permeability layer (H), K I represents the permeability of the target permeability parameter of the low permeability layer (I), A H represents the seepage area in the target permeability parameter of the high permeability layer (H), A I represents the seepage area in the target seepage parameter of the low permeability layer (I), ΔL represents the seepage distance in the target seepage parameter, μ w represents the preset formation water viscosity, and Δp represents the preset production pressure difference.
[0080] In the embodiment of the present invention, the ratio calculation of the target permeability parameter, the gas volume coefficient under the current formation condition and the water-sealed gas volume to obtain the pore volume ratio includes:
[0081] The pore volume ratio is calculated using the following formula:
[0082]
[0083]
[0084] Wherein, V represents the pore volume ratio, V I represents the pore volume of the low permeability layer (I) in the target permeability parameter, V H represents the pore volume of the high permeability layer (H) in the target permeability parameter, G F Represents the water seal air volume, W I represents the water storage volume of the low permeability layer (I) in the target permeability parameter, W H represents the water storage volume of the high permeability layer (H) in the target permeability parameter, B g Indicates the gas volume coefficient under the current formation conditions.
[0085] In the embodiment of the present invention, the expression of the pore volume ratio is converted to obtain the expression of the water-sealed air volume, including:
[0086] The expression of the water seal air volume is expressed as:
[0087]
[0088] Among them, G F Indicates the water seal air volume, B g represents the gas volume coefficient under the current formation conditions, V I represents the pore volume of the low permeability layer (I) in the target permeability parameter, W H represents the pore volume of the high permeability layer (H) in the target permeability parameter, W I represents the water storage volume of the low permeability layer (I) in the target permeability parameter, W H It represents the water storage volume of the high permeability layer (H) in the target permeability parameter.
[0089] In the embodiment of the present invention, the target permeability parameter and the gas reservoir produced water volume are added and calculated to obtain the gas reservoir water invasion volume, including:
[0090] The water invasion of the gas reservoir is obtained by adding and calculating using the following formula:
[0091] W e =W H +W I +W p B a
[0092] Among them, We represents the water intrusion of the gas reservoir, W I represents the water storage volume of the low permeability layer (I) in the target permeability parameter, W H represents the water storage volume of the high permeability layer (H) in the target permeability parameter, W p represents the produced water volume of the gas reservoir, B a Indicates the preset formation water volume coefficient.
[0093] In the embodiment of the present invention, the water seepage rate, the pore volume ratio, the expression of the water-sealed gas volume, the gas reservoir water intrusion and the permeability ratio are combined to obtain the target equation of the water-sealed gas volume, including:
[0094]
[0095]
[0096]
[0097]
[0098]
[0099]
[0100] The target equation is expressed as:
[0101]
[0102] Among them, G F Represents the water seal air volume, W I represents the water storage volume of the low permeability layer (I) in the target permeability parameter, W H represents the water storage volume of the high permeability layer (H) in the target permeability parameter, ω H represents the water seepage rate of the high permeability layer (H), ω I represents the water seepage of the low permeability layer (I), K H represents the permeability of the target permeability parameter of the high permeability layer (H), K I represents the permeability of the target permeability parameter of the low permeability layer (I), A H represents the seepage area in the target permeability parameter of the high permeability layer (H), A I represents the seepage area in the target permeability parameter of the low permeability layer (I), V represents the pore volume ratio, K represents the permeability ratio, W e represents the water intrusion of the gas reservoir, W p represents the produced water volume of the gas reservoir, B arepresents the formation water volume coefficient, B g Indicates the gas volume coefficient under the current formation conditions.
[0103] In the embodiment of the present invention, the gas reservoir equation corresponding to the gas reservoir water storage volume coefficient is calculated according to the gas volume coefficient, the dynamic reserves of the gas reservoir and the water storage volume coefficient of the gas reservoir under the original formation conditions, including:
[0104] The gas reservoir equation corresponding to the gas reservoir water volume coefficient is expressed as:
[0105]
[0106] Wherein, t represents the water volume coefficient of the gas reservoir, G represents the dynamic reserves of the gas reservoir, and B g0 It represents the gas volume coefficient under the original formation conditions, W e represents the water intrusion of the gas reservoir, W p represents the produced water volume of the gas reservoir, B a represents the formation water volume coefficient.
[0107] In an embodiment of the present invention, the step of generating a water-seal gas volume calculation model according to the target equation for the water-seal gas volume and the gas reservoir equation includes:
[0108] According to the target equation of the water-sealed gas volume, a reservoir characteristic constant is set, and the reservoir characteristic constant, the target equation and the gas reservoir equation are combined to obtain a water-sealed gas volume calculation model;
[0109] The reservoir characteristic constant is expressed as:
[0110]
[0111] Wherein, A represents the reservoir characteristic constant, V represents the pore volume ratio, and K represents the permeability ratio;
[0112] The water seal gas volume calculation model is expressed as:
[0113]
[0114] Among them, G F represents the water-sealed gas volume, G represents the dynamic reserves of the gas reservoir, and B g0 represents the gas volume coefficient under the original formation conditions, B g represents the gas volume coefficient under the current formation conditions, A represents the reservoir characteristic constant, and t represents the water volume coefficient of the gas reservoir.
[0115] S2, performing data conversion on the water-seal gas volume calculation model to obtain a material balance equation for a water-drive gas reservoir;
[0116] In the embodiment of the present invention, the data conversion of the water-seal gas volume calculation model to obtain the material balance equation of the water-drive gas reservoir includes:
[0117] Obtaining a gas volume coefficient in the water-seal gas volume calculation model, and updating the gas volume coefficient according to a preset formation pressure to obtain an updated gas volume coefficient;
[0118] A dimensionless pseudo-pressure is generated according to the updated gas volume coefficient, and a material balance equation of a water-drive gas reservoir is generated according to the dimensionless pseudo-pressure and the water-seal gas volume calculation model.
[0119] In the embodiment of the present invention, the gas volume coefficient is replaced according to the preset formation pressure to obtain an updated gas volume coefficient, including:
[0120] The formation pressure formula is expressed as:
[0121] B′ g =P1 / Z1
[0122] B g1 =P / Z
[0123] Among them, B′ g represents the updated gas volume coefficient under the current formation conditions, B g1 represents the updated gas volume coefficient under original formation conditions, P1 represents the original formation pressure, Z1 represents the gas deviation factor of the original formation pressure, P represents the current formation pressure, and Z represents the gas deviation factor of the current formation pressure.
[0124] In the embodiment of the present invention, generating dimensionless pseudo-pressure according to the updated gas volume coefficient includes:
[0125] The dimensionless pseudo-pressure is expressed as:
[0126]
[0127]
[0128] in, represents the dimensionless pseudo-pressure, B′ g represents the updated gas volume coefficient under the current formation conditions, B g1 represents the updated gas volume coefficient under original formation conditions, P1 represents the original formation pressure, Z1 represents the gas deviation factor of the original formation pressure, P represents the current formation pressure, and Z represents the gas deviation factor of the current formation pressure.
[0129] In the embodiment of the present invention, the material balance equation of the water-drive gas reservoir is generated according to the dimensionless pseudo-pressure and the water-seal gas volume calculation model, including:
[0130] Obtaining water invasion characteristic constants, current cumulative gas production and gas reservoir production degree, extracting gas reservoir water storage volume coefficient and gas reservoir dynamic reserves in the water-sealed gas volume calculation model, and establishing an equation based on the water invasion characteristic constants, the current cumulative gas production, the gas reservoir production degree, the gas reservoir water storage volume coefficient and the gas reservoir dynamic reserves;
[0131] A material balance equation for a water-driven gas reservoir is established based on the equation, the dimensionless pseudo-pressure and the water-seal gas volume calculation model.
[0132] In the embodiment of the present invention, the equation is expressed as:
[0133] ln t=B ln R
[0134]
[0135] Wherein, t represents the water volume coefficient of the gas reservoir, B represents the water invasion characteristic constant, R represents the gas reservoir recovery degree, G b represents the current cumulative gas production, and G represents the dynamic reserves of the gas reservoir;
[0136] The material balance equation of the water drive gas reservoir is expressed as:
[0137] GB g1 =(GG F -G b )B′ g +(W e -W p B a )
[0138]
[0139]
[0140]
[0141] in, represents the dimensionless pseudo-pressure, G represents the dynamic reserves of the gas reservoir, G F Represents the water seal gas volume, G b Indicates the current cumulative gas production, B′ g represents the updated gas volume coefficient under the current formation conditions, B g1 represents the renewal gas volume coefficient under the original formation conditions, W e represents the water intrusion of the gas reservoir, W p represents the produced water volume of the gas reservoir, B arepresents the formation water volume coefficient, R represents the gas reservoir recovery degree, A represents the reservoir characteristic constant, B represents the water invasion characteristic constant, p1 represents the original formation pressure, Z1 represents the gas deviation factor of the original formation pressure, P represents the current formation pressure, Z represents the gas deviation factor of the current formation pressure, R B The gas reservoir recovery degree representing the water invasion characteristic constant.
[0142] In the embodiment of the present invention, the material balance equation of the water-drive gas reservoir is generated according to the material balance principle through the equation, the dimensionless pseudo-pressure and the water-seal gas volume calculation model, and the condition that part of the gas is sealed is taken into consideration, so that the obtained material balance equation of the water-drive gas reservoir is more accurate and the credibility of the material balance equation of the water-drive gas reservoir is improved.
[0143] S3, drawing a theoretical map according to the material balance equation of the water-drive gas reservoir, dividing the theoretical map to obtain energy regions;
[0144] In the embodiment of the present invention, the step of drawing a theoretical diagram according to the material balance equation of the water drive gas reservoir includes:
[0145] Performing parameter analysis on the material balance equation of the water drive gas reservoir to obtain relevant parameters;
[0146] The dimensionless pseudo-pressure and gas reservoir recovery degree in the relevant parameters are obtained, and based on the material balance equation of the water-driven gas reservoir, a theoretical chart is generated by using a preset production decline analysis method with the dimensionless pseudo-pressure as a column vector and the gas reservoir recovery degree as a row vector.
[0147] In an embodiment of the present invention, the production decline analysis method includes but is not limited to the Blasingame method; according to the material balance equation of the water-drive gas reservoir, the curve change trend in the theoretical plate under different parameter conditions can be obtained, and the theoretical plate refers to the water invasion characteristic curve plate of the water-drive gas reservoir under water-sealed gas conditions.
[0148] In the embodiment of the present invention, the dividing of the theoretical map to obtain energy regions includes:
[0149] Extracting a reservoir characteristic constant and a water invasion characteristic constant from the material balance equation of the water drive gas reservoir, and performing a numerical judgment on the reservoir characteristic constant when the water invasion characteristic constant is not 0;
[0150] When the reservoir characteristic constant is 0, the theoretical map is divided into an energy replenishment zone;
[0151] When the reservoir characteristic constant is greater than 0 and less than 1, the theoretical map is divided into an energy replenishment dominant area;
[0152] When the reservoir characteristic constant is 1, the theoretical chart is divided on the balance line;
[0153] When the reservoir characteristic constant is greater than 1, the theoretical chart is divided into an energy isolation dominant area;
[0154] The energy replenishment area, the energy replenishment dominant area, the balance line and the energy isolation dominant area together form an energy area.
[0155] In the embodiment of the present invention, as Figure 3 shown, the theoretical chart refers to the water invasion characteristic curve chart of a water drive gas reservoir under the condition of considering water sealing gas. Under different water invasion and water sealing intensities, the theoretical chart is mainly divided into four energy areas, namely the energy replenishment area, the energy replenishment dominant area, the balance line and the energy isolation dominant area, which are divided according to the reservoir characteristic constant (A) and the water invasion characteristic constant (B). For example, in the energy replenishment area: A = 0, B ≠ 0; in the energy replenishment dominant area: 0 < A < 1, B ≠ 0; on the balance line: A = 1, B ≠ 0, and the balance line includes but is not limited to the replenishment line and the isolation line; in the energy isolation dominant area: A > 1, B ≠ 0.
[0156] S4. Determine the water invasion and water sealing intensity of the water-bearing gas reservoir in the water invasion area according to the dynamic characteristic parameters in the energy area.
[0157] In the embodiment of the present invention, the dynamic characteristic parameters refer to the reservoir characteristic constant (A) and the water invasion characteristic constant (B). The trend of the water invasion characteristic curve is judged according to the dynamic characteristic parameters. When the water invasion characteristic curve shows an upward trend, the water invasion intensity is greater; when the water invasion characteristic curve shows a downward trend, the water sealing intensity is greater, so as to determine the water invasion and water sealing intensity of the water invasion reservoir in the water invasion area. For example, when A = 0, B ≠ 0, the gas reservoir is a homogeneous gas reservoir, the water invasion characteristic is uniform water invasion, the water invasion characteristic curve is "concave upward", and the more severe the upward warping, the greater the water invasion intensity; when 0 < A < 1, B ≠ 0, water sealing gas is formed in the water invasion area, indicating that water quickly invades along the high-permeability layer, and part of the natural gas in the low-permeability layer is sealed before it can be driven, weakening the replenishment of water invasion energy to a certain extent, and the water invasion characteristic curve gradually shows a "convex downward" trend; when A = 1, B ≠ 0, the balance line includes but is not limited to the replenishment line and the isolation line, and the water invasion energy replenishment and the water sealing energy isolation reach balance; when A > 1, B ≠ 0, the gas reservoir has strong heterogeneity. After a certain degree of exploitation, a strong water sealing gas phenomenon appears. At that time, a large part of the gas is sealed off, the energy is weakened quickly, and the water invasion characteristic curve shows an "concave downward" shape, and the faster the downward drop, the more obvious the water sealing gas phenomenon.
[0158] Embodiment 2
[0159] To understand the present invention more clearly, the water invasion and water-sealing intensity of a water-bearing gas reservoir in the water invasion area will be further explained through a second embodiment below in the case of determining according to the dynamic characteristic parameters in the energy area.
[0160] In the embodiment of the present invention, determining the water invasion and water-sealing intensity of a water-bearing gas reservoir in the water invasion area according to the dynamic characteristic parameters in the energy area includes:
[0161] Determining the trend of the water invasion characteristic curve in the energy area according to the dynamic characteristic parameters;
[0162] When the water invasion characteristic curve is upward, the water invasion intensity of the water-bearing gas reservoir in the water invasion area is greater than the water-sealing intensity of the water invasion area;
[0163] When the water invasion characteristic curve is downward, the water-sealing intensity of the water invasion area is greater than the water invasion intensity of the water-bearing gas reservoir in the water invasion area;
[0164] When the water invasion characteristic curve is at the equilibrium line, the water invasion and water-sealing intensity of the water-bearing gas reservoir in the water invasion area are the same.
[0165] In the embodiment of the present invention, the dynamic characteristic parameters refer to the reservoir characteristic constant (A) and the water invasion characteristic constant (B). According to the dynamic characteristic parameters, the trend of the water invasion characteristic curve is judged. When the water invasion characteristic curve shows an upward trend, the water invasion intensity is greater; when the water invasion characteristic curve shows a downward trend, the water-sealing intensity is greater, so as to judge the water invasion and water-sealing intensity of the water invasion reservoir in the water invasion area. For example, when A = 0 and B ≠ 0, the gas reservoir is a homogeneous gas reservoir, the water invasion characteristic is uniform water invasion, the water invasion characteristic curve is "concave upward", and the more severe the upward warping, the greater the water invasion intensity; when 0 < A < 1 and B ≠ 0, water-sealed gas is formed in the water invasion area, indicating that water quickly invades along the high-permeability layer, and part of the natural gas in the low-permeability layer is sealed before being driven, weakening the replenishment of water invasion energy to a certain extent, and the water invasion characteristic curve gradually shows a "convex downward" trend; when A = 1 and B ≠ 0, the equilibrium line includes but is not limited to the replenishment line and the sealing line, and the balance between water invasion energy replenishment and water-sealing energy sealing is achieved; when A > 1 and B ≠ 0, the heterogeneity of the gas reservoir is relatively strong. After mining to a certain extent, a strong water-sealed gas phenomenon appears. At that time, a large part of the gas is sealed off, the energy is weakened quickly, the water invasion characteristic curve shows a "concave downward", and the faster the downward drop, the more obvious the water-sealed gas phenomenon.
[0166] Embodiment III
[0167] As Figure 4 shown, this embodiment also provides a functional module diagram of a system for determining the water invasion and water-sealing intensity of a water-bearing gas reservoir.
[0168] The water invasion and water seal strength determination system 400 for water-containing gas reservoirs described in this embodiment can be installed in an electronic device. According to the functions to be implemented, the water invasion and water seal strength determination system 400 for water-containing gas reservoirs can include a calculation model construction module 401, a data conversion module 402, an energy region generation module 403 and a strength determination module 404. The module described in the present invention can also be called a unit, which refers to a series of computer program segments that can be executed by an electronic device processor and can complete fixed functions, which are stored in the memory of the electronic device.
[0169] In this embodiment, the functions of each module / unit are as follows:
[0170] The calculation model building module 401 is used to obtain the permeability parameters of the water intrusion area and use the permeability parameters to build a water-sealed air volume calculation model;
[0171] The data conversion module 402 is used to perform data conversion on the water-seal gas volume calculation model to obtain a material balance equation for a water-drive gas reservoir;
[0172] The energy region generation module 403 is used to draw a theoretical map according to the material balance equation of the water drive gas reservoir, divide the theoretical map, and obtain energy regions;
[0173] The strength determination module 404 is used to determine the water invasion and water sealing strength of the water-containing gas reservoir in the water invasion area according to the dynamic characteristic parameters in the energy area.
[0174] In detail, the modules described in the water invasion and water seal strength determination system 400 for water-containing gas reservoirs described in the embodiment of the present invention adopt the same technical means as the water invasion and water seal strength determination method for water-containing gas reservoirs described in Example 1 and Example 2 when used, and can produce the same technical effects, which will not be repeated here.
[0175] Embodiment 4
[0176] like Figure 5 As shown, this embodiment also provides a computer electronic device, and the electronic device 500 may include a processor 501, a memory 502, a communication bus 503 and a communication interface 504, and may also include a computer program stored in the memory 502 and executable on the processor 501, such as a water invasion and water seal strength determination program for a water gas reservoir.
[0177] In some embodiments, the processor 501 may be composed of an integrated circuit, for example, a single packaged integrated circuit, or a plurality of packaged integrated circuits with the same or different functions, including one or more central processing units (CPUs), microprocessors, digital processing chips, graphics processors, and combinations of various control chips, etc. The processor 501 is the control core (ControlUnit) of the electronic device, and uses various interfaces and lines to connect various components of the entire electronic device, and executes or executes programs or modules stored in the memory 502 (for example, executing water invasion and water seal strength determination programs for water-containing gas reservoirs, etc.), and calls data stored in the memory 502 to execute various functions of the electronic device and process data.
[0178] The memory 502 includes at least one type of readable storage medium, and the readable storage medium includes a flash memory, a mobile hard disk, a multimedia card, a card-type memory (for example, SD or DX memory, etc.), a magnetic memory, a disk, an optical disk, etc. In some embodiments, the memory 502 may be an internal storage unit of an electronic device, such as a mobile hard disk of the electronic device. In other embodiments, the memory 502 may also be an external storage device of an electronic device, such as a plug-in mobile hard disk, a smart memory card (Smart Media Card, SMC), a secure digital (Secure Digital, SD) card, a flash card (Flash Card), etc. equipped on the electronic device. Further, the memory 502 may also include both an internal storage unit of the electronic device and an external storage device. The memory 502 may not only be used to store application software and various types of data installed in the electronic device, such as the code of the water invasion and water seal strength determination program of the water gas reservoir, but also be used to temporarily store data that has been output or is to be output.
[0179] The communication bus 503 may be a peripheral component interconnect (PCI) bus or an extended industry standard architecture (EISA) bus, etc. The bus may be divided into an address bus, a data bus, a control bus, etc. The bus is configured to realize connection and communication between the memory 502 and at least one processor 501, etc.
[0180] The communication interface 504 is used for communication between the above-mentioned electronic device and other devices, including a network interface and a user interface. Optionally, the network interface may include a wired interface and / or a wireless interface (such as a WI-FI interface, a Bluetooth interface, etc.), which is generally used to establish a communication connection between the electronic device and other electronic devices. The user interface may be a display (Display), an input unit (such as a keyboard (Keyboard)), and optionally, the user interface may also be a standard wired interface, a wireless interface. Optionally, in some embodiments, the display may be an LED display, a liquid crystal display, a touch-sensitive liquid crystal display, and an OLED (Organic Light-Emitting Diode, organic light-emitting diode) touch device, etc. Among them, the display may also be appropriately referred to as a display screen or a display unit, which is used to display information processed in the electronic device and to display a visual user interface.
[0181] The figure only shows an electronic device 500 with components. Those skilled in the art will understand that the structure shown in the figure does not constitute a limitation on the electronic device, and may include fewer or more components than shown in the figure, or combine certain components, or arrange the components differently.
[0182] For example, although not shown, the electronic device may also include a power source (such as a battery) for supplying power to each component. Preferably, the power source may be logically connected to the at least one processor 501 through a power management device, so that the power management device can realize functions such as charging management, discharging management, and power consumption management. The power source may also include any components such as one or more DC or AC power sources, recharging devices, power failure detection circuits, power converters or inverters, power status indicators, etc. The electronic device may also include a variety of sensors, Bluetooth modules, Wi-Fi modules, etc., which will not be repeated here.
[0183] It should be understood that the embodiment is for illustration only and the scope of the patent application is not limited to this structure.
[0184] The water invasion and water seal strength determination program for water-containing gas reservoirs stored in the memory 502 in the electronic device is a combination of multiple instructions. When running in the processor 501, it can achieve:
[0185] Obtaining permeability parameters of the water intrusion zone, and using the permeability parameters to construct a water-sealed air volume calculation model;
[0186] Performing data conversion on the water-seal gas volume calculation model to obtain a material balance equation for a water-drive gas reservoir;
[0187] Draw a theoretical map according to the material balance equation of the water-drive gas reservoir, divide the theoretical map into energy regions;
[0188] The water invasion and water sealing strength of the water-containing gas reservoir in the water invasion zone are determined according to the dynamic characteristic parameters in the energy region.
[0189] Specifically, the specific implementation method of the processor 501 for the above instructions can refer to the description of the relevant steps in the corresponding embodiment of the accompanying drawings, which will not be repeated here.
[0190] Furthermore, if the module / unit integrated in the electronic device is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. The computer-readable storage medium can be volatile or non-volatile. For example, the computer-readable medium may include: any entity capable of carrying the computer program code, a recording medium, a USB flash drive, a mobile hard disk, a magnetic disk, an optical disk, a computer memory, and a read-only memory (ROM).
[0191] Embodiment 5
[0192] This embodiment provides a storage medium storing a computer program. When the computer program is executed by a processor, the steps of the method for determining water invasion and water seal strength of a water-containing gas reservoir as described above are implemented.
[0193] These program codes can also be loaded onto a computer or other programmable data processing device so that a series of operation steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process in the computer or other programmable device. Figure 1 The steps of a specified function in a process or multiple processes.
[0194] Storage media include permanent and non-permanent, removable and non-removable media, and can be implemented by any method or technology to store information. Information can be computer-readable instructions, data structures, modules of programs or other data. Examples of storage media can include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, magnetic disk storage or other magnetic storage devices or any other non-transmission medium that can be used to store information that can be accessed by a computing device.
[0195] It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit the exemplary embodiments according to the present application. When the terms "include" and / or "include" are used in this specification, they indicate the presence of features, steps, operations, devices, components and / or their combinations.
[0196] It is to be understood that the terms so used are interchangeable under appropriate circumstances such that the embodiments of the application described herein are, for example, capable of operation in sequences other than those illustrated or described herein.
[0197] In the several embodiments provided by the present invention, it should be understood that the disclosed devices, systems and methods can be implemented in other ways. For example, the device embodiments described above are only illustrative, for example, the division of the modules is only a logical function division, and there may be other division methods in actual implementation.
[0198] The modules described as separate components may or may not be physically separated, and the components shown as modules may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
[0199] In addition, each functional module in each embodiment of the present invention may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of hardware plus software functional modules.
[0200] It is obvious to those skilled in the art that the present invention is not limited to the details of the above exemplary embodiments, and that the present invention can be implemented in other specific forms without departing from the spirit or essential characteristics of the present invention.
[0201] Therefore, no matter from which point of view, the embodiments should be regarded as illustrative and non-restrictive, and the scope of the present invention is limited by the appended claims rather than the above description, so it is intended that all changes falling within the meaning and scope of the equivalent elements of the claims are included in the present invention. Any attached figure mark in the claims should not be regarded as limiting the claims involved.
[0202] In addition, it is clear that the word "comprising" does not exclude other units or steps, and the singular does not exclude the plural. Multiple units or devices stated in the system claim can also be implemented by one unit or device through software or hardware. The words first, second, etc. are used to indicate names, and do not indicate any particular order.
[0203] Finally, it should be noted that the above embodiments are only used to illustrate the technical solution of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that the technical solution of the present invention can be modified or replaced by equivalents without departing from the spirit and scope of the technical solution of the present invention.
Claims
1. A method for determining water intrusion and water seal strength in a water-containing gas reservoir, characterized in that: The method comprises: Obtaining permeability parameters of the water intrusion zone, and using the permeability parameters to construct a water-sealed air volume calculation model; Performing data conversion on the water-seal gas volume calculation model to obtain a material balance equation for a water-drive gas reservoir; Draw a theoretical map according to the material balance equation of the water-drive gas reservoir, divide the theoretical map into energy regions; The water invasion and water sealing strength of the water-containing gas reservoir in the water invasion zone are determined according to the dynamic characteristic parameters in the energy region.
2. The method for determining water invasion and water seal strength of a water-containing gas reservoir according to claim 1, characterized in that: The method of constructing a water-sealed air volume calculation model using the permeability parameters includes: Selecting target permeability parameters of a high permeability layer and a low permeability layer from the permeability parameters, calculating the target permeability parameters, and obtaining the water body seepage rate; Obtaining the water-sealed gas volume of the low permeability layer and the gas volume coefficient under the current formation conditions, performing ratio calculation on the target permeability parameter, the gas volume coefficient under the current formation conditions and the water-sealed gas volume to obtain the pore volume ratio; Convert the expression of the pore volume ratio to obtain the expression of the water-sealed gas volume; Obtaining the produced water volume of the gas reservoir, adding and calculating the target permeability parameter and the produced water volume of the gas reservoir to obtain the water invasion volume of the gas reservoir; Setting the permeability ratio of the high permeability layer and the low permeability layer, combining the water body seepage, the pore volume ratio, the expression of the water-sealed gas volume, the gas reservoir water intrusion and the permeability ratio to obtain a target equation for the water-sealed gas volume; Obtaining the gas volume coefficient, the dynamic reserves of the gas reservoir, and the water storage volume coefficient of the gas reservoir under the original formation conditions, and calculating the gas reservoir equation corresponding to the water storage volume coefficient of the gas reservoir according to the gas volume coefficient, the dynamic reserves of the gas reservoir, and the water storage volume coefficient of the gas reservoir under the original formation conditions; A water-sealed gas volume calculation model is generated according to the target equation for the water-sealed gas volume and the gas reservoir equation.
3. The method for determining water invasion and water seal strength of a water-containing gas reservoir according to claim 2, characterized in that: The step of generating a water-seal gas volume calculation model according to the target equation for the water-seal gas volume and the gas reservoir equation comprises: According to the target equation of the water-sealed gas volume, a reservoir characteristic constant is set, and the reservoir characteristic constant, the target equation and the gas reservoir equation are combined to obtain a water-sealed gas volume calculation model; The reservoir characteristic constant is expressed as: Wherein, A represents the reservoir characteristic constant, V represents the pore volume ratio, and K represents the permeability ratio; The water seal gas volume calculation model is expressed as: Among them, G F represents the water-sealed gas volume, G represents the dynamic reserves of the gas reservoir, and B g0 represents the gas volume coefficient under the original formation conditions, B g represents the gas volume coefficient under the current formation conditions, A represents the reservoir characteristic constant, and t represents the water volume coefficient of the gas reservoir.
4. The method for determining water invasion and water seal strength of a water-containing gas reservoir according to claim 1, characterized in that: The data conversion of the water-seal gas volume calculation model is performed to obtain the material balance equation of the water-drive gas reservoir, including: Obtaining a gas volume coefficient in the water-seal gas volume calculation model, and updating the gas volume coefficient according to a preset formation pressure to obtain an updated gas volume coefficient; A dimensionless pseudo-pressure is generated according to the updated gas volume coefficient, and a material balance equation of a water-drive gas reservoir is generated according to the dimensionless pseudo-pressure and the water-seal gas volume calculation model.
5. The method for determining water invasion and water seal strength of a water-containing gas reservoir according to claim 4, characterized in that: The method of generating a material balance equation for a water-drive gas reservoir according to the dimensionless pseudo-pressure and the water-seal gas volume calculation model comprises: Obtaining water invasion characteristic constants, current cumulative gas production and gas reservoir production degree, extracting gas reservoir water storage volume coefficient and gas reservoir dynamic reserves in the water-sealed gas volume calculation model, and establishing an equation based on the water invasion characteristic constants, the current cumulative gas production, the gas reservoir production degree, the gas reservoir water storage volume coefficient and the gas reservoir dynamic reserves; A material balance equation for a water-driven gas reservoir is established based on the equation, the dimensionless pseudo-pressure and the water-seal gas volume calculation model.
6. The method for determining water invasion and water seal strength of a water-containing gas reservoir according to claim 1, characterized in that: The method of drawing a theoretical diagram according to the material balance equation of the water drive gas reservoir comprises: Performing parameter analysis on the material balance equation of the water drive gas reservoir to obtain relevant parameters; The dimensionless pseudo-pressure and gas reservoir recovery degree in the relevant parameters are obtained, and based on the material balance equation of the water-driven gas reservoir, a theoretical chart is generated by using a preset production decline analysis method with the dimensionless pseudo-pressure as a column vector and the gas reservoir recovery degree as a row vector.
7. The method for determining water invasion and water seal strength of a water-containing gas reservoir according to claim 1, characterized in that: The theoretical diagram is divided to obtain energy regions, including: Extracting a reservoir characteristic constant and a water invasion characteristic constant from the material balance equation of the water drive gas reservoir, and performing a numerical judgment on the reservoir characteristic constant when the water invasion characteristic constant is not 0; When the reservoir characteristic constant is 0, the theoretical map is divided into an energy replenishment zone; When the reservoir characteristic constant is greater than 0 and less than 1, the theoretical map is divided into an energy replenishment dominant area; When the reservoir characteristic constant is 1, the theoretical map is divided on the equilibrium line; When the reservoir characteristic constant is greater than 1, the theoretical map is divided into an energy isolation dominant area; The energy replenishment area, the energy replenishment leading area, the balance line and the energy isolation leading area together constitute an energy area.
8. A water invasion and water seal strength determination system for water-containing gas reservoirs, characterized in that: The system comprises: A calculation model building module, used to obtain the permeability parameters of the water intrusion area and build a water-sealed air volume calculation model using the permeability parameters; A data conversion module, used for performing data conversion on the water-seal gas volume calculation model to obtain a material balance equation for a water-drive gas reservoir; An energy region generation module is used to draw a theoretical map according to the material balance equation of the water-drive gas reservoir, divide the theoretical map, and obtain energy regions; The strength determination module is used to determine the water invasion and water sealing strength of the water-containing gas reservoir in the water invasion area according to the dynamic characteristic parameters in the energy area.
9. An electronic device, comprising: processor; a memory for storing instructions executable by the processor; The processor is configured to execute the instructions to implement the method for determining water invasion and water seal strength in a water-containing gas reservoir as described in any one of claims 1 to 7.
10. A computer-readable storage medium, characterized in that: A computer program is stored thereon, and when the program is executed by a processor, the method for determining water invasion and water seal strength of a water-containing gas reservoir as described in any one of claims 1 to 7 is implemented.