• 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊
  • Scopus数据库收录期刊

双孔隙压实膨润土的细观导热模型

周殷康, 阎长虹, 郑军, 谢胜华, 项国圣

周殷康, 阎长虹, 郑军, 谢胜华, 项国圣. 双孔隙压实膨润土的细观导热模型[J]. 岩土工程学报, 2021, 43(7): 1352-1359. DOI: 10.11779/CJGE202107022
引用本文: 周殷康, 阎长虹, 郑军, 谢胜华, 项国圣. 双孔隙压实膨润土的细观导热模型[J]. 岩土工程学报, 2021, 43(7): 1352-1359. DOI: 10.11779/CJGE202107022
ZHOU Yin-kang, YAN Chang-hong, ZHENG Jun, XIE Sheng-hua, XIANG Guo-sheng. Mesoscale model for thermal conductivity of compacted dual-porosity bentonite[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(7): 1352-1359. DOI: 10.11779/CJGE202107022
Citation: ZHOU Yin-kang, YAN Chang-hong, ZHENG Jun, XIE Sheng-hua, XIANG Guo-sheng. Mesoscale model for thermal conductivity of compacted dual-porosity bentonite[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(7): 1352-1359. DOI: 10.11779/CJGE202107022

双孔隙压实膨润土的细观导热模型  English Version

基金项目: 

江苏省重点研发计划项目 BE2019705

详细信息
    作者简介:

    周殷康(1988— ),男,博士,安徽太湖人,主要从事岩土体传热相关的试验和模拟工作。E-mail:Zhouyk@ahut.edu.cn

    通讯作者:

    阎长虹, E-mail:yanchh@nju.edu.cn

  • 中图分类号: TU45

Mesoscale model for thermal conductivity of compacted dual-porosity bentonite

  • 摘要: 压实膨润土导热系数是评价核废料处置库散热安全以及计算温度场分布的关键指标,准确预测不同条件下压实膨润土导热系数对于处置库的安全设计具有重要意义。数值模拟是预测土体导热系数的有效途径之一,然而,由于压实膨润土双孔隙结构的复杂性,基于细观数值模拟的压实膨润土导热模型较少。建立了一种适用于双孔隙压实膨润土导热系数的细观模型,选取了不同含水率和孔隙比下5种典型压实膨润土160个实测导热数据,并与模型预测导热系数值进行比较。结果显示膨润土导热系数的模型预测值基本在实测值±20%范围内,有效验证了模型的可靠性,表明本文模型在压实膨润土导热系数预测方面的潜力。为复杂双孔隙压实膨润土导热系数的精确预测提供了一种新的思路和方法。
    Abstract: The thermal conductivity of compacted bentonite is very crucial in evaluating the safety of heat dissipation and in calculating the temperature distribution of radioactive waste repository. An accurate predictive model is of great significance to the calculation of their thermal conductivity under various conditions. A mesoscale model is proposed for the thermal conductivity of compacted dual-porosity bentonite by numerical simulation. The thermal conductivity of 160 measured samples is selected and compared with their simulated one for five typical kinds of compacted bentonite with different dry densities and water contents. The results show that the most predictive values are within ±20% range of the measured ones, which verifies the reliability of the model effectively. The results in this study provide a new approach and/or strategy to predict the thermal conductivity of the compacted dual-porosity bentonite.
  • 在许多国家的高放废物(HLW)深地质处置库概念中,高压实膨润土块体等缓冲/回填材料被用来填充处置容器与围岩之间的空隙,构成一道重要的工程屏障,以延缓地下水入渗、阻滞核素迁移、传递核素衰变热和维护处置库结构稳定性[1]。然而,以高压实膨润土块体砌筑的工程屏障将不可避免地存在各种各样的施工接缝[2],包括块体与处置容器之间、块体与块体之间及块体与围岩之间的缝隙。这些施工接缝将成为地下水入渗和核素外泄的优势通道,大大降低了工程屏障的缓冲性能,直接威胁到处置库的长期安全运营[3]。为此,Salo等[4]首次提出采用高密度膨润土颗粒及其与膨润土粉末的混合物来填充施工接缝(图1(a)),或者取代膨润土块体作为主体缓冲/回填材料(图1(b))。

    图  1  膨润土颗粒在高放废物处置库中的应用
    Figure  1.  The application of bentonite pellets in HLW repository

    膨润土颗粒是一种由天然膨润土粉末制成的高干密度(通常大于1.80 g/cm3)颗粒状土材料。由于其具有高膨胀性、低渗透性和良好的核素吸附性能,又便于机械化制造、袋装化运输和自动化填充,且能够有效填充处置库中的各种施工接缝或空隙,膨润土颗粒已被瑞典SKB、芬兰Posiva、加拿大NWMO、瑞士NAGRA、法国ANDRA、比利时SCK/CEN及其它国家(如日本、中国、印度、韩国和乌克兰等)的处置库概念选为候选缓冲/回填材料[5]。许多学者也针对膨润土颗粒及其混合物的堆积性质、热传导性质、水力–力学特性和结构演化规律开展了大量研究工作[6]。遗憾的是,这些研究成果大多零散地分布在SKB、Posiva和NAGRA等机构的研究报告中,缺乏必要的整理与综述。

    本文拟对膨润土颗粒的制备方法、填充技术与堆积性质、热传导特性、水力特性、结构演化规律及力学特性等6个方面的研究成果和最新进展进行简要回顾与介绍,并分别指出各方面值得进一步深入研究的几个课题,以期为膨润土颗粒相关课题研究和工程实践提供参考与借鉴。

    按照制备工艺原理,膨润土颗粒的制备方法主要可分为5种:挤压法、辊压法、压实法、压实–破碎法和湿–干–破碎法。其中,挤压法和辊压法在国际上应用最为普遍,压实–破碎法在国内应用最为广泛。各种方法的原理和优缺点如表1所示。

    表  1  膨润土颗粒的制备方法比较
    Table  1.  Overview of manufacture techniques for bentonite pellets
    方法原理颗粒性质优点缺点
    挤压法机械压实形状规则大小统一工序简单,效率高,机械化程度高需要特定的制样模具,粒径范围有限
    辊压法
    压实法
    压实–破碎法机械压实机械破碎形状各异大小不一可制备各种粒径的颗粒工序多,效率低,难以控制机械破碎的初始粒径
    湿–干–破碎法吸力固结机械破碎
    下载: 导出CSV 
    | 显示表格

    挤压法采用一对碾轮将膨润土粉末压入多孔不锈钢碾盘中,成型的膨润土“条”从碾盘下部挤出,可根据需要截取特定长度的膨润土颗粒,颗粒呈杆状或圆柱状(图2[7-8]

    图  2  挤压法及其制备的颗粒[7-8]
    Figure  2.  Extrusion method and produced pellets[7-8]

    辊压法采用一组相向滚动的辊轮将膨润土粉末压入辊轮表面的凹“坑”中,形成枕状或杏仁状的颗粒(图3[7,9-11]

    图  3  辊压法及其制备的颗粒[7,11]
    Figure  3.  Roller compaction method and produced pellets[7,11]

    压实法是直接将膨润土粉末快速压入特定规格的模具中来制备高密度颗粒,颗粒呈规则的凸顶圆柱体状(图4(a)[12]

    图  4  压实法与压实–破碎法制备的颗粒[12-13]
    Figure  4.  Pellets manufactured by compaction and compaction-crushing methods[12-13]

    压实–破碎法先采用压力机将膨润土粉末压实成高密度块体,然后采用破碎机将块体破碎成不同粒径的颗粒,最后筛分成不同粒组以配制不同粒径级配的颗粒混合物(图4(b))。文献[6,13,1419]分别采用该方法制备了不同粒径的高庙子(GMZ)膨润土颗粒。

    湿–干–破碎法是先将膨润土粉末湿化成膏状,待其风干固结后,用破碎机械将其破碎成细小的颗粒,最后筛分成不同粒组以配制不同粒径级配的颗粒混合物[20]

    实际上,各种方法的颗粒产品质量(干密度、强度及微裂缝发育程度等)主要取决于膨润土粉末的矿物成分、初始含水率以及制样模具与制造工艺的精细化程度,而与制造原理本身关系不大。各种颗粒制备方法各有利弊,应根据工程实际需要来加以选择。条件允许时,可以联合采用两种或多种方法以提高颗粒制备效率和达到特定的颗粒规格要求。

    在地下原位填充试验和一些室内大型模拟填充试验中,常采用机械设备将颗粒混合物填充到膨润土块体与围岩的缝隙中,或者用于填充整个巷道。按照机械设备的工作原理,颗粒填充技术主要包括带式输送法、螺旋输送法、气动喷射法和人工填充法(图5),其中螺旋输送法和气动喷射法应用最为广泛[21-22]。各种方法的优缺点和适用性如表2所示。

    图  5  带式输送法、螺旋输送法和气动喷射法示意图
    Figure  5.  Schematic diagrams of pellets emplaced by belt, auger and pneumatic conveying
    表  2  膨润土颗粒的填充方法比较
    Table  2.  Overview of filling techniques for bentonite pellets
    方法主要优点主要缺点适用性
    带式输送效率高粉尘多
    易离析
    设备笨重
    填充较大的空隙或回填巷道
    螺旋输送效率较高
    粉尘较少
    均匀性较好
    设备较笨重填充较大的接缝、空隙或回填巷道
    气动喷射效率较高
    堆积密度大
    粉尘多
    均匀性差
    可填充各种接缝、空隙或回填巷道
    人工填充适应性强粉尘多
    效率低
    可填充各种接缝、空隙或回填巷道
    下载: 导出CSV 
    | 显示表格

    各种颗粒填充技术各具特色和局限性,应根据实际情况加以选择利用。评价颗粒填充效果需关注3项指标:堆积干密度、堆积均匀性和粉尘产生量。前两者直接影响到颗粒填充效果,进而关系到整个工程屏障的缓冲性能;而大量的粉尘将危害人体健康和机械设备的性能,降低工作面能见度,不利于准确操控机械设备,将间接地影响到颗粒填充效果。因此,提高填充干密度和均匀性、降低粉尘产生量是今后升级颗粒填充技术的关键突破点。

    堆积干密度和均匀性是影响颗粒混合物水力–力学特性的重要因素[10],是处置库概念设计的关键指标之一,通常通过颗粒填充试验来加以研究。按照填充区域的不同,颗粒填充试验可分为4类:①废物罐–围岩空隙填充试验;②矩形槽填充试验;③环形槽填充试验和④圆柱筒填充试验。

    废物罐–围岩空隙填充试验是将颗粒混合物填充到真实的或模拟的废物罐–围岩空隙区域内,然后取样测试堆积干密度和均匀性。在瑞士Mont Terri地下实验室开展的EB试验和HE-E试验中,都采用螺旋输送法将膨润土颗粒填充到废物罐与围岩之间的空隙区域,实测平均堆积干密度分别为1.36,1.457 g/cm3,且都观测到了不同程度的颗粒离析现象[23-24]。文献[25]分析认为粗细颗粒沿堆积坡面的滚落速率和距离不同、细颗粒填充粗颗粒间的大孔隙是导致颗粒离析的主要原因。

    矩形槽与环形槽填充试验是采用自由倾倒法或喷射法将颗粒混合物填充到长方体透明箱(图6(a))或圆环形透明槽缝(图6(b))中,通过测算颗粒总质量和填充总体积来计算堆积干密度,借助图像对比来分析堆积均匀性,重点考察了填充技术和级配的影响。研究表明[26-27]:①单一粒组、两粒组和满足Fuller级配的多粒组颗粒混合物的堆积干密度依次增大;②采用气动喷射、锤击和振动等手段能显著提高颗粒堆积干密度;③堆积均匀性由好到次依次为自由倾倒并振实、自由倾倒不振实和喷射法。

    图  6  矩形槽和环形槽填充试验[26-27]
    Figure  6.  Pellet filling tests in rectangle and annular gaps[26-27]

    圆柱筒填充试验是将颗粒混合物倒入圆柱形透明容器中,经反复多次拍击、敦实、摇晃、翻转和振动容器,分别测算松散和致密堆积干密度,并通过对比试验前后的图像或级配来分析堆积均匀性。研究表明[11,14-17]:①敦实、摇晃、翻转和振动等可显著提高颗粒堆积干密度;②随着细颗粒含量从0增大至100%,两粒组混合物的堆积干密度先增大后减小;③随着级配模量从0.1增大至0.9,多粒组颗粒混合物的堆积干密度先增大后减小;④峰值堆积干密度随最小粒径与最大粒径之比的增大而减小;⑤细颗粒的填充效应和疏松效应、粗颗粒的挤占效应和壁效应以及颗粒间的楔效应是影响颗粒堆积干密度的主要机制;⑥颗粒混合物的堆积干密度可由非线性堆积模型预测[18];⑦颗粒混合物的堆积均匀性随粒径特征指数(与粗细粒径比和质量比有关)的增大而减小[6,19]。此外,堆积方式也对堆积干密度和均匀性具有重要影响[12]

    综上所述,颗粒混合物的堆积干密度和均匀性与填充技术、级配和堆积方式等因素有关。大多数研究主要侧重于堆积干密度,而对堆积均匀性研究较少且缺乏定量化研究。实际上,由于不同粒径的颗粒具有不同的运动特性,堆积混合物中必然出现颗粒离析现象和局部颗粒间大孔隙,进而可能引发管涌、侵蚀和失稳等工程灾害。因此,研究颗粒离析的物理力学机制、定量描述堆积均匀性和提出改善堆积均匀性的措施,是今后需要进一步开展的工作[6,19]

    图1所示,处置库运营过程中核废物产生的核素衰变热将经废物罐和缓冲/回填材料向围岩传递。作为一种缓冲/回填材料,颗粒混合物是介于废物罐与围岩之间的热传递介质,其热传导特性决定了处置库近场和远场的温度分布。因此,热传导特性是颗粒混合物的一项重要研究内容。

    颗粒混合物的热传导系数可通过热探针法、热平板法、热流计法和热箱法4种途径进行测试(图7[8,28],其装置和原理分别如图7(a)~(d)所示。

    图  7  颗粒混合物的热传导系数测试技术[8,28]
    Figure  7.  Measuring techniques of thermal conductivity for pellet mixtures[8,28]

    热探针法和热平板法属于瞬态法,通过探针(或平板)以恒定功率释放热量并监测邻近颗粒温度随时间的变化来计算热传导系数;热流计法和热箱法属于稳态法,通过施加恒定温度梯度并监测流过单位面积试样的热流量来计算热传导系数。

    热探针法和热平板法仪器便携易用、测试效率高,但测试结果受探针(或平板)与颗粒的局部接触状态影响较大。热流计法和热箱法将颗粒混合物整体作为测试对象,克服了与颗粒的接触问题,测试结果较为可靠,但设备笨重不便携,测试效率低且成本较高(单次试验需要耗费大量的颗粒材料)。为了既高效又可靠地获得颗粒混合物在不同水力–力学状态下的热传导系数,应基于瞬态法和稳态法的基本原理,借鉴松散颗粒介质的热传导测试技术,同时考虑膨润土颗粒的特殊性(不均匀性、膨胀性和结构随吸力演化等),研发一套适用于膨润土颗粒混合物、颗粒与膨润土块体界面、颗粒与围岩界面的热传导测试技术和设备。

    利用不同的测试技术,不同学者对不同种类的颗粒混合物的热传导系数进行了试验研究。研究表明,颗粒混合物的热传导系数随干密度的增大而增大,随含水率的增大而增大,随温度的升高而增大(图8[11,25,28-30]

    图  8  热传导系数随干密度的变化[11,25]
    Figure  8.  Evolution of thermal conductivity with dry density[11,25]

    综上所述,目前对颗粒混合物热传导系数的研究还很不充分,现有文献报导中仅考虑了干密度、含水率和温度等影响因素。实际上,在处置库运营过程中,颗粒混合物将逐渐吸水饱和,同时伴随着强烈的结构特征演化过程。在这一过程中,颗粒混合物热传导系数的变化规律对处置库温度场预测和安全性能评价至关重要。因此,全饱和度(吸力、含水率)范围内颗粒混合物热传导系数的变化规律、基于结构特征演化的热传导系数变化机理以及改善颗粒混合物热传导性能的技术措施,是今后值得研究的几个课题。

    由于地下水的入渗过程极其缓慢,处置库中的颗粒混合物将长期处于非饱和状态。研究非饱和颗粒混合物的水力性质(持水特性、渗透特性、管涌与侵蚀等),是开展处置库安全评价的重要基础。

    研究表明,随着吸力的降低,单个颗粒的总孔隙比持续增大,但饱和度先增大后基本保持不变[31-33]。在恒体积条件下的高吸力范围内,不同干密度颗粒混合物的持水曲线基本重合(图9),表明水主要存在于颗粒内部孔隙中,持水性能主要取决于膨润土的矿物成分,而与干密度无关;在低吸力范围内,颗粒内部孔隙达到饱和状态,一部分水进入颗粒间大孔隙,此时持水曲线与干密度有关[10]

    图  9  不同干密度FEBEX膨润土颗粒混合物的持水曲线[10]
    Figure  9.  Curves of water retention of FEBEX bentonite pellets with different dry densities[10]

    实际处置库运营期间,膨润土颗粒可能会遭受温度循环和干湿循环作用[13,34],对持水曲线产生温度效应和滞回效应。温度对膨润土颗粒的持水特性及其滞回性的影响,仍是有待研究的课题。

    渗透特性包括饱和渗透与非饱和渗透两个方面。研究表明,膨润土颗粒混合物的饱和渗透系数随干密度的增大而增大[10,32,35]。在相同干密度条件下,颗粒混合物的饱和渗透系数与压实块体近似相等[6]。然而,与压实膨润土相比,水化过程中颗粒混合物的孔隙结构变化更加剧烈,导致渗透系数随水化时间持续发生显著变化。如图10所示,渗透试验开始时,颗粒之间存在大量相互连通的大孔隙,故渗透系数较大;随着试验的进行,颗粒吸水膨胀导致颗粒间大孔隙急剧减小甚至被堵塞,因此渗透系数迅速降低至稳定值。当饱和渗透系数达到稳定状态后,尽管试样外观看起来已经均一化,但干密度和含水率分布均表明试样仍未达到均一化[6]

    图  10  FEBEX膨润土颗粒混合物渗透系数随时间的演化[10]
    Figure  10.  Evolution of hydraulic conductivity with time for FEBEX bentonite pellets[10]

    颗粒混合物的非饱和渗透系数尚鲜见研究报导。实际上,文献[34,36]已经对MX80膨润土颗粒混合物开展了小型水–力耦合模型试验,分别监测了试样上下端同时通水和仅底端通水过程中不同截面的相对湿度变化,但未给出相应的非饱和渗透系数函数。文献[6,37]采用瞬时截面法测试或预测了不同条件(温度、干密度等)下颗粒混合物的非饱和渗透系数,结果表明,颗粒混合物的非饱和渗透系数随吸力的降低而先减小后增大,随温度的升高而增大,随干密度的增大而减小,随颗粒粒径的增大而增大。然而,地下水化学对非饱和渗透特性的影响尚待研究。

    除了渗流与结构演化,颗粒混合物水化过程中的另一个重要性质是管涌与侵蚀。许多学者通过圆柱管注水试验[7]、矩形槽注水试验[7]、环形槽注水试验[27]和模拟巷道接缝注水试验[38]研究发现,开始注水时水流优先沿着相互连通的颗粒间大孔隙流动,并侵蚀沿程细颗粒和膨润土胶体。随着颗粒吸水膨胀和胶体阻塞,颗粒间大孔隙逐渐闭合,管涌和侵蚀现象逐渐减弱直至消失。

    圆柱管注水试验[7]结果表明,管涌持续时间与侵蚀速率主要同注水速率有关:注水速率越大,则管涌持续时间越长,侵蚀速率也越大;在双对数坐标下,累计侵蚀量与累计水流量呈线性关系。此外,膨润土类型、颗粒形状与尺寸、细颗粒的分布情况、溶液盐度也对管涌和侵蚀具有显著影响,但由于缺乏系统的研究,这些因素的影响规律仍不明确。矩形槽[7]、环形槽[27]和模拟巷道接缝[38]注水试验从更大尺度上进一步研究了注水速率、溶液浓度、颗粒形状和大小及膨润土类型等因素对管涌行为的影响。但由于试验成果较少且不够系统,尚不足以形成完整的结论。

    不同于非膨胀性颗粒土,膨润土颗粒混合物的管涌和侵蚀行为是一个更加复杂的过程。管涌和侵蚀过程中,颗粒混合物中的细小颗粒或胶体被优势渗流侵蚀流失,导致颗粒混合物的整体干密度降低,进而引起膨胀力降低而饱和渗透系数增大,使其缓冲性能大幅衰减。在已有认知的基础上,开展管涌和侵蚀的影响因素及其机理研究并提出抑制管涌和侵蚀的技术措施,是今后针对颗粒混合物的重要研究方向之一。

    颗粒及其混合物的热–水–力特性与其结构特征密切相关。目前,关于颗粒及其混合物的结构演化规律的研究成果较少,仅部分学者采用MIP和μ-CT两种手段分别观测和分析了单个颗粒和颗粒混合物在水化过程中的结构演化规律。

    采用μ-CT扫描和MIP试验观测表明[31],吸力113 MPa的由压实法制备的单个颗粒内部尤其是靠近边缘部分存在明显的裂隙,其孔径分布呈双峰形态(图11(a)和(e))。随着吸力的降低,由于颗粒内部不同性质的物质吸水发生不均匀膨胀,导致颗粒内部的裂隙数量增多,裂隙连通性增强,颗粒宏观体积增大(图11(b)~(d));集合体内孔隙数量减少,集合体间孔隙数量增多且平均孔径增大(图11(f)~(h))。

    图  11  不同吸力下单个颗粒的结构特征[31]
    Figure  11.  Structural characteristics of single pellet at different suctions[31]

    在通水水化过程中,颗粒混合物的结构特征将发生显著改变。文献[39]采用μ-CT技术观测了FoCa膨润土颗粒/粉末(50/50)混合物在仅底端通水情况下的结构特征演化过程。结果显示(图12(a)),随着水流自下而上进入试样,膨润土颗粒由于吸水膨胀而密度逐渐降低,粉末由于受颗粒膨胀的压缩作用而密度逐渐增大,颗粒轮廓逐渐模糊,最终经历165 d后试样各处的干密度近似相等,基本达到了均一化状态。

    图  12  水化过程中颗粒混合物的结构演化[39-40]
    Figure  12.  Structural evolution of pellet mixtures during hydration[39-40]

    文献[31,40]采用μ-CT技术观测了MX80膨润土颗粒/粉末(80/20)在上下两端同时通水水化过程中的结构演化过程。结果表明,试样上、下两端的颗粒迅速水化膨胀,而试样中部的颗粒水化膨胀过程极其缓慢,及至水化56 d时部分颗粒和颗粒间孔隙仍然清晰可辨(图12(b))。最终经历100 d后,试样各处的干密度近似相等,基本达到了均一化状态。

    综上所述,目前主要采用μ-CT和MIP两种手段分别观测了单个颗粒在不同吸力状态下和颗粒混合物在通水水化过程中的结构特征演化规律,而对于颗粒混合物在不同吸力状态下的结构特征演化规律尚缺乏研究。此外,温度、孔隙水化学、外荷载等因素对颗粒混合物结构特征的影响也有待研究。

    颗粒应具有足够的强度和耐磨性,以减小在滚动、撞击和压实过程中的破损程度。对于单个颗粒而言,强度和耐磨性是评价其工程性能的重要力学指标。颗粒的强度常采用恒应变速率液压试验机(图13(a))和手动压力仪(图13(b))来测定,颗粒的耐磨性常采用耐磨试验机(图13(c))来测试[8]。研究表明[8,11],在一定含水率范围内,颗粒强度随含水率的增大而增大,且同一含水率下粒径大的颗粒强度更高;颗粒强度越高,耐磨性越好。

    图  13  颗粒强度和耐磨性试验设备[8]
    Figure  13.  Devices used in pellet strength and abrasion tests[8]

    目前,主要针对室温条件下、含水率低和形状规则的颗粒开展了简单的强度和耐磨性试验。今后应针对不同温度、不同含水率和不同形状的颗粒开展无侧限压缩试验、抗剪强度试验和多种形式的耐磨性试验(不同倾角溜滑、不同高度倾倒、不同转速翻滚、不同频率振动等),综合评价其强度和耐磨性。

    膨润土颗粒混合物的主要矿物成分为蒙脱石,具有显著的吸水膨胀特性。在恒体积条件下,蒙脱石吸水膨胀在宏观上表现为膨胀力,在非恒体积条件下则表现为膨胀变形。研究表明,颗粒混合物的最终膨胀力随试样干密度的增大而增大,随掺砂率的增大而减小,但膨胀力发展过程受试样尺寸的影响较大[9-10,35]。与压实膨润土块体类似[41],不同干密度和不同尺寸条件下,FoCa颗粒/粉末混合物(50/50)试样的膨胀力时程曲线均呈双峰形态(图14):膨胀力随水化时间先迅速增大至一个峰值,随后回落至一个谷值,然后以较小的速率再次增大至稳定值[9]。对于满足Fuller级配的FEBEX膨润土颗粒混合物(干密度1.36 g/cm3),其膨胀力时程曲线在达到峰值后基本不变,然后再逐渐增长至稳定值[42]图14)。由此可见,颗粒混合物的膨胀力发展过程与粒径级配有关。此外,随着吸力从300 MPa逐级降低至3 MPa,颗粒混合物的膨胀力先增大后轻微减小,膨胀力减小与应力路径贴近LC屈服线有关[10,43]

    图  14  不同条件下颗粒混合物的膨胀力时程曲线[9,42]
    Figure  14.  Evolution of swelling pressure with time for pellet mixtures under different conditions[9,42]

    一些学者研究了温度对颗粒混合物膨胀力的影响,但研究结果不尽相同。对于MX80颗粒混合物,文献[44]发现膨胀力随温度的升高而降低;文献[34]的研究则表明膨胀力随温度的升高而增大,而文献[35]发现膨胀力受温度的影响不明显。

    在恒定荷载下通水水化,颗粒混合物试样膨胀稳定时的变形与干密度和竖向荷载有关[10]。在相同干密度条件下,颗粒混合物的膨胀变形与竖向荷载的对数呈线性关系;在同一竖向荷载条件下,干密度越大,膨胀变形越大(图15(a))。在恒定荷载下逐级降低吸力(气相法),低干密度颗粒混合物的变形先增大后轻微减小,而高干密度时则持续增大[10]图15(b))。

    图  15  膨胀变形与竖向荷载及吸力的关系[10]
    Figure  15.  Evolution of swelling strain with vertical stress and suction[10]

    此外,通水水化时试样先小幅坍塌然后膨胀,膨胀变形达到峰值后又再次坍塌;气相法逐渐减低吸力时,不存在初始坍塌现象,试样持续膨胀到峰值后坍塌至稳定值[45]。造成这种区别的原因在于两种水化过程中的颗粒–水相互作用机理不同:通水水化过程中,颗粒间孔隙迅速被水饱和,导致颗粒间相互作用力减弱,在竖向荷载作用下颗粒间错动而引发坍塌;气相法降低吸力过程中,颗粒间孔隙一直处于非饱和状态,颗粒内、外的孔隙水势差异很小,毛细吸力引起的颗粒错动较小,而颗粒吸水膨胀相对显著,故不存在初始坍塌。

    尽管对颗粒混合物的膨胀特性已经开展了不少研究工作,但温度、级配及孔隙水化学等因素对膨胀力和膨胀变形的影响及膨胀力和膨胀变形的形成机理与发展过程仍有待进一步研究。

    文献[46]对7种不同膨润土颗粒混合物分别开展了3种不同级配下的击实试验,结果表明:压实性质主要受膨润土类型影响,而受级配影响较小;最优含水率随塑限的增大而线性增大,最大干密度随塑限的增大而减小。文献[20,29]分别对不同级配的巩义膨润土颗粒混合物和高庙子膨润土颗粒混合物开展了静力压实试验,结果表明级配对最优含水率的影响较小,但最大干密度随细颗粒含量的增多而增大。

    压缩试验结果表明[10,29],颗粒混合物的压缩系数随竖向应力的增大而先增大后减小,随含水率的增大而增大;弹性压缩系数κ和塑性压缩系数λ均随吸力的降低而增大(图16);级配对压缩指数与回弹系数等压缩性能具有一定影响,但主要表现在含水率较低的情况下。

    图  16  弹、塑性压缩系数随吸力的变化[10]
    Figure  16.  Evolution of elastoplastic compressibility coefficients (κ and λ) with suction[10]

    基于双孔结构理论,文献[45,47]采用宏观结构吸力和微观结构吸力分别描述颗粒间孔隙水势和颗粒内孔隙水势,采用BBM模型和BEXM模型分别描述微、宏观结构变形,分别建立了适用于描述颗粒混合物水力–力学行为的本构模型。所不同的是,前者将颗粒混合物简化为一维二元介质模型,并假设微、宏观相对渗透系数与相应的饱和度呈幂函数关系,后者假设颗粒间与颗粒内的孔隙水交换呈线性关系且仅由吸力驱动。

    目前,关于颗粒混合物的力学性质研究还处于宏观描述阶段,而缺乏对其微观力学行为(如颗粒运动、破碎、膨胀与胶结融合等)的研究。在今后的研究中,应综合借鉴非膨胀颗粒土和压实膨胀性土的研究经验,借助μ-CT、NMRI等成像技术,揭示颗粒混合物在不同应力状态下的微观力学行为机理,进而建立颗粒混合物的宏–细–微三重尺度耦合理论及本构模型。

    膨润土颗粒混合物是一种重要的缓冲/回填材料。20世纪80年代末以来,国内外学者针对膨润土颗粒的制备方法、填充技术与堆积性质、热传导特性、水力特性、结构演化规律及力学特性等方面开展了大量研究,主要成果如下:

    (1)膨润土颗粒可由挤压法、辊压法、压实法、压实–破碎法和湿–干–破碎法制备,可通过带式输送法、螺旋输送法、气动喷射法和人工填充法进行填充。

    (2)颗粒混合物的堆积干密度和均匀性与填充技术、级配、堆积方式等因素有关,其热传导系数与干密度、含水率和温度等因素有关,其水力–力学特性与级配、干密度及温度等因素有关。

    (3)水化过程中颗粒混合物由初始松散结构逐渐转变为胶结融合结构,及至水化饱和后基本达到宏观上的均一化结构,但微观层次的均一化过程仍将持续漫长的时间。

    纵观国内外文献报导,目前对膨润土颗粒混合物的研究范围较为广泛,但研究深度仍相对较浅。考虑到处置库实际运营工况的复杂性,科学高效的颗粒混合物填充技术、多场(热–水–化–力)耦合条件下的颗粒混合物水力–力学特性及结构演化规律是今后值得深入探索的研究方向。

  • 图  1   压实膨润土双孔隙结构示意图

    Figure  1.   Schematic view of compacted dual-porosity bentonite

    图  2   单元八方向生长示意图

    Figure  2.   Schematic view of eight growing directions for each cell

    图  3   不同孔隙比下压实膨润土双孔隙建模(黑色:固相颗粒;白色:孔隙气)

    Figure  3.   Modeling of dual-porosity for compacted bentonite with different void ratios (black: solid particle; white: pore air)

    图  4   双孔隙中有偏水分布建模(黑色:固相颗粒;灰色:孔隙水;白色:孔隙气)

    Figure  4.   Modeling of preferential water distribution within dual-porosity (black: solid particle; grey: pore water; white: pore air)

    图  5   稳态导热模拟的边界条件

    Figure  5.   Boundary conditions of steady heat conduction

    图  6   单向网格数优化分析

    Figure  6.   Optimization analysis of grid number

    图  7   基于Mx-80膨润土的其他三类膨润土em估计

    Figure  7.   Estimation of em for other three compacted bentonites based on Mx-80 bentonite

    图  8   5种膨润土导热系数模拟值与实测值对比图

    Figure  8.   Comparison between experimental and modeled thermal conductivities of five compacted bentonites

    图  9   Chen等模型与本文模型水分布比较

    Figure  9.   Comparison of water distribution between Chen et al. model and proposed model

    表  1   不同相对误差下压实膨润土导热系数预测的置信度

    Table  1   Confidence degrees of predictive thermal conductivity of compacted bentonite under different relative errors

    膨润土Mx-80FebexKyeongjuKunigel-V1GMZ-07
    样本个数2237711416
    D≤5%41%24%21%29%25%
    D≤15%91%81%53%71%81%
    D≤25%95%97%86%79%88%
    下载: 导出CSV
  • [1] 王驹. 中国高放废物地质处置21世纪进展[J]. 原子能科学技术, 2019, 53(10): 2072-2082. https://www.cnki.com.cn/Article/CJFDTOTAL-YZJS201910036.htm

    WANG Ju. Progress of geological disposal of high-level radioactive waste in China in the 21st century[J]. Atomic Energy Science and Technology, 2019, 53(10): 2072-2082. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YZJS201910036.htm

    [2] 陈宝, 张会新, 陈萍. 高碱溶液对高庙子膨润土侵蚀作用的研究[J]. 岩土工程学报, 2013, 35(1): 181-196. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201301022.htm

    CHEN Bao, ZHANG Hui-xin, CHEN Ping. Erosion effect of hyper-alkaline solution on gaomiaozi bentonite[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(1): 181-186. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201301022.htm

    [3] 刘月妙, 蔡美峰, 王驹. 内蒙古高庙子钠基膨润土导热性能研究. 全国岩土与工程学术大会[C]//全国岩土与工程学术大会. 2006, 武汉.

    LIU Yue-miao, CAI Mei-feng, WANG Ju. On the thermal conductivity of GMZ bentonite[C]//Proceedings of Rock and Engineering Conference. 2006, Wuhan. (in Chinese)

    [4] 叶为民, 王琼, 潘虹, 等. 高压实高庙子膨润土的热传导性能[J]. 岩土工程学报, 2010, 32(6): 821-826. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201006003.htm

    YE Wei-min, WANG Qiong, PAN Hong, et al. Thermal conductivity of compacted GMZ01 bentonite[J]. Chinese Journal of Geotechnical Engineering, 2010, 32(6): 821-826. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201006003.htm

    [5]

    KAHR G, MÜLLER-VONMOOS M. Wärmeleitfähigkeit von Bentonite MX80 und von Montigel nach der Heizdrahtmethode[R]. Zurich: Swiss National Cooperative for the Disposal of Radioactive Waste, NAGRA Technisher Bericht, 1982.

    [6]

    KNUTSSON S. On the Thermal Conductivity and Thermal Diffusivity of Highly Compacted Bentonite[R]. Stockholm: Swedish Nuclear Fuel and Waste Management Co, 1983.

    [7]

    SAKASHITA H, KUMADA T. Heat transfer model for predicting thermal conductivity of highly compacted bentonite[J]. Journal of the Atomic Energy Society of Japan, 1998, 40(3): 235-240.

    [8]

    ROMERO E. A microstructural insight into compacted clayey soils and their hydraulic properties[J]. Engineering Geology, 2013, 165(20): 3-19.

    [9] 叶为民, 赖小玲, 刘毅, 等. 高庙子膨润土微观结构时效性试验研究[J]. 岩土工程学报, 2013, 35(12): 2255-2261. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201312020.htm

    YE Wei-min, LAI Xiao-ling, LIU Yi, et al. Experimental study on ageing effects on microstructure of unsaturated GMZ01 bentonite[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(12): 2255-2261. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201312020.htm

    [10]

    LLORET A, VILLAR M V. Advances on the knowledge of the thermo-hydro-mechanical behaviour of heavily compacted “FEBEX” bentonite[J]. Physics and Chemistry of the Earth, 2007, 32(8): 701-715.

    [11]

    TANG A M, CUI Y J, LE T T. A study on the thermal conductivity of compacted bentonites[J]. Applied Clay Science, 2008, 41(3/4): 181-189.

    [12]

    CHEN Y F, WANG M, ZHOU S, et al. An effective thermal conductivity model for unsaturated compacted bentonites with consideration of bimodal shape of pore size distribution[J]. Science China Technological Sciences, 2015, 58(2): 369-380. doi: 10.1007/s11431-014-5738-3

    [13]

    CAI S S, ZHANG B X, CUI T F, et al. Mesoscopic study of the effective thermal conductivity of dry and moist soil[J]. International Journal of Refrigeration, 2019, 98(6): 171-181.

    [14] 周殷康, 阎长虹, 谢胜华, 等. 基于细观模拟的软土导热系数数值预测模型[J]. 工程地质学报, 2019, 27(5): 1070-1078. doi: 10.13544/j.cnki.jeg.2019078

    ZHOU Yin-kang, YAN Chang-hong, XIE Sheng-hua, et al. A numerical model for thermal conductivity of soft soils based on mesoscopic simulation[J]. Journal of Engineering Geology, 2019, 27(5): 1070-1078. (in Chinese) doi: 10.13544/j.cnki.jeg.2019078

    [15]

    JIANG Z, DIJKE M I J, WU K, et al. Stochastic pore network generation from 3D rock images[J]. Transport in Porous Media, 2012, 94(2): 571-593. doi: 10.1007/s11242-011-9792-z

    [16]

    JIANG Z, VAN DIJKE M I J, SORBIE K S, et al. Representation of multiscale heterogeneity via multiscale pore networks[J]. Water Resources Research, 2013, 49(9): 5437-5449. doi: 10.1002/wrcr.20304

    [17] 徐云山, 孙德安, 曾召田, 等. 膨润土热传导性能时效性试验研究[J]. 岩土力学, 2019, 40(11): 4324-4330. doi: 10.16285/j.rsm.2018.1707

    XU Yun-shan, SUN De-an, ZENG Zhao-tian, et al. Experimental study on ageing effect on bentonite thermal conductivity[J]. Rock and Soil Mechanics, 2019, 40(11): 4324-4330. (in Chinese) doi: 10.16285/j.rsm.2018.1707

    [18]

    DELAGE P, MARCIAL D, CUI Y J, et al. Ageing effects in a compacted bentonite: a microstructure approach[J]. Géotechnique, 2006, 56(5): 291-304. doi: 10.1680/geot.2006.56.5.291

    [19]

    CHEN Y F, ZHOU S, HU R, et al. A homogenization-based model for estimating effective thermal conductivity of unsaturated compacted bentonites[J]. International Journal of Heat and Mass Transfer, 2015, 83: 731-740. doi: 10.1016/j.ijheatmasstransfer.2014.12.053

    [20]

    WANG M R, WANG J K, PAN N, et al. Mesoscopic predictions of the effective thermal conductivity for microscale random porous media[J]. Physical Review E, 2007, 75(3): 036702. doi: 10.1103/PhysRevE.75.036702

    [21]

    WANG M R, PAN N. Predictions of effective physical properties of complex multiphase materials[J]. Materials Science and Engineering: R: Reports, 2008, 63(1): 1-30. doi: 10.1016/j.mser.2008.07.001

    [22]

    VILLAR M V. Caracterización termo-hidro-mecánica de una bentonita de Cabo de Gata[D]. Madrid: Universidad Complutense de Madrid, 2000.

    [23]

    LEE J O, CHOI H, LEE J Y. Thermal conductivity of compacted bentonite as a buffer material for a high-level radioactive waste repository[J]. Annals of Nuclear Energy, 2016, 94: 848-855. doi: 10.1016/j.anucene.2016.04.053

    [24]

    OULD-LAHOUCINE C, SAKASHITA H, KUMADA T. Measurement of thermal conductivity of buffer materials and evaluation of existing correlations predicting it[J]. Nuclear Engineering and Design, 2002, 216(1/2/3): 1-11.

    [25]

    NAKASHIMA Y. Nuclear magnetic resonance properties of water-rich gels of kunigel-V1 bentonite[J]. Journal of Nuclear Science and Technology, 2004, 41(10): 981-992. doi: 10.1080/18811248.2004.9726321

    [26]

    TANG A M, CUI Y J. Modelling the thermo-mechanical volume change behaviour of compacted expansive clays[J]. Géotechnique, 2009, 59(3): 185-195. doi: 10.1680/geot.2009.59.3.185

    [27] 徐云山, 曾召田, 吕海波, 等. 高温下红黏土热导率的变化规律试验研究[J]. 工程地质学报, 2017, 25(6): 1465-1473. doi: 10.13544/j.cnki.jeg.2017.06.009

    XU Yun-shan, ZENG Zhao-tian, LÜ Hai-bo, et al. Experimental study on variation of thermal conductivity of red clay at high temperature[J]. Journal of Engineering Geology, 2017, 25(6): 1465-1473. (in Chinese) doi: 10.13544/j.cnki.jeg.2017.06.009

  • 期刊类型引用(7)

    1. 刘桃根,李玲,王伟,谢敬礼,刘造保. 黏土-砂混合物三轴压缩超声波特性研究. 地下空间与工程学报. 2024(S1): 119-127 . 百度学术
    2. 陈贺新,胡立强,项国圣. 盐溶液作用下压实膨润土体变性能研究. 工程勘察. 2023(04): 1-6+35 . 百度学术
    3. 张小小,陈永贵,李昆鹏,叶为民,王琼. 前期化学循环作用下石墨烯改性GMZ膨润土持水特性. 岩土工程学报. 2023(07): 1490-1497 . 本站查看
    4. 胡志杰,项国圣,付文青,王浩,李华健. 荷载-溶液作用下膨润土压缩渗透性能研究. 地下空间与工程学报. 2023(06): 1851-1858 . 百度学术
    5. 项国圣,吕立勇,葛磊,周殷康,谢胜华. 温度对GMZ膨润土的膨胀性能影响研究. 岩土工程学报. 2021(01): 77-84 . 本站查看
    6. 胡云世,徐云山,孙德安,陈波,曾召田. 温度对颗粒膨润土热传导特性的影响. 岩土力学. 2021(07): 1774-1782 . 百度学术
    7. 刘俊新,唐伟,李军润,张建新,郭招群,陈龙,刘育田. 高温及碱性条件对高庙子钠基膨润土膨胀力的影响. 岩土力学. 2021(08): 2160-2172+2184 . 百度学术

    其他类型引用(7)

图(9)  /  表(1)
计量
  • 文章访问数:  260
  • HTML全文浏览量:  28
  • PDF下载量:  116
  • 被引次数: 14
出版历程
  • 收稿日期:  2020-07-16
  • 网络出版日期:  2022-12-02
  • 刊出日期:  2021-06-30

目录

/

返回文章
返回