CN109581481A - A kind of portable high frequency controlled source seismic signal harmonic wave interference removing method - Google Patents
A kind of portable high frequency controlled source seismic signal harmonic wave interference removing method Download PDFInfo
- Publication number
- CN109581481A CN109581481A CN201910018313.2A CN201910018313A CN109581481A CN 109581481 A CN109581481 A CN 109581481A CN 201910018313 A CN201910018313 A CN 201910018313A CN 109581481 A CN109581481 A CN 109581481A
- Authority
- CN
- China
- Prior art keywords
- signal
- seismic
- high frequency
- wave
- arrival
- 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
- 238000000034 method Methods 0.000 title claims abstract description 15
- 239000000758 substrate Substances 0.000 claims abstract description 20
- 230000005284 excitation Effects 0.000 claims abstract description 13
- 239000000284 extract Substances 0.000 claims abstract description 7
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 6
- 230000003213 activating effect Effects 0.000 claims description 4
- 230000021615 conjugation Effects 0.000 claims description 3
- 230000035939 shock Effects 0.000 claims 1
- 230000008030 elimination Effects 0.000 abstract description 5
- 238000003379 elimination reaction Methods 0.000 abstract description 5
- 238000001514 detection method Methods 0.000 abstract description 4
- 238000000605 extraction Methods 0.000 abstract 1
- 230000008901 benefit Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005674 electromagnetic induction Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
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/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
-
- 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention relates to a kind of portable high frequency controlled source seismic signal harmonic wave interference removing methods, construct according to the signal excitation and collecting flowchart of conventional Seismic Exploration with Vibrator, away within focus substrate 1m along line direction arrange one of wave detector.After epicenter excitation, from the signal that the near earthquake source substrate wave detector especially arranged acquires, according to the then different and capacity volume variance of first arrival signal and other seismic waves, identifies first arrival signal and extract.Cross correlation process is carried out using the first arrival signal of extraction as reference factor and the earthquake record of acquisition, the earthquake record of wavelet form is obtained, harmonic wave interference can be effectively eliminated.The present invention is suitable for eliminating the harmonic wave interference of portable high frequency controlled source seismic signal, the axis in the same direction of cross-correlation wavelet is more continuous, and seismic phase is clear, and signal-to-noise ratio is obviously improved, particularly with complex geologic conditions, substrate and the Harmonics elimination in the bad situation of ground-coupled are more effective.The present invention provides effective technical support for the Shallow High Resolutions exploration such as urban underground space detection.
Description
Technical field
The present invention relates to field of seismic exploration, especially a kind of portable high frequency controlled source seismic signal harmonic wave interference disappears
Except method.The harmonic wave interference of portable controlled source seismic signal can be eliminated using the present invention, particularly with complex geologic conditions,
Substrate and the Harmonics elimination in the bad situation of ground-coupled are more effective.
Background technique
Portable high frequency controlled source is a kind of seismic prospecting based on electromagnetic induction principle driving vibration excitor high-frequency vibration
Excitaton source.Compared with the current fluid pressure type controlled source routinely applied, it is wide in range that feature shows as excitation signal frequency band, usually may be used
It realizes the vibration for being higher than 1000Hz, while lightweight and portable in volume is mobile flexibly, is widely used in engineering geophysics, urban underground space detection
Equal Shallow High Resolutions seismic prospecting.
For the earthquake record of portable high frequency controlled source, wavelet similar with Pulse Source in order to obtain, take by
Earthquake record and reference factor carry out the preprocess method of cross-correlation.Current conventional method is by focus scanning signal or nearly substrate
Geophone signal is verified as reference factor by experimental data, and two methods have different degrees of harmonic wave interference.
Especially in the case where the regional substrate of complex geologic conditions and ground-coupled are out of condition, harmonic wave interference is especially prominent, serious to drop
The low resolution ratio of portable high frequency controlled source seismic signal, limits further applying for instrument.
Summary of the invention
The purpose of the present invention is to solve defect of the existing technology, a kind of portable high frequency controlled source is provided
Shake signal harmonic interference elimination method.
To achieve the goals above, the technical solution adopted by the present invention is that:
A kind of portable high frequency controlled source seismic signal harmonic wave interference removing method, includes the following steps:
A, in survey area, conventionally the requirement of seismic exploration designs activating system, the focus ginseng of each excitation point
Number is consistent, and scanning signal is s (t), and focus substrate and ground is kept to form good coupling；
B, conventionally the requirement of seismic exploration designs observation system, if along arrangement of measuring-line arterial highway wave detector, each road detection
The distance between device is identical, and wave detector and ground is kept to form good coupling；
C, when each focal point excites, one of wave detector is arranged along line direction within away from focus substrate 1m, referred to as closely
Focus substrate wave detector；
D, by portable high frequency controlled source earthquake-wave-exciting, the wave detector through arranging on survey line acquires seismic signal,
And saved by seismic detector, obtain original non-relative earthquake record.
Further, the invention also includes following steps:
E, the non-relative earthquake obtained for i excitation point records, and extracts the 1st signal, i.e. near earthquake source substrate wave detector
Signal gi(t).By gi(t) computing cross-correlation is carried out with focus scanning signal s (t), obtainedWhereinSymbol
Represent computing cross-correlation；
F, from yji(t) in, according to first arrival signal identification first arrival signal then different from other seismic waves.Wherein, by yji
(t) the crest value point in is considered first arrival signal then, and other seismic waves include back wave, refracted wave, surface wave and make an uproar at random
Sound；
G, with the first arrival signal of f step identification then for midpoint, according to the autocorrelation wavelet pulsewidth of scanning signal s (t), yji
(t) waveform morphology and the requirement of energy accounting, left and right respectively takes the window of same time length, extracts first arrival signal, obtains yfi
(t).Wherein, it is desirable that contain only first arrival signal in window and do not include other signals, while signal energy accounts for entirely just in window
To 90% or more of signal energy；
H, yf at this timei(t) it is obtained for cross-correlation wavelet form according to the relationship of computing cross-correlation and convolutionWherein, gfiIt (t) is the first arrival signal of FM signal form, symbol * represents convolution algorithm,It represents to letter
Number seek conjugation；
I, according to convolution theorem, Fourier transformation is done to above formula, is obtainedThenWherein F [] indicates to carry out Fourier transformation to signal, then carries out Fourier inversion,
Obtain gfi(t)=F-1[F[gfi(t)]], wherein F-1[] indicates to carry out Fourier inversion to signal；
J, by gfi(t) it is used as reference factor, the non-relative earthquake record obtained with i excitation point carries out computing cross-correlation,
The earthquake record of the wavelet form finally required.
The invention has the benefit that the harmonic wave that the present invention is suitable for eliminating portable high frequency controlled source seismic signal is dry
It disturbs, compared with existing conventional preprocess method, the wavelet of computing cross-correlation result axis in the same direction is more continuous, and seismic phase is clear, signal-to-noise ratio
It is obviously improved, particularly with complex geologic conditions, substrate and the Harmonics elimination in the bad situation of ground-coupled are more effective.This hair
It is bright to provide effective technical support for the Shallow High Resolutions exploration such as urban underground space detection.
Detailed description of the invention
Fig. 1 a is the portable high frequency controlled source earthquake record of conventional method processing；
Fig. 1 b is the time-frequency curve of the 10th signal of Fig. 1 a；
Fig. 2 a is the portable high frequency controlled source earthquake record of the method for the present invention processing；
Fig. 2 b is the time-frequency curve of the 10th signal of Fig. 2 a.
Specific embodiment
It is described in further detail with reference to the accompanying drawings and examples:
Such as Fig. 1 a, shown in Fig. 1 b, Fig. 2 a, Fig. 2 b, for northern China urban applications portable high frequency controlled source into
For the earthquake record that the row underground space detects, illustrate portable high frequency controlled source seismic signal related to the present invention
Harmonic wave interference removing method, including following order and step:
A, in survey area, conventionally the requirement of seismic exploration designs activating system, distance between each excitation point
6m, the scanning signal that focus is arranged is s (t), 20-200Hz of frequency range, sweep time 10s, sample frequency 1KHz, when excitation
Focus substrate and ground is kept to form good coupling；
B, conventionally the requirement of seismic exploration designs observation system, and along 24 wave detector of arrangement of measuring-line, offset distance is
10m, the distance between each road wave detector are 2m, it is desirable that wave detector and ground form good coupling；
C, special, when every focal point excites, at away from focus substrate 0.5m, along line direction arrangement together nearly substrate inspection
Wave device；
D, after activating system and observation system are provided with, earthquake is excited to the earth by portable high frequency controlled source
Wave, the wave detector through arranging on survey line acquires seismic signal, and is stored by seismic detector by the sample rate of 4K, obtains original
Non- relative earthquake record；
E, the non-relative earthquake obtained for No. 1 excitation point records, and extracts the 1st signal, i.e. near earthquake source substrate wave detector
Signal g1(t).By g1(t) computing cross-correlation is carried out with focus scanning signal s (t), obtainedWhereinSymbol
Represent computing cross-correlation；
F, from yj1(t) in, according to first arrival signal and other seismic waves (such as back wave, refracted wave, surface wave and random noise
Equal signals) then different identification first arrival signals.Wherein, by yj1(t) the crest value point 0.52ms in is considered first arrival signal
Then；
G, using the first arrival signal of f step identification, then 0.52ms is midpoint, according to the autocorrelation wavelet arteries and veins of scanning signal s (t)
Wide 8.1ms, and according to waveform morphology, left and right respectively takes the window of 4.5ms length, makes to contain only first arrival signal in window and not include
Other signals extract first arrival signal, obtain yf1(t).Meanwhile the signal energy extracted in calculation window, account for entire first arrival signal
The 91.36% of energy meets requirement of the invention；
H, the yf extracted at this time1(t) it is obtained for cross-correlation wavelet form according to the relationship of computing cross-correlation and convolution
It arrivesWherein, gfiIt (t) is the first arrival signal of FM signal form, symbol * represents convolution algorithm,Representative pair
Signal seeks conjugation；
I, according to convolution theorem, Fourier transformation is done to above formula, is obtainedThenWherein F [] indicates to carry out Fourier transformation to signal.Fourier inversion is carried out again,
Obtain gf1(t)=F-1[F[gf1(t)]], wherein F-1[] indicates to carry out Fourier inversion to signal.
J, by gf1(t) it is used as reference factor, the non-relative earthquake record obtained with No. 1 excitation point carries out computing cross-correlation,
Obtain the earthquake record for the wavelet form that the present invention finally requires.
The present invention is suitable for eliminating the harmonic wave interference of portable high frequency controlled source seismic signal, with existing conventional pretreatment
Method compares, and the wavelet of computing cross-correlation result axis in the same direction is more continuous, and seismic phase is clear, and signal-to-noise ratio is obviously improved, particularly with ground
Matter complicated condition, substrate and the Harmonics elimination in the bad situation of ground-coupled are more effective.The present invention is urban underground space spy
The exploration of the Shallow High Resolutions such as survey provides effective technical support.
The basic principles, main features and advantages of the present invention have been shown and described above.The technology of the industry
Personnel are it should be appreciated that the present invention is not limited to the above embodiments, and what is described in the above embodiment and the description is only the present invention
Principle, various changes and improvements may be made to the invention without departing from the spirit and scope of the present invention, these variation and
Improvement is both fallen in the range of claimed invention.The present invention claims protection scope by appended claims and its
Equivalent defines.
Claims (2)
1. a kind of portable high frequency controlled source seismic signal harmonic wave interference removing method, which comprises the steps of:
A, in survey area, conventionally the requirement of seismic exploration designs activating system, the focal shock parameter one of each excitation point
It causes, scanning signal is s (t), and focus substrate and ground is kept to form good coupling；
B, conventionally seismic exploration requirement design observation system, if along arrangement of measuring-line arterial highway wave detector, each road wave detector it
Between distance it is identical, and wave detector and ground is kept to form good coupling；
C, when each focal point excites, one of wave detector, referred to as near earthquake source are arranged along line direction within away from focus substrate 1m
Substrate wave detector；
D, by portable high frequency controlled source earthquake-wave-exciting, the wave detector acquisition seismic signal through being arranged on survey line, and by
Seismic detector saves, and obtains original non-relative earthquake record.
2. a kind of portable high frequency controlled source seismic signal harmonic wave interference removing method according to claim 1, special
Sign is, further includes following steps:
E, the non-relative earthquake obtained for i excitation point records, and extracts the 1st signal, i.e. near earthquake source substrate geophone signal gi
(t).By gi(t) computing cross-correlation is carried out with focus scanning signal s (t), obtainedWhereinSymbology is mutual
Related operation；
F, from yji(t) in, according to first arrival signal identification first arrival signal then different from other seismic waves.Wherein, by yji(t)
In crest value point be considered first arrival signal then, other seismic waves include back wave, refracted wave, surface wave and random noise；
G, with the first arrival signal of f step identification then for midpoint, according to the autocorrelation wavelet pulsewidth of scanning signal s (t), yji(t)
Waveform morphology and the requirement of energy accounting, left and right respectively take the window of same time length, extract first arrival signal, obtain yfi(t).Its
In, it is desirable that first arrival signal is contained only in window and does not include other signals, while signal energy accounts for entire first arrival signal in window
90% or more of energy；
H, yf at this timei(t) it is obtained for cross-correlation wavelet form according to the relationship of computing cross-correlation and convolution
Wherein, gfiIt (t) is the first arrival signal of FM signal form, symbol * represents convolution algorithm,Representative seeks conjugation to signal；
I, according to convolution theorem, Fourier transformation is done to above formula, is obtainedThenWherein F [] indicates to carry out Fourier transformation to signal, then carries out Fourier inversion,
Obtain gfi(t)=F-1[F[gfi(t)]], wherein F-1[] indicates to carry out Fourier inversion to signal；
J, by gfi(t) it is used as reference factor, the non-relative earthquake record obtained with i excitation point carries out computing cross-correlation, obtains
The earthquake record of the wavelet form finally required.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910018313.2A CN109581481B (en) | 2019-01-09 | 2019-01-09 | Portable high-frequency controllable seismic source seismic signal harmonic interference elimination method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910018313.2A CN109581481B (en) | 2019-01-09 | 2019-01-09 | Portable high-frequency controllable seismic source seismic signal harmonic interference elimination method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109581481A true CN109581481A (en) | 2019-04-05 |
| CN109581481B CN109581481B (en) | 2021-06-11 |
Family
ID=65916401
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910018313.2A Active CN109581481B (en) | 2019-01-09 | 2019-01-09 | Portable high-frequency controllable seismic source seismic signal harmonic interference elimination method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109581481B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111487682A (en) * | 2020-06-03 | 2020-08-04 | 东营市震点石油科技有限公司 | Surface layer investigation method based on node seismograph |
| CN111948703A (en) * | 2019-05-17 | 2020-11-17 | 中国石油天然气集团有限公司 | Seismic exploration method and device excited by mixed seismic source |
| CN113050166A (en) * | 2019-12-27 | 2021-06-29 | 中国石油天然气集团有限公司 | Seismic excitation signal correction method and device |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100157733A1 (en) * | 2008-12-23 | 2010-06-24 | Mark Daly | Compensating seismic data for source variations |
| CN103823243A (en) * | 2012-11-15 | 2014-05-28 | 地球物理集团公司 | Process for separating data recorded during continuous data acquisition seismic survey |
| US20140249757A1 (en) * | 2013-03-04 | 2014-09-04 | Bruno Gratacos | Apparatus and method for determination of far-field signature from variable-depth seismic data |
| CN105277983A (en) * | 2014-06-26 | 2016-01-27 | 中石化石油工程地球物理有限公司胜利分公司 | Novel method for suppressing vibroseis slide scanning harmonic interference |
| CN105319587A (en) * | 2014-07-04 | 2016-02-10 | 中国石油化工股份有限公司 | Vibroseis shaping filtering method based on first arrival wavelet extraction |
| CN108957551A (en) * | 2018-07-03 | 2018-12-07 | 吉林大学 | Controlled source harmonic suppression method based on reconstruct ground force signal |
-
2019
- 2019-01-09 CN CN201910018313.2A patent/CN109581481B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100157733A1 (en) * | 2008-12-23 | 2010-06-24 | Mark Daly | Compensating seismic data for source variations |
| CN103823243A (en) * | 2012-11-15 | 2014-05-28 | 地球物理集团公司 | Process for separating data recorded during continuous data acquisition seismic survey |
| US20140249757A1 (en) * | 2013-03-04 | 2014-09-04 | Bruno Gratacos | Apparatus and method for determination of far-field signature from variable-depth seismic data |
| CN105277983A (en) * | 2014-06-26 | 2016-01-27 | 中石化石油工程地球物理有限公司胜利分公司 | Novel method for suppressing vibroseis slide scanning harmonic interference |
| CN105319587A (en) * | 2014-07-04 | 2016-02-10 | 中国石油化工股份有限公司 | Vibroseis shaping filtering method based on first arrival wavelet extraction |
| CN108957551A (en) * | 2018-07-03 | 2018-12-07 | 吉林大学 | Controlled source harmonic suppression method based on reconstruct ground force signal |
Non-Patent Citations (1)
| Title |
|---|
| 杨志超: ""电磁式可控震源地震信号检测关键技术研究"", 《中国博士学位论文全文数据库•基础科学辑》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111948703A (en) * | 2019-05-17 | 2020-11-17 | 中国石油天然气集团有限公司 | Seismic exploration method and device excited by mixed seismic source |
| CN111948703B (en) * | 2019-05-17 | 2023-07-25 | 中国石油天然气集团有限公司 | Seismic exploration method and device for mixed seismic source excitation |
| CN113050166A (en) * | 2019-12-27 | 2021-06-29 | 中国石油天然气集团有限公司 | Seismic excitation signal correction method and device |
| CN111487682A (en) * | 2020-06-03 | 2020-08-04 | 东营市震点石油科技有限公司 | Surface layer investigation method based on node seismograph |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109581481B (en) | 2021-06-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Schimmel et al. | Using instantaneous phase coherence for signal extraction from ambient noise data at a local to a global scale | |
| Li et al. | A method for low-frequency noise suppression based on mathematical morphology in microseismic monitoring | |
| Tian et al. | Variable-eccentricity hyperbolic-trace TFPF for seismic random noise attenuation | |
| CN101598812B (en) | Method for removing abnormal noise in single-point reception of seismic record by digital detector | |
| CN109581481A (en) | A kind of portable high frequency controlled source seismic signal harmonic wave interference removing method | |
| CN110529087B (en) | Method and device for evaluating hydraulic fracturing effect of stratum | |
| Astuti et al. | Adaptive Short Time Fourier Transform (STFT) Analysis of seismic electric signal (SES): A comparison of Hamming and rectangular window | |
| CN102636807B (en) | Electromagnetic-type vibroseis seismic signal detection method | |
| CN116520419B (en) | Hot fluid crack channel identification method | |
| CN102879825A (en) | High power pulse noise detection and suppression method for vibroseis seismic data | |
| Taipodia et al. | Influence of data acquisition and signal preprocessing parameters on the resolution of dispersion image from active MASW survey | |
| Iqbal et al. | Iterative interferometry-based method for picking microseismic events | |
| CN100412569C (en) | Method for conducting earthquake signal high frequency compensation utilizing earthquake micro metering | |
| CN105301634A (en) | Method for increasing correlation quality of controlled seismic source single shot by utilizing record of detector near to seismic source | |
| CN106019377A (en) | Two-dimensional seismic exploration noise removing method based on time-space-domain frequency reduction model | |
| CN114623921B (en) | Demodulation method for distributed optical fiber vibration sensing | |
| Cheng et al. | A new method for estimating the correlation of seismic waveforms based on the NTFT | |
| CN109164490A (en) | A method of improving the weak catoptric imaging precision of desert area | |
| Qiao et al. | Ground Penetrating Radar Weak Signals Denoising via Semi-soft Threshold Empirical Wavelet Transform. | |
| CN102073065B (en) | Method for eliminating single-frequency interference of earthquake data | |
| Du et al. | Research and application of Rayleigh wave extraction method based on microtremors signal analysis | |
| CN111691876A (en) | Method and device for imaging adjacent well by using acoustic logging and storage medium | |
| Han et al. | A three-dimensional microseismic downhole noise suppression based on polarization filtering method in Shearlet transform | |
| Gui et al. | An Efficient Wavelet Transform for Hierarchical Echo Data Processing of Ground-penetrating Radar | |
| Tiwari et al. | Multifractal analysis for seismic wave in Kathmandu valley after Gorkha Earthquake-2015, Nepal |
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 |