CN106597534B - A kind of method and device detaching different order Free Surface related multiples - Google Patents
A kind of method and device detaching different order Free Surface related multiples Download PDFInfo
- Publication number
- CN106597534B CN106597534B CN201611076127.7A CN201611076127A CN106597534B CN 106597534 B CN106597534 B CN 106597534B CN 201611076127 A CN201611076127 A CN 201611076127A CN 106597534 B CN106597534 B CN 106597534B
- Authority
- CN
- China
- Prior art keywords
- free surface
- composite wave
- wave
- rank
- indicate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 239000002131 composite material Substances 0.000 claims abstract description 139
- 238000000926 separation method Methods 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 7
- 238000003860 storage Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 240000006829 Ficus sundaica Species 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The embodiment of the present application provides a kind of method and device detaching different order Free Surface related multiples, wherein the method includes：Obtain seismic data；According to the seismic data, primary reflection and the first composite wave are obtained；According to the primary reflection and first composite wave, the second composite wave is obtained；According to second composite wave and first composite wave, 1 rank Free Surface related multiple is obtained；Since i=1, according to i+1 composite wave and i rank Free Surface related multiples, the i-th+2 composite wave is obtained, described i-th+2 composite wave be from i+2 ranks to N rank Free Surface related multiples, the i+1 composite wave be from i+1 ranks to N rank Free Surface related multiples, wherein, i=1,2 ... N；According to the described i-th+2 composite wave and the i+1 composite wave, i+1 rank Free Surface related multiples are obtained.The embodiment of the present application can separate the wave of different orders in Free Surface related multiple.
Description
Technical field
This application involves offshore earthquake data process fields, related more particularly, to a kind of different order Free Surfaces of separation
The method and device of multiple wave.
Background technology
Ocean multiple wave is that interference the most serious, ocean multiple wave can be divided by its propagation path in offshore earthquake data
Involve interbed multiple for ghosting multiple wave, seabed first multiple, with the relevant micro- Qu Duoci of water layer.Wherein, seabed with
The Free Surface related multiple generated between sea is a kind of most commonly seen multiple wave.With the rapid development of computer technology
With the increase of seafari difficulty, multiple wave carries the utilization of information as the research emphasis of scholars, how using more than once
Wave carries information and has become the new breakthrough point for improving offshore earthquake data image quality.It is exactly to believe from earthquake using multiple wave premise
Multiple wave is effectively extracted in number.
The sides SRME (Surface Related-Multiple Elimination) that Verschuur et al. (1992) is proposed
Method is the Free Surface related multiple damped system that current industrial quarters is generally approved and applied, and this method can be disposable
The Free Surface related multiple of all orders is predicted, but when this method extraction multiple wave, the multiple wave of different orders
It is mingled in one piece, cannot be distinguished.If the multiple wave of different orders can be isolated, both can targetedly be applied,
It can individually predict and interfere illusion caused by different order multiple waves, to push application and the hair of multiple wave imaging technique
Exhibition.Regrettably, there is presently no the separation methods for exclusively carrying out different order Free Surface related multiples.
Invention content
The method for being designed to provide a kind of different order Free Surface related multiples of separation of the embodiment of the present application and
Device separates the wave of different orders in Free Surface related multiple.
In order to achieve the above objectives, the embodiment of the present application provides a kind of different order Free Surface related multiples of separation
Method, the method includes：
Obtain seismic data；
According to the seismic data, primary reflection and the first composite wave are obtained, first composite wave is from 1 rank to N
Rank Free Surface related multiple；
According to the primary reflection and first composite wave, the second composite wave is obtained, second composite wave is from 2
Rank is to N rank Free Surface related multiples；
According to second composite wave and first composite wave, 1 rank Free Surface related multiple is obtained；
According to i+1 composite wave and i rank Free Surface related multiples, the i-th+2 composite wave, the described i-th+2 combination are obtained
Wave is from i+2 ranks to N rank Free Surface related multiples, and the i+1 composite wave is related from i+1 ranks to N ranks Free Surface
Multiple wave, wherein i=1,2 ... N；
According to the described i-th+2 composite wave and the i+1 composite wave, i+1 rank Free Surface related multiples are obtained.
The embodiment of the present application also provides a kind of device of the different order Free Surface related multiples of separation, described devices
Including：
Acquisition module, for obtaining seismic data；
Reflection wave groups multiplex obtains module, for according to the seismic data, obtaining primary reflection and the first composite wave,
First composite wave is from 1 rank to N rank Free Surface related multiples；
Composite wave obtains module one, for according to the primary reflection and first composite wave, obtaining the second combination
Wave, second composite wave are from 2 ranks to N rank Free Surface related multiples；
1 rank wave obtains module, for according to second composite wave and first composite wave, obtaining 1 rank Free Surface
Related multiple；
Composite wave obtains module two, for according to i+1 composite wave and i rank Free Surface related multiples, obtaining i-th+2
Composite wave, the described i-th+2 composite wave are from i+2 ranks to N rank Free Surface related multiples, and the i+1 composite wave is from i+1
Rank is to N rank Free Surface related multiples, wherein i=1,2 ... N；
Single-order wave obtains module, for according to the described i-th+2 composite wave and the i+1 composite wave, it is free to obtain i+1 ranks
Surface-related multiple.
The technical solution provided by above the embodiment of the present application is it is found that the embodiment of the present application is obtained according to seismic data first
Primary reflection is arrived and from 1 rank to N rank Free Surface related multiples, then according to primary reflection and from 1 rank to N
Rank Free Surface related multiple has been obtained from 2 ranks to N rank Free Surface related multiples, further according to from 2 ranks to N rank Free Surfaces
Face related multiple and 1 rank Free Surface related multiple is obtained from 1 rank to N rank Free Surface related multiples, this just succeeds
1 rank Free Surface related multiple is isolated；It is separately separated out different order Free Surfaces further according to above method successively recursion
Face related multiple.
Description of the drawings
Attached drawing described herein is used for providing further understanding the embodiment of the present application, constitutes the embodiment of the present application
A part does not constitute the restriction to the embodiment of the present application.In the accompanying drawings：
Fig. 1 is a kind of method flow signal of the different order Free Surface related multiples of separation of the embodiment of the present application
Figure；
Fig. 2 a are the slave i ranks of the embodiment of the present application to 4 rank Free Surface related multiple schematic diagrames；
Fig. 2 b are the slave i+1 ranks of the embodiment of the present application to 4 rank Free Surface related multiple schematic diagrames；
Fig. 2 c are the i+1 rank Free Surface related multiple schematic diagrames of the embodiment of the present application；
Fig. 3 is the test model velocity field of the embodiment of the present application；
Fig. 4 a are the forward modeling single shot record of the embodiment of the present application；
Fig. 4 b are the embodiment of the present application from 1 rank to 4 rank order Free Surface related multiple schematic diagrames；
Fig. 4 c are the primary reflection schematic diagram of the embodiment of the present application；
Fig. 5 a are the embodiment of the present application from 2 ranks to 4 rank Free Surface related multiple schematic diagrames；
Fig. 5 b are that the 1 rank Free Surface of the embodiment of the present application first closes multiple wave；
Fig. 6 a are the embodiment of the present application from 3 ranks to 4 rank Free Surface related multiple schematic diagrames；
Fig. 6 b are that the 2 rank Free Surfaces of the embodiment of the present application first close multiple wave；
Fig. 7 a are that the 4 rank Free Surfaces of the embodiment of the present application first close multiple wave；
Fig. 7 b are that the 3 rank Free Surfaces of the embodiment of the present application first close multiple wave；
Fig. 8 is a kind of schematic device of the different order Free Surface related multiples of separation of the embodiment of the present application.
Specific implementation mode
To make the purpose, technical scheme and advantage of the embodiment of the present application be more clearly understood, with reference to embodiment and attached
Figure, is described in further details the embodiment of the present application.Here, the illustrative embodiments and their description of the embodiment of the present application are used for
It explains the embodiment of the present application, but is not intended as the restriction to the embodiment of the present application.
Below in conjunction with the accompanying drawings, the specific implementation mode of the embodiment of the present application is described in further detail.
As shown in Figure 1, a kind of side for the different order Free Surface related multiples of separation that the embodiment of the present application is provided
Method may comprise steps of.
S101 obtains seismic data.
The seismic data can be that filtered pretreated single shot record or big gun collection record, and described that treated is single
It can be comprising primary reflection and from 1 rank to N rank Free Surface related multiples in big gun record or big gun collection record.Wherein, described
Big gun collection record can be the set of multiple single shot records.The acquisition earthquake record can be that reception is filtered pretreated
Single shot record or big gun collection record.
S102 obtains primary reflection and the first composite wave according to the seismic data, and first composite wave is from 1
Rank is to N rank Free Surface related multiples.
In a kind of embodiment of the application, SRME methods may be used, the seismic data is handled, obtain
Primary reflection and the first composite wave.
S103 obtains the second composite wave, second composite wave according to the primary reflection and first composite wave
For the Free Surface related multiple from 2 ranks to N ranks.
In a kind of embodiment of the application, can by frequency domain primary reflection and first composite wave substitute into
The second composite wave is calculated in following formula in the way of matrix multiple,
In formula, M0Indicate primary reflection in frequency domain.M1,NIndicate the first composite wave in frequency domain.M2,NIt indicates second in frequency domain
Composite wave.xsIndicate focus in Free Surface position.xrIndicate geophone station in Free Surface position.xkIndicate multiple wave in freedom
The reflection position on surface, x when reality calculateskChoose xsWith xrBetween integer value.ω indicates frequency domain, i.e., is accordingly calculated as frequency domain
The calculating of interior corresponding wave.
S104 obtains 1 rank Free Surface related multiple according to second composite wave and first composite wave.
In a kind of embodiment of the application, can by second composite wave and first composite wave, substitute into
1 rank Free Surface related multiple is calculated in lower formula,
M1(xr,xs；ω)=M1,N(xr,xs；ω)-α·M2,N(xr,xs；ω)
In formula, M1Indicate that 1 rank Free Surface related multiple of frequency domain, α indicate that matching attribute, ω indicate that corresponding wave is frequency
Result in domain.Wherein, the matching attribute is to ensure the quantity between second composite wave and first composite wave
Grade is consistent.In one embodiment of the application, matching attribute can be obtained according to least square Wiener filtering.
S105 obtains the i-th+2 composite wave according to i+1 composite wave and i rank Free Surface related multiples.
Described i-th+2 composite wave is from i+2 ranks to N rank Free Surface related multiples, and the i+1 composite wave is from i+
1 rank is to N rank Free Surface related multiples, wherein i=1,2 ... N.
In a kind of embodiment of the application, can by by frequency domain i+1 composite wave and the i ranks Free Surface
Related multiple substitutes into following formula, and the i-th+2 composite wave is calculated in the way of matrix multiple,
In formula, Mi+2,NIndicate the i-th+2 composite wave in frequency domain, Mi+1,NIndicate i+1 composite wave, MiIndicate i rank Free Surfaces
Face related multiple.
It is related repeatedly to obtain i+1 ranks Free Surface according to the described i-th+2 composite wave and the i+1 composite wave by S106
Wave.
In a kind of embodiment of the application, the described i-th+2 composite wave and the i+1 composite wave can be substituted into
I+1 rank Free Surface related multiples are calculated in following formula,
Mi+1(xr,xs；ω)=Mi+1,N(xr,xs；ω)-α·Mi+2,N(xr,xs；ω)
In the application one embodiment, matching attribute can be obtained according to least square Wiener filtering, matching attribute
Effect is to ensure that the i-th+2 composite wave and the i+1 composite wave are the same order of magnitude.
Fig. 2 a~Fig. 2 c are the principle schematic of embodiment illustrated in fig. 1, M in figureiIndicate that i ranks Free Surface is related repeatedly
Wave.Fig. 2 a indicate that from i ranks to N rank Free Surface related multiples, Fig. 2 b indicate related repeatedly from i+1 ranks to N ranks Free Surface
Fig. 2 a are matched under the same order of magnitude to Fig. 2 b and are wanted to subtract that can be obtained by i rank Free Surfaces as shown in Figure 2 c related by wave
Multiple wave.
The technical solution that flow chart is provided as shown in Figure 1 is it is found that the embodiment of the present application is obtained according to seismic data first
Primary reflection and from 1 rank to N rank Free Surface related multiples, then according to primary reflection and from 1 rank to N ranks
Free Surface related multiple has been obtained from 2 ranks to N rank Free Surface related multiples, further according to from 2 ranks to N rank Free Surfaces
Related multiple and 1 rank Free Surface related multiple is obtained from 1 rank to N rank Free Surface related multiples, this just successfully divides
1 rank Free Surface related multiple is separated out；It is related it to be separately separated out different order Free Surfaces further according to above method recursion
Multiple wave.
In one embodiment of the application, acquired seismic data be by pretreated single shot record, and
The single shot record is obtained according to test model velocity field inverting as shown in Figure 3.From the figure 3, it may be seen that the test model velocity field
It is divided into two layers, lateral maximum distance is 10km, and the speed of depth capacity 2.4km, first layer (depth is 0~0.6km) are
The speed of 2000m/s, the second layer (depth is 0.6~2.4km) are 3000m/s.Fig. 4 a be in the test model velocity field just
The single shot record drilled, which includes primary reflection P, single order Free Surface related multiple M1, second order Free Surface
Related multiple M2, three rank Free Surface related multiple M3, quadravalence Free Surface related multiple M4.It will as shown in fig. 4 a
Different order Free Surface related multiple separation can follow the steps below in single shot record.
(1) primary reflection P (as illustrated in fig. 4 c) in the single shot record is calculated based on SRME and from 1 rank to 4
The Free Surface related multiple M of rank1,4(as shown in Figure 4 b).
(2) by the primary reflection P and M in frequency domain1,4Following formula is substituted into, is calculated in the way of matrix multiple
From the Free Surface related multiple M of 2 ranks to 4 ranks2,4(as shown in Figure 5 a).
In formula, xsIndicate focus in Free Surface position, xrIndicate geophone station in Free Surface position, xkIndicate multiple wave
In the reflection position of Free Surface.
(3) by M1,4And M2,4Following formula is substituted into, single order Free Surface related multiple M is calculated1, such as Fig. 5 b institutes
Show.
M1(xr,xs；ω)=M1,4(xr,xs；ω)-α·M2,4(xr,xs；ω)
In formula, matching attribute α can be obtained by least square Wiener filtering.
(4) by M2,4And M1Following formula is substituted into, the Free Surface related multiple M from 3 ranks to 4 ranks is calculated3,4, such as
Shown in Fig. 6 a,
(5) by M2,4And M3,4Following formula is substituted into, second order Free Surface related multiple M is calculated2, such as Fig. 6 b institutes
Show.
M2(xr,xs；ω)=M1,4(xr,xs；ω)-α·M2,4(xr,xs；ω)
In formula, matching attribute α can be obtained by least square Wiener filtering.
(6) recursion in the way of above-mentioned steps (4) and (5), can isolate three rank Free Surface related multiple M3
(as shown in Figure 7b) and quadravalence Free Surface related multiple M4(as shown in Figure 7a).
In above-described embodiment, primary reflection is obtained according to seismic data first and from 1 rank to 4 rank Free Surface phases
Multiple wave is closed, has then been obtained according to primary reflection and from 1 rank to 4 rank Free Surface related multiples from 2 ranks to 4 ranks certainly
By surface-related multiple, further according to related more to from 1 rank to 4 rank Free Surfaces from 2 ranks to 4 rank Free Surface related multiples
Subwave obtains 1 rank Free Surface related multiple, this has just been successfully separated out 1 rank Free Surface related multiple；Further according to
Upper method recursion is separately separated out different order Free Surface related multiples.
A kind of device detaching different order Free Surface related multiples is additionally provided in the embodiment of the present application, it is such as following
Embodiment described in.The principle solved the problems, such as due to the device and a kind of side detaching different order Free Surface related multiples
Method is similar, therefore the implementation of the device may refer to a kind of method implementation detaching different order Free Surface related multiples,
Overlaps will not be repeated.
As shown in figure 8, a kind of device detaching different order Free Surface related multiples may include with lower module.
Acquisition module 801, for obtaining seismic data.
Reflection wave groups multiplex obtains module 802, for according to the seismic data, obtaining primary reflection and the first combination
Wave, first composite wave are from 1 rank to N rank Free Surface related multiples.
Composite wave obtains module 1, for according to the primary reflection and first composite wave, obtaining second group
Multiplex, second composite wave are the Free Surface related multiple from 2 ranks to N ranks.
1 rank wave obtains module 804, for according to second composite wave and first composite wave, obtaining 1 rank Free Surface
Face related multiple.
Composite wave obtains module 2 805, for according to i+1 composite wave and i rank Free Surface related multiples, obtaining the
I+2 composite waves, the described i-th+2 composite wave be from i+2 ranks to N rank Free Surface related multiples, the i+1 composite wave be from
I+1 ranks are to N rank Free Surface related multiples, wherein i=1,2 ... N.
Single-order wave obtains module 806, for according to the described i-th+2 composite wave and the i+1 composite wave, obtaining i+1 ranks
Free Surface related multiple.
The technical solution provided by above-mentioned the application device embodiment is it is found that the embodiment of the present application provided according to earthquake first
Material has obtained primary reflection and from 1 rank to N rank Free Surface related multiples, then according to primary reflection and from 1
Rank to N rank Free Surface related multiples have been obtained from 2 ranks to N rank Free Surface related multiples, further according to from 2 ranks to N ranks
Free Surface related multiple and 1 rank Free Surface related multiple is obtained from 1 rank to N rank Free Surface related multiples, this
Just it has been successfully separated out 1 rank Free Surface related multiple；It is free it to be separately separated out different orders further according to above method recursion
Surface-related multiple.
The step of method or algorithm described in the embodiment of the present application can be directly embedded into hardware, processor execute it is soft
The combination of part module or the two.Software module can be stored in RAM memory, flash memory, ROM memory, EPROM storages
Other any form of storaging mediums in device, eeprom memory, register, hard disk, moveable magnetic disc, CD-ROM or this field
In.Illustratively, storaging medium can be connect with processor, so that processor can read information from storaging medium, and
It can be to storaging medium stored and written information.Optionally, storaging medium can also be integrated into processor.Processor and storaging medium can
To be set in ASIC, ASIC can be set in user terminal.Optionally, processor and storaging medium can also be set to use
In different components in the terminal of family.
In one or more illustrative designs, above-mentioned function described in the embodiment of the present application can be in hardware, soft
Part, firmware or the arbitrary of this three combine to realize.If realized in software, these functions can store and computer-readable
On medium, or with one or more instruction or code form be transmitted on the medium of computer-readable.Computer readable medium includes electricity
Brain storaging medium and convenient for allow computer program to be transferred to from a place telecommunication media in other places.Storaging medium can be with
It is that any general or special computer can be with the useable medium of access.For example, such computer readable media may include but
It is not limited to RAM, ROM, EEPROM, CD-ROM or other optical disc storage, disk storage or other magnetic storage devices or other
What can be used for carry or store with instruct or data structure and it is other can be by general or special computer or general or specially treated
The medium of the program code of device reading form.In addition, any connection can be properly termed computer readable medium, example
Such as, if software is to pass through a coaxial cable, fiber optic cables, double from a web-site, server or other remote resources
Twisted wire, Digital Subscriber Line (DSL) are defined with being also contained in for the wireless way for transmitting such as example infrared, wireless and microwave
In computer readable medium.The disk (disk) and disk (disc) includes compress disk, radium-shine disk, CD, DVD, floppy disk
And Blu-ray Disc, disk is usually with magnetic duplication data, and disk usually carries out optical reproduction data with laser.Combinations of the above
It can also be included in computer readable medium.
Particular embodiments described above has carried out further in detail the purpose, technical solution and advantageous effect of the application
It describes in detail bright, it should be understood that the foregoing is merely the specific embodiment of the embodiment of the present application, is not used to limit this Shen
Protection domain please, within the spirit and principles of this application, any modification, equivalent substitution, improvement and etc. done should all wrap
It is contained within the protection domain of the application.
Claims (10)
1. a kind of method detaching different order Free Surface related multiples, which is characterized in that include the following steps：
Obtain seismic data；
According to the seismic data, primary reflection and the first composite wave are obtained, first composite wave is for oneself from 1 rank to N ranks
By surface-related multiple；
According to the primary reflection and first composite wave, obtain the second composite wave, second composite wave be from 2 ranks to
N rank Free Surface related multiples；
According to second composite wave and first composite wave, 1 rank Free Surface related multiple is obtained；
Since i=1, according to i+1 composite wave and i rank Free Surface related multiples, the i-th+2 composite wave is obtained, described i-th
+ 2 composite waves are from i+2 ranks to N rank Free Surface related multiples, and the i+1 composite wave is from i+1 ranks to N rank Free Surfaces
Face related multiple, wherein i=1,2 ... N；
According to the described i-th+2 composite wave and the i+1 composite wave, i+1 rank Free Surface related multiples are obtained；
Wherein, described that i-th+2 composite wave is obtained according to i+1 composite wave and i rank Free Surface related multiples, specifically,
I+1 composite wave and i rank Free Surface related multiples are calculated i-th+2 in frequency domain according to following formula
Composite wave,
In formula, Mi+2,NIndicate the i-th+2 composite wave, Mi+1,NIndicate i+1 composite wave, MiIndicate i rank Free Surface related multiples,
ω indicates that corresponding wave is that frequency domain is interior as a result, xsIndicate focus in Free Surface position, xrGeophone station is indicated in Free Surface position,
xkReflection position of the expression multiple wave in Free Surface；
It is described according to the described i-th+2 composite wave and the i+1 composite wave, obtain i+1 rank Free Surface related multiples, specifically
For,
By the described i-th+2 composite wave and the i+1 composite wave, substitutes into following formula and i+1 ranks Free Surface correlation is calculated
Multiple wave,
Mi+1(xr,xs；ω)=Mi+1,N(xr,xs；ω)-α·Mi+2,N(xr,xs；ω)
In formula, Mi+1Indicate that i+1 rank Free Surface related multiples in frequency domain, α indicate matching attribute.
2. the method as described in claim 1, which is characterized in that it is described according to the seismic data, obtain primary reflection and
First composite wave, specifically,
According to the seismic data, primary reflection and the first composite wave are obtained using SRME methods.
3. method as claimed in claim 2, which is characterized in that described to be combined with described first according to the primary reflection
Wave obtains the second composite wave, specifically,
The second composite wave is calculated according to following formula in frequency domain in primary reflection and first composite wave,
In formula, M0Indicate primary reflection in frequency domain, M1,NIndicate the first composite wave in frequency domain, M2,NIndicate the second combination in frequency domain
Wave, xsIndicate focus in Free Surface position, xrIndicate geophone station in Free Surface position, xkIndicate multiple wave in Free Surface
Reflection position, ω indicate frequency domain.
4. method as claimed in claim 3, which is characterized in that described to be combined with described first according to second composite wave
Wave obtains 1 rank Free Surface related multiple, specifically,
By second composite wave and first composite wave, substitutes into following formula and 1 rank Free Surface correlation is calculated repeatedly
Wave,
M1(xr,xs；ω)=M1,N(xr,xs；ω)-α·M2,N(xr,xs；ω)
In formula, M1Indicate that 1 rank Free Surface related multiple of frequency domain, α indicate that matching attribute, ω indicate that corresponding wave is in frequency domain
As a result.
5. the method as described in claim 1, which is characterized in that the matching attribute is obtained according to least square Wiener filtering.
6. a kind of device of the different order Free Surface related multiples of separation, which is characterized in that described device includes：
Acquisition module, for obtaining seismic data；
Reflection wave groups multiplex obtains module, described for according to the seismic data, obtaining primary reflection and the first composite wave
First composite wave is from 1 rank to N rank Free Surface related multiples；
Composite wave obtains module one, for according to the primary reflection and first composite wave, obtaining the second composite wave, institute
It is from 2 ranks to N rank Free Surface related multiples to state the second composite wave；
1 rank wave obtains module, for according to second composite wave and first composite wave, it is related to obtain 1 rank Free Surface
Multiple wave；
Composite wave obtains module two, for according to i+1 composite wave and i rank Free Surface related multiples, obtaining the i-th+2 combination
Wave, the described i-th+2 composite wave be from i+2 ranks to N rank Free Surface related multiples, the i+1 composite wave be from i+1 ranks to
N rank Free Surface related multiples, wherein i=1,2 ... N；
Single-order wave obtains module, for according to the described i-th+2 composite wave and the i+1 composite wave, obtaining i+1 rank Free Surfaces
Related multiple；Wherein,
The composite wave obtains module two and is specifically used for, by i+1 composite wave and i rank Free Surface related multiples in frequency domain
It is interior that the i-th+2 composite wave is calculated according to following formula,
In formula, Mi+2,NIndicate the i-th+2 composite wave, Mi+1,NIndicate i+1 composite wave, MiIndicate i rank Free Surface related multiples,
ω indicates that corresponding wave is that frequency domain is interior as a result, xsIndicate focus in Free Surface position, xrGeophone station is indicated in Free Surface position,
xkReflection position of the expression multiple wave in Free Surface；
The single-order wave obtains module and is specifically used for, and the described i-th+2 composite wave and the i+1 composite wave substitute into following public
I+1 rank Free Surface related multiples are calculated in formula,
Mi+1(xr,xs；ω)=Mi+1,N(xr,xs；ω)-α·Mi+2,N(xr,xs；ω)
In formula, Mi+1Indicate that i+1 rank Free Surface related multiples in frequency domain, α indicate matching attribute.
7. device as claimed in claim 6, which is characterized in that the reflection wave groups multiplex obtains module and is specifically used for according to institute
Seismic data is stated, primary reflection and the first composite wave are obtained using SRME methods.
8. device as claimed in claim 7, which is characterized in that the composite wave obtains module one and is specifically used for, will be once anti-
In frequency domain, according to following formula, the second composite wave is calculated in ejected wave and first composite wave,
In formula, M0Indicate primary reflection in frequency domain, M1,NIndicate the first composite wave in frequency domain, M2,NIndicate the second combination in frequency domain
Wave, xsIndicate focus in Free Surface position, xrIndicate geophone station in Free Surface position, xkIndicate multiple wave in Free Surface
Reflection position, ω indicate frequency domain.
9. device as claimed in claim 8, which is characterized in that the 1 rank wave obtains module and is specifically used for, by described second group
Multiplex and first composite wave substitute into following formula and 1 rank Free Surface related multiple are calculated,
M1(xr,xs；ω)=M1,N(xr,xs；ω)-α·M2,N(xr,xs；ω)
In formula, M1Indicate that 1 rank Free Surface related multiple of frequency domain, α indicate that matching attribute, ω indicate that corresponding wave is in frequency domain
As a result.
10. device as claimed in claim 6, which is characterized in that the matching attribute is obtained according to least square Wiener filtering.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611076127.7A CN106597534B (en) | 2016-11-30 | 2016-11-30 | A kind of method and device detaching different order Free Surface related multiples |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611076127.7A CN106597534B (en) | 2016-11-30 | 2016-11-30 | A kind of method and device detaching different order Free Surface related multiples |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106597534A CN106597534A (en) | 2017-04-26 |
| CN106597534B true CN106597534B (en) | 2018-10-16 |
Family
ID=58593847
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611076127.7A Active CN106597534B (en) | 2016-11-30 | 2016-11-30 | A kind of method and device detaching different order Free Surface related multiples |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106597534B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108957546A (en) * | 2018-06-14 | 2018-12-07 | 中国石油天然气股份有限公司 | Prediction technique and system is imaged in multiple wave |
| CN109358363B (en) * | 2018-11-16 | 2019-05-24 | 中国海洋大学 | Three-dimensional singing tracking and prediction technique based on multiple reflections equivalent inter-phase model |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105334537A (en) * | 2015-10-26 | 2016-02-17 | 中国石油大学(华东) | Primary wave and multiple wave separation method based on alternative splitting Bregman iterative algorithm |
| CN105510976A (en) * | 2016-01-11 | 2016-04-20 | 中国海洋大学 | Multiple-wave combined self-adaptive attenuation method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8077544B2 (en) * | 2008-03-28 | 2011-12-13 | Westerngeco L.L.C. | Dual-wavefield multiple attenuation |
-
2016
- 2016-11-30 CN CN201611076127.7A patent/CN106597534B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105334537A (en) * | 2015-10-26 | 2016-02-17 | 中国石油大学(华东) | Primary wave and multiple wave separation method based on alternative splitting Bregman iterative algorithm |
| CN105510976A (en) * | 2016-01-11 | 2016-04-20 | 中国海洋大学 | Multiple-wave combined self-adaptive attenuation method |
Non-Patent Citations (2)
| Title |
|---|
| 分离多次波的独立分量分析法;耿建华 等;《同济大学学报（自然科学版）》;20130331;第41卷(第3期);第449-457页 * |
| 基于单程波偏移算子的地表相关多次波成像;叶月明 等;《地球物理学报》;20140430;第57卷(第4期);第1241-1250页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106597534A (en) | 2017-04-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Maggi et al. | An automated time-window selection algorithm for seismic tomography | |
| Wapenaar et al. | Marchenko imaging | |
| CN104536044B (en) | The interpolation denoising method and system of a kind of geological data | |
| Zhang et al. | Free-surface and internal multiple elimination in one step without adaptive subtraction | |
| Wang et al. | Reverse time migration of multiples: Eliminating migration artifacts in angle domain common image gathers | |
| WO2012021218A2 (en) | Attenuating internal multiples from seismic data | |
| Rocha et al. | Isotropic elastic wavefield imaging using the energy norm | |
| CN106597534B (en) | A kind of method and device detaching different order Free Surface related multiples | |
| CN103308941B (en) | A kind of formation method based on any wide angle wave equation and device | |
| CN107193045A (en) | A kind of seismic data processing technique and device | |
| Guo et al. | Target-oriented waveform redatuming and high-resolution inversion: Role of the overburden | |
| WO2015082419A1 (en) | Multi-dimensional deconvolution using exact boundary conditions | |
| CN103018777B (en) | A kind of formation method and device that can eliminate reverse-time migration noise | |
| BR102013000835A2 (en) | VARIABLE DEPTH STREAMER SRME | |
| CN104932011B (en) | A kind of method and device for picking up seismic first breaks | |
| CN104580183B (en) | A kind of method of data synchronization and device across cluster | |
| CN104570114B (en) | A kind of reverse-time migration Noise Elimination method based on wavefield decomposition | |
| CN105093328A (en) | Slip-sweep harmonic suppression method and device | |
| AU2015200555B2 (en) | Correction of sea surface state | |
| CA2985738A1 (en) | Method for generating multiple free seismic images | |
| CN104199105A (en) | Pickup method and device of first arrivals in high-density three-dimensional exploration | |
| CN105403919B (en) | A kind of reverse-time migration imaging method and device | |
| CN109358363B (en) | Three-dimensional singing tracking and prediction technique based on multiple reflections equivalent inter-phase model | |
| CN103954994B (en) | Seismic signal Enhancement Method and device based on continuous wavelet transform | |
| Zhang et al. | Anisotropic viscoacoustic wave modelling in VTI media using frequency-dependent complex velocity |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |