颗粒级配和含水率对含团粒泥炭土导热特性的影响

Effect of Grading and Water Content on Thermal Conductivity of Natural Peat Aggregated Soils

摘要

地表泥炭土是天然保温材料，对多年冻土保护具有重要意义.自然状态下，泥炭土通常含有粗团粒，且级配具有空间变异性.然而，级配对其导热特性的影响鲜有报道.采用非稳态探针法，对泥炭土开展了66组导热系数试验，以探究粗粒含量C_c、粒径比R和含水率w等因素对导热特性的影响.研究表明：在同一相对密度下，泥炭土的导热系数K随C_c的增加先增大后减小，在C_c=50%时取到最大值；粗粒主导时，随R的增大而增大，细粒主导时，R变化对其影响不显著.在相同孔隙比下，由于颗粒间接触热阻和配位数随C_c和R的变化，导热系数随着C_c的增加而增大；细粒主导时随R的增大而增大；粗粒主导时，则随R增大而降低.此外，泥炭土的导热系数随w增加而显著增大，其增长率随试样中团粒内部和团粒间孔隙中的水分分布而变化.

Abstract

Peat soil is a natural insulation material, which has a significant effect on the underlying frozen soil layer. Natural peat soil usually exists in the form of aggregates, and its particle size distribution has a spatial variability. However, the effect of gradings on the thermal conductivity is rarely reported. In this work, 66 laboratory tests were conducted on the peat soil to investigate the effects of coarse content, particle size ratio and water content on its thermal conductivity by a non-steady-state probe method.The results show that under the same relative density, the thermal conductivity of peat soil first increases and then decreases with the increase of coarse content, and reaches the peak value when the coarse content reaches 50%. For the peat soil with coarse dominated structure, the thermal conductivity increases distinctly with the increase of particle size ratio. However, the effect of particle size ratio reduces in the fine dominated structure. For a given void ratio, the thermal conductivity increases with the rising coarse content due to the thermal resistance and the particle contacts which are affected by particle size distribution; it increases with increasing particle size ratio in fine grain dominated structure but decreases in coarse grain dominated structure. In addition, the thermal conductivity of peat soil increases significantly with increasing water content, and the growth rate varies with the internal pores and the water distribution in the pores.

Graphical abstract

关键词

泥炭土 / 导热系数 / 团粒 / 粒径分布 / 含水率 / 工程地质学.

Key words

peat soil / thermal conductivity / aggregate particle / particle size distribution / moisture content / engineering geology

0 引言

在气候变暖和工程活动的共同影响下，全球冻土退化日益加剧（刘桂民等，2023；汪涛等， 2023；罗栋梁等，2024）.多年冻土地区表层的富有机质泥炭土因其低导热性和冻融条件下的导热差异性而成为下卧冻土的天然保温材料（程国栋和何平，2001； 陈槐等，2021；Bian et al.， 2024）.因此在冻土工程施工过程中，对地表植被和泥炭土层的扰动和破坏，使下卧冻土在夏季失去了良好的绝缘隔热层，从而导致其夏季温度升高，活动层加深，进而引起冻土融区面积扩大，冻土强度降低（Wu et al.， 2016），给寒区环境和冻土工程带来了众多负面影响（如地表植被退化、路基的不均匀沉降、基础塌陷和山体滑坡等）（Ding et al.， 2021；Che et al.， 2023；Qin and Zhou， 2023）.由此可见探究泥炭土的导热特性对于减缓冻土退化和寒区工程防灾减灾具有重要意义.

以往对泥炭土导热系数的研究主要集中在分析其受有机质含量、矿物含量、干密度、含水量、温度、冻融循环次数等的影响.He et al.（2022）通过室内试验发现土体的导热系数随有机质含量增加而显著减小.Mustamo et al.（2019）对按不同比例混合的泥炭土与矿物颗粒混合物进行导热系数测定，发现混合物导热系数随矿物含量的增加逐渐增大.冷毅飞（2011）和桂跃等（2023）将干密度和含水率作为土样变量，研究发现土样的导热系数随干密度和含水率增加而增大.杨宗维（2018）和Xu et al.（2020）研究了温度对泥炭土导热系数的影响，结果显示泥炭土导热系数在正温度范围内随温度降低呈线性下降，在负温度范围内呈指数上升.然而，已有研究表明，在基质吸力和有机质胶结的共同作用下，泥炭土常以一定大小的团粒形式存在（Schweizer et al.， 2019）.此外，由于团粒的双孔隙结构（内部孔隙和粒间孔隙）使团粒在气候与植被因素及外力作用更易发生破碎，从而使团粒的粒径分布产生空间变异性（Shi et al.， 2018； Xiao et al.， 2018；南富森等，2023；谭文峰等，2023）.而大量研究结果表明粒间结构变化（如孔隙比、颗粒形状、相对密度和粒径分布）是引起颗粒材料导热系数变化的主要原因（Lee et al.， 2017； 南博文， 2019）.同时，团粒作为双孔隙结构，其孔隙分布对土体导热系数的影响也受到了研究人员的广泛研究（Carminati et al.， 2008； Peth et al.， 2010； 蔡国庆等， 2020）.然而含团粒泥炭土粒径分布和孔隙分布对其导热系数的影响尚不明确.

本研究通过将天然状态的泥炭土进行充分的烘干和颗粒筛分处理，获得了不同粒径的团粒.在干燥状态下，对16种不同级配的团粒混合料进行了导热系数测试，以探究颗粒级配和孔隙分布对泥炭土导热系数的影响及相关作用机制.此外，本文还探究了在不同含水率下，泥炭土导热系数的变化规律.

1 试验材料和仪器

1.1　试验材料

本研究所用的材料是取自中国内蒙古呼伦贝尔根河市永久冻土区的表层泥炭土，原状泥炭土含水率高，土质均匀，呈黑色或棕褐色.为获取用于试验的泥炭土团粒，本研究首先将高含水率泥炭土进行自然风干.随着含水率的减小，大颗粒土块和团粒含量逐渐增加.为保持团粒的结构完整，将风干土样过2 mm筛网，去除粒径大于2 mm的土块，并将过筛颗粒继续烘干形成干燥状态的泥炭土团粒，烘干过程中为了防止泥炭土中有机质的分解，烘箱温度设置为65℃.如图1所示，干燥状态下泥炭土呈黑色或深灰色，且具有明显的宽级配颗粒材料特征（Shi et al.， 2021）（如图2中混合料级配曲线所示）.采用四分法取代表性土样配制悬浊液过0.075 mm洗筛以确定天然泥炭土材料的净砂含量，结果表明：土样中粒径大于0.075 mm的净砂含量低于总质量的7%.因此，用于试验的泥炭土主要由黏土颗粒在基质吸力和有机质胶结的共同作用下形成的团粒组成（Schweizer et al.， 2019）.泥炭土的基本物理特性参照GB/T 50123-2019标准进行测定，试验结果如表1所示.参照GB50021-2001标准，土样的有机质含量为17.42%，可将土体归类为弱泥炭质土；参照GB50007-2011标准，土样塑性指数I_P为15.29，可将土体归类为粉质黏土.

图2中泥炭土混合料粒径分布曲线表明：团粒粒径在0.075~0.250 mm范围内含量相对较低.同时，团粒结构的稳定性在此范围两侧表现出显著差异（刘亚龙等， 2023），因此在本文中以0.2 mm粒径为分界点，将泥炭土团粒体划分为粗粒组和细粒组.将泥炭土团粒体过孔径为0.2 mm筛网进行筛分，利用干筛法（刘亚龙等， 2023）分别确定如图2中细粒组粗粒组粒径分布曲线.粗粒和细粒的基本物理特性指标如表2所示.

1.2　试验仪器与原理

本研究采用非稳态探针（NSSP）法（即瞬态法）测量泥炭土的导热系数.与稳态方法相比，该方法的测试时间更短且更加准确（Lin et al.， 2023）.导热系数的测定装置如图3a所示，导热系数测定设备选用Hukseflux公司生产的TP08小型瞬态导热系数探针，探针针长70 mm，直径为1.2 mm，包括一个热电偶（热电偶一结点位于距针尖约15 mm处，另一结点位于底部）和一条加热丝，测量范围为0.1~6.0 W·m^-1·K^-1，在20 ℃下的测试精度为±（3%+0.02）W·m^-1·K^-1.其余设备包括：一个施加恒流的电源系统，一个数据记录仪，一台计算机.在每次测试中，用绝缘胶布将顶盖与试样筒之间的缝隙密封，将探针从顶盖预留的孔内插入试样中（图3b），并通过电源系统对加热丝施加直流电以产生热量，使针头的温度升高，随后热量从针头传递到周围的颗粒.在测试过程中，通过热电偶监测针内温度，并使用数据采集器连接探针的传感器部分，记录传感器信号，以确定时间和温差.

图4给出了温度-时间关系的典型结果，热导率K计算公式为：

K=\frac{Q}{\mathrm{4}\mathrm{\pi }}\frac{\mathrm{\Delta }\mathrm{l}\mathrm{n}t}{\mathrm{\Delta }T},（1）

式中：Q为单位长度内输入的线热源功率W/m；t是从加热开始的持续时间（s）；T是针头在相应时间的温度（K）.K是测量介质的热导率W·m^-1·K^-1.

由公式（1）分析可知，试样的热导率越大，热耗散速率将会越高.因此，热电偶检测到的温度升高速率会越低.温度时间轴的初始阶段会受针-土样耦合的影响，而在长时间加热后，试样边界会对温度变化产生扰动.根据ASTM D5334-14标准，热导率是从温度对数时间图的线性中心部分获得，如图4中Δlnt部分所示.测试程序和计算方法均参考ASTM D5334-14执行.

2 试验方案与试验步骤

2.1　试验方案

本文主要对含团粒泥炭土进行两个系列的试验研究.第一系列试验主要研究干燥条件下泥炭土颗粒级配和粒径比变化对导热系数的影响.

首先，为了得到不同级配的泥炭土试样，将上述筛分得到的粗粒组与细粒组团粒按不同比例均匀混合，其中混合料中粗颗粒含量（以下简称为粗粒含量C_c）是表征颗粒级配的参数指标.已有的试验结果（吴琪等， 2018； Basiri Parsa and Maleki， 2023）表明，对于具有明显颗粒粒径差异的粗细颗粒混合料，必定存在一个阈值粗粒含量，总结已有研究发现，其值一般高于20%，因此本文在等间距线性分布粗粒含量取值（0%，25%，50%，75%，100%）的基础上，进行局部内插加密，以获得粗细粒混合料不同的颗粒接触状态及阈值粗粒含量.粗粒含量设置为0%，25%，40%，50%，65%，75%，85%，100%，试样的粒径分布曲线如图5a所示.

其次，在研究粒径比（R=粗粒组d₅₀/细粒组d₅₀）变化对泥炭土导热系数的影响中，将细粒组、粗粒组代表性土样分别通过0.15、0.125、0.074、0.01 mm标准筛组和2、1、0.71、0.50、0.45、0.30、0.20 mm标准筛组得到如图5b、5c中粒径比R=10.2的细粒组、粗粒组团粒粒径分布曲线.在细粒组颗粒主导的结构中通过改变粗粒组团粒的最大粒径来改变粗粒组d₅₀，从而改变试样的粒径比（R），如图5b所示.为了使试验结果更显著，本文选用粗粒含量为40%的土样.与此相对应，在粗粒组颗粒主导的结构中，选用粗粒含量为75%的土样，通过改变细粒组团粒的最小粒径来改变粗粒组d₅₀，从而改变试样的粒径比（R），如图5c所示.具体的试验方案及试样的基本物理参数如表4所示.干土试样的最大、最小干密度参考GB/T 50123-2019标准进行测定.通过试验测定试样的颗粒破碎指数，结果表明：粗粒含量分别为75%，85%，100%的试样在相对密度为75%的情况下进行导热系数测定后，颗粒破碎指数分别为3.69%，5.32%，9.00%.因此在本研究中分别测定试样在相对密度（D_r）为30%，50%，75%情况下的导热系数.

在第二系列的研究中，以粗粒含量为50%的泥炭土击实试验确定的最大干密度（1.098 g/cm³）为基准，在压实度为100%，90%和80%对应的3个干密度，研究随含水率增加，泥炭土导热系数的变化规律.具体试验方案如表5所示.

2.2　试验步骤

在制备土样前，为了确保泥炭土处于干燥状态，将土体在65℃烘箱中烘干24 h后，密封冷却至室温并进行筛分.土样制备按以下步骤进行：

（1）将细粒与粗粒组团粒按照目标级配充分混合，以达到视觉上的均匀性.

（2）为配置满足目标含水率的土样，首先采用烘干法确定代表性土样的初始含水率.

（3）称取适量均匀混合的目标级配土体，基于初始含水率，选用喷雾器加水拌合至目标含水率，闷料12 h后测定土样的实际含水率（当制备干燥状态试样时，可忽略步骤2、3）.

（4）根据试样的目标干密度和含水率称取适量土体制样，为了保证试样的均匀性，采用欠压法（Lin et al.， 2023）将土体分5层（D_r=30%，50%，75%）填充在高100 mm，内部直径为60 mm测试容器中，以获得所需的相对密度，前一层的干密度比后一层的干密度小1%.

（5）试样制备完成后，对于非饱和试样，为了使试样中水分分布均匀，将试样闷料12 h后进行导热系数测定.

所有测试均在室内条件下进行，为了控制恒温环境，试验过程中关闭试验室门窗，防止空气对流影响传感器冷端的温度，并根据温度传感器的灵敏度，对热传导率进行校正，从而将环境温度的影响降至最低.同时，内部直径为60 mm的测试容器满足ASTM D5334-14要求，避免试验过程中边界效应的发生.导热系数测定按以下步骤进行：

（1）将TP08探针从容器顶盖垂直插入试样中心位置.

（2）连接试验设备，设置相关程序，开始试验.每次测试的时间为3 min，每隔15 min以上进行重复测试（Xu et al.， 2022；Yang et al.， 2022），以确保试样的温度恢复到初始值.每个试样进行3次测试，试样导热系数的最终值为重复测试的平均值，且测量值的标准差在3%以内.

（3）测试完成后，将试样用调土刀逐层从试筒中取出，并从试样中下部取两份15~30 g土料，平行测定试样含水率.试验结果表明：干燥状态试样在试验结束后其含水率均低于0.3%，因此水的影响不在此系列的研究范围内；不同含水率的试样在试验结束后其含水率与目标含水率误差均小于0.5%.

3 试验结果与分析

3.1　级配对最大和最小孔隙比的影响

图6给出了试样的最大孔隙比（e_max）和最小孔隙比（e_min）随粗粒含量的变化关系：当粗粒含量较低时，粗粒分散在细粒基质中，如图7a所示，试样的最大和最小孔隙比（e_max，e_min）随粗粒含量增加逐渐减小；当最大及最小孔隙比取得最小值时，粗粒彼此接触，粒间孔隙完全被细粒填充，如图7b所示，此时粗粒含量称为“过渡粗粒含量”（吴琪等， 2018）；在此之后，细粒分散在粗粒基质中，如图7c所示，试样的最大和最小孔隙比（e_max，e_min）逐渐增大.因此，基于试样的最大、最小孔隙比曲线可以确定泥炭土粗细粒混合料的过渡粗粒含量为50%.

同时，由于含团粒试样的孔隙比包括粒间孔隙与团粒内部孔隙，且团粒的内部孔隙随粒径的增加而逐渐增大（Usowicz et al.， 2013； He et al.， 2022）.对比试样最大孔隙比随粗粒含量变化规律与团粒内部孔隙的变化规律可以发现：试样最大孔隙比随粗粒含量变化主要受团粒间孔隙分布的影响.

基于粗细粒混合料粒间接触理论，当粗粒含量较低时，如图7a所示，细粒基质彼此接触形成土体骨架，粗颗粒完全被包裹在其中，此时土体为细粒主导结构（Shi et al.， 2020）；随着其含量增加，粗粒继续分散在细粒基质中，但由于颗粒粒径的增加，部分细粒及其粒间孔隙逐渐由粗粒取代，从而使最大和最小孔隙比（e_max，e_min）逐渐减小.直到粗粒含量达到50%的过渡粗粒含量时，如图7b所示，大部分粗颗粒彼此接触，部分细粒与粗粒接触共同形成土体骨架，同时粗粒间孔隙完全由细粒填充，此时最大最小孔隙比（e_max，e_min）取得最小值.随着粗粒含量的继续增加，如图7c所示，土体骨架完全由粗颗粒构成，细粒填充在粗粒间孔隙中，此时土体为粗粒主导结构（Shi et al.， 2020）.随着粗粒含量的继续增加，细颗粒无法充分填充粗粒间的孔隙，而使试样的最大和最小孔隙比逐渐增大.然而通过对击实试验完成后原粗粒含量为75%，85%，100%的土样进行筛分，试验结果表明：不同初始状态下颗粒破碎指数分别达到了4.73%，6.38%，11.18%.正因如此，在击实过程中，粗颗粒在外力作用发生颗粒破碎，破碎后的颗粒更好地填充了粗粒形成的孔隙，从而降低了粗粒含量增加对最小孔隙比的影响（施勇等， 2023）.

图8分别给出了在细粒和粗粒主导结构中，试样最大和最小孔隙比（e_max，e_min）随粒径比变化的关系.由图8a可知，在细粒主导结构中，粒径比的变化对试样最大和最小孔隙比无显著影响.这可能与细粒基质的粒径分布和孔隙结构有关.如前所述，在细粒主导结构中，粗粒分散在细粒基质中而彼此无接触，此时混合料的力学特性完全由细粒控制（Shi et al.， 2020）.因此，不改变粗粒含量而增大粗颗粒最大粒径从而增大粒径比不会对试样的最大和最小孔隙比产生显著影响.通过试验数据分析可知，当粒径比从3.8增加到8时，最大孔隙比逐渐减小了0.059；当粒径比从8增加到10时，最大孔隙比略微增大了0.041.这主要因为，随着粒径比的增大，粗粒组最大粒径逐渐增大，粗粒内部孔隙逐渐增加（Usowicz et al.， 2013）.同时，增大的粗粒逐渐取代细粒及粒间孔隙.在其双重作用下，试样孔隙比略微降低.而随着粒径比的继续增大，如图9a~9b团粒粒间孔隙变化图表明：当粗粒之间彼此接触并与细粒基质之间发生分层时，粗粒间产生的封闭孔隙增大了试样的孔隙比（Shi and Zhao， 2020）.

如图8b所示，在粗粒主导结构中，最大孔隙比和最小孔隙比随粒径比的增大而逐渐减小.这是因为：在粗粒主导结构中，细粒被限制在粗粒基质形成的孔隙中，如图10a所示，减小细粒的最小粒径，使细粒能够更好地填充粗粒间的大孔隙，如图10b所示，从而使试样的孔隙比减小.

3.2　粗粒含量对导热系数的影响

不同相对密度下土体导热系数随粗粒含量变化规律如图11所示.泥炭土试样的导热系数随着相对密度的增大而增大.这是因为随着相对密度的增大，试样的孔隙比减小，由于土颗粒的导热系数远大于空气的导热系数（唐盼盼等， 2019），孔隙比的减小使得颗粒间的有效传热面积增大，传热路径数量增加，从而使得试样的导热系数增大（Xiao et al.， 2018； Zhang et al.， 2020； He et al.， 2022）.在一定的相对密度下，泥炭土的导热系数随着粗粒含量的增加先逐渐增大；当粗粒含量接近50%（过渡粗粒含量）时，导热系数达到最大值；在此之后，导热系数随粗粒含量的增加逐渐减小.如图12中给出了相 对密度为50%的条件下，随着粗粒含量增加，试样孔隙比和导热系数的变化规律.对比结果表明：泥炭土导热系数随粗粒含量的变化与孔隙比密切相关，这表明孔隙比的变化是引起泥炭土导热系数随粗粒含量变化的关键因素.

正如之前的研究人员所指出的（Chen， 2008； Xiao et al.， 2018； Zhang et al.， 2020； He et al.， 2022）：级配对导热系数的影响还归因于颗粒之间的传热机制，该机制与孔隙率无关.因此，本节分析了在给定孔隙比下泥炭土随粗粒含量变化的导热系数变化规律，如图13所示，可以发现：在相同孔隙比下，随着粗粒含量的增加，导热系数逐渐增大.这是因为：如图14所示，在相同孔隙比下，当热流通过被测介质时，沿相同的传热路径，粒间接触点数目随颗粒粒径的增加而减少，这有效降低了颗粒接触产生的热阻（Zhang et al.， 2020； He et al.， 2022），从而提高了土体的传热性能.同时，在粗粒含量低于75%时，试样的不均匀系数随着粗粒含量的增加而增大，细粒逐渐填充大颗粒孔隙，增大了颗粒间的配位数，进而增加了试样的传热路径，使导热系数增大（Xiao et al.， 2018）.同时，随孔隙比的增大，泥炭土的导热系数减小（Chen， 2008； He et al.， 2022），二者呈近似的线性关系，R²均大于0.99，具有较高的拟合度.以上研究表明：泥炭土导热系数受粗粒含量的影响主要与孔隙比和颗粒间传热机制的变化有关.

3.3　粒径比对导热系数的影响

3.3.1　细粒主导结构中粒径比对导热系数的影响

为了分析粒径比对泥炭土导热系数的影响，图15给出了细粒主导结构中，不同相对密度下土体导热系数随粒径比的变化规律.结果表明：相同粒径比下，泥炭土的导热系数随着相对密度的增加而增大，这与之前的研究结果一致.在相同的相对密度下，细粒主导结构中，泥炭土导热系数随粒径比的变化而略有变化，这表明对于细粒主导结构，粒径比的变化对导热系数的影响较小.但在总体趋势上表现出随粒径比的增大先增大后减小的规律：当粒径比从3.8增加到8的过程中，导热系数缓慢增大；当粒径比从8增加到10时，导热系数略微减小.这可以解释为：在细粒主导结构中，粗粒最大粒径的增大在引入更多的颗粒内部孔隙的同时取代了部分细粒基质及粒间孔隙，使试样的孔隙比减小，传热路径增加，导热系数增大.当粗粒最大粒径继续增大使彼此接触形成粒间封闭空隙（Shi and Zhao， 2020）时，如图9所示，试样孔隙比增大，传热路径减少，导热系数减小.同时，由图16可知，相同粒径比下，导热系数随孔隙比的增大而减小.这表明在细粒主导结构中，孔隙比的变化是引起泥炭土导热系数随粒径比变化的关键因素.

由图16也可以发现，在细粒主导结构中，在一定孔隙比下，导热系数随着粒径比的增大逐渐增大；这可以解释为：随着粒径比的增加，试样的不均匀系数逐渐增大，土体颗粒间的接触面积和配位数随着均匀性的降低而增加（Xiao et al.，2018），如图9所示，使土体导热系数增大.同时，颗粒间传热产生的热阻受颗粒粒径的影响（Zhang et al.， 2020； He et al.， 2022）：在细粒主导结构中，随粒径比的增加，粗粒粒径增大有助于降低颗粒间热阻，从而使导热系数增大.

3.3.2　粗粒主导结构中粒径比对导热系数的影响

图17给出了粗粒主导结构中，不同相对密度下土体导热系数随粒径比的变化规律.由图可知：在粗粒主导结构中，导热系数随粒径比增大逐渐增大.这可以解释为：在粗粒主导结构中，随着粒径比的增大，细粒的粒径逐渐减小，更好地填充了粗粒间的孔隙，如图10所示，使试样的孔隙比减小，增加了试样的有效传热路径，从而使导热系数增加（Xiao et al.， 2018）.此外，由图18可知，粗粒主导结构中，在相同孔隙比下，泥炭土导热系数随粒径比的增大而减小.其主要原因是：如前所述，在粗粒主导结构中，粗颗粒彼此接触形成土体骨架，细粒填充在粗粒间的孔隙中，随着粒径比的增大，细粒的粒径变小，填充相同孔隙时所需的颗粒数目增加，如图10所示，在传热路径中因颗粒接触而产生的热阻增大，则试样的导热系数越小（南博文， 2019）.

3.4　含水率对导热系数的影响

图19给出了不同干密度下，泥炭土团粒混合料导热系数与含水率的关系曲线.由图可知，泥炭土导热系数随含水率的增加而显著增大.其主要原因在于水的导热系数为0.6 W·m^-1·K^-1，这远大于孔隙中空气的导热系数（0.023 W·m^-1·K^-1）（唐盼盼等， 2019）.同时，在一定的干密度下，随含水率的增加，导热系数增长率不断改变.这主要与试样中团粒内部孔隙和团粒间孔隙中的水分分布有关（Usowicz et al.， 2013）.试验结果表明：当含水率低于10%时，泥炭土导热系数随含水率的变化较小.这主要是因为泥炭土中大量的有机质储存在团粒内部的黏土颗粒的矿物表面，当含水率较低时，孔隙水首先进入团粒内部孔隙中，在黏土颗粒表面发生水化作用，随后在黏土颗粒接触处形成水膜，并随着含水率的增加，水膜相互重叠，团聚体间孔隙逐渐饱和，如图20中步骤1所示.此时，团粒内部导热系数由于孔隙水的增加而逐渐增大，但在团粒间孔隙区域的有效传热路径较少，因此导热系数增长率较低（Côté and Konrad， 2005； Peth et al.， 2010； 褚召祥等， 2020）.随着含水率的不断增加，孔隙水逐渐在团粒间接触点处以水膜的形式存在，如图20中步骤2所示，水膜的出现极大地增加了粒间孔隙中的有效传热路径，从而使导热系数迅速增加.当含水率继续增加，团粒间孔隙内水膜相互重叠，如图20中步骤3所示，导热系数进一步增加，但增长速率有所降低.当土体接近饱和时，如图20中步骤4所示，导热系数逐渐趋于稳定值（Carminati et al.， 2008； Lu and Dong， 2015）.此外，需要指出的是：在含水率逐渐增加的过程中，泥炭土中的团粒结构由于水力作用而发生崩解破坏或因基质吸力的减小和可溶性有机质的溶解而逐渐发生分解（Carminati et al.， 2008； Peth et al.， 2010），使试样中的孔隙分布趋于均匀（蔡国庆等， 2020），同样会引起导热系数增长率的改变.

4 结论

本研究在干燥状态下，将泥炭土中自然存在的团粒体划分为粗细两个粒组，通过控制泥炭土中初始级配的粗粒含量和粒径比，研究团粒的粒径分布和孔隙分布对泥炭土导热系数的影响.此外在泥炭土粗粒含量为50%的初始级配下，通过控制土体的含水率和干密度，进一步研究了泥炭土导热系数随含水率的变化规律.得出以下结论：

（1）含团粒泥炭土的粒间接触状态随粗粒含量的增加而变化：随着粗粒含量的增加，泥炭土由粗粒分散在细粒基质中的细粒主导结构逐渐转变为细粒填充在粗粒间孔隙中的粗粒主导结构.区分细粒与粗粒主导结构的阈值被称为过渡粗粒含量，此时试样在各相对密度下达到最密实状态.本研究中泥炭土团粒混合料的过渡粗粒含量为50%.

（2）泥炭土导热系数在试样的孔隙比和颗粒间传热机制的共同影响下随级配分布的变化而变化：在相同相对密度下，泥炭土的导热系数随粗粒含量的增加先增大后减小，当粗粒含量达到50%时取到最大值.在粗粒主导结构中，泥炭土的导热系数随粒径比的增大而增大，但在细粒主导结构中，粒径比变化对导热系数的影响不显著.在相同孔隙比下，随着粗粒含量的增加导热系数逐渐增加，在细粒主导结构中导热系数随粒径比的增大而增加；而在粗粒主导结构中，导热系数随粒径比增大而降低，

（3）在相同干密度下，泥炭土的导热系数随含水率增加而显著增大，其增长率随试样中团粒内部孔隙和团粒间孔隙中的水分分布而逐渐变化：随着含水率的增加，孔隙水首先填充团粒内部孔隙，导热系数增长率较小；当团粒内部孔隙达到饱和后，孔隙水在土颗粒形成水膜，粒间传热路径显著增加，从而使导热系数迅速增大；随着含水率的继续增加，水膜逐渐重叠，粒间孔隙逐渐饱和，导热系数增长率降低，数值逐渐趋于稳定.

参考文献

[1] Basiri Parsa, S., Maleki, M., 2023. Factors Affecting Small-Strain Shear Modulus of Sand-Silt Mixture Considering Different Moisture Contents. Iranian Journal of Science and Technology, Transactions of Civil Engineering, 47(1): 479-490. https://doi.org/10.1007/s40996-022-01018-1

[2] Bian, X., Ren, Z. L., Zeng, L. L., et al.,2024. Effects of Biochar on the Compressibility of Soil with High Water Content. Journal of Cleaner Production, 434: 140032. https://doi.org/10.1016/j.jclepro.2023.140032

[3] Cai, G.Q., Wu, T.C., Wang, Y.N., et al., 2020. Model of the Microstructure Evolution of Unsaturated Compacted Soils with Double-Pore Structure. Rock and Soil Mechanics, 41(11): 3583-3590 (in Chinese with English abstract).

[4] Carminati, A., Kaestner, A., Lehmann, P., et al., 2008. Unsaturated Water Flow across Soil Aggregate Contacts. Advances in Water Resources, 31(9): 1221-1232. https://doi.org/10.1016/j.advwatres.2008.01.008

[5] Che, L. N., Zhang, H. H., Wan, L. H., 2023. Spatial Distribution of Permafrost Degradation and Its Impact on Vegetation Phenology from 2000 to 2020. Science of the Total Environment, 877: 162889. https://doi.org/10.1016/j.scitotenv.2023.162889

[6] Chen, H., Wu, N., Wang, Y.F., et al., 2021. A Historical Overview about Basic Issues and Studies of Mires. Scientia Sinica (Terrae), 51(1): 15-26 (in Chinese).

[7] Chen, S. X., 2008. Thermal Conductivity of Sands. Heat and Mass Transfer, 44(10): 1241-1246. https://doi.org/10.1007/s00231-007-0357-1

[8] Cheng, G.D., He, P., 2001. Linearity Engineering in Permafrost Areas. Journal of Glaciolgy and Geocryology, 23(3): 213-217 (in Chinese with English abstract).

[9] Chu, Z.X., Zhou, G.Q., Rao, Z.H., et al., 2020. Predicting Correlation and Evolution Mechanisms of the Effective Thermal Conductivity of Granular Geomaterials. Chinese Journal of Rock Mechanics and Engineering, 39(2): 384-397 (in Chinese with English abstract).

[10] Côté, J., Konrad, J. M., 2005. A Generalized Thermal Conductivity Model for Soils and Construction Materials. Canadian Geotechnical Journal, 42(2): 443-458. https://doi.org/10.1139/t04-106

[11] Ding, Y. J., Mu, C. C., Wu, T. H., et al., 2021. Increasing Cryospheric Hazards in a Warming Climate. Earth-Science Reviews, 213: 103500. https://doi.org/10.1016/j.earscirev.2020.103500

[12] Gui, Y., Xie, Z.P., Gao, Y.F., 2023. Influence and Mechanism of Organic Matter on Thermal Conductivity of Cohesive Soil. Rock and Soil Mechanics, 44(S1): 154-162(in Chinese with English abstract).

[13] He, H., He, Y., Cai, G. J., et al., 2022. Influence of Particle Size and Packing on the Thermal Conductivity of Carbonate Sand. Granular Matter, 24(4): 117. https://doi.org/10.1007/s10035-022-01277

[14] Lee, C., Suh, H. S., Yoon, B., et al., 2017. Particle Shape Effect on Thermal Conductivity and Shear Wave Velocity in Sands. Acta Geotechnica, 12(3): 615-625. https://doi.org/10.1007/s11440-017-0524-6

[15] Leng, Y.F., 2011. Experimental Study on Physical and Mechanical Properties and Numerical Analysis of Temperature Field of Permafrost in Sino-Russian Oil Pipeline(Dissertation).Jilin University, Changchun (in Chinese with English abstract).

[16] Lin, J. Z., Shi, X. S., Zeng, Y. W., et al., 2023. Estimating the Thermal Conductivity of Granular Soils Based on a Simplified Homogenization Method. Cold Regions Science and Technology, 211: 103855. https://doi.org/10.1016/j.coldregions.2023.103855

[17] Liu, G.M., Zhang, B., Wang, L., et al., 2023. Permafrost Region and Permafrost Area in Globe and China. Earth Science, 48(12): 4689-4698 (in Chinese with English abstract).

[18] Liu, Y.L., Wang, P., Wang, J.K., 2023. Formation and Stability Mechanism of Soil Aggregates: Progress and Prospect. Acta Pedologica Sinica, 60(3): 627-643 (in Chinese with English abstract).

[19] Lu, N., Dong, Y., 2015. Closed-Form Equation for Thermal Conductivity of Unsaturated Soils at Room Temperature. Journal of Geotechnical and Geoenvironmental Engineering, 141(6): 04015016. https://doi.org/10.1061/(asce)gt.1943-5606.0001295

[20] Luo, D. L., Liu，J., Chen, F. F., et al., 2024. Research Progress and Prospect of Transition Zone in Permafrost. Earth Science, 49(11): 4063-4081 (in Chinese with English abstract).

[21] Mustamo, P., Ronkanen, A. K., Berglund, Ö., et al., 2019. Thermal Conductivity of Unfrozen and Partially Frozen Managed Peat Soils. Soil and Tillage Research, 191: 245-255. https://doi.org/10.1016/j.still.2019.02.017

[22] Nan, B.W., 2019. Experimental Study on the Influence of Particle Gradation and Shape and Thermal Characteristics of Particles on Thermal Conductivity of Soil (Dissertation).Chongqing University, Chongqing (in Chinese with English abstract).

[23] Nan F. S., Li Z. X., Zhang X. P., et al., 2023. Particle Size Fractal Characteristics of Soils in Desert-Steppe Transition Zone along the Northern Bank of Yellow River Basin in Lanzhou. Earth Science, 48(3): 1195-1204 (in Chinese with English abstract).

[24] Peth, S., Nellesen, J., Fischer, G., et al., 2010. Non-Invasive 3D Analysis of Local Soil Deformation under Mechanical and Hydraulic Stresses by μCT and Digital Image Correlation. Soil and Tillage Research, 111(1): 3-18. https://doi.org/10.1016/j.still.2010.02.007

[25] Qin, C., Zhou, J., 2023. On the Seismic Stability of Soil Slopes Containing Dual Weak Layers: True Failure Load Assessment by Finite-Element Limit-Analysis. Acta Geotechnica, 18(6): 3153-3175. https://doi.org/10.1007/s11440-022-01730-2

[26] Schweizer, S. A., Bucka, F. B., Graf-Rosenfellner, M., et al., 2019. Soil Microaggregate Size Composition and Organic Matter Distribution as Affected by Clay Content. Geoderma, 355: 113901. https://doi.org/10.1016/j.geoderma.2019.113901

[27] Shi, X. S., Herle, I., Muir Wood, D., 2018. A Consolidation Model for Lumpy Composite Soils in Open-Pit Mining. Géotechnique, 68(3): 189-204. https://doi.org/10.1680/jgeot.16.p.054

[28] Shi, X. S., Nie, J. Y., Zhao, J. D., et al., 2020. A Homogenization Equation for the Small Strain Stiffness of Gap-Graded Granular Materials. Computers and Geotechnics, 121: 103440. https://doi.org/10.1016/j.compgeo.2020.103440

[29] Shi, X. S., Zhao, J. D., 2020. Practical Estimation of Compression Behavior of Clayey/Silty Sands Using Equivalent Void-Ratio Concept. Journal of Geotechnical and Geoenvironmental Engineering, 146(6): 04020046. https://doi.org/10.1061/(asce)gt.1943-5606.0002267

[30] Shi, X. S., Zhao, J. D., Gao, Y. F., 2021. A Homogenization-Based State-Dependent Model for Gap-Graded Granular Materials with Fine-Dominated Structure. International Journal for Numerical and Analytical Methods in Geomechanics, 45(8): 1007-1028. https://doi.org/10.1002/nag.3189

[31] Shi, Y., Jia, X.L., Lü, G.E., et al., 2023. Determination of Maximum and Minimum Void Ratios of Calcareous Sand Considering Various Influence Factors. Journal of Ground Improvement, 5(4): 293-298 (in Chinese with English abstract).

[32] Tan, W.F., Xu, Y., Shi, Z.H., et al., 2023. The Formation Process and Stabilization Mechanism of Soil Aggregates Driven by Binding Materials. Acta Pedologica Sinica, 60(5): 1297-1308 (in Chinese with English abstract).

[33] Tang, P.P., Xu, J., Lu, Y.H., 2019. Experimental Study on Effects of Water Content and Temperature on Thermal Conductivity of Unsaturated Soils. Journal of Disaster Prevention and Mitigation Engineering, 39(4): 678-683 (in Chinese with English abstract).

[34] Usowicz, B., Lipiec, J., Usowicz, J. B., et al., 2013. Effects of Aggregate Size on Soil Thermal Conductivity: Comparison of Measured and Model-Predicted Data. International Journal of Heat and Mass Transfer, 57(2): 536-541. https://doi.org/10.1016/j.ijheatmasstransfer.2012.10.067

[35] Wang, T.,Wang, X.Y., Liu, D., et al., 2023. The Current and Future of Terrestrial Carbon Balance over the Tibetan Plateau. Scientia Sinica (Terrae), 53(7): 1506-1516 (in Chinese).

[36] Wu, Q. B., Zhang, Z. Q., Gao, S. R., et al., 2016. Thermal Impacts of Engineering Activities and Vegetation Layer on Permafrostin Different Alpine Ecosystems of the Qinghai–Tibet Plateau, China. The Cryosphere, 10(4): 1695-1706. https://doi.org/10.5194/tc-10-1695-2016

[37] Wu, Q., Chen, G.X., Zhou, Z.L., et al., 2018. Experimental Investigation on Liquefaction Resistance of Fine-Coarse-Grained Soil Mixtures Based on Theory of Intergrain Contact State. Chinese Journal of Geotechnical Engineering, 40(3): 475-485 (in Chinese with English abstract).

[38] Xiao, Y., Liu, H. L., Nan, B. W., et al., 2018. Gradation-Dependent Thermal Conductivity of Sands. Journal of Geotechnical and Geoenvironmental Engineering, 144(9): 06018010. https://doi.org/10.1061/(asce)gt.1943-5606.0001943

[39] Xu, X. T., Zhang, W. D., Fan, C. X., et al., 2020. Effects of Temperature, Dry Density and Water Content on the Thermal Conductivity of Genhe Silty Clay. Results in Physics, 16: 102830. https://doi.org/10.1016/j.rinp.2019.102830

[40] Xu, Y. S., Zhou, X. Y., Sun, D. A., et al., 2022. Thermal Properties of GMZ Bentonite Pellet Mixtures Subjected to Different Temperatures for High-Level Radioactive Waste Repository. Acta Geotechnica, 17(3): 981-992. https://doi.org/10.1007/s11440-021-01244-3

[41] Yang, Y. L., Zhang, T., Reddy, K. R., et al., 2022. Thermal Conductivity of Scrap Tire Rubber-Sand Composite as Insulating Material: Experimental Investigation and Predictive Modeling. Construction and Building Materials, 332: 127387. https://doi.org/10.1016/j.conbuildmat.2022.127387

[42] Yang, Z.W., 2018. Experimental Study on Thermal Conductivity of Typical Soils in Northeastern Inner Mongolia during Freezing/Thawing Process (Dissertation).Inner Mongolia University, Hohhot (in Chinese with English abstract).

[43] Zhang, X. R., Kong, G. Q., Wang, L. H., et al., 2020. Measurement and Prediction on Thermal Conductivity of Fused Quartz. Scientific Reports, 10(1): 6559. https://doi.org/10.1038/s41598-020-62299-y

基金资助

国家自然科学基金项目(52278346)

深圳市科技计划项目(KCXFZ20211020163816023)

河套-深圳-香港科技创新合作区项目(HZQB-KCZYB-2020083)