CN105807315B - Elastic vector reverse-time migration imaging method - Google Patents
Elastic vector reverse-time migration imaging method Download PDFInfo
- Publication number
- CN105807315B CN105807315B CN201610143873.7A CN201610143873A CN105807315B CN 105807315 B CN105807315 B CN 105807315B CN 201610143873 A CN201610143873 A CN 201610143873A CN 105807315 B CN105807315 B CN 105807315B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msubsup
- msub
- vector
- wave
- 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
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. analysis, for interpretation, for correction
Abstract
The invention belongs to exploration geophysics field, in particular it relates to a kind of elastic vector reverse-time migration imaging method.The present invention carries out wave field extrapolation using multi-component seismic data as input, using elastic wave propagation operator, builds elastic vector seismic wave field, and utilize the longitudinal and transverse wave field for decoupling the acquisition decoupling of extension equation method；Elastic wave cross-correlation image-forming condition is built using the vector dot product of the longitudinal and transverse wave field with focus illumination compensation, obtains multi-component seismic data elasticity vector reverse-time migration imaging results；Low frequency noise compacting is carried out to imaging results, obtains final elastic wave reverse-time migration imaging results.The present invention solves the scalar imaging problem of the longitudinal and transverse wave field of vector, improves imaging precision of the elastic reverse-time migration in the case of complicated structure.
Description
Technical field
The invention belongs to exploration geophysics field, in particular it relates to a kind of elastic vector reverse-time migration imaging method.
Background technology
It can be divided into two systems currently for the migration imaging of multi-component seismic data：One is scalar offset imaging body
System, under the system, multi-component seismic data is first broken down into compressional wave data and shear wave data, then carries out PP and PS respectively
Migration imaging；The second is multi -components combine elastic reverse-time migration imaging system, the system directly combines multi-component seismic data
As the boundary value condition of vector elastic wave field, traveling-wave field reconstruct is entered based on Time Migration of Elastic Wave Equation, it is final to obtain component imaging results
Or elastic wave field imaging results, it is a kind of imaging method accurate in theory.Normal conditions, using obtain elastic wave field into picture value as
The multi -components elastic reverse time migration imaging method by combining of final result application image-forming condition before using Helmholtz decompose into
Traveling-wave field is separated, but thus obtained compressional wave and shear wave are become relative to the amplitude and phase for being originally inputted wave field
Change so that wave field separation result protects width deficiency.
In the last few years, decoupling extension equation method was developing progressively as a kind of method of the practical longitudinal and transverse wave field of acquisition.Decoupling
Extension equation method can obtain the compressional wave and shear wave of the amplitude decoupling consistent with phase, but the longitudinal and transverse wave field after decoupling is
Vector field, the imaging results of scalar can not be directly obtained using these vector wave fields, this brings difficulty to following explanations work.
The content of the invention
In order to solve the problems, such as that vector wave field can not directly obtain scalar imaging results, the present invention provides the elastic vector inverse time
Offset imaging method, based on the elastic vector reverse-time migration imaging of elastic vector wave site product cross-correlation, scalar shape can be obtained
The elastic reverse-time migration imaging results of formula.
To achieve the above object, the present invention uses following technical proposals：
Elastic vector reverse-time migration imaging method, step are as follows：
Step 1：According to dielectric model, given source wavelet, realize that wave field just pushes out, using decoupling extension equation method
Obtain the longitudinal and transverse wave field of decoupling；
Step 2：According to dielectric model, using multi-component seismic data as boundary value condition, realize outside the inverse time of seismic wave field
Push away, utilize the longitudinal and transverse wave field for decoupling the acquisition decoupling of extension equation method；
Step 3：It is imaged, is marked using the vector wave site product cross-correlation image-forming condition of source wavefield illumination compensation
The imaging results of amount；
Step 4：Low frequency noise compacting is carried out to the imaging results of acquisition, obtains final imaging results.
It is as follows relative to prior art, beneficial effects of the present invention：It is real to provide a kind of elastic vector reverse-time migration imaging
Applying method, the elastic wave cross-correlation image-forming condition based on vector dot product structure can obtain scalar imaging knot using vector wave field
Fruit, avoid the imaging section as caused by being imaged component and be excessively difficult to the difficulty being explained further.
Brief description of the drawings
Fig. 1 is the schematic flow sheet of elastic vector reverse-time migration imaging method；
Fig. 2-1 is the velocity of longitudinal wave field according to the Marmosi2 elastic fluid models built；
Fig. 2-2 is the shear wave velocity field according to the Marmosi2 elastic fluid models built；
Fig. 3-1 is that forward direction is extrapolated to the 1s moment using decoupling continuation side in 5.8 kms to 12.4 km regions in the x direction
The component of the vector compressional wave of source wavefield after the separation of journey method in the x-direction；
Fig. 3-2 is that forward direction is extrapolated to the 1s moment using decoupling continuation side in 5.8 kms to 12.4 km regions in the x direction
The component of the vector compressional wave of source wavefield after the separation of journey method in the z-direction；
Fig. 3-3 is that forward direction is extrapolated to the 1s moment using decoupling continuation side in 5.8 kms to 12.4 km regions in the x direction
The component of the vector shear wave of source wavefield after the separation of journey method in the x-direction；
Fig. 3-4 is that forward direction is extrapolated to the 1s moment using decoupling continuation side in 5.8 kms to 12.4 km regions in the x direction
The component of the vector shear wave of source wavefield after the separation of journey method in the z-direction；
Fig. 4-1 is that inverse time continuation to the 1s moment uses decoupling continuation side in 5.8 kms to 12.4 km regions in the x direction
The component of the vector compressional wave of detection wave field after the separation of journey method in the x-direction；
Fig. 4-2 is that inverse time continuation to the 1s moment uses decoupling continuation side in 5.8 kms to 12.4 km regions in the x direction
The component of the vector compressional wave of detection wave field after the separation of journey method in the z-direction；
Fig. 4-3 is that inverse time continuation to the 1s moment uses decoupling continuation side in 5.8 kms to 12.4 km regions in the x direction
The component of the vector shear wave of detection wave field after the separation of journey method in the x-direction；
Fig. 4-4 is that inverse time continuation to the 1s moment uses decoupling continuation side in 5.8 kms to 12.4 km regions in the x direction
The component of the vector shear wave of detection wave field after the separation of journey method in the z-direction；
Fig. 5-1 be using vector dot product structure elastic wave cross-correlation image-forming condition obtain 5.8 kms in the x direction extremely
PP imaging sections in 12.4 km regions；
Fig. 5-2 be using vector dot product structure elastic wave cross-correlation image-forming condition obtain 5.8 kms in the x direction extremely
PS imaging sections in 12.4 km regions；
Fig. 5-3 be using vector dot product structure elastic wave cross-correlation image-forming condition obtain 5.8 kms in the x direction extremely
SP imaging sections in 12.4 km regions；
Fig. 5-4 be using vector dot product structure elastic wave cross-correlation image-forming condition obtain 5.8 kms in the x direction extremely
SS imaging sections in 12.4 km regions；
Fig. 6-1 is after final superposition and carries out the PP imaging sections of Laplce's filtering；
Fig. 6-2 is after final superposition and carries out the PS imaging sections of Laplce's filtering；
Fig. 6-3 is after final superposition and carries out the SP imaging sections of Laplce's filtering；
Fig. 6-4 is after final superposition and carries out the SS imaging sections of Laplce's filtering.
Embodiment
As shown in figure 1, elastic vector reverse-time migration method, step are as follows：
Step 1：According to dielectric model, given source wavelet, realize that wave field just pushes out, using decoupling extension equation method
The longitudinal and transverse wave field of decoupling is obtained, specific method is as follows：
(1), according to dielectric model set in advance, source wavelet is loaded, using elastic wave propagation operator, realizes seismic wave
The positive time continuation of field, build source vector wave field：
Utilize following one-order velocity-stress equations for elastic waves
Structure source wavefield is just pushing out staggering mesh finite-difference operator, and wherein ρ is density, and C is elastic fluid rigidity
Matrix, v=(vx,vy,vz)TParticle Vibration Velocity vector field is represented, superscript notation " T " represents transposition, vxRepresent particle vibration speed
Spend the component of vector field in the x-direction, vyRepresent component of the Particle Vibration Velocity vector field along y, vzRepresent Particle Vibration Velocity vector
Component of the field along z, τ=(σxx,σyy,σzz,τyz,τxz,τxy)TIt is stress tensor, σxx、σyyAnd σzzIt is direct stress, τyz、τxzAnd τxy
It is shearing stress,Derivative of the Particle Vibration Velocity vector field on time orientation is represented,Represent stress tensor on time orientation
Derivative, L is differential matrix,
Wherein, lx, lyAnd lzRepresent respectively along the derivative in x, y and z directionss, the stiffness matrix C tables in isotropic medium
It is shown as
Wherein λ and μ is Lame Coefficient.Discretization is carried out to above-mentioned wave equation and obtains following elastic wave continuation operator：
Wherein, τSRepresent the discrete stress field of source wavefield, vSThe discrete Particle Vibration Velocity field of source wavefield is represented, η is border
Absorption coefficient, absorption coefficient η=0 in target area, the absorption coefficient η=200 (0.5-0.5cos (π in the absorption region of border
R/R)), r=1,2 ..., R, R be absorbed layer thickness, π represents pi, and Δ t be time sampling interval, when n Δs t expressions are whole
Between point, (n+1/2) Δ t is half timing node, n=1,2 ..., N, T0The earthquake record that=N Δs t represents total receives duration, DfWith
DbHigh-order staggering mesh finite-difference matrix operator is represented respectively, and expression is
With
WhereinWithForwardly and rearwardly Difference Schemes with Staggered in the x-direction is represented,WithRepresent in the y-direction
Forwardly and rearwardly Difference Schemes with Staggered, whereinWithRepresent forwardly and rearwardly staggered-mesh difference lattice in the z-direction
Formula, specific form are as follows：
Wherein, Δ x, Δ y and Δ z be respectively along the sampling interval in x, y and z directionss,It is 2M
Rank staggering mesh finite-difference coefficient, f (i, j, k) are space networks lattice point (i, j, k) smooth functions.
(2), during source wavefield forward direction time continuation, the source wavefield for decoupling the acquisition decoupling of extension equation method is utilized
P-wave And S data：
Using the P-wave And S field data for decoupling the acquisition decoupling of extension equation method during wave field extrapolation, pass through public affairs first
Formula (3) calculates source wavefield compressional wave stress wave field：
WhereinSource wavefield compressional wave stress field,Derivative of the compressional wave stress wave field on time orientation is represented,Represent
Divergence operator, then the Particle Vibration Velocity field by compressional wave stress field acquisition focus compressional wave：
WhereinSource wavefield longitudinal-wave particle vibration velocity field is represented,Represent source wavefield compressional wave
Derivative of the Particle Vibration Velocity field on time orientation,Represent gradient operator,Represent the particle vibration of source wavefield compressional wave
The component of velocity vector field in the x-direction,Component of the Particle Vibration Velocity vector field of source wavefield compressional wave along y is represented,Table
Show component of the Particle Vibration Velocity vector field of source wavefield compressional wave along z, finally from the total wave field v of focusSIn subtract source wavefield
Longitudinal wave fieldSource wavefield shear wave field is obtained, i.e.,
WhereinThe Particle Vibration Velocity field of source wavefield shear wave is represented,Represent that source wavefield is horizontal
The component of the Particle Vibration Velocity vector field of ripple in the x-direction,Represent the Particle Vibration Velocity vector field of source wavefield shear wave along y
Component,Represent component of the Particle Vibration Velocity vector field of source wavefield shear wave along z.Outside specification word description,
The present invention can also be intuitively embodied in Figure of description, accompanying drawing 3-1 is the speed given to 2-1 and Fig. 2-2 with reference to the accompanying drawings
The focus longitudinal wave field obtained after the separation of the model elastic wave field that 5.8 kms to 12.4 km regions calculate in the x direction is in the x-direction
Component；Accompanying drawing 3-2 is the rate pattern given to 2-1 and Fig. 2-2 with reference to the accompanying drawings 5.8 kms to 12.4 kms in the x direction
The component of the focus longitudinal wave field obtained after the elastic wave field separation that region calculates in the z-direction；Accompanying drawing 3-3 is to 2-1 with reference to the accompanying drawings
Obtained after being separated with the given rate pattern elastic wave fields that 5.8 kms calculate to 12.4 km regions in the x direction of Fig. 2-2
The component of focus shear wave field in the x-direction, accompanying drawing 3-4 be to 2-1 and Fig. 2-2 with reference to the accompanying drawings give rate pattern in the x direction
The component of the focus shear wave field obtained after the elastic wave field separation that 5.8 kms to 12.4 km regions calculate in the z-direction.From accompanying drawing
3-1 to accompanying drawing 3-4 can clearly be seen that wave field is separated and come.
Step 2：According to dielectric model, using multi-component seismic data as boundary value condition, realize outside the inverse time of seismic wave field
Push away, it is as follows using the longitudinal and transverse wave field for decoupling the acquisition decoupling of extension equation method, specific method：
(1), according to dielectric model set in advance, using the multi-component seismic data of surface seismic records as boundary value condition, profit
With elastic wave propagation operator, the inverse time continuation of seismic wave field is realized, builds detection vector wave field：
By multi-component earthquake data as boundary value condition, it is loaded into following inverse time extrapolation equation
Wherein, τRRepresent the discrete stress field of detection wave field, vRThe discrete Particle Vibration Velocity field of detection wave field is represented, η is border
Absorption coefficient, absorption coefficient η=0 in target area, the absorption coefficient η=200 (0.5-0.5cos (π in the absorption region of border
R/R)), r=1,2 ..., R, R be absorbed layer thickness, π represents pi, and Δ t be time sampling interval, when n Δs t expressions are whole
Between point, (n+1/2) Δ t is half timing node, n=1,2 ..., N, T0The earthquake record that=N Δs t represents total receives duration, DfWith
DbHigh-order staggering mesh finite-difference matrix operator is represented respectively, and expression is
With
WhereinWithForwardly and rearwardly Difference Schemes with Staggered in the x-direction is represented,WithRepresent along y side
To forwardly and rearwardly Difference Schemes with Staggered, whereinWithRepresent forwardly and rearwardly staggered-mesh difference in the z-direction
Form, specific form are as follows：
Wherein, Δ x, Δ y and Δ z be respectively along the sampling interval in x, y and z directionss,It is 2M
Rank staggering mesh finite-difference coefficient, f (i, j, k) are space networks lattice point (i, j, k) smooth functions.
(2), during detection wave field inverse time continuation, the detection wave field for decoupling the acquisition decoupling of extension equation method is utilized
P-wave And S data：
Using the P-wave And S field data for decoupling the acquisition decoupling of extension equation method during wave field extrapolation, pass through public affairs first
Formula (7) calculates detection wave field compressional wave stress field：
WhereinDetection wave field compressional wave stress field,Derivative of the detection wave field compressional wave stress field on time orientation is represented,Divergence operator is represented, then the Particle Vibration Velocity field of detection wave field compressional wave is obtained by detection wave field compressional wave stress field：
WhereinDetection wave field longitudinal-wave particle vibration velocity field is represented,Represent detection wave field compressional wave
Derivative of the Particle Vibration Velocity field on time orientation,Represent gradient operator,Represent the particle vibration of detection wave field compressional wave
The component of velocity vector field in the x-direction,Component of the Particle Vibration Velocity vector field of detection wave field compressional wave along y is represented,Table
Show component of the Particle Vibration Velocity vector field of detection wave field compressional wave along z, finally from the total wave field v of detection wave fieldRIn subtract detection
Wave field longitudinal wave fieldDetection wave field shear wave field is obtained, i.e.,
WhereinThe Particle Vibration Velocity field of detection wave field shear wave is represented,Represent that detection wave field is horizontal
The component of the Particle Vibration Velocity vector field of ripple in the x-direction,Represent the Particle Vibration Velocity vector field edge of detection wave field shear wave
Y component,Represent component of the Particle Vibration Velocity vector field of detection wave field shear wave along z.Accompanying drawing 4-1 is that detection wave field passes through
The component of the vector longitudinal wave field that solution decoupling extension equation method obtains in the x-direction；Accompanying drawing 4-2 is detection wave field by decoupling continuation side
The component of the vector longitudinal wave field that journey method obtains in the z-direction；Accompanying drawing 4-3 is the arrow that detection wave field is obtained by decoupling extension equation method
Measure the component of shear wave field in the x-direction, accompanying drawing 4-4 is detection wave field by decoupling vector shear wave field that extension equation method obtains along z
The component in direction.It can clearly be seen that wave field is separated from accompanying drawing 4-1 to accompanying drawing 4-4 to come.
Step 3：It is imaged, is marked using the vector wave site product cross-correlation image-forming condition of source wavefield illumination compensation
The imaging results of amount, specific method are as follows：
The detection wave field decoupling P-wave And S that the source wavefield decoupling P-wave And S obtained using step 1 is obtained with step 2 enters
Go as follows
Vector dot product cross-correlation is imaged, and carries out focus illumination compensation, obtains the multi-component seismic data elasticity vector inverse time
Offset scalar imaging results.Wherein IPPRepresent PP imaging results, IPSIt is PS imaging results, ISPIt is SP imaging results, ISSSS into
As result, t represents time, T0Represent that earthquake record receives duration,Represent that source wavefield longitudinal-wave particle is shaken
Dynamic velocity field,The component of the Particle Vibration Velocity vector field of source wavefield compressional wave in the x-direction is represented,Represent source wavefield
Component of the Particle Vibration Velocity vector field of compressional wave along y,Represent the Particle Vibration Velocity vector field of source wavefield compressional wave along z's
Component,Detection wave field longitudinal-wave particle vibration velocity field is represented,Represent that the particle of detection wave field compressional wave shakes
The component of dynamic velocity vector field in the x-direction,Component of the Particle Vibration Velocity vector field of detection wave field compressional wave along y is represented,
Component of the Particle Vibration Velocity vector field of detection wave field compressional wave along z is represented,Represent source wavefield shear wave
Particle Vibration Velocity field,The component of the Particle Vibration Velocity vector field of source wavefield shear wave in the x-direction is represented,Represent
Component of the Particle Vibration Velocity vector field of source wavefield shear wave along y,Represent the Particle Vibration Velocity arrow of source wavefield shear wave
Component of the field along z is measured,The Particle Vibration Velocity field of detection wave field shear wave is represented,Represent detection wave field
The component of the Particle Vibration Velocity vector field of shear wave in the x-direction,Represent the Particle Vibration Velocity vector field of detection wave field shear wave
Along y component,Component of the Particle Vibration Velocity vector field of detection wave field shear wave along z is represented, symbol represents vector dot product.
Outside specification word description, the present invention can also be intuitively embodied in Figure of description, Fig. 5-1 is utilized based on illumination
The PP imaging sections that the elastic wave cross-correlation image-forming condition of the vector dot product structure of compensation obtains；Fig. 5-2 is utilized based on illumination
The PS imaging sections that the elastic wave cross-correlation image-forming condition of the vector dot product structure of compensation obtains；Fig. 5-3 is utilized based on illumination
The SP imaging sections that the elastic wave cross-correlation image-forming condition of the vector dot product structure of compensation obtains；Fig. 5-4 is utilized based on illumination
The SS imaging sections that the elastic wave cross-correlation image-forming condition of the vector dot product structure of compensation obtains.
Step 4：Low frequency noise compacting is carried out to the imaging results of acquisition, specific implementation method is as follows：
Utilize following Laplce's filtering algorithm
Processing is filtered to imaging section, obtains final imaging results.WhereinRepresent Laplace operator,WithIt is to carry out filtered PP, PS, SP and SS imaging section of Laplce respectively.Fig. 6-1 is most
The superposition at end simultaneously carries out the PP imaging sections of Laplce's filtering.Fig. 6-2 is final superposition and carries out Laplce's filtering
PS imaging sections, Fig. 6-3 is final superposition and carries out the SP imaging sections of Laplce's filtering, and Fig. 6-4 is final superposition
And carry out the SS imaging sections of Laplce's filtering.
Claims (1)
1. elastic vector reverse-time migration imaging method, it is characterised in that step is as follows：
Step 1：According to dielectric model, given source wavelet, realize that wave field just pushes out, obtained using extension equation method is decoupled
The longitudinal and transverse wave field of vector of decoupling；Concrete methods of realizing is as follows：
(1), according to dielectric model set in advance, source wavelet is loaded, using elastic wave propagation operator, realizes seismic wave field
Positive time continuation, build source vector wave field；
(2), during source wavefield forward direction time continuation, the vector for decoupling the acquisition source wavefield decoupling of extension equation method is utilized
P-wave And S field data；
Step 2：According to dielectric model, using multi-component seismic data as boundary value condition, the inverse time extrapolation of seismic wave field, profit are realized
The longitudinal and transverse wave field of vector of decoupling is obtained with decoupling extension equation method；Concrete methods of realizing is as follows：
(1), according to dielectric model set in advance, using the multi-component seismic data of surface seismic records as boundary value condition, bullet is utilized
Property ripple propagation operator, realize the inverse time continuation of seismic wave field, build detection vector wave field；
(2), during detection wave field inverse time continuation, indulged using the vector for decoupling the acquisition detection wave field decoupling of extension equation method,
Shear wave field data；
Step 3：It is imaged using the vector wave site product cross-correlation image-forming condition of source wavefield illumination compensation, obtains scalar
Imaging results；Concrete methods of realizing is that the longitudinal and transverse wave field of vector of the source wavefield decoupling obtained using step 1 is obtained with step 2
The longitudinal and transverse wave field of vector of detection wave field decoupling carry out vector dot product computing cross-correlation, obtain multi-component seismic data elasticity arrow
Measure reverse-time migration scalar imaging results；The following institute of vector wave site product cross-correlation image-forming condition with source wavefield illumination compensation
Show：
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>P</mi>
<mi>P</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<msub>
<mi>T</mi>
<mn>0</mn>
</msub>
</munderover>
<msubsup>
<mi>v</mi>
<mi>P</mi>
<mi>S</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>v</mi>
<mi>P</mi>
<mi>R</mi>
</msubsup>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<msub>
<mi>T</mi>
<mn>0</mn>
</msub>
</munderover>
<msubsup>
<mi>v</mi>
<mi>P</mi>
<mi>S</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>v</mi>
<mi>P</mi>
<mi>S</mi>
</msubsup>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>P</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<msub>
<mi>T</mi>
<mn>0</mn>
</msub>
</munderover>
<msubsup>
<mi>v</mi>
<mi>P</mi>
<mi>S</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>v</mi>
<mi>S</mi>
<mi>R</mi>
</msubsup>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<msub>
<mi>T</mi>
<mn>0</mn>
</msub>
</munderover>
<msubsup>
<mi>v</mi>
<mi>P</mi>
<mi>S</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>v</mi>
<mi>P</mi>
<mi>S</mi>
</msubsup>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>S</mi>
<mi>P</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<msub>
<mi>T</mi>
<mn>0</mn>
</msub>
</munderover>
<msubsup>
<mi>v</mi>
<mi>S</mi>
<mi>S</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>v</mi>
<mi>P</mi>
<mi>R</mi>
</msubsup>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<msub>
<mi>T</mi>
<mn>0</mn>
</msub>
</munderover>
<msubsup>
<mi>v</mi>
<mi>S</mi>
<mi>S</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>v</mi>
<mi>S</mi>
<mi>S</mi>
</msubsup>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<msub>
<mi>T</mi>
<mn>0</mn>
</msub>
</munderover>
<msubsup>
<mi>v</mi>
<mi>S</mi>
<mi>S</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>v</mi>
<mi>S</mi>
<mi>R</mi>
</msubsup>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<msub>
<mi>T</mi>
<mn>0</mn>
</msub>
</munderover>
<msubsup>
<mi>v</mi>
<mi>S</mi>
<mi>S</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>v</mi>
<mi>S</mi>
<mi>S</mi>
</msubsup>
</mrow>
</mfrac>
</mrow>
Wherein：IPPRepresent PP imaging results, IPSIt is PS imaging results, ISPIt is SP imaging results, ISSIt is SS imaging results, t is represented
Time, T0Represent that earthquake record receives duration,Source wavefield longitudinal-wave particle vibration velocity field is represented,Represent that detection wave field is indulged
Ripple Particle Vibration Velocity field,Source wavefield shear wave Particle Vibration Velocity field is represented,Represent detection wave field shear wave particle vibration
Velocity field, symbol represent vector dot product；
Step 4：Low frequency noise compacting is carried out to the imaging results of acquisition, obtains final imaging results；Specific implementation is using drawing
This filtering algorithm of pula is filtered processing to imaging section, obtains final imaging results.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610143873.7A CN105807315B (en) | 2016-03-14 | 2016-03-14 | Elastic vector reverse-time migration imaging method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610143873.7A CN105807315B (en) | 2016-03-14 | 2016-03-14 | Elastic vector reverse-time migration imaging method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105807315A CN105807315A (en) | 2016-07-27 |
| CN105807315B true CN105807315B (en) | 2018-03-20 |
Family
ID=56467364
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610143873.7A Active CN105807315B (en) | 2016-03-14 | 2016-03-14 | Elastic vector reverse-time migration imaging method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105807315B (en) |
Families Citing this family (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106842300B (en) * | 2016-12-21 | 2018-10-30 | 中国石油大学(华东) | A kind of high efficiency multi-component seismic data true amplitude migration imaging method |
| CN106772585B (en) * | 2017-01-26 | 2018-11-09 | 中国科学院地质与地球物理研究所 | A kind of quasi- analysis method and device of the optimization decoupling equation based on elastic wave |
| CN106970416B (en) * | 2017-03-17 | 2018-12-04 | 中国地质科学院地球物理地球化学勘查研究所 | Elastic wave least square reverse-time migration system and method based on wave field separation |
| CN106950598B (en) * | 2017-03-27 | 2019-01-29 | 中国科学院地质与地球物理研究所 | A kind of migration velocity field method for evaluating reliability |
| CN106896403B (en) * | 2017-05-05 | 2019-09-13 | 中国石油天然气集团有限公司 | Elastic Gaussian beam offset imaging method and system |
| CN106932820A (en) * | 2017-05-08 | 2017-07-07 | 厦门大学 | ACOUSTIC WAVE EQUATION reverse-time migration imaging method based on time domain puppet spectral method |
| CN107179551B (en) * | 2017-06-19 | 2018-04-06 | 吉林大学 | A kind of method using microseism record to subsurface structure direct imaging |
| CN108037526B (en) * | 2017-11-23 | 2018-12-07 | 中国石油大学(华东) | Reverse-time migration method based on all-wave wave field VSP/RVSP seismic data |
| CN109100784B (en) * | 2018-06-08 | 2020-04-28 | 恒泰艾普（北京）能源科技研究院有限公司 | Three-dimensional VSP source detection interchange full-wave-field imaging method |
| CN109143339B (en) * | 2018-08-14 | 2020-06-09 | 中国石油天然气集团有限公司 | Elastic reverse time migration imaging method and device based on transverse wave stress invariant |
| CN109581484B (en) * | 2018-11-01 | 2020-06-09 | 中国石油天然气集团有限公司 | Method and device for analyzing total illumination index of seismic longitudinal wave component multiples |
| US11086036B2 (en) | 2019-01-22 | 2021-08-10 | Saudi Arabian Oil Company | AVO imaging condition in elastic reverse time migration |
| CN111221037B (en) * | 2020-01-21 | 2021-07-23 | 中国石油大学(华东) | Decoupling elastic reverse time migration imaging method and device |
| CN111239804B (en) * | 2020-02-12 | 2021-07-02 | 中国石油大学(华东) | Elastic energy reverse time migration imaging method, device, equipment and system |
| CN111948710A (en) * | 2020-08-05 | 2020-11-17 | 河海大学 | Method for performing elastic vector wave imaging by using seabed multi-component seismic record |
| CN112327359B (en) * | 2020-10-14 | 2022-06-14 | 山东省科学院海洋仪器仪表研究所 | Elastic reverse time migration method based on imaging energy flow vector |
| CN112904426B (en) * | 2021-03-27 | 2022-09-30 | 中国石油大学(华东) | Decoupling elastic wave reverse time migration method, system and application |
| CN113406698B (en) * | 2021-05-24 | 2023-04-21 | 中国石油大学(华东) | Double-phase medium elastic wave reverse time migration imaging method based on longitudinal and transverse wave decoupling |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102156296A (en) * | 2011-04-19 | 2011-08-17 | 中国石油大学(华东) | Elastic reverse time migration imaging method by combining seismic multi-component |
| WO2013033651A1 (en) * | 2011-08-31 | 2013-03-07 | Spectraseis Ag | Full elastic wave equation for 3d data processing on gpgpu |
| CN103149585B (en) * | 2013-01-30 | 2016-02-17 | 中国石油天然气集团公司 | A kind of resilient bias seismic wave field construction method and device |
-
2016
- 2016-03-14 CN CN201610143873.7A patent/CN105807315B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105807315A (en) | 2016-07-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105807315B (en) | Elastic vector reverse-time migration imaging method | |
| US8953411B2 (en) | Apparatus and method for imaging a subsurface using frequency-domain elastic reverse-time migration | |
| Ladopoulos | Petroleum and gas reserves exploration by real-time expert seismology and non-linear seismic wave motion' | |
| CN103149585B (en) | A kind of resilient bias seismic wave field construction method and device | |
| CN112327358B (en) | Forward modeling method for acoustic seismic data in viscous medium | |
| CN107678062A (en) | The integrated forecasting deconvolution of hyperbolic Radon domains and feedback loop methodology multiple suppression model building method | |
| CN105242305A (en) | Longitudinal wave and transverse wave separation method and system | |
| CN109143340B (en) | Viscoelastic medium seismic wave simulation method and system based on constant Q model | |
| Fang et al. | Elastic full-waveform inversion based on GPU accelerated temporal fourth-order finite-difference approximation | |
| Trinh et al. | Efficient 3D elastic FWI using a spectral-element method on Cartesian-based mesh | |
| CN105652319A (en) | Estimation method of near-surface stratum Q value of complex mediums | |
| Chen et al. | Full waveform inversion based on time-integral-damping wavefield | |
| CN111505714B (en) | Elastic wave direct envelope inversion method based on rock physical constraint | |
| Zhang et al. | A depth variant seismic wavelet extraction method for inversion of poststack depth-domain seismic data | |
| CN111257930B (en) | Visco-elastic anisotropic double-phase medium area variable grid solving operator | |
| CN105182414A (en) | Direct wave removing method based on wave equation forward modeling | |
| Zhao et al. | Double plane wave least squares reverse time migration | |
| Xu et al. | A pure viscoacoustic equation for VTI media applied in anisotropic RTM | |
| Manning | Techniques to enhance the accuracy and efficiency of finite-difference modelling for the propagation of elastic waves | |
| Bhaumik et al. | Higher-order thin layer method (HTLM) based wavefield modeling approach | |
| Doyon et al. | Practical aspects of 2.5 D frequency-domain finite-difference modelling of viscoelastic waves | |
| CN108680957B (en) | Local cross-correlation time-frequency domain Phase-retrieval method based on weighting | |
| Chen et al. | Robust iterative deblending of simultaneous-source data | |
| Li et al. | Separation of primaries and different-order multiples based on focal transform and SRME | |
| Xie et al. | Full waveform inversion based on hybrid gradient |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20190401 Address after: 266580 No. 66 Changjiang West Road, Huangdao District, Qingdao, Shandong. Co-patentee after: Large department (Beijing) Technology Co. Ltd. Patentee after: China Petroleum University (East China) Address before: 266580 No. 66 Changjiang West Road, Huangdao economic and Technological Development Zone, Qingdao, Shandong Patentee before: China Petroleum University (East China) |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20211110 Address after: 266580 No. 66, Changjiang West Road, Huangdao District, Qingdao City, Shandong Province Patentee after: China University of Petroleum (East China) Address before: 266580 No. 66, Changjiang West Road, Huangdao District, Qingdao City, Shandong Province Patentee before: China University of Petroleum (East China) Patentee before: Dida Huineng (Beijing) Technology Co., Ltd |