双飞秒激光频率梳光谱测量技术研究进展

高宇炜, 方守龙, 武腾飞, 张磊

高宇炜, 方守龙, 武腾飞, 张磊. 双飞秒激光频率梳光谱测量技术研究进展[J]. 应用光学, 2021, 42(1): 157-175. DOI: 10.5768/JAO202142.0107003
引用本文: 高宇炜, 方守龙, 武腾飞, 张磊. 双飞秒激光频率梳光谱测量技术研究进展[J]. 应用光学, 2021, 42(1): 157-175. DOI: 10.5768/JAO202142.0107003
GAO Yuwei, FANG Shoulong, WU Tengfei, ZHANG Lei. Research progress of double femtosecond laser frequency comb spectroscopy measurement technology[J]. Journal of Applied Optics, 2021, 42(1): 157-175. DOI: 10.5768/JAO202142.0107003
Citation: GAO Yuwei, FANG Shoulong, WU Tengfei, ZHANG Lei. Research progress of double femtosecond laser frequency comb spectroscopy measurement technology[J]. Journal of Applied Optics, 2021, 42(1): 157-175. DOI: 10.5768/JAO202142.0107003

双飞秒激光频率梳光谱测量技术研究进展

基金项目: 航空基金(KJJ82004580)
详细信息
    作者简介:

    高宇炜(1997−),男,硕士研究生,主要从事双光梳光谱测量方面的研究。E-mail:1289725423@qq.com

  • 中图分类号: TN249;O433.1

Research progress of double femtosecond laser frequency comb spectroscopy measurement technology

  • 摘要: 双光梳光谱分析技术近年来凭借高分辨率、高灵敏度、宽光谱覆盖和快速测量的优势,在分子和原子光谱分析领域得到了快速发展。从双光梳光谱分析的原理出发,以使用泵浦光源数量为分类依据,对近年来国内外基于不同实现方案的多种双光梳光谱分析技术的原理进行了简要介绍,叙述了在不同数量泵浦光源条件下,不同实验小组提出的多种实验方案,介绍了他们进行的光谱实验结果。随着研究的不断深入,减少使用泵浦光源数量是双光梳光谱仪实现设备小型化的发展方向之一,增多光源则是从一维走向多维的重要发展思路,经典的双光源配置在测量精度和二维化方面起到了重要作用。
    Abstract: With the advantages of high resolution, high sensitivity, wide spectral coverage and rapid measurement, the dual-comb spectroscopy technology has developed rapidly in the field of molecular and atomic spectroscopy. According to the principle of dual-comb spectroscopy and taken the number of pump light sources used as the classification criteria, the principles of multiple dual-comb spectroscopy techniques are briefly introduced basen on different implementation schemes at home and abroad in recent years, the various experimental schemes were described proposed by different experimental groups with different numbers of pump light sources, and the results of their spectral experiments were introduced. With the deepening of research, reducing the number of pump light sources used is one of the important directions for the development of dual-comb spectrometers to achieve miniaturization of equipment, and the increase is an important development idea from one-dimensional to multi-dimensional. The classic dual light source configuration has made important progress in accuracy and two-dimensionalization.
  • 核电站燃料组件受水下工作环境、装配应力和热应力等因素的影响,在工作中产生局部变形等问题,因此,核电站水下核燃料棒的高精度检测与维护是其安全运行的重要保障[1]。由于水下检测条件恶劣,水下核燃料棒高精度三维测量面临着巨大的挑战。基于线激光扫描的单目视觉测量方法采用单目摄像机在被测物表面做左右方向的扫描运动,获取线激光条纹信息,通过提取激光条纹中心,将其转换为三维点云坐标,实现被测物尺寸的高精度测量。水下检测中的激光条纹照射质量、水下检测环境和待测物表面特性都会对摄像机采集的光条图像造成影响[2],从而给中心线的提取带来一系列问题。因此,如何精确快速地提取光条中心线,是保证整个检测系统测量精度的关键一步。

    近年来,国内外学者提出了许多条纹中心提取方法[3-4],主要有几何法[5]、极值法[6]、骨架提取法[7]、高斯拟合法[8]、方向模板法[9]、灰度重心法[10]等。极值法运算速度快,但若光条出现多个极值,定位精度会明显降低;骨架提取法方法简单,但其精度受环境影响大。田庆国[11]提出基于轮廓线多边形提取方法,可快速定位光带横截面方向并进行中心提取;方向模板法稳定性较弱,中心提取效果一般,王泽浩[12]提出自适应方向模板法,提高了方向模板的提取稳定性;高世一[13]提出变边限高斯拟合提取激光条纹中心线法,可适应条纹宽度进行高斯拟合;吴家勇[14]针对灰度梯度特性,提出了梯度重心法,有效提取了光条中心;Steger[15]基于Hessian矩阵得到光条横截面上一阶导数零点位置,求得光条纹中心线的亚像素位置,其考虑了光条几何特性,定位精度高,但无法解决光条中心多极值的问题,且运算速度慢,无法满足实时测量需求;灰度重心法考虑了光条灰度特性,运算快,熊会元[16]提出应用灰度梯度计算光条法线方向,在法线方向上应用灰度重心法精确提取光条中心的方法,解决了被测物体表面曲率变化大的问题;程鑫[17]结合Hessian矩阵对光条的中心进行定位,利用高斯曲线准确拟合出激光条的中心。以上算法提高了光条中心提取的速度及鲁棒性,但对于提取水下高反光表面光条中心的问题未能很好解决。

    基于上述分析,本文提出了一种反光表面条纹自适应中心提取算法,基于光条几何特性自适应修正高反光区域轮廓,结合图像灰度卷积处理方法自适应修正高反光区域灰度分布,经灰度重心法确定中心点的亚像素坐标,实现了水下核燃料棒表面激光条纹中心的高精度提取。该方法可以有效地减少水下噪声及核燃料棒反光影响,准确、快速地提取激光条纹中心,控制三维重建精度在0.2 mm以内,满足核燃料棒变形测量精度0.3 mm的需求。

    当激光器投影到物体表面,将形成一定宽度的激光条纹,其光强从条纹中心展开,满足高斯分布。其条纹截面灰度值分布曲线的表达式为

    $$ G\left(x\right) = \dfrac{S}{\sqrt{2{\text{π}} \sigma }}·\exp\left[-\dfrac{{\left(d-\mu \right)}^{2}}{2{\sigma }^{2}}\right]$$ (1)

    式中:S为灰度值的幅值;σ为标准差;μ为光条截面坐标;d为光条中心坐标。

    因激光器自身质量、测量环境及被测量物体材质的不同,单目摄像机所采集的激光条纹灰度分布差异性较大,为此,截取核燃料棒原始图像中激光条纹部分,绘制光条左部和中部截面的灰度分布曲线,见图1所示。

    图  1  原始激光条纹及其截面灰度分布
    Figure  1.  Original laser stripes and section gray distribution

    根据图1分析可以得到原始激光条纹图像具有以下特征:1) 光条边界附近存在水下环境噪声干扰;2) 光条截面灰度极值不同,截面宽度存在差异;3) 光条中间部分存在高反光区域,轮廓改变,灰度分布不满足近似高斯分布。

    通过以上特征可以得到,受水下测量环境影响,获取的激光条纹图像存在大量噪声,影响激光条纹分离提取;受物体表面反光特性影响,光条中部截面处存在多灰度极值的情况,靠近光条中心灰度值迅速增加;受金属镜面反光特性影响,条纹中心出现垂直方向的反光噪点,通过调整图片曝光时间及拍摄角度可以减轻;通过传统提取中心算法难以获得精确激光条纹中心,中心线提取依赖光条中灰度分布符合高斯分布的截面信息。

    文中出现较多物理量,主要分为3类。第1类是像素灰度值,其中ω表示像素点经预处理的灰度值;ω′表示像素点经自适应阈值提取的灰度值;ω″表示像素点经过灰度修正的灰度值。第2类是灰度分布,其中Tn表示灰度修正卷积模板;Gij)表示光条所在矩形内的点集。第3类是光条几何特征,其中wn表示激光条纹宽度;kv表示光条轮廓的切线方向;kn表示光条的截面方向。

    实际图像中因其他因素产生的噪声,需要通过图像去噪将其消除。通过对图像中噪声研究发现,水下核燃料图像中的噪声主要分为两种:一是孤立噪点或者由几个噪点组成的假光点;二是沿核燃料棒上下边缘的线状反光噪声。通过滤波方法及卷积运算消除噪声,保证对激光条纹区域提取的精确性。

    1) 首先利用中值滤波处理图像,通过在噪点周围定义一个5×5的邻域窗口,对窗口中像素点灰度值进行升序排列,取灰度值的中值作为滤波结果,可以去除水下核燃料棒激光条纹图像中的椒盐噪点。

    2) 对原始图像进行全局阈值的二值化处理,根据原始图像选择合适的灰度阈值Ct,计算实现图像二值化:

    $$\omega (\mu ,\nu ) = \left\{ {\begin{array}{*{20}{c}} {1,\omega \left( {\mu ,\nu } \right) \geqslant {C_{{t}}}} \\ {0,\omega \left( {\mu ,\nu } \right) \leqslant {C_{{t}}}} \end{array}} \right.$$ (2)

    式中ωμν)表示二值化后的光条在坐标(μν)处的灰度值。

    在灰度值非零的点上利用3×3模板去除噪点,模板如图2所示,a1a8分别对应图像中的每个像素的位置。

    图  2  卷积模板
    Figure  2.  Convolution template

    将该卷积模板在图像上移动,根据(3)式计算判定该点是否为噪点:

    $$A = \sum\limits_{r = 0}^8 {{{\left( {1 - {a_r}} \right)}^2}} $$ (3)
    $$\omega (\mu ,\nu ) = \left\{ {\begin{array}{*{20}{c}} {1,A \leqslant 4} \\ {0,A > 4} \end{array}} \right.$$ (4)

    通过对点(μν)的A值进行判断,去除图像中的孤立噪点和沿核燃料棒上下边缘的反光噪点。将去噪图像与原始图像各像素灰度值相乘,获得具有原始灰度信息的去噪激光图像,最后使用高斯滤波平滑图像。

    3) 根据图像灰度值,实现激光条纹的自适应阈值提取。从垂直激光平面方向开始扫描,寻找包含像素点的灰度值为255的行。当第一次检测到时,记该点为(iminjmin),最后一次检测到时,记该点为(imaxjmax),框选出包含光条位置的矩形,矩阵宽度为wd。则光条所在的矩形内点集Gij)为

    $$ \begin{split} \Bigg\{ G(i,j)\Bigg|& {i_{\min }} \leqslant i \leqslant {i_{\max }},\dfrac{{{j_{\min }} + {j_{\max }} - {w_{\rm{d}}}}}{2} <j < \\ & \dfrac{{{j_{\min }} + {j_{\max }} + {w_{\rm{d}}}}}{2} \Bigg\} \end{split} $$ (5)

    根据每列像素灰度值对该光条截面的不同影响程度确定阈值,选择使用基于梯度的自适应行阈值法在矩形内计算每列的阈值。如图3所示,将激光条纹图像中任一列的灰度分布抽象为两条一次函数:

    图  3  自适应阈值原理
    Figure  3.  Principle of self-adaptive threshold

    计算两像素点间灰度梯度Vmn

    $${V_{mn}} = \left| {{\omega _{mn}} - {\omega _{(m - 1)n}}} \right|$$ (6)

    式中ωmn)是点(μν)在图像第mn列的灰度值。

    计算每行灰度梯度的平均值Bn

    $${B_n} = \dfrac{{\displaystyle\sum\limits_{m = 2}^M {{V_{mn}}} }}{{D - 1}}$$ (7)

    式中D为光条大致宽度。

    在进行测量时,D为25个像素。通过下式判断点(μν)是否属于光条有效点:

    $${g_n} = 1 - \dfrac{{{B_n} \times D}}{{{G_{\max }}}}$$ (8)
    $$ {\omega' }(\mu ,\nu )=\left\{\begin{array}{c}\omega (\mu ,\nu ),\omega (\mu ,\nu )\geqslant {g}_{n}·{G}_{\max}\\ 0,\omega (\mu ,\nu )<{g}_{n}·{G}_{\max}\end{array}\right.$$ (9)

    式中:gn表示阈值界限系数;Gmax表示每行灰度的最大值;ω'μν)表示自适应阈值提取的光条在点(μν)处的灰度值。

    该方法能使光条的提取阈值根据列像素自适应变化,快速将激光从环境背景中提取,图4(a)为经自适应分割的光条图像。通过图4(b)分析,光条附近噪声被去除,提取出的光条图像边界明显,保证了提取精度。

    图  4  自适应阈值分割及其激光截面灰度分布图
    Figure  4.  Self-adaptive threshold segmentation and gray distribution of laser section

    提取的光条已从图像背景分割,但由于反光现象,激光条纹在被测物体表面的轮廓发生改变,此时直接进行光条中心线提取无法达到精度要求。为解决这一问题,选取灰度分布满足高斯分布的光条截面端点,采用3层BP神经网络对所选坐标点进行曲线拟合,拟合出两条激光条纹轮廓,对光条进行进一步的修正提取。

    图5所示,BP神经网络设定输入层结点为X,隐藏层节点为h,输出层节点为Y,输入层和输出层节点数为1,隐藏层的节点数C为2。设定wic为输入层到隐层第c个节点之间的连接权值,ψc表示隐层第c个节点的阈值;woc为隐层第c个节点到输出层之间的连接权值,θ表示输出层节点的阈值;隐层输入变量为hicc=1,2,隐层输出变量为hocc=1,2。

    图  5  神经网络的基本原理
    Figure  5.  Basic principle of neural network

    隐层激活函数fh选用(10)式Sigmoid函数,输出层激活函数fo则选用purelin函数。

    $${f_{{h}}}(x) = \dfrac{1}{{1 + \exp \left( { - x} \right)}}$$ (10)

    隐层的第c个节点的输入为

    $${h_{ic}} = {w_{ic}}X - {\psi _c}$$ (11)

    隐层的第c个节点的输出为

    $${h_{oc}} = {f_{{h}}}\left( {{w_{ic}}X - {\psi _c}} \right)$$ (12)

    输出层节点的输入为

    $${y_i} = {w_{oc}}{h_{oc}} - \theta = {w_{oc}}{f_{{h}}}\left( {{w_{ic}}X - {\psi _c}} \right) - \theta $$ (13)

    输出层节点的输出为

    $${y_{{o}}} = {f_{{o}}}\left( {{y_i}} \right) = {f_{{o}}}\left[ {{w_{oc}}{f_{{h}}}\left( {{w_{ic}}X - {\psi _c}} \right) - \theta } \right]$$ (14)

    设定Yk)为第k个点的纵坐标,对图像中提取的总数为K的像素点拟合轮廓,在网络反向传播时神经网络得到的轮廓位置误差为

    $${E^{\left( t \right)}} = \dfrac{1}{2}\sum\limits_{k = 1}^K {{{\left( {{y_{{o}}}^{\left( t \right)}\left( k \right) - Y\left( k \right)} \right)}^2}} $$ (15)

    式中yotk)为神经网络权值在第t次调整后得到的第k个点的纵坐标。

    根据Levenberg-Marquardt最优化算法调整,得到最终神经网络权值。训练直到网络输出的误差减少到给定误差,或进行到预先设定的训练次数时,可以停止训练。设定目标误差为0.01,设定最大训练次数为1 000次,训练数据选取各行两端的非高反光区域像素坐标点,确定其为曲线轮廓拟合的真实值,基于条纹轮廓真实值进行预测,最终得到激光条纹的拟合曲线如图6所示。

    图  6  轮廓拟合
    Figure  6.  Contour fitting

    分别提取光条两边有效像素点进行轮廓拟合,可以获得经过轮廓修正的激光条纹图像,如图7所示。

    图  7  轮廓修正后的激光条纹
    Figure  7.  Laser stripe after contour correction

    图8所示,两组数据分别为轮廓修正前后的激光条纹各行宽度,因反光干扰,修正前条纹各行宽度变化较大。通过BP神经网络对激光条纹的轮廓进行拟合,修正表面反光造成的光条宽度变化的部分,解决了激光条纹轮廓变形问题。

    图  8  轮廓修正前后的光条宽度
    Figure  8.  Width of light bar before and after contour correction

    激光条纹反光部分轮廓已经修正,但各像素灰度值仍未改变。此时若对光条进行中心线提取,反光区域将存在多极值点的情况。本文方法通过设定自适应卷积模板修正光条灰度分布,利用条纹宽度信息生成模板对图像进行卷积运算,获得经灰度分布正常的激光条纹图像。因为线激光照射在被测物体表面后其灰度分布近似符合高斯分布,所以灰度修正系数将通过高斯函数设定。同一列中,像素修正系数由中部向列两边依次递减,激光条纹图像中任意一点p灰度修正系数表示为

    $$\begin{split}& {f}_{p}\left({w}_{n},{m}_{p},{m}_{\min},\sigma \right) = \dfrac{\delta }{\sqrt{2{\text{π}} \sigma }}\cdot \\ &\quad\quad\exp\left\{-\dfrac{{\left[{m}_{p}-\left({m}_{\min}+\dfrac{{w}_{{n}}}{2}\right)\right]}^{2}}{2{\sigma }^{2}}\right\}\end{split}$$ (16)

    式中:wnn=0,1,…,N−1)为经过轮廓修正的激光条纹图像中点p所在第n列的非零灰度值像素个数,即光条带宽;N为含有非零灰度值的总列数,表示激光条纹长度;mp为点p所在行数;mmin为点p所在的n列中第1个非零像素点的行数;б为高斯函数标准差,取值为4;δ为固定系数,取值为0.95。

    因为点p的灰度修正系数已经求得,所以激光条纹上的点p修正灰度值如(17)式所示:

    $$ {\omega'' }\left({\mu }_{p},{\nu }_{p}\right)={f}_{p}\left({w}_{n},{m}_{p},{m}_{\min},\sigma \right)·{\omega '}\left({\mu }_{p},{\nu }_{p}\right)$$ (17)

    式中:ω″(μpνp)为点(μν)经过灰度修正后的灰度值。

    已知任意点的灰度修正系数后,可设计激光条纹图像中任意列的灰度修正卷积模板。卷积模板的行数设为激光条纹宽度wn,列数设为1,此时激光条纹的每个截面都可生成与其截面宽度相适应的灰度修正卷积模板Tnn= 0,1,…,N−1)。由(16)式知fp的范围在(0~1)之内,且根据光条截面宽度wnfp在模板中间时取最高值,然后依次向两边递减。激光条纹图像中任意列n的灰度修正卷积模板为

    $${T_n} = {\left( \begin{array}{c} {f_1} \\ {f_2} \\ \vdots \\ {f_{{M_n} - 1}} \\ {f_{{M_n}}} \\ \end{array} \right)_{{M_n} \times 1}}$$ (18)

    式中Mn为激光条纹第n列含有非零灰度值的总行数。

    根据激光条纹宽度生成自适应灰度修正卷积模板,在光条原始灰度分布矩阵Iβα)上自左向右进行滑动运算,可以得到激光光条修正后的灰度分布:

    $$ {{I'}}(\beta ,\alpha ) = {\left( {\begin{array}{*{20}{c}} 0&\cdots& 0&{{f_{m1}} \cdot G(m,1)}&0&\cdots&0\\ \vdots& & & & & &\vdots \\ 0& &\vdots & \vdots & \vdots &&0\\ {{f_{11}} \cdot G(1,1)}&& & \cdots & & &{{f_{\alpha 1}} \cdot G(\alpha ,1)}\\ \vdots && \cdots &{{f_{mn}} \cdot G(m,n)}& \cdots && \vdots \\ {{f_{1{\beta _1}}} \cdot G(1,{\beta _1})}&& & \cdots & & &{{f_{\alpha {\beta _\alpha }}} \cdot G(\alpha ,{\beta _\alpha })}\\ 0&& \vdots & \vdots & \vdots&&0\\ \vdots& & & & && \vdots \\ 0&\cdots&0&{{f_{m{\beta _m}}} \cdot G(m,{\beta _m})}&0&\cdots&0 \end{array}} \right)_{\alpha \times \beta }} $$ (19)

    式中:α为光条原始图像中含有非零灰度值的总列数;βm为第m列中含有非零灰度值的总行数,且1≤mαβ为激光条纹各列中含有非零灰度值总行数的最大值。

    对经轮廓修正的水下核燃料棒扫描图像进行卷积灰度修正,其结果如图9(a)所示。

    图  9  灰度修正后的激光条纹及其截面灰度分布
    Figure  9.  Laser stripe after gray correction and section gray distribution

    激光条纹反光部分灰度值过高,灰度分布改变,因此通过自适应卷积模板处理激光条纹图像灰度,对中间高反光区域进行灰度值降低使其满足高斯分布规律,如图9(b)所示,与图4(b)对比所见,获得了良好效果。

    采用灰度重心法对完成修正的光条图像进行激光条纹亚像素中心提取。因为核燃料棒是圆柱体,所以线激光平面与核燃料棒表面相交形成的激光条纹轮廓通过单目摄像机采集后成曲线轮廓。对于非直线激光条纹的中心提取,需要在其光条截面方向进行。如图10所示,因为光条轮廓线已确定,所以可以通过计算其轮廓线上任一像素点的切线kv,确定光条的截面方向kn

    图  10  激光截面中心
    Figure  10.  Center of laser section

    取轮廓线上一点为q,通过光条轮廓线在q点的一阶导数可求得该点在光条轮廓上的切线斜率kv。此时易得出经q点的光条截面方向斜率kn

    $${k_n} = - \dfrac{1}{{{k_v}}}$$ (20)

    光条截面方向线条经过q点后向前后方向延长,直至抵达光条边界像素点,自适应光条宽度。通过判断像素点所在的行向量与列向量是否与斜率kn同时相交,确定其是否被光条截面方向线条经过。将符合条件的像素点记为一组像素点集合Xq={xh| xh=(μhνh),h=1,2,…,H},其中H为光条截面方向斜率kn所经过的像素点总数。q点所在截面的亚像素中心通过在该像素点集合Xq内使用灰度重心法求解:

    $$ {q_c}\left( {{x_q},{y_q}} \right) = \left\{ {\begin{array}{*{20}{c}} {{x_q} = \dfrac{{\displaystyle\sum\limits_{h = 1}^H {{\mu _h} \cdot \omega ''({\mu _h},{\nu _h})} }}{{\displaystyle\sum\limits_{h = 1}^H {\omega ''({\mu _h},{\nu _h})} }}}\\ {{y_q} = \dfrac{{\displaystyle\sum\limits_{h = 1}^H {{\nu _h} \cdot \omega ''({\mu _h},{\nu _h})} }}{{\displaystyle\sum\limits_{h = 1}^H {\omega ''({\mu _h},{\nu _h})} }}} \end{array}} \right. $$ (21)

    式中:xqq点所在光条截面的亚像素中心横坐标;yq为纵坐标。

    采用Visio studio 2013作为软件开发平台、结合OpenCV 3.0库及PCL 1.8.0点云库实现文中方法。如图11所示,采用的水下核燃料棒线结构光测量平台,搭载Mako G-507 PoE 工业相机,其分辨率为2 464(H)× 2 056(V)像素,像素尺寸为3.45 μm×3.45 μm;激光器波长为520 nm。通过经水下标定的水下核燃料棒测量系统对核燃料棒尺寸进行测量。

    图  11  水下核燃料棒测量系统
    Figure  11.  Underwater nuclear fuel rod detection system

    图12(a)12(c)分别选用传统灰度重心法、基于最小二乘拟合中心提取算法[18]和文中算法进行激光条纹中心提取。从图12(a)可以看出,传统灰度重心法提取的光条中心线条受反光影响,光条中心点存在前后错开的现象,线条不连续;图12(b)中基于最小二乘拟合中心提取算法对光条进行提取,光条中心线条连续性增加,但由于其提取的中心点与拟合中心点重合不足,存在线条中心不唯一的情况,容易产生噪点;图12(c)为文中算法,算法充分利用了光条的有效数据,光条中心线条连续性较好,更客观地反映了光条中心分布情况。

    图  12  不同激光条纹中心提取方法对比
    Figure  12.  Comparison of different laser stripe center extraction methods

    使用3种算法提取激光条纹中心线后生成点云进行对比。图13(a)点云为传统灰度重心法生成,图13(b)点云为基于最小二乘拟合中心提取算法生成,图13(c)点云为文中算法生成。结构光中心线提取不准确,将导致通过图像坐标点求解出的三维坐标发生较大偏差,产生点云离群现象。可以看出,通过传统灰度重心法和基于最小二乘拟合中心提取算法得到的点云含有较多的离群噪点,被测物体各部分的反光导致多数截面光条中心线条出现断点和噪点,这些会影响点云测量的准确性;通过文中算法解决因反光造成的结构光中心线提取不准的现象,可以明显看到所获得的点云大量噪点被消除,圆柱点云表面更加光滑。

    图  13  点云计算结果
    Figure  13.  Calculation results of point cloud

    通过对10组测量提取的激光条纹中心线生成的点云进行直径测量,与核燃料棒人工采用游标卡尺测量的结果进行比对。表1所示为10组拟合直径统计结果。针对水下核燃料棒检测要求,测量数据应保持高精度与稳定性,传统灰度重心法在水下测量中稳定性较好,但是其测量数据误差较大,基于最小二乘拟合中心提取算法的测量数据的最大误差和平均误差有所减少,但是其稳定性不足。文中算法相较于传统灰度重心法和基于最小二乘拟合中心提取算法平均误差明显降低,分别为97%和89%,测量稳定性分别提高9%和33%,测量效果良好,拟合圆柱的误差在0.2 mm以内。可以看出,文中提出的算法能够有效减少核燃料棒水下测量的噪声,具有较好的稳定性和精确性。

    表  1  核燃料棒测量结果
    Table  1.  Measurement results of nuclear fuel rod mm
    方法灰度重
    心法
    基于最小二乘拟合
    中心提取算法
    文中算法
    测量值9.3979.4439.528
    9.3819.4829.822
    9.7319.5449.636
    9.4899.7639.732
    9.4789.7619.548
    9.6189.6329.586
    9.5829.8879.736
    9.7219.4649.622
    9.4919.8349.612
    9.5109.5529.841
    最大误差0.28900.21700.1710
    平均误差0.13020.03380.0037
    标准差0.12140.16330.1102
    下载: 导出CSV 
    | 显示表格

    针对水下核燃料棒检测过程中激光条纹受水下环境噪声及反光影响的问题,本文根据激光条纹特征提出了结合BP神经网络轮廓拟合与卷积运算的自适应激光条纹亚像素中心提取方法。文中算法通过对核燃料棒反光部分的几何特性与灰度特征进行自适应修正,实现了水下核燃料棒表面激光条纹中心线的高精度提取。实验证明,使用该方法实现了在水下测量650 mm处燃料棒尺寸误差控制在0.3 mm以内。与其他算法进行比较,平均误差分别降低97%和89%,稳定性分别提高9%和33%,验证了该算法测量的精确性和鲁棒性。文中算法对于核燃料棒尺寸的准确测量,保证核燃料组件安全运行具有重要意义。

  • 图  1   双光梳光谱测量原理

    Figure  1.   Measuring principle of dual-comb spectrum device

    图  2   单光梳扫频调谐型DCS方案

    Figure  2.   Structure and measuring principle of single-comb comb frequency-tuning DCS scheme

    图  3   自适应双偏振型双太赫兹梳光谱装置

    Figure  3.   Self-adaptive dual-polarization dual terahertz comb spectrum device

    图  4   $ {\mathrm{C}\mathrm{H}}_{3}\mathrm{C}\mathrm{N} $气体的吸收光谱图

    Figure  4.   Absorption spectrogram for $ {\mathrm{C}\mathrm{H}}_{3}\mathrm{C}\mathrm{N} $

    图  5   双波长锁模型DCS

    Figure  5.   Dual-wavelength mode-locked DCS

    图  6   与HITRAN数据库对比的吸收光谱图

    Figure  6.   Comparison of absorption spectrogram and HITRAN database

    图  7   腔内双向传输型DCS和吸收谱线对比图

    Figure  7.   In-cavity bidirectional transmission DCS and comparison with database

    图  8   芯片化微腔光梳系统

    Figure  8.   Microcavity comb DCS

    图  9   并联型双芯片DCS的光频及对应拍频吸收强度

    Figure  9.   Optical frequency of parallel dual chip DCS and corresponding beat frequency absorption intensity

    图  10   光梳对1和2测得的光谱数据与HITRAN对比

    Figure  10.   Comparison of optical data of optical comb pair 1 and 2 with HITRAN

    图  11   芯片化双微腔光梳系统

    Figure  11.   Microchip double microcavity comb system

    图  12   二氯甲烷($ {\mathrm{C}\mathrm{H}}_{2}{\mathrm{C}\mathrm{l}}_{2}) $的液相吸收光谱实验

    Figure  12.   Dichloromethane absorption spectrum device

    图  13   马-曾调制器DCS方案及其输出干涉谱线图

    Figure  13.   DCS scheme of Mach-Zehnder modulator and its output interference spectrum

    图  14   光频参考锁定光梳结构示意图

    Figure  14.   Configuration scheme of OFC locked with optical frequency reference

    图  15   典型双光梳光谱测量系统

    Figure  15.   Typical double comb spectrometric system

    图  16   采用光纤激光器的双光梳光谱测量系统测试结果

    Figure  16.   Test results of double comb spectrum measurement system with fiber laser

    图  17   双光梳锁定偏置频率方案示意图

    Figure  17.   Dual-comb offset frequency locking scheme

    图  18   偏频锁定方案吸收谱线图

    Figure  18.   Absorption spectrum of offset frequency locking scheme

    图  19   自适应DCS装置示意图(a)及其吸收谱线图(b)

    Figure  19.   Schematic diagram of self-adaptive DCS device and its absorption spectrum

    图  20   自适应DCS示意图(a)及其吸收谱线与HITRAN数据库对比图

    Figure  20.   Schematic diagram of self-adaptive DCS and comparison of absorption line and HITRAN database

    图  21   自适应双光梳光谱实验装置图

    Figure  21.   Self-adaptive dual optical comb spectrum experiment device

    图  22   二维双光梳光谱原理图

    Figure  22.   Two-dimensional dual-comb spectroscopy principle

    图  23   FWM信号频谱图

    Figure  23.   FWM signal spectrum

    图  24   自适应二维双梳方案示意图

    Figure  24.   Schematic diagram of self-adaptive two-dimensional dual-comb scheme

    图  25   三光梳二维相干光谱方案图

    Figure  25.   Tri-comb two-dimensional coherent spectrum scheme diagram

    图  26   光频域和射频域频率分布图

    Figure  26.   Frequency distribution diagram of optical frequency domain and RF domain

    图  27   铷原子不同同位素的谱线测量结果

    Figure  27.   Spectral line measurement results of different isotopes of rubidium atom

    表  1   各实验小组采用方案指标对比

    Table  1   Comparison table of scheme indicators adopted by each experimental group

    采用方案类型光谱分辨率信噪比/dB光谱宽度参考文献
    扫频调谐型 0.5 GHz 350(2 500 s) 8 THz [26]
    双偏振型 5 MHz 50(200 ns) 1.5 THz [27]
    双波长型 188 MHz 100(166 ms) 15 nm [28]
    双向输出型 93 GHz 40 nm [29]
    并联微腔 319±4 kHz 28(20 μs) 2.3 THz [30]
    并联微腔 10.453 GHz 50(1 s) 1.7 THz*2 [31]
    串联微腔 1.4±0.8 GHz 40(500 μs) 51 THz [32]
    马曾调制器 4 MHz 25 GHz [33]
    双精密锁定 1.1 kHz 35 3 THz [34]
    偏频锁定光频 0.086 cm−1 200(62.5 ms) 16.5 GHz [36]
    自适应双梳 1.1 GHz 20(467 μs) 14.5 THz [14]
    自适应双梳 1.09 GHz 350(3 s) 20 nm [38]
    下载: 导出CSV
  • [1] 魏志义, 张炜, 韩海年. 光学频率梳的研究进展和应用[C]//2008年全国频率控制技术年会论文集. 北京: 中国宇航学会, 中国电子元件行业协会, 2008: 108-112.

    WEI Zhiyi, ZHANG Wei, HAN Hainian. Research progress and application of optical frequency comb[C]//Proceedings of the 2008 National Annual Conference on Frequency Control Technology. Beijin: Chinese Aerospace Society, Chinese Electronic Components Industry Association, 2008: 108-112.

    [2]

    UDEM T, HOLZWARTH R, HÄNSCH T W. Optical frequency metrology[J]. Encyclopedia of Materials Science & Technology,2005,416(6877):1-5.

    [3] 孙青, 杨奕, 孟飞, 等. 基于频率梳的太赫兹频率精密测量方法研究[J]. 光学学报,2016, 36(4): 133-138.

    SUN Qing, YANG Yi, MENG Fei, et al. Research on precise measurement method of terahertz frequency based on frequency comb[J]. Acta Optics Sinica, 2016,36 (4): 133-138.

    [4]

    CHOU C W, HUME D B, ROSENBAND T, et al. Optical clocks and relativity[J]. Science, 2010,329(5999): 1630-1633.

    [5]

    HOLZWARTH R, UDEM T, HAENSCH T, et al. Optical frequency synthesizer for precision spectroscopy[J]. Phys. Rev. Lett.,2000,85(11):2264-2267. doi: 10.1103/PhysRevLett.85.2264

    [6]

    CODDINGTON I, SWANN W C, NENADOVIC L, et al. Rapid and precise absolute distance measurements at long range[J]. Nature Photonics,2009,3(6):351-356. doi: 10.1038/nphoton.2009.94

    [7]

    MANCE J G, La LONE B M, DOLAN D H, et al. Time-stretched photonic Doppler velocimetry[J]. Optics Express,2019,27(18):25022-25030. doi: 10.1364/OE.27.025022

    [8]

    NEWBURY N R, CODDINGTON I, SWANN W C. Optical frequency comb-based coherent LIDAR, US: 8558993[P]. 2013-10-15.

    [9]

    ADLER F, MASŁOWSKI P, FOLTYNOWICZ A, et al. Mid-infrared Fourier transform spectroscopy with a broadband frequency comb[J]. Optics Express,2010,18(21):21861. doi: 10.1364/OE.18.021861

    [10]

    GOHLE C, STEIN B, SCHLIESSER A, et al. Frequency comb vernier spectroscopy for broadband, high-resolution, high-sensitivity absorption and dispersion spectra[J]. Physical Review Letters,2007,99(26):263902-1-4.

    [11]

    CODDINGTON I, NEWBURY N, SWANN W. Dual-comb spectroscopy[J]. Optica,2016,3(4):414. doi: 10.1364/OPTICA.3.000414

    [12]

    SCHILLER S. Spectrometry with frequency combs[J]. Optics Letters,2002,27(9):766-768. doi: 10.1364/OL.27.000766

    [13]

    KEILMANN F, GOHLE C, HOLZWARTH R. Time-domain mid-infrared frequency-comb spectrometer[J]. Optics Letters,2004,29(13):1542-1544. doi: 10.1364/OL.29.001542

    [14]

    IDEGUCHI T, POISSON A, GUELACHVILI G, et al. Adaptive real-time dual-comb spectroscopy[J]. Nature Communications,2014(5):3375.

    [15]

    ROY J, DESCHÊNES J D, POTVIN S, et al. Continuous real-time correction and averaging for frequency comb interferometry[J]. Optics Express,2012,20(20):21932-21939. doi: 10.1364/OE.20.021932

    [16]

    CODDINGTON I, SWANN W C, NEWBURY N R. Coherent dual-comb spectroscopy at high signal-to-noise ratio[J]. Physical Review A,2010,82(4):43817. doi: 10.1103/PhysRevA.82.043817

    [17] 赵俊杰, 刘迪, 乔建翔, 等. 关于汽车尾气成分动态检测技术的分析[J]. 内燃机与配件,2019,280(4):137-138.

    ZHAO Junjie, LIU Di, QIAO Jianxiang, et al. Analysis on the dynamic detection technology of automobile exhaust components[J]. Internal Combustion Engine and Parts,2019,280(4):137-138.

    [18] 齐森. 废弃印刷电路板热拆解过程尾气成分分析及处理方法研究[D]. 合肥: 合肥工业大学, 2015.

    QI Sen. Research on components analysis and treatment methods of exhaust gas during heating and dismantling waste printed circuit boards[D]. Hefei: Hefei University of Technology, 2015.

    [19] 张文强. 恶臭气体的检测方法与技术研究[D]. 天津: 河北工业大学, 2014.

    ZHANG Wenqiang. Malodorous detection method and technology research[D]. Tianjin: Hebei University of Technology, 2014.

    [20] 何启欣. 基于红外激光吸收光谱技术的气体检测系统研究[D]. 长春: 吉林大学, 2018.

    HE Qixin. Research on gas detection system based on infrared laser absorption spectroscopy technique[D]. Changchun: Jilin University, 2018.

    [21] 刘文清, 崔志成, 刘建国, 等. 大气痕量气体测量的光谱学和化学技术[J]. 量子电子学报,2004,21(2):202-210. doi: 10.3969/j.issn.1007-5461.2004.02.011

    LIU Wenqing, CUI Zhicheng, LIU Jianguo, et al. Spectroscopy and chemical technology for atmospheric trace gas measurement[J]. Chinese Journal of Quantum Electronics,2004,21(2):202-210. doi: 10.3969/j.issn.1007-5461.2004.02.011

    [22]

    TAKURO I. Dual-comb spectroscopy[J]. Optics & Photonics News,2017,28(1):32.

    [23] 路桥, 时雷, 毛庆和. 双光梳光谱技术研究进展[J]. 中国激光,2018,45(4):7-28.

    LU Qiao, SHI Lei, MAO Qinghe. Research advances in dual-comb spectroscopy[J]. Chinese Journal of Lasers,2018,45(4):7-28.

    [24]

    CODDINGTON I, SWANN W C, NEWBURY N R. Coherent linear optical sampling at 15 bits of resolution[J]. Optics Letters, 2009,34(14): 2153-2155.

    [25] 张宏飞, 苏波, 何敬锁, 等. 超快太赫兹时域光谱系统[J]. 应用光学,2019,40(2):229-232.

    ZHANG Hongfei, SU Bo, HE Jingsuo, et al. Ultra-fast terahertz time domain spectroscopy system[J]. Journal of Applied Optics,2019,40(2):229-232.

    [26]

    LEE K, LEE J, JANG Y S, et al. Fourier-transform spectroscopy using an Er-doped fiber femtosecond laser by sweeping the pulse repetition rate[J]. Scientific Reports,2015(5):15726.

    [27]

    CHEN J, NITTA K, ZHAO X, et al. Comb-mode-resolved adaptive sampling terahertz dual-comb spectroscopy with a free-running single-cavity fiber laser[EB/OL]. [2019-08-03]. https://arxiv.org/abs/1908.01764.

    [28]

    ZHAO X, HU G Q, ZHAO B F, et al. Picometer-resolution dual-comb spectroscopy with a free-running fiber laser[J]. Optics Express,2016,24(19):21833-21845. doi: 10.1364/OE.24.021833

    [29]

    IDEGUCHI T, NAKAMURA T, KOBAYASHI Y, et al. Kerr-lens mode-locked bidirectional dual-comb ring laser for broadband dual-comb spectroscopy[J]. Optica,2016,3(7):748-753. doi: 10.1364/OPTICA.3.000748

    [30]

    LIN T, DUTT A, JOSHI C, et al. Broadband ultrahigh-resolution chip-scale scanning soliton dual-comb spectroscopy[EB/OL]. [2020-01-13]. https://arxiv.org/abs/2001.00869.

    [31]

    SHAMS-ANSARI A, YU M J, CHEN Z J, et al. An integrated lithium-niobate electro-optic platform for spectrally tailored dual-comb spectroscopy[EB/OL]. [2020-03-11]. https://arxiv.org/abs/2003.04533.

    [32]

    DUTT A, JOSHI C, JI X C, et al. On-chip dual-comb source for spectroscopy[J]. Science Advances,2018,4(3):e1701858. doi: 10.1126/sciadv.1701858

    [33]

    DENIEL L, WECKENMANN E, GALACHO D P, et al. Frequency-tuning dual-comb spectroscopy using silicon Mach-Zehnder modulators[J]. Optics Express,2020,28(8):10888-10898. doi: 10.1364/OE.390041

    [34]

    CODDINGTON I, SWANN W C, NEWBURY N R. Coherent, multi-heterodyne spectroscopy using stabilized optical frequency combs[J]. Physical Review Letters,2008,100(1):013902. doi: 10.1103/PhysRevLett.100.013902

    [35]

    COBURN S, ALDEN C B, WRIGHT R, et al. Regional trace-gas source attribution using a field-deployed dual frequency comb spectrometer[J]. Optica,2018,5(4):320. doi: 10.1364/OPTICA.5.000320

    [36] 张伟鹏, 杨宏雷, 陈馨怡, 等. 光频链接的双光梳气体吸收光谱测量[J]. 物理学报,2018,67(9):65-74.

    ZHANG Weipeng, YANG Honglei, CHEN Xinyi, et al. Optical frequency linked dual-comb absorption spectrum measurement[J]. Acta Phys. Sin,,2018,67(9):65-74.

    [37]

    GIACCARI P, DESCHÊNES J D, SAUCIER P, et al. Active Fourier-transform spectroscopy combining the direct RF beating of two fiber-based mode-locked lasers with a novel referencing method[J]. Optics Express,2008,16(6):4347-4365. doi: 10.1364/OE.16.004347

    [38] 杨力, 沈旭玲, 杨康文, 等. 自适应双光梳光谱原理分析与实现[J]. 光学学报,2018,38(5):156-162.

    YANG Li, SHEN Xuling, YANG Kangwen, et al. Analysis and realization of adaptive dual-comb spectroscopy[J]. Acta Optics Sinica,2018,38(5):156-162.

    [39]

    LOMSADZE B, CUNDIFF S T. Multi-heterodyne two dimensional coherent spectroscopy using frequency combs[J]. Scientific Reports,2017,7(1):14018. doi: 10.1038/s41598-017-14537-z

    [40]

    LOMSADZE B, SMITH B C, CUNDIFF S T. Tri-comb spectroscopy[J]. Nature Photonics,2018,12(11):676-680. doi: 10.1038/s41566-018-0267-4

  • 期刊类型引用(5)

    1. 周晟刚. 通风机械仪表盘在复杂背景环境中视觉故障检测与定位研究. 计算机测量与控制. 2024(03): 106-111+145 . 百度学术
    2. 吴来新,王志. 基于机器学习的眼位测定仪视线跟踪误差自动复合补偿技术. 自动化与仪器仪表. 2024(10): 370-374 . 百度学术
    3. 洪祥,张海越,宋骐. 基于图像识别的AH36钢激光焊缝节点定位技术研究. 计算机测量与控制. 2023(11): 299-305+314 . 百度学术
    4. 付晓艳,白雪,魏博凝. 基于机器学习的物流自动分拣机器人定位精度检测算法. 自动化与仪器仪表. 2022(08): 245-248+252 . 百度学术
    5. 刘翰林,张荣福. PCB焊点及芯线端头识别. 光学仪器. 2022(05): 42-52 . 百度学术

    其他类型引用(3)

图(29)  /  表(1)
计量
  • 文章访问数:  2779
  • HTML全文浏览量:  1728
  • PDF下载量:  137
  • 被引次数: 8
出版历程
  • 收稿日期:  2020-04-22
  • 修回日期:  2020-06-09
  • 网络出版日期:  2020-09-07
  • 刊出日期:  2021-01-14

目录

/

返回文章
返回