靜態CT行包安檢系統關鍵技術研究

唐利華，張國光，張 健，閆 平，張 雨

(中國原子能科學研究院 核技術應用研究所，北京 102413)

傳統的螺旋CT是通過旋轉滑環從各角度掃描物體，以獲得物體的透射信息，其設計中的一個關鍵限制是繞物體旋轉時可施加到X射線源上的最大加速度。高壓絕緣子上的機械壓力和流過陽極的冷卻油均受此影響。受環上過高的離心加速度及高壓絕緣材料的應力限制，螺旋CT正逐漸接近其轉速的極限[1]。為解決該問題，可考慮使用另一種系統結構。Zhang等[2]提出了一種不同于螺旋CT的掃描方法，這種掃描方法提出了一種在真空室內集成多個線性陣列陰極的新型X射線源。若以適當的掃描順序對該多焦點X射線管中的各X射線源進行脈沖控制，則無需任何物理運動即可獲得投影數據，這種掃描稱為靜態掃描。

靜態CT通過射線源高速切換，可實現相當于每秒幾十轉等效轉速，是螺旋CT轉速的幾倍甚至幾十倍，這樣不但減少了噪聲和振動，也避免了過高的向心加速度，使整體系統更為安全、可靠。另外，靜態CT的能耗只有傳統螺旋CT的一半左右。靜態CT在高速安檢的情況下，結合人工智能大數據，可智能識別安檢違禁品，極大地提高了安檢效率。本文介紹安檢用靜態CT技術的基本原理，對基于超順排碳納米管技術的冷陰極X射線源、基于光子計數技術的多能譜X射線成像技術和安檢用靜態CT系統圖像重建算法等進行研究。

1 基于超順排碳納米管技術的冷陰極X射線源

2002年，文獻[3]發現了碳納米管的超順排陣列，與普通陣列相比，超順排碳納米管陣列具有更直的碳納米管、更整齊的排列、較好的場發射性能和穩定的結構(圖1)。使用超順排碳納米管陣列制作的場發射器件易于操作，且可輕松集成到真空電子設備中[4]。本文研制一款以超順排碳納米管陣列作為冷陰極的X射線管，并進行測試。當陽極電壓低于130 kV時，研制的冷陰極X射線管可穩定工作。當陰極和柵極之間電勢差為940 V時，此冷陰極X射線管達到發射電流1.0 mA。陰極場發射的伏安特性曲線如圖2所示。通過小孔成像原理，測出了該X射線管的焦點尺寸，大小為0.6 mm×1.6 mm，符合設計目標。

圖1 碳納米管粉末(a)、普通陣列(b)和超順排陣列(c)的SEM圖像[5]

圖2 陰極場發射的伏安特性曲線

研制的冷陰極X射線管是一款場致發射X射線管，超順排碳納米管陣列陰極在柵極施加的電場引發下，發射出電子。電子被陽極電場加速后，撞擊到鎢靶上產生張角為80°×20°的扇形束X射線。冷陰極X射線管和高壓系統全部浸入絕緣油中以實現高壓絕緣和冷卻。

該X射線源使用兩組高壓電源系統，一組電源是+160 kV單極性電源，可在0～160 kV之間連續調節。另一組電源是-3 kV單極性電源，可使用直流模式也可使用脈沖模式。該冷陰極X射線源可瞬間開關，開關時間在μs級別，脈沖模式可在5～100 ms之間連續可調?？刂葡到y通過管電流的反饋，可自動調節柵極電壓，以保持管電流恒定。

2 基于光子計數技術的多能譜X射線成像技術

在電荷積分式探測器的成像系統中，高能光子對投影數據的貢獻更大，但系統所得到的圖像對比度很低，導致無法區分密度差異很小的物質。在基于光子計數技術的多能譜X射線成像技術中，不同的光子能量可具有不同的權重，從而有效地提高了圖像對比度[6]。探測器可設置探測的能量閾值，以有效地識別和消除噪聲，相應的成像方法稱為能量加權成像[7]，它在改善安全檢查成像的圖像質量中起著重要作用，因此，可期望其在多種應用中實現更好的圖像質量并改善材料的辨別能力[8-10]。本文研發一款基于硅酸釔(YSO：Ce)晶體和硅光電倍增管(SiPM)的基于光子計數技術的多能譜X射線探測系統，并進行測試。

2.1 探測器陣列和數據獲取系統

探測器是布置成線性陣列的光子計數型三能譜X射線探測器，其由YSO：Ce、SiPM和數據采集(DAQ)系統組成。YSO：Ce晶體使用硅油與SiPM進行光學耦合。YSO：Ce晶體將X射線光子轉換為可見光，從而在SiPM中生成電子脈沖，來自SiPM的電脈沖連接到前置放大器的輸入。DAQ系統根據預設的3個能量閾值對放大的脈沖數進行分類和計數，并將3個能量段的X射線強度傳輸到主機，所有能量閾值均可編程。

1個標準探測器模塊板上裝有4個線性排列的YSO：Ce探測器陣列，該模塊還包含電源電路、數據采集控制電路、讀取和傳輸電路及用于芯片配置的控制器。圖3為標準探測器模塊板(配備4個YSO：Ce探測器陣列，4×16像素，間距為1.2 mm×1.6 mm)和主板，最多可將10塊模板放在一起以構建掃描臂，掃描臂的有效長度是102.6 mm的倍數。

圖3 標準的探測器模塊板(a)和主板(b)

利用X射線照射尺寸為1.2 mm×1.6 mm×5.0 mm的線性YSO：Ce探測器，測試其計數率，測試結果如圖4所示。當X射線管電流在0.2～0.5 mA之間時，輸出的計數率線性增大，計數率在2.5×106s-1時達到飽和。

溫度漂移(溫漂)通常是指環境溫度變化時引起晶體管參數的變化。這可能導致靜態工作點和電路動態參數不穩定，甚至會阻止電路正常工作。該探測器系統具有溫漂修正系統，溫漂修正系統會根據環境溫度的變化自動調整偏差值，以確保探測器工作的穩定性。

圖4 計數率隨X射線管電流的變化和線性回歸

2.2 三能譜X射線成像技術

傳統電荷積分式探測器是通過能量積分獲取透射數據的，而三能譜X射線成像技術是通過1次掃描，根據預設的3個能量閾值對放大的脈沖數進行分類和計數，并將3個能量段的X射線強度傳輸至主機。主機依據3個能量段獲取的數據重構出3個能量段的圖像。圖5為光子計數型探測器與傳統高低能探測器的比較。由于物質對不同能量X射線衰減系數不同，在不同能量下得到的圖像會有不同的光密度分布。根據能量加權成像方法，經計算可分別給出被攝物品近似的等效原子序數。

采用管電壓為160 kV、管電流為1 mA的X射線，對該探測系統成像效果進行測試，得出3幅不同能量段的行李箱圖像(圖6)。圖6a、b、c分別為行李箱的低、中、高能段X射線透射圖像。在低能段的透射圖像中，水和塑料瓶等低密度物質可清晰地顯示；在高能段的圖像中，低密度物質成像較模糊。根據能量加權成像方法，經計算可分別給出被攝物品近似的等效原子序數。

a——傳統高低能探測器；b——光子計數型探測器

圖6 行李箱的三能譜成像

3 靜態CT圖像重建技術

靜態CT系統的掃描方式實質上為螺旋掃描，重建的第1步是通過修正得到與旋轉軸垂直的切片投影數據，可采用4π插值方法獲得該數據，然后采取等角扇形束重建方法進行重建圖像，則：

f(x,y)=

其中：h″(γ)為濾波函數；γ為射線與中心射線形成的角度，即探測器角度；β為中心射線與y軸的夾角，即投影角度(或視角)；L為射線源到重建點(γ，φ)的距離；D為射線源到中心的距離；q(γ,β)為扇形束投影[11]。

在CT系統中，圖像重建算法通常需完備的投影數據，在很多實際應用中，由于受數據采集時間或成像系統掃描的幾何位置約束，只能在有限角度范圍或在較少的投影角度得到數據，這些都屬于數據不完全問題[12]，本系統即屬于這種情況。利用基于壓縮傳感理論的稀疏角度重建算法，可在原有掃描方式上直接減少掃描角度數量，從而有效降低掃描過程總的照射劑量和掃描時間，也有利于在快速成像應用中獲得更好的成像效果[12-15]。在稀疏角度重建問題中，選擇使用基于TV最速下降和凸集投影約束相結合的ASD-POCS算法[16-18]。

使用Matlab軟件模擬擁有不同數量球管的靜態CT系統成像效果。設定每個射線源的扇形束射線張角為60°、射線源環的直徑為1 m、有效成像視野直徑為0.5 m，則使用20個球管的靜態CT系統，重建圖像的形狀基本可分辨出，球管數量越多，重建圖像空間分辨率越高(圖7)。

圖7 不同數量球管的shepp-logan頭模型重建圖像

4 結論

本文研制了一款基于超順排碳納米管技術的冷陰極X射線管，在130 kV以下，該X射線管可穩定工作。設計了一款直流和脈沖雙模式冷陰極X射線源，開關時間在μs級別，脈沖模式可在5～100 ms之間連續可調。研發了一款由YSO：Ce、SiPM和數據采集系統組成的三能譜光子計數型X射線探測系統，計數率在2.5×106s-1時達到飽和，得到3個能譜的成像結果。討論了靜態CT圖像重建的算法，提出了一種可行的解決途徑。使用Matlab軟件模擬了擁有不同數量球管的靜態CT系統成像效果，靜態CT系統關鍵技術的研究為最終研制靜態CT系統提供了技術基礎。