JP2022052684A - 干渉信号での鉱山による地震のｐ波を識別するための方法 - Google Patents

干渉信号での鉱山による地震のｐ波を識別するための方法 Download PDF

Info

Publication number
JP2022052684A
JP2022052684A JP2020197538A JP2020197538A JP2022052684A JP 2022052684 A JP2022052684 A JP 2022052684A JP 2020197538 A JP2020197538 A JP 2020197538A JP 2020197538 A JP2020197538 A JP 2020197538A JP 2022052684 A JP2022052684 A JP 2022052684A
Authority
JP
Japan
Prior art keywords
channel
wave
mine
signal
seismic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2020197538A
Other languages
English (en)
Other versions
JP7114849B2 (ja
Inventor
思園 ▲鞏▼
Siyuan Gong
汝佩 張
Rupei Zhang
志文 劉
Zhiwen Liu
威 瀋
Wei Shen
再君 薛
Zaijun Xue
佳勝 牛
Jiasheng Niu
金正 白
Jinzheng Bai
強 陸
Lu Qiang
志鋒 馬
Zhifeng Ma
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huating coal industry group co Ltd
Xuzhou Hongyi Technology Development Co Ltd
China University of Mining and Technology CUMT
Huaiyin Institute of Technology
Huating Coal Group Co Ltd
Original Assignee
Huating coal industry group co Ltd
Xuzhou Hongyi Technology Development Co Ltd
China University of Mining and Technology CUMT
Huaiyin Institute of Technology
Huating Coal Group Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huating coal industry group co Ltd, Xuzhou Hongyi Technology Development Co Ltd, China University of Mining and Technology CUMT, Huaiyin Institute of Technology, Huating Coal Group Co Ltd filed Critical Huating coal industry group co Ltd
Publication of JP2022052684A publication Critical patent/JP2022052684A/ja
Application granted granted Critical
Publication of JP7114849B2 publication Critical patent/JP7114849B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/01—Measuring or predicting earthquakes
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/288—Event detection in seismic signals, e.g. microseismics
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Acoustics & Sound (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Business, Economics & Management (AREA)
  • Emergency Management (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

【課題】干渉信号を識別し除去して各チャネルの効果的な鉱山による地震波信号を取得する、干渉信号での鉱山による地震のP波を識別するための方法を提供する。【解決手段】まず干渉信号を含む可能性のある鉱山による地震波信号を収集して記憶し、記憶された各チャネルの信号のうちの鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて、対応する最適なP波の到達時刻を選択し、初期に識別される全チャネルのP波の到達時刻ベクトルTを形成し、Tを相関判定式に代入して２つのチャネルによって識別された信号の発生源が同じか否かを判断し、最適かつ効果的なチャネルに基づく識別方法により、効果的な鉱山による地震波信号の波形を含むチャネルｋを決定し、チャネルｋを基準として、補正が必要な干渉信号チャネルを判断し、そして該チャネルの鉱山による地震波信号のP波の到達時刻の有効識別範囲を決定し、最適なP波の到達時刻を識別する。【選択図】図１

Description

本発明は、干渉信号での鉱山による地震のP波を識別するための方法に関し、炭鉱安全採掘の技術分野に属する。
現在、微小地震監視は、地への衝撃圧力による危険性の分析と早期警告を行うための重要な方法であり、なお、鉱山による地震波信号のうちのP波の識別精度は、鉱山による地震の測位とエネルギー解析の正確さに大きく影響する。識別精度が悪い場合、微小地震監視方法は、分析と早期警告において適切な役割を果たすことができなくなる可能性がある。P波を識別するための従来の方法は、例えば、エネルギー分析法、偏光分析法、高次統計法、フラクタル法、AIC法、スペクトル分析法、又は人工ニューラルネットワーク法などがあり、環境の背景ノイズが弱い場合、P波の到達時刻を正確に識別できるが、しかし、炭鉱坑内の作業環境が騒がしく、掘削工事、車両での輸送、大型電気設備の運転等の様々な干渉源で発生した干渉信号は、上記の識別方法の正確さを深刻に制約しており、その結果、干渉信号を鉱山による地震波信号として誤って識別しまったため、記憶されている鉱山による地震波信号が有効であるか否かを正確に識別できず、結局として、鉱山による地震の測位とエネルギー解析の正確さを確保できない。
本発明は、従来技術における上記不足を鑑み、鉱山による地震波信号における干渉信号を正確に識別できることにより、正確かつ効果的な鉱山による地震波信号を取得し、後段の鉱山による地震の測位とエネルギー解析の正確さを確保する、干渉信号での鉱山による地震のP波を識別するための方法を提供する。
上記目的を達成するために、本発明の技術的解決手段として、干渉信号での鉱山による地震のP波を識別するための方法は、
微小地震監視システムにより、干渉信号を含む可能性のある鉱山による地震波信号を収集して記憶するステップ（１）であって、
炭鉱坑内のトンネルの床に取り付けられたn個の地震波受信チャネルで構成される監視ネットワークにより、特定のサンプリング周波数SFで、切羽に対する作業とトンネルの掘削の過程中における、石炭含有岩内で発生した鉱山による地震によって放出された地震波をリアルタイムで収集し、そして、地上受信ユニットを介して、各受信チャネルからフィードバックされた鉱山による地震波信号をリアルタイムで記憶し、
ここで、掘削工事、採掘車両の通過、スイッチによる励起、作業者の操作などの影響により、これら記憶されている鉱山による地震波形信号の一部は、除去不可能な干渉信号が含まれるステップ（１）と、
移動時間窓と識別判断基準を設定して各チャネルのP波の到達時刻ベクトルTPを決定するステップ（２）であって、
いずれかのチャネルをチャネルi(i=1..n)とし、チャネルiに対応して記憶されている鉱山による地震波信号x(t)に対して、フロントタイム窓WB i(t)、バックタイム窓WA i(t) 、及び時間遅延窓WD i(t)という３つの移動時間ステップが1/SFの移動時間窓を一定の長さで設定し、以下のチャネルiの識別判断基準式を決定し、
R１ i≧H１ i(t)且つR２ i(t)≧H２ i(t)且つR３ i(t)≧H３ i(t)
但し、R1 i(t)は移動時間窓WB i(t)内の信号振幅の絶対値の平均値であり、H１ i(t)は第１の閾値であり、R２ i(t)はWA i(t)内の信号振幅の絶対値の平均値とWB i(t)内の信号振幅の絶対値の平均値との比であり、Hi ２(t)は第２の閾値であり、R３ i(t)はWD i(t)内の信号振幅の絶対値の平均値とWB i(t)内の信号振幅の絶対値の平均値との比であり、Hi ３(t)は第３の閾値であり、
上記識別判断基準式によれば、チャネルiのP波の到達時刻ベクトルTP iを取得し、更に、n個のチャネルに対してそれぞれn個の異なるP波の到達時刻ベクトルTP i(i=1..n)を形成できるステップ（２）と、
鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて、各チャネルのP波の到達時刻ベクトルから、対応する最適なP波の到達時刻を選択して、初期に識別される全チャネルの最適なP波の到達時刻ベクトルTを構成するステップ（３）であって、
一般に、鉱山による地震波信号の波形のエネルギーが記録されている波形の大部分であるという、鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて、R２ i(t)のうちの最大値に対応する時間t max R２ i(t)を、チャネルiの最適なP波の到達時刻を確認する最適な時刻として決定し、つまり、TP iのうちのt max R２ i(t)に最も近い時刻をチャネルiの最適なP波の到達時刻(TP i)Cとして選択し、n個のチャネルに対して繰り返し処理を掛けることにより、全チャネルの最適なP波の到達時刻ベクトルT=[(TP i)C…(TP n)C]を取得するステップ（３）と、
全チャネルの最適なP波の到達時刻ベクトルを相関判定式に代入して、到達時刻の相関行列Rを算出し、任意の２つのチャネルによって識別された信号が同一の鉱山による地震の地震源によって励起された信号であるか否かを判断するために使用するステップ（４）であって、
西安８０座標系では、各地震波受信チャネルの空間座標[x y z]を測定し、同時にトンネルの床における各チャネルのP波の速度値VPを、既知の震源位置での発破信号に基づいて逆推定し、任意のチャネルiとチャネルjの、ステップ（３）で得られた全チャネルの最適なP波の到達時刻ベクトルを、下記のチャネル相関算出式に代入して算出を行い、相関行列Rを形成し、
Figure 2022052684000002
Rij=1である場合、チャネルiとチャネルjによって記憶されている鉱山による地震信号は、同一の鉱山による地震の地震源によって励起された信号と決定され、Rij=0である場合、チャネルiとチャネルjの信号が相関していないことを意味し、その結果、得られたデータにより相関行列Rを構成するステップ（４）と、
相関行列Rに対して最適かつ効果的なチャネルの識別方法により、効果的な鉱山による地震波信号（干渉信号が含まれない鉱山による地震波信号）を含むチャネルkを決定するステップ（５）であって、
最適かつ効果的なチャネルの識別方法は、相関行列Rの列ごとに和をとり、
Figure 2022052684000003
効果的な鉱山による地震波信号チャネルkが、行ベクトルRcのうちの最大値max(Rc)が位置するチャネルであることを決定し、最大値が複数の場合、各チャネルの到達時刻においてチャネルkの最適なP波の到達時刻が最も古いと決定する方法であるステップ（５）と、
チャネルｋを基準として、補正が必要な干渉信号チャネルを判断し、及び、該チャネルの、鉱山による地震波信号の最適なP波の到達時刻の有効識別範囲を決定するステップ（６）であって、
チャネルｋを基準として、行ベクトルRcのうちのmax(Rc)未満の各チャネルiのそれぞれに対して、範囲拡大係数がαである場合における、各チャネルの鉱山による地震波信号の、最適なP波の到達時刻の有効識別範囲を決定し、
Figure 2022052684000004
次に、それぞれの有効識別範囲に基づいて、各チャネルiの最適なP波の到達時刻(TP i)Cが該範囲内にあるか否かを判断し、到達時刻が該範囲内にない場合、該チャネルを補正が必要な干渉信号チャネルとして記憶するステップ（６）と、
鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて、補正が必要な干渉信号チャネルの、対応するTPが有効識別範囲内にある時の最適なP波の到達時刻を識別するステップ（７）であって、
TP iにステップ（６）で決定された有効識別範囲に含まれる到達時刻がない場合、(TP i)Cは存在せず、該チャネルを削除しており、有効識別範囲に含まれる到達時刻がある場合、鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて、ステップ（３）に従って、有効識別範囲内において、チャネルiの最適なP波の到達時刻(TP i)Cを再決定し、干渉信号チャネルを補正するステップを完了させるステップ（７）と、
を含む。
更に、前記の第１の閾値、第２の閾値、及び第３の閾値は、フロントタイム窓WB i(t)をｐ個のサンプリングポイントだけシフトバックし、Ei(t)にヒルベルト変換し、変換された窓内の信号の平均値Ei avと標準偏Ei sdを算出し、第１の閾値がH1 i(t)= Ei av+3Ei sdとし、第２の閾値、及び第３の閾値が信号対ノイズ比SNRの倍数とする方法によって決定される。
更に、前記の地震波受信チャネルは地震波ピックアップである。
本発明は、従来技術と比較すると、まず、既存の微小地震監視システムにより干渉信号を含む可能性のあるn個のチャネルの鉱山による地震波信号を収集して記憶し、長さの異なる３つの移動時間窓を設定し、識別判断基準式を満たすチャネルiのような単一チャネルのP波の到達時刻ベクトルTP iを算出し、各チャネルに記憶された信号のうち、鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいてTP iから対応する最適なP波の到達時刻(TP i)Cを選択し、初期に識別される全チャネルのP波の到達時刻ベクトルT=[(TP i)C…(TP n)C]を構成し、Tを相関判定式に代入して到達時刻の相関行列Rを算出することにより２つのチャネルによって識別された信号の発生源が同じであるか否かを判断し、最適かつ効果的なチャネルの識別方法により、到達時刻の相関行列Rに対して、効果的な鉱山による地震波信号の波形を含むチャネルｋを決定し、チャネルｋを基準として補正が必要な干渉信号チャネルを判断し、及び該チャネルの鉱山による地震波信号のP波の到達時刻の有効識別範囲を決定し、鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて補正が必要なチャネルの、対応するTPが有効識別範囲内にある時の最適なP波の到達時刻を識別し、その結果、干渉信号を識別して除去し、各チャネルの効果的な鉱山による地震波信号を取得する。したがって、本発明は、鉱山による地震波信号における干渉信号を正確に識別でき、干渉信号によるP波の到達時刻の誤識別という問題を解決し、それによって正確かつ効果的な鉱山による地震波信号を取得し、後段の鉱山による地震の測位とエネルギー解析の正確さを確保する。また、本発明は、実用性が高く、結果の信頼性が高く、適用範囲が広く、且つコンピュータプログラミングが容易にするという利点を有する。
本発明の識別フローチャートである。 実施例において収集された、干渉信号を含む鉱山による地震波信号の波形図である。 実施例における、チャネル８のP波の到達時刻ベクトルTP 8を算出する概略図である。 実施例における、鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいてチャネル８の tmax(R2 8(t))を決定する概略図である。 実施例における、鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて、有効識別範囲内でチャネル８のtmax(R2 8(t))を決定する概略図である。
以下、本発明については、更に説明する。
＜実施例＞
図１に示すように、具体的に、
微小地震監視システムにより、干渉信号を含む可能性のある鉱山による地震波信号を収集して記憶するステップ（１）であって、
図２に示すように、炭鉱坑内のトンネルの床に取り付けられた８個の地震波受信チャネルで構成される監視ネットワークにより、サンプリング周波数５００Hzで、切羽に対する作業とトンネルの掘削の過程中における、石炭含有岩内で発生した鉱山による地震によって放出された地震波をリアルタイムで収集し、且つ、地上受信ユニットを介して、各受信チャネルからフィードバックされた鉱山による地震波信号をリアルタイムで記憶し、ここで、掘削工事、採掘車両の通過、スイッチによる励起、作業者の操作などの影響により、これら記憶されている鉱山による地震波形信号の一部は、除去不可能な干渉信号が含まれるステップ（１）と、
移動時間窓と識別判断基準を設定して各チャネルのP波の到達時刻ベクトルTPを決定するステップ（２）であって、
いずれかのチャネルをチャネルi(i=1..n)とし、チャネルiに対応して記憶されている鉱山による地震波信号x(t)に対して、フロントタイム窓WB i(t)、バックタイム窓WA i(t) 、及び時間遅延窓WD i(t)という３つの移動時間ステップが1/SFの移動時間窓を一定の長さで設定し、以下のチャネルiの識別判断基準式を決定し、
R１ i≧H１ i(t)且つR２ i(t)≧H２ i(t)且つR３ i(t)≧H３ i(t)
但し、R1 i(t)は移動時間窓WB i(t)内の信号振幅の絶対値の平均値であり、H１ i(t)は第１の閾値であり、R２ i(t)はWA i(t)内の信号振幅の絶対値の平均値とWB i(t)内の信号振幅の絶対値の平均値との比であり、Hi ２(t)は第２の閾値であり、R３ i(t)はWD i(t)内の信号振幅の絶対値の平均値とWB i(t)内の信号振幅の絶対値の平均値との比であり、Hi ３(t)は第３の閾値であり、上記の第１の閾値、第２の閾値、及び第３の閾値は、フロントタイム窓WB i(t)を２０個のサンプリングポイントだけシフトバックし、第１の閾値、第２の閾値、及び第３の閾値をいずれも定数とし、ここで、定数＝４であるとする方法によって決定され、
図３に示すように、上記識別判断基準式によれば、チャネルiのP波の到達時刻ベクトルTP iを取得し、更に、８個のチャネルのそれぞれに対して、TP 1 =[１１.２７４ １１.３２ １１.３２６ １１.３３６]；TP 2 =[１１.２７]； TP 3 =[１１.３１２ １１.３２６]；TP 4 =[１１.３２２ １１.３６６]；TP 5 =[１.８１６ １.８２ １.８３２ ２.２８８ ３.０１８ ３.０２６ ３.０３ ３.０３４ ３.０３８ ３.０５２ ３.９５２ ３.９７ ３.９７４ ４.８１２]；TP 6 =[１１.３１ １１.３４２]；TP 7 =[１１.３７２ １１.３８２ １１.３９４]；TP 8 =[１０.４３４ １０.４５８ １０.４６６ １１.３７８ １１.３８６]という８個の異なるP波の到達時刻ベクトルを形成することができるステップ（２）と、
鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて、各チャネルのP波の到達時刻ベクトルから、対応する最適なP波の到達時刻を選択して、初期に識別される全チャネルの最適なP波の到達時刻ベクトルTを構成するステップ（３）であって、
鉱山による地震波信号のエネルギーが最大であるメカニズムとは、一般に、鉱山による地震の波形のエネルギーが、記憶されている波形の大部分であることであり、このメカニズムによれば、R２ i(t)のうちの最大値に対応する時間tmaxR２ i(t)を、チャネルiの最適なP波の到達時刻を確認する最適な時刻として決定し、つまり、TP iのうちのtmaxR２ i(t)に最も近い時刻をチャネルiの最適なP波の到達時刻(TP i)Cとして選択し、ここで、(TP i)Cはそれぞれ[１１.２９８ １１.２９８ １１.３１４ １１.３２６ １.８１４ １１.３１２ １１.３７２ １０.４５８]であり、そして、全チャネルの最適なP波の到達時刻ベクトルT=[１１.３２ １１.２７ １１.３１２ １１.３２２ １.８１６ １１.３１ １１.３７２ １０.４５８]を取得するステップ（３）と、
全チャネルの最適なP波の到達時刻ベクトルを相関判定式に代入して到達時刻の相関行列Rを算出し、任意の２つのチャネルによって識別された信号が同一の鉱山による地震の地震源によって励起された信号であるか否かを判断するために使用するステップ（４）であって、
西安８０座標系では、算出の便宜上、現在の鉱山掘削工事計画の左下隅にある境界点を原点として再決定し、座標系を並進移動し、各地震波受信チャネルの空間座標[x y z]を測定し、同時にトンネルの床における各チャネルのP波の速度値VPを、既知の震源位置での発破信号に基づいて逆推定し、具体的には、以下の表のとおりであり、
Figure 2022052684000005
任意のチャネルiとチャネルjの、ステップ（３）で得られた全チャネルの最適なP波の到達時刻ベクトルを、下記のチャネル相関算出式に代入して算出を行い、相関行列Rを形成し、
Figure 2022052684000006
Figure 2022052684000007
Rij=1である場合、チャネルiとチャネルjによって記憶されている鉱山による地震信号は、同一の鉱山による地震の地震源によって励起された信号と決定され、Rij=0である場合、チャネルiとチャネルjの信号が相関していないことを意味し、その結果、得られたデータにより相関行列Rを構成するステップ（４）と、
Figure 2022052684000008
最適かつ効果的なチャネルの識別方法により、相関行列Rに対して、効果的な鉱山による地震波信号（干渉信号が含まれない鉱山による地震波信号）を含むチャネルを決定するステップ（５）であって、
最適かつ効果的なチャネルの識別方法は、具体的に、相関行列Rの列ごとに和をとり、行ベクトルRC= [５ ４ ６ ５ １ ６ ６ １]を取得し、効果的な鉱山による地震波信号チャネルkが、行ベクトルRCのうちの最大値max(RC)が位置するチャネルであり、つまり、チャネル３、６、７は和が６であることを確認し、３つのチャネルのうちのチャネル６の最適なP波の到達時刻が最も古いため、最適かつ効果的なチャネルがチャネル６であると決定する方法ステップ（５）と、
チャネル６を基準として、補正が必要な干渉信号チャネルを判断し、及び、該チャネルの、鉱山による地震波信号の最適なP波の到達時刻の有効識別範囲を決定するステップ（６）であって、
チャネル６を基準として、行ベクトルRCのうちのmax(RC)未満の各チャネルに、チャネル１、２、４、５、８の、鉱山による地震波信号は範囲拡大係数がα=1.05である時の最適なP波の到達時刻の有効識別範囲がそれぞれ[１０.６９,１１.９６]、[１０.７２６,１１.９２５]、[１０.７２,１１.９４]、[１０.７５,１１.９]、及び[１０.７１,１１.９４]であることを決定し、
次に、上記有効識別範囲に基づいて、対応する(T1 P)C=１１.３２、(T2 P)C =１１.２７、(T4 P)C =１１.３２２、(T5 P)C =１.８１６、(T8 P)C =１０.４５８がその範囲内にあるか否かを判断し、その結果、チャネル５及び８が該範囲内にないため、補正が必要な干渉信号チャネルとして記憶されるステップ（６）と、

鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて、補正が必要な干渉信号チャネルの、対応するTPが有効識別範囲内にある時の最適なP波の到達時刻を識別するステップ（７）であって、
チャネル５である場合、TP 5に識別範囲[１０.７５,１１.９]に含まれるP波の到達時刻がないため、(TP 5)cは存在せず、そのチャネルで受信した鉱山による地震波信号を削除しており、チャネル８である場合、TP 8に識別範囲[１０.７１,１１.９４]に含まれるP波の到達時刻があるため、識別範囲内においてtmaxR２ i(t)が１１.３５であることを決定でき、これにより、チャネル８のP波の到達時刻ベクトルから、１１.３７８を対応する最適なP波の到達時刻として選択し、その結果、干渉信号チャネルの補正のステップを完了させるステップ（７）と、を含む。
更に、上記地震波受信チャネルは地震波ピックアップである。
上記地震波受信チャネルで構成される監視ネットワーク及び地上受信ユニットは、いずれも炭鉱坑内の既存の設備である。

Claims (3)

  1. 微小地震監視システムにより、干渉信号を含む可能性のある鉱山による地震波信号を収集して記憶するステップ（１）であって、
    炭鉱坑内のトンネルの床に取り付けられたn個の地震波受信チャネルで構成される監視ネットワークにより、特定のサンプリング周波数SFで、切羽に対する作業とトンネルの掘削の過程中における、石炭含有岩内で発生した鉱山による地震によって放出された地震波をリアルタイムで収集し、そして、地上受信ユニットを介して、各受信チャネルからフィードバックされた鉱山による地震波信号をリアルタイムで記憶するステップ（１）と、
    Figure 2022052684000009
    Figure 2022052684000010
    鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて、各チャネルのP波の到達時刻ベクトルから、対応する最適なP波の到達時刻を選択して、初期に識別される全チャネルの最適なP波の到達時刻ベクトルTを構成するステップ（３）であって、
    Figure 2022052684000011
    全チャネルの最適なP波の到達時刻ベクトルを相関判定式に代入して到達時刻の相関行列Rを算出し、任意の２つのチャネルによって識別された信号が同一の鉱山による地震の地震源によって励起された信号であるか否かを判断するために使用するステップ（４）であって、
    西安８０座標系では、各地震波受信チャネルの空間座標[x y z]を測定し、同時にトンネルの床における各チャネルのP波の速度値VP、既知の震源位置での発破信号に基づいて逆推定し、任意のチャネルiとチャネルjの、ステップ（３）で得られた全チャネルの最適なP波の到達時刻ベクトルを、下記のチャネル相関算出式に代入して算出を行い、相関行列Rを形成し、
    Figure 2022052684000012
    Rij=1である場合、チャネルiとチャネルjによって記憶されている鉱山による地震信号は、同一の鉱山による地震の地震源によって励起された信号と決定され、Rij=0である場合、チャネルiとチャネルjの信号が相関していないことを意味し、その結果、得られたデータにより相関行列Rを構成するステップ（４）と、
    相関行列Rに対して最適かつ効果的なチャネルの識別方法により、効果的な鉱山による地震波信号を含むチャネルkを決定するステップ（５）であって、
    最適かつ効果的なチャネルの識別方法は、相関行列Rの列ごとに和をとり、
    Figure 2022052684000013
    効果的な鉱山による地震波信号チャネルkが、行ベクトルRcのうちの最大値max(Rc)が位置するチャネルであることを決定し、最大値が複数の場合、各チャネルの到達時刻においてチャネルkの最適なP波の到達時刻が最も古いと決定する方法であるステップ（５）と、
    チャネルｋを基準として、補正が必要な干渉信号チャネルを判断し、及び、該チャネルの、鉱山による地震波信号の最適なP波の到達時刻の有効識別範囲を決定するステップ（６）であって、
    チャネルｋを基準として、行ベクトルRcのうちのmax(Rc)未満の各チャネルiのそれぞれに対して、範囲拡大係数がαである場合における、各チャネルの鉱山による地震波信号の、最適なP波の到達時刻の有効識別範囲を決定し、
    Figure 2022052684000014
    次に、各有効識別範囲に基づいて、各チャネルiの最適なP波の到達時刻(TP i)Cが該範囲内にあるか否かを判断し、到達時刻が該範囲内にない場合、該チャネルを補正が必要な干渉信号チャネルとして記憶するステップ（６）と、
    鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて、補正が必要な干渉信号チャネルの、対応するTPが有効識別範囲内にある時の最適なP波の到達時刻を識別するステップ（７）であって、
    TP iにステップ（６）で決定された有効識別範囲に含まれる到達時刻がない場合、(TP i)Cは存在せず、該チャネルを削除しており、有効識別範囲に含まれる到達時刻がある場合、鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて、ステップ（３）に従って、有効識別範囲内において、チャネルiの最適なP波の到達時刻(TP i)Cを再決定し、そして、干渉信号チャネルを補正するステップを完了させるステップ（７）と、を含むことを特徴とする干渉信号での鉱山による地震のP波を識別するための方法。
  2. 前記第１の閾値、第２の閾値、及び第３の閾値は、フロントタイム窓Wi B (t)をｐ個のサンプリングポイントだけシフトバックし、Ei(t)にヒルベルト変換し、変換された窓内の信号の平均値Ei avと標準偏差Ei sdを算出し、第１の閾値がH1 i(t)= Ei av+3Ei sdとし、第２の閾値、及び第３の閾値が信号対ノイズ比SNRの倍数とする方法によって決定される、ことを特徴とする請求項１に記載の干渉信号での鉱山による地震のP波を識別するための方法。
  3. 前記地震波受信チャネルは地震波ピックアップであることを特徴とする請求項１に記載の干渉信号での鉱山による地震のP波を識別するための方法。
JP2020197538A 2020-09-23 2020-11-27 干渉信号での鉱山による地震のｐ波を識別するための方法 Active JP7114849B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202011007422.3 2020-09-23
CN202011007422.3A CN112180430B (zh) 2020-09-23 2020-09-23 一种存在干扰信号下的矿震p波初至识别方法

Publications (2)

Publication Number Publication Date
JP2022052684A true JP2022052684A (ja) 2022-04-04
JP7114849B2 JP7114849B2 (ja) 2022-08-09

Family

ID=73957109

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020197538A Active JP7114849B2 (ja) 2020-09-23 2020-11-27 干渉信号での鉱山による地震のｐ波を識別するための方法

Country Status (2)

Country Link
JP (1) JP7114849B2 (ja)
CN (1) CN112180430B (ja)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS52140374U (ja) * 1976-04-19 1977-10-24
JP2003504737A (ja) * 1999-07-01 2003-02-04 コミツサリア タ レネルジー アトミーク モデルによって生成される信号および判別式パラメータに基づくイベントまたは実体または状況の分類のための人工知能のシステム
JP2010085294A (ja) * 2008-10-01 2010-04-15 Tokyo Metropolitan Univ 緊急地震警報を発するための方法および装置

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9103935B2 (en) * 2013-02-04 2015-08-11 King Fahd University Of Petroleum And Minerals Method of first arrival picking of seismic refraction data
CN106772595A (zh) * 2016-12-05 2017-05-31 中国矿业大学 一种用于煤矿震动波波速反演中雷管延时消除的方法
CN107797147B (zh) * 2017-10-25 2019-07-05 中国石油化工股份有限公司 一种地震初至异常值快速剔除方法
CN108535773B (zh) * 2018-07-13 2019-07-05 广东石油化工学院 一种微地震事件检测方法及系统
CN110297271B (zh) * 2019-06-26 2020-09-11 中国矿业大学 一种用于矿震报警的单分量探头p波初至到时修正方法
CN110568499B (zh) * 2019-08-20 2021-06-04 长江大学 一种vsp地震资料的初至波时差校正方法及装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS52140374U (ja) * 1976-04-19 1977-10-24
JP2003504737A (ja) * 1999-07-01 2003-02-04 コミツサリア タ レネルジー アトミーク モデルによって生成される信号および判別式パラメータに基づくイベントまたは実体または状況の分類のための人工知能のシステム
JP2010085294A (ja) * 2008-10-01 2010-04-15 Tokyo Metropolitan Univ 緊急地震警報を発するための方法および装置

Also Published As

Publication number Publication date
JP7114849B2 (ja) 2022-08-09
CN112180430B (zh) 2021-08-20
CN112180430A (zh) 2021-01-05

Similar Documents

Publication Publication Date Title
US7675813B2 (en) Methods and systems for detecting arrivals of interest
CN110988981B (zh) 一种适用于钻爆法隧道的相控阵声波超前预报系统及方法
NO318907B1 (no) Fremgangsmate til a analysere Stoneley-bolgeformer og a karakterisere undergrunnsformasjoner.
CN110297271B (zh) 一种用于矿震报警的单分量探头p波初至到时修正方法
CN107678064B (zh) 一种声波时差实时提取方法
CN106646609B (zh) 多次扫描的微地震多参数联合快速反演方法
US20070081422A1 (en) PP/PS event matching (registration)
Dong et al. Arrival-time detection with multiscale wavelet analysis and source location of acoustic emission in rock
JP2022052684A (ja) 干渉信号での鉱山による地震のｐ波を識別するための方法
Arifullah et al. A preliminary result of automatic earthquakes localization from seismic temporary network in Northern Sumatra
GB2450163A (en) Detecting the location of seismic events without picking events in received seismic wave data
Chen et al. An approach to improve the performance of the earthquake early warning system for the 2018 Hualien earthquake in Taiwan
RU2010131268A (ru) Способ оценки основных характеристик ожидаемого сильного цунамигенного землетрясения и система для его осуществления
US10101485B2 (en) Method of coalescence microseismic mapping including model's uncertainty
CN115793056A (zh) 一种隧道施工超前地质探测方法及系统
CN109521467A (zh) 一种基于突出煤层巷道的超前探测方法
CN110967751A (zh) 基于地面浅井监测的微地震事件的定位方法及存储介质
JP2001083257A (ja) 地震早期検知警報・制御システム
Zollo et al. Source characterization for earthquake early warning
Soliman et al. Evaluation of Saudi National Seismic Network (SNSN) detectability
CN112200069A (zh) 时频域谱减法和经验模态分解联合的隧道滤波方法及系统
CN109507721B (zh) 井中微地震事件质量评价方法及系统
Shahvar et al. Strong motion records in Sarpol-e Zahab earthquake
CN114859411A (zh) 一种微地震震级的确定方法及系统
Manea et al. Evaluation of site parameters to inform seismic site characterization in New Zealand

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20201127

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210128

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220104

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220330

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20220621

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220707

R150 Certificate of patent or registration of utility model

Ref document number: 7114849

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150