Method, device and equipment for determining safety of carbon dioxide burying and storage medium

Abstract

The invention discloses a carbon dioxide burying safety determination method, device and equipment and a storage medium, and belongs to the technical field of oil and gas field development. For carbon dioxide burying in a target area, on the basis of establishing a geologic model and a numerical simulation model, sensitivity analysis is carried out on parameters influencing the carbon dioxide leakage rate, that is, the influence degree of each parameter on the carbon dioxide leakage rate is determined; carbon dioxide leakage risk evaluation indexes at different spatial positions of the target area are obtained, and the carbon dioxide leakage risk is determined according to the carbon dioxide leakage risk so as to guide on-site carbon dioxide leakage monitoring, for example, a carbon dioxide leakage monitoring well is deployed at a spatial position with high carbon dioxide leakage risk; effective monitoring of the carbon dioxide leakage condition is achieved, it is ensured that remedial measures are taken in time when the leakage condition occurs, and therefore the safety of carbon dioxide burying is guaranteed.

Description

Technical Field

[0001] The present application relates to the technical field of oil and gas field development, and in particular to a method, device, equipment and storage medium for determining the safety of carbon dioxide storage. Background Art

[0002] Carbon dioxide storage is a technology that buries carbon dioxide captured in industrial or energy production processes underground to reduce the concentration of carbon dioxide in the atmosphere. Currently, carbon dioxide storage is considered one of the important means to reduce greenhouse gas emissions.

[0003] The storage sites of carbon dioxide usually include oil, gas and water reservoirs. When the concentration of carbon dioxide in the air exceeds 2%, it can cause damage to the human respiratory organs. More seriously, high concentrations of carbon dioxide can even cause death due to lack of oxygen and suffocation. Therefore, before injecting carbon dioxide into oil, gas and water reservoirs for storage, it is necessary to determine the safety of carbon dioxide storage, that is, to evaluate the safety of carbon dioxide storage, so as to avoid large-scale leakage after carbon dioxide storage, which will cause casualties on the ground and bring economic losses. Summary of the invention

[0004] The embodiment of the present application provides a method, device, equipment and storage medium for determining the safety of carbon dioxide storage. The technical solution is as follows:

[0005] In one aspect, a method for determining safety of carbon dioxide sequestration is provided, the method comprising:

[0006] Establishing a three-dimensional geological model of the target area and performing grid division; wherein the three-dimensional geological model includes a cap layer or multiple interlayers;

[0007] According to the type of the target area and the carbon dioxide storage mechanism, a numerical simulation model based on the three-dimensional geological model is established;

[0008] Based on the numerical simulation model, a process of injecting carbon dioxide into the target area for storage is numerically simulated to obtain a parameter set that affects the carbon dioxide leakage of each grid; wherein the parameter set includes at least: the porosity, permeability and water saturation of the reservoir in the target area, and the carbon dioxide injection pressure and carbon dioxide injection rate;

[0009] For any grid, obtain the weight value of each parameter in the parameter set of the grid that affects the carbon dioxide leakage; wherein the weight value is used to characterize the degree of influence of the corresponding parameter on the carbon dioxide leakage;

[0010] Determining a carbon sequestration risk coefficient of the grid based on a value of each parameter in a parameter set of the grid and a weight value corresponding to each parameter;

[0011] When the carbon sequestration risk coefficient of the grid is higher than a preset threshold, it is determined that the grid has a carbon dioxide leakage risk.

[0012] In a possible implementation, obtaining the weight value of each parameter in the parameter set of the grid affecting the carbon dioxide leakage includes:

[0013] Performing a first normalization process on a first type of parameters in a parameter set of the grid;

[0014] performing a second normalization process on the second type of parameters in the parameter set of the grid;

[0015] Based on the normalized data, a weight value of each parameter in the parameter set of the grid affecting the carbon dioxide leakage is obtained;

[0016] The first type of parameters include porosity, water saturation, carbon dioxide injection pressure and carbon dioxide injection rate, and the second type of parameters include permeability.

[0017] In a possible implementation manner, performing a first normalization process on the first type of parameters in the parameter set of the grid includes:

[0018] For any parameter in the first category of parameters, obtain a value set of the parameter, a maximum value and a minimum value in the value set;

[0019] Based on the value set of the parameter, the maximum value and the minimum value in the value set, a normalized value set of the parameter is obtained.

[0020] In a possible implementation manner, performing a second normalization process on the second type of parameters in the parameter set of the grid includes:

[0021] For any parameter in the second category of parameters, obtaining a logarithmic value set of the parameter, and a maximum value and a minimum value in the logarithmic value set;

[0022] Based on the logarithmic value set of the parameter, the maximum value and the minimum value in the logarithmic value set, a normalized value set of the parameter is obtained.

[0023] In a possible implementation, obtaining the weight value of each parameter in the parameter set of the grid affecting the carbon dioxide leakage based on the normalized parameter value includes:

[0024] For any parameter, obtaining a value set of carbon dioxide leakage when the parameter is a first value;

[0025] Determine a maximum carbon dioxide leakage value and an average carbon dioxide leakage value in the value set;

[0026] Obtaining a derivative value set of carbon dioxide leakage when the parameter value is a normalized value set;

[0027] Based on the average carbon dioxide leakage amount, the maximum carbon dioxide leakage amount and the derivative value set, a weight value of the parameter affecting the carbon dioxide leakage amount is obtained.

[0028] In a possible implementation, establishing a numerical simulation model based on the three-dimensional geological model according to the type of the target area and the carbon dioxide storage mechanism includes:

[0029] When the target area is an oil reservoir, a numerical simulation model based on the three-dimensional geological model is established according to the structural storage mechanism of carbon dioxide, the dissolution storage mechanism of carbon dioxide, the bound storage mechanism of carbon dioxide and the mineralization storage mechanism of carbon dioxide; wherein the numerical simulation model includes a component model for obtaining crude oil components.

[0030] In a possible implementation, the method further includes:

[0031] After determining that the grid has a risk of carbon dioxide leakage, outputting a safety prompt;

[0032] The safety prompt is used to instruct to monitor the carbon dioxide leakage in the grid.

[0033] In another aspect, a device for determining safety of carbon dioxide sequestration is provided, the device comprising:

[0034] A first building module is configured to build a three-dimensional geological model of the target area and perform grid division; wherein the three-dimensional geological model includes a cap layer or multiple interlayers;

[0035] A second establishment module is configured to establish a numerical simulation model based on the three-dimensional geological model according to the type of the target area and the carbon dioxide storage mechanism;

[0036] A numerical simulation module is configured to perform numerical simulation on the process of injecting carbon dioxide into the target area for storage based on the numerical simulation model, and obtain a parameter set that affects the carbon dioxide leakage of each grid; wherein the parameter set includes at least: the porosity, permeability and water saturation of the reservoir in the target area, and the carbon dioxide injection pressure and carbon dioxide injection rate;

[0037] An acquisition module is configured to acquire, for any grid, a weight value of each parameter in a parameter set of the grid affecting the carbon dioxide leakage; wherein the weight value is used to characterize the degree of influence of the corresponding parameter on the carbon dioxide leakage;

[0038] A determination module is configured to determine the carbon storage risk coefficient of the grid based on the value of each parameter in the parameter set of the grid and the weight value corresponding to each parameter; when the carbon storage risk coefficient of the grid is higher than a preset threshold, it is determined that the grid has a risk of carbon dioxide leakage.

[0039] In a possible implementation, the acquisition module is configured to:

[0040] Performing a first normalization process on a first type of parameters in a parameter set of the grid;

[0041] performing a second normalization process on the second type of parameters in the parameter set of the grid;

[0042] Based on the normalized data, a weight value of each parameter in the parameter set of the grid affecting the carbon dioxide leakage is obtained;

[0043] The first type of parameters include porosity, water saturation, carbon dioxide injection pressure and carbon dioxide injection rate, and the second type of parameters include permeability.

[0044] In a possible implementation, the acquisition module is configured to:

[0045] For any parameter in the first category of parameters, obtain a value set of the parameter, a maximum value and a minimum value in the value set;

[0046] Based on the value set of the parameter, the maximum value and the minimum value in the value set, a normalized value set of the parameter is obtained.

[0047] In a possible implementation, the acquisition module is configured to:

[0048] For any parameter in the second category of parameters, obtaining a logarithmic value set of the parameter, and a maximum value and a minimum value in the logarithmic value set;

[0049] Based on the logarithmic value set of the parameter, the maximum value and the minimum value in the logarithmic value set, a normalized value set of the parameter is obtained.

[0050] In a possible implementation, the acquisition module is configured to:

[0051] For any parameter, obtaining a value set of carbon dioxide leakage when the parameter is a first value;

[0052] Determine a maximum carbon dioxide leakage value and an average carbon dioxide leakage value in the value set;

[0053] Obtaining a derivative value set of carbon dioxide leakage when the parameter value is a normalized value set;

[0054] Based on the average carbon dioxide leakage amount, the maximum carbon dioxide leakage amount and the derivative value set, a weight value of the parameter affecting the carbon dioxide leakage amount is obtained.

[0055] In a possible implementation, the second model building module is configured to:

[0056] When the target area is an oil reservoir, a numerical simulation model based on the three-dimensional geological model is established according to the structural storage mechanism of carbon dioxide, the dissolution storage mechanism of carbon dioxide, the bound storage mechanism of carbon dioxide and the mineralization storage mechanism of carbon dioxide; wherein the numerical simulation model includes a component model for obtaining crude oil components.

[0057] In a possible implementation manner, the device further includes:

[0058] The output module is configured to output a safety prompt after determining that there is a risk of carbon dioxide leakage in the grid; wherein the safety prompt is used to instruct to monitor the carbon dioxide leakage in the grid.

[0059] On the other hand, a computer device is provided, which includes a processor and a memory, wherein the memory stores at least one program code, and the at least one program code is loaded and executed by the processor to implement the above-mentioned method for determining the safety of carbon dioxide sequestration.

[0060] On the other hand, a computer-readable storage medium is provided, wherein at least one program code is stored in the storage medium, and the at least one program code is loaded and executed by a processor to implement the above-mentioned method for determining the safety of carbon dioxide sequestration.

[0061] On the other hand, a computer program product or a computer program is provided, which includes a computer program code, which is stored in a computer-readable storage medium. A processor of a computer device reads the computer program code from the computer-readable storage medium, and the processor executes the computer program code, so that the computer device executes the above-mentioned method for determining the safety of carbon dioxide sequestration.

[0062] For the carbon dioxide storage in the target area, the embodiment of the present application will, on the basis of establishing a geological model and a numerical simulation model, also conduct a sensitivity analysis on the parameters that affect the carbon dioxide leakage. That is, the embodiment of the present application will determine the degree of influence of various parameters on the carbon dioxide leakage, and then obtain the carbon dioxide leakage risk evaluation index (also called the carbon storage risk coefficient) at different spatial locations in the target area, and determine the carbon dioxide leakage risk accordingly to guide on-site carbon dioxide leakage monitoring. For example, carbon dioxide leakage monitoring wells may be deployed at spatial locations with a high risk of carbon dioxide leakage to achieve effective monitoring of carbon dioxide leakage and ensure that remedial measures are taken in a timely manner when leakage occurs, thereby ensuring the safety of carbon dioxide storage. BRIEF DESCRIPTION OF THE DRAWINGS

[0063] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0064] Figure 1 It is a flow chart of a method for determining the safety of carbon dioxide sequestration provided in an embodiment of the present application;

[0065] Figure 2 is a schematic diagram of a geological model of a target area provided in an embodiment of the present application;

[0066] Figure 3 This is an analysis result of the effect of porosity on carbon dioxide leakage provided in an embodiment of the present application;

[0067] Figure 4 This is another analysis result of the effect of porosity on carbon dioxide leakage provided in an embodiment of the present application;

[0068] Figure 5 This is an analysis result of the effect of permeability on carbon dioxide leakage provided by an embodiment of the present application;

[0069] Figure 6 This is another analysis result of the effect of permeability on carbon dioxide leakage provided by an embodiment of the present application;

[0070] Figure 7 This is an analysis result of the effect of water saturation on carbon dioxide leakage provided in an embodiment of the present application;

[0071] Figure 8 This is another analysis result of the effect of water saturation on carbon dioxide leakage provided in an embodiment of the present application;

[0072] Fig. 9This is an analysis result of the effect of carbon dioxide injection pressure on carbon dioxide leakage provided in an embodiment of the present application;

[0073] Fig.10 This is another analysis result of the effect of carbon dioxide injection pressure on carbon dioxide leakage provided in an embodiment of the present application;

[0074] Fig.11 This is an analysis result of the effect of carbon dioxide injection rate on carbon dioxide leakage provided in an embodiment of the present application;

[0075] Fig.12 This is another analysis result of the effect of carbon dioxide injection rate on carbon dioxide leakage provided in an embodiment of the present application;

[0076] Fig.13 is a schematic diagram of the distribution of carbon dioxide leakage risks at different spatial locations in a target area provided by an embodiment of the present application;

[0077] Fig.14 It is a schematic diagram of the deployment location of a suggested carbon dioxide leakage monitoring well provided in an embodiment of the present application;

[0078] Fig.15 It is a structural schematic diagram of a safety determination device for carbon dioxide storage provided in an embodiment of the present application;

[0079] Fig.16 It is a structural diagram of a computer device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0080] In order to make the objectives, technical solutions and advantages of the present application clearer, the implementation methods of the present application will be further described in detail below with reference to the accompanying drawings.

[0081] In this application, the terms "first", "second", etc. are used to distinguish the same or similar items with substantially the same role and function. It should be understood that there is no logical or temporal dependency between "first", "second", and "nth", nor is there a limitation on the quantity and execution order. It should also be understood that although the following description uses the terms first, second, etc. to describe various elements, these elements should not be limited by the terms.

[0082] These terms are only used to distinguish one element from another element. For example, without departing from the scope of various examples, a first element can be referred to as a second element, and similarly, a second element can also be referred to as a first element. Both the first element and the second element can be elements, and in some cases, can be separate and different elements.

[0083] Here, at least one means one or more than one, for example, at least one element can be one element, two elements, three elements, or any other integer greater than or equal to one. And multiple means two or more than two, for example, multiple elements can be two elements, three elements, or any other integer greater than or equal to two.

[0084] The "and / or" mentioned in this article indicates that there can be three relationships. For example, A and / or B can mean: A exists alone, A and B exist at the same time, and B exists alone. The character " / " generally indicates that the objects before and after are in an "or" relationship.

[0085] It should be noted that the information (including but not limited to user device information, user personal information, etc.), data (including but not limited to data used for analysis, stored data, displayed data, etc.) and signals involved in this application are all authorized by the user or fully authorized by all parties, and the collection, use and processing of relevant data must comply with relevant laws, regulations and standards in the relevant regions.

[0086] The embodiment of the present application provides a safety determination scheme for carbon dioxide storage, which provides a method for evaluating the safety of carbon dioxide storage during carbon dioxide injection. The execution subject of the scheme is a computer device. Exemplarily, the types of the above-mentioned computer devices include but are not limited to: desktop computers, laptops, smart phones, and tablet computers, etc., which are not limited in the present application.

[0087] Among them, the embodiment of the present application evaluates the carbon storage safety of the target area (such as oil, gas and water reservoirs) where carbon dioxide injection is implemented, and can solve the following two problems: on the one hand, it solves the problem that the related technology uses instruments and equipment to monitor carbon dioxide leakage, which is not only costly but also unable to achieve long-term monitoring; on the other hand, it at least solves the problem that the related technology cannot quantitatively characterize the safety of carbon dioxide storage.

[0088] For example, the carbon storage safety assessment method provided in the embodiment of the present application is mainly applied to oil, gas and water reservoirs for carbon dioxide injection and storage. By calculating the carbon dioxide leakage risk assessment index at different spatial positions in the target area, a quantitative assessment of the carbon dioxide leakage risk is achieved, thereby guiding the monitoring of carbon dioxide leakage in the target area during the on-site carbon dioxide injection process, for example, by deploying carbon dioxide leakage monitoring wells to monitor carbon dioxide leakage, thereby ensuring the safety of carbon dioxide injection.

[0089] The carbon sequestration safety evaluation method provided in the embodiment of the present application is described in detail through the following implementation method.

[0090] Figure 11 is a flow chart of a method for determining the safety of carbon dioxide storage provided by an embodiment of the present application. The execution subject of the method is a computer device. Figure 1 The method flow includes:

[0091] 101. Computer equipment establishes a three-dimensional geological model of the target area and performs grid division; wherein the three-dimensional geological model includes a cap layer or multiple interlayers.

[0092] In the embodiment of the present application, the computer equipment uses geological modeling software to establish a geological model of the target area. For example, the target area for carbon dioxide injection is X reservoir, see Figure 2 The number of grids in the geological model of the X reservoir is 84*38*16, and the average grid size is 65*65*10 meters; the longest part of the geological model is 5815 meters, and the widest part is 2550 meters; in addition, the geological model contains two cap rock areas and two reservoir areas.

[0093] In addition, the embodiments of the present application are Figure 2 The numerical axis shown in the lower right corner represents the pressure value of each grid in the geological model, in kPa.

[0094] 102. The computer equipment establishes a numerical simulation model based on the three-dimensional geological model according to the type of target area and the carbon dioxide storage mechanism.

[0095] Exemplarily, the types of the above-mentioned target areas include but are not limited to oil, gas and water reservoirs, which are not limited in this application. In the embodiment of this application, the computer equipment uses oil reservoir numerical simulation software to simulate the process of injecting carbon dioxide into the target area. In other words, the numerical simulation model lays the foundation for accurately simulating and predicting carbon dioxide leakage and evaluating the safety of carbon storage.

[0096] The first point that needs to be explained is that in order to accurately simulate the carbon dioxide injection process in the target area, the carbon dioxide tectonic storage mechanism, carbon dioxide dissolution storage mechanism (including the dissolution of carbon dioxide in oil and water), carbon dioxide bound storage and carbon dioxide mineralization storage mechanism are also comprehensively considered in the process of establishing the numerical simulation model.

[0097] The second point to be explained is that since different components of crude oil will react with carbon dioxide, it is also necessary to describe this interaction. For this reason, if the target area is an oil reservoir, a component model for obtaining crude oil components must be established. In other words, if the target area is an oil reservoir, the numerical simulation model includes a component model for obtaining crude oil components.

[0098] In summary, according to the type of target area and the carbon dioxide storage mechanism, a numerical simulation model based on a three-dimensional geological model is established, including but not limited to the following methods:

[0099] When the target area is an oil reservoir, a numerical simulation model based on a three-dimensional geological model is established according to the structural storage mechanism of carbon dioxide, the dissolution storage mechanism of carbon dioxide, the bound storage mechanism of carbon dioxide and the mineralization storage mechanism of carbon dioxide; wherein the numerical simulation model includes a component model for obtaining crude oil components.

[0100] The following is an introduction to the structural storage mechanism of carbon dioxide, the dissolution storage mechanism of carbon dioxide, the bound storage mechanism of carbon dioxide and the mineralization storage mechanism of carbon dioxide.

[0101] Regarding the structural storage of carbon dioxide, it is a storage mechanism that prevents the further migration of carbon dioxide based on the structural closure characteristics of the geological model.

[0102] Regarding the dissolution of carbon dioxide in water, the fugacity of carbon dioxide in water is calculated based on Henry's law, and the amount of carbon dioxide dissolved in water is further calculated based on phase equilibrium.

[0103] In view of the dissolution of carbon dioxide in oil, a fluid pseudo-component model was established based on the PR equation of state, and the proportion of carbon dioxide in crude oil components was further calculated.

[0104] For the bound storage of carbon dioxide, based on the Carlson two-phase hysteresis model theory, the storage amount of carbon dioxide bound in the pores of rock particles due to changes in saturation during carbon dioxide migration is simulated.

[0105] Regarding the mineralization and storage of carbon dioxide, the mineralization of carbon dioxide includes two steps, carbonization (combination with rocks) and hydrolysis (reaction with water). Among them, the ion reaction of carbon dioxide in water is based on the transient reaction model, and the dissolution and precipitation reaction of minerals is based on the transition state theory.

[0106] 103. Based on the established numerical simulation model, the computer equipment numerically simulates the process of injecting carbon dioxide into the target area for storage, and obtains a set of parameters that affect the carbon dioxide leakage of each grid; wherein the parameter set includes at least: the porosity, permeability and water saturation of the reservoir in the target area, as well as the carbon dioxide injection pressure and carbon dioxide injection rate.

[0107] In the embodiment of the present application, in order to obtain the degree of influence of different parameters on the carbon dioxide leakage, the parameters for parameter sensitivity analysis (step 104 below) include at least: the porosity, permeability, water saturation of the reservoir, and the injection pressure and injection rate during the carbon dioxide injection process.

[0108] 104. For any grid, the computer device obtains the weight value of each parameter in the parameter set of the grid that affects the carbon dioxide leakage; wherein, for any parameter, the weight value of the parameter is used to characterize the degree of influence of the parameter on the carbon dioxide leakage.

[0109] In the process of parameter sensitivity analysis, in order to avoid the influence of the difference in the values ​​of different parameters on the carbon storage safety evaluation results, the embodiment of the present application will normalize the parameters in the above parameter set. In the embodiment of the present application, the porosity, water saturation, carbon dioxide injection pressure and carbon dioxide injection rate are subjected to the first normalization, and the permeability is subjected to the second normalization. In other words, the weight value of each parameter in the parameter set of the grid affecting the carbon dioxide leakage is obtained, which is achieved in the following way:

[0110] A first normalization process is performed on the first type of parameters in the parameter set of the grid; a second normalization process is performed on the second type of parameters in the parameter set of the grid; based on the normalized data, a weight value of each parameter in the parameter set of the grid affecting the carbon dioxide leakage is obtained; wherein the first type of parameters includes porosity, water saturation, carbon dioxide injection pressure and carbon dioxide injection rate, and the second type of parameters includes permeability.

[0111] Taking the above-mentioned first normalization process as maximum and minimum normalization as an example, the above-mentioned first normalization process is performed on the first type of parameters in the parameter set of the grid, including:

[0112] For any parameter in the first category of parameters, obtain the value set of the parameter, the maximum value and the minimum value in the value set of the parameter; based on the value set of the parameter, the maximum value and the minimum value in the value set of the parameter, obtain the normalized value set of the parameter.

[0113] The above process can be described mathematically as the following formula (1).

[0114]

[0115] Wherein, a is a parameter set that affects the amount of carbon dioxide leakage, and the parameter set includes at least the porosity, permeability, water saturation, carbon dioxide injection pressure and carbon dioxide injection rate of the reservoir in the target area; x a is the value set of parameter a; x norm，a is the normalized value set of parameter a; min(x a ) is the minimum value in the set of values ​​of parameter a; max(x a ) is the maximum value in the set of values ​​of parameter a.

[0116] Taking the second normalization process as logarithmic normalization as an example, the second normalization process is performed on the second type of parameters in the parameter set of the grid of the parameter, including:

[0117] For any parameter in the second category, obtain the logarithmic value set of the parameter, the maximum value and the minimum value in the logarithmic value set of the parameter; based on the logarithmic value set of the parameter, the maximum value and the minimum value in the logarithmic value set of the parameter, obtain the normalized value set of the parameter.

[0118] The above process can be described mathematically as the following formula (2).

[0119]

[0120] Among them, lg(x a ) is the logarithmic value set of parameter a; min(lg(x a )) is the minimum value in the logarithmic value set of parameter a; max(lg(x a )) is the maximum value in the set of logarithmic values ​​of parameter a.

[0121] In addition, the embodiments of the present application are Figures 3 to 12 The analytical results of sensitivity analysis on porosity, permeability, water saturation, CO2 injection pressure and CO2 injection rate are shown respectively.

[0122] In a possible implementation, the present application embodiment selects five parameters of reservoir porosity, permeability, water saturation, carbon dioxide injection pressure and carbon dioxide injection rate to perform sensitivity analysis on the carbon dioxide leakage. In other words, the weight value of each parameter in the parameter set of the grid affecting the carbon dioxide leakage is obtained, which is implemented in the following way:

[0123] For any parameter, obtain a value set of carbon dioxide leakage when the parameter is a first value; determine the maximum carbon dioxide leakage and the average carbon dioxide leakage in the value set; obtain a derivative value set of carbon dioxide leakage when the parameter value is the above normalized value set; based on the above average carbon dioxide leakage, the above maximum carbon dioxide leakage and the above derivative value set of carbon dioxide leakage, obtain the weight value of the parameter affecting the carbon dioxide leakage.

[0124] The above process can be described mathematically as the following formula (3).

[0125]

[0126] Where w(a) is the weight value set of parameter a affecting the carbon dioxide leakage; L(x a ) is parameter a with value xa The value set of carbon dioxide leakage at the time; max(L(x a )) is the parameter a with the value x a The maximum value of carbon dioxide leakage in the value set of carbon dioxide leakage at that time; The value of parameter a is x a The average value of carbon dioxide leakage in the set of carbon dioxide leakage values; L′ norm (x norm,a ) is parameter a with value x norm，a The derivative value set of carbon dioxide leakage at time .

[0127] Exemplarily, the weight values ​​of the five parameters of porosity, permeability, water saturation, carbon dioxide injection pressure and carbon dioxide injection rate affecting carbon dioxide leakage are calculated by the above formulas (1) to (3) to be 0.0947, 0.2212, 0.0092, 1.3316 and 0.4355 respectively.

[0128] 105. The computer device determines the carbon sequestration risk coefficient of the grid based on the value of each parameter in the parameter set of the grid and the weight value corresponding to each parameter; when the carbon sequestration risk coefficient of the grid is higher than a preset threshold, it is determined that the grid has a risk of carbon dioxide leakage.

[0129] In the embodiment of the present application, based on the value of each parameter in the parameter set at different spatial locations of the target area and the weight value corresponding to each parameter, a weighted method is used to calculate the carbon dioxide leakage risk assessment index at different spatial locations, that is, the carbon storage risk coefficient. The larger the value, the higher the carbon dioxide leakage risk. The calculation formula of the carbon dioxide leakage risk assessment index is as follows:

[0130]

[0131] Among them, I is the carbon dioxide leakage risk assessment index, that is, the carbon storage risk coefficient; n is the number of parameters that affect the carbon dioxide leakage.

[0132] For example, the calculation result of the carbon dioxide leakage risk assessment index is as follows: Fig.13 That is, Fig.13 The distribution of carbon dioxide leakage risk at different spatial locations in the target area is shown. Fig.13 The numerical axis shown in the lower right corner represents the value of the carbon dioxide leakage risk assessment index of each grid in the geological model.

[0133] In one possible implementation, after determining that a certain spatial location within the target area has a risk of carbon dioxide leakage, the computer device will output a safety prompt; wherein the safety prompt is used to instruct the monitoring of carbon dioxide leakage in the grid. That is, the computer device will guide the on-site monitoring of carbon dioxide leakage based on the risk of carbon dioxide leakage. For example, carbon dioxide leakage monitoring wells are deployed at spatial locations with a high risk of carbon dioxide leakage to achieve effective monitoring of carbon dioxide leakage and ensure that remedial measures are taken in a timely manner when leakage occurs, thereby ensuring the safety of carbon dioxide storage. Among them, the recommended deployment locations of carbon dioxide leakage monitoring wells are as follows: Fig.14 shown.

[0134] The embodiment of the present application provides the industry with a low-cost quantitative evaluation scheme for the safety of carbon dioxide storage that takes into account the mechanism of carbon dioxide storage. In detail, for the carbon dioxide storage in the target area, the scheme will conduct a sensitivity analysis on the parameters that affect the carbon dioxide leakage on the basis of establishing a geological model and a numerical simulation model, that is, the scheme will determine the degree of influence of each parameter on the carbon dioxide leakage, and then obtain the carbon dioxide leakage risk evaluation index at different spatial locations in the target area, and determine the carbon dioxide leakage risk accordingly, so as to guide the on-site carbon dioxide leakage monitoring, such as deploying carbon dioxide leakage monitoring wells at spatial locations with high risks of carbon dioxide leakage, to achieve effective monitoring of carbon dioxide leakage, and ensure that timely remedial measures are taken when leakage occurs, thereby ensuring the safety of carbon dioxide storage.

[0135] Fig.15 Schematic diagram of a safety determination device for carbon dioxide storage provided in an embodiment of the present application. Fig.15 , the device comprises:

[0136] The first building module 1501 is configured to build a three-dimensional geological model of the target area and perform grid division; wherein the three-dimensional geological model includes a cap layer or multiple interlayers;

[0137] The second establishing module 1502 is configured to establish a numerical simulation model based on the three-dimensional geological model according to the type of the target area and the carbon dioxide sequestration mechanism;

[0138] The numerical simulation module 1503 is configured to perform numerical simulation on the process of injecting carbon dioxide into the target area for storage based on the numerical simulation model, and obtain a parameter set that affects the carbon dioxide leakage of each grid; wherein the parameter set at least includes: the porosity, permeability and water saturation of the reservoir in the target area, and the carbon dioxide injection pressure and carbon dioxide injection rate;

[0139] The acquisition module 1504 is configured to acquire, for any grid, a weight value of each parameter in the parameter set of the grid affecting the carbon dioxide leakage; wherein the weight value is used to represent the degree of influence of the corresponding parameter on the carbon dioxide leakage;

[0140] Determination module 1505 is configured to determine the carbon storage risk coefficient of the grid based on the value of each parameter in the parameter set of the grid and the weight value corresponding to each parameter; when the carbon storage risk coefficient of the grid is higher than a preset threshold, it is determined that the grid has a risk of carbon dioxide leakage.

[0141] The embodiment of the present application provides the industry with a low-cost quantitative evaluation scheme for the safety of carbon dioxide storage that takes into account the mechanism of carbon dioxide storage. In detail, for the carbon dioxide storage in the target area, the scheme will conduct a sensitivity analysis on the parameters that affect the carbon dioxide leakage on the basis of establishing a geological model and a numerical simulation model, that is, the scheme will determine the degree of influence of each parameter on the carbon dioxide leakage, and then obtain the carbon dioxide leakage risk evaluation index at different spatial locations in the target area, and determine the carbon dioxide leakage risk accordingly, so as to guide the on-site carbon dioxide leakage monitoring, such as deploying carbon dioxide leakage monitoring wells at spatial locations with high risks of carbon dioxide leakage, to achieve effective monitoring of carbon dioxide leakage, and ensure that timely remedial measures are taken when leakage occurs, thereby ensuring the safety of carbon dioxide storage.

[0142] In a possible implementation, the acquisition module is configured to:

[0143] Performing a first normalization process on a first type of parameters in a parameter set of the grid;

[0144] performing a second normalization process on the second type of parameters in the parameter set of the grid;

[0145] Based on the normalized data, a weight value of each parameter in the parameter set of the grid affecting the carbon dioxide leakage is obtained;

[0146] The first type of parameters include porosity, water saturation, carbon dioxide injection pressure and carbon dioxide injection rate, and the second type of parameters include permeability.

[0147] In a possible implementation, the acquisition module is configured to:

[0148] For any parameter in the first category of parameters, obtain a value set of the parameter, a maximum value and a minimum value in the value set;

[0149] Based on the value set of the parameter, the maximum value and the minimum value in the value set, a normalized value set of the parameter is obtained.

[0150] In a possible implementation, the acquisition module is configured to:

[0151] For any parameter in the second category of parameters, obtaining a logarithmic value set of the parameter, and a maximum value and a minimum value in the logarithmic value set;

[0152] Based on the logarithmic value set of the parameter, the maximum value and the minimum value in the logarithmic value set, a normalized value set of the parameter is obtained.

[0153] In a possible implementation, the acquisition module is configured to:

[0154] For any parameter, obtaining a value set of carbon dioxide leakage when the parameter is a first value;

[0155] Determine a maximum carbon dioxide leakage value and an average carbon dioxide leakage value in the value set;

[0156] Obtaining a derivative value set of carbon dioxide leakage when the parameter value is a normalized value set;

[0157] Based on the average carbon dioxide leakage amount, the maximum carbon dioxide leakage amount and the derivative value set, a weight value of the parameter affecting the carbon dioxide leakage amount is obtained.

[0158] In a possible implementation, the second model building module is configured to:

[0159] When the target area is an oil reservoir, a numerical simulation model based on the three-dimensional geological model is established according to the structural storage mechanism of carbon dioxide, the dissolution storage mechanism of carbon dioxide, the bound storage mechanism of carbon dioxide and the mineralization storage mechanism of carbon dioxide; wherein the numerical simulation model includes a component model for obtaining crude oil components.

[0160] In a possible implementation, the device further includes:

[0161] The output module is configured to output a safety prompt after determining that there is a risk of carbon dioxide leakage in the grid; wherein the safety prompt is used to instruct to monitor the carbon dioxide leakage in the grid.

[0162] All the above optional technical solutions can be arbitrarily combined to form optional embodiments of the present application, which will not be described one by one here.

[0163] It should be noted that: when the safety determination device for carbon dioxide storage provided in the above embodiment determines the safety of carbon dioxide storage during the injection of carbon dioxide, 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 device is divided into different functional modules to complete all or part of the functions described above. In addition, the safety determination device for carbon dioxide storage provided in the above embodiment and the safety determination method embodiment of carbon dioxide storage belong to the same concept. The specific implementation process is detailed in the method embodiment and will not be repeated here.

[0164] Fig.16 It is a structural diagram of a computer device 1600 provided in an embodiment of the present application.

[0165] Typically, the computer device 1600 includes a processor 1601 and a memory 1602 .

[0166] The processor 1601 includes one or more processing cores, such as a 4-core processor, an 8-core processor, etc. The processor 1601 is implemented in at least one hardware form of DSP (Digital Signal Processing), FPGA (Field-Programmable Gate Array), and PLA (Programmable Logic Array). Alternatively, the processor 1601 includes a main processor and a coprocessor. The main processor is a processor for processing data in the awake state, also known as a CPU (Central Processing Unit); the coprocessor is a low-power processor for processing data in the standby state. In one possible implementation, the processor 1601 is integrated with a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content to be displayed on the display screen. In one possible implementation, the processor 1601 also includes an AI (Artificial Intelligence) processor, which is used to process computing operations related to machine learning.

[0167] The memory 1602 includes one or more computer-readable storage media, which are non-transitory. The memory 1602 also includes a high-speed random access memory and a non-volatile memory, such as one or more disk storage devices and flash memory storage devices. In one possible implementation, the non-transitory computer-readable storage medium in the memory 1602 is used to store at least one program code, which is used to be executed by the processor 1601 to implement the safety determination method for carbon dioxide sequestration provided in the method embodiment of the present application.

[0168] In a possible implementation, the computer device 1600 further includes: a peripheral device interface 1603 and at least one peripheral device. The processor 1601, the memory 1602 and the peripheral device interface 1603 are connected via a bus or a signal line. Each peripheral device is connected to the peripheral device interface 1603 via a bus, a signal line or a circuit board. The peripheral device includes: at least one of a radio frequency circuit 1604, a display screen 1605, a camera assembly 1606, an audio circuit 1607, a positioning assembly 1608 and a power supply 1609.

[0169] The peripheral device interface 1603 is used to connect at least one peripheral device related to I / O (Input / Output) to the processor 1601 and the memory 1602. In one possible implementation, the processor 1601, the memory 1602, and the peripheral device interface 1603 are integrated on the same chip or circuit board; in another possible implementation, any one or two of the processor 1601, the memory 1602, and the peripheral device interface 1603 are implemented on a separate chip or circuit board, which is not limited in this application.

[0170] The radio frequency circuit 1604 is used to receive and transmit RF (Radio Frequency) signals, also known as electromagnetic signals. The radio frequency circuit 1604 communicates with the communication network and other communication devices through electromagnetic signals. The radio frequency circuit 1604 converts the electrical signal into an electromagnetic signal for transmission, or converts the received electromagnetic signal into an electrical signal. In one possible implementation, the radio frequency circuit 1604 includes: an antenna system, an RF transceiver, one or more amplifiers, a tuner, an oscillator, a digital signal processor, a codec chipset, a user identity module card, and the like. The radio frequency circuit 1604 communicates with other terminals through at least one wireless communication protocol. The wireless communication protocol includes, but is not limited to: the World Wide Web, a metropolitan area network, an intranet, various generations of mobile communication networks (2G, 3G, 4G and 5G), a wireless local area network and / or a WiFi (Wireless Fidelity) network. In one possible implementation, the radio frequency circuit 1604 also includes circuits related to NFC (Near Field Communication), which is not limited in this application.

[0171] The display screen 1605 is used to display a UI (User Interface). The UI includes graphics, text, icons, videos, and any combination thereof. In the case where the display screen 1605 is a touch display screen, the display screen 1605 also has the ability to collect touch signals on the surface or above the surface of the display screen 1605. The touch signal is input to the processor 1601 as a control signal for processing. At this time, the display screen 1605 is also used to provide virtual buttons and / or virtual keyboards, also known as soft buttons and / or soft keyboards. In one possible implementation, the display screen 1605 is one, which is arranged on the front panel of the computer device 1600; in another possible implementation, the display screen 1605 is at least two, which are arranged on different surfaces of the computer device 1600 or are folded; in another possible implementation, the display screen 1605 is a flexible display screen, which is arranged on a curved surface or a folded surface of the computer device 1600. Alternatively, the display screen 1605 is arranged in a non-rectangular irregular shape, that is, a special-shaped screen. The display screen 1605 is made of materials such as LCD (Liquid Crystal Display) and OLED (Organic Light-Emitting Diode).

[0172] The camera assembly 1606 is used to capture images or videos. In one possible implementation, the camera assembly 1606 includes a front camera and a rear camera. Typically, the front camera is disposed on the front panel of the terminal, and the rear camera is disposed on the back of the terminal. In one possible implementation, there are at least two rear cameras, which are any one of a main camera, a depth of field camera, a wide-angle camera, and a telephoto camera, so as to realize the fusion of the main camera and the depth of field camera to realize the background blur function, the fusion of the main camera and the wide-angle camera to realize panoramic shooting and VR (Virtual Reality) shooting function or other fusion shooting functions. In another possible implementation, the camera assembly 1606 also includes a flash. The flash is a monochrome temperature flash or a dual-color temperature flash. A dual-color temperature flash refers to a combination of a warm light flash and a cold light flash, which is used for light compensation at different color temperatures.

[0173] The audio circuit 1607 includes a microphone and a speaker. The microphone is used to collect sound waves from the user and the environment, and convert the sound waves into electrical signals and input them into the processor 1601 for processing, or input them into the radio frequency circuit 1604 to achieve voice communication. For the purpose of stereo acquisition or noise reduction, there are multiple microphones, which are respectively arranged at different parts of the computer device 1600. Alternatively, the microphone is an array microphone or an omnidirectional acquisition microphone. The speaker is used to convert the electrical signal from the processor 1601 or the radio frequency circuit 1604 into sound waves. The speaker is a traditional film speaker or a piezoelectric ceramic speaker. When the speaker is a piezoelectric ceramic speaker, not only can the electrical signal be converted into sound waves audible to humans, but also the electrical signal can be converted into sound waves inaudible to humans for purposes such as ranging. In a possible implementation, the audio circuit 1607 also includes a headphone jack.

[0174] The positioning component 1608 is used to locate the current geographic location of the computer device 1600 to implement navigation or LBS (Location Based Service). The positioning component 1608 can be a positioning component based on the US GPS (Global Positioning System), China's Beidou system, Russia's Grenas system, or the European Union's Galileo system.

[0175] The power supply 1609 is used to power various components in the computer device 1600. The power supply 1609 is an alternating current, a direct current, a disposable battery, or a rechargeable battery. In the case where the power supply 1609 includes a rechargeable battery, the rechargeable battery is a wired rechargeable battery or a wireless rechargeable battery. A wired rechargeable battery is a battery charged through a wired line, and a wireless rechargeable battery is a battery charged through a wireless coil. The rechargeable battery is also used to support fast charging technology.

[0176] Those skilled in the art will understand that Fig.16 The structure shown in the figure does not constitute a limitation on the computer device 1600, and the computer device 1600 may include more or less components than shown in the figure, or combine some components, or adopt a different arrangement of components.

[0177] In an exemplary embodiment, a computer-readable storage medium is also provided, such as a memory including a program code, and the program code can be executed by a processor in a computer device to complete the safety determination method of carbon dioxide sequestration in the above embodiment. For example, the computer-readable storage medium can be a read-only memory (ROM), a random access memory (RAM), a compact disc read-only memory (CD-ROM), a magnetic tape, a floppy disk, an optical data storage device, etc.

[0178] In an exemplary embodiment, a computer program product or a computer program is also provided, which includes a computer program code, and the computer program code is stored in a computer-readable storage medium. A processor of a computer device reads the computer program code from the computer-readable storage medium, and the processor executes the computer program code, so that the computer device executes the above-mentioned carbon dioxide storage safety determination method.

[0179] A person skilled in the art will understand that all or part of the steps to implement the above embodiments may be accomplished by hardware or by instructing related hardware through a program, and the program may be stored in a computer-readable storage medium, and the above-mentioned storage medium may be a read-only memory, a disk or an optical disk, etc.

[0180] The above description is only an optional embodiment of the present application and is not intended to limit the present application. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present application shall be included in the protection scope of the present application.

Claims

1. A method for determining the safety of carbon dioxide storage, characterized in that: The method comprises: Establishing a three-dimensional geological model of the target area and performing grid division; wherein the three-dimensional geological model includes a cap layer or multiple interlayers; According to the type of the target area and the carbon dioxide storage mechanism, a numerical simulation model based on the three-dimensional geological model is established; Based on the numerical simulation model, a process of injecting carbon dioxide into the target area for storage is numerically simulated to obtain a parameter set that affects the carbon dioxide leakage of each grid; wherein the parameter set includes at least: the porosity, permeability and water saturation of the reservoir in the target area, and the carbon dioxide injection pressure and carbon dioxide injection rate; For any grid, obtain the weight value of each parameter in the parameter set of the grid that affects the carbon dioxide leakage; wherein the weight value is used to characterize the degree of influence of the corresponding parameter on the carbon dioxide leakage; Determining a carbon sequestration risk coefficient of the grid based on a value of each parameter in a parameter set of the grid and a weight value corresponding to each parameter; When the carbon sequestration risk coefficient of the grid is higher than a preset threshold, it is determined that the grid has a carbon dioxide leakage risk.

2. The method according to claim 1, characterized in that The obtaining of the weight value of each parameter in the parameter set of the grid affecting the carbon dioxide leakage includes: Performing a first normalization process on a first type of parameters in a parameter set of the grid; performing a second normalization process on the second type of parameters in the parameter set of the grid; Based on the normalized data, a weight value of each parameter in the parameter set of the grid affecting the carbon dioxide leakage is obtained; The first type of parameters include porosity, water saturation, carbon dioxide injection pressure and carbon dioxide injection rate, and the second type of parameters include permeability.

3. The method according to claim 2, characterized in that The performing a first normalization process on the first type of parameters in the parameter set of the grid includes: For any parameter in the first category of parameters, obtain a value set of the parameter, a maximum value and a minimum value in the value set; Based on the value set of the parameter, the maximum value and the minimum value in the value set, a normalized value set of the parameter is obtained.

4. The method according to claim 2, characterized in that: The performing a second normalization process on the second type of parameters in the parameter set of the grid includes: For any parameter in the second category of parameters, obtaining a logarithmic value set of the parameter, and a maximum value and a minimum value in the logarithmic value set; Based on the logarithmic value set of the parameter, the maximum value and the minimum value in the logarithmic value set, a normalized value set of the parameter is obtained.

5. The method according to claim 2, characterized in that: The step of obtaining the weight value of each parameter in the grid parameter set affecting the carbon dioxide leakage based on the normalized parameter value includes: For any parameter, obtaining a value set of carbon dioxide leakage when the parameter is a first value; Determine a maximum carbon dioxide leakage value and an average carbon dioxide leakage value in the value set; Obtaining a derivative value set of carbon dioxide leakage when the parameter value is a normalized value set; Based on the average carbon dioxide leakage amount, the maximum carbon dioxide leakage amount and the derivative value set, a weight value of the parameter affecting the carbon dioxide leakage amount is obtained.

6. The method according to claim 1, characterized in that The step of establishing a numerical simulation model based on the three-dimensional geological model according to the type of the target area and the carbon dioxide storage mechanism includes: When the target area is an oil reservoir, a numerical simulation model based on the three-dimensional geological model is established according to the structural storage mechanism of carbon dioxide, the dissolution storage mechanism of carbon dioxide, the bound storage mechanism of carbon dioxide and the mineralization storage mechanism of carbon dioxide; wherein the numerical simulation model includes a component model for obtaining crude oil components.

7. The method according to any one of claims 1 to 6, characterized in that The method further comprises: After determining that the grid has a risk of carbon dioxide leakage, outputting a safety prompt; The safety prompt is used to instruct to monitor the carbon dioxide leakage in the grid.

8. A safety determination device for carbon dioxide storage, characterized in that: The device comprises: A first building module is configured to build a three-dimensional geological model of the target area and perform grid division; wherein the three-dimensional geological model includes a cap layer or multiple interlayers; A second establishment module is configured to establish a numerical simulation model based on the three-dimensional geological model according to the type of the target area and the carbon dioxide storage mechanism; A numerical simulation module is configured to perform numerical simulation on the process of injecting carbon dioxide into the target area for storage based on the numerical simulation model, and obtain a parameter set that affects the carbon dioxide leakage of each grid; wherein the parameter set includes at least: the porosity, permeability and water saturation of the reservoir in the target area, and the carbon dioxide injection pressure and carbon dioxide injection rate; An acquisition module is configured to acquire, for any grid, a weight value of each parameter in a parameter set of the grid affecting the carbon dioxide leakage; wherein the weight value is used to characterize the degree of influence of the corresponding parameter on the carbon dioxide leakage; A determination module is configured to determine the carbon storage risk coefficient of the grid based on the value of each parameter in the parameter set of the grid and the weight value corresponding to each parameter; when the carbon storage risk coefficient of the grid is higher than a preset threshold, it is determined that the grid has a risk of carbon dioxide leakage.

9. A computer device, characterized in that: The device comprises a processor and a memory, wherein the memory stores at least one program code, and the at least one program code is loaded and executed by the processor to implement the safety determination method for carbon dioxide sequestration as claimed in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: At least one program code is stored in the storage medium, and the at least one program code is loaded and executed by the processor to implement the safety determination method for carbon dioxide sequestration as claimed in any one of claims 1 to 7.

Images