基于最小平方法的微震震級與能量關系研究

王元杰，霍永金，張宗文，王傳朋，趙成利

(1.天地科技股份有限公司開采設計事業部，北京100013;2.煤炭科學研究總院開采設計研究分院，北京100013;3.中天合創能源有限責任公司，內蒙古鄂爾多斯017000;4.新汶礦業集團有限責任公司華豐煤礦，山東泰安271413)

微震監測技術 (MS)是通過監測煤巖體受力變形和破壞后發射出的震動波來進行監測工程巖體穩定性的技術方法［1］。微震監測技術作為一種區域性、實時性的監測技術，與其他傳統的監測手段如鉆屑法等相比，具有遠距離、三維、實時、動態監測的特點，并且微震監測具有很高的定位精度，已經成為目前監測煤巖體動力災害的主要技術手段［2］。目前，國內應用的微震監測系統多為引進國外或自發研制，如:新汶華豐煤礦與天地科技股份有限公司合作引進的波蘭ARAMIS M/E微震監測系統;冬瓜山煤礦與唐禮忠等合作引進的南非ISS微震監測系統;凡口鉛鋅礦與李庶林等合作研制的微震監測系統;徐建新和劉英利合作研制的微震監測系統等［3］。

1 微震震級與能量的關系

微震事件對煤巖體穩定性的影響程度大小取決于微震事件能量釋放多少，微震監測系統監測得到的只是各傳感器接收到的由震動機械能轉化成的電壓值，再把監測到的電壓值轉化為能量值，根據此能量值，在考慮地震波在傳播過程中因為反射、折射、透射等造成的衰減基礎上，得到震源能量和距離之間的能量關系式，從而推算出震源能量值［4－5］。

ARAMIS M/E微震監測系統采用綜合能量計算法計算微震事件的能量，能量的計算公式為:

其中，

式中，Ei為第i個拾震器計算得到的能量值;ODi為第i個拾震器距震源的距離;Vi為由理論公式得到的波速;f為極限頻率;Q為介質阻尼;CzTi為第i個拾震器的線性靈敏度;leb為震源深度 (Z坐標+測量處的標高);Ai，k為第i個拾震器的信號振幅;Npi為第i個拾震器初始信號振幅;Nki為第i個拾震器結束信號振幅;Zv為按速率表示的測量范圍;Zb為按bits表示的測量范圍;Takt為每秒取樣次數;G為介質密度。

該方法解決了對微震事件能量估計的問題，計算結果以焦耳為單位。但是，如果信號過強，某個拾振器信號里達到最大振幅的信號過多，通過這種信號得到的能量值準確性將受影響，用該拾振器計算得到的能量值不用于參與計算能量的均值。

微震震級的計算方法源于地震震級的計算方法，地震震級是通過測量地震波中的某個震相的振幅和周期，并考慮地震波傳播過程中的衰減計算得到的一個衡量地震相對大小的量，最早由美國的CharlesF.Riehter和日本的 Kiyoo Wadati在20世紀30年代提出的［6］，根據人對地震的感覺和地震的破壞程度將里氏震級劃分9個級別，3級以下的地震屬于弱震，人們一般不易覺察;6級以上屬于強震，根據震源的不同深度會造成嚴重的破壞和損壞［7］。地震的能量是根據地震儀記錄到的地震波振幅和周期計算得到的，地震的震級由能量推算得出，地震釋放出的能量與震級成正比，能量越大，震級越大［8］。震級 (M)和震源能量 (E)的關系為:

式中，a，b為常數，不同煤礦震級－能量關系式的a，b參數是不一樣的，因而結合煤礦具體的地質條件、開采狀況，很有必要確定不同煤礦的震級－能量關系式參數。

2 最小平方法優化參數

在科學研究中為了找出某些相關量之間的規律，通常采用數據擬合的方法，主要有傳統的三角函數逼近法、Pade逼近法、插值法、最佳一致逼近多項式法、最佳平方逼近、最小平方擬合等，以及現代的模糊逼近法、神經網絡逼近法、小波理論法、支持向量機函數逼近法等［9］。

最小平方擬合是一種最基本的計算方法，廣泛應用于建模中，解決問題時簡明而清晰，特別在大數據分析研究中，具有非常重要的作用和地位。最小平方法通過計算最小化誤差的平方和來求解數據的最佳函數。

假設直角坐標系中的一組數據 (x1，y1)，(x2，y2)，…，(xi，yi)，采用數學優化技術可以找到這些數據之間的函數關系。如果坐標系中的這些數據圖像近似于一條直線，就可以用直線方程:y=a+bx來表達這組數據的函數關系。最小平方法就是找到各種觀測數據之間函數關系的一種數學優化基本計算方法，它的原理是當最小誤差平方和S=Σ [yi－(a+bxi)]2為最小時，即

解得系數a，b為

3 應用驗證

華豐煤礦于2006年引進并裝備了波蘭ARAMIS M/E微震監測系統，截至目前已獲得了大量的監測數據，作為后續分析的基礎，如何保證監測數據的可靠性至關重要［10］。最初ARAMIS M/E微震監測系統計算得到震級對比當地地震臺提供的里氏震級誤差很大，微震監測系統計算得到的震級普遍高于地震臺提供的里氏震級，兩者差值最大的達到1.1級。針對微震系統計算得到震級誤差較大的問題，需要優化震級－能量關系式參數，使得系統計算得到的震級和當地地震臺提供數據更加吻合，提高系統數據的可靠性。華豐煤礦微震監測系統最初震級－能量關系式參數系數a和b，由波蘭地球物理學家Marek Dworak結合其在波蘭多年工作經驗設定，分別為a=2.2，b=1.9。

根據震級和震源發出的總能量之間關系lgE=a+bM，震源發出的總能量E的常用對數與震級M成線性關系，采用多個地震臺提供的震級數據作為計算震級－能量對應關系的樣本點，通過最小化誤差的平方和尋找數據的最佳直線匹配。如表1所示，設地震臺震級為xi，對應的微震系統計算能量的常用對數為yi，運用最小平方法生成震級－能量對數對應直線lgE=3.484+2.123M，如圖1所示，圖中的點代表樣本點，直線為震級－能量對數對應直線。

為驗證新的直線關系對于系統計算震級的誤差，選取多個地震臺數據與系統數據，并通過新的震級－能量關系式參數a=3.484，b=2.123計算得到新的震級，如表2所示。結果顯示采用最小平方法減小了系統震級的誤差，提高了微震監測系統監測數據的可靠性。

表1 震級調整計算對比結果

圖1 震級－能量對數對應直線

表2 震級調整計算對比結果

4 總結

采用最小平方法得出微震監測系統的震級－能量關系式參數:a=3.484，b=2.123。通過該參數得出的系統震級與當地地震臺震級的誤差明顯地減小。震級－能量關系式參數的優化取得了良好的效果，提高了華豐煤礦微震監測數據的可靠性，并為類似礦井微震監測系統的建設提供了借鑒和指導意義。

［1］齊慶新，竇林名.沖擊地壓理論與技術［M］.徐州:中國礦業大學出版社，2008.

［2］竇林名，何學秋.沖擊地壓理論與防治技術［M］.徐州:中國礦業大學出版社，2001.

［3］唐禮忠，潘長良，楊承祥，等.冬瓜山銅礦微震監測系統及其應用研究 ［J］.金屬礦山，2006(10):41－44，86.

［4］張宗文，王元杰，趙成利，等.微震和地音綜合監測在沖擊地壓防治中的應用［J］.煤炭科學技術，2011，39(1):44－47.

［5］王元杰，鄧志剛，王傳朋.提高深部開采微震事件定位精度的研究 ［J］.中國煤炭，2011，37(12):60－63，81.

［6］Aleksander Cianciara，Bogdan Cianciara.Method of predicting tremors on the basis of seismic emission registered in exploitation workings［C］.Tectonophysics，456(2008).

［7］劉英利，徐建新.強震與微震觀測系統的研究［J］.地震工程與工程振動，2006，26(3):264－266.

［8］李庶林，尹賢剛，鄭文達，等.凡口鉛鋅礦多通道微震監測系統及其應用研究 ［J］.巖石力學與工程學報，2005，24(12):2048－2053.

［9］魯振華，張連城.門頭溝礦微震的近場監測效能評估［J］.地震，1989(10):32－39.

［10］陳炳瑞，馮夏庭，李庶林.基于粒子群算法的巖體微震源分層定位方法 ［J］.巖石力學與工程學報，2009，28(4):740－749.