CN105717021B - A kind of method for quantitatively determining crack rock anisotropy of permeability - Google Patents
A kind of method for quantitatively determining crack rock anisotropy of permeability Download PDFInfo
- Publication number
- CN105717021B CN105717021B CN201610065036.7A CN201610065036A CN105717021B CN 105717021 B CN105717021 B CN 105717021B CN 201610065036 A CN201610065036 A CN 201610065036A CN 105717021 B CN105717021 B CN 105717021B
- Authority
- CN
- China
- Prior art keywords
- permeability
- crack rock
- angle
- anisotropy
- 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.)
- Expired - Fee Related
Links
- 230000035699 permeability Effects 0.000 title claims abstract description 94
- 239000011435 rock Substances 0.000 title claims abstract description 92
- 238000000034 method Methods 0.000 title claims abstract description 30
- 239000000463 material Substances 0.000 claims abstract description 39
- 238000009795 derivation Methods 0.000 claims abstract description 5
- 238000010586 diagram Methods 0.000 claims description 9
- 239000000843 powder Substances 0.000 claims description 4
- 239000002994 raw material Substances 0.000 claims description 4
- 238000012067 mathematical method Methods 0.000 abstract description 4
- 238000004088 simulation Methods 0.000 abstract 1
- 238000012360 testing method Methods 0.000 description 7
- 238000011065 in-situ storage Methods 0.000 description 3
- 238000001764 infiltration Methods 0.000 description 3
- 230000008595 infiltration Effects 0.000 description 3
- 230000035515 penetration Effects 0.000 description 3
- 238000010146 3D printing Methods 0.000 description 2
- 238000009533 lab test Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000012466 permeate Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000011083 cement mortar Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 238000013517 stratification Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 230000009466 transformation 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, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
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)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The present invention discloses a kind of method for quantitatively determining crack rock anisotropy of permeability, includes the following steps：The actual crack rock of simulation, prepares the analog material in different angle with horizontal direction；The permeability of different angle analog material is measured, determines the functional relation at permeability and inclination angle；To the functional relation derivation, the derived function of the functional relation is obtained, which represents change rate of the permeability to inclination angle；Derived function is standardized, obtains the relational expression at unit permeation rate variable quantity and inclination angle caused by unit change of pitch angle, which is the quantitative expression of actual crack rock anisotropy of permeability to be measured.The present invention determines the permeability variation rate of crack rock different directions with Mathematical Method, and the change rate at permeability and inclination angle is standardized, the variable quantity of unit permeation rate caused by unit change of pitch angle at different angle is obtained, can accurately represent the anisotropy of permeability of the larger crack rock of permeability contrast.
Description
Technical field
The present invention relates to a kind of method of definite crack rock penetration property, more particularly to one kind quantitatively determines crack rock
The method of anisotropy of permeability.
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 of structural plane and intact rock mostly, due to discontinuous structural plane complicated in crack rock
Presence so that the penetration property of crack rock is more complicated than general homogeneous, continuous, isotropic material penetration property
It is more, significant anisotropy of permeability is shown, there is significant impact to engineering safety, therefore crack rock is permeated each to different
The research of property has important theory value and engineering significance.
Although the anisotropy of permeability of crack rock significantly affects the seepage flow property and engineering safety of rock mass, mesh
The research that the preceding anisotropy of permeability on crack rock is done is seldom, this is because crack rock has generally comprised multigroup crack
And rock mass volume is larger, limited by field condition, testing equipment and testing expenses, carry out in-situ test measurement crack at the scene
The anisotropy of permeability of rock mass is very difficult even not possible with therefore can only being measured by the method for laboratory test.But
The natural rock sample taken by laboratory test sample size usually comprises only one group or not comprising crack, can not fully reflect rock mass
Anisotropy of permeability, and crack rock has the variability of geometric properties, and different rock crannies differ greatly, and are derived from part rock
The result of the test of the natural rock sample of body does not have adequate representation, and lacks the effective crack rock that quantitatively determines at present and permeate
Anisotropic method.
A kind of Chinese invention patent application " fracture-type reservoir permeability tensor and each of Application No. 201410604481.7
Anisotropy quantitative forecasting technique " proposes a kind of carrys out table by the ratio between maxima and minima of crack rock different directions permeability
Show the method for anisotropy of permeability size, but this method is only capable of the difference that reflection crack rock different directions permeability is most worth, no
It can reflect the difference of remaining directional permeability of crack rock, the value for the anisotropy of permeability tried to achieve does not have adequate representation, because
There is an urgent need to a kind of anisotropy of permeability that effective method quantitatively determines crack rock for this.
The content of the invention
Goal of the invention：The purpose of the present invention is directed to problems of the prior art, there is provided one kind quantitatively determines crack rock
The method of body anisotropy of permeability.
Technical solution：A kind of method for quantitatively determining crack rock anisotropy of permeability of the present invention, including it is as follows
Step：
Step 1, actual crack rock is simulated, prepares the analog material in different angle with horizontal direction；
Step 2, the permeability of different angle analog material is measured, determines the functional relation at permeability and inclination angle；
Step 3, to the functional relation derivation, the derived function of the functional relation is obtained, which represents permeability
To the change rate at inclination angle；
Step 4, derived function is standardized, obtain unit permeation rate variable quantity caused by unit change of pitch angle with
The relational expression at inclination angle, the relational expression are the quantitative expression of actual crack rock anisotropy of permeability to be measured.
By the actual crack rock of equivalent material simulating, crack rock is quantitatively determined in different directions with Mathematical Method
Permeability variation, and the change rate at permeability and inclination angle is standardized, obtained expression formula reflection is different
The variable quantity of unit permeation rate caused by unit change of pitch angle on inclination direction, can accurately represent that permeability contrast is larger
The anisotropy of permeability of crack rock.
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.
Above-mentioned inclination angle is 0 °, 15 °, 30 °, 45 °, 60 °, 75 ° or 90 °.
In above-mentioned steps 2, according to the permeability of the different angle analog material measured, made in reference axis permeability with
The scatter diagram of counter-bevel, with iterative method draw with the highest functional relation of scatter diagram degree of fitting, be permeability with inclining
The functional relation at angle.
In above-mentioned steps 4, it is to the method that derived function is standardized：, should by derived function divided by original permeability
Original permeability is the permeability for the analog material that inclination angle is 0 °.Due to inclination angle be 0 ° when permeability represent rock mass vertical direction
Permeability, which can obtain unit by derived function divided by original permeability and oozes as the original permeability of rock mass
The changing value of saturating 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, be the infiltration that can determine that crack rock by indoor model test
Anisotropy, saves a large amount of material resources and financial resources；Meanwhile crack rock is quantitatively determined in different directions with Mathematical Method
Permeability variation, and the change rate at permeability and inclination angle is standardized, obtained expression formula reflection is that difference is inclined
The variable quantity of unit permeation rate caused by unit change of pitch angle on angular direction, therefore can accurately represent that permeability contrast is larger
Crack rock anisotropy of permeability；(2) analog material of the invention can be prepared by 3D printing method, can accurately be reacted
The spatial shape in each crack and distribution, gained analog material can fully simulate the property of actual crack rock, and it is fixed to use it for
During amount analysis anisotropy of permeability, experimental result is representative.
Brief description of the drawings
Fig. 1 is flow chart of the method for the present invention.
Fig. 2 is the scatter diagram for the different angle analog material and inclination angle relation made in the embodiment of the present invention.
Embodiment
Technical scheme is described further below in conjunction with the accompanying drawings.
As shown in Figure 1, a kind of method for quantitatively determining crack rock anisotropy of permeability of the present invention, including following step
Suddenly：
Step 1, actual crack rock is simulated, prepares some analog materials 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, is 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 prints the analog material with horizontal direction in different angle for raw material.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, the permeability of different angle analog material is measured, determines the functional relation at permeability and inclination angle；
According to the permeability of the different angle analog material measured, dissipating for permeability and counter-bevel is made in reference axis
Point diagram, is drawn and the highest functional relation of scatter diagram degree of fitting, as permeability and the functional relation at inclination angle with iterative method
Formula.
Step 3, to the functional relation derivation, the derived function of the functional relation is obtained, which represents permeability
To the change rate at inclination angle；
Step 4, derived function is standardized, obtain unit permeation rate variable quantity caused by unit change of pitch angle with
The relational expression at inclination angle, the relational expression are the quantitative expression of actual crack rock anisotropy of permeability to be measured.
Permeability when inclination angle is 0 ° represents the permeability of rock mass vertical direction, the infiltration using inclination angle as 0 ° of analog material
Rate is original permeability, by derived function divided by the original permeability, as the quantitative expression of crack rock anisotropy of permeability,
The expression formula represents the variable quantity of unit permeation rate caused by unit change of pitch angle on different angle direction, characterizes permeability pair
The size of the anisotropy of permeability of the sensitivity of change of pitch angle, i.e. crack rock.
The present invention simulates actual crack rock by preparing analog material, it is not necessary to carries out in-situ test, passes through indoor mould
Type experiment can determine that the anisotropy of permeability of crack rock, save a large amount of material resources and financial resources.Meanwhile with Mathematical Method
Permeability variation of the crack rock in different directions is quantitatively determined, and place is standardized to the change rate of permeability and inclination angle
Reason, obtains the variable quantity of unit permeation rate caused by unit change of pitch angle on different angle direction, can accurately represent to permeate
The anisotropy of permeability for the crack rock that rate differs greatly.
Embodiment
To quantitatively determine at some hydropower station riverbed dam foundation exemplified by the anisotropy of permeability of crack rock, the method for the present invention includes
Following steps：
(1) split according to the spatial shape in crack, quantity and distribution situation, making in the volume of actual crack rock and rock mass
The three-dimensional digital model of gap rock mass, 3D printer is inputted by three-dimensional digital model, rotates the three-dimensional digital model of crack rock, with
The powder of actual crack rock prints respectively for raw material to incline from the horizontal by 0 °, 15 °, 30 °, 45 °, 60 °, 75 ° or 90 °
The actual crack rock analog material at angle.
(2) using the permeability of measuring gas permebility of rock device measurement different angle analog material, according to what is measured
The permeability of different angle analog material, using inclination angle as x-axis, permeability is y-axis, makes the infiltration of the crack rock analog material
The scatter diagram of rate and inclination angle relation, such as Fig. 2, is calculated and the highest functional relation y=of the scatter diagram degree of fitting with iterative method
3.84x0.43, characterize crack rock permeability and the relation at inclination angle.
(3) derivation, gained derived function y=1.65x are carried out to functional relation obtained by step (2)-0.57As surveyed crack
Expression formula of the rock mass permeability to the change rate at inclination angle.
(5) permeability is standardized the expression formula of the change rate at inclination angle, using inclination angle as 0 ° of similar material
The permeability k of material0For original permeability, k0=1.9411, derived function y=1.65x-0.57Divided by 1.9411 obtain standardization
Expression formula y=0.85x of the crack rock permeability afterwards to the change rate at inclination angle-0.57, as crack rock anisotropy of permeability
Quantitative expression, which represents the variable quantity of unit permeation rate caused by unit change of pitch angle on different angle direction,
Characterize sensitivity of the permeability to change of pitch angle, the i.e. size of the anisotropy of permeability of crack rock.
Claims (3)
- A kind of 1. method for quantitatively determining crack rock anisotropy of permeability, it is characterised in that include the following steps：Step 1, actual crack rock is simulated, prepares the analog material in different angle with horizontal direction；Step 2, the permeability of different angle analog material is measured, determines the functional relation at permeability and inclination angle；According to measuring Different angle analog material permeability, the scatter diagram of permeability and counter-bevel is made in reference axis, is obtained with iterative method Go out and the highest functional relation of scatter diagram degree of fitting, the as functional relation at permeability and inclination angle；Step 3, to the functional relation derivation, the derived function of the functional relation is obtained, which represents permeability to inclining The change rate at angle；Step 4, the derived function is standardized, obtain unit permeation rate variable quantity caused by unit change of pitch angle with The relational expression at inclination angle, the relational expression are the quantitative expression of actual crack rock anisotropy of permeability to be measured；Wherein, it is described right The method that derived function is standardized is：By the derived function divided by original permeability, which is that inclination angle is The permeability of 0 ° of analog material.
- 2. the method according to claim 1 for quantitatively determining crack rock anisotropy of permeability, it is characterised in that step 1 In, 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. the method according to claim 1 or 2 for quantitatively determining crack rock anisotropy of permeability, it is characterised in that institute Inclination angle is stated as 0 °, 15 °, 30 °, 45 °, 60 °, 75 ° or 90 °.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610065036.7A CN105717021B (en) | 2016-01-29 | 2016-01-29 | A kind of method for quantitatively determining crack rock anisotropy of permeability |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610065036.7A CN105717021B (en) | 2016-01-29 | 2016-01-29 | A kind of method for quantitatively determining crack rock anisotropy of permeability |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105717021A CN105717021A (en) | 2016-06-29 |
| CN105717021B true CN105717021B (en) | 2018-05-11 |
Family
ID=56154384
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610065036.7A Expired - Fee Related CN105717021B (en) | 2016-01-29 | 2016-01-29 | A kind of method for quantitatively determining crack rock anisotropy of permeability |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105717021B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107806339B (en) * | 2016-09-08 | 2020-10-23 | 中国石油化工股份有限公司 | Fracturing fracture conductivity experiment method |
| CN111398116B (en) * | 2020-03-04 | 2020-12-29 | 中国地质大学（武汉） | Method for representing shale anisotropy by mercury porosimetry in specific direction |
| CN112906312A (en) * | 2021-01-19 | 2021-06-04 | 中国矿业大学(北京) | Joint stiffness-variable fracture coal rock fluid-solid coupling simulation parameter selection method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101625354A (en) * | 2009-08-05 | 2010-01-13 | 山东大学 | Similar material for fluid-solid coupling model test |
| CN103454127A (en) * | 2013-07-29 | 2013-12-18 | 山东科技大学 | Similar material for small and medium solid-fluid coupled model test and preparation method thereof |
| CN104500050A (en) * | 2014-10-31 | 2015-04-08 | 中国石油大学(华东) | Fractured reservoir permeability tensor and anisotropy quantitative prediction method |
| CN104515696A (en) * | 2014-12-09 | 2015-04-15 | 河海大学 | Method for preparation of columnar jointed rock mass similar material sample by 3D printing technology |
-
2016
- 2016-01-29 CN CN201610065036.7A patent/CN105717021B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101625354A (en) * | 2009-08-05 | 2010-01-13 | 山东大学 | Similar material for fluid-solid coupling model test |
| CN103454127A (en) * | 2013-07-29 | 2013-12-18 | 山东科技大学 | Similar material for small and medium solid-fluid coupled model test and preparation method thereof |
| CN104500050A (en) * | 2014-10-31 | 2015-04-08 | 中国石油大学(华东) | Fractured reservoir permeability tensor and anisotropy quantitative prediction method |
| CN104515696A (en) * | 2014-12-09 | 2015-04-15 | 河海大学 | Method for preparation of columnar jointed rock mass similar material sample by 3D printing technology |
Non-Patent Citations (1)
| Title |
|---|
| 采动裂隙岩体应力恢复及其渗透性演化;王文学;《中国博士学位论文全文数据库 工程科技Ⅰ辑》;20150415(第4期);第100-101、114、126-131页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105717021A (en) | 2016-06-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105510207B (en) | Determine method of the crack rock different directions permeability to porosity sensitivity | |
| CN104374619B (en) | A kind of preparation method of irregular prismatical joint Fracture Networks model core sample | |
| CN103954511B (en) | A kind of Fracture Networks rock shearing-seepage flow coupling experiment method | |
| CN102289980B (en) | Method for preparing preset cracks in hydrofracture simulation specimen | |
| CN103674658B (en) | A kind of preparation method of random crack test model | |
| CN105606513B (en) | Determine the method for crack rock permeability and porosity to effective stress sensitivity | |
| CN104833555B (en) | The preparation method of interrupted filling fracture rock sample model | |
| Xia et al. | Parametric study of smooth joint parameters on the mechanical behavior of transversely isotropic rocks and research on calibration method | |
| CN105717021B (en) | A kind of method for quantitatively determining crack rock anisotropy of permeability | |
| CN105547831A (en) | Model test method for measuring deformation property of rock mass containing complex structural faces | |
| CN111099856B (en) | Similar simulation experiment material for dam body of underground reservoir and determination method for proportion of similar simulation experiment material | |
| CN105223080B (en) | Evaluation method for jointed rock performance and bolting and grouting effects under compression-shear condition | |
| CN106769753B (en) | The trace structure column joint samples preparation method containing chisel for seepage tests | |
| CN105486555A (en) | Method for preparing columnar jointed rock mass model containing blind joints and interlayer shear belts | |
| Niu et al. | Experimental study of split grouting reinforcement mechanism in filling medium and effect evaluation | |
| CN105241802B (en) | Prismatical joint rock mass permeability tensor method of testing and anisotropy quantitative evaluating method | |
| Wang et al. | Experimental and numerical study on failure modes and shear strength parameters of rock-like specimens containing two infilled flaws | |
| CN105388047A (en) | Preparation method of columnar joint fissure network model rock core samples | |
| CN104537674B (en) | A kind of detection method of epoxy asphalt concrete aggregate grading | |
| Tong et al. | DEM simulation of biaxial compression experiments of inherently anisotropic granular materials and the boundary effects | |
| CN208398855U (en) | Shield tunnel construction stress deformation indoor model test device | |
| QI et al. | Fracturing mechanism of rock-like specimens with different joint densities based on DIC technology | |
| CN114626215A (en) | Rock slope modeling method based on discrete element and synthetic rock mass technology | |
| CN105651964B (en) | A kind of method for determining Investigation of Representative Elementary Volume For Fractured Rock Mass product | |
| CN102915665A (en) | Coalbed methane earthquake physical model and making method thereof |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180511 |