宇宙射線μ子探測高Z材料的仿真研究

龐洪超，劉森林，王紅艷，徐勇軍，楊宏偉，王列民

（中國原子能科學研究院 輻射安全研究部，北京 102413）

2003年美國洛斯阿拉莫斯實驗室首先報道可利用宇宙射線μ子對重核物質輻射成像［1］，引起國內外科學家的廣泛關注［2－5］。天然宇宙射線μ子能量高，穿透能力強，不會對環境和人員造成額外輻射，在探測重核物質方面具有很大優勢。本課題組也對μ子探測重核物質進行了前期研究并取得一定成果［6－7］。本文在前期研究的基礎上，采用Geant4建立宇宙射線μ子探測高Z材料的仿真系統，研究漂移管探測器位置分辨率及系統成像時間對成像效果的影響，同時為漂移管的設計加工提供理論依據。

1 仿真模型的建立

圖1為宇宙射線μ子成像檢測系統仿真模型結構圖，上下各兩組探測器，中間區域為有效探測區域。當μ子穿過上兩組探測器時，能準確確定μ子的入射位置坐標，當μ子穿過待測物并從下面兩組探測器射出時，通過同樣方法可確定μ子的出射位置坐標。通過探測器記錄的數據可計算出μ子的位置和角度偏移信息，從而進一步確定探測區域是否含有高Z物質。

圖1 Geant4構建μ子探測系統仿真模型Fig.1 Model of muons detector system using Geant4simulation

探測器選擇位置靈敏探測器漂移管（DT），DT由外直徑30mm、壁厚0.6mm、長1 000mm的鋁管（陰極），端頭和直徑為50μm的鍍金鎢絲（陽極）組成，管內充氣壓為300kPa的氬和二氧化碳（體積比為93∶7）的混合氣體。

為增加模擬真實性，編寫了宇宙射線μ子產生器，使μ子的角分布和能譜分布與真實海平面μ子一致［7－8］。

探測空間為1m×1m×1m的立方體，內裝兩個物體，其中1個為直徑10cm的鐵球，另1個為直徑8cm的鈾球，其中鈾球的外面屏蔽一厚為1cm的鐵殼。上兩層探測器間距離與下兩層探測器間距離均為20cm。采用最大似然法［4－5］進行模擬計算。

2 漂移管位置分辨率對成像結果的影響

μ子產生器產生1min的入射粒子，最大似然算法迭代100次，探測器位置分辨率分別取0、200、300、400、500、1 000μm。成像顯示以像素的散射密度λ進行劃分，像素大小為10cm×10cm×10cm。當5＞λ≥1時，表示為低Z物質，用淺灰色顯示；15＞λ≥5，表示中Z物質，用深灰色顯示；λ≥15，表示為高Z物質，用黑色顯示（圖2）。由圖2可看出，隨探測器位置分辨率的增大，系統的成像效果逐漸變差。位置分辨率小于200μm時，在鈾球的位置出現黑色，鐵球的位置出現深灰色，并出現幾個淺灰色干擾點，成像效果理想。當探測器的位置分辨率取400μm時，成像區域出現較大的深灰色干擾點，成像效果變差，當探測器的位置分辨率為1mm時，雖然鈾球的位置還能出現黑色但深灰色干擾點較多，已無法區分物體。可看出，探測器的位置分辨率優于200μm時可得到理想的成像結果。

3 成像時間對成像結果的影響

選取探測器的位置分辨率為200μm，算法迭代次數取100，μ子的入射時間分別為30s、1min、2min，得到的成像結果示于圖3。由圖3可看出，在探測器的位置分辨率及算法迭代次數一定時，成像效果隨入射μ子數增加逐漸變好。當系統探測時間為30s時，鈾球位置出現黑色，而鐵球位置未出現，其它像素出現干擾點，成像效果差；當探測時間為1min時，鈾球和鐵球均出現，且出現幾個小散射密度點；當時間為2min時，效果非常理想，出現了鈾、鐵散射點，且干擾散射點的散射密度均小于1。可見，系統成像時間在1～2min內能得到理想的成像結果。

圖2 成像效果與探測器位置分辨率的關系Fig.2 Relationship between imaging and detector resolution

圖3 成像效果與系統探測時間的關系Fig.3 Relationship between imaging results and detecting time

4 結論

1）利用宇宙射線μ子實現高Z材料的快速探測是可行的。在漂移管的位置分辨率為200μm左右、成像時間1～2min內可得到理想的成像結果。

2）探測μ子的漂移管需具有200μm以內的位置分辨率，高于此值，辨別物體就存在困難，即本文的研究也為探測器的設計加工提供了理論依據。

3）本文的研究是基于蒙特卡羅方法的模擬研究，且探測區域情況較為理想，要將此方法進行實際應用，需做大量的實驗工作。

［1］BOROZDIN K，HOGAN G，MORRIS C，et al.Radiographic imaging with cosmic－ray muons［J］.Nature，2003，422（6929）：277.

［2］SCHULTZ L J，BLANPIED G S，BOROZDIN K N，et al.Statistical reconstructions for cosmic ray muon tomography［J］.IEEE Trans Imag Proce，2007，16（8）：1 985－1 993.

［3］WANG Guobao，QI Jinyi.Statistical image reconstruction for muon tomography using Gaussian scale mixture model［C］∥IEEE International Conference on Image Processing.［S.l.］：［s.n.］，2008：2 948－2 951.

［4］SCHULTZ L J.Cosmic ray muon radiography［D］.Portland：Portland State University，2003.

［5］劉圓圓，趙自然，陳志強，等.用于宇宙射線μ子成像的MLS－EM重建算法加速研究［J］.CT理論與應用研究，2007，16（3）：1－5.LIU Yuanyuan，ZHAO Ziran，CHEN Zhiqiang，et al.Research on accelerating MLS－EM reconstruction algorithm in cosmic ray muon radiography［J］.CT Theory and Applications，2007，16（3）：1－5（in Chinese）.

［6］羅煒，王紅艷，陳凌，等.用于核材料檢測的兩種宇宙射線μ子位置靈敏探測器［J］.核電子學與探測技術，2009，29（2）：353－356.LUO Wei，WANG Hongyan，CHEN Ling，et al.Two types of position－sensitive detectors about universe radial of muon apply detecting nucleus material［J］.Nuc1ear Electronics ＆ Detection Technology，2009，29（2）：353－356（in Chinese）.

［7］譚新建，馮朝陽，謝一岡，等.宇宙線μ用于反應堆燃料元件成像的嘗試性模擬研究［J］.原子能科學技術，2010，44（3）：344－347.TAN Xinjian，FENG Zhaoyang，XIE Yigang，et al.Simulation study on cosmic rayμtomography for reactor fuel elements［J］.At Energy Sci Technol，2010，44（3）：344－347（in Chinese）.

［8］GASSER T K，STANEV T.Review of particle physics［J］.Phys Lett B，2004，592：230.