CN112198548A - Method for establishing two-dimensional unsteady convolution filtering model - Google Patents
Method for establishing two-dimensional unsteady convolution filtering model Download PDFInfo
- Publication number
- CN112198548A CN112198548A CN202011070698.6A CN202011070698A CN112198548A CN 112198548 A CN112198548 A CN 112198548A CN 202011070698 A CN202011070698 A CN 202011070698A CN 112198548 A CN112198548 A CN 112198548A
- Authority
- CN
- China
- Prior art keywords
- dimensional unsteady
- filter
- unsteady
- filter model
- convolution
- 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
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000001914 filtration Methods 0.000 title claims description 19
- 238000013461 design Methods 0.000 claims abstract description 4
- 230000001629 suppression Effects 0.000 abstract description 6
- 238000012545 processing Methods 0.000 description 6
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000002939 conjugate gradient method Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 235000012771 pancakes Nutrition 0.000 description 1
- 229910052704 radon Inorganic materials 0.000 description 1
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
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/282—Application of seismic models, synthetic seismograms
-
- 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 invention discloses a two-dimensional unsteady convolution filter model establishing method, which adopts an unsteady filter rule noise suppression method through the two-dimensional unsteady convolution filter model establishing method to remove rule interference from seismic data, improve the signal-to-noise ratio of seismic records, effectively recover underground reflection signals and enhance the capability of the seismic signals to reflect underground structures. The invention designs the filter by utilizing the local characteristic of the effective reflected wave instead of aiming at certain specific regular interference, so that the prior knowledge of the regular interference is not needed, and any regular interference with different apparent speeds with the effective wave can be removed in a self-adaptive manner.
Description
Technical Field
The invention relates to the field of seismic data processing in oil and gas geophysical exploration, in particular to a method for establishing a two-dimensional unsteady convolution filter model.
Background
The seismic records collected in the field contain not only effective waves capable of reflecting underground geological information, but also regular interference generated by factors such as excitation, transmission and reception of seismic waves, near-surface conditions and the like, such as linear interference waves, surface waves, multiple waves and the like. Regular noise seriously affects the quality of a single shot, adversely affects subsequent data processing interpretation, and can cause structural artifacts on the superimposed profile when severe. Therefore, in order to meet the requirements of high-resolution seismic exploration, regular noise must be removed, and the data processing quality must be improved. The regular noise suppression process is a process for separating the effective signal from the linear noise, and simultaneously, the effective signal is ensured not to be damaged.
Conventional regular noise suppression methods perform noise suppression and signal enhancement, such as cut-off, band-pass filtering, f-k filtering, Radon transform, wavelet transform, etc., depending on the difference between the signal and the noise in frequency, apparent velocity, and propagation trajectory. Implicit assumptions of such denoising methods are: (1) the noise is stable in space, and the time distance curve of the noise meets a certain propagation track, such as linearity, hyperbola, parabola and the like; (2) the waveform characteristics and frequency characteristics of the noise are stable. However, in practice, the trajectory shape and frequency characteristics of regular perturbations, whether they are raw single shot records or traces after dynamic correction, are not completely steady-state. In contrast, the geometric and dynamic features of noise are point-by-point in spatial location, and appear as non-stationary features in space. This unsteady nature of seismic signals has heretofore made no accurate signal separation method, i.e., elimination of regular interference without compromising the effective signal.
At present, two-dimensional filtering is the main method for suppressing regular noise, and performs signal separation according to the regular noise and the difference of effective signal apparent speed and frequency. However, the greatest defect of the filtering method based on the two-dimensional fourier transform is that the filter cannot change along with time and space, which results in that seismic data processing personnel cannot suppress noise according to local characteristics of seismic data. For example, the conventional f-k filtering cannot adapt to unsteady characteristics of the seismic signals in time and space when regular noise is suppressed, so that the original effective reflected waves are inevitably damaged, and the application of the conventional f-k filtering to seismic signal processing is greatly reduced.
Disclosure of Invention
In order to better suppress regular interference, unsteady characteristics of seismic signals need to be considered, and an unsteady filter is designed to separate signals from regular noise. Therefore, the invention provides a method for establishing a two-dimensional unsteady convolution filter model, which allows the filter to change point by point along with time and space without limiting the type and the change mode of the two-dimensional filter. The method avoids the problem of effective signal damage caused by the overlapping of the effective signal and the regular interference in the f-k domain, and improves the effect of noise suppression.
In order to achieve the purpose, the invention provides the following technical scheme:
the method for establishing the two-dimensional unsteady convolution filter model comprises the following steps:
s1: calculating root mean square speed information H of the underground sound waves;
s2: the apparent velocity V of any point of the same-phase axis on the common shot gather is obtained by using H obtained from S1X；
S3: according to VXDesign of unstable state sector filter hP；
S4: h according to S3PDesigning self-adaptive band-pass cake-cutting filter hX；
S5: and calculating by a two-dimensional unsteady convolution filtering formula to obtain the seismic data yX without the regular noise.
In particular, said S1 is preceded by acquiring seismic data.
Specifically, the acquired seismic data is described as si(t), i ═ 1, 2, … n, where n is the number of seismic traces.
Specifically, the root mean square velocity information H is described as v (tj), where j is 1, 2,3, …, N; wherein N is the number of seismic traces.
Specifically, the visual velocity VX is described as v (t, x),
in particular, the non-stationary sector filter hP is denoted hP (t, x),
wherein f isNrIs the cut-off frequency, vrCut off apparent velocity.
Specifically, the band pass notch filter hX is described as h (t, x),
wherein f isNl、fNhRespectively representing low and high cut-off frequencies, vl～vhRepresenting the effective signal apparent velocity passband of the filter.
Specifically, the seismic data yX is described as y (t, x),
compared with the prior art, the invention has the following beneficial effects:
the two-dimensional unsteady convolution filtering model building method adopts an unsteady filtering rule noise suppression method to remove rule interference from seismic data, improve the signal-to-noise ratio of seismic records, effectively recover underground reflection signals and enhance the capacity of the seismic signals for reflecting underground structures. The invention designs the filter by utilizing the local characteristic of the effective reflected wave instead of aiming at certain specific regular interference, so that the prior knowledge of the regular interference is not needed, and any regular interference with different apparent speeds with the effective wave can be removed in a self-adaptive manner.
Drawings
FIG. 1: process flow diagram of the invention
FIG. 2: seismic records with regular noise
FIG. 3: seismic recording after regular noise removal by conventional f-k filtering
FIG. 4: seismic record with regular noise removal using two-dimensional non-stationary filtering
FIG. 5: regular noise removal by conventional f-k filtering
FIG. 6: regular noise removal using two-dimensional non-stationary filtering
FIG. 7: seismic record of actual seismic exploration acquisition
FIG. 8: the invention removes the result after regular noise
FIG. 9: the invention removes regular noise.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
(1) Acquisition of seismic data s using seismic exploration techniquesi(t), i is 1, 2, … n, where n is the number of seismic traces, as shown in fig. 2, which contains 5 reflected effective signals and 6 regular jammers;
(2) obtaining root mean square velocity information v (t) of the groundj)，j＝1，...，N；
(3) Using formulas
Solving the apparent velocity v (t, x) of any point of the same-phase axis on the common shot gather, and quickly solving by using a conjugate gradient method;
(4) designing an unsteady fan filter according to the obtained apparent velocity v (t, x):
wherein f isNrIs the cut-off frequency, vrCut off apparent velocity.
(5) Designing a self-adaptive band-pass cake-cutting filter according to the designed unsteady fan-shaped filter
Wherein f isNl、fNhRespectively representing low and high cut-off frequencies, vl～vhRepresenting the apparent velocity pass-band of the filter.
(6) The seismic data after the regular noise is removed is calculated by a two-dimensional unsteady convolution filtering formula, as shown in fig. 4:
fig. 3 and 4 are the results of denoising the model in fig. 2 by using the conventional two-dimensional steady-state f-k filtering method and the method, respectively, and fig. 5 and 6 are the noises obtained by using the conventional two-dimensional steady-state f-k filtering method and the method. By contrast, the present invention makes full use of the local characteristics of the effective signal, and adaptively selects different filters for each point in the space, so that the effective signal is retained to the maximum extent while regular interference is removed, as shown in fig. 6. However, in the conventional method, the same filter is used for the whole section, so that the apparent velocities of different spatial positions cannot be distinguished, and since the slope of the in-phase axis of the reflected wave at the offset of the cannon is greater than the minimum value of the slope of the in-phase axis of the linear interference, the apparent velocities of the effective wave and the linear interference are partially overlapped, so that the effective signal is inevitably damaged when the interference is removed, as shown in fig. 5.
This embodiment is a method for testing the feasibility of the new invention by using pre-stack seismic data of certain eastern oilfield blocks. The original prestack section is shown in FIG. 7 with a spatial sampling interval of 50m, a minimum offset of 40m, and a temporal sampling rate of 2 ms. In the single shot record, since noise is mainly a regular interference such as a surface wave and a linear interference, a filter at each point is a pancake filter. The high and low cut-off frequencies of the filter are 150Hz and 17Hz, respectively, and the pass-and-play band of the apparent velocity is (0.8 v-1.2 v), where v denotes the apparent velocity of the filter at any point.
As shown in FIG. 8, after the processing by the method, the surface wave in the original seismic data is completely eliminated, the effective signal is highlighted, and the event of the reflected wave becomes clear. In addition, the method does not damage the effective signal while removing the regular noise, and no obvious effective signal residue exists in the removed regular interference, as shown in fig. 9.
Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art will understand that various changes, modifications and substitutions can be made without departing from the spirit and scope of the invention as defined by the appended claims. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (8)
1. The method for establishing the two-dimensional unsteady convolution filter model is characterized by comprising the following steps of: the method comprises the following steps:
s1: calculating root mean square speed information H of the underground sound waves;
s2: the apparent velocity V of any point of the same-phase axis on the common shot gather is obtained by using H obtained from S1X；
S3: according to VXDesign of unstable state sector filter hP；
S4: h according to S3PDesigning self-adaptive band-pass cake-cutting filter hX；
S5: calculating to obtain seismic data y with regular noise removed by a two-dimensional unsteady convolution filtering formulaX。
2. The method of building a two-dimensional unsteady convolution filter model according to claim 1, wherein: said S1 is preceded by acquiring seismic data.
3. The two-dimensional unsteady convolution filter model building method of claim 2, characterized in that: acquired seismic data logIs represented by si(t), i ═ 1, 2, … n, where n is the number of seismic traces.
4. A method for building a two-dimensional unsteady-state convolution filter model according to claim 3, characterized in that: the root mean square velocity information H is described as v (tj), where j is 1, 2,3, …, N; wherein N is the number of seismic traces.
5. The two-dimensional unsteady convolution filter model building method of claim 4, characterized in that: the apparent velocity VXIs described as v (t, x),
6. the two-dimensional unsteady convolution filter model building method of claim 5, characterized in that: the unstable state sector filter hPIs described as hP(t,x),
Wherein f isNrIs the cut-off frequency, vrCut off apparent velocity.
7. The method of building a two-dimensional unsteady-state convolution filter model according to claim 6, wherein: the band-pass cake-cutting type filter hXIs described as h (t, x),
wherein f isNl、fNhRespectively representing low and high cut-off frequencies, vl～vhRepresenting the effective signal apparent velocity passband of the filter.
8. The two-dimensional unsteady-state convolution filter model building method of claim 7, wherein: the seismic data yXDescribed as y (t, x);
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011070698.6A CN112198548B (en) | 2020-10-09 | 2020-10-09 | Two-dimensional unsteady convolution filter model building method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011070698.6A CN112198548B (en) | 2020-10-09 | 2020-10-09 | Two-dimensional unsteady convolution filter model building method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112198548A true CN112198548A (en) | 2021-01-08 |
| CN112198548B CN112198548B (en) | 2024-01-30 |
Family
ID=74014323
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011070698.6A Active CN112198548B (en) | 2020-10-09 | 2020-10-09 | Two-dimensional unsteady convolution filter model building method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112198548B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4254480A (en) * | 1978-09-11 | 1981-03-03 | Standard Oil Company (Indiana) | Frequency independent directionally sensitive array in seismic surveying |
| CN1306621A (en) * | 1998-05-20 | 2001-08-01 | 施鲁博格控股有限公司 | Adaptive seismic noise and interference attenuation method |
| US20100114495A1 (en) * | 2008-10-31 | 2010-05-06 | Saudi Arabian Oil Company | Seismic Image Filtering Machine To Generate A Filtered Seismic Image, Program Products, And Related Methods |
-
2020
- 2020-10-09 CN CN202011070698.6A patent/CN112198548B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4254480A (en) * | 1978-09-11 | 1981-03-03 | Standard Oil Company (Indiana) | Frequency independent directionally sensitive array in seismic surveying |
| CN1306621A (en) * | 1998-05-20 | 2001-08-01 | 施鲁博格控股有限公司 | Adaptive seismic noise and interference attenuation method |
| US20100114495A1 (en) * | 2008-10-31 | 2010-05-06 | Saudi Arabian Oil Company | Seismic Image Filtering Machine To Generate A Filtered Seismic Image, Program Products, And Related Methods |
Non-Patent Citations (2)
| Title |
|---|
| HAOQI ZHAO 等: "Two-dimensional nonstationary convolutional filtering and an adaptive linear interference suppression method", SEG INTERNATIONAL EXPOSITION AND ANNUAL MEETING, pages 4710 - 4714 * |
| 李克沛 等: "组合F-K滤波器", 石油地球物理勘探, vol. 23, no. 3, pages 359 - 364 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112198548B (en) | 2024-01-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US4992995A (en) | Methods for attenuating noise in seismic data | |
| AU2009229124B2 (en) | Method for performing constrained polarization filtering | |
| CN112835103B (en) | Adaptive ghost wave removing and broadband quasi-zero phase deconvolution combined processing method and system | |
| CN110208856B (en) | Desert complex noise suppression method based on manifold partition 2D-VMD | |
| CN110658557B (en) | Seismic data surface wave suppression method based on generation of countermeasure network | |
| NO20130915A1 (en) | Removal of noise from a seismic paint | |
| CN110261910A (en) | Seismic data surface wave minimizing technology based on adaptive sparse S-transformation | |
| US20080232193A1 (en) | Methods for Noise Removal and/or Attenuation from Seismic Data by Wavelet Selection | |
| CN111708087A (en) | Method for suppressing seismic data noise based on DnCNN neural network | |
| CN110865410A (en) | Method for suppressing random noise of seismic exploration based on NAR-TFPF | |
| CN112198548B (en) | Two-dimensional unsteady convolution filter model building method | |
| CN113221746A (en) | Microseismic signal denoising method based on improved wavelet threshold function | |
| CN111257938A (en) | Time-lapse seismic virtual source wave field reconstruction method and system based on wavelet cross-correlation | |
| CN102338884B (en) | Elliptic window direction band-pass amplitude-preserved filtering data processing method in geophysical prospecting | |
| CN102338890B (en) | Round window band-pass amplitude preservation filtering data processing method in geophysical exploration | |
| CN111239814A (en) | Shallow profile data mechanical interference suppression method based on same-phase axis frequency division tracking smoothing | |
| CN111257931A (en) | Method for removing ship-passing interference noise in marine seismic exploration | |
| CN112213785B (en) | Seismic data desert noise suppression method based on feature-enhanced denoising network | |
| CN112764108B (en) | Novel seismic data noise suppression algorithm based on improved empirical wavelet transformation | |
| CN110596756B (en) | Desert seismic exploration noise suppression method based on self-adaptive mixed complex diffusion model | |
| Zhao et al. | Two-dimensional nonstationary convolutional filtering and an adaptive linear interference suppression method | |
| CN102368094B (en) | Rectangular-window band-pass amplitude-preserved filtering data processing method used in geophysical exploration | |
| US4972382A (en) | Time variant suppression of coherent noise signals in a seismic section | |
| NL2029800B1 (en) | Method for suppressing adjacent shot interference, system and computer readable medium | |
| Huaien et al. | Attenuation of marine coherent noise |
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 |