CN107367762A - A kind of method and device for determining reservoir parameter - Google Patents
A kind of method and device for determining reservoir parameter Download PDFInfo
- Publication number
- CN107367762A CN107367762A CN201710501254.5A CN201710501254A CN107367762A CN 107367762 A CN107367762 A CN 107367762A CN 201710501254 A CN201710501254 A CN 201710501254A CN 107367762 A CN107367762 A CN 107367762A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- rock
- mrow
- wave
- reservoir parameter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Abstract
The embodiment of the present application discloses a kind of method and device for determining reservoir parameter.Methods described includes：Obtain the attribute information of drilling well opening position rock in the purpose reservoir；Based on the attribute information of the rock, the target rock physicses template associated with reservoir parameter, P-S wave velocity ratio, the p-wave impedance of drilling well opening position in the purpose reservoir is established；The earthquake data before superposition of the purpose reservoir is obtained, based on the earthquake data before superposition, determines the inversion result data volume of the purpose reservoir；The inversion result data volume includes multiple data points；Based on the target rock physicses template, seed number strong point is determined from the inversion result data volume；Based on the target rock physicses template and the seed number strong point, reservoir parameter value corresponding to each data point in the inversion result data volume is determined.The technical scheme that the application provides, the degree of accuracy of identified reservoir parameter can be improved.
Description
Technical field
The application is related to technical field of geophysical exploration, more particularly to a kind of method and device for determining reservoir parameter.
Background technology
The parameters such as the porosity of reservoir, water saturation are for oil gas storage capacity and the oil and gas reserves progress to reservoir
The important parameter of evaluation.During seismic prospecting, by the geological information that is richly stored with the geological data that is gathered, institute
With researcher gathered geological data would generally be utilized to be predicted these reservoir parameters.
The main process of the method for the determination reservoir parameter of generally use is in the prior art：First, purpose reservoir is obtained
The attribute information of middle rock, for example, the modulus of elasticity of each mineral constituent and percent by volume, the porosity of rock become in rock
Change velocity of longitudinal wave, shear wave velocity and density log curve of scope and rock etc.；According to the attribute information of rock, mesh is established
The reservoir parameter of reservoir and the rock physics relations of elastic parameter, wherein, elastic parameter include velocity of longitudinal wave, shear wave velocity and
Density, and make rock physicses template using the rock physics relations；For example, according to the velocity of longitudinal wave of rock, shear wave velocity and
Density, the speed ratio and wave impedance of rock is calculated, coordinate system is established as transverse and longitudinal coordinate using speed ratio and wave impedance respectively, profit
With the rock physics relations established, isometric percentage curve and isoporosity curve is drawn on the coordinate system, these are bent
Line forms the rock physicses template；Then, the geological data of purpose reservoir is obtained, the geological data is carried out at prestack inversion
Reason, obtains the elastic parameter at multiple position of stratum in purpose reservoir；Finally, based on multiple fief layers in obtained purpose reservoir
The elastic parameter of opening position, using by the way of tabling look-up from rock physicses template determine purpose reservoir in correspond to stratum opening position
Reservoir parameter.
Inventor has found that at least there are the following problems in the prior art：Storage is waited when each in the rock physicses template drawn
When the spacing of layer parameter curve is larger, for example, isoporosity curve etc., according to a certain position of stratum in resulting purpose reservoir
The elastic parameter at place, data point corresponding to the speed ratio and wave impedance being calculated may the reservoir parameter curve such as each it
Between, not on reservoir parameter curve is waited, now, more difficult it may obtain the position of stratum from the rock physicses template exactly
The reservoir parameter at place.
The content of the invention
The purpose of the embodiment of the present application is to provide a kind of method and device for determining reservoir parameter, with storage determined by raising
The degree of accuracy of layer parameter.
In order to solve the above technical problems, it is so that the embodiment of the present application, which provides a kind of method and device for determining reservoir parameter,
Realize：
A kind of method for determining reservoir parameter, including：
Obtain the attribute information of drilling well opening position rock in purpose reservoir；Rock in the purpose reservoir includes at least one
Kind mineralogical composition；The attribute information is used for the feature for characterizing the rock；
Based on the attribute information of the rock, establish with the reservoir parameter of drilling well opening position in the purpose reservoir, in length and breadth
The wave velocity target rock physicses template more associated than, p-wave impedance；
The earthquake data before superposition of the purpose reservoir is obtained, based on the earthquake data before superposition, determines the purpose reservoir
Inversion result data volume；The inversion result data volume includes：Multiple data points；One data point is included at a position of stratum
P-S wave velocity ratio and p-wave impedance；
Based on the target rock physicses template, seed number strong point is determined from the inversion result data volume；
Based on the target rock physicses template and the seed number strong point, determine each in the inversion result data volume
Reservoir parameter value corresponding to individual data point.
In preferred scheme, the attribute information includes：Reservoir parameter change scope, containing fluid saturation, well logging compressional wave speed
Write music line, well logging shear wave velocity curve, well logging density curve, with least one mineral constituent distinguish corresponding initial elasticity
Modulus, the volume basis with the corresponding initial hole length-width ratio of at least one mineral constituent difference and various mineral constituents
Compare excursion.
It is described to be based on the rock physicses template and the seed number strong point in preferred scheme, determine the inversion result
Reservoir parameter value corresponding to each data point in data volume, including：
Based on the rock physicses template, it is determined that distinguishing corresponding reservoir parameter value with each seed number strong point, indulging
Shear wave velocity ratio and p-wave impedance；
Based on each seed number strong point corresponding reservoir parameter value, P-S wave velocity ratio and p-wave impedance respectively,
And P-S wave velocity ratio and p-wave impedance corresponding to a data point in the inversion result data volume, determine the data point pair
The reservoir parameter value answered.
It is described to be based on distinguishing corresponding reservoir parameter value, in length and breadth velocity of wave with each seed number strong point in preferred scheme
P-S wave velocity ratio corresponding to a data point and compressional wave resistance in degree ratio and p-wave impedance, and the inversion result data volume
It is anti-, reservoir parameter value corresponding to the data point is determined, including：
Based on each seed number strong point corresponding reservoir parameter value, P-S wave velocity ratio and p-wave impedance respectively,
And P-S wave velocity ratio and p-wave impedance corresponding to a data point in the inversion result data volume, using Power Interpolation
Method determines reservoir parameter value corresponding to the data point.
In preferred scheme, determine that reservoir corresponding to a data point is joined in the inversion result data volume using following formula
Numerical value：
In formula, ZjRepresent reservoir parameter value, z corresponding to j-th of data point in the inversion result data volumeiRepresent institute
State the reservoir parameter value at i-th of seed number strong point in seed number strong point, λijIt is described to represent that i-th of seed number strong point acts on
Weight in j-th of data point, M represent the number at the seed number strong point.
In preferred scheme, characterize i-th of seed number strong point using following formula and act in j-th of data point
Weight：
In formula, c represents scale factor, xiRepresent p-wave impedance corresponding to i-th of seed number strong point, yiRepresent institute
State P-S wave velocity ratio corresponding to i-th of seed number strong point, xjRepresent p-wave impedance corresponding to j-th of seed number strong point, yj
Represent P-S wave velocity ratio corresponding to j-th of seed number strong point.
It is described to be based on the earthquake data before superposition in preferred scheme, the inversion result data volume of the purpose reservoir is determined,
Including：
Inversion procedure is carried out to the earthquake data before superposition using the method for prestack Simultaneous Inversion, obtains the purpose reservoir
Formation at target locations opening position velocity of longitudinal wave, shear wave velocity and density；The formation at target locations position is appointed in the purpose reservoir
One position of stratum；
According to velocity of longitudinal wave, shear wave velocity and the density of the formation at target locations opening position, the formation at target locations position is determined
The P-S wave velocity ratio and p-wave impedance at place；
According to the P-S wave velocity ratio and p-wave impedance of the formation at target locations opening position, the inverting of the purpose reservoir is determined
Result data body.
It is described to be based on the rock physicses template in preferred scheme, determine seed number from the inversion result data volume
Strong point, including：
Determine isometric percentage curve of the rock physicses template with waiting reservoir ginseng in the inversion result data volume
Number of targets strong point corresponding to the P-S wave velocity ratio and p-wave impedance of the point of intersection of number curve, and using the number of targets strong point as
The seed number strong point.
A kind of device for determining reservoir parameter, described device include：Attribute information acquisition module, rock physicses template are established
Module, inversion result determining module, seed number strong point determining module and reservoir parameter value determining module；Wherein,
The attribute information acquisition module, for obtaining the attribute information of drilling well opening position rock in purpose reservoir；It is described
Rock in purpose reservoir includes at least one mineralogical composition；The attribute information includes：Reservoir parameter change scope, containing fluid
Saturation degree, well logging velocity of longitudinal wave curve, well logging shear wave velocity curve, well logging density curve and at least one mineral constituent
Initial elastic modulus corresponding to respectively, with least one mineral constituent corresponding initially hole length-width ratio and various respectively
The percent by volume excursion of mineral constituent；
The rock physicses template establishes module, for the attribute information based on the rock, establishes and is stored up with the purpose
The associated target rock physicses template of the reservoir parameter of drilling well opening position, P-S wave velocity ratio, p-wave impedance in layer；
The inversion result determining module, for obtaining the earthquake data before superposition of the purpose reservoir, based on the prestack
Geological data, determine the inversion result data volume of the purpose reservoir；The inversion result data volume includes：Multiple data points；
One data point includes the P-S wave velocity ratio and p-wave impedance at a position of stratum；
Seed number strong point determining module, for based on the rock physicses template, from the inversion result data volume
Middle determination seed number strong point；
The reservoir parameter value determining module, for based on the rock physicses template and the seed number strong point, it is determined that
Reservoir parameter value corresponding to each data point in the inversion result data volume.
In preferred scheme, the reservoir parameter value determining module includes：Seed information determining module and parameter value determine mould
Block；Wherein,
The seed information determining module, for based on the rock physicses template, it is determined that with each seed data
Reservoir parameter value, P-S wave velocity ratio and p-wave impedance corresponding to point difference；
The parameter value determining module, for based on each seed number strong point respectively corresponding reservoir parameter value,
In P-S wave velocity ratio and p-wave impedance, and the inversion result data volume P-S wave velocity ratio corresponding to a data point and
P-wave impedance, determine reservoir parameter value corresponding to the data point.
A kind of method and device for determining reservoir parameter of the embodiment of the present application, can be with based on the attribute information of the rock
Establish the target rock associated with reservoir parameter, P-S wave velocity ratio, the p-wave impedance of drilling well opening position in the purpose reservoir
Stone physical template；The earthquake data before superposition of the purpose reservoir can be obtained, based on the earthquake data before superposition, it may be determined that institute
State the inversion result data volume of purpose reservoir；, can be from the inversion result data volume really based on the rock physicses template
Determine seed number strong point；Based on the rock physicses template and the seed number strong point, it may be determined that the inversion result data volume
In reservoir parameter value corresponding to each data point.Reservoir parameter curve is not waited for data point in the rock physicses template
On situation, still can relatively accurately determine reservoir parameter value corresponding to the data point.
Brief description of the drawings
, below will be to embodiment or existing in order to illustrate more clearly of the embodiment of the present application or technical scheme of the prior art
There is the required accompanying drawing used in technology description to be briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments described in application, for those of ordinary skill in the art, do not paying the premise of creative labor
Under, other accompanying drawings can also be obtained according to these accompanying drawings.
Fig. 1 is a kind of flow chart for the embodiment of the method for determining reservoir parameter of the application；
Fig. 2 is the schematic diagram of rock physicses template in the embodiment of the present application；
Fig. 3 is the porosity diagrammatic cross-section of the purpose reservoir determined in the embodiment of the present application using the application method；
Fig. 4 is the composition structure chart for the device embodiment that the application determines reservoir parameter；
Fig. 5 is the composition structure chart that the application determines reservoir parameter value determining module in the device embodiment of reservoir parameter.
Embodiment
The embodiment of the present application provides a kind of method and device for determining reservoir parameter.
In order that those skilled in the art more fully understand the technical scheme in the application, it is real below in conjunction with the application
The accompanying drawing in example is applied, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described implementation
Example only some embodiments of the present application, rather than whole embodiments.It is common based on the embodiment in the application, this area
The every other embodiment that technical staff is obtained under the premise of creative work is not made, it should all belong to the application protection
Scope.
Fig. 1 is a kind of flow chart for the embodiment of the method for determining reservoir parameter of the application.As shown in figure 1, described determine storage
The method of layer parameter, comprises the following steps.
Step S101：Obtain the attribute information of drilling well opening position rock in purpose reservoir.
In one embodiment, the purpose reservoir can be the still undetermined reservoir of reservoir parameter.The purpose storage
Drilling well position in layer can refer to the drilling well by having been opened up in the purpose reservoir, be able to detect that the attribute letter of rock
The position of breath.The rock of drilling well opening position can include at least one mineralogical composition in the purpose reservoir.For example, the purpose
The rock of drilling well opening position can be calcite dolomite in reservoir, and it can include two kinds of mineral constituents, respectively dolomite and
Decompose stone.
The attribute information can be used for the feature for characterizing the rock.The attribute information can include：Reservoir parameter
Excursion, containing fluid saturation, well logging velocity of longitudinal wave curve, well logging shear wave velocity curve, well logging density curve, with it is described extremely
Lack initial elastic modulus corresponding to a kind of mineral constituent respectively, distinguish corresponding initial hole with least one mineral constituent
Length-width ratio and the percent by volume excursion of various mineral constituents.The initial elastic modulus can include：Initial volume
Modulus and original shear modulus.For example, the initial volume modulus of dolomite and original shear modulus can be 95 respectively in rock
Gigapascal (GPa) and 45GPa, the initial volume modulus of stone and original shear modulus are decomposed in rock can be respectively
76.8GPa and 32GPa.Initial hole length-width ratio corresponding to dolomite can be 0.08 in rock, be decomposed in rock corresponding to stone
Initial hole length-width ratio can be 0.06.The percent by volume excursion that stone is decomposed in rock is 0~100 percentage.It is described
Reservoir parameter can be porosity, and the reservoir parameter change may range from 2~20 percentages.
In one embodiment, can by way of drilling well and data acquisition mode, obtain the purpose reservoir
Rock core and log data.The log data can include：The WELL LITHOLOGY curve of drilling well opening position in the purpose reservoir,
Well logging porosity curve, well logging resistivity curve, well logging velocity of longitudinal wave curve, well logging shear wave velocity curve and well logging density curve
Deng.Physical Property Analysis processing is carried out to the rock core, can obtain with least one mineral constituent dividing in the attribute information
Not corresponding initial elastic modulus, and distinguish corresponding initial hole length-width ratio with least one mineral constituent.It is based on
The WELL LITHOLOGY curve, well logging porosity curve and well logging resistivity curve, it may be determined that reservoir is joined in the attribute information
Number excursion.Based on the WELL LITHOLOGY curve, it may be determined that the volume basis of various mineral constituents in the attribute information
Compare excursion.Based on well logging resistivity curve and well logging porosity curve, it may be determined that satisfy in the attribute information containing fluid
And degree.The fluid saturation that contains can be water saturation, oil saturation or gas saturation.
Step S102：Based on the attribute information of the rock, the reservoir with drilling well opening position in the purpose reservoir is established
The associated target rock physicses template of parameter, P-S wave velocity ratio, p-wave impedance.
Distinguish corresponding initial elastic modulus with least one mineral constituent in attribute information based on the rock,
With at least one mineral constituent corresponding initial hole length-width ratio, and containing fluid saturation, can establishing described respectively
The initial association of the elastic parameter of drilling well opening position and reservoir parameter, each mineral constituent percent by volume is closed in purpose reservoir
System.The elastic parameter can include：Velocity of longitudinal wave, shear wave velocity and density.Reservoir in attribute information based on the rock
Parameter variation range, the percent by volume excursion of each mineral constituent, and the initial association relation, it may be determined that institute
State the initial elasticity parameter curve for opening position saturated rock of being logged well in purpose reservoir.The initial elasticity parameter curve can wrap
Include：Initial velocity of longitudinal wave curve, initial shear wave velocity curve and initial density curve.In this way, with least one mineral group
A point initial elastic modulus corresponding to a difference and corresponding initial void is grown respectively with least one mineral constituent
Wide ratio, it can correspond to and obtain an initial elasticity parameter curve.It can calculate respectively initial in the initial elasticity parameter curve
Initial shear wave in first diversity factor of velocity of longitudinal wave curve and the well logging velocity of longitudinal wave curve, the initial elasticity parameter curve
Second diversity factor of rate curve and the well logging shear wave velocity curve, and initial density in the initial elasticity parameter curve
3rd diversity factor of curve and the well logging density curve.When first diversity factor, second diversity factor and the described 3rd
When diversity factor is respectively less than default diversity factor threshold value, can using initial elastic modulus corresponding to the initial elasticity parameter curve as
Object elastic modulus, and using initial hole length-width ratio corresponding to the initial elasticity parameter curve as target hole length and width
Than.In this way, can obtain the rock of the purpose reservoir with least one mineral constituent corresponding object elastic respectively
Modulus, and distinguish corresponding target hole length-width ratio with least one mineral constituent.
In one embodiment, it is right respectively with least one mineral constituent in the attribute information based on the rock
The initial elastic modulus answered, distinguish corresponding initial hole length-width ratio with least one mineral constituent, and satisfy containing fluid
And degree, establish the elastic parameter of drilling well opening position and reservoir parameter, each mineral constituent percent by volume in the purpose reservoir
Initial association relation.It can specifically include, divide in the attribute information based on the rock with least one mineral constituent
Not corresponding initial elastic modulus, it may be determined that the modulus of elasticity of Rock Matrix associates with each mineral constituent percent by volume
Relation.It is for instance possible to use following formula represent the modulus of elasticity of Rock Matrix and the pass of each mineral constituent percent by volume
Connection relation：
In formula (1),N represents the kind number of the mineral constituent of rock, Vk
Represent the percent by volume of kth kind mineral constituent in the rock, MkThe initial elasticity of kth kind mineral constituent in rock is stated in expression
Modulus, MsRepresent the modulus of elasticity of the Rock Matrix.Wherein, when modulus of elasticity is bulk modulus, M can be replaced with K；When
When modulus of elasticity is modulus of shearing, M can be replaced with μ.Based in the purpose reservoir each mineral constituent of rock it is default
Density, it may be determined that the incidence relation of the density of Rock Matrix and each mineral constituent percent by volume.Under it is for instance possible to use
State formula and represent the density of Rock Matrix and the incidence relation of each mineral constituent percent by volume：
In formula (2), ρsThe density of the Rock Matrix is represented, N represents the kind number of the mineral constituent of rock, VkRepresent institute
State the percent by volume of kth kind mineral constituent in rock, ρkRepresent the pre-set density of kth kind mineral constituent in rock.It can adopt
The incidence relation of the modulus of elasticity of rock matrix and the modulus of elasticity of Rock Matrix is represented with following formula：
In formula (3), KdAnd μdThe bulk modulus and modulus of shearing of the rock matrix are represented respectively, and N represents the ore deposit of rock
The kind number of thing component, fkRepresent to preset pore volume, P in the rock corresponding to the mineral constituent of kth kindskAnd QskInstitute is represented respectively
State bulk modulus polarization factor and modulus of shearing polarization factor corresponding to the mineral constituent of kth kind in rock, bulk modulus polarization because
Son initial hole length-width ratio corresponding with kth kind mineral constituent, Rock Matrix bulk modulus are associated, modulus of shearing polarization factor
Initial hole length-width ratio corresponding with kth kind mineral constituent, Rock Matrix modulus of shearing are associated, KsAnd μsDescribed in representing respectively
The bulk modulus and modulus of shearing of Rock Matrix, KkAnd μkVolume mould corresponding to the mineral constituent of kth kind in the rock is represented respectively
Amount and modulus of shearing,Following formula can be used to represent the modulus of elasticity and rock matrix of saturated rock
Modulus of elasticity, in the purpose reservoir reservoir parameter of drilling well opening position incidence relation：
In formula (4), KsatAnd μsatThe bulk modulus and modulus of shearing of saturated rock, K are represented respectivelydAnd μdRepresent respectively
The bulk modulus and modulus of shearing of the rock matrix, KsThe bulk modulus of the Rock Matrix is represented, φ represents the reservoir
Parameter, it can be porosity,Wherein, KhAnd KwThe bulk modulus of oil and water, S are represented respectivelyw
Water saturation in the attribute information is represented, e represents default index, and its span can be 1~5.It can use following
Formula represents the pass of the reservoir parameter of drilling well opening position in the density of saturated rock and the density of Rock Matrix, the purpose reservoir
Connection relation：
ρsat=(1- φ) ρs+φρh(1-Sw)+φρwSw (5)
In formula (5), ρsatRepresent the density of the saturated rock, ρsThe density of the Rock Matrix, ρhAnd ρwTable respectively
Show the density of oil and water, SwWater saturation in the attribute information is represented, φ represents the reservoir parameter, and it can be hole
Degree.Following formula can be used to represent the velocity of longitudinal wave of saturated rock and density, the modulus of elasticity of saturated rock of saturated rock
Incidence relation：
In formula (6), vp,satRepresent the velocity of longitudinal wave of saturated rock, KsatAnd μsatThe volume mould of saturated rock is represented respectively
Amount and modulus of shearing, ρsatRepresent the density of the saturated rock.Following formula can be used to represent the shear wave velocity of saturated rock
The incidence relation of the modulus of elasticity of density, saturated rock with saturated rock：
In formula (7), vs,satRepresent the shear wave velocity of saturated rock, μsatRepresent the modulus of shearing of saturated rock, ρsatTable
Show the density of the saturated rock.Based on above-mentioned formula (1)~(7), drilling well opening position in the purpose reservoir can be established
Elastic parameter and reservoir parameter, the initial association relation of each mineral constituent percent by volume.
Based on the object elastic modulus and target space length-width ratio, and it is described contain fluid saturation, according to upper
Step is stated, the elastic parameter of drilling well opening position and reservoir parameter, each mineral constituent volume in the purpose reservoir can be established
The target association relation of percentage.Reservoir parameter change scope in attribute information based on the rock, each mineral constituent
Percent by volume excursion, and the target association relation, it may be determined that opening position rock of being logged well in the purpose reservoir
Object elastic parameter value, i.e. target velocity of longitudinal wave, target shear wave velocity and target density.Based on the well logging opening position rock
The object elastic parameter value of stone, it may be determined that the P-S wave velocity ratio and p-wave impedance of the well logging opening position.In this way, the survey
The P-S wave velocity ratio of well opening position can be with reservoir parameter change scope, each mineral constituent in the attribute information of the rock
Percent by volume excursion be associated.It is described well logging opening position p-wave impedance can with the attribute information of the rock
Reservoir parameter change scope, the percent by volume excursion of each mineral constituent are associated.It can establish with p-and s-wave velocity
Than for ordinate and coordinate system that p-wave impedance is abscissa.It is reservoir parameter change scope in the attribute information of the rock, each
The percent by volume excursion of individual mineral constituent, a plurality of isometric percentage curve and more can be drawn on the coordinate system
The reservoir parameter curve such as bar.By a plurality of isometric percentage curve of the drafting and a plurality of it can wait reservoir parameter curve composition
Template as the target rock physicses template.For example, Fig. 2 is the schematic diagram of rock physicses template in the embodiment of the present application.
Abscissa and ordinate are respectively p-wave impedance and P-S wave velocity ratio in Fig. 2.It is a plurality of in Fig. 2 to wait calcite percent by volume bent
Line, such as V2=20%, V2=40%, V2=60%, V2=80%, V2=100% isometric percentage curve, with the hole such as a plurality of
Porosity curve, such as φ=2%, φ=6%, φ=10%, φ=14%, the isoporosity curve of φ=16%, the template of composition
For the target rock physicses template.Wherein, adjacent two percents by volume waited between calcite percent by volume curve are poor
It is worth identical, the porosity difference between adjacent two isoporosity curves is identical.
Step S103：The earthquake data before superposition of the purpose reservoir is obtained, based on the earthquake data before superposition, it is determined that described
The inversion result data volume of purpose reservoir.
The inversion result data volume can include：Multiple data points.One data point can include at a position of stratum
P-S wave velocity ratio and p-wave impedance.
By way of seismic prospecting and data acquisition, the earthquake data before superposition of the purpose reservoir can be obtained.Can be with
Inversion procedure is carried out to the earthquake data before superposition using the method for prestack Simultaneous Inversion, with obtaining the target of the purpose reservoir
Velocity of longitudinal wave, shear wave velocity and the density of layer opening position.The formation at target locations position can be in the purpose reservoir anyly
Layer position.According to velocity of longitudinal wave, shear wave velocity and the density of the formation at target locations opening position, it may be determined that the formation at target locations position
Put the P-S wave velocity ratio and p-wave impedance at place.According to the P-S wave velocity ratio and p-wave impedance of the formation at target locations opening position,
The inversion result data volume of the purpose reservoir can be determined.For example, black round dot represents the inversion result data volume in Fig. 2
In data point, i.e. the wave impedance of indicated inverting, speed are than data volume sampling point in Fig. 2.
Step S104：Based on the rock physicses template, seed number strong point is determined from the inversion result data volume.
Specifically, it can determine that isometric percentage of the rock physicses template is bent in the inversion result data volume
The P-S wave velocity ratio of point of intersection of the line with waiting reservoir parameter curve and the corresponding number of targets strong point of p-wave impedance, and by the mesh
Data point is marked as the seed number strong point for example, grey round dot S in Fig. 21For a seed number strong point, i.e. seed in Fig. 2
Point.
For example, table 1 is identified multiple seed number strong points in the embodiment of the present application, and divide with multiple seed number strong points
Not corresponding p-wave impedance, P-S wave velocity ratio and porosity.
P-wave impedance, P-S wave velocity ratio and porosity corresponding to the seed number strong point of table 1
| Seed number strong point | P-wave impedance | P-S wave velocity ratio | Porosity |
| S1 | 13710.4 | 1.67956 | 0.1 |
| S2 | 14743.6 | 1.70074 | 0.08 |
| S3 | 16011.7 | 1.72744 | 0.16 |
| S4 | 17491.2 | 1.75599 | 0.04 |
| S5 | 19064.6 | 1.79098 | 0.02 |
| S6 | 20778.9 | 1.82689 | 0 |
| S7 | 18782.8 | 1.75783 | 0.02 |
| S8 | 16998 | 1.69337 | 0.04 |
Step S105：Based on the rock physicses template and the seed number strong point, the inversion result data volume is determined
In reservoir parameter value corresponding to each data point.
Specifically, based on the rock physicses template, it may be determined that with the corresponding storage respectively of each seed number strong point
Layer parameter value, P-S wave velocity ratio and p-wave impedance.Based on each seed number strong point respectively corresponding reservoir parameter value,
In P-S wave velocity ratio and p-wave impedance, and the inversion result data volume P-S wave velocity ratio corresponding to a data point and
P-wave impedance, it may be determined that reservoir parameter value corresponding to the data point.Repeat reservoir parameter corresponding to described determination data point
The step of value, it may be determined that reservoir parameter value corresponding to each data point in the inversion result data volume.
In one embodiment, based on each seed number strong point corresponding reservoir parameter value, in length and breadth ripple respectively
P-S wave velocity ratio corresponding to a data point and compressional wave resistance in speed ratio and p-wave impedance, and the inversion result data volume
It is anti-, the method for Power Interpolation can be used to determine reservoir parameter value corresponding to the data point.Specifically, following formula can be used
Determine reservoir parameter value corresponding to a data point in the inversion result data volume：
Z in formula (1)jRepresent reservoir parameter value, z corresponding to j-th of data point in the inversion result data volumeiRepresent
The reservoir parameter value at i-th of seed number strong point, λ in the seed number strong pointijRepresent that i-th of seed number strong point acts on institute
The weight in j-th of data point is stated,M is represented
The number at seed number strong point, c represent scale factor,xiRepresent corresponding to i-th of seed number strong point
P-wave impedance, yiRepresent P-S wave velocity ratio corresponding to i-th of seed number strong point, xjRepresent j-th of seed number strong point
Corresponding p-wave impedance, yjRepresent P-S wave velocity ratio corresponding to j-th of seed number strong point.
For example, for a data point P1 in the inversion result data volume, its corresponding p-wave impedance is 15000,
P-S wave velocity ratio is 1.9.Table 2 is as 8 seed number strong points are respectively acting on the weight on data point P1 in table 1.Utilize public affairs
Formula (8), it is 8.4 percentages that can obtain porosity corresponding to data point P1.
The weight that the seed number strong point of table 2 is acted on data point P1
| Seed number strong point | λi,P1 |
| S1 | 0.148 |
| S2 | 0.240 |
| S3 | 0.241 |
| S4 | 0.115 |
| S5 | 0.054 |
| S6 | 0.028 |
| S7 | 0.058 |
| S8 | 0.116 |
Fig. 3 is the porosity diagrammatic cross-section of the purpose reservoir determined in the embodiment of the present application using the application method.Fig. 3
Middle gray value represents porosity, and two black lines represent the position of actual A wells and B wells respectively.Reality is explored result and shown, A
The daily gas of well epimere and hypomere is respectively 880,000 sides and 3.1 ten thousand sides, and the daily gas of B wells epimere and hypomere is respectively 5.9 ten thousand sides
With 2.2 ten thousand sides.A well epimere daily gas is most, while the porosity highest in Fig. 3 square frames, and about 0.12 or so.In this way, table
The degree of accuracy of the reservoir parameter of the bright purpose reservoir determined using the application method is higher.
The embodiment of the method for determining reservoir parameter, based on the attribute information of the rock, can be established and the mesh
Reservoir in the associated target rock physicses template of the reservoir parameter of drilling well opening position, P-S wave velocity ratio, p-wave impedance；Can
To obtain the earthquake data before superposition of the purpose reservoir, based on the earthquake data before superposition, it may be determined that the purpose reservoir
Inversion result data volume；Based on the rock physicses template, seed number strong point can be determined from the inversion result data volume；
Based on the rock physicses template and the seed number strong point, it may be determined that each data point in the inversion result data volume
Corresponding reservoir parameter value.In the case of grade on reservoir parameter curve, still may not be used in the rock physicses template for data point
Relatively accurately to determine reservoir parameter value corresponding to the data point.
Fig. 4 is the composition structure chart for the device embodiment that the application determines reservoir parameter.As shown in figure 4, described determine storage
The device of layer parameter can include：Attribute information acquisition module 100, rock physicses template establish module 200, inversion result determines
Module 300, seed number strong point determining module 400 and reservoir parameter value determining module 500.
The attribute information acquisition module 100, it can be used for the category for obtaining drilling well opening position rock in the purpose reservoir
Property information；Rock in the purpose reservoir can include at least one mineralogical composition；The attribute information can include：Reservoir
Parameter variation range, containing fluid saturation, well logging velocity of longitudinal wave curve, well logging shear wave velocity curve, well logging density curve and institute
State initial elastic modulus corresponding at least one mineral constituent respectively, with least one mineral constituent respectively it is corresponding initially
Hole length-width ratio and the percent by volume excursion of various mineral constituents.
The rock physicses template establishes module 200, can be used for the attribute information based on the rock, establish with it is described
The associated target rock physicses template of the reservoir parameter of drilling well opening position, P-S wave velocity ratio, p-wave impedance in purpose reservoir.
The inversion result determining module 300, it can be used for the earthquake data before superposition for obtaining the purpose reservoir, based on institute
Earthquake data before superposition is stated, determines the inversion result data volume of the purpose reservoir；The inversion result data volume can include：It is more
Individual data point；One data point can include the P-S wave velocity ratio and p-wave impedance at a position of stratum.
Seed number strong point determining module 400, it can be used for being based on the rock physicses template, from the inversion result
Seed number strong point is determined in data volume.
The reservoir parameter value determining module 500, it can be used for being based on the rock physicses template and the seed data
Point, determine reservoir parameter value corresponding to each data point in the inversion result data volume.
Fig. 5 is the composition structure chart that the application determines reservoir parameter value determining module in the device embodiment of reservoir parameter.
As shown in figure 5, reservoir parameter value determining module 500 includes in Fig. 4：Seed information determining module 510 and parameter value determining module
520。
The seed information determining module 510, can be used for being based on the rock physicses template, it is determined that with each kind
Reservoir parameter value, P-S wave velocity ratio and p-wave impedance corresponding to the difference of subnumber strong point.
The parameter value determining module 520, it can be used for being based on distinguishing corresponding reservoir with each seed number strong point
Ripple in length and breadth corresponding to a data point in parameter value, P-S wave velocity ratio and p-wave impedance, and the inversion result data volume
Speed ratio and p-wave impedance, determine reservoir parameter value corresponding to the data point.
The device embodiment for determining reservoir parameter is corresponding with the embodiment of the method for determining reservoir parameter, can be with
Realize the present processes embodiment, and the technique effect of adquisitiones embodiment.
In the 1990s, the improvement for a technology can clearly distinguish be on hardware improvement (for example,
Improvement to circuit structures such as diode, transistor, switches) or software on improvement (improvement for method flow).So
And as the development of technology, the improvement of current many method flows can be considered as directly improving for hardware circuit.
Designer nearly all obtains corresponding hardware circuit by the way that improved method flow is programmed into hardware circuit.Cause
This, it cannot be said that the improvement of a method flow cannot be realized with hardware entities module.For example, PLD
(Programmable Logic Device, PLD) (such as field programmable gate array (Field Programmable Gate
Array, FPGA)) it is exactly such a integrated circuit, its logic function is determined by user to device programming.By designer
Voluntarily programming comes a digital display circuit " integrated " on a piece of PLD, without asking chip maker to design and make
Special IC chip.Moreover, nowadays, substitution manually makes IC chip, this programming is also used instead mostly " patrols
Volume compiler (logic compiler) " software realizes that software compiler used is similar when it writes with program development,
And the source code before compiling also write by handy specific programming language, this is referred to as hardware description language
(Hardware Description Language, HDL), and HDL is also not only a kind of, but have many kinds, such as ABEL
(Advanced Boolean Expression Language)、AHDL(Altera Hardware Description
Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL
(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby
Hardware Description Language) etc., VHDL (Very-High-Speed are most generally used at present
Integrated Circuit Hardware Description Language) and Verilog2.Those skilled in the art
It will be apparent to the skilled artisan that only need method flow slightly programming in logic and being programmed into integrated circuit with above-mentioned several hardware description languages
In, it is possible to it is readily available the hardware circuit for realizing the logical method flow.
It is also known in the art that in addition to 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 controller with gate, switch, application specific integrated circuit, may be programmed
The form of logic controller and embedded microcontroller etc. realizes identical function.Therefore this controller is considered one kind
Hardware component, and it is used to realize that the device of various functions can also to be considered as the structure in hardware component to what is included in it.Or
Even, it not only can be able to will be the software module of implementation method for realizing that the device of various functions is considered as but also can be Hardware Subdivision
Structure in part.
Device that above-described embodiment illustrates, module, it can specifically be realized by computer chip or entity, or by with certain
The product of kind of function is realized.
For convenience of description, it is divided into various modules during description apparatus above with function to describe respectively.Certainly, this is being implemented
The function of each module can be realized in same or multiple softwares and/or hardware during application.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
Realized by the mode of software plus required general hardware platform.Based on such understanding, the technical scheme essence of the application
On the part that is contributed in other words to prior art can be embodied in the form of software product, in a typical configuration
In, computing device includes one or more processors (CPU), input/output interface, network interface and internal memory.The computer is soft
Part product can include some instructions make it that a computer equipment (can be personal computer, server, or network
Equipment etc.) perform method described in some parts of each embodiment of the application or embodiment.The computer software product can
To be stored in internal memory, internal memory may include the volatile memory in computer-readable medium, random access memory
(RAM) and/or the form such as Nonvolatile memory, such as read-only storage (ROM) or flash memory (flash RAM).Internal memory is computer
The example of computer-readable recording medium.Computer-readable medium includes permanent and non-permanent, removable and non-removable media can be by
Any method or technique come realize information store.Information can be computer-readable instruction, data structure, the module of program or its
His data.The example of the storage medium of computer includes, but are not limited to phase transition internal memory (PRAM), static RAM
(SRAM), dynamic random access memory (DRAM), other kinds of random access memory (RAM), read-only storage
(ROM), Electrically Erasable Read Only Memory (EEPROM), fast flash memory bank or other memory techniques, read-only optical disc are read-only
Memory (CD-ROM), digital versatile disc (DVD) or other optical storages, magnetic cassette tape, tape magnetic rigid disk storage or
Other magnetic storage apparatus or any other non-transmission medium, the information that can be accessed by a computing device available for storage.According to
Herein defines, and computer-readable medium does not include of short duration computer readable media (transitory media), such as modulation
Data-signal and carrier wave.
Each embodiment in this specification is described by the way of progressive, identical similar portion between each embodiment
Divide mutually referring to what each embodiment stressed is the difference with other embodiment.It is real especially for device
For applying example, because it is substantially similar to embodiment of the method, so description is fairly simple, related part is referring to embodiment of the method
Part explanation.
The application can be used in numerous general or special purpose computing system environments or configuration.Such as：Personal computer, clothes
Business device computer, handheld device or portable set, laptop device, multicomputer system, the system based on microprocessor, put
Top box, programmable consumer-elcetronics devices, network PC, minicom, mainframe computer including any of the above system or equipment
DCE etc..
The application can be described in the general context of computer executable instructions, such as program
Module.Usually, program module includes performing particular task or realizes routine, program, object, the group of particular abstract data type
Part, data structure etc..The application can also be put into practice in a distributed computing environment, in these DCEs, by
Task is performed and connected remote processing devices by communication network.In a distributed computing environment, program module can be with
In the local and remote computer-readable storage medium including storage device.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have it is many deformation and
Change is without departing from spirit herein, it is desirable to which appended claim includes these deformations and changed without departing from the application's
Spirit.
Claims (10)
- A kind of 1. method for determining reservoir parameter, it is characterised in that including：Obtain the attribute information of drilling well opening position rock in purpose reservoir；Rock in the purpose reservoir includes at least one ore deposit Thing composition；The attribute information is used for the feature for characterizing the rock；Based on the attribute information of the rock, establish and the reservoir parameter of drilling well opening position, in length and breadth velocity of wave in the purpose reservoir The associated target rock physicses template of degree ratio, p-wave impedance；The earthquake data before superposition of the purpose reservoir is obtained, based on the earthquake data before superposition, determines the anti-of the purpose reservoir Drill result data body；The inversion result data volume includes：Multiple data points；One data point includes vertical at a position of stratum Shear wave velocity ratio and p-wave impedance；Based on the target rock physicses template, seed number strong point is determined from the inversion result data volume；Based on the target rock physicses template and the seed number strong point, each number in the inversion result data volume is determined Reservoir parameter value corresponding to strong point.
- 2. a kind of method for determining reservoir parameter according to claim 1, it is characterised in that the attribute information includes： Reservoir parameter change scope, containing fluid saturation, well logging velocity of longitudinal wave curve, well logging shear wave velocity curve, well logging density curve, With at least one mineral constituent respectively corresponding initial elastic modulus, distinguish with least one mineral constituent it is corresponding Initial hole length-width ratio and the percent by volume excursion of various mineral constituents.
- 3. a kind of method for determining reservoir parameter according to claim 2, it is characterised in that described to be based on the rock thing Template and the seed number strong point are managed, determines reservoir parameter value corresponding to each data point in the inversion result data volume, Including：Based on the rock physicses template, it is determined that distinguishing corresponding reservoir parameter value, in length and breadth ripple with each seed number strong point Speed ratio and p-wave impedance；Based on each seed number strong point corresponding reservoir parameter value, P-S wave velocity ratio and p-wave impedance respectively, and P-S wave velocity ratio and p-wave impedance corresponding to a data point, are determined corresponding to the data point in the inversion result data volume Reservoir parameter value.
- 4. it is according to claim 3 it is a kind of determine reservoir parameter method, it is characterised in that it is described be based on it is each described Reservoir parameter value, P-S wave velocity ratio and p-wave impedance corresponding to the difference of seed number strong point, and the inversion result data volume In P-S wave velocity ratio and p-wave impedance corresponding to a data point, determine reservoir parameter value corresponding to the data point, including：Based on each seed number strong point corresponding reservoir parameter value, P-S wave velocity ratio and p-wave impedance respectively, and P-S wave velocity ratio and p-wave impedance corresponding to a data point in the inversion result data volume, using the method for Power Interpolation Determine reservoir parameter value corresponding to the data point.
- 5. a kind of method for determining reservoir parameter according to claim 4, it is characterised in that determine institute using following formula State reservoir parameter value corresponding to a data point in inversion result data volume：<mrow> <msub> <mi>Z</mi> <mi>j</mi> </msub> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <msub> <mi>&lambda;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <msub> <mi>z</mi> <mi>i</mi> </msub> </mrow>In formula, ZjRepresent reservoir parameter value, z corresponding to j-th of data point in the inversion result data volumeiRepresent the kind The reservoir parameter value at i-th of seed number strong point, λ in subnumber strong pointijRepresent that i-th of seed number strong point acts on the jth Weight in individual data point, M represent the number at the seed number strong point.
- 6. a kind of method for determining reservoir parameter according to claim 5, it is characterised in that characterize institute using following formula State the weight that i-th of seed number strong point is acted in j-th of data point：<mrow> <msub> <mi>&lambda;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>x</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>c</mi> <mn>2</mn> </msup> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>y</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> <mo>/</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <mrow> <mo>(</mo> <mfrac> <mn>1</mn> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>x</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>c</mi> <mn>2</mn> </msup> <msup> <mrow> <mo>(</mo> <mi>z</mi> <mo>-</mo> <msub> <mi>y</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow><mrow> <mi>c</mi> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>/</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <msub> <mi>y</mi> <mi>i</mi> </msub> </mrow>In formula, c represents scale factor, xiRepresent p-wave impedance corresponding to i-th of seed number strong point, yiRepresent described i-th P-S wave velocity ratio corresponding to individual seed number strong point, xjRepresent p-wave impedance corresponding to j-th of seed number strong point, yjRepresent P-S wave velocity ratio corresponding to j-th of seed number strong point.
- 7. a kind of method for determining reservoir parameter according to claim 1, it is characterised in that described based on the prestack Data are shaken, determine the inversion result data volume of the purpose reservoir, including：Inversion procedure is carried out to the earthquake data before superposition using the method for prestack Simultaneous Inversion, obtains the mesh of the purpose reservoir Mark velocity of longitudinal wave, shear wave velocity and the density at position of stratum；The formation at target locations position be in the purpose reservoir anyly Layer position；According to velocity of longitudinal wave, shear wave velocity and the density of the formation at target locations opening position, the formation at target locations opening position is determined P-S wave velocity ratio and p-wave impedance；According to the P-S wave velocity ratio and p-wave impedance of the formation at target locations opening position, the inversion result of the purpose reservoir is determined Data volume.
- 8. a kind of method for determining reservoir parameter according to claim 1, it is characterised in that described to be based on the rock thing Template is managed, seed number strong point is determined from the inversion result data volume, including：Determine isometric percentage curve of the rock physicses template with waiting reservoir parameter bent in the inversion result data volume Number of targets strong point corresponding to the P-S wave velocity ratio and p-wave impedance of the point of intersection of line, and using the number of targets strong point as described in Seed number strong point.
- 9. a kind of device for determining reservoir parameter, it is characterised in that described device includes：Attribute information acquisition module, rock thing Reason template establishes module, inversion result determining module, seed number strong point determining module and reservoir parameter value determining module；Wherein,The attribute information acquisition module, for obtaining the attribute information of drilling well opening position rock in purpose reservoir；The purpose Rock in reservoir includes at least one mineralogical composition；The attribute information includes：Reservoir parameter change scope, containing saturated with fluid Degree, well logging velocity of longitudinal wave curve, well logging shear wave velocity curve, well logging density curve and at least one mineral constituent are distinguished Corresponding initial elastic modulus, with least one mineral constituent corresponding initial hole length-width ratio and various mineral respectively The percent by volume excursion of component；The rock physicses template establishes module, for the attribute information based on the rock, establish with the purpose reservoir The target rock physicses template that reservoir parameter, P-S wave velocity ratio, the p-wave impedance of drilling well opening position are associated；The inversion result determining module, for obtaining the earthquake data before superposition of the purpose reservoir, based on the pre-stack seismic Data, determine the inversion result data volume of the purpose reservoir；The inversion result data volume includes：Multiple data points；One Data point includes the P-S wave velocity ratio and p-wave impedance at a position of stratum；Seed number strong point determining module, for based on the rock physicses template, from the inversion result data volume really Determine seed number strong point；The reservoir parameter value determining module, for based on the rock physicses template and the seed number strong point, it is determined that described Reservoir parameter value corresponding to each data point in inversion result data volume.
- 10. a kind of device for determining reservoir parameter according to claim 9, it is characterised in that the reservoir parameter value is true Cover half block includes：Seed information determining module and parameter value determining module；Wherein,The seed information determining module, for based on the rock physicses template, it is determined that with each seed number strong point point Not corresponding reservoir parameter value, P-S wave velocity ratio and p-wave impedance；The parameter value determining module, for based on each seed number strong point corresponding reservoir parameter value, in length and breadth respectively P-S wave velocity ratio and compressional wave corresponding to a data point in wave velocity ratio and p-wave impedance, and the inversion result data volume Impedance, determine reservoir parameter value corresponding to the data point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710501254.5A CN107367762B (en) | 2017-06-27 | 2017-06-27 | A kind of method and device of determining reservoir parameter |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710501254.5A CN107367762B (en) | 2017-06-27 | 2017-06-27 | A kind of method and device of determining reservoir parameter |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107367762A true CN107367762A (en) | 2017-11-21 |
| CN107367762B CN107367762B (en) | 2019-06-11 |
Family
ID=60305576
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710501254.5A Active CN107367762B (en) | 2017-06-27 | 2017-06-27 | A kind of method and device of determining reservoir parameter |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107367762B (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108761529A (en) * | 2018-05-03 | 2018-11-06 | 中国石油天然气股份有限公司 | The method and device of reservoir parameter is determined based on three-dimensional bullet electricity rock physics amount version |
| CN108873065A (en) * | 2018-05-10 | 2018-11-23 | 中国石油天然气股份有限公司 | Sandstone High-quality Reservoir prediction technique and device |
| CN109100802A (en) * | 2018-07-09 | 2018-12-28 | 中国石油天然气集团有限公司 | Reservoir determines method and apparatus |
| CN110426752A (en) * | 2019-08-20 | 2019-11-08 | 赛哲尔能源科技(北京)有限公司 | A kind of reservoir parameter inversion method and system based on petrophysical model |
| CN110609328A (en) * | 2019-09-11 | 2019-12-24 | 大庆油田有限责任公司 | Oil saturation prediction method based on well-seismic combination |
| CN111060964A (en) * | 2018-10-16 | 2020-04-24 | 中国石油天然气股份有限公司 | Method and device for determining stratum elasticity parameters |
| CN111077589A (en) * | 2018-10-19 | 2020-04-28 | 中国石油天然气股份有限公司 | Method and device for determining longitudinal-transverse wave velocity ratio of dry rock |
| WO2021096984A1 (en) * | 2019-11-15 | 2021-05-20 | Saudi Arabian Oil Company | Well logging to identify low resistivity pay zones in a subterranean formation using elastic attributes |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100191471A1 (en) * | 2009-01-23 | 2010-07-29 | Varel International Ind., L.P. | Method to determine rock properties from drilling logs |
| CN101887132A (en) * | 2009-05-15 | 2010-11-17 | 中国石油天然气股份有限公司 | Method for quantificationally predicting sandstone reservoir fluid saturation by combining well and seism |
| CN103376463A (en) * | 2012-04-13 | 2013-10-30 | 中国石油化工股份有限公司 | Inversion modeling method based on fault control |
| CN103412335A (en) * | 2013-08-20 | 2013-11-27 | 成都晶石石油科技有限公司 | Method for predicting reservoir stratums through seismic phase bodies |
| CN103760600A (en) * | 2014-01-07 | 2014-04-30 | 中国石油天然气股份有限公司 | Gas saturation inversion method |
| CN104391342A (en) * | 2014-12-09 | 2015-03-04 | 中国石油天然气股份有限公司 | Method and device for quantitatively forecasting reservoir permeability based on fracture porosity inversion |
| CN106054248A (en) * | 2016-07-15 | 2016-10-26 | 河海大学 | Earthquake rock physical inversion method based on large area tight reservoir |
| CN106597542A (en) * | 2016-11-23 | 2017-04-26 | 中国石油天然气股份有限公司 | Method and device for predicting reservoir characteristic parameters |
-
2017
- 2017-06-27 CN CN201710501254.5A patent/CN107367762B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100191471A1 (en) * | 2009-01-23 | 2010-07-29 | Varel International Ind., L.P. | Method to determine rock properties from drilling logs |
| CN101887132A (en) * | 2009-05-15 | 2010-11-17 | 中国石油天然气股份有限公司 | Method for quantificationally predicting sandstone reservoir fluid saturation by combining well and seism |
| CN103376463A (en) * | 2012-04-13 | 2013-10-30 | 中国石油化工股份有限公司 | Inversion modeling method based on fault control |
| CN103412335A (en) * | 2013-08-20 | 2013-11-27 | 成都晶石石油科技有限公司 | Method for predicting reservoir stratums through seismic phase bodies |
| CN103760600A (en) * | 2014-01-07 | 2014-04-30 | 中国石油天然气股份有限公司 | Gas saturation inversion method |
| CN104391342A (en) * | 2014-12-09 | 2015-03-04 | 中国石油天然气股份有限公司 | Method and device for quantitatively forecasting reservoir permeability based on fracture porosity inversion |
| CN106054248A (en) * | 2016-07-15 | 2016-10-26 | 河海大学 | Earthquake rock physical inversion method based on large area tight reservoir |
| CN106597542A (en) * | 2016-11-23 | 2017-04-26 | 中国石油天然气股份有限公司 | Method and device for predicting reservoir characteristic parameters |
Non-Patent Citations (2)
| Title |
|---|
| 李凌高 等: "面向叠前储层预测和油气检测的岩石物理分析新方法", 《内蒙古石油化工》 * |
| 闫玲玲 等: "叠后地质统计学反演在碳酸盐岩储层预测中的应用：以哈拉哈塘油田新垦区块为例", 《地学前缘》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108761529B (en) * | 2018-05-03 | 2019-10-11 | 中国石油天然气股份有限公司 | The method and device of reservoir parameter is determined based on three-dimensional bullet electricity rock physics amount version |
| CN108761529A (en) * | 2018-05-03 | 2018-11-06 | 中国石油天然气股份有限公司 | The method and device of reservoir parameter is determined based on three-dimensional bullet electricity rock physics amount version |
| CN108873065A (en) * | 2018-05-10 | 2018-11-23 | 中国石油天然气股份有限公司 | Sandstone High-quality Reservoir prediction technique and device |
| CN109100802B (en) * | 2018-07-09 | 2020-04-10 | 中国石油天然气集团有限公司 | Reservoir determination method and device |
| CN109100802A (en) * | 2018-07-09 | 2018-12-28 | 中国石油天然气集团有限公司 | Reservoir determines method and apparatus |
| CN111060964A (en) * | 2018-10-16 | 2020-04-24 | 中国石油天然气股份有限公司 | Method and device for determining stratum elasticity parameters |
| CN111060964B (en) * | 2018-10-16 | 2022-05-10 | 中国石油天然气股份有限公司 | Method and device for determining stratum elasticity parameters |
| CN111077589A (en) * | 2018-10-19 | 2020-04-28 | 中国石油天然气股份有限公司 | Method and device for determining longitudinal-transverse wave velocity ratio of dry rock |
| CN110426752A (en) * | 2019-08-20 | 2019-11-08 | 赛哲尔能源科技(北京)有限公司 | A kind of reservoir parameter inversion method and system based on petrophysical model |
| CN110609328A (en) * | 2019-09-11 | 2019-12-24 | 大庆油田有限责任公司 | Oil saturation prediction method based on well-seismic combination |
| CN110609328B (en) * | 2019-09-11 | 2022-03-22 | 大庆油田有限责任公司 | Oil saturation prediction method based on well-seismic combination |
| WO2021096984A1 (en) * | 2019-11-15 | 2021-05-20 | Saudi Arabian Oil Company | Well logging to identify low resistivity pay zones in a subterranean formation using elastic attributes |
| US11709287B2 (en) | 2019-11-15 | 2023-07-25 | Saudi Arabian Oil Company | Well logging to identify low resistivity pay zones in a subterranean formation using elastic attributes |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107367762B (en) | 2019-06-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107367762B (en) | A kind of method and device of determining reservoir parameter | |
| Cannon | Reservoir modelling: A practical guide | |
| Khan et al. | Machine learning derived correlation to determine water saturation in complex lithologies | |
| Burchette | Carbonate rocks and petroleum reservoirs: a geological perspective from the industry | |
| Deutsch et al. | The application of simulated annealing to stochastic reservoir modeling | |
| CN107817535A (en) | The determination method and apparatus of short lap | |
| Lia et al. | Uncertainties in reservoir production forecasts | |
| CN107942404B (en) | A kind of method and device in determining crack | |
| Hosseini et al. | Geostatistical modeling and spatial distribution analysis of porosity and permeability in the Shurijeh-B reservoir of Khangiran gas field in Iran | |
| CN106896406A (en) | The method and apparatus that High-quality Reservoir is predicted based on impedance Domain Properties | |
| Fachri et al. | Sensitivity of fluid flow to deformation-band damage zone heterogeneities: A study using fault facies and truncated Gaussian simulation | |
| CN109446735A (en) | A kind of generation method, equipment and the system of modeling logging data | |
| Kazemi et al. | Seismic history matching of nelson using time-lapse seismic data: An investigation of 4D signature normalization | |
| US11566518B2 (en) | Method and system for diagenesis-based rock classification | |
| CN105277979B (en) | The optimization method and device of a kind of seismic properties | |
| Bentley et al. | Future directions in reservoir modelling: new tools and ‘fit-for-purpose’workflows | |
| CN108035709A (en) | A kind of definite method and device of shale reservoir quality | |
| Alameedy et al. | Evaluating machine learning techniques for carbonate formation permeability prediction using well log data | |
| Chitralekha et al. | Application of the ensemble Kalman filter for characterization and history matching of unconventional oil reservoirs | |
| Suman et al. | Sensitivity study of rock-physics parameters for modeling time-lapse seismic response of Norne field | |
| CN106970423A (en) | A kind of method and device for determining shale gas stratigraphic structure model | |
| Sultan | New artificial neural network model for predicting the TOC from well logs | |
| Tariq | An intelligent functional network approach to develop mathematical model to predict sonic waves travel time for carbonate rocks | |
| Jalalalhosseini et al. | Predicting porosity by using seismic multi-attributes and well data and combining these available data by geostatistical methods in a South Iranian oil field | |
| Barman et al. | Fractured reservoir characterization using streamline-based inverse modeling and artificial intelligence tools |
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 |