CN109143351A - Prestack anisotropic character parameter inversion method and computer readable storage medium - Google Patents
Prestack anisotropic character parameter inversion method and computer readable storage medium Download PDFInfo
- Publication number
- CN109143351A CN109143351A CN201710507808.2A CN201710507808A CN109143351A CN 109143351 A CN109143351 A CN 109143351A CN 201710507808 A CN201710507808 A CN 201710507808A CN 109143351 A CN109143351 A CN 109143351A
- Authority
- CN
- China
- Prior art keywords
- velocity
- indicate
- longitudinal wave
- wave
- prestack
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000004590 computer program Methods 0.000 claims description 2
- 230000000149 penetrating effect Effects 0.000 claims 1
- 208000010392 Bone Fractures Diseases 0.000 abstract description 15
- 206010017076 Fracture Diseases 0.000 abstract description 15
- 230000015572 biosynthetic process Effects 0.000 abstract description 7
- 238000010586 diagram Methods 0.000 description 12
- 238000005755 formation reaction Methods 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 5
- 235000021185 dessert Nutrition 0.000 description 4
- 239000000284 extract Substances 0.000 description 4
- 238000000605 extraction Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 235000019738 Limestone Nutrition 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 239000003989 dielectric material Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000006028 limestone Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
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. for interpretation or for event detection
- 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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
-
- 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
- G01V2210/624—Reservoir parameters
-
- 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
- G01V2210/626—Physical property of subsurface with anisotropy
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Disclose a kind of prestack anisotropic character parameter inversion method and computer readable storage medium.This method may include: to obtain reflection coefficient based on prestack seismic gather；Prestack Simultaneous Inversion is carried out for earthquake data before superposition, obtains velocity of longitudinal wave, shear wave velocity and the density of top dielectric and layer dielectric respectively；Velocity of longitudinal wave, shear wave velocity and density based on reflection coefficient, top dielectric and layer dielectric carry out inverting according to Anisotropic parameters inversion formula, obtain the anisotropic parameters of top dielectric and layer dielectric.It converts on the basis of Ruger formula, potential is contacted with reflect between formation anisotropy parameter and fracture development.
Description
Technical field
The present invention relates to seismic data interpretation fields, more particularly, to a kind of prestack anisotropic character parametric inversion
Method and computer readable storage medium.
Background technique
Predict that crack is current relatively inexpensive, feasible side based on the azimuthal anisotropy AVO feature of P wave data
One of method.Theoretical basis of the Zoeppritz equation as AVO technology gives the anti-of plane wave propagation in HORIZONTAL LAYERED MEDIUM WITH HIGH ACCURACY
It penetrates and law of refraction, but due to its form complexity, it is difficult to determining physical significance is provided, accordingly, it is difficult to visual representations amplitude system
Relationship between several and medium parameter.Therefore, many experts and scholars carry out Zoeppritz equation approximate simplified, obtain
A variety of approximate formulas such as Aki, Shuey, Hilterman, but these are all premised on isotropic medium.With landing surface respectively to
Anisotropic problem is increasingly valued by people, and is discussed that influence of the anisotropic parameters to reflection coefficient becomes and inevitable is become
Gesture.
1977, Dely and Hron discussed interface reflection coefficient and transmission coefficient between the double-deck Method in Transverse Isotropic Medium
Accurate expression.1987, Wright analyzed in Method in Transverse Isotropic Medium anisotropy to reflecting and roll at interface
Penetrate generated influence.The same year, Banik discuss in Method in Transverse Isotropic Medium the physical significance of anisotropic parameters and each
Influence of the anisotropy parameter to reflection and transmission coefficients at interface.1993, it was anisotropy page that Kim, which establishes a upper layer,
Rock, the model that lower layer is isotropism gas sand, analyze influence of the anisotropy to P wave reflection coefficient at interface.
The approximate expression of longitudinal wave reflection coefficient in the Method in Transverse Isotropic Medium that Ruger was deduced weak anisotropy in 1996.
1998, Vavrycuk V and Psencik I gave in forward modeling in arbitrary anisotropic media longitudinal wave reflection coefficient with incidence angle, orientation
The approximate expression of angle variation.Previous research all shows that the anisotropy of medium can seriously affect isotropic medium model
AVO analysis under assuming that.
Facilitated to predict fracture strike according to azimuthal anisotropy AVO characteristic.1998, Ruger gave based on Ruger
The derivation of anisotropic parameters and fracture azimuth is intended in formula inverting, characterizes fracture development density with quasi- anisotropic parameters, from
And realize FRACTURE PREDICTION (R ü ger A.Variation of P-wave reflectivity with offset and
azimuth in anisotropic media[J].Geophysics,1998,63(3):935-947).2004, Gray gave
The application example for having gone out to carry out FRACTURE PREDICTION according to prestack orientation angle gathers data, provides and utilizes seismic location anisotropy
The thinking and techniqueflow in AVO feature prediction crack.2005, Zhu Zhaolin et al. have studied the longitudinal wave orientation of fracture medium respectively to
The inversion method of anisotropic AVO has derived the inversion formula based on Ruger formula and has been verified with model.2009, Xiao Pengfei
Et al. inquired into fracture spacing formula in TTI medium, further analyze influence of the fracture dip to inversion of Density.2010,
Downton analyzes the limitation of present azimuth AVO technology, and proposes a kind of new azimuthal anisotropy elastic inversion method
To overcome these limitations.In the inversion method that Ruger in 1998 is provided, inversion result is quasi- anisotropic parameters and side
Parallactic angle, rather than inverting obtains Thomsen anisotropic parameters.Therefore, it is necessary to which it is special to develop a kind of pre-stack seismic anisotropy
Parameter inversion method is levied, potential is contacted with reflect between formation anisotropy parameter and fracture development.
The information for being disclosed in background of invention part is merely intended to deepen the reason to general background technique of the invention
Solution, and it is known to those skilled in the art existing to be not construed as recognizing or imply that the information is constituted in any form
Technology.
Summary of the invention
The invention proposes a kind of prestack anisotropic character parameter inversion method and computer readable storage medium, energy
It is enough to be converted on the basis of Ruger formula, it potential is contacted with reflect between formation anisotropy parameter and fracture development.
According to an aspect of the invention, it is proposed that a kind of prestack anisotropic character parameter inversion method.The method can
To include: to obtain reflection coefficient based on prestack seismic gather；Prestack Simultaneous Inversion is carried out for earthquake data before superposition, is obtained respectively
Obtain velocity of longitudinal wave, shear wave velocity and the density of top dielectric and layer dielectric；Based on the reflection coefficient, the top dielectric with
Velocity of longitudinal wave, shear wave velocity and the density of layer dielectric carry out inverting according to Anisotropic parameters inversion formula, obtain on described
The anisotropic parameters of layer medium and layer dielectric.
Be preferably based on prestack seismic gather, obtain reflection coefficient include: the prestack seismic gather is divided into it is multiple
Orientation angle gathers；The multiple orientation angle gathers are converted into multiple incident angle gathers based on rate pattern；It is directed to along destination layer
Each incidence angle gathers extract reflection coefficient.
Preferably, described to carry out prestack Simultaneous Inversion for earthquake data before superposition, top dielectric is obtained respectively and lower layer is situated between
Velocity of longitudinal wave, shear wave velocity and the density of matter include: to carry out prestack Simultaneous Inversion to the earthquake data before superposition, obtain longitudinal wave speed
Spend body, shear wave velocity body and density body；Destination layer is divided into the top dielectric and the layer dielectric, is situated between on the upper layer
Window when selecting one in matter and the layer dielectric respectively, velocity of longitudinal wave body, shear wave velocity body corresponding to window when extracting described
With density body, and the root mean square of the corresponding velocity of longitudinal wave body, shear wave velocity body and density body is calculated separately, on described
Velocity of longitudinal wave, shear wave velocity and the density of layer medium and the layer dielectric.
Preferably, the Anisotropic parameters inversion formula are as follows:
Wherein, RpIndicate that reflected amplitude, i indicate incidence angle, φ indicates azimuth, and δ indicates longitudinal wave variation coefficient, Δ δ(V)
Indicate the difference of the longitudinal wave variation coefficient of upper and lower level medium, ε indicates longitudinal wave anisotropic parameters, Δ ε(V)Indicate upper and lower level medium
The difference of longitudinal wave anisotropic parameters, γ indicate that shear wave anisotropic parameters, Δ γ indicate the shear wave anisotropy of upper and lower level medium
The difference of parameter,Indicate that velocity of longitudinal wave, Δ α indicate the difference of the velocity of longitudinal wave of upper and lower level medium,Indicate shear wave velocity,It indicates
Wave impedance when longitudinal wave vertical incidence, Δ Z indicate difference of impedance when longitudinal wave vertical incidence,Indicate the tangential modulus of shear wave,
Δ G indicates that the tangential modulus of shear wave is poor.
Preferably, further includes: the anisotropic parameters based on the top dielectric and layer dielectric calculate the upper layer and are situated between
The difference of the anisotropic parameters of matter and layer dielectric.
According to another aspect of the invention, it is proposed that a kind of computer readable storage medium, is stored thereon with computer journey
Sequence, wherein performed the steps of when described program is executed by processor based on prestack seismic gather, obtain reflection coefficient；For
Earthquake data before superposition carries out prestack Simultaneous Inversion, obtain respectively the velocity of longitudinal wave of top dielectric and layer dielectric, shear wave velocity with
Density；Based on the reflection coefficient, the velocity of longitudinal wave of the top dielectric and layer dielectric, shear wave velocity and density, according to each
Anisotropy parametric inversion formula carries out inverting, obtains the anisotropic parameters of the top dielectric and layer dielectric.
Be preferably based on prestack seismic gather, obtain reflection coefficient include: the prestack seismic gather is divided into it is multiple
Orientation angle gathers；The multiple orientation angle gathers are converted into multiple incident angle gathers based on rate pattern；It is directed to along destination layer
Each incidence angle gathers extract reflection coefficient.
Preferably, described to carry out prestack Simultaneous Inversion for earthquake data before superposition, top dielectric is obtained respectively and lower layer is situated between
Velocity of longitudinal wave, shear wave velocity and the density of matter include: to carry out prestack Simultaneous Inversion to the earthquake data before superposition, obtain longitudinal wave speed
Spend body, shear wave velocity body and density body；Destination layer is divided into the top dielectric and the layer dielectric, is situated between on the upper layer
Window when selecting one in matter and the layer dielectric respectively, velocity of longitudinal wave body, shear wave velocity body corresponding to window when extracting described
With density body, and the root mean square of the corresponding velocity of longitudinal wave body, shear wave velocity body and density body is calculated separately, on described
Velocity of longitudinal wave, shear wave velocity and the density of layer medium and the layer dielectric.
Preferably, the Anisotropic parameters inversion formula are as follows:
Wherein, RpIndicate that reflected amplitude, i indicate incidence angle, φ indicates azimuth, and δ indicates longitudinal wave variation coefficient, Δ δ(V)
Indicate the difference of the longitudinal wave variation coefficient of upper and lower level medium, ε indicates longitudinal wave anisotropic parameters, Δ ε(V)Indicate upper and lower level medium
The difference of longitudinal wave anisotropic parameters, γ indicate that shear wave anisotropic parameters, Δ γ indicate the shear wave anisotropy of upper and lower level medium
The difference of parameter,Indicate that velocity of longitudinal wave, Δ α indicate the difference of the velocity of longitudinal wave of upper and lower level medium,Indicate shear wave velocity,It indicates
Wave impedance when longitudinal wave vertical incidence, Δ Z indicate difference of impedance when longitudinal wave vertical incidence,Indicate the tangential modulus of shear wave,
Δ G indicates that the tangential modulus of shear wave is poor.
The beneficial effects of the present invention are: using the present invention can steadily inverting obtain the variation of anisotropic parameters
Rate, and the anisotropic parameters distribution of reservoir is further obtained, the power of anisotropic parameters corresponds to the fracture development of reservoir
Intensity, anisotropy is stronger, and crack is more developed, the intensity of anisotropic parameters is stronger.Therefore, for unconventional equal slit formations oil
Gas reservoir can qualitatively indicate the advantageous dessert area of reservoir using anisotropic parameters.
Methods and apparatus of the present invention has other characteristics and advantages, these characteristics and advantages are attached from what is be incorporated herein
It will be apparent in figure and subsequent specific embodiment, or will be in the attached drawing being incorporated herein and subsequent specific reality
It applies in mode and is stated in detail, the drawings and the detailed description together serve to explain specific principles of the invention.
Detailed description of the invention
Exemplary embodiment of the present is described in more detail in conjunction with the accompanying drawings, of the invention is above-mentioned and other
Purpose, feature and advantage will be apparent, wherein in exemplary embodiments of the present invention, identical reference label is usual
Represent same parts.
Fig. 1 shows the flow chart of the step of prestack anisotropic character parameter inversion method according to the present invention.
Fig. 2 shows the schematic diagrames of prestack trace gather according to an embodiment of the invention.
Fig. 3 a, Fig. 3 b, Fig. 3 c, Fig. 3 d, Fig. 3 e respectively illustrate 0-50 degree according to an embodiment of the invention, 20-
70 degree, the schematic diagram of the orientation trace gather of 60-110 degree, 80-130 degree and 130-180 degree.
Fig. 4 shows the schematic diagram of rate pattern according to an embodiment of the invention.
Fig. 5 a, Fig. 5 b, Fig. 5 c, Fig. 5 d, Fig. 5 e respectively illustrate 0-50 degree according to an embodiment of the invention, 20-
70 degree, the schematic diagrames of the incident angle gathers of the orientation trace gather of 60-110 degree, 80-130 degree and 130-180 degree.
Fig. 6 shows the schematic diagram of reflected amplitude slice according to an embodiment of the invention.
Fig. 7 a, Fig. 7 b respectively illustrate the signal of the density slice of upper and lower level medium according to an embodiment of the invention
Figure.
Fig. 8 a, Fig. 8 b respectively illustrate the signal of the velocity of longitudinal wave of upper and lower level medium according to an embodiment of the invention
Figure.
Fig. 9 a, Fig. 9 b, Fig. 9 c respectively illustrate anisotropic parameters delta attribute according to an embodiment of the invention and cut
The schematic diagram that piece, anisotropic parameters ε attribute slice, anisotropic parameters γ attribute are sliced.
Specific embodiment
The present invention will be described in more detail below with reference to accompanying drawings.Although showing the preferred embodiment of the present invention in attached drawing,
However, it is to be appreciated that may be realized in various forms the present invention and should not be limited by the embodiments set forth herein.On the contrary, providing
These embodiments are of the invention more thorough and complete in order to make, and can will fully convey the scope of the invention to ability
The technical staff in domain.
Fig. 1 shows the flow chart of the step of prestack anisotropic character parameter inversion method according to the present invention.
In this embodiment, prestack anisotropic character parameter inversion method according to the present invention may include:
Step 101, it is based on prestack seismic gather, obtains reflection coefficient；In one example, it is based on prestack seismic gather,
Obtaining reflection coefficient includes: that prestack seismic gather is divided into multiple orientation angle gathers, and multiple orientation angle gathers are overlapped
Obtain channel set；Using root mean sequare velocity as rate pattern, angle analysis is done to channel set, obtains ranges of incidence angles, by incidence angle
Multiple small patches are averagely divided into, superposition is done respectively and obtains multiple incident angle gathers；Along destination layer for each incident angle gathers
Reflection coefficient is extracted, those skilled in the art, which can divide small patches as the case may be and be superimposed, obtains incident angle gathers.
Specifically, include incidence angle and azimuthal reflection coefficient to obtain, need prestack trace gather being divided into difference
Narrow orientation trace gather data incident angle gathers data are converted to based on rate pattern, finally along purpose to different direction angle gathers
Layer extracts obtained reflection coefficient and just contains incidence angle, azimuth information.
Step 102, prestack Simultaneous Inversion is carried out for earthquake data before superposition, obtains top dielectric and layer dielectric respectively
Velocity of longitudinal wave, shear wave velocity and density；In one example, prestack Simultaneous Inversion is carried out for earthquake data before superposition, obtained respectively
Velocity of longitudinal wave, shear wave velocity and the density for obtaining top dielectric and layer dielectric include: to carry out prestack simultaneously to earthquake data before superposition
Inverting obtains velocity of longitudinal wave body, shear wave velocity body and density body；It is boundary with the destination layer position that significant change occurs for wave impedance,
Such as the boundary of coal seam and layer of sand, the boundary of sandstone and limestone cause wave impedance that significant change occurs due to medium difference,
Boundary is the above are top dielectric, and the following are layer dielectrics for boundary, window when selecting one respectively in top dielectric and layer dielectric,
In conjunction with the practical geological condition in work area, choose include several groups of adjacent seismic events when window, those skilled in the art can be with
The when window comprising seismic event is chosen as the case may be.When extraction velocity of longitudinal wave body, shear wave velocity body corresponding to window with
Density body, and the root mean square of corresponding velocity of longitudinal wave body, shear wave velocity body and density body is calculated separately, as top dielectric under
Velocity of longitudinal wave, shear wave velocity and the density of layer medium.
Step 103, velocity of longitudinal wave, shear wave velocity and density based on reflection coefficient, top dielectric and layer dielectric, according to
Anisotropic parameters inversion formula carries out inverting, obtains the anisotropic parameters of top dielectric and layer dielectric.
In one example, Anisotropic parameters inversion formula are as follows:
Wherein, RpIndicate that reflected amplitude, i indicate incidence angle, φ indicates azimuth, and δ indicates longitudinal wave variation coefficient, Δ δ(V)
Indicate the difference of the longitudinal wave variation coefficient of upper and lower level medium, ε indicates longitudinal wave anisotropic parameters, Δ ε(V)Indicate upper and lower level medium
The difference of longitudinal wave anisotropic parameters, γ indicate that shear wave anisotropic parameters, Δ γ indicate the shear wave anisotropy of upper and lower level medium
The difference of parameter,Indicate that velocity of longitudinal wave, Δ α indicate the difference of the velocity of longitudinal wave of upper and lower level medium,Indicate shear wave velocity,It indicates
Wave impedance when longitudinal wave vertical incidence, Δ Z indicate difference of impedance when longitudinal wave vertical incidence,Indicate the tangential modulus of shear wave,
Δ G indicates that the tangential modulus of shear wave is poor.
Specifically, Anisotropic parameters inversion formula be formula (1), the road N incidence angle gathers earthquake record reflection coefficient with
Linear relationship between incidence angle, azimuth can use following systems of linear equations, i.e. formula (2) indicates:
Ax=b (2)
Wherein,
In formula (2), anisotropic parameters to be asked are indicated with matrix x, and coefficient matrices A is each incident angle gathers
Azimuth, incidence angle information, matrix b include the constant matrices of reflectance information.It, can be in the hope of by solution formula (2)
Anisotropic parameters matrix x.
In actual treatment, formula (2) is mostly overdetermined equation, can be solved with least square method, is solved as x=(ATA)-1ATB,
Work as ATWhen A is ill-condition matrix, the generalized inverse of singular value decomposition method solution matrix A can be used to seek vector x.By reflection coefficient, on
Layer medium and the velocity of longitudinal wave of layer dielectric, shear wave velocity and density substitute into formula (2), seek top dielectric and layer dielectric
Anisotropic parameters.
In one example, further includes: the anisotropic parameters based on top dielectric and layer dielectric, and then calculate upper layer
The difference of medium and the anisotropic parameters of layer dielectric.
Using the present invention can steadily inverting obtain the change rate of anisotropic parameters, and further obtain each of reservoir
Anisotropy parameter distribution, the power of anisotropic parameters correspond to the fracture development intensity of reservoir, and anisotropy is stronger, and crack is got over
Development, the intensity of anisotropic parameters are stronger.Therefore, for unconventional equal slit formations oil and gas reservoir, anisotropic parameters are utilized
It can qualitatively indicate the advantageous dessert area of reservoir.
Using example
A concrete application example is given below in the scheme and its effect of the embodiment of the present invention for ease of understanding.This field
It should be understood to the one skilled in the art that the example is only for the purposes of understanding the present invention, any detail is not intended to be limited in any way
The system present invention.
Fig. 2 shows the schematic diagrames of prestack trace gather according to an embodiment of the invention.
Fig. 3 a, Fig. 3 b, Fig. 3 c, Fig. 3 d, Fig. 3 e respectively illustrate 0-50 degree according to an embodiment of the invention, 20-
70 degree, the schematic diagram of the orientation trace gather of 60-110 degree, 80-130 degree and 130-180 degree.
Fig. 4 shows the schematic diagram of rate pattern according to an embodiment of the invention.
Fig. 5 a, Fig. 5 b, Fig. 5 c, Fig. 5 d, Fig. 5 e respectively illustrate 0-50 degree according to an embodiment of the invention, 20-
70 degree, the schematic diagrames of the incident angle gathers of the orientation trace gather of 60-110 degree, 80-130 degree and 130-180 degree.
Fig. 6 shows the schematic diagram of reflected amplitude slice according to an embodiment of the invention.
Fig. 7 a, Fig. 7 b respectively illustrate the signal of the density slice of upper and lower level medium according to an embodiment of the invention
Figure.
Fig. 8 a, Fig. 8 b respectively illustrate the signal of the velocity of longitudinal wave of upper and lower level medium according to an embodiment of the invention
Figure.
Fig. 9 a, Fig. 9 b, Fig. 9 c respectively illustrate anisotropic parameters delta attribute according to an embodiment of the invention and cut
The schematic diagram that piece, anisotropic parameters ε attribute slice, anisotropic parameters γ attribute are sliced.
As shown in Fig. 2, prestack trace gather to be divided into different narrow orientation trace gather data, 0-50 degree, 20-70 have been extracted respectively
The orientation trace gather of degree, 60-110 degree, 80-130 degree and 130-180 degree, as shown in Fig. 3 a, Fig. 3 b, Fig. 3 c, Fig. 3 d, Fig. 3 e, to not
Common-azimuth trace gather is converted to incident angle gathers data based on rate pattern as shown in Figure 4, such as Fig. 5 a, Fig. 5 b, Fig. 5 c, figure
Shown in 5d, Fig. 5 e, the seismic reflection lineups of purpose position layer by layer change with azimuthal different, then to different direction trace gather
Incident angle gathers data are extracted to obtain reflected amplitude slice, as shown in fig. 6, being to obtain to contain incidence angle, azimuth information
Reflection coefficient.
Prestack Simultaneous Inversion is carried out to earthquake data before superposition, obtains velocity of longitudinal wave body, shear wave velocity body and density body；By mesh
Mark layer is divided into top dielectric and layer dielectric, the window in top dielectric and window when selecting one respectively in layer dielectric, extraction
Corresponding velocity of longitudinal wave body, shear wave velocity body and density body, and calculate separately corresponding velocity of longitudinal wave body, shear wave velocity body with
The root mean square of density body, as the velocity of longitudinal wave of top dielectric and layer dielectric, shear wave velocity and density, such as Fig. 7 a, Fig. 7 b, figure
Shown in 8a, Fig. 8 b, it can be seen that the physical difference between upper and lower level is obvious.
The velocity of longitudinal wave, shear wave velocity and density of reflection coefficient, top dielectric and layer dielectric are substituted into formula (2), asked
The anisotropic parameters for taking top dielectric and layer dielectric, as shown in Fig. 9 a, Fig. 9 b, Fig. 9 c, anisotropic parameters δ is that longitudinal wave becomes
Different coefficient indicates the speed degree that longitudinal wave changes in vertical direction anisotropy, and anisotropic parameters ε is longitudinal wave anisotropy,
It is the parameter for measuring quasi-longitudinal wave strength of anisotropy, anisotropic parameters γ is shear wave anisotropy, measures quasi-transverse wave respectively to different
The parameter of property or shear wave splitting intensity.It can be found that the distribution of attribute abnormal value is substantially consistent with tomography trend, it may be speculated that
The anisotropic parameters difference larger part of upper and lower bed boundary, may crannied development.
In conclusion using the present invention can steadily inverting obtain the change rate of anisotropic parameters, and further
Anisotropic parameters to reservoir are distributed, and the power of anisotropic parameters corresponds to the fracture development intensity of reservoir, anisotropy
Stronger, crack is more developed, the intensity of anisotropic parameters is stronger.Therefore, for unconventional equal slit formations oil and gas reservoir, using each
Anisotropy parameter can qualitatively indicate the advantageous dessert area of reservoir.
It will be understood by those skilled in the art that above to the purpose of the description of the embodiment of the present invention only for illustratively saying
The beneficial effect of bright the embodiment of the present invention is not intended to limit embodiments of the invention to given any example.
According to an embodiment of the invention, providing a kind of computer readable storage medium, it is stored thereon with computer program,
Wherein, it is performed the steps of when program is executed by processor based on prestack seismic gather, obtains reflection coefficient；For prestack
It shakes data and carries out prestack Simultaneous Inversion, obtain velocity of longitudinal wave, shear wave velocity and the density of top dielectric and layer dielectric respectively；Base
In the velocity of longitudinal wave, shear wave velocity and density of reflection coefficient, top dielectric and layer dielectric, according to Anisotropic parameters inversion public affairs
Formula carries out inverting, obtains the anisotropic parameters of top dielectric and layer dielectric.
In one example, be based on prestack seismic gather, obtain reflection coefficient include: prestack seismic gather is divided into it is more
A orientation angle gathers；Multiple orientation angle gathers are converted into multiple incident angle gathers based on rate pattern；Along destination layer for every
A incidence angle gathers extract reflection coefficient.
In one example, prestack Simultaneous Inversion is carried out for earthquake data before superposition, obtains top dielectric and lower layer respectively
Velocity of longitudinal wave, shear wave velocity and the density of medium include: to carry out prestack Simultaneous Inversion to earthquake data before superposition, obtain velocity of longitudinal wave
Body, shear wave velocity body and density body；Destination layer is divided into top dielectric and layer dielectric, in top dielectric and layer dielectric
Window when selecting one respectively, velocity of longitudinal wave body, shear wave velocity body and density body corresponding to window when extraction, and calculate separately correspondence
Velocity of longitudinal wave body, shear wave velocity body and density body root mean square, as the velocity of longitudinal wave of top dielectric and layer dielectric, shear wave
Speed and density.
In one example, Anisotropic parameters inversion formula are as follows:
Wherein, RpIndicate that reflected amplitude, i indicate incidence angle, φ indicates azimuth, and δ indicates longitudinal wave variation coefficient, Δ δ(V)
Indicate the difference of the longitudinal wave variation coefficient of upper and lower level medium, ε indicates longitudinal wave anisotropic parameters, Δ ε(V)Indicate upper and lower level medium
The difference of longitudinal wave anisotropic parameters, γ indicate that shear wave anisotropic parameters, Δ γ indicate the shear wave anisotropy of upper and lower level medium
The difference of parameter,Indicate that velocity of longitudinal wave, Δ α indicate the difference of the velocity of longitudinal wave of upper and lower level medium,Indicate shear wave velocity,It indicates
Wave impedance when longitudinal wave vertical incidence, Δ Z indicate difference of impedance when longitudinal wave vertical incidence,Indicate the tangential modulus of shear wave,
Δ G indicates that the tangential modulus of shear wave is poor.
In one example, further includes: the anisotropic parameters based on top dielectric and layer dielectric, and then calculate upper layer
The difference of medium and the anisotropic parameters of layer dielectric.
Using the present invention can steadily inverting obtain the change rate of anisotropic parameters, and further obtain each of reservoir
Anisotropy parameter distribution, the power of anisotropic parameters correspond to the fracture development intensity of reservoir, and anisotropy is stronger, and crack is got over
Development, the intensity of anisotropic parameters are stronger.Therefore, for unconventional equal slit formations oil and gas reservoir, anisotropic parameters are utilized
It can qualitatively indicate the advantageous dessert area of reservoir.
Various embodiments of the present invention are described above, above description is exemplary, and non-exclusive, and
It is not limited to disclosed each embodiment.Without departing from the scope and spirit of illustrated each embodiment, for this skill
Many modifications and changes are obvious for the those of ordinary skill in art field.The selection of term used herein, purport
In the principle, practical application or improvement to the technology in market for best explaining each embodiment, or make the art
Other those of ordinary skill can understand each embodiment disclosed herein.
Claims (10)
1. a kind of prestack anisotropic character parameter inversion method, comprising:
Based on prestack seismic gather, reflection coefficient is obtained；
Prestack Simultaneous Inversion is carried out for earthquake data before superposition, obtains velocity of longitudinal wave, the cross of top dielectric and layer dielectric respectively
Wave velocity and density；
Based on the reflection coefficient, the velocity of longitudinal wave of the top dielectric and layer dielectric, shear wave velocity and density, according to it is each to
Anisotropic parameter inversion formula carries out inverting, obtains the anisotropic parameters of the top dielectric and layer dielectric.
2. prestack anisotropic character parameter inversion method according to claim 1, wherein it is based on prestack seismic gather,
Obtaining reflection coefficient includes:
The prestack seismic gather is divided into multiple orientation angle gathers；
The multiple orientation angle gathers are converted into multiple incident angle gathers based on rate pattern；
Reflection coefficient is extracted for each incident angle gathers along destination layer.
3. prestack anisotropic character parameter inversion method according to claim 1, wherein described to be directed to pre-stack seismic number
According to prestack Simultaneous Inversion is carried out, top dielectric is obtained respectively and velocity of longitudinal wave, shear wave velocity and the density of layer dielectric include:
Prestack Simultaneous Inversion is carried out to the earthquake data before superposition, obtains velocity of longitudinal wave body, shear wave velocity body and density body；
Destination layer is divided into the top dielectric and the layer dielectric, is divided in the top dielectric and the layer dielectric
Window when Xuan Ding not be one, velocity of longitudinal wave body, shear wave velocity body corresponding to window and density body when extracting described, and calculate separately institute
The root mean square for stating corresponding velocity of longitudinal wave body, shear wave velocity body and density body, as the top dielectric and the layer dielectric
Velocity of longitudinal wave, shear wave velocity and density.
4. prestack anisotropic character parameter inversion method according to claim 1, wherein the anisotropic parameters are anti-
Drill formula are as follows:
Wherein, RpIndicate that reflected amplitude, i indicate incidence angle, φ indicates azimuth, and δ indicates longitudinal wave variation coefficient, Δ δ(V)In expression
The difference of the longitudinal wave variation coefficient of layer dielectric, ε indicate longitudinal wave anisotropic parameters ,+ε(V)Indicate upper and lower level medium longitudinal wave respectively to
The difference of Anisotropic parameter, γ indicate that shear wave anisotropic parameters, Δ γ indicate the difference of the shear wave anisotropic parameters of upper and lower level medium,Indicate that velocity of longitudinal wave, Δ α indicate the difference of the velocity of longitudinal wave of upper and lower level medium,Indicate shear wave velocity,Indicate that longitudinal wave vertically enters
Wave impedance when penetrating, Δ Z indicate difference of impedance when longitudinal wave vertical incidence,Indicate that the tangential modulus of shear wave, Δ G indicate shear wave
Tangential modulus it is poor.
5. prestack anisotropic character parameter inversion method according to claim 1, wherein further include: on described
The anisotropic parameters of layer medium and layer dielectric, calculate the difference of the anisotropic parameters of the top dielectric and layer dielectric.
6. a kind of computer readable storage medium, is stored thereon with computer program, wherein when described program is executed by processor
It performs the steps of
Based on prestack seismic gather, reflection coefficient is obtained；
Prestack Simultaneous Inversion is carried out for earthquake data before superposition, obtains velocity of longitudinal wave, the cross of top dielectric and layer dielectric respectively
Wave velocity and density；
Based on the reflection coefficient, the velocity of longitudinal wave of the top dielectric and layer dielectric, shear wave velocity and density, according to it is each to
Anisotropic parameter inversion formula carries out inverting, obtains the anisotropic parameters of the top dielectric and layer dielectric.
7. computer readable storage medium according to claim 6, wherein be based on prestack seismic gather, obtain reflection system
Number includes:
The prestack seismic gather is divided into multiple orientation angle gathers；
The multiple orientation angle gathers are converted into multiple incident angle gathers based on rate pattern；
Reflection coefficient is extracted for each incident angle gathers along destination layer.
8. computer readable storage medium according to claim 6, wherein described to carry out prestack for earthquake data before superposition
Simultaneous Inversion, obtains top dielectric respectively and velocity of longitudinal wave, shear wave velocity and the density of layer dielectric include:
Prestack Simultaneous Inversion is carried out to the earthquake data before superposition, obtains velocity of longitudinal wave body, shear wave velocity body and density body；
Destination layer is divided into the top dielectric and the layer dielectric, is divided in the top dielectric and the layer dielectric
Window when Xuan Ding not be one, velocity of longitudinal wave body, shear wave velocity body corresponding to window and density body when extracting described, and calculate separately institute
The root mean square for stating corresponding velocity of longitudinal wave body, shear wave velocity body and density body, as the top dielectric and the layer dielectric
Velocity of longitudinal wave, shear wave velocity and density.
9. computer readable storage medium according to claim 6, wherein the Anisotropic parameters inversion formula are as follows:
Wherein, RpIndicate that reflected amplitude, i indicate incidence angle, φ indicates azimuth, and δ indicates longitudinal wave variation coefficient, Δ δ(V)In expression
The difference of the longitudinal wave variation coefficient of layer dielectric, ε indicate longitudinal wave anisotropic parameters, Δ ε(V)Indicate that the longitudinal wave of upper and lower level medium is each
The difference of anisotropy parameter, γ indicate shear wave anisotropic parameters, Δ γ indicate upper and lower level medium shear wave anisotropic parameters it
Difference,Indicate that velocity of longitudinal wave, Δ α indicate the difference of the velocity of longitudinal wave of upper and lower level medium,Indicate shear wave velocity,Indicate that longitudinal wave hangs down
Wave impedance when straight incident, Δ Z indicate difference of impedance when longitudinal wave vertical incidence,Indicate that the tangential modulus of shear wave, Δ G indicate
The tangential modulus of shear wave is poor.
10. computer readable storage medium according to claim 6, wherein further include: based on the top dielectric under
The anisotropic parameters of layer medium, calculate the difference of the anisotropic parameters of the top dielectric and layer dielectric.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710507808.2A CN109143351B (en) | 2017-06-28 | 2017-06-28 | Pre-stack anisotropy characteristic parameter inversion method and computer readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710507808.2A CN109143351B (en) | 2017-06-28 | 2017-06-28 | Pre-stack anisotropy characteristic parameter inversion method and computer readable storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109143351A true CN109143351A (en) | 2019-01-04 |
| CN109143351B CN109143351B (en) | 2020-07-24 |
Family
ID=64803068
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710507808.2A Active CN109143351B (en) | 2017-06-28 | 2017-06-28 | Pre-stack anisotropy characteristic parameter inversion method and computer readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109143351B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111551990A (en) * | 2020-04-07 | 2020-08-18 | 西安科技大学 | Method and system for calculating seismic wave reflection coefficient of HTI (HTI) coal seam |
| CN112379437A (en) * | 2020-11-02 | 2021-02-19 | 中国石油天然气集团有限公司 | Shale reservoir anisotropic parameter solving method and device |
| CN113568045A (en) * | 2021-07-26 | 2021-10-29 | 中南大学 | Inversion method of TTI medium crack inclination angle |
| CN114185092A (en) * | 2020-09-15 | 2022-03-15 | 中国石油化工股份有限公司 | Method and device for evaluating development degree of horizontal seam of reservoir, electronic equipment and medium |
| CN116299706A (en) * | 2023-03-23 | 2023-06-23 | 福瑞升(成都)科技有限公司 | Mixed impedance sand body identification method based on pre-stack quasi-transverse wave reflectivity attribute |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104407378A (en) * | 2014-11-25 | 2015-03-11 | 中国石油天然气股份有限公司 | Method and device for inversing anisotropy parameters |
| US20160054464A1 (en) * | 2013-04-02 | 2016-02-25 | Halliburton Energy Services, Inc. | Anisotropy analysis using direct and reflected arrivals in seismic survey data |
| CN105629303A (en) * | 2015-12-28 | 2016-06-01 | 中国石油大学(北京) | Prestack crack quantitative forecast method and system based on rock physics |
| WO2017048285A1 (en) * | 2015-09-18 | 2017-03-23 | Halliburton Energy Services, Inc. | Global inversion based estimation of anisotropy parameters for orthorhombic media |
-
2017
- 2017-06-28 CN CN201710507808.2A patent/CN109143351B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160054464A1 (en) * | 2013-04-02 | 2016-02-25 | Halliburton Energy Services, Inc. | Anisotropy analysis using direct and reflected arrivals in seismic survey data |
| CN104407378A (en) * | 2014-11-25 | 2015-03-11 | 中国石油天然气股份有限公司 | Method and device for inversing anisotropy parameters |
| WO2017048285A1 (en) * | 2015-09-18 | 2017-03-23 | Halliburton Energy Services, Inc. | Global inversion based estimation of anisotropy parameters for orthorhombic media |
| CN105629303A (en) * | 2015-12-28 | 2016-06-01 | 中国石油大学(北京) | Prestack crack quantitative forecast method and system based on rock physics |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111551990A (en) * | 2020-04-07 | 2020-08-18 | 西安科技大学 | Method and system for calculating seismic wave reflection coefficient of HTI (HTI) coal seam |
| CN111551990B (en) * | 2020-04-07 | 2023-05-23 | 西安科技大学 | Acquisition system for HTI type coal seam seismic wave reflection coefficient |
| CN114185092A (en) * | 2020-09-15 | 2022-03-15 | 中国石油化工股份有限公司 | Method and device for evaluating development degree of horizontal seam of reservoir, electronic equipment and medium |
| CN114185092B (en) * | 2020-09-15 | 2024-05-07 | 中国石油化工股份有限公司 | Reservoir horizontal seam development degree evaluation method and device, electronic equipment and medium |
| CN112379437A (en) * | 2020-11-02 | 2021-02-19 | 中国石油天然气集团有限公司 | Shale reservoir anisotropic parameter solving method and device |
| CN112379437B (en) * | 2020-11-02 | 2024-03-26 | 中国石油天然气集团有限公司 | Shale reservoir anisotropy parameter solving method and device |
| CN113568045A (en) * | 2021-07-26 | 2021-10-29 | 中南大学 | Inversion method of TTI medium crack inclination angle |
| CN116299706A (en) * | 2023-03-23 | 2023-06-23 | 福瑞升(成都)科技有限公司 | Mixed impedance sand body identification method based on pre-stack quasi-transverse wave reflectivity attribute |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109143351B (en) | 2020-07-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109143351A (en) | Prestack anisotropic character parameter inversion method and computer readable storage medium | |
| CN103733089B (en) | For including the system and method that the underground of uncertainty estimation characterizes | |
| CN103257363B (en) | A kind of method of fracture dip in Underground fracture-type reservoir | |
| Shen et al. | Azimuthal offset-dependent attributes applied to fracture detection in a carbonate reservoir | |
| Refunjol et al. | Inversion and attribute-assisted hydraulically induced microseismic fracture characterization in the North Texas Barnett Shale | |
| CN104005760A (en) | Azimuthal anisotropic elastic impedance based crack detection method | |
| US20140343859A1 (en) | Time-lapse 4d scattering for imaging hydraulically induced fractures | |
| Yang et al. | Integrated application of 3D seismic and microseismic data in the development of tight gas reservoirs | |
| CN104297800B (en) | A kind of from phased prestack inversion method | |
| Norton et al. | Surface seismic to microseismic: An integrated case study from exploration to completion in the Montney shale of NE British Columbia, Canada | |
| CN106501872B (en) | A kind of calculation method and device of fracture reservoir ground stress characteristics | |
| Zhang et al. | Velocity modeling and inversion techniques for locating microseismic events in unconventional reservoirs | |
| Deng et al. | The application of pure shear wave seismic data for gas reservoir delineation | |
| CN106353807A (en) | Fracture identification method and device | |
| Sirazhiev | Seismic attribute analysis of the mississippian chert at the wellington field, south-central kansas | |
| Xu et al. | Minimum semblance weighted stacking with polarity correction for surface microseismic data processing | |
| Parra et al. | Interwell seismic logging for formation continuity at the Gypsy test site, Oklahoma | |
| Parra et al. | Characterization of fractured low Q zones at the Buena Vista Hills reservoir, California | |
| Herwanger et al. | Applying time-lapse seismic methods to reservoir management and field development planning at South Arne, Danish North Sea | |
| Gupta | Using direct S-wave seismic modes for reservoir characterization in Wellington Field, Kansas | |
| Zhonghua et al. | Application of VSP survey for shale gas exploration | |
| Xiyuan et al. | Research and application of gas detection techniques using full-wave attributes in southwest China | |
| Sumantri et al. | Seismic Qualitative and Quantitative Analysis for Reservoir Characterization of Globigerina BH Gas Field, North East Java Basin. | |
| Lei et al. | Study on Application of Oil & Gas Reservoirs Prediction Using Multi-wave Seismic Technology | |
| Yang et al. | Inversion of Fracture Weakness Parameters Based on the 3CVSP Data? Part I: HTI Media Composed of A Single Fracture Set in the Isotropic Background Media |
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 |