Method and device for calculating CO2 structure burying and storage amount of depleted condensate gas reservoir, electronic equipment, storage medium and computer product

Abstract

The invention relates to the technical field of CO2 burying and storing amount calculation, in particular to an exhausted condensate gas reservoir CO2 structure burying and storing amount calculation method and device, electronic equipment, a storage medium and a computer product. The method comprises the following steps: constructing a depleted condensate gas reservoir CO2 buried substance balance equation; deducing the depleted condensate gas reservoir CO2 burying material balance equation, and constructing a depleted condensate gas reservoir CO2 structure burying amount calculation model; based on the depleted condensate gas reservoir CO2 structure burying stock calculation model, the depleted condensate gas reservoir CO2 structure burying stock is obtained. According to the method, the change of a mixed gas deviation factor after CO2 injection, the change of the saturation degree of residual stratum condensate oil in the CO2 injection process and the influence of deformation of rocks and bound water are considered. The calculation result of an example shows that the CO2 structure burying amount under different burying pressures in the new method is basically matched with the calculation result of a numerical simulation method, and the error is kept within 6.42%.

Description

Technical Field

[0001] This invention relates to the field of CO2 storage capacity calculation technology, and in particular to a method, apparatus, electronic device, storage medium and computer product for calculating CO2 structural storage capacity in depleted condensate gas reservoirs. Background Technology

[0002] Depleted condensate gas reservoirs are natural sites for underground CO2 storage. Accurate and rapid assessment of structural storage volume is fundamental to evaluating CO2 storage potential and is crucial for quantifying the environmental, economic, and social benefits of CO2 storage. Current methods for evaluating CO2 storage potential mainly include the area method, volume method, capacity coefficient method, analogy method, solubility method, mass balance method, and numerical simulation method. Among these, the area method, volume method, and analogy method are relatively simple in calculation, but their accuracy depends on the value of the storage correction coefficient. Furthermore, they do not adequately consider CO2 storage mechanisms, changes in rock fluid properties, and reservoir heterogeneity, often leading to results that deviate from reality. Numerical simulation is an effective means of assessing the CO2 storage potential of various reservoirs; however, high-precision numerical modeling is difficult, and the processes of production history fitting and numerical simulation are time-consuming, posing challenges for practical applications.

[0003] A literature review on CO2 storage calculation methods revealed that the mass balance method has broad application prospects. For example, the paper "A Method for Determining the CO2 Dissolution and Storage Potential of Saline Aquifers" published by China University of Petroleum (East China) solves the problem of CO2 storage in aquifers. The paper "A Method for Evaluating the CO2 Storage Potential of Gas Reservoirs Based on Mass Balance" published by Southwest Petroleum University establishes a method for calculating CO2 storage in dry gas reservoirs. Chinese patent CN115034489B discloses a method for predicting the CO2 storage potential of gas reservoirs considering dissolution. This patent improves upon traditional dry gas reservoir storage calculation methods by considering the effect of gas dissolution, and its research is also based on dry gas reservoirs. Chinese patent CN116502756A discloses a rapid prediction method for the theoretical CO2 storage in constant-volume closed oil reservoirs based on well testing theory. Existing CO2 storage calculation methods mainly focus on saline aquifers, oil reservoirs, and dry gas reservoirs. None of the aforementioned papers and patents can accurately calculate the structural CO2 storage in condensate gas reservoirs.

[0004] CO2 capture, utilization, and storage (CCUS) technology is considered a safety net for achieving carbon neutrality. Currently, coal seams, brine aquifers, depleted oil and gas reservoirs, and the deep sea are the main geological sites for CO2 storage. Compared to coal seams, brine aquifers, and the deep sea, depleted oil and gas reservoirs have a higher level of exploration and development, better geological understanding, and reduced leakage risk after CO2 injection. Furthermore, the existing injection and production well networks and surface pipelines within the reservoirs can effectively reduce the investment cost of CO2 geological storage. However, current CO2 storage projects are mostly concentrated in aquifers and depleted oil reservoirs, while research on CO2 storage in condensate gas reservoirs is still in the exploratory stage. During CO2 storage in condensate gas reservoirs, the gas-liquid balance of formation fluids changes. Conventional methods for calculating CO2 structural storage in dry gas reservoirs lack consideration for changes in reverse condensation and the compressibility factor after gas mixing, often leading to significant errors in calculation results in practical field applications. Therefore, how to scientifically and accurately evaluate the CO2 storage capacity of such condensate gas reservoirs has become a pressing technical challenge for the industry. Summary of the Invention

[0005] To address the aforementioned problems, this invention provides a method, apparatus, electronic device, storage medium, and computer product for calculating the CO2 structural reserves of depleted condensate gas reservoirs, thereby improving the prediction efficiency and calculation accuracy of CO2 structural reserves in such condensate gas reservoirs. The invention provides the following technical solution:

[0006] In a first aspect of the present invention, a method for calculating the structurally buried CO2 reserves of depleted condensate gas reservoirs is provided, the method comprising:

[0007] Based on the amount of original condensate gas in the depleted condensate gas reservoir, the amount of CO2 injected gas during the reservoir development stage or CO2 storage stage, the amount of produced condensate gas, the amount of produced condensate oil, the amount of currently remaining condensate gas, and the amount of currently remaining condensate oil, a CO2 storage material balance equation for the depleted condensate gas reservoir is constructed.

[0008] The material balance equation for CO2 burial in the depleted condensate gas reservoir is derived, and a calculation model for the structural burial amount of CO2 in the depleted condensate gas reservoir is constructed.

[0009] Based on the CO2 structural burial volume calculation model of the depleted condensate gas reservoir, the CO2 structural burial volume of the depleted condensate gas reservoir is obtained.

[0010] Specifically, the mass balance equation for CO2 burial in the depleted condensate gas reservoir is expressed as follows:

[0011]

[0012] In the formula, n ig The amount of original condensate gas in a depleted condensate gas reservoir; The amount of gaseous material injected during the gas reservoir development or CO2 storage phase; n pg The amount of condensate produced; n po The amount of condensate oil produced; n ro n represents the amount of remaining condensate oil. rg This represents the amount of remaining condensate gas.

[0013] Specifically, the calculation model for CO2 structural burial volume in depleted condensate gas reservoirs includes a calculation model for the CO2 structural burial volume in depleted condensate gas reservoirs and / or a calculation model for the CO2 structural burial mass in depleted condensate gas reservoirs.

[0014] Specifically, the calculation model for the structural burial volume of CO2 in depleted condensate gas reservoirs is expressed as follows:

[0015]

[0016] In the formula, For the CO2 structural burial volume of a depleted condensate gas reservoir; G p V represents the volume of condensate gas extracted from a depleted condensate gas reservoir. po The amount of condensate oil extracted; ρ o The density of the extracted condensate oil; M o V represents the molecular weight of condensate oil extracted under surface conditions. HC S represents the pore volume of hydrocarbon fluids in the gas reservoir. oc ρ represents the saturation of the remaining condensate oil under the current formation conditions. oc M represents the current density of the remaining condensate oil. oc The molecular weight of the remaining condensate oil; p i Z represents the original formation pressure; gi is the condensate gas deviation factor under the original formation conditions; R is the universal gas constant; T is the formation temperature; p c S represents the current formation pressure. w C represents the water saturation of a depleted condensate gas reservoir. w Z is the formation water compressibility coefficient; c This represents the deviation factor for the mixture of condensate gas and CO2 under current formation conditions; C f Δp is the compressibility coefficient of the formation rock; Δp is the pressure difference; p sc Pressure under standard conditions; Z sc T is the gas deviation factor under standard conditions. sc Temperature under standard conditions.

[0017] Specifically, the calculation model for the structural burial mass of CO2 in depleted condensate gas reservoirs is expressed as follows:

[0018]

[0019] In the formula, For the CO2 structural burial quality of depleted condensate gas reservoirs; G p The volume of condensate gas extracted from the depleted condensate gas reservoir; V is the density of CO2 under standard conditions. po The amount of condensate oil extracted; ρ o The density of the extracted condensate oil; M o V represents the molecular weight of condensate oil extracted under surface conditions. HC S represents the pore volume of hydrocarbon fluids in a depleted condensate gas reservoir. oc ρ represents the saturation of the remaining condensate oil under the current formation conditions. oc M represents the current density of the remaining condensate oil. oc The molecular weight of the remaining condensate oil; p i Z represents the original formation pressure; gi is the condensate gas deviation factor under the original formation conditions; R is the universal gas constant; T is the formation temperature; p c S represents the current formation pressure. w C represents the water saturation of a depleted condensate gas reservoir. w Z is the formation water compressibility coefficient; c This represents the deviation factor for the mixture of condensate gas and CO2 under current formation conditions; C f Δp is the compressibility coefficient of the formation rock; Δp is the pressure difference; p sc Pressure under standard conditions; Z sc T is the gas deviation factor under standard conditions. sc Temperature under standard conditions.

[0020] Specifically, the deviation factor of the mixture of condensate gas and CO2 under the current formation conditions is obtained by calculation using the equation of state;

[0021] The saturation of the remaining condensate oil under the aforementioned formation conditions was obtained through a PVT test.

[0022] In a second aspect of the invention, a device for calculating the structural burial volume of CO2 in a depleted condensate gas reservoir is provided, the device comprising:

[0023] The first building unit is to construct a CO2 burial material balance equation for a depleted condensate gas reservoir based on the amount of original condensate gas material, the amount of CO2 injected gas material during the reservoir development stage or CO2 burial stage, the amount of produced condensate gas material, the amount of produced condensate oil material, the amount of currently remaining condensate gas material, and the amount of currently remaining condensate oil material.

[0024] The second construction unit is used to derive the material balance equation for CO2 burial in the depleted condensate gas reservoir and construct a calculation model for the structural burial amount of CO2 in the depleted condensate gas reservoir.

[0025] The acquisition unit is used to acquire the CO2 structural burial volume of the depleted condensate gas reservoir based on the CO2 structural burial volume calculation model of the depleted condensate gas reservoir.

[0026] Specifically, the calculation model for the structural burial volume of CO2 in depleted condensate gas reservoirs is expressed as follows:

[0027]

[0028] In the formula, For the CO2 structural burial volume of a depleted condensate gas reservoir; G p The volume of condensate gas extracted from a depleted condensate gas reservoir; V po The amount of condensate oil extracted; ρ o The density of the extracted condensate oil; M o V represents the molecular weight of condensate oil extracted under surface conditions. HC S represents the pore volume of hydrocarbon fluids in a depleted condensate gas reservoir. oc ρ represents the saturation of the remaining condensate oil under the current formation conditions. oc M represents the current density of the remaining condensate oil. oc The molecular weight of the remaining condensate oil; p i Z represents the original formation pressure; gi is the condensate gas deviation factor under the original formation conditions; R is the universal gas constant; T is the formation temperature; p c S represents the current formation pressure. w C represents the water saturation of a depleted condensate gas reservoir. w Z is the formation water compressibility coefficient; c This represents the deviation factor for the mixture of condensate gas and CO2 under current formation conditions; C f Δp is the compressibility coefficient of the formation rock; Δp is the pressure difference; p sc Pressure under standard conditions; Z sc T is the gas deviation factor under standard conditions. sc Temperature under standard conditions.

[0029] Specifically, the calculation model for the structural burial volume of CO2 in depleted condensate gas reservoirs is also expressed as:

[0030]

[0031] In the formula, For the CO2 structural burial quality of depleted condensate gas reservoirs; G p The volume of condensate gas extracted from the depleted condensate gas reservoir; V is the density of CO2 under standard conditions. po The amount of condensate oil extracted; ρ o The density of the extracted condensate oil; M o V represents the molecular weight of condensate oil extracted under surface conditions. HCS represents the pore volume of hydrocarbon fluids in a depleted condensate gas reservoir. oc ρ represents the saturation of the remaining condensate oil under the current formation conditions. oc M represents the current density of the remaining condensate oil. oc The molecular weight of the remaining condensate oil; p i Z represents the original formation pressure; gi is the condensate gas deviation factor under the original formation conditions; R is the universal gas constant; T is the formation temperature; p c S represents the current formation pressure. w C represents the water saturation of a depleted condensate gas reservoir. w Z is the formation water compressibility coefficient; c This represents the deviation factor for the mixture of condensate gas and CO2 under current formation conditions; C f Δp is the compressibility coefficient of the formation rock; Δp is the pressure difference; p sc Pressure under standard conditions; Z sc T is the gas deviation factor under standard conditions. sc Temperature under standard conditions.

[0032] In a third aspect of the invention, an electronic device is provided, the electronic device comprising at least one processor and at least one memory, the memory being data-connected to the processor, wherein...

[0033] The memory stores instructions that can be executed by the at least one processor, which, when executed by the at least one processor, enables the at least one processor to perform the method described above.

[0034] In a fourth aspect of the invention, a computer-storeable medium is provided, characterized in that the storage medium stores computer instructions, which, when executed by a processor, specifically perform the steps in the method described above.

[0035] In a fifth aspect of the invention, a computer program product is provided, comprising computer instructions, characterized in that, when the computer instructions are executed by a processor, they specifically perform the steps in the method described above.

[0036] The technical effects and advantages of this invention are as follows:

[0037] The proposed method for calculating CO2 structural burial reserves in condensate gas reservoirs comprehensively considers numerous influencing factors during the CO2 burial process, which is of great significance for optimizing CO2 structural burial reserve calculations. Compared to traditional mass balance methods, this method focuses on the changes in the mixed gas deviation factor after CO2 injection, the changes in the saturation of residual condensate oil in the formation during CO2 injection, and the influence of rock and bound water deformation. Example calculations show that the CO2 structural burial reserves calculated by the new method under different burial pressures are in good agreement with those calculated by numerical simulation methods, with the error remaining within 6.42%.

[0038] The method of this invention can accurately and reliably assess the CO2 burial quantity in condensate gas reservoirs at any pressure, making it suitable for the effective assessment of CO2 structural burial quantities in condensate gas reservoirs in mining areas. Furthermore, ignoring variations in the mixed gas deviation factor leads to an overestimation of the calculated CO2 structural burial quantity, with an error reaching 18.61%. Ignoring variations in reverse condensation results in an underestimation of the calculated CO2 structural burial quantity, with an error as high as 9.93% in the low burial pressure range. The CO2 structural burial quantity calculation method that considers reverse condensation variations and corrects for the mixed gas deviation factor yields more accurate prediction results.

[0039] Other features and advantages of the invention will be set forth in the following description, and will be apparent in part from the description, or may be learned by practicing the invention. The objects and other advantages of the invention can be realized and obtained by means of the structures pointed out in the description and the drawings. Attached Figure Description

[0040] Figure 1 This is a flowchart of the method for calculating the structural burial volume of CO2 in depleted condensate gas reservoirs provided in the embodiments of this application;

[0041] Figure 2 This is a schematic diagram of the material balance method for CO2 burial process provided in an embodiment of this application;

[0042] Figure 3 This is a diagram of a device for calculating the structural burial volume of CO2 in a depleted condensate gas reservoir, provided in an embodiment of this application.

[0043] Figure 4 This is a block diagram of the electronic device structure according to an embodiment of this application;

[0044] Figure 5 A diagram showing the condensate saturation after CO2 injection into the remaining formation fluids, provided in an embodiment of this application.

[0045] Figure 6 A comparison chart of CO2 tectonic burial volume and numerical simulation calculation results provided in the embodiments of this application;

[0046] Figure 7A comparison chart of the calculation results of the burial volume before and after considering the change of the mixed gas deviation factor is provided for the embodiments of this application;

[0047] Figure 8 A comparison chart of the calculated CO2 tectonic stock volume before and after considering the changes in anti-condensation, provided for embodiments of this application. Detailed Implementation

[0048] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0049] To address the shortcomings of existing technologies, this invention discloses a method for calculating the structural burial volume of CO2 in depleted condensate gas reservoirs, such as... Figure 1 As shown, the method includes,

[0050] Step 1: Based on the amount of original condensate gas in the depleted condensate gas reservoir, the amount of CO2 injected gas during the reservoir development stage or CO2 storage stage, the amount of produced condensate gas, the amount of produced condensate oil, the amount of currently remaining condensate gas, and the amount of currently remaining condensate oil, construct the CO2 storage material balance equation for the depleted condensate gas reservoir.

[0051] Step 2: Derive the material balance equation for CO2 storage in the depleted condensate gas reservoir and construct a calculation model for the structural storage of CO2 in the depleted condensate gas reservoir;

[0052] Step 3: Based on the CO2 structural burial volume calculation model of the depleted condensate gas reservoir, obtain the CO2 structural burial volume of the depleted condensate gas reservoir.

[0053] In one specific embodiment of the present invention, for step 1, combined with Figure 2 As shown in the schematic diagram of the CO2 burial process material balance method of this invention, the original material of the condensate gas reservoir is mainly condensate gas. With the depressurization development of the reservoir, condensate oil is released, at which point the reservoir material mainly consists of condensate gas and condensate oil. During the burial process, CO2 gas is injected, and during this process, some condensate oil evaporates and transforms into the gas phase. Therefore, based on the amount of original material in the depleted condensate gas reservoir, the amount of CO2 injected, the amount of produced material, and the amount of remaining material in the depleted condensate gas reservoir, a material balance equation for CO2 burial of the depleted condensate gas reservoir is constructed, expressed as:

[0054]

[0055] In the formula, n igThe amount of original condensate gas in the depleted condensate gas reservoir, in kmol; The amount of gaseous substance injected during the gas reservoir development or CO2 storage phase, expressed in kmol; n pg The amount of condensate produced, in kmol; n po The amount of condensate produced is expressed in kmol; n ro n represents the current amount of remaining condensate oil, in kmol; rg The amount of remaining condensate gas, in kmol.

[0056] In a specific embodiment of the present invention, for step 2: further deriving the material balance equation for CO2 burial in the depleted condensate gas reservoir and constructing a calculation model for the structural burial volume of CO2 in the depleted condensate gas reservoir, including:

[0057] The amount of condensate gas under the original formation pressure can be expressed as:

[0058]

[0059] In the formula: p i The original formation pressure is MPa; V HC The pore volume of hydrocarbon fluids in the gas reservoir is m. 3 Z gi R is the condensate gas deviation factor under the original formation conditions; T is the general gas constant; and T is the formation temperature in K.

[0060] The amount of condensate gas separated from the surface is:

[0061]

[0062] Where: G p m represents the amount of natural gas produced from a condensate gas reservoir. 3 .

[0063] The amount of oil separated from the extracted condensate gas is:

[0064]

[0065] In the formula: V po To determine the amount of condensate oil extracted, m 3 ;ρ o To determine the density of the extracted condensate oil, kg / m³ 3 M o The molecular weight of the condensate oil produced under surface conditions is given in kg / kmol.

[0066] Amount of CO2 gas injected:

[0067]

[0068] In the formula: The amount of CO2 injected into the gas reservoir, m 3 .

[0069] The amount of remaining mixed gas substance n rg for:

[0070]

[0071] In the formula: p c Z represents the current formation pressure, in MPa; c C is the gas mixture deviation factor under current formation conditions. w C is the formation water compressibility coefficient. f ρ is the compressibility coefficient of the formation rock; Δp is the pressure difference, MPa. oc This represents the saturation of condensate oil under the current formation conditions.

[0072] The amount of residual condensate oil n ro for:

[0073]

[0074] In the formula: S oc ρ represents the saturation of condensate oil under the current formation conditions. oc Density of residual condensate oil, kg / m³ 3 M oc The value represents the molecular weight of the remaining condensate oil liquid, in kg / kmol.

[0075] Substituting equations (2) to (7) into equation (1) and rearranging, we obtain the following material balance equation for CO2 burial in condensate gas reservoirs, considering condensate gas anti-condensation and rock deformation:

[0076]

[0077] eliminate The equation for calculating the structural burial volume of CO2 in a depleted condensate gas reservoir, derived from the deformation, is as follows:

[0078]

[0079] According to formula (9), the equation for calculating the CO2 structural burial mass of depleted condensate gas reservoirs can be obtained as follows:

[0080]

[0081] In the formula, For the CO2 structural burial quality of depleted condensate gas reservoirs; G p The amount of natural gas extracted from a depleted condensate gas reservoir; V is the density of CO2 under standard conditions. po The amount of condensate oil extracted; ρo The density of the extracted condensate oil; M o To determine the molecular weight of condensate oil under surface conditions; V HC S represents the pore volume of hydrocarbon fluids in the gas reservoir. oc ρ represents the saturation of condensate oil under the current formation conditions. oc M represents the density of the remaining condensate oil. oc The molecular weight of the remaining condensate oil liquid; p i V represents the original formation pressure. Hc Z represents the pore volume of hydrocarbon fluids in the gas reservoir. gi is the condensate gas deviation factor under the original formation conditions; R is the universal gas constant; T is the formation temperature; p c S represents the current formation pressure. w C represents the water saturation of a depleted condensate gas reservoir. w Z is the formation water compressibility coefficient; c C is the gas mixture deviation factor under current formation conditions. f Δp is the compressibility coefficient of the formation rock; Δp is the pressure difference; p s1 Pressure under standard conditions; Z sc T is the gas deviation factor under standard conditions. sc Temperature under standard conditions.

[0082] Formulas (9) and (10) are formulas for calculating the CO2 structural burial volume of different types of condensate gas reservoirs. Except for Z, the parameters involved in the formulas are... c and S oc In addition, all of these can be obtained directly or indirectly from gas reservoir data.

[0083] In a specific embodiment of the present invention, the unknown parameter Z c and S oc It can be obtained through equation of state calculations and PVT (Pressure-Volume-Temperature) tests.

[0084] The process of calculating the mixed gas deviation factor under the current formation conditions using the equation of state is as follows:

[0085] The specific PR state equation is as follows:

[0086]

[0087] α=[1+m(1-T c 0.5 )] 2 (13)

[0088] m=0.37464+1.542260ω-0.26992ω 2 (14)

[0089] In the formula, P represents pressure; V m ρ is the specific volume; a is the gravitational constant; b is the repulsive constant; α is the temperature correction term; m is the correction term; p c The critical pressure; T c ω is the critical temperature; ω is the eccentricity factor.

[0090] Equation (11) can be transformed into a compression factor representation:

[0091] Z 3 -(1-B)Z 2 +(A-3B 2 -2B)Z-(AB-B 2 -B 3 )=0 (15)

[0092] in:

[0093]

[0094] In the formula, A and B are calculation coefficients; Z represents the deviation factor.

[0095] The saturation of condensate oil under the current formation conditions was obtained through PVT tests, including: conducting relevant PVT tests on the precipitation and evaporation of condensate oil during CO2 burial in condensate gas reservoirs; obtaining the variation of condensate oil saturation with formation pressure after mixing of formation fluids under different CO2 injection ratios; and calculating the condensate oil saturation S under different conditions using burial pressure and other parameters. oc .

[0096] In a second aspect, the present invention also provides a device for calculating the structural burial volume of CO2 in depleted condensate gas reservoirs, such as... Figure 3 As shown, the device includes,

[0097] The first building unit is used to construct the CO2 burial material balance equation of the depleted condensate gas reservoir based on the original material storage, CO2 injection, produced material storage, and current material storage of the depleted condensate gas reservoir.

[0098] The second construction unit is used to derive the CO2 burial material balance equation of the depleted condensate gas reservoir and construct a calculation model for the CO2 structural burial volume of the depleted condensate gas reservoir.

[0099] The acquisition unit is used to acquire the CO2 structural burial volume of the depleted condensate gas reservoir based on the CO2 structural burial volume calculation model of the depleted condensate gas reservoir.

[0100] The calculation model for CO2 structural burial reserves in depleted condensate gas reservoirs is expressed as follows:

[0101]

[0102] In the formula, For the CO2 structural burial quality of depleted condensate gas reservoirs; G p The amount of natural gas extracted from a depleted condensate gas reservoir; V is the density of CO2 under standard conditions. po The amount of condensate oil extracted; ρ o The density of the extracted condensate oil; M o To determine the molecular weight of condensate oil under surface conditions; V HC S represents the pore volume of hydrocarbon fluids in the gas reservoir. oc ρ represents the saturation of condensate oil under the current formation conditions. oc M represents the density of the remaining condensate oil. oc The molecular weight of the remaining condensate oil liquid; p i V represents the original formation pressure. Hc Z represents the pore volume of hydrocarbon fluids in the gas reservoir. gi is the condensate gas deviation factor under the original formation conditions; R is the universal gas constant; T is the formation temperature; p c S represents the current formation pressure. w C represents the water saturation of a depleted condensate gas reservoir. w Z is the formation water compressibility coefficient; c C is the gas mixture deviation factor under current formation conditions. f Δp is the compressibility coefficient of the formation rock; Δp is the pressure difference; p s1 Pressure under standard conditions; Z sc T is the gas deviation factor under standard conditions. sc Temperature under standard conditions.

[0103] The calculation model for CO2 structural burial reserves in depleted condensate gas reservoirs is also expressed as:

[0104]

[0105] In the formula, For the CO2 structural burial volume of a depleted condensate gas reservoir; G p V represents the amount of natural gas already produced from a depleted condensate gas reservoir. po The amount of condensate oil extracted; ρ o The density of the extracted condensate oil; M o To determine the molecular weight of condensate oil under surface conditions; V HC S represents the pore volume of hydrocarbon fluids in the gas reservoir. oc ρ represents the saturation of condensate oil under the current formation conditions. oc M represents the density of the remaining condensate oil. oc The molecular weight of the remaining condensate oil liquid; p i V represents the original formation pressure.Hc Z represents the pore volume of hydrocarbon fluids in the gas reservoir. gi is the condensate gas deviation factor under the original formation conditions; R is the universal gas constant; T is the formation temperature; p c S represents the current formation pressure. w C represents the water saturation of a depleted condensate gas reservoir. w Z is the formation water compressibility coefficient; c C is the gas mixture deviation factor under current formation conditions. f Δp is the compressibility coefficient of the formation rock; Δp is the pressure difference; p sc Pressure under standard conditions; Z sc T is the gas deviation factor under standard conditions. sc Temperature under standard conditions.

[0106] Regarding the apparatus in the above embodiments, the specific manner in which each unit performs its operation has been described in detail in the embodiments related to the method, and will not be elaborated upon here.

[0107] Based on the above disclosure, the present invention also provides an electronic device. For example... Figure 4 As shown, the electronic device of this disclosure includes at least one processor electrically connected to the present invention and at least one memory electrically connected to the processor, wherein the memory stores instructions executable by the at least one processor, the instructions being executed by the at least one processor to enable the at least one processor to perform the method steps as executed by the controller above.

[0108] The readable storage medium provided in this application is a computer-readable storage medium, which stores computer-readable program instructions (i.e., a computer program) for performing the above-described methods. Compared with the prior art, the beneficial effects of the computer-readable storage medium provided in this application are the same as the beneficial effects of the methods provided in the above embodiments, and will not be repeated here.

[0109] This application also provides a computer program product, including a computer program that, when executed by a processor, implements the steps of the method described above. Compared with the prior art, the beneficial effects of the computer program product provided in this application are the same as those of the methods provided in the above embodiments, and will not be repeated here.

[0110] The technical solution of the present invention will be further described below with reference to specific embodiments.

[0111] Example 1:

[0112] The specific steps for calculating the structural burial volume of CO2 in a depleted condensate gas reservoir are as follows:

[0113] Step 1: Obtain reservoir physical property parameters, gas reservoir production dynamic data, and fluid high-pressure physical property parameters.

[0114] In some embodiments, the reservoir physical properties may include parameters indicating the reservoir's physical properties, such as original formation pressure, reservoir temperature, and rock and formation water compressibility coefficients. Production dynamics data may include cumulative production, formation pressure, and other data. High-pressure reservoir fluid physical properties may include parameters such as original gas reservoir fluid properties, deviation factors, and condensate saturation variations.

[0115] Table 1 shows the calculation parameters for the structural burial volume of CO2 in condensate gas reservoirs as follows:

[0116]

[0117]

[0118] Step 2: Determine key calculation parameters

[0119] Unknown parameter Z c and S oc The deviation factors of condensate gas and CO2 under different pressure conditions can be calculated using the PR equation of state and PVT experiments, as detailed below:

[0120] The unknown parameter Z is obtained by calculating according to formula 11-17. c .

[0121] The saturation S of condensate oil under different conditions was obtained through PVT experiments. oc See Figure 5 .

[0122] Step 3: Calculate the CO2 burial volume under different burial pressures using formula (10).

[0123] With a storage pressure of 10 MPa, the calculated CO2 structural storage capacity of the target condensate gas reservoir is: 10,000 tons;

[0124] With a storage pressure of 15 MPa, the calculated CO2 structural storage capacity of the target condensate gas reservoir is: 10,000 tons;

[0125] With a storage pressure of 20 MPa, the calculated CO2 structural storage capacity of the target condensate gas reservoir is: 10,000 tons;

[0126] With a storage pressure of 25 MPa, the calculated CO2 structural storage capacity of the target condensate gas reservoir is: 10,000 tons;

[0127] With a storage pressure of 30 MPa, the calculated CO2 structural storage capacity of the target condensate gas reservoir is: 10,000 tons.

[0128] Comparative Example 1:

[0129] A numerical simulation model of block A was established to simulate the depletion development and CO2 burial process. The model grid step size was 10×20×4, the total number of grids was 5000, the formation dip angle was 15°, a gas well was set at the top of the model, and the CO2 injection rate was 30000 cubic meters / day to obtain the CO2 burial volume under different pressures.

[0130] The CO2 burial volume calculated by the numerical simulation method in Comparative Example 1 and the novel method proposed in this invention under different pressures is as follows: Figure 6 As shown. Obviously, the CO2 burial amount calculated by the method of the present invention is in good agreement with the numerical simulation results across the entire pressure range, with the error remaining within 6.42%, indicating that the model can accurately and reliably assess the CO2 burial amount in condensate gas reservoirs under any pressure.

[0131] Comparative Example 2:

[0132] Example of CO2 tectonic stockpile volume ignoring variations in mixed gas deviation factors:

[0133] According to the CO2 tectonic stock assessment method proposed in this invention, the CO2 tectonic stock is calculated considering the variation of the mixed gas deviation factor (Example 1) and ignoring the variation of the mixed gas deviation factor (Comparative Example 2), respectively. The calculation of ignoring the mixed gas deviation factor is performed according to Formula 10, and the parameter value is Zc = Z CO2 The result is as follows Figure 7 As shown in the figure, compared with considering the variation of the mixed gas deviation factor, the CO2 tectonic burial volume calculated by ignoring the mixed gas deviation factor is too large, with a maximum error of 18.61%.

[0134] Comparative Example 3:

[0135] Example of CO2 tectonic storage capacity ignoring changes in condensate oil:

[0136] According to the CO2 tectonics storage assessment method proposed by this invention, the CO2 tectonics storage is calculated considering changes in anti-condensate oil and neglecting changes in anti-condensate oil, respectively. The CO2 tectonics storage under the condition of neglecting anti-condensate changes is also calculated according to Formula 10. The calculation process lets S... OC =0, the result is as follows Figure 8 As shown in the figure, compared with considering the changes in anti-condensate oil, ignoring the changes in anti-condensate oil results in a larger calculated CO2 structural burial amount in the low burial pressure range, with an error as high as 9.93%. In the high burial pressure range, the calculated results of the two methods are basically the same.

[0137] Finally, it should be noted that the above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art can still modify the technical solutions described in the foregoing embodiments or make equivalent substitutions for some of the technical features. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A method for calculating the structural burial volume of CO2 in depleted condensate gas reservoirs, characterized in that, The method includes, Based on the amount of original condensate gas in the depleted condensate gas reservoir, the amount of CO2 injected gas during the reservoir development stage or CO2 storage stage, the amount of produced condensate gas, the amount of produced condensate oil, the amount of currently remaining condensate gas, and the amount of currently remaining condensate oil, a CO2 storage material balance equation for the depleted condensate gas reservoir is constructed. The material balance equation for CO2 burial in the depleted condensate gas reservoir is derived, and a calculation model for the structural burial amount of CO2 in the depleted condensate gas reservoir is constructed. Based on the CO2 structural burial volume calculation model of the depleted condensate gas reservoir, the CO2 structural burial volume of the depleted condensate gas reservoir is obtained.

2. The method for calculating the structural burial volume of CO2 in depleted condensate gas reservoirs according to claim 1, characterized in that, The mass balance equation for CO2 burial in the depleted condensate gas reservoir is expressed as follows: In the formula, n ig The amount of original condensate gas in a depleted condensate gas reservoir; The amount of gaseous material injected during the gas reservoir development or CO2 storage phase; n pg The amount of condensate produced; n po The amount of condensate oil produced; n ro n represents the amount of remaining condensate oil. rg This represents the amount of remaining condensate gas.

3. The method for calculating the structural burial volume of CO2 in depleted condensate gas reservoirs according to claim 1, characterized in that, The calculation model for CO2 structural burial volume in depleted condensate gas reservoirs includes a calculation model for the CO2 structural burial volume in depleted condensate gas reservoirs and / or a calculation model for the CO2 structural burial mass in depleted condensate gas reservoirs.

4. The method for calculating the structural burial volume of CO2 in depleted condensate gas reservoirs according to claim 3, characterized in that, The calculation model for the structural burial volume of CO2 in a depleted condensate gas reservoir is expressed as follows: In the formula, For the structural burial volume of CO2 in a depleted condensate gas reservoir; g p The volume of condensate gas extracted from a depleted condensate gas reservoir; v po The amount of condensate oil extracted; ρ o The density of the extracted condensate oil; m o The molecular weight of condensate oil extracted under surface conditions; v HC The pore volume of hydrocarbon fluids in the gas reservoir; s oc ρ represents the saturation of the remaining condensate oil under the current formation conditions. oc M represents the current density of the remaining condensate oil. oc The molecular weight of the remaining condensate oil; p i Z represents the original formation pressure; gi is the condensate gas deviation factor under the original formation conditions; R is the universal gas constant; T is the formation temperature; p c The current formation pressure; s w C represents the water saturation of a depleted condensate gas reservoir. w The coefficient of formation water compressibility; Z c This represents the deviation factor for the mixture of condensate gas and CO2 under current formation conditions; C f The coefficient of rock compression in the formation; Δp is the pressure difference; p sc Pressure under standard conditions; Z sc T is the gas deviation factor under standard conditions. sc Temperature under standard conditions.

5. The method for calculating the structural burial volume of CO2 in depleted condensate gas reservoirs according to claim 3, characterized in that, The calculation model for the structural burial mass of CO2 in depleted condensate gas reservoirs is expressed as follows: In the formula, For the CO2 structural burial quality of depleted condensate gas reservoirs; g p The volume of condensate gas extracted from the depleted condensate gas reservoir; v is the density of CO2 under standard conditions; po The amount of condensate oil extracted; ρ o The density of the extracted condensate oil; M o V represents the molecular weight of condensate oil extracted under surface conditions. HC The pore volume of hydrocarbon fluids in a depleted condensate gas reservoir; s oc ρ represents the saturation of the remaining condensate oil under the current formation conditions. oc M represents the current density of the remaining condensate oil. oc The molecular weight of the remaining condensate oil; p i Z represents the original formation pressure; gi is the condensate gas deviation factor under the original formation conditions; R is the universal gas constant; T is the formation temperature; p c S represents the current formation pressure. w C represents the water saturation of a depleted condensate gas reservoir. w The coefficient of formation water compressibility; Z c This represents the deviation factor for the mixture of condensate gas and CO2 under current formation conditions; C f The coefficient of rock compression in the formation; Δp is the pressure difference; p sc Pressure under standard conditions; Z sc T is the gas deviation factor under standard conditions. sc Temperature under standard conditions.

6. The method for calculating the structural burial volume of CO2 in depleted condensate gas reservoirs according to claim 4 or 5, characterized in that, The deviation factor of the mixture of condensate gas and CO2 under the current formation conditions is obtained by calculation through the equation of state; The saturation of the remaining condensate oil under the aforementioned formation conditions was obtained through a PVT test.

7. A device for calculating the structural burial volume of CO2 in a depleted condensate gas reservoir, characterized in that, The device includes, The first building unit is to construct a CO2 burial material balance equation for a depleted condensate gas reservoir based on the amount of original condensate gas material, the amount of CO2 injected gas material during the reservoir development stage or CO2 burial stage, the amount of produced condensate gas material, the amount of produced condensate oil material, the amount of currently remaining condensate gas material, and the amount of currently remaining condensate oil material. The second construction unit is used to derive the material balance equation for CO2 burial in the depleted condensate gas reservoir and construct a calculation model for the structural burial amount of CO2 in the depleted condensate gas reservoir. The acquisition unit is used to acquire the CO2 structural burial volume of the depleted condensate gas reservoir based on the CO2 structural burial volume calculation model of the depleted condensate gas reservoir.

8. The device for calculating the CO2 structural burial volume of depleted condensate gas reservoirs according to claim 7, characterized in that, The calculation model for CO2 structural burial reserves in depleted condensate gas reservoirs is expressed as follows: In the formula, For the CO2 structural burial volume of a depleted condensate gas reservoir; G p The volume of condensate gas extracted from a depleted condensate gas reservoir; V po The amount of condensate oil extracted; ρ o The density of the extracted condensate oil; M o V represents the molecular weight of condensate oil extracted under surface conditions. HC S represents the pore volume of hydrocarbon fluids in a depleted condensate gas reservoir. oc ρ represents the saturation of the remaining condensate oil under the current formation conditions. oc M represents the current density of the remaining condensate oil. oc The molecular weight of the remaining condensate oil; p i Z represents the original formation pressure; gi is the condensate gas deviation factor under the original formation conditions; R is the universal gas constant; T is the formation temperature; p c S represents the current formation pressure. w C represents the water saturation of a depleted condensate gas reservoir. w The coefficient of formation water compressibility; Z c This represents the deviation factor for the mixture of condensate gas and CO2 under current formation conditions; C f The coefficient of rock compression in the formation; Δp is the pressure difference; p sc Pressure under standard conditions; Z sc T is the gas deviation factor under standard conditions. sc Temperature under standard conditions.

9. The device for calculating the CO2 structural burial volume of depleted condensate gas reservoirs according to claim 7, characterized in that, The calculation model for CO2 structural burial reserves in depleted condensate gas reservoirs is also expressed as: In the formula, For the CO2 structural burial quality of depleted condensate gas reservoirs; G p The volume of condensate gas extracted from the depleted condensate gas reservoir; V is the density of CO2 under standard conditions. po The amount of condensate oil extracted; ρ o The density of the extracted condensate oil; M o V represents the molecular weight of condensate oil extracted under surface conditions. HC S represents the pore volume of hydrocarbon fluids in a depleted condensate gas reservoir. oc ρ represents the saturation of the remaining condensate oil under the current formation conditions. oc M represents the current density of the remaining condensate oil. oc The molecular weight of the remaining condensate oil; p i Z represents the original formation pressure; gi is the condensate gas deviation factor under the original formation conditions; R is the universal gas constant; T is the formation temperature; p c S represents the current formation pressure. w C represents the water saturation of a depleted condensate gas reservoir. w The coefficient of formation water compressibility; Z c This represents the deviation factor for the mixture of condensate gas and CO2 under current formation conditions; C f The coefficient of rock compression in the formation; Δp is the pressure difference; p sc Pressure under standard conditions; Z sc T is the gas deviation factor under standard conditions. sc Temperature under standard conditions.

10. An electronic device comprising at least one processor and at least one memory, the memory being data-connected to the processor, wherein, The memory stores instructions that can be executed by the at least one processor to enable the at least one processor to perform the method of any one of claims 1-6.

11. A computer-storable medium, characterized in that, The storable medium stores computer instructions, which, when executed by a processor, specifically perform the steps of the method as described in any one of claims 1-6.

12. A computer program product comprising computer instructions, characterized in that, When the computer instructions are executed by the processor, they specifically perform the steps in the method as described in any one of claims 1-6.

Images