CN113253343B - Method for identifying fault activity of underground gas storage based on microseism monitoring technology - Google Patents
Method for identifying fault activity of underground gas storage based on microseism monitoring technology Download PDFInfo
- Publication number
- CN113253343B CN113253343B CN202110518177.0A CN202110518177A CN113253343B CN 113253343 B CN113253343 B CN 113253343B CN 202110518177 A CN202110518177 A CN 202110518177A CN 113253343 B CN113253343 B CN 113253343B
- Authority
- CN
- China
- Prior art keywords
- microseism
- events
- value
- gas storage
- cluster
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000000694 effects Effects 0.000 title claims abstract description 39
- 238000003860 storage Methods 0.000 title claims abstract description 34
- 238000012544 monitoring process Methods 0.000 title claims abstract description 24
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000005516 engineering process Methods 0.000 title claims abstract description 10
- 238000004220 aggregation Methods 0.000 claims abstract description 6
- 230000002776 aggregation Effects 0.000 claims abstract description 6
- 238000007789 sealing Methods 0.000 abstract description 7
- 238000002347 injection Methods 0.000 abstract description 4
- 239000007924 injection Substances 0.000 abstract description 4
- 238000004519 manufacturing process Methods 0.000 abstract description 2
- 239000007789 gas Substances 0.000 description 29
- 238000004364 calculation method Methods 0.000 description 6
- 238000009826 distribution Methods 0.000 description 3
- 239000000700 radioactive tracer Substances 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/288—Event detection in seismic signals, e.g. microseismics
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/642—Faults
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a method for identifying fault activity of an underground gas storage based on a microseism monitoring technology, which comprises the following steps: 1) inputting a microseism positioning result and attribute information; 2) dividing the microseism events into R clusters according to the aggregation degree of the microseism events in the step 1); 3) calculating the b value of each cluster of microseism events according to the clustering condition in the step 2); 4) calculating the D value of each cluster of microseism events according to the clustering condition in the step 2); 5) and calculating the fault activity index F according to the result of the b value in the step 3) and the result of the D value in the step 4). The method is based on the microseism monitoring result, calculates the fault activity index of different cluster microseism events by using the attribute information of the microseism events, effectively identifies the fault activity in the injection and production operation process of the underground gas storage, can effectively prevent the fault activity sealing failure risk of the gas storage, and ensures the safe operation of the gas storage.
Description
Technical Field
The invention relates to a method for processing and explaining microseism data in gas storage safety monitoring and geophysical exploration, in particular to a method for identifying fault activity of an underground gas storage based on a microseism monitoring technology.
Background
The underground gas storage has no replaceable function and obvious advantages in the aspects of energy strategy and winter peak regulation to ensure the safe gas utilization of the civil life. The complex fault development of the geological structure of the domestic gas storage, and the gas storage has the characteristics of multi-cycle forced injection and forced extraction, alternating stress change and the like, and fault activity can be induced in the long-term injection and extraction process to cause the fault sealing failure of the gas storage and natural gas leakage accidents.
The fault tightness of the gas storage is a key factor for safe operation of the gas storage. Therefore, the gas storage construction strength is accelerated, and simultaneously, a gas storage fault monitoring system is continuously perfected, so that the safe construction and operation of the gas storage are guaranteed. At present, the safety monitoring of the fault of the gas storage is mainly carried out by two monitoring modes: firstly, monitoring is carried out by utilizing parameters such as temperature, pressure and fluid components, monitoring wells are arranged on two sides of a fault to monitor the parameters such as the temperature, the pressure and the fluid components, and when the parameters on the two sides of the fault are more and more the same, the fault sealing failure risk is considered to exist. Secondly, monitoring by using a tracer, injecting the tracer into the gas storage on one side of the fault, and if the same tracer is identified in the monitoring well on the other side of the fault outside the gas storage, determining that the fault sealing performance is invalid. The two monitoring modes have problems, namely the monitoring range is small, the fault tightness around the monitoring well can be monitored only, and when the monitored parameters reflect the problems, the fault tightness fails.
Disclosure of Invention
The invention aims to identify fault activities of an underground gas storage by using a microseism monitoring technology, avoid failure of fault tightness in the operation process of the gas storage, ensure safe operation of the gas storage and provide a method for identifying fault activities by using the microseism monitoring technology. The method utilizes the microseism monitoring technology to monitor the activity of the underground city of the gas storage in real time, accurately identifies the fault activity according to the distribution position and the related attributes of the microseism monitoring event positioning result, and can timely take related measures according to the fault activity condition in the later period to avoid fault sealing failure and ensure the safe operation of the gas storage.
The specific technical scheme is as follows:
the method for identifying the fault activity of the underground gas storage based on the microseism monitoring technology comprises the following steps:
1) inputting a microseism positioning result and attribute information;
2) dividing the microseism events into R clusters according to the aggregation degree of the microseism events in the step 1);
3) calculating the b value of each cluster of microseism events according to the clustering condition in the step 2);
4) calculating the D value of each cluster of microseism events according to the clustering condition in the step 2);
5) and calculating the fault activity index F according to the result of the b value in the step 3) and the result of the D value in the step 4).
Further, the method comprises the following steps:
1) inputting the microseism positioning result and the attribute information (x)i,yi,zi,Ti,Mi);
Wherein (x)i,yi,zi) Inverting spatial location coordinates, T, for microseismic eventsiGenerating time for microseismic events, MiFor the magnitude of the microseismic event, i belongs to [1, N ∈]N is the positioning number of the microseism events;
2) dividing the microseism events into R clusters according to the aggregation degree of the microseism events in the step 1); r is the number of clusters, SjFor the number of microseism events in each cluster, j belongs to [1, R ]],
3) Calculating the b value of each cluster of microseism events according to the clustering condition in the step 2), bjFor the j-th micro-seismic event b value, j belongs to [1, R ]];
log N=a-bM
In the formula: m is the microseism magnitude, N is the number of the microseism events of which the microseism magnitude is larger than M; a is a constant; b is the slope of magnitude and log of events at different magnitudes;
4) calculating the D value and D value of each cluster of microseism events according to the clustering condition of the step 2)jFor the j-th micro-seismic event D value, j belongs to [1, R ]];
C(r)∝rD
Wherein R is the microseism event distance, N is the number of microseism events when R is less than R, and D is the D value of the microseism event;
5) calculating a fault activity index F according to the result of the value b in the step 3) and the result of the value D in the step 4);
Fj=(|bj-1|+|Dj-2|)/2
wherein FjIs the fault activity index of the j-th micro-seismic event, when FjAnd when the number of the micro-seismic events is less than a certain constant c, the j clusters of micro-seismic events are considered to be fault activities.
The fault activity detection method can effectively identify the fault activity of the gas storage, avoid the fault sealing failure risk of the gas storage and ensure the safe operation of the gas storage.
The invention provides a method for identifying fault activities of an underground gas storage based on a microseism monitoring technology, which is based on microseism monitoring results, utilizes microseism event attribute information to calculate fault activity indexes of different cluster groups of microseism events, effectively identifies fault activities in the injection and production operation process of the underground gas storage, can effectively prevent the fault activity sealing failure risk of the underground gas storage and ensures the safe operation of the gas storage.
Description of the drawings:
FIG. 1 is a diagram of microseismic location and clustering results; the abscissa is the model X coordinate (unit: m); the ordinate is the model Y coordinate (unit: m), the size of the sphere represents the magnitude of the microseismic event, and the microseismic event is divided into 1 cluster because of more concentrated distribution.
FIG. 2 is a b-value calculation of a microseismic event; the magnitude of the abscissa (unit: Lei) and the number of microseismic events (unit: one) are the ordinate.
FIG. 3 is a D-value calculation of a microseismic event; the transverse distance r is logarithmic (unit: none), and the ordinate is the microseismic event distribution rule C (r) logarithmic (unit: none).
FIG. 4 is a graph of microseismic event fault activity index calculations; the abscissa is the cluster number (no units) and the ordinate is the microseismic event fault activity index F value (no units). The dashed black line is the value of a constant c, and when F < c, the cluster of microseismic events is identified as fault activity.
Detailed Description
The process of the invention comprises the following steps:
1) inputting the microseism positioning result and the attribute information (x)i,yi,zi,Ti,Mi);
Wherein (x)i,yi,zi) Inverting spatial location coordinates, T, for microseismic eventsiGenerating time for microseismic events, MiFor the magnitude of the microseismic event, i belongs to [1, N ∈]N is the positioning number of the microseism events;
2) dividing the microseism events into R clusters according to the aggregation degree of the microseism events in the step 1); r is the number of clusters, SjFor the number of microseism events in each cluster, j belongs to [1, R ]],
FIG. 1 shows the positioning and clustering results of microseismic events. The microseismic events in the area are more concentrated, so the microseismic events are divided into a cluster for analysis and research.
3) Calculating the b value of each cluster of microseism events according to the clustering condition in the step 2), bjFor the j-th micro-seismic event b value, j belongs to [1, R ]];
log N=a-bM
In the formula: m is the microseism magnitude, N is the number of the microseism events with the microseism magnitude larger than M; a is a constant; b is the slope of magnitude and log of events of different magnitude;
fig. 2 shows the calculation result of the b value of the input micro-seismic event in step 1), and the b value of the cluster of micro-seismic events in this embodiment is 1.2.
4) Calculating the D value and D value of each cluster of microseism events according to the clustering condition of the step 2)jFor the j-th micro-seismic event D value, j belongs to [1, R ]];
C(r)∝rD
Wherein R is the microseism event distance, N is the number of microseism events when R is less than R, and D is the D value of the microseism event;
fig. 3 is the calculation result of the D value of the input micro-seismic event in step 1), and the size of the D value of the micro-seismic event cluster is 2.02 in this embodiment.
5) Calculating a fault activity index F according to the result of the value b in the step 3) and the result of the value D in the step 4);
Fj=(|bj-1|+|Dj-2|)/2
wherein FjIs the fault activity index of the j-th micro-seismic event, when FjAnd when the number of the micro-seismic events is less than a certain constant c, the j clusters of micro-seismic events are considered to be fault activities.
FIG. 4 shows the calculation results of the F values of the fault activity indexes of different clusters. According to the formula, the fault activity index of the input micro-seismic event is 0.11, the value of the constant c is 0.25, and the input micro-seismic event cluster is judged to be fault activity.
Claims (1)
1. The method for identifying the fault activity of the underground gas storage based on the microseism monitoring technology is characterized by comprising the following steps of:
1) inputting a microseism positioning result and attribute information;
2) dividing the microseism events into R clusters according to the aggregation degree of the microseism events in the step 1);
3) calculating the b value of each cluster of microseism events according to the clustering condition in the step 2);
4) calculating the D value of each cluster of microseism events according to the clustering condition in the step 2);
5) calculating a fault activity index F according to the result of the value b in the step 3) and the result of the value D in the step 4);
the method specifically comprises the following steps:
1) inputting the microseism positioning result and the attribute information (x)i,yi,zi,Ti,Mi);
Wherein (x)i,yi,zi) Inverting spatial location coordinates, T, for microseismic eventsiGenerating time for microseismic events, MiFor the magnitude of the microseismic event, i belongs to [1, N ∈]N is the positioning number of the microseism events;
2) dividing the microseism events into R clusters according to the aggregation degree of the microseism events in the step 1); r is the number of clusters, SjFor the number of microseism events in each cluster, j belongs to [1, R ]],
3) Calculating the b value of each cluster of microseism events according to the clustering condition in the step 2), bjFor the j-th micro-seismic event b value, j belongs to [1, R ]];
log N=a-bM
In the formula: m is the microseism magnitude, N is the number of the microseism events with the microseism magnitude larger than M; a is a constant; b is the slope of magnitude and log of events of different magnitude;
4) calculating the D value and D value of each cluster of microseism events according to the clustering condition of the step 2)jFor the j-th micro-seismic event D value, j belongs to [1, R ]];
C(r)∝rD
Wherein R is the microseism event distance, N is the number of microseism events when R is less than R, and D is the D value of the microseism event;
5) calculating a fault activity index F according to the result of the value b in the step 3) and the result of the value D in the step 4);
Fj=(|bj-1|+|Dj-2|)/2
wherein FjIs the fault activity index of the j-th micro-seismic event, when FjAnd when the number of the micro-seismic events is less than a certain constant c, the j clusters of micro-seismic events are considered to be fault activities.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110518177.0A CN113253343B (en) | 2021-05-12 | 2021-05-12 | Method for identifying fault activity of underground gas storage based on microseism monitoring technology |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110518177.0A CN113253343B (en) | 2021-05-12 | 2021-05-12 | Method for identifying fault activity of underground gas storage based on microseism monitoring technology |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113253343A CN113253343A (en) | 2021-08-13 |
CN113253343B true CN113253343B (en) | 2022-05-31 |
Family
ID=77223083
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110518177.0A Active CN113253343B (en) | 2021-05-12 | 2021-05-12 | Method for identifying fault activity of underground gas storage based on microseism monitoring technology |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113253343B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117233835B (en) * | 2023-09-18 | 2024-05-28 | 北京戎彩科技有限公司 | Method for optimizing operation pressure interval of underground gas storage by microseism monitoring technology |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103760597A (en) * | 2013-12-27 | 2014-04-30 | 淮南万泰电子股份有限公司 | Automatic mine fault identification method |
FR3012623A1 (en) * | 2013-10-31 | 2015-05-01 | Cgg Services Sa | |
WO2017020461A1 (en) * | 2015-08-05 | 2017-02-09 | 深圳朝伟达科技有限公司 | Interpretation system of hydraulic fracturing micro-seismic event |
CN106483556A (en) * | 2016-10-09 | 2017-03-08 | 华北科技学院 | A kind of lasting earthquake magnitude based on Mine Earthquakes monitoring system and Richter scale conversion method |
CN107728200A (en) * | 2017-09-29 | 2018-02-23 | 中国石油化工股份有限公司 | Ground micro-seismic fracturing fracture dynamic spread method of real-time |
CN211954250U (en) * | 2020-04-30 | 2020-11-17 | 王娇华 | Gas storage monitoring system with microseism monitoring function |
CN112578457A (en) * | 2020-11-24 | 2021-03-30 | 中油奥博(成都)科技有限公司 | Optical fiber borehole earthquake monitoring method for carbon dioxide driven coal bed gas |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8213261B2 (en) * | 2008-05-22 | 2012-07-03 | Exxonmobil Upstream Research Company | Method for geophysical and geological interpretation of seismic volumes in the domains of depth, time, and age |
US8649980B2 (en) * | 2010-03-05 | 2014-02-11 | Vialogy Llc | Active noise injection computations for improved predictability in oil and gas reservoir characterization and microseismic event analysis |
CN111158045B (en) * | 2020-01-06 | 2022-02-22 | 中国石油化工股份有限公司 | Reservoir transformation microseism event scattered point clustering analysis method and system |
CN112464143A (en) * | 2020-10-23 | 2021-03-09 | 中国石油天然气集团有限公司 | Method and device for identifying underground coal in-situ gasification boundary |
-
2021
- 2021-05-12 CN CN202110518177.0A patent/CN113253343B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR3012623A1 (en) * | 2013-10-31 | 2015-05-01 | Cgg Services Sa | |
CN103760597A (en) * | 2013-12-27 | 2014-04-30 | 淮南万泰电子股份有限公司 | Automatic mine fault identification method |
WO2017020461A1 (en) * | 2015-08-05 | 2017-02-09 | 深圳朝伟达科技有限公司 | Interpretation system of hydraulic fracturing micro-seismic event |
CN106483556A (en) * | 2016-10-09 | 2017-03-08 | 华北科技学院 | A kind of lasting earthquake magnitude based on Mine Earthquakes monitoring system and Richter scale conversion method |
CN107728200A (en) * | 2017-09-29 | 2018-02-23 | 中国石油化工股份有限公司 | Ground micro-seismic fracturing fracture dynamic spread method of real-time |
CN211954250U (en) * | 2020-04-30 | 2020-11-17 | 王娇华 | Gas storage monitoring system with microseism monitoring function |
CN112578457A (en) * | 2020-11-24 | 2021-03-30 | 中油奥博(成都)科技有限公司 | Optical fiber borehole earthquake monitoring method for carbon dioxide driven coal bed gas |
Non-Patent Citations (2)
Title |
---|
A Time Picking Method for Microseismic Data Based on LLE and Improved PSO Clustering Algorithm;Haitao Ma;《IEEE Geoscience and Remote Sensing Letters》;20180725;全文 * |
页岩油水平井体积压裂及微地震监测技术实践;刘博 等;《岩性油气藏》;20200727;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113253343A (en) | 2021-08-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhu et al. | Natural gas pipeline valve leakage rate estimation via factor and cluster analysis of acoustic emissions | |
CN103268420B (en) | A kind of method for evaluating hazard of high rock slope | |
CN110298107B (en) | Working face impact risk evaluation method based on incremental stacking | |
CN113253344B (en) | Method for realizing pressure raising early warning of underground gas storage based on microseism monitoring technology | |
CN105389467A (en) | Method and apparatus of acquiring inter-well communication relationship | |
CN109254219B (en) | A kind of distribution transforming transfer learning method for diagnosing faults considering multiple factors Situation Evolution | |
CN113253343B (en) | Method for identifying fault activity of underground gas storage based on microseism monitoring technology | |
CN112364422B (en) | MIC-LSTM-based dynamic prediction method for shield construction earth surface deformation | |
CN201229289Y (en) | Corrosion predicting device | |
Jin et al. | Structural damage recognition based on filtered feature selection and convolutional neural network | |
Hao et al. | Quantification of margins and uncertainties for the risk of water inrush in a karst tunnel: representations of epistemic uncertainty with probability | |
Zhang et al. | Wireless monitoring–based real-time analysis and early-warning safety system for deep and large underground caverns | |
CN106246226A (en) | The recognition methods that a kind of Mine Gas Gushing is abnormal | |
CN103334740B (en) | Consider the method for the determination drainage front of free-boundary problem | |
Liu et al. | Study of roof water inrush forecasting based on EM-FAHP two-factor model | |
Guo et al. | Application of an improved cloud model and distance discrimination to evaluate slope stability | |
Jun et al. | A new dynamic assessment for multi-parameters information of water inrush in coal mine | |
Hong et al. | Evaluation of disaster-bearing capacity for natural gas pipeline under third-party damage based on optimized probabilistic neural network | |
CN117233835B (en) | Method for optimizing operation pressure interval of underground gas storage by microseism monitoring technology | |
CN116307706B (en) | Timing sequence model-based gas storage operation geological risk early warning and gas injection and production scheme optimization method | |
Tian et al. | [Retracted] Prediction of Coal Mining Subsidence Based on Machine Learning Probability Theory | |
CN104881536B (en) | A kind of rock mass discontinuity goodness of fit coefficient measuring method | |
Han et al. | Novel Long Short-Term Memory Model Based on the Attention Mechanism for the Leakage Detection of Water Supply Processes | |
Mohaghegh | Top-down, intelligent reservoir model | |
Lian et al. | Design and implementation of mine water hazard monitoring and early warning platform |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |