RISTニュース

No．67（2021）

環境中での空間線量率 3 次元分布計算システム（3D-ADRES）の研究開発の現状と今後の展望－福島復興に向けた市街地・森林等の環境中空間線量率分布の推定－

Present Status and Future Perspective of R&D for 3D-ADRES(3-Dimensional Air Dose Rate Evaluation System)- Evaluation Scheme for Air Dose Rate Distributions inCity and Forest Areas towards Fukushima’s revitalization 金敏植、Alex Malins、町田

日本原子力研究開発機構
昌彦、吉村 和也、斎藤 公明
東北大学大学院薬学研究科
吉田 浩子

柳

高度情報科学技術研究機構
秀明、吉田 亨、長谷川 幸弘

キーワード：福島の復興、空間線量率 3 次元分布、福島第一原子力発電所事故、放射性セシウ
ム、3D-ADRES、PHITS、家屋、市街地、森林
福島第一原子力発電所事故により環境中に放出された放射性セシウムの一部は陸域に降下した。放射性セシウムは土壌表層に留まる性質を有するため、土壌表面は主たる放射線源となった
他、道路、家屋の壁・屋根、森林内の樹木樹皮等にも一部残存し、市街地や森林等環境中の放射線量は複雑な分布を示すことが知られている。その上、それらの放射線源から発する放射線は複雑な環境構造物に吸収・散乱されるため、放射線量の詳細な分布を知ることは容易ではない。しかし、その環境中にて生活する住民にとって、大半の時間を過ごす家屋及びその近縁、職場等の生活空間における空間線量率の分布は重要な関心事である。本報告では、この環境中での複雑な放射線量の 3 次元的分布を可能な限り高精度に推定するため、研究開発を行ってきた3D-ADRESの現状を報告すると同時に、3D-ADRESを活用することで得られる成果の一端を示し今後の展望についても記す。
1 ．はじめにに強く吸着される性質を有するため、降着量2011年に発生した東京電力ホールディンの多い地域では、高い放射線量が計測され、グス（株）・ 福 島 第 一 原 子 力 発 電 所（以下多くの住民が避難する事態となった。避難を1F）事故により、放射性セシウムが環境中余儀なくされた住民にとって、かつて居住しに放出され、その一部は陸地へと降着した。ていた地域の空間線量率の分布や経時変化は

放射性セシウムは、土壌表層の粘土鉱物粒子重大な関心事である他、放射線量が比較的低－3－

RISTニュース

No．67（2021）

く、居住を続ける住民にとっても、周囲の放率の 3 次元分布をシミュレーションするに射線量の分布とその変化は、生活上の被ばくは、先ず事故初期における放射性核種の沈着を自ら知り、管理するための欠かせない情報分布を把握する必要がある。次に、「環境面」源となる。と定義する放射性物質が沈着した環境中の表

1F事故後、国・自治体は、居住する住民の面において、放射性核種がどのような動きを被ばくを低減するため、各地で除染を行い、示し、空間線量率がどのように変化するかをその結果、放射線量は低減したが、その実施環境動態の観点から理解する必要がある。

は平地からなる市街地や家屋周辺等に限定さこの観点に対し、吉村らは 文献［ 4 ］、1F事れ、福島県の凡そ 7 割を占める山林の多くは故後初期に福島県内を広域でかつ高頻度に環除染未実施のまま残されている文献［ 1 ］。
従って、境モニタリングを実施した。その結果、空間除染実施の有無、降着量の違い、そして環境線量率を137Cs沈着量で除した正規化線量率

中での放射性セシウムの移動等により、空間の減少速度は放射性崩壊よりも速やかであ線量率の分布は不均一となり、場所により異り、その要因として、放射性核種の深度方向なる分布を示すことが知られている文献［ 2 ］。

への移行や水平方向への流失が考えられると一方、空間線量率の経時変化に目を向ける報告している。さらに、その減少速度は土地と、現在、事故後10年余りが経過し、環境利用形態により、異なると報告している。

中の空間線量率の減少は予想を上回り、特にまた、1F事故後、大規模にCsの分布を調

人間の生活に密接に関係する環境の空間線量査するため、国の方針として放射線量等分布率は顕著に速く減少している文献［ 3 ］。このよマップ 文献［ 5 ，6 ］ の作成が行われ、航空機にうな環境中での変化を理解するには、放射性よる上空のモニタリングの他、地上でのモセシウムの動態を十分に理解する必要がある。

ニタリングが実施された。原子力機構では、以上の背景の下、住民帰還や福島の復興に他研究機関との協力を得て1Fから半径80km向け、重要となる情報は、生活圏内の空間圏内及び福島県内にある土壌面に対して、線量率の分布及びその経時変化となる。こ2011年 6 月から数多くのモニタリングを行の情報を提供すべく、日本原子力研究開発

い、沈着土壌濃度マップを作成した（2011

機構（JAEA）と高度情報科学技術研究機構

年 6 月から約2200地点、2012年から2016年（RIST）は、最新の計算科学技術を活用し、にかけては、測定地点を約6500地点に増や

図 1 に示すように福島の現地にある家屋・森

し沈着量及び空間線量率のモニタリングを実

林環境・線源分布を可能な限りリアルにモデ施）。この結果を基に、三上ら文献［ 7 ］は、1Fル化した上で、そのモデルを用いて空間線量事故後の初期から中期にかけてのモニタリ率の 3 次元分布の推定を可能にするシステ

ング結果から、2016年の平均空間線量率は、ム（ 以 下、3D-ADRES（3-Dimensional Air

セシウム134（134Cs、半減期：2.06年）の崩Dose Rate Evaluation System）） の 研 究 開

壊と共に放射性セシウムが土壌に浸透する影発を進めてきた。本報告では、3D-ADRES

響で、初期の約30%程度の減少が見られたと研究開発の現状と、それを用いた成果の一端報告した。また、未除染地域における土壌でを報告すると同時に、今後の研究開発の展望は、放射性崩壊から予想される速度で134Csについても報告する。

及び137Csの平均沈着量は減少し、放射性セシウムの水平方向の移動は比較的小さい一2 ．福島での空間線量率の分布と経時変化方、除染等が実施された地域のモニタリング

1F事故後の福島における詳細な空間線量結果をみると、定量的な分析は難しいものの、－4－

RISTニュース

No．67（2021）

平均空間線量率は約20%程度減少し、除染等かにするため、放射線輸送シミュレーションの効果が明らかになったと報告している。こを行った。その結果、事故後初期の林内の空れに加え、中間ら 文献［ 8 ］ による長期モニタ

間線量率は、樹冠・リター層からの影響が比

リングデータの結果でも、空間線量率は放射較的大きい一方、2015年以降は、殆どが土性セシウムの放射性崩壊による減衰速度以上

壌の上部 5 ㎝以内に達した放射性物質によるに早く減少し、特に非舗装面より、アスファものと報告している。

ルト道路の舗装面の方が早く減少したと報告

以上、1F事故後の初期・中期・そして最している。一方、森林においては、減少速度近の放射性セシウムの沈着状況と、さらにはは小さく、市街地からの距離が近くなるほど、土地利用形態や除染、人間活動等が、その場減少速度は増大するとの報告があった。これの空間線量率の分布に強く影響することが分らの結果を踏まえ、Malinsら 文献［ 9 ］ は森林かる文献［10］。このような現状を踏まえ、放射内の空間線量率を決める要因について、放射性セシウムの複雑な分布による 3 次元の空間性物質の存在位置と空間線量率の関係を明ら線量率分布の定量的推定とその可視化を目指

図1

空間線量率 3 次元分布計算システム（3D-ADRES）の概要と計算の流れ。市街地や森林等において、地理空間情報データから地形・建物・樹木モデルを作成し、線源分布を設定した後、PHITSを用いたモンテカルロシミュレーションを行い、空間線量率 3 次元分布を評価する。
－5－

RISTニュース

No．67（2021）

し、私たちは、3D-ADRESをこれまでに開築し、放射線場を効率的にシミュレーション発してきた。 3 章では、上記の背景を踏まえ、可能とすることを目標に、3D-ADRESを開開発を進めてきた内容の中でも重要なポイン発してきた。以下では、既存手法の限界等のトについて、その概略を紹介する。

課題を挙げ、改良を重ねてきた点を中心に報告する。

3 ．3D-ADRESの開発コンセプトと開発の現状

3．
1．家屋を中心とするリアルな屋内・屋外環境空間を再現するモデルの開発

3D-ADRESでは、環境における空間線量率の分布を解析するため、放射線の挙動を核

前述のように、3D-ADRESの開発により、反応モデルや核データ等を用いて模擬するモ

環境中に存在する複雑な地形・家屋・樹木等

ンテカルロ計算コード：Particle and Heavy

が 3 次元モデルとして再現可能となった。しIon Transport code System（PHITS）文献［11］かし、「解析対象」や「解析目的」に応じて、を用いて計算する。PHITSは、現在、地面・線量率の分布を評価するには、分布の細部に家屋・樹木・線源分布等の複雑な環境を簡渡る再現性という観点で限界があった。例え

易にモデル化するGraphical User Interface

ば、図 2 に示すように、帰還住民等が最も長（GUI）を備えたツールを持たないため、環い時間滞在する家屋の室内空間における空間境モデルを作成するには大きな手間がかか線量率をシミュレーションするには、①対象る。この状況を克服するため、JAEA・シス

家屋の近場にある線源が重要となるため、近テム計算科学センターは、RISTと連携し、場の土地利用状況を細かく設定し、再現性JAEA・福島研究開発部門の支援の下、実環を高める必要がある文献［13］。②さらには、建境における空間線量率の 3 次元分布が計算可

物直下の土壌が汚染されていない効果文献［14］能なシステム3D-ADRES文献［12］ を開発して（家屋直下に線源がないことが重要な効果をきた。3D-ADRESを用いると、地形・樹木・もたらす、un-contaminated effect）が重要建物等のリアルな 3 次元モデルの構築と不均となるため、家屋の形状に合わせて家屋真下質な放射性セシウムの線源分布の取り込みが等の未汚染形状の再現性を高める必要があ可能となり、福島県における市街地や森林等る。しかしながら、①の場合、計算コスト等のような複雑な環境モデルを再現できる。の問題で既存の地面格子は10ｍが限界であ

しかし、市街地、森林・農地等では、土地り、比較的、線源密度の高い庭及び周囲の草利用が複雑に入り組んだ環境となり、それに地等の存在は表現されず、家屋内外の空間線応じて沈着の度合いも変化する一方、除染等量率分布を詳細にシミュレートするには適切を含めた人間活動等を考慮し、且つ、家屋内とは言えない限界があった。②の場合は、当への各環境面からの影響を評価するには、地初、採用した家屋の形状が四角形のみで、縮面モデルの解像度の向上、実家屋モデルの構尺だけを考慮し、縦横高さのスケール変換の築、家屋室内空間のモデル化、各環境面におみで対応したため、家屋の真下等の未汚染形ける線源分布の細かな設定というモデル構築状の設定に限界があった。また、住民にとっ上の限界を広げる必要があった。更には、数て、家屋内の空間線量率の分布は重要だが、多くの複雑な環境モデルが設定されると、計その空間線量率がどの屋外の「環境面」から算負荷が高まる等の課題があり、効率的計算影響を受けているかを詳細に知りたいというを行うための工夫も必要となっていた。

ニーズがあった。これは、測定結果だけから以上、より詳細でリアルな環境モデルを構推定することは難しく、「環境面」単位での－6－

RISTニュースシミュレーションが必要となる。

No．67（2021）

が10m程度であり、且つモデル作成上の簡単

以 上 の 課 題 や ニ ー ズ を 受 け、 現 在 のさのため、規則格子を用いると、実際の福島

3D-ADRESでは、環境モデルのリアリティーの家屋周辺の土地利用分布を詳細に再現するを向上させ、PHITSシミュレーションを通ことが困難となるケースが多いことが分かっして屋内空間における空間線量率の由来がどた。特に家屋周囲の軒下や駐車場のような整のような家屋周辺の放射性物質の分布に起因備された舗装面と家屋周囲の庭や畑等の非舗するかを明らかにするため、以下のような研装面が交わる境界部分や、家屋に繋がる細い究開発を行い改良を加えた。

私道のような舗装面の表現は難しく、線源分布を再現性よく設定するには限界があった。

3．1．1．不規則メッシュの開発

これらを考慮するため、新しい3D-ADRES

福島県における様々な測定結果を見ると、では、既存の規則メッシュに加えて不規則三空間線量率の低減効果は、土地利用により大角メッシュを開発した。尚、3D-ADRESのきく異なっている。例えば、舗装面の場合、GUIでは、直接、不規則三角メッシュを作成

非舗装面よりも比較的、速やかに減少するこ

できないが、GISプログラムを利用することとが確認されている文献［15］。従って、正確なで、その不規則三角メッシュを3D-ADRES

空間線量率分布を解析するには、地面メッに入力可能とすることで課題を解決した。図シュを土地利用毎に詳細に分け作成すること

3 は、規則三角メッシュと不規則三角メッが望ましい。しかし、通常利用されている

シュを用いて家屋周囲のメッシュを作成した

GISデータや既存の地面メッシュは、解像度場合の違いを示した。図 3 より、明らかに不図2

周囲に林が隣接する典型的な福島の家屋環境：屋内の中心位置（床上 1 ｍ）に対し、周囲の環境から来るガンマ線の挙動イメージ図－7－

RISTニュース

No．67（2021）

規則三角メッシュの方が規則三角メッシュときいかを調べることも可能となる（図 3 ）。

比べ、現実に沿う形で詳細に土地利用分布をこの機能については、3.1.3.シナリオ機能の再現できることが分かる。また、不規則三角開発で追記する。

メッシュの特徴として、メッシュ一つ一つの3．
1．
2．形状リアリティーを再現可能とする

サイズを変えることも可能であり、図 3 の家
屋周囲の細い私道等を凡そ正確にモデル化で建物モデルの開発きることも分かる。図 4 に不規則三角メッ①多様な建物モデルの形状の再現

シュを使用し、作成した家屋環境モデルの例

建物の形状は単純ではなく建築工法も様々を示す。さらに、不規則三角メッシュは設定である上に、家屋周囲の線源分布は不均質でする位置も自由に設定できるため、家屋からあるため、建物内外の空間線量率分布は複雑

予め距離・環境面毎に、不規則三角メッシュ

な 3 次元分布を示す。中でも、建物内の空間を作成することで、調査対象とする家屋の室線量率や建物周囲の空間線量率は、当然建物内において、どの環境面からの影響が最も大周りの放射性セシウムの沈着分布に強く依

図3

家屋（灰色で示したエリア）とその周囲を規則三角メッシュでモデル化した場合（a）と不規則三角メッシュを用いた場合（b）の違い

図4

不規則三角メッシュを使用した家屋環境モデルの作成例（現行3D-ADRESにより作成した環境モデルをPHITSにより出力）－8－

RISTニュース

No．67（2021）存する。松田らによると文献［14］、建物がその以上、上記のニーズの下、建物形状を四面場の放射線場に及ぼす影響として、建物によ体メッシュ体系で構成し、多様な形状の建物る遮蔽効果の他、建物直下が汚染されていなモデル作成を可能とした（図 5 ）。四面体メッい効果があると報告している。従って、建シュの活用により、図 5 で示したコ型、セッ

物直下が汚染されていない効果を考慮したリ

トバック型、雁行型、Ｌ型の他、ほぼ任意形

アルな線源分布を再現するには、建物の立つ状の建物モデルが作成可能となった。

地面をリアルに再現することは勿論、建物②屋内空間の再現形状も反映させる必要がある。これまでの木造家屋のような薄い木材を材料とした仕3D-ADRESでは、既存の古田ら文献［16］ が用

切り壁は、家屋内の空間線量率分布に対し、いた四角形の家屋を参考にした簡易モデルを大きな影響を及ぼさないと推定されている。

用いてきたが、実環境を反映し多様な形状のしかし、工場や学校、大型施設等のコンクリー建物モデルを構築可能とすることが望まし

ト壁をもつ建物では、その影響は無視できない。その一方、PHITSでは、複雑な三次元い。また、環境省の委託研究「被ばく線量に体系を容易に構築するため、ポリゴン体系で影響を及ぼす住家内外の要因に関わる研究」

ある四面体メッシュ体系の作成機能が、最近（主任研究者：吉田浩子）に参加し、室内空導入された 文献［17，18］。四面体メッシュ体系間の空間線量率を実際に測定した所、家屋内はシンプルで且つ自由に複雑な三次元体系がにおいても空間線量率の分布があること、更作成可能である。こうして、柔軟に様々な形には、一室の室内においてさえも分布がある状表現が可能になる四面体メッシュ体系を用ことが判明している。従って、屋内のリアルいることで、複雑な家屋モデルの作成が可能な空間を作成可能とするツールを追加し、屋となった。

内空間の再現性を高めることが必要であるこ

図5

多様な建物モデルの形状の再現：（左：既存のボックス型モデル、右：四面体メッシュを用いた多様な建物モデルの例）－9－

RISTニュース

No．67（2021）とが分かった（図 6 参照）。現行3D-ADRES（b）地面除染の寄与が明確化されることで、では、上記ニーズに従い、屋内空間における除染効率等の予測が可能となる。更には、（c）

仕切り壁の作成や仕切り壁の材質の種類や厚家屋撤去の効果も推定可能となり、様々な計さ等を任意に設定可能とした。この改良によ画の妥当性が事前に検討可能となる。

り、屋内空間での線量率分布の解析精度が向3．
2．適合細分化格子による自動ネスト機能上することが期待される。の開発（集落環境を対象とするモデル開発）
3．1．3．シナリオ機能の開発

福島現地での空間線量率の分布を観測する

家屋室内における空間線量率分布を分析すと、集落毎に分布状況が異なる一方、ある一る際、家屋周りの環境面に分布する放射性セ定の特徴的分布を示すこと等が凡そ判明してシウムが、「環境面」毎に、どれ程の割合でいる。従って、対象とする集落をシミュレー

屋内に影響を及ぼすかを定量的に分析するこションするには、計算精度を可能な限り向上とが求められる。もし、このような分析が可させると同時に、計算コストを低減させるこ能になると、効率的な除染作業等の計画が立とを目的とし、解析の対象となる指定集落を案可能となる他、家屋の解体や樹木の伐採等最適な細分化格子で区分するための自動ネスによる効果の事前評価も可能となり、住民やト機能を開発した。一例として、格子作成前国・自治体に対し、意思決定に当たり重要な後の画像を図 8 に示す。図 8 にて分かるよう情報源の提供が可能となる。に、中心位置と中心からの距離の設定、距このような想定されるシナリオを基に、シ離毎の格子サイズの指定を可能としたことミュレーションを実施し、その影響を評価可

で、中心座標の設定の下、半径50m、100m、能とするシナリオ機能を3D-ADRESにて実150m、200m等の領域に円を描画後、範囲毎装した。この機能を使用し、実際の場面で分に指定したサイズで地形格子を生成すること析可能となる例を図 7 に示した。例えば、①が可能となった。これにより、局所的に解像室内の空間線量率に対し、（ 1 ）家屋に付着度の高い格子が作成可能となり、解析対象とした線源が及ぼす寄与、（ 2 ）周囲の樹木線する集落を中心に空間線量率分布の詳細計算源が及ぼす寄与、（ 3 ）地面線源が及ぼす寄と各線源からの寄与率等の分析が可能となった。

与が、各々、明確化する。次に、②除染やそ図 8 に示した適合細分化格子を用いた試計の他の人間活動に対し、（a）樹木伐採の寄与、算の例として、細分化格子を任意の領域に配

図6

屋内空間のモデリングのリアリティーを向上させる：上記のように任意に仕切り壁（材質も選択可）等の導入を可能とした。

－10－

RISTニュース①室内の空間線量率に対する周囲環境の影響評価（（ 1 ）～（ 3 ））②除染やその他の人間活動に対する影響評価（（a）～（c））

図73D-ADRESに実装したシナリオ機能の適用例のイメージ－11－

No．67（2021）

RISTニュース

No．67（2021）

置させた環境モデルを作成し、PHITSによかる。その他、PHITSにより動作確認を行っり空間線量率の計算を行った結果を図 9 に示た結果、線源分布（図 9（c））や空間線量率す。図 9 のモデル作成では、中心からの半分布の計算結果は、（図 9（d））が示す通り

径20mを 3 mメッシュ、半径40mを 5 mメッとなり、複雑な分布の再現が十分に可能とシュ、半径60mを10mメッシュ、半径80mをなっていることが分かる。

20mメッシュ、半径80m以降を40mメッシュで作成し、Surface及びCellを合わせて2.7万

4 ．3D-ADRES開発の成果と今後の展望程度のメッシュ数とし、線源を分布させ計算

4．
1．線量低減係数と除染・人間活動を考慮した空間線量率の将来予測を行った試計算の結果を示す。試計算では、一目で区別可能とする為、土壌の舗装面には

今 般、 新 し く 開 発 さ れ た 機 能 を も つ

非舗装面の10倍程度の線源を設置した（デ3D-ADRESは、二つの環境省委託事業であモのため非現実的な設定を実施）。計算結果

る放射線健康管理・健康不安対策事業等（①より、適合細分化格子による自動ネスト機能「被ばく線量に影響を及ぼす住家内外の要因を用い、計算領域の地面モデルを作成するこに関わる研究」主任研究者：吉田浩子、②「市とで（図 9（a））、計算領域の中心を詳細に街地の空間線量率の経時変化に対する人為的再現しながら、中心から離れるにつれ、メッな影響に関する研究」主任研究者：吉村和シュを大きくすることで効率的な計算が可能也）で活用されつつある。一例として、①のとなった。なお、道路等の比較的、幅が狭く、事業では、木造家屋をターゲットとし、屋内既存の規則メッシュでは再現が不十分であっ空間における空間線量率が、家屋周りのどのた部分も改善され、環境モデルでの現場の環「環境面」から影響を受けるのか？という課境の再現性が向上した（図 9（b））ことが分題について評価することを目指しているが、（a）格子作成前

図8（b）格子作成後最適な細分化格子の自動作成機能の画面（3D-ADRESのGUI画面）－12－

RISTニュース

No．67（2021）

不規則三角メッシュの開発より、家屋周りのするための研究が行われている。これについ

土地利用状況や線源分布が細かく再現可能とても、上記のように、集落単位での細分化モなり、精度の高い評価が可能となった。これ

デルの構築が可能となっており、着目地点毎らの機能拡張により、家屋周りの除染効果をの詳細モデリングが有効に機能すると考えら予測し、住民にとって、重要な情報である家れる。

屋線量低減係数文献［13］ を、現地の情報に併
4．
2．今後の課題せて各家屋単位で分析することが可能になると考えている。一方、②の事業では、被覆面以上、福島の実環境を可能な限り、リアリで異なる放射性セシウムの動態が、市街地にティーを有するモデルで再現し、より正確でおける空間線量率の減少に及ぼす影響を定量効率的な計算を可能とするための研究開発を的に評価することを目指し、人間活動を考慮行ってきたが、未だ実証及び検証が十分に行した空間線量率の減少要因の評価や将来予測われていない。今後は、上記で述べた研究課（a）計算領域（b）環境モデル（d）空間線量率分布の計算例（c）線源分布の設定

図9

適合細分化格子を用いた一つの集落に対する試計算の例－13－RISTニュースNo．67（2021）

題等を通して、実測データに基づいた解析や村和也）において実施したものを含みます。

モデルの実証及び検証を行い、3D-ADRES

研究開発においては、原子力機構・福島研究の有効性を示すことを予定している。
その上、開発部門の宮原要所長、飯島和毅ディビジョ

シミュレーションの大きなメリットである予ン長、眞田幸尚GL、新里忠史GL、操上広志

測評価や実測だけでは困難な定量的分析を行GL、佐久間一幸研究員、阿部智久技術員にい、3D-ADRESの活用範囲を広げることをは、多くのご教示をいただきました。ここに考えている。

感謝致します。システム計算科学センターからは、シミュレーションに当たり計算資源等5 ．まとめと結論で支援いただきました。武宮博センター長を福島の復興はいち早く実現すべき大きな日始めとするスタッフの皆様に感謝いたします

本社会の課題であるが、その実現に向け、多（特に、システム開発初期に参加した奥村雅くの科学的知見が蓄積されてきたことは極め

彦氏とGISデータ・衛星データ処理に関するて貴重である。本報告で記した3D-ADRES

岩田亜矢子氏の支援に感謝致します）。最後は、福島事故後の環境中での空間線量率分布にPHITS活用に際してはRISTの仁井田浩二を分析・考察する際、これまでに得られてき氏に感謝致します。

た様々な知見を基に、シミュレーションモデ
ルを構築する際のエッセンスが凝集されてお参考文献

り、家屋内の空間線量率に対する周囲の環境［ 1 ］ Kato, H. et al., Temporal changes ofからの影響を推定評価することが可能となっthe ambient dose rate in the forestている他、集落単位での解析も可能となってenvironments of Fukushima

おり、様々な活用が想定される。現在、空間Prefecture following the Fukushima

線量率の比較的低い地域での除染等による線reactor accident, Journal of

量率低減は一段落し、今後は帰還困難区域のEnvironmental Radioactivity,復興へと進んでいくと考えられているが、そvols.193-194, 2018, p.20-26の際、家屋や集落単位で、リアルなモデリン［ 2 ］ 眞田幸尚他，避難指示区域の解除にグと的確な除染等の計画が費用対効果も含め向けた特定復興再生拠点の放射線モて事前にシミュレーション可能となれば、復ニタリングと被ばく評価，日本原子興への貢献が期待できると考えている。

力学会和文論文誌，Vol.20,（2021），
No. 2 ，p.62-73［ 3 ］ Andoh, M. et al., Evaluation of

謝辞
3D-ADRESの開発に当たり、検証データdecreasing trend in air dose rateの取得等のため、福島での環境調査の際、住and ecological half-life within an

民及び自治体関係者に多大なご協力を頂き

80 km range from Fukushima Dai-ました。ここに、深く感謝の意を表します。

ichi Nuclear Power Plant, using

尚、この研究の一部は環境省委託事業「平成car-borne survey data measured

31年度 放射線健康管理・健康不安対策事業

by KURAMA systems up to 2018,（放射線の健康影響に係る研究調査事業）（主」

Journal of Nuclear Science and任研究者：吉田浩子）、「令和 2 年度 放射線Technology, vol.57, 2020, Issue 12

健康管理・健康不安対策事業（放射線の健康［ 4 ］ Yoshimura, K. et al., Initial Decrease

影響に係る研究調査事業）」（主任研究者：吉in the Ambient Dose Equivalent－14－

RISTニュース

No．67（2021）

Rate after the Fukushima Accident

Radioactivity, vol.210, 2019, 105878

and Its Difference from Chernobyl,［11］ Sato, T. et al., Features of Particle

Scientific Reports, vol.10, 2020,and Heavy Ion Transport code

p.3859-1-3859-9

System (PHITS) version 3.02, J.
Nucl. Sci. Technol. 55, 684-690 (2018)［ 5 ］ 放射線量等分布マップ拡大サイト
https://ramap.jmc.or.jp/map/［12］ Kim, M. et al., Simulation study of［ 6 ］ 放射性物質モニタリングデータの情

the effects of buildings, trees and報公開サイトhttps://emdb.jaea.go.jp/

paved surfaces on ambient dose

emdb/

equivalent rates outdoors at three［ 7 ］ Mikami, S. et al., The Deposition

suburban sites near Fukushima

Densities of Radiocesium and theDai-ichi, J. Environ. Radioact., 210,Air Dose Rates in Undisturbed(2019), 105803

Fields around the Fukushima［13］ Yoshida-Ohuchi, H. et al., Review

Dai-ichi Nuclear Power Plant;of reduction factors by buildings for

Their Temporal Changes for Five

gamma radiation from radiocaesium

Years after the Accident, Journal

deposited on the ground due toof Environmental Radioactivity,fallout, J. Environ. Radioact., 210,vol.210, 2019, p.105941-1-105941-12(2019), 105803［ 8 ］ Nakama, S. et al., Temporal［14］ Matsuda, N. et al., Measurements ofDecrease in Air Dose Rate in

air dose rates in and around houses

the Sub-Urban Area Affected byin the Fukushima Prefecture inthe Fukushima Dai-ichi Nuclear

Japan after the Fukushima accident,Power Plant Accident during Four

J. Environ. Radioact., 166 (2017), pp.

Years after Decontamination

427-435

Works, Journal of Environmental［15］ Yoshimura, K. et al., Distribution ofRadioactivity, vols.208-209, 2019,

137Cs on components in urban area

p.106013-1-106013-8

four years after the Fukushima Dai-

［ 9 ］ Malins et al., Calculations forichi Nuclear Power Plant accident, J.

ambient dose equivalent rates inEnviron. Radioact., 178-179 (2017),

nine forests in eastern Japan frompp.48-54

134Cs and 137Cs radioactivity［16］ Furuta, T. et al., JAEA-Researchmeasurements, Journal ofEnvironmental Radioactivity,

2014-003［17］ 古田

vol.226, 2021, p.106456-1-106456-12

琢哉，PHITSへの四面体メッシュ体系の導入，第13回PHITS研究［10］ Saito, K. et al., Summary of temporal会（2017）changes in air dose rates and［18］ Furuta, T. et al., Implementation

radionuclide deposition densities inof tetrahedral-mesh geometry in

the 80 km zone over ﬁve years afterMonte Carlo radiation transportthe Fukushima Nuclear Power Plant

code PHITS, Phys Med Biol. 2017

accident, Journal of Environmental

Jun 21;62(12):4798-4810－15－