CN106353816A - Seismic acquisition footprint noise suppression method and system - Google Patents
Seismic acquisition footprint noise suppression method and system Download PDFInfo
- Publication number
- CN106353816A CN106353816A CN201610648072.6A CN201610648072A CN106353816A CN 106353816 A CN106353816 A CN 106353816A CN 201610648072 A CN201610648072 A CN 201610648072A CN 106353816 A CN106353816 A CN 106353816A
- Authority
- CN
- China
- Prior art keywords
- wave
- seismic data
- geological data
- data
- stack seismic
- 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 53
- 230000001629 suppression Effects 0.000 title abstract 3
- 238000001228 spectrum Methods 0.000 claims abstract description 92
- 230000008030 elimination Effects 0.000 claims description 63
- 238000003379 elimination reaction Methods 0.000 claims description 63
- 238000001914 filtration Methods 0.000 claims description 41
- 239000011159 matrix material Substances 0.000 claims description 13
- 230000008859 change Effects 0.000 description 11
- 230000006870 function Effects 0.000 description 7
- 230000006872 improvement Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000009466 transformation Effects 0.000 description 5
- 238000013459 approach Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000000737 periodic effect Effects 0.000 description 3
- 230000017105 transposition Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013481 data capture Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 229920001296 polysiloxane Polymers 0.000 description 1
- 239000010979 ruby Substances 0.000 description 1
- 229910001750 ruby Inorganic materials 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/20—Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
-
- 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
An embodiment of the application discloses a seismic acquisition footprint noise suppression method and system; the method comprises: acquiring original post-stack seismic data; adjusting ranking mode of the original post-stack seismic data to obtain standard original post-stack seismic data; acquiring seismic data of a first time slice in the standard original post-stack seismic data, and subjecting the seismic data of the first time slice to Laplace transform to obtain Laplace transformed seismic data; subjecting the Laplace transformed seismic data to two-dimensional Fourier transform to obtain Laplace transformed seismic data wave-number spectrum, and acquiring a filter factor of the first time slice; subjecting the seismic data of the first time slice to two-dimensional Fourier transform to obtain seismic data wave-number spectrum of the first time slice; using the filter factor of the first time slice to perform acquisition footprint noise suppression on the seismic data wave-number spectrum of the first time slice. The ability to automatically recognize and suppress the acquisition footprint noise and maintain effective geological signals can be improved.
Description
Technical field
The application is related to the seismic data processing technology field in oil seismic exploration, particularly to a kind of earthquake-capturing foot
Print Noise Elimination method and system.
Background technology
Observation system used in the gatherer process of geological data can produce interference to the geological data collecting, can
Leave footprint to be visually called this interference in the seismic data, i.e. acquisition footprint.This footprint shows as folded
On the isochronous surface of geological data, regular amplitude variations occur afterwards, be a kind of systematicness being modulated on the reflected signal of stratum
Noise, it has severely impacted the degree of accuracy of seismic properties inverting.
Generally there are three kinds of approach decay acquisition footprint noises, the first approach is to make difference by acquisition system parameter setting
The road number change of offset distance is minimum；Second approach is to make difference between road collection to be superimposed minimum by pre-stack processing；The third
Approach is to suppress acquisition footprint by post-stack data processing.Seismic data process example shows, by optimizing acquisition parameter and prestack
After treatment measures, acquisition footprint is still developed.Accordingly, it would be desirable to Noise Elimination be carried out to post-stack seismic data, to improve earthquake number
According to signal to noise ratio and prominent effective reflected signal.
It is currently based on the method mainly frequency wave number (f-kx-ky) that poststack seismic exploration data suppresses acquisition footprint noise
Filtering method, this method is to wave-number domain by two-dimensional Fourier transform by original earthquake data by space field transformation, obtains
Original earthquake data wave-number spectrum, designs filtering factor according to acquisition footprint periodic characteristics, using described filtering factor to original
Geological data wave spectrum number is acquired footprint Noise Elimination, obtains seismic target earthquakes data wave spectrum number, by two-dimentional Fourier's contravariant
Described seismic target earthquakes data wave spectrum number of changing commanders is converted into spatial domain by wave-number domain, completes acquisition footprint Noise Elimination.
Inventor finds that in prior art, at least there are the following problems:
Part acquisition footprint, for example marine acquisition footprint noise, there is the feature in time with spatial variations, entirely
There is no obvious periodic distribution feature in shake data.Because prior art medium frequency wavenumber filtering method is according to acquisition footprint
Periodic characteristics design filtering factor to complete acquisition footprint Noise Elimination, so, the method is difficult to effectively identify and suppress
The described acquisition footprint with spatial variations in time, poor to geologic signals holding capacity.
Content of the invention
The purpose of the embodiment of the present application is to provide a kind of earthquake-capturing footprint Noise Elimination method and system, to improve automatically
The ability identifying and suppressing acquisition footprint noise and keep effective geologic signals.
For solving above-mentioned technical problem, the embodiment of the present application provides a kind of earthquake-capturing footprint Noise Elimination method and system
It is achieved in that
A kind of earthquake-capturing footprint Noise Elimination method, comprising:
Obtain original post-stack seismic data；
Adjust the arrangement mode of described original post-stack seismic data, obtain the original post-stack seismic data of standard；
Obtain the geological data of very first time section in the original post-stack seismic data of described standard, to very first time section
Geological data carries out Laplace transform, obtains Laplace transform geological data；
Described Laplace transform geological data is carried out two-dimensional Fourier transform, obtains Laplace transform geological data
Wave-number spectrum, according to described Laplace transform geological data wave-number spectrum, obtains the filtering factor of very first time section；
The geological data that the very first time is cut into slices carries out two-dimensional Fourier transform, obtains the geological data of very first time section
Wave-number spectrum；
Using the very first time, the filtering factor of section is acquired footprint to the geological data wave-number spectrum that the very first time cuts into slices
Noise Elimination.
In preferred version, the geological data ripple the very first time cut into slices in the filtering factor cut into slices using the described very first time
After number spectrum is acquired footprint Noise Elimination, methods described also includes:
The geological data wave-number spectrum of the very first time section after acquisition footprint Noise Elimination is transformed into space from wave-number domain
Domain, obtains middle post-stack seismic data；
Adjust the arrangement mode of described middle post-stack seismic data, obtain target post-stack seismic data.
The arrangement mode of described middle post-stack seismic data in preferred version, is adjusted by matrix transpose.
In preferred version, described according to described Laplace transform geological data wave-number spectrum, obtain very first time section
Filtering factor, comprising:
According to described Laplace transform geological data wave-number spectrum, obtain described Laplace transform geological data wave-number spectrum
In maximum wave number value and smallest wavenumber value；
According to preset rules and described maximum wave number value and described smallest wavenumber value, adjustment described Laplace transform ground
Wave number value in shake data wave-number spectrum；
The filtering factor that Laplace transform geological data wave-number spectrum after described adjustment was cut into slices as the very first time.
In preferred version, described according to preset rules and described maximum wave number value and described smallest wavenumber value, adjust institute
State the wave number value in Laplace transform geological data wave-number spectrum, comprising:
Realized using following formula:
In formula, aiRepresent i-th wave number value in the Laplace transform geological data wave-number spectrum after described adjustment；a0max
Represent the maximum wave number value in the Laplace transform geological data wave-number spectrum before described adjustment；a0minBefore representing described adjustment
Smallest wavenumber value in Laplace transform geological data wave-number spectrum；a0iRepresent the Laplace transform earthquake number before described adjustment
According to i-th wave number value in wave-number spectrum.
In preferred version, the geological data that the described filtering factor using the section of the described very first time was cut into slices to the very first time
Wave-number spectrum is acquired footprint Noise Elimination, comprising: the earthquake that the filtering factor that the very first time is cut into slices was cut into slices with the very first time
Data wave-number spectrum is multiplied in wave-number domain.
In preferred version, described geological data includes: time, amplitude, main profile number and cross-track number.
In preferred version, the described arrangement mode adjusting original post-stack seismic data, comprising: the described original poststack ground of adjustment
The arrangement mode of at least two dimensions of shake data；
Specifically, the time dimension at least two dimensions of described original post-stack seismic data is set to tie up the most slowly
Degree；At least two dimensions of described original post-stack seismic data include: described time corresponding time dimension, described main profile number
Corresponding inline dimension and described cross-track number corresponding cross-track dimension.
The arrangement mode of described original post-stack seismic data in preferred version, is adjusted by matrix transpose.
A kind of earthquake-capturing footprint Noise Elimination system, comprising: original post-stack seismic data acquiring unit, arrangement mode are adjusted
Whole unit, Laplace transform unit, filtering factor acquiring unit, Fourier transform unit and Noise Elimination unit；Wherein,
Described original post-stack seismic data acquiring unit, for obtaining original post-stack seismic data；
Described arrangement mode adjustment unit, for adjusting the arrangement mode of described original post-stack seismic data, obtains standard
Original post-stack seismic data；
Described Laplace transform unit, for obtaining very first time section in the original post-stack seismic data of described standard
Geological data, to the very first time, the geological data of section carries out Laplace transform, obtains Laplace transform geological data；
Described filtering factor acquiring unit, for carrying out two-dimentional Fourier's change by described Laplace transform geological data
Change, obtain Laplace transform geological data wave-number spectrum, according to described Laplace transform geological data wave-number spectrum, obtain first
The filtering factor of isochronous surface；
Described Fourier transform unit, carries out two-dimensional Fourier transform for the geological data that the very first time is cut into slices, obtains
Geological data wave-number spectrum to very first time section；
Described Noise Elimination unit, the earthquake number that the filtering factor for being cut into slices using the very first time was cut into slices to the very first time
It is acquired footprint Noise Elimination according to wave-number spectrum.
In preferred version, described system also includes: middle post-stack seismic data acquiring unit and target post-stack seismic data
Acquiring unit；Wherein,
Described middle post-stack seismic data acquiring unit, for by after acquisition footprint Noise Elimination the very first time section
Geological data wave-number spectrum is transformed into spatial domain from wave-number domain, obtains middle post-stack seismic data；
Seismic data acquisition cell after described object deck, for adjusting the arrangement mode of described middle post-stack seismic data,
Obtain target post-stack seismic data.
This application provides a kind of earthquake-capturing footprint Noise Elimination method and system, based on isochronous surface geological data
Laplace transform and the mode of wavenumber filtering combination, methods described is acquired footprint compacting to post-stack seismic data, can
The feature being changed with room and time according to marine streamer seismic data acquisition footprint, automatic identification is simultaneously suppressed acquisition footprint and is made an uproar
Sound, and can effectively keep geologic signals.
Brief description
In order to be illustrated more clearly that the embodiment of the present application or technical scheme of the prior art, below will be to embodiment or existing
Have technology description in required use accompanying drawing be briefly described it should be apparent that, drawings in the following description be only this
Some embodiments described in application, for those of ordinary skill in the art, in the premise not paying creative labor
Under, other accompanying drawings can also be obtained according to these accompanying drawings.
Fig. 1 is a kind of flow chart of the application earthquake-capturing footprint Noise Elimination embodiment of the method；
Fig. 2 is the schematic diagram of the geological data of isochronous surface at the 1000th millisecond of the application and its wave-number spectrum；
Fig. 3 is that the Laplce after the geological data Laplace transform of isochronous surface at the 1000th millisecond of the application becomes
Change geological data and the schematic diagram of Laplace transform geological data wave-number spectrum；
Fig. 4 be the application acquisition footprint Noise Elimination after the 1000th millisecond at the geological data of isochronous surface and its ripple
The schematic diagram of number spectrum；
Fig. 5 be the application acquisition footprint Noise Elimination before and after the 1000th millisecond at isochronous surface geological data difference and
The schematic diagram of its wave-number spectrum；
Fig. 6 is a kind of composition structure chart of the application earthquake-capturing footprint Noise Elimination system embodiment.
Specific embodiment
The embodiment of the present application provides a kind of earthquake-capturing footprint Noise Elimination method and system.
In order that those skilled in the art more fully understand the technical scheme in the application, real below in conjunction with the application
Apply the accompanying drawing in example, the enforcement it is clear that described is clearly and completely described to the technical scheme in the embodiment of the present application
Example is only some embodiments of the present application, rather than whole embodiments.Based on the embodiment in the application, this area is common
The every other embodiment that technical staff is obtained under the premise of not making creative work, all should belong to the application protection
Scope.
Fig. 1 is a kind of flow chart of the application earthquake-capturing footprint Noise Elimination embodiment of the method.As shown in figure 1, it is described
Earthquake-capturing footprint Noise Elimination method, comprises the following steps.
Step s101: obtain original post-stack seismic data.
Specifically, original post-stack seismic data can be read in calculator memory.Wherein, by earth's surface receive from ground
Under the geological data of same pip be overlapped after dynamic(al) correction, the signal to noise ratio of geological data can be effectively improved, superposition
Geological data afterwards, referred to as post-stack seismic data.Described original post-stack seismic data may include that time, amplitude, main profile
Number and cross-track (xline) number (inline).Described original post-stack seismic data is 3D seismic data, when may include that
Between direction, inline direction and xline direction.Wherein, inline direction and xline direction in described original post-stack seismic data
It is mutually perpendicular to, time orientation is all vertical with inline direction and xline direction.
For example, the amplitude in described original post-stack seismic data can be inline in described original post-stack seismic data
No. inline 5574 on direction and No. xline 3221 on xline direction are right at 1000 milliseconds of time on time orientation
The amplitude of the time sampling point answered.
Step s102: adjust the arrangement mode of described original post-stack seismic data, obtain the original post-stack seismic data of standard.
Specifically, the dimension of described original post-stack seismic data may include that corresponding time dimension of described time, described
The corresponding inline dimension of main profile number and described cross-track number corresponding cross-track dimension.Described original poststack earthquake number
According to being to be arranged according to the sequence of inline dimension-xline dimension-time dimension.Because acquisition footprint is along isochronous surface
There is the regularity feature such as ribbon, by matrix transpose, the arrangement mode of described original post-stack seismic data can be adjusted to
The sequence of time dimension-inline dimension-xline dimension.Described time dimension can be set to the slowest dimension.Can be by
Original post-stack seismic data after described arrangement mode adjustment is as the original post-stack seismic data of standard.
For example, the process of described matrix transpose can be realized using following manner:
Definition matrix a transposition is n × m rank matrix b, meets b=a (j, i), i.e. b (i, j)=a (j, i) (i-th row jth of b
Column element is jth row i-th column element of a), remember at=b.Two-dimensional matrixTranspositionThree-dimensional array in the same manner
Transposition, keeps one-dimensional constant, in addition bidimensional array realizes transposition according to two-dimensional matrix transposed way.
Step s103: obtain the geological data of very first time section in the original post-stack seismic data of described standard, to first
The geological data of isochronous surface carries out Laplace transform, obtains Laplace transform geological data.
Specifically, original for described standard post-stack seismic data can be transformed into isochronous surface domain by time domain, can obtain
The geological data of corresponding very first time section at the very first time.The geological data of the very first time section obtaining can be entered
Row Laplace transform, obtains Laplace transform geological data, can make adopting in the original post-stack seismic data of described standard
Collection footprint noise highlights in Laplace transform geological data.Wherein, the geological data of very first time section can be
The geological data of corresponding random time section at random time in the original post-stack seismic data of described standard.
For example, as shown in Fig. 2 in Fig. 2 left figure be the geological data of the isochronous surface at the 1000th millisecond.As shown in figure 3,
In Fig. 3, left figure is the Laplace transform after the geological data Laplace transform of isochronous surface at the 1000th millisecond in Fig. 2
Geological data.The geological data of the isochronous surface at can be seen that the 1000th millisecond in left figure from Fig. 3 becomes through Laplce
After changing, banded acquisition footprint noise highlights in described laplace transform geological data, so that collection foot
Print noise obtains Efficient Characterization.
Step s104: described Laplace transform geological data is carried out two-dimensional Fourier transform, obtains Laplce and become
Change geological data wave-number spectrum, according to described Laplace transform geological data wave-number spectrum, obtain the filtering of very first time section because
Son.
Specifically, Laplace transform geological data in left figure in Fig. 3 is carried out two-dimensional Fourier transform, obtain La Pula
This conversion earthquake data wave spectrum number.Laplace transform geological data wave-number spectrum according to right figure in Fig. 3, can obtain institute
State maximum wave number value and the smallest wavenumber value in Laplace transform geological data wave-number spectrum.According to preset rules and maximum ripple
Numerical value and smallest wavenumber value, can adjust the wave number value in described Laplace transform geological data wave-number spectrum.Described default rule
Then can be realized using following formula:
In formula, aiRepresent i-th wave number value in the Laplace transform geological data wave-number spectrum after described adjustment；a0max
Represent the maximum wave number value in the Laplace transform geological data wave-number spectrum before described adjustment；a0minBefore representing described adjustment
Smallest wavenumber value in Laplace transform geological data wave-number spectrum；a0iRepresent the Laplace transform earthquake number before described adjustment
According to i-th wave number value in wave-number spectrum.According to above-mentioned formula, the wave number value of described laplace transform geological data wave-number spectrum is entered
Row adjustment, the filtering factor that the Laplace transform geological data wave-number spectrum after adjustment was cut into slices as the very first time.Described
To the very first time section filtering factor be based on the very first time section geological data obtain.Therefore, each time
The filtering factor of section can be with the collection foot in the geological data of each isochronous surface of identification after Laplace transform
Print the change of noise and change, the acquisition footprint noise in the geological data of each isochronous surface effectively can be pressed
System.
For example: as shown in figure 3, right figure is that in Fig. 3, described in left figure, Laplace transform geological data passes through two dimension in Fig. 3
After Fourier transformation, the Laplace transform geological data wave spectrum number that obtains.According to described formula to the La Pu in right figure in Fig. 3
Lars conversion earthquake data wave spectrum number carry out wave number value adjustment, using adjustment after Laplace transform geological data wave-number spectrum as
The filtering factor of the isochronous surface at the 1000th millisecond.
Step s105: the geological data that the very first time is cut into slices carries out two-dimensional Fourier transform, obtains very first time section
Geological data wave-number spectrum.
Specifically, the geological data of very first time section in original for described standard post-stack seismic data is carried out in two-dimentional Fu
Leaf transformation, obtains the geological data wave-number spectrum of very first time section, so that using the filtering factor in above-mentioned steps in wave-number domain
Footprint compacting is acquired to the geological data of very first time section.
For example, the geological data of the isochronous surface in left figure in Fig. 2 the 1000th millisecond is carried out two-dimensional Fourier transform,
Obtain the geological data wave-number spectrum of the isochronous surface at the 1000th millisecond.In Fig. 2, right figure is after two-dimensional Fourier transform
The wave spectrum number of the geological data of isochronous surface at the 1000th millisecond.
It should be noted that step s105 can be before or after step s103 and/or step s104, the application is to this
Do not define.
Step s106: the filtering factor of section is carried out to the geological data wave-number spectrum that the very first time cuts into slices using the very first time
Acquisition footprint Noise Elimination.
Specifically, the geological data wave-number spectrum that the filtering factor that can be cut into slices the very first time was cut into slices with the very first time is in ripple
Number field is multiplied, and is acquired footprint compacting to the geological data wave-number spectrum of very first time section, obtains described acquisition footprint compacting
The geological data wave-number spectrum of very first time section afterwards.
For example, the geological data of the isochronous surface at the 1000th millisecond after Fig. 4 is the application acquisition footprint Noise Elimination
And its schematic diagram of wave-number spectrum.As shown in figure 4, in Fig. 4 left figure be acquisition footprint Noise Elimination after the 1000th millisecond at when
Between section geological data, in Fig. 4, right figure is the 1000 after the acquisition footprint Noise Elimination after two-dimensional Fourier transform
The geological data wave-number spectrum of the isochronous surface at millisecond.In left figure and Fig. 4 in contrast Fig. 2, left figure is it can be seen that banded adopt
Collection footprint noise can effectively be suppressed.Low, high wave number class acquisition footprint noise is can be seen that in right figure from Fig. 4
To be suppressed.
Fig. 5 be the application acquisition footprint Noise Elimination before and after the 1000th millisecond at isochronous surface geological data difference and
The schematic diagram of its wave-number spectrum.As shown in figure 5, in Fig. 5 left figure be acquisition footprint Noise Elimination before and after the 1000th millisecond at when
Between section geological data poor, in Fig. 5, right figure is the before and after the acquisition footprint Noise Elimination after two-dimensional Fourier transform
The geological data difference wave-number spectrum of the isochronous surface at 1000 milliseconds.From fig. 5, it can be seen that geological information before and after acquisition footprint compacting
It is positively maintained.
In another embodiment, can also be by the geological data of the very first time section after acquisition footprint Noise Elimination
Wave-number spectrum is transformed into spatial domain from wave-number domain, obtains middle post-stack seismic data.Described middle post-stack seismic data can be adjusted
Arrangement mode, obtain target post-stack seismic data.
Specifically, can using above-mentioned embodiment in described standard post-stack seismic data remove the 1000th millisecond at when
Between cut into slices beyond at other times the geological data of corresponding isochronous surface be acquired footprint compacting, obtain acquisition footprint noise
The geological data wave-number spectrum of all isochronous surfaces in described standard post-stack seismic data after compacting.Using two-dimentional Fourier's contravariant
Change, can be by the geological data of all isochronous surfaces in the described standard post-stack seismic data after described acquisition footprint Noise Elimination
Wave-number spectrum is transformed into isochronous surface domain from wave-number domain, then by isochronous surface domain, it is transformed into time domain again, obtains middle folded
Geological data afterwards.Can by matrix transpose by the arrangement mode of described middle post-stack seismic data be adjusted to inline dimension-
The sequence of xline dimension-time dimension, obtains target post-stack seismic data.It is achieved thereby that whole original post-stack seismic data
Acquisition footprint noise compacting.Wherein, in this step in the method for matrix transpose and step s102 matrix transpose method
Identical.
Described earthquake-capturing footprint Noise Elimination embodiment of the method, the Laplace transform based on isochronous surface geological data
The mode combining with wavenumber filtering, is acquired footprint compacting to post-stack seismic data, can be according to marine streamer seismic data
The feature that acquisition footprint changes with room and time, automatic identification simultaneously suppresses acquisition footprint noise, and can effectively retentively
Matter signal.
Fig. 6 is a kind of composition structure chart of the application earthquake-capturing footprint Noise Elimination system embodiment.As shown in fig. 6,
Described earthquake-capturing footprint Noise Elimination system may include that original post-stack seismic data acquiring unit 100, arrangement mode adjustment
Unit 200, Laplace transform unit 300, filtering factor acquiring unit 400, Fourier transform unit 500 and Noise Elimination list
Unit 600.
Described original post-stack seismic data acquiring unit 100, can be used for obtaining original post-stack seismic data.
Described arrangement mode adjustment unit 200, can be used for adjusting the arrangement mode of described original post-stack seismic data, obtains
To the original post-stack seismic data of standard.
Described Laplace transform unit 300, when can be used for obtaining first in the original post-stack seismic data of described standard
Between section geological data.Laplace transform can be carried out to the geological data of very first time section, obtain Laplce and become
Change geological data.
Described filtering factor acquiring unit 400, can be used for for described Laplace transform geological data carrying out two-dimentional Fu
In leaf transformation, obtain Laplace transform geological data wave-number spectrum.According to described Laplace transform geological data wave-number spectrum, can
To obtain the filtering factor of very first time section.
Described second Fourier transform unit 500, can be used for the geological data that the very first time cut into slices and carries out two-dimentional Fu
In leaf transformation, obtain the very first time section geological data wave-number spectrum.
Described Noise Elimination unit 600, can be used for using the filtering factor of very first time section, the very first time being cut into slices
Geological data wave-number spectrum be acquired footprint Noise Elimination.
In another embodiment, described earthquake-capturing footprint Noise Elimination system can also include: middle poststack ground
Seismic data acquisition cell 800 after shake data capture unit 700 and object deck.
Described middle post-stack seismic data acquiring unit 700, can be used for by after acquisition footprint Noise Elimination first when
Between section geological data wave-number spectrum be transformed into spatial domain from wave-number domain, obtain middle post-stack seismic data.
Seismic data acquisition cell 800 after described object deck, can be used for adjusting the row of described middle post-stack seismic data
Row mode, obtains target post-stack seismic data.
Described earthquake-capturing footprint Noise Elimination system embodiment is implemented with described earthquake-capturing footprint Noise Elimination method
Example is corresponding, it is possible to achieve post-stack seismic data is acquired with footprint compacting, can be gathered according to marine streamer seismic data
The feature that footprint changes with room and time, automatic identification simultaneously suppresses acquisition footprint noise, and can effectively keep geology to believe
Number.
In the nineties in 20th century, for a technology improvement can clearly distinguish be improvement on hardware (for example,
Improvement to circuit structures such as diode, transistor, switches) or software on improvement (for the improvement of method flow).So
And, with the development of technology, the improvement of current a lot of method flows can be considered as directly improving of hardware circuit.
Designer is nearly all by being programmed into obtain corresponding hardware circuit in hardware circuit improved method flow.Cause
This is it cannot be said that the improvement of a method flow cannot be realized with hardware entities module.For example, PLD
(programmable logic device, pld) (such as field programmable gate array (field programmable gate
Array, fpga)) it is exactly such a integrated circuit, its logic function to be determined to device programming by user.By designer
Voluntarily programming comes a digital display circuit " integrated " on a piece of pld, without chip maker please designing and make
Special IC chip 2.And, nowadays, replace and manually make IC chip, this programming is also used instead mostly
Realizing, when it is write with program development, software compiler used is similar for " logic compiler (logic compiler) " software
Seemingly, the source code and before will compiling also handy specific programming language writing, this is referred to as hardware description language
(hardware description language, hdl), and hdl is also not only a kind of, but have many kinds, such as abel
(advanced boolean expression language)、ahdl(altera hardware description
language)、confluence、cupl(cornell university programming language)、hdcal、jhdl
(java hardware description language)、lava、lola、myhdl、palasm、rhdl(ruby
Hardware description language) etc., that commonly use most at present is vhdl (very-high-speed
Integrated circuit hardware description language) and verilog2.Those skilled in the art
Should understand it is only necessary to slightly make programming in logic and be programmed into integrated circuit method flow with above-mentioned several hardware description languages
In it is possible to be readily available the hardware circuit realizing this logical method flow process.
Controller can be implemented in any suitable manner, and for example, controller can take such as microprocessor or process
Device and storage can by the computer of the computer readable program code (such as software or firmware) of this (micro-) computing device
Read medium, gate, switch, special IC (application specific integrated circuit,
Asic), the form of programmable logic controller (PLC) and embedded microcontroller, the example of controller includes but is not limited to following microcontroller
Device: arc 625d, atmel at91sam, microchip pic18f26k20 and silicone labs c8051f320, deposit
Memory controller is also implemented as a part for the control logic of memorizer.
It is also known in the art that in addition to controller is realized in pure computer readable program code mode, complete
Full controller can be made with gate, switch, special IC, programmable by method and step is carried out programming in logic
The form of logic controller and embedded microcontroller etc. is realizing identical function.Therefore this controller is considered one kind
Hardware component, and the structure in hardware component can also be considered as to the device for realizing various functions including in it.Or
Even, the device being used for realizing various functions can be considered as not only being the software module of implementation method but also can be Hardware Subdivision
Structure in part.
System, device, module or unit that above-described embodiment illustrates, specifically can be realized by computer chip or entity,
Or to be realized by the product with certain function.
For convenience of description, it is divided into various units to be respectively described with function when describing apparatus above.Certainly, implementing this
The function of each unit can be realized in same or multiple softwares and/or hardware during application.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
Mode by software plus necessary general hardware platform to be realized.Based on such understanding, the technical scheme essence of the application
On in other words prior art is contributed partly can be embodied in the form of software product, typical configure at one
In, computing device includes one or more processors (cpu), input/output interface, network interface and internal memory.This computer is soft
Part product can include some instructions with so that a computer equipment (can be personal computer, server, or network
Equipment etc.) execution each embodiment of the application or embodiment some partly described methods.This computer software product can
To be stored in internal memory, internal memory potentially includes the volatile memory in computer-readable medium, random access memory
(ram) and/or the form such as Nonvolatile memory, such as read only memory (rom) or flash memory (flash ram).Internal memory is computer
The example of computer-readable recording medium.Computer-readable medium includes permanent and non-permanent, removable and non-removable media can be by
Any method or technique is realizing information Store.Information can be computer-readable instruction, data structure, the module of program or its
His data.The example of the storage medium of computer includes, but are not limited to phase transition internal memory (pram), static RAM
(sram), dynamic random access memory (dram), other kinds of random access memory (ram), read only memory
(rom), Electrically Erasable Read Only Memory (eeprom), fast flash memory bank or other memory techniques, read-only optical disc are read-only
Memorizer (cd-rom), digital versatile disc (dvd) or other optical storage, magnetic cassette tape, tape magnetic rigid disk storage or
Other magnetic storage apparatus or any other non-transmission medium, can be used for storing the information that can be accessed by a computing device.According to
Herein defines, and computer-readable medium does not include of short duration computer readable media (transitory media), as modulation
Data signal and carrier wave.
Each embodiment in this specification is all described by the way of going forward one by one, identical similar portion between each embodiment
Divide mutually referring to what each embodiment stressed is the difference with other embodiment.Real especially for system
For applying example, because it is substantially similar to embodiment of the method, so description is fairly simple, referring to embodiment of the method in place of correlation
Part illustrate.
The application can be used in numerous general or special purpose computing system environments or configuration.For example: personal computer, clothes
Business device computer, handheld device or portable set, laptop device, multicomputer system, the system based on microprocessor, put
Top box, programmable consumer-elcetronics devices, network pc, minicomputer, mainframe computer, inclusion any of the above system or equipment
Distributed computing environment etc..
The application can be described in the general context of computer executable instructions, such as program
Module.Usually, program module includes execution particular task or the routine realizing particular abstract data type, program, object, group
Part, data structure etc..The application can also be put into practice in a distributed computing environment, in these distributed computing environment, by
The remote processing devices connected by communication network are executing task.In a distributed computing environment, program module is permissible
It is located in the local and remote computer-readable storage medium including storage device.
Although the application is depicted by embodiment, it will be appreciated by the skilled addressee that the application have many deformation and
Change is without deviating from spirit herein it is desirable to appended claim includes these deformation and change without deviating from the application's
Spirit.
Claims (11)
1. a kind of earthquake-capturing footprint Noise Elimination method is it is characterised in that include:
Obtain original post-stack seismic data；
Adjust the arrangement mode of described original post-stack seismic data, obtain the original post-stack seismic data of standard；
Obtain the geological data of very first time section in the original post-stack seismic data of described standard, the earthquake to very first time section
Data carries out Laplace transform, obtains Laplace transform geological data；
Described Laplace transform geological data is carried out two-dimensional Fourier transform, obtains Laplace transform geological data wave number
Spectrum, according to described Laplace transform geological data wave-number spectrum, obtains the filtering factor of very first time section；
The geological data that the very first time is cut into slices carries out two-dimensional Fourier transform, obtains the geological data wave number of very first time section
Spectrum；
Using the very first time, the filtering factor of section is acquired footprint noise to the geological data wave-number spectrum that the very first time cuts into slices
Compacting.
2. according to claim 1 a kind of earthquake-capturing footprint Noise Elimination method it is characterised in that using described first
After the filtering factor of isochronous surface is acquired footprint Noise Elimination to the geological data wave-number spectrum that the very first time cuts into slices, described
Method also includes:
The geological data wave-number spectrum of the very first time section after acquisition footprint Noise Elimination is transformed into spatial domain from wave-number domain, obtains
To middle post-stack seismic data；
Adjust the arrangement mode of described middle post-stack seismic data, obtain target post-stack seismic data.
3. according to claim 2 a kind of earthquake-capturing footprint Noise Elimination method it is characterised in that being adjusted by matrix transpose
The arrangement mode of whole described middle post-stack seismic data.
4. according to claim 1 a kind of earthquake-capturing footprint Noise Elimination method it is characterised in that drawing described in described basis
Laplace transform geological data wave-number spectrum, obtains the filtering factor of very first time section, comprising:
According to described Laplace transform geological data wave-number spectrum, obtain in described Laplace transform geological data wave-number spectrum
Maximum wave number value and smallest wavenumber value；
According to preset rules and described maximum wave number value and described smallest wavenumber value, adjust described Laplace transform earthquake number
According to the wave number value in wave-number spectrum；
The filtering factor that Laplace transform geological data wave-number spectrum after described adjustment was cut into slices as the very first time.
5. according to claim 4 a kind of earthquake-capturing footprint Noise Elimination method it is characterised in that described basis preset rule
Then and described maximum wave number value and described smallest wavenumber value, adjust the ripple in described Laplace transform geological data wave-number spectrum
Numerical value, comprising:
Realized using following formula:
In formula, aiRepresent i-th wave number value in the Laplace transform geological data wave-number spectrum after described adjustment；a0maxRepresent
The maximum wave number value in Laplace transform geological data wave-number spectrum before described adjustment；a0minRepresent the La Pu before described adjustment
Lars converts the smallest wavenumber value in earthquake data wave number spectrum；a0iRepresent the Laplace transform geological data ripple before described adjustment
I-th wave number value in number spectrum.
6. according to claim 1 a kind of earthquake-capturing footprint Noise Elimination method it is characterised in that described using described
The filtering factor of one isochronous surface is acquired footprint Noise Elimination to the geological data wave-number spectrum that the very first time cuts into slices, comprising:
The geological data wave-number spectrum that the filtering factor that the very first time is cut into slices was cut into slices with the very first time is multiplied in wave-number domain.
7. according to claim 1 a kind of earthquake-capturing footprint Noise Elimination method it is characterised in that described geological data bag
Include: time, amplitude, main profile number and cross-track number.
8. according to claim 7 a kind of earthquake-capturing footprint Noise Elimination method it is characterised in that described adjustment is original folded
The arrangement mode of geological data afterwards, comprising: adjust the arrangement mode of at least two dimensions of described original post-stack seismic data；
Specifically, the time dimension at least two dimensions of described original post-stack seismic data is set to the slowest dimension；Institute
At least two dimensions stating original post-stack seismic data include: described time corresponding time dimension, described main profile correspondence
Inline dimension and described cross-track number corresponding cross-track dimension.
9. according to claim 1 a kind of earthquake-capturing footprint Noise Elimination method it is characterised in that being adjusted by matrix transpose
The arrangement mode of whole described original post-stack seismic data.
10. a kind of earthquake-capturing footprint Noise Elimination system is it is characterised in that include: original post-stack seismic data acquiring unit,
Arrangement mode adjustment unit, Laplace transform unit, filtering factor acquiring unit, Fourier transform unit and Noise Elimination list
Unit；Wherein,
Described original post-stack seismic data acquiring unit, for obtaining original post-stack seismic data；
Described arrangement mode adjustment unit, for adjusting the arrangement mode of described original post-stack seismic data, obtains standard original
Post-stack seismic data；
Described Laplace transform unit, for obtaining the earthquake of very first time section in the original post-stack seismic data of described standard
Data, to the very first time, the geological data of section carries out Laplace transform, obtains Laplace transform geological data；
Described filtering factor acquiring unit, for described Laplace transform geological data is carried out two-dimensional Fourier transform, obtains
To Laplace transform geological data wave-number spectrum, according to described Laplace transform geological data wave-number spectrum, obtain the very first time
The filtering factor of section；
Described Fourier transform unit, carries out two-dimensional Fourier transform for the geological data that the very first time is cut into slices, and obtains
The geological data wave-number spectrum of one isochronous surface；
Described Noise Elimination unit, the geological data ripple that the filtering factor for being cut into slices using the very first time was cut into slices to the very first time
Number spectrum is acquired footprint Noise Elimination.
11. according to claim 10 a kind of earthquake-capturing footprint Noise Elimination system it is characterised in that described system is also wrapped
Include: seismic data acquisition cell after middle post-stack seismic data acquiring unit and object deck；Wherein,
Described middle post-stack seismic data acquiring unit, for the earthquake that the very first time after acquisition footprint Noise Elimination cuts into slices
Data wave-number spectrum is transformed into spatial domain from wave-number domain, obtains middle post-stack seismic data；
Seismic data acquisition cell after described object deck, for adjusting the arrangement mode of described middle post-stack seismic data, obtains
Target post-stack seismic data.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610648072.6A CN106353816B (en) | 2016-08-09 | 2016-08-09 | A kind of earthquake-capturing footprint Noise Elimination method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610648072.6A CN106353816B (en) | 2016-08-09 | 2016-08-09 | A kind of earthquake-capturing footprint Noise Elimination method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106353816A true CN106353816A (en) | 2017-01-25 |
| CN106353816B CN106353816B (en) | 2018-10-16 |
Family
ID=57844746
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610648072.6A Active CN106353816B (en) | 2016-08-09 | 2016-08-09 | A kind of earthquake-capturing footprint Noise Elimination method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106353816B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108445540A (en) * | 2018-03-13 | 2018-08-24 | 林腾飞 | A kind of auxiliary foot print drawing method |
| CN111599234A (en) * | 2020-05-19 | 2020-08-28 | 黑龙江工业学院 | Automatic English spoken language scoring system based on voice recognition |
| CN112255684A (en) * | 2020-10-26 | 2021-01-22 | 中国石油天然气集团有限公司 | Processing method and device for collecting footprint of seismic data |
| CN113640882A (en) * | 2021-08-10 | 2021-11-12 | 南方海洋科学与工程广东省实验室(湛江) | Method for removing noise of collected footprint, electronic device and computer readable storage medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101852866A (en) * | 2010-05-13 | 2010-10-06 | 中国石油天然气股份有限公司 | Post-stack seismic data filtering method |
| US20130345985A1 (en) * | 2012-06-25 | 2013-12-26 | Schlumberger Technology Corporation | Seismic orthogonal decomposition attribute |
| CN103558636A (en) * | 2013-11-08 | 2014-02-05 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for collecting attenuation of footprints from post-stack seismic data |
| CN105425294A (en) * | 2015-11-30 | 2016-03-23 | 中国石油天然气集团公司 | Determination method and apparatus for acquisition footprint |
-
2016
- 2016-08-09 CN CN201610648072.6A patent/CN106353816B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101852866A (en) * | 2010-05-13 | 2010-10-06 | 中国石油天然气股份有限公司 | Post-stack seismic data filtering method |
| US20130345985A1 (en) * | 2012-06-25 | 2013-12-26 | Schlumberger Technology Corporation | Seismic orthogonal decomposition attribute |
| CN103558636A (en) * | 2013-11-08 | 2014-02-05 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for collecting attenuation of footprints from post-stack seismic data |
| CN105425294A (en) * | 2015-11-30 | 2016-03-23 | 中国石油天然气集团公司 | Determination method and apparatus for acquisition footprint |
Non-Patent Citations (2)
| Title |
|---|
| ZOU ZHEN EL.: ""Acquisition footprint suppression of towed streamer 4D seismic data"", 《2015 SEG NEW ORLEANS ANNUAL MEETING》 * |
| 麻三怀 等: ""采集脚印分析和处理方法综述"", 《地球物理学进展》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108445540A (en) * | 2018-03-13 | 2018-08-24 | 林腾飞 | A kind of auxiliary foot print drawing method |
| CN111599234A (en) * | 2020-05-19 | 2020-08-28 | 黑龙江工业学院 | Automatic English spoken language scoring system based on voice recognition |
| CN112255684A (en) * | 2020-10-26 | 2021-01-22 | 中国石油天然气集团有限公司 | Processing method and device for collecting footprint of seismic data |
| CN112255684B (en) * | 2020-10-26 | 2024-03-01 | 中国石油天然气集团有限公司 | Processing method and device for acquisition footprint of seismic data |
| CN113640882A (en) * | 2021-08-10 | 2021-11-12 | 南方海洋科学与工程广东省实验室(湛江) | Method for removing noise of collected footprint, electronic device and computer readable storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106353816B (en) | 2018-10-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Chen et al. | Empirical low-rank approximation for seismic noise attenuation | |
| Zhang et al. | Hybrid rank-sparsity constraint model for simultaneous reconstruction and denoising of 3D seismic data | |
| Chen et al. | Robust reduced-rank filtering for erratic seismic noise attenuation | |
| Kreimer et al. | Tensor completion based on nuclear norm minimization for 5D seismic data reconstruction | |
| Trorey | Diffractions for arbitrary source-receiver locations | |
| Margrave et al. | Gabor deconvolution: Estimating reflectivity by nonstationary deconvolution of seismic data | |
| Gao et al. | Parallel matrix factorization algorithm and its application to 5D seismic reconstruction and denoising | |
| CN106353816B (en) | A kind of earthquake-capturing footprint Noise Elimination method and system | |
| Krischer et al. | Large‐scale seismic inversion framework | |
| CN107193045B (en) | A kind of seismic data processing technique and device | |
| Liu et al. | OC-seislet: Seislet transform construction with differential offset continuation | |
| Gholami | Deconvolutive radon transform | |
| Mansour et al. | Improved wavefield reconstruction from randomized sampling via weighted one-norm minimization | |
| Carozzi et al. | Robust tensor-completion algorithm for 5D seismic-data reconstruction | |
| Elison et al. | Data-driven control over short-period internal multiples in media with a horizontally layered overburden | |
| Kaplan et al. | Data reconstruction with shot-profile least-squares migration | |
| CN111239808A (en) | Seismic inversion spectrum decomposition method and device, electronic equipment and storage medium | |
| CN106226812A (en) | The method and device of seismic data resolution is improved based on time-frequency Quadratic Spectrum | |
| CN104199083A (en) | Method and device for determining wide-line seismological observation system | |
| CN106353818A (en) | Seismic data time difference correcting method and system | |
| Innanen | Time-and offset-domain internal multiple prediction with nonstationary parameters | |
| CN108897036B (en) | Seismic data processing method and device | |
| CN107561588B (en) | Seismic data noise suppression method and device | |
| CN106291691A (en) | A kind of seismic migration imaging method and device | |
| Oboué et al. | Enhanced low-rank matrix estimation for simultaneous denoising and reconstruction of 5D seismic data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | 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 |