Influence of flow velocity on conservative solute transport in sets of parallel fractures with fracture skin

Analysis of contaminant transport through fractured crystalline rocks has received considerable attention, particularly with regard to subsurface nuclear waste repositories. Most of the studies have employed the dual continuum approach, with the fractures and the rock matrix as the two continuums, assuming that fractures control the overall conductivity of the rock and the porous matrix just provides storage. However, field observations of rock fractures have shown that the real situation can be very complex. Based on some recent investigations, it has been reported that the portion of the rock matrix adjacent to many open fractures is physically and chemically altered. These alterations, referred to as the fracture skin, can have different sorption and diffusion properties compared to those of the undisturbed rock matrix and this may influence the transport of solutes through such formations. In the present study, a numerical model is developed to simulate conservative solute transport in a fractured crystalline rock formation using the triple continuum approach ?? with the fracture, fracture skin and the rock matrix as the three continuums. The model is solved using a fully implicit finite difference scheme. Contaminant migration in the fractured formation with and without skin has been simulated. It is observed that contaminant penetration along the fracture is enhanced at large flow velocities. The effect of flow velocity on conservative solute transport is investigated for different fracture apertures and fracture skin thicknesses. The influence of flow velocity on contaminant transport is demonstrated to be more with change in fracture aperture than with change in skin thickness.     

Investigations on Diffusion Characteristics of Granite and Chalk Rock Mass

Contaminant transport through fractured rock mass is predominated by diffusion. This is due to the continuous interaction of the mobile water present in the fracture network and relatively immobile pore water, which is adsorbed on the surface and in the rock matrix itself. Even though the advective flow through the fracture network is high, besides sorption of rock mass, the diffusive exchange into the rock mass leads to significant retardation of contaminant transport. Hence, for describing contaminant transport in fractured rock mass, more precisely, the effect of retardation attributed to the matrix diffusion must be taken in account. With this in view, a methodology, which can be employed for determination of the diffusion characteristics of the rock mass, has been developed and its details are presented in this paper. Validation of the methodology has been demonstrated with the help of Archie’s law.     

Permeability and stress in crystalline rocks

Groundwater from crystalline rocks is a significant resource in many areas of the world. It is also an important medium for contaminant transport from, for example, deep nuclear waste repositories. Stress distributions in fractured rocks are important in controlling groundwater flow in several ways: (i) palaeostress fields are responsible for the evolution of fracture systems which transmit groundwater; (ii) current in situ stress fields will influence the shape and aperture of fractures; (iii) humans can influence the natural stress field in a rock mass to enhance fracture flows. The significance of stresses for groundwater flow can be investigated by field techniques (hydraulic fracturing), laboratory techniques (stress cells) or by numerical modelling.     

Constraining the uncertainty in fracture geometry using tracer tests

In fractured-rock aquifers, the geometric and hydraulic properties of the fractures commonly have a dominant influence on transport. Tracer tests are often used to estimate directly the gross transport properties of a fractured rock mass. The prospects for understanding characteristics of the heterogeneities in a fractured porous medium were explored from evidence provided by tracer experiments. The approach was to simulate flow and transport on a large set of prescribed fracture networks in a two-dimensional homogeneous permeable medium, thus generating synthetic tracer test data. The fracture orientation, aperture, spacing and network geometry were systematically altered from one case to the next. A classification scheme was devised for the tracer breakthrough curves using principal component analysis and this classification was linked to the fracture pattern properties. Even under highly simplified and controlled conditions, quite different fracture patterns can produce very similar breakthrough curves. The classification scheme thus demonstrates that a single breakthrough curve cannot reveal the fracture geometry with any precision. However, the scheme provided a methodology for rejecting geometric properties that do not belong to the fracture pattern under investigation, thus reducing the uncertainty in fracture geometry.     

Hybrid analytical and finite element numerical modeling of mass and heat transport in fractured rocks with matrix diffusion

Quantification of mass and heat transport in fractured porous rocks is important to areas such as contaminant transport, storage and release in fractured rock aquifers, the migration and sorption of radioactive nuclides from waste depositories, and the characterization of engineered heat exchangers in the context of enhanced geothermal systems. The large difference between flow and transport characteristics in fractures and in the surrounding matrix rock means models of such systems are forced to make a number of simplifications. Analytical approaches assume a homogeneous system, numerical approaches address the scale at which a process is operating, but may lose individual important processes due to averaging considerations. Numerical stability criteria limit the contrasts possible in defining material properties. Here, a hybrid analytical–numerical method for transport modeling in fractured media is presented. This method combines a numerical model for flow and transport in a heterogeneous fracture and an analytical solution for matrix diffusion. By linking the two types of model, the advantages of both methods can be combined. The methodology as well as the mathematical background are developed, verified for simple geometries, and applied to fractures representing experimental field conditions in the Grimsel rock laboratory.     

Modeling of subsidence and stress-dependent hydraulic conductivity for intact and fractured porous media

Summary This study investigates the changes in deformation and stress dependent hydraulic conductivities that occur as a result of underground mining in intact and fractured porous media. The intact porous medium is assumed to be comprised of regularly packed spherical grains of uniform size. The variation in grain size or pore space due to the effect of changing intergranular stresses results in a change in rock hydraulic conductivity. A model is developed to describe the sensitivity of hydraulic conductivity to effective stresses through Hertzian contact of spherical grains. The fractured porous medium is approximated as an equivalent fracture network in which a single fracture is idealized as a planar opening having a constant equivalent thickness or aperture. Changes in fracture aperture as a result of changes in elastic deformation control the variation of hydraulic conductivity. A model is presented to illustrate the coupling between strain and hydraulic conductivity. Subsidence induced deformations that result from mining induced changes in hydraulic conductivity in both intact and fractured media. These changes are examined and compared with results from a mining case study.     

Boundary element analysis of contaminant transport in fractured porous media

A two-dimensional boundary integral method to analyse the flow of contaminant in fractured media having a two- or three-dimensional orthogonal fracture network is presented. The method assumes that the fractures provide the paths of least resistance for transport of contaminants while the matrix, because of its low permeability, acts as ‘storage blocks’ into which the contaminant diffuses. Laplace transform is used to eliminate the time variable in the governing equation in order to facilitate the formulation of a boundary integral equation in the Laplace transform space. Conventional boundary element techniques are applied to solve for the contaminant concentrations at specified locations in the spatial domain. The concentration in the time domain is then obtained by using an efficient inversion technique developed by Talbot. The method is able to analyse the behaviour of waste repositories which have diminishing concentration due to the mass transport of the contaminant into the surrounding fractured media.     

非平衡-非线性吸附情况下填埋场污染物运移分析

采用Langmuir等温吸附线方程描述非线性吸附性能, 基于改进的混合元方法，通过数值计算与分析，探讨了非平衡吸附条件下污染物运移过程及其机理。计算结果表明: 当考虑非平衡、非线性吸附性能时, 污染物穿透曲线即浓度的时程变化曲线尖锐而狭窄、峰值点前移, “拖长尾”现象不明显, 由此说明，土颗粒对污染物的非平衡、非线性吸附使得污染物的穿透能力增强, 滞留能力下降。进一步的变动参数比较分析表明：Langmuir等温线方程中的参数B、压实粘土衬里的渗透性及地下水渗流速度对污染物运移过程具有显著的影响。     

Numerical modeling of stress effects on solute transport in fractured rocks

The effects of stress/deformation on fluid flow and contaminant transport in fractured rocks is one of the major concerns for performance and safety assessments of many subsurface engineering problems, especially radioactive waste disposal and oil/gas reservoir fields. However, very little progress has been made to study this issue due to difficulties in both experiments and numerical modeling. The objective of this study is to systematically investigate the influence of stress on solute transport in fractured rocks for the first time, considering different stress and hydraulic pressure conditions. A hybrid approach combining discrete element method (DEM) for stress-flow simulations and a particle tracking algorithm is developed. The impact of matrix diffusion (diffusion of molecular size solutes in and out of the rock matrix, and sorption onto the surface of micropores in rock matrix) is also included. The numerical results show that stress not only significantly changes the solute residence time through the fracture networks, but also changes the solute travel paths. Matrix diffusion plays a dominant role in solute transport when the hydraulic gradient is small, which is often encountered in practice.     

Application of multiple-point geostatistics on modelling groundwater flow and transport in a cross-bedded aquifer (Belgium)

Sedimentological processes often result in complex three-dimensional subsurface heterogeneity of hydrogeological parameter values. Variogram-based stochastic approaches are often not able to describe heterogeneity in such complex geological environments. This work shows how multiple-point geostatistics can be applied in a realistic hydrogeological application to determine the impact of complex geological heterogeneity on groundwater flow and transport. The approach is applied to a real aquifer in Belgium that exhibits a complex sedimentary heterogeneity and anisotropy. A training image is constructed based on geological and hydrogeological field data. Multiple-point statistics are borrowed from this training image to simulate hydrofacies occurrence, while intrafacies permeability variability is simulated using conventional variogram-based geostatistical methods. The simulated hydraulic conductivity realizations are used as input to a groundwater flow and transport model to investigate the effect of small-scale sedimentary heterogeneity on contaminant plume migration. Results show that small-scale sedimentary heterogeneity has a significant effect on contaminant transport in the studied aquifer. The uncertainty on the spatial facies distribution and intrafacies hydraulic conductivity distribution results in a significant uncertainty on the calculated concentration distribution. Comparison with standard variogram-based techniques shows that multiple-point geostatistics allow better reproduction of irregularly shaped low-permeability clay drapes that influence solute transport.     

Numerical Modeling of Natural Fracture Distributions in Shale

The production efficiency of shale gas is affected by the interaction between hydraulic and natural fractures. This study presents a simulation of natural fractures in shale reservoirs, based on a discrete fracture network (DFN) method for hydraulic fracturing engineering. Fracture properties of the model are calculated from core fracture data, according to statistical mathematical analysis. The calculation results make full use of the quantitative information of core fracture orientation, density, opening and length, which constitute the direct and extensive data of mining engineering. The reliability and applicability of the model are analyzed with regard to model size and density, a calculation method for dominant size and density being proposed. Then, finite element analysis is applied to a hydraulic fracturing numerical simulation of a shale fractured reservoir in southeastern Chongqing. The hydraulic pressure distribution, fracture propagation, acoustic emission information and in situ stress changes during fracturing are analyzed. The results show the application of fracture statistics in fracture modeling and the influence of fracture distribution on hydraulic fracturing engineering. The present analysis may provide a reference for shale gas exploitation.     

Permeability tensor and representative elementary volume of fractured rock masses

Based on a simulation of three-dimensional fracture networks and a superposition principle of liquid dissipation energy for fractured rock masses, a model of the fracture permeability tensor is proposed. An elastic constitutive model of rock fractures, considering fracture closure and dilation during shearing, is also proposed, based on the dilation angle of the fracture. Algorithms of flow-path searching and calculation of the effective flow coefficients for fracture networks are presented, together with a discussion on the influence of geometric parameters of the fractures (trace length, spacing, aperture, orientation and the number of fracture sets) on magnitude, anisotropy of hydraulic permeability and the size of a representative elementary volume (REV). The anisotropy of hydraulic permeability of fractured rock masses is mainly affected by orientation and the number of fracture sets, and the REV size is mainly influenced by trace length, spacing and the number of fracture sets. The results of studies on REV size and the influence of in-situ stress on hydraulic conductivity of the rock mass on the slope of Jinping-I hydropower station, China, are presented using the developed models and methods. The simulation results agreed well with the results obtained from field water-pressure measurements, with an error of less than 10 %.     

Modelling coupled heat and contaminant transport in groundwater

This paper presents the results of a research programme conducted on the geotechnical centrifuge at The University of Western Australia to investigate coupled heat and contaminant transport in the soil surrounding a buried waste source. The phenomena which govern heat and contaminant transport through porous media are discussed, the principles of geotechnical centrifuge modelling are outlined, and relevant scaling laws that govern the relationship between a centrifuge model and the prototype, with respect to the problem of coupled waste transport, are presented. A model test, simulating two-dimensional migration from a buried heat and contaminant source, is described, and the results from four model tests are presented. The experimental data show that hydraulic instability is responsible for the transport of contaminant in the soil around the source and that the mode of instability is determined by the magnitude of the effective Rayleigh number.     

运移参数对GCL中污染物运移的影响研究

分析了运移参数的变化对GCL中污染物运移的影响。从衬垫底部浓度比、对流和扩散通量、污染物累积运移量等方面,对污染物在GCL中的运移进行了计算。结果表明,减小渗透系数和扩散系数可以有效减缓衬垫底部浓度比的增加速度。渗透系数对对流通量的大小、扩散通量的变化速度的影响较大;扩散系数主要影响最大扩散通量。当渗透系数较大时,可以忽略扩散累积运移量,用对流累积运移量代替总累积运移量引起的误差很小。当渗透系数较小时,忽略扩散运移量可能产生较大的误差,尤其是在扩散系数也较大的情况下误差更加明显。      

裂隙岩体渗透系数确定方法研究

裂隙岩体渗透系数以及渗透主方向的确定对研究岩体渗透性大小及各向异性具有重要意义。高放废物地质处置库介质岩体的渗透性能将直接影响其使用安全性。本文运用离散裂隙网络模拟的方法对我国高放废物处置库甘肃北山预选区3#钻孔附近裂隙岩体进行了渗透性质分析。通过对3#钻孔171.5~178.0m段压水试验数据的反演,标定了离散裂隙网络渗流模型中的裂隙渗透参数(导水系数T)。利用标定的离散裂隙网络模型对场区裂隙岩体进行了渗流模拟,确定了该区域裂隙岩体的渗流表征单元体(REV)的尺寸大小以及渗透主值和主渗透方向。运用离散裂隙网络模型计算得出的渗透主值的几何均值与现场压水试验计算结果较接近,证明了计算结果的有效性。     

岩体水力学基础(二)─—岩体水力学的基础理论

单裂隙水流立方定律是岩体水力学中最基本的．也是最重要的定理。本文从不变形平直等宽单裂隙水流运动定理开始，讨论了受应力作用岩体单裂隙、一组平行裂隙、一组正交裂隙以及岩体裂隙系统中水流运动规律，分析了单裂隙中水流渗透压力与裂隙变形之关系。运用水力学原理推导了岩溶单管道流定理，为岩体水力学模型研究奠定了基础。     

裂隙岩体渗透性影响因素的研究

采用UDEC离散单元法中关于裂隙岩体开挖模拟及水力全耦合分析模型,分析裂隙岩体洞室开挖后,因围岩应力与水力耦合作用导致裂隙隙宽变化及渗流变化的过程。为了更直观地了解耦合作用对裂隙岩体渗透特性的影响,以隧洞开挖为例,用开挖后隧洞内总涌水量来表征岩体的渗透特性。利用数值试验的方法,研究了块体边界大小、初始应力比、裂隙隙宽和裂隙夹角对开挖后隧洞内涌水量变化的影响,进而可以看出它们对裂隙岩体渗透性的影响。并得出如下结论：随着块体尺寸和初始应力比的增大,隧洞内总涌水量减少;随初始隙宽的增大涌水量增加并当达到某一固定值时保持不变;隧洞涌水量在θ2/θ1=3.5,其中θ1=30°,即两组节理的夹角为75°处达到最大。     

Sensitivity analysis of hydrological parameters in modeling flow and transport in the unsaturated zone of Yucca Mountain, Nevada, USA

The unsaturated fractured volcanic deposits at Yucca Mountain in Nevada, USA, have been intensively investigated as a possible repository site for storing high-level radioactive waste. Field studies at the site have revealed that there exist large variabilities in hydrological parameters over the spatial domain of the mountain. Systematic analyses of hydrological parameters using a site-scale three-dimensional unsaturated zone (UZ) flow model have been undertaken. The main objective of the sensitivity analyses was to evaluate the effects of uncertainties in hydrologic parameters on modeled UZ flow and contaminant transport results. Sensitivity analyses were carried out relative to fracture and matrix permeability and capillary strength (van Genuchten α) through variation of these parameter values by one standard deviation from the base-case values. The parameter variation resulted in eight parameter sets. Modeling results for the eight UZ flow sensitivity cases have been compared with field observed data and simulation results from the base-case model. The effects of parameter uncertainties on the flow fields were evaluated through comparison of results for flow and transport. In general, this study shows that uncertainties in matrix parameters cause larger uncertainty in simulated moisture flux than corresponding uncertainties in fracture properties for unsaturated flow through heterogeneous fractured rock.     

微生物灌浆加固裂隙岩体的渗流特性分析

裂隙岩体的防渗加固是影响工程稳定和安全的关键技术问题之一。在以往微生物加固砂土技术基础上,将微生物加固技术应用于裂隙岩体灌浆加固,研究结果表明,(1) 加固前围压对裂隙岩样的水力开度影响很大,加固后裂隙被碳酸钙沉积物充填,岩样的渗流由裂隙渗流转变为孔隙渗流,单位渗流量和渗透系数由围压和渗透水压共同控制;(2) 微生物灌浆加固后,不同围压和渗透水压力作用下裂隙岩样的单位时间渗流量减小了80.12%～90.04%,渗透系数可达到10–6 cm/s数量级;(3) 裂隙岩样灌浆加固过程中微生物诱导产生的CaCO3沉积物具有较好胶结作用,将劈裂岩样胶结为一个整体,达到了加固防渗效果。研究成果可为裂隙岩体的灌浆加固提供较好的参考。     

离散裂隙渗流方法与裂隙化渗透介质建模

流体渗流模拟的连续介质方法通常适用于多孔地质体，并不一定适用于裂隙岩体，由于裂隙分布及其特征与孔隙差异较大。若流体渗流主要受裂隙的控制，对于一定尺寸的裂隙岩体，多孔介质假设则较难刻划裂隙岩体的渗流特征。离散裂隙渗流方法不但可直接用于模拟裂隙岩体非均质性和各向异性等渗流特征，而且可用其确定所研究的裂隙岩体典型单元体及其水力传导（渗透）张量大小。主要讨论了以下问题：（1）饱和裂隙介质中一般的离散流体渗流模拟；（2）裂隙岩体中的REV（典型单元体）及其水力传导（渗透）张量的确定；（3）利用离散裂隙网络流体渗流模型研究裂隙方向几何参数对水力传导系数和REV的影响；（4）在二维和三维离散裂隙流体渗流模型中对区域大裂隙和局部小裂隙的处理方法。调查结果显示离散裂隙流体渗流数学模型可用来评价不同尺度上的裂隙岩体的水力特征，以及裂隙方向对裂隙化岩体的水力特征有着不可忽视的影响。同时，局部小裂隙、区域大裂隙应当区别对待，以便据其所起的作用及水力特征，建立裂隙化岩体相应的流体渗流模型。