Mechanism of seepage control of pulsating grouting in completely weathered granite stratum
-
摘要: 针对全风化花岗岩地层的基础防渗加固,传统灌浆工艺限制于结构松散难以成孔起压,加剧了地基处理的难度。提出采用“钻灌一体,脉动灌浆”技术工艺；并利用数值模拟实现了不同脉动周期压力、钻灌时间控制参数下浆液渗透扩散规律及脉动周期压力下地层应力-应变的过程分析,通过脉动和稳压灌浆的浆液扩散与防渗影响范围对比,揭示了脉动灌浆防渗控制机理,得出了灌浆参数；结合黏土水泥浆灌浆材料进行了现场工程试验,现场试验验证了控制参数的合理性。试验表明,1.5～2.0 MPa脉动灌浆压力下,单排孔布设间距1 m,常规吕荣试验检查孔段次统计,透水率<2 Lu的占70%；疲劳吕荣试验透水率<3 Lu的比例为97%,透水率稳定区间为0.5～2.4 Lu；破坏吕荣试验透水率稳定区间为2～8.5 Lu。灌浆后全风化、强风化地层的岩体完整性均有不同程度的提高,声波提升幅度范围为17.3%～52.5%。检查孔取出的芯样较完整,芯样抗压强度平均达7.3 MPa,且灌浆过程地层抬升小。研究成果对于全风化花岗岩以及同类地层具有较强的应用性,为此类地层防渗处理提供了一个可借鉴的工程案例。Abstract: For the foundation reinforcement of the completely weathered granite stratum, the traditional grouting technology is limited to the loose structure, and it is difficult to form holes and to increase pressures, which aggravates the difficulty of foundation treatment. The technology of "drilling and grouting simultaneously, pulsating grouting" is put forward. The numerical simulation is carried out under different pulsating pressures and drilling and grouting time. The analysis of stress-strain of the stratum under pulsating pressures is realized. The control mechanism of pulsating grouting reinforcement is revealed by comparing the influence ranges of slurry diffusion of pulsating grouting and stabilized grouting. The field tests are carried put on the clay-cement slurry to verify the rationality of the control parameters. The numerical simulation results show that under the pulsating grouting pressures of 1.5~2.0 MPa and the spacing of holes of 1 m, according to the statistics of the test holes examined by the routine Lugeon tests, 70% of the sections with the permeability is smaller than 2 Lu. By the Lugeon tests on fatigue, the proportion of permeability smaller than 3 Lu is 97%, and the stable permeability is 0.5 ~ 2.4 Lu. The stable permeability by the Lugeon tests on damage is 2 ~ 8.5 Lu. After grouting, the strata integrity is improved to different degrees, and the acoustic wave elevation range is 17.3%~52.5%. The compressive strength of the core samples reaches 7.3 MPa on average, and the stratum uplift is small during grouting. The results of the study have a better application to the completely weathered granite and the similar strata, and provide a reference case for the reinforcement treatment of such strata.
-
Keywords:
- pulsating grouting /
- granite /
- Lugeon test /
- reinforcement mechanism /
- complete weathering
-
0. 引言
持水特征曲线指土体持水状态（可用体积含水率、质量含水率或饱和度表达）与土中水的能量状态（可用吸力、势能或相对湿度表达）之间的关系[1],描述了土中水（液相）与土（固相）间的能量平衡状态[2]。
土体持水的物理化学机制有毛细、吸附与渗透。毛细机制与弯液面上的表面张力有关,受Young- Laplace方程控制（土体孔径越小其毛细吸力越大）,当土中水能量状态低至其空化点以下,毛细水汽化,相应毛细机制消失。吸附机制指（黏）土颗粒表面由于van der Waals力、双电层力、阳离子水化等造成的短程吸附作用导致的土中水能量降低,可用Kelvin方程描述。通常将毛细与吸附机制导致的土中水相对自由水状态而言势能的降低称为基质势（用负值表达）；将对应土体固相的持水/吸水能力（用正值表达）称为基质吸力；二者数值相等,符号相反。此外,渗透机制指溶解盐会降低土中水的能量状态,其能量降低程度与溶质类型和浓度相关。
SWRC表征土体持水能力,原本用于农业领域,随后的研究表明非饱和状态下的渗流、强度与变形等土力学行为[1]乃至土的成分与结构等土质学特征[2]均与SWRC密切相关,因而成为描述非饱和状态下土体性状的关键指标与工具,在非饱和土力学研究与工程应用中具有基础的、特殊的、关键的作用与位置,其研究成果亦非常丰富[1-3],但大部分SWRC[1-3]由轴平移方法获得,受陶土板进气值限制,其基质吸力目前仅能控制在1500 kPa以内,根据Kelvin方程（式（1））可知其对应的相对湿度RH为98.9%以上（环境温度25℃时）。然而,自然界中RH变化范围远比轴平移方法所能控制的RH范围大；例如中国主要城市各月RH统计[4]显示,最低14%（拉萨,2月）,最高90%（长沙,10月与11月；海口,1月）；因而高吸力（RH低于95%或吸力高于7.1 MPa）下的土体持水特征研究具有特别的现实意义。
高吸力下土体持水特征测试常用方法主要有：滤纸法[5-8]（适用于10～100 MPa）、冷镜湿度计法[9]（适用于0～300 MPa）、饱和盐溶液等压湿度控制法[10-12]（适用于3～400 MPa）、干湿空气混合湿度控制法[13]（适用于7～600 MPa）等。前两种为吸力量测方法,难以控制吸力；后两种为吸力控制方法,适于测定完整的SWRC,但蒸汽平衡历时漫长,如文献[14]采用饱和盐溶液等压湿度控制法测定3.3～359.1 MPa吸力范围的膨胀土SWRC耗时319 d,阻碍了高吸力下土体持水特征的研究与应用。
近年来出现的动态露点等温线方法[15-16]使得上述情况有所改观,该方法能够在较短时间（几天内）获得相对湿度在3%～95%范围的SWRC。为了给高吸力下相关土力学工程实践提供依据和参考,本文采用该方法：①对从黏土质砂到极高塑性黏土的9种土样开展高吸力下吸脱湿过程持水特征测试；②采用RH = 95%下的吸湿质量含水率w95量化高吸力下持水能力,采用吸脱湿SWRC间滞回区域面积HHA量化高吸力下SWRC滞回效应强弱,采用同一吸力下吸脱湿SWRC含水率之差最大值Δwmax量化水力滞回导致的含水率偏差范围；③探讨w95、HHA、Δwmax与液限wL、塑性指数Ip、阳离子交换量CEC、比表面积SSA等土性参数间量化关系；④构建修正的Fredlund-Xing模型描述高吸力下吸脱湿过程SWRC。
1. 研究方案
1.1 试验土样
试验采用Wyoming膨润土、宁明膨胀土、荆门黄褐色膨胀土、Denver黏土岩、荆门棕褐色膨胀土、武汉黏土、三门峡粉质黏土、郑州粉土和开封黄河砂共9种土样,其物理性质指标测试结果见表1。
表 1 试验用土的物理性质指标Table 1. Physical property indexes of soils土样名称 颗粒相对质量密度[17] 颗粒所占比例[17]/% 液限[17]/% 塑限[17]/% 塑性指数 塑性[18] USCS 定名[19] 比表面积EGME[20]/(m2·g-1) CEC(NH4+)[17]/( mmol·kg-1) 自由膨胀率[21]/% 膨胀潜势[21] >0.075 mm 0.005~0.075 mm 0.002~0.005 mm <0.002 mm Wyoming膨润土 2.70 5.1 14.2 3.8 77.0 217.7 33.1 184.5 极高 CH 500.1 621 775 强 宁明膨胀土 2.73 5.8 2.9 0.3 91.1 164.6 37.6 127.0 极高 CH 393.0 508 655 强 荆门黄褐色膨胀土 2.75 2.1 47.4 21.5 29.0 62.9 25.5 37.4 高 CH 236.5 309 75 中 Denver黏土岩 2.72 8.5 44.2 20.6 26.7 46.0 23.2 22.8 中等 CL 88.7 160 38 — 荆门棕褐色膨胀土 2.72 11.6 42.4 26.7 19.3 41.8 20.7 21.1 中等 CL 110.4 254 40 弱 武汉黏土 2.73 5.5 63.9 9.6 21.0 40.3 18.2 22.1 中等 CL 112.5 154 39 — 三门峡粉质黏土 2.72 0.5 82.8 5.7 11.1 35.2 19.1 16.1 中等 CL 108.5 149 30 — 郑州粉土 2.70 8.7 77.4 2.5 11.4 26.1 14.8 11.3 低 CL 27.7 50 9 — 开封黄河砂 2.68 62.6 36.0 0.2 1.3 — — — — SC 7.5 16 — — 1.2 蒸汽吸附分析仪工作原理
美国Meter公司生产的蒸汽吸附分析仪（Vapor Sorption Analyzer, VSA）见图1,用于测试动态的或平衡的水分吸附等温线（moisture sorption isotherms）,即给定温度下的SWRC。
VSA采用动态露点等温线（dynamic dew point isotherm, DDI）方法获得动态水分吸附等温线；采用控制湿度下测试样质量的方法（dynamic vapor sorption, DVS）获得静态（平衡）水分吸附等温线。其工作原理是：DDI方法既不控制含水率亦不控制RH,而是在用湿气流湿化/干气流脱湿过程中通过高精度天平测试样质量变化、通过冷镜露点法测RH得到动态水分吸附等温线。DVS方法是在控制RH条件下跟踪试样质量变化以获得静态（平衡）水分吸附等温线,每级湿度下的稳定（平衡）标准可用预设时间间隔或预设试样质量变化幅度。Likos等[15]与Arthur等[22]试验证实了DDI方法的精度与可靠性,系统试验研究表明DDI方法与DVS方法获得的SWRC无明显偏差。
VSA相对湿度控制范围是3%～95%（25℃下对应吸力范围为7.1～482.9 MPa）,试样质量建议范围为500～5000 mg,试样温度可控范围为15℃～60 ℃。
1.3 高吸力下持水特征试验方案
采用VSA测定9种土样高吸力下SWRC。将土样风干碾散过0.5 mm筛,置于烘箱中烘至衡量（105 ℃下烘干48 h）,取出置于干燥器内冷却至室温,每次取约1000 mg试验。试验方案与步骤如下：
（1）仪器准备
将VSA在工作台上调平,安装干燥管,水箱中注入去离子水,检查天平,标定相对湿度。
（2）设定SWRC量测步骤
鉴于采用DDI方法与DVS方法测土样SWRC具有很好的一致性[22],且前者测试历时约为后者的1/25[15],因此采用DDI方法。由于试样是完全烘干后的,质量含水率w为0,第1步设置初始RH为3%,最终RH为95%,用DDI方法采集吸湿过程w数据RH,即土样SWRC吸湿边界线；第2步设置初始RH为95%、最终RH为3%,目的是获得土样SWRC的脱湿曲线；由于脱湿曲线起始点为RH=95%对应点,所以并非脱湿边界线。此外,设定样品室温度为25℃。
（3）测定SWRC
将试样置入样品杯,放入试样室开始测试,测试由VSA自动执行。试验结束后,取出试样杯称试样质量,烘干后复测干土质量,并在计算机终端下载SWRC试验数据。
2. 试验结果分析
由于Kelvin方程（式(1)）确定了吸力ψ与RH间关系,VSA试验结果既可用w-RH关系亦可用w-ψ关系表达；鉴于非饱和土力学中常用w-ψ关系,因此将9种土样w-ψ关系表达的SWRC绘制在图2中,将每种土样试验历时,试验数据点数,吸湿起始点、吸湿结束点（脱湿起始点）、脱湿结束点对应的RH,ψ,w值列在表2中。
表 2 高吸力下SWRC的吸湿/脱湿起始点与结束点Table 2. Start point/end point of adsorption/desorption of SWRCs under high suctions土样名称 试验历时/min 试验数据点数 吸湿起始点 吸湿结束点（脱湿起始点） 脱湿结束点 RH ψ/kPa w/% RH ψ/kPa w95/% RH ψ/kPa w/% Wyoming膨润土 3135 203 0.0375 451647 0.67 0.9457 7686 21.46 0.0350 461106 0.87 宁明膨胀土 3777 205 0.0324 471755 0.76 0.9472 7468 19.45 0.0343 463909 1.02 荆门黄褐色膨胀土 2982 204 0.0297 484119 1.28 0.9456 7700 13.40 0.0331 468823 1.96 Denver黏土岩 2113 201 0.0313 476506 0.79 0.9476 7410 6.79 0.0338 465943 1.14 荆门棕褐色膨胀土 2523 212 0.0299 483240 0.89 0.9474 7439 9.63 0.0245 510367 1.28 武汉黏土 2343 202 0.0294 485114 0.44 0.9455 7715 7.37 0.0326 470929 0.91 三门峡粉质黏土 2186 203 0.0244 510892 0.43 0.9478 7380 6.19 0.0327 470498 0.97 郑州粉土 1787 197 0.0246 509791 0.42 0.9500 7061 2.34 0.0245 510347 0.52 开封黄河砂 1295 148 0.0245 510383 0.65 0.9519 6786 1.28 0.0247 509237 0.67 ψ=−RTνwlnRH, (1) 式中 ψ为吸力（kPa）；R = 8.314 J/(K·mol),为通用气体常数；T为热力学温度,样品室温度25℃对应T = 298.15 K；
νw = 1.8×10-5 m³/mol,为水的摩尔体积；RH为相对湿度,量纲为1。由表2可见,9种土样中宁明膨胀土测试历时最长,为3777 min,即2.6 d,且有205个数据点。一方面说明试验历时为文献[14]中饱和盐溶液等压湿度控制法测试历时（319 d）的1/122,极大缩短了测试时间；另一方面说明数据点丰富,能够更好描述土样的持水特征。
从以下3方面探讨高吸力下的SWRC：①反映土样亲水性强弱的持水能力特征点w95；②描述水力滞回效应强弱的SWRC滞回区域面积HHA；③反映水力滞回导致的含水率偏差范围的SWRC滞回效应最大值Δwmax。
2.1 高吸力下持水能力特征点w95
由图2与表2可见：9种土样从完全干燥状态吸湿结束时（控制RH为95%）对应的质量含水率w95差别相当大,高至Wyoming膨润土的21.46%,宁明膨胀土的19.45%,表明持水能力/亲水性非常强,低至开封黄河砂的1.28%,表明其持水能力/亲水性非常弱。由此可见,w95作为RH = 95%下的吸湿质量含水率,是表征土样持水（吸水）能力的特征点,其数值反映了土样亲水性强弱。
将9种土样的w95与表1给出的CEC对应绘制在图3中,可见w95与CEC具有非常高的线性相关性,其线性拟合结果见式（2）,决定系数R2高达0.99。
w95= 0.035 CEC + 1.3。 (2) 鉴于《膨胀土地区建筑技术规范》[21]采用CEC分别为340,260,170 mmol/kg作为判断土样膨胀潜势强（CEC > 340 mmol/kg）、中（CEC =260～340 mmol/kg）、弱（CEC = 170～260 mmol/kg）的界限,将其代入式（2）得到采用w95判别土样膨胀潜势的相应界限分别为13.0%,10.2%,7.1%。即：w95 > 13.0%,判别为强膨胀土；w95=10.2%～13.0%,判别为中膨胀土；w95=7.1%～10.2%,判别为弱膨胀土；w95<7.1%,判别为非膨胀土。
下面对上述方法的可行性进行验证：表3为根据《膨胀土地区建筑技术规范》[21]测定的Denver膨润土、信阳黏土、商用高岭土的CEC、自由膨胀率及确定的膨胀潜势。采用VSA获得的w95亦列在表3中,Denver膨润土w95为22.47%,大于13.0%,根据w95判别为强膨胀土；信阳黏土w95为8.77%,在7.1%～10.2%,根据w95判别为弱膨胀土；商用高岭土w95为1.03%,小于7.1 %,根据w95判别为非膨胀土。表明采用w95作为膨胀潜势判别标准与采用CEC的判别结果具有很好的一致性,至于如何提高其判别精度,则应收集更多的试验结果,确定更为可靠的强、中、弱膨胀潜势的w95界限值。
表 3 3种土样膨胀潜势的判别指标Table 3. Eestimation indexes of swelling potential of 3 kinds of soils2.2 高吸力下SWRC滞回区域面积HHA
由图2可见,高吸力段(RH = 3%～95%),9种土样SWRC均呈现或强或弱的水力滞回行为,即相同吸力下脱湿路径上的w高于吸湿路径上的w。
SWRC水力滞回效应强弱可用脱湿与吸湿曲线间滞回区域面积HHA (hydraulic hysteresis area)描述。根据图2中的脱湿与吸湿曲线数据在Origin软件中数值积分获得9种土样滞回区域面积见表4。
表 4 高吸力下SWRC滞回效应特征参数（RH = 3%~95%）Table 4. Hysteresis properties of SWRCs under high suctions土样名称 滞回区域面积HHA/kPa HHA 相对值 Δwmax/% ψ(Δwmax)/kPa RH(Δwmax) wd(Δwmax)/% ww(Δwmax)/% w相对误差/% Δwmaxwd Δwmaxww Wyoming膨润土 13641 1.00 6.27 68737 0.6069 13.08 6.81 47.9 92.1 宁明膨胀土 10717 0.79 4.18 72744 0.5895 11.00 6.81 38.0 61.4 荆门黄褐色膨胀土 4512 0.33 1.77 37559 0.7612 9.87 8.10 17.9 21.9 Denver黏土岩 2367 0.17 0.69 37184 0.7633 4.80 4.11 14.4 16.8 荆门棕褐色膨胀土 3520 0.26 1.38 30993 0.7984 7.41 6.03 18.6 22.9 武汉黏土 3412 0.25 1.15 31077 0.7979 5.13 3.98 22.4 28.9 三门峡粉质黏土 3439 0.25 1.04 89267 0.5228 3.55 2.51 29.3 41.4 郑州粉土 1002 0.07 0.37 155335 0.3235 1.22 0.85 30.3 43.5 开封黄河砂 196 0.01 0.16 221700 0.1999 0.88 0.72 18.2 22.2 由图2与表4可见：9种土样HHA差别显著,对于强膨胀土,Wyoming膨润土高达13641 kPa,宁明膨胀土亦高达10717 kPa；对于中膨胀土,荆门黄褐色膨胀土为4512 kPa；对于弱膨胀土,荆门棕褐色膨胀土为3520 kPa；武汉黏土、三门峡粉质黏土亦达到3412,3439 kPa；而郑州粉土低至1002 kPa；开封黄河砂仅有196 kPa,是Wyoming膨润土的1/70,滞回效应可忽略不计。
以上数据表明：HHA与土样膨胀潜势正相关,将9种土样的HHA与表1给出的CEC、SSA对应绘制在双y轴坐标系中,见图4,可见HHA与CEC、SSA均具有非常高的线性相关性,其过坐标原点的线性拟合结果为
CEC = 0.049 HHA, (3) SSA = 0.037 HHA, (4) 式中,HHA与CEC线性相关的R2达到0.94,HHA与SSA线性相关的R2高达0.97。
CEC与SSA是强烈相关但有明显区别的土性指标。CEC指单位质量土样吸附的可交换阳离子总量,反映由于颗粒带负电荷导致的土样吸附阳离子的能力,更大的CEC意味着更强的膨胀潜势。SSA指单位质量土样具有的表面积,反映土样的颗粒级配与矿物成分,更大的SSA意味着更大的塑性,更强的膨胀潜势,更强的亲水性。HHA与CEC、SSA均正相关表明,滞回效应大小既受土样吸附阳离子能力的强烈影响,亦受土样比表面积控制,且从R2上看,其受土样比表面积影响更大。
将8种细粒土HHA与其液限wL、塑性指数Ip绘制在双y轴坐标系中,见图5。
![]()
图 5 8种细粒土的wL-Ip-HHA关系Figure 5. wL-Ip-HHA relationship of 8 kinds of fine-grained soils可见HHA与wL、Ip均正相关,过坐标原点的线性拟合结果分别为
wL= 0.015 HHA, (5) Ip= 0.012 HHA。 (6) 式中,HHA与wL线性相关的R2高达0.97,HHA与Ip线性相关的R2亦达到0.92,表明细粒土HHA与wL具有更高的线性相关性。由于wL、Ip代表土样塑性大小,表明HHA受土样塑性大小控制。
2.3 高吸力下SWRC滞回效应最大值Δwmax
SWRC滞回效应指水力滞回（脱/吸湿路径）对SWRC（此处是w-ψ关系）的影响,此处SWRC滞回效应最大值定义为相同ψ下,脱湿曲线对应含水率wd与吸湿曲线对应含水率ww之差的最大值Δwmax。根据该定义,在图2中搜索到每种土的Δwmax,及其对应的吸力ψ、RH、脱湿曲线含水率wd、吸湿曲线含水率ww,列入表4中。
由表4与图2可见,9种土样Δwmax存在较大差异：强膨胀土Δwmax相当大,如Wyoming膨润土为6.27%,宁明膨胀土为4.18%；中膨胀土,荆门黄褐色膨胀土为1.77 %；弱膨胀土,荆门棕褐色膨胀土为1.38%；中等塑性黏土,武汉黏土为1.15%,三门峡粉质黏土为1.04%,Denver黏土岩为0.69%；低塑性黏土,郑州粉土为0.37%；开封黄河砂仅为0.16%。对比表1可见,土样Δwmax值与wL,Ip,CEC,SSA均正相关,反映了水力滞回导致的含水率偏差范围。取Δwmax与相应wd、ww的比值作为SWRC的含水率相对误差,计算9种土样Δwmax/wd、Δwmax/ww结果列于表4与图6。
![]()
图 6 由于水力滞回导致的含水率相对误差Figure 6. Relative errors of water content due to hydraulic hysteresis表4与图6中,同一土样ww小于wd,因此采用吸湿曲线含水率ww作为分母的相对误差比采用脱湿曲线含水率wd作为分母的相对误差要大。由相对误差计算结果可见：在高吸力范围,不考虑水力滞回会导致持水特征曲线SWRC上质量含水率w最高92.1 %（Wyoming膨润土）的相对误差。
另一方面,如果工程中需要用当前含水率估算吸力,则同一w下,脱湿路径上吸力与吸湿路径上吸力的比值最高可达2.71（图2中Wyoming膨润土,w为6%,脱湿路径上吸力205959 kPa,吸湿路径上吸力76119 kPa,比值为2.71）,印证了高吸力下强膨胀土SWRC滞回效应之强烈。
2.4 高吸力下SWRC滞回行为机制探讨
由图2及2.2节、2.3节对SWRC滞回区域面积HHA与滞回效应最大值Δwmax的分析,可知高吸力下,9种土样中,高塑性（液限高、塑性指数大）、阳离子交换量大、比表面积大的土样滞回效应大,在此对高吸力下SWRC滞回行为机制展开探讨。
水力滞回行为机制解释主要有墨水瓶效应[1]与接触角效应[13],前者由于土中孔隙尺寸与形状（可用孔径r概化）的不均匀性所致,后者由于吸湿接触角
αw 通常大于脱湿接触角αd 所致[23]。墨水瓶效应与接触角效应均属毛细范畴,受Young-Laplace方程控制：ψc=ua−uw=2Tscosαr。 (7) 式中 ψc为毛细吸力,可用孔隙气压力ua与孔隙水压力uw之差表达；Ts为水的表面张力；
α 为接触角；r为毛细半径。当吸力增大、土中水能量状态低到其空化点以下,毛细水汽化,（土中毛细水比自由水空化点更低,但其对应的ψc或RH尚无定论,一般而言RH低于80%,毛细水难以存在[2]）,墨水瓶效应与接触角效应难以解释高吸力下的滞回行为。Lu等[24]给出了高吸力下SWRC滞回机制的一种解释：非膨胀黏土,脱湿与吸湿主要发生在颗粒表面,因而滞回效应较小。膨胀性土,脱湿与吸湿亦可发生于存在可交换阳离子的晶层之间,吸湿时阳离子水化吸附极性水分子,水分子需克服较大的,van der Waals力方能进入层间；脱湿时由于层间距离较大van der Waals力较弱,因此相同ψ下,吸湿路径对应较小的w,脱湿路径对应较大的w。在此结合文中数据探讨该解释的合理性：
高吸力下Wyoming膨润土HHA=13641 kPa,Δwmax=6.27 %,滞回效应最大；自由膨胀率为775 %,膨胀性最大；SSA高达500.1 m2/g,表明其包含了较大的层间比表面积[25],CEC高达621 mmol/kg,表明其可交换阳离子总量大；印证了吸湿状况下膨胀土层间阳离子水化受van der Waals力阻碍致其含水率较低的机制解释。
相对而言,郑州粉土,SSA为27.7 m2/g,CEC为50 mmol/kg,比表面积小,可交换阳离子总量小。一方面其自由膨胀率（9%）为Wyoming膨润土的1.2 %,表明其为非膨胀土；另一方面,其HHA（1002 kPa）为Wyoming膨润土的7%,Δwmax（0.37%）为Wyoming膨润土的6%,滞回效应很小。印证了非膨胀土脱湿吸湿主要发生在颗粒表面（而非层间）致其滞回效应较小的机制解释。
以上讨论明晰了2.2节中的论述：高吸力下滞回效应既受CEC的强烈影响,亦受SSA控制,且从R2上看,其受SSA影响更大。至于是否考虑滞回效应,应视土样与工况等具体情况而定,鉴于Wyoming膨润土具有最大的HHA,可采用其作为标准,其它土样的HHA相对值亦列入表4中,可见从开封黄河砂的0.01到宁明膨胀土的0.79,HHA相对值可作为实际工况中是否考虑滞回效应的依据。
3. 高吸力下持水特征的本构描述
基于对以上9种土样高吸力下SWRC规律性的认知,经尝试,建议了Fredlund-Xing（1994）模型[1]的修正形式描述高吸力下SWRC：
w(ψ)=C(ψ)wu{ln[exp(1)+((ψ−ψb)a)n]}m=[1−ln(1+ψ/ψr)ln(1+106/ψr)]wu{ln[exp(1)+((ψ−ψb)a)n]}m。 (8) 式中 w为质量含水率（%）；ψ为吸力（kPa）；w(ψ)为质量含水率表达的SWRC；C(ψ)为修正系数函数,能够迫使SWRC通过零含水率对应吸力106 kPa的点。4个参数a,n,m,ψr为Fredlund-Xing(1994)模型已有参数,a为与拐点位置相关的土性参数（kPa）；n为与拐点处斜率相关的土性参数,量纲为1,m为与曲率相关的土性参数,量纲为1,ψr为残余吸力（kPa）。
为描述高吸力下SWRC,将Fredlund-Xing模型中饱和质量含水率w0用wu替代,wu为高吸力下SWRC质量含水率上边界（最大值）,相应吸力为ψb。增加土性参数ψb是式（8）与Fredlund-Xing模型形式上的主要区别,ψb单位为kPa,物理意义是高吸力下SWRC吸力下边界。由此可见,修正的Fredlund-Xing模型共有6个参数：wu,ψb,ψr,a,n,m。为区分吸湿与脱湿过程,相应参数增加上标w（wetting,表示吸湿）或d（drying,表示脱湿）。
采用式（8）分别对9种土样高吸力下SWRC进行拟合,流程如下：
（1）采用表2给出的吸湿结束点（脱湿起始点）含水率w95值作为wu值,对应吸力值作为ψb值；wu与ψb为吸脱湿过程共用参数值。
（2）对于吸湿过程,采用表2给出的吸湿起始点吸力值作为ψr值；至此,吸湿过程所需确定的参数仅余a,n,m。依据式（8）在Excel中通过规划求解功能对吸湿过程SWRC数据进行非线性拟合获得a,n,m最优解以及拟合结果的决定系数R2。
（3）对于脱湿过程,采用表2给出的脱湿结束点吸力值作为ψr值；依据式（8）在Excel中通过规划求解功能对脱湿过程SWRC数据进行非线性拟合获得a,n,m最优解以及R2。
前8种土样采用上述流程获得的模型参数与R2见表5,可见8种土样,无论脱湿吸湿,修正的Fredlund-Xing模型再现结果的R2均在0.99以上。为节约篇幅,其修正Fredlund-Xing模型数值再现结果亦放在图2中,可见式（8）相当好地再现了高吸力下8种土样脱湿与吸湿过程SWRC。
表 5 高吸力下修正Fredlund-Xing模型参数Table 5. Modified Fredlund-Xing model parameters for 9 kinds of soils under high suctions土样名称 共用参数 吸湿过程 脱湿过程 wu/% ψb/kPa ψrw/kPa aw/kPa nw mw R2 ψrd/kPa aw/kPa nd md R2 Wyoming膨润土 21.46 7686 451647 248967 0.94 10.54 1.00 461106 90593 1.06 2.09 0.99 宁明膨胀土 19.45 7468 471755 245967 0.82 8.02 1.00 463909 90230 0.94 2.13 1.00 荆门黄褐色膨胀土 13.40 7700 484119 237692 0.61 4.34 1.00 468823 91452 0.82 1.90 1.00 Denver黏土岩 6.79 7410 476506 267853 0.57 4.14 1.00 465943 91075 0.66 1.79 1.00 荆门棕褐色膨胀土 9.63 7439 483240 87263 0.64 2.81 1.00 510367 86362 0.80 1.96 0.99 武汉黏土 7.37 7715 485114 81572 0.67 3.77 1.00 470929 80590 0.71 2.39 0.99 三门峡粉质黏土 6.19 7380 510892 213682 0.67 4.50 1.00 470498 91368 0.82 1.56 1.00 郑州粉土 2.34 7061 509791 30152 0.78 1.17 0.99 510347 31576 0.58 0.86 1.00 开封黄河砂 1.28 6786 — 76542 0.66 1.04 0.96 — 76542 0.66 1.04 0.96 对第9种土样开封黄河砂,由图2可见其SWRC滞回行为非常小,因此采用同一组模型参数描述其SWRC,此外不考虑式（8）中修正系数函数C(ψ),即C(ψ) = 1,对吸脱湿SWRC数据进行非线性拟合获得a,n,m最优解见表5,同时描述其吸脱湿SWRC的R2达到0.96,由图2的再现结果可见式（8）非常好地对其SWRC进行了本构描述。
此外,对比表5中吸脱湿过程a,n,m值可见,对同种土样,相对脱湿过程,吸湿过程a值较大说明其拐点位置靠后,n值较低说明其拐点处斜率相对较小,m值较大说明其曲率较大。
由此可见,文中建议的Fredlund-Xing模型修正形式有6个参数：wu,ψb,ψr值直接从实测SWRC数据确定,a,n,m值通过非线性拟合获得。Fredlund- Xing(1994)模型中4个参数（a,n,m,ψr）需通过非线性拟合获得,还有一个参数饱和质量含水率w0需试验测定（从实测SWRC数据确定）,实际上有5个参数。文中建议的修正形式仅增加了一个参数ψb,且ψb与wu、ψr直接从实测SWRC数据确定,只需非线性拟合3个参数（a,n,m）,且该修正形式能够很好地分别描述吸湿与脱湿过程中SWRC强烈的非线性行为。
4. 结论
（1）9种土样高吸力下吸湿与脱湿过程持水特征测试最长历时为2.6 d,有205个数据点,说明动态露点法极大缩短了高吸力下SWRC测试时间,且数据点丰富,能够更好描述土样持水特征。
（2）w95作为RH = 95%下的吸湿质量含水率,是表征土样持水能力的特征点,反映土样亲水性强弱,w95与CEC具有非常高的线性相关性,w95作为膨胀潜势的判别标准是可行的。
（3）SWRC水力滞回效应强弱可用吸脱湿曲线间滞回区域面积HHA描述,HHA与CEC、SSA均有非常高的线性相关性,表明滞回效应强弱既受土样吸附阳离子能力的强烈影响,亦受土样SSA控制,且从R2看,其受SSA影响更大。HHA与wL,Ip均正相关,表明HHA受土样塑性大小控制。
（4）SWRC滞回效应最大值可用同一ψ下脱湿与吸湿曲线对应w之差最大值Δwmax量化,土样Δwmax值与wL,Ip,CEC,SSA均正相关,反映了水力滞回导致的w偏差范围。高吸力范围,不考虑水力滞回会导致SWRC上w最高92.1 %的相对误差。如工程中需要用当前w估算ψ,则同一w下,脱湿路径上ψ与吸湿路径上ψ的比值最高可达2.71。
（5）建议采用Fredlund-Xing模型的修正形式描述高吸力下SWRC,方程有6个参数,wu,ψb,ψr值直接从实测SWRC数据确定,a,n,m值通过非线性拟合获得。对9种土样SWRC的修正Fredlund -Xing模型数值再现结果表明该方程可分别很好地描述吸湿与脱湿过程中SWRC强烈的非线性行为。
-
![]()
图 2 钻灌一体,脉动灌浆工艺流程图
Figure 2. Flow chart of “drilling and grouting simultaneously, pulsating grouting” technology
![]()
图 3 地层网格划分与随机孔隙率地层模型
Figure 3. Stratigraphic meshing and random porosity stratigraphic model
![]()
图 4 浆液相对质量密度1.5和1.6的黏度方程拟合
Figure 4. Fitting of slurry with specific gravity of 1.5 and 1.6
![]()
图 5 1 MPa脉动压力浆液扩散随时间变化图
Figure 5. Variation of slurry diffusion with time under pulsating grouting pressure of 1 MPa
![]()
图 6 1.5 MPa脉动压力浆液扩散随时间变化图
Figure 6. 1.5 MPa pulsating grouting pressure of slurry diffusion over time
![]()
图 7 2 MPa脉动压力浆液扩散随时间变化图
Figure 7. Variation of slurry diffusion with time under pulsating grouting pressure of 2 MPa
![]()
图 8 1.5 MPa稳压灌浆浆液扩散随时间变化图
Figure 8. Variation of slurry diffusion with time under stabilized grouting pressure of 1.5 MPa
![]()
图 11 全风化花岗岩粒度曲线
Figure 11. Grain-size distribution curve of completely weathered granite
![]()
图 15 检查孔J3疲劳吕荣透水率与检查孔J4破坏吕荣透水率
Figure 15. Permeabilities by Lugeon tests on fatigue of test hole J3 and Lugeon tests on damage of test hole J4
表 1 水泥黏土稳定浆液黏度随时间变化关系
Table 1 Variation of viscosity of clay-cement-stabilized slurry with time
(s) 相对质量密度 0 min 10 min 20 min 30 min 40 min 1.5 20 22 24 27 29 1.6 22 24 27 28 30 
下载: 导出CSV
表 2 全风化花岗岩基本物理参数
Table 2 Basic physical parameters of completely weathered granite
密度/(g·cm-3) 干密度/(g·cm-3) 含水率/% 液限/% 塑限/% 孔隙率/% 1.93 1.71 7.6 8.1 17.9 36
下载: 导出CSV
表 3 各水固比浆液性能参数
Table 3 Performance parameters of slurry with different water-solid ratios
水固比 黏土∶水泥 相对质量密度 初凝时间/min 析水率/% 流动度/cm 黏度/s 28 d抗压强度/MPa 0.8 1∶1 1.55 67 0 13.2 24 21.8 1.0 2∶1 1.48 92 0 19.7 21 5.2 1.2 2∶0.75 1.45 108 1 24.3 19 3.3 1.5 2∶0.5 1.41 121 1 27.0 18 1.5
下载: 导出CSV
表 4 声波检测数据记录表vs
Table 4 Record sheet of acoustic detection data vs
(km/s) 灌浆前后 岩石 J4 J1 J2 J3 灌浆后 全风化 1.10 1.18 1.67 1.35 强风化 2.15 2.78 2.72 2.29 灌浆前 全风化 0.84 0.97 0.79 1.10 强风化 1.66 2.31 1.28 1.87
下载: 导出CSV
-
[1] 王凯, 李术才, 杨磊, 等. 全风化花岗岩加固特性注浆模拟试验[J]. 天津大学学报(自然科学与工程技术版), 2017, 50(11): 1199-1209. https://www.cnki.com.cn/Article/CJFDTOTAL-TJDX201711013.htm WANG Kai, LI Shu-cai, YANG Lei, et al. Grouting simulation experiment on reinforcement characteristics of completely decomposed granite[J]. Journal of Tianjin University (Science and Technology), 2017, 50(11): 1199-1209. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TJDX201711013.htm
[2] WANG S, CHAN D, LAM K C. Experimental study of the effect of fines content on dynamic compaction grouting in completely decomposed granite of Hong Kong[J]. Construction and Building Materials, 2009, 23(3): 1249-1264. doi: 10.1016/j.conbuildmat.2008.08.002
[3] ARIAS D, PANDO L, LÓPEZ-FERNÁNDEZ C, et al. Deep weathering of granitic rocks: a case of tunnelling in NW Spain[J]. CATENA, 2016, 137: 572-580. doi: 10.1016/j.catena.2015.10.026
[4] 杨磊, 林荣峰, 李召峰, 等. 水泥浆液黏度对全风化花岗岩注浆加固效果的影响[J]. 中国公路学报, 2018, 31(10): 246-254. doi: 10.3969/j.issn.1001-7372.2018.10.024 YANG Lei, LIN Rong-feng, LI Zhao-feng, et al. Influence of grout viscosity on the grouting reinforcement effect of completely weathered granite[J]. China Journal Highway Transportation, 2018, 31(10): 246-254. (in Chinese) doi: 10.3969/j.issn.1001-7372.2018.10.024
[5] 袁敬强, 陈卫忠, 黄世武, 等. 全风化花岗岩注浆加固特性试验研究[J]. 岩石力学与工程学报, 2016, 35(增刊1): 2876-2882. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2016S1033.htm YUAN Jing-qiang, CHEN Wei-zhong, HUANG Shi-wu, et al. Experimental study on physico- mechanical properties of grouted completely weathered granite[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(S1): 2876-2882. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2016S1033.htm
[6] ZHANG D, FANG Q, LOU H. Grouting techniques for the unfavorable geological conditions of Xiang'an subsea tunnel in China[J]. Journal of Rock Mechanics and Geotechnical Engineering, 2014, 6(5): 438-446. doi: 10.1016/j.jrmge.2014.07.005
[7] 张聪, 梁经纬, 阳军生, 等. 堤坝脉动注浆浆液扩散机制及应用研究[J]. 岩土力学, 2019, 40(4): 1507-1514. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201904033.htm ZHANG Cong, LIANG Jing-wei, YANG Jun-sheng, et al. Research on the diffusion mechanism and application of pulsate grouting in embankment and dam. Rock and Soil Mechanics, 2019, 40(4): 1507-1514. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201904033.htm
[8] 张贵金, 梁经纬, 杨东升, 等. 松软地层脉动灌浆封孔浆体止浆机制初步研究[J]. 岩土工程学报, 2016, 38(2): 210-218. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201602004.htm ZHANG Gui-jin, LIANG Jing-wei, YANG Dong-sheng, et al. Mechanism of stemming materials for pulsating grouting in soft strata[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(2): 210-218. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201602004.htm
[9] SHI Y, YANG J, WANG S. Sub-horizontal reinforcement of weathered granite before tunneling beneath a spillway[J]. International Journal of Rock Mechanics and Mining Sciences, 2014, 72: 283-293. doi: 10.1016/j.ijrmms.2014.08.005
[10] 吕晓光, 王德发, 姜洪福. 储层地质模型及随机建模技术[J]. 大庆石油地质与开发, 2000(1): 10-13, 16-52. https://www.cnki.com.cn/Article/CJFDTOTAL-DQSK200001003.htm LU Xiao-guang, WANG De-fa, JIANG Hong-fu. Reservoir geological model and stochastic modeling technique[J]. Petroleum Geology and Development in Daqing, 2000(1): 10-13, 16-52. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DQSK200001003.htm
[11] 张聪, 梁经纬, 张箭, 等. 基于脉动注浆的宾汉流体渗透扩散机制研究[J]. 岩土力学, 2018, 39(8): 2740-2746, 2754. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201808004.htm ZHANG Cong, LIANG Jing-wei, ZHANG Jian, et al. Mechanism of Bingham fluid permeation and diffusion based on pulse injection[J]. Rock and Soil Mechanics, 2018, 39(8): 2740-2746. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201808004.htm
[12] 张聪, 阳军生, 谢亦朋, 等. 非均质软弱围岩隧道注浆加固圈分布特性[J]. 交通运输工程学报, 2019, 19(3): 58-70. https://www.cnki.com.cn/Article/CJFDTOTAL-JYGC201903007.htm ZHANG Cong, YANG Jun-sheng, XIE Yi-peng, et al. Distribution characteristics of grouting reinforcemen tring fortunnel with heterogeneous weak surrounding rock[J]. Journal of Traffic and Transportation Engineering, 2019, 19(3): 58-70. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JYGC201903007.htm
[13] 欧阳进武, 张贵金, 刘杰. 劈裂灌浆扩散机理研究[J]. 岩土工程学报, 2018, 40(7): 1328-1335. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201807026.htm OUYANG Jin-wu, ZHANG Gui-jin, LIU Jie. Diffusion mechanism of split grouting[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(7): 1328-1335. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201807026.htm
[14] ZHANG G, LIU J, LI Y, et al. A pasty clay-cement grouting material for soft and loose ground under groundwater conditions[J]. Advances in Cement Research, 2017, 29(2): 54-62.
[15] LIU J, LI Y, ZHANG G, et al. Effects of cementitious grout components on rheological properties[J]. Construction and Building Materials, 2019, 27(2): 1-13.
[16] LI Y, ZHANG G, JIANG H, et al. Performance assessment of a newly developed and highly stable sandy cementitious grout for karst aquifers in China[J]. Environmental Earth Sciences, 2020, 79(6): 1-15.
-
期刊类型引用(1)
1. 张辉，吴曙光，杨凯丞. 预留土支护基坑旋转破坏模式下的极限抗力上限解. 土木与环境工程学报(中英文). 2023(05): 116-124 . 
百度学术
其他类型引用(2)
计量
- 文章访问数: 282
- HTML全文浏览量: 20
- PDF下载量: 154
- 被引次数: 3
