吳建星，劉 佳

（武漢科技大學(xué)資源與環(huán)境工程學(xué)院，湖北 武漢，430081）

在采礦工程實(shí)踐中，由于大量人為采掘，以致造成巖土體變形、地面沉降、井巷塌陷以及泥石流等地質(zhì)現(xiàn)象而引發(fā)微震。微震不僅破壞采礦工程設(shè)備和礦區(qū)資源環(huán)境，危及人類(lèi)生命和財(cái)產(chǎn)安全，而且還嚴(yán)重地影響采礦生產(chǎn)。微震以震動(dòng)波的形式突然釋放，其量度為1～3級(jí)。這一類(lèi)地震人們無(wú)法感覺(jué)，僅靠?jī)x器測(cè)出。在礦山微震實(shí)時(shí)監(jiān)測(cè)和震源定位研究方面，邴紹丹等［1］采用小波分析濾波過(guò)的信號(hào)來(lái)確定礦震震源位置，發(fā)現(xiàn)其系統(tǒng)的監(jiān)測(cè)結(jié)果與實(shí)際礦震相吻合；趙龍［2］采用模擬微震試驗(yàn)方法采集礦震波信號(hào)，得到微震理論位置、發(fā)生時(shí)間及其系統(tǒng)定位、定時(shí)的誤差；朱超等［3］研發(fā)了一套用于礦山微震實(shí)時(shí)監(jiān)測(cè)系統(tǒng)，并介紹了該系統(tǒng)的功能和用途。為了預(yù)防和控制礦山微震事件的發(fā)生，本文在文獻(xiàn)［3］的基礎(chǔ)上研制一套礦山微震多通道監(jiān)測(cè)定位系統(tǒng)，并根據(jù)牛頓迭代目標(biāo)定位原理，應(yīng)用該系統(tǒng)對(duì)礦山微震定位進(jìn)行研究，以期為預(yù)報(bào)和控制礦山地質(zhì)災(zāi)害的發(fā)生提供依據(jù)。

1 監(jiān)測(cè)系統(tǒng)的描述

根據(jù)某鐵礦的實(shí)際地質(zhì)情況，本文研制的礦山微震監(jiān)測(cè)定位系統(tǒng)由傳感器、轉(zhuǎn)換器、數(shù)據(jù)傳輸系統(tǒng)、數(shù)據(jù)處理中心等部分組成，如圖1所示。該系統(tǒng)能實(shí)現(xiàn)實(shí)時(shí)在線監(jiān)測(cè)，可作出三維圖形，且成本低廉。

1.1 硬件組成

圖1 監(jiān)測(cè)系統(tǒng)結(jié)構(gòu)圖Fig.1 Monitoring system diagram

微震監(jiān)測(cè)系統(tǒng)硬件由防爆計(jì)算機(jī)、數(shù)據(jù)采集卡、采集卡電源、前置放大器、傳感器、光纖收發(fā)機(jī)和光纜等部分所組成。微震傳感器共設(shè)置了16個(gè)采集通道，所采集的微震信號(hào)由前置放大器進(jìn)行前置放大，再通過(guò)電纜傳輸?shù)綌?shù)據(jù)采集卡，經(jīng)過(guò)光纜后傳至防爆計(jì)算機(jī)對(duì)微震信號(hào)進(jìn)行處理分析，從而實(shí)現(xiàn)對(duì)礦山微震的監(jiān)測(cè)和定位。

1.2 軟件組成

本研究微震監(jiān)測(cè)定位系統(tǒng)采用LabVIEW軟件，該軟件由NI有限公司開(kāi)發(fā)。它使用的是圖形化編輯語(yǔ)言G編寫(xiě)程序，其程序?yàn)榭驁D形式。在微震監(jiān)測(cè)定位過(guò)程中，由于LabVIEW軟件可充分發(fā)揮計(jì)算機(jī)強(qiáng)大的數(shù)據(jù)處理功能而進(jìn)行數(shù)據(jù)采集、數(shù)據(jù)分析、數(shù)據(jù)顯示以及數(shù)據(jù)存儲(chǔ)等。LabVIEW軟件程序框圖使編程更加直觀，處理數(shù)據(jù)更加便捷，并在一個(gè)硬件的情況下，通過(guò)改變其軟件可實(shí)現(xiàn)不同儀器儀表的功能，而且Lab－VIEW的擴(kuò)展性使它能夠?qū)崿F(xiàn)與其他軟件如Matlab的聯(lián)合使用。

1.3 定位原理［4］

如果在礦區(qū)某一區(qū)域發(fā)生微震現(xiàn)象，則只象征著有潛在的巖體破壞危險(xiǎn)，并存在著較大的隨機(jī)性和復(fù)雜性。設(shè)微震源的空間坐標(biāo)為x、y、z，微震事件的起始時(shí)間為t，P波在巖體中傳播的平均速度為c，則第i個(gè)傳感器與微震源之間的非線性到時(shí)方程為

式中：xi、yi、zi分別為第i個(gè)傳感器的空間坐標(biāo)；ti為微震源到第i個(gè)傳感器的到時(shí)，ms；n為接受信號(hào)的傳感器個(gè)數(shù)；x、y、z、t分別為所要求的微震源時(shí)空參數(shù)。

取i個(gè)傳感器中的4個(gè)，與微震源組成的方程組為

圖2為微震定位示意圖。將已知傳感器坐標(biāo)和到時(shí)信息代入式（2）進(jìn)行迭代計(jì)算，計(jì)算結(jié)果為微震源的一次定位坐標(biāo)。經(jīng)過(guò)多次定位計(jì)算后取其平均值，就可得到更為精準(zhǔn)的微震源坐標(biāo)。

圖2 微震定位示意圖Fig.2 Microseismic positioning schematic diagram

2 實(shí)例應(yīng)用

某鐵礦礦體賦存條件復(fù)雜，礦巖受構(gòu)造與蝕變的影響，其巖性差異較大。礦區(qū)內(nèi)有許多巖種具有強(qiáng)烈的水理特性，如矽卡巖、角頁(yè)巖、泥質(zhì)角巖以及粉礦，當(dāng)?shù)V巖經(jīng)過(guò)水的作用后，會(huì)發(fā)生崩解、膨脹、軟化等變化，使其穩(wěn)定性大大降低。為此，應(yīng)用一套多通道微震監(jiān)測(cè)定位系統(tǒng)，通過(guò)合理布置傳感器以及對(duì)所采集的信號(hào)波形進(jìn)行震源定位計(jì)算，就可達(dá)到精準(zhǔn)定位微震源的目的。

2.1 傳感器布置

為了充分發(fā)揮監(jiān)測(cè)系統(tǒng)的作用，必須合理地對(duì)傳感器進(jìn)行空間布置和分配，以滿(mǎn)足在微震監(jiān)測(cè)工程技術(shù)指標(biāo)的條件下，使傳感器陣列監(jiān)測(cè)范圍擴(kuò)大，以充分發(fā)揮監(jiān)測(cè)系統(tǒng)的作用。本研究所用傳感器為PS－60B型傳感器，其監(jiān)測(cè)靈敏度達(dá)到85%以上，符合實(shí)際要求。按照震源定位誤差不大于10m的指標(biāo)要求，建立被測(cè)區(qū)域的物理模型，利用被測(cè)區(qū)域內(nèi)聲傳波速度的相關(guān)信息，對(duì)16個(gè)傳感器陣列內(nèi)的監(jiān)測(cè)范圍進(jìn)行模擬分析，并在模擬過(guò)程中不斷調(diào)整傳感器的位置，使監(jiān)測(cè)范圍內(nèi)的微震定位精度滿(mǎn)足技術(shù)要求。圖3為微震傳感器布置圖。根據(jù)地質(zhì)結(jié)構(gòu)特點(diǎn)，采用分布法來(lái)布置傳感器，使采集到的數(shù)據(jù)精度更高。

圖3 微震傳感器布置圖Fig.3 Microseismic sensor layout drawing

在礦區(qū)內(nèi)，將采集系統(tǒng)放入水平－340m的值班室，24h連續(xù)不斷實(shí)時(shí)在線采集微震數(shù)據(jù)，并將其數(shù)據(jù)分時(shí)段存儲(chǔ)（隔4h存儲(chǔ)1次），以方便數(shù)據(jù)的傳遞與閱讀。16個(gè)不同位置的傳感器將微震信號(hào)采集后，傳遞到采集卡轉(zhuǎn)換為數(shù)字信號(hào)，再通過(guò)光纖轉(zhuǎn)換器轉(zhuǎn)換為光纖信號(hào)，使信號(hào)傳遞時(shí)所受干擾和衰減均較小，以方便遠(yuǎn)程傳輸。采集系統(tǒng)接收信號(hào)時(shí)，利用光線轉(zhuǎn)換器將光纖信號(hào)還原成數(shù)字信號(hào)，就能被工控機(jī)識(shí)別，并通過(guò)相關(guān)軟件顯示出來(lái)。

2.2 震源定位計(jì)算

通過(guò)傳感器的合理布置，現(xiàn)場(chǎng)能采集到大量的實(shí)時(shí)在線微震信號(hào)，應(yīng)用相關(guān)軟件對(duì)這些微震信號(hào)進(jìn)行識(shí)別分析。表1為每個(gè)傳感器的坐標(biāo)和到時(shí)信息等參數(shù)，圖4為計(jì)算震源位置的程序框圖。

表1 傳感器坐標(biāo)和到時(shí)參數(shù)Table1 Sensor coordinate and arrival time

在監(jiān)測(cè)范圍內(nèi)的不同位置布置了16個(gè)傳感器，若選擇每4個(gè)組成一個(gè)方程組計(jì)算震源位置，則有255種不同的組合方法。根據(jù)到時(shí)的長(zhǎng)短，先把到時(shí)較短的組合在一起進(jìn)行定位計(jì)算，再把到時(shí)較長(zhǎng)的組合在一起進(jìn)行定位計(jì)算。對(duì)全部定位結(jié)果取其空間幾何中心值，計(jì)算所有定位結(jié)果到中心點(diǎn)的距離。計(jì)算出其平均距離后，再計(jì)算各個(gè)距離到平均距離的方差。取均方差或取均方差的一個(gè)倍數(shù)，舍去方差大于這個(gè)系數(shù)的定位結(jié)果。通過(guò)取舍，留有10組較精準(zhǔn)的定位結(jié)果如表2所示。把所有計(jì)算結(jié)果相加除10得到平均值大小，就是經(jīng)計(jì)算得到的震源位置。實(shí)踐應(yīng)用表明，震源位置產(chǎn)生的誤差小于8.8m，所得的定位結(jié)果較為精確。

圖4 震源定位程序框圖Fig.4 Positioning block diagram

表2 定位誤差和震源位置Table2 The positioning error and source location

由表2可看出，通過(guò)計(jì)算得到震源位置x為5044、y為7670、z為－342m，與實(shí)際震源位置x為5048、y為7663、z為－340m相比，其計(jì)算結(jié)果x、y、z的誤差分別為－4、7、－2m。

由于巖石密度不同，微震波在不同方位的傳播速度也不同，再加上傳感器不靈敏，導(dǎo)致計(jì)算精度誤差約為10m，但這個(gè)誤差在微震監(jiān)測(cè)系統(tǒng)所允許的范圍內(nèi)，表明本監(jiān)測(cè)系統(tǒng)仍具有較高的實(shí)時(shí)監(jiān)測(cè)能力和定位精度。

3 結(jié)論

（1）引入牛頓迭代法對(duì)礦山微震源進(jìn)行定位計(jì)算，在迭代過(guò)程中只要迭代幾次就可得到精確的解。

（2）對(duì)16個(gè)傳感器進(jìn)行選擇性組合，按照到時(shí)的長(zhǎng)短，將到時(shí)較短的組合在一起，到時(shí)較長(zhǎng)的另組合在一起。通過(guò)進(jìn)一步計(jì)算、取舍和結(jié)合得到震源位置，其產(chǎn)生的誤差在微震監(jiān)測(cè)系統(tǒng)所允許范圍內(nèi)。

（3）本研究研制一套實(shí)時(shí)在線監(jiān)測(cè)定位系統(tǒng)，以爆破地點(diǎn)為微震源，監(jiān)測(cè)微震傳播過(guò)程，并采集、分析微震信號(hào)，確定了精準(zhǔn)的微震源位置，為避免礦山地質(zhì)災(zāi)害的發(fā)生尋求到一條有效途徑。

［1］邴紹丹，潘一山.礦山微震定位方法及應(yīng)用研究［J］.煤礦開(kāi)采，2007，12（5）：1－4.

［2］趙龍.煤礦微震監(jiān)測(cè)定位系統(tǒng)誤差分析［J］.實(shí)驗(yàn)科學(xué)與技術(shù)，2010，8（5）：33－34.

［3］朱超 吳建星，趙智，等.基于LabVIEW的微震實(shí)時(shí)監(jiān)測(cè)系統(tǒng)［J］.工業(yè)安全與環(huán)保，2011，37（10）：47－48.

［4］吳治濤，駱循，李仕雄.聯(lián)合小波變換與偏振分析自動(dòng)拾取微地震P波到時(shí)［J］.地球物理學(xué)進(jìn)展，2012，27（1）：131－136.