核燃料芯块的表面裂纹检测算法研究 下载: 936次
1 引言
截至2015年,中国核能发电量为1689.93亿千瓦时,约占全国发电总量的3.01%,而核电利用率最高的法国,约占73.3%。可见,我国核电领域仍有很大的发展空间[1]。核燃料芯块作为反应堆的基础核心元件,会经历压制、烧结、研磨等复杂的工艺过程,在装棒之前必须确保燃料芯块外观无瑕疵,以杜绝芯块和包层发生相互作用[2]。燃料芯块破损可能引起反应堆结构安全、热工安全、辐射屏蔽系统承受能力以及工作人员辐射安全等问题[3]。目前主要采用的无损检测手段包括中子照相[4]、DR/CT成像[5]等,虽然能够对芯块内部和外部进行整体的质量检测,但当检测任务仅仅针对芯块表面时[6],存在成本高、技术难度大、检测效率低的问题。因此研究更高效的燃料芯块的表面缺损检测方法,特别是裂纹检测方法,极为重要。
在机器视觉领域,针对其他物体的裂纹检测算法已较为成熟。如文献[
7]中将张量投票算法用于裂纹片段的连接,有效克服了复杂背景下裂纹检出不连续的问题,但该算法运行速度较慢且易受伪裂纹影响;文献[
8]中提取图像的LBP特征,结合SVM对图像是否含有裂纹进行分类,该方法在图像背景较为复杂的情况下效果不好,且提取图像局部二值模式(LBP)特征耗时较长;文献[
9]中利用卷积神经网络(CNN)对道路表面病害进行识别,克服了路面裂纹检测中复杂多变的光照条件,但仅适用于尺寸较大的裂纹;文献[
10]中将Beamlet变换用于图像中线性特征的提取,能够在强噪声的背景中提取线特征,但算法计算速度较慢,无法用于实时检测。同时,燃料元件为黑灰色圆柱体,且尺寸较小,标准规格为
针对上述情况,本文提出基于CNN和Beamlet变换的算法。该算法采用CNN和窗口滑动技术识别图像中的含裂纹区域,去除伪裂纹和背景的干扰;针对含裂纹区域采用Beamlet算法提取裂纹特征;根据形态学判断去除干扰,得到裂纹检测结果。由于真实芯块具有放射性,不便于在实验室长期使用。在系统算法设计和验证阶段,采用中国原子能科学研究院提供的高仿真类陶瓷仿体作为研究和检测对象。实验结果表明,相比于文献[
7]中方法,本文算法的
2 算法实现
本文算法流程如
2.1 图像预处理
待检燃料仿体端面图像如
图 2. 预处理过程。 (a)原始图像; (b)掩模区域; (c)感兴趣区域图像
Fig. 2. Pretreatment process. (a) Original image; (b) mask region; (c) region of interest
得到待处理图像后,按照60 pixel×60 pixel的窗口、30 pixel步长对其进行分割,如
最后,按照1∶1的比例,从所有分割图像块中挑选出无裂纹图像块和含裂纹图像块,将无裂纹图像块标记为0,含裂纹图像块标记为1,作为CrackCNN的数据集,共得到12000组标记数据。为了使CNN能更好地识别微裂纹特征,提高识别准确率,将数据集中所有图像块的大小变换为120 pixel×120 pixel。
2.2 CNN的构建与实现
在裂纹识别中,含裂纹图像块和无裂纹图像块在灰度、形貌、纹理等方面存在差异,可以用CNN对其进行分类,达到裂纹识别的目的。因此根据建立的数据集,按照4∶1的比例将其分为训练集和测试集两种,用于训练和测试搭建CrackCNN模型。CrackCNN的结构设计如
表 1. CrackCNN 模型参数配置
Table 1. Parameter configurations of CrackCNN
|
1) 通过卷积层对输入图像进行卷积运算,提取图像特征,网络中所有卷积核采用高斯分布初始化,标准差设为0.05。卷积计算公式为
式中:
该激活函数能有效解决模型梯度消失和梯度爆炸的问题。
2) 对于卷积层得到的特征图层,为了降低特征维度,提升训练效率,一般采用池化层对其进行降维。最大池化使网络具有一定的平移不变性,可以更好地保留图像裂纹特征,因此选用最大池化进行池化操作。
3) 通过多个卷积池化层,将图像的底层特征转化为高层特征,并经过全连接层对其逐层融合,组成图像的全局特征表述,有效减少了特征信息的损失。由于裂纹分类的互斥性(有无裂纹),网络最后一层神经元数目设定为2,并采用softmax函数输出分类结果的概率分布。同时,为了防止网络过拟合,在全连接层之间采用Dropout方法[12-13]。
采用交叉熵损失函数作为损失函数,通过最小化损失函数,寻找模型最优参数。假定训练集合
式中:
为了使损失函数最小化,采用梯度下降和反向传播算法对模型参数进行更新[14-15],权值矩阵
式中:
在训练过程中,学习率初始值为0.01,每5000次迭代变为原来的0.1倍,bitch_size大小设置为64。网络训练过程中对应的损失函数和测试集上的准确率如
2.3 Beamlet裂纹检测算法
Beamlet算法是由 Donoho等[16]提出的具有局部性、方向性和多尺度特性的图像几何分析工具。该方法在强噪声图像中线特征提取效果较好,但计算量较大。因此本文只针对CrackCNN识别得到的含裂纹区域进行裂纹检测,极大地减少了Beamlet算法的计算量,检测效率大幅度提升。Beamlet算法中有两个重要概念:Beamlet基和Beamlet变换。
对于一幅
图 6. 不同尺度下的 Beamlet 基。(a)尺度为 0;(b)尺度为1;(c)尺度为2;(d)尺度为3
Fig. 6. Beamlets at different scales. (a) Scale is 0; (b) scale is 1; (c) scale is 2; (d) scale is 3
为了方便计算机处理,一个Beamlet基的离散Beamlet变换就是这个基上所有像素点的灰度值加权和。整幅图像中离散Beamlet变换[17]定义为
式中:
提取裂纹特征时,计算每个基的能量统计,然后将能量统计小于阈值的基置0并作为背景,大于阈值的基组成裂纹特征。能量统计公式为
式中:
3 实验结果与分析
核燃料芯块表面裂纹检测中,检测精度和检测效率都非常重要。高精度的检测能够有效降低漏判带来的安全隐患,避免大量误判带来的经济损失;高效率的检测能提高芯块产能。因此,随机挑选30个含裂纹仿体,采用Basler Aca 1600-20gm相机进行端面图像采集,并在Intel(R)Xeon (R) E5-2637 V3@3.5 GHz CPU和NVIDIA Quadro K 2200 GPU平台上,对这30幅图像进行实验,对本文算法性能进行评估。评估指标主要有以下3个[18],查准率(
式中:
3.1 Beamlet变换中本文阈值方法和固定阈值方法结果对比
在选取Beamlet线段时,阈值的设定对检测结果的影响很大。阈值过大可能会引起漏检,过小则可能导致误检,且不同图片之间理想阈值差异较大。为了获得最佳的阈值效果,本文选取一系列固定阈值对30幅含裂纹图像进行处理,即每一个阈值均对30幅图像进行裂纹检测,并计算裂纹检测结果的平均
3.2 裂纹检测结果分析
将本文算法与双重阈值、张量投票算法[7]、传统Beamlet算法[10]进行对比。实验过程中,CrackCNN模型的建立和测试使用的数据集共12000个样本,其中训练集为9600个,测试集为2400个,该数据集中不含用于算法性能评估的30幅图像。模型在10000步后测试集准确率达到97%并趋于稳定,选取此时模型参数用于算法性能分析。采用的Beamlet算法的尺度因子
表 2. 不同方法对比
Table 2. Comparison of different methods
|
从
4 结论
针对传统检测方法中伪裂纹误检以及噪声干扰影响较大的问题,提出CNN与Beamlet相结合的算法,并在燃料芯块仿体上对算法性能进行测试分析。实验结果表明,本文算法有效克服了伪裂纹区域和噪声干扰的问题,具有较高的裂纹检出率,大幅提高了Beamlet算法的检测效率,且在仿真生产线上进行验证,满足相关技术指标,为真实核燃料的表面裂纹检测提供了重要参考。但采用Beamlet算法对裂纹区域检测时,由于在一个子图像块中可能提取多个基,检出裂纹的宽度往往略大于真实裂纹。后续需要对该问题进行深入研究,进一步提高算法性能。
[1] 许建明. 中国核电发展现状与国际化发展对策研究[J]. 电力设备管理, 2018, 2: 1-2.
Xu J M. Research on the status of China's nuclear power development and international development countermeasures[J]. Power Equipment Management, 2018, 2: 1-2.
[2] Piro M HA, SunderlandD, LivingstoneS, et al. A review of pellet-clad interaction behavior in zirconium alloy fuel cladding[M] ∥Sereni J G R. Reference module in materials science and materials engineering. Amsterdam: Elsevier, 2017: 1- 68.
[3] 林晓玲. 建立反应堆燃料元件破损运行判据的思考[J]. 核动力工程, 2013, 34(2): 126-128.
Lin X L. A scheme for establishing of criterions for reactor safe operation in condition of fuel clad failure[J]. Nuclear Power Engineering, 2013, 34(2): 126-128.
[4] 魏国海, 陈东风, 韩松柏, 等. 核燃料元件三维中子成像定量无损检测方法研究[J]. 核动力工程, 2017, 38(S1): 71-73.
Wei G H, Chen D F, Han S B, et al. Study on quantitative non-destructive test method for nuclear fuel rods by three dimensions neutron images[J]. Nuclear Power Engineering, 2017, 38(S1): 71-73.
[5] 涂修佳. 高温气冷堆球形燃料元件质量DR/CT检测与系统开发[D]. 北京: 北京航空航天大学, 2013: 32- 49.
Tu XJ. DR/CT testing system development for spherical fuel elements used in high temperature gas-cooled reactor[D]. Beijing:Beihang University, 2013: 32- 49.
[6] de Oliveira A F, Lapa C M F, et al. . Automatic defect identification on PWR nuclear power station fuel pellets[J]. Nuclear Engineering and Design, 2012, 245: 62-77.
[7] 李慧娴, 张斌, 刘丹, 等. 基于双重阈值和张量投票的表面裂纹检测算法[J]. 激光与光电子学进展, 2018, 55(5): 051010.
[8] Chen F C, Jahanshahi M R, Wu R T, et al. A texture-based video processing methodology using Bayesian data fusion for autonomous crack detection on metallic surfaces[J]. Computer-Aided Civil and Infrastructure Engineering, 2017, 32(4): 271-287.
[9] 沙爱民, 童峥, 高杰. 基于卷积神经网络的路表病害识别与测量[J]. 中国公路学报, 2018, 31(1): 1-10.
Sha A M, Tong Z, Gao J. Recognition and measurement of pavement disasters based on convolutional neural networks[J]. China Journal of Highway and Transport, 2018, 31(1): 1-10.
[10] 李庆嵘. 一种基于Beamlet变换的线特征提取方法[J]. 弹箭与制导学报, 2009, 29(3): 96-98.
Li Q R. A method of line feature extraction based on Beamlet transform[J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2009, 29(3): 96-98.
[11] ZhangL, YangF, Daniel Zhang Y M, et al. Road crack detection using deep convolutional neural network[C]∥2016 IEEE International Conference on Image Processing (ICIP), September 25-28, 2016, Phoenix, AZ, USA. New York: IEEE, 2016: 3708- 3712.
[12] Sutskever I, Hinton G, Krizhevsky A, et al. Dropout: a simple way to prevent neural networks from overfitting[J]. Journal of Machine Learning Research, 2014, 15: 1929-1958.
[13] 曲磊, 王康如, 陈利利, 等. 基于RGBD图像和卷积神经网络的快速道路检测[J]. 光学学报, 2017, 37(10): 1010003.
[14] Cha Y J, Choi W, Büyüköztürk O. Deep learning-based crack damage detection using convolutional neural networks[J]. Computer-Aided Civil and Infrastructure Engineering, 2017, 32(5): 361-378.
[15] 王民, 刘可心, 刘利, 等. 基于优化卷积神经网络的图像超分辨率重建[J]. 激光与光电子学进展, 2017, 54(11): 111005.
[16] Donoho DL, Huo XM. Beamlets and multiscale image analysis[M] ∥Barth T J, Chan T, Haimes R. Multiscale and multiresolution methods. Lecture notes in computational science and engineering. Berlin, Heidelberg: Springer, 2002, 20: 149- 196.
[17] DonohoD, Huo XM. Applications of beamlets to detection and extraction of lines, curves and objects in very noisy images[C]∥Proceedings of Nonlinear Signal and Image Processing(NSIP), Baltimore. [S.l. : s. n.], 2001.
[18] 徐威, 唐振民, 吕建勇. 基于图像显著性的路面裂缝检测[J]. 中国图象图形学报, 2013, 18(1): 69-77.
Xu W, Tang Z M, Lü J Y. Pavement crack detection based on image saliency[J]. Journal of Image and Graphics, 2013, 18(1): 69-77.
Article Outline
宋文豪, 张斌, 李峰宇, 杨腾达, 李建宁, 杨小会. 核燃料芯块的表面裂纹检测算法研究[J]. 激光与光电子学进展, 2019, 56(16): 161008. Wenhao Song, Bin Zhang, Fengyu Li, Tengda Yang, Jianning Li, Xiaohui Yang. Surface Crack Detection Algorithm for Nuclear Fuel Pellets[J]. Laser & Optoelectronics Progress, 2019, 56(16): 161008.