孫 闖，敖云鶴，張家鳴

(遼寧工程技術(shù)大學(xué) 土木工程學(xué)院，遼寧 阜新 123000)

0 引言

我國(guó)公路隧道建設(shè)迅速發(fā)展，小凈距公路隧道因經(jīng)濟(jì)、社會(huì)和環(huán)境效益顯著，逐漸成為山區(qū)狹窄地形條件下修建隧道的主要結(jié)構(gòu)形式之一。在弱節(jié)理復(fù)雜地質(zhì)條件下對(duì)小凈距隧道進(jìn)行設(shè)計(jì)、施工及安全穩(wěn)定性研究，具有十分重要的現(xiàn)實(shí)意義[1-2]。

在弱節(jié)理小凈距隧道研究中，關(guān)鍵科學(xué)問(wèn)題是合理確定弱節(jié)理圍巖的力學(xué)參數(shù)，提出準(zhǔn)確的計(jì)算分析方法并采取可行的支護(hù)和加固措施。針對(duì)以上問(wèn)題，國(guó)內(nèi)外學(xué)者已經(jīng)開(kāi)展了系列研究。孫闖等[3]采用超前地質(zhì)預(yù)報(bào)和圍巖評(píng)級(jí)系統(tǒng)相結(jié)合的方法，合理確定了隧道節(jié)理巖體的GSI并獲得符合工程實(shí)際的圍巖力學(xué)參數(shù)；在數(shù)值計(jì)算方面，王康[4]修正了已有的摩爾-庫(kù)倫屈服接近度，并基于彈塑性有限差分?jǐn)?shù)值模擬，對(duì)不同圍巖等級(jí)的隧道給出了合理的施工方案建議；O.Jenck[5]等采用數(shù)值分析方法分析了淺埋小凈距隧道開(kāi)挖過(guò)程中對(duì)既有線襯砌結(jié)構(gòu)的變形和內(nèi)力的影響。在中夾巖柱研究方面，周佳媚等[6]采用三維數(shù)值模擬的手段，得出小凈距隧道中夾巖柱應(yīng)力狀態(tài)與其凈距的關(guān)系；應(yīng)國(guó)剛[7]等通過(guò)預(yù)留支護(hù)界面上的正應(yīng)力和剪應(yīng)力分布表達(dá)式，計(jì)算中夾巖柱的平均應(yīng)力和附加彎矩從而得到其應(yīng)力分布規(guī)律。目前，在小凈距隧道相關(guān)研究方面已經(jīng)獲得了一定的成果，但在弱節(jié)理小凈距隧道合理凈距及圍巖穩(wěn)定性方面仍待深入研究。

本研究在前人的基礎(chǔ)上，以大荒溝小凈距隧道工程為研究背景，基于量化GSI圍巖評(píng)級(jí)系統(tǒng)，確定圍巖力學(xué)參數(shù)，基于Hoek-Brown強(qiáng)度準(zhǔn)則的應(yīng)變軟化模型，采用 FLAC3D數(shù)值模擬得到不同凈距隧道中夾巖柱塑性區(qū)分布，確定最小合理凈距，理論計(jì)算合理凈距條件下的圍巖壓力，并進(jìn)行支護(hù)結(jié)構(gòu)設(shè)計(jì)，對(duì)比分析不同開(kāi)挖方式下小凈距隧道圍巖變形量，確定合理的開(kāi)挖及支護(hù)方案。

1 弱節(jié)理隧道圍巖參數(shù)及力學(xué)模型

1.1 工程概況

大荒溝隧道位于丹東市寬甸縣青山溝鄉(xiāng)彎溝村與本溪市桓仁縣向陽(yáng)鄉(xiāng)和平村之間，走向南西212°～238°左右。隧道最大埋深330 m，隧道長(zhǎng)度1 130 m。通過(guò)現(xiàn)場(chǎng)地質(zhì)勘察及隧道工作面觀察可知，大荒溝隧道小凈距段巖石主要為砂質(zhì)泥巖，節(jié)理裂隙發(fā)育，屬于弱節(jié)理巖層，巖體整體強(qiáng)度較差，巖層節(jié)理分布特征如圖1所示。

圖1 構(gòu)造破碎帶及節(jié)理Fig.1 Structural fracture zone and joint

1.2 量化GSI圍巖評(píng)級(jí)系統(tǒng)

合理確定弱節(jié)理小凈距隧道圍巖的力學(xué)參數(shù)，在工程設(shè)計(jì)中起著至關(guān)重要的作用。E.Hoek和Brown在大量的工程實(shí)踐中發(fā)現(xiàn)，對(duì)于軟弱、破碎的巖體結(jié)構(gòu)，GSI圍巖評(píng)級(jí)系統(tǒng)能夠較為理想的反應(yīng)節(jié)理巖體特征。隨著GSI圍巖評(píng)級(jí)系統(tǒng)的廣泛使用，M.Cai等[8]提出了量化GSI 系統(tǒng)的方法，其中包括結(jié)構(gòu)面條件因子Jc和塊體尺寸Vb，量化的GSI方法可以讓工程師在現(xiàn)場(chǎng)更客觀的對(duì)節(jié)理巖體參數(shù)的進(jìn)行取值。具體方法如圖2所示。

圖2 量化GSI圍巖評(píng)級(jí)系統(tǒng)Fig.2 Quantitative GSI surrounding rock rating system

結(jié)合現(xiàn)場(chǎng)勘查、工作面觀測(cè)結(jié)果，估算該區(qū)段弱節(jié)理圍巖峰值GSI在52～55之間，峰后殘余GSI在31～34之間，如圖3中峰值區(qū)域及殘余強(qiáng)度區(qū)域所示。現(xiàn)場(chǎng)對(duì)隧道內(nèi)砂質(zhì)泥巖進(jìn)行取樣，室內(nèi)采用朝陽(yáng)GAW-2000微機(jī)控制電液伺服剛性壓力試驗(yàn)機(jī)開(kāi)展試驗(yàn)，通過(guò)現(xiàn)場(chǎng)取樣獲得隧道工作面砂質(zhì)泥巖巖樣，實(shí)驗(yàn)獲得砂質(zhì)泥巖的單軸抗壓強(qiáng)度σci=15 MPa，彈性模量為E=1.1 GPa，泊松比υ=0.3。

1.3 應(yīng)變軟化力學(xué)模型

基于Hoek-Brown(H-B)屈服準(zhǔn)則[9]，即：

(1)

式中，σ1，σ3分別為巖體破壞時(shí)最大、最小主應(yīng)力；σci為巖塊單軸抗壓強(qiáng)度；mb，s，a均為常數(shù)，其中mb為mi的折算值，s反映巖體破壞程度，其取值范圍在0～1之間。該準(zhǔn)則將擾動(dòng)系數(shù)D引入到巖體強(qiáng)度計(jì)算當(dāng)中，各參數(shù)可由下式表示：

(2)

若假設(shè)H-B常數(shù)中的mb，s隨η值線性衰減，則可得到H-B應(yīng)變軟化模型，其關(guān)系式為：

(3)

式中，ωp為峰值參數(shù);ωr為殘余參數(shù)；ω可以代替H-B模型中的mb；s常數(shù)。

將軟化參數(shù)η定義為塑性剪切應(yīng)變，通過(guò)最大主塑性應(yīng)變和最小主塑性應(yīng)變的差值獲得，即

(4)

根據(jù)現(xiàn)場(chǎng)掌子面節(jié)理巖體特征及觀測(cè)結(jié)果，取GSIp=52，GSIr=32，確定基于量化GSI系統(tǒng)的圍巖力學(xué)參數(shù)如表1所示。隧道弱節(jié)理圍巖的彈性模量根據(jù)Hoek提出的方法確定[10]，即

表1 弱節(jié)理巖體力學(xué)參數(shù)Tab.1 Mechanical parameters of weak jointed rock mass

2 弱節(jié)理小凈距隧道合理凈距的確定

本研究以FLAC3D有限差分軟件為工具，基于Hoek-Brown強(qiáng)度準(zhǔn)則的應(yīng)變軟化模型，建立三維小凈距隧道數(shù)值模型，數(shù)值分析小凈距隧道中夾巖柱的塑性區(qū)分布特征及塑性區(qū)的突變規(guī)律，綜合確定隧道合理凈距，計(jì)算工況分別取小凈距為0.9B～1.7B等9種，其中B為隧道的跨度。

2.1 弱節(jié)理小凈距隧道數(shù)值計(jì)算模型

大荒溝小凈距隧道建筑限界寬10.25 m，凈高5 m，雙車道隧道。內(nèi)輪廓采用曲墻三心圓形式，寬10.84 m，高8.9 m。計(jì)算模型力學(xué)邊界條件采用位移邊界條件，上邊界為自由邊界，左右兩邊邊界采用水平位移+應(yīng)力約束，更加準(zhǔn)確地還原隧道的原巖應(yīng)力狀態(tài)，固定縱向方向的位移，底部邊界為固定邊界，根據(jù)公式P=λγh可知，P1=6.8 MPa，P2=3.4 MPa，側(cè)壓力系數(shù)λ為0.5。

圖3 隧道數(shù)值模型及邊界條件Fig.3 Tunnel numerical model and boundary conditions

2.2 數(shù)值計(jì)算結(jié)果分析

本研究以FLAC3D數(shù)值模擬軟件計(jì)算0.9B～1.7B凈距條件下，雙洞開(kāi)挖后圍巖塑性區(qū)分布情況。對(duì)小凈距隧道合理凈距的判定條件，主要包括兩個(gè)方面，一是小凈距隧道中夾巖柱不出現(xiàn)塑性區(qū)貫通[11-12]，二是塑性區(qū)頂部不能出現(xiàn)突變情況，若出現(xiàn)上部塑性區(qū)突變則說(shuō)明中夾巖柱易產(chǎn)生塌方隱患，通過(guò)計(jì)算獲得不同凈距隧道中夾巖柱的塑性區(qū)分布如圖4所示。由數(shù)值計(jì)算結(jié)果可知，當(dāng)隧道小凈距D取 0.9B～1.1B時(shí)，塑性區(qū)在中間巖柱位置全部處于貫通狀態(tài)，說(shuō)明在這種條件下不滿足小凈距隧道的施工要求；當(dāng)隧道小凈距D取 1.2B～1.4B時(shí)，隧道中夾巖柱的塑性區(qū)并沒(méi)有貫通，圍巖具有一定的自穩(wěn)能力，但是隧道的肩部出現(xiàn)明顯突變現(xiàn)象，說(shuō)明在這種條件下開(kāi)挖時(shí)，隧道中夾巖柱仍然存在塌方隱患；當(dāng)D=1.5B時(shí)，塑性區(qū)頂部突變現(xiàn)象明顯減小，且兩硐室塑性區(qū)處于分離狀態(tài)，隧道圍巖趨于穩(wěn)定。根據(jù)綜合分析大荒溝隧道地質(zhì)條件和數(shù)值計(jì)算結(jié)果，將其最小合理凈距取值確定為D=1.5B較為適宜。

圖4 不同凈距隧道中夾巖柱的塑性區(qū)分布Fig.4 Distribution of plastic zone of middle rock pillar in tunnels with different clear distances

3 弱節(jié)理小凈距隧道支護(hù)結(jié)構(gòu)設(shè)計(jì)

小凈距隧道雙洞內(nèi)側(cè)應(yīng)力狀態(tài)較單洞情況復(fù)雜[13-15]，其影響因素不僅包括隧道尺寸和埋深等工程條件，還包括圍巖參數(shù)等地質(zhì)條件，本研究以松散介質(zhì)平衡理論(普式理論)為基礎(chǔ)，計(jì)算深埋隧道開(kāi)挖后圍巖的垂直壓力及兩側(cè)水平壓力[16]。

3.1 小凈距隧道圍巖壓力

(1)垂直壓力

深埋隧道垂直壓力由基本松散壓力q1和附加松散壓力q2組成，可按下式計(jì)算：

外側(cè)：q外=q1+q2=γ(hq1+hq2)=68.37 kPa,

(6)

內(nèi)側(cè)：q內(nèi)=q1+q′2=γ(hq1+h′q2)=102.92 kPa,

(7)

式中，q1為小凈距隧道的基本垂直壓力荷載，單側(cè)洞室形成的穩(wěn)定平衡拱下部的土壓力, 假定其為均布荷載；q2為小凈距隧道的附加垂直壓力荷載。左右洞室共同形成的極限平衡拱下部松散土壓力減去基本松散土壓力及中夾巖柱體上部土壓力荷載后的荷載, 假定為梯形分布荷載,如圖5所示。hq1為基本荷載高度；hq2為附加荷載高度；γ為拱頂附近巖體的平均容重。

圖5 小凈距隧道荷載分布計(jì)算示意圖Fig.5 Schematic diagram of load distribution calculation of small clear distance tunnel

(2)兩側(cè)水平壓力

(8)

(9)

3.2 支護(hù)結(jié)構(gòu)設(shè)計(jì)

(1)錨桿支護(hù)設(shè)計(jì)

系統(tǒng)錨桿長(zhǎng)度和間距應(yīng)根據(jù)圍巖的破碎程度、隧道寬度，采用工程類比法確定，經(jīng)驗(yàn)公式：

① 錨桿長(zhǎng)度：L=1/5Bt～1/3Bt=2.17～3.61，取L=3 m；

② 錨桿間距：P=0.4L～0.7L=1.2～2.1，取P=1.2 m。

式中，L為錨桿長(zhǎng)度；Bt隧道開(kāi)挖寬度，Bt=10.84 m；P為錨桿的設(shè)置間距。使用HRB400直徑為25 mm的錨桿，矩形排列縱向間距和排距為1.2 m 和HRB235的鋼墊。根據(jù)文獻(xiàn)[17]可知，錨桿所能提供的最大承載力為0.139 MPa。錨桿參數(shù)如表2所示。

表2 錨桿參數(shù)Tab.2 Parameters of rock bolts

(2)噴射混凝土支護(hù)設(shè)計(jì)

選取C30普通硅酸鹽水泥，支護(hù)噴射混凝土厚度為20 cm。根據(jù)文獻(xiàn)可知，混凝土所能提供的最大承載為1.071 MPa?；炷羺?shù)如表3所示。

表3 混凝土參數(shù)Tab.3 Parameters of concrete

(3)中夾巖柱的支護(hù)設(shè)計(jì)

本隧道凈距為D=16.26 m，根據(jù)現(xiàn)場(chǎng)的地質(zhì)條件選取中夾巖柱加固措施，即使用直徑為25 mm的HRB400預(yù)應(yīng)力錨桿在中夾巖柱區(qū)域連接兩個(gè)洞室，以達(dá)到對(duì)中巖墻的加固，選用C30混凝土對(duì)上盤巖、下盤巖進(jìn)行注漿加固。

4 圍巖開(kāi)挖方案優(yōu)化及穩(wěn)定性分析

采用FLAC3D軟件分別對(duì)支護(hù)條件下大荒溝小凈距隧道采用二臺(tái)階法、單側(cè)壁導(dǎo)坑法、雙側(cè)壁導(dǎo)坑法、CD法、CRD法、環(huán)形開(kāi)挖預(yù)留核心土法進(jìn)行數(shù)值模擬，分析不同工況下圍巖拱頂位移及水平位移的分布特征，選擇合理的施工方法，開(kāi)挖方法如圖6所示。其中二臺(tái)階法采用短臺(tái)階開(kāi)挖，臺(tái)階步距為0.9B，單側(cè)壁導(dǎo)坑法、雙側(cè)壁導(dǎo)坑法的開(kāi)挖步距為0.6B，CD法、CRD法的開(kāi)挖步距為0.5B，環(huán)形開(kāi)挖預(yù)留核心土法的上部開(kāi)挖步距為0.1B，下部臺(tái)階開(kāi)挖步距為0.4B。

圖6 小凈距隧道不同開(kāi)挖方法圖Fig.6 Different excavation methods for small clear distance tunnel

從圖7中可以看出不同的開(kāi)挖方法對(duì)小凈距隧道拱頂沉降量有一定的影響，環(huán)形開(kāi)挖預(yù)留核心土法對(duì)拱頂沉降量控制最為明顯。而隧道左洞、右洞之間的沉降量無(wú)較大差別，保持為1 mm左右。

從表4中可以看出，拱頂最終沉降量環(huán)形開(kāi)挖預(yù)留核心土法最小，然后分別是雙側(cè)壁導(dǎo)坑法、單側(cè)壁導(dǎo)坑法、CD法、CRD法、二臺(tái)階法。

圖7 不同開(kāi)挖方法時(shí)后方拱頂沉降曲線Fig.7 Settlement curves of rear vault using different excavation methods

圖8為不同工況下中間巖柱水平位移變化曲線(向左為負(fù)值，向右為正值)，由文獻(xiàn)[18]中提出的小凈距隧道中間巖柱的分布及加固法可知，中間巖柱的穩(wěn)定性對(duì)小凈距隧道施工起著至關(guān)重要的作用，對(duì)中間巖柱進(jìn)行加固可減小隧道圍巖變形，控制圍巖塑性區(qū)發(fā)展，提高隧道圍巖的穩(wěn)定性。無(wú)論是中間巖柱的左側(cè)還是右側(cè)，環(huán)形開(kāi)挖預(yù)留核心土法施工時(shí)中間巖柱的水平位移最小，二臺(tái)階法水平位移最大，由此可知在應(yīng)用環(huán)形開(kāi)挖預(yù)留核心土法施工時(shí)對(duì)隧道圍巖造成的擾動(dòng)最小。對(duì)比巖柱向左和向右移動(dòng)的水平位移可知，隧道開(kāi)挖時(shí)，中間巖柱向左洞移動(dòng)趨勢(shì)明顯大于右洞，圍巖有先向左洞移動(dòng)的趨勢(shì)。

表4 隧道拱頂最終沉降量(單位：mm)Tab.4 Final settlement of tunnel vault (unit: mm)

圖8 不同工況中間巖柱的水平位移曲線Fig.8 Horizontal displacement curves of middle rock column under different working conditions

5 小凈距隧道施工方案及監(jiān)測(cè)分析

5.1 施工方案的選擇

通過(guò)對(duì)不同工況的數(shù)值模擬分析可知，環(huán)形開(kāi)挖預(yù)留核心土法、雙側(cè)壁導(dǎo)坑法、CRD法引起的隧道圍巖擾動(dòng)較小。從整體上來(lái)看，環(huán)形開(kāi)挖預(yù)留核心土法在控制圍巖變形方面優(yōu)于后兩者。而環(huán)形開(kāi)挖留有核心土支撐著開(kāi)挖面，能夠及時(shí)進(jìn)行初期支護(hù)，所以開(kāi)挖面穩(wěn)定性好，其次核心土和下部開(kāi)挖是在初期支護(hù)保護(hù)下進(jìn)行的，施工安全性好。故選擇環(huán)形開(kāi)挖預(yù)留核心土法進(jìn)行隧道施工。

5.2 現(xiàn)場(chǎng)監(jiān)測(cè)分析

為驗(yàn)證數(shù)值模擬的可靠性，對(duì)大荒溝小凈距隧道采用環(huán)形開(kāi)挖預(yù)留核心土法施工，采用國(guó)產(chǎn)XB338-B型智能數(shù)顯滑動(dòng)式沉降儀對(duì)拱頂沉降進(jìn)行量測(cè)，采用激光斷面儀對(duì)水平收斂進(jìn)行量測(cè)，觀測(cè)點(diǎn)布置如圖9所示。

圖9 隧道現(xiàn)場(chǎng)監(jiān)測(cè)點(diǎn)布置圖Fig.9 Layout of tunnel field monitoring points

圖10 現(xiàn)場(chǎng)監(jiān)測(cè)變形量變化曲線Fig.10 Curves of field monitoring deformation

圖10為DK424+205～+255段拱頂沉降及水平收斂位移隨測(cè)點(diǎn)距掌子面距離的變化曲線。從圖10中可以看出，隧道拱頂沉降量在初始階段迅速增大，隨著隧道的掘進(jìn)，拱頂沉降量逐漸趨于穩(wěn)定，在距掌子面35 m處趨于穩(wěn)定，達(dá)到8.02 mm，小于數(shù)值分析結(jié)果。從圖10中還可以看出，水平收斂量隨著隧道掘進(jìn)逐漸增大，在距掌子面35 m左右趨于穩(wěn)定，達(dá)到3.71 mm?，F(xiàn)場(chǎng)監(jiān)測(cè)數(shù)據(jù)值與數(shù)值模擬分析的最終結(jié)果基本吻合，從圖中還可以看出，隧道開(kāi)挖造成的擾動(dòng)范圍大約為1.5B洞跨。

通過(guò)對(duì)現(xiàn)場(chǎng)監(jiān)測(cè)結(jié)果及數(shù)值模擬的分析可知，環(huán)形開(kāi)挖預(yù)留核心土法在施工安全方面具有可行性，環(huán)形開(kāi)挖預(yù)留核心土法能夠更加有效的控制隧道圍巖的穩(wěn)定性，結(jié)合實(shí)際工程，選擇1.5B凈距及環(huán)形開(kāi)挖預(yù)留核心土法施工更加合理。

6 結(jié)論

(1)采用量化GSI圍巖評(píng)級(jí)系統(tǒng)，能夠確定合理的節(jié)理巖體力學(xué)參數(shù)，基于Hoek-Brown強(qiáng)度準(zhǔn)則的應(yīng)變軟化模型，采用 FLAC3D數(shù)值模擬軟件計(jì)算確定最小合理凈距，合理凈距的判定條件為小凈距隧道中夾巖柱不出現(xiàn)塑性區(qū)貫通，塑性區(qū)頂部不能出現(xiàn)突變情況。

(2)按照規(guī)范理論計(jì)算圍巖壓力后，進(jìn)行支護(hù)結(jié)構(gòu)設(shè)計(jì)，采用FLAC3D數(shù)值模擬軟件對(duì)小凈距隧道不同開(kāi)挖方式進(jìn)行計(jì)算，確定了環(huán)形開(kāi)挖法為理想的施工方法；得出錨噴支護(hù)與中夾巖柱長(zhǎng)錨桿共同作用，能夠有效控制小凈距隧道中夾巖柱松動(dòng)圈的擴(kuò)展。

(3)通過(guò)現(xiàn)場(chǎng)監(jiān)測(cè)驗(yàn)證了施工方案的合理性，在計(jì)算小凈距隧道穩(wěn)定性方面，未考慮巖體的時(shí)效變形特性，可能導(dǎo)致計(jì)算結(jié)果存在偏差，這方面內(nèi)容需進(jìn)一步研究。