田朝陽 蘭恒星② 張 寧 許博聞

(①長安大學(xué)地質(zhì)工程與測繪學(xué)院， 西安 710054， 中國)

(②中國科學(xué)院地理科學(xué)與資源研究所， 資源與環(huán)境信息系統(tǒng)國家重點實驗室， 北京 100101， 中國)

0 引 言

CZ鐵路是我國藏區(qū)綜合交通運輸體系的重要組成部分，對維護我國國家長治久安、促進西藏自治區(qū)經(jīng)濟社會發(fā)展具有重大意義(彭建兵等， 2020)。CZ鐵路沿線地貌陡峻，內(nèi)外地質(zhì)作用異常活躍，多場環(huán)境耦合作用強烈，地質(zhì)環(huán)境條件極其復(fù)雜。全線橋隧占比超過80%，存在多條特長深埋隧道，在極端復(fù)雜的工程環(huán)境及強烈的工程擾動影響下，這些深埋長隧建設(shè)過程中工程災(zāi)害問題更加凸顯(郭長寶等， 2017； 薛翊國等， 2020； 趙勇等， 2021)。其中：色季拉山隧道位于西藏自治區(qū)林芝市境內(nèi)，全長37.965km，最大埋深1696m，為新建CZ鐵路雅安林芝段先期開工隧道和全線重點控制性工程。隧道地處青藏高原東南緣的喜馬拉雅東構(gòu)造節(jié)附近，屬于印度洋板塊向歐亞板塊俯沖、楔入的最前緣，是現(xiàn)今全球構(gòu)造活動最為劇烈的地區(qū)之一(常利軍等， 2015； 蘭恒星等， 2021)。強烈構(gòu)造背景下形成的構(gòu)造應(yīng)力和深埋應(yīng)力耦合成的高地應(yīng)力特征，加之隧道沿線高度發(fā)育花崗巖、閃長巖和片麻巖等硬質(zhì)巖，使得色季拉山隧道在建設(shè)過程中，將面臨著嚴重的高地應(yīng)力巖爆問題。

巖爆是在工程開挖過程中，圍巖發(fā)生應(yīng)力釋放而產(chǎn)生脆性破壞的一種巖體失穩(wěn)過程，具有突發(fā)性、高危害性等特點，是高地應(yīng)力地區(qū)地下隧道建設(shè)過程中所面對的主要風(fēng)險之一(靳曉光等， 2013； 蔣鈺峰等，2019)。巖爆的發(fā)生取決于兩個條件：一是巖體能夠儲存大量能量，同時在破壞過程中表現(xiàn)出強烈的彈射傾向； 二是巖體所在區(qū)域存在形成較高的應(yīng)力集中的地應(yīng)力環(huán)境(王開洋等， 2014； Pan et al.，2021)。

準確的預(yù)測預(yù)判巖爆發(fā)生的位置以及強烈程度對高地應(yīng)力區(qū)域地下工程的修建具有重要指導(dǎo)意義，也是巖爆研究中的熱點和難點問題。Neyma于1972年提出了能夠反映巖石彈性變形能儲存能力的巖爆傾向性指數(shù)Wet(Homand et al.，1990)，其后諸多學(xué)者采用此指標進行巖爆傾向性評估，并根據(jù)不同的適用條件、提出了相應(yīng)的巖爆傾向性評價標準(Kidybinski， 1981； Singh， 1988)。然而，除了圍巖的自身性質(zhì)之外，巖爆的發(fā)生還取決于圍巖所處區(qū)域的地應(yīng)力環(huán)境(李鵬翔等， 2019)。Barton(1974)研究挪威隧道圍巖穩(wěn)定性過程中發(fā)現(xiàn)，巖石單軸抗壓強度與最大主應(yīng)力的比值與巖爆之間具有一定的關(guān)系，并根據(jù)強度應(yīng)力比指標劃分了相應(yīng)的巖爆等級分級標準。此后，國內(nèi)外學(xué)者根據(jù)工程實踐經(jīng)驗以及理論推導(dǎo)，先后提出了Russense判據(jù)(Russense,1974)、Hoek判據(jù)(Hoek et al., 1980)、陶振宇判據(jù)(張鏡劍等， 2008)、二郎山判據(jù)(徐林生等， 1999)、國標法(中華人民共和國國家標準編寫組，2008)等一系列應(yīng)力指標型(應(yīng)力與強度或者強度與應(yīng)力的比值)巖爆判據(jù)。隨著新技術(shù)理論的發(fā)展，以及對巖爆問題探索的深入，涌現(xiàn)了一批如模糊綜合預(yù)測、BP神經(jīng)網(wǎng)絡(luò)預(yù)測、聲發(fā)射現(xiàn)場檢測預(yù)報等新的巖爆預(yù)測方法與理論(李庶林， 2009； 陳鵬宇等， 2014； Roohollah et al.，2019)。

上述方法中，應(yīng)力指標判據(jù)法同時考慮了巖爆發(fā)生的兩大條件：圍巖自身性質(zhì)以及外部地應(yīng)力條件，并分別將其用地應(yīng)力量值、巖石強度這兩個相對易獲取的指標參數(shù)進行表征。在實際應(yīng)用中，通過有限的實測地應(yīng)力數(shù)據(jù)進行反演，結(jié)合地質(zhì)資料中巖石的強度參數(shù)，即可方便快捷地得到隧道的巖爆預(yù)測結(jié)果，是目前應(yīng)用最廣泛的巖爆預(yù)測方法。然而，現(xiàn)有強度應(yīng)力判據(jù)采用巖石強度參數(shù)多通過完整巖塊的室內(nèi)試驗獲取，忽視了原位巖體強度的結(jié)構(gòu)損傷效應(yīng)。此外，這些判據(jù)多為巖爆等級的定性預(yù)測，難以給出巖爆的范圍和程度的定量預(yù)測。而量化研究巖爆對于硐室圍巖的破壞程度和影響深度，對于選取恰當?shù)募庸谭绞揭约爸ёo參數(shù)，預(yù)防并降低巖爆對于隧道建設(shè)的不利影響、確保工程的順利開展具有重要意義。

為此，本文收集了大量國內(nèi)外典型巖爆隧道的特征參數(shù)并進行統(tǒng)計分析，建立了一種通過應(yīng)力強度比來預(yù)測巖爆爆坑深度的巖爆定量化預(yù)測模型。同時，通過室內(nèi)力學(xué)試驗和巖體強度損傷理論獲取了隧道沿線巖體強度參數(shù)，以多元線性回歸法反演了隧道沿線地應(yīng)力參數(shù)，最終對CZ鐵路色季拉山隧道的巖爆風(fēng)險進行了定量化預(yù)測和分析，并與同屬雅魯藏布縫合帶影響區(qū)、地質(zhì)條件相近的已貫通CZ鐵路巴玉隧道巖爆風(fēng)險進行了對比評價，成果可以為色季拉山隧道施工與巖爆風(fēng)險防控提供重要參考。

1 巖爆定量化預(yù)測模型與巖體強度損傷準則

1.1 脆性破壞定量化預(yù)測模型

巖爆是多種因素共同耦合的結(jié)果，其中應(yīng)力和圍巖性質(zhì)(主要包括強度和完整性)是導(dǎo)致脆性破壞發(fā)生的最顯著的兩個條件(Martin et al.，1994)。對于裂隙網(wǎng)絡(luò)不發(fā)育的完整巖體，當應(yīng)力強度比較高時，開挖界面附近的巖體發(fā)生應(yīng)力集中而產(chǎn)生脆性破壞，且隧道圍巖的破壞深度隨著應(yīng)力強度比的增加而逐漸增加，即硐室的破壞深度與其應(yīng)力強度比之間具有一定的函數(shù)關(guān)系(Martin et al.，1999)?；谶@一發(fā)現(xiàn)，Martin et al. (1999)根據(jù)實測數(shù)據(jù)對應(yīng)力強度比和硐室圍巖的破壞深度進行了擬合(圖1)，建立了應(yīng)力強度比與硐室破壞深度的關(guān)系表達式：

圖1 硐室圍巖破壞深度與最大切應(yīng)力關(guān)系圖(Martin et al.， 1999)Fig. 1 Relationship between depth of failure and the maximum tangential stress at the boundary of the opening

(1)

式中：Df為硐室破壞影響深度，等于隧道破壞深度與隧道半徑的和；a為隧道半徑；σc為完整巖石的單軸抗壓強度；σmax=3σ1-σ3，為硐室圍巖的最大切向應(yīng)力，其中，σ1為最大主應(yīng)力，σ3為最小主應(yīng)力。

由圖1可以看出，硐室圍巖破壞深度與最大切向應(yīng)力之間具有良好的線性關(guān)系，通過式(1)反算，得到圍巖發(fā)生脆性破壞的臨界應(yīng)力為0.4σc(Df/a=1)。該方法針對應(yīng)力引起的完整或少量微裂隙的脆性圍巖破壞建立，將傳統(tǒng)的應(yīng)力強度比與圍巖破壞深度聯(lián)系起來，可以實現(xiàn)量化預(yù)測應(yīng)力集中下的圍巖脆性破壞深度。

1.2 巖爆定量化預(yù)測模型

基于上述思路，通過收集66個典型巖爆隧道的實例數(shù)據(jù)，對包括烏茲別克斯坦甘姆奇克隧道、二郎山隧道、錦屏二級水電站引水隧洞、拉林鐵路桑珠嶺隧道、花瓶山隧道等典型巖爆隧道的巖性、地應(yīng)力、抗壓強度、隧道跨度以及最大爆坑深度等參數(shù)進行了統(tǒng)計。典型巖爆隧道參數(shù)統(tǒng)計見表1。

表1 典型巖爆隧道統(tǒng)計Table 1 Statistics of typical rockburst tunnels

表中σ1為最大主應(yīng)力；σ3為最小主應(yīng)力；σc為巖石單軸抗壓強度

對收集的巖爆隧道的應(yīng)力強度比以及爆坑深度數(shù)據(jù)分析發(fā)現(xiàn)，應(yīng)力強度比與隧道巖爆的爆坑深度具有良好的線性關(guān)系。通過建立隧道巖爆爆坑深度與應(yīng)力強度比之間的關(guān)系，對收集的66組數(shù)據(jù)進行回歸分析，得到擬合公式：

(2)

擬合關(guān)系曲線如圖2所示?？梢钥闯觯佑绊懮疃菵f與應(yīng)力強度比σmax/σc之間具有良好的線性關(guān)系，同時，由式(2)推導(dǎo)出的巖爆起始應(yīng)力為0.32σc(Df/a=1)，與Martin et al. (1999)提出的臨界脆性破裂應(yīng)力范圍0.4σc接近。

圖2 巖爆量化預(yù)測模型Fig. 2 A universal quantitative rockburst prediction model

Hoek et al. (1980)根據(jù)南非金礦工程經(jīng)驗，提出采用圍巖切向應(yīng)力與巖石單軸抗壓強度比值σmax/σc來對巖爆進行評價并給出了巖爆等級劃分界限以及相應(yīng)的支護建議，其后得到廣泛的應(yīng)用。因此，綜合隧道爆坑深度和Hoek et al. (1980)提出的巖爆等級劃分依據(jù)，確定出本文預(yù)測模型中巖爆等級劃分依據(jù)如表2所示。

表2 巖爆等級分級標準Table 2 Classification of rockburst grades

本文建立預(yù)測模型所用參數(shù)包括巖石的單軸抗壓強度，最大切向應(yīng)力和爆坑影響深度等參數(shù)。最大切向應(yīng)力反映了對隧道圍巖脆性破壞發(fā)生起重要作用的原位應(yīng)力狀態(tài)，而采用爆坑影響深度，很好地考慮到巖爆對于整個隧道圍巖破壞程度的量化值。因此，基于該預(yù)測模型，通過室內(nèi)試驗獲得完整巖石的單軸抗壓強度，并結(jié)合現(xiàn)場原位地應(yīng)力測試，即可實現(xiàn)巖爆風(fēng)險的定量化預(yù)測。

1.3 結(jié)構(gòu)巖體強度損傷準則

Hoek et al. (2002)考慮到受結(jié)構(gòu)裂隙、風(fēng)化以及開挖過程中的工程擾動等因素的影響，巖體的強度相比于實驗室測得的巖塊強度將有大幅度的下降，并根據(jù)完整巖石的物理學(xué)試驗和裂隙巖體的模型研究結(jié)果，建立了一種考慮結(jié)構(gòu)損傷效應(yīng)的巖體強度估算理論。該方法的原理為：通過引入相關(guān)折減系數(shù)，對完整巖石的基本物理參數(shù)和力學(xué)參數(shù)進行折減，以獲取含結(jié)構(gòu)裂隙巖體的強度參數(shù)。該理論目前是較為完善的巖體強度估算方法之一(謝和平等， 2022)，在對巖體強度進行合理估計的同時，大大節(jié)約了巖體強度測試的成本和不方便性，在實際工程中得到廣泛的應(yīng)用。其最新修正表達式為(Hoek et al.，2019)：

(3)

式中：σ1、σ3分別為巖體破壞時的最大、最小主應(yīng)力； σc為完整巖石單軸抗壓強度；mb為巖石材料常數(shù)mi(范圍： 0.0000001～25)的折減系數(shù)，可根據(jù)式(4)確定；

(4)

式中：s和α為巖體常數(shù)，可分別由式(5)和式(6)確定：

(5)

(6)

式(4)～式(6)中：GSI為巖石的地質(zhì)強度指標；D為擾動常數(shù)，反映巖體受開挖、爆破等外界因素所影響的程度，介于0～1之間。

在開挖過程中，靠近隧道邊界的巖體處于臨空狀態(tài)，因此可認為σ3=0，此時巖體的強度可由式(3)簡化為：

σcm=σcisα

(7)

式中：σcm為巖體強度。

對于質(zhì)量控制良好的隧道掘進機開挖等方式，其對隧道圍巖的擾動可以控制在較小的水平，D可取值為0(Hoek et al.，2019)。此時，地質(zhì)強度指標GSI的大小將直接決定結(jié)構(gòu)巖體強度損傷的程度，即可根據(jù)室內(nèi)完整巖塊強度參數(shù)和GSI指標求取原位巖體強度。

2 工程實例應(yīng)用

2.1 工程概況及工程地質(zhì)條件

圖3 色季拉山隧道地質(zhì)縱斷面Fig. 3 Geological profile of Sejila Tunnel

表3 色季拉山隧道圍巖等級統(tǒng)計Table 3 Statistics of surrounding rock grade of Sejila Tunnel

2.2 工程區(qū)巖石力學(xué)性質(zhì)

將工程區(qū)取回的鉆孔巖芯加工成標準樣，分別進行單軸抗壓強度試驗和巴西劈裂實驗，得到工程區(qū)巖石力學(xué)參數(shù)統(tǒng)計結(jié)果見表4?？梢钥闯?，隧址區(qū)巖石抗壓強度較高，一般大于100MPa。巖石自身是否具有發(fā)生巖爆的本質(zhì)屬性，通常以巖爆傾向性指標進行評價。目前，反映巖石硬脆程度的脆性系數(shù)B，是目前應(yīng)用最簡便、廣泛的巖爆性傾向指標之一(王開洋等， 2014)，其計算公式為：

表4 隧道巖石基本力學(xué)參數(shù)Table 4 Statistics of basic mechanical parameters of rock of Sejila Tunnel

B=σc/σt

(8)

式中：σt為巖石的抗壓強度。

根據(jù)式(8)以及巖石強度參數(shù)測試結(jié)果，計算得到隧址區(qū)花崗巖、閃長巖以及片麻巖的脆性系數(shù)B分別為20.07、21.56以及16。徐林生等(1999)提出脆性系數(shù)B判斷巖爆傾向性的界限：B≥18時，強烈?guī)r爆； 18>B≥14，中等巖爆； 14>B≥10，輕微巖爆； 10>B無巖爆。因此，由式(8)可判斷出，除片麻巖為中等巖爆傾向外，閃長巖和花崗巖均具有強烈?guī)r爆傾向性。

隧址區(qū)巖體以弱風(fēng)化為主，巖體完整性總體較好，以完整、較完整巖體為主，揭露巖芯柱可長達0.5～1m，部分段落存在斷層破碎帶及節(jié)理密集帶?；谑覂?nèi)巖石物理力學(xué)試驗參數(shù)、野外鉆孔巖芯描述、工程物探及初步巖體分級等資料，采用定性與定量相結(jié)合的方法確定出隧道沿線巖體質(zhì)量指標GSI分布范圍，如圖4所示。利用前述Hoek-Brown巖體強度損傷準則及軟件RocLab對隧道軸線巖體強度參數(shù)進行分析估算，如表5所示。

圖4 隧道沿線巖體質(zhì)量GSI分布范圍Fig. 4 Distribution of GSI along the tunnel

表5 利用Hoek-Brown強度準則估算隧道軸線巖體強度Table 5 Strength of rock mass around tunnel section based on Hoek-Brown Criterion

2.3 隧址地應(yīng)力特征

隧址區(qū)位于喜馬拉雅東構(gòu)造結(jié)附近，區(qū)域內(nèi)新生代以來地殼急劇抬升、板塊強烈擠壓，造成了地應(yīng)力高且主應(yīng)力方向復(fù)雜的初始地應(yīng)力場特征。在前期勘察過程中，沿隧道沿線共布置地應(yīng)力測量孔10個，采用水壓致裂法對代表性測段進行測試，為了減小山坡卸荷和溝谷應(yīng)力集中等干擾因素，測段深度一般超過100m，最大孔深達1450m(張玉璽， 2021)。地應(yīng)力測試結(jié)果見表6。

表6 隧址區(qū)地應(yīng)力測量結(jié)果Table 6 In situ stress measurement results

為進一步分析隧址區(qū)地應(yīng)力狀態(tài)，根據(jù)水力壓裂測試結(jié)果分別計算出隧址區(qū)地應(yīng)力特征參數(shù)Kav、KHv和Khv，其表達式如下：

Kav=(SH+Sh)/2Sv

(9)

KHv=SH/Sv

(10)

Khv=Sh/Sv

(11)

式中：Kav為側(cè)壓力系數(shù)；KHv為最大水平應(yīng)力系數(shù)；Khv為最小水平應(yīng)力系數(shù)； 可以反映水平應(yīng)力相對于垂向自重應(yīng)力的關(guān)系。

隧址區(qū)地應(yīng)力特征參數(shù)隨深度變化的規(guī)律如圖5所示?？梢钥闯觯纠剿淼赖貞?yīng)力特征參數(shù)分布特征表現(xiàn)出相似的趨勢，可通過雙曲線進行擬合。在埋深較淺的近地表附近，各地應(yīng)力特征參數(shù)基本大于1，三向主應(yīng)力之間關(guān)系為SH>Sh>Sv，此時，反映出隧址區(qū)受水平構(gòu)造作用影響較大。同時，由于受地形切割、卸荷以及風(fēng)化剝蝕等作用的影響，淺部地應(yīng)力特征參數(shù)離散程度較大。隨著深度增大，隧道地應(yīng)力特征參數(shù)及其波動范圍逐漸減小。在埋深大于500m之后，Kav趨于1.14±0.2，KHv趨于1.36±0.26，Khv趨于0.97±0.19，三向主應(yīng)力之間關(guān)系逐漸轉(zhuǎn)變?yōu)镾H>Sv>Sh，隧址區(qū)的應(yīng)力場仍表現(xiàn)為以水平應(yīng)力為主的構(gòu)造控制作用，但自重應(yīng)力的影響顯著加強。

圖5 色季拉山隧道地應(yīng)力特征參數(shù)Fig. 5 Characteristic parameters of in-situ stress of the Sejila Tunnela. 側(cè)壓力系數(shù)Kav； b. 最大水平應(yīng)力系數(shù)KHv； c. 最小水平應(yīng)力系數(shù)Khv

2.4 初始地應(yīng)力場反演分析

地應(yīng)力是隧道工程圍巖發(fā)生變形破壞的關(guān)鍵性因素和重要工程荷載來源，直接影響隧道工程的施工安全(李天斌等， 2019)。由于場地、經(jīng)費等諸多因素的影響，難以做到對隧道全線地應(yīng)力進行水力壓裂測試(黃藝丹等， 2021)。因此，結(jié)合實地工程地質(zhì)條件及地應(yīng)力測試數(shù)據(jù)，進行地應(yīng)力數(shù)值模擬與反演，對于獲取分析隧址區(qū)及隧道沿線初始地應(yīng)力場具有重要意義。目前常用的地應(yīng)力場反演分析方法包括多元線性回歸分析法、邊界荷載調(diào)整法、應(yīng)力函數(shù)法等(王金安等， 2015； 李天斌等， 2016； 王慶武等， 2018)。其中：多元線性回歸法地應(yīng)力場反演概念清晰，在許多工程中已得到了成功的應(yīng)用和推廣(裴啟濤等， 2012)。因此，本文采用多元線性回歸法進行地應(yīng)力反演。

在Rhino中建立三維數(shù)值計算模型并剖分地層，模型構(gòu)建過程中，采用250m薄層來模擬斷層破碎帶及其影響帶。模型范圍如下：以隧道為中心區(qū)域，垂直于隧道走向為x軸，長42km； 以隧道走向為y軸， 16km； 豎直向上為z軸，模型基底以隧道出口點高程為基準向下延伸2000m。

結(jié)合研究區(qū)的工程條件，本文選取自重應(yīng)力狀態(tài)、南北向水平擠壓構(gòu)造運動、東西向水平擠壓構(gòu)造運動、水平面內(nèi)均勻剪切構(gòu)造運動等4個因素作為回歸巖體初始應(yīng)力場的基本要素，回歸初始地應(yīng)力場表達式為：

σ地=b1σ構(gòu)1+b2σ構(gòu)2+b3σ構(gòu)3+b4σ自4

(12)

式中：σ構(gòu)1、σ構(gòu)2分別為x、y方向施加的均勻擠壓應(yīng)力；σ構(gòu)3為xy方向施加的均勻應(yīng)力；σ自4為自重應(yīng)力，b1～b4為回歸系數(shù)。

將模型導(dǎo)入FLAC3D軟件中進行計算。模型所使用的隧址區(qū)巖體物理力學(xué)參數(shù)綜合現(xiàn)場地質(zhì)勘察資料、室內(nèi)試驗結(jié)果以及Hoek-Brown經(jīng)驗公式(Hoek et al.，2002)確定，如表7所示。采用彈性本構(gòu)模型對隧址區(qū)地應(yīng)力進行反演，以ZK6鉆孔實測地應(yīng)力值為回歸目標值，分別對式(12)確定的4種工況進行模擬計算，求解回歸方程。分別得到隧址區(qū)回歸系數(shù)b1、b2、b3、b4為1.89、0.14、1.05、0.45。

表7 隧址區(qū)巖體力學(xué)參數(shù)Table 7 Mechanical parameters of rock mass

經(jīng)過應(yīng)力場平衡計算，得到隧址區(qū)最大水平主應(yīng)力云圖如圖6所示，可以看出，隧址區(qū)核心部最大主應(yīng)力明顯大于隧道進出口段。

圖6 隧址區(qū)最大主應(yīng)力云圖Fig. 6 Maximum horizontal principal stress at tunnel site

考慮到靠近地表地應(yīng)力實測數(shù)據(jù)較離散，因此，選取ZK6埋深大于400m的地應(yīng)力實測數(shù)據(jù)與地應(yīng)力場反演結(jié)果進行對比，如圖7所示?？梢钥闯?，地應(yīng)力回歸數(shù)據(jù)與鉆孔實測地應(yīng)力數(shù)據(jù)比較接近，且兩者隨深度的變化具有較好的一致性。各測點實測值與回歸模擬值相對誤差計算結(jié)果見表8?？梢钥闯?，地應(yīng)力回歸模擬值與實測值的誤差基本控制在10%以內(nèi)，表現(xiàn)出較高的準確率。

圖7 地應(yīng)力實測與回歸值對比Fig. 7 Comparison between measured and simulated results of Sejila Tunnel

表8 地應(yīng)力回歸值誤差Table 8 Comparison between measured and simulated results of Sejila Tunnel

數(shù)值模擬得到的色季拉山隧道沿線地應(yīng)力分布如圖8所示?？梢钥闯鏊淼姥鼐€最大水平主應(yīng)力為6.53～51.24MPa，最小水平主應(yīng)力為4.8～36.12MPa，垂向主應(yīng)力為2.26～43.62MPa。全線最大水平主應(yīng)力量值水平較高，在斷層帶附近有明顯下降。在淺埋深段水平主應(yīng)力大于鉛垂應(yīng)力，三向主應(yīng)力關(guān)系為SH>Sh>Sv，表現(xiàn)出明顯的構(gòu)造應(yīng)力為主，當埋深超過500m，垂向應(yīng)力大于最小水平主應(yīng)力，三向主應(yīng)力關(guān)系轉(zhuǎn)變?yōu)镾H>Sv>Sh。

圖8 色季拉山隧道沿線地應(yīng)力Fig. 8 Principal stress value of Sejila Tunnel

2.5 色季拉山隧道巖爆量化預(yù)測

根據(jù)色季拉山隧道沿線巖體抗壓強度估算結(jié)果和地應(yīng)力反演得到的隧道軸線地應(yīng)力參數(shù)，采用式(2)對色季拉山隧道巖爆風(fēng)險進行量化預(yù)測。為此，首先需要進行隧道最大切向應(yīng)力的計算。色季拉山隧道開挖斷面可近似于圓形，故可將隧道橫截面假定為平面應(yīng)變模型(張重遠等， 2022)，此時，隧道截面方向受到水平向正應(yīng)力σn和上覆巖石自重作用引起的垂向應(yīng)力σv，根據(jù)線彈性公式，水平向正應(yīng)力計算公式為：

(13)

式中：θ為最大水平主應(yīng)力SH與隧道軸向法線的夾角，根據(jù)水力壓裂測試結(jié)果，色季拉山隧道θ為64°。最終，可以得到色季拉山隧道最大切向應(yīng)力計算公式：

(14)

對于斷層破碎帶等不具備彈性應(yīng)變儲存條件的巖體強烈破碎區(qū)域，在較高的應(yīng)力下可能發(fā)生大變形而非巖爆，因此在預(yù)測及統(tǒng)計隧道沿線巖爆時對該段落進行去除。最終計算得到色季拉山隧道沿線巖爆預(yù)測結(jié)果如圖9所示。

圖9 色季拉山隧道最大爆坑深度預(yù)測結(jié)果Fig. 9 Prediction results of maximum explosion pit depth of Sejila Tunnel

可以看出，除斷層影響地帶和埋深較淺的溝谷、進出口地帶，隧道沿線均可能有不同程度的巖爆出現(xiàn)。其中：在埋深較大的地段，花崗巖、閃長巖地層將出現(xiàn)強烈等級程度的巖爆，片麻巖地層中將有中等程度巖爆發(fā)生。

同時對沿線巖爆段落進行統(tǒng)計，得到色季拉山隧道巖爆評價結(jié)果，如表9所示?？梢钥闯?，發(fā)生巖爆的段落普遍具有較高的地應(yīng)力，最大水平主應(yīng)力SH普遍高于20MPa，部分段落地應(yīng)力可達40MPa以上，有強烈等級的巖爆發(fā)生。對不同等級巖爆段落長度進行統(tǒng)計，結(jié)果顯示，色季拉山隧道預(yù)測巖爆段落達23254m，占隧道全長的61.25%，其中：中等巖爆段落達12188m，占隧道全長32.1%，強烈等級巖爆段落達3502m，占隧道全長的9.23%。預(yù)測隧道爆坑最大深度比Df/a可達1.67(1.46+0.21)，按色季拉山隧道當量半徑5.1m計算，預(yù)測最大爆坑深度可達3.42m。

表9 色季拉山隧道巖爆評價結(jié)果Table 9 Statistical chart of rockburst prediction results

2.6 色季拉山巖爆風(fēng)險可控性分析

為評估色季拉山隧道巖爆風(fēng)險的可控性，選取已成功貫通的CZ鐵路拉林段巴玉隧道進行相似工程案例風(fēng)險的定量對比分析。工程地質(zhì)類比法屬于定性評價的一種，它的實質(zhì)是將已有的隧道風(fēng)險進行總結(jié)歸納后，應(yīng)用到地質(zhì)條件相似的隧道風(fēng)險的評價研究中，是一種應(yīng)用廣泛的經(jīng)驗方法(張寧等， 2022)。

巴玉隧道位于桑日至加查下游段，與色季拉山隧道同屬于雅魯藏布縫合帶影響區(qū)，隧道全長13.047km，最大埋深約2073m。巴玉隧道巖爆問題極為突出，據(jù)13047m正洞和8131m平導(dǎo)的統(tǒng)計結(jié)果顯示，隧道巖爆里程段高達94%(嚴健等， 2019)。

表10分別列出了兩座隧道的基本概況和基本地質(zhì)條件，同時也分別給出了巴玉隧道實測最大爆坑深度和色季拉山隧道預(yù)測最大爆坑深度?？梢钥闯觯瑑勺淼谰哂袠O其相似的工程地質(zhì)條件，兩座隧道地層巖性均以花崗巖為主，且?guī)r石強度較為接近。此外，色季拉山隧道最大主應(yīng)力為51.24MPa，略大于巴玉隧道最大主應(yīng)力45.77MPa。根據(jù)預(yù)測結(jié)果，色季拉山隧道巖爆的最大爆坑深度可達3.42m，略小于巴玉隧道實測最大爆坑深度3.5m，同時，色季拉山隧道預(yù)測巖爆里程段占比為61.25%，遠小于巴玉隧道實際巖爆段占比94%。因此，可以認為，在合理的工程措施下，色季拉山隧道巖爆風(fēng)險總體可控。

表10 典型隧道巖爆情況對比Table 10 Comparison of rockburst conditions of typical tunnels

3 結(jié) 論

本文在對66個巖爆隧道的特征參數(shù)統(tǒng)計分析基礎(chǔ)上，建立了一種通過應(yīng)力強度比來預(yù)測巖爆爆坑深度的巖爆定量化預(yù)測模型，并通過室內(nèi)力學(xué)試驗、巖體強度損傷準則和地應(yīng)力場多元線性回歸反演，分別得到模型計算所需要的隧道沿線巖體強度及地應(yīng)力參數(shù)，最終對新建CZ鐵路色季拉山隧道的巖爆風(fēng)險進行了定量化預(yù)測和相似工程案例對比分析，得到如下結(jié)論：

(1)本文根據(jù)巖爆隧道爆坑深度與應(yīng)力強度比線性統(tǒng)計關(guān)系所建立的巖爆風(fēng)險預(yù)測模型，可以預(yù)測巖爆發(fā)生時的最大爆坑影響深度，進而實現(xiàn)對巖爆風(fēng)險的定量化預(yù)測及評估。

(2)色季拉山隧道沿線地應(yīng)力以水平構(gòu)造應(yīng)力為主導(dǎo)，反演出最大水平主應(yīng)力可達51.24MPa，且地應(yīng)力普遍較高，同時，硬脆性巖層所占比例較大且均具有較高的巖爆傾向性，具有較高的高地應(yīng)力巖爆風(fēng)險。

(3)預(yù)測色季拉山隧道可能發(fā)生巖爆段落達23254m，占隧道全長的61.25%，其中：中等巖爆風(fēng)險段落12188m，占隧道全長32.1%，強烈?guī)r爆風(fēng)險段落3502m，占隧道全長的9.23%。巖爆將主要發(fā)生在弱風(fēng)化花崗巖、閃長巖以及埋深較大的片麻巖段落中，且最大爆坑深度普遍大于1m。

(4)相似工程案例對比分析表明，色季拉山隧道預(yù)測最大爆坑深度為3.42m，小于已貫通的巴玉隧道實測最大爆坑深度(3.5m)，在合理的施工和防控措施下可認為風(fēng)險總體可控。