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

脱水泥饼堆填体的压缩及渗透特性试验研究

王子昂, 詹良通, 郭晓刚, 刘超洋, 尧俊

王子昂, 詹良通, 郭晓刚, 刘超洋, 尧俊. 脱水泥饼堆填体的压缩及渗透特性试验研究[J]. 岩土工程学报, 2021, 43(10): 1915-1923. DOI: 10.11779/CJGE202110018
引用本文: 王子昂, 詹良通, 郭晓刚, 刘超洋, 尧俊. 脱水泥饼堆填体的压缩及渗透特性试验研究[J]. 岩土工程学报, 2021, 43(10): 1915-1923. DOI: 10.11779/CJGE202110018
WANG Zi-ang, ZHAN Liang-tong, GUO Xiao-gang, LIU Chao-yang, YAO Jun. Experimental study on compressibility and permeability characteristics of dewatered slurry cake backfill[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(10): 1915-1923. DOI: 10.11779/CJGE202110018
Citation: WANG Zi-ang, ZHAN Liang-tong, GUO Xiao-gang, LIU Chao-yang, YAO Jun. Experimental study on compressibility and permeability characteristics of dewatered slurry cake backfill[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(10): 1915-1923. DOI: 10.11779/CJGE202110018

脱水泥饼堆填体的压缩及渗透特性试验研究  English Version

基金项目: 

浙江省重点研发计划项目 2019C03107

苏州市民生科技项目—科技示范工程项目 SS201804-01

详细信息
    作者简介:

    王子昂(1997— ),男,硕士研究生,主要从事环境岩土工程方面研究工作。E-mail:ziangwang1226@zju.edu.cn

    通讯作者:

    詹良通, E-mail:zhanlt@zju.edu.cn

  • 中图分类号: TU441

Experimental study on compressibility and permeability characteristics of dewatered slurry cake backfill

  • 摘要: 脱水泥饼是工程泥浆经压滤脱水后形成的块状产物,经装卸、转运等处理环节破碎成尺寸不一的碎散体,其堆填体具有明显的双孔隙特征,由块体间的大孔隙和块体内的小孔隙组成。为了指导泥饼堆填作业,对泥饼堆填体的压缩性及渗透性开展研究。利用大尺寸固体废弃物压缩-渗透装置模拟了泥饼堆填过程中上覆荷载增加,泥饼堆积单元体逐渐被压密的过程,并利用CT扫描对压缩后泥饼块体间大孔隙的变化情况进行了分析;同时测试了泥饼堆积单元体在压缩过程中的渗透性变化规律。试验结果表明:泥饼堆积单元体试样随着上覆荷载增加逐渐被压密,泥饼间孔隙与泥饼内孔隙的尺度差异逐渐减小,双孔隙特征减弱。在上覆荷载为295 kPa时,试样空气孔隙率为5%,此时试样已成为均质块状泥饼,对应压实度为0.87,试样透水性已由初始阶段(0 kPa)的中等透水等级转变为弱透水等级。考虑可操作性和经济性,建议将0.85作为泥饼堆填作业过程中的压实度控制指标。
    Abstract: After the slurry is dewatered by pressure filtration, the dewatered slurry cake is formed. After loading and unloading, transportation and other processing links, the dewatered slurry cake is broken into pieces of different sizes. Dewatered slurry cake backfill has obvious dual-porosity characteristics, including macro-pores between cakes and micro-pores inside cakes. In order to guide the filling operation of the dewatered slurry cake, the compressibility and permeability of the dewatered slurry cake backfill are studied in this research. The compaction of the dewatered slurry cake as the load increases during the backfilling is simulated by using the large-scale solid waste compression-permeation experimental system. The change of pore structure in the sample during compression is analyzed based on the CT scan. And the change of permeability of dewatered slurry cake backfill in the process of compression is studied at the same time. The results show that the sample of the dewatered slurry cake is compacted gradually with the increase of external load. The scale difference between the macro-pores and the micro-pores decreases gradually, and the dual-porosity characteristics weakens with the process of compression. When the overlying load is 295 kPa, the volume ratio of air of the sample is 5%. At this time, the sample becomes a homogeneous lump slurry cake, and the corresponding compaction degree is 0.87. In addition, the permeability of the sample also changes from a medium permeability grade at the initial stage (0 kPa) to a weak permeability grade. Considering the operability and economy, it is suggested that the compaction degree should be set as 0.85 during the backfilling of the dewatered slurry cake.
  • 在高放射性废物深地质处置中,压实膨润土因具良好的吸附性、膨胀性以及低渗透性被确定为首选缓冲/回填材料[1-3]。花岗岩等结晶岩因力学强度高,完整性较好而作为处置库围岩备选之一,如中国处置库预选区甘肃北山新场围岩为花岗岩[4-5]。然而,处置库建设过程中开挖扰动产生的围岩裂隙[6],以及岩体内部存在大量剪切形成的连通微裂隙,成为地下水的优势渗流通道[7]。地下水侵蚀作用下,处于裂隙围岩交界处的压实膨润土会释放1~1000 nm粒径的膨润土胶体(BC)到地下水中,这些胶体颗粒的主要矿物成分为蒙脱石,比表面积较大,对放射性核素具有较强的吸附和交换能力[8]。在有利地下水环境中,稳定、移动的膨润土胶体可通过吸附作用与放射性核素结合,沿导流裂隙快速迁移,从而威胁处置安全[9]图1)。

    图  1  膨润土胶体形成及核素迁移[9]
    Figure  1.  Colloid formation and radionuclide migration [9]

    通常,足够浓度的稳定胶体能吸附放射性核素,而胶体移动且吸附不可逆时,则会对核素迁移产生重要影响[10]。为此,国外学者开展了大量与膨润土胶体相关的试验和模拟研究。例如,欧盟关于工程屏障长期性能和核素迁移影响项目(BELBaR),重点关注了膨润土侵蚀和胶体的产生机理及其稳定性[11];瑞士GTS(Grimsel Test Site)试验场进行的原位胶体迁移和阻滞项目(CRR)[10]、胶体形成和迁移项目(CFM)[12],重点考察了膨润土胶体相关的放射性核素迁移行为。目前,中国学者主要针对高庙子(GMZ)膨润土胶体的稳定性[13-15]、吸附性和迁移特性以及胶体和单一核素(U(Ⅵ)、Eu(Ⅲ))的共迁移特性开展了相关试验研究[16-17],揭示各因素影响下的作用规律和机理,但缺乏原位试验和模拟研究。

    本文全面归纳和总结处置库膨润土胶体的吸附迁移特性,阐述了膨润土胶体-核素共迁移的试验和模拟研究成果,并进一步揭示了胶体和核素的共迁移作用机理,以期为中国高放废物处置库工程屏障系统的设计和长期性能评估提供依据。

    膨润土胶体吸附、迁移特性影响胶体与核素共迁移作用,已有研究揭示了胶体浓度、地下水离子强度和pH等对胶体吸附迁移特性的影响规律和作用机理。

    膨润土胶体主要矿物成分为蒙脱石,粒径小,比表面积大,表现出对核素离子较强的吸附性。由批式吸附试验得到的吸附率q和吸附分配系数Kdsor分析发现,膨润土胶体对核素吸附性能受初始胶体浓度、核素浓度、地下水化学条件和温度影响明显,受胶体粒径影响不大[17-19]

    研究表明,胶体对核素吸附率随胶体浓度增加而增加。图2(a)中,GMZ膨润土胶体在1 mM NaCl电解质中对Eu(Ⅲ)的吸附率、吸附分配系数均随胶体浓度增加先呈指数增长,后保持不变;B3378膨润土胶体在离子强度为2 mM的模拟雨水中对Cs(Ⅰ)和U(Ⅵ)的吸附率随胶体浓度增加呈线性增长,而吸附分配系数随胶体浓度增加呈指数衰减。胶体浓度决定胶体表面的吸附位点数量,随浓度增加,有效吸附位点增加直至饱和[17-18]。核素初始浓度对吸附性影响也与胶体浓度相关,图2(b)中,当胶体浓度较低时,Cs(Ⅰ)和U(Ⅵ)的吸附率、吸附分配系数随核素初始浓度增加而下降;而胶体浓度较高时,两者随核素初始浓度增加先增加后减小。此外,Cs(Ⅰ),U(Ⅵ)共存时,竞争吸附会降低核素的吸附率[18]

    图  2  不同因素影响下的胶体吸附性
    Figure  2.  Adsorption properties of bentonite colloids under different influences

    地下水的离子强度、阳离子价态、pH等水化学条件对膨润土胶体吸附性影响显著[17-18]。当地下水离子强度较高或存在二价阳离子时,因离子竞争吸附或络合反应,吸附率明显下降。在高离子强度(169 mM)、较高浓度Na+,Ca2+的地下水中,B3378膨润土胶体对Cs(Ⅰ)和U(Ⅵ)的吸附率都低于60%[18]图2(c)中,当pH值小于7时,GMZ膨润土胶体对Eu(Ⅲ)的吸附率随离子强度和Eu(Ⅲ)浓度的增加呈明显下降趋势。膨润土胶体对Eu(Ⅲ)的吸附率在pH值小于7时,随pH值增加逐渐提高至90%,而后随pH值增加吸附率基本不变。这是因为随pH值增大,膨润土胶体对Eu(Ⅲ)的吸附由外表面络合转变为内表面络合,而内表面络合基本不受离子强度和pH影响[17-18]

    此外,膨润土胶体对Eu(Ⅲ)的吸附为自发吸热反应[17],升温可提高金属离子在吸附剂表面的亲和力或增加吸附剂表面的电荷和电势,从而提高吸附性[20]

    膨润土胶体在低离子强度、弱碱性地下水中呈现快速迁移的特性,表现为胶体先于示踪剂HTO达到穿透曲线的峰值,如图3所示[21]

    图  3  膨润土胶体的快速迁移特征[21]
    Figure  3.  Rapid migration characteristics of bentonite colloids [21]

    胶体在尺寸排阻和电荷排阻效应的耦合作用下,迁移速率可达平均水流速率的1.4倍。膨润土胶体相对于溶质离子粒径更大,在岩石基质中的扩散作用可忽略不计;在尺寸排阻效应作用下,胶体在裂隙中迁移时更容易进入水力传导性更强、孔隙更大的区域。此外,膨润土胶体和裂隙岩石表面在中性、碱性条件下均带负电荷,在电荷排阻效应作用下胶体被驱动着远离介质表面向流速更快的区域运动,进一步加快胶体迁移[22]。胶体的迁移速率可采用阻滞因子Rf定性评价[21]

    Rf=VtracVcol=1+ρKdθ, (1)

    式中,Vtrac为示踪剂流速,Vcol为胶体流速,Kd为吸附分配系数(cm),ρ为柱填充材料体积密度(g/mL),θ为孔隙度。

    胶体的迁移穿透能力采用最大穿透率(Max C/ C0)和回收率(质量回收率)Re进行评价[23-24]。穿透试验表明,同流速下测得的膨润土胶体回收率要明显低于示踪剂,是胶体过滤所致[21, 25]

    膨润土胶体的迁移穿透能力受胶体浓度以及地下水流速、离子强度和pH影响显著,如图4。胶体浓度较高时,黏度增加导致颗粒间作用增强,产生熟化过滤效应,使胶体回收率下降[26]。地下水流速越高,胶体在介质中迁移时间越短,水动力剪切作用越大,胶体和介质形成的联结容易被剪断,胶体不容易被介质过滤而滞留在裂隙表面,随流速增大回收率提高[21, 26]。地下水离子强度和pH通过改变胶体的粒径和电荷特性,使稳定性变化,随离子强度升高或酸性增强,分散、稳定的胶体形成大的絮凝体沉积下来,使胶体最大穿透率和回收率下降[15, 17, 25]。此外,当地下水由高离子强度转变为低离子强度时,过滤的胶体难以恢复迁移能力,而二价阳离子会大大降低胶体迁移性[17, 25]

    图  4  不同因素影响下的胶体迁移性[17, 25]
    Figure  4.  Migration properties of bentonite colloids under different influences [17, 25]

    膨润土胶体和核素共迁移试验主要采用动态柱和偶极子流场试验,关注移动和过滤胶体对核素迁移的作用。

    动态柱试验能够有效控制试验变量,常用于实验室研究迁移规律。装置见图5,由进液瓶、蠕动泵、层析柱和自动部分收集器组成,通过在岩块间、岩块和套筒之间预留接缝构造单一、平行裂隙[27-28],或在层析柱中填充花岗岩碎屑、石英砂等构造多孔介质来模拟岩石裂隙系统[16-17, 29],但难以反映真实处置库中复杂的岩石裂隙系统。试验过程包括预平衡、穿透和淋洗3个阶段,当背景电解质平衡好层析柱后开始穿透实验,注入一定孔隙体积的样品(浓度为C0)并间隔收集流出液(浓度为C)至浓度稳定后,改用背景电解质淋洗层析柱,直至流出液检测不到样品[16-17]

    图  5  柱试验装置示意图[17, 27-28]
    Figure  5.  Schematic diagram of column test devices [17, 27-28]

    动态柱试验结果(表1)表明,膨润土胶体和核素的共迁移作用依赖于地下水离子强度和pH。低离子强度、碱性地下水中,胶体稳定、可移动,可促进核素迁移[16, 21, 29];反之,高离子强度、存在二价阳离子、酸性地下水中,胶体不稳定,反而阻滞核素迁移[17, 30]

    表  1  离子强度和pH对共迁移影响
    Table  1.  Effect of ionic strength and pH on co-migration
    核素-胶体电解质(Ⅰ; pH)胶体状态穿/%效果文献
    Cs(Ⅰ)DW稳定、可移动0.08促进文献[29]
    Cs+BC 1.89
    Sr(Ⅱ)NaClO4 50 mM; 9.5稳定、可移动72促进文献[21]
    Sr+BC  100
    U(Ⅵ)NaCl 1 mM; 7稳定、可移动45促进文献[16]
    U+BC  68
    Cs(Ⅰ)碳酸盐微咸水170 mM; 7.6不稳定0.03~0.29促进文献[30]
    Cs+BC0.27~0.31
    Ce(Ⅲ)17~41阻滞
    Ce+BC0.8~1.4
    U(Ⅵ)39~67阻滞
    U+BC23~40
    Eu(Ⅲ)NaCl 1 mM; 6.5稳定、可移动34促进文献[17]
    Eu+BC  78
    Eu(Ⅲ)NaCl 10 mM; 6.5不稳定68阻滞
    Eu+BC  18
    Eu(Ⅲ)NaCl 1 mM; 3.5不稳定58阻滞
    Eu+BC  39
    下载: 导出CSV 
    | 显示表格

    低离子强度、碱性电解质中,膨润土胶体促进了Cs(Ⅰ),U(Ⅵ)在饱和花岗岩碎屑和Sr(Ⅱ)在裂隙中迁移(图6(a))。去离子水中,Cs(Ⅰ)单独迁移时测得的活性度极低,而加入膨润土胶体后,部分Cs(Ⅰ)不可逆地和胶体吸附位点结合后由胶体快速迁移,活性度提高至2%。对胶体、Cs(Ⅰ)共迁移时的流出液分析,发现98%的Cs(Ⅰ)从胶体表面解吸出来滞留在裂隙表面[29]。离子强度50 mM、pH9.5的NaClO4电解质中,85Sr(Ⅱ) 单独迁移穿透曲线呈单峰,而和FEBEX膨润土胶体共迁移的穿透曲线呈现3个主峰,分别对应不同的阻滞因子和回收率,其中第一个峰值反映Sr吸附在移动胶体上不受阻滞地迁移,且吸附可逆(Rf~1;Re=1%);第二个峰值反映Sr以溶质态迁移,且和岩面相互作用(Rf~6;Re=8%);第三个峰值反映了Sr被岩面胶体所吸附(Rf~20;Re=40%)。Sr(Ⅱ)总回收率由单独迁移72%增至共迁移100% [21]。类似地,离子强度1 mM、pH 7的NaCl电解质中,GMZ膨润土胶体促进了U(Ⅵ)的迁移,最大穿透率由45%提高至68%[16]

    图  6  动态柱胶体-核素共迁移试验结果
    Figure  6.  Results of dynamic column co-migration tests

    高离子强度或酸性电解质中,膨润土胶体阻滞了Ce(Ⅲ)在碳酸盐裂隙和Eu(Ⅲ)在饱和石英砂中的迁移(图6(b))。在170 mM、pH 7.6、含高浓度Ca2+地下水中,膨润土胶体不稳定,Cs(Ⅰ)在胶体表面的吸附性、迁移性均较差,胶体加入对Cs(Ⅰ)回收率影响较小;失去迁移性的胶体会抑制Ce(Ⅲ)自身形成的真胶体进入和穿透裂隙,并与Ce(Ⅲ)胶体絮凝形成大的异质胶体集合沉积下来,使Ce(Ⅲ)的回收率大幅下降;U(Ⅵ)与Ca2+,CO32-形成络合物,抑制其吸附在胶体和裂隙表面,胶体对其迁移影响较小[30]。类似地,低离子强度、弱酸性电解质中,GMZ膨润土胶体促进了Eu(Ⅲ)迁移,最大穿透率提高,而高离子强度或酸性电解质中胶体则阻滞Eu(Ⅲ)迁移[17]。此外,pH和流速会改变迁移介质对核素的吸附阻滞特性,从而影响胶体和核素共迁移作用。如pH为8时,MX80膨润土胶体促进Np(Ⅴ) 迁移,而pH为10时,高流速(0.8 mL/min)条件下胶体促进Np(V)迁移,低流速(0.3 mL/min)则阻滞迁移。花岗岩对Np(V)的吸附性随pH增大而增强,流速越大胶体越不容易被介质过滤[28]

    偶极子流场试验能够模拟复杂岩石裂隙系统,更接近真实工况,多用于原位试验。学者在GTS试验场的花岗岩剪切区开展了CRR、CFM原位偶极子流场试验,研究了多种核素共存时,膨润土胶体对其在裂隙剪切带迁移的影响规律。

    CRR、CFM原位偶极子流场试验如图7(a),CRR试验中注射流量为10 mL/min,提取流量为150 mL/min,迁移时间约1 h,而CFM试验通过减小提取流量(25 mL/min)来提高迁移时间[11, 31]。其中,原位偶极子流场试验中的水流速度比天然地下水流速大几个数量级,可能对试验结果带来一定影响[32]。因此,有学者设计了室内偶极子流场试验(图7(b)),在裂隙岩块中钻取若干和水平主裂隙相交的钻孔,其中注射钻孔和高效液相色谱注射泵相连,提取钻孔和馏分收集器相连,注射流量为0.25 mL/min [33]。偶极子流场试验过程及结束条件和柱试验基本一致。

    图  7  偶极子流场试验原理
    Figure  7.  Schematic diagram of dipole flow field tests

    CRR偶极子流场试验混合核素单独迁移以及和FEBEX膨润土胶体共迁移结果见表2图8[34]。由表2的回收率分析可知,在离子强度为1.2 mM,pH9.6的花岗岩地下水中,膨润土胶体对Am(Ⅲ)、Pu(Ⅳ)、U(Ⅵ))的迁移起到促进作用,对Np(Ⅴ)的迁移起到阻滞作用,而对Sr(Ⅱ)迁移影响较小。由表2的胶体分数分析可知,Am(Ⅲ),Pu(Ⅳ),Th(Ⅳ)单独迁移时其自身会水解成真胶体而加快迁移,穿透曲线峰值位置出现时间略早于示踪剂I。Np(Ⅴ)和U(Ⅵ)的穿透曲线基本一致,均存在明显的拖尾现象。Cs(Ⅰ)穿透曲线出现两个峰值,第一个峰值反映少量的Cs(Ⅰ)由胶体携带运输,阻滞因子较小,而第二个峰值反映大部分Cs(Ⅰ)的迁移受到迁移介质的强烈阻滞[10, 34]。膨润土胶体对Sr(Ⅱ)迁移影响很小,这是因为胶体对Sr(Ⅱ)的吸附可逆,容易被裂隙表面吸附。此外,室内偶极子流场试验结果同样表明,Sr(Ⅱ)和膨润土胶体在共迁移过程中会从胶体表面解吸出来,在高浓度胶体过滤作用下Sr(Ⅱ)的回收率降低,迁移被阻滞[33]

    表  2  偶极子流场胶体-核素共迁移结果[34]
    Table  2.  Co-migration results in dipole flow field [34]
    试验内容胶体分数/%阻滞因子回收率/%
    1单独迁移I(Ⅰ)-13101100
    Sr(Ⅱ)-850387
    Am(Ⅲ)-2436~580.8830
    Pu(Ⅳ)-2425~580.8718
    Th(IV)-23220~300.86
    Np(Ⅴ)-2370~10197
    U(Ⅵ)-2380~12153
    2共迁移I(Ⅰ)-131+BC0192
    Sr(II)-85+BC02.588
    Am(Ⅲ)-241+BC990.955
    Pu(IV)-244+BC840.977
    Th(IV)-232 + BC940.955
    Np(Ⅴ)-237 + BC0~1178
    U(Ⅵ)-233 + BC6198
    Cs-137+BC80.81;12170
    下载: 导出CSV 
    | 显示表格
    图  8  原位偶极子流场共迁移试验结果[34]
    Figure  8.  Results of dipole flow field migration [34]

    膨润土胶体-核素的共迁移作用取决于胶体的吸附-解吸作用以及胶体的过滤作用,并与迁移介质表面特性有关。

    溶解性和移动性较差的核素在胶体吸附作用下成为放射性假胶体,提高了其在水中的有效浓度[35];同时,尺寸排阻、电荷排阻效应以及水动力弥散作用使得核素的基质扩散和孔隙表面对其吸附反应等阻滞作用被减弱,核素由胶体携带着快速迁移[22]图9(a))。然而,只有胶体表面存在核素的强吸附且不可逆吸附位点,即核素在胶体表面的解吸作用较弱时,携带核素的胶体才有可能在岩石裂隙中进行长距离迁移[36]

    图  9  胶体-核素共迁移作用机理[22, 35]
    Figure  9.  Mechanism of colloid and nuclide co-migration[22, 35]

    膨润土胶体对不同核素具有不同的吸附-解吸特性[37-38],见图10。FEBEX膨润土胶体对Cs(Ⅰ)的第一次解吸分配系数小于吸附分配系数且基本不随时间变化,而随解吸次数增加,解吸分配系数增大且随接触时间增加而提高。分析认为,部分Cs(Ⅰ)在胶体表面的吸附可逆,位于阳离子交换位点上,容易进入胶体表面,为弱吸附。胶体对U(Ⅵ)的第一次解吸分配系数略大于吸附分配系数,且随解吸次数增加而增大,但基本不随时间变化,反映出胶体对U(Ⅵ)吸附不是单纯阳离子交换,且U(Ⅵ)在胶体上不存在动态结合位点。由此可见,Cs(Ⅰ)和U(Ⅵ)在膨润土胶体上的吸附不完全可逆[37]。研究发现,在胶体和Cs(Ⅰ)共迁移中,Cs(Ⅰ)极易从胶体表面解吸出来滞留在裂隙表面,导致胶体对迁移的促进作用较小[39-40],从胶体表面解吸出来的核素有可能进一步吸附在裂隙填充物上。地下水同时存在胶体、核素、裂隙填充物时,胶体对Am(Ⅲ),Pu(Ⅳ),Th(Ⅳ)呈现强吸附特性,而对Np(Ⅴ),U(Ⅵ)和Tc(Ⅶ)呈现弱吸附特性。由式(2)计算解吸速率,发现Pu(Ⅳ)和Th(Ⅳ)的解吸速率明显小于Am(Ⅲ)[38]

    图  10  核素在膨润土胶体的吸附-解吸动力学[37]
    Figure  10.  Adsorption-desorption kinetics of bentonite colloids [37]
    C(t)=C0ekt (2)

    式中,C(t)为t时刻核素浓度(mol/L),C0为核素初始浓度(mol/L),k为解吸速率(h-1),t为时间(h)。

    悬浮移动的胶体受介质环境影响,发生絮凝、过滤或沉降,从水相中去除。这种胶体移动性降低,由移动液相转变为不可移动固相的过程称为过滤作用。胶体过滤作用取决于胶体电荷和粒径特性、迁移介质固体表面化学特性以及地下水流速和化学条件[21]

    膨润土胶体在裂隙中迁移,受地下水流速、离子强度以及pH影响显著,是因为胶体稳定性和移动性受地下水动力和化学条件影响。一方面,地下水流速降低,胶体迁移滞留时间增加,胶体和裂隙面碰撞几率增大,导致介质表面对胶体吸附作用增强。另一方面,地下水离子强度、pH以及阳离子都会改变胶体表面电荷分布及粒径[13, 15],进而改变稳定性和移动性。当地下水离子强度较高且超过胶体聚沉值时,胶体絮凝形成较大集合体,在孔隙收缩处形成多颗粒架桥现象,或发生重力沉降附着在裂隙表面。酸性地下水中,胶体边缘电荷由负电荷转变为正电荷,附着在带负电荷的介质表面。阳离子吸附在胶体表面,进而对胶体表面电荷起到屏蔽作用,增加胶体和介质表面的静电吸引作用。吸附核素的胶体因过滤作用滞留在介质表面,从而阻滞核素迁移[17]图9(b))。

    介质阻塞效应多指多孔介质中迁移的胶体,在表面过滤、附着以及内阻塞等机理作用下[41],滞留在介质中并阻塞孔隙,导致介质渗透性下降[35]。当较大浓度胶体和核素共迁移时,会引发介质阻塞效应。高浓度核素离子会促使胶体絮凝沉积在介质表面,同时高浓度胶体则会在孔隙处聚集,絮凝成大颗粒,从而堵塞迁移通道[42-43]。阻塞效应还受到胶体和迁移介质的相对电荷特性[44]、胶体粒径和介质孔径的相对大小[45]以及地下水流速[43]、离子强度和pH[42]等影响,本质是胶体稳定性和移动性的改变。阻塞效应在裂隙介质中的膨润土胶体和核素共迁移研究中鲜有报道。

    膨润土胶体和核素共迁移穿透曲线可采用双重渗透率和双重孔隙介质模型模拟[46-47],描述特定裂隙-孔隙介质下膨润土胶体对核素迁移的促进现象。

    双重渗透率模型是由Ibaraki等[46]基于胶体促进水相污染物在离散断裂多孔介质的迁移机理建立的二维数值模型,能够描述对流-弥散水相溶质在岩石基质和裂隙中的迁移,胶体在裂隙中的迁移以及溶质在裂隙、移动和过滤胶体表面的吸附反应,主要控制方程为[27]

    θRfCt+qiCxixiθτDdCxj+θRfλdC=0, (3)
    qi=Kijhxij, (4)

    式中,θ为基质孔隙度,C为基质中溶质浓度,qi为基质中达西流量,τDd为有效弥散系数,λd为衰变常数,Kij为水力传导系数,h为水头。

    分布在水中、裂隙表面、移动胶体和过滤胶体表面的核素质量,通过下式描述各相间的交换过程:

    2b(Cft+qfCfxxDfCfx+λdCf+tSm)+2(tSs+tSσ)=0, (5)

    式中,2b为裂隙宽度,Cf为核素在裂隙中的浓度,qf为沿裂隙的达西通量,Df为核素弥散系数,Ss,Sm,Sσ为吸附在裂隙表面、移动和过滤胶体表面的核素质量。Sm,Sσ分别采用下式描述:

    2b(Smt+qmSmxxDmSmx+λdSmtRcfSm+λfqmSm)=0, (6)
    2b(1bSσt+1bλdSσ1btRCfSσλfqmSm)=0, (7)

    式中,Dm为胶体的弥散系数,λf为胶体过滤系数(m-1),qm为胶体在裂隙中的流速。

    采用平衡或动力学吸附反应方程式描述核素在移动、过滤胶体表面的吸附反应:

    Sm=Kdm×Cf Sσ=Kdσ×Cf, (8)
    Smt=βm(Kdm×CfSm) Sσt=βσ(K×CfSσ), (9)

    式中,Kdm,K为移动、过滤胶体的吸附分配系数,βm,βσ为移动胶体、过滤胶体的吸附动力学速率常数。

    由式(8),(9)可知,双重渗透率模型模拟胶体和核素共迁移参数包括:①裂隙系统几何参数,如地下水流动通道、裂隙宽度、孔隙度、渗透系数等;②吸附、迁移相关参数,如吸附分配系数、动力学反应速率和胶体移动速度等[31]。该模型模拟室内Kunipia-F膨润土胶体和Cs(Ⅰ)的共迁移穿透曲线见图11,发现仅考虑平衡吸附或动力学吸附反应时,模拟结果与试验值偏差较大,而同时考虑动力学吸附和胶体过滤的结果比较准确[27]。此外,对CRR现场试验中FEBEX膨润土胶体和Am(Ⅲ)的共迁移模拟发现,同时考虑动力学吸附以及胶体尺寸、电荷排阻效应时,双重渗透率模型计算值和试验数据较吻合[31]

    图  11  双重渗透率模型模拟结果[27]
    Figure  11.  Simulated results of dual permeability model[27]

    双重孔隙介质模型由Baek和Pitt建立,以描述胶体促进核素在饱和岩石基质和单一平面裂隙系统中的迁移行为,可同时考虑胶体过滤和再活化作用[47-48]。模型中核素和移动、过滤胶体间的交换反应Qm,Qσ

    Qm=km[CcmKdmCfCm], (10)
    Qσ=Kσ[CcσKdσCfCσ] (11)

    而胶体和裂隙表面的交换反应,包括过滤和再活化作用为

    Qc=CcσtλfucCcmbRmbCcσ (12)

    胶体促进的核素迁移过程为

    Cft+uCfzDf2Cfz2=Qm1bQfr1bQσ1bQfpλdCf=nλdnCfn (13)

    式中Ccf,Ccm,C为吸附在裂隙、移动和过滤胶体的浓度(mol/m3);Cf,Cm,Cσ为溶解在水中、吸附在移动和过滤胶体上的核素浓度(mol/m3);km,kσ为核素在移动和过滤胶体上的吸附动力学反应速率(a-1);Rmb为胶体的再活化系数(a-1);uc,u为胶体、核素的平均速度(m/a);Qfr,Qfp为核素和岩石裂隙、基质孔隙的交换过程;n为核素相关指数。

    双重孔隙介质模型中,吸附分配系数Kd、吸附动力学反应速率k、过滤系数λf和再活化系数Rmb都对模拟结果准确性起着关键作用。图12为模拟CFM项目GTS试验场偶极子流场试验FEBEX膨润土胶体和核素同系物共迁移的穿透曲线,结果表明,核素从胶体上解吸后快速吸附到裂隙填充材料上,胶体在裂隙表面的过滤、再活化等动态过程可解释试验结果[48]。此外,岩石也可吸附核素,因此考虑胶体和岩石表面吸附作用的LANL模型用于CFM胶体和核素共迁移模拟,可考虑不同吸附、解吸速率以及胶体核素结合的时间效应,更好地模拟试验结果[12]

    图  12  双重孔隙介质模型模拟结果[48]
    Figure  12.  Simulated results of double-porosity model [48]

    近年来,各国学者从工程实际出发,在高放射性废物深地质处置库膨润土胶体的吸附迁移特性以及膨润土胶体和核素的共迁移试验、模拟、共迁移作用机理等方面都进行了一定的探索研究,可得到3点结论。

    (1)膨润土胶体的吸附迁移特性受胶体浓度、地下水离子强度和pH影响显著。一定胶体浓度范围内,随胶体浓度增加,胶体吸附性提高而迁移性下降;随离子强度升高和pH下降,胶体的吸附迁移性下降。然而已有成果难以评价膨润土胶体对核素的吸附性及胶体的迁移能力。

    (2)动态柱和偶极子流场试验研究同时关注了移动胶体对核素迁移的促进作用以及过滤胶体对核素迁移的阻滞作用。然而,胶体的阻滞作用依赖于地下水化学环境,不利于胶体稳定的地下水环境中是否存在较高浓度的胶体是值得关注的问题,目前尚缺乏原位膨润土胶体携带核素远距离迁移的研究成果。

    (3)膨润土胶体和核素共迁移效果依赖胶体的吸附-解吸作用以及胶体的过滤作用,裂隙介质的孔隙、电荷特性的变化值得进一步关注。

    (4)双重渗透率模型和双重孔隙介质模型能够模拟一定试验条件下膨润土胶体和某些核素的共迁移穿透曲线,研究成果集中在国外FEBEX膨润土胶体和核素的共迁移模拟,模拟的裂隙系统多为单一平行裂隙,鲜少考虑核素的竞争吸附效应。国内已有试验成果集中在饱和多孔介质,尚未开展膨润土胶体-核素在真实岩石裂隙中的室内及原位迁移试验,缺乏模拟所需的相关参数,如裂隙几何参数、胶体解吸动力学参数及过滤系数等,因此未见关于国内GMZ膨润土胶体和核素的共迁移穿透模拟成果的报道。

    总体而言,对处置库膨润土胶体和核素迁移的研究尚处于试验研究阶段,理论研究成果也不够成熟。试验方面,需要深入开展膨润土胶体对多种核素的竞争吸附-解吸试验,考察胶体迁移过程中对裂隙岩石介质潜在的阻塞效应;同时,开展大尺度原位迁移试验,评估胶体-核素的远距离共迁移能力。理论方面,进一步考虑多核素共存时胶体表面的竞争吸附效应;可建立复杂裂隙网络系统并合理建模,改进已有模型,综合评价胶体对核素迁移的作用。

  • 图  1   脱水泥饼实物图

    Figure  1.   Physical picture of dewatered slurry cake

    图  2   大尺寸固体废弃物压缩-渗透装置

    Figure  2.   Large-scale solid waste compression-permeation experimental system

    图  3   泥饼压缩-渗透系列试验方案

    Figure  3.   Experimental scheme for compression-permeation of dewatered slurry cake

    图  4   试样分级加载时压缩量、孔隙比及饱和度变化曲线

    Figure  4.   Change of deformation, void ratio and saturation of samples during loading

    图  5   A试样

    Figure  5.   Sample A

    图  6   A试样在不同荷载条件下孔隙分布示意图

    Figure  6.   Pore distribution of Sample A under different loading conditions

    图  7   试样空气孔隙率和压实度与荷载的关系

    Figure  7.   Change of volume ratio of air and degree of compaction of samples against external loads

    图  8   试样孔隙比和压实度与荷载的关系

    Figure  8.   Curves of void ratio and degree of compaction of samples against external loads

    图  9   渗气系数与荷载的关系

    Figure  9.   Curve of gas permeability coefficient against external loads

    图  10   渗气系数与ea的关系

    Figure  10.   Curve of gas permeability coefficient against ea

    图  11   渗透系数和压实度与荷载的关系

    Figure  11.   Curves of permeability coefficient and degree of compaction of the sample against external loads

    表  1   试验所用泥饼的基本物理性质

    Table  1   Basic physical properties of dewatered slurry cake

    相对质量密度Gs塑性指数含水率/%初始干密度/(g·cm-3)最大干密度/(g·cm-3)最优含水率/%粒径组成/%
    d>0.0750.075≥d>0.005d≤0.005
    2.6613.5245.161.191.3829.840.9568.530.55
    下载: 导出CSV

    表  2   试验所用泥饼的矿物组分及含量

    Table  2   Mineral composition and content of dewatered slurry cake

    矿物成分石英云母方解石钙长石高岭石
    含量33%31%10%22%4%
    下载: 导出CSV

    表  3   分级加载过程中试样的固结系数

    Table  3   Consolidation coefficient of samples during graded loading

    荷载/kPa固结系数cv*/(mm2·s-1)
    A试样B试样
    50609.31806.78
    100140.85497.62
    200108.36115.81
    40038.6032.80
    8002.332.02
    下载: 导出CSV
  • [1] 郭卫社, 王百泉, 李沿宗, 等. 盾构渣土无害化处理、资源化利用现状与展望[J]. 隧道建设(中英文), 2020, 40(8): 1101-1112. https://www.cnki.com.cn/Article/CJFDTOTAL-JSSD202008004.htm

    GUO Wei-she, WANG Bai-quan, LI Yan-zong, et al. Status quo and prospect of harmless disposal and reclamation of shield muck in China[J]. Tunnel Construction, 2020, 40(8): 1101-1112. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JSSD202008004.htm

    [2] 丁建文, 洪振舜, 刘松玉. 疏浚淤泥流动固化处理与流动性试验研究[J]. 岩土力学, 2011, 32(增刊1): 280-284. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2011S1051.htm

    DING Jian-wen, HONG Zhen-shun, LIU Song-yu. Study of flow-solidification method and fluidity test of dredged clays[J]. Rock and Soil Mechanics, 2011, 32(S1): 280-284. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2011S1051.htm

    [3]

    BARENBLATT G I, ZHELTOV I P, KOCHINA I N. Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks[J]. Journal of Applied Mathematics, 1960, 24: 1286-1303.

    [4]

    DELAGE P, LEFEBVRE G. Study of the structure of a sensitive champlain clay and of its evolution during consolidation[J]. Canadian Geotechnical Journal, 1984, 21(1): 21-35. doi: 10.1139/t84-003

    [5]

    JUANG C H, HOLTZ R D. Fabric, pore size distribution, and permeability of sandy soils[J]. Journal of Geotechnical Engineering, 1986, 112(9): 855-868. doi: 10.1061/(ASCE)0733-9410(1986)112:9(855)

    [6]

    ROMERO E, DELLA VECCHIA G, JOMMI C. An insight into the water retention properties of compacted clayey soils[J]. Géotechnique, 2011, 61(4): 313-328. doi: 10.1680/geot.2011.61.4.313

    [7]

    ALONSO E E, VAUNAT J, GENS A. Modelling the mechanical behaviour of expansive clays[J]. Engineering Geology, 1999, 54: 173-183. doi: 10.1016/S0013-7952(99)00079-4

    [8]

    GENS A, ALONSO E E. A framework for the behaviour of unsaturated expansive clays[J]. Canadian Geotechnical Journal, 1992, 29(6): 1013-1032. doi: 10.1139/t92-120

    [9] 李光耀. 毛细阻滞型覆盖层微观-宏观水气传导特性及服役性能[D]. 杭州: 浙江大学, 2020.

    LI Guang-yao. Micro- Macro-Scale Water and Gas Transport Properties and Practical Performance of Capillary Barrier Covers[D]. Hangzhou: Zhejiang University, 2020. (in Chinese)

    [10]

    ROSQVIST H, DESTOUNI G. Solute transport through preferential pathways in municipal solid waste[J]. Journal of Contaminant Hydrology, 2000, 46(1/2): 39-60.

    [11]

    MCCREANOR P T, REINHART D R. Mathematical modeling of leachate routing in a leachate recirculating landfill[J]. Water Research, 2000, 34(4): 1285-1295. doi: 10.1016/S0043-1354(99)00243-2

    [12]

    GERKE H H, GENUCHTEN M T. A dual-porosity model for simulating the preferential movement of water and solutes in structured porous media[J]. Water Resources Research, 1993, 29(2): 305-319. doi: 10.1029/92WR02339

    [13] 杨益彪, 詹良通, 王顺玉, 等. 封顶覆盖黄土气体渗透特性测试及填埋气一维运移分析[J]. 岩土力学, 2015, 36(7): 1973-1980. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201507026.htm

    YANG Yi-biao, ZHAN Liang-tong, WANG Shun-yu, et al. Measurement of gas permeability of cover loess and one-dimensional analysis of landfill gas emission[J]. Rock and Soil Mechanics, 2015, 36(7): 1973-1980. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201507026.htm

    [14]

    TULI A, HOPMANS J W, ROLSTON D E, et al. Comparison of air and water permeability between disturbed and undisturbed soils[J]. Soil Science Society of America Journal, 2005, 69(5): 1361-1371. doi: 10.2136/sssaj2004.0332

    [15] 王卫华, 王全九, 樊军. 原状土与扰动土导气率、导水率与含水率的关系[J]. 农业工程学报, 2008, 24(8): 25-29. https://www.cnki.com.cn/Article/CJFDTOTAL-NYGU200808007.htm

    WANG Wei-hua, WANG Quan-jiu, FAN Jun. Relationship between air permeability, water conductivity and water content for undisturbed and disturbed soils[J]. Transactions of the Chinese Society of Agricultural Engineering, 2008, 24(8): 25-29. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-NYGU200808007.htm

    [16] 刘锦阳, 李喜安, 简涛, 等. 马兰黄土渗气率与饱和渗透系数的关系研究[J]. 水文地质工程地质, 2017, 44(6): 154-162. https://www.cnki.com.cn/Article/CJFDTOTAL-SWDG201706025.htm

    LIU Jin-yang, LI Xi-an, JIAN Tao, et al. A study of the relationship between gas permeability and saturated coefficient of permeability of the Malan loess[J]. Hydrogeology and Engineering Geology, 2017, 44(6): 154-162. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SWDG201706025.htm

    [17] 彭尔兴, 孙文博, 章定文, 等. 非饱和含砂细粒土的气体渗透特性研究[J]. 岩土力学, 2016, 37(5): 1301-1306, 1316. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201605012.htm

    PENG Er-xing, SUN Wen-bo, ZHANG Ding-wen, et al. Air permeability of unsaturated fine sandy soil[J]. Rock and Soil Mechanics, 2016, 37(5): 1301-1306, 1316. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201605012.htm

    [18] 苗强强, 陈正汉, 张磊, 等. 非饱和黏土质砂的渗气规律试验研究[J]. 岩土力学, 2010, 31(12): 3746-3750, 3757. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201012010.htm

    MIAO Qiang-qiang, CHEN Zheng-han, ZHANG Lei, et al. Experimental study of gas permeability of unsaturated clayey sand[J]. Rock and Soil Mechanics, 2010, 31(12): 3746-3750, 3757. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201012010.htm

    [19] 刘锦阳, 李喜安, 张凯旋, 等. 马兰黄土渗气性试验研究[J]. 水土保持学报, 2018, 32(1): 141-146. https://www.cnki.com.cn/Article/CJFDTOTAL-TRQS201801023.htm

    LIU Jin-yang, LI Xi-an, ZHANG Kai-xuan, et al. Experiment study on air permeability of malan loess[J]. Journal of Soil and Water Conservation, 2018, 32(1): 141-146. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TRQS201801023.htm

    [20] 土工试验方法标准:GB/T 50123—2019[S]. 2019.

    Standard for Geotechnical Testing Method: GB/T 50123—2019[S]. 2019. (in Chinese)

    [21] 土的工程分类标准:GB/T 50145—2007[S]. 2007.

    Standard for Engineering Classification of Soil: GB/T50145-2007[S]. 2007. (in Chinese)

    [22] 陈存礼, 张登飞, 张洁, 等. 等向应力条件下原状Q3黄土的渗气特性研究[J]. 岩土工程学报, 2017, 39(2): 287-294. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201702015.htm

    CHEN Cun-li, ZHANG Deng-fei, ZHANG Jie, et al. Gas permeability of intact Q3 loess under isotropic stresses[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(2): 287-294. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201702015.htm

    [23] 陈正汉, 谢定义, 王永胜. 非饱和土的水气运动规律及其工程性质研究[J]. 岩土工程学报, 1993, 15(3): 9-20. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC199303001.htm

    CHEN Zheng-han, XIE Ding-yi, WANG Yong-sheng. Experimental studies of laws of fluid motion, suction and pore pressures in unsaturated soil[J]. Chinese Journal of Geotechnical Engineering, 1993, 15(3): 9-20. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC199303001.htm

    [24] 姚志华, 陈正汉, 黄雪峰, 等. 非饱和Q3黄土渗气特性试验研究[J]. 岩石力学与工程学报, 2012, 31(6): 1264-1273. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201206024.htm

    YAO Zhi-hua, CHEN Zheng-han, HUANG Xue-feng, et al. Experimental research on gas permeability of unsaturated Q3 loess[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(6): 1264-1273. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201206024.htm

    [25] 谢定义. 非饱和土土力学[M]. 北京: 高等教育出版社, 2015.

    XIE Ding-yi. Soil Mechanics for Unsaturated Soils[M]. Beijing: Higher Education Press, 2015. (in Chinese)

    [26]

    FREDLUND D G, RAHARDJO H. Soil Mechanics for Unsaturated Soils[M]. New York: John Wiley & Sons, 1993.

    [27] 万智. 细粒土的压实特性与潮湿地区公路路基填筑控制技术研究[D]. 长沙: 中南大学, 2010.

    WAN Zhi. Research on Compaction Characteristics of Fine-Grained Soil and Compaction Control Technology in Wet Area[D]. Changsha: Central South University, 2010. (in Chinese)

    [28]

    Great Britain, Department of Transport, GBW Office. Specification for Road and Bridge Works[M]. 5th ed. London: HMSO, 1976.

    [29] 林鹏, 许镇鸿, 徐鹏, 等. 软土压缩过程中固结系数的研究[J]. 岩土力学, 2003, 24(1): 106-108, 112. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200301025.htm

    LIN Peng, XU Zhen-hong, XU Peng, et al. Research on coefficient of consolidation of soft clay under compression[J]. Rock and Soil Mechanics, 2003, 24(1): 106-108, 112. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200301025.htm

    [30] 朱前林, 李小春, 魏宁, 等. 多孔介质中气泡尺寸对流动阻力的影响[J]. 岩土力学, 2012, 33(3): 913-918. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201203042.htm

    ZHU Qian-lin, LI Xiao-chun, WEI Ning, et al. Influence of bubble size on fluid resistance in porous media[J]. Rock and Soil Mechanics, 2012, 33(3): 913-918. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201203042.htm

    [31]

    FREDLUND D G, RAHARDJO H, FREDLUND M D. Air Flow through Unsaturated Soils[M]. New York: John Wiley & Sons, Ltd, 2012.

    [32]

    BLACKWELL P S, RINGROSE-VOASE A J, JAYAWARDANE N S, et al. The use of air-filled porosity and intrinsic permeability to air to characterize structure of macropore space and saturated hydraulic conductivity of clay soils[J]. Journal of Soil Science, 1990, 41(2): 215-228.

    [33]

    LOLL P, MOLDRUP P, SCHJONNING P, et al. Predicting saturated hydraulic conductivity from air permeability: application in stochastic water infiltration modeling[J]. Water Resources Research, 1999, 35(8): 2387-2400.

    [34] 水利水电工程地质勘察规范:GB 50487—2008[S]. 2008.

    Code for Engineering Geological Invcstigation of Water Resources and Hydropower: GB 50487—2008[S]. 2008. (in Chinese)

  • 期刊类型引用(5)

    1. 唐清枫,徐真,孙亚楼,张奇,谭琪,潘多强,吴王锁. 黏土矿物-生物复合胶体影响放射性核素迁移行为与机理研究进展. 核化学与放射化学. 2025(01): 1-16 . 百度学术
    2. 吴鹏,王驹,凌辉,周志超,段佳欣,李南,段先哲. 甘肃北山新场深部地下水中铀的赋存形态及其影响因素的地球化学模拟研究. 原子能科学技术. 2024(05): 1007-1016 . 百度学术
    3. 田云婷,谭凯旋,李咏梅,李春光,唐治鹏,李小杰. 还原条件下膨润土胶体稳定性研究. 南华大学学报(自然科学版). 2024(03): 47-52 . 百度学术
    4. 邹威燕,周书葵,段毅. 不同地层中放射性核素迁移模型. 有色金属(冶炼部分). 2022(01): 79-87 . 百度学术
    5. 朱帅润,李绍红,何博,吴礼舟. 改进的Picard法在非饱和土渗流中的应用研究. 岩土工程学报. 2022(04): 712-720 . 本站查看

    其他类型引用(7)

图(11)  /  表(3)
计量
  • 文章访问数: 
  • HTML全文浏览量:  0
  • PDF下载量: 
  • 被引次数: 12
出版历程
  • 收稿日期:  2021-01-27
  • 网络出版日期:  2022-12-02
  • 刊出日期:  2021-09-30

目录

/

返回文章
返回