CN109374497A - A kind of rock micropore structure test method - Google Patents
A kind of rock micropore structure test method Download PDFInfo
- Publication number
- CN109374497A CN109374497A CN201811269514.1A CN201811269514A CN109374497A CN 109374497 A CN109374497 A CN 109374497A CN 201811269514 A CN201811269514 A CN 201811269514A CN 109374497 A CN109374497 A CN 109374497A
- Authority
- CN
- China
- Prior art keywords
- rock
- modulus
- stiff
- fracture
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000011435 rock Substances 0.000 title claims abstract description 109
- 238000010998 test method Methods 0.000 title claims abstract description 13
- 230000003068 static effect Effects 0.000 claims abstract description 31
- 238000009825 accumulation Methods 0.000 claims abstract description 22
- 238000000034 method Methods 0.000 claims abstract description 17
- 238000011545 laboratory measurement Methods 0.000 claims description 17
- 238000010008 shearing Methods 0.000 claims description 8
- 238000005259 measurement Methods 0.000 claims description 4
- 230000009471 action Effects 0.000 claims description 2
- 238000004088 simulation Methods 0.000 claims description 2
- 230000008859 change Effects 0.000 abstract description 9
- 238000012512 characterization method Methods 0.000 abstract description 2
- 239000011148 porous material Substances 0.000 description 13
- 238000012360 testing method Methods 0.000 description 8
- 239000011800 void material Substances 0.000 description 6
- 238000011160 research Methods 0.000 description 4
- 239000004575 stone Substances 0.000 description 4
- 230000001419 dependent effect Effects 0.000 description 3
- 239000002245 particle Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000008187 granular material Substances 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000002706 hydrostatic effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000002922 simulated annealing Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000000547 structure data Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The present invention provides a kind of rock micropore structure test methods, this method comprises: establishing the relationship of characterization elastic modulus of rock and micropore structure；Estimate the aspect ratio of hard hole in rock；Calculate the accumulation fracture density of rock under each effective pressure；Establish the relationship between effective pressure and fracture density；Distribution characteristics based on static modulus of elasticity estimation fracture porosity；Calculate the distribution characteristics of static buildup fracture density and fracture porosity.The present invention fully considers the static modulus of elasticity feature more sensitive to pressure change, using the static modulus of elasticity replacement dynamic modulus of elasticity come the micropore structure distribution characteristics inside prediction rock, prediction result is more close to the truth of rock, it can be achieved that rock interior micropore structure is accurately portrayed.
Description
Technical field
The present invention relates to a kind of test method of rock interior fracture porosity more particularly to a kind of elasticity using rock
Modulus tests the characterization test method of its microscopic void distribution characteristics.
Background technique
Carry out impact analysis of the rock micropore structure to elastic modulus of rock, this has to subsurface reservoir feature is portrayed
Important directive significance (Tran et al, 2008).Ultrasonic velocity (dynamic elasticity mould of numerous scholars based on laboratory measurement
Amount) and stress-strain relation (static modulus of elasticity) come study the influence (King, 1966；Nur and Simmons, 1969；
Coyner,1984；Baud et al,2017).By research discovery dynamic modulus of elasticity in low effective pressure (effective pressure p=
Pc-Pp, wherein PcAnd PPConfining pressure and pore pressure respectively) under change it is obvious, this is because (crack is in length and breadth for aspect ratio
Than the ratio for being defined as fracture width half and length half) lesser crack is gradually closed with pressure increase, and then leads to rock
The resilient nature variation of stone is obvious.In essence, pressure-dependent rock elastic property and the close phase of micropore structure
It closes.Walsh (1965) has derived dry rock (pore pressure PPBe zero) bulk modulus (bulk modulus be elasticity modulus one
Kind, it is defined as the ratio of hydrostatic pressure and body strain) with confining pressure PcThe expression formula of variation finds crack (soft hole) to elasticity modulus
Influence it is more more significant than intergranular pore (hard hole), and estimated the fracture porosity of rock sample.Cheng and(1979) using KT equivalent medium mode (Kuster and1974) microcosmic hole in variety classes rock has been estimated
The porosity distribution characteristics of gap, i.e., with the hole of different aperture aspect ratio features, corresponding different porosity.Due to KT model
Can not estimate the elasticity modulus of high porosity rock, Tran etc. (2008) by introduce DEM equivalent medium mode (Berryman,
1992) improve Cheng andMethod.However, final calculated result is not unique, and the reliability of result
Prior model dependent on input.To solve this problem, Izumotani and Onozuka (2013) utilizes fast simulated annealing
Algorithm (a kind of algorithm based on probability) has estimated the porosity distribution characteristics of microscopic void.
Eberhart-Phillips etc. (1989) has derived seimic wave velocity and effective based on linear and exponential term combination
Rule-of-thumb relation between pressure.Zimmerman (1991) is pressed using three kinds of sample of sandstone of the relational expression approximate calculation
Contracting coefficient, and estimated the porosity distribution characteristics of microscopic void.Method based on Zimmerman, David and Zimmerman
(2012) using Mori-Tanaka it is theoretical (Mori and Tanaka, 1973；Benveniste, 1987) from the earthquake of dry rock
The porosity distribution characteristics (hereinafter referred to as DZ model) of microscopic void is obtained in wave velocity, and predicts to satisfy by Gassmann equation
With the p-and s-wave velocity of rock with the variation tendency of effective pressure.Though these methods easily realize, estimation result usually with reality
Border data misfit, research find based on the dynamic modulus of elasticity predict result it is often less than normal (Cheng and1979；
Pervukhina etc., 2010).
It is to sum up told, there are still following problems in existing technical research:
(1) dynamic modulus of elasticity does not fully consider the sensibility of pressure；
(2) true feelings of the accumulative fracture density, fracture porosity of dynamic modulus of elasticity prediction well below rock interior
Condition.
Summary of the invention
Goal of the invention: in view of the above problems, the present invention proposes a kind of rock micropore structure test method, this method base
The pore structure distribution characteristics of rock interior microscopic void is extracted in static modulus of elasticity, improves accumulative fracture density, crack
The test accuracy of the pore structures such as porosity.
Technical solution: the technical scheme adopted by the invention is that: a kind of rock micropore structure test method, including with
Lower step:
S1: establishing the relationship of elastic modulus of rock and micropore structure based on Mori-Tanaka theory, which includes
The porosity φ of hard hole in elastic modulus of rock shown in formula (1), formula (2) and rockstiffRelationship:
K in formulastiffAnd GstiffRespectively contain only the equivalent volume modulus and equivalent shear modulus of hard porous rocks, K0And G0
The bulk modulus and modulus of shearing of porous rocks are respectively free of, assumes that the hard hole in rock is elliposoidal, φ in formulastiffFor
The porosity of hard hole, P, Q are respectively the form factor of hard hole in rock, are had with the aspect ratio α of the hard hole of elliposoidal corresponding
Relationship；
And formula (3), (4) indicate elastic modulus of rock and rock in accumulate fracture density Γ relationship:
K in formulaeffAnd GeffThe respectively equivalent volume modulus and equivalent shear modulus of rock, vstiffIt is to contain only hard hole rock
The Poisson's ratio v of stonestiff=(3Kstiff-2Gstiff)/(6Kstiff+2Gstiff), Γ is accumulation fracture density；
Wherein, formula (1), P, Q in formula (2) is defined as:
Wherein, v=(3K0-2G0)/(6K0+2G0), g is only related with the aspect ratio α of hard hole, expression formula are as follows:
S2: stress-strain relation is passed through using laboratory routineThe method for measuring bulk modulus,
Measurement corresponding rock Equivalent Static bulk modulus when crack is closed under the action of effective pressure in rockAnd pass through
The aspect ratio of hard hole is calculated in formula (1)
S3: based on corresponding equivalent volume modulus K under each effective pressure of laboratory measurementeffOr shear modulus Geff, benefit
It is calculated with formula (3) or (4), corresponding accumulation fracture density Γ under each effective pressure p can be obtainedp,st(α)；Specifically include with
Lower process:
S31: based on corresponding equivalent volume modulus K under each effective pressure of laboratory measurementeffOr shear modulus Geff；
S32: by the aspect ratio of the hard hole obtained in S2 stepIt is calculated using formula (1) and formula (2)
Contain only the equivalent volume modulus K of hard porous rocksstiffWith equivalent shear modulus Gstiff；
S33: Poisson's ratio v is calculatedstiff=(3Kstiff-2Gstiff)/(6Kstiff+2Gstiff)；
S34: corresponding accumulation fracture density under each effective pressure p is calculated using formula (3) or formula (4)
Γp,st(α)。
S4: using fracture density with the variation exponentially decay law of effective pressureWhereinIt is
Initial fracture density when effective pressure is zero,It is one and pressure p with the pressure-constant of the order of magnitude, it is resulting according to S3
Corresponding accumulation fracture density Γ under each effective pressure pp,stParameter in (α) digital simulation relational expression, to obtain
The fracture density of rock with the variation of effective pressure matched curve or relational expression；
S5: pass through the relational expression of fracture porosityIt establishes between fracture porosity and aspect ratio
Relationship, to obtain the distribution characteristics of rock interior fracture porosity；Specifically include following procedure:
S51: being divided into n equal portions for laboratory measurement pressure range first, utilizes the relationship between Crack aspect ratio and pressure
Formula obtains corresponding Crack aspect ratio under each pressure；Wherein, the relational expression between the Crack aspect ratio and pressure are as follows:
In formulaIt is the Equivalent Static Young's modulus under high effective pressure, it can be by formula
Be calculated, Poisson's ratio byIt is calculated.
S52: obtaining pressure and accumulative fracture density relational expression secondly by S4, seeks corresponding accumulation crack under each pressure
Density；
S53: finally by the relational expression of fracture porosityEstablish fracture porosity and aspect ratio it
Between relationship, to obtain the distribution characteristics of rock interior fracture porosity.
S6: the distribution characteristics based on rock interior fracture porosity obtained by S5, to all fracture porosity and density
It sums, is calculated in rock and accumulates fracture porosity and accumulative fracture density curve, then the asymptotic value on the two curve
The fracture porosity and fracture density of rock interior respectively to be measured.
The utility model has the advantages that the present invention designs rock micropore structure extracting method, dynamic is replaced using static modulus of elasticity
Elasticity modulus estimates the porosity distribution characteristics of microscopic void.The present invention fully considers static modulus of elasticity to pressure change more
Sensitive feature replaces the dynamic modulus of elasticity using static modulus of elasticity to test the micropore structure of rock interior and be distributed spy
Sign, test result is more close to the truth of rock, it can be achieved that rock interior micropore structure is accurately portrayed.
Detailed description of the invention
Fig. 1 is that the present invention is based on the flow diagrams of static modulus of elasticity test rock micropore structure method；
Fig. 2 is the stress-strain curve of Navajo and Weber sandstone；
Fig. 3 is change curve of the Navajo sandstone p-and s-wave velocity with effective pressure；
Fig. 4 is change curve of the Weber sandstone p-and s-wave velocity with effective pressure；
Fig. 5 is the distribution characteristics of Navajo and Weber sandstone fracture porosity；
Fig. 6 is the accumulation fracture porosity distribution characteristics of Navajo and Weber sandstone；
Fig. 7 is the accumulation fracture density distribution characteristics of Navajo and Weber sandstone；
Fig. 8 is hard hole, crack and the total porosity of Navajo and Weber sandstone with the change curve of effective pressure.
Specific embodiment
Specific embodiments of the present invention will be described in further detail with reference to the accompanying drawings and examples.
The present invention provides a kind of rock micropore structure test methods, as shown in Figure 1, specifically includes the following steps:
Step S1: the relationship of elastic modulus of rock and micropore structure is established based on Mori-Tanaka theory, then rock
The equivalent elastic modulus of medium can be expressed as:
Wherein KstiffAnd GstiffRespectively contain only the equivalent volume and modulus of shearing of hard porous rocks, K0And G0Respectively rock
The volume and modulus of shearing of stone particle, φstiffFor the porosity of hole hard in rock, P, Q are respectively the form factor of hard hole,
It is related with the aspect ratio α of elliposoidal hole and the Poisson's ratio v of rock particles, it is defined as:
Wherein, v=(3K0-2G0)/(6K0+2G0),
Using the medium comprising hard hole as background phase, consider influence of the crack to rock elastic property, then rock medium
Equivalent elastic modulus can be expressed as:
Wherein, Keff、GeffThe respectively equivalent volume and modulus of shearing of rock, vstiff=(3Kstiff-2Gstiff)/
(6Kstiff+2Gstiff) it is the Poisson's ratio for containing only hard porous rocks, Γ is fracture density, it can be obtained by following formula:
Wherein, a is the mean radius in crack, and N is built-in the sum of the crack in cell cube (volume V), and bracket indicates
Average value.
Step S2: in view of Walsh (1965) give the clossing pressure relational expression p of holeclose=π E0α0/(4(1-
(v)2)), wherein E0=3K0(1-2v) is the Young's modulus of rock particles, α0Pore components when being p=0.Then in typical sand
Rock (E0=50GPa) in, the clossing pressure for the hole that aspect ratio is 0.01 is 500MPa, this is considerably beyond in laboratory measurement
Used pressure.Therefore, David and Zimmerman (2012) think the aspect ratio of hard hole in the range of 0.01 < α < 1,
The aspect ratio of soft hole is less than 0.01.
Only have hard hole since hole (crack) soft under high effective pressure is almost closed, in rock to exist.Therefore, may be used
The Equivalent Static bulk modulus of rock under the high effective pressure as measured by laboratoryBased on stress-strain relation
It calculates:
Wherein subscript " hp " is high effective pressure, and subscript " st " indicates that the value is the calculated value based on static modulus of elasticity,
Δ V is the variable quantity of cell cube volume V.Finally, be based on formula (1), using least square method, that is, the bulk modulus sought with it is quiet
State bulk modulusBetween error square be minimum when, corresponding most identical aspect ratio α, the aspect ratio of as hard hole
High effective pressure described in this step, for Navajo sandstone, range is generally taken as 80~100MPa,
Weber sandstone is then generally taken as 60~100MPa, for any one rock to be tested, can pass through fracture porosity
With the change curve of pressure, select when fracture porosity is near 0% corresponding effective pressure value as in this step,
Can directly take experiment test when institute attainable pressure maximum value, such as in this experiment the value be 100MPa, then utilize
It is finally fitted obtained fracture porosity and verifies whether the pressure value meets the requirements with the change curve of pressure.
Step S3: due to the pressure dependency and fracture density close relation of rock volume modulus, when known to fracture density
When, the elasticity modulus of rock can be obtained by formula (5) and (6), otherwise the elasticity modulus that may be based on laboratory measurement is estimated
Calculate the fracture density of rock.It therefore, can be by carrying out least square method calculating to the elasticity modulus measured under each effective pressure
Accumulation fracture density Γ under each pressurep,st(α).Specifically include following procedure:
S31: based on corresponding equivalent volume modulus K under each effective pressure of laboratory measurementeffOr shear modulus Geff；
S32: by the aspect ratio of the hard hole obtained in S2 stepIt is calculated using formula (1) and formula (2)
To the equivalent volume modulus K for containing only hard porous rocksstiffWith equivalent shear modulus Gstiff；
S33: Poisson's ratio v is calculatedstiff=(3Kstiff-2Gstiff)/(6Kstiff+2Gstiff)；
S34: being calculated using formula (3) or formula (4), using least square method, when utilization formula (3) or formula (4)
When square-error minimum between the equivalent volume modulus or modulus of shearing and experimental measurements of calculating, obtains corresponding accumulation and split
Gap density Γp,st(α)。
Step S4: domestic and foreign scholars propose (e.g., fracture density follows exponential damping law with the variation of effective pressure
Shapiro, 2003):
Wherein,It is the initial fracture density when effective pressure is zero,It is one and pressure p with the pressure of the order of magnitude
Force constant.Based on the fracture density gone out under each effective pressure p given in step S3, the fitting parameter in formula (9) can be obtained.
Step S5: in view of the crack (soft hole) in rock, when content is in Spectral structure in a certain range in length and breadth, because
This, when crack is closed with pressure increase, the initial aspect ratio of the minimum of each effective pressure p not closed all slits is
Wherein, Kst(Γp) it is rock Equivalent Static bulk modulus under effective pressure p, it can be calculated by formula (5).
Formula (9) is updated in formula (10):
Γ is carried out to formula (11)iTo ΓpIntegral between (α) can be obtained
Wherein,It is the Equivalent Static Poisson's ratio under high effective pressure.
In conjunction with formula (9) and (12), available pore componentsRelational expression between pressure p:
Wherein,It is the Equivalent Static Young's modulus under high effective pressure, is defined as
Poisson's ratio byIt is calculated.
When being gradually increased in view of effective pressure p, have the crack of different features in length and breadth with the variable quantity of pressure difference increment dp
It is identical.Therefore, when variable quantity, that is, dp of pressure is sufficiently small, the reduction amount of fracture density is attributable to aspect ratio and is less thanCrack closure caused by.Therefore rock interior fracture porosity can be obtained based on the static volume modulus of laboratory measurement
Distribution characteristics.Including following procedure:
S51: being divided into n equal portions for the pressure range of laboratory measurement first, utilizes the pass between pore components and pressure
It is formula, obtains corresponding aspect ratio under each pressure；
S52: obtaining pressure and accumulative fracture density relational expression secondly by S4, seeks corresponding accumulation crack under each pressure
Density；
S53: finally by the relational expression of fracture porosityEstablish fracture porosity and aspect ratio it
Between relationship, to obtain the distribution characteristics of rock interior fracture porosity.
Step S6: the stress-strain data based on sandstone laboratory measurement is calculated in rock and accumulates fracture porosity
With accumulative fracture density curve, then the asymptotic value on the two curve is fracture porosity and the crack for being respectively rock interior to be measured
Density.
Embodiment
It is further illustrated using Navajo and Weber sandstone as embodiment, the ess-strain used and ultrasonic velocity are real
Data are tested, from Navajo the and Weber sandstone data of Coyner (1984) laboratory measurement.Test measuring condition are as follows: dry
Dry state, effective pressure are up to 100MPa.The porosity of Navajo sandstone is 0.118, and the grain density of dry rock is
2316kg/m3, the volume and modulus of shearing of matrix granule are respectively 30GPa and 33GPa.The porosity of Weber sandstone is
0.095, the grain density of dry rock is 2392kg/m3, the bulk modulus of matrix granule and modulus of shearing be respectively 42GPa with
30GPa。
Fig. 2 is the stress-strain curve of Navajo and Weber sandstone, the stress proposed using Zimmerman (1991)
Fit correlation formulaIt is fitted.For Navajo sandstone, can be obtained
Δ V/V=0.0489p+0.0007 (1-e-p/0.014), wherein the unit of pressure p is GPa, goodness of fit R2=0.9998.And it is right
In Weber sandstone, then matched curve relational expression is Δ V/V=0.0563p+0.0059 (1-e-p/0.0111), goodness of fit R2=
0.9999.Experimental data is shown: due to the presence in crack, the load-deformation curve of two kinds of sandstone all shows nonlinear characteristic.
Fig. 3 and 4 is respectively change curve of the Navajo and Weber sandstone p-and s-wave velocity with effective pressure.It can find
When effect pressure is smaller, with the closure in crack, velocity variations are violent, and the research discovery of the trend and numerous scholars are consistent (such as
Carcione etc., 2003；Agersborg etc., 2008；Han etc., 2011；Asefet et al., 2013；Yurikov etc., 2017).
As step S1 and S2 is told, the aspect ratio of hard hole is estimated using the laboratory measurements under high effective pressure.Static bullet
The hard pore components for the Navajo and Weber sandstone that property modulus is estimated are respectively 0.18 and 0.08, and the dynamic modulus of elasticity
Estimate that hard pore components are respectively 0.28 and 0.11.
The fit correlation of pressure and fracture density can be obtained according to step S3 and S4.For Navajo sandstone, based on static state
The fit correlation formula that elasticity modulus obtains is Γp,st=0.8842e-p/0.0089, goodness of fit R2=0.9509；And based on dynamic
The relational expression that elasticity modulus obtains is Γp,dy=0.2007e-p/0.0237, goodness of fit R2=0.9894.For Weber sandstone,
The relational expression obtained based on static modulus is Γp,st=5.1001e-p/0.0089, goodness of fit R2=0.9887；Based on dynamic analog
The relational expression measured is Γp,dy=1.4596e-p/0.0154, goodness of fit R2=0.9951.
The distribution characteristics of Navajo and Weber sandstone fracture porosity can be obtained according to step S5.Using the prior art, base
In the ultrasonic velocity data of laboratory measurement, the distribution that pore components in rock are calculated based on the dynamic modulus of elasticity is special
Sign, accumulative fracture density and fracture porosity, and compared and analyzed with the calculated result of step S5.As shown in figure 5, real in figure
Line and double-crossed are respectively indicated based on the static calculated result with the dynamic modulus of elasticity.As figure shows, it is based on static volume modulus
The fracture porosity estimated is higher than the prediction of the dynamic modulus of elasticity.It estimates to obtain in Navajo sandstone based on static modulus of elasticity
The peak value of fracture porosity is 0.00035%, and corresponding aspect ratio is 0.00023；It is respectively 0.0023% He in Weber sandstone
0.00026.Likewise, the peak value for the Navajo sandstone fracture porosity estimated based on the dynamic modulus of elasticity and score in length and breadth
It is not 7 × 10-5% and 0.0006, Weber sandstone are respectively 6 × 10-4% and 0.00045.
Accumulation fracture porosity and the accumulative fracture density distribution spy of Navajo and Weber sandstone can be obtained according to step S6
Sign.According to DZ model, when accumulating fracture porosity or accumulative fracture density reaches asymptotic value, the asymptotic value of the two is respectively rock
Fracture porosity and fracture density inside stone.Using the prior art, the ultrasonic velocity data based on laboratory measurement are based on
The dynamic modulus of elasticity accumulates fracture porosity and accumulative fracture density to calculate in rock, and carries out with the calculated result of step S6
Comparative analysis, as shown in Figures 6 and 7.As shown, the accumulation fracture pore of the Navajo sandstone based on the estimation of static volume modulus
Degree is 0.084%, and accumulation fracture density is 0.88, and is 0.047% based on the fracture porosity of dynamic volume modulus, and crack is close
Degree is 0.2.For Weber sandstone, the accumulation fracture porosity based on the estimation of static volume modulus is 0.55%, and accumulation crack is close
Degree is 5.1, and is 0.25% based on the fracture porosity of dynamic volume modulus, fracture density 1.5.The result shows that based on quiet
The porosity and density in the accumulation crack of state bulk modulus estimation are all higher than the value that dynamic volume modulus obtains.
To illustrate that the rock pore structure data tested based on static modulus of elasticity more meet the true feelings of rock interior
Condition will be tested the pressure-dependent curve of resulting rock porosity and be compared with experimental measurements, and be based on dynamic bullet
The corresponding data that property modulus test obtains are compared, as shown in Figure 8.The porosity of hard hole is that have effect based on height in figure
The linear fit relationship of data is estimated to obtain (Pervukhina etc., 2010) under power.For Navajo sandstone, it is based on static elastic
The related coefficient of the estimated total porosity of Moduli data is 0.9933, and the related coefficient that the dynamic modulus of elasticity is estimated
0.9849.In addition, related coefficient is respectively 0.9820 and 0.9038 for Weber sandstone.The result shows that being based on static volume
The fracture porosity and experimental data of modulus estimation are more identical.
Claims (5)
1. a kind of rock micropore structure test method, which comprises the following steps:
S1: establishing the relationship of elastic modulus of rock and micropore structure based on Mori-Tanaka theory, which includes formula
(1), the porosity φ of elastic modulus of rock shown in formula (2) and hole hard in rockstiffRelationship:
K in formulastiffAnd GstiffRespectively contain only the equivalent volume modulus and equivalent shear modulus of hard porous rocks, K0And G0Respectively
Assume that the hard hole in rock is elliposoidal, φ for bulk modulus and modulus of shearing without porous rocks, in formulastiffFor rock
In hard hole porosity, P, Q are respectively the form factor of hard hole, have corresponding relationship with the aspect ratio α of the hard hole of elliposoidal；
And formula (3), (4) indicate elastic modulus of rock and rock in accumulate fracture density Γ relationship:
K in formulaeffAnd GeffThe respectively equivalent volume modulus and equivalent shear modulus of rock, vstiffContain only hard porous rocks
Poisson's ratio vstiff=(3Kstiff-2Gstiff)/(6Kstiff+2Gstiff), Γ is accumulation fracture density；
S2: stress-strain relation is passed through using laboratory routineThe method for measuring bulk modulus, measurement
The corresponding rock Equivalent Static bulk modulus when crack is closed under the action of effective pressure in rockAnd pass through formula
(1) aspect ratio of hard hole is calculated
S3: based on corresponding equivalent volume modulus K under each effective pressure of laboratory measurementeffOr shear modulus Geff, utilize public affairs
Formula (3) or (4) calculate, and corresponding accumulation fracture density Γ under each effective pressure p can be obtainedp,st(α)；
S4: using fracture density with the variation exponentially decay law of effective pressureWhereinIt is to have
Initial fracture density when pressure is zero is imitated,It is one and pressure p with the pressure-constant of the order of magnitude, it is resulting each according to S3
Corresponding accumulation fracture density Γ under effective pressure pp,stParameter in (α) digital simulation relational expression, to obtain rock
Fracture density with effective pressure variation matched curve or relational expression；
S5: pass through the relational expression of fracture porosityThe relationship between fracture porosity and aspect ratio is established,
To obtain the distribution characteristics of rock interior fracture porosity；
S6: the distribution characteristics based on rock interior fracture porosity obtained by S5 carries out all fracture porosity and density
Summation is calculated in rock and accumulates fracture porosity and accumulative fracture density curve, then the asymptotic value difference on the two curve
For the fracture porosity and fracture density of rock interior to be measured.
2. rock micropore structure test method according to claim 1, it is characterised in that: public affairs described in step S1
P, Q in formula (1), formula (2) is defined as:
Wherein, v=(3K0-2G0)/(6K0+2G0), g is only related with the aspect ratio α of hard hole, expression formula are as follows:
3. rock micropore structure test method according to claim 1, it is characterised in that: step S3 includes following mistake
Journey:
S31: based on corresponding equivalent volume modulus K under each effective pressure of laboratory measurementeffOr shear modulus Geff；
S32: by the aspect ratio of the hard hole obtained in S2 stepIt is calculated and contains only using formula (1) and formula (2)
The equivalent volume modulus K of hard porous rocksstiffWith equivalent shear modulus Gstiff；
S33: Poisson's ratio v is calculatedstiff=(3Kstiff-2Gstiff)/(6Kstiff+2Gstiff)；
S34: corresponding accumulation fracture density Γ under each effective pressure p is calculated using formula (3) or formula (4)p,st
(α)。
4. rock micropore structure test method according to claim 1, it is characterised in that: step S5 includes following mistake
Journey:
S51: being divided into n equal portions for laboratory measurement pressure range first, using the relational expression between Crack aspect ratio and pressure, obtains
Corresponding Crack aspect ratio under to each pressure；
S52: obtaining pressure and accumulative fracture density relational expression secondly by S4, and it is close to seek corresponding accumulation crack under each pressure
Degree；
S53: finally by the relational expression of fracture porosityIt establishes between fracture porosity and aspect ratio
Relationship, to obtain the distribution characteristics of rock interior fracture porosity.
5. rock micropore structure test method according to claim 4, it is characterised in that: split described in step S51
Relational expression between gap aspect ratio and pressure are as follows:
Wherein,It is the Equivalent Static Young's modulus under high effective pressure, it can be by formulaBe calculated, Poisson's ratio byIt calculates
It arrives.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811269514.1A CN109374497B (en) | 2018-10-29 | 2018-10-29 | Rock micro-pore structure testing method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811269514.1A CN109374497B (en) | 2018-10-29 | 2018-10-29 | Rock micro-pore structure testing method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109374497A true CN109374497A (en) | 2019-02-22 |
| CN109374497B CN109374497B (en) | 2020-12-11 |
Family
ID=65390484
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811269514.1A Active CN109374497B (en) | 2018-10-29 | 2018-10-29 | Rock micro-pore structure testing method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109374497B (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109827881A (en) * | 2019-02-26 | 2019-05-31 | 西南石油大学 | Rock hydration degree characterizing method and system |
| CN111208566A (en) * | 2020-03-04 | 2020-05-29 | 中国石油大学(北京) | Hole seam parameter inversion method and device based on SCA model and storage medium |
| CN111595752A (en) * | 2020-06-22 | 2020-08-28 | 中国科学技术大学 | Method for determining effective porosity of rock |
| CN113009565A (en) * | 2021-03-24 | 2021-06-22 | 中国石油大学(北京) | Method, device and equipment for determining seismic wave velocity parameters based on SCA model |
| CN113887150A (en) * | 2021-09-18 | 2022-01-04 | 河海大学 | Method for estimating length of characteristic jet flow of compact sandstone |
| CN113916738A (en) * | 2020-07-07 | 2022-01-11 | 中国石油化工股份有限公司 | Method for measuring pore compression coefficient of crack medium |
| CN115575505A (en) * | 2022-10-10 | 2023-01-06 | 四川大学 | Method for calculating longitudinal wave velocity and attenuation of rock under stress action condition |
Citations (24)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN2570772Y (en) * | 2002-10-11 | 2003-09-03 | 石油大学(北京) | Apparatus for measuring fine and close rock pore |
| CN201730615U (en) * | 2010-06-26 | 2011-02-02 | 大庆油田有限责任公司 | Evacuating saturation device for testing rock porosity |
| CN102096107A (en) * | 2009-12-09 | 2011-06-15 | 中国石油天然气股份有限公司 | Method for evaluating permeability of reservoir layer according to interval transit time and density inversed pore flat degree |
| CN202041426U (en) * | 2011-04-27 | 2011-11-16 | 浙江农林大学 | Measuring device for longitudinal shear elasticity modulus of wood |
| CN102445708A (en) * | 2010-10-14 | 2012-05-09 | 中国石油大学(北京) | Rich argillaceous sand stone three-dimension equivalent speed prediction model |
| CN102608013A (en) * | 2012-03-02 | 2012-07-25 | 河海大学 | Method for measuring porosity in piping development process |
| EP2490140A1 (en) * | 2010-08-19 | 2012-08-22 | PART Engineering GmbH | Computer-implemented method for determining the elasto-plastic mechanical behaviour of inhomogeneous anisotropic materials |
| CN102830051A (en) * | 2012-06-08 | 2012-12-19 | 河海大学 | Rock mass fracture occurrence recognizing system in oscillation testing system |
| CN102841046A (en) * | 2012-09-20 | 2012-12-26 | 中国石油大港油田勘探开发研究院 | Method for measuring rock porosity and measuring device |
| CN102854531A (en) * | 2012-09-11 | 2013-01-02 | 中国石油天然气股份有限公司 | Multi-scale rock physical charting method and device for detecting reservoir hydrocarbon |
| CN103091352A (en) * | 2013-01-28 | 2013-05-08 | 河海大学 | Multiscale prediction method of coefficients of thermal expansion of common cement paste in early stages |
| CN103105485A (en) * | 2013-01-28 | 2013-05-15 | 河海大学 | Hardened ordinary cement paste thermal expansion coefficient multiscale predication method |
| CN103207137A (en) * | 2013-04-08 | 2013-07-17 | 河海大学 | Device for fully automatically measuring penetration rate and void ratio of compact rock under dynamic confining pressure |
| CN103207138A (en) * | 2013-04-08 | 2013-07-17 | 河海大学 | Method for combined measurement of permeability and porosity of dense rock under dynamic confining pressure |
| CN103235338A (en) * | 2013-03-25 | 2013-08-07 | 中国石油大学(华东) | Rock fracture parameter inversion method |
| CN103576195A (en) * | 2013-10-28 | 2014-02-12 | 西北大学 | Method for forecasting fissured medium transverse wave velocity varying with pressure |
| CN103792581A (en) * | 2013-12-31 | 2014-05-14 | 刘致水 | Shale gas reservoir crustal stress logging prediction method based on rock physics model |
| CN104655501A (en) * | 2015-01-23 | 2015-05-27 | 同济大学 | Method for acquiring mechanical property of electrochemical deposition product |
| CN105510207A (en) * | 2016-01-29 | 2016-04-20 | 河海大学 | Method for determining sensitivity of crack rock mass permeability in all directions to porosity |
| US20160139084A1 (en) * | 2014-11-18 | 2016-05-19 | Kabushiki Kaisha Toshiba | Signal processing apparatus, server, detection system, and signal processing method |
| CN106054248A (en) * | 2016-07-15 | 2016-10-26 | 河海大学 | Earthquake rock physical inversion method based on large area tight reservoir |
| CN107832513A (en) * | 2017-11-01 | 2018-03-23 | 中科吉奥（北京）地质工程咨询有限责任公司 | Rock mass engineering project calculation method of parameters and system based on statistical mechanics of rockmass |
| CN107921178A (en) * | 2015-05-05 | 2018-04-17 | 哈佛大学的校长及成员们 | Tubular tissue constructs and Method of printing |
| CN108387931A (en) * | 2018-01-22 | 2018-08-10 | 中国石油天然气股份有限公司 | The determination method and apparatus of reservoir pore pressure |
-
2018
- 2018-10-29 CN CN201811269514.1A patent/CN109374497B/en active Active
Patent Citations (24)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN2570772Y (en) * | 2002-10-11 | 2003-09-03 | 石油大学(北京) | Apparatus for measuring fine and close rock pore |
| CN102096107A (en) * | 2009-12-09 | 2011-06-15 | 中国石油天然气股份有限公司 | Method for evaluating permeability of reservoir layer according to interval transit time and density inversed pore flat degree |
| CN201730615U (en) * | 2010-06-26 | 2011-02-02 | 大庆油田有限责任公司 | Evacuating saturation device for testing rock porosity |
| EP2490140A1 (en) * | 2010-08-19 | 2012-08-22 | PART Engineering GmbH | Computer-implemented method for determining the elasto-plastic mechanical behaviour of inhomogeneous anisotropic materials |
| CN102445708A (en) * | 2010-10-14 | 2012-05-09 | 中国石油大学(北京) | Rich argillaceous sand stone three-dimension equivalent speed prediction model |
| CN202041426U (en) * | 2011-04-27 | 2011-11-16 | 浙江农林大学 | Measuring device for longitudinal shear elasticity modulus of wood |
| CN102608013A (en) * | 2012-03-02 | 2012-07-25 | 河海大学 | Method for measuring porosity in piping development process |
| CN102830051A (en) * | 2012-06-08 | 2012-12-19 | 河海大学 | Rock mass fracture occurrence recognizing system in oscillation testing system |
| CN102854531A (en) * | 2012-09-11 | 2013-01-02 | 中国石油天然气股份有限公司 | Multi-scale rock physical charting method and device for detecting reservoir hydrocarbon |
| CN102841046A (en) * | 2012-09-20 | 2012-12-26 | 中国石油大港油田勘探开发研究院 | Method for measuring rock porosity and measuring device |
| CN103091352A (en) * | 2013-01-28 | 2013-05-08 | 河海大学 | Multiscale prediction method of coefficients of thermal expansion of common cement paste in early stages |
| CN103105485A (en) * | 2013-01-28 | 2013-05-15 | 河海大学 | Hardened ordinary cement paste thermal expansion coefficient multiscale predication method |
| CN103235338A (en) * | 2013-03-25 | 2013-08-07 | 中国石油大学(华东) | Rock fracture parameter inversion method |
| CN103207137A (en) * | 2013-04-08 | 2013-07-17 | 河海大学 | Device for fully automatically measuring penetration rate and void ratio of compact rock under dynamic confining pressure |
| CN103207138A (en) * | 2013-04-08 | 2013-07-17 | 河海大学 | Method for combined measurement of permeability and porosity of dense rock under dynamic confining pressure |
| CN103576195A (en) * | 2013-10-28 | 2014-02-12 | 西北大学 | Method for forecasting fissured medium transverse wave velocity varying with pressure |
| CN103792581A (en) * | 2013-12-31 | 2014-05-14 | 刘致水 | Shale gas reservoir crustal stress logging prediction method based on rock physics model |
| US20160139084A1 (en) * | 2014-11-18 | 2016-05-19 | Kabushiki Kaisha Toshiba | Signal processing apparatus, server, detection system, and signal processing method |
| CN104655501A (en) * | 2015-01-23 | 2015-05-27 | 同济大学 | Method for acquiring mechanical property of electrochemical deposition product |
| CN107921178A (en) * | 2015-05-05 | 2018-04-17 | 哈佛大学的校长及成员们 | Tubular tissue constructs and Method of printing |
| CN105510207A (en) * | 2016-01-29 | 2016-04-20 | 河海大学 | Method for determining sensitivity of crack rock mass permeability in all directions to porosity |
| CN106054248A (en) * | 2016-07-15 | 2016-10-26 | 河海大学 | Earthquake rock physical inversion method based on large area tight reservoir |
| CN107832513A (en) * | 2017-11-01 | 2018-03-23 | 中科吉奥（北京）地质工程咨询有限责任公司 | Rock mass engineering project calculation method of parameters and system based on statistical mechanics of rockmass |
| CN108387931A (en) * | 2018-01-22 | 2018-08-10 | 中国石油天然气股份有限公司 | The determination method and apparatus of reservoir pore pressure |
Non-Patent Citations (3)
| Title |
|---|
| BA JIN 等: "Compressional wave dispersion due to rock matrix stiffening by clay squirt flow", 《GEOPHYSICAL RESEARCAL LETTERS》 * |
| O"CONNELL R J 等: "Seismic velocities in dry and saturated cracked solids", 《JOURNAL OF GEOPHYSICAL RESEARCH》 * |
| 朱祖扬 等: "利用有限差分法模拟双孔介质波场传播", 《2016中国地球科学联合学术年会论文集（二十三）——专题46：地震波传播与成像》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109827881A (en) * | 2019-02-26 | 2019-05-31 | 西南石油大学 | Rock hydration degree characterizing method and system |
| CN109827881B (en) * | 2019-02-26 | 2022-05-31 | 西南石油大学 | Rock hydration degree characterization method and system |
| CN111208566A (en) * | 2020-03-04 | 2020-05-29 | 中国石油大学(北京) | Hole seam parameter inversion method and device based on SCA model and storage medium |
| CN111595752A (en) * | 2020-06-22 | 2020-08-28 | 中国科学技术大学 | Method for determining effective porosity of rock |
| CN111595752B (en) * | 2020-06-22 | 2021-07-06 | 中国科学技术大学 | Method for determining effective porosity of rock |
| CN113916738A (en) * | 2020-07-07 | 2022-01-11 | 中国石油化工股份有限公司 | Method for measuring pore compression coefficient of crack medium |
| CN113009565A (en) * | 2021-03-24 | 2021-06-22 | 中国石油大学(北京) | Method, device and equipment for determining seismic wave velocity parameters based on SCA model |
| CN113887150A (en) * | 2021-09-18 | 2022-01-04 | 河海大学 | Method for estimating length of characteristic jet flow of compact sandstone |
| CN113887150B (en) * | 2021-09-18 | 2022-05-06 | 河海大学 | Method for estimating length of characteristic jet flow of compact sandstone |
| CN115575505A (en) * | 2022-10-10 | 2023-01-06 | 四川大学 | Method for calculating longitudinal wave velocity and attenuation of rock under stress action condition |
| CN115575505B (en) * | 2022-10-10 | 2024-02-13 | 四川大学 | Calculation method for rock longitudinal wave velocity and attenuation under stress condition |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109374497B (en) | 2020-12-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109374497A (en) | A kind of rock micropore structure test method | |
| Shapiro | Elastic piezosensitivity of porous and fractured rocks | |
| David et al. | Pore structure model for elastic wave velocities in fluid‐saturated sandstones | |
| Lade et al. | Modelling the elastic behaviour of granular materials | |
| Prevost et al. | Constitutive relations for soil materials | |
| WO2008103986B1 (en) | Method and system for creating upscaled reservoir models | |
| CN108399270B (en) | Method for determining anisotropic shale proportion in shale stratum | |
| CN111399044A (en) | Reservoir permeability prediction method and device and storage medium | |
| Tarrahi et al. | Integration of microseismic monitoring data into coupled flow and geomechanical models with ensemble Kalman filter | |
| WO2013158873A2 (en) | System and method for calibrating permeability for use in reservoir modeling | |
| Polavarapu et al. | Some challenges of middle atmosphere data assimilation | |
| Sayers et al. | The effect of pore fluid on the stress-dependent elastic wave velocities in sandstones | |
| Oskay et al. | A survey of geotechnical system identification techniques | |
| Lissa et al. | Squirt flow in cracks with rough walls | |
| Edwards et al. | Characterizing the sediment bed in terms of resistance to motion: Toward an improved model of saltation thresholds for aeolian transport | |
| Jia et al. | Dynamic and cyclic properties of soils | |
| Sergeant et al. | Numerical modeling of iceberg capsizing responsible for glacial earthquakes | |
| Lin et al. | Two-stage covariance-based multisensing damage detection method | |
| Son et al. | Shear wave propagation in a layered poroelastic structure | |
| Pak et al. | Experimental synthesis of seismic horizontal free-field motion of soil in finite-domain simulations with absorbing boundary | |
| Gao et al. | Pressure-dependent seismic velocities based on effective compliance theory and an asperity deformation model | |
| CN109001804B (en) | Method, device and system for determining effective force based on three-dimensional seismic data | |
| Seyedi-Viand et al. | An alternative empirical function to predict air–water mixture bulk modulus for numerical modeling of liquefaction behavior of induced partially saturated sands | |
| CN109490988A (en) | Establish the method for being suitable for the rock physics new model of hard rock | |
| Soudkhah et al. | Wave absorbing-boundary method in seismic centrifuge simulation of vertical free-field ground motion |
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 |