单向与多层互剪搅拌桩性能模型试验对比研究

葛春巍; 刘钟; 兰伟; 杨宁晔; 文磊; 周洁

doi:10.11779/CJGE20230791

单向与多层互剪搅拌桩性能模型试验对比研究 English Version

葛春巍^{1, 2,},
刘钟^{1, 2, 3, ,},
兰伟^{1, 2},
杨宁晔^{1, 2},
文磊^{1, 2},
周洁^{4, 5}

1.
浙江坤德创新岩土工程有限公司, 浙江宁波 315100
2.
坤德智慧岩土技术研究院, 浙江宁波 315100
3.
中冶集团中冶建筑研究总院有限公司, 北京 100088
4.
同济大学地下建筑与工程系, 上海 200092
5.
同济大学岩土及地下工程教育部重点实验室, 上海 200092

基金项目:

宁波市高新区重大科技专项 2023CX050004

详细信息

作者简介:
葛春巍(1988—)，男，博士，工程师，主要从事地基处理装备研发及工艺工法研究。Email: gexxx063@alum.mit.edu

通讯作者:
刘钟, E-mail: zzliu8@163.com

中图分类号: TU443
计量
- 文章访问数: 486
- HTML全文浏览量: 41
- PDF下载量: 99
出版历程
- 收稿日期: 2023-08-15
- 网络出版日期: 2024-03-24
- 刊出日期: 2024-10-31

Model tests on performances of DSM and CS-DSM piles

GE Chunwei^{1, 2,},
LIU Zhong^{1, 2, 3, ,},
LAN Wei^{1, 2},
YANG Ningye^{1, 2},
WEN Lei^{1, 2},
ZHOU Jie^{4, 5}

1.
ZheJiang Kunde Innovate Geotechnical Engineering Co., Ltd., Ningbo 315100, China
2.
Kunde Research Institute of Intelligent Geotechnical Technology, Ningbo 315100, China
3.
Central Research Institute of Building and Construction, MCC Group, Beijing 100088, China
4.
Department of Geotechnical Engineering, College of Civil Engineering, Tongji University, Shanghai 200092, China
5.
Key Laboratory of Geotechnical and Underground Engineering (Tongji University), Ministry of Education, Shanghai 200092, China

摘要

摘要: 利用自研单向与多层互剪搅拌桩模型钻机，对传统DSM桩与新型CS-DSM桩的工艺因素及成桩质量进行了模型试验对比研究，探索了水泥掺量和单位桩长搅拌次数T对两类搅拌桩成桩均匀性与桩身强度的影响。试验结果证实CS-DSM桩的表观均匀性和桩身连续性明显好于DSM桩，在T值相同条件下，两者平均桩身强度比R_s为1.41~6.4，且R_s随着搅拌次数T增加而呈指数趋势降低，24组模型试验结果还揭示出DSM桩与CS-DSM桩在UCS-T之间的本质联系。提出的DSM桩与CS-DSM桩的T值和UCS值的计算方法，可以指导搅拌桩施工参数合理选取以实现桩身设计强度目标。通过DSM工法与CS-DSM工法的试验对比研究获得的两类桩的桩身质量差异性结果，能够为高质量的CS-DSM桩工艺控制原则和质量保障体系提供试验依据。
- 搅拌桩 /
- 多层互剪搅拌桩 /
- CS-DSM工法 /
- 搅拌次数 /
- 桩身均匀性和抗压强度
Abstract: Using the small-scale model rig, a comparative experimental study is conducted to investigate the process control and quality control in the deep soil mixing (DSM) method and contra-rotational shear deep soil mixing (CS-DSM) method. The research emphasis is placed on the effects of cement content, blade rotation number T on the uniformity and unconfined compression strength (UCS) of mixing piles. A notably improved uniformity and continuity of the CS-DSM piles over the DSM piles is confirmed by the test results. The strength ratio R_s ranges from 1.41 to 6.4, and it shows an exponential decrease with the increasing T. Furthermore, the results from 24 model piles provide insights into the fundamental relationship and distinctions between UCS and T. On this basis, the construction parameters can be optimized to ensure the target design strength for both types of piles. The model test results clearly demonstrate the differences in the construction quality between the DSM method and the CS-DSM method. The technical basis established in this study serves as a cornerstone for the process control and quality assurance in installing high-quality CS-DSM piles.
- DSM pile /
- CS-DSM pile /
- CS-DSM method /
- blade rotation number /
- pile uniformity and strength

HTML全文

0. 引言

随着中国西部大开发战略的持续推进，西部山岭地区隧道工程建设发展迅速。然而，隧道建设不可避免穿越活动断裂带，将会对隧道结构安全构成严重的威胁。例如，2008年汶川地震中，穿断层破碎带区域的隧道结构发生了严重的衬砌垮塌^[1]；2022年门源地震中，由于冷龙岭断裂带的剧烈走滑错动，使得大梁铁路隧道在断层错动面处发生超过2 m的水平位错，从而造成了极为严重的结构破坏^[2]。

目前国内外学者采用数值模拟手段针对穿活动断裂带隧道衬砌的地震响应特征开展了大量研究工作。Shahidi等^[3]基于二维有限元模型和拟静力方法，通过在断层区域一侧施加强制错动位移，研究了隧道结构沿纵向的受力分布规律；Anastasopoulos等^[4]基于有限元模拟对比分析了正断层错动和近断层强震动作用下的隧道结构非线性响应行为，发现断层错动作用对隧道的影响远大于近断层强地震动作用；甘星球等^[5]与Liu等^[6]同时考虑了断层错动与地震动的共同作用，研究表明活动断层错动作用是导致隧道结构破坏的主要因素。值得注意的是，目前研究主要针对给定断层错动量或地震动强度作用下的隧道结构响应。然而，实际活动断裂带错动作用具有显著的随机性分布特征，且历史地震观测表明，近断层地震动强度与活动断层错动量存在一定的相关性^[7-9]，可见目前研究仍无法反映活动断裂带运动的不确定性，因而无法在断层潜在风险下的工程决策提供概率性风险评估指导。

地震易损性分析是评估结构抗震性能的重要手段，其中易损性是指在一定强度地震作用下结构损伤状态等于或超过某极限状态的条件概率。目前基于性能的地震工程（performance-based earthquake engineering, PBEE）框架，并结合数值计算方法获得易损性曲线已成为该领域的主流趋势。Argyroudis等^[10-11]基于增量动力分析（incremental dynamic analysis, IDA）方法提出了水平横向地震动作用下存在腐蚀作用的浅埋地铁隧道结构易损性分析方法；Hu等^[12]利用等效线性化方法来考虑土层的动力非线性特征，进而提出了一种适用于软土隧道的地震易损性分析方法；Osmi等^[13]基于考虑岩石隧道非线性力学特征的数值分析模型，提出了含锚杆加固区岩石隧道的地震易损性分析方法。可见目前隧道工程抗震易损性主要关注地震动对隧道结构的影响，尚缺乏针对穿断裂带隧道结构抗错动性能的评估方法，因此难以满足跨断层隧道工程的实际建设需求。

本文旨在建立面向穿活动断裂带隧道结构抗错动性能评估的易损性分析方法。首先引入概率性断层位错灾害分析（PFDHA）方法，明确矩震级与活动断层位错量的概率关系；其次建立围岩-断层带-隧道结构耦合体系的三维弹塑性有限元模型，并进行拟静力推覆弹塑性分析；最后给出活动断裂带错动作用下隧道抗错动易损性曲线。基于该方法，通过参数化分析，进一步探究隧道选址、断层倾角等因素对穿活动断裂带隧道抗震易损性的影响规律。

1. 隧道抗错动易损性分析流程

本文建立的穿活动断裂带隧道抗错动易损性基本分析流程如图 1所示，具体如下：

图 1 穿活动断裂带隧道抗错动易损性分析流程

Figure 1. Analysis process of seismic fragility analysis method for evaluation of dislocation resistance of tunnels crossing active fault zones

下载: 全尺寸图片幻灯片

（1）建立断裂带平均位错（D_ave）与地震动强度（intensity measure, IM）的概率关系^[8]、位错相对值（D_norm）与破裂段相对位置（l/RL）的概率关系（D_norm为破裂带任意位置处位错量D与D_ave的比值^[14-15]），从而给出隧址区断层位错概率模型。

（2）建立围岩-断层-隧道相互作用体系的弹塑性分析数值模型，通过计算分析得出损伤指标（damage index，DI）与断层位错量D的量化关系。

（3）依据步骤（1）建立的隧址区断层位错概率公式，计算结构在不同性能水平下的失效概率，即可得到穿活动断裂带隧道抗错动易损性。

2. 断层位错超越概率计算方法与地震动强度参数的确定

2.1 断层位错超越概率计算流程

目前已有学者^[7]提出概率断层位错危险性分析方法（probabilistic fault displacement hazard analysis，PFDHA），旨在通过研究断层运动规律与概率关系以评估穿断层结构面临的危险性。

PFDHA方法一般简化断层几何结构模型如图 2所示。图中l为沿断层走向距左侧破裂端部的位置；RL为该破裂长度（rapture length）。

图 2 断层几何结构模型

Figure 2. Geometrical model for surface fault

下载: 全尺寸图片幻灯片

该方法认为断层位错发生在破裂段上，且破裂段平均位错量D_ave与地震动强度IM相关^[8]（目前PFDHA方法均以矩震级M_W作为IM）、破裂段上任意位置处的相对位错大小D_norm^[14-15]与破裂段相对位置l/RL相关。

可见，断层位错超越概率分析需要事先确定概率密度函数f_Dave(D_ave|IM)和f_Dnorm(D_norm|l/SRL)。为此，给出隧址区断层位错的超越概率计算流程如下：

（1）建立D_ave-IM概率关系以及D_norm-l/RL概率关系。

（2）计算断层破裂段任意位置处在不同地震动强度下的断层位错超越概率。

（3）最后考虑隧道穿越断层破裂段位置的不确定性，计算隧址区断层位错量的超越概率。

2.2 矩震级与断层平均位错关系

地震动强度参数（IM）是地震易损性分析的关键指标之一。在传统隧道抗震易损性分析中，通常将峰值地表加速度（PGA）作为IM指标^[16]。虽然目前研究仍无法给出活动断裂带错动作用下PGA与断层位错量之间的量化表征关系，但是关于断层参数与矩震级（M_W）关系的研究相对较多^{[8, 17-19]}，矩震级概念最早由Kanamori^[18]提出，定义如下式所示^[20]：

${M_0} = \mu \cdot {\text{RL}} \cdot {\text{RW}} \cdot {D_{{\text{ave}}}} \text{，}$

(1)

${M_{\text{W}}} = \frac{2}{3}\log \left[ {\frac{{{M_0}}}{{{M_0}(0)}}} \right] 。$

(2)

式中：M₀为地震矩；RL为该地震矩下断层的破裂长度；RW为该地震矩下断层的破裂宽度；D_ave为该地震矩下断层的平均位错量；μ为断层的剪切模量；M₀(0)为矩震级为0的地震矩，其值为10^9.1 N·m。

矩震级作为地震矩的一种映射变换，能够与断裂带破裂过程的物理表征直接关联，且能够合理反映地震作用下断层错动的物理特性和本征表达，故本文选用M_W作为IM来表征地震动强度。

由式（1），（2）可知，断层平均位错量D_ave与矩震级M_W服从对数关系，可假设给定M_W下的D_ave服从对数正态分布。若lg(D_ave)的均值和标准差分别为 $\lg ({\mu _{{\text{AD}}}})$ 和σ_AD，则随机变量D_ave|M_W的概率密度函数可表示为

$\lg ({\mu _{{\text{AD}}}}) = a \cdot {M_{\text{W}}} - b \text{，}$

(3)

$\begin{aligned} & f_{D_{\mathrm{avc}}}\left(D_{\mathrm{ave}} \mid M_{\mathrm{W}}\right) \\ & \quad=\frac{1}{\ln 10 \sqrt{2 \pi} D_{\mathrm{ave}} \sigma_{\mathrm{AD}}} \cdot \exp \left[-\frac{\left(\lg D_{\mathrm{ave}}-\lg \mu_{\mathrm{AD}}\right)^2}{2 \sigma_{\mathrm{AD}}^2}\right] \end{aligned} 。$

(4)

式中：a和b为经验系数，可根据不同的数据集拟合求得。

2.3 位错相对值与破裂带相对位置关系

目前普遍采用的断层位错分布模型通过基于观测的经验分析来描绘D_norm与l/RL的关系。记D_norm的均值为μ_ND，通过拟合可以得到μ_ND-l/RL函数关系^[14]。

根据拟合函数的差异性，目前已有针对 ${\mu _{{\text{ND}}}}$ -l/RL曲线形状的不同分布模型，例如对称椭圆分布^[15]、对称双折线分布^[7]、非对称三角分布^[19]等，可以统一表示为

${\mu _{{\text{ND}}}} = g(l/{\text{RL}}) 。$

(5)

式中：g(·)为拟合函数形式（椭圆形、三角形等）。

通常假设随机变量D_norm|l/RL服从对数正态分布，以拟合结果残差标准差作为lg(D_norm)的标准差 ${\sigma _{{\text{ND}}}}$ ，则其概率密度函数可以表示为

$\begin{aligned} & f_{D_{\text {norm }}}\left(D_{\text {norm }} \mid l / \mathrm{RL}\right) \\ & \quad=\frac{1}{\ln 10 \sqrt{2 \pi} D_{\text {norm }} \sigma_{\mathrm{ND}}} \cdot \exp \left[-\frac{\left(\lg D_{\text {norm }}-\lg \mu_{\mathrm{ND}}\right)^2}{2 \sigma_{\mathrm{ND}}^2}\right] 。 \end{aligned}$

(6)

2.4 隧址区断层位错超越概率计算公式

若已知 ${f_{{D_{{\text{ave}}}}}}$ (D_ave|M_W)和 ${f_{{D_{{\text{norm}}}}}}$ (D_norm|l/RL)，断层位错概率关系可采用联合条件概率密度函数表示为

$\begin{aligned} & f_{\mathrm{D}}\left(D \mid M_{\mathrm{w}}, l / \mathrm{RL}\right)=f_{\mathrm{D}}\left(D_{\text {nom }} \cdot D_{\text {ave }} \mid M_{\mathrm{w}}, l / \mathrm{RL}\right) \\ & \quad=\int_0^{+\infty} \frac{1}{D_{\text {ave }}} f_{D_{\text {ave }}}\left(D_{\text {ave }} \mid M_{\mathrm{w}}\right) f_{D_{\text {ave }}}\left(\left.\frac{D}{D_{\text {ave }}} \right\rvert\, l / \mathrm{RL}\right) \mathrm{d} D_{\text {ave }} 。 \end{aligned}$

(7)

则在某地震动作用下穿破裂带任意位置处隧道结构的断层位错超越概率为

$P(D > {D_{\text{o}}}|{M_{\text{W}}},l/{\text{RL}}) = \int\limits_{{D_0}}^{ + \infty } {{f_{\text{D}}}(D|{M_{\text{W}}},l/{\text{RL}})} 。$

(8)

另外，若隧道选址位置仍存在不确定性，可假设隧道选址在破裂段服从均匀分布，则隧址区在某地震动作用下超越一定位错的概率可表示为

$\begin{aligned} P(D & \left.>D_0 \mid M_{\mathrm{w}}\right) \\ & =\int_0^1 P\left(D>D_0 \mid M_{\mathrm{w}}, l / \mathrm{RL}\right) \mathrm{d}(l / \mathrm{RL}) \quad 。 \end{aligned}$

(9)

3. 数值模型与损伤指标确定

3.1 数值模型建立

本文重点关注活动断层的错动作用，忽略断层产生的地震动及惯性力效应，采用拟静力方法分析。

本节模型以走滑断层为例，图 3为穿活动断裂带隧道的围岩-断层-结构三维模型示意图。模型宽度取隧道直径的5倍以上，可忽略侧向边界的影响^[21]。图中阴影区域为活动断裂带分布范围，宽度为W。假设断层错动面位于断裂带中心断面，错动面倾角为β，且错动面处的力学行为服从莫尔-库仑摩擦接触定律。基于拟静力分析方法，活动断层带两侧地盘约束如下：一侧地盘视为固定；另一侧通过施加强制位移来模拟活动断层的错动作用。

图 3 模型示意图

Figure 3. Schematic diagram of model

下载: 全尺寸图片幻灯片

隧道结构抗错动基本分析流程如下：

步骤1：围岩初始地应力平衡；

步骤2：采用收敛约束法^[22]进行隧道开挖模拟；

步骤3：采用静力形式逐级输入断层位错。

3.2 结构损伤指标

结构损伤指标（DI）是易损性分析所需的关键参数之一，用于定量描述隧道结构的损伤状态。

基于有限元模型，考虑到隧道结构在断层错动作用下的三维非线性响应特征，选用混凝土损伤塑性（CDP）模型对衬砌结构损伤状态进行表征，该模型能够较好地表征隧道衬砌的损伤状态以及复杂受力状态下混凝土材料的开裂及压碎行为^[23]。CDP中混凝土的受压和受拉损伤采用归一化指标d_c和d_t来表征。本文针对隧道横断面定义某断面的损伤指标。

首先选取结构每个单元质心位置处的损伤因子（压碎损伤指标d_ci与开裂损伤指标d_ti最大值）来表征该单元的损伤大小。那么隧道不同横断面位置处的结构损伤可由该断面所有单元损伤对单元体积加权平均得到，即损伤体积比，如式（10）。该指标可定量表征该断面结构的整体损伤程度。

${\text{DI}} = \frac{{\Sigma _{i = 1}^n{d_i}{V_i}}}{{\Sigma _{i = 1}^n{V_i}}},{d_i} = {\text{Max}}\left\{ {{d_{{\text{c}}i}},{d_{{\text{t}}i}}} \right\} 。$

(10)

式中：d_i为第i个单元质心位置处的损伤因子大小，V_i为第i个单元的体积。

以隧道结构的弹塑性损伤发展状况为安全状态划分依据，将各衬砌截面按照损伤指标DI大小划分为安全、轻微损伤、中等损伤以及严重损伤4个损伤状态（damage status，ds），如表 1所示^[24]。

表 1 损伤状态

Table 1. Damage states

${\text{d}}{{\text{s}}_i}$	损伤状态	DI	描述
${\text{d}}{{\text{s}}_1}$	安全	0≤DI＜0.1	主要发生弹性变形
${\text{d}}{{\text{s}}_2}$	轻微损伤	0.1≤DI＜0.4	结构出现较大的塑性变形，但强度和功能性无实质性减弱
${\text{d}}{{\text{s}}_3}$	中等损坏	0.4≤DI＜0.7	局部或整体塑性变形超出功能性限值，结构仍保持完整性
${\text{d}}{{\text{s}}_4}$	严重损伤	DI≥0.7	严重损伤，位移和塑性变形发散，结构发生失稳、压溃、断裂、失去完整性

下载: 导出CSV

| 显示表格

4. 分析与讨论

4.1 分析模型及参数

本节旨在通过算例分析说明本文方法的具体实施流程及可行性。假设某隧道工程穿越走滑活动断层，埋深95 m，其轴线与断层走向正交。该断层倾角76°，断层破碎带宽度36 m，断层错动面位于破碎带中心。隧道截面形式为圆形，内径和外径分别为8.6，10 m。隧道和围岩的力学参数如表 2所示。

表 2 围岩力学参数

Table 2. Mechanical parameters of surrounding rock

名称	密度 ρ/(g·cm^-2)	杨氏模量 E/MPa	泊松比 $\nu$	摩擦角φ/(°)	屈服应力 ${\sigma _{\text{y}}}$ /MPa
破碎带	2.0	300	0.35	30	5
岩体	2.5	2500	0.30	45	30
C30 混凝土	2.4	24500	0.20	—	—

下载: 导出CSV

| 显示表格

依据本文建立的易损性计算方法对该穿断层隧道结构的抗错动性能进行分析，建立有限元模型如图 4所示，模型整体尺寸为400 m×200 m×200 m。隧道与围岩均采用六面体实体单元模拟，两者间设置绑定接触。以摩擦接触模拟断层错动面，摩擦系数取0.7。模拟断层错动时，位错施加在上盘围岩两侧。岩体塑性本构选用Drucker-Prager屈服准则，混凝土选用CDP塑性损伤本构模型，本文C30混凝土对应的CDP本构曲线标定如图 5所示。

图 4 有限元分析模型

Figure 4. Schematic diagram of FEM model

下载: 全尺寸图片幻灯片

图 5 C30混凝土应力应变及损伤关系

Figure 5. Relationship among stress, plastic strain and damage of C30

下载: 全尺寸图片幻灯片

4.2 断层位错的概率密度函数

对于D_ave-IM概率关系，基于Anderson等^[9]共55次地震数据进行拟合，得到走滑断层平均位错量与矩震级关系，如所示，其中标准差 ${\sigma _{{\text{AD}}}}$ = 0.251，平均位错量均值 ${\mu _{{\text{AD}}}}$ 关于M_W的拟合函数表达为

图 6 断层破裂带平均位错关于矩震级拟合关系

Figure 6. Regression diagram of rupture average dislocation with respect to moment magnitude

下载: 全尺寸图片幻灯片

$\lg ({\mu _{{\text{AD}}}}) = 0.598{M_{\text{W}}} - 4.055 。$

(11)

关于D_norm-l/RL概率关系，本文基于Wesnousky^[25]共9次地震的位错分布数据，采用椭圆分布进行拟合，得到如所示的走滑断层位错相对值与破裂带相对位置的关系，其中标准差 ${\sigma _{{\text{ND}}}}$ = 0.379，位错相对值均值 ${\mu _{{\text{ND}}}}$ 关于破裂带相对位置l/RL的拟合函数表示为

图 7 断层位错相对值与破裂带相对位置拟合关系

Figure 7. Fitting relationship between relative value of dislocation and relative position of rupture

下载: 全尺寸图片幻灯片

$\lg ({\mu _{{\text{ND}}}}) = 1.699\sqrt {1 - 4{{(l/{\text{RL}} - 0.5)}^2}} - 1.742 。$

(12)

将式（10），（11）代入式（4），（6），即可得到 ${f_{{D_{{\text{ave}}}}}}$ (D_ave|IM)与 ${f_{{D_{{\text{norm}}}}}}$ (D_norm|l/RL)，再根据数值模拟确定各损伤状态的临界位错量，代入式（8）即可求得隧址区在某地震动水准下超越该临界位错的概率，即隧道结构的抗活断层错动易损性。

4.3 易损性分析

根据3.2节定义的损伤指标，重点针对隧道结构3个关键断面进行易损性分析，即错动面与隧道相交断面、上下盘-断层破碎带边界面与隧道相交断面，如图 8所示。

图 8 关键断面位置示意图

Figure 8. Section diagram of three key tunnel locations

下载: 全尺寸图片幻灯片

以错动面处隧道断面损伤情况为例，图 9为隧道衬砌损伤指标DI随断层错动量D的演化关系曲线。根据表 1定义的隧道结构损伤状态分级，从图中可以看出：该断面在断层错动量达到0.2 m时进入轻微损伤状态；断层错动量达到0.32 m时，结构进入中等损坏状态；断层错动量达到0.5 m时达到严重损伤状态。

图 9 断层错动面处衬砌段损伤指标与位错量关系

Figure 9. Relationship between damage index of lining section at dislocation plane and slip-displacement

下载: 全尺寸图片幻灯片

图 10展示了上述4个阶段下该隧道断面的损伤分布情况。从损伤演化特征可以看出，该断面结构损伤首先在拱顶和仰拱处开始发展，随着断层错动量的增加，结构损伤逐渐扩展至整个截面，最终导致该衬砌截面的整体失效。将断面各点损伤以体积进行加权平均可得该断面的损伤指标。

图 10 不同错动量下断层错动面处衬砌段损伤开展分布情况

Figure 10. Damage development and distribution of lining section at dislocation surface

下载: 全尺寸图片幻灯片

选取穿越断层破碎带的上下盘边界面以及断层错动面的3处关键隧道断面，各断面的易损性曲线如图 11所示。图中曲线分别为隧道3个关键断面在不同矩震级下达到或超越规定损伤状态的概率。

图 11 关键断面位置处隧道易损性曲线

Figure 11. Fragility curves of lining section at key section

下载: 全尺寸图片幻灯片

由图 11可知，在同一矩震级下，断层错动面处隧道结构的损伤概率最高，上盘与断层交界处隧道断面次之，下盘与断层交界处隧道断面最低。例如，对于中度损伤状态，在矩震级为7的情况下，断层错动面处隧道的易损性概率为0.67，上盘与断层交界处隧道失效概率为0.51，下盘与断层交界处仅为0.47。这是由于断层错动面以强制剪切作用为主，变形最大且损伤最严重；而上下盘与断层交界处的地层存在差异，会发生差异变形，但这种变形程度不及在错动面直接发生的剪切作用，因此断层错动面处对隧道更为不利。对比上下盘与断层交界处隧道断面的易损性曲线可以看出，在同一震级下，上盘位置处隧道达到同等级损伤状态的概率高于下盘位置处隧道，即上盘效应^[26]。

此外，还可以发现在震级为4~5级和8~9级的范围内，易损性曲线变化较为平缓，特别是在4~5级震级范围内，隧道的损伤概率相对较低且变化较小，可以不需要采取额外的抗震措施。然而，在5~8级震级范围内，易损性曲线变化非常显著，表明对震级变化的敏感性较强，即使是小幅度的震级变化，也会导致隧道损伤状态出现显著的波动。因此，在工程设计阶段，需要精确评估地震动水平并提供足够的抗震性能，以应对这种较大的随机性和波动性。

4.4 参数化分析

（1）断层倾角对隧道抗错动易损性的影响

在4.1节分析模型及参数的基础上，为进一步研究断层倾角对隧道抗错动易损性的影响，考虑到实际断层倾角的分布情况，模型中断层倾角分别取为63°，76°，90°，计算模型及参数同上节。

以中度损伤状态为例，图 12给出了不同断层倾角下3个关键断面处隧道结构的抗错动易损性曲线。由图可知，断层倾角越大，隧道易损性相对越高。如矩震级为7时，断层错动面处隧道的易损性概率分别为0.62（断层倾角为63°）、0.67（断层倾角为76°）和0.71（断层倾角为90°）。这是由于随着断层倾角的增加，断层与隧道衬砌相交范围减小，隧道衬砌段遭受的断层位错更为集中，损伤也更加严重。因此，应对穿越倾角较大走滑断层的隧道加以重视，应在选线时尽量避开或者采用合理的抗错断措施。

图 12 穿不同倾角断层隧道结构中等损伤状态的易损性曲线

Figure 12. Moderate damage fragility curves of tunnel crossing different dip angles fault at key section

下载: 全尺寸图片幻灯片

此外，通常认为上下盘在断层处存在的地应力突变是导致上盘效应的关键因素^[26]，而且断层倾角越小，地应力突变也就越大，上盘效应也越为显著。如图 12所示，当倾角为90°时，上下盘与断层交界处隧道抗错动易损性曲线基本重合，但对于倾角为76°及63°的工况而言，下盘与断层交界处隧道的抗错动易损性低于上盘，且断层倾角越小此现象越为明显。

（2）穿断层隧道选址对结构抗错动易损性的影响

在4.1节分析模型及参数的基础上，为进一步研究穿断层隧道选址对结构抗错动易损性的影响机制，假设隧道穿破裂带位置l/RL分别取为0，0.1，0.2，0.3，0.4，0.5，即涵盖了破裂带端部至中部不同位置，计算模型及参数同上节。

以断层错动面处隧道结构的中度损伤状态为例，图 13给出了不同选址处的抗错动易损性曲线。由图 13可知，选址越靠近破裂带端部，隧道结构的易损性相对越低。如矩震级为7时，隧道选址l/RL为0，0.1，0.2相对l/RL=0.5时的易损性可分别降低99.3%，54.6%，22.4%。然而当选址位置进一步靠近破裂带中心，如l/RL = 0.4，易损性变化幅度明显减小，仅降低1.6%。若要求该断面在矩震级7级时达到中等破坏的概率低于0.4，则隧道选址需位于距破裂带端部1/10的总破裂带长度的区间内。分析表明即使无法完全避免隧道穿越破裂带，如选择尽量靠近破裂带端部位置穿越则可显著降低隧道结构的受损概率。因此，建议隧道选址应充分考虑断层破裂带与隧道线路的位置关系，最大程度提高隧道结构的抗错动安全性。

图 13 穿不同位置破裂带隧道中等损伤状态易损性曲线

Figure 13. Fragility curves of moderate damage state in tunnels crossing different rupture positions

下载: 全尺寸图片幻灯片

5. 结论

本文提出了面向穿活动断裂带隧道结构的抗错动易损性分析方法，首先建立了PFDHA概率分析模型，并基于该模型选用矩震级作为地震动强度指标，依据混凝土损伤塑性模型建立了隧道结构损伤指标，进而通过矩震级大小与活动断层位错量的概率统计回归关系，结合基于三维弹塑性有限元模型的拟静力推覆分析方法，建立了穿断层隧道抗错动易损性程度的表征方法，从而为穿活动断裂带隧道抗震风险评估及结构抗震设计提供了有效的分析手段。基于该方法对某穿走滑断裂带隧道进行应用分析，并研究了断层倾角对隧道抗错动易损性的影响规律，主要得到以下3点结论。

（1）断层错动面处的衬砌段易损性最高，断层破碎带与上下盘交界面处的衬砌段易损性次之。相比而言，穿走滑断裂带隧道上盘易损性高于下盘，即表现出明显的上盘效应，而且断层倾角越小，错动面处隧道易损性越小，上盘效应也更为显著。

（2）在矩震级为4~5级和8~9级范围内，隧道易损性曲线变化较为平缓，而在5~8级矩震级范围内，隧道的易损性变化显著，对震级变化较为敏感。因此在工程设计阶段，需要准确评估地震动水平以应对较大概率的波动变化。

（3）隧道选址对结构抗错动易损性影响显著，如尽量靠近破裂带端部位置穿越可显著降低隧道结构的受损概率。建议隧道工程线路规划时应充分考虑穿越破裂带的不同位置关系，以最大程度提高隧道结构的抗错动安全性。

需要说明的是，本文建立的隧道易损性分析方法适用于任何活动断层类型，文中采用的概率密度函数和有限元模型边界条件仅适用于走滑断层，但对于其他类型断层，仅需调整有限元模型并依据统计关系修改概率密度函数。针对具体工程，可依据工程场地的历史震害数据建立概率分析模型，进而应用本方法为穿活动断裂带隧道结构抗震设计与抗错动性能评估提供科学依据。

图 1 DSM桩、DDM桩、CS-DSM桩的施工钻头对比

Figure 1. Comparison between three different drill bits

下载: 全尺寸图片幻灯片

图 2 SXJ-110-D型多层互剪搅拌桩施工钻机装备

Figure 2. CS-DSM drilling rig model No. SXJ-110-D

下载: 全尺寸图片幻灯片

图 3 智能控制ACS系统框图

Figure 3. Layout of autonomous control system

下载: 全尺寸图片幻灯片

图 4 模型箱、模型钻机与ACS系统

Figure 4. Model case, drilling rig and ACS

下载: 全尺寸图片幻灯片

图 5 模型钻机动力头

Figure 5. Rotary power heads for drilling rig

下载: 全尺寸图片幻灯片

图 6 搅拌桩钻具

Figure 6. Illustration of drill bits for DSM and CS-DSM piles

下载: 全尺寸图片幻灯片

图 7 模型桩位布置图

Figure 7. Position of piles in model foundation

下载: 全尺寸图片幻灯片

图 8 色砂试验的表观均匀性对比照片(T = 400)

Figure 8. Results of uniformity tests for DSM and CS-DSM piles

下载: 全尺寸图片幻灯片

图 9 两类搅拌桩的试块均匀性对比照片(T = 640~650)

Figure 9. Pictures of soilcrete cubes for DSM and CS-DSM piles

下载: 全尺寸图片幻灯片

图 10 水泥掺量与桩身强度的关系曲线

Figure 10. Relationship between UCS and cement content

下载: 全尺寸图片幻灯片

图 11 DSM桩的桩身强度与搅拌次数的关系曲线

Figure 11. Relationship between UCS and T for DSM piles

下载: 全尺寸图片幻灯片

图 12 CS-DSM桩的桩身强度与搅拌次数的相关关系

Figure 12. Relationship between UCS and T for CS-DSM piles

下载: 全尺寸图片幻灯片

图 13 DSM桩与CS-DSM桩的桩身强度对比

Figure 13. Comparison of UCS between DSM and CS-DSM piles

下载: 全尺寸图片幻灯片

图 14 强度比R_s与单位桩长搅拌次数T的相关性

Figure 14. Strength ratios between DSM and CS-DSM piles

下载: 全尺寸图片幻灯片

图 15 两类桩身强度变异系数与搅拌次数的相关关系

Figure 15. Relationship between COV and T

下载: 全尺寸图片幻灯片

表 1 模拟地基的物理力学指标

Table 1 Physico-mechanical parameters of model soil

密度/(g·cm^-3)	含水率/%	液限/%	塑限/%	黏聚力/kPa	内摩擦角/(°)	塑性指数
1.84	30	34.1	20.5	8.4	13.7	13.6

下载: 导出CSV

表 2 桩身强度与水泥掺量关系的试验方案

Table 2 Test plans for studying cement content and UCS

桩号	水泥掺量/%	阶段	转速/ rpm		升降速度/(m·min^-1)	浆液流量/(L·min^-1)	搅拌次数/(rev·m^-1)
桩号	水泥掺量/%	阶段	内杆	外杆	升降速度/(m·min^-1)	浆液流量/(L·min^-1)	搅拌次数/(rev·m^-1)
1	8	下/上	单向	24	0.28/0.28	0.33/ —	686
2	13	下/上	单向	24	0.28/0.28	0.54/ —	686
3	18	下/上	单向	24	0.28/0.28	0.75/ —	686
CS1	8	下/上	10	6	0.28/0.28	0.33/ —	686
CS2	13	下/上	10	6	0.28/0.28	0.54/ —	686
CS3	18	下/上	10	6	0.28/0.28	0.75/ —	686
注：平均养护温度为17℃。

下载: 导出CSV

表 3 搅拌次数对比试验方案：DSM桩试验参数

Table 3 Test plans for studying blade rotation number for DSM piles

桩号	阶段	转速/ rpm		升降速度/（m·min^-1）	浆液流量/（L·min^-1）	搅拌次数/（rev·m^-1）
桩号	阶段	下钻	提钻	升降速度/（m·min^-1）	浆液流量/（L·min^-1）	搅拌次数/（rev·m^-1）
4	下/上	9	14	0.25/0.36	0.48/ —	300
5	下/上	11	20	0.25/0.36	0.48/ —	398
6	下/上	13	22	0.25/0.36	0.48/ —	452
7	下/上	15	23	0.25/0.36	0.48/ —	496
8	下/上	17	25	0.25/0.36	0.48/ —	550
9	下/上	18	28	0.25/0.36	0.48/ —	599
10	下/上	21	28	0.25/0.36	0.48/ —	647
11	下/上	22	31	0.25/0.36	0.48/ —	696
12	下/上	22	32	0.25/0.32	0.48/ —	752
13	下/上	23	33	0.25/0.30	0.48/ —	808
14	下/上	25	33	0.25/0.30	0.48/ —	853
15	下/上	25	35	0.25/0.28	0.48/ —	900
16	下/上	25	36	0.25/0.24	0.48/ —	1000
注：平均养护温度为32 ℃。

下载: 导出CSV

表 4 搅拌次数对比试验方案：CS-DSM桩试验参数

Table 4 Test plans for studying blade rotation number for CS-DSM piles

桩号	阶段	转速/ rpm		升降速度/(m·min^-1)	浆液流量/(L·min^-1)	搅拌次数/(rev·m^-1)
桩号	阶段	内杆	外杆	升降速度/(m·min^-1)	浆液流量/(L·min^-1)	搅拌次数/(rev·m^-1)
CS4	下/上	6	4	0.28/0.28	0.54/ —	429
CS5	下/上	9	6	0.28/0.28	0.54/ —	643
CS6	下/上	11	7	0.28/0.28	0.54/ —	771
CS7	下/上	13	9	0.28/0.28	0.54/ —	943
CS8	下/上	19	12	0.28/0.28	0.54/ —	1329
注：平均养护温度为32℃。

下载: 导出CSV

参考文献(26)

[1]	COLLE E R. Mixed in place pile: USA, US3270511[P]. 1966-09-06.
[2]	PORBAHA A. State of the art in deep mixing technology: part Ⅰ Basic concepts and overview[J]. Proceedings of the Institution of Civil Engineers-Ground Improvement, 1998, 2(2): 81-92. doi: 10.1680/gi.1998.020204
[3]	DENIES N, HUYBRECHTS N. Deep Mixing Method[M]// Ground Improvement Case Histories. Amsterdam: Elsevier, 2015: 311-350.
[4]	FILZ G, BRUCE D. Innovation and collaboration in deep mixing[C]// Grouting 2017. Honolulu, 2017: 336-353.
[5]	PORBAHA A, SHIBUYA S, KISHIDA T. State of the art in deep mixing technology. Part Ⅲ: geomaterial characterization[J]. Proceedings of the Institution of Civil Engineers - Ground Improvement, 2000, 4(3): 91-110. doi: 10.1680/grim.2000.4.3.91
[6]	KITAZUME M. Quality Control and Assurance of the Deep Mixing Method[M]. London: CRC Press, 2021.
[7]	O'ROURKE T D, MCGINN A J. Case history of deep mixing soil stabilization for Boston central artery[C]//Geotechnical Engineering for Transportation Projects. Los Angeles, 2004: 77-136.
[8]	SEIJI M, SHUBUYA S, TETSUJI M, et al. Soil improvement device: Japan, JP08-199556[P]. 1996-08-06.
[9]	木付拓磨, 澤口宏, 今井正, など. 大口径大深度深層混合処理工法の適用におけるリアルタイム管理システムの導入[J]. 日本材料科学学会雑誌, 2018, 67(1): 93-98. KIZUKI T, SAWAGUCHI H, IMAI T, et al. An introduction of real-time management system for applying deep-mixing method with large diameter and improvement depth (DCS method)[J]. Journal of the Society of Materials Science, Japan, 2018, 67(1): 93-98. (in Japanese)
[10]	鈴木孝, 齋藤邦夫, 原満生, など. 複合相対攪拌翼を用いた深層混合処理工法の改良原理と適用事例[J]. 日本材料科学学会雑誌, 2010, 59(1): 32-37. SUZUKI K, SAITOH K, HARA M, et al. Mixing mechanism and case study of deep stabilization method using contra-rotational mixing head[J]. Journal of the Society of Materials Science, 2010, 59(1): 32-37. (in Japanese)
[11]	MORI K, UKAJI N, MIYAKAWA M. Invited Lecture: recent trends in the development of deep mixing methods in Japan[C]// Geotechnics for Sustainable Infrastructure Development. Singapore, 2020: 541-554.
[12]	島野嵐. 大口径相対攪拌工法の概要と施工事例: KS-S·MIX工法(特集基礎工, 地盤改良)建設機械施工一般社団法人[J]. 日本建設機械施工協会誌, 2017, 69(7): 59-63. SHIMANO A. Case study: large diameter contra-rotational soil mixing pils using KS-S MIX method[J]. JCMA, 2017, 69(7): 59-63. (in Japanese)
[13]	刘松玉, 宫能和, 冯锦林, 等. 双向水泥土搅拌桩机: CN1632232A[P]. 2005-06-29. LIU Songyu, GONG Nenghe, FEN Jinling, et al. Bi-directional cement soil mixing pile machine: CN1632232A[P]. 2005-06-29. (in Chinese)
[14]	刘松玉, 易耀林, 朱志铎. 双向搅拌桩加固高速公路软土地基现场对比试验研究[J]. 岩石力学与工程学报, 2008, 27(11): 2272-2280. doi: 10.3321/j.issn:1000-6915.2008.11.014 LIU Songyu, YI Yaolin, ZHU Zhiduo. Comparison tests on field bidirectional mixing column for soft ground improvement in expressway[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(11): 2272-2280. (in Chinese) doi: 10.3321/j.issn:1000-6915.2008.11.014
[15]	刘松玉. 新型搅拌桩复合地基理论与技术[M]. 南京: 东南大学出版社, 2014. LIU Songyu. Innovative Deep Mixing Method-Theory and Technology[M]. Nanjing: Southeast University Press, 2014. (in Chinese)
[16]	刘钟, 陈天雄, 杨宁晔, 等. 一种具有双向旋搅机构的智能钻机装备: CN218405399U[P]. 2023-01-31. LIU Zhong, CHEN Tianxiong, YANG Ningye, et al. Intelligent drilling machine equipment with bidirectional rotary stirring mechanism: CN218405399U[P]. 2023-01-31. (in Chinese)
[17]	刘钟, 李国民, 王占丑, 等. 一种搅拌桩机用同心三管三通道钻杆结构和组合注浆钻具: CN114319334A[P]. 2022-04-12. LIU Zhong, LI Guomin, WANG Zhanchou, et al. Concentric three-pipe three-channel drill rod structure for mixing pile machine and combined grouting drilling tool: CN114319334A[P]. 2022-04-12. (in Chinese)
[18]	陈湘生, 洪成雨, 苏栋. 智能岩土工程初探[J]. 岩土工程学报, 2022, 44(12): 2151-2159. doi: 10.11779/CJGE202212001 CHEN Xiangsheng, HONG Chengyu, SU Dong. Intelligent geotechnical engineering[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(12): 2151-2159. (in Chinese) doi: 10.11779/CJGE202212001
[19]	ANTSAKLIS P J, PASSINO K M, WANG S J. An introduction to autonomous control systems[J]. IEEE Control Systems Magazine, 1991, 11(4): 5-13. doi: 10.1109/37.88585
[20]	NAKAO K, INAZUMI S, TAKAUE T, et al. Visual evaluation of relative deep mixing method type of ground-improvement method[J]. Results in Engineering, 2021, 10: 100233. doi: 10.1016/j.rineng.2021.100233
[21]	NAKAO K, INAZUMI S, TAKAUE T, et al. Evaluation of discharging surplus soils for relative stirred deep mixing methods by MPS-CAE analysis[J]. Sustainability, 2021, 14(1): 58. doi: 10.3390/su14010058
[22]	ANDROMALOS K B, BAHNER E W. The application of various deep mixing methods for excavation support systems[C]// Grouting and Ground Treatment, New Orleans, 2003: 515-526.
[23]	SHEN S L, MIURA N, KOGA H. Interaction mechanism between deep mixing column and surrounding clay during installation[J]. Canadian Geotechnical Journal, 2003, 40(2): 293-307.
[24]	深层搅拌法地基处理技术规范: DL/T5425—2018[S]. 北京: 中国电力出版社, 2018. Technical Specification for Deep Mixing Ground Treatment: DL/T5425—2018[S]. Beijing: China Electric Power Press, 2018. (in Chinese)
[25]	水泥土配合比设计规程: JGJ/T 233—2011[S]. 北京: 中国建筑工业出版社, 2011. Specification for Mix Proportion Design of Cement Soil: JGJ/T 233—2011[S]. Beijing: China Architecture & Building Press, 2011. (in Chinese)
[26]	NAVIN M P, FILZ G M. Statistical analysis of strength data from ground improved with DMM columns[C]//Deep Mixing '05: Proc Int Conf on Deep Mixing Best Practice and Recent Advances, Stockholm, 2005.