CN111562613B - Method for analyzing seismic wave frequency-dependent reflection coefficient of thin reservoir or mutual reservoir model - Google Patents

Method for analyzing seismic wave frequency-dependent reflection coefficient of thin reservoir or mutual reservoir model Download PDF

Info

Publication number
CN111562613B
CN111562613B CN202010294977.4A CN202010294977A CN111562613B CN 111562613 B CN111562613 B CN 111562613B CN 202010294977 A CN202010294977 A CN 202010294977A CN 111562613 B CN111562613 B CN 111562613B
Authority
CN
China
Prior art keywords
reservoir
thin
seismic
wave
mutual
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
CN202010294977.4A
Other languages
Chinese (zh)
Other versions
CN111562613A (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.)
East China Institute of Technology
Original Assignee
East China Institute of Technology
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 East China Institute of Technology filed Critical East China Institute of Technology
Priority to CN202010294977.4A priority Critical patent/CN111562613B/en
Publication of CN111562613A publication Critical patent/CN111562613A/en
Application granted granted Critical
Publication of CN111562613B publication Critical patent/CN111562613B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/282Application of seismic models, synthetic seismograms
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30Analysis
    • G01V1/301Analysis for determining seismic cross-sections or geostructures

Abstract

The invention relates to a method for analyzing a seismic wave frequency-dependent reflection coefficient of a thin reservoir or mutual reservoir model, which comprises the following steps: step 1: describing an actual oil and gas reservoir by using a modified Biot-jet flow pore medium model, and describing a top plate layer and a bottom plate layer of the oil and gas reservoir by using a common single-phase medium model; step 2: solving and correcting a Biot-jet flow pore medium wave equation according to stress, strain and fluid pressure continuous conditions of a stratum interface to generate a seismic wave reflection coefficient of a thin (mutual) layer reservoir; and step 3: generating the reflection of the seismic waves at the reservoir interface along with the change of the frequency according to the seismic wave reflection coefficient of the thin (mutual) layer reservoir, and the relation between the transmission coefficient and the seismic wave incident angle; the method overcomes the defect that the traditional analysis method does not strictly consider the condition that the reservoir is a typical multiphase pore medium, can be used for analyzing the change of the reflection seismic amplitude of the thin (mutual) layer reservoir along with the frequency and the incidence angle, and provides method support for the seismic prediction of the oil and gas reservoir.

Description

Method for analyzing seismic wave frequency-dependent reflection coefficient of thin reservoir or mutual reservoir model
Technical Field
The invention relates to an amplitude interpretation technology in oil and gas reservoir reflection seismic exploration, in particular to a method for analyzing seismic wave frequency-dependent reflection coefficients of a thin reservoir or mutual reservoir model for correcting Biot-jet flow.
Background
In the development of oil gas exploration and development for more than 150 years, the demand of social development on energy is greatly met, and inexhaustible power is provided for the sustainable development of economy and society of various countries in the world. Reflection seismic amplitude interpretation technology (AVO technology and the like) is an important method for reservoir seismic prediction, and the core of the method is to study the change of the reflection coefficient of seismic reflection waves at a reservoir interface along with the incidence angle. The modern oil and gas seismic exploration develops towards reservoir fine description and micro-pore structure research, and a thin (mutual) layer reservoir prediction technology is a difficult technology. The reservoir is a typical pore medium, and the reservoir seismic prediction technical research based on the pore medium wave propagation model is more consistent with the actual medium form and receives more and more attention.
A plurality of evidences indicate that seismic waves have a high frequency dispersion effect in a reservoir, the phase velocity and the attenuation change along with the frequency, and further, the reflection coefficient also has a frequency variation effect. The frequency-dependent reflection coefficient provides more analyzable attributes for reservoir seismic prediction. Therefore, the method for analyzing and researching the frequency-dependent reflection coefficient of the seismic wave of the thin (mutual) reservoir model based on the pore medium wave propagation theory has important scientific research and application value for oil and gas seismic exploration and prediction.
For the reservoir seismic wave reflection coefficient analysis technology, a common method is to use a Zoeppritz equation based on a single-phase medium to analyze the change of the reflection coefficient along with the incident angle. However, the Zoeppritz equation is complex and not practical. In 1961, bortfeld proposed a simplified formula, but the relationship between Bortfeld and lithology was not clear and was not popularized. In 1983, hilterman divides the longitudinal wave reflection coefficient in the Bortfeld formula into a rigid factor and a fluid factor, gives the wave reflection coefficients of longitudinal waves and transverse waves, and has a complex equation. In 1985, shuey simplified the longitudinal wave reflection coefficient in Aki and Richards formulas, and the longitudinal wave reflection coefficient is decomposed into a plurality of units and is easy to apply on a computer. Many oil and gas reservoirs exist in the form of thin (mutual) layers, and the analysis of the seismic reflection coefficient of the thin (mutual) layers is always the difficulty of the prediction of the reflection seismic reservoirs.
Many researchers have developed thin (inter) layer seismic reflections based on Zeoppritz's simplified equations for single phase media (Kallweit et al, 1982, liu et al, 2003, lu et al, 2019). With the increased awareness of reservoirs, hydrocarbon reservoirs are increasingly recognized as typical pore media. The reflection coefficient of seismic waves should develop based on the pore medium theory. 1990. In the years, wanshanghai studied the propagation of elastic waves in two-phase media, and derived the equations for the reflection coefficient and transmission coefficient of elastic waves in two-phase media. In 1992, the reflection and transmission of sound waves at the interface of a two-phase medium were known from a series of studies by Johngxiao et al. In 2006, yong school and so on studied the two-phase medium AVO equation and parameter simplification thereof based on the Biot model, and the practicability of the two-phase medium AVA equation is enhanced through parameter simplification in the process.
In 2015, lihongxing and the like derived AVA equations of wave reflection and transmission at an interface containing a multiphase medium based on an improved BISQ multiphase medium theory. With the development of petrophysics, people realize that the frequency range of seismic waves is also in dispersion in size, and particularly, the dispersion is more obvious when the seismic waves propagate in oil and gas reservoirs. Chapman (2005), guo chile et al (2016) performed frequency-dependent AVO analyses based on a viscoelastic model. As the pore media continues to evolve theoretically, researchers have further realized that fluid flow at the mesoscopic scale is a significant cause of wave attenuation and dispersion. The seismic reflection wave amplitude reservoir prediction based on the pore medium theory can accurately and directly describe the fluid in the oil and gas reservoir.
In summary, the following problems mainly exist in the research of the existing method: the vast majority of reservoir seismic reflection coefficient analysis is based on a simplified Zeoppritz equation of a single-phase medium, and the characteristics of multiphase and pore media of the reservoir are not considered. Meanwhile, the frequency-dependent reflection coefficient analysis is also based on the equivalent medium thought, and the fluid characteristics of the oil and gas reservoir are difficult to accurately describe. Furthermore, no seismic reflection coefficient analysis method studies based on thin (inter) layers of porous media have been found.
Disclosure of Invention
The invention aims to solve the defects in the prior art and provides a frequency-dependent reflection coefficient analysis method for seismic waves of a thin reservoir or mutual reservoir model for correcting a Biot-jet flow, so as to analyze the change rule of the reflection coefficient of the seismic waves of seismic exploration reflected by the thin (mutual) reservoir along with the frequency and the incident angle.
In order to achieve the aim, the invention provides a method for analyzing the frequency-dependent reflection coefficient of seismic waves of a thin reservoir or interburden model of a modified Biot-jet flow, which comprises the following steps:
step 1: describing an actual oil and gas reservoir by using a modified Biot-jet flow pore medium model, and describing a top plate layer and a bottom plate layer of the oil and gas reservoir by using a common single-phase medium model;
step 2: solving and correcting a Biot-jet flow pore medium wave equation according to stress, strain and fluid pressure continuous conditions of a stratum interface to generate a seismic wave reflection coefficient of a thin (mutual) layer reservoir;
and step 3: and generating the reflection of the seismic waves at the reservoir interface along with the change of the frequency according to the thin (mutual) layer reservoir seismic wave reflection coefficient, and the relation between the transmission coefficient and the seismic wave incidence angle.
Further, in step 1, the modified Biot-jet pore media model used to describe the actual hydrocarbon reservoir is as follows:
Figure BDA0002451843900000021
Figure BDA0002451843900000022
Figure BDA0002451843900000023
ρ 1 =(1-φ)ρ s
ρ 2 =φρ f
ρ 12 =-ρ a
ρ 11 =ρ 112
ρ 22 =ρ 212
wherein the content of the first and second substances,
λ -pore medium framework elastic modulus;
μ — pore medium framework shear modulus;
α 0 -the pore elastic coefficient;
b-dissipation factor;
u-solid phase displacement vector;
u-liquid phase displacement vector
F-Biot flow coefficient;
s-jet flow coefficient;
ρ s -solid phase particle density;
ρ f -pore fluid density;
ρ a -solid-flow coupling density;
phi-porosity;
p-fluid pressure.
Further, in the step 2, the stress, strain and fluid pressure continuous conditions of the formation interface are as follows:
Figure BDA0002451843900000031
P 1 =P 2 | z=0 ,
Figure BDA0002451843900000032
P 2 =P 3 | z=h , />
Figure BDA0002451843900000033
wherein the content of the first and second substances,
σ -stress;
h-thin (inter) layer thickness;
subscript x, z-the direction of the vector components of displacement, stress, etc.;
upper/ subscript 1,2,3 — reservoir overburden media, reservoir pore media, reservoir underburden media.
Substituting the boundary conditions into a modified Biot-jet flow pore medium wave equation to generate a thin-layer reservoir seismic wave reflection coefficient solving equation: [ BD ] 1 -D 2 ]R T And =0, generating a solution equation of seismic wave reflection coefficients of the thin interbed reservoir:
Figure BDA0002451843900000041
wherein:
B=MWM -1
Figure BDA0002451843900000042
Figure BDA0002451843900000043
/>
Figure BDA0002451843900000044
Figure BDA0002451843900000051
Figure BDA0002451843900000052
E lm =A l +2N l cos 2 α lm +(Q 1 +R 1 )v lm +Q l
G lm =N l sin 2α lm
H lm =N l cos 2α lm
S lm =φ l (1-v lm )cosα lm
Figure BDA0002451843900000053
Figure BDA0002451843900000054
Figure BDA0002451843900000055
wherein the content of the first and second substances,
k is the complex number of waves;
alpha-angle;
a, N-elastic parameters of the framework without water drainage;
q, R-solid-fluid coupling elastic parameter;
v-fluid-solid displacement ratio;
l =1,2,3-reservoir overburden media, reservoir pore media, reservoir underburden media;
m =1,2,3-fast longitudinal wave, slow longitudinal wave, converted transverse wave;
Figure BDA0002451843900000056
-reflection wave coefficients;
Figure BDA0002451843900000057
-the transmitted wave coefficient;
p1, p2, s-fast longitudinal wave, slow longitudinal wave, and converted transverse wave.
The embodiment of the invention has the beneficial effects that:
the invention provides a novel reflection coefficient (amplitude) analysis method in reflection seismic exploration of a thin (mutual) layer oil and gas reservoir, and overcomes the problem that the traditional analysis does not strictly consider the current situation that the reservoir is a typical multiphase pore medium.
Aiming at the difficulty of analyzing the reflection seismic amplitude of the thin (mutual) layer reservoir, the method analyzes the reflection coefficient changing along with the frequency and the incident angle based on the multiphase pore medium model, has good applicability, can be used for analyzing the reflection seismic amplitude of the thin (mutual) layer reservoir changing along with the frequency and the incident angle, and provides method support for the seismic prediction of the oil and gas reservoir.
Drawings
FIG. 1 is a flow chart of a seismic wave frequency-dependent reflection coefficient analysis method for a modified Biot-jet flow thin (mutual) reservoir model according to an embodiment of the present invention;
FIG. 2 is a model explanatory diagram of a thin (inter) layer hydrocarbon reservoir in accordance with an embodiment of the present invention;
FIG. 3 is a flow chart of calculating seismic reflection coefficients of a thin (inter) reservoir in accordance with an embodiment of the present invention;
FIG. 4 is a reservoir parameter of an embodiment of the present invention;
FIG. 5 is a graph of reflection coefficient as a function of frequency and angle of incidence for (a) thick reservoirs, (b) thin reservoirs, according to an embodiment of the present invention;
FIG. 6 is a graph showing the variation of reflection coefficient with frequency and incident angle for different thicknesses of thin layers in accordance with an embodiment of the present invention;
FIG. 7 is a graph of the slow longitudinal wave velocity and inverse quality factor analysis of the model at different lacunarity;
FIG. 8 is a graph of the analysis of the slow longitudinal wave velocity and inverse quality factor of the model at different permeabilities;
FIG. 9 is a graph of the analysis of the slow compressional velocity and inverse quality factor of the model at different gas saturations.
Detailed Description
In order to make the purpose, technical solution and advantages of the embodiments of the present invention more apparent, the embodiments of the present invention are further described in detail below with reference to the accompanying drawings. The present invention is described in detail with reference to the accompanying drawings, which are incorporated in and constitute a part of this specification.
Referring to fig. 1 to 9, the present embodiment provides a method for seismic frequency-dependent reflection coefficient analysis of a modified Biot-jet thin reservoir or mutual reservoir model, as shown in fig. 1, the method including:
step S1: and describing an actual oil and gas reservoir by using a modified Biot-jet flow pore medium model, and describing a top plate layer and a bottom plate layer of the oil and gas reservoir by using a common single-phase medium model.
Fig. 2 is an explanatory diagram of a thin (inter) hydrocarbon reservoir model, which is described using a modified Biot-jet multiphase pore medium model:
Figure BDA0002451843900000061
Figure BDA0002451843900000062
Figure BDA0002451843900000071
ρ 1 =(1-φ)ρ s
ρ 2 =φρ f
ρ 12 =-ρ a
ρ 11 =ρ 112
ρ 22 =ρ 212
wherein the content of the first and second substances,
lambda is the pore medium framework elastic modulus;
μ -pore medium framework shear modulus;
α 0 -the pore elastic coefficient;
b-dissipation factor;
u-solid phase displacement vector;
u-liquid phase displacement vector
F-Biot flow coefficient;
s-jet flow coefficient;
ρ s -solid phase particle density;
ρ f -pore fluid density;
ρ a -solid-flow coupling density;
phi-porosity;
p-fluid pressure.
Step S2: and solving and correcting a Biot-jet flow pore medium wave equation according to the stress, strain and fluid pressure continuous conditions of the stratum interface to generate the seismic wave reflection coefficient of the thin (mutual) layer reservoir.
Fig. 3 is a process of calculating the seismic wave reflection coefficient of the thin (inter) reservoir in the invention, and in step S2, the continuous conditions of stress, strain and fluid pressure at the formation interface are as follows:
Figure BDA0002451843900000081
P 1 =P 2 | z=0 , />
Figure BDA0002451843900000082
P 2 =P 3 | z=h , />
Figure BDA0002451843900000083
wherein σ is stress;
h-thin (inter) layer thickness;
subscripts x, z-the direction of vector components of displacement, stress, etc.;
upper/ subscript 1,2,3 — reservoir overburden media, reservoir pore media, reservoir underburden media.
Substituting the boundary conditions into a modified Biot-jet flow pore medium wave equation to generate a thin-layer reservoir seismic wave reflection coefficient solving equation: [ BD ] 1 -D 2 ]R T And =0, generating a solution equation of seismic wave reflection coefficients of the thin interbed reservoir:
Figure BDA0002451843900000084
wherein:
B=MWM -1
Figure BDA0002451843900000085
Figure BDA0002451843900000086
Figure BDA0002451843900000091
Figure BDA0002451843900000092
Figure BDA0002451843900000093
E lm =A l +2N l cos 2 α lm +(Q 1 +R 1 )v lm +Q l
G lm =N l sin 2α lm
H lm =N l cos 2α lm
S lm =φl l (1-v lm )cosα lm
Figure BDA0002451843900000094
Figure BDA0002451843900000095
Figure BDA0002451843900000096
/>
wherein the content of the first and second substances,
k is the complex number of waves;
alpha-angle;
a, N-elastic parameters of the framework without water drainage;
q, R-solid-fluid coupling elastic parameter;
v-fluid-solid displacement ratio;
l =1,2,3-reservoir overburden media, reservoir pore media, reservoir underburden media;
m =1,2,3-fast longitudinal wave, slow longitudinal wave, converted transverse wave;
Figure BDA0002451843900000101
-reflection wave coefficients;
Figure BDA0002451843900000102
-a transmitted wave coefficient;
p1, p2, s-fast longitudinal wave, slow longitudinal wave, transverse wave.
The present invention is illustrated by the following specific examples.
The reservoir parameters are shown in figure 4.
The thickness of the thin (inter) layer reservoir is 5 m, which is less than the seismic wave wavelength.
Fig. 5 shows the variation of the reflection coefficient of the reflected seismic wave with frequency and incidence angle for a thin reservoir, and the variation of the reflection coefficient of the reflected seismic wave with frequency and incidence angle for a thick reservoir (the thickness is much larger than the seismic wavelength).
Fig. 6 shows the variation of reflection coefficient of the reflected seismic wave with frequency and incident angle when the thickness of the thin layer is different.
From the above, the invention provides a novel reflection coefficient (amplitude) analysis method in reflection seismic exploration of thin (mutual) layer oil and gas reservoirs, and overcomes the problem that the traditional analysis does not strictly consider the current situation that the reservoirs are typical multiphase pore media. The method is used for analyzing the reflection seismic amplitude variation of the thin (mutual) reservoir stratum along with the frequency and the incidence angle based on the multiphase pore medium model aiming at the difficulty of the reflection seismic amplitude analysis of the thin (mutual) reservoir stratum, has good applicability, can be used for analyzing the reflection seismic amplitude variation of the thin (mutual) reservoir stratum along with the frequency and the incidence angle, and provides method support for the seismic prediction of the oil and gas reservoir stratum.
The above-mentioned embodiments, objects, technical solutions and advantages of the present invention are further described in detail, it should be understood that the above-mentioned embodiments are only examples of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.

Claims (1)

1. A method for analyzing seismic wave frequency-dependent reflection coefficients of a thin reservoir or mutual reservoir model is characterized by comprising the following steps:
step 1: describing an actual oil and gas reservoir by using a modified Biot-jet flow pore medium model, and describing a top plate layer and a bottom plate layer of the oil and gas reservoir by using a common single-phase medium model;
step 2: solving and correcting a Biot-jet flow pore medium wave equation according to stress, strain and fluid pressure continuous conditions of a stratum interface to generate a seismic wave reflection coefficient of a thin reservoir or a mutual reservoir;
and step 3: generating the reflection of the seismic waves at the reservoir interface along with the change of the frequency according to the reflection coefficient of the seismic waves of the thin reservoir or the mutual reservoir, and the relation between the transmission coefficient and the incident angle of the seismic waves;
in step 1, the modified Biot-jet pore media model used to describe the actual hydrocarbon reservoir is as follows:
Figure DEST_PATH_IMAGE001
Figure DEST_PATH_IMAGE002
Figure DEST_PATH_IMAGE003
Figure DEST_PATH_IMAGE004
Figure DEST_PATH_IMAGE005
Figure DEST_PATH_IMAGE006
Figure DEST_PATH_IMAGE007
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE008
-pore medium backbone elastic modulus;
Figure DEST_PATH_IMAGE009
-pore medium framework shear modulus;
Figure DEST_PATH_IMAGE010
-the pore elastic coefficient;
Figure DEST_PATH_IMAGE011
-a dissipation factor;
Figure DEST_PATH_IMAGE012
-a solid phase displacement vector;
Figure DEST_PATH_IMAGE013
-liquid phase displacement vector
Figure DEST_PATH_IMAGE014
-Biot flow coefficient;
Figure DEST_PATH_IMAGE015
-the jet flow coefficient;
Figure DEST_PATH_IMAGE016
-solid phase particle density;
Figure DEST_PATH_IMAGE017
-pore fluid density;
Figure DEST_PATH_IMAGE018
-solid-flow coupling density;
Figure DEST_PATH_IMAGE019
-porosity;
Figure DEST_PATH_IMAGE020
-the fluid pressure; />
In the step 2, the continuous conditions of stress, strain and fluid pressure of the formation interface are as follows:
Figure DEST_PATH_IMAGE021
,/>
Figure DEST_PATH_IMAGE022
,/>
Figure DEST_PATH_IMAGE023
, />
Figure DEST_PATH_IMAGE024
, />
Figure DEST_PATH_IMAGE025
, />
Figure DEST_PATH_IMAGE026
,
Figure DEST_PATH_IMAGE027
,/>
Figure DEST_PATH_IMAGE028
,/>
Figure DEST_PATH_IMAGE029
, />
Figure DEST_PATH_IMAGE030
, />
Figure DEST_PATH_IMAGE031
, />
Figure DEST_PATH_IMAGE032
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE033
-stress;
Figure DEST_PATH_IMAGE034
-thin reservoir thickness or interbed thickness;
subscript
Figure DEST_PATH_IMAGE035
-direction of vector components of displacement, stress, etc.;
upper/lower 1,2,3-reservoir overburden media, reservoir pore media, reservoir underburden media;
substituting the boundary conditions into a modified Biot-jet flow pore medium wave equation to generate a thin-layer reservoir seismic wave reflection coefficient solving equation:
Figure DEST_PATH_IMAGE036
generating a solution equation of the seismic wave reflection coefficient of the thin interbed reservoir:
Figure DEST_PATH_IMAGE037
wherein:
Figure DEST_PATH_IMAGE038
Figure DEST_PATH_IMAGE039
Figure DEST_PATH_IMAGE040
Figure DEST_PATH_IMAGE041
/>
Figure DEST_PATH_IMAGE042
Figure DEST_PATH_IMAGE043
Figure DEST_PATH_IMAGE044
Figure DEST_PATH_IMAGE045
Figure DEST_PATH_IMAGE046
Figure DEST_PATH_IMAGE047
Figure DEST_PATH_IMAGE048
Figure DEST_PATH_IMAGE049
Figure DEST_PATH_IMAGE050
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE051
-a complex number of wave numbers;
Figure DEST_PATH_IMAGE052
-an angle;
Figure DEST_PATH_IMAGE053
,/>
Figure DEST_PATH_IMAGE054
-elasticity parameters of the non-draining carcass;
Figure DEST_PATH_IMAGE055
,/>
Figure DEST_PATH_IMAGE056
-a solid-flow coupling elasticity parameter;
Figure DEST_PATH_IMAGE057
-the flow-solid displacement ratio;
Figure DEST_PATH_IMAGE058
-reservoir overburden media, reservoir pore media, reservoir underburden media;
Figure DEST_PATH_IMAGE059
-fast longitudinal waves, slow longitudinal waves, transverse waves are converted;
Figure DEST_PATH_IMAGE060
-reflection wave coefficients;
Figure DEST_PATH_IMAGE061
-the transmitted wave coefficient;
Figure DEST_PATH_IMAGE062
, />
Figure DEST_PATH_IMAGE063
,/>
Figure DEST_PATH_IMAGE064
fast longitudinal waves, slow longitudinal waves, transverse waves. />
CN202010294977.4A 2020-04-15 2020-04-15 Method for analyzing seismic wave frequency-dependent reflection coefficient of thin reservoir or mutual reservoir model Active CN111562613B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010294977.4A CN111562613B (en) 2020-04-15 2020-04-15 Method for analyzing seismic wave frequency-dependent reflection coefficient of thin reservoir or mutual reservoir model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010294977.4A CN111562613B (en) 2020-04-15 2020-04-15 Method for analyzing seismic wave frequency-dependent reflection coefficient of thin reservoir or mutual reservoir model

Publications (2)

Publication Number Publication Date
CN111562613A CN111562613A (en) 2020-08-21
CN111562613B true CN111562613B (en) 2023-03-31

Family

ID=72070676

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010294977.4A Active CN111562613B (en) 2020-04-15 2020-04-15 Method for analyzing seismic wave frequency-dependent reflection coefficient of thin reservoir or mutual reservoir model

Country Status (1)

Country Link
CN (1) CN111562613B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113009571B (en) * 2021-02-18 2022-03-08 中国矿业大学(北京) Method for determining reflection coefficient and transmission coefficient of horizontal crack in two-phase medium
CN113253346B (en) * 2021-04-20 2022-04-15 中国石油大学(北京) Prediction method for horizontal fractures of thin interbed reservoir

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102508296A (en) * 2011-11-14 2012-06-20 中国石油天然气股份有限公司 Method and device for analyzing dispersion and attenuation of unsaturated double-porosity medium earthquake waves
CN104155693A (en) * 2014-08-29 2014-11-19 成都理工大学 Angle gather seismic response numerical computation method of reservoir fluid fluidity
CN104570084A (en) * 2013-10-29 2015-04-29 中国石油化工股份有限公司 Cross-scale seismic rock physical attenuation model and method for predicating attenuation and dispersion
CN109298443A (en) * 2017-07-25 2019-02-01 中国石油化工股份有限公司 Two-fluid injection stream model elastic wave evaluation method and computer readable storage medium

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102508296A (en) * 2011-11-14 2012-06-20 中国石油天然气股份有限公司 Method and device for analyzing dispersion and attenuation of unsaturated double-porosity medium earthquake waves
CN104570084A (en) * 2013-10-29 2015-04-29 中国石油化工股份有限公司 Cross-scale seismic rock physical attenuation model and method for predicating attenuation and dispersion
CN104155693A (en) * 2014-08-29 2014-11-19 成都理工大学 Angle gather seismic response numerical computation method of reservoir fluid fluidity
CN109298443A (en) * 2017-07-25 2019-02-01 中国石油化工股份有限公司 Two-fluid injection stream model elastic wave evaluation method and computer readable storage medium

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
nfluence of the squirt flow on reflection and refraction of elastic waves at a fluid/fluid一saturated poroelastic solid interface;Zhi-Wen Cui,Ke-Xie Wang;《International Journal of Engineering Science》;20031231;全文 *
基于改进BISQ模型的AVA方程;李红星,桂磊峰,陶春辉等;《石油学报》;20150930;第36卷(第9期);全文 *

Also Published As

Publication number Publication date
CN111562613A (en) 2020-08-21

Similar Documents

Publication Publication Date Title
Tang A unified theory for elastic wave propagation through porous media containing cracks—An extension of Biot’s poroelastic wave theory
Carcione et al. Bottom-simulating reflectors: Seismic velocities and AVO effects
CN111562613B (en) Method for analyzing seismic wave frequency-dependent reflection coefficient of thin reservoir or mutual reservoir model
CN111856560B (en) Natural gas hydrate reservoir information evaluation method and application thereof
CN112255680A (en) Seismic response simulation method for heterogeneous hydrate reservoir
Gajewski et al. Vector wavefields for weakly attenuating anisotropic media by the ray method
CN104007466B (en) The reservoir that a kind of no restriction from borehole data prestack inversion based on P-wave amplitude realizes and fluid prediction method
CN109031413B (en) vector wave reverse time migration system and method based on fluctuating submarine cable data
CN106249294A (en) A kind of reservoir detecting method of hydrocarbon
Guo et al. Seismic signatures of reservoir permeability based on the patchy-saturation model
Daley et al. Tube-wave suppression in single-well seismic acquisition
CN109407148A (en) One kind being based on the formation information prediction technique of " V " type AVO feature
Wapenaar Amplitude-variation-with-angle behavior of self-similar interfaces
Parra et al. Characterization of fractured low Q zones at the Buena Vista Hills reservoir, California
Goloshubin et al. Using frequency-dependent seismic attributes in imaging of a fractured reservoir zone
CN116819616B (en) Method for determining thickness of ultrathin high-quality shale reservoir
CN117316329B (en) Joint intelligent inversion method and system for natural gas hydrate saturation acoustic-electric logging
Qiu et al. Effect of viscosity on pseudo-Scholte wave propagation at liquid/porous medium interface
SUN A two-parameter model of elastic wave velocities in rocks and numerical avo modeling
Wang et al. DIRECT PREDICTION METHOD OF FRACTURING ABILITY IN SHALE FORMATIONS BASED ON PRE-STACK SEISMIC INVERSION
Li et al. The effect of attenuation on frequency-dependent AVO based on porous media theory
Li et al. Reflection and transmission of waves at a porous-medium/water interface
Santos et al. Analysis of mesoscopic loss effects in anisotropic poroelastic media using harmonic finite element simulations
Amundsen et al. The plane-wave primary reflection response from an impedance gradient interface
Xu et al. Seismic response analysis of sand-shale thin interbeds based on propagation matrix theory

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