安慶銅礦充填體內破壞區和非破壞區地震波衰減規律研究

袁本勝，史曉鵬，郭利杰，董凱程

（1.北京礦冶研究總院，北京 100160；2.金屬礦山智能開采技術北京市重點實驗室，北京 100160）

安慶銅礦充填體內破壞區和非破壞區地震波衰減規律研究

袁本勝1，2，史曉鵬1，郭利杰1，董凱程1

（1.北京礦冶研究總院，北京 100160；2.金屬礦山智能開采技術北京市重點實驗室，北京 100160）

本文首先采用爆破震動測試法對充填體內的地震波進行測試，然后根據數理統計原理運用最小二乘法進行方程回歸擬合，獲得測試區段充填體內破壞區和非破壞區的地震波衰減公式，進而通過解方程求算充填體內安全允許震速。然后采用聲波測試法根據聲波傳播速度的變化分析爆破引起的充填體破壞狀況驗證所求得的充填體地震波衰減公式和安全允許震速的正確性。

地震波；衰減規律；充填體；爆破

隨著膠結充填等充填新工藝的逐漸完善，采用充填采礦法的礦山比例逐年增加，本文所述礦山即采用了充填采礦法開采。隨著采礦作業的不斷開展，該礦形成了兩面充填體中礦柱回采、三面充填體下礦房回采以及充填體下水平礦柱回采的局面。為了在回采過程中有效保護充填體，首先需要了解充填體的震動破壞特性。本文首先采用爆破震動測試法回歸擬合充填體破壞區和非破壞區的地震波衰減公式，然后通過解方程求得充填體內安全允許震速；然后采用聲波測試法根據聲波傳播速度變化分析爆破引起的充填體破壞狀況以驗證所得充填體內破壞區和非破壞區的地震波衰減規律的正確性。

1 充填體概述

被測礦區地層主要為三迭系中上統灰巖和粉矽巖，出露于礦區的東部與南側，露頭發育程度中等，地層較為完整。礦巖主要成分為大理巖、矽卡巖和閃長巖。礦體傾向西南，傾角較陡。礦區構造基本為單斜層，地層走向為北西－南東，傾向東北，傾角中等至平緩。礦區褶皺主要有西馬鞍山背斜和龜形山背斜。斷層和褶皺帶較發育，對礦床影響明顯，構造發育程度中等。

測試場地位于－400m水平7＃采場充填體中間部位，在－400m至－460m天井聯絡道中－416m位置延伸原來的充填體調查橫巷，橫巷深20m，寬1.5m，測點和爆源點均布置在橫巷頂板。

對1＃礦體7＃、9＃礦房－488m、－471m 和－416m水平充填體做了開挖調查，并對取得的150個試樣進行了測試，其平均強度為0.78MPa，其中在本次試驗點位置所取試樣的平均強度為0.72MPa。

礦山生產過程中，1∶10配比的充填料漿平均濃度為70.38%，28天和60天強度平均值分別為1.14MPa、1.44MPa；1∶10配比濃度為70%～73%的充填體試塊28天抗壓強度試驗室試驗結果為1.23M～1.34MPa。

2 測試系統搭建

爆破震動測試系統由傳感器、震動測試儀和微型計算機等組成。傳感器為磁電式垂直型震動速度傳感器，它為彈簧－質量系統，當物體震動時，由于慣性導致線圈與磁鋼相對運動，切割磁力線人而在線圈內產生感應電壓。震動測試儀采用高速微控制器，將速度傳感器輸出的電壓量進行處理，然后由12位高速A/D轉換器對電壓量進行量化并將量化結果保存到存儲器內。通過該儀器與計算機的RS－232串行接口通信可將測試的數據結果存入計算機磁盤，之后利用配套軟件即可進行分析處理。

聲波速度測試采用數字存儲示波器和計算機采集和存儲波形。根據雙通道波形起始點判讀直接得出兩測點間聲波傳播時間，由此和兩點間距離即可求出傳播速度。當充填體破壞時，波速會有不同程度的降低；反過來，根據同一測試區間的波速變化可分析充填體的破壞狀況。

測點包括震動速度測點和聲波速度測點，均布置在巷道頂板，見圖1。震動速度測點共8個，編號為1＃～8＃。聲波速度測點有2個，編號為a和b。圖中A、B、C、D、E為5個爆破用鉆孔，垂直向上，孔深0.5 m。

圖1 測試點位置示意圖

圖2 第1次測試獲得的波形圖

3 爆破震動衰減規律分析

圖2為第1次測試獲得的波形（限于篇幅，其他波形略）。根據5次測試獲得的波形，經過計算得出比例藥量和峰值震速結果見表1。

比例藥量：ρ＝Q1/3/R，式中 Q 表示炸藥量，齊發爆破時為總裝藥量，微差爆破時為最大一段裝藥量，單位kg；R表示測點與爆破中心的距離，m。

峰值震速：V＝震動測試儀測得的峰值電壓/傳感器靈敏度。

從表1可以看出，在5次測試中，比例藥量取值范圍為0.0232～0.2656kg1/3/m，峰值質點震動速度最小值為0.88cm/s，最大值為14.85cm/s。峰值質點震動速度與比例藥量在雙對數坐標系中的散點圖參見圖3。從散點圖可以看出，當比例藥量由小逐漸增大時，震速的變化趨勢也是由小變大，但當比例藥量較大時，震速增長趨勢變慢，大約在0.1時出現拐點。此拐點可以理解為充填體破壞區和非破壞區的分界點，它的值即是本文需推算的臨界比例藥量及與之對應的臨界震速，而關鍵之所在就是如何確定此拐點。

表1 5次爆破震動測試的比例藥量和峰值震速

具體推算方法為：按比例藥量由小到大的順序逐漸增加回歸點并進行回歸，同時作相關系數和F檢驗。當檢驗值出現拐點時，即以此前檢驗的最佳值對應的回歸公式作為正常狀態下（也即非破壞狀態下）充填體的地震波衰減公式。再按比例藥量由大到小的順序逐漸增加回歸點并進行回歸，同時作相關系數和F檢驗以得出非正常狀態下（也即破壞狀態下）充填體的衰減公式。兩條回歸線的交點即是所求的臨界比例藥量和臨界震速。

3.1 非破壞狀態下充填體中地震波衰減規律

從比例藥量小于0.07的所有點開始回歸，然后按比例藥量由小到大的順序逐點加入進行回歸（因篇幅所限，回歸分析數據表略）。

回歸曲線相關系數R范圍為0.722～0.907，若在0.05水平檢驗，查相關系數檢驗表得出其范圍為0.304～0.381，實際回歸得到的R均超出該范圍，表明地震波在充填體內的衰減規律按公式V＝K（Q1/3/R）α擬合的效果較好。對 F 值在 0.05 和0.01水平檢驗，查F分布表可得其范圍分別為4.08～4.24和7.31～7.77，而上表中的 F 值在25.048～147.725之間，全部超出查表值，可見擬合的曲線關系高度顯著。比較各次回歸的剩余標準差s，以s值最小為優，也就是說以s為最小值時回歸得出的K、α值作為未被破壞的充填體中地震波衰減公式的K、α值。從表1中數據可求出，未被破壞的充填體中地震波衰減公式見式（1）。

3.2 破壞狀態下充填體中地震波衰減規律

對于充填體破壞區的地震波衰減公式回歸，數據選取按比例藥量從大到小的順序取至其稍小于0.1（未被破壞的充填體中地震波衰減公式回歸取點的最大比例藥量為0.1027）即可。

根據前述非破壞區充填體同樣的計算方法得到破壞區充填體中地震波衰減公式（式（2））。

最終獲得的破壞區和非破壞區充填體中地震波衰減公式的回歸線見圖3。

將破壞區和非破壞區充填體中地震波衰減公式（1）、式（2）聯立解方程得到安全允許震速 V＝9.22cm/s，對應的比例藥量ρ＝0.0955kg1/3/m。

圖3 破壞區與非破壞區充填體中震動數據散點與回歸線

4 利用聲波測試法驗證衰減規律

聲波速度測試點為2個固定接收點，通過錘擊人為產生聲源。錘擊點與接收點同在一條直線上。兩接收點的距離為0.71m。在5次測試爆破之前和前4次爆破之后都進行了聲波測試，雙通道數據采集器記錄的波形分別見圖4。取得的聲波傳播時間和經計算得出相應的聲波傳播速度見表2。

表2 充填體聲波測試數據

由表2可知，從第1次爆破后波速就開始減小，以后各次比前一次也都減小。由此可見，第1次爆破就已經對兩測試點的充填體產生破壞，且隨著爆源的逐漸靠近，破壞程度越來越嚴重。

圖4 聲波測試波形

按式（1）計算，在第1次爆破時，藥量為0.075kg，爆心與第a和b兩測試點的距離分別為3.8 6 m和4.6 m，求得此二點的峰值震速分別為11.27cm/s和8.67cm/s。兩測試點所求得的峰值震速在前文所求得的安全允許震速9.22cm/s的兩邊，說明充填體被破壞的位置在兩測試點之間。

上述由爆破震動和聲波速度測試分析兩種方法所得結論一致，相互驗證了正確性，充分說明本研究對具體充填體條件進行實測得出的衰減公式和安全允許震速是可信的。

5 結論

本文對同一次爆破采用爆破震動測試數值擬合法和聲波測試法兩種方法對爆破震動在充填體內的衰減規律進行測試，并根據數理統計原理運用最小二乘法進行方程回歸擬合，取得測試區段充填體內破壞區和非破壞區的地震波衰減規律，并輔以聲波速度測試反映充填體破壞狀況，驗證了前述結果的正確性。獲得了如下結論。

1）按比例藥量逐點加入回歸以剩余標準差s判優，并利用相關系數的顯著性檢驗和F檢驗在未知破壞區和非破壞區的情況下找出了此二區兩條回歸線的回歸系數，據此劃定兩區界線，從而確定充填體內安全允許震動速度。

2）對于本測試的充填體，灰砂比為1∶10時，充填體非破壞區的震動衰減公式為V＝307.912（Q1/3/R）1.494。

3）對于本測試的充填體，灰砂比為1∶10時，充填體破壞區的震動衰減公式為 V＝26.972（Q1/3/R）0.457。其衰減指數明顯低于非破壞區的衰減指數，即爆破震動在破壞區的衰減慢。場地系數也因充填體被破壞而明顯減小。

4）對于本測試的充填體，灰砂比為1∶10時，安全允許震動速度為9.22cm/s。

5）可通過測量充填體內兩個點的聲波速速來判斷充填體的破壞狀況。

［1］趙剛，張銀平.DSVM系列測振儀在爆破振動測試中的應用［J］.工程爆破，1995，1（4）：61－65.

［2］胡建華，雷濤，周科平，等.充填采礦環墇境下爆破振動效應預測與控制［J］.2009，26（4）：81－84.

［3］呂濤，石永強，黃誠，等.非線性回歸法求解爆破振動速度衰減公式參數［J］.巖土力學，2007，28（9）：1871－1878.

［4］李洪濤，舒大強.爆破震動衰減規律的影響因素［J］.武漢大學學報：工學版，2005，38（1）：79－82.

［5］LI Hong－tao，SHU Da－qiang.Influential factors analysis of blasting vibration attenuation law［J］.Engineering Journal of Wuhan University，2005，38（1）：79－82.

［6］ZANGXiu－ping，WANG Hong.Preliminary exploration of the velocity attenuation equation form of blast vibration ［J］.Metal Mine，2002（4）：36－37.

［7］臧秀平，劉升寬，董濤.露天爆破振動參數衰減模型研究［J］.有色金屬：礦山部分，2005，57（5）：31－33.

［8］張繼春，彭瓊芳.巖體爆破地震波衰減規律的現場試驗與分析［J］.遼寧工程技術大學學報：自然科學版，2001，20（4）：399－401.

［9］葉洲元，周志華.爆破震動安全距離的優化計算［J］.中國安全科學學報，2005，15（3）：57－60.

Study on blasting vibration attenuation law for the failure zone and non－failure zone of backfilling in Anqing copper mine

YUAN Ben－sheng1，2，SHI Xiao－peng1，GUO Li－jie1，DONG Kai－cheng1
（1.Beijing General Research Institute of Mining and Metallurgy，Beijing 100160，China；2.Beijing Key Laboratory of Nonferrous Intelligent Mining Technology，Beijing 100160，China）

At first this paper use the method of the blasting vibration to test the seismic wave on the filling body，then apply the principle of mathematical statistics and the least squares to fit of regression equation，and obtain the attenuation formula of the seismic wave at the failure zone and non－failure zone of the filling body，also by solving the equation to calculate the allowable safety vibration velocity of filling body.At last it use the sound wave test method with its features of that the change of the seisimic wave velocity can indicate the destroys state of the filling body，to validate the attenuation formula of the seismic wave at the failure zone and non－failure zone of the filling body.

seismic wave；attenuation law；backfilling；blasting

袁本勝（1982－），男，漢族，工程師，碩士，主要從事礦山安全監測方向的工作。E－mail：yuanbsh＠163.com。

P631

A

1004－4051（2014）S2－0309－04

2014－10－17

“十二五”國家科技支撐計劃項目資助（編號：2013BAB－02B02；2012BAB01B04）；國際合作項目資助 （編號：2011DFA－71990）；北京礦冶研究總院院基金課題資助（編號：YJ－2012－01）；礦山主被動結合高精度微震監測技術與裝備研發項目資助（編號：2012BAK09B03）