高压实膨润土的非饱和渗透膨胀模型

高放废弃物深地质处置库中,作为缓冲/回填材料的高压实膨润土的设计功能是,遇水吸湿导致土体膨胀变形,以密封高压实膨润土块体砌置时形成的块体与块体之间和块体与围岩之间的施工缝隙以及围岩中因处置库开挖卸载引起的裂缝,形成阻障围岩中地下水渗入内库引起核素迁移,与库内高放废物的辐射扩散的人工屏障。高压实膨润土在自由膨胀条件下的膨胀应变和非饱和渗透系数由渗透试验测得;孔隙比和吸力的关系由压汞试验测得。吸力和含水率的函数关系可由土-水特征曲线测试试验得到。因此,自由膨胀条件下膨润土的非饱和渗透膨胀模型由膨胀应变、非饱和渗透系数和孔隙比与吸力的函数关系建立。结果表明：渗透系数和体积应变都随吸力减小而增加,因为吸力减小会导致孔隙比增加。这一发现有助于设计高放废物的最终处置方案。     

高放废物地质处置系统核素迁移模型研究

高放废物地质处置系统核素迁移研究是高放废物安全处理和处置的重要内容。从整个系统的角度出发,把高放废物地质处置系统核素迁移过程分为三种模式:工程屏障中的释放—扩散模型;地质屏障中的对流—弥散模型和在生物圈中的分区传递模型。在分区传递模型中利用转移系数描述核素在各分区之间的迁移。通过对每种模型核素迁移机理的分析得出概化模型,并给出了相应的数学描述。     

中国高放废物处置库选址中灵敏地质特征的初步确认

本文首先根据中国、瑞士和日本3国高放废物处置概念的相似性，建立中国第一处置库系统性能评价模型；然后通过对9个地质变量在8个事件中的灵敏度分析，初步确认出3个灵敏地质物征。此项研究对中国目前高放废物处置库选址工作有一定的指导意义。     

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

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

The effect of heterogeneity of diffusion parameters on chloride transport in low-permeability argillites

Understanding flow and transport in low-permeability media is very important in the context of nuclear waste disposal, oil and gas reservoirs and long term evolution of groundwater systems. In low-permeability media, transport by diffusion is often the most important mass transport process. This study investigates the effect of the heterogeneity of diffusion parameters on mass transport in low-permeability media. A geostatistical approach for integrating heterogeneity of diffusion parameters in groundwater flow and transport models is proposed and applied to the Toarcian argillites in France which are studied in the framework of feasibility of storing radioactive waste in deep clayey massifs. Stochastic fields of the diffusion parameters of the Toarcian argillites (France) are generated based on 64 measured values of diffusion coefficient and diffusion accessible porosity and used as input for a 3D local-scale groundwater flow and transport model. The chloride concentrations computed by these heterogeneous models are compared to the measured chloride concentrations and to concentrations calculated with a model in which the Toarcian argillites are subdivided into several homogeneous zones. The heterogeneous simulations result in a slightly better correspondence between measured and calculated values and have the additional advantage that the measured diffusion coefficient values in the Toarcian are perfectly honored in the model. This study shows that small-scale variability of diffusion parameters has a significant effect on solute concentrations and omitting this heterogeneity may be a problem in transport calculations in low-permeability media, depending on the specific setting and objectives of the study.     

System reliability analysis for near-surface radioactive waste disposal facilities

The low- and intermediate-level radioactive wastes are generally disposed in near-surface disposal facility (NSDF). The NSDF is composed of engineered barriers. The probabilistic safety assessment model has been developed to analyse the performance of NSDF. The endpoints of the assessment are the concentration of the radionuclide in the groundwater and the corresponding dose rate of the radionuclide. The barrier system can have multiple failure modes but practically the possible failure modes are failure of top cover, failure of waste container, degradation of waste form and failure of bottom cover, which are usually considered as independent failure events. Through a sensitivity analysis, the most critical parameters affecting the design reliability for failure criteria are identified as the groundwater velocity and distribution coefficient. The study shows that for the NSDF considered, there is a high degree of dependence between the failure modes, and demonstrates the probability of the simultaneous occurrence of failures. Thus, the need to consider the system reliability in the NSDF is highlighted. The study also advocates the use of optimisation techniques to evaluate the probability of failure, which provides a better estimate of the probability of failure, as validated from the results obtained from Monte Carlo simulations.     

Stochastic analysis of the effect of heterogeneity and fractures on radionuclide transport in a low-permeability clay layer

Deep low-permeability clay layers are considered as safe environments for disposal of high-level radioactive waste. In Belgium, the Boom Clay is a candidate host rock for deep geological disposal. In this study, we analyze the effects of fractures and spatially variable hydraulic conductivity on radionuclide migration through the clay. Fracture geometry and properties are simulated with Monte Carlo simulation. The heterogeneity of hydraulic conductivity is simulated by direct sequential co-simulation using measurements of hydraulic conductivity and four types of secondary variables. The hydraulic conductivity and fracture simulations are used as input for a transport model. Radionuclide fluxes computed with this heterogeneous model are compared with fluxes obtained with a homogeneous model. The output fluxes of the heterogeneous model differ at most 8% from the homogeneous model. The main safety function of the Boom Clay is thus not affected by the fractures and the spatial variability of hydraulic conductivity.     

西北某干旱区多级储水洼地地下水数值模拟

西北某放射性废物处置场预选区,区域地下水系统包含多级独立第四系储水洼地,洼地出口以泉水排泄地下水,继而回渗补给下级洼地。根据含水层底板起伏特征,利用GMS模拟软件中的排水沟模块与溪流模块概化了泉水,通过泉流量校准与监测孔水位拟合,校正了研究区渗透系数、给水度、和储水率等水文地质参数,计算了地下水流速、流向以及地下水资源量。结果显示,研究区地下水由南向北东径流,东北部为最终排泄洼地,地下水流速缓慢,水资源相对匮乏,有利于放射性废物处置的安全。     

黏土岩时效变形特性试验与理论研究

高放废物地质处置岩体材料选取中,黏土岩因具有低渗透性、损伤自修复特性、对放射性核素具有良好的吸附作用等优点,被认为是一种合理的高放废物地质处置屏障。以黏土岩为研究对象,从黏土岩短期、长期力学特性等方面开展研究工作,主要内容如下：（1）进行黏土岩短期自然固结试验,确定了黏土岩的基本物理力学参数,主要包括前期固结压力、压缩指数等;（2）进行流-固耦合固结试验,试验结果表明盐水作用对黏土岩力学特性具有重要影响,且流-固耦合过程中的膨胀现象与黏土岩的黏土矿物含量和类型密切相关;（3）通过对黏土岩进行固结流变试验,研究黏土岩渗流-应力耦合作用下的长期力学特性,试验表明黏土岩具有明显的流变特性,且流变现象与载荷密切相关;（4）根据固结流变试验建立黏土岩一维流变本构模型,同时将模型计算结果与试验结果进行对比,分析表明,该模型能够很好地反映黏土岩一维固结流变特性。该研究对我国未来黏土岩高放废物处置库的规划、设计、选址和运营等具有重要的参考意义。     

核废物地质处置模型试验近场热-水-应力耦合数值模拟

笔者描述了所建立的饱和-非饱和介质中热-水-应力耦合现象的控制方程。针对国际合作项目DECOVALEX II 中Task-2的核废物地质处置轴对称模型试验及计算参数,使用所研制的有限元程序进行了热-水-应力耦合过程模拟。对缓冲层及岩体中若干部位的温度、含水量及接触应力的分布、变化进行了分析, 并将其与量测值和使用THAMES程序的计算结果作了比较,看到了三者的异、同之处,并指出了可能的影响原因。通过与量测值和THAMES程序的计算结果的对比,显示了笔者所研制的有限元程序在热-水-应力耦合数值分析中是可靠和实用的。     

高放废物深地质处置场近场剂量场-温度场-水势场-应力场的耦合计算

回顾了高放废物处置概况, 给出了剂量场-温度场-水势场-应力场耦合研究的意义, 探讨了剂量场计算公式与热-湿-力耦合方程及其求解, 最后指出了这一领域的若干优先研究方向.      

Application of uncertainty analysis to groundwater pollution modeling

Prediction and evaluation of pollution of the subsurface environment and planning remedial actions at existing sites may be useful for siting and designing new land-based waste treatment or disposal facilities. Most models used to make such predictions assume that the system behaves deterministically. A variety of factors, however, introduce uncertainty into the model predictions. The factors include model and pollution transport parameters and geometric uncertainty. The Monte Carlo technique is applied to evaluate the uncertainty, as illustrated by applying three analytical groundwater pollution transport models. The uncertainty analysis provides estimates of statistical reliability in model outputs of pollution concentration and arrival time. Examples are provided that demonstrate: (a) confidence limits around predicted values of concentration and arrival time can be obtained, (b) the selection of probability distributions for input parameters affects the output variables, and (c) the probability distribution of the output variables can be different from that of the input variables, even when all input parameters have the same probability distribution     

基于BQ系统的阿拉善巴彦诺日公NRG01号钻孔岩体质量评价

高放废物深部地质处置目前受到世界各国的高度重视,处置库围岩岩体质量评价是高放废物处置库选址的关键问题之一。文章基于BQ系统对我国高放废物地质处置库阿拉善预选区巴彦诺日公NRG01号钻孔进行了岩体质量评价。NRG01号孔钻孔资料丰富,采用岩芯编录、波速测井、水压致裂地应力测量及岩石波速测试、单轴压缩试验等多种方法手段相结合的方式获取了BQ评价系统中的各个参数。评价过程中,综合考虑了岩体坚硬程度、完整程度及地下水、软弱结构面、初始地应力等因素对围岩稳定的影响。根据计算的[BQ]值对岩体质量进行评价并与RQD评价结果进行对比。评价结果表明,NRG01号孔岩体质量较好,Ⅰ级和Ⅱ级岩体占90%,410~500 m的岩体范围可作为处置库建设的目标岩体。[BQ]值与RQD值所得的岩体质量评价结果一致性良好,但由于BQ系统考虑的影响因素更为全面,评价结果更为精细。本次评价可为巴彦诺日公地段高放废物地质处置库场址比选及处置库目标岩体选择提供依据。     

孔隙气压对核废物地质处置热-流-固耦合影响的有限元分析

为考察孔隙气体压力对高放废物地质处置中的热-流-固耦合过程的影响,借用Leiws等建立的可变形孔隙介质中非等温空气流和水流模型,在其水连续性方程中加入了温度梯度引起的水分扩散项,研制出相应的热-水-气-应力耦合弹塑性二维有限元程序。针对一个假定的高放废物地质处置库模型,在相同的初始温度、孔隙水压力和岩体应力条件下,取3种缓冲层中的初始孔隙气体压力,通过数值模拟考察了处置库近场的主应力、饱和度、气和水的流速、温度和孔隙压力的分布与变化。结果显示,当缓冲层中的初始孔隙气压力较高时,其对围岩中应力影响较大。     

高放废物处置库北山预选区区域水文地质概念模型

地下水流数值模拟是高放废物处置库选址中水文地质评价的重要内容.水文地质概念模型是地下水流数值模拟的基础和前提.在北山预选区水文地质分析基础上,利用DEM数据生成了研究区地表水系,据此确定了模拟范围;借助GMS软件,建立了研究区三维水文地质结构模型;再通过边界条件、流场特征、地下水均衡项和水文地质参数的分析和描述,建立了北山地区水文地质概念模型.     

Hydrogeological modeling of radionuclide transport in low permeability media: a comparison between Boom Clay and Ypresian Clay

Deep low-permeability clay layers are considered as suitable environments for disposal of high-level radioactive waste. In Belgium, the Boom Clay is the reference host formation and the Ypresian Clay an alternative host formation for research and safety and feasibility assessment of deep disposal of nuclear waste. In this study, two hydrogeological models are built to calculate the radionuclide fluxes that would migrate from a potential repository through these two clay formations. Transport parameter heterogeneity is incorporated in the models using geostatistical co-simulations of hydraulic conductivity, diffusion coefficient and diffusion accessible porosity. The calculated radionuclide fluxes in the two clay formations are compared. The results show that in the Ypresian Clay larger differences between the fluxes through the lower and the upper clay boundary occur, larger total output radionuclide amounts are calculated and a larger effect of parameter heterogeneity on the calculated fluxes is observed, compared to the Boom Clay.     

中国高放废物地质处置:现状和展望

本文阐述了中国的高放废物处置政策,回顾了中国高放废物地质处置在选址和场址评价、缓冲回填材料等方面的进展,提出了中国高放废物地质处置的初步规划,指出我国计划在2015年确定地下实验室场址.还展望了国际合作前景.     

基于离散裂隙网络模型的核素粒子迁移数值模拟研究

高放废物地质处置过程中涉及的核素在围岩裂隙地下水中的迁移问题已引起广泛关注,数值模拟是研究核素粒子运移的重要方法。目前裂隙介质中渗流模型主要是等效连续介质模型、双重介质模型和离散裂隙网络模型。对于岩体尺度裂隙地下水的流动,离散裂隙网络模型能充分表现裂隙介质的各向异性、不连续性等特征。因此,针对裂隙介质准确概化及核素迁移模拟等难点,文章结合Monte Carlo随机生成裂隙方法、裂隙渗流有限元算法和高放射性核素衰变方程等方法,依据花岗岩深钻孔裂隙统计数据,采用离散裂隙网络模型对内蒙古阿拉善高放废物地质处置预选区展开了核素粒子迁移数值模拟研究,并讨论了实例预测分析结果。结果显示:针对设定的地质模型,核素粒子从中心运移到边界的迁移路径长度平均为1293.35 m,粒子运移到边界耗费的时间平均为1.70E+11 d。     

考虑膨胀效应UH模型的有限元实现

高庙子膨润土已被确定为我国高放废物地质处置的首选缓冲/回填材料,进行其力学特性及应用模拟研究对合理设计处置库具有十分重要的意义。考虑膨胀效应的UH模型是在超固结非饱和土本构模型（非饱和土UH模型）的基础上提出的,因引入了土体团粒的吸水膨胀作用,使其适用于非饱和膨润土。基于有限元软件ABAQUS的二次开发平台,利用Fortran语言编制模型的UMAT（user-difined material）材料子程序,进行了考虑膨胀效应UH模型的有限元实现。然后针对高放废物处置概念模型,利用编制的模型子程序对高放废物处置系统进行了三维有限元模拟,并对渗流场和应力场的演化规律进行了初步探究。模拟结果显示处置系统中膨润土的内应力增加,土体的稳定性较好,说明选择高庙子膨润土作为缓冲/回填材料的合理性,也同时验证了程序的正确性。     

Study on the residence time of deep groundwater for high-level radioactive waste geological disposal

Residence time of deep groundwater is one of the most important parameters in safety and performance assessment for high-level radioactive waste geological disposal. In this study, we collected the deep groundwater samples of Jijicao in Gansu Beishan pre-selected region. The deep groundwater residence time at two depths estimated by Helium-4 accumulation method were 3.8 ka and 5.0 ka respectively upon measurement and calculation, which indicates that the deep groundwater is not derived from the deep crust circulation process. Hence, deep groundwater is featured with long residence time as well as slow circulation and update rate, and such features are conductive to the safe disposal of high-level radioactive waste.