岩体粗糙单裂隙渗流及溶质运移研究

与多孔介质相比,裂隙介质中的渗流、溶质运移要复杂得多,许多学者在这方面进行了深入的研究,他们在模型的建立、求解、验证等方面都取得很大突破,但仍有许多难点没有解决。着眼于粗糙单裂隙开度变化的实质特性,将原有裂隙介质渗流、溶质运移方程进行修改,建立粗糙单裂隙渗流、溶质运移数学模型,并进行数值模拟。     

裂隙岩体地下水溶质运移的尺度问题

核废料的安全贮存、地下污染的预测控制等目前面临的许多重大问题要求对大尺度裂隙岩体地下水溶质运移问题进行研究,然而裂隙岩体地下水溶质运移问题的尺度效应给这项研究带来了极大的挑战.本文介绍了裂隙岩体地下水溶质运移的多尺度概念模型及其物理和数学模拟方法,并探讨了模拟尺度与观测尺度、预测尺度、裂隙介质尺度之间的关系及尺度效应的分形特征.     

岩溶区地下水数值模拟研究进展

岩溶含水介质的不均一性导致岩溶地下水流动、溶质运移和热量迁移的数学模拟研究成为地下水模拟的难点。本文综述了岩溶区地下水流模拟的几种方法,重点阐述了等效多孔介质法、双重连续介质法和三重介质法的定义、发展过程和适用范围,并回顾了这几种方法的研究成果。从等效多孔介质法到三重介质法,模拟精度不断提高,适用范围也逐渐由大区域实际问题向小区域理论研究过渡。介绍了溶质运移模拟和热迁移模拟的研究方法及实例。溶质运移模拟以对流弥散方程为基础,其中尺度效应是溶质运移模拟的重点研究问题;热量迁移模拟应考虑地下热水密度变化对地下热水运动的影响。溶质运移模拟和热量迁移模拟往往是将迁移模型和已经调试成功的地下水流动模型相耦合,从而达到模拟溶质及热量迁移的目的。由于溶质运移和热量迁移的复杂性,现阶段水流模型多数处于等效多孔介质模型阶段。综合理论及实际应用,指出精确刻画裂隙及管道和注重基础数学算法是岩溶水数值模拟进步的关键。      

岩溶管道与裂隙介质间溶质暂态存储机制

溶质暂态存储是岩溶地下水溶质运移过程中的普遍现象。为揭示岩溶管道与裂隙介质间溶质暂态存储机制,本文构建室内管道-裂隙物理模型,开展集中补给条件下的定量示踪试验,运用双区对流弥散模型实现溶质运移过程模拟。研究表明：随着集中补给水动力条件的增强,裂隙暂态存储水量呈线性增加趋势,溶质穿透曲线由单峰型向双峰型转变;管道和裂隙中的平均流速呈负相关关系,溶质在管道和裂隙中的滞留时间差决定了穿透曲线的形态;溶质暂态存储引发了穿透曲线的拖尾效应和双峰现象,对岩溶地下水溶质运移过程具有重要的控制作用。     

煤矿区放水试验中示踪剂的运移规律

以地下水动力学的相关理论为基础,提出应用溶质运移模型对放水和示踪联合试验中获取的水量、水位及示踪剂浓度时间序列数据进行定量研究的方法。对裂隙介质溶质运移的基本规律、数学模型及其数值解法进行了系统的总结,对比各种解法的优缺点及其适用条件,并给出了联合试验溶质运移规律模拟分析的流程。     

地下多组分反应溶质运移数值模拟:地质资源和环境研究的新方法

地下多组分反应溶质运移数值模拟(RTM)是解释地球系统中的耦合过程和不同时空尺度对其影响的重要工具。RTM是研究地球科学基础理论、地质资源和环境等复杂地球化学过程的一个新方法,可用于如废物处置安全性评估、地下水污染研究、二氧化碳地质储存、金属矿床的地浸开采等的研究中。笔者首先回顾了反应溶质溶质运移模拟的发展历史,然后总结了反应溶质运移模拟的发展现状,再从耦合过程、空间尺度、裂隙和非均质介质处理角度说明了反应溶质运移模拟所面临的挑战。结合地下水质的演化、生物降解、CO2地质储存等具体实例讨论了反应溶质运移模拟的广泛应用前景,探讨了反应溶质运移模拟的未来发展方向。     

裂隙网络管道模型弥散试验

为了探求不同裂隙几何参数对裂隙网络溶质运移的影响,基于离散裂隙网络思想和优势流、沟槽流理论,建立裂隙管道网络概念模型,搭建不同管径、不同连通方式的管道网络试验装置,进行渗流和溶质运移实验。运用应用广泛的模拟软件CHEMFLO-2000建立等效多孔介质模型,拟合不同几何参数下等效弥散度,定量刻画不同管道网络几何参数对溶质运移的影响,讨论了不同管径、连通方式等与等效弥散度之间的关系。通过进一步分析得出:在连通方式相同的情况下,不同管径的管道网络等效弥散度存在差异,但是差异不大;溶质在小管径中的穿透时间短于大管径,穿透曲线缓和程度相差不大;管道网络连通方式越复杂,其等效弥散度越大、对溶质运移的影响越大、穿透曲线越缓和、路径越长,等效弥散度越大;用这种等效弥散度的方法表征管道网络对溶质运移的影响,与多孔介质弥散度具有相似性;管道数目、管道面数目与等效弥散度成正相关关系,且等效弥散度随尺度的增大而增加。     

地表地下耦合的水量及溶质运移模拟研究进展

地表水与地下水之间存在强烈的相互作用,构建地表水和地下水耦合的水量和溶质运移模型可以揭示地表水和地下水之间的相互作用机制,并提高地表地下水量水质预测的精度。介绍了现有的地表地下耦合的水量及溶质运移模型,对比了各模型的求解方法和适用领域,总结了这些模型的应用现状。其中水量模型中全耦合模型精度较高,是未来模型发展的主流方向,溶质运移模型中溶质反应过程模拟较为欠缺。未来应在优化求解算法、提高溶质反应模拟和加强裂隙介质模拟等方面深入开展相关研究。     

岩体裂隙网络矿物溶解-沉淀-迁移数值模拟

为了研究岩体裂隙网络中矿物溶解-沉淀的变化规律及其对溶质运移的影响,联合矿物溶解-沉淀动力学模型、渗流模型和溶质运移模型建立了岩体裂隙网络矿物溶解-沉淀-迁移模型,经过水溶液络合物计算和矿物成分分析,对研究区域方解石溶解-沉淀作用及溶质运移进行了数值模拟。得出的结论主要有:(1)受裂隙网络分布的影响,溶质运移分布极不均匀,溶质主要通过连通的裂隙从上游向下游迁移;(2)考虑方解石的溶解-沉淀作用之后,在溶解区域溶质浓度增加,在沉淀区域溶质浓度减少;(3)通过对比分析,可知方解石溶解-沉淀受水溶液中CO23-浓度影响较大。     

裂隙岩体溶质运移模型研究

核废料的地下深部贮存、垃圾填埋造成的污水下渗、海水入侵、输油管道老化引起渗漏等诸多应用领域的发展,要求从机理、试验和模型等方面发展裂隙岩体溶质运移理论。基于裂隙岩体系统溶质运移的多尺度概念模型,分析了目前描述裂隙岩体系统中溶质运移各种数学模型的适用性和优缺点,为选取合理的数学模型求解具体的问题提供了参考依据。最后提出了需要进一步研究的问题。     

Experimental investigation of contaminant migration in filled fracture

Experiment equipments involved in the single tube fracture and double tube fracture models are designed to research the characteristics of groundwater flow and solute transport in filled fracture. During the experiment, the state of groundwater flow can be characterized as linear flow, and satisfies Darcy’s law. Therefore, based on the pipe flow of hydraulics and Darcy’s law, the flow rate and water flow velocity can be calculated. Also, dispersion parameters were calculated with the fitting of observed data and analytical solution in the single tube fracture model. Furthermore, effects of some factors on solute transport are involved in the double tube fracture model, and length of branch fracture, particles’ diameter and flow rate in water inlet have been discussed. Results show that the arrival time of concentration peak value in the single tube fracture model is faster than that in the double tube fracture model, and two concentration peak values exist in the double tube fracture model. Arrival time of concentration peak value is faster with the increase of branch fracture length. Furthermore, if the branch fracture is longer, arrival time of the first concentration peak value is faster, while arrival time of the second concentration peak value is slower, relative to short branch fracture.     

裂隙岩体中非饱和渗流与运移的概念模型及数值模拟

探讨了裂隙岩体中非饱和地下水渗流与溶质运移的几种概念模型的构造及数值模拟问题 ,如裂隙网络模型、连续体模型、等效连续体模型、双孔隙度 (单渗透率 )模型、双渗透率模型、多组份连续体模型等。在裂隙岩体中 ,非饱和地下水的渗流可能只局限于岩体中的岩石组份、或裂隙网络 ,也可能在裂隙和岩石中同时发生 ;对前一种情形只需考虑单一连续体中的流动 ,而后一种情况则需要包括地下水在岩石和裂隙之间的交换。岩体中的裂隙网络往往是溶质运移的主要通道 ;但当溶质在裂隙与岩石之间的渗透和扩散是重要的运移机制时 ,就需要考虑岩石与裂隙界面处的溶质交换。为了模拟岩石与裂隙之间地下水和溶质的交换 ,就需要了解岩石与裂隙之间相互作用的模式和范围 ,使得这类问题的概念模型较单一连续体模型多了一层不确定性、其数值模拟也变得更为困难。因为在实际问题中不易、甚至根本不能判别非饱和渗流的实际形态 ,具体采用哪种模型主要取决于分析的目的和对现场数据的掌握程度。不论哪种模型都会受到模型及参数不确定性的影响 ,因此必须考虑与其他辅助模型的比较.     

Fracture step structure: geometrical characterization and effects on fluid flow and breakthrough curve

Groundwater flow and solute transport through fractured rock is highly responsive to the hydraulic anisotropy and heterogeneity that are specific to every major fracture. A major fracture is modeled as the combination of some primal master fractures and several splay fractures that branch out from primal master fractures: step structures (or jog parts). Step structures are commonly observed along a major fracture on various scales. Master fractures were formed and developed by shear movement while some splay fractures were formed by extension normal to their wall. This difference in fracturing process may lead to a permeability difference between master fractures and splay fractures which seems to be one of the major factors controlling flow and solute transport through the fracture networks due to its hydraulic anisotropic and heterogeneous features. This study is composed of two major components: (1) identification and characterization of a step structure from borehole data; (2) evaluation of effect of some idealized step structures on breakthrough curve by numerical simulations. The fracture data of four 1000-m boreholes were used to make clear fracture patterns in the Tono area of Japan. Some major fractures were identified using stereographic projection technique. On the basis of these results, several idealized models of a major fracture having a step was constructed for the numerical study. The obtained results from numerical simulations clearly imply that geometry of step structure plays an important role in flow and transport through the fracture networks.     

Review: The state-of-art of sparse channel models and their applicability to performance assessment of radioactive waste repositories in fractured crystalline formations

Laboratory and field experiments done on fractured rock show that flow and solute transport often occur along flow channels. ‘Sparse channels’ refers to the case where these channels are characterised by flow in long flow paths separated from each other by large spacings relative to the size of flow domain. A literature study is presented that brings together information useful to assess whether a sparse-channel network concept is an appropriate representation of the flow system in tight fractured rock of low transmissivity, such as that around a nuclear waste repository in deep crystalline rocks. A number of observations are made in this review. First, conventional fracture network models may lead to inaccurate results for flow and solute transport in tight fractured rocks. Secondly, a flow dimension of 1, as determined by the analysis of pressure data in well testing, may be indicative of channelised flow, but such interpretation is not unique or definitive. Thirdly, in sparse channels, the percolation may be more influenced by the fracture shape than the fracture size and orientation but further studies are needed. Fourthly, the migration of radionuclides from a waste canister in a repository to the biosphere may be strongly influenced by the type of model used (e.g. discrete fracture network, channel model). Fifthly, the determination of appropriateness of representing an in situ flow system by a sparse-channel network model needs parameters usually neglected in site characterisation, such as the density of channels or fracture intersections.     

Connectivity and single/dual domain transport models: tests on a point-bar/channel aquifer analogue

In porous aquifers, groundwater flow and solute transport strongly depend on the sedimentary facies distribution at fine scale, which determines the heterogeneity of the conductivity field; in particular, connected permeable sediments could form preferential flow paths. Therefore, properly defined statistics, e.g. total and intrinsic facies connectivity, should be correlated with transport features. In order to improve the assessment of the relevance of this relationship, some tests are conducted on two ensembles of equiprobable realizations, obtained with two different geostatistical simulation methods—sequential indicator simulation and multiple point simulation (MPS)—from the same dataset, which refers to an aquifer analogue of sediments deposited in a fluvial point-bar/channel association. The ensembles show different features; simulations with MPS are more structured and characterised by preferential flow paths. This is confirmed by the analysis of transport connectivities and by the interpretation of data from numerical experiments of conservative solute transport with single and dual domain models. The use of two ensembles permits (1) previous results obtained for single realizations to be consolidated on a more firm statistical basis and (2) the application of principal component analysis to assess which quantities are statistically the most relevant for the relationship between connectivity indicators and flow and transport properties.     

反应性溶质在不同质地饱和土柱中运移的数值模拟

盐渍土壤中的物理化学作用对溶质运移具有重要影响.吸附和离子交换作用是土壤中常见的反应.利用室内土柱出流实验对这两种作用下的单组分和多组分溶质运移进行了探讨,用CXTFIT软件模拟了只考虑对流-弥散的常规溶质运移;用水文地球化学模拟软件PHREEQC进行了耦合吸附和离子交换反应的模拟.结果表明,土壤质地对单组分溶质的运移具有重要影响,而在多组分溶质运移中,组分之间的相互作用对溶质运移具有更为重要的影响,并且耦合物理化学作用的模拟精度更高.     

http://dx.doi.org/10.1016/j.gsf.2015.11.001

A new method, the characteristic finite element method(CFEM), was developed to simulate solute transport in a cross-fracture. The solution of this mathematical model for solute transport considered that the contribution of convection and dispersion terms was deduced using the single-step, trace-back method and routine finite element method(FEM). Also, experimental models were designed to verify the reliability and validity of the CFEM. Results showed that experimental data from a single fracture model agreed with numerical simulations obtained from the use of the CFEM. However, routine FEM caused numerical oscillation and dispersion during the calculation of solute concentration. Furthermore, in this cross-fracture model, CFEM simulation results predicted that the arrival time of concentration peak values decreased with increasing flux. Also, the second concentration peak value was obvious with the decrease of flux, which may have resulted from the convergence of solute concentrations from main, and branch, fractures.