CN109061735A - A kind of adaptive controlled source slip scan harmonic suppression method - Google Patents
A kind of adaptive controlled source slip scan harmonic suppression method Download PDFInfo
- Publication number
- CN109061735A CN109061735A CN201810552006.8A CN201810552006A CN109061735A CN 109061735 A CN109061735 A CN 109061735A CN 201810552006 A CN201810552006 A CN 201810552006A CN 109061735 A CN109061735 A CN 109061735A
- Authority
- CN
- China
- Prior art keywords
- big gun
- harmonic wave
- wave interference
- harmonic
- 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.)
- Granted
Links
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. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
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)
- Radar Systems Or Details Thereof (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a kind of adaptive controlled source slip scan harmonic suppression methods, blow out time sequencing according to SEGD trace header record, slip scan single shot record is ranked up according to time sequencing is blown out；System loads are observed to slip scan single shot record data, obtain the shot point and geophone station coordinate information of each single shot record；According to the shot point of each single shot record and geophone station coordinate and successive time sequencing, the (n+1)th big gun is calculated in the position of the n-th big gun；The spatial dimension size of the position analysis harmonic wave interference distribution of the n-th big gun is appeared in the (n+1)th big gun, calculates the radius r of harmonic wave interference coverage；The coverage of the (n+1)th big gun in time to the n-th big gun is calculated according to the temporal regularity of distribution calculation formula of harmonic wave interference and sliding time；Using space-time variation frequency median filter method in the n-th big gun by progress harmonic wave compacting in the (n+1)th big gun harmonic wave coverage.Compared with prior art, the present invention can greatly improve the speed and effect of harmonic wave compacting.
Description
Technical field
The present invention relates to the controlled source slip scan field of seismic prospecting, especially a kind of adaptive controlled source sliding
Scan harmonic suppression method.
Background technique
With the appearance of Seismic Exploration with Vibrator, people have recognized harmonic wave very early, before occurring due to slip scan
Scanning mode in, influence very little of the harmonic wave interference to data quality can almost be ignored, so the development of harmonic wave compact technique is slow
Slowly.In order to improve the production efficiency, there is the efficient acquisition technique such as slip scan, since people further investigate acquisition speed,
Sliding time is smaller and smaller, and harmonic wave interference influences technology continuous renewal that is also just increasing, therefore suppressing harmonic wave interference, especially
It is foreign countries, proposed a series of methods of compacting harmonic wave interferences in recent years.Such as: the frequency conversion changeable phases generally used in field acquisition
Scanning technique, the technology are the polar characters based on scanning signal and harmonic wave distorted signal, are swept in collection process by changing
The initial vibration phase of signal is retouched, harmonic wave interference can effectively be suppressed by the way of being superimposed after first related later,
But this method specific location for occurring for harmonic wave selects suitable parameter to suppress harmonic wave, implements relatively difficult, matter
Amount monitoring is also difficult to carry out.Li et al. people (1994) is based on linear frequency sweep mode, proposes pure phase-shift filtering (PPSM) method
Harmonic wave interference is suppressed, this method is simple and efficient intuitive stabilization, especially suitable for VSP data, but is not suitable for slip scan note
Record.In this regard, Huang Jianping et al. (2012) makes improvements, a kind of phase shift filter that can be adapted for slip sweep is proposed
Wave method, this method can effectively suppress the interference of this big gun, but not be fine for the pressing result of adjacent big gun interference.
Therefore, harmonic wave interference is mainly interfered present in current controlled source highly effective gathering data.It is directed to harmonic wave at present
Compacting usually suppresses harmonic wave interference with force signal using female record before correlation, but is not suitable for quickly answering in production
With.It can be suppressed using median filtering according to the feature of the energy exception of harmonic wave, but median filtering cannot identify harmonic wave
The position of interference generally requires the disturbance suppression for carrying out whole big gun, reduces the fidelity of computational efficiency and denoising, and denoising
The significant wave that may be injured in the process.
Summary of the invention
The invention aims to solve the deficiencies in the prior art, a kind of adaptive controlled source sliding is provided and is swept
Harmonic suppression method is retouched, automatically the range of harmonic wave interference can be identified, improves the specific aim and effect of harmonic wave compacting.
In order to achieve the above objectives, the present invention is implemented according to following technical scheme:
A kind of adaptive controlled source slip scan harmonic suppression method, the adaptive controlled source slip scan harmonic wave compacting side
Method the following steps are included:
Step 1, according to SEGD trace header record blow out time sequencing, by slip scan single shot record according to blow out time sequencing into
Row sequence；
Step 2, system loads are observed to slip scan single shot record data, obtain shot point and the detection of each single shot record
Point coordinate information；
Step 3, according to the shot point of each single shot record and geophone station coordinate and successive time sequencing, the (n+1)th big gun is calculated in the n-th big gun
Position；
Step 4, the spatial dimension size that the position analysis harmonic wave interference distribution of the n-th big gun is appeared in the (n+1)th big gun, it is dry to calculate harmonic wave
Disturb the radius r of coverage；
Step 5, the (n+1)th big gun is calculated to n-th according to the temporal regularity of distribution calculation formula of harmonic wave interference and sliding time
The coverage of big gun in time；
Step 6, using space-time variation frequency median filter method the n-th big gun by the (n+1)th big gun harmonic wave coverage carry out harmonic wave pressure
System.
Further technical solution is, in step 2, the SPS document definition observation system provided first according to field, and so
Afterwards according to the corresponding relationship of document number and SPS observation system, shot point and the detection that load obtains each single-shot are carried out to trace header
Point coordinate information.
Further technical solution is in step 3, to calculate the (n+1)th big gun in the position of the n-th big gun, calculation formula are as follows:
Taoist monastic name where harmonic wave interference are as follows:
Harm_c=int(x_sn- x_sn_c1)/r_dis+ int(x_sn+1-x_sn)/r_dis；
Arrangement number where harmonic wave interference are as follows:
Harm_rline=int(y_sn+1-y_sn_c1)/rline_dis+first_rline_no；
X_s in formulanFor the x coordinate of the n-th big gun, x_sn_c1For first of x coordinate of the n-th big gun, r_dis is track pitch from x_sn+1It is
The x coordinate of n+1 big gun；y_sn+1For the y-coordinate of the (n+1)th big gun, y_sn_c1For first of y-coordinate of the n-th big gun, rline_dis is to receive
Linear distance, first_rline_no are the wire size of first reception line.
Further technical solution is in step 4, to do in the calculated arrangement rline_no of step 3 first each
The RMS amplitude in road counts, and calculation formula is as follows:
, wherein f is time domain function, and n is sampling point number, obtains the RMS amplitude of per pass
Value En, and seek the average root-mean-square amplitude E in the adjacent road mm,En/Em=R takes the Taoist monastic name of maximum two R values to be denoted as R1, R2,
It is the range of harmonic wave interference between R1 and R2.
Further technical solution is in steps of 5, to calculate the time effects range of harmonic wave interference according to the following formula:
,,
Wherein: T1For the initial time of harmonic wave interference, T2For the termination time of harmonic wave interference, T is sweep length, and k is harmonic wave rank
It is secondary, flTo scan initial frequency, fuFrequency is terminated for scanning, W is scan blindness.
Further technical solution is the harmonic wave interference coverage calculated in step 6 according to step 4 and step 5,
Harmonic wave is suppressed using median filter method within the scope of this, median filtering calculation formula is as follows:
G(x, y)=med { f (x-k, y-l), (k, l ∈ W) }, wherein after f (x, y), g (x, y) they are respectively original image and handle
Image.
Compared with prior art, the present invention can automatically calculate the harmonic wave that the (n+1)th big gun generates in slip scan acquisition
The coverage to the n-th big gun is interfered, and accurately used time space-variant frequency dividing median filter method is done in the coverage
Compacting is disturbed, the speed and effect of harmonic wave compacting can be greatly improved.
Detailed description of the invention
Fig. 1 is the specific implementation flow chart of adaptive controlled source slip scan harmonic suppression method of the invention；
It is distribution figure of the harmonic wave interference in single shot record that Fig. 2, which is in specific embodiments of the present invention,；
Fig. 3 is to pass through the pressed single shot record figure of adaptive harmonic wave in specific embodiments of the present invention；
Fig. 4 is that harmonic wave suppresses preceding stacked profile map in specific embodiments of the present invention；
Fig. 5 is stacked profile map after harmonic wave compacting in specific embodiments of the present invention.
Specific embodiment
The invention will be further described combined with specific embodiments below, in the illustrative examples and explanation of the invention
For explaining the present invention, but it is not as a limitation of the invention.
As shown in Figure 1, a kind of adaptive controlled source slip scan harmonic suppression method, specific flow chart are as follows:
Step 101, time sequencing is blown out according to SEGD trace header record, slip scan single shot record is arranged according to the time
Sequence；
Step 102, the SPS document definition observation system provided first according to field, then according to document number and SPS observation system
Corresponding relationship, to trace header carry out load obtain each single-shot shot point and geophone station coordinate information；
Step 103, adjacent two big gun is calculated, i.e., the (n+1)th big gun is in the position of the n-th big gun, calculation formula are as follows:
Taoist monastic name where harmonic wave interference are as follows:
Harm_c=int(x_sn-x_sn_c1)/r_dis+ int(x_sn+1-x_sn) /r_dis；
Arrangement number where harmonic wave interference are as follows:
Harm_rline=int(y_sn+1-y_sn_c1)/ rline_dis+first_rline_no；
X_s in formulanFor the x coordinate of the n-th big gun, x_sn_c1For first of x coordinate of the n-th big gun, r_dis is track pitch from x_sn+1It is
The x coordinate of n+1 big gun；y_sn+1For the y-coordinate of the (n+1)th big gun, y_sn_c1For first of y-coordinate of the n-th big gun, rline_dis is to receive
Linear distance, first_rline_no are the wire size of first reception line.
It is 526 that position where calculating harmonic wave interference in the present embodiment, which is Taoist monastic name, and arrangement number is 1327.
Step 104, it is done in the calculated arrangement rline_no of step 103 first per RMS amplitude statistics together,
Calculation formula is as follows:
, f is time domain function, and n is sampling point number, obtains the RMS amplitude value E of per passn,
And seek the average root-mean-square amplitude E in the adjacent road mm,En/Em=R, takes the Taoist monastic name of maximum two R values, R1, R2, R1 and R2 it
Between be harmonic wave interference range, carried out between R1 and R2 targeting harmonic wave compacting.
It is the amplitude for counting calculating per pass in the 1327th line in this example, then according to radium computing method, calculates
Obtained radius is 15.
Step 105, the time effects range of harmonic wave interference is calculated according to the following formula:
,
Wherein: T1For the initial time of harmonic wave interference, T2For the termination time of harmonic wave interference, T is sweep length, and k is harmonic wave rank
It is secondary, flTo scan initial frequency, fuFrequency is terminated for scanning, W is scan blindness.
The Annual distribution range for the 2 subharmonic interference being calculated in the present embodiment is 2000-18000.
As shown in Fig. 2, being distribution of the harmonic wave interference in single shot record.
Step 106, the distribution of the harmonic wave calculated according to step 104 and step 105, is filtered using intermediate value in the range
Wave method suppresses harmonic wave, and median filtering calculation formula is as follows:
G(x, y)=med { f (x-k, y-l), (k, l ∈ W) },
Wherein, f (x, y), g (x, y) are respectively image after original image and processing.
By the pressed single shot record of adaptive harmonic wave as shown in figure 3, as shown in Figure 4, Figure 5, suppressing front and back for harmonic wave
Stacked section comparison, as can be seen that adaptive harmonic suppression method through the invention, can automatically count from above-mentioned figure
The coverage of harmonic wave interference is calculated, and can carry out targeting accurate harmonic wave compacting to be interior in the coverage of harmonic wave, is improved
The fidelity of harmonic wave pressing result and denoising.
The limitation that technical solution of the present invention is not limited to the above specific embodiments, it is all to do according to the technique and scheme of the present invention
Technology deformation out, falls within the scope of protection of the present invention.
Claims (6)
1. a kind of adaptive controlled source slip scan harmonic suppression method, which is characterized in that the adaptive controlled source sliding
Scan harmonic suppression method the following steps are included:
Step 1, according to SEGD trace header record blow out time sequencing, by slip scan single shot record according to blow out time sequencing into
Row sequence；
Step 2, system loads are observed to slip scan single shot record data, obtain shot point and the detection of each single shot record
Point coordinate information；
Step 3, according to the shot point of each single shot record and geophone station coordinate and successive time sequencing, the (n+1)th big gun is calculated in the n-th big gun
Position；
Step 4, the spatial dimension size that the position analysis harmonic wave interference distribution of the n-th big gun is appeared in the (n+1)th big gun, it is dry to calculate harmonic wave
Disturb the radius r of coverage；
Step 5, the (n+1)th big gun is calculated to n-th according to the temporal regularity of distribution calculation formula of harmonic wave interference and sliding time
The coverage of big gun in time；
Step 6, using space-time variation frequency median filter method the n-th big gun by the (n+1)th big gun harmonic wave coverage carry out harmonic wave pressure
System.
2. adaptive controlled source slip scan harmonic suppression method according to claim 1, which is characterized in that in step
In 2, then the SPS document definition observation system provided first according to field is closed according to document number is corresponding with SPS observation system
System carries out the shot point and geophone station coordinate information that load obtains each single-shot to trace header.
3. adaptive controlled source slip scan harmonic suppression method according to claim 1, which is characterized in that in step
In 3, the (n+1)th big gun is calculated in the position of the n-th big gun, calculation formula are as follows:
Taoist monastic name where harmonic wave interference are as follows:
Harm_c=int(x_sn- x_sn_c1)/r_dis+ int(x_sn+1-x_sn)/r_dis；
Arrangement number where harmonic wave interference are as follows:
Harm_rline=int(y_sn+1-y_sn_c1)/rline_dis+first_rline_no；
X_s in formulanFor the x coordinate of the n-th big gun, x_sn_c1For first of x coordinate of the n-th big gun, r_dis is track pitch from x_sn+1It is
The x coordinate of n+1 big gun；y_sn+1For the y-coordinate of the (n+1)th big gun, y_sn_c1For first of y-coordinate of the n-th big gun, rline_dis is to receive
Linear distance, first_rline_no are the wire size of first reception line.
4. adaptive controlled source slip scan harmonic suppression method according to claim 3, which is characterized in that in step
In 4, done in the calculated arrangement rline_no of step 3 first per RMS amplitude statistics together, calculation formula is as follows:
, wherein f is time domain function, and n is sampling point number, obtains the RMS amplitude value of per pass
En, and seek the average root-mean-square amplitude E in the adjacent road mm,En/Em=R takes the Taoist monastic name of maximum two R values to be denoted as R1, R2, R1
Between R2 be harmonic wave interference range.
5. adaptive controlled source slip scan harmonic suppression method according to claim 4, which is characterized in that in step
In 5, the time effects range of harmonic wave interference is calculated according to the following formula:
,,
Wherein: T1For the initial time of harmonic wave interference, T2For the termination time of harmonic wave interference, T is sweep length, and k is harmonic wave rank
It is secondary, flTo scan initial frequency, fuFrequency is terminated for scanning, W is scan blindness.
6. adaptive controlled source slip scan harmonic suppression method according to claim 5, which is characterized in that in step
The harmonic wave interference coverage calculated in 6 according to step 4 and step 5, in the range using median filter method to harmonic wave into
Row compacting, median filtering calculation formula are as follows:
G(x, y)=med { f (x-k, y-l), (k, l ∈ W) }, wherein after f (x, y), g (x, y) they are respectively original image and handle
Image.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810552006.8A CN109061735B (en) | 2018-05-31 | 2018-05-31 | Self-adaptive vibroseis sliding scanning harmonic suppression method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810552006.8A CN109061735B (en) | 2018-05-31 | 2018-05-31 | Self-adaptive vibroseis sliding scanning harmonic suppression method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109061735A true CN109061735A (en) | 2018-12-21 |
| CN109061735B CN109061735B (en) | 2021-07-16 |
Family
ID=64819831
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810552006.8A Active CN109061735B (en) | 2018-05-31 | 2018-05-31 | Self-adaptive vibroseis sliding scanning harmonic suppression method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109061735B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004063771A1 (en) * | 2003-01-08 | 2004-07-29 | Westerngeco, L.L.C. | Method of harmonic noise attenuation in correlated sweep data |
| CN104111477A (en) * | 2014-07-04 | 2014-10-22 | 中国石油天然气集团公司 | Method and system for acquiring sliding time for harmonic suppression |
| CN106842323A (en) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | A kind of slip scan harmonic wave disturbance suppression method based on scaling down processing |
-
2018
- 2018-05-31 CN CN201810552006.8A patent/CN109061735B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004063771A1 (en) * | 2003-01-08 | 2004-07-29 | Westerngeco, L.L.C. | Method of harmonic noise attenuation in correlated sweep data |
| CN104111477A (en) * | 2014-07-04 | 2014-10-22 | 中国石油天然气集团公司 | Method and system for acquiring sliding time for harmonic suppression |
| CN106842323A (en) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | A kind of slip scan harmonic wave disturbance suppression method based on scaling down processing |
Non-Patent Citations (3)
| Title |
|---|
| CHARLES SICKING ET AL.: "Slip Sweep Harmonic Noise Rejection on Correlated Shot Data", 《SEG HOUSTON 2009 INTERNATIONAL EXPOSITION AND ANNUAL MEETING》 * |
| 李振春等: "可控震源两种谐波产生机理与特征研究", 《石油物探》 * |
| 林娟等: "可控震源滑动扫描谐波干扰压制方法", 《石油地球物理勘探》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109061735B (en) | 2021-07-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106443802B (en) | A kind of quantitatively characterizing method and system of the Electrical imaging reservoir fracture hole body based on morphologic filtering | |
| CN102681014B (en) | Regular linear interference suppressing method based on polynomial fitting | |
| CN102819043B (en) | Array signal random noise adaptive model denoising method | |
| CN101551465B (en) | Method for adaptively recognizing and eliminating seismic exploration single-frequency interference | |
| CN103376464A (en) | Inversion method for stratigraphic quality factor | |
| CN104007469A (en) | Weak seismic signal reconstruction method based on curvelet transform | |
| CN105510975B (en) | Improve the method and device of geological data signal to noise ratio | |
| CN108983284A (en) | A kind of domain f-p ghost reflection drawing method suitable for marine tiltedly cable data | |
| CN105242309A (en) | Earthquake data regular interference suppressing method and device | |
| CN103344989B (en) | The analytical approach of impulse noise interference in vibroseis seismologic record | |
| CN107144879A (en) | A kind of seismic wave noise-reduction method combined based on adaptive-filtering with wavelet transformation | |
| CN109633752B (en) | Offshore towing cable data self-adaptive ghost wave compression method based on three-dimensional fast Radon transformation | |
| CN105866852A (en) | Remote reference magnetotelluric impedance calculation method based on correlation detection | |
| CN109885903A (en) | A kind of ground nuclear magnetic resonance signal peaks noise remove method based on model | |
| CN105259579A (en) | A high-amplitude shielding layer rejecting method based on seismic data instantaneous attributes | |
| CN102798894A (en) | Sliding scanning harmonic suppression method based on sparse inversion | |
| CN104133248B (en) | A kind of high fidelity sound wave interference drawing method | |
| CN105911590B (en) | A kind of compression process method and device of harmonic noise | |
| CN109581481B (en) | Portable high-frequency controllable seismic source seismic signal harmonic interference elimination method | |
| CN103076626B (en) | A kind of wave field purification treatment method | |
| CN106094033A (en) | The orientation seismic beam forming method of singular value decomposition | |
| CN106842323A (en) | A kind of slip scan harmonic wave disturbance suppression method based on scaling down processing | |
| CN109061735A (en) | A kind of adaptive controlled source slip scan harmonic suppression method | |
| CN110737022B (en) | Suppression method for vibroseis to excite noise of seismic data black triangle area | |
| CN104570118B (en) | A kind of automatic identification based on dual factors and the method for removing industry disturbance |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |