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 PDF

Info

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
Application number
CN201910018313.2A
Other languages
Chinese (zh)
Other versions
CN109581481B (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.)
East China Institute of Technology
Original Assignee
East China Institute of Technology
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 East China Institute of Technology filed Critical East China Institute of Technology
Priority to CN201910018313.2A priority Critical patent/CN109581481B/en
Publication of CN109581481A publication Critical patent/CN109581481A/en
Application granted granted Critical
Publication of CN109581481B publication Critical patent/CN109581481B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/16Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/30Noise handling
    • G01V2210/32Noise 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

A kind of portable high frequency controlled source seismic signal harmonic wave interference removing method
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.
CN201910018313.2A 2019-01-09 2019-01-09 Portable high-frequency controllable seismic source seismic signal harmonic interference elimination method Active CN109581481B (en)

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)

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

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

Patent Citations (6)

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

* Cited by examiner, † Cited by third party
Title
杨志超: ""电磁式可控震源地震信号检测关键技术研究"", 《中国博士学位论文全文数据库•基础科学辑》 *

Cited By (4)

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