石廣斌 肖 清 魏娟盆 張 雯 洪 勇

(1.西安建筑科技大學(xué)資源工程學(xué)院,陜西 西安 710055;2.陜西省巖土與地下空間重點(diǎn)實(shí)驗(yàn)室,陜西 西安 710055)

隨著我國(guó)經(jīng)濟(jì)的快速發(fā)展,深部資源得以開(kāi)發(fā),隨之而來(lái)的深部巖體所具有的特殊工程地質(zhì)問(wèn)題也更加地突出[1-4]。由于地下硐室群建在深部巖體中,開(kāi)挖可能引起圍巖釋放壓力,導(dǎo)致大型硐室圍巖變形。當(dāng)圍巖的變形超過(guò)周?chē)鷰r體的承受能力時(shí),硐室群巖體發(fā)生滑動(dòng)和破壞[5-8]。目前對(duì)大型地下硐室群圍巖破壞的研究方法主要有三維數(shù)值模擬法,如肖益蓋等[9]采用三維數(shù)值仿真模擬計(jì)算方法對(duì)某大型地下巖體開(kāi)挖工程進(jìn)行整體性穩(wěn)定評(píng)價(jià);代碧波等[10]依據(jù)礦山的開(kāi)采計(jì)劃,建立了數(shù)值模型,分析了采動(dòng)影響下巷道圍巖的應(yīng)力、塑性區(qū)和位移變化規(guī)律。根據(jù)前期地質(zhì)勘探資料顯示,某大型地下硐室群存在不良地質(zhì)體,即PH擠壓帶為發(fā)育在T3Z2[3]巖組內(nèi)的2條規(guī)模較大的碳質(zhì)千枚巖帶,其物理力學(xué)指標(biāo)在Ⅳ2類(lèi)及以下,直接影響著大跨度地下硐室群圍巖穩(wěn)定。因此,本研究將通過(guò)三維彈塑性有限元方法模擬硐室群巖體開(kāi)挖的力學(xué)過(guò)程,基于塑性區(qū)分布建立力學(xué)模型。根據(jù)毛洞開(kāi)挖圍巖變形提出永久性和連續(xù)性措施,為硐室群圍巖支護(hù)措施的制定提供科學(xué)的依據(jù)。

1 工程地質(zhì)概況

某大型地下硐室群位于左岸下壩線稍靠下游山體內(nèi),主要由主廠房和主變室等平行布置,硐室群垂直最大埋深為250 m,主廠房開(kāi)挖尺寸(長(zhǎng)×寬×高)為183.0 m×25.8 m×65.25 m,主變室開(kāi)挖尺寸(長(zhǎng)×寬×高)為117.2 m×16.5 m×30.8 m,硐室群間距為39.2 m。

硐室群出露地層巖性為三疊系上統(tǒng)雜谷腦組上段(T3Z2)薄～厚層狀變質(zhì)細(xì)砂巖夾碳質(zhì)千枚巖。其外圍構(gòu)造簡(jiǎn)單,無(wú)大斷層通過(guò),根據(jù)前期地質(zhì)勘探揭露該硐室群存在不良地質(zhì)體,主要為層間擠壓錯(cuò)動(dòng)帶(PH條帶)(圖1)。PH條帶的物理力學(xué)指標(biāo)等同于Ⅳ2類(lèi)圍巖,在此地質(zhì)條件下對(duì)大跨度地下硐室群開(kāi)挖將面臨圍巖變形較大,從而造成圍巖變形失穩(wěn)或破壞?？紤]邊界條件、荷載情況等,力學(xué)參數(shù)值見(jiàn)表1。

圖1 工程地質(zhì)圖Fig.1 Engineering geological map

表1 巖體質(zhì)量分級(jí)及力學(xué)參數(shù)Table 1 Rock mass quality classification and mechanical parameters

2 圍巖力學(xué)本構(gòu)模型

硐室群巖體開(kāi)挖過(guò)程中,圍巖初始應(yīng)力因破壞而重新分布,當(dāng)釋放的應(yīng)力大于圍巖強(qiáng)度時(shí),圍巖產(chǎn)生的屈服破壞轉(zhuǎn)變成塑性變形,為防止硐室出現(xiàn)坍塌,硐室群圍巖的破壞問(wèn)題多采用Mohr-Coulomb屈服準(zhǔn)則。因此為便于分析,現(xiàn)假設(shè)硐室群圍巖為均質(zhì)各向同性的連續(xù)介質(zhì),本研究在眾多研究者的基礎(chǔ)上[11、12],基于 Mohr-Coulomb屈服準(zhǔn)則,考慮圍巖開(kāi)挖—支護(hù)過(guò)程中錨固效應(yīng),對(duì)圍巖進(jìn)行彈塑性分析。

3 不良地質(zhì)體條帶寬度對(duì)圍巖穩(wěn)定影響

3.1 數(shù)值分析模型

根據(jù)硐室群實(shí)際尺寸建立地下硐室群三維整體模型。模型坐標(biāo)系采用右手法則,X軸平行于硐室群縱軸線,Y軸垂直于硐室群縱軸線,Z軸指向上為正。模型兩邊為Ⅲ類(lèi)圍巖,中間為Ⅳ類(lèi)圍巖(代表PH條帶),通過(guò)改變Ⅳ類(lèi)圍巖的寬度,即PH條帶寬度,研究不良地質(zhì)體寬度對(duì)圍巖穩(wěn)定性的影響。計(jì)算模型邊界取開(kāi)挖高度的3倍,模型四周采用法向約束,屈服準(zhǔn)則采用莫爾—庫(kù)倫準(zhǔn)則,硐室群開(kāi)挖體模型單元示意見(jiàn)圖2,計(jì)算分析方案見(jiàn)表2。

圖2 開(kāi)挖體單元模型示意Fig.2 Schematic diagram of excavated body unit model

表2 計(jì)算分析方案Table 2 Calculation and analysis plans

3.2 開(kāi)挖分期組合

根據(jù)地下硐室群的布置,結(jié)合硐室群岔洞的結(jié)構(gòu)特點(diǎn)和施工條件,對(duì)地下硐室群進(jìn)行分層開(kāi)挖,分期開(kāi)挖組合見(jiàn)表3,分層高度示意見(jiàn)圖3。

表3 洞室群分期開(kāi)挖順序Table 3 Staged excavation sequence of cavern group

圖3 地下硐室群分層高度示意(標(biāo)高單位:m)Fig.3 Schematic diagram of the layered height of the underground chamber group

3.3 計(jì)算結(jié)果與分析

3.3.1 硐室群塑性區(qū)分布

根據(jù)計(jì)算分析方案,由三維數(shù)值模擬計(jì)算可得出,圍巖塑形應(yīng)變見(jiàn)表4,上下游邊墻和頂拱塑性區(qū)分布最大深度與PH條帶寬度間的關(guān)系見(jiàn)圖4,部分計(jì)算分析方案的塑性區(qū)分布云圖見(jiàn)圖5。

表4 圍巖最大塑性區(qū)深度Table 4 Maximum plastic zone depth of surrounding rock

從圖4、圖5可看出,PH條帶寬度越大,圍巖塑性區(qū)變化越大,由于硐室間的相互作用影響,隨著地下硐室下部開(kāi)挖,邊墻的塑性區(qū)不斷向深部發(fā)展,當(dāng)全部圍巖類(lèi)別為Ⅲ類(lèi)時(shí),主廠房的下游邊墻與主變室的上游邊墻之間的塑性區(qū)并未全部貫通。當(dāng)有不良地質(zhì)體存在時(shí),硐室開(kāi)挖結(jié)束后,主廠房下游邊墻與主變室的上游邊墻間的塑性區(qū)基本貫通,塑性區(qū)最大深度為40 m,此時(shí)硐室群圍巖極不穩(wěn)定。

圖4 硐室群塑性區(qū)深度與PH條帶寬度間的關(guān)系Fig.4 The relationship between the plastic zone depth of chamber group and PH strip width

圖5 圍巖塑性區(qū)分布云圖Fig.5 Cloud map of the distribution of plastic zone of surrounding rock

總體上,隨著開(kāi)挖次序增加,頂拱塑性區(qū)深度變化不大,基本處于穩(wěn)定狀態(tài),主廠房的上游邊墻塑性區(qū)分布最大深度為34.1 m,主變室下游邊墻塑性區(qū)分布最大深度為24 m。

3.3.2 圍巖變形

隨著不良地質(zhì)體PH條帶寬度的增加,硐室群上游邊墻的變形要比下游邊墻的變形大,這主要受地應(yīng)力傾向和相應(yīng)洞室開(kāi)挖的影響,部分開(kāi)挖方案的圍巖位移變形圖見(jiàn)圖6。硐室位移相對(duì)收斂變化率如圖7所示,參照國(guó)家標(biāo)準(zhǔn)《巖土錨桿與噴射混凝土支護(hù)工程技術(shù)規(guī)范》(GB 50086—2015)[13],Ⅲ類(lèi)圍巖允許相對(duì)位移收斂率在0.2%～0.50%,Ⅳ類(lèi)圍巖允許相對(duì)位移收斂率在0.40%～1.20%。通過(guò)圖7可看出,隨著PH條帶寬度的增加,圍巖相對(duì)位移值隨之增大,當(dāng)PH條帶寬度超過(guò)15 m時(shí),圍巖允許相對(duì)位移收斂率超過(guò)1.2%,由此得出圍巖變形程度較大,必須在一定的支護(hù)措施下控制圍巖的穩(wěn)定性。

圖6 硐室開(kāi)挖圍巖變形圖Fig.6 Deformation diagram of surrounding rock in chamber excavation

圖7 圍巖位移相對(duì)收斂變化率Fig.7 Relative convergence rate of change of surrounding rock displacement

4 巖體變形控制措施及效應(yīng)分析

4.1 支護(hù)參數(shù)設(shè)計(jì)

在第3節(jié)毛洞開(kāi)挖圍巖塑性區(qū)和變形分析基礎(chǔ)上,依據(jù)國(guó)內(nèi)外已建的相關(guān)工程實(shí)例和現(xiàn)行規(guī)范《巖土錨桿與噴射混凝土支護(hù)工程技術(shù)規(guī)范》(GB 50086—2015)[13]的基礎(chǔ)上,結(jié)合該地下洞室群工程的地質(zhì)條件及樞紐特點(diǎn),并參照最終確定的支護(hù)設(shè)計(jì)方案如表5。限于篇幅,只對(duì)PH條帶寬15m時(shí)的硐室群進(jìn)行支護(hù)研究。

表5 地下硐室群系統(tǒng)支護(hù)參數(shù)Table 5 Supporting parameters of underground chamber group system

4.2 計(jì)算結(jié)果與分析

(1)計(jì)算方法及模型。計(jì)算模型在原無(wú)支護(hù)圍巖穩(wěn)定計(jì)算模型基礎(chǔ)上增加錨桿、錨索及噴混凝土等支護(hù)措施,施工順序與無(wú)支護(hù)時(shí)相同,圍巖系統(tǒng)支護(hù)見(jiàn)圖8。

圖8 硐室群圍巖系統(tǒng)支護(hù)Fig.8 Supporting diagram of surrounding rock system of cave room group

(2)圍巖塑性區(qū)分布。開(kāi)挖完成后的有無(wú)支護(hù)圍巖塑性區(qū)最大深度比較見(jiàn)圖9,有支護(hù)情況下的塑性區(qū)分布規(guī)律與無(wú)支護(hù)基本相同。硐室群圍巖支護(hù)后,主廠房上游邊墻塑性區(qū)分布最大深度為4 m,主變室下游邊墻塑性區(qū)分布最大深度為1.4 m,塑性區(qū)最大深度均未超過(guò)錨桿長(zhǎng)度的70%,有效抑制了硐周破壞區(qū)發(fā)展。從計(jì)算結(jié)果分析,錨固支護(hù)有效縮減了塑性區(qū)深度變化,支護(hù)效果顯著,從支護(hù)后塑性區(qū)的分布情況看,各硐室圍巖穩(wěn)定性得到安全保證。

圖9 硐室群開(kāi)挖圍巖支護(hù)前后塑性區(qū)變化Fig.9 Changes of plastic zone before and after support of the cavern group surrounding rock excavation

(3)圍巖位移分布。為進(jìn)一步對(duì)錨固效果進(jìn)行對(duì)比,對(duì)有、無(wú)支護(hù)時(shí)的圍巖位移變化進(jìn)行分析統(tǒng)計(jì),見(jiàn)圖10。硐室群支護(hù)后,圍巖最大位移變形分布在上游墻,主廠房圍巖最大位移為360 mm,主變室圍巖最大位移為210 mm,與無(wú)支護(hù)相比,有支護(hù)時(shí)的圍巖最大位移明顯小于無(wú)支護(hù)時(shí)圍巖最大位移,硐周?chē)鷰r允許相對(duì)位移值為1.0%,說(shuō)明錨固支護(hù)有效抑制了硐室群邊墻、頂拱的變形。

圖10 硐室群開(kāi)挖圍巖支護(hù)前后位移變形Fig.10 Displacement and deformation before and after support of surrounding rock excavation of chamber group

5 結(jié) 論

以某大型地下硐室群巖體開(kāi)挖為工程背景,通過(guò)前期的地質(zhì)勘探得來(lái)的地質(zhì)資料,利用三維有限元數(shù)值模擬分析軟件MIDAS-GTS/NX,基于莫爾—庫(kù)倫準(zhǔn)則,對(duì)硐室群圍巖穩(wěn)定進(jìn)行彈塑性分析,采用不同的模擬開(kāi)挖方案,對(duì)圍巖穩(wěn)定進(jìn)行研究分析,得出了以下結(jié)論:

(1)硐室群開(kāi)挖完成后,主廠房的下游邊墻與主變室的上游邊墻塑性區(qū)基本貫通,延伸長(zhǎng)度為40 m,但連接的范圍較小;由于邊墻效應(yīng),隨著開(kāi)挖次序增加,硐室群的邊墻變形要比頂拱變形大。

(2)通過(guò)對(duì)不同開(kāi)挖方案的計(jì)算分析,可得出圍巖穩(wěn)定性與不良地質(zhì)體寬度有關(guān)系,當(dāng)不良地質(zhì)體的寬度越大時(shí),圍巖變形越大,圍巖越不穩(wěn)定。

(3)采用不同的支護(hù)措施,根據(jù)地下硐室群圍巖特征,通過(guò)改變錨索粘結(jié)力計(jì)算可得出,不同的錨索粘結(jié)力對(duì)圍巖穩(wěn)定性影響基本相同,有支護(hù)后的圍巖變形要比無(wú)支護(hù)時(shí)的圍巖變形要小?？傊?通過(guò)系統(tǒng)地采取支護(hù)措施,有效減少了圍巖的變形,保證了硐室群開(kāi)挖過(guò)程中的圍巖穩(wěn)定要求,可為后續(xù)的施工提供安全可靠的技術(shù)支撐。