• 中国科技论文统计源期刊
  • 中文核心期刊
  • 中国科技核心期刊
  • 中国机械工程学会无损检测分会会刊
高级检索

锅炉膜式水冷壁漏磁检测磁化结构设计及分析

陈涛, 孙晗, 杨象岳, 刘延雷, 金斌戈, 熊伟东

陈涛, 孙晗, 杨象岳, 刘延雷, 金斌戈, 熊伟东. 锅炉膜式水冷壁漏磁检测磁化结构设计及分析[J]. 无损检测, 2025, 47(3): 1-5. DOI: 10.11973/wsjc240306
引用本文: 陈涛, 孙晗, 杨象岳, 刘延雷, 金斌戈, 熊伟东. 锅炉膜式水冷壁漏磁检测磁化结构设计及分析[J]. 无损检测, 2025, 47(3): 1-5. DOI: 10.11973/wsjc240306
CHEN Tao, SUN Han, YANG Xiangyue, LIU Yanlei, JIN Binge, XIONG Weidong. Design and analysis of magnetization structure for leakage magnetic detection of boiler membrane waterwall[J]. Nondestructive Testing, 2025, 47(3): 1-5. DOI: 10.11973/wsjc240306
Citation: CHEN Tao, SUN Han, YANG Xiangyue, LIU Yanlei, JIN Binge, XIONG Weidong. Design and analysis of magnetization structure for leakage magnetic detection of boiler membrane waterwall[J]. Nondestructive Testing, 2025, 47(3): 1-5. DOI: 10.11973/wsjc240306

锅炉膜式水冷壁漏磁检测磁化结构设计及分析

基金项目: 

国家市场监督管理总局科技计划项目 2022MK162

详细信息
    作者简介:

    陈涛(1986—),男,博士,高级工程师,主要从事特种设备安全评价及无损检测技术研究工作

    通讯作者:

    熊伟东,chentao_zju@163.com

  • 中图分类号: TH49;TG115.28

Design and analysis of magnetization structure for leakage magnetic detection of boiler membrane waterwall

  • 摘要:

    在采用漏磁检测技术对锅炉水冷壁管进行缺陷检测时,受水冷壁鳍片结构的影响,水冷壁管缺陷处的漏磁信号强度会被削弱,从而导致检测灵敏度下降。针对此问题,提出了增加磁化单元,对待检水冷壁管两侧管道同时进行磁化的改进方案,并且对该方案的可行性进行了深入研究。建立了并列三个磁化单元的有限元模型,对比分析了增加磁化单元前后水冷壁管缺陷处的漏磁信号强度,最后通过试验验证了数值计算结果的准确性和改进方案的可行性,为锅炉水冷壁管漏磁检测设备的开发和应用提供参考。

    Abstract:

    When using magnetic flux leakage detection technology to detect the defect of boiler water wall tube, due to the influence of the fin structure of the water-cooled wall, the magnetic flux leakage signal strength at the defect location of the water-cooled wall tube will be weakened, leading to a decrease in detection sensitivity. In response to the above issues, the author of this article proposed an improved scheme of adding a magnetization unit to simultaneously magnetize the pipes on both sides of the water-cooled wall pipe to be inspected. Furthermore, the author conducted in-depth research on the feasibility of this scheme. A finite element model of the magnetization structure of the boiler water-cooled wall tube containing three magnetization units was established. The magnetic leakage signal strength at the defects of the water-cooled wall tube before and after adding magnetization units was compared and analyzed. Finally, the accuracy of the numerical calculation results and the feasibility of the improvement plan were verified through experiments. The research content of this article could provide reference for the development and application of magnetic leakage detection equipment for boiler water-cooled wall tubes.

  • 随着我国电力工业的高速发展,绝缘子的使用量大幅增加。硅橡胶复合绝缘子因耐污闪性能高、比强度高、性价比高等特点,较其他类型绝缘子的使用更为广泛[1]。然而,复合绝缘子由于材料、结构以及制作工艺的特殊性,在生产过程中容易产生端部密封不良、芯棒-护套黏接不良等缺陷;在使用过程中,受到酸雨等自然因素的影响,其内部缺陷会进一步扩大,最终导致断串及芯棒击穿故障[2-4],对输电线路的安全稳定运行构成严重威胁。因此,有必要在生产端对复合绝缘子的黏接质量进行检测评估,防患于未然。

    目前,相关学者针对复合绝缘子内部缺陷的检测提出了诸多方法,如太赫兹无损检测法[5]、微波致热检测法[6]、超声导波检测法[7]、错位散斑干涉检测法[8]等。但上述检测方法普遍检测效率低且实际应用的难度较大,而相控阵超声检测技术具有安全性高、灵敏度高、检测速度快、成本较低等突出优点,较为适用于廉价且用量大的复合绝缘子检测。谢从珍等[9]利用相控阵超声波检测仪对4种复合绝缘子典型内部缺陷进行了定性识别。陈海燕等[10]提出了一种基于水囊柔性耦合的相控阵检测方法,该方法适用于带电现场检测。梁进祥等[11]采用了仿真模拟与试验相结合的方式对复合绝缘子的内部破损和贯穿缺陷进行了定性分析。然而,上述研究均未对脱黏缺陷信号进行深入研究。

    在实际检测中发现,复合绝缘子的直径较小,相控阵换能器难以精确对准绝缘子中心,且声能会被曲面部分反射,同时硅橡胶材料声能的高衰减性会导致回波信号幅度大幅下降[12]。因此,为了提升分层界面的回波幅值,需将采集增益调到最大才能分辨出有效信号,但也产生了较大的固有电路噪声和结构噪声,严重影响信号分析及成像。另外,由于复合绝缘子的黏接结构特点,一旦涂胶过量,硅橡胶/玻璃纤维双层结构会转变为硅橡胶/黏接胶层双层结构,从而导致C扫描成像中出现异常高亮区域,严重影响对缺陷的判断[13]。故亟需一种信号处理方法消除这些影响。

    曹弘毅等[14]利用相控阵超声技术对复合材料层压板的分层缺陷进行了C扫描成像和定量评估,但未提出C扫描图像的优化方法。常青等[15]为了提高对钛合金扩散焊接界面中微小缺陷的识别,提出基于小波变换的C扫描图像融合算法,提高了图像的对比度,但不适用于黏接界面异常信号的消除。吴杰等[16]为提升玻璃纤维复合材料的C扫描成像效果,提出一种基于VMD的相关性分析成像算法,提升了成像质量,但是VMD参数需要人工试验,且其相关性分析只用到了信号合成,并不能解决异常高亮图像的问题。

    综上所述,文章采用相控阵超声检测方法对复合绝缘子的脱黏缺陷进行试验,针对C扫描图像中存在异常高亮区域、对比度低的问题,提出一种基于SSA-VMD的DTW相关性分析方法。

    由于异质材料的声学性质不同,超声波P0从一种物质透射进入另一种物质时,除了产生透射声波Pt外,还会在异质界面处返回一种界面回波Pr,其简化模型如图1所示。

    图  1  超声传播简化模型

    声压反射率r和声压透射率T与介质阻抗的关系[17]

    {r=Z2-Z1Z2+Z1T=2Z2Z2+Z1 (1)

    式中:Z1为介质1的阻抗;Z2为介质2的阻抗。

    查阅资料可知,常温下空气的声阻抗为0.000 41 g/(cm2·s);硅橡胶的声阻抗为1.4 g/(cm2·s);玻璃纤维的声阻抗为6.04 g/(cm2·s)。根据式(1)可计算出从橡胶到空气的声压反射率约为−1,从橡胶到玻璃纤维的声压反射率约为0.62,所以脱黏处界面回波的相位应与黏接良好处界面回波的相位相反并且幅值更大。

    变分模态分解(Variational mode decomposition,VMD)是在2014年由DRAGOMIRETSKIY等[18]提出的一种信号处理算法。

    首先构造约束变分模型,经典约束变分模型为

    min{uk},{ωk}{kt[(δ(t)+jπt)uk(t)]e-jωkt22} (2)
    s.t.   k=1Kuk=f (3)

    式中:f为原始输入信号;K为信号分解的模态个数;{uk}={u1,,uk},{ωk}={ω1,,ωk}分别为各模态分量及其中心频率;t为时间;δ(t)为狄拉克函数;∂tuk(t)的单边频谱。

    引入拉格朗日乘子λ和二次惩罚项α,原始模型转换为非约束问题,得到增广拉格朗日表达式

    L({uk},{ωk},λ)=αkt[(δ(t)+jπt)uk(t)]e-jωkt22+f(t)-kuk(t)22+λ(t),f(t)-kuk(t) (4)

    通过交替方向乘子算法(ADMM)求解式的鞍点,即可求得各模态及其中心频率,结果为

    λ^n+1(ω)=λ^n(ω)+γ(f^(ω)-ku^kn+1(ω)) (5)
    u^kn+1(ω)=f^(ω)-iku^i(ω)+λ^(ω)21+2α(ω-ωk)2 (6)
    ωkn+1=0ω|u^k(ω)|2dω0|u^k(ω)|2dω (7)

    式中:n为迭代收敛次数;u^k(ω),f^(ω),λ^(ω)分别为u^k(t),f^(t),λ^(t)的傅里叶变换;γ为噪声容忍度。

    将式(6)结果进行傅里叶反变换,即可求得各IMF分量。

    麻雀搜索算法(Sparrow search algorithm,SSA)是在2020年由薛建凯等[19]提出的一种元启发式算法。该算法模仿自然界中麻雀的群体觅食和反觅食行为,将觅食群体划分为探索者和追随者,探索者为能够寻找较好食物的麻雀,其余麻雀为追随者。受到探索者影响,反觅食群体会进行侦察并根据觅食者动作做出相应行为。其算法模型大致如下。

    群体位置初始化

    X=lb+rand*(ub-lb) (8)

    式中:ublb分别为种群上下位置边界;rand为(0,1)之间的随机数。

    d维空间的探索者位置更新表达式为

    Xi,jm+1={Xi,jm · exp(-iαMmax)R2<STXi,jm+Q · LR2ST (9)

    式中:m为当前迭代次数;Mmax为最大迭代次数;Xi,j为当前迭代中第i只麻雀在第j维的位置信息;α为(0,1]中的随机数;STR2分别为安全值和预警值参数;Q为均值为0、方差为1的随机数;L为元素均为1的1×d矩阵。

    追随者的位置更新表达式为

    Xi,jm+1={Qexp(Xworst,jm-Xi,jmi2)i>Npop/2Xp,jm+1+|Xi,jm-Xp,jm+1| · A+ · LiNpop/2 (10)

    式中:Xp,j为发现者的最佳位置;Xworst,j为当前迭代中的全局最差麻雀位置;Npop为麻雀种群规模;A为1×d的矩阵,且A+=AT(AAT)-1

    侦察者的位置更新表达式为

    Xi,jm+1={Xbest,jm+β · |Xi,jm-Xbest,jm|fifgXi,jm+K · (Xi,jm-Xworst jm(fi-fw)+ε)fi=fg (11)

    式中:Xbest,j为当前迭代中的最佳麻雀位置;β为步长控制参数;K′为[-1,1]中的随机数;ε为极小的常数;fi为预警者的适应度;fgfw分别为当前麻雀个体的最优和最差适应度。

    VMD算法的困难性在于Kα的确定,不合适的参数会导致出现模态混叠问题。由此引入SSA算法进行VMD分解参数的寻优,同时采用排列熵作为SSA的适应度函数。设有一长度为L的时间序列X={x(1),x(2),x(3), ,x(L)},其排列熵计算公式为

    H(m)=-j=1KPjlnPj (12)

    式中;Pj为第j种排列模式在整个时间序列中出现的概率。

    SSA-VMD算法具体的寻优步骤为:①设置SSA算法参数,包括寻优的上下限值,种群规模和最大迭代数;②调用排列熵作为适应度函数,得到初始化后的排列熵和对应的最佳参数;③以初始化的参数组合为基础,更新麻雀的位置和速度;④在不同麻雀的位置下对信号进行VMD分解,计算当前位置下的排列熵;⑤对比各个排列熵大小,确定当下最小排列熵,并更新个体极值和全局极值;⑥更新麻雀的位置和速度,重复步骤④,⑤。循环迭代直到迭代数到达预定最大值,跳出循环得到最小排列熵和最佳参数组合。

    动态时间规整(Dynamic time warping,DTW)是一种衡量两个时间序列之间相似度的方法,主要应用于比较语音信号的相似性。其最大的优点是能够摆脱序列长度的影响,主要比较两段序列的趋势相似性。文章引入DTW距离作为脱黏信号和非脱黏信号的相似性判断方法,设要计算相似度的两个时间序列为AB,采用递归算法计算最短路径长度。两个序列各点之间的距离矩阵M(i,j)为

    M(i,j)=|A(i)-B(j)| (13)

    式中:A(i)和B(j)分别为时间序列对应点的值。

    DTW距离矩阵为

    Lmin(1, 1)=M(1, 1) (14)
    L(i,j)=min{Lmin(i,j-1),Lmin(i-1,j),Lmin(i-1,j-1)}+M(i,j) (15)

    式中:L(i,j)为两个时间序列的规整路径距离,即DTW距离。

    通过对比不同信号和参考信号的DTW距离,设定合适的阈值,即可实现对不同信号的相似度评估。文章算法的总体流程图如图2所示。

    图  2  文章算法的整体流程图

    为确定合适的聚焦法则,使用CIVA软件构建水/硅橡胶双层介质模型,模拟使用频率为10.0,5.0,2.5 MHz的相控阵换能器(阵元中心距为0.3 mm)。在32阵元,16阵元,8阵元的聚焦条件下,分别距离硅橡胶表面10,20,30,40,50 mm处进行声场模拟,虚拟聚焦点均设置在硅橡胶5 mm深处。得到换能器参数、聚焦阵元个数与预设聚焦点处横向分辨率的关系曲线如图3所示。文章试验要求精度为1 mm,为了兼顾计算机性能,故,采用频率为10 MHz的相控阵换能器进行8阵元聚焦的线性扫查较为合适。

    图  3  换能器参数与预设点处横向分辨率的关系曲线

    试验对象为不带伞裙的复合绝缘子试件(直径为34 mm),其硅橡胶护套厚度为5 mm,在其中埋入一个直径约为10 mm的圆形缺陷。整体的试验装置如图4所示,该装置主要由上位机控制系统、机械臂运动控制装置和信号采集装置组成。试验使用的信号采集设备为M2M公司生产的Multi 2000型便携式相控阵超声采集卡,使用的相控阵换能器型号为10L64-0.3×10,试验主机的CPU型号为i9-13900HX,GPU型号为RTX4070,内存大小为32 G。共进行了3组试验,分别在距绝缘子表面10,30,50 mm处采集信号。

    图  4  整体试验装置

    试验所获得的典型信号如图5所示,可见,脱黏信号之间的波形都极为相似,且与界面信号波形差距较大;异常高亮信号之间的波形也极为相似,且其波形几乎为界面信号的放大版。

    图  5  距绝缘子表面10,30,50 mm处的典型信号

    图5中原始信号的噪声较大,很可能会对信号的互相关性判断造成影响,因此,采用SSA-VMD算法对10 mm处的界面信号进行模态分解,得到的VMD分解图和SSA算法适应度曲线分别如图6,7所示。

    图  6  界面信号各分解模态及其频谱
    图  7  SSA算法适应函数收敛曲线

    图6可知,原始信号被分解为两个模态,IMF1模态为噪声的主要存在模态,频率集中在8 MHz,IMF2模态为有效信号的主要存在模态,频率集中在2 MHz。

    图7可知,算法在12次迭代时适应度达到最小,说明当K=2,α=6 892时,分解效果最好。为验证参数普适性,对脱黏信号和异常高亮信号进行了测试,得到模态分解图如图8,9所示。

    图  8  脱黏信号各分解模态及其频谱
    图  9  异常高亮信号各分解模态及其频谱

    图8,9可知,脱黏信号被分解为IMF1和IMF2,IMF2为所需模态;异常高亮信号被分解为IMF1和IMF2,IMF2为所需模态。综上所述,K=2,α=6 892为较佳的VMD参数。

    将脱黏信号、界面信号和异常高亮信号的有效回波进行归一化后的对比如图10所示,可以看到界面信号和异常高亮信号的形状走势极为相似,而脱黏信号与界面信号波形差距较大。

    图  10  3种典型信号的IMF2模态对比

    为有效判断3种信号的互相关程度,分别在异常高亮区域和脱黏区域采集5组信号,与正常界面信号的DTW距离如表1所示。将阈值设定为5,大于5即判定为脱黏区域,相关系数矩阵计算公式为

    CM=sgn(DTW-5) (16)

    式中:CM为相关系数矩阵;sgn(·)为符号函数;DTW为DTW距离。

    Table  1.  5组信号的DTW距离
    编号异常高亮区域信号与正常界面信号的DTW距离脱黏区域信号与正常界面信号的DTW距离
    11.514 48.812 6
    21.741 88.236 4
    31.652 48.127 8
    41.977 78.897 7
    51.687 37.987 8
    下载: 导出CSV 
    | 显示表格

    对整幅图像的A扫描信号进行相关性分析,得到相关系数矩阵,再与原始C扫描图进行加权成像,得到优化后的C扫描图像,将其二值化后发现,原始图像的异常高亮区域已经被完全消除。优化前后及二值化后的C扫描图如图11所示。

    图  11  优化前后及二值化后的C扫描图对比

    文章提出一种基于SSA-VMD信号预处理和DTW相关性分析的C扫描成像方法。试验结果表明,该方法能够有效去除C扫描图像中的异常高亮信号,提高图像对比度,为后续的自动化缺陷判别奠定基础,同时也为相似构成工件的脱黏缺陷检测提供一种成像思路。

  • 图  1   水冷壁磁化结构示意

    图  2   有无鳍片水冷壁的漏磁场磁通量密度对比

    图  3   三磁化单元并列布置的磁化结构模型及有限元网格

    图  4   磁化结构模型的三维体磁通密度模

    图  5   缺陷二维漏磁场分布云图

    图  6   漏磁场信号提取路径示意

    图  7   3种磁化结构下,水冷壁管缺陷处漏磁场磁通量密度的径向分量和轴向分量

    图  8   试验装置布置现场

    图  9   3种不同检测方案的霍尔电压数据对比

    图  10   3种不同检测方案得到的缺陷特征参数对比

    Table  1   膜式水冷壁磁化结构模型参数

    参数数值
    水冷壁管外径60
    水冷壁管壁厚4
    水冷壁鳍片厚度5
    水冷壁管中心距95
    磁铁尺寸(长×宽×高,下同)90×60×25
    衔铁尺寸90×200×25
    磁化器提离值4
    传感器提离值1
    缺陷尺寸(长×宽×深)3×3×3
    下载: 导出CSV

    Table  2   3种不同检测方案的漏磁信号波峰波谷差值

    缺陷尺寸/mm无鳍片电压/V有鳍片电压/V有鳍片三磁化头电压/V
    2×2×30.1860.1080.180
    3×3×30.2680.1510.257
    4×4×30.3280.2000.320
    5×5×30.5090.3990.496
    下载: 导出CSV
  • [1] 杜佼伟,李霄,吕品正,等. 660 MW锅炉水冷壁管道泄漏分析[J]. 金属加工:热加工,2024(6):45-49.
    [2] 郭海霞,张欣耀,袁晓冬. 锅炉水冷壁管失效分析[J]. 热处理,2024,39(2):55-57,60.
    [3] 赵文宇. 电站锅炉水冷壁管的泄漏原因及应对措施研究[J]. 中国设备工程,2023(1):134-136.
    [4] 黄斌鑫. 热电厂锅炉水冷壁管失效原因探究[J]. 全面腐蚀控制,2022,36(9):36-38,46.
    [5] 周公文. 600 MW超临界锅炉水冷壁管爆管原因分析[J]. 化工装备技术,2022,43(2):26-28.
    [6] 周爽. 浅谈某电厂2号锅炉水冷壁管泄漏原因[J]. 锅炉制造,2022(4):19-20,23.
    [7] 王红轩,李芳. 一种煤粉锅炉水冷壁管爆管原因分析[J]. 中国科技信息,2022(12):66-68.
    [8] 吕运容,段志宏,连玮琦,等. 循环流化床锅炉水冷壁管泄漏原因分析[J]. 热加工工艺,2023,52(17):139-142.
    [9] 陈申乾,黄松岭,赵伟,等. 燃煤锅炉热力管道缺陷检测技术[J]. 无损检测,2017,39(1):42-48.
    [10] 张重阳. 漏磁技术在蒸汽锅炉水冷壁检测中的应用[J]. 现代电子技术,2010,33(19):139-140,143.
    [11] 陈世旺. 在役锅炉水冷壁管腐蚀检测技术研究[J]. 质量技术监督研究,2012(6):17-19.
    [12] 刘安仓,叶鑫,孙伟生,等. 火电厂锅炉水冷壁管缺陷检测技术研究进展[J]. 工业安全与环保,2023,49(8):64-69.
    [13] 陈涛,钮建良,廖晓玲,等. 锅炉水冷壁管鳍片对漏磁检测影响的分析[J]. 无损检测,2024,46(1):1-5.
图(10)  /  表(2)
计量
  • 文章访问数:  19
  • HTML全文浏览量:  5
  • PDF下载量:  17
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-07-09
  • 刊出日期:  2025-03-09

目录

/

返回文章
返回