CN109655936A - A kind of the elastic parameter calculation method and system of the replacement of clastic rock lithology - Google Patents
A kind of the elastic parameter calculation method and system of the replacement of clastic rock lithology Download PDFInfo
- Publication number
- CN109655936A CN109655936A CN201710942437.0A CN201710942437A CN109655936A CN 109655936 A CN109655936 A CN 109655936A CN 201710942437 A CN201710942437 A CN 201710942437A CN 109655936 A CN109655936 A CN 109655936A
- Authority
- CN
- China
- Prior art keywords
- rock
- lithology
- new
- porosity
- modulus
- 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
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
- G01V2210/6242—Elastic parameters, e.g. Young, Lamé or Poisson
-
- 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/66—Subsurface modeling
- G01V2210/661—Model from sedimentation process modeling, e.g. from first principles
Abstract
The invention belongs to earthquake rock physics research fields, specifically propose the elastic parameter calculation method and system of a kind of clastic rock lithology replacement, this method comprises: step 1: acquisition destination layer rock specimens carry out geology thin section identification, obtain rock basic parameter；Step 2: in conjunction with core data, log data and lithofacies, establishing the empirical relation of lithofacies, mineral content and porosity；Step 3: being demarcated and logged well with rock core, the lithologic character of destination layer is determined using well-log information；Step 4: according to rock basic parameter, porosity parameter and deposition characteristics, establishing petrophysical model；Step 5: the Rock Elastic Parameters that lithology replacement calculates new lithology are carried out, to calculate P- and S-wave velocity.More meet practical rock using the formation of elastic parameters that this method is extrapolated, improves the elastic parameter computational accuracy of the different lithology in extrapolation, the simulation problem of variation of lithological in solution reservoir quantitative forecast, reservoir description.
Description
Technical field
The invention belongs to petrophysics research fields, and in particular to a kind of method of Clastic Stratum of Country Rocks lithology replacement and be
System derives the displacement simulation rock elasticity of clastic rock pores degree, variation of lithological from deposition characteristics and rock physics theory
Calculation method of parameters.
Background technique
Earthquake rock physics are the bases of oil-gas exploration and development.Earthquake rock physics research mainly by rock physics modeling and
Rock physical analysis composition.Rock physical analysis is to analyze structure feature, the fluid of rock based on suitable petrophysical model
Influence to macroscopical response characteristic.In geophysics oil-gas exploration and development field, the petrophysical effect of earthquake is more and more important.
It is general to carry out fluid replacement analysis, it is assumed that reservoir porosity in Seismic Reservoir Prediction, reservoir description and monitoring process
It is constant, the fluid and its type of filling hole are changed, variation and its seismic facies Ying Te of Earthquake Resilient property are simulated
Sign establishes foundation for Seismic Reservoir Prediction and fluid identification.Under actual clastic rock reservoir condition, the porosity of reservoir and make
Rock mineral content is controlled by diagenesis, and porosity, rock matrix mineral constituent are related with lithofacies.To the reservoir characteristics of fixed well
Detailed research is all carried out, for the cross directional variations of reservoir, often can only be extrapolated and be predicted by seismic data.According to
Know that well data is extrapolated, analysis fluid, pressure variation caused by seismic signature change, the result and reality many times predicted
Border falls far short, and being primarily due to the difference of depositional environment, the porosity of reservoir, the charges in hole is horizontally variation
, the variation of storage and collection performance is caused, is not accounted in extrapolation process.For example, an alluvial fan sand body, middle part, root in fan body
Portion, alar part, rock matrix component and porosity are all different, if not considering the variation of rock constituents and porosity, cause outer
The seismic signature pushed away does not conform to the actual conditions, equally, same for mud shale stratum.
Rock physics are the strong means for carrying out reservoir extrapolation, still, at present to pore-fluid and its Change of types
(i.e. fluid replacement) pays much attention, and does to the analog study of variation of lithological few.There is scholar to carry out porosity variation
Rock physical modeling research, it is common practice to, do not consider skeleton variation, simply assume hole variation and its hole in
The variation of fluid, this is not inconsistent with actual reservoir situation.
In Rock physical analysis, the variation of porosity is not only considered, it is also necessary to consider the rock in porosity variation
Skeleton variation.
Summary of the invention
In order to solve problems in the prior art, the present invention studies a kind of rock that skeleton also changes while porosity variation
Property replacement calculation method.
The present invention be propose it is a kind of based on the elastic parameter that sedimentary principle and rock physics theory carry out lithology replacement
Calculation method, the opening relationships between reservoir parameter and seismic response, instructs Seismic Reservoir Prediction and fluid identification.
An invention according to the present invention provides a kind of elastic parameter calculation method of clastic rock lithology replacement, this method
Include:
Step 1: acquisition destination layer rock specimens carry out geology thin section identification and test, obtain rock basic parameter；Rock
Basic parameter includes rock forming mineral component and its content, sedimentary facies belt, porosity；
Step 2: in conjunction with core data, log data and lithofacies, establish lithofacies, mineral content and porosity through customs examination
System；
Step 3: being demarcated and logged well with rock core, the lithologic character of destination layer is determined using well-log information；
Step 4: according to rock basic parameter, porosity parameter and deposition characteristics, establishing petrophysical model；
Step 5: the Rock Elastic Parameters that lithology replacement calculates new lithology are carried out, to calculate P- and S-wave velocity.
Further, it is demarcated and is logged well with rock core, be utilized respectively well-log information and determine that reliable lithologic character includes rock mine
Object content, porosity, characteristic of fluid and its parameter.
Further, the empirical relation for establishing lithofacies, mineral content and porosity includes:
Establish quartz content variable quantity dQ, the clay content variable quantity dC of different location rock and the relationship of porosity d φ:
DQ=f (d φ, a) dC=f (d φ, b)
Wherein a, b are coefficient.
Further, establishing petrophysical model includes:
1. utilizing the concrete moduli and equivalent density of Voigt-Reuss-Hill average computation skeleton；
2. utilizing Wood model solution pore-fluid bulk modulus Kfl；
3. calculating dry Modulus of Rocks Kdry、μdry；
4. constructing fluid saturated rocks model using Gsssman equation, the elastic modulus of rock of saturation fluid is calculated.
Further, rock matrix concrete moduli are as follows:
Wherein, KV、μVIt is the bulk modulus and modulus of shearing that Voigt is obtained, K respectivelyR、μRIt is the body that Reuss is obtained respectively
Product module amount and modulus of shearing, Km、μmIt is the equivalent volume modulus and modulus of shearing that Hill is averagely obtained respectively；
Rock matrix equivalent density are as follows:
Wherein, fi、ρiThe respectively volume content of i-th of mineral constituent, density, N are the mineral constituent quantity of rock.
Further, with Batzle-Wang equation calculation gas, elasticity of fluid parameter, Wood model solution is recycled to obtain
Obtain pore-fluid bulk modulus Kfl, Wood model are as follows:
Wherein, Kw, Ko、KgThe respectively bulk modulus of water, oil, gas, Sw, So、SgRespectively water, oil, gas volume content, and
Sw+So+Sg=1.
Further, dry Modulus of Rocks K is determined by experimental determination datadry、μdry, or use following formula meter
It calculates:
Wherein, φ is rock porosity, and φ c is Critical porosity, Km、μmRespectively the bulk modulus of rock matrix and cut
Shear modulu.
Further, it is calculated by the following formula the elastic modulus of rock of saturation fluid:
μsat=μdry
ρsat=ρm(1-φ)+ρflφ
Wherein, K, μ are respectively bulk modulus and modulus of shearing, and ρ is density, and subscript sat, dry, fl, m respectively represent fluid
Saturation, dry rock, fluid and skeleton, φ are porosity.
Further, carry out lithology replacement include: change formation porosity, and according to established lithofacies, mineral content
The basic parameter that new lithology is determined with the empirical relation of porosity calculates the elasticity of new lithology according to the petrophysical model of step 4
Modulus and density.
Further, porosity variation is d φ=φnew- φ, φ are the original hole in stratum, φnewFor new lithology hole
Degree；
According to the relationship of quartz content variable quantity and porosity, the increment dQ of quartz content is calculated, according to clay content
Variable quantity and porosity relationship, calculate the calculating increment dC of clay content, obtain the quartz content Q of new lithologynew, clay contains
Measure Cnew, other mineral constituents remain unchanged；
According to new rock skeleton parameter and original fluid parameter, repeat step 4 1. 3. 4., calculate
The Rock Elastic Parameters K of new lithologynew、μnew、ρnew, finally calculate P- and S-wave velocity:
ρnew=ρm(1-φnew)+ρflφnew
Wherein, ρnew、Vpnew、VsnewDensity, velocity of longitudinal wave and the shear wave velocity of respectively new lithology, φnewFor new lithology
Porosity, ρflFor fluid density, ρmRepresent skeletal density, Knew、μnewFor the bulk modulus and modulus of shearing of new lithology.
Lithology replacement is completed, has extrapolated porosity and formation of elastic parameters when fluid all changes.
According to another aspect of the present invention, a kind of elastic parameter computing system of clastic rock lithology replacement, the system are provided
Include:
Memory is stored with computer executable instructions；
Processor, the processor run the computer executable instructions in the memory, execute following steps:
Step 1: acquisition destination layer rock specimens carry out geology thin section identification and test, obtain rock basic parameter, rock
Basic parameter includes rock forming mineral component and its content, sedimentary facies belt (microfacies), porosity；
Step 2: in conjunction with core data, log data and lithofacies, establish lithofacies, mineral content and porosity through customs examination
System；
Step 3: being demarcated and logged well with rock core, the lithologic character of destination layer is determined using well-log information；
Step 4: according to rock basic parameter, porosity parameter and deposition characteristics, establishing petrophysical model；
Step 5: the Rock Elastic Parameters that lithology replacement calculates new lithology are carried out, to calculate P- and S-wave velocity.
The present invention is to introduce the pass of sedimentary principle, lithology (component) and porosity relative to the innovative point of prior art
System, variation of lithological more meet practical rock, improve the elastic parameter computational accuracy of the different lithology in extrapolation, it is fixed to solve reservoir
The simulation problem of variation of lithological in amount prediction, reservoir description.
Detailed description of the invention
Disclosure illustrative embodiments are described in more detail in conjunction with the accompanying drawings, the disclosure above-mentioned and its
Its purpose, feature and advantage will be apparent, wherein in disclosure illustrative embodiments, identical reference label
Typically represent same parts.
Fig. 1 shows flow diagram according to the method for the embodiment of the present invention.
Fig. 2 shows that the sandstone layer porosity of S1 well according to an embodiment of the present invention is replaced with clay content curve and lithology
The P- and S-wave velocity and density curve of front and back are changed, wherein 1 (dotted line)-represents=8%, 2 (solid lines)-represent primary porosity, 3
(chain-dotted line)-represents=22%.
Specific embodiment
The preferred embodiment of the disclosure is more fully described below with reference to accompanying drawings.Although showing the disclosure in attached drawing
Preferred embodiment, however, it is to be appreciated that may be realized in various forms the disclosure without the embodiment party that should be illustrated here
Formula is limited.On the contrary, these embodiments are provided so that this disclosure will be more thorough and complete, and can be by the disclosure
Range is completely communicated to those skilled in the art.
The present invention be propose it is a kind of based on the elastic parameter that sedimentary principle and rock physics theory carry out lithology replacement
Calculation method belongs to earthquake rock physics research field.Method of the invention can comprise the following steps that
(1) destination layer rock specimens geology thin section identification and test, obtain rock basic parameter, and rock basic parameter includes
Rock forming mineral component and its content, sedimentary facies belt (microfacies), porosity；
(2) core data, log data and lithofacies are combined to establish lithofacies, quartz, clay content and porosity relationship；
(3) rock core calibration well logging, is utilized respectively well-log information and determines reliable lithologic character, including rock forming mineral content,
Porosity, characteristic of fluid and other parameters；
(4) suitable petrophysical model is constructed according to rock basic parameter, porosity basic parameter and deposition characteristics；
(5) relationship established according to second step changes lithology, porosity, carries out lithology replacement, obtains the bullet of different lithology
Property parameter.
More meet practical rock using the formation of elastic parameters that this method is extrapolated, improves the bullet of the different lithology in extrapolation
Property parameter computational accuracy, solve reservoir quantitative forecast, in reservoir description variation of lithological simulation problem.
According to embodiment of the present invention, a kind of elastic parameter calculation method of clastic rock lithology replacement is provided, it should
Method includes:
Step 1: acquisition destination layer rock specimens carry out geology thin section identification and test, obtain rock basic parameter；Rock
Basic parameter includes rock forming mineral component and its content, sedimentary facies belt (microfacies), porosity；
Step 2: in conjunction with core data, log data and lithofacies, establish lithofacies, mineral content and porosity through customs examination
System；
Step 3: being demarcated and logged well with rock core, the lithologic character of destination layer is determined using well-log information；
Step 4: according to rock basic parameter, porosity parameter and deposition characteristics, establishing petrophysical model；
Step 5: the Rock Elastic Parameters that lithology replacement calculates new lithology are carried out, to calculate P- and S-wave velocity.
Specifically, as shown in Figure 1, can be realized according to the method for the embodiment of the present invention by following steps:
(1) acquisition destination layer rock specimens carry out geology thin section identification and test, obtain rock basic parameter.
Rock basic parameter includes rock forming mineral component and its content, sedimentary facies belt (microfacies), porosity.
(2) core data, log data and lithofacies are combined to establish lithofacies, mineral content and porosity empirical relation.
Different location rock mineral constituent quartz, clay content variable quantity dQ, dC and porosity d φ relationship are as follows:
DQ=f (d φ, a) dC=f (d φ, b)
Wherein, a, b are coefficient.
(3) with rock core demarcate log well, be utilized respectively well-log information determine reliable lithologic character include rock forming mineral content,
Porosity, characteristic of fluid and its parameter.
(4) according to reservoir sandstone basic parameter, porosity etc., petrophysical model is established, is specifically included:
1. utilizing the concrete moduli K of Voigt-Reuss-Hill average computation skeletonMAnd μMWith equivalent density ρm；
2. utilizing Wood model solution pore-fluid bulk modulus Kfl；
3. calculating dry Modulus of Rocks Kdry、μdry；
4. constructing fluid saturated rocks model using Gsssman equation, the elastic modulus of rock of saturation fluid is calculated.
1. utilizing the concrete moduli and equivalent density of Voigt-Reuss-Hill average computation skeleton:
The upper bound K of concrete modulivWith lower bound KR:
Wherein, fi、KiItem and μiIt is volume content, bulk modulus and the modulus of shearing of i-th of mineral constituent respectively.So
Rock matrix concrete moduli are as follows:
Rock matrix equivalent density are as follows:
Wherein, fi、ρiThe respectively volume content of i-th of mineral constituent, density, N are the mineral constituent quantity of rock.
2. pore-fluid bulk modulus Kfl, Batzle-Wang equation calculation gas, elasticity of fluid parameter are generally used, then press
Wood model solution obtains.
Wood model are as follows:
Wherein, φ is rock porosity, Kw, Ko、KgThe respectively bulk modulus of water, oil, gas.Sw, So、SgRespectively water,
Oil, air volume content, and Sw+So+Sg=1.
The density of fluid are as follows:
ρfl=SWρw+SOρO+Sgρg
3. calculating dry Modulus of Rocks Kdry、μdry。
If there is experimental determination data, so that it may directly obtain these parameters.If general using such as without experimental data
Lower formula calculates:
Wherein, φ is rock porosity, φcFor Critical porosity, Km、μmRespectively the bulk modulus of rock matrix and cut
Shear modulu.
4. constructing fluid saturated rocks model with Gsssman equation, the elastic modulus of rock of saturation fluid is obtained.
μsat=μdry
ρsat=ρm(1-φ)+ρflφ
Wherein, K, μ are respectively bulk modulus and modulus of shearing, and ρ is density, and subscript sat, dry, fl, m respectively represent fluid
Saturation, dry rock, fluid and skeleton, φ are porosity.(5) carrying out lithology replacement includes: change formation porosity, and according to institute
The empirical relation of the lithofacies of foundation, mineral content and porosity determines the basic parameter of new lithology, according to the rock mould of step 4
Type calculates the elasticity modulus and density of new lithology.
Firstly, changing lithology (sedimentary micro) according to the relationship that (2) step is established, the original hole in stratum is φ, new lithology
Porosity is φnew, then porosity variation is d φ=φnew- φ, and rock porosity and rock matrix component (quartz content Q,
Clay content C) it synchronizes and changes, according to the relationship of quartz content variable quantity and porosity, calculate the increment of quartz content
DQ calculates the calculating increment dC of clay content according to clay content variable quantity and porosity relationship, with original (at wellhole)
Quartz content is added with clay content, obtains the quartz content Q of new lithologynew, clay content Cnew, other mineral constituents holdings are not
Become.
According to new rock skeleton parameter and original fluid parameter, repeat step 4 1. 3. 4., calculate new lithology
(pore-fluid bulk modulus KflRemain unchanged) Rock Elastic Parameters Knew、μnew、ρnew, finally calculate P- and S-wave velocity:
ρnew=ρm(1-φnew)+ρflφnew
Wherein, ρnew、Vpnew、VsnewDensity, velocity of longitudinal wave and the shear wave velocity of respectively new lithology, φnewFor new lithology
Porosity, ρflFor fluid density, ρmRepresent skeletal density, Knew、μnewFor the bulk modulus and modulus of shearing of new lithology.
Lithology replacement is completed, has extrapolated porosity and formation of elastic parameters when fluid all changes.
Another embodiment according to the present invention provides a kind of elastic parameter computing system of clastic rock lithology replacement, should
System includes:
Memory is stored with computer executable instructions；
Processor, the processor run the computer executable instructions in the memory, execute following steps:
Step 1: acquisition destination layer rock specimens carry out geology thin section identification, obtain rock basic parameter；
Step 2: in conjunction with core data, log data and lithofacies, establish lithofacies, mineral content and porosity through customs examination
System；
Step 3: being demarcated and logged well with rock core, the lithologic character of destination layer is determined using well-log information；
Step 4: according to rock basic parameter, porosity parameter and deposition characteristics, establishing petrophysical model；
Step 5: the Rock Elastic Parameters that lithology replacement calculates new lithology are carried out, to calculate P- and S-wave velocity.
A concrete application example is given below in the scheme and its effect of the embodiment of the present invention for ease of understanding.This field
It should be understood to the one skilled in the art that the example is only for the purposes of understanding the present invention, any detail is not intended to be limited in any way
The system present invention.
As shown in Fig. 2, the target zone of S1 well is 7.5 meters, porosity is 9%-12% or so, clay content 7%-18%
(left figure of Fig. 2).Assuming that porosity becomes 8% (dotted line) and 22% (chain-dotted line) respectively, increase according to quartz content, clay content
The relationship of amount and porosity, clay content increase dC=1.5d φ+0.366 on the basis of the original, and same quartz content is in original
Increase dQ, dQ=0.634-1.5d φ on the basis of coming, the stratum for obtaining that porosity is respectively 8% and 22% is replaced by lithology
(left side number of Fig. 2 plays the 2nd for velocity of longitudinal wave (left side number of Fig. 2 plays the 3rd figure), shear wave velocity (the rightmost figure of Fig. 2) and density
A figure).
More meet practical rock using the formation of elastic parameters that this method is extrapolated, improves the bullet of the different lithology in extrapolation
Property parameter computational accuracy, solve reservoir quantitative forecast, in reservoir description variation of lithological simulation problem.
The presently disclosed embodiments is described above, above description is exemplary, and non-exclusive, and
It is not limited to disclosed each embodiment.Without departing from the scope and spirit of illustrated each embodiment, for this skill
Many modifications and changes are obvious for the those of ordinary skill in art field.The selection of term used herein, purport
In principle, the practical application or to the technological improvement in market for best explaining each embodiment, or make the art its
Its those of ordinary skill can understand each embodiment disclosed herein.
Claims (10)
1. a kind of elastic parameter calculation method of clastic rock lithology replacement, which is characterized in that this method comprises:
Step 1: acquisition destination layer rock specimens carry out geology thin section identification and test, obtain rock basic parameter, rock is basic
Parameter includes rock forming mineral component and its content, sedimentary facies belt, porosity；
Step 2: in conjunction with core data, log data and lithofacies, establishing the empirical relation of lithofacies, mineral content and porosity；
Step 3: being demarcated and logged well with rock core, the lithologic character of destination layer is determined using well-log information；
Step 4: according to rock basic parameter, porosity parameter and deposition characteristics, establishing petrophysical model；
Step 5: the Rock Elastic Parameters that lithology replacement calculates new lithology are carried out, to calculate P- and S-wave velocity.
2. the elastic parameter calculation method of clastic rock lithology replacement according to claim 1, which is characterized in that establish rock
The empirical relation of phase, mineral content and porosity includes:
Establish quartz content variable quantity dQ, the clay content variable quantity dC of different location rock and the relationship of porosity d φ:
DQ=f (d φ, a) dC=f (d φ, b)
Wherein, a, b are coefficient.
3. the elastic parameter calculation method of clastic rock lithology replacement according to claim 1, which is characterized in that establish rock
Physical model includes:
1. utilizing the concrete moduli K of Voigt-Reuss-Hill average computation skeletonmAnd μmWith equivalent density ρm；
2. utilizing Wood model solution pore-fluid bulk modulus Kfl；
3. calculating dry Modulus of Rocks Kdry、μdry；
4. constructing fluid saturated rocks model using Gsssman equation, the elastic modulus of rock of saturation fluid is calculated.
4. the elastic parameter calculation method of clastic rock lithology replacement according to claim 3, which is characterized in that rock matrix
Concrete moduli are as follows:
Wherein, KV、μVIt is the bulk modulus and modulus of shearing that Voigt is obtained, K respectivelyR、μRIt is the volume mould that Reuss is obtained respectively
Amount and modulus of shearing, Km、μmIt is the equivalent volume modulus and modulus of shearing that Hill is averagely obtained respectively；
Rock matrix equivalent density are as follows:
Wherein, fi、ρiThe respectively volume content of i-th of mineral constituent, density, N are the mineral constituent quantity of rock.
5. the elastic parameter calculation method of clastic rock lithology replacement according to claim 3, which is characterized in that use
Batzle-Wang equation calculation gas, elasticity of fluid parameter recycle Wood model solution to obtain pore-fluid bulk modulus
Kfl, Wood model are as follows:
Wherein, Kw, Ko、KgThe respectively bulk modulus of water, oil, gas, Sw, So、SgRespectively water, oil, gas volume content, and Sw+So
+Sg=1.
6. the elastic parameter calculation method of clastic rock lithology replacement according to claim 3, which is characterized in that pass through experiment
Room determination data determines dry rock volume modulus Kdry, modulus of shearing μdry, or be calculated using the following equation:
Wherein, φ is rock porosity, φcFor Critical porosity, Km、μmThe respectively bulk modulus and shearing mould of rock matrix
Amount.
7. the elastic parameter calculation method of clastic rock lithology replacement according to claim 3, which is characterized in that by following
Formula calculates the elastic modulus of rock of saturation fluid:
μsat=μdry
ρsat=ρm(1-φ)+ρflφ
Wherein, K, μ are respectively bulk modulus and modulus of shearing, and ρ is density, and it is full that subscript sat, dry, fl, m respectively represent fluid
With dry rock, fluid and skeleton, φ is porosity.
8. the elastic parameter calculation method of clastic rock lithology replacement according to claim 1, which is characterized in that carry out lithology
Replacement includes: change formation porosity, and determines new rock according to the empirical relation of lithofacies, mineral content and the porosity established
The basic parameter of property calculates the elasticity modulus and density of new lithology according to the petrophysical model of step 4.
9. the elastic parameter calculation method of clastic rock lithology replacement according to claim 8, which is characterized in that porosity becomes
Turn to d φ=φnew- φ, φ are the original hole in stratum, φnewFor new lithology porosity；
According to the relationship of quartz content variable quantity and porosity, the increment dQ of quartz content is calculated, is changed according to clay content
Amount and porosity relationship calculates the calculating increment dC of clay content, obtains the quartz content Q of new lithologynew, clay content
Cnew, other mineral constituents remain unchanged；
According to new rock skeleton parameter and original fluid parameter, repeat step 4 1. 3. 4., calculate the rock of new lithology
Elastic parameter Knew、μnew、ρnew, finally calculate P- and S-wave velocity:
ρnew=ρm(1-φnew)+ρflφnew
Wherein, ρnew、Vpnew、VsnewDensity, velocity of longitudinal wave and the shear wave velocity of respectively new lithology, φnewFor the hole of new lithology
Degree, ρflFor fluid density, ρmRepresent skeletal density, Knew、μnewFor the bulk modulus and modulus of shearing of new lithology.
10. a kind of elastic parameter computing system of clastic rock lithology replacement, which is characterized in that the system includes:
Memory is stored with computer executable instructions；
Processor, the processor run the computer executable instructions in the memory, execute following steps:
Step 1: acquisition destination layer rock specimens carry out geology thin section identification and test, obtain rock basic parameter, rock is basic
Parameter includes rock forming mineral component and its content, sedimentary facies belt, porosity；
Step 2: in conjunction with core data, log data and lithofacies, establishing the empirical relation of lithofacies, mineral content and porosity；
Step 3: being demarcated and logged well with rock core, the lithologic character of destination layer is determined using well-log information；
Step 4: according to rock basic parameter, porosity parameter and deposition characteristics, establishing petrophysical model；
Step 5: the Rock Elastic Parameters that lithology replacement calculates new lithology are carried out, to calculate P- and S-wave velocity.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710942437.0A CN109655936B (en) | 2017-10-11 | 2017-10-11 | Elastic parameter calculation method and system for clastic rock lithology replacement |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710942437.0A CN109655936B (en) | 2017-10-11 | 2017-10-11 | Elastic parameter calculation method and system for clastic rock lithology replacement |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109655936A true CN109655936A (en) | 2019-04-19 |
| CN109655936B CN109655936B (en) | 2021-06-29 |
Family
ID=66108278
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710942437.0A Active CN109655936B (en) | 2017-10-11 | 2017-10-11 | Elastic parameter calculation method and system for clastic rock lithology replacement |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109655936B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112649870A (en) * | 2019-10-12 | 2021-04-13 | 中国石油化工股份有限公司 | Method and system for determining mineral elastic parameters in rock physical modeling |
| CN113866824A (en) * | 2020-06-30 | 2021-12-31 | 中国石油化工股份有限公司 | Method for determining elastic parameters of stratum between salts, storage medium and computer equipment |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102156297A (en) * | 2011-05-16 | 2011-08-17 | 中国石油大学(北京) | Fluid substitution method based on sandstone reservoir post-stack seismic data |
| US20150160368A1 (en) * | 2012-07-10 | 2015-06-11 | Statoil Petroleum As | Anisotropy parameter estimation |
| CN106569210A (en) * | 2016-10-18 | 2017-04-19 | 长安大学 | Soil moisture inversion method and soil moisture inversion device based on remote sensing |
| US10379255B2 (en) * | 2010-07-27 | 2019-08-13 | Exxonmobil Upstream Research Company | Inverting geophysical data for geological parameters or lithology |
-
2017
- 2017-10-11 CN CN201710942437.0A patent/CN109655936B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10379255B2 (en) * | 2010-07-27 | 2019-08-13 | Exxonmobil Upstream Research Company | Inverting geophysical data for geological parameters or lithology |
| CN102156297A (en) * | 2011-05-16 | 2011-08-17 | 中国石油大学(北京) | Fluid substitution method based on sandstone reservoir post-stack seismic data |
| US20150160368A1 (en) * | 2012-07-10 | 2015-06-11 | Statoil Petroleum As | Anisotropy parameter estimation |
| CN106569210A (en) * | 2016-10-18 | 2017-04-19 | 长安大学 | Soil moisture inversion method and soil moisture inversion device based on remote sensing |
Non-Patent Citations (1)
| Title |
|---|
| 马中高: "碎屑岩地震岩石物理学特征研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112649870A (en) * | 2019-10-12 | 2021-04-13 | 中国石油化工股份有限公司 | Method and system for determining mineral elastic parameters in rock physical modeling |
| CN113866824A (en) * | 2020-06-30 | 2021-12-31 | 中国石油化工股份有限公司 | Method for determining elastic parameters of stratum between salts, storage medium and computer equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109655936B (en) | 2021-06-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3571532B1 (en) | Systematic evaluation of shale plays | |
| US10572611B2 (en) | Method and system for characterizing fractures in a subsurface region | |
| US10846447B2 (en) | Method and system for stacking fracture prediction | |
| EP2960680B1 (en) | Fracturing and reactivated fracture volumes | |
| US10007015B2 (en) | Methods, systems and devices for predicting reservoir properties | |
| CN103258091B (en) | Unconventionaloil pool hides the method and device that net horizontal section three-dimensional mechanical models for rock mass is set up | |
| CN105089615B (en) | A kind of log data historical regression processing method based on reservoir model | |
| CA2913827C (en) | Methods, systems and devices for predicting reservoir properties | |
| CN104898161B (en) | Effective sandstone predicting method based on logging response simulator | |
| Angus et al. | Integrated hydro-mechanical and seismic modelling of the Valhall reservoir: a case study of predicting subsidence, AVOA and microseismicity | |
| US20110015912A1 (en) | Transport Property Data Calculated From Derivative Seismic Rock Property Data For Transport Modeling | |
| CN103256046A (en) | Unconventional oil and gas reservoir horizontal well section full-fracture-length fracturing parameter analog method and device | |
| CN105445791A (en) | Stratum aperture pressure prediction method based on variety earthquake attributes | |
| CN103914620B (en) | Method for computing distribution of opening spaces of fractures of fault rupture zones | |
| CN110501744A (en) | Hydrocarbon source rock organic carbon geophysics quantitative forecasting technique, device, equipment and storage medium | |
| Kim et al. | Petrophysical approach for S-wave velocity prediction based on brittleness index and total organic carbon of shale gas reservoir: a case study from Horn River Basin, Canada | |
| Kianoush et al. | Determining the drilling mud window by integration of geostatistics, intelligent, and conditional programming models in an oilfield of SW Iran | |
| CN109655936A (en) | A kind of the elastic parameter calculation method and system of the replacement of clastic rock lithology | |
| Dadashpour et al. | Porosity and permeability estimation by integration of production and time-lapse near and far offset seismic data | |
| Akhter et al. | Integrated interpretation with Gassmann fluid substitution for optimum field development of Sanghar area, Pakistan: a case study | |
| Sanei et al. | Building 1D and 3D static reservoir geomechanical properties models in the oil field | |
| Amjad et al. | Sembar Formation as an Unconventional Prospect: New Insights in Evaluating Shale Gas Potential Combined with Deep Learning | |
| CN109490988A (en) | Establish the method for being suitable for the rock physics new model of hard rock | |
| Duffaut et al. | Stress and fluid sensitivity in two North Sea oil fields—Comparing rock physics models with seismic observations | |
| Ali et al. | An integrated analysis of petrophysics, cross-plots and Gassmann fluid substitution for characterization of Fimkassar area, Pakistan: A case study |
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 |