JP2002207002A - 大気環境シミュレーションシステム - Google Patents
大気環境シミュレーションシステムInfo
- Publication number
- JP2002207002A JP2002207002A JP2001002458A JP2001002458A JP2002207002A JP 2002207002 A JP2002207002 A JP 2002207002A JP 2001002458 A JP2001002458 A JP 2001002458A JP 2001002458 A JP2001002458 A JP 2001002458A JP 2002207002 A JP2002207002 A JP 2002207002A
- Authority
- JP
- Japan
- Prior art keywords
- simulation
- data
- source
- numerical
- input
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Instructional Devices (AREA)
- Processing Or Creating Images (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
(57)【要約】
【課題】大気汚染物質の測定値（実測データ）と、汚染
物質の拡散分布状況を数値シミュレーションした計算結
果とを組合わせて汚染物質発生源位置を推定することが
できる大気環境シミュレーションシステムを提供するこ
と。 【解決手段】数値計算処理部１５は選定されたシミュレ
ーションモデル、固定パラメータ、解析条件および調整
パラメータに基づいて得られた数値計算データを用いて
汚染物質の拡散分布状況のシミュレーションを実行す
る。その結果と地図を表示する一連の解析を発生源の個
数だけ繰り返し実施し、それらの結果と大気汚染物質の
測定値とを統計的手法により大気汚染物質の発生源位置
を発生源推定手段１６で推定する。
物質の拡散分布状況を数値シミュレーションした計算結
果とを組合わせて汚染物質発生源位置を推定することが
できる大気環境シミュレーションシステムを提供するこ
と。 【解決手段】数値計算処理部１５は選定されたシミュレ
ーションモデル、固定パラメータ、解析条件および調整
パラメータに基づいて得られた数値計算データを用いて
汚染物質の拡散分布状況のシミュレーションを実行す
る。その結果と地図を表示する一連の解析を発生源の個
数だけ繰り返し実施し、それらの結果と大気汚染物質の
測定値とを統計的手法により大気汚染物質の発生源位置
を発生源推定手段１６で推定する。
Description
【０００１】
【発明の属する技術分野】本発明は大気汚染物質の発生
源による大気汚染物質の拡散分布状態を解析した結果と
大気汚染濃度の実測値（測定データ）とを組合わせて、
統計的手法を用いたデータ解析を行うことにより、実測
データに寄与している発生源の位置等を推定し、対策検
討のための情報を提供する大気環境シミュレーションシ
ステムに関する。
源による大気汚染物質の拡散分布状態を解析した結果と
大気汚染濃度の実測値（測定データ）とを組合わせて、
統計的手法を用いたデータ解析を行うことにより、実測
データに寄与している発生源の位置等を推定し、対策検
討のための情報を提供する大気環境シミュレーションシ
ステムに関する。
【０００２】
【従来の技術】近年、地球環境保全の観点から環境アセ
スメント法の制定やISO14000シリーズの発行など、大気
環境の改善や大気汚染対策のため環境影響評価が重要に
なってきている。特に大気環境に汚染物質を放出する工
場や発電所、燃焼施設では計画・立案段階は勿論のこ
と、運転稼働中も環境影響を定期的に把握する必要があ
る。
スメント法の制定やISO14000シリーズの発行など、大気
環境の改善や大気汚染対策のため環境影響評価が重要に
なってきている。特に大気環境に汚染物質を放出する工
場や発電所、燃焼施設では計画・立案段階は勿論のこ
と、運転稼働中も環境影響を定期的に把握する必要があ
る。
【０００３】従来、大気汚染物質を放出する工場や発電
所、燃焼施設等からの汚染物質の拡散予測を行い、予測
結果を表示装置上に表示するシステムとしては、例えば
特開平１１−１４５３１号公報が挙げられる。このよう
な予測システムは発生源からの汚染物質の拡散分布を計
算するものである。
所、燃焼施設等からの汚染物質の拡散予測を行い、予測
結果を表示装置上に表示するシステムとしては、例えば
特開平１１−１４５３１号公報が挙げられる。このよう
な予測システムは発生源からの汚染物質の拡散分布を計
算するものである。
【０００４】
【発明が解決しようとする課題】しかしながら、例えば
実測された大気汚染物質が複数の発生源のどれから放出
されたものか、すなわち汚染物質濃度の測定値（実測デ
ータ）から汚染物質発生源の位置等を推定することは困
難であった。
実測された大気汚染物質が複数の発生源のどれから放出
されたものか、すなわち汚染物質濃度の測定値（実測デ
ータ）から汚染物質発生源の位置等を推定することは困
難であった。
【０００５】本発明の目的は、大気汚染物質の測定値
（実測データ）と、汚染物質の拡散分布状況を数値シミ
ュレーションした計算結果とを組合わせて汚染物質発生
源の位置を推定することができる大気環境シミュレーシ
ョンシステムを提供することにある。
（実測データ）と、汚染物質の拡散分布状況を数値シミ
ュレーションした計算結果とを組合わせて汚染物質発生
源の位置を推定することができる大気環境シミュレーシ
ョンシステムを提供することにある。
【０００６】
【課題を解決するための手段】本発明は上記目的を達成
するために、大気汚染物質の発生源位置を表示画面上に
格子状に規則正しく複数個設定し、各発生源から放出さ
れる汚染物質の拡散分布状況をシミュレーションするシ
ミュレーション条件データに基づきこのシミュレーショ
ンに適したシミュレーションモデル、固定パラメータ及
び解析条件を選定するとともに、シミュレーション条件
データにより拡散分布状況の数値シミュレーションを行
い類似している環境における実測値（測定値）との誤差
が小さくなるように調整パラメータを選定し、選定され
たシミュレーションモデル、固定パラメータ、解析条件
および調整パラメータに基づいて得られた数値計算デー
タを用いて汚染物質の拡散分布状況の数値シミュレーシ
ョンを実行してその結果と地図を表示する一連の解析を
発生源の個数だけ繰り返し実施し、それらの結果と大気
汚染物質の測定値とを統計的手法を用いて計算すること
により大気汚染物質の発生源位置を推定して表示装置に
表示するようにしたことを特徴とする。
するために、大気汚染物質の発生源位置を表示画面上に
格子状に規則正しく複数個設定し、各発生源から放出さ
れる汚染物質の拡散分布状況をシミュレーションするシ
ミュレーション条件データに基づきこのシミュレーショ
ンに適したシミュレーションモデル、固定パラメータ及
び解析条件を選定するとともに、シミュレーション条件
データにより拡散分布状況の数値シミュレーションを行
い類似している環境における実測値（測定値）との誤差
が小さくなるように調整パラメータを選定し、選定され
たシミュレーションモデル、固定パラメータ、解析条件
および調整パラメータに基づいて得られた数値計算デー
タを用いて汚染物質の拡散分布状況の数値シミュレーシ
ョンを実行してその結果と地図を表示する一連の解析を
発生源の個数だけ繰り返し実施し、それらの結果と大気
汚染物質の測定値とを統計的手法を用いて計算すること
により大気汚染物質の発生源位置を推定して表示装置に
表示するようにしたことを特徴とする。
【０００７】
【発明の実施の形態】以下、本発明の実施の形態を、図
面を参照して説明する。
面を参照して説明する。
【０００８】図１に本発明の一実施例に係る大気環境シ
ミュレーションシステムの構成図を示す。
ミュレーションシステムの構成図を示す。
【０００９】図1において、大気環境シミュレーション
システム１０は、入力データ（シミュレーション条件デ
ータ）を入力する入力装置１と解析結果を表示，出力す
るディスプレイやプリンタ等の出力装置２とを結合する
もので、モデル選定部１１、データベース部１２、入力
パラメータチューニング部１３、入力作成部１４、数値
計算処理部１５、発生源推定処理部１６、表示合成部１
７、データ格納部１８、可視化データ選定処理部１９、
表示用処理部２０から構成されている。
システム１０は、入力データ（シミュレーション条件デ
ータ）を入力する入力装置１と解析結果を表示，出力す
るディスプレイやプリンタ等の出力装置２とを結合する
もので、モデル選定部１１、データベース部１２、入力
パラメータチューニング部１３、入力作成部１４、数値
計算処理部１５、発生源推定処理部１６、表示合成部１
７、データ格納部１８、可視化データ選定処理部１９、
表示用処理部２０から構成されている。
【００１０】モデル選定部１１には入力装置１から汚染
物質の種類や放出量、放出位置、大気条件などの入力条
件３が入力される。モデル選定部１１は入力された入力
条件３に基づき汚染物質の種類や形状，放出位置，放出
条件，評価時点での大気条件，評価地点の地形状況、等
の観点から分類を行い、データベース部１２の中に格納
されている解析事例データベース１２１、法令・指針デ
ータベース１２２およびシミュレーションモデルデータ
ベース１２３を参照することによって、入力条件選定部
（Ａ）４において入力条件３に最適なシミュレーション
モデル４１と固定パラメータ４２と解析条件４３を選定
する。
物質の種類や放出量、放出位置、大気条件などの入力条
件３が入力される。モデル選定部１１は入力された入力
条件３に基づき汚染物質の種類や形状，放出位置，放出
条件，評価時点での大気条件，評価地点の地形状況、等
の観点から分類を行い、データベース部１２の中に格納
されている解析事例データベース１２１、法令・指針デ
ータベース１２２およびシミュレーションモデルデータ
ベース１２３を参照することによって、入力条件選定部
（Ａ）４において入力条件３に最適なシミュレーション
モデル４１と固定パラメータ４２と解析条件４３を選定
する。
【００１１】入力条件選定部（Ｂ）５は入力条件３に基
づく解析結果の精度を向上させ、実測値（測定値）と近
づけるため、入力パラメータチューニング部１３でデー
タベース部１２の中に格納されている測定値データベー
ス１２４を参照して調整パラメータ５１を選定し、入力
作成部１４に与える。入力作成部１４は選定した条件を
基に、数値計算処理部１５のシミュレーション計算に必
要な数値計算入力データ６を作成する。
づく解析結果の精度を向上させ、実測値（測定値）と近
づけるため、入力パラメータチューニング部１３でデー
タベース部１２の中に格納されている測定値データベー
ス１２４を参照して調整パラメータ５１を選定し、入力
作成部１４に与える。入力作成部１４は選定した条件を
基に、数値計算処理部１５のシミュレーション計算に必
要な数値計算入力データ６を作成する。
【００１２】数値計算処理部１５は、入力作成部１４で
作成された数値計算入力データを用い、内蔵した物理モ
デルから解析対象領域全範囲にわたる物理現象を並列計
算等の高速化手法を用いて数値解析を実施し、その対象
領域全範囲の高精度数値解析結果１５１を表示合成部１
７のデータ格納部１８にある解析結果データベース１８
２の中に分類整理して格納する。
作成された数値計算入力データを用い、内蔵した物理モ
デルから解析対象領域全範囲にわたる物理現象を並列計
算等の高速化手法を用いて数値解析を実施し、その対象
領域全範囲の高精度数値解析結果１５１を表示合成部１
７のデータ格納部１８にある解析結果データベース１８
２の中に分類整理して格納する。
【００１３】表示合成部１７は、可視化に必要なデータ
を格納しておくデータ格納部１８と可視化に必要なデー
タを選定する可視化データ選定処理部１９と表示用処理
部２０から構成される。
を格納しておくデータ格納部１８と可視化に必要なデー
タを選定する可視化データ選定処理部１９と表示用処理
部２０から構成される。
【００１４】データ格納部１８には解析結果を格納する
解析結果データベース１８２の他に解析結果表示の際に
背景データとして使用する地図データを格納した、地図
データベース１８１を有する。地図データベース１８１
は解析対象領域全範囲に関する多種類の精度の地図デー
タが格納されている。
解析結果データベース１８２の他に解析結果表示の際に
背景データとして使用する地図データを格納した、地図
データベース１８１を有する。地図データベース１８１
は解析対象領域全範囲に関する多種類の精度の地図デー
タが格納されている。
【００１５】可視化データ選定処理部１９は入力装置１
から入力された表示条件７に基づきデータ格納部１８の
データベース１８１、１８２を参照し、表示用処理部２
０へ表示条件に適した表示精度地図データ２０１と表示
精度可視化用物理量分布データ２０２を選定し指示す
る。表示用処理部２０は可視化データ選定処理部１９で
指示された表示精度地図データ２０１と表示精度可視化
用物理量分布データ２０２を可視化処理を行い、合成し
て出力装置２に描画する物理量分布可視化データ２０３
を作成する。
から入力された表示条件７に基づきデータ格納部１８の
データベース１８１、１８２を参照し、表示用処理部２
０へ表示条件に適した表示精度地図データ２０１と表示
精度可視化用物理量分布データ２０２を選定し指示す
る。表示用処理部２０は可視化データ選定処理部１９で
指示された表示精度地図データ２０１と表示精度可視化
用物理量分布データ２０２を可視化処理を行い、合成し
て出力装置２に描画する物理量分布可視化データ２０３
を作成する。
【００１６】これまでの一連の操作を、入力条件３にお
いて汚染物質の発生源位置のみを格子状に規則正しく複
数個設定し、各発生源条件ごとに数値計算処理部１５１
において複数回繰り返し計算させ、複数回の計算結果を
データ格納部１８の解析結果データベース１８２に格納
する。汚染物質発生源位置が異なる前記複数の計算結果
および入力条件３であらかじめ入力された実測値（測定
値）に基づき、発生源推定処理部１６において統計的手
法を用いて計算処理を行い、汚染物質の発生源として最
も確からしい位置を推定し、推定結果を再び解析結果デ
ータベース１８２に格納する。
いて汚染物質の発生源位置のみを格子状に規則正しく複
数個設定し、各発生源条件ごとに数値計算処理部１５１
において複数回繰り返し計算させ、複数回の計算結果を
データ格納部１８の解析結果データベース１８２に格納
する。汚染物質発生源位置が異なる前記複数の計算結果
および入力条件３であらかじめ入力された実測値（測定
値）に基づき、発生源推定処理部１６において統計的手
法を用いて計算処理を行い、汚染物質の発生源として最
も確からしい位置を推定し、推定結果を再び解析結果デ
ータベース１８２に格納する。
【００１７】このように、入力装置１から格子状に規則
正しく配置した汚染物質発生源位置、汚染物質濃度が実
測されている位置及び濃度、並びにその他の入力条件を
入力条件３として入力することにより、複数の発生源位
置に対する汚染物質の拡散分布状況を数値計算処理部１
５において計算し、複数の計算結果から汚染物質濃度が
実測されている位置における計算値を算出し、当該位置
における汚染物質濃度の実測値（測定データ）と比較す
ることにより、どの位置の汚染物質発生源から放出され
た可能性が高いのかを容易に推定することができる。
正しく配置した汚染物質発生源位置、汚染物質濃度が実
測されている位置及び濃度、並びにその他の入力条件を
入力条件３として入力することにより、複数の発生源位
置に対する汚染物質の拡散分布状況を数値計算処理部１
５において計算し、複数の計算結果から汚染物質濃度が
実測されている位置における計算値を算出し、当該位置
における汚染物質濃度の実測値（測定データ）と比較す
ることにより、どの位置の汚染物質発生源から放出され
た可能性が高いのかを容易に推定することができる。
【００１８】さらに、表示条件７を可視化データ選定処
理部１９で表示範囲に応じた表示精度可視化用物理量分
布データ２０２と表示精度地図データ２０１を、また発
生源位置座標データ２０４と表示精度地図データ２０１
を、をそれぞれ合成して表示することができる。
理部１９で表示範囲に応じた表示精度可視化用物理量分
布データ２０２と表示精度地図データ２０１を、また発
生源位置座標データ２０４と表示精度地図データ２０１
を、をそれぞれ合成して表示することができる。
【００１９】次にデータベース部に格納されるデータに
ついて説明する。
ついて説明する。
【００２０】図２に法令・指針データベース１２２に格
納されているデータの一例を示す。
納されているデータの一例を示す。
【００２１】図２（ａ）は法令・指針データベース１２
２の例を示す。法令・指針データベース１２２には、法
令・指針名１２２１と環境項目１２２２、解析対象分野
１２２３、地域１２２４、モデル１２２５に分類された
情報があり、これらの関係を示す対応表１２２６から構
成される。対応表１２２６には、それぞれの情報の関係
だけでなくその重要性も数値化して持っている。この事
例では重要性を二重丸、丸、三角の記号で表している。
２の例を示す。法令・指針データベース１２２には、法
令・指針名１２２１と環境項目１２２２、解析対象分野
１２２３、地域１２２４、モデル１２２５に分類された
情報があり、これらの関係を示す対応表１２２６から構
成される。対応表１２２６には、それぞれの情報の関係
だけでなくその重要性も数値化して持っている。この事
例では重要性を二重丸、丸、三角の記号で表している。
【００２２】図２（ｂ）は解析事例データベース１２１
の一例を示す。解析事例データベース１２１には、排出
源情報１２１１、モデル１２１２、拡散パラメータ１２
１３、有効煙突高さ１２１４などの解析に用いた各種条
件を分類した情報と、これらの関係を示す対応表１２１
５から構成される。
の一例を示す。解析事例データベース１２１には、排出
源情報１２１１、モデル１２１２、拡散パラメータ１２
１３、有効煙突高さ１２１４などの解析に用いた各種条
件を分類した情報と、これらの関係を示す対応表１２１
５から構成される。
【００２３】大気環境シミュレーションシステム１０の
モデル選定部１１はこれらデータベース１２１〜１２３
を用いることによって、入力条件３に応じた環境項目や
解析対象から、参考とする法令・指針やモデルの選定を
行い、過去の計算事例から有効煙突高さや、拡散パラメ
ータなどの値を選定する。
モデル選定部１１はこれらデータベース１２１〜１２３
を用いることによって、入力条件３に応じた環境項目や
解析対象から、参考とする法令・指針やモデルの選定を
行い、過去の計算事例から有効煙突高さや、拡散パラメ
ータなどの値を選定する。
【００２４】さて、次に入力パラメータチューニング部
１３の処理動作をモデル選定部１１とデータベース部１
２と関連付けて図３のフロー図を用いて説明する。
１３の処理動作をモデル選定部１１とデータベース部１
２と関連付けて図３のフロー図を用いて説明する。
【００２５】入力装置１から入力された入力条件３はモ
デル選定部１１に与えられ、入力条件Ａが選定される
（ステップＳ１）。選定された入力条件Ａは入力パラメ
ータチューニング部１３へ加えられる。入力パラメータ
チューニング部１３はこの条件Ａに基づく数値シミュレ
ーションの準備を行ない（ステップＳ２）、入力データ
作成から数値シミュレーションまでの処理を入力作成部
１４と数値計算処理部１５へ必要な情報を伝達して実施
する（処理１３１）。この処理１３１を図３では入力デ
ータ作成、数値シミュレーション処理と称している。
デル選定部１１に与えられ、入力条件Ａが選定される
（ステップＳ１）。選定された入力条件Ａは入力パラメ
ータチューニング部１３へ加えられる。入力パラメータ
チューニング部１３はこの条件Ａに基づく数値シミュレ
ーションの準備を行ない（ステップＳ２）、入力データ
作成から数値シミュレーションまでの処理を入力作成部
１４と数値計算処理部１５へ必要な情報を伝達して実施
する（処理１３１）。この処理１３１を図３では入力デ
ータ作成、数値シミュレーション処理と称している。
【００２６】この処理１３１の詳細を図４のフロー図を
用いて説明する。
用いて説明する。
【００２７】入力データ作成、数値シミュレーション部
１３１では選定した入力条件Ａを入力作成部１４に指示
し（ステップＳ３）、数値シミュレーション用の入力デ
ータを作成する（ステップＳ４）。この数値シミュレー
ション用データは、分析者が入力した入力条件３とは異
なり、数値シミュレーションで扱える物理データ、例え
ば大気状態は大気安定度という気象状況を表す指標から
水平方向と鉛直方向の拡散パラメータに変換されたもの
である。
１３１では選定した入力条件Ａを入力作成部１４に指示
し（ステップＳ３）、数値シミュレーション用の入力デ
ータを作成する（ステップＳ４）。この数値シミュレー
ション用データは、分析者が入力した入力条件３とは異
なり、数値シミュレーションで扱える物理データ、例え
ば大気状態は大気安定度という気象状況を表す指標から
水平方向と鉛直方向の拡散パラメータに変換されたもの
である。
【００２８】次に数値計算処理部１５へその入力データ
を指示し（ステップＳ５）、シミュレーション計算を行
い（ステップＳ６）、シミュレーション結果を選定した
入力条件結果として出力する。
を指示し（ステップＳ５）、シミュレーション計算を行
い（ステップＳ６）、シミュレーション結果を選定した
入力条件結果として出力する。
【００２９】図３に戻り入力データチューニング部１３
の動作を説明する。
の動作を説明する。
【００３０】図４で説明したように入力データ作成、数
値シミュレーション処理部１３１の処理結果として、入
力条件Aシミュレーション解析結果が得られる。また、
入力パラメータチューニング部１３には入力条件３が加
えられ、データベース部１２の測定値データベース１２
４を検索して入力条件Ａに応じた測定値、測定条件を検
索し（ステップＳ８）、測定値および測定条件１３３の
データを作成する。
値シミュレーション処理部１３１の処理結果として、入
力条件Aシミュレーション解析結果が得られる。また、
入力パラメータチューニング部１３には入力条件３が加
えられ、データベース部１２の測定値データベース１２
４を検索して入力条件Ａに応じた測定値、測定条件を検
索し（ステップＳ８）、測定値および測定条件１３３の
データを作成する。
【００３１】入力パラメータチューニング部１３では先
に求めた入力条件Aシミュレーション解析結果１３２と
測定値、測定条件１３３より、解析結果と測定値との比
較を行い（ステップＳ９）、感度解析のためのパラメー
タを選定の上チューニングする（ステップＳ１０）。感
度解析のパラメータには、例えば水平方向拡散パラメー
タσｙや鉛直方向拡散パラメータσｚ及び修正係数など
がある。
に求めた入力条件Aシミュレーション解析結果１３２と
測定値、測定条件１３３より、解析結果と測定値との比
較を行い（ステップＳ９）、感度解析のためのパラメー
タを選定の上チューニングする（ステップＳ１０）。感
度解析のパラメータには、例えば水平方向拡散パラメー
タσｙや鉛直方向拡散パラメータσｚ及び修正係数など
がある。
【００３２】次にそのパラメータを変更して入力データ
作成、数値シミュレーション（処理１３１）を行い、そ
の結果から誤差を求め（ステップＳ１１）、誤差が最小
かどうか最小２乗法などを用いて判定する（ステップＳ
１２）。その結果、最小でない場合は、パラメータの値
をチューニングし（ステップＳ１０）、同様にステップ
Ｓ１０からステップＳ１２までの処理を繰り返し実行す
る。誤差が最小になった場合は、そのパラメータを入力
条件（B）５の処理結果として入力データ作成部に与え
る（ステップＳ１３）。
作成、数値シミュレーション（処理１３１）を行い、そ
の結果から誤差を求め（ステップＳ１１）、誤差が最小
かどうか最小２乗法などを用いて判定する（ステップＳ
１２）。その結果、最小でない場合は、パラメータの値
をチューニングし（ステップＳ１０）、同様にステップ
Ｓ１０からステップＳ１２までの処理を繰り返し実行す
る。誤差が最小になった場合は、そのパラメータを入力
条件（B）５の処理結果として入力データ作成部に与え
る（ステップＳ１３）。
【００３３】このようにして入力パラメータチューニン
グ部１３で類似している環境における実測値と一致する
ように調整入力パラメータを選定することができ、より
精度の高い数値シミュレーションを行うことができる。
グ部１３で類似している環境における実測値と一致する
ように調整入力パラメータを選定することができ、より
精度の高い数値シミュレーションを行うことができる。
【００３４】次に発生源推定処理部１６の処理動作を入
力条件３、数値計算処理部１５および解析結果データベ
ース１８２と関連付けて図５のフロー図を用いて説明す
る。
力条件３、数値計算処理部１５および解析結果データベ
ース１８２と関連付けて図５のフロー図を用いて説明す
る。
【００３５】入力装置１から、入力条件３として格子状
に配列したｎ個の発生源位置（座標）、および１箇所ま
たは複数箇所における汚染物質の測定値（実測データ）
を与える（ステップＳ１）。入力条件３に基づき数値計
算処理部１５でｎ個の各発生源について数値シミュレー
ションを行い（ステップＳ２）、解析結果を解析結果デ
ータベース１８２に格納する。
に配列したｎ個の発生源位置（座標）、および１箇所ま
たは複数箇所における汚染物質の測定値（実測データ）
を与える（ステップＳ１）。入力条件３に基づき数値計
算処理部１５でｎ個の各発生源について数値シミュレー
ションを行い（ステップＳ２）、解析結果を解析結果デ
ータベース１８２に格納する。
【００３６】ｎ個の発生源位置ごとの数値シミュレーシ
ョン結果から測定値（実測データ）のある位置における
解析による汚染物質濃度を計算し（ステップＳ３）、測
定値（実測データ）がある位置における解析結果と測定
値の差の２乗和が最小となる格子点の座標を求めると、
これが汚染物質発生源と推定される。次に、推定された
汚染物質発生源位置は、表示用処理部２０の発生源位置
可視化データに変換し、表示精度地図データと合わせて
出力装置２に表示する。
ョン結果から測定値（実測データ）のある位置における
解析による汚染物質濃度を計算し（ステップＳ３）、測
定値（実測データ）がある位置における解析結果と測定
値の差の２乗和が最小となる格子点の座標を求めると、
これが汚染物質発生源と推定される。次に、推定された
汚染物質発生源位置は、表示用処理部２０の発生源位置
可視化データに変換し、表示精度地図データと合わせて
出力装置２に表示する。
【００３７】
【発明の効果】本発明によれば、解析条件に適したモデ
ルや入力パラメータを自動的に設定することができ、解
析者に高度な知識や経験がなくても数値シミュレーショ
ンを行い、さらに測定値（実測データ）がどの位置の発
生源から放出された汚染物質によるものかを推定するこ
とができる。
ルや入力パラメータを自動的に設定することができ、解
析者に高度な知識や経験がなくても数値シミュレーショ
ンを行い、さらに測定値（実測データ）がどの位置の発
生源から放出された汚染物質によるものかを推定するこ
とができる。
【図１】本発明の一実施例を示す構成図である。
【図２】データベースの形式を示す一例である。
【図３】入力パラメータチューニング部の動作フロー図
である。
である。
【図４】入力データ作成、数値シミュレーション処理の
動作フロー図である。
動作フロー図である。
【図５】発生源推定処理の動作フロー図である。
１…入力装置、２…出力装置、３…入力条件、４…入力
条件選定結果Ａ、５…入力条件選定結果Ｂ、６…数値計
算入力データ、７…表示条件１０…大気環境シミュレー
ションシステム、１１…モデル選定部、１２…データベ
ース部、１２１…解析事例データベース、１２１１…排
出源情報、１２１２…モデル、１２１３…拡散パラメー
タ、１２１４…有効煙突高さ、１２１５…対応表、１２
２…法令・指針データベース、１２２１…法令・指針
名、１２２２…環境項目、１２２３…解析対象分野、１
２２４…地域、１２２５…モデル、１２２６…対応表、
１２３…シミュレーションモデルデータベース、１２４
…測定値データベース、１３…入力パラメータチューニ
ング部、１３１…入力データ作成、数値シミュレーショ
ンの処理、１３２…入力条件Aシミュレーション計算結
果、１３３…測定値、測定条件、１４…入力作成部、１
５…数値計算処理部、１５１…対象領域全範囲の数値解
析結果、１６…発生源推定処理部、１６１…複数発生源
の数値解析結果の統計処理、１７…表示合成部、１８…
データ格納部、１８１…地図データベース、１８２…解
析結果データベース １９…可視化データ選定処理部、２０…表示用処理部、
２０１…表示精度地図データ、２０２…表示精度可視化
用物理量分布データ、２０３…物理量分布可視化デー
タ、２０４…発生源位置座標データ、２０５…発生源位
置座標可視化データ
条件選定結果Ａ、５…入力条件選定結果Ｂ、６…数値計
算入力データ、７…表示条件１０…大気環境シミュレー
ションシステム、１１…モデル選定部、１２…データベ
ース部、１２１…解析事例データベース、１２１１…排
出源情報、１２１２…モデル、１２１３…拡散パラメー
タ、１２１４…有効煙突高さ、１２１５…対応表、１２
２…法令・指針データベース、１２２１…法令・指針
名、１２２２…環境項目、１２２３…解析対象分野、１
２２４…地域、１２２５…モデル、１２２６…対応表、
１２３…シミュレーションモデルデータベース、１２４
…測定値データベース、１３…入力パラメータチューニ
ング部、１３１…入力データ作成、数値シミュレーショ
ンの処理、１３２…入力条件Aシミュレーション計算結
果、１３３…測定値、測定条件、１４…入力作成部、１
５…数値計算処理部、１５１…対象領域全範囲の数値解
析結果、１６…発生源推定処理部、１６１…複数発生源
の数値解析結果の統計処理、１７…表示合成部、１８…
データ格納部、１８１…地図データベース、１８２…解
析結果データベース １９…可視化データ選定処理部、２０…表示用処理部、
２０１…表示精度地図データ、２０２…表示精度可視化
用物理量分布データ、２０３…物理量分布可視化デー
タ、２０４…発生源位置座標データ、２０５…発生源位
置座標可視化データ
───────────────────────────────────────────────────── フロントページの続き (72)発明者 三好 由美子 茨城県日立市幸町三丁目２番１号 日立エ ンジニアリング株式会社内 Ｆターム(参考） 2C032 HB03 HB05 HB15 HC27 5B050 AA00 BA17 CA07 EA19 FA02
Claims (2)
- 【請求項１】大気汚染物質の発生源から放出される汚染
物質の拡散分布状況をシミュレーションするシミュレー
ション条件データを入力する入力手段と、解析事例、法
令・指針、シミュレーションモデルおよび測定値のそれ
ぞれのデータを有するデータベースと、前記シミュレー
ション条件データを入力してこのシミュレーションに適
したシミュレーションモデル、固定パラメータおよび解
析条件を前記データベースから選定するモデル選定手段
と、前記シミュレーション条件データを入力して拡散分
布状況の数値シミュレーションを行い、類似している環
境における実測値を格納している前記データベースの測
定値との誤差が小さくなるように調整パラメータを選定
するパラメータチューニング手段と、前記モデル選定手
段により選定されたシミュレーションモデル、固定パラ
メータ、解析条件と前記パラメータチューニング手段に
より選定された調整パラメータに基づいて得られた数値
計算データを用いて前記汚染物質の拡散分布状況の数値
シミュレーションを実行する数値計算処理手段と、前記
数値計算処理手段による解析結果データおよび地図デー
タを格納するデータ格納手段と、前記数値計算処理手段
による計算結果と前記汚染物質の実測値データに基づい
て少なくとも発生源の位置を含む発生源情報を推定する
発生源推定手段と、前記数値計算処理手段による計算結
果、前記地図データによる地図および発生源位置を表示
装置に表示する表示処理手段とを具備したことを特徴と
する大気拡散シミュレーションシステム。 - 【請求項２】請求項１において、前記発生源推定手段
は、前記入力手段により汚染物質発生源の位置、放出時
間、放出量を入力し、前記数値計算処理手段を用いて前
記汚染物質の拡散分布状況を数値シミュレーションした
計算結果を汚染物質発生源の位置をパラメータとして複
数回計算し、汚染物質濃度実測値と合わせて統計的手法
を用いたデータ解析を行い汚染物質発生源位置を推定す
ることを特徴とする大気環境シミュレーションシステ
ム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001002458A JP2002207002A (ja) | 2001-01-10 | 2001-01-10 | 大気環境シミュレーションシステム |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001002458A JP2002207002A (ja) | 2001-01-10 | 2001-01-10 | 大気環境シミュレーションシステム |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2002207002A true JP2002207002A (ja) | 2002-07-26 |
Family
ID=18870962
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001002458A Pending JP2002207002A (ja) | 2001-01-10 | 2001-01-10 | 大気環境シミュレーションシステム |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2002207002A (ja) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100851520B1 (ko) * | 2005-12-28 | 2008-08-11 | 한국건설기술연구원 | 통계적 분석기법을 이용한 지하 공간의 대기 환경 예측장치 및 방법 |
| JP2013185988A (ja) * | 2012-03-08 | 2013-09-19 | Toshiba Corp | 情報処理装置及び放出点推定方法 |
| WO2014020747A1 (ja) * | 2012-08-02 | 2014-02-06 | 株式会社日立製作所 | 物理量推定システム及び物理量推定方法 |
| JP2015172494A (ja) * | 2014-03-11 | 2015-10-01 | 新日鐵住金株式会社 | 環境負荷分子発生源評価方法、環境負荷分子発生源評価システム、及びコンピュータプログラム |
| US9326969B2 (en) | 2012-09-11 | 2016-05-03 | Norgine Bv | Compositions |
| CN111443015A (zh) * | 2020-03-04 | 2020-07-24 | 平安国际智慧城市科技股份有限公司 | 大气污染物的溯源方法、装置、计算机设备及存储介质 |
| CN112986072A (zh) * | 2021-02-09 | 2021-06-18 | 上海英凡环保科技有限公司 | 污染源臭气在线监测系统及方法 |
-
2001
- 2001-01-10 JP JP2001002458A patent/JP2002207002A/ja active Pending
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100851520B1 (ko) * | 2005-12-28 | 2008-08-11 | 한국건설기술연구원 | 통계적 분석기법을 이용한 지하 공간의 대기 환경 예측장치 및 방법 |
| JP2013185988A (ja) * | 2012-03-08 | 2013-09-19 | Toshiba Corp | 情報処理装置及び放出点推定方法 |
| WO2014020747A1 (ja) * | 2012-08-02 | 2014-02-06 | 株式会社日立製作所 | 物理量推定システム及び物理量推定方法 |
| JPWO2014020747A1 (ja) * | 2012-08-02 | 2016-07-11 | 株式会社日立製作所 | 物理量推定システム及び物理量推定方法 |
| US9326969B2 (en) | 2012-09-11 | 2016-05-03 | Norgine Bv | Compositions |
| JP2015172494A (ja) * | 2014-03-11 | 2015-10-01 | 新日鐵住金株式会社 | 環境負荷分子発生源評価方法、環境負荷分子発生源評価システム、及びコンピュータプログラム |
| CN111443015A (zh) * | 2020-03-04 | 2020-07-24 | 平安国际智慧城市科技股份有限公司 | 大气污染物的溯源方法、装置、计算机设备及存储介质 |
| CN111443015B (zh) * | 2020-03-04 | 2024-09-13 | 平安国际智慧城市科技股份有限公司 | 大气污染物的溯源方法、装置、计算机设备及存储介质 |
| CN112986072A (zh) * | 2021-02-09 | 2021-06-18 | 上海英凡环保科技有限公司 | 污染源臭气在线监测系统及方法 |
| CN112986072B (zh) * | 2021-02-09 | 2022-06-03 | 上海英凡环保科技有限公司 | 污染源臭气在线监测系统及方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Haylock et al. | Downscaling heavy precipitation over the United Kingdom: a comparison of dynamical and statistical methods and their future scenarios | |
| US8903646B2 (en) | Method and determining system for automatically determining emission locations, and method and traffic control system based thereon for immission-dependent traffic control | |
| KR101668753B1 (ko) | 대기 오염 관리 시스템 | |
| CN112308292A (zh) | 一种火灾风险等级分布图的绘制方法 | |
| Zhou et al. | Ensemble transform Kalman filter (ETKF) for large-scale wildland fire spread simulation using FARSITE tool and state estimation method | |
| JP5743930B2 (ja) | 大気拡散物質発生源探査装置、大気拡散物質発生源探査システムおよび大気拡散物質発生源探査方法 | |
| CN110457829A (zh) | 一种基于集成大气扩散模型的源项释放反演和扩散预测方法 | |
| Lesieur et al. | Data assimilation for urban noise mapping with a meta-model | |
| US8665270B2 (en) | Method of determining implicit hidden features of phenomena which can be represented by a point distribution in a space | |
| JP4206024B2 (ja) | 大気物質指標分布解析装置 | |
| JP2000194742A (ja) | 大気環境シミュレ―ションシステム | |
| CN111461163B (zh) | 城市内部pm2.5浓度模拟和人口暴露度评估方法和装置 | |
| Crenna | An introduction to WindTrax | |
| CN112989153B (zh) | 数据处理方法、装置及计算机设备 | |
| Mileu et al. | Development of a QGIS plugin to dasymetric mapping | |
| JP4141537B2 (ja) | 建築用ｃａｄシステム | |
| JPWO2010092680A1 (ja) | 地図情報処理装置、地図情報処理方法、地図情報処理プログラム、及び記録媒体 | |
| Jakobsson | Data quality and quality management–examples of quality evaluation procedures and quality management in European national mapping agencies | |
| KR102660454B1 (ko) | 대기질 예측 모델에 의한 예측 데이터 및 관측 데이터의 관계를 제공하는 방법 | |
| KR20050004316A (ko) | 복수 기상관측자료를 이용한 대기오염 확산 예측 방법 | |
| Kedrin et al. | System for visualization the map of air pollution by vehicles | |
| CN112598763A (zh) | 一种区域地表温度场的时间过程模拟方法及系统 | |
| JP5433559B2 (ja) | 拡散物質の発生源推定装置および発生源推定方法 | |
| CN117910701A (zh) | 基于ftir的大气污染物溯源和臭氧生成潜势评估方法及系统 | |
| Heyns | Reduced target-resolution strategy for rapid multi-observer site location optimisation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20040622 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20041102 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20041224 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050329 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050517 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20051108 |