基于动态多步流动法的非饱和土持水特性试验研究

2021-11-16冯巧云毕庆涛马佳兵梁世春

水利与建筑工程学报 2021年5期

冯巧云，毕庆涛，马佳兵，梁世春，徐迪

(1.华北水利水电大学，河南郑州 450046；2.山西省交通规划勘察设计院有限公司, 山西太原 030012)

非饱和土是一种由固相、液相和气相组成的三相体。与饱和土相比，气相的存在使非饱和土的性质更加复杂，这给非饱和土工程特性的研究带来了许多困难。随着我国现代化建设的推进，工程中遇到的非饱和土问题日益突出，由于目前缺乏合理的非饱和土力学理论来指导工程设计和施工，引发了很多工程事故[1]。其中，基质吸力是影响非饱和土工程性质的重要参数。通常采用土-水特征曲线(简写为SWCC曲线)描述基质吸力与饱和度之间的关系[2]。SWCC曲线作为非饱土和力学中最基础的本构关系，它对非饱和土工程特性，如渗透系数[3-5]，抗剪强度[6-8]、应力和应变[9-10]等的研究具有重要意义。

国内外很多学者开展了非饱和土水力特性研究，并对SWCC曲线进行了深入的探讨。Tarantino[11]认为，通过改变孔隙比大小使土的密实程度发生变化，进而导致SWCC曲线的形态随之改变；李同录等[12]用压汞试验、滤纸法测SWCC曲线，并用扫描电镜获得其微观结构图像，进行了击实黄土孔隙结构对土水特征曲线的影响分析；郑方等[13]利用GCTS土-水特征曲线仪对三种不同粒度的土进行增减湿循环试验，分析粒度对非饱和土土-水特性的影响；Wang等[14]通过分析压力板试验原理和土样的微观结构，发现尺寸较小的圆形试样可以缩短SWCC曲线测试的时间；李幻等[15]在改进的压力板仪测试系统上对粉土试样开展动态多步流动实验，最终得到平衡状态下SWCC脱湿曲线，与常规平衡态测试方法相比可以大大节省时间。

目前关于非饱和土SWCC曲线的研究大多数都是基于传统测试方法，即在平衡状态下测试，一条SWCC曲线的测试时间为数周甚至数月，耗时较长。动态多步流动法[16]的提出为缩短试验时间做出了贡献。该方法是在试样含水率与基质吸力未达到平衡时，就施加下一级的基质吸力；通过理论分析，得到平衡状态下试样的SWCC曲线。

本文以山西太原非饱和黄土为对象，在改进的压力板仪上利用动态多步流动法测定原状土和重塑土的土-水特征曲线，并对比分析两者差异，为非饱和土持水特性的研究提供参考。

1 仪器改进和试验方案

1.1 试验仪器

试验利用GEO-Experts压力板仪测定非饱和黄土原状样与重塑样脱湿过程的SWCC曲线。此压力板仪存在不能实时监测溢出水的质量，测试时间较长，无法准确测量陶土板底部气泡体积等问题，导致动态多步流动法试验无法进行，因此需要改进。改进后的压力板仪，增加了数据自动采集系统，气泡体积测量系统和储水冲刷系统，如图1所示。

图1 改进的压力板仪测试系统示意图

改进后的压力板仪可以人为设定记录溢出水质量的时间间隔(本试验设置为每30 s采集一次)，避免了试验过程中将试样拿出时由于土体掉落产生的误差；在重力的作用下，打开储水冲刷系统的开关后可将陶土板下产生的气泡排除，气泡体积测量系统可以精确的测量气泡体积，减小试验数据的误差。

1.2 试验方案

试验土样取自山西太原某工程的非饱和原状黄土。根据《土工试验方法标准》[17](GB/T 50123—2019)，测得土样的基本物理性质指标如表1所示，颗粒级配曲线见图2。

表1 试验黄土基本物理性质指标

图2 试验黄土颗粒级配曲线

制备原状土与重塑土环刀试样(直径61.8 mm，高20 mm)。制备原状样时，把环刀放在水平土样层面，将环刀以外的土切掉，缓慢旋转环刀至土样完全充满环刀，过程中应尽量避免环刀内土体结构受到扰动；利用原状土制备重塑样时，将土样烘干(105℃)，碾碎，利用制样器与千斤顶制备与原状样相同含水率与干密度的重塑样。最后对试样进行真空抽气饱和后开始试验。

开始试验前，将底座中的陶土板和整个系统线路进行饱和。然后将土样放入压力室，转动调压阀施加一小吸力至10 kPa，使得试样从饱和状态变为非饱和状态，即不再有水溢出后，开始动态多步流动试验。各级吸力的加载时间人为确定，这里参考伊盼盼等[16]使用的加载方案，设定为20 kPa→30 kPa→40 kPa→60 kPa→90 kPa→130 kPa→170 kPa→210 kPa→290 kPa，加载所用时间见表2。利用数据采集系统记录溢出水的质量。试验过程中，在每级吸力加载结束后要对聚集在陶土板底部的气泡体积进行冲刷和测量，用以修正各级吸力作用过程中溢出水的质量。

表2 动态多步流动试验加载时间

2 试验数据分析

2.1 溢水量

改进后的压力板仪可以采集到饱和后土样溢出水累计质量的数据，但由于试验过程中陶土板底部气泡的存在导致采集系统监测的数据与实际溢水量相比偏大，需将排除的气泡体积换算成等体积水的质量，即多余溢水质量，然后对原始采集数据进行修正，可得到相对精确的试验数据。

对修正后的溢水量进行处理，获得了试样溢水量和基质吸力与时间的关系，如图3和图4所示。从图中可以看出，试验前期，原状样和重塑样的溢水量和溢水速度均较小，曲线较平缓，这是因为施加的吸力较小，水不易从试样中排出。试验中期，对原状样而言，当基质吸力达到40 kPa～60 kPa时，曲线趋势开始变陡，斜率较大，此时溢水速率较大；而重塑样的溢水速率开始变大时，对应的基质吸力则为60 kPa～90 kPa。试验后期，随着基质吸力的变大，原状样与重塑样的溢水速率开始逐渐变小，曲线趋于平缓，此时试样将要达到残余状态；且与重塑样相比，原状样的溢水量较大，达到残余状态所需时间较短。

图3 原状样溢水量和基质吸力与时间关系

图4 重塑样溢水量和基质吸力与时间关系

2.2 饱和度方程

由图3、图4可以看出，试验在各级吸力作用下未达到平衡状态，为此伊盼盼等[16]结合动态多步流动法原理，推导出非饱和土的饱和度演化方程：

(1)

方程(1)考虑到了各级基质吸力的增量对饱和度的影响，可将这种影响结果称为各级吸力下的“遗传效应”。

饱和度方程对于平衡态与非平衡态均适用，这为非平衡态下试验数据确定平衡态下的土水特征曲线奠定了基础。方程(1)中Ci与τi为未知量，其余量均为试验实测量。

可以利用Origin 软件建立饱和度随时间的演化方程模型，编译成功后，把试验数据导入软件中进行参数拟合，得到每级基质吸力对应的未知参数Ci和τi。

根据上述方法进行拟合，得到原状样与重塑样在各级基质吸力下的饱和度方程参数容水率Ci和特征时间τi，如表3所示。

表3 饱和度方程参数表

利用饱和度方程反算非平衡态试验结果，并与实测结果对比，如图5、图6 所示；可以看出实测结果和基于饱和度时间演化方程模型的拟合结果曲线几乎重合，各级基质吸力对应的曲线拟合相关系数均达到0.95，拟合度较高。

3 土-水特征曲线

土-水特征曲线是指平衡状态下饱和度与基质吸力的关系曲线。由上节饱和度方程得到每级基质吸力增量下的平衡态土水特征曲线斜率C，从而可以计算出第k级基质吸力下对应的平衡态饱和度，即

(2)

(3)

根据式(2)、式(3)，便可绘制土-水特征曲线。由试验结果，原状样与重塑样的土-水特征曲线如图7所示。

图5 原状样饱和度、基质吸力与时间的关系

图6 重塑样饱和度、基质吸力与时间的关系

图7 土-水特征曲线

由图7可以看出，原状土的土-水特征曲线位于重塑土的下方，即在相同基质吸力作用下，原状土的饱和度均小于重塑土，二者差异明显。饱和度的降低速率体现了土体内部水量减少的快慢，这与施加的吸力大小，土样孔径及内部含水率有关。原状土与重塑土的饱和度降低速率均随着基质吸力的逐级增加呈现先增长后减小，最终曲线趋于平缓的规律。

通过观察容水率C值的大小和土水特征曲线趋势可以看出原状土的进气值介于40 kPa～60 kPa，而重塑土的进气值在 60 kPa～90 kPa之间。由土的应力历史可知，原状土受上部土体长期作用后，内部颗粒按一定方向和次序排列，会形成具有方向性的孔隙结构，气体首先进入尺寸较大的通道，因此进气值较小；重塑土样为静压制作，孔隙分布无方向性，水流路径较少，导致气体较难进入土体，表现出较高的进气值。

试验后期，继续增加基质吸力，土-液接触面之间的毛细作用力增大，试样中孔隙中的水分会持续排出，但是排水速率逐渐减小。当基质吸力达到130 kPa时，原状土排水速率开始低于重塑土，二者的土-水特征曲线开始并行，直到毛细作用力与基质吸力达到平衡时试样饱和度不再随基质吸力的增大而减小，达到残余饱和度状态，将此时的饱和度转换为对应的体积含水率即为残余含水率。重塑土的残余含水率高于原状土，重塑土的持水性能优于原状土。