CN104375185A - Method and device for removing surface waves from seismic records - Google Patents
Method and device for removing surface waves from seismic records Download PDFInfo
- Publication number
- CN104375185A CN104375185A CN201410612211.0A CN201410612211A CN104375185A CN 104375185 A CN104375185 A CN 104375185A CN 201410612211 A CN201410612211 A CN 201410612211A CN 104375185 A CN104375185 A CN 104375185A
- Authority
- CN
- China
- Prior art keywords
- seismologic record
- ground roll
- frequency range
- record
- seismologic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The embodiment of the invention discloses a method and device for removing surface waves from seismic records. The method includes the steps that frequency division scanning is conducted on the seismic records to be processed, and the distribution range of surface wave frequency bands is acquired; according to the distribution range of the surface wave frequency bands, wavelet frequency division is conducted on the seismic records, and the seismic records on different frequency bands are acquired; according to the distribution range of the surface wave frequency bands, the seismic records containing the surface waves are recognized from the seismic records on the different frequency bands, and radial way conversion is conducted on the seismic records containing the surface waves; according to the distribution range of the surface wave frequency bands, the high-pass filtering frequency is determined, and high-pass filtering is conducted on the seismic records acquired after radial way conversion in a radial way domain; Reverse radial way conversion is conducted on the seismic records acquired after high-pass filtering; wavelet reconstruction is conducted on the seismic records acquired after reverse radial way conversion and the seismic records on the different frequency bands and not containing the surface waves, and the seismic records with the surface waves removed are acquired.
Description
Technical field
The application relates to petroleum prospecting Seismic Data Processing Technique field, particularly relates to a kind of method and apparatus removing ground roll from seismologic record.
Background technology
Because ground roll exists extensively, energy is strong, have a strong impact on the usable reflection of mid-deep strata, reduce the signal to noise ratio (S/N ratio) of seismic data, so surface wave eliminations becomes an important process in seism processing, for the feature of ground roll, the method for multiple removal ground roll can be taked, as the method such as radial trace transform, wavelet transformation, but the method for these single elimination ground rolls does not often reach desirable effect.
Wavelet transformation is that seismologic record is carried out wavelet decomposition on multiple different yardstick, the wavelet conversion coefficient that each yardstick obtains through wavelet decomposition just represents original place and shakes the information be recorded in different time or frequency, so just can utilize ground roll and the significant wave difference in frequency, ground roll and significant wave are made a distinction, and then removal ground roll, but the ground roll utilizing wavelet transformation to remove has residual, and noise removal capability is weak.
Radial trace transform is a kind of simple geometry seismic trace mapping algorithm, because ground roll and seismologic record exist huge difference in speed, the method is by the conversion of coordinate system, ground roll and effective seismologic record are made a distinction, and then removal ground roll, but have the effective seismologic record of part and be also removed.
Summary of the invention
The object of the application is to provide a kind of method and apparatus removing ground roll from seismologic record, in order to when reducing seismologic record loss, effectively removes ground roll.
To achieve these goals, this application provides a kind of method removing ground roll from seismologic record, the method comprises:
S101: carry out fraction frequency to pending seismologic record, obtains the distribution range of ground roll frequency range;
S102: according to the distribution range of described ground roll frequency range, small echo frequency division is carried out to described seismologic record, obtain the seismologic record on different frequency range, the seismologic record on described different frequency range comprises the seismologic record comprising ground roll and the seismologic record not comprising ground roll section;
S103: according to the distribution range of described ground roll frequency range, identify the seismologic record comprising ground roll from the seismologic record described different frequency range, carries out radial trace transform to the described seismologic record comprising ground roll;
S104: according to the frequency of the distribution range determination high-pass filtering of described ground roll frequency range, in territory, radial road, high-pass filtering is carried out to the seismologic record after radial trace transform in S103;
S105: anti-radial trace transform is carried out to the seismologic record after high-pass filtering in S104;
S106: the seismologic record after radial trace transform anti-in S105 and the seismologic record not comprising ground roll described in S102 are carried out wavelet reconstruction, obtains the seismologic record removing ground roll.
In a preferred embodiment, described seismologic record have recorded multiple tracks sampled point, and per pass has multiple sampled point, and described sampled point have recorded the amplitude of signal.
In a preferred embodiment, describedly carry out fraction frequency to described seismologic record, the distribution range obtaining ground roll frequency range comprises:
Spectrum analysis is carried out to described seismologic record, obtains the band limits of seismologic record;
Described seismologic record is carried out frequency range division, obtains the seismologic record of different frequency ranges;
In each band limits, frequency filtering is carried out to described seismologic record, obtains the seismologic record of the time domain of different frequency range;
According to the seismologic record determination ground roll of the time domain of described different frequency range, obtain the distribution range of ground roll frequency range.
In a preferred embodiment, the described distribution range according to described ground roll frequency range, carries out small echo frequency division to described seismologic record, and the seismologic record obtained on different frequency range comprises:
The division range of described small echo frequency division is set according to the distribution range of ground roll frequency range;
Horse traction spy (Mallat) algorithm is adopted to carry out frequency division to described seismologic record, seismologic record is resolved into seismologic record wavelet details and seismologic record wavelets approximation, and the seismologic record wavelets approximation decomposed the last time is as inputting next time, obtain the seismologic record on different frequency range, formula is as follows:
In above formula, High, Low are pair of orthogonal mirror filter; J decomposes number of times, j=0,1 ... j-1, j=0 represent the 1st decomposition; 2
janalyze yardstick,
jth+2 secondary earthquake record wavelet details,
jth+1 secondary earthquake record wavelets approximation,
it is jth+2 secondary earthquake record wavelets approximation.
In a preferred embodiment, describedly radial trace transform carried out to the described seismologic record comprising ground roll comprise:
The seismologic record comprising ground roll described in t-x territory is transformed to τ-v territory, and utilizes cubic spline interpolation to realize radial trace transform.As follows from the mapping relations transforming to τ-v territory in t-x territory:
t→τ＝t-t
0
x→v＝(x-x
0)/(t-t
0)
In above formula, t is two-way time; X is offset distance; x
0for the origin position of radial trace transform; t
0for the two-way time of the origin position of radial trace transform; τ is Delay time; V is speed.
In a preferred embodiment, describedly anti-radial trace transform carried out to the seismologic record after high-pass filtering in S104 comprise:
Seismologic record after the described high-pass filtering in τ-v territory is transformed to t-x territory, as follows from the mapping relations transforming to τ-v territory in t-x territory:
τ→t＝τ+t
0
v→x＝v(t-t
0)+x
0
In above formula, t is two-way time; X is offset distance; x
0for the origin position of radial trace transform; t
0for the two-way time of the origin position of radial trace transform; τ is Delay time; V is speed.
The application also provides a kind of device removing ground roll from seismologic record on the other hand, and this device comprises:
Ground roll frequency range determining unit, for carrying out fraction frequency to pending seismologic record, obtains the distribution range of ground roll frequency range;
Frequency unit, for the distribution range according to described ground roll frequency range, carry out small echo frequency division to described seismologic record, obtain the seismologic record on different frequency range, the seismologic record on described different frequency range comprises the seismologic record comprising ground roll and the seismologic record not comprising ground roll section;
Radial trace transform unit, for the distribution range according to described ground roll frequency range, identifies the seismologic record comprising ground roll from the seismologic record described different frequency range, carries out radial trace transform to the described seismologic record comprising ground roll;
Filter unit, for the frequency of the distribution range determination high-pass filtering according to described ground roll frequency range, the seismologic record after the radial trace transform obtain radial road converter unit carries out high-pass filtering in territory, radial road;
Anti-radial trace transform unit, carries out anti-radial trace transform for the seismologic record after the high-pass filtering that obtains filter unit;
Wavelet reconstruction unit, carries out wavelet reconstruction for the seismologic record not comprising ground roll described in the seismologic record after the anti-radial trace transform obtained in anti-radial trace transform unit and frequency unit being obtained, and obtains the seismologic record removing ground roll.
The technical scheme provided from above the embodiment of the present application, the application is by carrying out fraction frequency to pending seismologic record, obtain the distribution range of ground roll frequency range, according to the distribution range of ground roll frequency range, small echo frequency division is carried out to pending seismologic record, obtain the seismologic record of different frequency range, the seismologic record comprising ground roll is identified from the seismologic record of different frequency range, carry out radial trace transform, ground roll and effective earthquake is utilized to be recorded in difference in speed, High frequency filter is carried out in territory, radial road, again filtered seismologic record is carried out radial road inverse transformation, the seismologic record not comprising ground roll in the seismologic record of the different frequency range obtained after the seismologic record after radial road inverse transformation and above-mentioned small echo frequency division is carried out wavelet reconstruction, obtain the seismologic record removing ground roll.Compared with prior art, not removing only ground roll, and greatly reduce the loss of effective seismologic record.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present application or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, the accompanying drawing that the following describes is only some embodiments recorded in the application, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.
Fig. 1 is a kind of process flow diagram removing the method for ground roll from seismologic record that the embodiment of the present application provides;
Fig. 2 is pending seismologic record sectional view;
Fig. 3 is the seismologic record sectional view that the seismologic record in Fig. 2 is obtained by the method that small echo frequency division removes ground roll;
Fig. 4 is the seismologic record sectional view that the seismologic record in Fig. 2 is obtained by the method that radial trace transform removes ground roll;
Fig. 5 is the seismologic record sectional view that the seismologic record in Fig. 2 provides a kind of method removing ground roll from seismologic record to obtain by the embodiment of the present application;
Fig. 6 is a kind of schematic diagram removing the device of ground roll from seismologic record that the embodiment of the present application provides.
Embodiment
Technical scheme in the application is understood better in order to make those skilled in the art person, below in conjunction with the accompanying drawing in the embodiment of the present application, technical scheme in the embodiment of the present application is clearly and completely described, obviously, described embodiment is only some embodiments of the present application, instead of whole embodiments.Based on the embodiment in the application, those of ordinary skill in the art are not making the every other embodiment obtained under creative work prerequisite, all should belong to the scope of the application's protection.
The specific implementation of the embodiment of the present application is described below in detail with concrete example.
As shown in Figure 1, a kind of method removing ground roll from seismologic record that the embodiment of the present application provides, the method comprises:
S201: carry out fraction frequency to pending seismologic record, obtains the distribution range of ground roll frequency range.
In practice, obtaining seismologic record needs to excite by blowing out on ground, and surface geophone receives the seismologic record reflected, and this seismologic record is as pending seismologic record.Above-mentioned seismologic record have recorded multiple tracks sampled point, and per pass has multiple sampled point, and described sampled point have recorded the amplitude of signal.
The seismologic record collected first is carried out Fourier transform, according to the frequency spectrum of seismologic record after Fourier transform, obtains the band limits of seismologic record; Then, frequency range division is carried out to seismologic record, obtain the seismologic record of different frequency ranges; And according to the division of frequency range, a wave filter is designed to each frequency range; Finally, the seismologic record of each frequency range of carrying out after Fourier transform is multiplied with each wave filter, then inversefouriertransform is carried out to each multiplied result; The seismologic record of the time domain of different frequency range can be obtained, from the seismologic record of the time domain of different frequency range, determine ground roll place frequency range, obtain the distribution range of ground roll frequency range.
S202: according to the distribution range of above-mentioned ground roll frequency range, carry out small echo frequency division to seismologic record, obtain the seismologic record on different frequency range, the seismologic record on different frequency range comprises the seismologic record comprising ground roll and the seismologic record not comprising ground roll.
As shown in Figure 2 be pending seismologic record sectional view, wherein, horizontal ordinate is number of channels, and ordinate is the time, and unit is ms.The seismologic record sectional view of the part in the 1220th road, road to 1440 is have recorded in figure.
In practice, the division range of described small echo frequency division is set according to the distribution range of ground roll frequency range; Such as when the distribution range of ground roll frequency range is 6-15HZ, can directly by the setting of division range when carrying out small echo frequency division to seismologic record, the part seismologic record of 6-15HZ frequency range is branched away, as the seismologic record comprising ground roll, and other be the seismologic record not comprising ground roll.
Small echo frequency division is carried out by horse traction spy (Mallat) algorithm realization to seismologic record, is equivalent to seismologic record to resolve into seismologic record wavelet details and seismologic record wavelets approximation.Wavelet details, the namely wavelet transformation of our certain usually said yardstick, the result of its in fact high-pass filtering; Wavelets approximation, it is actually the result of low-pass filtering.The seismologic record wavelets approximation decomposed last time, as inputting next time, removes the seismologic record wavelet details and the seismologic record wavelets approximation that obtain thicker yardstick.Like this two points several times, until remaining low frequency signal void value obtains and can divide, and then obtain the seismologic record on different frequency range, formula is as follows:
In above formula, High, Low are pair of orthogonal mirror filter; J decomposes number of times, j=0,1 ... j-1, j=0 represent the 1st decomposition; 2
janalyze yardstick,
jth+2 secondary earthquake record wavelet details,
jth+1 secondary earthquake record wavelets approximation,
it is jth+2 secondary earthquake record wavelets approximation.
It will be apparent to those skilled in the art that and to be not limited only to by a kind of mode of horse traction spy (Mallat) algorithm the frequency division mode of seismologic record for realization, can also be other forms in practice, the embodiment of the present application be as limit.
S203: according to the distribution range of ground roll frequency range, identifies the seismologic record comprising ground roll frequency range from the seismologic record different frequency range, carries out radial trace transform to the seismologic record comprising ground roll frequency range.
In practice, owing to when carrying out small echo frequency division to seismologic record being the distribution range according to ground roll frequency range, so the band limits comprising the seismologic record of ground roll is exactly the distribution range of ground roll frequency range, therefore can directly according to the distribution range of ground roll frequency range, identify the seismologic record comprising ground roll in the seismologic record of the different frequency range obtained from S202, this seismologic record comprising ground roll comprises ground roll and effective seismologic record.
Owing to there is larger difference in ground roll and effective seismologic record speed, face wave propagation velocity is less, and the formation velocity specific surface wave propagation velocity that seismologic record records is much larger, so carry out radial trace transform to the seismologic record identified, the seismologic record comprising ground roll in t-x territory is transformed to τ-v territory, ground roll and effective earthquake recording areas are separated.Owing to being over-sampling in the position close to t-x territory initial point, and be down-sampling away from the position of initial point, therefore radial trace transform needs interpolation.The embodiment of the present application utilizes cubic spline interpolation to realize radial trace transform.
As follows from the mapping relations transforming to τ-v territory in t-x territory:
t→τ＝t-t
0
x→v＝(x-x
0)/(t-t
0)
In above formula, t is two-way time; X is offset distance; x
0for the origin position of radial trace transform; t
0for the two-way time of the origin position of radial trace transform; τ is Delay time; V is speed.
It will be apparent to those skilled in the art that the interpolation method for radial trace transform is not limited only to be realized by a kind of mode of cubic spline interpolation, can also be other forms in practice, the embodiment of the present application be as limit.
S204: according to the frequency of the distribution range determination high-pass filtering of ground roll frequency range, in territory, radial road, high-pass filtering is carried out to the seismologic record after radial trace transform in S203.
In practice, high-pass filtering to be carried out to the seismologic record after radial trace transform in territory, radial road, remove ground roll, the frequency of high-pass filtering need be first set.The frequency of high-pass filtering is exactly the highest frequency of ground roll frequency range, the distribution range of such as ground roll frequency range is 6-15HZ, then the frequency of high-pass filtering can be set to 15HZ, due to the speed of ground roll and effective seismologic record having larger difference, therefore carry out high-pass filtering in territory, radial road after the seismologic record comprising ground roll is done radial trace transform, can ground roll be removed.
S205: anti-radial trace transform is carried out to the seismologic record after high-pass filtering in S204.
Seismologic record after high-pass filtering is carried out anti-radial trace transform, the seismologic record after the described high-pass filtering in τ-v territory is transformed to t-x territory, as follows from the mapping relations transforming to τ-v territory in t-x territory:
τ→t＝τ+t
0
v→x＝v(t-t
0)+x
0
In above formula, t is two-way time; X is offset distance; x
0for the origin position of radial trace transform; t
0for the two-way time of the origin position of radial trace transform; τ is Delay time; V is speed.
S206: the seismologic record after radial trace transform anti-in S205 and the seismologic record not comprising ground roll in S202 are carried out wavelet reconstruction, obtains the seismologic record removing ground roll.
Seismologic record in above-mentioned S205 after anti-radial trace transform is the seismologic record after the seismologic record removal ground roll comprising ground roll in the seismologic record of the different frequency range obtained in above-mentioned S202, the seismologic record not comprising ground roll in the seismologic record of the different frequency range obtained in the seismologic record after this anti-radial trace transform and above-mentioned S202 is carried out wavelet reconstruction, finally obtains the seismologic record removing ground roll.Namely the process of wavelet reconstruction is passed through
with
reconstruct obtains
also namely reconstructed by jth+2 secondary earthquake record wavelet details in wavelet decomposition process and jth+2 secondary earthquake record wavelets approximation and obtain jth+1 secondary earthquake record wavelets approximation.
Below by way of the scheme of comparative illustration the embodiment of the present application relative to the technique effect of prior art:
As shown in Figure 3 be the seismologic record sectional view that method that seismologic record in Fig. 2 directly removes ground roll by small echo frequency division obtains, wherein, horizontal ordinate is number of channels, and ordinate is the time, and unit is ms.Find to eliminate part ground roll through too small frequency division from figure, but there is frequency dispersion due to ground roll, in nearly geophone offset there also residual.
As shown in Figure 4 be the seismologic record sectional view that method that seismologic record in Fig. 2 directly removes ground roll by radial trace transform obtains, wherein, horizontal ordinate is number of channels, and ordinate is the time, and unit is ms.Find effectively to eliminate ground roll through radial trace transform from figure, but also have lost effective seismologic record of deep layer.
As shown in Figure 5 be the seismologic record sectional view that the seismologic record in Fig. 2 provides a kind of method removing ground roll from seismologic record to obtain by the embodiment of the present application, wherein, horizontal ordinate is number of channels, and ordinate is the time, and unit is ms.From figure, find the seismologic record providing a kind of method removing ground roll from seismologic record to obtain by the embodiment of the present application, not only effectively eliminate ground roll, and effectively reduce the loss of seismologic record.
As can be seen here, a kind of method removing ground roll from seismologic record that the embodiment of the present application provides, by carrying out fraction frequency to pending seismologic record, obtain the distribution range of ground roll frequency range, according to the distribution range of ground roll frequency range, small echo frequency division is carried out to pending seismologic record, obtain the seismologic record of different frequency range, the seismologic record comprising ground roll is identified from the seismologic record of different frequency range, carry out radial trace transform, ground roll and effective earthquake is utilized to be recorded in difference in speed, High frequency filter is carried out in territory, radial road, again filtered seismologic record is carried out radial road inverse transformation, the seismologic record not comprising ground roll in the seismologic record of the different frequency range obtained after the seismologic record after radial road inverse transformation and above-mentioned small echo frequency division is carried out wavelet reconstruction, obtain the seismologic record removing ground roll.Compared with prior art, not removing only ground roll, and greatly reduce the loss of effective seismologic record.
The embodiment of the present application is also corresponding provides a kind of device 300 removing ground roll from seismologic record, as shown in Figure 6, this device 300 comprises: ground roll frequency range determining unit 301, frequency unit 302, radial trace transform unit 303, filter unit 304, anti-radial trace transform unit 305 and wavelet reconstruction unit 306.
Wherein, ground roll frequency range determining unit 301, for carrying out fraction frequency to pending seismologic record, obtains the distribution range of ground roll frequency range; Frequency unit 302, for the distribution range according to ground roll frequency range, carries out small echo frequency division to seismologic record, and obtain the seismologic record on different frequency range, the seismologic record on this different frequency range comprises the seismologic record comprising ground roll and the seismologic record not comprising ground roll section; Radial trace transform unit 303, for the distribution range according to above-mentioned ground roll frequency range, identifies the seismologic record comprising ground roll from the seismologic record above-mentioned different frequency range, carries out radial trace transform to the seismologic record comprising ground roll; Filter unit 304 is for the frequency of the distribution range determination high-pass filtering according to above-mentioned ground roll frequency range, and the seismologic record after the radial trace transform obtain radial road converter unit 303 carries out high-pass filtering in territory, radial road; Anti-radial trace transform unit 305 carries out anti-radial trace transform for the seismologic record after the high-pass filtering that obtains filter unit 304; Wavelet reconstruction unit 306, for the seismologic record after the anti-radial trace transform obtained in anti-radial trace transform unit 305 and frequency unit 302 are obtained the seismologic record not comprising ground roll carry out wavelet reconstruction, obtain removing the seismologic record of ground roll.
From with, a kind of device removing ground roll from seismologic record that the embodiment of the present application provides, by carrying out fraction frequency to pending seismologic record, obtain the distribution range of ground roll frequency range, according to the distribution range of ground roll frequency range, small echo frequency division is carried out to pending seismologic record, obtain the seismologic record of different frequency range, the seismologic record comprising ground roll is identified from the seismologic record of different frequency range, carry out radial trace transform, ground roll and effective earthquake is utilized to be recorded in difference in speed, High frequency filter is carried out in territory, radial road, again filtered seismologic record is carried out radial road inverse transformation, the seismologic record not comprising ground roll in the seismologic record of the different frequency range obtained after the seismologic record after radial road inverse transformation and above-mentioned small echo frequency division is carried out wavelet reconstruction, obtain the seismologic record removing ground roll.Compared with prior art, not removing only ground roll, and greatly reduce the loss of effective seismologic record.
Each embodiment in this instructions all adopts the mode of going forward one by one to describe, and what each embodiment stressed is the difference with other embodiments, between each embodiment identical similar part mutually see.Especially, for system embodiment, because it is substantially similar to embodiment of the method, so description is fairly simple, relevant part illustrates see the part of embodiment of the method.
Although depict the application by embodiment, those of ordinary skill in the art know, the application has many distortion and change and do not depart from the spirit of the application, and the claim appended by wishing comprises these distortion and change and do not depart from the spirit of the application.
Claims (7)
1. from seismologic record, remove a method for ground roll, it is characterized in that, the method comprises:
S1: carry out fraction frequency to pending seismologic record, obtains the distribution range of ground roll frequency range;
S2: according to the distribution range of described ground roll frequency range, small echo frequency division is carried out to described seismologic record, obtain the seismologic record on different frequency range, the seismologic record on described different frequency range comprises the seismologic record comprising ground roll and the seismologic record not comprising ground roll section;
S3: according to the distribution range of described ground roll frequency range, identify the seismologic record comprising ground roll from the seismologic record described different frequency range, carries out radial trace transform to the described seismologic record comprising ground roll;
S4: according to the frequency of the distribution range determination high-pass filtering of described ground roll frequency range, in territory, radial road, high-pass filtering is carried out to the seismologic record after radial trace transform in S3;
S5: anti-radial trace transform is carried out to the seismologic record after high-pass filtering in S4;
S6: the seismologic record after radial trace transform anti-in S5 and the seismologic record not comprising ground roll described in S2 are carried out wavelet reconstruction, obtains the seismologic record removing ground roll.
2. the method removing ground roll from seismologic record according to claim 1, is characterized in that, described seismologic record have recorded multiple tracks sampled point, and per pass has multiple sampled point, and described sampled point have recorded the amplitude of signal.
3. the method removing ground roll from seismologic record according to claim 1, is characterized in that, describedly carries out fraction frequency to described seismologic record, and the distribution range obtaining ground roll frequency range comprises:
Spectrum analysis is carried out to described seismologic record, obtains the band limits of seismologic record;
Described seismologic record is carried out frequency range division, obtains the seismologic record of different frequency ranges;
In each band limits, frequency filtering is carried out to described seismologic record, obtains the seismologic record of the time domain of different frequency range;
According to the seismologic record determination ground roll of the time domain of described different frequency range, obtain the distribution range of ground roll frequency range.
4. the method removing ground roll from seismologic record according to claim 1, is characterized in that, the described distribution range according to described ground roll frequency range, carries out small echo frequency division to described seismologic record, and the seismologic record obtained on different frequency range comprises:
The division range of described small echo frequency division is set according to the distribution range of ground roll frequency range;
Horse traction spy (Mallat) algorithm is adopted to carry out frequency division to described seismologic record, seismologic record is resolved into seismologic record wavelet details and seismologic record wavelets approximation, and the seismologic record wavelets approximation decomposed the last time is as inputting next time, obtain the seismologic record on different frequency range, formula is as follows:
In above formula, High, Low are pair of orthogonal mirror filter; J decomposes number of times, j=0,1 ... j-1, j=0 represent the 1st decomposition; 2
janalyze yardstick,
jth+2 secondary earthquake record wavelet details,
jth+1 secondary earthquake record wavelets approximation,
it is jth+2 secondary earthquake record wavelets approximation.
5. the method removing ground roll from seismologic record according to claim 4, is characterized in that, describedly carries out radial trace transform to the described seismologic record comprising ground roll and comprises:
The seismologic record comprising ground roll described in t-x territory is transformed to τ-v territory, and utilizes cubic spline interpolation to realize radial trace transform; As follows from the mapping relations transforming to τ-v territory in t-x territory:
t→τ＝t-t
0
x→v＝(x-x
0)/(t-t
0)
In above formula, t is two-way time; X is offset distance; x
0for the origin position of radial trace transform; t
0for the two-way time of the origin position of radial trace transform; τ is Delay time; V is speed.
6. the method removing ground roll from seismologic record according to claim 5, is characterized in that, describedly carries out anti-radial trace transform to the seismologic record after high-pass filtering in S4 and comprises:
Seismologic record after the described high-pass filtering in τ-v territory is transformed to t-x territory, as follows from the mapping relations transforming to τ-v territory in t-x territory:
τ→t＝τ+t
0
v→x＝v(t-t
0)+x
0
In above formula, t is two-way time; X is offset distance; x
0for the origin position of radial trace transform; t
0for the two-way time of the origin position of radial trace transform; τ is Delay time; V is speed.
7. from seismologic record, remove a device for ground roll, it is characterized in that, this device comprises:
Ground roll frequency range determining unit, for carrying out fraction frequency to pending seismologic record, obtains the distribution range of ground roll frequency range;
Frequency unit, for the distribution range according to described ground roll frequency range, carry out small echo frequency division to described seismologic record, obtain the seismologic record on different frequency range, the seismologic record on described different frequency range comprises the seismologic record comprising ground roll and the seismologic record not comprising ground roll section;
Radial trace transform unit, for the distribution range according to described ground roll frequency range, identifies the seismologic record comprising ground roll from the seismologic record described different frequency range, carries out radial trace transform to the described seismologic record comprising ground roll;
Filter unit, for the frequency of the distribution range determination high-pass filtering according to described ground roll frequency range, the seismologic record after the radial trace transform obtain radial road converter unit carries out high-pass filtering in territory, radial road;
Anti-radial trace transform unit, carries out anti-radial trace transform for the seismologic record after the high-pass filtering that obtains filter unit;
Wavelet reconstruction unit, carries out wavelet reconstruction for the seismologic record not comprising ground roll described in the seismologic record after the anti-radial trace transform obtained in anti-radial trace transform unit and frequency unit being obtained, and obtains the seismologic record removing ground roll.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410612211.0A CN104375185A (en) | 2014-11-04 | 2014-11-04 | Method and device for removing surface waves from seismic records |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410612211.0A CN104375185A (en) | 2014-11-04 | 2014-11-04 | Method and device for removing surface waves from seismic records |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN104375185A true CN104375185A (en) | 2015-02-25 |
Family
ID=52554213
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410612211.0A Pending CN104375185A (en) | 2014-11-04 | 2014-11-04 | Method and device for removing surface waves from seismic records |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104375185A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106950600A (en) * | 2017-02-16 | 2017-07-14 | 中国石油大学(华东) | A kind of minimizing technology of near surface scattering surface ripple |
| CN109239780A (en) * | 2017-07-10 | 2019-01-18 | 中国石油化工股份有限公司 | Based on the synchronous method for squeezing wavelet transformation removal surface wave |
| CN113687420A (en) * | 2020-05-18 | 2021-11-23 | 中国石油天然气股份有限公司 | Method and device for removing surface wave noise of seismic data |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101915939A (en) * | 2010-08-17 | 2010-12-15 | 中国科学院地质与地球物理研究所 | Surface wave pressing method |
| US20140016436A1 (en) * | 2012-07-10 | 2014-01-16 | Pgs Geophysical As | Methods and systems for reconstruction of low frequency particle velocity wavefields and deghosting of seismic streamer data |
-
2014
- 2014-11-04 CN CN201410612211.0A patent/CN104375185A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101915939A (en) * | 2010-08-17 | 2010-12-15 | 中国科学院地质与地球物理研究所 | Surface wave pressing method |
| US20140016436A1 (en) * | 2012-07-10 | 2014-01-16 | Pgs Geophysical As | Methods and systems for reconstruction of low frequency particle velocity wavefields and deghosting of seismic streamer data |
Non-Patent Citations (4)
| Title |
|---|
| 刘志鹏等: "径向道变换压制相干噪声方法研究", 《地球物理学进展》 * |
| 夏洪瑞等: "小波变换及其在去噪中的应用", 《石油地球物理勘探》 * |
| 张军华等: "地震资料去噪方法技术综合评述", 《地球物理学进展》 * |
| 张少梅: "地震资料面波正演模拟及去噪技术研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106950600A (en) * | 2017-02-16 | 2017-07-14 | 中国石油大学(华东) | A kind of minimizing technology of near surface scattering surface ripple |
| CN106950600B (en) * | 2017-02-16 | 2019-02-19 | 中国石油大学(华东) | A kind of minimizing technology of near surface scattering surface wave |
| CN109239780A (en) * | 2017-07-10 | 2019-01-18 | 中国石油化工股份有限公司 | Based on the synchronous method for squeezing wavelet transformation removal surface wave |
| CN113687420A (en) * | 2020-05-18 | 2021-11-23 | 中国石油天然气股份有限公司 | Method and device for removing surface wave noise of seismic data |
| CN113687420B (en) * | 2020-05-18 | 2023-08-22 | 中国石油天然气股份有限公司 | Method and device for removing surface wave noise of seismic data |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Mousavi et al. | Automatic noise-removal/signal-removal based on general cross-validation thresholding in synchrosqueezed domain and its application on earthquake data | |
| CN102193108B (en) | Method for improving SNR (signal-to-noise ratio) of oil exploration data processing | |
| US11880011B2 (en) | Surface wave prediction and removal from seismic data | |
| CN104360393A (en) | Seismic data reconstruction method | |
| CN103399348A (en) | Denoising method for seismic signal based on Shearlet transform | |
| CN105700020A (en) | Random noise suppression method and apparatus for seismic data | |
| CN104614769A (en) | Beam-forming filtering method for suppressing seismic surface waves | |
| CN110208856B (en) | Desert complex noise suppression method based on manifold partition 2D-VMD | |
| CN110261910A (en) | Seismic data surface wave minimizing technology based on adaptive sparse S-transformation | |
| CN104375185A (en) | Method and device for removing surface waves from seismic records | |
| CN104849757A (en) | System and method for eliminating random noise in seismic signals | |
| CN104635264B (en) | The processing method of earthquake data before superposition and equipment | |
| CN102841380B (en) | A kind of surface relief complex structure Earthquakes grouping of data is concerned with noise attenuation method | |
| CN110244360B (en) | Seismic data separation method and system based on effective frequency wave number domain de-aliasing | |
| CN105093282A (en) | Energy replacement surface wave suppressing method based on frequency constraint | |
| CN102323618A (en) | Coherent noise suppression method based on fractional order Fourier transformation | |
| CN103645504A (en) | Weak earthquake signal processing method based on generalized instantaneous phase and P norm negative norm | |
| CN105277987A (en) | Vibroseis harmonic wave suppressing method based on predictive filtering method and pure phase shifting method | |
| CN109959964B (en) | Broadband background noise suppression method for seismic signals of high-speed rail seismic source | |
| US20110075514A1 (en) | Apparatus and methods for attenuating seismic noise associated with atmospheric pressure fluctuations | |
| CN102338884B (en) | Elliptic window direction band-pass amplitude-preserved filtering data processing method in geophysical prospecting | |
| CN106291680A (en) | A kind of data low frequency continuation method | |
| CN101907726B (en) | Method for automatically identifying and eliminating industrial electrical interference in earthquake exploration | |
| CN112764108B (en) | Novel seismic data noise suppression algorithm based on improved empirical wavelet transformation | |
| CN112213785B (en) | Seismic data desert noise suppression method based on feature-enhanced denoising network |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150225 |