,,,,

(中国兵器科学研究院宁波分院，宁波 315103)

工业CT检测中主要工艺参数定量取值方法

张维国,倪培君,王晓艳,齐子诚,唐盛明

(中国兵器科学研究院宁波分院，宁波 315103)

工业CT作为一种先进的无损检测技术，应用范围已越来越广泛。目前在用的工业CT系统已有数百套，这些工业CT因为硬件配备的不同，工艺参数也有很大差别。现行的工业CT标准中对主要工艺参数的取值给出了指导性原则，但没有给出定量取值方法而导致操作不便。从检测需求出发，依据采样理论分别给出了主要参数的取值公式，方便了CT检测工艺的编写，并用实例给出了操作过程，并对方法的正确性进行了验证。

工业CT; 工艺参数;定量

工业CT作为一种先进的无损检测技术，已被广泛应用于航天、航空、兵器、汽车、考古、地质等多个行业领域[1-2]。据不完全统计，国内已安装的工业CT系统已超过200多套，而且在以每年5～10套的速度增加。目前安装的系统中，从探测器和扫描装置的不同可以分为：线阵三代、线阵二代以及三维锥束面阵。从射线焦点及能量的不同可以分为：纳米CT、微焦CT、常规X射线源CT以及加速器工业CT系统。上述系统的生产厂家、型号、使用的探测器种类、扫描装置及射线源种类各不相同，其系统的工艺参数差别也就很大。现行的通用CT标准中，对于工业CT参数的取值都给出了一些指导性原则，但很少给出定量取值方法。笔者从检测需求出发，依据采样理论分别给出CT检测中常用工艺参数的定量取值方法，为工艺卡的制定提供帮助[3-8]。

1 工业CT扫描工艺参数取值分析

1.1 基本原则

工业CT检测中，CT图像质量与检测时间及系统开销是相关的。因此检测过程必须要以检测需求为基础，制定其他相关的工艺参数。

在检测前，需要了解以下的信息。

检测设备信息：射线源焦点尺寸a，探测器通道尺寸d，以及系统的最佳空间分辨率L。

试样信息：试样的材料，试样的最大直径D，试样的重量。

另外需要确定检测目的，常见的检测目的有缺陷检测、尺寸测量、密度表征、结构分析等。对于缺陷检测，要确定需要检出的最小缺陷尺寸Defmin。

1.2 转台位置

对于转台位置可以沿射线平行方向移动的系统，源到探测器距离SDD不变，源到转台中心距离SOD是可变的。转台位置确定了有效射束宽度，转台的位置应该在最佳放大倍数Mopt位置，见式(1)。

从式(1)可以看出，d≪a时，Mopt≈1，此时转台的最佳位置越靠近探测器越好，其主要受转台尺寸限制。d≫a时，Mopt很大，此时转台的位置越靠近射线源，分辨率越高，其主要受转台尺寸、试样尺寸和面板尺寸限制，目前微焦面阵探测器CT多工作在此模式下。需要注意的是微焦射线源的焦点尺寸通常随射线源能量的增加而增大，随着射线源能量的增加，最佳放大倍数会减小。

1.3 重建矩阵和重建范围

重建范围(Re)受到试样外形尺寸的影响，通常重建范围可以用式(2)确定。

其即为试样外形最大尺寸的1.5倍。检测时应防止试样未摆放到转台中心，造成CT图像没有完全包含产品。

图1 三代模式下扫描示意

重建矩阵(H×H)用于以一定的精度来量化重建范围，如图1所示(图中圆形区域是重建范围，圆形区域中的小方格表示重建矩阵)，重建范围固定后，重建矩阵越大，分割试样空间越精细，能描述的空间分辨率越高。

像元尺寸p用公式(3)计算，重建矩阵的取值满足式(4)的最接近值。

根据采样定理，通常要求像元尺寸p小于要求检出最小缺陷Defmin的1/2，以保证最小缺陷可以正常地显示。缺陷是否能显示出来，还取决于系统的空间分辨率是否能够优于最小缺陷。

重建矩阵的大小决定了CT图像文件的规模，在达到系统极限空间分辨率图像可正常显示的条件后，即使再增加图像矩阵，也不能再增加细节的分辨率，而且还需要更多的硬盘和内存开销。

另外显示器的物理分辨率也在一定程度上限制了图像细节的显示，当原始图像像元尺寸小于显示器的物理分辨率时，图像显示时会被放大，小的噪声信号也可能被放大到人眼可识别的程度，造成CT图像背景噪声增大，而干扰小缺陷的识别[9-11]；当显示器完全显示所有图像信息时，像素的大小被显示器的物理分辨率替代，细节特征可能会被丢失。

1.4 焦点尺寸

对于焦点尺寸可以选择的X射线源，小焦点具有更高的空间分辨率，但小焦点因为散热问题，功率较小，射线强度弱，探测器采样需要更长的积分时间。对于微焦射线源，为了保证试样的穿透性，需要提高扫描电压，焦点同时也变大。因此应在保证穿透性和射线强度的前提下，尽量使用小焦点。

1.5 扫描电压和扫描电流

扫描电压决定射线的能量，ISO/DIS 15708-1Non-destructiveTesting-RadiationMethods-ComputedTomography-Part1:Principle,EquipmentandSamples中从信噪比的角度给出了穿透试样的射线强度占入射射线强度的10%～20%时，信噪比最好的结论。当穿透射线强度太低，会造成噪声增加，信噪比降低；当穿透射线强度太高时，信号对比度降低，信噪比也降低。也就是说大约穿透2～3个半值层时，信噪比最佳。

扫描电流确定射线强度，在探测器的线形阶段，射线强度越高，信噪比越大。因此在保证探测器不饱和状态下，扫描电流应尽可能大一些。

1.6 采样幅数

对于三代扫描模式，在采样过程中，探测器通道间隔代表了试样空间上的径向采样频率，而采样幅数n表示了圆周方向的采样频率。按照采样定理，n应满足式(5)，以避免最大圆周方向上因采样不足造成噪声增大的问题。

式中：M为放大倍数。

从公式(5)可推导出公式(6)。

而探测器通道间隔受到了探测器硬件成本的限制，因此很多CT系统通过机械系统微动，增加扫描次数来增加径向采样频率，提高空间分辨率。 其是用增加检测时间来换取检测精度的方法。

1.7 切片厚度

对于线阵探测器，切片厚度h由后准直器水平宽度确定。切片厚度决定了z轴方向的分辨率，也影响密度分辨率，切片厚度越小，信噪比越小，图像噪声增大，但对薄型缺陷更容易检出，z方向分辨率越高。从缺陷定量角度分析，切片厚度h满足公式(7)。

1.8 积分时间

积分时间的设置主要与射线扫描电流(射线强度)以及扫描时间相关，积分时间的设置要保证探测器采样信号不饱和。适当增加积分时间可以提高信噪比，但在射线强度严重不足时，仅增加积分时间会导致噪声的增加。因此积分时间的确定要综合考虑探测器响应参数、射线强度等因素。

2 实例分析

2.1 试样信息

采用小孔试片进行实例分析。其材料为不锈钢；试样的规格为：外径50 mm，厚度1.2 mm，其上加工有直径分别为0.1，0.3，0.5，0.7，0.9，1.1 mm的6类人工通孔，每类沿径向分布3个。

2.2 检测设备

采用固鸿科技生产的IPT6110 6 MeV高能工业CT系统进行检测。

设备的加速器能量为6 MeV；焦点尺寸为标称2 mm；探测器608通道，垂直准直器开口0.3 mm，通道间隔1.3 mm，水平准直器0.25～5 mm可调。源到探测器距离(SDD)为3 421.6 mm；源到转台中心距离(SOD)为2 639.4～3 019.4 mm。

2.3 检测要求

可发现0.5 mm以上人工孔。

2.4 工艺参数确定

(1) 将参数代入式(1)，得到Mopt=1.03；受机械系统限制，选择SOD=3 019.4 mm，此时M=1.13。

(2) 将参数代入式(2)，重建范围Re=75 mm。

(3) 将参数代入式(4)，重建矩阵H>320，因系统中最小选项为1 024,故重建矩阵选择1 024×1 024。

(4) 该系统为单焦点，焦点尺寸为2.0 mm。

(5) 扫描电压及扫描电流：该系统射线源为6 MeV加速器系统，因此射线能量只能为6 MeV。

(6) 将参数代入式(6)，得到采样幅数n>1 135.4，因此采样幅数选择满足要求的最小值2 048。

(7) 按式(7)，切片厚度h设置为0.5 mm。

(8) 微动次数是IPT6110设备专用工艺参数，其用扫描时间补充了探测器数量。系统设置选项为2,5,10次，5次微动后可以保证采样频率小于0.2 mm[由(1.3-0.3)/5计算得出],故选择5次微动。

(9) 转台转速为IPT6110设备专用工艺参数，其与采样幅数及加速器触发频率有关，在加速器触发频率固定时，采样幅数与转台转速成反比。此次扫描加速器触发频率依据设备当前条件，设置为170 Hz，因此转台转速为5 r·min-1。

(10) 扫描部位通过在DR图像上的实测来确定，扫描层数为1层。

2.5 检测结果

按以上工艺参数，使用IPT6110设备对试样进行三代CT扫描，检测结果如图2(a)所示，在其上可清晰发现3个0.5 mm孔。

2.6 讨论

CT图像质量通常采用空间分辨率，对比度灵敏度和伪像3个指标来评价。文章使用小孔模拟缺陷的讨论主要针对空间分辨率指标。

图2 试样CT扫描图像及其增加采样幅数和微动次数后的CT图像

实际上，缺陷尤其是小缺陷能否检测出来，还很大程度地依赖于系统的对比度灵敏度。通过增加采样幅数、微动次数、切片厚度、扫描电流等降低噪声、增加对比度灵敏度的手段对小缺陷的检出是比较有利的。图2(b)是增加采样幅数和微动次数后的CT图像，在其上可清晰发现0.3 mm人工孔，但检测时间增加了5倍，可见该检测质量是靠牺牲检测效率来得到的。

3 结论

从检测结果图像可以清晰地发现0.5 mm人工模拟缺陷，验证了检测工艺的正确性，证明了基于采样理论推导出的工艺参数取值公式是可行的。

[1] 张祥春,张鹭,王俊涛.工业CT技术在航空发动机单晶叶片壁厚测量中的应用[J].无损检测,2015,37(2):20-22.

[2] 蒋绍青,栾传彬,满月娥,等. 工业CT在大型复杂机匣检测中的应用[J].无损检测2017,39(2):18-21.

[3] ASTM E 1441-2011 Standard guide for computed tomography (CT) imaging[S].

[4] ASTM E 1570-2011 Standard praactice for computed tomography (CT) examination [S].

[5] ISO/DIS 15708-1-2016 Non-destructive testing-radiation methods-computed tomography-part1: principle,equipment and samples[S].

[6] GJB 5312-2004 工业射线层析成像(CT)检测[S].

[7] GB/T 29070-2012 无损检测 工业CT 检测通用要求[S].

[8] GB/T 29071-2012 无损检测 火工装置的工业CT 检测方法[S].

[9] 何凯，晏磊, 赵红颖.基于人眼视觉分辨力的数字影像质量评价方法研究[J].北京邮电大学学报,2005,28(S1):127-131.

[10] 曾明，张建勋,王湘晖,等.基于视觉特性和复杂度加权处理的图像增强新算法[J].光电子激光,2005,16(3):363-367.

[11] KARUNSEKERA S A, KINGSBURY N K. A distortion measure for blocking artifacts in image based on human visual sensitivity[J].IEEE Transactions on Image Processing,1995,4(6):173-176.

QuantitativeMethodofMainProcessParametersinIndustrialCTDetection

ZHANG Weiguo, NI Peijun, WANG Xiaoyan, QI Zicheng, TANG Shengmin

(The Ningbo Branch of Ordnance Science Institute of China, Ningbo 315103, China)

As an advanced nondestructive testing technology, the application of industrial CT has become increasingly widespread. There are hundreds of industrial CT systems currently in use. If the industrial CT hardware equipment is different, the corresponding process parameters are also very different. Although the current industrial CT standards in the main process parameters provide guide lines for selection principles, they do not give a quantitative selection method and hence the operation is inconvenient. Based on the sampling theory, this paper gives the formula of the main parameters, which greatly facilitates the preparation of the CT detection process. The operation procedure is given by examples and the correctness of the method is verified.

ICT; process parameter; quantitative

TG115.28

A

1000-6656(2017)12-0007-03

2017-06-01

国家自然科学基金资助项目(61471411)；浙江省自然科学基金资助项目(LQ15E010003)

张维国(1972-)，男，研究员，主要从事工业CT等无损检测方法和研究工作

张维国，zwg1972@163.com

10.11973/wsjc201712002