産業大気汚染源 大気拡散モデリング とは、大気汚染物質が 周囲の大気中 でどのように拡散するかを数学的にシミュレーションすること です。これは、汚染物質の拡散を支配する数式を解くアルゴリズム を含むコンピュータプログラムを用いて行われます。拡散モデルは 、工業プラント、自動車交通、偶発的な化学物質放出などの発生源から排出される大気汚染物質や毒素の風下側の大気中濃度 を推定するために使用されます。また、特定のシナリオ(例えば、排出源の変化)の下での将来の濃度を予測するためにも使用できます。そのため、大気質政策の策定において最も広く用いられているモデルです。拡散モデルは、広範囲に拡散し、大気中で反応する可能性のある汚染物質に最も有効です。空間的・時間的に非常に変動が大きい汚染物質(例えば、ブラックカーボン のように発生源からの距離による減衰が非常に急峻な汚染物質)や疫学 研究においては、統計的な土地利用回帰モデルも使用されます。
拡散モデルは、大気質の 保護と管理を担当する政府機関にとって重要です。これらのモデルは通常、米国およびその他の国の既存のまたは計画中の新しい工業施設が 国家大気質基準 (NAAQS)に準拠しているかどうか、または準拠する予定であるかどうかを判断するために使用されます。これらのモデルはまた、有害な大気汚染物質の排出を 削減するための効果的な制御戦略の設計にも役立ちます。1960 年代後半、米国 EPA の大気汚染管理局は、都市計画担当者や交通計画担当者が使用するモデルの開発につながる研究プロジェクトを開始しました。[ 1 ] このような研究から生まれた道路拡散モデルの主要かつ重要な応用例は、 1971 年にカナダのスパディナ高速道路 に適用されました。
大気拡散モデルは、公共安全対応者や緊急事態管理担当者が偶発的な化学物質放出の緊急計画策定にも使用されます。モデルは、有害物質や有毒物質の偶発的な放出による影響を判断するために使用されます。偶発的な放出は、化学物質や放射性核種などの有害物質が関わる火災、流出、爆発を引き起こす可能性があります。最悪のケースの偶発的放出源の用語 と気象条件を使用した拡散モデリングの結果は、影響を受ける場所、環境濃度を推定し、放出が発生した場合に適切な保護措置を決定するために使用できます。適切な保護措置には、風下方向の人々の避難や屋内待避 が含まれる場合があります。産業施設では、この種の影響評価または緊急計画策定は、米国連邦規則集第40巻第68部に規定されている米国 大気浄化法 (CAA)によって義務付けられています。
分散モデルは、モデル開発に使用される数学的手法によって異なりますが、いずれも以下のようなデータの入力が必要です。
風速や風向、大気 乱流 の量( 「安定度クラス」 と呼ばれるもので特徴づけられる)、外気温、上空に存在する可能性のある逆転 層の底までの高さ、雲量、日射量などの気象条件。発生源項(排出または偶発的放出の発生源項 における毒素の濃度または量)および物質の温度排出または放出パラメータには、発生源の位置と高さ、発生源の種類(火災、プール、または換気煙突など)、排出速度 、排出温度、質量流量 または放出速度などが含まれます。 発生源地点および近隣の住宅、学校、企業、病院などの受波地点における地形標高。 放出されたガスプルームの経路にある障害物(建物やその他の構造物など)の位置、高さ、幅、表面の粗さ、またはより一般的なパラメータである「農村部」または「都市部」の地形の使用。 最新の高度な拡散モデリングプログラムの多くは、気象データなどの入力を行うための前処理 モジュールを備えており、出力データをグラフ化したり、大気汚染物質の影響を受けた地域を地図上にプロットし たりするための後処理モジュールも備えています。影響を受けた地域のプロットには、濃度が最小から高濃度までの範囲を示す等値線が 含まれる場合もあり、これにより健康リスクが最も高い地域を特定できます。等値線プロットは、一般市民や対応要員に対する保護措置を決定する際に役立ちます。
大気拡散モデルは、大気拡散モデル、空気拡散モデル、大気質モデル、大気汚染拡散モデルとも呼ばれます。
大気層 大気中の汚染物質が大気中でどのように拡散するかを理解するには、地球の大気 層について議論する必要があります。地球の表面に最も近い層は対流圏 と呼ばれています。対流圏は海面から約18 km(11 マイル) の高さまで広がり、大気全体の質量の約80パーセントを占めています。成層圏は 次の層で、18 km(11 マイル) から約50 km(31 マイル) まで広がっています。3番目の層は中間圏で 、 50 km(31 マイル) から約80 km(50 マイル) まで広がっています。80 kmより上にも他の層がありますが、大気拡散モデルに関しては重要ではありません。
対流圏の最下層は惑星境界層 (PBL)または大気境界層 と呼ばれます。PBLの気温は高度の上昇とともに低下し、キャッピングインバージョン と呼ばれる逆転層に達します。キャッピングインバージョンとは、暖かい空気が冷たい空気よりも大気中で高い位置にある逆転層 の一種です。PBLのキャッピングインバージョンより下の領域を対流惑星境界層 と呼び、その高さは通常1.5~2 km(0.93~1.24 マイル) です。対流圏の上部(つまり、逆転層より上)は自由対流圏と呼ばれ、 対流圏界面 (地球大気における対流圏と成層圏の境界)まで広がっています。熱帯および中緯度地域では、日中の自由対流層は 対流圏全体を占めることがあり、熱帯収束帯 ではその高さは最大で10~18 km (6.2~11.2 マイル) になります。
大気境界層(PBL)は、大気汚染物質の輸送と拡散において重要な役割を果たします。これは、風の乱流力が地表付近で最も強いためです。PBLのうち、地表と逆転層の底の間にある部分は混合層と呼ばれます。大気中に放出された大気汚染物質のほぼすべては、この混合層内で輸送・拡散されます。一部の汚染物質は逆転層を貫通し、PBLより上の自由対流圏へと流れ込みます。
要約すると、地表から上空にかけての地球大気の層は、混合層と逆転層からなる境界層(PBL)、自由対流圏、成層圏、中間圏などです。多くの大気拡散モデルは、主に大気境界層(ABL)内での大気汚染物質の拡散をモデル化するため、境界層モデル と呼ばれています。混乱を避けるため、メソスケールモデル と呼ばれるモデルは、水平方向に数百キロメートルまで拡散をモデル化する機能を備えています。これは、中間圏での拡散をモデル化するという意味ではありません。
ガウス型大気汚染物質拡散方程式 大気汚染拡散に関する技術文献は非常に豊富で、1930年代以前にまで遡ります。初期の大気汚染プルーム拡散方程式の1つは、BosanquetとPearsonによって導出されました。[ 2 ] 彼らの方程式はガウス分布を 仮定しておらず、汚染プルームの地表反射の影響も含まれていませんでした。
サー・グラハム・サットンは1947年に大気汚染物質のプルーム拡散方程式を導出した[ 3 ]。 この方程式には、プルームの垂直方向および横風方向の拡散に対するガウス分布の仮定が含まれており、またプルームの地表反射の影響も含まれていた。
厳格な環境規制 の出現によってもたらされた刺激により、1960 年代後半から今日までの間に大気汚染物質プルーム拡散計算の使用が飛躍的に増加しました。その期間中に、大気汚染物質排出の拡散を計算するための多数のコンピュータ プログラムが開発され、それらは「大気拡散モデル」と呼ばれました。これらのモデルのほとんどの基礎は、以下に示す連続浮力大気汚染プルーム のガウス拡散モデリングの完全な方程式です。 [ 4 ] [ 5 ]
C = Q u ⋅ f σ y 2 π ⋅ g 1 + g 2 + g 3 σ z 2 π {\displaystyle C={\frac {\;Q}{u}}\cdot {\frac {\;f}{\sigma _{y}{\sqrt {2\pi }}}}\;\cdot {\frac {\;g_{1}+g_{2}+g_{3}}{\sigma _{z}{\sqrt {2\pi }}}}}
上記の式には、地面からの上向き反射だけでなく、大気中に存在する逆転層の底からの下向き反射も含まれています。
4 つの指数項の合計g 3 {\displaystyle g_{3}} 最終値に非常に速やかに収束します。ほとんどの場合、 m = 1、m = 2、m = 3の級数の和で十分な解が得られます。
σ z \displaystyle \sigma _{z}} そしてσ y \displaystyle \sigma _{y}} これらは、大気安定度クラス(すなわち、周囲の大気の乱流の度合いを示す指標)と、受容体までの風下距離の関数である。汚染物質排出拡散度合いに影響を与える最も重要な2つの変数は、排出源の高さと大気乱流の度合いである。乱流が大きいほど、拡散度合いは高くなる。
式[ 6 ] [ 7 ] σ y \displaystyle \sigma _{y}} そしてσ z \displaystyle \sigma _{z}} は:
σ y \displaystyle \sigma _{y}} (x) = exp(I y + J y ln(x) + K y [ln(x)] 2 )
σ z \displaystyle \sigma _{z}} (x) = exp(I z + J z ln(x) + K z [ln(x)] 2 )
(単位)σ z \displaystyle \sigma _{z}} 、 そしてσ y \displaystyle \sigma _{y}} (xはメートル単位)
安定性クラスの分類はF. Pasquillによって提案された。[ 8 ] 6つの安定性クラスは以下の通りである:A-極めて不安定、B-やや不安定、C-やや不安定、D-中立、E-やや安定、F-やや安定
大気汚染物質濃度 の計算結果は、調査対象地域における汚染物質濃度の空間的変動を示すために、しばしば大気汚染物質 濃度等高線図 として表現される。このようにして等高線は、感受性の高い受容体の位置と重なり合い、大気汚染物質と対象地域との空間的な関係を明らかにすることができる。
従来のモデルでは、 σ y \displaystyle \sigma _{y}} そしてσ z \displaystyle \sigma _{z}} より最近のモデルでは、これらのパラメータを導出するために、モニン・オブホフの相似理論 にますます依存するようになっている。
ブリッグスの煙の上昇方程式 (上記で説明した)ガウス大気汚染物質拡散方程式では、汚染物質プルームの中心線の高さであるHを入力として必要とします。Hは、 H s (汚染物質プルームの排出源の実際の物理的な高さ)とΔ H (プルームの浮力によるプルームの上昇)の合計です。
浮力を持つガウス型大気汚染物質拡散プルームの可視化 Δ H を決定するために、1960 年代後半から 2000 年代初頭にかけて開発された大気拡散モデルの多く、あるいはほとんどは、ブリッグス方程式として知られるものを使用していた。GA ブリッグスは、1965 年に初めて煙の上昇に関する観測と比較を発表した。[ 9 ] 1968 年、CONCAWE (オランダの組織) が主催するシンポジウムで、彼は当時文献で入手可能だった多くの煙の上昇モデルを比較した。[ 10 ] 同年、ブリッグスは、スレードが編集した出版物[ 11 ] の、煙の上昇モデルの比較分析に関するセクションも執筆した。1969 年に、彼は煙の上昇に関する文献全体の古典的な批判的レビューを発表し[ 12 ] 、その中で「ブリッグス方程式」として広く知られるようになった一連の煙の上昇方程式を提案した。その後、ブリッグスは1971年と1972年に1969年の煙の上昇方程式を修正した。[ 13 ] [ 14 ]
ブリッグスは、大気汚染物質のプルームを以下の4つの一般的なカテゴリーに分類した。
穏やかな外気条件下における冷たいジェット噴流 風の強い外気条件下における冷たいジェット噴流 穏やかな周囲の空気中で発生する、高温で浮力のある噴煙 風の強い周囲の空気中で発生する、高温で浮力のある噴煙 ブリッグスは、冷たいジェットプルームの軌跡は初期速度運動量によって支配され、高温で浮力のあるプルームの軌跡は浮力運動量によって支配され、初期速度運動量は比較的重要ではないと考えた。ブリッグスは上記の各プルームカテゴリーについてプルーム上昇方程式を提案したが、広く使用されるようになった「ブリッグスの方程式」は、彼が曲がった高温の浮力のあるプルームに対して提案したものであることを強調することが重要である。
一般的に、ブリッグスの、曲がった高温の浮力のある煙柱に関する方程式は、大規模発電所で化石燃料 を燃焼させる蒸気発生ボイラーの煙突など、典型的な燃焼源からの煙 柱に関する観測とデータに基づいています。したがって、煙突出口の速度はおそらく20~100フィート/秒(6~30メートル/秒)の範囲で、出口温度は250~500 °F(120~260 °C)の範囲であったと考えられます。
ブリッグス方程式[ 4 ] を用いて、曲がった浮力のある噴煙の上昇軌跡を求めるための論理図を以下に示します。
ブリッグスの方程式で使用される上記のパラメータについては、ベイチョクの著書で説明されています。[ 4 ]
参考文献 ↑ Fensterstock, JC 他、「環境計画による大気汚染ポテンシャルの低減」、JAPCA、第21巻、第7号、1971年。 ↑ Bosanquet, CH および Pearson, JL、「煙突からの煙とガスの拡散」、Trans. Faraday Soc.、32:1249、1936 ↑ サットン、OG、「下層大気における拡散の問題」、QJRMS、73:257、1947年、および「工場の煙突からの大気汚染の理論的分布」、QJRMS、73:426、1947年 1 2 3 ベイチョック、ミルトン R. (2005).スタックガス分散の基礎 (第 4 版). 著者出版. ISBN 0-9644588-0-2 。 ↑ Turner, DB (1994). 大気拡散推定ワークブック:拡散モデリング入門 (第2 版). CRC Press. ISBN 1-56670-023-X 。↑ サインフェルド、ジョン・H. (2006). 「第18章」. 大気化学と物理学:大気汚染から気候変動まで . ワイリー. ISBN 9780471720171 。↑ ハンナ、スティーブン (1982)。 「大気拡散ハンドブック」 。 米国エネルギー省報告書 。 ↑ W、クルーグ(1984年4月)。 大気拡散(第3版)。F. パスクイル、FB スミス。エリス・ホーウッド、(ジョン・ワイリー・アンド・サンズ)チチェスター、1983年 (第3 版)。ニューヨーク:王立気象学会季刊誌。 ↑ Briggs, GA、「観測結果と比較した煙の上昇モデル」、JAPCA、15:433–438、1965年 ↑ Briggs, GA、「CONCAWE会議:異なる煙の上昇式による影響の比較に関する議論」、Atmos. Envir.、2:228–232、1968年 ↑ スレード、DH(編)「気象学と原子力 1968」、米国商務省大気資源研究所、1968年 ↑ ブリッグス、GA、「煙の上昇」、米国原子力委員会批評シリーズ、1969年 ↑ Briggs, GA、「煙の上昇観測に関する最近の分析」、第2回国際大気浄化会議議事録、Academic Press、ニューヨーク、1971年 ↑ Briggs, GA、「考察:中立かつ安定した環境における煙突の煙」、大気環境、6:507–510、1972年
さらに読む
本 入門編 ベイチョック、ミルトン R. (2005).煙突ガス分散の基礎 (第4 版). 著者出版. ISBN 0-9644588-0-2 。 化学プロセス安全センター(1999)。化学プロセス定量的リスク分析ガイドライン (第2 版)。米国化学工学会、ニューヨーク州ニューヨーク。ISBN 978-0-8169-0720-5 。 化学プロセス安全センター(1996)。蒸気雲および発生源拡散モデルの使用に関するガイドライン、演習例付き (第2 版)。米国化学工学会、ニューヨーク州ニューヨーク。ISBN 978-0-8169-0702-1 。 Schnelle, Karl B. & Dey, Partha R. (1999).大気拡散モデリング準拠ガイド (第1 版). McGraw-Hill Professional. ISBN 0-07-058059-6 。 Turner, DB (1994).大気拡散推定ワークブック:拡散モデリング入門 (第2 版). CRC Press. ISBN 1-56670-023-X 。 高度な Arya, S. Pal (1998).大気汚染の気象学と拡散 (第1 版). オックスフォード大学出版局. ISBN 0-19-507398-3 。 ロッド、バラット (2001)。大気分散モデリング (第 1 版)。アーススキャン出版。ISBN 1-85383-642-7 。 コルズ、ジェレミー(2002)。大気汚染 (第1 版)。スポン・プレス(英国)。ISBN 0-415-25565-1 。 Cooper JR、Randle K、Sokh RG (2003)。環境における放射性物質の放出 ( 第1 版)。John Wiley & Sons。ISBN 0-471-89924-0 。 欧州プロセス安全センター(1999)。大気拡散 (第1 版)。ラグビー:化学工学会。ISBN 0-85295-404-2 。 ゴディッシュ、サッド(2003)。大気質 (第4 版)。CRCプレス。ISBN 1-56670-586-X 。 Hanna, SR & Drivas, DG (1996).水蒸気雲拡散モデルの使用に関するガイドライン (第2 版). Wiley-American Institute of Chemical Engineers. ISBN 0-8169-0702-1 。 Hanna, SR & Strimaitis, DG (1989).蒸気雲源拡散モデルのテストケースワークブック (第1 版)。米国化学工学会化学プロセス安全センター。ISBN 0-8169-0455-3 。 Hanna , SR & Britter, RE (2002).工業地帯および都市部における風の流れと水蒸気雲の拡散 (第1 版)。Wiley-American Institute of Chemical Engineers。ISBN 0-8169-0863-X 。Perianez, Raul (2005).海洋環境における放射性核種の拡散モデリング :入門 (第1 版). Springer. ISBN 3-540-24875-7 。 ピールケ、ロジャー A. (2001).メソスケールモデリング (第2 版). エルゼビア. ISBN 0-12-554766-8 。 Zannetti, P. (1990).大気汚染モデリング :理論、計算方法、および利用可能なソフトウェア 。Van Nostrand Reinhold。ISBN 0-442-30805-1 。 Forago I、Georgiev K、Havasi A 編 (2004)。環境安全保障のための大気汚染モデリングの進歩 (NATO ワークショップ) 。Springer、2005 年。ISSN 0957-4352。 Kretzschmar JG、Cosemans G 編 (1996)。規制目的のための大気拡散モデリングにおける調和 (第 4 回ワークショップ) 。International Journal of Environment and Pollution、第 8 巻、第 3-6 号、Interscience Enterprises、1997 年。ISSN 0957-4352 。 Bartzis, J G. 編 (1998).規制目的のための大気拡散モデリングにおける調和 (第 5 回ワークショップ) . International Journal of Environment and Pollution、第 14 巻、第 1-6 号、Interscience Enterprises、2000 年。ISSN 0957-4352 . Coppalle, A. 編 (1999).規制目的のための大気拡散モデリングにおける調和 (第 6 回ワークショップ) . International Journal of Environment and Pollution、第 16 巻、第 1-6 号、Inderscience Enterprises、2001 年。ISSN 0957-4352 。 Batchvarova, E. 編 (2002).規制目的のための大気拡散モデリングにおける調和 (第 8 回ワークショップ) . International Journal of Environment and Pollution, vol. 20 no. 1–6, Inderscience Enterprises, 2003. ISSN 0957-4352 . Suppan, P. 編 (2004).規制目的のための大気拡散モデリングにおける調和 (第 8 回ワークショップ) . International Journal of Environment and Pollution、第 24 巻、第 1-6 号および第 25 巻、第 1-6 号、Inderscience Enterprises、2005 年。ISSN 0957-4352 。 Zannetti, P. 編 (1993).国際大気汚染会議 (第 1 回、メキシコシティ) . Computational Mechanics、1993 年。ISBN 1-56252-146-2 。 De Wispelaere, C. 編 (1980).大気汚染モデリングとその応用に関する国際技術会議 (第 11 回) . Plenum Press、1981 年。ISBN 0-306-40820-1 。 De Wispelaere, C. 編 (1982).大気汚染モデリングとその応用に関する国際技術会議 (第 13 回) . NATO 現代社会の課題に関する委員会 [by] Plenum Press、1984 年。ISBN 0-306-41491-0 。 Gryning, S.、Schiermeir, FA 編 (1995).大気汚染モデリングとその応用に関する国際技術会議 (第 21 回) . NATO 現代社会の課題に関する委員会 [by] Plenum Press、1996 年。ISBN 0-306-45381-9 。 Gryning, S.、Chaumerliac, N. 編 (1997).大気汚染モデリングとその応用に関する国際技術会議 (第 22 回) . NATO 現代社会の課題に関する委員会 [by] Plenum Press、1998 年。ISBN 0-306-45821-7 。 Gryning, S.、Batchvarova, E. 編 (1998).大気汚染モデリングとその応用に関する国際技術会議 (第 23 回) . NATO 現代社会の課題に関する委員会 [by] Kluwer Academic/Plenum Press、2000 年。ISBN 0-306-46188-9 。 Gryning, S.、Schiermeir, FA 編 (2000)。大気汚染モデリングとその応用に関する国際技術会議 (第 24 回) 。NATO 現代社会の課題に関する委員会 [by] Kluwer Academic、2001 年。ISBN 0-306-46534-5 。 Borrego, C.、Schayes, G. 編 (2000)。大気汚染モデリングとその応用に関する国際技術会議 (第 25 回) 。NATO 現代社会の課題に関する委員会 [by] Kluwer Academic、2002 年。ISBN 0-306-47294-5 。 Borrego, C.、Incecik, S. 編 (2003).大気汚染モデリングとその応用に関する国際技術会議 (第 26 回) . NATO 現代社会の課題に関する委員会 [by] Kluwer Academic/Plenum Press、2004 年。ISBN 0-306-48464-1 。 米国科学アカデミー出版局、有害物質放出の大気拡散に関する委員会編(2002年)。有害物質放出の大気拡散の追跡と予測(ワークショップ) 。2003年。ISBN 0-309-08926-3 。
ガイダンス Hanna, SR; Briggs, GA & Hosker, RP (1982).大気拡散ハンドブック 。米国エネルギー省技術情報センター。DOE/TIC- 11223。OSTI 5591108 。 米国環境保護庁(1993年)。有害物質/有毒物質の大気放出に対する改良拡散モデルの適用に関するガイダンス 。大気質計画基準局、EPA-454/R-93-002。 米国環境保護庁(1999)。オフサイト影響分析に関するリスク管理プログラムガイダンス(付録) (PDF) 。固体廃棄物・緊急対応局、EPA 550-B-99-009。2010年4月17日にオリジナル(PDF)からアーカイブ。 2010年4月9日 に取得 。 米国環境保護庁(1999年)。無水アンモニア、水酸化アンモニア、塩素、二酸化硫黄のオフサイト影響分析に関する技術的背景文書 (PDF) 。化学物質緊急事態対策・予防局。2008年3月28日に オリジナル(PDF) からアーカイブ。 米国環境保護庁(2009年)。第4章:オフサイト影響分析。化学事故防止のためのリスク管理プログラムに関する一般ガイダンス(40 CFR Part 68) (PDF) 。固体廃棄物・緊急対応局、EPA 555-B-04-001。2008年3月28日にオリジナル(PDF)からアーカイブ。
外部リンク EPA(米国環境保護庁)の大気規制モデリング支援センター 米国環境保護庁(EPA)の大気質モデリンググループ(AQMG) NOAA(米国海洋大気庁)の大気資源研究所(ARL) 英国大気拡散モデリング連絡委員会ウェブサイト 英国分散モデリング局のウェブサイト 大気化学輸送モデル LOTOS-EUROS 運用優先物質モデル(OPS )(オランダ語) HAMS-GPS分散モデリング 大気拡散モデリングに関するWiki。大気拡散モデリングの国際的なコミュニティ(主に研究者、そしてモデルの利用者)を対象としています。その目的は、拡散モデリング担当者が業務を通じて得た経験を集約することです。