Method, device, equipment and storage medium for calculating storage capacity of aquifer gas storage

Abstract

The present invention discloses a method, device, equipment and storage medium for calculating the storage capacity of an aquifer gas storage reservoir. The method is used for a gas storage reservoir with an open aquifer water body type, and includes the following steps: obtaining geological data and test data of a groundwater layer to determine the ultimate pore volume and ultimate storage capacity of the aquifer-type gas storage reservoir; obtaining water property data of the aquifer and constructing a first functional relationship group for the gas-water interface of the gas storage reservoir changing with aquifer pressure; obtaining a second functional relationship group for the two-phase region of a phase permeability curve changing with displacement cycles; obtaining the gas storage pore volume corresponding to irreducible water saturation and the movable gas pore volume corresponding to movable gas saturation in different displacement cycles; and calculating the storage capacity of the aquifer-type gas storage reservoir under the displacement cycles based on the first and second functional relationship groups and the gas storage pore volume and movable gas pore volume corresponding to the displacement cycles. The present invention improves the accuracy of calculating the storage capacity of the aquifer-type gas storage reservoir.

Description

Technical Field

[0001] The present invention relates to the field of petroleum and natural gas, and in particular to a storage capacity calculation method, device, equipment and storage medium of an aquifer gas storage reservoir. Background Art

[0002] Conventional natural gas storage methods include gaseous storage and liquid storage. Gaseous storage includes ground tank storage, pipeline storage, and underground gas storage (USG). Ground tank storage and pipeline storage can only be used as measures to eliminate the imbalance in gas consumption between day and night. To solve the problem of seasonal gas consumption imbalance, the fundamental solution is to build underground gas storage.

[0003] The storage capacity of a gas storage facility is an important indicator for measuring its peak-shaving capacity. Verifying the accurate storage capacity helps optimize the working gas volume of the gas storage facility and assists in making decisions on the supply and demand balance of the natural gas pipeline network.

[0004] After research, the inventors found that the conventional gas storage capacity determination method in the prior art is not suitable for calculating the storage capacity of a gas storage converted from an underground aquifer, and often has the defect that the calculated gas storage capacity is too large.

[0005] The information disclosed in this background technology section is only intended to enhance understanding of the overall background of the invention and should not be regarded as an admission or any form of suggestion that the information constitutes the prior art already known to a person skilled in the art. Summary of the Invention

[0006] The object of the present invention is to improve the accuracy of calculating the storage capacity of an aquifer-type gas storage reservoir.

[0007] The present invention provides a method for calculating the storage capacity of an aquifer gas storage reservoir, which is used for a gas storage reservoir whose aquifer water body type is an open water body, comprising the steps of:

[0008] S11. Acquire geological data and test data of the groundwater layer to determine the ultimate pore volume and ultimate storage capacity of the aquifer-type gas storage reservoir; the geological data include trap characteristic data and structural characteristic data for determining the ultimate pore volume of the aquifer-type gas storage reservoir; the test data include pressure data, injection and production data, and core flooding test data;

[0009] S12. Acquire water property data of the aquifer, and construct a first functional relationship group of the gas-water interface of the gas storage reservoir and the change of the aquifer pressure based on the water property data; the water property data includes water pressure and water intrusion coefficient;

[0010] S13, obtaining multiple rounds of gas-water mutual displacement experimental data, and obtaining a second functional relationship group of the two-phase region of the phase permeability curve changing with the displacement rounds; the gas-water mutual displacement experimental data includes displacement pressure, displacement time, cumulative gas production, and cumulative water production;

[0011] S14, obtaining the gas storage pore volume corresponding to the irreducible water saturation and the movable gas pore volume corresponding to the movable gas saturation in different displacement cycles;

[0012] S15 calculates the storage capacity of the aquifer-type gas storage reservoir under the displacement cycle based on the first functional relationship group, the second functional relationship group, and the gas storage pore volume and movable gas pore volume corresponding to the displacement cycle; the storage capacity includes dynamic storage capacity and movable gas storage capacity.

[0013] In the present invention, the calculation formula of the limit pore volume includes:

[0014]

[0015] Among them, V max is the limiting pore volume, is the average porosity, A max is the limiting enclosed area, h max is the trap closure height.

[0016] In the present invention, the calculation formula of the limit storage capacity includes:

[0017] G max =V max ·(1-S wc );

[0018] Among them, G max is the limit storage capacity, S wc The irreducible water saturation is obtained from the gas-water displacement experiment in the core flooding experiment.

[0019] In the present invention, the method of obtaining a first functional relationship group of the gas-water interface of the gas storage reservoir changing with the aquifer pressure includes:

[0020] The formula for obtaining the pore volume of the aquifer corresponding to the gas-water interface at time t is: and,

[0021] According to the formula The relationship between water intrusion and aquifer pressure W e =f(P);

[0022] Among them, V t is the pore volume, is the average porosity, A is the gas-bearing area, is the air-water interface at time t, Cs is the water intrusion coefficient, W e is the water intrusion, and f(P) is the functional relationship between the water intrusion and the aquifer pressure.

[0023] In the present invention, the step of obtaining multiple rounds of gas-water mutual displacement experimental data and obtaining a second functional relationship of the two-phase region of the phase permeability curve changing with the displacement rounds includes:

[0024] The gas-water mutual displacement experimental data is used to obtain the gas-water relative permeability curve; wherein, the relative permeability curve processing process of the gas-water displacement process includes:

[0025] S21. Correct the cumulative liquid production value measured at the rock sample outlet pressure to the value at the average rock sample pressure using the following formula:

[0026]

[0027] Among them, V i is the cumulative liquid production value at time i, ΔV wi is the water production from time i-1 to time i, V i-1 is the cumulative liquid production value at time i-1, P a is atmospheric pressure, ΔP is displacement pressure difference, ΔV gi It is the gas increment value measured at a certain time interval under atmospheric pressure;

[0028] S22, correct the cumulative liquid production value, and use the non-steady-state gas-water relative permeability method to calculate the end-face gas saturation S at different times of gas-water displacement. ge , water phase relative permeability K rw and the gas phase relative permeability K rg , the formula is as follows:

[0029]

[0030]

[0031]

[0032] Among them, f w (S g ) is the moisture content, Cumulative water production, is the cumulative liquid production, K rw is the relative permeability of water phase, k rg is the gas phase relative permeability, μ g is the viscosity of natural gas under formation conditions, μ w is the formation water viscosity under formation conditions, I is the flow capacity ratio, Q t is the liquid volume at the core outlet end at time t, Q0 is the water production at the core outlet end at the initial time, Sge is the gas saturation at the core outlet end face;

[0033] The second functional relationship group includes:

[0034] According to the formula S wc =f1(n) is the functional relationship between irreducible water saturation and displacement rounds;

[0035] According to the formula S gr =f2(n) The function relationship between the residual gas saturation and the displacement round number;

[0036] According to the formula S glive =f3(n) The functional relationship between the movable gas saturation and the displacement round number;

[0037] Among them, S wc is the irreducible water saturation, n is the number of displacement cycles, f1(n) is the functional relationship between irreducible water saturation and displacement cycles, S gr is the residual gas saturation, f2(n) is the functional relationship between the residual gas saturation and the displacement rounds, S glive is the movable gas saturation, and f3(n) is the functional relationship between the movable gas saturation and the displacement rounds.

[0038] In the present invention, the gas storage pore volume corresponding to the irreducible water saturation is obtained using the following formula:

[0039] V wc =V t (1-S wc )

[0040] Among them, V wc is the gas storage pore volume of aquifer-type gas storage;

[0041] The movable air pore volume corresponding to the movable water saturation is obtained using the following formula:

[0042] V glive =V t S glive

[0043] Among them, V glive is the movable gas pore volume of aquifer-type gas storage;

[0044] In the present invention, the storage capacity of the aquifer-type gas storage reservoir under the displacement rounds is calculated, including:

[0045] S31. Calculate the dynamic storage capacity G of aquifer-type gas storage g ,include:

[0046]

[0047] Among them, V wc is the gas storage pore volume of aquifer-type gas storage; B g is the volume coefficient of natural gas at formation pressure;

[0048] S32. Calculate the maximum value G of the dynamic library capacity gmax ,include:

[0049] The dynamic storage capacity G g Constrained by the formation upper limit pressure, when the injection and production capacity is expanded to the overflow point, if the formation pressure does not break through the formation upper limit pressure, the dynamic storage capacity G g The maximum value of G gmax The calculation formula includes:

[0050] G gmax =G max

[0051] If the formation pressure exceeds the formation upper limit pressure, the dynamic storage capacity G g The maximum value of G gmax The calculation formula includes:

[0052]

[0053] Among them, V t,limit is the pore volume of the aquifer corresponding to the formation upper limit pressure during the injection and production expansion process, and the calculation formula is:

[0054]

[0055] in, The gas-water interface corresponding to the formation upper limit pressure during the injection-production expansion process;

[0056] S33. Calculate the movable gas storage capacity G of aquifer-type gas storage glive The formulas include:

[0057]

[0058] Among them, B g V is the volume coefficient of natural gas under formation pressure; glive is the movable gas pore volume of the aquifer-type gas storage.

[0059] In another aspect of the present invention, a storage capacity calculation device for an aquifer gas storage reservoir is provided, which is used for a gas storage reservoir whose aquifer water body type is an open water body, comprising:

[0060] a data acquisition unit, configured to acquire geological data and test data of the groundwater layer to determine the ultimate pore volume and ultimate storage capacity of the aquifer-type gas storage reservoir; the geological data includes trap characteristic data and structural characteristic data, which are used to determine the ultimate pore volume of the aquifer-type gas storage reservoir; the test data includes pressure data, injection and production data, and core flooding experiment data;

[0061] a first function acquisition unit configured to acquire water property data of the aquifer and construct a first functional relationship group of changes in the gas-water interface of the gas storage reservoir with the aquifer pressure based on the water property data; the water property data including the water pressure and the water intrusion coefficient;

[0062] A second function acquisition unit is used to acquire multiple rounds of gas-water mutual displacement experimental data and obtain a second function relationship group of the two-phase region of the phase permeability curve changing with the displacement rounds; the gas-water mutual displacement experimental data includes displacement pressure, displacement time, cumulative gas production and cumulative water production;

[0063] A pore volume calculation unit is used to obtain the gas storage pore volume corresponding to the irreducible water saturation in different displacement cycles, and the movable gas pore volume corresponding to the movable gas saturation;

[0064] A storage capacity calculation unit is used to calculate the storage capacity of the aquifer-type gas storage reservoir under the displacement cycle based on the first functional relationship group, the second functional relationship group, and the gas storage pore volume and movable gas pore volume corresponding to the displacement cycle; the storage capacity includes dynamic storage capacity and movable gas storage capacity.

[0065] In another embodiment of the present invention, a gas storage capacity calculation device is provided, including:

[0066] memory for storing computer programs;

[0067] A processor is used to call and execute the computer program to implement each step of the gas storage capacity calculation method as described in any one of the above items.

[0068] In another aspect of the embodiment of the present invention, a storage medium is further provided, on which a computer program is stored. When the computer program is executed by a processor, each step of the gas storage capacity calculation method as described in any one of the above items is implemented.

[0069] The gas storage reservoir capacity calculation device includes a computer program stored on a medium, and the computer program includes program instructions. When the program instructions are executed by a computer, the computer executes the methods described in the above aspects and achieves the same technical effects.

[0070] Compared with the prior art, the present invention has the following beneficial effects:

[0071] It can be seen from the above scheme that the storage capacity calculation method, device, equipment and storage medium for aquifer gas storage provided by the present invention take into account the water type and properties of the underground aquifer and the influence of multi-cycle seepage of the gas storage when calculating the dynamic storage capacity and movable gas storage capacity of the aquifer-type gas storage converted from an open groundwater layer. The natural reduction of the storage capacity caused by periodic injection and production is corrected through a dynamic function relationship, avoiding the defect of the gas storage capacity calculation value being too large that is easy to occur in the prior art. As a result, the accuracy of the calculation results in the present invention is greatly improved compared with the prior art. Moreover, as the injection and production operation cycle of the gas storage is gradually extended, the calculation accuracy will be further improved, which can provide accurate storage capacity data for the construction and operation of aquifer gas storage.

[0072] The above description is only an overview of the technical solution of the present invention. In order to more clearly understand the technical means of the present invention and to implement it according to the contents of the specification, and to make the above and other purposes, technical features and advantages of the present invention easier to understand, one or more preferred embodiments are listed below and described in detail with reference to the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0073] In order to more clearly illustrate the technical solution of the present invention, the following briefly introduces the drawings required for use in the embodiments. 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 any creative work.

[0074] Figure 1 It is a step diagram of the method for calculating the storage capacity of a gas storage reservoir described in the present invention;

[0075] Figure 2 Schematic diagram of the gas-water relative permeability curve of the present invention;

[0076] Figure 3 It is a structural schematic diagram of the gas storage capacity calculation device of the present invention;

[0077] Figure 4 It is a structural diagram of the gas storage capacity calculation device described in the present invention. DETAILED DESCRIPTION

[0078] The specific embodiments of the present invention are described in detail below with reference to the accompanying drawings, but it should be understood that the protection scope of the present invention is not limited by the specific embodiments.

[0079] Unless expressly stated otherwise, throughout the specification and claims, the term “comprise” or variations such as “include” or “comprising” will be understood to include the stated elements or components but not to exclude other elements or components.

[0080] In this document, the terms "first", "second", etc. are used to distinguish two different elements or parts, and are not used to limit specific positions or relative relationships. In other words, in some embodiments, the terms "first", "second", etc. can also be interchangeable with each other.

[0081] Example 1

[0082] In order to improve the accuracy of calculating the storage capacity of aquifer-type gas storage, Figure 1 As shown, in an embodiment of the present invention, a method for calculating the storage capacity of an aquifer gas storage reservoir is provided, which is used for a gas storage reservoir whose aquifer water body type is an open water body, comprising the steps of:

[0083] S11. Acquire geological data and test data of the groundwater layer to determine the ultimate pore volume and ultimate storage capacity of the aquifer-type gas storage reservoir; the geological data include trap characteristic data and structural characteristic data for determining the ultimate pore volume of the aquifer-type gas storage reservoir; the test data include pressure data, injection and production data, and core flooding test data;

[0084] In the embodiment of the present invention, the pressure data is used to determine the water body type and obtain water body property data, and can also be used to determine the first functional relationship group. The core flooding experimental data can be used to construct the second functional relationship group of the phase permeability curve two-phase region changes with the flooding cycle;

[0085] In practical applications, the formula for calculating the limiting pore volume can be as follows based on the aquifer trap area, aquifer trap closure height, and average porosity:

[0086]

[0087] Among them, V max is the limiting pore volume, is the average porosity, A max is the limiting enclosed area, h max is the trap closure height, the trap area is determined based on the trap characteristics, and the structural overflow point is determined based on the structural characteristics. The aquifer trap closure height is determined based on the structural overflow point and the trap characteristics, and the structural overflow point is determined based on the structural characteristics.

[0088] The test data in the embodiment of the present invention may include pressure data, injection and production data, and core flooding test data. Based on the ultimate pore volume of the reconstructed aquifer-type gas storage and the irreducible water saturation obtained from the first round of gas-displacement water experiments in the core flooding experiment, the formula for calculating the ultimate storage capacity of the reconstructed aquifer-type gas storage may be as follows:

[0089] G max =Vmax ·(1-S wc );

[0090] Among them, G max is the limit storage capacity, S wc The irreducible water saturation is obtained from the gas-water displacement experiment in the core flooding experiment.

[0091] S12. Acquire water property data of the aquifer, and construct a first functional relationship group of the gas-water interface of the gas storage reservoir and the change of the aquifer pressure based on the water property data; the water property data includes water pressure and water intrusion coefficient;

[0092] The application scenario of the embodiment of the present invention is to calculate the storage capacity of an aquifer-type gas storage reservoir converted from an open water body groundwater layer. The aquifer water type is an open water body. The water body properties in the embodiment of the present invention may specifically include water body pressure and water intrusion coefficient.

[0093] Water pressure and water intrusion coefficient are the decision parameters that affect the conversion of open water body type into aquifer type gas storage. Water pressure is the original water layer pressure, and water intrusion coefficient C is the initial pressure. s The formula for determination can be as follows:

[0094]

[0095] Among them, C s is the water intrusion coefficient, W e is the water intrusion, P t is the air-water boundary pressure at time t;

[0096] The water intrusion volume is calculated based on the material balance equation of the gas injection expansion process, and the formula is as follows:

[0097] W p B w =W e +(G i -G p )B g

[0098] Among them, W p is the volume of produced water, B w is the formation water coefficient, G i is the cumulative gas injection volume, G p is the cumulative gas production, B g is the volume coefficient of natural gas at formation pressure.

[0099] The functional relationship between the gas-water interface and the aquifer pressure in the embodiment of the present invention (i.e., the first functional relationship group) includes:

[0100] According to the formula Obtain the pore volume of the aquifer corresponding to the gas-water interface at time t;

[0101] According to the formula Obtain the relationship between water intrusion and aquifer pressure W e =f(P);

[0102] Among them, V t is the pore volume, is the average porosity, A is the gas-bearing area, is the air-water interface at time t, W e is the water intrusion, and f(P) is the functional relationship between the water intrusion and the aquifer pressure.

[0103] In an embodiment of the present invention, the gas-water interface is used to delineate the gas-bearing range, calculate the pore volume of the delineated gas-bearing range, and establish a functional relationship between water intrusion and aquifer pressure (i.e., the first functional relationship group) to evaluate the smoothness of the gas-water interface movement during the injection and production process of the gas storage reservoir.

[0104] S13, obtaining multiple rounds of gas-water mutual displacement experimental data, and obtaining a second functional relationship group of the two-phase region of the phase permeability curve changing with the displacement rounds; the gas-water mutual displacement experimental data includes displacement pressure, displacement time, cumulative gas production, and cumulative water production;

[0105] The multiple rounds of gas-water mutual drive experimental data obtained in the embodiment of the present invention may specifically include displacement pressure, displacement time, cumulative gas production, and cumulative water production; the gas-water mutual drive experimental data is used to obtain a gas-water relative permeability curve;

[0106] K rw is the relative permeability of water phase, K rg is the gas phase relative permeability

[0107] The following takes the gas-water displacement process as an example to illustrate the processing of the phase permeability curve:

[0108] S21. Correct the cumulative liquid production value measured at the rock sample outlet pressure to the value at the average rock sample pressure using the following formula:

[0109]

[0110] Among them, V i is the cumulative liquid production value at time i, ΔV wi is the water production from time i-1 to time i, V i-1 is the cumulative liquid production value at time i-1, P a is atmospheric pressure, ΔP is displacement pressure difference, ΔV gi It is the gas increment value measured at a certain time interval under atmospheric pressure;

[0111] S22. Correct the cumulative liquid production value and use the non-steady-state gas-water relative permeability method to calculate the end-face gas saturation, water phase relative permeability, and gas phase relative permeability at different times of gas-water displacement. The formula is as follows:

[0112]

[0113]

[0114] Among them, f w (S g ) is the moisture content, Cumulative water production, is the cumulative liquid production, K rw is the relative permeability of water phase, K rg is the gas phase relative permeability, μ g is the viscosity of natural gas under formation conditions, μ w is the formation water viscosity under formation conditions, I is the flow capacity ratio, Q t is the liquid volume at the core outlet end at time t, Q0 is the water production at the core outlet end at the initial time, S ge is the gas saturation at the core outlet end face;

[0115] When generating displacement phase permeability curves based on multiple rounds (taking 4 rounds as an example) of gas-water mutual displacement experimental data, the obtained gas-water phase permeability curves can be as follows: Figure 2 As shown in the figure, krw-round1, krw-round2, krw-round3 and krw-round4 are used to represent the relative permeability K of the water phase in the four rounds respectively. rw ; krg-round1, krg-round2, krg-round3 and krg-round4 are used to represent the gas phase relative permeability K of the four rounds respectively rg .

[0116] In the present invention, a trap refers to a place where oil and gas can be prevented from continuing to migrate and accumulate. The trap consists of three parts: the reservoir, the cap rock, and the obstruction that prevents the continued migration of oil and gas and causes the accumulation of oil and gas. The residual gas saturation refers to the percentage of the pore volume of the reservoir occupied by the residual gas, expressed in Figure 2 For example, the residual gas saturation is the intersection of 1 and the horizontal axis of the gas phase relative permeability curvature. Irreducible water saturation refers to the percentage of the pore volume of the reservoir occupied by residual water, expressed in Figure 2 For example, the irreducible water saturation is the intersection of the water phase relative permeability curve and the abscissa.

[0117] In an embodiment of the present invention, multiple rounds of gas-water mutual displacement experiments can be used to observe the distribution and flow patterns of gas-water two-phase fluids inside the core, and to construct a functional relationship (i.e., a second functional relationship group) between irreducible water saturation, residual gas saturation, and recoverable gas saturation and displacement cycles.

[0118] The functional relationship of the two-phase region of the phase permeability curve changing with the displacement cycle (i.e., the second functional relationship group) includes:

[0119] According to the formula S wc =f1(n) to obtain the functional relationship between irreducible water saturation and displacement rounds;

[0120] According to the formula S gr =f2(n) to obtain the functional relationship between residual gas saturation and displacement rounds;

[0121] According to the formula S glive =f3(n) to obtain the functional relationship between movable gas saturation and displacement rounds;

[0122] Among them, S wc is the irreducible water saturation, n is the number of displacement cycles, f1(n) is the functional relationship between irreducible water saturation and displacement cycles, S gr is the residual gas saturation, f2(n) is the functional relationship between the residual gas saturation and the displacement rounds, S gilve is the movable gas saturation, and f3(n) is the functional relationship between the movable gas saturation and the displacement rounds.

[0123] In each round of gas-water mutual displacement experiment, the bound water saturation S wc and residual gas saturation S gr are all constant values, but the relative permeability K of the water phase in different rounds of gas-water mutual flooding experiments is different. rw and phase relative permeability K rg will change. Specifically, as the number of rounds increases, the bound water saturation S wc and / or residual gas saturation S gr Will gradually increase.

[0124] S14, obtaining the gas storage pore volume corresponding to the irreducible water saturation and the movable gas pore volume corresponding to the movable gas saturation in different displacement cycles;

[0125] In the embodiment of the present invention, the calculation formula for the gas storage pore volume of the aquifer-type gas storage corresponding to the irreducible water saturation can be as follows:

[0126] V wc =V t (1-S wc );

[0127] Among them, Vwc is the gas storage pore volume of aquifer-type gas storage;

[0128] The calculation formula for the movable pore volume of an aquifer-type gas storage reservoir corresponding to the movable water saturation can be as follows:

[0129] V glive =V t S glive ; Among them, V glive is the movable gas pore volume of the aquifer-type gas storage.

[0130] S15. Calculate the storage capacity of the aquifer-type gas storage reservoir under the displacement cycle based on the first functional relationship group, the second functional relationship group, and the gas storage pore volume and movable gas pore volume corresponding to the displacement cycle; the storage capacity includes the dynamic storage capacity and the movable gas storage capacity;

[0131] S31. Calculate the dynamic storage capacity G of aquifer-type gas storage g ,include:

[0132] The calculation formula for the dynamic storage capacity of aquifer-type gas storage can be as follows:

[0133]

[0134] Among them, G g is the dynamic storage capacity of aquifer-type gas storage; V wc is the gas storage pore volume of aquifer-type gas storage; B g is the volume coefficient of natural gas at formation pressure;

[0135] S32. Calculate the maximum value G of the dynamic library capacity gmax ,include:

[0136] In the embodiment of the present invention, the dynamic reservoir capacity is restricted by the formation upper limit pressure. When the injection and production capacity is expanded to the overflow point, if the formation pressure does not exceed the formation upper limit pressure, the maximum value of the dynamic reservoir capacity G gmax for:

[0137] G gmax =G max

[0138] If the formation pressure exceeds the formation upper limit pressure, the maximum value of the dynamic reservoir capacity G gmax for:

[0139]

[0140] Among them, V t,limit The pore volume of the aquifer corresponding to the formation upper limit pressure during the injection and production expansion process is calculated as follows:

[0141]

[0142] in, The gas-water interface corresponding to the formation upper limit pressure during the injection-production expansion process;

[0143] S33. Calculate the movable gas storage capacity G of aquifer-type gas storage glive The formulas include:

[0144] The movable gas storage capacity G of the aquifer-type gas storage reservoir glive The calculation is as follows:

[0145]

[0146] Among them, B g V is the volume coefficient of natural gas under formation pressure; glive is the movable gas pore volume of the aquifer-type gas storage.

[0147] In summary, the storage capacity calculation method for an aquifer gas storage reservoir provided in an embodiment of the present invention takes into account the water type and properties of the underground aquifer and the influence of multi-cycle seepage of the gas storage reservoir when calculating the dynamic storage capacity and movable gas storage capacity of an aquifer-type gas storage reservoir converted from an open groundwater layer, and corrects the natural reduction in storage capacity caused by periodic injection and production by establishing a dynamic function relationship; avoids the defect of the existing technology that the calculated value of the gas storage reservoir storage capacity is too large, thereby significantly improving the accuracy of the calculation results in the present invention compared with the existing technology, and as the injection and production operation cycle of the gas storage reservoir is successively extended, the calculation accuracy will be further improved, and accurate storage capacity data can be provided for the construction and operation of aquifer gas storage reservoirs.

[0148] Example 2

[0149] Corresponding to the method embodiment, another aspect of the present invention provides a storage capacity calculation device for an aquifer gas storage reservoir. Figure 3 The schematic diagram of the structure of the storage capacity calculation device for aquifer gas storage provided by the embodiment of the present invention is shown. The storage capacity calculation device for aquifer gas storage is Figure 1 The device corresponding to the storage capacity calculation method for aquifer gas storage in the corresponding embodiment, that is, the virtual device is used to realize Figure 1 In the corresponding embodiment of the storage capacity calculation method for an aquifer gas storage reservoir, each virtual module constituting the storage capacity calculation device for an aquifer gas storage reservoir can be executed by an electronic device, such as a network device, a terminal device, or a server. Specifically, the storage capacity calculation device for an aquifer gas storage reservoir in the embodiment of the present invention includes:

[0150] The storage capacity calculation device for an aquifer gas storage reservoir in an embodiment of the present invention is used for a gas storage reservoir whose aquifer water body type is an open water body, and includes:

[0151] The data acquisition unit 01 is used to acquire geological data and test data of the groundwater layer to determine the ultimate pore volume and ultimate storage capacity of the aquifer-type gas storage reservoir. The geological data includes trap characteristic data and structural characteristic data, which are used to determine the ultimate pore volume of the aquifer-type gas storage reservoir. The test data includes pressure data, injection and production data, and core flooding test data.

[0152] The first function acquisition unit 02 is used to acquire water property data of the aquifer and construct a first functional relationship group of the gas-water interface of the gas storage reservoir and the aquifer pressure based on the water property data; the water property data includes water pressure and water intrusion coefficient;

[0153] The second function acquisition unit 03 is used to acquire multiple rounds of gas-water mutual displacement experimental data, and obtain a second functional relationship group of the two-phase region of the phase permeability curve changing with the displacement rounds; the gas-water mutual displacement experimental data includes displacement pressure, displacement time, cumulative gas production and cumulative water production;

[0154] The pore volume calculation unit 04 is used to obtain the gas storage pore volume corresponding to the irreducible water saturation in different displacement cycles, and the movable gas pore volume corresponding to the movable gas saturation;

[0155] The storage capacity calculation unit 05 is used to calculate the storage capacity of the aquifer-type gas storage reservoir under the displacement cycle based on the first functional relationship group, the second functional relationship group, and the gas storage pore volume and movable gas pore volume corresponding to the displacement cycle; the storage capacity includes the dynamic storage capacity and the movable gas storage capacity that is corrected for the natural reduction in storage capacity caused by periodic injection and production by establishing a dynamic functional relationship.

[0156] It should be noted that the specific implementation and technical effects of the storage capacity calculation device for aquifer gas storage in the embodiment of the present invention can be referred to Figure 1 The corresponding storage capacity calculation method for aquifer gas storage will not be described here.

[0157] Example 3

[0158] Corresponding to the method embodiment, the present invention also provides a storage capacity calculation device for an aquifer gas storage reservoir, such as a terminal, a server, etc. The server can be an independent physical server, a server cluster or a distributed system composed of multiple physical servers, or a cloud server that provides basic cloud computing services such as cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communications, middleware services, domain name services, security services, CDN, and big data and artificial intelligence platforms. The terminal can be a smartphone, tablet computer, laptop computer, desktop computer, etc., but is not limited to these.

[0159] An example of a hardware structure block diagram of a storage capacity calculation device for an aquifer gas storage reservoir provided in an embodiment of the present application is shown in FIG. Figure 4 As shown, this may include:

[0160] Processor 1, communication interface 2, memory 3 and communication bus 4;

[0161] The processor 1, the communication interface 2, and the memory 3 communicate with each other via the communication bus 4;

[0162] Optionally, the communication interface 2 may be an interface of a communication module, such as an interface of a GSM module;

[0163] The processor 1 may be a central processing unit (CPU), or an application-specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of the present application.

[0164] The memory 3 may include a high-speed RAM memory, and may also include a non-volatile memory (non-volatile memory), such as at least one disk memory.

[0165] The processor 1 is specifically configured to execute the computer program stored in the memory 3 to perform the following steps:

[0166] The storage capacity calculation method for aquifer gas storage is used for a gas storage with an open aquifer water body type, including the following steps:

[0167] S11. Acquire geological data and test data of the groundwater layer to determine the ultimate pore volume and ultimate storage capacity of the aquifer-type gas storage reservoir; the geological data include trap characteristic data and structural characteristic data for determining the ultimate pore volume of the aquifer-type gas storage reservoir; the test data include pressure data, injection and production data, and core flooding test data;

[0168] S12. Acquire water property data of the aquifer, and construct a first functional relationship group of the gas-water interface of the gas storage reservoir and the change of the aquifer pressure based on the water property data; the water property data includes water pressure and water intrusion coefficient;

[0169] S13, obtaining multiple rounds of gas-water mutual displacement experimental data, and obtaining a second functional relationship group of the two-phase region of the phase permeability curve changing with the displacement rounds; the gas-water mutual displacement experimental data includes displacement pressure, displacement time, cumulative gas production, and cumulative water production;

[0170] S14, obtaining the gas storage pore volume corresponding to the irreducible water saturation and the movable gas pore volume corresponding to the movable gas saturation in different displacement cycles;

[0171] S15 calculates the storage capacity of the aquifer-type gas storage reservoir under the displacement cycle based on the first functional relationship group, the second functional relationship group, and the gas storage pore volume and movable gas pore volume corresponding to the displacement cycle; the storage capacity includes dynamic storage capacity and movable gas storage capacity.

[0172] The above-mentioned product can execute the method provided in the embodiment of the present invention and has the corresponding functional modules and beneficial effects of the execution method. For technical details not fully described in this embodiment, please refer to the storage capacity calculation method for aquifer gas storage provided in the embodiment of the present invention.

[0173] Example 4

[0174] In an embodiment of the present invention, a storage medium is further provided. The storage medium may store a program suitable for execution by a processor, wherein the program is used to:

[0175] A gas storage facility for an open aquifer water body includes the following steps:

[0176] S11. Acquire geological data and test data of the groundwater layer to determine the ultimate pore volume and ultimate storage capacity of the aquifer-type gas storage reservoir; the geological data include trap characteristic data and structural characteristic data for determining the ultimate pore volume of the aquifer-type gas storage reservoir; the test data include pressure data, injection and production data, and core flooding test data;

[0177] S12. Acquire water property data of the aquifer, and construct a first functional relationship group of the gas-water interface of the gas storage reservoir and the change of the aquifer pressure based on the water property data; the water property data includes water pressure and water intrusion coefficient;

[0178] S13, obtaining multiple rounds of gas-water mutual displacement experimental data, and obtaining a second functional relationship group of the two-phase region of the phase permeability curve changing with the displacement rounds; the gas-water mutual displacement experimental data includes displacement pressure, displacement time, cumulative gas production, and cumulative water production;

[0179] S14, obtaining the gas storage pore volume corresponding to the irreducible water saturation and the movable gas pore volume corresponding to the movable gas saturation in different displacement cycles;

[0180] S15. Calculate the storage capacity of the aquifer-type gas storage reservoir under the displacement cycle based on the first functional relationship group, the second functional relationship group, and the gas storage pore volume and movable gas pore volume corresponding to the displacement cycle; the storage capacity includes the dynamic storage capacity and the movable gas storage capacity.

[0181] Optionally, the detailed functions and extended functions of the program may refer to the above description.

[0182] The above-mentioned product can execute the method provided by the embodiment of the present invention, and has the functional modules and beneficial effects corresponding to the execution method. For technical details not fully described in this embodiment, please refer to the methods provided by other embodiments of the present invention.

[0183] The above-mentioned product can execute the method provided by the embodiment of the present invention, and has the functional modules and beneficial effects corresponding to the execution method. For technical details not fully described in this embodiment, please refer to the method provided by the embodiment of the present invention.

[0184] Those skilled in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

[0185] In the several embodiments provided herein, it should be understood that the disclosed systems, devices, and methods may be implemented in other ways. Furthermore, the couplings or direct couplings or communication connections shown or discussed may be indirect couplings or communication connections through interfaces, devices, or units, and may be electrical, mechanical, or other forms.

[0186] The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed across multiple network units. Some or all of these units may be selected to achieve the purpose of this embodiment according to actual needs.

[0187] In addition, each functional unit in each embodiment of the present application 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.

[0188] It should be understood that in the embodiments of the present application, the various embodiments and features can be combined with each other to solve the aforementioned technical problems.

[0189] If the functions are implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present application, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes several instructions for enabling a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the method described in each embodiment of the present application. The aforementioned storage medium includes various media that can store program codes, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk.

[0190] The above description of the disclosed embodiments is intended to enable one skilled in the art to implement or use the present application. Various modifications to these embodiments will be readily apparent to one skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the present application. Therefore, the present application is not limited to the embodiments shown herein, but is intended to conform to the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A method for calculating the storage capacity of an aquifer gas storage reservoir, which is used for a gas storage reservoir with an open aquifer water body, characterized in that: Including steps: S11. Acquire geological data and test data of the groundwater layer to determine the ultimate pore volume and ultimate storage capacity of the aquifer-type gas storage reservoir; the geological data include trap characteristic data and structural characteristic data for determining the ultimate pore volume of the aquifer-type gas storage reservoir; the test data include pressure data, injection and production data, and core flooding test data; S12. Acquire water property data of the aquifer, and construct a first functional relationship group of the gas-water interface of the gas storage reservoir and the change of the aquifer pressure based on the water property data; the water property data includes water pressure and water intrusion coefficient; S13, obtaining multiple rounds of gas-water mutual displacement experimental data, and obtaining a second functional relationship group of the two-phase region of the phase permeability curve changing with the displacement rounds; the gas-water mutual displacement experimental data includes displacement pressure, displacement time, cumulative gas production, and cumulative water production; S14, obtaining the gas storage pore volume corresponding to the irreducible water saturation and the movable gas pore volume corresponding to the movable gas saturation in different displacement cycles; S15. Calculate the storage capacity of the aquifer-type gas storage reservoir under the displacement cycle based on the first functional relationship group, the second functional relationship group, and the gas storage pore volume and movable gas pore volume corresponding to the displacement cycle; the storage capacity includes the dynamic storage capacity and the movable gas storage capacity.

2. The storage capacity calculation method for aquifer gas storage according to claim 1, characterized in that: The calculation formula of the limiting pore volume includes: Among them, V max is the limiting pore volume, is the average porosity, A max is the limiting enclosed area, h max is the trap closure height.

3. The storage capacity calculation method for aquifer gas storage according to claim 2, characterized in that: The calculation formula of the limit storage capacity is include: G max =V max ·(1-S wc ); Among them, G max is the limit storage capacity, S wc The irreducible water saturation is obtained from the gas-water displacement experiment in the core flooding experiment.

4. The storage capacity calculation method for aquifer gas storage according to claim 3, characterized in that: The first functional relationship group of the gas-water interface of the gas storage reservoir and the change of the aquifer pressure includes: The formula for obtaining the pore volume of the aquifer corresponding to the gas-water interface at time t is: and, According to the formula The relationship between water intrusion and aquifer pressure W e =f(P); Among them, V t is the pore volume of the aquifer, is the average porosity, A is the gas-bearing area, is the air-water interface at time t, C s is the water intrusion coefficient, W e is the water intrusion, and f(P) is the functional relationship between the water intrusion and the aquifer pressure.

5. The storage capacity calculation method for aquifer gas storage according to claim 4, characterized in that: The method of obtaining multiple rounds of gas-water mutual displacement experimental data and obtaining a second functional relationship between the two-phase region of the phase permeability curve and the displacement rounds includes: The gas-water mutual displacement experimental data is used to obtain the gas-water relative permeability curve; wherein, the relative permeability curve processing process of the gas-water displacement process includes: S21. Correct the cumulative liquid production value measured at the rock sample outlet pressure to the value at the average rock sample pressure using the following formula: Among them, V i is the cumulative liquid production value at time i, ΔV wi is the water production from time i-1 to time i, V i-1 is the cumulative liquid production value at time i-1, P a is atmospheric pressure, ΔP is displacement pressure difference, ΔV gi It is the gas increment value measured at a certain time interval under atmospheric pressure; S22, correct the cumulative liquid production value, and use the non-steady-state gas-water relative permeability method to calculate the end-face gas saturation S at different times of gas-water displacement. ge , water phase relative permeability K rw and the gas phase relative permeability K rg , the formula is as follows: Among them, f w (S g ) is the moisture content, Cumulative water production, is the cumulative liquid production, K rw is the relative permeability of water phase, K rg is the gas phase relative permeability, μ g is the viscosity of natural gas under formation conditions, μ w is the formation water viscosity under formation conditions, I is the flow capacity ratio, Q t is the liquid volume at the core outlet end at time t, Q0 is the water production at the core outlet end at the initial time, S ge The gas saturation of the end face at different times of gas-driven water displacement; The second functional relationship group includes: According to the formula S wc =f1(n) is the functional relationship between irreducible water saturation and displacement rounds; According to the formula S gr =f2(n) The function relationship between the residual gas saturation and the displacement round number; According to the formula S glive =f3(n) The functional relationship between the movable gas saturation and the displacement round number; Among them, S wc is the irreducible water saturation, n is the number of displacement cycles, f1(n) is the functional relationship between irreducible water saturation and displacement cycles, S gr is the residual gas saturation, f2(n) is the functional relationship between the residual gas saturation and the displacement rounds, S glive is the movable gas saturation, and f3(n) is the functional relationship between the movable gas saturation and the displacement rounds.

6. The storage capacity calculation method for aquifer gas storage according to claim 5, characterized in that: The formula for obtaining the gas storage pore volume corresponding to the irreducible water saturation is as follows: V wc =V t (1-S wc ) Among them, V wc is the gas storage pore volume of aquifer-type gas storage; The movable air pore volume corresponding to the movable water saturation is obtained using the following formula: V glive =V t S glive Among them, V glive is the movable gas pore volume of the aquifer-type gas storage.

7. The storage capacity calculation method for aquifer gas storage according to claim 6, characterized in that: Calculating the storage capacity of the aquifer-type gas storage reservoir under the displacement rounds, including: S31. Calculate the dynamic storage capacity G of aquifer-type gas storage g ,include: Among them, V wc is the gas storage pore volume of aquifer-type gas storage; B g is the volume coefficient of natural gas at formation pressure; S32. Calculate the maximum value G of the dynamic library capacity gmax ,include: The dynamic storage capacity G g Constrained by the formation upper limit pressure, when the injection and production capacity is expanded to the overflow point, if the formation pressure does not break through the formation upper limit pressure, the dynamic storage capacity G g The maximum value of G gmax The calculation formula includes: G gmax =G max If the formation pressure exceeds the formation upper limit pressure, the dynamic storage capacity G g The maximum value of G gmax The calculation formula includes: Among them, V t,limit is the pore volume of the aquifer corresponding to the formation upper limit pressure during the injection and production expansion process, and the calculation formula is: in, The gas-water interface corresponding to the formation upper limit pressure during the injection-production expansion process; S33. Calculate the movable gas storage capacity G of aquifer-type gas storage glive The formulas include: Among them, B g V is the volume coefficient of natural gas under formation pressure; glive is the movable gas pore volume of the aquifer-type gas storage.

8. A storage capacity calculation device for an aquifer gas storage reservoir, which is used for a gas storage reservoir with an open aquifer water body, characterized in that: include: a data acquisition unit, configured to acquire geological data and test data of the groundwater layer to determine the ultimate pore volume and ultimate storage capacity of the aquifer-type gas storage reservoir; the geological data includes trap characteristic data and structural characteristic data, which are used to determine the ultimate pore volume of the aquifer-type gas storage reservoir; the test data includes pressure data, injection and production data, and core flooding experiment data; a first function acquisition unit configured to acquire water property data of the aquifer and construct a first functional relationship group of changes in the gas-water interface of the gas storage reservoir with the aquifer pressure based on the water property data; the water property data including the water pressure and the water intrusion coefficient; A second function acquisition unit is used to acquire multiple rounds of gas-water mutual displacement experimental data and obtain a second function relationship group of the two-phase region of the phase permeability curve changing with the displacement rounds; the gas-water mutual displacement experimental data includes displacement pressure, displacement time, cumulative gas production and cumulative water production; A pore volume calculation unit is used to obtain the gas storage pore volume corresponding to the irreducible water saturation in different displacement cycles, and the movable gas pore volume corresponding to the movable gas saturation; A storage capacity calculation unit is used to calculate the storage capacity of the aquifer-type gas storage reservoir under the displacement cycle based on the first functional relationship group, the second functional relationship group, and the gas storage pore volume and movable gas pore volume corresponding to the displacement cycle; the storage capacity includes dynamic storage capacity and movable gas storage capacity.

9. A storage capacity calculation device for an aquifer gas storage reservoir, characterized in that: include: memory for storing computer programs; A processor is used to call and execute the computer program to implement the steps of the storage capacity calculation method for an aquifer gas storage reservoir as described in any one of claims 1 to 7.

10. A storage medium, characterized in that: The method comprises a software program adapted to execute, by a processor, the steps of the storage capacity calculation method for an aquifer gas storage reservoir as claimed in any one of claims 1 to 7.

Images