黏土岩水力耦合特性试验研究

对于高放废物地质处置工程,地质屏障系统是放射性有害物质进入环境的最后一道屏障,也是一条重要防线。黏土岩由于其低渗透性、渗透损伤自修复以及较强的吸附能力等特性,被认为是一种合理的高放废物处置地质屏障。结合比利时正开展的黏土岩高放废物地质处置相关研究课题,通过一系列室内试验研究黏土岩的水力耦合机制及长期蠕变特性。三轴压缩试验表明,黏土岩压缩强度、超孔隙水压力与围压正相关。渗透试验表明,黏土岩渗透性呈现显著的各向异性特征,围压增大使黏土岩渗透性显著降低。蠕变试验表明,黏土岩蠕变变形、蠕变变形速率与载荷密切相关,即：载荷越大,黏土岩蠕变变形越显著,蠕变变形速率达到稳定所需的时间越长,且相应的稳态蠕变速率越大;根据应力阈值和等时曲线法初步确定该黏土岩长期强度在1.0～1.2 MPa。研究结果将为我国未来黏土岩高放废物处置库的选址和安全性评估提供重要科学依据。     

考虑开挖损伤的高放废物地质处置库温度-渗流-应力耦合数值模拟方法

高放废物地质处置库处于温度-渗流-应力（THM）多场耦合环境中,对高放废物处置库进行安全评估时,需进行多场耦合分析。然而高放废物处置库开挖引起硐壁附近围岩应力重分布,产生损伤,导致围岩热学参数（T）、渗流参数（H）和力学参数（M）发生变化,且在空间上分布不均匀,这将会对运营期处置库THM耦合演化过程产生显著影响。通过分析高放废物处置库温度-渗流-应力三场的耦合原理和处置库围岩损伤的分布和演化规律,定义了损伤变量和损伤演化准则,并将损伤变量与热学参数、渗流参数、力学参数以及多场耦合参数（Biot系数、Biot模量和温度排水系数）建立联系,将围岩损伤与温度-渗流-应力建立联系,形成了一个弹塑性损伤温度-渗流-应力多场耦合数值模型,然后利用建立的模型对瑞士Mont Terri高放废物地质处置库围岩加热试验进行模拟,对比了模拟值和试验值,比较了考虑开挖损伤和不考虑开挖损伤对高放废物地质处置库温度-渗流-应力的影响,并分析了在多场耦合作用下开挖损伤的演化规律。     

中国高放废物地质处置缓冲材料大型试验台架和热-水-力-化学耦合性能研究

缓冲材料作为高放废物深地质处置库中一道重要的人工屏障,与高放废物容器和处置库围岩直接接触,在高放废物衰变热、辐射作用和地下水等影响下产生复杂的热-水-力-化学耦合作用,为了验证缓冲材料是否能长期有效地发挥其屏障材料的作用,核工业北京地质研究院利用高庙子钠基膨润土组装并运行了模拟中国高放废物地质处置室 尺寸的大型缓冲材料膨润土试验台架（China-Mock-Up）。建立了缓冲材料试验台架的安装和试验方法,依据实测数据和理论分析,揭示了热-水-力-化学耦合作用条件下膨润土中的相对湿度是在加热器的热效应和外部供水的湿效应共同作用下发生变化的,压实膨润土中应力的变化主要是由于膨润土遇水膨胀和加热器的热效应引起的,试验验证了模拟高放废物地质处置室内加热器（废物罐）运行初期的位移过程,为缓冲材料和高放废物地质处置库的设计提供了重要的工程参数和理论依据。     

考虑开挖损伤的高放废物地质处置库温度-渗流-应力耦合数值模拟方法

高放废物地质处置库处于温度?渗流?应力（THM）多场耦合环境中,对高放废物处置库进行安全评估时,需进行多场耦合分析。然而,高放废物处置库开挖引起硐壁附近围岩应力重分布,产生损伤,导致围岩热学参数（T）、渗流参数（H）和力学参数（M）发生变化,且在空间上分布不均匀,这将会对运营期处置库THM耦合演化过程产生显著影响。通过分析高放废物处置库温度?渗流?应力三场的耦合原理和处置库围岩损伤的分布和演化规律,定义了损伤变量和损伤演化准则,并将损伤变量与热学参数、渗流参数、力学参数以及多场耦合参数（Biot系数、Biot模量和温度排水系数）建立联系,将围岩损伤与温度?渗流?应力建立联系,形成了一个弹塑性损伤温度?渗流?应力多场耦合数值模型,然后利用建立的模型对瑞士Mont Terri高放废物地质处置库围岩加热试验进行模拟,对比了模拟值和试验值,比较了考虑开挖损伤和不考虑开挖损伤对高放废物地质处置库温度?渗流?应力的影响,并分析了在多场耦合作用下开挖损伤的演化规律。     

黏土岩温度-渗流-应力耦合特性试验与本构模型研究进展

高放废物处置库、垃圾填埋场等工程中常常涉及到温度场（T）、渗流场（H）和应力场（M）的耦合作用的问题。从试验和理论模型两个角度综述国内外黏土岩温度-渗流-应力耦合特性的研究进展,主要包括其传热特性、温度影响下的渗流特性、变形、强度、蠕变特性。在此基础上,重点分析了黏土岩水-热迁移模型以及热-力耦合本构模型的适应性。基于上述认识,通过试验研究了比利时Boom clay在温度作用下的强度、渗透性、蠕变性等特征。结果表明：随着温度升高,Boom clay的强度有所降低,渗透性显著增强,蠕变速率明显加快。提出了适用于Boom clay的THM耦合弹塑性损伤模型,计算结果验证了模型能合理反映温度的影响。最后,探讨了黏土岩THM耦合机理研究的不足和今后的研究方向。     

超固结土热黏弹塑性本构关系

热黏弹塑性本构模型是描述土在温度（热）和时间（黏）耦合作用下的应力-应变关系的本构模型。在一些新型岩土工程诸如高放核废料地质处置、地热资源开发与贮存的建设中,需要同时考虑温度和时间对土的影响,所以建立一个热黏弹塑性本构模型具有理论和实际意义。将温度变化对黏土体积和强度参数的影响引入笔者之前提出的超固结土等向应力-应变-时间关系,建立了一个等向应力条件下的应力-应变-时间-温度关系。随后,基于该关系推导了屈服面硬化定律,并将其与超固结土统一硬化模型的屈服方程和流动法则结合,建立了超固结土的热黏弹塑性本构模型。最后,使用新模型预测室内试验,证明新模型能够反映时间和温度对土体积、一维压缩曲线和前期固结压力的耦合影响。     

考虑细观等效热学参数的高放废物处置库围岩应力−渗流−温度耦合模拟方法

高放废物处置库中乏燃料持续释放的热量对围岩的应力场和渗流场及其长期稳定性具有重要影响。围岩的热学参数依赖于岩石矿物组成、孔隙率和孔隙流体等因素,准确取值是进行高放废物地质处置库多场耦合分析的前提。通过细观力学分析,建立了围岩等效热学参数（热容、热传导系数、热膨胀系数）取值方法,并基于Biot孔隙介质理论,建立应力?温度?渗流三场耦合模型,进而提出了高放废物处置库围岩应力?渗流?温度耦合数值模拟方法。最后通过COMSOL Multiphysics多场耦合软件,利用瑞士Mont Terri高放废物地下试验室围岩温度?渗流?应力多场耦合现场试验数据对数值模拟方法进行验证,并探讨了温度?渗流?应力耦合过程的演化规律。研究表明,模拟结果和试验值吻合良好。研究结果可为我国高放废物处置库的安全评估和选址提供科学依据。     

高放废物地质处置黏土岩处置库围岩研究现状

世界上很多国家都对处置库的可能围岩进行了详细研究。通过对比,认为花岗岩、黏土岩、岩盐比较适合作为处置库围岩,而黏土岩由于具有自封闭性、渗透率低等其他岩石类型不可比拟的优点,因而将黏土岩作为高放废物地质处置库围岩越来越受到各国的关注。文章同时介绍了瑞士、法国、比利时等国家在黏土岩中所进行的大量研究,均认为在黏土岩中处置高放废物和乏燃料是安全的。文章还对黏土岩处置库概念设计、黏土岩处置库围岩地下实验室研究,以及我国开展黏土岩处置库研究的意义等进行了综述。     

高放废物地质处置数据资源集成开发进展

数据资源是高放废物地质处置库选址、安全评价、工程建设等研发的重要基础,该领域数据具有典型的大数据特性,其开发利用是一项长期、系统的科学工程。文章对我国高放废物地质处置数据资源集成开发策略、方法和现状进行了简要介绍,阐述了高放废物地质处置研发过程中产生的各类科学数据的具体特性,提出建立我国高放废物地质处置大数据平台的规划设想,重点介绍了处置库预选区地学信息集成开发现状。     

一个描述软黏土时效特性的一维弹黏塑性模型

为了描述软黏土一维应力-应变关系的时效特性,基于Bjerrum的等时间线体系,提出等黏塑性应变率线概念,推导了黏塑性应变率与黏塑性应变增量的关系,建立了软黏土的一维弹黏塑性模型;从理论上分析了新建模型与3种典型的一维弹黏塑性时效本构模型的内在联系,表明新建模型与其他3种模型在本质上是等效的,且形式更简洁,物理意义更明确;利用新建模型对软黏土的固结-蠕变耦合效应、应变率效应、应力松弛效应等时效特性进行了理论分析,并得到了相应的解析解;结合宁波软黏土的一维固结试验,阐述了模型参数的确定方法,并用新建模型对宁波软黏土的固结-蠕变试验、温州软黏土的一维多级等应变率试验、香港海相软黏土的一维应力松弛试验进行模拟,验证了新建模型的有效性。研究结果表明,新建模型能很好地模拟软黏土的一维时效特性。     

黏土岩饱和过程中水分运移规律试验研究

黏土岩作为高放废物地质处置库的备选介质,目前得到世界各国的高度重视。黏土岩地质处置库巷道施工过程中,一方面,围岩因开挖损伤生成裂隙使得渗透性增强,对核素的阻滞作用降低;另一方面,在应力和水的耦合作用下,黏土岩良好的裂隙渗透损伤自修复能力使得围岩的渗透性逐渐恢复接近于原始状态。基于电阻率测试,首先开展了黏土岩试样在不同条件下的饱和过程试验研究,得到了黏土岩试样饱和过程中等效电阻率的变化规律,分析了不同损伤程度试样、盐溶液对等效电阻率的影响,进而揭示黏土岩饱和过程中水分运移规律。试验结果表明：（1）等效电阻率随着含水率增加而逐渐减小,并逐渐趋于一个稳定值;（2）等效电阻率的大小不仅与含水率有关,试样内部裂隙的存在也会影响等效电阻率分布,这一发现为电阻率法可以探测试样中裂隙的存在提供了依据;（3）水流在黏土岩中扩散,内部裂隙成为优先通道,水流在裂隙中的快速扩散加快了黏土岩的饱和速度。同时,随着黏土岩中水分与黏土矿物的水化膨胀反应,内部裂隙有一定程度闭合,加深对裂隙闭合机制认识,通过电阻率测试可以有效地揭示这一过程。     

热−水−力耦合条件下黏土岩蠕变特性研究

黏土岩作为放射性核废料处置的备选介质,长期处于热-水-力耦合的复杂条件下。为研究围岩的长期稳定性,开展了一系列不同围压和偏应力下的黏土岩加温-降温排水蠕变试验。得出以下结论：温度的升高会提升黏土岩的蠕变速率,延长衰减蠕变阶段的时间,但是降温过程中试样主要以冷缩变形为主,并未见明显的蠕变变形;围压的降低以及偏应力的增加会提高黏土岩的蠕变速率,并且这种影响会随温度的升高明显加剧。基于试验结果,在Perzyna过应力理论的基础上引入蠕变硬化、蠕变损伤以及热损伤,建立了黏土岩热-水-力耦合蠕变模型。通过ABAQUS软件及其子程序对该模型进行数值实现,模拟结果与试验结果的对比初步表明该模型可以很好地描述热-水-力耦合条件下黏土岩的蠕变特性。     

The influence of an excavation damaged zone on the thermal-mechanical and hydro-mechanical behaviors of an underground excavation

In Korea, a reference disposal system, KRS, was proposed in 2006 after 10 years of research and development. In the KRS, the high-level radioactive waste repository is considered to be located in a crystalline rock likes granite. For a validation of the feasibility, safety, and stability of the KRS, an underground research tunnel, KURT was constructed in Nov. 2006. During the construction of KURT by a controlled blasting, the size and characteristics of an excavation damaged zone(EDZ) were investigated by in situ as well as laboratory tests. The possible influences of an EDZ around a tunnel on the thermal, hydraulic and mechanical behaviors of the near field were investigated by using hydro-mechanical and thermo-mechanical coupling analyses. From this study, it was found that the existence of an EDZ can influence the thermal, hydraulic, and mechanical behaviors of the near field and it was recommended that an EDZ should be considered as an important parameter during the design of underground repositories.     

Laboratory Investigations of the Hydro-Mechanical–Chemical Coupling Behaviour of Sandstone in CO2 Storage in Aquifers

This paper is devoted to experimental investigations of the hydro-mechanical–chemical coupling behaviour of sandstone in the context of CO₂ storage in aquifers. We focused on the evolution of creep strain, the transport properties and the elastic modulus of sandstone under the effect of CO₂–brine or CO₂ alone. A summary of previous laboratory results is first presented, including mechanical, poromechanical and hydro-mechanical–chemical coupling properties. Tests were then performed to investigate the evolution of the creep strain and permeability during the injection of CO₂–brine or CO₂ alone. After the injection of CO₂–brine or CO₂ alone, an instantaneous volumetric dilatancy was observed due to the decrease in the effective confining stress. However, CO₂ alone had a significant influence on the creep strain and permeability compared to the small influence of CO₂–brine. This phenomenon can be attributed to the acceleration of the CO₂–brine–rock reaction by the generation of carbonic acid induced by the dissolution of CO₂ into the brine. The original indentation tests on samples after the CO₂–brine–rock reaction were also performed and indicated that the elastic modulus decreased with an increasing reaction time. The present laboratory results can advance our knowledge of the hydro-mechanical–chemical coupling behaviour of sandstone in CO₂ storage in aquifers.     

Fully coupled hydro-mechanical numerical manifold modeling of porous rock with dominant fractures

Coupled hydro-mechanical (HM) processes are significant in geological engineering such as oil and gas extraction, geothermal energy, nuclear waste disposal and for the safety assessment of dam foundations and rock slopes, where the geological media usually consist of fractured rock masses. In this study, we developed a model for the analysis of coupled hydro-mechanical processes in porous rock containing dominant fractures, by using the numerical manifold method (NMM). In the current model, the fractures are regarded as different material domains from surrounding rock, i.e., finite-thickness fracture zones as porous media. Compared with the rock matrix, these fractured porous media are characterized with nonlinear behavior of hydraulic and mechanical properties, involving not only direct (poroelastic) coupling but also indirect (property change) coupling. By combining the potential energy associated with mechanical responses, fluid flow and solid–fluid interactions, a new formulation for direct HM coupling in porous media is established. For indirect coupling associated with fracture opening/closure, we developed a new approach implicitly considering the nonlinear properties by directly assembling the corresponding strain energy. Compared with traditional methods with approximation of the nonlinear constitutive equations, this new formulation achieves a more accurate representation of the nonlinear behavior. We implemented the new model for coupled HM analysis in NMM, which has fixed mathematical grid and accurate integration, and developed a new computer code. We tested the code for direct coupling on two classical poroelastic problems with coarse mesh and compared the results with the analytical solutions, achieving excellent agreement, respectively. Finally, we tested for indirect coupling on models with a single dominant fracture and obtained reasonable results. The current poroelastic NNM model with a continuous finite-thickness fracture zone will be further developed considering thin fractures in a discontinuous approach for a comprehensive model for HM analysis in fractured porous rock masses.     

The use of two-part Hooke’s model (TPHM) to model the mine-by test at Mont Terri Site,Switzerland

The full-scale mine-by (MB) test conducted in 2008 in the Mont Terri underground rock laboratory (Switzerland) investigated the deformation and the coupled hydro-mechanical behavior of the Opalinus Clay in response to tunnel excavations. The Opalinus Clay is currently under investigation in Switzerland as a potential host rock for geologic disposal of high-level radioactive waste. To further improve the understanding and modeling of the coupled processes and their impact on the performance of a geologic repository in Opalinus Clay, a newly developed two-part Hooke’s model (TPHM) was implemented into a geomechanical simulator. A three-dimensional simulation model based on the TPHM was then developed to predict the deformation and pore pressure responses in the near field of the MB Niche 2 test at the Mont Terri Site. The usefulness and validity of the TPHM are demonstrated by the consistency between simulation results and field observations. Simulation results show that the pore pressure disturbance becomes visible at about 11 m ahead of the mine-by excavation advancing face (along the longitudinal direction of the MB Niche). The results also demonstrate that there exists a good correlation between the excavation damage zone (EDZ) and the pore-pressure evolution, which may have important practical implications for monitoring EDZ evolution with pore-pressure sensors. The simulation results, which are sensitive to the constitutive relationships used in the model, capture both the observed displacements and the size of the damage zone, whereas the approach based on the conventional Hooke’s law underestimates both. The comparison between simulated and observed results also indicates that laboratory-measured mechanical properties can be used to accurately predict field-scale mechanical deformations, as long as valid constitutive relationships are employed.     

Study of the Anisotropic Properties of Argillite Under Moisture and Mechanical Loads

Due to various factors, such as sedimentation, layered morphology of clay minerals, in situ stress, etc., argillite rocks often exhibit anisotropic behavior. In order to study the anisotropic properties of the Callovo-Oxfordian (COx) argillite of the Meuse–Haute-Marne site in France considered as a possible host rock for high-level radioactive nuclear waste repository, a series of tests including uniaxial compression and dehydration and hydration at different constant applied stress levels are carried out. In this study, a specific setup combining moisture and mechanical loading with optical observation is used and it allows to continuously capture surface images from which the full-field strains are determined by using Digital Image Correlation techniques. The results show evidence of the mechanical and hydric anisotropy of the material. The anisotropy parameters are identified, assuming the studied argillite as transversely isotropic. The shrinkage and swelling depend on the applied stress and the angle with respect to the vertical direction of the mechanical load and the stratification plane, and this dependence is quantified. The non-linearity and the hysteresis observed during dehydration and hydration cycles are discussed.