CN106556866A - A kind of method of utilization rock physicses template to reservoir quantitative interpretation - Google Patents
A kind of method of utilization rock physicses template to reservoir quantitative interpretation Download PDFInfo
- Publication number
- CN106556866A CN106556866A CN201510618211.6A CN201510618211A CN106556866A CN 106556866 A CN106556866 A CN 106556866A CN 201510618211 A CN201510618211 A CN 201510618211A CN 106556866 A CN106556866 A CN 106556866A
- Authority
- CN
- China
- Prior art keywords
- rock
- template
- physicses
- wave
- 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.)
- Granted
Links
- 239000011435 rock Substances 0.000 title claims abstract description 169
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000004458 analytical method Methods 0.000 claims abstract description 15
- 238000009826 distribution Methods 0.000 claims abstract description 15
- 238000005457 optimization Methods 0.000 claims abstract description 15
- 238000012360 testing method Methods 0.000 claims abstract description 13
- 238000001615 p wave Methods 0.000 claims description 31
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 26
- 239000004575 stone Substances 0.000 claims description 22
- 239000012530 fluid Substances 0.000 claims description 19
- 239000011148 porous material Substances 0.000 claims description 19
- 241001274660 Modulus Species 0.000 claims description 14
- 239000004576 sand Substances 0.000 claims description 14
- 230000008859 change Effects 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 12
- 238000010008 shearing Methods 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 229920006395 saturated elastomer Polymers 0.000 claims description 8
- 238000013519 translation Methods 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 6
- 239000000758 substrate Substances 0.000 claims description 6
- 239000000470 constituent Substances 0.000 claims description 3
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 2
- 239000011707 mineral Substances 0.000 claims description 2
- 238000005553 drilling Methods 0.000 claims 1
- 239000000203 mixture Substances 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 239000007787 solid Substances 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/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
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
The present invention proposes a kind of method of utilization rock physicses template to reservoir quantitative interpretation, including：Rock sample test analysis；Calculate rock physicses template initial model；Optimization rock physicses template initial model；Pre-stack seismic inversion result is explained using the rock physicses template initial model of optimization；Obtain oily window.The method of the present invention can be used for the explanation of interpretation of logging data and pre-stack seismic inversion result.The rock physicses template method for building up that the present invention is provided can quickly set up the rock physicses template for well logging and seismic interpretation, analyze oil-gas possibility by the distribution to data in rock physicses template, can reduce the multi-solution of seismic data interpretation.
Description
Technical field
The present invention relates to petrophysics and geophysical well logging technology field, and in particular to one kind utilizes rock
Method of the physical template to reservoir quantitative interpretation.
Background technology
Petrophysics is the bridge for connecting petrophysical parameter, fluid properties parameter and seismic elastic parameter,
Quantitative interpretation for geological data provides solid foundation, greatly reduces the multi-solution of seismic interpretation.In recent years
Come, significant progress is achieved in terms of earthquake rock physicses basic research both at home and abroad, while also more successful
Ground achievement in research is applied to into the aspects such as Lithology Discrimination, fluid detection and porosity inversion, seismic prospecting,
The fields such as oil reservoir development played an important role.
With deepening continuously for oil-gas exploration and development, oil gas reservoir body becomes increasingly complex, and exploration emphasis is by original
Structural deposit to lithologic deposit, subtle pool shift.One side exploration targets is gradually deepened, is become
Little, difficulties in exploration is increased, and the required precision to reservoir prediction is higher.On the other hand, reservoir prediction is final
Purpose is except judging the existence of reservoir, predicting reservoir development degree ground quality, while in needing to judge reservoir
The property of contained fluid.
Substantial amounts of to scientific investigations showed that, China possesses abundant Stratigraphic and subtle reservoirs resource, has huge surveying
Visit potentiality to be exploited.And seismic inversion is the important technical for obtaining lithologic character stratum elastic parameter, by poststack
Inverting can obtain p-wave impedance data volume, by prestack inversion can obtain p-wave impedance, S-wave impedance,
The parameter such as density or velocity of longitudinal wave, shear wave velocity, density.These elastic parameters are widely used in oil-containing
Gas reservoir predicts, achieves many prediction effects well, but also has the example of failure, does not reach differentiation rock
Property and identification fluid purpose.In terms of the reservoir information that contained of seismic inversion data is excavated, have scholar and
Scientific research personnel has done substantial amounts of research and production application work, the such as fluid identification such as Poisson impedance, speed ratio
The appearance of the factor.Explanation to these elastic parameters and its conversion at present is carried out in qualitatively degree,
The multi-solution that a kind of parameter of single application carries out reservoir prediction work is stronger, brings to the explanation of seismic data
Risk, limits the application of seismic inversion and explains work.
The content of the invention
Explanation means and the shortcoming being of limited application for seismic inversion, the inventive method are based on rock
Stone physical model, establishes the rock physicses template for seismic data quantitative interpretation.By multiparameter and rock
Physical template carries out intersection analysis, analyzes the distribution of each parameter in rock physicses template, and by purpose
The elastic parameter in region is sketched the contours in section, section or three dimensions, realizes the explanation to target geological body
And prediction work, reduce explaining the multi-solution of work.
The present invention provides a kind of method of utilization rock physicses template to reservoir quantitative interpretation, for instructing earthquake
Inverting and the explanation of inversion result, realize oil and gas reservoir forecast of distribution and fluid detection.The method of the present invention,
Including：Rock sample test analysis；Calculate rock physicses template initial model；Optimization rock physicses template introductory die
Type；Pre-stack seismic inversion result is explained using the rock physicses template initial model of optimization；Contained
Oil gas window.
Further, in rock sample test analysis step, by carrying out test analysis to reservoir rock sample, obtain
The mineralogical composition and parameter of pore structure of composition rock.
Alternatively, in sand shale, elastic modelling quantity, density including acquisition sandstone and mud stone, hole ellipticity,
Described elastic modelling quantity includes bulk moduluses, modulus of shearing；Wherein the elastic modelling quantity of sandstone, density are symbol
It is expressed as Ks,μs,ρs, the elastic modelling quantity of mud stone, density are Kc,μc,ρc；Hole ellipticity is αs,αc；Parameter
(αs,αc,Ks,μs,ρs,Kc,μc,ρc) for the input data of rock physicses forward model.
Alternatively, rock physicses template initial model is calculated, including approximate using dry rock for tight sand
The petrophysical model of tight sand, and calculating speed, density are set up with Gassmann formula with blowhole
The conversion curve of degree, water saturation and shale content, realizes rock physicses forward modelling.
Further, described dry rock is approximately on the basis of given Rock Matrix elastic modelling quantity, it is considered to hole
Gap structure parameter come the method that calculates dry elastic modulus of rock, particular by porosity and parameter of pore structure
To be modified to initial substrate elastic modelling quantity, obtain rock in porosity be φ when, pore structure be α when
Dry elastic modulus of rock：
K (φ)=K0(1-φ)p, μ (φ)=μ0(1-φ)q (1)
K0And μ0It is the elastic modelling quantity of substrate, p and q is the function of volumn concentration and pore structure, K (φ)
With μ (φ) be porosity be φ when dry elastic modulus of rock, vlThe volumn concentration of sandstone and mud stone is represented,
Respectively vsAnd vc, Tiijj(αl) and F (αl) it is parameter of pore structure.
Alternatively, described Gassmann formula can calculate fluid saturated rocks by dry elastic modulus of rock
Elastic modelling quantity, thus calculates the velocity of longitudinal wave and shear wave velocity of isotropic medium：
ρsat=φ (Sw·ρw+(1-Sw·ρg))+(1-φ)·ρma (5)
Wherein, KfFor the bulk moduluses of hole saturation fluid, KsatAnd μsatFor fluid saturated rocks bulk moduluses
And modulus of shearing, SwFor water saturation, ρsatIt is the density of fluid saturated rocks, ρw、ρg、ρmaIt is water
Density, the density of gas and Rock Matrix density, KdryAnd μdryIt is the bulk moduluses and shearing mould of dry rock
Amount, Vρ estAnd Vs estIt is the velocity of longitudinal wave and shear wave velocity of isotropic medium respectively.
Further, shale content is set to into definite value, gives water saturation, calculating speed and density are with hole
The change curve of porosity；Line translation is entered to above-mentioned curve, p-wave impedance and speed ratio is obtained, then will conversion
It is p-wave impedance that curve afterwards is displayed in abscissa, during vertical coordinate is for the cross plot of speed ratio, obtains rock thing
Reason template initial model, wherein p-wave impedance are the product of velocity of longitudinal wave and density, and speed ratio is velocity of longitudinal wave
With the ratio of shear wave velocity.
Further, rock physicses template initial model is optimized using actual measurement shear wave logging data, is obtained
Rock physicses template to after optimization：
Wherein,WithFor the compressional wave and shear wave velocity of actual measurement, VpAnd VsFor the compressional wave just drilled and shear wave velocity,
λ is shear wave weight coefficient, and span is between 0 to 1.Wherein by adjust rock constituents in sandstone and
The elastic modelling quantity and hole ellipticity of mud stone so that the p-and s-wave velocity just drilled and the p-and s-wave velocity residual error surveyed
Minimum, the elastic modelling quantity and hole ellipticity (α of sandstone and mud stone after now adjustings,αc,Ks,μs,ρs,Kc,μc,ρc)
It is exactly the rock physicses template parameter after optimizing.
Further, the p-and s-wave velocity that prestack inversion obtained, density enter line translation, obtain p-wave impedance
And P-S wave velocity ratio, the distribution of the p-wave impedance and speed ratio of prestack inversion in rock physicses template is analyzed,
Oil-gas possibility is analyzed by the porosity of inversion result region, water saturation, oil-containing transom window is obtained
Mouthful.Wherein p-wave impedance is the product of velocity of longitudinal wave and density, and speed ratio is velocity of longitudinal wave and shear wave velocity
Ratio.
The method of the present invention can be used for the explanation of interpretation of logging data and pre-stack seismic inversion result.The present invention
The rock physicses template method for building up of offer can quickly set up the rock thing for well logging and seismic interpretation
Reason template, analyzes oil-gas possibility by the distribution to data in rock physicses template, can reduce earthquake
The multi-solution of data interpretation.
Description of the drawings
Disclosure illustrative embodiments are described in more detail by combining accompanying drawing, the disclosure it is above-mentioned
And other purposes, feature and advantage will be apparent from, wherein, in disclosure illustrative embodiments
In, identical reference number typically represents same parts.
Fig. 1 is the flow chart of the present invention.
Fig. 2 is the speed ratio-p-wave impedance rock physicses template in embodiment 1.
Fig. 3 is the oily window in the speed ratio-p-wave impedance rock physicses template in embodiment 1, will circle
The data of choosing are shown with particular color in slice of data shows, so as to sketch the contours of the distribution in especially potentiality region
Scope and form.
Specific embodiment
The preferred implementation of the disclosure is more fully described below with reference to accompanying drawings.Although showing in accompanying drawing
The preferred implementation of the disclosure, however, it is to be appreciated that may be realized in various forms the disclosure and should not be by
Embodiments set forth herein is limited.Conversely, thesing embodiments are provided so that the disclosure is more saturating
It is thorough and complete, and the scope of the present disclosure intactly can be conveyed to those skilled in the art.
The invention provides a kind of method of utilization rock physicses template to reservoir quantitative interpretation, can be used to log well
The explanation of data interpretation and pre-stack seismic inversion result.The rock physicses template method for building up that the present invention is provided
Quickly can set up for well logging and seismic interpretation rock physicses template, by data in rock thing
Distribution in reason template can reduce the multi-solution of seismic data interpretation analyzing oil-gas possibility.
In order to set up the rock physicses template for Sandstone Gas Reservoir quantitative interpretation, the present invention is first to reservoir
Rock sample carries out test analysis and obtains constituting the mineralogical composition and parameter of pore structure of rock, then using rock thing
Reason model calculates Rock Elastic Parameters with the change of porosity, water saturation and shale content, obtains rock
Physical template initial model, recycles actual measurement shear wave logging data to carry out rock physicses template initial model excellent
Change, the rock physicses template after being optimized, the same rock physicses of the elastic parameter that finally prestack inversion is obtained
Template is simultaneously displayed in cross plot, by distribution of the elastic parameter of prestack inversion in rock physicses 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 of utilization rock physicses template to reservoir quantitative interpretation,
Belong to petrophysics and geophysical well logging technology field.Methods described includes：Rock sample test analysis；Meter
Calculate rock physicses template initial model；Optimization rock physicses template initial model；Using the rock physicses of optimization
Template initial model is explained to pre-stack seismic inversion result；Obtain oily window.
Further, the method for the present invention includes：Carry out test analysis to obtain constituting rock to reservoir rock sample
Mineralogical composition and parameter of pore structure；Using petrophysical model calculate Rock Elastic Parameters with porosity, contain
The change of water saturation and shale content, obtains rock physicses template initial model；Using actual measurement SWAL
Data is optimized to rock physicses template initial model, the rock physicses template after being optimized；By prestack
The elastic parameter that inverting is obtained is simultaneously displayed in cross plot with rock physicses template, by the bullet of prestack inversion
Property distribution of the parameter in rock physicses template is carrying out quantitative interpretation to reservoir and oil-gas possibility.
For ease of understanding the scheme and its effect of the embodiment of the present invention, a concrete application example given below.
It will be understood by those skilled in the art that the example is only for the purposes of understanding the present invention, its any detail is not
It is intended to limit the present invention by any way.
Used as an embodiment, the method for the present invention is comprised the following steps：
Test analysis are carried out to rock sample first.Carry out test analysis to obtain constituting rock to Sandstone Gas Reservoir rock sample
The mineralogical composition and parameter of pore structure of stone.In sand shale, elastic modelling quantity including sandstone and mud stone and close
The parameters such as the hole ellipticity of degree, sandstone and mud stone.Described elastic modelling quantity includes bulk moduluses, modulus of shearing；
The elastic modelling quantity of described sandstone and mud stone, density and hole ellipticity are the defeated of follow-up rock physicses forward model
Enter data.The elastic modelling quantity of sandstone, density are expressed as K for symbols,μs,ρs,；The elastic modelling quantity of mud stone, density
For Kc,μc,ρc；Hole ellipticity is αs,αc.This group of parameter (αs,αc,Ks,μs,ρs,Kc,μc,ρc) for rock physicses just
Drill the input data of model.
For tight sand using dry rock is approximate and Gassmann formula carry out the rock physicses mould of tight sand
Type is set up, and calculating speed, density are with the conversion curve of rock porosity, water saturation and shale content,
Realize rock physicses forward modelling.
Dry rock in the content is approximately on the basis of given Rock Matrix elastic modelling quantity, it is considered to which hole is tied
Structure parameter comes right particular by porosity and parameter of pore structure come the method for calculating dry elastic modulus of rock
Initial substrate elastic modelling quantity is modified, obtain rock in porosity be φ when, pore structure be α when it is dry
Elastic modulus of rock.Concrete form is shown in equation below (1)：
K (φ)=K0(1-φ)p, μ (φ)=μ0(1-φ)q (1)
K0And μ0It is the elastic modelling quantity of substrate (skeleton), p and q is the letter of volumn concentration and pore structure
Number,It is that porosity is with μ (φ)When dry elastic modulus of rock.vlRepresent the volume hundred of sandstone and mud stone
Divide content, respectively vsAnd vc, obtained by shale content.Tiijj(αl) and F (αl) it is parameter of pore structure.
Gassmann formula can be calculated the elastic modelling quantity of fluid saturated rocks by dry elastic modulus of rock, thus
The velocity of longitudinal wave and shear wave velocity of isotropic medium are calculated, equation below (3), (4), (5) are specifically shown in：
ρsat=φ (Sw·ρw+(1-Sw·ρg))+(1-φ)·ρma (5)
Wherein, KfIt is full of the bulk moduluses of fluid, K for holesatAnd μsatTo be full of fluid rock volume modulus
And modulus of shearing.
Wherein, SwFor water saturation.ρsatIt is the density of fluid saturated rocks, ρw、ρg、ρmaIt is water
The density of density, the density of gas and rock matrix, is known.KdryAnd μdryIt is the bulk moduluses of dry rock
And modulus of shearing, Vρ estAnd Vs estIt is the velocity of longitudinal wave and shear wave velocity of isotropic medium respectively.
Above-mentioned steps are established by given rock forming mineral component and parameter of pore structure, porosity, shale content
Calculate and each just drill process to rock speed in the same direction.
Next, shale content is set to definite value, give water saturation, can calculating speed and density with hole
The change curve of porosity.As one embodiment, calculate successively water saturation be 0,0.2,0.4,0.6,
0.8th, 1 when speed, density with porosity change curve；Water saturation is set to into definite value, is calculated successively
Change curve of the speed, density when calculating shale content is 0,0.2,0.4,0.6,0.8,1 with porosity.
Line translation is entered to above-mentioned curve, p-wave impedance and speed ratio is obtained, is then included the curve after conversion in horizontal stroke
Coordinate is p-wave impedance, during vertical coordinate is for the cross plot of speed ratio, obtains rock physicses template 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.
Next, being optimized to rock physicses template initial model using actual measurement shear wave logging data, obtain
Rock physicses template after optimization, is shown in equation below (6)：
In the formula (6),WithFor the compressional wave and shear wave velocity of actual measurement, VpAnd VsFor the compressional wave just drilled
And shear wave velocity.λ is shear wave weight coefficient, and span is between 0 to 1.By adjusting in rock constituents
Sandstone and mud stone elastic modelling quantity and hole ellipticity so that the ripple in length and breadth of the p-and s-wave velocity just drilled and actual measurement
Speed residual error is minimum.The elastic modelling quantity and hole ellipticity of sandstone and mud stone after now adjusting
(αs,αc,Ks,μs,ρs,Kc,μc,ρc) be exactly optimization after rock physicses template parameter, thus just carrying out rock physicses
Drill p-and s-wave velocity when calculating shale content, water saturation, porosity change.
Finally, the p-and s-wave velocity that prestack inversion obtained, density enter line translation, obtain p-wave impedance and vertical
Shear wave velocity ratio, p-wave impedance and speed ratio are shown in cross plot, analyze the p-wave impedance of prestack inversion
With distribution of the speed ratio in rock physicses template, by the porosity of inversion result region, aqueous full
Oil-gas possibility is analyzed with degree.Wherein p-wave impedance is the product of velocity of longitudinal wave and density, and speed ratio is compressional wave
The ratio of speed and shear wave velocity.
The method of the present invention establishes the rock physicses explained for Prestack seismic data based on petrophysical model
Template, by the optimization to rock physicses template is realized with the comparison of actual measurement shear wave velocity, oily is stored up
Layer prediction is from qualitative to one step of fixed quanting pushing.The rock physicses template set up is to instruct seismic inversion and anti-
The explanation of result is drilled, the purpose of oil and gas reservoir forecast of distribution and fluid detection is realized.
Below by taking the real data in somewhere as an example, rock physicses template is set up using the inventive method, and then
Illustrate the effect of the present invention.
Fig. 2 is the speed ratio-p-wave impedance rock physicses template in embodiment 1.During Fig. 3 is embodiment 1
Oily window in speed ratio-p-wave impedance rock physicses template, the data of circle choosing are shown in slice of data
Middle particular color shows, so as to sketch the contours of the distribution and form in especially potentiality region.
Table 1 is conventional elastic parameter, is obtained by petrophysics experiment.
| Lithology and elastic parameter | Numerical value |
| The bulk moduluses of sandstone | 38GPa |
| The modulus of shearing of sandstone | 44GPa |
| The bulk moduluses of mud stone | 20.9GPa |
| The modulus of shearing of mud stone | 6.85GPa |
| The bulk moduluses of water | GPa |
The conventional elastic parameter of table 1
Embodiment 1
The present embodiment is enforcement reality of the present invention for the practical logging data and prestack inversion data in somewhere
Example, method of the present invention step are as follows：
1st, test analysis are carried out to Sandstone Gas Reservoir rock sample first, obtains the input of rock physicses forward model
Data.Elastic modelling quantity, density and hole ellipticity including sandstone and mud stone.
2nd, for tight sand using dry rock is approximate and Gassmann formula carry out the rock thing of tight sand
Reason model is set up, and calculating speed, density are with the conversion of rock porosity, water saturation and shale content
Curve, realizes rock physicses forward modelling.
3rd, on the basis of step 2, shale content is set to into definite value first 1., gives water saturation, can count
Calculate the change curve of speed and density with porosity.Successively calculate water saturation be 0,0.2,0.4,0.6,
0.8th, 1 when speed, density with porosity change curve.2. water saturation is set to into definite value, calculate according to
The secondary speed calculated when shale content is 0,0.2,0.4,0.6,0.8,1, density are bent with the change of porosity
Line.3. line translation is entered to above-mentioned curve, obtain p-wave impedance and speed ratio (p-wave impedance be velocity of longitudinal wave with
The product of density is obtained, and speed ratio is obtained divided by shear wave velocity for velocity of longitudinal wave), then by the curve after conversion
It is p-wave impedance to be displayed in abscissa, and vertical coordinate is in the cross plot of speed ratio.Obtain at the beginning of rock physicses template
Beginning model.
4th, rock physicses template initial model is optimized using actual measurement shear wave logging data, after being optimized
Rock physicses template.Using the rock physicses template parameter after optimization recalculate different water cut saturation,
Compressional wave, shear wave velocity and density in the case of shale content.I.e. using the parameter after optimization, above-mentioned step is carried out
Rapid 2 job content.
5th, the p-and s-wave velocity that obtains prestack inversion, density enter line translation, obtain p-wave impedance and ripple in length and breadth
Speed ratio, p-wave impedance and speed ratio are shown in cross plot, analyze the p-wave impedance and speed of prestack inversion
Distribution of the degree ratio in rock physicses template, by the porosity of inversion result region, water saturation
To analyze oil-gas possibility.
The method of the present invention can be used for the explanation of interpretation of logging data and pre-stack seismic inversion result.The present invention
The rock physicses template method for building up of offer can quickly set up the rock thing for well logging and seismic interpretation
Reason template, analyzes oil-gas possibility by the distribution to data in rock physicses template, can reduce earthquake
The multi-solution of data interpretation.
It is described above the presently disclosed embodiments, described above is exemplary, and non-exclusive,
And it is also not necessarily limited to disclosed each embodiment.In the scope and spirit without departing from illustrated each embodiment
In the case of, many modifications and changes will be apparent from for those skilled in the art.
The selection of term used herein, it is intended to best explain the principle of each embodiment, practical application or to market
In technology technological improvement, or other those of ordinary skill of the art is understood that
Each embodiment.
Claims (10)
1. method of a kind of utilization rock physicses template to reservoir quantitative interpretation, including：
Rock sample test analysis；
Calculate rock physicses template initial model；
Optimization rock physicses template initial model；
Pre-stack seismic inversion result is explained using the rock physicses template initial model of optimization；
Obtain oily window.
2. method of the utilization rock physicses template according to claim 1 to reservoir quantitative interpretation, wherein
In rock sample test analysis step, by carrying out test analysis to reservoir rock sample, obtain constitute rock mineral into
Divide and parameter of pore structure.
3. method of the utilization rock physicses template according to claim 2 to reservoir quantitative interpretation, wherein
In sand shale, elastic modelling quantity, density including acquisition sandstone and mud stone, hole ellipticity, described elasticity
Modulus includes bulk moduluses, modulus of shearing；
Wherein the elastic modelling quantity of sandstone, density are expressed as K for symbols,μs,ρs,
The elastic modelling quantity of mud stone, density are Kc,μc,ρc；
Hole ellipticity is αs,αc；
Parameter (αs,αc,Ks,μs,ρs,Kc,μc,ρc) for the input data of rock physicses forward model.
4. method of the utilization rock physicses template according to claim 3 to reservoir quantitative interpretation, wherein
Rock physicses template initial model is calculated, dry rock is approximate and Gassmann is public including utilizing for tight sand
Formula sets up the petrophysical model of tight sand, and calculating speed, density with rock porosity, containing water saturation
The conversion curve of degree and shale content, realizes rock physicses forward modelling.
5. method of the utilization rock physicses template according to claim 4 to reservoir quantitative interpretation, wherein
Described dry rock is approximately on the basis of given Rock Matrix elastic modelling quantity, it is considered to which parameter of pore structure is counting
The method for calculating dry elastic modulus of rock, comes to initial substrate bullet particular by porosity and parameter of pore structure
Property modulus is modified, and obtaining rock in porosity isWhen, pore structure be α when dry rock elasticity mould
Amount：
K (φ)=K0(1-φ)p, μ (φ)=μ0(1-φ)q (1)
K0And μ0It is the elastic modelling quantity of substrate, p and q is the function of volumn concentration and pore structure,
It is that porosity is with μ (φ)When dry elastic modulus of rock, vlThe volumn concentration of sandstone and mud stone is represented,
Respectively vsAnd vc, Tiijj(αl) and F (αl) it is parameter of pore structure.
6. method of the utilization rock physicses template according to claim 5 to reservoir quantitative interpretation, wherein
Described Gassmann formula can be calculated the elastic modelling quantity of fluid saturated rocks by dry elastic modulus of rock, by
This calculates the velocity of longitudinal wave and shear wave velocity of isotropic medium：
ρsat=φ (Sw·ρw+(1-Sw·ρg))+(1-φ)·ρma (5)
Wherein, KfFor the bulk moduluses of hole saturation fluid, KsatAnd μsatFor fluid saturated rocks bulk moduluses
And modulus of shearing, SwFor water saturation, ρsatIt is the density of fluid saturated rocks, ρw、ρg、ρmaIt is water
Density, the density of gas and Rock Matrix density, KdryAnd μdryIt is the bulk moduluses and shearing mould of dry rock
Amount, Vρ estAnd Vs estIt is the velocity of longitudinal wave and shear wave velocity of isotropic medium respectively.
7. method of the utilization rock physicses template according to claim 6 to reservoir quantitative interpretation, wherein
Shale content is set to into definite value, water saturation, the change curve of calculating speed and density with porosity is given；
Line translation is entered to above-mentioned curve, p-wave impedance and speed ratio is obtained, then the curve after conversion is shown
It is p-wave impedance in abscissa, during vertical coordinate is for the cross plot of speed ratio, obtains rock physicses template introductory die
Type, wherein p-wave impedance are the product of velocity of longitudinal wave and density, and speed ratio is velocity of longitudinal wave and shear wave velocity
Ratio.
8. method of the utilization rock physicses template according to claim 1 to reservoir quantitative interpretation, wherein
Rock physicses template initial model is optimized using actual measurement shear wave logging data, the rock after being optimized
Physical template：
Wherein,WithFor the compressional wave and shear wave velocity of actual measurement, VpAnd VsFor the compressional wave just drilled and shear wave velocity,
λ is shear wave weight coefficient, and span is between 0 to 1.
9. method of the utilization rock physicses template according to claim 1 to reservoir quantitative interpretation, wherein
P-and s-wave velocity that prestack inversion is obtained, density enter line translation, obtain p-wave impedance and P-S wave velocity ratio,
Distribution of the p-wave impedance and speed ratio of analysis prestack inversion in rock physicses template, by inversion result institute
Porosity, water saturation in region obtains oily window analyzing oil-gas possibility；
Wherein p-wave impedance is the product of velocity of longitudinal wave and density, and speed ratio is velocity of longitudinal wave and shear wave velocity
Ratio.
10. method of the utilization rock physicses template according to claim 8 to reservoir quantitative interpretation, its
In by adjust rock constituents in sandstone and mud stone elastic modelling quantity and hole ellipticity so that just drilling in length and breadth
Wave velocity is minimum with the p-and s-wave velocity residual error of actual measurement, the elastic modelling quantity of sandstone and mud stone after now adjusting and
Hole ellipticity (αs,αc,Ks,μs,ρs,Kc,μc,ρc) be exactly optimization after rock physicses template parameter.
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 true CN106556866A (en) | 2017-04-05 |
| CN106556866B 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) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109115987A (en) * | 2018-07-20 | 2019-01-01 | 中国石油天然气股份有限公司 | A kind of evaluation method and device of the fluid factor based on petrophysical model |
| CN109521463A (en) * | 2017-09-20 | 2019-03-26 | 中国石油化工股份有限公司 | Determine the method and system of the best earthquake excitation li thology of igneous rock near surface |
| CN109709608A (en) * | 2018-12-07 | 2019-05-03 | 中国石油天然气股份有限公司 | The mixed product hydrocarbon source rock TOC content of one kind and lithology compound mensuration method and application |
| CN110426752A (en) * | 2019-08-20 | 2019-11-08 | 赛哲尔能源科技(北京)有限公司 | A kind of reservoir parameter inversion method and system based on petrophysical model |
| CN110658555A (en) * | 2018-06-29 | 2020-01-07 | 中国石油化工股份有限公司 | Reservoir fluid factor construction method and reservoir fluid identification method |
| CN110945384A (en) * | 2017-05-11 | 2020-03-31 | 艾奎诺能源公司 | Fluid displacement |
| CN111856583A (en) * | 2020-06-12 | 2020-10-30 | 中国石油天然气股份有限公司 | Reservoir water saturation quantitative prediction method based on porosity constraint condition |
| CN112147689A (en) * | 2020-08-27 | 2020-12-29 | 中国石油天然气股份有限公司 | Carbonate reservoir interpretation method and interpretation quantity version establishment method |
| CN112255688A (en) * | 2020-10-27 | 2021-01-22 | 中国海洋石油集团有限公司 | 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 (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5265016A (en) * | 1989-07-19 | 1993-11-23 | Amoco Corporation | Method of shear wave velocity estimation |
| US20100326669A1 (en) * | 2008-04-09 | 2010-12-30 | Yaping Zhu | Method for generating anisotropic resistivity volumes from seismic and log data using a rock physics model |
| 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 |
-
2015
- 2015-09-24 CN CN201510618211.6A patent/CN106556866B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5265016A (en) * | 1989-07-19 | 1993-11-23 | Amoco Corporation | Method of shear wave velocity estimation |
| US20100326669A1 (en) * | 2008-04-09 | 2010-12-30 | Yaping Zhu | Method for generating anisotropic resistivity volumes from seismic and log data using a rock physics model |
| 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模型横波速度预测的误差分析", 《地球物理学报》 * |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110945384B (en) * | 2017-05-11 | 2022-04-05 | 艾奎诺能源公司 | Fluid displacement |
| CN110945384A (en) * | 2017-05-11 | 2020-03-31 | 艾奎诺能源公司 | Fluid displacement |
| CN109521463A (en) * | 2017-09-20 | 2019-03-26 | 中国石油化工股份有限公司 | Determine the method and system of the best earthquake excitation li thology of igneous rock near surface |
| CN110658555A (en) * | 2018-06-29 | 2020-01-07 | 中国石油化工股份有限公司 | Reservoir fluid factor construction method and reservoir fluid identification method |
| CN109115987A (en) * | 2018-07-20 | 2019-01-01 | 中国石油天然气股份有限公司 | A kind of evaluation method and device of the fluid factor based on petrophysical model |
| CN109709608A (en) * | 2018-12-07 | 2019-05-03 | 中国石油天然气股份有限公司 | The mixed product hydrocarbon source rock TOC content of one kind and lithology compound mensuration method and application |
| CN110426752A (en) * | 2019-08-20 | 2019-11-08 | 赛哲尔能源科技(北京)有限公司 | A kind of reservoir parameter inversion method and system based on petrophysical model |
| CN111856583A (en) * | 2020-06-12 | 2020-10-30 | 中国石油天然气股份有限公司 | Reservoir water saturation quantitative prediction method based on porosity constraint condition |
| CN111856583B (en) * | 2020-06-12 | 2022-08-05 | 中国石油天然气股份有限公司 | Reservoir water saturation quantitative prediction method based on porosity constraint condition |
| CN112147689A (en) * | 2020-08-27 | 2020-12-29 | 中国石油天然气股份有限公司 | Carbonate reservoir interpretation method and interpretation quantity version establishment method |
| CN112147689B (en) * | 2020-08-27 | 2024-01-30 | 中国石油天然气股份有限公司 | Carbonate reservoir interpretation method and interpretation volume setting method |
| CN112255688A (en) * | 2020-10-27 | 2021-01-22 | 中国海洋石油集团有限公司 | 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 |
| CN112651171B (en) * | 2020-12-15 | 2024-05-14 | 清华大学 | Data-driven rock physical modeling method and system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106556866B (en) | 2019-03-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106556866A (en) | A kind of method of utilization rock physicses template to reservoir quantitative interpretation | |
| CN104950331B (en) | A kind of porosity of sand mud reservoir and the earthquake prediction method of shale content | |
| Radwan | Three-dimensional gas property geological modeling and simulation | |
| 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 | |
| Ali et al. | 3D integrated structural, facies and petrophysical static modeling approach for complex sandstone reservoirs: A case study from the coniacian–santonian matulla formation, july oilfield, gulf of suez, Egypt | |
| CN105527653B (en) | A kind of virtual log construction method based on geological information | |
| CN105469159A (en) | Method capable of realizing quantitative prediction on favorable oil gas accumulation area | |
| CN101236257A (en) | Oil well location determination technical method | |
| CN102906600B (en) | Method of prospecting directly for free gas in stratum | |
| CN103675907A (en) | AVO inversion hydrocarbon detection method based on petrographic constraints | |
| Kianoush et al. | Compressional and shear interval velocity modeling to determine formation pressures in an oilfield of SW Iran | |
| CN107515957A (en) | Mud shale SEQUENCE STRATIGRAPHIC DIVISION method | |
| Ali et al. | Prediction of Cretaceous reservoir zone through petrophysical modeling: Insights from Kadanwari gas field, Middle Indus Basin | |
| Ganguli | Integrated reservoir studies for CO2-enhanced oil recovery and sequestration: Application to an Indian mature oil field | |
| Kianoush et al. | Application of pressure-volume (PV) fractal models in modeling formation pressure and drilling fluid determination in an oilfield of SW Iran | |
| Khadem et al. | Integration of rock physics and seismic inversion for rock typing and flow unit analysis: a case study | |
| CN112505754B (en) | Method for collaborative partitioning sedimentary microfacies by well-seismic based on high-precision sequence grid model | |
| Guo et al. | Integrated geophysical technologies for unconventional reservoirs and case study within Fuling shale gas field, Sichuan Basin, China | |
| CN109283574B (en) | Low-frequency model construction method and computer-readable storage medium | |
| CN114910953B (en) | Method for explaining irregular section of complex fracture zone | |
| Mabrey | Rock quality index for Niobrara horizontal well drilling and completion optimization, Wattenberg field, Colorado | |
| CN103278852B (en) | Utilize the method for seismic data volume waveform configuration characteristic model predicting oil/gas | |
| Naji et al. | Prediction of sonic shear wave using artificial neural network | |
| Kimple et al. | Stratigraphy and Porosity Modeling of South-Central Illinois (USA) chester (upper mississippian) series sandstones using petrel | |
| Gunnarsson | 3D modeling in Petrel of geological CO2 storage site |
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 |