來(lái)興平 單鵬飛 任奮華 曹建濤 苗勝軍 張?zhí)厝A

(1.西安科技大學(xué)能源學(xué)院;2.教育部西部礦井開(kāi)采及災(zāi)害防治重點(diǎn)實(shí)驗(yàn)室;3.北京科技大學(xué)土木與環(huán)境工程學(xué)院)

礦山高陡邊坡穩(wěn)定性分析包括工程地質(zhì)分析、極限平衡法、極限分析法、數(shù)值模擬和物理相似模擬等方法［1-2］。相比而言，作為實(shí)驗(yàn)巖石力學(xué)主要研究手段之一的物理相似模擬實(shí)驗(yàn)在研究邊坡穩(wěn)定性問(wèn)題中優(yōu)勢(shì)獨(dú)特［3-8］，尤其是隨著邊坡逐步向深部發(fā)展，所形成的高陡邊坡極易誘發(fā)大面積動(dòng)力學(xué)失穩(wěn)情況下，三維物理相似模擬實(shí)驗(yàn)平臺(tái)的科學(xué)構(gòu)建能夠準(zhǔn)確獲取邊坡與采場(chǎng)圍巖變形破裂特征，服務(wù)于高陡邊坡動(dòng)力失穩(wěn)災(zāi)害預(yù)報(bào)及開(kāi)挖的科學(xué)設(shè)計(jì)。

加強(qiáng)對(duì)模型材料力學(xué)性能的認(rèn)識(shí)，實(shí)現(xiàn)模型材料配比定量化選擇，有利于提高物理相似模擬實(shí)驗(yàn)的準(zhǔn)確性［9］。在物理相似模擬實(shí)驗(yàn)發(fā)展的幾十年間，國(guó)內(nèi)外的研究機(jī)構(gòu)相繼開(kāi)發(fā)了多種配比的模型材料，它們具有不同的力學(xué)性能，并且在實(shí)際應(yīng)用中都取得了一定的成果。例如意大利貝加莫結(jié)構(gòu)與模型實(shí)驗(yàn)研究所［10］采用氧化鉛模擬巖石基本成分，張強(qiáng)勇［11-12］、李仲奎［13］和韓伯鯉［14］分別研制 NIOS、IBSCM和MIB模型材料，研究手段是制作標(biāo)準(zhǔn)巖樣并獲取力學(xué)參數(shù)，模型具有與原型極為相似的力學(xué)性能［15］。

水廠鐵礦北山采區(qū)西部9～33勘探線露天開(kāi)采境界線將向西擴(kuò)展1 400 m，露天開(kāi)采底標(biāo)高將從原開(kāi)采設(shè)計(jì)的－350 m降至－440 m，巖體最大強(qiáng)度達(dá)160 MPa。三維物理相似模擬實(shí)驗(yàn)(3.6 m(長(zhǎng))×2.0 m(寬)×2.5 m(高))的幾何相似比為1∶400，為小比例模擬，模擬難度大［16］。但上述各種模型材料有關(guān)小比例變質(zhì)巖為主的礦山高陡邊坡的模型材料定量化研究相對(duì)較少。為了克服上述問(wèn)題，西安科技大學(xué)教育部西部礦井開(kāi)采及災(zāi)害防治重點(diǎn)實(shí)驗(yàn)室研制了PRCM模型材料，以水廠鐵礦北山采區(qū)西部9～33勘探線間高陡邊坡穩(wěn)定性研究為工程背景，借助標(biāo)準(zhǔn)巖樣制備、力學(xué)參數(shù)測(cè)試，得出了PRCM模型材料配比與單軸抗壓強(qiáng)度間的定量化關(guān)系及水對(duì)PRCM模型材料力學(xué)性能的影響，驗(yàn)證了PRCM模型材料模擬礦山高陡邊坡變質(zhì)巖體的可行性，為物理模擬實(shí)驗(yàn)提供了科學(xué)依據(jù)。

1 PRCM模型材料

1.1 PRCM模型材料的制備

PRCM模型材料(多孔復(fù)合介質(zhì)類巖材料，Porous Rock－like Composite Material)選擇河砂作為骨料(粒徑為0.5～1.0 mm占3.9%，0.25～0.5 mm占29.1%，0.1～0.25 mm占67.0%)，應(yīng)視為多孔松散顆粒介質(zhì)材料;石膏(乙級(jí))和大白粉均為膠結(jié)劑;水為塑性影響劑，且對(duì)實(shí)現(xiàn)配比定量化有重要影響。通過(guò)增減骨料控制PRCM模型材料的力學(xué)強(qiáng)度，大白粉與石膏相對(duì)比例不變(均為5個(gè)單位)。本次測(cè)試共8組，第1組材料配比為河砂∶石膏∶大白粉(下同)=50∶5∶5，其他組均較前1組河砂比例增加5個(gè)單位;此外為了分析水對(duì)PRCM模型材料力學(xué)性能的影響，每組均制備9個(gè)試件，加水20、30和40 mL各3個(gè)。將材料按理論值量取混合整體夯壓成型(?50 mm×100 mm)。制備完成后晾曬3 d進(jìn)行單軸壓縮測(cè)試。

1.2 水對(duì)PRCM模型材料的非線性耦合作用

水對(duì)PRCM模型材料的非線性耦合作用顯著［17-18］。松散顆粒組成的PRCM模型材料各向異性，抗壓強(qiáng)度大小取決于水對(duì)PRCM模型材料的非線性耦合作用:一方面是水與各組分間的相互耦合，另一方面水對(duì)材料的非線性作用(σa為有效應(yīng)力)是影響材料強(qiáng)度的主因，表現(xiàn)在孔隙靜水壓力(μ)與動(dòng)水壓力(τd)的疊加作用。塑性影響劑添加量不同，使得PRCM模型材料有效應(yīng)力表現(xiàn)為以下3種可能。

(1)PRCM模型材料未達(dá)到水飽和狀態(tài):外載荷作用時(shí)，由于孔隙水無(wú)法提供有效內(nèi)部應(yīng)力(σ)，使得內(nèi)部微裂隙產(chǎn)生并發(fā)生擴(kuò)容現(xiàn)象降低材料強(qiáng)度;但此時(shí)水的增多會(huì)導(dǎo)致靜水壓力無(wú)益性逐漸減小(公式(1))，此時(shí)PRCM模型材料強(qiáng)度處于“穩(wěn)中有升”，PRCM模型材料的有效應(yīng)力處于穩(wěn)增變化階段。

(2)PRCM模型材料達(dá)到水飽和狀態(tài):外載荷作用下導(dǎo)致內(nèi)部孔隙間經(jīng)歷一個(gè)排水平衡過(guò)程且結(jié)構(gòu)趨于穩(wěn)定狀態(tài)，骨架和水體共同承擔(dān)外荷載，孔隙水將對(duì)骨架施加有益靜水壓力(公式(2))逐漸達(dá)到峰值，使得PRCM模型材料的有效應(yīng)力增至最大，PRCM模型材料的有效應(yīng)力處于峰值變化階段。

(3)PRCM模型材料達(dá)到水超飽和狀態(tài):孔隙將充滿流動(dòng)水體，將對(duì)骨架施加孔隙動(dòng)水壓力(公式(3)中R、J分別為水的容重及水力坡度)產(chǎn)生切向應(yīng)力，剪應(yīng)力效應(yīng)導(dǎo)致骨架切向變形加大、微裂隙繼續(xù)擴(kuò)張甚至出現(xiàn)模型完全破壞，即出現(xiàn)動(dòng)力學(xué)失穩(wěn)現(xiàn)象，PRCM模型材料力學(xué)特性出現(xiàn)“不升反降”，有效應(yīng)力處于快速降低階段。

2 PRCM模型材料力學(xué)性能分析

2.1 單軸抗壓強(qiáng)度測(cè)試結(jié)果

PRCM模型材料單軸壓縮測(cè)試系統(tǒng)采用無(wú)測(cè)限壓力機(jī)壓縮所制備的標(biāo)準(zhǔn)巖樣［19］，測(cè)試系統(tǒng)主要由3部分組成:

(1)實(shí)時(shí)記錄單軸應(yīng)力(σc)、應(yīng)變(ε)等數(shù)據(jù)，繪制巖樣全程應(yīng)力－應(yīng)變曲線，并在應(yīng)力峰值的40%處確定PRCM模型材料的彈性模量等基本力學(xué)參數(shù)［20］;

(2)采用聲發(fā)射(AE)對(duì)單軸壓縮狀態(tài)下的標(biāo)準(zhǔn)巖樣進(jìn)行內(nèi)部損傷演化規(guī)律分析，確定應(yīng)力累積－事件總數(shù)等參數(shù)間的定量化關(guān)系;

(3)采用光學(xué)攝像(OC)對(duì)巖樣宏觀裂隙發(fā)育情況進(jìn)行描述，得出實(shí)時(shí)動(dòng)態(tài)變形信息。

測(cè)試結(jié)果表明(如圖1(a)～(h)及表1所示)，PRCM模型材料的彈性模量可在比較大的范圍內(nèi)調(diào)整，通過(guò)改變骨料組分含量，可以使PRCM模型材料的彈性模量在2～190 MPa，單軸抗壓強(qiáng)度在0.118～0.378 MPa，滿足小比例變質(zhì)巖為主的礦山高陡邊坡三維物理相似模擬實(shí)驗(yàn)要求［16］。此外，各組內(nèi)標(biāo)準(zhǔn)巖樣應(yīng)力峰值均隨著水量的增多而持續(xù)性增加，說(shuō)明材料并未達(dá)到水超飽和狀態(tài)。

當(dāng)配比=50∶5∶5時(shí)(如圖1(a)所示)，水量的改變使應(yīng)力－應(yīng)變呈現(xiàn)不同特征:水量為20 mL時(shí)試樣的峰應(yīng)力值與最大軸向變形均最小，分別為0.118 MPa與1.6%，彈性變形階段顯著;而水量為30 mL時(shí)，較前者保持應(yīng)力峰值能力明顯增加，表現(xiàn)出良好的塑性特征;當(dāng)水量進(jìn)一步加大到40 mL時(shí)，應(yīng)力峰值也隨之加大至0.301 MPa，穩(wěn)定發(fā)展階段楊氏模量明顯增加，表明內(nèi)部破壞速率加快。如圖1(b)所示，當(dāng)配比55∶5∶5時(shí)，水量為20與30mL下的試樣應(yīng)力峰值較前組均有所增加，其中水量20 mL時(shí)增量最大，而水量為40mL時(shí)的峰值保持與前組一致且在峰值附近出現(xiàn)應(yīng)力躍進(jìn)，為突發(fā)能量釋放所致;水量30 mL時(shí)，仍舊具有良好的保持峰值應(yīng)力強(qiáng)度的能力。

當(dāng)河砂比例加至60(如圖1(c))時(shí)，峰值強(qiáng)度明顯降低，水量為20與30 mL的應(yīng)力峰值與40 mL的差距進(jìn)一步減小，且均有較大的殘余強(qiáng)度。當(dāng)配比65∶5∶5時(shí)(圖1(d))，應(yīng)力峰值整體應(yīng)力值均偏小。而當(dāng)河砂加至70時(shí)(圖1(e))，峰值強(qiáng)度顯著提高，水量40 mL時(shí)的峰值應(yīng)力至0.378 MPa，為最大值，水量30 mL時(shí)破壞前應(yīng)力強(qiáng)化顯著，而水量20 mL時(shí)，破壞前的曲線所圍面積遠(yuǎn)大于峰后面積，表明在此類情況下動(dòng)力失穩(wěn)現(xiàn)象極易出現(xiàn)。隨著骨料的再次增加(圖1(f))，水量為20 mL時(shí)的峰前曲線面積遠(yuǎn)小于峰后面積，與前組恰好相反，說(shuō)明在此情況下試樣的破壞過(guò)程穩(wěn)定;當(dāng)配比為80∶5∶5與85∶5∶5時(shí)，骨料的增加使得試件的出現(xiàn)明顯的延性特征。

圖1 外載荷作用下不同配比PRCM模型材料的破壞過(guò)程

2.2 水對(duì)PRCM模型材料特性的影響

水對(duì)PRCM模型材料的力學(xué)性能影響至關(guān)重要，表現(xiàn)在以下2個(gè)方面。

(1)水與骨料、膠結(jié)劑等物質(zhì)間的相互耦合(圖2)不同于連續(xù)固體巖石材料的線性關(guān)聯(lián)，材料配比與單軸抗壓強(qiáng)度間存在明顯的非線性定量化關(guān)系，通過(guò)對(duì)測(cè)試結(jié)果擬合表明:水為20、30及40 mL時(shí)，PRCM模型材料配比與單軸抗壓強(qiáng)度間的函數(shù)關(guān)系滿足“雙馬鞍”模型，擬合后均出現(xiàn)3個(gè)極大與2個(gè)極小值，在相同條件下大體在配比70∶5∶5時(shí)材料的單軸抗壓強(qiáng)度達(dá)到最大，而在配比63∶5∶5為強(qiáng)度最小，配比53∶5∶5和80∶5∶5處單軸強(qiáng)度為第2與第3極大值位置，在配比77∶5∶5附近大體出現(xiàn)第2極小值;另一方面，水的逐漸增多使得測(cè)試結(jié)果波動(dòng)性愈大，其中加水40 mL的波動(dòng)性最大，測(cè)試范圍為0.118至0.378 MPa，可為物理相似模擬實(shí)驗(yàn)科學(xué)配比選取提供參考數(shù)據(jù)。

表1 單軸壓縮主要力學(xué)參數(shù)測(cè)試結(jié)果

圖2 不同配比下單軸抗壓強(qiáng)度變化趨勢(shì)

(2)根據(jù)單軸壓縮測(cè)試結(jié)果(圖3)，在測(cè)試范圍內(nèi)，隨著水量逐漸增加，PRCM模型材料單軸抗壓強(qiáng)度增大，均在40mL時(shí)取得極大值，表明測(cè)試各組類巖材料的有效應(yīng)力變化均未處于快速降低階段，μ仍有增大可能性，即孔隙動(dòng)水壓力τd為0;此外，不同配比下曲線變化趨勢(shì)各異，配比為60∶5∶5及70∶5∶5時(shí)曲線為線性變化，而配比為65∶5∶5與85∶5∶5為雙曲線，其他均為拋物線變化趨勢(shì)。

圖3 不同加水量單軸抗壓強(qiáng)度峰值變化趨勢(shì)

3 在礦山高陡邊坡穩(wěn)定性模擬實(shí)驗(yàn)中應(yīng)用

水廠鐵礦是首都鋼鐵公司自有礦石資源的主要礦山和重要鐵礦石基地之一，且又是遷安鐵礦區(qū)中儲(chǔ)量規(guī)模最大的鐵礦床，實(shí)際采礦規(guī)模為1 100萬(wàn)t/a。按照擴(kuò)幫工程方案設(shè)計(jì)，北山采區(qū)西部9～33勘探線長(zhǎng)約1 200 m的露天開(kāi)采境界線將向西擴(kuò)展200 m，露天開(kāi)采底標(biāo)高將從原開(kāi)采設(shè)計(jì)的－350 m降至－440 m，深度為840 m，形成深凹露天礦山邊坡［21］，將對(duì)邊坡穩(wěn)定性造成嚴(yán)重影響，并有誘發(fā)大面積邊坡動(dòng)力學(xué)失穩(wěn)的可能性。依據(jù)深凹露天礦高陡邊坡失穩(wěn)的自身特點(diǎn)與規(guī)律，進(jìn)行礦山高陡邊坡穩(wěn)定性三維物理相似模擬實(shí)驗(yàn)。水廠鐵礦礦床以變質(zhì)巖為主，通過(guò)PRCM模型材料實(shí)驗(yàn)確定主要地層包括片麻巖、花崗巖、礫巖、泥巖和堆積物等原型地層的配比［16］。圖4為鋪裝后的三維物理相似模擬模型照片。本次模擬實(shí)驗(yàn)監(jiān)測(cè)過(guò)程采用專業(yè)型高速地質(zhì)雷達(dá)、光學(xué)鉆孔攝像儀、聲發(fā)射、數(shù)字化三維立體攝影測(cè)量等，成功解決了邊坡動(dòng)力失穩(wěn)災(zāi)害預(yù)報(bào)及開(kāi)挖的科學(xué)設(shè)計(jì)等問(wèn)題。

圖4 水廠鐵礦三維物理相似模擬模型

4 結(jié)論

(1)PRCM模型材料適用于小比例變質(zhì)巖為主的礦山高陡邊坡穩(wěn)定性模擬，其彈性模量和抗壓強(qiáng)度等主要力學(xué)參數(shù)滿足物理相似模擬實(shí)驗(yàn)的要求且制備簡(jiǎn)便。

(2)水對(duì)PRCM模型材料的非線性耦合作用顯著，使得PRCM模型材料有效應(yīng)力表現(xiàn)為3種可能(PRCM模型材料未達(dá)到水飽和狀態(tài)、達(dá)到水飽和狀態(tài)、達(dá)到水超飽和狀態(tài)下的應(yīng)力表現(xiàn))。

(3)PRCM模型材料配比與應(yīng)力峰值在不同水量下為非線性定量化關(guān)系，即擬合后均表現(xiàn)出“雙馬鞍”變化趨勢(shì)，且強(qiáng)度均隨水量增多而增大且都在40 mL時(shí)取得極大。

［1］ 楊天鴻，張鋒春，于慶磊，等.露天礦高陡邊坡穩(wěn)定性研究現(xiàn)狀及發(fā)展趨勢(shì)［J］.巖土力學(xué)，2011，32(5):1437-1472.

［2］ 蔡美峰，何滿朝，劉東燕.巖石力學(xué)與地下工程［M］.北京:科學(xué)出版社，2002.

［3］ LU Ai-hong，MAO Xian-biao，LIU Hai-shun.Physical simulation of rock burst induced by stresswaves［J］.Journal of China University of Mining and Technology，2008，18(3):0401-0405.

［4］ 李天斌，王湘鋒，孟陸波.巖爆的相似材料物理模擬研究［J］.巖石力學(xué)與工程學(xué)報(bào)，2011，30(增1):2610-2616.

［5］ 王克忠，李仲奎.深埋長(zhǎng)大引水隧洞三維物理模型滲透性試驗(yàn)研究［J］.巖石力學(xué)與工程學(xué)報(bào)，2009，28(4):725-737.

［6］ 肖江，高喜才.復(fù)雜地質(zhì)條件下露天礦邊坡穩(wěn)定的相似模型研究［J］.巖石力學(xué)與工程學(xué)報(bào)，2006，25(增2):3661-3666.

［7］ 黃慶享.淺埋煤層保水開(kāi)采隔水層穩(wěn)定性的模擬研究［J］.巖石力學(xué)與工程學(xué)報(bào)，2009，28(5):988-992.

［8］ 潘一山，章夢(mèng)濤，王來(lái)貴，等.地下硐室?guī)r爆的相似材料模擬試驗(yàn)研究［J］.巖土工程學(xué)報(bào)，1997，19(4):49-56.

［9］ 林韻梅.實(shí)驗(yàn)巖石力學(xué)［M］.北京:煤炭工業(yè)出版社，1984.

［10］ 姜小蘭，陳進(jìn)，孫紹文，等.錦屏工程地質(zhì)力學(xué)模型壩基巖體新型材料研究［J］.長(zhǎng)江科學(xué)院院報(bào)，2009，26(6):40-43.

［11］ 張強(qiáng)勇，李術(shù)才，郭小紅，等.鐵晶砂膠結(jié)新型巖土相似材料的研制及其應(yīng)用［J］.巖土力學(xué)，2008，29(8):2126-2130.

［12］ 王漢鵬，李術(shù)才，張強(qiáng)勇，等.新型地質(zhì)力學(xué)模型試驗(yàn)相似材料的研制［J］.巖石力學(xué)與工程學(xué)報(bào)，2006，25(9):1842-1847.

［13］ 馬芳平，李仲奎，羅光福.NIOS模型材料及其在地質(zhì)力學(xué)相似模擬試驗(yàn)中的應(yīng)用［J］.水力發(fā)電學(xué)報(bào)，2004，23(1):48-51.

［14］ 韓伯鯉，陳霞齡，宋一樂(lè).巖體相似材料的研究［J］.武漢水利電力大學(xué)學(xué)報(bào)，1997，30(2):6-9.

［15］ 吳綿拔.加載速率對(duì)巖石抗壓和抗拉強(qiáng)度的影響［J］.巖土工程學(xué)報(bào)，1982，4(2):97-106.

［16］ 來(lái)興平，單鵬飛，任奮華，等.礦山高陡邊坡多介質(zhì)復(fù)合相似材料強(qiáng)度對(duì)比試驗(yàn)［J］.西安科技大學(xué)學(xué)報(bào)，2012，32(2):1-5.

［17］ 韓 濤，楊維好，楊志江，等.多孔介質(zhì)固液耦合相似材料的研制［J］.巖土力學(xué)，2011，32(5):1411-1417.

［18］ 邢義川，謝定義，李 振.非飽和土的應(yīng)力傳遞機(jī)理與有效應(yīng)力原理［J］.巖土工程學(xué)報(bào)，2001，23(1):53-56.

［19］ 王學(xué)濱.巖樣單軸壓縮軸向及側(cè)向變形耗散能量及穩(wěn)定性分析［J］.巖石力學(xué)與工程學(xué)報(bào)，2005，24(5):846-853.

［20］ JO S，Chandler N A.Innovative laboratory testing［J］.International Journal of Rock Mechanics and Mining Science，2004，41 (8):1427-1445.

［21］ 首鋼礦業(yè)公司地質(zhì)研究所.水廠鐵礦修改設(shè)計(jì)補(bǔ)充勘探地質(zhì)報(bào)告［R］.遷安:首鋼礦業(yè)公司地質(zhì)研究所，2009.