Method for testing aquifer reservoir geological body, electronic device, and program product

Abstract

This application discloses a testing method, electronic equipment, and program product for aquifer gas storage geological bodies, relating to the field of underground gas storage design and construction technology. It aims to provide a scheme for testing geological bodies used in the construction of gas storage facilities utilizing aquifers. The method includes: acquiring location data of injection wells and monitoring wells at the aquifer; the first temperature of the injection fluid injected into the aquifer at the injection well; the second temperature of the aquifer; the effective porosity of the aquifer; and the viscosity of the injection fluid at the first temperature; acquiring background monitoring data; acquiring monitoring data during injection testing at the injection well; and retrieving the permeability, vertical dispersion, and lateral dispersion of the aquifer in three dimensions based on the tracer breakthrough curve, the first temperature, the second temperature, the effective porosity, the hydrostatic pressure difference between the injection well and the monitoring well, and based on a water-heat-mass coupling solute transport analytical model.

Description

Technical Field

[0001] This application relates to the field of underground gas storage design and construction technology, and in particular to a geological testing method, electronic equipment and program product for aquifer gas storage. Background Technology

[0002] Underground gas storage facilities play an important role in ensuring energy security, optimizing pipeline operation, supporting renewable energy, improving economic efficiency, and protecting the environment.

[0003] Current gas storage facilities are mainly oil and gas reservoir type and salt cavern type gas storage facilities. The structural resources can no longer meet the requirements for gas storage facility construction. In order to further expand the site selection range of gas storage facilities, it is necessary to further explore the structure of aquifers and other structures.

[0004] Therefore, it became our need to test the geological bodies of aquifers. Summary of the Invention

[0005] The purpose of this application is to provide a method, electronic equipment and program product for testing geological bodies of aquifer gas storage facilities, aiming to provide a solution for testing geological bodies used to construct gas storage facilities using aquifers.

[0006] In a first aspect, this application provides a method for testing geological bodies in an aquifer gas reservoir. The method includes: acquiring location data of injection wells and monitoring wells at the aquifer; a first temperature of the injection fluid injected into the aquifer at the injection well; a second temperature of the aquifer; the effective porosity of the aquifer; and the viscosity of the injection fluid at the first temperature; the injection fluid includes a tracer; acquiring background monitoring data; the background monitoring data consists of baseline pressure and temperature data of the aquifer collected by the monitoring well, assuming both the injection well and the monitoring well are shut in and the pressure and temperature of the aquifer remain stable for a preset time period; acquiring monitoring data during the injection test at the injection well; the monitoring data includes pump pressure at the injection well, injection operation start and stop times, and aquifer pressure curve, caprock overburden permeability pressure curve, and tracer breakthrough curve at the monitoring well; and determining the appropriate parameters by comparing the aquifer pressure curve with the injection operation start and stop times and the pressure baseline data. Aquifer connectivity is determined as follows: If the aquifer pressure curve shows fluctuations deviating from the pressure baseline at the injection operation start and stop times, the aquifer is considered connected; if the aquifer pressure curve remains consistent with the pressure baseline at the injection operation start and stop times, the aquifer is considered disconnected. Caprock sealing is determined by comparing the caprock overburden permeability pressure curve with the injection operation start and stop times and the pressure baseline: If the caprock overburden permeability pressure curve remains consistent with the pressure baseline at the injection operation start and stop times, the caprock is considered sealed; if the caprock overburden permeability pressure curve shows fluctuations deviating from the pressure baseline at the injection operation start and stop times, the caprock is considered leaking. Based on the tracer breakthrough curve, first temperature, second temperature, effective porosity, hydrostatic pressure difference between the injection well and the monitoring well, and based on the water-heat-mass coupling solute transport analytical model, the aquifer's permeability, vertical dispersion, and lateral dispersion in three dimensions are inverted.

[0007] Optionally, based on the tracer breakthrough curve, first temperature, second temperature, effective porosity, hydrostatic pressure difference between the injection well and the monitoring well, and based on the water-thermal-mass coupled solute transport analytical model, the permeability, vertical dispersion, and lateral dispersion of the aquifer in three dimensions are inverted. This includes: determining the three-dimensional velocity components of the injected fluid in the aquifer based on the peak arrival time and location data of the tracer, and calculating the synthetic average velocity; calculating the vertical and lateral spatial variances of the injected fluid based on the tracer breakthrough curve; and using the synthetic average velocity, vertical spatial variance, lateral spatial variance, and peak arrival time, inverting the initial values ​​of the vertical dispersion and lateral dispersion through moment analysis, and explicitly subtracting the values ​​obtained from the first... Temperature-corrected molecular diffusion terms determined by a first temperature and a second temperature; initial values ​​of three-dimensional permeability are inverted using three-dimensional velocity components, hydrostatic pressure difference, and fluid viscosity corrected for viscosity and temperature coefficient of viscosity; based on an analytical model of water-heat-mass coupling solute transport including temperature-corrected molecular diffusion terms, the three-dimensional velocity components are used as input to the analytical model, and multiple sets of parameter samples are generated based on initial values ​​of longitudinal dispersion, lateral dispersion, and three-dimensional permeability, and corresponding breakthrough simulation curves are calculated. The peak concentration, peak arrival time, and at least one curve broadening characteristic time of the breakthrough simulation curve and the tracer breakthrough curve are used to construct the fitting error. The Monte Carlo method is used to search for the parameter combination that minimizes the fitting error, resulting in optimized longitudinal dispersion, lateral dispersion, and three-dimensional permeability.

[0008] Optionally, the analytical model for solute transport based on water-heat-mass coupling satisfies the following relationship: ; in, Indicates solute concentration; Indicates time; , , Represents three-dimensional spatial coordinates; Indicates the delay factor; Indicates the injected fluid in Diffusion coefficient in the direction of travel; Indicates the injected fluid in Diffusion coefficient in the direction of travel; Indicates the injected fluid in Diffusion coefficient in the direction of travel; Indicates the injected fluid in Flow velocity component in the direction; Indicates the injected fluid in Flow velocity component in the direction; Indicates the injected fluid in Flow velocity component in the direction; Indicates source and sink terms; delay factor Determined by the following formula: ; in, This represents the solid-phase medium volume density of the aquifer; Indicates the linear adsorption partition coefficient of the solute; Indicates effective porosity.

[0009] Optionally, the diffusion coefficient satisfies the following relationship: ; in, Indicates lateral dispersion. Represents the temperature-corrected molecular diffusion term; ; in, This represents the diffusion coefficient of the injected fluid at the first temperature; This indicates the viscosity of the injected fluid at the first temperature; This indicates the viscosity of the injected fluid at the second temperature; ; in, This represents the viscosity temperature coefficient.

[0010] Optionally, the solute concentration satisfies the following relationship: ; in, Represents the unit impulse response function; Indicates time The varying injection concentration rate; This indicates the upper limit of the simulation duration.

[0011] Optionally, the longitudinal dispersion and the lateral dispersion are determined according to the following formulas: ; ; in, Indicates longitudinal dispersion; Indicates the composite average flow rate; Indicates lateral dispersion; Represents the longitudinal spatial variance; Indicates the lateral spatial variance; This indicates the peak arrival time.

[0012] Alternatively, the three-dimensional permeability can be inverted according to the following formula: ; in, Indicates that the aquifer is in Permeability in the direction of travel; Indicates that the aquifer is in Permeability in the direction of travel; Indicates that the aquifer is in Permeability in the direction of travel; Indicates the hydrostatic pressure difference of the fluid; Indicates the relationship between the injection well and the monitoring well. Coordinate difference in direction; Indicates the relationship between the injection well and the monitoring well. Coordinate difference in direction; Indicates the relationship between the injection well and the monitoring well. The difference in coordinates along the direction.

[0013] Optionally, multiple sets of parameter samples are generated according to the following formula: ; in, No. Parameters in the group parameter sample ; Indicates parameters The prior minimum value; Indicates parameters The prior maximum value; This represents a random number that follows a triangular distribution; The fitting error satisfies the following relationship: ; in, Indicates the first Group parameter sample parameters The corresponding fitting error; Represents the normalization coefficient; This indicates the theoretical peak concentration that breaks through the simulated curve; Indicates the peak concentration of the tracer breakthrough curve; This indicates the theoretical peak arrival time determined based on the breakthrough simulation curve; Indicates the peak arrival time determined based on the tracer breakthrough curve; This indicates the broadening characteristic time of the theoretical curve extracted from the breakthrough simulation curve; This represents the curve broadening characteristic time extracted from the tracer breakthrough curve.

[0014] This application provides a geological testing method for aquifer gas storage facilities, specifically addressing the current situation where existing gas storage facilities are mainly oil and gas reservoirs and salt cavern types, and structural resources are increasingly scarce. This method offers a geological testing approach for gas storage facilities built using aquifers. By acquiring background pressure / temperature baselines for injection and monitoring wells, reservoir pressure curves during injection testing, overburden permeability pressure curves, and tracer breakthrough curves, and combining this with a water-thermal-mass coupled solute transport analytical model, this method achieves a single test and multi-parameter inversion of aquifer connectivity, caprock sealing, and reservoir three-dimensional permeability, vertical dispersion, and lateral dispersion. Compared to conventional isothermal tracer testing, this application fully considers the influence of the temperature difference between the injected fluid and the gas reservoir on fluid viscosity and molecular diffusion coefficient during the inversion process, effectively overcoming the parameter inversion distortion problem caused by traditional isothermal models in scenarios with significant temperature differences. This provides a more reliable basis for key parameters in the layout of injection and production wells, capacity assessment, and sealing evaluation of aquifer-converted gas storage facilities.

[0015] Secondly, this application provides a geological testing device for aquifer gas storage reservoirs, the device comprising various functional modules for the method described in the first aspect above.

[0016] Thirdly, this application provides an electronic device comprising: a processor and a memory; the memory storing instructions executable by the processor; the processor being configured to, when executing the instructions, cause the electronic device to perform the method provided in the first aspect above.

[0017] Fourthly, this application provides a computer-readable storage medium comprising: computer software instructions; when the computer software instructions are executed in an electronic device, the electronic device causes the electronic device to implement the method provided in the first aspect above.

[0018] Fifthly, this application provides a computer program product including computer instructions that, when executed on an electronic device, cause the electronic device to perform the method provided in the first aspect above.

[0019] The beneficial effects of aspects two through five above can be referred to in aspect one, and will not be repeated here. Attached Figure Description

[0020] To more clearly illustrate the technical solutions of the embodiments of this application, the drawings used in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0021] Figure 1A schematic diagram illustrating the application scenario of the geological body testing method for aquifer gas storage provided in this application embodiment; Figure 2 A schematic flowchart illustrating the geological testing method for aquifer gas storage provided in this application embodiment; Figure 3 A schematic diagram of the composition of the geological testing device for aquifer gas storage provided in this application embodiment; Figure 4 This is a schematic diagram illustrating the composition of an electronic device provided in an embodiment of this application. Detailed Implementation

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

[0023] The terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of technical features indicated. Therefore, a feature defined as "first" or "second" may explicitly or implicitly include one or more of that feature. In the description of this application, unless otherwise stated, "a plurality of" means two or more.

[0024] In embodiments of this application, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, article, or apparatus that includes that element.

[0025] In the embodiments of this application, the words "exemplarily" or "for example" are used to indicate examples, illustrations, or explanations. Any embodiment or design described as "exemplarily" or "for example" in the embodiments of this application should not be construed as being more preferred or advantageous than other embodiments or design solutions. Specifically, the use of the words "exemplarily" or "for example" is intended to present the relevant concepts in a specific manner.

[0026] Underground gas storage facilities play an important role in ensuring energy security, optimizing pipeline operation, supporting renewable energy, improving economic efficiency, and protecting the environment.

[0027] Current gas storage facilities are mainly oil and gas reservoir type and salt cavern type gas storage facilities. The structural resources can no longer meet the requirements for gas storage facility construction. In order to further expand the site selection range of gas storage facilities, it is necessary to further explore the structure of aquifers and other structures.

[0028] Therefore, it became our need to test the geological bodies of aquifers.

[0029] Based on this, embodiments of this application provide a method, electronic equipment, and program product for testing geological bodies in aquifer gas storage facilities, aiming to provide a solution for testing geological bodies in aquifers.

[0030] The following description is provided in conjunction with the accompanying drawings.

[0031] Figure 1 This is a schematic diagram illustrating an application scenario for the geological testing method for aquifer gas storage provided in this application embodiment. For example... Figure 1 As shown, the vertical wellbore on the left is the injection well, with fiber optic cables running throughout. Electronic and fiber optic pressure gauges are installed at the bottom of the wellbore, and bidirectional water / gas arrows at the bottom of the formation indicate the injection of a water-gas mixture containing tracers into the formation. The wellbore on the right is the monitoring well, also with fiber optic cables running throughout. Multiple packers are arranged vertically at intervals within the well to achieve formation layer isolation. Electronic pressure gauges are installed between the packers, and a direct-reading fiber optic pressure gauge is installed at the bottom of the wellbore. The formation is divided from top to bottom into overlying permeable layer 3, overlying permeable layer 2, overlying permeable layer 1, and aquifer. The entire system can simultaneously acquire pump injection pressure, layered caprock pressure time-series curves, global temperature, and fluid acoustic vibration signals. Combined with inter-well tracer migration data, it can qualitatively determine reservoir connectivity and caprock sealing, and quantitatively invert formation seepage parameters.

[0032] It should be noted that the execution entity of the aquifer gas reservoir geological body testing method provided in this application embodiment is an aquifer gas reservoir geological body testing device (hereinafter referred to as the testing device). This testing device can be an electronic device with computing processing capabilities, such as a computer or server. Optionally, the testing device can also be a processor (e.g., a central processing unit, CPU) in the aforementioned electronic device; or, the testing device can also be an application (APP) with aquifer gas reservoir geological body testing function installed in the aforementioned electronic device; or, the testing device can also be a software system or platform deployed in the aforementioned electronic device; or, the testing device can also be a functional module in the aforementioned electronic device used to execute the aquifer gas reservoir geological body testing method, etc. This application embodiment does not impose any limitations on these aspects.

[0033] For simplicity, the following description will use the testing device as the execution subject of the aquifer gas reservoir geological body testing method provided in the embodiments of this application as an example.

[0034] Figure 2 This is a schematic flowchart illustrating the geological testing method for aquifer gas storage reservoirs provided in this application embodiment. Figure 2 As shown, the method may include the following steps: S101. Obtain the location data of the injection well and monitoring well at the aquifer, the first temperature of the injection fluid injected into the aquifer through the injection well, the second temperature of the aquifer, the effective porosity of the aquifer, and the viscosity of the injection fluid at the first temperature.

[0035] The injected fluid includes a tracer.

[0036] S102. Obtain background monitoring data.

[0037] Among them, the background monitoring data are the pressure baseline data and temperature baseline data of the aquifer collected by the monitoring well when both the injection well and the monitoring well are shut in and the pressure and temperature of the aquifer remain stable within a preset time.

[0038] As an example, maintaining temperature can mean that the fluctuation value of pressure or temperature is less than or equal to a preset threshold.

[0039] S103. Acquire monitoring data during the injection test in the injection well; the monitoring data includes the pump pressure at the injection well, the start and stop times of the injection operation, and the aquifer pressure curve, caprock permeability pressure curve, and tracer breakthrough curve at the monitoring well.

[0040] S104. By comparing the aquifer pressure curve with the injection operation start and stop time points and pressure baseline data, the connectivity of the aquifer is determined: if the aquifer pressure curve shows fluctuations deviating from the pressure baseline data at the injection operation start and stop time points, the aquifer is determined to be connected; if the aquifer pressure curve remains consistent with the pressure baseline data at the injection operation start and stop time points, the aquifer is determined to be disconnected.

[0041] S105. By comparing the overlying permeation pressure curve with the injection operation start and stop time points and pressure baseline data, the sealing performance of the caprock is determined: if the overlying permeation pressure curve remains consistent with the pressure baseline data at the injection operation start and stop time points, the caprock is confirmed to be sealed; if the overlying permeation pressure curve shows fluctuations deviating from the pressure baseline data at the injection operation start and stop time points, the caprock is confirmed to be leaking.

[0042] S106. Based on the tracer breakthrough curve, first temperature, second temperature, effective porosity, hydrostatic pressure difference between the injection well and the monitoring well, and based on the water-heat-mass coupling solute transport analytical model, the permeability, vertical dispersion and horizontal dispersion of the aquifer in three dimensions are inverted.

[0043] In some embodiments, S106 may specifically include the following steps: Step 1: Based on the peak arrival time and location data of the tracer, determine the three-dimensional velocity components of the injected fluid in the aquifer and calculate the synthetic average velocity.

[0044] Step 2: Calculate the longitudinal and transverse spatial variances of the injected fluid based on the tracer breakthrough curve.

[0045] Step 3: Using the synthetic average flow rate, longitudinal spatial variance, lateral spatial variance, and peak arrival time, the initial values ​​of longitudinal and lateral diffusion are inverted through moment analysis, and the temperature-corrected molecular diffusion term determined by the first and second temperatures is explicitly subtracted in the moment analysis.

[0046] Step 4: Using the three-dimensional velocity components, the hydrostatic pressure difference of the fluid, and the fluid viscosity corrected according to the viscosity and viscosity temperature coefficient, the initial value of the three-dimensional permeability is inverted.

[0047] Step 5: Based on the water-thermal-mass coupling solute transport analytical model that includes a temperature-corrected molecular diffusion term, the three-dimensional flow velocity component is used as the input to the water-thermal-mass coupling solute transport analytical model. Based on the initial values ​​of longitudinal dispersion, lateral dispersion, and three-dimensional permeability, multiple sets of parameter samples are generated and the corresponding breakthrough simulation curves are calculated. The peak concentration, peak arrival time, and at least one curve broadening characteristic time of the breakthrough simulation curve and the tracer breakthrough curve are used to construct the fitting error. The Monte Carlo method is used to search for the parameter combination that minimizes the fitting error, and the optimized longitudinal dispersion, lateral dispersion, and three-dimensional permeability are obtained.

[0048] As an example, the analytical model of solute transport based on water-heat-mass coupling satisfies the following relationship: Formula (1) in, Indicates solute concentration; Indicates time; , , Represents three-dimensional spatial coordinates; Indicates the delay factor; Indicates the injected fluid in Diffusion coefficient in the direction of travel; Indicates the injected fluid in Diffusion coefficient in the direction of travel; Indicates the injected fluid in Diffusion coefficient in the direction of travel; Indicates the injected fluid in Flow velocity component in the direction; Indicates the injected fluid in Flow velocity component in the direction; Indicates the injected fluid in Flow velocity component in the direction; Indicates source and sink terms; delay factor Determined by the following formula: Formula (2) in, This represents the solid-phase medium volume density of the aquifer; Indicates the linear adsorption partition coefficient of the solute; Indicates effective porosity.

[0049] As an example, the velocity components satisfy the following relationship: Formula (3) in, This indicates that the injection well and the monitoring well are in Coordinate difference in direction; This indicates that the injection well and the monitoring well are in Coordinate difference in direction; This indicates that the injection well and the monitoring well are in Coordinate difference in direction; This indicates the peak arrival time.

[0050] As an example, the diffusion coefficient satisfies the following relationship: Formula (4) in, Indicates lateral dispersion. This represents the temperature-corrected molecular diffusion term.

[0051] Formula (5) in, This represents the diffusion coefficient of the injected fluid at the first temperature; This indicates the viscosity of the injected fluid at the first temperature; This indicates the viscosity of the injected fluid at the second temperature.

[0052] Formula (6) in, This represents the viscosity temperature coefficient.

[0053] As an example, the solute concentration satisfies the following relationship: Formula (7) in, Represents the unit impulse response function; Indicates time The varying injection concentration rate; This indicates the upper limit of the simulation duration.

[0054] As an example, longitudinal diffusion and lateral diffusion are determined according to the following formulas: Formula (8) Formula (9) in, Indicates longitudinal dispersion; Indicates the composite average flow rate; Indicates lateral dispersion; Represents the longitudinal spatial variance; Indicates the lateral spatial variance; This indicates the peak arrival time.

[0055] As an example, three-dimensional permeability is inverted according to the following formula: Formula (10) in, Indicates that the aquifer is in Permeability in the direction of travel; Indicates that the aquifer is in Permeability in the direction of travel; Indicates that the aquifer is in Permeability in the direction of travel; Indicates the hydrostatic pressure difference of the fluid; Indicates the relationship between the injection well and the monitoring well. Coordinate difference in direction; Indicates the relationship between the injection well and the monitoring well. Coordinate difference in direction; Indicates the relationship between the injection well and the monitoring well. The difference in coordinates along the direction.

[0056] As an example, multiple sets of parameter samples are generated according to the following formula: Formula (11) in, No. Parameters in the group parameter sample ; Indicates parameters The prior minimum value; Indicates parameters The prior maximum value; This represents a random number that follows a triangular distribution.

[0057] The fitting error satisfies the following relationship: Formula (12) in, Indicates the first Group parameter sample parameters The corresponding fitting error; Represents the normalization coefficient; This indicates the theoretical peak concentration that breaks through the simulated curve; Indicates the peak concentration of the tracer breakthrough curve; This indicates the theoretical peak arrival time determined based on the breakthrough simulation curve; Indicates the peak arrival time determined based on the tracer breakthrough curve; This indicates the broadening characteristic time of the theoretical curve extracted from the breakthrough simulation curve; This represents the curve broadening characteristic time extracted from the tracer breakthrough curve.

[0058] For example, the parameter sample may include three-dimensional permeability, longitudinal diffusion, and lateral diffusion.

[0059] The geological testing method for aquifer gas storage provided in this application addresses the current situation where existing gas storage construction is mainly of oil and gas reservoir and salt cavern type, and structural resources are increasingly scarce. This method specifically provides a geological testing approach for constructing blow-up storage facilities using aquifers. By acquiring the background pressure / temperature baselines of the injection well and monitoring well, the reservoir pressure curve during injection testing, the overburden permeability pressure curve, and the tracer breakthrough curve, and combining this with a water-thermal-mass coupled solute transport analytical model, this method achieves a single test and multi-parameter inversion of aquifer connectivity, caprock sealing, and reservoir three-dimensional permeability, vertical dispersion, and lateral dispersion. Compared with conventional isothermal tracer testing, this application fully considers the influence of the temperature difference between the injected fluid and the gas reservoir on fluid viscosity and molecular diffusion coefficient during the inversion process. This effectively overcomes the parameter inversion distortion problem caused by traditional isothermal models in scenarios with significant temperature differences, providing a more reliable key parameter basis for injection and production well layout, storage capacity assessment, and sealing evaluation of aquifer-converted gas storage facilities.

[0060] The foregoing primarily describes the solutions provided by the embodiments of this application from a methodological perspective. To achieve the aforementioned functions, each device, such as a testing apparatus, includes corresponding hardware structures and / or software modules for executing each function. Those skilled in the art should readily recognize that, in conjunction with the algorithmic steps of the examples described in the embodiments disclosed herein, this application can be implemented in hardware or a combination of hardware and computer software. Whether a function is executed in hardware or by computer software driving hardware depends on the specific application and design constraints of the technical solution. Experts may 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.

[0061] This application embodiment can divide the testing device into functional modules according to the above method embodiment. For example, each function can be divided into a separate functional module, or two or more functions can be integrated into one functional module. The integrated module can be implemented in hardware or software. It should be noted that the module division in this application embodiment is illustrative and only represents one logical functional division. In actual implementation, there may be other division methods. The following description uses the example of dividing each functional module according to each function.

[0062] In an exemplary embodiment, this application provides a geological testing device for aquifer gas storage reservoir. Figure 3 A schematic diagram illustrating the composition of the geological testing device for an aquifer gas storage reservoir provided in this embodiment of the application. Figure 3 As shown, the geological testing device for the aquifer gas reservoir includes: an acquisition module 301 and a processing module 302.

[0063] The acquisition module 301 is used to acquire the location data of the injection well and the monitoring well at the aquifer, the first temperature of the injection fluid injected into the aquifer at the injection well, the second temperature of the aquifer, the effective porosity of the aquifer, and the viscosity of the injection fluid at the first temperature; the injection fluid includes a tracer; acquire background monitoring data; the background monitoring data consists of the pressure baseline data and temperature baseline data of the aquifer collected by the monitoring well when both the injection well and the monitoring well are in a shut-in state and the pressure and temperature of the aquifer remain stable within a preset time period; acquire monitoring data during the injection test at the injection well; the monitoring data includes the pump pressure at the injection well, the start and stop times of the injection operation, and the aquifer pressure curve, the overburden permeability pressure curve, and the tracer breakthrough curve at the monitoring well; Processing module 302 is used to determine the connectivity of the aquifer by comparing the aquifer pressure curve with the injection operation start and stop time points and pressure baseline data: if the aquifer pressure curve fluctuates and deviates from the pressure baseline data at the injection operation start and stop time points, the aquifer is determined to be connected; if the aquifer pressure curve remains consistent with the pressure baseline data at the injection operation start and stop time points, the aquifer is determined to be disconnected. It also determines the caprock sealing by comparing the overburden permeability pressure curve with the injection operation start and stop time points and pressure baseline data: if the overburden permeability pressure curve remains consistent with the pressure baseline data at the injection operation start and stop time points, the caprock is determined to be sealed; if the overburden permeability pressure curve fluctuates and deviates from the pressure baseline data at the injection operation start and stop time points, the caprock is determined to be leaking. Based on the tracer breakthrough curve, first temperature, second temperature, effective porosity, hydrostatic pressure difference between the injection well and the monitoring well, and based on the water-heat-mass coupling solute transport analytical model, it inversely calculates the aquifer's permeability, vertical dispersion, and lateral dispersion in three dimensions.

[0064] In some embodiments, the processing module 302 is specifically used to determine the three-dimensional velocity components of the injected fluid in the aquifer based on the peak arrival time and location data of the tracer, and to calculate the synthetic average velocity; to calculate the longitudinal spatial variance and lateral spatial variance of the injected fluid based on the tracer breakthrough curve; to use the synthetic average velocity, longitudinal spatial variance, lateral spatial variance, and peak arrival time to invert the initial values ​​of longitudinal dispersion and lateral dispersion through moment analysis, and to explicitly subtract the temperature-corrected molecular diffusion term determined by the first temperature and the second temperature in the moment analysis; and to use the three-dimensional velocity components, the hydrostatic pressure difference, and the root... The initial value of three-dimensional permeability is inverted based on the viscosity and viscosity-temperature coefficient-corrected fluid viscosity. Using a water-thermal-mass coupling solute transport analytical model incorporating a temperature-corrected molecular diffusion term, the three-dimensional velocity component is taken as input. Based on the initial values ​​of longitudinal dispersion, lateral dispersion, and three-dimensional permeability, multiple sets of parameter samples are generated, and corresponding breakthrough simulation curves are calculated. The peak concentration, peak arrival time, and at least one curve broadening characteristic time of the breakthrough simulation curve and the tracer breakthrough curve are used to construct the fitting error. The Monte Carlo method is then used to search for the parameter combination that minimizes the fitting error, resulting in optimized longitudinal dispersion, lateral dispersion, and three-dimensional permeability.

[0065] In some embodiments, the analytical model for solute transport based on water-heat-mass coupling satisfies the following relationship: ; in, Indicates solute concentration; Indicates time; , , Represents three-dimensional spatial coordinates; Indicates the delay factor; Indicates the injected fluid in Diffusion coefficient in the direction of travel; Indicates the injected fluid in Diffusion coefficient in the direction of travel; Indicates the injected fluid in Diffusion coefficient in the direction of travel; Indicates the injected fluid in Flow velocity component in the direction; Indicates the injected fluid in Flow velocity component in the direction; Indicates the injected fluid in Flow velocity component in the direction; Indicates source and sink terms; delay factor Determined by the following formula: ; in, This represents the solid-phase medium volume density of the aquifer; Indicates the linear adsorption partition coefficient of the solute; Indicates effective porosity.

[0066] In some embodiments, the diffusion coefficient satisfies the following relationship: ; in, Indicates lateral dispersion. Represents the temperature-corrected molecular diffusion term; ; in, This represents the diffusion coefficient of the injected fluid at the first temperature; This indicates the viscosity of the injected fluid at the first temperature; This indicates the viscosity of the injected fluid at the second temperature; ; in, This represents the viscosity temperature coefficient.

[0067] In some embodiments, the solute concentration satisfies the following relationship: ; in, Represents the unit impulse response function; Indicates time The varying injection concentration rate; This indicates the upper limit of the simulation duration.

[0068] In some embodiments, longitudinal dispersion and lateral dispersion are determined according to the following formula: ; ; in, Indicates longitudinal dispersion; Indicates the composite average flow rate; Indicates lateral dispersion; Represents the longitudinal spatial variance; Indicates the lateral spatial variance; This indicates the peak arrival time.

[0069] In some embodiments, the three-dimensional permeability is inverted according to the following formula: ; in, Indicates that the aquifer is in Permeability in the direction of travel; Indicates that the aquifer is in Permeability in the direction of travel; Indicates that the aquifer is in Permeability in the direction of travel; Indicates the hydrostatic pressure difference of the fluid; Indicates the relationship between the injection well and the monitoring well. Coordinate difference in direction; Indicates the relationship between the injection well and the monitoring well. Coordinate difference in direction; Indicates the relationship between the injection well and the monitoring well. The difference in coordinates along the direction.

[0070] In some embodiments, multiple sets of parameter samples are generated according to the following formula: ; in, No. Parameters in the group parameter sample ; Indicates parameters The prior minimum value; Indicates parameters The prior maximum value; This represents a random number that follows a triangular distribution; The fitting error satisfies the following relationship: ; in, Indicates the first Group parameter sample parameters The corresponding fitting error; Represents the normalization coefficient; This indicates the theoretical peak concentration that breaks through the simulated curve; Indicates the peak concentration of the tracer breakthrough curve; This indicates the theoretical peak arrival time determined based on the breakthrough simulation curve; Indicates the peak arrival time determined based on the tracer breakthrough curve; This indicates the broadening characteristic time of the theoretical curve extracted from the breakthrough simulation curve; This represents the curve broadening characteristic time extracted from the tracer breakthrough curve.

[0071] It should be noted that the above Figure 3 Modules in a module can also be called units; for example, a processing module can be called a processing unit. Additionally, in... Figure 3 In the embodiments shown, the names of the modules may not be the same as those shown in the figure. For example, the acquisition module may also be called the transceiver module or the communication module.

[0072] Figure 3If the various modules in the application are implemented as software functional modules and sold or used as independent products, they can be stored in a readable storage medium. Based on this understanding, the technical solutions of the embodiments of this application, in essence, or the parts that contribute to the prior art, or all or part of the technical solutions, can be embodied in the form of software products. These software products are stored in a storage medium and include several instructions to cause an electronic device (which may be a mobile phone, personal computer, server, or network device, etc.) or processor to execute all or part of the steps of the methods of the various embodiments of this application. Storage media for storing computer software products include various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0073] In an exemplary embodiment, this application also provides an electronic device. Figure 4 This is a schematic diagram illustrating the composition of an electronic device provided in an embodiment of this application. For example... Figure 4 As shown, the electronic device includes a processor 402, a communication interface 403, and a bus 404. As an example, the electronic device may also include a memory 401.

[0074] Processor 402 may implement or execute various exemplary logic blocks, modules, and circuits described in conjunction with the disclosure of this application. Processor 402 may be a central processing unit, a general-purpose processor, a digital signal processor, an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or other programmable logic devices, transistor logic devices, hardware components, or any combination thereof. It may implement or execute various exemplary logic blocks, modules, and circuits described in conjunction with the disclosure of this application. Processor 402 may also be a combination that implements computing functions, such as including one or more microprocessor combinations, a combination of a DSP and a microprocessor, etc.

[0075] Communication interface 403 is used to connect to other devices via a communication network. This communication network can be Ethernet, wireless access network, wireless local area network (WLAN), etc.

[0076] The memory 401 may be a read-only memory (ROM) or other type of static storage device capable of storing static information and instructions, random access memory (RAM) or other type of dynamic storage device capable of storing information and instructions, or electrically erasable programmable read-only memory (EEPROM), disk storage medium or other magnetic storage device, or any other medium capable of carrying or storing desired program code in the form of instructions or data structures and accessible by a computer, but is not limited thereto.

[0077] As one possible implementation, the memory 401 can exist independently of the processor 402. The memory 401 can be connected to the processor 402 via a bus 404 and is used to store instructions or program code. When the processor 402 calls and executes the instructions or program code stored in the memory 401, it can implement the geological body testing method for aquifer gas storage provided in this application embodiment.

[0078] In another possible implementation, the memory 401 can also be integrated with the processor 402.

[0079] Bus 404 can be an extended industry standard architecture (EISA) bus, etc. Bus 404 can be divided into address bus, data bus, control bus, etc. For ease of representation, Figure 4 The bus is represented by a single thick line, but this does not mean that there is only one bus or one type of bus.

[0080] Through the above description of the implementation methods, those skilled in the art can clearly understand that, for the sake of convenience and brevity, only the division of the above functional modules is used as an example. In actual applications, the above functions can be assigned to different functional modules as needed, that is, the internal structure of the aquifer gas storage geological body testing device can be divided into different functional modules to complete all or part of the functions described above.

[0081] In exemplary embodiments, this application also provides a computer-readable storage medium including computer software instructions that, when executed in an electronic device, cause the electronic device to implement the methods described in the above embodiments. The readable storage medium can also be an external storage device of the electronic device, such as a plug-in hard drive, smart media card (SMC), secure digital (SD) card, flash card, etc., equipped on the electronic device. Further, the readable storage medium can include both internal storage units and external storage devices of the electronic device. The readable storage medium is used to store the software instructions and other programs and data required by the electronic device. The readable storage medium can also be used to temporarily store data that has been output or will be output.

[0082] In an exemplary embodiment, this application also provides a computer program product including computer instructions that, when executed on an electronic device, cause the electronic device to perform the methods described in the above method embodiments.

[0083] In the above embodiments, implementation can be achieved, in whole or in part, through software, hardware, firmware, or any combination thereof. When implemented using software programs, implementation can be, in whole or in part, in the form of a computer program product. This computer program product includes one or more computer-executable instructions. When these computer-executable instructions are loaded and executed on a computer, all or part of the flow or function according to the embodiments of this application is generated. The computer can be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer-executable instructions can be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another. For example, computer-executable instructions can be transmitted from one website, computer, server, or data center to another website, computer, server, or data center via wired (e.g., coaxial cable, fiber optic, digital subscriber line (DSL)) or wireless (e.g., infrared, Bluetooth, microwave, etc.) means. The computer-readable storage medium can be any available medium accessible to a computer or a data storage device containing one or more servers, data centers, etc., that can be integrated with the medium. The available medium can be a magnetic medium (e.g., floppy disk, hard disk, magnetic tape) or an optical medium (e.g., DVD), etc.

[0084] Although this application has been described herein in conjunction with various embodiments, those skilled in the art, by reviewing the accompanying drawings, disclosure, and appended claims, will understand and implement other variations of the disclosed embodiments in carrying out the claimed application. In the claims, the word "comprising" does not exclude other components or steps, and "a" or "an" does not exclude multiple instances. A single processor or other unit can implement several functions listed in the claims. While different dependent claims may recite certain measures, this does not mean that these measures cannot be combined to produce good results.

[0085] Although this application has been described in conjunction with specific features and embodiments, it is obvious that various modifications and combinations can be made thereto without departing from the spirit and scope of this application. Accordingly, this specification and drawings are merely exemplary illustrations of this application as defined by the appended claims, and are considered to cover any and all modifications, variations, combinations, or equivalents within the scope of this application. Clearly, those skilled in the art can make various alterations and modifications to this application without departing from the spirit and scope of this application. Thus, if such modifications and modifications of this application fall within the scope of the claims of this application and their equivalents, this application is also intended to include such modifications and modifications.

[0086] The above are merely specific embodiments of this application, but the scope of protection of this application is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the scope of the technology disclosed in this application should be included within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A method for testing geological bodies in an aquifer gas reservoir, characterized in that, The method includes: The system acquires location data of injection wells and monitoring wells at the aquifer, a first temperature of the injection fluid injected into the aquifer through the injection well, a second temperature of the aquifer, the effective porosity of the aquifer, and the viscosity of the injection fluid at the first temperature; the injection fluid includes a tracer. Acquire background monitoring data; the background monitoring data refers to the pressure baseline data and temperature baseline data of the aquifer collected by the monitoring well when both the injection well and the monitoring well are in a shut-in state and the pressure and temperature of the aquifer remain stable within a preset time period. Acquire monitoring data during the injection test at the injection well; the monitoring data includes the pumping pressure at the injection well, the start and stop times of the injection operation, and the aquifer pressure curve, caprock overburden permeability pressure curve, and tracer breakthrough curve at the monitoring well; The connectivity of the aquifer is determined by comparing the aquifer pressure curve with the injection operation start and stop time points and the pressure baseline data: if the aquifer pressure curve fluctuates at the injection operation start and stop time points, deviating from the pressure baseline data, the aquifer is determined to be connected; if the aquifer pressure curve remains consistent with the pressure baseline data at the injection operation start and stop time points, the aquifer is determined to be disconnected. The capping layer sealing performance is determined by comparing the permeation pressure curve on the capping layer with the injection operation start and stop time points and the pressure baseline data: if the permeation pressure curve on the capping layer remains consistent with the pressure baseline data at the injection operation start and stop time points, the capping layer is confirmed to be sealed; if the permeation pressure curve on the capping layer shows fluctuations that deviate from the pressure baseline data at the injection operation start and stop time points, the capping layer is confirmed to be leaking. Based on the tracer breakthrough curve, the first temperature, the second temperature, the effective porosity, the hydrostatic pressure difference between the injection well and the monitoring well, and the water-heat-mass coupled solute transport analytical model, the permeability, vertical dispersion, and horizontal dispersion of the aquifer in three dimensions are inverted.

2. The method according to claim 1, characterized in that, The method of retrieving the permeability, vertical dispersion, and lateral dispersion of the aquifer in three dimensions based on the tracer breakthrough curve, the first temperature, the second temperature, the effective porosity, the hydrostatic pressure difference between the injection well and the monitoring well, and based on the water-thermal-mass coupled solute transport analytical model, includes: Based on the peak arrival time of the tracer and the location data, the three-dimensional velocity components of the injected fluid in the aquifer are determined, and the composite average velocity is calculated. The longitudinal and lateral spatial variances of the injected fluid are calculated based on the tracer breakthrough curve. Using the synthetic average flow rate, the longitudinal spatial variance, the transverse spatial variance, and the peak arrival time, the initial values ​​of the longitudinal dispersion and the transverse dispersion are inverted through moment analysis, and the temperature-corrected molecular diffusion term determined by the first temperature and the second temperature is explicitly subtracted in the moment analysis. The initial value of the three-dimensional permeability is inverted using the three-dimensional velocity components, the hydrostatic pressure difference, and the fluid viscosity corrected according to the viscosity and the viscosity temperature coefficient. Based on the water-thermal-mass coupling solute transport analytical model that includes the temperature-corrected molecular diffusion term, the three-dimensional flow velocity component is used as the input of the water-thermal-mass coupling solute transport analytical model. Based on the initial values ​​of the longitudinal dispersion, the lateral dispersion, and the three-dimensional permeability, multiple sets of parameter samples are generated and corresponding breakthrough simulation curves are calculated. The peak concentration, peak arrival time, and at least one curve broadening characteristic time of the breakthrough simulation curve and the tracer breakthrough curve are used to construct the fitting error. The Monte Carlo method is used to search for the parameter combination that minimizes the fitting error, and the optimized longitudinal dispersion, lateral dispersion, and three-dimensional permeability are obtained.

3. The method according to claim 2, characterized in that, The analytical model for solute transport based on water-heat-mass coupling satisfies the following relationship: ； in, Indicates solute concentration; Indicates time; , , Represents three-dimensional spatial coordinates; Indicates the delay factor; Indicates that the injected fluid is in Diffusion coefficient in the direction of travel; Indicates that the injected fluid is in Diffusion coefficient in the direction of travel; Indicates that the injected fluid is in Diffusion coefficient in the direction of travel; Indicates that the injected fluid is in Flow velocity component in the direction; Indicates that the injected fluid is in Flow velocity component in the direction; Indicates that the injected fluid is in Flow velocity component in the direction; Represents source and sink terms; the delay factor Determined by the following formula: ； in, This represents the solid-phase medium volume density of the aquifer; Indicates the linear adsorption partition coefficient of the solute; Indicates effective porosity.

4. The method according to claim 3, characterized in that, The diffusion coefficient satisfies the following relationship: ； in, Indicates lateral dispersion. This represents the temperature-corrected molecular diffusion term; ； in, This represents the diffusion coefficient of the injected fluid at the first temperature; This indicates the viscosity of the injected fluid at the first temperature; This indicates the viscosity of the injected fluid at the second temperature; ； in, This represents the viscosity temperature coefficient.

5. The method according to claim 3, characterized in that, The solute concentration satisfies the following relationship: ； in, Represents the unit impulse response function; Indicates time The varying injection concentration rate; This indicates the upper limit of the simulation duration.

6. The method according to claim 3, characterized in that, Longitudinal and lateral dispersion are determined according to the following formulas: ； ； in, Indicates longitudinal dispersion; This represents the synthesized average flow rate; Indicates lateral dispersion; This represents the longitudinal spatial variance; This represents the lateral spatial variance; This indicates the arrival time of the peak value.

7. The method according to claim 3, characterized in that, The three-dimensional permeability is inverted according to the following formula: ； in, This indicates that the aquifer is in Permeability in the direction of travel; This indicates that the aquifer is in Permeability in the direction; This indicates that the aquifer is in Permeability in the direction of travel; This represents the hydrostatic pressure difference of the fluid; This indicates that the injection well and the monitoring well are in Coordinate difference in direction; This indicates that the injection well and the monitoring well are in Coordinate difference in direction; This indicates that the injection well and the monitoring well are in The difference in coordinates along the direction.

8. The method according to claim 3, characterized in that, The multiple sets of parameter samples are generated according to the following formula: ； in, No. Parameters in the group parameter sample ; Indicates parameters The prior minimum value; Indicates parameters The prior maximum value; This represents a random number that follows a triangular distribution. The fitting error satisfies the following relationship: ； in, Indicates the first Group parameter sample parameters The corresponding fitting error; Represents the normalization coefficient; This represents the theoretical peak concentration of the breakthrough simulation curve; This indicates the peak concentration of the tracer in the breakthrough curve; This indicates the theoretical peak arrival time determined based on the breakthrough simulation curve; This indicates the peak arrival time determined based on the tracer breakthrough curve; This represents the broadening characteristic time of the theoretical curve extracted from the breakthrough simulation curve; This represents the curve broadening characteristic time extracted from the tracer breakthrough curve.

9. An electronic device, characterized in that, include: Memory and processor; The memory stores instructions that the processor can execute; When the processor is configured to execute the instructions, the electronic device performs the method as described in any one of claims 1-8.

10. A computer program product, characterized in that, include: Computer instructions; When the computer instructions are executed in an electronic device, the electronic device causes the electronic device to perform the method as described in any one of claims 1-8.

Images