福島原発事故シミュレーション
amazon：Excelで操る! ここまでできる科学技術計算 第2版 楽天： Excelで操る! ここまでできる科学技術計算 [ 神足 史人 ] |
図3 放射線物質の大気拡散分布図(ガウス分布)
計算方法
| 大気汚染計算は、通常、汚染物質の排出量を入力とし、拡散方程式を適用して地上の汚染物質濃度を求めるものであるが、ここでは逆に各地の放射能レベルを入力とし、拡散方程式を逆に適用することによって、1日当たりに放出されている放射線物質の量を推定する。 |
観測された各地の放射能レベル
放射線レベルの連続的監視が行われている地点のうち、福島第一原子力発電所から約４０ｋｍの距離にある７地点を選んで、3月24日、25日の放射線レベルを示すと、図１になる。汚染は放出点から一様に広がっているのではなく、風道にあたるところに高い濃度が現れる。飯舘村はそのような風道の上にある。原発と飯舘村を結ぶ風道の直線から１０ｋｍ程度しか離れていない南相馬では飯舘村に比べての濃度が極端に下がっていることからみて、40km地点での風道の幅は４〜６ｋｍと見られる。 |
シーベルトからベクレルに換算
飯舘村の放射線影響強度は１０μSv/hである。これを大気中の放射線濃度(Bq/m3)に換算する。 計算結果は130Bq/m3となる。  図３ Excelを用いて、シーベルトをベクレルに換算 |
適用する拡散モデルはプルームモデル
拡散モデルにはパフモデルとプルーフモデルがあるが、図１の汚染レベルの分布から見て、風速２〜５ｍ/secのプルームモデルを適応する。(文献１）プルームとは煙のことである。プルームモデルでは、煙は風下の方向（ｘ軸）に向かって幅（ｙ軸）と高さ（ｚ軸）を増しながら広がる。風速と気象状況により表１ 日本式パスキル安定度階級分類法に基づき安定度階級（A〜Ｇ）を決定し、図２のグラフより、風下距離(x)の位置のプルーフの幅（σy）と高さ（σZ）を求める。プルーフ内の濃度は(1)式のガウス分布で近似できる。(図３参照) exp(-(z/σz)2-(y/σy)2) ・・・・・・ (1) 放出量(Q)は、プルームの高さ２σｚと幅２σｙと風速vと平均濃度Ｃavをかけて求められる。 Q = 4σz・σy・v・Cav ・・・・・・・ (2) 平均濃度はプルームの中心から高さσzと幅σyのところに現れる。 ３月23日、24日の気象条件は風速２〜３m/s、曇り小雨で、安定度はＤとみなせる。図２より、 σz=200m σy=1000m とする。 |
観測値から放出量の推定
まず、排出煙元の高さＨを推定する必要があるが、原発から立ち昇る水蒸気の湯気から判断すると、原発建屋の5〜6になる。原発建屋の高さは約５０ｍなのでＨは300mと推定する。そして、飯舘村の300m上空をプルームの中心が通過したと仮定する。プルーム濃度がガウス分布することにより、地上の計測値から平均濃度を算出することが出来る。平均濃度はプルーム中心から下にσz(200m)のところの濃度である。Ｅｘｃｅｌで1日当たりの放出量を算出すると、６３ＴＢｑ/dとなる。(T(テラ）は1012)  図４ Ｅｘｃｅｌを用いて、観測値から放出量の推定 |
参考文献１ 環境科学辞典
参考文献２ 現代化学 2011年5月号
理論物理計算が示す原発事故の真相
放射性物質の海域への放出量計算
EXCELを用いた科学技術計算
EXCELのソルバーを用いた科学技術計算
EXCELを用いた構造力学(材料力学)計算
〒６７３−００３６ 兵庫県明石市松江６２−１４
(有)ゴッドフット企画
感想をお聞かせ下さい
Copyright (C) GODFOOT 2013. All Rights Reserved.