JP2002207002A - Atmospheric environment simulation system - Google Patents

Atmospheric environment simulation system

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
Application number
JP2001002458A
Other languages
Japanese (ja)
Inventor
Tetsuaki Matsuura
哲明 松浦
Asao Yamamoto
朝男 山本
Yumiko Miyoshi
由美子 三好
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Engineering Co Ltd
Original Assignee
Hitachi Engineering Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hitachi Engineering Co Ltd filed Critical Hitachi Engineering Co Ltd
Priority to JP2001002458A priority Critical patent/JP2002207002A/en
Publication of JP2002207002A publication Critical patent/JP2002207002A/en
Pending legal-status Critical Current

Links

Landscapes

  • Instructional Devices (AREA)
  • Processing Or Creating Images (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

PROBLEM TO BE SOLVED: To provide an atmospheric environment simulation system capable of estimating the position of a pollutant generation source by combining the measured value (actually measured data) of an atmospheric pollutant and a calculation result obtained by the numerical simulation of the diffusion distribution state of the pollutants. SOLUTION: A numerical value calculation processing part 15 performs the simulation of the diffusion distribution state of the pollutants using numerical value calculation data obtained on the basis of a selected simulation model, a fixed parameter, an analyzing condition and an adjusting parameter. A series of analyses displaying the simulation result and a map are repeatedly performed by the number of pollutant generation sources, and these results and the measured value of the atmospheric pollutants are used to estimate the position of the atmospheric polluting substance generation source according to a statistic technique by a generation source estimation means 16.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は大気汚染物質の発生
源による大気汚染物質の拡散分布状態を解析した結果と
大気汚染濃度の実測値(測定データ)とを組合わせて、
統計的手法を用いたデータ解析を行うことにより、実測
データに寄与している発生源の位置等を推定し、対策検
討のための情報を提供する大気環境シミュレーションシ
ステムに関する。
BACKGROUND OF THE INVENTION The present invention combines the results of analyzing the distribution distribution of air pollutants due to the source of air pollutants with the actual measured values (measurement data) of the air pollutant concentration.
The present invention relates to an atmospheric environment simulation system that estimates the position of a source contributing to actual measurement data by performing data analysis using a statistical method and provides information for studying countermeasures.

【0002】[0002]

【従来の技術】近年、地球環境保全の観点から環境アセ
スメント法の制定やISO14000シリーズの発行など、大気
環境の改善や大気汚染対策のため環境影響評価が重要に
なってきている。特に大気環境に汚染物質を放出する工
場や発電所、燃焼施設では計画・立案段階は勿論のこ
と、運転稼働中も環境影響を定期的に把握する必要があ
る。
2. Description of the Related Art In recent years, environmental impact assessment has become important for improving the air environment and taking measures against air pollution, such as the establishment of the Environmental Assessment Law and the issuance of the ISO14000 series from the viewpoint of global environmental protection. Particularly in factories, power plants and combustion facilities that emit pollutants into the air environment, it is necessary to periodically monitor the environmental impact during operation, as well as during the planning and planning stages.

【0003】従来、大気汚染物質を放出する工場や発電
所、燃焼施設等からの汚染物質の拡散予測を行い、予測
結果を表示装置上に表示するシステムとしては、例えば
特開平11−14531号公報が挙げられる。このよう
な予測システムは発生源からの汚染物質の拡散分布を計
算するものである。
[0003] Conventionally, as a system for predicting the diffusion of pollutants from factories, power plants, combustion facilities, etc. that emit air pollutants, and displaying the prediction results on a display device, for example, JP-A-11-14531 is known. Is mentioned. Such prediction systems calculate the diffusion distribution of contaminants from the source.

【0004】[0004]

【発明が解決しようとする課題】しかしながら、例えば
実測された大気汚染物質が複数の発生源のどれから放出
されたものか、すなわち汚染物質濃度の測定値(実測デ
ータ)から汚染物質発生源の位置等を推定することは困
難であった。
However, for example, from which of a plurality of sources the actually measured air pollutants are emitted, that is, from the measured value (actually measured data) of the pollutant concentration, the position of the pollutant source is determined. It was difficult to estimate such.

【0005】本発明の目的は、大気汚染物質の測定値
(実測データ)と、汚染物質の拡散分布状況を数値シミ
ュレーションした計算結果とを組合わせて汚染物質発生
源の位置を推定することができる大気環境シミュレーシ
ョンシステムを提供することにある。
An object of the present invention is to estimate the position of a pollutant source by combining a measured value (actually measured data) of an air pollutant and a calculation result obtained by numerically simulating the state of diffusion distribution of the pollutant. An object of the present invention is to provide an atmospheric environment simulation system.

【0006】[0006]

【課題を解決するための手段】本発明は上記目的を達成
するために、大気汚染物質の発生源位置を表示画面上に
格子状に規則正しく複数個設定し、各発生源から放出さ
れる汚染物質の拡散分布状況をシミュレーションするシ
ミュレーション条件データに基づきこのシミュレーショ
ンに適したシミュレーションモデル、固定パラメータ及
び解析条件を選定するとともに、シミュレーション条件
データにより拡散分布状況の数値シミュレーションを行
い類似している環境における実測値(測定値)との誤差
が小さくなるように調整パラメータを選定し、選定され
たシミュレーションモデル、固定パラメータ、解析条件
および調整パラメータに基づいて得られた数値計算デー
タを用いて汚染物質の拡散分布状況の数値シミュレーシ
ョンを実行してその結果と地図を表示する一連の解析を
発生源の個数だけ繰り返し実施し、それらの結果と大気
汚染物質の測定値とを統計的手法を用いて計算すること
により大気汚染物質の発生源位置を推定して表示装置に
表示するようにしたことを特徴とする。
According to the present invention, in order to achieve the above object, a plurality of sources of air pollutants are regularly set in a grid pattern on a display screen, and contaminants emitted from each source are set. A simulation model, fixed parameters and analysis conditions suitable for this simulation are selected based on the simulation condition data that simulates the diffusion distribution situation of the simulation, and a numerical simulation of the diffusion distribution situation is performed using the simulation condition data, and the actual measurement values in a similar environment Adjustment parameters are selected so that the error from (measured values) is small, and the diffusion distribution of contaminants is calculated using the selected simulation model, fixed parameters, analysis conditions, and numerical calculation data obtained based on the adjustment parameters. Run a numerical simulation of Estimate the source location of air pollutants by repeating a series of analyzes that display the results and maps for the number of sources and calculating the results and measured values of air pollutants using statistical methods And displaying it on a display device.

【0007】[0007]

【発明の実施の形態】以下、本発明の実施の形態を、図
面を参照して説明する。
Embodiments of the present invention will be described below with reference to the drawings.

【0008】図1に本発明の一実施例に係る大気環境シ
ミュレーションシステムの構成図を示す。
FIG. 1 shows a configuration diagram of an atmospheric environment simulation system according to one embodiment of the present invention.

【0009】図1において、大気環境シミュレーション
システム10は、入力データ(シミュレーション条件デ
ータ)を入力する入力装置1と解析結果を表示,出力す
るディスプレイやプリンタ等の出力装置2とを結合する
もので、モデル選定部11、データベース部12、入力
パラメータチューニング部13、入力作成部14、数値
計算処理部15、発生源推定処理部16、表示合成部1
7、データ格納部18、可視化データ選定処理部19、
表示用処理部20から構成されている。
In FIG. 1, an atmospheric environment simulation system 10 is a combination of an input device 1 for inputting input data (simulation condition data) and an output device 2 such as a display or a printer for displaying and outputting analysis results. Model selection unit 11, database unit 12, input parameter tuning unit 13, input creation unit 14, numerical calculation processing unit 15, generation source estimation processing unit 16, display synthesis unit 1
7, data storage unit 18, visualization data selection processing unit 19,
It comprises a display processing unit 20.

【0010】モデル選定部11には入力装置1から汚染
物質の種類や放出量、放出位置、大気条件などの入力条
件3が入力される。モデル選定部11は入力された入力
条件3に基づき汚染物質の種類や形状,放出位置,放出
条件,評価時点での大気条件,評価地点の地形状況、等
の観点から分類を行い、データベース部12の中に格納
されている解析事例データベース121、法令・指針デ
ータベース122およびシミュレーションモデルデータ
ベース123を参照することによって、入力条件選定部
(A)4において入力条件3に最適なシミュレーション
モデル41と固定パラメータ42と解析条件43を選定
する。
The input condition 1 such as the type and release amount of the pollutant, the release position, the atmospheric conditions, etc. is input from the input device 1 to the model selection unit 11. The model selection unit 11 performs classification based on the input condition 3 based on the type and shape of the pollutant, the release position, the release condition, the atmospheric condition at the time of evaluation, the topographical condition of the evaluation point, and the like, and the database unit 12 By referring to the analysis case database 121, the law / guideline database 122, and the simulation model database 123 stored in, a simulation model 41 and a fixed parameter 42 that are optimal for the input condition 3 in the input condition selecting unit (A) 4. And the analysis condition 43 are selected.

【0011】入力条件選定部(B)5は入力条件3に基
づく解析結果の精度を向上させ、実測値(測定値)と近
づけるため、入力パラメータチューニング部13でデー
タベース部12の中に格納されている測定値データベー
ス124を参照して調整パラメータ51を選定し、入力
作成部14に与える。入力作成部14は選定した条件を
基に、数値計算処理部15のシミュレーション計算に必
要な数値計算入力データ6を作成する。
The input condition selection unit (B) 5 is stored in the database unit 12 by the input parameter tuning unit 13 in order to improve the accuracy of the analysis result based on the input condition 3 and bring it closer to the actually measured value (measured value). The adjustment parameter 51 is selected with reference to the measured value database 124 and given to the input creation unit 14. The input creation unit 14 creates the numerical calculation input data 6 necessary for the simulation calculation of the numerical calculation processing unit 15 based on the selected conditions.

【0012】数値計算処理部15は、入力作成部14で
作成された数値計算入力データを用い、内蔵した物理モ
デルから解析対象領域全範囲にわたる物理現象を並列計
算等の高速化手法を用いて数値解析を実施し、その対象
領域全範囲の高精度数値解析結果151を表示合成部1
7のデータ格納部18にある解析結果データベース18
2の中に分類整理して格納する。
The numerical calculation processing unit 15 uses the numerical calculation input data created by the input creating unit 14 to numerically analyze physical phenomena over the entire range of the analysis target area from a built-in physical model by using a parallel calculation or other high-speed technique. Analysis is performed, and a high-precision numerical analysis result 151 of the entire target area is displayed and displayed by the display combining unit 1.
7 analysis result database 18 in the data storage unit 18
2 and sort and arrange them.

【0013】表示合成部17は、可視化に必要なデータ
を格納しておくデータ格納部18と可視化に必要なデー
タを選定する可視化データ選定処理部19と表示用処理
部20から構成される。
The display synthesizing section 17 comprises a data storage section 18 for storing data necessary for visualization, a visualization data selection processing section 19 for selecting data necessary for visualization, and a display processing section 20.

【0014】データ格納部18には解析結果を格納する
解析結果データベース182の他に解析結果表示の際に
背景データとして使用する地図データを格納した、地図
データベース181を有する。地図データベース181
は解析対象領域全範囲に関する多種類の精度の地図デー
タが格納されている。
The data storage unit 18 has a map database 181 storing map data used as background data when displaying the analysis results, in addition to an analysis result database 182 storing the analysis results. Map database 181
Stores map data of various kinds of accuracy with respect to the entire range of the analysis target area.

【0015】可視化データ選定処理部19は入力装置1
から入力された表示条件7に基づきデータ格納部18の
データベース181、182を参照し、表示用処理部2
0へ表示条件に適した表示精度地図データ201と表示
精度可視化用物理量分布データ202を選定し指示す
る。表示用処理部20は可視化データ選定処理部19で
指示された表示精度地図データ201と表示精度可視化
用物理量分布データ202を可視化処理を行い、合成し
て出力装置2に描画する物理量分布可視化データ203
を作成する。
The visualization data selection processing unit 19 includes the input device 1
The display processing unit 2 refers to the databases 181 and 182 of the data storage unit 18 based on the display condition 7 input from the
To 0, the display accuracy map data 201 and the display accuracy visualization physical quantity distribution data 202 suitable for the display conditions are selected and instructed. The display processing unit 20 performs a visualization process on the display accuracy map data 201 and the display accuracy visualization physical quantity distribution data 202 specified by the visualization data selection processing unit 19, and combines and displays the physical quantity distribution visualization data 203 to be drawn on the output device 2.
Create

【0016】これまでの一連の操作を、入力条件3にお
いて汚染物質の発生源位置のみを格子状に規則正しく複
数個設定し、各発生源条件ごとに数値計算処理部151
において複数回繰り返し計算させ、複数回の計算結果を
データ格納部18の解析結果データベース182に格納
する。汚染物質発生源位置が異なる前記複数の計算結果
および入力条件3であらかじめ入力された実測値(測定
値)に基づき、発生源推定処理部16において統計的手
法を用いて計算処理を行い、汚染物質の発生源として最
も確からしい位置を推定し、推定結果を再び解析結果デ
ータベース182に格納する。
In a series of operations up to this point, only a plurality of source positions of contaminants are regularly set in a grid in input condition 3, and a numerical calculation processing unit 151 is set for each source condition.
The calculation is repeatedly performed a plurality of times, and the calculation results of the plurality of times are stored in the analysis result database 182 of the data storage unit 18. Based on the plurality of calculation results at different contaminant source positions and the actually measured values (measured values) input in advance under the input condition 3, the source estimation processing unit 16 performs a calculation process using a statistical method, The most probable position is estimated as the generation source of, and the estimation result is stored in the analysis result database 182 again.

【0017】このように、入力装置1から格子状に規則
正しく配置した汚染物質発生源位置、汚染物質濃度が実
測されている位置及び濃度、並びにその他の入力条件を
入力条件3として入力することにより、複数の発生源位
置に対する汚染物質の拡散分布状況を数値計算処理部1
5において計算し、複数の計算結果から汚染物質濃度が
実測されている位置における計算値を算出し、当該位置
における汚染物質濃度の実測値(測定データ)と比較す
ることにより、どの位置の汚染物質発生源から放出され
た可能性が高いのかを容易に推定することができる。
As described above, by inputting, from the input device 1, the positions of the contaminant generating sources regularly arranged in a grid pattern, the positions and concentrations where the contaminant concentrations are actually measured, and other input conditions, the input conditions 3 are used. Numerical calculation processing unit 1
5 to calculate the calculated value at the position where the pollutant concentration is actually measured from the plurality of calculation results, and compare the calculated value with the actually measured value (measurement data) of the pollutant concentration at the position. It can be easily estimated that the possibility of release from the source is high.

【0018】さらに、表示条件7を可視化データ選定処
理部19で表示範囲に応じた表示精度可視化用物理量分
布データ202と表示精度地図データ201を、また発
生源位置座標データ204と表示精度地図データ201
を、をそれぞれ合成して表示することができる。
Further, the display condition 7 is displayed by the visualization data selection processing unit 19 on the display accuracy visualization physical quantity distribution data 202 and the display accuracy map data 201 according to the display range, and on the source position coordinate data 204 and the display accuracy map data 201.
Can be combined and displayed.

【0019】次にデータベース部に格納されるデータに
ついて説明する。
Next, data stored in the database unit will be described.

【0020】図2に法令・指針データベース122に格
納されているデータの一例を示す。
FIG. 2 shows an example of data stored in the law / guideline database 122.

【0021】図2(a)は法令・指針データベース12
2の例を示す。法令・指針データベース122には、法
令・指針名1221と環境項目1222、解析対象分野
1223、地域1224、モデル1225に分類された
情報があり、これらの関係を示す対応表1226から構
成される。対応表1226には、それぞれの情報の関係
だけでなくその重要性も数値化して持っている。この事
例では重要性を二重丸、丸、三角の記号で表している。
FIG. 2A shows a law / guideline database 12.
2 shows an example. The law / guideline database 122 has information classified into a law / guideline name 1221, an environmental item 1222, an analysis target field 1223, a region 1224, and a model 1225, and is configured from a correspondence table 1226 indicating a relationship between these. The correspondence table 1226 quantifies not only the relationship between the respective information but also its importance. In this case, the significance is represented by double circles, circles, and triangles.

【0022】図2(b)は解析事例データベース121
の一例を示す。解析事例データベース121には、排出
源情報1211、モデル1212、拡散パラメータ12
13、有効煙突高さ1214などの解析に用いた各種条
件を分類した情報と、これらの関係を示す対応表121
5から構成される。
FIG. 2B shows an analysis case database 121.
An example is shown below. The analysis case database 121 includes emission source information 1211, a model 1212, a diffusion parameter 12
13. Correspondence table 121 showing the information obtained by classifying various conditions used for the analysis of the effective chimney height 1214 and the like, and the relationship between them.
5 is comprised.

【0023】大気環境シミュレーションシステム10の
モデル選定部11はこれらデータベース121〜123
を用いることによって、入力条件3に応じた環境項目や
解析対象から、参考とする法令・指針やモデルの選定を
行い、過去の計算事例から有効煙突高さや、拡散パラメ
ータなどの値を選定する。
The model selecting section 11 of the atmospheric environment simulation system 10 stores these databases 121 to 123
By using the above, a law, a guideline or a model to be referred to is selected from environmental items and analysis targets according to the input condition 3, and values such as an effective stack height and a diffusion parameter are selected from past calculation examples.

【0024】さて、次に入力パラメータチューニング部
13の処理動作をモデル選定部11とデータベース部1
2と関連付けて図3のフロー図を用いて説明する。
Next, the processing operation of the input parameter tuning unit 13 will be described with reference to the model selecting unit 11 and the database unit 1.
2 will be described with reference to the flowchart of FIG.

【0025】入力装置1から入力された入力条件3はモ
デル選定部11に与えられ、入力条件Aが選定される
(ステップS1)。選定された入力条件Aは入力パラメ
ータチューニング部13へ加えられる。入力パラメータ
チューニング部13はこの条件Aに基づく数値シミュレ
ーションの準備を行ない(ステップS2)、入力データ
作成から数値シミュレーションまでの処理を入力作成部
14と数値計算処理部15へ必要な情報を伝達して実施
する(処理131)。この処理131を図3では入力デ
ータ作成、数値シミュレーション処理と称している。
The input condition 3 input from the input device 1 is given to the model selecting section 11, and the input condition A is selected (step S1). The selected input condition A is applied to the input parameter tuning unit 13. The input parameter tuning unit 13 prepares for a numerical simulation based on the condition A (step S2), and transmits necessary information to the input creating unit 14 and the numerical calculation processing unit 15 for processing from input data creation to numerical simulation. (Step 131). This processing 131 is called input data creation and numerical simulation processing in FIG.

【0026】この処理131の詳細を図4のフロー図を
用いて説明する。
The details of the process 131 will be described with reference to the flowchart of FIG.

【0027】入力データ作成、数値シミュレーション部
131では選定した入力条件Aを入力作成部14に指示
し(ステップS3)、数値シミュレーション用の入力デ
ータを作成する(ステップS4)。この数値シミュレー
ション用データは、分析者が入力した入力条件3とは異
なり、数値シミュレーションで扱える物理データ、例え
ば大気状態は大気安定度という気象状況を表す指標から
水平方向と鉛直方向の拡散パラメータに変換されたもの
である。
The input data creation and numerical simulation section 131 instructs the input creation section 14 on the selected input condition A (step S3), and creates input data for a numerical simulation (step S4). This data for numerical simulation is different from the input condition 3 input by the analyst, and physical data that can be handled by numerical simulation, for example, the atmospheric state is converted from an index indicating a weather condition such as atmospheric stability into diffusion parameters in the horizontal and vertical directions. It was done.

【0028】次に数値計算処理部15へその入力データ
を指示し(ステップS5)、シミュレーション計算を行
い(ステップS6)、シミュレーション結果を選定した
入力条件結果として出力する。
Next, the input data is instructed to the numerical calculation processing unit 15 (step S5), a simulation calculation is performed (step S6), and the simulation result is output as a selected input condition result.

【0029】図3に戻り入力データチューニング部13
の動作を説明する。
Returning to FIG. 3, the input data tuning unit 13
Will be described.

【0030】図4で説明したように入力データ作成、数
値シミュレーション処理部131の処理結果として、入
力条件Aシミュレーション解析結果が得られる。また、
入力パラメータチューニング部13には入力条件3が加
えられ、データベース部12の測定値データベース12
4を検索して入力条件Aに応じた測定値、測定条件を検
索し(ステップS8)、測定値および測定条件133の
データを作成する。
As described with reference to FIG. 4, an input condition A simulation analysis result is obtained as a processing result of the input data creation and numerical simulation processing section 131. Also,
The input condition tuning unit 13 receives the input condition 3 and the measured value database 12 of the database unit 12.
4 to retrieve a measurement value and a measurement condition corresponding to the input condition A (step S8), and create data of the measurement value and the measurement condition 133.

【0031】入力パラメータチューニング部13では先
に求めた入力条件Aシミュレーション解析結果132と
測定値、測定条件133より、解析結果と測定値との比
較を行い(ステップS9)、感度解析のためのパラメー
タを選定の上チューニングする(ステップS10)。感
度解析のパラメータには、例えば水平方向拡散パラメー
タσyや鉛直方向拡散パラメータσz及び修正係数など
がある。
The input parameter tuning unit 13 compares the analysis result with the measured value based on the previously obtained input condition A simulation analysis result 132 and the measured value and the measurement condition 133 (step S9), and sets parameters for sensitivity analysis. Is selected and tuned (step S10). The parameters for the sensitivity analysis include, for example, a horizontal diffusion parameter σy, a vertical diffusion parameter σz, and a correction coefficient.

【0032】次にそのパラメータを変更して入力データ
作成、数値シミュレーション(処理131)を行い、そ
の結果から誤差を求め(ステップS11)、誤差が最小
かどうか最小2乗法などを用いて判定する(ステップS
12)。その結果、最小でない場合は、パラメータの値
をチューニングし(ステップS10)、同様にステップ
S10からステップS12までの処理を繰り返し実行す
る。誤差が最小になった場合は、そのパラメータを入力
条件(B)5の処理結果として入力データ作成部に与え
る(ステップS13)。
Next, the parameters are changed, input data is created, and a numerical simulation (process 131) is performed. An error is obtained from the result (step S11), and it is determined whether or not the error is minimum by using a least square method or the like (step S11). Step S
12). As a result, if it is not the minimum, the value of the parameter is tuned (step S10), and the processing from step S10 to step S12 is similarly executed repeatedly. When the error is minimized, the parameter is given to the input data creating unit as the processing result of the input condition (B) 5 (step S13).

【0033】このようにして入力パラメータチューニン
グ部13で類似している環境における実測値と一致する
ように調整入力パラメータを選定することができ、より
精度の高い数値シミュレーションを行うことができる。
As described above, the input parameter tuning unit 13 can select the adjustment input parameter so as to match the actually measured value in a similar environment, and a more accurate numerical simulation can be performed.

【0034】次に発生源推定処理部16の処理動作を入
力条件3、数値計算処理部15および解析結果データベ
ース182と関連付けて図5のフロー図を用いて説明す
る。
Next, the processing operation of the source estimation processing unit 16 will be described with reference to the flowchart of FIG. 5 in association with the input condition 3, the numerical calculation processing unit 15, and the analysis result database 182.

【0035】入力装置1から、入力条件3として格子状
に配列したn個の発生源位置(座標)、および1箇所ま
たは複数箇所における汚染物質の測定値(実測データ)
を与える(ステップS1)。入力条件3に基づき数値計
算処理部15でn個の各発生源について数値シミュレー
ションを行い(ステップS2)、解析結果を解析結果デ
ータベース182に格納する。
From the input device 1, n source positions (coordinates) arranged in a grid as input conditions 3 and measured values (actually measured data) of contaminants at one or a plurality of locations
(Step S1). Based on the input condition 3, the numerical calculation processing unit 15 performs a numerical simulation on each of the n sources (step S2), and stores the analysis result in the analysis result database 182.

【0036】n個の発生源位置ごとの数値シミュレーシ
ョン結果から測定値(実測データ)のある位置における
解析による汚染物質濃度を計算し(ステップS3)、測
定値(実測データ)がある位置における解析結果と測定
値の差の2乗和が最小となる格子点の座標を求めると、
これが汚染物質発生源と推定される。次に、推定された
汚染物質発生源位置は、表示用処理部20の発生源位置
可視化データに変換し、表示精度地図データと合わせて
出力装置2に表示する。
From the numerical simulation results for each of the n source positions, the contaminant concentration at the position where the measured value (actually measured data) is present is calculated (step S3), and the analysis result at the position where the measured value (actually measured data) is present Finding the coordinates of the grid point that minimizes the sum of squares of the difference between
This is estimated to be a source of pollutants. Next, the estimated contaminant source position is converted into source position visualization data of the display processing unit 20, and displayed on the output device 2 together with the display accuracy map data.

【0037】[0037]

【発明の効果】本発明によれば、解析条件に適したモデ
ルや入力パラメータを自動的に設定することができ、解
析者に高度な知識や経験がなくても数値シミュレーショ
ンを行い、さらに測定値(実測データ)がどの位置の発
生源から放出された汚染物質によるものかを推定するこ
とができる。
According to the present invention, a model and an input parameter suitable for an analysis condition can be automatically set, a numerical simulation can be performed without an advanced knowledge and experience of an analyst, and a measured value can be obtained. It is possible to estimate from which source the (measured data) is due to the pollutants released.

【図面の簡単な説明】[Brief description of the drawings]

【図1】本発明の一実施例を示す構成図である。FIG. 1 is a configuration diagram showing one embodiment of the present invention.

【図2】データベースの形式を示す一例である。FIG. 2 is an example showing a format of a database.

【図3】入力パラメータチューニング部の動作フロー図
である。
FIG. 3 is an operation flowchart of an input parameter tuning unit.

【図4】入力データ作成、数値シミュレーション処理の
動作フロー図である。
FIG. 4 is an operation flowchart of input data creation and numerical simulation processing.

【図5】発生源推定処理の動作フロー図である。FIG. 5 is an operation flowchart of a generation source estimation process.

【符号の説明】[Explanation of symbols]

1…入力装置、2…出力装置、3…入力条件、4…入力
条件選定結果A、5…入力条件選定結果B、6…数値計
算入力データ、7…表示条件10…大気環境シミュレー
ションシステム、11…モデル選定部、12…データベ
ース部、121…解析事例データベース、1211…排
出源情報、1212…モデル、1213…拡散パラメー
タ、1214…有効煙突高さ、1215…対応表、12
2…法令・指針データベース、1221…法令・指針
名、1222…環境項目、1223…解析対象分野、1
224…地域、1225…モデル、1226…対応表、
123…シミュレーションモデルデータベース、124
…測定値データベース、13…入力パラメータチューニ
ング部、131…入力データ作成、数値シミュレーショ
ンの処理、132…入力条件Aシミュレーション計算結
果、133…測定値、測定条件、14…入力作成部、1
5…数値計算処理部、151…対象領域全範囲の数値解
析結果、16…発生源推定処理部、161…複数発生源
の数値解析結果の統計処理、17…表示合成部、18…
データ格納部、181…地図データベース、182…解
析結果データベース 19…可視化データ選定処理部、20…表示用処理部、
201…表示精度地図データ、202…表示精度可視化
用物理量分布データ、203…物理量分布可視化デー
タ、204…発生源位置座標データ、205…発生源位
置座標可視化データ
DESCRIPTION OF SYMBOLS 1 ... Input device, 2 ... Output device, 3 ... Input condition, 4 ... Input condition selection result A, 5 ... Input condition selection result B, 6 ... Numerical calculation input data, 7 ... Display condition 10 ... Atmospheric environment simulation system, 11 ... Model selection part, 12 ... Database part, 121 ... Analysis case database, 1211 ... Emission source information, 1212 ... Model, 1213 ... Diffusion parameter, 1214 ... Effective stack height, 1215 ... Correspondence table, 12
2 Law / guideline database, 1221 Law / guide name, 1222 environmental item, 1223 analysis target field, 1
224 area, 1225 model, 1226 correspondence table,
123: simulation model database, 124
... Measured value database, 13 ... Input parameter tuning unit, 131 ... Input data creation, numerical simulation processing, 132 ... Input condition A simulation calculation result, 133 ... Measured value, measurement condition, 14 ... Input creation unit, 1
5 Numerical calculation processing unit, 151 Numerical analysis result of entire target area, 16 Source generation estimation unit, 161 Statistical processing of numerical analysis results of multiple sources, 17 Display synthesis unit, 18
Data storage unit, 181: Map database, 182: Analysis result database 19: Visualization data selection processing unit, 20: Display processing unit,
201: display accuracy map data, 202: physical quantity distribution data for visualizing display accuracy, 203: physical quantity distribution visualization data, 204: source position coordinate data, 205: source position coordinate visualization data

───────────────────────────────────────────────────── フロントページの続き (72)発明者 三好 由美子 茨城県日立市幸町三丁目2番1号 日立エ ンジニアリング株式会社内 Fターム(参考) 2C032 HB03 HB05 HB15 HC27 5B050 AA00 BA17 CA07 EA19 FA02 ────────────────────────────────────────────────── ─── Continuing on the front page (72) Inventor Yumiko Miyoshi 3-2-1 Sachimachi, Hitachi-shi, Ibaraki F-term within Hitachi Engineering Co., Ltd. (reference) 2C032 HB03 HB05 HB15 HC27 5B050 AA00 BA17 CA07 EA19 FA02

Claims (2)

【特許請求の範囲】[Claims] 【請求項1】大気汚染物質の発生源から放出される汚染
物質の拡散分布状況をシミュレーションするシミュレー
ション条件データを入力する入力手段と、解析事例、法
令・指針、シミュレーションモデルおよび測定値のそれ
ぞれのデータを有するデータベースと、前記シミュレー
ション条件データを入力してこのシミュレーションに適
したシミュレーションモデル、固定パラメータおよび解
析条件を前記データベースから選定するモデル選定手段
と、前記シミュレーション条件データを入力して拡散分
布状況の数値シミュレーションを行い、類似している環
境における実測値を格納している前記データベースの測
定値との誤差が小さくなるように調整パラメータを選定
するパラメータチューニング手段と、前記モデル選定手
段により選定されたシミュレーションモデル、固定パラ
メータ、解析条件と前記パラメータチューニング手段に
より選定された調整パラメータに基づいて得られた数値
計算データを用いて前記汚染物質の拡散分布状況の数値
シミュレーションを実行する数値計算処理手段と、前記
数値計算処理手段による解析結果データおよび地図デー
タを格納するデータ格納手段と、前記数値計算処理手段
による計算結果と前記汚染物質の実測値データに基づい
て少なくとも発生源の位置を含む発生源情報を推定する
発生源推定手段と、前記数値計算処理手段による計算結
果、前記地図データによる地図および発生源位置を表示
装置に表示する表示処理手段とを具備したことを特徴と
する大気拡散シミュレーションシステム。
1. An input means for inputting simulation condition data for simulating a diffusion distribution state of pollutants emitted from a source of air pollutants, and respective data of analysis examples, laws and guidelines, simulation models and measured values. And a model selecting means for inputting the simulation condition data and selecting a simulation model suitable for the simulation, fixed parameters and analysis conditions from the database, and inputting the simulation condition data to obtain a numerical value of a diffusion distribution state. A parameter tuning means for performing a simulation and selecting an adjustment parameter such that an error from a measured value of the database storing actual measured values in a similar environment is reduced; and a model selected by the model selecting means. Numerical calculation processing means for performing a numerical simulation of the diffusion distribution state of the contaminant using a simulation model, fixed parameters, analysis conditions and numerical calculation data obtained based on the adjustment parameters selected by the parameter tuning means, Data storage means for storing analysis result data and map data by the numerical calculation processing means; and source information including at least the location of the source based on the calculation results by the numerical calculation processing means and the actually measured value data of the pollutant. An atmospheric dispersion simulation system comprising: a source estimating unit for estimating; and a display processing unit for displaying a calculation result by the numerical calculation processing unit, a map based on the map data, and a source position on a display device.
【請求項2】請求項1において、前記発生源推定手段
は、前記入力手段により汚染物質発生源の位置、放出時
間、放出量を入力し、前記数値計算処理手段を用いて前
記汚染物質の拡散分布状況を数値シミュレーションした
計算結果を汚染物質発生源の位置をパラメータとして複
数回計算し、汚染物質濃度実測値と合わせて統計的手法
を用いたデータ解析を行い汚染物質発生源位置を推定す
ることを特徴とする大気環境シミュレーションシステ
ム。
2. A method according to claim 1, wherein said source estimating means inputs the position, release time, and release amount of the pollutant source by said input means, and diffuses said contaminant using said numerical calculation processing means. Estimate the location of pollutant sources by calculating the results of a numerical simulation of the distribution situation multiple times using the location of the pollutant source as a parameter, and using the statistical method to analyze the data together with the actual measured values of the pollutant concentration. An atmospheric environment simulation system characterized by the following.
JP2001002458A 2001-01-10 2001-01-10 Atmospheric environment simulation system Pending JP2002207002A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2001002458A JP2002207002A (en) 2001-01-10 2001-01-10 Atmospheric environment simulation system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2001002458A JP2002207002A (en) 2001-01-10 2001-01-10 Atmospheric environment simulation system

Publications (1)

Publication Number Publication Date
JP2002207002A true JP2002207002A (en) 2002-07-26

Family

ID=18870962

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2001002458A Pending JP2002207002A (en) 2001-01-10 2001-01-10 Atmospheric environment simulation system

Country Status (1)

Country Link
JP (1) JP2002207002A (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100851520B1 (en) * 2005-12-28 2008-08-11 한국건설기술연구원 Apparatus and method for inferring atmosphere environment of underground space using statistical analysis
JP2013185988A (en) * 2012-03-08 2013-09-19 Toshiba Corp Information processing device and emission point estimation method
WO2014020747A1 (en) * 2012-08-02 2014-02-06 株式会社日立製作所 Physical quantity estimation system and physical quantity estimation method
JP2015172494A (en) * 2014-03-11 2015-10-01 新日鐵住金株式会社 Environment load molecule source evaluation method, environment load molecule source evaluation system, and computer program
US9326969B2 (en) 2012-09-11 2016-05-03 Norgine Bv Compositions
CN111443015A (en) * 2020-03-04 2020-07-24 平安国际智慧城市科技股份有限公司 Atmospheric pollutant tracing method and device, computer equipment and storage medium
CN112986072A (en) * 2021-02-09 2021-06-18 上海英凡环保科技有限公司 Pollution source odor online monitoring system and method

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100851520B1 (en) * 2005-12-28 2008-08-11 한국건설기술연구원 Apparatus and method for inferring atmosphere environment of underground space using statistical analysis
JP2013185988A (en) * 2012-03-08 2013-09-19 Toshiba Corp Information processing device and emission point estimation method
WO2014020747A1 (en) * 2012-08-02 2014-02-06 株式会社日立製作所 Physical quantity estimation system and physical quantity estimation method
JP5864753B2 (en) * 2012-08-02 2016-02-17 株式会社日立製作所 Physical quantity estimation system and physical quantity estimation method
JPWO2014020747A1 (en) * 2012-08-02 2016-07-11 株式会社日立製作所 Physical quantity estimation system and physical quantity estimation method
US9326969B2 (en) 2012-09-11 2016-05-03 Norgine Bv Compositions
JP2015172494A (en) * 2014-03-11 2015-10-01 新日鐵住金株式会社 Environment load molecule source evaluation method, environment load molecule source evaluation system, and computer program
CN111443015A (en) * 2020-03-04 2020-07-24 平安国际智慧城市科技股份有限公司 Atmospheric pollutant tracing method and device, computer equipment and storage medium
CN112986072A (en) * 2021-02-09 2021-06-18 上海英凡环保科技有限公司 Pollution source odor online monitoring system and method
CN112986072B (en) * 2021-02-09 2022-06-03 上海英凡环保科技有限公司 Pollution source odor online monitoring system and method

Similar Documents

Publication Publication Date Title
US8903646B2 (en) Method and determining system for automatically determining emission locations, and method and traffic control system based thereon for immission-dependent traffic control
CN110348746B (en) Air quality influence assessment method and device based on single pollution source
JP4206024B2 (en) Atmospheric substance index distribution analyzer
JP2002207002A (en) Atmospheric environment simulation system
WO2017221856A1 (en) Analysis device, analysis method, and recording medium
US7487070B2 (en) Method for planning a security array of sensor units
CN107870858A (en) Version contrast test method and system
KR102326602B1 (en) Fine spatial resolution map generation method for air pollutants using data assimilation and hybrid models and iformatio provison method
JP2000194742A (en) Ambient air quality simulation system
Mileu et al. Development of a QGIS plugin to dasymetric mapping
Stanislawski et al. Estimating positional accuracy of data layers within a GIS through error propagation
JP4141537B2 (en) Architectural CAD system
CN116702477B (en) Leakage source positioning method and system based on atmospheric diffusion model
JP2006259940A (en) Statistical information presentation device
Kedrin et al. System for Visualization the Map of Air Pollution by Vehicles
Williams et al. Pollution source direction identification: embedding dispersion models to solve an inverse problem
Jakobsson Data quality and quality management–examples of quality evaluation procedures and quality management in European national mapping agencies
Crenna An introduction to WindTrax
JP2013065208A (en) Substance discharge estimation device, method for the same and program
CN111461163B (en) Urban interior PM2.5 concentration simulation and population exposure evaluation method and device
JP3940834B2 (en) Analysis work support CAE device using numerical map information data
CN104424515A (en) Novel BeiDou technology and GIS technology-based tour group real-time management system
Böttinger et al. Visualization of 2D Uncertainty in Decadal Climate Predictions.
JP2013195170A (en) Air flow state prediction device, air flow state prediction method, diffusion state prediction device and diffusion state prediction method
JP2012132824A (en) Source estimation device and source estimation method for diffusion material

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 Written amendment

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 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20050517

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20051108