JP2012132824A - 拡散物質の発生源推定装置および発生源推定方法 - Google Patents

Abstract

【課題】より柔軟且つ簡便に発生源を推定し得る拡散物質の発生源推定装置および発生源推定方法を提供する。
【解決手段】観測器５−ｉ（ｉ＝１〜ｎ）の位置情報および計測した濃度情報を観測情報入手部２１によって入手し、仮想格子設定部２２により仮想格子上で仮想放出地点Ｐｏｊ（ｊ＝１〜ｍ）を設定し、影響関数算出部２４により影響関数Ｄ_ｉｊを算出し、残差ノルム算出部２５により、仮想放出地点Ｐｏｊ毎に、各観測器５−ｉの濃度情報と、各観測器５−ｉの該仮想放出地点Ｐｏｊに対する影響関数Ｄ_ｉｊと該仮想放出地点Ｐｏｊの放出強度ｑ_ｊとの積と、の差の平方和である残差ノルムＲ（ｑ_ｊ）を算出し、推定部２６により、算出された全仮想放出地点の残差ノルムの内、残差ノルムが最小となる仮想放出地点を放出地点と推定する。
【選択図】図１

Description

本発明は、拡散物質の発生源を推定する拡散物質の発生源推定装置および発生源推定方法に関するものである。

プラント施設（火力発電所、ゴミ焼却施設、化学プラント等）での事故等による汚染物質の放出や、テロ等による毒性ガス等の放出に対して、その発生源を同定して即座に対処する必要性から、事故または事件の現場情報（濃度計測値等）から発生源情報（放出地点の位置および発生量）を推定する拡散物質の発生源推定装置および発生源推定方法について各種技法が提案されている。

例えば、非特許文献１の「センシング情報を用いた応答係数法に基づく環境影響物質の発生源同定法」では、仮想の放出点および時刻等の情報に基づき観測位置での影響を評価し、変分原理によって誤差を最小とする放出点を求め、該放出点における時刻および放出量を同定する手法が提案されている。

石田義洋，加藤信介，樋山恭助，"センシング情報を用いた応答係数法に基づく環境影響物質の発生源同定法（第２報）"，平成２１年度空気調和・衛生工学学術講演会講演論文集，（２００９年９月）

しかしながら、上述した非特許文献１に開示された技術においては、観測点の数が仮想の放出点の数より多くなければならないという制約があり、より柔軟に拡散物質の発生源を推定し得る手法が望まれていた。

本発明は、上記従来の事情に鑑みてなされたものであって、より柔軟且つ簡便に発生源を推定し得る拡散物質の発生源推定装置および発生源推定方法を提供することを目的としている。

上記課題を解決するため、本発明は以下の手段を採用する。
本発明に係る拡散物質の発生源推定装置は、複数個の観測器からの情報に基づき、ガス発生源情報を推定する拡散物質の発生源推定装置であって、前記観測器の位置情報および計測した濃度情報を入手する観測情報入手手段と、相互の離間距離が一定の仮想格子上で各格子線が交差する位置を仮想放出地点として設定する仮想格子設定手段と、拡散モデルを用いて、前記観測器と前記仮想放出地点との相対位置、並びに相対時刻に応じて定まる影響関数を算出する影響関数算出手段と、前記仮想放出地点毎に、各観測器の濃度情報と、各観測器の該仮想放出地点に対する影響関数と該仮想放出地点の放出強度との積と、の差の平方和である残差ノルムを算出する残差ノルム算出手段と、算出された全仮想放出地点の残差ノルムの内、残差ノルムが最小となる仮想放出地点を放出地点と推定する推定手段と、を備えることを特徴とする。

本発明によれば、設定した仮想放出地点毎の残差ノルムを評価して、該残差ノルムが最小となる放出強度を求め、該放出強度の仮想放出地点を放出位置とし、また該放出強度を該放出位置からの放出量として推定するので、「観測地点の数≧仮想放出地点の数」という制約を持つことなく放出地点を推定することができ、より柔軟且つ簡便に発生源を推定し得る拡散物質の発生源推定装置を実現することができる。

また、本発明は、上記記載の拡散物質の発生源推定装置において、前記影響関数算出手段は、数値拡散計算により前記影響関数を算出することを特徴とする。

例えば、平地一様流場では拡散モデルを用いて影響関数を算出し、また複雑気流場では数値拡散（シミュレーション）計算により影響関数を算出することで、様々な地形における拡散に対してより精度良く発生源を推定することができる。

また、本発明は、上記記載の拡散物質の発生源推定装置において、仮想格子設定手段は、前記推定手段で推定された放出地点を含み、且つ相互の離間距離がより短い仮想格子上で各格子線が交差する位置を仮想放出地点として再設定することを特徴とする。

本発明によれば、格子線の相互の離間距離が相対的に長い仮想格子から相対的に短い仮想格子へと絞り込みながら仮想放出地点を再設定して発生源を推定するようにすれば、離間距離を最短とした１面の仮想格子上で仮想放出地点を設定して発生源を推定する場合と比べて、１面の仮想格子上での仮想放出地点数を格段に小さくすることができ、全体的な処理の計算量を抑制してより速い発生源の推定を行うことができる。

また、本発明は、上記記載の拡散物質の発生源推定装置において、仮想放出時刻を設定する仮想放出時刻設定手段、を備え、前記残差ノルム算出手段は、仮想放出時刻別に各仮想放出地点毎の残差ノルムを算出し、前記推定手段は、仮想放出時刻別に算出された全仮想放出地点の残差ノルムの内、残差ノルムが最小となる仮想放出時刻および仮想放出地点をそれぞれ放出時刻および放出地点と推定することを特徴とする。

本発明によれば、放出時刻が不明な場合でも仮想放出時刻設定手段によって仮想放出時刻を設定することで、より柔軟に発生源を推定し得る拡散物質の発生源推定装置を実現することができる。

また、本発明は、上記記載の拡散物質の発生源推定装置において、前記影響関数算出手段は、想定される前記観測器と前記仮想放出地点との相対位置およびまたは相対時刻に応じた影響関数を予め算出してデータベースに保存しておくことを特徴とする。

本発明によれば、データベースを参照して予め算出した影響関数を処理に用いることにより、影響関数を算出する処理分の計算量を排除することができ、全体的な処理の計算量を低減して計算時間を短縮し、より速く発生源の推定を行うことができる。

また、本発明に係る拡散物質の発生源推定方法は、複数個の観測器からの情報に基づき、ガス発生源情報を推定する拡散物質の発生源推定方法であって、前記観測器の位置情報および計測した濃度情報を入手する観測情報入手ステップと、相互の離間距離が一定の仮想格子上で各格子線が交差する位置を仮想放出地点として設定する仮想格子設定ステップと、拡散モデルを用いて、前記観測器と前記仮想放出地点との相対位置、並びに相対時刻に応じて定まる影響関数を算出する影響関数算出ステップと、前記仮想放出地点毎に、各観測器の濃度情報と、各観測器の該仮想放出地点に対する影響関数と該仮想放出地点の放出強度との積と、の差の平方和である残差ノルムを算出する残差ノルム算出ステップと、算出された全仮想放出地点の残差ノルムの内、残差ノルムが最小となる仮想放出地点を放出地点と推定する推定ステップと、を備えることを特徴とする。

本発明によれば、設定した仮想放出地点毎の残差ノルムを評価して、該残差ノルムが最小となる放出強度を求め、該放出強度の仮想放出地点を放出位置とし、また該放出強度を該放出位置からの放出量として推定するので、「観測地点の数≧仮想放出地点の数」という制約を持つことなく放出地点を推定することができ、より柔軟且つ簡便に発生源を推定し得る拡散物質の発生源推定方法を実現することができる。

また、本発明は、上記記載の拡散物質の発生源推定方法において、前記影響関数算出ステップは、数値拡散計算により前記影響関数を算出することを特徴とする。

例えば、平地一様流場では拡散モデルを用いて影響関数を算出し、また複雑気流場では数値拡散（シミュレーション）計算により影響関数を算出することで、様々な地形における拡散に対してより精度良く発生源を推定することができる。

また、本発明は、上記記載の拡散物質の発生源推定方法において、記仮想格子設定ステップは、前記推定ステップで推定された放出地点を含み、且つ相互の離間距離がより短い仮想格子上で各格子線が交差する位置を仮想放出地点として再設定することを特徴とする。

本発明によれば、格子線の相互の離間距離が相対的に長い仮想格子から相対的に短い仮想格子へと絞り込みながら仮想放出地点を再設定して発生源を推定するようにすれば、離間距離を最短とした１面の仮想格子上で仮想放出地点を設定して発生源を推定する場合と比べて、１面の仮想格子上での仮想放出地点数を格段に小さくすることができ、全体的な処理の計算量を抑制してより速い発生源の推定を行うことができる。

また、本発明は、上記記載の拡散物質の発生源推定方法において、仮想放出時刻を設定する仮想放出時刻設定ステップ、を備え、前記残差ノルム算出ステップは、仮想放出時刻別に各仮想放出地点毎の残差ノルムを算出し、前記推定ステップは、仮想放出時刻別に算出された全仮想放出地点の残差ノルムの内、残差ノルムが最小となる仮想放出時刻および仮想放出地点をそれぞれ放出時刻および放出地点と推定することを特徴とする。

本発明によれば、放出時刻が不明な場合でも仮想放出時刻設定ステップによって仮想放出時刻を設定することで、より柔軟に発生源を推定し得る拡散物質の発生源推定方法を実現することができる。

さらに、本発明は、上記記載の拡散物質の発生源推定方法において、前記影響関数算出ステップは、想定される前記観測器と前記仮想放出地点との相対位置およびまたは相対時刻に応じた影響関数を予め算出してデータベースに保存しておくことを特徴とする。

本発明によれば、データベースを参照して予め算出した影響関数を処理に用いることにより、影響関数を算出する処理分の計算量を排除することができ、全体的な処理の計算量を低減して計算時間を短縮し、より速く発生源の推定を行うことができる。

本発明によれば、観測地点数に制約を持つことなく放出地点を推定することができ、より柔軟且つ簡便に発生源を推定し得る拡散物質の発生源推定装置および発生源推定方法を実現することができるという効果を奏する。

本発明の第１実施形態に係る拡散物質の発生源推定装置の構成図である。 拡散現象の線形性を説明する説明図である。 従来法による残差ノルムの算出例を例示する説明図である。 本発明における残差ノルムの算出例を例示する説明図である。 第１実施形態の拡散物質の発生源推定方法を説明するフローチャートである。 第２実施形態の拡散物質の発生源推定方法を説明するフローチャートである。 Ｎ段目の仮想格子を用いて残差ノルムを算出する例を例示する説明図である。 Ｎ＋１段目の仮想格子を用いて残差ノルムを算出する例を例示する説明図である。

以下、本発明の拡散物質の発生源推定装置および発生源推定方法の実施形態について、第１実施形態、第２実施形態の順に図面を参照して詳細に説明する。

〔第１実施形態〕
図１は本発明の第１実施形態に係る拡散物質の発生源推定装置の構成図である。
同図において、本実施形態の拡散物質の発生源推定装置３は、通信インタフェース１１，入力部１３、発生源推定処理部１５、記憶部１７および出力部１９を備えて構成されている。すなわち、発生源推定装置３はいわゆるコンピュータシステムの構成であり、発生源推定処理部１５はＭＰＵ（マイクロプロセッサ）やＤＳＰ（ディジタルシグナルプロセッサ）等のプロセッサで具現される。また、記憶部１７はＲＡＭ（Random Access Memory）、ＲＯＭ（Read Only Memory）、ＨＤＤ（Hard Disk Drive）等の記憶装置で、入力部１３はキーボード、マウス等の入力装置で、出力部１９はディスプレイ、プリンタ等の出力装置でそれぞれ具現される。

また、図１において、発生源推定装置３は、ｎ個の観測器５−１〜５−ｎ（ｎは正整数）からの情報、即ち該観測器の位置情報並びに該観測器が計測した濃度情報および計測時刻情報を、通信インタフェース１１を介して入手する構成である。観測器５−１〜５−ｎは、少なくともその設置場所における大気中の所望ガスの濃度を計測する機能を備え、定期的に該観測器の位置情報並びに該観測器が計測した濃度情報および計測時刻情報を発信する機能も備える。なお、観測器５−１〜５−ｎからリアルタイムに情報入手が可能な場合には、観測器の位置情報および濃度情報のみを受信して、発生源推定装置３側での受信タイミングで計測時刻を代用するようにしても良い。

また、観測器５−１〜５−ｎは固定設置または可動設置の何れでも良く、可動設置の場合には、ＧＰＳ機能を備えて、常時、該観測器の位置情報を取得可能な構成とすれば良い。またさらに、固定設置の場合には、観測器の位置情報を固有の識別コード（機器番号等）で代用することも可能である。この場合、発生源推定装置３側で、観測器の識別コードに対応したテーブル等に基づき観測器の位置情報を導出することになる。

また、観測地点から情報を入手する方法は、図１に例示した構成に限定されない。例えば、民間または公共機関等によるガス濃度観測システムが存在して、各観測地点における濃度データを逐次蓄積していくデータベース機能を備えたものである場合には、通信インタフェース１１を介しインターネット等のネットワーク上に設置されたデータベースにアクセスして、各観測地点の位置データ並びに該観測器が計測した濃度データおよび観測時刻データを入手するようにしても良い。なお、ガス放出に対してその発生源を即座に同定する観点からは、本実施形態の構成（図１）が望ましい。

また、発生源推定装置３の発生源推定処理部１５には、観測情報入手部２１、仮想格子設定部２２、仮想放出時刻設定部２３、影響関数算出部２４、残差ノルム算出部２５および推定部２６を備えている。これら各構成要素は、ＭＰＵやＤＳＰ等のプロセッサ上で実行されるプログラムの機能的まとまりとして具現されるものである。

ここで、発生源推定処理部１５の各構成要素の具体的な機能を説明する前に、本発明における拡散物質の発生源を推定する手法の基本的な考え方について、図２〜図４を参照して説明する。図２は拡散現象の線形性を説明する説明図であり、図３は従来法（非特許文献１）による残差ノルムの算出例を例示する説明図であり、図４は本発明における残差ノルムの算出例を例示する説明図である。

まず、拡散現象の基本的な性質である線形性について説明する。単純な例として、図２（ａ）に示すように、２箇所の放出地点Ｐｏ１およびＰｏ２からの放出を評価地点（観測地点）Ｐｖで観測する場合について考える。なお、放出地点近傍でｘ方向に一様な風が吹いていると仮定し、風に直角な方向をｙ方向（、風に鉛直な方向をｚ方向）とする。

このとき、評価地点Ｐｖにおける濃度は、図２（ｂ）に示すような放出地点Ｐｏ１からの放出による影響と、図２（ｃ）に示すような放出地点Ｐｏ２からの放出による影響と、の和で表される。すなわち、放出地点Ｐｏ１およびＰｏ２の放出強度をそれぞれｑ_１，ｑ_２とし、放出地点Ｐｏ１およびＰｏ２の放出に対する影響関数をそれぞれＤ_１，Ｄ_２とすると、評価地点Ｐｖにおける濃度Ｄは「Ｄ＝ｑ_１・Ｄ_１＋ｑ_２・Ｄ_２」として表すことができる。

このような拡散現象の放出強度に対する線形性により、放出地点が複数（ｍ箇所；ｍは正整数）あるとき、任意の時刻（ｔ）での評価位置（ｘ，ｙ）における濃度Ｄ（ｘ，ｙ，ｔ）は、各放出地点からの放出による影響の和で表され、次式が成立する。

ここに、ｑ_ｊは放出位置（ｘ_ｊ，ｙ_ｊ）の放出強度であり、Ｄ_ｊ（ｘ−ｘ_ｊ，ｙ−ｙ_ｊ，ｔ−ｔ_ｊ）は評価位置と放出位置との相対位置および相対時刻によって決まる影響関数である。

また、複数の観測地点（即ち、ｎ箇所（ｎは正整数）の評価位置）における濃度（Ｄ（ｘ_ｉ，ｙ_ｉ，ｔ）；ｉ＝１，２，…，ｎ）が計測されていれば、次式が成り立つ。

また、放出地点の数ｍが観測地点の数ｎ以下ならば、連立方程式（２）から各放出位置（ｘ_ｊ，ｙ_ｊ）における放出強度ｑ_ｊが得られる。具体的には、式（２）の左辺と右辺の差の平方和である残差ノルムが最小となるように放出強度ｑを決める。この残差ノルムは次式で表される。

また、残差ノルムを最小にする放出強度ｑ_ｊは変分法により次式で表される。

このとき、残差ノルムは次式となる。

ここで、本発明の手法と対比するために、従来（非特許文献１）の手法について、図３を参照して説明する。後で詳述するように仮想格子上での各格子線の交差点を仮想放出地点として、９箇所（ｍ＝９）の仮想放出地点Ｐｏ１〜Ｐｏ９があり、仮想放出地点数（ｍ＝９）よりも大きい１０箇所（ｎ＝１０）の観測地点Ｐｖ１〜Ｐｖ１０があるものとする。また、各仮想放出地点Ｐｏｊ（ｊ＝１〜９）の放出強度はそれぞれｑ_ｊ（ｊ＝１〜９）とする。

この場合、従来手法では、式（３）の残差ノルムは同図中に示した式に展開され、該式について残差ノルムＲを最小にする放出強度ｑ_ｊを変分法等によって求める。また、計算において放出強度ｑ_ｊが負値となった場合には、その仮想放出地点を除去して再度計算を行う。こうして得られた放出強度ｑ_ｊ（ｊ＝１〜９）の内、最大となる放出強度ｑ_ｊの仮想放出地点Ｐｏｊを放出位置として推定する。

しかしながら、式（４）の示すところは、式（４）を満たす全ての放出強度ｑ_ｊ（ｊ＝１〜ｍ）を適用したときに初めて式（３）の残差ノルムを最小化することができる、というのが本来の意味であり、計算で最大となった放出強度ｑ_ｊの仮想放出地点Ｐｏｊのみを適用した式（３）の残差ノルムが、他の仮想放出地点を適用して評価した残差ノルムよりも小さいという保証はない。これは、従来手法が全仮想放出地点からの影響を想定して残差ノルムを計算するものであり、複数の放出源の可能性を想定した式を用いているにもかかわらず、最終的に１点に絞り込んでいるためと考えられる。

これに対して、本発明では、式（４）を解くことはせずに、各放出源（即ち、各仮想放出地点Ｐｏｊ（ｊ＝１〜ｍ））に対する残差ノルムを評価して、該残差ノルムが最小となる放出強度ｑ_ｊを求め、該放出強度ｑ_ｊの仮想放出地点Ｐｏｊを放出位置とし、また該放出強度ｑ_ｊを該放出位置からの放出量として推定する。

具体的に、各仮想放出地点Ｐｏｊ（ｊ＝１〜ｍ）に対する残差ノルムは次式で表される。

また、この残差ノルムを最小にする放出強度ｑ_ｊ並びにそのときの残差ノルムＲ_ｊは、次式で表される。

ここで、図４に例示した本発明における残差ノルムの算出例を参照して説明する。図４では、従来の手法と対比するために、仮想放出地点数および観測地点数を図３と同数としている。また、本発明では仮想放出地点Ｐｏｊ（ｊ＝１〜９）毎に残差ノルムを計算するが、図４では代表例として仮想放出地点Ｐｏ１に対する残差ノルムの算出について例示している。この場合、式（６）の残差ノルムは同図中に示した式に展開される。同様に、他の仮想放出地点Ｐｏｊ（ｊ＝２〜９）に対する残差ノルムを算出して、該残差ノルムが最小となる放出強度ｑ_ｊを求め、該放出強度ｑ_ｊの仮想放出地点Ｐｏｊを放出位置とし、また該放出強度ｑ_ｊを該放出位置からの放出量として推定する。

次に、以上説明した基本的な理論を踏まえて、発生源推定処理部１５の各構成要素（即ち、観測情報入手部２１、仮想格子設定部２２、仮想放出時刻設定部２３、影響関数算出部２４、残差ノルム算出部２５および推定部２６）の具体的な機能について説明する。

まず、観測情報入手部２１は、観測器５−１〜５−ｎからの情報、即ち該観測器の位置情報、該観測器が計測した濃度情報および計測時刻情報を、通信インタフェース１１を介して入手する。入手した各情報は関連づけて記憶部１７の所定領域に保存しておくのが望ましい。

また、仮想格子設定部２２は、相互の離間距離が一定の仮想格子を想定して、該仮想格子の領域上で各格子線が交差する位置を仮想放出地点Ｐｏｊ（ｊ＝１〜ｍ）として設定する。具体的に、図４に示した例では、仮想格子上で９箇所（ｍ＝９＝３×３）の仮想放出地点Ｐｏ１〜Ｐｏ９を設定した。この仮想放出地点数ｍが多ければ多いほど放出地点の推定精度は向上し、そのためには仮想格子の離間距離をより小さくして交差点数を増やせば良いが、その分だけ計算量は増大することになる。したがって、要求される処理時間（観測地点から情報を入手してから放出地点を推定するまでの時間）に応じて、予め発生源推定処理部１５を具現するプロセッサの処理性能等を勘案して仮想放出地点数ｍを概算しておくのが望ましい。また、経験則等に基づいて注目すべき地域を判断し、その注目地域に仮想格子の領域を設定することが望ましい。

また、仮想放出時刻設定部２３は、放出時刻が不明であるとき仮想放出時刻を設定する。つまり、放出時刻が不明の場合には、仮想放出時刻を所定時間刻みで複数設定して式（７）を評価することとしている。

また、影響関数算出部２４は、拡散モデルを用いて影響関数を算出する。上述のように、影響関数Ｄ_ｉｊは評価地点（観測器５−ｉ（ｉ＝１〜ｎ）の位置）と仮想放出地点Ｐｏｊ（ｊ＝１〜ｍ）との相対位置、並びに（放出時刻と観測器５−ｉの計測時刻との）相対時刻に応じて定まる関数であり、ｎ×ｍ個の影響関数Ｄ_ｉｊを算出することになる。

本実施形態では、平地一様流場（拡散域の地形が平地であり、風が一様に流れている状態）を想定する場合には、拡散モデルとしてパフモデルを使用する。パフモデルによれば、風速をＵ［m/sec］としたとき、拡散係数Ｄ_ｉｊは次式で与えられる。

ここに、σ_ｘ，σ_ｙ，σ_ｚはそれぞれ濃度分布のｘ，ｙ，ｚ方向の拡散パラメータ［m］であり、パスキル・ギフォード線図や経験式等に基づいて得られる。なお、拡散モデルとしては、パフモデルに限定されることなく、例えばプルームモデル等の他のモデルを使用することも可能である。

また、一般の市街地などでは、建屋、地形等の影響で一様流でないケースが多いが、このような複雑気流場の場合には、数値拡散計算により影響関数を算出する。すなわち、各種シミュレーションモデルを用いて、想定する放出地点から単位強度の放出をした場合の評価地点での濃度（即ち、影響係数）を求めるものである。このシミュレーションモデルとしては、例えば、修正プルームモデル、ポテンシャル流モデル、粘性流モデルが挙げられる。なお、これらの詳細については、例えば、三菱重工技報vol.21 No.5抜刷（1984年9月）「排煙拡散数値シミュレーションモデルの開発」pp.1-8 にて、開示されている。また、環境庁大気保全局大気規制課編「窒素酸化物総量規制マニュアル」にも詳述されているプルーム・パフモデルや、同じく周知のモデルであるセル内粒子法、ラグランジュ型粒子モデルを使用しても求めることができる。なお、この数値拡散計算は、仮想放出地点数ｍが例えば２５（＝５×５）点程度であれば、計算時間の負荷は小さく全体の計算量に与える影響は少ない。

なお、放出時刻が不明で、仮想放出時刻設定部２３によって複数（ｒ個）の仮想放出時刻が設定されている場合には、ｎ×ｍ×ｒ個（ｒは正整数）の影響関数Ｄ_ｉｊを算出する必要がある。したがって、複雑気流場の場合、仮想放出地点数ｍおよび仮想放出時刻数ｒの値によっては、影響関数算出処理の計算量が全体に与える影響が大きくなってくることも考えられる。そのような場合には、観測情報を入手してから逐一算出処理を行うのではなく、想定され得る地域毎の仮想格子（仮想放出地点Ｐｏｊ）について予め相対時刻毎に影響関数Ｄ_ｉｊを算出し、影響関数データベースに登録しておくようにして良い。ここで、影響関数データベースのデータは記憶部１７内に保持される。これにより、影響関数算出処理の全体計算量に与える影響を殆ど無くすことができる。

また、残差ノルム算出部２５は、仮想放出地点Ｐｏｊ毎に、各観測器５−ｉ（ｉ＝１〜ｎ）の濃度情報と、各観測器５−ｉの該仮想放出地点Ｐｏｊに対する影響関数Ｄ_ｉｊと該仮想放出地点の放出強度ｑ_ｊとの積と、の差の平方和である残差ノルムＲ_ｊを算出する。すなわち、仮想格子設定部２２によって設定された各仮想放出地点Ｐｏｊについて、式（６）に基づき残差ノルムＲ_ｊを算出する。なお、放出時刻が不明で、仮想放出時刻設定部２３によって仮想放出時刻が設定されている場合には、設定されている複数の仮想放出時刻に対してそれぞれ残差ノルムを算出することとなる。

また、推定部２６は、算出された全仮想放出地点Ｐｏｊについての残差ノルムの内、残差ノルムが最小となる放出強度ｑ_ｊを求め、該放出強度ｑ_ｊの仮想放出地点Ｐｏｊを放出位置とし、また該放出強度ｑ_ｊを該放出位置からの放出量として推定する。なお、放出時刻が不明で、仮想放出時刻設定部２３によって仮想放出時刻が設定されている場合には、複数の仮想放出時刻に対してそれぞれ最小となる残差ノルムの中でさらに最小の残差ノルムの放出強度ｑ_ｊを求め、該放出強度ｑ_ｊの仮想放出地点を放出位置とし、また該放出強度ｑ_ｊを該放出位置からの放出量として推定すると共に、対応する仮想放出時刻を放出時刻として推定する。

次に、以上のような構成要素を備えた発生源推定装置３における拡散物質の発生源推定方法について、図５を参照して説明する。ここで、図５は本実施形態に係る拡散物質の発生源推定方法を説明するフローチャートである。

まず、ステップＳ１０１では、観測情報入手部２１において、観測器５−ｉ（ｉ＝１〜ｎ）の位置情報、該観測器５−ｉが計測した濃度情報および計測時刻情報を、通信インタフェース１１を介して入手する。

次に、ステップＳ１０２では、仮想格子設定部２２により、相互の離間距離が一定の仮想格子を想定して、該仮想格子の領域上で各格子線が交差する位置を仮想放出地点Ｐｏｊ（ｊ＝１〜ｍ）として設定する。なお、ここで、放出時刻が不明であるときには、仮想放出時刻設定部２３により仮想放出時刻を設定する。

次に、ステップＳ１０３では、影響関数算出部２４により、拡散モデル（例えば式（８））を用いて影響関数Ｄ_ｉｊを算出する。なお、仮想放出時刻設定部２３により仮想放出時刻が設定されている場合には、各仮想放出時刻に応じた相対時刻毎に影響関数Ｄ_ｉｊを算出することとなる。また、上述したように予め影響関数データベースに登録されている場合には、影響関数データベースを参照して影響関数Ｄ_ｉｊを獲得するようにしても良い。

次に、ステップＳ１０４では、残差ノルム算出部２５により、仮想格子設定部２２で設定された各仮想放出地点Ｐｏｊについて、式（６）に基づき残差ノルムＲ_ｊを算出する。なお、仮想放出時刻設定部２３により仮想放出時刻が設定されている場合には、設定されている複数の仮想放出時刻に対してそれぞれ残差ノルムＲ_ｊを算出する。

また、ステップＳ１０４では、推定部２６により、算出された全仮想放出地点Ｐｏｊについての残差ノルムの内、残差ノルムが最小となる放出強度ｑ_ｊを求め、該放出強度ｑ_ｊの仮想放出地点Ｐｏｊを放出位置とし、また該放出強度ｑ_ｊを該放出位置からの放出量として推定する。なお、仮想放出時刻設定部２３によって仮想放出時刻が設定されている場合には、複数の仮想放出時刻に対してそれぞれ最小となる残差ノルムの中でさらに最小の残差ノルムの放出強度ｑ_ｊを求め、該放出強度ｑ_ｊの仮想放出地点を放出位置とし、また該放出強度ｑ_ｊを該放出位置からの放出量として推定すると共に、対応する仮想放出時刻を放出時刻として推定する。

以上説明したように、本実施形態の拡散物質の発生源推定装置および発生源推定方法では、観測器５−ｉ（ｉ＝１〜ｎ）の位置情報および計測した濃度情報を観測情報入手部２１によって入手し、仮想格子設定部２２により相互の離間距離が一定の仮想格子上で各格子線が交差する位置を仮想放出地点Ｐｏｊ（ｊ＝１〜ｍ）として設定し、影響関数算出部２４により、拡散モデルを用いて観測器５−ｉと仮想放出地点Ｐｏｊとの相対位置、並びに相対時刻に応じて定まる影響関数Ｄ_ｉｊを算出し、残差ノルム算出部２５により、仮想放出地点Ｐｏｊ毎に、各観測器５−ｉの濃度情報と、各観測器５−ｉの該仮想放出地点Ｐｏｊに対する影響関数Ｄ_ｉｊと該仮想放出地点Ｐｏｊの放出強度ｑ_ｊとの積と、の差の平方和である残差ノルムＲ_ｊを算出し、推定部２６により、算出された全仮想放出地点の残差ノルムの内、残差ノルムが最小となる仮想放出地点を放出地点と推定する。

このように、設定した仮想放出地点Ｐｏｊ毎の残差ノルムＲ（ｑ_ｊ）を評価して、該残差ノルムＲ（ｑ_ｊ）が最小となる放出強度を求め、該放出強度の仮想放出地点を放出位置とし、また該放出強度を該放出位置からの放出量として推定するので、従来のように「観測地点の数ｎ≧仮想放出地点の数ｍ」という制約を持つことなく放出地点を推定することができ、より柔軟且つ簡便に発生源を推定し得る拡散物質の発生源推定装置および発生源推定方法を実現することができる。

また、本実施形態の拡散物質の発生源推定装置および発生源推定方法では、平地一様流場ではパフモデル等の拡散モデルを用いて影響関数を算出し、また複雑気流場では数値拡散（シミュレーション）計算により想定する放出地点から単位強度の放出をした場合の評価地点での濃度として影響関数を算出するので、様々な地形における拡散に対してより精度良く発生源を推定することができる。

また、本実施形態の拡散物質の発生源推定装置および発生源推定方法では、放出時刻が不明な場合には、仮想放出時刻設定部２３によって仮想放出時刻を設定することとしたので、より柔軟に発生源を推定し得る拡散物質の発生源推定装置および発生源推定方法を実現することができる。

さらに、本実施形態の拡散物質の発生源推定装置および発生源推定方法では、影響関数算出部２４において影響関数データベースを参照して予め算出した影響関数を獲得し、その後の処理に用いる構成（手順）とすることにより、影響関数を算出する処理分の計算量を排除することができ、全体的な処理の計算量を低減してより速く発生源の推定を行うことができる。

〔第２実施形態〕
次に、本発明の第２実施形態に係る拡散物質の発生源推定装置および発生源推定方法について説明する。本実施形態における拡散物質の発生源推定装置の構成は第１実施形態の構成（図１参照）と同等であり、各構成要素についての詳細な説明は省略する。

但し、第１実施形態では１面の仮想格子上で仮想放出地点Ｐｏｊ（ｊ＝１〜ｍ）を設定したのに対し、本実施形態では、格子線の相互の離間距離が段階的に異なるｓ面（ｓは正整数）の仮想格子（一面の仮想格子における格子線の相互の離間距離は一定）上で、それぞれ仮想放出地点ＰｏＮｊ（Ｎ＝１〜ｓ，ｊ＝１〜ｐ；ｐはｐ＜ｍで正整数）を設定する点が異なる。なお、仮想格子の面の適用は、格子線の相互の離間距離が相対的に長い仮想格子から相対的に短い仮想格子へ段階的に使用していくものとし、ｓ段目の仮想格子の面が（離間距離が最短の）最小グリッドであるものとする。

また、本実施形態では、想定され得る地域毎の仮想格子（仮想放出地点ＰｏＮｊ）について、予め影響関数算出部２４によって相対時刻毎に影響関数Ｄ_ｉｊが算出され、影響関数データベース１８に登録されて、データが記憶部１７内に保持されているものとする。

次に、本実施形態における拡散物質の発生源推定方法について、図６〜図８を参照して説明する。ここで、図６は本実施形態に係る拡散物質の発生源推定方法を説明するフローチャートであり、図７はＮ段目の仮想格子を用いて残差ノルムを算出する例を例示する説明図であり、図８はＮ＋１段目の仮想格子を用いて残差ノルムを算出する例を例示する説明図である。

まず、ステップＳ２０１では、観測情報入手部２１において、観測器５−ｉ（ｉ＝１〜ｎ）の位置情報、該観測器５−ｉが計測した濃度情報および計測時刻情報を、通信インタフェース１１を介して入手する。

次に、ステップＳ２０２では、仮想格子設定部２２により、Ｎ段目の仮想格子を想定して、該仮想格子の領域上で各格子線が交差する位置を仮想放出地点ＰｏＮｊ（Ｎ＝１〜ｓ，ｊ＝１〜ｐ）として設定する。なお、ここで、放出時刻が不明であるときには、仮想放出時刻設定部２３により仮想放出時刻を設定する。

次に、ステップＳ２０３では、影響関数算出部２４により、影響関数データベース１８を参照して影響関数Ｄ_ｉｊを獲得する。

次に、ステップＳ２０４では、残差ノルム算出部２５により、仮想格子設定部２２で設定された各仮想放出地点ＰｏＮｊについて、式（６）に基づき残差ノルムＲ_ｊを算出する。なお、仮想放出時刻設定部２３により仮想放出時刻が設定されている場合には、設定されている複数の仮想放出時刻に対してそれぞれ残差ノルムＲ_ｊを算出する。

また、ステップＳ２０４では、推定部２６により、算出された全仮想放出地点ＰｏＮｊについての残差ノルムの内、残差ノルムが最小となる放出強度ｑ_ｊを求め、該放出強度ｑ_ｊの仮想放出地点ＰｏＮｊを放出候補地点の位置とし、また該放出強度ｑ_ｊを該放出候補地点からの放出量として推定する。なお、仮想放出時刻設定部２３によって仮想放出時刻が設定されている場合には、複数の仮想放出時刻に対してそれぞれ最小となる残差ノルムの中でさらに最小の残差ノルムの放出強度ｑ_ｊを求め、該放出強度ｑ_ｊの仮想放出地点を放出候補地点の位置とし、また該放出強度ｑ_ｊを該放出候補地点からの放出量として推定すると共に、対応する仮想放出時刻を放出時刻として推定する。

次に、ステップＳ２０５では、現在適用している仮想格子の面が最小グリッド（ｓ段目で）であるか否かを判断し、最小グリッドである場合には終了し、また最小グリッドでなければ、ステップＳ２０６に進んでＮをインクリメントした後、ステップＳ２０２に戻る。

ここで、図７および図８に示すように、仮想格子の面として３×３の仮想格子を使用し、観測器５−ｉによる観測地点が１０箇所（Ｐｖ１〜Ｐｖ１０）あるとした場合について説明する。まず、Ｎ段目（Ｎ＝１）の仮想格子を用いて残差ノルムを算出する場合には、図７に示す如く、仮想格子設定部２２により９個の仮想放出地点Ｐｏ１１〜Ｐｏ１９が設定され（ステップＳ２０２）、推定部２６により、全仮想放出地点Ｐｏ１１〜Ｐｏ１９についての残差ノルムを評価したところ、放出強度ｑ_４で残差ノルムが最小となり、該放出強度ｑ_４の仮想放出地点Ｐｏ１４が放出候補地点となった。

次に、Ｎがインクリメントされて（ステップＳ２０６）、Ｎ＋１段目（Ｎ＋１＝２）の仮想格子を用いて残差ノルムを算出する場合には、図８に示す如く、仮想格子設定部２２により、９個の仮想放出地点Ｐｏ２１〜Ｐｏ２９が設定され（ステップＳ２０２）、推定部２６により、全仮想放出地点Ｐｏ２１〜Ｐｏ２９についての残差ノルムを評価したところ、放出強度ｑ_２で残差ノルムが最小となり、該放出強度ｑ_２の仮想放出地点Ｐｏ２２が放出候補地点となった。なお、Ｎ＋１段目の仮想格子は、前段（Ｎ段目）の仮想格子上で放出候補地点とされたＰｏ１４を、その仮想格子面内に含むように設定される。図８では、放出候補地点Ｐｏ１４が面の中心となるように、Ｎ＋１段目の仮想格子上で仮想放出地点Ｐｏ２５が重複して設定されているが、必ずしも重複設定する必要はなく、仮想格子面内に含んでいれば良い。

このように、格子線の相互の離間距離が相対的に長い仮想格子から相対的に短い仮想格子へと、段階的に、（離間距離が最短で）最小グリッドとなるｓ段目の仮想格子の面まで適用していき、ｓ段目の仮想格子の面の適用において、推定部２６により、全仮想放出地点ＰｏＮｊについての残差ノルムの内、残差ノルムが最小となる放出強度ｑ_ｊが求められ、該放出強度ｑ_ｊの仮想放出地点ＰｏＮｊが放出候補地点と推定されたとき、該放出候補地点が放出地点として推定されることとなる。

以上説明したように、本実施形態の拡散物質の発生源推定装置および発生源推定方法では、仮想格子設定部２２において、前段の仮想格子面の適用で推定部２６によって推定された放出地点を含み、且つ相互の離間距離がより短い仮想格子上で各格子線が交差する位置を仮想放出地点として再設定することとし、格子線の相互の離間距離が相対的に長い仮想格子から相対的に短い仮想格子へと絞り込みながら仮想放出地点を再設定して発生源を推定するようにしている。これにより、第１実施形態（離間距離を最短とした１面の仮想格子上で仮想放出地点を設定して発生源を推定する場合）と比べて、１面の仮想格子上での仮想放出地点数を格段に小さくすることができ、全体的な処理の計算量を抑制して計算時間を短縮し、より速く発生源の推定を行うことができる。

以上、本発明の実施形態について図面を参照して詳述したが、本発明はこれら実施形態に限定されるものではなく、本発明の要旨を逸脱しない範囲の設計変更等があっても本発明に含まれる。

３ 拡散物質の発生源推定装置
５−１〜５−ｎ 観測器
１１ 通信インタフェース
１３ 入力部
１５ 発生源推定処理部
１７ 記憶部
１８ 影響関数データベース
１９ 出力部
２１ 観測情報入手部
２２ 仮想格子設定部
２３ 仮想放出時刻設定部
２４ 影響関数算出部
２５ 残差ノルム算出部
２６ 推定部

Claims (10)

1. 複数個の観測器からの情報に基づき、ガス発生源情報を推定する拡散物質の発生源推定装置であって、
   前記観測器の位置情報および計測した濃度情報を入手する観測情報入手手段と、
   相互の離間距離が一定の仮想格子上で各格子線が交差する位置を仮想放出地点として設定する仮想格子設定手段と、
   拡散モデルを用いて、前記観測器と前記仮想放出地点との相対位置、並びに相対時刻に応じて定まる影響関数を算出する影響関数算出手段と、
   前記仮想放出地点毎に、各観測器の濃度情報と、各観測器の該仮想放出地点に対する影響関数と該仮想放出地点の放出強度との積と、の差の平方和である残差ノルムを算出する残差ノルム算出手段と、
   算出された全仮想放出地点の残差ノルムの内、残差ノルムが最小となる仮想放出地点を放出地点と推定する推定手段と、
   を有することを特徴とする拡散物質の発生源推定装置。
2. 前記影響関数算出手段は、数値拡散計算により前記影響関数を算出することを特徴とする請求項１に記載の拡散物質の発生源推定装置。
3. 前記仮想格子設定手段は、前記推定手段で推定された放出地点を含み、且つ相互の離間距離がより短い仮想格子上で各格子線が交差する位置を仮想放出地点として再設定することを特徴とする請求項１または請求項２に記載の拡散物質の発生源推定装置。
4. 仮想放出時刻を設定する仮想放出時刻設定手段、を有し、
   前記残差ノルム算出手段は、仮想放出時刻別に各仮想放出地点毎の残差ノルムを算出し、
   前記推定手段は、仮想放出時刻別に算出された全仮想放出地点の残差ノルムの内、残差ノルムが最小となる仮想放出時刻および仮想放出地点をそれぞれ放出時刻および放出地点と推定することを特徴とする請求項１〜請求項３の何れか１項に記載の拡散物質の発生源推定装置。
5. 前記影響関数算出手段は、想定される前記観測器と前記仮想放出地点との相対位置およびまたは相対時刻に応じた影響関数を予め算出してデータベースに保存しておくことを特徴とする請求項１〜請求項４の何れか１項に記載の拡散物質の発生源推定装置。
6. 複数個の観測器からの情報に基づき、ガス発生源情報を推定する拡散物質の発生源推定方法であって、
   前記観測器の位置情報および計測した濃度情報を入手する観測情報入手ステップと、
   相互の離間距離が一定の仮想格子上で各格子線が交差する位置を仮想放出地点として設定する仮想格子設定ステップと、
   拡散モデルを用いて、前記観測器と前記仮想放出地点との相対位置、並びに相対時刻に応じて定まる影響関数を算出する影響関数算出ステップと、
   前記仮想放出地点毎に、各観測器の濃度情報と、各観測器の該仮想放出地点に対する影響関数と該仮想放出地点の放出強度との積と、の差の平方和である残差ノルムを算出する残差ノルム算出ステップと、
   算出された全仮想放出地点の残差ノルムの内、残差ノルムが最小となる仮想放出地点を放出地点と推定する推定ステップと、
   を有することを特徴とする拡散物質の発生源推定方法。
7. 前記影響関数算出ステップは、数値拡散計算により前記影響関数を算出することを特徴とする請求項６に記載の拡散物質の発生源推定方法。
8. 前記仮想格子設定ステップは、前記推定ステップで推定された放出地点を含み、且つ相互の離間距離がより短い仮想格子上で各格子線が交差する位置を仮想放出地点として再設定することを特徴とする請求項６または請求項７に記載の拡散物質の発生源推定方法。
9. 仮想放出時刻を設定する仮想放出時刻設定ステップ、を有し、
   前記残差ノルム算出ステップは、仮想放出時刻別に各仮想放出地点毎の残差ノルムを算出し、
   前記推定ステップは、仮想放出時刻別に算出された全仮想放出地点の残差ノルムの内、残差ノルムが最小となる仮想放出時刻および仮想放出地点をそれぞれ放出時刻および放出地点と推定することを特徴とする請求項６〜請求項８の何れか１項に記載の拡散物質の発生源推定方法。
10. 前記影響関数算出ステップは、想定される前記観測器と前記仮想放出地点との相対位置およびまたは相対時刻に応じた影響関数を予め算出してデータベースに保存しておくことを特徴とする請求項６〜請求項９の何れか１項に記載の拡散物質の発生源推定方法。

Images