基于二次正交組合試驗(yàn)的圍巖松動(dòng)圈影響因素分析

劉亞?wèn)|

?

基于二次正交組合試驗(yàn)的圍巖松動(dòng)圈影響因素分析

劉亞?wèn)|

(中鐵十九局集團(tuán)公司云桂鐵路四標(biāo)段項(xiàng)目部, 廣西南寧, 530012)

圍巖松動(dòng)圈的形成實(shí)質(zhì)是松動(dòng)圈內(nèi)圍巖產(chǎn)生碎脹變形, 碎脹變形是圍巖卸荷擾動(dòng)產(chǎn)生的松動(dòng)與破壞, 其影響因素包括地下水含量、地應(yīng)力以及其它地質(zhì)力學(xué)指標(biāo)如霍克布朗參數(shù)等。采用二次正交組合試驗(yàn)的4因素5水平試驗(yàn)方案, 并結(jié)合地質(zhì)雷達(dá), 按照組合設(shè)計(jì)表進(jìn)行了17次野外試驗(yàn)。經(jīng)過(guò)計(jì)算得出影響因素的顯著性由高至低分別為: 霍克布朗參數(shù)、、和的交互作用及隧道埋深。

霍克布朗參數(shù); 圍巖松動(dòng)圈; 二次正交組合試驗(yàn); 隧道

圍巖松動(dòng)圈理論是分析隧道穩(wěn)定性的一個(gè)重要手段。圍巖松動(dòng)圈的分類(lèi)方法和所確定的支護(hù)形式和參數(shù)與現(xiàn)場(chǎng)實(shí)際符合, 取得了顯著的技術(shù)、經(jīng)濟(jì)與社會(huì)效益。工程實(shí)踐證明[1–2], 該理論抓住了支護(hù)的主要對(duì)象, 即松動(dòng)圈的形成實(shí)質(zhì)是松動(dòng)圈內(nèi)圍巖產(chǎn)生碎脹變形。然而松動(dòng)圈厚度值多使用超聲波測(cè)試技術(shù)等手段獲得, 不僅費(fèi)時(shí), 而且費(fèi)用高昂, 限制了松動(dòng)圈理論的進(jìn)一步推廣應(yīng)用。因此, 研究如何快速準(zhǔn)確地獲取松動(dòng)圈厚度值是勢(shì)在必行的工作。楊志法[3]等基于反分析原理, 由實(shí)際監(jiān)測(cè)值反饋, 根據(jù)開(kāi)挖的推進(jìn)與監(jiān)測(cè)應(yīng)力和位移的數(shù)據(jù)變化, 進(jìn)一步驗(yàn)證了圍巖松動(dòng)圈的實(shí)際存在。高瑋等[4]運(yùn)用人工神經(jīng)網(wǎng)絡(luò)來(lái)預(yù)測(cè)圍巖松動(dòng)圈的厚度, 避免了傳統(tǒng)神經(jīng)網(wǎng)絡(luò)存在網(wǎng)絡(luò)結(jié)構(gòu)及計(jì)算參數(shù)等不易確定的問(wèn)題, 但沒(méi)有研究網(wǎng)絡(luò)權(quán)值進(jìn)化的優(yōu)化問(wèn)題。李順才等[5]采用有限元分析軟件ANSYS對(duì)不同地應(yīng)力、圍巖強(qiáng)度和跨度的巷道進(jìn)行了數(shù)值模擬, 得出最顯著的影響因素是地應(yīng)力, 其次為圍巖強(qiáng)度, 其跨度影響最小, 但沒(méi)有考慮地質(zhì)力學(xué)指標(biāo); 陳慶發(fā)等[6]采用聲波對(duì)不同采深等變量相異的巷道進(jìn)行了實(shí)測(cè), 得出圍巖松動(dòng)圈不同的影響因素, 并給出了相對(duì)應(yīng)的經(jīng)驗(yàn)公式; 陳秋紅等[7–8]采用FLAC3D有限元分析軟件, 通過(guò)模擬計(jì)算黏聚力強(qiáng)度()、摩擦角()、圍巖壓力()和洞室跨高比()這4個(gè)不同條件下得出的不同松動(dòng)圈厚度值, 并用回歸法建立圍巖松動(dòng)圈=(,,,)的關(guān)系函數(shù), 但未考慮地下水含量等因素對(duì)圍巖松動(dòng)圈的影響。陳建功等[9]采用動(dòng)靜力學(xué)分析得出圍巖松動(dòng)圈的半徑; 李政林[10]提出基于損傷理論來(lái)確定隧道圍巖松動(dòng)圈的大小。已有的研究很少考慮巖體賦存的地質(zhì)情況以及開(kāi)挖前后被擾動(dòng)程度。地下工程由于工程區(qū)域的地質(zhì)構(gòu)造、圍巖產(chǎn)狀和巖性、巖體力學(xué)參數(shù)、地應(yīng)力場(chǎng)構(gòu)成、地下水分布以及開(kāi)挖支護(hù)等不確定因素的存在, 使得每一個(gè)地下工程都具有特殊性。松動(dòng)圈的支護(hù)對(duì)象是松動(dòng)圈內(nèi)圍巖的碎脹變形, 碎脹變形實(shí)質(zhì)就是圍巖卸荷擾動(dòng)產(chǎn)生的松動(dòng)與破壞, 其影響因素包括地下水含量、巖體強(qiáng)度、地應(yīng)力以及其它地質(zhì)力學(xué)指標(biāo)。而Hoek-Brown準(zhǔn)則中的參數(shù)和全面考慮了巖體強(qiáng)度、巖體風(fēng)化程度、節(jié)理間距和巖體擾動(dòng)程度等一系列因素[11–14], 因而可以考慮采用Hoek-Brown準(zhǔn)則中的參數(shù)進(jìn)行松動(dòng)圈的影響因素分析。松動(dòng)圈的影響因素可能是單個(gè)或多個(gè)因素, 也可能是兩兩因素之間的交互作用。二次正交組合設(shè)計(jì)實(shí)驗(yàn)考慮因素間交互作用, 具有試驗(yàn)次數(shù)少、結(jié)果可靠的特點(diǎn)[15]。本文采用二次正交組合試驗(yàn), 考慮地下水含量、隧道埋深(表征地應(yīng)力)和Hoek-Brown參數(shù)來(lái)分析圍巖松動(dòng)圈的影響因素。

1 二次回歸正交組合設(shè)計(jì)試驗(yàn)因素的確定

董方庭[1]提出的圍巖松動(dòng)圈支護(hù)理論, 認(rèn)為影響圍巖松動(dòng)圈的因素有原巖應(yīng)力和圍巖強(qiáng)度, 而隧道的幾何尺寸影響不大, 但很多學(xué)者認(rèn)為圍巖松動(dòng)圈理論的影響因素還存在爆破震動(dòng)、節(jié)理裂隙發(fā)育程度和施工方法等。董方庭建立了方程(,),為松動(dòng)圈厚度,為地應(yīng)力,為圍巖強(qiáng)度。靖洪文等[2]采用實(shí)測(cè)數(shù)據(jù)和最小二乘法擬合得到

= 1.293(/-0.379)1/2。 (1)

式(1)表明: 在巖石強(qiáng)度不變情況下, 原巖應(yīng)力與松動(dòng)圈厚度正相關(guān); 在原巖應(yīng)力不變時(shí), 巖石強(qiáng)度與松動(dòng)圈厚度值反相關(guān)。閆長(zhǎng)斌等[15]對(duì)地下洞室開(kāi)挖爆破對(duì)圍巖松動(dòng)圈的影響進(jìn)行了數(shù)值分析, 得出判定巖體開(kāi)始松動(dòng)的損傷系數(shù)閾值推導(dǎo)公式, 說(shuō)明圍巖松動(dòng)圈同樣受到爆破震動(dòng)的影響。圍巖卸荷后, 隧道周?chē)膽?yīng)力釋放, 導(dǎo)致局部應(yīng)力集中, 圍巖強(qiáng)度隨之降低, 從而由臨空面向里形成一定范圍的松動(dòng)圈, 因而圍巖松動(dòng)圈必須考慮地應(yīng)力這一因素。對(duì)于以自重應(yīng)力為主的地下工程, 以埋深為參考。圍巖強(qiáng)度大小可由圍巖的力學(xué)參數(shù)得到, 因而圍巖的破壞過(guò)程與力學(xué)參數(shù)相關(guān), 所以巖體力學(xué)參數(shù)是圍巖松動(dòng)圈影響因素之一。而地下水的賦存多少也影響圍巖破壞程度, 地下水與巖體之間的交互作用, 相互改變著兩者的物理、化學(xué)組分及力學(xué)性質(zhì), 因此地下水的賦存情況也是圍巖松動(dòng)圈考慮因素之一。圍巖無(wú)論采用爆破開(kāi)挖或機(jī)械開(kāi)挖方式, 圍巖都會(huì)受到擾動(dòng)。Hoek-Brown準(zhǔn)則是Hoek E和Brown E T在修正Griffith理論的基礎(chǔ)上, 并通過(guò)對(duì)大量巖石三軸試驗(yàn)及巖體現(xiàn)場(chǎng)試驗(yàn)成果的統(tǒng)計(jì)分析, 采用試錯(cuò)法導(dǎo)出巖體破壞時(shí)極限主應(yīng)力和巖塊之間的關(guān)系式, 其成功之處在于將巖體強(qiáng)度與地質(zhì)描述聯(lián)系在了一起。因此, 可由Hoek-Brown參數(shù),來(lái)表征圍巖擾動(dòng)程度大小和地質(zhì)強(qiáng)度變化, 地質(zhì)強(qiáng)度指標(biāo)GSI和擾動(dòng)參數(shù)相當(dāng)于M-C準(zhǔn)則中的黏聚力強(qiáng)度()、摩擦角()。綜上所述, 考慮以地下水含量、埋深、Hoek-Brown參數(shù)和4個(gè)因素作為表征圍巖松動(dòng)圈的影響因素。

試驗(yàn)時(shí), 可參考施工現(xiàn)場(chǎng)超前地質(zhì)預(yù)報(bào)資料和紅外探水方法獲取所測(cè)隧道的含水情況, 地下水條件水平值可參考表1。Hoek-Brown準(zhǔn)則以室內(nèi)巖石力學(xué)試驗(yàn)為基礎(chǔ), 綜合考慮巖體中節(jié)理風(fēng)化程度、節(jié)理間距、巖體擾動(dòng)程度等。Hoek-Brown參數(shù)和取值可參考巖體力學(xué)與工程[16]。

表1 地下水條件參考

2 二次回歸正交組合設(shè)計(jì)

2.1 回歸方程的建立

計(jì)算回歸系數(shù), 建立含規(guī)范變量的回歸方程?；貧w系數(shù)的計(jì)算公式為: 常數(shù)項(xiàng); 一次項(xiàng)偏回歸系數(shù); 交互項(xiàng)偏回歸系數(shù);二次項(xiàng)偏回歸系數(shù)。

2.2 回歸方程顯著性檢驗(yàn)

3 工程概況及二次正交試驗(yàn)

3.1 工程簡(jiǎn)介

長(zhǎng)嶺崗隧道是云桂鐵路的重要隧道之一, 隧道處于云貴高原, 喀斯特地貌廣泛分布, 存在地形起伏大、地層繁多、巖性紛雜、斷裂構(gòu)造等工程地質(zhì)問(wèn)題。根據(jù)二次正交法, 得到如表2所示的各影響因素。和分別代表地下水含量、埋深。根據(jù)文獻(xiàn)[11], 得到表2。其中因變量為所測(cè)得圍巖松動(dòng)圈厚度值, 用雷達(dá)測(cè)得。與表2中試驗(yàn)號(hào)對(duì)應(yīng)的分別為1= 0.84,2= 1.43,3= 1.65,4= 1.58,5= 1.78,6= 1.47,7= 0.52,8= 2.11,9= 0.53,10= 1.22,11= 1.37,12= 2.14,13= 0.66,14= 1.64,15= 0.74,16= 1.68,17= 1.46。

表2 4因素二次回歸正交組合設(shè)計(jì)計(jì)算表

3.2 試驗(yàn)分析

3.2.1 因素水平確定

由因素?cái)?shù)= 4, 取0= 1, 根據(jù)星號(hào)臂長(zhǎng)的計(jì)算公式得= 1.353。

依據(jù)確定的圍巖松動(dòng)圈影響因素, 每10 m長(zhǎng)的隧道涌量(1)為10～100 L/min, 上限可確定為1γ= 100, 下限為1-= 10, 所以零水平10= 55。計(jì)算得到變化間距1= (100-55)/1.353 = 33.259, 同理可得、和的編碼, 如表3所示。

表3 因素水平的編碼

3.2.2 試驗(yàn)方案的實(shí)施

因變量為圍巖松動(dòng)圈的厚度值。采用LTD-2200探地雷達(dá)對(duì)隧道進(jìn)行探測(cè), 考慮分辨率與探測(cè)深度, 采用900 MHz可滿(mǎn)足要求。根據(jù)巖石力學(xué)參數(shù)計(jì)算松動(dòng)圈值與實(shí)測(cè)松動(dòng)范圍的對(duì)比結(jié)果見(jiàn)表2。部分處理后雷達(dá)波形圖見(jiàn)圖1, 可解譯讀出深度標(biāo)尺即松動(dòng)圈大小。

圖1 雷達(dá)現(xiàn)場(chǎng)測(cè)試及波形圖(所標(biāo)記白線(xiàn)為松動(dòng)圈實(shí)測(cè)值)

3.2.3 回歸方程的建立

3.2.4 顯著性判斷

表4 方差分析表

注: *越多代表顯著性越高, 反之越低;0.05(1, 3) = 10.13,0.01(1, 3) = 34.12,0.01(14, 3) = 25.93。

4 結(jié)論

(1)二次正交組合設(shè)計(jì)試驗(yàn)與正交試驗(yàn)相比較, 能有效減少試驗(yàn)次數(shù), 考慮了試驗(yàn)因素之間的交互作用, 得到的數(shù)據(jù)也更為全面。

(2) 建立了松動(dòng)圈影響因素的4因素5水平組合設(shè)計(jì)試驗(yàn)方案。4因素分別為埋深、地下水含量、霍克布朗參數(shù)和霍克布朗參數(shù); 5水平即上星號(hào)水平、上水平、零水平、下水平和下星號(hào)水平及各因素和水平對(duì)應(yīng)的取值。

(3) 對(duì)于松動(dòng)圈的影響因素, 以往學(xué)者只是定性分析單個(gè)因素, 忽略了兩兩因素的交互作用對(duì)松動(dòng)圈大小的影響。確定了圍巖松動(dòng)圈的影響因素試驗(yàn)方案, 并建立了回歸方程, 再由回歸方程的顯著性判別影響因素的重要程度。結(jié)果表明, 影響圍巖松動(dòng)圈的顯著性大小依次為:霍克布朗參數(shù)、、和的交互作用及隧道埋深。

參考文獻(xiàn)：

[1] 董方庭, 宋宏偉, 郭志宏, 等. 圍巖松動(dòng)圈巷道支護(hù)理論的基本觀(guān)點(diǎn)[J]. 建井技術(shù), 1994(4–5): 3–9.

[2] 郭志宏, 周希圣, 何勝春, 等. 巷道支護(hù)的圍巖松動(dòng)圈分類(lèi)方法[J]. 建井技術(shù), 1994(Z1): 11–13.

[3] 楊志法, 熊順成, 王存玉, 等. 關(guān)于位移反分析法的某些考慮[J]. 巖石力學(xué)與工程學(xué)報(bào), 1999, 18(l): 70–74.

[4] 高瑋, 鄭穎人. 巷道圍巖松動(dòng)圈預(yù)測(cè)的進(jìn)化神經(jīng)網(wǎng)絡(luò)法[J]. 巖石力學(xué)與工程學(xué)報(bào), 2002, 21(5): 658–661.

[5] 劉剛, 宋宏偉. 圍巖松動(dòng)圈影響因素的數(shù)值模擬[J]. 礦冶工程, 2003, 23(1): 1–3.

[6] 陳慶發(fā), 周科平, 胡建華,等. 緩傾薄層弱結(jié)構(gòu)松動(dòng)圈聲波測(cè)試時(shí)測(cè)孔布置的理論依據(jù)與驗(yàn)證[J]. 中南大學(xué)學(xué)報(bào): 自然科學(xué)版, 2009, 41(5): 1 406–1 410.

[7] 陳秋紅, 李仲奎, 張志增. 松動(dòng)圈分區(qū)模型及其在地下工程反饋分析中的應(yīng)用[J]. 巖石力學(xué)與工程學(xué)報(bào), 2010, 29(Z1): 3 216–3 220.

[8] 蔣樹(shù)屏, 蔡志偉, 林志,等. 考慮松動(dòng)圈的卡爾曼濾波與有限元耦合反分析法及其在圍巖穩(wěn)定性分析中的應(yīng)用[J]. 巖土力學(xué), 2009, 30(08): 2 529–2 534.

[9] 陳建功, 賀虎, 張永興. 隧道圍巖松動(dòng)圈形成機(jī)理的動(dòng)靜力學(xué)解析[J]. 巖土工程學(xué)報(bào), 2011, 33(12): 1 964–1 967.

[10] 李政林. 基于損傷理論的隧道圍巖松動(dòng)圈確定方法[J]. 地下空間與工程學(xué)報(bào), 2011, 7(6): 1 060–1 064.

[11] Hoek E, Carranza-Torres, Corkum. Hoek-Brown failure criterion [C]//Proe of the 5th North American Rock Mechanics Symposium and 17th Tunneling Association of Canada Conference. Toronto: University of Toronto Press, 2002: 267–271.

[12] Sonmez H, Ulusay R. Modifications to the geological strength index (GSI) and their applicability to stability of slopes assessment of the in-situ shear strength of rock masses and discontinuities [J]. International Journal of Rock Mechanics and Mining Sciences, 1999, 36(5): 743–760.

[13] 閆長(zhǎng)斌, 徐國(guó)元. 對(duì)Hoek-Brown公式的改進(jìn)及其工程應(yīng)用[J]. 巖石力學(xué)與工程學(xué)報(bào), 2005, 24(22): 4 030–4 035.

[14] 閆長(zhǎng)斌, 徐國(guó)元, 楊飛. 爆破動(dòng)荷載作用下圍巖累積損傷效應(yīng)聲波測(cè)試研究[J]. 巖土工程學(xué)報(bào), 2007, 29(1): 88–93.

[15] 李云雁, 胡傳榮. 試驗(yàn)設(shè)計(jì)與數(shù)據(jù)處理[M]. 北京: 化學(xué)工業(yè)出版社, 2005.

[16] 工程巖體分級(jí)標(biāo)準(zhǔn): GB50218-1994 [S].

(責(zé)任編校: 劉剛毅)

Analysis of influence factors for the broken rock zone based on quadratic orthogonal method

Liu Yadong

(Yungui Railway Department,China Railway 19th Bureau Group Co Ltd, Nanning 530012, China)

The essence of the broken rock zone is the formation of the hulking deformation of surrounding rock, while the hulking deformation of surrounding rock is produced by loosening or damaging after the unloading of tunnels, the influence factors include not only the groundwater and the ground stress, but also the other geological mechanic indicators, such as Hoek-Brown parameters and so on. The quadratic orthogonal combination experiment with 4 factors 5 level scheme is used, and with the help of Ground Penetrating Radar, 17 times of the field tests are carried out in accordance with the combination design table. The tests results from high to low are listed as follows: the Hoek-Brown parameter,, interaction ofandand the buried depth of tunnels.

Hoek-Brown parameters; the broken rock zone; quadratic orthogonal method; tunnel

10.3969/j.issn.1672–6146.2015.03.017

U 458

1672–6146(2015)03–0073–05

劉亞?wèn)|, xiuxiu128@qq.com。

2015–05–08

國(guó)家自然科學(xué)基金(51478065)。