褚吉祥，于慶磊，曹永勝，孫 鵬

(1.安徽太平礦業(yè)有限公司； 2.東北大學(xué)巖石破裂與失穩(wěn)研究所)

引 言

矽卡巖型礦床是中國最常見的鐵礦床類型之一，其賦存形態(tài)復(fù)雜多變，礦巖軟弱破碎。矽卡巖體在形成過程中因變質(zhì)和蝕變作用，導(dǎo)致其物理力學(xué)性質(zhì)差異很大，巖石強(qiáng)度低，自穩(wěn)性差，水理性強(qiáng)，遇水易膨脹、崩解、碎化，甚至在無淋水條件下，矽卡巖體吸收空氣中的水分，也不斷膨脹、崩解[1]。長期生產(chǎn)實(shí)踐表明，矽卡巖體常是井巷變形破壞及其他形式地壓顯現(xiàn)的突發(fā)點(diǎn)和關(guān)鍵因素[2]。因此，確定矽卡巖體力學(xué)參數(shù)，對(duì)于評(píng)估巖體力學(xué)行為，指導(dǎo)該類型礦山生產(chǎn)設(shè)計(jì)、施工、維護(hù)、評(píng)價(jià)圍巖穩(wěn)定性、預(yù)防工程事故等有十分重要的意義。

矽卡巖體由于吸水易崩解的特性，難以制成標(biāo)準(zhǔn)的巖石試樣，遂無法通過試驗(yàn)手段獲取其完整的力學(xué)參數(shù)。近年來，利用圍巖各種力學(xué)響應(yīng)反推模型參數(shù)的反分析方法，成為當(dāng)前的熱點(diǎn)研究問題，目前已形成完整的分析體系。反分析方法主要包括解析法和數(shù)值法，解析法由于其局限性難以適用于復(fù)雜的計(jì)算模型。數(shù)值法因其計(jì)算效率高、求解難度小等優(yōu)點(diǎn)得到了廣泛推廣。

目前，最常用的反分析方法是位移反分析法，即將現(xiàn)場(chǎng)巖體的變形監(jiān)測(cè)數(shù)據(jù)作為反演依據(jù)，調(diào)整數(shù)值模型中的基本巖體參數(shù)，得出最優(yōu)解。張晨明等[3]利用增量位移反分析法對(duì)硐室圍巖彈性模量進(jìn)行了反演。朱珍德等[4]利用優(yōu)化后的PSO-BP隧道位移反分析系統(tǒng)，基于現(xiàn)場(chǎng)監(jiān)測(cè)的變形位移數(shù)據(jù)，對(duì)隧道圍巖參數(shù)進(jìn)行了反演分析，并通過與實(shí)測(cè)數(shù)據(jù)對(duì)比驗(yàn)證了該方法的可行性。張社榮等[5]綜合考慮大型地下硐室群施工現(xiàn)場(chǎng)的開挖進(jìn)度信息、支護(hù)進(jìn)度信息及新出露的地質(zhì)信息，提出適用于大型地下硐室群施工期圍巖力學(xué)參數(shù)實(shí)時(shí)動(dòng)態(tài)反演分析方法。劉勇等[6]基于隧洞施工過程中的動(dòng)態(tài)監(jiān)測(cè)數(shù)據(jù)，利用正交設(shè)計(jì)反分析方法對(duì)松散土層、注漿加固土層和初期支護(hù)的彈性模量進(jìn)行了反演分析。全面的測(cè)量信息是地下工程巖體力學(xué)參數(shù)反分析中的關(guān)鍵一環(huán)[7]。圍巖松動(dòng)圈厚度分布規(guī)律亦是反映擾動(dòng)條件下圍巖力學(xué)反饋的重要參數(shù)。近年來，考慮松動(dòng)圈的反分析方法已有不少學(xué)者研究，如陳秋紅等[8]建立了反映松動(dòng)圈差異性的反饋模型，并提出了數(shù)值模擬中對(duì)松動(dòng)圈的處理方法，在錦屏一級(jí)水電站反饋分析中取得了良好的效果。張飛等[9]建立了基于松動(dòng)圈范圍的參數(shù)弱化基本模型，并通過后驗(yàn)差法驗(yàn)證了該方法的適用性。張繼勛等[10]利用遺傳算法，結(jié)合松動(dòng)圈深度信息，通過反演發(fā)現(xiàn)隧道設(shè)計(jì)采用的參數(shù)取值偏安全。綜上所述，巷道力學(xué)行為是巷道圍巖在多種動(dòng)態(tài)變化因素?cái)_動(dòng)下的直觀反映，但目前巷道圍巖力學(xué)參數(shù)的反演研究大多局限于單獨(dú)采用位移監(jiān)測(cè)結(jié)果或松動(dòng)圈測(cè)量結(jié)果，反演的力學(xué)參數(shù)可能不夠全面[11]。

本文以圍巖具有蝕變特征的安徽太平礦業(yè)有限公司-480 m中段運(yùn)輸巷道為研究對(duì)象，以數(shù)值模擬為研究工具，采用黃金分割法，結(jié)合現(xiàn)場(chǎng)位移監(jiān)測(cè)結(jié)果和塑性松動(dòng)區(qū)兩類測(cè)量指標(biāo)，系統(tǒng)地反演分析了矽卡巖體力學(xué)參數(shù)，為該巷道后續(xù)的穩(wěn)定性分析提供數(shù)據(jù)支撐。

1 反演分析方法

1.1 反演參數(shù)確定

巖體力學(xué)參數(shù)主要包括變形參數(shù)(彈性模量(E)及泊松比(μ))與強(qiáng)度參數(shù)(內(nèi)聚力(C)、內(nèi)摩擦角(φ)、抗拉強(qiáng)度(σt))。理論上,當(dāng)模型中包含的反演參數(shù)越多，得到的結(jié)果就與現(xiàn)場(chǎng)巖體實(shí)際參數(shù)越符合；但過多的反演參數(shù)不僅會(huì)增加反演分析的難度，還會(huì)降低最終結(jié)果的準(zhǔn)確性。依據(jù)工程經(jīng)驗(yàn)，變形參數(shù)中彈性模量(Ε)對(duì)巖體變形的影響程度高于泊松比(μ)，強(qiáng)度參數(shù)中內(nèi)聚力(C)對(duì)巖體的影響高于內(nèi)摩擦角(φ)和抗拉強(qiáng)度(σt)。因此，本文選用巷道模型的彈性模量(E)和內(nèi)聚力(C)作為反演參數(shù)。其他影響相對(duì)較小的參數(shù)依據(jù)經(jīng)驗(yàn)參數(shù)作為已知項(xiàng)代入數(shù)值計(jì)算。在本文中，強(qiáng)度參數(shù)的反演依據(jù)現(xiàn)場(chǎng)測(cè)試圍巖松動(dòng)圈范圍，變形參數(shù)的反演則依據(jù)現(xiàn)場(chǎng)監(jiān)測(cè)的圍巖位移結(jié)果。

1.2 反演過程

為了評(píng)價(jià)反演參數(shù)的準(zhǔn)確性，構(gòu)造位移及松動(dòng)圈反分析的目標(biāo)函數(shù)來評(píng)價(jià)對(duì)比每次反演模擬的準(zhǔn)確度。設(shè)現(xiàn)場(chǎng)位移及松動(dòng)圈的實(shí)測(cè)值為X0，再令數(shù)值模擬得到的計(jì)算值Xi=(Ei，Ci，μ，φ，σt)，i=1，2，…，n，作為第i次數(shù)值模擬計(jì)算的參數(shù)值。然后將測(cè)量值與模擬計(jì)算值的差值記為:

(1)

式中：R為目標(biāo)函數(shù);X0為測(cè)量值;Xi為模擬計(jì)算值。

選取不同的計(jì)算參數(shù)則得到的目標(biāo)函數(shù)也不同。運(yùn)用黃金分割法計(jì)算位移與測(cè)量位移越接近，則R的值越小。根據(jù)礦山已有報(bào)告中巷道圍巖資料，結(jié)合測(cè)量結(jié)果，優(yōu)先確定巖體變形參數(shù)，再反演巖體強(qiáng)度參數(shù)。

2 巷道圍巖參數(shù)反演分析

2.1 工程概況

安徽太平礦業(yè)有限公司-480 m中段運(yùn)輸巷道是該礦山深部開采的明豎井-盲豎井開拓的轉(zhuǎn)運(yùn)中段，其穩(wěn)定性直接關(guān)乎深部資源的開采與利用，具有十分重要的戰(zhàn)略意義。受成礦期構(gòu)造斷裂影響，礦區(qū)內(nèi)地層變動(dòng)強(qiáng)烈，地層產(chǎn)狀變化顯著，圍巖蝕變作用(主要為矽卡巖化)顯著。蝕變后的矽卡巖體(蝕變巖)常為礦體的直接頂、底板，其形態(tài)及范圍受接觸構(gòu)造控制，暴露后，遇水易膨脹、崩解、碎化，甚至在無淋水條件下，矽卡巖體吸收空氣中的水分，也不斷崩解、膨脹、完全失去強(qiáng)度，巖體穩(wěn)固性極差，給巷道掘進(jìn)及后期支護(hù)工作帶來極大困擾。盡管巷道采用了保守的錨網(wǎng)噴支護(hù)和U形鋼拱架支護(hù)，但是巷道仍會(huì)發(fā)生底鼓、側(cè)幫膨出變形；持續(xù)的膨出變形導(dǎo)致巷道圍巖片幫塌落，進(jìn)而影響正常的采礦生產(chǎn)和礦石運(yùn)輸，甚至造成安全生產(chǎn)事故。因此，反演分析矽卡巖體巷道的力學(xué)參數(shù)對(duì)評(píng)價(jià)圍巖穩(wěn)定性、設(shè)計(jì)維護(hù)方案等具有重要意義。

2.2 巷道變形量值

巷道地層位移監(jiān)測(cè)能夠直觀反映巷道開挖及支護(hù)后圍巖的各種動(dòng)態(tài)響應(yīng)。采用多點(diǎn)位移計(jì)監(jiān)測(cè)-480 m 中段運(yùn)輸巷道1#測(cè)點(diǎn)處(見圖1)的變形情況。在巷道的左幫和拱頂分別設(shè)置埋深為1 m、2 m、3 m、4 m 的4個(gè)多點(diǎn)位移計(jì)以測(cè)量巷道圍巖位移。對(duì)各測(cè)點(diǎn)進(jìn)行為期2個(gè)月的實(shí)時(shí)監(jiān)測(cè)，監(jiān)測(cè)頻率為1 min/次，巷道圍巖的累積變形量見圖2。

圖1 測(cè)點(diǎn)位置示意圖

由圖2可知：左幫矽卡巖體受開挖擾動(dòng)影響較大，變形經(jīng)歷了30多d的急劇增長，最高變化率達(dá)1 mm/d，此后變形逐漸穩(wěn)定，最大位移為48 mm。拱頂位移雖變化較小，最大位移只有2.6 mm，但拱頂變形持續(xù)時(shí)間較長，48 d以后趨于平穩(wěn)。并且，隨孔深增加，巷道圍巖穩(wěn)定性逐漸增加，各測(cè)點(diǎn)變形量逐漸減小。

2.3 松動(dòng)圈測(cè)試

通過聲波測(cè)試得出的松動(dòng)圈大小，直接反映巖體的完整性程度，進(jìn)而更深層次地體現(xiàn)圍巖強(qiáng)度參數(shù)(內(nèi)聚力)的強(qiáng)弱。因此，松動(dòng)圈大小和分布在判定圍巖穩(wěn)定性的同時(shí)也可作為圍巖參數(shù)反演的重要依據(jù)。采用意大利FASTWAVE高速高精度探地雷達(dá)，在-480 m中段運(yùn)輸巷道處沿巷道走向進(jìn)行了掃描探測(cè)。左幫在巷道掘進(jìn)過程中發(fā)生了局部坍塌，地質(zhì)雷達(dá)探測(cè)結(jié)果顯示該段巷道松動(dòng)圈大部分區(qū)域松動(dòng)圈厚度約2.00 m，局部超過2.50 m(見圖3)。經(jīng)過統(tǒng)計(jì)計(jì)算，矽卡巖體巷道圍巖松動(dòng)圈厚度約為2.25 m。

2.4 計(jì)算模型建立

以安徽太平礦業(yè)有限公司-480 m中段運(yùn)輸巷道為研究對(duì)象，采用Flac3D軟件建立三維數(shù)值模型(見圖4)，模型尺寸為60.0 m×40.0 m×40.0 m，巷道處于模型中心，斷面尺寸為2.5 m×2.5 m。模型四周施加法向約束邊界條件，底部為全約束邊界條件，為模擬模型上部賦巖的重量，在頂部施加16.6 MPa均布載荷。模型共劃分單元數(shù)523 156個(gè)，節(jié)點(diǎn)數(shù)94 447個(gè)。采用更能反映巖體強(qiáng)度隨塑性變形變化的應(yīng)變軟化模型，反演計(jì)算巷道周圍矽卡巖體的力學(xué)參數(shù)。

圖2 測(cè)點(diǎn)位移變化曲線

圖3 巷道松動(dòng)圈范圍

2.5 反演計(jì)算分析

利用黃金分割法進(jìn)行迭代計(jì)算，根據(jù)現(xiàn)場(chǎng)地質(zhì)報(bào)告中的圍巖參數(shù)資料，確定待反演參數(shù)E、C的取值(見表1、表2)分別為0～17.53 GPa和0～4.23 MPa。

圖4 巷道三維數(shù)值模型及尺寸

計(jì)算過程中其他參數(shù)保持不變，泊松比取0.3，內(nèi)摩擦角為33°，抗拉強(qiáng)度為0.5 MPa，密度為2 630 kg/m3。當(dāng)目標(biāo)函數(shù)取值小于0.01時(shí)，即得到反演結(jié)果。由表1可知：經(jīng)過7次黃金分割，10次數(shù)值計(jì)算后，矽卡巖體的彈性模量目標(biāo)函數(shù)符合目標(biāo)值。此時(shí)，彈性模量為2.1 GPa，模擬計(jì)算得到巷道左幫位移為48.1 mm，與巷道監(jiān)測(cè)到的最終變形量(48 mm)相吻合，等于設(shè)置的目標(biāo)函數(shù)值0.01。由表2可知：當(dāng)選取的內(nèi)聚力為1.51 MPa時(shí)，模擬計(jì)算得到的松動(dòng)圈范圍為2.25 m，與實(shí)際測(cè)試深度吻合，小于目標(biāo)函數(shù)值(0.01)。綜上所述，矽卡巖體的彈性模量取2.1 GPa，內(nèi)聚力取1.51 MPa時(shí)，為該參數(shù)組合的最優(yōu)解。

3 結(jié) 論

本文聯(lián)合位移-松動(dòng)圈2種測(cè)量指標(biāo)，通過Flac3D軟件建立三維模型，基于黃金分割法，對(duì)安徽太平礦業(yè)有限公司-480 m中段運(yùn)輸巷道的力學(xué)參數(shù)進(jìn)行了反演分析。主要結(jié)論如下：

1)反演過程中選取圍巖的彈性模量與內(nèi)聚力作為反演參數(shù)，而將其他影響程度較小的參數(shù)當(dāng)成已知參數(shù)代入數(shù)值計(jì)算，可以精確計(jì)算出拱頂沉降位移及圍巖松動(dòng)圈范圍，節(jié)約迭代次數(shù)，減少計(jì)算量。

表1 巖體參數(shù)E反演結(jié)果

表2 巖體參數(shù)C反演結(jié)果

2)本次工程圍巖力學(xué)參數(shù)反演結(jié)果為：彈性模量2.1 GPa，內(nèi)聚力1.51 MPa，計(jì)算結(jié)果與現(xiàn)場(chǎng)測(cè)量數(shù)據(jù)吻合，說明了反演分析方法的可行性。

3)基于位移-松動(dòng)圈2種測(cè)量指標(biāo)的反演分析方法具有很高的可信度，本文的方法也可以用于類似工程的參數(shù)分析。