[論文レビュー] Modeling of X-ray attenuation via photon statistics evolution

Marc Guirao, Sergey Leble|arXiv (Cornell University)|Nov 29, 2014

Medical Imaging Techniques and Applications参考文献 10被引用数 4

ひとこと要約

本稿は、コルモゴロフ前方方程式の1回散乱近似を用いて、ベリリウム膜を通ったX線の減衰をモデル化し、時間分解強度および遅延計算による光子統計の進化を追跡する。有限点検出器における強度について閉形式の式を導出し、散乱角および経路長に明示的な依存性を示しながら、散乱によるパルスの広がりと時間遅延を示す。

ABSTRACT

We consider a formulation of Cauchy problem for Kolmogorov equation which corresponds a localized source of particles to be scattered by medium with given scattering amplitude density. The multiple scattering amplitudes are introduced and the corresponding series solution of the equation is constructed. We investigate the one-fold integral representation for the series terms, its estimations and values of photon number of a finite and point receivers. An application to X-ray beam scattering for orthogonal and inclined to a layer is considered.

研究の動機と目的

• コルモゴロフ前方方程式に基づく確率的アプローチを用いて金属層を通過するX線の伝播をモデル化すること。
• 電磁波の力学を一定速度伝播に簡略化しながら、ベリリウム層内の弾性散乱および吸収を考慮すること。
• 1回の散乱イベント後に発生する有限点検出器に到達する散乱X線パルスの強度について解析的表現を導出すること。
• 散乱幾何学および経路長の変動がパルスの時間遅延および時間幅広がりに与える影響を定量化すること。
• 今後の研究において、複数回散乱や不均質媒体へのモデル拡張の基盤を提供すること。

提案手法

• 全断面積および微分断面積を組み込んだ1粒子分布関数に対するコルモゴロフ前方方程式を定式化する。
• 点パルス源を初期条件とし、一定速度で自由空間およびベリリウム層内を通過する光子の軌道をモデル化する。
• 球座標を用い、コサイン法則を用いて散乱角関係を記述する：$\cos\gamma = \cos\theta\cos\theta' + \sin\theta\sin\theta'\cos(\phi - \phi')$。
• 経路長および飛行時間依存性を組み込んだ1重積分表現として散乱強度を導出する。
• 時間関数としての強度を、$I/I_0 = \frac{1}{t_2^2} \exp[-\sigma_t(t_3 + t_4)] \sigma(\cos\gamma, z)$ を用いて評価する。ここで、$t_3$ および $t_4$ は層通過前後における時間成分である。
• 特定のパラメータ（$\vartheta_0 = 0.1$, $z_0 = 0.01$, $\Delta = 0.001$）を用いた数値プロット（SWP）により、強度対時間の関係を可視化する。

実験結果

リサーチクエスチョン

• RQ1 コルモゴロフ方程式の1回散乱近似は、ベリリウム膜を通過した後の検出器におけるX線強度分布をどのようにモデル化できるか？
• RQ2 散乱X線パルスの時間的進化はどのように変化し、散乱角および経路幾何学にどのように依存するか？
• RQ3 散乱および吸収断面積は、検出パルスの強度および遅延にどのように影響するか？
• RQ4 前方散乱幾何学における点受信機について、強度の時間関数としての解析的形はどのようなものか？
• RQ5 このモデルは、複数回散乱や不均質媒体を含めるためにどのように拡張可能か？

主な発見

• 検出器における強度は $I/I_0 = \frac{1}{t_2^2} \exp[-\sigma_t(t_3 + t_4)] \sigma(\cos\gamma, z)$ で与えられ、$t_2$ は散乱点から検出器までの時間を表す。
• パルスの到着時刻は遅れ、広がりを示す。選択したパrametersでは、開始時刻 $t_0 = 3 \times 10^{-2}$、終了時刻 $t_e = 7.5 \times 10^{-2}$ となる。
• 飛行時間は散乱点の高さ $z$ に依存し、$z = \frac{1}{2} \cos\vartheta_0 (t - z_0) \frac{t + z_0}{t - z_0 \cos\vartheta_0}$ として表される。
• 経路長および吸収の増加に伴い相対的強度が減少し、指数減衰項 $\exp[-\sigma_t(t_3 + t_4)]$ で捉えられる。
• モデルは、有限で遅延され、広がったパルスが検出器に到達することを予測しており、ピーク強度は最小経路長に対応する時刻に現れる。
• 結果は、散乱が測定可能な時間遅延および強度低下を引き起こすことを示しており、導出された解析的表現を用いて定量的に評価可能である。