王璽 李施慶 李秋濤 程力 劉煥新

摘要：為全面了解不同震源機制下地震波傳播和分布規(guī)律及其對巷道圍巖的影響，基于現(xiàn)場監(jiān)測數(shù)據(jù)研究了不同震源參數(shù)下PGV的分布規(guī)律。研究結果表明，震源距離、矩震級及斷層滑移角都會對PGV的分布產(chǎn)生不同程度的影響。通過比較經(jīng)驗公式計算PGV值與實測值間的擬合相關系數(shù)，判斷各預測方法適用性。研究發(fā)現(xiàn)Sadov預測公式對爆破事件的PGV預測值與實測結果相關度更高。依據(jù)研究結果，計算了動態(tài)震源機制下PGV空間分布，并依此給出了不同震源位置和震級條件下的圍巖支護建議。

關鍵詞：金屬礦山;震源;震動監(jiān)測;PGV;預測公式;擬合相關系數(shù)

中圖分類號：TD853.34文章編號：1001-1277（2020）10-0029-08

文獻標志碼：Adoi：10.11792/hj20201006

引 言

隨著地表淺部資源的開采殆盡，金屬礦開采開始轉向深部。而由于地下金屬礦復雜的開采技術條件及多變的礦巖受力狀態(tài)，其采場內(nèi)圍巖會產(chǎn)生不同程度的開采擾動，同時也導致地下金屬礦采場中存在多種高危的安全隱患[1]。在地下金屬礦開采過程中，若對礦山井下強擾動的影響范圍、具體位置、能量大小等缺乏準確的了解和認知，礦山工人與生產(chǎn)設備將持續(xù)處于極其危險且頻率高發(fā)的災害威脅中，并有可能引發(fā)重大的財產(chǎn)損失與人員傷亡，最終帶來難以估量的損失。

在礦山震動的研究工作中，預測震動強度對震動安全的判據(jù)有著重要意義。震動強度（I）的預測方法依據(jù)表征震動強度物理量的不同而對應不同的形式，通常各類預測方法都符合以下函數(shù)形式[2]：

I=f（m，n，k）（1）

式中：m為震源;n為傳播途徑;k為儀器特性。

通過公式（1）發(fā)現(xiàn)，影響震動強度的因素多樣且復雜，在計算過程中需要對多因素的影響進行綜合考量[3-4]。由于通過多變量確認函數(shù)形式的難度很高，該過程一般會選擇某些獨立變量作為常數(shù)，在此基礎上對其他變量進行研究，找到它們之間存在的函數(shù)關系。運用公式（1）及現(xiàn)場監(jiān)測所得的震動相關數(shù)據(jù)能夠對質點峰值震速、加速度、位移等參數(shù)進行預測，再由預測結果對震動強度進行評判。實際工作中運用最廣泛的是對PGV（Peak Ground Velocity，PGV）進行預測。但是，在運用PGV經(jīng)驗公式對實測數(shù)據(jù)進行回歸分析時，計算得出的回歸系數(shù)離散性往往較大。所以經(jīng)驗公式在實際運用時會出現(xiàn)較多問題和局限性。為此，研究者結合現(xiàn)有經(jīng)驗將公式進行變形，得出了很多適用于各類情況的相關公式。例如：陳壽如等[5]在計算露天礦爆破產(chǎn)生的震動模型時加入了高差這一變量，得到了新的變形公式;韓子榮[6]基于對主頻與結構自震頻率影響的考慮得出了折合速度公式。在預測PGV時除運用傳統(tǒng)線性回歸方法以外，近些年研究者研發(fā)出了不少新的預測方法，比較典型的有：王民壽等[7]通過對雙隨機變量進行回歸分析來預測震速;徐全軍等[8-9]通過神經(jīng)網(wǎng)絡對震動峰值進行預測;黃光球等[10]通過遺傳規(guī)劃法對震動峰值進行預測。

綜合上述研究，針對震動事件產(chǎn)生的應力波，采礦中經(jīng)常使用PGV描述圍巖在應力波作用下的能量平衡狀態(tài)（能量需求與能量吸收能力）和動態(tài)變化，由此對圍巖的穩(wěn)定性和完整性作出評估[11]。因此，本次研究基于現(xiàn)場震動監(jiān)測數(shù)據(jù)，對礦震數(shù)據(jù)在震源參數(shù)影響下的分布特征進行研究;并選取某礦震事件的震源參數(shù)，運用國內(nèi)外已有的PGV計算公式進行預測計算;通過研究計算結果與實測值間的擬合相關系數(shù)，判斷不同的計算公式在該事件PGV預測的適用性;最后運用優(yōu)選計算公式的預測結果給出巖體支護建議。

1 圍巖震動測定與分析

1.1 震源參數(shù)與PGV測定

山東金洲礦業(yè)集團有限公司金青頂?shù)V區(qū)（下稱“金青頂?shù)V區(qū)”）坐落在黃壘河南岸的金青頂上，位于山東省乳山市下初鎮(zhèn)。金青頂?shù)V區(qū)整體屬于硬巖礦山，最大開采深度已超過1 000 m，大部分采場和巷道布置集中、開挖空間較小，且眾多采場與井巷工程的開挖空間結構參數(shù)和圍巖性質不盡相同，這些因素導致該礦區(qū)的震動事件錯綜復雜，嚴重影響礦區(qū)范圍內(nèi)礦震監(jiān)測和分析工作。

為解決該問題，需要對PGV在各項震源參數(shù)下的分布情況進行研究。利用局域微震監(jiān)測系統(tǒng)（見圖1）對該礦區(qū)-1 145 m中段E8、E9、E10及E11穿脈內(nèi)的活動事件進行了跟蹤記錄，記錄信息包括事件時間、矩震級、實測PGV、傳感器距震源距離及相對位置等參數(shù)。將記錄的監(jiān)測數(shù)據(jù)與震源特征參數(shù)繪制成圖，即可得到PGV與相關震源參數(shù)的關系圖。本文分別討論了震源距離、矩震級及斷層滑移角對PGV分布的影響情況。

1.2 震源參數(shù)對PGV分布影響

1）震源距離對PGV的影響。隨機抽取8次震級不一的事件，繪制不同矩震級（Mw）下震源距離（d）與PGV（v）的關系，見圖2。

從圖2可以看出：隨著震源距離的增大，震動的峰值速度都呈指數(shù)型下降趨勢，說明傳播過程中，巖石介質將對PGV產(chǎn)生明顯的削弱作用。此外，PGV的衰減程度與震源距離也有著密切的聯(lián)系，當震源距離d<300 m時，事件8（Mw=0.542 397 9）PGV隨震源距離的衰減率為k1=4.776 784×10-6;當300 m≤d≤500 m時，衰減率k2=2.772 07×10-7;當d>500 m時，衰減率k3=5.883 29×10-8。由此可見，整體上PGV在傳播過程中的衰減作用受震源距離的影響，在距震源距離較小時，震源距離對PGV的衰減影響較大;隨著震源距離增大，PGV的衰減受震源距離的影響作用逐漸弱化。

2）矩震級對PGV的影響。矩震級是利用地震矩的大小確定震級，它由地震斷層的破裂面積、平均錯動量及巖石剪切模量的乘積來確定，是一個描述震動發(fā)生時力學強度的物理量?；诒O(jiān)測數(shù)據(jù)，取震源距離為100～200 m、200～300 m、300～400 m、400～500 m及500 m以上的5種情況下矩震級和PGV數(shù)值，分別繪制PGV（v）相對于矩震級（Mw）大小的散點分布圖，見圖3。

從圖3可以看出：矩震級增大時，PGV值高的散點數(shù)量增多，整體上PGV呈現(xiàn)增大趨勢;分析發(fā)現(xiàn)其原因可能是，隨著矩震級的增大，斷層的長度等尺寸參數(shù)變大，使得PGV隨矩震級增大也顯著變大[12]。但是，由PGV-Mw變化趨勢公式（見表1）可知：在不同震源距離下，PGV隨矩震級增大的程度有所差異，在近震源情況下，PGV受矩震級影響較大，隨矩震級呈近指數(shù)趨勢增長;但當監(jiān)測點逐漸遠離震源時，PGV受矩震級的影響減小，增大趨勢逐漸減弱。由此說明，矩震級對PGV的分布有較大影響，且該作用效果隨震源距離的由近至遠逐漸衰減。

3）斷層滑移角對PGV的影響。在金屬礦山中，斷層滑移角通常指礦體上下盤相對運動的方向（斷層滑移方向）與斷層的走向之間的夾角。為探究不

同斷層滑移角對PGV分布的影響，在監(jiān)測所得的數(shù)據(jù)中取滑移角（αr）為0°、45°、90°、-45°和-90°的5次震動事件，分別繪制PGV關于震源距離變化的曲線，見圖4。

從圖4可以看出：不同斷層滑移角下隨著震源距離的增大，所監(jiān)測到的PGV值都在逐漸減小，整體上呈現(xiàn)出隨震源距離增大而衰減的趨勢，但在相同震源距離下，不同滑移角下的PGV值大小存在差異。

（1）整體上看，走滑斷層（αr=0°）下的PGV值最大，且分布較為規(guī)律，其次是正斷層（αr=45°）、垂直斷層（αr=±90°），而逆斷層（αr=-45°）下產(chǎn)生的PGV值最小。

（2）從曲線的走向和趨勢上看，αr=-45°、αr=90°和αr=-90°的3條曲線較為吻合，呈現(xiàn)的規(guī)律性相似，特別是當d>600 m時，3條曲線幾乎重合，說明逆斷層和垂直斷層PGV隨震源距離分布具有一定的相似性，且這種相似性隨與震源距離的增加而愈發(fā)明顯。

2 PGV預測方法

對于爆破作業(yè)引起的震動事件，薩道夫斯基的經(jīng)驗公式[13]，又稱“Sadov公式”，是目前最為公認、應用面最廣的經(jīng)驗公式，其具體表達式如下：

vs=K（Q1/3R）α（2）

式中：vs為測點最大震動速度（cm/s）;K為與地形、爆破方式等因素相關的系數(shù)，堅硬巖石50～150，中等硬度巖石150～250，軟巖250～350;Q為裝藥量（kg）;R為測點距爆心的距離（m）;α為爆破地震波衰減指數(shù)。

董隴軍等[14]根據(jù)Sadov公式的原理，發(fā)現(xiàn)裝藥量對爆破威力的影響與地震動峰值速度反映出的爆破荷載的能量幅值具有相同趨勢，基于上述理論，用震源PGV替代裝藥量Q作為計算參數(shù)，對Sadov公式進行了變形，從而得到新的圍巖爆破震動峰值速度計算公式：

vb=m1v1/30/lm2（3）

式中：vb為質點震動峰值速度（m/s）;v0為震源震動峰值速度（m/s）;l為測點距爆心的距離（m）;m1為受地形、爆破方式影響的相關系數(shù);m2為震動波衰減系數(shù)。

對于非爆破事件引起的礦山震動事件，目前已有的PGV計算通用方程式大多是基于McGarr等人先前的工作：

lg Rvmax=AM+1（4）

式（4）變形可得：

vmax=C10AM/R（5）

式中：vmax為質點峰值速度（m/s）;M為地震時間的震級;A和C為礦山特征參數(shù)。

1996年Kaiser對來自Brunswick礦山、El Teniente礦山和Creighton礦山的地震數(shù)據(jù)總體進行了95 %置信區(qū)間下的回歸分析，綜合McGarr 1984年的研究結果，確定參數(shù)A=0.5和C=0.25。因此，有：

vmax=C10a（Mw+1.5）R≈1.4×10（Mw+1.5）R（6）

式中：a為與礦山相關的特征參數(shù)。

由上述內(nèi)容可知，目前國內(nèi)外對PGV計算方法的研究已經(jīng)較為成熟。但是，對經(jīng)驗公式在不同震動事件下適用性的討論研究尚少，工程和研究人員在選擇計算公式時缺少理論依據(jù)，所以需要結合現(xiàn)場數(shù)據(jù)對上述公式進行計算討論，并將理論計算結果與實際監(jiān)測值進行擬合，根據(jù)擬合的相關系數(shù)來確定不同公式的適用范圍和精確度[15-16]。

為了獲得礦山特征參數(shù)，取Mw=0.053 067 5事件中的震源距離和實測PGV作為數(shù)據(jù)來源（見表2），對公式（5）進行迭代擬合運算，從而獲得適用于金青頂?shù)V區(qū)的特征參數(shù)：C=0.260 2和A=0.355 1。

此外，由于PGV和矩震級由震源能量釋放決定，同時能量與速度的平方呈正相關，所以可將公式（3）

中的震源v0替換為矩震級Mw，得到公式（3）的變形公式：

vb=m1M1/6w/Rym2（7）

式中：Ry為測點距震源的距離（m）。

同樣采用Mw=0.334 718 2事件中的數(shù)據(jù)對公式（7）進行迭代擬合，求得m1、m2的取值為18和1.9。 公式（7）、公式（5）和公式（6）的實際監(jiān)測值、公式計算值與震源距離的關系曲線見圖5（為表述方便，后文統(tǒng)稱公式（?。?、公式（ⅱ）、公式（ⅲ））。

對公式（ⅰ）、公式（ⅱ）、公式（ⅲ）的計算值與其對應的實測值進行相關系數(shù)擬合計算，獲得3個公式的相關系數(shù)分別為r1=0.859、r2=0.798和r3=0.756。 結合實測值、計算值同震源距離的關系圖及二者間的擬合相關系數(shù)，對3種PGV計算公式的擬合情況和適用性進行分析：

1）在3種公式的關系圖中，PGV隨震源距離的變化趨勢皆與擬合曲線相同，說明3種公式中都體現(xiàn)了PGV隨震源距離的衰減規(guī)律。

2）通過比較實測值的擬合曲線與3條計算曲線可發(fā)現(xiàn)，3條曲線在與擬合曲線的吻合度上存在差異，公式（?。┯嬎闱€在d>500 m的區(qū)間內(nèi)與擬合曲線有較高的擬合度，公式（ⅱ）和公式（ⅲ）整體上與擬合曲線有較大差異，說明3種PGV計算公式在當前情況下適用性不一。

3）比較三者的相關系數(shù)：r1=0.859>r2=0.798>r3=0.756，發(fā)現(xiàn)公式（?。┑挠嬎闱€和擬合曲線的吻合程度較高[17]，證明公式（?。┑暮瘮?shù)曲線和PGV的實際分布規(guī)律高度相關，而公式（ⅱ）、公式（ⅲ）的擬合度較低，說明對于金青頂?shù)V區(qū)的震動事件Mw=0.053 067 5而言，公式（?。┑倪m用性優(yōu)于另2種公式。

4）3種計算公式對該事件PGV計算所得結果仍存在較大誤差。

針對3種公式在該事件中體現(xiàn)出的適用性差異，結合公式原理和震動事件本身的特征對誤差進行分析：

（1）基于對原始數(shù)據(jù)的分析，了解到Mw=0.053 067 5 事件的震源以爆破作用為主導，而公式（?。┦怯蒘adov 公式推導變形產(chǎn)生，Sadov公式是以爆破信號為震源的PGV計算公式，其計算原理和有關參數(shù)是基于爆破震動的形成機理和傳播規(guī)律，所以理論上其變形公式更適用于圍巖爆破震動峰值速度的計算[18];而公式（ⅱ）、公式（ⅲ）是基于McGarr經(jīng)驗公式的改進，其本身的參數(shù)A、C僅與礦山特征相關，而與爆破震動特征相關性較低，所以在計算非剪切作用為主導的震動事件時體現(xiàn)出較低的適用性。

（2）公式（ⅲ）中的礦山特征參數(shù)A、C源自于Kaiser對Brunswick礦山、El Teniente礦山和Creighton礦山的地震數(shù)據(jù)進行95 %置信區(qū)間下回歸分析所得的結果，但理論上金屬礦山間一般都存在較大的差異性，因此公式（ⅲ）中的特征參數(shù)缺乏普適性。

（3）因為深部金屬礦在地質條件、采場結構和作業(yè)情況等方面都極為復雜，所以通過迭代擬合計算和總結出的爆破相關系數(shù)m1、衰減系數(shù)m2及礦山特征參數(shù)C、A可能與礦山實際情況不符，從而導致公式（?。┘肮剑áⅲ┑挠嬎憬Y果與實測值間存在一定差異。

3 基于監(jiān)測數(shù)據(jù)與預測結果的支護建議

根據(jù)金青頂?shù)V區(qū)的震動監(jiān)測數(shù)據(jù)和分析結果，可以發(fā)現(xiàn)不同的震源參數(shù)會對采場及巷道圍巖內(nèi)的PGV分布產(chǎn)生明顯的影響，且在同一震源條件下，PGV的分布在空間上表現(xiàn)出較強的不均勻性。因此，若在支護工作中不考慮震源參數(shù)對PGV分布規(guī)律的影響，則易導致出現(xiàn)支護手段與圍巖穩(wěn)定狀態(tài)不匹配的現(xiàn)象[19-20]。根據(jù)Ju Ma等[12]提出的在考慮震源機制下的圍巖支護設計建議，應基于該地區(qū)震動數(shù)據(jù)監(jiān)測記錄，對地震事件進行震源機制的反演分析，獲得斷層面解和矩震級的分布規(guī)律，根據(jù)所獲取的斷層面解和矩震級分布，合成該地區(qū)在可能的震源模型下及可能的最大矩震級風險下的PGV分布圖，結合基于動態(tài)荷載和巖體質量等級的支護建議（見表3），采取適宜的支護手段。

由于公式（?。┰谠摰V區(qū)條件下的適用性優(yōu)于另外2種，所以選用Sadov公式的變形公式對矩震級Mw=1.65和Mw=0.99的2次震動事件中監(jiān)測所得數(shù)據(jù)進行2種震源假設位置下的PGV分布預測分析。

3.1 震源位于采場M8-3～M8-5區(qū)域

震源位于采場M8-3～M8-5區(qū)域內(nèi)時所獲得的E8、E9、E10和E11穿脈內(nèi)圍巖PGV預測結果分布情況見圖6。

1）由圖6-a）中的預測結果可知：在Mw=1.65震源條件下，E8、E9、E10和E11穿脈內(nèi)圍巖PGV預測最大值分別為0.004 46 m/s、0.004 22 m/s、0.002 89 m/s和0.002 47 m/s，根據(jù)表3對各穿脈內(nèi)的支護作業(yè)給出如下建議：

（1）E8、E9穿脈內(nèi)預測結果較大，建議對2條穿脈靠近采場位置的圍巖進行噴射混凝土+普通錨桿支護，對存在節(jié)理裂隙的區(qū)域可利用錨網(wǎng)進行補充。

（2）E10、E11穿脈內(nèi)圍巖PGV預測值偏低，故整體上采用噴射混凝土作為支護手段。但是，E10、E11穿脈與采場的交叉區(qū)域內(nèi)巖石暴露面積較大，應對該區(qū)域的圍巖進行噴射混凝土+吸能螺栓+錨網(wǎng)的預支護。

2）由圖6-b）中的預測結果可知：在Mw=0.99震源條件下，E8、E9、E10和E11穿脈內(nèi)圍巖PGV預測最大值分別為0.001 48 m/s、0.000 96 m/s、0.000 44 m/s 和0.000 19 m/s。通過與圖6-a）所示結果比較可知：

（1）Mw=0.99震源條件下PGV隨震源距離的變化趨勢與Mw=1.65大致相同，但臨近穿脈或采場間的PGV分布差異性較小，整體上PGV分布會更加均勻。

（2）根據(jù)上述現(xiàn)象，可將噴射混凝土作為各穿脈的通用支護手段，同時針對作業(yè)點附近的圍巖或存在局部破碎的圍巖，可考慮在噴射混凝土的基礎上安裝普通錨桿作為補充支護。

3.2 震源位于E10-1鉆孔硐室內(nèi)

震源位于E10-1鉆孔硐室內(nèi)時所獲得的E8、E9、E10和E11穿脈內(nèi)圍巖PGV預測結果分布情況見圖7。

1）由圖7-a）中的預測結果可知：在該震源條件下，E8、E9、E10和E11穿脈內(nèi)圍巖PGV預測最大值分別為0.001 84 m/s、0.002 71 m/s、0.004 49 m/s和0.003 22 m/s，根據(jù)表3對各穿脈內(nèi)的支護作業(yè)給出如下建議：

（1）E10穿脈內(nèi)預測結果最大，且由于穿脈與采場及鉆孔硐室相接，圍巖的暴露面積大。故建議對E10穿脈靠近采場位置的圍巖進行噴射混凝土+普通錨桿+錨網(wǎng)支護，對靠近主巷道的圍巖進行噴射混凝土支護。

（2）E11穿脈內(nèi)圍巖PGV預測值偏低，建議對E11穿脈采用噴射混凝土的方式作為支護手段。由圖7可知，E11、E10穿脈間采場內(nèi)的預測值偏大，且該采場內(nèi)頂板圍巖破碎并已設置木樁支護，所以建議在木樁支護的基礎上添加噴射混凝土+普通錨桿作為補充支護。

（3）E8、E9穿脈內(nèi)PGV預測值普遍低于0.002 5 m/s，故對E8、E9穿脈圍巖采取噴射混凝土支護即可。

2）由圖7-b）中的預測結果可知：在Mw=0.99震源條件下，E8、E9、E10和E11穿脈內(nèi)圍巖PGV預測最大值分別為0.000 12 m/s、0.000 67 m/s、0.001 50 m/s 和0.001 48 m/s，根據(jù)表3對各穿脈內(nèi)的支護作業(yè)給出如下建議：

（1）圖7-b）中各穿脈內(nèi)PGV值都較小，說明該情況下震動對圍巖穩(wěn)定性的影響較弱，且各穿脈間PGV分布的差異性較低，所以建議對E8、E9、E10和E11穿脈內(nèi)的圍巖采用噴射混凝土支護。

（2）由于E10穿脈與采場及鉆孔硐室相接，其圍巖的暴露面積大。故建議對E10穿脈靠近采場位置的圍巖進行噴射混凝土+普通錨桿支護。

4 結 論

本文在研究國內(nèi)外震動機理和相關理論的基礎上，采用理論分析、現(xiàn)場監(jiān)測數(shù)據(jù)分析和擬合相關系數(shù)計算等手段對不同震源參數(shù)下PGV分布規(guī)律和預測技術進行了深入研究和分析，得出以下結論：

1）震源距離對PGV分布產(chǎn)生明顯作用;且PGV隨震源距離的衰減程度與震源距離之間存在密切的聯(lián)系，當距震源較遠時，震源距離對PGV的衰減作用逐漸弱化。

2）相同震源距離下的PGV會隨矩震級增大顯著變大，且這種影響會根據(jù)距震源遠近逐漸衰減。

3）不同的斷層滑移角下的PGV大小和分布情況存在較大差異，一般情況下走滑斷層（αr=0°）>正斷層（αr=45°）>垂直斷層（αr=±90°）≈逆斷層（αr=-45°）。

4）Sadov公式的變形公式在爆破條件下PGV計算中體現(xiàn)出良好的適用性，同時也說明運用與震動事件的類型和機理相符的公式，所得結果更能反映PGV的實際分布規(guī)律。

5）通過迭代擬合獲得的McGarr經(jīng)驗公式中的特征值A、C不能完全適合該礦區(qū)，且Kaiser回歸分析獲得的特征值在特定礦山下缺乏適用性，導致計算結果與實測值間的相關性較低。

6）在制定采場或巷道的支護方式時，應結合不同震動事件下，支護區(qū)域及周邊圍巖的PGV分布情況，根據(jù)不同震源機制下的圍巖支護設計建議選擇更加合理和時效化的支護方式。

[參 考 文 獻]

[1] MORRISION R G K.Theory and the practical problem of rock bursts[J].Engineering and Mining Journal，1948，149（3）：66-72.

[2] 馮叔瑜，呂毅，顧毅成.城市控制爆破[M].北京：中國鐵道出版社，1987.

[3] 何姣云.礦山采動災害監(jiān)測及控制技術研究[D].武漢：武漢理工大學，2007.

[4] 熊仁欽.頂板大面積來壓破壞機理的研究[J].煤炭學報，1995（增刊1）：38-41.

[5] 陳壽如，宋光明，史秀志，等.近河堤采礦爆破的震動監(jiān)測與控制[J].中國有色金屬學報，2000，10（1）：136-139.

[6] 韓子榮.金川礦區(qū)露天地下聯(lián)合開采的爆破震動安全性評價[J].礦冶工程，1985，5（1）：6-11.

[7] 王民壽，郭慶海.用雙隨機變量回歸改進爆破震速回歸分析[J].爆炸與沖擊，1998，18（3）：283-288.

[8] 徐全軍，劉強，聶渝軍，等.爆破地震峰值預報神經(jīng)網(wǎng)絡研究[J].爆炸與沖擊，1999，19（2）：133-138.

[9] 徐全軍，張慶明，揮壽榕.爆破地震峰值的神經(jīng)網(wǎng)絡預報模型[J].北京理工大學學報，1998，18（14）：472-475.

[10] 黃光球，桂中岳.確定爆破工程中真實經(jīng)驗公式的遺傳規(guī)劃方法[J].工程爆破，1997，3（3）：15-22.

[11] 熊代余，顧毅成.巖石爆破理論與技術新進展[M].北京：冶金工業(yè)出版社，2002.

[12] MA J，DONG L J，ZHAO G Y，et al.Qualitative method and case study for ground vibration of tunnels induced by faultslip in underground mine[J].Rock Mechanics and Rock Engineering，2018，52（3）：1-15.

[13] 田蜜.爆破擾動誘發(fā)煤礦沖擊地壓的數(shù)值模擬研究[D].阜新：遼寧工程技術大學，2014.

[14] 董隴軍，王鈞暉，馬舉.不同微震震源機制下地下硐室圍巖響應及支護建議[J].隧道與地下工程災害防治，2019，1（3）：68-76.

[15] 肖和平.煤礦構造礦震機理[J].湖南地質，1999，18（2/3）：141-146.

[16] 肖和平.煤礦礦震應力窗口效應[J].華南地震，1999，19（1）：85-90.

[17] 董隴軍，李夕兵，唐禮忠，等.無需預先測速的微震震源定位的數(shù)學形式及震源參數(shù)確定[J].巖石力學與工程學報，2011，30（10）：2 057-2 067.

[18] 鄒景波.爆破地震波作用下結構的動力響應及安全評價研究[D].青島：青島理工大學，2012.

[19] 毛暉.建筑物爆破震動的安全控制技術研究[D].長沙：中南大學，2004.

[20] 許紅濤，盧文波.幾種爆破震動安全判據(jù)[J].爆破，2002，19（1）：8-10.

Study on the distribution and prediction of PGV in surrounding

rock under different seismic source parameters

Wang Xi1，Li Shiqing2，Li Qiutao2，Cheng Li1，Liu Huanxin1

（1.Deep Mining Laboratory Subsidiary of Shandong Gold Mining Technology Co.，Ltd.;

2.Shandong Jinzhou Mining Group Co.，Ltd.）

Abstract：In order to totally understand the law of seismic wave propagation and distribution in the surrounding rock under different seismic source mechanisms，the paper studies the distribution of PGV under different seismic source parameters based on field monitoring data.The results show that the distance between stations and the source，moment magnitude scale，and the rake of fault slip have different effects on the distribution of PGV.The applicability of each prediction methods is analyzed by comparing the fitting correlation coefficient between PGV calculated from empirical formulas and others from field monitoring.The study finds that the predicted value of PGV calculated from the Sadovs prediction formula has a higher correlation to the measured result with blasting vibration event.Based on the above research results，the spatial distribution of PGV under dynamic seismic mechanisms is calculated before suggestions on the surrounding rock support under different conditions of the seismic source location and magnitude are given.

Keywords：metal mine;seismic source;vibration monitoring;PGV;prediction formula;fitting correlation coefficient