CN102798894B - Sliding scanning harmonic suppression method based on sparse inversion - Google Patents
Sliding scanning harmonic suppression method based on sparse inversion Download PDFInfo
- Publication number
- CN102798894B CN102798894B CN201110140361.2A CN201110140361A CN102798894B CN 102798894 B CN102798894 B CN 102798894B CN 201110140361 A CN201110140361 A CN 201110140361A CN 102798894 B CN102798894 B CN 102798894B
- Authority
- CN
- China
- Prior art keywords
- record
- harmonic
- relevant
- arrival
- strong
- 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
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a sliding scanning harmonic suppression method based on sparse inversion. Scanning of subsequent scanning in a sliding scanning continuous record after hearing time is considered as a record which is not polluted by the harmonic and serves as a starting shot for harmonic suppression processing; and a relevant after-normal single-shot record after the harmonic of a front shot is suppressed is obtained through inversion processing along the actual acquisition time. By the method, the main harmonic pollution components can be extracted and obviously suppressed; compared with the quality of the conventional single shot, the quality of the suppressed single-shot record is obviously improved, and the strong high-frequency harmonic pollution of a middle deep layer is eliminated.
Description
Technical field
The present invention relates to controlled source slip scan acquisition technique field in oil exploration, is a kind of slip scan harmonic suppression method based on sparse inversion for removing the harmonic pollution existed in controlled source slip scan image data.
Background technology
Adopt controlled source when gathering by land, due to reasons such as controlled source and the Non-linear coupling on earth's surface, the non-linear factor of controlled source itself and insufficient feedback systems, often can observe force signal and receiving record all exists harmonic distortion.If adopt down sweep mode, the fundametal compoment of force signal and receiving record cross-correlation, what humorous wave interference was present in record continues to position, and the humorous wave interference of first arrival so by force often affects the identification of back wave; And adopting upsweep mode, the fundametal compoment of force signal and receiving record cross-correlation, humorous wave interference is present in the front end of record, and the humorous wave interference of so strong first arrival would not affect back wave.When therefore constructing, often adopt raising frequency linear scanning mode in the wild.
The people such as A.J.Seriff in 1970 propose the essence & characteristics of harmonic wave in slip scan record and the trial of various compacting harmonic wave, and wherein major part processes for gross data.Slip scan collection is exactly before the data record end of last group of controlled source scanning, scanning is started after the time interval (sliding time) that rear one group of ready controlled source sets at instrument, wherein sliding time can not be less than record length, instrument processes the original mother's record received, first to be syncopated as the original record of every big gun according to the record length of every big gun according to the TB time, then relevant to the original record of each big gun by reference signal, export the single shot record after being correlated with.Therefore, it is that the controlled source of group more than a kind of utilization carries out continuously or the acquisition technique of overlapping scan, improves collecting efficiency because it can shorten the acquisition time of single shake frequency.
But, controlled source conventional Christmas method can weaken the impact of harmonic pollution to a great extent by upsweep, but the harmonic pollution impact that follow up scan produces front one scan then must be faced when the operation of controlled source slip scan, be forced in seeking balance between collecting efficiency and data quality---extend sliding time on the one hand, reduce the impact of humorous wave interference, but reduce collecting efficiency; Then tolerate harmonic pollution to a certain degree on the other hand, can production efficiency be significantly improved, but harmonic pollution makes data quality be affected.Thus, the bottleneck problem affecting its extensive use is become.
Summary of the invention
The object of the invention is to carry and a kind ofly solves the harmonic pollution problems existed in controlled source slip scan record, improves the slip scan harmonic suppression method based on sparse inversion of data quality.
The present invention is realized by following concrete steps:
1) controlled source acquiring seismic data, data and force signal before slip scan is relevant, first process from last big gun, short time discrete Fourier transform is utilized to be separated at time-frequency domain to the force signal of last big gun, isolate first-harmonic, second harmonic, triple-frequency harmonics composition, by isolated fundametal compoment to relevant before record and makes cross-correlation, record x after obtaining being correlated with
0(t);
2) for one, take absolute value maximum x from relevant rear record
max, search from sequence | x
0(t) | > α x
max, formed with x
maxcentered by one section of region, determine the strong first arrival of controlled source;
Wherein: x
maxfor maximum absolute value, α is the constant of < 1;
3) utilize isolated fundametal compoment from force signal to carry out auto-correlation and obtain zero-phase wavelet, utilize this wavelet to carry out match tracing in relevant rear record, obtain the rarefaction reflectance factor of first arrival, obtain the fundametal compoment recording strong first arrival before being correlated with;
4) fundametal compoment deducting strong first arrival from relevant front record obtains remaining record
by 1 order harmonics s of itself and force signal
1t () carries out cross-correlation, obtain relevant after record x
1(t);
5) after relevant, x is recorded
1(t), the wavelet w obtained by 1 order harmonics auto-correlation
1t () adopts match tracing method in strong first arrival region, obtain the first arrival coefficient of an order harmonics, thus obtain strong first arrival record before correlation in 1 order harmonic component;
6) from step 4) the residue record that obtains
in deduct step 5) strong first arrival one order harmonic component that obtains, obtain another residue record, according to from step 4) to step 5) extract the method for an order harmonic component, extract strong first arrival record before correlation in the more high-order harmonic wave component such as 2 rank, 3 rank;
7) adopt step 3) obtain relevant before record fundametal compoment and the step 5 of strong first arrival) and step 6) obtain 1 rank, the 2 contour order harmonic component in rank, according to the sliding time gathered, from forward gun, deduct these harmonic components in relevant position, obtain the relevant front record not having harmonic pollution;
8) step 7) obtain do not have harmonic pollution relevant before record and the fundametal compoment that extracts from forward gun force signal carry out cross-correlation, normal single shot record after just to obtain after forward gun compacting harmonic wave relevant.
The present invention suppresses the degree very nearly the same without Harmonic Theory data reaching and simulate by harmonic wave, can well main harmonic pollution constituents extraction out, and well mute, from time frequency analysis result, the harmonic components of 1 time, the 2 times high frequencies within 5 times is all by obviously compacting (shown in Fig. 2).Single shot record (shown in Fig. 4) after compression process significantly improves with original single big gun specific mass quality, and the strong high-frequency harmonic eliminating mid-deep strata pollutes, and these harmonic pollutions all concentrate near migration range scope (as shown in Figure 5).
Accompanying drawing explanation
Fig. 1 is last Bao Mou road time-frequency record, the time-domain record in Shi Gai road above, and here is corresponding time frequency analysis;
Fig. 2 is last big gun detuning Bo Houmou road time-frequency, removing only the harmonic wave within 5 times herein, the time-domain record in Shi Gai road above, and here is corresponding time frequency analysis;
Fig. 3 removes humorous wavefront single shot record, the harmonic pollution of the visible high frequency of near migration range mid-deep strata;
Fig. 4 is single shot record after removal harmonic wave, within removing 5 times after harmonic pollution;
Fig. 5 is the harmonic wave display of removing, the display of harmonic wave on single big gun within 5 times that remove.
Detailed description of the invention
The present invention is described in detail below in conjunction with accompanying drawing.
Harmonic noise in slip scan record is rendered as the pollution of follow up scan to preceding scans record, therefore, in slip scan can being recorded continuously, the scanning of follow up scan after listening the time is thought not by the record of harmonic pollution, as the initial big gun of harmonic wave compression process, along the counter movement process of actual acquisition time.The present invention decomposes the harmonic wave in the non-relative recording of slip scan and removes.
The invention process flow process is as follows:
(1) the controlled source slip scan geological data gathered comprises: slip scan is correlated with front data and force signal corresponding to every big gun.First process from last big gun, force signal is positioned at the force signal of auxiliary channel to last big gun and utilizes short time discrete Fourier transform to be separated at time-frequency domain, isolate first-harmonic, second harmonic, triple-frequency harmonics composition, and isolated fundametal compoment is made cross-correlation to relevant front record, obtain the record x after being correlated with
0(t);
(2) for one, take absolute value maximum x from relevant rear record
max, search from sequence | x
0(t) | > α x
max, determine first arrival scope;
Wherein: x
maxfor maximum absolute value, α is the constant of < 1, can get the size of α, such as 0.9 according to data characteristics; Region after searching for like this defines with x
maxcentered by one section of region, this region is just the strong first arrival region in this Bao Zhe road.
(3) to utilize in the first step isolated fundametal compoment from force signal to carry out auto-correlation and obtain zero-phase wavelet, this wavelet is utilized to carry out match tracing in relevant rear record, obtain the rarefaction reflectance factor of first arrival, concrete grammar is: make cross-correlation with zero-phase wavelet and relevant rear trace record, find relevant maximum point, extract this sample value, and then search for the 2nd relevant maximum point, extract its sample value, so, extract the sample value of series of discrete successively, such sample value sequence is exactly the first arrival reflection coefficient sequence of the first-harmonic of rarefaction.There is reflectance factor, with first-harmonic composition and the reflectance factor convolution of the force signal obtained in the first step, just obtained the fundametal compoment recording strong first arrival before being correlated with;
(4) from relevant front record, deduct the fundametal compoment of the strong first arrival that the 3rd step obtains, obtain remaining record
, will
with 1 order harmonics s of force signal
1t () carries out cross-correlation, obtain relevant after record x
1(t);
(5) after relevant, x is recorded
1t (), carries out auto-correlation with 1 order harmonics obtained in the first step and obtains wavelet w
1t (), utilizes w
1t () records x after relevant
1t () obtains first arrival coefficient of first-harmonic with the 3rd step in strong first arrival region the same, adopts match tracing method, obtains the first arrival coefficient of an order harmonics, thus obtain strong first arrival record before correlation in 1 order harmonic component;
(6) the residue record obtained from the 4th step
in deduct strong first arrival one order harmonic component that (5) step obtains, obtain another residue record, according to from (4)-extract the method for an order harmonic component to (5) step, extract strong first arrival record before correlation in the more high-order harmonic wave component such as 2 rank, 3 rank;
(7) had (3) step to obtain relevant before record 1 rank that the fundametal compoment of strong first arrival and the (5) (6) step obtain, the 2 contour order harmonic component in rank, according to the sliding time gathered, from forward gun, deduct these harmonic components in relevant position, obtain the relevant front record not having harmonic pollution;
(8) the front fundametal compoment recorded and extract from forward gun force signal of being correlated with of harmonic pollution that do not have that (7) step obtains is carried out cross-correlation, just obtain the relevant rear normally single shot record (Fig. 4 is shown) after forward gun compacting harmonic wave.
The application of the present invention in experimental data, main acquisition parameter is as shown in the table:
Acquisition parameter table
Focus group number | 4 groups |
Sweep length | 12s |
Sliding time | 7s |
Record length | 7s |
Scan frequency | 10--72Hz |
Record length | 7s |
Track pitch | 25m |
Big gun distance | 25m |
Number of channels | 297 roads (281 recorded traces, 16 auxiliary channels) |
Data type | Non-related data (SEGD form) |
To the data analysis gathered, the auxiliary channel in the data that collection in worksite obtains comprises following 4 contents: reference signal, force signal, weight signal, dull and stereotyped signal; Phase selection close before nearly shot point road carry out time frequency analysis (as Fig. 1), as seen from the harmonic pollution of a big gun below.After relevant there is very strong harmonic energy (as Fig. 3) in single big gun mid-deep strata, because high-frequency harmonic decay is very fast, so harmonic pollution mainly concentrates within the scope of nearly geophone offset.
After implementing the present invention, from time frequency analysis result, the harmonic components of 1 time, the 2 times high frequencies within 5 times is all by obviously compacting (shown in Fig. 2).Single shot record (shown in Fig. 4) after compression process significantly improves with original single big gun specific mass quality, and the strong high-frequency harmonic eliminating mid-deep strata pollutes, and these harmonic pollutions all concentrate near migration range scope (as shown in Figure 5).
The degree very nearly the same without Harmonic Theory data reaching and simulate suppressed by harmonic wave of the present invention, can well main harmonic pollution constituents extraction out, and well mute.
Claims (1)
1., based on a slip scan harmonic suppression method for sparse inversion, feature is realized by following concrete steps:
1) controlled source acquiring seismic data, data and force signal before slip scan is relevant;
Before described slip scan is relevant, data and force signal first process from last big gun, short time discrete Fourier transform is utilized to be separated at time-frequency domain to the force signal of last big gun, isolate first-harmonic, first harmonic, second harmonic, triple-frequency harmonics composition, by isolated first-harmonic composition to relevant before record and makes cross-correlation, record x after obtaining being correlated with
0(t);
2) for together, after relevant, x is recorded
0t take absolute value in () maximum x
max, search from sequence | x
0(t) | > α x
max, formed with x
maxcentered by one section of region, determine the strong first arrival of controlled source;
Wherein: x
maxfor maximum absolute value, α is the constant of <1;
3) utilize isolated first-harmonic composition from force signal to carry out auto-correlation and obtain zero-phase wavelet, utilize this wavelet to record x after relevant
0carry out match tracing in (t), obtain the rarefaction reflectance factor of first arrival, obtain the fundametal compoment recording strong first arrival before being correlated with;
4) fundametal compoment deducting strong first arrival from relevant front record obtains remaining record
by 1 order harmonics s of itself and force signal
1t () carries out cross-correlation, obtain relevant after record x
1(t);
5) after relevant, x is recorded
1t wavelet w that () obtains by 1 order harmonics auto-correlation
1t (), adopts match tracing method in strong first arrival region, obtain the first arrival coefficient of 1 order harmonics, thus obtain strong first arrival record before correlation in 1 order harmonic component;
6) from step 4) the residue record that obtains
in deduct step 5) strong first arrival 1 order harmonic component that obtains, obtain another residue record, according to from step 4) to step 5) extract the method for 1 order harmonic component, extract strong first arrival record before correlation in 2 rank, the above high-order harmonic wave component in 3 rank or 3 rank;
7) adopt step 3) obtain relevant before record fundametal compoment and the step 5 of strong first arrival) and step 6) obtain 1 rank, the above high-order harmonic wave component in 2 rank or 3 rank, according to the sliding time gathered, from forward gun, deduct these harmonic components in relevant position, obtain the relevant front record not having harmonic pollution;
8) step 7) obtain do not have harmonic pollution relevant before record and the fundametal compoment that extracts from forward gun force signal carry out cross-correlation, normal single shot record after just to obtain after forward gun compacting harmonic wave relevant.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110140361.2A CN102798894B (en) | 2011-05-27 | 2011-05-27 | Sliding scanning harmonic suppression method based on sparse inversion |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110140361.2A CN102798894B (en) | 2011-05-27 | 2011-05-27 | Sliding scanning harmonic suppression method based on sparse inversion |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102798894A CN102798894A (en) | 2012-11-28 |
CN102798894B true CN102798894B (en) | 2015-04-29 |
Family
ID=47198049
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110140361.2A Active CN102798894B (en) | 2011-05-27 | 2011-05-27 | Sliding scanning harmonic suppression method based on sparse inversion |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102798894B (en) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103885085B (en) * | 2012-12-21 | 2017-02-08 | 中国石油集团东方地球物理勘探有限责任公司 | Method for suppressing controllable epicenter harmonic wave interference |
CN104422954A (en) * | 2013-08-20 | 2015-03-18 | 中国石油化工股份有限公司 | Sparse sampling method in seismic data regularization |
CN105277983B (en) * | 2014-06-26 | 2018-04-17 | 中石化石油工程地球物理有限公司胜利分公司 | The method of suppressing controllable seismic-source slip scan harmonic wave interference |
CN104111477B (en) * | 2014-07-04 | 2017-01-11 | 中国石油天然气集团公司 | Method and system for acquiring sliding time for harmonic suppression |
CN106199691B (en) * | 2015-05-26 | 2018-09-14 | 中石化石油工程地球物理有限公司胜利分公司 | Parallel the Forward Modeling based on controlled source slip scan method |
CN106199715B (en) * | 2015-05-30 | 2019-01-08 | 中国石油化工股份有限公司 | Strong reflection separation method and device based on model |
CN106094024A (en) * | 2016-06-21 | 2016-11-09 | 中国石油天然气集团公司 | A kind of controlled source slip scan data acquisition treatment method and device |
CN107255835B (en) * | 2017-06-14 | 2019-01-01 | 中国石油化工股份有限公司 | A method of for pressing machinery factor harmonic wave interference |
CN108957551B (en) * | 2018-07-03 | 2020-06-19 | 吉林大学 | Vibroseis harmonic suppression method based on reconstructed ground force signal |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2359363B (en) * | 2000-02-15 | 2002-04-03 | Geco Prakla | Processing simultaneous vibratory seismic data |
US7260021B1 (en) * | 2003-01-08 | 2007-08-21 | Westerngeco L.L.C. | Method of harmonic noise attenuation in correlated sweep data |
-
2011
- 2011-05-27 CN CN201110140361.2A patent/CN102798894B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN102798894A (en) | 2012-11-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102798894B (en) | Sliding scanning harmonic suppression method based on sparse inversion | |
CN101334483B (en) | Method for attenuating rayleigh wave scattered noise in earthquake data-handling | |
CN102176056B (en) | Seismic effective wave multi-domain energy compensating method based on time-frequency analysis | |
CN109669212B (en) | Seismic data processing method, stratum quality factor estimation method and device | |
CN110058305A (en) | A kind of DAS seismic data noise-reduction method based on convolutional neural networks | |
CN105911585B (en) | A kind of extracting method and device of earthquake record regular interference | |
CN101598809A (en) | A kind of self-adaptation is eliminated the method for linear programming noise and multiple reflection interference | |
CN109738951B (en) | Time-varying deconvolution method based on seismic event sub-spectrum | |
CN105572723B (en) | The design method of controlled source scanning signal based on autocorrelation wavelet | |
CN103675899A (en) | Method for expanding post-stack seismic data frequency band based on wavelet compression | |
CN105182417B (en) | A kind of surface wave separation method and system based on anatomic element analysis | |
CN102262243B (en) | Method for suppressing harmonic interference in seismic data of controlled source by filtering | |
CN110031899A (en) | Compressed sensing based weak signal extraction algorithm | |
CN106680874A (en) | Harmonic noise suppression method based on waveform morphology sparse modeling | |
CN105445801A (en) | Processing method for eliminating random noises of two dimensional seismic data | |
CN108845357A (en) | A method of the equivalent quality factor in stratum is estimated based on the synchronous wavelet transformation that squeezes | |
CN109946739A (en) | A kind of seismic profile Enhancement Method based on compressive sensing theory | |
CN106842324A (en) | A kind of optical cable coupled noise drawing method based on the sparse optimization of time-frequency domain | |
CN109164492B (en) | Method for extracting sound wave velocity of cased well stratum | |
CN109581481B (en) | Portable high-frequency controllable seismic source seismic signal harmonic interference elimination method | |
CN103076626A (en) | Wave field purification treatment method | |
CN105093282A (en) | Energy replacement surface wave suppressing method based on frequency constraint | |
CN110398778B (en) | Aviation electromagnetic shallow data correlation identification method based on equivalent sampling | |
CN102749646B (en) | Method for depth-frequency analysis of Rayleigh surface waves | |
CN102759749B (en) | Velocity analysis method for Rayleigh surface waves |
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 |