Method and device for determining site-specific efficiency factors for saline aquifer carbon dioxide storage

Abstract

This disclosure provides a method and apparatus for determining the carbon dioxide sequestration efficiency factor (CSP) of a site-level saline aquifer, belonging to the field of carbon dioxide sequestration technology. The method for determining the CSP efficiency factor of a site-level saline aquifer includes: obtaining the geological parameter coefficients, carbon dioxide gravity number, carbon dioxide mobility ratio, and bound water saturation of the saline aquifer in the target area; determining a CSP efficiency factor calculation model applicable to the target area based on the carbon dioxide gravity number of the saline aquifer; and inputting the geological parameter coefficients, carbon dioxide gravity number, carbon dioxide mobility ratio, and bound water saturation of the saline aquifer into the CSP efficiency factor calculation model to determine the CSP efficiency factor of the target area. By determining a CSP efficiency factor calculation model applicable to the target area, the embodiments of this disclosure can accurately calculate the CSP efficiency factor of the target area, providing a theoretical basis for the planning and implementation of carbon dioxide sequestration projects.

Description

Technical Field

[0001] This disclosure relates to the field of carbon dioxide sequestration technology, and specifically to a method and apparatus for determining the efficiency factor of carbon dioxide sequestration in a site-level saline aquifer. Background Technology

[0002] Site-level saline aquifer sequestration sites refer to specific locations within particular geological structures used for carbon dioxide (CO2) sequestration. These sites typically consist of geological elements such as traps, reservoirs, and caprocks, forming a complete sequestration unit. Research on these sites requires relatively high precision, generally within a range of 10–20 km. 2 Within this scope, as our understanding of CO2 sequestration potential deepens, research focus has gradually shifted from macro-scale to more refined site and infusion scales, aiming to more accurately assess the effectiveness of CO2 sequestration.

[0003] Due to the complexity of CO2 geological sequestration mechanisms in saline aquifers and the numerous influencing factors, the sequestration efficiency factor is widely used as a key indicator to measure the combined impact of multiple factors. Currently, mainstream research methods for sequestration efficiency factors, such as those proposed by the U.S. Department of Energy (DOE), the Carbon Capture and Storage Leaders Forum (CSLF), the U.S. Geological Survey (USGS), the International Energy Agency (IEA), and various methods proposed by academia, typically rely on numerical simulations. Parameter values ​​are assigned using a random normal distribution, and the sequestration efficiency factor values ​​under different confidence intervals are ultimately obtained through numerical simulations.

[0004] However, existing methods for studying the storage efficiency factor have certain limitations at the site scale. On the one hand, the uncertainty of parameter values ​​at the site scale leads to significant differences in calculation results among different methods, making it difficult to provide accurate and consistent predictions. On the other hand, the application scope of different methods is not clearly defined, making it difficult to directly apply model predictions to specific sites, resulting in limited applicability. Summary of the Invention

[0005] The purpose of this disclosure is to provide a method and apparatus for determining the carbon dioxide sequestration efficiency factor of a site-level saline aquifer, aiming to partially or completely solve the technical problems in the background art.

[0006] To achieve the above objectives, in a first aspect, embodiments of this disclosure provide a method for determining the carbon dioxide sequestration efficiency factor of a site-level saline aquifer, comprising: acquiring the geological parameter coefficients, carbon dioxide gravity number, carbon dioxide mobility ratio, and bound water saturation of the saline aquifer in the target area; determining a sequestration efficiency factor calculation model applicable to the target area based on the carbon dioxide gravity number of the saline aquifer; and determining the sequestration efficiency factor of the target area based on the geological parameter coefficients, carbon dioxide gravity number, carbon dioxide mobility ratio, bound water saturation, and the sequestration efficiency factor calculation model of the saline aquifer.

[0007] In some embodiments, based on the carbon dioxide gravity number of the saline aquifer, a storage efficiency factor calculation model applicable to the target area is determined, including: when the carbon dioxide gravity number of the saline aquifer is less than a first threshold, the storage efficiency factor calculation model applicable to the target area is determined as a first storage efficiency factor calculation model; when the carbon dioxide gravity number of the saline aquifer is greater than or equal to the first threshold, the storage efficiency factor calculation model applicable to the target area is determined as a second storage efficiency factor calculation model.

[0008] In some embodiments, the storage efficiency factor of the target area is determined based on the geological parameter coefficients of the saline aquifer, the carbon dioxide gravity number, the carbon dioxide mobility ratio, the bound water saturation, and a storage efficiency factor calculation model. This includes: determining the storage efficiency factor based on the geological parameter coefficients of the saline aquifer, the carbon dioxide mobility ratio, the bound water saturation, and a first storage efficiency factor calculation model; wherein, the first storage efficiency factor calculation model is: In the formula, E B For the storage efficiency factor, E geol Geological parameter coefficients, S w For bound water saturation, (1-S w ) av λ represents the average gas saturation and λ represents the carbon dioxide mobility ratio.

[0009] In some embodiments, the storage efficiency factor of the target area is determined based on the geological parameter coefficients of the saline aquifer, the carbon dioxide gravity number, the carbon dioxide mobility ratio, the bound water saturation, and a storage efficiency factor calculation model. This includes: determining the storage efficiency factor based on the geological parameter coefficients of the saline aquifer, the carbon dioxide gravity number, the carbon dioxide mobility ratio, the bound water saturation, and a second storage efficiency factor calculation model; wherein, the second storage efficiency factor calculation model is:

[0010] In the formula, E B For the storage efficiency factor, E geol Γ is the geological parameter coefficient, Γ is the carbon dioxide gravity number, and S is the gravimetric coefficient. w For bound water saturation, (1-S w ) avλ represents the average gas saturation and λ represents the carbon dioxide mobility ratio.

[0011] In some embodiments, the geological parameter coefficient is: E geol =E An / At ×E hn / hg ×E φe / φtot In the formula, E geol Geological parameter coefficients, E An / At The ratio of effective carbon dioxide storage area to total area, E hn / hg The ratio of effective carbon dioxide sequestration thickness to total thickness, E φe / φtot The effective carbon dioxide sequestration porosity is the ratio of the total porosity to the total carbon dioxide porosity; the carbon dioxide gravimetric number is: In the formula, Γ is the gravimetric number of carbon dioxide, Δρ is the difference between the density of brine and the density of supercritical carbon dioxide, k is the permeability of the saline aquifer, g is the acceleration due to gravity, and λ is the gravitational acceleration. b H is the brine mobility, H is the brine layer thickness, and Q is the injection velocity; the carbon dioxide mobility ratio is: In the formula, λ is the carbon dioxide mobility ratio, For supercritical carbon dioxide mobility, λ b This refers to the saline fluidity.

[0012] In some embodiments, supercritical carbon dioxide mobility for: In the formula, The relative permeability of supercritical carbon dioxide, μ CO2 The viscosity of supercritical carbon dioxide; the brine mobility λ b for: In the formula, k rb The relative permeability of saline solution, μ b This refers to the viscosity of the salt water.

[0013] In some embodiments, the method further includes: obtaining the actual storage efficiency factor of the target area; and optimizing the storage efficiency factor calculation model based on the actual storage efficiency factor.

[0014] Secondly, embodiments of this disclosure provide an apparatus for determining the carbon dioxide sequestration efficiency factor of a site-level saline aquifer, comprising: an acquisition unit for acquiring geological parameter coefficients, carbon dioxide gravity number, carbon dioxide mobility ratio, and bound water saturation of the saline aquifer in a target area; a first determination unit for determining a calculation model for the sequestration efficiency factor applicable to the target area based on the carbon dioxide gravity number of the saline aquifer; and a second determination unit for determining the sequestration efficiency factor of the target area based on the geological parameter coefficients, carbon dioxide gravity number, carbon dioxide mobility ratio, bound water saturation, and the sequestration efficiency factor calculation model of the saline aquifer.

[0015] Thirdly, embodiments of this disclosure provide a computer device comprising: a memory configured to store instructions; and a processor configured to retrieve instructions from the memory and, when executing the instructions, to implement the method for determining the site-level saline aquifer carbon dioxide sequestration efficiency factor provided in the first aspect or any embodiment of the first aspect.

[0016] Fourthly, embodiments of this disclosure provide a machine-readable storage medium storing instructions for causing a machine to execute the method for determining the site-level saline aquifer carbon dioxide sequestration efficiency factor provided in the first aspect or any embodiment of the first aspect.

[0017] Through the above technical solutions, the method for determining the site-level saline aquifer carbon dioxide sequestration efficiency factor provided in this disclosure can accurately calculate the sequestration efficiency factor of the target area by determining a sequestration efficiency factor calculation model applicable to the target area, thus providing a theoretical basis for the planning and implementation of carbon dioxide sequestration projects.

[0018] Other features and advantages of the embodiments disclosed herein will be described in detail in the following detailed description section. Attached Figure Description

[0019] The accompanying drawings are provided to further illustrate the embodiments of this disclosure and form part of the specification. They are used together with the following detailed description to explain the embodiments of this disclosure, but do not constitute a limitation thereof. In the drawings:

[0020] Figure 1 This is a flowchart illustrating a method for determining the carbon dioxide sequestration efficiency factor of a site-level saline aquifer according to Embodiment 1 of this disclosure.

[0021] Figure 2 This is a schematic diagram of the interosmosis relationship curve between carbon dioxide and brine according to Embodiment 1 of this disclosure.

[0022] Figure 3 This is a flowchart illustrating a method for determining the efficiency factor of carbon dioxide sequestration in a site-level saline aquifer according to Embodiment 2 of this disclosure.

[0023] Figure 4 This is a schematic diagram of a device for determining the efficiency factor of carbon dioxide sequestration in a site-level saline aquifer, according to an embodiment of this disclosure. Detailed Implementation

[0024] To make the objectives, technical solutions, and advantages of the embodiments of this disclosure clearer, the technical solutions of the embodiments of this disclosure will be clearly and completely described below with reference to the accompanying drawings. It should be understood that the specific embodiments described herein are only for illustration and explanation of the embodiments of this disclosure and are not intended to limit the embodiments of this disclosure. All other embodiments obtained by those skilled in the art based on the embodiments of this disclosure without creative effort are within the scope of protection of this disclosure.

[0025] It should be noted that the acquisition, transmission, storage, use, and processing of data in this disclosed technical solution comply with the relevant provisions of national laws and regulations. In the embodiments of this disclosure, certain existing industry solutions such as software, components, and models may be mentioned. These should be considered exemplary, intended only to illustrate the feasibility of implementing the technical solution of this disclosure, and do not imply that the applicant has already used or necessarily used such solutions.

[0026] Furthermore, if the embodiments of this disclosure involve descriptions such as "first" or "second," these descriptions are for descriptive purposes only and should not be construed as indicating or implying their relative importance or implicitly specifying the number of technical features indicated. Therefore, a feature defined with "first" or "second" may explicitly or implicitly include at least one of those features. Additionally, the technical solutions of the various embodiments can be combined with each other, but this must be based on the ability of those skilled in the art to implement them. If the combination of technical solutions is contradictory or impossible to implement, it should be considered that such a combination of technical solutions does not exist and is not within the scope of protection claimed by this disclosure.

[0027] Figure 1 This is a flowchart illustrating a method for determining the carbon dioxide sequestration efficiency factor of a site-level saline aquifer according to Embodiment 1 of this disclosure.

[0028] like Figure 1 As shown, the method for determining the carbon dioxide sequestration efficiency factor of the site-level saline aquifer may include steps S101 to S103.

[0029] Step S101: Obtain the geological parameter coefficients, carbon dioxide gravity number, carbon dioxide mobility ratio, and bound water saturation of the saline aquifer in the target area.

[0030] The site-level saline aquifer sequestration efficiency factor is influenced by a combination of geological, engineering, and fluid factors. This disclosure mainly considers the influence of three aspects: the influence of geological factors, the influence of macroscopic displacement effects, and the influence of microscopic displacement efficiency.

[0031] Regarding the influence of geological factors, since geological parameter coefficients can represent the overall geological characteristics of the saline aquifer, they are used to partially quantify the impact of geological factors on the sequestration efficiency factor. Geological parameter coefficient E geol Characterized as: E geol =E An / At ×E hn / hg ×E φe / φtot , of which E An / At The ratio of effective carbon dioxide storage area to total area, E hn / hg The ratio of effective carbon dioxide sequestration thickness to total thickness, E φe / φtot It is the ratio of effective carbon dioxide sequestration porosity to total porosity. That is, the geological parameter coefficient is the product of the ratio of effective carbon dioxide sequestration area to total area, the ratio of effective carbon dioxide sequestration thickness to total thickness, and the ratio of effective carbon dioxide sequestration porosity to total porosity.

[0032] The impact on macroscopic displacement efficiency includes the combined effects of heterogeneity within the plume, buoyancy, and the differences in density and viscosity between carbon dioxide and the saline water layer. Within the plume, the macroscopic displacement efficiency E... H Defined as the ratio of the volume of carbon dioxide to the total reservoir volume at the maximum plume radius, it is related to the gravity number. Since the carbon dioxide gravity number can represent the vertical transport capacity of carbon dioxide in a saline aquifer, the carbon dioxide mobility ratio can represent the lateral transport capacity of carbon dioxide in a saline aquifer, and the bound water saturation can represent the effective storage space of carbon dioxide in a saline aquifer, the carbon dioxide gravity number, carbon dioxide mobility ratio, and bound water saturation can be used together to partially quantify the impact of macroscopic displacement efficiency on the storage efficiency factor.

[0033] It is understandable that when the gravity number is small, viscous forces dominate, and lateral transport of carbon dioxide occurs at the bottom of the saline layer, thereby improving the storage efficiency. When the gravity number is large, the flow of carbon dioxide mainly occurs at the top of the saline layer, resulting in a lower storage efficiency factor. Therefore, when calculating the macroscopic displacement efficiency, it is necessary to further distinguish based on the gravity number.

[0034] Specifically, when the gravitational number is less than 0.5, the macroscopic displacement efficiency E H for: In the formula, (1-S w ) av For average gas saturation, S w For bound water saturation, (1-S w ) av Let be the average gas saturation, and λ be the carbon dioxide mobility ratio. When the gravitational number is greater than or equal to 0.5, the macroscopic displacement efficiency E... H for: In the formula, (1-S w )av For average water saturation, S w For bound water saturation, (1-S w ) av Γ represents the average gas saturation, Γ represents the carbon dioxide gravity number, and λ represents the carbon dioxide mobility ratio.

[0035] The impact on micro-displacement efficiency, due to the micro-displacement efficiency E W On the plume scale: E W =(1-S w ) av ≈1-S w In the formula, (1-S w ) av For average gas saturation, S w The bound water saturation is used to represent the effective storage space for carbon dioxide in a saline aquifer. Therefore, the bound water saturation can be used to partially quantify the impact of microscopic displacement efficiency on the storage efficiency factor.

[0036] Therefore, this embodiment of the disclosure requires obtaining the geological parameter coefficients, carbon dioxide gravity number, carbon dioxide mobility ratio, and bound water saturation of the saline aquifer in the target area.

[0037] Step S102: Based on the carbon dioxide gravity number of the saline aquifer, determine the appropriate calculation model for the storage efficiency factor in the target area.

[0038] Based on the above geological factors, macroscopic displacement efficiency, and microscopic displacement efficiency, according to formula E B =E geol ×E H ×E W A pre-established model for the sealing efficiency factor is given, where E is the factor for sealing efficiency. geol Due to the influence of geological factors (geological parameter coefficients), E H For the impact of macroscopic displacement efficiency (macroscopic displacement efficiency), E W The impact on micro-displacement efficiency (micro-displacement efficiency).

[0039] In some embodiments, a model for calculating the storage efficiency factor applicable to the target area is determined based on the carbon dioxide gravity number of the saline aquifer, including:

[0040] When the carbon dioxide gravity number of the saline aquifer is less than the first threshold, the first storage efficiency factor calculation model is determined to be the applicable storage efficiency factor calculation model for the target area.

[0041] When the carbon dioxide gravity number of the saline aquifer is greater than or equal to the first threshold, the storage efficiency factor calculation model applicable to the target area is determined to be the second storage efficiency factor calculation model.

[0042] It should be noted that when the carbon dioxide gravity number is low, viscous forces dominate, and lateral migration of carbon dioxide occurs at the bottom of the saline layer, thereby improving the storage efficiency. When the carbon dioxide gravity number is high, the flow of carbon dioxide mainly occurs at the top of the saline layer, resulting in a lower storage efficiency factor. Therefore, the embodiments of this disclosure adjust the storage efficiency factor calculation model according to the range of carbon dioxide gravity number, and different storage efficiency factor calculation models are applicable to different carbon dioxide gravity number ranges.

[0043] Step S103: Based on the geological parameter coefficients of the saline aquifer, the carbon dioxide gravity number, the carbon dioxide mobility ratio, the bound water saturation, and the storage efficiency factor calculation model, determine the storage efficiency factor of the target area.

[0044] In some embodiments, when the carbon dioxide gravity number of the saline aquifer is less than a first threshold, the storage efficiency factor is determined based on the geological parameter coefficients of the saline aquifer, the carbon dioxide mobility ratio, the bound water saturation, and a calculation model of the first storage efficiency factor.

[0045] The first model for calculating the storage efficiency factor (applicable to gravity numbers less than the first threshold) is as follows:

[0046] In the formula, E B For the storage efficiency factor, E geol Geological parameter coefficients, S w For bound water saturation, (1-S w ) av λ represents the average gas saturation and λ represents the carbon dioxide mobility ratio.

[0047] In some embodiments, when the carbon dioxide gravity number of the saline aquifer is greater than or equal to a first threshold, the storage efficiency factor is determined based on the geological parameter coefficients of the saline aquifer, the carbon dioxide gravity number, the carbon dioxide mobility ratio, the bound water saturation, and a second storage efficiency factor calculation model.

[0048] The second sealing efficiency factor calculation model (applicable to gravity numbers greater than or equal to the first threshold) is as follows:

[0049] In the formula, E B For the storage efficiency factor, E geol Γ is the geological parameter coefficient, Γ is the carbon dioxide gravity number, and S is the gravimetric coefficient. w For bound water saturation, (1-S w ) av λ represents the average gas saturation and λ represents the carbon dioxide mobility ratio.

[0050] The following provides a further explanation of the parameter calculation and value selection in step S101 of the present disclosure embodiments.

[0051] In some embodiments, the geological parameters, carbon dioxide gravity number, carbon dioxide mobility ratio, and bound water saturation of the saline aquifer in the target area can be calculated from data of the target area. Before calculating these parameters, it is necessary to obtain the data required for the calculation. Typically, the required data may include single-well data, seismic data, fluid data, and injection data of the target area.

[0052] Specifically, the acquired single-well data may include raw data, logging data, and logging interpretation data from all wells within the target area, as well as other single-well related data. The acquired seismic data may include geophysical prediction data such as target layer amplitude, frequency, and wave impedance. The acquired fluid data may include fluid-related data such as relative permeability, brine viscosity, and bound water saturation within the target area. The acquired injection data may include injection-related data such as carbon dioxide injection rate. It should be noted that, in this embodiment, the injection data primarily considers the carbon dioxide injection rate and temporarily disregards the impact of other engineering injection parameters on the sealing efficiency factor; however, this is not a further limitation of this embodiment.

[0053] In some embodiments, the geological parameter coefficients of the saline aquifer can be calculated using the acquired seismic data and single-well data. The geological parameter coefficients are: E geol =E An / At ×E hn / hg ×E φe / φtot The ratio of effective carbon dioxide storage area to total area, E. An / At The value of is assigned using seismic properties that characterize the planar distribution of sand and mud. The ratio E of effective carbon dioxide sequestration thickness to total thickness is... hn / hg The value of this parameter is determined by a weighted average of the sand-to-sludge ratio data from all single wells within the target area. The value of E represents the ratio of effective carbon dioxide storage porosity to total porosity. φe / φtot The value of this parameter is determined by a weighted average of the logging interpretation data of all single wells within the target area.

[0054] The weighted average method in the above value selection process can be: In the formula, P is the final parameter value, n is the total number of wells in the target area, and P i For the parameter values ​​of the i-th well, h i Let be the thickness value of the i-th well.

[0055] In some embodiments, the gravitational number of carbon dioxide is: In the formula, Γ is the gravimetric number of carbon dioxide, Δρ is the difference between the density of brine and the density of supercritical carbon dioxide, k is the permeability of the saline aquifer, g is the acceleration due to gravity, and λ is the gravitational acceleration. b H is the brine mobility, H is the brine layer thickness, and Q is the carbon dioxide injection rate.

[0056] The relevant parameters for calculating the gravimetric number of carbon dioxide and the methods for determining their values ​​are shown in Table 1 below.

[0057] Table 1. Relevant parameters and methods for calculating the gravimetric number of carbon dioxide.

[0058]

[0059] In some embodiments, the carbon dioxide mobility ratio is: In the formula, λ is the carbon dioxide mobility ratio, For supercritical carbon dioxide mobility, λ b This refers to the brine mobility. Among them, the supercritical carbon dioxide mobility... for: In the formula, The relative permeability of supercritical carbon dioxide, The viscosity is the supercritical carbon dioxide viscosity. The brine mobility λ b for: In the formula, k rb The relative permeability of saline solution, μ b This refers to the viscosity of the brine. It should be noted that supercritical carbon dioxide refers to the state of carbon dioxide when it is at or above its critical temperature and critical pressure.

[0060] The values ​​of supercritical carbon dioxide mobility and brine mobility can be further determined based on the relative permeability curves of carbon dioxide and brine in the actual target area. If there are no relevant relative permeability curves in the target area, the relative permeability curves given in this disclosure can be used for calculation (see reference). Figure 2 ).like Figure 2 As shown, the red curve represents the relative permeability of the brine, and the blue curve represents the relative permeability of carbon dioxide. With increasing water saturation, the brine content increases, leading to an increase in the relative permeability of the brine; simultaneously, the relative permeability of carbon dioxide decreases.

[0061] In some embodiments, the value of bound water saturation can be determined during actual engineering injection monitoring. If the bound water saturation value is not measured in the saline aquifer of the target area, a value of 0.5 is given as an approximate value for the parameter. The average gas saturation value is also obtained through monitoring data, and the average water saturation value actually reflects the bound water saturation value.

[0062] For example, taking the calculation of the carbon dioxide sequestration efficiency factor of a site-level saline aquifer as an example, the specific steps of the method for determining the carbon dioxide sequestration efficiency factor of a site-level saline aquifer provided in this disclosure embodiment will be described in general. The target area of ​​site-level A is 8.5 km². 2 .

[0063] Step 1: Based on the seismic data and single-well data of the target area, calculate the ratio E of the effective carbon dioxide storage area to the total area. An / At Given a value of 0.77, calculate the ratio E of the effective carbon dioxide sequestration thickness to the total thickness. hn / hg The value is 0.74; calculate the ratio E of effective carbon dioxide storage porosity to total porosity. φe / φtot It is 0.71. Based on the geological parameter coefficient E... geol Definition E geol =E An / At ×E hn / hg ×E φe / φtot By substituting relevant parameters, the geological parameter coefficient E is obtained. geol It is 0.41.

[0064] Step 2: Based on the relevant parameter values ​​in the target area that are related to the calculation of the carbon dioxide gravimetric number, and according to the definition of the carbon dioxide gravimetric number... The carbon dioxide gravity number of the target area was determined to be 0.43 by substituting relevant parameters. The relevant parameter values ​​for the target area are shown in Table 2 below.

[0065] Table 2. Target Area Relevant Parameter Values

[0066] parameter Value unit saline aquifer temperature 75 ℃ saline aquifer pressure 24 MPa <![CDATA[Supercritical CO2 density]]> 720 <![CDATA[kg·m -3 ]]> Salt water density 1043 <![CDATA[kg·m -3 ]]> gravitational acceleration 10 <![CDATA[m·s -2 ]]> saline aquifer permeability <![CDATA[101×10 -3 ]]> <![CDATA[μm 2 ]]> Salt water viscosity 0.69 mPa·s <![CDATA[Viscosity of Supercritical CO2]]> 0.08 mPa·s Saltwater layer thickness 20 m <![CDATA[CO2 injection rate]]> <![CDATA[2.4×10 6 ]]> <![CDATA[m 3 ·d -1 ]]>

[0067] Step 3: Since the carbon dioxide gravity number in the target area is 0.43 < 0.5, the first storage efficiency factor calculation model is selected.

[0068] Step 4: Calculate the carbon dioxide mobility ratio in the target area. (Reference) Figure 2 As shown, when the water saturation is 0.57, the relative permeability of both carbon dioxide and brine is 0.05. Based on the data values ​​of brine viscosity and supercritical carbon dioxide viscosity in Table 2, and the definition of carbon dioxide mobility ratio λ, The carbon dioxide mobility ratio λ in the target area was found to be 8.6.

[0069] It should be noted that the execution order of the first, second, and fourth steps in this embodiment is not limited, except that the third step must be performed after the second step. The first four steps in this embodiment must be completed before the fifth step.

[0070] Step 5: Calculate the sealing efficiency factor for the target area. Based on the calculated geological parameter coefficient E... geol The carbon dioxide mobility ratio λ is taken as the bound water saturation of 0.46, and the average water saturation (1-S) is determined. w ) av Approximately 1-S w The value is 0.54, based on the formula of the first sealing efficiency factor calculation model. Substituting the relevant parameters, the storage efficiency factor of the target area is calculated to be 0.0139, which is equivalent to 1.39%.

[0071] The method for determining the site-level carbon dioxide sequestration efficiency factor of a saline aquifer, provided in Embodiment 1 of this disclosure, selects a sequestration efficiency factor calculation model suitable for the target area by calculating the carbon dioxide gravimetric number, and then calculates the sequestration efficiency factor suitable for the target area based on the model. This method can calculate the site-level sequestration efficiency factor based on relevant geological and fluid data using a theoretical calculation model without performing specific numerical simulations. This disclosure has advantages such as strong targeting, simple application, and high prediction accuracy, making it more suitable for calculating the site-level saline aquifer sequestration efficiency factor, and also providing important guidance for the feasibility of engineering implementation.

[0072] Figure 3 This is a flowchart illustrating a method for determining the efficiency factor of carbon dioxide sequestration in a site-level saline aquifer, according to Embodiment 2 of this disclosure.

[0073] like Figure 3 As shown, the method for determining the carbon dioxide sequestration efficiency factor of the site-level saline aquifer includes steps S101 to S103, and may also include steps S104 and S105.

[0074] Step S101: Obtain the geological parameter coefficients, carbon dioxide gravity number, carbon dioxide mobility ratio, and bound water saturation of the saline aquifer in the target area.

[0075] Step S102: Based on the carbon dioxide gravity number of the saline aquifer, determine the appropriate calculation model for the storage efficiency factor in the target area.

[0076] Step S103: Based on the geological parameter coefficients of the saline aquifer, the carbon dioxide gravity number, the carbon dioxide mobility ratio, the bound water saturation, and the storage efficiency factor calculation model, determine the storage efficiency factor of the target area.

[0077] Other technical features in steps S101 to S103 can be found in the description of steps S101 and S103 in Embodiment 1, and will not be repeated here.

[0078] Step S104: Obtain the actual sealing efficiency factor of the target area.

[0079] The process of obtaining the actual carbon dioxide sequestration efficiency factor involves calculating the effective carbon dioxide sequestration amount and the theoretical carbon dioxide sequestration amount. The effective carbon dioxide sequestration amount refers to the actual amount of carbon dioxide injected into the underground reservoir, minus any amount that may escape into the atmosphere or other formations. This can be obtained by monitoring the injection process and using geophysical exploration techniques (such as seismic imaging). The theoretical carbon dioxide sequestration amount refers to the maximum possible amount assuming all injected carbon dioxide can be safely stored in the target geological body without loss. This can be predicted using mathematical models based on parameters such as reservoir physical properties (e.g., porosity, permeability), reservoir volume, and the phase behavior of carbon dioxide (e.g., solubility, density). The actual sequestration efficiency factor is the ratio between the effective carbon dioxide sequestration amount and the theoretical carbon dioxide sequestration amount; this ratio reflects the actual efficiency of the sequestration process.

[0080] By obtaining the actual storage efficiency factor, we can provide data support for subsequent model optimization, making the storage efficiency factor calculation model closer to actual conditions.

[0081] Step S105: Correct the sealing efficiency factor based on the actual sealing efficiency factor.

[0082] After obtaining the actual storage efficiency factor, the model can be adjusted and optimized by comparing the deviation between the predicted and actual storage efficiency factors, thereby improving its predictive accuracy. The optimized model can be used to correct the storage efficiency factor, more effectively guiding the planning and execution of storage projects, reducing resource waste, and improving storage efficiency. Furthermore, optimizing the model helps ensure that storage projects are both economical and environmentally friendly, achieving more sustainable carbon dioxide storage by reducing unnecessary storage costs and environmental risks.

[0083] The method for determining the site-level saline aquifer carbon dioxide sequestration efficiency factor provided in Embodiment 2 of this disclosure, based on Embodiment 1, further enhances the accuracy of the sequestration efficiency factor calculation by introducing the measurement of the actual sequestration efficiency factor and model optimization. This helps to improve the planning and execution capabilities of carbon dioxide sequestration projects, ensuring that they are both economical and environmentally friendly.

[0084] The other beneficial effects of the method for determining the carbon dioxide sequestration efficiency factor of a site-level saline aquifer provided in Embodiment 2 of this disclosure are the same as those of the method for determining the carbon dioxide sequestration efficiency factor of a site-level saline aquifer provided in Embodiment 1 above, and will not be repeated here.

[0085] Figure 4 This is a schematic diagram of a device for determining the efficiency factor of carbon dioxide sequestration in a site-level saline aquifer, according to an embodiment of this disclosure.

[0086] like Figure 4 As shown, the device 100 for determining the carbon dioxide sequestration efficiency factor of the site-level saline aquifer includes an acquisition unit 110, a first determination unit 120, and a second determination unit 130.

[0087] The acquisition unit 110 is used to acquire the geological parameter coefficients, carbon dioxide gravity number, carbon dioxide mobility ratio and bound water saturation of the saline aquifer in the target area.

[0088] The first determining unit 120 is used to determine a calculation model for the storage efficiency factor applicable to the target area based on the carbon dioxide gravity number of the saline aquifer.

[0089] The second determining unit 130 is used to determine the storage efficiency factor of the target area based on the geological parameter coefficients of the saline aquifer, the carbon dioxide gravity number, the carbon dioxide mobility ratio, the bound water saturation, and the storage efficiency factor calculation model.

[0090] The apparatus for determining the carbon dioxide sequestration efficiency factor of a site-level saline aquifer provided in this disclosure adopts the method for determining the carbon dioxide sequestration efficiency factor of a site-level saline aquifer provided in the above embodiments, and can solve the technical problems mentioned in the background art.

[0091] The beneficial effects of the apparatus for determining the carbon dioxide sequestration efficiency factor of a site-level saline aquifer provided in this disclosure are the same as those of the method for determining the carbon dioxide sequestration efficiency factor of a site-level saline aquifer provided in the above embodiments. Furthermore, the other technical features of the apparatus for determining the carbon dioxide sequestration efficiency factor of a site-level saline aquifer are the same as those disclosed in the method for determining the carbon dioxide sequestration efficiency factor of a site-level saline aquifer, and will not be repeated here.

[0092] This disclosure also provides a computer device, comprising: a memory configured to store instructions; and a processor configured to retrieve instructions from the memory and, when executing the instructions, to implement the method for determining the site-level saline aquifer carbon dioxide sequestration efficiency factor provided in any of the above embodiments.

[0093] The beneficial effects of the computer equipment provided in this embodiment are the same as those of the method for determining the site-level saline aquifer carbon dioxide sequestration efficiency factor provided in the above embodiments, and will not be repeated here.

[0094] This disclosure also provides a machine-readable storage medium storing instructions for causing a machine to execute the method for determining the site-level saline aquifer carbon dioxide sequestration efficiency factor provided in any of the above embodiments.

[0095] The beneficial effects of the machine-readable storage medium provided in this disclosure are the same as those of the method for determining the site-level saline aquifer carbon dioxide sequestration efficiency factor provided in the above embodiments, and will not be repeated here.

[0096] Those skilled in the art will understand that embodiments of this disclosure can be provided as methods, systems, or computer program products. Therefore, this disclosure can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, this disclosure can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0097] This disclosure is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of this disclosure. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create a machine for implementing the flowchart illustrations. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0098] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0099] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0100] In a typical configuration, a computing device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0101] Memory may include non-persistent memory in computer-readable media, such as random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash RAM. Memory is an example of computer-readable media.

[0102] Computer-readable media includes both permanent and non-permanent, removable and non-removable media that can store information using any method or technology. Information can be computer-readable instructions, data structures, modules of programs, or other data. Examples of computer storage media include, but are not limited to, phase-change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, CD-ROM, digital versatile optical disc (DVD) or other optical storage, magnetic tape, disk storage or other magnetic storage devices, or any other non-transferable medium that can be used to store information accessible by a computing device. As defined herein, computer-readable media does not include transient computer-readable media, such as modulated data signals and carrier waves.

[0103] It should also be noted that the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, 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 process, method, article, or apparatus. Unless otherwise specified, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes that element.

[0104] The above are merely embodiments of this disclosure and are not intended to limit the scope of this disclosure. Various modifications and variations can be made to this disclosure by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this disclosure should be included within the scope of the claims of this disclosure.

Claims

1. A method for determining the carbon dioxide sequestration efficiency factor of a site-level saline aquifer, characterized in that, include: Obtain geological parameters, carbon dioxide gravity number, carbon dioxide mobility ratio, and bound water saturation of the saline aquifer in the target area; Based on the carbon dioxide gravity number of the saline aquifer, a calculation model for the storage efficiency factor applicable to the target area is determined; Based on the geological parameter coefficients of the saline aquifer, the carbon dioxide gravity number, the carbon dioxide mobility ratio, the bound water saturation, and the calculation model of the sequestration efficiency factor, the sequestration efficiency factor of the target area is determined.

2. The method of claim 1, wherein the site-scale saline aquifer CO2 storage efficiency factor is determined by: The calculation model for determining the storage efficiency factor applicable to the target area based on the carbon dioxide gravity number of the saline aquifer includes: When the carbon dioxide gravity number of the saline aquifer is less than the first threshold, the storage efficiency factor calculation model applicable to the target area is determined to be the first storage efficiency factor calculation model. When the carbon dioxide gravity number of the saline aquifer is greater than or equal to the first threshold, the storage efficiency factor calculation model applicable to the target area is determined to be the second storage efficiency factor calculation model.

3. The method of claim 2, wherein the site-scale saline aquifer CO2 storage efficiency factor is determined by: The calculation model based on the geological parameter coefficients of the saline aquifer, the carbon dioxide gravity number, the carbon dioxide mobility ratio, the bound water saturation, and the storage efficiency factor determines the storage efficiency factor of the target area, including: Based on the geological parameter coefficients of the saline aquifer, the carbon dioxide mobility ratio, the bound water saturation, and the calculation model of the first storage efficiency factor, the storage efficiency factor is determined. The calculation model for the first sealing efficiency factor is as follows: In the formula, E B For the storage efficiency factor, E geol Geological parameter coefficients, S w For bound water saturation, (1-S w ) av λ represents the average gas saturation and λ represents the carbon dioxide mobility ratio.

4. The method of claim 2, wherein the site-scale saline aquifer CO2 storage efficiency factor is determined by: The calculation model based on the geological parameter coefficients of the saline aquifer, the carbon dioxide gravity number, the carbon dioxide mobility ratio, the bound water saturation, and the storage efficiency factor determines the storage efficiency factor of the target area, including: Based on the geological parameter coefficients of the saline aquifer, the carbon dioxide gravity number, the carbon dioxide mobility ratio, the bound water saturation, and the calculation model of the second storage efficiency factor, the storage efficiency factor is determined. The calculation model for the second sealing efficiency factor is as follows: In the formula, E B For the storage efficiency factor, E geol Γ is the geological parameter coefficient, Γ is the carbon dioxide gravity number, and S is the gravimetric coefficient. w For bound water saturation, (1-S w ) av λ represents the average gas saturation and λ represents the carbon dioxide mobility ratio.

5. The method for determining the site-level carbon dioxide sequestration efficiency factor of a saline aquifer according to any one of claims 1-4, characterized in that, The geological parameter coefficient is: E geol =E An / At ×E hn / hg ×E φe / φtot In the formula, E geol Geological parameter coefficients, E An / At The ratio of effective carbon dioxide storage area to total area, E hn / hg The ratio of effective carbon dioxide sequestration thickness to total thickness, E φe / φtot It is the ratio of effective carbon dioxide storage porosity to total porosity; The gravitational number of carbon dioxide is: In the formula, Γ is the gravimetric number of carbon dioxide, Δρ is the difference between the density of brine and the density of supercritical carbon dioxide, k is the permeability of the saline aquifer, g is the acceleration due to gravity, and λ is the gravitational acceleration. b H is the brine mobility, H is the brine layer thickness, and Q is the injection velocity; The carbon dioxide mobility ratio is: In the formula, λ is the carbon dioxide mobility ratio, For supercritical carbon dioxide mobility, λ b This refers to the saline fluidity.

6. The method for determining the carbon dioxide sequestration efficiency factor of a site-level saline aquifer according to claim 5, characterized in that, the supercritical carbon dioxide flow rate is: wherein, is the supercritical carbon dioxide relative permeability, is the supercritical carbon dioxide viscosity; The brine mobility λ b is: where k rb is the brine relative permeability, μ b is the brine viscosity.

7. The method of claim 1, wherein: Also includes: Obtain the actual sealing efficiency factor of the target area; Based on the actual storage efficiency factor, the calculation model of the storage efficiency factor is optimized.

8. A device for determining the carbon dioxide sequestration efficiency factor of a site-level saline aquifer, characterized in that, include: The acquisition unit is used to acquire geological parameter coefficients, carbon dioxide gravity number, carbon dioxide mobility ratio, and bound water saturation of the saline aquifer in the target area. The first determining unit is used to determine a storage efficiency factor calculation model applicable to the target area based on the carbon dioxide gravity number of the saline aquifer. The second determining unit is used to determine the storage efficiency factor of the target area based on the geological parameter coefficients of the saline aquifer, the carbon dioxide gravity number, the carbon dioxide mobility ratio, the bound water saturation, and the storage efficiency factor calculation model.

9. A computer device, comprising: include: The memory is configured to store instructions; as well as The processor is configured to retrieve the instructions from the memory and, when executing the instructions, to implement the method for determining the site-level saline aquifer carbon dioxide sequestration efficiency factor according to any one of claims 1 to 7.

10. A machine-readable storage medium, characterized in that, The machine-readable storage medium stores instructions for causing the machine to perform the method for determining the site-level saline aquifer carbon dioxide sequestration efficiency factor as described in any one of claims 1 to 7.

Images