趙伯明，趙天次，周玉書

(1.北京交通大學(xué) 城市地下工程教育部重點實驗室, 北京 100044；2.北京交通大學(xué) 土木建筑工程學(xué)院, 北京 100044)

活斷層在我國分布廣泛，突發(fā)性的斷層錯動是產(chǎn)生地震的主要根源，長期以來給人類的生命財產(chǎn)和結(jié)構(gòu)工程的安全性帶來了很大威脅。大量研究表明，活斷層錯動導(dǎo)致的同震地表破裂、土體變形、強地面震動及其他次生地質(zhì)災(zāi)害均會給穿越活斷層隧道帶來危害[1-11]。然而在大型工程的選線和設(shè)計中不可避免的會穿越或臨近活斷層，因此跨活斷層隧道的抗位錯研究具有重要的理論意義和工程價值。

目前針對穿越活斷層隧道抗位錯分析國內(nèi)外開展了諸多研究，研究方法包括現(xiàn)場調(diào)研、模型試驗、解析計算和數(shù)值模擬等。

現(xiàn)場調(diào)研可以直接得到最真實的一手資料，通過震后短時間內(nèi)的儀器監(jiān)測可得到斷層附近的地表同震位移[12]，運用回歸分析手段將實地調(diào)研數(shù)據(jù)進(jìn)行處理，建立地表同震位移與震級、上覆土層厚度等指標(biāo)的定量關(guān)系[13-17]。但該方法受限于觀測點的數(shù)量和密度，不同區(qū)域的地質(zhì)構(gòu)造和土層性質(zhì)差異也會導(dǎo)致回歸關(guān)系式的差異[18]，且大多數(shù)研究重點關(guān)注了地表的同震位移，對隧道等地下工程而言，土體內(nèi)部的變形對結(jié)構(gòu)受力影響更大，斷層錯動導(dǎo)致的圍巖狀態(tài)變化同樣值得關(guān)注。

部分學(xué)者針對斷層錯動對隧道結(jié)構(gòu)的影響進(jìn)行了模型實驗研究[19-24]，但在實驗中由于模型尺寸的限制，對斷層錯動的模擬通常簡化為單側(cè)底板的抬升，并不能模擬真實的土體變形情況。

在基于彈性半空間的解析解法方面，學(xué)者們[25-27]應(yīng)用簡化的斷層模型和位錯理論，給出了半無限空間中矩形滑動斷層的地震位移場的解析解。該方法推導(dǎo)清晰，操作便捷，其中應(yīng)用最為廣泛的是Okada提出的均勻半無限空間解析解法，但其缺點是對于多層介質(zhì)模型或三維復(fù)雜情況并不適用。

隨著地質(zhì)資料和勘測手段的進(jìn)步，Okada的解析解法已經(jīng)略顯粗糙，在均勻半無限空間解析解基礎(chǔ)上進(jìn)一步發(fā)展而來的層狀半空間中的靜態(tài)彈性變形計算方法考慮了地球的層狀構(gòu)造，位移矢量需要在柱坐標(biāo)系下利用柱諧函數(shù)來表達(dá)，并通過傳播矩陣的算法求解[28-29]，近年來已經(jīng)逐漸成為主流趨勢。

對于地球?qū)訝顦?gòu)造對同震變形的影響程度，文獻(xiàn)[30-31]的研究表明：在地殼內(nèi)發(fā)生的地震，在震源距不大于100 km的情況下，層狀構(gòu)造的影響可達(dá)到20%。

因此，隨著計算理論的發(fā)展，人們對斷層位錯發(fā)生后土體的變形了解逐步加深，有限元模型中的輸入也應(yīng)盡量接近真實的構(gòu)造環(huán)境，以提高計算結(jié)果的可靠性，更加精細(xì)化的位錯輸入方法有待提出。

本文首先對現(xiàn)有的跨斷層隧道抗位錯計算方法進(jìn)行誤差分析，宏觀上分析了土體層狀構(gòu)造和斷層破裂面上位錯分布導(dǎo)致的同震變形場差異，微觀上分析了不同輸入方法導(dǎo)致的土體變形差異，揭示了現(xiàn)有隧道抗位錯計算中同震變形輸入方法的不足。然后基于正交化格林函數(shù)方法編制程序，對層狀半無限空間中的靜態(tài)彈性變形進(jìn)行計算，結(jié)合Abaqus的input文件將斷層錯動導(dǎo)致的土體變形場以位移邊界的形式施加到有限元模型，并采用Matlab建立了輸入流程的GUI界面，使建模流程大為簡化。最后將該方法應(yīng)用于成蘭鐵路某隧道的三維抗位錯分析中，實現(xiàn)了隧道工程抗位錯分析中斷層錯動的精細(xì)化輸入。

1 現(xiàn)有跨斷層隧道抗位錯機理研究

我國西部地區(qū)地質(zhì)條件復(fù)雜，活斷層分布密集，隨著我國西部地區(qū)鐵路隧道逐漸增多，部分隧道不可避免的需要穿越眾多的活斷層。當(dāng)隧道工程處于斷層附近時往往需要進(jìn)行抗位錯分析，以得到隧道結(jié)構(gòu)在斷層錯動作用下變形效應(yīng)和破壞機理，隧道工程與活斷層之間的的空間關(guān)系見圖1，黑色單箭頭代表斷層的錯動。

圖1 跨斷層隧道空間位置關(guān)系

在隧道工程的抗位錯分析中，通常截取有限區(qū)域為計算區(qū)域，并設(shè)置人工邊界來模擬活斷層的位錯[10,32]，但現(xiàn)有的建模過程中往往存在著過度的簡化，簡化分析模型見圖2。

圖2 現(xiàn)有研究中的簡化模型[32]

由圖2可知，上述簡化建模方法中將斷層破裂面直接置于隧道中部，通過斷層上盤底面的整體抬升模擬逆斷層的錯動。而實際地層環(huán)境中大部分?jǐn)鄬拥钠屏衙娌⑽粗苯咏財嗨淼?，對于埋深不為零的隱伏斷層或破裂面未直達(dá)地表的斷層而言，該模擬方法顯然是不適用的。此外，上述研究中將斷層破裂面上的位錯量設(shè)置為均勻分布，與實際斷層面上的位錯分布存在較大的差別。通過分析可知建模過程中的誤差來源如下：①土體的層狀結(jié)構(gòu)；②破裂面上的位錯分布；③斷層錯動的輸入方法。本節(jié)以典型地震為例，通過對比分析得到各誤差源的影響大小，說明了建立精細(xì)化位錯輸入方法的必要性。

1.1 土體層狀結(jié)構(gòu)對地表同震變形的影響

土體的層狀構(gòu)造在同震變形場的計算中存在影響，為量化這一誤差，本節(jié)選取典型地震2013年Mw6.7級蘆山地震作為算例，對成層土體模型和均勻土體模型下的地表同震變形進(jìn)行對比分析，研究土體層狀結(jié)構(gòu)對同震變形場的影響。

由王衛(wèi)民等[33]給出的蘆山地震反演結(jié)果可知，發(fā)震斷層為一條埋深2 km的隱伏逆斷層，斷層破裂面并未延伸到地表，斷層參數(shù)見表1。斷層參數(shù)及坐標(biāo)系如圖3所示。其中斷層長度為L，寬度為W，傾角為β，斷層面上位錯滑動角為α，滑動量為u。

圖3 斷層空間坐標(biāo)系

表1 蘆山地震震源參數(shù)

成層土體模型的土層劃分參考參考蘆山地震震源區(qū)地殼地震波速度模型[34]及前人的研究成果[35-41]，綜合給出了蘆山震區(qū)的成層土體模型，自地表向下共劃分了4個不均勻波速層，土層參數(shù)見表2。均勻土體模型與成層土體模型的最上層土采用相同的力學(xué)參數(shù)。

表2 蘆山地震震區(qū)地殼結(jié)構(gòu)模型

首先計算均勻土體模型中蘆山地震導(dǎo)致的地表同震變形，豎向位移等值線和水平位移矢量分別見圖4，其中虛線矩形框表示斷層的地表投影。

圖4 均勻位錯模型地表豎向、水平位移矢量圖

地表位移的分布受斷層產(chǎn)狀影響，豎向位移的數(shù)值總體上大于水平位移。由圖4(a)可知，地表隆起主要分布在斷層上盤區(qū)域即Y軸正半軸，沉降集中下斷層下盤區(qū)域，同震變形在斷層上斷線的地表投影附近變化最為劇烈。豎向位移的隆起峰值為0.802 m，沉降峰值為0.062 m，這是由于斷層傾角較小，斷層面上的位錯量經(jīng)過土層傳導(dǎo)出現(xiàn)了一定的偏移。

由圖4(b)可知，上盤和下盤的土體均向斷層上斷線的地表投影擠壓，隨著與上斷線地表投影線距離的減小，水平位移值逐漸增大，峰值為0.338 m。

基于分層土體模型的蘆山地震地表同震變形計算結(jié)果見圖5，從結(jié)果來看，地表豎向位移和水平位移的分布和數(shù)值都與GPS監(jiān)測數(shù)據(jù)[42-43]給出的結(jié)果相近。

圖5 分層位錯模型地表豎向、水平位移矢量圖

對比圖5兩種模型的地表同震變形計算結(jié)果，考慮地殼層狀構(gòu)造之后地表的豎向位移和水平位移在分布規(guī)律上絕大部分都保持一致，豎向位移隆起峰值為0.775 m，沉降峰值為0.046 m，水平位移峰值為0.395 m，但是在峰值位置上有所差異。均勻土體模型的隆起峰值和沉降峰值都要大于分層土體模型，但水平位移峰值小于分層土體模型，這一差異是由均勻土體模型較分層土體模型更軟的土體力學(xué)參數(shù)導(dǎo)致的。為清晰的展示兩者差異，將兩種模型的計算結(jié)果取差值，見圖6。

圖6 分層模型和均勻模型的地表豎向、水平位移差值

由圖6(a)可知，豎向位移的差異主要集中在上斷線投影附近，地表隆起存在一條沿斷層走向的0.1 m以上分布帶，地表沉降則對稱分布在上斷線投影，量值上僅0.05 m左右。這說明均勻模型具有更大的隆起和沉降峰值，豎向位移最大差值為0.308 m，約占分層土體模型豎向位移最大值的39.7%。

由圖6(b)可知，水平位移的差異同樣體現(xiàn)在上斷線投影處，差異帶寬度10 km左右，最大差值0.236 m，達(dá)到分層土體模型水平位移最大值的59.8%。

由蘆山地震震例的計算結(jié)果可見，土體的分層結(jié)構(gòu)對地表同震位移的影響是不可忽略的，位移差異達(dá)40%～60%，尤其是在靠近斷層上斷線地表投影區(qū)域，需要考慮地震發(fā)生地區(qū)地殼層狀構(gòu)造的影響。

1.2 破裂面上位錯量分布對地表同震變形的影響

在真實的地震中，斷層破裂面的幾何形狀與位錯量分布都很復(fù)雜。在以往的跨斷層隧道計算中，往往假設(shè)斷層面位錯是均勻分布的，但實際斷層面上的位錯分布是非一致的。吳忠良[44]通過對多個地震斷層面上的滑動量調(diào)查統(tǒng)計發(fā)現(xiàn)斷層面上的位錯分布是高度不均勻的；Banerjee等[45]利用41個連續(xù)觀測的遠(yuǎn)場GPS資料對Sumatra-Andaman地震斷層的滑動量進(jìn)行計算，發(fā)現(xiàn)斷層上的最大滑動量出現(xiàn)在斷層中部附近。

在經(jīng)典半無限空間的位錯理論中，Okada均勻土體模型的有限矩形位錯源公式是由點源公式在斷層面上積分所得，分層土體模型中有限斷層源是由離散的點源格林函數(shù)線性疊加而得，兩種模型均可對斷層進(jìn)行離散化，使破裂面上的位錯量接近真實的位錯分布，再求得對應(yīng)的同震變形場。

為分析斷層破裂面上位錯量分布對同震變形場的影響，選取典型地震1989年美國Mw6.95級Loma地震分別計算位錯量均勻分布和真實分布兩種情況的地表同震變形，Loma地震發(fā)震斷層參數(shù)見表3，Loma地區(qū)的土體參數(shù)見表4[47]。

表3 Loma地震震源參數(shù)

表4 Loma地震震區(qū)地殼結(jié)構(gòu)模型

沿走向和傾向?qū)l(fā)震斷層劃分為2 km×2 km的子斷層，每個子斷層上的位錯量反演結(jié)果見圖7[46]。

圖7 Loma地震斷層破裂面位錯量分布(單位：cm)

圖7中每個子斷層上位錯量的精確反演結(jié)果反映了斷層面上位錯的真實分布情況，再求取各個子斷層面上位錯量的平均值作為均勻分布模型的位錯量，采用分層土體模型分別計算Loma地震兩種位錯分布模式下的地表同震位移。

真實位錯分布模式下地表三個方向的位移等值線分別見圖8。

圖8 真實位錯分布下豎向、X向、Y向地表位移等值線(單位：m)

由圖8可知，從結(jié)果來看地表位移分布符合斜滑斷層的運動機制。在斷層傾向抬升作用和走向擠壓作用的疊加影響下，豎向位移有整體向左移動的趨勢，豎向位移峰值在斷層中部左側(cè)約7 km的位置，峰值量約為0.7 m。X方向位移主要受走滑分量的控制，表現(xiàn)為上盤區(qū)域沿走向向左移動，下盤區(qū)域沿走向向右移動，X向位移峰值在斷層中部左側(cè)約12 km的位置，峰值量約為0.4 m。Y方向位移受到走滑分量和傾滑分量的共同作用，表現(xiàn)為大部分區(qū)域向Y正方向移動，只有斷層右側(cè)小部分區(qū)域向Y負(fù)方向移動，Y向位移峰值在斷層中部左側(cè)約8 km的位置，峰值量約為0.3 m。

為了量化位錯分布模式對地表同震位移的影響，將均勻分布模式下的結(jié)果與真實分布模式的結(jié)果取差值進(jìn)行分析，三個方向上的位移差值分布分別見圖9。由圖9可知，均勻分布下的三方向位移在整個斷層區(qū)域內(nèi)都與真實分布有較大的差別，誤差最大的位置對應(yīng)位移峰值所在位置。

圖9 位錯均勻分布與真實分布的豎向、X向、Y向地表位移差(單位：m)

考察斷層走向中線剖面位置處的地表同震位移見圖10。由圖10可知，三個方向的位移值都隨著遠(yuǎn)離斷層而逐漸減小，隨著與斷層距離的減小，均勻分布模式會產(chǎn)生較大的誤差，具體數(shù)值見表5。

圖10 斷層走向中線剖面的地表同震位移對比

表5 不同位錯分布模式的位移峰值誤差

由表5可見，均勻分布模式在水平和豎向都會給地表同震位移帶來誤差，水平同震位移峰值誤差在X方向達(dá)到36.6%，在Y方向達(dá)到8.4%；豎向同震位移峰值誤差在Z方向達(dá)到9.6%，這些誤差放大了斷層錯動的破壞效應(yīng)。

綜上所述，將斷層面上的位錯量按均勻分布考慮會給計算結(jié)果帶來誤差，在計算土體的同震位移場時，采用反演得到的真實位錯分布具有一定的必要性。

1.3 斷層錯動輸入方法對地表同震變形的影響

在求得斷層錯動導(dǎo)致的土體同震變形場之后，將該同震變形場準(zhǔn)確合理的輸入截取的計算區(qū)域是保證跨斷層隧道抗位錯分析準(zhǔn)確性的重要一環(huán)。現(xiàn)有研究中通常采用斷層上盤底面的整體抬升模擬逆斷層的錯動，見圖2，這是一種過度簡化的手段，顯然會帶來一定的誤差。本小節(jié)選取矩震級Mw7.0級典型地震對斷層錯動導(dǎo)致的模型底面同震變形進(jìn)行計算，并將模型底面的真實變形與簡化方法的輸入值進(jìn)行對比，量化了簡化方法的誤差，說明了選擇合理斷層錯動輸入方法的必要性。

為使對比結(jié)果更為清晰，忽略土體分層和位錯分布的影響，按均勻土體模型且位錯均勻分布計算，土層參數(shù)設(shè)置為表2中的第一層土，選取逆斷層進(jìn)行分析，由斷層破裂參數(shù)和矩震級Mw之間的統(tǒng)計公式[46,48-49]得Mw7.0級典型地震的震源參數(shù)見表6。

表6 Mw7.0級典型地震震源參數(shù)

選取斷層上斷線中軸面為觀察剖面，斷層模型及觀測剖面位置如圖11，分別計算深度為0、50、100、150、200 m處觀察剖面的土體水平和豎向同震變形值，將各深度真實的土體同震變形值和簡化輸入變形值對比，見圖12。

圖11 斷層模型及觀測剖面

由圖12(a)可知，豎直方向土體的真實同震變形與簡化方法中輸入的土體變形值趨勢一致，隨著深度的增加簡化值與真實值的差值逐漸減小，上盤區(qū)域差值約0.2 m.下盤區(qū)域差值約0.25 m，誤差為真實變形值的25%左右；由圖12(b)可知，而水平方向土體的真實同震變形與簡化輸入變形值差異顯著，簡化方法中認(rèn)為斷層上盤固定不動僅下盤抬升的輸入方式顯然是不準(zhǔn)確的，上盤和下盤的輸入差值都為0.4 m左右，水平變形輸入誤差達(dá)80%左右。

2 跨斷層隧道抗位錯精細(xì)化計算

由上節(jié)分析可知，土體層狀結(jié)構(gòu)、破裂面上的位錯量分布及錯動輸入方法均會導(dǎo)致地表同震變形的計算出現(xiàn)誤差，其中影響最為顯著的是錯動輸入方法，誤差達(dá)25%～80%；土體層狀結(jié)構(gòu)影響次之，誤差為40%～60%；破裂面上的位錯量分布影響最小，僅10%左右，該誤差在位錯量反演不明確的條件下可忽略。

因此，本節(jié)提出了一套精細(xì)化的位錯輸入方法，首先由斷層位錯出發(fā)，求得斷層錯動導(dǎo)致的模型邊界上的同震變形場，再將該同震變形場以位移人工邊界的形式施加至有限元模型的邊界上，以實現(xiàn)斷層錯動的精確輸入。以成蘭鐵路跨越龍門山破裂帶某隧道為例，探究隧道在斷層錯動作用下的破壞機理及力學(xué)特性，為隧道工程抗位錯措施的建立提供指導(dǎo)。

2.1 精細(xì)化位錯輸入方法的實現(xiàn)過程

早期研究中將斷層上盤或下盤的移動簡化為模型單側(cè)底面的抬升，這是較不合理的一種簡化輸入方法。本文提出的精細(xì)化斷層錯動輸入方法建模思路見圖13。

圖13 精細(xì)化輸入方法建模思路

計算流程如下：首先輸入斷層參數(shù)和土層參數(shù)，求得斷層錯動導(dǎo)致的有限元模型區(qū)域的同震變形場，再對計算結(jié)果進(jìn)行提取，獲得模型邊界處的變形值并寫入input文件，最后通過input文件將同震變形場施加到Abaqus有限元模型，以實現(xiàn)斷層錯動的精確輸入，自編程序計算流程見圖14。

圖14 計算流程圖

2.2 跨活斷層隧道數(shù)值計算模型

以我國西部成蘭鐵路穿越北川—映秀活斷層某隧道為例進(jìn)行抗位錯分析，北川—映秀是龍門山中央斷裂帶的主斷裂，也是汶川地震的發(fā)震構(gòu)造之一，本文將其簡化為逆沖斷層，斷層參數(shù)見表6。

由圖12可知，斷層位錯在土體產(chǎn)生的位錯分布主要集中在上斷線投影前后500 m范圍內(nèi)，在這個范圍外的土體位移值雖然很大，但是位移分布一致，位移差值很小，隧道結(jié)構(gòu)的破壞主要是由土體位移差值造成的，所以有限元模型的規(guī)模取縱向長為500 m?？紤]圍巖影響范圍為3～5倍洞徑，模型寬和高分別取為100 m、100 m，Abaqus模型見圖15，坐標(biāo)系原點為圖中模型左下角點，模型整體位于第一象限。

圖15 Abaqus有限元模型

為了驗證數(shù)值模型的準(zhǔn)確性，首先對不加隧道結(jié)構(gòu)只有土體的數(shù)值模型進(jìn)行計算，對該模型施加斷層位錯作用產(chǎn)生的土體同震變形場，提取地表沿縱向中線的水平和豎向位移值，與Okada解析解進(jìn)行對比，結(jié)果見圖16。由圖16可知，從計算結(jié)果可以看出，本文精細(xì)化輸入方法的數(shù)值解與使用Okada位錯理論計算得到的解析解相差極小，驗證了本文方法的準(zhǔn)確性。

圖16 有限元模型計算結(jié)果驗證

下一步加入隧道結(jié)構(gòu)形成土體-隧道有限元模型，隧道頂板覆土厚度為20 m，隧道斷面選取圓形截面，尺寸參數(shù)見圖17，模型中各材料的參數(shù)見表7。由于本文重點為為計算方法展示，為節(jié)約計算時間隧道襯砌結(jié)構(gòu)采用彈性本構(gòu)，若有更高的計算需求可替換為各類彈塑性本構(gòu)，本文方法同樣適用。土體采用Abaqus中的實體單元C3D20R模擬，隧道襯砌采用殼單元S4R模擬，土體-隧道有限元模型網(wǎng)格見圖18，模型X軸為隧道橫向?qū)挾确较?，Y軸為隧道縱向長度方向，Z軸為土體的豎向高度方向。

圖17 隧道結(jié)構(gòu)斷面尺寸參數(shù)(單位：cm)

表7 有限元模型材料參數(shù)

圖18 有限元模型網(wǎng)格劃分

為計算斷層錯動而引發(fā)的隧道襯砌結(jié)構(gòu)的受力與變形效應(yīng)，將上述基于分層土體位錯模型計算得到的土體同震變形場作為模型的位移邊界施加在模型上來模擬斷層的錯動作用，得到隧道襯砌結(jié)構(gòu)的破壞發(fā)展過程，斷層錯動前后的隧道變形絕對值云圖見圖19。

圖19 斷層錯動前后的隧道變形絕對值

由圖19可知，變形絕對值沿隧道縱向先減小后增大，隧道中部變形值的變化速率最大，為進(jìn)一步分析隧道襯砌的受力分布情況，選取典型觀測位置見圖20，繪制隧道各處主應(yīng)力隨隧道縱向距離變化曲線見圖21。

圖20 隧道橫截面觀測位置分布

圖21 各觀測位置主應(yīng)力隨隧道縱向距離變化

由圖21可知，隧道的拱腰和拱腳各處整體呈受壓狀態(tài)，最大壓應(yīng)力出現(xiàn)在拱腰附近，拱腳處也具有較大的壓應(yīng)力峰值；而拱頂和仰拱在斷層上盤靠近上斷線投影處呈受拉狀態(tài)，這與斷層上盤在該處的大幅抬升密切相關(guān)，最大拉應(yīng)力出現(xiàn)在仰拱附近。由汶川地震的震害調(diào)研可知[24,50-51]，隧道的仰拱和拱頂主要出現(xiàn)抗拉破壞，拱腰和拱腳主要出現(xiàn)抗壓破壞，拱肩損傷較小，這與計算結(jié)果相符。因此，本文基于分層土體模型提出的精細(xì)化斷層錯動輸入方法可為工程中隧道的抗位錯分析提供一定的指導(dǎo)作用。

3 結(jié)論

通過理論分析和數(shù)值計算，探究了土體層狀構(gòu)造、斷層破裂面上位錯分布和錯動輸入方法在同震變形場計算中帶來的誤差，揭示了現(xiàn)有隧道抗位錯計算中同震變形輸入方法的不足，提出了精細(xì)化的位錯輸入方法，并將該方法應(yīng)用于成蘭鐵路某隧道的三維抗位錯分析，得出以下結(jié)論：

(1)在斷層錯動導(dǎo)致土體同震變形場的計算中，土體分層結(jié)構(gòu)的影響是不可忽略的，地表同震位移差異達(dá)40%～60%。尤其是在靠近斷層上斷線地表投影區(qū)域，位移差異更為顯著，需要考慮地震發(fā)生地區(qū)地殼層狀構(gòu)造的影響。

(2)斷層面上的位錯量的分布模式同樣給同震變形場的計算結(jié)果帶來誤差，但量值僅為10%左右，該誤差在震源位錯量反演不明確的條件下可忽略。

(3)錯動輸入方法導(dǎo)致的地表同震變形誤差最大，傳統(tǒng)簡化方法的誤差可達(dá)25%～80%，將斷層上盤固定不動僅下盤抬升的輸入方式顯然是不準(zhǔn)確的。

(4)將本文方法的數(shù)值解與Okada解析解對比，驗證了本文精細(xì)化輸入方法的準(zhǔn)確性，并以成蘭鐵路某隧道為例分析了其在典型斷層錯動下的破壞狀態(tài)，隧道的仰拱和拱頂主要出現(xiàn)抗拉破壞，拱腰和拱腳主要出現(xiàn)抗壓破壞，拱肩損傷較小，可為跨斷層隧道的抗位錯分析提供參考。