文献参数: 白一兵, 张紫昭, 徐元鹏, 等.冻融循环条件下伊犁地区滑坡黄土湿陷特性及微观机理[J].水土保持通报,2026,46(2):32-42. Citation:Bai Yibing, Zhang Zizhao, Xu Yuanpeng, et al. Collapsibility characteristics and microscopic mechanisms of loess in Ili region under freeze-thaw cycles [J]. Bulletin of Soil and Water Conservation,2026,46(2):32-42.
黄土是一种兼具多孔隙与弱胶结特性的特殊性沉积物
[1],具有强烈的水敏性和湿陷性。在周期性冻融作用下,黄土经历多次冻结-融化过程,造成土体孔隙中液态水与固态水的周期性相变,进而引发土体颗粒间胶结作用减弱、孔隙结构重塑等,最终导致黄土边坡出现非均匀变形特征,严重影响边坡的安全性,极易引发工程灾害
[2]。同时,冻融循环会破坏黄土的多孔结构,导致融雪、降雨等大量地表水渗入坡体,使土体含水率升高,进而显著增强黄土的湿陷性
[3]。因此,开展黄土在冻融循环条件下湿陷特性的研究对于边坡安全稳定性具有重要的理论和实际意义。
目前,有多名学者针对冻融循环因素研究黄土的湿陷特性,并基于扫描电镜等仪器对黄土微观结构特征进行研究。冻融循环是一种对于土体的强风化作用。郑郧等
[4]、宁俊等
[5]通过试验得出冻融循环极大地改变了土体的湿陷性及微观结构。肖东辉等
[6]研究表明,冻融循环会促使黄土孔隙中的大颗粒逐步分解为小颗粒,不仅改变了黄土的孔隙率,还对其湿陷性产生了间接影响。李国玉等
[7]试验表明,在反复冻融作用下,黄土呈现出干密度减小,孔隙比增大的特征,同时水分会持续增加并向中上区域聚集,这些因素共同导致了湿陷变形量的上升。在湿陷特性方面,早在20世纪30年代,学术界便针对黄土湿陷展开系统性探索,研究重点集中在湿陷假说构建、湿陷系数与起始压力的关联、湿陷变形量测算及湿陷性评价方法等,这些研究为工程建设的科学规划、地质灾害的风险评估以及相关防治措施的实践应用提供了可靠的理论基础与技术保障,在相关领域实践中发挥了不可替代的关键作用
[8-10]。在理论层面,针对黄土湿陷机制的研究已形成多元体系,从早期的毛管假说、溶盐假说,到后续的胶体不足说、水膜楔入说,再到欠压密理论等,各类学说相互补充,不断完善,共同推动了对黄土湿陷本质认知的深化;江美英等
[11-14]提出水与外力是导致黄土湿陷性的外部因素,而黄土自身独特的架空孔隙结构及物质组成,则是影响其湿陷性的内部因素。杨鑫
[15]对原状黄土的增湿变形特性规律展开分析后发现,随着增湿含水率的升高,其增湿湿陷变形呈现逐渐减小的趋势。高英等
[16]研究了不同增湿含水率条件下微结构特征与湿陷变形的关联,结果表明当上覆压力增大时,强湿陷性黄土湿陷系数的增长斜率更为陡峭,且初始含水率较低的强湿陷性黄土,其变形发展速度更快。Liu Zhen等
[17],Li Ping等
[18], Li Xian等
[19]通过定量分析湿陷过程中颗粒形态、颗粒定向性及孔隙大小、孔隙面积等参数的变化规律,并结合物质成分分析结果,提出了黄土湿陷机理。随着扫描电镜等仪器出现,人们开展了对黄土微结构的研究,土体的微观结构对更好地了解土体性质具有重要意义。学界普遍认为,黄土的湿陷机制可从其自身结构特征切入解释,包括骨架颗粒的形态、接触关系与排列方式、胶结状态的差异以及孔隙的尺寸大小与形态分布等,这一认知逐步发展,最终形成了阐释黄土湿陷机制的结构学说
[20]。赵鲁庆等
[21]从微观角度出发,研究了土粒形态、连接方式、排列形式及土粒排列、形态、尺度分布等微观结构参数与冻融次数的关系。徐璐等
[22]研究发现,黄土湿陷过程会导致孔隙分形维数升高,同时黄土的孔隙形态也变得更为复杂。陈果等
[23]对比了黄土湿陷前后的微结构特征,结果显示,湿陷作用使黄土骨架颗粒的接触形态产生明显改变。本文以伊犁地区喀拉海依苏滑坡黄土为研究对象,基于其基本物理性质,开展冻融循环下的室内湿陷试验,分析湿陷系数与不同初始含水率以及冻融次数间的关系;结合扫描电镜试验,从定性及定量角度分析土体在不同含水率、不同循环次数条件下湿陷前、后黄土微观结构的变化特征和变化规律,揭示黄土湿陷特性,为滑坡灾害防治提供理论依据。
1 材料与方法
1.1 研究区概况
研究区位于新疆伊犁地区新源县阿热勒托别镇喀拉海依苏村,属典型的温带大陆气候。该区域属于大型深层黄土滑坡,平均厚度约40 m,滑坡隐患体所处斜坡地貌单元为剥蚀堆积的圆顶低山,形态近似半圆形,其上发育有多级拉张裂缝,裂缝多呈弧形近东西向排列,高程1 276~1 470 m,斜坡体坡向320°,相对高差194 m,纵向长度300.4 m,横向宽度50.4 m,面积9.80×104 m2。研究区大面积被第四系所覆盖,按时代由老至新划分为第四系上更新统风积层(Q3eol)、第四系滑坡堆积层(Q4dl)及第四系全新统冲洪积层(Q4apl)。
1.2 区内黄土基本特性
试验土样采自坡顶位置,土样呈黄褐色且质地相对均匀。取样深度为3 m,取样时间为2025年3月。为避免土样在运输环节因外界干扰发生结构破坏,采集后采用保鲜膜与泡沫膜进行多层包裹处理,再转运至实验室。随后,对采集的原状土样依次进行风干、粉碎处理,并通过室内基本物理性质试验测定相关参数,其主要物理指标详见
表1。研究区内土种为粉质黄土,其中,粉土占比为79.1%,砂土和黏土占比分别为10.08%和10.82%。土体中主要矿物有石英、白云母及绿泥石等(
表2)。
1.3 方法
1.3.1 冻融循环试验
试验以粉质黄土最优含水率(17.2%)为基准,设置6组梯度含水率(14.2%,15.2%,16.2%,18.2%,19.2%,20.2%),并采用最大干密度(1.83 g/cm³)的环刀试样。试验于2025年3—5月在新疆大学地质与矿业工程学院国家重点实验室进行。土样先过0.5 mm土工筛,烘干后通过喷壶少量多次补水,经充分搅拌均匀后装入保湿皿密封养护48 h。随后使用轻型击实法制备环刀试样,直径61.8 mm,高度20.00 mm。将成型后的试样置于恒温恒湿试验箱中,进行多次冻融循环测试。为防止土样水分流失,环刀内土样采用保鲜膜密封。设定冻结温度-20 ℃,融化温度20 ℃。冻结和融化阶段的持续时间均为12 h,以此还原土壤的冻结与融解过程;此外,温度波动度控制在±0.5 ℃,温度均匀度控制在±2 ℃,湿度维持在20%~98%,湿度偏差不超过±2.5%。冻融循环结束后,首先对环刀试样进行质量测定,要求质量误差不超过0.1%,同时含水率变化控制在0.3%以内,符合上述条件后方可进入下一轮冻融。试验共设置9次循环。为探究循环次数对湿陷特性的影响,分别选取0,1,3,5,7,9次循环作为典型节点,用于开展后续宏微观试验;每个典型次数均制备3个平行试样,保障数据可靠性与重复性。
1.3.2 湿陷试验
采用十六联固结仪开展湿陷试验(
表3)。采用双线法,将冻融循环完成后的试样,一个试样经历冻融循环试验后进行分级加荷,另一个冻融循环试样需要浸水饱和,将2个试样安装在十六联固结仪中,分级加载至
p=25 kPa,
p=50 kPa,
p=100 kPa,
p=200 kPa,
p=300 kPa,
p=500 kPa,
p=800 kPa,
p=1 000 kPa,待下沉稳定后,2个试样同时加至最后一级压力,下沉稳定变形量标准为0.01 mm/h。每间隔1 h读取并记录试验数据,采用饱和试样在浸水后的总变形量与非饱和试样在浸水前的压缩变形量之差表示湿陷变形。计算公式如下:
式中:为湿陷系数; 是在压力(kPa)作用下,浸水变形稳定后的土样高度与稳定前高度之差(mm); 为试样在未受压力作用下的原始高度(mm)。浸水变形前后高度差的增量 代表土壤湿陷量。
1.3.3 微观SEM试验
对冻融循环、湿陷后的试样进行扫描电镜试验(scanning electron microscope,简称SEM)。选取冻融循环0,1,9次以及含水率为14.2%,20.2%的湿陷前后试样,首先,将湿陷试验前后的试样自然风干,通过控制水分缓慢蒸发,最大限度减少对土体内部结构的扰动,同时确保试样满足扫描电镜试验所需的真空环境要求;随后,在土样的特定部位进行切样操作,并完成喷金处理,用场发射扫描电子显微镜Sigma 300提取相关微观参数。湿陷前后孔径分布、丰度、定向频率、分形维数计算公式见
表4。
2 结果与分析
2.1 冻融循环对压缩特性的影响
由
图1可知,在竖向荷载施加初期,土体孔隙体积快速减小,曲线表现出较大的斜率,表明该阶段土粒结构的不稳定性较高,颗粒发生快速重排引起变形。随荷载继续增加,孔隙内的水分与气体逐步被排出,孔隙体积缩小,斜率不断减小,压缩曲线呈现非线性递减趋势。相同冻融循环次数条件下,不同含水率粉土的压缩曲线总体变化趋势类似,表现为孔隙比随荷载增加而逐渐降低,主要是由于外加竖向荷载使土体内部微结构逐渐密实,进而孔隙比不断降低。在相同竖向荷载条件下,高含水率土样的压缩曲线始终高于低含水率土样,说明其初始孔隙比更大,含水率升高使土体内孔隙水含量增加,颗粒周围水膜增厚,从而削弱了颗粒间的摩擦作用,更易发生相对滑移,最终表现出更大的压缩变形。相反,低含水率土样中孔隙水较少,颗粒间摩擦阻力较大,限制了颗粒的运动,因此压缩变形较小,整体压缩性较低,变形程度也相对有限。
2.2 冻融循环对湿陷特性的影响
研究结果表明,在经历相同冻融循环次数后,不同含水率黄土的湿陷系数随垂直荷载增大均呈先上升后下降的趋势(
图2)。当轴向压力由初始值增加至300 kPa时,湿陷系数迅速增长,
-
P曲线斜率较大,表明该阶段湿陷变形发展显著;在300 kPa时湿陷系数达到峰值,说明试样在此荷载下湿陷性最为明显。随后轴向压力继续增大至1 000 kPa,湿陷系数逐渐减小,表明湿陷变形的增长幅度持续减缓,
-
P曲线逐渐趋于平缓。这种变化特征主要是由于水与荷载共同作用下黄土粒间孔隙发生塌陷,随后孔隙结构逐渐向镶嵌型转变,在持续压力作用下土体不断压密、孔隙被压缩,导致湿陷变形逐渐减弱。
由
图2可知,垂直荷载在300 kPa时,不同含水率的黄土在冻融循环过程中,湿陷系数峰值的下降幅度有所不同。含水率为14.2%时,在冻融循环9次后的湿陷系数峰值从0.094降至0.084,下降了10.6%;含水率为20.2%时,在冻融循环9次后的湿陷系数峰值从0.079降至0.076,下降了3.8%,表明含水率较低的黄土受冻融循环的影响更为显著,其湿陷性的降低幅度更大。同时,随着初始含水率的增加,黄土的湿陷系数峰值整体上表现出逐渐减小的趋势。在冻融循环为0次的初始状态时,湿陷系数由14.2%含水率下的0.094逐渐降低到20.2%含水率下的0.079(
表5),说明随着初始含水率的增加,黄土的初始湿陷性降低。
2.3 黄土湿陷微观机制研究
2.3.1 黄土微观图像定性分析
在初始含水量为14.2%的条件下,土样浸水前的SEM图像显示颗粒排列紧密,颗粒边界清晰,孔隙较小,颗粒间接触方式以紧密接触或少量黏土与胶结物质连接为主,表现出较为稳定的致密结构(
图3)。经浸水处理后,该土样颗粒之间的排列明显趋于松散,颗粒间隙显著增大,颗粒边界逐渐模糊,原有黏土与胶结物质发生软化乃至部分溶解,使得颗粒间原有的紧密连接遭到明显削弱甚至丧失,进而降低了土体的结构稳定性,表现出明显的膨胀与结构松散特征。
随着冻融循环次数的增加,粉质黄土的微观结构发生显著变化。含水率为14.2%的黄土,在冻融1次的样品中,土体中存在较多的大孔隙,骨架由较大的粉黏团粒和颗粒构成,颗粒间接触方式主要为点-点接触,结构松散且不稳定,此时土粒主要表现为较大的团粒状或片状结构,颗粒表面有较少细小碎屑,接触界面简单,易受水的作用而发生明显的湿陷变形(
图4a)。当冻融循环增加至9次时,原本的大孔隙已明显转化为中、小孔隙,大粒径的团粒状骨架颗粒受到反复冻胀挤压作用而破碎成中、小粒径颗粒,微观结构中碎屑和片状颗粒的含量显著增多,这些细小碎屑附着在粗骨架颗粒表面,使颗粒表面的粗糙度增加(
图4b)。同时,多次冻融循环使颗粒重新排列组合,颗粒接触更多表现为点-点接触或点-边接触,孔隙结构趋于细化,土体内可压缩的空间逐渐减小。这种微观结构的演变导致黄土的湿陷性明显减弱,湿陷系数随着循环次数的增加而逐渐减小。湿陷试验后,粉质黄土的微观结构进一步压密,与试验前进行对比,其骨架结构发生了显著的改变。冻融1次且经湿陷后的样品中,原先团粒状的骨架颗粒已经被水载破坏得更小,颗粒间的胶结物质被溶蚀或破坏,颗粒表面的碎屑物质增多,片状或次棱角状的粉粒可见,土体结构较未湿陷前更加密实(
图4c)。经过9次冻融循环并进行湿陷试验后,土颗粒间的排列更加紧凑,边-边接触数量显著增加,部分棱角较明显的粉粒进一步磨圆为偏圆状,孔隙进一步收缩为细小且分散的形态,虽然仍有少量团聚结构与架空孔隙残留,但颗粒的进一步重组压缩空间已十分有限,整体结构趋于稳定(
图4d)。多次冻融循环后,黄土的结构性湿陷特征几乎消失,湿陷系数在反复循环作用下逐渐趋于稳定,不再大幅变化。
2.3.2 黄土微观图像定量分析
(1) 孔径分布。孔径分布能够直接反映土体中不同尺度孔隙的数量与占比,是决定土体压缩性与湿陷潜力的关键结构特征。本研究根据黄土、软土孔径分类标准,将土壤孔隙直径大小分为4个区间:微孔隙(直径<2 μm),小孔隙(2~5 μm之间),中等孔隙(5~20 μm之间),大孔隙(直径>20 μm)。
由
图5可知,在同一含水率不同冻融循环次数条件下,黄土湿陷前后均以小孔隙最为丰富,其次为中孔隙,随着冻融循环次数的增加,黄土中不同孔径类型的孔隙数量占比发生变化。中孔隙数量从冻融循环1次湿陷前至冻融循环9次湿陷后,降低幅度达68.3%;大孔隙占比降幅为25.0%。相反,在经历冻融循环与湿陷作用后,小孔隙数量相对较高,下降幅度为18.4%。因此,研究区黄土湿陷变形的机理可概括为颗粒间黏粒胶结作用减弱,致使原有架空结构稳定性丧失而发生崩解;随之颗粒重新排列并充填孔隙,细颗粒逐渐进入大孔隙内,导致架空孔隙数量递减,孔隙类型也相应转为以镶嵌孔隙和胶结孔隙为主。
(2) 孔隙方向频率。颗粒或孔隙的长轴与同一平面法线方向的夹角称为方向角,方向频率的定义是在0°~180°的角度区间中,每个10°分段内方向角数量的占比
[29]。由
图6可知,在相同含水率条件下,不同冻融循环次数对研究区黄土的孔隙方向分布与数量产生显著影响。冻融循环1次时,湿陷前孔隙主要集中在0°~40°与130°~180°,表现出明显的方向性结构;而湿陷后,孔隙数量整体减少,且在0°,90°及180°减幅最为显著,表明水平方向与竖直方向的孔隙更易在湿陷过程中发生压缩闭合。随着冻融循环次数增加至9次,孔隙方向分布趋于均匀,方向性减弱。
(3) 孔隙丰度。土体孔隙的短轴与长轴的比值即为孔隙丰度,它反映了孔隙形态的几何特征。若孔隙丰度值较小,则孔隙的形状趋于长条形而非圆形;相反,如果孔隙丰度值较大,则孔隙形态趋向于圆形。
图7表明,14.2%含水率条件下,湿陷后孔隙丰度呈现增加趋势,冻融循环1次湿陷前后,孔隙丰度增加28%,冻融循环9次湿陷前后,孔隙丰度增加23%,意味着一部分长条形孔隙可能向圆形状转化,同时也有一小部分孔隙由长条状转为圆形。孔隙丰度增加表明孔隙形状趋向于圆形,形状更加规则。整体拟合曲线呈单调上升趋势,表明随着平均孔隙面积的增加,孔隙的方向丰度亦逐步增强,二者呈显著正相关关系,同时,反映出孔隙分布由无序向有序的方向增强。
(4) 孔隙分形维数。颗粒或孔隙的分形维数能够综合度量孔隙或颗粒边界轮廓的复杂程度,其值越大表明越复杂
[30]。本研究计算了土体中所有孔隙的分形维数,得到了14.2%含水率条件下不同冻融循环次数和湿陷前后的平均分形维数值(
图8)。
可以看出,分形维数随着冻融循环次数和湿陷作用的共同作用呈现显著上升趋势。A点的分形维数最低,表明初始状态下孔隙结构较为简单;B点相较A点有明显提升,说明1次冻融及其后的湿陷使颗粒结构破碎,孔隙分布趋于复杂。继续冻融至9次后,C点的分形维数进一步提高,说明多次冻融循环累积效应强化了孔隙微观扰动与细化趋势;D点分形维数达到最大值,表明经过反复冻融作用后,土体结构趋于碎化,孔隙形态高度不规则,该趋势表明冻融与湿陷过程在微观上显著改变了孔隙形态及其分布特征。
(5) 宏微观参数相关性。研究结果表明,当含水率为14.2%,冻融循环次数分别为0,1,9次时,湿陷系数与分形维数之间存在极强的负相关性(
表6),表明随着冻融循环次数的增加,土体内部颗粒结构趋于复杂化与不规则化,这一微观结构的重塑过程直接导致了宏观湿陷性的显著降低。
湿陷系数与大孔隙占比之间存在极强的正相关性,其p值为0.053,略高于0.05的显著性水平,可能是由于样本量有限所致,但也表明了强烈的相关趋势。该结果从定量上支持了经典的湿陷机理,即大孔隙起到了对黄土湿陷变形的主导作用。冻融循环作用通过破坏原有的架空孔隙结构,导致大孔隙数量减少,从而从根本上削弱了土体的湿陷性,使得宏观湿陷系数下降。
上述相关性分析表明,宏微观参数之间存在明确的内在联系,由此构建了冻融循环到微观结构重塑再到宏观湿陷性减弱的完整关联。
2.3.3 黄土微观湿陷机制
(1) 湿陷前微观机理。在含水率为14.2%的初始条件下,经过1次冻融循环处理的粉质黄土内部具有较多的大孔隙,土体骨架主要由较大的粉黏粒团聚体和颗粒构成。颗粒之间以点-点接触为主,整体结构松散且稳定性较差。此时土粒常呈较大的团粒状或片状,颗粒表面附着的细小碎屑很少,接触界面简单平直。这样的松散结构在遇水时极易发生湿陷变形,即水分渗入后削弱了颗粒间本已有限的胶结作用,导致土体骨架迅速塌陷压缩。与之相比,经过9次冻融循环后的粉质黄土微观结构明显趋于密实。反复的冻胀和融化使原本的大孔隙大部分转化为中、小孔隙;原先的大粒径团聚骨架在冻融应力反复挤压下破碎成中、小粒径的颗粒。显微结构中细碎的残片和片状粉粒含量显著增多,这些细小碎屑附着在粗骨架颗粒表面,增大了颗粒表面的粗糙度(
图9a)。多次冻融循环还使颗粒重新排列组合,颗粒间接触方式更多地表现为点接触或点-边接触,孔隙结构明显细化,可压缩的孔隙空间逐渐缩小。宏观上,由于微观结构的这种演变,黄土的湿陷性被显著削弱,随着冻融循环次数的增加,湿陷系数逐渐减小,土体抗湿陷能力提高。
(2) 湿陷后微观机理。在经过湿陷试验后,1次冻融循环的样品其微观结构发生了变化。原先团粒状的土骨架在水的冲击和荷载作用下被破坏得更细小,颗粒之间的胶结物质被大量溶蚀或破坏。破坏后的骨架颗粒表面碎屑明显增多,土粒形态上出现片状或次棱角状的粉粒。与湿陷前相比,此时土体结构更为密实紧致,孔隙体积减少,说明湿陷过程使得1次冻融后的黄土内部结构受到额外改造,削除了部分残余的大孔隙和不稳定结构。
经历9次冻融循环并进行湿陷试验后的样品,其微观结构表现出更强的稳定性和紧密性。湿陷过程使已经重组多次的颗粒进一步挤紧排列,颗粒间的边-边接触数量显著增加,架空孔隙和松散结构几乎被压垮殆尽(
图9b)。一些原本棱角较为锋利的粉粒在反复摩擦中被磨圆,变为偏圆形,从而改善了颗粒的嵌合程度和接触界面。试验后孔隙进一步收缩为细小且分散的形式,尽管仍可见极少量残留的团聚结构或局部架空孔隙,但颗粒进一步重组压缩的空间已非常有限,整体结构趋近稳定。由于多次冻融循环已大幅破坏并重构了土的初始结构,黄土的结构性湿陷特征在反复冻融后几乎消失。宏观上表现为湿陷系数在多次循环作用下趋于稳定,不再出现大的变化,说明土体经过反复冻融与湿陷过程后已达到相对稳定的状态。
3 结 论
(1) 冻融循环作用使黄土反复发生冻胀和融沉,随着冻融循环次数的增加,各个区间含水率的黄土孔隙比逐渐降低,土体在荷载作用下结构逐渐紧实。湿陷系数随着垂直荷载的增加,在300 kPa的压力下达到峰值。在低含水率的情况下,湿陷系数大于高含水率黄土,并随着冻融循环次数的增加逐渐减小。
(2) 随着冻融循环次数的增加,黄土的骨架颗粒发生偏转、变形,导致浸水前、后土粒形态大小发生变化,颗粒之间的连接方式由原有的点-点接触逐渐变为边-边接触,颗粒间的架空孔隙逐渐减小;颗粒轮廓从不规则形状逐渐边缘化、圆形化,颗粒从混乱、无序的分布向具有一定的方向性滑移,孔径的变化规律与定性分析结论一致。
(3) 冻融循环作为一种结构重组作用,通过减少主导湿陷变形的大孔隙空间,并增加结构的几何复杂性,从而驱动土体向更稳定的状态转化,宏观上表现为湿陷性的显著减弱。
(4) 黄土的架空孔隙结构崩塌以及颗粒间的排列方式变化是黄土湿陷变形的本质,为湿陷提供主要空间。随冻融循环次数增加,细碎的残片和片状粉粒含量显著增多,这些细小碎屑附着在粗骨架颗粒表面,增大了颗粒表面的粗糙度,导致了黄土的湿陷变形。
国家自然科学基金项目“伊犁谷地冰雪融水叠加降雨入渗黄土滑坡形成机理与预警判据”(42367021)
新疆维吾尔自治区天山英才培养计划“天山北坡矿山开采地质灾害-水文地质-生态环境耦合互馈机制及其工程地质意义”(2023TSYCCX0010)