CN106353805A - Intuitive microearthquake event classification method - Google Patents
Intuitive microearthquake event classification method Download PDFInfo
- Publication number
- CN106353805A CN106353805A CN201510422714.6A CN201510422714A CN106353805A CN 106353805 A CN106353805 A CN 106353805A CN 201510422714 A CN201510422714 A CN 201510422714A CN 106353805 A CN106353805 A CN 106353805A
- Authority
- CN
- China
- Prior art keywords
- micro
- arrival
- seismic event
- event
- microearthquake
- 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.)
- Granted
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses an intuitive microearthquake event classification method. The technical scheme includes: firstly, performing de-noising processing on original observation records and increasing signal-to-noise ratio of the records; finding out possible microearthquake events by a long time window and short time window proportion method; collecting first arrivals of the microearthquake events according to continuity of synthetic shafts, performing arrangement processing on the acquired first arrivals and finally superimposing all the first arrivals of the microearthquake events to display. According to distribution of the superimposed first arrivals, the microearthquake events similar in first arrival forms are defined as the same microearthquake event group, unified selection of parameters in post-processing is facilitated, and the events with the curve forms different from the overall trend are defined as stray events or invalid microearthquake events. The intuitive microearthquake event classification method has the advantages that classification of the microearthquake events is convenient, intuitive and reliable.
Description
Technical field
The present invention relates to fracturing micro-seismic monitoring field, specifically one kind intuitively microseism
Event classification method.
Background technology
Hydraulic fracturing technology is the main path improving the unconventionaloil pool reservoir recovery ratios such as shale oil gas, micro-ly
Shake monitoring is by gathering the signal that fracturing causes catalase to produce, rupture location being positioned, thus
Instruct pressing crack construction and pressure break result is evaluated.Typically pick up in existing microseism data processor
Directly positioned after taking micro-seismic event first arrival, then the distribution situation according to micro-seismic event positioning result
Carry out the division of microseismic event.There is following point in this handling process: one is not have priori to microseism
Event is classified, and would generally select identical positional parameter in position fixing process, so often affects micro-
The precision of seismic events positioning;Two is that micro-seismic event to all pickups positions, so can be greatly
Waste computing resource, it is particularly evident that this embodies during real-time monitoring.
Content of the invention
The purpose of the present invention is for solving problem present in above-mentioned tradition microseism data handling procedure, improving
The location efficiency of micro-seismic event and positioning precision, provide a kind of convenient, intuitively micro-seismic event classification
Method.
The present invention can by the following method measure realizing:
One kind intuitively micro-seismic event sorting technique, comprising:
1) protocol is carried out with denoising, improves the signal to noise ratio of record;
2) micro-seismic event is found out according to energy ratio between the long window and the short window method;
3) determine the first arrival of micro-seismic event;
4) the micro-seismic event first arrival in step 3 is carried out regular process;
5) by step 4) in microseism first arrival be plotted in the same coordinate system;
6) according to step 5) in the form of coordinate system first arrival micro-seismic event is classified.
Such scheme further includes:
The first arrival of the determination micro-seismic event described in described step 3 comprises the steps of
A) according to micro-seismic event according to the seriality of axle in the same direction, first break picking;
B) the corresponding crest of first break picking method choice pickup, trough or three kinds of modes of zero crossing, will
Corresponding sampling point position is multiplied by the sampling interval to be recorded as first arrival.
The micro-seismic event first arrival of described step 4 carries out regular comprising the steps of
C) set first most [t1 t2 ... the tn] of micro-seismic event, in formula, n is cymoscope series, meter
Calculate the average of pickup micro-seismic event first arrival, then average utilizes following formula to calculate
D) the first arrival unification of micro-seismic event is deducted the average calculating in step c), advised
First arrival [t1 t2 ... tn]-t after wholeave.
Described step 6 is classified to micro-seismic event, the following sorting criterion of application:
E) using the micro-seismic event with shooting point first arrival plesiomorphism as validity event;
F) first arrival plesiomorphism micro-seismic event is divided into same event group;
G) the microseism leave of absence larger for first arrival morphological differencess is defined as the event of peeling off or anomalous event.
One kind of the present invention intuitively micro-seismic event sorting technique, by original observational record
Carry out denoising, improve the signal to noise ratio of record, then find out possibility using long short time-window than method
Micro-seismic event, pick up the first arrival of micro-seismic event according to the seriality of axle in the same direction, will obtain
First arrival carry out regular process, finally all of micro-seismic event first arrival is overlapped display. pass through
The first arrival distribution situation of the micro-seismic event of superposition, the micro-seismic event of first arrival plesiomorphism is fixed
Justice is same micro-seismic event group, to facilitate the unified selection entering line parameter in last handling process,
The event definition simultaneously tracing pattern being different from overall trend is the event that peels off or invalid micro-
Shake event.The feature of the method is convenient, directly perceived, reliable for the classification of micro-seismic event.
Brief description
Accompanying drawing 1 is this invention flow chart.
Accompanying drawing 2 is original three-component borehole microseismic logout.
Record section after accompanying drawing 3 protocol denoising.Wherein grey lines are according to axle in the same direction
The first arrival of seriality pickup.
Accompanying drawing 4 renders to display in the same coordinate system for regular rear microseism first arrival.Wherein: 1 is
The corresponding first arrival of micro-seismic event in Fig. 3,2,3 is two micro-seismic event groups, and 4,5 are two
Peel off event.
Specific embodiment
Above and other objects, features and advantages for enabling the present invention become apparent, hereafter
Especially exemplified by going out preferred embodiment, and coordinate institute's accompanying drawings, be described in detail below.
One kind intuitively micro-seismic event sorting technique, comprising:
1) protocol is carried out with denoising, improves the signal to noise ratio of record;
2) micro-seismic event is found out according to energy ratio between the long window and the short window method;
3) determine the first arrival of micro-seismic event;
4) the micro-seismic event first arrival in step 3 is carried out regular process;
5) by step 4) in microseism first arrival be plotted in the same coordinate system;
6) according to step 5) in the form of coordinate system first arrival micro-seismic event is classified.
Such scheme further includes:
The first arrival of the determination micro-seismic event described in described step 3 comprises the steps of
A) according to micro-seismic event according to the seriality of axle in the same direction, first break picking;
B) the corresponding crest of first break picking method choice pickup, trough or three kinds of modes of zero crossing, will
Corresponding sampling point position is multiplied by the sampling interval to be recorded as first arrival.
The micro-seismic event first arrival of described step 4 carries out regular comprising the steps of
C) set first most [t1 t2 ... the tn] of micro-seismic event, in formula, n is cymoscope series, meter
Calculate the average of pickup micro-seismic event first arrival, then average utilizes following formula to calculate
D) the first arrival unification of micro-seismic event is deducted the average calculating in step c), advised
First arrival [t1 t2 ... tn]-t after wholeave.
Described step 6 is classified to micro-seismic event, the following sorting criterion of application:
E) using the micro-seismic event with shooting point first arrival plesiomorphism as validity event;
F) first arrival plesiomorphism micro-seismic event is divided into same event group;
G) the microseism leave of absence larger for first arrival morphological differencess is defined as the event of peeling off or anomalous event.
1-4 is described further to above-mentioned embodiment below in conjunction with the accompanying drawings.
First, original three-component microseismograms (as shown in Figure 2) is carried out by direct current, goes
The denoisings such as exceptional value, bandpass filtering, obtain the record shown in Fig. 3.
Second, it is as shown in Figure 3 to apply energy ratio between the long window and the short window method to find out micro-seismic event.
Third, picking up the first arrival of outgoing event according to the consistent principle of crest, (as Fig. 3 gray line institute
Show), calculate the true first arrival of this micro-seismic event according to the sampling interval simultaneously.
Fourth, by step 3 pickup first arrival carry out regular process and project to coordinate system, such as
In Fig. 4 shown in mark 1.
Fifth, repeat step 3-4, the first arrival of all of regular micro-seismic event is thrown the most at last
It is put in coordinate as shown in Figure 4.
Finally, according to the first arrival form according to micro-seismic event in Fig. 4 coordinate system, according to invention
Judgment criterion, can by micro-seismic event be divided into 2 microseismic event groups (in Fig. 4 mark 2,
3) and 2 events that peel off (mark 4,5) in Fig. 4, it is achieved thereby that micro-seismic event point
Class.
Claims (4)
1. one kind intuitively micro-seismic event sorting technique it is characterised in that include,
1) protocol is carried out with denoising, improves the signal to noise ratio of record;
2) micro-seismic event is found out according to energy ratio between the long window and the short window method;
3) determine the first arrival of micro-seismic event;
4) the micro-seismic event first arrival in step 3 is carried out regular process;
5) by step 4) in microseism first arrival be plotted in the same coordinate system;
6) according to step 5) in the form of coordinate system first arrival micro-seismic event is classified.
2. the intuitively micro-seismic event sorting technique according to claim 1, its feature
It is that the first arrival of the determination micro-seismic event described in described step 3 comprises the steps of
A) according to micro-seismic event according to the seriality of axle in the same direction, first break picking;
B) the corresponding crest of first break picking method choice pickup, trough or three kinds of modes of zero crossing, will
Corresponding sampling point position is multiplied by the sampling interval to be recorded as first arrival.
3. the intuitively micro-seismic event sorting technique according to claim 2, its feature exists
Micro-seismic event first arrival in described step 4 carries out regular comprising the steps of
C) set first most [t1 t2 ... the tn] of micro-seismic event, in formula, n is cymoscope series, meter
Calculate the average of pickup micro-seismic event first arrival, then average utilizes following formula to calculate
D) the first arrival unification of micro-seismic event is deducted the average calculating in step c), advised
First arrival [t1 t2 ... tn]-t after wholeave.
4. the intuitively micro-seismic event sorting technique according to claim 3, its feature exists
In described step 6, micro-seismic event is classified, the following sorting criterion of application:
E) using the micro-seismic event with shooting point first arrival plesiomorphism as validity event;
F) first arrival plesiomorphism micro-seismic event is divided into same event group;
G) the microseism leave of absence larger for first arrival morphological differencess is defined as the event of peeling off or abnormal thing
Part.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510422714.6A CN106353805B (en) | 2015-07-17 | 2015-07-17 | Intuitive microseism event classification method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510422714.6A CN106353805B (en) | 2015-07-17 | 2015-07-17 | Intuitive microseism event classification method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106353805A true CN106353805A (en) | 2017-01-25 |
CN106353805B CN106353805B (en) | 2019-12-20 |
Family
ID=57842296
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510422714.6A Active CN106353805B (en) | 2015-07-17 | 2015-07-17 | Intuitive microseism event classification method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106353805B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109143364A (en) * | 2017-06-27 | 2019-01-04 | 中国石油化工股份有限公司 | First arrival time pick-up method |
CN110286408A (en) * | 2019-06-28 | 2019-09-27 | 广东石油化工学院 | A kind of microseismic event detection method and system |
CN110886599A (en) * | 2018-09-07 | 2020-03-17 | 中国石油化工股份有限公司 | Non-fracturing event identification method and system based on fracture speed |
CN111881299A (en) * | 2020-08-07 | 2020-11-03 | 哈尔滨商业大学 | Outlier event detection and identification method based on duplicate neural network |
CN113514878A (en) * | 2021-07-09 | 2021-10-19 | 长沙矿山研究院有限责任公司 | Mine microseismic waveform type automatic identification method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103837891A (en) * | 2014-03-26 | 2014-06-04 | 宋维琪 | High-precision pickup method of first arrival of microseism |
CN104216008A (en) * | 2013-06-05 | 2014-12-17 | 中国石油天然气集团公司 | Downhole fracturing microseismic event identification method |
CN104730579A (en) * | 2013-12-18 | 2015-06-24 | 中国石油化工股份有限公司 | P-wave and S-wave combined static correction method based on surface S-wave velocity inversion |
-
2015
- 2015-07-17 CN CN201510422714.6A patent/CN106353805B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104216008A (en) * | 2013-06-05 | 2014-12-17 | 中国石油天然气集团公司 | Downhole fracturing microseismic event identification method |
CN104730579A (en) * | 2013-12-18 | 2015-06-24 | 中国石油化工股份有限公司 | P-wave and S-wave combined static correction method based on surface S-wave velocity inversion |
CN103837891A (en) * | 2014-03-26 | 2014-06-04 | 宋维琪 | High-precision pickup method of first arrival of microseism |
Non-Patent Citations (1)
Title |
---|
魏路路: "水力压裂微地震事件自动识别与拾取", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109143364A (en) * | 2017-06-27 | 2019-01-04 | 中国石油化工股份有限公司 | First arrival time pick-up method |
CN109143364B (en) * | 2017-06-27 | 2020-04-14 | 中国石油化工股份有限公司 | First arrival time picking method |
CN110886599A (en) * | 2018-09-07 | 2020-03-17 | 中国石油化工股份有限公司 | Non-fracturing event identification method and system based on fracture speed |
CN110286408A (en) * | 2019-06-28 | 2019-09-27 | 广东石油化工学院 | A kind of microseismic event detection method and system |
CN111881299A (en) * | 2020-08-07 | 2020-11-03 | 哈尔滨商业大学 | Outlier event detection and identification method based on duplicate neural network |
CN111881299B (en) * | 2020-08-07 | 2023-06-30 | 哈尔滨商业大学 | Outlier event detection and identification method based on replicated neural network |
CN113514878A (en) * | 2021-07-09 | 2021-10-19 | 长沙矿山研究院有限责任公司 | Mine microseismic waveform type automatic identification method |
CN113514878B (en) * | 2021-07-09 | 2022-06-24 | 长沙矿山研究院有限责任公司 | Mine microseismic waveform type automatic identification method |
Also Published As
Publication number | Publication date |
---|---|
CN106353805B (en) | 2019-12-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106353805A (en) | Intuitive microearthquake event classification method | |
CN110680302B (en) | Automatic identification method for electrocardiosignal characteristic wave | |
CN108664841A (en) | A kind of sound state object recognition methods and device based on laser point cloud | |
CN105785439B (en) | The Forecasting Methodology and device of small scale heterogeneous geologic body spatial distribution position | |
CN104266894A (en) | Mine microearthquake signal preliminary wave moment extracting method based on correlation analysis | |
CN104459779A (en) | Method and device for automatically identifying abnormal seismic traces | |
CN113092972B (en) | Partial discharge type identification method and device based on non-contact ultrasonic sensor | |
CN104851301B (en) | Vehicle parameter identification method based on deceleration strip sound analysis | |
CN105022027B (en) | A kind of data filtering method and device | |
CN107402405B (en) | Quiet phase virtual source trace gather construction method | |
CN109696609A (en) | A kind of acquisition methods and its acquisition system of novel bow net current carrying quality | |
CN104714251B (en) | Oblique stacking peak amplitude portion edge detecting method for event automatic picking | |
CN102254321B (en) | Method for automatically identifying reversed polarity channel based on first-motion wave | |
CN102590359B (en) | Method and system for identifying redundancy signals | |
CN104597486A (en) | Method for quickly detecting shot gather data through first arrival time | |
CN109088762A (en) | A kind of service health monitor method based on micro services | |
CN106646592B (en) | The real-time stacking method of single shot record | |
CN104932742A (en) | Method and system for obtaining multi-point operation tracks of touch screen device | |
CN116466396A (en) | Nonlinear optimal positioning method for microseism focus | |
CN108919068B (en) | Intermittent defect signal identification method for power equipment | |
CN112394402A (en) | Method and system for detecting microseism signals based on synchronous extrusion wavelet transform | |
CN102998699A (en) | Casing-harmonic-containing vertical seismic profile data interval velocity inversion method | |
CN104570118A (en) | Method for automatically identifying and removing industrial interface based on dual factors | |
CN113217109B (en) | Waveform completion method for rockburst mine microseismic monitoring system | |
CN112764109B (en) | Separation and extraction method and device for dipole transverse wave reflected waves |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |