The permeability of intact marble samples collected from the depth of 1.6 km in southwestern China is investigated under moderate confining pressures and temperatures. No microcracks initiate or propagate during the tests, and the variation of permeability is due to the change of aperture of microcracks. Test results show a considerable decrease of permeability along with confining pressure increase from 10 to 30 MPa and temperature increase from 15 to 40°C. The thermal effect on the permeability is notable in comparison with the influence of the stress. A simple permeability evolution law is developed to correlate the permeability and the porosity in the compressive regime based on the microphysical geometric linkage model. Using this law, the permeability in the compressive regime for crystalline rock can be predicted from the volumetric strain curve of mechanical tests.
Chinese Fundamental Research (973) Program2013CB036001. Introduction
The transport properties of rocks are of significant importance for many topics of earth sciences, for example, exploration and production of hydrocarbons, oil/gas storage in deep cavern, and high-level radioactive nuclear waste geological disposal [1–4]; they are also vital to understand the fundamental geologic processes such as heat and mass transfer, earthquake, and metamorphism [5–9]. Permeability is a key parameter to describe the transport properties of materials and it can be measured using different methods, for example, steady method and pulse-test [10]. Due to Klinkenberg effect, the measured apparent gas permeability for tight rock is often different from the intrinsic permeability, which is related solely to the pore geometry (e.g., porosity, pore shape, and pore-size distribution) of the rock and is independent of the properties of the fluid. In the following texts, the term permeability is used to represent intrinsic permeability for brevity.
Several factors, for example, stress level, chemical effects, and temperature, can influence rock permeability by changing the pore geometry of the material [11–17]. These factors can significantly change the permeability of rock sample. For instance, Souley et al. [18], Schulze et al. [19], Jiang et al. [20], and Chen et al. [4] investigated the permeability change due to mechanical loading, and they found that the microcrack growth can result in an increase of permeability for 3~5 orders of magnitudes. Many studies focused on the thermal effects on rock permeability subjected to high temperatures (>100°C) [13, 21–24]. For instance, Zharikov et al. [23] studied the permeability of amphibolite collected from Kola surface, Kola superdeep borehole (at depth of 11.4 km), and KTB borehole (at depth of 3.8 km) at temperatures up to 600°C. It was found that the permeability decreased first and then re-increased while the temperature exceeded 300°C during heating at low confining pressure (30–80 MPa). Microstructure investigations and analysis showed that the increase of rock matrix compressibility reduces the microcracks apertures from 100 to 300°C, and it leads to the decrease of the permeability. While temperature was higher than 300°C, the intensive microcrack generated by heat cracking at mineral grain boundaries was responsible to the increase of the permeability.
Compared with the extensive studies on permeability evolution in the dilatant regime and high temperatures which can result in generation and propagation of microcracks, few permeability tests and analysis have been performed under moderate temperatures in the compressive regime. However, the permeability variation of rock material under moderate pressures and temperatures is of great importance in some cases, for example, fluid movement at the shallow depth of the crust [25], long-term performance of an underground nuclear waste repository constructed in crystalline rock [3], weathering action of building claddings, and historical monuments due to seasonal temperature changes [26]. In this study, the transport properties of intact marble collected at depth of 1.6 km are measured under moderate temperatures and confining pressures. To explain the experimental results, a simple formula correlating permeability with porosity is proposed based on a penny-shaped crack linkage model. The applicability of the formula is verified against the measured results.
2. Tested Material and Experimental Procedure2.1. Tested Material
The tested material was taken from the diversion tunnel at the depth of 1600 m of the JinPing hydropower station in southwestern China. To measure the permeability under moderate temperature, two cylindrical samples (D20 and D30) with a diameter of 49.1 mm and a length of 80.0 mm were drilled from the same block. There are no preferred fabric textures of the tested material. The sample axes are parallel to the vertical axis of the ground. The density of the dry sample is 2.76 kg/m3, Young’s modulus is between 45 GPa and 60 GPa at high confining pressures (40 MPa to 60 MPa), and the uniaxial compressive strength is 103 MPa. The results of mineral X-ray diffraction show that the material consists of calcite (8.6% in weight), dolomite (91.0%), and mica (0.4%) [27].
The microstructure of the material has been investigated using optical microscopy and scanning electron microscopy (SEM). From the optical images examined under plane polarized light (Figure 1(a)), it can be seen that there is no preferred orientation of the minerals. The mineral sizes are mostly among 0.2 to 0.5 mm, and the transport path of fluid is mainly the interface of the grains. The parallel lines inside grains are the cleavage of calcite. The SEM images (Figure 1(b)) show that the rock is extremely compact. The grain boundaries in the SEM picture can only be distinguished by the variation of cleavage direction. The pore size and porosity of the marble are evaluated using mercury intrusion porosimetry. The porosities of the two samples are 0.17% and 0.14%, respectively. The pore size of the material mainly varies between 40 nm and 400 nm (Figure 2). It is close to the results of Italy and Sweden marble, in which the pore-size distributions cover a range of 2–200 nm [28].
Microstructure of the Jinping marble investigated by optical microscopes and SEM.
Pictures under optical microscope
Pictures under SEM
Pore size distribution of Jinping marble.
2.2. Experimental Procedure
The steady-state gas flow method is employed to investigate the transport property of the intact marble. The inlet gas pressure is maintained during the permeability test and the outlet pressure is the atmosphere pressure. A high-precision bubble gas flow meter is used to measure the flow rate of the outlet gas. The temperature is controlled using a heating system mounted around the sample cell. To seal the sample, a thin inner casing (nitrile butadiene rubber, thickness 0.3 mm) and a thick outer casing (thick Viton jacket, thickness 7 mm) are coated to prevent the oil from leaking into the sample. The likelihood of bypass flow that could take place at the interface between the sample and the inner casing may enlarge the flow rate substantially for tight rock [26]. In this study, liquid silica rubber (LSR) is used to stick the sample and inner casing together to eliminate the bypass flow.
Before permeability tests, the specimens were dried at 105°C in a vacuum oven for 24 hr. The sample D20 was tested under the hydrostatic pressures of 10, 20, and 30 MPa at room temperature (16°C), while the sample D30 under the constant pressure of 20 MPa was tested at the successively increasing temperature, that is, 15, 30, and 40°C. Several gas pressure levels were chosen to estimate the possible Klinkenberg effect.
When the gas permeability measurement is finished, the sample D20 was then put in MTS to measure the volume strain under the same mechanical loading conditions. Note that the low hydrostatic stress will not induce any damage of the sample.
3. Experimental Results
The steady gas flow velocity and the corresponding test conditions of the samples D20 and D30 are listed in Tables 1 and 2, respectively.
Test results on flow rate of sample D20 in different conditions.
Due to small pore size, Klinkenberg effect could not be ignored during gas permeability measurement. Klinkenberg in 1941 found that the permeability of a medium to gas is higher than that to liquid, and he attributed this phenomenon to “slip flow” between gas molecules and the pore wall surfaces. In Darcy flow, the collision between molecules and pore walls is negligible, and the flow velocity of the liquid is approximated as zero on the pore walls. However, the additional flux due to gas flow at the wall surfaces effectively increases as the pore radius approaches the mean free path of the gas molecules [2, 29]. This effect is expressed as follows: (1)Kg=K1+bp,where p is the average pore pressure and b is the Klinkenberg coefficient, which is related to pore geometry only.
The average pore pressure is often considered as the mean value of the inlet and outlet pressures and (1) can be expressed as follows:(2)Kg=K1+bp0+pL/2,where p0 and pL are the inlet and outlet pressures, respectively.
The use of (p0+pL)/2 in (2) is a rough estimation of the average gas pressure in the sample. Wu et al. [30] developed an analytical solution to estimate the permeability that incorporates the Klinkenberg effect. This solution is applicable when the sample temperature equals to the ambient air temperature. For tests performed on the sample D30, the sample temperature is different from the ambient air temperature, and the solution of Wu et al. [30] can be extended into(3)Qv,rprLμAp0-pLToTr=Kb+Kp0+pL2,where L (m) and A (m2) are the length and cross-sectional area of the test sample. Tr and Qv,r (m3/s) are the room temperature and the volumetric flux under room conditions. μ is the gas viscosity (Pa·s). To is the sample temperature.
During testing, the outlet pressure pL is kept constant while p0 is varied and Qv,r is measured. Then, K and b are evaluated by fitting results according to (3). Figure 3 illustrates the fitting curves in which Y=Qv,rprLμTo/ATr(p0-pL) is plotted against X=(p0+pL)/2. K equals the slope of the fitting line and Kb is the value of y intercept. The permeability, Klinkenberg coefficient, and correlation coefficient R2 are listed in Table 3. The results show that the Klinkenberg coefficient (b) increases as the permeability (K) decreases. It follows the same trend as the results from Tanikawa and Shimamoto [31].
Fitting results of permeability and Klinkenberg coefficient of samples D20 and D30.
Sample
Pressure or temperature
K/(10−21 m2)
b/MPa
R2
D20
10 MPa
21.1
1.64
0.993
20 MPa
11.0
3.02
0.999
30 MPa
6.92
3.55
0.994
D30
15°C
21.3
0.47
0.993
30°C
12.6
0.67
0.998
40°C
3.68
0.93
0.999
Fitting curve of test results for samples D20 and D30 by (3) (in the figure, X=(p0+pL)/2 and Y=Qv,rprLμTo/ATr(p0-pL)).
D20
D30
The permeability of Jinping marble decreases as the pressure and temperature increase. The permeability decreases from 21.1 × 10−21 m2 to 6.92 × 10−21 m2 as the pressure increases from 10 to 30 MPa; the decrease range is 67.2%. And it decreases from 21.3 × 10−21 m2 to 3.68 × 10−21 m2 as the temperature increases from 15 to 40°C; the decrease range is 82.7%. Comparing with the influence of the moderate pressure on the permeability, the effect of temperature on the hydraulic property of tight rock is substantial.
4. Microphysical Model and Test Results Analysis4.1. Microphysical Model
A significant number of studies have focused on developing the permeability evolution law in the dilatant regime based on percolation theory [32–34] or linear elastic fracture mechanics [35]. Little attention has been paid to permeability evolution in the compressive regime. In this study, the measured outflow of gas indicates the preexistence of linkage network of cracks in the samples. The variation of permeability during compression is due to the change of crack apertures.
Dienes [32] employed the continuum-averaging method and Poiseuille’s law for laminar flow between parallel plates to develop an analytical result for the permeability tensor of rock intersected by isotropically distributed penny-shaped cracks. The permeability tensor is isotropic and is expressed as follows:(4)K=8π215θA¯3N0c¯5,where θ accounts for the deviations of the crack shape from a uniform lamina, containing the drag effects of crack shape, crack end, and the roughness, A¯ is the mean aspect ratio (w¯/c¯, where w¯ is the mean crack half aperture (m) and c¯ is the mean crack radius (m)) of the cracks, and N0 denotes the number of cracks per unit volume that are not isolated.
Suppose the porosity (ϕ) of rock sample is composed of penny-shaped cracks which were due to rock damage during tectonic processes, sample collection, preparation, and so forth. It can be estimated from the mean crack aperture, radius, and spacing as follows [36]:(5)ϕ=2παc¯2w¯l¯3,where l¯ is the mean crack spacing (m) and α is the volumetric shape factor giving the deviation in volume from the ideal penny-shaped crack, defined as (real crack volume)/(ideal penny-shaped crack volume).
The connection between the permeability and the porosity (ϕ) is easily obtained from (4) and (5):(6)K=θN015πα3·l¯3c¯23·ϕ3.
In the compressive regime, the values of N0, l¯, c¯, and α in (6) remain unchanged because no new microcrack initiates and propagates. The value of θ may change during compression. However, considering the averaging effect of many cracks, θ is assumed to be constant here. Therefore, the item θN0/15πα3·l¯3/c¯23 is considered as a constant during compression, and it is defined as β, which has no relations to the apertures of cracks. So, (6) can be written as(7)K=βϕ3.
Equation (7) establishes the relationship between the permeability and the porosity (ϕ). It should be noted here that when crack propagates or new crack initiates, the value of β is no longer a constant.
4.2. Test Results Analysis
It is impossible to measure the porosity directly during testing. However, the value of porosity (ϕ) can be estimated from the initial porosity ϕ0, volumetric strain εv, and mineral compressibility Cs:(8)ϕ=ϕ0-εv+CsP.
So, (7) can be written as(9)K=βϕ0-εv+CsP3.
It can be seen from (9) that the intrinsic permeability decreases following a cubic relationship as the compressive volumetric strain increases.
The solid line in Figure 4 presents the volumetric strain of sample D20 as pressure increases. In (9), the compressibility (Cs) of dolomite can be obtained from its bulk modulus [37], and its value is 1.06×10-5/MPa. The compressive volumetric strains (εv) at 10 MPa, 20 MPa, and 30 MPa are obtained from the volumetric strain curve in Figure 4, and they are εv10 = 0.00071, εv20 = 0.0011, and εv30 = 0.0013, respectively. Combining the permeability test results (K in Table 3) at 10 MPa, 20 MPa, and 30 MPa confining pressures (P), the parameters of β and ϕ0 can be fitted from (9). (10)β=0.94×10-11m2,ϕ0=0.19×10-2.
Permeability variation with volumetric strain during compression.
The initial porosity ϕ0 determined from (9) is close to the intrusion test results (0.17% and 0.14%). Equation (9) can be used to predict variation of rock permeability from the variation of volumetric strain. The estimated permeability curve for D20 is plotted against volumetric strain in Figure 4 as the dot line. It should be noted that the linear Hook’s law should not be applied to calculate the volumetric strain because the permeability variation law is based on the variation of crack aperture other than the elastic deformation of matrix. Replacing β in (7), the porosity can be estimated and it decreases from 0.13% to 0.09% as the pressure increases from 10 MPa to 30 MPa; the decrease range is 31%.
For D30, the permeability variation with temperature can be estimated by replacing the term CsP in (9) with α(T-T0), where α is the thermal expansion coefficient of the mineral; T and T0 are the current temperature and initial temperature. In the current test, the volumetric strains are not measured at different temperatures. So, the parameters β cannot be obtained from the test results of D30. Supposing that samples D30 and D20 have the same value of β, which means that cracks in both specimen follow the same distribution, the variation of porosity ϕ with temperatures can be roughly estimated from (7), as shown in Figure 5. The results show that the porosity decreases from 0.13% to 0.07% as the temperature increases from 15°C to 40°C; the decrease range is 44%. Equation (5) shows that the mean crack aperture (2w¯) changes linearly with the porosity. Thus, the mean aperture also decreases 44% as temperature increases from 15°C to 40°C. The porosity change due to the thermal effect is also found by Sun et al. [38], who found that the porosity of granite decreases (from 0.88% to 0.75%) as temperature increases (from 25°C to 50°C) in the moderate temperature range. Comparing with the influence of the moderate pressure on the permeability presented above, the thermal effect on the permeability is notable and it should be taken into account into the estimation of the transport properties of intact rock.
Variation of porosity with temperature under hydrostatic stress (20 MPa) for Jinping marble.
5. Conclusions
The steady state method is used to study the permeability of Jinping marble under moderate pressures and temperatures in the compressive regime. The experimental results show that the permeability of Jinping marble decreases from 21.1 × 10−21 m2 to 6.92 × 10−21 m2 as the pressure increases from 10 to 30 MPa and decreases from 21.3 × 10−21 m2 to 3.68 × 10−21 m2 as the temperature increases from 15 to 40°C. The results imply that the influence of temperature on marble permeability is as remarkable as stress.
The permeability evolution law in the compressive regime is studied based on a geometric model of the penny-shaped crack distribution. The analytical result shows that the permeability decreases following a cubic law. The initial porosity of the marble estimated from the permeability tests is 0.19%, which is close to the results of the mercury intrusion porosimetry (0.14% and 0.17%). The analytical results show that the porosity decreases from 0.13% to 0.09% as the pressure increases from 10 MPa to 30 MPa, and from 0.13% to 0.07% as the temperature increases from 15°C to 40°C. From the permeability evolution law, the permeability for crystalline rock in the compressive regime due to moderate pressure and temperature can be estimated from the volumetric strain curve. The proposed permeability variation law will be furtherly verified to investigate the effect of moderate stress and temperature on the transport properties of host rock of underground nuclear waste repository in the future.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors gratefully acknowledge the support of the Chinese Fundamental Research (973) Program through Grant no. 2013CB03600.
CosenzaP.GhoreychiM.Bazargan-SabetB.de MarsilyG.In situ rock salt permeability measurement for long term safety assessment of storage199936450952610.1016/s0148-9062(99)00017-02-s2.0-0033150498YangD.WangW.ChenW.WangS.WangX.Experimental investigation on the coupled effect of effective stress and gas slippage on the permeability of shale201772-s2.0-8501574090310.1038/srep4469644696TsangC.-F.BernierF.DaviesC.Geohydromechanical processes in the excavation damaged zone in crystalline rock, rock salt, and indurated and plastic clays—in the context of radioactive waste disposal200542110912510.1016/j.ijrmms.2004.08.0032-s2.0-10944234687ChenY.HuS.WeiK.HuR.ZhouC.JingL.Experimental characterization and micromechanical modeling of damage-induced permeability variation in Beishan granite20147164762-s2.0-8490538720810.1016/j.ijrmms.2014.07.002WitherspoonP. A.GaleJ. E.Mechanical and hydraulic properties of rocks related to induced seismicity197711123552-s2.0-001747110310.1016/0013-7952(77)90018-7EtheridgeM. A.WallV. J.CoxS. F.VernonR. H.High fluid pressures during regional metamorphism and deformation: implications for mass transport and deformation mechanisms1984896434443582-s2.0-0020138311AgueJ. J.Models of permeability contrasts in subduction zone mélange: implications for gradients in fluid fluxes, Syros and Tinos Islands, Greece20072393-421722710.1016/j.chemgeo.2006.08.0122-s2.0-34147137153FarrellN. J. C.HealyD.TaylorC. W.Anisotropy of permeability in faulted porous sandstones20146350672-s2.0-8489626355310.1016/j.jsg.2014.02.008KolzenburgS.RussellJ. K.Welding of pyroclastic conduit infill: a mechanism for cyclical explosive eruptions201411975305532310.1002/2013JB0109312-s2.0-84906089319BoulinP. F.BretonnierP.GlandN.LombardJ. M.Contribution of the steady state method to water permeability measurement in very low permeability porous media20126733874012-s2.0-8486396471310.2516/ogst/2011169BoltonE. W.LasagaA. C.RyeD. M.Long-term flow/chemistry feedback in a porous medium with heterogenous permeability: kinetic control of dissolution and precipitation1999299116810.2475/ajs.299.1.12-s2.0-0033408194DongJ.-J.HsuJ.-Y.WuW.-J.ShimamotoT.HungJ.-H.YehE.-C.WuY.-H.SoneH.Stress-dependence of the permeability and porosity of sandstone and shale from TCDP Hole-A2010477114111572-s2.0-7795665260710.1016/j.ijrmms.2010.06.019KožušníkováA.KonečnýP.Influence of temperature on the permeability of rocks201161121081108510.1680/geot.8.T.0342-s2.0-80955158220ChenJ.YangX.DuanQ.PeachC.Integrated measurements of permeability, effective porosity, and specific storage of core samples using water as the pore fluid20157955622-s2.0-8494004631010.1016/j.ijrmms.2015.08.007RathnaweeraT. D.RanjithP. G.PereraM. S. A.Effect of salinity on effective CO2 permeability in reservoir rock determined by pressure transient methods: an experimental study on hawkesbury sandstone20154852093211010.1007/s00603-014-0671-02-s2.0-84939563138XuP.YangS.-Q.Permeability evolution of sandstone under short-term and long-term triaxial compression2016851521642-s2.0-8496191402010.1016/j.ijrmms.2016.03.016WangL. L.ZhangG. Q.HallaisS.TanguyA.YangD. S.The swelling of shales: a multi-scale experimental investigation201710.1021/acs.energyfuels.7b01223SouleyM.HomandF.PepaS.HoxhaD.Damage-induced permeability changes in granite: a case example at the URL in Canada200138229731010.1016/s1365-1609(01)00002-82-s2.0-0035251433SchulzeO.PoppT.KernH.Development of damage and permeability in deforming rock salt2001612-31631802-s2.0-003542502910.1016/S0013-7952(01)00051-5JiangT.ShaoJ. F.XuW. Y.ZhouC. B.Experimental investigation and micromechanical analysis of damage and permeability variation in brittle rocks20104757037132-s2.0-7795403486210.1016/j.ijrmms.2010.05.003ChakiS.TakarliM.AgbodjanW.Influence of thermal damage on physical properties of a granite rock: porosity, permeability and ultrasonic wave evolutions20082271456146110.1016/j.conbuildmat.2007.04.0022-s2.0-43049112449ShmonovV. M.VitovtovaV. M.ZharikovA. V.Experimental study of seismic oscillation effect on rock permeability under high temperature and pressure19993634054122-s2.0-003278914310.1016/S0148-9062(99)00020-0ZharikovA. V.VitovtovaV. M.ShmonovV. M.GrafchikovA. A.Permeability of the rocks from the Kola superdeep borehole at high temperature and pressure: implication to fluid dynamics in the continental crust20033701-417719110.1016/S0040-1951(03)00185-92-s2.0-0042472824ChenS.YangC.WangG.Evolution of thermal damage and permeability of Beishan granite2017110153315422-s2.0-8498795662710.1016/j.applthermaleng.2016.09.075ManningC. E.IngebritsenS. E.Permeability of the continental crust: implications of geothermal data and metamorphic systems199937112715010.1029/1998rg9000022-s2.0-0033010497SelvaduraiA. P. S.CarnaffanP.A transient pressure pulse method for the measurement of permeability of a cement grout19972434895022-s2.0-003070860610.1139/l96-132ZhaoX.ZuoJ.PeiJ.Meso-experimental study of fracture mechanism of bedded marble in Jinping20123135345422-s2.0-84859384579Malaga-StarzecK.LindqvistJ. E.SchouenborgB.Experimental study on the variation in porosity of marble as a function of temperature20022051818810.1144/GSL.SP.2002.205.01.072-s2.0-0037004573DullienF. A. C.19922ndNew York, NY, USAAcademic PressWuY.-S.PruessK.PersoffP.Gas flow in porous media with klinkenberg effects199832111713710.1023/A:10065352116842-s2.0-0345601543TanikawaW.ShimamotoT.Comparison of Klinkenberg-corrected gas permeability and water permeability in sedimentary rocks20094622292382-s2.0-5834910500710.1016/j.ijrmms.2008.03.004DienesJ. K.GoodmanR. E.HeuzeF. E.Permeability, percolation and statistical crack mechanics1982New York, NY, USAAmerican Institute of Mining, Metallurgical and Petroleum Engineers8694GueguenY.DienesJ.Transport properties of rocks from statistics and percolation19892111132-s2.0-000180067910.1007/BF00897237ZhuW.WongT.-F.Network modeling of the evolution of permeability and dilatancy in compact rock199910422963297110.1029/1998JB9000622-s2.0-0033540636SimpsonG.GuéguenY.SchneiderF.Permeability enhancement due to microcrack dilatancy in the damage regime200110633999401610.1029/2000JB9001942-s2.0-0034952225PeachC. J.SpiersC. J.Influence of crystal plastic deformation on dilatancy and permeability development in synthetic salt rock19962561-41011282-s2.0-002975179010.1016/0040-1951(95)00170-0AhrensT. J.1995Washington, DC, USAAmer Geophysical UnionSunQ.ZhangW.XueL.ZhangZ.SuT.Thermal damage pattern and thresholds of granite2015743234123492-s2.0-8493751401410.1007/s12665-015-4234-9