CN109558613A - The inversion method and system of anisotropic rock physical model - Google Patents
The inversion method and system of anisotropic rock physical model Download PDFInfo
- Publication number
- CN109558613A CN109558613A CN201710892661.3A CN201710892661A CN109558613A CN 109558613 A CN109558613 A CN 109558613A CN 201710892661 A CN201710892661 A CN 201710892661A CN 109558613 A CN109558613 A CN 109558613A
- Authority
- CN
- China
- Prior art keywords
- parameter
- anisotropic
- clay
- formula
- inversion method
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention discloses the inversion method and system of a kind of anisotropic rock physical model, which includes: to establish the petrophysical model with transverse isotropy elastic characteristic；Based on bayesian theory, Posterior probability distribution is obtained；Based on petrophysical model and Posterior probability distribution, the objective function of inverting is established；The optimal solution of fracture intensity parameter is obtained based on objective function.The present invention has the advantages that proposing a kind of anisotropic rock Method of Physical Modeling, anisotropic clay and the microcrack that aligns can be comprehensively considered on the anisotropic influence of shale, it and is that constraint carries out rock physics inversion with prior information and log parameter based on Bayesian frame, it obtains the elastic parameter of clay mineral and aligns the density in crack, to further discuss its anisotropic character.
Description
Technical field
The present invention relates to oil gas technical field of physical geography, more particularly, to a kind of anisotropic rock physical model
Inversion method and system.
Background technique
Shale rock physics model the relationship for being intended to establish shale component, structure and fluid-filling Yu its elastic parameter,
And further analyze its resilient nature.Shale possesses more complicated mineral constituent and structure, this makes the research work to shale
Make relative difficulty.To the research of shale petrophysical model early stage can trace back to Vernik and Nur (1992), Vernik and
Liu (1997) etc. is based on anisotropy Backus average theory, has studied the North America for belonging to rich organic matter black shale type
The rock physics modeling problem of Bakken shale.The same period, Hornby et al. (1994) are approximate (SCA) using anisotropy self-compatibility
With the theoretical progress shale rock physics modeling of difference EFFECTIVE MEDIUM (DEM).
In recent years, many scholars utilize rock physics modeling method, discuss the composition and its fine structure of shale, and divide
Analyse its anisotropy attribute: Mba and Prasad (2010) analyzes the relationship of shale Mineralss and elastic anisotropy；
Spikes (2011) based on statistical method have studied fracture intensity parameter uncertainty and pore morphology to shale elasticity
Parameter influences；Jiang (2013) establishes petrophysical model to Wei Er pages of Haynes using self-compatibility model and Chapman theory
Rock has carried out parameter Estimation.Cause the anisotropic factor of shale complicated, existing rock physics modeling means can't consider
The influence of whole factors, therefore comprehensively consider that the anisotropic factor of shale is caused to be that current shale rock physics are built as far as possible
One important development direction of mould.
Existing technology just causes shale anisotropic some or two factors to be discussed mostly, each such as clay
Anisotropy, the distribution of organic matter, horizontal bedding and crack etc..Paradoxically with consideration influence factor as much as possible, constraint
Condition is insufficient, that is, lacks corresponding data.In addition, clay mineral is composition shale important component, in certain areas, its content is very
To can achieve 50% or more, therefore decide the quality of rock physics modeling to the estimation of its elastic parameter.Clay mineral is gathered around
There is complicated component, and there are the characteristics that aligning, this has differences its elastic parameter of the clay of different regions, existing
Research is mostly based on the measured clay elastic parameter of forefathers' experiment and studies, and has ignored its compositional difference and aligns institute
Caused by influence.
Therefore, it is necessary to develop a kind of anisotropic rock Method of Physical Modeling and system, can comprehensively consider each to different
Property clay and the microcrack that aligns do not used directly on the anisotropic influence of shale, and during considering clay mineral
The elastic parameter that forefathers measure, and direct inversion obtains the elastic parameter of clay mineral.
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 the inversion methods and system of a kind of anisotropic rock physical model, can be examined by synthesis
Anisotropic clay and the microcrack aligned are considered on the anisotropic influence of shale, and based on Bayesian frame with priori
Information and log parameter are that constraint carries out rock physics inversion, obtain the elastic parameter of clay mineral and align the close of crack
Degree, to further discuss its anisotropic character.
According to an aspect of the invention, it is proposed that a kind of inversion method of anisotropic rock physical model, the inverting
Method includes:
Establish the petrophysical model with transverse isotropy elastic characteristic；
Based on bayesian theory, Posterior probability distribution is obtained；
Based on the petrophysical model and the Posterior probability distribution, the objective function of inverting is established；
The optimal solution of fracture intensity parameter is obtained based on the objective function.
Preferably, establishing petrophysical model includes:
It is theoretical based on Hashi-Shtrikman boundary, calculate non-stick great soil group mineral and the kerogenic isotropism of organic matter
Elastic parameter；
Clay mineral is mixed with above-mentioned non-stick great soil group mineral, kerogen based on Backus average theory, calculating has VTI
The elasticity modulus of anisotropic solid matrix, and consider the heterogeneity of clay mineral elastic parameter；
The elastic parameter of fluid-mixing in hole-fracture spaces is sought based on Wood formula；
Multiple dimensioned hole-Fracture System is introduced into solid matrix based on Chapman theory, is obtained with VTI anisotropy
Shale reservoir elastic parameter；
It is converted based on Bond, in the case of calculating pendage, the TTI elastic anisotropy parameter of shale reservoir.
Preferably, the TTI elastic anisotropy parameter matrix of the shale reservoir are as follows:
In formula,For the elastic anisotropy parameter coefficient matrix of TTI model；
vpclayFor the velocity of longitudinal wave of clay mineral；
vsclayFor the shear wave velocity of clay mineral；
ε is the hole of parallel bedding or the density of microcrack.
Preferably, it is defined according to rock physics, in the situation known to well logging porosity φ, the hole of parallel bedding or micro-
The density ε in crack are as follows:
Then formula (1) are as follows:
In formula, α is pore components, is denary logarithm.
Preferably, the Posterior probability distribution are as follows:
P(B∣A)∝P(A∣B)P(B) (4)
In formula, P (B) is the prior probability that event B occurs；
P (A ∣ B) is likelihood function；
P (B ∣ A) is posterior probability.
Preferably, to three parameter vp of invertingclay, vsclayMeet ternary Gaussian Profile with α:
In formula, x be parameter composed by clay mineral P- and S-wave velocity and pore components denary logarithm to
Amount；
μ1Vector composed by mean value for parameter；
Σ1For the covariance matrix of parameter vector.
Preferably, likelihood function meets dimensional gaussian distribution:
In formula, y is vector composed by the vertical P- and S-wave velocity of TTI model；
μ2For the vector of well logging P- and S-wave velocity composition；
Σ2For covariance matrix, the degrees of tolerance for deviateing model solution practical logging speed decide；
Then formula (4) are as follows:
Preferably, the objective function obtained according to formula (7) are as follows:
In formula, JiFor the objective function of i-th of logging point；
yiVector is formed for the P- and S-wave velocity that the rock physics of i-th of logging point model；
For the vector of the actual measurement P- and S-wave velocity composition of i-th of logging point；
xiFor vector composed by i-th of logging point clay P- and S-wave velocity and pore components；
For vector composed by the priori mean value of i-th of logging point clay P- and S-wave velocity and pore components, yiFor
xiFunction.
Preferably, the optimal solution of the fracture intensity parameter is obtained when the objective function minimum.
According to another aspect of the invention, it is proposed that a kind of Inversion System of anisotropic rock physical model, the system
System includes:
Memory is stored with computer executable instructions；
Processor, the processor run the computer executable instructions in the memory, execute following steps:
Establish the petrophysical model with transverse isotropy elastic characteristic；
Based on bayesian theory, Posterior probability distribution is obtained；
Based on the petrophysical model and the Posterior probability distribution, the objective function of inverting is established；
Obtain the optimal solution of the objective function.
The inversion method and system of a kind of anisotropic rock physical model according to the present invention, the advantage is that: proposing
A kind of anisotropic rock Method of Physical Modeling can comprehensively consider anisotropic clay and the microcrack that aligns to shale
Anisotropic influence, and be that constraint carries out rock physics inversion with prior information and log parameter based on Bayesian frame,
It obtains the elastic parameter of clay mineral and aligns the density in crack, to further discuss its anisotropic character.
Method and system of the invention have 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 and subsequent specific implementation being incorporated herein
It is stated in detail in example, these the drawings and specific embodiments are used together to explain specific principle 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.
Fig. 1 shows a kind of inverting of anisotropic rock physical model of an exemplary embodiment according to the present invention
The flow chart of the step of method.
Fig. 2 shows the schematic diagrames for establishing petrophysical model of an exemplary embodiment according to the present invention.
Fig. 3 a and Fig. 3 b respectively illustrate the TTI model longitudinal wave of an exemplary embodiment according to the present invention, shear wave speed
The contrast schematic diagram of degree and actual speed log.
The crack that the petrophysical model inverting that Fig. 4 shows an exemplary embodiment according to the present invention obtains is in length and breadth
Than the schematic diagram of curve.
The petrophysical model inverting that Fig. 5 a and Fig. 5 b respectively illustrate an exemplary embodiment according to the present invention obtains
The velocity of longitudinal wave of the clay mineral obtained and the schematic diagram of shear wave velocity.
Fig. 6 shows the schematic diagram of the fracture spacing curve of an exemplary embodiment according to the present invention.
The longitudinal wave anisotropic parameters ε and shear wave that Fig. 7 shows an exemplary embodiment according to the present invention are respectively to different
The schematic diagram of property parameter gamma curve.
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.
The present invention provides a kind of inversion method of anisotropic rock physical model, which includes:
Establish the petrophysical model with transverse isotropy elastic characteristic；
Based on bayesian theory, Posterior probability distribution is obtained；
Based on petrophysical model and Posterior probability distribution, the objective function of inverting is established；
The optimal solution of fracture intensity parameter is obtained based on objective function.
Petrophysical model is established by the shale reservoir with transverse isotropy parameter attribute.
Preferably, establishing petrophysical model includes:
It is theoretical based on Hashi-Shtrikman boundary, calculate non-stick great soil group mineral (quartz, dolomite, limestone, Huang Tie
Mine) and the kerogenic isotropic elasticity parameter of organic matter, the research based on forefathers, it is assumed that kerogen is spatially in random
Lumps distribution；
Clay mineral is mixed with above-mentioned non-stick great soil group mineral, kerogen based on Backus average theory, calculating has VTI
The elasticity modulus of (Transversely isotropy with a vertical axis) anisotropic solid matrix, and
Consider the heterogeneity of clay mineral elastic parameter；
The elastic parameter of fluid-mixing in hole-fracture spaces (gas and water) is sought based on Wood formula；
Multiple dimensioned hole-Fracture System is introduced into solid matrix based on Chapman theory, is obtained with VTI anisotropy
Shale reservoir elastic parameter；
It is converted based on Bond, in the case of calculating pendage, the TTI elastic anisotropy parameter of shale reservoir.
It proposes a kind of anisotropic rock Method of Physical Modeling, anisotropic clay can be comprehensively considered and aligns
Microcrack is that constraint carries out rock with prior information and log parameter on the anisotropic influence of shale, and based on Bayesian frame
Stone physics inverting obtains the elastic parameter of clay mineral and aligns the density in crack, to further discuss it respectively to different
Property feature.
In above-mentioned petrophysical model, clay mineral velocity of longitudinal wave, shear wave velocity and hole or fracture spacing are undetermined
Parameter, other parameters can be determined or be estimated by log data, therefore design the inversion method of petrophysical model, page for convenience
The TTI elastic anisotropy parameter matrix of rock reservoir are as follows:
In formula,For the elastic anisotropy parameter coefficient matrix of TTI model；
vpclayFor the velocity of longitudinal wave of clay mineral；
vsclayFor the shear wave velocity of clay mineral；
ε is the hole of parallel bedding or the density of microcrack.
Preferably, it is defined according to rock physics, in the situation known to well logging porosity φ, the hole of parallel bedding
The density ε of gap or microcrack are as follows:
Then formula (1) are as follows:
In formula, α is pore components, is denary logarithm.
It by log data calculating clay mineral velocity of longitudinal wave, shear wave velocity and is these three parameters of pore components, if
Count the inversion algorithm of petrophysical model.
Using clay velocity of longitudinal wave, clay shear wave velocity and the pore components in modeling process as fitting parameter, with
Velocity of longitudinal wave and shear wave velocity that well logging obtains carry out the inverting of petrophysical model as constraint.
Preferably, Posterior probability distribution are as follows:
P(B∣A)∝P(A∣B)P(B) (4)
In formula, P (B) is the prior probability that event B occurs；
P (A ∣ B) is likelihood function；
P (B ∣ A) is posterior probability.
Assume that the velocity of longitudinal wave of clay mineral and shear wave velocity meet Gaussian Profile, velocity of longitudinal wave and shear wave speed in the present invention
The mean value of degree takes the measurement result of forefathers, and pore components α is denary logarithm, also meets Gaussian Profile, hole is in length and breadth
Than the changing greatly with lithology in different regions, biggish variance is chosen.Meet a ternary to three parameters of inverting in this way
Gaussian Profile, the formula met are as follows:
In formula, x be parameter composed by clay mineral P- and S-wave velocity and pore components denary logarithm to
Amount；
μ1Vector composed by mean value for parameter；
Σ1For the covariance matrix of parameter vector.
Wherein, covariance matrix not only reflects that each variable deviates the degree of mean value, but also can reflect between different variables
Correlation.Present invention assumes that clay mineral velocity of longitudinal wave and shear wave velocity are positive correlation, i.e. velocity of longitudinal wave becomes larger, shear wave
Speed also becomes larger, and the clay speed that covariance is measured by forefathers counts to obtain, while assuming pore components and clay
Velocity of longitudinal wave, shear wave velocity do not have correlation, i.e. covariance is 0.
In order to establish the objective function of inverting, it is assumed that forward modeling result and Log-constrained parameter (the longitudinal wave speed measured of model
Degree and shear wave velocity) between residual error, i.e. likelihood function also meets dimensional gaussian distribution.
Preferably, likelihood function meets dimensional gaussian distribution:
In formula, y is vector composed by the vertical P- and S-wave velocity of TTI model；
μ2For the vector of well logging P- and S-wave velocity composition；
Σ2For covariance matrix, the degrees of tolerance for deviateing model solution practical logging speed decide；
Then formula (4) are as follows:
Wherein, the maximum point of posterior probability density corresponds to optimal solution, therefore, available objective function are as follows:
In formula, JiFor the objective function of i-th of logging point；
yiVector is formed for the P- and S-wave velocity that the rock physics of i-th of logging point model；
For the vector of the actual measurement P- and S-wave velocity composition of i-th of logging point；
xiFor vector composed by i-th of logging point clay P- and S-wave velocity and pore components；
For vector composed by the priori mean value of i-th of logging point clay P- and S-wave velocity and pore components, yiFor
xiFunction.
Preferably, the optimal solution of fracture intensity parameter is obtained when objective function minimum.
It can establish rock constituents, porosity, pore morphology and fluidity using the inversion method of petrophysical model
Contacting between matter and rock integral, flexible parameter, further analyzes its anisotropic character, has to shale oil-gas reservoir is found
Important reference significance.Meanwhile the crack in shale is important hydrocarbon migration channel and reservoir space, fracture spacing bigger generation
Table higher reservoir space, the present invention by the fracture intensity parameter that rock physics inversion obtains be evaluation reservoir storage capacity
One of important parameter.
The present invention also provides a kind of Inversion Systems of anisotropic rock physical model, which is characterized in that system includes:
Memory is stored with computer executable instructions；
Processor, the computer executable instructions in processor run memory execute following steps:
Establish the petrophysical model with transverse isotropy elastic characteristic；
Based on bayesian theory, Posterior probability distribution is obtained；
Based on petrophysical model and Posterior probability distribution, the objective function of inverting is established；
Obtain the optimal solution of objective function.
Embodiment
Fig. 1 shows a kind of inverting of anisotropic rock physical model of an exemplary embodiment according to the present invention
The flow chart of the step of method.
As shown in Figure 1, a kind of inversion method of anisotropy petrophysical model of the present embodiment, the inversion method include:
Establish the petrophysical model with transverse isotropy elastic characteristic；
Based on bayesian theory, Posterior probability distribution is obtained；
Based on petrophysical model and Posterior probability distribution, the objective function of inverting is established；
The optimal solution of fracture intensity parameter is obtained based on objective function.
Fig. 2 shows the schematic diagrames for establishing petrophysical model of an exemplary embodiment according to the present invention.
As shown in Fig. 2, establishing petrophysical model and including:
It is theoretical based on Hashi-Shtrikman boundary, calculate non-stick great soil group mineral (quartz, dolomite, limestone, Huang Tie
Mine) and the kerogenic isotropic elasticity parameter of organic matter, the research based on forefathers, it is assumed that kerogen is spatially in random
Lumps distribution；
Clay mineral is mixed with above-mentioned non-stick great soil group mineral, kerogen based on Backus average theory, calculating has VTI
The elasticity modulus of (Transversely isotropy with a vertical axis) anisotropic solid matrix, and
Consider the heterogeneity of clay mineral elastic parameter；
The elastic parameter of fluid-mixing in hole-fracture spaces (gas and water) is sought based on Wood formula；
Multiple dimensioned hole-Fracture System is introduced into solid matrix based on Chapman theory, is obtained with VTI anisotropy
Shale reservoir elastic parameter；
It is converted based on Bond, in the case of calculating pendage, the TTI elastic anisotropy parameter of shale reservoir.
At home certain well log well one by one shop using the above method carry out petrophysical model inverting.
Fig. 3 a and Fig. 3 b respectively illustrate the TTI model longitudinal wave of an exemplary embodiment according to the present invention, shear wave speed
The contrast schematic diagram of degree and actual speed log.
As shown in Figure 3a and Figure 3b shows, dark line segment form show the vertical velocity of longitudinal wave of TTI petrophysical model (described in Fig. 3 a) and
Shear wave velocity (shown in Fig. 3 b), light line segment form show the velocity of longitudinal wave and shear wave velocity of well logging actual measurement.It can be seen that utilizing this hair
Bright objective function carries out constraint and population inversion method, and log has obtained preferable fitting, and regression criterion is smaller, says
Bright objective function of the invention has stronger restraining force.
Wherein, the purpose of population inverting is the optimal solution of model parameter when obtaining objective function minimum.
The hole that the petrophysical model inverting that Fig. 4 shows an exemplary embodiment according to the present invention obtains is in length and breadth
Than the schematic diagram of curve.
As shown in figure 4, the pore components of shale interval concentrate between 0.1~0.3, before overall pore components meet
Statistical result of the people to this area's shale pore components.
The petrophysical model inverting that Fig. 5 a and Fig. 5 b respectively illustrate an exemplary embodiment according to the present invention obtains
The velocity of longitudinal wave of the clay mineral obtained and the schematic diagram of shear wave velocity.
As shown in figure 5 a and 5b, the clay velocity of longitudinal wave (shown in Fig. 5 a) and shear wave velocity that inverting obtains (shown in Fig. 5 b)
Greater than forefathers' as a result, the reason of causing the result can be clay that the clay component in the shale is studied compared to forefathers its
In illite content it is higher, and illite P- and S-wave velocity with higher.This inversion result further illustrates this hair
Bright proposed refutation strategy is able to reflect the influence of change of component bring speed in clay mineral.
Fig. 6 shows the schematic diagram of the fracture spacing curve of an exemplary embodiment according to the present invention.Fig. 7 is shown
The signal of the longitudinal wave anisotropic parameters ε and shear wave anisotropic parameters gamma curve of an exemplary embodiment according to the present invention
Figure.
As shown in fig. 7, be parallel to stratification horizontal fracture enhancing stratum VTI anisotropy, it illustrates longitudinal wave respectively to
Anisotropic parameter ε and shear wave anisotropic parameters gamma curve can be used for analyzing the anisotropic character of this area's shale, parameter ε generation
Horizontally and vertically velocity of longitudinal wave difference, parameter γ represent level side caused by anisotropy caused by table anisotropy
To with vertical direction velocity of longitudinal wave difference.
The embodiment of the present invention is described above, above description is exemplary, and non-exclusive, and also not
It is limited to disclosed embodiment.Without departing from the scope and spirit of embodiment described, for the art
Those of ordinary skill for many modifications and changes are obvious.The selection of term used herein, it is intended to best
Principle, practical application or the improvement to the technology in market of embodiment are explained in ground, or make the other common of the art
Technical staff can understand embodiments disclosed herein.
Claims (10)
1. a kind of inversion method of anisotropic rock physical model, which is characterized in that the inversion method includes:
Establish the petrophysical model with transverse isotropy elastic characteristic；
Based on bayesian theory, Posterior probability distribution is obtained；
Based on the petrophysical model and the Posterior probability distribution, the objective function of inverting is established；
The optimal solution of fracture intensity parameter is obtained based on the objective function.
2. the inversion method of anisotropic rock physical model according to claim 1, wherein establish petrophysical model
Include:
It is theoretical based on Hashi-Shtrikman boundary, calculate non-stick great soil group mineral and the kerogenic isotropic elasticity of organic matter
Parameter；
Clay mineral is mixed with above-mentioned non-stick great soil group mineral, kerogen based on Backus average theory, calculate have VTI respectively to
The elasticity modulus of anisotropic solid matrix, and consider the heterogeneity of clay mineral elastic parameter；
The elastic parameter of fluid-mixing in hole-fracture spaces is sought based on Wood formula；
Multiple dimensioned hole-Fracture System is introduced into solid matrix based on Chapman theory, is obtained with the anisotropic page of VTI
Rock reservoir elastic parameter；
It is converted based on Bond, in the case of calculating pendage, the TTI elastic anisotropy parameter of shale reservoir.
3. the inversion method of anisotropic rock physical model according to claim 2, wherein the shale reservoir
TTI elastic anisotropy parameter matrix are as follows:
In formula,For the elastic anisotropy parameter coefficient matrix of TTI model；
vpclayFor the velocity of longitudinal wave of clay mineral；
vsclayFor the shear wave velocity of clay mineral；
ε is the hole of parallel bedding or the density of microcrack.
4. the inversion method of anisotropic rock physical model according to claim 3, wherein fixed according to rock physics
Justice, in the situation known to well logging porosity φ, the hole of parallel bedding or the density ε of microcrack are as follows:
Then formula (1) are as follows:
In formula, α is pore components, is denary logarithm.
5. the inversion method of anisotropic rock physical model according to claim 4, wherein the Posterior probability distribution
Are as follows:
P(B∣A)∝P(A∣B)P(B) (4)
In formula, P (B) is the prior probability that event B occurs；
P (A ∣ B) is likelihood function；
P (B ∣ A) is posterior probability.
6. the inversion method of anisotropic rock physical model according to claim 5, wherein three parameters to inverting
vpclay, vsclayMeet ternary Gaussian Profile with α:
In formula, x is parameter vector composed by clay mineral P- and S-wave velocity and pore components denary logarithm；
μ1Vector composed by mean value for parameter；
Σ1For the covariance matrix of parameter vector.
7. the inversion method of anisotropic rock physical model according to claim 6, wherein likelihood function meets two dimension
Gaussian Profile:
In formula, y is vector composed by the vertical P- and S-wave velocity of TTI model；
μ2For the vector of well logging P- and S-wave velocity composition；
Σ2For covariance matrix, the degrees of tolerance for deviateing model solution practical logging speed decide；
Then formula (4) are as follows:
8. the inversion method of anisotropic rock physical model according to claim 7, wherein obtained according to formula (7)
The objective function are as follows:
In formula, JiFor the objective function of i-th of logging point；
yiVector is formed for the P- and S-wave velocity that the rock physics of i-th of logging point model；
For the vector of the actual measurement P- and S-wave velocity composition of i-th of logging point；
xiFor vector composed by i-th of logging point clay P- and S-wave velocity and pore components；
For vector composed by the priori mean value of i-th of logging point clay P- and S-wave velocity and pore components, yiFor xi's
Function.
9. the inversion method of anisotropic rock physical model according to claim 8, wherein the objective function is minimum
When obtain the optimal solution of the fracture intensity parameter.
10. a kind of Inversion System of anisotropic rock physical model, which is characterized in that the system comprises:
Memory is stored with computer executable instructions；
Processor, the processor run the computer executable instructions in the memory, execute following steps:
Establish the petrophysical model with transverse isotropy elastic characteristic；
Based on bayesian theory, Posterior probability distribution is obtained；
Based on the petrophysical model and the Posterior probability distribution, the objective function of inverting is established；
Obtain the optimal solution of the objective function.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710892661.3A CN109558613B (en) | 2017-09-27 | 2017-09-27 | Inversion method and system of anisotropic rock physical model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710892661.3A CN109558613B (en) | 2017-09-27 | 2017-09-27 | Inversion method and system of anisotropic rock physical model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109558613A true CN109558613A (en) | 2019-04-02 |
| CN109558613B CN109558613B (en) | 2021-11-23 |
Family
ID=65864276
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710892661.3A Active CN109558613B (en) | 2017-09-27 | 2017-09-27 | Inversion method and system of anisotropic rock physical model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109558613B (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110135112A (en) * | 2019-06-04 | 2019-08-16 | 南京信息工程大学 | A kind of shale anisotropy estimation method based on particle filter algorithm |
| CN110276094A (en) * | 2019-04-30 | 2019-09-24 | 清华大学 | Current elements 3-d inversion method based on Bayes's elastic network(s) regularization method |
| CN110512622A (en) * | 2019-08-23 | 2019-11-29 | 郑州安源工程技术有限公司 | Face the vertical fracture filling grouting method of slope surface rock mass |
| CN111175824A (en) * | 2020-01-06 | 2020-05-19 | 中国石油化工股份有限公司 | Time-frequency joint domain seismic inversion method under lithofacies driving |
| CN111751886A (en) * | 2020-08-12 | 2020-10-09 | 重庆地质矿产研究院 | Shale gas reservoir fracture modeling method based on microseism monitoring data |
| CN113687412A (en) * | 2020-05-18 | 2021-11-23 | 中国石油化工股份有限公司 | Method and device for predicting pressure of stratum between salts, electronic equipment and medium |
| CN114280670A (en) * | 2021-12-13 | 2022-04-05 | 华能煤炭技术研究有限公司 | Multi-physical-model well logging shear wave velocity curve reconstruction method and system and electronic equipment |
| CN117420202A (en) * | 2023-10-18 | 2024-01-19 | 中国地质大学（北京） | Shale type determining method, shale type determining device, shale type determining terminal and shale type storing medium |
| CN117420202B (en) * | 2023-10-18 | 2024-06-21 | 中国地质大学（北京） | Shale type determining method, shale type determining device, shale type determining terminal and shale type storing medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050090986A1 (en) * | 2003-10-27 | 2005-04-28 | Fugro-Jason | Method and device for the generation and application of anisotropic elastic parameters |
| US20110292766A1 (en) * | 2010-05-26 | 2011-12-01 | Ran Bachrach | Estimating Anisotropic Parameters |
| CN105203739A (en) * | 2015-10-29 | 2015-12-30 | 中国石油大学(北京) | Rock physical model-based shale stratum anisotropy parameter prediction method |
| CN107203005A (en) * | 2016-03-18 | 2017-09-26 | 中国石油化工股份有限公司 | A kind of method that quantification calculates crack characterising parameter |
-
2017
- 2017-09-27 CN CN201710892661.3A patent/CN109558613B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050090986A1 (en) * | 2003-10-27 | 2005-04-28 | Fugro-Jason | Method and device for the generation and application of anisotropic elastic parameters |
| US20110292766A1 (en) * | 2010-05-26 | 2011-12-01 | Ran Bachrach | Estimating Anisotropic Parameters |
| CN105203739A (en) * | 2015-10-29 | 2015-12-30 | 中国石油大学(北京) | Rock physical model-based shale stratum anisotropy parameter prediction method |
| CN107203005A (en) * | 2016-03-18 | 2017-09-26 | 中国石油化工股份有限公司 | A kind of method that quantification calculates crack characterising parameter |
Non-Patent Citations (1)
| Title |
|---|
| 刘喜武，董宁，刘宇巍: ""裂缝性孔隙介质频变AVAZ反演方法研究进展"", 《石油物探》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110276094A (en) * | 2019-04-30 | 2019-09-24 | 清华大学 | Current elements 3-d inversion method based on Bayes's elastic network(s) regularization method |
| CN110135112A (en) * | 2019-06-04 | 2019-08-16 | 南京信息工程大学 | A kind of shale anisotropy estimation method based on particle filter algorithm |
| CN110512622A (en) * | 2019-08-23 | 2019-11-29 | 郑州安源工程技术有限公司 | Face the vertical fracture filling grouting method of slope surface rock mass |
| CN111175824A (en) * | 2020-01-06 | 2020-05-19 | 中国石油化工股份有限公司 | Time-frequency joint domain seismic inversion method under lithofacies driving |
| CN113687412A (en) * | 2020-05-18 | 2021-11-23 | 中国石油化工股份有限公司 | Method and device for predicting pressure of stratum between salts, electronic equipment and medium |
| CN113687412B (en) * | 2020-05-18 | 2024-03-26 | 中国石油化工股份有限公司 | Method and device for predicting formation pressure between salts, electronic equipment and medium |
| CN111751886A (en) * | 2020-08-12 | 2020-10-09 | 重庆地质矿产研究院 | Shale gas reservoir fracture modeling method based on microseism monitoring data |
| CN114280670A (en) * | 2021-12-13 | 2022-04-05 | 华能煤炭技术研究有限公司 | Multi-physical-model well logging shear wave velocity curve reconstruction method and system and electronic equipment |
| CN117420202A (en) * | 2023-10-18 | 2024-01-19 | 中国地质大学（北京） | Shale type determining method, shale type determining device, shale type determining terminal and shale type storing medium |
| CN117420202B (en) * | 2023-10-18 | 2024-06-21 | 中国地质大学（北京） | Shale type determining method, shale type determining device, shale type determining terminal and shale type storing medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109558613B (en) | 2021-11-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109558613A (en) | The inversion method and system of anisotropic rock physical model | |
| CN107203005B (en) | Method for quantitatively calculating crack description parameters | |
| EP2583121B1 (en) | Production simulator for simulating a mature hydrocarbon field | |
| CN109655903B (en) | Shale layer transverse wave velocity prediction method and system | |
| Hummel et al. | Nonlinear diffusion-based interpretation of induced microseismicity: A Barnett Shale hydraulic fracturing case study | |
| CN103258091A (en) | Unconventional oil and gas reservoir horizontal well section three-dimensional rock mass mechanics modeling method and device | |
| CN104297785A (en) | Lithofacies constrained reservoir physical property parameter inversion method and device | |
| CN105974474A (en) | Method for evaluating reservoir pore type based on pore replacement | |
| CN103163553A (en) | Earthquake hydrocarbon detection method and detection device based on multiple pore medium model | |
| CN104570065B (en) | Method for quantitatively inverting porosity by using seismic wave impedance | |
| CN106199712B (en) | A kind of method and device in definite pressure break casing deformation region | |
| Gupta et al. | Residual Oil Zone Recovery Evaluation and Forecast Methodology: A Wasson Field Case Study | |
| CN110320575A (en) | Method and device is determined based on the shale content of organic matter of petrophysical model | |
| Bachrach et al. | A nonlinear and anisotropic data-driven litho-petroelastic inversion for single-loop subsurface characterization: Theory, uncertainties, and a case study | |
| CN109459791A (en) | A kind of method and system determining river location using log | |
| Shchipanov et al. | A new approach to deformable fractured reservoir characterization: case study of the Ekofisk field | |
| CN115964842A (en) | Method and system for determining natural gas hydrate resource amount | |
| CN114460639B (en) | Shale oil reservoir permeability prediction method and device | |
| Pettersson et al. | Dynamic estimates of extreme-case CO2 storage capacity for basin-scale heterogeneous systems under geological uncertainty | |
| CN105137480A (en) | Method and apparatus for determining pre-stack time migration velocity of VTI medium | |
| CN112379437B (en) | Shale reservoir anisotropy parameter solving method and device | |
| CN114185110B (en) | Method, device and medium for determining relative content of spatial structure in tight reservoir | |
| Wang et al. | Study on a New Rock Fracability Evaluation Model of Shale Gas Reservoir | |
| CN114114390B (en) | Oil and gas exploration method, device, equipment and storage medium | |
| Jervis et al. | Deep learning applied to seismic facies classification |
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 |