CN108563802B - Method for improving numerical simulation precision of seismic converted waves - Google Patents
Method for improving numerical simulation precision of seismic converted waves Download PDFInfo
- Publication number
- CN108563802B CN108563802B CN201711483283.XA CN201711483283A CN108563802B CN 108563802 B CN108563802 B CN 108563802B CN 201711483283 A CN201711483283 A CN 201711483283A CN 108563802 B CN108563802 B CN 108563802B
- Authority
- CN
- China
- Prior art keywords
- wave
- equation
- formula
- field
- elastic 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Abstract
The invention discloses a method for improving the numerical simulation precision of seismic converted waves, which comprises the following steps: setting an original elastic wave fluctuation equation; introducing transverse wave components Ss and Sp to obtain a modified elastic wave fluctuation equation; and (3) solving the corrected elastic wave equation after equation equivalence is proved through mathematical formula derivation to obtain an Ss wave field and an Sp wave field, namely a converted wave field and a longitudinal wave field. The invention obtains P, S wave field completely separated, the wave field is clean, and there is no interference between longitudinal wave and transverse wave. The amplitude and the phase of the P, S wave after separation are kept unchanged, which is beneficial to further research on P, S propagation characteristics.
Description
Technical Field
The invention belongs to the technical field of seismic numerical simulation, and particularly relates to a method for improving the numerical simulation precision of seismic converted waves.
Background
The traditional elastic wave field obtained based on the numerical simulation of the elastic wave fluctuation equation contains longitudinal waves (P waves) and transverse waves (S waves), the two wave fields are separated before the P waves or the S waves are processed, and the process is difficult to ensure that the wave fields of the P waves and the S waves are not distorted and can also hardly ensure that the two wave fields can be completely separated. In recent years, a seabed node technology characterized by four-component reception is gradually developed in the field of earthquake, especially, the application of converted waves is more and more concerned by geophysicists, and the converted waves have very important significance for searching fractured hydrocarbon reservoirs, estimating lithology and the like.
At present, the elastic wave simulation can only obtain a mixed elastic wave field, including P waves and S waves, and the separation of longitudinal waves and transverse waves by a traditional gradient rotation operator is difficult to ensure that the wave field after separation has no amplitude and phase distortion.
Disclosure of Invention
In view of the above, the present invention provides a method for improving the accuracy of numerical simulation of seismic converted waves.
In order to solve the technical problem, the invention discloses a method for improving the numerical simulation precision of seismic converted waves, which comprises the following steps:
step 1: setting an original elastic wave fluctuation equation which is composed of a formula (1-1), a formula (1-2) and a formula (1-3),
wherein x, y and z represent three directions of space, u, v and w represent components of the elastic wave field in the three directions of x, y and z, and t represents propagation time; vp represents a longitudinal wave propagation velocity, and Vs represents a transverse wave propagation velocity;
step 2: introducing transverse wave components Ss and Sp into the original elastic wave fluctuation equation in the step 1 to obtain a modified elastic wave fluctuation equation;
introduction of variables Sp, Ss
Sp=(UP、VP、WP)(2)
Ss=(Us,Vs,Ws)(3)
Wherein, Up, Vp, Wp represent the components of P wave in X, Y, Z directions, Us, Vs, Ws represent the components of S wave in X, Y, Z directions;
substituting the formulas (2) and (3) into the original elastic wave fluctuation equation in the step one to obtain a corrected elastic wave fluctuation equation; the corrected elastic wave fluctuation equation consists of a formula (4-1) and a formula (4-9); the formula (4-1) -the formula (4-9) forms an equation set (4);
U=UP+US(4-1)
V=VP+VS(4-2)
W=WP+WS(4-3)
and step 3: and (3) solving the corrected elastic wave equation in the step (2) after equation equivalence is proved through mathematical formula derivation to obtain expressions of the Ss wave field and the Sp wave field, and realizing the expressions of the Ss wave field and the Sp wave field, namely a converted wave field and a longitudinal wave field, through a program.
Further, the derivation of the mathematical formula in step 3 proves that the equation is equivalent to:
substituting the (4-4) th formula and the (4-7) th formula in the modified elastic wave fluctuation equation into the (4-1) th formula,
the above equation is formula (5), and formula (5) is fully equivalent to formula (1-1) in equation set (1);
similarly, an equation set obtained by substituting the (4-5) th equation and the (4-8) th equation in the modified elastic wave fluctuation equation into the (4-2) th equation is completely equivalent to the (1-2) th equation in the original elastic wave fluctuation equation: substituting the (4-6) th formula and the (4-9) th formula in the modified elastic wave fluctuation equation into the (4-2) th formula to obtain an equation set which is completely equivalent to the (1-1) th formula in the original elastic wave fluctuation equation, namely the modified elastic wave fluctuation equation is equivalent to the original elastic wave fluctuation equation;
then, it turns out that the introduced variables Sp and Ss can represent the compressional wave field and the shear wave field, respectively:
because the second derivative of the Sp rotation is 0, the first derivative of the Sp rotation is deduced to be a constant, the Sp rotation is a linear function or a constant, and the Sp rotation is obtained according to the condition of the initial zero solution of the wave equation, wherein the Sp rotation is 0, namely Sp is a non-rotation field, and Sp is a longitudinal wave; likewise, Ss is proven to be a shear wave.
Further, the step 3 of solving the corrected elastic wave equation to obtain the Ss wave field and the Sp wave field specifically includes: solving equations (4-1) - (4-9) by finite difference method to obtain the expressions of P wave and S wave, space 6 order:
wherein ui represents a P wave, uj represents an S wave, and DH represents a grid length; and realizing the expression of P wave and S wave by a program to obtain a converted wave shear wave field and a converted wave longitudinal wave field.
Compared with the prior art, the invention can obtain the following technical effects:
1) the invention obtains P, S wave field completely separated, the wave field is clean, and there is no interference between longitudinal wave and transverse wave.
2) The amplitude and the phase of the P, S wave after separation are kept unchanged, which is beneficial to further research on P, S propagation characteristics.
3) The invention accelerates based on the GPU and improves the simulation efficiency.
Of course, it is not necessary for any one product in which the invention is practiced to achieve all of the above-described technical effects simultaneously.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the invention and not to limit the invention. In the drawings:
FIG. 1 is a diagram of a longitudinal wave velocity model according to the present invention;
FIG. 2 is a diagram of a shear wave velocity model according to the present invention;
FIG. 3 is a snapshot of the longitudinal wave field at time 0.5s according to the present invention;
FIG. 4 is a fast graph of the converted wave field at 0.5s according to the present invention;
FIG. 5 is a P-wave single shot record diagram of the present invention;
fig. 6 is a diagram of an S-wave single shot record of the present invention.
Detailed Description
The following embodiments are described in detail with reference to the accompanying drawings, so that how to implement the technical features of the present invention to solve the technical problems and achieve the technical effects can be fully understood and implemented.
The idea of the invention is that: the invention provides an equivalent elastic wave equation, which takes P wave and S wave as variables, obtains a completely separated converted wave field by solving the wave field equation and simultaneously obtains a longitudinal wave field. The method can play an important role in the aspects of elastic wave theory research, converted wave field characteristic analysis, converted wave processing explanation and the like.
The theoretical basis of the converted wave numerical simulation is that through the replacement of variables, transverse wave components Ss and Sp are introduced into an original elastic wave fluctuation equation, and after the equation equivalence is proved through the derivation of a mathematical formula, a corrected elastic fluctuation equation is solved to obtain an Ss wave field and an Sp wave field, namely a converted wave field and a longitudinal wave field. The specific mathematical derivation process is as follows:
wherein x, y and z represent three directions of space, u, v and w represent components of the elastic wave field in the three directions of x, y and z, and t represents propagation time; vp represents a longitudinal wave propagation velocity, and Vs represents a transverse wave propagation velocity.
Formula (1) represents the original elastic wave fluctuation equation, including formula (1-1), formula (1-2), and formula (1-3);
setting a longitudinal wave velocity model, wherein the velocity in the velocity model represents Vp, and as shown in fig. 1, the velocity model is divided into 3 layers, from top to bottom, the first layer is seawater with a velocity of 1500m/s, the second layer is 2000m/s, and the third layer is 2500 m/s; a wave detection point is on the surface, a shot point is on the surface, and an actual seismic data acquisition mode is simulated;
setting a transverse wave velocity model, wherein the velocity in the velocity model represents Vs, and as shown in FIG. 2, the transverse wave velocity of each layer is obtained by converting longitudinal wave velocity through a geological empirical formula, the velocity model is divided into 3 layers, from top to bottom, the velocity of the first layer is 866m/s, the velocity of the second layer is 1155m/s, and the velocity of the third layer is 1443 m/s; the demodulator probe is on the surface, and the shot point is on the surface;
step 2, introducing transverse wave components Ss and Sp to obtain a modified elastic wave fluctuation equation which is an equation set (4):
SP=(UP,VP,WP)(2)
Ss=(Us,Vs,Ws)(3)
wherein, Up, Vp, Wp represent the components of P wave in x, y, z directions, respectively, and Us, Vs, Ws represent the components of S wave in x, y, z directions, respectively.
Bringing the formulas (2) and (3) into the formula (1) to obtain the following equation set (4), wherein the equation (4) includes the formula (4-1) to the formula (4-9):
U=UP+US(4-1)
V=VP+VS(4-2)
W=WP+WS(4-3)
wherein, Up, Vp, Wp represent the components of P wave in x, y, z directions, Us, Vs, Ws represent the components of S wave in x, y, z directions; x, y and z represent three directions of space, u, v and w represent components of the elastic wave field in the three directions of x, y and z, and t represents propagation time; vp represents a longitudinal wave propagation velocity, and Vs represents a transverse wave propagation velocity.
Step 3, solving the corrected elastic wave equation after equation equivalence is proved through mathematical formula derivation to obtain an Ss wave field and an Sp wave field, namely a converted wave field and a longitudinal wave field:
first, equation (1) and equation (4) prove to be equivalent.
Substituting the (4-4) th and (4-7) th equations in equation (4) into the (4-1) th equation, we obtain:
wherein x, y and z represent three directions of space, u, v and w represent components of the elastic wave field in the three directions of x, y and z, and t represents propagation time. Vp represents a longitudinal wave propagation velocity, and Vs represents a transverse wave propagation velocity. Up and Us respectively represent the longitudinal wave field separation and the transverse wave field component on the x component;
it is easy to find that the formula (5) is completely equivalent to the (1-1) th formula in the equation set (1). The equation set (1) herein refers to equation (1), which represents the original elastic wave fluctuation equation.
Similarly, the equation set obtained by substituting the (4-5) th and (4-8) th equations in the equation set (4) into the (4-2) th equation is completely equivalent to the (1-2) th equation in the equation set (1): the equation set obtained by substituting the (4-6) th expression and the (4-9) th expression in the equation set (4) into the (4-2) th expression is completely equivalent to the (1-1) th expression in the equation set (1), that is, the equation set (4) is equivalent to the equation (1).
Then, it turns out that the introduced variables Sp and Ss can represent the compressional wave field and the shear wave field, respectively:
wherein, i, j, k respectively represent coordinates of three-dimensional space grid points in x, y, z directions and represent rotation operators
Since the second derivative of the Sp rotation is 0, it can be deduced that the first derivative of the Sp rotation is constant, and the Sp rotation is a linear function or constant, which can be obtained according to the condition of the initial zero solution of the wave equation, the Sp rotation is 0, i.e., Sp is a non-rotation field, and Sp is a longitudinal wave. Likewise, Ss may also be demonstrated as a shear wave.
Example 1 takes a two-dimensional case as an example, and solves the problem by a finite difference method.
A method for improving the numerical simulation precision of seismic converted waves is carried out according to the following steps:
U=UP+US(4-1)
W=WP+WS(4-3)
step 2: solving the equation by a finite difference method to obtain expressions (6-order space) of P wave and S wave:
wherein i, j represents coordinates in x, y directions of spatial grid points involved in calculation, corresponding to i, j in the context, and a (n) represents a difference coefficient; ui denotes a P wave, uj denotes an S wave, and DH denotes a distance of the adjacent mesh involved in the calculation.
Step 3, realizing the discrete expression in the step 2 through a program to obtain a converted wave transverse wave field and a converted wave longitudinal wave field, and specifically comprising the following steps:
step 3.1, for the elastic wave equation in two-dimensional isotropic medium
Where x, y, and z represent three directions in space, λ and μ represent lame constants, τ zz, τ xx, and τ xz represent components of stress, and vx and vz represent components of velocity components.
Second order approximation of time
vx, vz denote the components of the velocity component, Δ t denotes the sampling rate, and m denotes the order taken after taylor expansion.
Subtracting the two equations to obtain an approximate difference equation of 2M order
vx represents the component of the velocity component and Δ t represents the sampling rate.
In the above n-pair equation, the equal-sign sides of the nth pair equation are simultaneously multiplied by the coefficient cn (n), and then subtracted, and all the equations are summed to obtain:
only a term is reserved at the right end of the formula, the coefficient is 1, and the other coefficients are 0, so that the following equation set is obtained
Substituting the above formula into (16) to obtain
Substituted into (15) to obtain
The same can be said for all discrete formats of stress-displacement equations:
wherein i, j, k respectively represent coordinates of three-dimensional space grid points in x, y, z directions, x, y, z represent three directions of space, λ μ represents a Lamei constant, τ zz, τ xx, τ xz represent components of stress, vx, vz represent components of velocity components, Δ t represents a sampling rate, N represents a difference order, Δ x, Δ z represent grid steps in the horizontal and vertical directions, and cn (N) is a coefficient.
Thus, a simulation result of the elastic wave was obtained.
And 3.2, extracting longitudinal and transverse waves from vx and vz components obtained by elastic wave simulation through formulas (7) and (8), wherein the derivation processes of the formulas (7) and (8) are as above.
Where i, j represents the coordinates of a certain grid point, a (n) represents the differential system ui represents the P wave, uj represents the S wave, and DH represents the grid length, i.e., Δ x, Δ z as mentioned in the above formula.
The converted wave extraction method in the current market obtains a longitudinal and transverse wave field by separating an elastic wave field by using a divergence rotation operator, the method hardly causes the amplitude and phase information of the wave field to be distorted, and is unfavorable for later data research, so that the invention provides a completely separated longitudinal and transverse wave data simulation method. Compared with the traditional divergence curl operator for separating longitudinal waves and transverse waves, the longitudinal waves and the transverse waves obtained by the method are completely separated from each other and do not interfere with each other, and the obtained longitudinal wave and transverse wave fields retain original amplitude and phase information, so that the method has great benefits on converted wave characteristic research, converted wave data processing and inversion method research.
Fig. 3 is a wave field snapshot of longitudinal waves, fig. 4 is a wave field snapshot of converted transverse waves, which is obtained from step 2 and step 3 of the embodiment, and through comparison of the two graphs, the longitudinal waves and the transverse waves can be found to be completely separated from each other and not to interfere with each other, thereby verifying the effectiveness of the invention.
The longitudinal wave record of fig. 5, fig. 6 is a shear wave record, which is obtained from step 2 and step 3 of the specific embodiment, fig. 3 and fig. 4 are by wave field snapshot comparison, fig. 5 and fig. 6 are by record comparison finally obtained, the P wave and the S wave are completely separated from each other, no converted wave is generated at the collar offset, the polarity of the converted wave is opposite on both sides of the shot point, which are the same as the actual phenomenon, and the effectiveness of the invention is proved again.
As can be seen from fig. 3 to 6, the two wave fields are completely independent wave fields, and there are no shear waves on the compressional wave recordings and no compressional waves on the compressional wave recordings. The single shot record is generally the final result of forward modeling, and from the single shot, the longitudinal and transverse waves are completely independent wave fields, so that the phenomenon that the transverse waves exist on the longitudinal wave field or the longitudinal waves exist on the transverse wave record due to incomplete separation of the longitudinal waves and the transverse waves in the traditional method can be avoided.
While the foregoing description shows and describes several preferred embodiments of the invention, it is to be understood, as noted above, that the invention is not limited to the forms disclosed herein, but is not to be construed as excluding other embodiments and is capable of use in various other combinations, modifications, and environments and is capable of changes within the scope of the inventive concept as expressed herein, commensurate with the above teachings, or the skill or knowledge of the relevant art. And that modifications and variations may be effected by those skilled in the art without departing from the spirit and scope of the invention as defined by the appended claims.
Claims (3)
1. A method for improving numerical simulation precision of seismic converted waves is characterized by comprising the following steps:
step 1: setting an original elastic wave fluctuation equation which is composed of a formula (1-1), a formula (1-2) and a formula (1-3),
wherein x, y and z represent three directions of space, u, v and w represent components of the elastic wave field in the three directions of x, y and z, and t represents propagation time; vp represents a longitudinal wave propagation velocity, and Vs represents a transverse wave propagation velocity;
step 2: introducing transverse wave components Ss and Sp into the original elastic wave fluctuation equation in the step 1 to obtain a modified elastic wave fluctuation equation;
introduction of variables Sp, Ss
Sp=(UP、VP、WP)(2)
Ss=(Us,Vs,Ws)(3)
Wherein, Up, Vp, Wp represent the components of P wave in X, Y, Z directions, Us, Vs, Ws represent the components of S wave in X, Y, Z directions;
substituting the formulas (2) and (3) into the original elastic wave fluctuation equation in the step one to obtain a corrected elastic wave fluctuation equation; the corrected elastic wave fluctuation equation consists of a formula (4-1) and a formula (4-9); the formula (4-1) -the formula (4-9) forms an equation set (4);
U=UP+US(4-1)
V=VP+VS(4-2)
W=WP+WS(4-3)
and step 3: and (3) solving the corrected elastic wave equation in the step (2) after equation equivalence is proved through mathematical formula derivation to obtain expressions of the Ss wave field and the Sp wave field, and realizing the expressions of the Ss wave field and the Sp wave field, namely a converted wave field and a longitudinal wave field, through a program.
2. The method for improving the numerical simulation accuracy of the seismic converted waves according to claim 1, wherein the equation equivalence proved by the derivation of the mathematical formula in the step 3 is as follows:
substituting the (4-4) th formula and the (4-7) th formula in the modified elastic wave fluctuation equation into the (4-1) th formula,
the above equation is formula (5), and formula (5) is fully equivalent to formula (1-1) in equation set (1);
similarly, an equation set obtained by substituting the (4-5) th equation and the (4-8) th equation in the modified elastic wave fluctuation equation into the (4-2) th equation is completely equivalent to the (1-2) th equation in the original elastic wave fluctuation equation: substituting the (4-6) th formula and the (4-9) th formula in the modified elastic wave fluctuation equation into the (4-2) th formula to obtain an equation set which is completely equivalent to the (1-1) th formula in the original elastic wave fluctuation equation, namely the modified elastic wave fluctuation equation is equivalent to the original elastic wave fluctuation equation;
then, it turns out that the introduced variables Sp and Ss can represent the compressional wave field and the shear wave field, respectively:
wherein: i, j and k respectively represent coordinates of three-dimensional space grid points in x, y and z directions and represent rotation operators; because the second derivative of the Sp rotation is 0, the first derivative of the Sp rotation is deduced to be a constant, the Sp rotation is a linear function or a constant, and the Sp rotation is obtained according to the condition of the initial zero solution of the wave equation, wherein the Sp rotation is 0, namely Sp is a non-rotation field, and Sp is a longitudinal wave; likewise, Ss is proven to be a shear wave.
3. The method for improving the accuracy of numerical simulation of seismic converted waves according to claim 1, wherein the step 3 of solving the modified elastic wave equation to obtain the Ss wave field and the Sp wave field specifically comprises: solving equations (4-1) - (4-9) by finite difference method to obtain the expressions of P wave and S wave, space 6 order:
wherein i, j represents coordinates in x, y directions of spatial grid points involved in calculation, corresponding to i, j in the context, and a (n) represents a difference coefficient; ui represents a P wave, uj represents an S wave, and DH represents a grid length; and realizing the expression of P wave and S wave by a program to obtain a converted wave shear wave field and a converted wave longitudinal wave field.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711483283.XA CN108563802B (en) | 2017-12-29 | 2017-12-29 | Method for improving numerical simulation precision of seismic converted waves |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711483283.XA CN108563802B (en) | 2017-12-29 | 2017-12-29 | Method for improving numerical simulation precision of seismic converted waves |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108563802A CN108563802A (en) | 2018-09-21 |
CN108563802B true CN108563802B (en) | 2021-12-17 |
Family
ID=63529547
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711483283.XA Active CN108563802B (en) | 2017-12-29 | 2017-12-29 | Method for improving numerical simulation precision of seismic converted waves |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108563802B (en) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104122585A (en) * | 2014-08-08 | 2014-10-29 | 中国石油大学(华东) | Seismic forward modeling method based on elastic wave field vector decomposition and low-rank decomposition |
CN104330822A (en) * | 2014-10-23 | 2015-02-04 | 中国石油天然气股份有限公司 | Method and device for determining remaining oil and gas distribution through coupling four-dimensional seismic inversion |
CN104635260A (en) * | 2013-11-08 | 2015-05-20 | 北京瑞拉爱堡地质勘探技术有限公司 | Seismic data prestack imaging method |
CN105242305A (en) * | 2015-09-06 | 2016-01-13 | 中国科学院地质与地球物理研究所 | Longitudinal wave and transverse wave separation method and system |
CN105467443A (en) * | 2015-12-09 | 2016-04-06 | 中国科学院地质与地球物理研究所 | A three-dimensional anisotropy elastic wave numerical simulation method and system |
CN105891884A (en) * | 2016-06-28 | 2016-08-24 | 中国石油大学(北京) | Micro-earthquake focus mechanism inversion method and micro-earthquake focus mechanism inversion device |
CN106680870A (en) * | 2017-03-16 | 2017-05-17 | 中国海洋大学 | High-precision seismic travel time ray tracing method |
CN106970416A (en) * | 2017-03-17 | 2017-07-21 | 中国地质科学院地球物理地球化学勘查研究所 | Elastic wave least square reverse-time migration system and method based on wave field separation |
CN107462922A (en) * | 2016-06-06 | 2017-12-12 | 中国石油化工股份有限公司 | full waveform inversion method and device |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150081223A1 (en) * | 2013-09-19 | 2015-03-19 | Schlumberger Technology Corporation | Microseismic survey |
EP3209859B1 (en) * | 2014-10-24 | 2021-04-28 | Schlumberger Technology B.V. | Travel-time objective function for full waveform inversion |
-
2017
- 2017-12-29 CN CN201711483283.XA patent/CN108563802B/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104635260A (en) * | 2013-11-08 | 2015-05-20 | 北京瑞拉爱堡地质勘探技术有限公司 | Seismic data prestack imaging method |
CN104122585A (en) * | 2014-08-08 | 2014-10-29 | 中国石油大学(华东) | Seismic forward modeling method based on elastic wave field vector decomposition and low-rank decomposition |
CN104330822A (en) * | 2014-10-23 | 2015-02-04 | 中国石油天然气股份有限公司 | Method and device for determining remaining oil and gas distribution through coupling four-dimensional seismic inversion |
CN105242305A (en) * | 2015-09-06 | 2016-01-13 | 中国科学院地质与地球物理研究所 | Longitudinal wave and transverse wave separation method and system |
CN105467443A (en) * | 2015-12-09 | 2016-04-06 | 中国科学院地质与地球物理研究所 | A three-dimensional anisotropy elastic wave numerical simulation method and system |
CN107462922A (en) * | 2016-06-06 | 2017-12-12 | 中国石油化工股份有限公司 | full waveform inversion method and device |
CN105891884A (en) * | 2016-06-28 | 2016-08-24 | 中国石油大学(北京) | Micro-earthquake focus mechanism inversion method and micro-earthquake focus mechanism inversion device |
CN106680870A (en) * | 2017-03-16 | 2017-05-17 | 中国海洋大学 | High-precision seismic travel time ray tracing method |
CN106970416A (en) * | 2017-03-17 | 2017-07-21 | 中国地质科学院地球物理地球化学勘查研究所 | Elastic wave least square reverse-time migration system and method based on wave field separation |
Non-Patent Citations (3)
Title |
---|
The oceanic nature of the African slab subducted under Peloponnesus: thin-layer resolution from multiscale analysis of teleseismic P-to-S converted waves;A. Gesret 等;《 Geophysical Journal International 》;20101002;833 - 849 * |
单程波动方程地震数值模拟新方法研究;熊晓军;《中国博士学位论文全文数据库》;20070401;A011-37 * |
滩浅海高精度地震资料处理研究;王福海,童思友 等;《青岛海洋大学学报》;20031230;32-36 * |
Also Published As
Publication number | Publication date |
---|---|
CN108563802A (en) | 2018-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108873066B (en) | Elastic medium wave equation reflected wave travel time inversion method | |
CN105467444B (en) | A kind of elastic wave full waveform inversion method and device | |
WO2012170201A2 (en) | System and method for seismic data inversion by non-linear model update | |
CN106133556A (en) | Determine the component of wave field | |
CN110187382B (en) | Traveling time inversion method for wave equation of reverse wave and reflected wave | |
CN105549068A (en) | 3D anisotropic micro seismic interference inverse positioning method and 3D anisotropic micro seismic interference inverse positioning system | |
CN103412328A (en) | Wave number field amplitude preservation wave field separation method based on staggered mesh finite difference algorithm | |
Araujo et al. | Symplectic scheme and the Poynting vector in reverse-time migration | |
Gao et al. | Full waveform inversion using deterministic source encoding | |
CN111045077A (en) | Full waveform inversion method of land seismic data | |
Maharramov et al. | Robust joint full-waveform inversion of time-lapse seismic data sets with total-variation regularization | |
DU et al. | Hybrid PS/FD numerical simulation and stability analysis of pure P-wave propagation in VTI media | |
CN104597489A (en) | Seismic source wavelet optimal setting method and device | |
NO20190489A1 (en) | Seismic modeling | |
Albertin et al. | 3D orthorhombic elastic full-waveform inversion in the reflection domain from hydrophone data | |
CN106054242A (en) | Three dimensional anisotropy attenuation medium wave field simulation method | |
Sharan et al. | Sparsity-promoting joint microseismic source collocation and source-time function estimation | |
CN111505714B (en) | Elastic wave direct envelope inversion method based on rock physical constraint | |
Zhang et al. | Massive 3D seismic data compression and inversion with hierarchical Tucker | |
CN108563802B (en) | Method for improving numerical simulation precision of seismic converted waves | |
Huang et al. | Waveform-based source location method using a source parameter isolation strategy | |
Raknes et al. | Challenges and solutions for performing 3D time-domain elastic full-waveform inversion | |
Jiang et al. | Full waveform inversion based on inversion network reparameterized velocity | |
Masmoudi et al. | Waveform inversion in acoustic orthorhombic media with a practical set of parameters | |
CN113866827A (en) | Method, system, medium and device for explanatory velocity modeling seismic imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |