CN106556866B - A method of using rock physics template to reservoir quantitative interpretation - Google Patents
A method of using rock physics template to reservoir quantitative interpretation Download PDFInfo
- Publication number
- CN106556866B CN106556866B CN201510618211.6A CN201510618211A CN106556866B CN 106556866 B CN106556866 B CN 106556866B CN 201510618211 A CN201510618211 A CN 201510618211A CN 106556866 B CN106556866 B CN 106556866B
- Authority
- CN
- China
- Prior art keywords
- rock
- rock physics
- modulus
- physics template
- density
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- 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/624—Reservoir parameters
- G01V2210/6244—Porosity
Abstract
The invention proposes a kind of using rock physics template to the method for reservoir quantitative interpretation, comprising: rock sample test analysis；Calculate rock physics template initial model；Optimize rock physics template initial model；Pre-stack seismic inversion result is explained using the rock physics template initial model of optimization；Obtain oily window.Method of the invention can be used for the explanation of interpretation of logging data and pre-stack seismic inversion result.Rock physics template method for building up provided by the invention can quickly be set up for logging well and the rock physics template of seismic interpretation, by analyzing oil-gas possibility to distribution of the data in rock physics template, can reduce the multi-solution of seismic data interpretation.
Description
Technical field
The present invention relates to petrophysics and geophysical well logging technology fields, and in particular to a kind of to utilize rock physics mould
Method of the plate to reservoir quantitative interpretation.
Background technique
Petrophysics is the bridge for connecting petrophysical parameter, fluid properties parameter and seismic elastic parameter, is earthquake
The quantitative interpretation of data provides solid foundation, greatly reduces the multi-solution of seismic interpretation.In recent years, both at home and abroad in earthquake rock
Stone physical basis research aspect achieve significant progress, while also more successfully by research achievement be applied to Lithology Discrimination,
Fluid detection and porosity inversion etc. played an important role in fields such as seismic prospecting, oil reservoir developments.
With deepening continuously for oil-gas exploration and development, oil gas reservoir body becomes increasingly complex, and explores emphasis by original construction
Oil-gas reservoir is shifted to lithologic deposit, subtle pool.One side exploration targets gradually deepens, becomes smaller, and difficulties in exploration increases, right
The required precision of reservoir prediction is higher.On the other hand, existence of the final purpose of reservoir prediction in addition to judging reservoir, prediction storage
Layer development degree it is fine or not, while needing to determine the property of contained fluid in reservoir.
Largely scientific investigations showed that, China possesses Stratigraphic and subtle reservoirs resource abundant, there is huge exploration and development
Potentiality.And seismic inversion is the important technical for obtaining lithologic character stratum elastic parameter, passes through the available longitudinal wave of post-stack inversion
Impedance data body, by the available p-wave impedance of prestack inversion, S-wave impedance, density or velocity of longitudinal wave, shear wave velocity, close
The parameters such as degree.These elastic parameters are widely used in oil-bearing reservoir prediction, achieve many good prediction effects, but
There is the example of failure, the purpose distinguished lithology and identify fluid is not achieved.Believe excavating the reservoir that seismic inversion data are contained
Breath aspect has scholar and scientific research personnel to do a large amount of research and production application work, such as the fluids such as Poisson impedance, speed ratio
The appearance of recognition factor.It to the explanation of these elastic parameters and its transformation is carried out in qualitative degree at present, it is single to answer
The multi-solution for carrying out reservoir prediction work with a kind of parameter is stronger, brings risk to the explanation of seismic data, limits earthquake
The application and explanation work of inversion result.
Summary of the invention
Explanation means for seismic inversion and the shortcomings that be of limited application, the method for the present invention is based on rock physics
Model establishes the rock physics template for seismic data quantitative interpretation.Multi-parameter is intersected with rock physics template
Analysis, analyzes the distribution of each parameter in rock physics template, and by the elastic parameter of destination region section, slice or
Three-dimensional space is sketched the contours, and realizes explanation and prediction work to target geological body, reduces the multi-solution for explaining work.
The present invention provides a kind of method using rock physics template to reservoir quantitative interpretation, for instruct seismic inversion and
Oil and gas reservoir forecast of distribution and fluid detection are realized in the explanation of inversion result.Method of the invention, comprising: rock sample test analysis；
Calculate rock physics template initial model；Optimize rock physics template initial model；It is initial using the rock physics template of optimization
Model explains pre-stack seismic inversion result；Obtain oily window.
Further, in rock sample test analysis step, by carrying out test analysis to reservoir rock sample, composition rock is obtained
Mineralogical composition and parameter of pore structure.
Optionally, in sand shale, including obtaining sandstone and the elasticity modulus of mud stone, density, hole ellipticity, the bullet
Property modulus includes bulk modulus, modulus of shearing；Wherein the elasticity modulus of sandstone, density are that symbol is expressed as Ks,μs,ρs, mud stone
Elasticity modulus, density Kc,μc,ρc；Hole ellipticity is αs,αc；Parameter (αs,αc,Ks,μs,ρs,Kc,μc,ρc) be rock physics just
Drill the input data of model.
Optionally, calculate rock physics template initial model, including for tight sand using dry rock approximation and
Gassmann formula establishes the petrophysical model of tight sand, and calculating speed, density are with rock porosity, water saturation
And the conversion curve of shale content, realize rock physics forward modelling.
Further, the dry rock is approximately to consider pore structure on the basis of given Rock Matrix elasticity modulus
Parameter is come particular by porosity and parameter of pore structure to initial substrate elasticity come the method for calculating dry elastic modulus of rock
Modulus is modified, and is obtained rock and is in dry elastic modulus of rock of porosity when being φ, when pore structure is α:
K (φ)=K0(1-φ)p, μ (φ)=μ0(1-φ)q (1)
K0And μ0It is the elasticity modulus of matrix, p and q are the function of volumn concentration and pore structure, K (φ) and μ (φ)
Be porosity be φ when dry elastic modulus of rock, vlRepresent the volumn concentration of sandstone and mud stone, respectively vsAnd vc,
Tiijj(αl) and F (αl) it is parameter of pore structure.
Optionally, the Gassmann formula can be calculated the springform of fluid saturated rocks by dry elastic modulus of rock
Thus amount calculates the velocity of longitudinal wave and shear wave velocity of isotropic medium:
ρsat=φ (Sw·ρw+(1-Sw·ρg))+(1-φ)·ρma (5)
Wherein, KfThe bulk modulus of fluid, K are saturated for holesatAnd μsatFor fluid saturated rocks bulk modulus and shearing
Modulus, SwFor water saturation, ρsatIt is the density of fluid saturated rocks, ρw、ρg、ρmaIt is the density of water, the density of gas and rock
The density of matrix, KdryAnd μdryIt is the bulk modulus and modulus of shearing of dry rock, Vρ estAnd Vs estIt is isotropic medium respectively
Velocity of longitudinal wave and shear wave velocity.
Further, shale content is set as definite value, gives water saturation, calculating speed and density with the change of porosity
Change curve；Above-mentioned curve is converted, p-wave impedance and speed ratio are obtained, then shows transformed curve in abscissa
For p-wave impedance, ordinate is to obtain rock physics template initial model, wherein p-wave impedance is vertical in the cross plot of speed ratio
The product of wave velocity and density, speed ratio are the ratio of velocity of longitudinal wave and shear wave velocity.
Further, rock physics template initial model is optimized using actual measurement shear wave logging data, is optimized
Rock physics template afterwards:
Wherein,WithFor the longitudinal wave and shear wave velocity of actual measurement, VpAnd VsFor the longitudinal wave and shear wave velocity of forward modeling, λ is cross
Wave weight coefficient, value range are between 0 to 1.Wherein by adjusting sandstone and mud stone in rock constituents elasticity modulus and
Hole ellipticity, so that the p-and s-wave velocity residual error of the p-and s-wave velocity of forward modeling and actual measurement is minimum, sandstone and mud adjusted at this time
The elasticity modulus and hole ellipticity (α of rocks,αc,Ks,μs,ρs,Kc,μc,ρc) be exactly optimization after rock physics template parameter.
Further, p-and s-wave velocity that prestack inversion obtains, density are converted, obtains p-wave impedance and longitudinal and shear wave
Speed ratio analyzes the distribution of the p-wave impedance and speed ratio of prestack inversion in rock physics template, by where inversion result
The porosity in region, water saturation analyze oil-gas possibility, obtain oily window.Wherein p-wave impedance be velocity of longitudinal wave with
The product of density, speed ratio are the ratio of velocity of longitudinal wave and shear wave velocity.
Method of the invention can be used for the explanation of interpretation of logging data and pre-stack seismic inversion result.It is provided by the invention
Rock physics template method for building up can quickly be set up for logging well and the rock physics template of seismic interpretation, and logarithm is passed through
Oil-gas possibility is analyzed according to the distribution in rock physics template, can reduce the multi-solution of seismic data interpretation.
Detailed description of the invention
Disclosure illustrative embodiments are described in more detail in conjunction with the accompanying drawings, the disclosure above-mentioned and its
Its purpose, feature and advantage will be apparent, wherein in disclosure illustrative embodiments, identical reference label
Typically represent same parts.
Fig. 1 is flow chart of the invention.
Fig. 2 is speed ratio-p-wave impedance rock physics template in embodiment 1.
Fig. 3 is the oily window in speed ratio-p-wave impedance rock physics template in embodiment 1, by the number of circle choosing
It is shown according in slice of data is shown with particular color, to sketch the contours of the distribution and form in especially potentiality region.
Specific embodiment
The preferred embodiment of the disclosure is more fully described below with reference to accompanying drawings.Although showing the disclosure in attached drawing
Preferred embodiment, however, it is to be appreciated that may be realized in various forms the disclosure without the embodiment party that should be illustrated here
Formula is limited.On the contrary, these embodiments are provided so that this disclosure will be more thorough and complete, and can be by the disclosure
Range is completely communicated to those skilled in the art.
The present invention provides a kind of using rock physics template to the method for reservoir quantitative interpretation, can be used for well-log information solution
It releases and the explanation of pre-stack seismic inversion result.Rock physics template method for building up provided by the invention can quickly be set up
For logging well and the rock physics template of seismic interpretation, by analyzing oily to distribution of the data in rock physics template
Property, it can reduce the multi-solution of seismic data interpretation.
In order to establish the rock physics template for tight sandstone reservoir quantitative interpretation, the present invention first to reservoir rock sample into
Row test analysis obtains the mineralogical composition and parameter of pore structure of composition rock, then calculates rock bullet using petrophysical model
Property parameter with the variation of porosity, water saturation and shale content, obtain rock physics template initial model, recycle actual measurement
Shear wave logging data optimizes rock physics template initial model, and the rock physics template after being optimized will finally fold
The elastic parameter that preceding inverting obtains is simultaneously displayed in cross plot with rock physics template, is existed by the elastic parameter of prestack inversion
Distribution in rock physics template to carry out quantitative interpretation to reservoir and oil-gas possibility.
As shown in Figure 1, the present invention discloses a kind of method using rock physics template to reservoir quantitative interpretation, belong to rock
Physics and geophysical well logging technology field.The described method includes: rock sample test analysis；Calculate rock physics template introductory die
Type；Optimize rock physics template initial model；Using the rock physics template initial model of optimization to pre-stack seismic inversion result
It explains；Obtain oily window.
Further, the method comprise the steps that reservoir rock sample carry out test analysis obtain composition rock mineral at
Point and parameter of pore structure；Rock Elastic Parameters are calculated using petrophysical model to contain with porosity, water saturation and shale
The variation of amount obtains rock physics template initial model；Using actual measurement shear wave logging data to rock physics template initial model
It optimizes, the rock physics template after being optimized；The elastic parameter that prestack inversion is obtained with rock physics template simultaneously
It is shown in cross plot, by distribution of the elastic parameter of prestack inversion in rock physics template come to reservoir and oil-gas possibility
Carry out quantitative interpretation.
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.
As an embodiment, the method for the present invention includes the following steps:
Test analysis is carried out to rock sample first.Test analysis is carried out to tight sandstone reservoir rock sample and obtains the mine of composition rock
Object ingredient and parameter of pore structure.In sand shale, the hole of elasticity modulus and density, sandstone and mud stone including sandstone and mud stone
The parameters such as gap ellipticity.The elasticity modulus includes bulk modulus, modulus of shearing；The elasticity modulus of the sandstone and mud stone,
Density and hole ellipticity are the input datas of subsequent rock physics forward model.The elasticity modulus of sandstone, density are symbol expression
For Ks,μs,ρs,；Elasticity modulus, the density K of mud stonec,μc,ρc；Hole ellipticity is αs,αc.This group of parameter (αs,αc,Ks,μs,
ρs,Kc,μc,ρc) be rock physics forward model input data.
It is built for tight sand using the petrophysical model that dry rock approximation and Gassmann formula carry out tight sand
It is vertical, and calculating speed, density are realizing rock physics just with the conversion curve of rock porosity, water saturation and shale content
Drill calculating.
Dry rock in the content is approximately to consider parameter of pore structure on the basis of given Rock Matrix elasticity modulus
Method to calculate dry elastic modulus of rock, comes particular by porosity and parameter of pore structure to initial substrate elasticity modulus
It is modified, obtains rock and be in dry elastic modulus of rock of porosity when being φ, when pore structure is α.Concrete form see as
Lower formula (1):
K (φ)=K0(1-φ)p, μ (φ)=μ0(1-φ)q (1)
K0And μ0It is the elasticity modulus of matrix (skeleton), p and q are the functions of volumn concentration and pore structure,
It is that porosity is with μ (φ)When dry elastic modulus of rock.vlRepresent the volumn concentration of sandstone and mud stone, respectively vs
And vc, obtained by shale content.Tiijj(αl) and F (αl) it is parameter of pore structure.
Gassmann formula can be calculated the elasticity modulus of fluid saturated rocks by dry elastic modulus of rock, thus be calculated each
To the velocity of longitudinal wave and shear wave velocity of same sex medium, it is specifically shown in following formula (3), (4), (5):
ρsat=φ (Sw·ρw+(1-Sw·ρg))+(1-φ)·ρma (5)
Wherein, KfIt is full of the bulk modulus of fluid, K for holesatAnd μsatTo be full of fluid rock volume modulus and shearing
Modulus.
Wherein, SwFor water saturation.ρsatIt is the density of fluid saturated rocks, ρw、ρg、ρmaBe the density of water, gas it is close
The density of degree and rock matrix, is known.KdryAnd μdryIt is the bulk modulus and modulus of shearing of dry rock, Vρ estAnd Vs estPoint
It is not the velocity of longitudinal wave and shear wave velocity of isotropic medium.
Above-mentioned steps are established to be calculated respectively by given rock forming mineral component and parameter of pore structure, porosity, shale content
To the forward modeling process of rock speed in the same direction.
Next, shale content is set as definite value, give water saturation, can calculating speed and density with porosity change
Change curve.As one embodiment, speed when water saturation is 0,0.2,0.4,0.6,0.8,1, density are successively calculated with hole
The change curve of porosity；Water saturation is set as definite value, calculating and successively calculating shale content is 0,0.2,0.4,0.6,0.8,1
When speed, density with porosity change curve.Above-mentioned curve is converted, obtains p-wave impedance and speed ratio, then
Transformed curve is shown that in abscissa be p-wave impedance, ordinate is to obtain rock physics mould in the cross plot of speed ratio
Plate initial model.Wherein p-wave impedance is the product of velocity of longitudinal wave and density, and speed ratio is the ratio of velocity of longitudinal wave and shear wave velocity
Value.
Next, being optimized using actual measurement shear wave logging data to rock physics template initial model, after obtaining optimization
Rock physics template, see following formula (6):
In the formula (6),WithFor the longitudinal wave and shear wave velocity of actual measurement, VpAnd VsFor longitudinal wave and the shear wave speed of forward modeling
Degree.λ is shear wave weight coefficient, and value range is between 0 to 1.By adjusting the springform of sandstone and mud stone in rock constituents
Amount and hole ellipticity, so that the p-and s-wave velocity of forward modeling and the p-and s-wave velocity residual error of actual measurement are minimum.Sandstone adjusted at this time
With the elasticity modulus and hole ellipticity (α of mud stones,αc,Ks,μs,ρs,Kc,μc,ρc) be exactly optimization after rock physics template parameter,
Thus p-and s-wave velocity when rock physics forward modelling shale content, water saturation, porosity variation is carried out.
Finally, p-and s-wave velocity that prestack inversion obtains, density are converted, p-wave impedance and p-and s-wave velocity are obtained
Than p-wave impedance and speed ratio being shown in cross plot, the p-wave impedance and speed ratio for analyzing prestack inversion are in rock physics
Distribution in template analyzes oil-gas possibility by the porosity of inversion result region, water saturation.Wherein longitudinal wave hinders
Resist the product for velocity of longitudinal wave and density, speed ratio is the ratio of velocity of longitudinal wave and shear wave velocity.
Method of the invention establishes the rock physics template explained for Prestack seismic data based on petrophysical model,
By realizing the optimization to rock physics template compared with surveying shear wave velocity, by oil-bearing reservoir prediction from qualitative to fixed
Amount pushes away further.Explanation of the rock physics template established to instruct seismic inversion and inversion result realizes oil gas storage
The purpose of layer forecast of distribution and fluid detection.
Below by taking the real data in somewhere as an example, rock physics template is established using the method for the present invention, further relates to this
The effect of invention.
Fig. 2 is speed ratio-p-wave impedance rock physics template in embodiment 1.Fig. 3 is that the speed ratio-in embodiment 1 is vertical
Oily window in wave impedance rock physics template will enclose the data selected and be shown in slice of data is shown with particular color,
To sketch the contours of the distribution and form in especially potentiality region.
Table 1 is common elastic parameter, is obtained by petrophysics experiment.
| Lithology and elastic parameter | Numerical value |
| The bulk modulus of sandstone | 38GPa |
| The modulus of shearing of sandstone | 44GPa |
| The bulk modulus of mud stone | 20.9GPa |
| The modulus of shearing of mud stone | 6.85GPa |
| The bulk modulus of water | GPa |
The common elastic parameter of table 1
Embodiment 1
The present embodiment is the present invention for the practical logging data in somewhere and the embodiment of prestack inversion data, this hair
Bright method and step is as follows:
1, test analysis is carried out to tight sandstone reservoir rock sample first, obtains the input data of rock physics forward model.
Elasticity modulus, density and hole ellipticity including sandstone and mud stone.
2, the petrophysical model of tight sand is carried out using dry rock approximation and Gassmann formula for tight sand
It establishes, and calculating speed, density realize rock physics with the conversion curve of rock porosity, water saturation and shale content
Forward modelling.
3, on the basis of step 2, shale content is 1. set as definite value first, give water saturation, can calculating speed and
Density with porosity change curve.Speed when water saturation is 0,0.2,0.4,0.6,0.8,1, density are successively calculated with hole
The change curve of porosity.2. water saturation is set as definite value, calculate successively calculate shale content be 0,0.2,0.4,0.6,
0.8,1 when speed, density with porosity change curve.3. converting to above-mentioned curve, p-wave impedance and speed are obtained
Than (p-wave impedance obtains for the product of velocity of longitudinal wave and density, and speed ratio is that velocity of longitudinal wave is obtained divided by shear wave velocity), then
Transformed curve is shown that in abscissa be p-wave impedance, ordinate is in the cross plot of speed ratio.Obtain rock physics mould
Plate initial model.
4, rock physics template initial model is optimized using actual measurement shear wave logging data, the rock after being optimized
Physical template.In the case of recalculating different water cut saturation degree, shale content using the rock physics template parameter after optimization
Longitudinal wave, shear wave velocity and density.I.e. using the parameter after optimization, the action of above-mentioned steps 2 is carried out.
5, p-and s-wave velocity that prestack inversion obtains, density are converted, obtain p-wave impedance and P-S wave velocity ratio,
P-wave impedance and speed ratio are shown in cross plot, the p-wave impedance and speed ratio for analyzing prestack inversion are in rock physics template
In distribution, oil-gas possibility is analyzed by the porosity of inversion result region, water saturation.
Method of the invention can be used for the explanation of interpretation of logging data and pre-stack seismic inversion result.It is provided by the invention
Rock physics template method for building up can quickly be set up for logging well and the rock physics template of seismic interpretation, and logarithm is passed through
Oil-gas possibility is analyzed according to the distribution in rock physics template, can reduce the multi-solution of seismic data interpretation.
The presently disclosed embodiments is 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 technological improvement to the technology in market for best explaining each embodiment, or lead this technology
Other those of ordinary skill in domain can understand each embodiment disclosed herein.
Claims (6)
1. it is a kind of using rock physics template to the method for reservoir quantitative interpretation, comprising:
Rock sample test analysis；
Calculate rock physics template initial model；
Optimize rock physics template initial model；
Pre-stack seismic inversion result is explained using the rock physics template initial model of optimization；
Obtain oily window；
Wherein in rock sample test analysis step, by carrying out test analysis to reservoir rock sample, the mineralogical composition of composition rock is obtained
And parameter of pore structure；
Wherein in sand shale, including obtaining sandstone and the elasticity modulus of mud stone, density, hole ellipticity, the elasticity modulus
Including bulk modulus, modulus of shearing；
Wherein, the bulk modulus, modulus of shearing of sandstone, the symbol of density are expressed as Ks,μs,ρs,
The bulk modulus of mud stone, modulus of shearing, density symbol be expressed as Kc,μc,ρc；
Sandstone and the symbol of mud stone hole ellipticity are expressed as αs,αc；
Parameter (αs,αc,Ks,μs,ρs,Kc,μc,ρc) be rock physics forward model input data；
Wherein, the step of calculating rock physics template initial model include: for tight sand, using dry rock approximation and
Gassmann formula establishes the petrophysical model of tight sand, and calculating speed, density are with rock porosity, water saturation
And the conversion curve of shale content, realize rock physics forward modelling；
Shale content is wherein set as definite value, gives water saturation, calculating speed and density with the change curve of porosity；
Above-mentioned change curve is converted, p-wave impedance and speed ratio are obtained, then shows transformed change curve
Abscissa is p-wave impedance, and ordinate is to obtain rock physics template initial model in the cross plot of speed ratio, and wherein longitudinal wave hinders
Resist the product for velocity of longitudinal wave and density, speed ratio is the ratio of velocity of longitudinal wave and shear wave velocity.
2. it is according to claim 1 using rock physics template to the method for reservoir quantitative interpretation, wherein the dry rock
Stone is approximately to consider parameter of pore structure on the basis of the matrix elasticity modulus of given rock to calculate dry elastic modulus of rock
Method is modified initial substrate elasticity modulus particular by porosity and parameter of pore structure, obtains rock and be in
Dry elastic modulus of rock when porosity is φ, when pore structure is α:
K (φ)=K0(1-φ)p, μ (φ)=μ0(1-φ)q (1)
K0And μ0It is matrix elasticity modulus, p and q are the functions of volumn concentration and pore structure, and K (φ) and μ (φ) are holes
Dry elastic modulus of rock when degree is φ, vlRepresent the volumn concentration of sandstone and mud stone, respectively vsAnd vc, Tiijj(αl) and
F(αl) it is parameter of pore structure.
3. it is according to claim 2 using rock physics template to the method for reservoir quantitative interpretation, wherein described
Gassmann formula can be calculated the elasticity modulus of fluid saturated rocks by dry elastic modulus of rock, thus calculated isotropism and be situated between
The velocity of longitudinal wave and shear wave velocity of matter:
ρsat=φ (Sw·ρw+(1-Sw·ρg))+(1-φ)·ρma (5)
Wherein, KfThe bulk modulus of fluid, K are saturated for holesatAnd μsatFor fluid saturated rocks bulk modulus and modulus of shearing,
SwFor water saturation, ρsatIt is the density of fluid saturated rocks, ρw、ρg、ρmaIt is the density of water, the density of gas and Rock Matrix
Density, KdryAnd μdryIt is the bulk modulus and modulus of shearing of dry rock, Vρ estAnd Vs estIt is the longitudinal wave speed of isotropic medium respectively
Degree and shear wave velocity.
4. it is according to claim 1 using rock physics template to the method for reservoir quantitative interpretation, wherein horizontal using actual measurement
Wave well-log information optimizes rock physics template initial model, the rock physics template after being optimized:
Wherein,WithFor the longitudinal wave and shear wave velocity of actual measurement, VpAnd VsFor the longitudinal wave and shear wave velocity of forward modeling, λ is shear wave power
Weight coefficient, value range 0-1.
5. it is according to claim 1 using rock physics template to the method for reservoir quantitative interpretation, wherein including by prestack
P-and s-wave velocity that inverting obtains, density are converted, and p-wave impedance and P-S wave velocity ratio are obtained；It analyzes p-wave impedance and indulges
Shear wave velocity is than the distribution in rock physics template；Divided by the porosity of inversion result region, water saturation
Oil-gas possibility is analysed, oily window is obtained.
6. it is according to claim 4 using rock physics template to the method for reservoir quantitative interpretation, wherein including adjust rock
The elasticity modulus and hole ellipticity of sandstone and mud stone in stone component, so that the velocity of wave in length and breadth of the p-and s-wave velocity of forward modeling and actual measurement
Spend the rock physics template parameter that residual error is minimum, after the as optimization of the elasticity modulus and hole ellipticity of sandstone and mud stone at this time.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510618211.6A CN106556866B (en) | 2015-09-24 | 2015-09-24 | A method of using rock physics template to reservoir quantitative interpretation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510618211.6A CN106556866B (en) | 2015-09-24 | 2015-09-24 | A method of using rock physics template to reservoir quantitative interpretation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106556866A CN106556866A (en) | 2017-04-05 |
| CN106556866B true CN106556866B (en) | 2019-03-12 |
Family
ID=58415526
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510618211.6A Active CN106556866B (en) | 2015-09-24 | 2015-09-24 | A method of using rock physics template to reservoir quantitative interpretation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106556866B (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2562285B (en) * | 2017-05-11 | 2022-03-02 | Equinor Energy As | Fluid substitution |
| CN109521463B (en) * | 2017-09-20 | 2020-09-29 | 中国石油化工股份有限公司 | Method and system for determining optimal seismic excitation lithology of igneous rock near surface |
| CN110658555A (en) * | 2018-06-29 | 2020-01-07 | 中国石油化工股份有限公司 | Reservoir fluid factor construction method and reservoir fluid identification method |
| CN109115987B (en) * | 2018-07-20 | 2021-01-29 | 中国石油天然气股份有限公司 | Rock physical model-based fluid factor evaluation method and device |
| CN109709608B (en) * | 2018-12-07 | 2020-08-11 | 中国石油天然气股份有限公司 | Method for measuring TOC content and lithologic components of mixed accumulated hydrocarbon source rock and application |
| CN110426752B (en) * | 2019-08-20 | 2021-05-25 | 赛哲尔能源科技(北京)有限公司 | Reservoir parameter inversion method and system based on rock physical model |
| CN111856583B (en) * | 2020-06-12 | 2022-08-05 | 中国石油天然气股份有限公司 | Reservoir water saturation quantitative prediction method based on porosity constraint condition |
| CN112147689B (en) * | 2020-08-27 | 2024-01-30 | 中国石油天然气股份有限公司 | Carbonate reservoir interpretation method and interpretation volume setting method |
| CN112255688B (en) * | 2020-10-27 | 2022-08-02 | 中国海洋石油集团有限公司 | Method for inverting formation pressure by three-dimensional earthquake based on rock physics theory |
| CN112651171A (en) * | 2020-12-15 | 2021-04-13 | 清华大学 | Data-driven rock physical modeling method and system |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5265016A (en) * | 1989-07-19 | 1993-11-23 | Amoco Corporation | Method of shear wave velocity estimation |
| CN102353989A (en) * | 2011-08-24 | 2012-02-15 | 成都理工大学 | Method for estimating velocity of transverse waves based on inversion of equivalent elastic modulus for self-adapting matrix minerals |
| CN102854531A (en) * | 2012-09-11 | 2013-01-02 | 中国石油天然气股份有限公司 | Multi-scale rock physical charting method and device for detecting reservoir hydrocarbon |
| CN103135135A (en) * | 2013-02-26 | 2013-06-05 | 中国石油天然气股份有限公司 | Method and device used for hydrocarbon quantitative prediction and based on unconsolidated sandstone model |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2716196C (en) * | 2008-04-09 | 2016-11-01 | Exxonmobil Upstream Research Company | Method for generating anisotropic resistivity volumes from seismic and log data using a rock physics model |
-
2015
- 2015-09-24 CN CN201510618211.6A patent/CN106556866B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5265016A (en) * | 1989-07-19 | 1993-11-23 | Amoco Corporation | Method of shear wave velocity estimation |
| CN102353989A (en) * | 2011-08-24 | 2012-02-15 | 成都理工大学 | Method for estimating velocity of transverse waves based on inversion of equivalent elastic modulus for self-adapting matrix minerals |
| CN102854531A (en) * | 2012-09-11 | 2013-01-02 | 中国石油天然气股份有限公司 | Multi-scale rock physical charting method and device for detecting reservoir hydrocarbon |
| CN103135135A (en) * | 2013-02-26 | 2013-06-05 | 中国石油天然气股份有限公司 | Method and device used for hydrocarbon quantitative prediction and based on unconsolidated sandstone model |
Non-Patent Citations (1)
| Title |
|---|
| 基于Xu-White模型横波速度预测的误差分析;白俊雨 等;《地球物理学报》;20130228;第55卷(第2期);第589-594页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106556866A (en) | 2017-04-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106556866B (en) | A method of using rock physics template to reservoir quantitative interpretation | |
| CN111257926B (en) | Method for predicting ancient valley uranium reservoir by using old seismic data | |
| CN106054248A (en) | Earthquake rock physical inversion method based on large area tight reservoir | |
| Bodziak et al. | The role of seismic attributes in understanding the hydraulically fracturable limits and reservoir performance in shale reservoirs: An example from the Eagle Ford Shale, south Texas | |
| CN105093306B (en) | Reservoir is explained and thickness acquiring method automatically in a kind of geophysical exploration | |
| CN104570127B (en) | A kind of method of utilization seimic wave velocity Simultaneous Inversion porosity and shale content | |
| CN101236257A (en) | Oil well location determination technical method | |
| CN110275202B (en) | Method for predicting brittleness of compact oil reservoir | |
| Mahmood et al. | Total organic carbon content and total porosity estimation in unconventional resource play using integrated approach through seismic inversion and well logs analysis within the Talhar Shale, Pakistan | |
| CN104834003B (en) | Phase-controlled compression coefficient earthquake prediction method of unconventional tight gas reservoir | |
| CN106597543A (en) | Stratigraphic sedimentary facies division method | |
| CN105842733A (en) | Shale reservoir earthquake identification method | |
| CN109270578A (en) | Beach shallow sea law quality reservoirs rolling development method | |
| CN110456412B (en) | Method for identifying fluid saturation of carbonate reservoir based on post-stack seismic data | |
| CN107728205A (en) | A kind of Formation pressure prediction method | |
| Shi et al. | Characterization of hydraulic fracture configuration based on complex in situ stress field of a tight oil reservoir in Junggar Basin, Northwest China | |
| CN113393335A (en) | Reservoir oil and gas prediction method and device based on multi-seismic attribute optimization | |
| Gaddipati et al. | An Integrated Reservoir Modeling Case Study to Simulate Multi-Stage Hydraulically Fractured Horizontal Wells, based on Seismic, Petrophysical and Geological data for Pinedale Tight Gas Fluvial Reservoir | |
| Kimple et al. | Stratigraphy and Porosity Modeling of South-Central Illinois (USA) chester (upper mississippian) series sandstones using petrel | |
| Alam et al. | Synthetic core from conventional logs: A new method of interpretation for identification of key reservoir properties | |
| Adabnezhad et al. | Three-dimensional modeling of geomechanical units using acoustic impedance in one of the gas fields in South of Iran | |
| Siddiqui et al. | Facies characteristics and static reservoir connectivity of some siliciclastic tertiary outcrop successions in Bintulu and Miri, Sarawak, East Malaysia | |
| Clemons et al. | Seismic attributes: Exploiting seismic data to understand heterogeneous reservoir performance in the Eagle Ford Shale, south Texas, USA | |
| Yu | Facies-controlled modeling for permeability of tight gas reservoir based on hydrodynamic and geophysics characteristics | |
| CN109736791A (en) | Second interpretation technical method towards reservoir in a kind of oil and gas development |
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 |