(57)【要約】【課題】少数の観測データを用いて震源地決定を行う従来の方法では、（１）精度が低い、（２）解が出るまでに時間がかかる、（３）同時に異なる地点で発生した地震があった時に同一の地震と見なされる、などの問題があった。【解決手段】双曲面Ci1の当該片業が、ｔ＝ｔ1の時にはａ＝０で、Ci1は平面（X＝０）となり、震源が点S1、Siの対称面（X＝０）の点S1側の領域（Di1）にあることが言える。任意の観測点Siに関しても、同じような制限が課されるので、Di1の共通領域たる角柱の内部にあることになる。この場合には深さに関しては、地震活動から知られる制限を使うことになる。Δｔが大きくなると、Ci1はより鋭い曲面となり、点S1に向かって移動すると共に収縮し、震源の存在できる領域Dilは小さくなり、推定値はより精度が高くなる。【選択図】        図３
2006105862.tif 000002
【特許請求の範囲】
【請求項１】
地震波観測データによる震源の時刻および位置の決定方法であって、１個以上の地震波観測データを用い、高度化着未着法によって震源要素および、その精度を推定することを特徴とするリアルタイム地震危険度予測の方法。
【請求項２】
地震波観測データによる震源の時刻および位置の決定方法であって、前記１個以上の地震波観測データを用い、前記高度化着未着法とユレダス法を融合して震源を決定することを特徴とするリアルタイム地震危険度予測の方法。
【請求項３】
地震波観測データによる震源の時刻および位置の決定方法であって、ほぼ同時に２カ所以上で発生した可能性のある地震を、前記高度化着未着法によって識別することを特徴とするリアルタイム地震危険度予測の方法。
【請求項４】
地震波観測データによる震源の時刻および位置の決定方法を使い、特定地区での利用に供するために、高精度で、特定サイトの地震危険度を推定することを特徴とするリアルタイム地震危険度予測の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
地震観測網が一定の水準にあって、データが即時的（リアルタイム）あるいは準リアルタイムで収集されている環境で、地震の発生後かつ、地震波主要動が到達するまでにその情報をユーザーに伝え、地震の早期警報あるいは自動防災対応措置を行うためのリアルタイム地震危険度予測の方法に関する。
【背景技術】
【０００２】
全国的または地域的に整備された、一定規模の地震観測網からのデータに基づく、リアルタイム震源要素決定の方法については、代表的にユレダス法（鉄道総研報告１９巻３号、ｐｐ１〜１２）および、着未着法（特願２００１−３０８２１８号公報参照）が提案され、実証試験も行われている。前者の特徴は、観測点が一点から震源の情報が求まることであり、後者では、２ケ所以上の観測点データから相当に高い精度で震源パラメータが求まることである。
【０００３】
【特許文献１】特願２００１−３０８２１８号公報
【非特許文献１】鉄道総研報告１９巻３号、ｐｐ１〜１２
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかし、このようなリアルタイム地震危険度予測の方法において、特定のユーザーがその地震危険度を十分に早くかつ精確に推定し、当該の対象物の保全に供するためには、従来の方法では十分とは言えない面があり、その改善の方法が求められている。 
地震の発生と共に、一定の密度をもった地震観測網からリアルタイム地震情報が発信され始める。これ等の情報を主に用いて、高信頼度が要求される緊急対応自動的制御などに供するため、特定のユーザーがその地震危険度を、より早く、十分な精度で、かつ、精度の範囲を明示して推定できるシステムが必要とされている。
【０００５】
ユレダス法の特徴は、検測点が一点の地震波形から震源の情報が求める。この方法では、精度がやや低いこと、震源の深さが求まらないという課題がある。また、着未着法では、迅速性に優れており、２観測点以上のデータを用いて相当に高い精度で震源パラメータが求められるが、２個以上のデータの到達を待つ必要があること、推定解の精度が不明であることが課題である。
【０００６】
また、着未着法により当初それほど多くない複数の観測データを用いて震源地決定を行う場合、同時に異なる地点で発生した地震があった時には、同一の地震に属する地震波と見なされ、推定パラメータに大きな誤差が入る可能性がある。
本発明では、これらの問題を解決する方法を提示する。すなわち、リアルタイム地震危険度予測の方法の課題を解決し、実用性を向上させるため、次の２項を提供する。
１）震源情報の高精度・高速度決定の方法
２）ユーザー側でのリアルタイム高精度地震危険度推定の方法（以下「高度化着未着法」という）
【課題を解決するための手段】
【０００７】
本発明では、新たな方法の提示により、以上に説明した課題に応えるものである。
この目的を達成するため、本発明においては、特許請求の範囲に記載するように構成している。すなわち、本発明は、請求項１に記載のように、
地震波観測データによる震源の時刻および位置の決定方法であって、１個以上の地震波観測データを用い、高度化着未着法によって震源要素および、その精度を推定することを特徴とするリアルタイム地震危険度予測の方法を構成する。
また、本発明は、請求項２に記載のように、
地震波観測データによる震源の時刻および位置の決定方法であって、前記１個以上の地震波観測データを用い、前記高度化着未着法とユレダス法を融合して震源を決定することを特徴とするリアルタイム地震危険度予測の方法を構成する。
また、本発明は、請求項３に記載のように、
地震波観測データによる震源の時刻および位置の決定方法であって、ほぼ同時に２カ所以上で発生した可能性のある地震を、前記高度化着未着法によって識別することを特徴とするリアルタイム地震危険度予測の方法を構成する。
また、本発明は、請求項４に記載のように、
地震波観測データによる震源の時刻および位置の決定方法を使い、特定地区での利用に供するために高精度で、特定サイトの地震危険度を推定することを特徴とするリアルタイム地震危険度予測の方を構成する。
【発明の効果】
【０００８】
本発明に係るリアルタイム地震危険度予測の方法によって得られる効果を請求項ごとに説明する。
まず、請求項１の発明においては、高度化着未着法によって一観測点データを用いて、迅速かつ、かなりの高精度で震源時刻および位置が、その信頼限界と共に求まる。
また、請求項２の発明においては、ユレダス法と高度化着未着法による融合法によって、一つの観測データを用いて、高精度に震源の位置が、その信頼限界と共に決定できる。
また、請求項３の発明においては、許容の値を基準として、ほぼ同時に２ケ所以上で観測された地震波が、同一の地震によるものであるか否かを判定することが可能である。
また、請求項４の発明においては、特定サイトの地震危険度を、高精度に推定することが可能であり、緊急地震速報の利活用に供される。
【発明を実施するための最良の形態】
【０００９】
以下に実施例1〜４を示す。
【実施例１】
【００１０】
図１は、本発明に係る地震観測網の例を示す図である。図１に例示するように、日本全国に観測点数Nの観測網が一定のレベルに整備されているとする。例として高感度地震観測網（Hi-net）や気象庁の津波地震早期検知システム、多機能型地震観測網を想定する。Hi-netでは、現在約800点の観測点がある。間隔は約20kmであるが、北海道と本州の間に津軽海峡があり、そこでは測点の密度は小さいなど観測点密度は一様ではない。又、海の測点数は、現在、若干であるに過ぎない。
【００１１】
地震観測データは連続的に一定時間幅（Hi-netでは１秒）のパケットで伝送されてくるので、図２に例示するように、パケット単位で地震の到達した測点の到着時間が読み取られ、震源パラメータ（位置、震源時）が推定される。又、振幅から地震の規模を示すマグニチュードが推定される。
【００１２】
この発明では、小数点のデータから、震源パラメータを推定する新たな方法を提示する。なお、パケット毎の伝送データの場合の取り扱いは、連続伝送の一つの特殊ケースとして扱えるので、ここでは連続的に伝送されているとする。
【００１３】
地震波が検知された時刻をｔ1、ｔ2、ｔ3、・・・ｔNとし、その測点をS1、S2、S3・・・SNとする（図２）。Nは測点の数である。S1、S2、S3・・・の位置ベクトルをｒ1、ｒ2、ｒ3・・・とする。簡単の為に観測点が水平地表面上にあると仮定し、座標系は３次元直角座標とする。水平面上にｘ、ｙ軸をとり、鉛直上方にｚ軸の正方向をとる。
【００１４】
震源の位置をR、震源時をｔ0とする。地震波の速度をP波に対してはVp、S波に対してはVsとする。ここでは簡単のため一様媒質と仮定し、概ねVp を６km/s、Vsを3.6km/sとする。
（１）１点検測（ｔ1≦ｔ＜ｔ2）の場合
一定の基準で当該イベントが地震によるものと判定されたとする。この判定方法そのものが一つの重要な課題であるが、本稿では取り扱わない。時間がｔ1から第２の観測点に検知される時刻ｔ2までの区間で震源のパラメータが、どのように推定されるかを調べる。
任意の時刻ｔ∈（ｔ1、ｔ2）では、震源をRとして以下の条件が成り立つ。
【００１５】
｜R −r1｜／Vp  ＝  (ｔ1 − ｔ0)                         (1-1)
｜R −r2｜／Vp  ＞  (ｔ − ｔ0)                          (1-2)
・・・・・・・・・・・・・・・・
｜R −ri｜／Vp  ＞  (ｔ − ｔ0)                          (1-i)
・・・・・・・・・・・・・・・・
式(1-1)、(1-i)から、
｜R  ー  rj｜  ー  ｜R  −  r1｜  ＞  Δt・Vp                (2)
２  ≦  ｉ  ≦  N
となる。すなわち、測点Siと測点S1を焦点とする双曲面Ci1の、測点S1を内部に含む双曲面
【００１６】
【数１】
2006105862.tif 000003
【００１７】
の内側に震源があることになる（図３参照）。ここで座標軸を、点SiとS1を結ぶ線をXi軸とし、それに直交して鉛直線上にZ軸、さらに右手示になるようにYi軸をとる。原点Oiは、点S1と点Siを結ぶ線の中点とする。
【００１８】
この場合、式(3)のa、bは、  ｔ  −  ｔ1  ＝  Δtとして
｜ri  −  r1｜  ＝  ２ｌi1
a  ＝  Vp  ・  Δt                                    (4)
【００１９】
【数２】
2006105862.tif 000004
【００２０】
となる。双曲面Ci1の当該片業がX軸と交わる点の座標は、（−ａ、０、０）である。ｔ＝ｔ1の時にはａ＝０で、Ci1は平面（X＝０）となり、震源が点S1、Siの対称面（X＝０）の点S1側の領域（Di1）にあることが言える。任意の観測点Siに関しても、同じような制限が課されるので、震源の存在可能領域はDi1（i=1,2,3,,,,）の共通領域である角柱の内部にあることになる。この場合、深さに関しては、地震活動から知られる制限を使うことになる。
Δｔが大きくなると、Ci1はより鋭い曲面となり、点S1に向かって移動すると共に収縮し、震源の存在できる領域Dilは小さくなる。
【００２１】
【数３】
2006105862.tif 000005
【００２２】
すなわち、震源位置Rの推定精度が高くなる。
Di1(t)（ｉ＝２、・・・、N）の共通領域をD1 (t)とする
【００２３】
【数４】
2006105862.tif 000006
【００２４】
地震発生域をＤと考える。Ｄは鉛直方向には、地表面Z＝０と、最大震源深度ｄmとで挟まれる領域で、水平方向は、地震発生帯で区切られるとして良い。従って全測点（i＝２、・・・、N）の検測結果を使って、震源の存在領域D1(t)は、
【００２５】
【数５】
2006105862.tif 000007
【００２６】
と表される、
【００２７】
【数６】
2006105862.tif 000008
【００２８】
D1(t)は、３次元の立体領域となる。D1を十分な精度で推定するには、点S1を囲み、最も小さな多角形をなす測点群を考えれば十分である。
震源位置の推定値を、立方領域D1(t)の重心と定める。
【００２９】
【数７】
2006105862.tif 000009
【００３０】
信頼限界  σ    ＝  （σx、σy、σz）は、
【００３１】
【数８】
2006105862.tif 000010
【００３２】
【数９】
2006105862.tif 000011
【００３３】
【数１０】
2006105862.tif 000012
【００３４】
と定めることができる。
【００３５】
地震活動度の空間分布S(ｒ)がわかっているとして、それに比例する重みを使って、
より合理的な推定値を定義することもできる。
【００３６】
【数１１】
2006105862.tif 000013
【００３７】
信頼限界は、
【００３８】
【数１２】
2006105862.tif 000014
【００３９】
【数１３】
2006105862.tif 000015
【００４０】
【数１４】
2006105862.tif 000016
【００４１】
と定義できる。
一方、読み取り誤差Δt1による震源位置推定誤差は以下のように見積られる。
【００４２】
【数１５】
2006105862.tif 000017
【００４３】
いま、D1+(t)  ＝  D（ｔ｜ｔ1＋Δｔ1）とすると、Δｔ1＞０の時には、
【００４４】
【数１６】
2006105862.tif 000018
【００４５】
となる。このときの震源位置の推定値をR+(t)とすると、
【００４６】
【数１７】
2006105862.tif 000019
【００４７】
となる。RとR＋の差の絶対値で読み取り誤差によって推定誤差ΔRを定める、
【００４８】
【数１８】
2006105862.tif 000020
【００４９】
【数１９】
2006105862.tif 000021
【００５０】
【数２０】
2006105862.tif 000022
【００５１】
となる。
（２）２点検測（ｔ2≦ｔ＜ｔ3）の場合
ｔが大きくなり２番目の観測点S2にも地震動が到達した後、３番目の点S3に到達する前までの時間帯を対象にする。
【００５２】
この時間帯では、震源は、
【００５３】
【数２１】
2006105862.tif 000023
【００５４】
なる曲面C12上にある。
推定値を求める前に、ｔ＝ｔ2で測点S2にＰ波（初微動）が検出されたことの妥当性の検証を行う。点S2がないとして、点S1での検測、S3、S4・・・では未到着であるとすると、震源の存在できる領域は
【００５５】
【数２２】
2006105862.tif 000024
【００５６】
であり、点S1がないとして、点S2のみの検測とした場合の震源の存在できる領域は、
【００５７】
【数２３】
2006105862.tif 000025
【００５８】
となる。
【００５９】
【数２４】
2006105862.tif 000026
【００６０】
この時には、S1またはS2のＰ波観測に誤りがあったことになる。どれが誤りであるかは不明である。
【００６１】
判定する方法として、地震動による粒子運動の方向と許容領域と測点との相対関係から予想される粒子運動面との整合性を調べることが考えられる。また、Ｐ波か否かについても、粒子運動の極性を使うことが出来る。なお、この粒子運動対照法の確実度を定量化しておくことも必要である。
【００６２】
【数２５】
2006105862.tif 000027
【００６３】
この場合には、震源の存在可能な領域が、
【００６４】
【数２６】
2006105862.tif 000028
【００６５】
と推測される。震源位置の推定値として、領域D12の重心と定めることができる。信頼限界として、領域D12の広がりで定義できる。
（３）３点検測の場合
前項２で述べた方法を順次拡張して行くことで、３点S1、S2、S3検測の条件での、解の推定ができる。３点以上になると条検測の妥当性の判定に多数決法も導入できる。
（４）データ収集が断続的に行われる場合
地震のデータが一定の時間Δtごとのフレーム伝送によって収集される場合には、任意の時間での決定ではなく、t＝Δt、２Δt、３Δt・・・・・という時間ごとに決定することになるのみで、これまでの導き方に変更はない。
（５）簡便な計算法
実用に供するために簡易な計算の方法を呈示する。点S0を囲む４つの測点からなる四辺形S′1S′2S′3S′4の近似長方形をS1S2S3S4とする（図４参照）。図４は、平面図である。ここで近似長方形Ｓ0とＳ1のデータから、曲面C01の左葉のｘ軸負側に震源が制限される。C01を表現する式は、
【００６６】
【数２７】
2006105862.tif 000029
【００６７】
また、C02を表現する式は、
【００６８】
【数２８】
2006105862.tif 000030
【００６９】
となる。ここで、
【００７０】
【数２９】
2006105862.tif 000031
【００７１】
【数３０】
2006105862.tif 000032
【００７２】
【数３１】
2006105862.tif 000033
【００７３】
二つの双曲面C01、C02のZ＝０（地表）上の交点P＋012、P−012の座標を求める（図5参照）。
【００７４】
【数３２】
2006105862.tif 000034
【００７５】
【数３３】
2006105862.tif 000035
【００７６】
但し、二つの解の内
−ｌ02  ＜  X  ＜  ｌ01
の方を選ぶ。式（18a）、（18b）から
【００７７】
【数３４】
2006105862.tif 000036
【００７８】
ここで、ｌ01  ＝  ｌ02  の時には、
X  ＝  ０
となる。又b＝b′であるので、双曲面C01、C02は合同となる。この時は、C01、C02で囲まれる領域の重心のX方向はS0のX座標と一致する。
【００７９】
【数３５】
2006105862.tif 000037
【００８０】
【数３６】
2006105862.tif 000038
【００８１】
となる。他の解は、X＜−ｌ02でこれは使わない。式(18a)にX02を入れて求まるｙの値を±Y02とする。XZ面で対応する点は（X012、−Y02）である。一般の場合には、式(19)の二つの解X1、X2（X1＜X2）の内、X2を使う。このX2を式(18a)に代入して求められるｙをY2とする。この領域C01C02重心のX座標をXとする。
【００８２】
【数３７】
2006105862.tif 000039
【００８３】
積分は台形公式（台形数３ー４個）で求める。範囲は、
【００８４】
【数３８】
2006105862.tif 000040
【００８５】
同じようにS3、S4に対する曲面C03C04から決まる領域C03C04の重心のY座標をInv.(Y)とする。
【００８６】
【数３９】
2006105862.tif 000041
【００８７】
【数４０】
2006105862.tif 000042
【００８８】
（Inv.(X)、Inv.(Y)）を領域C01C02C03C04の重心の水平方向の座標と近似できる。
（ΔZ）については、
【００８９】
【数４１】
2006105862.tif 000043
【００９０】
【数４２】
2006105862.tif 000044
【００９１】
として、絶対値の小さい方をとる。
【００９２】
【数４３】
2006105862.tif 000045
【００９３】
小さい値を与える領域をC0iC0i＋１とすると、この間の最深点座標がY0i（i＋1）であり、Zの範囲を
【００９４】
【数４４】
2006105862.tif 000046
【００９５】
とすることが出来る。Y0i(i＋1)が地震活動度深度ｈＳを越えるときは、不定扱いとする。
境界点の場合
点S0の周囲点４つの内１つまたは２つが欠けている場合である。この時には欠けた点の方向での震源の存在領域の制限に当たっては、地震活動から知られる範囲の知識を使う。
以上のように、震源位置の最大値および、分散を求めることで、一観測点データによる、震源時刻および位置が求まる（高度化着未着法）。その結果、リアルタイム地震危険度の予測の高速度・高精度化が可能となる。
【実施例２】
【００９６】
ユレダス法と高度化着未着法のハイブリッド（融合）法による一つの観測データを用いた震源パラメータの推定方法を説明する。ユレダス法では、地震波の震央到達方向θと、震央距離Δが求まる。これ等の推定誤差をｄθ，ｄΔとする。これから計算される震源の位置をRｕ、確率分布をProb（Rｕ）とする。一方、高度化着未着法による震源の推定位置をR aとし、確率分布をProb（R a）とする。なお、この分布は、ユレダス法では、測定におけるデータのＳ／Ｎ（信号対雑音）比などによって決まり、少数着未着法では、主として観測点数および、場所の関数となる。
ハイブリット法では、ユレダス法と高度化着未着法との結果を重みつき平均で推定する。すなわち、
R h =  λR a＋（1−λ）R u                        （２６）
ここで、正数λは、両方法に対する重みであり、確率分布から求められる値である。
ユレダス法と高度化着未着法との確率分布の分散を、それぞれΔｕ，ΔＡとすると、重みλは、例えば、以下のように定める。
λ＝Δｕ−１／（Δｕ−１＋ΔＡ−１）                  （２７）
震源時間についても、同様な方法で推定する。
以上説明したように、ユレダス法と高度化着未着法によるハイブリッド法によって、一つの観測データを用いて、高精度に震源の位置が決定できる。
【実施例３】
【００９７】
次に、同一の地震によるものであるか否かの判定方法について説明する。
高度化着未着法では、二つ以上の観測データまたは、通常の地震観測では、３点以上の地震観測であるから、その際、データが同一の地震に対応するか否かを、確かめておくことが必要となる。ここでは、その方法を提示する。
同一の地震によるものであるか否かの判定は、第１の実施形態、あるいは第２の実施形態の方法を、地震波を新たな観測点で検知する度に行う。すなわち、第一の測点データによる震源位置、震源時間を（R １，ｔ01）、第二の測点による震源位置、時間を（R ２，
ｔ0２）とする。二つのベクトルの差を
｜ΔR ｜＝｜R １−R ２｜＋Ｖｐ｜ｔ01−ｔ0２｜            （２８）
と定め、許容の値｜ΔR ｜Ｓと比較して、
｜ΔR ｜≦｜ΔR｜Ｓ；同一地震                        （２９ａ）
｜ΔR｜＞｜ΔR ｜Ｓ；異なる地震                      （２９ｂ）
と判定する。｜ΔR ｜Ｓは、それぞれの方法による誤差を基準として、式（２８）になら
って、決めるものとする。統計的検定（帰無仮説）を行って判定してもよい。
３点以上になったときは、多数法を占めるグループに入るか否かを判定する。
以上説明したように、許容の値｜ΔR ｜Ｓを基準として、ほぼ同時に２ケ所で検測された、同一の地震によるものであるか否かを判定することが可能である。
【実施例４】
【００９８】
次に、リアルタイムで特定サイトの地震危険度を推定する方法を説明する。
４．１  危険度推定方法
図６は、危険度推定方法の説明図である。図６に示す例では、地震の震央Ｓ（符号１）とユーザーＵ（符号６（Ｕ））が、ある程度離れている場合（観測網の平均的間隔をＬとして、３Ｌ〜４Ｌ）について説明する。
【００９９】
ある防災緊急対応を行うのに必要な最小時間をＴa（allowed time）とし、準備に要する時間をＴｐ（preparatory time）とする。Ｔａはほぼ１〜数秒、Ｔｐは１０秒前後と考えてよい。というよりこのような特性を有する機器がリアルタイム地震情報の活用に適すると思われる。ユーザーは、主要動の到達前Tａに緊急制御を行い、Ｔｐより前に予備操作を開始する。予備操作は何段階かに分かれていても良い。但し、直近で地震が発生する場合には、到達前情報が、最小時間Ｔａより小さいことがあることに注意するべきである。このような場合には、ユーザーサイド（特定サイト）の地震計データを参照し緊急処置を取ることとなる。
以下、処理方法を説明する。但し全ての観測データが融合して扱える場合を例としてあげるが、ユーザーサイドの地震計データが独立に扱う必要がある時も、ほぼ同様な考え方で処理できる。
【０１００】
図７は、危険度（震源値、時間）の算定アルゴリズムを示す図である。図７に示すように、観測データ、震源情報が逐次送達されてくるものとする。ユーザーゾーン（注意ゾーンＰＺ、動作ゾーンＲＺ）にあるデータをユーザー領域データＤ（ｕ）として、以下、図７の危険度（震源値、時間）の算定アルゴリズムを４種類（Ａ〜Ｄ）について説明する。なお、それがネットワークのものであることもあり、又ユーザーサイドのものであることもあるが、マージされているものとする。こうすることによって公用のデータと、ユーザーの固有のデータを統一的に活用できる。
Ａ．センターからｉ番目の震源情報Ｐｉが出されたとする。その確度がＰｒｏｂ（Ｐｉ）である。その時の危険度（地震強度、余裕時間）を算定する。それらをまとめてＲｉとし、確立分布をProb（Ｒｉ）とする。
Ｂ．Ｐ波がＰＺに到達する予定の時刻まで：
発信される情報のチェックを送信側のアルゴリズムと独立な算出方法で行う。すなわち、広域データＤ（ｒ）とユーザー領域データＤ（ｕ）を使って、推定震源パラメーターの妥当性を検定する。
【０１０１】
（１）変更が必要でないとき：
このときは、確かにまだP波がＰＺに到達していないとして良いので、準備的行動を取るか否か待機する体制となる。但し、安全サイドのリスクマネージメントでは、準備的行動を取る選択もある。
【０１０２】
（２）変更が必要となったとき：
Ｐ波が既にＰＺに到達しているか、または当該地震以外の原因で当該測点にトリガーがかかったことになる。当該トリガーが妥当であるか否かを次の方法で判定する。
(1)近傍の測点に既に地震P波が到達しているか否か
(2)予想される地震強度か否か
(3)見かけ速度に矛盾がないか
(4)波形データが有る場合に質点の運動方向が震源方向と矛盾がないか
Ｃ．Ｐ波がPＺに入った時：
パラメータをチェック、補正する（例、(特願２００１−２５７７６５号公報参照)、発明者（藤縄等））。なお、特願２００１−２５７７６５号では、ユーザー側のデータを用いて補正するとしているが、本発明ではユーザー側と広域観測データとを統一的に扱ってより正しい補正を期す。前提として、測定条件を把握しておく必要がある。
（ａ）観測が正しく行われているか否か。
（ｂ）バックグランドノイズレベルの特性を常時調べる。
この上で、Ｐ波がＰＺに入ったことの真偽を検定し、その確度を求める。予め定めた確度の、しきい値以上で真であるとすれば、Ｐ波がＰＺに入ったことと判定する。
【０１０３】
（ｃ−１）ＰＺ内一点検測のとき：
推定P波振幅との比較：
推定値との整合性が真と判定されれば、この検測を真とする。偽となった場合には、
PＺに検測されたイベントと、当該地震は、同一でないとする。
【０１０４】
（ｃ−２）ＰＺ内2点検測のとき：
(1)推定P波振幅との比較
(2)前述の整合性の判定
これらの判定を双方とも使って総合的に行う。（ｃ−１）の場合と同様に、イベントの同一性がないと判定された場合には、別扱いとなる。ＰＺ内の地震データから求められる危険度Ｒiが一定の閾値以上なら準備操作信号を出す。この操作は、補正震源パラメータを用いて、Ｐ波が動作ゾーンＲＺに入ると思われる時間まで行う。
【０１０５】
（ｃ−３）別イベントの時：
同時に２つ以上の地震が発生したことになる。ＰＺ内の地震データのみを用い、先に提示した少数データによる震源パラメータの推定を行なって、危険度Ｒiが一定の閾値以上なら余裕時間を勘案して、準備操作信号単独または緊急制御も殆ど同時に出す。
【０１０６】
Ｄ．Ｐ波が危険領域Ｒｚ（図６の７）に入った後：        
Ｒｚ内の二つ以上測点で検知され、地震危険度及び確度が一定レベル以上の時には、制御開始信号を出す。
以上、方法Ａ〜ＤをＳ波がＲｚに到達するまで繰り返す。
ユーザーサイドの地震観測網
以上の結果を用いて、ユーザー地点Ｕに当該施設があり、この施設を地震から守るための地震観測網を設計する。ユーザー地点からＲＺまでの距離をρｒ，ＰＺまでの距離を
ρｐとする。警戒時間をＴｒ，準備時間をＴｐとすると、
ρｒ＝ＶｓｘＴｒ                                  （３０a）
ρp＝ＶｓｘＴｒ                                  （３０b）
と定める。
大体の値は、Vs＝３．6ｋｍ／ｓ、Ｔｒ＝１ｓとして、ρｒは３．１ｋｍとなる。
【０１０７】
そして、ＰＺに一定の地震計を配置して、公式発表データの検証を行うこととする。地震計Ｕを中心とする円形にかつ、方位角θで均等に並べるとする。正ｎ角形の頂点に配置すると、
ρｒ／ρｒ'＝ｃｏｓθ                              （３１）
θ＝２π／ｎ
となる。
たとえば、θ＝６０°とすると、ｎ＝６、ρｒ／ρp＝１／２となり、ρp≒７．２ｋｍとなり、ほぼ合理的な距離に設置できる。ちなみに、ｎ＝4では、ρp＝∞と、ほとんど非実用的な距離になる。
このようにすると、任意の到来方向の地震到来に対して、少なくとも２点でのデータを使用して準備操作ができる。この場合の１点検測による予備時間Ｔｐ（１）は、
√3／２ｓｅｃ≦Ｔｐ（１）≦１ｓｅｃ
また、２点検測による準備時間Ｔｐ（２）は、
０≦Ｔｐ（２）≦√3／２ｓｅｃ
の範囲となる。
さらに、Ｔｐ（２）をＴｐ'としたい場合には、
Δρｒ＝Ｔｐ'ｘＶＳ
だけ広げる。これに応じて、ρp はΔ２ρｒだけ大きくする。なお、対象をシナリオ地震に限る場合には、当該地震の方向をカバーすればよいので、その分、必要な地震計の数は少なくなる。
４．３  地震危険度推定値の評価
推定値が妥当か否かの判定が必要になることが多い。それの基本となる考え方を説明する。図８は、地震の強度（通常は計測震度が使われるが、対象に応じて、パラメータが選択される）の推定の評価グラフを示す。図８に示す例では、計測震度を使っている。先の広いロート状の領域（Oabcd）が合格の範囲としている。被害発生最小震度Iｔｍ以下の震度では、被害が無いという結果になればよいので、推定値がOABCのなかに入れば合格と言うことになるが、最大無被害震度などは、対象によって異なること、被害が震度に依存することを考慮すれば、先の広いロート状の領域を合格範囲とする合理性はあるであろう。
一方、図９は、地震危険度の予測値の評価グラフのうち、地震危険度Ｒの今ひとつの指標である余裕時間に関するものである。緊急防災対応を行うのに必要な最小時間（Ｔｒ）、準備行動を取るのに必要な時間（Ｔｐ）を使って、余裕時間の合否範囲の例を示している。
【０１０８】
正しい到達時間をＴ（ｔ）とし、推定値をＴ（e）、誤差をΔＴとすると、
Ｔ（e）＝Ｔ（ｔ）＋ΔＴ                                （３２）
となり、場合として3つが考えられる。
１）Ｔｒ≧Ｔ（ｔ）
この場合には、主要動が到達するまでに終えたい操作を時間内に完了させることが出来ない。この場合には、対処として、対象システムの特性にあわせて、何も行わない、あるいは、完了しないまでも、制御を開始することが選択される。
２）Ｔｐ≧Ｔ（ｔ）≧Ｔｒ： 
このケースでは、準備操作には間に合わないが、緊急操作を行うには、時間の余裕が有ることになる。推定の結果が閾値以上で有れば、緊急対応操作を行う。
３）Ｔ（ｔ）≧Ｔp:
このケースでは、準備操作も行える。推定の結果がそれぞれ閾値以上で有れば、準備操作および緊急対応操作を行う。
以上説明したように、入力された震源情報をリアルタイムに処理して、特定サイトの地震危険度を推定することが可能である。
【図面の簡単な説明】
【０１０９】
【図１】本発明に係る地震観測網を示す図である。
【図２】図１に示した地震観測網による地震波の説明図である。
【図３】図２に示した地震波が一点で観測され時に、震源が存在できる範囲を求める方法の説明図である。
【図４】図２に示した地震波が一点（S0）で観測され時に、震源が存在できる範囲を近似的に求める為に擬似的な観測点（S１、S２、S３、S４）による詳細な説明図である。
【図５】図４に示した近似観測点配置の場合の、震源存在領域推定の説明図である。
【図６】ユーザー側での地震危険度高精度推定方法についての詳細な説明図である。
【図７】図６に示した地震危険度の予測値の評価方法を示す説明図である。
【図８】地震の強度推定の評価グラフを示す説明図である。
【図９】図８に示した地震強度推定の評価グラフの余裕時間に関する説明図である。
【符号の説明】
【０１１０】
１  震央
２〜５  震央近傍の地震波観測点
６（ｕ）  ユーザーの場所
７  緊急指令を出す時点の主要動の波面の包絡線（Ｒｚ）
８  準備指令を出す時点の主要動の波面の包絡線（Ｏｚ）
１０〜１５  ユーザー近傍の地震波観測点
２１  地震波観測点（最初の地震波到達点）
２２  ある時間の地震波が最初の観測点２１に到達した時の波面の位置
２３〜３２  地震波が未到達の観測点
【図１】
2006105862.tif 000047
【図２】
2006105862.tif 000048
【図３】
2006105862.tif 000049
【図４】
2006105862.tif 000050
【図５】
2006105862.tif 000051
【図６】
2006105862.tif 000052
【図７】
2006105862.tif 000053
【図８】
2006105862.tif 000054
【図９】
2006105862.tif 000055