一种基于ORB和RANSAC的油相流速测量方法技术

本发明专利技术属于油田生产技术领域,具体涉及一种基于ORB和RANSAC的油相流速测量方法，1、根据FAST算法处理图像的数据信息，检测图像中的特征点；2、确定特征点的方向θ，根据改进的BRIEF算法描述步骤二中检测的特征点，使其描述子具有稳定性；3、根据Hamming距离对采集的多帧图像中特征点进行匹配，初步提取匹配点对；4、利用改进的RANSAC算法对初步提取的匹配点对进行筛选，提取最终的匹配点对5、根据最终提取到的匹配点对的位置信息和图像间采集的时间差t计算所有特征粒子的平均运动速度

A method of oil phase velocity measurement based on orb and RANSAC

【技术实现步骤摘要】
一种基于ORB和RANSAC的油相流速测量方法
：本专利技术属于油田生产
,具体涉及一种基于ORB和RANSAC的油相流速测量方法。
技术介绍
：在油田生产和化工领域，油水两相流是一种最为常见的两相流。目前油水两相流的主要特点是低流速、高含水，流型较为复杂，使其流动特性与其他条件下流体的流动特性不同，加大了流速等流动参数测量的难度。油田产液油的总量估算和油田设备的安全操控与流动参数的准确测量至关重要。粒子追踪技术(PTV)是一种常见的多相流流速测量技术，其关键是追踪单个粒子的运动轨迹，即对单个粒子进行匹配。粒子匹配的准确性直接影响最终多相流流速等流动参数的准确测量。在油水两相流中，水是无色透明的，低流速情况下无法观测到水的流动特征，因此，粒子的匹配即是针对油相进行匹配的，获取最终的平均速度亦为油相运动的平均速度。ORB特征提取与匹配方法是目前特征匹配方法中计算速度较快的一种匹配方法，其后续计算过程中所需的存储空间也相对较少。RANSAC算法可以从存在错误数据中提取正确数据，但其提炼的数据较少。若结合ORB特征提取及匹配和RANSAC算法，并对其进行改进，提高算法的计算速度，增多提炼的正确数据，并将其算法应用到PTV算法过程中，用于测量油水两相流油相流速，使其油相速度更为精确。
技术实现思路
：本专利技术的目的是为解决低流速高含水油水两相流油相流速测量的过程中，由于图像粒子匹配受图像模糊、亮度不均等影响流速测量精度的问题，提供一种基于改进的ORB和RANSAC算法的油水两相流油相流速测量方法。本专利技术采用的技术方案为：一种基于ORB和RANSAC的油相流速测量方法，所述测量方法包括以下步骤：步骤一：根据FAST算法处理图像的数据信息，检测图像中的特征点；所述FAST算法计算方法为：判断某一像素点O(i,j)是否为特征点，将以点O为中心，半径为3的周围像素点对应的像素值大小与点O处的像素值大小进行比较，若超过12个像素点对应的像素值与点O的像素值差的绝对值大于设定的阈值则判定点O为特征点，否则点O不是特征点；判断点O(i,j)是否为特征点的表达式为式中，若E为1时，判定点O为特征点，否则点O不是特征点，其中Dij表达式为式中TN为设定的阈值，amn为矩阵A的元素，矩阵A的表达式为I(i,j)是指点O(i,j)处的像素值大小；步骤二：确定特征点的方向θ，根据改进的BRIEF算法描述步骤二中检测的特征点，使其描述子具有稳定性；所述的确定特征点方向的方法为：若经步骤一判定点O(i,j)为特征点，则确定以点O(i,j)为中心，半径为r的所有像素点像素值的质心M(xm，ym)，质心M(xm,ym)的计算公式为式中x,y分别表示以点O为中心，半径为r内像素点对应的横纵坐标值；I(x,y)为点(x,y)处的像素值大小，特征点O的方向θ计算公式为所述改进的BRIEF算法的计算方法为：在点O周围的31×31区域内随机抽取l+1个随机点，将其定义为矩阵S，矩阵S的表达式为为使矩阵S对于特征点固定且唯一，将矩阵S进行旋转θ得Sθ，Sθ的计算公式为式中Rθ为特征点方向θ对应的旋转矩阵，Pl是指随机点旋转后的位置；比较矩阵Sθ中Pk和Pk+1(1≤k≤l)代表的点像素值的大小，记录对应的二值码，构成的二值码串即为点O的描述子，点O描述子fO的表达式为式中k是描述子中一二值码对应的序号，fO的初始值为0，I(Pk)是指矩阵Pk代表点处的像素值大小；步骤三：根据Hamming距离对采集的多帧图像中特征点进行匹配，初步提取匹配点对；所述匹配点对粗提取的计算方法为式中G表示特征点是否匹配成功，若G＝1,则表示匹配成功，否则表示匹配失败；fO1、fO2分别表示第一帧和第二帧图像中特征点O的描述子，q表示图像中特征点的序号，s表示图像中检测到特征点的个数；粗提取后的匹配点对用矩阵H表示，H表示为式中z表示粗提取匹配点对的个数，x1z,y1z分别是到第一帧图像中经粗提取得到的特征点的横纵坐标，x2z,y2z分别是到第二帧图像中经粗提取得到的特征点横纵坐标，H1z和H2z分别表示为步骤四：利用改进的RANSAC算法对初步提取的匹配点对进行筛选，提取最终的匹配点对；所述改进的RANSAC算法的计算方法为流程1：设置Hamming距离一个阈值，提取低于阈值的匹配点对得到精提取的匹配点对，标记为矩阵I，I表示为式中z′表示为筛选后匹配点对的个数。x1z′,y1z′分别是到第一帧图像中经粗提取得到的特征点的横纵坐标，x2z′,y2z′分别是到第二帧图像中经粗提取得到的特征点横纵坐标；I1z′和I2z′分别表示为改进RANSAC算法的观测数据和模型计算数据分别为矩阵H和矩阵I；流程2：匹配点对的坐标H和I均需进行归一化，以减少噪声的干扰，计算公式为式中Prv、Mrv(r,v＝1,2)分别为第v帧归一化后粗提取(r＝1)或精提取(r＝2)的特征点的坐标和归一化转换矩阵；Mrv(r,v＝1,2)计算公式为式中分别为第v帧归一化前粗提取(r＝1)或精提取(r＝2)的特征点的横纵坐标的平均值，λrv(r,v＝1,2)分别为第v帧粗提取(r＝1)或精提取(r＝2)的特征点的归一化伸缩比，λrv(r,v＝1,2)的计算公式为式中u,u′分别是指归一化前粗提取(r＝1)或精提取(r＝2)的图像中特征点的序号；流程3：根据精提取后的匹配点对P21、P22估计模型，获取基础矩阵F′，基础矩阵F′的计算公式为式中是P22的转置，对基础矩阵F′进行反归一化处理获取符合观测数据的基础矩阵F：式中表示先对矩阵M12进行求逆，后进行转置；是M22的转置，表示矩阵M11的逆；流程4：计算观测数据H到对极线的距离，低于设定的阈值TN判定为内点，记录内点数β，β的计算公式为式中P11x、P11y、P12x、P12y分别为P11、P12的第一行和第二行，w是矩阵P11和P12列的序号，P11w、P12w分别为P11、P12的第w列。P11和P12的改写方式为流程5：重复流程3和流程4，选取最大内点数βmax以及对应的内点，内点即为最终的匹配数据。内点的表达公式为式中Pin1、Pin2分别是第一帧和第二帧图像已成功匹配的特征点。x1′、y1′分别为第一帧图像已成功匹配特征点的横纵坐标，x1″、y1″分别为特征点(x1′,y1′)对应第二帧图像中特征点的横纵坐标；步骤五：根据最终提取到的匹配点对的位置信息和图像间采集的时间差t计算所有特征粒子的平均运动速度即为油相的平均运动速度所述油相的平均运动速度计算公式为：式中h表示最终匹配点对的序号，t是指两本文档来自技高网...

【技术保护点】
1.一种基于ORB和RANSAC的油相流速测量方法，其特征在于：所述测量方法包括以下步骤：/n步骤一：根据FAST算法处理图像的数据信息，检测图像中的特征点；/n所述FAST算法计算方法为：/n判断某一像素点O(i,j)是否为特征点，将以点O为中心，半径为3的周围像素点对应的像素值大小与点O处的像素值大小进行比较，若超过12个像素点对应的像素值与点O的像素值差的绝对值大于设定的阈值则判定点O为特征点，否则点O不是特征点；/n判断点O(i,j)是否为特征点的表达式为/n

【技术特征摘要】
1.一种基于ORB和RANSAC的油相流速测量方法，其特征在于：所述测量方法包括以下步骤：
步骤一：根据FAST算法处理图像的数据信息，检测图像中的特征点；
所述FAST算法计算方法为：
判断某一像素点O(i,j)是否为特征点，将以点O为中心，半径为3的周围像素点对应的像素值大小与点O处的像素值大小进行比较，若超过12个像素点对应的像素值与点O的像素值差的绝对值大于设定的阈值则判定点O为特征点，否则点O不是特征点；
判断点O(i,j)是否为特征点的表达式为



式中，若E为1时，判定点O为特征点，否则点O不是特征点，其中Dij表达式为



式中TN为设定的阈值，amn为矩阵A的元素，矩阵A的表达式为I(i,j)是指点O(i,j)处的像素值大小；
步骤二：确定特征点的方向q，根据改进的BRIEF算法描述步骤二中检测的特征点，使其描述子具有稳定性；
所述的确定特征点方向的方法为：
若经步骤一判定点O(i,j)为特征点，则确定以点O(i,j)为中心，半径为r的所有像素点像素值的质心M(xm，ym)，质心M(xm,ym)的计算公式为



式中x,y分别表示以点O为中心，半径为r内像素点对应的横纵坐标值；I(x,y)为点(x,y)处的像素值大小，特征点O的方向q计算公式为



所述改进的BRIEF算法的计算方法为：
在点O周围的31×31区域内随机抽取l+1个随机点，将其定义为矩阵S，矩阵S的表达式为



为使矩阵S对于特征点固定且唯一，将矩阵S进行旋转q得Sq，Sq的计算公式为



式中Rq为特征点方向q对应的旋转矩阵，Pl是指随机点旋转后的位置；比较矩阵Sq中Pk和Pk+1(1≤k≤l)代表的点像素值的大小，记录对应的二值码，构成的二值码串即为点O的描述子，点O描述子fO的表达式为



式中k是描述子中一二值码对应的序号，fO的初始值为0，I(Pk)是指矩阵Pk代表点处的像素值大小；
步骤三：根据Hamming距离对采集的多帧图像中特征点进行匹配，初步提取匹配点对；
所述匹配点对粗提取的计算方法为



式中G表示特征点是否匹配成功，若G＝1,则表示匹配成功，否则表示匹配失败；fO1、fO2分别表示第一帧和第二帧图像中特征点O的描述子，q表示图像中特征点的序号，s表示图像中检测到特征点的个数；粗提取后的匹配点对用矩阵H表示，H表示为



式中z表示粗提取匹配点对的个数，x1z,y1z分别是到第一帧图像中经粗提取得到的特征点的横纵坐标，x2z,y2z分别是到第二帧图像中经粗提取得到的特征点横纵坐标，H1z和H2z分别表示为



步骤四：利用改...

【专利技术属性】
技术研发人员：韩连福，王海霞，付长凤，刘兴斌，谢荣华，刘超，
申请(专利权)人：东北石油大学，
类型：发明
国别省市：黑龙;23