Earthquake monitoring method, device and equipment based on structural parameters of gas storage and medium

Abstract

The embodiment of the invention provides an earthquake monitoring method, device and equipment based on gas storage structure parameters and a medium, and relates to the field of safety analysis technology and earthquake risk assessment based on the gas storage structure parameters. The method comprises the following steps: collecting seismic geological data, engineering geology and structural design information of a gas storage building area, generating a plurality of groups of random seismic oscillation time history records as input conditions, building a simulation model and carrying out dynamic response calculation to obtain structural mechanical response information; obtaining an earthquake risk cost by adopting a risk quantitative evaluation model, creating a corresponding risk control measure, and obtaining a disposal method library under different intensity seismic oscillation conditions; constructing an earthquake risk assessment model, and determining risk control measures; in the operation stage of the gas storage, the operation state of the gas storage is monitored in real time according to the three-dimensional geologic model obtained based on the structural parameters of the gas storage, potential earthquake risks are found and processed in time, and long-term safe operation of the gas storage is ensured.

Description

Technical Field

[0001] Embodiments of the present invention relate to the technical field of earthquake monitoring, and in particular to an earthquake monitoring method, device, equipment and medium based on gas storage reservoir structural parameters. Background Art

[0002] The Northwest region has abundant wind and solar energy resources, which account for a significant proportion of the country and is an important base for the development of new energy. The large-scale grid connection of new energy power generation has an urgent need for large-capacity energy storage power stations, while the Northwest region lacks water resources. Compressed air energy storage power stations are the most promising large-capacity energy storage methods in the Northwest. Considering the geological conditions in the Northwest, it is suitable to build a compressed air energy storage power station with supporting artificial underground gas storage. At the same time, the Northwest is located in China's strong earthquake zone and belongs to a high earthquake intensity area. Its history of seismic activity and geological structural characteristics determine that the region has a high earthquake risk. Artificial underground gas storage is a newly excavated building. At the same time, the sealing method of artificial underground gas storage is developing towards a flexible sealing system, which has the advantages of low cost and easy repair. However, the flexible sealing system has higher requirements for structural stability and airtightness. Its structural stability and airtightness determine the long-term operation safety of compressed air energy storage power stations. It is necessary to form an earthquake safety risk assessment method for artificial underground gas storage to provide support for risk management plans.

[0003] In the existing technology, the prediction of earthquake hazards mainly focuses on determining the earthquake source, calculating the relationship between the earthquake amplitude and frequency of the earthquake source, calculating the relationship between the ground peak acceleration of the underground gas storage and the earthquake frequency, etc., but there is no feasible and effective implementation method for quantitative evaluation of post-disaster risks after the actual earthquake and the designation of corresponding risk management plans. In addition, when an earthquake occurs, it is impossible to quickly assess the risk cost of the disaster according to the scale of the earthquake in the first time and make corresponding post-earthquake disposal responses. Summary of the invention

[0004] Based on the above technical problems, the embodiments of the present invention provide an earthquake monitoring method, device, equipment and medium based on the structural parameters of a gas storage reservoir, aiming to evaluate the earthquake risk cost of the gas storage reservoir through an earthquake risk assessment model and output corresponding risk control measures.

[0005] A first aspect of an embodiment of the present invention provides a method for earthquake monitoring based on gas storage reservoir structural parameters, the method comprising: Collecting seismic geological data, engineering geological and structural design information of the gas storage construction area, wherein the seismic geological data include: historical regional earthquake records and seismic motion parameters of the gas storage construction area, and the engineering geological and structural design information include: regional rock and soil layer distribution characteristics, groundwater conditions, geostress field, rock and soil physical and mechanical parameters, gas storage body design parameters, and gas storage structure parameters; Based on the historical regional earthquake records and seismic parameters of the gas storage construction area, multiple groups of random seismic time history records with different exceedance probabilities are generated; Based on the seismic geological data and the engineering geological and structural design information, a three-dimensional numerical simulation model is established; The generated multiple groups of random earthquake time history records under different exceedance probabilities are used as input conditions to perform dynamic response calculation of the three-dimensional numerical simulation model to obtain damage degree information, wherein the damage degree information includes: the damage range of the sealing layer and the damage range of the lining; Using a risk quantitative assessment model, according to the exceedance probability in the input conditions and the damage degree information under the earthquake conditions, the earthquake risk cost under the input conditions is obtained; Based on the earthquake risk cost, corresponding risk control measures are created to obtain a library of disposal methods under different intensities of earthquake motion conditions; Based on multiple sets of random seismic time history records with different exceedance probabilities and the corresponding earthquake risk costs, an earthquake risk assessment model is constructed, and risk control measures matching different earthquake risk costs are determined in the disposal method library; During the operation of the current gas storage, monitoring the earthquake motion in the area where the current gas storage is located and the stress level of the gas storage structure; When an earthquake occurs, the earthquake risk cost of the current gas storage facility is quickly assessed through the earthquake risk assessment model, and corresponding risk control measures are output.

[0006] Optionally, the generating multiple groups of random earthquake time history records under different exceedance probabilities are used as input conditions to perform dynamic response calculation of the three-dimensional numerical simulation model to obtain damage degree information, including: Based on the generated multiple groups of random earthquake time history records under different exceedance probabilities as input conditions, dynamic response calculation of the three-dimensional numerical simulation model is performed, and the stress of the sealing layer and the lining under each group of the input conditions is recorded; Among them, the allowable stress of the sealing layer material is σ Rs , the allowable stress of the lining material is σ Rc In each group of calculation results, the sealing layer stress is greater than the allowable stress σ Rs The total volume of the region V s , is determined as the damage range of the sealing layer; the lining stress is greater than σ Rc The total volume of the region V c , determined as the damage range of the lining.

[0007] Optionally, the earthquake risk cost includes: repair cost and production loss; using a risk quantitative evaluation model, according to the exceedance probability in the input conditions and the damage degree information under the earthquake conditions, the earthquake risk cost under the input conditions is obtained, including: Determine the repair cost C according to the damage range of the sealing layer and the damage range of the lining. repair , the repair cost C repair The amount of work used for reaction repair works; According to the number of days of shutdown and the daily benefit of the energy storage power station, the loss of production reduction C is determined. loss , the number of days of suspension of production is related to the amount of work of the repair project; P f ×(C repair +C loss ) is determined as the earthquake risk cost, where P f is the probability of occurrence of earthquake intensity at the site.

[0008] Optionally, an earthquake risk assessment model is constructed according to the multiple groups of random earthquake motion time history records under different exceedance probabilities and corresponding earthquake risk costs, including: The multiple groups of random seismic time history records under different exceedance probabilities and the corresponding earthquake risk costs are used as training samples to train the machine learning model to obtain the earthquake risk assessment model.

[0009] Optionally, based on historical regional earthquake records and seismic parameters in the gas storage construction area, multiple groups of random seismic time history records with different exceedance probabilities are generated, including: Based on the historical regional earthquake records of the gas storage construction area, determine the upper and lower limits of the ground motion intensity PGA; For the target PGA, based on the Monte Carlo simulation method, the multiple groups of random seismic time history records under different exceedance probabilities are generated with PGA as the target value.

[0010] Optionally, based on the earthquake risk cost, corresponding risk control measures are created to obtain a library of disposal methods under earthquake motion conditions of different intensities, including: Based on the earthquake risk cost, risk control measures for corresponding engineering quantities are created, and the risk control measures include: strengthening structures, sealing and repairing, and increasing emergency reserves; Risk control measures created for different earthquake risk costs form a library of disposal methods under earthquake motion conditions of different intensities.

[0011] A second aspect of an embodiment of the present invention provides an earthquake monitoring device based on gas storage reservoir structural parameters, the device comprising: A data collection module is used to collect seismic geological data, engineering geological and structural design information of the gas storage construction area, wherein the seismic geological data include: historical regional earthquake records and seismic motion parameters of the gas storage construction area, and the engineering geological and structural design information include: regional rock and soil layer distribution characteristics, groundwater conditions, ground stress field, rock and soil body physical and mechanical parameters, gas storage body design parameters, and gas storage structure parameters; An input condition generation module is used to generate multiple groups of random seismic time history records under different exceedance probabilities based on historical regional earthquake records and seismic parameters in the gas storage construction area; A simulation model building module, used to build a three-dimensional numerical simulation model based on the seismic geological data and the engineering geological and structural design information; A damage degree determination module is used to use the generated multiple groups of random seismic time history records under different exceedance probabilities as input conditions to perform dynamic response calculation of the three-dimensional numerical simulation model to obtain damage degree information, wherein the damage degree information includes: the damage range of the sealing layer and the damage range of the lining; A risk quantitative evaluation module, used to obtain the earthquake risk cost under the input conditions according to the exceedance probability in the input conditions and the damage degree information under the earthquake conditions; The disposal method library establishment module is used to create corresponding risk control measures based on earthquake risk costs and obtain a disposal method library under different intensities of earthquake motion conditions; The earthquake risk assessment module is used to construct an earthquake risk assessment model based on multiple sets of random earthquake motion time history records with different exceedance probabilities and the corresponding earthquake risk costs, and to determine the risk control measures matching different earthquake risk costs in the disposal method library; A monitoring module, used to monitor the seismic motion in the area where the current gas storage is located and the stress level of the gas storage structure during the operation of the current gas storage; The risk control measures output module is used to evaluate the earthquake risk cost of the current gas storage reservoir through an earthquake risk assessment model and output corresponding risk control measures.

[0012] Optionally, the generated multiple groups of random earthquake time history records under different exceedance probabilities are used as input conditions to perform dynamic response calculation of the three-dimensional numerical simulation model to obtain damage degree information; The damage degree determination module comprises: A dynamic response calculation submodule is used to perform dynamic response calculation of a three-dimensional numerical simulation model based on the generated multiple groups of random seismic time history records under different exceedance probabilities as input conditions, and record the stress of the sealing layer and the lining under each group of input conditions; The damage range determination submodule is used to determine the damage range of the sealing layer stress greater than the allowable stress σ in each group of calculation results. Rs The total volume of the region V s , is determined as the damage range of the sealing layer; the lining stress is greater than σ Rc The total volume of the region V c , determined as the damage range of the lining.

[0013] Optionally, the earthquake risk cost includes: repair cost and production loss; a risk quantitative evaluation model is used to obtain the earthquake risk cost under the input conditions according to the exceedance probability in the input conditions and the damage degree information under the earthquake conditions; The risk quantitative assessment module includes: The repair cost determination submodule is used to determine the repair cost C according to the damage range of the sealing layer and the damage range of the lining. repair , the repair cost C repair The amount of work used for reaction repair works; The production reduction loss determination submodule is used to determine the production reduction loss C according to the number of days of shutdown and the daily benefit of the energy storage power station. loss , the number of days of suspension of production is related to the amount of work of the repair project; The risk cost determination submodule converts P f ×(C repair +C loss ) is determined as the earthquake risk cost, where P f is the probability of occurrence of earthquake intensity at the site.

[0014] Optionally, an earthquake risk assessment model is constructed according to the multiple groups of random earthquake motion time history records under different exceedance probabilities and corresponding earthquake risk costs; The earthquake risk assessment module comprises: The assessment model training submodule is used to use the multiple groups of random seismic time history records under different exceedance probabilities and the corresponding earthquake risk costs as training samples to train the machine learning model to obtain the earthquake risk assessment model.

[0015] Optionally, based on historical regional earthquake records and seismic parameters in the gas storage construction area, multiple groups of random seismic time history records with different exceedance probabilities are generated; Optionally, the input condition generating module includes: A historical record confirmation submodule, used to determine the upper and lower limits of the ground motion intensity PGA based on the historical regional earthquake records in the gas storage construction area; The model generation submodule is used to generate the multiple groups of random seismic time history records under different exceedance probabilities for the target PGA based on the Monte Carlo simulation method and taking PGA as the target value.

[0016] Optionally, based on the earthquake risk cost, corresponding risk control measures are created to obtain a library of disposal methods under earthquake motion conditions of different intensities; The disposal method library establishment module includes: A risk control measure determination submodule is used to create risk control measures for corresponding engineering quantities based on the earthquake risk cost, wherein the risk control measures include: strengthening structures, sealing repairs, and increasing emergency reserves; The method library forms a sub-module, which is used to create risk control measures for different earthquake risk costs, forming a disposal method library under different intensities of earthquake motion conditions.

[0017] A third aspect of an embodiment of the present invention provides an electronic device, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the computer program is executed by the processor, an earthquake monitoring method based on gas storage reservoir structural parameters as in the first aspect of an embodiment of the present invention is implemented.

[0018] A fourth aspect of an embodiment of the present invention provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the seismic monitoring method based on gas storage reservoir structural parameters of the first aspect of the embodiment of the present invention is implemented.

[0019] Through the seismic monitoring method based on gas storage reservoir structural parameters of the embodiment of the present invention, seismic geological data, engineering geological and structural design information of the gas storage reservoir construction area are collected, the seismic geological data include: historical regional seismic records and seismic motion parameters of the gas storage reservoir construction area, and the engineering geological and structural design information include: regional rock and soil layer distribution characteristics, groundwater conditions, ground stress field, rock and soil body physical and mechanical parameters, gas storage reservoir body design parameters, gas storage reservoir structural parameters; based on the historical regional seismic records and seismic motion parameters of the gas storage reservoir construction area, multiple groups of random seismic motion time history records with different exceedance probabilities are generated; based on the seismic geological data and the engineering geological and structural design information, a three-dimensional numerical simulation model is established; using the generated multiple groups of random seismic motion time history records with different exceedance probabilities as input conditions, dynamic response calculation of the three-dimensional numerical simulation model is performed to obtain the damage The damage degree information includes: the damage range of the sealing layer and the damage range of the lining; using a risk quantitative evaluation model, according to the exceedance probability in the input conditions and the damage degree information under the earthquake conditions, the earthquake risk cost under the input conditions is obtained; based on the earthquake risk cost, corresponding risk control measures are created to obtain a disposal method library under earthquake motion conditions of different intensities; according to multiple groups of random earthquake motion time history records under different exceedance probabilities and the corresponding earthquake risk costs, an earthquake risk assessment model is constructed, and risk control measures matching different earthquake risk costs are determined in the disposal method library; during the operation of the current gas storage, the earthquake motion in the area where the current gas storage is located and the stress level of the gas storage structure are monitored; when an earthquake occurs, the earthquake risk cost of the current gas storage is quickly assessed through the earthquake risk assessment model, and the corresponding risk control measures are output.

[0020] In this embodiment, based on the historical regional earthquake records and seismic parameters of the gas storage construction area, multiple groups of random seismic time history records under different exceedance probabilities are generated, and then based on the seismic geological data and the engineering geological and structural design information, a three-dimensional numerical simulation model is established, and then the multiple groups of random seismic time history records under different exceedance probabilities are used as input conditions to perform dynamic response calculations of the three-dimensional numerical simulation model to obtain damage degree information, and then obtain the seismic risk cost under the input conditions, so as to achieve the purpose of quantitatively evaluating the risk cost brought by the earthquake; finally, based on the multiple groups of random seismic time history records under different exceedance probabilities and the corresponding seismic risk costs, an earthquake risk assessment model is constructed, and a disposal method library corresponding to seismic conditions of different intensities is used to assess the seismic risk cost of the current gas storage, and output corresponding risk control measures, so as to achieve the purpose of assessing the risk cost in the first time when an earthquake occurs and making corresponding risk management. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0022] Figure 1 is a flow chart of a method for earthquake monitoring based on gas storage structure parameters according to an embodiment of the present invention; Figure 2 It is a structural block diagram of an earthquake monitoring device based on gas storage reservoir structural parameters provided by an embodiment of the present invention; Figure 3 It is a schematic diagram of an electronic device shown in an embodiment of the present invention. DETAILED DESCRIPTION

[0023] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0024] Please refer to Figure 1 , Figure 1 FIG. 1 is a flow chart of an earthquake monitoring method based on gas storage reservoir structural parameters according to an embodiment of the present invention. Figure 1 As shown, the method may include steps S101 to S109: Step S101: Collecting seismic geological data, engineering geological and structural design information of the gas storage construction area, wherein the seismic geological data include: historical regional earthquake records and seismic motion parameters of the gas storage construction area, and the engineering geological and structural design information include: regional rock and soil layer distribution characteristics, groundwater conditions, geostress field, rock and soil physical and mechanical parameters, gas storage body design parameters, and gas storage structure parameters.

[0025] In this embodiment, since earthquakes are low-probability events, a project may not encounter an earthquake during its 50-year operation period. Therefore, it is necessary to collect historical earthquake geological data and combine them with engineering geology and structural design information to simulate the possible impacts on the gas storage facility under different earthquake conditions.

[0026] For example, seismic geology data, engineering geology, and structural design information can be obtained through open source database downloads, drone measurements, and purchases.

[0027] Step S102: Based on the historical regional earthquake records and seismic parameters of the gas storage construction area, multiple groups of random seismic time history records with different exceedance probabilities are generated.

[0028] In this embodiment, multiple groups of random seismic time history records with different exceedance probabilities can be obtained by using, but not limited to, the Monte Carlo simulation method. The exceedance probability is the probability that the gas storage facility may encounter an earthquake intensity value greater than or equal to a given seismic intensity value or seismic parameter value within a certain period of time. The seismic time history describes the process of changes in seismic acceleration, velocity or displacement over time, and includes three key elements, namely, the intensity, spectrum and duration of the seismic motion. The intensity is usually characterized by the maximum amplitude (such as the peak acceleration); the spectral characteristics reflect the intensity distribution of different frequency components; and the duration affects the final degree of damage to the structure, especially the structural response in the nonlinear stage.

[0029] By using multiple sets of random seismic time history records with different random exceedance probabilities, the impact of earthquakes of different magnitudes when they actually occur can be simulated.

[0030] Step S103: establishing a three-dimensional numerical simulation model based on the seismic geological data and the engineering geological and structural design information.

[0031] In this embodiment, the three-dimensional numerical simulation model of the gas storage reservoir can be established by, but not limited to, using a finite element structural model establishment method.

[0032] In the present application, a simulation model can be established for an existing gas storage facility, a risk assessment of the gas storage facility can be conducted, and suggestions can be given for its corresponding earthquake risk management. At the same time, a simulation model can also be established for a gas storage facility that has not yet been built based on its design data, and a risk assessment of the gas storage facility can be conducted to estimate the earthquake risk that it may subsequently suffer.

[0033] Step S104: using the generated multiple groups of random seismic time history records under different exceedance probabilities as input conditions, performing dynamic response calculation of the three-dimensional numerical simulation model to obtain damage degree information, wherein the damage degree information includes: the damage range of the sealing layer and the damage range of the lining.

[0034] In this embodiment, multiple groups of random seismic time history records with different exceedance probabilities are generated as input conditions to simulate earthquake conditions of different degrees. The process of calculating the dynamic response of the three-dimensional numerical simulation model can be carried out by but not limited to the method of using dynamic analysis model calculation to finally obtain the specific impact of different earthquake conditions on the gas storage, that is, the damage degree information.

[0035] In this application, since the gas storage has higher requirements on structural stability and airtightness, the damage degree information of the gas storage is determined by the damage range of the sealing layer and the lining.

[0036] Step S105: using a risk quantitative assessment model, according to the exceedance probability in the input conditions and the damage degree information under the earthquake conditions, to obtain the earthquake risk cost under the input conditions.

[0037] In this embodiment, a monetary risk quantitative evaluation model is used to evaluate the impact of earthquakes on gas storage operations based on the probability of earthquake occurrence and the degree of damage under the earthquake conditions. The earthquake risk cost includes the repair cost, that is, the earthquake risk cost under a certain earthquake intensity, and the production loss, that is, the time cost impact of repairing the losses caused by the earthquake.

[0038] Step S106: Based on the earthquake risk cost, corresponding risk control measures are created to obtain a library of disposal methods under earthquake motion conditions of different intensities.

[0039] In this embodiment, based on the earthquake risk cost, the repair measures to be taken for the gas storage when it is affected by earthquakes of different degrees are integrated and classified. The collection of different response and treatment measures is the disposal method library.

[0040] Different methods should be adopted according to the different degrees of damage suffered by the gas storage facility. For example, if the damage is minor, only filling and repair can be considered. If there is extensive damage, structural reinforcement or overall replacement will need to be considered.

[0041] Step S107: construct an earthquake risk assessment model based on multiple groups of random earthquake motion time history records under different exceedance probabilities and corresponding earthquake risk costs, and determine risk control measures plans that match different earthquake risk costs in the disposal method library.

[0042] In this embodiment, multiple sets of earthquake risk cost calculation results can be integrated through but not limited to the use of machine learning methods to construct an earthquake risk assessment model for predicting the disasters that may be caused by subsequent earthquakes, and matching it with a disposal method library to obtain risk control measures that match different earthquake risk costs.

[0043] Step S108: During the operation of the current gas storage, the seismic motion in the area where the current gas storage is located and the stress level of the gas storage structure are monitored.

[0044] In this embodiment, during the actual operation of the gas storage facility, real-time seismic monitoring should be carried out in the surrounding areas to provide timely warning and response to possible earthquake disasters. At the same time, the structural stress level of the gas storage facility itself should be monitored. Since the structural stability of the gas storage facility can be reflected by the structural stress level, monitoring the structural stress level of the gas storage facility itself can ensure that problems are discovered and handled in a timely manner, which is beneficial to the maintenance and use of the gas storage facility.

[0045] Step S109: Assess the earthquake risk cost of the current gas storage reservoir through an earthquake risk assessment model, and output corresponding risk control measures.

[0046] In this embodiment, after the earthquake risk assessment model is established, the earthquake risk cost of the gas storage facility can be judged at any time according to the seismic conditions in the current area. Therefore, when an actual earthquake occurs, the earthquake risk can be timely assessed through the earthquake assessment model associated with the disposal method library, and corresponding response plans can be made.

[0047] In combination with the above embodiments, in one implementation, the present invention further provides a seismic monitoring method based on gas storage reservoir structural parameters. In the method, the above step S104 may specifically include steps S201 to S202: Step S201: Based on the generated multiple groups of random seismic time history records under different exceedance probabilities as input conditions, dynamic response calculation of the three-dimensional numerical simulation model is performed, and the stress of the sealing layer and the lining under each group of the input conditions is recorded.

[0048] In this embodiment, the damage degree information includes at least: the damage range of the sealing layer and the damage range of the lining. Among them, the generated multiple groups of random seismic time history records under different exceedance probabilities are used as input conditions to simulate earthquake conditions of different degrees. The process of dynamic response calculation of the three-dimensional numerical simulation model can be achieved by, but not limited to, using the method of dynamic analysis model calculation, and finally obtain the specific impact of different earthquake conditions on the gas storage. The specific manifestation of the damage to the sealing layer and the lining is the stress received by the sealing layer and the lining.

[0049] Step S202: The allowable stress of the sealing layer material is σ Rs , the allowable stress of the lining material is σ Rc In each group of calculation results, the sealing layer stress is greater than the allowable stress σ Rs The total volume of the region V s , is determined as the damage range of the sealing layer; the lining stress is greater than σ Rc The total volume of the region V c , determined as the damage range of the lining.

[0050] In this embodiment, the damage of the sealing layer and lining can be judged by calculating the stress borne by the sealing layer and lining. The allowable stress threshold of the sealing layer and lining of the gas storage reservoir can be obtained through the gas storage reservoir design documents and actual construction information. When the sealing layer or lining is subjected to stress exceeding the threshold, it can be regarded as damaged. By calculating the total volume of the area subjected to stress exceeding the allowable stress, the damage range of the sealing layer or lining can be obtained.

[0051] In combination with the above embodiments, in one implementation, the present invention further provides a seismic monitoring method based on gas storage reservoir structural parameters. In the method, the above step S105 may specifically include steps S301 to S303: Step S301: Determine the repair cost C according to the damage range of the sealing layer and the damage range of the lining. repair , the repair cost C repair The amount of work used to reflect the repair work.

[0052] In this embodiment, the damage range of the sealing layer and the damage range of the lining are used to determine the repair cost C repair According to the extent of the damage, the amount of work required for the repair project will also vary, and the economic cost of the repair project will also vary accordingly. The repair of the gas storage facility includes the repair of the sealing layer and the repair of the lining.

[0053] Step S302: Determine the production reduction loss C according to the number of days of shutdown and the daily benefit of the energy storage power station. loss The number of days of suspension of production is related to the workload of the repair project.

[0054] In this embodiment, the earthquake risk cost also includes the production loss C loss After suffering an earthquake disaster, the energy storage power station supplied by the gas storage reservoir cannot continue to work during the repair process. Therefore, the earthquake risk cost includes not only the repair cost of the gas storage reservoir, but also the economic loss caused by the suspension of the energy storage power station. The number of days of suspension depends on the degree of damage to the gas storage reservoir. The greater the degree of damage, the greater the amount of repair work and the more days of suspension. The smaller the degree of damage, the smaller the amount of repair work and the fewer days of suspension. The loss of production reduction C loss It is also related to the daily benefits of the energy storage power station. The product of the daily economic benefits of the energy storage power station and the number of days of shutdown is the production loss C. loss .

[0055] Step S303: P f ×(C repair +C loss ) is determined as the earthquake risk cost, where P f is the probability of occurrence of earthquake intensity at the site.

[0056] In this example, the earthquake risk cost is not only affected by the repair cost C repair and production loss C loss The two aspects of influence also depend on the probability of occurrence of earthquake intensity P at the gas storage site. f The probability of occurrence includes the probability of earthquake occurrence and the probability that the earthquake intensity will cause damage to the gas storage. Taking the three factors into consideration, the earthquake risk cost P of the earthquake risk assessment model is constructed. f ×(C repair +C loss ).

[0057] In combination with the above embodiments, in one implementation, the present invention further provides a seismic monitoring method based on gas storage reservoir structural parameters. In the method, the above step S107 may specifically include step S401: Step S401: using the multiple groups of random seismic time history records under different exceedance probabilities and the corresponding earthquake risk costs as training samples, training the machine learning model to obtain the earthquake risk assessment model.

[0058] In this embodiment, first, for each sample seismic time history record, the sample seismic time history record is marked according to whether the sample seismic time history record indicates that the earthquake will cause damage to the three-dimensional numerical simulation model. The marking process can be completed manually, that is, manually judging whether the sample seismic time history record poses a risk to the three-dimensional numerical simulation model, and then marking the sample seismic time history record, so that the sample seismic time history record carries the mark.

[0059] Each sample seismic time history record is marked according to the above process, thereby obtaining a plurality of sample seismic time history records carrying marks.

[0060] Then, features are extracted from each labeled sample seismic time history record, and the extracted features are input into a classifier (e.g., softmax, Bayesian, SVM, xgboost) for training until the first preset model after training has the function of predicting the probability of a three-dimensional numerical simulation model suffering from earthquake risk represented by a single seismic time history record. The first preset model after training is the earthquake risk assessment model.

[0061] Finally, after obtaining the earthquake risk assessment model, for a single seismic motion time history record, for example, the seismic motion time history record of an actual earthquake, the single seismic motion time history record can be input into the earthquake risk assessment model, and the earthquake risk assessment model can be used to predict the probability that the gas storage facility may be subject to earthquake risk represented by the seismic motion time history record.

[0062] In combination with the above embodiments, in one implementation, the present invention further provides a seismic monitoring method based on gas storage reservoir structural parameters. In the method, the above step S102 may specifically include steps S501 to S502: Step S501: Determine the upper and lower limits of the ground motion intensity PGA based on the historical regional earthquake records in the gas storage construction area.

[0063] In this example, due to the regional nature of earthquake risk assessment, for the northwestern region where earthquakes are more frequent, establishing upper and lower limits of earthquake intensity can more specifically simulate possible earthquake situations, and the random seismic motion time history records used as input are more realistic.

[0064] Step S502: For the target PGA, based on the Monte Carlo simulation method, with PGA as the target value, generate the multiple groups of random seismic time history records under different exceedance probabilities.

[0065] In this example, multiple groups of random seismic time history records with different exceedance probabilities can be obtained by, but not limited to, using the Monte Carlo simulation method. The Monte Carlo simulation method establishes a probability model or a random process so that its parameters or numerical characteristics are equal to the solution of the problem: these parameters or numerical characteristics are then calculated through observation or sampling tests of the model or process, and finally the approximate value of the solution is given. The error of this method is independent of the dimension of the problem. It can directly solve problems with statistical properties, and there is no need to discretize continuous problems. In the application of random seismic time history record simulation, the randomness of earthquake occurrence can be simulated more accurately, and a large amount of data input can be provided for the establishment of a reliable earthquake risk assessment model.

[0066] In combination with the above embodiments, in one implementation, the present invention further provides a seismic monitoring method based on gas storage reservoir structural parameters. In the method, the above step S106 may specifically include steps S601 to S602: Step S601: Based on the earthquake risk cost, risk control measures for corresponding engineering quantities are created, wherein the risk control measures include: strengthening structures, sealing and repairing, and increasing emergency reserves.

[0067] In this example, based on the earthquake risk cost, targeted risk control measures are proposed, such as strengthening structures, sealing and repairing, and increasing emergency reserves (the engineering volume of the emergency measures plan is totaled according to the earthquake risk cost). According to the specific quantified risk cost obtained by the risk assessment model, corresponding emergency reserves can be made to quickly respond to earthquake damage when an earthquake disaster occurs.

[0068] Step S602: Risk control measures are created for different earthquake risk costs to form a disposal method library under earthquake motion conditions of different intensities.

[0069] In this example, due to the randomness of earthquake occurrence, different degrees of response measures should be taken for different degrees of earthquake damage. Therefore, the risk assessment model can be used to estimate the loss of gas storage caused by the earthquake, evaluate the amount of post-disaster repair work, and form a disposal plan based on the amount of work. A library of disposal methods under different intensities of seismic motion conditions can be formed. After an extreme disaster occurs, the speed of risk response can be improved to provide support for the risk management plan.

[0070] It should be noted that, for the sake of simplicity, the method embodiments are described as a series of action combinations, but those skilled in the art should be aware that the embodiments of the present invention are not limited by the order of the actions described, because according to the embodiments of the present invention, certain steps can be performed in other orders or simultaneously. Secondly, those skilled in the art should also be aware that the embodiments described in the specification are all preferred embodiments, and the actions involved are not necessarily required by the embodiments of the present invention.

[0071] Based on the same inventive concept, an embodiment of the present invention provides an earthquake monitoring device 700 based on gas storage reservoir structural parameters. Figure 2 , Figure 2 FIG. 1 is a structural block diagram of an earthquake monitoring device based on gas storage reservoir structural parameters provided by an embodiment of the present invention. Figure 2 As shown, the device 700 includes: The data collection module 701 is used to collect seismic geological data, engineering geological and structural design information of the gas storage construction area, wherein the seismic geological data includes: historical regional earthquake records and seismic motion parameters of the gas storage construction area, and the engineering geological and structural design information includes: regional rock and soil layer distribution characteristics, groundwater conditions, ground stress field, rock and soil body physical and mechanical parameters, gas storage body design parameters, and gas storage structure parameters; An input condition generation module 702 is used to generate multiple groups of random seismic time history records under different exceedance probabilities based on historical regional earthquake records and seismic parameters in the gas storage construction area; A simulation model building module 703 is used to build a three-dimensional numerical simulation model based on the seismic geological data and the engineering geological and structural design information; The damage degree determination module 704 is used to use the generated multiple groups of random earthquake time history records under different exceedance probabilities as input conditions to perform dynamic response calculation of the three-dimensional numerical simulation model to obtain damage degree information, wherein the damage degree information includes: the damage range of the sealing layer and the damage range of the lining; The risk quantitative evaluation module 705 is used to obtain the earthquake risk cost under the input conditions according to the exceedance probability in the input conditions and the damage degree information under the earthquake conditions; The disposal method library establishment module 706 is used to create corresponding risk control measures based on earthquake risk costs, and obtain a disposal method library under different intensities of earthquake motion conditions; The earthquake risk assessment module 707 is used to construct an earthquake risk assessment model based on multiple groups of random earthquake motion time history records under different exceedance probabilities and corresponding earthquake risk costs, and determine risk control measures plans that match different earthquake risk costs in the disposal method library; A monitoring module 708 is used to monitor the earthquake motion in the area where the current gas storage is located and the stress level of the gas storage structure during the operation of the current gas storage; The risk control measures output module 709 is used to evaluate the earthquake risk cost of the current gas storage reservoir through an earthquake risk assessment model and output corresponding risk control measures.

[0072] Optionally, the generated multiple groups of random earthquake time history records under different exceedance probabilities are used as input conditions to perform dynamic response calculation of the three-dimensional numerical simulation model to obtain damage degree information; The damage degree determination module 704 includes: A dynamic response calculation submodule is used to perform dynamic response calculation of a three-dimensional numerical simulation model based on the generated multiple groups of random seismic time history records under different exceedance probabilities as input conditions, and record the stress of the sealing layer and the lining under each group of input conditions; The damage range determination submodule is used to determine the damage range of the sealing layer stress greater than the allowable stress σ in each group of calculation results. Rs The total volume of the region V s , is determined as the damage range of the sealing layer; the lining stress is greater than σ Rc The total volume of the region V c , determined as the damage range of the lining.

[0073] Optionally, the earthquake risk cost includes: repair cost and production loss; a risk quantitative evaluation model is used to obtain the earthquake risk cost under the input conditions according to the exceedance probability in the input conditions and the damage degree information under the earthquake conditions; The risk quantitative assessment module 705 includes: The repair cost determination submodule is used to determine the repair cost C according to the damage range of the sealing layer and the damage range of the lining. repair , the repair cost C repair The amount of work used for reaction repair works; The production reduction loss determination submodule is used to determine the production reduction loss C according to the number of days of shutdown and the daily benefit of the energy storage power station. loss , the number of days of suspension of production is related to the amount of work of the repair project; The risk cost determination submodule converts P f ×(C repair +C loss ) is determined as the earthquake risk cost, where P f is the probability of occurrence of earthquake intensity at the site.

[0074] Optionally, an earthquake risk assessment model is constructed according to the multiple groups of random earthquake motion time history records under different exceedance probabilities and corresponding earthquake risk costs; The earthquake risk assessment module 707 includes: The assessment model training submodule is used to use the multiple groups of random seismic time history records under different exceedance probabilities and the corresponding earthquake risk costs as training samples to train the machine learning model to obtain the earthquake risk assessment model.

[0075] Optionally, based on historical regional earthquake records and seismic parameters in the gas storage construction area, multiple groups of random seismic time history records with different exceedance probabilities are generated; Optionally, the input condition generating module 702 includes: A historical record confirmation submodule, used to determine the upper and lower limits of the ground motion intensity PGA based on the historical regional earthquake records in the gas storage construction area; The model generation submodule is used to generate the multiple groups of random seismic time history records under different exceedance probabilities for the target PGA based on the Monte Carlo simulation method and taking PGA as the target value.

[0076] Optionally, based on the earthquake risk cost, corresponding risk control measures are created to obtain a library of disposal methods under earthquake motion conditions of different intensities; The disposal method library establishment module 706 includes: A risk control measure determination submodule is used to create risk control measures for corresponding engineering quantities based on the earthquake risk cost, wherein the risk control measures include: strengthening structures, sealing repairs, and increasing emergency reserves; The method library forms a sub-module, which is used to create risk control measures for different earthquake risk costs, forming a disposal method library under different intensities of earthquake motion conditions.

[0077] Based on the same inventive concept, another embodiment of the present invention provides an electronic device 800, such as Figure 3 shown. Figure 3801 and 802. The electronic device includes a memory 802, a processor 801 and a computer program stored in the memory and executable on the processor, and the processor implements the steps of the seismic monitoring method based on gas storage reservoir structural parameters described in any of the above embodiments of the present invention when executing the computer program.

[0078] As for the device embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment.

[0079] The various embodiments in this specification are described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the various embodiments can be referenced to each other.

[0080] Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, devices, or computer program products. Therefore, embodiments of the present invention may take the form of complete hardware embodiments, complete software embodiments, or embodiments combining software and hardware. Furthermore, embodiments of the present invention may take the form of computer program products implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes.

[0081] The embodiments of the present invention are described with reference to the flowcharts and / or block diagrams of the methods, terminal devices (systems), and computer program products according to the embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of the processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing terminal device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing terminal device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0082] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing terminal device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce a manufactured product including an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0083] These computer program instructions can also be loaded onto a computer or other programmable data processing terminal device so that a series of operating steps are executed on the computer or other programmable terminal device to produce a computer-implemented process, thereby providing instructions for executing on the computer or other programmable terminal device to implement the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0084] Although the preferred embodiments of the present invention have been described, those skilled in the art may make additional changes and modifications to these embodiments once they have learned the basic creative concept. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments and all changes and modifications that fall within the scope of the embodiments of the present invention.

[0085] Finally, it should be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or terminal device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or terminal device. In the absence of further restrictions, the elements defined by the sentence "comprise a ..." do not exclude the existence of other identical elements in the process, method, article or terminal device including the elements.

[0086] The above is a detailed introduction to the seismic monitoring method, device, equipment and medium based on the structural parameters of the gas storage reservoir provided by the present invention. Specific examples are used in this article to illustrate the principles and implementation methods of the present invention. The description of the above embodiments is only used to help understand the method of the present invention and its core idea; at the same time, for those skilled in the art, according to the idea of ​​the present invention, there will be changes in the specific implementation method and application scope. In summary, the content of this specification should not be understood as a limitation on the present invention.

Claims

1. A seismic monitoring method based on gas storage reservoir structural parameters, characterized in that: The method comprises: Collecting seismic geological data, engineering geological and structural design information of the gas storage construction area, wherein the seismic geological data include: historical regional earthquake records and seismic motion parameters of the gas storage construction area, and the engineering geological and structural design information include: regional rock and soil layer distribution characteristics, groundwater conditions, ground stress field, rock and soil physical and mechanical parameters, gas storage body design parameters, and gas storage structure parameters; Based on the historical regional earthquake records and seismic parameters of the gas storage construction area, multiple groups of random seismic time history records with different exceedance probabilities are generated; Based on the seismic geological data and the engineering geological and structural design information, a three-dimensional numerical simulation model is established; The generated multiple groups of random earthquake time history records under different exceedance probabilities are used as input conditions to calculate the dynamic response of the three-dimensional numerical simulation model to obtain damage degree information, which includes: the damage range of the sealing layer and the damage range of the lining; Using a risk quantitative assessment model, according to the exceedance probability in the input conditions and the damage degree information under the earthquake conditions, the earthquake risk cost under the input conditions is obtained; Based on the earthquake risk cost, corresponding risk control measures are created to obtain a library of disposal methods under different intensities of earthquake motion conditions; Based on multiple sets of random seismic time history records with different exceedance probabilities and the corresponding earthquake risk costs, an earthquake risk assessment model is constructed, and risk control measures matching different earthquake risk costs are determined in the disposal method library; During the operation of the current gas storage, monitoring the earthquake motion in the area where the current gas storage is located and the stress level of the gas storage structure; The earthquake risk cost of the current gas storage facility is evaluated through an earthquake risk assessment model, and corresponding risk control measures are output.

2. The earthquake monitoring method based on gas storage reservoir structural parameters according to claim 1 is characterized in that: The method uses the generated multiple groups of random earthquake time history records under different exceedance probabilities as input conditions to perform dynamic response calculation of the three-dimensional numerical simulation model to obtain damage degree information, including: Based on the generated multiple groups of random earthquake time history records under different exceedance probabilities as input conditions, dynamic response calculation of the three-dimensional numerical simulation model is performed, and the stress of the sealing layer and the lining under each group of the input conditions is recorded; Among them, the allowable stress of the sealing layer material is σ Rs , the allowable stress of the lining material is σ Rc In each group of calculation results, the sealing layer stress is greater than the allowable stress σ Rs The total volume of the region V s , is determined as the damage range of the sealing layer; the lining stress is greater than σ Rc The total volume of the region V c , determined as the damage range of the lining.

3. The earthquake monitoring method based on gas storage reservoir structural parameters according to claim 1 is characterized in that: The earthquake risk cost includes: repair cost and production loss; the risk quantitative evaluation model is used to obtain the earthquake risk cost under the input conditions according to the exceedance probability in the input conditions and the damage degree information under the earthquake conditions, including: Determine the repair cost C according to the damage range of the sealing layer and the damage range of the lining. repair , the repair cost C repair The amount of work used for reaction repair works; According to the number of days of shutdown and the daily benefit of the energy storage power station, the loss of production reduction C is determined. loss , the number of days of suspension of production is related to the amount of work of the repair project; P f ×(C repair +C loss ) is determined as the earthquake risk cost, where P f is the probability of occurrence of earthquake intensity at the site.

4. The earthquake monitoring method based on gas storage reservoir structural parameters according to claim 1 is characterized in that: According to the multiple groups of random earthquake motion time history records under different exceedance probabilities and the corresponding earthquake risk costs, an earthquake risk assessment model is constructed, including: The multiple groups of random seismic time history records under different exceedance probabilities and the corresponding earthquake risk costs are used as training samples to train the machine learning model to obtain the earthquake risk assessment model.

5. The earthquake monitoring method based on gas storage reservoir structural parameters according to claim 1 is characterized in that: Based on the historical regional earthquake records and seismic parameters of the gas storage construction area, multiple groups of random seismic time history records with different exceedance probabilities are generated, including: Based on the historical regional earthquake records of the gas storage construction area, determine the upper and lower limits of the ground motion intensity PGA; For the target PGA, based on the Monte Carlo simulation method, the multiple groups of random seismic time history records under different exceedance probabilities are generated with PGA as the target value.

6. The earthquake monitoring method based on gas storage reservoir structural parameters according to claim 3 is characterized in that: Based on the earthquake risk cost, corresponding risk control measures are created to obtain a library of disposal methods under different intensities of earthquake motion conditions, including: Based on the earthquake risk cost, risk control measures for corresponding engineering quantities are created, and the risk control measures include: strengthening structures, sealing and repairing, and increasing emergency reserves; Risk control measures created for different earthquake risk costs form a library of disposal methods under earthquake motion conditions of different intensities.

7. An earthquake monitoring device based on gas storage reservoir structural parameters, characterized in that: The device comprises: A data collection module is used to collect seismic geological data, engineering geological and structural design information of the gas storage construction area, wherein the seismic geological data include: historical regional earthquake records and seismic motion parameters of the gas storage construction area, and the engineering geological and structural design information include: regional rock and soil layer distribution characteristics, groundwater conditions, ground stress field, rock and soil body physical and mechanical parameters, gas storage body design parameters, and gas storage structure parameters; An input condition generation module is used to generate multiple groups of random seismic time history records under different exceedance probabilities based on historical regional earthquake records and seismic parameters in the gas storage construction area; A simulation model building module, used to build a three-dimensional numerical simulation model based on the seismic geological data and the engineering geological and structural design information; A damage degree determination module is used to use the generated multiple groups of random seismic time history records under different exceedance probabilities as input conditions to perform dynamic response calculation of the three-dimensional numerical simulation model to obtain damage degree information, wherein the damage degree information includes: the damage range of the sealing layer and the damage range of the lining; A risk quantitative evaluation module, used to obtain the earthquake risk cost under the input conditions according to the exceedance probability in the input conditions and the damage degree information under the earthquake conditions; The disposal method library establishment module is used to create corresponding risk control measures based on earthquake risk costs and obtain a disposal method library under different intensities of earthquake motion conditions; The earthquake risk assessment module is used to construct an earthquake risk assessment model based on multiple sets of random earthquake motion time history records with different exceedance probabilities and the corresponding earthquake risk costs, and to determine the risk control measures matching different earthquake risk costs in the disposal method library; A monitoring module, used to monitor the seismic motion in the area where the current gas storage is located and the stress level of the gas storage structure during the operation of the current gas storage; The risk control measures output module is used to evaluate the earthquake risk cost of the current gas storage reservoir through an earthquake risk assessment model and output corresponding risk control measures.

8. The earthquake monitoring device based on gas storage reservoir structural parameters according to claim 7, characterized in that: The damage degree determination module comprises: A dynamic response calculation submodule is used to perform dynamic response calculation of a three-dimensional numerical simulation model based on the generated multiple groups of random seismic time history records under different exceedance probabilities as input conditions, and record the stress of the sealing layer and the lining under each group of the input conditions; The damage range determination submodule is used to determine the sealing layer stress greater than the allowable stress σ of the sealing layer material in each group of calculation results. Rs The total volume of the region V s , is determined as the damage range of the sealing layer; the lining stress is greater than the allowable stress σ of the lining material Rc The total volume of the region V c , determined as the damage range of the lining.

9. The earthquake monitoring device based on gas storage reservoir structural parameters according to claim 7, characterized in that: The risk quantitative assessment module includes: The repair cost determination submodule is used to determine the repair cost C according to the damage range of the sealing layer and the damage range of the lining. repair , the repair cost C repair The amount of work used for reaction repair works; The production reduction loss determination submodule is used to determine the production reduction loss C according to the number of days of shutdown and the daily benefit of the energy storage power station. loss , the number of days of suspension of production is related to the amount of work of the repair project; The risk cost determination submodule converts P f ×(C repair +C loss ) is determined as the earthquake risk cost, where P f is the probability of occurrence of earthquake intensity at the site.

10. The earthquake monitoring device based on gas storage reservoir structural parameters according to claim 7, characterized in that: The earthquake risk assessment module comprises: The assessment model training submodule is used to use the multiple groups of random seismic time history records under different exceedance probabilities and the corresponding earthquake risk costs as training samples to train the machine learning model to obtain the earthquake risk assessment model.

11. The earthquake monitoring device based on gas storage reservoir structural parameters according to claim 7, characterized in that: The input condition generation module comprises: A historical record confirmation submodule, used to determine the upper and lower limits of the ground motion intensity PGA based on the historical regional earthquake records in the gas storage construction area; The model generation submodule is used to generate the multiple groups of random seismic time history records under different exceedance probabilities for the target PGA based on the Monte Carlo simulation method and taking PGA as the target value.

12. The earthquake monitoring device based on gas storage reservoir structural parameters according to claim 9, characterized in that: The disposal method library establishment module includes: A risk control measure determination submodule is used to create risk control measures for corresponding engineering quantities based on the earthquake risk cost, wherein the risk control measures include: strengthening structures, sealing repairs, and increasing emergency reserves; The method library forms a sub-module, which is used to create risk control measures for different earthquake risk costs, forming a disposal method library under different intensities of earthquake motion conditions.

13. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the computer program is executed by the processor, the seismic monitoring method based on gas storage reservoir structural parameters as described in any one of claims 1 to 6 is implemented.

14. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the seismic monitoring method based on gas storage reservoir structural parameters as described in any one of claims 1 to 6 is implemented.

15. A computer program product comprising a computer program, characterized in that When the computer program is executed by a processor, the seismic monitoring method based on gas storage reservoir structural parameters as described in any one of claims 1 to 6 is implemented.

Images