CN108680959A - A kind of adaptive ship disturbance suppression method - Google Patents
A kind of adaptive ship disturbance suppression method Download PDFInfo
- Publication number
- CN108680959A CN108680959A CN201810469197.1A CN201810469197A CN108680959A CN 108680959 A CN108680959 A CN 108680959A CN 201810469197 A CN201810469197 A CN 201810469197A CN 108680959 A CN108680959 A CN 108680959A
- Authority
- CN
- China
- Prior art keywords
- ship
- disturbance suppression
- interference
- corrected
- correlation
- 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
- 230000001629 suppression Effects 0.000 title claims abstract description 24
- 238000000034 method Methods 0.000 title claims abstract description 23
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 12
- 230000002452 interceptive effect Effects 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims abstract description 3
- 239000013535 sea water Substances 0.000 claims description 3
- 241000208340 Araliaceae Species 0.000 claims description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 2
- 235000008434 ginseng Nutrition 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000003825 pressing Methods 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/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)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention relates to a kind of adaptive ship disturbance suppression methods, it is characterised in that includes the following steps：1) based on the seismic data interfered containing ship, the location parameter of interference ship is determined by the way of cross-correlation；2) distance for interfering ship to each geophone station is calculated, and is library track by the corresponding wave detector of minimum range, be corrected using the time delay between other roads of cross-correlation calculation and the library track, and by the time delay, each seismic channel set after being corrected；3) each seismic channel set after correction is filtered, it is then counter to be corrected into disturbance suppression of navigating, obtain the seismic channel set after final ship disturbance suppression.The present invention can effectively suppress ship interference, and enhance significant wave, thus can be widely applied in marine seismic prospectiong field.
Description
Technical field
The invention belongs to marine seismic prospectiong field, ship interference, enhancing significant wave can be suppressed more particularly to a kind of
Adaptive ship disturbance suppression method.
Background technology
The interference that external source is suffered from marine seismic prospectiong to later stage seism processing and explains the huge shadow of generation
It rings, wherein ship interference is common one of external source interference, effectively removes the influence of ship interference, can greatly improve latter earthquake money
The processing accuracy of material, and then the reliability of seismic interpretation is improved, it has a very important significance.Traditional ship disturbance suppression technology
It is to be transformed into the shot gather data interfered containing ship in the trace gather of other arrangements, in the trace gather, ship interference shows as having non-
The noise of coherence, may be used random noise compacting mode suppressed, but such method to lineups vertex near
Noise energy pressing result it is poor, and be easy to damage useful signal.Another kind of common method is by human-computer interaction or outer
Source scanning etc. modes come determine ship interference track, then suppressed using linear filter method, but such method is often
Manual intervention is needed, or is needed to scan larger spatial dimension and related prior information to determine the position of interference ship, it is logical
Normal needs take a substantial amount of time and human and material resources.
Invention content
In view of the above-mentioned problems, the object of the present invention is to provide a kind of adaptive ship disturbance suppression method, in efficiently removal ship
While interference, it is effectively protected useful signal.
To achieve the above object, the present invention takes following technical scheme：A kind of adaptive ship disturbance suppression method, feature
It is to include the following steps：1) based on the seismic data interfered containing ship, the position ginseng of interference ship is determined by the way of cross-correlation
Number；2) distance for interfering ship to each geophone station is calculated, and is library track by the corresponding wave detector of minimum range, utilize cross-correlation meter
The time delay between other roads and the library track is calculated, and the time delay is corrected, each seismic channel after being corrected
Collection；3) each seismic channel set after correction is filtered, it is then counter to be corrected into disturbance suppression of navigating, obtain final ship
Seismic channel set after disturbance suppression.
In the step 1), use the mode of cross-correlation determine interference ship location parameter method for：First from containing ship
In the seismic data of interference, three signals are randomly selected；Then it is library track to take wherein 1, is determined by way of cross-correlation
In addition the time delay between twice signal and library track signal, establishes equation group；Finally the equation group of foundation is solved,
Obtain the location parameter of interference ship.
The calculation formula of location parameter of the interference ship is：
In formula, v is the speed of seawater, and (x, y) is the location parameter for interfering ship, (x0,y0) it is library track R0Corresponding geophone station
Location parameter, (xi,yi) and (xj,yj) it is respectively other twice signal Ri,RjThe location parameter of corresponding geophone station.
In the step 2), the distance of interference ship to k-th of geophone station is：
In formula, k is geophone station serial number, and k=1,2 ..., N, N are the total number of wave detector.
The invention adopts the above technical scheme, which has the following advantages：1, the present invention is true by the way of cross-correlation
Surely the position for interfering ship, without flat cable it is assumed that and need not assume to interfere ship within the scope of cable that there is well adapting to property.
2, the present invention determines the time delay of each seismic channel set by cross-correlation in the way of, be corrected, filter, it is counter correct, be easy to manage
Solution, and quickly and effectively, significant wave can be protected while removing ship interference.Thus, the present invention can be widely applied to ocean
Ship disturbance suppression field in seismic prospecting.
Description of the drawings
Fig. 1 is adaptive ship disturbance suppression method flow diagram of the invention；
Fig. 2 is the schematic diagram of ship interference；
Fig. 3 (a)~Fig. 3 (d) is the treatment effect figure of the present invention, and wherein Fig. 3 (a) is the earthquake record interfered without ship,
Fig. 3 (b) is the earthquake record containing simulation ship interference, and Fig. 3 (c) is the ship interference muted, and Fig. 3 (d) is after ship disturbance suppression
Earthquake record.
Specific implementation mode
The present invention is described in detail below with reference to the accompanying drawings and embodiments.
As shown in Figure 1, the present invention provides a kind of adaptive ship disturbance suppression method, include the following steps：
(1) in the seismic data interfered containing ship, three signal R are randomly selected0,Ri,Rj, and it is library track to take wherein 1
R0, other twice signal R is determined by way of cross-correlationi,RjWith library track signal R0Between time delay Δ ti,Δtj,
It establishes equation group and solves, obtain the position of interference ship.
As shown in Fig. 2, containing in ship interference earthquake record, it is assumed that library track R0The location parameter of corresponding geophone station is (x0,
y0), other twice Ri,RjThe location parameter of corresponding geophone station is respectively (xi,yi) and (xj,yj), by the way of cross-correlation pair
The formula that is calculated of location parameter of interference ship is：
In formula, v is the speed of seawater, and about 1500m/s, (x, y) is the location parameter for interfering ship.
(2) distance that interference ship arrives each geophone station is calculated, and is library track by the corresponding wave detector of minimum range, using mutual
Time delay between other roads of correlation computations and the library track, and the time delay is corrected, i.e., ship is interfered to corresponding to
Curve even up.Wherein, it is conventional method to the corrected method of time delay, details are not described herein by the present invention.
Interference ship to k-th of geophone station distance be：
In formula, k is geophone station serial number, and k=1,2 ..., N, N are the total number of wave detector.Remember minimum range d*It is corresponding
Wave detector be library track R*。
(3) trace gather after correction is filtered, it is then counter to be corrected into disturbance suppression of navigating, obtain final ship
Seismic channel set after disturbance suppression.Wherein, the trace gather after correction is filtered and the anti-method corrected is conventional side
Method, details are not described herein by the present invention.
Embodiment one
As shown in Fig. 3 (a)~3 (d), the present embodiment verifies the method for the present invention by taking Analog seismic data as an example.Its
In, Fig. 3 (a) is the earthquake record interfered without ship, and Fig. 3 (b) is the earthquake record containing simulation ship interference, it can be seen that ship
Interference is presented as a string of parallel curves from shallow to deep, can be determined using seismic signal cross-correlation method each road and library track it
Between time delay, and then be corrected and ship interference compacting.Fig. 3 (c) is the ship interference muted, it can be seen that pressure
Trace without significant wave in the ship interference that system is fallen.Fig. 3 (d) is the earthquake record after ship disturbance suppression, it can be seen that ship interferes
Effective compacting is arrived, significant wave has obtained effective enhancing, can improve the precision of follow-up data imaging, be conducive to seismic data
Reasonable dismissal.The present invention can be effectively removed the influence of ship interference, prominent significant wave, and can be effectively protected effective letter
Breath, can effectively improve the precision of follow-up datum processing.
The various embodiments described above are merely to illustrate the present invention, wherein the structure of each component, connection type and manufacture craft etc. are all
It can be varied from, every equivalents carried out based on the technical solution of the present invention and improvement should not exclude
Except protection scope of the present invention.
Claims (4)
1. a kind of adaptive ship disturbance suppression method, it is characterised in that include the following steps：
1) based on the seismic data interfered containing ship, the location parameter of interference ship is determined by the way of cross-correlation；
2) distance for interfering ship to each geophone station is calculated, and is library track by the corresponding wave detector of minimum range, utilize cross-correlation
The time delay between other roads and the library track is calculated, and the time delay is corrected, each earthquake after being corrected
Trace gather；
3) each seismic channel set after correction is filtered, it is then counter to be corrected into disturbance suppression of navigating, obtain final ship
Seismic channel set after disturbance suppression.
2. a kind of adaptive ship disturbance suppression method as described in claim 1, it is characterised in that：In the step 1), use
The mode of cross-correlation determines that the method for interfering the location parameter of ship is：
First from the seismic data interfered containing ship, three signals are randomly selected；
Then it is library track to take wherein 1, is determined by way of cross-correlation between other twice signal and library track signal
Time delay establishes equation group；
Finally the equation group of foundation is solved, obtains the location parameter of interference ship.
3. a kind of adaptive ship disturbance suppression method as claimed in claim 2, it is characterised in that：The position ginseng of the interference ship
Several calculation formula are：
In formula, v is the speed of seawater, and (x, y) is the location parameter for interfering ship, (x0,y0) it is library track R0The position of corresponding geophone station
Parameter is set, (xi,yi) and (xj,yj) it is respectively other twice signal Ri,RjThe location parameter of corresponding geophone station.
4. a kind of adaptive ship disturbance suppression method as described in claim 1, it is characterised in that：In the step 2), interference
The distance of ship to k-th of geophone station is：
In formula, k is geophone station serial number, and k=1,2 ..., N, N are the total number of wave detector.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810469197.1A CN108680959A (en) | 2018-05-16 | 2018-05-16 | A kind of adaptive ship disturbance suppression method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810469197.1A CN108680959A (en) | 2018-05-16 | 2018-05-16 | A kind of adaptive ship disturbance suppression method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108680959A true CN108680959A (en) | 2018-10-19 |
Family
ID=63806678
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810469197.1A Pending CN108680959A (en) | 2018-05-16 | 2018-05-16 | A kind of adaptive ship disturbance suppression method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108680959A (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0381367A2 (en) * | 1989-01-30 | 1990-08-08 | Halliburton Geophysical Services, Inc. | Marine vibrator tuneable array |
| CN1306621A (en) * | 1998-05-20 | 2001-08-01 | 施鲁博格控股有限公司 | Adaptive seismic noise and interference attenuation method |
| CN102890291A (en) * | 2012-10-09 | 2013-01-23 | 中国海洋石油总公司 | Self-adapting suppressing method for external source interference in seismic data |
| CN103852790A (en) * | 2014-03-31 | 2014-06-11 | 中国石油大学(北京) | Morphological filtering method and system of offshore earthquake data |
| CN104777015A (en) * | 2014-01-13 | 2015-07-15 | 中国石油天然气股份有限公司 | Pressing method and pressing apparatus for side reflection waves in marine earthquake data |
-
2018
- 2018-05-16 CN CN201810469197.1A patent/CN108680959A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0381367A2 (en) * | 1989-01-30 | 1990-08-08 | Halliburton Geophysical Services, Inc. | Marine vibrator tuneable array |
| CN1306621A (en) * | 1998-05-20 | 2001-08-01 | 施鲁博格控股有限公司 | Adaptive seismic noise and interference attenuation method |
| CN102890291A (en) * | 2012-10-09 | 2013-01-23 | 中国海洋石油总公司 | Self-adapting suppressing method for external source interference in seismic data |
| CN104777015A (en) * | 2014-01-13 | 2015-07-15 | 中国石油天然气股份有限公司 | Pressing method and pressing apparatus for side reflection waves in marine earthquake data |
| CN103852790A (en) * | 2014-03-31 | 2014-06-11 | 中国石油大学(北京) | Morphological filtering method and system of offshore earthquake data |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106526677B (en) | A kind of wideband reverse-time migration imaging method of marine adaptive compacting ghost reflection | |
| CN109669212B (en) | Seismic data processing method, stratum quality factor estimation method and device | |
| CN113625337B (en) | Ultra-shallow water high-precision seismic data rapid imaging method | |
| WO2009105964A1 (en) | Method of pre-stack two-dimension-like transformation of three-dimensional seismic record | |
| CN108983284A (en) | A kind of domain f-p ghost reflection drawing method suitable for marine tiltedly cable data | |
| Yao et al. | DnResNeXt network for desert seismic data denoising | |
| CN107229071B (en) | A kind of subsurface structure inversion imaging method | |
| CN104570116A (en) | Geological marker bed-based time difference analyzing and correcting method | |
| CN110389377B (en) | Microseism offset imaging positioning method based on waveform cross-correlation coefficient multiplication | |
| CN112099082B (en) | Seismic folding wave travel time inversion method for coplanar element common azimuth gather | |
| CN109975868A (en) | A kind of wave equation ghost reflection drawing method based on Taylor expansion | |
| CN116520419B (en) | Hot fluid crack channel identification method | |
| CN110780341B (en) | Anisotropic seismic imaging method | |
| CN108303741A (en) | The double inspections of earthquake based on water inspection integral merge singing drawing method | |
| CN108919345B (en) | Submarine cable land detection noise attenuation method | |
| CN112183407A (en) | Tunnel seismic wave data denoising method and system based on time-frequency domain spectral subtraction | |
| CN108680959A (en) | A kind of adaptive ship disturbance suppression method | |
| CN109212609A (en) | Near surface Noise Elimination method based on wave equation continuation | |
| CN111538082B (en) | Automatic first arrival picking method for seismic wave time-frequency domain | |
| Liu et al. | Research on a matching detection method for magnetic anomaly of underwater targets | |
| CN112394393B (en) | CRP gather data volume reconstruction method | |
| CN112327361B (en) | Inclination interference elimination method based on linear same-phase axis iterative tracking attenuation | |
| CN116660979B (en) | OBN data full waveform inversion method based on Kaiser time window integration | |
| CN114114422B (en) | Prestack seismic data noise elimination method based on directional multi-scale decomposition | |
| Yao et al. | Weak Signals Recovery of Downhole DAS with Scale-Weighted Non-Local Selective Attention |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181019 |