特許 7114849 干渉信号での鉱山による地震のＰ波を識別するための方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2022-08-01

(45)【発行日】2022-08-09

(54)【発明の名称】干渉信号での鉱山による地震のＰ波を識別するための方法

(51)【国際特許分類】

   G01V 1/30 20060101AFI20220802BHJP

   G01H 1/00 20060101ALI20220802BHJP

【ＦＩ】

G01V1/30

G01H1/00 E

【請求項の数】 3

(21)【出願番号】P 2020197538

(22)【出願日】2020-11-27

(65)【公開番号】P2022052684

(43)【公開日】2022-04-04

【審査請求日】2020-11-27

(31)【優先権主張番号】202011007422.3

(32)【優先日】2020-09-23

(33)【優先権主張国・地域又は機関】CN

(73)【特許権者】

【識別番号】518431602

【氏名又は名称】中国砿業大学

【氏名又は名称原語表記】ＣＨＩＮＡ ＵＮＩＶＥＲＳＩＴＹ ＯＦ ＭＩＮＩＮＧ ＡＮＤ ＴＥＣＨＮＯＬＯＧＹ

【住所又は居所原語表記】Ｎｏ．１，Ｄａｘｕｅ Ｒｏａｄ，Ｔｏｎｇｓｈａｎ，Ｘｕｚｈｏｕ，Ｊｉａｎｇｓｕ，２２１１１６，Ｃｈｉｎａ

(73)【特許権者】

【識別番号】520468690

【氏名又は名称】華亭煤業集団有限責任公司

【氏名又は名称原語表記】Ｈｕａｔｉｎｇ Ｃｏａｌ Ｉｎｄｕｓｔｒｙ Ｇｒｏｕｐ Ｃｏ．， Ｌｔｄ．

【住所又は居所原語表記】Ｎｏ． １０９， Ｌｉａｎｈｕ Ｒｏａｄ， Ｓｈａｎｇｔｉｎｇ Ｃｏｍｍｕｎｉｔｙ， Ｘｉｈｕａ Ｔｏｗｎ， Ｈｕａｔｉｎｇ Ｃｉｔｙ， Ｐｉｎｇｌｉａｎｇ Ｃｉｔｙ， Ｇａｎｓｕ Ｐｒｏｖｉｎｃｅ，Ｃｈｉｎａ

(73)【特許権者】

【識別番号】520468704

【氏名又は名称】淮陰工学院

【氏名又は名称原語表記】ＨＵＡＩＹＩＮ ＩＮＳＴＩＴＵＴＥ ＯＦ ＴＥＣＨＮＯＬＯＧＹ

【住所又は居所原語表記】Ｍｅｉｃｈｅｎｇ Ｒｏａｄ Ｃａｍｐｕｓ， Ｈｕａｉｙｉｎ Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｔｅｃｈｎｏｌｏｇｙ， Ｎｏ． １， Ｍｅｉｃｈｅｎｇ Ｅａｓｔ Ｒｏａｄ， Ｈｕａｉａｎ Ｃｉｔｙ， Ｊｉａｎｇｓｕ Ｐｒｏｖｉｎｃｅ，Ｃｈｉｎａ

(73)【特許権者】

【識別番号】520468715

【氏名又は名称】徐州弘毅科技発展有限公司

【氏名又は名称原語表記】Ｘｕｚｈｏｕ Ｈｏｎｇｙｉ Ｔｅｃｈｎｏｌｏｇｙ Ｄｅｖｅｌｏｐｍｅｎｔ Ｃｏ．， Ｌｔｄ．

【住所又は居所原語表記】Ｎｏ． ２０－１－６１１， Ｚｈｏｎｇｎｅｎｇ Ｓｃｉｅｎｃｅ ａｎｄ Ｔｅｃｈｎｏｌｏｇｙ Ｐａｒｋ， Ｔｈｉｒｄ Ｒｉｎｇ Ｓｏｕｔｈ Ｒｏａｄ， Ｑｕａｎｓｈａｎ Ｄｉｓｔｒｉｃｔ， Ｘｕｚｈｏｕ Ｃｉｔｙ，Ｊｉａｎｇｓｕ， Ｃｈｉｎａ

(74)【代理人】

【識別番号】100146374

【弁理士】

【氏名又は名称】有馬 百子

(72)【発明者】

【氏名】▲鞏▼ 思園

(72)【発明者】

【氏名】張 汝佩

(72)【発明者】

【氏名】劉 志文

(72)【発明者】

【氏名】瀋 威

(72)【発明者】

【氏名】薛 再君

(72)【発明者】

【氏名】牛 佳勝

(72)【発明者】

【氏名】白 金正

(72)【発明者】

【氏名】陸 強

(72)【発明者】

【氏名】馬 志鋒

【審査官】森口 正治

(56)【参考文献】

【文献】特開２０１０－８５２９４（ＪＰ，Ａ）

【文献】特表２００３－５０４７３７（ＪＰ，Ａ）

【文献】実開昭５２－１４０３７４（ＪＰ，Ｕ）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｇ０１Ｖ １／００－９９／００

Ｇ０１Ｈ １／００－１７／００

(57)【特許請求の範囲】

【請求項1】

微小地震監視システムにより、干渉信号を含む可能性のある鉱山による地震波信号を収集して記憶するステップ（１）であって、
炭鉱坑内のトンネルの床に取り付けられたｎ個の地震波受信チャネルで構成される監視ネットワークにより、特定のサンプリング周波数ＳＦで、切羽に対する作業とトンネルの掘削の過程中における、石炭含有岩内で発生した鉱山による地震によって放出された地震波をリアルタイムで収集し、そして、地上受信ユニットを介して、各受信チャネルからフィードバックされた鉱山による地震波信号をリアルタイムで記憶するステップ（１）と、

鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて、各チャネルのＰ波の到達時刻ベクトルから、対応する最適なＰ波の到達時刻を選択して、初期に識別される全チャネルの最適なＰ波の到達時刻ベクトルＴを構成するステップ（３）であって、

全チャネルの最適なＰ波の到達時刻ベクトルを相関判定式に代入して到達時刻の相関行列Ｒを算出し、任意の２つのチャネルによって識別された信号が同一の鉱山による地震の地震源によって励起された信号であるか否かを判断するために使用するステップ（４）であって、
西安８０座標系では、各地震波受信チャネルの空間座標［ｘ ｙ ｚ］を測定し、同時にトンネルの床における各チャネルのＰ波の速度値Ｖ_Ｐ、既知の震源位置での発破信号に基づいて逆推定し、任意のチャネルｉとチャネルｊの、ステップ（３）で得られた全チャネルの最適なＰ波の到達時刻ベクトルを、下記のチャネル相関算出式に代入して算出を行い、相関行列Ｒを形成し、

Ｒ_ｉｊ＝１である場合、チャネルｉとチャネルｊによって記憶されている鉱山による地震信号は、同一の鉱山による地震の地震源によって励起された信号と決定され、Ｒ_ｉｊ＝０である場合、チャネルｉとチャネルｊの信号が相関していないことを意味し、その結果、得られたデータにより相関行列Ｒを構成するステップ（４）と、
相関行列Ｒに対して最適かつ効果的なチャネルの識別方法により、効果的な鉱山による地震波信号を含むチャネルｋを決定するステップ（５）であって、
最適かつ効果的なチャネルの識別方法は、相関行列Ｒの列ごとに成分の和をとり、

効果的な鉱山による地震波信号チャネルｋは、前記行ベクトルＲ_ｃの成分のうちの最大値ｍａｘ（Ｒ_ｃ）である成分のｊの値のチャネルｋ（ｋは、ｊの値）であると決定し、前記行ベクトルＲ _ｃ の成分のうち前記最大値ｍａｘ（Ｒ _ｃ ）である成分のチャネルが複数ある場合、前記複数のチャネルのうち前記最適なＰ波の到達時刻が最も古いチャネルをチャネルｋと決定するステップ（５）と、
前記チャネルｋを基準として、補正が必要な干渉信号チャネルを判断し、及び、該チャネルの、鉱山による地震波信号の最適なＰ波の到達時刻の有効識別範囲を決定するステップ（６）であって、
前記チャネルｋを基準として、行ベクトルＲ_ｃの成分のうちの前記最大値ｍａｘ（Ｒ_ｃ）未満の成分のｊの値の各チャネルのそれぞれに対して、範囲拡大係数がαである場合における、前記各チャネルの鉱山による地震波信号の、最適なＰ波の到達時刻の有効識別範囲を下記式により算出し、

（なお、上記式において、チャネルｉのｉの値は、前記行ベクトルＲ _ｃ の成分のうちの前記最大値ｍａｘ（Ｒ _ｃ ）未満の成分のｊの値を指し、以下のチャネルｉも同様である）
次に、前記各チャネルについて算出した前記有効識別範囲に基づいて、各チャネルｉの最適なＰ波の到達時刻（Ｔ_Ｐ ^ｉ）_Ｃが前記有効識別範囲内にあるか否かを判断し、到達時刻（Ｔ _Ｐ ^ｉ ） _Ｃ が前記有効識別範囲内にない場合、該チャネルを補正が必要な干渉信号チャネルとして記憶するステップ（６）と、
鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて、前記補正が必要な干渉信号チャネルの、対応するＴ_Ｐが有効識別範囲内にある時の最適なＰ波の到達時刻を識別するステップ（７）であって、
前記チャネルｉのＰ波の到達時刻ベクトルＴ_Ｐ ^ｉにステップ（６）で決定された前記有効識別範囲に含まれる到達時刻がない場合、前記チャネルｉの最適なＰ波の到達時刻（Ｔ_Ｐ ^ｉ）_Ｃは存在せず、該チャネルｉを削除し、前記有効識別範囲に含まれる到達時刻がある場合、鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて、ステップ（３）に従って、前記有効識別範囲内の前記チャネルｉの最適なＰ波の到達時刻（Ｔ_Ｐ ^ｉ）_Ｃを再決定することにより、前記干渉信号チャネルを補正するステップ（７）と、を含むことを特徴とする干渉信号での鉱山による地震のＰ波を識別するための方法。

【請求項2】

前記第１の閾値、第２の閾値、及び第３の閾値は、フロントタイム窓Ｗ^ｉ _Ｂ （ｔ）をｐ個のサンプリングポイントだけシフトバックし、Ｅ_ｉ（ｔ）にヒルベルト変換し、変換された窓内の信号の平均値Ｅ^ｉ _ａｖと標準偏差Ｅ^ｉ _ｓｄを算出し、第１の閾値がＨ^１ _ｉ（ｔ）＝ Ｅ^ｉ _ａｖ＋３Ｅ^ｉ _ｓｄとし、第２の閾値、及び第３の閾値が信号対ノイズ比ＳＮＲの倍数とする方法によって決定される、ことを特徴とする請求項１に記載の干渉信号での鉱山による地震のＰ波を識別するための方法。

【請求項3】

前記地震波受信チャネルは地震波ピックアップであることを特徴とする請求項１に記載の干渉信号での鉱山による地震のＰ波を識別するための方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、干渉信号での鉱山による地震のP波を識別するための方法に関し、炭鉱安全採掘の技術分野に属する。

【背景技術】

【0002】

現在、微小地震監視は、地への衝撃圧力による危険性の分析と早期警告を行うための重要な方法であり、なお、鉱山による地震波信号のうちのP波の識別精度は、鉱山による地震の測位とエネルギー解析の正確さに大きく影響する。識別精度が悪い場合、微小地震監視方法は、分析と早期警告において適切な役割を果たすことができなくなる可能性がある。P波を識別するための従来の方法は、例えば、エネルギー分析法、偏光分析法、高次統計法、フラクタル法、AIC法、スペクトル分析法、又は人工ニューラルネットワーク法などがあり、環境の背景ノイズが弱い場合、P波の到達時刻を正確に識別できるが、しかし、炭鉱坑内の作業環境が騒がしく、掘削工事、車両での輸送、大型電気設備の運転等の様々な干渉源で発生した干渉信号は、上記の識別方法の正確さを深刻に制約しており、その結果、干渉信号を鉱山による地震波信号として誤って識別しまったため、記憶されている鉱山による地震波信号が有効であるか否かを正確に識別できず、結局として、鉱山による地震の測位とエネルギー解析の正確さを確保できない。

【発明の概要】

【発明が解決しようとする課題】

【0003】

本発明は、従来技術における上記不足を鑑み、鉱山による地震波信号における干渉信号を正確に識別できることにより、正確かつ効果的な鉱山による地震波信号を取得し、後段の鉱山による地震の測位とエネルギー解析の正確さを確保する、干渉信号での鉱山による地震のP波を識別するための方法を提供する。

【課題を解決するための手段】

【0004】

上記目的を達成するために、本発明の技術的解決手段として、干渉信号での鉱山による地震のP波を識別するための方法は、
微小地震監視システムにより、干渉信号を含む可能性のある鉱山による地震波信号を収集して記憶するステップ（１）であって、

【0005】

炭鉱坑内のトンネルの床に取り付けられたn個の地震波受信チャネルで構成される監視ネットワークにより、特定のサンプリング周波数SFで、切羽に対する作業とトンネルの掘削の過程中における、石炭含有岩内で発生した鉱山による地震によって放出された地震波をリアルタイムで収集し、そして、地上受信ユニットを介して、各受信チャネルからフィードバックされた鉱山による地震波信号をリアルタイムで記憶し、

【0006】

ここで、掘削工事、採掘車両の通過、スイッチによる励起、作業者の操作などの影響により、これら記憶されている鉱山による地震波形信号の一部は、除去不可能な干渉信号が含まれるステップ（１）と、
移動時間窓と識別判断基準を設定して各チャネルのP波の到達時刻ベクトルT_Pを決定するステップ（２）であって、
いずれかのチャネルをチャネルi(i=1..n)とし、チャネルiに対応して記憶されている鉱山による地震波信号x(t)に対して、フロントタイム窓W_B ⁱ(t)、バックタイム窓W_A ⁱ(t) 、及び時間遅延窓W_D ⁱ(t)という３つの移動時間ステップが1/SFの移動時間窓を一定の長さで設定し、以下のチャネルiの識別判断基準式を決定し、
R_１ ⁱ≧H_１ ⁱ(t)且つR_２ ⁱ(t)≧H_２ ⁱ(t)且つR_３ ⁱ(t)≧H_３ ⁱ(t)
但し、R₁ ⁱ(t)は移動時間窓W_B ⁱ(t)内の信号振幅の絶対値の平均値であり、H_１ ⁱ(t)は第１の閾値であり、R_２ ⁱ(t)はW_A ⁱ(t)内の信号振幅の絶対値の平均値とW_B ⁱ(t)内の信号振幅の絶対値の平均値との比であり、Hⁱ _２(t)は第２の閾値であり、R_３ ⁱ(t)はW_D ⁱ(t)内の信号振幅の絶対値の平均値とW_B ⁱ(t)内の信号振幅の絶対値の平均値との比であり、Hⁱ _３(t)は第３の閾値であり、
上記識別判断基準式によれば、チャネルiのP波の到達時刻ベクトルT_P ⁱを取得し、更に、n個のチャネルに対してそれぞれn個の異なるP波の到達時刻ベクトルT_P ⁱ(i=1..n)を形成できるステップ（２）と、

【0007】

鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて、各チャネルのP波の到達時刻ベクトルから、対応する最適なP波の到達時刻を選択して、初期に識別される全チャネルの最適なP波の到達時刻ベクトルTを構成するステップ（３）であって、
一般に、鉱山による地震波信号の波形のエネルギーが記録されている波形の大部分であるという、鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて、R_２ ⁱ(t)のうちの最大値に対応する時間^t max R_２ ⁱ(t)を、チャネルiの最適なP波の到達時刻を確認する最適な時刻として決定し、つまり、T_P ⁱのうちの^t max R_２ ⁱ(t)に最も近い時刻をチャネルiの最適なP波の到達時刻(T_P ⁱ)_Cとして選択し、n個のチャネルに対して繰り返し処理を掛けることにより、全チャネルの最適なP波の到達時刻ベクトルT=[(T_P ⁱ)_C…(T_P ⁿ)_C]を取得するステップ（３）と、

【0008】

全チャネルの最適なP波の到達時刻ベクトルを相関判定式に代入して、到達時刻の相関行列Rを算出し、任意の２つのチャネルによって識別された信号が同一の鉱山による地震の地震源によって励起された信号であるか否かを判断するために使用するステップ（４）であって、

【0009】

西安８０座標系では、各地震波受信チャネルの空間座標[x y z]を測定し、同時にトンネルの床における各チャネルのP波の速度値V_Pを、既知の震源位置での発破信号に基づいて逆推定し、任意のチャネルiとチャネルjの、ステップ（３）で得られた全チャネルの最適なP波の到達時刻ベクトルを、下記のチャネル相関算出式に代入して算出を行い、相関行列Rを形成し、

【0010】

R_ij=1である場合、チャネルiとチャネルjによって記憶されている鉱山による地震信号は、同一の鉱山による地震の地震源によって励起された信号と決定され、R_ij=0である場合、チャネルiとチャネルjの信号が相関していないことを意味し、その結果、得られたデータにより相関行列Rを構成するステップ（４）と、

【0011】

相関行列Rに対して最適かつ効果的なチャネルの識別方法により、効果的な鉱山による地震波信号（干渉信号が含まれない鉱山による地震波信号）を含むチャネルkを決定するステップ（５）であって、
最適かつ効果的なチャネルの識別方法は、相関行列Rの列ごとに和をとり、

【0012】

効果的な鉱山による地震波信号チャネルkが、行ベクトルR_cのうちの最大値max(R_c)が位置するチャネルであることを決定し、最大値が複数の場合、各チャネルの到達時刻においてチャネルkの最適なP波の到達時刻が最も古いと決定する方法であるステップ（５）と、

【0013】

チャネルｋを基準として、補正が必要な干渉信号チャネルを判断し、及び、該チャネルの、鉱山による地震波信号の最適なP波の到達時刻の有効識別範囲を決定するステップ（６）であって、

【0014】

チャネルｋを基準として、行ベクトルR_cのうちのmax(R_c)未満の各チャネルiのそれぞれに対して、範囲拡大係数がαである場合における、各チャネルの鉱山による地震波信号の、最適なP波の到達時刻の有効識別範囲を決定し、

【0015】

次に、それぞれの有効識別範囲に基づいて、各チャネルiの最適なP波の到達時刻(T_P ⁱ)_Cが該範囲内にあるか否かを判断し、到達時刻が該範囲内にない場合、該チャネルを補正が必要な干渉信号チャネルとして記憶するステップ（６）と、

【0016】

鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて、補正が必要な干渉信号チャネルの、対応するT_Pが有効識別範囲内にある時の最適なP波の到達時刻を識別するステップ（７）であって、

【0017】

T_P ⁱにステップ（６）で決定された有効識別範囲に含まれる到達時刻がない場合、(T_P ⁱ)_Cは存在せず、該チャネルを削除しており、有効識別範囲に含まれる到達時刻がある場合、鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて、ステップ（３）に従って、有効識別範囲内において、チャネルiの最適なP波の到達時刻(T_P ⁱ)_Cを再決定し、干渉信号チャネルを補正するステップを完了させるステップ（７）と、
を含む。
更に、前記の第１の閾値、第２の閾値、及び第３の閾値は、フロントタイム窓W_B ⁱ(t)をｐ個のサンプリングポイントだけシフトバックし、E_i(t)にヒルベルト変換し、変換された窓内の信号の平均値Eⁱ _avと標準偏Eⁱ _sdを算出し、第１の閾値がH¹ _i(t)= Eⁱ _av+3Eⁱ _sdとし、第２の閾値、及び第３の閾値が信号対ノイズ比SNRの倍数とする方法によって決定される。
更に、前記の地震波受信チャネルは地震波ピックアップである。

【発明の効果】

【0018】

本発明は、従来技術と比較すると、まず、既存の微小地震監視システムにより干渉信号を含む可能性のあるn個のチャネルの鉱山による地震波信号を収集して記憶し、長さの異なる３つの移動時間窓を設定し、識別判断基準式を満たすチャネルiのような単一チャネルのP波の到達時刻ベクトルT_P ⁱを算出し、各チャネルに記憶された信号のうち、鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいてT_P ⁱから対応する最適なP波の到達時刻(T_P ⁱ)_Cを選択し、初期に識別される全チャネルのP波の到達時刻ベクトルT=[(T_P ⁱ)_C…(T_P ⁿ)_C]を構成し、Tを相関判定式に代入して到達時刻の相関行列Rを算出することにより２つのチャネルによって識別された信号の発生源が同じであるか否かを判断し、最適かつ効果的なチャネルの識別方法により、到達時刻の相関行列Rに対して、効果的な鉱山による地震波信号の波形を含むチャネルｋを決定し、チャネルｋを基準として補正が必要な干渉信号チャネルを判断し、及び該チャネルの鉱山による地震波信号のP波の到達時刻の有効識別範囲を決定し、鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて補正が必要なチャネルの、対応するT_Pが有効識別範囲内にある時の最適なP波の到達時刻を識別し、その結果、干渉信号を識別して除去し、各チャネルの効果的な鉱山による地震波信号を取得する。したがって、本発明は、鉱山による地震波信号における干渉信号を正確に識別でき、干渉信号によるP波の到達時刻の誤識別という問題を解決し、それによって正確かつ効果的な鉱山による地震波信号を取得し、後段の鉱山による地震の測位とエネルギー解析の正確さを確保する。また、本発明は、実用性が高く、結果の信頼性が高く、適用範囲が広く、且つコンピュータプログラミングが容易にするという利点を有する。

【図面の簡単な説明】

【0019】

【図1】本発明の識別フローチャートである。

【図2】実施例において収集された、干渉信号を含む鉱山による地震波信号の波形図である。

【図3】実施例における、チャネル８のP波の到達時刻ベクトルT_P ⁸を算出する概略図である。

【図4】実施例における、鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいてチャネル８の t_max(R₂ ⁸(t))を決定する概略図である。

【図5】実施例における、鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて、有効識別範囲内でチャネル８のt_max(R₂ ⁸(t))を決定する概略図である。

【発明を実施するための形態】

【0020】

以下、本発明については、更に説明する。
＜実施例＞
図１に示すように、具体的に、
微小地震監視システムにより、干渉信号を含む可能性のある鉱山による地震波信号を収集して記憶するステップ（１）であって、

【0021】

図２に示すように、炭鉱坑内のトンネルの床に取り付けられた８個の地震波受信チャネルで構成される監視ネットワークにより、サンプリング周波数５００Hzで、切羽に対する作業とトンネルの掘削の過程中における、石炭含有岩内で発生した鉱山による地震によって放出された地震波をリアルタイムで収集し、且つ、地上受信ユニットを介して、各受信チャネルからフィードバックされた鉱山による地震波信号をリアルタイムで記憶し、ここで、掘削工事、採掘車両の通過、スイッチによる励起、作業者の操作などの影響により、これら記憶されている鉱山による地震波形信号の一部は、除去不可能な干渉信号が含まれるステップ（１）と、
移動時間窓と識別判断基準を設定して各チャネルのP波の到達時刻ベクトルT_Pを決定するステップ（２）であって、
いずれかのチャネルをチャネルi(i=1..n)とし、チャネルiに対応して記憶されている鉱山による地震波信号x(t)に対して、フロントタイム窓W_B ⁱ(t)、バックタイム窓W_A ⁱ(t) 、及び時間遅延窓W_D ⁱ(t)という３つの移動時間ステップが1/SFの移動時間窓を一定の長さで設定し、以下のチャネルiの識別判断基準式を決定し、
R_１ ⁱ≧H_１ ⁱ(t)且つR_２ ⁱ(t)≧H_２ ⁱ(t)且つR_３ ⁱ(t)≧H_３ ⁱ(t)
但し、R₁ ⁱ(t)は移動時間窓W_B ⁱ(t)内の信号振幅の絶対値の平均値であり、H_１ ⁱ(t)は第１の閾値であり、R_２ ⁱ(t)はW_A ⁱ(t)内の信号振幅の絶対値の平均値とW_B ⁱ(t)内の信号振幅の絶対値の平均値との比であり、Hⁱ _２(t)は第２の閾値であり、R_３ ⁱ(t)はW_D ⁱ(t)内の信号振幅の絶対値の平均値とW_B ⁱ(t)内の信号振幅の絶対値の平均値との比であり、Hⁱ _３(t)は第３の閾値であり、上記の第１の閾値、第２の閾値、及び第３の閾値は、フロントタイム窓W_B ⁱ(t)を２０個のサンプリングポイントだけシフトバックし、第１の閾値、第２の閾値、及び第３の閾値をいずれも定数とし、ここで、定数＝４であるとする方法によって決定され、

【0022】

図３に示すように、上記識別判断基準式によれば、チャネルiのP波の到達時刻ベクトルT_P ⁱを取得し、更に、８個のチャネルのそれぞれに対して、T_P ¹ =[１１.２７４ １１.３２ １１.３２６ １１.３３６]；T_P ² =[１１.２７]； T_P ³ =[１１.３１２ １１.３２６]；T_P ⁴ =[１１.３２２ １１.３６６]；T_P ⁵ =[１.８１６ １.８２ １.８３２ ２.２８８ ３.０１８ ３.０２６ ３.０３ ３.０３４ ３.０３８ ３.０５２ ３.９５２ ３.９７ ３.９７４ ４.８１２]；T_P ⁶ =[１１.３１ １１.３４２]；T_P ⁷ =[１１.３７２ １１.３８２ １１.３９４]；T_P ⁸ =[１０.４３４ １０.４５８ １０.４６６ １１.３７８ １１.３８６]という８個の異なるP波の到達時刻ベクトルを形成することができるステップ（２）と、

【0023】

鉱山による地震波信号のエネルギーが最大であるメカニズムに基づいて、各チャネルのP波の到達時刻ベクトルから、対応する最適なP波の到達時刻を選択して、初期に識別される全チャネルの最適なP波の到達時刻ベクトルTを構成するステップ（３）であって、
鉱山による地震波信号のエネルギーが最大であるメカニズムとは、一般に、鉱山による地震の波形のエネルギーが、記憶されている波形の大部分であることであり、このメカニズムによれば、R_２ ⁱ(t)のうちの最大値に対応する時間^tmaxR_２ ⁱ(t)を、チャネルiの最適なP波の到達時刻を確認する最適な時刻として決定し、つまり、T_P ⁱのうちの^tmaxR_２ ⁱ(t)に最も近い時刻をチャネルiの最適なP波の到達時刻(T_P ⁱ)_Cとして選択し、ここで、(T_P ⁱ)_Cはそれぞれ[１１.２９８ １１.２９８ １１.３１４ １１.３２６ １.８１４ １１.３１２ １１.３７２ １０.４５８]であり、そして、全チャネルの最適なP波の到達時刻ベクトルT=[１１.３２ １１.２７ １１.３１２ １１.３２２ １.８１６ １１.３１ １１.３７２ １０.４５８]を取得するステップ（３）と、

【0024】

全チャネルの最適なP波の到達時刻ベクトルを相関判定式に代入して到達時刻の相関行列Rを算出し、任意の２つのチャネルによって識別された信号が同一の鉱山による地震の地震源によって励起された信号であるか否かを判断するために使用するステップ（４）であって、

【0025】

西安８０座標系では、算出の便宜上、現在の鉱山掘削工事計画の左下隅にある境界点を原点として再決定し、座標系を並進移動し、各地震波受信チャネルの空間座標[x y z]を測定し、同時にトンネルの床における各チャネルのP波の速度値V_Pを、既知の震源位置での発破信号に基づいて逆推定し、具体的には、以下の表のとおりであり、

【0026】

任意のチャネルiとチャネルjの、ステップ（３）で得られた全チャネルの最適なP波の到達時刻ベクトルを、下記のチャネル相関算出式に代入して算出を行い、相関行列Rを形成し、

R_ij=1である場合、チャネルiとチャネルjによって記憶されている鉱山による地震信号は、同一の鉱山による地震の地震源によって励起された信号と決定され、R_ij=0である場合、チャネルiとチャネルjの信号が相関していないことを意味し、その結果、得られたデータにより相関行列Rを構成するステップ（４）と、

【0027】

最適かつ効果的なチャネルの識別方法により、相関行列Rに対して、効果的な鉱山による地震波信号（干渉信号が含まれない鉱山による地震波信号）を含むチャネルを決定するステップ（５）であって、

【0028】

最適かつ効果的なチャネルの識別方法は、具体的に、相関行列Rの列ごとに和をとり、行ベクトルR_C= [５ ４ ６ ５ １ ６ ６ １]を取得し、効果的な鉱山による地震波信号チャネルkが、行ベクトルR_Cのうちの最大値max(R_C)が位置するチャネルであり、つまり、チャネル３、６、７は和が６であることを確認し、３つのチャネルのうちのチャネル６の最適なP波の到達時刻が最も古いため、最適かつ効果的なチャネルがチャネル６であると決定する方法ステップ（５）と、

【0029】

チャネル６を基準として、補正が必要な干渉信号チャネルを判断し、及び、該チャネルの、鉱山による地震波信号の最適なP波の到達時刻の有効識別範囲を決定するステップ（６）であって、

【0030】

チャネル６を基準として、行ベクトルR_Cのうちのmax(R_C)未満の各チャネルに、チャネル１、２、４、５、８の、鉱山による地震波信号は範囲拡大係数がα=1.05である時の最適なP波の到達時刻の有効識別範囲がそれぞれ[１０.６９,１１.９６]、[１０.７２６,１１.９２５]、[１０.７２,１１.９４]、[１０.７５,１１.９]、及び[１０.７１,１１.９４]であることを決定し、
次に、上記有効識別範囲に基づいて、対応する(T¹ _P)_C=１１.３２、(T² _P)_C =１１.２７、(T⁴ _P)_C =１１.３２２、(T⁵ _P)_C =１.８１６、(T⁸ _P)_C =１０.４５８がその範囲内にあるか否かを判断し、その結果、チャネル５及び８が該範囲内にないため、補正が必要な干渉信号チャネルとして記憶されるステップ（６）と、

【0031】

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

【図1】

【図2】

【図3】

【図4】

【図5】