逆断层错动作用下地铁隧道结构损伤分析

安韶 陶连金 边金 韩学川 吴晓娲

摘   要：以乌鲁木齐地铁隧道穿越西山活动逆断层工程为例，建立三维弹塑性有限元模型. 首先模拟分析了逆断层错动作用下隧道二次衬砌塑性应变发展过程，拉压损伤因子、剪切应变的横向及纵向分布规律，计算了混凝土的裂缝宽度;其次研究了不同错动位移、隧道底部距围岩交界面不同垂直距离及不同破碎带宽度的结构损伤规律，最后进行了设置柔性接头的减灾效果研究. 结果表明：二次衬砌结构破坏首先出现在拱顶;然后是拱底，最后在拱腰处累积. 破裂面附近拱腰处发生拉压剪的共同破坏;远离破裂面上盘拱顶，破碎带拱底处发生受拉破坏;远离破裂面上盘拱底，破碎带拱顶处发生受压破坏. 基于混凝土裂缝得到隧道拉裂破坏的严重与轻微受损区分别为10 m和30 m. 错动位移越大，结构受损越严重;隧道底部距围岩交界面垂直距离越大，土层越厚，耗散能量越多，结构受损越轻;破碎带宽度越大，隧道破坏越严重，当破碎带宽度达到26 m时，破碎带宽度对隧道的影响基本保持稳定. 设置柔性接头可以显著降低结构的损伤，基本满足在设防错动位移下的设计要求.

關键词：逆断层;地铁隧道;损伤分析;影响因素;柔性接头

中图类分号：TU443                                文献标志码：A

Abstract：DamageAnalysis on Subway Tunnel Structure under Effect of Reverse Fault DislocationAN Shao1，TAO Lianjin1，BIAN Jin2，HAN Xuechuan1，WU Xiaowa1 （1. Beijing Key Laboratory of Earthquake Engineering and Structural Retrofit，Beijing University of Technology，Beijing 100124，China;2. Guang Dong Ocean University，Zhanjiang 524008，China）Abstract：Taking the subway tunnel in Urumqi crossing Xishan reverse active fault as an example，a three- dimensional elastic-plastic finite element model of soil-tunnel structure was established to simulate and analyze plastic strain development process of secondary lining，transverse and longitudinal distribution of tensile damage factor，compressive damage factor and shear strain，calculate the crack width of concrete under the action of reverse fault dislocation. Then the damage law of the tunnel structure with different dislocation displacements，different vertical distances from tunnel bottom to the interface of surrounding rocks and different widths of fault zone was studied. Finally，the disaster reduction effect of the flexible joint was studied. Results show that the damage of secondary lining first appeared in the vault，then the bottom，and finally accumulated at the waist. The shear，tensile and compressive damage all appear at the tunnel waist of the rupture surface. When the tunnel structure is away from rupture surface，the tensile damage appears at the tunnel vault of hanging wall and invert of fault zone，the compressive damage appears at the tunnel invert of hanging wall and vault of fault zone. The severe and slightly tensile-crack damaged length are 10 m and 30 m respectively based on the crack width of concrete. The larger the fault displacement is，the more serious the structural damage is. The larger the vertical distances from tunnel bottom to the interface of surrounding rocks，the thicker the soil is，resulting in that more energy is dissipated and the structural damage is lighter. The damage to tunnel increases with the increase of width of fault zone，but when the width increases to 26 m，the influence of width of fault zone on tunnel damage becomes stable. Setting flexible joints can significantly reduce structural damage and substantially meet the design requirements for fault displacement.

Key words：reverse fault;subway tunnel;damage analysis;influence factor;flexible joint

目前，我国很多城市正在进行大规模的地铁建设. 但如北京、乌鲁木齐、天津及太原等城市位于高烈度区且存在多条活动断裂带，由于地铁工程走向往往取决于城市交通功能的需求，避让原则[1]常常无法实现. 穿越活动断裂带的隧道在断层错动作用下会发生严重破坏，如5.12汶川地震白云顶隧道距映秀端洞口60 m处二衬错台近40 cm，隧道破坏长20 m，局部发生垮塌，路面破损、开裂起拱 [2]. 因此研究断层错动作用对城市地铁隧道结构的影响，是保障隧道工程安全的关键性课题之一.

跨断层隧道结构的影响研究方法包括震后调查、数值模拟及模型试验法，其中震后调查受诸多条件限制，难以大规模开展，因此主要采用后两种方法. 在模型试验方面，Kontogianni等[3]通过模型试验研究了在不同断层倾角的逆断层、走滑断层作用下隧道衬砌管片环向接头部位应力应变规律. Lin等[4]通过提升底板试验模拟了逆冲断层作用下砂土中盾构隧道变形破坏特点. 陈俊岭等[5]通过振动台试验研究了近断层滑冲效应下风力发电塔动力响应，结果表明近断层地震滑冲效应脉冲使得结构响应明显增大. 刘学增等[6]通过模型试验研究了断层倾角为75°的逆断层错动下公路隧道衬砌的反应规律，给出了受拉受压区间以及最大拉压应变发生的位置. 黄强兵等[7]通过地裂缝活动模型试验，研究了西安地铁2号线隧道正交穿越地裂缝带的设防参数. 胡辉[8]进行了相似模型试验，研究了跨断层隧道变形缝设置的有效性. 孙风伯等[9]等推导了隧道最大节段长度并通过实验得以验证. 王明年，崔光耀等[10]展开了断裂黏滑错动作用下模型试验并研究了减震缝与减震层的减震效果，为实际隧道工程的建设提供了参考和依據. 在数值模拟方面，Mohammad等[11]通过数值模拟研究了逆断层错动作用下距错动处不同位置的隧道附近土层剪切带出现的位置及传播路线. 王琼[12]  采用有限元方法研究了断层错动对长大山岭隧道的非线性反应，得到了应力屈服破坏顺序、屈服范围和产生高应力区的位置. 张维庆[13]采用数值模拟研究了穿越断层隧道在断层错动和地震力分别作用下的震害机理. 赵颖等[14]采用数值模拟研究了走滑断层位错作用下城市地铁隧道衬砌结构的损伤发展，重点分析了结构损伤的开始部位、发展过程以及最终的损伤程度，并建立了能够估计结构损伤范围及破坏最严重位置的统计关系式. Shahidi等[15]针对伊朗Koohrang-III输水隧道跨活动断裂的情况，提出断层带部分衬砌采用铰接设计方法，并计算确定了适宜的衬砌节段长度和柔性连接宽度. 赵坤等[16]采用有限元法通过正交试验研究了衬砌节段长度、柔性连接宽度、柔性连接材料强度3种因素对衬砌结构内力的影响. 刘国庆等[17]对跨断层水工隧洞建立三维有限元模型，考虑动接触后在地震作用下围岩出现了明显的错动位移，断层破碎带的存在造成了衬砌位移和应力的增加.

以上研究成果对于穿越活动断层的隧道工程均具有一定的指导意义，但仍然存在一些不足. 目前的研究对象主要集中在深埋山岭隧道，对城市浅埋地铁隧道研究较少，且受条件所限，现有模型试验尺寸都比较小，缺乏对隧道纵向尺度方向的考虑;数值模拟试验结构大多采用弹性本构，结构塑性损伤分析方面研究较少.

本文以乌鲁木齐地铁2号线隧道穿越西山断裂带工程为依托，考虑活动断层错动作用下土层与隧道结构的相互作用，采用ABAQUS建立三维有限元弹塑性模型，研究了隧道二次衬砌塑性应变发展过程、拉压损伤因子及剪切应变的分布，给出了隧道拉裂严重受损、轻微受损及无受损的区间范围，分析了不同错动位移、隧道底部距围岩交界面不同垂直距离时隧道及不同断层破碎带宽度的损伤破坏规律，进行了设置柔性接头的减灾效果研究，以期为实际工程提供参考.

1   工程背景

乌鲁木齐轨道交通2号线马料地站-平川路站区间，由西外环的马料地站起，沿马料地街向西，再向北转到平川路北园春市场平川路站，中间段穿越西山断层，全长约873.1 m. 拟建区间底板埋深17.6～26.5 m，标高831～842.6 m，地质剖面简图见图1. 西山断层形成于中更新世中晚期，最新活动时间为晚更新世晚期，属晚更新世活动断层，走向N45°-75°E，倾向N，倾角44°～83°，具逆冲性质. 穿越断层段前后200 m范围内隧道区间拟采用矿山法施工，其余地段拟采用盾构法施工. 根据乌鲁木齐轨道交通2号线详勘报告，考虑未来100年有突发事件的可能性，建议设防垂直位错hv = 0.165 m.

2   模型建立及参数选取

取断层破碎带前后各100 m范围内建立计算模型，模型长×宽×高为252 m × 68 m × 60 m，隧道埋深为12 m，断层倾角为45°，隧道直径为6.7 m，模型网格划分见图2，隧道结构采用复合式衬砌，纵向划分126份网格，网格纵向间距为2 m. 土体假设为理想弹塑性材料，采用摩尔库伦屈服准则，各土层参数见表1，衬砌结构采用ABAQUS中的塑性损伤模型，依据标准GB 50010—2010《混凝土结构设计规范》[18]的规定，定义衬砌混凝土弹塑性应力应变关系，初期支护与二次衬砌结构参数见表2. 初期支护厚0.35 m，混凝土强度等级为C45;二次衬砌厚0.25 m，混凝土强度等级为C50. 考虑隧道与土体之间的摩擦，其接触面设置罚摩擦，摩擦系数取为0.4[19]. 模拟分析分为三步：初始地应力平衡-隧道开挖-施加断层位移，在前两个步骤中，模型底部及侧边界施加法向约束，上部边界为自由边界，在最后一个步骤中，由于逆断层错动为上盘相对下盘作向上运动，因此为实现逆断层错动过程，释放上盘底部及侧边界法向约束，并在相应的位置上施加位移荷载以模拟断层错动，位移加载边界条件见图3.

图20和图21为隧道沿纵向各个横截面上最大受压和受拉损伤因子的纵向分布曲线图. 当破碎带宽度分别为2 m、26 m和52 m时，最大受压损伤因子分别为0.159 43、0.866 27及0.901 48，最大受拉损伤因子分别为0.902 12、0.933 2及0.954 3. 由图可知，当破碎带宽度由2 m变为26 m时，最大受压损伤因子增加了443.35%，最大受拉损伤因子增加了3.45%;当破碎带宽度由26 m变为52 m时，最大受压损伤因子增加了4.06%，最大受拉损伤因子增加了2.26%. 结合以上分析可知，当断层破碎带宽度由2 m变为26 m时，隧道衬砌破坏显著增加，当破碎带宽度由26 m变为52 m时，衬砌破坏会出现很小程度的增加，基本保持稳定. 产生这种现象的原因依据是泥岩与煤层为Ⅴ级围岩，粉质黏土与圆砾为Ⅵ级围岩，断层错动位置发生在泥岩与煤层接触面处，当破碎带（煤层）宽度较小时，下盘内的粉质黏土发挥主要作用，其围岩性质较软，有利于隧道的变形，衬砌破坏较轻. 随着破碎带宽度的逐渐增加，破碎带围岩对隧道的影响越来越大，二次衬砌破坏逐渐增加，当破碎带宽度达到26 m时，破碎带围岩对隧道的影响基本保持稳定.

5   柔性接头减灾效果研究

断层错动会对地下结构造成严重的破坏，因此研究有效的抗断设计方法十分必要且具实际意义. 目前常用的抗断减灾设计方法主要包括超挖设计、柔性接头、设置减震层及注浆加固4种，该工程采用设置柔性接头的方法. 为确定衬砌节段长度，在设防错动位移作用下二次衬砌的弯矩纵向分布见图22，如图所示，在错动面两侧32 m处隧道二次衬砌形成反弯点，依据Shahidi[15]等人的研究，衬砌节段长度最大不应大于32 m，考虑到混凝土台车的宽度一般为6 ～ 12 m，同时考虑到防水要求，衬砌节段长度不宜设置太短，因此节段衬砌长度最终确定为12 m，设置范围分布在断层破碎带及其附近共96 m，柔性接头长度0.15 m，按不利工况考虑，假定断层错动的位置发生在节段衬砌中点处，模拟中柔性接头被简化为理想弹性材料，密度为1 000 kg/m3，弹性模量为7.8 MPa，泊松比0.47，柔性接头设置示意图见图23.

施加断层错动位移16.5 cm，二次衬砌的塑性应变、拉压损伤因子及剪应变云图见图24. 对比第3节相应的分析可知，二次衬砌塑性应变降低了97.8%，最大受压损伤因子降低了90.7%，最大受拉损伤因子降低了88.8%，最大剪应变降低了98.5%. 值得注意的是，设置柔性接头以后，二次衬砌发生破坏的位置主要集中在跨断层错动面附近的柔性接头处，因此应对此处的结构进行重点加固. 综合以上分析可知，设置柔性接头可以显著降低断层错动变形对结构的损害，但由于柔性接头弹性模量低，较为“柔软”，因此成为了抗断薄弱部位，结构损伤最大的位置也往往出现在柔性接头附近. 以上分析体现出了跨断层地下结构的“以适应地层变形为主，局部结构加固为辅”的设计理念.

6   结   论

1）逆断层错动作用下，错动位移不大于3 cm时，二次衬砌塑性应变为0，结构未发生破坏. 当错动位移进一步加大时，塑性应变首先出现在拱顶，然后是拱底，最后在拱腰处累积. 当错动位移达到16.5 cm时，错动面及其附近的拱腰处产生拉压剪的共同破坏，远离破裂面上盘内拱顶处受拉，拱底处受压，在受损最为严重的拱腰处指向下盘方向的破碎带内拱顶处受压，拱底处受拉.

2）基于拉伸损伤因子计算得到了二次衬砌的混凝土裂缝宽度，并依此为依据将隧道二次衬砌受损区间划分为拉裂严重受损区、轻微受損区与无受损区. 其中，拉裂严重受损区集中在破碎带内，发生在纵向区间105 m至115 m处，总长共10 m;轻微受损区出现在严重受损区两侧，沿破碎带指向上盘为16 m，指向下盘为14 m，总长为30 m.

3）断层错动位移越大，结构受损越严重. 隧道底部距围岩交界面垂直距离越大，土层越厚，错动位移传播过程中消耗能量越多，结构受损越轻. 断层破碎带宽度越大，隧道破坏越严重，当破碎带宽度达到一定值时，破碎带对隧道的影响基本保持稳定.

4）节段衬砌之间设置柔性接头可以显著降低错动位移对结构的危害. 但由于柔性接头较为柔软，是抗断薄弱部位，设置柔性接头以后受损最为严重的区域往往出现在变形缝及其附近，因此应对该处结构进行重点加固.

参考文献

[1]    GB 20909—2014   城市轨道交通结构抗震设计规范[S]. 北京：中国计划出版社，2014：56—57.

GB 50909—2014  Code for seismic design of urban rail transit structures[S]. Beijing：China Planning Press，2014：56—57. （In Chinese）

[2]    孙风伯. 穿越活动断裂山岭隧道抗位错机理与方法研究[D]. 北京：北京交通大学土木建筑工程学院，2018：27—29.

SUN F B. Failure mechanism and method of anti-dislocation of mountain tunnel through active fault[D]. Beijing：School of Civil Engineering，Beijing Jiaotong University，2018：27—29. （In Chinese）

[3]    KONTOGIANNI V A，STIROS S C. Earthquakes and seismic faulting：effects on tunnels[J]. Turkish Journal of Earth Sciences，2003，12（1）：153—156.

[4]    LIN M L，CHUNG C F，JENG F S，et al. The deformation of overburden soil induced by thrust faulting and its impact on underground tunnels[J]. Engineering Geology，2007，92 （3/4）：110—132.

[5]    陈俊岭，阳荣昌，马人乐. 近断层地震滑冲效应下风力发电塔动力响应和振动控制试验研究[J]. 湖南大学学报（自然科学版），2013，40（8）：28—33.

CHEN J L，YANG R C，MA R L. Experimental study of the seismic response and vibration control of wind turbine tower considering near-fault ground motion with fling step[J]. Journal of Hunan University （Natural Sciences），2013，40（8）：28—33. （In Chinese）

[6]    刘学增，林亮伦. 75°倾角逆断层黏滑错动对公路隧道影响的模型试验研究[J]. 岩石力学与工程学报，2011，30（12）：2523—2530.

LIU X Z，LIN L L. Research on model experiment of effect of thrust fault with 75° dip angle stick-slip dislocation on highway tunnel[J]. Chinese Journal of Rock Mechanics and Engineering，2011，30（12）：2523—2530. （In Chinese）

[7]    黄强兵，彭建兵，邓亚虹，等. 西安地铁2号线隧道穿越地裂缝带的设防参数[J]. 岩土力学，2010，31（9）：2882—2888.

HUANG Q B，PENG J B，DENG Y H，et al. Design parameters of Xi'an Metro Line 2 Tunnel passing through active ground fissure zones [J]. Rock and Soil Mechancs，2010，31（9）：2882—2888. （In Chinese）

[8]    胡辉. 穿越活动断层的隧道减震结构研究[D]. 成都：西南交通大学土木工程学院，2013：68—78.

HU H. Study on aseismic structure for tunnel passing through active fault[D]. Chengdu：School of Civil Engineering，Southwest Jiaotong University，2013：68—78. （In Chinese）

[9]    孙风伯，赵伯明，杨清源，等. 穿越活动断层隧道组合抗震缝定量设置的计算公式及试验验证[J]. 中国铁道科学，2018，39（2）：61—70.

SUN F B，ZHAO B M，YANG Q Y，et al. Calculation formula and test verification for quantitative setting of combined seismic joint for tunnel through active fault[J]. China Railway Science，2018，39（2）：61—70. （In Chinese）

[10]  王明年，崔光耀. 高烈度地震区隧道减震模型的建立及其减震效果模型试验研究[J]. 岩土力学，2010，31（6）：1885—1890.

WANG M N，CUI G Y. Establishment of tunnel damping model and research on damping effect with model test in highly seismic area[J]. Rock and Soil Mechancs，2010，31（6）：1885—1890. （In Chinese）

[11]  MOHAMMAD H B，ALI N，RONAK M，et al. Evaluation of underground tunnel response to reverse fault rupture using numerical approach[J]. Soil Dynamics and Earthquake Engineering，2016，83：1—17.

[12]  王瓊. 跨断层隧道地震反应分析[D]. 哈尔滨：中国地震局工程力学研究所，2011：43—70.

WANG Q. Seismic response analysis of crossing-fault tunnel[D]. Harbin：Institute of Engineering Mechanics，China Earthquake Administration，2011：43—70. （In Chinese）

[13]  张维庆. 穿越断层隧道震害机理以及抗减震技术研究[D]. 成都：西南交通大学土木工程学院，2012：61—71.

ZHANG W Q. Study on the seismic damage mechanism of tunnels through fault and the anti-seismic and damping measures[D]. Chengdu：School of Civil Engineering，Southwest Jiaotong University，2012：61—71. （In Chinese）

[14]  赵颖，郭恩栋，刘智，等. 走滑断层位错作用下城市地铁隧道损伤分析[J]. 岩土力学，2014，35（S2）：467—473.

ZHAO Y，GUO E D，LIU Z，et al. Damage analysis of urban metro tunnel under strike-slip fault[J]. Rock and Soil Mechancs，2014，35（S2）：467—473. （In Chinese）

[15]  SHAHIDI A R，VAFAEIAN M. Analysis of longitudinal profile of the tunnels in the active faulted zone and designing the flexible lining[J]. Tunnelling and Underground Space Technology，2005，20（3）：213—221.

[16]  赵坤，陈卫忠，赵武胜，等. 逆断层错动作用下隧道衬砌铰接设计参数研究[J]. 岩石力学与工程学报，2018，37（S1）：3412—3421.

ZHAO K，CHEN W Z，ZHAO W S，et al. Study on parameters of articulated design of tunnel lining under reverse fault dislocation[J]. Chinese Journal of Rock Mechanics and Engineering，2018，37（S1）：3412—3421. （In Chinese）

[17]  刘国庆，肖明，杨阳，等. 跨断层水工隧洞地震响应数值模拟分析方法[J]. 湖南大学学报（自然科学版），2018，45（11）：141—148.

LIU G Q，XIAO M，YANG Y，et al. Numerical simulation and analysis method of seismic response for hydraulic tunnel across fault[J]. Journal of Hunan University （Natural Sciences），2018，45（11）：141—148. （In Chinese）

[18]  GB 50010—2010   混凝土结构设计规范[S]. 北京：中国建筑工业出版社，2015：19—20.

GB 50010—2010   Code for design of concrete structures[S]. Beijing： China Architecture & Building Press， 2015：19—20. （In Chinese）

[19]  HUO H. Seismic design and analysis of rectangular underground structures[D]. West Lafayette：School of Civil Engineering，Purdue University，2005：23—28.

[20]  施士昇. 混凝土的抗剪强度、剪切模量和弹性模量[J].土木工程学报，1994，32（2）：47—51.

SHI S S. Shear strength，modulus of rigidity and Young's modulus of concrete[J]. China Civil Engineering Journal. 1994，32（2）：47—51. （In Chinese）

[21]  PAUL Z. Torsional strength of prestressed concrete member[J]. ACI Journal，1961，32（10）：1337—1359.

[22]  CHEN G M，TENG J G，ASCE M，et al. Finite-element modeling of intermediate crack debonding in FRP-Plated RC beams[J]. Journal of Composites for Construction. 2011，15（3）：339—353.

[23] ROTS J G. Computational modeling of concrete fracture[D]. Netherlands：School of Civil Engineering，Delft University of Technology，1988：49—50.

[24]  SL 191—2008   水工混凝土結构设计规范[S]. 北京：中国水利水电出版社，2008：20—21.

SL 191—2008   Design code for hydraulic concrete structures[S]. Beijing：China Water & Power Press，2008：20—21. （In Chinese）