JPH10508109A

JPH10508109A - 地震モニタリング

Info

Publication number: JPH10508109A
Application number: JP9504834A
Authority: JP
Inventors: ヨヘンツシャオ、
Original assignee: ゲオフォルシュンクスツェントルム・ポツダム
Priority date: 1995-07-05
Filing date: 1996-07-05
Publication date: 1998-08-04
Also published as: US5890094A; WO1997002499A1; DE19680510B4

Abstract

(57)【要約】ある調査領域の所定の調査ポイントにおける地震活動のその空間的および時間的における変化をモニタリングする方法であって、特に地震予測の改善のための方法においては、比較的に低い強さ(マイクロ地震)の振動のロケーションとそのタイムが前記領域において時間的・空間的に感知され、そしてこのマイクロ地震の地震活動は、当該調査ポイントおよび調査タイムにおいて少なくとも１つのパラメータ(ＳＥＩＳＭＯＬＡＰ)手段によって数量化され、１つの確率値が発生パラメータとして確定され、タイムにおける確率値の発生形成がＳＥＩＳＭＯＬＡＰの変則的異常を検出するためにモニタされ、その活動のクラスタおよび／または相対的な地震の不活動フェーズによって、誘導される更に激しい地震が発生するロケーションを正確に高い信頼性で提供する。

Description

【発明の詳細な説明】地震モニタリング記本発明は、別紙の請求項１前段に記載のようなある１つの調査地域における地震活動パラメータのモニタリングのための一方法に係わり、特に地震の予報( 予測)の可能性を改善するものに関する。今日までのところ、予想される地震の発生時刻、場所(震央)および強さ(マグニチュード)の確かな精度での予報は不可能であり、よって、対応する安全な先行的警告、詳しくはその人口に留意した保護のための警告がとられ得る。従来知られている方法は、地震警報をこれに一致してセッティングするに必要な信頼性を発揮することはできなかった。また、例えば大規模な避難を実施することの判断をするには信頼性が乏しかった。ＳＥＩＳＭＯＬＡＰ(ＳＥＩＳＭｉｃＯｖｅｒＬＡＰｐｉｎｇ)という名称で本出願人により提案されている方法は、地震予測においてかなり改良された一方法である。この方法は、″GeoForschungszentrum Potsdam″の１９９２/９３の２年レポートに記述されている。そこでこの文献の概要についてここに示す。この公知のＳＥＩＳＭＯＬＡＰ方法は、２つの現象因子、つまり名称的には地震のクラスタリングおよび地震の不活動(不活性)と呼べるものの記録および数量化から構成されたものである。地震クラスタリングは、その後に続く激しい地震の揺れを、時間的および空間的にもマイナーなマイクロ地震について注目し予測することにより解る。この現象は、いわゆるＳＥＩＳＭＯＬＡＰパラメータを形成することにより数量化される。このＳＥＩＳＭＯＬＡＰパラメータ(以下に記すＳ)は、次の様に確定される (図２参照)。調査されたある地域においては、観測ロケーションまたはポイント(地点)の１つのネットワークが提供され、発生したマイクロ地震のテーブル(表)またはカタログ(目録)が編成されており、各地震のタイム(時間・時刻)のポイント(時点)や、そのロケーション(場所・位置)およびそのマグニチュードが一覧表にされている。公知の方法により、各マイクロ地震は１つの２次元〜４次元のボディ(体系) が割り当てられている。最も単純なケースにおいてこれは１つの二乗により形成され、これは、震央における１つの三乗によってその地震の深さを含む場合であり、更に１つの４次元構成により、観測報告の発生における時間的インターバル（タイムウインドウ）を考慮することも含んでいた。後者は、次元「ｋｍ³・日」を有し、時間軸が未来ではなく過去だけに向かっている。調査された地域においては、グリッドライン(碁盤の目状)のネットワークが形成され、この中間点には２次元〜４次元ボディが割り当てられ、例えば２次元の場合は、それがこのマイクロ地震に割り当てられた二乗と同じ大きさを有す１つの二乗によって形成されている。通常、これら地震または中間点に割り当てられた２次元〜４次元ボディは、大きさにおいて等しく、その固有の条件、更に詳しくは探求され又は予測される主要な地震の強さに基づいて選択される。パラメータＳ₁は、２次元〜４次元ボディのすべての重複するものの合計から得られる結果であり、それらは夫々、式（１）に基づいて地震または中間点（観測ポイント）に割り当てられる。この２次元のケースは図２に示されている。地震のロケーションはアスタリスクでシンボライズされている。また、黒い領域は重複、即ちパラメータＳ₁を生成する合計を示している。 (Ｘ_j，Ｙ_j，Ｚ_j)は、時刻Ｔ_jにこのポイントで発生した地震の座標であり、同時に、(Ｘ，Ｙ，Ｚ)は調査ポイントの座標である。Ｄ_1-3 ＝ＤＸは、スペースウインドウの次元であり、Ｄ₄ ＝ＤＴはタイムウインドウを表わす。これらのパラメータにより２次元〜４次元ボディの「ボリューム(量)」が定義され、一方、ノーマライジングファクタとして式（１）に表現される。このように、オリジナルのバージョンにおいては、ＳＥＩＳＭＯＬＡＰ手法は、マイクロ地震の空間的な集中の一手法を単に確立し、過去に知られたバージョンにおいても、一時的な周期がＳＥＩＳＭＯＬＡＰパラメータに含まれていた。この目的のために４次元の構成配置は、３つの空間的座標から離れて第４の次元としてのタイム(時間時刻)を含んでスーパーインポーズ（重ね合せ）されていた。これは実際に、マイクロ地震が地球の地殻における被調査ポイントから遠く離れていることを意味するかまたは、時間的に遥かに逆かのぼったことを意味し、僅かに寄与するだけか又は、調査されたロケーションおよび時間的所定ポイントのためのＳＥＩＳＭＯＬＡＰパラメータに対してはまったく寄与しない。対比してみれば、例えば調査ロケーションについての１つの空間的ウインドウ内において、これらはこのポイント付近に位置されており、そして時間的にも遥かに逆かのぼった日付けではない。すなわち、考慮された時間内におけるポイントについての例えばタイムウインドウ内において主要な貢献(寄与)を果たす。よって、調査された地殻ロケーションについてのマイクロ地震の集中が密になればなるほど、これらが更に頻繁に起こり、より大きくなるのがこのＳＥＩＳＭＯＬＡＰパラメータである。このＳＥＩＳＭＯＬＡＰパラメータは、調査の特定なポイントや広範な領域のためにも決定され得る。また、マイクロ地震の強さは任意に重み付けされてもよい。マイクロ地震が強くなればなるほど、選ばれた対応する４次元的な配置構成が更に拡大する。この公知な方法の効果は、サパンカ(Sapanca)地震（1988年12月9日、オープンリヒタスケールにおけるＭ＝４.２）で、アプリケーション（「ポストキャスティング」）によってデモンストレートされた。この地震の２〜３日前には、震央におけるＳＥＩＳＭＯＬＡＰパラメータはゼロ基線から力強く増加し、そして、地震の起こる日の夕刻には特に強くその予兆を知らせている（参照：、図１０の前述した1992/93の２年レポート）。このＳＥＩＳＭＯＬＡＰパラメータの主要な地震の揺れの直接的なタイムウインドウにおける変化の動きは、ＳＥＩＳＭＯＬＡＰの偏差(変則性)に収束する。ＳＥＩＳＭＯＬＡＰ＿１は地震のクラスタリングの直接的な尺度であり、この地震の不活動さは、このパラメータ(ＳＥＩＳＭＯＬＡＰ＿２又は、Ｓ₂)の逆数によって計量される。次式（２）を参照。ＳＥＩＳＭＯＬＡＰ＿２＝Ｓ₂(Ｘ，Ｙ，Ｚ，Ｔ) ＝ (Ｓ₁)^-1 （２）このようにして、活動フェーズのみならず不活動フェーズであることを判別できる可能性が、ＳＥＩＳＭＯＬＡＰ＿１の逆数値を考慮することによってわかっていた。この式（２）に係わり、ＳＥＩＳＭＯＬＡＰ＿１の逆数値は、調査ポイントについての空間における更に大きなウインドウに対応して形成されてもよい。これは、パラメータＳ₁のインデックスＧＲによって示される。Ｓ₁およびＳ₂の両パラメータのリニアな組合せは、Ａ・Ｓ₁ − Ｂ・１／Ｓ_1,GR ＝Ａ・Ｓ₁ − Ｂ・Ｓ_2,GR であり、これは、活動および不活動のこれら両フェーズの逐次的な判別を制限し、インデックスＧＲはこの調査ポイントについての更に広いエリアを示し、Ｓ_1, _GR はこのエリアのために確定されたＳＥＩＳＭＯＬＡＰ＿１のパラメータである。 1988年12月のサパンカ地震に例えばこれを適用すると、明白な不活動フェーズは、地震が起こる活動フェーズのおよそ２週間から１〜２日以前から続いていた。他のケースでは地震のやや長い不活動フェーズが具現され、これにより不活動フェーズは部分的に主要な地震発生の数ケ月から２〜３年以前には既に始まっていた。この遅延は、その後に到来する激しい主要な地震に関係するものである。ＳＥＩＳＭＯＬＡＰ＿２のタイムデベロップメントにおいては、このような不活動フェーズは、特に良好にハイライト(強調)され得る。このような観測報告は、ＳＥＩＳＭＯＬＡＰ＿２パラメータの計算および記録に注目された調査の結果得られたものであり、これは、地震発生の２〜３週間または数ケ月前または、極めて激しい地震の場合では、何年かそれ以前に予測するためのポジションにある。ここに採用された方法は、Ｓ₁またはＳ₂の値が計算やモニタされたロケーションを考慮に入れないで引用されたものであり、異なる基本的な活動形態を呈してもよい。地震の基本的活動が高いあるロケーションにおいては、低い基本的な不活動を呈するロケーションで記録された相対的なものと比べても、どれとも異なる意味を有している。従来この方法はこれらロケーションの相違を考慮することができないので、誤った予測が行われてしまいう場合もあった。それ故に、あらゆるケースにおいて予測が信頼できるものではなかった。本発明の目的は、調査領域における地震活動のパラメータをモニタリングするための改善された方法および、この方法を実施するための装置を提供することにあり、この方法は特に地震予測における信頼性の増大を図るものである。さらに詳しくは、ロケーション的に記録されたＳＥＩＳＭＯＬＡＰ＿２のパラメータを比較する改良された手段を備えることのできる地震活動パラメータをモニタする一方法を定義することが本発明の目的である。この目的は請求項１および請求項１１の要旨によって達成される。また本発明の好適な実施例はこれらの従属項に提示されている。本発明によれば上記目的が達成でき、詳しくは、ＳＥＩＳＭＯＬＡＰ＿１またはＳＥＩＳＭＯＬＡＰ＿２の値が、切迫した上述の主要な地震発生（地震の強さＭ＞４）の１つのインジケータとしてはもはや調査されない。しかしながら、その代わりとしてＳＥＩＳＭＯＬＡＰの値の妥当性は使用される。ＳＥＩＳＭＯＬＡＰ＿２の極大値は、例えば地震の不活動ケースにおいて観察され、そして、極めて低い確率を伴う警告可能な値として発生する。これは、マイクロ地震の存在する場合にはロケーション的な値は比較しても等しく成り得ることを確実にする。この確率は、例えば「パーソン(Pearson)タイプ３の確率分布」をＳＥＩＳＭＯＬＡＰ値の頻度分布に適応するような方法によって定義され得る。この確率分布は、極値である水位のモデリングを行う水管理にたびたび利用されるものである。この方法は、それぞれのロケーションのためや１つの予測が行われるため、適応可能な頻度分布が１つ存在する。即ち、時間的にも十分に長い期間にわたって記録されたＳＥＩＳＭＯＬＡＰの値のロケーションのために、このＳＥＩＳＭＯＬＡＰの値が供給され、そしてこれらＳＥＩＳＭＯＬＡＰの値はパーソンタイプ３確率分布の１つのプロット(点)によって適合され得る。本発明は、次に添付された図面を参照することにより更に詳しく記述される。すなわち、図１は、調査の２つの異なる地点のためのＳＥＩＳＭＯＬＡＰ＿１パラメータの確率分布をプロットしたものである。図２は、従来技術から知られているものとして、ＳＥＩＳＭＯＬＡＰ＿１パラメータを示したものである。本発明によれば、パラメータＳＥＩＳＭＯＬＡＰ＿１又はＳＥＩＳＭＯＬＡＰ＿２の少なくとも１つのパラメータが、調査領域、特に興味深いロケーションにおいては所定時間を超過して最初にまず確定され、その後、対応する頻度分布がプロットされる。この目的のためには、ＳＥＩＳＭＯＬＡＰの値の計算に使用されるタイムウインドウのｎ回（ｎ＞１）の期間が選択される。通常このｎは３に等しく、例えば１００日の１つのタイムウインドウのためには、３００日のＳＥＩＳＭＯＬＡＰのデータがその頻度分布のプロットやパーソンタイプ３の確率曲線の適応のために必要となる。つづいて、その対応する発生のためには、１つの確率分布をその頻度分布に適合することで１つの確率値が各パラメータから確定される。好ましくは、強い非対称を示す確率分布が選択され、相対的に希な発生の極値を適用することを許すものが選択される。これは、例えばパーソンタイプ３確率分布によって満足されるものである。次式(3)によれば、Ｐ(ε) ＝（１/Γ(ｍ)）∫₀ ^εα^m-1ｅ^- ^αｄα （３） ε ＝（２/Ｃss）｛(S−μs)/σs＋２/Ｃss｝＞０ｍ＝（４/Ｃ²ss）＞０ μs ＝ (１/ι)Σ_i=0 ^i-1Ｓi σs＝（１/(ι−1)）｛Σ_i=0 ^i-1（Ｓi-μs）²｝^1/2 Ｃss ＝｛ι^1/2Σ_i-0 ^i-1（Ｓi-μs）³｝/｛（ι-１）^3/2σs³｝ここで、ｓ＝Ｓ₂(Ｔ_k)は、Ｔ＝Ｔ_kの時点でのＳＥＩＳＭＯＬＡＰ＿２の即値であり、確率値に対してこの値が割り当てられる。ｓ_i＝Ｓ₂(Ｔ_k)は、固定的なタイムウインドウにおけるＴ＝Ｔ_k-iの時点でのＳＥＩＳＭＯＬＡＰ＿２の値であり、１はこのウインドウにおけるＳＥＩＳＭＯＬＡＰ＿２の値の番号であり、Γ(ｍ)はガンマ関数である。この式（３）において、Ｐ(ε)は、超過した値εの確率である。もし、Ｃssが正数であれば、Ｐ(ε)は超過したＳＥＩＳＭＯＬＡＰ＿２の値ｓに等しい。この場合、１/Ｐ(ε)（１/確率に対応する式）は、地震の不活動の直接的な基準尺度である。よって、更に高いものは１/Ｐ(ε)であり、それより低いものは極めて低レベルの地震活動で平常な状態であるとして考慮され得る。一方、もしＣssが負数ならば、１/(１−Ｐ(ε))がその代わりに使われる。それぞれの確率値から、更なるパラメータＳＥＩＳＭＯＬＡＰ＿３が確定される。すなわち、１/Ｐ(ε) ＳＥＩＳＭＯＬＡＰ＿３＝Ｓ₃(Ｘ，Ｙ，Ｚ，Ｔ)＝（４）１/(１−Ｐ(ε)) 少なくとも１つのパラメータ(ＳＥＩＳＭＯＬＡＰ＿１により定義され得る確率値ＳＥＩＳＭＯＬＡＰ＿２として相対的に記述されるものとしての)この確率値は次に、興味深い調査ポイント(Ｘ，Ｙ，Ｚ)のため、調査タイム(Ｔ_k)毎に格納される。図１において、２つの異なるロケーションのパラメータＳＥＩＳＭＯＬＡＰ＿１の確率分布は、例として曲線Ａおよび曲線Ｂ(破線)によって表わされている。ＳＥＩＳＭＯＬＡＰ＿２は、ＳＥＩＳＭＯＬＡＰ＿１の逆数値を主に表わすので、対応する確率曲線はＳＥＩＳＭＯＬＡＰ＿２のためにもプロットされ得る。また、図１から次のことが明らかである。すなわち、プロットされたＡに代表されるロケーションは、比較的低いレベルの地震活動を示す。何故ならばこれは、ＳＥＩＳＭＯＬＡＰ＿１の比較的低いレベルにあるプロットされたＢに比べて、その確率分布の最大を示す故である。したがって、相対的な地震の不活動のフェーズはロケーションＡにおいては更に大きな確率にあると考えられ、ＳＥＩＳＭＯＬＡＰ＿１の低い値からは、まだ差し迫った活動を結論づけるような判断は不要である。ここで、プロットＢの場合においては、確率分布の最大値はＳＥＩＳＭＯＬＡＰ＿１のより高い値にシフトされる。そしてロケーションＢは、地震活動の相対的な基準レベルを示す。また、ロケーションＢにおいては、地震活動の相対的レベルのフェーズが、ロケーションＡのそれよりもかなり大きい確率にあり、よってここでも、ＳＥＩＳＭＯＬＡＰ＿１の高い値から、まだ差し迫った活動を結論づけるような判断は不要である。この様にこの手法は、興味深いそれぞれのロケーションのための確立されたＳＥＩＳＭＯＬＡＰ値の時間における開発技術ではないが、その確率分布から確定され得るそれぞれの確率ではある。この確率値はロケーション的に比較可能であり、この予測方法の比較可能性はすべてにロケーショナルに改善され、よって、変則性のある空間的ポジションがさらに正確に決定され得る。結論的には、地震を引き起こす活動の発生はそれぞれ記録され、そして／または、調査ポイントの連続に対応する確率値が著しい変化（降下又は上昇傾向）を示す時期においては、マイクロ地震のクラスタおよび／または、調査ポイントの地震的不活動の形態で警告する。この警告は、時間単位毎の定義された変化率のしきい値に対し超過する場合において、その調査領域に対応する確率が所定のしきい値を侵害したら直ちに行われるものであり、これは光学的出力または音響的警報または他の表示警告というような適宜なものであればよい。この発明による方法によれば、この方法は空間的に比較された結果を許容するという利点からは離れ、上述のようにロケーションの更に良好なピンポイント化を有し、また、信号対ノイズ率が改善されることもまたこの利点である。これに加えて、その結果は、選択された空間ウインドウにおいて起こる地震発生の平均回数の多少に係わらず達成される。本発明に基づく方法は、「ポストキャスティング(postcasting)」な１００件以上の地震から得られるデータで継続的にテストされ、被災地である例えば、スピタック(アルメニア、１９８８年)、ローマプリエタ(ＵＳ、１９８９年)、ランダース(ＵＳ、１９９２年)、北海道(日本、１９９３年)、ノースリッジ(米国、１９９４年)、神戸(日本、１９９５年)およびアイゴン(ギリシャ、１９９５年 )を含んでいる。本発明のさらに特殊な改良や変形例としては、更なるウインドウが５次元のために与えられてもよく、これは詳しくはマイクロ地震のマグニチュードであり、その他４つの次元（空間的および時間的）のための方法と同様に応用し得る。この変形例は、地震がたとえ所定のマグニチュード範囲だけで発生した場合であっても、地震を示すＳＥＩＳＭＯＬＡＰの変則的異常を検知することを可能にする。この目的のため、調査ポイントは中心的なマグニチュードに関する１つのマグニチュードウインドウに割り当てられる。同様に、各マイクロ地震はその送られたマグニチュードに係わる１つのウインドウに割り当てられ、通常、同じサイズのウインドウは両種類のために選択される。２つのウインドウの重なりは次に、他の４つの次元と同様な方法で、ＳＥＩＳＭＯＬＡＰパラメータに対して、マグニチュードの所定帯域（又はインターバルまたはレンジ）の寄与分担を決定する。したがって、この処理は、許容できる更に低いマグニチュードのしきい値であっても有利に実行でき、例えば極めて激しい地震の発生を誘導する地震活動（マイクロ地震のマグニチュードでやや大きく現れるもの）は、これら発生からは別個な独立したものであり、（マイクロ地震のマグニチュードでやや小さく単に現れること）はそれほど激しいものではない。また、それぞれの分離・判別はこの方法（タイムウインドウおよび空間ウインドウ）の他の任意な２つのパラメータの適宜な選択によって同様に成されてもよい。本発明に基づく方法を実施するための装置は、例えば、１５個の地震計のネットワークを使用し、これら地震計は調査地域内に適切に配置されたものである。その調査地域は、例えばおよそ２００・１００kmのサイズを有している。これら地震計は地震の振動を感知して、そのマイクロ地震の時間、位置および強さを決定するシステムにこの振動値を供給し、次には前述したような方法を具現化する。このシステムは、本発明に基づく方法の処理ステップのそれぞれを実行するための構成手段としては特にコンピュータまたはメモリユニットであってもよいし、また、適宜な警告信号手段であってもよい。

Claims

【特許請求の範囲】１．一調査地域において、所定の調査ポイント(Ｘ，Ｙ，Ｚ)および、調査タイム(Ｔ_k)での地震活動の空間的および時間的な変化をモニタリングする方法であって、この方法においては、相対的に低い強さ(マイクロ地震)の揺れのロケーションおよびタイムが前記調査地域で感知され、その後このマイクロ地震の地震活動は、前記調査ポイント(Ｘ，Ｙ，Ｚ)および前記調査タイム(Ｔ_k)の少なくとも１つのパラメータ(Ｓ₁，Ｓ₂)の手段によって数量化され、前記パラメータ(Ｓ₁，Ｓ₂) のそれぞれの少なくとも１つは、前記マイクロ地震の空間的および時間的な蓄積のための時間的変化の可能な尺度であり、および／または前記調査ポイント(Ｘ，Ｙ，Ｚ)での当該地震の不活動であって、ａ）所定の期間にわたって少なくとも１つのパラメータ(Ｓ₁，Ｓ₂)の値を感知し、その後、パラメータ(Ｓ₁，Ｓ₂)のそれぞれの少なくとも１つの頻度分布を生成するステップと、ｂ）前記頻度分布に１つの確率分布を適合することにより、各パラメータの発生に対応してそれぞれ１つの確率値を確定するステップと、ｃ）更なる各調査タイム(Ｔ_k)における前記調査ポイント(Ｘ，Ｙ，Ｚ)のための少なくとも１つの前記確率値を確定および格納するステップと、ｄ）調査ポイントの連続に対応する確率値が著しい変化を示す時期において、マイクロ地震のクラスタおよび／または、この調査ポイントの地震的不活動の形態で、地震を引き起こす活動の発生を記録および／または信号を送るステップと、から成ることを特徴とする地震モニタリングの方法。２．少なくとも１つのパラメータは、第１パラメータ(SEISMOLAP＿1，Ｓ₁) および／または第２パラメータ(SEISMOLAP＿2，Ｓ₂)によって形成され、前記第２パラメータ(SEISMOLAP＿2，Ｓ₂)は、前記第１パラメータ(SEIS MOLAP＿1，Ｓ₁)の逆数であることを特徴とする、請求項１に記載の方法。３．少なくとも１つのパラメータは、第１パラメータ(SEISMOLAP＿1，Ｓ₁) および／または第２パラメータ(SEISMOLAP＿2，Ｓ₂)によって形成され、後者が調査地域に比べてより大きい基準地域のために確定された場合においては、前記第２パラメータ(SEISMOLAP＿2，Ｓ₂)は、前記第１パラメータ(SEISMOLAP＿1，Ｓ₁ )に対応するパラメータ(SEISMOLAP＿1_G，Ｓ_1,G)の逆数であることを特徴とする、請求項１に記載の方法。４．前記第１パラメータ(SEISMOLAP＿1，Ｓ₁)は、各マイクロ地震に対して空間的および時間的に１つのウインドウが割り当てられるように決定され、その後、各マイクロ地震のこれらウインドウの重なりは計算されて、前記調査ポイントにおける空間的および時間的に定義された１つのウインドウを伴って合算されることを特徴とする、請求項２又は請求項３に記載の方法。５．この方法の前記ステップａ）における前記の期間は倍数であり、詳しくは、前記のタイムウインドウの長さの３倍であることを特徴とする、請求項４に記載の方法。６．前記確率分布は強い非対称な確率分布であり、これは相対的に希に発生する極値に対して適用可能であることを特徴とする、請求項１に記載の方法。７．前記確率分布は、パーソン(Pearson)タイプ３の確率分布であることを特徴とする、請求項６に記載の方法。８．前記の空間的ウインドウのマグニチュードは、前記マイクロ地震の強度に応答して選択されることを特徴とする、請求項４に記載の方法。９．前記調査ポイントおよび、前記マイクロ地震のそれぞれは、１つのウインドウに付加的に割り当てられ、前記調査ポイントの前記マグニチュードの前記重なりは、前記第１のパラメータの計算に、前記マイクロ地震の各マグニチュードウインドウと共に入力されることを特徴とする、請求項４に記載の方法。１０．前記調査ポイントは、更に低いマグニチュードしきい値が割り当てられ、マイクロ地震では、前記しきい値の以上にあるマグニチュードのものだけが計算されることを特徴とする、請求項１に記載の方法。１１．調査された前記領域における地震活動を検出するための装置の１つのネットワークであり、前記の少なくとも１つのパラメータを計算するためおよび、前記確率分布を前記頻度分布に適合させるための１つのコンピュータシステムによるものであることを特徴とする、前述した何れかの請求項に記載の方法を実施するための装置。