A method and system for evaluating the stability of salt cavern carbon pools based on multi-source data

Abstract

This invention discloses a method and system for evaluating the stability of salt cavern carbon storage ponds based on multi-source data, belonging to the field of stability evaluation technology. The method involves collecting geological survey data of the area where the salt cavern carbon storage pond is located and constructing a geological structure model; outputting temperature and pressure variation parameter ranges based on historical meteorological data; extracting abnormal operating pressure data of the salt cavern; combining the temperature and pressure variation parameter ranges with the abnormal operating pressure data, and incorporating the cavity burial depth data to set the duration and depth range of the impact; defining a disturbance input vector and constructing a response evolution rate index; constructing a response evolution intensity based on the response evolution rate index; setting a dynamic risk value calculation function based on the response evolution intensity; and pre-setting a threshold range to analyze and provide stability early warning. This method not only achieves a comprehensive evaluation of salt cavern stability from static structural analysis to dynamic evolution response but also significantly improves risk identification capabilities and decision-making efficiency under complex environments and extreme operating conditions.

Description

Technical Field

[0001] This invention relates to the field of stability evaluation technology, specifically to a method and system for evaluating the stability of salt cavern carbon storage pools based on multi-source data. Background Technology

[0002] Currently, salt caverns are mainly used for underground storage of energy sources such as natural gas and hydrogen. With the rise of carbon storage technology, their application potential in carbon dioxide sequestration is becoming increasingly apparent. Meanwhile, rock salt strata possess unique mechanical characteristics such as creep, plasticity, and low permeability, making the structural stability of salt caverns a critical issue for engineering safety under disturbances such as high-pressure injection and extraction. Related research is gradually focusing on constructing refined geological models, simulating multi-source disturbance scenarios, and evaluating rock mass response mechanisms. Preliminary evaluation frameworks based on finite element simulation, discrete element analysis, and risk matrices have been established. For example, the Analytic Hierarchy Process (AHP) combined with geological parameters (such as salt mine core data and interlayer distribution) is used for stability evaluation, or FLAC3D software is used for long-term creep simulation of gas storage in salt caverns (such as the stability analysis based on the Mohr-Coulomb model in the Zhaoji Salt Mine case), providing certain technical support for engineering practice.

[0003] However, existing methods for evaluating the carbon storage stability of salt caverns generally suffer from problems such as limited data dimensions, idealized boundary settings, and coarse characterization of disturbance response mechanisms. For example, while the analytic hierarchy process (AHP) can achieve a qualitative-to-quantitative transformation through multi-factor weight allocation, it lacks the ability to respond in real time to dynamic disturbances (such as sudden temperature and pressure changes, and variations in injection and production frequency). Numerical tools such as FLAC3D can simulate long-term creep behavior, but they struggle to integrate multi-source data (such as meteorological changes and operational anomalies) for dynamic risk modeling. In most studies, only a single dimension of formation structure or injection and production parameters is considered for analysis, lacking the ability to integrate multi-source information (such as meteorological data, historical operational data, and geological exploration data) for modeling, making it difficult to accurately reflect the nonlinear response behavior of salt caverns under multiple disturbances in complex environments. Furthermore, existing methods often use ideal operating conditions to simulate sudden temperature and pressure changes, lacking boundary construction mechanisms based on historical measured data, and failing to dynamically evaluate the coupling relationship between disturbance intensity and structural evolution. More critically, most current model evaluation indicators fail to incorporate attenuation control for factors such as disturbance depth and duration, resulting in limited accuracy in risk warning. Summary of the Invention

[0004] The purpose of this invention is to provide a method and system for evaluating the stability of salt cavern carbon storage pools based on multi-source data, so as to solve the problems mentioned in the background art.

[0005] To solve the above-mentioned technical problems, the present invention provides the following technical solution:

[0006] A method for evaluating the stability of salt cavern carbon reservoirs based on multi-source data is disclosed. The method includes the following steps: Step S1: Collect geological survey data of the area where the salt cavern carbon reservoir is located and construct a geological structure model; Step S2: Output the temperature and pressure variation parameter ranges based on historical meteorological data; extract abnormal operating pressure data of the salt caverns; combine the temperature and pressure variation parameter ranges with the abnormal operating pressure data, and combine this with the cavity burial depth data to set the duration and depth range of the action; Step S3: Define a disturbance input vector and construct a response evolution rate index; construct the response evolution intensity based on the response evolution rate index; Step S4: Based on the response evolution intensity, set a dynamic risk value calculation function; preset a threshold range, analyze and issue a stability warning.

[0007] As a preferred embodiment of the multi-source data-based method for evaluating the stability of salt cavern carbon reservoirs described in this invention, geological survey data of the area where the salt cavern carbon reservoir is located is collected. This geological survey data includes salt mine core data, stratigraphic structure data, interlayer distribution data, cavity depth data, and stratigraphic dip angle data. Based on this geological survey data, a geological structure model is constructed. This geological structure model is used to reflect the mechanical differences and stress sensitivity characteristics of the rock salt layer, as detailed below:

[0008] Based on the salt mine core data, record the lithology types at different depths and output preliminary vertical lithology distribution information. Based on the vertical lithology distribution information and combined with stratigraphic structure data, output the stratigraphic structure profile.

[0009] In the rock strata structure profile, based on the interlayer distribution data, a rock strata structure model containing interlayer information is output; based on the rock strata structure model and cavity burial depth data, the depth range where the cavity is located is marked, and the relationship between the cavity embedding position and the surrounding rock composition is output;

[0010] Based on the relationship between the cavity embedding position and the surrounding rock composition, and combined with the stratum dip angle data, a three-dimensional spatial structure diagram is constructed; in the three-dimensional spatial structure diagram, corresponding lithology and structural parameters are set for different rock layers and interlayers to construct a geological structure model.

[0011] As a preferred embodiment of the multi-source data-based method for evaluating the stability of salt cavern carbon storage as described in this invention, based on historical meteorological data, temperature and pressure records with abrupt changes are extracted, the magnitude and rate of change are recorded, and the range of temperature and pressure change parameters is output; in the historical data of salt cavern operation, abnormal operating pressure data that occurred in a short period of time are extracted, including sudden increases in operating pressure or changes in injection and extraction frequency, and the corresponding operating condition change magnitude and duration are output.

[0012] By combining the temperature and pressure change parameter ranges with the abnormal operating pressure data, several sets of simulated boundary conditions are generated to reflect different types of temperature and pressure change scenarios. Combined with the cavity burial depth data, the duration and depth range of action for various temperature and pressure change conditions are set as the boundary control basis for subsequent simulations. The temperature change, pressure change, action time, and action depth are uniformly organized into a set of input parameters.

[0013] As a preferred embodiment of the salt cavern carbon pool stability evaluation method based on multi-source data described in this invention, a perturbation input vector is defined as follows: in, ΔT represents the comprehensive characteristics of the i-th perturbation event. i Let ΔP represent the temperature change caused by the i-th disturbance event. i τ represents the pressure change caused by the i-th disturbance event. i d represents the duration of the i-th disturbance event. i This indicates the depth of effect of the i-th perturbation event;

[0014] The response evolution rate metric is constructed as follows:

[0015]

[0016] Where α and β represent preset weighting factors for temperature and pressure changes, γ1 and γ2 represent preset sensitivity factors (controlling the degree of nonlinear enhancement (such as the sensitivity of crack propagation to temperature or pressure changes)), λ represents preset application time weighting factor, and μ represents preset depth attenuation coefficient.

[0017] It should be noted that the first term of this formula represents the combined effect of thermal and pressure disturbances, and instead of using linear weighting, it employs an exponential form for nonlinear amplification. This reflects the characteristics that in high-temperature / high-pressure zones, effects such as material crack propagation and interface instability are more likely to occur and grow faster. The intermediate term simulates that the longer the disturbance lasts, the stronger the destructive effect on the rock mass, but it will not increase indefinitely. This reflects the physical law that the cumulative response has a saturation upper limit, such as the gradual weakening of the time effect of material creep or damage. The last term is based on the mechanism of gradual attenuation of thermal and pressure disturbances propagating downwards in actual underground structures. The response at depth will be much smaller than that at shallow depths, which is consistent with the basic laws of thermal diffusion and stress wave attenuation. The parameters can be flexibly adjusted according to different salt cavern geological structures, burial depths, and pressure / temperature histories, possessing good versatility and adaptability, and can be extended to different types of salt cavern carbon reservoir evaluation tasks.

[0018] Based on the aforementioned response evolution rate index, the response evolution intensity is constructed as follows:

[0019]

[0020] Among them, w i Let represent the influence coefficient of the i-th disturbance event, and Λ represent the response evolution intensity under the current disturbance history.

[0021] As a preferred embodiment of the salt cavern carbon pool stability evaluation method based on multi-source data described in this invention, a dynamic risk value calculation function is set based on the response evolution intensity Λ, as follows:

[0022] Ω=κ1·Λ+κ2·ΔT max +κ3·ΔP max

[0023] Where Λ represents the intensity of the response evolution under the current disturbance history, ΔT max ΔP represents the maximum temperature change in the current period. max κ1, κ2, and κ3 represent the maximum pressure change in the current cycle, κ1, κ2, and κ3 represent the preset empirical control parameters, characterizing the contribution of each factor to stability, and Ω represents the dynamic risk value.

[0024] A preset risk threshold range is defined. If the dynamic risk value Ω is less than the minimum value within the risk threshold range, the salt cavern is determined to have excellent structural stability under continuous disturbance.

[0025] If the dynamic risk value Ω is within the risk threshold range, it is determined that the salt cavern is unstable under continuous disturbance. If the dynamic risk value Ω is greater than the maximum value within the risk threshold range, it is determined that the structural stability of the salt cavern is poor under continuous disturbance, and a stability warning is issued.

[0026] A stability evaluation system for salt cavern carbon storage pools based on multi-source data. The system includes: a data acquisition and model building module, a data acquisition and processing module, an index and intensity calculation module, and a risk calculation, analysis, and early warning module.

[0027] The data acquisition and model building module collects geological survey data of the area where the salt cavern carbon storage pond is located and builds a geological structure model.

[0028] The data acquisition and processing module outputs the temperature and pressure change parameter range based on historical meteorological data; extracts abnormal operating pressure data of the salt cavern; combines the temperature and pressure change parameter range with the abnormal operating pressure data, and, in conjunction with the cavity burial depth data, sets the duration and depth range of action.

[0029] The index and intensity calculation module: defines the perturbation input vector and constructs the response evolution rate index; based on the response evolution rate index, it constructs the response evolution intensity;

[0030] The risk calculation, analysis, and early warning module: based on the response evolution intensity, sets a dynamic risk value calculation function; presets a threshold range, analyzes and provides stability early warning.

[0031] Furthermore, the data acquisition and model building module includes a data acquisition and model building unit;

[0032] The data acquisition and model building unit collects geological survey data of the area where the salt cavern carbon reservoir is located. This geological survey data includes salt mine core data, stratigraphic structure data, interlayer distribution data, cavity depth data, and stratigraphic dip angle data. Based on this geological survey data, a geological structure model is constructed. This geological structure model is used to reflect the mechanical differences and stress sensitivity characteristics of the rock salt layer, as detailed below:

[0033] Based on the salt mine core data, record the lithology types at different depths and output preliminary vertical lithology distribution information. Based on the vertical lithology distribution information and combined with stratigraphic structure data, output the stratigraphic structure profile.

[0034] In the rock strata structure profile, based on the interlayer distribution data, a rock strata structure model containing interlayer information is output; based on the rock strata structure model and cavity burial depth data, the depth range where the cavity is located is marked, and the relationship between the cavity embedding position and the surrounding rock composition is output;

[0035] Based on the relationship between the cavity embedding position and the surrounding rock composition, and combined with the stratum dip angle data, a three-dimensional spatial structure diagram is constructed; in the three-dimensional spatial structure diagram, corresponding lithology and structural parameters are set for different rock layers and interlayers to construct a geological structure model.

[0036] Furthermore, the data acquisition and processing module includes a data acquisition and processing unit;

[0037] The data acquisition and processing unit: Based on historical meteorological data, extracts temperature and pressure records with abrupt changes, records the magnitude and rate of change, and outputs the range of temperature and pressure change parameters; In the historical data of salt cavern operation, extracts abnormal operating pressure data that occurred in a short period of time, including sudden increases in operating pressure or changes in injection and extraction frequency, and outputs the corresponding operating condition change magnitude and duration.

[0038] By combining the temperature and pressure change parameter ranges with the abnormal operating pressure data, several sets of simulated boundary conditions are generated; by combining the cavity burial depth data, the duration and depth range of action are set for various temperature and pressure abrupt change conditions; and the temperature change, pressure change, action time, and action depth are uniformly organized into a set of input parameters.

[0039] Furthermore, the index and intensity calculation module includes an index calculation unit and an intensity calculation unit;

[0040] The index calculation unit defines a disturbance input vector, denoted as... in, ΔT represents the comprehensive characteristics of the i-th perturbation event. i Let ΔP represent the temperature change caused by the i-th disturbance event. i τ represents the pressure change caused by the i-th disturbance event. i d represents the duration of the i-th disturbance event. i Indicate the depth of influence of the i-th perturbation event; construct a response evolution rate index;

[0041] The intensity calculation unit: constructs the response evolution intensity Λ based on the response evolution rate index.

[0042] Furthermore, the risk calculation and analysis early warning module includes a risk calculation unit and an analysis early warning unit;

[0043] The risk calculation unit: sets a dynamic risk value calculation function based on the response evolution intensity Λ;

[0044] The analysis and early warning unit has a preset risk threshold range. If the dynamic risk value Ω is less than the minimum value within the risk threshold range, the salt cavern is determined to have good structural stability under continuous disturbance. If the dynamic risk value Ω is within the risk threshold range, the salt cavern is determined to be unstable under continuous disturbance. If the dynamic risk value Ω is greater than the maximum value within the risk threshold range, the salt cavern is determined to have poor structural stability under continuous disturbance, and a stability warning is issued.

[0045] Compared with existing technologies, the beneficial effects achieved by this invention are as follows: The method and system for evaluating the stability of salt cavern carbon reservoirs based on multi-source data provided by this invention first constructs an accurate geological structure model, providing a physical basis for subsequent stability simulation; by integrating historical meteorological data and salt cavern operational anomaly data, various disturbance boundary conditions are set, and the time and depth of action are considered, achieving comprehensive simulation of complex operating conditions; further, disturbance input vectors and response evolution indices are introduced to quantify disturbance intensity and capture its spatiotemporal evolution characteristics, establishing a cumulative characterization system for the historical impact of disturbances; finally, a dynamic stability early warning mechanism is constructed through a risk value calculation function and graded threshold intervals. The overall technical solution not only achieves a comprehensive evaluation of salt cavern stability from static structural analysis to dynamic evolution response, but also significantly improves risk identification capabilities and decision-making efficiency under complex environments and extreme operating conditions, possessing good engineering adaptability and promotional value. Attached Figure Description

[0046] The accompanying drawings are provided to further illustrate the invention and form part of the specification. They are used together with the embodiments of the invention to explain the invention and do not constitute a limitation thereof.

[0047] Figure 1 This is a schematic diagram illustrating the steps of a method for evaluating the stability of salt cavern carbon storage pools based on multi-source data according to the present invention.

[0048] Figure 2 This is a schematic diagram of the structure of a salt cavern carbon storage stability evaluation system based on multi-source data according to the present invention. Detailed Implementation

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

[0050] Please see Figure 1 In this first embodiment, a method for evaluating the stability of salt cavern carbon storage pools based on multi-source data is provided. This method includes the following steps:

[0051] Step S1: Collect geological survey data of the area where the salt cavern carbon storage pond is located and construct a geological structure model.

[0052] Specifically, geological survey data of the area where the salt cavern carbon reservoir is located is collected. This data includes salt mine core data, stratigraphic structure data, interlayer distribution data, cavity depth data, and stratigraphic dip angle data. Based on this data, a geological structure model is constructed to reflect the mechanical differences and stress-sensitive characteristics of the rock salt layers, as detailed below:

[0053] Based on salt mine core data, lithological types at different depths are recorded, and preliminary vertical lithological distribution information is output. Based on the vertical lithological distribution information and combined with stratigraphic structure data, a rock strata structure profile is output. In the actual construction of salt cavern carbon storage ponds, lithology at different depths has a significant impact on the stability of the salt cavern. For example, in some areas, the roof and interlayers of salt mines are composed of silty mudstone, which poses a risk of leakage when storing high-pressure CO2 gas. Therefore, accurately grasping the lithological distribution is crucial.

[0054] In the rock strata structure profile, based on the interlayer distribution data, a rock strata structure model containing interlayer information is output; based on the rock strata structure model and cavity burial depth data, the depth range of the cavity is marked, and the relationship between the cavity embedment location and the surrounding rock is output; for example, in the design process of a salt cavern gas storage facility, the relationship between the cavity and the surrounding rock was determined by analyzing the interlayer distribution and cavity burial depth, providing a basis for subsequent stability evaluation.

[0055] Based on the relationship between the cavity embedding position and the surrounding rock composition, and combined with the stratum dip angle data, a three-dimensional spatial structure diagram is constructed; in the three-dimensional spatial structure diagram, corresponding lithology and structural parameters are set for different rock layers and interlayers to construct a geological structure model.

[0056] In this invention, by collecting multi-dimensional geological exploration data, including salt mine core data, stratigraphic structure data, interlayer distribution data, cavity burial depth data, and stratigraphic dip angle data, a detailed geological structure model is constructed. This model comprehensively depicts the structure, lithological characteristics, and mechanical differences of the rock-salt layers surrounding the salt cavern carbon reservoir. Furthermore, this model can be used to analyze the contact relationship and stress state between the cavity and the surrounding rock, and is particularly instructive for analyzing the stability and stress transmission path of the salt cavern. This step provides a physical basis for subsequent thermo-compressional disturbance response and risk analysis, improves the accuracy and reliability of the simulation boundary, and thus enhances the scientific rigor of the overall stability assessment.

[0057] Step S2: Based on historical meteorological data, output the temperature and pressure change parameter range; extract abnormal operating pressure data of the salt cavern; combine the temperature and pressure change parameter range with the abnormal operating pressure data, and combine the cavity burial depth data to set the duration and depth range of action.

[0058] Specifically, based on historical meteorological data, temperature and pressure records with abrupt changes are extracted, the magnitude and rate of change are recorded, and the range of temperature and pressure change parameters is output. In the historical data of salt cavern operation, abnormal operating pressure data that occurred in a short period of time are extracted. The abnormal operating pressure data includes sudden increases in operating pressure or changes in injection and extraction frequency, and the corresponding magnitude and duration of the operating condition changes are output. For example, in the actual operation of salt caverns, there may be sudden increases in operating pressure due to equipment failure, or changes in injection and extraction frequency due to adjustments in production plans. These abnormal situations need to be considered in the evaluation.

[0059] By combining the temperature and pressure variation parameter ranges with abnormal operating pressure data, several sets of simulated boundary conditions are generated to reflect different types of temperature and pressure abrupt change scenarios. Based on the cavity burial depth data, the duration and depth range of each temperature and pressure abrupt change condition are set as the boundary control basis for subsequent simulations. Temperature changes, pressure changes, duration, and depth are uniformly organized into a unified set of input parameters. In the actual operation of salt cavern carbon storage tanks, cavities at different depths respond differently to temperature and pressure changes. For example, beyond a depth of 800m, the CO2 density does not change significantly with depth, but the impact of temperature and pressure changes on cavity stability still exists. This step, by setting boundary conditions in conjunction with the cavity burial depth, better reflects the actual situation.

[0060] In this invention, multiple simulated boundary conditions are generated by combining external natural meteorological data (such as records of sudden temperature and pressure changes) with internal operating conditions (such as sudden pressure increases or changes in injection and extraction frequencies), achieving comprehensive modeling of various real-world disturbance scenarios. Furthermore, by combining cavity burial depth data to define the temporal and spatial range of the disturbance's impact, the simulation input possesses spatial hierarchy and temporal evolution. This step enables the model to not only reflect stability under normal operating conditions but also evaluate dynamic stability responses under extreme conditions (such as heavy rain, extreme heat, and sudden operational changes), thereby significantly improving the adaptability and foresight of risk analysis.

[0061] Step S3: Define the perturbation input vector and construct the response evolution rate index; based on the response evolution rate index, construct the response evolution intensity.

[0062] Specifically, the perturbation input vector is defined as . in, ΔT represents the comprehensive characteristics of the i-th perturbation event. i Let ΔP represent the temperature change caused by the i-th disturbance event. i τ represents the pressure change caused by the i-th disturbance event. i d represents the duration of the i-th disturbance event. i This indicates the depth of effect of the i-th perturbation event;

[0063] The response evolution rate metric is constructed as follows:

[0064]

[0065] Where α and β represent the preset weighting factors for temperature and pressure changes, γ1 and γ2 represent the preset sensitivity factors (controlling the degree of nonlinear enhancement (such as the sensitivity of crack propagation to temperature or pressure changes)), λ represents the preset weighting factor for the duration of action, and μ represents the preset depth attenuation coefficient. In the actual process of carbon storage in salt caverns, temperature and pressure changes will cause changes in the mechanical properties of salt rocks. For example, the creep characteristics of salt rocks are significantly affected by temperature and pressure. This formula comprehensively considers the influence of these factors on the stability of salt caverns by introducing weighting factors and sensitivity coefficients.

[0066] Based on the aforementioned response evolution rate index, the response evolution intensity is constructed as follows:

[0067]

[0068] Among them, w i Let represent the influence coefficient of the i-th disturbance event, and Λ represent the response evolution intensity under the current disturbance history.

[0069] In this invention, this step quantifies the perturbation process into a perturbation input vector (including temperature change, pressure change, duration of action, and depth of action), and introduces weighting factors and sensitivity coefficients to construct a response evolution rate index for nonlinear enhancement characteristics. By integrating temperature and pressure factors with spatiotemporal dimensions, a mathematical expression that can measure the degree of perturbation impact is formed. Furthermore, the response evolution rate is integrated into the historical impact intensity of the perturbation, constituting an index system reflecting the cumulative effect of multiple rounds of perturbations. This method innovatively links transient perturbations with long-term risks, providing quantitative support for risk trend prediction, thereby realizing the transformation from single-point analysis to evolutionary analysis and enhancing the ability to judge the long-term stability of salt caverns.

[0070] Step S4: Based on the response evolution intensity, set the dynamic risk value calculation function; preset the threshold range, analyze and issue stability warnings.

[0071] Specifically, based on the response evolution intensity Λ, a dynamic risk value calculation function is defined as follows:

[0072] Ω=κ1·Λ+κ2·ΔT max +κ3·ΔP max

[0073] Where Λ represents the intensity of the response evolution under the current disturbance history, ΔT max ΔP represents the maximum temperature change in the current period. max κ1, κ2, and κ3 represent the maximum pressure change in the current cycle, κ1, κ2, and κ3 represent the preset empirical control parameters, characterizing the contribution of each factor to stability, and Ω represents the dynamic risk value. In practical applications, different salt cavern carbon storage ponds have different stability requirements due to different geological conditions, cavity structures, and other factors. By setting empirical control parameters κ1, κ2, and κ3, the evaluation needs of different salt caverns can be adapted more flexibly.

[0074] A preset risk threshold range is defined. If the dynamic risk value Ω is less than the minimum value within the risk threshold range, the salt cavern is determined to have excellent structural stability under continuous disturbance.

[0075] If the dynamic risk value Ω is within the risk threshold range, it is determined that the salt cavern is unstable under continuous disturbance. If the dynamic risk value Ω is greater than the maximum value within the risk threshold range, it is determined that the structural stability of the salt cavern is poor under continuous disturbance, and a stability warning is issued.

[0076] It should be noted that, based on the existing evolution intensity, a dynamic risk value function containing multiple control parameters is set and compared with a preset risk threshold range to achieve the judgment and early warning of the stability of the salt cavern structure under the current disturbance state. By integrating the maximum temperature and pressure disturbance parameter and the evolution intensity index, the dynamic risk value comprehensively reflects the risk level of the current state. This step not only realizes the quantitative output of risk, but also improves the real-time performance and decision-making of the risk response of the salt cavern carbon storage system by introducing a threshold range, thus ensuring the safe and controllable operation of the carbon sequestration system.

[0077] Please see Figure 2 In this second embodiment: a stability evaluation system for salt cavern carbon storage pools based on multi-source data is provided. The system includes: a data acquisition and model building module, a data acquisition and processing module, an index and intensity calculation module, and a risk calculation, analysis, and early warning module.

[0078] The data acquisition and model building module collects geological survey data of the area where the salt cavern carbon storage pond is located and builds a geological structure model.

[0079] The data acquisition and processing module outputs the temperature and pressure change parameter range based on historical meteorological data; extracts abnormal operating pressure data of the salt cavern; combines the temperature and pressure change parameter range with the abnormal operating pressure data, and, in conjunction with the cavity burial depth data, sets the duration and depth range of action.

[0080] The index and intensity calculation module: defines the perturbation input vector and constructs the response evolution rate index; based on the response evolution rate index, it constructs the response evolution intensity;

[0081] The risk calculation, analysis, and early warning module: based on the response evolution intensity, sets a dynamic risk value calculation function; presets a threshold range, analyzes and provides stability early warning.

[0082] Furthermore, the data acquisition and model building module includes a data acquisition and model building unit;

[0083] The data acquisition and model building unit collects geological survey data of the area where the salt cavern carbon reservoir is located. This geological survey data includes salt mine core data, stratigraphic structure data, interlayer distribution data, cavity depth data, and stratigraphic dip angle data. Based on this geological survey data, a geological structure model is constructed. This geological structure model is used to reflect the mechanical differences and stress sensitivity characteristics of the rock salt layer, as detailed below:

[0084] Based on the salt mine core data, record the lithology types at different depths and output preliminary vertical lithology distribution information. Based on the vertical lithology distribution information and combined with stratigraphic structure data, output the stratigraphic structure profile.

[0085] In the rock strata structure profile, based on the interlayer distribution data, a rock strata structure model containing interlayer information is output; based on the rock strata structure model and cavity burial depth data, the depth range where the cavity is located is marked, and the relationship between the cavity embedding position and the surrounding rock composition is output;

[0086] Based on the relationship between the cavity embedding position and the surrounding rock composition, and combined with the stratum dip angle data, a three-dimensional spatial structure diagram is constructed; in the three-dimensional spatial structure diagram, corresponding lithology and structural parameters are set for different rock layers and interlayers to construct a geological structure model.

[0087] Furthermore, the data acquisition and processing module includes a data acquisition and processing unit;

[0088] The data acquisition and processing unit: Based on historical meteorological data, extracts temperature and pressure records with abrupt changes, records the magnitude and rate of change, and outputs the range of temperature and pressure change parameters; In the historical data of salt cavern operation, extracts abnormal operating pressure data that occurred in a short period of time, including sudden increases in operating pressure or changes in injection and extraction frequency, and outputs the corresponding operating condition change magnitude and duration.

[0089] By combining the temperature and pressure change parameter ranges with the abnormal operating pressure data, several sets of simulated boundary conditions are generated; by combining the cavity burial depth data, the duration and depth range of action are set for various temperature and pressure abrupt change conditions; and the temperature change, pressure change, action time, and action depth are uniformly organized into a set of input parameters.

[0090] Furthermore, the index and intensity calculation module includes an index calculation unit and an intensity calculation unit;

[0091] The index calculation unit defines a disturbance input vector, denoted as... in, ΔT represents the comprehensive characteristics of the i-th perturbation event. i Let ΔP represent the temperature change caused by the i-th disturbance event. i τ represents the pressure change caused by the i-th disturbance event. i d represents the duration of the i-th disturbance event. i Indicate the depth of influence of the i-th perturbation event; construct a response evolution rate index;

[0092] The intensity calculation unit: constructs the response evolution intensity Λ based on the response evolution rate index.

[0093] Furthermore, the risk calculation and analysis early warning module includes a risk calculation unit and an analysis early warning unit;

[0094] The risk calculation unit: sets a dynamic risk value calculation function based on the response evolution intensity Λ;

[0095] The analysis and early warning unit has a preset risk threshold range. If the dynamic risk value Ω is less than the minimum value within the risk threshold range, the salt cavern is determined to have good structural stability under continuous disturbance. If the dynamic risk value Ω is within the risk threshold range, the salt cavern is determined to be unstable under continuous disturbance. If the dynamic risk value Ω is greater than the maximum value within the risk threshold range, the salt cavern is determined to have poor structural stability under continuous disturbance, and a stability warning is issued.

[0096] It should be noted that, in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such process, method, article, or apparatus.

[0097] Finally, it should be noted that the above descriptions are merely preferred embodiments of the present invention and are not intended to limit the present invention. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art can still modify the technical solutions described in the foregoing embodiments or make equivalent substitutions for some of the technical features. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A method for evaluating stability of a salt cavern carbon storage based on multi-source data, characterized in that, The method comprises the following steps: Step S1: Collecting geological survey data of the region where the salt cavern carbon storage is located, and constructing a geological structure model; Step S2: According to historical meteorological data, output the temperature and pressure change parameter range; extract the abnormal data of the operating pressure of the salt cavern; combine the temperature and pressure change parameter range with the abnormal data of the operating pressure, and set the action time length and the action depth range in combination with the cavity burial depth data; Step S3: Defining a disturbance input vector and constructing a response evolution rate index; based on the response evolution rate index, constructing a response evolution intensity; Step S4: Based on the response evolution intensity, setting a dynamic risk value calculation function; presetting a threshold interval, analyzing and performing stability warning; The specific implementation process of step S3 comprises: The disturbance input vector is defined as wherein, represents the comprehensive feature of the i-th disturbance action event, ΔT i represents the temperature change amount of the i-th disturbance action event, ΔP i represents the pressure change amount of the i-th disturbance action event, τ i represents the action time of the i-th disturbance action event, d i represents the action depth of the i-th disturbance action event; The response evolution rate index is constructed as follows: Wherein, α and β represent the preset weight factors of temperature change and pressure change, γ1 and γ2 represent the preset sensitivity factors, λ represents the preset action time weight factor, and μ represents the preset depth attenuation coefficient; Based on the response evolution rate index, the response evolution intensity is constructed as follows: where w i represents the influence coefficient of the ith disturbance action event, and represents the response evolution intensity under the current disturbance history. The specific implementation process of step S4 comprises: Based on the response evolution intensity Λ, the dynamic risk value calculation function is set as follows: Ω = κ1 · Λ + κ2 · ΔT + κ3 · ΔP max max ​ wherein, Λ represents the response evolution intensity under the current disturbance history, ΔT max represents the maximum temperature change in the current cycle, ΔP max represents the maximum pressure change in the current cycle, κ1, κ2 and κ3 represent preset empirical control parameters, which represent the contribution degree of each factor to stability, and Ω represents a dynamic risk value. The risk threshold interval is preset, if the dynamic risk value Ω is less than the minimum value in the risk threshold interval, it is determined that the structural stability of the salt cavern under continuous disturbance is good; If the dynamic risk value Ω exists in the risk threshold interval, it is determined that the salt cavern under continuous disturbance is unstable, and if the dynamic risk value Ω is greater than the maximum value in the risk threshold interval, it is determined that the structural stability of the salt cavern under continuous disturbance is poor, and then a stability warning is issued.

2. The salt cavern carbon repository stability evaluation method based on multi-source data according to claim 1, characterized in that, The specific implementation process of step S1 comprises: Collecting the geological survey data of the region where the salt cavern carbon storage is located, the geological survey data comprising salt mine core data, stratum structure data, interlayer distribution data, cavity burial depth data and stratum inclination data, and based on the geological survey data, constructing a geological structure model, the geological structure model being used to reflect the mechanical difference and stress sensitivity characteristic of rock salt layer, and the specific implementation process being as follows: According to the salt mine core data, the lithology type at different depth positions is recorded, and the preliminary lithology vertical distribution information is outputted, and based on the lithology vertical distribution information, the stratum structure profile is outputted in combination with the stratum structure data; In the stratum structure profile, the stratum structure model containing interlayer information is outputted according to the interlayer distribution data; based on the stratum structure model and the cavity burial depth data, the depth interval where the cavity is located is marked, and the cavity embedding position and surrounding rock constituting relationship are outputted; Based on the cavity embedding position and surrounding rock constituting relationship, the three-dimensional space structure diagram is constructed in combination with the stratum inclination data; in the three-dimensional space structure diagram, the corresponding lithology and structure parameters of different strata and interlayers are set, and the geological structure model is constructed.

3. The method for evaluating the stability of salt cavern carbon storage according to claim 2, wherein, The specific implementation process of step S2 comprises: According to historical meteorological data, temperature and pressure records with mutation characteristics are extracted, the change amplitude and change rate are recorded, and the temperature and pressure change parameter range is output; in the operation history data of the salt cave, operation pressure abnormal data occurring in a short time are extracted, the operation pressure abnormal data include operation pressure sudden rise or injection-production frequency change, and the corresponding working condition change amplitude and duration are output; The temperature and pressure change parameter range and the operation pressure abnormal data are combined to generate a plurality of groups of simulated boundary conditions; the temperature change amount, the pressure change amount, the action time, and the action depth are unified and arranged as an input parameter set in combination with the cavity burial depth data.

4. A salt cavern carbon storage stability evaluation system based on multi-source data, which executes a salt cavern carbon storage stability evaluation method based on multi-source data according to any one of claims 1-3, characterized in that, The system comprises a data acquisition and model construction module, a data acquisition and processing module, an index and intensity calculation module, and a risk calculation and analysis warning module; The data acquisition and model construction module collects geological survey data of the region where the salt cave carbon storage is located, and constructs a geological structure model; The data acquisition and processing module outputs the temperature and pressure change parameter range according to historical meteorological data; extracts the salt cave operation pressure abnormal data; combines the temperature and pressure change parameter range and the operation pressure abnormal data, and sets the action time length and the action depth range in combination with the cavity burial depth data; The index and intensity calculation module defines a disturbance input vector and constructs a response evolution rate index; based on the response evolution rate index, a response evolution intensity is constructed; The risk calculation and analysis warning module sets a dynamic risk value calculation function based on the response evolution intensity; a threshold interval is preset, and stability warning is analyzed and performed.

5. The salt cavern carbon repository stability evaluation system based on multi-source data according to claim 4, characterized in that: The data acquisition and model construction module comprises a data acquisition and model construction unit; The data acquisition and model construction unit collects geological survey data of the region where the salt cave carbon storage is located, and constructs a geological structure model based on the geological survey data, wherein the geological survey data comprises salt mine core data, stratum structure data, interlayer distribution data, cavity burial depth data, and stratum inclination data, and the geological structure model is used to reflect the mechanical difference and stress sensitivity characteristics of rock salt layer, and specifically as follows: According to the salt mine core data, the lithology type at different depth positions is recorded, and the preliminary lithology longitudinal distribution information is output; based on the lithology longitudinal distribution information, the stratum structure profile is output in combination with the stratum structure data; In the stratum structure profile, the stratum structure model containing interlayer information is output according to the interlayer distribution data; based on the stratum structure model and the cavity burial depth data, the depth interval where the cavity is located is marked, and the cavity embedding position and surrounding rock composition relationship are output; Based on the cavity embedding position and surrounding rock composition relationship, a three-dimensional space structure diagram is constructed in combination with the stratum inclination data; In the three-dimensional space structure diagram, corresponding lithology and structure parameters are set for different rock layers and interlayers, and the geological structure model is constructed.

6. The salt cavern carbon repository stability evaluation system based on multi-source data according to claim 5, characterized in that: The data acquisition and processing module comprises a data acquisition and processing unit; The data acquisition and processing unit: according to historical meteorological data, extract the temperature and pressure records with mutation characteristics, record the change amplitude and change rate, output the temperature and pressure change parameter range; in the operation history data of salt cavern, extract the operation pressure abnormal data occurring in a short time, the operation pressure abnormal data includes operation pressure sudden rise or injection production frequency change, output the corresponding working condition change amplitude and duration; Combine the temperature and pressure change parameter range with the operation pressure abnormal data to generate a plurality of groups of simulation boundary conditions; combine the cavity burial depth data to set the action time length and action depth range for various temperature and pressure mutation conditions; unify the temperature change amount, pressure change amount, action time and action depth into an input parameter set.

7. The salt cavern carbon repository stability evaluation system based on multi-source data according to claim 6, characterized in that: The index and intensity calculation module includes an index calculation unit and an intensity calculation unit; The index calculation unit: define the disturbance input vector, denoted as Wherein, The comprehensive characteristics of the i th disturbance event, ΔT i The temperature change of the i th disturbance event, ΔP i The pressure change of the i th disturbance event, τ i The action time of the i th disturbance event, d i The action depth of the i th disturbance event; Construct the response evolution rate index; The intensity calculation unit: based on the response evolution rate index, construct the response evolution intensity Lambda.

8. The salt cavern carbon repository stability evaluation system based on multi-source data according to claim 7, characterized in that: The risk calculation and analysis warning module includes a risk calculation unit and an analysis warning unit; The risk calculation unit: based on the response evolution intensity Lambda, set a dynamic risk value calculation function; The analysis warning unit: preset a risk threshold interval, if the dynamic risk value Omega is less than the minimum value in the risk threshold interval, it is determined that the structural stability of the salt cavern under continuous disturbance is good; If the dynamic risk value Omega exists in the risk threshold interval, it is determined that the salt cavern under continuous disturbance exists instability, if the dynamic risk value Omega is greater than the maximum value in the risk threshold interval, it is determined that the structural stability of the salt cavern under continuous disturbance is poor, then the stability warning is issued.

Images