JP7463187B2 - 震度推定システム及び震度推定方法 - Google Patents

Description

本発明は、震度推定システム及び震度推定方法に関するものである。

特許文献１には、地震計に関する技術が開示されている。この地震計は、計測部、記録部及び演算部を備える。計測部は、観測地点において地動を観測し、地動の地震動を検出する。記録部には、過去に発生した地震データのＳ波の加速度の最大値と、Ｐ波が到達してから連続する複数の予測期間における地震データのＰ波の加速度の最大値との相関関係が、各予測期間ごとに記録されている。演算部は、各予測期間における相関関係と、計測部により計測された地震動のＰ波の加速度の最大値とを用いて、観測地点に到達する地震動のＳ波の加速度の最大値を各予測期間が経過するごとに推定する。

特許文献２には、地震報知システムに関する技術が開示されている。この地震報知システムは、震度計及びＦＭ多重表示装置を備える。震度計は、地震波のうち初期微動を検知することにより、地震の予測情報を演算し、演算された予測情報に基づいて地震警報信号を出力する。ＦＭ多重表示装置は、震度計に接続され、ＦＭ多重放送の文字情報をスクロール表示する。ＦＭ多重表示装置は、地震警報信号を入力したときに、地震に関する地震文字情報をＦＭ多重放送の文字情報より優先して表示させる。

特開２０１８－０４４７８４号公報 特開２０１１－０６５４９６号公報

近年の地震検知システムにおいては、地震が発生すると、分散配置された複数の観測地点におけるＰ波の振幅および到達時刻から震源位置及び地震規模（マグニチュード）を算出し、その震源位置を基準とする距離減衰の経験式に基づいて、各地の震度を地震発生直後に推定する。しかしながら、このような方式では、震源位置からの距離が遠い地点ほど推定誤差が大きくなるという問題がある。実際、多くの地点において震源位置からの距離は数十キロないし数百キロに及び、推定誤差を無視できない場合がある。また、地震計を設置している観測地点においては、地震計により検出されたＰ波の加速度の最大値に基づく経験式からＳ波の加速度の最大値を推定することも可能であるが（特許文献１を参照）、観測地点ではない地点における正確な震度を推定することが求められる場合もある。例えば、鉄道設備においては、地震が発生した際、鉄道の運行を停止するか否かを正確に決定する為に、各駅における震度、及び駅間の線路上における震度を精度良く推定することが望まれる。通常、これらの鉄道設備においては鉄道の走行による地面の揺れが頻繁に発生するので、これらの鉄道設備は地震計の設置に適さない。

本発明は、このような問題点に鑑みてなされたものであり、観測地点ではない地点における震度を精度良く推定することが可能な震度推定システム及び震度推定方法を提供することを目的とする。

上述した課題を解決するために、本発明の震度推定システムは、互いに離れた複数の観測地点のそれぞれにおいて地動を観測し、地動の地震動を検出する複数の地震計と、Ｐ波最大加速度値を含む地震動に関する情報を複数の地震計から受ける上位装置と、を備える。上位装置においては、各観測地点を中心とする複数の地域区画が予め設定されている。複数の地域区画のうち二以上の地域区画のそれぞれには、震度を予測すべき対象施設が含まれている。上位装置は、各観測地点にＳ波が到達するまでの間、Ｐ波最大加速度値を各地震計から周期的に受ける。そして、上位装置は、Ｐ波最大加速度値に基づいて各観測地点における工学基盤の震度又はＳ波最大加速度値を推測演算し、各地域区画における工学基盤の震度又はＳ波最大加速度値が個々の地域区画内において均一であるとの仮定の下、各地域区画内における表層地盤増幅率の分布に基づいて、地表面の震度を各地域区画内において面的に推定する。

上述した課題を解決するために、本発明の震度推定方法は、互いに離れた複数の観測地点のそれぞれを中心とする複数の地域区画を予め設定する設定ステップと、各観測地点において地動を観測し、地動の地震動を検出する検出ステップと、各観測地点にＳ波が到達するまでの間、各観測地点のＰ波最大加速度値に基づいて各観測地点における工学基盤の震度又はＳ波最大加速度値を推測演算し、各地域区画における工学基盤の震度又はＳ波最大加速度値が個々の地域区画内において均一であるとの仮定の下、各地域区画内における表層地盤増幅率の分布に基づいて、地表面の震度を各地域区画内において面的に推定する推定ステップと、を含む。各地域区画は、隣り合う少なくとも一つの地域区画と接するか又は縁部同士において重なり合っている。複数の地域区画のうち二以上の地域区画のそれぞれには、震度を予測すべき対象施設が含まれている。

これらの震度推定システム及び震度推定方法によれば、地表面の震度を各地域区画内において面的に推定するので、観測地点の周辺における観測地点ではない地点においても震度を精度良く推定することが可能であり、各地域区画に含まれる予測対象施設の震度を精度良く推定することができる。また、各地域区画内の単一の観測地点における観測結果（Ｐ波最大加速度値）に基づいて当該地域区画内の震度を推定するので、複数の観測地点における観測結果に基づいて震源位置及び地震規模を算出してから各地の震度を推定する方式と比較して、地震が発生してから震度を推定するまでの所要時間を短くすることができる。

上記の震度推定システム及び震度推定方法において、各地域区画は、隣り合う少なくとも一つの地域区画と接するか又は縁部同士において重なり合っていてもよい。この場合、連続して延在する広い地域において震度を面的に推定できるので、例えば鉄道設備といった、長い距離にわたって連続する設備の各地点における震度を精度良く推定することができる。

上記の震度推定システムは、複数の地域区画における地表面の震度の分布を視覚的に表示する一又は複数の表示部を更に備えてもよい。この場合、各地域区画内において面的に推定された地表面の震度を、視覚情報によって容易に且つ短時間で理解することができる。

上記の震度推定システムは、複数の地域区画における地表面の震度の分布に関するデータを記録する記録部を更に備えてもよい。この場合、各地域区画内において面的に推定された地表面の震度に関するデータを、地震の解析に用いることができる。

上記の震度推定システムは、予測震度別に対象施設をグループ分けしたデータを出力する出力部を更に備えてもよい。この場合、推定震度が大きい（すなわち地震被害が大きいと推定される）対象施設を容易に且つ短時間で判別することができる。

上記の震度推定システム及び震度推定方法において、複数の地域区画の形状は四角形であってもよい。日本の地震ハザードステーションが提供する表層地盤増幅率のデータの区画は四角形状であるので、地域区画の形状が四角形であることによって演算を容易にできる。

上記の震度推定システム及び震度推定方法において、複数の地域区画の面積は互いに等しくてもよい。この場合、各地域区画内における震度の推定時間及び推定精度を複数の地域区画間で均一に近づけることができる。

上記の震度推定システム及び震度推定方法において、隣り合う地域区画の縁部同士が重なり合っている地域においては、隣り合う地域区画それぞれにおいて推定された地表面の震度のうち値が最も大きいものを用いてもよい。この場合、当該地域において地震被害を過小に推定することを回避できる。

上記の震度推定システムにおいて、上位装置は、各地震計から周期的に受けるＰ波最大加速度値から各観測地点におけるＳ波最大加速度値を推測演算し、Ｓ波最大加速度値から各観測地点における震度の推定を行い、該震度に基づいて各観測地点における工学基盤の震度を推測演算してもよい。同様に、上記の震度推定方法の推定ステップでは、各観測地点から周期的に受けるＰ波最大加速度値から各観測地点におけるＳ波最大加速度値を推測演算し、Ｓ波最大加速度値から各観測地点における震度の推定を行い、該震度に基づいて各観測地点における工学基盤の震度を推測演算してもよい。この場合、Ｐ波最大加速度値に基づく工学基盤の震度の推測演算を精度良く行うことができる。

上記の震度推定システム及び震度推定方法において、対象施設は鉄道設備であってもよい。上述したように、本発明の震度推定システム及び震度推定方法によれば、鉄道設備のような広い範囲に延在する設備であっても、対象施設の震度を精度良く推定することができる。したがって、地震が発生した際、鉄道の運行を停止するか否かの正確な決定に寄与できる。

本発明による震度推定システム及び震度推定方法によれば、観測地点ではない地点における震度を精度良く推定することができる。

本発明の一実施形態にかかる震度推定システムの構成を示すブロック図である。 上位装置の演算部において予め設定されている、各観測地点を中心とする複数の地域区画を示す図である。 Ｓ波最大加速度値と、複数の予測期間におけるＰ波最大加速度値との相関関係を示す散布図である。 Ｓ波最大加速度値と、複数の予測期間におけるＰ波最大加速度値との相関関係を示す散布図である。 Ｓ波最大加速度値と、複数の予測期間におけるＰ波最大加速度値との相関関係を示す散布図である。 地震動のＳ波最大加速度値を推定する方法を説明するためのグラフである。 或る地域区画における、地表面と工学基盤とを模式的に示す図である。 表示部における表示例を示す図である。 印刷部から出力される印刷物の内容例を示す。 一実施形態による震度推定方法を示すフローチャートである。

以下、添付図面を参照しながら本発明による震度推定システム及び震度推定方法の実施の形態を詳細に説明する。なお、図面の説明において同一の要素には同一の符号を付し、重複する説明を省略する。

図１は、本発明の一実施形態にかかる震度推定システム１の構成を示すブロック図である。図１に示されるように、本実施形態の震度推定システム１は、複数の地震計２と、上位装置３と、情報端末４と、通信回路５とを備えている。

複数の地震計２は、互いに離れた複数の観測地点のそれぞれにおいて地動を観測し、地動の地震動を検出する。本実施形態の地震計２は、観測部２１と、処理部２２とを有する。観測部２１は、３つの加速度計２３と、Ａ／Ｄ変換器を含む計測基板２４とを含んで構成されている。加速度計２３は、観測地点の表層において地動を観測し、地動の地震動を検出する。具体的には、３つの加速度計２３が、東西（Ｅ－Ｗ）方向、南北（Ｎ－Ｓ）方向、上下（Ｕ－Ｄ）方向の地震動の加速度を計測する。３つの加速度計２３は、上記３方向の地震動の加速度の大きさを示すアナログ信号を、所定のサンプリング周期に同期させて計測基板２４に出力する。加速度計２３は、サーボ加速度計であってもよく、或いは、サーボ加速度計に代えてＭＥＭＳ（Micro Electro Mechanical System）加速度センサであってもよい。なお、サーボ加速度計或いはＭＥＭＳ加速度センサを用いることにより、上記３方向の地震動の加速度について所定の精度が得られることが望ましい。また、加速度計２３は、上記３方向の加速度の大きさを計測できれば、４つ以上設けられてもよい。

計測基板２４は、３つの加速度計２３に電気的に接続されている。計測基板２４のＡ／Ｄ変換器は、アナログ信号をデジタル信号（加速度データ）に変換する。Ａ／Ｄ変換器は、所定の周期で加速度データを生成する。Ａ／Ｄ変換器は、加速度データを所定のデータに細分化したのち処理部２２に出力する。

処理部２２は、計測基板２４と電気的に接続されている。処理部２２には、計測基板２４のＡ／Ｄ変換器から出力された加速度データが入力される。処理部２２は、地震動の発生を監視する。具体的には、処理部２２は、観測データ処理部２５を含んで構成されている。観測データ処理部２５は、加速度データから地震動のＰ波及びＳ波の到達を判定する。観測データ処理部２５は、地震動のＰ波又はＳ波が到達したと判定した後、加速度データの解析処理を行う。観測データ処理部２５は、地震動のＰ波が到達したと判定した時刻以降の加速度データから、地震動のＰ波の加速度の最大値を、トリガ時刻以降の連続する複数の予測期間ごとに抽出する。また、観測データ処理部２５は、地震動のＳ波が到達したと判定した時刻（トリガ時刻）以降の加速度データから、地震動のＳ波の加速度の最大値を、トリガ時刻以降の連続する複数の計測期間ごとに抽出する。処理部２２は、観測データ処理部２５において抽出したＰ波最大加速度値又はＳ波最大加速度値を含む地震動に関する情報（以下、地震情報という）を、通信回路５を介して上位装置３へ出力する。なお、観測データ処理部２５がＰ波及びＳ波のいずれでもないと判定した場合には、処理部２２は、上位装置３への地震情報の出力を行わない。処理部２２が地震情報を出力する周期は、例えば１秒以上３秒以下であり、一実施例では１秒である。

上位装置３は、演算部３１と、入出力部３２と、表示部３３と、操作部３４と、時刻校正部３５とを有する。演算部３１は、通信回路５と接続されており、通信回路５を介して複数の地震計２から地震情報を周期的に受ける。演算部３１は、各観測地点にＳ波が到達するまでの間、Ｐ波最大加速度値に基づいて、各観測地点及びその周辺の地表面のＳ波最大加速度値及び震度を推定する。また、演算部３１は、各観測地点にＳ波が到達して以降、Ｓ波最大加速度値に基づいて、各観測地点及びその周辺の地表面の震度を推定する。演算部３１は、例えばＣＰＵ及び記憶装置（ＲＯＭ、ハードディスク等）を含むコンピュータにより構成され得る。演算部３１の各機能は、記憶装置に記憶されたプログラムをＣＰＵが読み込んで実行することにより実現される。

ここで、図２は、上位装置３の演算部３１において予め設定されている、各観測地点Ｑを中心とする複数の地域区画Ｒを示す図である。図２には、震度を予測すべき対象施設としての鉄道設備（鉄道駅Ｔａ、鉄道線路Ｔｂおよび変電所Ｔｃ）が併せて示されている。また、図２には、図の中心を震源と仮定したときの、Ｐ波の１秒ごとの平均的な広がりを示す同心円Ｃ１～Ｃ３が併せて示されている。複数の鉄道駅Ｔａのうち、２つの路線が結合する鉄道駅Ｔａ１は、特に重要な対象施設である。他の鉄道駅Ｔａ２は、鉄道駅Ｔａ１に次いで重要な対象施設である。各観測地点Ｑには、地震計２が設置される。観測地点Ｑは、当該地域区画Ｒ内の対象施設に近接する（但し、一致しない）位置に設定される。

各地域区画Ｒは、例えば多角形または円形といった任意の様々な形状とされ得る。一例では、各地域区画Ｒは四角形であり、図示例の各地域区画Ｒは正方形である。複数の地域区画Ｒの面積は互いに等しいことが好ましいが、互いに異なっていてもよい。地域区画Ｒの形状が正方形である場合、一辺の長さは例えば０．２５ｋｍ以上１ｋｍ以下である。複数の地域区画Ｒのうち二以上の地域区画Ｒのそれぞれには、震度を予測すべき対象施設である鉄道駅Ｔａ及び鉄道線路Ｔｂが含まれている。なお、図示例では全ての地域区画Ｒに鉄道駅Ｔａ及び鉄道線路Ｔｂの双方が含まれているが、一部の地域区画Ｒには鉄道駅Ｔａ及び鉄道線路Ｔｂのうち鉄道線路Ｔｂのみが含まれてもよい。各地域区画Ｒは、隣り合う少なくとも一つの地域区画Ｒと接するか、又は、隣り合う少なくとも一つの地域区画Ｒと縁部同士において重なり合っている。

また、図２に示されるように、複数の観測地点Ｑを一本の線Ｌにて結んだとき、その線Ｌ上において互いに隣り合う観測地点Ｑ間の距離は全て等しい。言い換えると、線Ｌ上において互いに隣り合う観測地点Ｑ同士を線分で結んだとき、全ての線分の長さは互いに等しい。線Ｌ上において互いに隣り合う観測地点Ｑ間の距離は、例えば１ｋｍ以上６ｋｍ以下である。対象施設が鉄道設備である場合、その距離は例えば５ｋｍである。なお、複数の観測地点Ｑを一本の線Ｌにて結ぶ場合、その結び方（観測地点Ｑの結合順序）は複数通りあるが、そのうちの少なくとも一通りの結び方においてこの条件を満たしていればよい。線Ｌ上には、通信回線が敷設されてもよい。

演算部３１は、各観測地点ＱにＳ波が到達するまでの間、各地震計２から周期的に受けるＰ波最大加速度値に基づいて、まず各観測地点ＱにおけるＳ波最大加速度値を推測演算する。次に、演算部３１は、各観測地点ＱにおけるＳ波最大加速度値に基づいて、各観測地点Ｑにおける震度の推定を行う。そして、演算部３１は、推定した震度に基づいて、各観測地点Ｑの直下における工学基盤の震度を推測演算する。或いは、演算部３１は、各観測地点ＱにおけるＳ波最大加速度値に基づいて、各観測地点Ｑの直下における工学基盤のＳ波最大加速度値を推測演算してもよい。演算部３１は、各地域区画Ｒにおける工学基盤の震度又はＳ波最大加速度値が個々の地域区画Ｒ内において均一であるとの仮定の下、各地域区画Ｒ内における表層地盤増幅率（ＡＲＶ）の分布に基づいて、地表面の震度を各地域区画Ｒ内において面的に推定する。以下、このような演算部３１の処理について詳細に説明する。

まず、演算部３１は、各観測地点ＱにおけるＰ波最大加速度値に基づいて、当該観測地点ＱにおけるＳ波最大加速度値の推測演算を行う。この処理において、演算部３１は、各地震計２から周期的に受けるＰ波最大加速度値、及び、過去に発生した地震データのＰ波最大加速度値とＳ波最大加速度値との相関関係を示す関係式を用いて、Ｓ波最大加速度値を各周期毎に推定する。この関係式は、演算部３１が有する記憶装置に記憶されている。或いは、この関係式は、入出力部３２が有する記録部３２ａ（後述）に記憶されてもよい。

この相関関係は、過去に発生した複数の地震における、Ｓ波最大加速度値と、Ｐ波が到達してから連続する複数の予測期間におけるＰ波最大加速度値との相関関係である。なお、過去に発生した複数の地震に関するデータとしては、例えば国立研究開発法人防災科学技術研究所（ＮＩＥＤ）が所有する国内地震の地震データを用いることができる。

図３（ａ）、図３（ｂ）、図４（ａ）、図４（ｂ）、図５（ａ）、及び図５（ｂ）は、Ｓ波最大加速度値と、複数の予測期間Ｔ_１～Ｔ_６におけるＰ波最大加速度値との相関関係を示す散布図である。図３（ａ）はＰ波の到達後から０．５０秒経過時までの予測期間Ｔ_１、図３（ｂ）はＰ波の到達後０．５１秒経過時から１．００秒経過時までの予測期間Ｔ_２、図４（ａ）はＰ波の到達後１．０１秒経過時から２．００秒経過時までの予測期間Ｔ_３、図４（ｂ）はＰ波の到達後２．０１秒経過時から３．００秒経過時までの予測期間Ｔ_４、図５（ａ）はＰ波の到達後３．０１秒経過時から４．００秒経過時までの予測期間Ｔ_５、図５（ｂ）はＰ波の到達後４．０１秒経過時から５．００秒経過時までの予測期間Ｔ_６における関係を示す。これらの散布図において、横軸はＰ波最大加速度値（単位：Ｇａｌ、ただし１Ｇａｌ＝０．０１ｍ／ｓ^２）の対数値を表し、縦軸はＳ波最大加速度値（単位：Ｇａｌ）の対数値を表している。グラフＧ１０～Ｇ１５は、これらの散布図の直線近似式である。これらの散布図に示されるように、Ｓ波最大加速度値Ｙと、複数の予測期間Ｔ_１～Ｔ_６におけるＰ波最大加速度値Ｘとの間には、有意な相関関係が存在する。

相関関係を示す関係式は、例えばこれらの散布図において、グラフＧ１０～Ｇ１５に示される直線近似式である。関係式は、次の数式のように表される。

係数ａ及びｂは、Ｓ波最大加速度値ＹとＰ波最大加速度値Ｘとを直線近似したときに得られる係数である。なお、これらの散布図に示されるように、地震データのＰ波が到達してから、予測期間Ｔ_６に近づくほど、相関関係は収束している。すなわち、地震データのＰ波が到達してから時間が経過するほど、Ｓ波最大加速度値Ｙを精度良く推定することができる。演算部３１は、各地震計２から提供されたＰ波最大加速度値を、（１）式のＰ波最大加速度値Ｘに当て嵌めることにより、のちに観測地点Ｑに到達する地震動のＳ波最大加速度値を推定することができる。

ここで、地震動のＳ波最大加速度値を推定する具体的な方法について説明する。図６は、地震動のＳ波最大加速度値を推定する方法を説明するためのグラフである。図６において、横軸は時間（単位：秒）を表し、縦軸は地震動のＰ波及びＳ波の加速度（単位：Ｇａｌ）を表している。Ｔ_Ｐは、地震動のＰ波が到達した時刻であり、Ｔ_Ｓは、地震動のＳ波が到達した時刻である。Ｓ_ＭＡＸは、地震動のＳ波最大加速度値を示す。図６において、時刻Ｔ_Ｐから時刻Ｔ_Ｓの間には、連続する複数の予測期間Ｔ_１～Ｔ_４が含まれる。すなわち、各予測期間Ｔ_１～Ｔ_４は、地震動のＰ波が到達してから地震動のＳ波が到達するまでのＰ波の継続期間に含まれる。図６に示されるように、地震動のＳ波最大加速度値Ｓ_ＭＡＸに対して、各予測期間Ｔ_１～Ｔ_４における地震動のＰ波最大加速度値は、それぞれａ１倍、ｂ１倍、ｃ１倍、ｄ１倍である。ａ１、ｂ１、ｃ１、ｄ１は互いに異なる。したがって、各予測期間Ｔ_１～Ｔ_４に対応する（１）式をそれぞれ用いることにより、地震動のＳ波最大加速度値を推定することができる。具体的には、各予測期間Ｔ_１～Ｔ_４に対応する（１）式に、各予測期間Ｔ_１～Ｔ_４における地震動のＰ波最大加速度値Ｘをそれぞれ代入することによって、前述したように各予測期間Ｔ_１～Ｔ_４ごとに地震動の時間が経過するほど、Ｓ波最大加速度値Ｙを精度良く推定することができる。

次に、演算部３１は、各観測地点Ｑに関する推定したＳ波最大加速度値に基づいて、各観測地点Ｑにおける震度の推定を行う。Ｓ波最大加速度値に基づく震度の算出は、震度計算アルゴリズムを用いて行うことができる。震度計算アルゴリズムは、気象庁により制定された、震度計が満たすべき性能の技術基準に基づいて定められている。

続いて、演算部３１は、推定した震度に基づいて、各観測地点Ｑの直下における工学基盤の震度を推測演算する。ここで、図７は、或る地域区画Ｒにおける、地表面Ｇａと、工学基盤Ｇｂとを模式的に示す図である。図７には、更に、地表面Ｇａにおいて一つの地域区画Ｒを分割して成る複数の単位区画Ｒａと、各単位区画Ｒａに対して個別に与えられている表層３０ｍの平均Ｓ波速度データ（ＡＶＳ３０、防災科学技術研究所提供の地震ハザードステーション（Ｊ－ＳＨＩＳ）の観測点データとして公開）Ｄａとが模式的に示されている。各単位区画Ｒａの平面形状は正方形であり、各単位区画Ｒａの一辺の長さは例えば２５０ｍである。地域区画Ｒにおいて、単位区画Ｒａは複数行及び複数列の二次元状に配列されている。各単位区画Ｒａにおける表層地盤増幅率ＡＲＶは、ＡＶＳ３０に基づいて、経験的に、下記の数式（２）により求められる。

また、工学基盤Ｇｂにおける震度Ｉｒｂは、地表面Ｇａにおける震度Ｉｒａと表層地盤増幅率ＡＲＶとに基づいて、経験的に、下記の数式（３）により求められる。

演算部３１は、数式（２）を用いて、観測地点Ｑにおける表層地盤増幅率ＡＲＶを算出する。そして、演算部３１は、既に推定した観測地点Ｑにおける震度Ｉｒａを数式（３）に当て嵌めることにより、観測地点Ｑの直下の地点Ｑｂにおける震度Ｉｒｂを推定する（図中の矢印Ａ１）。

或いは、演算部３１は、観測地点Ｑにおける推定したＳ波最大加速度値に基づいて、各観測地点Ｑの直下における工学基盤ＧｂのＳ波最大加速度値を推測演算してもよい。この演算もまた、上記数式（３）に準じた関係によって行われる。この場合、各観測地点Ｑに関する推定したＳ波最大加速度値に基づく、各観測地点Ｑにおける震度の推定は不要である。

続いて、演算部３１は、各地域区画Ｒにおける工学基盤Ｇｂの震度Ｉｒｂ又はＳ波最大加速度値が個々の地域区画Ｒ内において均一であるとの仮定の下、各地域区画Ｒ内における表層地盤増幅率ＡＲＶの分布に基づいて、地表面Ｇａの震度Ｉｒａを各地域区画Ｒ内において面的に推定する。震源からの地震波は、地震基盤及び工学基盤Ｇｂを経由して地表面Ｇａの観測地点Ｑに到達する。工学基盤Ｇｂまでは比較的均質であるので、工学基盤Ｇｂでの震度Ｉｒｂ及びＳ波最大加速度値は、或る限定された区画（地域区画Ｒ）の範囲内ではほぼ均一であると仮定することができる。また、面的に推定するとは、地域区画Ｒを分割して成る複数の単位区画Ｒａのそれぞれにおいて震度Ｉｒａを推定することにより、地域区画Ｒの全域において震度Ｉｒａを推定することをいう。

すなわち、演算部３１は、数式（２）を用いて、観測地点Ｑを含む単位区画Ｒａを除く他の単位区画Ｒａにおける表層地盤増幅率ＡＲＶを算出する。これにより、地域区画Ｒ内における表層地盤増幅率ＡＲＶの分布が得られる。そして、演算部３１は、算出した各単位区画Ｒａの表層地盤増幅率ＡＲＶと、工学基盤Ｇｂの震度Ｉｒｂとを数式（３）に当て嵌めることにより（すなわち地域区画Ｒ内における表層地盤増幅率ＡＲＶの分布に基づいて）、各単位区画Ｒａの震度Ｉｒａを推定する（図中の矢印Ａ２）。或いは、演算部３１は、算出した各単位区画Ｒａの表層地盤増幅率ＡＲＶと、工学基盤ＧｂのＳ波最大加速度値とを数式（３）に準ずる数式に当て嵌めることにより、各単位区画Ｒａの地表面ＧａにおけるＳ波最大加速度値を推定し、そのＳ波最大加速度値から各単位区画Ｒａの震度Ｉｒａを求める。演算部３１は、このような推定演算を、全ての地域区画Ｒについて行う。図２に示された対象設備（鉄道駅Ｔａ、鉄道線路Ｔｂおよび変電所Ｔｃ）はいずれかの単位区画Ｒａに必ず含まれるので、このような推定演算により、対象設備における震度Ｉｒａを推定することができる。

なお、各観測地点にＳ波が到達した後においては、演算部３１は、各地震計２から周期的に提供されるＳ波最大加速度値に基づいて、各観測地点Ｑにおける震度Ｉｒａを算出する。そして、演算部３１は、算出した震度Ｉｒａに基づいて（又は、各観測地点ＱにおけるＳ波最大加速度値に基づいて）、各観測地点Ｑの直下における工学基盤の震度Ｉｒｂ（又はＳ波最大加速度値）を演算する。演算部３１は、各地域区画Ｒにおける工学基盤の震度Ｉｒｂ（又はＳ波最大加速度値）が個々の地域区画Ｒ内において均一であるとの仮定の下、各地域区画Ｒ内における表層地盤増幅率ＡＲＶの分布に基づいて、各単位区画Ｒａの震度Ｉｒａを推定（すなわち各地域区画Ｒ内において震度Ｉｒａを面的に推定）する。

また、図２において隣り合う地域区画Ｒの縁部同士が重なり合っている地域においては、隣り合う地域区画Ｒそれぞれにおいて推定された地表面Ｇａの震度Ｉｒａのうち値が最も大きいものを、その地域の震度とする。

再び図１を参照する。入出力部３２は、通信回路５を介して複数の情報端末４と通信を行うための通信インターフェイスである。また、入出力部３２は、上位装置３が有する演算部３１、表示部３３、操作部３４、及び時刻校正部３５と電気的に（例えばデータバスを介して）接続されている。更に、入出力部３２は、警報部６及び外部表示部７と電気的に（又は無線通信手段を通じて）接続されている。入出力部３２は、これらの構成要素間の通信に用いられる通信プロトコルの変換、及びこれらの構成要素に入出力する信号を論理レベル変換する機能を有する。入出力部３２は、演算部３１が推定した震度Ｉｒａに関するデータを演算部３１から受け取る。そして、入出力部３２は、震度Ｉｒａに関するデータを、通信回路５を介して複数の情報端末４へ出力し、更に表示部３３及び外部表示部７へ出力する。

入出力部３２は、記録部３２ａを含む。記録部３２ａとしては、例えばＣＦカードなどの記憶媒体が用いられる。記録部３２ａは、演算部３１により推測された、複数の地域区画Ｒにおける地表面Ｇａの震度Ｉｒａの分布に関するデータを記録する。なお、記録部３２ａに記録されるデータは、例えばパソコンを用いて閲覧が可能なフォーマットにて記録される。このように、震度分布に関するデータを記録することにより、当該データの二次的利用や、震度に基づく地域または施設のグループ分けが容易となる。

表示部３３は、複数の地域区画Ｒにおける地表面Ｇａの震度Ｉｒａの分布を、例えば時間変化を含む動画画像として視覚的に表示することができる。この表示画像により、地震動の広がりや各施設の被害状況などの推測が視覚に基づいて可能となる。操作部３４は、使用者が表示部３３の表示画面の操作を行う為に設けられている。具体的には、保守運用の操作、警報解除信号の出力、画面表示のクリア、及び、例えば推定データや実測データを含む地震情報の一覧の表示といった操作が行われる。

時刻校正部３５には、時刻を管理する高精度時計が内蔵されている。屋外に設置されているＧＰＳレシーバから１秒ごとに出力されるパルス信号によって、高精度時計の校正処理が自動的に行われる。３つの加速度計２３のサンプリング周期は、高精度時計に同期している。サンプリング周期の同期は、３つの加速度計２３に対して同時に行われる。

警報部６は、回転灯及び警報音により警報を発生する機能を有する。また、警報部６は、例えば利用者の携帯機器であってもよく、その場合、携帯機器を無線通信により入出力部３２に接続して、警報を発生する機能を有する。警報部６は、各予測期間が経過するごとに、入出力部３２から入力される警報発生信号或いは警報解除信号に従って、警報を発生するか、或いは警報の発生を解除する。警報部６は、警報発生信号が入力された場合に警報を発生し、警報解除信号が入力された場合に警報を解除する。

外部表示部７は、演算部により推定された震度等の地震情報を表示する機能を有する。外部表示部７は、上位装置３が設置されている場所とは異なる場所に設置され、地震動のＳ波が到達するまで予測データを表示する早期地震表示盤であってもよく、或いは、観測地点Ｑから離れた場所に設置され、地震情報を表示する遠隔表示盤であってもよい。

各情報端末４は、例えば主要施設（指令・保守区など）に配置される。各情報端末４は、対象施設の管理者が携帯可能な大きさを有してもよい。各情報端末４は、通信回路５を介して上位装置３と接続されており、上位装置３において推測された震度等の連絡情報を上位装置３から受け取る。各情報端末４は、入出力部４１と、表示部４２と、印刷部４３とを有する。入出力部４１は、通信回路５を介して上位装置３と通信を行うための通信インターフェイスである。また、入出力部４１は、当該情報端末４が有する表示部４２及び印刷部４３と電気的に（例えばデータバスを介して）接続されている。

表示部４２は、上位装置３の表示部３３と同様の表示を行う。すなわち、表示部４２は、複数の地域区画Ｒにおける地表面Ｇａの震度Ｉｒａの分布を、例えば時間変化を含む動画画像として視覚的に表示することができる。また、表示部４２は、例えば推定データや実測データを含む地震情報の一覧の表示を行うこともできる。一例では、表示部３３及び４２はカラー液晶ディスプレイとタッチパネルとを含む。タッチパネルは、カラー液晶ディスプレイの表示画面を視認可能な状態にて、カラー液晶ディスプレイに取り付けられる。タッチパネルは、利用者のタッチ操作を受け付ける。タッチパネルがタッチ操作されることにより、カラー液晶ディスプレイの表示画面の操作が行われる。

図８は、表示部３３及び４２における表示例を示す図である。この例では、対象施設（鉄道駅Ｔａ、鉄道線路Ｔｂおよび変電所Ｔｃ）を含む地図が示され、各対象施設の近傍に、その対象施設における推定震度を施設名称（Ａ～Ｉ及びＸ）とともに示す表示領域が設けられている。推定震度の表示領域は、推定震度の大きさに応じて色分けして示される。また、推定震度の表示領域は、タッチパネルにおける震度表示選択部Ｂ１を兼ねる。震度表示選択部Ｂ１は、画面中の部分Ｄにおいて詳細な地震情報を表示する施設を切り替えるための入力を受け付ける。

また、タッチパネルは、路線選択部Ｂ２及びＢ３、画面表示クリア部Ｂ４、地震履歴切り替え部Ｂ５、及び印刷操作部Ｂ６を含む。路線選択部Ｂ２及びＢ３は、カラー液晶ディスプレイに表示する路線を切り替えるための入力を受け付ける。画面表示クリア部Ｂ４は、画面表示をクリア（初期化）するための入力を受け付ける。地震履歴切り替え部Ｂ５は、カラー液晶ディスプレイに表示する情報として地図情報と地震履歴情報とを相互に切り替えるための入力を受け付ける。印刷操作部Ｂ６は、印刷部４３による印刷を行うための入力を受け付ける。

前述した平均Ｓ波速度データＤａには、観測点データと、緯度及び経度のデータとが含まれる。そこで、表示部３３及び４２には、対象施設を含む所定範囲における地図情報の緯度及び経度と、各単位区画Ｒａの緯度及び経度とを一致させて、推定震度を示す情報を地図に重ね合わせて表示してもよい。また、地図上において震度分布を色分けして表示してもよい。その場合、所定範囲の推定震度を視覚情報として容易に認識することができる。

印刷部４３は、本実施形態における出力部の例であり、例えば推定震度別に対象施設をグループ分けしたデータ等、推定震度に関するデータを印刷して出力することができる。図９は、印刷部４３から出力される印刷物の内容例を示す。一例として、印刷部４３は、震度情報（地震発生を最初に観測した地震計の設置場所、日時、ガル値）、推定震度、重要施設に対する規制（停止、警戒、又は規制なし）、及び区域別規制内容（重要施設と対応して示す）を印刷する。なお、各情報端末４は、表示部４２及び印刷部４３のうち一方のみを有してもよい。各情報端末４が表示部４２を有しない場合、表示部は上位装置３にのみ設けられてもよい。

通信回路５は、例えばインターネット、電話回線等の情報通信網である。通信回路５は、有線及び無線のいずれであってもよい。

図１０は、本実施形態による震度推定方法を示すフローチャートである。この震度推定方法は、例えば上述した震度推定システム１を用いて実施することができる。この震度推定方法では、まず、上位装置３において、互いに離れた複数の観測地点Ｑのそれぞれを中心とする複数の地域区画Ｒを設定する（設定ステップＳＴ１）。次に、複数の地震計２において地動を観測し、地動の地震動を検出する（検出ステップＳＴ２）。この検出ステップＳＴ２において、各地震計２は、３つの加速度計２３から得られる加速度データから、地震動のＰ波の到達を判定する。各地震計２がＰ波の到達を判定した場合（検出ステップＳＴ２：ＹＥＳ）、上位装置３は、これらの地震計２に対応する観測地点Ｑの地表面ＧａにおけるＳ波最大加速度値を推測演算する（第１の推定ステップＳＴ３）。このとき、上位装置３は、必要に応じて、観測地点Ｑの地表面ＧａにおけるＳ波最大加速度値から、観測地点Ｑの地表面Ｇａにおける震度を求める。また、上位装置３は、Ｊ－ＳＨＩＳから予め取得した、各観測地点Ｑにおける表層３０ｍの平均Ｓ波速度データ（ＡＶＳ３０）を読み出す（ステップＳＴ４）。そして、上位装置３は、第１の推定ステップＳＴ３において演算した各観測地点Ｑにおける震度（又はＳ波最大加速度値）と、ステップＳＴ４において読み出した各観測地点ＱのＡＶＳ３０とに基づいて、各観測地点Ｑの直下における工学基盤Ｇｂの震度（又はＳ波最大加速度値）を推測演算する（第２の推定ステップＳＴ５）。

続いて、上位装置３は、Ｊ－ＳＨＩＳから予め取得した、各観測地点Ｑを除く他の単位区画ＲａのＡＶＳ３０を読み出す（ステップＳＴ６）。そして、上位装置３は、各地域区画Ｒにおける工学基盤Ｇｂの震度（又はＳ波最大加速度値）が個々の地域区画Ｒ内において均一であるとの仮定の下、第２の推定ステップＳＴ５において演算した工学基盤Ｇｂの震度（又はＳ波最大加速度値）と、ステップＳＴ６において読み出した各単位区画ＲａのＡＶＳ３０とに基づいて、観測地点Ｑを除く他の単位区画Ｒａの地表面Ｇａの震度を推定する（第３の推定ステップＳＴ７）。そして、上位装置３は、推定した震度等の連絡情報を、情報端末４へ発信する（ステップＳＴ８）。また、上位装置３は、所定地域の面的震度分布を作成し、記録する（ステップＳＴ９）。上位装置３は、以上のステップＳＴ３～ＳＴ９を、所定の予測時間が経過する毎に繰り返す（ステップＳＴ１０）。

以上に説明した本実施形態による震度推定システム１及び震度推定方法によって得られる効果について説明する。地震被害をもたらす地震波には、Ｐ波とＳ波とがある。Ｐ波の伝搬速度は約６ｋｍ／秒であり、Ｓ波の伝搬速度は約３．５ｋｍ／秒である。これらの地震波をリアルタイムで捉えて強震動の予測及び情報発信を行うためには、観測網における観測点間の距離を地震波の伝搬速度と同等もしくはそれ以下とし、データ取得の時間間隔を１秒以下とすることが望ましい。また、地震被害は、主にＳ波及び表面波によりもたらされる。一方、Ｐ波の振幅とＳ波の振幅との間には相関関係が存在し、Ｐ波の振幅に基づいて、Ｓ波の振幅ひいては震度を予測することができる。Ｐ波を用いた震度予測は、震央からの距離および震源の深さが増すほど、時間的な余裕をもって揺れに備えることができるので有用である。また、地殻内にて発生するマグニチュード７程度かそれ以上の地震の場合、その震源の深さは少なくとも１０ｋｍ程度あり、Ｐ波とＳ波との到達時間差は１．７秒以上ある。加えて、破壊の成長時間による強震動の到達の遅れを考慮すると、Ｐ波の到達から強震動の到達までは３～４秒の時間がある。この時間は、回避行動のためには過小ではあるが、最も被害が大きい震源地付近に対しても強震動の到達前に予測の伝達が可能であることを意味する。例えば、鉄道の自動管制その他の稼働施設のシャットダウンのために予測情報を用いることにより、被害を軽減することが可能となる。

本実施形態の震度推定システム１及び震度推定方法によれば、地表面Ｇａの震度を各地域区画Ｒ内において面的に推定するので、観測地点Ｑの周辺における観測地点Ｑではない地点においても震度を精度良く推定することが可能であり、各地域区画Ｒに含まれる対象施設の震度を精度良く推定することができる。また、各地域区画Ｒ内の単一の観測地点Ｑにおける観測結果（Ｐ波最大加速度値）に基づいて当該地域区画Ｒ内の震度を推定するので、複数の観測地点Ｑにおける観測結果に基づいて震源位置及び地震規模（マグニチュード）を算出してから各地の震度を推定する方式と比較して、地震が発生してから震度を推定するまでの所要時間を短くすることができる。よって、被害状況の早期把握が可能となる。

また、震源位置及び地震規模を算出してから各地の震度を推定する方式では、３点以上の観測地点における観測結果に基づいて震源位置及び地震規模を算出するので、震源及びその周辺地域は、震度予測が不可能な、いわゆる“ブラインドゾーン”となる。すなわち、震源位置及び地震規模を算出するには約６ｋｍ／秒で伝播するＰ波を３つ以上の観測地点で捉える必要があるので、例えば観測地点間の平均間隔が２０ｋｍである場合、震源を中心とする半径３０ｋｍ～４０ｋｍの地域では震度予測が不可能となる。震源及びその周辺の地域は被害が集中する破壊域であるので、この地域の震度予測が不可能であることは大きな問題である。これに対し、本実施形態では、一つの観測地点における観測結果に基づいて、該観測地点の周辺の震度分布を推定するので、このようなブラインドゾーンの問題は生じない。したがって、震源付近の震度予測も可能となる。

また、従来の震度推定方式では、震度の過小評価を回避して対象施設の安全を優先するために、震度を比較的高めに推定する傾向がある。対象施設が鉄道設備等である場合、推定した震度を高めとすると運行再開判断のランクも高くなり、運転中止や保守点検を必要する範囲が広くなってしまう。これに対し、本実施形態では、対象施設の近傍に位置する観測地点ＱにおけるＰ波最大加速度値の実測データを利用して対象施設の震度を推定するので、震度をより正確に推定することができ、運転中止や保守点検を必要する範囲が無駄に広がることを抑制できる。

また、本実施形態によれば、観測地点Ｑに設置される地震計２の機能は少なくともＰ波到達判断とＰ波最大振幅値の計測とが可能であれば足りるので、システム全体のコスト削減を図ることができる。

本実施形態のように、各地域区画Ｒは、隣り合う少なくとも一つの地域区画Ｒと接するか又は縁部同士において重なり合っていてもよい。この場合、連続して延在する広い地域において震度を面的に推定できるので、例えば鉄道設備といった、長い距離にわたって連続する設備の各地点における震度を精度良く推定することができる。

本実施形態のように、震度推定システム１は、複数の地域区画Ｒにおける地表面Ｇａの震度の分布を視覚的に表示する表示部３３及び４２を備えてもよい。この場合、各地域区画Ｒ内において面的に推定された地表面Ｇａの震度を、視覚情報によって容易に且つ短時間で理解することができる。そして、表示部３３及び４２に表示された情報に基づいて、対象施設の規制や地震防災対策を迅速に実施することが可能となる。

本実施形態のように、震度推定システム１は、複数の地域区画Ｒにおける地表面Ｇａの震度の分布に関するデータを記録する記録部３２ａを備えてもよい。この場合、各地域区画Ｒ内において面的に推定された地表面Ｇａの震度に関するデータを、地震の解析に用いることができる。

本実施形態のように、震度推定システム１は、予測震度別に対象施設をグループ分けしたデータを出力する出力部（印刷部４３）を備えてもよい。この場合、推定震度が大きい（すなわち地震被害が大きいと推定される）対象施設を容易に且つ短時間で判別することができる。そして、出力されたデータに基づいて、対象施設の規制や地震防災対策を迅速に実施することが可能となる。

本実施形態のように、複数の地域区画Ｒの形状は四角形であってもよい。Ｊ－ＳＨＩＳが提供するＡＶＳ３０の区画は四角形状であるので、地域区画Ｒの形状が四角形であることによって演算を容易にできる。

本実施形態のように、複数の地域区画Ｒの面積は互いに等しくてもよい。この場合、各地域区画Ｒ内における震度の推定時間及び推定精度を複数の地域区画Ｒ間で均一に近づけることができる。

本実施形態のように、隣り合う地域区画Ｒの縁部同士が重なり合っている地域においては、隣り合う地域区画Ｒそれぞれにおいて推定された地表面Ｇａの震度のうち値が最も大きいものを用いてもよい。この場合、当該地域において地震被害を過小に推定することを回避できる。

本実施形態のように、上位装置３は（推定ステップＳＴ３，ＳＴ５，ＳＴ７では）、各地震計２から周期的に受けるＰ波最大加速度値から各観測地点ＱにおけるＳ波最大加速度値を推測演算し、Ｓ波最大加速度値から各観測地点Ｑにおける震度の推定を行い、該震度に基づいて各観測地点Ｑにおける工学基盤Ｇｂの震度を推測演算してもよい。この場合、Ｐ波最大加速度値に基づく工学基盤Ｇｂの震度の推測演算を精度良く行うことができる。

本実施形態のように、対象施設は鉄道設備（鉄道駅Ｔａ、鉄道線路Ｔｂおよび変電所Ｔｃ）であってもよい。上述したように、本実施形態の震度推定システム１及び震度推定方法によれば、鉄道設備のような広い範囲に延在する設備であっても、対象施設の震度を精度良く推定することができる。したがって、地震が発生した際、鉄道の運行を停止するか否かの正確な決定に寄与でき、また、運行再開判断の支援を行うことができる。なお、対象施設は鉄道設備に限られず、例えば鉄道以外の建造物や道路といった様々な施設を対象施設とすることができる。

本発明による震度推定システム及び震度推定方法は、上述した実施形態に限られるものではなく、他に様々な変形が可能である。例えば、上記実施形態では、Ｐ波最大加速度値に基づく震度又はＳ波最大加速度値の推定を上位装置３が行っているが、震度又はＳ波最大加速度値を推定するための各演算の少なくとも一部を各地震計２が行ってもよい。言い換えると、上位装置３の演算部３１の機能の少なくとも一部を各地震計２が担ってもよい。

１…震度推定システム、２…地震計、３…上位装置、４…情報端末、５…通信回路、６…警報部、７…外部表示部、２１…観測部、２２…処理部、２３…加速度計、２４…計測基板、２５…観測データ処理部、３１…演算部、３２…入出力部、３２ａ…記録部、３３…表示部、３４…操作部、３５…時刻校正部、４１…入出力部、４２…表示部、４３…印刷部、Ｂ１…震度表示選択部、Ｂ２…路線選択部、Ｂ３…路線選択部、Ｂ４…画面表示クリア部、Ｂ５…地震履歴切り替え部、Ｂ６…印刷操作部、Ｄａ…平均Ｓ波速度データ、Ｇａ…地表面、Ｇｂ…工学基盤、Ｌ…線、Ｑ…観測地点、Ｑｂ…観測地点の直下の地点、Ｒ…地域区画、Ｒａ…単位区画、Ｔ_１～Ｔ_６…予測期間、Ｔ_Ｐ，Ｔ_Ｓ…時刻、Ｔａ，Ｔａ１，Ｔａ２…鉄道駅、Ｔｂ…鉄道線路、Ｔｃ…変電所。

Claims (9)

1. 互いに離れた複数の観測地点のそれぞれにおいて地動を観測し、前記地動の地震動を検出する複数の地震計と、
   Ｐ波最大加速度値を含む前記地震動に関する情報を前記複数の地震計から受ける上位装置と、
   を備え、
   前記上位装置において、各観測地点を中心とする複数の地域区画が予め設定されており、
   前記複数の地域区画のうち二以上の前記地域区画のそれぞれには、震度を予測すべき対象施設が含まれており、
   前記上位装置は、各観測地点にＳ波が到達するまでの間、前記Ｐ波最大加速度値を各地震計から周期的に受け、前記Ｐ波最大加速度値に基づいて各観測地点における工学基盤の震度又はＳ波最大加速度値を推測演算し、各地域区画における工学基盤の震度又はＳ波最大加速度値が個々の前記地域区画内において均一であるとの仮定の下、各地域区画内における表層地盤増幅率の分布に基づいて、地表面の震度を各地域区画内において面的に推定し、
   各地域区画は、隣り合う少なくとも一つの前記地域区画と接するか又は縁部同士において重なり合っており、
   前記複数の観測地点を一本の線にて結んだとき、その線上において互いに隣り合う観測地点間の距離は全て等しく、
   前記複数の地域区画における地表面の震度の分布を視覚的に表示する一又は複数の表示部を更に備える、震度推定システム。
2. 前記複数の地域区画における地表面の震度の分布に関するデータを記録する記録部を更に備える、請求項１に記載の震度推定システム。
3. 予測震度別に前記対象施設をグループ分けしたデータを出力する出力部を更に備える、請求項１または２に記載の震度推定システム。
4. 前記複数の地域区画の形状は四角形である、請求項１～３のいずれか１項に記載の震度推定システム。
5. 前記複数の地域区画の面積は互いに等しい、請求項１～４のいずれか１項に記載の震度推定システム。
6. 隣り合う前記地域区画の縁部同士が重なり合っている地域においては、前記隣り合う地域区画それぞれにおいて推定された前記地表面の震度のうち値が最も大きいものを用いる、請求項１～５のいずれか１項に記載の震度推定システム。
7. 前記上位装置は、各地震計から周期的に受ける前記Ｐ波最大加速度値から各観測地点におけるＳ波最大加速度値を推測演算し、前記Ｓ波最大加速度値から各観測地点における震度の推定を行い、該震度に基づいて各観測地点における工学基盤の震度を推測演算する、請求項１～６のいずれか１項に記載の震度推定システム。
8. 前記対象施設は鉄道設備である、請求項１～７のいずれか１項に記載の震度推定システム。
9. 互いに離れた複数の観測地点のそれぞれを中心とする複数の地域区画を予め設定する設定ステップと、
   各観測地点において地動を観測し、前記地動の地震動を検出する検出ステップと、
   各観測地点にＳ波が到達するまでの間、各観測地点のＰ波最大加速度値に基づいて各観測地点における工学基盤の震度又はＳ波最大加速度値を推測演算し、各地域区画における工学基盤の震度又はＳ波最大加速度値が個々の前記地域区画内において均一であるとの仮定の下、各地域区画内における表層地盤増幅率の分布に基づいて、地表面の震度を各地域区画内において面的に推定する推定ステップと、
   前記複数の地域区画における地表面の震度の分布を視覚的に表示するステップと、
   を含み、
   各地域区画は、隣り合う少なくとも一つの前記地域区画と接するか又は縁部同士において重なり合っており、
   前記複数の地域区画のうち二以上の前記地域区画のそれぞれには、震度を予測すべき対象施設が含まれており、
   前記複数の観測地点を一本の線にて結んだとき、その線上において互いに隣り合う観測地点間の距離は全て等しい、震度推定方法。

Images