冷凝器中污垢沉积率和剥蚀率的在线测量方法
【专利摘要】一种冷凝器中污垢沉积率和剥蚀率的在线测量方法，其步骤为：S1：获得管壁结垢过程的图像序列；S2：识别沉积与剥蚀颗粒；对采集到的图像序列进行分析，识别出污垢颗粒；进行帧间污垢颗粒的匹配，完成对污垢颗粒运动轨迹的实时追踪；S3：建立颗粒的三维几何模型；根据污垢颗粒三维运动轨迹和方向梯度，求得污垢颗粒在管壁的位置和分布，对污垢颗粒及其运动轨迹进行三维重构，建立颗粒的三维几何模型；S4：求得沉积颗粒的体积和与剥蚀颗粒的体积和，进而得到污垢沉积率和剥蚀率。本发明具有实时性好、准确率高、能够直观反映冷凝器管内污垢动态变化过程等优点。
【专利说明】
冷凝器中污垢沉积率和剥蚀率的在线测量方法
技术领域
[0001] 本发明主要设及到冷凝器领域，特指一种应用冷凝器中的污垢沉积率和剥蚀率的 在线测量方法。
【背景技术】
[0002] 目前，冷凝器在电力、化工、冶金、机械等各大领域得到广泛应用。然而，冷凝器中 由于污垢的存在，大大降低了换热面的传热能力，使得汽轮机排汽压力和排汽溫度升高，输 出功率减少，热耗率增加，从而导致机组经济效益下降。此外，污垢的积聚，也常常引起冷凝 管的局部腐蚀乃至穿透，严重威胁到冷凝器的安全运行。正因如此，必须采取有效的污垢控 制措施。
[0003] 冷凝器结垢过程是污垢沉积与污垢剥蚀相互作用的结果，沉积使污垢增加，剥蚀 使污垢减少，沉积率与剥蚀率之差即为污垢结垢速率，污垢控制的本质在于控制沉积率、剥 蚀率，进而达到控制结垢速率的目的。因此，对于冷凝器污垢控制，正确的方法应当是根据 污垢沉积率、剥蚀率的动态变化自适应地调节控制量。然而，由于很难实时、精确地把握管 内污垢的动态变化，且冷凝器污垢是多种颗粒污垢的混合体，从而使得对污垢沉积率、剥蚀 率在线测量的方法的研究难度加大，进展缓慢，至今都没有很好的方法被提出。
【发明内容】
[0004] 本发明要解决的技术问题就在于:针对现有技术存在的技术问题，本发明提供一 种实时性好、准确率高、能够直观反映冷凝器管内污垢动态变化过程的冷凝器中污垢沉积 率和剥蚀率的在线测量方法。
[0005] 为解决上述技术问题，本发明采用W下技术方案：
[0006] -种冷凝器中污垢沉积率和剥蚀率的在线测量方法，其步骤为：
[0007] SI:获得管壁结垢过程的图像序列；
[000引S2:识别沉积与剥蚀颗粒;对采集到的图像序列进行分析，识别出污垢颗粒;进行 帖间污垢颗粒的匹配，完成对污垢颗粒运动轨迹的实时追踪；
[0009] S3:建立颗粒的=维几何模型;根据污垢颗粒=维运动轨迹和方向梯度，求得污垢 颗粒在管壁的位置和分布，对污垢颗粒及其运动轨迹进行=维重构，建立颗粒的=维几何 模型；
[0010] S4:求得沉积颗粒的体积和与剥蚀颗粒的体积和，进而得到污垢沉积率和剥蚀率。
[0011] 作为本发明方法的进一步改进:在所述步骤S4中，根据污垢颗粒的=维几何模型， 求得沉积颗粒的体积和剥蚀颗粒的体积，污垢沉积率、剥蚀率则为单位时间内污垢的沉积 量与剥蚀量，即沉积颗粒的体积和与剥蚀颗粒的体积和。
[0012] 作为本发明方法的进一步改进:所述步骤Sl中是通过图像采集装置实时采集管壁 内的图像，得到管壁结垢过程的图像序列。
[0013] 作为本发明方法的进一步改进:所述步骤S2中，识别颗粒的方法为:先是采用二维 高斯函数描述污垢颗粒的灰度分布，通过对颗粒灰度分布进行函数拟合获得颗粒几何形状 参数，根据颗粒的几何形状参数来识别颗粒。
[0014] 作为本发明方法的进一步改进:所述步骤S2中，颗粒匹配的方法为基于几何特征 相似度和归一化转动惯量特征相似度结合的方法。
[0015] 作为本发明方法的进一步改进:所述步骤S3中，污垢颗粒和污垢运动轨迹的=维 重构方法为:先是基于双目视觉，计算出颗粒的空间位置，实现颗粒的=维重构；之后结合 前面颗粒在各帖间的运动轨迹信息，进行颗粒运动轨迹的=维重构。
[0016] 与现有技术相比，本发明的优点在于：
[0017] 1、本发明本发明的冷凝器中污垢沉积率和剥蚀率的在线测量方法，原理简单、操 作简便、实时性高，能够真正实现对污垢沉积率、剥蚀率的在线测量。
[0018] 2、本发明的冷凝器中污垢沉积率和剥蚀率的在线测量方法，准确度高、能够直观 的反应冷凝器管内污垢的动态变化过程，从而准确得出污垢颗粒在管壁的位置和分布。
[0019] 3、本发明的冷凝器中污垢沉积率和剥蚀率的在线测量方法，易实现，精准可靠，所 需成本低，并能应用于复杂环境当中，有利于管内污垢的控制，提高了运行效率，易推广，增 大了社会效益。
【附图说明】
[0020] 图1是本发明方法的流程示意图。
[0021] 图2是本发明在具体应用实例中进行粒子匹配时的原理示意图。
[0022] 图3是本发明在具体应用实例中对粒子进行=维空间坐标重构的示意图。
[0023] 图4是本发明在具体应用实例中对污垢颗粒进行匹配时帖间捜寻流程的原理示意 图。
[0024] 图5是本发明在具体应用实例中所采用图像采集装置的示意图。
【具体实施方式】
[0025] W下将结合说明书附图和具体实施例对本发明做进一步详细说明。
[0026] 如图1所示，本发明的冷凝器中污垢沉积率和剥蚀率的在线测量方法，步骤为：
[0027] SI:获得管壁结垢过程的图像序列；
[0028] 通过图像采集装置(如:摄像头)实时采集管壁内的图像，得到管壁结垢过程的图 像序列；
[0029] S2:识别沉积与剥蚀颗粒；
[0030] 对采集到的图像序列进行分析，识别出污垢颗粒;进行帖间污垢颗粒的匹配，从而 实现对污垢颗粒运动轨迹的实时追踪。
[0031] S3:建立颗粒的S维几何模型；
[0032] 根据污垢颗粒=维运动轨迹和方向梯度，可求得污垢颗粒在管壁的位置和分布； 良P，对污垢颗粒及其运动轨迹进行=维重构，建立颗粒的=维几何模型；
[0033] S4:求得沉积颗粒的体积和与剥蚀颗粒的体积和，进而得到污垢沉积率和剥蚀率；
[0034] 根据污垢颗粒的=维几何模型，求得沉积颗粒的体积和与剥蚀颗粒的体积，得到 污垢沉积率和剥蚀率。其中，污垢沉积率、剥蚀率是指单位时间内污垢的沉积量与剥蚀量， 即沉积颗粒的体积和与剥蚀颗粒的体积和，根据建立的污垢=维几何模型可求得颗粒体 积。
[0035] 利用双目视觉传输的图像信息计算出颗粒的空间位置，结合前面对颗粒的实时追 踪(二维追踪轨迹），实现颗粒运动轨迹的=维重构。建立颗粒的=维几何模型，求得颗粒体 积，从而实现污垢沉积率、剥蚀率的在线测量。
[0036] 在上述步骤S2中，污垢颗粒运动轨迹实时追踪关键在于颗粒的识别、匹配。本发明 中先是采用二维高斯函数描述污垢颗粒的灰度分布，通过对颗粒灰度分布进行函数拟合获 得颗粒几何形状参数，根据颗粒的几何形状参数来识别颗粒，而实现颗粒匹配的算法是提 出一种基于几何特征相似度和归一化转动惯量特征相似度结合的算法。
[0037] 进一步，本发明针对在颗粒识别、匹配过程中遇到的颗粒重叠、被遮挡的问题，提 出一种帖间边缘预测算法、帖间灰度相关性修正补偿算法。至于存在的颗粒絮凝变形现象， 则采用聚类预估方法。
[0038] 在上述步骤S3中，污垢颗粒和污垢运动轨迹的=维重构，先是基于双目视觉，计算 出颗粒的空间位置，实现颗粒的=维重构。之后结合前面颗粒在各帖间的运动轨迹信息，进 行颗粒运动轨迹的=维持=维重构。
[0039] 在具体应用实例中，本发明的冷凝器中污垢沉积率和剥蚀率的在线测量方法，详 细步骤为：
[0040] SlOO:获得管壁结垢过程的图像序列；
[0041] 如，通过摄像头采集图像，从而获得管壁结垢过程的图像序列。
[0042] S200:利用图像序列，实现颗粒运动轨迹的S维重构；
[0043] 目P :对图像序列进行处理，提取出每帖图像中的污垢颗粒，并对每一帖用双目摄像 机得到污垢颗粒图像进行=维重构W获得其空间位置。
[0044] S300:对图像序列中污垢点进行实时追踪，从而获得图像序列中污垢颗粒的运动 路径。
[0045] S400:定义污垢沉积率与剥蚀率为单位时间内污垢的沉积量和剥蚀量。
[0046] 参见图2~图4,在具体应用时，步骤S200的具体流程如下：
[0047] S201:先对污垢颗粒的位置进行=维重构，建立污垢颗粒空间位置所依赖的世界 坐标系，其同左相机坐标系的旋转矩阵和平移矩阵分别为化和化，其同右相机坐标系的旋转 矩阵和平移矩阵分别为化和Tr,运些矩阵都是固定数值矩阵，由标定可得。
[004引 S202:设左右相机中的粒子像P '和P是一对匹配粒子，其粒子=维空间坐标=维重 构几何原理如图3所示，其中I和I'分别为左右相机的相机中屯、，粒子P所在的=维空间位置 即为两直线PI和p'r的交点。
[0049] 设空间粒子G在所定义的世界坐标系中的空间坐标为扣节，是），而其在左右相 机中的粒子像分别为P和P'，且它们在世界坐标系中的坐标分别设为和 (巧，皆，巧）。因易求得粒子像巧日P'分别在左右相机坐标系的坐标，分别设为(Xp，Yp，Zp) 和(Xp'，押'，Zp'），
[0050] 故可得到如下计算式： 陶]
a)
[00对
樹
[0化3] 又易知左右相机中屯、I和I '在左右相机坐标系的坐标分别为（0，0，Zi )和（0，0， Zr),将其转换为上述的世界坐标，分别设它们在世界坐标系的坐标为〇￥7，叹.马）和 (乃，F/-，巧)，则坐标转换公式为：
[0化4] (3)
[0055] (4)
[0056] 由上述计算式已经将点全部转换到定义的世界坐标系中，易知，直线PI和GI是同 一直线，在世界坐标系中的表达式相同。同样，p'r和GI'也是同一直线，拥有相同的表达 式。故空间粒子G的坐标可由如下式(5)至(9)计算式所得：
[0化7]
[0化引 [0化9]
[0060]
[0061] 觸
[0062] 役）
[0063] 在具体应用时，步骤S300的具体流程如下：
[0064] 由图4所示，对上述步骤S2中获得的图像序列中污垢颗粒位置及提取出的污垢颗 粒图像进行处理，对图像序列中第一帖图像检测污垢颗粒及其位置点，在下一帖中尝试跟 踪运些颗粒，即找到运些污垢颗粒在新的图像帖中的位置，如此跟踪完整的图像序列，如图 4所示，图像Ii和12是相邻帖。同时，对无匹配的污垢颗粒依据其所处位置及运动趋势进行标 记，区分为沉积粒子和剥蚀粒子。
[0065] 根据上述操作完成后，根据实际需要，可W选择将颗粒区分为沉积颗粒、剥蚀颗 粒、逃逸颗粒。其中，逃逸颗粒指所有帖中都有匹配颗粒且无沉积趋势。由于污垢颗粒点会 随水流移动，并且发生旋转等现象。因此，为了在帖间匹配污垢粒子，本发明采取如下技术 方案：
[0066] 首先，根据工况环境确定在下一帖中的查询区域，应用初始帖中颗粒特征参数进 行相似匹配，其中特征参数包括颜色、几何形状、图像质量、图像转动惯量、正投影面积W及 角点的相对位置。当完成由图4所示的流程进行颗粒特征匹配部分时，需要根据新一帖的匹 配颗粒对颗粒的特征参数进行更新，已获得污垢颗粒的完整特征。
[0067] 其中，图像质量及图像转动惯量方案如下：
[0068] 将双摄像头得到的两幅图像进行灰度拉伸处理W提取污垢颗粒的二值化图像。如 图5所示，在具体应用实例中，图像采集装置包括光源1、计算机3和双目相机4,其中光源1具 有如图一般的横截面2,上述图像采集装置用来对冷凝管5进行实时监测。其中，像素的灰度 值应与距离光源的远近负相关，即距离光源越远，灰度值越低。设提取的污垢颗粒的图像像
[0070] (10) 素区域范围为l《x《M，l《y《N，令f(x，y)表示对应像素灰度值。[0069] 定义灰度平均值为图像质量m:
[0071] 痒方式如下：
[0072] (Ii)
[007；3] (城
[0074]
[0075] (13)
[0076] 在整个追踪污垢颗粒运动路径中，针对颗粒间的重叠、被遮挡现象，拟用帖间灰度 相关性修正补偿算法，即通过颗粒运动的序列图像，对前一帖图像的像素进行移位和插值 来补偿当前帖的边缘缺损。而对于颗粒在运动中出现絮凝变形现象，本案例采用聚类预估 的方法，W追踪的粒子为聚类中屯、，析晶颗粒为聚类数据，选取合适的距离度量，采用k均值 聚类算法，获得颗粒子群，W此颗粒子群作为对被追踪颗粒絮团形状的预估，得到颗粒子群 的形状参数，并与实际获得的絮团图像进行几何相似度特征比较，从而解决絮凝造成的颗 粒匹配问题。
[0077] 在具体应用时，步骤S400的具体流程如下：
[0078] 通过上述的步骤，摄像头记录了管壁从光滑到结垢的全过程，An和Bn分别表示第n 帖图像管壁区域的沉积颗粒和剥蚀颗粒，在前面的粒子追踪中对运两种趋势的粒子分别作 了两个不同的表识。经过=维重构后，对沉积颗粒和剥蚀颗粒作不同的标记(如用红色标记 沉积颗粒，蓝色标记剥蚀颗粒），加 W区分。设VAn和VBn分别表示第n帖图像=维重构后管壁 区域沉积颗粒和剥蚀颗粒的体积，利用步骤S3中得到的污垢粒子的特征参数，利用其角点 相对位置重构其几何模型，进而计算得到其体积。
[0079] 又令T表示摄像机采集两帖图像的相差时间，Van和Vbn表示帖与帖间的沉积率和剥 蚀率。则右々n下管井.
[0080] (14)
[0081] (15)
[0082] 计时，经过时间t后的沉积率Va和剥蚀率Vb有：
[00削 狂目）
[0084] (切
[0085] W上仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例， 凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的 普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，应视为本发明的保护 范围。
【主权项】
1. 一种冷凝器中污垢沉积率和剥蚀率的在线测量方法，其特征在于，步骤为： Sl:获得管壁结垢过程的图像序列； S2:识别沉积与剥蚀颗粒;对采集到的图像序列进行分析，识别出污垢颗粒;进行帧间 污垢颗粒的匹配，完成对污垢颗粒运动轨迹的实时追踪； S3:建立颗粒的三维几何模型；根据污垢颗粒三维运动轨迹和方向梯度，求得污垢颗粒 在管壁的位置和分布，对污垢颗粒及其运动轨迹进行三维重构，建立颗粒的三维几何模型； S4:求得沉积颗粒的体积和与剥蚀颗粒的体积和，进而得到污垢沉积率和剥蚀率。2. 根据权利要求1所述的冷凝器中污垢沉积率和剥蚀率的在线测量方法，其特征在于， 在所述步骤S4中，根据污垢颗粒的三维几何模型，求得沉积颗粒的体积和剥蚀颗粒的体积， 污垢沉积率、剥蚀率则为单位时间内污垢的沉积量与剥蚀量，即沉积颗粒的体积和与剥蚀 颗粒的体积和。3. 根据权利要求1所述的冷凝器中污垢沉积率和剥蚀率的在线测量方法，其特征在于， 所述步骤Sl中是通过图像采集装置实时采集管壁内的图像，得到管壁结垢过程的图像序 列。4. 根据权利要求1或2或3所述的冷凝器中污垢沉积率和剥蚀率的在线测量方法，其特 征在于，所述步骤S2中，识别颗粒的方法为:先是采用二维高斯函数描述污垢颗粒的灰度分 布，通过对颗粒灰度分布进行函数拟合获得颗粒几何形状参数，根据颗粒的几何形状参数 来识别颗粒。5. 根据权利要求1或2或3所述的冷凝器中污垢沉积率和剥蚀率的在线测量方法，其特 征在于，所述步骤S2中，颗粒匹配的方法为基于几何特征相似度和归一化转动惯量特征相 似度结合的方法。6. 根据权利要求1或2或3所述的冷凝器中污垢沉积率和剥蚀率的在线测量方法，其特 征在于，所述步骤S3中，污垢颗粒和污垢运动轨迹的三维重构方法为:先是基于双目视觉， 计算出颗粒的空间位置，实现颗粒的三维重构;之后结合前面颗粒在各帧间的运动轨迹信 息，进行颗粒运动轨迹的三维重构。
【文档编号】G01N15/04GK105954163SQ201610244912
【公开日】2016年9月21日
【申请日】2016年4月19日
【发明人】樊绍胜, 李寒松, 叶紫权, 甘彬
【申请人】长沙理工大学