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 PDF

Info

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
Application number
CN201610064275.0A
Other languages
Chinese (zh)
Other versions
CN105510207A (en
Inventor
巢志明
王环玲
徐卫亚
贾朝军
闫龙
杨兰兰
赵恺
夏季
滕志强
官明开
李栋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hohai University HHU
Original Assignee
Hohai University HHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hohai University HHU filed Critical Hohai University HHU
Priority to CN201610064275.0A priority Critical patent/CN105510207B/en
Publication of CN105510207A publication Critical patent/CN105510207A/en
Application granted granted Critical
Publication of CN105510207B publication Critical patent/CN105510207B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, 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 counter-bevel direction analog material respectively;Each derived function is standardized, obtains the variable quantity of unit permeation rate caused by the unit porosity change of counter-bevel direction, the variable quantity is the sensitivity that the permeability in actual crack rock counter-bevel direction to be measured changes to porosity.

Description

Determine method of the crack rock different directions permeability to porosity sensitivity
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 in-situ 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 counter-bevel 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 counter-bevel direction The variable quantity of position permeability, the variable quantity are that the permeability in actual crack rock counter-bevel 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 three-dimensional digital model of rock mass;
(2) three-dimensional digital model is inputted into 3D printer, the three-dimensional 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 above-mentioned 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 above-mentioned steps 2, apply 3MPa, 8MPa, 15Mpa, 25MPa, 30MPa and 35MPa confined pressure step by step.
In above-mentioned 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 counter-bevel 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 in-situ 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 three-dimensional digital model of crack rock;The three-dimensional digital model is inputted into 3D printer, the three-dimensional 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 counter-bevel 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 counter-bevel direction The variable quantity of position permeability, the variable quantity are that the permeability in actual crack rock counter-bevel 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 in-situ 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 three-dimensional digital model of the crack rock, and three-dimensional 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 counter-bevel 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 counter-bevel direction unit porosity caused by unit ooze The variable quantity of saturating rate, the variable quantity be the permeability in actual crack rock counter-bevel 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 counter-bevel 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 Three-dimensional digital model;
(2) by the three-dimensional digital model input 3D printer, rotate the three-dimensional 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.
CN201610064275.0A 2016-01-29 2016-01-29 Determine method of the crack rock different directions permeability to porosity sensitivity Active CN105510207B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610064275.0A CN105510207B (en) 2016-01-29 2016-01-29 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) 2016-01-29 2016-01-29 Determine method of the crack rock different directions permeability to porosity sensitivity

Publications (2)

Publication Number Publication Date
CN105510207A CN105510207A (en) 2016-04-20
CN105510207B true CN105510207B (en) 2018-02-06

Family

ID=55718348

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610064275.0A Active CN105510207B (en) 2016-01-29 2016-01-29 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)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107917865B (en) * 2016-10-11 2020-01-31 中国石油化工股份有限公司 compact sandstone reservoir multi-parameter permeability prediction method
CN108875096B (en) * 2017-05-09 2021-11-23 中国石油化工股份有限公司 Method for constructing simulation model of core scale
CN109752302B (en) * 2017-11-02 2021-09-28 中国石油天然气股份有限公司 Method for predicting porosity of tight sandstone reservoir in upward inclination direction
CN108760598B (en) * 2018-04-12 2019-12-10 河海大学 Method for determining minimum cyclic loading and unloading times required by filling joint rock compaction
CN108681793B (en) * 2018-05-18 2022-02-01 中国石油天然气集团有限公司 Deep reservoir oil extraction index prediction method and device
CN109063334B (en) * 2018-08-03 2020-10-02 河海大学 Fluid-solid coupling model construction method of low-permeability porous medium material
CN109374497B (en) * 2018-10-29 2020-12-11 河海大学 Rock micro-pore structure testing method

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101248531B1 (en) * 2011-12-05 2013-04-03 한국지질자원연구원 Apparatus and method for measuring porosity and permeability of dioxide carbon underground storage medium
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
CN103528934B (en) * 2013-10-29 2016-05-25 中国石油大学(北京) Measure the cross-correlation method of Oil in Super-low Permeability rock permeability stress sensitivity
CN104515696A (en) * 2014-12-09 2015-04-15 河海大学 Method for preparation of columnar jointed rock mass similar material sample by 3D printing technology
CN105241802B (en) * 2015-11-03 2018-01-19 河海大学 Prismatical joint rock mass permeability tensor method of testing and anisotropy quantitative evaluating method

Also Published As

Publication number Publication date
CN105510207A (en) 2016-04-20

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 gas-water 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) Full-automatic 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 non-Darcy 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 Thick-Walled 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