侯敬民李 巖劉偉東郭英杰于慶波劉少虹

(1.內(nèi)蒙古平莊能源股份有限公司,內(nèi)蒙古自治區(qū)赤峰市,024076; 2.天地科技股份有限公司開采設(shè)計事業(yè)部,北京市朝陽區(qū),100013; 3.煤炭科學研究總院開采研究分院,北京市朝陽區(qū),100013; 4.內(nèi)蒙古平莊能源股份有限公司古山煤礦,內(nèi)蒙古自治區(qū)赤峰市,024076)

基于奇異值分解法的古山礦微震系統(tǒng)監(jiān)測能力優(yōu)化?

侯敬民1李 巖2,3劉偉東4郭英杰1于慶波1劉少虹2,3

(1.內(nèi)蒙古平莊能源股份有限公司,內(nèi)蒙古自治區(qū)赤峰市,024076; 2.天地科技股份有限公司開采設(shè)計事業(yè)部,北京市朝陽區(qū),100013; 3.煤炭科學研究總院開采研究分院,北京市朝陽區(qū),100013; 4.內(nèi)蒙古平莊能源股份有限公司古山煤礦,內(nèi)蒙古自治區(qū)赤峰市,024076)

微震監(jiān)測系統(tǒng)已成為古山煤礦礦山動力災害的主要監(jiān)測手段.為了保證微震系統(tǒng)定位和能量計算的準確性,采用奇異值分解法對古山煤礦微震監(jiān)測系統(tǒng)的現(xiàn)用臺網(wǎng)監(jiān)測能力進行評價,從而發(fā)現(xiàn)了微震系統(tǒng)定位誤差較大區(qū)域,由此制定監(jiān)測系統(tǒng)優(yōu)化方案,并對微震系統(tǒng)的臺網(wǎng)布置進行優(yōu)化.最后采用現(xiàn)場爆破試驗對優(yōu)化結(jié)果進行校驗,結(jié)果表明微震系統(tǒng)的監(jiān)測誤差明顯下降.

采煤工作面 動力災害 微震系統(tǒng) 奇異值分解法 定位誤差 強礦壓顯現(xiàn)臺網(wǎng)布置優(yōu)化

微震監(jiān)測系統(tǒng)的定位精度是對沖擊地壓等動力災害預測預報的基礎(chǔ).微震事件的定位是礦山微震研究的第一步,也是最重要的一步.微震事件定位的精確性取決于多種因素,比如臺站分布、速度模型、震相讀取誤差、走時區(qū)域異常、定位算法、設(shè)備運行狀態(tài)和環(huán)境噪音、定位算法等,其中臺網(wǎng)布置和走時區(qū)域異常是現(xiàn)實中可以人為優(yōu)化改進的.合理的臺網(wǎng)布置可以提高定位精度,增強監(jiān)測效果,降低設(shè)備維護和數(shù)據(jù)處理成本.可以說,臺網(wǎng)的優(yōu)化布置是利用好所有臺站、使網(wǎng)內(nèi)所有臺站有效配合發(fā)揮最大作用的關(guān)鍵.

微震系統(tǒng)臺網(wǎng)定位誤差的計算方法主要包括D值優(yōu)化法、蒙特卡洛法、動態(tài)姿態(tài)網(wǎng)組合法等.本文在上述方法的基礎(chǔ)上,采用奇異值分解法,通過震源參數(shù)的線性相關(guān)性評價震源定位結(jié)果的優(yōu)劣,從而對微震臺網(wǎng)定位精度進行分析.

古山煤礦三井已出現(xiàn)強礦壓顯現(xiàn),而作為礦山動力災害的主要監(jiān)測手段,微震監(jiān)測系統(tǒng)已在古山煤礦取得較好的應用效果,為保證礦震定位和能量計算的準確,提高預測預報礦井沖擊地壓的可行性,采用奇異值分解法對古山煤礦微震監(jiān)測系統(tǒng)的現(xiàn)用臺網(wǎng)的監(jiān)測能力進行評價,從而發(fā)現(xiàn)微震系統(tǒng)定位誤差較大區(qū)域,由此對系統(tǒng)進行優(yōu)化.

1 微震臺網(wǎng)監(jiān)測能力評價原理

假設(shè)微震臺網(wǎng)包含n個臺站,各臺站坐標β＝(xi,yi,zi),i＝1,2,…,n.震源坐標為β＝(x0,y0,z0),微震事件震源的發(fā)震時刻為t0.震源參數(shù)θ＝(x0,y0,z0,t0).當發(fā)生一個微震事件時,該臺網(wǎng)記錄到了m個觀測到時tj,j＝1,2,…,m,第i個臺站理論走時Ti和到時殘差ri.

在震源參數(shù)近似值θ?＝(x?0,y?0,z?0,t?0)處對走時Ti應用一階Taylor展開式,并結(jié)合觀測到時ti可以得到到時殘差ri的表達式.使用矩陣表示如下:

式中:A——m×n維微震波走時偏導數(shù)矩陣;

b——常量;

m——收到信號的拾震器的個數(shù);

x——待求參數(shù),包括震源坐標(x0,y0,z0)和發(fā)震時刻t0.

因為可以進行微震事件震源定位的前提條件是必須最少有4個臺站可以判讀震相初至,因此此處n≥4.

對矩陣A進行奇異值分解,得到:

式中:U——m×m維矩陣,其列向量為AAT的m個正交歸一本征矢量;

V——n×n維矩陣,其列向量為AAT的正交歸一本征矢量;

S——m×n維對角奇異值矩陣,元素包括λ1,λ2,…λn;

AT——A的轉(zhuǎn)置;

T——矩陣轉(zhuǎn)置運算.

再由估計參量空間的協(xié)方差矩陣推出觀測值的協(xié)方差矩陣C(r),主對角元素是各觀測值的到時方差,可得到估計參量空間的協(xié)方差矩陣H.

H的主對角元素是震源參數(shù)標準誤差的估計值.震中和震源深度的標準誤差分別為:

式中:H11、H22、H33——矩陣H的主對角元素;

σxy——震中的標準誤差;

σz——震源深度的標準誤差.

協(xié)方差矩陣H的非對角元素可以用于估計震源參數(shù)間的線性相關(guān)性,而較高的震源參數(shù)線性相關(guān)性會導致震源定位迭代過程失穩(wěn),對計算結(jié)果和收斂速度影響很大.因此,依據(jù)式(4)的計算結(jié)果繪制云圖,對微震臺網(wǎng)的定位誤差進行評估,從而確定微震臺網(wǎng)的監(jiān)測能力及其優(yōu)劣性.

2 古山礦三井現(xiàn)用微震臺網(wǎng)監(jiān)測能力評價

首先對古山礦三井現(xiàn)用微震臺網(wǎng)監(jiān)測能力進行評價,微震臺網(wǎng)評價取P波平均波速為4000 m/s,到時讀取方差為0.05 s,以能量為100 J的微震事件的定位精度作為評價基準.

依據(jù)式(4)繪制現(xiàn)有微震監(jiān)測臺網(wǎng)定位誤差云圖,如圖1所示.三井井下共安裝9臺拾震器.其中,S1、S7拾震器位于東068-1工作面切眼附近,S3、S6拾震器分別位于東068-1回風巷和運輸巷中距離石門較近的位置,其余S12、S4、S5、S8、S2拾震器皆位于三井開拓大巷內(nèi).由圖1可知,現(xiàn)有微震臺網(wǎng)存在兩個問題.

(1)東068-1工作面切眼處存在定位誤差劇烈變化帶這是由于切眼處于微震臺網(wǎng)包圍區(qū)域外側(cè),且該處S1和S7拾震器距離太近,二者相距僅約為110 m;與之對比,S1和S2拾震器的距離約為310 m,其外側(cè)區(qū)域就沒有定位誤差的劇烈變化帶.

(2)西068掘進工作面附近區(qū)域,尤其是S12拾震器左側(cè)區(qū)域定位誤差較大且存在誤差劇烈變化帶.這是由于該區(qū)域拾震器布置較少,微震臺網(wǎng)對其監(jiān)測能力較弱,導致定位誤差出現(xiàn)了明顯的增加.

鑒于上述問題,目前微震臺網(wǎng)將不能滿足古山三井的監(jiān)測需求,由此,將對古山礦三井微震臺網(wǎng)布置進行優(yōu)化.

圖1 現(xiàn)用臺網(wǎng)震中定位誤差云圖

3 微震臺網(wǎng)優(yōu)化

3.1 優(yōu)化方案及其定位效果

根據(jù)目前礦井采掘情況及現(xiàn)用微震臺網(wǎng)存在的問題,提出微震臺網(wǎng)布置優(yōu)化方案.

將S4拾震器移至+320西運輸大巷與回風下山交叉處,移除S5,將S6移至東068-1運輸巷巷道口處.優(yōu)化方案需要移動和去除的拾震器詳細信息見表1,表中拾震器坐標還要根據(jù)實際安裝完成后通過井下實測結(jié)果進一步修正.

表1 優(yōu)化方案中需要操作的拾震器信息

依據(jù)式(4)繪制優(yōu)化方案的微震監(jiān)測臺網(wǎng)定位誤差云圖,如圖2所示.

圖2 優(yōu)化方案震中定位誤差云圖

3.2 優(yōu)化方案與現(xiàn)用臺網(wǎng)定位效果對比

對比圖1和圖2可知,將S4移至+320西運輸大巷與回風下山交叉處后,S12左側(cè)區(qū)域的定位誤差明顯下降,原本S12左側(cè)區(qū)域的定位誤差最低為70 m,優(yōu)化后降低至40 m以下,并且誤差劇變情況大大減弱,誤差劇變帶基本消失,這對西068工作面微震事件的監(jiān)測是有利的.

將拾震器S6移至東068-1運輸巷口后,定位誤差小于30 m的區(qū)域隨之向左側(cè)增加,擴展至拾震器S12上,對西068掘進工作面上下兩巷的監(jiān)測精度進一步提高.

4 現(xiàn)場爆破校驗

在工作面回采過程中,采用煤層卸壓爆破和頂板預裂爆破作為沖擊危險區(qū)域的解危措施,微震監(jiān)測系統(tǒng)記錄到這些爆破事件的波形.統(tǒng)計10個能量大于103J的典型放炮事件,對微震監(jiān)測系統(tǒng)優(yōu)化前后實際監(jiān)測誤差與計算誤差進行比較,微震系統(tǒng)的計算誤差由式(4)計算,如表2所示.首先,對比實際的震中和震源誤差與計算得出的震中和震源誤差,已知微震監(jiān)測系統(tǒng)震中誤差允許范圍為± 20 m,震源誤差允許范圍為±50 m,則震中吻合率為75%,震源吻合率為90%.可知本文對定位精度的評價方法可信度較高,可以為現(xiàn)場提供可靠參考.然后,對比微震監(jiān)測系統(tǒng)臺網(wǎng)優(yōu)化前后監(jiān)測的實際誤差,可知通過微震監(jiān)測系統(tǒng)臺網(wǎng)優(yōu)化監(jiān)測誤差下降了8%～19%.

利用奇異值分解技術(shù)相結(jié)合的方法對微震監(jiān)測臺網(wǎng)進行定位誤差計算,由此對微震臺網(wǎng)布置進行優(yōu)化,對保證微震系統(tǒng)的監(jiān)測精度以及對沖擊地壓災害的預測預報具有重要的指導意義.

表2 系統(tǒng)優(yōu)化前后能量大于103J的放炮事件監(jiān)測與計算誤差對比

5 結(jié)論

(1)微震臺網(wǎng)的布置和拾震器安裝均對微震定位的精確性起著重要作用.首先,微震臺網(wǎng)布置優(yōu)良,將會大幅減少甚至消除定位誤差較大區(qū)域,提高定位精度和速度.其次,拾震器安裝基礎(chǔ)越好,噪音越低,接收的震動波越清晰,將有利于確定P波初至時刻,同時還可增加拾震器監(jiān)測范圍.

(2)以微震事件的定位精度作為評價基準,對古山礦三井現(xiàn)用微震臺網(wǎng)的監(jiān)測能力進行評價,發(fā)現(xiàn)了古山礦三井已用微震臺網(wǎng)存在的問題.

(3)針對古山礦三井現(xiàn)用微震臺網(wǎng)存在的問題,基于微震臺網(wǎng)的優(yōu)化布置原則,對現(xiàn)用臺網(wǎng)進行優(yōu)化.通過優(yōu)化,微震監(jiān)測誤差下降了8%～19%.

(4)通過井下爆破對本文的計算結(jié)果進行校驗,由校驗結(jié)果可知本文的評價方法可信度較高,可以為現(xiàn)場提供可靠參考.

[1] 王元杰,鄧志剛,王傳朋.提高深部開采微震事件定位精度的研究[J].中國煤炭,2011(12)

[2] 印海南.高精度三維地震勘探在山西晉城礦區(qū)的應用[J].中國煤炭,2011(2)

[3] 焦遠碧,吳開統(tǒng)等.我國地震臺網(wǎng)監(jiān)測能力及臺網(wǎng)觀測條件質(zhì)量評定[J].中國地震,1990(4)

[4] 朱元清,趙仲和.提高地震定位精度新方法的研究[J].地震地磁觀測與研究,1997(5)

[5] 姜長寧.江蘇數(shù)字地震臺網(wǎng)地震定位能力分析[J].地震學刊,2001(2)

[6] 郭飆.首都圈數(shù)字地震臺網(wǎng)的微震定位實驗[J].地震地質(zhì),2002(3)

[7] 張玲,梁向軍,董春麗等.山西數(shù)字地震臺網(wǎng)監(jiān)測能力分析[J].山西地震,2010(4)

[8] 鞏思園,竇林名,曹安業(yè)等.煤礦微震監(jiān)測臺網(wǎng)優(yōu)化布設(shè)研究[J].地球物理學報,2010(2)

[9] Sato Y,Skoko D.Opimum distribution of seismic observation pointⅡ[J].Bull of earthquake Res Inst Tokyo Univ,1965(43)

[10] 潘俊鋒,秦子晗,馮美華等.巖漿巖床下伏短壁綜放面集中靜載荷型沖擊啟動原理[J].巖土力學, 2015(9)

[11] 劉少虹.動載沖擊地壓機理分析與防治實踐[D].北京:煤炭科學研究總院,2013

[12] 李宏鑒,斯圖爾特.微震臺網(wǎng)的原理及應用[M].北京:地震出版社,1984

[13] 朱介壽.地震學中的計算方法[M].北京:地震出版社,1988

[14] 田玥,陳曉非.地震定位研究綜述[J].地球物理學進展,2002(1)

[15] 陳祥熊.地震精確定位問題綜述[J].福建地震, 2007(1～2)

[16] 陳法兵,王穎,任文濤等.基于奇異值分解法的微震子臺網(wǎng)監(jiān)測能力分析[J].地震地磁觀測與研究,2015(1)

[17] 楊志國,于潤滄,郭然等.基于微震監(jiān)測技術(shù)的礦山高應力區(qū)采動研究[J].巖土力學與工程學報, 2009(2)

[18] 黃慶國,高潤平.特厚煤層綜放面高精度微震監(jiān)測技術(shù)的應用[J].中國煤炭,2008(8)

Optimization of monitoring ability of microseismic system in Gushan Mine based on singular value decomposition method

Hou Jingmin1,Li Yan2,3,Liu Weidong4,Guo Yingjie1,Yu Qingbo1,Liu Shaohong2,3
(1.Inner Mongolia Pingzhuang Energy Co.,Ltd.,Chifeng,Inner Mongolia 024076,China; 2.Coal Mining and Designing Department,Tiandi Science and Technology Co.,Ltd.,Chaoyang,Beijing 100013,China; 3.Coal Mining Branch,China Coal Research Institute,Chaoyang,Beijing 100013,China; 4.Gushan Mine,Inner Mongolia Pingzhuang Energy Co.,Ltd.,Chifeng,Inner Mongolia 024076,China)

Microseismic monitoring system had been the main monitoring method of mine dynamic disaster in Gushan Coal Mine.To ensure the system localization and energy calculation, singular value decomposition was adopted to evaluate the monitoring ability of current network of microseismic system,which had found out the area with larger location error,so monitoring system optimum proposal for network layout was worked out.Field blasting tests for verifying the optimizing results showed that the monitoring error of microseismic system declined obviously.

coal mining face,dynamic disaster,microseismic system,singular valued composition method,location error,strong strata behaviors,network layout optimization

TD326

A

侯敬民(1974-),男,赤峰平莊人,漢族,高工,從事采礦工程與煤礦安全科研及管理工作。

(責任編輯 張毅玲)

中國煤炭科工集團有限公司科技創(chuàng)新基金項目(KJ-2015-TDKC-05、KJ-2015-TDKC-11、KJ-2016-TDKC-01)