CN109425897B - Method and system for eliminating seismic data outlier interference - Google Patents
Method and system for eliminating seismic data outlier interference Download PDFInfo
- Publication number
- CN109425897B CN109425897B CN201710766068.4A CN201710766068A CN109425897B CN 109425897 B CN109425897 B CN 109425897B CN 201710766068 A CN201710766068 A CN 201710766068A CN 109425897 B CN109425897 B CN 109425897B
- Authority
- CN
- China
- Prior art keywords
- time window
- seismic
- value
- channel
- outlier
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
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. analysis, for interpretation, for correction
- G01V1/30—Analysis
-
- 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
-
- 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
- G01V2210/324—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
A method and system for eliminating the wild value interference of seismic data are disclosed. The method can comprise the following steps: picking up outlier channels based on the original shot gather data, and determining the visual speed of the top of the first wave band and the time window length of the first wave band; calculating the top time of the head-wave-band time window corresponding to each seismic channel, determining the head-wave-band time window corresponding to each seismic channel, and obtaining a head-wave horizontal time window and a reflected-wave horizontal time window corresponding to each seismic channel; calculating the mean value of the maximum absolute value amplitude of the outlier channel and the maximum absolute value amplitude of the adjacent normal channel in the time window, calculating an amplitude correction factor, obtaining an amplitude correction value, and obtaining the outlier channel after amplitude correction; and carrying out spatial median filtering to obtain a median filtering value, replacing the amplitude correction value with the median filtering value to obtain a substitute seismic channel, and obtaining the seismic channel for eliminating the wild value interference. The invention can effectively eliminate the influence of outlier interference, does not transform normal seismic channels and avoids the data distortion of the normal seismic channels.
Description
Technical Field
The invention relates to the technical field of seismic data processing, in particular to a method and a system for eliminating seismic data outlier interference.
Background
In the process of processing seismic data, the interference of outliers often has a serious influence on the final processing quality of the seismic data. In general, it is difficult for a processing person to eliminate the influence of outlier interference according to a normal processing flow during a processing procedure, and the influence of outlier interference may be eliminated only by a special processing technique for outlier interference. Human handling of outlier interference also goes through two stages: firstly, when processing seismic data, generally removing seismic channels with outlier interference as abnormal channels (namely, filling amplitude values of the seismic channels with outlier interference to zero), later finding that the seismic channels with outlier interference also contain useful seismic information, considering to keep the seismic channels with outlier interference, and considering how to suppress outlier interference, thereby achieving the purpose of changing the seismic channels with outlier interference into normal seismic channels.
Although Liu Shi Cheng et al have filed a patent for removing outliers of seismic data, there are few documents on this aspect in the field of seismic data processing, and many technical documents on how to remove outliers in data processing in the field of data processing such as communication, signal processing, automatic control, and fault diagnosis. For example, a method for eliminating outlier interference related to document 1 (summer satellite, dawn east. fastest tracking differentiator filter characteristic research based on robust filter theory. electro-optic and control, 2010,17(12):78-81) belongs to the technical field of active disturbance rejection control; document 2 (Cao Yi. research on a wild value elimination algorithm for improving the phase difference measurement accuracy. signal and information processing 2006,36(11):31-33) discusses how to select a method for eliminating the wild value in an interferometer direction-finding system to improve the phase difference accuracy.
The reason why the articles related to the outlier rejection method in the technical field of seismic data processing are few is that the method for eliminating outlier interference used in the current practical work is in an unsatisfactory place. As is known, the method for eliminating outlier interference at present is used to eliminate outlier interference and to reconstruct normal seismic data, for example, the OMEGA processing software, which is the most frequently used seismic data processing software system in the conventional seismic data processing technology field, utilizes an abnormal amplitude noise attenuation processing module (AAA processing module) to eliminate outlier interference, and the method used by the module is to convert pre-stack seismic data in a time domain to a frequency domain, then utilize median filtering in a space domain to achieve the purpose of eliminating outlier interference, where an amplitude value of a frequency band is determined by a given middle amplitude value of a threshold value, and finally convert the filtered seismic data in the frequency domain back to the original time domain. The method for eliminating outlier interference easily causes people to eliminate outlier interference and influence normal seismic data, so that seismic data distortion is caused to a certain extent, and adverse influence is caused on the seismic attribute parameter extraction by using the pre-stack seismic data in the future. Therefore, if the outlier interference area in the actual seismic data is large, the outlier interference can be eliminated by using a processing module such as AAA, but if the outlier interference in the actual seismic data is relatively small, such a processing method is not suitable. Therefore, there is a need to develop a method and system for eliminating outlier interference of seismic data, which can only eliminate outlier interference without affecting normal seismic traces.
The information disclosed in this background section is only for enhancement of understanding of the general background of the invention and should not be taken as an acknowledgement or any form of suggestion that this information forms the prior art already known to a person skilled in the art.
Disclosure of Invention
The invention provides a method and a system for eliminating outlier interference of seismic data, which can effectively eliminate the influence of outlier interference by calculating an amplitude correction factor, so that outlier channels become normal seismic channels, the normal seismic channels cannot be transformed, and the data distortion of the normal seismic channels is avoided.
According to one aspect of the invention, a method for eliminating outlier interference from seismic data is provided. The method may include: picking up field value channels based on original shot gather data, and determining the visual speed of the top of a first wave band and the time window length of the first wave band according to the position of the first wave band in the original shot gather; calculating the top time of a first wave band time window corresponding to each seismic channel according to the visual velocity and the offset distance of other seismic channels in the original shot gather data, determining the first wave band time window corresponding to each seismic channel by combining the length of the first wave band time window, and further obtaining a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel; respectively calculating the maximum value of the absolute value amplitude of the outlier channel and the mean value of the maximum absolute value amplitudes of N adjacent normal channels in the time window in a head wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel, further calculating an amplitude correction factor, correcting the amplitude value of the outlier channel through the amplitude correction factor, obtaining an amplitude correction value, and further obtaining the outlier channel after amplitude correction; and respectively carrying out spatial median filtering in a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel to obtain a median filtering value of the field value channel after amplitude correction in the time window, replacing the amplitude correction value with the median filtering value in the field value channel after amplitude correction to obtain a substitute seismic channel, and further obtaining the seismic channel for eliminating field value interference.
Preferably, the calculating the amplitude correction factor includes: calculating the average value of the maximum absolute value amplitudes of the N adjacent normal tracks based on the maximum value of the absolute value amplitude of each track of the N adjacent normal tracks in the time window; and dividing the average value of the maximum absolute value amplitudes of the N adjacent normal tracks by the maximum value of the absolute value amplitude of the outlier track to obtain the amplitude correction factor.
Preferably, the amplitude correction value is obtained by multiplying the amplitude correction factor by the amplitude value of the outlier channel.
Preferably, the head wave horizontal time window and the reflected wave horizontal time window corresponding to each seismic channel are obtained through the following steps: and carrying out zero filling processing on the amplitude value above the first wave band time window corresponding to each seismic channel, leveling the seismic channel along the top of the first wave band, and obtaining a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel.
Preferably, the alternative seismic channels are corrected to be in a state before being flattened, and the seismic channels with the elimination of the outlier interference are obtained.
According to another aspect of the present invention, a system for eliminating outlier interference in seismic data is provided, which may include: a memory storing computer-executable instructions; a processor executing computer executable instructions in the memory to perform the steps of: picking up field value channels based on original shot gather data, and determining the visual speed of the top of a first wave band and the time window length of the first wave band according to the position of the first wave band in the original shot gather; calculating the top time of a first wave band time window corresponding to each seismic channel according to the visual velocity and the offset distance of other seismic channels in the original shot gather data, determining the first wave band time window corresponding to each seismic channel by combining the length of the first wave band time window, and further obtaining a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel; respectively calculating the maximum value of the absolute value amplitude of the outlier channel and the mean value of the maximum absolute value amplitudes of N adjacent normal channels in the time window in a head wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel, further calculating an amplitude correction factor, correcting the amplitude value of the outlier channel through the amplitude correction factor, obtaining an amplitude correction value, and further obtaining the outlier channel after amplitude correction; and respectively carrying out spatial median filtering in a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel to obtain a median filtering value of the field value channel after amplitude correction in the time window, replacing the amplitude correction value with the median filtering value in the field value channel after amplitude correction to obtain a substitute seismic channel, and further obtaining the seismic channel for eliminating field value interference.
Preferably, the calculating the amplitude correction factor includes: calculating the average value of the maximum absolute value amplitudes of the N adjacent normal tracks based on the maximum value of the absolute value amplitude of each track of the N adjacent normal tracks in the time window; and dividing the average value of the maximum absolute value amplitudes of the N adjacent normal tracks by the maximum value of the absolute value amplitude of the outlier track to obtain the amplitude correction factor.
Preferably, the amplitude correction value is obtained by multiplying the amplitude correction factor by the amplitude value of the outlier channel.
Preferably, the head wave horizontal time window and the reflected wave horizontal time window corresponding to each seismic channel are obtained through the following steps: and carrying out zero filling processing on the amplitude value above the first wave band time window corresponding to each seismic channel, leveling the seismic channel along the top of the first wave band, and obtaining a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel.
Preferably, the alternative seismic channels are corrected to be in a state before being flattened, and the seismic channels with the elimination of the outlier interference are obtained.
The method and apparatus of the present invention have other features and advantages which will be apparent from or are set forth in detail in the accompanying drawings and the following detailed description, which are incorporated herein, and which together serve to explain certain principles of the invention.
Drawings
The above and other objects, features and advantages of the present invention will become more apparent by describing in more detail exemplary embodiments thereof with reference to the attached drawings, in which like reference numerals generally represent like parts.
FIG. 1 shows a flow chart of the steps of a method of canceling seismic data outlier interference according to the present invention.
FIGS. 2a and 2b show schematic diagrams of raw shot data traces containing outlier interference and seismic traces with outlier interference removed, respectively, according to an embodiment of the invention.
Detailed Description
The invention will be described in more detail below with reference to the accompanying drawings. While the preferred embodiments of the present invention are shown in the drawings, it should be understood that the present invention may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
FIG. 1 shows a flow chart of the steps of a method of canceling seismic data outlier interference according to the present invention.
In this embodiment, the method for eliminating seismic data outlier interference according to the invention can include:
102, calculating the top time of a first wave band time window corresponding to each seismic channel according to the apparent velocity and the offset distance of other seismic channels in the original shot gather data, determining the first wave band time window corresponding to each seismic channel by combining the length of the first wave band time window, and further obtaining a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel; in one example, the head wave horizontal time window and the reflected wave horizontal time window corresponding to each seismic trace are obtained by the following steps: and carrying out zero filling processing on the amplitude value above the time window of the head wave band corresponding to each seismic channel, leveling the seismic channel along the top of the head wave band, and obtaining the horizontal time window of the head wave and the horizontal time window of the reflected wave corresponding to each seismic channel.
103, calculating the maximum value of the absolute value amplitude of the outlier channel and the mean value of the maximum absolute value amplitudes of N adjacent normal channels in the time window in a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel, so as to calculate an amplitude correction factor, correcting the amplitude value of the outlier channel through the amplitude correction factor, obtaining an amplitude correction value, and further obtaining the outlier channel after amplitude correction; in one example, calculating the amplitude correction factor includes: calculating the average value of the maximum absolute value amplitudes of the N adjacent normal tracks based on the maximum value of the absolute value amplitude of each track of the N adjacent normal tracks in the time window; and dividing the average value of the maximum absolute value amplitudes of the N adjacent normal tracks by the maximum value of the absolute value amplitude of the outlier track to obtain the amplitude correction factor. In one example, the amplitude correction value is obtained by multiplying the amplitude value of the outlier track by an amplitude correction factor.
And 104, respectively carrying out spatial median filtering in a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel to obtain a median filtering value of the field value channel after amplitude correction in the time window, replacing the amplitude correction value with the median filtering value in the field value channel after amplitude correction to obtain a substitute seismic channel, and further obtaining the seismic channel for eliminating field value interference. In one example, the surrogate seismic traces are corrected back to a pre-flattening state to obtain seismic traces with outlier interference removed.
Specifically, on shot gather data for which outlier interference needs to be eliminated, outlier tracks are picked up, and file numbers and track numbers corresponding to the outlier tracks are recorded on files. According to the position of the initial wave band (including direct wave number, refracted wave and multiple waves near the earth surface) in the original shot set, determining the visual velocity of the top of the initial wave band, the time window length of the initial wave band and further determining the time window range of the initial wave band, wherein the reflected wave band is below the time window of the initial wave band, and therefore the reflected wave band is divided into the initial wave band and the reflected wave band because the intensity of the reflected wave is often much weaker than that of the above waves.
Calculating the top time of a first wave band time window corresponding to each seismic channel according to the apparent velocity and the offset distance of other seismic channels in the original shot gather data, determining the first wave band time window corresponding to each seismic channel by combining the length of the first wave band time window, carrying out zero filling processing on an amplitude value above the first wave band time window corresponding to each seismic channel, leveling the seismic channels along the top of the first wave band, and obtaining a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel.
Calculating the mean value of the maximum absolute value amplitudes of N adjacent normal channels based on the maximum value of the absolute value amplitudes of the N adjacent normal channels in a head wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel; and dividing the average value of the maximum absolute value amplitudes of the N adjacent normal tracks by the maximum value of the absolute value amplitude of the field value track to obtain an amplitude correction factor, multiplying the amplitude correction factor by the amplitude value of the field value track to obtain an amplitude correction value, and further obtaining the field value track after amplitude correction. And respectively carrying out spatial median filtering in a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel to obtain a median filtering value of the field value channel after amplitude correction in the time window, replacing an amplitude correction value with the median filtering value in the field value channel after amplitude correction to obtain a substitute seismic channel, correcting the substitute seismic channel to a state before leveling, and obtaining the seismic channel for eliminating field value interference.
By calculating the amplitude correction factor, the method can effectively eliminate the influence of outlier interference, change outlier channels into normal seismic channels, and avoid the data distortion of the normal seismic channels without transforming the normal seismic channels.
Application example
To facilitate understanding of the solution of the embodiments of the present invention and the effects thereof, a specific application example is given below. It will be understood by those skilled in the art that this example is merely for the purpose of facilitating an understanding of the present invention and that any specific details thereof are not intended to limit the invention in any way.
FIGS. 2a and 2b show schematic diagrams of raw shot data traces containing outlier interference and seismic traces with outlier interference removed, respectively, according to an embodiment of the invention.
And selecting shot gather data of a certain place for test processing, wherein the shot gather sampling interval is 1ms, and the recording length is 6 s.
And picking up the outlier tracks on the original shot gather data, and recording the file numbers and track numbers corresponding to the outlier tracks to the text file bigvalue. According to the position of a first wave band in an original shot gather, determining that the apparent velocity of the top of the first wave band is 2000m/s, the length of a time window of the first wave band is 600ms, then calculating the time corresponding to the top of the time window of the first wave band corresponding to each track of a shot gather data track according to the offset value in each track head, and recording the shot number, the track number and the top time corresponding to each track in a text file TopTime.
Calculating the top time of a first wave band time window corresponding to each seismic channel according to the apparent velocity and the offset distance of other seismic channels in the original shot gather data, determining the first wave band time window corresponding to each seismic channel by combining the length of the first wave band time window, carrying out zero filling processing on an amplitude value above the first wave band time window corresponding to each seismic channel, leveling the seismic channel to a position of 300ms along the top of the first wave band, determining that the bottom time of the first wave time window is 900ms according to the length of the first wave band time window, and then taking 900ms-6000ms as a horizontal time window of reflected waves.
Calculating the mean value of the maximum absolute value amplitudes of N adjacent normal channels based on the maximum value of the absolute value amplitudes of the N adjacent normal channels in a head wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel; and dividing the average value of the maximum absolute value amplitudes of the N adjacent normal tracks by the maximum value of the absolute value amplitude of the field value track to obtain an amplitude correction factor, multiplying the amplitude correction factor by the amplitude value of the field value track to obtain an amplitude correction value, and further obtaining the field value track after amplitude correction. And respectively carrying out spatial median filtering in a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel to obtain a median filtering value of the field value channel after amplitude correction in the time window, replacing an amplitude correction value with the median filtering value in the field value channel after amplitude correction to obtain a substitute seismic channel, correcting the substitute seismic channel to a state before leveling, and obtaining the seismic channel for eliminating field value interference, wherein the state is shown in fig. 2 b.
In conclusion, the invention can effectively eliminate the influence of outlier interference by calculating the amplitude correction factor, so that outlier channels become normal seismic channels, normal seismic channels cannot be transformed, and the data distortion of the normal seismic channels is avoided.
It will be appreciated by persons skilled in the art that the above description of embodiments of the invention is intended only to illustrate the benefits of embodiments of the invention and is not intended to limit embodiments of the invention to any examples given.
According to an embodiment of the present invention, a system for eliminating seismic data outlier interference is provided, which may include: a memory storing computer-executable instructions; a processor executing computer executable instructions in the memory to perform the steps of: based on the original shot gather data, picking up a field value channel, and determining the apparent velocity of the top of a first wave band and the time window length of the first wave band according to the position of the first wave band in the original shot gather; calculating the top time of a first wave band time window corresponding to each seismic channel according to the apparent velocity and the offset distance of other seismic channels in the original shot gather data, determining the first wave band time window corresponding to each seismic channel by combining the length of the first wave band time window, and further obtaining a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel; respectively calculating the maximum value of the absolute value amplitude of the outlier channel and the mean value of the maximum absolute value amplitudes of N adjacent normal channels in the time window in a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel, further calculating an amplitude correction factor, correcting the amplitude value of the outlier channel through the amplitude correction factor, obtaining an amplitude correction value, and further obtaining the outlier channel after amplitude correction; and respectively carrying out spatial median filtering in a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel to obtain a median filtering value of the field value channel after amplitude correction in the time window, replacing the amplitude correction value with the median filtering value in the field value channel after amplitude correction to obtain a substitute seismic channel, and further obtaining the seismic channel for eliminating field value interference.
In one example, calculating the amplitude correction factor includes: calculating the average value of the maximum absolute value amplitudes of the N adjacent normal tracks based on the maximum value of the absolute value amplitude of each track of the N adjacent normal tracks in the time window; and dividing the average value of the maximum absolute value amplitudes of the N adjacent normal tracks by the maximum value of the absolute value amplitude of the outlier track to obtain the amplitude correction factor.
In one example, the amplitude correction value is obtained by multiplying the amplitude value of the outlier track by an amplitude correction factor.
In one example, the head wave horizontal time window and the reflected wave horizontal time window corresponding to each seismic trace are obtained by the following steps: and carrying out zero filling processing on the amplitude value above the time window of the head wave band corresponding to each seismic channel, leveling the seismic channel along the top of the head wave band, and obtaining the horizontal time window of the head wave and the horizontal time window of the reflected wave corresponding to each seismic channel.
In one example, the surrogate seismic traces are corrected back to a pre-flattening state to obtain seismic traces with outlier interference removed.
By calculating the amplitude correction factor, the invention can effectively eliminate the influence of outlier interference, change outlier channels into normal seismic channels, and can not transform the normal seismic channels, thereby avoiding the data distortion of the normal seismic channels.
It will be appreciated by persons skilled in the art that the above description of embodiments of the invention is intended only to illustrate the benefits of embodiments of the invention and is not intended to limit embodiments of the invention to any examples given.
Having described embodiments of the present invention, the foregoing description is intended to be exemplary, not exhaustive, and not limited to the embodiments disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments.
Claims (10)
1. A method of eliminating outlier interference from seismic data, comprising:
picking up field value channels based on original shot gather data, and determining the visual speed of the top of a first wave band and the time window length of the first wave band according to the position of the first wave band in the original shot gather;
calculating the top time of a first wave band time window corresponding to each seismic channel according to the visual velocity and the offset distance of other seismic channels in the original shot gather data, determining the first wave band time window corresponding to each seismic channel by combining the length of the first wave band time window, and further obtaining a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel;
respectively calculating the maximum value of the absolute value amplitude of the outlier channel in a head wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel, calculating the mean value of the maximum absolute value amplitudes of N adjacent normal channels in the time window, further calculating an amplitude correction factor, correcting the amplitude value of the outlier channel through the amplitude correction factor, obtaining an amplitude correction value, and further obtaining the outlier channel after amplitude correction;
and respectively carrying out spatial median filtering in a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel to obtain a median filtering value of the field value channel after amplitude correction in the time window, replacing the amplitude correction value with the median filtering value in the field value channel after amplitude correction to obtain a substitute seismic channel, and further obtaining the seismic channel for eliminating field value interference.
2. The method of canceling seismic outlier interference of claim 1, wherein said calculating an amplitude correction factor comprises:
calculating the average value of the maximum absolute value amplitudes of the N adjacent normal tracks based on the maximum value of the absolute value amplitude of each track of the N adjacent normal tracks in the time window;
and dividing the average value of the maximum absolute value amplitudes of the N adjacent normal tracks by the maximum value of the absolute value amplitude of the outlier track to obtain the amplitude correction factor.
3. The method of canceling seismic data outlier interference of claim 2, wherein said amplitude correction value is obtained by multiplying said amplitude correction factor by said outlier trace amplitude value.
4. The method of eliminating seismic data outlier interference of claim 1, wherein the head wave horizontal time window and the reflected wave horizontal time window corresponding to each seismic trace are obtained by:
and carrying out zero filling processing on the amplitude value above the first wave band time window corresponding to each seismic channel, leveling the seismic channel along the top of the first wave band, and obtaining a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel.
5. The method of eliminating seismic data outlier interference of claim 4, wherein said alternate seismic traces are corrected back to a pre-leveling state to obtain said outlier interference eliminated seismic traces.
6. A system for eliminating outlier interference from seismic data, the system comprising:
a memory storing computer-executable instructions;
a processor executing computer executable instructions in the memory to perform the steps of:
picking up field value channels based on original shot gather data, and determining the visual speed of the top of a first wave band and the time window length of the first wave band according to the position of the first wave band in the original shot gather;
calculating the top time of a first wave band time window corresponding to each seismic channel according to the visual velocity and the offset distance of other seismic channels in the original shot gather data, determining the first wave band time window corresponding to each seismic channel by combining the length of the first wave band time window, and further obtaining a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel;
respectively calculating the maximum value of the absolute value amplitude of the outlier channel in a head wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel, calculating the mean value of the maximum absolute value amplitudes of N adjacent normal channels in the time window, further calculating an amplitude correction factor, correcting the amplitude value of the outlier channel through the amplitude correction factor, obtaining an amplitude correction value, and further obtaining the outlier channel after amplitude correction;
and respectively carrying out spatial median filtering in a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel to obtain a median filtering value of the amplitude correction outlier channel in the time window, wherein the amplitude correction outlier channel replaces the amplitude correction value with the median filtering value to obtain a substitute seismic channel, and further obtain the seismic channel for eliminating outlier interference.
7. The system for canceling seismic outlier interference of claim 6, wherein the calculating the amplitude correction factor comprises:
calculating the average value of the maximum absolute value amplitudes of the N adjacent normal tracks based on the maximum value of the absolute value amplitude of each track of the N adjacent normal tracks in the time window;
and dividing the average value of the maximum absolute value amplitudes of the N adjacent normal tracks by the maximum value of the absolute value amplitude of the outlier track to obtain the amplitude correction factor.
8. The system for canceling seismic data outlier interference of claim 7, wherein said amplitude correction is obtained by multiplying said amplitude correction factor by an amplitude value of said outlier trace.
9. The system for eliminating outlier interference of seismic data of claim 6, wherein the head wave horizontal time window and the reflected wave horizontal time window corresponding to each seismic trace are obtained by:
and carrying out zero filling processing on the amplitude value above the first wave band time window corresponding to each seismic channel, leveling the seismic channel along the top of the first wave band, and obtaining a first wave horizontal time window and a reflected wave horizontal time window corresponding to each seismic channel.
10. The system for canceling seismic data outlier interference of claim 9, wherein said alternate seismic traces are corrected back to a pre-leveling state to obtain said outlier canceling seismic traces.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710766068.4A CN109425897B (en) | 2017-08-30 | 2017-08-30 | Method and system for eliminating seismic data outlier interference |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710766068.4A CN109425897B (en) | 2017-08-30 | 2017-08-30 | Method and system for eliminating seismic data outlier interference |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109425897A CN109425897A (en) | 2019-03-05 |
| CN109425897B true CN109425897B (en) | 2020-04-07 |
Family
ID=65504201
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710766068.4A Active CN109425897B (en) | 2017-08-30 | 2017-08-30 | Method and system for eliminating seismic data outlier interference |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109425897B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112147690A (en) * | 2019-06-28 | 2020-12-29 | 中国石油化工股份有限公司 | Method and system for improving shallow weak reflection signal-to-noise ratio in first wave band |
| CN111736224B (en) * | 2020-07-14 | 2021-04-20 | 西安交通大学 | Method, storage medium and equipment for suppressing linear interference of pre-stack seismic data |
| CN112379418B (en) * | 2020-10-16 | 2022-06-17 | 山东大学 | Method and system for calculating wave velocity of seismic direct arrival wave |
| CN112925024B (en) * | 2021-01-26 | 2023-10-20 | 中国石油化工股份有限公司 | Square wave pressing method for seismic recording |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103308943A (en) * | 2013-05-10 | 2013-09-18 | 中国石油天然气股份有限公司 | Method and device for attenuating interbed multiples during process of processing marine seismic data |
| CN105277974A (en) * | 2014-07-23 | 2016-01-27 | 中国石油化工股份有限公司 | Stratum data interpolation method |
| CN106772601A (en) * | 2016-12-21 | 2017-05-31 | 成都理工大学 | A kind of fracture reservoir method for quick predicting based on earthquake data before superposition |
| US10444393B2 (en) * | 2015-03-27 | 2019-10-15 | Halliburton Energy Services, Inc. | Coherent noise attenuation using statistical methods |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9274239B2 (en) * | 2012-01-13 | 2016-03-01 | Westerngeco L.L.C. | Wavefield deghosting |
-
2017
- 2017-08-30 CN CN201710766068.4A patent/CN109425897B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103308943A (en) * | 2013-05-10 | 2013-09-18 | 中国石油天然气股份有限公司 | Method and device for attenuating interbed multiples during process of processing marine seismic data |
| CN105277974A (en) * | 2014-07-23 | 2016-01-27 | 中国石油化工股份有限公司 | Stratum data interpolation method |
| US10444393B2 (en) * | 2015-03-27 | 2019-10-15 | Halliburton Energy Services, Inc. | Coherent noise attenuation using statistical methods |
| CN106772601A (en) * | 2016-12-21 | 2017-05-31 | 成都理工大学 | A kind of fracture reservoir method for quick predicting based on earthquake data before superposition |
Non-Patent Citations (2)
| Title |
|---|
| 低信噪比资料的去噪技术;徐颖等;《油气地球物理实用新技术—中国石化石油勘探开发研究院南京石油物探研究所2004年学术交流会论文集》;20050930;第271-276页 * |
| 泌阳凹陷三维地震资料叠前噪声压制方法研究;吴长江;《石油地质与工程》;20091130;第23卷(第6期);第41-44页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109425897A (en) | 2019-03-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109425897B (en) | Method and system for eliminating seismic data outlier interference | |
| CN107272066B (en) | A kind of noisy seismic signal first-arrival traveltime pick-up method and device | |
| US11215700B2 (en) | Method and system for real-time motion artifact handling and noise removal for ToF sensor images | |
| US20110010098A1 (en) | Method of pre-stack two-dimension-like transformation of three-dimensional seismic record | |
| CN104793253A (en) | Airborne electromagnetic data denoising method based on mathematical morphology | |
| CN114994754B (en) | Seismic source mechanism joint inversion method based on direct wave and depth seismic phase initial motion polarity | |
| CN107367759A (en) | A kind of geological data based on architectural feature protects side denoising method | |
| CN108957553B (en) | Method and device for dynamic correction without stretching distortion by recursion correction of dynamic correction value | |
| CN110007346A (en) | Adaptive intermediate value very noisy based on wavelet decomposition interferes seismic data noise attenuation method | |
| CN103076626A (en) | Wave field purification treatment method | |
| CN111257938A (en) | Time-lapse seismic virtual source wave field reconstruction method and system based on wavelet cross-correlation | |
| CN111781644B (en) | Linear interference attenuation method and device for shallow and medium stratum seismic data | |
| CN109212609A (en) | Near surface Noise Elimination method based on wave equation continuation | |
| CN105301655B (en) | A kind of common imaging gather line noise minimizing technology and device | |
| CN111833368B (en) | Speech restoration method based on phase consistency edge detection | |
| EP3764542A1 (en) | Signal enhancer and method for weak and blind signal denoising and extracting | |
| CN112198548B (en) | Two-dimensional unsteady convolution filter model building method | |
| CN112649873B (en) | Seismic data frequency division motion correction processing method and system | |
| CN112444853B (en) | Ground microseism noise suppression method | |
| CN110737019B (en) | Processing method and system for improving velocity spectrum precision | |
| CN111562616A (en) | Seismic data scattering noise suppression method and device | |
| CN111965704A (en) | Optical cable secondary oscillation noise suppression method and system for borehole seismic data | |
| CN109061735B (en) | Self-adaptive vibroseis sliding scanning harmonic suppression method | |
| CN109991659A (en) | A kind of post-stack seismic data high frequency sequence method for automatic tracking | |
| CN108305220B (en) | Airborne infrared degraded image correction method |
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 |