CN106908839B - A kind of seismic wave accurate reflection coefficient analysis method and device - Google Patents

A kind of seismic wave accurate reflection coefficient analysis method and device Download PDF

Info

Publication number
CN106908839B
CN106908839B CN201710312357.7A CN201710312357A CN106908839B CN 106908839 B CN106908839 B CN 106908839B CN 201710312357 A CN201710312357 A CN 201710312357A CN 106908839 B CN106908839 B CN 106908839B
Authority
CN
China
Prior art keywords
wave
velocity
reflection coefficient
medium
longitudinal
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
Application number
CN201710312357.7A
Other languages
Chinese (zh)
Other versions
CN106908839A (en
Inventor
宗兆云
印兴耀
程广森
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum East China
Original Assignee
China University of Petroleum East China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN201710312357.7A priority Critical patent/CN106908839B/en
Publication of CN106908839A publication Critical patent/CN106908839A/en
Application granted granted Critical
Publication of CN106908839B publication Critical patent/CN106908839B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/67—Wave propagation modeling

Abstract

The invention discloses a kind of seismic wave accurate reflection coefficient analysis method and device, wherein this method includes：Obtain incident compressional angle θp1, longitudinal wave angle of transmission θp2, transverse wave reflection angle θs1, shear wave angle of transmission θs2, velocity of longitudinal wave vp, shear wave velocity vs, longitudinal wave reflection coefficient rp, transverse wave reflection coefficient rs and density reflectance factor rd；Seismic wave reflectance factor is parsed according to the following formula：Realize that the parsing of seismic reflection coefficient has the characteristics that precision is higher according to technical solution disclosed in the present application, further, it uses and establishes the processes of AVO forward simulations and can avoid the application range of limitation to(for) tri- parametric inversions of conventional AVO, to solve problems of the prior art.

Description

A kind of seismic wave accurate reflection coefficient analysis method and device
Technical field
The present invention relates to post-stack inversion technical fields, are parsed more specifically to a kind of seismic wave accurate reflection coefficient Method and device.
Background technology
Post-stack inversion utilizes poststack seismic data, main inverting stratum longitudinal wave information.Prestack seismic data contains many Information not available for poststack data can also estimate that formation shear, density, Modulus of Rocks, fluid are quick in addition to inverting longitudinal wave information The information such as parameter are felt, so pre-stack seismic inversion is by extensively should be in reservoir prediction and fluid identification.
Tri- parametric inversions of conventional AVO usually utilize Zoeppritz equations approximate expression (Aki-Richards approximations and Fatti approximations etc.) establish the process of AVO forward simulations, however in P wave incidence angles cross critical angle and elastic parameter variation is violent In the case of, approximate formula precision is limited, limits the application range of tri- parametric inversions of conventional AVO.
It is current in conclusion how to provide a kind of higher technical solution of precision to establish the process of AVO forward simulations Those skilled in the art's urgent problem to be solved.
Invention content
It is anti-to improve seismic wave the object of the present invention is to provide a kind of seismic wave accurate reflection coefficient analysis method and device Penetrate the precision of coefficient parsing.
To achieve the goals above, the present invention provides the following technical solutions：
A kind of seismic wave accurate reflection coefficient analytic method, including：
Obtain incident compressional angle θp1, longitudinal wave angle of transmission θp2, transverse wave reflection angle θs1, shear wave angle of transmission θs2, velocity of longitudinal wave vp、 Shear wave velocity vs, longitudinal wave reflection coefficient rp, transverse wave reflection coefficient rs and density reflectance factor rd；
Seismic wave reflectance factor is parsed according to the following formula：
Wherein,
Preferably, velocity of longitudinal wave v is obtainedpAnd shear wave velocity vs, including：
The velocity of longitudinal wave v of first medium is obtained respectivelyp1, shear wave velocity vs1, the velocity of longitudinal wave v of second mediump2, shear wave speed Spend vs2And the velocity of longitudinal wave difference Δ v between first medium and second mediump, shear wave velocity difference Δ vs；
Velocity of longitudinal wave v is calculated separately according to the following formulapAnd shear wave velocity vs：
Preferably, longitudinal wave reflection coefficient rp, transverse wave reflection coefficient rs and density reflectance factor rd are obtained, including：
The density p of the first medium is obtained respectively1, second medium density p1Between first medium and second medium Density contrast Δ ρ；
Longitudinal wave reflection coefficient rp, transverse wave reflection coefficient rs and density reflectance factor rd are calculated separately according to the following formula：
Wherein,
A kind of seismic wave accurate reflection coefficient resolver, including：
Acquisition module is used for：Obtain incident compressional angle θp1, longitudinal wave angle of transmission θp2, transverse wave reflection angle θs1, shear wave angle of transmission θs2, velocity of longitudinal wave vp, shear wave velocity vs, longitudinal wave reflection coefficient rp, transverse wave reflection coefficient rs and density reflectance factor rd；
Parsing module is used for：Seismic wave reflectance factor is parsed according to the following formula：
Wherein,
Preferably, the acquisition module includes：
First acquisition unit is used for：The velocity of longitudinal wave v of first medium is obtained respectivelyp1, shear wave velocity vs1, second medium Velocity of longitudinal wave vp2, shear wave velocity vs2And the velocity of longitudinal wave difference Δ v between first medium and second mediump, shear wave velocity difference Δ vs；Velocity of longitudinal wave v is calculated separately according to the following formulapAnd shear wave velocity vs：
Preferably, the acquisition module includes：
Second acquisition unit is used for：The density p of the first medium is obtained respectively1, second medium density p1With first Density contrast Δ ρ between medium and second medium；Longitudinal wave reflection coefficient rp, transverse wave reflection coefficient are calculated separately according to the following formula Rs and density reflectance factor rd：
Wherein,
The present invention provides a kind of seismic wave accurate reflection coefficient analysis method and devices, it is demonstrated experimentally that according to the application Disclosed technical solution realizes that the parsing of seismic reflection coefficient has the characteristics that precision is higher, further, uses foundation The process of AVO forward simulations can avoid the limitation of the application range for tri- parametric inversions of conventional AVO, existing to solve The problem of technology.
Description of the drawings
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below There is attached drawing needed in technology description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this The embodiment of invention for those of ordinary skill in the art without creative efforts, can also basis The attached drawing of offer obtains other attached drawings.
Fig. 1 is a kind of flow chart of seismic wave accurate reflection coefficient analytic method provided in an embodiment of the present invention；
Fig. 2 is Goodway models essence in a kind of seismic wave accurate reflection coefficient analytic method provided in an embodiment of the present invention Spend inspection result schematic diagram；
Fig. 3 is Ostrander models in a kind of seismic wave accurate reflection coefficient analytic method provided in an embodiment of the present invention Accuracy test result schematic diagram；
Fig. 4 is a kind of structural schematic diagram of seismic wave accurate reflection coefficient resolver provided in an embodiment of the present invention.
Specific implementation mode
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete Site preparation describes, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other Embodiment shall fall within the protection scope of the present invention.
Referring to Fig. 1, it illustrates a kind of seismic wave accurate reflection coefficient analytic methods provided in an embodiment of the present invention Flow chart may include：
S11：Obtain incident compressional angle θp1, longitudinal wave angle of transmission θp2, transverse wave reflection angle θs1, shear wave angle of transmission θs2, velocity of longitudinal wave vp, shear wave velocity vs, longitudinal wave reflection coefficient rp, transverse wave reflection coefficient rs and density reflectance factor rd.
It should be noted that above-mentioned each concept is consistent with the meaning for corresponding to concept in the prior art, details are not described herein.
S12：Seismic wave reflectance factor is parsed according to the following formula：
Wherein,
It is demonstrated experimentally that realizing that the parsing of seismic reflection coefficient has precision according to above-mentioned technical proposal disclosed in the present application Higher feature further uses and establishes the processes of AVO forward simulations and can avoid for tri- parametric inversions of conventional AVO Application range limitation, to solve problems of the prior art.
In addition in above-mentioned technical proposal disclosed in the present application for realizing seismic reflection coefficient parsing formula by six Independent seismologic parameter is reduced to four, i.e., respectively velocity of longitudinal wave reflectance factor, shear wave velocity reflectance factor, density reflection are Number and P-S wave velocity ratio, wherein P-S wave velocity ratio may be configured as constant to greatly facilitate accurate Zoeppritz equations Application and practice.
With reference to the essence of classical Goodway models and Ostrander model testings above-mentioned formula disclosed in the present application Exactness (it is assumed that known to P-S wave velocity ratio), above-mentioned formula precision disclosed in the present application is higher it can be seen from Fig. 2 and Fig. 3, right The method for the above-mentioned realization seismic reflection coefficient parsing disclosed in the present application answered has correctness.
In a kind of seismic wave accurate reflection coefficient analytic method provided in an embodiment of the present invention, velocity of longitudinal wave v is obtainedpAnd it is horizontal Wave velocity vs, may include：
The velocity of longitudinal wave v of first medium is obtained respectivelyp1, shear wave velocity vs1, the velocity of longitudinal wave v of second mediump2, shear wave speed Spend vs2And the velocity of longitudinal wave difference Δ v between first medium and second mediump, shear wave velocity difference Δ vs；
Velocity of longitudinal wave v is calculated separately according to the following formulapAnd shear wave velocity vs：
A kind of seismic wave accurate reflection coefficient analytic method provided in an embodiment of the present invention, acquisition longitudinal wave reflection coefficient rp, Transverse wave reflection coefficient rs and density reflectance factor rd may include：
The density p of first medium is obtained respectively1, second medium density p1It is close between first medium and second medium The poor Δ ρ of degree；
Longitudinal wave reflection coefficient rp, transverse wave reflection coefficient rs and density reflectance factor rd are calculated separately according to the following formula：
Wherein,
Corresponding parameter can be simply and rapidly got through the above way, ensure that on provided in an embodiment of the present invention State the realization efficiency of technical solution.
Specifically, the application to the adjustment of accurate Zoeppritz equations (1) first by obtaining a new equation (2)：
Rpp=(E-F+-G+H-p2)/D (1)
Wherein,
D=E+F++G-H-p2, a=ρ2(1-2vs2 2p2)-ρ1(1-2vs1 2p2)
B=ρ2(1-2vs2 2p2)+2ρ1vs1 2p2, c=ρ1(1-2vs1 2p2)+2ρ2vs2 2p2
D=2 ρ2vs2 2-2ρ1vs1 2,
Wherein,
Wherein, θp1,θp2,θs1,θs2Incident compressional angle, longitudinal wave angle of transmission, transverse wave reflection angle and shear wave transmission are indicated respectively Angle, vp1,vs1,ρ1,vp2,vs2,ρ2It indicates respectively in velocity of longitudinal wave, shear wave velocity and density and the second medium in first medium Velocity of longitudinal wave, shear wave velocity and density.Wherein first medium and second medium are in neighbouring two layer medium.
Pass throughRelationship can obtain formula (3)：
Formula (4)-(10) can similarly be obtained：
Wherein, Δ vp,Δvs, Δ ρ is the poor velocity of longitudinal wave of lower upper two layer medium, shear wave velocity difference and density contrast respectively.
It is longitudinal wave reflection coefficient, transverse wave reflection coefficient, density reflection system respectively Number and P-S wave velocity ratio.
Formula (3)-(10), which are updated to, can obtain a new Exact Equation (11) in equation (2)
Wherein,
The embodiment of the present invention additionally provides a kind of seismic wave accurate reflection coefficient resolver, as shown in figure 4, can wrap It includes：
Acquisition module 11, is used for：Obtain incident compressional angle θp1, longitudinal wave angle of transmission θp2, transverse wave reflection angle θs1, shear wave transmission Angle θs2, velocity of longitudinal wave vp, shear wave velocity vs, longitudinal wave reflection coefficient rp, transverse wave reflection coefficient rs and density reflectance factor rd；
Parsing module 12, is used for：Seismic wave reflectance factor is parsed according to the following formula：
Wherein,
A kind of seismic wave accurate reflection coefficient resolver provided in an embodiment of the present invention, acquisition module may include：
First acquisition unit is used for：The velocity of longitudinal wave v of first medium is obtained respectivelyp1, shear wave velocity vs1, second medium Velocity of longitudinal wave vp2, shear wave velocity vs2And the velocity of longitudinal wave difference Δ v between first medium and second mediump, shear wave velocity difference Δ vs；Velocity of longitudinal wave v is calculated separately according to the following formulapAnd shear wave velocity vs：
A kind of seismic wave accurate reflection coefficient resolver provided in an embodiment of the present invention, acquisition module may include：
Second acquisition unit is used for：The density p of first medium is obtained respectively1, second medium density p1And first medium Density contrast Δ ρ between second medium；Calculate separately according to the following formula longitudinal wave reflection coefficient rp, transverse wave reflection coefficient rs and Density reflectance factor rd：
Wherein,
The explanation of relevant portion please join in a kind of seismic wave accurate reflection coefficient resolver provided in an embodiment of the present invention See the detailed description of corresponding part in a kind of seismic wave accurate reflection coefficient analytic method provided in an embodiment of the present invention, herein not It repeats again.
The foregoing description of the disclosed embodiments enables those skilled in the art to realize or use the present invention.To this A variety of modifications of a little embodiments will be apparent for a person skilled in the art, and the general principles defined herein can Without departing from the spirit or scope of the present invention, to realize in other embodiments.Therefore, the present invention will not be limited It is formed on the embodiments shown herein, and is to fit to consistent with the principles and novel features disclosed in this article widest Range.

Claims (6)

1. a kind of seismic wave accurate reflection coefficient analytic method, which is characterized in that including：
Obtain incident compressional angle θp1, longitudinal wave angle of transmission θp2, transverse wave reflection angle θs1, shear wave angle of transmission θs2, velocity of longitudinal wave vp, shear wave Speed vs, longitudinal wave reflection coefficient rp, transverse wave reflection coefficient rs and density reflectance factor rd；
Seismic wave reflectance factor is parsed according to the following formula：
Wherein,
2. according to the method described in claim 1, it is characterized in that, obtaining velocity of longitudinal wave vpAnd shear wave velocity vs, including：
The velocity of longitudinal wave v of first medium is obtained respectivelyp1, shear wave velocity vs1, the velocity of longitudinal wave v of second mediump2, shear wave velocity vs2, And the velocity of longitudinal wave difference Δ v between first medium and second mediump, shear wave velocity difference Δ vs；
Velocity of longitudinal wave v is calculated separately according to the following formulapAnd shear wave velocity vs：
3. according to the method described in claim 2, it is characterized in that, obtain longitudinal wave reflection coefficient rp, transverse wave reflection coefficient rs and Density reflectance factor rd, including：
The density p of the first medium is obtained respectively1, second medium density p2It is close between first medium and second medium The poor Δ ρ of degree；
Longitudinal wave reflection coefficient rp, transverse wave reflection coefficient rs and density reflectance factor rd are calculated separately according to the following formula：
Wherein,
4. a kind of seismic wave accurate reflection coefficient resolver, which is characterized in that including：
Acquisition module is used for：Obtain incident compressional angle θp1, longitudinal wave angle of transmission θp2, transverse wave reflection angle θs1, shear wave angle of transmission θs2, it is vertical Wave velocity vp, shear wave velocity vs, longitudinal wave reflection coefficient rp, transverse wave reflection coefficient rs and density reflectance factor rd；
Parsing module is used for：Seismic wave reflectance factor is parsed according to the following formula：
Wherein,
5. device according to claim 4, which is characterized in that the acquisition module includes：
First acquisition unit is used for：The velocity of longitudinal wave v of first medium is obtained respectivelyp1, shear wave velocity vs1, the longitudinal wave of second medium Speed vp2, shear wave velocity vs2And the velocity of longitudinal wave difference Δ v between first medium and second mediump, shear wave velocity difference Δ vs；It presses Velocity of longitudinal wave v is calculated separately according to following equationpAnd shear wave velocity vs：
6. device according to claim 5, which is characterized in that the acquisition module includes：
Second acquisition unit is used for：The density p of the first medium is obtained respectively1, second medium density p2And first medium Density contrast Δ ρ between second medium；Calculate separately according to the following formula longitudinal wave reflection coefficient rp, transverse wave reflection coefficient rs and Density reflectance factor rd：
Wherein,
CN201710312357.7A 2017-05-05 2017-05-05 A kind of seismic wave accurate reflection coefficient analysis method and device Active CN106908839B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710312357.7A CN106908839B (en) 2017-05-05 2017-05-05 A kind of seismic wave accurate reflection coefficient analysis method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710312357.7A CN106908839B (en) 2017-05-05 2017-05-05 A kind of seismic wave accurate reflection coefficient analysis method and device

Publications (2)

Publication Number Publication Date
CN106908839A CN106908839A (en) 2017-06-30
CN106908839B true CN106908839B (en) 2018-09-28

Family

ID=59210322

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710312357.7A Active CN106908839B (en) 2017-05-05 2017-05-05 A kind of seismic wave accurate reflection coefficient analysis method and device

Country Status (1)

Country Link
CN (1) CN106908839B (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101957454A (en) * 2010-09-17 2011-01-26 中国石油大学(华东) AVO (Amplitude Versus Offset) approximate formula based on incidence angle and attribute extraction method
CN102590863A (en) * 2012-01-29 2012-07-18 中国石油天然气股份有限公司 Methods and devices for determining reflection coefficient dispersion and reflection energy and identifying gas layer
CN103257361A (en) * 2013-05-24 2013-08-21 中国石油天然气集团公司 Petroleum-gas prediction method and system based on Zoeppritz equation approximate expression
CN103675911A (en) * 2014-01-06 2014-03-26 中国石油大学(华东) Method for combined inversion of intercept and gradient based on compressional waves and converted shear waves
CN104570072A (en) * 2013-10-16 2015-04-29 中国石油化工股份有限公司 Method for modeling reflection coefficient of spherical PP wave in viscoelastic medium
CN104965224A (en) * 2015-06-03 2015-10-07 北京多分量地震技术研究院 PP wave and PS wave joint AVO inversion method using average angle of incidence gather
CN104977610A (en) * 2015-06-23 2015-10-14 中国石油大学(华东) Method of attribute extraction through AVO approximate formula based on incidence angle

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2920008C (en) * 2013-07-29 2017-06-13 Cgg Services Sa Method and device for the generation and application of anisotropic elastic parameters in horizontal transverse isotropic (hti) media

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101957454A (en) * 2010-09-17 2011-01-26 中国石油大学(华东) AVO (Amplitude Versus Offset) approximate formula based on incidence angle and attribute extraction method
CN102590863A (en) * 2012-01-29 2012-07-18 中国石油天然气股份有限公司 Methods and devices for determining reflection coefficient dispersion and reflection energy and identifying gas layer
CN103257361A (en) * 2013-05-24 2013-08-21 中国石油天然气集团公司 Petroleum-gas prediction method and system based on Zoeppritz equation approximate expression
CN104570072A (en) * 2013-10-16 2015-04-29 中国石油化工股份有限公司 Method for modeling reflection coefficient of spherical PP wave in viscoelastic medium
CN103675911A (en) * 2014-01-06 2014-03-26 中国石油大学(华东) Method for combined inversion of intercept and gradient based on compressional waves and converted shear waves
CN104965224A (en) * 2015-06-03 2015-10-07 北京多分量地震技术研究院 PP wave and PS wave joint AVO inversion method using average angle of incidence gather
CN104977610A (en) * 2015-06-23 2015-10-14 中国石油大学(华东) Method of attribute extraction through AVO approximate formula based on incidence angle

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"一种计算反射/传输系数的快速精确算法";法林 等;《中国科学》;20081231;第38卷(第7期);第873-895页 *

Also Published As

Publication number Publication date
CN106908839A (en) 2017-06-30

Similar Documents

Publication Publication Date Title
US10578757B2 (en) Omnidirectional vector seismic data processing method and apparatus, computer-readable storage medium and device
CN103149586B (en) Wave field the Forward Modeling in a kind of inclination stratified viscoelastic media
EP3028071B1 (en) Method and device for the generation and application of anisotropic elastic parameters in horizontal transverse isotropic (hti) media
US20170299745A1 (en) Prestack egs migration method for seismic wave multi-component data
CN104570072A (en) Method for modeling reflection coefficient of spherical PP wave in viscoelastic medium
CN105182408A (en) Manufacturing method and device for synthesizing earthquake record
CN1013226B (en) Method for depth imaging multicomponent seismic data
CN105093292A (en) Data processing method and device for earthquake imaging
CN102156297B (en) Fluid substitution method based on sandstone reservoir post-stack seismic data
CN105068117A (en) AVO (Amplitude Versus Offset) retrieval method, device and equipment for fractured medium
US20120010820A1 (en) Fresnel Zone Fat Ray Tomography
CN107589443A (en) Method and system based on elastic wave least square reverse-time migration imaging
CN102866426B (en) A kind of method utilizing AVO wide-angle road set analysis rock mass hydrocarbon information
CN102105900A (en) Method for propagating pseudo acoustic quasi-P waves in anisotropic media
CN106199704B (en) A kind of Three-dimendimal fusion submarine cable seismic data velocity modeling method
CN105242313B (en) A kind of bearing calibration of elastic wave reverse-time migration polarity inversion and system
CN106324665A (en) Method and system of inverting fracture density
CN109143351B (en) Pre-stack anisotropy characteristic parameter inversion method and computer readable storage medium
CN106324666A (en) Extraction method and extraction device for extracting converted wave angle gather of transversely isotropic media
CN104977609A (en) Prestack longitudinal wave and transverse wave combined inversion method based on rapid simulated annealing
CN106908839B (en) A kind of seismic wave accurate reflection coefficient analysis method and device
Gao et al. An Efficient Multiscale Finite‐Element Method for Frequency‐Domain Seismic Wave Propagation
CN104133249B (en) A kind of micro-logging data and the method and device of Sonic Logging Data joint interpretation
Chen et al. Anisotropic layering and seismic body waves: Deformation gradients, initial S-polarizations, and converted-wave birefringence
CN104898162B (en) Crack detection method in geological exploration

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