CN109425897A - Eliminate the method and system of seismic data outlier interference - Google Patents
Eliminate the method and system of seismic data outlier interference Download PDFInfo
- Publication number
- CN109425897A CN109425897A CN201710766068.4A CN201710766068A CN109425897A CN 109425897 A CN109425897 A CN 109425897A CN 201710766068 A CN201710766068 A CN 201710766068A CN 109425897 A CN109425897 A CN 109425897A
- Authority
- CN
- China
- Prior art keywords
- amplitude
- window
- outlier
- road
- seismic channel
- 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 36
- 238000012937 correction Methods 0.000 claims abstract description 81
- 238000001914 filtration Methods 0.000 claims abstract description 30
- 238000006467 substitution reaction Methods 0.000 claims abstract description 12
- 238000012545 processing Methods 0.000 claims description 22
- 239000012141 concentrate Substances 0.000 claims description 5
- 230000008030 elimination Effects 0.000 abstract description 10
- 238000003379 elimination reaction Methods 0.000 abstract description 10
- 230000002159 abnormal effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000002547 anomalous effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012360 testing method 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/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. 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
-
- 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
Disclose a kind of method and system of elimination seismic data outlier interference.This method may include: to pick up outlier road, and determine the time window length of the apparent velocity at the top of Mintrop wave band, Mintrop wave band based on original shot gather data；The top time of window when calculating each seismic channel corresponding Mintrop wave band, window when determining each seismic channel corresponding Mintrop wave band, window and window when back wave level when obtaining each seismic channel corresponding Mintrop wave level；The maximum value and the mean value of the maximum value amplitude of adjacent wrap in window when this, calculated amplitude modifying factor for calculating the absolute amplitude in outlier road obtain amplitude correction values, the outlier road after obtaining correction of amplitude；Spatial median filtering is carried out, median filtering value is obtained, replaces amplitude correction values with median filtering value, substitution seismic channel is obtained, obtains the seismic channel for eliminating outlier interference.The present invention can effectively eliminate the influence of outlier interference, and normal seismic channel will not be transformed, and avoid the distortion of normal earthquake track data.
Description
Technical field
The present invention relates to Seismic Data Processing Technique fields, interfere more particularly, to a kind of elimination seismic data outlier
Method and system.
Background technique
In process of seismic data processing, the interference of outlier usually causes serious shadow to the final process quality of seismic data
It rings.In general, treatment people is the shadow for being difficult to eliminate outlier interference according to common regular process flow during processing
Loud, it is necessary to be possible to eliminate the influence of outlier interference by the special processing technique for outlier interference.People treat
The processing of outlier interference also experienced two stages: originally, when handling seismic data, outlier interference generally will be present
Seismic channel is weeded out and (the amplitude for the seismic channel for having outlier to interfere is zeroized) as abnormal trace, is found to have later
The seismic channel of outlier interference also contains useful earthquake information, considers to retain the seismic channel interfered containing outlier, then considers such as
What compacting outlier interference, so that the seismic channel for having outlier to interfere be achieved the purpose that become normal earthquake road.
Currently, applied for a patent in terms of seismic data unruly-value rejecting in spite of Liu Zhicheng et al., but generally speaking,
In seism processing field, the document in relation to this respect is seldom, and examines in communication, signal processing, automatic control and failure
The data processing fields such as disconnected are many relating to how to the technical literature of rejecting outlier in data handling.Such as document 1 (Xia Weixing,
Steepest Nonlinear Tracking Differentiator filtering characteristic research electric light and control of the Yang Xiaodong based on robust filtering theory, 2010,17 (12):
The method that outlier interference is eliminated involved in 78-81) belongs to Auto Disturbances Rejection Control Technique field；(Cao Ying improves phase difference to document 2
The Research on Outliers Eliminating Algorithm Signal and Information Processing of measurement accuracy, 2006,36 (11): 31-33) it has inquired into and has been surveyed in interferometer
The method for rejecting outlier how is selected to improve phase difference precision into system.
Cause Seismic Data Processing Technique field in relation to the less reason of the article in terms of method of abnormal value removing and correction largely
On be used in current real work the method for elimination outlier interference there are places not fully up to expectations.It is well known that mesh
The method of preceding elimination outlier interference used is also transformed normal seismic data while eliminating outlier interference,
Such as the current most commonly used Seismic data processing software system of common seismic data processing technology field --- OMEGA processing is soft
Part is eliminated outlier using anomalous amplitude sound attenuation processing module (AAA processing module) and interfered, used in the module
Method is that the earthquake data before superposition of time-domain is first transformed to frequency domain, is then disappeared using the median filtering of spatial domain to reach
Except the purpose of outlier interference, the amplitude of midband is determined by the intermediate amplitude of given threshold value, will finally be passed through
The seismic data of filtered frequency domain switches back to original time-domain.Method more than utilization eliminates outlier interference and is easy to lead
It causes people while eliminating outlier interference, influence also is produced on normal seismic data, to cause to a certain extent
The distortion of seismic data can cause adverse effect to Seismic Attribute Parameters are extracted using earthquake data before superposition later in this way.Institute
With if outlier confusion area is larger in actual seismic data, it is dry to eliminate outlier to can use the processing module of AAA etc
It disturbs, but if outlier interference is relatively fewer in actual seismic data, such processing method is with regard to improper.Therefore, it is necessary to
The method and system for developing a kind of elimination seismic data outlier interference, can only need except outlier is interfered without making to normal earthquake road
At influence.
The information for being disclosed in background of invention part is merely intended to deepen the reason to general background technique of the invention
Solution, and it is known to those skilled in the art existing to be not construed as recognizing or imply that the information is constituted in any form
Technology.
Summary of the invention
The invention proposes a kind of method and system of elimination seismic data outlier interference, pass through calculated amplitude modifying factor
Son can effectively eliminate the influence of outlier interference, so that outlier road is become normal seismic channel, and normal earthquake will not be transformed
Road avoids the distortion of normal earthquake track data.
According to an aspect of the invention, it is proposed that a kind of method for eliminating the interference of seismic data outlier.The method can be with
Include: outlier road to be picked up based on original shot gather data, and concentrate the position where Mintrop wave band according to the original big gun, determines institute
State the apparent velocity at the top of Mintrop wave band, the time window length of the Mintrop wave band；According to the apparent velocity and the original shot gather data
The offset distance of other interior each seismic channels, the top time of window when calculating each seismic channel corresponding Mintrop wave band, in conjunction with the Mintrop wave band
Time window length, window when determining each seismic channel corresponding Mintrop wave band, so when obtaining each seismic channel corresponding Mintrop wave level window with
Window when back wave level；In each seismic channel corresponding Mintrop wave level, window calculates the open country in window when back wave level respectively
It is worth the maximum value and the mean value of the maximum value amplitude of N number of adjacent wrap in window when this of the absolute amplitude in road, Jin Erji
The amplitude correc-tion factor is calculated, and by the amplitude in outlier road described in the amplitude correc-tion factor correction, obtains amplitude correction values, into
And obtain the outlier road after correction of amplitude；The window and when back wave level in window in each seismic channel corresponding Mintrop wave level respectively,
Carry out spatial median filtering, the median filtering value of outlier road window at this after obtaining the correction of amplitude, in the amplitude school
The amplitude correction values are replaced with the median filtering value in outlier road after just, obtain substitution seismic channel, and then it is wild to obtain elimination
It is worth the seismic channel of interference.
Preferably, the calculated amplitude modifying factor include: based on when this in window N number of adjacent wrap per pass absolute value
The maximum value of amplitude calculates the mean value of the maximum value amplitude of N number of adjacent wrap；With N number of adjacent wrap
Maximum value amplitude mean value divided by the outlier road absolute amplitude maximum value, obtain the amplitude correc-tion because
Son.
Preferably, the amplitude correction values are obtained multiplied by the amplitude in the outlier road with the amplitude correc-tion factor.
Preferably, window and window when stating back wave level when obtaining each seismic channel corresponding Mintrop wave level by following steps:
Amplitude when Mintrop wave band corresponding to each seismic channel more than window carries out the processing that zeroizes, and will be described along the top of the Mintrop wave band
Seismic channel is evened up, window and window when back wave level when obtaining each seismic channel corresponding Mintrop wave level.
Preferably, by the substitution earthquake channel correction pull back it is flat before state, obtain the earthquake for eliminating outlier interference
Road.
According to another aspect of the invention, it is proposed that a kind of system for eliminating the interference of seismic data outlier, may include: deposit
Reservoir is stored with computer executable instructions；Processor, the processor run the executable finger of the computer in the memory
It enables, executes following steps: being based on original shot gather data, pick up outlier road, and according to where the original big gun concentration Mintrop wave band
Position determines the apparent velocity at the top of the Mintrop wave band, the time window length of the Mintrop wave band；According to the apparent velocity and the original
The offset distance of other each seismic channels in beginning shot gather data, the top time of window when calculating each seismic channel corresponding Mintrop wave band, in conjunction with
The time window length of the Mintrop wave band, window when determining each seismic channel corresponding Mintrop wave band, and then obtain the corresponding Mintrop wave of each seismic channel
Window and window when back wave level when horizontal；The window and when back wave level in window in each seismic channel corresponding Mintrop wave level respectively,
Calculate the maximum value of the absolute amplitude in the outlier road with when this in window the maximum value amplitude of N number of adjacent wrap it is equal
Value, and then calculated amplitude modifying factor, and by the amplitude in outlier road described in the amplitude correc-tion factor correction, obtain amplitude
Corrected value, and then obtain the outlier road after correction of amplitude；Window and back wave water in each seismic channel corresponding Mintrop wave level respectively
Usually in window, spatial median filtering, the median filtering value of outlier road window at this after obtaining the correction of amplitude, in institute are carried out
The amplitude correction values are replaced with the median filtering value in outlier road after stating correction of amplitude, obtain substitution seismic channel, and then obtain
The seismic channel of outlier interference must be eliminated.
Preferably, the calculated amplitude modifying factor include: based on when this in window N number of adjacent wrap per pass absolute value
The maximum value of amplitude calculates the mean value of the maximum value amplitude of N number of adjacent wrap；With N number of adjacent wrap
Maximum value amplitude mean value divided by the outlier road absolute amplitude maximum value, obtain the amplitude correc-tion because
Son.
Preferably, the amplitude correction values are obtained multiplied by the amplitude in the outlier road with the amplitude correc-tion factor.
Preferably, window and window when back wave level when obtaining each seismic channel corresponding Mintrop wave level by following steps: right
Amplitude when the corresponding Mintrop wave band of each seismic channel more than window carries out the processing that zeroizes, and will describedly along the top of the Mintrop wave band
Shake road is evened up, window and window when back wave level when obtaining each seismic channel corresponding Mintrop wave level.
Preferably, by the substitution earthquake channel correction pull back it is flat before state, obtain the earthquake for eliminating outlier interference
Road.
Methods and apparatus of the present invention has other characteristics and advantages, these characteristics and advantages are attached from what is be incorporated herein
It will be apparent in figure and subsequent specific embodiment, or will be in the attached drawing being incorporated herein and subsequent specific reality
It applies in mode and is stated in detail, the drawings and the detailed description together serve to explain specific principles of the invention.
Detailed description of the invention
Exemplary embodiment of the present is described in more detail in conjunction with the accompanying drawings, of the invention is above-mentioned and other
Purpose, feature and advantage will be apparent, wherein in exemplary embodiments of the present invention, identical reference label is usual
Represent same parts.
Fig. 1 shows the flow chart of the step of method of elimination seismic data outlier interference according to the present invention.
Fig. 2 a and Fig. 2 b respectively illustrate the original big gun collection number according to an embodiment of the invention containing outlier interference
According to road and the schematic diagram for the seismic channel for eliminating outlier interference.
Specific embodiment
The present invention will be described in more detail below with reference to accompanying drawings.Although showing the preferred embodiment of the present invention in attached drawing,
However, it is to be appreciated that may be realized in various forms the present invention and should not be limited by the embodiments set forth herein.On the contrary, providing
These embodiments are of the invention more thorough and complete in order to make, and can will fully convey the scope of the invention to ability
The technical staff in domain.
Fig. 1 shows the flow chart of the step of method of elimination seismic data outlier interference according to the present invention.
In this embodiment, the method according to the present invention for eliminating the interference of seismic data outlier may include:
Step 101, it is based on original shot gather data, picks up outlier road, and concentrates the position where Mintrop wave band according to original big gun,
Determine the apparent velocity at the top of Mintrop wave band, the time window length of Mintrop wave band.
Step 102, according to the offset distance of other each seismic channels in apparent velocity and original shot gather data, each seismic channel pair is calculated
The top time of window when the Mintrop wave band answered, in conjunction with the time window length of Mintrop wave band, window when determining each seismic channel corresponding Mintrop wave band, into
And window and window when back wave level when obtaining each seismic channel corresponding Mintrop wave level；In one example, it is obtained by following steps
Each seismic channel corresponding Mintrop wave level when window with window when stating back wave level: when Mintrop wave band corresponding to each seismic channel more than window
Amplitude carry out the processing that zeroizes, and seismic channel is evened up along the top of Mintrop wave band, it is horizontal to obtain the corresponding Mintrop wave of each seismic channel
When window and window when back wave level.
Step 103, in each seismic channel corresponding Mintrop wave level, window calculates outlier road in window when back wave level respectively
Absolute amplitude maximum value and this when window in N number of adjacent wrap maximum value amplitude mean value, and then calculate vibration
Width modifying factor, and pass through the amplitude in amplitude correc-tion factor correction outlier road, amplitude correction values are obtained, and then obtain amplitude school
Outlier road after just；In one example, calculated amplitude modifying factor includes: based on N number of adjacent wrap per pass in window when this
Absolute amplitude maximum value, calculate the mean value of the maximum value amplitude of N number of adjacent wrap；With N number of adjacent wrap
Maximum value amplitude mean value divided by outlier road absolute amplitude maximum value, obtain the amplitude correc-tion factor.At one
In example, with the amplitude correc-tion factor multiplied by the amplitude in outlier road, amplitude correction values are obtained.
Step 104, in each seismic channel corresponding Mintrop wave level, window carries out in space in window when back wave level respectively
Value filtering, median filtering value of the outlier road at this in window after obtaining correction of amplitude, the outlier road after correction of amplitude, in
Value filtering value replaces amplitude correction values, obtains substitution seismic channel, and then obtain the seismic channel for eliminating outlier interference.In an example
In, will substitute earthquake channel correction pull back it is flat before state, obtain the seismic channel for eliminating outlier interference.
Specifically, on the shot gather data for needing to eliminate outlier interference, outlier road is picked up, and by text corresponding to outlier road
Piece number and Taoist monastic name are recorded on file.Concentrate Mintrop wave band (including the more of direct wave number, refracted wave and adjacent ground surface according to original big gun
Subwave) where position, determine the apparent velocity at the top of Mintrop wave band, the time window length of Mintrop wave band, so determine Mintrop wave band when
Window range below window is exactly reflection band when Mintrop wave band, why to be divided into Mintrop wave band and reflection band, be because of back wave
Intensity often it is weaker than the intensity of those above wave very much.
According to the offset distance of other each seismic channels in apparent velocity and original shot gather data, the corresponding Mintrop wave of each seismic channel is calculated
With when window top time, in conjunction with the time window length of Mintrop wave band, window when determining each seismic channel corresponding Mintrop wave band, to each seismic channel
Amplitude when corresponding Mintrop wave band more than window carries out the processing that zeroizes, and evens up seismic channel along the top of Mintrop wave band, obtains each
Window and window when back wave level when seismic channel corresponding Mintrop wave level.
The window and when back wave level in window in each seismic channel corresponding Mintrop wave level respectively is based on when this N number of phase in window
The maximum value of the absolute amplitude of adjacent wrap calculates the mean value of the maximum value amplitude of N number of adjacent wrap；With N number of phase
The mean value of the maximum value amplitude of adjacent wrap divided by the absolute amplitude in outlier road maximum value, obtain amplitude correc-tion because
Son obtains amplitude correction values multiplied by the amplitude in outlier road with the amplitude correc-tion factor, and then obtains the outlier after correction of amplitude
Road.In each seismic channel corresponding Mintrop wave level, window carries out spatial median filtering, is shaken in window when back wave level respectively
The median filtering value of outlier road window at this after width correction in outlier road after correction of amplitude, is replaced with median filtering value
Amplitude correction values obtain substitution seismic channel, will substitute earthquake channel correction pull back it is flat before state, obtain the ground for eliminating outlier interference
Shake road.
This method can effectively eliminate the influence of outlier interference, become outlier road just by calculated amplitude modifying factor
Normal seismic channel, and normal seismic channel will not be transformed, avoid the distortion of normal earthquake track data.
Using example
A concrete application example is given below in the scheme and its effect of the embodiment of the present invention for ease of understanding.This field
It should be understood to the one skilled in the art that the example is only for the purposes of understanding the present invention, any detail is not intended to be limited in any way
The system present invention.
Fig. 2 a and Fig. 2 b respectively illustrate the original big gun collection number according to an embodiment of the invention containing outlier interference
According to road and the schematic diagram for the seismic channel for eliminating outlier interference.
The big gun set information in somewhere is selected to carry out test process, which integrates the sampling interval as 1ms, record length 6s.
Outlier road is picked up on original shot gather data, and text text is recorded in document number corresponding to outlier road and Taoist monastic name
On part bigvalue.txt.The position where Mintrop wave band is concentrated according to original big gun, determines that the apparent velocity at the top of Mintrop wave band is
The time window length of 2000m/s, Mintrop wave band calculate shot gather data road then according to the offset distance value in per pass trace header for 600ms
The time corresponding to window top when Mintrop wave band, big gun number corresponding to each road, Taoist monastic name and top time are remembered corresponding to each road
Record is in text file TopTime.txt.
According to the offset distance of other each seismic channels in apparent velocity and original shot gather data, the corresponding Mintrop wave of each seismic channel is calculated
With when window top time, in conjunction with the time window length of Mintrop wave band, window when determining each seismic channel corresponding Mintrop wave band, to each seismic channel
Amplitude when corresponding Mintrop wave band more than window carries out the processing that zeroizes, and seismic channel is evened up 300ms along the top of Mintrop wave band
Position, according to Mintrop wave band time window length determine Mintrop wave when window bottom time be 900ms, then 900ms-6000ms be back wave water
Usually window.
The window and when back wave level in window in each seismic channel corresponding Mintrop wave level respectively is based on when this N number of phase in window
The maximum value of the absolute amplitude of adjacent wrap calculates the mean value of the maximum value amplitude of N number of adjacent wrap；With N number of phase
The mean value of the maximum value amplitude of adjacent wrap divided by the absolute amplitude in outlier road maximum value, obtain amplitude correc-tion because
Son obtains amplitude correction values multiplied by the amplitude in outlier road with the amplitude correc-tion factor, and then obtains the outlier after correction of amplitude
Road.In each seismic channel corresponding Mintrop wave level, window carries out spatial median filtering, is shaken in window when back wave level respectively
The median filtering value of outlier road window at this after width correction in outlier road after correction of amplitude, is replaced with median filtering value
Amplitude correction values obtain substitution seismic channel, will substitute earthquake channel correction pull back it is flat before state, obtain the ground for eliminating outlier interference
Road is shaken, as shown in Figure 2 b.
In conclusion the present invention can effectively eliminate the influence of outlier interference, make outlier by calculated amplitude modifying factor
Road becomes normal seismic channel, and normal seismic channel will not be transformed, and avoids the distortion of normal earthquake track data.
It will be understood by those skilled in the art that above to the purpose of the description of the embodiment of the present invention only for illustratively saying
The beneficial effect of bright the embodiment of the present invention is not intended to limit embodiments of the invention to given any example.
It may include: storage according to an embodiment of the invention, providing a kind of system of elimination seismic data outlier interference
Device is stored with computer executable instructions；Processor, the computer executable instructions in processor run memory, execute with
Lower step: being based on original shot gather data, picks up outlier road, and the position where Mintrop wave band is concentrated according to original big gun, determines Mintrop wave
The time window length of the apparent velocity at the top of band, Mintrop wave band；According in apparent velocity and original shot gather data other each seismic channels it is inclined
It moves away from the top time of window determines each seismic channel in conjunction with the time window length of Mintrop wave band when calculating each seismic channel corresponding Mintrop wave band
Window when corresponding Mintrop wave band, and then window and window when back wave level when obtaining each seismic channel corresponding Mintrop wave level；Respectively each
Window and when back wave level in window when seismic channel corresponding Mintrop wave level, calculate the maximum value of the absolute amplitude in outlier road with should
When window in N number of adjacent wrap maximum value amplitude mean value, and then calculated amplitude modifying factor, and pass through amplitude correc-tion
The amplitude in factor correction outlier road obtains amplitude correction values, and then obtains the outlier road after correction of amplitude；Respectively in each earthquake
Window and when back wave level in window when the corresponding Mintrop wave level in road, carries out spatial median filtering, the outlier after obtaining correction of amplitude
The median filtering value of road window at this, the outlier road after correction of amplitude are replaced amplitude correction values with median filtering value, are replaced
For seismic channel, and then obtain the seismic channel for eliminating outlier interference.
In one example, calculated amplitude modifying factor include: based on when this in window N number of adjacent wrap per pass it is absolute
It is worth the maximum value of amplitude, calculates the mean value of the maximum value amplitude of N number of adjacent wrap；With the maximum of N number of adjacent wrap
The mean value of absolute amplitude obtains the amplitude correc-tion factor divided by the maximum value of the absolute amplitude in outlier road.
In one example, amplitude correction values are obtained multiplied by the amplitude in outlier road with the amplitude correc-tion factor.
In one example, window and back wave level is stated when obtaining each seismic channel corresponding Mintrop wave level by following steps
When window: amplitude when Mintrop wave band corresponding to each seismic channel more than window carries out the processing that zeroizes, and along the top of Mintrop wave band by ground
Shake road is evened up, window and window when back wave level when obtaining each seismic channel corresponding Mintrop wave level.
In one example, will substitute earthquake channel correction pull back it is flat before state, obtain the seismic channel for eliminating outlier interference.
The present invention can effectively eliminate the influence of outlier interference, become outlier road just by calculated amplitude modifying factor
Normal seismic channel, and normal seismic channel will not be transformed, avoid the distortion of normal earthquake track data.
It will be understood by those skilled in the art that above to the purpose of the description of the embodiment of the present invention only for illustratively saying
The beneficial effect of bright the embodiment of the present invention is not intended to limit embodiments of the invention to given any example.
Various embodiments of the present invention are described above, above description is exemplary, and non-exclusive, and
It is not limited to disclosed each embodiment.Without departing from the scope and spirit of illustrated each embodiment, for this skill
Many modifications and changes are obvious for the those of ordinary skill in art field.
Claims (10)
1. a kind of method for eliminating the interference of seismic data outlier, comprising:
Based on original shot gather data, outlier road is picked up, and concentrates the position where Mintrop wave band according to the original big gun, described in determination
The time window length of the apparent velocity at the top of Mintrop wave band, the Mintrop wave band；
According to the offset distance of other each seismic channels in the apparent velocity and the original shot gather data, it is corresponding to calculate each seismic channel
The top time of window when Mintrop wave band, in conjunction with the time window length of the Mintrop wave band, window when determining each seismic channel corresponding Mintrop wave band, into
And window and window when back wave level when obtaining each seismic channel corresponding Mintrop wave level；
In each seismic channel corresponding Mintrop wave level, window calculates the absolute value in the outlier road in window when back wave level respectively
The mean value of the maximum value of amplitude and the maximum value amplitude of N number of adjacent wrap in window when this, and then calculated amplitude modifying factor
Son, and by the amplitude in outlier road described in the amplitude correc-tion factor correction, amplitude correction values are obtained, and then obtain amplitude school
Outlier road after just；
In each seismic channel corresponding Mintrop wave level, window carries out spatial median filtering, obtains in window when back wave level respectively
The median filtering value of outlier road window at this after the correction of amplitude, the outlier road after the correction of amplitude, in described
Value filtering value replaces the amplitude correction values, obtains substitution seismic channel, and then obtain the seismic channel for eliminating outlier interference.
2. the method according to claim 1 for eliminating the interference of seismic data outlier, wherein the calculated amplitude modifying factor
Include:
Based on when this in window the absolute amplitude of N number of adjacent wrap per pass maximum value, calculate N number of adjacent wrap
The mean value of maximum value amplitude；
With the mean value of the maximum value amplitude of N number of adjacent wrap divided by the maximum of the absolute amplitude in the outlier road
Value, obtains the amplitude correc-tion factor.
3. the method according to claim 2 for eliminating the interference of seismic data outlier, wherein multiplied with the amplitude correction factor
With the amplitude in the outlier road, the amplitude correction values are obtained.
4. the method according to claim 1 for eliminating the interference of seismic data outlier, wherein obtain various regions by following steps
Window and window when the back wave level when shaking the corresponding Mintrop wave level in road:
Amplitude when Mintrop wave band corresponding to each seismic channel more than window carries out the processing that zeroizes, and will along the top of the Mintrop wave band
The seismic channel is evened up, window and window when back wave level when obtaining each seismic channel corresponding Mintrop wave level.
5. the method according to claim 4 for eliminating the interference of seismic data outlier, wherein by the substitution earthquake channel correction
The state pulled back before putting down obtains the seismic channel for eliminating outlier interference.
6. a kind of system for eliminating the interference of seismic data outlier, which is characterized in that the system includes:
Memory is stored with computer executable instructions；
Processor, the processor run the computer executable instructions in the memory, execute following steps:
Based on original shot gather data, outlier road is picked up, and concentrates the position where Mintrop wave band according to the original big gun, described in determination
The time window length of the apparent velocity at the top of Mintrop wave band, the Mintrop wave band；
According to the offset distance of other each seismic channels in the apparent velocity and the original shot gather data, it is corresponding to calculate each seismic channel
The top time of window when Mintrop wave band, in conjunction with the time window length of the Mintrop wave band, window when determining each seismic channel corresponding Mintrop wave band, into
And window and window when back wave level when obtaining each seismic channel corresponding Mintrop wave level；
In each seismic channel corresponding Mintrop wave level, window calculates the absolute value in the outlier road in window when back wave level respectively
The mean value of the maximum value of amplitude and the maximum value amplitude of N number of adjacent wrap in window when this, and then calculated amplitude modifying factor
Son, and by the amplitude in outlier road described in the amplitude correc-tion factor correction, amplitude correction values are obtained, and then obtain amplitude school
Outlier road after just；
In each seismic channel corresponding Mintrop wave level, window carries out spatial median filtering, obtains in window when back wave level respectively
The median filtering value of correction of amplitude outlier road window at this, correction of amplitude outlier road are replaced with the median filtering value
The amplitude correction values are changed, obtain substitution seismic channel, and then obtain the seismic channel for eliminating outlier interference.
7. the system according to claim 6 for eliminating the interference of seismic data outlier, wherein the calculated amplitude modifying factor
Include:
Based on when this in window the absolute amplitude of N number of adjacent wrap per pass maximum value, calculate N number of adjacent wrap
The mean value of maximum value amplitude；
With the mean value of the maximum value amplitude of N number of adjacent wrap divided by the maximum of the absolute amplitude in the outlier road
Value, obtains the amplitude correc-tion factor.
8. the system according to claim 7 for eliminating the interference of seismic data outlier, wherein multiplied with the amplitude correc-tion factor
With the amplitude in the outlier road, the amplitude correction values are obtained.
9. the system according to claim 6 for eliminating the interference of seismic data outlier, wherein obtain various regions by following steps
Window and window when the back wave level when shaking the corresponding Mintrop wave level in road:
Amplitude when Mintrop wave band corresponding to each seismic channel more than window carries out the processing that zeroizes, and will along the top of the Mintrop wave band
The seismic channel is evened up, window and window when back wave level when obtaining each seismic channel corresponding Mintrop wave level.
10. the system according to claim 9 for eliminating the interference of seismic data outlier, wherein by the substitution seismic channel school
The state just pulled back before putting down obtains the seismic channel for eliminating outlier interference.
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 true CN109425897A (en) | 2019-03-05 |
| CN109425897B 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) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111736224A (en) * | 2020-07-14 | 2020-10-02 | 西安交通大学 | Method, storage medium and equipment for suppressing linear interference of pre-stack seismic data |
| CN112147690A (en) * | 2019-06-28 | 2020-12-29 | 中国石油化工股份有限公司 | Method and system for improving shallow weak reflection signal-to-noise ratio in first wave band |
| CN112379418A (en) * | 2020-10-16 | 2021-02-19 | 山东大学 | Method and system for calculating wave velocity of seismic direct arrival wave |
| CN112925024A (en) * | 2021-01-26 | 2021-06-08 | 中国石油化工股份有限公司 | Square wave pressing method for seismic record |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130182536A1 (en) * | 2012-01-13 | 2013-07-18 | Westerngeco L.L.C. | Wavefield deghosting |
| 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 |
-
2017
- 2017-08-30 CN CN201710766068.4A patent/CN109425897B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130182536A1 (en) * | 2012-01-13 | 2013-07-18 | Westerngeco L.L.C. | Wavefield deghosting |
| 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年学术交流会论文集》 * |
Cited By (7)
| 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 |
| CN111736224A (en) * | 2020-07-14 | 2020-10-02 | 西安交通大学 | Method, storage medium and equipment for suppressing linear interference of pre-stack seismic data |
| CN111736224B (en) * | 2020-07-14 | 2021-04-20 | 西安交通大学 | Method, storage medium and equipment for suppressing linear interference of pre-stack seismic data |
| CN112379418A (en) * | 2020-10-16 | 2021-02-19 | 山东大学 | Method and system for calculating wave velocity of seismic direct arrival wave |
| CN112379418B (en) * | 2020-10-16 | 2022-06-17 | 山东大学 | Method and system for calculating wave velocity of seismic direct arrival wave |
| CN112925024A (en) * | 2021-01-26 | 2021-06-08 | 中国石油化工股份有限公司 | Square wave pressing method for seismic record |
| CN112925024B (en) * | 2021-01-26 | 2023-10-20 | 中国石油化工股份有限公司 | Square wave pressing method for seismic recording |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109425897B (en) | 2020-04-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107272066B (en) | A kind of noisy seismic signal first-arrival traveltime pick-up method and device | |
| CN109425897A (en) | Eliminate the method and system of seismic data outlier interference | |
| CN106462957B (en) | The minimizing technology and system of fringes noise in a kind of infrared image | |
| CN109310329A (en) | System and method for reducing noise caused by the stimulation artifact in the nerve signal received as neuromodulation device | |
| CN111885275B (en) | Echo cancellation method and device for voice signal, storage medium and electronic device | |
| CN108665423B (en) | Infrared image non-uniformity correction method based on guide filtering and high-pass filtering | |
| CN109743473A (en) | Video image 3 D noise-reduction method, computer installation and computer readable storage medium | |
| CN109885903B (en) | Model-based ground nuclear magnetic resonance signal spike noise removing method | |
| CN108540754A (en) | Methods, devices and systems for more video-splicings in video monitoring | |
| CN104133247B (en) | Suppressing method and device for sleeve waves in data of vertical seismic profile | |
| CN106052852A (en) | Pulse sound signal detection method and device | |
| CN102990225A (en) | Method for detecting laser welding quality in real time | |
| CN109682326A (en) | Pot seedling upright degree detection device and detection method based on depth image | |
| WO2010142262A1 (en) | Method for determining an averaged frequency-dependent transmission function for a disturbed linear time-invariant system, evaluation device and computer program product | |
| CN105654959B (en) | Adaptive filtering coefficient updating method and device | |
| CN106714067A (en) | Automatic detection method on production line and automatic detection device thereof | |
| CN106814394B (en) | Station formula wave detector combines low frequency restoration method with conventional wave detector | |
| CN108037533B (en) | A kind of seismic prospecting noise drawing method recombinating expectation log-likelihood based on block | |
| JPS63194477A (en) | Background picture extracting method | |
| CN107918155B (en) | Inverse migration analogue data TEC time error correction method and system | |
| CN100479497C (en) | Method and device of vedio noise reduction | |
| CN111419211B (en) | Method and device for suppressing pseudo Gibbs effect in ECG signal | |
| CN106707333A (en) | Prestack de-noising method and device of digital wave detector | |
| CN106361328B (en) | A kind of EEG signals extracting method under magnetic resonance environment | |
| CN109583309A (en) | Signal de-noising method, apparatus, computer equipment and storage medium |
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 |