Method for evaluating the safety of co2 storage in saline aquifers based on dimensionless numbers ca, gr and bo

Abstract

The application discloses a kind of based on dimensionless number Ca, Gr and Bo's saltwater layer CO2 storage safety evaluation method, it is related to oil reservoir and underground water numerical simulation technical field, comprising: (1) the establishment of saltwater layer CO2 storage geological model, including target reservoir, interlayer and overlying saltwater layer;(2) crack model and monitoring point are established in interlayer;(3) the geological model is initialized;(4) research CO2 leakage dynamics, calculate force and dimensionless number;(5) calculate storage safety factor;(6) change initialization condition, repeat steps (4)-(5), obtain three groups of dimensionless number and storage safety factor;(7) establish the linear regression model of storage safety factor and dimensionless number, calculate saltwater layer CO2 storage safety factor.The method of the application can adapt to various schemes under different properties reservoir, different pressure gradient, the main driving force and leakage law in the process of CO2 leakage through crack are clarified, with strong adaptability, high calculation efficiency advantage.

Description

TECHNICAL FIELD

[0001] The application relates to the technical field of oil reservoir and underground water numerical simulation, and particularly relates to a brine layer CO2 storage safety evaluation method based on dimensionless numbers Ca, Gr and Bo. BACKGROUND

[0002] The emission and treatment of greenhouse gases are becoming a topic of increasing concern for countries, and China has clearly proposed the "double carbon target" of striving to reach the peak of CO2 emissions before 2030 and striving to achieve carbon neutrality before 2060. Separating CO2 from industrial emission sources, transporting it and utilizing or geologically storing it (CCUS) is one of the effective ways to achieve "net zero" CO2 emissions. Potential CO2 storage geological bodies include abandoned oil and gas reservoirs, coal seams, deep brine layers and the like, and the storage capacity of the brine layer is hundreds of times that of the oil and gas field. China has a considerable deep brine layer block in the Songliao Basin, the Tarim Basin, the Bohai Bay Basin, the Subei Basin and the Ordos Basin and the like, and the CO2 storage potential is about 2420 billion tons.

[0003] After CO2 is injected into the brine layer, leakage along abandoned wells, faults / fractures and the like is a major challenge of the project. Although there are many studies on CO2 leakage at present, various CO2 leakage risk evaluation methods have been proposed, but the results do not consider the risk evaluation of in-situ leakage in the underground, the statistics of the CO2 leakage rate is not real-time, and cannot be applied to general fractured brine layers; the definition of the CO2 storage safety factor of the geological body is not clear, and the nature of the force in the leakage process is not considered, and there is a lack of unified understanding of the operation factors and the geological factors. SUMMARY

[0004] To solve the above technical problems, the application discloses a brine layer CO2 storage safety evaluation method based on dimensionless numbers Ca, Gr and Bo, a fractured brine layer CO2 geological model is established, the CO2 leakage rule in the long-term storage process is studied, the relationship between the storage safety factor and the dimensionless number is clarified, and the fast and accurate prediction of the storage safety factor is realized.

[0005] To achieve the above purpose, the application adopts the following technical scheme:

[0006] The application discloses a brine layer CO2 storage safety evaluation method based on dimensionless numbers Ca, Gr and Bo, and the method comprises the following steps:

[0007] (1) a brine layer CO2 storage geological model is established, including a target reservoir, a barrier interlayer and an overlying brine layer;

[0008] (2) a fracture model and a monitoring point are established in the barrier interlayer;

[0009] (3) initializing the geological model;

[0010] (4) researching CO2 leakage dynamics, calculating force and dimensionless number;

[0011] (5) calculating the storage safety factor;

[0012] (6) changing the initialization condition, repeating steps (4) to (5) to obtain three groups of dimensionless numbers and storage safety factors;

[0013] (7) establishing a linear regression model of the storage safety factor and the dimensionless number, and calculating the CO2 storage safety factor of the saline aquifer.

[0014] Optionally, in step (4), the step of researching CO2 leakage dynamics, calculating force and dimensionless number is specifically: researching the dynamics of CO2 plume passing through the fracture and leaking upward by a numerical simulation method, calculating the force and dimensionless number in the CO2 leakage process according to the data of the simulation monitoring point, and the data of the simulation monitoring point includes the density of CO2 and salt water, the viscosity of CO2, the flow velocity of CO2, the capillary pressure at the monitoring point and the permeability.

[0015] Optionally, in step (5), the storage safety factor is the proportion of CO2 leaked from the target reservoir to the overlying saline aquifer existing in the form of dissolution storage and mineralization storage, as shown in formula (1):

[0016]

[0017] In the formula, y is the storage safety factor; M l is the total mass of CO2 leaked from the target reservoir to the overlying saline aquifer; M d is the mass of CO2 leaked in the form of dissolution storage; M m is the mass of CO2 leaked in the form of mineralization storage.

[0018] Optionally, in step (6), the dimensionless number is:

[0019]

[0020]

[0021]

[0022] In the formula, Ca is the capillary number; Gr is the gravity number; Bo is the Bond number; F B , F V , F CP are buoyancy, viscous force and capillary force respectively; ρ w , ρ grespectively, g is the acceleration of gravity, and l is the characteristic length. cr is the characteristic length. is the CO2 viscosity. is the CO2 flow rate, k is the permeability, and p is the pressure. cr is the capillary pressure.

[0023] Optionally, in step (7), the storage safety factor is linearly regressed with the dimensionless number as follows:

[0024] y = A·Ca + B·Bo + C (5)

[0025] wherein A, B and C are regression coefficients.

[0026] The second aspect of the present application proposes a system based on the method of the first aspect, comprising:

[0027] a geological model construction module for establishing a geological model containing fractures and initializing model parameters;

[0028] a calculation module for calculating the dimensionless number and the storage safety factor;

[0029] a storage safety factor regression model construction module for establishing a linear regression model of the storage safety factor and the dimensionless number;

[0030] The establishment of the geological model containing fractures and the initialization of the model parameters comprises: establishing a CO2 storage geological model of a saline aquifer, including a target reservoir, a barrier interlayer and an overlying saline aquifer; establishing a fracture model and a monitoring point in the barrier interlayer; initializing the geological model;

[0031] The calculation of the dimensionless number and the storage safety factor comprises: studying CO2 leakage dynamics, calculating the acting force and the dimensionless number; calculating the storage safety factor; changing the initialization conditions to obtain three groups of dimensionless numbers and storage safety factors;

[0032] The establishment of the linear regression model of the storage safety factor and the dimensionless number comprises: establishing a linear regression model of the storage safety factor and the dimensionless number, and calculating the CO2 storage safety factor of the saline aquifer.

[0033] The third aspect of the present application proposes a computer readable storage medium having a program stored thereon, which, when executed by a processor, implements the steps of the CO2 storage safety evaluation method of the saline aquifer based on the dimensionless numbers Ca, Gr and Bo according to the first aspect of the present application.

[0034] The fourth aspect of the present application provides an electronic device, comprising a memory, a processor and a program stored in the memory and executable on the processor, wherein the processor implements the steps in the method for evaluating the safety of CO2 storage in a saline aquifer based on dimensionless numbers Ca, Gr and Bo according to the first aspect of the present application when executing the program.

[0035] The present application has the advantages that the method for evaluating the safety factor of CO2 storage in a general fractured saline aquifer according to the present application can be applied to various schemes under different reservoir properties and pressure gradients, and the main driving force and leakage rule in the process of CO2 leakage through fractures are determined, so that the method has the advantages of strong adaptability and high calculation efficiency. The method can effectively evaluate the leakage risk of CO2 storage in a fractured saline aquifer, and provide a reference for the site selection of a target area and horizon for CO2 storage in a saline aquifer. BRIEF DESCRIPTION OF DRAWINGS

[0036] Figure 1 The flowchart of the method for evaluating the safety of CO2 storage in a saline aquifer based on dimensionless numbers Ca, Gr and Bo according to the present application is shown in FIG. 1.

[0037] Figure 2 The geological model of CO2 storage in a fractured saline aquifer and the grid division diagram according to the embodiment of the present application are shown in FIG. 2.

[0038] Figure 3 The stress analysis diagram in the process of CO2 leakage along fractures according to the embodiment of the present application is shown in FIG. 3.

[0039] Figure 4 The average dimensionless number in the process of CO2 leakage under different pressure gradients according to the embodiment of the present application is shown in FIG. 4.

[0040] Figure 5 The graph of the change of CO2 leakage rate with time under different pressure gradients according to the embodiment of the present application is shown in FIG. 5.

[0041] Figure 6 The multiple linear regression diagram of the safety factor according to the embodiment of the present application is shown in FIG. 6. DETAILED DESCRIPTION

[0042] The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative efforts fall within the scope of the present application.

[0043] Embodiment 1

[0044] A method for evaluating the safety of CO2 storage in a saline aquifer based on dimensionless numbers Ca, Gr and Bo, as shown inFigure 1 The method comprises the following steps as shown in the figure:

[0045] (1) Establish a geological model of CO2 storage in saline aquifer according to literature data, including target reservoir, interlayer and overlying saline aquifer;

[0046] (2) Establish a fracture model in the interlayer by local grid refinement method to connect the target reservoir and the overlying saline aquifer, and set simulation monitoring points in the fracture area;

[0047] (3) Initialize the geological model according to the actual block and literature research data, set the lower target reservoir as a high-pressure CO2 saturated area, and simulate the conditions after CO2 injection into the saline aquifer;

[0048] (4) Study the dynamics of CO2 plume passing through the fracture and leaking upward by numerical simulation method, and calculate the force and dimensionless number in the CO2 leakage process according to the data of the simulation monitoring points; the data of the simulation monitoring points include the density of CO2 and saline water, the viscosity of CO2, the flow velocity of CO2, the capillary pressure at the monitoring point and the permeability;

[0049] (5) Statistically analyze the leakage rate, cumulative leakage amount and proportion of CO2 in various storage modes during the simulation process, and calculate the corresponding storage safety factor; the storage safety factor is the proportion of CO2 existing in the form of dissolution storage and mineralization storage from the target reservoir to the overlying saline aquifer, as shown in formula (1):

[0050]

[0051] In the formula, y is the storage safety factor; M l is the total mass of CO2 leaked from the target reservoir to the overlying saline aquifer; M d is the mass of CO2 existing in the form of dissolution storage in the leaked CO2; M m is the mass of CO2 existing in the form of mineralization storage in the leaked CO2;

[0052] (6) Change the initial pressure gradient of the target reservoir and the overlying saline aquifer in step (3), repeat steps (4) to (5), and obtain three groups of dimensionless numbers and storage safety factors; the dimensionless numbers are:

[0053]

[0054]

[0055]

[0056] In the formula, Ca is the capillary number; Gr is the gravity number; Bo is the Bond number; F B , F V , FCP respectively buoyancy, viscous force and capillary force; p w , p g respectively salt water and CO2 density; g is the acceleration of gravity; l cr is the characteristic length; is the CO2 viscosity; is the CO2 flow rate; k is the permeability; p cr is the capillary pressure;

[0057] (7) A linear regression model of the sealing safety factor and dimensionless number is established to calculate the CO2 sealing safety factor of the salt water layer; the linear regression model of the sealing safety factor and dimensionless number is as follows:

[0058] y=A·Ca+B·Bo+C (5)

[0059] wherein A, B and C are all regression coefficients.

[0060] Application Example

[0061] A two-dimensional cylindrical coordinate geological model is established to study the migration law of the leaked CO2 from the target reservoir through the fracture, and a safety evaluation method of the CO2 sealing of the salt water layer based on the dimensionless numbers Ca, Gr and Bo is proposed, combined with Figures 2 to 6 , and the specific steps are as follows:

[0062] (1) A two-dimensional cylindrical coordinate geological model is established by using the numerical simulation software CMG-GEM, which is divided into an overlying salt water layer, a barrier interlayer and a target reservoir, as shown in Figure 2 The top depth of the overlying salt water layer is 1615.44 m, the radial length of the model is 1219.20 m, and the total thickness of the model is 60.96 m. Among them, the thicknesses of the overlying salt water layer and the barrier interlayer are both 15.24 m, and the thickness of the target reservoir is 30.48 m, which is divided into 100x45 (r x z) grid numbers in a radial increasing manner;

[0063] (2) A fracture model is established in the barrier interlayer by using the method of local grid densification, which connects the target reservoir and the overlying salt water layer, and a simulation monitoring point A is set in the fracture area, as shown in Figure 2 The fracture porosity is 0.9, the fracture permeability is 2000 mD, and the permeability of the remaining barrier layer is 0 mD;

[0064] (3) The geological model is initialized, the lower target reservoir is set as a high-pressure CO2 saturated area, the initial pressure is set to 20 MPa, the corresponding initial pressure gradient of the reservoir and the overlying salt water layer is 143.65 kPa / m, and the initial CO2 saturation is 0.8. The top pressure of the overlying salt water layer is 15.62 MPa, and the temperature is 60℃. The overlying salt water layer and the barrier interlayer are initially saturated with formation water;

[0065] (4) The dynamics of the CO2 plume passing through the fracture and leaking upward in the target reservoir were simulated using CMG-GEM, and the forces acting on the CO2 during this process were analyzed, such as... Figure 3 As shown, it can be seen that during the CO2 leakage process, the buoyancy generated by the density difference is the main driving force for leakage, while the viscous force and capillary force are the resistances to the leakage process. The capillary number, gravity number, and Bond number can be calculated using formulas (2) to (4). Figure 4 As shown;

[0066]

[0067]

[0068]

[0069] In the formula, Ca is the capillary number; Gr is the gravitational number; Bo is the bond number; F B F V F CP These are buoyancy, viscous force, and capillary force, respectively; ρ w ρ g The densities of saltwater and CO2 are respectively; g is the acceleration due to gravity; l cr The characteristic length; CO2 viscosity; CO2 flow rate; k is permeability; p cr Capillary pressure;

[0070] (5) The simulation time was set to 1000 years, and the CO2 leakage rate during this process was statistically analyzed, such as... Figure 5 As shown, the CO2 leakage rate generally exhibits a gradual decreasing characteristic, which is due to the continuous decrease in reservoir pressure and the weakening of the driving force for leakage. The mass of each occurrence form of the leaked CO2 is statistically analyzed, and the corresponding sequestration safety factor is calculated using formula (1).

[0071]

[0072] In the formula, y is the sealing safety factor; M l M represents the total mass of CO2 leaking from the target reservoir to the overlying saline aquifer. d M represents the mass of the leaked CO2 that is dissolved and sequestered. m The mass of the leaked CO2 that exists in the form of mineralized sequestration;

[0073] (6) change the initial pressure of the target reservoir in step (3), respectively set to 25MPa and 30MPa, so that the initial pressure gradient of the corresponding target reservoir and the overlying saline aquifer is 307.64kPa / m and 471.63kPa / m respectively; repeat steps (4)-(5) to obtain three groups of dimensionless numbers and storage safety factors.

[0074] (7) using a multiple linear regression tool, combined with the above data, a linear regression model of the storage safety factor and the dimensionless number can be established, as shown in formula (5), the regression coefficients A=-0.00553, B=11.73, C=17.81 are obtained, and the linear regression model of the storage safety factor and the dimensionless number is: Figure 6

[0075] y=-0.00553*Ca+11.73*Bo+17.81 (5)

[0076] Using the model, the reservoir safety factor of the saline aquifer for long-term storage of CO2 can be quickly calculated.

[0077] At present, in the safety research of the saline aquifer for storing CO2, various methods for evaluating the risk of CO2 leakage are proposed, but the essence of the force in the leakage process is not considered, and the storage safety factor established cannot well adapt to different geological models, and there is a lack of unified understanding of operation factors and geological factors. The application provides an evaluation method for the storage safety factor of the general fractured saline aquifer for storing CO2, which can adapt to various schemes under different properties of the reservoir and different pressure gradients, clearly defines the main driving force and leakage rule in the leakage process of CO2 through the fractures, has the advantages of strong adaptability and high calculation efficiency, and can effectively evaluate the risk of CO2 leakage in the fractured saline aquifer for storing CO2, thereby providing a reference for site selection of the target area and horizon of the saline aquifer for storing CO2.

[0078] Embodiment 2

[0079] A system based on the method of the first aspect, comprising:

[0080] a geological model construction module for establishing a geological model containing fractures and initializing model parameters;

[0081] a calculation module for calculating dimensionless numbers and storage safety factors;

[0082] a storage safety factor regression model construction module for establishing a linear regression model of the storage safety factor and the dimensionless number;

[0083] the establishment of the geological model containing fractures and the initialization of the model parameters comprises: establishing a geological model of the saline aquifer for storing CO2, including a target reservoir, a barrier interlayer and an overlying saline aquifer; establishing a fracture model and a monitoring point in the barrier interlayer; initializing the geological model; ​

[0084] The calculation of the dimensionless number and the sealing safety coefficient comprises: researching CO2 leakage dynamics, calculating force and dimensionless number; calculating the sealing safety coefficient; changing the initial condition to obtain three groups of dimensionless numbers and sealing safety coefficients;

[0085] The linear regression model of the sealing safety coefficient and the dimensionless number is established, and the CO2 sealing safety coefficient of the saline aquifer is calculated.

[0086] Embodiment 3

[0087] A computer readable storage medium has a program stored thereon, and the program is executed by a processor to implement the steps in the CO2 sealing safety evaluation method of the saline aquifer based on the dimensionless numbers Ca, Gr and Bo according to the first aspect of the present application, comprising:

[0088] (1) Establish a CO2 geological model of the saline aquifer, including a target reservoir, a barrier interlayer and an overlying saline aquifer;

[0089] (2) Establish a fracture model and a monitoring point in the barrier interlayer;

[0090] (3) Initialize the geological model;

[0091] (4) Research CO2 leakage dynamics, calculate force and dimensionless number;

[0092] (5) Calculate the sealing safety coefficient;

[0093] (6) Change the initial condition, repeat steps (4) to (5) to obtain three groups of dimensionless numbers and sealing safety coefficients;

[0094] (7) Establish a linear regression model of the sealing safety coefficient and the dimensionless number, and calculate the CO2 sealing safety coefficient of the saline aquifer.

[0095] The detailed steps of the method implemented by the above program are the same as the CO2 sealing safety evaluation method of the saline aquifer based on the dimensionless numbers Ca, Gr and Bo provided in Embodiment 1, and will not be repeated here.

[0096] Embodiment 4

[0097] An electronic device comprises a memory, a processor and a program stored on the memory and executable on the processor, and the processor implements the steps in the CO2 sealing safety evaluation method of the saline aquifer based on the dimensionless numbers Ca, Gr and Bo according to the first aspect of the present application when executing the program, comprising:

[0098] (1) Establish a CO2 geological model of the saline aquifer, including a target reservoir, a barrier interlayer and an overlying saline aquifer;

[0099] (2) Establishing a fracture model and monitoring points in the interlayer;

[0100] (3) Initializing the geological model;

[0101] (4) Studying CO2 leakage dynamics, calculating force and dimensionless number;

[0102] (5) Calculating the storage safety factor;

[0103] (6) Changing the initialization conditions, repeating steps (4) to (5) to obtain three groups of dimensionless numbers and storage safety factors;

[0104] (7) Establishing a linear regression model of the storage safety factor and the dimensionless number, and calculating the CO2 storage safety factor of the saline aquifer.

[0105] The detailed steps of the method realized by the above procedure are the same as the method for evaluating the CO2 storage safety of the saline aquifer based on the dimensionless numbers Ca, Gr and Bo provided in Embodiment 1, and thus will not be repeated here.

[0106] Those skilled in the art will understand that the embodiments disclosed in the present application can be provided as a method, a system, or a computer program product. Therefore, the present application can take the form of a hardware embodiment, a software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present application can take the form of a computer program product implemented on one or more computer-usable storage media (including, but not limited to, disk storage and optical storage, etc.) containing computer-usable program code.

[0107] The present application is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present disclosure. It should be understood that each flow and / or block in the flowcharts and / or block diagrams, as well as combinations of flows and / or blocks in the flowcharts 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, a special-purpose computer, an embedded processor, or other programmable data processing apparatus to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing apparatus produce a device that implements the functions specified in the flowcharts and / or block diagrams. Figure 1 The functions specified in a flow or multiple flows and / or blocks Figure 1 The functions specified in a flow or multiple flows and / or blocks

[0108] These computer program instructions can also be stored in a computer-readable memory that can direct the computer or other programmable data processing apparatus to work in a specific manner, so that the instructions stored in the computer-readable memory produce a manufactured product including instruction devices that implement the functions specified in the flowcharts and / or block diagrams. Figure 1 The functions specified in a flow or multiple flows and / or blocks Figure 1 The functions specified in a flow or multiple flows and / or blocks

[0109] 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.

[0110] Those skilled in the art will understand that all or part of the processes in the above embodiments can be implemented by a computer program instructing related hardware. The program can be stored in a computer-readable storage medium, and when executed, it can include the processes of the embodiments of the above methods. The storage medium can be a magnetic disk, optical disk, read-only memory (ROM), or random access memory (RAM), etc.

[0111] Of course, the above description is not intended to limit the present invention, and the present invention is not limited to the examples given above. Any changes, modifications, additions or substitutions made by those skilled in the art within the scope of the present invention should also fall within the protection scope of the present invention.

Claims

1. A method for assessing the safety of CO2 sequestration in saline aquifers based on dimensionless numbers Ca, Gr, and Bo, characterized in that, Includes the following steps: (1) Establish a geological model for CO2 sequestration in a saline aquifer, including the target reservoir, interlayers and overlying saline aquifer; (2) Establish crack models and monitoring points in the interlayer; (3) Initialize the geological model; (4) Study the dynamics of CO2 leakage and calculate the forces and dimensionless numbers; (5) Calculate the sealing safety factor; (6) Change the initialization conditions and repeat steps (4) to (5) to obtain three sets of dimensionless numbers and sealing safety factors; (7) Establish a linear regression model between the sequestration safety factor and the dimensionless number, and calculate the CO2 sequestration safety factor of the saline aquifer; In step (4), the steps of studying the dynamics of CO2 leakage and calculating the forces and dimensionless numbers are as follows: by using numerical simulation, the dynamics of CO2 plume passing through cracks and leaking upwards are studied, and the forces and dimensionless numbers in the CO2 leakage process are calculated based on the data of the simulated monitoring points. The data of the simulated monitoring points include the density of CO2 and salt water, CO2 viscosity, CO2 flow velocity, capillary pressure and permeability at the monitoring points. In step (5), the sequestration safety factor is the ratio of CO2 leaked from the target reservoir to the overlying saline aquifer in the form of dissolution sequestration and mineralization sequestration, as shown in equation (1): In the formula, y is the sealing safety factor; M l M represents the total mass of CO2 leaking from the target reservoir to the overlying saline aquifer. d M represents the mass of the leaked CO2 that is dissolved and sequestered. m The mass of CO2 that is contained in the leaked CO2 in the form of mineralized sequestration; In step (6), the dimensionless number is: In the formula, Ca is the capillary number; Gr is the gravitational number; Bo is the bond number; F B F V F CP These are buoyancy, viscous force, and capillary force, respectively; ρ w ρ g The densities of saltwater and CO2 are respectively; g is the acceleration due to gravity; l cr The characteristic length; CO2 viscosity; CO2 flow rate; k is permeability; p cr Capillary pressure; In step (7), the linear regression model of the sealing safety factor and the dimensionless number is as follows: y = A·Ca + B·Bo + C(5) Where A, B, and C are all regression coefficients.

2. A safety assessment system for CO2 sequestration in saline aquifers based on dimensionless numbers Ca, Gr, and Bo, wherein the method described in claim 1 is characterized in that, include: The geological model building module is used to create a geological model containing fractures and initialize the model parameters; The calculation module is used to calculate dimensionless numbers and the sealing safety factor; The module for constructing a regression model for the storage safety factor is used to establish a linear regression model between the storage safety factor and the dimensionless number. The process of establishing a fractured geological model and initializing model parameters includes: establishing a saline aquifer CO2 sequestration geological model, including the target reservoir, interlayers, and overlying saline aquifers; establishing fracture models and monitoring points in the interlayers; and initializing the geological model. The calculation of dimensionless numbers and storage safety factors includes: studying the dynamics of CO2 leakage, calculating the forces and dimensionless numbers; calculating the storage safety factor; changing the initialization conditions to obtain three sets of dimensionless numbers and storage safety factors; The establishment of the linear regression model between the sequestration safety factor and the dimensionless number includes: establishing the linear regression model between the sequestration safety factor and the dimensionless number, and calculating the sequestration safety factor of CO2 in the saline aquifer.

3. A computer-readable storage medium having a program stored thereon, characterized in that, When executed by the processor, the program implements the steps in the method for assessing the safety of CO2 sequestration in saline aquifers based on dimensionless numbers Ca, Gr, and Bo as described in claim 1.

4. An electronic device comprising a memory, a processor, and a program stored in the memory and executable on the processor, characterized in that, When the processor executes the program, it implements the steps in the method for assessing the safety of CO2 sequestration in saline aquifers based on dimensionless numbers Ca, Gr, and Bo as described in claim 1.

Images