CN108761529B - Method and device for determining reservoir parameters based on three-dimensional elastic-electric rock physical quantity plate - Google Patents
Method and device for determining reservoir parameters based on three-dimensional elastic-electric rock physical quantity plate Download PDFInfo
- Publication number
- CN108761529B CN108761529B CN201810412688.2A CN201810412688A CN108761529B CN 108761529 B CN108761529 B CN 108761529B CN 201810412688 A CN201810412688 A CN 201810412688A CN 108761529 B CN108761529 B CN 108761529B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- parameter
- data
- reservoir
- rock
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 192
- 238000000034 method Methods 0.000 title claims abstract description 74
- 238000012545 processing Methods 0.000 claims abstract description 24
- 230000005611 electricity Effects 0.000 claims description 120
- 230000015572 biosynthetic process Effects 0.000 claims description 49
- 238000010008 shearing Methods 0.000 claims description 18
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 16
- 238000003860 storage Methods 0.000 claims description 15
- 239000002245 particle Substances 0.000 claims description 13
- 238000012937 correction Methods 0.000 claims description 10
- 239000012530 fluid Substances 0.000 claims description 9
- 239000004575 stone Substances 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 239000004927 clay Substances 0.000 claims description 5
- 229920006395 saturated elastomer Polymers 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 3
- 239000008398 formation water Substances 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 239000011148 porous material Substances 0.000 claims description 3
- 238000010206 sensitivity analysis Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 10
- 238000001514 detection method Methods 0.000 abstract description 7
- 238000010276 construction Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 16
- 230000006870 function Effects 0.000 description 12
- 230000008569 process Effects 0.000 description 9
- 238000004590 computer program Methods 0.000 description 7
- 239000000203 mixture Substances 0.000 description 6
- 238000001615 p wave Methods 0.000 description 6
- 230000000153 supplemental effect Effects 0.000 description 5
- 241000208340 Araliaceae Species 0.000 description 4
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 4
- 235000003140 Panax quinquefolius Nutrition 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 235000008434 ginseng Nutrition 0.000 description 4
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- 230000000750 progressive effect Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000012216 screening Methods 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000011017 operating method Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000002023 wood Substances 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 241001274660 Modulus Species 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 229910021389 graphene Inorganic materials 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 238000004088 simulation Methods 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
-
- 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
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 embodiment of the application discloses a method and a device for determining reservoir parameters based on a three-dimensional elastic-electric rock physical quantity plate, wherein the method comprises the steps of obtaining three-dimensional elasticity and electric parameter actual measurement data of a target stratum; projecting the three-dimensional elasticity and electric parameter actual measurement data to a pre-constructed three-dimensional elastic-electric rock physical quantity plate to obtain a projection point, wherein the construction of the three-dimensional elastic-electric rock physical quantity plate comprises the following steps: performing gridding subdivision processing on the three-dimensional reservoir parameter space, and calculating elasticity and electrical parameter data corresponding to grid nodes to obtain a three-dimensional elasto-electric rock physical quantity version; calculating the distance between the projection point and each grid node, and determining the reservoir parameter data of the projection point according to the reservoir parameter data of the grid node corresponding to the minimum distance; and determining the reservoir parameter data of the target stratum according to the reservoir parameter data of the projection points. By utilizing the embodiments of the application, the precision of the reservoir parameters for oil and gas detection can be improved.
Description
Technical field
The present invention relates to petroleum gas seismic exploration technique fields, particularly, are related to a kind of based on three-dimensional bullet electricity rock
Physical quantity version determines the method and device of reservoir parameter.
Background technique
Rock physics template is by establishing petrophysical model appropriate, and simulation reservoir physical parameter and seismic response are special
The quantitative interpretation relationship between reservoir parameter and seismic properties is established in connection between sign.The establishment of this relationship is to reservoir object
Property and its state research be it is extremely important, can be enhanced the reliability that reservoir prediction is result, improve reservoir fluid and lithology is known
Other ability reduces the risk of quantitative interpretation.
On the basis of rock physics theoretical research, domestic and international some scholars to use rock physics template carry out reservoir ginseng
Several and lithology prediction has carried out exploration and research.Currently, based on two-dimentional rock physics template come predicting reservoir parameter or identification
Lithology comparative maturity, but two ginsengs of the two-dimentional rock physics template in prediction porosity, hydrocarbon saturation and shale content
It needs to fix one of parameter when number, can just obtain the prediction result of other two parameter in this way.Shale content is compared
Stable area is relatively applicable in, and changes the hydrocarbon saturation and shale content of bigger place prediction for shale content
It is larger with actual prediction resultant error.Therefore, need in the industry it is a kind of can more accurate predicting reservoir parameter method.
Summary of the invention
Being designed to provide for the embodiment of the present application a kind of determines reservoir parameter based on three-dimensional bullet electricity rock physics amount version
Method and device, can more accurate quantitative forecast reservoir parameter, further increase the physical parameter for oil and gas detection
Precision.
It is provided by the present application a kind of to determine that the method and device of reservoir parameter is logical based on three-dimensional bullet electricity rock physics amount version
It crosses including following manner realization:
A method of reservoir parameter is determined based on three-dimensional bullet electricity rock physics amount version, comprising:
Obtain the three dimensional elasticity and electrical parameter measured data of formation at target locations；
By the three dimensional elasticity and electrical parameter measured data to the three-dimensional bullet electricity rock physics amount version constructed in advance into
Row projection, obtains subpoint, wherein the building of three-dimensional bullet electricity rock physics amount version includes: to carry out three-dimensional reservoir parameter space
The processing of gridding subdivision calculates the corresponding elasticity of grid node and electrical parameter data, obtains three-dimensional bullet electricity rock physics amount
Version；
The subpoint is calculated at a distance from each grid node, is joined according to the reservoir apart from the corresponding grid node of minimum value
Number data determine the reservoir parameter data of the subpoint；
The reservoir parameter data of the formation at target locations are determined according to the reservoir parameter data of the subpoint.
The method for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version of the embodiment of the present application, it is described to store up three-dimensional
Layer parameter space carries out the processing of gridding subdivision, comprising:
The reservoir parameter data value interval in target work area is determined according to log data；
The reservoir parameter data value interval in three-dimensional reservoir parameter space is based on default accuracy value to divide,
Obtain gridding subdivision treated three-dimensional reservoir parameter space.
The method for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version of the embodiment of the present application, the calculating grid
The corresponding elastic parameter data of node, comprising:
Elastic rock computation model is constructed based on ellipsoid particle contact model, is calculated according to the elastic rock computation model
Determine the corresponding elastic parameter data of grid node.
The method for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version of the embodiment of the present application, the calculating grid
The corresponding electrical parameter data of node, comprising:
The corresponding electrical parameter data of grid node are calculated according to electrical rock computation model, wherein the electrical property rock
Computation model includes improved Xi Mendu formula:
Wherein, Rw、Rsh、RtRespectively indicate formation water resistivity, clay resistivity, formation resistivity, φ, Sw、VshRespectively
Indicate that porosity, water saturation, shale content, a, m and n are constant.
The method for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version of the embodiment of the present application, the calculating determine
The corresponding elastic parameter data of grid node, comprising:
The elastic parameter includes velocity of longitudinal wave；
The bulk modulus and modulus of shearing of dry rock are calculated based on ellipsoid particle contact model；
According to the rock volume modulus and shearing under the conditions of the bulk modulus and modulus of shearing, saturated with fluid of the dry rock
Modulus, rock density, which calculate, obtains velocity of longitudinal wave data.
The method for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version of the embodiment of the present application, it is described to be based on ellipsoid
Particle contact model calculates the bulk modulus and modulus of shearing of dry rock, comprising:
Wherein, Kdry、GdryRespectively indicate the bulk modulus and modulus of shearing of dry rock, F1=1- [α-0.0903 -1]1.456,
F2=1- [α-0.1026 -1]1.531,ξ is the correction constant for considering shearing stress direction, and α is vertical for hole
Horizontal ratio, n are ligancy, νma、GmaRespectively indicate matrix Poisson's ratio, modulus of shearing, effective pressure Peff=(ρg-ρw) gh, ρg、ρwPoint
Not Biao Shi density of matrix, fluid density, g indicate acceleration of gravity, h indicate depth of stratum.
The method for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version of the embodiment of the present application, the calculating grid
The corresponding elasticity of node and electrical parameter data, comprising:
According to log data or core data to the elastic rock computation model and electrical rock computation model into
Row correction；
Elastic rock computation model and electrical rock computation model are stated using after correction, it is corresponding to calculate grid node
Elasticity and electrical parameter data.
The method for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version of the embodiment of the present application, it is described according to distance
The reservoir parameter data of the corresponding grid node of minimum value determine the reservoir parameter data of the subpoint, comprising:
Filter out the grid node that the subpoint is less than preset threshold at a distance from each grid node；
The grid node that distance is less than preset threshold is ranked up according to the size of distance, is obtained corresponding apart from minimum value
Grid node.The method for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version of the embodiment of the present application, the acquisition
The three dimensional elasticity and electrical parameter measured data of formation at target locations, comprising:
Obtain the earthquake data before superposition and electromagnetic data of formation at target locations；
Elastic parameter measured data is obtained according to the prestack earthquake data inversion, and according to the electromagnetic data inverting
Obtain electrical parameter measured data.
The method for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version of the embodiment of the present application, the acquisition target
The three dimensional elasticity and electrical parameter measured data on stratum, comprising:
The elastic parameter and electrical parameter for obtaining preset kind are joined according to the elastic parameter to preset kind and electrically
The sensitivity analysis of the opposite reservoir parameter of number is as a result, it is preferred that go out sensibility elasticity parameter and sensitive electrical parameter；
Using the sensibility elasticity parameter and sensitive electrical parameter as the three dimensional elasticity of formation at target locations and electrical parameter,
Obtain the three dimensional elasticity and electrical parameter measured data of formation at target locations.
On the other hand, the embodiment of the present application, which also provides, a kind of determines reservoir parameter based on three-dimensional bullet electricity rock physics amount version
Device, comprising:
Module is obtained, for obtaining the three dimensional elasticity and electrical parameter measured data of formation at target locations；
Projection module, for by the three dimensional elasticity and electrical parameter measured data to the three-dimensional bullet electricity rock constructed in advance
Stone physical quantity version is projected, and obtains subpoint, wherein the building of three-dimensional bullet electricity rock physics amount version includes: by three-dimensional reservoir
Parameter space carries out the processing of gridding subdivision, calculates the corresponding elasticity of grid node and electrical parameter data, obtains three-dimensional bullet
Electric rock physics amount version；
Reservoir parameter determining module, for calculating the subpoint at a distance from each grid node, according to apart from minimum value
The reservoir parameter data of corresponding grid node determine the reservoir parameter data of the subpoint, according to the reservoir of the subpoint
Supplemental characteristic determines the reservoir parameter data of the formation at target locations.
The device that reservoir parameter is determined based on three-dimensional bullet electricity rock physics amount version of the embodiment of the present application, including processor and
For the memory of storage processor executable instruction, when described instruction is executed by the processor realize the following steps are included:
Obtain the three dimensional elasticity and electrical parameter measured data of formation at target locations；
By the three dimensional elasticity and electrical parameter measured data to the three-dimensional bullet electricity rock physics amount version constructed in advance into
Row projection, obtains subpoint, wherein the building of three-dimensional bullet electricity rock physics amount version includes: to carry out three-dimensional reservoir parameter space
The processing of gridding subdivision calculates the corresponding elasticity of grid node and electrical parameter data, obtains three-dimensional bullet electricity rock physics amount
Version；
The subpoint is calculated at a distance from each grid node, is joined according to the reservoir apart from the corresponding grid node of minimum value
Number data determine the reservoir parameter data of the subpoint, with determining the target according to the reservoir parameter data of the subpoint
The reservoir parameter data of layer.
One kind that this specification one or more embodiment provides determines that reservoir is joined based on three-dimensional bullet electricity rock physics amount version
Several method and devices can determine reservoir parameter by constructing three-dimensional bullet electricity rock physics amount version in advance in three dimensions
With the corresponding relationship between bullet-electrical parameter；Then the three dimensional elasticity of formation at target locations and electrical parameter measured data are projected to
In the three-dimensional bullet electricity rock physics amount version constructed in advance, to quantitatively determine the corresponding reservoir parameter data of formation at target locations.Utilize this
Apply for each embodiment, bullet-accurately can be portrayed between electrical parameter and reservoir parameter according to three-dimensional bullet electricity rock physics amount version
3-D quantitative explain relationship, be used in three-dimensional bullet-electrical parameter data inverting and obtain reservoir parameter data, to further mention
Precision of the height for the physical parameter of oil and gas detection.
Detailed description of the invention
In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
The some embodiments recorded in application, for those of ordinary skill in the art, in the premise of not making the creative labor property
Under, it is also possible to obtain other drawings based on these drawings.In the accompanying drawings:
Fig. 1 is that a kind of method for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version that this specification provides is implemented
The flow diagram of example；
Stream based on three-dimensional bullet electricity physics template quantitative forecast reservoir parameter in the example that Fig. 2 provides for this specification
Cheng Tu；
Three-dimensional based on the building of velocity of longitudinal wave, density and logarithmic resistance rate in the example that Fig. 3 provides for this specification
The schematic diagram of bullet electricity physics template prediction gas saturation；
Three-dimensional based on the building of velocity of longitudinal wave, density and logarithmic resistance rate in the example that Fig. 4 provides for this specification
The schematic diagram of bullet electricity physics template prediction porosity；
Three-dimensional based on the building of velocity of longitudinal wave, density and logarithmic resistance rate in the example that Fig. 5 provides for this specification
The schematic diagram of bullet electricity physics template prediction shale content；
P-wave impedance, density and the logarithmic resistance rate number of log measurement in another example that Fig. 6 provides for this specification
According to schematic diagram；
P-wave impedance, density and logarithmic resistance based on log measurement in another example that Fig. 7 provides for this specification
Gas saturation, porosity and the shale content schematic diagram data of rate data prediction；
Fig. 8 is a kind of device implementation that reservoir parameter is determined based on three-dimensional bullet electricity rock physics amount version that this specification provides
The modular structure schematic diagram of example.
Specific embodiment
In order to make those skilled in the art more fully understand the technical solution in this specification, below in conjunction with this explanation
Attached drawing in book one or more embodiment carries out the technical solution in this specification one or more embodiment clear, complete
Site preparation description, it is clear that described embodiment is only specification a part of the embodiment, instead of all the embodiments.Based on saying
Bright book one or more embodiment, it is obtained by those of ordinary skill in the art without making creative efforts all
The range of this specification example scheme protection all should belong in other embodiments.
When for current two-dimentional rock physics amount version come predicting reservoir parameter or identification lithology, two-dimensional measurement can only be utilized
The problem of data carry out quantitative forecast two dimension reservoir parameter, can not accurately complete the prediction work of multidimensional reservoir parameter.The application
A kind of three-dimensional bullet electricity rock physics amount version is provided quantitatively to portray oil-bearing reservoir.For example, can be by constructing three-dimensional bullet electricity
Rock physics amount version, which is established, a kind of can accurately portray elastic parameter, electrical parameter and reservoir parameter (such as porosity-saturation degree-mud
Matter content) between 3-D quantitative explain relationship, be used in three-dimensional bullet-electrical parameter data inverting and obtain reservoir parameter data.From
And it can be further improved the precision of the physical parameter for oil and gas detection using scheme provided by the present application.
In each implementation of the application, the type of elastic parameter, electrical parameter and reservoir parameter to be predicted can basis
Actual conditions determine.For example, the reservoir parameter may include reservoir porosity, gas saturation, shale content etc., the bullet
Property parameter may include longitudinal wave (or shear wave) speed, density, longitudinal wave (shear wave) impedance, Poisson's ratio, P-S wave velocity ratio etc., institute
Stating electrical parameter may include resistivity, logarithmic resistance rate etc..
Fig. 1 is a kind of method for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version that this specification provides
Embodiment flow diagram.Although present description provides as the following examples or method operating procedure shown in the drawings or device
Structure, but may include more in the method or device or part merging based on routine or without creative labor
Less operating procedure or modular unit afterwards.In the step of there is no necessary causalities in logicality or structure, these steps
The modular structure of rapid execution sequence or device is not limited to this specification embodiment or execution shown in the drawings sequence or module knot
Structure.Device in practice, server or the end product of the method or modular structure are in application, can be according to embodiment
Either method or modular structure carry out sequence execution shown in the drawings or parallel execute (such as parallel processor or multithreading
The environment of processing, the even implementation environment including distributed treatment, server cluster).
Specific one embodiment is as shown in Figure 1, this specification offer determines reservoir based on three-dimensional bullet electricity rock amount version
In one embodiment of the method for parameter, the method may include:
S2, the three dimensional elasticity and electrical parameter measured data for obtaining formation at target locations.
In the present embodiment, it can be directed to the different target stratum in target work area, obtain corresponding three dimensional elasticity and electrical property
Parameter measured data.Wherein, elastic parameter may include longitudinal wave (or shear wave) speed, rock density, longitudinal wave (shear wave) impedance,
Poisson's ratio, P-S wave velocity ratio etc., the electrical parameter may include resistivity or the logarithmic resistance rate on stratum etc..It is specific real
Shi Shi, can according in the three-dimensional bullet electricity rock physics amount version constructed in advance elastic parameter and electrical parameter type determine need
The three dimensional elasticity and electrical parameter measured data to be obtained.
For example, elastic parameter and electrical parameter type difference in the three-dimensional bullet electricity rock physics amount version constructed in advance
Are as follows: velocity of longitudinal wave, rock density and logarithmic resistance rate, correspondingly, the velocity of longitudinal wave of available formation at target locations, rock density
And logarithmic resistance rate measured data, form the three dimensional elasticity and electrical parameter measured data of formation at target locations.This specification
In one embodiment, then the earthquake data before superposition and electromagnetic data in available target work area utilize earthquake data before superposition
Inverting obtains elastic parameter data, and obtains electrical parameter data using electromagnetic data inverting, to improve the elasticity of acquisition
And the accuracy of electrical parameter data.
S4, by the three dimensional elasticity and electrical parameter measured data to the three-dimensional bullet electricity rock physics amount version constructed in advance
It is projected, obtains subpoint.
It, can be by three dimensions of the elasticity obtained in above-mentioned steps and electrical parameter measured data composition in the present embodiment
According to being projected to the three-dimensional bullet electricity rock physics amount version constructed in advance, obtain the corresponding subpoint of the three-dimensional data.For example,
Assuming that above-mentioned three-dimensional data is (vp,ρ,Rt), wherein vpIndicate that velocity of longitudinal wave, ρ indicate rock density, RtIndicate resistivity.Then
It can be by (vp,ρ,Rt) projected to three-dimensional bullet electricity rock physics amount version, obtain (vp,ρ,Rt) three-dimensional bullet electricity rock object
Subpoint position in reason amount version space.
In one embodiment of this specification, the building of the three-dimensional bullet electricity rock physics amount version may include: will be three-dimensional
Reservoir parameter space carries out the processing of gridding subdivision, then calculates the corresponding elasticity of each grid node and electrical parameter data,
To which the corresponding reservoir parameter data of grid node each in three-dimensional reservoir parameter space and three-dimensional to be played to the bullet-in electrical parameter space
Electrical parameter data associates one by one, obtains three-dimensional bullet electricity rock physics amount version.
Further, in one embodiment of this specification, it can determine that the reservoir parameter data in target work area become first
Change range, each reservoir parameter data value interval is determined according to the reservoir parameter data variation range；Then, by three-dimensional reservoir
The reservoir parameter data value interval in parameter space is based on default accuracy value and is divided, thus by three-dimensional reservoir parameter
Space carries out gridding subdivision, obtains gridding subdivision treated three-dimensional reservoir parameter space.When it is implemented, described default
Accuracy value can be preset according to the actual situation in the case where taking into account precision with efficiency.The reservoir in the target work area
Supplemental characteristic variation range can be determined according to measured datas such as such as log datas.Then, computational gridding subdivision is handled
The corresponding elastic and electrical parameter data of each grid node in three-dimensional reservoir parameter space afterwards, to construct three-dimensional bullet electricity rock
Stone physical quantity version.
In the present embodiment, by the data variation model for determining each reservoir parameter according to the measured data in target work area first
It encloses, then, the processing of gridding subdivision is carried out to three-dimensional reservoir parameter space in the range, the effect of gridding processing can be improved
Rate；At the same time it can also reduce grid node number, the data calculation amount of the corresponding bullet-electrical parameter of grid node is reduced, is significantly mentioned
The efficiency of high data processing.Further, data screening range can also be effectively limited using measured data, reduced screening and missed
Difference improves the accuracy of the subsequent reservoir parameter data that target work area is determined based on bullet electricity rock physics amount version.
For example, in one or more embodiment of this specification, it is assumed that the reservoir parameter type include: porosity,
Water saturation and shale content.The hole in target work area can be then determined according to the data that rock core and well log interpretation obtain first
The variation range of three porosity, water saturation and shale content reservoir parameters, and according to the reservoir parameter data variation
Range determines each reservoir parameter data value interval.Then, the reservoir parameter data in three-dimensional reservoir parameter space are taken
Value section is based on default accuracy value accordingly and is divided.Correspondingly, it is hole that porosity, which corresponds to the equal point of numerical value to be linked to be line,
The corresponding equipotential lines of porosity, similarly, it is aqueous full that water saturation and shale content, which correspond to the equal point of numerical value to be linked to be line,
Equipotential lines corresponding with degree and shale content.The intersection point of each parameter equipotential lines is grid node, to be formed by porosity, be contained
The three-dimensional grid space of water saturation and shale content composition, to obtain gridding subdivision treated three-dimensional reservoir parameter
Space.
It is then possible to the corresponding elasticity of each grid node and electrically in three-dimensional reservoir parameter space after computational gridding
Supplemental characteristic.It, can be according to the elastic rock computation model and electrical property constructed in advance in some embodiments of this specification
Rock computation model calculates the corresponding elasticity of each grid node and electrical parameter data.In one embodiment of this specification,
Elastic rock computation model can be constructed based on ellipsoid particle contact model, be calculated and determined according to the elastic rock computation model
Elastic parameter data.Wherein, the ellipsoid particle contact model includes: that formation rock can be equivalent to ellipsoid particle, i.e.,
Formation pore aspect ratio is not equal to 1.Formation rock is equivalent to ellipsoid particle and more meets stratum actual conditions, thus based on ellipse
Ball particle contact model can be further improved the accuracy of elastic parameter data calculating.
In one or more embodiment of this specification, it is assumed that the corresponding elasticity of three-dimensional bullet electricity rock physics amount version and
Electrical parameter is respectively vp,ρ,Rt, reservoir parameter to be predicted is porosity φ, water saturation SwWith shale content Vsh.Then exist
In the above-mentioned three-dimensional grid space being made of porosity, water saturation and shale content, for each grid node φ,
Sw、VshThree-dimensional data, Voigt-Roy Si-Ilyushin (Voigt-Reuss- can be utilized by rock mineral composition and its content
Hill) the bulk modulus and modulus of shearing of formula average computation Rock Matrix.It is then possible to based on ellipsoid particle contact model
Calculate the bulk modulus K of dry rockdryAnd shear modulus Gdry, wherein corresponding calculation formula can be expressed as:
In formula, F1=1- [α-0.0903 -1]1.456, F2=1- [α-0.1026 -1]1.531,ξ is to examine
Consider the correction constant in shearing stress direction, α is pore components, and n is ligancy, ν ma, GmaRespectively indicate matrix Poisson's ratio, shearing
Modulus, effective pressure Peff=(ρg-ρw) gh, ρg、ρwDensity of matrix, fluid density are respectively indicated, g indicates acceleration of gravity, h table
Show depth of stratum.
It is then possible to calculate the equivalent volume modulus of fluid-mixing using Wood (Wood) formula:
In formula, KflFor the bulk modulus of fluid, KgFor the bulk modulus of oil gas, KwFor the bulk modulus of water.
Later, it using the rock volume modulus in the case of Jia Siman (Gassmann) formula calculating saturated with fluid and can cut
Shear modulu, formula are as follows:
Gsat=Gdry (5)
In formula, KsatAnd GsatThe respectively bulk modulus and modulus of shearing of saturated with fluid.
Later, the density that rock can be calculated by rock physics volume-based model, then calculates velocity of longitudinal wave again:
In formula, ρ is the density of saturated rock.
In another embodiment of this specification, resistivity can be calculated using improved Xi Mendu formula::
In formula, Rw、RshAnd RtRespectively formation water resistivity, clay resistivity and formation resistivity, a, m and n be Ah
The constant of your odd formula.
It repeats the above process, the elastic parameter data and electrical parameter data for obtaining all grid nodes can be calculated, it will
The corresponding three-dimensional space grid of reservoir parameter three-dimensional space corresponding with bullet-electrical parameter associates one by one, forms three-dimensional bullet electricity
Rock physics amount version.
It, can also be according to log data or core data to the elastic rock meter in one embodiment of this specification
It calculates model and electrical rock computation model is corrected, when it is implemented, can be to elastic rock computation model and electrical property
Constant data in rock computation model is corrected.Then, the elastic rock computation model and electrical rock after correction are utilized
Stone computation model calculates the corresponding elasticity of each grid node and electrical parameter data, to further increase three-dimensional bullet electricity rock
The accuracy of stone physical quantity version building.
S6, the subpoint is calculated at a distance from each grid node, according to the storage apart from the corresponding grid node of minimum value
Layer parameter data determine the reservoir parameter data of the subpoint, determine the mesh according to the reservoir parameter data of the subpoint
Mark the reservoir parameter data on stratum.
It, can be in predetermined reservoir parameter data variation range for the subpoint obtained in S4 in the present embodiment
It is interior, calculate the distance between each grid node (such as Euclidean distance) in subpoint and three-dimensional bullet electricity rock physics amount version.Then,
Subpoint corresponding grid node of minimum value at a distance from each grid node is filtered out, is obtained described apart from the corresponding net of minimum value
Reservoir parameter number of the reservoir parameter data (such as: porosity, water saturation and shale content) of lattice node as the subpoint
According to.And subpoint and the actual measurement bullet of each formation at target locations-electrical parameter data correspond, it is determined that the reservoir parameter data of subpoint
Namely the reservoir parameter data of each formation at target locations are determined.To predicted in target work area respectively using above-mentioned steps are quantitative
The reservoir parameter data on stratum.
In one embodiment of this specification, subpoint can be filtered out first at a distance from each grid node less than default
Then the grid node that distance is less than preset threshold is ranked up, thus really by the grid node of threshold value according to the size of distance
Make grid node corresponding apart from minimum value.And then obtain the reservoir ginseng apart from the corresponding grid node of minimum value
Reservoir parameter data of the number data as the subpoint.Wherein, the size of the preset threshold can carry out according to the actual situation
Adjustment, for example, the size of preset threshold can be reduced if the grid node filtered out is more；, whereas if not screening
Grid node out, then can the size appropriate that expand preset threshold.Using the present embodiment above scheme, can be further reduced
Operand improves data-handling efficiency.
It repeats the above steps, can quantitatively set the goal the corresponding reservoir parameter data in each stratum in work area really.Meanwhile this theory
The bright each embodiment of book, which also passes through, to be introduced electrical parameter (such as resistivity) information to compensate for elastic parameter insensitive to saturation degree
Problem, to further increase the precision of prediction of hydrocarbon saturation.
It, can also be before constructing three-dimensional bullet electricity rock physics amount version, according to not in some embodiments of this specification
With target work area, preferably corresponding sensibility elasticity parameter and sensitive electrical parameter out, or be adapted in the target work area
Reservoir parameter to be predicted.Then, it is more suitable for using sensitive bullet (electricity) the property parameter and with the building of predicting reservoir parameter
The three-dimensional bullet electricity rock physics amount version in the target work area, so that the three-dimensional bullet electricity rock physics amount version of building more adapts to difference
The complicated landform situation in area.
In one embodiment of this specification, the elastic parameter and electrical parameter of available preset kind are then right
The elastic parameter and electrical parameter of preset kind carry out the sensitivity analysis of reservoir parameter relatively to be predicted, based on the analysis results
It is preferred that going out sensibility elasticity or electrical parameter.The sensibility elasticity or electrical parameter that preferably go out are used to construct three-dimensional bullet electricity rock
Physical quantity version, and the sensibility elasticity of formation at target locations or the measured data of electrical parameter are obtained, for determining target
The reservoir parameter data of layer.For example, when can be changed with bullet (electricity) the property supplemental characteristic of the multiple types of comparative analysis, accordingly to pre-
Survey the variation degree of reservoir parameter.Using reservoir parameter change to be predicted it is relatively obvious when corresponding bullet (electricity) property parameter type as
Sensitive bullet (electricity) property parameter.Using the scheme of the present embodiment, the bullet electricity rock physics amount version of foundation can be improved to formation at target locations
Adaptability, to further increase the accuracy of reservoir parameter forecast.
Certainly, can also be according to the actual conditions of objective area in some embodiments of this specification, selection can be compared with
The elastic parameter or electrical parameter easily obtained constructs three-dimensional bullet electricity rock physics amount version.Then, for three-dimensional bullet electricity rock object
Elastic parameter and electrical parameter in reason amount version obtain the three dimensional elasticity and electrical parameter measured data of formation at target locations, into
One step determines the reservoir parameter of the formation at target locations.
So that the scheme in above-described embodiment, can be applied not only to conventional siltstone, selection can also be passed through
Suitable elasticity, electrical parameter and reservoir parameter construct different three-dimensional bullet electricity rock physics amount versions, to meet complicated reservoirs
The application environment of parameter prediction.Also, the scheme that each embodiment of this specification provides can also be applied to multiple scale (rocks
The heart, well logging and earthquake-electromagnetism) data predicting reservoir parameter, to provide a kind of very simple reality for quantitative forecast reservoir
Method.
In order to further verify the practicability and feasibility of this method, this specification additionally provides one using above scheme
A specific example, Fig. 2 indicates the flow chart using three-dimensional bullet electricity rock physics amount version quantitative forecast reservoir parameter, according to target work
The actual conditions in area construct elastic parameter and electrical parameter that three-dimensional bullet electricity rock physics amount version uses in this example are as follows: vertical
Wave velocity, rock density and logarithmic resistance rate, reservoir parameter to be predicted are as follows: gas saturation, porosity and shale contain
Amount.As shown in Fig. 2, specific embodiment can be such that
(1) elastic parameter and electrical parameter measured data are obtained.
Earthquake data before superposition, the electromagnetic data in target work area are obtained, and inverting obtains elastic parameter (velocity of longitudinal wave vp, rock
Stone density p) and electrical parameter (logarithmic resistance rate log10 (Rt)) measured data.
(2) measured data is projected to the three-dimensional bullet electricity rock physics amount version constructed in advance.
Fig. 3, Fig. 4, Fig. 5 are respectively indicated using the three-dimensional bullet electricity rock physics amount version prediction porosity φ constructed in advance, are contained
Water saturation SwAnd shale content VshSchematic diagram, by velocity of longitudinal wave (P-wave velocity), rock density
(Density) and logarithmic resistance rate (log10 (Rt)) composition three-dimensional measured data into three-dimensional bullet electricity rock physics amount version
Projection, as shown in Figure 3-Figure 5, the point in figure indicate corresponding subpoint.Wherein, the unit of velocity of longitudinal wave are as follows: thousand meter per second (Km/
S), the unit of rock density are as follows: gram/cc (g/cm3), the unit of logarithmic resistance rate are as follows: ohm meter (ohm-m).
Correspondingly, three-dimensional bullet electricity rock physics amount version can be constructed in advance according to following steps:
It is possible, firstly, to determine the change of porosity, water saturation and shale content according to rock core and well log interpretation achievement
Change range, porosity, water saturation and the shale content in corresponding data variation range are carried out according to preset requirement
Equal part.To which three-dimensional reservoir parameter space is carried out gridding subdivision, formed by porosity, water saturation and shale content
The three-dimensional space grid of composition.
As shown in Figure 3-Figure 5, the variation range of porosity (Porosity) is 0.1-0.3, is carried out etc. with 0.1 for step-length
Point；The variation range of gas saturation (gas saturation) is 0.1-0.6, carries out equal part with 0.1 for step-length；Shale content
The variation range of (Clay content) is 0-0.8, carries out equal part with 0.4 for step-length.Wherein, porosity 0.1,0.2,0.3
Line be the corresponding equipotential lines of porosity, the line of gas saturation 0.1,0.2,0.3,0.4,0.5,0.6 is that gassiness is full
The corresponding equipotential lines with degree, the line of shale content 0,0.4,0.8 is the corresponding equipotential lines of shale content.Porosity, gassiness
The intersection point that saturation degree and each equipotential lines of shale content are formed is grid node.To be formed by porosity, aqueous full
With the three-dimensional space grid of degree and shale content composition.
Secondly, porosity corresponding for each grid node, water saturation and shale contain in three-dimensional space grid
Amount is based on ellipsoid particle contact model, can use formula (1)-(6) and calculates determining velocity of longitudinal wave, and according to rock physics
Volume-based model, which calculates, determines rock density；And calculated using improved Xi Mendu formula (7) and determine resistivity, to obtain logarithm
Resistivity data.To reservoir parameter data (φ, S corresponding for each grid nodew,Vsh), it can use preparatory building
The calculating of elastic and electrical rock computation model determine corresponding three-dimensional bullet-electrical parameter data (vp,ρ,log10(Rt)).From
And the corresponding three-dimensional space grid of reservoir parameter three-dimensional space corresponding with bullet-electrical parameter associates one by one, it obtains three-dimensional
Bullet electricity rock physics amount version.
(3) reservoir parameter of formation at target locations is determined.
To the subpoint of acquisition, subpoint and three can be calculated within the scope of predetermined reservoir parameter data variation
Tie up the distance between each grid node in bullet electricity rock physics amount version.And the subpoint is calculated at a distance from each grid node,
The reservoir parameter data of the subpoint are determined according to the reservoir parameter data apart from the corresponding grid node of minimum value.Subpoint
It is corresponded with the actual measurement bullet of each formation at target locations-electrical parameter data, it is determined that the reservoir parameter data of subpoint namely determined
The reservoir parameter data of each formation at target locations.To the reservoir that has predicted target work area Zhong Ge stratum quantitative using above-mentioned steps
Supplemental characteristic.
Fig. 6 and Fig. 7, which is respectively indicated, does not consolidate reservoir log data based on three-dimensional bullet electricity rock physics amount version quantitative forecast
Example schematic.Wherein, elastic parameter are as follows: p-wave impedance (P-wave Impedance), rock density (Density), electricity
Property parameter are as follows: logarithmic resistance rate (log10 (Rt)), wherein the unit of p-wave impedance are as follows: thousand * grams/cc of meter per second (km/
s*g/cm3).Reservoir parameter to be predicted are as follows: gas saturation (gas saturation), porosity (Porosity) and mud
Matter content (Clay content).
Fig. 6 (a), Fig. 6 (b), Fig. 6 (c) respectively indicate the p-wave impedance that actual measurement and inverting obtain, rock density and right
Number resistivity is with depth of stratum (Depth), and wherein solid line indicates p-wave impedance, rock density and the logarithmic resistance rate number of actual measurement
According to dotted line indicates p-wave impedance, the rock obtained according to above-mentioned elastic rock computation model and electrical rock computation model inverting
Stone density and logarithmic resistance rate data.Fig. 7 (a), Fig. 7 (b), Fig. 7 (c) respectively indicate the gas saturation of prediction, porosity
And the data of shale content.Fig. 7 (a), Fig. 7 (b), the dot in Fig. 7 (c) indicate to log well actually measured gas saturation,
Porosity and shale content data, solid line indicate the gas saturation using three-dimensional bullet electricity rock physics amount version prediction, hole
Degree and shale content data.Wherein, the unit of depth of stratum are as follows: the downward depth in seabed, unit are rice (mbsf).
Analysis chart 6 (a), Fig. 6 (b), Fig. 6 (c) be not it is found that only the dotted line in Fig. 6 (b) and solid line are completely coincident.It is i.e. sharp
The p-wave impedance and logarithmic resistance rate data and measured data obtained with the computation model inverting after the correction of this specification embodiment
Substantially it coincide, only there are some deviations for rock density actual measurement and inverting data.Thus elasticity or electricity after showing based on correction
Property rock computation model inverting obtain electricity-elastic data accuracy it is higher, it is pre- to further improve subsequent reservoir parameter data
The accuracy of survey.
Analysis chart 7 (a), Fig. 7 (b), Fig. 7 (c) it is found that based on Fig. 6 (a), Fig. 6 (b), Fig. 6 (c) inverting bullet-electrical parameter
Porosity, gas saturation and the shale content and result of log interpretation that data predict coincide preferable.To show to utilize this
The scheme that specification said one or multiple embodiments provide, can greatly improve reservoir parameter (especially water saturation)
Precision of prediction, and can preferably complete the prediction work of multi-source data and multidimensional reservoir parameter.
All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment
Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Specifically it is referred to
The description of aforementioned relevant treatment related embodiment, does not do repeat one by one herein.
It is above-mentioned that this specification specific embodiment is described.Other embodiments are in the scope of the appended claims
It is interior.In some cases, the movement recorded in detail in the claims or step can be come according to the sequence being different from embodiment
It executes and desired result still may be implemented.In addition, process depicted in the drawing not necessarily require show it is specific suitable
Sequence or consecutive order are just able to achieve desired result.In some embodiments, multitasking and parallel processing be also can
With or may be advantageous.
One kind that this specification one or more embodiment provides determines that reservoir is joined based on three-dimensional bullet electricity rock physics amount version
Several methods can determine reservoir parameter and bullet-electricity by constructing three-dimensional bullet electricity rock physics amount version in advance in three dimensions
Corresponding relationship between parameter；Then the three dimensional elasticity of formation at target locations and electrical parameter measured data are projected into preparatory building
Three-dimensional bullet electricity rock physics amount version in, to quantitatively determine the corresponding reservoir parameter data of formation at target locations.It is each using the application
It is fixed can accurately to portray three-dimensional of the bullet-between electrical parameter and reservoir parameter according to three-dimensional bullet electricity rock physics amount version for embodiment
Explanation relationship is measured, inverting in three-dimensional bullet-electrical parameter data is used for and obtains reservoir parameter data, to further increase for oil
The precision of the physical parameter of gas detection.
Based on the method described above for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version, this specification one
Or multiple embodiments also provide a kind of device that reservoir parameter is determined based on three-dimensional bullet electricity rock physics amount version.The device can
With system, software (application), module, component, server for include the use of this specification embodiment the method etc. and combine
The necessary device for implementing hardware.Based on same innovation thinking, in one or more embodiments that this specification embodiment provides
Device as described in the following examples.Since the implementation that device solves the problems, such as is similar to method, this specification is real
The implementation for applying the specific device of example may refer to the implementation of preceding method, and overlaps will not be repeated.It is used below, term
The combination of the software and/or hardware of predetermined function may be implemented in " unit " or " module ".Although described in following embodiment
Device is preferably realized with software, but the realization of the combination of hardware or software and hardware is also that may and be contemplated.
Specifically, a kind of Installation practice that reservoir parameter is determined based on three-dimensional bullet electricity rock physics amount version that Fig. 8 specification provides
Modular structure schematic diagram, such as Fig. 8, the apparatus may include:
Module 102 is obtained, can be used for obtaining the three dimensional elasticity and electrical parameter measured data of formation at target locations；
Projection module 104 can be used for the three dimensional elasticity and electrical parameter measured data to three constructed in advance
Dimension bullet electricity rock physics amount version projected, obtain subpoint, wherein the building of three-dimensional bullet electricity rock physics amount version include: by
Three-dimensional reservoir parameter space carries out the processing of gridding subdivision, calculates the corresponding elasticity of grid node and electrical parameter data, obtains
Obtain three-dimensional bullet electricity rock physics amount version；
Reservoir parameter determining module 106 can be used for calculating the subpoint at a distance from each grid node, according to distance
The reservoir parameter data of the corresponding grid node of minimum value determine the reservoir parameter data of the subpoint；According to the subpoint
Reservoir parameter data determine the reservoir parameter data of the formation at target locations.
It should be noted that device described above can also include other embodiment party according to the description of embodiment of the method
Formula.Concrete implementation mode is referred to the description of related method embodiment, does not repeat one by one herein.
One kind that this specification one or more embodiment provides determines that reservoir is joined based on three-dimensional bullet electricity rock physics amount version
Several devices can determine reservoir parameter and bullet-electricity by constructing three-dimensional bullet electricity rock physics amount version in advance in three dimensions
Corresponding relationship between parameter；Then the three dimensional elasticity of formation at target locations and electrical parameter measured data are projected into preparatory building
Three-dimensional bullet electricity rock physics amount version in, to quantitatively determine the corresponding reservoir parameter data of formation at target locations.It is each using the application
It is fixed can accurately to portray three-dimensional of the bullet-between electrical parameter and reservoir parameter according to three-dimensional bullet electricity rock physics amount version for embodiment
Explanation relationship is measured, inverting in three-dimensional bullet-electrical parameter data is used for and obtains reservoir parameter data, to further increase for oil
The precision of the physical parameter of gas detection.
Method or apparatus described in above-described embodiment that this specification provides can realize that business is patrolled by computer program
It collects and records on a storage medium, the storage medium can be read and be executed with computer, realize this specification embodiment institute
The effect of description scheme.Therefore, this specification, which also provides, a kind of determines reservoir parameter based on three-dimensional bullet electricity rock physics amount version
Device, the memory including processor and storage processor executable instruction, described instruction are realized when being executed by the processor
The following steps are included:
Obtain the three dimensional elasticity and electrical parameter measured data of formation at target locations；
By the three dimensional elasticity and electrical parameter measured data to the three-dimensional bullet electricity rock physics amount version constructed in advance into
Row projection, obtains subpoint, wherein the building of three-dimensional bullet electricity rock physics amount version includes: to carry out three-dimensional reservoir parameter space
The processing of gridding subdivision calculates the corresponding elasticity of grid node and electrical parameter data, obtains three-dimensional bullet electricity rock physics amount
Version；
The subpoint is calculated at a distance from each grid node, is joined according to the reservoir apart from the corresponding grid node of minimum value
Number data determine the reservoir parameter data of the subpoint；With determining the target according to the reservoir parameter data of the subpoint
The reservoir parameter data of layer.
The storage medium may include the physical unit for storing information, usually by after information digitalization again with benefit
The media of the modes such as electricity consumption, magnetic or optics are stored.It may include: that letter is stored in the way of electric energy that the storage medium, which has,
The device of breath such as, various memory, such as RAM, ROM；The device of information is stored in the way of magnetic energy such as, hard disk, floppy disk, magnetic
Band, core memory, magnetic bubble memory, USB flash disk；Using optical mode storage information device such as, CD or DVD.Certainly, there are also it
Readable storage medium storing program for executing of his mode, such as quantum memory, graphene memory etc..
It should be noted that device described above can also include other embodiment party according to the description of embodiment of the method
Formula.Concrete implementation mode is referred to the description of related method embodiment, does not repeat one by one herein.
A kind of device determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version, Ke Yitong described in above-described embodiment
It crosses and constructs three-dimensional bullet electricity rock physics amount version in advance, determine in three dimensions corresponding between reservoir parameter and bullet-electrical parameter
Relationship；Then the three dimensional elasticity of formation at target locations and electrical parameter measured data are projected to the three-dimensional bullet electricity rock constructed in advance
In physical quantity version, to quantitatively determine the corresponding reservoir parameter data of formation at target locations.It, can basis using each embodiment of the application
Three-dimensional bullet electricity rock physics amount version accurately portrays 3-D quantitative of the bullet-between electrical parameter and reservoir parameter and explains relationship, uses it
Inverting obtains reservoir parameter data in three-dimensional bullet-electrical parameter data, to further increase the physical property ginseng for oil and gas detection
Several precision.
It should be noted that this specification device described above can also include according to the description of related method embodiment
Other embodiments, concrete implementation mode are referred to the description of embodiment of the method, do not repeat one by one herein.This explanation
Various embodiments are described in a progressive manner in book, and same and similar part refers to each other i.e. between each embodiment
Can, each embodiment focuses on the differences from other embodiments.It is situated between especially for hardware+program class, storage
For matter+program embodiment, since it is substantially similar to the method embodiment, so be described relatively simple, related place referring to
The part of embodiment of the method illustrates.
It is above-mentioned that this specification specific embodiment is described.Other embodiments are in the scope of the appended claims
It is interior.In some cases, the movement recorded in detail in the claims or step can be come according to the sequence being different from embodiment
It executes and desired result still may be implemented.In addition, process depicted in the drawing not necessarily require show it is specific suitable
Sequence or consecutive order are just able to achieve desired result.In some embodiments, multitasking and parallel processing be also can
With or may be advantageous.
System, device, module or the unit that above-described embodiment illustrates can specifically realize by computer chip or entity,
Or it is realized by the product with certain function.It is a kind of typically to realize that equipment is computer.Specifically, computer for example may be used
That thinks in personal computer, laptop computer, vehicle-mounted human-computer interaction device, tablet computer or these equipment any sets
Standby combination.
For convenience of description, it is divided into various modules when description apparatus above with function to describe respectively.Certainly, implementing this
The function of each module can be realized in the same or multiple software and or hardware when specification one or more, it can also be with
The module for realizing same function is realized by the combination of multiple submodule or subelement etc..Installation practice described above is only
It is only illustrative, for example, in addition the division of the unit, only a kind of logical function partition can have in actual implementation
Division mode, such as multiple units or components can be combined or can be integrated into another system or some features can be with
Ignore, or does not execute.Another point, shown or discussed mutual coupling, direct-coupling or communication connection can be logical
Some interfaces are crossed, the indirect coupling or communication connection of device or unit can be electrical property, mechanical or other forms.
It is also known in the art that other than realizing controller in a manner of pure computer readable program code, it is complete
Entirely can by by method and step carry out programming in logic come so that controller with logic gate, switch, specific integrated circuit, programmable
Logic controller realizes identical function with the form for being embedded in microcontroller etc..Therefore this controller is considered one kind
Hardware component, and the structure that the device for realizing various functions that its inside includes can also be considered as in hardware component.Or
Person even, can will be considered as realizing the device of various functions either the software module of implementation method can be hardware again
Structure in component.
The present invention be referring to according to the method for the embodiment of the present invention with the flowchart and/or the block diagram of computer program product
Come what is described.It should be understood that can be realized by computer program instructions each flow and/or block in flowchart and/or the block diagram,
And the combination of the process and/or box in flowchart and/or the block diagram.These computer program instructions be can provide to general meter
Calculation machine, special purpose computer, Embedded Processor or other programmable data processing devices processor to generate a machine, make
It obtains and is generated by the instruction that computer or the processor of other programmable data processing devices execute for realizing in flow chart one
The device for the function of being specified in a process or multiple processes and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
In a typical configuration, calculating equipment includes one or more processors (CPU), input/output interface, net
Network interface and memory.
It should also be noted that, the terms "include", "comprise" or its any other variant are intended to nonexcludability
It include so that the process, method, commodity or the equipment that include a series of elements not only include those elements, but also to wrap
Include other elements that are not explicitly listed, or further include for this process, method, commodity or equipment intrinsic want
Element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that including described want
There is also other identical elements in the process, method or equipment of element.
It will be understood by those skilled in the art that this specification one or more embodiment can provide as method, system or calculating
Machine program product.Therefore, this specification one or more embodiment can be used complete hardware embodiment, complete software embodiment or
The form of embodiment combining software and hardware aspects.Moreover, this specification one or more embodiment can be used at one or
It is multiple wherein include computer usable program code computer-usable storage medium (including but not limited to magnetic disk storage,
CD-ROM, optical memory etc.) on the form of computer program product implemented.
This specification one or more embodiment can computer executable instructions it is general on
It hereinafter describes, such as program module.Generally, program module includes executing particular task or realization particular abstract data type
Routine, programs, objects, component, data structure etc..This this specification one can also be practiced in a distributed computing environment
Or multiple embodiments, in these distributed computing environments, by being held by the connected remote processing devices of communication network
Row task.In a distributed computing environment, program module can be located at the local and remote computer including storage equipment
In storage medium.
All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment
Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for system reality
For applying example, since it is substantially similar to the method embodiment, so being described relatively simple, related place is referring to embodiment of the method
Part explanation.In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ",
The description of " specific example " or " some examples " etc. means specific features described in conjunction with this embodiment or example, structure, material
Or feature is contained at least one embodiment or example of this specification.In the present specification, to the signal of above-mentioned term
Property statement must not necessarily be directed to identical embodiment or example.Moreover, specific features, structure, material or the spy of description
Point may be combined in any suitable manner in any one or more of the embodiments or examples.In addition, without conflicting with each other,
Those skilled in the art can be by different embodiments or examples described in this specification and different embodiments or examples
Feature is combined.
The above description is only an example of the present application, is not intended to limit this application.For those skilled in the art
For, various changes and changes are possible in this application.All any modifications made within the spirit and principles of the present application are equal
Replacement, improvement etc., should be included within the scope of the claims of this application.
Claims (12)
1. a kind of method for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version characterized by comprising
Obtain the three dimensional elasticity and electrical parameter measured data of formation at target locations；
The three dimensional elasticity and electrical parameter measured data are thrown to the three-dimensional bullet electricity rock physics amount version constructed in advance
Shadow obtains subpoint, wherein the building of three-dimensional bullet electricity rock physics amount version includes: that three-dimensional reservoir parameter space is carried out grid
Change subdivision processing, calculates the corresponding elasticity of grid node and electrical parameter data, obtain three-dimensional bullet electricity rock physics amount version；
The subpoint is calculated at a distance from each grid node, according to the reservoir parameter number apart from the corresponding grid node of minimum value
According to the reservoir parameter data of the determination subpoint；
The reservoir parameter data of the formation at target locations are determined according to the reservoir parameter data of the subpoint.
2. the method according to claim 1 for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version, feature exist
In described that three-dimensional reservoir parameter space is carried out the processing of gridding subdivision, comprising:
The reservoir parameter data value interval in target work area is determined according to log data；
The reservoir parameter data value interval in three-dimensional reservoir parameter space is based on default accuracy value to divide, is obtained
Gridding subdivision treated three-dimensional reservoir parameter space.
3. the method according to claim 1 for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version, feature exist
In the corresponding elastic parameter data of the calculating grid node, comprising:
Elastic rock computation model is constructed based on ellipsoid particle contact model, is calculated and is determined according to the elastic rock computation model
The corresponding elastic parameter data of grid node.
4. the method according to claim 3 for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version, feature exist
In the corresponding electrical parameter data of the calculating grid node, comprising:
The corresponding electrical parameter data of grid node are calculated according to electrical rock computation model, wherein the electrical property rock calculates
Model includes improved Xi Mendu formula:
Wherein, Rw、Rsh、RtRespectively indicate formation water resistivity, clay resistivity, formation resistivity, φ, Sw、VshRespectively indicate hole
Porosity, water saturation, shale content, a, m and n are constant.
5. the method according to claim 1 for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version, feature exist
In the calculating determines the corresponding elastic parameter data of grid node, comprising:
The elastic parameter includes velocity of longitudinal wave；
The bulk modulus and modulus of shearing of dry rock are calculated based on ellipsoid particle contact model；
According to the rock volume modulus and shearing mould under the conditions of the bulk modulus and modulus of shearing, saturated with fluid of the dry rock
Amount, rock density, which calculate, obtains velocity of longitudinal wave data.
6. the method according to claim 5 for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version, feature exist
In the bulk modulus and modulus of shearing for calculating dry rock based on ellipsoid particle contact model, comprising:
Wherein, Kdry、GdryRespectively indicate the bulk modulus and modulus of shearing of dry rock, F1=1- [α-0.0903-1]1.456, F2=1-
[α-0.1026-1]1.531,ξ is the correction constant for considering shearing stress direction, and α is pore components, and n is
Ligancy, νma、GmaRespectively indicate matrix Poisson's ratio, modulus of shearing, effective pressure Peff=(ρg-ρw) gh, ρg、ρwRespectively indicate base
Matter density, fluid density, g indicate that acceleration of gravity, h indicate depth of stratum.
7. the method according to claim 4 for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version, feature exist
In, the method also includes:
School is carried out to the elastic rock computation model and electrical rock computation model according to log data or core data
Just；
Correspondingly, described calculated according to the elastic rock computation model determines that the corresponding elastic parameter data of grid node includes
The corresponding elastic parameter data of grid node is calculated using the elastic rock computation model after correction, it is described according to electrical rock
Stone computation model calculates the corresponding electrical parameter data of grid node including the use of the electrical rock computation model after correction
Calculate the corresponding electrical parameter data of grid node.
8. the method according to claim 1 for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version, feature exist
In the reservoir parameter data according to apart from the corresponding grid node of minimum value determine the reservoir parameter number of the subpoint
According to, comprising:
Filter out the grid node that the subpoint is less than preset threshold at a distance from each grid node；
The grid node that distance is less than preset threshold is ranked up according to the size of distance, is obtained apart from the corresponding net of minimum value
Lattice node.
9. the method according to claim 1 for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version, feature exist
In the three dimensional elasticity and electrical parameter measured data for obtaining formation at target locations, comprising:
Obtain the earthquake data before superposition and electromagnetic data of formation at target locations；
Elastic parameter measured data is obtained according to the prestack earthquake data inversion, and is obtained according to the electromagnetic data inverting
Electrical parameter measured data.
10. the side according to claim 1 to 8 for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version
Method, which is characterized in that the three dimensional elasticity and electrical parameter measured data for obtaining formation at target locations, comprising:
The elastic parameter and electrical parameter for obtaining preset kind, according to the elastic parameter and electrical parameter phase to preset kind
Sensitivity analysis to reservoir parameter is as a result, it is preferred that go out sensibility elasticity parameter and sensitive electrical parameter；
Using the sensibility elasticity parameter and sensitive electrical parameter as the three dimensional elasticity of formation at target locations and electrical parameter, obtain
The three dimensional elasticity and electrical parameter measured data of formation at target locations.
11. a kind of device for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version characterized by comprising
Module is obtained, for obtaining the three dimensional elasticity and electrical parameter measured data of formation at target locations；
Projection module, for by the three dimensional elasticity and electrical parameter measured data to the three-dimensional bullet electricity rock object constructed in advance
Reason amount version is projected, and obtains subpoint, wherein the building of three-dimensional bullet electricity rock physics amount version includes: by three-dimensional reservoir parameter
Space carries out the processing of gridding subdivision, calculates the corresponding elasticity of grid node and electrical parameter data, obtains three-dimensional bullet electricity rock
Stone physical quantity version；
Reservoir parameter determining module, for calculating the subpoint at a distance from each grid node, according to corresponding apart from minimum value
The reservoir parameter data of grid node determine the reservoir parameter data of the subpoint, according to the reservoir parameter of the subpoint
Data determine the reservoir parameter data of the formation at target locations.
12. a kind of device for determining reservoir parameter based on three-dimensional bullet electricity rock physics amount version, which is characterized in that including processor and
For the memory of storage processor executable instruction, when described instruction is executed by the processor realize the following steps are included:
Obtain the three dimensional elasticity and electrical parameter measured data of formation at target locations；
The three dimensional elasticity and electrical parameter measured data are thrown to the three-dimensional bullet electricity rock physics amount version constructed in advance
Shadow obtains subpoint, wherein the building of three-dimensional bullet electricity rock physics amount version includes: that three-dimensional reservoir parameter space is carried out grid
Change subdivision processing, calculates the corresponding elasticity of grid node and electrical parameter data, obtain three-dimensional bullet electricity rock physics amount version；
The subpoint is calculated at a distance from each grid node, according to the reservoir parameter number apart from the corresponding grid node of minimum value
According to the reservoir parameter data of the determination subpoint；
The reservoir parameter data of the formation at target locations are determined according to the reservoir parameter data of the subpoint.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810412688.2A CN108761529B (en) | 2018-05-03 | 2018-05-03 | Method and device for determining reservoir parameters based on three-dimensional elastic-electric rock physical quantity plate |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810412688.2A CN108761529B (en) | 2018-05-03 | 2018-05-03 | Method and device for determining reservoir parameters based on three-dimensional elastic-electric rock physical quantity plate |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108761529A CN108761529A (en) | 2018-11-06 |
| CN108761529B true CN108761529B (en) | 2019-10-11 |
Family
ID=64009434
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810412688.2A Active CN108761529B (en) | 2018-05-03 | 2018-05-03 | Method and device for determining reservoir parameters based on three-dimensional elastic-electric rock physical quantity plate |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108761529B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10928536B2 (en) * | 2017-12-07 | 2021-02-23 | Saudi Arabian Oil Company | Mapping chemostratigraphic signatures of a reservoir with rock physics and seismic inversion |
| CN109613624B (en) * | 2019-01-09 | 2020-03-13 | 中国石油大学(华东) | Reservoir rock acoustoelectric property combined simulation method |
| CN111487685B (en) * | 2019-01-25 | 2022-12-02 | 中国石油天然气股份有限公司 | Method and device for calculating resistivity of rock reservoir and storage medium |
| CN110599594B (en) * | 2019-07-29 | 2021-07-20 | 成都理工大学 | Three-dimensional modeling method for rock physical structure |
| CN116593295B (en) * | 2023-07-19 | 2023-10-03 | 北京科技大学 | Method and device for improving acoustic emission positioning precision by utilizing rock anisotropic wave velocity |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102854531A (en) * | 2012-09-11 | 2013-01-02 | 中国石油天然气股份有限公司 | Multi-scale rock physical plate method and device for reservoir hydrocarbon detection |
| CN104516017A (en) * | 2013-09-29 | 2015-04-15 | 中国石油化工股份有限公司 | Carbonate rock physical parameter seismic inversion method |
| CN106597542A (en) * | 2016-11-23 | 2017-04-26 | 中国石油天然气股份有限公司 | Reservoir characteristic parameter prediction method and device |
| CN107367762A (en) * | 2017-06-27 | 2017-11-21 | 中国石油天然气股份有限公司 | Method and device for determining reservoir parameters |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3203270B1 (en) * | 2016-02-03 | 2024-04-10 | CGG Services SAS | Complex pore geometry modeling by continuously varying inclusions (ci) method for elastic parameter prediction using inclusion models |
-
2018
- 2018-05-03 CN CN201810412688.2A patent/CN108761529B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102854531A (en) * | 2012-09-11 | 2013-01-02 | 中国石油天然气股份有限公司 | Multi-scale rock physical plate method and device for reservoir hydrocarbon detection |
| CN104516017A (en) * | 2013-09-29 | 2015-04-15 | 中国石油化工股份有限公司 | Carbonate rock physical parameter seismic inversion method |
| CN106597542A (en) * | 2016-11-23 | 2017-04-26 | 中国石油天然气股份有限公司 | Reservoir characteristic parameter prediction method and device |
| CN107367762A (en) * | 2017-06-27 | 2017-11-21 | 中国石油天然气股份有限公司 | Method and device for determining reservoir parameters |
Non-Patent Citations (2)
| Title |
|---|
| 《岩石物理模型分析在致密砂岩含气储层预测中的应用》;王永刚等;《2015年物探技术研讨会》;20151231;全文 * |
| 《精细刻画南海古近系油藏的叠前地质统计学反演》;薛怀艳等;《石油天然气学报（江汉石油学院学报）》;20130228;第35卷(第2期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108761529A (en) | 2018-11-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108761529B (en) | Method and device for determining reservoir parameters based on three-dimensional elastic-electric rock physical quantity plate | |
| Torabi et al. | Fault zone architecture and its scaling laws: where does the damage zone start and stop? | |
| Matthäi et al. | Control-volume finite-element two-phase flow experiments with fractured rock represented by unstructured 3D hybrid meshes | |
| US20170316128A1 (en) | Method and system for characterizing fractures in a subsurface region | |
| US8452580B2 (en) | Method and system for using multiple-point statistics simulation to model reservoir property trends | |
| CN103282747B (en) | For producing the system and method for the renewal of geological model | |
| CN104297785B (en) | Lithofacies constrained reservoir physical property parameter inversion method and device | |
| US20180031719A1 (en) | Method and System for Generating a Subsurface Model | |
| Durlofsky et al. | Uncertainty quantification for subsurface flow problems using coarse-scale models | |
| CN104755965A (en) | Methods and systems for determining minimum porosity for presence of clathrates in sediment | |
| CN107942404B (en) | Method and device for determining cracks | |
| US20210149067A1 (en) | Methods and Systems of Generating Simulation Grids for Different Structural Interpretations | |
| CN108983285A (en) | A kind of a variety of source wavefield analogy methods and device based on moment tensor | |
| US11409024B2 (en) | Methods and systems for generating simulation grids via zone by zone mapping from design space | |
| CN109856674A (en) | Engineering dessert evaluating method and device | |
| US10718188B2 (en) | Method for evaluation of fluid transport properties in heterogenous geological formation | |
| CN109870719A (en) | A kind of carbonate rock thin tight reservoir distribution determination method, apparatus and system | |
| Doveton | Theory and applications of vertical variability measures from Markov chain analysis | |
| CN109709602A (en) | Remote detection sound wave offset imaging method, device and system | |
| McKenna et al. | Fractured continuum approach to stochastic permeability modeling | |
| Taleghani et al. | Post-treatment assessment of induced fracture networks | |
| CN108710149A (en) | Thin reservoir prediction method and device | |
| CN109001804B (en) | Method, device and system for determining effective force based on three-dimensional seismic data | |
| CN109358364B (en) | Method, device and system for establishing underground river reservoir body geological model | |
| Jang et al. | The oil production performance analysis using discrete fracture network model with simulated annealing inverse method |
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 |