CN106842323A - A kind of slip scan harmonic wave disturbance suppression method based on scaling down processing - Google Patents
A kind of slip scan harmonic wave disturbance suppression method based on scaling down processing Download PDFInfo
- Publication number
- CN106842323A CN106842323A CN201510883681.5A CN201510883681A CN106842323A CN 106842323 A CN106842323 A CN 106842323A CN 201510883681 A CN201510883681 A CN 201510883681A CN 106842323 A CN106842323 A CN 106842323A
- Authority
- CN
- China
- Prior art keywords
- data
- window
- harmonic wave
- domain
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000012545 processing Methods 0.000 title claims abstract description 17
- 230000001629 suppression Effects 0.000 title claims abstract description 16
- 238000011282 treatment Methods 0.000 claims description 26
- 230000002159 abnormal effect Effects 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- 238000009826 distribution Methods 0.000 abstract description 5
- 238000013461 design Methods 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000013508 migration Methods 0.000 abstract description 2
- 230000005012 migration Effects 0.000 abstract description 2
- 238000004519 manufacturing process Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000002547 anomalous effect Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (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
It is a kind of slip scan harmonic wave disturbance suppression method based on scaling down processing the invention belongs to geophysical exploration technology.Obtain seismic data after, depending on data volume number, by whole district's packet, can one group of one group of a bunch or multiwire line, geological data is sorted into CMP domain again from big gun domain.Because the interference of slip scan harmonic wave has strong energy feature, integrated distribution is near migration range in single-shot, therefore the identification and compacting in big gun domain are more difficult, but after data are sorted into other data fields, harmonic wave interference profile relative distribution, can design corresponding incorrect energy identification and drawing method is suppressed it according to its strong energy feature.The present invention is changed by data field using the strong energy feature of slip scan harmonic wave interference, design corresponding noise identification and drawing method, and simple to operate, operation efficiency is high, production is may be directly applied to, be a kind of simple efficient slip scan harmonic wave disturbance suppression method.
Description
Technical field
The present invention relates to Processing Seismic Data field, belong to prestack denoising technology category,
Specifically a kind of slip scan harmonic wave disturbance suppression method based on scaling down processing.
Background technology
In due to controlled source slip scan acquisition mode, previous big gun focus vibrations are not over, latter big gun
Focus vibrations just have begun to, and front and rear two big guns vibrations there occurs overlap in time, therefore slip scan is gathered
There is very serious harmonic wave in mode to disturb.Harmonic wave interfering energy is strong, is mainly distributed near migration range and ground
The deep layer of shake data, and distribution concentrate very much, removal difficulty it is very big.Therefore slip scan harmonic wave interference effect
The signal to noise ratio of pre-stack data, it is necessary to suppressed it before skew.
Conventional slip scan harmonic wave disturbance suppression method, sets up harmonic prediction and calculates generally by various methods
Son, is subtracted dry to suppress harmonic wave in each order harmonicses model tried to achieve by harmonic prediction operator from initial data
Disturb.But slip scan typically uses controlled source separate unit single-shot, collection data signal to noise ratio is very low, is based on
The computational efficiency of the harmonic prediction operator of design data is low and easy distortion, so as to cause harmonic-model predicted distortion.
Estimate that true force signal can be obtained so as to obtain harmonic prediction operator based on focus backplane signal and weight signal
Certain effect, but with a distance from still having from practical application.
The content of the invention
Present invention aim at regarding to the issue above, a kind of slip scan harmonic wave based on scaling down processing is proposed
Disturbance suppression method.The method does not need the first break information of substantial amounts of artificial pickup perforation data, keeps away
The error in the artificial pick process caused because data signal to noise ratio is relatively low is exempted from, has been using three-dimensional
Geological data asks for a kind of static correcting method of ground micro-seismic static correction value.
The present invention can be realized by following technological means：
A kind of slip scan harmonic wave disturbance suppression method based on scaling down processing, including：
1) the slip scan seismic data original data packet that will be gathered；
2) data in group are sorted again；
The first-selected CMP domain of sorting, common receiver domain, common offset are alternately；
3) the data field application in sorting is based on the harmonic suppression method of split-band；I.e.：
Data field road collection first in sorting divides time-domain window, harmonic wave is carried out in window and is done
The identification and compacting disturbed, when V ≠ 0, the computing formula at window top bottom is：
(formula 1)
Wherein, X is offset distance, and T (X) is the time at offset distance X, and DT is the delay of window
Time, T0It it is time of the window at 0 offset distance, V is the velocity amplitude at window time；
Data track transformed into frequency domain in each window and split-band treatment is carried out；
The mean amplitude of tide of each frequency band Nei Gedao is calculated, computing formula is：
(formula 2)
Wherein, EftkIt is the mean amplitude of tide of frequency band Nei Gedao, nftkIt is per pass sampled point number,For
The amplitude of each sampled point；
Then application space medium filtering recognizes abnormal harmonics amplitude in the window, will identify that
Harmonic wave is disturbed into row interpolation or attenuation processing；
When all frequency band treatment are completed all frequency band Data Synthesis Whole frequency band data later in window
Time-domain is transformed into so as to complete a treatment for time-domain window；
After the completion of all window treatments, that is, complete a treatment for road collection；
4) after the treatment of whole road collections is completed, the treatment of data in group is just completed；
5) repeat step 2 to 4, until total data all treatment completions；
6) data are sorted into big gun domain again；
7) geological data is exported.
Such scheme is further included：
Divide time-domain window and be divided into hyperbolic form using offset distance and velocity information, window has
Certain length, and window has certain overlap in time orientation up and down；Data track is transformed into frequency
Rate domain and carry out split-band be according to 4Hz step-lengths divide.
Time window length selects 200-300ms, window to overlap selection up and down and overlap 50%.
To group in data sorted again when selection road collection on below 3s and 3000 meter of offset distance with
Interior data carry out subsequent treatment.
Energy and characteristic distributions that the present invention is disturbed using data field conversion and slip scan harmonic wave, using base
In the split-band harmonic wave compact technique of Fourier transformation, it is not necessary to calculate each order harmonicses predictive operator, can keep away
Exempt from harmonic prediction operator predicted distortion, it is not required that record controlled source base plate and weight signal.The method has
There are strong stability, efficiency high and calculate easy advantage, process for producing can be effectively applied to.
Brief description of the drawings
Fig. 1 is a kind of slip scan harmonic wave disturbance suppression method based on scaling down processing of the present invention
Flow chart；
Fig. 2 is the single-shot data sectional drawing comprising harmonic wave interference；
Fig. 3 is the stacked section sectional drawing comprising harmonic wave interference；
Fig. 4 is to be sorted into common midpoint gather data sectional drawing comprising design window；
Fig. 5 is frequency dividing tape handling schematic diagram；
Fig. 6 is using median filter identification harmonic wave interference schematic diagram；
Fig. 7 is rejecting harmonic wave interference way application interpolation reconstruction interference way schematic diagram；
Fig. 8 is the single-shot sectional drawing after being disturbed using this method removal harmonic wave；
Fig. 9 is the stacked section sectional drawing after being disturbed using this method removal harmonic wave.
Specific embodiment
To enable above and other objects, features and advantages of the invention to become apparent, it is cited below particularly go out compared with
Good embodiment, and coordinate institute's accompanying drawings, elaborating, specific flow is as shown in Figure 1.
A kind of slip scan harmonic wave disturbance suppression method based on scaling down processing
1) the slip scan seismic data according to collection is by original data packet；Packet Main Basiss
Data volume size and computer process ability；
2) data in group are sorted again；
The first-selected CMP domain of sorting, common receiver domain, common offset also can reach point in addition
Dissipate the effect that harmonic wave disturbs integrated distribution；
3) the data field application in sorting is based on the harmonic wave compact technique of split-band；
Concrete implementation method：Time-domain window is divided in new data field road collection first, in window
The intraoral identification and compacting for carrying out harmonic wave interference.Window is divided to be drawn using offset distance and velocity information
It is divided into hyperbolic form, it is desirable to which window has certain length, and window has one in time orientation up and down
It is fixed to overlap.When V ≠ 0, window top, the computing formula at bottom are：
(formula 1)
Wherein, X is offset distance, and T (X) is the time at offset distance X, and DT is the delay of window
Time, T0It it is time of the window at 0 offset distance, V is the velocity amplitude at window time.
Then data track transformed into frequency domain in each window and carries out split-band treatment.
Typically scaling down processing can be carried out according to 4Hz step-lengths, such as be divided into 0-4Hz, 4-8Hz, 8-12Hz,
12-16Hz……
The mean amplitude of tide of each frequency band Nei Gedao is calculated, computing formula is：
(formula 2)
Wherein, EftkIt is the mean amplitude of tide of frequency band Nei Gedao, nftkIt is per pass sampled point number,For
The amplitude of each sampled point.
Then application space medium filtering recognizes abnormal harmonics amplitude in the window, will identify that
Harmonic wave is disturbed into row interpolation or attenuation processing.When later all frequency band treatment complete in window
All frequency band Data Synthesis Whole frequency band data are transformed into time-domain so as to complete a time-domain window
Treatment.After the completion of all window treatments, that is, complete a treatment for road collection；
4) after the treatment of whole road collections is completed, the treatment of data in group is just completed；
5) repeat step 2 to 4, until total data all treatment completions；
6) data are sorted into big gun domain again；
7) geological data is conventionally exported.
Below by Application Example, the present invention will be further described again：
1) data of the XX slip scans construction of XX oil fields construction is chosen, single-shot is present in work area
More serious harmonic wave interference, harmonic wave interference is still very serious after conventional method denoising, such as schemes
More serious harmonic wave interference is there is also shown in 2, on stacked section, as shown in Figure 3；
2) 1-20 beam datas are chosen, data are sorted into CMP domain, obtain common midpoint gather,
3) time-domain window as shown in Figure 4 is designed in each road collection, window longitudinal length 300ms,
Upper and lower windows overlay 50%, divides altogether 50 windows；
4) data in each window are carried out into Fourier transformation, then divides tape handling, as shown in Figure 5；
5) ask in each frequency band per mean amplitude of tide together, anomalous amplitude recognized using median filter,
As shown in Figure 6；
6) the anomalous amplitude road that will identify that is reconstructed by the way of interpolation or decay, as shown in Figure 7.
7) repeat step 3 to 6, until all windows are all calculated and finished in road collection；
8) repeat step 3- to 7 is until all road collections are all calculated and finished.
9) common midpoint gather after harmonic wave interference will be suppressed and is sorted into big gun domain again.It is illustrated in figure 8
Except the single-shot after the interference of slip scan harmonic wave, Fig. 9 is to remove the stacked section after harmonic wave interference；
Comparison diagram 2 and Fig. 8, and Fig. 3 and Fig. 9 be can see, and the interference of slip scan harmonic wave is obtained
Compacting well, data SNR is significantly improved.
10) data output.
Claims (4)
1. a kind of slip scan harmonic wave disturbance suppression method based on scaling down processing, it is characterised in that
Including：
1) the slip scan seismic data original data packet that will be gathered；
2) data in group are sorted again；
The first-selected CMP domain of sorting, common receiver domain, common offset are alternately；
3) the data field application in sorting is based on the harmonic suppression method of split-band；I.e.：
Data field road collection first in sorting divides time-domain window, harmonic wave is carried out in window and is done
The identification and compacting disturbed, when V ≠ 0, the computing formula at window top bottom is：
Wherein, X is offset distance, and T (X) is the time at offset distance X, and DT is the delay of window
Time, T0It it is time of the window at 0 offset distance, V is the velocity amplitude at window time；
Data track transformed into frequency domain in each window and split-band treatment is carried out；
The mean amplitude of tide of each frequency band Nei Gedao is calculated, computing formula is：
Wherein, EftkIt is the mean amplitude of tide of frequency band Nei Gedao, nftkIt is per pass sampled point number,For
The amplitude of each sampled point；
Then application space medium filtering recognizes abnormal harmonics amplitude in the window, will identify that
Harmonic wave is disturbed into row interpolation or attenuation processing；
When all frequency band treatment are completed all frequency band Data Synthesis Whole frequency band data later in window
Time-domain is transformed into so as to complete a treatment for time-domain window；
After the completion of all window treatments, that is, complete a treatment for road collection；
4) after the treatment of whole road collections is completed, the treatment of data in group is just completed；
5) repeat step 2 to 4, until total data all treatment completions；
6) data are sorted into big gun domain again；
7) geological data is exported.
2. the slip scan harmonic wave disturbance suppression based on scaling down processing according to claim 1
Method, it is characterised in that divide time-domain window and divide in pairs using offset distance and velocity information
Curved state, window has certain length, and window has certain overlap in time orientation up and down；Will
Data track transforms to frequency domain and carries out split-band to be divided according to 4Hz step-lengths.
3. the slip scan harmonic wave disturbance suppression method based on scaling down processing according to claim 2, its
It is characterised by, time window length selection 200-300ms, window overlaps selection and overlaps 50% up and down.
4. the slip scan harmonic wave disturbance suppression side based on scaling down processing according to claim 1,2 or 3
Method, it is characterised in that when being sorted again to data in group below 3s and 3000 meter in selection road collection
Data carry out subsequent treatment within offset distance.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510883681.5A CN106842323A (en) | 2015-12-04 | 2015-12-04 | A kind of slip scan harmonic wave disturbance suppression method based on scaling down processing |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510883681.5A CN106842323A (en) | 2015-12-04 | 2015-12-04 | A kind of slip scan harmonic wave disturbance suppression method based on scaling down processing |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106842323A true CN106842323A (en) | 2017-06-13 |
Family
ID=59149758
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510883681.5A Pending CN106842323A (en) | 2015-12-04 | 2015-12-04 | A kind of slip scan harmonic wave disturbance suppression method based on scaling down processing |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106842323A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108957551A (en) * | 2018-07-03 | 2018-12-07 | 吉林大学 | Controlled source harmonic suppression method based on reconstruct ground force signal |
| CN109061735A (en) * | 2018-05-31 | 2018-12-21 | 中石化石油工程技术服务有限公司 | A kind of adaptive controlled source slip scan harmonic suppression method |
| CN109655915A (en) * | 2017-10-11 | 2019-04-19 | 中国石油化工股份有限公司 | A kind of frequency dividing prestack time migration method and system |
| CN111239814A (en) * | 2020-01-19 | 2020-06-05 | 中国海洋大学 | Shallow profile data mechanical interference suppression method based on same-phase axis frequency division tracking smoothing |
| CN114740530A (en) * | 2021-01-07 | 2022-07-12 | 中国石油天然气股份有限公司 | Medium-high frequency quasi-linear noise suppression method and device based on hyperbolic time window constraint |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102478671A (en) * | 2010-11-23 | 2012-05-30 | 中国石油天然气集团公司 | Method for suppressing controllable seismic-source harmonic-wave interference |
| CN102692649A (en) * | 2011-03-23 | 2012-09-26 | 中国石油天然气集团公司 | Method for eliminating harmonic interference in controllable earthquake source earthquake exploration |
| CN105093328A (en) * | 2015-08-17 | 2015-11-25 | 中国石油天然气集团公司 | Slip-sweep harmonic suppression method and device |
-
2015
- 2015-12-04 CN CN201510883681.5A patent/CN106842323A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102478671A (en) * | 2010-11-23 | 2012-05-30 | 中国石油天然气集团公司 | Method for suppressing controllable seismic-source harmonic-wave interference |
| CN102692649A (en) * | 2011-03-23 | 2012-09-26 | 中国石油天然气集团公司 | Method for eliminating harmonic interference in controllable earthquake source earthquake exploration |
| CN105093328A (en) * | 2015-08-17 | 2015-11-25 | 中国石油天然气集团公司 | Slip-sweep harmonic suppression method and device |
Non-Patent Citations (4)
| Title |
|---|
| ZHANG JUNHU ET AL.: "The S/N evaluation of high-density seismic data after bin-divisible processing", 《CPS/SEG BEIJING 2009 INTERNATIONAL GEOPHYSICAL CONFERENCE & EXPOSITION》 * |
| 刘宏杰等: "基于能量差异的滑动扫描数据谐波压制方法", 《新疆石油地质》 * |
| 林娟等: "可控震源滑动扫描谐波干扰压制方法", 《石油地球物理勘探》 * |
| 赵邦六等: "《低渗透薄储层地震勘探关键技术》", 30 June 2013 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109655915A (en) * | 2017-10-11 | 2019-04-19 | 中国石油化工股份有限公司 | A kind of frequency dividing prestack time migration method and system |
| CN109061735A (en) * | 2018-05-31 | 2018-12-21 | 中石化石油工程技术服务有限公司 | A kind of adaptive controlled source slip scan harmonic suppression method |
| CN108957551A (en) * | 2018-07-03 | 2018-12-07 | 吉林大学 | Controlled source harmonic suppression method based on reconstruct ground force signal |
| CN111239814A (en) * | 2020-01-19 | 2020-06-05 | 中国海洋大学 | Shallow profile data mechanical interference suppression method based on same-phase axis frequency division tracking smoothing |
| CN114740530A (en) * | 2021-01-07 | 2022-07-12 | 中国石油天然气股份有限公司 | Medium-high frequency quasi-linear noise suppression method and device based on hyperbolic time window constraint |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106842323A (en) | A kind of slip scan harmonic wave disturbance suppression method based on scaling down processing | |
| CN106597532A (en) | Pre-stack seismic data frequency band expanding method of combining well information and horizon information | |
| CN103376464B (en) | A kind of inversion method for stratigraphic quality factor | |
| CN102998706B (en) | Method and system for attenuating seismic data random noise | |
| CN103412329B (en) | A kind of method improving seismic data resolution | |
| CN107144879A (en) | A kind of seismic wave noise-reduction method combined based on adaptive-filtering with wavelet transformation | |
| CN102681014A (en) | Regular linear interference suppressing method based on polynomial fitting | |
| CN104849756A (en) | Method for improving resolution ratio of seismic data and enhancing energy of valid weak signals | |
| CN102288994A (en) | Method for regularizing high-dimensional seismic data under constraint of Radon spectrum | |
| CN102819043A (en) | Array signal random noise adaptive model denoising method | |
| CN105676291A (en) | Multiple wave matching attenuation method based on optimized phase axis tracking | |
| CN104330826A (en) | A method for removing various noises under the condition of complex surface | |
| CN105652322A (en) | T-f-k field polarization filtering method for multi-component seismic data | |
| CN104133248B (en) | A kind of high fidelity sound wave interference drawing method | |
| CN116520419B (en) | Hot fluid crack channel identification method | |
| CN103954998B (en) | Residual amplitude compensating method based on AVO | |
| CN106019376A (en) | Seismic wave compensation method constructed through frequency driving space-variant Q value model | |
| CN104635264B (en) | Pre-stack seismic data processing method and device | |
| CN103728661B (en) | A kind of high accuracy inverse Q filtering Processing Seismic Data | |
| CN105093282B (en) | Energy replacement surface wave pressing method based on frequency constraint | |
| CN104614767A (en) | Method for correcting seismic time-varying wavelet phase based on sectional prolongation | |
| CN106019377B (en) | A kind of two-dimensional seismic survey noise remove method based on time-space domain frequency reducing model | |
| CN105068119B (en) | The damped system and its device of face ripple in low frequency seismic data | |
| CN104502977B (en) | A kind of well control protects panel height resolution seismic data processing method | |
| CN106950597B (en) | Mixing source data separation method based on the filtering of three sides |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170613 |
|
| WD01 | Invention patent application deemed withdrawn after publication |