特開 2018-54339 震源情報可視化システム及び震源情報可視化方法

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

(12)【公報種別】公開特許公報(A)

(11)【公開番号】特開2018-54339(P2018-54339A)

(43)【公開日】2018年4月5日

(54)【発明の名称】震源情報可視化システム及び震源情報可視化方法

(51)【国際特許分類】

   G01V 1/34 20060101AFI20180309BHJP

   G06T 17/05 20110101ALI20180309BHJP

   G06T 11/20 20060101ALI20180309BHJP

【ＦＩ】

   G01V1/34

   G06T17/05

   G06T11/20 600

【審査請求】未請求

【請求項の数】7

【出願形態】ＯＬ

【全頁数】18

(21)【出願番号】特願2016-187518(P2016-187518)

(22)【出願日】2016年9月26日

(71)【出願人】

【識別番号】504137912

【氏名又は名称】国立大学法人 東京大学

(71)【出願人】

【識別番号】515000742

【氏名又は名称】株式会社ライブ・アース

(74)【代理人】

【識別番号】110000925

【氏名又は名称】特許業務法人信友国際特許事務所

(72)【発明者】

【氏名】河合 研志

(72)【発明者】

【氏名】庄司 真史

(72)【発明者】

【氏名】小林 佑介

【テーマコード（参考）】

2G105

5B050

5B080

【Ｆターム（参考）】

2G105AA03

2G105MM00

5B050AA06

5B050BA06

5B050BA18

5B050CA07

5B050EA09

5B050EA12

5B050EA27

5B050EA28

5B050EA30

5B050FA02

5B050FA05

5B050FA08

5B050GA08

5B080AA19

5B080BA08

5B080FA02

5B080GA00

(57)【要約】

【課題】地震の位置及び規模の情報から、地震発生の状況の正確な把握ができると共に、地下の断層面の適切な推定が行えるようにする。
【解決手段】一定の期間内に発生した地震について、規模、震源位置、発生日時及びその震源位置の地表面からの深さの情報を、震源情報として取得する。そして、取得した震源情報から、指定された座標範囲の震源情報を読出し、それぞれの震源情報の震源位置と深さに対応した座標位置を、３次元空間内にプロットした３次元震源分布画像を生成する。その３次元震源分布画像の表示などで、利用者に提示する。震源情報を選ぶ期間や座標位置の範囲などを選択することで、様々な条件での地震の発生状況を表示できるようになる。
【選択図】図１０

【特許請求の範囲】

【請求項1】

地震情報のデータベースから地震情報を収集する地震情報収集部と、
前記地震情報収集部が収集した地震情報を記憶する地震情報蓄積部と、
前記地震情報蓄積部に記憶された地震情報を検索する地震情報検索部と、
前記地震情報検索部に地震情報の検索を指示し、その検索結果を震源情報として取得する地震情報取得部と、
前記地震情報取得部が取得する震源情報の取得条件を設定する読出し範囲設定部と、
前記地震情報取得部が取得した震源情報に含まれるそれぞれの震源の位置と深さに対応した座標位置を、３次元空間の内部にプロットした３次元震源分布画像を生成する３次元画像生成部と、
前記３次元画像生成部が生成する３次元震源分布画像の表示条件を設定する表示条件設定部と、
前記３次元画像生成部が生成した前記３次元震源分布画像を提示する提示部とを備えた
震源情報可視化システム。

【請求項2】

前記３次元画像生成部は、前記３次元空間として、前記表示条件設定部により設定された深さの地表面と最深面との間に形成される３次元空間とし、前記３次元震源分布画像内に、地表面及び最深面を表示するようにした
請求項１に記載の震源情報可視化システム。

【請求項3】

前記震源情報から得たそれぞれの地震の規模に応じて、前記３次元震源分布画像内の各震源位置を示す印のサイズ又は表示形態を変更するようにした
請求項１又は２に記載の震源情報可視化システム。

【請求項4】

前記震源情報から得たそれぞれの地震の震源位置の深さに応じて、前記３次元震源分布画像内の各震源位置を示す印の表示形態を変更するようにした
請求項１又は２に記載の震源情報可視化システム。

【請求項5】

前記震源情報から得たそれぞれの地震の発生日又は時間に応じて、前記３次元震源分布画像内の各震源位置を示す印の表示形態を変更するようにした
請求項１又は２に記載の震源情報可視化システム。

【請求項6】

前記３次元画像生成部は、３次元空間内にプロットした複数の震源位置から求めた近似平面を、３次元震源分布画像に加えるようにした
請求項１〜５のいずれか１項に記載の震源情報可視化システム。

【請求項7】

地震情報のデータベースから地震情報を収集する地震情報収集処理と、
前記地震情報収集処理が収集した地震情報を記憶する地震情報蓄積処理と、
前記地震情報蓄積処理により記憶された地震情報を検索する地震情報検索処理と、
前記地震情報検索処理に地震情報の検索を指示し、その検索結果を震源情報として取得する地震情報取得処理と、
前記地震情報取得処理により取得した震源情報が取得する震源情報の取得条件を設定する読出し範囲設定処理と、
前記地震情報取得処理により取得した震源情報に含まれるそれぞれの震源の位置と深さに対応した座標位置を、３次元空間内にプロットした３次元震源分布画像を生成する３次元画像生成処理と、
前記３次元画像生成処理により生成する３次元震源分布画像の表示条件を設定する表示条件設定処理と、
前記３次元画像生成処理により生成した前記３次元震源分布画像を提示する提示処理とを含む
震源情報可視化方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、地震の震源に関する情報を可視化する震源情報可視化システム及び震源情報可視化方法に関する。

【背景技術】

【0002】

地震は、地下の断層で起きる破壊を伴うすべり運動（破壊すべり）である。このすべり運動による断層が地表面で生じた場合は、比較的容易に特定することができる。地上に現れる断層面は、長年の地質踏査による活断層マッピング等により広く知られたものとなっており、防災分野で活用されている。
特許文献１には、航空機で撮影した地上の画像から、地表面の断層を検出する手法が記載されている。

【先行技術文献】

【特許文献】

【0003】

【特許文献1】特開２０１５−１１１２４８号公報

【発明の概要】

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

【0004】

ところが、地下に存在する断層面は、地質踏査では特定するのが困難であり、地球物理学的手法が必要であった。例えば、人工的に起こした地震波などの断層面での反射波を用いて地下断層面を特定する手法や、地図上での余震の分布から本震のすべり面を特定する手法などが知られている。反射波を用いて地下断層面を特定するためには、大規模な物理探査が必要であり、簡単に地下断層面を特定することはできないという問題があった。

【0005】

また、地図上での余震の分布から本震のすべり面を特定する場合には、規模や地質情報を用いて予め余震が発生したと想定される期間を定めて、その期間内に設定した地域で発生した地震の情報を使って、１つの地下断層面を特定することが行われる。ここで、時間的にどこまで後の地震を余震とするか、どこまでの距離の地震がすべり面の影響によるものであるか、などの条件の設定が難しいという問題があった。
例えば、本震と思われる地震から特性時間（たとえば、１週間後）までに発生した地震の情報からすべり面を特定した場合と、本震と思われる地震から特性時間の２倍（２週間後）までに発生した地震の情報からすべり面を特定した場合とでは、すべり面の位置が相違する可能性がある。従来は、解析を行う者が、解析に使用する地震情報の時間的な範囲や地域的な範囲を決めて、その条件ですべり面を求めるようにしていた。しかしながら、実際にはどの地震が本当に余震であるのかを判断することは非常に難しく、決めた条件が本当に正しいかどうかの検証は、容易にはできないという問題があった。
また、地震の発生状況によっては、複数のすべり面の影響で地震が発生することもあるが、そのような複数のすべり面が影響しているのか否かを、簡単に検証することも困難であった。

【0006】

本発明の目的は、地震の位置、規模及び発生日時の情報から、視覚的に地震発生の状況の正確な把握ができると共に、地下の断層面の適切な推定が行えるようにすることにある。

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

【0007】

本発明の震源情報可視化システムは、地震情報収集部と、地震情報蓄積部と、地震情報検索部と、地震情報取得部と、読出し範囲設定部と、３次元画像生成部と、表示条件設定部と、提示部とを備える。
地震情報収集部は、地震情報のデータベースから、定期的に地震情報を収集する。地震情報のデータベースは、例えばインターネット上で公開されているデータベースを利用する。利用者が作成した、または指定したデータベースを追加で入力する場合も、ここで収集する。
地震情報蓄積部は、地震情報収集部が収集した地震情報を記憶する。
地震情報検索部は、地震情報取得部から指定される条件に従って、地震情報蓄積部に記憶された地震情報を検索する。
地震情報取得部は、読出し範囲設定部により指定される条件に従って、地震情報検索部に地震情報の検索を指示し、その検索結果を震源情報として取得する。
読出し範囲設定部は、地震の発生日時、地震の震源の座標位置、震源の深さ、及び地震の規模といった、地震情報取得部が取得する震源情報の取得条件を設定する。
３次元画像生成部は、地震情報取得部が取得した震源情報に含まれるそれぞれの震源の位置と深さに対応した座標位置を、３次元空間内にプロットした３次元震源分布画像を生成する。
表示条件設定部は、３次元画像生成部が生成する３次元震源分布画像の表示条件を設定する。
提示部は、３次元画像生成部が生成した３次元震源分布画像を提示する。

【0008】

また本発明の震源情報可視化方法は、以下の(a)〜(h)の処理を行う工程を含む。
(a)インターネット上で公開されている地震情報のデータベースから、定期的に地震情報を収集する地震情報収集処理および利用者指定データベースの地震情報収集処理
(b)地震情報収集処理が収集した地震情報を記憶する地震情報蓄積処理
(c)地震情報取得処理から指定される条件に従って、地震情報蓄積処理により記憶された地震情報を検索する地震情報検索処理
(d)読出し範囲設定処理により指定される条件に従って、地震情報検索処理に地震情報の検索を指示し、その検索結果を震源情報として取得する地震情報取得処理
(e)地震の発生日時、地震の震源の座標位置、震源の深さ、及び地震の規模といった、地震情報取得処理が取得する震源情報の取得条件を設定する読出し範囲設定処理
(f)地震情報取得処理により取得した震源情報に含まれるそれぞれの震源の位置と深さに対応した座標位置を、３次元空間内にプロットした３次元震源分布画像を生成する３次元画像生成処理
(g)３次元画像生成処理により生成する３次元震源分布画像の表示条件を設定する表示条件設定処理
(h)３次元画像生成処理により生成した３次元震源分布画像を提示する提示処理

【発明の効果】

【0009】

本発明によると、本震とその本震の前後に発生した地震を３次元空間内にプロットした３次元震源分布画像が得られ、地震発生位置の付近の地下での地震の発生状況を３次元震源分布画像から的確に評価することができる。例えば、３次元震源分布画像内での本震と余震のプロット位置から、本震を発生した地下断層面を的確かつ簡単に推定できるようになる。また、画像内に表示する震源情報の発生期間や座標位置の範囲などの条件を変更することで、異なった条件での地震の発生状況を簡単に表示できるため、様々な条件での地震の評価が簡単にできるようになる。

【図面の簡単な説明】

【0010】

【図1】本発明の一実施の形態例によるシステム構成例を示すブロック図である。

【図2】本発明の一実施の形態例のシステムに適用されるハードウェア構成例を示すブロック図である。

【図3】本発明の一実施の形態例による３次元震源分布画像の生成処理例を示すフローチャートである。

【図4】本発明の一実施の形態例による設定画面の例を示す説明図である。

【図5】本発明の一実施の形態例による表示形態の設定例を示す説明図である。

【図6】本発明の一実施の形態例による３次元震源分布画像の例（例１）を示す説明図である。

【図7】本発明の一実施の形態例による３次元震源分布画像の例（例２）を示す説明図である。

【図8】本発明の一実施の形態例による３次元震源分布画像の例（例３）を示す説明図である。

【図9】本発明の一実施の形態例による３次元震源分布画像の例（例４）を示す説明図である。

【図10】本発明の一実施の形態例による３次元震源分布画像の例（例５）を示す説明図である。

【図11】本発明の一実施の形態例による３次元震源分布画像の例（例６）を示す説明図である。

【図12】本発明の一実施の形態例による３次元震源分布画像の例（例７）を示す説明図である。

【図13】本発明の一実施の形態例による３次元震源分布画像の例（例８）を示す説明図である。

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

【0011】

以下、本発明の一実施の形態例（以下、本例と称する）を、添付図面を参照して説明する。
［１．システム構成例］
図１は、本例の震源情報可視化システム１００の機能的な構成を示す機能ブロック図である。
震源情報可視化システム１００は、サーバー装置１１０とクライアント装置１２０とで構成される。サーバー装置１１０は、地震情報収集部１１１を備え、ネットワークＮを経由して、定期的に、例えば１日１回、地震情報データベース１０にアクセスして、最新の地震情報を取得する地震情報収集処理を行う。定期的に最新情報を取得するのは一例であり、一括してデータを収集してもよい。
地震情報データベース１０としては、例えば気象庁が過去の地震情報を蓄積して、インターネット上で公開しているデータベースを利用する。また、地震情報収集部１１１は、利用者指定の特定のデータベースを使用したり、あるいは利用者自身が入力したデータセットをデータベースとして用意して、利用してもよい。

【0012】

地震情報収集部１１１は、取得した地震情報を、地震情報蓄積部１１２に記憶する地震情報蓄積処理を行う。ここで地震情報収集部１１１が地震情報データベース１０より取得（収集）する地震情報には、地震が発生した日時、地震の震源の座標位置、震源の深さ、及び地震の規模（マグニチュード）の情報が含まれる。なお、震源の深さは、例えば平均海水面（標高０ｍ）からの深さを示す。
地震情報収集部１１１は、取得したそれぞれの地震情報を、地震情報蓄積部１１２に記憶する。本例の震源情報可視化システム１００では、この地震情報蓄積部１１２に記憶された地震情報で示される震源を表示して利用者に提示する処理が行われる。

【0013】

クライアント装置１２０の地震情報取得部１２１は、地震情報検索部１１３に対し、利用者により読出し範囲設定部１２２に入力された条件に合致する地震情報の取得を要求する。読出し範囲設定部１２２に入力される条件としては、例えば後述する設定画面（図４）で示されるような、地震が発生した日付（又は時間）の範囲、震源位置の緯度及び経度の範囲、震源の深さの範囲、地震の規模の範囲などがある。

【0014】

地震情報検索部１１３は、地震情報蓄積部１１２に蓄積された地震情報のうち、指定された条件に合致する地震情報を検索して抽出する地震情報検索処理を行い、地震情報取得部１２１に返却する。地震情報取得部１２１は、前項にて指定された地震情報の検索条件に従って地震情報検索部１１３により抽出された地震情報を取得する地震情報取得処理を行う。

【0015】

地震情報取得部１２１が地震情報検索部１１３より取得した地震情報は、３次元画像生成部１２３に渡される。さらに、３次元画像生成部１２３は、その仮想的に生成された３次元空間の内部に、地震情報取得部１２１より渡された地震情報の各震源位置をプロットした３次元震源分布画像を生成する。震源位置をプロットする３次元空間は、震源の深さの範囲の条件で設定した地表面及び最深面を、上面及び下面とした３次元空間であり、３次元震源分布画像では、この地表面と最深面が表示される。このとき、地表面には、該当する緯度及び経度の範囲内の地図又は地形が表示される。
なお、表示モードの設定によっては、３次元震源分布画像を表示する際に、距離や深さなどの値、緯度値、経度値を表示するようにしてもよい（後述する図６，図７参照）。
また、ここでは上面を地表面と称するが、後述する設定画面は、深さの表示範囲として、上面を標高０ｍ以外の面としたときに、画像内に表示される地表面が標高０ｍ以外の面となる場合もある。

【0016】

３次元画像生成部１２３が３次元震源分布画像を生成する際には、３次元震源分布画像をどのアングルから見た画像とするか等の表示条件が、表示条件設定部１２４から指示される。この表示アングルなどの指示は、例えばコンピューター装置の場合にはマウスなどの操作に基づいて行われ、タブレット端末やスマートフォンなどのタッチパネルを備えた装置の場合には、画像を表示した画面上のタッチ操作などに基づいて行われる。
震源位置を示す印の表示形態についても、表示条件設定部１２４からの指示で設定される。本例の場合には、３次元画像生成部１２３が３次元空間の内部に震源位置をプロットする場合、例えばそれぞれの震源を球体（円形）の印で示す。この場合、円形の球体の印のサイズや表示色が、表示条件設定部１２４で設定した条件により決まる。例えば、地震のマグニチュードが大きいほど、大きい球体の印とする。あるいは、地震が発生した日付や期間ごとに異なる色の球体の印とする。あるいはまた、震源位置の深さによって、異なる色の球体の印とする。

【0017】

さらに、３次元震源分布画像を作成する際には、３次元空間の内部の震源の位置から求まる近似平面を表示することもできる。ここでの近似平面とは、３次元空間に配置した各震源の座標位置からの距離が最小となる平面である。すなわち、それぞれの震源の座標位置が、３軸ｘ，ｙ，ｚ上の点列(x₀,y₀,z₀)、(x₁,y₁,z₁)、・・・、(x_n,y_n,z_n)で示されるとき、その点列(x₀,y₀,z₀)、(x₁,y₁,z₁)、・・・、(x_n,y_n,z_n)を最小二乗近似する平面ｚ＝ａｘ＋ｂｙ＋ｃ（ここで、ａ，ｂ，ｃは係数）が近似平面となる。

【0018】

この近似平面を求める処理をより詳細に説明すると、点列(x₀,y₀,z₀),(x₁,y₁,z₁),・・・,(x_n,y_n,z_n)について、点(x_i,y_i,z_i)と平面ｚ＝ａｘ＋ｂｙ＋ｃのｚ方向の距離ｄ_ｉの２乗和Ｑ＝Σd_i²＝Σ(ax_i+by_i+c-z_i)²を最小にする係数(a,b,c)を決定し代入した平面ｚが、近似平面である。
この近似平面を求めることは、地震により発生したすべり面（断層面）を推定することに相当する。これらの表示を行う条件についても、設定画面から利用者が設定することが可能である。
さらにまた、３次元震源分布画像内の特定の震源位置の印（球体）が、利用者による画面のクリックなどで選択されたときには、その座標位置の震源についての詳細情報（マグニチュードの値、発生日時、緯度、経度、深さなど）を文字や数字で直接表示することができる。
これらの表示形態の具体な例については、図６以降の表示例で後述する。

【0019】

３次元画像生成部１２３が生成した３次元震源分布画像は、表示部１２５により表示される。この表示部１２５により、利用者に３次元震源分布画像が提示される。また、３次元画像生成部１２３が生成した３次元震源分布画像は、出力部１２６から外部に出力するようにしてもよい。例えば、震源情報可視化システム１００に外部の端末からのアクセスで３次元震源分布画像の作成指示がある場合には、出力部１２６が、外部の端末に対して３次元震源分布画像を送信して、外部の端末の利用者に３次元震源分布画像を提示する。また、出力部１２６として、プリンターを用意し、３次元震源分布画像の印刷物を出力するようにしてもよい。

【0020】

［２．震源情報可視化システムのハードウェア構成例］
図２は、本例の震源情報可視化システム１００に適用されるコンピューター装置のハードウェア構成例を示す。
図１に示す震源情報可視化システム１００のサーバー装置１１０及びクライアント装置１２０は、それぞれコンピューター装置Ｃにより構成される。クライアント装置１２０を構成するコンピューター装置としては、スマートフォンやタブレット端末などの、様々な形態の情報処理装置が含まれる。

【0021】

図２に示すコンピューター装置Ｃは、バスラインＣ４にそれぞれ接続されたＣＰＵ（Central Processing Unit：中央処理装置）Ｃ１、ＲＯＭ（Read Only Memory）Ｃ２、及びＲＡＭ（Random Access Memory）Ｃ３を備える。さらに、コンピューター装置Ｃは、不揮発性ストレージＣ５、入出力部Ｃ６、及びディスプレイＣ７を備える。

【0022】

ＣＰＵ Ｃ１は、本例の震源情報可視化システム１００が備える各機能を実現するソフトウェアのプログラムコードをＲＯＭ Ｃ２から読み出して実行する。ＲＡＭ Ｃ３には、演算処理の途中に発生した変数やパラメータ等が一時的に書き込まれる。本例の場合、ＣＰＵ Ｃ１がＲＯＭ Ｃ２に記憶されているプログラムを読み出すことで、３次元震源分布画像の生成処理が実行される。

【0023】

不揮発性ストレージＣ５としては、例えば、ＨＤＤ（Hard disk drive）、ＳＳＤ（Solid State Drive）、フレキシブルディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＣＤ−Ｒ、磁気テープ、不揮発性のメモリ等が用いられる。この不揮発性ストレージＣ５には、ＯＳ（Operating System）、各種のパラメータの他に、コンピューター装置Ｃを震源情報可視化システム１００として機能させるためのプログラムが記録されている。

【0024】

入出力部Ｃ６は、例えばネットワークに接続するためのインターフェイス装置などが用いられ、各種のデータを外部との間で送受信することが可能である。本例の場合、入出力部Ｃ６は、３次元震源分布画像のデータを外部の端末に送信する。あるいは、プリンターに印刷用の画像データを出力するようにしてもよい。
ディスプレイＣ７は、コンピューター装置Ｃにより実行した処理結果を表示する。ここでは、ＣＰＵ Ｃ１の制御で作成された３次元震源分布画像がディスプレイＣ７に表示される。

【0025】

［３．３次元震源分布画像を作成する処理例］
図３は、本例の震源情報可視化システム１００が３次元震源分布画像を作成する際の処理の流れを示すフローチャートである。
まず、震源情報可視化システム１００の読出し範囲設定部１２２が、３次元震源分布画像を作成する際に使用する地震情報を読出す条件として、表示日時範囲を設定する（ステップＳ１１）。
次に、読出し範囲設定部１２２は、３次元震源分布画像として震源を表示する範囲（座標位置）を設定する読み出し範囲設定処理を行う（ステップＳ１２）。

【0026】

さらに、読出し範囲設定部１２２は、震源の深さの表示範囲の設定を行うとともに（ステップＳ１３）、地震のマグニチュードの表示範囲の設定を行う（ステップＳ１４）。これらのステップＳ１１〜Ｓ１４での表示条件設定処理は、設定画面上での利用者による入力に基づいて行われる。
そして、表示条件設定部１２４は、表示色閾値の設定（ステップＳ１５）、及びグラデーションの設定（ステップＳ１６）を、それぞれ設定画面上での利用者による入力に基づいて行う。
次に、３次元画像生成部１２３は、読出し範囲設定部１２２により設定された条件の地震情報を地震情報蓄積部１１２から読出す。そして、表示条件設定部１２４は、設定された条件を表示形態とする指示を３次元画像生成部１２３に対して行う。なお、ステップＳ１１からステップＳ１６の処理については、上述した順序に限らない。つまり、いずれの順序で処理してもよく、図３のフローチャートに示す順序はあくまでも一例を示したに過ぎない。

【0027】

ステップＳ１５における表示色閾値の設定は、マグニチュードの表示色を変更する閾値を設定するものである。この表示色閾値を設定する際には、設定した閾値以上のマグニチュードの地震について、他の地震と区別した特別な色（例えば目立つ色）で表示する処理が行われる。例えば、この表示色閾値の機能を使って、本震よりも若干低いマグニチュードを表示色閾値に設定することで、３次元震源分布画像で、本震だけを特別な色で表示することができる。

【0028】

また、ステップＳ１６におけるグラデーションの設定は、複数の地震を色分け表示する条件を設定するものである。例えば、マグニチュードによる色分け、日付による色分け、及び深さによる色分けの３つのモードが用意され、設定画面において利用者により表示されるモードの設定が行われる。

【0029】

次に、表示条件設定部１２４は、震源位置の近似平面を表示する設定がなされているか否かを判断する（ステップＳ１７）。このステップＳ１７で、近似平面を表示する設定がなされている場合には（ステップＳ１７のＹＥＳ）、表示条件設定部１２４は、３次元震源分布画像として近似平面を表示する設定を行う（ステップＳ１８）。なお、この近似平面は、震源位置の集合を最もよく説明した平面であり、３次元空間の内部にプロットした各震源位置からの距離が最小となるように設定した平面である。また、このステップＳ１７からステップ１８の処理については、前述のステップＳ１１からステップＳ１６の各処理の間のどこかにおいて行われてもよく、図３のフローチャートに示す順序はあくまでも一例を示したに過ぎない。

【0030】

ステップＳ１７で近似平面を表示する設定がなされていない場合（ステップＳ１７のＮＯ）、及びステップＳ１８での近似平面の設定が行われた場合には、地震情報取得部１２１が設定された条件で地震情報の検索指示を地震情報検索部１１３に対して行い、地震情報検索部１１３は指示された条件に基づいて地震情報蓄積部１１２より地震情報を検索する（ステップＳ１９）。そして、その検索された地震情報を使って、３次元画像生成部１２３は、３次元空間の内部にそれぞれの地震情報をプロットする３次元画像生成処理を行う（ステップＳ２０）。また、３次元画像生成部１２３は、表示条件設定部１２４から指示されたアングルで３次元空間の内部に地震情報がプロットされた３次元震源分布画像の描画処理を行う（ステップＳ２１）。なお、３次元画像生成部１２３が３次元震源分布画像を描画する際には、必要により距離，深さ，緯度，経度などの数値を画像中に付加する。また、いずれかの震源の詳細表示を行う指示がある場合には、３次元画像生成部１２３は、該当する震源の詳細情報（マグニチュードの値、発生日時、緯度、経度、深さなど）を文字や数字で加える。
このようにして描画された３次元震源分布画像は、表示部１２５に表示されるとともに、出力部１２６からプリントなどの出力処理が行われ、利用者に提示される。

【0031】

［４．設定画面の例］
図４は、３次元震源分布画像を生成する際の各項目の設定画面の例である。
この設定画面は、表示部１２５に表示される画面であり、画面内の各設定項目の入力は、利用者の操作により行われる。
図４に示す設定画面には、上から順に、日付の範囲の入力箇所、緯度の範囲の入力箇所、経度の範囲の入力箇所、深さの範囲の入力箇所、マグニチュードの入力箇所、表示色閾値の値の入力箇所、近似平面の表示の有無の入力箇所、及びグラデーションのモードの選択箇所が設けられている。

【0032】

日付の範囲の入力箇所は、３次元震源分布画像内にプロットされる地震の発生日の範囲を指示する。図４では、２０１６年４月１４日から２０１６年４月２０日を範囲として設定した例を示す。
緯度及び経度の範囲の入力箇所は、３次元震源分布画像内にプロットされ地震の緯度及び経度の範囲を指示する。図４では、緯度（北緯）として３１°から３５°の範囲、経度（東経）として１２９°から１３３°の範囲を設定した例を示す。

【0033】

深さの範囲の入力箇所は、３次元震源分布画像内にプロットされる地震の地表面から最深面までの範囲を指示する。図４では、地表面０ｋｍから最深面２０ｋｍまでの範囲を設定した例を示す。
マグニチュードの範囲の入力箇所は、三次元空間の内部にプロットされる地震のマグニチュードの範囲を指示する。図４では、マグニチュードの範囲の指示を行わない例を示し、この場合には、地震情報として得られた全てのマグニチュードの地震が表示される。

【0034】

表示色閾値の値の入力箇所には、閾値を超えたときに特定の色（表示形態）で表示する際の、その閾値のマグニチュードの値が入力される。図４では、マグニチュード６以上を設定した例を示す。図６以降の表示例では、閾値を超えた地震の印は、網掛け表示としてある。
近似平面の表示の有無の入力箇所は、近似平面を表示する指示があるとき、チェック有りになり、近似平面を表示する指示がないときには、チェック無しになる。
グラデーションのモードの選択箇所は、日付による色分けモード、マグニチュードによる色分けモード、深さによる色分けモードのいずれか１つのモードでの表示が選択される。図４では、日付による色分けモードを設定した例を示している。
図４には示さないが、設定画面において、距離，深さ，緯度，経度などの数値の表示を行うか否かの選択や、地図の表示を行うか否かの選択、地図の種類の選択、緯度経度深さのスケールを合わせるか否かの選択等を行うようにしてもよい。

【0035】

［５．３次元震源分布画像の表示形態の設定例］
３次元画像生成部１２３は、表示条件設定部１２４からの指示に基づいて、３次元空間を見るアングルを任意の方向に設定することにより、３次元震源分布画像を生成することができる。
図５Ａ，Ｂ，Ｃは、それぞれ別のアングルで３次元空間を見た場合の例を示す。なお、図５では、３次元空間の表示形態を示し、震源を示す印はプロットされていない。
図５Ａは、３次元空間の地表面Ｐ１を真上から見たアングルである。地表面Ｐ１には、該当する範囲の地図ｍが表示されている。また、南Ｓと北Ｎの方向を示すラインと、東Ｅと西Ｗの方向を示すラインが、地表面Ｐ１内に表示される。なお、図５Ａでは、地表面Ｐ１は最深面と重なっている。

【0036】

この図５Ａに示す状態から、表示条件設定部１２４が、地表面Ｐ１を水平方向（矢印θ１の方向）に回転させる指示を行うと、３次元画像生成部１２３はその指示に対応した回転位置の画像を描画する。また、表示条件設定部１２４が、東Ｅと西Ｗの方向を示すラインを軸として、その軸の周囲を矢印θ２の方向に回転させる指示を行うと、３次元画像生成部１２３はその指示に対応したアングルの画像を描画する。図５Ｂ及び図５Ｃは、図５Ａに示す状態から、矢印θ２の方向の回転指示を行ったアングルでの表示形態を示している。

【0037】

図５Ｂは、３次元空間を斜め上から見たアングルの例を示す。
この状態では、地表面Ｐ１と最深面Ｐ２とが、上下に離れてほぼ同じ形状で表示され、その地表面Ｐ１と最深面Ｐ２との間に、地震情報で示される震源位置がプロットされる。地表面Ｐ１には、南Ｓと北Ｎの方向を示すライン、東Ｅと西Ｗの方向を示すライン、及び地図ｍが表示される。
また、地表面Ｐ１の中心位置と最深面Ｐ２の中心位置を貫通する鉛直軸Ｚが表示される。
なお、図５Ｂは、近似平面Ｐ３の一例を示した図である。この近似平面Ｐ３は、画像中に表示された全ての震源の座標位置からの距離が最小となる位置とした、震源の集合を最も良く表した平面であり、図５Ｂに示す表示位置はあくまでも一例である。

【0038】

図５Ｃは、３次元空間を横から見たアングルの例を示す。
この状態では、地表面Ｐ１と最深面Ｐ２とが、相互に反転した台形形状で表示され、その地表面Ｐ１と最深面Ｐ２との間に、地震情報で示される震源位置がプロットされる。地表面Ｐ１には、南Ｓと北Ｎの方向を示すライン、東Ｅと西Ｗの方向を示すライン、鉛直軸Ｚ、及び地図ｍが表示される点は、図５Ｂのアングルと同じである。この図５Ｃの例も、近似平面Ｐ３の一例である。

【0039】

［６．実際の３次元震源分布画像の例］
図６〜図１３は、実際の震源分布を示す３次元震源分布画像の例を示した図である。図６〜図１３の３次元震源分布画像２００ａ〜２００ｈでは、ライン２０１は、図５に示した東Ｅと西Ｗの方向を示すラインであり、ライン２０２は、図５に示した南Ｓと北Ｎの方向を示すラインである。

【0040】

なお、図６及び図７は、距離，深さ，緯度，経度の値を表示させた例を示す。すなわち、図６及び図７に示すように、ライン２０１，２０２の中心から５０ｋｍの地点に、距離の値［５０ｋｍ］を表示している。また、図６に示すように、ライン２０１の緯度［３３°］や、ライン２０２の経度［１３１°］を表示している。さらに、図７に示すように、深さを示すライン２０３に、５０ｋｍごとの深さの値を表示させてある。図８以降ではこれらの数値表示は省略するが、同様に表示するようにしてもよい。図６や図７で距離や深さを５０ｋｍ単位で表示するようにしたのは一例であり、より細かい距離や深さを表示してもよい。
また、図６〜図１３の例では、各画像２００ａ〜２００ｈの右隅に、震源を表示する日時の範囲を設定するスライダ２１０と、そのスライダ２１０内の日時の始点のノブ２１１及び終点のノブ２１２を示している。

【0041】

図６〜図１３の３次元震源分布画像２００ａ〜２００ｈでは、各震源の位置を球体（円形）で示し、その球体の大きさをマグニチュードの値によって、複数段階に変化させている。
なお、図６〜図１３の３次元震源分布画像２００ａ〜２００ｈでは、地表面Ｐ１での地図（地形）の表示形態として、地形を薄く表示した形態としたが、より明確に地形や地図が分かる表示形態としてもよい。この地図（地形図）の表示を行う際には、既に知られた活断層マップなどから取得した断層の位置を表示するようにしてもよい。

【0042】

この図６〜図１３に示した３次元震源分布画像２００ａ〜２００ｈは、九州の熊本地方で、２０１６年４月１６日１時２５分に本震が発生した地震（いわゆる平成２８年熊本地震）についての、本震と余震の３次元震源分布画像である。
平成２８年熊本地震の経緯について簡単に説明すると、２０１６年４月１４日２１時２６分に、深さ１１ｋｍ，マグニチュード６．５の前震と解釈できる地震があり、その後の２０１６年４月１６日１時２５分に、深さ１２ｋｍ，マグニチュード７．３の本震があった。

【0043】

図６〜図１３から分かるように、最初の前震と解釈できる地震の発生から長期間に亘って余震が発生している。余震の中では、２０１６年４月１５日０時３分にマグニチュード６．４の余震があり、他の余震はマグニチュード６．０未満であった。なお、ここでのマグニチュードの値は気象庁が発表したものである。
なお、図６〜図１０では、日付による色分けモードを設定した状態を示し、地震が発生した日付ごとに異なる色（図面上では異なる濃さ）で、各球体の印を表示させている。

【0044】

図６の３次元震源分布画像２００ａは、３次元空間を真上から見たアングルの例を示し、表示日時の範囲の始点及び終点のノブ２１１，２１２は、４月１４日及び４月２１日の全範囲とした例である。真上から見たアングルであるため、図６の例は３次元表示ではなく平面的な表示形態になっている。すなわち、図６に示す震源分布画像２００ａの表示形態は、従来から知られた、２次元の地図中に震源位置をプロットした震源分布図と同じである。平成２８年熊本地震の場合、図６に示すように、多くの余震がほぼ直線状に並び、この直線状に余震が並んだ箇所に、断層があることが推測される。

【0045】

図７に示す３次元震源分布画像２００ｂは、３次元空間をほぼ横から見たものである。この図７の３次元震源分布画像２００ｂは、表示範囲として、４月１４日の２１時から４月１５日の２１時まで、つまり前震が発生してから本震が発生する直前での震源分布を示す。この期間には、マグニチュード６を超える地震が２回あり、そのマグニチュード６を超える地震の震源位置の円形の印ｅ１，ｅ２が、他の印とは異なる表示形態（図では網掛け状の表示形態）で示される。
また、図７に示す３次元震源分布画像２００ｂは、利用者によるクリックなどの操作で、１つの印ｅ１が選択されたときの震源の詳細の表示例を示す。すなわち、印ｅ１が選択された場合、３次元震源分布画像２００ｂは、印ｅ１に該当する震源の詳細（マグニチュードの値、発生日付及び時刻、緯度及び経度、深さ等）を、文字及び数字で印ｅ１の脇に表示している。この震源の詳細表示は、画像中のいずれの震源についても可能である。また、図７以外の他の表示（図６、図８−図１３に示す表示）を行った場合にも、震源の選択操作で、該当する詳細表示が行われる。

【0046】

図８に示す３次元震源分布画像２００ｃは、図７の３次元震源分布画像２００ｂから表示範囲の終点（ノブ２１２の位置）をシフトさせて、４月１４日の２１時から４月１７日の２時まで、つまり前震が発生してから本震が発生した直後までの震源分布を示している。
図８で新たに加わった印ｅ３が、本震の位置を示すものである。図８の画像２００ｃと図７の画像２００ｂを比較すると分かるように、本震の直後には多数の余震が発生していることが分かる。

【0047】

図９に示す３次元震源分布画像２００ｄは、図８の３次元震源分布画像２００ｃからさらに表示範囲の終点（ノブ２１２の位置）をシフトさせて、４月１４日の２１時から４月１８日の２時まで、つまり前震が発生してから本震の約１日後までの震源分布を示したものである。
このように本震発生から約１日の経過で、さらに余震が増えていく状態が、図９に示す３次元震源分布画像２００ｄから分かる。

【0048】

図１０に示す３次元震源分布画像２００ｅは、図９の３次元震源分布画像２００ｄからさらに表示範囲の終点（ノブ２１２の位置）をシフトさせて、４月１４日の２１時から４月２１日の０時までの約１週間の震源分布を示す。
この図１０から分かるように、マグニチュード６を超える３つの地震の印ｅ１，ｅ２，ｅ３の周囲に、多数の余震が発生した状況が分かる。ここで、発生した余震の日付ごとに印の色（濃さ）が異なると共に、その印の球体の大きさが地震の規模（マグニチュード）に対応しているため、１日ごとにどの程度の余震がどの位置や深さで発生したかが評価できるようになる。

【0049】

図１０に示す３次元震源分布画像２００ｅから各地震の３次元分布を見ると、地震の発生位置として、グループＧａと、グループＧｂの２つに分かれていることが分かる。図１０は１つのアングルから見た状態であるが、図５に示す水平方向θ１の回転で他の方位から見ることで、２つのグループＧａ，Ｇｂに分かれていることが、より明確に分かる。
地震の発生位置が、２つのグループＧａ，Ｇｂに分かれているということは、この２つのグループＧａ，Ｇｂが、別の断層で起きた破壊による誘発地震である可能性が高いと評価することができる。

【0050】

つまり、図１０は、２つの断層に分かれて地震が発生した状態を示しており、本震の印ｅ３を含む断層の破壊に伴って、別の断層の破壊も生じたことが、図１０に示す３次元震源分布画像２００ｅから明確に理解できる。従来の平面的な震源分布図を使った解析では、１つの断層が連続したものであるのか、あるいは別の２つの断層であるのか、簡単には評価できない場合が多々あった。これに対して、本例のように立体的に示した場合には、画像から簡単に地震の発生状況等を適切に評価することができる。特に、本例の場合には、図７〜図１０に示すように、表示する震源の時間範囲を様々に変更して表示させることで、震源位置の時間的な変化を三次元的に評価することができ、断層位置などを適切に推定できるようになる。
また、各地震の発生状況から、断層破壊がどの程度の深さでどの程度生じたかについても評価できるようになる。

【0051】

図１１に示す３次元震源分布画像２００ｆは、震源の印の表示形態を、深さによる色分けモードとした例である。この例は、深さが深い程、濃い色で表示させた例である。但し、表示色閾値として設定されたマグニチュードを超えた地震の印ｅ１，ｅ２，ｅ３については、網掛け表示として、この３つの地震の印ｅ１，ｅ２，ｅ３が、区別できようにしている。
この図１１に示す３次元震源分布画像２００ｆから、各地震の深さを評価することができる。

【0052】

図１２及び図１３に示す３次元震源分布画像２００ｆ及び２００ｇは、近似平面Ｐ３を表示させた例である。図１２の震源分布画像２００ｆは真上から見た平面図であり、図１３の３次元震源分布画像２００ｇは、立体的に近似平面Ｐ３を見た図である。
この近似平面Ｐ３は、各震源の位置を近似して求めた平面であり、図１２のように真上から見た平面図の震源分布画像２００ｆでは、近似平面Ｐ３は直線として表示される。一方、立体的に近似平面Ｐ３を見た場合には、近似平面Ｐ３の地下空間内での配置状態が分かるようになる。この近似平面Ｐ３から、地震により生じた断層の評価ができるようになる。
なお、図１２及び図１３の例は、表示範囲全体から近似平面Ｐ３を求めた例であるが、例えば図１０で説明したそれぞれのグループＧａ、Ｇｂごとの地震の表示範囲を設定して、表示範囲ごとの近似平面Ｐ３を表示させる処理を行うこともできる。このように表示範囲ごとの近似平面Ｐ３を表示することにより、より実際の断層の位置に近い近似平面Ｐ３を表示させることができる。
したがって、本例によると、地震により生じた地下の断層を適切に評価できるようになる。なお、画面上で、それぞれ異なる条件で得た複数の近似平面を同時に表示するようにして、それぞれの近似平面の比較や評価ができるようにしてもよい。

【0053】

［７．変形例］
なお、図６から図１３のそれぞれの図に示した３次元震源分布画像の表示例はあくまでも一例であり、その他の表示形態としてもよい。すなわち、図６〜図１３の例では、各震源の位置に、マグニチュードの値に対応した球体を印として配置した画像としたが、その他の印を配置した画像としてもよい。マグニチュードの値に対応して球体の大きさを変化させる点についても一例であり、その他の表示形態で、マグニチュードなどの地震の規模が分かるようにしてもよい。

【0054】

また、ここまで説明した例では、地震情報として、震源のマグニチュードの値と座標位置と深さの値を得るようにしたが、その他の情報を取得して、３次元震源分布画像を生成する際に、その取得した情報を使った表示を行うようにしてもよい。例えば、地震波が地中を伝わる速度の情報を取得して、３次元分布画像として、震源情報と重ね合わせ表示を行うようにしてもよい。

【0055】

また、近似平面を算出する際に、各震源からの距離が最小となる平面を算出する方法は、好適な一例を示すものであり、その他の演算処理を用いて複数の震源の集合状態を表した面を算出して近似平面としてもよい。また、演算処理で求まる近似面の少なくとも一部が曲面である場合には、その曲面を近似面として表示してもよい。また、３次元空間内に配置された震源のいくつかを利用者が選択し、その選択された震源を元にひとつあるいは複数の近似平面を表示するモードを用意してもよい。

【0056】

また、図４に示す設定画面は、ひとつの画面ですべての地震情報取得／表示条件を入力する例としたが、いくつかの条件ごとに設定画面をわけるようにしてもよい。

【0057】

また、図１に示す震源情報可視化システム１００は、サーバー装置１１０とクライアント装置１２０で構成した例とした。これに対して、１つのコンピューター装置内で地震情報の収集から３次元震源分布画像まで全ての処理を行うようにしてもよい。図１に示すようにクライアント装置１２０がサーバー装置１１０にアクセスする場合にも、３次元震源分布画像の生成までをサーバー装置１１０が実行して、クライアント装置１２０側では、サーバー装置１１０から伝送される画像の表示やプリントなどを行うようにしてもよい。

【0058】

また、上述した実施の形態例では、作成した３次元震源分布画像を利用者に提示する提示部として、２次元の画像表示を行う表示部とした例を示したが、その他の形態で利用者に提示するようにしてもよい。例えば、立体画像表示が可能な表示装置を使用して、３次元震源分布画像そのものが立体的に表示されるようにしてもよい。また、表示部の代わりに、プリンターによる印刷などで３次元震源分布画像を利用者に提示してもよい。

【0059】

さらに、本発明の請求項に記載した構成や処理は、上述した実施の形態の例に限定されるものではない。本発明の要旨を逸脱しない限り、種々の改変、組み合わせ、他の実施の形態例が生じうることは、当業者にとって当然のことと理解される。

【符号の説明】

【0060】

１０…地震情報データベース、１００…地震情報可視化システム、１１０…サーバー装置、１１１…地震情報収集部、１１２…地震情報蓄積部、１１３…地震情報検索部、１２０…クライアント装置、１２１…地震情報取得部、１２２…読出し範囲設定部、１２３…３次元画像生成部、１２４…表示条件設定部、１２５…表示部、１２６…出力部、Ｃ…コンピューター装置、Ｃ１…ＣＰＵ、Ｃ２…ＲＯＭ、Ｃ３…ＲＡＭ、Ｃ４…バスライン、Ｃ５…ストレージ装置、Ｃ６…入出力部、Ｃ７…ディスプレイ

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】