摘要
大直径盾构隧道面临着复杂地层带来的挑战,传统的设计方法往往选择钻孔揭示的最不利地层断面或者采用平均化的地层,无法准确反映掘进面临的真实地层情况.因此,利用现有钻孔勘察信息建立耦合马尔可夫链模型,确定盾构掘进过程中最可能穿越的典型地层断面.采用不同设计参数展开数值模拟,以盾构隧道收敛变形相关指标的标准差作为鲁棒性指标,提出了一种基于地层变异性的大直径盾构隧道鲁棒性设计方法.研究结果表明,耦合马尔可夫链模型的最典型断面为盾构隧道提供了更准确的地层条件,使得鲁棒性设计更具适用性.最后,对现阶段常用的地层设计方法进行了优化,确保掘进过程中地层选择更接近实际情况,在有效降低隧道建设成本的同时,使得隧道具备更好的鲁棒性.
关键词
大直径盾构隧道穿越多重地层,导致其受力条件复
目前国内外众多学者对地层变异性开展了大量研究.Deng
已有研究结合不确定性分析和优化模型,开展鲁棒性设
为了解决盾构隧道穿越地层形态变化的问题,提出一种考虑地层变异性和岩土参数变异性的盾构隧道鲁棒性设计方法.结合佛山市季华路地铁线路的设计资料,并充分利用现有的钻孔信息.首先,建立基于耦合马尔可夫链的地层随机场模型,通过概率统计分析得出盾构隧道最可能穿越的典型地层断面.然后在此基础上,引入岩土参数的变异性,并在保证结构安全可靠的条件下,以隧道收敛变形的相关指标作为鲁棒性指标,评估地层变异性对盾构隧道鲁棒性的影响程度,从而证明地层模拟的可行性和适应性.最后,对盾构隧道实际设计方法进行优化,综合考虑掘进过程和钻孔信息,为大直径盾构隧道的设计提供具有可行性的建议.
1 基于钻孔条件的地层变异性模型
Krumbein于1968年提出马尔可夫链模型在模拟地层序列中的运用,Park
为更真实地模拟地层变异性,采用Park

图1 耦合马尔可夫链状态表示
Fig.1 Coupled Markov state representation
(1) |
式中:为单元(m,n)处于Sk状态的概率;Ny为最右边单元格的列号;假设单元(m-1,n)、(m,n-1)和(Nx,n-1)分别处于Sl、St和Sk状态;和分别为水平方向上(Nx-m)步长从Sk、Sf到Sq的概率.
Deng
(2) |
该定律指出,不同地层在垂直和水平方向的沉积顺序相同,因此不同类型地层的转换次数相同(=,m≠n).由于在相同的地质环境中共同沉积,所以假设地层的平均倾斜度都相同,令K=/, K是唯一未知的参数.所有单元从一个状态转移到其自身的概率总是高于其他状态.如果钻孔揭示三种地层,则垂直转移计数矩阵如
(3) |
耦合马尔可夫链模型的具体步骤如下:
1)利用钻孔信息计算垂直转移概率矩阵.
2)假设Nk个K值(K1,K2,…,)用于计算水平转移概率矩阵.
3)模拟条件由两侧距离最近的钻孔提供.
4)由每个单元状态概率的积表示钻孔岩土类型出现的概率Pl,即
(4) |
式中:P(Zm,n=)为单元(m,n)的状态转化为的概率;kn为单元在竖直方向上的个数.取Nk个结果中Pl最大时对应的K值,用于后续地层模拟.蒙特卡罗模拟地层时一般取1 500次即可满足精度要
2 基于随机场的盾构穿越地层模型
在考虑地层变异性的基础上得到最大穿越概率的典型横断面,引入一次二阶矩法来描述每一层岩土性质的内在变异性,真实反映盾构隧道穿越地层多维度性质的变化.
2.1 工程概况
采用佛山市季华路西延线隧道的勘察设计资料.综合考虑各种因素,确定盾构隧道的内径为13.7 m,管片厚度为0.65 m.在掘进起始方向的480 m范围内,共有12个钻孔,编号从SD059到SD128.
如

图2 纵断面钻孔图
Fig.2 Longitudinal borehole profile
地层 | 压缩模量范围/MPa | 均值/MPa | 变异系数 |
---|---|---|---|
中风化泥岩 | 4.50~9.29 | 6.54 | 0.55 |
中风化砂岩 | 2.72~24.70 | 12.39 | 0.33 |
中风化砂砾岩 | 2.72~23.70 | 13.95 | 0.44 |
中风化玄武岩 | 15.80~58.40 | 33.23 | 0.48 |
2.2 基于耦合马尔可夫链的地层随机场模型
利用耦合马尔可夫链模型和钻孔数据,结合蒙特卡罗模拟,预测地层在水平和竖直方向上的分布,并确定最佳参数以进行地层模拟.
2.2.1 钻孔信息转换
为了建立一个二维耦合马尔可夫链模型,用于描述盾构穿越的地层,将所有钻孔投影到盾构掘进方向上,并重新标记为D1~D12.以

图3 钻孔相对位置
Fig.3 Relative position of borehole
在所有的钻孔中,钻孔长度从-24.5 m到20 m不等.水平方向上,钻孔之间的最大间距为13 m.投影结果如
(5) |
2.2.2 耦合马尔可夫链模拟
根据
地层状态 | 黏土 | 砂岩 | 砾岩 |
---|---|---|---|
黏土 | 0.867 8 | 0.086 2 | 0.045 98 |
砂岩 | 0.130 4 | 0.826 1 | 0.043 48 |
砾岩 | 0 | 0 | 1 |
地层状态 | 黏土 | 砂岩 | 砾岩 |
---|---|---|---|
黏土 | 0.836 7 | 0.163 3 | 0 |
砂岩 | 0.399 8 | 0.600 2 | 0 |
砾岩 | 0.194 9 | 0.091 4 | 0.713 7 |
2.2.3 预测结果
K为8~12时,地层模拟效果较
经验证,K取8.4时Pl对应最大,地层模拟效果最好.其中一次模拟的地层如

(a) 1~110 m地层预测

(b) 111~240 m地层预测

(c) 241~360 m地层预测

(d) 361~480 m地层预测
图4 耦合马尔可夫链预测
Fig.4 Coupled Markov formation prediction
2.2.4 准确性验证
提取钻孔Jz-Ⅲ20-S18和SD116钻孔信息,验证模拟地层真实性,选取情况如
钻孔 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 |
---|---|---|---|---|---|---|---|---|---|---|---|---|
选取 | √ | √ | √ | √ | √ | √ | √ | √ | √ | √ |

(a) 钻孔Jz-Ⅲ20-S18预测

(b) 钻孔SD116预测
图5 耦合马尔可夫链模型验证
Fig.5 Coupled Markov formation prediction
2.2.5 典型穿越地层选择
通过勘察设计资料,获得盾构掘进过程中的埋深变化以及相应的掘进长度.在耦合马尔可夫链模拟的地层中,结合盾构穿越的区间进行统计分析.蒙特卡罗模拟生成的地层情况存在多种可能性,但盾构最常穿越的典型地层断面分布情况趋于一致.
如

图6 盾构穿越埋深对应
Fig.6 Shield crossing buried depth
3 基于模拟地层的鲁棒性设计模型
利用最常穿越地层和相应岩土参数的力学特性,为数值模拟提供更准确的地层条件.在确保盾构隧道安全的前提下,评估与收敛变形相关的鲁棒性指标,并使用多目标优化方法获取适宜的鲁棒性设计方案.
3.1 力学模型
Zhang

图7 有限元模型
Fig.7 Finite element model
3.2 鲁棒性指标
《盾构法隧道施工及验收规范》(GB 50446—2017
(6) |
大直径盾构隧道在穿越复合地层时不可避免地面临复杂多变的地质条件.为了有效评估衬砌结构对这种地层的适应能力,选取安全系数Fs的变异系数δ作为鲁棒性指标,如
(7) |
3.3 多目标优化模型
依据上述设计流程与背景,将优化的问题简化为二维优化问题.优化参数为管片厚度h、管片配筋率ρ和横向接头螺栓直径Db,结合设计资料,建立成本函数如
(8) |
引入各类成本参数,符号说明见
符号 | 意义 | 单位 |
---|---|---|
Pc | 衬砌混凝土单价 |
元/ |
Ps | 钢筋单价 | 元/kg |
Pb | 接头螺栓单价 | 元/kg |
Vc | 混凝土用量 |
|
ms | 钢筋用量 | kg/m |
mb | 接头螺栓总质量 | kg/m |
根据现场的设计概算资料,取衬砌混凝土单价1 200元/
(9) |
(10) |
结合定价和盾构设计情况,得出每延米所需管片、钢筋以及螺栓的单价,具体如
种类 | 数量 | 每延米数量 | 每延米费用/元 |
---|---|---|---|
管片混凝土 |
44 335.16 |
30.118 | 22 864.31 |
钢筋 | 7 282.02 t | 5.045 15 t | 12 836 |
8.8级M36螺栓 | 31 458套 | 21.37套 | 892.872 |
初始方案是季华路隧道项目的实际设计参数:管片配筋(6ϕ25+8ϕ22)/环,ρ=1.194%,h=0.65 m, Db=36 mm.盾构隧道段横断面设计综合考虑建筑限界、设备限界及综合误差等因素.管片厚度设计过小,则衬砌环变形量越大,对施工中的拼装和隧道防水越不利;管片厚度设计过大,则会影响建筑限界.综合考虑大直径盾构隧道断面的特点,选取各待优化设计参数的上下界如
设计参数 | 限值 |
---|---|
h/m | [0.40,0.60] |
ρ/% | [0.2,4.0] |
Db/mm |
[ |
结合数值模拟得到隧道在不同设计下的鲁棒性,对算例参数进行影响性分析和归一化处理后,得到设计与鲁棒性指标的表达式,如
(11) |
每一种设计都有对应的鲁棒性和成本,成本越高相应鲁棒性越好,但影响程度难以评估.引入双目标优化方法权衡鲁棒性和经济性,筛选不同Pareto组中的最优解,获得完整的Pareto前沿面.
3.4 不同地层选择方法
从工程实际设计方法出发,不考虑地层变异性.第一种方法选取钻孔统计得到的最不利地层进行设计.第二种对比方法采用线性插值法,直接得到地层.最后统一对比不同设计的鲁棒性,在现阶段设计基础上提出有效建议.
4 参数研究
4.1 基于地层模拟的鲁棒性分析的盾构最优设计
通过NSGA-Ⅱ算法获得的Pareto前沿如

图8 基于NSGA-Ⅱ优化的Pareto前沿
Fig.8 Pareto based on NSGA-Ⅱ optimization
鲁棒性指标从0.05增加到0.062时,衬砌费用迅速下降至8.606 3万元,而当鲁棒性指标从0.062增加到0.09时,衬砌费用变化不明显.故Knee Point具有较高的鲁棒性,并且与Zhang
设计点 | h/mm | ρ/% | Db/mm | Pm/万元 | δ |
---|---|---|---|---|---|
实际点 | 0.65 | 1.194 | 36 | 9.791 | 0.077 6 |
优化点1 | 0.60 | 1.036 | 36 | 9.791 | 0.049 1 |
优化点2 | 0.60 | 1.4 | 34 | 7.981 4 | 0.077 6 |
最优点 | 0.60 | 1.6 | 36 | 8.606 3 | 0.060 59 |
4.2 不同地层设计选择对鲁棒性的影响
使用和3.3节中一致的设计区间,结合工程设计中较常采用的两种地层选择方法,即基于钻孔所得的最不利地层模型和加权平均地层模型,对盾构隧道展开鲁棒性设计.
如

(a) 最不利地层模型

(b) 加权平均地层模型
图9 考虑不同地层条件下的Pareto前沿
Fig.9 Pareto frontier under stratigraphic conditions
如
设计点 | h/mm | ρ/% | Db/mm | Pm/万元 | δ |
---|---|---|---|---|---|
最不利地层 | 0.62 | 1.736 | 36 | 9.791 0 | 0.071 64 |
加权平均地层 | 0.55 | 1.62 | 34 | 8.706 4 | 0.068 82 |
4.3 基于工程设计的优化建议
蒙特卡罗模拟需要大量计算,为了提高效率,综合考虑工程中常采用的加权平均地层、地铁隧道穿越埋深情况.采用2.2.5节中提出的典型穿越断面选择方法,得到最易穿越的断面如

图10 基于实际设计的优化地层选择
Fig.10 Optimal formation selection based on actual design
该断面状态下隧道中心埋深21.6 m,上覆土为黏土,隧道穿越砂岩和砾岩,岩性转换边界对应的标高为-23.45 m,与耦合马尔可夫链模型模拟的结果相近.
如

图11 优化地层选择的Pareto前沿
Fig.11 Pareto frontier for optimal formation selection
5 结 论
通过对具体工程案例的分析,利用钻孔信息进行地层模拟,得到大直径盾构隧道最可能经过的典型断面.在考虑结构对岩土参数变异敏感性和大直径盾构隧道建设成本的基础上,形成一个约束目标优化问题,建立了一套大直径盾构隧道的鲁棒性设计框架.得到的主要结论如下:
1) 大直径盾构隧道所处地层环境复杂,通过结合钻孔勘探的信息,能更精确地预测盾构开挖阶段所经过的地层分布情况,更加准确地描述盾构周边地层条件.
2) 盾构隧道穿越区间的地层边界和内部岩土参数存在较大变异性,选取鲁棒性指标时考虑安全系数的响应,精确评价结构的鲁棒性,进而结合设计进行优化.
3) 盾构隧道横断面设计时,应引入穿越地层变化的理念,综合考虑掘进穿越过程和钻孔所得条件,提升盾构隧道的鲁棒性,并保证计算效率,为大直径盾构隧道的设计提供可靠的依据.
4) 为大直径盾构隧道设计提供了切实可行的鲁棒性设计框架,并对大直径盾构隧道地层分布做出了较准确的预测,为工程提供了有价值的参考.
参考文献
JUANG C H, ZHANG J, SHEN M F,et al.Probabilistic methods for unified treatment of geotechnical and geological uncertainties in a geotechnical analysis[J].Engineering Geology,2019, 249: 148-161. [百度学术]
WANG C H,WANG K,TANG D F,et al.Spatial random fields-based Bayesian method for calibrating geotechnical parameters with ground surface settlements induced by shield tunneling[J].Acta Geotechnica,2022,17(4):1503-1519. [百度学术]
DENG Z P,LI D Q,QI X H,et al.Reliability evaluation of slope considering geological uncertainty and inherent variability of soil parameters[J].Computers and Geotechnics, 2017,92:121-131. [百度学术]
CHESSA A G. A Markov chain model for subsurface characterization: theory and applications[J]. Mathematical Geology,2006,38(4):503-505. [百度学术]
HUANG H W,GONG W P,KHOSHNEVISAN S,et al.Simplified procedure for finite element analysis of the longitudinal performance of shield tunnels considering spatial soil variability in longitudinal direction[J].Computers and Geotechnics,2015, 64: 132-145. [百度学术]
GONG W P,ZHAO C,JUANG C H,et al.Stratigraphic uncertainty modelling with random field approach[J].Computers and Geotechnics,2020,125:103681. [百度学术]
ZHAO C,GONG W P,LI T Z,et al.Probabilistic characterization of subsurface stratigraphic configuration with modified random field approach[J].Engineering Geology,2021,288:106138. [百度学术]
DENG Z P,JIANG S H,NIU J T,et al.Stratigraphic uncertainty characterization using generalized coupled Markov chain[J].Bulletin of Engineering Geology and the Environment, 2020, 79(10): 5061-5078. [百度学术]
张东明, 代鉷锋, 王慧, 等.考虑地层变异的浅基础承载力分析[J].地下空间与工程学报, 2020, 16(5): 1412-1419. [百度学术]
ZHANG D M,DAI H F,WANG H,et al.Analysis on shallow foundation capacity considering geological heterogeneity[J].Chinese Journal of Underground Space and Engineering, 2020, 16(5): 1412-1419.(in Chinese) [百度学术]
ZHANG D M,DAI H F,WANG H,et al.Investigating the effect of geological heterogeneity of strata on the bearing capacity of shallow foundations using Markov random field[J]. ASCE-ASME Journal of Risk and Uncertainty in Engineering Systems,Part A:Civil Engineering, 2021, 7(4): 04021060. [百度学术]
赵森森, 张冬梅, 黄忠凯.大直径盾构隧道分布式凹凸榫受力特性[J].现代隧道技术,2022,59(5):54-62. [百度学术]
ZHAO S S,ZHANG D M,HUANG Z K.Mechanical characteristics of distributed mortise-and-tenon of large-diameter shield tunnels[J].Modern Tunnelling Technology,2022,59(5):54-62.(in Chinese) [百度学术]
张冬梅,刘杰,李保军,等.大直径盾构隧道斜螺栓环缝抗剪特性研究[J].中国公路学报,2020,33(12):142-153. [百度学术]
ZHANG D M,LIU J,LI B J,et al.Shearing behavior of circumferential joints with oblique bolts in large diameter shield tunnel[J].China Journal of Highway and Transport,2020, 33(12): 142-153.(in Chinese) [百度学术]
ZHANG D M,ZHOU W D,BU X H,et al.Failure mechanism and stiffness degradation of double lining with inner R/FRC lining subjected to internal water pressure[J]. Tunnelling and Underground Space Technology,2022,130:104737. [百度学术]
ZHANG D M,CHEN C C,ZHANG D M.Ground surface movement of shallow-buried large-sectional tunnel under full-ring pipe-jacking roof and ground freezing[J].Tunnelling and Underground Space Technology,2022,127:104600. [百度学术]
ZHANG D M,CHEN S,WANG R C,et al.Behaviour of a large-diameter shield tunnel through multi-layered strata[J].Tunnelling and Underground Space Technology,2021,116:104062. [百度学术]
李守巨,刘军豪,刘文怡.琼州海峡隧道盾构管片内力分析及截面设计[J].辽宁工程技术大学学报(自然科学版),2018,37(3):547-552. [百度学术]
LI S J,LIU J H,LIU W Y.Internal force analysis and section design for shield tunnel segments in Qiongzhou strait[J].Journal of Liaoning Technical University (Natural Science),2018,37(3):547-552.(in Chinese) [百度学术]
朱瑶宏,夏杨于雨,董子博,等.通用环错缝拼装盾构隧道结构设计计算参数研究[J].现代隧道技术,2018,55(4):114-121. [百度学术]
ZHU Y H,XIA Y Y Y,DONG Z B,et al.Calculation parameters of a shield tunnel structure with a staggered segmental tunnel lining[J].Modern Tunnelling Technology,2018,55(4):114-121.(in Chinese) [百度学术]
赵子蓬, 刘国彬, 白云.高水压大直径盾构隧道管片设计[J].公路交通技术, 2018, 34(增刊1):98-106. [百度学术]
ZHAO Z P, LIU G B, BAI Y.Design of the shield tunnel segments with large diameter in high water pressure[J].Technology of Highway and Transport, 2018, 34(Sup.1): 98-106.(in Chinese) [百度学术]
DODIGOVIĆ F, IVANDIĆ K, KOVAČEVIĆ M S, et al.Modified, reliability-based robust geotechnical design method, in accordance with eurocode 7[J]. Applied Sciences, 2021, 11(18):8423. [百度学术]
GAO K,DO D M, CHU S,et al.Robust topology optimization of structures under uncertain propagation of imprecise stochastic-based uncertain field[J].Thin-Walled Structures, 2022, 175:109238. [百度学术]
陈坤, 黄宏伟, 张东明,等. 基于约束多目标优化算法的盾构隧道设计方法[J]. 土木工程学报, 2020, 53(增刊1): 81-86. [百度学术]
CHEN K, HUANG H W, ZHANG D M, et al. Constrained multi-objective optimization algorithm based design method of shield tunnel[J]. China Civil Engineering Journal, 2020, 53(Sup.1): 81- 86. (in Chinese) [百度学术]
华雨杉, 黄宏伟, 陈坤, 等. 双层衬砌输水盾构隧道鲁棒性设计方法[J]. 湖南大学学报(自然科学版), 2022, 49(9):145-155. [百度学术]
HUA Y S,HUANG H W,CHEN K,et al.Robust design method of water conveyance shield tunnel with double lining[J].Journal of Hunan University (Natural Sciences), 2022, 49(9): 145-155.(in Chinese) [百度学术]
黄大维, 冯青松, 刘开富, 等. 基于弯矩最小化的地铁盾构隧道横断面优化设计分析[J]. 铁道学报, 2018, 40(9): 159-166. [百度学术]
HUANG D W, FENG Q S, LIU K F,et al.Analysis for optimal design of cross section of metro shield tunnel based on minimum bending moment[J]. Journal of the China Railway Society, 2018, 40(9): 159-166.(in Chinese) [百度学术]
李春良, 王方彦, 张立辉, 等.复杂条件下的地铁盾构隧道的力学行为研究[J].铁道学报, 2016, 38(5): 103-109. [百度学术]
LI C L,WANG F Y,ZHANG L H,et al.Study on mechanical behavior of metro shield tunnel under complex conditions[J].Journal of the China Railway Society,2016,38(5):103-109.(in Chinese) [百度学术]
PARK E,ELFEKI A,DEKKING M.Characterization of subsurface heterogeneity:integration of soft and hard information using multidimensional coupled Markov chain approach[J].Developments in Water Science,2005,52:193-202. [百度学术]
DO N A, DIAS D, ORESTE P, et al.2D numerical investigation of segmental tunnel lining behavior[J]. Tunnelling and Underground Space Technology, 2013, 37: 115-127. [百度学术]