CN104834010B - The poststack Forecasting Methodology and device of many subwaves between two interfaces - Google Patents
The poststack Forecasting Methodology and device of many subwaves between two interfaces Download PDFInfo
- Publication number
- CN104834010B CN104834010B CN201510240408.0A CN201510240408A CN104834010B CN 104834010 B CN104834010 B CN 104834010B CN 201510240408 A CN201510240408 A CN 201510240408A CN 104834010 B CN104834010 B CN 104834010B
- Authority
- CN
- China
- Prior art keywords
- mrow
- interface
- msub
- contact surface
- section
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a kind of poststack Forecasting Methodology of many subwaves between two interface and device, wherein, this method includes：Obtain it is deep during given first and when second it is deep, wherein, be more than deeply deep when first when second, deeply corresponding subsurface reflective boundary is the first interface when first, and deeply corresponding subsurface reflective boundary is second contact surface when second；Stacked section is divided into by Part I data and Part II data with the first interface；The energy for the seismic channel for obtaining many subwaves in the first interface is calculated according to the whole process and pegleg multiples at the first interface, the energy for the seismic channel for obtaining many subwaves of second contact surface is calculated according to the whole process and pegleg multiples of second contact surface, and calculating obtains normalized cross-correlation section；According to normalized related section and Part II data, the interbed multiple obtained between the first interface and second contact surface is calculated.The above embodiment of the present invention has reached the technique effect for effectively improving forecasting accuracy and forecasting efficiency.
Description
Technical field
The present invention relates to technical field of geophysical exploration, the poststack of many subwaves between more particularly to a kind of two interface is pre-
Survey method and apparatus.
Background technology
Interface between some special rock mass and country rock often produces stronger many subwaves or interbed multiple, example
Such as：The exploratory area of the developments such as carbonate rock, igneous rock and coal seam is just easy to produce stronger many subwaves or interbed multiple, this meeting
Large effect is compared in velocity analysis, imaging and inverting generation to seismic data.
In order to reduce this influence, generally require and many subwaves are suppressed, provided at present for the higher earthquake of signal to noise ratio
Material, multiple suppression typically uses prestack method, and for the low seismic data of signal to noise ratio, multiple suppression typically uses poststack
Method.Further, it is contemplated that method for marine seismic data signal to noise ratio is higher, land low signal-to-noise ratio seismic data signal to noise ratio is low, therefore
The many subwaves for suppressing method for marine seismic data typically use prestack method, and many subwaves of compacting land seismic data typically use poststack side
Method, poststack method is specifically described below：
The method that the method for poststack multiple suppression mainly uses wave field extrapolation, that is, give the one or more generations in underground many
The interface of subwave and corresponding stack velocity, then by a down going wave wave field extrapolation and a upgoing wave wave field extrapolation come pre-
Survey the first multiple and pegleg multiples relevant with the interface.It is theoretical for the interbed multiple prediction between two interfaces
On can be realized by six wave field extrapolations, however, the shortcoming of this implementation be need two interfaces between layer speed
Degree, and by stack velocity be converted into interval velocity process can so that there is larger error in the result of interval velocity, and, it is more multiple
Wave field extrapolation can cause computational efficiency relatively low.
The problem of for how efficiently and accurately to predict the interbed multiple between two interfaces, not yet propose at present effective
Solution.
The content of the invention
The embodiments of the invention provide a kind of poststack Forecasting Methodology of many subwaves between two interface, to solve prior art
In two interfaces interbed multiple prediction accuracy it is not high, the low technical problem of forecasting efficiency, this method includes：
Obtain stacked section and stacking velocity field corresponding with the stacked section；
Obtain it is deep during given first and when second it is deep, wherein, be more than deeply when second it is deep when first, it is deeply corresponding when first
Subsurface reflective boundary is the first interface, and deeply corresponding subsurface reflective boundary is second contact surface when second；
The stacked section is divided into by Part I data and Part II data with first interface, wherein, it is described
Part I data are the data of first interface above section, and the Part II data are first interface with bottom
The data divided；
Found out from the stacking velocity field when described first deep when deeply corresponding first stack velocity and described second
Corresponding second stack velocity；
The first depth for obtaining first interface is calculated according to deep and first stack velocity when described first, according to
Deep and second stack velocity calculates the second depth for obtaining the second contact surface when described second；
According to first depth and first stack velocity, using downgoing wave equation by the stacked section from benchmark
Face continuation obtains the first wave field to first interface, then using up-going wave equation, by first wave field from first boundary
Face continuation obtains the whole process and pegleg multiples at first interface to the reference plane；
According to second depth and second stack velocity, using downgoing wave equation by the stacked section from benchmark
Face continuation obtains the second wave field to the second contact surface, then using up-going wave equation, by second wave field from second boundary
Face continuation obtains the whole process and pegleg multiples of the second contact surface to the reference plane；
The seismic channel for obtaining many subwaves in the first interface is calculated according to the whole process and pegleg multiples at first interface
Energy, the seismic channel for obtaining the second contact surface many subwaves is calculated according to the whole process and pegleg multiples of the second contact surface
Energy；
According to the whole process and pegleg multiples and the whole process and pegleg multiples of the second contact surface at first interface, meter
Calculation obtains the cross-correlation section between first interface and many subwaves of the second contact surface；
According to the energy of the seismic channel of many subwaves in the first interface, the energy of the seismic channel of many subwaves of the second contact surface
With the cross-correlation section, calculating obtains normalized cross-correlation section；
According to the normalized related section and the Part II data, calculating obtains first interface and described
Interbed multiple between second contact surface.
In one embodiment, first depth and the second depth are calculated according to below equation：
Z1(x, y)=T1(x,y) V1(x,y)
Z2(x, y)=T2(x,y) V2(x,y)
Wherein, V1(x, y)=V (x, y, T1(x, y)), V2(x, y)=V (x, y, T2(x, y)), Z1Represent the first depth, Z2
Represent the second depth, T1Represent deep when first, T2Represent deep when second, V1Represent the first stack velocity, V2Represent the second superposition speed
Degree, V represents the stacking velocity field of the stacked section, and x represents the vertical line direction coordinate of stacked section, and y represents stacked section
Horizontal line direction coordinate.
In one embodiment, energy and the institute of the seismic channel of many subwaves in the first interface are calculated according to below equation
State the energy of the seismic channel of many subwaves of second contact surface：
Wherein, M1(x, y, t) represents the whole process and pegleg multiples at the first interface, M2(x, y, t) represents the complete of second contact surface
Journey and pegleg multiples, E1(x, y) represents the energy of the seismic channel of many subwaves in the first interface, E2(x, y) represents that second contact surface is multiple
The energy of the seismic channel of ripple, x represents the vertical line direction coordinate of stacked section, and y represents the horizontal line direction coordinate of stacked section,
T represents time coordinate, tmRepresent the dominant record time of stacked section.
In one embodiment, according to below equation calculate first interface and the second contact surface many subwaves it
Between cross-correlation section：
Wherein, R (x, y, t) represents the cross-correlation section between first interface and many subwaves of the second contact surface,
M1(x, y, t) represents the whole process and pegleg multiples at the first interface, M2(x, y, t) represents that the whole and micro- of second contact surface bends multiple
Ripple,Represent computing cross-correlation.
In one embodiment, normalized cross-correlation section is calculated according to below equation：
Wherein, Rn(x, y, t) represents normalized cross-correlation section.
In one embodiment, the interlayer between first interface and the second contact surface is calculated according to below equation
Many subwaves：
M3(x, y, t)=Rn(x,y,t)*D2(x,y,t)
Wherein, M3(x, y, t) represents the interbed multiple between first interface and the second contact surface, D2(x,y,t)
The Part II data are represented, * represents convolution operation.
The embodiment of the present invention additionally provides a kind of poststack prediction meanss of many subwaves between two interface, to solve existing skill
The accuracy of the prediction of the interbed multiple at two interfaces is not high in art, and the low technical problem of forecasting efficiency, the device includes：
Stacked section and stacking velocity field acquiring unit, for obtaining stacked section and corresponding with the stacked section folded
Acceleration field；
When deep acquiring unit, it is deep when given for obtaining first and when second it is deep, wherein, when being more than first deeply when second
Deep, deeply corresponding subsurface reflective boundary is the first interface when first, and deeply corresponding subsurface reflective boundary is second contact surface when second；
Stacked section division unit, for the stacked section to be divided into Part I data and with first interface
Two partial datas, wherein, the Part I data are the data of first interface above section, the Part II data
For the data of part below first interface；
Searching unit, for found out from the stacking velocity field when described first deeply corresponding first stack velocity and
Deeply corresponding second stack velocity when described second；
Depth calculation unit, first boundary is obtained for being calculated according to deep and first stack velocity when described first
First depth in face, it is deep according to obtain the second contact surface when described second with second stack velocity calculating deeply second
Degree；
First continuation unit, for according to first depth and first stack velocity, being incited somebody to action using downgoing wave equation
The stacked section obtains the first wave field from reference plane continuation to first interface, then using up-going wave equation, by described
One wave field obtains the whole process and pegleg multiples at first interface from first interface continuation to the reference plane；
Second continuation unit, for according to second depth and second stack velocity, being incited somebody to action using downgoing wave equation
The stacked section obtains the second wave field from reference plane continuation to the second contact surface, then using up-going wave equation, by described
Two wave fields obtain the whole process and pegleg multiples of the second contact surface from the second contact surface continuation to the reference plane；
The energy calculation unit of seismic channel, institute is obtained for being calculated according to the whole process and pegleg multiples at first interface
The energy of the seismic channel of many subwaves in the first interface is stated, is calculated according to the whole process and pegleg multiples of the second contact surface described in obtaining
The energy of the seismic channel of many subwaves of second contact surface；
Cross-correlation section computing unit, for according to the whole process at first interface and pegleg multiples and second boundary
The whole process and pegleg multiples in face, calculate the cross-correlation obtained between first interface and many subwaves of the second contact surface and cut open
Face；
Normalization unit, energy, the second contact surface for the seismic channel according to many subwaves in the first interface is multiple
The energy of the seismic channel of ripple and the cross-correlation section, calculating obtain normalized cross-correlation section；
Interbed multiple computing unit, for according to the normalized related section and the Part II data, meter
Calculation obtains the interbed multiple between first interface and the second contact surface.
In one embodiment, the depth calculation unit according to below equation specifically for calculating first depth
With the second depth：
Z1(x, y)=T1(x,y) V1(x,y)
Z2(x, y)=T2(x,y) V2(x,y)
Wherein, V1(x, y)=V (x, y, T1(x, y)), V2(x, y)=V (x, y, T2(x, y)), Z1Represent the first depth, Z2
Represent the second depth, T1Represent deep when first, T2Represent deep when second, V1Represent the first stack velocity, V2Represent the second superposition speed
Degree, V represents the stacking velocity field of the stacked section, and x represents the vertical line direction coordinate of stacked section, and y represents stacked section
Horizontal line direction coordinate.
In one embodiment, the energy calculation unit of the seismic channel according to below equation specifically for calculating described
The energy of the energy of the seismic channel of many subwaves in first interface and the seismic channel of many subwaves of the second contact surface：
Wherein, M1(x, y, t) represents the whole process and pegleg multiples at the first interface, M2(x, y, t) represents the complete of second contact surface
Journey and pegleg multiples, E1(x, y) represents the energy of the seismic channel of many subwaves in the first interface, E2(x, y) represents that second contact surface is multiple
The energy of the seismic channel of ripple, x represents the vertical line direction coordinate of stacked section, and y represents the horizontal line coordinates of stacked section, t tables
Show time coordinate, tmRepresent the dominant record time of stacked section.
In one embodiment, the cross-correlation section computing unit according to below equation specifically for calculating described the
Cross-correlation section between one interface and many subwaves of the second contact surface：
Wherein, R (x, y, t) represents the cross-correlation section between first interface and many subwaves of the second contact surface,
M1(x, y, t) represents the whole process and pegleg multiples at the first interface, M2(x, y, t) represents that the whole and micro- of second contact surface bends multiple
Ripple,Represent computing cross-correlation.
In embodiments of the present invention, it would be desirable to could be realized by two interfacial layer speed and six wave field extrapolations two
The prediction of the interbed multiple at interface, two interfaces for being reduced to only to need stack velocity and four wave field extrapolations just to realize
The prediction of interbed multiple, so that the accuracy for solving the prediction of the interbed multiple at two interfaces in the prior art is not high, in advance
The technical problem of inefficiency is surveyed, the technique effect for effectively improving forecasting accuracy and forecasting efficiency has been reached.
Brief description of the drawings
Accompanying drawing described herein is used for providing a further understanding of the present invention, constitutes the part of the application, not
Constitute limitation of the invention.In the accompanying drawings：
Fig. 1 is the method flow diagram of the poststack Forecasting Methodology of many subwaves between two interfaces according to embodiments of the present invention；
Fig. 2 is two ground lower bounds that the stacked section of Pluto generated datas according to embodiments of the present invention and user provide
Face schematic diagram；
Geological data when Fig. 3 is second after data according to embodiments of the present invention are split below depth；
Fig. 4 is the whole process and pegleg multiples ray path schematic diagram at a certain interface in underground according to embodiments of the present invention；
Fig. 5 is that the whole process and pegleg multiples at the first interface of use wave field extrapolation prediction according to embodiments of the present invention are shown
It is intended to；
Fig. 6 is that the whole process and pegleg multiples of the second contact surface of use wave field extrapolation prediction according to embodiments of the present invention are shown
It is intended to；
Fig. 7 is the normalizing between the first interface according to embodiments of the present invention and the whole process and pegleg multiples of second contact surface
Change cross-correlation diagrammatic cross-section；
Fig. 8 is the interbed multiple signal between the first interface of prediction according to embodiments of the present invention and second interface
Figure；
Fig. 9 is that the multiple wave self-adaption of use according to embodiments of the present invention subtracts each other, and suppresses the whole process at the first and first interface
Stacked section schematic diagram after interbed multiple between pegleg multiples and the first and second interfaces；
Figure 10 is the structured flowchart of the poststack prediction meanss of many subwaves between two interfaces according to embodiments of the present invention.
Embodiment
It is right with reference to embodiment and accompanying drawing for the object, technical solutions and advantages of the present invention are more clearly understood
The present invention is described in further details.Here, the exemplary embodiment of the present invention and its illustrating to be used to explain the present invention, but simultaneously
It is not as a limitation of the invention.
In view of in the prior art predict two interfaces many subwaves when, why exist resultant error than it is larger,
The problem of efficiency comparison is low, is primarily due to need to calculate interval velocity during prediction, and wave field extrapolation number of times is relatively more.For
This, inventors herein proposes a kind of poststack Forecasting Methodology of many subwaves between two interface, theoretical based on backscattering, it would be desirable to two
Interfacial layer speed and six wave field extrapolations are reduced to only need to using pre- come the interbed multiple method predicted between two interfaces
The whole process and pegleg multiples for two subsurface interfaces surveyed, it is possible to predict the mode of the interbed multiple between two interfaces, subtract
Lack computing cost, also improve computational accuracy.
As shown in figure 1, the poststack Forecasting Methodology of many subwaves between two interface comprises the following steps：
Step 101：Stacked section and stacking velocity field corresponding with the stacked section are obtained, is to obtain as shown in Figure 2
Pluto generated datas stacked section schematic diagram；
Step 102：Obtain it is deep during given first and when second it is deep, wherein, be more than deeply when second it is deep when first, when first
Deep corresponding subsurface reflective boundary is the first interface, and deeply corresponding subsurface reflective boundary is second contact surface when second；
Specifically, this when deeply can user give, two articles of white dashed lines as shown in Figure 2 are that user gives
A period of time depth T1With depth T when second2；
Step 103：The stacked section is divided into by Part I data and Part II data with first interface, its
In, the Part I data are the data of first interface above section, and the Part II data are first boundary
The data of part below face；It is illustrated in figure 3 the Part II data D for carrying out being obtained after data segmentation2(x, y, t) schematic diagram.
Step 104：Deeply corresponding first stack velocity is found out when described first from the stacking velocity field and described
Deeply corresponding second stack velocity when second；
Step 105：Calculated according to deep and first stack velocity when described first and obtain the first of first interface
Depth, the second depth for obtaining the second contact surface is calculated according to deep and second stack velocity when described second；
For example：When being got by above-mentioned steps after deep and stack velocity, first can be calculated according to below equation deep
Degree and the second depth：
Z1(x, y)=T1(x,y) V1(x,y)
Z2(x, y)=T2(x,y) V2(x,y)
Wherein, V1(x, y)=V (x, y, T1(x, y)), V2(x, y)=V (x, y, T2(x, y)), Z1Represent the first depth, Z2
Represent the second depth, T1Represent deep when first, T2Represent deep when second, V1Represent the first stack velocity, V2Represent the second superposition speed
Degree, V represents the stacking velocity field of the stacked section, and x represents the vertical line direction coordinate of stacked section, and y represents stacked section
Horizontal line direction coordinate.
Step 106：According to first depth and first stack velocity, the superposition is cutd open using downgoing wave equation
Face obtains the first wave field from reference plane continuation to first interface, then using up-going wave equation, by first wave field from institute
The first interface continuation is stated to the reference plane, the whole process and pegleg multiples at first interface is obtained；Similarly, according to described
Second depth and second stack velocity, using downgoing wave equation by the stacked section from reference plane continuation to described second
Interface, obtains the second wave field, then using up-going wave equation, by second wave field from the second contact surface continuation to the benchmark
Face, obtains the whole process and pegleg multiples of the second contact surface；
The continuation processing based on reference plane is carried out namely based on depth and stack velocity data, as shown in figure 4, being underground
Many subwaves produce the whole process and pegleg multiples ray path schematic diagram at interface, are illustrated in figure 5 in the way of in step 106
Predicted after continuation the whole process and pegleg multiples M at the first obtained interface1(x, y, t) schematic diagram, be illustrated in figure 6 according to
Mode in step 106 predicted after continuation the whole process and pegleg multiples M of obtained second contact surface2(x, y, t) schematic diagram.
Step 107：Calculated according to the whole process and pegleg multiples at first interface and obtain many subwaves in the first interface
Seismic channel energy, according to the second contact surface whole process and pegleg multiples calculate obtain many subwaves of the second contact surface
The energy of seismic channel；
When actually performing, the energy for obtaining seismic channel can be calculated according to below equation：
Wherein, M1(x, y, t) represents the whole process and pegleg multiples at the first interface, M2(x, y, t) represents the complete of second contact surface
Journey and pegleg multiples, E1(x, y) represents the energy of the seismic channel of many subwaves in the first interface, E2(x, y) represents that second contact surface is multiple
The energy of the seismic channel of ripple, x represents the vertical line direction coordinate of stacked section, and y represents the horizontal line direction coordinate of stacked section,
T represents time coordinate, tmRepresent the dominant record time of stacked section.
Step 108：Bent according to the whole process and pegleg multiples at first interface and the whole and micro- of the second contact surface
Many subwaves, calculate the cross-correlation section obtained between first interface and many subwaves of the second contact surface；
When implementing, cross-correlation section can be calculated according to below equation：
Wherein, R (x, y, t) represents the cross-correlation section between first interface and many subwaves of the second contact surface,
M1(x, y, t) represents the whole process and pegleg multiples at the first interface, M2(x, y, t) represents that the whole and micro- of second contact surface bends multiple
Ripple,Represent computing cross-correlation.
Step 109：According to the earthquake of the energy of the seismic channel of many subwaves in the first interface, many subwaves of the second contact surface
The energy in road and the cross-correlation section, calculating obtain normalized cross-correlation section；
Specifically, normalized cross-correlation section can be calculated according to below equation：
Wherein, Rn(x, y, t) represents normalized cross-correlation section, as shown in fig. 7, being the first interface and second contact surface
Cross-correlation section between many subwaves, and normalized cross-correlation diagrammatic cross-section.
Step 110：According to the normalized related section and the Part II data, calculating obtains first boundary
Interbed multiple between face and the second contact surface：
Specifically, based on normalized related section and Part II data are obtained in above-mentioned steps, can be according to following
Formula calculates the interbed multiple between first interface and the second contact surface：
M3(x, y, t)=Rn(x,y,t)*D2(x,y,t)
Wherein, M3(x, y, t) represents the interbed multiple between first interface and the second contact surface, D2(x,y,t)
The Part II data are represented, * represents convolution operation, be illustrated in figure 8 between the first interface of prediction and second interface
Interbed multiple schematic diagram, Fig. 9 is after being subtracted each other using multiple wave self-adaption after compacting many subwaves relevant with two interfaces
Stacked section schematic diagram.
Based on same inventive concept, the poststack that a kind of many subwaves between two interface are additionally provided in the embodiment of the present invention is pre-
Device is surveyed, as described in the following examples.Because the poststack prediction meanss of many subwaves between two interfaces solve the principle of problem
The poststack Forecasting Methodology of many subwaves between two interfaces is similar, therefore the poststack prediction meanss of many subwaves between two interfaces
Implementation may refer to the implementation of the poststack Forecasting Methodology of many subwaves between two interfaces, repeats part and repeats no more.It is following to be made
, term " unit " or " module " can realize the combination of the software and/or hardware of predetermined function.Although following examples
Described device is preferably realized with software, but hardware, or the combination of software and hardware realization be also may be simultaneously
It is contemplated.Figure 10 is a kind of structured flowchart of the poststack prediction meanss of many subwaves between two interfaces of the embodiment of the present invention,
As shown in Figure 10, including：Stacked section and stacking velocity field acquiring unit 1001, when deep acquiring unit 1002, stacked section draw
Subdivision 1003, searching unit 1004, depth calculation unit 1005, the first continuation unit 1006, the second continuation unit 1007,
Energy calculation unit 1008, cross-correlation section computing unit 1009, normalization unit 1010 and the interbed multiple for shaking road are calculated
Unit 1011, is illustrated to the structure below.
Stacked section and stacking velocity field acquiring unit 1001, for obtaining stacked section and corresponding with the stacked section
Stacking velocity field；
When deep acquiring unit 1002, it is deep when given for obtaining first and when second it is deep, wherein, be more than the when second deeply
Deep for the moment, deeply corresponding subsurface reflective boundary is the first interface when first, and deeply corresponding subsurface reflective boundary is second when second
Interface；
Stacked section division unit 1003, for the stacked section to be divided into Part I data with first interface
With Part II data, wherein, the Part I data be first interface above section data, the Part II
Data are the data of part below first interface；
Searching unit 1004, for finding out deeply corresponding first superposition speed when described first from the stacking velocity field
Degree and deeply corresponding second stack velocity when described second；
Depth calculation unit 1005, described the is obtained for being calculated according to deep when described first and first stack velocity
First depth at one interface, calculates according to deep and second stack velocity when described second and obtains the second of the second contact surface
Depth；
First continuation unit 1006, for according to first depth and first stack velocity, using down going wave side
The stacked section from reference plane continuation to first interface, is obtained the first wave field by journey, then using up-going wave equation, by institute
The first wave field is stated from first interface continuation to the reference plane, the whole process and pegleg multiples at first interface is obtained；
Second continuation unit 1007, for according to second depth and second stack velocity, using down going wave side
The stacked section from reference plane continuation to the second contact surface, is obtained the second wave field by journey, then using up-going wave equation, by institute
The second wave field is stated from the second contact surface continuation to the reference plane, the whole process and pegleg multiples of the second contact surface is obtained；
The energy calculation unit 1008 of seismic channel, is calculated for the whole process and pegleg multiples according to first interface
To the energy of the seismic channel of many subwaves in the first interface, calculated and obtained according to the whole process and pegleg multiples of the second contact surface
The energy of the seismic channel of many subwaves of second contact surface；
Cross-correlation section computing unit 1009, for the whole process and pegleg multiples according to first interface and described the
The whole process and pegleg multiples of second interface, calculating obtain mutual between first interface and many subwaves of the second contact surface
Close section；
Normalization unit 1010, energy, the second contact surface for the seismic channel according to many subwaves in the first interface
The energy of the seismic channel of many subwaves and the cross-correlation section, calculating obtain normalized cross-correlation section；
Interbed multiple computing unit 1011, for according to the normalized related section and the Part II number
According to calculating obtains the interbed multiple between first interface and the second contact surface.
In one embodiment, depth calculation unit 1005 can calculate first depth and the according to below equation
Two depth：
Z1(x, y)=T1(x,y) V1(x,y)
Z2(x, y)=T2(x,y) V2(x,y)
Wherein, V1(x, y)=V (x, y, T1(x, y)), V2(x, y)=V (x, y, T2(x, y)), Z1Represent the first depth, Z2
Represent the second depth, T1Represent deep when first, T2Represent deep when second, V1Represent the first stack velocity, V2Represent the second superposition speed
Degree, V represents the stacking velocity field of the stacked section, and x represents the vertical line direction coordinate of stacked section, and y represents stacked section
Horizontal line direction coordinate.
In one embodiment, the energy calculation unit 1008 of seismic channel can calculate described first according to below equation
The energy of the energy of the seismic channel of many subwaves in interface and the seismic channel of many subwaves of the second contact surface：
Wherein, M1(x, y, t) represents the whole process and pegleg multiples at the first interface, M2(x, y, t) represents the complete of second contact surface
Journey and pegleg multiples, E1(x, y) represents the energy of the seismic channel of many subwaves in the first interface, E2(x, y) represents that second contact surface is multiple
The energy of the seismic channel of ripple, x represents the vertical line direction coordinate of stacked section, and y represents the horizontal line direction coordinate of stacked section,
T represents time coordinate, tmRepresent the dominant record time of stacked section.
In one embodiment, cross-correlation section computing unit 1009 can calculate first boundary according to below equation
Cross-correlation section between face and many subwaves of the second contact surface：
Wherein, R (x, y, t) represents the cross-correlation section between first interface and many subwaves of the second contact surface,
M1(x, y, t) represents the whole process and pegleg multiples at the first interface, M2(x, y, t) represents that the whole and micro- of second contact surface bends multiple
Ripple,Represent computing cross-correlation.
In one embodiment, normalization unit 1010 can calculate normalized cross-correlation according to below equation and cut open
Face：
Wherein, Rn(x, y, t) represents normalized cross-correlation section.
In one embodiment, interbed multiple computing unit 1011 can calculate first boundary according to below equation
Interbed multiple between face and the second contact surface：
M3(x, y, t)=Rn(x,y,t)*D2(x,y,t)
Wherein, M3(x, y, t) represents the interbed multiple between first interface and the second contact surface, D2(x,y,t)
The Part II data are represented, * represents convolution operation.
In another embodiment, a kind of software is additionally provided, the software is used to perform above-described embodiment and preferred real
Apply the technical scheme described in mode.
In another embodiment, a kind of storage medium is additionally provided, be stored with above-mentioned software in the storage medium, should
Storage medium includes but is not limited to：CD, floppy disk, hard disk, scratch pad memory etc..
As can be seen from the above description, the embodiment of the present invention realizes following technique effect：It will need by two
The prediction of the interbed multiple at two interfaces that interfacial layer speed and six wave field extrapolations could be realized, is reduced to only need superposition
The prediction of the interbed multiple at two interfaces that speed and four wave field extrapolations can just be realized, so as to solve in the prior art two
The accuracy of the prediction of the interbed multiple at interface is not high, the low technical problem of forecasting efficiency, has reached and has effectively improved prediction
The technique effect of accuracy and forecasting efficiency.
Obviously, those skilled in the art should be understood that each module or each step of the above-mentioned embodiment of the present invention can be with
Realized with general computing device, they can be concentrated on single computing device, or be distributed in multiple computing devices
On the network constituted, alternatively, the program code that they can be can perform with computing device be realized, it is thus possible to by it
Store and performed in the storage device by computing device, and in some cases, can be to be held different from order herein
They, are either fabricated to each integrated circuit modules or will be multiple in them by the shown or described step of row respectively
Module or step are fabricated to single integrated circuit module to realize.So, the embodiment of the present invention is not restricted to any specific hard
Part and software are combined.
The preferred embodiments of the present invention are the foregoing is only, are not intended to limit the invention, for the skill of this area
For art personnel, the embodiment of the present invention can have various modifications and variations.Within the spirit and principles of the invention, made
Any modification, equivalent substitution and improvements etc., should be included in the scope of the protection.
Claims (10)
1. a kind of poststack Forecasting Methodology of many subwaves between two interface, it is characterised in that including：
Obtain stacked section and stacking velocity field corresponding with the stacked section；
Obtain it is deep during given first and when second it is deep, wherein, deep, deeply corresponding underground when first when being more than first when second deeply
Reflecting interface is the first interface, and deeply corresponding subsurface reflective boundary is second contact surface when second；
The stacked section is divided into by Part I data and Part II data with first interface, wherein, described first
Partial data is the data of first interface above section, and the Part II data are part below first interface
Data；
Corresponded to deeply when deeply corresponding first stack velocity and described second are found out when described first from the stacking velocity field
The second stack velocity；
The first depth for obtaining first interface is calculated according to deep and first stack velocity when described first, according to described
Deep and second stack velocity calculates the second depth for obtaining the second contact surface when second；
According to first depth and first stack velocity, the stacked section is prolonged from reference plane using downgoing wave equation
First interface is opened up, the first wave field is obtained, then using up-going wave equation, first wave field is prolonged from first interface
The reference plane is opened up, the whole process and pegleg multiples at first interface is obtained；
According to second depth and second stack velocity, the stacked section is prolonged from reference plane using downgoing wave equation
The second contact surface is opened up, the second wave field is obtained, then using up-going wave equation, second wave field is prolonged from the second contact surface
The reference plane is opened up, the whole process and pegleg multiples of the second contact surface is obtained；
The energy for the seismic channel for obtaining many subwaves in the first interface is calculated according to the whole process and pegleg multiples at first interface
Amount, the energy for the seismic channel for obtaining many subwaves of the second contact surface is calculated according to the whole process and pegleg multiples of the second contact surface
Amount；
According to the whole process and pegleg multiples and the whole process and pegleg multiples of the second contact surface at first interface, calculate
Cross-correlation section between first interface and many subwaves of the second contact surface；
According to the energy, the energy of the seismic channel of many subwaves of the second contact surface and institute of the seismic channel of many subwaves in the first interface
Cross-correlation section is stated, calculating obtains normalized cross-correlation section；
According to the normalized related section and the Part II data, calculating obtains first interface and described second
Interbed multiple between interface.
2. the method as described in claim 1, it is characterised in that calculate first depth and second deeply according to below equation
Degree：
Z1(x, y)=T1(x,y)V1(x,y)
Z2(x, y)=T2(x,y)V2(x,y)
Wherein, V1(x, y)=V (x, y, T1(x, y)), V2(x, y)=V (x, y, T2(x, y)), Z1Represent the first depth, Z2Represent
Second depth, T1Represent deep when first, T2Represent deep when second, V1Represent the first stack velocity, V2Represent the second stack velocity, V
The stacking velocity field of the stacked section is represented, x represents the vertical line direction coordinate of stacked section, and y represents the horizontal stroke of stacked section
Line direction coordinate.
3. method as claimed in claim 1 or 2, it is characterised in that calculate many subwaves in the first interface according to below equation
Seismic channel energy and many subwaves of the second contact surface seismic channel energy：
<mrow>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
</munderover>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
1
<mrow>
<msub>
<mi>E</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
</munderover>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
Wherein, M1(x, y, t) represents the whole process and pegleg multiples at the first interface, M2(x, y, t) represent second contact surface whole process and
Pegleg multiples, E1(x, y) represents the energy of the seismic channel of many subwaves in the first interface, E2(x, y) represents many subwaves of second contact surface
The energy of seismic channel, x represents the vertical line direction coordinate of stacked section, and y represents the horizontal line direction coordinate of stacked section, t tables
Show time coordinate, tmRepresent the dominant record time of stacked section.
4. method as claimed in claim 3, it is characterised in that calculate first interface and described second according to below equation
Cross-correlation section between many subwaves at interface：
<mrow>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&CircleTimes;</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
Wherein, R (x, y, t) represents the cross-correlation section between first interface and many subwaves of the second contact surface, M1(x,
Y, t) represent the first interface whole process and pegleg multiples, M2(x, y, t) represents the whole process and pegleg multiples of second contact surface,
Represent computing cross-correlation.
5. method as claimed in claim 4, it is characterised in that calculate normalized cross-correlation section according to below equation：
<mrow>
<msub>
<mi>R</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<msqrt>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>E</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
</msqrt>
</mfrac>
</mrow>
Wherein, Rn(x, y, t) represents normalized cross-correlation section.
6. method as claimed in claim 5, it is characterised in that calculate first interface and described second according to below equation
Interbed multiple between interface：
M3(x, y, t)=Rn(x,y,t)*D2(x,y,t)
Wherein, M3(x, y, t) represents the interbed multiple between first interface and the second contact surface, D2(x, y, t) is represented
The Part II data, * represents convolution operation.
7. a kind of poststack prediction meanss of many subwaves between two interface, it is characterised in that including：
Stacked section and stacking velocity field acquiring unit, for obtaining stacked section and superposition speed corresponding with the stacked section
Spend field；
When deep acquiring unit, it is deep when given for obtaining first and when second it is deep, wherein, be more than deeply when second it is deep when first,
Deeply corresponding subsurface reflective boundary is the first interface when first, and deeply corresponding subsurface reflective boundary is second contact surface when second；
Stacked section division unit, for the stacked section to be divided into Part I data and second with first interface
Divided data, wherein, the Part I data are the data of first interface above section, and the Part II data are institute
State the data of part below the first interface；
Searching unit, for finding out when described first deeply corresponding first stack velocity from the stacking velocity field and described
Deeply corresponding second stack velocity when second；
Depth calculation unit, first interface is obtained for being calculated according to deep and first stack velocity when described first
First depth, the second depth for obtaining the second contact surface is calculated according to deep and second stack velocity when described second；
First continuation unit, for according to first depth and first stack velocity, using downgoing wave equation by described in
Stacked section obtains the first wave field from reference plane continuation to first interface, then using up-going wave equation, by the first wave
Field obtains the whole process and pegleg multiples at first interface from first interface continuation to the reference plane；
Second continuation unit, for according to second depth and second stack velocity, using downgoing wave equation by described in
Stacked section obtains the second wave field from reference plane continuation to the second contact surface, then using up-going wave equation, by second ripple
Field obtains the whole process and pegleg multiples of the second contact surface from the second contact surface continuation to the reference plane；
The energy calculation unit of seismic channel, described the is obtained for being calculated according to the whole process and pegleg multiples at first interface
The energy of the seismic channel of many subwaves in one interface, calculates according to the whole process and pegleg multiples of the second contact surface and obtains described second
The energy of the seismic channel of many subwaves in interface；
Cross-correlation section computing unit, for the whole process according to first interface and pegleg multiples and the second contact surface
Whole and pegleg multiples, calculates the cross-correlation section obtained between first interface and many subwaves of the second contact surface；
Normalization unit, energy, many subwaves of the second contact surface for the seismic channel according to many subwaves in the first interface
The energy of seismic channel and the cross-correlation section, calculating obtain normalized cross-correlation section；
Interbed multiple computing unit, for according to the normalized related section and the Part II data, calculating
Interbed multiple between first interface and the second contact surface.
8. device as claimed in claim 7, it is characterised in that the depth calculation unit is specifically for according to below equation meter
Calculate first depth and the second depth：
Z1(x, y)=T1(x,y)V1(x,y)
Z2(x, y)=T2(x,y)V2(x,y)
Wherein, V1(x, y)=V (x, y, T1(x, y)), V2(x, y)=V (x, y, T2(x, y)), Z1Represent the first depth, Z2Represent
Second depth, T1Represent deep when first, T2Represent deep when second, V1Represent the first stack velocity, V2Represent the second stack velocity, V
The stacking velocity field of the stacked section is represented, x represents the vertical line direction coordinate of stacked section, and y represents the horizontal stroke of stacked section
Line direction coordinate.
9. device as claimed in claim 7 or 8, it is characterised in that the energy calculation unit of the seismic channel specifically for by
The seismic channel of energy and many subwaves of the second contact surface according to the seismic channel of the below equation calculating many subwaves in the first interface
Energy：
<mrow>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
</munderover>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>E</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
</munderover>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
Wherein, M1(x, y, t) represents the whole process and pegleg multiples at the first interface, M2(x, y, t) represent second contact surface whole process and
Pegleg multiples, E1(x, y) represents the energy of the seismic channel of many subwaves in the first interface, E2(x, y) represents many subwaves of second contact surface
The energy of seismic channel, x represents the vertical line direction coordinate of stacked section, and y represents the horizontal line coordinates of stacked section, when t is represented
Between coordinate, tmRepresent the dominant record time of stacked section.
10. device as claimed in claim 9, it is characterised in that the cross-correlation section computing unit specifically for according to
Lower formula calculates the cross-correlation section between first interface and many subwaves of the second contact surface：
<mrow>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&CircleTimes;</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
Wherein, R (x, y, t) represents the cross-correlation section between first interface and many subwaves of the second contact surface, M1(x,
Y, t) represent the first interface whole process and pegleg multiples, M2(x, y, t) represents the whole process and pegleg multiples of second contact surface,
Represent computing cross-correlation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510240408.0A CN104834010B (en) | 2015-05-13 | 2015-05-13 | The poststack Forecasting Methodology and device of many subwaves between two interfaces |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510240408.0A CN104834010B (en) | 2015-05-13 | 2015-05-13 | The poststack Forecasting Methodology and device of many subwaves between two interfaces |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104834010A CN104834010A (en) | 2015-08-12 |
| CN104834010B true CN104834010B (en) | 2017-09-12 |
Family
ID=53812006
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510240408.0A Active CN104834010B (en) | 2015-05-13 | 2015-05-13 | The poststack Forecasting Methodology and device of many subwaves between two interfaces |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104834010B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112099088B (en) * | 2020-09-16 | 2022-04-12 | 中油奥博(成都)科技有限公司 | Oil-gas indication and characterization method based on high-density optical fiber seismic data |
| CN114910965A (en) * | 2021-02-07 | 2022-08-16 | 中国石油天然气股份有限公司 | Multiple suppression method and system for post-stack seismic data |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6219621B1 (en) * | 1998-06-30 | 2001-04-17 | Exxonmobil Upstream Research Co. | Sparse hyperbolic inversion of seismic data |
| CN101630019B (en) * | 2008-07-16 | 2011-12-07 | 中国石油集团东方地球物理勘探有限责任公司 | Method for correcting predicated large time difference of multiple benthal waves |
| CN103376466B (en) * | 2012-04-13 | 2016-02-03 | 中国石油化工股份有限公司 | A kind of multiple reflection drawing method |
| CN104280777B (en) * | 2013-07-12 | 2017-02-22 | 中国石油天然气集团公司 | Method for suppressing interference of seismic data multiples on land |
| CN104502974A (en) * | 2014-12-31 | 2015-04-08 | 中国石油天然气集团公司 | Suppressing multiple reflection combining method and device |
-
2015
- 2015-05-13 CN CN201510240408.0A patent/CN104834010B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN104834010A (en) | 2015-08-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hale | Dynamic warping of seismic images | |
| Liu et al. | Least-squares reverse time migration using controlled-order multiple reflections | |
| Yilmaz et al. | Migration velocity analysis by wave-field extrapolation | |
| CN106932824B (en) | The adaptive interlayer multiple suppression method of dimensionality reduction of land seismic prospecting data | |
| US9128205B2 (en) | Process for creating image gathers | |
| CN1036682C (en) | Migration velocity analysis using limited aperture migration | |
| EP3507626B1 (en) | Attenuation of multiple reflections | |
| Wang et al. | Simultaneous reverse time migration of primaries and free-surface related multiples without multiple prediction | |
| CN102819039B (en) | Interval velocity model building method under complicated geological conditions | |
| Liu et al. | Fast least-squares reverse time migration of VSP free-surface multiples with dynamic phase-encoding schemes | |
| US6665618B1 (en) | Seismic survey design technique | |
| US10162070B2 (en) | Converting a first acquired data subset to a second acquired data subset | |
| Zhu et al. | Overthrust imaging with tomo-datuming: A case study | |
| Ma et al. | Wavefield extrapolation in pseudodepth domain | |
| Sahar et al. | Quantification of modification of ground motion due to urbanization in a 3D basin using viscoelastic finite-difference modelling | |
| CN103210323A (en) | System and method for processing seismic data | |
| CN104834010B (en) | The poststack Forecasting Methodology and device of many subwaves between two interfaces | |
| He et al. | Improving the least-squares image by using angle information to avoid cycle skipping | |
| CN104635267B (en) | A kind of seismic event sublevel reverse-time migration weighted stacking formation method | |
| CN107422375B (en) | The determination method and apparatus of the CRP gather of subsea node | |
| US20220342101A1 (en) | Automatic data enhancement for full waveform inversion in the midpoint-offset domain | |
| Dhakal et al. | Validation of the deep velocity structure of the Tokachi basin based on 3-D simulation of long-period ground motions | |
| Pylypenko et al. | Seismic migration in near-vertical and wide-angle reflection and refraction studies: towards a unified approach | |
| Jeong et al. | Effective suppression of spurious events when generating reflected P-and PS-wave data using seismic interferometry | |
| Sivaram | Simulation of Strong Ground Motion in the National Capital Region, India, from a Future 8.5 Magnitude Earthquake Using Two‐Step Empirical Green’s Function Method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |