(57)【要約】  設置されているセンサから十分な量のデータを取得できない場合であっても、高速に、信頼性の高い物理量の推定を行うために、本発明は、シミュレーションの条件として入力された値と、前記入力された値に基づく所定のシミュレーションの結果として算出された、一つ以上の位置における物理量と、を含むシナリオを複数保持するシナリオ保持部と、計測された物理量及び前記物理量が計測された位置に基づく所定の計算によって前記複数のシナリオに含まれる物理量を合成することによって、前記一つ以上の位置における物理量を推定するシナリオ合成部と、前記推定された物理量を保持する合成シナリオ保持部と、を有する。
2014020747.tif 000002
【特許請求の範囲】
【請求項１】
シミュレーションの条件として入力された値と、前記入力された値に基づく所定のシミュレーションの結果として算出された、一つ以上の位置における物理量と、を含むシナリオを複数保持するシナリオ保持部と、
計測された物理量及び前記物理量が計測された位置に基づく所定の計算によって前記複数のシナリオに含まれる物理量を合成することによって、前記一つ以上の位置における物理量を推定するシナリオ合成部と、
前記推定された物理量を保持する合成シナリオ保持部と、を有することを特徴とする物理量推定システム。
【請求項２】
請求項１に記載の物理量推定システムであって、
前記シナリオ合成部は、前記計測された物理量及び前記物理量が計測された位置に基づく所定の計算として、前記計測された物理量と、前記複数のシナリオに含まれる、前記物理量が計測された位置における物理量との類似度に従って定めた重みを、前記複数のシナリオに含まれる物理量に付けた重み付け演算を実行することを特徴とする物理量推定システム。
【請求項３】
請求項２に記載の物理量推定システムであって、
前記複数のシナリオを、前記各シナリオに含まれる物理量に基づくクラスタリングによって、各々が複数のシナリオを含む複数のシナリオクラスタに分類し、前記各シナリオクラスタに含まれる複数のシナリオの前記各位置における物理量の統計量を算出し、前記各シナリオクラスタに含まれる複数のシナリオの前記各位置における物理量のばらつきが所定の値より大きい場合、前記各位置を重要スポットとして特定し、前記各シナリオが属する前記シナリオクラスタを特定する情報、及び、前記各シナリオクラスタの前記重要スポットを特定する情報を前記シナリオ保持部に格納するシナリオ管理部をさらに有し、
前記シナリオ合成部は、
前記各シナリオクラスタについて算出された前記統計量と、前記計測された物理量とを比較することによって、前記計測された物理量と最も類似するシナリオクラスタを特定し、
前記特定されたシナリオクラスタの前記重要スポットにおける物理量を、前記重み付け演算によって推定し、
前記特定されたシナリオクラスタの前記重要スポット以外の位置における物理量を、前記重要スポットにおいて使用されるものとは異なる重みを付けた重み付け演算によって推定することを特徴とする物理量推定システム。
【請求項４】
請求項３に記載の物理量推定システムであって、
第１の種類の前記シミュレーションの結果に含まれる前記一つ以上の位置における第１の種類の前記物理量に基づいて、第２の種類の前記シミュレーションの条件として入力される値を生成するシミュレーション条件生成部をさらに有し、
前記各シナリオは、少なくとも、前記第１の種類のシミュレーションの条件、前記第１の種類のシミュレーションの結果として得られる前記第１の種類の物理量、及び、前記第２の種類のシミュレーションの結果として得られる第２の種類の物理量を含み、
前記シナリオ合成部は、前記各シナリオクラスタについて算出された前記統計量と、第１の種類のセンサによって計測された前記第１の種類の物理量と、第２の種類のセンサによって計測された前記第２の種類の物理量と、を比較することによって、前記計測された物理量と最も類似するシナリオクラスタを特定することを特徴とする物理量推定システム。
【請求項５】
請求項４に記載の物理量推定システムであって、
入力された前記第１の種類の物理量に基づいて、一つ以上の位置における前記第１の種類の物理量を算出する第１の種類のシミュレータ部と、入力された前記第２の種類の物理量に基づいて、一つ以上の位置における前記第２の種類の物理量を算出する第２の種類のシミュレータ部と、をさらに有することを特徴とする物理量推定システム。
【請求項６】
請求項３に記載の物理量推定システムであって、
一つの前記シミュレーション条件を入力されると、各々が、前記入力されたシミュレーション条件に含まれる値を含む範囲内の複数の値の各々を含む、複数のシミュレーション条件を生成するシミュレーション条件生成部をさらに有することを特徴とする物理量推定システム。
【請求項７】
請求項３に記載の物理量推定システムであって、
前記各シナリオクラスタについて算出された前記統計量は、前記各シナリオクラスタ内の前記複数のクラスタの各位置における前記物理量に同一の重みを付けて計算された平均値を含み、
前記シナリオ合成部は、
前記特定されたシナリオクラスタに含まれる前記複数のシナリオの前記重要スポットにおける物理量に、前記計測された物理量との類似度が高いほど大きい重みを付けて計算した平均値を、前記特定されたシナリオクラスタの前記重要スポットにおける物理量として推定し、
前記特定されたシナリオクラスタの前記重要スポット以外の位置について算出された前記統計量に含まれる平均値を、前記特定されたシナリオクラスタの前記重要スポット以外の位置における物理量として推定することを特徴とする物理量推定システム。
【請求項８】
請求項３に記載の物理量推定システムであって、
前記シナリオ合成部は、
前記特定されたシナリオクラスタに含まれる前記複数のシナリオの前記重要スポットにおける物理量に、前記計測された物理量との類似度が高いほど大きい重みを付けて計算した平均値を、前記特定されたシナリオクラスタの前記重要スポットにおける物理量として推定し、
前記特定されたシナリオクラスタに含まれる前記複数のシナリオの前記重要スポット以外の位置における物理量のうち、前記計測された物理量との類似度が最も高いものの重みとして０より大きい値を付け、それ以外の物理量の重みとして０を付けて計算した平均値を、前記特定されたシナリオクラスタの前記重要スポット以外の位置における物理量として推定することを特徴とする物理量推定システム。
【請求項９】
物理量推定システムが実行する物理量推定方法であって、
前記物理量推定システムは、プロセッサと、前記プロセッサに接続されるインターフェースと、前記プロセッサに接続される記憶装置と、を有し、
前記記憶装置は、シミュレーションの条件として入力された値と、前記入力された値に基づく所定のシミュレーションの結果として算出された、一つ以上の位置における物理量と、を含むシナリオを複数保持し、
前記物理量推定方法は、
前記インターフェースが、計測された物理量の入力を受ける第１手順と、
前記プロセッサが、前記計測された物理量及び前記物理量が計測された位置に基づく所定の計算によって前記複数のシナリオに含まれる物理量を合成することによって、前記一つ以上の位置における物理量を推定し、前記推定された物理量を前記記憶装置に格納する第２手順と、を含むことを特徴とする物理量推定方法。
【請求項１０】
請求項９に記載の物理量推定方法であって、
前記第２手順は、前記プロセッサが、前記計測された物理量と、前記複数のシナリオに含まれる、前記物理量が計測された位置における物理量との類似度に従って定めた重みを、前記複数のシナリオに含まれる物理量に付けた重み付け演算を実行する手順を含むことを特徴とする物理量推定方法。
【請求項１１】
請求項１０に記載の物理量推定方法であって、
前記プロセッサが、前記複数のシナリオを、前記各シナリオに含まれる物理量に基づくクラスタリングによって、各々が複数のシナリオを含む複数のシナリオクラスタに分類し、前記各シナリオクラスタに含まれる複数のシナリオの前記各位置における物理量の統計量を算出し、前記各シナリオクラスタに含まれる複数のシナリオの前記各位置における物理量のばらつきが所定の値より大きい場合、前記各位置を重要スポットとして特定し、前記各シナリオが属する前記シナリオクラスタを特定する情報、及び、前記各シナリオクラスタの前記重要スポットを特定する情報を前記記憶装置に格納する第３手順をさらに含み、
前記第２手順は、
前記プロセッサが、前記各シナリオクラスタについて算出された前記統計量と、前記計測された物理量とを比較することによって、前記計測された物理量と最も類似するシナリオクラスタを特定する第４手順と、
前記プロセッサが、前記特定されたシナリオクラスタの前記重要スポットにおける物理量を、前記重み付け演算によって推定する第５手順と、
前記プロセッサが、前記特定されたシナリオクラスタの前記重要スポット以外の位置における物理量を、前記重要スポットにおいて使用されるものとは異なる重みを付けた重み付け演算によって推定する第６手順と、を含むことを特徴とする物理量推定方法。
【請求項１２】
請求項１１に記載の物理量推定方法であって、
前記プロセッサが、第１の種類の前記シミュレーションの結果に含まれる前記一つ以上の位置における第１の種類の前記物理量に基づいて、第２の種類の前記シミュレーションの条件として入力される値を生成する第７手順をさらに含み、
前記各シナリオは、少なくとも、前記第１の種類のシミュレーションの条件、前記第１の種類のシミュレーションの結果として得られる前記第１の種類の物理量、及び、前記第２の種類のシミュレーションの結果として得られる第２の種類の物理量を含み、
前記第４手順は、前記各シナリオクラスタについて算出された前記統計量と、第１の種類のセンサによって計測された前記第１の種類の物理量と、第２の種類のセンサによって計測された前記第２の種類の物理量と、を比較する手順を含むことを特徴とする物理量推定方法。
【請求項１３】
請求項１２に記載の物理量推定方法であって、
前記プロセッサが、入力された前記第１の種類の物理量に基づいて、一つ以上の位置における前記第１の種類の物理量を算出し、入力された前記第２の種類の物理量に基づいて、一つ以上の位置における前記第２の種類の物理量を算出する第８手順をさらに含むことを特徴とする物理量推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明はシミュレーションを利用して物理量を推定する技術に関し、特に、複数種類のセンサによって取得された、位置情報及び時刻情報を伴う種々のデータを利用したシミュレーションによって任意の物理量を推定する技術に関する。
【背景技術】
【０００２】
センサが実際に計測したデータを初期条件として与えられると、所定の計算を行うことによって、現在の（又は任意の時刻の）物理量を推定して出力する、リアルタイムシミュレーション技術が知られている。【０００３】
例えば、特開２０１０−５４４６０号公報（特許文献１）には、「第一推定値の気象予報値を観測値で修正した大気解析データに含まれる風向及び風速データを基にデータ補間してＮ時の波浪予報データを推算すると、これをＮ+α時の波浪予測のための初期条件値として波浪推算プログラムに設定し、Ｎ+α時の波浪を予測する。本処理は大気解析データが更新される所定時間間隔（α）ごとに繰り返し実行される。また、地球規模での気象変化を予測計算した全球気象予報データや局地気象データの風データを波浪推算プログラムに設定して波浪予測を行う。これにより、時間経過に伴う波浪状態の変化予報データに誤差拡大がない高精度な波浪予報データをあらゆる地点について迅速に算出することができる」と記載されている。
【先行技術文献】
【特許文献】
【０００４】
特許文献１：特開２０１０−５４４６０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
防災又は商業利用等を目的として種々の物理量を計測し、計測した物理量を含むデータを提供するセンサが、人間の生活域（又はその周辺の地域）に設置される。その例をいくつか挙げるとすれば、地震計、雨量計、風速計、風向計、河川等の水位計、各家庭又は事業所の電力計、及び、道路に設置された交通量センサ等である。しかし、センサを設置・維持するコスト、災害によるセンサ自体又はセンサからのデータが経由する通信経路の障害、及び、当該通信経路における通信速度、等の種々の要因のため、目的に応じた十分な量のデータを得られないことが想定される。そこで、センサが実際に取得できない地点又は時刻における物理量を、シミュレーションによって推定することが考えられる。しかし、通常、そのようなシミュレーションの実行には多くの時間を要するため、実際にセンサから取得したリアルタイムデータを初期条件として使用したシミュレーションによって、現在の（又は未来の所望の時点における）所望の地点の物理量を推定することは難しい。
【課題を解決するための手段】
【０００６】
上記の課題を解決するために、本発明は、シミュレーションの条件として入力された値と、前記入力された値に基づく所定のシミュレーションの結果として算出された、一つ以上の位置における物理量と、を含むシナリオを複数保持するシナリオ保持部と、計測された物理量及び前記物理量が計測された位置に基づく所定の計算によって前記複数のシナリオに含まれる物理量を合成することによって、前記一つ以上の位置における物理量を推定するシナリオ合成部と、前記推定された物理量を保持する合成シナリオ保持部と、を有する。
【発明の効果】
【０００７】
本発明の一実施形態によれば、設置されているセンサから十分な量のデータを取得できない場合であっても、事前に実行されたシミュレーション結果を利用することによって、高速に、信頼性の高い物理量の推定を行うことができる。
【図面の簡単な説明】
【０００８】
【図１】本発明の第１の実施形態の計算機システムの構成を示す機能ブロック図である。
【図２】本発明の第１の実施形態の物理量推定システムのハードウェア構成を示すブロック図である。
【図３Ａ】本発明の第１の実施形態の物理量推定システムが実行する事前処理のフローチャートである。
【図３Ｂ】本発明の第１の実施形態の物理量推定システムが実行するリアルタイム処理のフローチャートである。
【図４】本発明の第１の実施形態の物理量推定システムが実行するシナリオＤＢ構築処理のフローチャートである。
【図５】本発明の第１の実施形態の物理量推定システムが実行するシミュレーション条件生成処理のフローチャートである。
【図６】本発明の第１の実施形態の物理量推定システムが実行するシミュレーション処理のフローチャートである。
【図７】本発明の第１の実施形態の物理量推定システムが実行するシミュレーション処理の一例である氾濫シミュレーションの説明図である。
【図８】本発明の第１の実施形態の物理量推定システムが実行するシナリオＤＢ構築処理のフローチャートである。
【図９Ａ】本発明の第１の実施形態の物理量推定システムによって作成されるシナリオの説明図である。
【図９Ｂ】本発明の第１の実施形態のシナリオＤＢに含まれるシナリオ管理テーブルの説明図である。
【図９Ｃ】本発明の第１の実施形態のシナリオＤＢに含まれるシナリオ物理量管理テーブルの説明図である。
【図１０】本発明の第１の実施形態の物理量推定システムが実行するシナリオクラスタ生成処理のフローチャートである。
【図１１Ａ】本発明の第１の実施形態のシナリオＤＢに含まれるシナリオクラスタリングテーブルの説明図である。
【図１１Ｂ】本発明の第１の実施形態のシナリオＤＢに含まれるシナリオクラスタ管理テーブルの説明図である。
【図１１Ｃ】本発明の第１の実施形態のシナリオＤＢに含まれるシナリオクラスタ管理テーブルの具体例の説明図である。
【図１１Ｄ】本発明の第１の実施形態のシナリオＤＢに含まれる重要スポット管理テーブルの説明図である。
【図１２】本発明の第１の実施形態の物理量推定システムが実行するシナリオ合成処理のフローチャートである。
【図１３】本発明の第１の実施形態の物理量推定システムが実行する実物理量入力処理のフローチャートである。
【図１４Ａ】本発明の第１の実施形態の物理量推定システムが現況入力システムから受け付けるデータの説明図である。
【図１４Ｂ】本発明の第１の実施形態の物理量推定システムがセンサから受け付けるデータの説明図である。
【図１５】本発明の第１の実施形態の物理量推定システムが実行するシナリオ合成処理のフローチャートである。
【図１６】本発明の第１の実施形態の物理量推定システムが実行するシナリオ合成処理の概念の説明図である。
【図１７】本発明の第１の実施形態のアプリケーションシステムが実行する誘導情報生成処理のフローチャートである。
【図１８】本発明の第２の実施形態の計算機システムの構成を示す機能ブロック図である。
【図１９Ａ】本発明の第２の実施形態の物理量推定システムが実行する事前処理のフローチャートである。
【図１９Ｂ】本発明の第２の実施形態の物理量推定システムが実行するリアルタイム処理のフローチャートである。
【図２０】本発明の第２の実施形態の物理量推定システムが実行するシナリオＤＢ構築処理のフローチャートである。
【図２１】本発明の第２の実施形態のシミュレータが実行するシミュレーション処理の一例である津波シミュレーションの説明図である。
【図２２Ａ】本発明の第２の実施形態の物理量推定システムによって作成されるシナリオの説明図である。
【図２２Ｂ】本発明の第２の実施形態のシナリオＤＢに含まれるシナリオ管理テーブルの説明図である。
【図２２Ｃ】本発明の第２の実施形態のシナリオＤＢに含まれる水深管理テーブルの説明図である。
【図２２Ｄ】本発明の第２の実施形態のシナリオＤＢに含まれるスマートメータ故障管理テーブルの説明図である。
【図２３Ａ】本発明の第２の実施形態の物理量推定システムが現況入力システムから受け付けるデータの説明図である。
【図２３Ｂ】本発明の第２の実施形態の物理量推定システムがセンサから受け付けるデータの説明図である。
【図２３Ｃ】本発明の第２の実施形態の物理量推定システムが別の種類のセンサから受け付けるデータの説明図である。
【図２４】本発明の第２の実施形態のアプリケーションシステムが実行する被害地域情報生成処理のフローチャートである。
【発明を実施するための形態】
【０００９】
以下、添付図面に基づいて、本発明の実施形態を説明する。
【００１０】
＜第１の実施形態＞
図１は、本発明の第１の実施形態の計算機システムの構成を示す機能ブロック図である。
【００１１】
本実施形態の計算機システムは、物理量推定システム１００、シナリオ例入力システム１１０、センサ１２０、現況入力システム１３０及びアプリケーションシステム１４０によって構成される。シナリオ例入力システム１１０はネットワーク１５０Ａを介して、センサ１２０及び現況入力システム１３０はネットワーク１５０Ｂを介して、アプリケーションシステム１４０はネットワーク１５０Ｃを介して、それぞれ物理量推定システム１００に接続される。
【００１２】
物理量推定システム１００は、センサ１２０及び現況入力システム１３０から取得したデータに基づいて物理量を推定する計算機であり、シミュレーション条件生成部１０１、シミュレータ１０２、シナリオ管理部１０３、シナリオデータベース（ＤＢ）１０４、入力受付部１０５、シナリオ合成部１０６、合成シナリオＤＢ１０７及びシナリオ配信部１０８を有する。これらの各部の詳細及び物理量推定システム１００のハードウェア構成については後述する。
【００１３】
なお、本実施形態において、「シミュレータ」とは、物理的な条件が入力されると、演繹的に他の物理量を計算する機能を有するモジュールである。すなわち、一般的にシミュレータと呼ばれるものは本実施形態のシミュレータに該当するが、一般的にシミュレータと呼ばれないものであっても、上記の機能を有するものであれば本実施形態のシミュレータに該当する。本実施形態の物理量推定システム１００は、種類の異なる複数のシミュレータ１０２（例えば、地震シミュレータ及び氾濫シミュレータ等）を有してもよい。
【００１４】
また、本実施形態において、「シナリオ」とは、実際に起こりうる物理的な現象の記録であり、その物理的な現象が発生する確率の情報を含んでもよい。具体的には、シナリオは、シミュレーションに入力される条件と、そのシミュレーションの結果（すなわち物理量）との集合である。例えば、震源地と地震の規模とがシミュレーションの条件として入力されると各地点の揺れの大きさを計算する地震シミュレータの場合、入力された震源地及び地震の規模と、計算された各時刻における各地点の揺れの大きさの値と、の集合が一つのシナリオである。あるいは、堤防の決壊地点とそこから流出する水量とがシミュレーションの条件として入力されると各地点の浸水深を計算する氾濫シミュレータの場合、入力された決壊地点及び水量と、計算された各時刻における各地点の浸水深の値と、の集合が一つのシナリオである。
【００１５】
シナリオ例入力システム１１０は、物理量推定システム１００の開発者、管理者又はユーザ等（以下、単にユーザ等とも記載する）がシナリオ例を入力するために使用するものであり、ネットワーク１５０Ａに接続された計算機によって実現される。例えば、シナリオ例入力システム１１０は、当該計算機上で動作するＷｅｂブラウザであってもよい。なお、ここで入力されるシナリオ例とは、例えば、複数のシミュレーションの間の関係を定義する情報であり、その詳細については後述する。
【００１６】
センサ１２０は、所定のタイミングで（例えば定期的に、又は物理量推定システム１００等からの要求に応じて）物理量を計測し、計測した物理量を含むデータ（すなわち計測した物理量そのもの又はそれに基づいて計算された値）を、ネットワーク１５０Ｂを介して物理量推定システム１００に送信する。本実施形態の図１には一つのセンサ１２０のみを図示しているが、実際には複数の種類の複数のセンサ１２０が各地点に設置される。各センサ１２０は、例えば、地震計、雨量計、河川に設置された水位計、各家庭又は事業所に設置された電力計（例えばいわゆるスマートメータ）、又は、道路に設置された交通量センサ等である。各センサ１２０から送信されるデータは、計測した物理量（例えば水位又は電力量等）に加えて、少なくともその物理量が計測された時刻の情報を含む。物理量推定システム１００は、各センサが設置された位置の情報（座標値等）を保持している。
【００１７】
なお、センサ１２０は、所定の位置に固定的に設置されたものではなく、移動可能なものであってもよい。その場合には、センサ１２０から送信されるデータに、センサ１２０の位置情報も含まれる。移動可能なセンサの一例は、走行する自動車に設置され、所定のタイミングでその自動車の位置情報を送信する車載端末（いわゆるプローブカーに設置されるプローブ端末）であり、このような車載端末からの情報に基づいて道路の交通量を計測することができる。このような場合、物理量推定システム１００は、センサ１２０から送信された位置情報を保持し、それに基づいてセンサ１２０の位置を特定する。
【００１８】
現況入力システム１３０は、物理量推定システム１００のユーザ等が、物理量推定システム１００による物理量の推定に利用される情報を入力するために使用するものであり、ネットワーク１５０Ｂに接続された計算機によって実現される。例えば、現況入力システム１３０は、当該計算機上で動作するＷｅｂブラウザであってもよい。例えば、地震が発生したときに、ユーザ等は、現況入力システム１３０を使用して、自身が体感した揺れの程度を示す情報を入力してもよい。現況入力システム１３０は、入力された情報を物理量推定システム１００に送信する。
【００１９】
この場合、物理量推定システム１００は、現況入力システム１３０から受信した情報を、複数のセンサ１２０のうち地震計から送信されたデータと同様に扱ってもよい。ただし、ユーザ等が体感震度を入力した場合、その精度は地震計の制度より低いと考えられる。ユーザ等が他の物理量を入力する場合の精度も、一般に、その物理量を計測するセンサの精度より低い。すなわち、物理量推定システム１００は、現況入力システム１３０からの情報を、精度の低いセンサ１２０からの情報として扱うことができる。
【００２０】
現況入力システム１３０を用いることによって、センサ１２０では感知困難な物理量、又は、センサ１２０が設置されていない地点における物理量を人手によって入力することができる。
【００２１】
アプリケーションシステム１４０は、物理量推定システム１００によって推定された物理量に基づく任意の処理を実行する。例えば、物理量推定システム１００が任意の時刻における震度を推定した場合、アプリケーションシステム１４０は、その震度に基づく地震警報、津波警報又は避難情報等を送信してもよい。物理量推定システム１００が任意の時刻における洪水による浸水深を推定した場合、アプリケーションシステム１４０は、浸水していない道路のみを経由する経路を探索してもよい。アプリケーションシステム１４０は、例えば、ネットワーク１５０Ｃに接続された計算機を含むシステムによって実現される。
【００２２】
ネットワーク１５０Ａ、１５０Ｂ及び１５０Ｃは、それらに接続されたシステム間の通信を可能にするものであれば、どのようなものであってもよい。例えば、ネットワーク１５０Ａ、１５０Ｂ及び１５０Ｃは、いわゆるインターネット又は電話網のような公衆ネットワークであってもよいし、物理量推定システム１００のための専用ネットワークであってもよいし、それらの組み合わせであってもよい。また、図１には三つのネットワーク１５０Ａ、１５０Ｂ及び１５０Ｃを示したが、これらが一つのネットワークによって実現されてもよい。以下の説明においてネットワーク１５０が参照された場合、それは、ネットワーク１５０Ａ、１５０Ｂ及び１５０Ｃのいずれか、又は、ネットワーク１５０Ａ、１５０Ｂ及び１５０Ｃを実現する一つのネットワークである。
【００２３】
図２は、本発明の第１の実施形態の物理量推定システム１００のハードウェア構成を示すブロック図である。
【００２４】
本実施形態の物理量推定システム１００は、相互に接続されたプロセッサ２０１、メインメモリ２０２、入力装置２０３、出力装置２０４、インターフェース（Ｉ／Ｆ）２０５及び記憶装置２０６を有する計算機である。
【００２５】
プロセッサ２０１は、メインメモリ２０２に格納されたプログラムを実行する。
【００２６】
メインメモリ２０２は、例えば半導体メモリであり、プロセッサ２０１によって実行されるプログラム及びプロセッサ２０１によって参照されるデータを格納する。具体的には、記憶装置２０６に格納されたプログラム及びデータの少なくとも一部が、必要に応じてメインメモリ２０２にコピーされる。
【００２７】
入力装置２０３は、物理量推定システム１００のユーザ等からの入力を受ける。入力装置２０３は、例えばキーボード及びマウス等を含んでもよい。
【００２８】
出力装置２０４は、例えば液晶表示装置のような画像表示装置である。
【００２９】
Ｉ／Ｆ２０５は、ネットワーク１５０に接続され、シナリオ例入力システム１１０、センサ１２０、現況入力システム１３０及びアプリケーションシステム１４０と通信するインターフェースである。ネットワーク１５０Ａ、１５０Ｂ及び１５０Ｃが独立している場合、物理量推定システム１００は、それぞれのネットワークに接続される複数のＩ／Ｆ２０５を有する。
【００３０】
記憶装置２０６は、例えばハードディスク装置（ＨＤＤ）又はフラッシュメモリのような不揮発性の記憶装置である。本実施形態の記憶装置２０６には、少なくとも、シミュレーション条件生成部１０１、シミュレータ１０２、シナリオ管理部１０３、入力受付部１０５、シナリオ合成部１０６、シナリオ配信部１０８、シナリオＤＢ１０４及び合成シナリオＤＢ１０７が格納される。
【００３１】
シナリオＤＢ１０４及び合成シナリオＤＢ１０７は、後述するように、それぞれシナリオ管理部１０３によって生成されたシナリオ及びシナリオ合成部１０６によって合成されたシナリオを保持する。
【００３２】
シミュレーション条件生成部１０１、シミュレータ１０２、シナリオ管理部１０３、入力受付部１０５、シナリオ合成部１０６及びシナリオ配信部１０８は、プロセッサ２０１によって実行されるプログラムである。以下の説明においてこれらの各部が実行する処理は、実際にはプロセッサ２０１によって実行される。なお、図２には一つのシミュレータ１０２のみを示しているが、本実施形態の物理量推定システム１００は互いに種類の異なる複数のシミュレータ１０２（例えば、地震シミュレータプログラム及び氾濫シミュレータプログラム等）を保持する。
【００３３】
シナリオ例入力システム１１０、現況入力システム１３０及びアプリケーションシステム１４０は、物理量推定システム１００と同様のハードウェアを有する計算機によって実現できるため、それらのハードウェア構成については図示を省略する。
【００３４】
シナリオ例入力システム１１０の記憶装置（図示省略）又はメインメモリ（図示省略）にはシナリオ例入力システムを実現するためのプログラム（例えばＷｅｂブラウザ、図示省略）が格納され、それをシナリオ例入力システム１１０のプロセッサ（図示省略）が実行する。
【００３５】
同様に、現況入力システム１３０の記憶装置（図示省略）等に格納された現況入力システムを実現するためのプログラム（例えばＷｅｂブラウザ、図示省略）を現況入力システム１３０のプロセッサ（図示省略）が実行する。アプリケーションシステム１４０の記憶装置（図示省略）等に格納されたアプリケーションプログラム（図示省略）をアプリケーションシステム１４０のプロセッサ（図示省略）が実行する。
【００３６】
次に、物理量推定システム１００が実行する処理について説明する。
【００３７】
物理量推定システム１００は、事前処理（図３Ａ）及びリアルタイム処理（図３Ｂ）を実行する。リアルタイム処理は、物理量推定システム１００がセンサ１２０等から取得した物理量に基づいて、現在又は未来の任意の時点の物理量を推定し、その推定結果を含むシナリオ（合成シナリオ）をアプリケーションシステム１４０に配信する処理である。事前処理は、リアルタイム処理を実行するための準備に相当し、物理量推定システム１００を導入してから最初にリアルタイム処理を実行する前に少なくとも１回実行する必要がある。
【００３８】
図３Ａは、本発明の第１の実施形態の物理量推定システム１００が実行する事前処理のフローチャートである。
【００３９】
物理量推定システム１００は、シナリオＤＢ構築処理を実行する（ステップ３０１）。これによって事前に想定されるシナリオが多数作成され、シナリオＤＢ１０４に蓄積される。シナリオＤＢ構築処理の詳細については図４等を参照して後述する。
【００４０】
図３Ｂは、本発明の第１の実施形態の物理量推定システム１００が実行するリアルタイム処理のフローチャートである。
【００４１】
物理量推定システム１００は、シナリオＤＢ１０４を参照してシナリオ合成処理を実行する（ステップ３０２）。これによって、リアルタイムに計測されたデータと、事前に作成されたシナリオとに基づいて新たなシナリオが合成され、その結果が合成シナリオＤＢ１０７に格納される。シナリオ合成処理の詳細については図１５等を参照して後述する。
【００４２】
アプリケーションシステム１４０は、合成シナリオＤＢ１０７を参照して種々のアプリケーションを実行することができる。例えば、アプリケーションシステム１４０は、合成シナリオＤＢ１０７を参照して誘導情報生成処理を実行し（ステップ３０３）、生成した誘導情報を、ネットワーク１５０Ｄを介して移動体端末３０４に送信してもよい。誘導情報生成処理の詳細については図１７を参照して後述する。
【００４３】
ここで、移動体端末３０４は、アプリケーションシステム１４０のユーザが所持する携帯電話機等の端末装置である。ネットワーク１５０Ｄは、例えば携帯電話網又は無線ＬＡＮ（Local Area Network）等であるが、アプリケーションシステム１４０と移動体端末３０４との間の通信を可能にするものであれば、どのようなものであってもよい。
【００４４】
ここで、事前処理（図３Ａ）の詳細について説明する。
【００４５】
図４は、本発明の第１の実施形態の物理量推定システム１００が実行するシナリオＤＢ構築処理のフローチャートである。この処理は、図３Ａのステップ３０１において実行される。
【００４６】
最初に、シミュレーション条件生成部１０１がシミュレーション条件生成処理を実行する（ステップ４０１）。シミュレーション条件生成処理の詳細については図５を参照して後述する。
【００４７】
次に、シミュレータ１０２が、ステップ４０１において生成されたシミュレーション条件に従ってシミュレーション処理を実行する（ステップ４０２）。シミュレーション処理の詳細については図６等を参照して後述する。
【００４８】
次に、シナリオ管理部１０３が、ステップ４０２において得られたシミュレーション結果を用いてシナリオＤＢ１０４を構築する（ステップ４０３）。さらに、シナリオ管理部１０３は、シナリオＤＢ１０４に格納されたシナリオを参照して、シナリオクラスタを生成し、その結果をシナリオＤＢ１０４に格納する（ステップ４０４）。これらの処理の詳細については図８及び図１０等を参照して後述する。
【００４９】
図５は、本発明の第１の実施形態の物理量推定システム１００が実行するシミュレーション条件生成処理のフローチャートである。この処理は、図４のステップ４０１においてシミュレーション条件生成部１０１によって実行される。
【００５０】
最初に、シミュレーション条件生成部１０１は、シミュレート条件受付処理を実行する（ステップ５０１）。具体的には、シミュレーション条件生成部１０１は、シナリオ例入力システム１１０を介して、物理量推定システム１００のユーザ等からシミュレート条件の入力を受け付ける。ユーザ等は、経験又は学術的な知見等に基づいて、実際に起こりうると考えられる事象をシミュレート条件として入力することができる。例えば、ユーザは、地震シミュレータに関するシミュレート条件として、発生する可能性がある地震の震源地及びその地震の規模をシミュレート条件として入力してもよい。あるいは、ユーザは、河川の氾濫時の浸水深を推定する氾濫シミュレータに関するシミュレート条件として、河川の堤防が決壊する可能性がある位置及びそこから流出する水の流量、又はある地点の降雨量をシミュレート条件として入力してもよい。
【００５１】
さらに、ユーザ等は、複数の種類のシミュレーションを関連付ける情報、より詳細には、あるシミュレータによるシミュレーションの結果に基づいて別のシミュレータのシミュレート条件を生成する方法を定義する情報を入力することができる。例えば、シミュレータ１０２の一つが上記の地震シミュレータであり、別の一つが、出火地点を入力されると延焼の範囲をシミュレートする火災シミュレータである場合、ユーザ等は、地震の揺れの大きさと出火の確率とを関連付ける情報を入力してもよい。シミュレータ１０２がさらに各家庭及び事業所等の使用電力量をシミュレートする電力シミュレータを含む場合、ユーザ等は、地震の揺れの大きさと停電の発生確率とを関連付ける情報、及び、火災の発生と停電の発生確率とを関連付ける情報を入力してもよい。
【００５２】
あるいは、シミュレータ１０２の一つが氾濫シミュレータ又は津波による浸水深を推定する津波シミュレータであり、別の一つがスマートメータによって計測される電力量を推定する電力シミュレータである場合、ユーザ等は、浸水深とスマートメータの故障の発生確率とを関連付ける情報を入力してもよい。このような例については本発明の第２の実施形態において説明する。
【００５３】
次に、シミュレーション条件生成部１０１は、境界条件抽出処理を実行する（ステップ５０２）。ここで、境界条件とは、シミュレータ１０２にシミュレート条件として入力される具体的な数値である。シミュレーション条件生成部１０１は、ステップ５０１において入力されたシミュレート条件から境界条件を抽出する。例えば、シミュレータ１０２に入力可能な具体的な数値がステップ５０１においてシミュレート条件として入力された場合、シミュレーション条件生成部１０１は、入力された数値をそのまま境界条件として抽出してもよい。ステップ５０１において曖昧さを含むシミュレート条件が入力された場合、そのシミュレート条件から具体的な数値を抽出する。
【００５４】
ここで、シミュレーション条件生成部１０１は、入力された一つのシミュレート条件から一つの境界条件のみを抽出してもよいが、通常は、一つのシミュレート条件から複数の境界条件を抽出する。例えば、シミュレータ１０２に入力可能な具体的な数値がシミュレート条件として入力された場合、シミュレーション条件生成部１０１は、入力された数値に加えて、その数値を少しずつ変更した複数の数値を複数の境界条件として抽出してもよい。曖昧さを含むシミュレート条件が入力された場合、そのシミュレート条件に含まれる（又はそのシミュレート条件に近い）複数の数値を複数の境界条件として抽出してもよい。
【００５５】
例えば、地震シミュレータに関するシミュレート条件として、一つの震源地の座標値と、その震源地において発生する地震の規模を示す一つの数値が入力された場合、シミュレーション条件生成部１０１は、その震源地を含む範囲内の複数の地点の座標値と、入力された地震の規模を示す数値を含む複数の数値と、の複数の組み合わせを複数の境界条件として抽出してもよい。あるいは、震源地のおおよその範囲と、地震のおおよその規模とがシミュレート条件として入力された場合、シミュレーション条件生成部１０１は、その範囲内の（又はその範囲を含む所定の範囲内の）複数の地点の座標値と、入力された規模の範囲内の（又はその範囲を含む所定の範囲内の）地震の規模を示す複数の数値と、の複数の組み合わせを複数の境界条件として抽出してもよい。このとき、シミュレーション条件生成部１０１は、入力されたシミュレート条件の中心を特定し、抽出された境界条件の当該シミュレート条件の中心から遠いほど小さくなるような重みを、抽出された境界条件に付与してもよい。
【００５６】
あるいは、確率的なシミュレート条件が入力された場合、シミュレーション条件生成部１０１は、その確率に従って境界条件を抽出してもよい。例えば震源地の確率分布と地震の規模を示す数値の確率分布とが入力された場合、シミュレーション条件生成部１０１は、それらの確率分布に従って、震源地の座標値と地震の規模を示す数値との複数の組み合わせを複数の境界条件として抽出してもよい。このとき、シミュレーション条件生成部１０１は、確率分布に応じた重み（例えば、分布の中心から遠いほど小さくなるような重み）を、抽出された境界条件に付与してもよい。
【００５７】
以上でシミュレーション条件生成処理が終了する。
【００５８】
図６は、本発明の第１の実施形態の物理量推定システム１００が実行するシミュレーション処理のフローチャートである。この処理は、図４のステップ４０２においてシミュレータ１０２によって実行される。
【００５９】
最初に、シミュレータ１０２は、境界条件入力処理を実行する（ステップ６０１）。具体的には、シミュレータ１０２は、ステップ５０２において抽出された境界条件を入力する。
【００６０】
次に、シミュレータ１０２は、時間発展処理を実行する（ステップ６０２）。具体的には、シミュレータ１０２は、ステップ６０１において入力された境界条件に基づいて、シミュレーションを実行することによって、各時刻における各地点の物理量を計算する。
【００６１】
ステップ５０２において複数の境界条件が抽出された場合、シミュレータ１０２は、それぞれの境界条件についてステップ６０１及びステップ６０２を実行する。
【００６２】
また、ステップ５０１において、シミュレーションの結果に基づいて別のシミュレータのシミュレート条件を生成する方法が定義された場合には、ステップ６０２が終了した後、その結果を用いてシミュレーション条件生成部１０１がステップ５０２を実行し、それによって得られた境界条件を用いてシミュレータ１０２が再びステップ６０１及びステップ６０２を実行する。
【００６３】
抽出されたすべての境界条件についてステップ６０１及びステップ６０２が実行されると、シミュレーション処理が終了する。
【００６４】
図７は、本発明の第１の実施形態の物理量推定システム１００が実行するシミュレーション処理の一例である氾濫シミュレーションの説明図である。
【００６５】
氾濫シミュレーションにおいて、シミュレーション対象の空間（すなわち河川の氾濫が発生する地域）が、所定の大きさの複数の格子（グリッド）７０１に分割され、各グリッドの各時刻における水深及び流速が計算される。図７の例では各グリッドは１辺の長さがｄの正方形である。
【００６６】
氾濫シミュレーションの処理結果７１０は、グリッドＩＤ７１１、時刻７１２、水深７１３及び流速７１４を含む。グリッドＩＤ７１１は各グリッドの識別子であり、各グリッドの座標に対応する。グリッドＩＤ７１１によって識別されるグリッドの、時刻７１２が示す時刻における浸水深及び水の流れる速度（具体的には例えば流速のｘ方向成分及びｙ方向成分）が、それぞれ水深７１３及び流速７１４として計算される。
【００６７】
このような処理結果７１０及びそれに対応するシミュレーションの境界条件がシナリオとしてシナリオＤＢ１０４に格納されてもよいが、シナリオ管理部１０３が処理結果７１０に基づいて所望のシナリオを算定し、そのシナリオをシナリオＤＢ１０４に格納してもよい（図８等参照）。なお、氾濫シミュレーションでは上記のように水深及び流速が物理量として計算されるが、他のシミュレーションでは他の物理量（例えば地震シミュレータでは各グリッドの各時刻における揺れの大きさ等）が計算され、シミュレーション処理結果として取得される。
【００６８】
なお、本実施形態において記載される「地点」及び「位置」は、緯度及び経度等の座標によって特定される点及びその位置であってもよいが、図７に示すようなグリッドＩＤによって特定されるグリッド（すなわち所定の大きさの範囲）及びその位置であってもよいし、後述する地物ＩＤ（図１１Ｂ参照）によって特定される地物及びその位置であってもよい。
【００６９】
図８は、本発明の第１の実施形態の物理量推定システム１００が実行するシナリオＤＢ構築処理のフローチャートである。この処理は、図４のステップ４０３においてシナリオ管理部１０３によって実行される。
【００７０】
最初に、シナリオ管理部１０３は、シミュレーション結果受付処理を実行する（ステップ８０１）。これによって、シミュレーション処理の結果（例えば図７に示す処理結果７１０等）が受け付けられる。
【００７１】
次に、シナリオ管理部１０３は、シナリオ算定処理を実行する（ステップ８０２）。そして、シナリオ管理部１０３は、算定されたシナリオをシナリオＤＢ１０４に格納するために、ＤＢ格納処理を実行する（ステップ８０３）。
【００７２】
以上でシナリオＤＢ構築処理が終了する。
【００７３】
上記のステップ８０２において算定され、ステップ８０３において格納されるシナリオの一例を、図９Ａ〜図９Ｃを参照して説明する。ここでは、シナリオの一例として、図７の処理結果７１０に基づいて構築される道路被災シナリオを説明する。
【００７４】
図９Ａは、本発明の第１の実施形態の物理量推定システム１００によって作成されるシナリオの説明図である。
【００７５】
図９Ａに示す複数のグリッド９０１は、図７に示すグリッド７０１と同等である。図９Ａでは、さらに、シミュレーション対象の空間に存在する道路９０２が図示される。シナリオ管理部１０３は、シミュレーション結果を用いたシナリオ算定処理（ステップ８０２）を実行することによって、道路９０２ごと及び時刻ごとの水深を計算する。
【００７６】
図９Ｂは、本発明の第１の実施形態のシナリオＤＢ１０４に含まれるシナリオ管理テーブル９１０の説明図である。
【００７７】
シナリオ管理テーブル９１０は、シナリオＩＤ９１１及び氾濫条件データ９１２を含む。シナリオＩＤ９１１は、各シナリオの識別情報（ＩＤ）である。氾濫条件データ９１２は、シナリオを構築する元になったシミュレーション処理のシミュレート条件（又は境界条件）に相当し、例えば、河川の堤防が決壊した、降雨量が極端に大きくなった、などの氾濫の初期条件に関するデータである。氾濫シミュレーション以外のシミュレーションに基づくシナリオの場合は、その元になったシミュレーションの初期条件が氾濫条件データ９１２の代わりに登録される。
【００７８】
図９Ｃは、本発明の第１の実施形態のシナリオＤＢ１０４に含まれるシナリオ物理量管理テーブル９２０の説明図である。
【００７９】
シナリオ物理量管理テーブル９２０は、シナリオＩＤ９２１、道路ＩＤ９２２、時刻９２３、水深９２４及び流速９２５を含む。
【００８０】
シナリオＩＤ９２１は、各シナリオの識別情報であり、図９ＢのシナリオＩＤ９１１に対応する。
【００８１】
道路ＩＤ９２２は、シミュレーション対象の空間内の各道路９０２の識別情報であり、各道路９０２の座標値を含んでもよい。
【００８２】
道路ＩＤ９２２によって識別される道路の、時刻９２３が示す時刻における浸水深及び水の流れる速度（具体的には例えば流速のｘ方向成分及びｙ方向成分）が、それぞれ水深９２４及び流速９２５として格納される。
【００８３】
シナリオ管理部１０３は、図８のステップ８０２において、図７の処理結果７１０に基づいて、シナリオ物理量管理テーブル９２０に登録される情報を算定する。これらの情報は、具体的には例えば次のように算定される。物理量推定システム１００は、少なくとも道路の位置に関する情報を含む地図情報（図示省略）を保持する。シナリオ管理部１０３は、地図情報に基づいて、各グリッド９０１（すなわち図７のグリッド７０１）と各道路９０２との対応関係、すなわち、各道路９０２がどのグリッド内にあるかを特定する。そして、シナリオ管理部１０３は、特定されたグリッド９０１の水深７１３及び流速７１４に基づいて、そのグリッドに対応する道路９０２の水深９２４及び流速９２５を算定する。
【００８４】
図９Ｃに示すシナリオ物理量管理テーブル９２０は、一つのシナリオの、一つの道路に関する、一つの時刻における水深及び流速を示す。ステップ８０２では、一つのシナリオについて、シミュレーション対象の地域内の各道路の各時刻における水深及び流速が算定され、それぞれの値を含む複数のシナリオ物理量管理テーブル９２０がシナリオＤＢ１０４に格納される。
【００８５】
上記のように、図９Ｃの例では、道路ＩＤ９２２及びそれによって識別される道路における水深９２４等がシナリオ物理量管理テーブル９２０に格納される。これは、後述するように、アプリケーションシステム１４０が誘導情報生成処理３０３を実行するために各道路が浸水しているか否かを示す情報を利用するからである。このようなアプリケーションシステム１４０の処理及びそれに適するシナリオ物理量管理テーブル９２０は一例であり、シナリオ管理部１０３は、他の処理に適するシナリオ物理量管理テーブル９２０を生成することもできる。例えば、アプリケーションシステム１４０が各グリッドにおける物理量を利用した処理を実行する場合には、シナリオ物理量管理テーブル９２０は道路ＩＤ９２２の代わりにグリッドＩＤを含んでもよい。アプリケーションシステム１４０が道路以外の任意の地物（例えば河川又は建物等）における物理量を利用した処理を実行する場合には、シナリオ物理量管理テーブル９２０は道路ＩＤ９２２の代わりに任意の地物ＩＤを含んでもよい。
【００８６】
図１０は、本発明の第１の実施形態の物理量推定システム１００が実行するシナリオクラスタ生成処理のフローチャートである。この処理は、図４のステップ４０４においてシナリオ管理部１０３によって実行される。
【００８７】
シナリオ管理部１０３は、最初に、シナリオ特徴量算定処理を実行し（ステップ１００１）、続いて、得られたシナリオの特徴量に基づいてクラスタリングを実行する（ステップ１００２）。
【００８８】
シナリオ特徴量算定処理（ステップ１００１）は、シナリオのデータからシナリオの特徴を表す量（特徴量）を抽出する処理である。抽出された特徴量に基づいてクラスタリングを行うことによって、クラスタリングの計算量を削減することができる。
【００８９】
図９Ｃに例示するように、一つのシナリオ（すなわち一つのシナリオＩＤによって識別されるシナリオ）は、各地点の各時刻における物理量（図９Ｃの例では、各道路の各時刻における水深及び流速）を含む。すなわち、各シナリオは、各地点の各時刻における物理量を要素として含む多変量ベクトルによって表現することができる。シナリオ管理部１０３は、主成分分析又は独立成分分析といった公知の方法によって、シナリオを表現するベクトルの次元数を削減することができる。
【００９０】
あるいは、シナリオ管理部１０３は、クラスタリングの結果への寄与が大きいと考えられる特定の地点の物理量を特徴量として選択してもよい。例えば、クラスタリング対象の複数のシナリオのいずれにおいても、ある地点（第１地点）の物理量（例えば水深）はほとんど同じであるが、別のある地点（第２地点）の物理量はシナリオごとに大きく異なる場合には、第２地点の物理量はクラスタリングの結果に大きく寄与するが、第１地点の物理量の寄与は比較的小さいと考えられる。このような場合、シナリオ管理部１０３は、第２地点の物理量を含み、第１地点の物理量を含まない特徴量を算定してもよい。
【００９１】
このように、シナリオの特徴量を抽出し、その特徴量に基づいてクラスタリングを行うことによって、クラスタリングの結果への寄与が小さい物理量をクラスタリングの計算から除外することができ、これによって計算量を削減できる。ただし、計算量を削減する必要がない場合には、シナリオ管理部１０３はステップ１００１を実行せずに（すなわち各シナリオを表現するベクトルをそのまま特徴量として使用することによって）ステップ１００２のクラスタリングを実行してもよい。
【００９２】
クラスタリング（ステップ１００２）は、ステップ１００１において抽出された特徴量を用いて、ｋ−ｍｅａｎｓ法、混合ガウス分布に基づく方法、階層クラスタリングなど、公知の方法によって実行することができる。クラスタリングによって、各シナリオクラスタに一つの又は互いに比較的類似する二つ以上のシナリオが含まれるように、複数のシナリオが複数のシナリオクラスタに分類される。
【００９３】
次に、シナリオ管理部１０３は、クラスタ統計量算出処理を実行する（ステップ１００３）。具体的には、シナリオ管理部１０３は、各シナリオクラスタに含まれるシナリオの（例えば各地点の各時刻の）物理量の平均値、分散値及び期待値といった統計量を算出する。
【００９４】
次に、シナリオ管理部１０３は、重要スポット抽出処理を実行する（ステップ１００４）。ここで、重要スポットとは、シナリオクラスタ内で特に値のばらつきが大きい地点（又は地点及び時刻の組）である。上記のように、クラスタリングによって二つ以上のシナリオが一つのシナリオクラスタに分類された場合、それらのシナリオの特徴量は互いに比較的類似している。すなわち、それらの二つ以上のシナリオの各地点の各時刻における物理量は互いに類似していることが多い。しかし、特定の地点の（又は、特定の地点の特定の時刻における）物理量は、一つのクラスタ内のシナリオ間でばらつきが大きい場合がある。シナリオ管理部１０３は、上記のような物理量のばらつきの大きい地点（又は地点及び時刻の組み合わせ）を、重要スポットとして抽出する。
【００９５】
例えば氾濫シミュレーションの場合、一つのシナリオクラスタに含まれるシナリオにおける堤防の決壊地点及び浸水域等はおおむね類似すると考えられる。しかし、シミュレーション対象の地域のうち、比較的標高の高い地点では、シナリオクラスタ内のいずれのシナリオでも全く浸水せず、比較的標高の低い地点ではすべてのシナリオで浸水し、中程度の標高の地点に関しては、シナリオクラスタに含まれる一部のシナリオでは浸水し、残りのシナリオでは浸水しない、という場合がある。
【００９６】
クラスタ内でのある地点の物理量（例えば浸水深）のばらつきが大きいことは、シミュレーション条件のわずかな相違に応じて、その地点で計算される物理量が大きく異なる可能性があること、すなわち、その地点の物理量を精度よく推定するためには、他の地点より精密な推定を行う必要があることを意味する。このため、シミュレーション結果のばらつきの大きい地点（上記の例では中程度の標高の地点）が、精密な推定を必要とする重要スポットとして抽出される。
【００９７】
具体的には、シナリオ管理部１０３は、シナリオクラスタ内のシナリオの物理量のばらつきを示す指標（例えば分散値）が所定の条件を満たす（例えば所定の値を超える）地点及び時刻を、重要スポットとして抽出してもよい。後述するシナリオ合成処理において、重要スポットの物理量は、実際にセンサ１２０等によって計測されたデータを用いてより精密に推定される。
【００９８】
シナリオクラスタにおいて、ある地点の物理量のばらつきが時刻にかかわらず大きい場合もあるが、特定地点の特定の時刻の物理量のみのばらつきが大きい場合もあり、後者の場合には、当該特定の地点及び特定の時刻の組み合わせが重要スポットとして抽出される。以下、重要スポットとして地点及び時刻の組み合わせが抽出された場合を例として説明するが、地点のみが抽出された場合（言い換えるとある地点と全時刻の範囲との組み合わせが重要スポットとして抽出された場合）にも以下の説明を適用することができる。
【００９９】
以上でシナリオクラスタ生成処理が終了する。シナリオ管理部１０３は、図１１Ａ〜図１１Ｄを参照して後述するように、シナリオクラスタ生成処理の結果をシナリオＤＢ１０４に格納する。
【０１００】
図１１Ａは、本発明の第１の実施形態のシナリオＤＢ１０４に含まれるシナリオクラスタリングテーブル１１００の説明図である。
【０１０１】
シナリオクラスタリングテーブル１１００は、シナリオＩＤ１１０１及びシナリオクラスタＩＤ１１０２を含む。シナリオＩＤ１１０１は、各シナリオの識別情報であり、図９ＢのシナリオＩＤ９１１等に対応する。シナリオクラスタＩＤ１１０２は、各シナリオが属するシナリオクラスタの識別情報であり、これがクラスタリング（ステップ１００２）の結果に相当する。
【０１０２】
図１１Ｂは、本発明の第１の実施形態のシナリオＤＢ１０４に含まれるシナリオクラスタ管理テーブル１１１０の説明図である。
【０１０３】
シナリオクラスタ管理テーブル１１１０は、シナリオクラスタＩＤ１１１１、地物ＩＤ１１１２、クラスタ内平均値１１１３及びクラスタ内分散値１１１４を含む。
【０１０４】
シナリオクラスタＩＤ１１１１は、各シナリオクラスタの識別情報であり、図１１ＡのシナリオクラスタＩＤ１１０２に対応する。
【０１０５】
地物ＩＤ１１１２は、各シナリオクラスタに属するシナリオにおいて、物理量と対応付けられる地物の識別情報である。例えば、図９Ｃの例のように水深等の物理量が道路に対応付けられる場合、地物ＩＤ１１１２は道路ＩＤである。また、シナリオにおいて物理量がグリッドに対応付けられる場合（例えば図７参照）には、地物ＩＤ１１１２はグリッドＩＤである。あるいは、地物ＩＤ１１１２は、建物の識別情報等であってもよい。
【０１０６】
クラスタ内平均値１１１３及びクラスタ内分散値１１１４は、それぞれ、各シナリオクラスタに属するシナリオの各地物に対応付けられた物理量の平均値及び分散値である。これらは、クラスタ統計量算出処理（ステップ１００３）によって算出された値である。
【０１０７】
なお、図１１Ｂは、シナリオクラスタ管理テーブル１１１０の代表的な形態を例示したものであり、実際には例えば図１１Ｃに示すように種々のバリエーションがあり得る。
【０１０８】
図１１Ｃは、本発明の第１の実施形態のシナリオＤＢ１０４に含まれるシナリオクラスタ管理テーブル１１１０の具体例の説明図である。
【０１０９】
図１１Ｃに示すシナリオクラスタ管理テーブル１１２０は、図９Ａ〜図９Ｃに示すシナリオをクラスタリングすることによって得られたシナリオクラスタに関するシナリオクラスタ管理テーブル１１１０の具体例であり、シナリオクラスタＩＤ１１２１、道路ＩＤ１１２２、時刻１１２３、平均水深１１２４及び平均流速１１２５を含む。
【０１１０】
シナリオクラスタＩＤ１１２１は、各シナリオクラスタの識別情報であり、図１１ＢのシナリオクラスタＩＤ１１１１に対応する。
【０１１１】
道路ＩＤ１１２２は、各道路９０２の識別情報であり、図１１Ｂの地物ＩＤ１１１２に対応する。
【０１１２】
時刻１１２３は、各シナリオに含まれる物理量に関連付けられる時刻であり、図９Ｃの時刻９２３に対応する。
【０１１３】
平均水深１１２４及び平均流速１１２５は、各道路９０２に関連付けられる各時刻の水深及び流速のシナリオクラスタ内の平均値であり、図１１Ｂのクラスタ内平均値１１１３に対応する。これらの値は、図９Ｃの水深９２４及び流速９２５から算出される。
【０１１４】
シナリオクラスタ管理テーブル１１２０は、さらに、水深及び流速の分散値を含んでもよい。
【０１１５】
図１１Ｄは、本発明の第１の実施形態のシナリオＤＢ１０４に含まれる重要スポット管理テーブル１１３０の説明図である。
【０１１６】
重要スポット管理テーブル１１３０は、シナリオクラスタＩＤ１１３１、道路ＩＤ１１３２及び時刻１１３３を含む。
【０１１７】
シナリオクラスタＩＤ１１３１は、各シナリオクラスタの識別情報であり、図１１ＡのシナリオクラスタＩＤ１１０２に対応する。
【０１１８】
道路ＩＤ１１３２及び時刻１１３３は、重要スポットとして抽出された地点及び時刻の組を識別する情報である。例えば、ある道路におけるある時刻の物理量（例えば水深又は流速）のシナリオクラスタ内の分散値が所定の値より大きい場合、その道路及び時刻の組が重要スポットとして抽出される。複数の重要スポットが抽出された場合には、複数の重要スポットを識別する道路ＩＤ及び時刻の複数の組が重要スポット管理テーブル１１３０に登録される。
【０１１９】
以上で事前処理（図３Ａ）が終了する。続いて、リアルタイム処理（図３Ｂ）の詳細について説明する。
【０１２０】
図１２は、本発明の第１の実施形態の物理量推定システム１００が実行するシナリオ合成処理のフローチャートである。この処理は、図３Ｂのステップ３０２において実行される。
【０１２１】
最初に、物理量推定システム１００の入力受付部１０５が実物理量入力処理を実行する（ステップ１２０１）。実物理量入力処理の詳細については図１３等を参照して後述する。
【０１２２】
次に、物理量推定システム１００のシナリオ合成部１０６がシナリオ合成処理を実行する（ステップ１２０２）。ステップ１２０２のシナリオ合成処理の詳細については図１５等を参照して後述する。合成されたシナリオは合成シナリオＤＢ１０７に格納される。
【０１２３】
次に、アプリケーションシステム１４０が誘導情報生成処理を実行する（ステップ１２０３）。誘導情報生成処理の詳細については図１７を参照して後述する。
【０１２４】
図１３は、本発明の第１の実施形態の物理量推定システム１００が実行する実物理量入力処理のフローチャートである。この処理は、図１２のステップ１２０１において実行される。
【０１２５】
入力受付部１０５は、目視情報受付処理（ステップ１３０１）及びセンサ情報受付処理（ステップ１３０２）を実行する。入力受付部１０５は、ステップ１３０１及び１３０２のどちらを先に実行してもよい。
【０１２６】
目視情報受付処理は、ユーザ等が現況入力システム１３０を介して入力したデータを、入力受付部１０５が、ネットワーク１５０Ｂ及びＩ／Ｆ２０５を介して受け付ける処理である。入力されるデータの例を図１４Ａに示す。
【０１２７】
センサ情報受付処理は、センサ１２０から送信されたデータを、入力受付部１０５が、ネットワーク１５０Ｂ及びＩ／Ｆ２０５を介して受け付ける処理である。センサ１２０から送信されるデータの例を図１４Ｂに示す。
【０１２８】
図１４Ａは、本発明の第１の実施形態の物理量推定システム１００が現況入力システム１３０から受け付けるデータの説明図である。
【０１２９】
物理量推定システム１００の入力受付部１０５が現況入力システム１３０から受け付ける目視情報のデータ１４００は、例えば、地点１４０１、時刻１４０２及び浸水程度１４０３を含む。
【０１３０】
地点１４０１は、ユーザ等が目視（又は体感）によって物理量を観測した地点であり、座標値等であってもよいが、町名及び番地のような情報であってもよい。町名及び番地が入力された場合、物理量推定システム１００が地図情報（図示省略）を参照してそれを座標値（又は座標値の範囲）に変換することができる。
【０１３１】
時刻１４０２は、ユーザ等が目視（又は体感）によって物理量を観測した時刻である。
【０１３２】
浸水程度１４０３は、ユーザ等が目視（又は体感）によって観測した、河川の氾濫等による浸水の程度（すなわち浸水深）を示す情報である。浸水程度１４０３は、例えば「５０センチメートル」のような具体的な数値であってもよいが、ユーザ等の目視による観測の精度は低いと考えられるため、例えば「浸水なし」「床下浸水」又は「床上浸水」の３段階のいずれかを示す情報であってもよい。
【０１３３】
上記のデータ１４００は一例であり、入力受付部１０５は、例えば地震の揺れの程度を示す情報等、ユーザ等の目視（又は体感）によって観測され、現況入力システム１３０から送信された任意の物理量に関するデータを受け付けることができる。例えば地震の場合、目視情報のデータ１４００は、浸水程度１４０３の代わりに、ユーザ等が体感した震度を含んでもよい。
【０１３４】
図１４Ｂは、本発明の第１の実施形態の物理量推定システム１００がセンサ１２０から受け付けるデータの説明図である。
【０１３５】
図１４Ｂには、物理量推定システム１００の入力受付部１０５がセンサ１２０から受け付けるデータの一例として、浸水センサ（すなわち液面検知センサ）から送付されるデータ１４１０を示す。データ１４１０は、センサＩＤ１４１１、時刻１４１２及び水深１４１３を含む。
【０１３６】
センサＩＤ１４１１は、データ１４１０を送信したセンサ１２０の識別情報である。物理量推定システム１００は各センサ１２０が設置された位置を特定する情報（座標値等）を保持しているため、センサＩＤ１４１１に基づいてセンシングが行われた地点を特定することができる。
【０１３７】
時刻１４１２は、センサ１２０によるセンシングが行われた時刻である。
【０１３８】
水深１４１３は、センサ１２０によってセンシングされた水深である。
【０１３９】
上記のデータ１４１０は一例であり、入力受付部１０５は、例えば震度等、センサ１２０から送信された任意の物理量に関するデータを受け付けることができる。例えば震度が送信される場合、データ１４１０は、水深１４１３の代わりに震度を含む。
【０１４０】
図１５は、本発明の第１の実施形態の物理量推定システム１００が実行するシナリオ合成処理のフローチャートである。この処理は、図１２のステップ１２０２において実行される。
【０１４１】
最初に、物理量推定システム１００のシナリオ合成部１０６は、現況シナリオ特徴量算定処理を実行する（ステップ１５０１）。具体的には、シナリオ合成部１０６は、実物理量入力処理（ステップ１２０１）において入力された物理量（すなわち実際にセンサ１２０によって計測された物理量及びユーザ等によって入力された物理量）の特徴量を計算する。例えば、シナリオ合成部１０６は、実物理量入力処理において入力された物理量を対象として、図１０のステップ１００１と同様の処理を行ってもよい。
【０１４２】
次に、シナリオ合成部１０６は、シナリオクラスタ判定処理を実行する（ステップ１５０２）。具体的には、シナリオ合成部１０６は、ステップ１５０１において抽出された特徴量と、シナリオＤＢ１０４に格納されたシナリオクラスタの特徴量との類似度を計算し、最も類似度の高いシナリオクラスタを選択する。ここで、シナリオクラスタの特徴量は、シナリオクラスタに含まれるシナリオの特徴量の代表値、例えば、図１０のステップ１００１において抽出された特徴量の各要素をシナリオクラスタ内で平均することによって得られる値を要素として含むベクトルであってもよい。
【０１４３】
なお、一般に、センサ１２０又は現況入力システム１３０による物理量の観測が行われる地点及び時刻の組み合わせの数は、シミュレーションによって物理量が計算される地点及び時刻の組み合わせの数よりずっと少ない。このため、シナリオ合成部１０６は、実物理量入力処理（ステップ１２０１）において入力された物理量が観測された地点及び時刻に対応する地点及び時刻の各シナリオの物理量に基づいてシナリオクラスタの特徴量を計算し、ステップ１５０２においてそのシナリオクラスタの特徴量と現況シナリオ特徴量との類似度を計算してもよい。
【０１４４】
ステップ１５０１において抽出された特徴量及びシナリオクラスタの特徴量はいずれもベクトルであるため、それらのベクトル間の距離の近さを類似度の高さとして扱うことができる。なお、ここで用いる距離としては、マハラノビス距離など、多変量解析の分野で知られている距離の定義を用いてもよい。
【０１４５】
次に、シナリオ合成部１０６は、重要スポット情報推定処理を実行する（ステップ１５０３）。この処理は、重要スポットの物理量を、ステップ１５０２において選択されたクラスタに含まれる全シナリオの物理量と、実物理量入力処理（ステップ１２０１）において入力された物理量と、に基づいて推定する処理である。具体的には、シナリオ合成部１０６は、ステップ１５０２において選択されたシナリオクラスタに含まれる全シナリオの重要スポットの物理量について、それらの物理量と、重要スポットの地点及び時刻に対応するステップ１２０１において入力された物理量と、の類似度に基づく重み付け演算（具体的には例えば重み付け平均）を行うことによって、重要スポットにおける物理量の期待値を計算する。
【０１４６】
ここで、類似度は、重要スポットにおけるシナリオの物理量とそれに対応するステップ１２０１において入力された物理量との距離（差）の逆数であってもよい。例えば、あるシナリオ（ここではシナリオＡと記載する）の重要スポットの物理量とステップ１２０１において入力された物理量との距離を、全シナリオの重要スポットの物理量とステップ１２０１において入力された物理量との距離の合計値で除算した値を、シナリオＡの重要スポットの物理量の重みとして計算してもよい。シナリオ合成部１０６は、同様の重みをシナリオクラスタ内の全シナリオの重要スポットの物理量について計算し、それによって重み付けされた物理量の平均値を、当該シナリオクラスタの重要スポットの物理量の期待値として計算することができる。
【０１４７】
次に、シナリオ合成部１０６は、物理量推定処理を実行する（ステップ１５０４）。具体的には、シナリオ合成部１０６は、重要スポット以外の地点及び時刻（以下、非重要スポットとも記載する）については、それぞれの地点及び時刻における、選択されたシナリオクラスタの代表値を、それぞれの地点及び時刻における物理量として推定し、重要スポットについては、ステップ１５０３における重み付け演算によって計算された期待値を、当該重要スポットの地点及び時刻における物理量として推定し、それらの推定した物理量を含むシナリオを新たに生成する。このようにして生成されたシナリオは、上記のように、重み付け演算等によって複数のシナリオを合成することによって得られた物理量を含むため、合成シナリオと記載される。シナリオ合成部１０６は、生成した合成シナリオを合成シナリオＤＢ１０７に格納する。
【０１４８】
ここで、重要スポット以外の地点及び時刻における、選択されたシナリオクラスタの代表値とは、その地点及び時刻における、選択されたシナリオクラスタに含まれる複数のシナリオの物理量を代表する値であり、具体的にはそれらの複数のシナリオのうちいずれか一つの物理量であってもよいし、複数のシナリオの物理量の統計量であってもよい。例えば、代表値は、選択されたシナリオクラスタ内の全シナリオの物理量の平均値であってもよいし、選択されたシナリオクラスタ内の全シナリオのうち、ステップ１５０１において抽出された特徴量との類似度が最も高いシナリオの物理量であってもよいし、選択されたシナリオクラスタ内の全シナリオのうち、ステップ１５０１において抽出された特徴量との類似度が所定の値以上の複数のシナリオの物理量の平均値であってもよい。また、統計量として、平均値以外の値（例えば中央値等）を使用してもよい。
【０１４９】
なお、上記の選択されたシナリオクラスタ内の全シナリオの物理量の平均値とは、シナリオクラスタ内の全シナリオの物理量に同一の重みを付けて重み付け平均を行った結果と言い換えることもできる。ステップ１５０１において抽出された特徴量との類似度が最も高いシナリオの物理量とは、ステップ１５０１において抽出された特徴量との類似度が最も高いシナリオの物理量に０以外の重みを付け、他の物理量に０の重みを付けて重み付け平均を行った結果と言い換えることもできる。ステップ１５０１において抽出された特徴量との類似度が所定の値以上の複数のシナリオの物理量の平均値とは、ステップ１５０１において抽出された特徴量との類似度が所定の値以上の複数のシナリオの物理量に０以外の同一の重みを付け、他の物理量に０の重みを付けて重み付け平均を行った結果と言い換えることもできる。このように、非重要スポットにおける物理量の推定も、重み付け演算によって求めることができるが、重みの値を決定する方法が重要スポットにおけるものとは異なる。
【０１５０】
シナリオ合成部１０６は、非重要スポットにおける物理量を推定する計算量をさらに削減するため、時刻の間引き又はシナリオの間引きを行ってもよい。前者を実行する場合、シナリオ合成部１０６は、非重要スポットのすべての時刻における代表値を計算する代わりに、一部の時刻（例えば所定の間隔で間引きされた一部の時刻）の代表値を計算し、その他の時刻の物理量は、計算された代表値の値を補完することによって推定する。後者を実行する場合、シナリオ合成部１０６は、上記のような類似度によらずに、シナリオクラスタに含まれる複数のシナリオの一部を選択し、選択されたシナリオの非重要スポットにおける物理量の平均値を代表値として計算してもよい。時刻の間引き及びシナリオの間引きを組み合わせてもよい。
【０１５１】
重要スポット以外の地点及び時刻における物理量のシナリオクラスタ内でのばらつきは小さいため、どのシナリオの物理量を代表値として選択しても、あるいは平均値を代表値として計算しても、推定の結果は大きく変わらない。これに対して、重要スポットにおける物理量のシナリオクラスタ内でのばらつきは大きいため、ステップ１２０１において入力された物理量（すなわちリアルタイムで実測された物理量）に基づく重み付け平均によって物理量が推定される。これによって、実測された物理量に整合する高精度の物理量推定を行うことができる。
【０１５２】
図１６は、本発明の第１の実施形態の物理量推定システム１００が実行するシナリオ合成処理の概念の説明図である。
【０１５３】
例えば、道路１６０１にいくつかのセンサ１６０２が設置されている場合に、シナリオ合成部１０６は、それらのセンサ１６０２から取得された物理量の特徴量を現況シナリオ特徴量として計算する（ステップ１５０１）。
【０１５４】
そして、シナリオ合成部１０６は、現況シナリオ特徴量との類似度が最も高い特徴量を有するシナリオクラスタを選択する（ステップ１５０２）。
【０１５５】
続いて、シナリオ合成部１０６は、そのシナリオクラスタの重要スポットについては、センサ１６０２から取得された物理量に基づく、選択されたシナリオクラスタ内の全シナリオ１６１０の物理量の重み付け平均によって期待値を計算する（ステップ１５０３）。
【０１５６】
一方、重要スポット以外の地点及び時刻については、シナリオ合成部１０６は、現況シナリオに最も類似するシナリオの物理量（又は、選択されたシナリオクラスタ内の全シナリオ１６１０の物理量の、重み付けをしない平均値）を物理量の期待値として取得する。
【０１５７】
そして、シナリオ合成部１０６は、重要スポット及びそれ以外の地点及び時刻のそれぞれの物理量の期待値（すなわち推定された物理量）を合成することによって合成シナリオ１６２０を生成し、それを合成シナリオＤＢ１０７に格納する（ステップ１５０４）。
【０１５８】
物理量推定システム１００のシナリオ配信部１０８は、必要に応じて（例えばアプリケーションシステム１４０からの要求に応じて）、合成シナリオＤＢ１０７に格納された合成シナリオ（例えば上記の合成シナリオ１６２０）をアプリケーションシステム１４０に配信する。
【０１５９】
図１７は、本発明の第１の実施形態のアプリケーションシステム１４０が実行する誘導情報生成処理のフローチャートである。
【０１６０】
図１７には、アプリケーションシステム１４０が提供するアプリケーションの一例として、河川等が氾濫した（又はこれから氾濫すると予想される）時に、避難等のためにユーザを誘導するアプリケーションを示す。このため、図１７に示す例において、アプリケーションシステム１４０は、公知のナビゲーションシステム等が有するものと同様の道路ネットワークデータ（図示省略）及びそれに基づく経路探索機能を有する。
【０１６１】
最初に、アプリケーションシステム１４０は、経路コスト算定処理を実行する（ステップ１７０１）。このとき、アプリケーションシステム１４０は、シナリオ配信部１０８から配信された合成シナリオ（例えば図１６の合成シナリオ１６２０）を参照する。図１６に例示する合成シナリオ１６２０は、氾濫シミュレーションの結果を含むため、推定された物理量として各道路の浸水深の値を含む。そこで、アプリケーションシステム１４０は、例えば、浸水深の値が０より大きい（又は所定の値より大きい）道路を、通行不可能な道路と判定して、その道路が経路探索によって選択されないように、道路ネットワークデータを変更する。
【０１６２】
例えば、アプリケーションシステム１４０は、通行不可能と判定された道路に、通行止め等の、通行を禁止する情報を設定してもよい。あるいは、アプリケーションシステム１４０は、通行不可能と判定された道路のコストを、その道路が経路探索結果として選択されない程度に十分に大きい値に変更してもよい。
【０１６３】
そして、アプリケーションシステム１４０は、上記のように変更された道路ネットワークデータを参照して、例えば、ユーザの現在地から所定の避難場所までの経路のコストを計算する。
【０１６４】
次に、アプリケーションシステム１４０は、ステップ１７０１において計算された経路コストに基づいて、経路探索処理を実行し、その結果を出力する（ステップ１７０２）。例えば、アプリケーションシステム１４０は、経路探索処理の結果を、ネットワーク１５０Ｄを介して、ユーザが保持する移動体端末３０４に送信してもよい。
【０１６５】
以上で誘導情報生成処理が終了する。なお、上記のような誘導情報生成処理は一例であり、アプリケーションシステム１４０は、合成シナリオを利用して種々の処理を実行することができる。
【０１６６】
以上のように、本発明の第１の実施形態の物理量推定システム１００は、種々の条件に基づいて物理量を推定するシミュレーションを実行し、各シミュレーション条件とそれに基づくシミュレーション結果との組をシナリオとして保持し、そのようにして得られた多数のシナリオをクラスタリングし、各クラスタ内の物理量のばらつきの大きい地点及び時刻を重要スポットとして特定しておく。そして、物理量推定システム１００は、センサ等によって観測された物理量を含むリアルタイムデータを取得すると、取得したリアルタイムデータの特徴量と類似するクラスタを特定し、そのクラスタの重要スポットにおける物理量をリアルタイムデータに基づいて推定し、重要スポットにおける推定された物理量と、重要スポット以外の地点及び時刻におけるそのクラスタの物理量の代表値（例えば平均値）と、からなる合成シナリオを生成する。この合成シナリオに含まれる物理量が、現在までにセンサ１２０等から取得された物理量に基づく現在及び未来の各時刻における各地点の推定された物理量である。
【０１６７】
上記の処理のうち、物理量推定システム１００がセンサ１２０等から取得したリアルタイムデータを用いなければ実行できないのは、クラスタの特定及び重要スポットにおける物理量の推定のみであり、シミュレーションの実行から重要スポットの特定までは、リアルタイムデータを取得する前に実行しておくことができる。これによって、リアルタイムデータの取得後に実行すべき計算の量を削減することができるため、高速に物理量を推定し、その物理量に基づいて、迅速に避難情報等を提供することができる。また、クラスタ内の物理量のばらつきが大きい重要スポットについては、センサ１２０等から取得された物理量に基づく推定を行うことによって、より尤もらしい推定結果を得ることができる。
【０１６８】
＜第２の実施形態＞
図１８は、本発明の第２の実施形態の計算機システムの構成を示す機能ブロック図である。
【０１６９】
本実施形態の計算機システムは、物理量推定システム１８００、シミュレータ１８１０、シナリオ例入力システム１１０、センサ１２０、現況入力システム１３０及びアプリケーションシステム１４０によって構成される。シミュレータ１８１０はネットワーク１５０Ｅを介して、シナリオ例入力システム１１０はネットワーク１５０Ａを介して、センサ１２０及び現況入力システム１３０はネットワーク１５０Ｂを介して、アプリケーションシステム１４０はネットワーク１５０Ｃを介して、それぞれ物理量推定システム１００に接続される。
【０１７０】
本実施形態の計算機システムのうち、物理量推定システム１８００、シミュレータ１８１０及びネットワーク１５０Ｅ以外の構成は、図１において同一の符号を付された構成と同一の機能を有するため、それらの説明は省略する。
【０１７１】
物理量推定システム１８００は、シミュレータ１０２を有さず、シミュレータ通信部１８０１及びシミュレータリスト１８０２を有する点を除いて、第１の実施形態の物理量推定システム１００と同一である。上記以外の物理量推定システム１８００の構成は、図１において同一の符号を付された構成と同一の機能を有するため、それらの説明は省略する。
【０１７２】
シミュレータ通信部１８０１は、シミュレータ１８１０との通信を実行する機能を有する。その詳細については、図２０を参照して後述する。
【０１７３】
シミュレータリスト１８０２は、所望のシミュレーションを実行するシミュレータを特定する情報を含む。具体的には、シミュレータリスト１８０２は、例えば、シミュレーション条件生成部によって生成された境界条件に基づくシミュレーションを実行できるシミュレータ１８１０を特定する情報（例えばそのシミュレータ１８１０のアドレス）等のリストであってもよい。
【０１７４】
シミュレータ１８１０は、第１の実施形態のシミュレータ１０２と同様の機能を有する。ただし、第１の実施形態のシミュレータ１０２は、物理量推定システム１００のプロセッサ２０１が物理量推定システム１００の記憶装置２０６に格納されたプログラムを実行することによって実現さたが、本実施形態のシミュレータ１８１０は、物理量推定システム１００以外の計算機によって実現される。本実施形態の計算機システムは、それぞれが異なる種類のシミュレーションを実行する複数のシミュレータ１８１０を有してもよい。
【０１７５】
ネットワーク１５０Ｅは、シミュレータ１８１０と物理量推定システム１８００との間の通信を可能にするものであれば、どのようなものであってもよい。例えば、ネットワーク１５０Ｅは、ネットワーク１５０Ａ、１５０Ｂ及び１５０Ｃのいずれかと同種のものであってもよい。また、図１には三つのネットワーク１５０Ａ、１５０Ｂ、１５０Ｃ及び１５０Ｅを示したが、これらが一つのネットワークによって実現されてもよい。以下の説明においてネットワーク１５０が参照された場合、それは、ネットワーク１５０Ａ、１５０Ｂ、１５０Ｃ及び１５０Ｅのいずれか、又は、ネットワーク１５０Ａ、１５０Ｂ、１５０Ｃ及び１５０Ｅを実現する一つのネットワークである。
【０１７６】
物理量推定システム１８００のハードウェア構成は、図２に示した第１の実施形態の物理量推定システム１００と同様であるため、図示を省略する。ただし、物理量推定システム１８００の記憶装置２０６には、シミュレータ１０２が格納されず、シミュレータ通信部１８０１を実現するプログラム及びシミュレータリスト１８０２が格納される。以下の説明においてシミュレータ通信部１８０１が実行する処理は、実際にはプロセッサ２０１によって実行される。
【０１７７】
次に、本実施形態の物理量推定システム１８００が実行する処理について説明する。
【０１７８】
物理量推定システム１８００は、事前処理（図１９Ａ）及びリアルタイム処理（図１９Ｂ）を実行する。これらの処理は、第１の実施形態の事前処理（図３Ａ）及びリアルタイム処理（図３Ｂ）と同様の目的で実行される。
【０１７９】
図１９Ａは、本発明の第２の実施形態の物理量推定システム１８００が実行する事前処理のフローチャートである。
【０１８０】
物理量推定システム１８００は、シナリオＤＢ構築処理を実行する（ステップ１９０１）。これによってシナリオＤＢ１０４が構築される。シナリオＤＢ構築処理の詳細については図２０等を参照して後述する。
【０１８１】
図１９Ｂは、本発明の第２の実施形態の物理量推定システム１８００が実行するリアルタイム処理のフローチャートである。
【０１８２】
物理量推定システム１８００は、シナリオＤＢ１０４を参照してシナリオ合成処理を実行し（ステップ１９０２）、その結果を合成シナリオＤＢ１０７に格納する。本実施形態のシナリオ合成処理は、第１の実施形態と同様である。
【０１８３】
アプリケーションシステム１４０は、合成シナリオＤＢ１０７を参照して種々のアプリケーションを実行することができる。例えば、アプリケーションシステム１４０は、合成シナリオＤＢ１０７を参照して被害推定表示処理を実行し（ステップ１９０３）、生成した誘導情報を、ネットワーク１５０Ｆを介して被害状況可視化端末１９０４に送信してもよい。被害推定表示処理の詳細については図２４を参照して後述する。
【０１８４】
ここで、被害状況可視化端末１９０４は、アプリケーションシステム１４０のユーザが所持するパーソナルコンピュータ等の端末装置である。ネットワーク１５０Ｆは、アプリケーションシステム１４０と移動体端末３０４との間の通信を可能にするものであれば、どのようなものであってもよい。
【０１８５】
ここで、事前処理（図１９Ａ）の詳細について説明する。
【０１８６】
図２０は、本発明の第２の実施形態の物理量推定システム１８００が実行するシナリオＤＢ構築処理のフローチャートである。この処理は、図１９Ａのステップ１９０１において実行される。
【０１８７】
最初に、シミュレーション条件生成部１０１がシミュレーション条件生成処理を実行する（ステップ４０１）。この処理は、第１の実施形態のシミュレーション条件生成部１０１が実行するものと同様であるため（図５参照）、ここでは説明を省略する。
【０１８８】
次に、シミュレータ通信部１８０１が、必要シミュレータ選択処理（ステップ２００１）及びシミュレーション実行依頼処理（ステップ２００２）を実行する。
【０１８９】
必要シミュレータ選択処理（ステップ２００１）は、ステップ４０１において生成された境界条件に基づくシミュレーションを実行するシミュレータ１８１０を選択する処理である。具体的には、シミュレータ通信部１８０１は、シミュレータリスト１８０２を参照して、生成された境界条件に適合するシミュレータ１８１０を選択する。
【０１９０】
シミュレーション実行依頼処理（ステップ２００２）は、選択されたシミュレータ１８１０にシミュレーション処理の実行を依頼する処理である。具体的には、シミュレータ通信部１８０１は、ステップ２００１において選択されたシミュレータ１８１０に、ステップ４０１において生成された境界条件を含むシミュレーション処理実行依頼を送信し、そのシミュレータ１８１０からシミュレーション処理結果を受信する。
【０１９１】
シミュレータ１８１０は、ステップ４０１において生成され、シミュレータ通信部１８０１から送信された境界条件に従ってシミュレーション処理を実行する（ステップ４０２）。第１の実施形態のシミュレータ１０２と同様に、本実施形態のシミュレータ１８１０は種々のシミュレーションを実行することができる。したがって、シミュレータ通信部１８０１は、図７に示すようなシミュレーション処理結果を受信してもよいが、本実施形態では別のシミュレーション処理結果（具体的には、シミュレータ１８１０の一つである津波シミュレータによる津波シミュレーション処理の結果）の例について後述する（図２１参照）。
【０１９２】
次に、シナリオ管理部１０３が、ステップ４０２において得られたシミュレーション結果を用いてシナリオＤＢ１０４を構築する（ステップ４０３）。さらに、シナリオ管理部１０３は、シナリオＤＢ１０４に格納されたシナリオを参照して、シナリオクラスタを生成し、その結果をシナリオＤＢ１０４に格納する（ステップ４０４）。これらの処理は、第１の実施形態のシナリオ管理部１０３が実行するものと同様であるため（図８及び図１０参照）、ここでは説明を省略する。したがって、これらの処理によって、シナリオＤＢに図９Ｂ、図９Ｃ及び図１１Ａ〜図１１Ｄに示すようなシナリオが格納されてもよいが、本実施形態では別のシナリオの例について後述する（図２２Ａ〜図２２Ｄ参照）。
【０１９３】
図２１は、本発明の第２の実施形態のシミュレータ１８１０が実行するシミュレーション処理の一例である津波シミュレーションの説明図である。
【０１９４】
図２１に示すグリッド２１０１は、図７に示すグリッド７０１と同様である。図２１の津波シミュレーションの処理結果２１１０に含まれるグリッドＩＤ２１１１、時刻２１１２、水深２１１３及び流速２１１４は、それぞれ、図７の氾濫シミュレーションの処理結果７１０に含まれるグリッドＩＤ７１１、時刻７１２、水深７１３及び流速７１４と同様である。このため、これらについての説明は省略する。ただし、図２１の水深２１１３及び流速２１１４は、河川等の氾濫によるものではなく、津波による浸水深及び水の流れる速度である。
【０１９５】
本実施形態のシナリオ管理部１０３が実行するシナリオＤＢ構築処理（ステップ４０３）及びシナリオクラスタ生成処理（ステップ４０４）の手順は、第１の実施形態のシナリオ管理部１０３によって実行されるものと同様である。シナリオＤＢ構築処理（ステップ４０３）及びシナリオクラスタ生成処理（ステップ４０４）の詳細は、それぞれ図８及び図１０に示した通りである。これによって、本実施形態の物理量推定システム１８００は、図９Ａ〜図９Ｃに示すものと同様のシナリオＤＢ１０４を構築し、図１１Ａ〜図１１Ｄに示すものと同様のシナリオクラスタの生成及び重要スポットの抽出をすることもできる。ただし、本実施形態では、第１の実施形態の道路被災シナリオとは別の例として、図２１の処理結果２１１０に基づいて構築される震災津波シナリオについて、図２２Ａ〜図２２Ｄを参照して説明する。
【０１９６】
図２２Ａは、本発明の第２の実施形態の物理量推定システム１８００によって作成されるシナリオの説明図である。図２２Ａに示す複数のグリッド２２０１は、図２１に示すグリッド２１０１と同等である。
【０１９７】
図２２Ｂは、本発明の第２の実施形態のシナリオＤＢ１０４に含まれるシナリオ管理テーブル２２１０の説明図である。
【０１９８】
シナリオ管理テーブル２２１０は、シナリオＩＤ２２１１及び地震条件データ２２１２を含む。シナリオＩＤ２２１１は、各シナリオの識別情報（ＩＤ）である。地震条件データ２２１２は、シナリオを構築する元になったシミュレーション処理のシミュレート条件（又は境界条件）に相当し、例えば、震源地及び地震の規模などの、地震（及びそれに起因する津波）の初期条件に関するデータである。津波シミュレーション以外のシミュレーションに基づくシナリオの場合は、その元になったシミュレーションの初期条件が地震条件データ２２１２の代わりに登録される。
【０１９９】
図２２Ｃは、本発明の第２の実施形態のシナリオＤＢ１０４に含まれる水深管理テーブル２２２０の説明図である。
【０２００】
水深管理テーブル２２２０は、シナリオＩＤ２２２１、グリッドＩＤ２２２２、時刻２２２３、水深２２２４及び流速２２２５を含む。
【０２０１】
シナリオＩＤ２２２１は、各シナリオの識別情報であり、図２２ＢのシナリオＩＤ２２１１に対応する。
【０２０２】
グリッドＩＤ２２２２は、シミュレーション対象の空間内の各グリッド２２０１の識別情報であり、図２１のグリッドＩＤ２１１１に対応する。
【０２０３】
グリッドＩＤ２２２２によって識別されるグリッドの、時刻２２２３が示す時刻における浸水深及び水の流れる速度（具体的には例えば流速のｘ方向成分及びｙ方向成分）が、それぞれ水深２２２４及び流速２２２５として格納される。
【０２０４】
シナリオ管理部１０３は、図２０のステップ４０３（より詳細には、ステップ４０３に含まれる図８のステップ８０２）において、図２１の処理結果２１１０に基づいて、水深管理テーブル２２２０に登録される情報を算定する。
【０２０５】
図２２Ｃに示す水深管理テーブル２２２０は、一つのシナリオの、一つのグリッドに関する、一つの時刻における水深及び流速を示す。ステップ８０２では、一つのシナリオについて、シミュレーション対象の地域内の各グリッドの各時刻における水深及び流速が算定され、それぞれの値を含む複数の水深管理テーブル２２２０がシナリオＤＢ１０４に格納される。
【０２０６】
図２２Ｄは、本発明の第２の実施形態のシナリオＤＢ１０４に含まれるスマートメータ故障管理テーブル２２３０の説明図である。
【０２０７】
本実施形態のシミュレータ通信部１８０１は、一つ（又は複数）のシミュレータ１８１０から取得したシミュレーション処理結果を、シナリオ管理部１０３だけでなく、シミュレーション条件生成部１０１に渡してもよい。あるシミュレータ１８１０によるシミュレーションの結果に基づいて別のシミュレータ１８１０のシミュレート条件を生成する方法を定義する情報がシミュレーション条件生成部１０１に入力されている場合、シミュレーション条件生成部１０１は、その情報と、シミュレータ通信部１８０１から渡されたシミュレーション処理結果と、に基づいて、新たな境界条件を生成し、それをシミュレータ通信部１８０１に入力する。
【０２０８】
シミュレーションの結果に基づいて別のシミュレート条件を生成する方法を定義する情報の一例は、河川の氾濫又は津波等による各地点の浸水深と、その地点に設置されたセンサ１２０の一種であるスマートメータの故障発生確率との関係である。スマートメータが設置された施設等が浸水した場合、スマートメータへの電力の供給が停止する、スマートメータに浸水して動作が停止する、又はスマートメータが流失する、といった故障のためにそのスマートメータがデータを送信できなくなる確率が上昇する。一般に、浸水深が深くなるほど、スマートメータの故障発生確率は高くなる。
【０２０９】
例えば、シミュレーション条件生成部１０１は、上記のような関係に基づいて生成した境界条件をシミュレータ通信部１８０１に入力する。シミュレータ通信部１８０１は、入力された境界条件を含むシミュレーション実行依頼を、適切に選択されたシミュレータ１８１０（例えば電力シミュレータ）に送信し、シミュレーション処理結果を受信する。このようにして受信したシミュレーション処理結果は、例えば、各地点においてスマートメータの故障が発生する時刻の情報を含んでもよい。
【０２１０】
図２２Ｄには、上記のようなシミュレーション処理結果に基づいて生成されるスマートメータ故障管理テーブル２２３０の例を示す。この例において、スマートメータ故障管理テーブル２２３０は、シナリオＩＤ２２３１、メータＩＤ２２３２及び故障時刻２２３３を含む。
【０２１１】
シナリオＩＤ２２３１は、各シナリオの識別情報であり、図２２ＢのシナリオＩＤ２２１１に対応する。
【０２１２】
メータＩＤ２２３２は、各スマートメータの識別情報である。
【０２１３】
故障時刻２２３３は、各スマートメータが故障すると思われる時刻であり、シミュレーション条件生成部１０１によって生成された境界条件に基づいてシミュレータ１８１０が実行したシミュレーションの処理結果として得られた値である。
【０２１４】
このようにして生成された各シナリオは、少なくとも、津波シミュレータに入力された境界条件と、その境界条件に基づく津波シミュレータによるシミュレーション結果と、そのシミュレーション結果に基づいて生成された境界条件が入力された電力シミュレータによるシミュレーション結果と、の組を含み、さらに、電力シミュレータの境界条件を含んでもよい。
【０２１５】
本実施形態のシナリオ管理部１０３は、図２２Ａ〜図２２Ｄに示すシナリオを対象として、図１０に示すシナリオクラスタ生成処理を実行する。これによって、スマートメータ故障管理テーブル２２３０における故障時刻２２３３を含むシナリオ特徴量が計算され（ステップ１００１）、その特徴量に基づくクラスタリングが行われ（ステップ１００２）、生成されたシナリオクラスタの統計量が算出され（ステップ１００３）、シナリオクラスタの重要スポットが抽出される（ステップ１００４）。
【０２１６】
本実施形態のシナリオ合成処理（ステップ１９０２）の手順は、第１の実施形態のシナリオ合成処理（ステップ３０２）と同様である。したがって、本実施形態の入力受付部１０５及びシナリオ合成部１０６は、図１２に示す実物理量入力処理（ステップ１２０１）及びシナリオ合成処理（ステップ１２０２）を順次実行する。実物理量入力処理及びシナリオ合成処理の詳細は、それぞれ図１３及び図１５に示した通りである。これによって、本実施形態の物理量推定システム１８００は、図１４Ａ、図１４Ｂ及び図１６に例示したものと同じ処理を実行することもできる。ただし、本実施形態では、上記のものとは別の例として、スマートメータから得られるデータも利用してシナリオを合成する処理の例を説明する。このため、本実施形態のステップ１２０１において入力されるデータとして、第１の実施形態と異なる例を図２３Ａ〜図２３Ｃに示す。
【０２１７】
図２３Ａは、本発明の第２の実施形態の物理量推定システム１８００が現況入力システム１３０から受け付けるデータの説明図である。
【０２１８】
物理量推定システム１８００の入力受付部１０５が現況入力システム１３０から受け付ける目視情報のデータ２３００は、例えば、地点２３０１、時刻２３０２及び浸水程度２３０３を含む。この目視情報のデータ２３００は、図１４Ａに示した目視情報のデータ１４００と同様であるため、説明を省略する。
【０２１９】
図２３Ｂは、本発明の第２の実施形態の物理量推定システム１８００がセンサ１２０から受け付けるデータの説明図である。
【０２２０】
図２３Ｂには、物理量推定システム１８００の入力受付部１０５がセンサ１２０から受け付けるデータの一例として、浸水センサ（すなわち液面検知センサ）から送付されるデータ２３１０を示す。データ２３１０は、センサＩＤ２３１１、時刻２３１２及び水深２３１３を含む。このデータ２３１０は、図１４Ｂに示した浸水センサから送付されるデータ１４１０と同様であるため、説明を省略する。
【０２２１】
図２３Ｃは、本発明の第２の実施形態の物理量推定システム１８００が別の種類のセンサ１２０から受け付けるデータの説明図である。
【０２２２】
図２３Ｃには、物理量推定システム１８００の入力受付部１０５がセンサ１２０から受け付けるデータの一例として、スマートメータから送付されるデータ２３２０を示す。データ２３２０は、メータＩＤ２３２１、時刻２３２２及び消費電力２３２３を含む。
【０２２３】
メータＩＤ２３２１は、データ２３２０を送信したセンサ１２０（すなわちスマートメータ）の識別情報である。物理量推定システム１８００は各スマートメータが設置された位置を特定する情報（座標値等）を保持しているため、メータＩＤ２３２１に基づいてセンシングが行われた地点を特定することができる。
【０２２４】
時刻２３２２は、スマートメータによるセンシングが行われた時刻である。
【０２２５】
消費電力２３２３は、スマートメータによってセンシングされた（すなわちスマートメータが設置された施設等の）消費電力量（例えばワット／時間）である。河川の氾濫又は津波等によって、スマートメータが設置された施設等が浸水した場合、スマートメータの故障のためにデータ２３２０を送信できなくなる確率が上昇する。物理量推定システム１８００は、例えば、あるスマートメータからデータ２３２０を受信できない時間が所定の長さを超えた場合、そのスマートメータが故障したと判定してもよい。
【０２２６】
本実施形態のシナリオ合成部１０６は、データ２３２０に基づくスマートメータの故障の判定結果を含む現況シナリオ特徴量を計算し（ステップ１５０１）、その現況シナリオ特徴量と、ステップ１００２で生成されたシナリオクラスタの特徴量とに基づいて、現況シナリオ特徴量に最も類似するシナリオクラスタを選択する。さらに、シナリオ合成部１０６は、選択したシナリオクラスタの重要スポットについて、第１の実施形態と同様の方法によって重要スポット情報推定処理を実行し（ステップ１５０３）、その結果と、ステップ１５０２において選択されたシナリオクラスタの代表値とを含む合成シナリオを生成し（ステップ１５０４）、生成した合成シナリオを合成シナリオＤＢ１０７に格納する。
【０２２７】
図２４は、本発明の第２の実施形態のアプリケーションシステム１４０が実行する被害地域情報生成処理のフローチャートである。
【０２２８】
図２４には、アプリケーションシステム１４０が提供するアプリケーションの一例として、津波による被害を受けた地域（又はこれから被害を受けると予想される地域）を表示するアプリケーションを示す。
【０２２９】
具体的には、アプリケーションシステム１４０は、水深可視化処理を実行する（ステップ２４０１）。この処理は、アプリケーションシステム１４０が、シナリオ配信部１０８から配信された合成シナリオの、各地点（グリッド）の、各時刻の浸水深を参照し、その浸水深をアプリケーションシステム１４０の出力装置２０４（例えば画像表示装置）から出力する処理である。例えば、浸水深は、数値によって表示されてもよいし、地図上の記号、図形又は色彩等によって表示されてもよい。
【０２３０】
以上で被害地域情報生成処理が終了する。なお、上記のような被害地域情報生成処理は一例であり、アプリケーションシステム１４０は、合成シナリオを利用して種々の処理（例えば第１の実施形態と同様の誘導情報生成処理）を実行することができる。
【０２３１】
以上の本発明の第２の実施形態によれば、物理量推定システム１８００がシミュレータ１０２を有しない場合であっても、外部のシミュレータ１８１０を利用して、第１の実施形態と同様の物理量推定を行うことができる。
【０２３２】
また、あるシミュレータによるシミュレーションの結果と、別のシミュレータによるシミュレーションの入力条件とを関連付ける情報をあらかじめ設定しておくことによって、本実施形態の物理量推定システム１８００は、センサ１２０等（例えばスマートメータ）から取得されたリアルタイムデータを、本来の目的とは異なる物理量の推定（例えば浸水深の推定）に利用することができる。
【０２３３】
より詳細には、ある地域においてある事象が発生した場合に、その事象が当該地域内の複数の種類のセンサ１２０が取得するデータに影響を与える場合がある。例えば、ある地域において津波が発生し、それが原因で停電した場合、水位計（潮位計）は高い水位を計測し、スマートメータによって計測される電力量はゼロになる（又はスマートメータが故障することによってデータの送信が途絶える）。このように、複数種類のセンサが一つの事象に関連するデータを取得する場合がある。しかし、各種のセンサ１２０は、例えば、津波の検知、電力供給の制御又は課金、等の、それぞれに固有の目的のために使用されるものであり、それぞれのセンサ１２０が取得したデータを、それぞれのセンサ１２０の本来の目的を超えて、それらが設置された地域に発生した事象に関連する物理量を推定するために利用することはできなかった。
【０２３４】
例えば、上記のように、水位計と、水位計が計測した値を条件として入力されると津波シミュレーションを実行する津波シミュレータと、スマートメータと、スマートメータが計測した値が条件として入力されると電力に関するシミュレーションを実行する電力シミュレータと、が用意されている場合に、スマートメータの計測結果を、津波シミュレーションによる浸水深の推定に利用することはできなかった。しかし、本実施形態によれば、上記のように、センサ１２０等から取得されたリアルタイムデータを、本来の目的とは異なる物理量の推定に利用することが可能になる。
【０２３５】
これによって、十分な数のセンサ１２０が設置されていない場合であっても、物理量の推定の精度を向上させることができる。例えば、水位計は、河川及び港湾等には設置されていても、人の居住地域に設置されていることは少ない。一方、スマートメータ（電力計）は、人の居住地域に豊富に存在することが期待できる。スマートメータの計測値を津波による浸水深の推定に利用することによって、水位計の乏しい居住地域における津波による浸水深の推定の精度を向上させることができる。
【０２３６】
なお、本発明の第１の実施形態と第２の実施形態の本質的な相違点は、シミュレータが物理量推定システムの内部にあるか否かだけである。このため、上記のようにセンサを本来の目的以外の目的に使用することは、第１の実施形態の物理量推定システム１００によっても実現することができる。
【０２３７】
上記の第１及び第２の実施形態ではシナリオのクラスタリング及び重要スポットの特定が行われるが、これらの実施形態の変形例として、クラスタリング等を省略する形態も考えられる。例えば、物理量を推定しようとする対象地域が比較的狭く、発生が予想される事象も限定されているような場合には、物理量推定システム１００（又は１８００）は、クラスタリングを行わずに（言い換えるとすべてのシナリオを一つのシナリオクラスタに属するものとみなして）物理量の推定を行うこともできる。このとき、物理量推定システム１００は、重要スポットの特定を実行してもよいが、これを実行せずに、すべての地点及び時刻について、類似度に基づく重み付け演算によるシナリオの合成を行ってもよい。
【０２３８】
上記の実施形態の各機能を実現するプログラム、テーブル、ファイル等の情報は、不揮発性半導体メモリ、ハードディスクドライブ、ＳＳＤ（Ｓｏｌｉｄ  Ｓｔａｔｅ  Ｄｒｉｖｅ）等の記憶デバイス、または、ＩＣカード、ＳＤカード、ＤＶＤ等の計算機読み取り可能な非一時的データ記憶媒体に格納することができる。
【０２３９】
本発明は上記した実施形態に限定されるものではなく、様々な変形例が含まれる。例えば、上記した実施形態は本発明を分かりやすく説明するために詳細に説明したものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。また、ある実施形態の構成の一部を他の実施形態の構成に置き換えることが可能であり、また、ある実施形態の構成に他の実施形態の構成を加えることも可能である。また、各実施形態の構成の一部について、他の構成の追加・削除・置換をすることが可能である。
【図１】
2014020747.tif 000003
【図２】
2014020747.tif 000004
【図３Ａ】
2014020747.tif 000005
【図３Ｂ】
2014020747.tif 000006
【図４】
2014020747.tif 000007
【図５】
2014020747.tif 000008
【図６】
2014020747.tif 000009
【図７】
2014020747.tif 000010
【図８】
2014020747.tif 000011
【図９Ａ】
2014020747.tif 000012
【図９Ｂ】
2014020747.tif 000013
【図９Ｃ】
2014020747.tif 000014
【図１０】
2014020747.tif 000015
【図１１Ａ】
2014020747.tif 000016
【図１１Ｂ】
2014020747.tif 000017
【図１１Ｃ】
2014020747.tif 000018
【図１１Ｄ】
2014020747.tif 000019
【図１２】
2014020747.tif 000020
【図１３】
2014020747.tif 000021
【図１４Ａ】
2014020747.tif 000022
【図１４Ｂ】
2014020747.tif 000023
【図１５】
2014020747.tif 000024
【図１６】
2014020747.tif 000025
【図１７】
2014020747.tif 000026
【図１８】
2014020747.tif 000027
【図１９Ａ】
2014020747.tif 000028
【図１９Ｂ】
2014020747.tif 000029
【図２０】
2014020747.tif 000030
【図２１】
2014020747.tif 000031
【図２２Ａ】
2014020747.tif 000032
【図２２Ｂ】
2014020747.tif 000033
【図２２Ｃ】
2014020747.tif 000034
【図２２Ｄ】
2014020747.tif 000035
【図２３Ａ】
2014020747.tif 000036
【図２３Ｂ】
2014020747.tif 000037
【図２３Ｃ】
2014020747.tif 000038
【図２４】
2014020747.tif 000039
【手続補正書】
【提出日】平成27年11月12日(2015.11.12)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
シミュレーションの条件として入力された値と、前記入力された値に基づく所定のシミュレーションの結果として算出された、一つ以上の位置における物理量と、を含むシナリオを複数保持するシナリオ保持部と、
計測された物理量及び前記物理量が計測された位置に基づく所定の計算によって前記複数のシナリオに含まれる物理量を合成することによって、前記一つ以上の位置における物理量を推定するシナリオ合成部と、
前記推定された物理量を保持する合成シナリオ保持部と、を有し、
前記シナリオ合成部は、前記計測された物理量及び前記物理量が計測された位置に基づく所定の計算として、前記計測された物理量と、前記複数のシナリオに含まれる、前記物理量が計測された位置における物理量との類似度に従って定めた重みを、前記複数のシナリオに含まれる物理量に付けた重み付け演算を実行することを特徴とする物理量推定システム。
【請求項２】
請求項１に記載の物理量推定システムであって、
前記複数のシナリオを、前記各シナリオに含まれる物理量に基づくクラスタリングによって、各々が複数のシナリオを含む複数のシナリオクラスタに分類し、前記各シナリオクラスタに含まれる複数のシナリオの前記各位置における物理量の統計量を算出し、前記各シナリオクラスタに含まれる複数のシナリオの前記各位置における物理量のばらつきが所定の値より大きい場合、前記各位置を重要スポットとして特定し、前記各シナリオが属する前記シナリオクラスタを特定する情報、及び、前記各シナリオクラスタの前記重要スポットを特定する情報を前記シナリオ保持部に格納するシナリオ管理部をさらに有し、
前記シナリオ合成部は、
前記各シナリオクラスタについて算出された前記統計量と、前記計測された物理量とを比較することによって、前記計測された物理量と最も類似するシナリオクラスタを特定し、
前記特定されたシナリオクラスタの前記重要スポットにおける物理量を、前記重み付け演算によって推定し、
前記特定されたシナリオクラスタの前記重要スポット以外の位置における物理量を、前記重要スポットにおいて使用されるものとは異なる重みを付けた重み付け演算によって推定することを特徴とする物理量推定システム。
【請求項３】
請求項２に記載の物理量推定システムであって、
第１の種類の前記シミュレーションの結果に含まれる前記一つ以上の位置における第１の種類の前記物理量に基づいて、第２の種類の前記シミュレーションの条件として入力される値を生成するシミュレーション条件生成部をさらに有し、
前記各シナリオは、少なくとも、前記第１の種類のシミュレーションの条件、前記第１の種類のシミュレーションの結果として得られる前記第１の種類の物理量、及び、前記第２の種類のシミュレーションの結果として得られる第２の種類の物理量を含み、
前記シナリオ合成部は、前記各シナリオクラスタについて算出された前記統計量と、第１の種類のセンサによって計測された前記第１の種類の物理量と、第２の種類のセンサによって計測された前記第２の種類の物理量と、を比較することによって、前記計測された物理量と最も類似するシナリオクラスタを特定することを特徴とする物理量推定システム。
【請求項４】
請求項３に記載の物理量推定システムであって、
入力された前記第１の種類の物理量に基づいて、一つ以上の位置における前記第１の種類の物理量を算出する第１の種類のシミュレータ部と、入力された前記第２の種類の物理量に基づいて、一つ以上の位置における前記第２の種類の物理量を算出する第２の種類のシミュレータ部と、をさらに有することを特徴とする物理量推定システム。
【請求項５】
請求項２に記載の物理量推定システムであって、
一つの前記シミュレーション条件を入力されると、各々が、前記入力されたシミュレーション条件に含まれる値を含む範囲内の複数の値の各々を含む、複数のシミュレーション条件を生成するシミュレーション条件生成部をさらに有することを特徴とする物理量推定システム。
【請求項６】
請求項２に記載の物理量推定システムであって、
前記各シナリオクラスタについて算出された前記統計量は、前記各シナリオクラスタ内の前記複数のシナリオの各位置における前記物理量に同一の重みを付けて計算された平均値を含み、
前記シナリオ合成部は、
前記特定されたシナリオクラスタに含まれる前記複数のシナリオの前記重要スポットにおける物理量に、前記計測された物理量との類似度が高いほど大きい重みを付けて計算した平均値を、前記特定されたシナリオクラスタの前記重要スポットにおける物理量として推定し、
前記特定されたシナリオクラスタの前記重要スポット以外の位置について算出された前記統計量に含まれる平均値を、前記特定されたシナリオクラスタの前記重要スポット以外の位置における物理量として推定することを特徴とする物理量推定システム。
【請求項７】
請求項２に記載の物理量推定システムであって、
前記シナリオ合成部は、
前記特定されたシナリオクラスタに含まれる前記複数のシナリオの前記重要スポットにおける物理量に、前記計測された物理量との類似度が高いほど大きい重みを付けて計算した平均値を、前記特定されたシナリオクラスタの前記重要スポットにおける物理量として推定し、
前記特定されたシナリオクラスタに含まれる前記複数のシナリオの前記重要スポット以外の位置における物理量のうち、前記計測された物理量との類似度が最も高いものの重みとして０より大きい値を付け、それ以外の物理量の重みとして０を付けて計算した平均値を、前記特定されたシナリオクラスタの前記重要スポット以外の位置における物理量として推定することを特徴とする物理量推定システム。
【請求項８】
物理量推定システムが実行する物理量推定方法であって、
前記物理量推定システムは、プロセッサと、前記プロセッサに接続されるインターフェースと、前記プロセッサに接続される記憶装置と、を有し、
前記記憶装置は、シミュレーションの条件として入力された値と、前記入力された値に基づく所定のシミュレーションの結果として算出された、一つ以上の位置における物理量と、を含むシナリオを複数保持し、
前記物理量推定方法は、
前記インターフェースが、計測された物理量の入力を受ける第１手順と、
前記プロセッサが、前記計測された物理量及び前記物理量が計測された位置に基づく所定の計算によって前記複数のシナリオに含まれる物理量を合成することによって、前記一つ以上の位置における物理量を推定し、前記推定された物理量を前記記憶装置に格納する第２手順と、を含み、
前記第２手順は、前記プロセッサが、前記計測された物理量と、前記複数のシナリオに含まれる、前記物理量が計測された位置における物理量との類似度に従って定めた重みを、前記複数のシナリオに含まれる物理量に付けた重み付け演算を実行する手順を含むことを特徴とする物理量推定方法。
【請求項９】
請求項８に記載の物理量推定方法であって、
前記プロセッサが、前記複数のシナリオを、前記各シナリオに含まれる物理量に基づくクラスタリングによって、各々が複数のシナリオを含む複数のシナリオクラスタに分類し、前記各シナリオクラスタに含まれる複数のシナリオの前記各位置における物理量の統計量を算出し、前記各シナリオクラスタに含まれる複数のシナリオの前記各位置における物理量のばらつきが所定の値より大きい場合、前記各位置を重要スポットとして特定し、前記各シナリオが属する前記シナリオクラスタを特定する情報、及び、前記各シナリオクラスタの前記重要スポットを特定する情報を前記記憶装置に格納する第３手順をさらに含み、
前記第２手順は、
前記プロセッサが、前記各シナリオクラスタについて算出された前記統計量と、前記計測された物理量とを比較することによって、前記計測された物理量と最も類似するシナリオクラスタを特定する第４手順と、
前記プロセッサが、前記特定されたシナリオクラスタの前記重要スポットにおける物理量を、前記重み付け演算によって推定する第５手順と、
前記プロセッサが、前記特定されたシナリオクラスタの前記重要スポット以外の位置における物理量を、前記重要スポットにおいて使用されるものとは異なる重みを付けた重み付け演算によって推定する第６手順と、を含むことを特徴とする物理量推定方法。
【請求項１０】
請求項９に記載の物理量推定方法であって、
前記プロセッサが、第１の種類の前記シミュレーションの結果に含まれる前記一つ以上の位置における第１の種類の前記物理量に基づいて、第２の種類の前記シミュレーションの条件として入力される値を生成する第７手順をさらに含み、
前記各シナリオは、少なくとも、前記第１の種類のシミュレーションの条件、前記第１の種類のシミュレーションの結果として得られる前記第１の種類の物理量、及び、前記第２の種類のシミュレーションの結果として得られる第２の種類の物理量を含み、
前記第４手順は、前記各シナリオクラスタについて算出された前記統計量と、第１の種類のセンサによって計測された前記第１の種類の物理量と、第２の種類のセンサによって計測された前記第２の種類の物理量と、を比較する手順を含むことを特徴とする物理量推定方法。
【請求項１１】
請求項１０に記載の物理量推定方法であって、
前記プロセッサが、入力された前記第１の種類の物理量に基づいて、一つ以上の位置における前記第１の種類の物理量を算出し、入力された前記第２の種類の物理量に基づいて、一つ以上の位置における前記第２の種類の物理量を算出する第８手順をさらに含むことを特徴とする物理量推定方法。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】０００６
【補正方法】変更
【補正の内容】
【０００６】
上記の課題を解決するために、本発明は、シミュレーションの条件として入力された値と、前記入力された値に基づく所定のシミュレーションの結果として算出された、一つ以上の位置における物理量と、を含むシナリオを複数保持するシナリオ保持部と、計測された物理量及び前記物理量が計測された位置に基づく所定の計算によって前記複数のシナリオに含まれる物理量を合成することによって、前記一つ以上の位置における物理量を推定するシナリオ合成部と、前記推定された物理量を保持する合成シナリオ保持部と、を有し、前記シナリオ合成部は、前記計測された物理量及び前記物理量が計測された位置に基づく所定の計算として、前記計測された物理量と、前記複数のシナリオに含まれる、前記物理量が計測された位置における物理量との類似度に従って定めた重みを、前記複数のシナリオに含まれる物理量に付けた重み付け演算を実行する。
【国際調査報告】
2014020747.tif 000040
2014020747.tif 000041