在中国东南部沿海地区有大量的吹填造地项目，以前吹填土方来源以砂为主，但是目前由于砂料不足，许多地区采用近海或河流底部的淤泥、淤泥质土为吹填材料，经水力吹填到场地内的这种超软土的物理力学性质很差，具有含水率高、黏粒含量高、压缩性高、承载力低等显著特点^[1-3].

为了加快这类超软土地基的固结，国内常使用无砂真空预压法^[4-5]，该方法具有施工方便、成本低廉、工期短等特点.真空预压法需要在超软土中设置竖向排水体，目前由芯板和包裹芯板起过滤作用的过滤膜组成的塑料排水板最为常用.大量超软土真空预压的工程实践表明，在预压过程中紧靠塑料排水板的周围会出现渗透性极低的“土桩”^[6-9]，这些“土桩”本身比周围的超软土具有更高的强度和密度，它的存在阻隔了“土桩”外围土体中水的流入，影响了这些区域土体的固结，有时甚至无法再继续固结而导致地基处理失败.针对这一“土桩”问题，近年来，一些学者进行了相应的研究：唐彤芝等^[7]通过现场试验认为吹填淤泥土颗粒细、结构性强度低、含水率高、流动性大是形成淤泥抱团(即“土桩”)的主要原因.陈平山等^[9]通过现场试验和数值分析，认为“土桩”是排水板周围土体细颗粒在渗透力作用下向排水板聚积的结果，其形状为倒锥形.武亚军等^[10]通过真空预压室内模型试验研究认为“土桩”是土颗粒在排水板处聚积所致.沈杰等^[11]通过室内模型试验，探讨真空荷载沿径向上的传递规律并分析造成排水板周围“土桩”现象的成因.以上关于“土桩”的研究，通常是采用常规测试技术获得真空预压前后排水板周围土体的颗粒分布、压缩及渗透性指标，并对这些指标的时空分布特点进行分析来推测“土桩”的大致形状及形成原因.但常规测试技术目前还存在一定的不足：常规测试技术在取样过程中会对土样产生扰动; 取样点是离散的，获取的数据不连续; 取样过程是不可逆的，导致每次取样所获得的并不是同一土样，这些不足会使试验结果产生一定的误差，因此，寻求一种更加直观的、非侵入式的研究手段尤为重要.

采用物理力学性质与天然土体相近的透明材料来模拟天然土，开展土体内部变形可视化物理模型试验研究逐渐成为一种趋势^[12].Iskander、Liu等^[13-14]采用无定形二氧化硅和与其折射率相匹配的孔隙液体制配出透明土，并对其进行了室内土工试验，发现其岩土工程性质与黏土相似.White等^[15]将PIV技术运用在岩土工程中，并通过试验对透明土体的变形进行了研究，表明PIV技术可以非常精确地观测土体内部变形.孔纲强等^[16]基于透明土材料和PIV技术，对扩底楔形桩、楔形桩和等截面桩的水平向受荷过程进行了模型试验研究，获得了其周围土体位移场的变化规律.宫全美等^[17]通过常规固结试验、直剪固快强度试验以及室内压板载荷试验验证了无定形硅粉类透明土在模型试验中模拟黏性土的可行性.夏元友等^[18]采用透明土材料和粒子图像测速法，通过物理模型试验研究了锚杆拔出机理.Wang等^[19]采用透明土进行模型试验，研究了基坑降水中防渗墙和抽油井的耦合效应.目前尚未见有学者采用透明土模拟超软土，实际上超软土的固结特性与黏土基本一致^[20]，因此，超软土在一定程度上可以理解为是含水率更高、黏粒含量更高、压缩性更高且强度更低的黏性土.

综上，将透明土试验技术应用于超软土真空预压模型试验，通过高清数码相机对真空预压全过程进行连续拍摄，得到真空预压过程中排水板周围土体变形的位移场，从而初步探讨“土桩”的形成机理.

1 真空预压模型试验 1.1 透明超软土制配及基本性质

本文所用透明超软土，采用无定形二氧化硅模拟“土”颗粒，采用与无定形二氧化硅折射率一致的混合矿物油模拟孔隙“水”.混合矿物油由正十二烷与二号白矿物油按照1: 11的质量比调配而成，其折射率为1.447 6.“土”颗粒的比重为2.11，粒径分布曲线如图 1所示.为了方便对比，图 1还给出了天津和温州超软土^[21-22]的粒径分布曲线.可以看出，透明超软土的粒径分布曲线大致位于两地超软土的粒径分布曲线之间，说明透明超软土的粒径分布在天然超软土粒径分布的变化范围内.以“土”颗粒、孔隙“水”按质量比1: 7的比例制配透明超软土，得到孔隙比e为2.47的试样(由无定形二氧化硅材料制配成的透明土的孔隙比计算方法, 可查阅相关研究成果^{[13-14, 17]}).借鉴Iskander等^[13]透明土制配方法，本文透明超软土试样的制配步骤为：1)将质量比为1: 7的“土”颗粒与孔隙“水”混合并置于密封桶中，使用固定转速的搅拌器搅拌5~10 min，至均匀黏稠状; 2)采用真空饱和装置，去除试样中由于搅动等因素产生的气泡，这一过程需要20~24 h; 3)对饱和后的试样，用手动搅拌器搅拌20~30 min，手动搅拌的优点是在保证试样整体均匀性的同时较大程度地避免在试样中引入气泡.固结后的透明超软土试样如图 2所示，透过4 cm厚的土样可以清晰地看到其下方的多孔板.

由无定形二氧化硅材料制配成的透明土的物性指标，可查阅相关文献研究成果^{[13-14, 17-18]}; 本文试验所用透明超软土的主要物性指标见表 1.为了方便对比，表 1还列出了温州和天津超软土的物性指标^[23-24].可以看出，透明超软土的初始孔隙比、压缩系数和压缩模量均在温州超软土对应指标的变化范围内; 透明超软土与天津超软土相比，其初始孔隙比较天津超软土的初始孔隙比最小值低0.01;压缩系数在天津超软土的变化范围内; 压缩模量比天津超软土的压缩模量最大值高0.11 MPa.

图 3为透明超软土的固结系数随荷载变化曲线，可以看出，透明超软土的固结系数随荷载的增大持续增大，这与张明等^[25]通过吹填超软土室内固结试验得到的规律基本一致.从数值上看，透明超软土的固结系数比吹填超软土的固结系数大，这与透明超软土本身较高的渗透性有关.图 4为透明超软土渗透系数k与孔隙比e的半对数关系曲线，曲线近似直线，这与Taylor等^[26]建立的反映黏性土渗透性的e-logk渗透模型一致.随着孔隙比的变化，透明超软土的渗透系数在3.3×10^-7~1.79×10^-5 cm/s变化，数值上约比相同孔隙比下吹填超软土的渗透系数高1~2个数量级^[27].

由上述可知，透明超软土的固结系数和渗透系数比天然超软土的高1~2个数量级，即透明土的排“水”固结速度较天然超软土更快，在荷载下完成变形的时间更短; 但是，透明超软土的初始孔隙比和压缩性指标均处于天然超软土对应指标的变化范围内，在相同外荷载作用下，透明超软土的压缩变形规律与天然超软土是相近的，因此，认为尽管透明土的诸多性质与天然超软土存在较大差异，利用透明超软土与天然超软土相似的压缩变形规律，从定性的层面上进行真空预压的“土桩”形成机理研究是可行的.

1.2 试验器材及步骤

模型试验器材包括模型箱(长50 cm×高40 cm×宽10 cm)、排水板(长15 cm×宽1 cm)、土工布及密封膜、抽滤瓶、空压机、真空发生器、高清数码相机及计算机等，具体如图 5所示.

将制配好的透明超软土置于模型箱中，静置48 h后排出上清液.为了利用PIV技术获取土体内部的位移场，在透明超软土中放入示踪颗粒(直径3 mm的塑料珠，可悬浮在透明超软土中)，使所有的示踪颗粒间隔均匀地布于同一个平面内(间隔2 cm)，该平面垂直于模型箱的宽度方向，与模型箱正面的距离为5 cm，称为示踪颗粒面.之后将排水板放置在示踪颗粒面的中间位置并垂直于示踪颗粒面，排水板上端露出透明超软土表面，再在其上铺设土工布和密封膜并在周边密封，之后连接排水板与真空系统开始抽真空.抽真空过程中对真空度、排“水”量进行监测，由于模型箱厚度较小，通过调压装置将真空度控制在60 kPa.以排“水”速率低于30 mL/h为卸载标准，该排水速率下模型箱内土体的表面沉降速率为1.44 mm/d，满足《港口工程地基规范》 (JTS 147-1—2010)的卸载标准，抽真空13 h后卸载.土体的变形过程用佳能D60数码相机记录，相机固定在三脚架上，通过电脑控制相机拍摄可得到真空预压过程中的一系列图像.通过PIV技术和数据可视化软件Tecplot对得到的图片进行处理，可获得真空预压过程中排水板周边土体的位移矢量图与位移等值线图.

2 试验结果与分析 2.1 排“水”速率

图 6为真空预压过程中的排“水”速率曲线.排“水”速率随着抽真空时间逐渐减小，大致可分为3个阶段，分别为高峰期(1 000 mL/h以上)、平稳期(1 000~400 mL/h)和缓慢期(400 mL/h以下)，3个阶段的持续时间分别为1，4和8 h.上述透明超软土排水过程分为3个阶段的试验现象与文献[22, 28]中基本一致，说明透明超软土的真空固结排水过程与天然超软土具有一定相似性.

2.2 排水板周围土体位移矢量图

对拍摄得到的一系列图片进行处理分析(为了方便分析，所有图片仅取右半部分)，发现排水板周围的土体位移在排“水”高峰期、平稳期和缓慢期呈现出不同的阶段性特征，将上述3个阶段内的土体位移用箭头矢量图来表示，列于图 7(a)、(b)和(c)，竖向排水板与土工布以黑色线表示，密封膜以黑色实线表示.由图 7可知，在排“水”高峰期，排水板周围约18 cm范围内的土体向着排水板斜向下位移，位移区域大致呈漏斗状，使得排水板周边土体的表面沉降在真空预压的初期高于稍远处.在排“水”平稳期，距排水板大约3 cm以内的土体位移以竖向向下为主，距排水板约3 cm以外的土体存在明显的斜向下位移，且位移的竖向分量随离开排水板距离的增加逐渐增大，使得土体表面产生沿着水平向上的不均匀沉降.在排“水”缓慢期，排水板周围24 cm范围内的土体均竖向向下位移，位移大小随着离开排水板水平距离的增大逐渐增大，导致距离排水板较远区域的土表面沉降比近处更明显.

2.3 排水板周围土体位移等值线图

为了更直观地表示排水板周围土体位移场，将结果绘制成等值线图，如图 8所示.考虑到边界效应，取距离排水板、密封膜最低点及模型箱边界各2 cm的方形区域作为观察区域.

图 8(a)和(d)分别为排“水”高峰期土体的水平和竖向位移等值线图，可以看出，土体的最大水平位移位于离开排水板约6 cm的浅层土体中，土体的竖向位移大致随离开排水板距离的增加逐渐减小.分析认为，透明超软土的初始含“水”率很高，抽真空开始后，大量自由“水”在压力差的作用下向排水板渗流排出，使得土体的浮容重变为有效容重后产生竖向压缩，靠近排水板位置土体的自由“水”含量降低得最多，因此，产生比排水板远端土体更大的竖向沉降; 孔隙“水”在渗流过程中，对“土”颗粒产生渗流力，使“土”颗粒沿着渗流方向移动，随着离开排水板距离的增加，渗流力逐渐衰减，“土”颗粒移动趋势也随之减弱，因此，排水板近端产生比远端更大的水平位移; 土颗粒会在自重下沉积，相互接触形成土骨架，浅层土体由于较低的上覆土重，形成的土骨架更松散，因此，在渗流力作用下产生较深层土体更加明显的渗流压密现象.

图 8(b)和(e)分别为排“水”平稳期土体的水平和竖向位移等值线图，土体的水平位移最大值位于距排水板约10 cm的中层土体中，土体的竖向位移随着离开排水板距离的增大逐渐增大.分析认为，随着抽真空过程的推进，排水板周边一定范围内的土体中已经形成了结构较为稳定的“土”骨架，渗流力得以向排水板远端的土体中传递，因此，较远处的土体产生渗流压密现象; 随着土中自由“水”含量的降低，膜上的“真空荷载”由初始的孔隙“水”承担转换为由“土”骨架承担，土体在自重和“真空荷载”的作用下产生竖向压缩，排水板远端的土体在渗流力的作用下趋于松散，因此，其竖向压缩量较排水板近端土体的竖向压缩量更大.

图 8(c)和(f)分别为排“水”缓慢期土体的水平和竖向位移等值线图，可以看出，土体的竖向位移随着离开排水板距离的增大逐渐增大，水平位移很小且无明显规律，可能是由土体竖向压缩过程中的水平向变形引起的.分析认为，在抽真空的后期，土体中的大部分孔隙“水”已经被排出，渗流力已不足以使土颗粒产生位移; 排水板远端土体较排水板近端土体仍然更为松散，因此，产生的竖向压缩量更大.

3 “土桩”形状及形成机理 3.1 “土桩”形状确定方法

PIV技术是将土体变形前后摄取的灰度图像分割成许多网格，每一网格称之为Interrogation块.将变形前任一Interrogation块与变形后灰度图像进行全场匹配或相关计算，根据峰值相关系数确定该Interrogation块在变形前后的位置，由此可以得到该块的平均位移，对变形前所有Interrogation块进行类似运算，就可得到整个位移场^[29].图像匹配的标准关联函数为

$ C(\Delta x, \Delta y) = \frac{1}{{MN}}\sum\limits_{m = 0}^{M - 1} {\sum\limits_{n = 0}^{N - 1} f } (m, n)g(m + \Delta x, n + \Delta y). $

式中：M，N为图像块的长、宽，f为t₁时刻图像中某图像块中心点坐标(m，n)处的灰度值分布函数，g为t₂时刻图像中该图像块中心点坐标(m+Δx，n+Δy)处的灰度值分布函数，Δx，Δy分别为x，y方向的位移增量.基于上述PIV技术的原理，可以求得“土桩”外包线，具体为：在某一深度下，对各相邻土体单元作Δx的差值，总能获得一个点，该点两侧单元的Δx差值为零，为方便讨论，将此类点称为位移差值零点，位移差值零点在深度方向的连线则称为位移差值零线.在位移差值零线以左的区域内，土体的Δx随着离开排水板距离的增大逐渐增大，即该区域的土体产生水平向的压缩，形成渗流压密区.在位移差值零线以右的区域内，土体的Δx随着离开排水板距离的增大逐渐减小，即该区域的土体在水平方向趋于松散，形成松散区.位移差值零线与排水板之间的土体实际上就是“土桩”，位移差值零线即为“土桩”的外包线.如图 9所示，“土桩”的半径由浅至深逐渐变小，最大半径约为11 cm，本文关于“土桩”形状的发现与陈平山等^[9]的研究结论相似.陈雷等^[30]通过真空预压室内模型试验，发现排水板周围土体的压缩系数较小，存在半径约为10 cm的密实区，这与本文关于“土桩”半径的发现接近.总体来说，透明土模型试验能较好地模拟真空预压过程中的“土桩”现象.

3.2 “土桩”形成机理解释

本文中的模型试验，笔者采用同样的试验设备和试验方法，在相同的环境下短期内进行过多次.每次试验结果均会出现“土桩”现象; 每次试验中排“水”速率随时间的变化规律、“土桩”形状随深度的变化规律都是一致的，因此，该模型试验的重复性和可靠性可以保证.

为了更好地解释“土桩”的形成机理，绘制一定深度处孔隙水的渗流示意图，如图 10所示.真空预压开始后，由于密封膜的存在，真空度向排水板中传递，引起排水板内部孔隙水压力的消散.密封膜使加固土体与大气隔绝，因此，排水板内的孔隙水压力p₁最大可降低至绝对压力零点^[30]，这使得土体内与排水体之间产生孔隙水压力差Δp₁，孔隙水在Δp₁的作用下向排水板方向渗流，降低渗流路径上的孔隙水压力，逐渐形成p₂、p₃和p₄.p₂、p₃和p₄的存在又形成了Δp₂和Δp₃，上述过程不断循环，在土体中逐渐形成一个相互连通的渗流网络，该渗流网络持续工作，降低整个土体内的孔隙水压力，达到土体固结的目的.在渗流过程中，孔隙水会对土颗粒施加渗流力，使土颗粒产生沿着渗流方向的位移，吹填超软土(尤其指新吹填超软土)中有大量细小颗粒处于悬浮状态，颗粒间结构强度尚未完全形成，在渗流力作用下很容易产生沿着渗流方向的位移.因此，大量的土颗粒聚积在排水板周边，并随着Δp由近至远向远处扩散，越来越多的土颗粒聚积在排水板周边.聚积在排水板周边的这部分土体，经过水平的渗流压密和竖向的固结压缩，逐渐形成“土桩”.随着离开排水板距离的增加，Δp逐渐衰减，土颗粒的移动趋势随之减弱，土体在水平方向上趋于松散，由此形成“土桩”范围以外的“软弱带”.由上所述，“土桩”范围内土体是经过聚积和渗流压密的，因此，其密度较“软弱带”范围内土体的更高，压缩性较“软弱带”范围内土体的更低.两块区域内的土体在自重和“真空荷载”作用下产生的压缩量不同，在土体表面形成“桩头”和“弧形凹槽”.图 11为真空预压试验过程中拍摄所得图片，图中的“桩头”和“弧形凹槽”现象明显.

试验结果表明，“土桩”半径随深度的增加有减小趋势，产生这种现象的原因主要有两个：其一，浅层土体比深层土体的细颗粒含量更高，且颗粒间结构强度尚未完全形成，在Δp作用下会产生比深层土体更加明显的沿着渗流方向的位移和渗流压密.其二，随着土层深度的增加，排水板内部的真空度存在衰减^{[4, 11]}，排水板内部的孔隙水压力消散程度也随之下降，即图 10中p₁的值随着深度增加有所提高，因此，孔隙水压力差Δp₁，Δp₂和Δp₃随深度增加有所下降.土颗粒的径向位移以孔隙水施加的渗流力为动力，真空度沿深度方向的损失导致了排水板径向上孔隙水压力差的减小，进一步导致了渗流力的衰减，因此，土颗粒的水平位移趋势随深度增加逐渐减弱，“土桩”的半径及土体的径向固结度也随之下降.实际工程中，深度每增加一米，排水板内的真空度衰减就可达10 kPa^{[4, 11]}，因此，真空度随深度的衰减，是“土桩”半径及土体径向固结度随深度发生变化的重要影响因素.

4 结论

1) 基于透明土试验技术和PIV技术，实现了真空预压过程土体内部位移的非侵入研究，相对常规测试技术，对土样不存在扰动且可以获得连续的位移场，具有经济合理、简便有效等优点.

2) 真空预压的影响范围随着抽真空时间逐渐扩大，模型箱内土体位移在真空预压过程中不同阶段呈现出不同的阶段性特征，导致土体表面沉降存在阶段性特征; 在排“水”高峰期，土体产生朝着排水板方向的斜向下位移; 在排“水”平稳期，土体的位移特征在离开排水板不同距离有较大差异; 在排“水”缓慢期，土体几乎仅有竖向位移，且位移随着离开排水板水平距离的增大逐渐增大.

3) 排水板内真空度的存在，使得周围土体内部形成一个相互连通的渗流网络，土颗粒在渗流力作用下不断向排水板处聚积; 土颗粒向排水板处聚积，并经过水平向渗流压密和竖向固结压缩，是形成“土桩”而外围土体强度过低的主要原因.

4)“土桩”侧面形状显示，“土桩”半径从上到下有减小的趋势，但本试验中模型箱底边界的存在，限制了底部土体的侧向位移，为了进一步研究“土桩”沿竖向的变化特征，需要在下一步研究中加大模型箱尺寸.

此外，透明土本身的性质限制了试验容器的大小，由此产生边界效应对试验结果的影响，需要进一步研究.