刘 力, 万 旭 升*, 王 智 猛, 邱 恩 喜, 龚 富 茂, 钟 昌 茂, 吴 松 波
( 1.西南石油大学 土木工程与测绘学院, 四川 成都 610500;2.中铁二院工程集团有限责任公司, 四川 成都 610031 )
特殊土在工程建设中比较常见,其特有的性质和危害给工程建设带来了不少困难.盐渍土作为一种特殊土,是指易溶盐含量大于0.3%,并具有溶陷、盐胀和腐蚀等工程特性的岩土,它在我国分布范围很广,尤其是西北干旱地区[1].盐渍土的盐胀融沉引发的基础工程问题严重制约了该地区经济发展.土的渗透性是土体评价的一项重要内容,渗透系数则是衡量土体渗透性能的一项重要指标.研究盐渍土的渗透系数对防治和改良其病害有着重要的意义.
土体渗透系数的影响因素有很多,王双等[2]采用神经网络研究了碎石土级配特征对渗透系数的影响,而Shafiee[3]、朱崇辉等[4]研究了混合黏土的级配对渗透系数的影响,发现颗粒的圆度越差,渗透系数越大.邵龙潭等[5]提出了一种渗透系数的测量方法并给出了经验预测公式.Zhou等[6]、邓友生等[7]、孙大松等[8]通过室内试验,得出了渗透系数与干容重、土颗粒大小和所受固结压力的关系,并推导了非饱和土渗透系数的理论表达式.然而,渗透系数的测定非常依赖于仪器设备,由不同仪器设备测量而得到的结果区别较大,目前还没有建立盐渍土渗透系数理论模型.在土的渗透改良研究方面,张登武等[9]发现硫酸盐水泥、石灰和粉煤灰不同比例的配合对盐渍土的改良效果比单独使用其中任何一种好,所以对其进行了冻融试验和渗透试验.用混合料改良后,盐渍土的渗透性很小.米海珍等[10]研究了石灰土、二灰土和水泥土3种改良土的渗透性能,得出了防渗土料掺入料的最佳配比.随着纳米材料的兴起,很多学者都尝试利用纳米材料进行改性研究,Ren等[11]、Rong等[12]、陈学军等[13]通过一系列试验研究,发现新型纳米材料可以改变土体的物理性质,纳米材料(纳米二氧化硅、纳米碳酸钙)可以提高黏土的单轴抗压强度和液塑限,减小其平均孔径,提高水泥抗压抗折强度并影响红黏土黏聚力、内摩擦角和抗剪强度,从而改变材料的强度特性.虽然微纳米材料在很多工程领域都有了应用成功的经验,但关于其对盐渍土渗透性能影响的研究并不多见.
本文针对盐渍土渗透模型不健全及渗透特性改良方面存在的问题,通过建立颗粒级配曲线与孔径分布的关系,借助土水特征曲线建立渗透系数计算模型.并通过室内渗透试验,研究添加微纳米材料(纳米二氧化硅、纳米碳酸钙和相变材料)后硫酸钠盐渍土渗透性能.进一步探究盐渍土中溶液迁移机理和渗透特性改良机理,以期为盐渍土地区工程实践提供一定的参考依据.
从微观孔隙角度出发,将土水特征曲线作为反映土体孔隙通道的间接指标,陶高梁等[14]建立了饱和土的渗透模型:
(1)
式中:ks为饱和土的渗透系数;ρw、g为水的密度和重力加速度(各取1 000 kg/m3、9.8 m/s2);σ为表面张力,25 ℃时取0.072 N/m;α为接触角(饱和土可取α=0°);τ为孔隙通道实际长度与土样长度的比值,即迂曲度,可视为常数2[15];μ为水的黏度,25 ℃时取0.893 7×10-3Pa·s;ψm为基质吸力(kPa),其表达式为
ψm=2σcosα/r
(2)
其中r为孔隙半径;θmax、θmin分别为最大(饱和)和最小体积含水率.
考虑到盐渍土中含有盐分,土体中不仅存在基质势,还会存在溶质势.当含盐量较高时,溶质势不可忽视,所以将式(2)中ψm项修正为土水势ψ,其表达式为
ψ=ψm+ψs
(3)
ψs为溶质势,在常温条件下,假定盐分充分溶解到孔隙水中,则溶质势可表示为[16]
ψs=ρwRTlnaw/Mw
(4)
其中R为气体常数,取8.314 J/(mol·K);T为温度(K);Mw为水的摩尔质量,取18 g/mol;aw为盐溶液中的水分活度[17-21].故饱和土的渗透模型可修正为
(5)
土体中孔隙分布特征是决定其渗透性的关键.利用分形理论研究多孔介质孔隙分布是一种好方法,但是对于土体而言,更简洁有效的方法是从粒度分布的角度去研究孔隙分布,采用孔径、含水量与粒度分布的关系[22]:
(6)
式中:ri为第i个粒度组所对应的孔隙半径;e为土体孔隙比;γ为尺度参数,取1.4[22];Ri为第i个粒度组的平均粒径;ni为第i个粒度组的颗粒数量,可用下式计算:
(7)
其中wi为颗粒曲线中第i个粒度组含量,ds为土粒相对密实度.对颗粒级配曲线划分微段,可建立wi与Ri的关系,具体如图1所示.
图1 土颗粒粒径分布曲线划分
(8)
式中:θi为小于第i个粒度组尺寸的累积体积含水率;Φ为土体孔隙率;Sw为饱和含水量测量值与孔隙率理论值的比值,取1.
对于特定孔隙半径ri而言,其土水势可由下式计算:
(9)
为了建立含水量与土水势的关系,首先通过颗粒级配曲线建立不同粒径与其对应含水量之间的关系,然后基于孔径与粒径的半经验公式建立不同孔径与对应含水量的关系,进而得到土水特征曲线.具体步骤如图2所示.
图2 变量关系图
试验用土取自甘肃兰州地区,取土深度为0.5~1.0 m,土性为粉质黏土,土样基本物理性质见表1,土样粒度分布如图3所示.掺加的微纳米材料为纳米二氧化硅、纳米碳酸钙和相变材料,3种材料外形均为白色粉末,粒径分别为(30±10) nm、40~80 nm和1~50 μm.纳米二氧化硅比表面积150~300 m2/g,含量≥99.5%,堆积密度0.057 g/cm3.纳米碳酸钙主含量>97%,pH 8.0~11.0,沉降体积>2.4 mL/g.相变材料为保温材料,主要成分为改性芯材、甲苯二异氰酸酯、聚乙二醇,相变热165 J/g.本次试验采用人工配制硫酸钠盐渍土,以消除试验差异并控制土样含盐量.试验用土经过脱盐、烘干、碾碎和过筛,最后进行密封保存.
表1 土样基本物理性质
图3 土样颗粒级配曲线
本次试验采用变水头渗透试验方法,试验装置为TST-55型渗透仪,如图4所示.渗透容器由环刀、透水石、套环、上盖和下盖组成,环刀内径61.8 mm、高40 mm,透水石的渗透系数大于0.001 cm/s;变水头装置由渗透容器、变水头管、供水瓶、进水管等组成,变水头管内径6.6 mm、长2.0 m,满足《土工试验方法标准》[23]的要求.
试验土样设置含盐量为0、1%、2%、3%和4%,由此探寻含盐量对土样渗透性的影响规律.
图4 TST-55型渗透仪
再取高含盐量(3%)土样,分别掺加不同含量(1.0%、1.5%、2.0%、2.5%、3.0%)的纳米二氧化硅、纳米碳酸钙和相变材料.文中涉及含量的数据均为质量分数.
制样时,控制土样含水量为18.5%,首先按计算数据配制硫酸钠盐溶液,需要掺加掺料时,先将称量好的土和掺料充分混合后,再均匀地喷入盐溶液,将盐溶液和土样充分拌和,装入密封袋并放置于25 ℃恒温箱中稳定24 h,进一步使盐溶液和土样均匀混合.装样时取试样干密度为1.6 g/cm3.记录试验数据时,需要同时记录时间、水柱高度和出水温度,当水柱下降一定高度(大于10 cm)后再次记录.最后改变水柱高度,对同一试样重复5次试验,以便计算得出渗透系数平均值,减小试验误差.其余具体试验方法和步骤均按照《土工试验方法标准》进行.试验装置简图如图5所示.
图5 试验装置简图
试验过程中盐溶液不断填充孔隙进而发生渗流,故在计算过程中可认为土样是饱和的.不同含盐量土样对应盐溶液浓度见表2.
由于土样孔隙在渗流过程中不断填充,可近似认为有水流流出时为饱和状态,可根据式(1)~(4)计算盐渍土渗透系数.
表2 不同含盐量土样对应盐溶液浓度
处理数据时按变水头渗透系数计算公式进行计算[18].
(10)
式中:a、A分别为变水头管和土样的断面面积(cm2);2.3为ln和lg的变换因数;L为渗径,即试样高度(cm);t1、t2分别为测读水头的起始和终止时间(s);H1、H2分别为起始和终止水头.
按照式(10)计算得出渗透系数,进行温度校正后对3个在允许误差值范围内的数据取平均值,作为一个试样在标准水温(20 ℃)下的渗透系数.
R2=0.999 6
(11)
为了获取高精度,需要土颗粒级配曲线拟合函数,对土颗粒粒径取对数之后再按间隔0.065分组,即
lgRi-1=lgRi-0.065
(12)
对土颗粒粒径取对数后再进行均匀分组,所得土颗粒含量分布与试验值非常接近,如图6所示.
图6 土颗粒含量拟合曲线
修正后土水势、溶质势与水分活度的关系如图7所示.由图7(a)可知,溶质势随着水分活度增大而增大;由图7(b)可知,土水势随着含水量增加而降低,当土体含水量较低时,溶质势对土水势的影响较小,当土体趋于饱和时,溶质势对土水势的影响较明显.
(a) 溶质势与水分活度的关系
(b) 不同水分活度土水势与含水量的关系
试验测试并计算得出最终渗透系数,取其平均值的数据点绘制成图,从图8中可以看出,在同一干密度下,随着含盐量的增加,盐渍土的渗透系数呈减小趋势,与邓友生等[7]所得结论保持一致.
试验中兰州粉质黏土渗透系数首先随着含盐量的增加大幅度降低,当含盐量大于2%后,其渗透系数随着含盐量的变化很小.通过计算值与试验值对比可知,本文修正的计算模型可较好地计算硫酸钠盐渍土渗透系数.土样渗透系数随含盐量增加而减小,是因为硫酸钠盐的存在降低了盐渍土土水势,而当含盐量较高时,盐分达到饱和状态致使其水分活度基本不变.但是高含盐量硫酸钠盐渍土中部分盐晶体的析出会造成土中孔隙减小,致使渗透系数微弱减小从而降低了土体渗透特性.
图8 含盐量对盐渍土渗透系数的影响
微纳米材料对盐渍土渗透系数影响如图9所示.纳米二氧化硅和相变材料可以降低盐渍土的渗透系数,并随着掺量增加,渗透系数降低更多,掺加相变材料的盐渍土其渗透系数变化幅度较大.低掺量时,掺加相变材料的盐渍土与无掺加土样的渗透系数相差范围很小,而掺加纳米二氧化硅后的土样渗透系数则比无掺加时明显减小.当掺量大于2.5%时,掺加两种材料后的渗透系数相差不大,均比无掺加时有明显减小.掺料为纳米碳酸钙时,随掺量增加,土样渗透系数总体呈上升趋势,掺量为1.0%时,其渗透系数比无掺加时有微量减小,其余掺量对土样渗透系数均为增加作用.
图9 各掺料对3%含盐量盐渍土渗透系数的影响
掺加纳米材料后,相当于降低了土样的初始含水量,在掺料适量添加后,纳米材料形成团聚体,起到较好的填充作用,如图10所示.由于纳米二氧化硅堆积密度较小,随着其掺量增加,土样密度减小,从而渗透系数变小减缓.纳米碳酸钙作为亲水材料在用量增大后,表面活化剂的浸润和增塑作用更加显著,加大了纳米粒子的流动分散[24],而且纳米碳酸钙纳米级的较细粒径,表面能高,具有较强的吸附作用,变相减小了试样含盐量,使得试样的渗透系数增大;随着渗流进行,在吸附作用影响下含盐量增大,使得试样渗透系数减小.当相变材料的掺量较低时,相变材料会释放热量,在一定程度上影响土样内部温度从而减小水的黏滞系数,使渗透系数增大.作为疏水材料的相变材料,在大量添加后起到了较好的填充作用,有效降低了渗透系数.
(a) 未添加纳米材料
(b) 添加纳米材料
但是,当土体含水量较低时,溶质势对土水势的影响较小,故在非饱和土渗透系数的计算中可忽略溶液的影响.本研究中未考虑环境温度变化对硫酸钠盐渍土渗透特性的影响,对硫酸钠盐渍土而言,在温度降低时会生成十水硫酸钠晶体引发土体盐胀,故盐胀变形引发的孔隙变化对土体渗透特性的影响会在以后的工作中继续探究.
(1)随着土体趋于饱和,溶质势对土水势的影响越来越大,故针对饱和盐渍土渗透系数进行土水势修正极为必要;推导得到的饱和盐渍土渗透系数计算模型的计算值与试验值接近,该模型符合饱和硫酸钠盐渍土渗透系数的变化规律.
(2)当含盐量在4%以内时,硫酸钠盐渍土中土水势随着含盐量的增加而逐渐减小.同样,渗透系数随着含盐量的增加逐渐减小,当含盐量大于2%后,渗透系数变化很小.
(3)对于高含盐量硫酸钠盐渍土,掺加纳米二氧化硅和相变材料可以降低渗透系数,随掺量增加渗透系数减小更多.当纳米碳酸钙掺量大于1%后,渗透系数随着其含量的增加而大幅度增加.硫酸钠盐渍土渗透系数随着相变材料的增加先增大后减小,当掺量大于2.5%后变化很小.