基于密集臺(tái)陣的遼寧老虎臺(tái)礦區(qū)微震事件重定位

摘要：

遼寧撫順老虎臺(tái)煤礦開(kāi)采已百年有余，是礦震災(zāi)害較為嚴(yán)重的礦區(qū)之一。然而，該地區(qū)地震監(jiān)測(cè)臺(tái)站較少，導(dǎo)致對(duì)該礦區(qū)的微震監(jiān)測(cè)能力有所不足。為了更好地監(jiān)測(cè)老虎臺(tái)礦區(qū)微震的活動(dòng)情況，本研究基于2019年11月3—25日在礦區(qū)布設(shè)的兩條線(xiàn)性密集臺(tái)陣，開(kāi)展礦區(qū)微震事件的雙差定位研究。重定位結(jié)果顯示：本次微震活動(dòng)的最大事件為發(fā)生于渾河斷裂附近的ML 2.4級(jí)地震，震源深度為0.5 km。在渾河ML 2.4級(jí)地震發(fā)生之前，ML gt;1.0的微震事件逐漸增多且有向渾河斷裂西端遷移的特征；在渾河ML 2.4級(jí)地震發(fā)生之后，微震活動(dòng)主要集中于礦區(qū)內(nèi)部的F1和F2斷層附近，震源深度主要集中在0.2～1.5 km之間?；谖⒄鹗录臅r(shí)空分布特征，推測(cè)老虎臺(tái)礦區(qū)內(nèi)部的微震活動(dòng)可能受到渾河ML 2.4級(jí)地震的影響，礦區(qū)內(nèi)部可能發(fā)生斷層活化現(xiàn)象。此外，礦區(qū)內(nèi)部的F2斷層可能是一條延伸至渾河斷裂的NNW向隱伏斷層。

關(guān)鍵詞：老虎臺(tái)礦區(qū)；密集臺(tái)陣；微震監(jiān)測(cè)；雙差定位；渾河斷裂；礦震災(zāi)害

doi：10.13278/j.cnki.jjuese.20230102

中圖分類(lèi)號(hào)：P315

文獻(xiàn)標(biāo)志碼：A

孔慶翰，呂子強(qiáng)，張廣偉. 基于密集臺(tái)陣的遼寧老虎臺(tái)礦區(qū)微震事件重定位. 吉林大學(xué)學(xué)報(bào)（地球科學(xué)版），2024，54（5）：16851695. doi：10.13278/j.cnki.jjuese.20230102.

Kong Qinghan， Lü Ziqiang， Zhang Guangwei. Microseismic Events Relocation of" Laohutai Coal Mine in Liaoning Province Based on Dense Seismic Array. Journal of Jilin University （Earth Science Edition）， 2024， 54 （5）： 16851695. doi：10.13278/j.cnki.jjuese.20230102.

收稿日期：20230419

作者簡(jiǎn)介：孔慶翰（1999—），男，碩士研究生，主要從事微震定位方面的研究，E-mail： kqh_lntu@163.com

通信作者：呂子強(qiáng)（1982—），男，副教授，主要從事深部結(jié)構(gòu)成像方面的研究，E-mail： ziqianglyu@sina.com

基金項(xiàng)目：國(guó)家自然科學(xué)基金項(xiàng)目（42274129）

Supported by the National Natural Science Foundation of China （42274129）

Microseismic Events Relocation of" Laohutai Coal Mine in Liaoning Province Based on Dense Seismic Array

Kong Qinghan1， Lü Ziqiang1， Zhang Guangwei2

1. College of Mining， Liaoning Technical University， Fuxin 123000， Liaoning， China

2." National Institute of Natural Hazards， Ministry of Emergency Management of China， Beijing 100085， China

Abstract：

Laohutai coal mine， located in Fushun， Liaoning Province， has been exploited for almost a century and is one of the mining areas with serious mine earthquake disasters. However， due to the lack of dense seismic monitoring stations， the mine’s microseismic event monitoring capacity is insufficient. In order to" monitor the microseismic activity better， this study carried out" double-difference location of microseismic events in the mining area based on two linear dense seismic arrays deployed in "Laohutai coal mine from November 3 to 25， 2019. The relocation results indicate that the largest event of this microseismic activity was the ML 2.4 earthquake that occurred near" Hunhe fault， with an epicenter depth of 0.5 km. Before the occurrence of the Hunhe ML 2.4 earthquake， the microseismic events of ML gt;1.0 gradually increased and migrated to the western end of" Hunhe fault. After the occurrence of the Hunhe ML 2.4 earthquake， the microseismic activities in the mining region were mainly concentrated near" F1 and F2 faults in the mining area， and the depth of the epicenter was mainly concentrated between 0.21.5 km. The spatial and temporal distribution characteristics of the microseismic events indicate that the microseismic activity in" Laohutai coal mine may be affected by the Hunhe ML 2.4 earthquake， and fault activation may occur in the mining area. Additionally， we speculate that the F2 fault may be a NNW-trending buried fault extending to" Hunhe fault based on the observed distribution of the microseismic relocation.

Key words：

Laohutai coal mine; dense seismic array; microseismic monitoring; double-difference location; Hunhe fault; mine earthquake disaster

0" 引言

老虎臺(tái)礦區(qū)位于遼寧省撫順市南部撫順煤田的中部，西臨撫順西露天礦，東接龍鳳礦，北至渾河斷裂，東西長(zhǎng)約4.95 km，南北寬約2 km，礦區(qū)面積約為10 km2。老虎臺(tái)煤礦屬新華夏構(gòu)造第二隆起帶，東部為山地丘陵，西南方向?yàn)閾犴樑璧?，西部融入遼河平原［1］。研究區(qū)內(nèi)多條斷層縱橫交錯(cuò)，嚴(yán)重破壞了巖土體的完整性和連續(xù)性。受到煤層開(kāi)采的影響，研究區(qū)內(nèi)的主要斷層長(zhǎng)期處于不穩(wěn)定狀態(tài)，導(dǎo)致研究區(qū)內(nèi)礦震災(zāi)害頻發(fā)，給礦區(qū)的安全生產(chǎn)帶來(lái)了極大的隱患［2］。如：2001年1月12日ML 2.8級(jí)礦震造成2人死亡、23人受傷；2002年1月26日ML 3.7級(jí)礦震造成井下900 m左右采空區(qū)破壞嚴(yán)重、10余人受傷；2003年5月3日由于沖擊地壓引發(fā)礦震，造成3人死亡、9人受傷；2009年12月11日ML 2.3級(jí)礦震造成1人死亡。老虎臺(tái)礦區(qū)的微震活動(dòng)是影響礦區(qū)安全生產(chǎn)的重要因素之一，通過(guò)分析微震的時(shí)空特征，可以確定礦區(qū)的應(yīng)力集中和主要發(fā)震位置，是一種有效的區(qū)域性礦震監(jiān)測(cè)和預(yù)警手段［34］。最新的研究［5］表明，微震有助于我們探測(cè)一些看不到的地下構(gòu)造。高精度微震事件的時(shí)空位置可以精確刻畫(huà)地下精細(xì)結(jié)構(gòu)，從而推斷地下應(yīng)力狀態(tài)和流體運(yùn)移情況，對(duì)區(qū)域地震危險(xiǎn)性做出研判［68］。目前，微震監(jiān)測(cè)已在多個(gè)方面展開(kāi)應(yīng)用，例如利用微震的空間位置精細(xì)刻畫(huà)地下斷層產(chǎn)狀，也可以利用微震監(jiān)測(cè)區(qū)域水壓致裂情況［910］。而在礦區(qū)，微震活動(dòng)可以很好地反映不同煤層深度的地質(zhì)構(gòu)造情況［1112］。1738年，英國(guó)南斯坦福煤田是世界上首次記錄到有礦震現(xiàn)象的礦區(qū)。1908年，世界上第一個(gè)監(jiān)測(cè)礦井活動(dòng)的地震臺(tái)在德國(guó)波鴻的魯爾區(qū)建立，一直運(yùn)行至20世紀(jì)30年代［13］。1933年，我國(guó)首次在撫順地區(qū)開(kāi)展了礦震的相關(guān)研究［14］。老虎臺(tái)煤礦是撫順地區(qū)現(xiàn)今開(kāi)采的重點(diǎn)煤礦之一，近年來(lái)隨著礦區(qū)開(kāi)采深度與開(kāi)采強(qiáng)度的增加，礦震也呈現(xiàn)出震級(jí)增大且活動(dòng)性增強(qiáng)的趨勢(shì)。礦震的頻發(fā)給礦區(qū)的安全生產(chǎn)帶來(lái)巨大威脅［15］，因此礦區(qū)地震活動(dòng)研究對(duì)有效預(yù)防礦震災(zāi)害，保障人類(lèi)生命財(cái)產(chǎn)安全具有重要意義。

2019年11月12日10時(shí)04分，遼寧撫順老虎臺(tái)礦區(qū)附近（41.86°N，123.94°E）發(fā)生了ML 2.4級(jí)地震。該地區(qū)地質(zhì)構(gòu)造復(fù)雜，斷層易于活動(dòng)［1617］，但周?chē)卣鹋_(tái)網(wǎng)相對(duì)稀疏，導(dǎo)致對(duì)該地區(qū)的地震監(jiān)測(cè)能力不足。我們于2019年11月3—25日在老虎臺(tái)礦區(qū)布置兩條線(xiàn)性密集臺(tái)陣［18］，處于此次微震活動(dòng)的發(fā)震區(qū)間，為深入分析微震的發(fā)震規(guī)律提供了重要的數(shù)據(jù)基礎(chǔ)。本研究基于此開(kāi)展礦區(qū)微震事件的雙差定位研究，以增強(qiáng)我們對(duì)礦區(qū)微震活動(dòng)規(guī)律和孕震風(fēng)險(xiǎn)的認(rèn)識(shí)，為礦區(qū)的安全生產(chǎn)提供參考依據(jù)。

1" 數(shù)據(jù)與方法

1.1" 臺(tái)站與數(shù)據(jù)

本研究于2019年11月3—25日在遼寧撫順老虎臺(tái)礦區(qū)布置兩條線(xiàn)性密集臺(tái)陣，每條測(cè)線(xiàn)各有30臺(tái)三分量短周期地震計(jì)，共計(jì)60臺(tái)。臺(tái)站間距為400 m，儀器記錄頻帶范圍為0.2～150.0 Hz，采樣率為100 Hz，兩條測(cè)線(xiàn)相距6 km，分別橫跨老虎臺(tái)礦區(qū)和龍鳳礦區(qū)，測(cè)線(xiàn)長(zhǎng)度約為12 km（圖1）。

模板匹配定位法利用模板波形與記錄的連續(xù)地震波形進(jìn)行波形互相關(guān)來(lái)識(shí)別和定位微震事件［1920］。我們以中國(guó)地震臺(tái)網(wǎng)渾河ML 2.4級(jí)地震為模板，以

模板事件為中心，沿經(jīng)度、緯度方向搜索范圍均為0.25°，搜索間距為0.001°，沿深度方向搜索范圍為5 km，搜索間距為0.01 km。通過(guò)模板匹配定位法對(duì)密集臺(tái)陣拾取到的23 d地震資料進(jìn)行掃描，共檢測(cè)到325個(gè)微震事件，其中ML 2.0級(jí)以上的地震2個(gè)，ML 1.0～2.0級(jí)的地震23個(gè)，ML 0.0～1.0級(jí)的地震235個(gè)，ML -0.5～0.0級(jí)的地震65個(gè)。圖2展示出檢測(cè)的ML 1.8級(jí)地震事件波形，其中平均相關(guān)系數(shù)閾值設(shè)為0.65，平均信噪比閾值設(shè)為15，檢測(cè)結(jié)果表明在連續(xù)地震波形中能有效拾取微震事件。

不同波形代表同一微震事件在不同地震臺(tái)站上的匹配結(jié)果?；疑ㄐ螢楸尘斑B續(xù)波形；藍(lán)色波形為檢測(cè)的微震事件波形。以2019年11月9日ML1.8級(jí)事件為示例。

1.2" 地震重定位

雙差地震定位法利用兩個(gè)事件對(duì)的走時(shí)殘差進(jìn)行定位，能有效減少地殼結(jié)構(gòu)引起的誤差，與傳統(tǒng)定位結(jié)果相比，其定位結(jié)果更為精確［2122］。采用雙差地震定位法對(duì)通過(guò)模板匹配定位法拾取的325個(gè)微震事件進(jìn)行了重定位。參數(shù)設(shè)置如下：事件對(duì)與臺(tái)站之間的最大距離為9 km，事件對(duì)之間的最大距離為10 km，每個(gè)事件對(duì)的鄰居上限不超過(guò)45個(gè)，每個(gè)鄰居所需要的最小震相對(duì)數(shù)目不少于15對(duì)，每個(gè)事件使用的最少震相不少于4個(gè)，同時(shí)為兼顧定位精度與計(jì)算效率，設(shè)置最多震相不多于35個(gè)。通過(guò)以上參數(shù)將拾取的地震縱（P）、橫（S）波到時(shí)數(shù)據(jù)轉(zhuǎn)變?yōu)槭录?duì)的走時(shí)數(shù)據(jù)。然后采用共軛梯度法進(jìn)行反演，得到295個(gè)重定位事件（圖3）。在重定位過(guò)程中，一維P波速度模型參考了該地區(qū)的人工地震剖面結(jié)果［23］，并依據(jù)煤層氣儲(chǔ)層地震波速度對(duì)該地區(qū)淺層P波速度進(jìn)行修正［24］，建立了重定位使用的速度模型（表1）。

2" 結(jié)果

重定位結(jié)果顯示，本次微震事件最大震級(jí)地震（ML 2.4級(jí)）位于渾河斷裂附近，震源深度0.5 km（圖3）。緯度、經(jīng)度和深度方向上經(jīng)過(guò)20次迭代的平均誤差分別為8.1、8.7和7.8 m，走時(shí)殘差由0.18 s降為0.01 s（圖4）。重定位后走時(shí)殘差明顯減小，震源相對(duì)位置的定位精度有較大提升。重定位后微震事件主要分布在渾河斷裂南部的老虎臺(tái)礦

區(qū)和F1、F2兩條斷層附近，而在渾河斷裂北側(cè)地震事件偏少（圖3a）。重定位后的震源深度統(tǒng)計(jì)結(jié)果表明，地震事件主要集中在0.2～1.5 km范圍內(nèi)（圖3d），這與老虎臺(tái)煤礦目前的開(kāi)采深度相似。

為了更好地展示微震事件的時(shí)序變化，以渾河ML 2.4級(jí)地震為基準(zhǔn)，按不同時(shí)間段展示微震事件的分布情況（圖5）。第一階段為渾河地震發(fā)生前2～10 d，微震事件在渾河斷裂兩側(cè)均有分布，但分布比較彌散、不集中。第二階段為渾河地震的臨近階段（前2 d），微震事件主要分布在距離渾河地震較

近的老虎臺(tái)礦區(qū)附近，呈現(xiàn)向渾河斷裂西部遷移的特征，且ML 1.0～2.0級(jí)的微震事件顯著增加。第三階段為渾河地震發(fā)生后2 d，該階段微震事件明顯較少，主要分布在距離渾河地震較近的位置，呈現(xiàn)逐漸向礦區(qū)移動(dòng)的趨勢(shì)。第四階段為渾河地震發(fā)生后2～14 d，微震事件主要分布在渾河斷裂南部的老虎臺(tái)礦區(qū)和龍鳳礦區(qū)，而在渾河斷裂北部基本沒(méi)有地

震事件。整體表現(xiàn)出微震事件隨時(shí)間逐漸向礦區(qū)內(nèi)部遷移的特征。

渾河ML 2.4級(jí)地震發(fā)生前后呈現(xiàn)出兩種不同的時(shí)空分布特征，渾河地震發(fā)生前微震事件震源深度普遍較深，渾河地震發(fā)生后微震事件震源深度較淺（圖6）。為了更清晰地觀(guān)察微震事件與渾河斷裂的關(guān)系，選取渾河斷裂兩側(cè)2 km范圍內(nèi)的微震事件，沿渾河斷裂走向方向繪制微震事件與渾河斷裂的位置關(guān)系圖（圖7），從圖7中可以看出，微震事件中震級(jí)較大的地震事件均分布在渾河斷裂附近，而震級(jí)較小的地震事件則集中分布在距離渾河斷裂較遠(yuǎn)的礦區(qū)位置，并呈現(xiàn)隨時(shí)間向渾河斷裂南部遷移的特征。

3" 討論

速度模型對(duì)雙差地震定位結(jié)果具有重要影響［25］。張冰等［18］基于區(qū)域一維P波速度模型的定位結(jié)果表明，沿緯度、經(jīng)度及深度方向相對(duì)誤差分別為44、46、61 m，走時(shí)殘差為0.31 s。我們依據(jù)小尺度速度模型的定位結(jié)果整體誤差均得到了明顯收斂，沿緯度、經(jīng)度及深度的相對(duì)誤差分別降低到8.1、8.7、7.8 m，走時(shí)殘差降低為0.01 s。以上兩種不同速度模型的定位誤差對(duì)比表明，本文所采用的小尺度速度模型使雙差定位精度得到進(jìn)一步提升。微震事件的重定位結(jié)果表明震源深度集中在0.2～

1.5 km之間（圖3），這可能與礦區(qū)開(kāi)采活動(dòng)具有一定的相關(guān)性。目前的煤礦開(kāi)采主要集中于老虎臺(tái)礦區(qū)，西南側(cè)的撫順西露天礦已停止開(kāi)采活動(dòng)。煤礦開(kāi)采活動(dòng)可能影響地下介質(zhì)的應(yīng)力狀態(tài)，使其介質(zhì)強(qiáng)度較低，而應(yīng)力更傾向于在介質(zhì)強(qiáng)度較低的區(qū)域積累，從而導(dǎo)致本次的微震活動(dòng)主要集中于正在開(kāi)采的老虎臺(tái)礦區(qū)。目前，老虎臺(tái)礦區(qū)的開(kāi)采深度達(dá)到1 000 m左右，與此次微震事件的優(yōu)勢(shì)發(fā)震深度一致［2627］。這種微震事件的空間相關(guān)性可能一定程度上反映了老虎臺(tái)礦區(qū)的開(kāi)采部位和開(kāi)采深度特征。

一般認(rèn)為微小斷層滑動(dòng)是引起礦震的重要因素之一［28］，與采礦過(guò)程中引起的震動(dòng)在極性、振幅和頻率上有所差異。斷層滑動(dòng)可能受到區(qū)域擠壓或拉張應(yīng)力的作用，初至波形在不同方位的極性可能向上或向下；而采礦過(guò)程中的微震主要是由采空區(qū)上

方受重力影響向下塌落引起，因此不同方位記錄到的極性一般為向下［29］。在振幅特征上，斷層滑動(dòng)導(dǎo)致的微震振幅通常較大，而采礦活動(dòng)引起的微震振幅則較小［30］。在頻率特征上，斷層滑動(dòng)導(dǎo)致的微震頻帶較寬，而采礦活動(dòng)引起的微震頻帶則較窄且主要以低頻信號(hào)為主［31］。通常認(rèn)為，0＜MLlt;1.5 的地震屬于采礦直接型礦震，而ML≥1.5的地震則屬于采礦間接型礦震［32］。本次微震事件的重定位結(jié)果表明，渾河ML 2.4級(jí)地震發(fā)生后，MLgt;1.5的微震事件顯著增加，并且呈現(xiàn)隨時(shí)間向礦區(qū)內(nèi)部遷移的特征（圖6）。推斷渾河地震的發(fā)生可能破壞了該區(qū)域原有的應(yīng)力平衡狀態(tài)，致使礦區(qū)內(nèi)部出現(xiàn)應(yīng)力擾動(dòng)，導(dǎo)致礦區(qū)出現(xiàn)斷層發(fā)生活化現(xiàn)象，進(jìn)而引起大量微震的產(chǎn)生。地質(zhì)調(diào)查顯示老虎臺(tái)礦區(qū)內(nèi)部存在多條微小斷裂，而微震事件一般在斷裂附近保持較高的活動(dòng)水平［33］，表明先存斷層更易引發(fā)剪切運(yùn)動(dòng)，從而誘發(fā)微震活動(dòng)［34］。震源機(jī)制解研究結(jié)果［18］表明，渾河ML 2.4級(jí)地震可能使礦區(qū)整體處于拉張狀態(tài)，為斷層滑動(dòng)提供了良好條件。先前的研究表明，采礦活動(dòng)積累的能量沿?cái)鄬觾?yōu)先釋放，在開(kāi)采過(guò)程中微小的應(yīng)力變化可能會(huì)導(dǎo)致現(xiàn)有斷層規(guī)模的擴(kuò)大或新裂縫的產(chǎn)生［3536］，從而使微震活動(dòng)維持在較高水平。結(jié)合本次微震事件的分布特征，推斷由渾河地震引起的礦區(qū)微小斷層滑動(dòng)可能是此次微震事件產(chǎn)生的主要因素。

渾河ML 2.4級(jí)微震事件的時(shí)序變化（圖5）表明：渾河地震發(fā)生前2～10 d微震事件在F1、F2斷層附近分布較少；渾河地震發(fā)生前2 d微震主要集中在渾河地震附近，可能暗示著渾河地震附近區(qū)域的應(yīng)力正在不斷積累；渾河地震發(fā)生后微震事件數(shù)量逐漸增加并呈現(xiàn)出向南部礦區(qū)遷移的特征，可能是由于渾河地震的發(fā)生，破壞了礦區(qū)內(nèi)部斷層原有的應(yīng)力平衡狀態(tài)［3738］，斷層失穩(wěn)產(chǎn)生相對(duì)滑移并伴隨能量的釋放，進(jìn)而誘發(fā)了大量微震的產(chǎn)生［3940］。

重定位后的地震空間分布特征與地表斷層及深

部構(gòu)造存在對(duì)應(yīng)關(guān)系［41］。為進(jìn)一步分析主要發(fā)震

剖面位置見(jiàn)圖3a。

斷層的微震演化規(guī)律，分別繪制了沿渾河斷裂、F1、F2走向的三個(gè)微震事件剖面（圖8）。渾河斷裂的剖

面（AA′）橫跨整個(gè)研究區(qū)，微震分布集中在東西兩端，在剖面中部地震事件較少。F1剖面（BB′）位于老虎臺(tái)礦區(qū)內(nèi)部，震源深度多分布在1 km左右，可能是由于老虎臺(tái)煤礦當(dāng)前的開(kāi)采深度對(duì)其的影響。微震事件在F2剖面（CC′）上呈現(xiàn)出很好的連續(xù)性，F(xiàn)2斷層處的小震活動(dòng)逐漸向北擴(kuò)展至渾河斷裂。然而，以往的勘察過(guò)程中并未在地表發(fā)現(xiàn)此斷層。結(jié)合最新微震重定位結(jié)果的展布特征，推斷F2斷層可能是延伸至渾河斷裂的NNW向隱伏斷層。

4" 結(jié)論

本研究基于2019年11月3—25日在遼寧撫順老虎臺(tái)礦區(qū)布設(shè)的密集臺(tái)陣數(shù)據(jù)，利用雙差地震定位方法對(duì)拾取到的325個(gè)地震事件進(jìn)行了重定位，主要獲得了以下認(rèn)識(shí)：

1）重定位結(jié)果表明，微震事件主要分布在渾河斷裂南側(cè)以及F1、F2兩條斷層附近，而在北側(cè)的地震事件偏少，整體微震事件主要集中在0.2～1.5 km范圍內(nèi)，與采礦活動(dòng)密切相關(guān)。

2）渾河地震發(fā)生之前，微震事件在渾河斷裂南北兩側(cè)均有分布，而在渾河地震發(fā)生后，微震事件主要分布在渾河斷裂南部的老虎臺(tái)礦區(qū)內(nèi)部，推斷渾河地震可能使礦區(qū)內(nèi)出現(xiàn)應(yīng)力擾動(dòng)，導(dǎo)致礦區(qū)內(nèi)斷層發(fā)生活化現(xiàn)象，進(jìn)而引起大量微震的產(chǎn)生。

3）礦區(qū)內(nèi)部F2斷層附近的微震活動(dòng)逐漸向北擴(kuò)展達(dá)到渾河斷裂，推斷F2斷層可能是一條延伸至渾河斷裂的NNW向隱伏斷層。

參考文獻(xiàn)（References）：

［1］" Chen D， Wang E Y， Li N. Study on the Source Parameters of the Micro-Earthquakes in Laohutai Ccoal Mine Based on Double Difference Relocation［J］. Soil Dynamics and Earthquake Engineering， 2021， 142： 106540.

［2］" 劉招君，孟慶濤，柳蓉，等. 撫順盆地始新統(tǒng)計(jì)軍屯組油頁(yè)巖地球化學(xué)特征及其地質(zhì)意義［J］.巖石學(xué)報(bào)，2009，25（10）：23402350.

Liu Zhaojun， Meng Qingtao， Liu Rong， et al. Geochemical Characteristics of Oil Shale of Eocene Jijuntun Formation and Its Geological Significance， Fushun Basin［J］. Acta Petrologica Sinica， 2009， 25（10）： 23402350.

［3］" Xuan D Y， Xu J L， Zhu W B. Dynamic Disaster Control Under a Massive Igneous Sill by Grouting from Surface Boreholes［J］. International Journal of Rock Mechanics and Mining Sciences， 2014， 71： 176187.

［4］" 竇林名，何江，鞏思園，等. 采空區(qū)突水動(dòng)力災(zāi)害的微震監(jiān)測(cè)案例研究［J］.中國(guó)礦業(yè)大學(xué)學(xué)報(bào)，2012，41（1）：2025.

Dou Linming， He Jiang， Gong Siyuan， et al.A Case Study of Micro-Seismic Monitoring： Goaf Water-Inrush Dynamic Hazards［J］. Journal of China University of Mining and Technology， 2012， 41（1）： 2025.

［5］" 楊宏偉，王霽川，孔慶豐，等. 井中地震粘聲逆時(shí)偏移成像影響因素分析［J］.物探與化探，2022，46（4）：877886.

Yang Hongwei， Wang Jichuan， Kong Qingfeng， et al. An Analysis of Influencing Factors of Visco-Acoustic Reverse Time Migration Imaging in Borehole Seismic［J］. Geophysical and Geochemical Exploration， 2022， 46（4）： 877886.

［6］" Ross Z E， Cochran E S， Trugman D T， et al. 3D Fault Architecture Controls the Dynamism of Earthquake Swarms［J］. Science， 2020， 368： 13571361.

［7］" 范廣偉，Wallace T C，趙大鵬，等. 應(yīng)用地震層析成象技術(shù)研究深部速度結(jié)構(gòu)［J］.CT理論與應(yīng)用研究，1999，8（1）：18.

Fan Guangwei， Wallace T C， Zhao Dapeng， et al. A Seismic Tomography Study of Deep Velocity Structure［J］. CT Theory and Application， 1999， 8（1）： 18.

［8］" 李洪麗，劉財(cái)，田有，等. 中國(guó)東北地殼結(jié)構(gòu)的地震層析成像［J］.吉林大學(xué)學(xué)報(bào)（地球科學(xué)版），2022，52（1）：270280.

Li Hongli， Liu Cai， Tian You， et al. Seismic Tomography of the Crustal Structure in Northeast China［J］. Journal of Jilin University （Earth Science Edition）， 2022， 52（1）： 270280.

［9］" Drew J， White R S， Tilmann F， et al. Coalescence Microseismic Mapping［J］. Geophysical Journal International， 2013， 195（3）： 17731785.

［10］" Maity D， Ciezobka J. Using Microseismic Frequency-Magnitude Distributions from Hydraulic Fracturing as an Incremental Tool for Fracture Completion Diagnostics［J］. Journal of Petroleum Science and Engineering， 2019， 176（C）： 11351151.

［11］" 黃宇奇，查華勝，高級(jí)，等. 基于密集臺(tái)陣地震背景噪聲成像預(yù)測(cè)煤礦瓦斯分布［J］.地球物理學(xué)報(bào)，2021，64（11）： 39974011.

Huang Yuqi，Zha Huasheng， Gao Ji， et al. Predicting the Distribution of Coalbed Methane by Ambient Noise Tomography with a Dense Seismic Array［J］. Chinese Journal of Geophysics， 2021， 64（11）： 39974011.

［12］" 張昭，張龍飛，段剛，等. 基于正演模擬的煤層采空區(qū)地震響應(yīng)特征［J］.CT理論與應(yīng)用研究，2021，30（3）：291300.

Zhang Zhao， Zhang Longfei， Duan Gang， et al.Seismic Response Characteristics of Coal Goaf Based on Forward Modeling［J］. CT Theory and Application， 2021， 30（3）： 291300.

［13］" Hassani H， Hlouek F， Alexandrakis C， et al. Migration-Based Microseismic Event Location in the Schlema-Alberoda Mining Area［J］. International Journal of Rock Mechanics and Mining Sciences， 2018， 110： 161167.

［14］" 潘一山，李忠華，章夢(mèng)濤. 我國(guó)沖擊地壓分布、類(lèi)型、機(jī)理及防治研究［J］.巖石力學(xué)與工程學(xué)報(bào)，2003，22（11）：18441851.

Pan Yishan， Li Zhonghua， Zhang Mengtao.Distribution， Type， Mechanism and Prevention of Rockbrust in China［J］. Chinese Journal of Rock Mechanics and Engineering， 2003， 22（11）： 18441851.

［15］" 李永靖，張向東，袁世君，等. 撫順老虎臺(tái)礦礦山地震活動(dòng)趨勢(shì)研究［J］.中國(guó)礦業(yè)，2004，13（11）：6567.

Li Yongjing， Zhang Xiangdong， Yuan Shijun， et al. Study of Fushun Laohutai Mineral Earthquake Activity Trend［J］. China Mining Magazine， 2004， 13（11）： 6567.

［16］" 劉財(cái)，鄭確，田有，等. 遼寧海城及其鄰區(qū)地震活動(dòng)性研究［J］.吉林大學(xué)學(xué)報(bào)（地球科學(xué)版），2014，44（4）：13481356.

Liu Cai， Zheng Que， Tian You， et al. Study on the Activity of the Earthquakes of Haicheng Region in Liaoning Province［J］. Journal of Jilin University （Earth Science Edition）， 2014， 44（4）： 13481356.

［17］" 張?zhí)锾?，楊為民，萬(wàn)飛鵬. 渾河斷裂帶地質(zhì)災(zāi)害發(fā)育特征及其成因機(jī)制［J］.吉林大學(xué)學(xué)報(bào)（地球科學(xué)版），2022，52（1）：149161.

Zhang Tiantian， Yang Weimin， Wan Feipeng. Characteristics and Formation Mechanism of Geohazards in Hunhe Fault Zone［J］. Journal of Jilin University （Earth Science Edition）， 2022， 52（1）： 149161.

［18］" 張冰，張廣偉，焦明若，等. 基于密集臺(tái)陣研究2019遼寧撫順M2.4礦震震源參數(shù)［J］.地球物理學(xué)報(bào)，2021，64（4）：12271235.

Zhang Bing， Zhang Guangwei， Jiao Mingruo， et al.Source Characterization of the 2019 M2.4 Earthquakes Induced by Fushun， Liaoning Mining Based on a Dense Seismic Array［J］. Chinese Journal of Geophysics， 2021， 64（4）： 12271235.

［19］" Zhang M， Wen L X. Earthquake Characteristics Before Eruptions of Japan’s Ontake Volcano in 2007 and 2014［J］. Geophysical Research Letters， 2015， 42（17）： 69826988.

［20］" 張淼. 地震定位和檢測(cè)［D］.合肥：中國(guó)科學(xué)技術(shù)大學(xué)，2015.

Zhang Miao. Earthquake Location and Detection［D］. Hefei： University of Science and Technology of China， 2015.

［21］" Waldhauser F， Ellsworth W L.A Double-Difference Earthquake Location Algorithm： Method and Application to the Northern Hayward Fault， California［J］. Bulletin of the Seismological Society of America， 2000， 90（6）： 13531368.

［22］" Waldhauser F， Ellsworth W L. Fault Structure and Mechanics of the Hayward Fault， California， from Double-Difference Earthquake Locations［J］. Journal of Geophysical Research： Solid Earth， 2002， 107（B3）： 115.

［23］" 王兆國(guó)，劉財(cái)，馮晅，等. 中國(guó)東北地區(qū)地震空間分布與主要斷裂帶、深部構(gòu)造及應(yīng)力場(chǎng)關(guān)系［J］.世界地質(zhì)，2009， 28（4）：513519.

Wang Zhaoguo， Liu Cai， Feng Xuan， et al.Earthquake Space Distribution and Its Relationships with Main Faults， Deep Structure and Stress Field in Northeast China［J］. World Geology， 2009， 28（4）： 513519.

［24］" 唐志遠(yuǎn)，楊延強(qiáng)，劉國(guó)昌. 基于雙重孔隙結(jié)構(gòu)的煤層氣儲(chǔ)層震波速度預(yù)測(cè)研究［J］.煤炭科學(xué)技術(shù)，2018，46（12）：187191.

Tang Zhiyuan， Yang Yanqiang， Liu Guochang. Study on Prediction on Seismic Wave Velocity of Coalbed Methane Reservoir Based on Dual Pore Structure［J］. Coal Science and Technology， 2018， 46 （12）： 187191.

［25］" 王勤彩，王中平，張金川，等. 2010年4月玉樹(shù)MS 7.3地震序列的斷層結(jié)構(gòu)［J］.地球物理學(xué)報(bào)，2015，58（6）：19311940.

Wang Qincai， Wang Zhongping， Zhang Jinchuan， et al.Fault Structure of Ms7.3 Yushu Earthquake Sequence in April， 2010［J］. Chinese Journal of Geophysics， 2015， 58（6）： 19311940.

［26］" Li T， Cai M F， Cai M. A Review of Mining-Induced Seismicity in China［J］. International Journal of Rock Mechanics and Mining Sciences， 2007， 44（8）： 11491171.

［27］" 李文健. 微震監(jiān)測(cè)技術(shù)在沖擊地壓礦井的應(yīng)用［J］.中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào)，2015，26（4）：116120.

Li Wenjian. Application of Microseismic Monitoring Technology in Rock Burst Mine［J］. The Chinese Journal of Geological Hazard and Control， 2015， 26（4）： 116120.

［28］" Cai M F， Ji H G， Wang J A. Study of the Time-Space-Strength Relation for Mining Seismicity at Laohutai Coal Mine and Its Prediction［J］. International Journal of Rock Mechanics amp; Mining Sciences， 2005， 42（1）： 145151．

［29］" Meng S J， Mu H L， Wang M， et al. Characteristics and Identification Method of Natural and Mine Earthquakes： A Case Study on the Hegang Mining Area［J］. Minerals， 2022， 12（10）： 1256.

［30］" 和雪松，李世愚，沈萍，等. 用小波包識(shí)別地震和礦震［J］.中國(guó)地震，2006，22（4）：425434.

He Xuesong， Li Shiyu， Shen Ping， et al. A Wavelet Packet Approach to Wave Classification of Earthquakes and Mining Shocks［J］. Earthquake Research in China， 2006， 22（4）： 425434.

［31］" 張思萌. 鶴崗地區(qū)地震、爆破與礦震記錄的識(shí)別［J］.地震地磁觀(guān)測(cè)與研究，2021，42（6）：8288.

Zhang Simeng. Identification of Earthquake， Blasting and Mining Seismic Records in the Hegang Region［J］. Seismological and Geomagnetic Observation and Research， 2021， 42（6）： 8288.

［32］" 李鐵，張建偉，呂毓國(guó)，等. 采掘活動(dòng)與礦震關(guān)系［J］.煤炭學(xué)報(bào)，2011，36（12）：21272132.

Li Tie， Zhang Jianwei， Lü Yuguo， et al. Relationship Between Mining and Mining-Induced Seismicty［J］. International Journal of Coal Science， 2011， 36（12）： 21272132.

［33］" 孫曉元，楊威，張飛，等. 老虎臺(tái)礦55002# 掘進(jìn)面通過(guò)F71斷層前后微震事件的時(shí)空特征分析［J］.中國(guó)安全生產(chǎn)科學(xué)技術(shù)，2011，7（4）：7174.

Sun Xiaoyuan， Yang Wei， Zhang Fei， et al.Space-Time Characteristics Analysis of Microearthquake Events Before and After F71 Fault on 55002 Heading Face in Laohutai Colliery［J］. China Occupational Safety and Health， 2011， 7（4）： 7174.

［34］" 唐禮忠，翦英驊，李地元，等. 基于微震矩張量的礦山圍巖破壞機(jī)制分析［J］.巖土力學(xué)，2017，38（5）：14361444.

Tang Lizhong， Jian Yinghua， Li Diyuan， et al. Analysis of Damage Mechanism for Surrounding Rock Based on Microseismic Moment Tensor［J］. Rock and Soil Mechanics， 2017， 38（5）： 14361444.

［35］" Liu J P， Liu Z S， Wang S Q， et al. Analysis of Microseismic Activity in Rock Mass Controlled by Fault in Deep Metal Mine［J］. International Journal of Mining Science and Technology， 2016， 26（2）： 235239.

［36］" Liu J P， Feng X T， Li Y H， et al. Studies on Temporal and Spatial Variation of Microseismic Activities in a Deep Metal Mine［J］. International Journal of Rock Mechanics and Mining Sciences， 2013， 60： 171179.

［37］" 郭曉強(qiáng)，竇林名，陸菜平，等. 采動(dòng)誘發(fā)斷層活化的微震活動(dòng)規(guī)律研究［J］.煤礦安全，2011，42（1）：2630.

Guo Xiaoqiang， Dou Linming， Lu Caiping， et al.Research on the Microseismic Activity of Fault Reaction Induced by Coal Mining［J］. Safety in Coal Mines， 2011， 42 （1）： 2630.

［38］" 余傳濤，劉鴻福，于艷梅，等. CSAMT法在煤礦隱伏斷層探測(cè)中的應(yīng)用［J］.CT理論與應(yīng)用研究，2010，19（1）：2833.

Yu Chuantao， Liu Hongfu， Yu Yanmei， et al. Application of CSAMT Method in Buried Fault Detecting in Coal Mine［J］. CT Theory and Application， 2010， 19（1）： 2833.

［39］" 胡剛，段寶平，何正勤，等. 利用井間地震層析成像方法探測(cè)隱伏斷層［J］.CT理論與應(yīng)用研究，2015，24（3）：345355.

Hu Gang， Duan Baoping， He Zhengqin， et al. Using Cross-Well Seismic Tomography to Detect Hidden Faults［J］. CT Theory and Application， 2015， 24（3）： 345355.

［40］" 呂苗苗，徐小明. 臺(tái)陣地震學(xué)方法及對(duì)間斷面結(jié)構(gòu)的研究進(jìn)展［J］.CT理論與應(yīng)用研究，2014，23（6）：10411050.

Lü Miaomiao， Xu Xiaoming. Array Seismological Methods and Its Progress in the Discontinuity Structure Research［J］. CT Theory and Application， 2014， 23（6）： 10411050.

［41］" 李欣蔚，張廣偉，謝卓娟，等. 2021年四川瀘縣M6.0地震發(fā)震機(jī)理及地震活動(dòng)時(shí)空演化特征［J］.地球物理學(xué)報(bào)，2022，65（11）：42844298.

Li Xinwei， Zhang Guangwei， Xie Zhuojuan， et al.Seismogenic Mechanism of the 2021 M6.0 Luxian Earthquake and Seismicity Spatio-Temporal Characteristics Around the Source Region［J］. Chinese Journal of Geophysics， 2022， 65（11）： 42844298.