医療画像における放射線線量計算の精度向上:ICRPメッシュ型ファントムによる電子SAFデータセット

解説対象論文: Electron specific absorbed fractions for the ICRP family of mesh-type reference computational phantoms
Physics in Medicine and Biology|被引用数: 0(2026年10月23日時点)
本記事では、国際放射線防護委員会(ICRP)が開発した最新のメッシュ型リファレンス計算ファントム(MRCPs)を用いた電子比吸収割合(SAF)の包括的なデータセット作成に関する研究について解説します。従来のボクセル型ファントム(VRCPs)では捉えきれなかったマイクロメートルスケールの放射線感受性組織や臓器内血液の含有をMRCPsが明示的にモデル化することで、内部被ばく線量評価の精度が大幅に向上します。PHITSモンテカルロコードとコンポーネントベースの方法を用いた計算結果は、Geant4による独立した検証と高い一致を示し、特に幼い対象者では年齢依存的なSAF値の増加が確認されました。このデータセットは、放射線防護基準の更新と、医療・職業被ばくにおける内部線量評価のゴールデンスタンダード確立に貢献する重要な基盤となります。
はじめに:放射線被ばく線量評価の進化
どうも、Beyond the Pixelです。医療画像インフォマティクスが日々進化する中、放射線被ばく線量の正確な評価は、診断や治療における安全性を確保するための重要な基盤となります。特に体内に入った放射性物質(放射性核種)からの内部被ばく線量の推定は、放射線防護の分野で中心的な役割を果たします。体内に入った放射性核種は、医療目的での投与、吸入、摂取、創傷汚染といった経路を通じて体内に分布し、特定の臓器や組織に保持されます。この放射性核種が時間とともに減少していく様子(生体内動態モデル)は、国際放射線防護委員会(ICRP)によって詳しく記述されています。内部被ばく線量の平均吸収線量(D)は、放射性核種が均一に分布した線源組織から標的組織に到達するエネルギーを考慮して計算されます。この計算において、S値(S係数)と比吸収割合(SAF)が鍵となります。S値は、放射性核種の放射線放出特性と、線源と標的の解剖学的構造の両方を考慮した値であり、SAFは、線源組織から放出されたエネルギーが標的組織に吸収される割合を示します。

この研究で用いられた線源領域のリストと、

標的領域のリストです。歴史的に、電子線量計算では、すべての電子およびベータ粒子のエネルギーが線源臓器内で局所的に吸収されるという簡略化された仮定が用いられてきました。これは、骨内線量測定の場合を除き、ICRP Publication 30(1979年)で正式に採用された方法でした。しかし、この仮定は電子エネルギーが約300~500 keVを超えると破綻し、高エネルギー電子が臓器の境界を越えて隣接する組織に無視できない量のエネルギーを沈着させることが明らかになりました。解剖学的にリアルな臓器形状を提供するボクセル型計算ファントム(VRCPs)の登場により、放射線輸送コードで荷電粒子を直接シミュレーションすることが可能になり、この局所沈着の仮定は見直されることになりました。ICRPは、VRCPsを用いて光子、電子、アルファ粒子、核分裂スペクトル中性子の包括的なSAFデータセットをPublication 133および155で作成しました。しかし、VRCPsにもいくつかの解剖学的・幾何学的な限界が指摘されていました。例えば、臓器体積は実質組織のみを定義し、臓器内の血液を含んでいなかったため、SAF値に事後補正が必要でした。また、ボクセル解像度では、呼吸器や消化管の粘膜、皮膚の基底層、膀胱壁といったマイクロメートルスケールの放射線感受性標的組織層を直接モデル化することができませんでした。これらの課題に対応するため、ICRPはメッシュ型リファレンス計算ファントム(MRCPs)を開発し、Publication 145および156で公開しました。

この図は、ICRPのメッシュ型リファレンス計算ファントム(MRCPs)を示しており、新生児から成人までの男女の参照対象者が含まれています。本研究は、ICRPのMRCPsファミリーすべてについて、拡張された線源・標的組織の組み合わせを持つ完全な電子SAFデータセットを初めて提供するものであり、内部線量測定の精度をさらに高めることを目的としています。
メッシュ型リファレンス計算ファントム(MRCPs)の力
ICRPのMRCPsは、本研究における電子SAF計算の解剖学的基盤として用いられました(ただし、骨内幾何学構造を除く)。MRCPsは、ICRPのVRCPsを四面体メッシュ形式に変換し、マイクロメートルスケールの放射線感受性標的組織の解剖学的にリアルなモデルを組み込み、VRCPsの既知の解剖学的欠点をいくつか修正することによって開発されました。これらの修正には、臓器内の血液の明示的な追加、骨格系、結腸、小腸、胸郭外気道、眼、歯、甲状腺のモデルの改訂が含まれます。電子線量測定の観点から見ると、四面体メッシュ形式は光子輸送に関連する利点に加えて、さらに優位性を提供します。メッシュ形状による滑らかな臓器境界は、エネルギー沈着の大部分が標的境界近くに集中する短距離粒子の線量計算におけるアーティファクトを低減します。MRCPs内に薄い放射線感受性標的層を明示的に定義することで、サロゲート組織や外部様式モデルを使用することなく、自己整合的な電子SAF計算が可能になります。さらに、MRCPsへの臓器内血液の組み込みは臓器の体積を増加させ、臓器境界間の適合関係を変化させます。
電子SAFの計算方法:PHITSとコンポーネントベースの手法
本研究では、ICRPの成人および小児のMRCPsを用いて、定義されたすべての線源および標的領域((表1)と(表2)参照)について電子SAFを計算しました。電子SAFデータセットを生成するために、PHITS(Particle and Heavy Ion Transport code System)バージョン3.35という汎用モンテカルロ粒子輸送シミュレーションコードが用いられました。電子SAFの総量は、衝突による寄与と制動X線による寄与を分離するコンポーネントベースの方法を用いて計算されました。衝突成分は、電子が物質と直接相互作用することでエネルギーを失うことによるもので、放射成分は、電子が物質中で減速する際に放出される制動X線(Bremsstrahlung X-rays)によるものです。後者の放射成分は、光子SAFデータセットから導出されました。

この図は、5歳女児のMRCPにおける肝臓、甲状腺、海綿骨体積(T-bone-V)、第2胸郭外領域結合部(ET2-bnd)の線源領域における様々な線源電子エネルギーで生成された制動光子数を示しています。

この図では、5歳女児のMRCPにおける様々な線源と標的の組み合わせにおける、放射成分、衝突成分、および総電子比吸収割合(SAF)が示されています。線源から標的までの距離が電子の飛程内にある場合、電子は衝突相互作用を通じて標的領域に多くのエネルギーを沈着させ、放射成分を凌駕することがわかります。骨格組織外の線源から赤色骨髄(R-marrow)および骨内膜表面(Endost-BS)への衝突電子SAFは、特定の式を用いて計算されました。また、呼吸器の複雑な構造内でのSAFについては、Shinらの先行研究(2025年)で発表されたデータセットが採用され、本研究のSAFデータセットに追加されました。骨格系内の線源からの電子SAFについては、VRCPでは捉えきれなかった微細構造モデルから導出されたICRP Publication 133および155のデータセットが、現在のところ最も完全なセットとして用いられました。
品質保証と低エネルギー補外
本研究で報告された電子SAF値が、使用された放射線輸送コードのアーティファクトではないことを確認するため、成人男性MRCPにおける線源・標的組織のサブセットについて、Geant4バージョン11.01.p02を用いて独立にシミュレーションを行いました。この品質保証計算では、PHITSの単一シミュレーション法、PHITSのコンポーネントベース法、およびGeant4の単一シミュレーション法の3つのモンテカルロ手法間で比較が行われました。

この図は、成人男性のメッシュ型リファレンス計算ファントム(MRCP)における、PHITS単一シミュレーション法とPHITSコンポーネントベース法、Geant4単一シミュレーション法とPHITS単一シミュレーション法、Geant4単一シミュレーション法とPHITSコンポーネントベース法における比吸収割合(SAF)値の比率の分布を示すボックスプロットです。各ボックスは16個のデータポイントを含み、いずれかの値がゼロのペアは除外されています。これら3つの方法は、全エネルギー範囲およびすべての線源・標的組織の組み合わせにおいて、おおむね一致を示しました。PHITSの2つの方法間の平均比率は1.10、中央値は1.00でした。また、Geant4とPHITSの単一シミュレーション間の平均比率は0.95、中央値は0.98でした。PHITSのコンポーネントベース法はGeant4の単一シミュレーション結果を平均比率1.06、中央値1.00で再現しました。ICRP Publication 133および155で確立された低エネルギーSAF挙動の解析的表現に従い、電子SAFの総値はゼロエネルギーへの制限値へと補外されました。非重複する線源・標的組織のペアでは、粒子エネルギーがゼロに近づくにつれてSAFはゼロに近づきます。自己吸収の場合、つまり線源と標的が同じ領域である場合、放出されたすべてのエネルギーは局所的に吸収され、SAFは標的質量(kg)の逆数に収束します。血液が線源領域である場合、総血液量の一部が各臓器内に存在するため、全身の標的組織と部分的に空間的に重なり合います。この場合のゼロエネルギー制限SAFは、標的組織内に含まれる総血液量の割合で決定されます。同様の部分的重なり条件は、消化管内の一部の線源・標的組織の組み合わせでも発生し、線源領域体積の定義された部分が標的体積と一致する場合、制限SAFは標的内の線源の体積比によって決定されます。
重要な結果:年齢依存性とVRCPとの比較
本研究で生成された電子SAFデータセットは、すべてのICRP MRCPsファミリーについて包括的なもので、5332の線源・標的組織の組み合わせと28のエネルギー点(0~10 MeV)をカバーしています。年齢依存性の傾向がMRCPsファミリー全体で観察され、年齢が減少するにつれてSAF値は増加する傾向にありました。これは、幼いファントムでは臓器質量が小さく、臓器間距離が短いため、より高い吸収割合(AF)が生じることを反映しています。

この図は、線源が肝臓である電子被ばくにおける比吸収割合(SAF)を、標的組織が肝臓、肺胞間質(AI)、直腸S状結腸の幹細胞層(RS-stem)、赤色骨髄(R-marrow)である場合について、新生児から成人までの各メッシュ型リファレンス計算ファントム(MRCPs)で計算した結果を示しています。

この図は、線源が甲状腺である電子被ばくにおける比吸収割合(SAF)を、標的組織が甲状腺、食道の基底層(Oesophag-bas)、小腸の幹細胞層(SI-stem)、第2胸郭外領域の基底層(ET2-bas)である場合について、新生児から成人までの各MRCPsで計算した結果を示しています。

この表は、自己吸収およびクロスファイアケースにおける小児対成人の電子比吸収割合(SAF)比の中央値、5パーセンタイル、95パーセンタイルを示しています。年齢とともにSAF比の中央値が1に近づく傾向が見られます。ICRP Publication 133および155で報告された値との比較では、自己吸収幾何学構造については広く一致が見られましたが、クロスファイア幾何学構造では相違が確認されました。この相違は、ファントムの解剖学的構造、臓器の体積と臓器間距離、制動X線光子生成結果、低エネルギー外挿法、および放射成分の計算に使用された基礎となる光子SAFデータセットの違いなど、複数の要因の組み合わせに起因すると考えられます。例えば、臓器内の血液を明示的に含めることで臓器の体積が増加し、線源と標的の分離が効果的に減少するため、特定のクロスファイアケースで顕著な比率の違いが生じました。

この図は、ICRP Publication 133および155の値に対するメッシュ型計算ファントム(MRCPs)の比吸収割合(SAF)の比率を示しており、血液量が多い臓器(肝臓と脾臓)および血液量が少ない臓器(筋肉と脳)を含む線源・標的ペアで表示されています。また、VRCPsの「階段状アーティファクト」によって不完全な骨の被覆が生じていた骨標的組織に対する筋肉線源の計算では、本研究の方が低エネルギーで低いSAF値を示しました。これは、筋肉内で発生した電子が皮質骨を通過せずに直接海綿骨や髄腔に到達できる露出領域がないためと考えられます。

この図は、ICRP Publication 133および155の値に対するMRCPsの比吸収割合(SAF)の比率を示しており、筋肉を線源とし、骨を標的とする場合について、VRCPsにおける階段状アーティファクトによる海綿骨および髄腔領域の不完全な被覆の影響を調査するために表示されています。さらに、MRCPsで明示的にモデル化された消化管内のマイクロメートルスケールの標的組織層の影響も観察されました。

この図は、ICRP Publication 133および155の値に対するMRCPsの比吸収割合(SAF)の比率を示しており、MRCPsに明示的に定義された消化管のマイクロメートルスケール標的層を含む線源・標的ペアで表示されています。MRCPs開発中にその幾何学が大幅に修正された臓器においても、臓器依存的な違いが確認されました。

この図は、ICRP Publication 133および155の値に対するMRCPsの比吸収割合(SAF)の比率を示しており、VRCPsからMRCPsへの開発中にその幾何学が修正された臓器および組織を含む線源・標的ペアで表示されています。

この図は、5歳女児における自己吸収幾何学とクロスファイア幾何学での電子エネルギーに対する比吸収割合(SAF)値の比率の分布を示すボックスプロットです。ボックスは、自己吸収幾何学で92組、クロスファイア幾何学で3142組の一致する線源・標的ペアで構成されています。特定の放射性核種(Y-90、I-131、Lu-177)について、5歳女児MRCPとVRCPを用いた肝臓線源からのS値計算も行われました。

この表は、5歳女児MRCPと5歳女児VRCPについて、選択された放射性核種と線源・標的組織の組み合わせで計算された放射性核種S値を示しています。自己吸収およびクロスファイアの両ケースで、MRCPとVRCPの値の間には比較的良好な一致が見られましたが、Y-90の肝臓線源から小腸の幹細胞層への場合が最も大きな差を示しました。
研究の限界と今後の展望
本研究の限界としては、包括的な放射性核種S値の比較評価がなされていない点が挙げられます。本研究では電子SAFの計算に焦点を当てており、MRCPsから計算された光子SAFと電子SAFを用いて、今後、包括的なS値が計算される予定です。これらの新しいS値データは、内部線量測定の精度をさらに高め、放射線防護の実践において重要な貢献を果たすことが期待されます。
用語集
- 比吸収割合 (SAF): 特定の線源組織から放出された放射線エネルギーが、特定の標的組織にどれだけ吸収されるかを示す割合です。
- S値: 放射性核種が特定の線源組織に存在する場合に、特定の標的組織が受ける吸収線量を示す係数です。放射線放出特性と解剖学的構造を考慮します。
- 国際放射線防護委員会 (ICRP): 放射線防護に関する勧告を行う国際的な非政府組織です。
- メッシュ型リファレンス計算ファントム (MRCPs): ICRPが開発した、メッシュ形式で解剖学的構造を詳細に表現した人体モデルで、従来のボクセル型よりも高精度な線量計算を可能にします。
- ボクセル型リファレンス計算ファントム (VRCPs): ICRPが以前使用していた、CT画像データに基づきボクセル(立方体の要素)で人体構造を表現した計算ファントムです。
- PHITS: 粒子と重イオンの輸送をシミュレーションするためのモンテカルロコードで、放射線輸送計算に用いられます。
- Geant4: 物質中の粒子の通過をシミュレーションするためのモンテカルロツールキットで、放射線物理学の分野で広く利用されます。
- モンテカルロ法: 乱数を用いて現象をシミュレーションし、その結果から物理量を推定する計算手法です。
- 制動X線: 電子などの荷電粒子が物質中で減速する際に放出されるX線のことです。
- 自己吸収幾何学: 放射線源となる組織と放射線エネルギーを吸収する標的組織が同一である場合の幾何学的配置を指します。
- クロスファイア幾何学: 放射線源となる組織と放射線エネルギーを吸収する標的組織が異なる場合の幾何学的配置を指します。
コメントを残す