任富強，常 遠，汪 東，常力月，夏圓圓

(1.中國礦業(yè)大學(北京) a.力學與建筑工程學院；b.深部巖土力學與地下工程國家重點實驗室，北京 100083；2.成都理工大學 地球科學學院,成都 610059)

?

集寧板巖巖爆碎屑分形特征分析

任富強1a,1b，常 遠1a，汪 東2，常力月1a，夏圓圓1a,1b

(1.中國礦業(yè)大學(北京) a.力學與建筑工程學院；b.深部巖土力學與地下工程國家重點實驗室，北京 100083；2.成都理工大學 地球科學學院,成都 610059)

為了補充板巖的巖爆試驗研究，利用深部巖爆模擬試驗系統(tǒng)對內(nèi)蒙古集寧板巖進行了室內(nèi)巖爆試驗，通過三向加載—分級卸載—垂向加載來模擬板巖巖爆的發(fā)生。收集巖爆碎屑進行篩分，對粒徑>10 mm的碎屑按不同的巖塊形態(tài)的衡量指標，分別計算其分形維數(shù)。結(jié)果表明：巖塊形態(tài)的衡量指標對板巖巖爆的碎屑分維值有一定的影響，分維值分布在1.62～2.73之間，且按長度/寬度/厚度-數(shù)量計算得到的分維值離散性最小，可作為計算板巖巖爆碎屑分維值的優(yōu)選衡量指標；厚度、長度、寬度3個方向的破碎難度依次增強；長寬比、長厚比和寬厚比與數(shù)量分維值可表征碎屑的條狀、片狀、塊狀特征，板巖巖爆碎屑以條狀、條塊狀和塊狀為主，片狀特征不明顯。

板巖；室內(nèi)巖爆試驗；碎屑；不同巖塊形態(tài)；分形

1 研究背景

隨著社會經(jīng)濟的發(fā)展，資源和空間逐漸減少，為滿足生活、生產(chǎn)的需要，必然要加大對地下資源以及空間的利用。近年來深部工程越來越多而且深度越來越大，工程巖體所處的地質(zhì)壞境更加復雜，并且伴隨著地應力的增高，使得工程擾動誘發(fā)的工程地質(zhì)災害(如巖爆)問題更加嚴重和突出[1]。研究者在巖爆產(chǎn)生的機理及防治等多個方面開展了較為詳細的研究，Cook等[2-3]利用剛度理論及能量理論分析了巖爆產(chǎn)生的機制，譚以安[4]闡述了巖爆的破壞過程及防治措施，何滿潮等[5-6]利用自主研發(fā)的試驗機結(jié)合應力采集儀、聲發(fā)射等設(shè)備對巖爆的發(fā)生條件、分類、破壞特征及控制對策進行了詳細的研究。汪洋等[7]利用數(shù)值分析的方法分析了巖體開挖卸荷效應對巖爆的影響；仝躍等[8]和郭慶清等[9]分別利用統(tǒng)計學的算法對巖爆進行了預測。謝和平等[10-11]把分形理論應用于巖爆的分形特征與機制研究，得出了巖爆時刻分形維數(shù)比較低的結(jié)論，并在細觀層面利用損傷理論及分形理論對巖石的破壞特征做了闡述。此外何滿潮等[12]還對室內(nèi)巖爆試驗所得的碎屑進行了詳細的分類并總結(jié)了碎屑分形特征的研究方法。

以往的巖爆研究多以花崗巖為試驗對象，并少量為砂巖、大理巖和玄武巖，針對板巖尚未開展系統(tǒng)的巖爆研究。本文利用室內(nèi)巖爆模擬試驗得到了板巖巖爆的碎屑，通過對板巖巖爆碎屑的質(zhì)量以及尺寸(包括長度、寬度、厚度3個方向)的測量，計算了3類衡量指標的分形維數(shù)，進而獲得了板巖巖爆碎屑的分形特征。

2 巖爆試驗

2.1 試驗樣品

試驗樣品為灰色板巖標準長方體試件，板狀結(jié)構(gòu)，取自內(nèi)蒙古集寧，高為149.22 mm，寬59.53 mm，厚29.34 mm，質(zhì)量為709.30 g，密度2.680 g/cm3，X衍射分析其礦物組成為石英45%、鉀長石19.3%、斜長石16.7%、黏土礦物19%(以高嶺石為主)。試件縱波速度為5 394.11 m/s，橫波速度為2 709.88 m/s，樣品較完整，無原生裂紋，左下角處有缺陷，用硅橡膠填充，在正面、背面有豎向裂紋并貫穿于試件底部。

2.2 試驗方法及過程

利用深部巖爆模擬試驗系統(tǒng)[5]進行試驗，該系統(tǒng)包括試驗主機(可實現(xiàn)對長方柱巖石試件三向六面獨立加載，單面突然卸載)、液壓控制系統(tǒng)和數(shù)據(jù)采集系統(tǒng)3大組成部分。數(shù)據(jù)采集系統(tǒng)由動態(tài)高速壓力、位移采集儀、聲發(fā)射測試系統(tǒng)以及高速攝影系統(tǒng)組成。

本次試驗采用三向加載—分級卸載—垂向加載的方式來模擬巖爆的發(fā)生。首先對試件JNB-1#施加三向不等的初始應力來模擬巖體的原巖應力狀態(tài)，即試驗過程中采用分級加載，加載速率為0.5 MPa/s，每級應力5 MPa，加載間隔約2.5 min，達到原巖應力值后保持30 min穩(wěn)定；然后分3次逐步卸載σ3方向單面荷載，分別從5 MPa→4 MPa→3 MPa→0，卸載的同時逐級增加σ1，每卸載1次各項應力均保持15 min，最后一次卸載的同時觀察暴露面現(xiàn)象，試驗過程中加卸載應力路徑曲線如圖1所示。

圖1 JNB-1#應力路徑曲線Fig.1 Stress path curves of specimen JNB-1#

圖2 JNB-1#試驗前、后對比Fig.2 Comparison of specimen JNB-1# before and after test

2.3 試驗結(jié)果

板巖樣品的巖爆模擬試驗在86.85 min時開始卸載，90.21 min時發(fā)生破壞。剛開始卸載時應力集中約70 s有微小聲響出現(xiàn)，大約25 s后有小顆粒掉落，68 s后右上角有小顆粒彈射，23 s大塊碎屑掉落，此過程持續(xù)了16 s，大塊碎屑從試件的右上部和左上部掉落，試件發(fā)生破壞。試驗過程中的試件三向加載和卸載時的應力值及聲音描述見表1。試件JNB-1#試驗前后6個表面的形態(tài)對比如圖2所示，試件側(cè)面有一條明顯的從頂部貫穿到中底部的豎向張拉裂紋。上頂面有明顯壓碎的痕跡，下底面有碎屑掉落后的小坑洞。

3 JNB-1#巖爆碎屑分形計算

分形維數(shù)計算是對粒徑>10 mm的JNB-1#巖爆碎屑(見圖3)進行的，按3類衡量指標研究其分形特征：(Ⅰ)碎屑特征粒度(即等效邊長)；(Ⅱ)碎屑長度、寬度和厚度等3個方向的最大值；(Ⅲ)碎屑長寬比、長厚比和寬厚比。

圖3 JNB-1#板巖巖爆直徑>10 mm的碎屑Fig.3 Fragments of diameter larger than 10 mm of specimen JNB-1# rockburst

圖4 粒度-數(shù)量的分形維數(shù)Fig.4 Fractal dimension of granularity-quantity

3.1 粒度-數(shù)量的分形

將巖爆碎屑看成是正方體，根據(jù)所測碎屑的長、寬、厚度值(測量的均為碎屑長度l、寬度w、厚度h值的最大值)換算出正方體的等效邊長Leq，按式(1)計算分維值。

(1)

式中：Leq=(l×w×h)1/3；N為所選尺度內(nèi)碎屑特征粒度(等效邊長Leq)不小于Leq的碎屑數(shù)量；N0為具有最大特征粒度Leq max對應的碎屑數(shù)量；D為分維值。

當以lg(Leq max/Leq)和lgN分別為橫、縱坐標時，直線的斜率即為分形維數(shù)。

集寧板巖巖爆碎屑粒度-數(shù)量分形維數(shù)的計算結(jié)果如圖4所示，通過對散點圖的線性擬合得到擬合曲線的斜率, 即其分維值為2.26。

3.2 粒度-質(zhì)量的分形

根據(jù)單曉云等[13]的研究， 通過理論推導得到碎屑的質(zhì)量-等效邊長分形維數(shù)D按式(2)進行計算。

(2)

(3)

式中：MLeq/M為碎屑特征粒度(等效邊長Leq)小于Leq的碎屑的累計百分含量;Leq為碎屑質(zhì)量為MLeq時所對應的等效邊長;M為計算尺度內(nèi)碎屑的總質(zhì)量。集寧板巖巖爆碎屑粒度-質(zhì)量分形維數(shù)的計算結(jié)果如圖5所示，通過對散點圖的線性擬合得到擬合直線的斜率α=0.27, 進而計算得其分維值為2.73。

圖5 粒度-質(zhì)量的分形維數(shù)Fig.5 Fractal dimension of granularity-mass

3.3 長度/寬度/厚度-數(shù)量的分形

根據(jù)李德建等[14]的研究，將所測量巖爆碎屑的長度、寬度和厚度(均為該方向上的最大值)，按照2.1節(jié)中所述的等效邊長-累計數(shù)量的分形維數(shù)計算方法，進行長度/寬度/厚度-數(shù)量的分形維數(shù)計算,計算結(jié)果如圖6所示(圖中s代表長度/寬度/厚度，單位為mm)，分別對長度、寬度和厚度方向得到的散點圖進行線性擬合得到斜率，即分維值分別為1.70,1.91,1.62，可見寬度、長度、厚度方向的分維值依次減小。

圖6 3種尺度確定分形維數(shù)Fig.6 Fractal dimensions of length/width/thickness vs.quantity

3.4 長寬比/長厚比/寬厚比-數(shù)量的分形

長寬比/長厚比/寬厚比-數(shù)量的分形維數(shù)計算方法與長度-數(shù)量的計算方法相同，JNB-1#碎屑長寬比/長厚比/寬厚比-數(shù)量分形維數(shù)的計算結(jié)果如圖7所示，分別對長寬比/長厚比/寬厚比-數(shù)量得到的散點圖進行線性擬合得到斜率，即分維值分別為2.00，1.64，1.73，從圖7中可看出長寬比和寬厚比的數(shù)據(jù)非常接近，但斜率卻差別較大，可見長寬比的值比較集中，即碎屑的長寬比集中在一定區(qū)間內(nèi)。

圖7 3種尺度比確定分形維數(shù)Fig.7 Fractal dimensions of ratios of length to width, length to thickness,and width to thickness vs. quantity

4 巖爆碎屑分形結(jié)果綜合分析

通過對板巖巖爆碎屑的分形計算和關(guān)系曲線的線性擬合，即可得到擬合直線斜率，進而得到分維值如表2所示(表中n表示累計數(shù)量;m表示累計質(zhì)量;l，w，h分別表示長度、寬度和厚度;等效邊長用Leq表示)。

表2 板巖巖爆碎屑的分形計算結(jié)果

由表2可知，表征巖塊形態(tài)的衡量指標不同，計算得到的分維值則不同，說明指標的選取對分形維數(shù)有一定影響，其中按粒度-質(zhì)量的分形維數(shù)最大，厚度-數(shù)量的最小，從分形維數(shù)的物理意義來看，其值大小與破碎的難度呈正相關(guān)性。對比3種尺度及3種尺度之比的分形維數(shù)可以發(fā)現(xiàn)厚度方向分維值最小，即厚度方向的破碎相對較容易，其次是長度方向，最后是寬度方向，這與板巖的板狀結(jié)構(gòu)有關(guān)，并與破壞之后的破壞形式(見圖2)相對應。其次長、寬、厚3個方向的分維值有明顯差異，且寬厚方向的分維值差值最大為0.29，長厚方向的分維值差值最小為0.08。從表2還可看出，分維值的均值第Ⅰ類為最大，達到2.50，最小為第Ⅱ類，為1.74。從方差來看第Ⅱ類的最小，即計算得到的分維值離散性最小，是3類方法中較好的一種。此外長厚比-數(shù)量的分形維數(shù)值可以表征巖爆碎屑的片狀特征，長寬比-數(shù)量和寬厚比-數(shù)量的分維值分別可以表征巖爆碎屑的條狀和塊狀特征，且長厚比、長寬比和寬厚比的分維值越大，片狀、條狀、塊狀特征越明顯，所以由表2的分形維數(shù)可知板巖巖爆碎屑的條狀、塊狀、片狀特征依次減弱，即碎屑大部分是條狀、條塊狀或塊狀組成，片狀特征不明顯。

5 結(jié) 論

分維值是巖爆破壞程度的一個重要衡量指標，研究板巖的巖爆破壞后的分形特征有一定的應用指導意義，本文根據(jù)對集寧板巖巖爆試驗后的碎屑進行量測，并通過不同的衡量指標計算并對比分維值后，得出以下3條具體結(jié)論，可供類似的工程參考。

(1) 表征巖塊形態(tài)的衡量指標對板巖巖爆碎屑分形維數(shù)的計算有一定程度的影響，分維值分布在1.62～2.73之間，且按長度/寬度/厚度-數(shù)量計算得到的分維值離散性最小，可作為計算板巖巖爆碎屑分維值的優(yōu)選衡量指標。

(2) 板巖巖爆碎屑的分維值在長、寬、厚3個方向上有明顯差異，最小和最大差值分別為0.08和0.29，且厚度、長度、寬度3個方向的破碎難度依次增強。

(3) 長寬比、寬厚比和長厚比與數(shù)量的分維值依次減小，即板巖巖爆碎屑的條狀、塊狀、片狀特征依次減弱。

[1] 齊慶新，歐陽振華，趙善坤，等.我國沖擊地壓礦井類型及防治方法研究[J].煤炭科學技術(shù)，2014，42(10)：1-5.

[2] COOK N G W. The Failure of Rock[J]. International Journal of Rock Mechanics and Mining Science and Geomechanics Abstracts, 1965, 2(4): 389-403.

[3] COOK N G W, PRETORIOUS H E, ORTLEPP J P G,etal. Rock Mechanics Applied to the Study of Rockbursts[J]. Journal of the South African Institute of Mining and Metallurgy, 1966, 66(10): 435-528.

[4] 譚以安. 巖爆形成機理研究[J]. 水文地質(zhì)工程地質(zhì)，1989，16(1):34-38.

[5] 何滿潮，苗金麗，李德建，等. 深部花崗巖試樣巖爆過程實驗研究[J]. 巖石力學與工程學報，2007，26(5)：865-876.

[6] 何滿潮，賈雪娜，苗金麗，等. 巖爆機制及其控制對策實驗研究[C]∥第十一次全國巖石力學與工程學術(shù)大會論文集. 武漢：中國電力出版社，2010：46-56.

[7] 汪 洋，尹健民，李永松，等. 基于巖體開挖卸荷效應的巖爆機理研究[J].長江科學院院報，2014,31(11):120-124.

[8] 仝 躍，陳 亮，黃宏偉. 基于PSO-SVM算法的高放廢物處置北山預選區(qū)巖爆預測[J].長江科學院院報，doi:10.11988/ckyyb.20160058.

[9] 郭慶清，劉磊磊，張紹和，等. 基于組合賦權(quán)法和聚類分析法的巖爆預測[J]. 長江科學院院報，2013,30(12):54-59.

[10]謝和平. 巖爆的分形特征和機制[J]. 巖石力學與工程學報，1993，12(1)：28-37.

[11]謝和平，高 峰，周宏偉，等. 巖石斷裂和破碎的分形研究[J]. 防災減災工程學報，2003，23(4)：1-9.

[12]何滿潮，楊國興，苗金麗，等. 巖爆實驗碎屑分類及其研究方法[J]. 巖石力學與工程學報，2009，28(8)：1521-1529.

[13]單曉云，李占金. 分形理論和巖石破碎的分形研究[J]. 河北理工學院學報，2003，25(2)：11-17.

[14]李德建，賈雪娜，苗金麗，等.花崗巖巖爆試驗碎屑分形特征分析[J]. 巖石力學與工程學報, 2010，29(1)：3280-3289.

(編輯：姜小蘭)

Fractal Characteristics of Fragments Generated from RockburstTest on Slate Taken from Jining

REN Fu-qiang1,2, CHANG Yuan1, WANG Dong3, CHANG Li-yue1, XIA Yuan-yuan1,2

(1.School of Mechanics and Civil Engineering, China University of Mining and Technology，Beijing 100083，China; 2.State Key Laboratory for Geomechanics & Deep Underground Engineering, China University of Mining and Technology， Beijing 100083, China; 3. College of Earth Sciences, Chengdu University of Technology, Chengdu 610059, China)

In order to replenish the study of slate’s rockburst test, we carried out indoor rockburst test on slates from Jining by loading in three directions, gradual unloading from one surface, and then vertical gradual loading tosimulate the rockburst of slate using deep rockburst simulation test system. Then we collected and screened the rockburst fragments, and calculated the fractal dimensions of fragments of particle size greater than 10mm according to different forms of rock. The results indicate that different forms of rock have influence on the fractal dimension which ranges from 1.62-2.73. The fractal dimensions of length, width, and thickness to quantity have the smallest discrete and can be taken as the preferred indicators. The degree of broken difficulty in thickness, length, width is enhanced gradually. In addition, the fractal dimensions of ratios of length to thickness, length to width and width to thickness could characterize strip, sheet, and block of fragments,respectively. The results of the test indicate that slate rockburst fragments are mainly strip, strip-block and block and the sheet feature of fragments is unobvious.

slate; indoor rockburst test; fragments; different forms of rock; fractal

2015-09-06；

2015-10-16

國家自然科學基金重大項目(51134005)

任富強(1992-),男,山西呂梁人,博士研究生，研究方向為巖爆機理及其控制，(電話)18813089558(電子信箱)564123634@qq.com。

10.11988/ckyyb.20150742

2016,33(10):102-105,110

TU45

A

1001-5485(2016)10-0102-04