CN106353805A - Intuitive microearthquake event classification method - Google Patents

Intuitive microearthquake event classification method Download PDF

Info

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
Application number
CN201510422714.6A
Other languages
Chinese (zh)
Other versions
CN106353805B (en
Inventor
陈雨茂
左建军
武小杰
魏国华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Geophysical Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Geophysical Research Institute of Sinopec Shengli Oilfield Co
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Geophysical Research Institute of Sinopec Shengli Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN201510422714.6A priority Critical patent/CN106353805B/en
Publication of CN106353805A publication Critical patent/CN106353805A/en
Application granted granted Critical
Publication of CN106353805B publication Critical patent/CN106353805B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

One kind intuitively micro-seismic event sorting technique
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
t a v e = σ i = 1 n t i / n - - - ( 1 )
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
t a v e = σ i = 1 n t i / n - - - ( 1 )
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
t a v e = σ i = 1 n t i / n - - - ( 1 )
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.
CN201510422714.6A 2015-07-17 2015-07-17 Intuitive microseism event classification method Active CN106353805B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
魏路路: "水力压裂微地震事件自动识别与拾取", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *

Cited By (8)

* Cited by examiner, † Cited by third party
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