CN105093275B - Useful signal based on rate pattern extracts ISS wave field separation methods - Google Patents

Useful signal based on rate pattern extracts ISS wave field separation methods Download PDF

Info

Publication number
CN105093275B
CN105093275B CN201410195660.XA CN201410195660A CN105093275B CN 105093275 B CN105093275 B CN 105093275B CN 201410195660 A CN201410195660 A CN 201410195660A CN 105093275 B CN105093275 B CN 105093275B
Authority
CN
China
Prior art keywords
data
wave field
useful signal
signal
time
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.)
Active
Application number
CN201410195660.XA
Other languages
Chinese (zh)
Other versions
CN105093275A (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
Sinopec Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
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, Sinopec Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201410195660.XA priority Critical patent/CN105093275B/en
Publication of CN105093275A publication Critical patent/CN105093275A/en
Application granted granted Critical
Publication of CN105093275B publication Critical patent/CN105093275B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

ISS wave field separation methods are extracted the invention provides a kind of useful signal based on rate pattern, belong to oil gas (energy) Exploration Domain.This method includes:(1) CDP trace gathers are generated:The data of many focus simultaneous shots are converted into common CDP trace gathers data;(2) speed spectrum analysis and t0Time picks up:Normal-moveout spectrum is done to the common CDP trace gathers data, and analysis is interacted to normal-moveout spectrum, the t of each reflecting layer position is tentatively picked up0Time;(3) the effective continuous signal of pickup is followed the trail of, useful signal is obtained:(4) other noise signals are removed:It will not be that other all information of useful signal regard noise as, and these noises directly subtracted from seismic channel data, there is no the signal for interfering big gun noise obtained from;(5) seismic channel data inverse transformation:Inverse transformation is carried out to the common CDP trace gathers data after step (4) processing, the shot gather data body after wave field separation is obtained.

Description

Useful signal based on rate pattern extracts ISS wave field separation methods
Technical field
The invention belongs to oil gas (energy) Exploration Domain, and in particular to a kind of useful signal based on rate pattern is extracted ISS wave field separation methods, the early stage for many focus simultaneous shots (controlled source synchronous scanning) seismic data in oil-gas exploration Processing.
Background technology
Many focus simultaneous shot seismic acquisition technologies are a kind of efficient seismic data acquisition technologies.It can be using routinely Shake data collection technology and gather time of a big gun data and gather two big guns or many big gun data, therefore can greatly improve collecting efficiency, Reduce acquisition cost.Seismic acquisition technology can be divided into a variety of, and independent synchronization scanning technique (ISS) is wherein efficiency highest Acquisition technique.But because the same time there are multiple focus exciting, cause wave field in earthquake record to interfere, data processing is One difficult point of the technology.
The treatment technology route for many vibroseis acquisition seismic datas that different methods is gathered is each variant, and ISS is provided The processing of material mainly has two technology paths, and one, be according to the characteristics of many focus seismic datas, to develop special data processing skill Art;2nd, it is to use wave field separation technology, many focus seismic datas is separated into single focus seismic data, then using existing Common seismic data treatment technology is handled.Need to develop a series of technology for former approach, at present still in grinding Study carefully the stage.Many focus seism processings at this stage mainly use later approach.
Wave field separation for ISS data is mainly randomness according to simultaneous shot focus on firing time, by ground Shake data are transformed in non-big gun trace gather (for example:Common midpoint gather, common offset trace gather or common detector gather etc.), now, Useful signal shows as being relevant, and interferes big gun noise to be random, using this species diversity, is made an uproar at random using various removals The method of sound can just eliminate interference big gun noise, so as to reach the effect of wave field separation.But the various sides for removing random noise Method all has the damage to useful signal, and this damage is often very fatal for lithology and reservoir exploration. In order to avoid this damage, when selecting denoising parameter, often more from protection useful signal so that noise is filtered Effect is undesirable, ultimately results in wave field separation undesirable.And how to realize many focus seismic datas high-purity separation be realize The key of its data high accuracy imaging.
The content of the invention
It is an object of the invention to solve problem present in above-mentioned prior art, there is provided a kind of having based on rate pattern Signal extraction ISS wave field separation methods are imitated, the high-purity separation to many focus seismic datas of ISS is realized, so that existing to utilize Common seismic data treatment technology the high Precision Processing of many focus seismic datas such as ISS is laid the foundation.
The present invention is achieved by the following technical solutions:
A kind of useful signal based on rate pattern extracts ISS wave field separation methods, including:
(1) CDP trace gathers are generated:The data of many focus simultaneous shots are converted into common CDP trace gathers data;
(2) speed spectrum analysis and t0Time picks up:The common CDP trace gathers data are done with normal-moveout spectrum, and normal-moveout spectrum is carried out Interaction analysis, tentatively picks up the t of each reflecting layer position0Time;
(3) the effective continuous signal of pickup is followed the trail of, useful signal is obtained:
(4) other noise signals are removed:To not be that other all information of useful signal regard noise as, and from seismic channel number These noises are directly subtracted in, do not have the signal (signal that namely we need) for interfering big gun noise obtained from;
(5) seismic channel data inverse transformation:Inverse transformation is carried out to the common CDP trace gathers data after step (4) processing, obtained Shot gather data body after wave field separation.
What the step (3) was realized in:
According to t0Time opens a hour window, using the uniformity in zero-offset road and its neighbouring a few road signals, using phase Algorithm is closed, the time corresponding to zero-offset road peak swing is calculated in hour window, and as the t of modification0Time, then According to the uniformity of useful signal, using related algorithm calculate the maximum data of adjacent track data zero-offset road correlation with The time difference tn in zero-offset road, and by its t with modification0Time is added, so as to obtain the travelling of adjacent one of reflected signal When, the reflected signal of the seismic channel is thus obtained, by that analogy, the reflected signal of each seismic channel is obtained.
In above-mentioned steps, step (2) speed spectrum analysis and t0In time pickup, it is desirable to take t0Time, just it is necessary to have speed Model is spent, and rate pattern is obtained by speed spectrum analysis, it can be seen that the present invention is the wave field separation based on rate pattern Method.
Compared with prior art, the beneficial effects of the invention are as follows:This method is used to gather the ISS that 2D and 3D is simulated and provided Material carries out wave field separation and achieves ideal effect, while this Weak Signal Extraction is when certain in window, Only it is tracked according to the continuity of useful signal, and to the feature of noise signal without being distinctly claimed, therefore this method is not only Suitable for the wave field separation to ISS data, for other kinds of highly effective gathering seismic data (such as:D4S wave field separation) It is effective.
Brief description of the drawings
Fig. 1 is the step block diagram of the inventive method.
Fig. 2-1 is the big gun record before 2D simulation ISS data wave field separations, and interference big gun noise is developed very much.
Fig. 2-2 is the big gun record after 2D simulation ISS data wave field separations, and interference big gun noise is substantially eliminated.
Fig. 3-1 is the big gun record before 3D simulation ISS data wave field separations, and interference big gun noise is developed very much.
Fig. 3-2 is the big gun record after 3D simulation ISS data wave field separations, and interference big gun noise is substantially eliminated.
Fig. 4-1 is the big gun record before field D4S collection data wave field separations, well-regulated interference big gun noise.
Fig. 4-2 is the big gun record after field D4S collection data wave field separations, and interference big gun noise is substantially eliminated.
Embodiment
The present invention is described in further detail below in conjunction with the accompanying drawings:
In the denoising of seismic data, the denoising of special low SNR data is a technological difficulties.Occur in recent years A kind of weak signal extracts noise-removed technology.Seismic data is transformed to CDP trace gathers by the technology first, and normal-moveout spectrum point is carried out to it Analysis, according to each reflecting layer position t in speed spectrum analysis data0Time opens the hour window of a certain length, then gathers related skill Art, with reference to the stack velocity at each layer position interface, is tracked according to the continuity of seismic signal between neighboring track to useful signal Pickup, it is other to be removed as noise.This denoising method avoids damage of the conventional noise-removed technology to useful signal, reaches More preferably protect the effect of useful signal.
Many focus seismic acquisition technologies gather data because the same time often has the excitation signal of multiple focus mutual Interfere and cause signal to noise ratio very low, therefore can be regarded as a kind of special low SNR data, can be extracted and gone with this weak signal Technology is made an uproar to eliminate interference big gun noise, so that the effect of wave field separation is reached, and also this noise-removed technology is lossless to useful signal Wound, therefore be a kind of high-precision ISS wave field separations technology.
Useful signal of the present invention based on rate pattern extracts the step block diagram of ISS wave field separation methods as shown in figure 1, bag Include:
1) CDP trace gathers are generated
The data of many focus simultaneous shots are exactly converted into common CDP trace gathers data.
2) speed spectrum analysis and t0Time picks up
Using existing business conventional treatment software, (processing system of seismic data is a lot, and all systems have this Function, such as:FOCUS processing systems, OMEGA processing systems etc.) normal-moveout spectrum is done to common CDP trace gathers, and normal-moveout spectrum data is carried out Interaction analysis, tentatively picks up the t of each reflecting layer position0Time.
3) effectively continuous signal follows the trail of pickup
According to t0Time opens a hour window, using the similitude in zero-offset road and its neighbouring a few road signals, using phase Stem algorithm, changes t0Time, coherent algorithm is recycled to pick up useful signal according to the continuity of useful signal.
4) other noise signals are removed
To not be that other all information of useful signal are regarded noise as and filtered out.
5) seismic channel data inverse transformation
Inverse transformation is carried out to common CDP trace gathers data, so that the shot gather data body formed after wave field separation.
One embodiment of this method is as follows:
This method is used for the wave field separation that ISS data is simulated to 2D first, before Fig. 2-1 and Fig. 2-2 are wave field separation The contrast of shot gather data afterwards.Interference big gun noise energy is very strong it can be seen from big gun record figure before its wave field separation, and Develop very much, it is difficult resolution to cause useful signal.After wave field separation, interference big gun noise is substantially completely eliminated, and illustrates this Method is a kind of precision very high wave field separation technology.Secondly this method is used for the wave field separation that ISS data is simulated to 3D, Fig. 3-1 and Fig. 3-2 is the contrast of shot gather data before and after 3D simulation ISS data wave field separations.Thus figure is it can be seen that this method Separation accuracy it is very high.The method of this invention is finally applied to the wave field separation of D4S field acquisition data, such as Fig. 4-1 With shown in Fig. 4-2.Thus it is also seen that the separating effect of the inventive method.
The wave field separation of many focus seismic datas is to realize to complete to many focus seismic datas using existing business software One of key technology of reason, is that a kind of weak signal extractive technique based on stack velocity model is used for wave field separation by this, so that The useful signal for being developed based on rate pattern extracts ISS wave field separation technologies, i.e., utilize normal-moveout spectrum solution to concentrating in common CDP The result released, useful signal is extracted in the capable and experienced hybrid recording for relating to big gun noise, so that elimination interference big gun noise is reached, it is real The purpose of existing wave field separation.
Above-mentioned technical proposal is one embodiment of the present invention, for those skilled in the art, at this On the basis of disclosure of the invention application process and principle, it is easy to make various types of improvement or deformation, this is not limited solely to Invent the method described by above-mentioned embodiment, therefore previously described mode is preferred, and and without limitation The meaning of property.

Claims (2)

1. a kind of useful signal based on rate pattern extracts ISS wave field separation methods, it is characterised in that:Methods described includes:
(1) CDP trace gathers are generated:The data of many focus simultaneous shots are converted into common CDP trace gathers data;
(2) speed spectrum analysis and t0Time picks up:Normal-moveout spectrum is done to the common CDP trace gathers data, and normal-moveout spectrum interacted point Analysis, tentatively picks up the t of each reflecting layer position0Time;
(3) the effective continuous signal of pickup is followed the trail of, useful signal is obtained:
(4) other noise signals are removed:To not be that other all information of useful signal regard noise as, and from seismic channel data These noises are directly subtracted, do not have the signal for interfering big gun noise obtained from;
(5) seismic channel data inverse transformation:Inverse transformation is carried out to the common CDP trace gathers data after step (4) processing, wave field is obtained Shot gather data body after separation.
2. the useful signal according to claim 1 based on rate pattern extracts ISS wave field separation methods, its feature exists In:What the step (3) was realized in:
According to t0Time opens a hour window, using the uniformity in zero-offset road and its neighbouring a few road signals, using correlation calculation Method, calculates the time corresponding to zero-offset road peak swing in hour window, and as the t of modification0Time, further according to The uniformity of useful signal, it is inclined using the maximum data of related algorithm calculating adjacent track data zero-offset road correlation and zero The time difference tn in Yi Ju roads, and by its t with modification0Time is added, so that when obtaining the travelling of adjacent one of reflected signal, by This obtains the reflected signal of the seismic channel, by that analogy, obtains the reflected signal of each seismic channel.
CN201410195660.XA 2014-05-09 2014-05-09 Useful signal based on rate pattern extracts ISS wave field separation methods Active CN105093275B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410195660.XA CN105093275B (en) 2014-05-09 2014-05-09 Useful signal based on rate pattern extracts ISS wave field separation methods

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410195660.XA CN105093275B (en) 2014-05-09 2014-05-09 Useful signal based on rate pattern extracts ISS wave field separation methods

Publications (2)

Publication Number Publication Date
CN105093275A CN105093275A (en) 2015-11-25
CN105093275B true CN105093275B (en) 2017-08-18

Family

ID=54574129

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410195660.XA Active CN105093275B (en) 2014-05-09 2014-05-09 Useful signal based on rate pattern extracts ISS wave field separation methods

Country Status (1)

Country Link
CN (1) CN105093275B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113552632B (en) * 2020-04-23 2022-11-01 中国石油大学(北京) Seismic cross-correlation signal pickup method and system based on wavelet domain Kalman filtering
CN111474584B (en) * 2020-05-29 2023-05-05 核工业北京地质研究院 Focusing superposition imaging method and system based on correlated seismic interference

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101893719A (en) * 2010-04-16 2010-11-24 西安石油大学 Method for separating and de-noising wave fields of longitudinal and transverse earthquake waves
CN103675903A (en) * 2012-09-07 2014-03-26 中国石油化工股份有限公司 Multi-source seismic wave field separating method by utilizing random denoising technology

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6738715B2 (en) * 2001-09-14 2004-05-18 Exxonmobil Upstream Research Company Method for attenuating noise in seismic data

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101893719A (en) * 2010-04-16 2010-11-24 西安石油大学 Method for separating and de-noising wave fields of longitudinal and transverse earthquake waves
CN103675903A (en) * 2012-09-07 2014-03-26 中国石油化工股份有限公司 Multi-source seismic wave field separating method by utilizing random denoising technology

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
FK滤波在三维VSP波场分离中的应用;吴永栓 等;《中国地球物理2009》;20091231;第470页 *
地震勘探波场分离与去噪处理思路;沈鸿雁;《西部探矿工程》;20121231(第3期);第87-90,93页 *

Also Published As

Publication number Publication date
CN105093275A (en) 2015-11-25

Similar Documents

Publication Publication Date Title
CN101630017B (en) Method for separating seismic wave fields of different types in two-dimensional vertical seismic profile
US9766357B2 (en) Seismic image dip decomposition estimation and recomposition
CN103926623B (en) Method for suppressing reverse time migration low frequency noise
CN106154319B (en) A kind of method for separating for being imaged trace gather
CN106154330B (en) Electromagnetic type controlled source parabolic Radon transform suppresses parallel focus crosstalk noise method
CN103675903A (en) Multi-source seismic wave field separating method by utilizing random denoising technology
CN104620132A (en) Signal enhancement using diversity shot stacking for reverse time migration (DeSSeRT)
CN105911585A (en) Method and device for extracting seismic record regular interference waves
CN110389377B (en) Microseism offset imaging positioning method based on waveform cross-correlation coefficient multiplication
CN107807393A (en) Separate unit station collection preliminary wave Enhancement Method based on seismic interference method
CN105319587A (en) Vibroseis shaping filtering method based on first arrival wavelet extraction
CN105629300B (en) The method for improving complicated structure offset data signal-to-noise ratio
CN104570116A (en) Geological marker bed-based time difference analyzing and correcting method
CN105093275B (en) Useful signal based on rate pattern extracts ISS wave field separation methods
CN104749623B (en) A kind of imaging of seismic data processing method
CN109581481B (en) Portable high-frequency controllable seismic source seismic signal harmonic interference elimination method
CN103076626A (en) Wave field purification treatment method
AU2006237348A1 (en) Method of processing seismic data for AVO or AVOA characterisation
CN105301634A (en) Method for increasing correlation quality of controlled seismic source single shot by utilizing record of detector near to seismic source
AU2015200555B2 (en) Correction of sea surface state
CN105093282A (en) Energy replacement surface wave suppressing method based on frequency constraint
CN110737022B (en) Suppression method for vibroseis to excite noise of seismic data black triangle area
CN110780346A (en) Separation method for advanced detection of complex seismic wave field in tunnel
CN105425300A (en) Residual static correction method
CN104536048A (en) Method for automatically cutting seismic wave arrival

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant