高温冻土固化前后结构性变化对融化压缩特性影响

孙杲辰，张建明，党迎生，张 虎，丁 聪，陈 新

(1.冻土工程国家重点实验室(中国科学院西北生态环境资源研究院)，兰州 730000；2.中国科学院大学 工程科学学院，北京 100049； 3.中国市政工程西北设计研究院有限公司，兰州 730000；4.石河子大学 水利与建筑工程学院，新疆 石河子 832000； 5.上海勘测设计研究院有限公司，上海 200434)

土的结构性是指土中颗粒或集合体以及它们之间孔隙的性状和排列形式[1].1925年，太沙基对土体结构性进行了研究，他所倡导的土体微观结构概念与思想是土体微观结构研究的开端，但是由于当时观察设备落后，难以对土体复杂的结构性进行系统的定性描述；20世纪50年代至60年代中期，随着科学技术的发展，众多学者开始注重土体结构要素的完整性，特别是土体结构单元体的定向排列和分布特征；20世纪60年代后期至90年代中期，随着电子计算机技术的引入和新方法的出现，土体结构性由定性进入了定量阶段[2].中国土体结构性研究起步于1996年，许多知名专家和学者都对土体结构性进行了大量的研究和探索，如沈珠江院士[3]认为土的结构性是影响土体力学特性最为重要的要素，是“20世纪土力学的核心问题”，谢定义等[4]发展了土体结构性研究的土力学方法，开辟了土体结构性与定量化参数研究的新途径.然而，上述土体结构性研究方法仅限于常规土中，对于冻土特别是压缩性大、未冻水含量高、抗剪强度低的高温冻土结构性研究尚属空白.

众所周知，高温冻土又称近相变区冻土，国外许多学者也将其表述为“frozen soil at temperature close to 0 ℃”，通常用是指温度区间变化范围在-1.5～0 ℃，相对温度较高的冻土.高温冻土由土颗粒、未冻水、土中气体以及固态冰4种组分构成，由于特定的形成条件其结构性特征有着自身的特点，这种结构性对高温冻土的工程性质有重要的影响[5-7].多年以来，许多专家都以主动降温[8-10]，被动保护[11-13]的方式来保护冻土自身结构不因外界环境的变化而遭到破坏，但是在全球气候变暖的背景下，这些措施必然存在一定的局限性.因此，本研究从高温冻土固化的角度出发，对比不同固化剂加入高温冻土后，固化前后土体物理性质及结构性的变化规律，并以孔隙定向分维数Df与融化压缩系数av为桥梁，建立了冻土固化后微观孔隙结构特征与宏观力学性质之间的关系，为高温冻土结构性研究提供新的途径与方法.

1 试 验

1.1 试验材料及固化土的配比设计

本研究以采自青藏高原的粉质黏土为研究对象，其中土体基本物理指标及颗粒分布如表1，2所示.土体固化剂主要以无机类固化剂为主.近年来，有部分学者[14-16]通过普通硅酸盐水泥(portland cement，以下简称GC)对冻土进行了固化，但是受低温环境对其反应温度的限制，往往需要加入众多外加剂来促使其发生水化反应，单独使用时反应速率低且强度增长缓慢[17-19].然而，众多研究发现，高性能硫铝酸盐水泥(sulphate aluminum cement，以下简称SC)具有反应温度低且低温水化反应速率快等特点[20-21].因此，本研究通过在高温冻土中加入普通硅酸盐水泥(GC)、高性能硫铝酸盐水泥(SC)以及(普通硅酸盐水泥+硫铝酸盐水泥)(以下简称SC+GC)3种固化剂对比其固化效果.其中普通硅酸盐水泥采用甘肃祁连山水泥有限公司生产的P·O42.5型号普通硅酸盐水泥，硫铝酸盐水泥采用唐山北极熊水泥有限公司生产的C52.5硫铝酸盐水泥，两种水泥颗粒分布如表3所示.根据中华人民共和国《水泥土配合比设计规程》(JGJ/T 233—2011)水泥土掺入基准可根据使用目的及当地经验，按工程要求的水泥土性能指标确定，并宜取3%～25%，根据前人经验[22]以及高温冻土本身的性质，选取GC、SC以及GC+SC固化剂掺入质量分数为干土重的15%，其中3类固化土样的固化剂具体掺量如表4所示.

表1 青藏高原粉质黏土基本物理指标

表2 土体的颗粒分布

表3 水泥颗粒分布

表4 固化剂具体掺量

注:固化剂掺量均为干土重的质量分数/%.

1.2 试样制备

由表1可知，土体液限WP=18.6%，塑限WL=36.7%，为了确保土样的均匀性，采用静压法制备土样，试样的密度为1.46 g/cm3，同时设置土样的初始含水率为30%，试样制备步骤如下：

1)将试验粉质黏土充分在阳光下暴晒风干，之后碾碎并过2 mm筛，将100 g风干粉质黏土放入环刀中，置于烘箱中烘干24 h测定其土样天然含水率为3%，并将过筛后土样放置于塑料密封袋中密封，置于-5 ℃低温试验室中预冷24 h以待使用.

2)为了模拟青藏高原的真实环境，本次制样过程均在-5 ℃低温实验室内进行，将固化剂按照表4中掺量与预冷后土样混合搅拌均匀，并与事先准备好的冰粉与蒸馏水(冰粉过2 mm筛，冰水比例为3∶1)再次混合搅拌均匀后装入涂有凡士林的铁质模具中(φ61.8 mm×h40.0 mm).用油压千斤顶压实，脱膜后将其用保鲜薄膜包裹并置于-1 ℃的恒温箱中养护7 d后进行相关试验.制样过程详见图1.

图1 固化土样制样过程

1.3 试验方法

烘干法测定固化土样总含水量；利用土壤水分测定仪(time domain reflectometry，简称TDR)观测不同固化土样固化前后未冻水含量；通过扫描电镜试验(scanning electron microscope，简称SEM)观测固化前后高温冻土的微观孔隙结构；利用TD3500-X射线衍射分析仪进行固化高温冻土物相分析；运用GZQ-1型全自动气压固结仪，测试固化前后土体不同荷载与不同温度下所对应的压缩系数.

1.3.1 高温冻土固化前后总含水量测定

为了研究3种固化剂GC、 SC以及 SC+GC加入高温冻土后，高温冻土内部总含水量的变化规律，将养护7 d后的3种固化土样静置于105 ℃ DGG-9140B型烘箱中，烘干24 h测定固化前后土样总含水量变化.

1.3.2 固化前后土样未冻水含量确定

将养护7 d固化土样从-1 ℃恒温箱中拿出，并迅速将土壤水分测定仪(TDR)插入用塑料瓶及保鲜薄膜装好的土样中心(图2)，置于预先预冷过的-1 ℃智能低温恒温槽中，测定不同固化土样未冻水含量变化.

图2 固化前后土体未冻水含量测定

1.3.3 SEM试验土样制备

对于制备固化前后高温冻土微结构试样而言，难点在于如何置换与去除试样中的孔隙水与冰，保证试样不变形从而获取平整新鲜的扫描电镜观测面.已有的阴干法、风干法以及烘干法在试样干燥过程会产生明显的体积收缩，从而导致观测到的试样微结构与其真实结构之间存在较大误差[23-24].因此，从冷冻干燥技术入手，利用上海豫明仪器制造厂生产的“FA-1A-50冷冻干燥机”(图3)对试样进行干燥后再将土样转移至扫描电镜下观测其微结构.

图3 冷冻干燥机(FA-1A-50)

1.3.4 融化压缩试验

融化压缩试验在GZQ-1型全自动气压固结仪上进行(图4)，将在-1 ℃下养护7 d的No cement、GC、SC以及(SC+GC)4种土样分别置于4个压缩盒中，设定恒温箱的温度为-1 ℃静置24 h，待其完成后，在同一温度下陆续加载0.1，0.2和0.3 MPa 3级荷载.加载完成后，调整恒温箱温度至-0.5和0 ℃(每次调节温度时间应当控制试样变形少于0.025 mm/h)，重复以上加载过程.

图4 高温冻土固化前后融化压缩试验

Fig.4 Melting and compression test for high warm and frozen soil before and after stabilization

2 试验结果分析

2.1 高温冻土固化前后含水量变化规律

2.1.1 总含水量变化

高温冻土固化前后土样总含水量变化如图5所示.可以看出，No cement土样含水量为30.06%，随着固化剂的加入逐渐与高温冻土中的自由水发生水化反应， GC土样养护 7 d后的含水量由30.06%降至26.14%，总含水量降低了3.92%；与GC土样相比，由于硫铝酸盐水泥优越的低温反应性能，SC土样总含水量由30.06%降低至24.27%，总含水量降低了5.79%；SC+GC固化土样总含水量由30.06%降至22.33%，总含水量降低了7.73%，普通硅酸盐水泥与硫铝酸盐水泥相混合后，两者水化反应的过程中达到了动力上的协调，硫铝酸盐水泥进一步刺激了硅酸盐水泥的反应，使得水化反应更加剧烈，土体总含水量大幅度降低.

图5 高温冻土固化前后总含水量柱状分布

2.1.2 未冻水含量变化

固化剂的加入虽然降低了高温冻土中的总含水量，但是由于其与土体中的自由水发生水化反应放出的热量融化了土颗粒周围的结晶冰，增加了土体内部未冻水含量.高温冻土固化前后未冻水含量分布如图6所示，与总含水量相比，固化土体内部未冻水含量变化与总含水量变化趋势正好相反，固化效果越好，土体内部未冻水含量越多，但是随着温度的降低土体内部未冻水含量基本稳定在10.0%左右.

图6 高温冻土固化前后未冻水含量变化

2.2 固化后土体结构单元体粒度定量分析

图7为高温冻土固化后结构单元体粒度定量分析结果.可以看出，固化后结构单元体的直径主要分布在2～5 μm与5～10 μm内，其他区间内分布较小.随着固化剂加入后固化效果不同，结构单元体粒度分布不同，其中在5～10 μm和10～20 μm反映了土体固化后水化反应胶结作用土体颗粒团聚的主要结果，SC+GC土样等效直径分布几率最大为52.85%，其次为SC土样分布几率为50.82%,而GC土样分布最小为46.68%.充分说明SC+GC固化剂加入土体后，土颗粒团聚效果最好，大颗粒增多.

图7 高温冻土固化后固化土样结构单元体粒度分布

Fig.7 Element size distribution histogram of stabilized soil structure

2.3 高温冻土固化前后孔隙面积变化规律

高温冻土固化前后孔隙面积在各孔径范围内面积分布的改变是其微观孔隙结构改变最直接的表象.由表5与图8可以看出，加入固化剂后高温冻土各孔径范围内孔隙面积均呈现减小规律，但是固化剂的种类不同，受低温条件与反应速率的限制，其孔隙面积改变规律有所不同，整体表现为：SC+GC土样，孔隙面积减小最多；SC固化土样孔隙面积次之；GC固化土样孔隙面积减少最小.

表5 固化前后各孔径对应孔隙面积分布 m2·g-1

图8 固化前后孔隙面积分布

2.4 固化前后土体孔隙丰度评价

固化前后土体孔隙的丰度值C为孔隙短轴B与长轴A的比值，C=B/A,该值代表了分析区域内孔隙的圆度，即土体孔隙展示的二维平面几何形状，反映了土体的均匀程度.丰度值C的取值范围在(0，1)，C值越小，表明孔隙平面形状趋向于长条形，越不均匀；C值越大，则说明孔隙平面形状趋向于圆形，土体越均匀[25].本次研究将丰度值划分为8个统计区间，各区间内丰度值分布如图9所示.结果表明，No cement固化土样的丰度值主要集中在0.2～0.3，土颗粒孔隙主要以狭长状为主；GC固化土样丰度值集中在0.3～0.4，说明普通硅酸盐水泥加入高温冻土后，随着水化产物对土体孔隙的填充，土体孔隙由狭长状向半椭圆状转化；而SC与SC+GC固化土样的丰度值主要集中在0.8～0.9，说明硫铝酸盐水泥加入后，土体孔隙由椭圆状向圆形转换，土体整体趋于均匀化.

2.5 固化高温冻土XRD试验结果分析

为了观察固化后，土样内部水化产物生成总量，将养护7 d后的固化高温冻土样品冷冻干燥后，经器具研磨过0.08 mm方孔筛，然后通过电子天平称取5 g试样，利用TD3500-X射线衍射分析仪进行固化高温冻土物相分析.扫描分析采用连续扫描测量法，扫描起始角为10°，终止角为70°，并利用JADE6.0软件对衍射结果进行统计.从图10分析结果可以看出，GC固化土样中钙矾石(AFt)与水化硅酸钙凝胶(C-S-H)仅有微小峰值，而SC和SC+GC固化土样的钙矾石(AFt)峰值较大，特别是硫铝酸盐水泥与普通硅酸盐水泥混合加入土体后，不仅提高了钙矾石(AFt)的峰值，而且促进了水化硅酸钙凝胶(C-S-H)的生成，充分证明硫铝酸盐水泥具有低温水化反应速率快的特点.

图9 固化前后孔隙丰度值

图10 固化高温冻土7 d X射线衍射图

Fig.10 X-ray diffraction of solidified high warm and frozen soil samples at 7 d

2.6 固化前后微观结构变化固化机理分析

图11为-1 ℃下养护7 d的5 000×倍扫描电镜照片，可以看出，SC+GC固化土样(图11(d))经过7 d的养护生成了大量的钙矾石(AFt)与水化硅酸钙(CSH)凝胶，AFt与CSH相互交错黏结在一起，不仅填充了孔隙空间，而且相互搭接形成了牢固的空间结构；SC固化土样(图11(c))与SC+GC固化土样相比，钙矾石(AFt)与水化硅酸钙(CSH)凝胶生成较少；GC固化土样(图11(b))基本无钙矾石(AFt)的生成.充分证明，与硫铝酸盐水泥相比，普通硅酸盐水泥负温下反应速率缓慢，单独加入硫铝酸盐水泥时，虽能生成一定的钙矾石(AFt)，但是水化产物量较少.当普通硅酸盐水泥与硫铝酸盐水泥混和时，不仅加快了钙矾石(AFt)的生成速率，减少了钙矾石(AFt)对凝胶类水化产物的破坏作用，而且修复了凝胶类水化物被破坏胶结结构，协调了两者反应的动力学关系，改变了土体内部的孔隙结构，使得土体内部孔隙更加致密.

图11 不同固化土扫描电镜图片(5 000×)

2.7 固化前后孔隙定向性分布特征

为了研究高温冻土固化前后土体孔隙分布状况，将通过CURVEEXTRACT预处理过的图片导入Image-ProPlus6.0(IPP)软件中统计固化土样的孔隙角度以及每个角度区间内的概率分布，从而计算土体孔隙的定向性.土体孔隙的定向性定义为孔隙长轴与水平x轴的夹角θ.因为孔隙方向在0°～180°和180°～360°具有对称性， 以0°～180°的玫瑰风向图，表示土体孔隙长轴方向的角度频率分布，其中共划分了14个孔隙数目统计区间.

孔隙分布概率计算如式(1)所示：

(1)

式中：ni，n分别为夹角θ落入(θi-1，θi)的土体孔隙数量和孔隙总数[26].

通过冻土物理学可知，黏土颗粒在不同应力作用下会促使其呈一定方向的排列.从图12固化前后高温冻土玫瑰风向图可以看出，对于固化前后高温冻土，孔隙定向性分布具有一定的相似性，孔隙长轴与水平x轴之间的夹角主要分布在90°以及180°左右2个区间内，表明土体在该区间内均表现出较大的定向性.其中不同固化土体90°以及180°区间内分布概率见表6.可以看出，SC+GC固化土样与SC土样在90°左右分布概率为23.25%和25.25%，而90°区间内长轴方向与自重应力方向相同，其在该区间内孔隙概率分布越小，从微观尺度上就能更好地反映其颗粒接触越充分，从而表现出压缩系数越小，变形越小的宏观工程性质.

图12 不同固化土孔隙长轴方向玫瑰风向图

表6 不同固化土体90°及180°孔隙分布概率

Tab.6 90°and 180°pore area probability before and after curing

2.8 固化前后孔隙定向性分维分析

众所周知，土体是具有一定自相似性的分形结构体.根据分形理论来定量描述高温冻土固化前后的微观结构特征，不仅表现方法直观化，而且可以定量描述固化前后土体力学性质[27].本研究使用Mandelbort提出的容量分维方法，计算固化前后土体孔隙的平面定向分维，计算方法如下：

N(r)=cr-Df.

(2)

式中：r代表格子的边长，N(r)表示所占的格子数.绘制高温冻土固化前后lnN(r)-lnr散点图，如图13所示，lnN(r)-lnr曲线均具有较好的线性关系.定向分维(Df)反映土体孔隙的定向程度，Df值越小，土的定向性越好，Df值越大，土的定向程度越差，混乱度越大.

Fig.13 Pore direction fractal dimension relationship curves of solidified high warm and frozen soil before and after curing

通过表7可以看出，固化土样的分维值在1.009 9～1.103 7变化，固化后土样的定向分维呈现减小的趋势，且SC+GC固化土样的定向分维值最小，SC固化土样定向分维值次之，GC固化土样定向分维值较大.充分说明该区间内定向分维值Df越小，孔隙分布概率F(θ)越小，这与前面所描述的孔隙定向频率相吻合.

3 讨 论

3.1 高温冻土固化前后融化压缩试验研究

土体压缩系数av是反映高温冻土宏观力学性质的一个重要指标.根据土工试验规(GB/T50123—1999)，压缩系数是指外部荷载0.1～0.2 MPa作用下变形量与荷载变化的比值，即

av=Δh/(H×ΔP).

(3)

式中：av为融化压缩系数(MPa-1)，Δh为土样在0.1～0.2 MPa时总变形量(mm)，H为试样初始高度(mm)，ΔP为荷载差值.

表7 孔隙定向性分维统计

图14为固化前后土体压缩系数变化.可以看出，No cement土样压缩系数受温度变化影响，且其变化范围较大.负温下压缩系数变化范围较小，主要因为土体中存在冰晶，对土体骨架起到了支撑作用，当温度升至0 ℃时，由于土中冰晶的融化，土体压缩系数显著增大.而对于3种固化土来说，压缩系数受温度变化影响较小，但是由于固化效果不同，不同固化土压缩系数相差较大，其中SC+GC固化土样压缩系数最小，即使温度升至0 ℃变形量也较小，充分证明硫铝酸盐水泥+普通硅酸盐水泥固化剂的加入，改变了土体内部的孔隙结构，水化产物替代了土体内部的冰晶成为冻土的主要支撑结构.

图14 固化前后土体压缩系数变化

3.2 高温冻土融化压缩系数与孔隙定向分维关系

通过以上分析可知，固化剂的加入彻底改变了土体的孔隙结构.固化剂与冻土中自由水反应产生的水化产物取代了冰晶，成为土体结构的主要支撑.为了更加真实地反应固化前后孔隙定向分维Df与压缩系数av之间的关系，选取0 ℃时土体压缩系数av与固化后土体孔隙定向分维数Df进行拟合，结果见图15.可以看出，孔隙定向分维数随着压缩系数的增大而增大，两者之间存在良好的线性关系，孔隙定向分维数越大，土体力学性质越好.

图15 孔隙定向分维Df与压缩系数av之间关系

Fig.15 Correlation between pore directional fractal dimensionDfand compressibility coefficientavof samples

3.3 高温冻土固化反应微观模型

黄新、宁建国等[28]通过普通硅酸盐水泥对常规粉质黏土进行固化，提出了水泥固化土微观机制模型，指出土体固化实质是水泥发生水化反应后，产生的水化产物对土颗粒进行包裹、胶结并填充的过程.根据上述高温冻土固化前后微观孔隙结构变化规律以及融化压缩性质结果，参照刘松玉等[29]碳化固化粉土模型，提出高性能水泥固化高温冻土微观结构反应模型.掺入硫铝酸盐水泥+普通硅酸盐水泥混合固化剂时，硫铝酸盐水泥以其具有反应温度低、与水化反应速率快的特点，快速与高温冻土中的大量自由水发生反应形成水化硅酸钙凝胶(CSH)包裹在(土-冰)团粒表面(见图16(a)).随着水化反应进程的加剧，水泥水化放出的热量逐渐开始融化团粒之间的冰，冰融化后产生的自由水进一步与剩余水泥反应生成CSH，将土颗粒胶结在一起(见图16(b)).同时，六棱柱状的钙矾石AFt迅速在CSH包裹土团粒表面生长，AFt晶体与纤维状的CSH凝胶纵横交替搭接成复杂的网状结构，通过填充与压实作用来减小土团粒之间的孔隙(见图16(c)).硫铝酸盐水泥与普通硅酸盐水泥混合后，加快了钙矾石的生成速率，增强了凝胶类水化产物对破坏胶结结构的修复作用，形成了水化过程中的动力学关系.

图16 高温冻土固化反应微观模型

4 结 论

1)固化剂加入高温冻土后，与土体内部的自由水发生反应，降低土体总含水量的同时，生成的水化产物将不同粒径的颗粒胶结在一起产生土体团聚作用，使得在5～10 μm和10～20 μm土颗粒等效直径分布几率增大，其中SC+GC固化剂固化效果最好，SC次之，GC最差.

2)高温冻土固化后冻土孔隙面积主要分布于r<2 500 nm区间，孔隙长轴与水平x轴之间的夹角主要集中于90°及180°左右.固化后土体孔隙面积大幅减小，在90°，180°左右表现出较大的定向性.其中SC+GC固化效果明显，孔隙面积减小最多，且90°，180°左右定向性表现突出.充分证明硫铝酸盐水泥与普通硅酸盐水泥混合后，加快了钙矾石的生成速率，增强了凝胶类水化产物对破坏胶结结构的修复作用.90°，180°左右定向性孔隙概率分布越小，从微观尺度上能更好地反映其颗粒接触越充分，压缩系数较小，变形较小的宏观工程性质.固化剂与高温冻土中未冻土发生水化反应替代冻土中冰晶成为冻土中主要的支撑结构.

3)固化后土体孔隙定向分维Df与压缩系数系数av之间存在较好的线性关系，且孔隙定向分维Df越大，压缩系数av越大，土体分布越不均匀，Df越小压缩系数av越小，土体越趋于均匀化.