0.   引言

中国现阶段正处于水下盾构隧道建设的高峰期，水下隧道的建设总体呈现出大断面、大埋深、高水压和地质条件日益复杂的趋势。不良地质段的沉降控制是海底隧道施工的核心技术之一，它不仅对施工安全至关重要，而且对隧道运营阶段的结构受力及长期安全性也是极其重要的^[1]。断层破碎带是海底盾构隧道建设过程中常见的不良地质段，在断层破碎带地层界面附近，地层特性和隧道上覆荷载均发生显著变化，极易造成隧道的不均匀沉降，从而导致环缝张开、错台、螺栓应力增大乃至接缝渗漏水等现象，严重影响盾构隧道的结构稳定和运营安全^[2]。

目前常用的隧道上覆荷载计算方法有全土柱法、太沙基理论^[3]、普氏理论^[4]和规范公式^[5]等。黎春林^[6]在Terzaghi^[7]基础上推导了一种可以考虑地层损失和管片刚度的松动土压力计算公式。王浩等^[8]基于普氏平衡拱理论，提出适用于3洞小净距隧道的围岩压力计算方法。颉永斌等^[9]用筒仓理论结合地层应力分布特征来求解断层破碎带内隧道的上覆荷载。由于断层破碎带形态各异，填充物复杂，目前断层破碎带地层隧道上覆荷载仍难以全面准确描述。

不同学者^[10-12]将既有盾构隧道简化为无限长弹性直梁或Euler-Bernoulli梁放置于Winkler地基模型或Pasternak地基模型上，研究既有隧道的力学响应。然而Euler-Bernoulli梁忽略了实际工程常发生的剪切变形。为更合理地描述隧道的实际变形模式，Wu等^[13]针对隧道的剪切变形，基于Liao等^[14]接头有一定影响范围思路同样引入剪切刚度影响系数推导隧道管片等效抗剪刚度，将隧道视为Timoshenko梁，从而考虑隧道纵向剪切效应。Zhang等^[15]在Wu等^[13]基础上将既有隧道简化为能够考虑地层的压缩及剪切效应的kerr三参数地基上Timoshenko梁，并依据高斯曲线推导出新隧道对既有隧道的纵向响应的解析解。目前Timoshenko梁在多段非均布复杂上覆荷载作用下难以直接求解，通常以有限差分法进行数值计算。

为此，本文将隧道上覆荷载概括为4种特征形态断层破碎带松动土压力的线性组合，提出了归一化上覆荷载理论公式，并将复杂上覆荷载在5%的误差内简化为三次泰勒多项式，结合Timoshenko梁理论推导建立了断层破碎带海底盾构隧道管片纵向沉降及内力系列解析解，并通过现场实测结果验证了本文理论模型和解析解的正确性。

1.   断层破碎带隧道纵向上覆荷载分析

1.1   断层破碎带形态分析

断层破碎带是海底隧道建设中经常遇到的不良地质条件，由于断层破碎带内岩体构造发育，渗透性较大，使得海底盾构隧道的上覆荷载比较复杂。如图 1所示，将隧道沿纵向分为4段：ab段为正常段，bc段为相交段，cd段为延伸段，de段为正常段。不同的相交段和延伸段组合可模拟不同形态的断层破碎带。

图  1  盾构隧道分段示意图

Figure  1.  Sectional diagram of shield tunnel

下载: 全尺寸图片 幻灯片

断层破碎带的规模根据地质构造的不同，从几米到几千米不等。对于任意断层破碎带形态都可以看成以下4种特征形态的线性组合，如图 2所示。其中相交段为形态一、三的线性组合，延伸段为形态二、四的线性组合。

图  2  断层破碎带地层典型特征形态

Figure  2.  Typical characteristics and morphologies of strata in fault fracture zone

下载: 全尺寸图片 幻灯片

1.2   隧道纵向上覆荷载计算

对断层破碎带范围内海底隧道任意横剖面分析，隧道的横向上覆荷载均可视为多层土体作用下隧道的松动土压力，力学模型如图 3所示。土体上作用水压力p₀，忽略水的渗流作用。从上至下3种土层厚度分别为h_a1，h_a2，h_a3。

图  3  隧道松动土压力计算模型

Figure  3.  Computational model for loosening earth pressure in tunnel

下载: 全尺寸图片 幻灯片

取第一层正常围岩松动区内深度y处宽为$2{B_t}$，厚度为dy的土条单元，根据土条竖向力的平衡，可得

式中：${\gamma _{\text{s}}}$为正常围岩的重度；${\sigma _y}$为深度y处的垂直方向应力；${\tau _{\text{s}}}$为作用于正常围岩地层土条上向上的摩擦应力；${B_t}$为隧道开挖影响宽度的一半，t=1，2，B₁为相交段范围内隧道开挖影响宽度的一半，B₂为其余段隧道开挖影响宽度的一半，${\tau _{\text{s}}}$和${B_{\text{t}}}$由下式计算：

式中：${\sigma _x}$为深度y+dy/2处的水平方向应力；c_s为正常围岩黏聚力；${\varphi _{\text{s}}}$为正常围岩内摩擦角；${\varphi _{\text{p}}}$为断层破碎带内摩擦角；r为隧道的外半径；k_s为正常围岩静止侧压力系数，按经验公式${k_{\text{0}}} = {\text{1}} - {\text{sin}}{\varphi _{\text{s}}}$进行计算^[16]。

将式(2)，(3)代入式(1)，结合边界条件$y=\text{0}，{\sigma }_{\text{y}}={p}_{\text{0}}$，并代入y=${h_{{\text{a1}}}}$可以得到第一层土体底部荷载：

同理，将上层土体底部荷载当作边界条件代入下层土体中进行计算，即可求得3层土体作用下隧道松动土压力：

式中：c_p为断层破碎带的黏聚力；${\gamma _{\text{p}}}$为断层破碎带的重度；k_p为断层破碎带静止侧压力系数，按经验公式${k_{\text{0}}} = {\text{1}} - {\text{sin}}{\varphi _{\text{p}}}$进行计算^[16]。

式(5)即为断层破碎带地层海底盾构隧道归一化上覆荷载理论公式。4种特征形态下海底盾构隧道纵向上覆荷载，如表 1所示。

表  1  4种特征形态下海底盾构隧道纵向上覆荷载

Table  1.  Longitudinal overlying loads of subsea shield tunnel for four characteristic forms

形态	纵向上覆荷载	${B_t}$取值
形态一	${q_{\text{1}}}{\text{(}}x{\text{)}} = q{\text{(}}x{\text{)}}\left| {{h_{{\text{a1}}}}{\text{ = }}{h_{\text{1}}}{\text{, }}} \right.{h_{{\text{a1}}}}{\text{ = }}{h_{\text{1}}}{\text{, }}{h_{{\text{a3}}}}{\text{ = 0}}$	${B_t} = {B_{\text{1}}}$
形态二	${q_{\text{2}}}{\text{(}}x{\text{)}} = q{\text{(}}x{\text{)}}\left| {{h_{{\text{a1}}}}{\text{ = 0, }}} \right.{h_{{\text{a2}}}}{\text{ = }}{h_3}{\text{, }}{h_{{\text{a3}}}}{\text{ = }}{h_4}$	${B_t} = {B_{\text{2}}}$
形态三	${q_3}{\text{(}}x{\text{)}} = q{\text{(}}x{\text{)}}\left| {{h_{{\text{a1}}}} = {h_{{\text{a3}}}} = 0, } \right.{h_{{\text{a2}}}} = {h_5}$	${B_t} = {B_{\text{1}}}$
形态四	${q_4}{\text{(}}x{\text{)}} = q{\text{(}}x{\text{)}}\left| {{h_{{\text{a1}}}} = {h_{\text{6}}}, {h_{{\text{a2}}}} = {h_7}, } \right.{h_{{\text{a3}}}} = {h_8}$	${B_t} = {B_{\text{2}}}$

下载: 导出CSV 

| 显示表格

对于正常段ab、de，海底盾构隧道的纵向上覆荷载q₀(x)：

结合正常段、相交段、延伸段的纵向上覆荷载，对于不同形态断层破碎带，纵向上覆荷载如表 2所示，其中H为隧道埋深，b_p为断层破碎带宽度，$\alpha$为断层破碎带与隧道轴线的夹角。

表  2  不同倾角断层破碎带海底隧道纵向上覆荷载

Table  2.  Longitudinal overlying loads of subsea tunnel in fault fracture zone with different inclination angles

倾角	上覆荷载
$\alpha = {90^ \circ }$	${q_0}(x{\text{)}}$$- \infty {\text{ < }}x{\text{ < 0}}$	${q_3}(x{\text{)}}$${\text{0}} \leqslant x{\text{ < }}{b_{\text{p}}}$	${q_0}(x{\text{)}}$${b_{\text{p}}} \leqslant x{\text{ < + }}\infty$
$\tan \alpha > \frac{H}{{{b_{\text{p}}}}}$	${q_0}(x{\text{)}}$$- \infty {\text{ < }}x{\text{ < 0}}$	${q_1}(x{\text{)}}$${\text{0}} \leqslant x{\text{ < }}\frac{H}{{{\text{tan}}\alpha }}$	${q_3}(x{\text{)}}$$\frac{H}{{{\text{tan}}\alpha }} \leqslant x \leqslant {b_{\text{p}}}$	${q_2}(x{\text{)}}$${b_{\text{p}}}{\text{ < }}x{\text{ < }}{b_{\text{p}}}{\text{ + }}\frac{H}{{{\text{tan}}\alpha }}$	${q_0}(x{\text{)}}$${b_{\text{p}}}{\text{ + }}\frac{H}{{{\text{tan}}\alpha }} \leqslant x < + \infty$
$\tan \alpha = \frac{H}{{{b_{\text{p}}}}}$	${q_0}(x{\text{)}}$$- \infty {\text{ < }}x{\text{ < 0}}$	${q_1}(x{\text{)}}$${\text{0}} \leqslant x{\text{ < }}{b_{\text{p}}}$	${q_3}(x{\text{)}}$$x = {b_{\text{p}}}$	${q_2}(x{\text{)}}$${b_{\text{p}}}{\text{ < }}x \leqslant {\text{2}}{b_{\text{p}}}$	${q_0}(x{\text{)}}$${\text{2}}{b_{\text{p}}}{\text{ < }}x{\text{ < + }}\infty$
$\tan \alpha < \frac{H}{{{b_{\text{p}}}}}$	${q_0}(x{\text{)}}$$- \infty {\text{ < }}x{\text{ < 0}}$	${q_1}(x{\text{)}}$${\text{0}} \leqslant x{\text{ < }}{b_{\text{p}}}$	${q_4}(x{\text{)}}$${b_{\text{p}}} \leqslant x < \frac{H}{{{\text{tan}}\alpha }}$	${q_2}(x{\text{)}}$$\frac{H}{{{\text{tan}}\alpha }} \leqslant x \leqslant {b_{\text{p}}}{\text{ + }}\frac{H}{{{\text{tan}}\alpha }}$	${q_0}(x{\text{)}}$${b_{\text{p}}}{\text{ + }}\frac{H}{{{\text{tan}}\alpha }}{\text{ < }}x{\text{ < }} + \infty$

下载: 导出CSV 

| 显示表格

2.   盾构隧道纵向沉降解析解

2.1   隧道纵向沉降理论模型

为了更加准确计算盾构隧道的纵向沉降与内力, 考虑管片接头对隧道整体剪切刚度的影响，将隧道管片结构视为置于Winkler地基上的Timoshenko梁，其纵向力学理论模型如图 4所示。其中f_j(x)为盾构隧道纵向上覆荷载，由表 3可知，隧道纵向上覆荷载f_j(x)共分为四类，即j=1~4。${k_{\text{2}}}$为断层破碎带地层地基基床系数；${k_{\text{1}}}$为正常地层地基基床系数。根据隧道纵向上覆荷载的分段及地层的变化，盾构隧道沿纵向可以分为多段Timoshenko梁。

图  4  盾构隧道管片结构纵向理论模型

Figure  4.  Longitudinal theoretical model for segment structure of shield tunnel

下载: 全尺寸图片 幻灯片

表  3  4种破碎带形态及其对应上覆荷载

Table  3.  Four types of fracture zones and their corresponding overlying loads

破碎带形态	$\alpha = {90^ \circ }$	$\tan \alpha > \frac{H}{{{b_{\text{p}}}}}$	$\tan \alpha = \frac{H}{{{b_{\text{p}}}}}$	$\tan \alpha < \frac{H}{{{b_{\text{p}}}}}$
上覆荷载	f1(x)	f2(x)	f3(x)	f4(x)

下载: 导出CSV 

| 显示表格

根据Timoshenko梁理论, 各段梁弯矩M_ij、剪力Q_ij、转角${\theta _{ij}}$满足下面的关系：

式中：(EI)_eq为隧道纵向等效抗弯刚度；(kGA)_eq为隧道纵向等效抗剪刚度；i为梁段号，${k_i}$为第i段梁的地基基床系数；M_ij为第j类上覆荷载作用在隧道上i段梁产生的弯矩。4类上覆荷载对应破碎带形态如表 3所示。

在纵向上覆荷载f_j(x)作用下，可以得到在Winkler地基上Timoshenko梁关于纵向沉降w_ij(x)的微分方程：

式(10)对应的齐次方程为

式中，${A_{i{\text{a}}}} = - \frac{{{\text{2}}{k_i}r}}{{{{{\text{(}}kGA{\text{)}}}_{{\text{eq}}}}}}$，${A_{i{\text{b}}}} = \frac{{{\text{2}}{k_i}r}}{{{{{\text{(}}EI{\text{)}}}_{{\text{eq}}}}}}$。

其齐次方程通解为

式中，${\alpha _i}{\text{ = }}\frac{{\sqrt {{\text{2}}\sqrt {{A_{i{\text{b}}}}} - {A_{i{\text{a}}}}} }}{{\text{2}}}$，${\beta _i} = \frac{{{\text{4}}{A_{i{\text{b}}}} - A_{i{\text{b}}}^{\text{2}}}}{{\sqrt {{\text{2}}\sqrt {{A_{i{\text{b}}}}} - {A_{i{\text{a}}}}} }}$，c_i1，c_i2，c_i3，c_i4为待定系数。

2.2   隧道纵向上覆荷载简化

由表 3可知，f_j(x)在其作用范围内任意开区间(x₀，x_n)具有直到n+1阶的导数。x_t为开区间(x₀，x_n)内的一点，对于任意x∈(x₀，x_n)，由泰勒公式得：

式中，$\xi$为x_t与x之间的某个值。当x→x_t时，误差$\left| {{R_n}{\text{(}}x{\text{)}}} \right|$是比${{\text{(}}x - {x_t}{\text{)}}^{n + 1}}$高阶的无穷小，x越趋近于x_t，误差就越小。利用泰勒公式这个性质，可以提出高精度泰勒公式分段展开法，即对不同的x进行不同x_t处高精度n次多项式展开。

本文将纵向上覆荷载f_j(x)展开为三次泰勒多项式，其与f_j(x)的误差为

对于有界函数f_j(x)，当x∈(x₀，x_n)时，$\mid f_j^{({\text{4}})}(x)\mid$ $\leqslant B$，则有下列不等式：

若本文最大误差取为5%，可以得到

联立求解式(16)，(17)，可以求出f_j(x)在x_t处最大误差为5%的展开区间(x₀，x_n)。即f_j(x)可以在5%误差内简化为三次多项式函数：

此时纵向沉降微分方程的特解为

即Timoshenko梁纵向沉降微分方程通解w_ij(x)为

假设隧道两端自由，可得边界条件：

假定隧道沉降曲线连续可微可导，在各段梁节点处隧道的弯矩、剪力、转角、纵向变形必然相等，可以列出以下方程：

式中，${x_{\text{p}}}$为各连接点横坐标。

由式(21)~(26)及式(7)~(9)，(12)可以求出通解中待定系数c_i1，c_i2，c_i3，c_i4。

2.3   相关参数的确定

隧道纵向等效抗弯刚度(EI)_eq及纵向等效剪切刚度(kGA)_eq是决定隧道变形及内力分布的重要参数，其等效纵向抗弯刚度^[17]表达式为

式中：E_c为管片结构弹性模量；I_c为管片环纵向惯性矩；$\phi$为中性轴位置的角度，由下式确定：

式中：k_b接头螺栓的平均线刚度，${k_{\text{b}}}{\text{ = }}{E_{\text{b}}}{A_{\text{b}}}{\text{/}}{l_{\text{b}}}$，E_b为螺栓的弹性模量，A_b为螺栓的横截面积, l_b为螺栓长度；n为螺栓个数；l_s为管片环宽；A_c为隧道管片横截面积。

盾构隧道等效剪切刚度(kGA)_eq公式^[13]为

式中：$\zeta$为修正系数，在本文其值取1；${\kappa _{\text{b}}}$及${\kappa _{\text{c}}}$分别为螺栓及管片环Timoshenko剪切系数，对于圆形截面螺栓，${\kappa _{\text{b}}}$取0.9，对于环形隧道管片环结构，${\kappa _{\text{c}}}$取0.5；G_b及G_c为分别螺栓剪切刚度及隧道管片剪切刚度，两者与各自弹性模量的关系为${G_{\text{b}}} = \frac{{{E_{\text{b}}}}}{{2(1 + {\nu _{\text{b}}})}}$及${G_{\text{c}}} = \frac{{{E_{\text{c}}}}}{{{\text{2(1 + }}{\nu _{\text{c}}}{\text{)}}}}$，其中，${\nu _{\text{b}}}$及${\nu _{\text{c}}}$分别为螺栓及管片的泊松比。

Winkler地基基床系数k表达式^[18]为

式中：E_s为正常围岩地层弹性模量；$\nu _{\text{s}}^{}$为正常围岩地层泊松比；E_p为断层破碎带地层弹性模量；$\nu _{\text{p}}^{}$为断层破碎带地层泊松比。

3.   工程实例

3.1   工程概况

珠三角水资源配置工程狮子洋段海底输水隧洞长2252 m，埋深40 m，采用泥水平衡盾构法施工，最大水深为27 m。管片采用C50钢筋混凝土管片错缝拼装，外径8.3 m，内径7.5 m，管片环宽度1.6 m。隧洞穿越f114断层破碎带，该断面水深10 m。断层宽度28.5 m，倾角约为40°。表 4给出断层破碎带与周边岩层的力学参数，其中$\gamma$为重度，$\varphi$为内摩擦角，c为黏聚力，E为弹性模量，$\nu$为泊松比，k为静止侧压力系数。断层破碎带的形态及其与盾构掘进面相对位置关系如图 6所示。

表  4  断层破碎带及围岩力学参数

Table  4.  Mechanical parameters of fault fracture zone and surrounding rock

类型	$\gamma {\text{/}}$${\text{(kN}} \cdot {{\text{m}}^{ - {\text{3}}}}{\text{)}}$	$\varphi$/(°)	$c{\text{/}}$kPa	$E{\text{/}}$GPa	$\nu$	k
正常围岩	26	30	100	1.6	0.33	0.49
断层破碎带	20	25	30	0.5	0.40	0.67

下载: 导出CSV 

| 显示表格

图  6  断层破碎带与隧道剖面图

Figure  6.  Section of fault fracture zone and tunneling position

下载: 全尺寸图片 幻灯片

由表 3可知，隧道上覆荷载为第四类f₄(x)。根据破碎带的上覆荷载及地层的变化将隧道分为5段Timoshenko梁，其范围分别为：(-∞, 0]，(0, b_p]，(b_p，H/tan$\alpha$]，(H/tan$\alpha$，b_p+H/tan$\alpha$]，(b_p+H/tan$\alpha$，+∞)。纵向沉降微分方程式(20)相关的系数保留两位小数后分别为：c₁₁=-3.49×10^-9；c₁₂=7.04×10^-9；c₁₃=0；c₁₄=0；c₂₁=4.01×10^-10；c₂₂=-2.54×10^-10；c₂₃=-14.76；c₂₄= 17.53；c₃₁=-1.26×10^-23；c₃₂=-3.22×10^-23；c₃₃=2.0×10⁹；c₃₄=8.75×10⁸；c₄₁=-1.14×10^-30；c₄₂=-8.73×10^-32；c₄₃=1.76×10¹¹；c₄₄=3.60×10¹¹；c₅₁=0；c₅₂=0；c₅₃=7.22×10¹⁸；c₅₄=-5.95×10¹⁷。

3.2   监测点布设

根据工程地质条件，选取5个代表管片，分别为处于断层破碎带范围内的第1529、1534和1538环，断层破碎带之外的第1551、1561环。采用全分布式光纤传感器对管片进行监测，如图 7所示。

图  7  光纤传感器分布示意图

Figure  7.  Distribution of optical fiber sensors

下载: 全尺寸图片 幻灯片

3.3   现场监测试验结果对比分析

采用偏压短柱法^[19]通过实测钢筋应变反算盾构管片内力。如图 8所示，将隧洞衬砌考虑成偏心受压短柱，建立力矩平衡方程，计算相关参数见表 5。

图  8  管片内力计算示意图

Figure  8.  Calculation of internal force of segment

下载: 全尺寸图片 幻灯片

表  5  螺栓及钢筋相关参数

Table  5.  Relevant parameters of bolt and reinforcement

a'/mm	Ee/GPa	n1	n2	Ec/GPa
50	210	20	10	3.15×104
nb	Eb/GPa	lb/mm	Ab/m2	${\alpha _{\text{b}}}{\text{/}}$(°)
28	210	730	7.069×10-4	35
注：$a'$为钢筋保护层厚度；Ee为钢筋弹性模量；n1为管片内侧环向主筋数量；n2为管片外侧环向主筋数量；Ec为管片弹性模量；nb为纵向螺栓个数；Eb为螺栓弹性模量；lb为螺栓长度，Ab为螺栓横截面面积；${\alpha _{\text{b}}}$为螺栓倾角。

下载: 导出CSV 

| 显示表格

式中：${P_{\text{1}}}$为管片外侧单根钢筋集中力；${P_{\text{2}}}$为管片内侧单根钢筋集中力；${\sigma _{{\text{c1}}}}$为管片内侧钢筋位置处混凝土应力；${\sigma _{{\text{c2}}}}$为管片外侧钢筋位置处混凝土应力；${n_{\text{c}}}$为管片内(外)侧环向主筋数量；b为管片厚度；$a'$为钢筋保护层厚度；A_c为管片横截面面积。

对第1534、1538环管片结构的内力进行分析，管片环宽度为1.6 m，同一环管片上理论计算弯矩值会在一定的范围波动，其与监测弯矩对比结果如图 9，10所示。

图  9  1529环管片实测弯矩与计算弯矩

Figure  9.  Measured and calculated bending moments of segments of No. 1529 ring

下载: 全尺寸图片 幻灯片

图  10  1538环管片实测弯矩与计算弯矩

Figure  10.  Measured and calculated bending moments of segments of No. 1538 ring

下载: 全尺寸图片 幻灯片

由图 9，10可知：现场实测得到的隧洞弯矩与理论计算较为接近。1529环根据现场实测钢筋应变反算管片弯矩值在212.2~228.5 kN，在监测后期管片的弯矩值基本稳定在220 kN·m左右。1538环根据现场实测钢筋应变反算管片弯矩值在175.4~193.6 kN·m，1538环管片弯矩值稳定在180 kN·m左右。

4.   参数敏感性分析

4.1   断层破碎带倾角

假定隧道埋深40 m，断层破碎带宽度为30 m，水深10 m，分别对断层破碎带倾角为20°，40°，60°和80°条件下隧道纵向沉降进行分析，如图 11所示。在断层破碎带区域隧道沉降先急剧增大，后缓慢增加、最后急剧减小，即曲线整体呈现三段变化规律：激增段、缓增段和骤减段。在断层破碎带与正常围岩界面附近隧道不均匀沉降急剧增加。

图  11  不同断层破碎带倾角隧道纵向沉降图

Figure  11.  Longitudinal settlements of tunnel with different inclination angles of fault fracture zone

下载: 全尺寸图片 幻灯片

4.2   断层破碎带宽度

假定隧道埋深40 m，断层破碎带倾角为60°，分别对断层破碎带宽度为20，30，40，50，60 m条件下隧道纵向沉降进行分析，如图 12所示。随着断层破碎带宽度的增大，断层破碎带范围内隧道的最大纵向沉降逐渐增大并趋于稳定，最大沉降出现位置逐渐远离断层破碎带中部。对于8.3 m外径海底盾构隧道，倾斜断层破碎带对延伸侧1倍破碎带宽度范围内隧道的沉降影响较大，当断层破碎带宽度超过50 m时，隧道的最大纵向沉降基本保持不变。

图  12  不同断层破碎带宽度隧道纵向沉降图

Figure  12.  Longitudinal settlements of tunnel with different widths of fault fracture zone

下载: 全尺寸图片 幻灯片

不同断层破碎带宽度对隧道纵向弯矩影响如图 13所示。随着断层破碎带宽度增加，在断层破碎带和正常围岩界面附近，隧道纵向弯矩基本不变，在断层破碎带内部隧道的纵向弯矩逐渐减小。对于8.3 m外径海底盾构隧道，当破碎带宽度超过40 m，断层破碎带中部几乎无弯矩。

图  13  不同断层破碎带宽度隧道弯矩图

Figure  13.  Bending moments of tunnel with different widths of fault fracture zone

下载: 全尺寸图片 幻灯片

4.3   断层破碎带弹性模量

假定隧道埋深40 m，断层破碎带宽度为30 m，倾角为60°，分别对断层破碎带与正常围岩弹性模量比值为1∶2，1∶4，1∶6，1∶8进行分析，如图 14所示。随着断层破碎带与正常围岩弹性模量比值的减小，隧道沉降趋势基本不变，但不均匀沉降迅速增大。对于外径为8.3 m海底盾构隧道，当宽度为30 m的断层破碎带与正常围岩弹性模量比值小于1∶4时，隧道的沉降量超过预警值10 mm^[20]。

图  14  不同断层破碎带弹性模量隧道纵向沉降图

Figure  14.  Longitudinal settlements of tunnel with different elastic moduli of fault fracture zone

下载: 全尺寸图片 幻灯片

5.   结论

本文针对断层破碎带海底盾构隧道纵向沉降问题，基于太沙基松动土压力理论，结合泰勒级数展开公式，对断层破碎带海底盾构隧道纵向沉降及内力进行分析，得到以下3点结论。

(1) 通过分析海底盾构隧道穿越断层破碎带地层全过程，提出了可考虑断层破碎带宽度和倾角的统一荷载模型和理论公式。

(2) 针对复杂上覆荷载条件下Timoshenko梁难以直接求解问题，利用泰勒级数展开公式提出了高精度泰勒公式分段展开法，结合Timoshenko梁理论，推导建立了可用于复杂荷载作用下海底盾构隧道管片纵向沉降及内力系列解析解。工程实例表明本文理论推导的系列解析解能够较为准确地计算管片结构内力。

(3) 计算结果表明，断层破碎带内隧道的沉降曲线表现为激增段、缓增段和骤减段。断层破碎带对延伸侧1倍破碎带宽度范围内隧道的沉降影响较大；当断层破碎带宽度超过6倍洞径时，隧道的沉降峰值基本保持不变。