CN105510207B  Determine method of the crack rock different directions permeability to porosity sensitivity  Google Patents
Determine method of the crack rock different directions permeability to porosity sensitivity Download PDFInfo
 Publication number
 CN105510207B CN105510207B CN201610064275.0A CN201610064275A CN105510207B CN 105510207 B CN105510207 B CN 105510207B CN 201610064275 A CN201610064275 A CN 201610064275A CN 105510207 B CN105510207 B CN 105510207B
 Authority
 CN
 China
 Prior art keywords
 porosity
 permeability
 analog material
 crack rock
 inclination angle
 Prior art date
 Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
 Active
Links
 230000035699 permeability Effects 0.000 title claims abstract description 96
 239000011435 rock Substances 0.000 title claims abstract description 89
 230000035945 sensitivity Effects 0.000 title claims abstract description 25
 239000000463 material Substances 0.000 claims abstract description 79
 238000009795 derivation Methods 0.000 claims abstract description 5
 238000010586 diagram Methods 0.000 claims description 15
 239000000843 powder Substances 0.000 claims description 4
 239000002994 raw material Substances 0.000 claims description 4
 230000000875 corresponding Effects 0.000 claims description 3
 238000009738 saturating Methods 0.000 claims description 2
 239000011148 porous material Substances 0.000 description 5
 230000035515 penetration Effects 0.000 description 4
 238000010192 crystallographic characterization Methods 0.000 description 3
 238000011065 insitu storage Methods 0.000 description 3
 238000001764 infiltration Methods 0.000 description 3
 238000010146 3D printing Methods 0.000 description 2
 230000000694 effects Effects 0.000 description 2
 238000009533 lab test Methods 0.000 description 2
 238000005259 measurement Methods 0.000 description 2
 238000004445 quantitative analysis Methods 0.000 description 2
 238000004458 analytical method Methods 0.000 description 1
 239000011083 cement mortar Substances 0.000 description 1
 230000018109 developmental process Effects 0.000 description 1
 238000005516 engineering process Methods 0.000 description 1
 238000004321 preservation Methods 0.000 description 1
 238000003325 tomography Methods 0.000 description 1
 230000001131 transforming Effects 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, porevolume, or surfacearea of porous materials
 G01N15/08—Investigating permeability, porevolume, or surface area of porous materials
Abstract
The present invention discloses a kind of determination crack rock different directions permeability to the method for porosity sensitivity, comprises the following steps：Actual crack rock is simulated, prepares the analog material in different angle with horizontal direction；Actual crack rock is simulated, prepares the analog material in different angle with horizontal direction；Apply confined pressure step by step to each inclination angle analog material, measure the permeability and porosity of each inclination angle analog material after every grade of confined pressure steady respectively, determine the permeability of each inclination direction analog material and the functional relation of porosity；To each functional relation derivation, gained derived function represents rate of change of the permeability to porosity of counterbevel direction analog material respectively；Each derived function is standardized, obtains the variable quantity of unit permeation rate caused by the unit porosity change of counterbevel direction, the variable quantity is the sensitivity that the permeability in actual crack rock counterbevel direction to be measured changes to porosity.
Description
Technical field
The present invention relates to a kind of method for determining crack rock penetration property, and in particular to one kind determines that crack rock is different
Method of the directional permeability to porosity sensitivity.
Background technology
It experienced very long geology into raw transformation effect and table life to develop, rock mass of the preservation in certain geology ambient stress
Internal all extensive developments the structural planes such as joint, stratification, contact band, shear band, tomography, therefore the rock run into Practical Project
The crack rock that body is all made up of structural plane and intact rock mostly, due to discontinuous structural plane complicated in crack rock
There is control action to the penetration property of rock mass so that the penetration property of crack rock is than in general homogeneous, continuous, isotropism material
Complicated more of the penetration property of material, have significant impact to engineering safety, therefore crack rock pore structure is permeated to it
Property Quality Research has important theory value and engineering significance.
The research done at present on influence of the pore structure to the seepage flow property of rock mass of crack rock is seldom, this be by
Multigroup crack has been generally comprised in crack rock and rock mass volume is larger, has been limited by field condition, testing equipment and testing expenses
System, at the scene carry out insitu test measurement crack rock permeability and porosity be it is very difficult even not possible with, because
This can only be measured by the method for laboratory test.But the natural rock sample taken by laboratory test sample size generally comprises only one group
Or not comprising crack, it can not fully reflect the influence of the pore structure of crack rock to the seepage flow property of rock mass, and crack rock
Body has the polytropy of geometric properties, and different rock crannies differ greatly, and is derived from the result of the test of the natural rock sample of part rock mass
Without adequate representation, and lack effective seepage flow property of the pore structure to rock mass for quantitatively determining crack rock at present
Influence method, but the pore structure of crack rock significantly affects on the permeability of rock mass, and a small amount of change of porosity is
Being widely varied for permeability can be caused, therefore there is an urgent need to one kind to quantitatively determine crack rock different directions permeability to hole
Spend the method for sensitivity
The content of the invention
Goal of the invention：The purpose of the present invention is to be directed to problems of the prior art, there is provided one kind determines crack rock
Method of the different directions permeability to porosity sensitivity.
Technical scheme：Side of the crack rock different directions permeability to porosity sensitivity is determined the present invention relates to a kind of
Method, comprise the following steps：
Step 1, actual crack rock is simulated, prepares the analog material in different angle with horizontal direction；
Step 2, apply confined pressure step by step to each inclination angle analog material, measure the similar material in each inclination angle after every grade of confined pressure steady respectively
The permeability and porosity of material, determine the permeability of each inclination direction analog material and the functional relation of porosity；
Step 3, the infiltration of counterbevel direction analog material is represented respectively to each functional relation derivation, gained derived function
Rate of change of the rate to porosity；
Step 4, each derived function is standardized, obtained single caused by the unit porosity change of counterbevel direction
The variable quantity of position permeability, the variable quantity are that the permeability in actual crack rock counterbevel direction to be measured changes to porosity
Sensitivity.
The rate of change of permeability on different angle direction to porosity is obtained by the method for the present invention, after standardization
Obtain on each inclination direction unit permeation rate variable quantity caused by the change of unit porosity, can accurate characterization permeability contrast it is larger
Anisotropy of permeability rock mass different directions permeability to the sensitivity of porosity.
Specifically, in step 1, prepare includes with horizontal direction in the method for the analog material of different angle：
(1) according to the form in crack and distribution in the geometric shape of actual crack rock and the rock mass, the crack is made
The threedimensional digital model of rock mass；
(2) threedimensional digital model is inputted into 3D printer, the threedimensional digital model is rotated, with the powder of actual crack rock
The analog material with horizontal direction in different angle is printed for raw material in end.
Preferably, inclination angle is 0 °, 15 °, 30 °, 45 °, 60 °, 75 ° or 90 °.
Preferably, in step 2, any inclination angle analog material is chosen, with the hole of the inclination angle analog material under every grade of confined pressure
Porosity represents porosity of all inclination angle analog materials under this grade of confined pressure.Due to different angle direction analog material porosity
Size and with confined pressure changing rule close to consistent, therefore the porosity of measurable a certain inclination angle analog material is all to represent
The porosity of analog material on inclination direction.
In abovementioned steps 2, according to permeability and porosity of each inclination angle analog material measured under every grade of confined pressure, with hole
Porosity is abscissa, permeability is ordinate, makes each inclination direction porosity and the scatter diagram of permeability, passes through iterative method point
Do not draw and each scatter diagram degree of fitting highest functional relation, the permeability and hole of as corresponding inclination direction analog material
The functional relation of degree.
In abovementioned steps 2, apply 3MPa, 8MPa, 15Mpa, 25MPa, 30MPa and 35MPa confined pressure step by step.
In abovementioned steps 4, the method being standardized to each derived function is：By the derived function obtained by step 3 divided by
The original permeability of counterbevel direction analog material, the original permeability are the infiltration of the analog material measured under initial confined pressure
Rate.
Beneficial effect：Compared with prior art, remarkable advantage of the invention is：(1) present invention is by preparing similar material
Expect to simulate actual crack rock, it is not necessary to carry out insitu test, can determine that anisotropy of permeability splits by indoor model test
Gap rock mass different directions permeability saves a large amount of material resources and financial resources to the sensitivity of porosity；Meanwhile with mathematical analysis side
Standard measure determines rate of change of the crack rock different directions permeability to porosity, and the rate of change expression formula is standardized
Processing, the variable quantity of unit permeation rate caused by the change of unit porosity on different angle direction is obtained, can be oozed with accurate characterization
Sensitivity of the crack rock different directions permeability that saturating rate differs greatly to porosity；(2) analog material of the invention can
Prepared by 3D printing method, can accurately react spatial shape and the distribution in each crack, gained analog material being capable of abundant mould
Intend the property of actual crack rock, use it for sensitive journey of the permeability to porosity of quantitative analysis anisotropy of permeability rock mass
When spending, experimental result is representative.
Brief description of the drawings
Fig. 1 is the porosity of analog material and the scatter diagram of permeability that inclination angle is 0 °；
Fig. 2 is the porosity of analog material and the scatter diagram of permeability that inclination angle is 15 °；
Fig. 3 is the porosity of analog material and the scatter diagram of permeability that inclination angle is 30 °；
Fig. 4 is the porosity of analog material and the scatter diagram of permeability that inclination angle is 45 °；
Fig. 5 is the porosity of analog material and the scatter diagram of permeability that inclination angle is 60 °；
Fig. 6 is the porosity of analog material and the scatter diagram of permeability that inclination angle is 75 °；
Fig. 7 is the porosity of analog material and the scatter diagram of permeability that inclination angle is 90 °.
Embodiment
Technical scheme is described further below in conjunction with the accompanying drawings.
A kind of determination crack rock different directions permeability of the present invention is including as follows to the method for porosity sensitivity
Step：
Step 1, actual crack rock is simulated, prepares the analog material in different angle with horizontal direction；
Inclination angle can be 0 °, 15 °, 30 °, 45 °, 60 °, 75 ° or 90 °；Specifically, prepare with horizontal direction in different angle
The method of analog material includes the following two kinds：
The first, according to the form in crack and distribution in the geometric shape of actual crack rock and the rock mass, makes respectively
Act on the mould for preparing different angle analog material；Cement mortar is configured according to the mechanical property of actual crack rock, poured into
Mould, maintenance, the demoulding, obtains the analog material of different angle.
Second, according to the form in crack and distribution in the geometric shape of actual crack rock and the rock mass, making should
The threedimensional digital model of crack rock；The threedimensional digital model is inputted into 3D printer, the threedimensional digital model is rotated, with reality
The powder of crack rock is that raw material prints the analog material in different angle with horizontal direction.Phase is prepared by 3D printing method
Like material, spatial shape and the distribution in each crack can be accurately reacted, gained analog material can fully simulate actual crack rock
The property of body, when using it for quantitative analysis anisotropy of permeability, experimental result is representative.
Step 2, apply confined pressure step by step to each inclination angle analog material, measure the similar material in each inclination angle after every grade of confined pressure steady respectively
The permeability and porosity of material, determine the permeability of each inclination direction analog material and the functional relation of porosity；
Size due to different angle direction analog material porosity and the changing rule with confined pressure are close consistent, therefore
Porosity of a certain inclination angle analog material under every grade of confined pressure can be only measured, it can represent all inclination angle analog materials and be enclosed in the level
The porosity of pressure.
3MPa, 8MPa, 15Mpa, 25MPa, 30MPa and 35MPa confined pressure can be applied step by step.
According to permeability and porosity of each inclination angle analog material measured under every grade of confined pressure, using porosity as horizontal seat
Mark, permeability are ordinate, make each inclination direction porosity and the scatter diagram of permeability respectively in reference axis, pass through iteration
Method draw respectively with each scatter diagram degree of fitting highest functional relation, the permeability of as corresponding inclination direction analog material with
The functional relation of porosity.
Step 3, the infiltration of counterbevel direction analog material is represented respectively to each functional relation derivation, gained derived function
Rate of change of the rate to porosity；
Step 4, each derived function is standardized, obtained single caused by the unit porosity change of counterbevel direction
The variable quantity of position permeability, the variable quantity are that the permeability in actual crack rock counterbevel direction to be measured changes to porosity
Sensitivity.
Original permeability using the permeability of each inclination direction analog material that is measured under initial confined pressure as the direction, by certain
Derived function on one inclination direction divided by original permeability in this direction, obtain caused by unit porosity change in this direction
The variable quantity of unit permeation rate, characterize the sensitivity that permeability changes to porosity in this direction.
Actual crack rock is simulated by preparing analog material, it is not necessary to carry out insitu test, pass through indoor model test
Can determine that sensitivity of the anisotropy of permeability crack rock different directions permeability to porosity, save a large amount of material resources and
Financial resources；Meanwhile rate of change of the crack rock different directions permeability to porosity is quantitatively determined with Mathematical Method, and to this
Rate of change expression formula is standardized, and obtains unit permeation rate caused by the change of unit porosity on different angle direction
Variable quantity, sensitivity that can be with the larger crack rock different directions permeability of accurate characterization permeability contrast to porosity.
Embodiment
To quantitatively determine at some hydropower station riverbed dam foundation anisotropy of permeability rock mass different directions permeability to stress sensitive
Exemplified by degree, determination crack rock permeability of the invention comprises the following steps to the method for porosity sensitivity：
(1) according to the spatial shape of structural plane in the volume of actual anisotropy of permeability crack rock and rock mass, quantity and
Distribution situation, makes the threedimensional digital model of the crack rock, and threedimensional digital model is inputted into 3D printer, rotates three dimension
Word model, using the powder of actual crack rock to be measured as raw material, print respectively from the horizontal by 0 °, 15 °, 30 °, 45 °,
The crack rock analog material at 60 °, 75 ° and 90 ° inclination angles.
(2) 3MPa, 8MPa, 15Mpa, 25MPa, 30MPa and 35MPa are applied step by step respectively to different angle analog material
Confined pressure, utilize the permeability of different angle analog material under measuring gas permebility of rock measurement device difference confined pressure；Due to not
Size with inclination direction analog material porosity and the changing rule with confined pressure are close consistent, therefore with the similar material in 0 ° of inclination angle
Expect that the porosity under different confined pressures represents the porosity of all inclination angle analog materials；Using porosity φ as abscissa, permeability k
For ordinate, the permeability of different angle direction analog material and the scatter diagram of porosity are made respectively, such as Fig. 1~Fig. 7；Then
Calculated respectively with iterative method with each inclination direction scatter diagram degree of fitting highest functional relation, as crack rock not Tongfang
To permeability and the functional relation of porosity, such as table 1.
The functional relation of the crack rock different directions permeability of table 1 and porosity
Inclination angle/°  The functional relation of permeability and porosity 
0  K=11.18 φ^{0}.^{69} 
15  K=130.71 φ^{0}.^{82} 
30  K=147.73 φ^{0}.^{57} 
45  K=131.19 φ^{0}.^{81} 
60  K=141.91 φ^{0}.^{66} 
75  K=117.51 φ^{0}.^{62} 
90  K=129.56 φ^{0}.^{7} 
(3) functional relation of the permeability to different angle direction analog material and porosity carries out derivation, and gained is led
Expression formula of the absolute value of function as crack rock different directions permeability to the rate of change of porosity, such as table 2.
Expression formula of the crack rock different directions permeability of table 2 to the rate of change of porosity
(4) crack rock different directions permeability is standardized to the expression formula of the rate of change of porosity, just
Original permeability during beginning confined pressure 3MPa at 0 °, 15 °, 30 °, 45 °, 60 °, 75 ° and 90 ° inclination direction is respectively 5.11e^{17}、
5.07×10^{16}、7.38×10^{16}、5.23×10^{16}、2.05×10^{15}、1.33×10^{16}、5.87×10^{16}, with correspondence direction
The permeability of crack rock takes absolute value to the expression formula divided by the original permeability of the rate of change of porosity, obtains crack
Quantitative expression of the anisotropy of permeability rock mass permeability to porosity sensitivity.
Such as table 3, the expression formula represents the variable quantity of unit permeation rate caused by the change of correspondence direction unit porosity respectively,
The value characterizes the sensitivity that permeability changes to porosity on correspondence direction.
Quantitative expression of the crack rock different directions permeability of table 3 to porosity sensitivity
Claims (4)
1. a kind of determine method of the crack rock different directions permeability to porosity sensitivity, it is characterised in that including such as
Lower step：
Step 1, actual crack rock is simulated, prepares the analog material in different angle with horizontal direction；
Step 2, apply confined pressure step by step to each inclination angle analog material, measure each inclination angle analog material after every grade of confined pressure steady respectively
Permeability and porosity, determine the permeability of each inclination direction analog material and the functional relation of porosity；Wherein, choose and appoint
One inclination angle analog material, all inclination angle analog materials are represented in the level with the porosity of the inclination angle analog material under every grade of confined pressure
Porosity under confined pressure；According to permeability and porosity of each inclination angle analog material measured under every grade of confined pressure, with porosity
It is ordinate for abscissa, permeability, makes each inclination direction porosity and the scatter diagram of permeability, by iterative method respectively
Go out with each scatter diagram degree of fitting highest functional relation, the permeability of as corresponding inclination direction analog material and porosity
Functional relation；
Step 3, the permeability pair of counterbevel direction analog material is represented respectively to each functional relation derivation, gained derived function
The rate of change of porosity；
Step 4, each derived function is standardized, obtain the change of counterbevel direction unit porosity caused by unit ooze
The variable quantity of saturating rate, the variable quantity be the permeability in actual crack rock counterbevel direction to be measured porosity is changed it is quick
Sense degree；The method that is standardized to each derived function is：By the derived function obtained by step 3 divided by counterbevel side
To the original permeability of analog material, the original permeability is the permeability of the analog material measured under initial confined pressure.
2. according to claim 1 determine method of the crack rock different directions permeability to porosity sensitivity, its
It is characterised by, in step 1, the preparation includes with horizontal direction in the method for the analog material of different angle：
(1) according to the form in crack and distribution in the geometric shape of actual crack rock and the rock mass, the crack rock is made
Threedimensional digital model；
(2) by the threedimensional digital model input 3D printer, rotate the threedimensional digital model, using the powder of actual crack rock as
It is in the analog material of different angle that raw material, which is printed with horizontal direction,.
3. according to claim 1 or 2 determine method of the crack rock different directions permeability to porosity sensitivity,
Characterized in that, the inclination angle is 0 °, 15 °, 30 °, 45 °, 60 °, 75 ° or 90 °.
4. according to claim 1 determine method of the crack rock different directions permeability to porosity sensitivity, its
It is characterised by, in step 2, applies 3MPa, 8MPa, 15Mpa, 25MPa, 30MPa and 35MPa confined pressure step by step.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201610064275.0A CN105510207B (en)  20160129  20160129  Determine method of the crack rock different directions permeability to porosity sensitivity 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201610064275.0A CN105510207B (en)  20160129  20160129  Determine method of the crack rock different directions permeability to porosity sensitivity 
Publications (2)
Publication Number  Publication Date 

CN105510207A CN105510207A (en)  20160420 
CN105510207B true CN105510207B (en)  20180206 
Family
ID=55718348
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201610064275.0A Active CN105510207B (en)  20160129  20160129  Determine method of the crack rock different directions permeability to porosity sensitivity 
Country Status (1)
Country  Link 

CN (1)  CN105510207B (en) 
Families Citing this family (7)
Publication number  Priority date  Publication date  Assignee  Title 

CN107917865B (en) *  20161011  20200131  中国石油化工股份有限公司  compact sandstone reservoir multiparameter permeability prediction method 
CN108875096B (en) *  20170509  20211123  中国石油化工股份有限公司  Method for constructing simulation model of core scale 
CN109752302B (en) *  20171102  20210928  中国石油天然气股份有限公司  Method for predicting porosity of tight sandstone reservoir in upward inclination direction 
CN108760598B (en) *  20180412  20191210  河海大学  Method for determining minimum cyclic loading and unloading times required by filling joint rock compaction 
CN108681793B (en) *  20180518  20220201  中国石油天然气集团有限公司  Deep reservoir oil extraction index prediction method and device 
CN109063334B (en) *  20180803  20201002  河海大学  Fluidsolid coupling model construction method of lowpermeability porous medium material 
CN109374497B (en) *  20181029  20201211  河海大学  Rock micropore structure testing method 
Family Cites Families (5)
Publication number  Priority date  Publication date  Assignee  Title 

KR101248531B1 (en) *  20111205  20130403  한국지질자원연구원  Apparatus and method for measuring porosity and permeability of dioxide carbon underground storage medium 
CN103207137A (en) *  20130408  20130717  河海大学  Device for fully automatically measuring penetration rate and void ratio of compact rock under dynamic confining pressure 
CN103528934B (en) *  20131029  20160525  中国石油大学(北京)  Measure the crosscorrelation method of Oil in Superlow Permeability rock permeability stress sensitivity 
CN104515696A (en) *  20141209  20150415  河海大学  Method for preparation of columnar jointed rock mass similar material sample by 3D printing technology 
CN105241802B (en) *  20151103  20180119  河海大学  Prismatical joint rock mass permeability tensor method of testing and anisotropy quantitative evaluating method 

2016
 20160129 CN CN201610064275.0A patent/CN105510207B/en active Active
Also Published As
Publication number  Publication date 

CN105510207A (en)  20160420 
Similar Documents
Publication  Publication Date  Title 

CN105510207B (en)  Determine method of the crack rock different directions permeability to porosity sensitivity  
CN105606513B (en)  Determine the method for crack rock permeability and porosity to effective stress sensitivity  
US10712252B2 (en)  Reservoir depletion/injection simulation under true triaxial stress conditions  
CN104568694B (en)  Method for testing gaswater relative permeability of dense core  
CN103575631B (en)  Rock permeability test macro and method of testing  
CN103760320B (en)  The method of testing of swelled ground water cut and bulging force relation under tunnel support condition  
CN105675469A (en)  Fullautomatic test system and measurement method for gas permeability of rock  
CN105241802B (en)  Prismatical joint rock mass permeability tensor method of testing and anisotropy quantitative evaluating method  
CN107622139B (en)  Calculation method of crack permeability  
CN203672756U (en)  Indoor constant head test device for anisotropic soil bidirectional permeability coefficient  
CN107462936B (en)  Utilize the method for pressure monitoring Data Inversion low permeability reservoir nonDarcy percolation law  
Yin et al.  Investigating the roles of included angle and loading condition on the critical hydraulic gradient of real rock fracture networks  
CN105424466A (en)  Method for evaluating resistance of sandy soil on embedded submarine pipeline under action of waves  
Qian et al.  Water seepage flow in concrete  
CN110161216A (en)  A kind of device and method measuring rock core irreducible water saturation  
Zhang et al.  A study on effect of seepage direction on permeability stress test  
CN105606261A (en)  Device and test method for testing three dimensional stress state in concrete  
CN105717021B (en)  A kind of method for quantitatively determining crack rock anisotropy of permeability  
CN109406362B (en)  Method for determining relative permeability of gas and water  
CN106841002B (en)  Determine the method that the blind joint density threshold of slippage effect is generated in rock test  
CN109083630B (en)  Method for evaluating plugging performance of drilling fluid  
CN108444886B (en)  A kind of method of determining rock relative permeability  
CN107907464B (en)  Device and method for measuring performance of permeable stone cement slurry for fracturing  
Santana et al.  Workflow on Incorporating ThickWalled Cylinder Test Results in Finite Element Models of near Wellbore for Sanding Prediction Studies  
CN105651964B (en)  A kind of method for determining Investigation of Representative Elementary Volume For Fractured Rock Mass product 
Legal Events
Date  Code  Title  Description 

C06  Publication  
PB01  Publication  
C10  Entry into substantive examination  
SE01  Entry into force of request for substantive examination  
GR01  Patent grant  
GR01  Patent grant 