技术领域
[0001] 本
发明属于
多相流测试领域，特别涉及一种管道内周向液膜特征的光学畸变校正和测量方法。
背景技术
[0002] 对于管道内
气液两相流，当气相占据管道截面中心，而液相附着于管道内壁流动时，称为气液环状流。环状流具有较宽的气、液相流速范围，而且液膜流动消耗少，热质交换能
力强，广泛存在于石油、化工、
水利、
能源、
冶金、核工业等国民经济重要的工业生产过程、系统和装备中。对于管道内液膜重要的特征参数，如液膜厚度及其分布、
波速、液膜与管壁间剪切
应力等的准确测量具有非常重要的意义。在高流速下，由于液雾以及气液相间的相互作用，液膜流动呈现非轴对称性，除了沿液膜流动方向的管道轴向液膜流动结构，解析管道圆周截面上的液膜流动特征和流动结构具有更为重要的研究意义。
[0003] 典型的管道内液膜流场测试方法，如基于电导探针的测量方法仅能探测其附近的
信号，且具有平均效应，在液膜厚度较大时，信号灵敏度会降低，即使采用多电导
电极阵列，对于随
时空高速变化的动态液膜而言，仍然难以获得足够高的时空
分辨率。平面激光诱导
荧光(PLIF)技术是一种具有原理优势的流场诊断方法，具有非侵入、高时空分辨率、场测量以及
可视化的突出优势，已成为环状流和液膜研究的前沿和热点。该方法是采用接近荧光物质最大吸收
波长的激光作为激发
光源，照射空间流场某一截面，诱导溶于液相中的荧光分子发出荧光，利用高速摄像机获取荧光图像并进行诊断的一种先进光学诊断方法。
[0004] 针对基于PLIF的液膜流动测试，由于管道以及气液相折射率的差异，以及管道弯曲形状所带来的影响，摄取到的液膜图像难免失真，采集计算得到的实际液膜厚度大于原始液膜厚度。目前为了克服管道内液膜光路畸变和失真带来的影响，使其易于诊断，一般会建立矩形管道，或通过折射率匹配的方法，即在圆管外安装方形管道，其间充满折射率相近的工作液体来进行光路补偿，这种基于外加结构的折射率匹配方法增加了测量的复杂性，在有些场合也难以实现。因此，如何实现管道内液膜特征的普适性、高
精度测量成为具有突破性的关键技术之一。
发明内容
[0005] 为克服
现有技术的不足，本发明旨在提出一种基于靶标校准的管道内周向液膜光学畸变校正和测量方法，不依赖于管道的材料、形状、尺寸以及工作液体的折射率,即使在圆管外不加方形水套也可以实现周向液膜特征的非侵入、高精度测量与分析。本发明采用的技术方案是，管道内周向液膜特征光学畸变校正与测量方法，步骤如下：
[0006] 1)、设计适用于管道内周向液膜测量的圆形校准靶标，采集液相工质条件下待测周向截面的靶标畸变图像，提取靶标畸变图像上的特征点；
[0007] 2)、将靶标畸变图像与理想靶标图像上的特征点进行匹配，基于多项式函数建立原始图像与畸变图像特征点之间的映射数学模型；
[0008] 3)、调整片状激光，使其水平照亮靶标平面处的管道横截面，激发溶解有特定荧光
溶剂的轴向液膜产生荧光，利用高速摄像机摄取圆管内周向液膜的双视
角荧光图像；
[0009] 4)、将采集的液膜原始图像进行灰度化、差影及滤波预处理，并基于映射数学模型对预处理后的液膜图像进行畸变校正，具体分为两步，一是对原畸变图像与畸变校正后图像进行映射变换，二是基于灰度插值确定畸变校正后图像的灰度值；
[0010] 5)、将畸变校正后的双视角图像中心
位置重合，并对摄取的左、右半周液膜图像有效部分进行融合，准确还原周向液膜的流动状态，以便提取周向液膜特征，并进一步分析管道内液膜的演化特性及
传热传质特性。
[0011] 靶标设计为圆形，直径为管道内径2倍，靶标特征由距离靶心不同远近的圆点组成，位于靶心的中心圆点直径为1mm，距离其最近且与竖直方向夹角分别为0°、90°、180°及270°处分别排列4个小圆点，直径为0.25mm，将管道中心作为投影基准，其余圆点按横向与纵向间距1.5mm依次排列，直径为0.5mm，，根据双视角周向液膜的分布特性，采用单台高速摄像机和折光分光光路虚拟两台高速摄像机，将靶标平面置于片状激光照亮的待测圆周截面位置上，在管道内加注液相；
[0012] 将靶标的圆点质心作为待提取特征点，根据光路分析，圆管内靶标中心区域的光路畸变最小，因此可以将中心区域处4个小圆点的位置关系作为已知条件，对靶标畸变图像进行灰度化、差影、滤波以及
阈值分割，提取各圆点特征的质心及面积，根据面积属性获取靶标中心
位姿关系，并利用已知4个小圆点与中心点之间的位置关系，由近及远依次推导出其余圆点的质心坐标，建立理想靶标图像，并得到其特征点坐标及其位姿关系，根据畸变靶标图像与理想图像特征点之间的畸变投影关系匹配特征点，基于多项式函数建模的方法得到畸变校正的数学模型，其原理如下：
[0013] 设p(u,v)为理想图像，p(x,y)为对应畸变的图像，两者的特征点坐标存在非线性对应关系F，即
[0014] (x,y)＝F(u,v) (1)
[0015] 上式可用多项式模型来近似，即
[0016]
[0017] 式中，aij和bij为待定系数，n为多项式次数。通过特征点匹配获得多组匹
配对，基于最小二乘拟合得到式(2)中多项式系数。
[0018] 进一步地，基于靶标校正后，利用激光照射该横截面流场，采用高速摄像机摄取并对畸变的液膜荧光图像进行校正处理，即对图像上每一坐标值(u,v)进行变换F，得到在畸变图像上的对应坐标(x,y)，在获得畸变图像与校正图像坐标之间的变换关系后，采用灰度插值法，确定校正后图像灰度值。
[0019] 本发明的特点及有益效果是：
[0020] 本发明方法通过基于靶标的畸变校正方法，结合平面激光诱导荧光技术，实现了圆管内周向液膜的光学畸变校正与高精度可视化测量。该方法突破了传统测量手段必须建立矩形管道，或通过折射率匹配(在圆管外安装方形透明管，其间充满折射率相近的工作液体)进行光路补偿的局限性，基于畸变校正数学模型实现了管道内周向液膜的高精度、普适性测量，为后续气液环状流中液膜流动结构及演化特性分析奠定了坚实
基础。
附图说明：
[0021] 图1为管道内周向液膜流动特征不同测量方法示意图(以竖直管道为例)。图中，(a)为传统折射率匹配方法示意图，(b)为本发明测量方法示意图。
[0022] 图2为设计的校准靶标与高速摄像机采集的靶标图像。图中，(a)为设计的靶标样式及局部放大图，(b)为实际采集的靶标图像。
[0023] 图3为管道内周向液膜畸变校正前与畸变校正后的图像。图中，(a)为畸变校正前的周向液膜图像(预处理后)，(b)为畸变校正后的液膜图像，(c)为左右视角融合后的液膜图像。
具体实施方式
[0024] 本方法适用于环状流圆管内周向液膜特征的光学畸变校正和测量，同样也适用于其他流型或降膜圆管内周向液膜特征的光学畸变校正和测量。
[0025] 针对圆管内周向液膜，本发明通过基于靶标的光学畸变校正，结合平面激光诱导荧光技术，实现了周向液膜特征的光学畸变校正与高精度测量。设计了新型圆形校准靶标，基于靶标建立畸变图像与真实图像之间的变换关系，并通过液膜图像预处理、畸变校正及图像融合，即使不采用方形管道或增加水套装置，仍然可以有效重建管道内液膜图像，进而高精度地提取周向液膜时空流动特征。
[0026] 具体校正及测量方法如下：
[0027] 1、设计适用于管道内周向液膜测量的圆形校准靶标，采集液相工质条件下待测周向截面的靶标畸变图像，提取靶标畸变图像上的特征点；
[0028] 2、将靶标畸变图像与理想靶标图像上的特征点进行匹配，基于多项式函数建立原始图像与畸变图像特征点之间的映射数学模型；
[0029] 3、调整片状激光，使其水平照亮靶标平面处的管道横截面，激发溶解有特定荧光溶剂的轴向液膜产生荧光，利用高速摄像机摄取圆管内周向液膜的双视角荧光图像；
[0030] 4、将采集的液膜原始图像进行灰度化、差影及滤波等预处理，并基于映射数学模型对预处理后的液膜图像进行畸变校正，具体分为两步，一是对原畸变图像与畸变校正后图像进行映射变换，二是基于灰度插值确定畸变校正后图像的灰度值；
[0031] 5、将畸变校正后的双视角图像中心位置重合，并对摄取的左、右半周液膜图像有效部分进行融合，准确还原周向液膜的流动状态，以便提取周向液膜特征，并进一步分析管道内液膜的演化特性及传热传质特性等。
[0032] 本发明一种基于靶标校准的管道内周向液膜流动参数光学畸变校正和测量方法，其独特之处在于：不依赖于管道的材料、形状、尺寸以及工作液体的折射率,对于圆形管道而言，即使不外加方形水套也可以实现周向液膜厚度、分布等流动特征的无扰、准确测量和分析，对于非轴对称气液环状流流动结构和流动特性的研究而言尤为重要。
[0033] 下面结合附图对本发明做进一步的说明：
[0034] 图1所示为管道内周向液膜的不同测量方法示意图。以竖直管道内环状流液膜为例，圆形管道内为被荧光剂(如罗丹明B)
染色后的液膜，激光平面水平照射管道内液膜流场，两台高速摄像机分别从两个视角获取液膜左、右半周液膜图像并进行处理。由于管道以及气液相折射率的差异，以及管道弯曲形状所带来的影响，管道内液膜会产生光路畸变和失真，传统的测量方法需要在圆管外安装方形透明管段，其间充满折射率相近的工作液体(如水，其折射率n＝1.33)来进行光路补偿，如图1(a)所示。本发明方法基于靶标校准建立畸变图像与真实图像之间的关系，不加装方形水套也可以实现周向液膜特征的准确测量，如图1(b)所示。
[0035] 图2所示为设计的靶标样式及高速摄像机采集得到的靶标图像。靶标设计为圆形，直径为管道内径2r。为了方便、快速、准确地提取特征进行投影变换，设计靶标如图2(a)所示。靶标特征由距离靶心不同远近的圆点组成，中心圆点直径为1mm，距离其最近且与竖直方向夹角分别为0°、90°、180°及270°处分别排列4个小圆点，直径为0.25mm，如图2(a)中放大图所示，由于管道中心处畸变较小，可将其作为投影基准，其余圆点按横向与纵向间距1.5mm依次排列，直径为0.5mm。根据双视角周向液膜的分布特性，为了严格同步及有效利用高速摄像机视场，采用单台高速摄像机和折光分光光路虚拟两台高速摄像机，将靶标平面置于片状激光照亮的待测圆周截面位置上，在管道内加注液相，摄取的双视角靶标图像如图2(b)所示。从图中可以看出，由于管道形状及不同折射率的影响，靶标图像产生了明显的畸变，中心区域畸变较小，靠近管道处畸变较大。
[0036] 将靶标的圆点质心作为待提取特征点，根据光路分析，圆管内靶标中心区域的光路畸变最小，因此可以将中心区域处4个小圆点的位置关系作为已知条件。对靶标畸变图像进行灰度化、差影、滤波以及阈值分割等，提取各圆点特征的质心及面积，根据面积属性获取靶标中心位姿关系，并利用已知4个小圆点与中心点之间的位置关系，由近及远依次推导出其余圆点的质心坐标。建立理想靶标图像，并得到其特征点坐标及其位姿关系，根据畸变靶标图像与理想图像特征点之间的畸变投影关系匹配特征点，基于多项式函数建模的方法得到畸变校正的数学模型，其原理如下：
[0037] 设p(u,v)为理想图像，p(x,y)为对应畸变的图像，两者的特征点坐标存在非线性对应关系F，即
[0038] (x,y)＝F(u,v) (1)
[0039] 上式可用多项式模型来近似，即
[0040]
[0041] 式中，aij和bij为待定系数，n为多项式次数。通过特征点匹配获得多组匹配对，基于最小二乘拟合得到式(2)中多项式系数，为了满足计算量与畸变校正精度的需求，多项式次数n一般取5。
[0042] 基于靶标校正后，利用激光照射该横截面流场，采用高速摄像机摄取并对畸变的液膜荧光图像进行校正处理，即对图像上每一坐标值(u,v)进行变换F，得到在畸变图像上的对应坐标(x,y)。在获得畸变图像与校正图像坐标之间的变换关系后，采用灰度插值法确定校正后图像灰度值，
[0043] 图3为采用本发明方法采集得到的周向液膜原始图像以及畸变校正后的图像。其中，图3(a)为基于PLIF摄取的环状流圆管内周向液膜双视角图像，并经过了灰度化、差影及中值滤波等图像预处理。从图中可以看出，由于光路畸变的影响，摄取的液膜图像产生了很大的畸变。基于畸变校正数学模型对图3(a)进行处理，得到畸变校正后的液膜图像，如图3(b)所示。从图中可以看出，基于靶标的光学畸变校正方法可以很好地对管道内液膜畸变进行逆变换，从而有效地去除由于圆形管道及不同折射率所带来的畸变影响。将畸变校正后的双视角图像中心位置重合，并对图像中液膜有效部分进行逻辑“与”运算，融合后的周向液膜图像如图3(c)所示，准确还原了周向液膜的流动结构。
[0044] 基于靶标校准后的周向液膜图像，可以有效还原管道内周向液膜流动分布，准确获取周向液膜厚度及分布，进而获得液膜时空统计特征及
波动特征。该发明为科学研究和工程应用领域中实现管道内液膜的准确、普适测量，深入认识液膜演化及传热传质特性，以及气液环状流流动特性等提供了新的技术途径，有力推动了两相/多相流动测试和流动理论的发展。
