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

Abstract

The application discloses a safety determination method and device for carbon dioxide storage, equipment and a storage medium, and belongs to the technical field of oil and gas field development. For carbon dioxide storage in a target area, the application performs sensitivity analysis on parameters affecting the carbon dioxide leakage amount on the basis of establishing a geological model and a numerical simulation model, that is, the application determines the influence degree of each parameter on the carbon dioxide leakage amount, and then obtains a carbon dioxide leakage risk evaluation index at different spatial positions of the target area, and determines the carbon dioxide leakage risk according to the index, so as to guide the on-site carbon dioxide leakage amount monitoring, for example, deploying a carbon dioxide leakage monitoring well at a spatial position with a high carbon dioxide leakage risk, realizing effective monitoring of the carbon dioxide leakage condition, ensuring that remedial measures are taken in time when leakage occurs, and thus the safety of the carbon dioxide storage is ensured.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of oil and gas field development, and particularly relates to a safety determination method and device for carbon dioxide storage, equipment and a storage medium. BACKGROUND

[0002] Carbon dioxide storage is a technology of storing carbon dioxide captured in industrial or energy production processes in the underground to reduce the concentration of carbon dioxide in the atmosphere. At present, carbon dioxide storage is considered as one of important means for reducing greenhouse gas emissions.

[0003] The storage site of carbon dioxide usually includes oil and gas reservoirs. When the concentration of carbon dioxide in the air exceeds 2%, it can cause damage to human respiratory organs, and more seriously, high-concentration carbon dioxide can even cause human suffocation and death due to oxygen deficiency. Therefore, before injecting carbon dioxide into the oil and gas reservoir for storage, the safety of carbon dioxide storage needs to be determined, that is, the safety of carbon dioxide storage needs to be evaluated to avoid large-scale leakage of carbon dioxide after storage, thereby causing casualties and economic losses on the ground. SUMMARY

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

[0005] In one aspect, a safety determination method for carbon dioxide storage is provided, and the method comprises:

[0006] establishing a three-dimensional geological model of a target region and performing grid division; wherein the three-dimensional geological model contains a cap rock or multiple interlayers;

[0007] establishing a numerical simulation model based on the three-dimensional geological model according to the type of the target region and the carbon dioxide storage mechanism;

[0008] performing numerical simulation on the process of injecting carbon dioxide into the target region for storage based on the numerical simulation model to obtain a parameter set of carbon dioxide leakage amount affecting each grid; wherein the parameter set at least includes the porosity, permeability and water saturation of the reservoir in the target region, and the carbon dioxide injection pressure and carbon dioxide injection speed;

[0009] for any one grid, obtaining a weight value of each parameter in the parameter set of the grid affecting the carbon dioxide leakage amount; wherein the weight value is used to represent the influence degree of the corresponding parameter on the carbon dioxide leakage amount;

[0010] determining a 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.

[0011] In a case where the carbon burial risk coefficient of the grid is higher than a preset threshold, it is determined that the grid has a risk of carbon dioxide leakage.

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

[0013] performing first normalization processing on a first type of parameter in the parameter set of the grid;

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

[0015] obtaining, based on the normalized data, the weight value of each parameter in the parameter set of the grid on the carbon dioxide leakage amount;

[0016] The first type of parameter comprises porosity, water saturation, carbon dioxide injection pressure, and carbon dioxide injection speed, and the second type of parameter comprises permeability.

[0017] In a possible implementation, the first normalization processing on the first type of parameter in the parameter set of the grid comprises:

[0018] For any one parameter in the first type of parameter, a value set of the parameter, a maximum value and a minimum value in the value set are obtained;

[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, the second normalization processing on the second type of parameter in the parameter set of the grid comprises:

[0021] For any one parameter in the second type of parameter, a logarithmic value set of the parameter, a maximum value and a minimum value in the logarithmic value set are obtained;

[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, the obtaining of the weight value of each parameter in the parameter set of the grid on the carbon dioxide leakage amount based on the normalized parameter value comprises:

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

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

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

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

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

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

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

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

[0032] The safety prompt is used to indicate that the grid is monitored for carbon dioxide leakage amount.

[0033] On the other hand, a carbon dioxide storage safety determination device is provided, comprising:

[0034] The first establishment module is configured to establish a three-dimensional geological model of a target region and perform grid division; wherein the three-dimensional geological model contains a cap rock or multiple interlayers;

[0035] The 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 region and the carbon dioxide storage mechanism;

[0036] The numerical simulation module is configured to perform numerical simulation on the process of carbon dioxide storage in the target region based on the numerical simulation model to obtain a parameter set affecting the carbon dioxide leakage amount of each grid; wherein the parameter set at least includes the porosity, permeability and water saturation of the reservoir in the target region, and the carbon dioxide injection pressure and the carbon dioxide injection speed;

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

[0038] The determination module is configured to determine a 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 the parameter, and determine that the grid has a carbon dioxide leakage risk in a case where the carbon sequestration risk coefficient of the grid is higher than a preset threshold.

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

[0040] perform first normalization processing on a first type of parameter in the parameter set of the grid;

[0041] perform second normalization processing on a second type of parameter in the parameter set of the grid;

[0042] acquire, based on the normalized data, the weight value of each parameter in the parameter set of the grid affecting the carbon dioxide leakage amount;

[0043] The first type of parameter includes porosity, water saturation, carbon dioxide injection pressure, and carbon dioxide injection speed, and the second type of parameter includes permeability.

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

[0045] For any one parameter in the first type of parameter, acquire a value set of the parameter, a maximum value and a minimum value in the value set;

[0046] acquire, 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.

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

[0048] For any one parameter in the second type of parameter, acquire a logarithmic value set of the parameter, a maximum value and a minimum value in the logarithmic value set;

[0049] acquire, 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.

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

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

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

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

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

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

[0056] In the case that the target area is an oil reservoir, a numerical simulation model based on the three-dimensional geological model is established according to carbon dioxide tectonic storage mechanism, carbon dioxide dissolution storage mechanism, carbon dioxide binding storage mechanism, and carbon dioxide mineralization storage mechanism; wherein the numerical simulation model includes a component model for obtaining a crude oil component.

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

[0058] The output module is configured to output a safety prompt after determining that the grid has a carbon dioxide leakage risk; wherein the safety prompt is used to indicate that the grid is monitored for carbon dioxide leakage amount.

[0059] In another aspect, a computer device is provided, the device including a processor and a memory, the memory having at least one program code stored therein, the at least one program code being loaded and executed by the processor to implement the above-mentioned carbon dioxide storage safety determination method.

[0060] In another aspect, a computer readable storage medium is provided, the storage medium having at least one program code stored therein, the at least one program code being loaded and executed by a processor to implement the above-mentioned carbon dioxide storage safety determination method.

[0061] In another aspect, a computer program product or computer program is provided, the computer program product or computer program including computer program code stored in a computer readable storage medium, the computer program code being read by a processor of a computer device from the computer readable storage medium, the processor executing the computer program code to cause the computer device to perform the above-mentioned carbon dioxide storage safety determination method.

[0062] For carbon dioxide storage in the target area, the embodiment of the application will also perform sensitivity analysis on the parameters affecting the carbon dioxide leakage amount on the basis of establishing the geological model and the numerical simulation model, that is, the embodiment of the application will determine the influence degree of each parameter on the carbon dioxide leakage amount, and then obtain the carbon dioxide leakage risk evaluation index (also referred to as the carbon storage risk coefficient) at different spatial positions of the target area, and determine the carbon dioxide leakage risk accordingly, so as to guide the on-site carbon dioxide leakage amount monitoring, such as deploying carbon dioxide leakage monitoring wells at spatial positions with high carbon dioxide leakage risk, to realize effective monitoring of the carbon dioxide leakage situation, ensure that remedial measures are taken in time when leakage occurs, and thus ensure the safety of carbon dioxide storage. BRIEF DESCRIPTION OF DRAWINGS

[0063] In order to more clearly illustrate the technical solutions in the embodiments of the application, the drawings needed in the embodiment description will be briefly introduced. Obviously, the drawings in the following description are only some embodiments of the application, and other drawings can be obtained by those skilled in the art without creative labor.

[0064] Figure 1 is a flowchart of a carbon dioxide storage safety determination method provided by the embodiment of the application;

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

[0066] Figure 3 is an analysis result of the influence of porosity on carbon dioxide leakage amount provided by the embodiment of the application;

[0067] Figure 4 is another analysis result of the influence of porosity on carbon dioxide leakage amount provided by the embodiment of the application;

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

[0069] Figure 6 is another analysis result of the influence of permeability on carbon dioxide leakage amount provided by the embodiment of the application;

[0070] Figure 7 is an analysis result of the influence of water saturation on carbon dioxide leakage amount provided by the embodiment of the application;

[0071] Figure 8 is another analysis result of the influence of water saturation on carbon dioxide leakage amount provided by the embodiment of the application;

[0072] Figure 9is an analysis result of the influence of carbon dioxide injection pressure on carbon dioxide leakage amount provided by an embodiment of the present application;

[0073] Figure 10 is another analysis result of the influence of carbon dioxide injection pressure on carbon dioxide leakage amount provided by an embodiment of the present application;

[0074] Figure 11 is an analysis result of the influence of carbon dioxide injection speed on carbon dioxide leakage amount provided by an embodiment of the present application;

[0075] Figure 12 is another analysis result of the influence of carbon dioxide injection speed on carbon dioxide leakage amount provided by an embodiment of the present application;

[0076] Figure 13 is a distribution diagram of carbon dioxide leakage risk at different spatial positions in a target area provided by an embodiment of the present application;

[0077] Figure 14 is a diagram of a suggested deployment position of a carbon dioxide leakage monitoring well provided by an embodiment of the present application;

[0078] Figure 15 is a structural diagram of a safety determination device for carbon dioxide storage provided by an embodiment of the present application;

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

[0080] In order to make the purpose, technical scheme and advantages of the present application clearer, the embodiments of the present application will be further described in detail below with reference to the drawings.

[0081] The terms "first", "second", and the like are used in the present application to distinguish between elements or items having substantially the same function and effect, and it should be understood that there is no logical or time sequence dependency between "first", "second", and "n", and the number and execution order are not limited. It should also be understood that although the following description uses the terms first, second, and the like 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. 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. The first element and the second element can both be elements, and in some cases, can be separate and distinct elements.

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

[0084] In this paper, "and / or" means that there can be three relationships, for example, A and / or B can mean: A exists alone, A and B exist together, and B exists alone. The character " / " generally represents that the associated objects before and after are 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 for analysis, stored data, displayed data, etc.) and signals involved in the present application are all authorized by the user or fully authorized by all parties, and the collection, use and processing of related data need to comply with relevant laws, regulations and standards in the relevant region.

[0086] The embodiment of the present application provides a carbon dioxide storage safety determination scheme, which gives a carbon dioxide storage safety evaluation method in the process of injecting carbon dioxide. The execution subject of the scheme is a computer device. Illustratively, the types of the above-mentioned computer device include but are not limited to desktop computers, notebook computers, smart phones and tablet computers, etc., which are not limited by the present application.

[0087] Among them, the embodiment of the present application evaluates the carbon storage safety of the target area (such as oil and gas reservoir) where carbon dioxide is injected, which can solve the following two problems: on the one hand, the related technology uses instrument equipment to monitor the leakage of carbon dioxide, which not only has high cost but also cannot realize long-term monitoring; on the other hand, at least solve the problem that the related technology cannot realize quantitative characterization of carbon dioxide storage safety.

[0088] Illustratively, the carbon storage safety evaluation method provided by the embodiment of the present application is mainly applied to the oil and gas reservoir where carbon dioxide is injected for storage. By calculating the carbon dioxide leakage risk evaluation index at different spatial positions of the target area, the quantitative evaluation of the carbon dioxide leakage risk is realized, and then the carbon dioxide leakage amount in the target area during the injection of carbon dioxide is monitored, for example, by deploying carbon dioxide leakage monitoring wells to monitor the carbon dioxide leakage amount, and then the safety of injecting carbon dioxide is ensured.

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

[0090] Figure 1is a flowchart of a method for determining the safety of carbon dioxide storage provided in an embodiment of the present application. The method is executed by a computer device. Referring to Figure 1 , the method flow includes:

[0091] 101. The computer device establishes a three-dimensional geological model of the target area and performs grid division; wherein the three-dimensional geological model contains a cap rock or multiple interlayers.

[0092] In an embodiment of the present application, the computer device uses a geological modeling software to establish a geological model of the target area. Taking the target area where carbon dioxide is injected as an X reservoir for example, referring to Figure 2 , the grid number of the geological model of the X reservoir is 84*38*16, and the average grid size is 65*65*10 meters; wherein 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 zones and two reservoir zones.

[0093] In addition, in the present embodiment, the pressure values of each grid in the geological model are represented by the number line shown in the lower right corner of Figure 2 , with the unit of kPa.

[0094] 102. The computer device establishes 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.

[0095] Exemplarily, the type of the above-mentioned target area includes but is not limited to oil and gas reservoirs, which is not limited in the present application. In an embodiment of the present application, the computer device uses a reservoir numerical simulation software to simulate the process of injecting carbon dioxide into the target area. In other words, the numerical simulation model lays a foundation for accurately simulating and predicting the leakage amount of carbon dioxide and evaluating the safety of carbon storage.

[0096] The first point to be explained is that in order to accurately simulate the process of injecting carbon dioxide into the target area, the carbon dioxide structural storage mechanism, the carbon dioxide dissolution storage mechanism (including the dissolution of carbon dioxide in oil and water), the carbon dioxide bound storage and the carbon dioxide mineralization storage mechanism are also 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, the interaction needs to be described. Therefore, if the target area is a reservoir, a component model for obtaining the components of crude oil also needs to be established. In other words, if the target area is a reservoir, the numerical simulation model includes a component model for obtaining the components of crude oil.

[0098] In summary, 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, including but not limited to the following ways:

[0099] In the case that the target region is an oil reservoir, a numerical simulation model based on a three-dimensional geological model is established according to a carbon dioxide structural storage mechanism, a carbon dioxide dissolution storage mechanism, a carbon dioxide binding storage mechanism and a carbon dioxide mineralization storage mechanism; wherein the numerical simulation model includes a component model for obtaining a crude oil component.

[0100] The carbon dioxide structural storage mechanism, the carbon dioxide dissolution storage mechanism, the carbon dioxide binding storage mechanism and the carbon dioxide mineralization storage mechanism are introduced as follows.

[0101] For the carbon dioxide structural storage, it is a storage mechanism based on the structural trap characteristics of the geological model to prevent further migration of the carbon dioxide.

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

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

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

[0105] For the carbon dioxide mineralization storage, the carbon dioxide mineralization includes two steps, carbonization (combination with the rock) and hydrolysis (reaction with water). Among them, the ionic reaction of the carbon dioxide in water is based on the instantaneous reaction model, and the dissolution and precipitation reaction of the mineral is based on the transition state theory.

[0106] 103. The computer device performs numerical simulation on the process of storing the injected carbon dioxide in the target region based on the established numerical simulation model, and obtains a parameter set of the carbon dioxide leakage amount affecting each grid; wherein the parameter set at least includes: the porosity, the permeability and the water saturation of the reservoir in the target region, and the injection pressure and the injection speed of the carbon dioxide.

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

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

[0109] In the process of performing the parameter sensitivity analysis, in order to avoid the influence of the value difference of different parameters on the carbon sequestration safety evaluation result, the embodiments of the present application will normalize each parameter in the parameter set. In the embodiments of the present application, the porosity, the water saturation, the carbon dioxide injection pressure and the carbon dioxide injection speed are subjected to first normalization processing, and the permeability is subjected to second normalization processing. In other words, the weight value of each parameter in the parameter set of the grid affecting the carbon dioxide leakage amount is obtained by using the following method:

[0110] The first type of parameters in the parameter set of the grid are subjected to first normalization processing, and the second type of parameters in the parameter set of the grid are subjected to second normalization processing; based on the normalized data, the weight value of each parameter in the parameter set of the grid affecting the carbon dioxide leakage amount is obtained; wherein the first type of parameters includes the porosity, the water saturation, the carbon dioxide injection pressure and the carbon dioxide injection speed, and the second type of parameters includes the permeability.

[0111] Taking the first normalization processing as the maximum-minimum normalization as an example, the first normalization processing of the first type of parameters in the parameter set of the grid includes:

[0112] For any one parameter in the first type of parameters, the value set of the parameter, the maximum value and the minimum value in the value set of the parameter are obtained; based on the value set of the parameter, the maximum value and the minimum value in the value set of the parameter, the normalized value set of the parameter is obtained.

[0113] The above process is described by mathematical expression as formula (1) below.

[0114]

[0115] Wherein, a is a parameter set affecting the carbon dioxide leakage amount, the parameter set at least including the porosity, the permeability, the water saturation, the carbon dioxide injection pressure and the carbon dioxide injection speed of the reservoir in the target region; x a is the value set of the parameter a; x norm，a is the normalized value set of the parameter a; min(x a ) is the minimum value in the value set of the parameter a; max(x a ) is the maximum value in the value set of the parameter a.

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

[0117] For any one of the second type of parameters, a logarithmic value set of the parameter, a maximum value and a minimum value in the logarithmic value set of the parameter are obtained; based on the logarithmic value set of the parameter, the maximum value and the minimum value in the logarithmic value set of the parameter, a normalized value set of the parameter is obtained.

[0118] The above process is described by a mathematical expression as follows.

[0119]

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

[0121] In addition, the embodiments of the present application obtain the weight value of each parameter in the parameter set of the grid by Figures 3 to 12 The analysis results of the sensitivity analysis on the porosity, the permeability, the water saturation, the carbon dioxide injection pressure and the carbon dioxide injection speed are shown respectively.

[0122] In a possible implementation, the embodiments of the present application select the porosity, the permeability, the water saturation, the carbon dioxide injection pressure and the carbon dioxide injection speed of the reservoir to perform the sensitivity analysis on the carbon dioxide leakage amount. In other words, the weight value of each parameter in the parameter set of the grid affecting the carbon dioxide leakage amount is obtained by using the following method:

[0123] For any one parameter, a value set of the carbon dioxide leakage amount when the parameter is the first value is obtained; the maximum value and the average value of the carbon dioxide leakage amount in the value set are determined; a derivative value set of the carbon dioxide leakage amount when the value of the parameter is the normalized value set is obtained; based on the average value of the carbon dioxide leakage amount, the maximum value of the carbon dioxide leakage amount and the derivative value set of the carbon dioxide leakage amount, the weight value of the parameter affecting the carbon dioxide leakage amount is obtained.

[0124] The above process is described by a mathematical expression as follows.

[0125]

[0126] wherein, w(a) is the weight value set of the parameter a affecting the carbon dioxide leakage amount; L(x a ) is the value of the parameter a taking xa a value of the carbon dioxide leakage amount when the parameter a takes the value x; max(L(x a )) is a maximum value of the carbon dioxide leakage amount in the value set of the carbon dioxide leakage amount when the parameter a takes the value x a a value of the carbon dioxide leakage amount when the parameter a takes the value x; L'(x a ) is an average value of the carbon dioxide leakage amount in the value set of the carbon dioxide leakage amount when the parameter a takes the value x norm a value of the carbon dioxide leakage amount when the parameter a takes the value x; L'(x norm,a ) is a derivative value of the carbon dioxide leakage amount when the parameter a takes the value x norm，a a value of the carbon dioxide leakage amount when the parameter a takes the value x; L'(x .

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

[0128] 105、The computer device determines a 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; and determines that the grid has a carbon dioxide leakage risk in a case where the carbon sequestration risk coefficient of the grid is higher than a preset threshold.

[0129] In the embodiment of the present application, the weighted method is used to calculate the carbon dioxide leakage risk evaluation index, i.e. the carbon sequestration risk coefficient, at different spatial positions of the target region based on the value of each parameter in the parameter set at different spatial positions of the target region and the weight value corresponding to each parameter, and the larger the value is, the higher the carbon dioxide leakage risk is. The calculation formula of the carbon dioxide leakage risk evaluation index is as follows:

[0130]

[0131] wherein, I is the carbon dioxide leakage risk evaluation index, i.e. the carbon sequestration risk coefficient; and n is the number of parameters affecting the carbon dioxide leakage amount.

[0132] Exemplarily, the calculation result of the carbon dioxide leakage risk evaluation index is as shown in Figure 13 . That is, Figure 13 shows the distribution of the carbon dioxide leakage risk at different spatial positions in the target region. In addition, the value of the carbon dioxide leakage risk evaluation index of each grid in the geological model is represented by the number axis shown in the lower right corner of Figure 13 .

[0133] In a possible implementation, after determining that a certain spatial position in the target region has a carbon dioxide leakage risk, the computer device outputs a safety prompt; wherein the safety prompt is used to instruct the grid to monitor the carbon dioxide leakage amount. That is, the computer device guides the field to monitor the carbon dioxide leakage amount according to the carbon dioxide leakage risk. For example, a carbon dioxide leakage monitoring well is deployed at a spatial position with a high carbon dioxide leakage risk, so as to effectively monitor the carbon dioxide leakage and ensure that remedial measures are taken in time when leakage occurs, thereby ensuring the safety of carbon dioxide storage. Figure 14

[0134] The embodiment of the present application provides a carbon dioxide storage safety quantitative evaluation scheme with low cost and considering the carbon dioxide storage mechanism. In detail, for carbon dioxide storage of a target region, the scheme further performs sensitivity analysis on parameters affecting the carbon dioxide leakage amount on the basis of establishing a geological model and a numerical simulation model, that is, the scheme determines the influence degree of each parameter on the carbon dioxide leakage amount, and then obtains a carbon dioxide leakage risk evaluation index of different spatial positions of the target region, and determines the carbon dioxide leakage risk according to the index, so as to guide the field to monitor the carbon dioxide leakage amount, for example, a carbon dioxide leakage monitoring well is deployed at a spatial position with a high carbon dioxide leakage risk, so as to effectively monitor the carbon dioxide leakage and ensure that remedial measures are taken in time when leakage occurs, thereby ensuring the safety of carbon dioxide storage.

[0135] Figure 15 FIG. 1 is a structural schematic diagram of a carbon dioxide storage safety determination device provided by the embodiment of the present application. Referring to FIG. 1, the device comprises: Figure 15

[0136] A first establishing module 1501 is configured to establish a three-dimensional geological model of a target region and perform grid division; wherein the three-dimensional geological model contains a cap rock or multiple sets of interlayers;

[0137] A 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 region and a carbon dioxide storage mechanism;

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

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

[0140] The determination module 1505 is configured to determine a 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 the parameter, and determine that the grid has a carbon dioxide leakage risk in a case where the carbon sequestration risk coefficient of the grid is higher than a preset threshold.

[0141] The embodiment of the present application provides a carbon dioxide sequestration safety quantitative evaluation scheme with low cost and considering the carbon dioxide sequestration mechanism in the industry. In detail, for carbon dioxide sequestration of a target region, the scheme will perform sensitivity analysis on parameters affecting the carbon dioxide leakage amount on the basis of establishing a geological model and a numerical simulation model, that is, the scheme will determine the influence degree of each parameter on the carbon dioxide leakage amount, and then obtain a carbon dioxide leakage risk evaluation index at different spatial positions of the target region, and determine the carbon dioxide leakage risk accordingly, to guide the on-site carbon dioxide leakage amount monitoring, such as deploying a carbon dioxide leakage monitoring well at a spatial position with a high carbon dioxide leakage risk, to realize effective monitoring of the carbon dioxide leakage condition, ensure that remedial measures are taken in time when leakage occurs, and thus ensure the safety of carbon dioxide sequestration.

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

[0143] perform first normalization processing on a first type of parameter in the parameter set of the grid;

[0144] perform second normalization processing on a second type of parameter in the parameter set of the grid;

[0145] acquire, based on the normalized data, a weight value of each parameter in the parameter set of the grid affecting the carbon dioxide leakage amount;

[0146] The first type of parameter includes porosity, water saturation, carbon dioxide injection pressure, and carbon dioxide injection speed, and the second type of parameter includes permeability.

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

[0148] For any one parameter in the first type of parameter, acquire a value set of the parameter, a maximum value and a minimum value in the value set;

[0149] Based on the parameter value set, 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 obtaining module is configured to:

[0151] For any one of the second type of parameters, a logarithmic value set of the parameter, a maximum value and a minimum value in the logarithmic value set are obtained.

[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 obtaining module is configured to:

[0154] For any one of the second type of parameters, a logarithmic value set of the parameter, a maximum value and a minimum value in the logarithmic value set are obtained.

[0155] A maximum value and an average value of the carbon dioxide leakage amount in the value set are determined.

[0156] A derivative value set of the carbon dioxide leakage amount when the value of the parameter is a normalized value set is obtained.

[0157] Based on the average value of the carbon dioxide leakage amount, the maximum value of the 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 establishing module is configured to:

[0159] In a case where the target region is an oil reservoir, a numerical simulation model based on the three-dimensional geological model is established according to a carbon dioxide structural storage mechanism, a carbon dioxide dissolution storage mechanism, a carbon dioxide binding storage mechanism and a carbon dioxide mineralization storage mechanism; wherein the numerical simulation model includes a component model for obtaining a crude oil component.

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

[0161] The output module is configured to output a safety prompt after determining that the grid has a carbon dioxide leakage risk; wherein the safety prompt is used to instruct to perform carbon dioxide leakage amount monitoring on the grid.

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

[0163] It should be noted that the above embodiment provides the safety determination device for carbon dioxide storage. In the process of injecting carbon dioxide, the safety determination of carbon dioxide storage is determined. Only the above-mentioned division of each functional module is used as an example for illustration. In actual application, the above-mentioned functions can be completed by different functional modules according to the needs, 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 and the safety determination method for carbon dioxide storage provided in the above embodiment belong to the same concept, and the specific implementation process is described in the method embodiment. Here, it is not repeated.

[0164] Figure 16 is a structural schematic diagram of a computer device 1600 provided by an embodiment of the present application.

[0165] Generally, 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 of a hardware form of a DSP (Digital Signal Processing), an FPGA (Field-Programmable Gate Array), and a 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 an awake state, also known as a CPU (Central Processing Unit). The coprocessor is a low-power processor for processing data in a standby state. In a possible implementation manner, the processor 1601 is integrated with a GPU (Graphics Processing Unit) that is responsible for rendering and drawing the content required to be displayed by the display screen. In a possible implementation manner, the processor 1601 further includes an AI (Artificial Intelligence) processor that is used to process machine learning-related computing operations.

[0167] The memory 1602 includes one or more computer-readable storage media. The computer-readable storage media is non-transitory. The memory 1602 includes high-speed random access memory and can also include non-volatile memory, such as one or more magnetic disk storage devices, optical disk storage devices, flash memory devices, or other non-volatile solid-state storage devices. In one possible implementation, the non-transitory computer-readable storage medium in the memory 1602 is used for storing the at least one program code for being executed by the processor 1601 to implement the method for determining the safety of carbon dioxide storage provided by the method embodiments in the present application.

[0168] In one 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 through a bus or a signal line. Each peripheral device is connected to the peripheral device interface 1603 through 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 component 1606, an audio circuit 1607, a positioning component 1608, and a power supply 1609.

[0169] The peripheral device interface 1603 is used for connecting the at least one peripheral device related to input / output (I / O) 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, and the present application does not limit this.

[0170] The radio frequency circuit 1604 is configured to receive and send RF (Radio Frequency) signals, also known as electromagnetic signals. The radio frequency circuit 1604 communicates with communication networks and other communication devices through electromagnetic signals. The radio frequency circuit 1604 converts electrical signals into electromagnetic signals for transmission, or converts electromagnetic signals received into electrical signals. In a 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 coding chipset, a subscriber 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 a possible implementation, the radio frequency circuit 1604 also includes NFC (Near Field Communication) related circuitry, which is not limited in the present application.

[0171] The display screen 1605 is configured to display a UI (User Interface). The UI includes graphics, text, icons, video, 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 or above the surface of the display screen 1605. The touch signals are input to the processor 1601 as control signals for processing. At this time, the display screen 1605 is also configured to provide virtual buttons and / or a virtual keyboard, also known as soft buttons and / or soft keyboard. In a possible implementation, the display screen 1605 is one, arranged on the front panel of the computer device 1600; in another possible implementation, the display screen 1605 is at least two, arranged on different surfaces of the computer device 1600 or in a folding design; in another possible implementation, the display screen 1605 is a flexible display screen, arranged on a curved surface or a folding surface of the computer device 1600. Alternatively, the display screen 1605 is arranged in an irregular shape other than a rectangle, i.e., a special-shaped screen. The display screen 1605 is made of materials such as LCD (Liquid Crystal Display), OLED (Organic Light-Emitting Diode), and the like.

[0172] The camera assembly 1606 is configured to capture images or videos. In a possible implementation, the camera assembly 1606 includes a front camera and a rear camera. Generally, the front camera is arranged on the front panel of the terminal, and the rear camera is arranged on the back of the terminal. In a possible implementation, the rear camera is at least two, which is any one of a main camera, a depth-of-field camera, a wide-angle camera, and a long-focus camera, to realize the background blur function by fusing the main camera and the depth-of-field camera, realize the panoramic shooting and VR (Virtual Reality) shooting function by fusing the main camera and the wide-angle camera, or realize other fusion shooting functions. In another possible implementation, the camera assembly 1606 further includes a flash. The flash is a single-color-temperature flash or a double-color-temperature flash. The double-color-temperature flash refers to a combination of a warm light flash and a cold light flash, which is used for light compensation under different color temperatures.

[0173] The audio circuit 1607 includes a microphone and a speaker. The microphone is configured to collect sound waves of a user and an environment, and convert the sound waves into an electrical signal input to the processor 1601 for processing or input to the radio frequency circuit 1604 to realize voice communication. For the purpose of stereo sound collection or noise reduction, the microphone is multiple, which is arranged at different parts of the computer device 1600. Alternatively, the microphone is an array microphone or an omnidirectional collection microphone. The speaker is configured to convert an electrical signal from the processor 1601 or the radio frequency circuit 1604 into sound waves. The speaker is a traditional diaphragm speaker or a piezoelectric ceramic speaker. In the case where the speaker is a piezoelectric ceramic speaker, not only can the electrical signal be converted into a sound wave audible to humans, but also can be converted into an inaudible sound wave to humans for ranging purposes. In a possible implementation, the audio circuit 1607 further includes a headphone jack.

[0174] The positioning component 1608 is configured to position the current geographical position of the computer device 1600 to realize navigation or LBS (Location Based Service). The positioning component 1608 can be a positioning component based on the GPS (Global Positioning System) of the United States, the Beidou system of China, the GLONASS system of Russia, or the Galileo system of the European Union.

[0175] The power supply 1609 is configured to supply power to each component 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 charging battery or a wireless charging battery. The wired charging battery is a battery charged through a wired line, and the wireless charging 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 can understand that, Figure 16 The structure shown in the figure does not constitute a limitation on the computer device 1600, and can include more or fewer components than the figure, or combine certain components, or adopt a different component arrangement.

[0177] In an exemplary embodiment, a computer readable storage medium, such as a memory including program code executable by a processor in a computer device to perform the above-described method for determining the safety of carbon dioxide storage is also provided. 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, and an optical data storage device, etc.

[0178] In an exemplary embodiment, a computer program product or computer program is also provided, which includes computer program code stored in a computer readable storage medium, and 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 to cause the computer device to perform the above-described method for determining the safety of carbon dioxide storage.

[0179] Those of ordinary skill in the art can understand that all or part of the steps of the above-described embodiments can be completed by hardware, or by a program instructing relevant hardware, and the program can be stored in a computer readable storage medium, such as a Read-Only Memory, a magnetic disk or an optical disk, etc.

[0180] The above is only an optional embodiment of the present application, and is not intended to limit the present application. Any modification, equivalent replacement, improvement, 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 by, The method comprises: establishing a three-dimensional geological model of a target area and performing grid division; wherein the three-dimensional geological model contains a cap rock or multiple interlayers; 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; based on the numerical simulation model, performing numerical simulation on the process of carbon dioxide injection and storage in the target area to obtain a parameter set of carbon dioxide leakage amounts 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 speed; for any one grid, obtaining a weight value of each parameter in the parameter set of the grid affecting the carbon dioxide leakage amount; wherein the weight value is used to represent the influence degree of the corresponding parameter on the carbon dioxide leakage amount; based on the value of each parameter in the parameter set of the grid and the weight value corresponding to each parameter, determining the carbon storage risk coefficient of the grid; in the case where the carbon storage risk coefficient of the grid is higher than a preset threshold, determining that the grid has a carbon dioxide leakage risk.

2. The method of claim 1, wherein, The method comprises: performing first normalization processing on the first type of parameters in the parameter set of the grid; performing second normalization processing on the second type of parameters in the parameter set of the grid; based on the normalized data, obtaining the weight value of each parameter in the parameter set of the grid affecting the carbon dioxide leakage amount; wherein the first type of parameters includes porosity, water saturation, carbon dioxide injection pressure and carbon dioxide injection speed, and the second type of parameters includes permeability.

3. The method of claim 2, wherein, The method comprises: for any one parameter in the first type of parameters, obtaining 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, obtaining a normalized value set of the parameter.

4. The method of claim 2, wherein, The method comprises: for any one parameter in the second type of parameters, obtaining a logarithmic value set of the parameter, 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, obtaining a normalized value set of the parameter.

5. The method of claim 2, wherein, The method comprises: for any one parameter, obtaining a value set of the carbon dioxide leakage amount when the parameter is a first value; determining the maximum value and the average value of the carbon dioxide leakage amount in the value set; obtaining a derivative value set of the carbon dioxide leakage amount when the value of the parameter is a normalized value set; Based on the average value of the carbon dioxide leakage amount, the maximum value of the 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 of claim 1, wherein, The numerical simulation model based on the three-dimensional geological model is established according to the type of the target region and the carbon dioxide storage mechanism, and includes: In the case that the target region is an oil reservoir, a numerical simulation model based on the three-dimensional geological model is established according to a carbon dioxide structural storage mechanism, a carbon dioxide dissolution storage mechanism, a carbon dioxide binding storage mechanism, and a carbon dioxide mineralization storage mechanism; wherein the numerical simulation model includes a component model for obtaining a crude oil component.

7. The method according to any one of claims 1 to 6, characterized in that, The method further includes: After determining that the grid has a carbon dioxide leakage risk, a safety prompt is outputted; The safety prompt is used to indicate that the grid is monitored for carbon dioxide leakage amount.

8. A safety determination device for carbon dioxide storage, characterized by comprising: The device includes: A first establishment module configured to establish a three-dimensional geological model of a target region and perform grid division; wherein the three-dimensional geological model contains a cap rock or multiple sets of interlayers; A second establishment module configured to establish a numerical simulation model based on the three-dimensional geological model according to the type of the target region and the carbon dioxide storage mechanism; A numerical simulation module configured to perform numerical simulation on the process of carbon dioxide storage in the target region based on the numerical simulation model to obtain a parameter set affecting the carbon dioxide leakage amount of each grid; wherein the parameter set at least includes the porosity, permeability, and water saturation of the reservoir in the target region, and the carbon dioxide injection pressure and carbon dioxide injection speed; An acquisition module configured to obtain, for any one grid, a weight value of each parameter in the parameter set of the grid affecting the carbon dioxide leakage amount; wherein the weight value is used to represent the influence degree of the corresponding parameter on the carbon dioxide leakage amount; A determination module configured to determine a 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; in the case that the carbon storage risk coefficient of the grid is higher than a preset threshold, it is determined that the grid has a carbon dioxide leakage risk.

9. A computer device, comprising: The device includes a processor and a memory, and the memory stores at least one program code, which is loaded and executed by the processor to implement the safety determination method for carbon dioxide storage according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, The storage medium stores at least one program code, which is loaded and executed by the processor to implement the safety determination method for carbon dioxide storage according to any one of claims 1 to 7.

Images