CN111239804B

CN111239804B - Elastic energy reverse time migration imaging method, device, equipment and system

Info

Publication number: CN111239804B
Application number: CN202010088034.6A
Authority: CN
Inventors: 杜启振; 张晓语
Original assignee: China University of Petroleum East China
Current assignee: China University of Petroleum East China
Priority date: 2020-02-12
Filing date: 2020-02-12
Publication date: 2021-07-02
Anticipated expiration: 2040-02-12
Also published as: CN111239804A

Abstract

The embodiment of the specification discloses an elastic energy reverse time migration imaging method, device, equipment and system. The method comprises the steps of decoupling a forward continuation operator based on a constructed mode to obtain vector field information of a seismic source end; constructing a seismic source end energy vector field by using seismic source end vector field information; based on the constructed mode decoupling reverse continuation operator, obtaining vector field information of a detection end; constructing a detection end energy vector field by using detection end vector field information; imaging a seismic source end energy vector field and a detection end energy vector field according to the constructed mode decoupling elastic energy cross-correlation imaging condition to obtain a single shot mode decoupling energy imaging result; and performing multi-shot superposition on the energy imaging result decoupled in the single shot mode to obtain an elastic energy reverse time migration imaging result. The embodiment of the specification can effectively reduce the problem of energy coupling in the imaging result, reduce the difficulty of interpretation work of the subsequent energy imaging result and improve the imaging efficiency.

Description

Elastic energy reverse time migration imaging method, device, equipment and system

Technical Field

The embodiment scheme of the specification belongs to the field of exploration geophysics, and particularly relates to an elastic energy reverse time migration imaging method, device, equipment and system.

Background

The elastic reverse time migration method can utilize multi-component seismic data to carry out structural imaging on the complex earth medium. Elastic reverse time migration methods based on elastic wave equations for wave field construction may require scalar imaging results to be extracted from the vector elastic wave field. In the early stage, most of elastic reverse time migration methods are scalar imaging result extraction based on a vector elastic particle vibration velocity field or a displacement field, and scalar imaging result extraction based on stress is developed subsequently.

In recent years, the elastic reverse time migration method is further developed to realize scalar quantization representation of an elastic vector field by taking the elastic energy density as a field function, and extract a scalar energy imaging result. However, the elastic wave field energy reverse time migration method developed based on the first-order elastic wave equation can only provide the total elastic energy imaging result, and the energy coupling problem exists in the total elastic energy imaging result, which may cause difficulty in the subsequent interpretation work by using the energy imaging result of the pure wave mode and the converted wave energy imaging result.

Disclosure of Invention

Embodiments of the present disclosure provide an elastic energy reverse time migration imaging method, apparatus, device, and system, which can effectively reduce the problem of energy coupling in the elastic wave field energy reverse time migration imaging result, reduce the difficulty in performing interpretation work by using the energy imaging result in the pure wave mode and the converted wave energy imaging result, and improve the imaging efficiency.

The elastic energy reverse time migration imaging method, device, equipment and system provided by the specification are realized by the following modes:

an elastic energy reverse time migration imaging method comprises the following steps:

based on the constructed mode decoupling forward continuation operator, obtaining vector field information of a seismic source end, wherein the vector field information comprises a particle vibration velocity vector field and a pseudo stress vector field of longitudinal waves, and a particle vibration velocity vector field and a pseudo stress vector field of transverse waves;

constructing a seismic source end energy vector field by using the seismic source end vector field information;

based on the constructed mode decoupling reverse continuation operator, obtaining vector field information of a detection end;

constructing a detection end energy vector field by using the detection end vector field information;

imaging the seismic source end energy vector field and the detection end energy vector field according to the constructed mode decoupling elastic energy cross-correlation imaging condition to obtain a single shot mode decoupling energy imaging result;

and performing multi-shot superposition on the energy imaging result decoupled in the single shot mode to obtain an elastic energy reverse time migration imaging result.

An elastic energy reverse time migration imaging device, comprising:

the seismic source end vector field information acquisition module is used for decoupling the forward continuation operator based on the constructed mode to acquire seismic source end vector field information, wherein the vector field information comprises a particle vibration velocity vector field and a pseudo stress vector field of longitudinal waves and a particle vibration velocity vector field and a pseudo stress vector field of transverse waves;

the seismic source end energy vector field construction module is used for constructing a seismic source end energy vector field by utilizing the seismic source end vector field information;

the detection end vector field information acquisition module is used for decoupling a reverse continuation operator based on the constructed mode to acquire detection end vector field information;

the detection end energy vector field construction module is used for constructing a detection end energy vector field by using the detection end vector field information;

the single-shot energy imaging module is used for imaging the seismic source end energy vector field and the detection end energy vector field according to the constructed mode decoupling elastic energy cross-correlation imaging condition to obtain a single-shot mode decoupling energy imaging result;

and the elastic energy reverse time migration imaging module is used for performing multi-shot superposition on the single-shot mode decoupled energy imaging result to obtain an elastic energy reverse time migration imaging result.

A resilient energy reverse time shift imaging device comprising a processor and a memory for storing processor-executable instructions which when executed by the processor implement steps comprising:

A resilient energy reverse time migration imaging system comprising at least one processor and a memory storing computer-executable instructions that, when executed by the processor, perform the steps of the method of any one of the method embodiments of the present specification.

The specification provides an elastic energy reverse time migration imaging method, device, equipment and system. In some embodiments, the mode decoupling forward continuation operator and the mode decoupling reverse continuation operator are constructed on the basis of decoupling continuation, so that mode decoupling of an elastic wave field particle vibration velocity field and a stress field can be realized while forward continuation and reverse continuation of the elastic wave field are realized, and the energy coupling problem in an elastic wave field energy reverse time migration imaging result can be effectively reduced. The constructed mode decoupling energy cross-correlation imaging condition is utilized to obtain the energy imaging result of the pure wave mode and the converted wave energy imaging result from the decoupled elastic energy field, so that the difficulty of interpretation work by utilizing the energy imaging result of the pure wave mode and the converted wave energy imaging result in the follow-up process can be reduced, and the imaging efficiency is improved. By adopting the implementation scheme provided by the specification, the problem of energy coupling in the elastic wave field energy reverse time migration imaging result can be effectively solved, the difficulty in subsequent interpretation work by utilizing the pure wave mode energy imaging result and the converted wave energy imaging result is reduced, and the imaging efficiency is improved.

Drawings

In order to more clearly illustrate the embodiments of the present specification or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, it is obvious that the drawings in the following description are only some embodiments described in the present specification, and for those skilled in the art, other drawings can be obtained according to the drawings without any creative effort.

FIG. 1 is a schematic flow chart diagram of one embodiment of an elastic energy reverse time shift imaging method provided herein;

FIG. 2a is a schematic diagram of a longitudinal wave velocity field of a concave model;

FIG. 2b is a schematic diagram of a shear velocity field of a dimple model;

FIG. 2c is a schematic illustration of a dimple pattern density field;

FIG. 3a is a schematic diagram of a seismic source end longitudinal wave pseudo normal stress field along the xx direction, which is constructed according to the model given in FIG. 2a, FIG. 2b and FIG. 2c by using a mode decoupling forward continuation operator to forward extrapolate to the time of 0.8 s;

FIG. 3b is a schematic diagram of a seismic source end shear wave pseudo-normal stress field in the xx direction constructed by forward extrapolation of a mode decoupling forward continuation operator to 0.8s according to the model given in FIG. 2a, FIG. 2b and FIG. 2 c;

FIG. 4a is a schematic illustration of the components of a 2.7s multicomponent seismic recording acquired at the surface along the x-direction;

FIG. 4b is a schematic component view in the z-direction of a 2.7s multicomponent seismic recording acquired at the surface;

FIG. 5a is a schematic cross-sectional view of a final PP wave energy imaging of a multi-shot stack;

FIG. 5b is a schematic cross-sectional view of the resulting superimposed PS wave energy imaging;

FIG. 5c is a schematic cross-sectional view of the resulting superimposed SP wave energy imaging;

FIG. 5d is a schematic cross-sectional view of the resulting superimposed SS wave energy imaging;

FIG. 6 is a block diagram illustrating an embodiment of an elastic energy reverse time migration imaging device provided herein;

fig. 7 is a block diagram of a hardware configuration of an embodiment of a resilient energy reverse time migration imaging server provided in the present specification.

Detailed Description

In order to make those skilled in the art better understand the technical solutions in the present specification, the technical solutions in the embodiments of the present specification will be clearly and completely described below with reference to the drawings in the embodiments of the present specification, and it is obvious that the described embodiments are only a part of the embodiments in the present specification, and not all of the embodiments. All other embodiments that can be obtained by a person skilled in the art on the basis of one or more embodiments of the present description without inventive step shall fall within the scope of protection of the embodiments of the present description.

The elastic wave field energy reverse time migration method developed based on the first-order elastic wave equation can only provide a total elastic energy imaging result, and the total elastic energy imaging result has an energy coupling problem, so that the difficulty is provided for the subsequent interpretation work by using the energy imaging result of a pure wave mode and the converted wave energy imaging result.

The specification provides an elastic energy reverse time migration imaging method, device, equipment and system. By constructing the mode decoupling forward continuation operator and the mode decoupling reverse continuation operator on the basis of decoupling continuation, the mode decoupling of the elastic wave field particle vibration velocity field and the stress field can be realized while the forward continuation and the reverse continuation of the elastic wave field are realized, and therefore the energy coupling problem in the elastic wave field energy reverse time migration imaging result can be effectively reduced. The constructed mode decoupling energy cross-correlation imaging condition is utilized to obtain the energy imaging result of the pure wave mode and the converted wave energy imaging result from the decoupled elastic energy field, so that the difficulty of interpretation work by utilizing the energy imaging result of the pure wave mode and the converted wave energy imaging result in the follow-up process can be reduced, and the imaging efficiency is improved.

The following describes an embodiment of the present disclosure with a specific application scenario as an example. Specifically, fig. 1 is a schematic flow chart diagram of an embodiment of an elastic energy reverse time shift imaging method provided in the present specification. Although the present specification provides the method steps or apparatus structures as shown in the following examples or figures, more or less steps or modules may be included in the method or apparatus structures based on conventional or non-inventive efforts. In the case of steps or structures which do not logically have the necessary cause and effect relationship, the execution order of the steps or the block structure of the apparatus is not limited to the execution order or the block structure shown in the embodiments or the drawings of the present specification. When the described method or module structure is applied to a device, a server or an end product in practice, the method or module structure according to the embodiment or the figures may be executed sequentially or in parallel (for example, in a parallel processor or multi-thread processing environment, or even in an implementation environment including distributed processing and server clustering).

It should be noted that the following description of the embodiments does not limit the technical solutions in other extensible application scenarios based on the present specification. In one embodiment, as shown in fig. 1, the present disclosure provides an elastic energy reverse time shift imaging method, which may include:

s0: and (3) decoupling a forward continuation operator based on the constructed mode to obtain vector field information of the seismic source end, wherein the vector field information comprises a particle vibration velocity vector field and a pseudo stress vector field of longitudinal waves and a particle vibration velocity vector field and a pseudo stress vector field of transverse waves.

The source-end vector field information may include a particle vibration velocity vector field and a pseudo stress vector field of a source-end longitudinal wave, and a particle vibration velocity vector field and a pseudo stress vector field of a source-end transverse wave. In the embodiment of the description, the mode decoupling prolongation operator can realize the separation of the vibration velocity of longitudinal and transverse wave particles and the separation of longitudinal and transverse wave stress.

In an embodiment of the present specification, a given medium model and a given seismic source wavelet may be utilized, and a forward continuation operator is decoupled according to a constructed mode, so as to realize forward extrapolation of an elastic wave field, and obtain a longitudinal wave particle vibration velocity vector field-pseudo stress vector field and a transverse wave particle vibration velocity vector field-pseudo stress vector field decoupled at a seismic source end.

In an embodiment of the present specification, constructing the pattern decoupling forward continuation operator may include: loading a seismic source wavelet by utilizing a preset medium model; loading the seismic source wavelet to the longitudinal wave body stress of an elastic wave field, and constructing a mode decoupling forward continuation equation; discretizing the mode decoupling forward continuation equation to obtain a mode decoupling forward continuation operator.

For example, in some implementations, a source wavelet may be loaded according to a given media model. Wherein the source wavelet may be represented as:

F_p＝h(t)Φ(x,y,z) (1)

wherein, F_pRepresents the source wavelet, h (t) represents the seismic wavelet function, and Φ (x, y, z) represents the source spatial attenuation function. In some implementations, h (t) may be a rake wavelet h (t) ═ 1-2(π f)_pt)²]exp[-(πf_pt)²]Other wavelets may be used, and this specification does not limit this, where pi represents a circumferential ratio, and f represents_pIndicating dominant frequency, t time, exp exponential function. In some implementations, Φ (x, y, z) ═ exp { - α²[(x-x_s)²+(y-y_s)²+(z-z_s)²]α denotes the attenuation parameter, (x)_s,y_s,z_s) Representing the source location.

In some implementation scenarios, the source wavelet may be loaded on the longitudinal wave body stress of the elastic wave field, and the mode decoupling forward continuation equation is constructed as follows:

wherein, tau_pDenotes the longitudinal wave bulk stress, the superscript "·" denotes the derivative of the respective field function in the time direction, the subscript "p" denotes the longitudinal wave mode, λ and μ denote the Lame coefficient of the background medium,. v ═ denotes the divergence operator, v ═ denotes (v ═ v-_x,v_y,v_z)^TRepresenting the particle vibration velocity vector field, v_x、v_yAnd v_zRepresenting the components of the particle vibration velocity vector field in the x, y and z directions, τ, respectively_s＝(τ_s,xx,τ_s,yy,τ_s,zz,τ_s,xy,τ_s,xz,τ_s,yz)^TShowing a crossWave stress vector field, τ_s,xx、τ_s,yyAnd τ_s,zzDenotes the positive stress of the transverse wave in the x-direction, in the y-direction and in the z-direction, respectively, τ_s,yz、τ_s,xzAnd τ_s,xyRespectively, shear stress of transverse waves in yz direction, xz direction and xy direction, the subscript "s" represents transverse wave mode, C_sRepresenting the stiffness matrix in the transverse wave mode of the background medium, L representing the differential matrix, L_x、l_yAnd l_zRepresenting the spatial derivatives in the x, y and z directions, respectively, the superscript symbol "T" representing transposition, τ_qp＝(τ_qp,xx,τ_qp,yy,τ_qp,zz,τ_qp,yz,τ_qp,xy,τ_qp,xz)^TRepresenting a longitudinal wave pseudo-stress vector field, tau_qp,xx、τ_qp,yyAnd τ_qp,zzDenotes the longitudinal wave pseudo-normal stress in the x-direction, in the y-direction and in the z-direction, respectively, τ_qp,yz、τ_qp,xzAnd τ_qp,xyRespectively representing the longitudinal wave pseudo-shear stress in the yz direction, in the xz direction and in the xy direction, the subscript "qp" representing the longitudinal wave pseudo-stress pattern, τ_qs＝(τ_qs,xx,τ_qs,yy,τ_qs,zz,τ_qs,yz,τ_qs,xy,τ_qs,xz)^TRepresenting a shear wave pseudo-stress vector field, τ_qs,xx、τ_qs,yyAnd τ_qs,zzDenotes the shear wave pseudo-normal stress in the x-direction, in the y-direction and in the z-direction, respectively, τ_qs,yz、τ_qs,xzAnd τ_qs,xyRespectively representing shear wave pseudo-shear stress in yz direction, xz direction and xy direction, the subscript "qs" representing shear wave pseudo-stress mode, C representing stiffness matrix of background medium, ρ representing density of background medium, v_p＝(v_p,x,v_p,y,v_p,z)^TRepresenting the vibration velocity vector field of longitudinal wave particles, v_p,x、v_p,yAnd v_p,zRepresents components of the vibration velocity of the longitudinal wave particle along the x, y and z directions, respectively,. represents a gradient operator,. v_s＝(v_s,x,v_s,y,v_s,z)^TRepresenting the vibration velocity vector field of shear wave particle, v_s,x、v_s,yAnd v_s,zRespectively showing the transverse wave qualityThe components of the point vibration velocity in the x, y and z directions. In some of the implementation scenarios, it is preferred that,

in some implementation scenarios, discretizing the mode decoupling forward continuation equation can obtain the following mode decoupling forward continuation operators:

where the superscript "S" represents the source end wavefield function, and, correspondingly,

showing the stress field of the longitudinal wave body at the seismic source end,

represents the transverse wave stress vector field at the seismic source end,

and

respectively representing the positive shear wave stress at the source end in the x direction, in the y direction and in the z direction,

and

respectively represents the shear stress of the source end along the yz direction, the xz direction and the xy direction,

representing a longitudinal wave pseudo stress vector field at the seismic source end,

and

respectively represents the source end longitudinal wave pseudo-normal stress along the x direction, the y direction and the z direction,

and

respectively represents the pseudo shear stress of the longitudinal wave at the seismic source end along the yz direction, the xz direction and the xy direction, the subscript mark 'qp' represents the pseudo stress mode of the longitudinal wave,

representing the source end shear wave pseudo stress vector field,

and

respectively representing source-end shear pseudo-normal stresses in the x-direction, in the y-direction and in the z-direction,

and

respectively represents the shear wave pseudo shear stress at the seismic source end along the yz direction, the xz direction and the xy direction, the subscript mark 'qs' represents the shear wave pseudo stress mode,

representing the vibration velocity vector field of longitudinal wave particles at the seismic source end,

and

respectively representing components of the vibration velocity of longitudinal wave particles at the seismic source end along the directions of x, y and z,

representing the vibration velocity vector field of the seismic source end shear wave particle,

and

respectively representing components of the seismic source end shear particle vibration velocity in the x, y and z directions,

representing the vibration velocity field of elastic wave particles at the seismic source end,

and

representing components of vibration velocity of elastic wave mass points at a seismic source end along the directions of x, y and z, eta represents a boundary absorption coefficient, delta T represents a time sampling interval, N delta T represents a whole time point corresponding to a discrete point N, (N +1/2) delta T represents a half-time node corresponding to a discrete point (N +1/2), N is 1,2₀N Δ t represents the total seismic recording reception time,

and

respectively, representing a high-order staggered grid finite difference matrix operator.

In some embodiments, the absorption coefficient η in the target region is 0, the absorption coefficient η in the boundary absorption region is 200(0.5-0.5cos (rr R)), R is 1, 2.

In some of the implementation scenarios, it is preferred that,

and

the specific expressions of (a) may be:

and

wherein the content of the first and second substances,

and

respectively representing the forward and backward staggered mesh differential format in the x-direction,

and

respectively representing the forward and backward staggered mesh differential format in the y-direction,

and

the forward and backward staggered mesh differential format along the z-direction is respectively represented, and the specific format is as follows:

where Δ x, Δ y, and Δ z represent sampling intervals in the x, y, and z directions, respectively,

represents a 2M-order interleaved grid finite difference coefficient, u (x)_i,y_j,z_k) Representing spatial grid points (x)_i,y_j,z_k) Elastic wave field function of.

In some embodiments, after obtaining the mode decoupling forward continuation operator, the source-side vector field information may be obtained based on the constructed mode decoupling forward continuation operator. For example, in some implementation scenarios, the forward continuation operator may be decoupled based on the constructed pattern, and the decoupled seismic source-side longitudinal wave particle vibration velocity, pseudo stress, and shear wave particle vibration velocity, pseudo stress may be obtained. In some embodiments, the elastic stress may be decomposed into a longitudinal wave bulk stress and a shear wave stress vector, but in some embodiments, the shear wave stress vector may contain a strong longitudinal wave component, and is referred to herein as a longitudinal wave pseudo stress vector and a shear wave pseudo stress vector in order to distinguish the longitudinal wave bulk stress from the shear wave stress vector.

For example, fig. 2a, 2b and 2c show given medium models, fig. 3a and 3b show the source-side vector field information obtained by the given medium models, wherein FIG. 2a is a schematic diagram of a longitudinal wave velocity field of a concave model, FIG. 2b is a schematic diagram of a transverse wave velocity field of a concave model, FIG. 2c is a schematic diagram of a density field of a sag model, FIG. 3a is a schematic diagram of a source-side longitudinal pseudo normal stress field in the xx direction constructed according to the longitudinal velocity model given in FIG. 2a, the shear velocity model given in FIG. 2b and the density model given in FIG. 2c by using a mode decoupling forward continuation operator to forward extrapolate to the time of 0.8s, FIG. 3b is a schematic diagram of a seismic source end shear pseudo normal stress field along the xx direction, which is constructed according to the longitudinal wave velocity model given in FIG. 2a, the shear wave velocity model given in FIG. 2b and the density model given in FIG. 2c by forward extrapolation of a mode decoupling forward prolongation operator to 0.8 s. The separation of the elastic stress field is evident from fig. 3a, 3 b.

S2: and constructing a seismic source end energy vector field by using the seismic source end vector field information.

In an embodiment of the present specification, the constructing a seismic source energy vector field by using the seismic source vector field information may include: and constructing a longitudinal wave energy vector field and a transverse wave energy vector field of the seismic source end by using the elastic energy density as a field function and using the vector field information of the seismic source end.

For example, in some implementation scenarios, the elastic energy density is used as a field function, and the obtained longitudinal wave particle vibration velocity vector field, the pseudo stress vector field, the transverse wave particle vibration velocity vector field, and the pseudo stress vector field of the seismic source at each moment can be used to construct the longitudinal wave energy vector field and the transverse wave energy vector field of the seismic source according to the following formulas:

wherein the content of the first and second substances,

and

respectively representing a longitudinal wave energy vector field at the seismic source end and a transverse wave energy vector field at the seismic source end,

and

respectively representing the longitudinal wave pseudo body stress of the seismic source end and the transverse wave pseudo body stress of the seismic source end,

in one embodiment of the present description, after the seismic source end energy vector field is constructed, the longitudinal wave energy vector field and the transverse wave energy vector field at each time point of the seismic source end can be saved.

S4: and (4) decoupling a reverse continuation operator based on the constructed mode to obtain vector field information of the detection end.

The vector field information may include a particle vibration velocity vector field and a pseudo stress vector field of a longitudinal wave, and a particle vibration velocity vector field and a pseudo stress vector field of a transverse wave. The pickup-side vector field information may include a particle vibration velocity vector field and a pseudo stress vector field of a pickup-side longitudinal wave, and a particle vibration velocity vector field and a pseudo stress vector field of a pickup-side transverse wave.

In an embodiment of the present specification, a given medium model and a multi-component seismic record can be used, and a constructed mode decoupling reverse continuation operator is used to realize elastic wave field reverse extrapolation, so as to obtain a longitudinal wave particle vibration velocity vector field-pseudo stress vector field and a transverse wave particle vibration velocity vector field-pseudo stress vector field decoupled at a detection end.

In an embodiment of the present specification, constructing the pattern decoupling reverse prolongation operator may include: constructing a mode decoupling reverse continuation equation by using a preset medium model and taking multi-component seismic records as boundary conditions; discretizing the mode decoupling reverse continuation equation to obtain a mode decoupling reverse continuation operator.

Fig. 4a and 4b show given multi-component seismic records, wherein fig. 4a is a schematic diagram of a 2.7s multi-component seismic record acquired from the surface along the x-direction, and fig. 4b is a schematic diagram of a 2.7s multi-component seismic record acquired from the surface along the z-direction. In some implementation scenarios, the elastic wave field backward time extension may be performed by using the multi-component seismic records predefined in fig. 4a and 4b as boundary conditions according to a given medium model, that is, the following mode decoupling backward time extension equation is used for performing the backward time extension:

wherein r ═ r (r)_x,r_y,r_z)^TRepresenting a predetermined multi-component seismic recording vector field, r_xRepresenting a component in the x-direction, r_yRepresenting a component in the y-direction, r_zRepresenting a component in the z direction, z₀Representing the depth of burial at the time of seismic record acquisition. It should be noted that other letters can refer to the above explanation of the letter in formula (2), and the description thereof is omitted.

In some implementation scenarios, discretizing the mode decoupling reverse continuation equation can obtain the following mode decoupling reverse continuation operators:

where the superscript "R" denotes the detection end wavefield function, and, accordingly,

and

respectively showing the stress field of longitudinal wave body and the stress vector field of transverse wave at the detection end,

and

respectively showing the positive strain of the detection-end transverse wave in the x direction, the y direction and the z direction,

and

respectively represents the shear stress of the detection end transverse wave in the yz direction, in the xz direction and in the xy direction,

represents a longitudinal wave pseudo stress vector field at the detection end,

and

respectively showing the detection end longitudinal wave pseudo-normal stress along the x direction, along the y direction and along the z direction,

and

respectively represents the detection end longitudinal wave pseudo-shear stress along the yz direction, the xz direction and the xy direction, the subscript mark 'qp' represents the longitudinal wave pseudo-stress mode,

represents a pseudo stress vector field of a transverse wave at the detection end,

and

respectively represents the wave detection end transverse wave pseudo-normal stress along the x direction, the y direction and the z direction,

and

respectively represent the detection end shear wave pseudo-shear stress in the yz direction, in the xz direction and in the xy direction, the subscript "qs" represents the shear wave pseudo-stress mode,

representing the vibration velocity vector field of longitudinal wave mass point at the detection end,

and

respectively representing components of the detection end longitudinal wave particle vibration velocity vector field along the x direction, the y direction and the z direction,

representing the vibration velocity vector field of transverse wave mass point at the detection end,

and

respectively representing the components of the vibration velocity of the transverse wave mass point at the detection end along the x direction, the y direction and the z direction,

representing the vibration velocity vector field of the elastic wave mass point at the detection end,

and

representing components of the vibration velocity of the elastic wave mass point at the detection end along the x direction, the y direction and the z direction respectively, eta represents a boundary absorption coefficient, delta T represents a time sampling interval, N delta T represents a whole time node corresponding to a discrete point N, (N +1/2) delta T represents a half time node corresponding to a discrete point (N +1/2), (N-1/2) delta T represents a half time node corresponding to a discrete point (N-1/2) delta T, N is 1,2₀N Δ T is gradually reduced to 0, T₀N Δ t represents the total seismic recording reception time,

and

In some embodiments, after obtaining the mode decoupling reverse continuation operator, the detection end vector field information may be obtained based on the constructed mode decoupling reverse continuation operator. For example, in some implementation scenarios, the constructed mode decoupling reverse continuation operator can be used to obtain a decoupled pickup end compressional particle vibration velocity vector field, a pseudo stress vector field, and shear particle vibration velocity vector field, and a pseudo stress vector field.

S6: and constructing a detection end energy vector field by using the detection end vector field information.

In an embodiment of the present specification, the constructing a detection-end energy vector field by using the detection-end vector field information may include: and constructing a longitudinal wave energy vector field and a transverse wave energy vector field of the seismic source end by using the detection end vector field information by taking the elastic energy density as a field function.

For example, in some implementation scenarios, the elastic energy density is used as a field function, and the obtained longitudinal wave particle vibration velocity vector field, pseudo stress vector field, and transverse wave particle vibration velocity vector field, pseudo stress vector field of the detector at each moment can be used to construct the longitudinal wave energy vector field and transverse wave energy vector field of the detector according to the following formulas according to the principle of forward direction consistency of the base vectors:

wherein the content of the first and second substances,

and

respectively representing a longitudinal wave energy vector field and a transverse wave energy vector field of the detection end,

and

respectively represents the stress of a longitudinal wave pseudo body at the detection end and the stress of a transverse wave pseudo body at the detection end, the sign "-" is a negative sign and represents the direction opposite to the time change,

it should be noted that the principle of consistency in the positive base vector direction can be understood that the positive base vector directions of the source-end energy vector field and the detection-end energy vector field are consistent.

In one embodiment of the present disclosure, after the detection end energy vector field is constructed, the detection end longitudinal and transverse wave energy vector fields at each time can be stored.

S8: and imaging the seismic source end energy vector field and the detection end energy vector field according to the constructed mode decoupling elastic energy cross-correlation imaging condition to obtain a single shot mode decoupling energy imaging result.

The energy imaging results of the single shot mode decoupling can comprise PP wave energy imaging results, PS wave energy imaging results, SP wave energy imaging results and SS wave energy imaging results.

In an embodiment of the present specification, the obtained source end longitudinal wave and transverse wave energy vector fields and the detection end longitudinal wave and transverse wave energy vector fields may be used to respectively image according to the constructed mode decoupling elastic energy cross-correlation imaging conditions, so as to obtain a mode decoupling energy imaging result. For example, in some implementation scenarios, for each continuation time, energy inner product operations may be performed on the obtained decoupled longitudinal and transverse wave energy vector fields at the source end of the corresponding time and the obtained decoupled longitudinal and transverse wave energy vector fields at the detection end of the corresponding time, and time superposition may be performed on the energy inner product results at all the continuation times, that is, imaging is performed under the constructed mode decoupling elastic energy cross-correlation imaging conditions, so as to obtain the mode decoupling energy imaging result.

In one embodiment of the present disclosure, the constructed mode-decoupled elastic energy cross-correlation imaging condition may be expressed as:

wherein, I_E,PPRepresenting PP wave energy imaging results, I_E,PSShows the PS wave energy imaging result, I_E,SPRepresenting SP wave energy imaging results, I_E,SSRepresenting SS wave energy imaging result, t representing time, e representing Young modulus of background medium, g representing Poisson's ratio of the background medium, p representing density of the background medium, and symbol<*，*>_EThe sign of the operation of the inner product of energy is expressed,

and

and

and

and

respectively showing the detecting end longitudinal wave pseudo body stress and the detecting end transverse wave pseudo body stress,

is a transverse wave particle vibration velocity vector field at the detection end,

and

and

respectively represents the pseudo shear stress of the source end longitudinal wave along the yz direction, the xz direction and the xy direction,

and

and

respectively represents the pseudo shear stress of the source end transverse wave along the yz direction, the xz direction and the xy direction,

and

and

respectively represents the pseudo-shear stress of the detection end longitudinal wave along the yz direction, the xz direction and the xy direction,

and

and

the method is characterized in that the method respectively represents the detection end transverse wave pseudo shear stress along the yz direction, the xz direction and the xy direction, the subscript symbol 'p' represents a longitudinal wave mode, the subscript symbol 'S' represents a transverse wave mode, the superscript 'S' represents a seismic source end wave field function, the superscript 'R' represents a detection end wave field function, the subscript symbol 'qs' represents a transverse wave pseudo stress mode, and the subscript symbol 'qp' represents a longitudinal wave pseudo stress mode. Wherein the energy inner product can be expressed as

S10: and performing multi-shot superposition on the energy imaging result decoupled in the single shot mode to obtain an elastic energy reverse time migration imaging result.

In one embodiment of the present description, multiple shot superposition may be performed on the obtained decoupled PP wave energy, PS wave energy, SP wave energy, and SS wave energy imaging results of a single shot to obtain the final imaging result.

The final imaging results are shown in fig. 5a, 5b, 5c and 5d, wherein fig. 5a is a schematic diagram of a multi-shot stacked final PP wave energy imaging profile, fig. 5b is a schematic diagram of a final stacked PS wave energy imaging profile, fig. 5c is a schematic diagram of a final stacked SP wave energy imaging profile, and fig. 5d is a schematic diagram of a final stacked SS wave energy imaging profile.

According to the elastic wave field energy reverse time migration imaging method, the mode decoupling forward continuation operator and the mode decoupling reverse continuation operator are constructed on the basis of decoupling continuation, so that mode decoupling of an elastic wave field particle vibration velocity field and a stress field can be realized while forward continuation and reverse continuation of an elastic wave field are realized, and the energy coupling problem in an elastic wave field energy reverse time migration imaging result can be effectively reduced. The constructed mode decoupling energy cross-correlation imaging condition is utilized to obtain the energy imaging result of the pure wave mode and the converted wave energy imaging result from the decoupled elastic energy field, so that the difficulty of interpretation work by utilizing the energy imaging result of the pure wave mode and the converted wave energy imaging result in the follow-up process can be reduced, and the imaging efficiency is improved.

In the present specification, each embodiment of the method is described in a progressive manner, and the same and similar parts in each embodiment may be joined together, and each embodiment focuses on the differences from the other embodiments. Relevant points can be obtained by referring to part of the description of the embodiment of the method.

Based on the elastic energy reverse time migration imaging method, one or more embodiments of the present disclosure further provide an elastic energy reverse time migration imaging apparatus. The apparatus may include systems (including distributed systems), software (applications), modules, components, servers, clients, etc. that use the methods described in the embodiments of the present specification in conjunction with any necessary apparatus to implement the hardware. Based on the same innovative conception, embodiments of the present specification provide an apparatus as described in the following embodiments. Since the implementation scheme of the apparatus for solving the problem is similar to that of the method, the specific implementation of the apparatus in the embodiment of the present specification may refer to the implementation of the foregoing method, and repeated details are not repeated. As used hereinafter, the term "unit" or "module" may be a combination of software and/or hardware that implements a predetermined function. Although the means described in the embodiments below are preferably implemented in software, an implementation in hardware, or a combination of software and hardware is also possible and contemplated.

Specifically, fig. 6 is a schematic block diagram of an embodiment of an elastic energy reverse time migration imaging device provided in the present specification, and as shown in fig. 6, the elastic energy reverse time migration imaging device provided in the present specification may include: the seismic source end vector field information acquisition module 120, the seismic source end energy vector field construction module 122, the detection end vector field information acquisition module 124, the detection end energy vector field construction module 126, the single shot energy imaging module 128 and the elastic energy reverse time migration imaging module 130.

The seismic source-side vector field information obtaining module 120 may be configured to decouple the forward continuation operator based on the constructed mode to obtain seismic source-side vector field information, where the vector field information includes a particle vibration velocity vector field of a longitudinal wave, a pseudo stress vector field, and a particle vibration velocity vector field and a pseudo stress vector field of a transverse wave;

a seismic source end energy vector field construction module 122, configured to construct a seismic source end energy vector field by using the seismic source end vector field information;

the detection end vector field information obtaining module 124 may be configured to decouple the reverse continuation operator based on the constructed mode to obtain detection end vector field information;

a detection end energy vector field constructing module 126, configured to construct a detection end energy vector field by using the detection end vector field information;

the single-shot energy imaging module 128 may be configured to image the seismic source end energy vector field and the detection end energy vector field according to a constructed mode decoupling elastic energy cross-correlation imaging condition, so as to obtain a single-shot mode decoupling energy imaging result;

the elastic energy reverse time migration imaging module 130 may be configured to perform multi-shot superposition on the single-shot mode decoupled energy imaging result to obtain an elastic energy reverse time migration imaging result.

According to the elastic energy reverse time migration imaging device provided by the specification, the mode decoupling forward continuation operator and the mode decoupling reverse continuation operator are constructed on the basis of decoupling continuation, so that the mode decoupling of an elastic wave field particle vibration velocity field and a stress field can be realized while the forward continuation and the reverse continuation of an elastic wave field are realized, and the energy coupling problem in an elastic wave field energy reverse time migration imaging result can be effectively reduced. The constructed mode decoupling energy cross-correlation imaging condition is utilized to obtain the energy imaging result of the pure wave mode and the converted wave energy imaging result from the decoupled elastic energy field, so that the difficulty of interpretation work by utilizing the energy imaging result of the pure wave mode and the converted wave energy imaging result in the follow-up process can be reduced, and the imaging efficiency is improved.

It should be noted that the above-mentioned description of the apparatus according to the method embodiment may also include other embodiments, and specific implementation manners may refer to the description of the related method embodiment, which is not described herein again.

The present specification also provides embodiments of a resilient energy reverse time migration imaging device comprising a processor and a memory for storing processor-executable instructions which, when executed by the processor, implement steps comprising:

It should be noted that the above description of the apparatus according to the method embodiment may also include other embodiments. The specific implementation manner may refer to the description of the related method embodiment, and is not described in detail herein.

The present specification also provides embodiments of a resilient energy reverse time migration imaging system, comprising at least one processor and a memory storing computer-executable instructions, which when executed by the processor, implement the steps of the method described in any one or more of the above embodiments, for example, comprising: based on the constructed mode decoupling forward continuation operator, obtaining vector field information of a seismic source end, wherein the vector field information comprises a particle vibration velocity vector field and a pseudo stress vector field of longitudinal waves, and a particle vibration velocity vector field and a pseudo stress vector field of transverse waves; constructing a seismic source end energy vector field by using the seismic source end vector field information; based on the constructed mode decoupling reverse continuation operator, obtaining vector field information of a detection end; constructing a detection end energy vector field by using the detection end vector field information; imaging the seismic source end energy vector field and the detection end energy vector field according to the constructed mode decoupling elastic energy cross-correlation imaging condition to obtain a single shot mode decoupling energy imaging result; and performing multi-shot superposition on the energy imaging result decoupled in the single shot mode to obtain an elastic energy reverse time migration imaging result. The system may be a single server, or may include a server cluster, a system (including a distributed system), software (applications), an actual operating device, a logic gate device, a quantum computer, etc. using one or more of the methods or one or more of the example devices of the present specification, in combination with a terminal device implementing hardware as necessary.

The method embodiments provided in the present specification may be executed in a mobile terminal, a computer terminal, a server or a similar computing device. Taking an example of the application on a server, fig. 7 is a hardware block diagram of an embodiment of a resilient energy reverse time migration imaging server provided in this specification, where the server may be a resilient energy reverse time migration imaging apparatus or a resilient energy reverse time migration imaging system in the above embodiment. As shown in fig. 7, the server 10 may include one or more (only one shown) processors 100 (the processors 100 may include, but are not limited to, a processing device such as a microprocessor MCU or a programmable logic device FPGA, etc.), a memory 200 for storing data, and a transmission module 300 for communication functions. It will be understood by those skilled in the art that the structure shown in fig. 7 is only an illustration and is not intended to limit the structure of the electronic device. For example, the server 10 may also include more or fewer components than shown in FIG. 7, and may also include other processing hardware, such as a database or multi-level cache, a GPU, or have a different configuration than shown in FIG. 7, for example.

The memory 200 may be used to store software programs and modules of application software, such as program instructions/modules corresponding to the elastic energy reverse time shift imaging method in the embodiment of the present specification, and the processor 100 executes various functional applications and data processing by executing the software programs and modules stored in the memory 200. Memory 200 may include high speed random access memory and may also include non-volatile memory, such as one or more magnetic storage devices, flash memory, or other non-volatile solid-state memory. In some examples, memory 200 may further include memory located remotely from processor 100, which may be connected to a computer terminal through a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.

The transmission module 300 is used for receiving or transmitting data via a network. Specific examples of the network described above may include a wireless network provided by a communication provider of the computer terminal. In one example, the transmission module 300 includes a Network adapter (NIC) that can be connected to other Network devices through a base station so as to communicate with the internet. In one example, the transmission module 300 may be a Radio Frequency (RF) module, which is used for communicating with the internet in a wireless manner.

The foregoing description has been directed to specific embodiments of this disclosure. Other embodiments are within the scope of the following claims. In some cases, the actions or steps recited in the claims may be performed in a different order than in the embodiments and still achieve desirable results. In addition, the processes depicted in the accompanying figures do not necessarily require the particular order shown, or sequential order, to achieve desirable results. In some embodiments, multitasking and parallel processing may also be possible or may be advantageous.

The method or apparatus provided by the present specification and described in the foregoing embodiments may implement service logic through a computer program and record the service logic on a storage medium, where the storage medium may be read and executed by a computer, so as to implement the effect of the solution described in the embodiments of the present specification.

The storage medium may include a physical device for storing information, and typically, the information is digitized and then stored using an electrical, magnetic, or optical media. The storage medium may include: devices that store information using electrical energy, such as various types of memory, e.g., RAM, ROM, etc.; devices that store information using magnetic energy, such as hard disks, floppy disks, tapes, core memories, bubble memories, and usb disks; devices that store information optically, such as CDs or DVDs. Of course, there are other ways of storing media that can be read, such as quantum memory, graphene memory, and so forth.

The above embodiments of the elastic energy reverse time migration imaging method or apparatus provided in this specification may be implemented in a computer by a processor executing corresponding program instructions, for example, implemented in a PC using a c + + language of a windows operating system, implemented in a linux system, or implemented in an intelligent terminal using, for example, android, iOS system programming languages, implemented in processing logic based on a quantum computer, and the like.

It should be noted that descriptions of the apparatus, the computer storage medium, and the system described above according to the related method embodiments may also include other embodiments, and specific implementations may refer to descriptions of corresponding method embodiments, which are not described in detail herein.

The embodiments in the present application are described in a progressive manner, and the same and similar parts among the embodiments can be referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the hardware + program class embodiment, since it is substantially similar to the method embodiment, the description is simple, and the relevant points can be referred to the partial description of the method embodiment.

The embodiments of this specification are not limited to what must be in compliance with industry communication standards, standard computer data processing and data storage rules, or the description of one or more embodiments of this specification. Certain industry standards, or implementations modified slightly from those described using custom modes or examples, may also achieve the same, equivalent, or similar, or other, contemplated implementations of the above-described examples. The embodiments using the modified or transformed data acquisition, storage, judgment, processing and the like can still fall within the scope of the alternative embodiments of the embodiments in this specification.

In the 90 s of the 20 th century, improvements in a technology could clearly distinguish between improvements in hardware (e.g., improvements in circuit structures such as diodes, transistors, switches, etc.) and improvements in software (improvements in process flow). However, as technology advances, many of today's process flow improvements have been seen as direct improvements in hardware circuit architecture. Designers almost always obtain the corresponding hardware circuit structure by programming an improved method flow into the hardware circuit. Thus, it cannot be said that an improvement in the process flow cannot be realized by hardware physical modules. For example, a Programmable Logic Device (PLD), such as a Field Programmable Gate Array (FPGA), is an integrated circuit whose Logic functions are determined by programming the Device by a user. A digital system is "integrated" on a PLD by the designer's own programming without requiring the chip manufacturer to design and fabricate application-specific integrated circuit chips. Furthermore, nowadays, instead of manually making an Integrated Circuit chip, such Programming is often implemented by "logic compiler" software, which is similar to a software compiler used in program development and writing, but the original code before compiling is also written by a specific Programming Language, which is called Hardware Description Language (HDL), and HDL is not only one but many, such as abel (advanced Boolean Expression Language), ahdl (alternate Hardware Description Language), traffic, pl (core universal Programming Language), HDCal (jhdware Description Language), lang, Lola, HDL, laspam, hardward Description Language (vhr Description Language), vhal (Hardware Description Language), and vhigh-Language, which are currently used in most common. It will also be apparent to those skilled in the art that hardware circuitry that implements the logical method flows can be readily obtained by merely slightly programming the method flows into an integrated circuit using the hardware description languages described above.

The controller may be implemented in any suitable manner, for example, the controller may take the form of, for example, a microprocessor or processor and a computer-readable medium storing computer-readable program code (e.g., software or firmware) executable by the (micro) processor, logic gates, switches, an Application Specific Integrated Circuit (ASIC), a programmable logic controller, and an embedded microcontroller, examples of which include, but are not limited to, the following microcontrollers: ARC 625D, Atmel AT91SAM, Microchip PIC18F26K20, and Silicone Labs C8051F320, the memory controller may also be implemented as part of the control logic for the memory. Those skilled in the art will also appreciate that, in addition to implementing the controller as pure computer readable program code, the same functionality can be implemented by logically programming method steps such that the controller is in the form of logic gates, switches, application specific integrated circuits, programmable logic controllers, embedded microcontrollers and the like. Such a controller may thus be considered a hardware component, and the means included therein for performing the various functions may also be considered as a structure within the hardware component. Or even means for performing the functions may be regarded as being both a software module for performing the method and a structure within a hardware component.

The systems, devices, modules or units illustrated in the above embodiments may be implemented by a computer chip or an entity, or by a product with certain functions. One typical implementation device is a computer. In particular, the computer may be, for example, a personal computer, a laptop computer, a vehicle-mounted human-computer interaction device, a cellular telephone, a camera phone, a smart phone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or a combination of any of these devices.

Although one or more embodiments of the present description provide method operational steps as described in the embodiments or flowcharts, more or fewer operational steps may be included based on conventional or non-inventive approaches. The order of steps recited in the embodiments is merely one manner of performing the steps in a multitude of orders and does not represent the only order of execution. When an actual apparatus or end product executes, it may execute sequentially or in parallel (e.g., parallel processors or multi-threaded environments, or even distributed data processing environments) according to the method shown in the embodiment or the figures. The terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, the presence of additional identical or equivalent elements in a process, method, article, or apparatus that comprises the recited elements is not excluded. The terms first, second, etc. are used to denote names, but not any particular order.

For convenience of description, the above devices are described as being divided into various modules by functions, and are described separately. Of course, when implementing one or more of the present description, the functions of each module may be implemented in one or more software and/or hardware, or a module implementing the same function may be implemented by a combination of multiple sub-modules or sub-units, etc. The above-described embodiments of the apparatus are merely illustrative, and for example, the division of the units is only one logical division, and other divisions may be realized in practice, for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may be in an electrical, mechanical or other form.

The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

In a typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.

The memory may include forms of volatile memory in a computer readable medium, Random Access Memory (RAM) and/or non-volatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). Memory is an example of a computer-readable medium.

Computer-readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), Static Random Access Memory (SRAM), Dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), Read Only Memory (ROM), Electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), Digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage, graphene storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, a computer readable medium does not include a transitory computer readable medium such as a modulated data signal and a carrier wave.

As will be appreciated by one skilled in the art, one or more embodiments of the present description may be provided as a method, system, or computer program product. Accordingly, one or more embodiments of the present description may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, one or more embodiments of the present description may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.

The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the system embodiment, since it is substantially similar to the method embodiment, the description is simple, and for the relevant points, reference may be made to the partial description of the method embodiment. In the description of the specification, reference to the description of the term "one embodiment," "some embodiments," "an example," "a specific example," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the specification. In this specification, the schematic representations of the terms used above are not necessarily intended to refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, various embodiments or examples and features of different embodiments or examples described in this specification can be combined and combined by one skilled in the art without contradiction.

The above description is merely exemplary of one or more embodiments of the present disclosure and is not intended to limit the scope of one or more embodiments of the present disclosure. Various modifications and alterations to one or more embodiments described herein will be apparent to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the scope of the claims.

Claims

1. An elastic energy reverse time migration imaging method, comprising:

imaging the seismic source end energy vector field and the detection end energy vector field according to the constructed mode decoupling elastic energy cross-correlation imaging condition to obtain a single shot mode decoupling energy imaging result; wherein the constructed mode-decoupled elastic energy cross-correlation imaging condition is expressed as:

wherein, I_E,PPRepresenting PP wave energy imaging results, I_E,PSShows the PS wave energy imaging result, I_E,SPRepresenting SP wave energy imaging results, I_E,SSRepresents SS wave energy imaging results, t represents time, e represents Young modulus of a background medium, g represents Poisson's ratio of the background medium, p represents density of the background medium,

and

and

and

and

and

and

and

and

and

and

and

and

respectively representing the pseudo shear stress of the detection end transverse wave along yz direction, xz direction and xy direction, the subscript symbol 'p' represents a longitudinal wave mode, the subscript symbol 'S' represents a transverse wave mode, the superscript 'S' represents a seismic source end wave field function, the superscript 'R' represents a detection end wave field function, the subscript symbol 'qs' represents a transverse wave pseudo stress mode, the subscript symbol 'qp' represents a longitudinal wave pseudo stress mode, the subscript symbol 'E' represents energy, T₀Representing a total seismic record reception duration;

2. The method of claim 1, wherein constructing the pattern decoupling forward prolongation operator comprises:

loading a seismic source wavelet by utilizing a preset medium model;

loading the seismic source wavelet to the longitudinal wave body stress of an elastic wave field, and constructing a mode decoupling forward continuation equation;

discretizing the mode decoupling forward continuation equation to obtain the following mode decoupling forward continuation operators:

wherein the content of the first and second substances,

and

respectively representing the stress field of longitudinal wave body and the stress vector field of transverse wave at the seismic source end, F_pA representation of the source wavelet is shown,

representing the vibration velocity vector field, v, of the seismic source-side transverse wave particle^SRepresenting the vibration velocity vector field of the elastic wave particle at the seismic source end,

and

respectively representing a longitudinal wave pseudo stress vector field and a transverse wave pseudo stress vector field at a seismic source end, eta represents a boundary absorption coefficient, delta t represents a time sampling interval, N delta t represents a whole time point corresponding to a discrete point N, (N +1/2) delta t represents a half-time node corresponding to a discrete point (N +1/2), N is 1,2Stiffness matrix of background Medium, C_sRepresenting the stiffness matrix in the background dielectric shear mode,

and

respectively representing a high-order staggered grid finite difference matrix operator, lambda and mu representing the Lame coefficient of the background medium, and rho representing the density of the background medium.

3. The method of claim 1, wherein said using said source end vector field information to construct a source end energy vector field comprises:

and constructing a longitudinal wave energy vector field and a transverse wave energy vector field of the seismic source end according to the following formula by taking the elastic energy density as a field function:

wherein the content of the first and second substances,

and

and

and

respectively representing a longitudinal wave pseudo stress vector field and a transverse wave pseudo stress vector field at the seismic source end.

4. The method of claim 1, wherein constructing the pattern decoupling reverse prolongation operator comprises:

constructing a mode decoupling reverse continuation equation by using a preset medium model and taking multi-component seismic records as boundary conditions;

discretizing the mode decoupling reverse continuation equation to obtain the following mode decoupling reverse continuation operators:

wherein the content of the first and second substances,

and

v represents the vibration velocity vector field of transverse wave particle at the detection end^RRepresenting the vibration velocity vector field of the elastic wave mass point at the detection end,

and

respectively representing a longitudinal wave pseudo-stress vector field and a transverse wave pseudo-stress vector field at a wave detection end, r representing a preset multi-component seismic record, eta representing a boundary absorption coefficient, delta t representing a time sampling interval, N delta t representing a whole time point corresponding to a discrete point N, (N +1/2) delta t representing a half-time node corresponding to the discrete point (N +1/2), (N-1/2) delta t representing a half-time node corresponding to the discrete point (N-1/2) delta t, N being 1,2, N, N representing discrete points corresponding to the total receiving time of the seismic record, C representing a stiffness matrix of a background medium, C representing a stiffness matrix of the background medium, C representing a total receiving time of the seismic record, and_srepresenting the stiffness matrix in the background dielectric shear mode,

and

5. The method of claim 1, wherein said using said detector end vector field information to construct a detector end energy vector field comprises:

and constructing a detection end longitudinal wave energy vector field and a detection end transverse wave energy vector field by taking the elastic energy density as a field function according to the following formula:

wherein the content of the first and second substances,

and

and

and

respectively showing a longitudinal wave pseudo stress vector field and a transverse wave pseudo stress vector field of the detection end.

6. An elastic energy reverse time migration imaging apparatus, comprising:

the single-shot energy imaging module is used for imaging the seismic source end energy vector field and the detection end energy vector field according to the constructed mode decoupling elastic energy cross-correlation imaging condition to obtain a single-shot mode decoupling energy imaging result; wherein the constructed mode-decoupled elastic energy cross-correlation imaging condition is expressed as:

wherein，I_E,PPRepresenting PP wave energy imaging results, I_E,PSShows the PS wave energy imaging result, I_E,SPRepresenting SP wave energy imaging results, I_E,SSRepresents SS wave energy imaging results, t represents time, e represents Young modulus of a background medium, g represents Poisson's ratio of the background medium, p represents density of the background medium,

and

and

and

and

and

and

and

and

and

and

and

and

7. A resilient energy reverse time migration imaging device comprising a processor and a memory for storing processor-executable instructions that when executed by the processor implement steps comprising:

and

and

and

and

and

and

and

and

and

and

and

and

respectively representing the pseudo shear stress of the detection end transverse wave along yz direction, xz direction and xy direction, the subscript symbol 'p' represents a longitudinal wave mode, the subscript symbol 'S' represents a transverse wave mode, the superscript 'S' represents a seismic source end wave field function, the superscript 'R' represents a detection end wave field function, the subscript symbol 'qs' represents a transverse wave pseudo stress mode, the subscript symbol 'qp' represents a longitudinal wave pseudo stress mode, and the subscript symbol 'E' tableEnergy indicating, T₀Representing a total seismic record reception duration;

8. A resilient energy reverse time migration imaging system comprising at least one processor and a memory storing computer executable instructions which when executed by the processor implement the steps of the method of any one of claims 1 to 5.