CN112558151A - AVA waveform simulation method and system in attenuation medium - Google Patents

AVA waveform simulation method and system in attenuation medium Download PDF

Info

Publication number
CN112558151A
CN112558151A CN202011530847.2A CN202011530847A CN112558151A CN 112558151 A CN112558151 A CN 112558151A CN 202011530847 A CN202011530847 A CN 202011530847A CN 112558151 A CN112558151 A CN 112558151A
Authority
CN
China
Prior art keywords
frequency
attenuation
seismic
reflection coefficient
ava
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.)
Pending
Application number
CN202011530847.2A
Other languages
Chinese (zh)
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.)
Tongji University
Original Assignee
Tongji University
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 Tongji University filed Critical Tongji University
Priority to CN202011530847.2A priority Critical patent/CN112558151A/en
Publication of CN112558151A publication Critical patent/CN112558151A/en
Pending legal-status Critical Current

Links

Images

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/282—Application of seismic models, synthetic seismograms
    • 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

Abstract

The invention relates to an AVA waveform simulation method and system in an attenuation medium, which comprises the following steps: constructing an attenuation medium rock physical model, and calculating the seismic wave velocity changing along with the frequency; acquiring a reflection coefficient which changes along with the frequency and the incident angle; acquiring a quality factor representing the propagation energy attenuation of the seismic waves, and constructing an attenuation function according to the quality factor; constructing an unsteady reflection coefficient according to the attenuation function and the reflection coefficient which changes along with the frequency and the incidence angle; and performing frequency decomposition on the seismic wavelet, performing convolution on each single-frequency component of the seismic wavelet and a single-frequency unsteady reflection coefficient, and superposing all frequencies to obtain AVA waveform analog record in the attenuation medium. In addition, the invention also discloses a system. Compared with the prior art, the method can simulate AVA waveform recording in the attenuation medium under the condition of the same calculation complexity, and the simulated seismic record simultaneously contains seismic wave propagation attenuation effect and reflection attenuation effect and can be accurately matched with the actual seismic waveform recording.

Description

AVA waveform simulation method and system in attenuation medium
Technical Field
The invention relates to the field of seismic wave forward modeling in exploration seismology, in particular to an AVA waveform simulation method and system in an attenuation medium.
Background
The synthetic seismic recording method based on the convolution model is the most widely applied seismic wave forward modeling method due to the relative simplicity of the calculation degree. At present, an AVA inversion method based on a convolution model is widely adopted in the industry to carry out oil and gas reservoir characterization, and the AVA waveform simulation method based on the convolution model is based on the AVA waveform simulation method, but the existing AVA seismic waveform simulation method based on the convolution model does not consider the seismic wave energy propagation attenuation effect and the seismic wave reflection attenuation effect at the same time, so that the simulated AVA waveform and the actual seismic record waveform have larger errors, and the reliability of AVA inversion is reduced.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide an AVA waveform simulation method and system in an attenuation medium.
The purpose of the invention can be realized by the following technical scheme:
a method for simulating an AVA waveform in an attenuation medium comprises the following steps:
s1: acquiring logging data, constructing an attenuation medium rock physical model, establishing a speed model changing along with frequency according to the attenuation medium rock model, and calculating the seismic wave speed changing along with frequency;
s2: acquiring a reflection coefficient which changes along with the frequency and the incident angle according to a speed model which changes along with the frequency and a pore elastic medium theory;
s3: acquiring a quality factor representing the propagation energy attenuation of the seismic waves according to the seismic wave speed changing along with the frequency, and constructing an attenuation function according to the quality factor;
s4: constructing an unsteady reflection coefficient according to the attenuation function and the reflection coefficient which changes along with the frequency and the incidence angle;
s5: and performing frequency decomposition on the seismic wavelet, performing convolution on each single-frequency component of the seismic wavelet and a single-frequency unsteady reflection coefficient, and then stacking all frequencies to obtain AVA waveform analog record in the attenuation medium.
Preferably, the rock physical model of the attenuation medium is as follows:
Q-1(ω)＝2Im(kp)/RE(kp)
wherein Q is a quality factor, kpIs the longitudinal wave number and omega is the seismic wave frequency.
Preferably, the velocity model with frequency variation is:
v(ω)＝ω/RE(kp)
where v (ω) is the seismic velocity, k, which varies with frequencypIs the longitudinal wave number and omega is the seismic wave frequency.
Preferably, the quality factor is calculated by:
Figure BDA0002852011340000021
wherein Q isave(t, ω) is the mean of the quality factors of the subsurface N layers, Qk(t, ω) is the k-th layer quality factor, Δ tkWhen the seismic wave propagation double-travel in the k-th underground layer is carried out, N is the total number of underground layers, and k is the number of underground layers.
Preferably, the attenuation function is:
α(t,ω)＝e-πωt/Q(t,ω)
wherein, α (t, ω) is an attenuation function, t is a seismic wave two-way travel time, ω is a seismic wave frequency, Q (t, ω) is a quality factor, and e is a natural constant.
Preferably, in S4, the attenuation function is multiplied by the reflection coefficient varying with frequency and incidence angle to obtain the unsteady reflection coefficient, where the unsteady reflection coefficient is:
Figure BDA0002852011340000022
wherein the content of the first and second substances,
Figure BDA0002852011340000023
for unsteady reflection coefficients, α (t, ω) is the attenuation function, and c (t, ω) is the reflection coefficient as a function of frequency and angle of incidence.
Preferably, in S5, the formula for obtaining the analog record of the AVA waveform in the attenuation medium is:
Figure BDA0002852011340000024
wherein d issFor AVA waveform analog recording in attenuating media, s is a seismic wavelet,
Figure BDA0002852011340000025
for unsteady reflection coefficients, ω is the frequency.
Preferably, the well log data comprises compressional and shear wave velocity, density, porosity, pore tortuosity, fluid water saturation, fluid density, fluid viscosity coefficient, rock mineral bulk modulus, rock shear modulus and rock permeability.
An AVA waveform simulation system in an attenuating medium, the system comprising:
a memory storing computer-executable instructions;
a processor executing computer executable instructions in the memory to perform the steps of:
acquiring logging data, constructing an attenuation medium rock physical model, establishing a speed model changing along with frequency according to the attenuation medium rock model, and calculating the seismic wave speed changing along with frequency;
acquiring a reflection coefficient which changes along with the frequency and the incident angle according to a speed model which changes along with the frequency and a pore elastic medium theory;
acquiring a quality factor representing the propagation energy attenuation of the seismic waves according to the seismic wave speed changing along with the frequency, and constructing an attenuation function according to the quality factor;
constructing an unsteady reflection coefficient according to the attenuation function and the reflection coefficient which changes along with the frequency and the incidence angle;
and performing frequency decomposition on the seismic wavelet, performing convolution on each single-frequency component of the seismic wavelet and a single-frequency unsteady reflection coefficient, and then stacking all frequencies to obtain AVA waveform analog record in the attenuation medium.
Compared with the prior art, the method is based on the rock physical model of the attenuation medium, the velocity dispersion relation is established, the change of the reflection coefficient along with the frequency is obtained, the attenuation effect and the reflection attenuation effect of seismic wave propagation energy are considered at the same time, the unsteady state reflection coefficient which takes the seismic wave propagation and the reflection attenuation effect along with the change of the frequency and the incidence angle is established, single-frequency AVA waveform simulation is realized in a frequency domain based on the convolution model, then all the frequencies are superposed, and the broadband AVA waveform simulation record is obtained. Compared with the traditional convolution model, the method can simulate AVA waveform recording in the attenuation medium under the condition of the same calculation complexity, the simulated seismic recording not only comprises seismic wave interface reflection attenuation effect, but also comprises seismic wave propagation attenuation effect, the seismic wave energy propagation attenuation effect and the seismic wave reflection attenuation effect are comprehensively considered, the actual seismic waveform recording is conveniently and accurately matched, the error is small, and the inversion reliability is high.
Drawings
FIG. 1 is a flow chart of a method of the present invention;
FIG. 2 is a schematic flow chart of the present invention;
FIG. 3 is a schematic diagram of the principles of the present invention;
figure 4 comparison of the simulation effect of the AVA waveform of the present invention.
Detailed Description
The invention is described in detail below with reference to the figures and specific embodiments. Note that the following description of the embodiments is merely a substantial example, and the present invention is not intended to be limited to the application or the use thereof, and is not limited to the following embodiments.
Examples
A method for simulating an AVA waveform in an attenuating medium, as shown in fig. 1 and 2, comprising the following steps:
s1: and acquiring logging data, constructing an attenuation medium rock physical model, establishing a speed model changing along with frequency according to the attenuation medium rock model, and calculating the seismic wave speed changing along with frequency.
Specifically, the attenuation medium petrophysical model is:
Q-1(ω)＝2Im(kp)/RE(kp)
wherein Q is a quality factor, kpIs the longitudinal wave number and omega is the seismic wave frequency.
And when specifically calculating:
Figure BDA0002852011340000041
in the above formula:
Figure BDA0002852011340000042
b0＝-β(Kd+4μ/3+α2/β)/α
Figure BDA0002852011340000043
c＝(α-bsρ/ρfb0)/(α+bs)
bs＝ρfθω2
Figure BDA0002852011340000044
α＝1-Kd/Ks
Figure BDA0002852011340000045
Figure BDA0002852011340000046
θ＝iκ(ω)/ηω
Figure BDA0002852011340000047
in the above formula, ω is the seismic frequency, K, KdRespectively, the bulk modulus, K, of the fluid saturated rock and the dry rocks、KfIs the bulk modulus of rock minerals and pore fluids, mu is the rock shear modulus,
Figure BDA0002852011340000048
is the porosity, ps、ρfRock mineral and pore fluid density, respectively, k0Is static permeability, tau is porosity tortuosity, eta is pore fluid viscosity coefficient, subscript +, -of parameter respectively represents fast and slow longitudinal wave, c, b0、bs、α、β、θ、kp0、
Figure BDA0002852011340000049
For calculating the parameter, i represents a complex number.
And, the velocity model with frequency variation is:
v(ω)＝ω/RE(kp)
where v (ω) is the seismic velocity, k, which varies with frequencypIs the longitudinal wave number and omega is the seismic wave frequency.
S2: and acquiring the reflection coefficient which changes along with the frequency and the incidence angle according to a speed model which changes along with the frequency and a pore elastic medium theory.
Specifically, in this embodiment, the velocity model varying with frequency is: c (t, ω) R11
Where t represents the two-way travel time of seismic wave propagation, R11Is calculated as follows：
PX＝B
Figure BDA0002852011340000051
Figure BDA0002852011340000052
x＝sin(α)
Figure BDA0002852011340000053
Figure BDA0002852011340000054
Figure BDA0002852011340000055
Figure BDA0002852011340000056
Γn(Z)＝1-2(nZ)2
Figure BDA0002852011340000057
R11In the calculation process of (a), β represents a Biot coefficient, λ represents a rock lame coefficient, M represents a fluid saturated rock longitudinal wave modulus, μ represents a rock shear modulus, subscripts M1 and 2 represent upper and lower media of a reflecting interface respectively, n 1,2 and 3 represent a fast longitudinal wave, a slow longitudinal wave and a transverse wave respectively, and other parameters are calculation parameters.
S3: and acquiring a quality factor representing the propagation energy attenuation of the seismic waves according to the seismic wave speed changing along with the frequency, and constructing an attenuation function according to the quality factor.
Specifically, the calculation formula of the quality factor is:
Figure BDA0002852011340000061
wherein Q isave(t, ω) is the mean of the quality factors of the subsurface N layers, Qk(t, ω) is the k-th layer quality factor, Δ tkWhen the seismic wave propagation double-travel in the k-th underground layer is carried out, N is the total number of underground layers, and k is the number of underground layers.
The constructed attenuation function is:
α(t,ω)＝e-πωt/Q(t,ω)
wherein, α (t, ω) is an attenuation function, t is a seismic wave two-way travel time, ω is a seismic wave frequency, Q (t, ω) is a quality factor, and e is a natural constant.
S4: an unsteady reflection coefficient is constructed from the attenuation function and the reflection coefficient as a function of frequency and angle of incidence.
In S4, the attenuation function is multiplied by the reflection coefficient varying with frequency and incidence angle to obtain the unsteady reflection coefficient, where the unsteady reflection coefficient is:
Figure BDA0002852011340000062
wherein the content of the first and second substances,
Figure BDA0002852011340000063
for unsteady reflection coefficients, α (t, ω) is the attenuation function, and c (t, ω) is the reflection coefficient as a function of frequency and angle of incidence.
S5: and performing frequency decomposition on the seismic wavelet, performing convolution on each single-frequency component of the seismic wavelet and a single-frequency unsteady reflection coefficient, and then stacking all frequencies to obtain AVA waveform analog record in the attenuation medium.
As shown in fig. 3, specifically, the formula for obtaining an analog record of the AVA waveform in an attenuated medium is:
Figure BDA0002852011340000064
wherein d issFor AVA waveform analog recording in attenuating media, s is a seismic wavelet,
Figure BDA0002852011340000065
for unsteady reflection coefficients, ω is the frequency.
In the present invention, the well log data includes compressional and shear wave velocity, density, porosity, pore tortuosity, fluid water saturation, fluid density, fluid viscosity coefficient, rock mineral bulk modulus, rock shear modulus, and rock permeability.
An AVA waveform simulation system in an attenuating medium, the system comprising:
a memory storing computer-executable instructions;
a processor executing computer executable instructions in the memory to perform the steps of:
acquiring logging data, constructing an attenuation medium rock physical model, establishing a speed model changing along with frequency according to the attenuation medium rock model, and calculating the seismic wave speed changing along with frequency;
acquiring a reflection coefficient which changes along with the frequency and the incident angle according to a speed model which changes along with the frequency and a pore elastic medium theory;
acquiring a quality factor representing the propagation energy attenuation of the seismic waves according to the seismic wave speed changing along with the frequency, and constructing an attenuation function according to the quality factor;
constructing an unsteady reflection coefficient according to the attenuation function and the reflection coefficient which changes along with the frequency and the incidence angle;
and performing frequency decomposition on the seismic wavelet, performing convolution on each single-frequency component of the seismic wavelet and a single-frequency unsteady reflection coefficient, and then stacking all frequencies to obtain AVA waveform analog record in the attenuation medium.
FIG. 4 is an AVA waveform simulation result in which a) an AVA waveform simulation record without any attenuation effects; b) analog recording of AVA waveform only containing reflection interface effect; c) the method comprises AVA waveform simulation recording of the reflection interface effect and the seismic wave propagation effect, and has the advantages of obviously smaller simulation effect error and high reliability.
The above embodiments are merely examples and do not limit the scope of the present invention. These embodiments may be implemented in other various manners, and various omissions, substitutions, and changes may be made without departing from the technical spirit of the present invention.

Claims (10)

1. A method for simulating an AVA waveform in an attenuation medium is characterized by comprising the following steps:
s1: acquiring logging data, constructing an attenuation medium rock physical model, establishing a speed model changing along with frequency according to the attenuation medium rock model, and calculating the seismic wave speed changing along with frequency;
s2: acquiring a reflection coefficient which changes along with the frequency and the incident angle according to a speed model which changes along with the frequency and a pore elastic medium theory;
s3: acquiring a quality factor representing the propagation energy attenuation of the seismic waves according to the seismic wave speed changing along with the frequency, and constructing an attenuation function according to the quality factor;
s4: constructing an unsteady reflection coefficient according to the attenuation function and the reflection coefficient which changes along with the frequency and the incidence angle;
s5: and performing frequency decomposition on the seismic wavelet, performing convolution on each single-frequency component of the seismic wavelet and a single-frequency unsteady reflection coefficient, and then stacking all frequencies to obtain AVA waveform analog record in the attenuation medium.
2. The method of claim 1, wherein the attenuation medium petrophysical model is:
Q-1(ω)＝2Im(kp)/RE(kp)
wherein Q is a quality factor, kpIs the longitudinal wave number and omega is the seismic wave frequency.
3. The method of claim 1 wherein the frequency dependent velocity model is:
v(ω)＝ω/RE(kp)
where v (ω) is the seismic velocity, k, which varies with frequencypIs the longitudinal wave number and omega is the seismic wave frequency.
4. The method of claim 1 wherein the figure of merit is calculated as:
Figure FDA0002852011330000011
wherein Q isave(t, ω) is the mean of the quality factors of the subsurface N layers, Qk(t, ω) is the k-th layer quality factor, Δ tkWhen the seismic wave propagation double-travel in the k-th underground layer is carried out, N is the total number of underground layers, and k is the number of underground layers.
5. The method of claim 1 wherein the decay function is:
α(t,ω)＝e-πωt/Q(t,ω)
wherein, α (t, ω) is an attenuation function, t is a seismic wave two-way travel time, ω is a seismic wave frequency, Q (t, ω) is a quality factor, and e is a natural constant.
6. The method of claim 1, wherein the attenuation function is multiplied by a reflection coefficient varying with frequency and incidence angle to obtain an unsteady reflection coefficient in S4, the unsteady reflection coefficient is:
Figure FDA0002852011330000021
wherein the content of the first and second substances,
Figure FDA0002852011330000022
for unsteady reflection coefficients, α (t, ω) is the attenuation function, and c (t, ω) is the reflection coefficient as a function of frequency and angle of incidence.
7. The method of claim 1, wherein in step S5, the formula for obtaining the analog record of the AVA waveform in the attenuation medium is:
Figure FDA0002852011330000023
wherein d issFor AVA waveform analog recording in attenuating media, s is a seismic wavelet,
Figure FDA0002852011330000024
for unsteady reflection coefficients, ω is the frequency.
8. The method of claim 1 wherein the well log data includes compressional-compressional velocity, density, porosity, pore tortuosity, fluid water saturation, fluid density, fluid viscosity coefficient, rock mineral bulk modulus, rock shear modulus, and rock permeability.
9. An AVA waveform simulation system in an attenuating medium, the system comprising:
a memory storing computer-executable instructions;
a processor executing computer executable instructions in the memory to perform the steps of:
acquiring logging data, constructing an attenuation medium rock physical model, establishing a speed model changing along with frequency according to the attenuation medium rock model, and calculating the seismic wave speed changing along with frequency;
acquiring a reflection coefficient which changes along with the frequency and the incident angle according to a speed model which changes along with the frequency and a pore elastic medium theory;
acquiring a quality factor representing the propagation energy attenuation of the seismic waves according to the seismic wave speed changing along with the frequency, and constructing an attenuation function according to the quality factor;
constructing an unsteady reflection coefficient according to the attenuation function and the reflection coefficient which changes along with the frequency and the incidence angle;
and performing frequency decomposition on the seismic wavelet, performing convolution on each single-frequency component of the seismic wavelet and a single-frequency unsteady reflection coefficient, and then stacking all frequencies to obtain AVA waveform analog record in the attenuation medium.
10. The system of claim 9, wherein the formula for obtaining the analog record of the AVA waveform in the attenuating medium is:
Figure FDA0002852011330000031
wherein d issFor AVA waveform analog recording in attenuating media, s is a seismic wavelet,
Figure FDA0002852011330000032
for unsteady reflection coefficients, ω is the frequency.
CN202011530847.2A 2020-12-22 2020-12-22 AVA waveform simulation method and system in attenuation medium Pending CN112558151A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011530847.2A CN112558151A (en) 2020-12-22 2020-12-22 AVA waveform simulation method and system in attenuation medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011530847.2A CN112558151A (en) 2020-12-22 2020-12-22 AVA waveform simulation method and system in attenuation medium

Publications (1)

Publication Number Publication Date
CN112558151A true CN112558151A (en) 2021-03-26

Family

ID=75030831

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011530847.2A Pending CN112558151A (en) 2020-12-22 2020-12-22 AVA waveform simulation method and system in attenuation medium

Country Status (1)

Country Link
CN (1) CN112558151A (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104570071A (en) * 2013-10-12 2015-04-29 中国石油化工股份有限公司 Bayesian linear AVA and AVF retrieval method of sticky sound medium
CN106249282A (en) * 2015-06-12 2016-12-21 中国石油化工股份有限公司 A kind of frequency domain seismic channel set creation method being applicable to AVAF inverting
CN108181654A (en) * 2018-01-18 2018-06-19 中国石油大学(北京) AVAF analogy methods and device based on multi-scale rock physical model
CN111123354A (en) * 2019-12-30 2020-05-08 中国石油大学(北京) Method and equipment for predicting dense gas layer based on frequency-dependent reflection amplitude attenuation

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104570071A (en) * 2013-10-12 2015-04-29 中国石油化工股份有限公司 Bayesian linear AVA and AVF retrieval method of sticky sound medium
CN106249282A (en) * 2015-06-12 2016-12-21 中国石油化工股份有限公司 A kind of frequency domain seismic channel set creation method being applicable to AVAF inverting
CN108181654A (en) * 2018-01-18 2018-06-19 中国石油大学(北京) AVAF analogy methods and device based on multi-scale rock physical model
CN111123354A (en) * 2019-12-30 2020-05-08 中国石油大学(北京) Method and equipment for predicting dense gas layer based on frequency-dependent reflection amplitude attenuation

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JINGKANG YANG 等: "A frequency-decomposed nonstationary convolutional model for amplitude-versus-angle-and-frequency forward waveform modeling in attenuative media", 《GEOPHYSICS》 *
孙卫涛 等: "孔隙介质弹性波频散-衰减理论模型", 《地球物理学进展》 *

Similar Documents

Publication Publication Date Title
KR101861060B1 (en) Multi-parameter inversion through offset dependent elastic fwi
Sams et al. The measurement of velocity dispersion and frequency-dependent intrinsic attenuation in sedimentary rocks
CN102466816B (en) Inversion method for stratum elasticity constant parameter of pre-stack seismic data
CN110133718B (en) Attenuation anisotropic fluid elastic impedance inversion method
CN108535775B (en) Non-stationary seismic data sound impedance inversion method and device
CN106125135B (en) Gas-bearing sandstone reservoir seismic response method for numerical simulation based on petrophysical model
CN106324675B (en) A kind of broad-band teleseismic wave impedance low-frequency information prediction technique and system
CN112987088B (en) Seepage medium seismic transverse wave numerical simulation and imaging method
CN112255680A (en) Seismic response simulation method for heterogeneous hydrate reservoir
CN110618450B (en) Intelligent gas-bearing property prediction method for tight reservoir based on rock physical modeling
Thiel et al. Comparison of acoustic and elastic full‐waveform inversion of 2D towed‐streamer data in the presence of salt
Kumar et al. Frequency-dependent AVO analysis using the scattering response of a layered reservoir
CN111427086B (en) Method for predicting maximum attenuation parameter of stratum by utilizing amplitude frequency component of longitudinal wave and transverse wave
CN112558151A (en) AVA waveform simulation method and system in attenuation medium
Xue et al. Seismic attenuation estimation using a complete ensemble empirical mode decomposition-based method
CN106249294A (en) A kind of reservoir detecting method of hydrocarbon
CN104714253A (en) AVO/AVA analysis method based on dispersion viscosity wave equation
CN102967884B (en) Reliability of wave impedance inversion data evaluation method and device
CN113176610B (en) Seismic data transmission loss compensation method based on unsteady state model
CN108196301B (en) Method and device for acquiring gather with amplitude changing along with offset
EP3887873A1 (en) Implicit property modeling
Marr et al. Exploration horizons from new seismic concepts of CDP and digital processing
Weibull et al. A comparison of the effects of layering and attenuation on the seismic response of gas reservoirs
CN107764697A (en) Gas potential detection method based on the progressive equation non-linear inversion of pore media
CN104316959B (en) Fluid identification based on equivalent fluid acoustic wave impedance

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20210326