(57)【要約】【課題】災害現場からの観測データに基づいて模擬規則を決定し、被災範囲予測を行う災害監視装置を得る。【解決手段】監視部２は、統合指示部１からの命令に従って災害現場から観測データを取得し、現状データ生成部３は観測データを基に災害状況を簡略化して表す現状データを生成する。模擬規則推定部４は、統合指示部１からの命令に従って現状データに基づいて推定規則を生成する。模擬部５は、統合指示部１からの命令に従って現状データに基づき、生成された推定規則を用いて被災範囲予測を行う。【選択図】        図１
2004295633.tif 000002
【特許請求の範囲】
【請求項１】
指定された地域、時刻および事象に基づいて、災害現場からの観測情報を取得する命令と、上記観測情報に基づく被災範囲の予測に用いる推定規則を生成する命令と、上記推定規則を用いて被災範囲の予測を行う命令を作成する統合指示部と、
上記統合指示部が作成した観測情報取得の命令に基づいて災害現場から観測機器を介して観測データを取得する監視部と、
上記観測データから災害状況を簡略化して表す現状データを生成する現状データ生成部と、
災害現場で発生する事象の情報を発生地点、発生時刻と共に格納する事象データベースから、与えられた条件に合致する事象の情報を抽出する事象管理部と、
上記統合指示部が作成した推定規則生成の命令に従い、上記現状データに上記事象管理部から取得した事象の情報を反映して推定規則を生成する模擬規則推定部と、
上記統合指示部が作成した被災範囲予測の命令に従い、上記現状データに上記事象管理部から取得した事象の情報を反映し、上記推定規則を用いて被災範囲予測を行う模擬部と、
上記模擬部による被災範囲予測の結果を上記現状データ等と共に表示する統合表示部とを備えた災害監視装置。
【請求項２】
模擬規則推定部により推定される模擬規則は、模擬対象地域を複数のセルの集合として表し、上記模擬規則に従ってセル間で状態変化を伝播させていくことにより災害の模擬データを作成する規則であり、
上記模擬規則推定部は、上記セル間での状態変化の伝播の発生の尤度が最大となるように上記模擬規則を生成することを特徴とする請求項１記載の災害監視装置。
【請求項３】
模擬規則推定部により推定される模擬規則は、推定対象変数として個々のセルの隣接セル範囲、各種類のセル間での伝播確率、各種類のセルの特性値およびセル間での各事象の伝播が阻害される確率を含み、
上記模擬規則推定部は、上記個々のセルの隣接セル範囲を、セル間での状態変化の伝播の発生の尤度が最大となるように算定する隣接範囲推定部と、
上記各種類のセル間での伝播確率を、セル間での状態変化の伝播の発生の尤度が最大となるように算定する伝播確率推定部と、
上記各種類のセルの特性値を、セル間での状態変化の伝播の発生の尤度が最大となるように算定するセル特性値推定部と、
上記セル間での各事象の伝播が阻害される確率を、セル間での状態変化の伝播の発生の尤度が最大となるように算定する伝播阻害確率推定部とを備えたことを特徴とする請求項２記載の災害監視装置。
【請求項４】
模擬規則推定部は、模擬対象地域を所与の数の模擬範囲に分割し、各セルが所属する模擬範囲を推定対象変数としてセル間での状態変化の伝播の発生の尤度が最大となるように算定する均一領域推定部を備えたことを特徴とする請求項２または請求項３記載の災害監視装置。
【請求項５】
模擬規則推定部は、模擬開始から模擬終了までの時間を所与の数の模擬時間に分割し、各セルが所属する模擬時間を推定対象変数としてセル間での状態変化の伝播の発生の尤度が最大となるように算定する均一時間帯推定部を備えたことを特徴とする請求項２または請求項３記載の災害監視装置。
【請求項６】
模擬規則推定部は、伝播の方向と大きさを推定対象変数としてセル間での状態変化の伝播の発生の尤度が最大となるように算定する風向風速推定部を備えたことを特徴とする請求項２または請求項３記載の災害監視装置。
【請求項７】
模擬規則推定部は、セルの種類の指定に誤りがあるときは、種類を修正するセルを、推定対象変数としてセル間での状態変化の伝播の発生の尤度が最大となるように決定するセル種類修正部を備えたことを特徴とする請求項２または請求項３記載の災害監視装置。
【請求項８】
模擬規則推定部は、伝播の時間単位を、推定対象変数としてセル間での状態変化の伝播の発生の尤度が最大となるように定める時間単位推定部を備えたことを特徴とする請求項２または請求項３記載の災害監視装置。
【請求項９】
模擬規則推定部は、与えられた事象について、該事象が模擬データ上で欠落していると仮定し、該事象の発生地点と発生時刻を推定対象変数としてセル間での状態変化の伝播の発生の尤度が最大となるように推定する事象推定部を備えたことを特徴とする請求項２または請求項３記載の災害監視装置。
【請求項１０】
模擬規則推定部は、模擬規則に含まれる推定対象変数の算定を行う共通推定部を備え、
上記共通推定部は、上記推定対象変数の初期値および伝播尤度の暫定最適値を設定する変数初期値決定部と、
通知された推定対象変数値に基づいて伝播尤度を算出する伝播尤度算出部と、
上記伝播尤度算出部が算出した伝播尤度と上記暫定最適値を比較し、該伝播尤度が上記暫定最適値よりも大きい場合には、対応する推定対象変数値と伝播尤度を新たな暫定最適値とすると共に変数変更部に推定対象変数の変更命令を通知し、一方、伝播尤度の暫定最適値が所与の回数の推定対象変数の変更を行っても変化しなくなった場合には、上記共通推定部の処理を終了する変数推定終了判定部と、
上記変数推定終了判定部からの推定対象変数の変更命令を受け、上記推定対象変数を所与の方式で変更し、上記伝播尤度算出部に供給する変数変更部を有することを特徴とする請求項２から請求項９のうちのいずれか１項記載の災害監視装置。
【請求項１１】
伝播尤度算出部は、各セルの隣接セルを抽出する隣接網作成部と、
状態が異なるセル同士の境界に位置するセルを外縁セルとして抽出する外縁抽出部と、
上記外縁セル間の伝播パスを列挙する外縁間パス算出部と、
上記外縁間パス算出部が列挙した各伝播パスの伝播尤度を算出し、それらの総和を求める伝播尤度算出部とを備えたことを特徴とする請求項１０記載の災害監視装置。
【請求項１２】
統合指示部は、指定された地域の未来の時刻における計画データ、該計画データに反映する事象および使用する推定規則に基づいて、災害現場からの観測情報を取得する命令、上記計画データを実現するための計画策定を行う命令を作成し、
監視部は、上記統合指示部が作成した観測情報取得の命令に基づいて災害現場から観測データを取得し、
現状データ生成部は、上記現状画像を基に災害状況を簡略化して表す現状データを生成し、
模擬規則推定部は、上記統合指示部が作成した計画策定を行う命令、上記現状データ、上記反映する事象および指定された推定規則に基づいて、事象推定部により該事象の発生地点と発生時刻を推定して事象計画を生成し、
統合表示部は、上記事象計画を表示し、事象管理部は、上記事象計画を事象データベースに登録することを特徴とする請求項９記載の災害監視装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
この発明は、災害による被害を予測し、対策決定を支援する災害監視装置に関するものである。
【０００２】
【従来の技術】
火災や洪水などの災害による被害の範囲を最小限に抑えるためには、災害現場の被災状況に即して適切な対策を迅速に実行することが重要である。このため、災害現場に設置された観測機器から得られる監視画像や観測データ等に基づいて、セルラーオートマトン等の原理を利用した模擬規則により模擬データを作成し、被災範囲の予測を行うシステムが多数提案されている。得られた予測結果を基に、予想される被害をできるだけ防ぐように対策を立てることができる。
【０００３】
例えば、特許文献１に開示された従来の環境シミュレータでは、セルラーオートマトンの原理に基づく所与の模擬規則によって、拡散現象を擬似的に解析することにより、被災状況の模擬データを生成している。
また、特許文献２に開示された従来の３次元火災延焼シミュレーションシステムでは、同じくセルラーオートマトンの原理に基づいて、所与のセル材質および延焼規則に従って火災延焼の模擬データを生成している。
【０００４】
【特許文献１】
特開２００１−１３３３５９号公報
【特許文献２】
特開平８−２４９３１３号公報
【０００５】
【発明が解決しようとする課題】
従来の災害監視装置においては、利用者が予め模擬規則を決定し、それに従って被災範囲の予測を行っていた。そのため災害現場から得られた監視画像や観測データに基づいて適切な模擬規則を決定しようとするとそのための時間を要し、実用的な時間内に被災範囲予測データを得られないという問題があった。
【０００６】
この発明は上記のような課題を解決するためになされたもので、災害現場から得られる監視データに基づいて自動的に模擬規則を決定し、被災範囲予測を行う災害監視装置を得ることを目的とする。
【０００７】
【課題を解決するための手段】
この発明に係る災害監視装置は、指定された地域、時刻および事象に基づいて、災害現場からの観測情報を取得する命令と、観測情報に基づく被災範囲の予測に用いる推定規則を生成する命令と、推定規則を用いて被災範囲の予測を行う命令を作成する統合指示部と、統合指示部が作成した観測情報取得の命令に基づいて災害現場から観測機器を介して観測データを取得する監視部と、観測データから災害状況を簡略化して表す現状データを生成する現状データ生成部と、災害現場で発生する事象の情報を発生地点、発生時刻と共に格納する事象データベースから、与えられた条件に合致する事象の情報を抽出する事象管理部と、統合指示部が作成した推定規則生成の命令に従い、現状データに事象管理部から取得した事象の情報を反映して推定規則を生成する模擬規則推定部と、統合指示部が作成した被災範囲予測の命令に従い、現状データに事象管理部から取得した事象の情報を反映し、推定規則を用いて被災範囲予測を行う模擬部と、模擬部による被災範囲予測の結果を現状データ等と共に表示する統合表示部とを備えたものである。
【０００８】
【発明の実施の形態】
以下、この発明の実施の形態を説明する。
実施の形態１．
図１は、この発明の実施の形態１による、災害監視装置１００の構成を示すブロック図である。図に示すように、災害監視装置１００は、統合指示部１、監視部２、現状データ生成部３、模擬規則推定部４、模擬部５、事象管理部６、および統合表示部７を備える。統合指示部１、監視部２、現状データ生成部３、模擬規則推定部４、模擬部５、事象管理部６および統合表示部７は、実際には災害監視装置１００のプロセッサの動作を制御するプログラムのモジュールに従って便宜的に分割したものである。
また、監視データベース８および事象データベース９は、リレーショナル型データベース等によって管理されるデータレコードの集合であり、記憶装置に格納される。監視データベース８および事象データベース９は、災害監視装置１００と通信回線を介して接続される。なお、監視データベース８および事象データベース９は、同一の記憶装置に格納されていてもよい。
監視部２は、例えば光学カメラ、赤外線カメラ、合成開口レーダ等の複数の観測機器と接続されており、それらの観測機器から災害現場の複数の監視画像や観測データを取得することができる。取得したデータは監視データベース８に格納される。
【０００９】
図２は模擬規則推定部４の詳細を示すブロック図である。図に示すように、模擬規則推定部４は、隣接範囲推定部４１、伝播確率推定部４２、セル特性値推定部４３、伝播阻害確率推定部４４、均一領域推定部４５１、均一時間帯推定部４５２、風向風速推定部４５３、セル種類修正部４５４、時間単位修正部４５５、事象推定部４５６、および共通推定部４６、終了判定部４７、模擬規則格納部４８を備える。隣接範囲推定部４１、伝播確率推定部４２、セル特性値推定部４３、伝播阻害確率推定部４４、均一領域推定部４５１、均一時間帯推定部４５２、風向風速推定部４５３、セル種類修正部４５４、時間単位修正部４５５、事象推定部４５６、共通推定部４６、終了判定部４７および模擬規則格納部４８は、実際には災害監視装置１００のプロセッサの動作を制御するプログラムのモジュールに従って便宜的に分割したものである。
【００１０】
また、図３は共通推定部４６の詳細を示すブロック図である。図に示すように、共通推定部４６は、変数初期値決定部４６１、伝播尤度算出部４６２、終了判定部（変数推定終了判定部）４６３、変数変更部４６４を備える。また、伝播尤度算出部４６２は、隣接網作成部４６２１、外縁抽出部４６２２、外縁間伝播尤度算出部（外縁間パス算出部、伝播尤度算出部）４６２３を備える。
【００１１】
次に、実施の形態１による災害監視装置１００の動作について説明する。
統合指示部１は、システム利用者が入力装置（図示せず）を介して入力した内容に従って、災害監視装置１００の各部に与える複数の命令を作成し、各部へ通知する。統合指示部１が作成する命令には、例えば、監視部２に災害現場の監視データの取得を命令する監視指示、模擬規則推定部４に被災範囲予測のための模擬規則を推定させる規則推定指示、模擬規則推定部４に災害対策のための計画を作成させる計画策定指示、模擬部５に模擬規則を用いて被災範囲予測を行わせる模擬指示がある。
【００１２】
図５に、システム利用者が入力操作を行うときに使用するディスプレイ（図示せず）に表示される入力画面１６１の例を示す。入力画面１６１は、例えばタッチパネル等のＧＵＩ（Ｇｒａｐｈｉｃａｌ  Ｕｓｅｒ  Ｉｎｔｅｒｆａｃｅ）を介して項目が指定できるようになっており、「表示」ボタン１６２、「模擬」ボタン１６３、「規則推定」ボタン１６４、「計画策定」ボタン１６５を選択することにより、作成する命令を指定することができる。システム利用者は、作成する命令を指定したら、地域指定領域１６６、時刻指定領域１６７、事象指定領域１６８に命令に含まれる具体的な内容を入力する。
【００１３】
システム利用者は「表示」ボタン１６２を指定することにより、監視部２に所望の監視画像あるいは観測データを取得させ、それらをディスプレイに表示させる命令を作成することができる。この場合、地域指定領域１６６で監視画像や観測データの表示対象地域を指定することができる。地域の指定方法は、例えば入力装置を使って緯度経度を指定してもよいし、画面に地図を表示させて指定地域を選択するようにしてもよい。その他、地域を正確に指定できる方法であればどのようなものであってもよい。時刻指定領域１６７では表示対象時刻を時刻列として指定することができる。表示対象時刻は、監視画像や観測データが取得された時刻であり、時刻列は複数の表示対象時刻から構成される。具体的には、例えば開始時刻、表示終了時刻、および表示間隔を入力することにより、複数の表示対象時刻が指定される。事象指定領域１６８では例えば観測機器の選択等を行うことができる。これにより、システム利用者は、指定した地域における災害の監視画像や観測データを一定の時間間隔で連続して取得することができる。
「模擬」ボタン１６３を選択した場合には、後述する被災範囲予測処理を行うための各種命令を作成することができる。この場合には、地域指定領域１６６で被災範囲予測を行う対象地域を指定し、時刻指定領域１６７で例えば被災範囲予測の開始時刻と終了時刻、および時間間隔を入力し、予測結果の表示を行う時刻の列（以下、模擬時刻列という。）を指定することができる。また、事象指定領域１６８で、模擬に反映させる事象を指定することができる。ここでの事象とは、例えば火災による被災範囲を予測するのであれば、出火、鎮火等が該当する。
「規則推定」ボタン１６４を選択した場合には、後述する模擬規則推定処理を行うための各種命令を作成することができる。この場合には、地域指定領域１６６で模擬規則推定の対象地域を指定する。時刻指定領域１６７で模擬規則推定に用いる監視画像や観測データを指定する時刻の列（以下、規則推定時刻列という。）を指定し、事象指定領域１６８で規則推定に反映させる事象を指定することができる。
「計画策定」ボタン１６５を選択した場合には、後述する計画策定処理を行うための各種命令を作成することができる。この場合には、地域指定領域１６６で計画データの作成を行い、時刻指定領域１６７で計画データの時刻指定を行い、事象指定領域１６８で計画したい事象の指定を行うことができる。
【００１４】
次に、図４のフローチャートを用いて、この発明の実施の形態１による模擬規則の自動推定処理、および推定された模擬規則に基づく被災範囲予測処理について説明する。
【００１５】
統合指示部１は、システム利用者が「規則推定」ボタン１６４を選択して指定した内容に基づいて模擬規則推定部４に通知する規則推定指示を作成する。また、システム利用者が「模擬」ボタン１６３を選択して指定した内容に基づいて模擬部５に通知する模擬指示を作成する。また、統合指示部１は監視部２に通知する監視指示を作成する（ステップＳＴ１０１）。
【００１６】
監視指示には、システム利用者が画面上で指定した模擬対象地域、規則推定時刻列、および模擬時刻列中の模擬開始時刻が含まれる。
規則推定指示には、模擬対象地域および規則推定時刻列が含まれる。また、模擬規則推定部４には、規則推定指示と共に事象獲得指示が通知される。事象獲得指示には、システム利用者が画面上で指定した反映事象が含まれる。
模擬指示には、システム利用者が画面上で指定した模擬対象地域および模擬時刻列が含まれる。模擬部５には、模擬指示と共に事象獲得指示が通知される。事象獲得指示には、システム利用者が画面上で指定した反映事象が含まれる。
【００１７】
統合指示部１は、作成した監視指示を監視部２に通知し、規則推定指示および事象獲得指示を模擬規則推定部４に通知し、模擬指示および事象獲得指示を模擬部５に通知する。
【００１８】
さらに、統合指示部１は、システム利用者が入力画面上で指定した反映事象を事象管理部６に通知し、事象管理部６は、該事象を事象データベース９にユーザ定義事象として登録する。
【００１９】
ここで、図６を用いて事象データベース９に格納される情報の内容について説明する。事象データベース９は、事象内容、事象の発生した時刻および事象の発生した位置の緯度経度によって定義される事象のレコードの集合である。事象内容には、例えば、出火、鎮火のようにある位置における災害の発生状況を表す事象と、現場からの被災・救出人数報告、消火・救援計画に従った消火要員・消火機材・救援要員の配置のように災害対策の目的で行われる事象がある。それぞれの事象は、事象定義によって分類される。事象定義は「規定」と「ユーザ定義」の２通りがある。「規定」とは、対応する事象内容が既に発生した事実であるか、あるいは発生することが計画されている事象であることを示す。一方、「ユーザ定義」とは、システム利用者によって模擬規則推定などで用いるために仮に設定された事象であることを表す。例えば、システム利用者が入力画面１６１の事象指定領域１６８で指定する規則推定に反映させる事象の情報は、ユーザ定義事象として事象データベース９に格納される。図６に示す例では、２種類の消火事象のうち消防車の消火は、規定事象として当該時刻に発生することが計画されており、ヘリコプターの消火は、ユーザが仮に設定した事象であることを示す。
【００２０】
監視部２は、統合指示部１から通知された監視指示に基づいて監視データベース８から観測情報を取得し現状画像を生成する。あるいは、観測機器を制御して監視画像や観測データを取得し、現状画像を生成する（ステップＳＴ１０２）。
【００２１】
ここで、図７を用いて監視データベース８に格納される情報の内容について説明する。監視データベース８には、観測情報が、観測時刻、観測範囲、観測手段と関連付けて記憶されている。観測情報は、光学カメラ、赤外線カメラ、合成開口レーダなどの観測機器によって得られた監視画像や観測データであり、監視データベース８には、それらのデータファイルへのリンク情報が保持される。観測時刻は観測機器によって観測情報が得られた時刻あるいは時刻範囲であり、観測範囲は観測された地域の範囲を緯度経度で示したものである。
【００２２】
監視部２は、監視データベース８の観測範囲と観測時刻が、統合指示部１から通知された監視指示に含まれる模擬対象地域、および規則推定時刻列に含まれる時刻または模擬開始時刻と合致する観測情報を取得する。監視部２は、取得した観測情報の画像データから現状画像を生成し、模擬地域および規則推定時刻または模擬開始時刻と共に現状データ生成部３に通知する。
【００２３】
現状データ生成部３は、監視部２から通知された規則推定時刻列の各時刻および模擬開始時刻における現状画像から、各時刻における現状データを生成し、模擬規則推定部４に通知する（ステップＳＴ１０３）。
監視部２から得られる現状画像の形式は、画素毎に強度情報や位相情報を有する多次元配列データ等である。一方、模擬規則推定部４における規則推定処理および模擬部５における被災範囲予測処理では、画素数よりも数が少なく種類も少ないセル集合を使用する。よって、現状画像から、セル集合に合わせた現状データへの変換が必要になる。例えば現状画像中の一定範囲の画素をまとめて１つのセルとする。セル毎に画素が示す強度情報や位相情報の平均を算出し、その値に基づいてセルを数種類に分類する。これにより、画素よりも数が少なく、種類数も少ないセルの集合からなる現状データに変換する。
このようにして得られる現状データの例を図８の４０１に示す。図に示すように、模擬対象地域が１１〜１５，２１〜２５，３１〜３５，４１〜４５で示す複数セルの集合として表される。各セルは３種類のセルに分類されており、種類に応じたセル特性値を有する。また、個々のセルは状態値を有する。各セルは、状態値の値に応じていくつかの状態に分類することができる。ここでは、燃焼状態と非燃焼状態の２状態に類別している。
【００２４】
模擬規則推定部４は、統合指示部１から規則推定指示および事象獲得指示の通知を受け、かつ現状データ生成部３から現状データを受け取ったら模擬規則推定処理を開始する（ステップＳＴ１０４）。
【００２５】
まず、模擬規則推定部４は、受け取った事象獲得指示を事象管理部６に通知する。
事象管理部６は、事象獲得指示を受け取ると、事象データベース９に格納されたレコードの中から、事象内容が、事象獲得指示に含まれる事象と合致するレコードを全て抽出し、模擬規則推定部４に事象集合として出力する。
【００２６】
模擬規則推定部４は、規則推定指示、規則推定時刻列に含まれる各時刻での現状データ、および事象集合に基づいて、模擬規則を生成する。模擬規則は、図８の４０１に示すようなセル集合に対する規則であり、その構成要素として模擬範囲、模擬時間、隣接範囲、セル種間伝播確率、セル特性値、事象毎の伝播阻害確率、共通バイアスを含む。
【００２７】
模擬範囲は、同一の模擬規則を適用する範囲であり、規則推定指示で与えられた模擬地域内の部分にあたる。図８の４０１，４０２では、セル１１〜１５、セル２１〜２５、セル３１〜３５，セル４１〜４５の２０個のセルで示される範囲を１つの模擬範囲としている。模擬時間は、同一の模擬規則を適用する時間である。図８の４０１，４０２では１つの模擬時間に属する２時刻での状態を表している。
隣接範囲は、隣接セルの範囲を規定する数値である。平面上に格子状に並んだセルにおいて、格子の単位長さ（セル幅）を１とすると、隣接範囲の値が１のときには上下左右のセルが隣接セルとなる。隣接範囲の値が１．５であれば、４つの辺で接する上下左右のセルに加え、４つの頂点で接するセルも隣接セルとなる。例えば４０１において、セル２２の隣接セルは、隣接範囲１のときは２１、１２、２３、３２であり、隣接範囲１．５のときには１１、１２、１３、２１、２３、３１、３２、３３の８個のセルが隣接セルとなる。
セル種間伝播確率は、例えば燃焼の場合を例にとって説明すると、あるセルが燃焼状態にある時に、そのセルの隣接セルに対し、単位時間あたり状態値に１単位の影響を及ぼす確率である。なお、ここでは状態値０を非燃焼状態とする。
セル特性値とは、セル種類によって決まる値であり、燃焼の場合を例にとって説明すると、セルの状態値が初期状態から何単位の影響を受けた時に燃焼状態になるかを表す。例えば、あるセルＡの隣接セルのうち燃焼状態にあるセルがセルＢ１つだった場合、セルＡのセル種別特性値をａ、伝播確率をｐとすると、セルＢからセルＡに延焼するには最低ａ単位時間、平均ａ／ｐ単位時間かかることになる。また、セルＢからセルＡにａ単位時間内に延焼する尤度は、ｐのａ乗となる。燃焼状態になったセルの燃焼状態値は、セル種類毎に所与の燃焼状態値に設定され、セル種類毎かつ燃焼状態値毎に所与の減少率に従って減少していく。また、複数の模擬範囲もしくは模擬時間に属するセルをまたがる伝播確率は、例えば伝播を受ける側のセルに適用された模擬規則を用いる。
伝播阻害確率とは、事象によって決まる値であって、該事象に属するセルの隣接セルの燃焼状態値に、単位時間あたり−１単位の影響を及ぼす確率である。例えば、あるセルＣが燃焼状態にあり、その隣接セルの中ではセルＤにだけ消火事象が設定されていたとする。セルＣの燃焼状態値をｃ、伝播阻害確率をｑとすると、セルＣがセルＤの効果で鎮火するまでに最低ｃ単位時間、平均ｃ／ｑ単位時間かかることになる。
共通バイアスとは、模擬範囲全体にわたって、例えば風によって発生する伝播確率の偏り（バイアス）を指す。共通バイアスは方向と大きさからなる。例えば東の方向に大きさｒの共通バイアスがある場合、セルＢからセルＡへの伝播確率をｐとすると、例えばセルＢから見てセルＡが東にあればｐ＋ｒ、西にあればｐ−ｒの伝播確率とする。ここでは、伝播確率が０を下回る場合には０、１を上回る場合には１とする。
【００２８】
隣接範囲推定部４１において隣接範囲が算出され、伝播確率推定部４２においてセル種間伝播確率が算出され、セル特性値推定部４３においてセル特性値が算出され、伝播阻害確率推定部４４において各事象の伝播阻害確率が算出される。上記各部は、共通推定部４６における共通推定処理よって各数値の算出を行う。
【００２９】
共通推定部４６における共通推定処理について説明する。
まず、変数初期値決定部４６１において、推定対象変数に初期値を設定し、伝播尤度の暫定最適値を設定する。初期値の設定方法は、例えば各変数のとりうる範囲内で乱数によって定めるようにしてもよい。あるいは、前回の推定処理で得られた値を保持しておき、初期値としてもよい。変数は、隣接範囲推定部４１における隣接範囲推定処理の場合には隣接範囲、伝播確率推定部４２におけるセル種間伝播確率推定処理の場合には各セル種間の伝播確率、セル特性値推定部４３におけるセル特性値推定処理の場合には各セル種のセル特性値、伝播阻害確率推定部４４における伝播阻害確率推定処理の場合には各事象の伝播阻害確率である。
【００３０】
次に、伝播尤度算出部４６２において、隣接網作成部４６２１による隣接網作成処理、外縁抽出部４６２２による外縁抽出処理、および外縁間伝播尤度算出部４６２３による外縁間伝播尤度算出処理が順次実行される。これらの各部の処理により、セル間での状態変化の伝播をベースとした模擬規則における伝播尤度を最大にする諸変数を決定する。以下、燃焼事象における模擬規則を例にとり説明する。
【００３１】
まず、隣接網作成部４６２１において、設定されている隣接範囲の値に基づいて、各セルの隣接セルを抽出し、それらを記述した隣接網を作成する。
【００３２】
次に、外縁抽出部４６２２において外縁のセルを抽出する。ここで外縁とは、例えば燃焼事象を反映した現状データの場合には、非燃焼セルに隣接した燃焼セルの集合（燃焼境界）と、燃焼セルに隣接した非燃焼セルの集合（非燃焼境界）を意味する。
模擬規則推定部４は、事象管理部６から取得した事象集合中の各事象を、各事象の発生時刻以降の現状データにおける該当する位置のセルに設定する。外縁抽出部４６２２は、各時刻における現状データに対して外縁抽出処理を実施する。
図８の４０２は、現状データ４０１の時刻の一定時間後の現状データであり、燃焼が進んだ様子を示している。現状データ４０２において、隣接範囲が１のときの燃焼境界は、３１、３２、３３、２３、１４である。また、非燃焼境界は、４１、４２、４３、３４、２４、１５である。
また、図９の５０２は、現状データ５０１の時刻の一定時間後の現状データであり、燃焼が進んだ様子を示している。現状データ５０１の時刻でセル１３に消火事象が設定されており、現状データ５０２の時刻においてセル１２が燃焼終了状態になっている。現状データ５０２において、隣接範囲１のとき、燃焼境界は２２、３１、３２、３３であり、非燃焼境界は、４１、４２、４３、３４、２３である。
【００３３】
次に、外縁間伝播尤度算出部４６２３において、外縁抽出部４６２２で抽出された各時刻における現状データの燃焼境界および非燃焼境界について、ある時刻の燃焼境界からその後の時刻の燃焼境界に至る各セルの各時刻における伝播尤度、ならびにある時刻の燃焼境界からその後の時刻の非燃焼境界への各セルの各時刻における伝播尤度を算出する。
セルｉの時刻ｔにおける伝播尤度をＺ（ｉ，ｔ）、セルｉの隣接セルの数をＮ（ｉ）、セルｉからセルｊへの伝播確率をＰ（ｉ，ｊ）、セルｊのセル特性値をＱ（ｊ）、セルｊの隣接セルからの伝播阻害確率をＡ（ｊ）とする。時刻範囲ｔ＋１からｔ＋Ｔまでの間の時刻τにおける、セルｊの隣接セルｉからの伝播パスの尤度は、
Ｚ（ｉ，ｔ）Ｐ（ｉ，ｊ）（１−Ａ（ｊ））τ−１
と表すことができる。このようなパスはＴＮ（ｊ）本存在する。よってセルｊの時刻ｔ＋Ｔにおける伝播尤度Ｚ（ｊ，ｔ＋Ｔ）は、上記の中からＱ（ｊ）本の伝播パスで伝播が起きる尤度を、ＴＮ（ｊ）本からＱ（ｊ）本選択する組合せ全部にわたって合計したものとすることができる。このように、各セルの各時刻における伝播尤度は、燃焼境界に属するセルから反復的に計算することができる。
【００３４】
ここで、図１０を用いて、伝播パスの概念を具体的に説明する。
図において、時刻０ではセル１１のみが燃焼状態であり、セル２１、２２、３１は非燃焼状態である。時刻１および時刻２においては、セル１１は燃焼状態であるが、その他のセルの状態は不明である。時刻３においては、セル１１およびセル２１が燃焼状態であり、セル２２およびセル３１が非燃焼状態である。また、セル１１と２１、セル２１とセル２２、セル２１とセル３１はそれぞれ隣接セルである。この場合、図中に矢印で示したような隣接セルへの延焼可能性が生じる。この延焼可能性の矢印と、同一セルを結んだ太線とによって、時刻０の燃焼セル１１から時刻３の燃焼セル２１への伝播パスとして３通りが考えられる。また、時刻０の燃焼セル１１から時刻３の２つの非燃焼セル２２，３１への伝播パス１本ずつをそれぞれ考えることができる。時刻０のセル１１の伝播尤度を１、セル１１からセル２１への伝播確率を０．３、セル２１のセル特性値を１とすると、セル１１からセル２１への３本の伝播パスによる時刻３までの伝播尤度は、１×０．３×３＝０．９
と計算できる。
また、セル特性値を２とした場合には、３本の伝播パスのうち２本によって伝播が成立することになり、その組合せは３通りなので、
１×０．３×０．３×３＝０．２７
となる。
また、セル１１に伝播阻害確率が設定されている場合は、セル１１から別の時刻でのセル１１への線についても確率が与えられていると考える。セル１１の伝播阻害確率を０．６とすると、セル１１からセル２１への３本の伝播パスによる時刻３までの伝播尤度は、２本で伝播の起きる尤度の総和は、
（０．３×０．３×０．４）＋（０．３×０．３×０．１６）
＋（０．３×０．４×０．３×０．１６）＝０．０５６１６
と計算できる。また、時刻３までにセル１１からセル３１へ伝播しない尤度は、時刻２におけるセル２１の伝播尤度と、セル２１からセル３１へ伝播しない確率との積となる。伝播しない確率は、１から伝播確率を引いた値となる。
【００３５】
図８に示す例を用いて、上述の方法を説明する。４０３は、現状データ４０１の時刻における燃焼境界セル１１およびセル１２から、現状データ４０２における燃焼境界のうちセル３１とセル３２への、伝播パスの例を示したものである。セル３１へは、１１→２１→３１、１２→２２→２１→３１、１２→２２→３２→３１、１２→１３→２３→３３→３２→３１といったパスが可能である。また、４０４は、現状データ４０２の時刻における燃焼境界から、同時刻における非燃焼境界への伝播パスの例を示したものである。
現状データ４０１の時刻の燃焼境界セル１１およびセル１２から現状データ４０２の時刻のセル３１への各伝播パスの伝播尤度は、上述したように、隣接セルからの伝播確率と自セルのセル特性値によって、現状データ４０２へ至るまでの間の時刻ｔにおける伝播尤度と、その後の時刻ｔ＋Ｔまでにさらに隣接セルへ伝播する尤度の積で表すことができる。また、現状データ４０１の時刻から現状データ４０２の時刻までに燃焼境界セル１１およびセル１２からセル４１に伝播しない尤度は、現状データ４０２へ至る途中の時刻ｔまでにセル３１へ伝播する尤度と、その後現状データ４０２の時刻までにセル３１からセル４１へ伝播しない尤度の積で表すことができる。
【００３６】
さらに、図９の例について説明する。５０３は、現状データ５０１の時刻における燃焼境界セル１１およびセル１２から、現状データ５０２における燃焼境界のうちセル３１への、伝播パスの例を示したものである。セル３１へは、１１→２１→３１、１２→２２→２１→３１、１２→２２→３２→３１といったパスが可能である。５０４は、現状データ５０２の時刻における燃焼境界から、同時刻における非燃焼境界への伝播パスの例を示したものである。
現状データ５０１の時刻から現状データ５０２の時刻までに燃焼境界セル１１およびセル１２からセル３１への各伝播パスの伝播尤度は、該パスに沿って、上記のように時刻ｔの伝播尤度、隣接セルからの伝播確率と自セルの伝播阻害確率およびセル特性値から伝播尤度を算出し、それを反復することによって計算できる。また、５０１の時刻から５０２の時刻までに燃焼境界セル１１およびセル１２からセル４１に伝播しない尤度は、３１へ伝播する尤度と３１から４１へ伝播しない尤度の積で表すことができる。
【００３７】
伝播尤度算出部４６２における伝播尤度算出が終了すると、終了判定部４６３は、推定対象変数の値および伝播尤度を記録する。各伝播パスの伝播尤度の総合値、例えば対数の総和が暫定最適値より大きい場合には、該伝播尤度ならびに該変数値を暫定最適値とする。終了判定部４６３は変数変更部４６４に推定対象変数の値を変更するように通知する。変更後の推定対象変数は伝播尤度算出部４６２に供給され、再度伝播尤度算出処理が行われる。
一方、終了判定部４６３は伝播尤度の暫定最適値が予め定められた回数の伝播尤度算出処理を行っても改善しなくなった場合は反復処理を停止し、共通推定部４６の処理を終了する。
【００３８】
変数変更部４６４は、終了判定部４６３から推定対象変数を変更するように通知されると、推定対象変数の値を変更し、再び伝播尤度算出部４６２において伝播尤度算出処理を行う。ここでの変数の変更方式は、推定対象変数の種類によって決まった方式を採用する。
例えば、隣接範囲の値を変数とする隣接範囲推定部４１は、変数初期値決定部４６１での初期値をセル幅と等しくし、変数変更部４６４での変更方式は、例えばセル幅の半分の刻みで増加させるようにして、共通推定部４６において尤度最大となる隣接範囲を算出する。
また、各セル種間の伝播確率を変数とする伝播確率推定部４２は、変数初期値決定部４６１での初期値を例えば０．５とし、変数変更部４６４での変更方式を０から１の範囲の二分探索として、共通推定部４６において尤度最大となる伝播確率を算出する。
また、各セル種のセル特性値を変数とするセル特性値推定部４３は、変数初期値決定部４６１での初期値をセル種類毎の所与の値とし、変数変更部４６４での変更方式を例えば幅１のランダム方向探索として、共通推定部４６において尤度最大となるセル特性値を算出する。
また、各事象内容の伝播阻害確率を変数とする伝播阻害確率推定部４４は、変数初期値決定部４６１での初期値を事象内容毎の所与の値とし、変数変更部４６４での変更方式を例えば幅１のランダム方向探索として、共通推定部４６において尤度最大となる伝播阻害確率を算出する。
【００３９】
次に、模擬規則推定の補助的な処理として、均一領域推定部４５１、均一時間帯推定部４５２、風向風速推定部４５３、セル種類修正部４５４、時間単位修正部４５５、および事象推定部４５６において、共通推定部４６を用いて各変数の推定を行う（ステップＳＴ１０５）。なお、上記各部の処理は模擬規則推定のための補助的な処理であるので、全て実行してもよいし、一部のみの実行でもよい。あるいは、いずれも実行しなくてもよい。
【００４０】
均一領域推定部４５１では、模擬地域を所与の数の模擬範囲に分割する。各模擬範囲内で模擬規則を同一とし、異なる模擬範囲における模擬規則は別の規則として推定する。各セルの所属する模擬範囲を各セルの変数とし、共通推定部４６において尤度最大となる所属模擬範囲を求める。変数初期値決定部４６１での初期値は、例えば模擬範囲の数だけランダムにセルを選び、各セルの初期値は一番距離の近いセルとする。変数変更部４６４での変更方式は、例えば隣接セルのいずれかをランダムに選択し、その模擬範囲と同一とする。このように、模擬地域の分割により、地域性の変化・時間経過に伴う状況の変化に応じた模擬規則を自動的に推定する。
【００４１】
均一時間帯推定部４５２では、模擬開始時刻から模擬停止時刻までの時間を、所与の数の模擬時間に分割する。各模擬時間内での模擬規則を同一とし、模擬時間の違う模擬規則は別の規則を推定する。模擬時間の境界時刻を変数とし、共通推定部４６において尤度最大となる境界時刻を求める。変数初期値決定部４６１での初期値は、例えば模擬時間の境界時刻数だけランダムに時刻を選ぶ。変数変更部４６４での変更方式は、例えば各境界時刻をランダムな幅だけ増減させる。このように、模擬時間帯の分割により、地域性の変化・時間経過に伴う状況の変化に応じた模擬規則を自動的に推定する。
【００４２】
風向風速推定部４５３では、共通バイアスの方向と大きさを変数として、共通推定部４６において尤度最大となる方向と大きさを求める。変数初期値決定部４６１での初期値は例えばランダムに選ぶ。変数変更部４６４での変更方式は、例えば方向および大きさをランダムな幅だけ増減させる。このように、共通バイアスを推定することにより、模擬地域・模擬時間帯にわたって均一な模擬規則成分を自動的に推定する。
【００４３】
セル種類修正部４５４では、所与の数のセルの種類指定に誤指定があると仮定したときに、どのセルが誤っているかを推定する。誤指定セルの位置を変数とし、共通推定部４６において尤度最大となる誤指定セルを求める。変数初期値決定部４６１での初期値は例えばランダムに選ぶ。変数変更部４６４での変更方式も例えばランダムな選択とする。このように、セル種類の誤りを推定することにより、セル種類が誤っていても自動的に模擬規則を推定する。
【００４４】
時間単位修正部４５５では、１回の模擬規則適用の時間単位を推定する。時間単位を変数とし、共通推定部４６において尤度最大となる時間単位を求める。変数初期値決定部４６１での初期値は所与の値とし、変数変更部４６４での変更方式は、例えばランダムな幅だけ増減させる。このように、模擬規則適用のための適正な時間単位を推定することにより、隣接範囲を小さくすることができ、模擬規則を高速に推定できる。
【００４５】
事象推定部４５６では、所与の事象内容について欠落があったと仮定してその位置および時刻を推定する。事象の位置および時刻を変数とし、共通推定部４６において伝播尤度最大となる位置および時刻を求める。変数初期値決定部４６１での初期値は例えばランダムに選ぶ。変数変更部４６４での変更方式も例えばランダムな選択とする。このように、事象データベース９にない事象を推定することにより、事象データベース９に欠落があっても自動的に模擬規則を推定することができる。
【００４６】
次に、終了判定部４７において、隣接範囲推定部４１、伝播確率推定部４２、セル特性値推定部４３、伝播阻害確率推定部４４、均一領域推定部４５１、均一時間帯推定部４５２、風向風速推定部４５３、セル種類修正部４５４および時間単位修正部４５５の各部における模擬規則の変更幅が所与の幅以内のとき、あるいは事象推定部４５６における事象情報の変更が所与の反復回数の間なかったとき、あるいは所与の実行時間を経過したときに模擬規則推定処理の終了と判定し、ステップＳＴ１０７へ進む。終了と判定されなかった場合にはステップＳＴ１０４に戻る（ステップＳＴ１０６）。
【００４７】
模擬規則格納部４８は、推定された模擬規則を災害監視装置１００内の記憶部（図示せず）あるいは外部記憶装置（図示せず）へ格納すると共に模擬部５へ通知する（ステップＳＴ１０７）。
【００４８】
模擬部５は、統合指示部１からの模擬指示および事象獲得指示を受け、現状データ生成部３から模擬指示に含まれる模擬対象地域の模擬開始時刻における観測情報に基づいて生成された現状データを受け、かつ模擬規則推定部４から模擬規則を受け取ったら、被災範囲予測処理を行う（ステップＳＴ１０８）。
【００４９】
模擬部５は、統合指示部１からの模擬指示に含まれる事象獲得指示を事象管理部６に送付し、その応答として事象集合を獲得する。現状データ生成部３から受け取った現状データと、模擬規則推定部４から受け取った模擬規則に基づき、事象管理部６から受け取った事象集合を参照して、模擬指示に含まれる模擬時刻列に含まれる時刻での被災範囲予測を行い、各時刻での模擬データを生成する。
【００５０】
ここで、計画データとは、現状データと同様の、模擬対象地域をセル集合によって表したものである。図１１に計画データの例を示す。計画データ７０１は、模擬開始時刻から模擬停止時刻までの間のある時刻における模擬範囲の例を示している。模擬データ７０１においては、模擬開始時刻における現状データが示され、さらに事象管理部６から取得した事象集合の例として発火事象がセル１１およびセル１２に、消火事象がセル１３に示されている。７０２はこの模擬範囲における模擬規則の例である。発火事象の伝播確率が３つのセル種類について示され、消火事象の伝播阻害確率が示されている。また、模擬データ７０３は、模擬データ７０１の時刻より後の模擬時刻計画に含まれる時刻における模擬データであり、模擬データ７０１を基準に模擬規則７０２に従って模擬を実施した結果得られたものである。模擬データ生成は、従来技術に関連して述べたセルラーオートマトンの原理に基づいて行われる。模擬データ７０３では、セル２１、２２、３１、３２、３３は燃焼状態にあり、セル１１および１２は燃焼終了状態にあることが示されている。
模擬部５は、被災範囲予測によって得られた各時刻での模擬データを統合表示部７に送付する。
【００５１】
統合表示部７は、模擬部５から受領した模擬データを、位置や時刻の対応が明確になるように加工して、地図、現状画像、事象などのその他の情報と共に時刻別に、あるいは時系列に沿って表示する（ステップＳＴ１０９）。
図１２に、統合表示部７による出力画面の例を示す。図に示すように、出力画面１７１は、現状画像を表示するための現状画像表示領域１７２、表示されている現状画像の時刻を表示する画像時刻表示領域１７３、模擬データを表示する模擬データ表示領域１７４、表示されている模擬データの時刻を表示する模擬時刻表示領域１７５、事象集合および事象計画を表示する事象位置表示領域１７６、事象時刻表示領域１７７を含む。
【００５２】
次に、図１３のフローチャートを用いて、この発明の実施の形態１による計画策定処理について説明する。計画策定処理では、システム利用者が災害地域の一定時間後のあるべき状況を計画データとして設定し、その状況を実現するために発生させる事象（以下、計画事象という。）の的確な発生位置および発生時刻を推定する。
【００５３】
例えば、ある地域で発生した火災を完全に鎮火させるための対策計画を策定する場合、システム利用者は計画事象として例えば消火を選択し、計画データとして、その地域が非燃焼状態もしくは燃焼終了状態であることを示すデータを作成する。なお。計画データは現状データ生成部３が生成する現状データと同一形式のデータである。
また、ガソリンスタンドや化学薬品工場など、危険度の高い地域を防災重点地域として火災から防護する計画を策定する場合には、計画事象として消火を選択し、それらの防災重点地域が非燃焼状態であることを示す計画データを作成する。この場合、例えば、火災に対する危険度の高い地点を防災重点地域として予めリストに保持しておき、火災が発生した時には自動的にそれらの地域の消火を実現する対策計画を策定するようにしてもよい。また、消火手段が複数ある場合には、それぞれ別の消火事象として登録しておき、適切な方法を選択するようにしてもよい。
【００５４】
システム利用者は、システム利用者が入力画面１６１上で、「計画策定」ボタン１６５を選択することにより、地域指定領域１６６に対策地域の計画データ、時刻指定領域１６７に計画データで示される状態を実現させる時刻、事象指定領域１６８に計画事象を入力することができる。また、計画データと共に、計画策定に用いる模擬規則を指定する。各領域への入力は、ＧＵＩ等により支援される。
【００５５】
統合指示部１は、監視部２に通知される監視指示と、模擬規則推定部４に通知される計画策定指示を作成する（ステップＳＴ２０１）。監視指示には、システム利用者の入力操作によって指定された対策地域および事象推定開始時刻が含まれる。計画策定指示には、対策地域と計画データ、計画事象、および指定模擬規則が含まれる。また、模擬規則推定部４へは、計画策定指示と共に、指定された計画事象を含む事象獲得指示が通知される。また、統合指示部１は、事象データベース９に計画事象をユーザ定義事象として登録するよう事象管理部６に通知する。
【００５６】
統合指示部１から監視部２へ監視指示が通知されると、監視部２は監視データベース８を参照し、レコードの観測範囲と観測時刻が監視指示に含まれる対策地域、事象推定開始時刻と一致する観測情報が取得される。あるいは、監視部２は観測機器を制御して、当該地点における当該時刻の観測情報を取得する。監視部２は、取得した観測情報から現状画像を生成し、現状データ生成部３に通知する（ステップＳＴ２０２）。
【００５７】
現状データ生成部３は、監視部２から通知された現状画像を模擬規則推定部４の要求する形式に加工し、事象推定開始時刻における現状データとして模擬規則推定部４に通知する（ステップＳＴ２０３）。
【００５８】
模擬規則推定部４は、統合指示部１から計画策定指示および事象獲得指示が通知され、現状データ生成部３から現状データが通知されると、計画策定処理を開始する（ステップＳＴ２０４）。
まず、模擬規則推定部４は、受け取った事象獲得指示を事象管理部６に通知する。事象管理部６は、模擬規則推定部４から事象獲得指示を受け取ると、事象データベース９に格納されたレコードの中から、事象内容が、事象獲得指示に含まれる事象と合致するレコードを全て抽出し、模擬規則推定部４に事象集合として出力する。
【００５９】
模擬規則推定部４は、災害監視装置１００内の記憶部等に記憶されている模擬規則の中から、計画策定指示に含まれる指定模擬規則に該当するものを取得する。上述した事象推定部４５６において、取得した模擬規則を用いて事象推定を実施し、指定された計画事象の発生地点および時刻を推定する（ステップＳＴ２０４）。
【００６０】
終了判定部４７は、事象情報の変更が所与の反復回数の間なかったとき、または所与の実行時間を経過したときには計画策定処理を終了し、ステップＳＴ２０６に進む。終了しない場合にはステップＳＴ２０４に戻る（ステップＳＴ２０５）。
【００６１】
模擬規則推定部４において推定した事象を、事象計画を事象管理部６および統合表示部７に送付する。事象管理部６は、事象データベース９に推定した事象の情報をシステム推定事象として登録する。また、事象計画を統合表示部７に送付する（ステップＳＴ２０６）。
【００６２】
統合表示部７において、地図、現状画像、事象など他の情報との、位置や時刻の対応が明確になるように、模擬規則推定部４から受領した事象計画を加工して表示する（ステップＳＴ２０７）。
【００６３】
以上のように、この実施の形態１によれば、統合指示部１がシステム利用者の入力に基づいて作成した命令に従って、模擬規則推定部４が観測装置が収集した災害現場の観測データに基づいて模擬規則を推定し、推定された規則を用いて模擬部５が被災範囲予測を行うようにしたので、災害現場の現状に即した模擬規則を迅速に取得することができ、より精度の高い被災範囲予測を行うことができるという効果がある。
【００６４】
また、被災地域を複数セルの集合として表し、セル間での事象の伝播に基づく被災範囲予測における模擬規則について、伝播発生の尤度を最大化するように模擬規則の諸変数を決めるようにしたので、セル間での事象の伝播に基づく模擬規則を自動推定することが可能である。
【００６５】
また、セル間での伝播をベースとした模擬における模擬規則について、伝播尤度を外縁間のパスを列挙することによって算出し、その最大化を基準として模擬規則の諸変数を決めることにより、セル間の伝播をベースとした模擬における模擬規則の、自動推定を可能にする。
【００６６】
また、ある地域について一定時間後にあるべき状況を設定し、その状況が実現するために発生させる事象の的確な位置および発生時刻を推定することにより、所望の結果を得るための災害対策計画を策定することができる。
【００６７】
【発明の効果】
以上のように、この発明によれば、模擬規則推定部が観測装置から得られる災害現場の観測データに基づいて模擬規則を推定し、推定された規則を用いて模擬部が被災範囲予測を行うようにしたので、災害現場から得られる監視データに基づいて自動的に模擬規則を決定し、被災範囲予測を行う災害監視装置を得られるという効果がある。
【図面の簡単な説明】
【図１】この発明の実施の形態１による、災害監視装置の構成を示すブロック図である。
【図２】この発明の実施の形態１による、模擬規則推定部の詳細を示すブロック図である。
【図３】この発明の実施の形態１による、共通推定部の詳細を示すブロック図である。
【図４】この発明の実施の形態１による、模擬規則推定処理、および被災範囲予測処理のフローチャートである。
【図５】この発明の実施の形態１による、入力画面の例を示す図である。
【図６】この発明の実施の形態１による、事象データベースの内容の例を示す図である。
【図７】この発明の実施の形態１による、監視データベースの内容の例を示す図である。
【図８】この発明の実施の形態１による、模擬規則の推定処理を説明する図である。
【図９】この発明の実施の形態１による、模擬規則の推定処理を説明する図である。
【図１０】この発明の実施の形態１による、伝播尤度算出処理を説明する図である。
【図１１】この発明の実施の形態１による、被災範囲予測処理を説明する図である。
【図１２】この発明の実施の形態１による、統合表示部による出力画面の例を示す図である。
【図１３】この発明の実施の形態１による、計画策定処理のフローチャートである。
【符号の説明】
１  統合指示部、２  監視部、３  現状データ生成部、４  模擬規則推定部、５  模擬部、６  事象管理部、７  統合表示部、８  監視データベース、９  事象データベース、４１  隣接範囲推定部、４２  伝播確率推定部、４３  セル特性値推定部、４４  伝播阻害確率推定部、４６  共通推定部、４７  終了判定部、４８  模擬規則格納部、１００  災害監視装置、４５１  均一領域推定部、４５２  均一時間帯推定部、４５３  風向風速推定部、４５４  セル種類修正部、４５５  時間単位修正部、４５６  事象推定部、４６１  変数初期値決定部、４６２  伝播尤度算出部、４６３  終了判定部（変数推定終了判定部）、４６４  変数変更部、４６２１  隣接網作成部、４６２２  外縁抽出部、４６２３  外縁間伝播尤度算出部（外縁間パス算出部、伝播尤度算出部）。
【図１】
2004295633.tif 000003
【図２】
2004295633.tif 000004
【図３】
2004295633.tif 000005
【図４】
2004295633.tif 000006
【図５】
2004295633.tif 000007
【図６】
2004295633.tif 000008
【図７】
2004295633.tif 000009
【図８】
2004295633.tif 000010
【図９】
2004295633.tif 000011
【図１０】
2004295633.tif 000012
【図１１】
2004295633.tif 000013
【図１２】
2004295633.tif 000014
【図１３】
2004295633.tif 000015