CN1673775A - Method for conducting earthquake signal high frequency compensation utilizing earthquake micro metering - Google Patents

Method for conducting earthquake signal high frequency compensation utilizing earthquake micro metering Download PDF

Info

Publication number
CN1673775A
CN1673775A CN 200510049340 CN200510049340A CN1673775A CN 1673775 A CN1673775 A CN 1673775A CN 200510049340 CN200510049340 CN 200510049340 CN 200510049340 A CN200510049340 A CN 200510049340A CN 1673775 A CN1673775 A CN 1673775A
Authority
CN
China
Prior art keywords
seismic
earthquake
depth
signal
frequency compensation
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
CN 200510049340
Other languages
Chinese (zh)
Other versions
CN100412569C (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CNB2005100493404A priority Critical patent/CN100412569C/en
Publication of CN1673775A publication Critical patent/CN1673775A/en
Application granted granted Critical
Publication of CN100412569C publication Critical patent/CN100412569C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

The present invention discloses HF compensating method for seismic signal with seismic microlog. The seismic waveform information is utilized in analyzing spectral characteristic of primary waveform in direct wave and finding out the attenuation features of seismic signal in different depths. Meanwhile, system input and output are used in analyzing absorbing attenuation response of near surface system and performing HF compensation of seismic data. The present invention applies seismic microlog technology in the research of attenuation law of near surface medium, and the waveform characteristic of serial seismic primary waves is used in searching the waveform spectrum variation in different depths, so that one reverse filter factor is found out for the HF compensation of superposed seismic data with artificial seism reflecting signal. The present invention has raised seismic signal distinguishing capacity.

Description

Utilize the earthquake micro logging to carry out seismic signal high-frequency compensation method
Technical field
The present invention relates to a kind of earthquake micro logging that utilizes and carry out seismic signal high-frequency compensation method, particularly utilize the attenuation law of micro logging observation signal to carry out high-frequency compensation to improve the effective way of seismic signal resolution.
Background technology
Earthquake micro logging technology mainly is used to investigate in the past the near surface geologic structure, it mainly is to utilize the different focus of a series of degree of depth to understand the degree of depth of the thickness and the water table of underground low velocity layer time of arrival to the differently seismic wave first arrival of surface geophone, and what it was paid close attention to is the time of arrival of earthquake primary wave.
The attenuation by absorption of near surface compensation often adopts the mode of Q inverse filtering to carry out in the past, but its effect is because the measuring accuracy of Q parameter limits is affected.
Summary of the invention
The purpose of this invention is to provide a kind of method of utilizing the earthquake micro logging to carry out the seismic signal high-frequency compensation.
It is to the resulting micro logging record of different depth focus, the editing and processing that keeps primary wave, its spectral amplitude and phase spectrum are analyzed, try to achieve of the system filter device response of this degree of depth with respect to primary wave on the minimum-depth, system filter device with the minimum-depth primary wave responds as desired output, filter response for the different depth primary wave, obtain the inverse filter that it absorbs earth-attenuation according to the principle of least square inverse filtering, after acting on post-stack seismic data, thereby realized eliminating the high-frequency compensation that near surface is decayed.
The present invention utilizes earthquake micro logging technology to study the attenuation law of near surface medium, it mainly is that the waveform frequency spectrum of utilizing the waveform character of a series of different earthquake primary waves to understand on the different depth changes, comprise amplitude Pu and phase spectrum, be equivalent to the near surface stratum as a system, the response of under the situation of known system input and output, obtaining system.Thereby obtain the high-frequency compensation that an inverse filtering factor pair post-stack seismic data carries out the artificial earthquake reflected signal.The advantage signal to noise ratio (S/N ratio) bandwidth that compensates later reflection wave seismic signal has obtained significantly widening, thereby has improved the resolution characteristic of seismic signal, for lithology seismic prospecting oil gas provides favourable instrument.
Description of drawings
Fig. 1 (a) is the individual well observed pattern synoptic diagram of earthquake micro logging;
Fig. 1 (b) is the twin-well observed pattern synoptic diagram of earthquake micro logging;
Fig. 2 is that the somewhere micro logging is resulting a series of through the later different depth seismic event first arrival signal waveforms of editing and processing;
Fig. 3 is the power spectrum of the corresponding direct wave of the different source depth in somewhere, and among the figure: a left side is 13 meters of the degree of depth, and the right side is 3 meters for the degree of depth;
Fig. 4 is the comparing result of earthquake stacked section before and after the somewhere high-frequency compensation is handled, and among the figure: a left side is for before handling, and is right for after handling;
The power spectrum that Fig. 5 handles the front and back superposition of data for the somewhere contrasts, and among the figure: a left side is for before handling, and the right side is after handling;
The frequency sweep result that Fig. 6 handles the front and back stacked section for the somewhere contrasts, and a left side is for before handling, and the right side is after handling.
Embodiment
The present invention keeps the editing and processing of primary wave to the resulting earthquake micro logging of different depth focus (being generally detonator) record, and its spectral amplitude and phase spectrum are analyzed, and tries to achieve the filter response of this degree of depth with respect to the primary wave system of minimum-depth.System filter device with the minimum-depth primary wave responds as desired output, filter response for the different depth primary wave, just can obtain the inverse filter that it absorbs earth-attenuation according to the principle of least square inverse filtering, after all post-stack seismic datas are handled, thus the high-frequency compensation of realization seismic signal.
Open-air earthquake micro logging observation technology of the present invention is provided by following diagram.Fig. 1 (a) is the individual well observed pattern, and Fig. 1 (b) is the twin-well observed pattern.At near surface lithology competitive list for the moment, individual well micro logging observed pattern can be adopted, and during near surface lithology more complicated, two well logging observed patterns should be adopted.
Fig. 2 has represented the signal of a series of micro logging seismic event first arrival, and their waveform has faint difference as can be seen.We utilize this difference just can carry out spectrum analysis to the earthquake primary wave signal on the different depth, thereby obtain their decay absorbing rule on spectral amplitude and phase spectrum.
The micro logging in past is often only noted the time of arrival of primary wave, in order to see the time of primary wave take-off clearly, thereby sometimes the gain of seismic instrument is transferred very high, also be enlargement factor give very big, make waveform subsequent distortion occur, this situation must prevent in the application of present technique in strictness.
Fig. 3 is the result of the seismic event first arrival signal amplitude analysis of spectrum of (13 meters and 3 meters) on two different depths in a certain micro logging, and both spectral amplitudes are different as can be seen, and particularly the radio-frequency component of deep signal has been attenuated significantly as can be seen.
It is as follows that the present invention compensates the inverse filter principle:
Excite if the well depth in micro logging is the position of i, the seismic signal that ground receives is x i(t), be 1 through thickness jThe stratum to another degree of depth j, then Ci Shi seismic signal is
x j(t)=q j(t)*x i(t) (1)
In the formula, q j(t) expression zone thickness 1 jTo the influence of fading of reflected signal, i.e. the response of system.So, inverse operator P j = q j - 1 = { p i ( t ) } Just can compensate stratum 1 jDecay to reflected signal frequency and amplitude has
x i(t)=p j(t)*x j(t) (2)
In time domain, equation (2) is to utilize the Levinson algorithm to obtain p jCan certainly find the solution in frequency field, Fourier transform is done on equation (2) both sides, have
X i(f)=P j(f)×X j(f) (3)
Can get by equation (3)
P j ( f ) = X i ( f ) / X i ( f ) = X i ( f ) × Xj ( f ) ‾ / | Xj ( f ) | 2 - - - ( 4 )
Consider noise, have
P j ( f ) = X i ( f ) / X j ( f ) = X i ( f ) × Xj ( f ) ‾ / [ | Xj ( f ) | 2 + α 2 ] - - - ( 5 )
A represents noise in the formula (5).We are with inverse operator p j(t) be called the absorption compensation inverse filter.When frequency field finally realizes inverse filtering, note the difference that time-sampling rate difference is caused.Our way is the inverse filtering factor length lengthening with time domain, so that both frequency intervals reach consistent.
The compensation result contrast:
As shown in Figure 4, left side figure is the stacked section before the high-frequency compensation, and right figure is the later stacked section of high-frequency compensation.The radio-frequency component showed increased of seismic signal on the section after handling as can be seen, quantity of information obviously increases, and a lot of detail sections are reflection to some extent all, and some interferes complex wave to obtain separation.
Fig. 5 is the spectrum analysis result to seismic prospecting zone of interest among Fig. 4, it has illustrated the high-frequency compensation treatment effect especially quantitatively, left side figure is the spectrum analysis result of zone of interest on the preceding seismic section of high-frequency compensation, and right figure is the spectrum analysis result of zone of interest on the later seismic section of high-frequency compensation.Advantage signal to noise ratio (S/N ratio) bandwidth before high-frequency compensation will be handled is as can be seen widened to 0~160 hertz from 0~70 hertz, the high band energy has obtained remarkable reinforcement, the basic simultaneously low-frequency component that keeps does not have to change, and so just makes the resolution characteristic of seismic signal be greatly improved.
Fig. 6 is to the frequency sweep analysis of two seismic sections among Fig. 4 (80~150 hertz) comparison diagram as a result, as can be seen in 80~150 hertz of frequency ranges, seismic section before the compensation deals has not been seen reflection line-ups, and the still visible later on a plurality of seismic reflection lineups of high-frequency compensation.

Claims (1)

1. one kind is utilized the earthquake micro logging to carry out seismic signal high-frequency compensation method, it is characterized in that, to the resulting micro logging record of different depth focus, the editing and processing that keeps primary wave, its spectral amplitude and phase spectrum are analyzed, try to achieve of the system filter device response of this degree of depth with respect to primary wave on the minimum-depth, system filter device with the minimum-depth primary wave responds as desired output, filter response for the different depth primary wave, obtain the inverse filter that it absorbs earth-attenuation according to the principle of least square inverse filtering, after acting on post-stack seismic data, thereby realized eliminating the high-frequency compensation that near surface is decayed.
CNB2005100493404A 2005-03-11 2005-03-11 Method for conducting earthquake signal high frequency compensation utilizing earthquake micro metering Expired - Fee Related CN100412569C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2005100493404A CN100412569C (en) 2005-03-11 2005-03-11 Method for conducting earthquake signal high frequency compensation utilizing earthquake micro metering

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2005100493404A CN100412569C (en) 2005-03-11 2005-03-11 Method for conducting earthquake signal high frequency compensation utilizing earthquake micro metering

Publications (2)

Publication Number Publication Date
CN1673775A true CN1673775A (en) 2005-09-28
CN100412569C CN100412569C (en) 2008-08-20

Family

ID=35046434

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2005100493404A Expired - Fee Related CN100412569C (en) 2005-03-11 2005-03-11 Method for conducting earthquake signal high frequency compensation utilizing earthquake micro metering

Country Status (1)

Country Link
CN (1) CN100412569C (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101929332A (en) * 2009-06-26 2010-12-29 中国石油集团东方地球物理勘探有限责任公司 Method for determining speed and buried depth of undrilled stratum by adopting long-array micrometer well
CN102323616A (en) * 2011-06-08 2012-01-18 浙江大学 Frequency division matching method for increasing seismic data resolution of limestone exposure area
CN102109612B (en) * 2009-12-23 2012-10-17 中国石油天然气集团公司 Seismic wave absorption and attenuation compensation method
CN111722275A (en) * 2019-03-21 2020-09-29 中石化石油工程技术服务有限公司 Broadband scanning signal design method based on absorption attenuation compensation
CN113552616A (en) * 2020-04-24 2021-10-26 江苏省地震局 Detector data compensation method and system

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102269822B (en) * 2010-06-02 2013-07-31 中国石油天然气集团公司 Mixed stratigraphic absorption compensation method

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU428331A1 (en) * 1971-07-14 1974-05-15 Институт геофизики DEVICE FOR ACOUSTIC WAGGES
US4168483A (en) * 1977-09-06 1979-09-18 The United States Of America As Represented By The Administrator Of The National Aeronautics & Space Administration System for detecting substructure microfractures and method therefor
US4319346A (en) * 1980-04-21 1982-03-09 Halliburton Company Acoustic frequency spectrum logging

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101929332A (en) * 2009-06-26 2010-12-29 中国石油集团东方地球物理勘探有限责任公司 Method for determining speed and buried depth of undrilled stratum by adopting long-array micrometer well
CN101929332B (en) * 2009-06-26 2012-12-12 中国石油集团东方地球物理勘探有限责任公司 Method for determining speed and buried depth of undrilled stratum by adopting long-array micrometer well
CN102109612B (en) * 2009-12-23 2012-10-17 中国石油天然气集团公司 Seismic wave absorption and attenuation compensation method
CN102323616A (en) * 2011-06-08 2012-01-18 浙江大学 Frequency division matching method for increasing seismic data resolution of limestone exposure area
CN111722275A (en) * 2019-03-21 2020-09-29 中石化石油工程技术服务有限公司 Broadband scanning signal design method based on absorption attenuation compensation
CN111722275B (en) * 2019-03-21 2023-03-10 中国石油化工集团有限公司 Broadband scanning signal design method based on absorption attenuation compensation
CN113552616A (en) * 2020-04-24 2021-10-26 江苏省地震局 Detector data compensation method and system
CN113552616B (en) * 2020-04-24 2024-03-15 江苏省地震局 Method and system for compensating detector data

Also Published As

Publication number Publication date
CN100412569C (en) 2008-08-20

Similar Documents

Publication Publication Date Title
Lu et al. Seismic spectral decomposition using deconvolutive short-time Fourier transform spectrogram
Xue et al. Application of the empirical mode decomposition and wavelet transform to seismic reflection frequency attenuation analysis
CN106896409B (en) A kind of varying depth cable ghost reflection drawing method based on wave equation boundary values inverting
CN1673775A (en) Method for conducting earthquake signal high frequency compensation utilizing earthquake micro metering
CN106526678A (en) Reflection acoustic logging wave field separation method and device
US20100286922A1 (en) Method for detecting and/or processing seismic signals
CN102998703B (en) Method and device for conducting reservoir prediction and based on earth surface consistency deconvolution
CN106483564B (en) A method of fluid identification is carried out using earthquake low-frequency information
Banjade et al. Earthquake accelerogram denoising by wavelet-based variational mode decomposition
CN104391324A (en) Seismic trace set dynamic correction stretching correction pre-processing technology before AVO inversion depending on frequency
CN103135133A (en) Vector noise reduction method and device for multi-component seismic data
Li et al. Wigner-Ville distribution and its application in seismic attenuation estimation
CN106324702B (en) A kind of quantitative evaluation method of seismic interference method imaging observation system design
CN107561578B (en) A kind of surface wave pressing method for wide-azimuth seismic data
Kuhn A numerical study of Lamb's problem
Ebadi Coherent and incoherent seismic noise attenuation using parabolic radon transform and its application in environmental geophysics
CN105572742B (en) It is a kind of to determine sea water advanced method and apparatus
Zhang et al. Data processing of a wide-azimuth, broadband, high-density 3D seismic survey using a low-frequency vibroseis: a case study from Northeast China
Baradello et al. Vibroseis deconvolution: A comparison of pre and post correlation vibroseis deconvolution data in real noisy data
Haney et al. Insight into the output of reverse-time migration: What do the amplitudes mean?
CN107024716A (en) A kind of seismic wave field absorption compensation imaging method and system
Dricker et al. Statistically optimal technique of simultaneous event location and focal mechanism determination of weak microseismicity using surface arrays
Du et al. Research and application of Rayleigh wave extraction method based on microtremors signal analysis
Gui et al. An Efficient Wavelet Transform for Hierarchical Echo Data Processing of Ground-penetrating Radar
Aniwetalu et al. Spectral analysis of Rayleigh waves in south-eastern parts of Niger delta, Nigeria

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20080820

Termination date: 20120311