技术领域
[0001] 本
发明属于气液两相流动测量技术领域，具体涉及一种基于光场成像技术的气液两相流中气泡三维测量系统及其方法。
背景技术
[0002] 气液两相流动作为两相流中一种常见的流动形式，广泛存在于化工、
能源、环保和轻工业等行业中。随着科学技术的进步，气液两相流在科学研究、环境保护、工业生产及人类生活中日益重要。气泡作为气液两相流中的分散相，其参数是表征气液两相流动的重要参数之一。气泡的测量，有助于进一步理解气液两相流的内部流动机理，也有助于揭示两相系统内部动量、
质量和
能量传递规律，最终对于指导两相流实验系统的设计、两相流设备的运行优化和有效控制各种化学反应过程产生重要意义。因此，气泡测量是两相流研究的一个重要课题。
[0003] 对气液两相流中气泡的测量有
接触式和非接触式两种方法，接触式测量方法的通常做法为将
传感器浸入到液体中，利用气泡穿过时传感器
探头产生的反馈
信号来计算气泡的尺寸和速度，由于接触式测量会
对流场本身产生干扰，并且传感器数量有限，所以此种测量方法的准确性受到了大大的限制。非接触式测量方法由于其无干扰，测量
精度高等优点，因而得到了迅速的发展，目前常用的方法有激光多普勒技术、
核磁共振技术、过程
层析成像技术和图像测量技术等。
[0004] 随着光学技术、计算机技术以及工业CCD的快速发展，图像测量技术已经成为气泡测量研究的热点，通过对相机拍摄得到的气泡图像进行后期处理，可以得到气泡的物理参数以及空间分布，并且具有无干扰、瞬态、全场
可视化等优点。但是，传统图像测量技术只能记录气泡二维信息，即单个截面的气泡信息，无法获得整个测量体内气泡三维信息。
[0005] 光场成像技术是一种全新的成像方法，通过在相机镜头和CCD传感器之间放置微透镜阵列，光场相机能同时记录光线的强度信息和
角度信息，克服了传统成像只能记录光线强度的缺点，光场成像技术目前已经应用到3D复原、火焰测量和流场测速等领域，但还无法应用于立体三维气泡的测量中。
发明内容
[0006] 本发明所要解决的技术问题是针对上述
现有技术存在的不足，而提供一种可实现气泡
三维重建，进而统计气泡群尺寸分布、空间
位置及含气率等参数的基于光场成像技术的气液两相流中气泡三维测量系统及其方法。
[0007] 本发明基于光场成像技术，通过单相机单次曝光将被测流场的光场数据都储存在一张原始光场图像上，结合计算成像技术获取被测流场的全聚焦图像以及不同深度位置的重聚焦
图像序列，通过系统标定获取重聚焦图像序列的深度以及放大率函数，对全聚焦图像进行处理获得气泡的深度方向上的投影图像，对重聚焦图像序列进行清晰度评价得到气泡深度信息，基于获取的气泡投影图像、深度以及放大率信息，计算得到气泡三维坐标以及当量直径，进而实现气泡三维重建，并进一步统计分析两相流场的气泡群尺寸分布、空间位置及含气率等参数。
[0008] 为解决上述技术问题，本发明采用的技术方案是：
[0009] 一种气液两相流中气泡测量系统，包括流场产生装置、LED平面
光源、光场相机以及数据分析储存装置；其特征在于，
[0010] 所述流场产生装置，用于产生气液两相流场，并可通过控制
阀门开度获取不同工况下的流场；
[0011] 所述LED平面光源，用于为气液两相流场提供照明，降低相机曝光时间，消除流场图像中气泡拖影等负面效果；
[0012] 所述光场相机，用于获取气液两相流场的光场信息，并将所采集的光场信息传输至数据分析储存装置；
[0013] 所述数据分析储存装置，对光场相机所记录的流场的光场信息进行处理，得到流场的全聚焦图像和聚焦在不同深度位置的重聚焦图像序列；对所得到的全聚焦图像进行二值化处理，得到气泡图像质心二维坐标以及当量直径；对重聚焦图像序列进行轴向清晰度比对结合系统标定，得到物空间气泡的质心深度及该深度处放大率，进而计算得到物空间气泡质心三维坐标以及当量直径。
[0014] 所述数据分析储存装置得到气泡图像质心二维坐标以及当量直径的方法是：将全聚焦图像转化为灰度图，进行二值化分割和
像素填充，得到气泡投影图像，基于已知的像素尺寸，计算气泡图像质心二维坐标和当量直径。
[0015] 所述数据分析储存装置得到物空间气泡的质心深度的方法是：
[0016] 将重聚焦图像序列中的每一幅图像划分成M×N个正方形网格，基于清晰度评价函数计算每个网格的清晰度评价值，对于M和N，有
[0017]
[0018]
[0019] 式中：l和m分别为为重聚焦图像高度和
水平方向像素数，np为单个窗口边长像素数，为保证清晰度评价的准确性，np取值为7～11。
[0020] 轴向对比所有重聚焦图像中同一处网格的清晰度评价值，最大值所在的位置便为该网格的深度；
[0021] 对所有网格进行运算，得到重聚焦图像的像素
深度图，将像素深度图与全聚焦二值化图取交集得到气泡边缘深度图，取每个气泡边缘深度的平均值作为该气泡的质心深度。
[0022] 一种气液两相流中气泡三维测量方法，其特征在于，包括以下步骤：
[0023] 步骤一：采用光场相机对气液两相流场进行拍摄，获取流场的光场信息，获取流场全聚焦图像在不同深度位置聚焦的流场图像序列；
[0024] 步骤二：对流场全聚焦图像进行二值化处理，得到像空间气泡质心二维坐标以及当量直径；
[0025] 步骤三：将重聚焦图像序列各对应像素进行轴向清晰度比对，获取全聚焦图像深度图；在全聚焦图像深度图中，计算气泡质心深度位置；
[0026] 步骤四：根据气泡三维重建模型，结合步骤二得到的像空间气泡质心二维坐标以及当量直径及步骤三得到的气泡质心深度位置和对应放大率，对气泡进行三维重建。
[0027] 还包括，
[0028] 步骤五、气泡进行三维重建后，根据重建结果，得到气泡群尺寸分布、空间位置及含气率数。
[0029] 所述步骤四中，所述气泡三维重建模型由物空间气泡参数与像空间气泡参数之间的关系建立，物空间气泡三维重建参数的表达式如下：
[0030]
[0031] M＝g(zc)
[0032] xc＝x′c/M
[0033] yc＝y′c/M
[0034] A＝A′/M2
[0035] de＝d′e/M
[0036] 式中：d(x，y)为物空间气泡投影面各处深度，xc、yc、zc为物空间气泡质心三维坐标，x′c、y′c为气泡全聚焦图像中气泡质心二维坐标；A为气泡投影面积，A′为全聚焦图像中气泡面积，de为气泡当量直径，d′e为全聚焦图像中气泡当量直径；M是放大率。
[0037] 所述步骤三中，基于点锐度清晰度评价函数进行清晰度比对，其中点锐度清晰度评价函数为：
[0038]
[0039] 式中：df表示图像灰度变化的幅值，dx表示像素之间的距离，df/dx计算时取该像素的8个邻域，I表示图像网格的像素数。
[0040] 步骤五：将重聚焦图像序列各对应像素进行轴向清晰度比对，获取全聚焦图像深度图，并根据步骤五中气泡投影图像计算气泡质心深度位置，得到物空间气泡质心三维坐标；
[0041] 所述像空间气泡位置及尺寸参数计算表达式如下：
[0042]
[0043]
[0044]
[0045]
[0046] 式中，x′c、y′c为像空间气泡二维坐标，A′为气泡图像面积，de为气泡图像等面积折合当量直径，xi，j、yi，j为CCD上像素中心坐标， 为气泡图像像素中心坐标平均值，a为单个像素面积。
[0047] 所述全聚焦像素深度图通过以下步骤获取：将重聚焦图像序列中的每一幅图像划分成469×313个网格，单个网格大小为7×7像素；计算每个网格的清晰度评价值，然后轴向对比所有重聚焦图像中同一处网格的清晰度评价值，最大值所在的位置便为该网格的深度；对所有网格进行运算，并以水箱内壁面作为z轴零点；
[0048] 根据流场全聚焦二值化图像获得气泡
边缘图像；根据全聚焦像素深度图获得气泡
深度图像；二者取交集获得单个气泡边缘深度图；气泡边缘深度平均值作为气泡质心深度。
[0049] 放大率M为深度z的函数，M＝g(z)，放大率函数M计算模型为：
[0050]
[0051]
[0052] u＝z+LOC
[0053]
[0054]
[0055]
[0056] 式中：u、v、f分别表示物距、像距、相机主镜头焦距，M(0)为z＝0处放大率，LOC为物空间
坐标系原点O到主透镜
光心C的距离。
[0057] 本发明首次将光场成像技术引入气泡测量领域中，利用单相机单次曝光来捕捉整个光场信息，通过后期
算法处理可以获得不同深度气泡图像，计算分析得到气泡位置及尺寸参数，进而重建气泡的三维分布。与传统成像比，无须多次曝光、改变相机位置或者布置多相机阵列，具有操作简单、实时性强、设备体积小和应用场合广等优点。
[0058] 有益效果：与现有技术对比有以下有点：本发明采用单个光场相机单次曝光即可获得流场光场信息，利用计算成像结合
图像处理技术可实现气泡三维重建及参数测量；测量系统设备简单，算法简洁，运算量小，适用场合广，且为非接触式测量，对流场不产生干扰。
附图说明
[0059] 图1测量系统结构图。
[0060] 图2数字重聚焦原理图
[0061] 图3重聚焦图像序列和全聚焦图像示意图。
[0062] 图4气泡三维重建模型图。
[0063] 图5重聚焦图像深度标定装置图。
[0064] 图6物空间深度与重聚焦图像序号标定结果
[0065] 图7放大率与深度关系曲线。
[0066] 图8气液两相流场全聚焦图像。
[0067] 图9气液两相流场二值化图像。
[0068] 图10全聚焦像素深度图。
[0069] 图11气泡三维重建效果图。
[0070] 图12气泡质心深度分布图。
[0071] 图13气泡当量直径分布图。
[0072] 图14截面含气率曲线。
具体实施方式
[0073] 下面结合附图和具体
实施例，进一步阐述本发明。应理解这些实施例仅用于说明本发明而不用于限制本发明的范围，在阅读本发明之后，本领域技术人员对本发明的各种等价形式的
修改落于本
申请所附
权利要求所限定的范围。本实施例中光场相机采用德国Raytrix公司的R29型光场相机。
[0074] 本发明一种基于光场成像技术的气液两相流中气泡三维测量方法，包括以下步骤：
[0075] 步骤一：利用光场相机拍摄气液两相流场，系统布置如图1所示，获取流场的光场信息。数字重聚焦技术过程如图2所示，其中，光线L(u，v，s，t)的方向由微透镜中心坐标(u，v)和CCD上像素坐标(s，t)唯一确定，光线强度值由像素记录，CCD与微透镜中心距离为D。基于获取的光线强度和方向信息，在空间任意位置放置成像面，即重聚焦面，重聚焦面与CCD的距离为d，在重聚焦面上对光线进行二重积分，得到重聚焦图像，积分表达式如式(1)所示。
[0076]
[0077] 式中：E(s′，t′)为重聚焦面坐标(s′，t′)处光线强度。
[0078] 利用数字重聚焦技术在不同位置成像可以得到重聚焦图像序列，过程如图3所示，[0079] 其中：z′1、z′2、z′3为重聚焦面示例位置，将所有重聚焦面上的清晰像素组合得到全聚焦图。
[0080] 步骤二：在图4中，B为物空间气泡，B的质心坐标为mc(xc，yc，zc)，面积为A，当量直径为de，B经过成像过程在像空间聚焦，得到全聚焦图像B′，放大率为M，B′的质心坐标为m′c(x′c，y′c，z′c)，面积为A′，当量直径为d′e，则物空间气泡参数的计算如式(2)-(7)所示。
[0081]
[0082] M＝g(zc) (3)
[0083] xc＝x′c/M (4)
[0084] yc＝y′c/M (5)
[0085] A＝A′/M2 (6)
[0086] de＝d′e/M (7)
[0087] 式中，d(x，y)为物空间气泡投影面各处深度，放大率M是物空间深度的函数，通过步骤四的标定工作确定，在像空间，则有
[0088]
[0089]
[0090] A′＝∫∫B′dx′dy′ (10)
[0091]
[0092] CCD记录的光线强度信息是离散的，式(8)-(11)的离散形式如式(12)-(15)所示，气泡图像质心坐标为各像素点中心坐标平均值，气泡图像面积为各像素点面积之和。
[0093]
[0094]
[0095]
[0096]
[0097] 式中，xi，j和yi，j为每个像素点中心坐标， 为气泡图像像素中心坐标平均值，a为每个像素的面积。
[0098] 步骤三：如图5所示，将标定尺沿深度方向放置，与水平面呈一定夹角θ，尺子底端与镜头前段距离为d0，用光场相机拍摄后重聚焦得到处于某一深度范围的图像序列，基于点锐度清晰度评价函数，确定每一张图片最清晰的刻度，清晰度评价函数表达式如式(16)所示。
[0099]
[0100] 式中：df表示图像灰度变化的幅值，dx表示像素之间的距离，df/dx计算时取该像素的8个邻域，I表示图像网格的像素数。Li值越大，代表图像越清晰。
[0101] 若某清晰点A刻度值为LA，则在物空间A点到镜头前段距离dA为：
[0102] dA＝d0+LAcosθ (17)
[0103] 由于标定尺最小刻度为1mm，因此标定的误差范围为(-0.5mm，0.5mm)，为减小标定误差，对同一毫米刻度的深度数据在±0.5mm的范围内进行线性插值，设
流体与水箱接触面深度为零，得到的实际深度与重聚焦图像序号的关系如图6所示。
[0104] 步骤四：为确定气泡实际尺寸与气泡像尺寸之间的关系，需要标定物空间不同深度的放大率M，M是深度z的函数，M＝g(z)。根据成像公式，
[0105]
[0106]
[0107] u＝z+LOC (20)
[0108] 其中，f为主透镜焦距，LOC为物空间坐标系原点O到主透镜光心C的距离，联立解得放大率M与物面深度z的关系式M＝g(z)如式(21)所示。
[0109]
[0110] 由于LOC未知，因此首先需要标定LOC的值。采用的标定方法为，将标定尺放置在z＝0的位置，即水箱内壁处，拍摄并由式(21)和(22)计算放大率M(0)及LOC。
[0111]
[0112]
[0113] 其中，L为标定尺实际长度，np为标定尺的像素长度，l为CCD单个像素长度。求得M(0)后代入式(22)，在z区间[0，20mm]内求得的M＝g(z)函数曲线如图7所示。
[0114] 步骤五：气液两相流场全聚焦图像如图8所示，通过
阈值分割获得其二值化图像，如图9所示，通过二值化图像可以获取气泡边缘数据，对二值化图像执行填充后，根据式(12)-(15)，计算像空间每一个气泡质心二维坐标x′c、y′c以及当量直径d′e。
[0115] 步骤六：将重聚焦图像序列中的每一幅图像划分网格，以单个网格尺寸7×7像素为例，计算每个网格的清晰度评价值，然后轴向对比所有重聚焦图像中同一处网格的清晰度评价值，最大值所在的位置便为该网格的深度。对所有网格进行相同运算，并以水箱内壁面作为z轴零点，得到如图10所示的像素深度图。将图9与图10取交集得到气泡边缘深度图，根据式(2)，取每个气泡边缘深度的平均值作为该气泡的质心深度zc，根据图7确定z＝zc处放大率M(zc)，并根据式(4)-(7)计算物空间气泡质心二维坐标xc、yc以及当量直径de。
[0116] 步骤七：基于获取的气泡质心三维坐标(xc，yc，zc)和当量直径de，对气泡进行三维重建。
[0117] 效果如图11所示，测量体的体含气率GVF可由式(24)计算。
[0118]
[0119] 式中，di是气泡当量直径，n为气泡个数，Vm为测量体体积。
[0120] 对气泡数据做进一步分析，可以分别得到所有气泡深度位置、当量直径分布和截面含气率曲线，如图12-14所示，气泡质心深度分布与截面含气率曲线一致，从而验证了测量方法的准确性。
