一种基于双目摄像头的气体泄漏量测量方法与流程

1.本发明属于工业气体管网领域，尤其涉及一种基于双目摄像头的气体泄漏量测量方法。


背景技术：

2.目前发电厂、化工厂、天然气管网等领域有大量输送工业气体的管道，部分管道使用年限较长、管道因腐蚀穿孔或连接处焊接质量差，导致泄漏事故时有发生。考虑到工业气体大多属于易燃易爆气体或高温高压气体，一旦发生泄漏，极易引发火灾、爆炸等重大安全事故，甚至造成人员伤亡，对公共安全造成极大威胁。因此对工业管道泄漏定位及泄漏速率测量显得尤为重要。由于工业管道周围环境复杂危险，存在管道处于高位、多条管道并排布置、管道周围高温、管道周围存在有毒有害气体等众多不利于测试人员靠近的情况。


技术实现要素：

3.本发明的目的是克服现有技术中的不足，提供一种基于双目摄像头的气体泄漏量测量方法。
4.这种基于双目摄像头的气体泄漏量测量方法，包括以下步骤：
5.步骤1、搭建气体泄漏量测量硬件平台，硬件平台包括双目摄像头、气体泄漏探测器和主机；双目摄像头为两台相同的相机l和相机r；
6.步骤2、基于双目摄像头测得测试点和管道待测泄漏点的距离；设待测泄漏点p在相机l和相机r上的成像点分别为pl和pr；相机l和相机r的焦距相同，均为f；相机l和相机r的之间的中心距为b；设测试点和待测泄漏点p的距离为z；
7.根据相似三角形原理：
[0008][0009]
推导得：
[0010][0011]
上式中，x
l
为相机l的x坐标，xr为相机r的x坐标；x
l-xr为视差，视差是同一个空间点在相机l和相机r成像中对应的x坐标的差值；当相机l和相机r的相对位置处于水平时，将相机l和相机r得到的rgb图编码成灰度图，来反映测试点和待测泄漏点p的距离z的远近，离镜头越近的在灰度图中越亮；因此通过双目摄像原理可以测得泄漏点离测量点的实际距离；
[0012]
步骤3、基于阵列测量的波束形成信号处理技术测得待测泄漏点p的噪声信号强度；
[0013]
步骤4、根据已知管道内气体成分相关参数和气体泄漏探测器测得管道泄漏点处的气体泄漏量，根据能量守恒方程和动量守恒方程定量描述绝热流动过程，来建立泄漏速
率模型，求得气体泄漏量。
[0014]
作为优选，步骤1中相机l和相机r集成在同一硬件平台上，相机l和相机r之间保持距离b；相机l和相机r用于远距离复杂环境中气体泄漏定位；气体泄漏探测器用于进行气体泄漏量测量；主机用于存储和显示相机l和相机r传输的图像数据。
[0015]
作为优选，步骤2中焦距f和相机l和相机r的中心距b通过标定得到。
[0016]
作为优选，步骤3具体包括以下步骤：
[0017]
步骤3.1、根据各个声音传声器接收声音信号的时间差以及声音传声器的位置来获取高度方向性波束特征，识别声源；
[0018]
步骤3.2、对声音传声器阵列中各阵元的输出进行延迟，补偿各阵元上的传播延时；然后通过加运算和等运算，在所有声音传声器感受到指定声源的同一个瞬时波前，产生一个空间响应极大值，将各阵元的输出信号聚焦到声源的真实位置；具体公式为：
[0019][0020][0021]
上式中，是聚焦位置为信号角频率为ω时声音传声器阵列中阵元的输出；m是声音传声器个数；pm(ω)是第m号声音传声器接收到的声音信号；为聚焦位置处的假设点声源辐射声波传播至第m号声音传声器时相对参考点的延时量；为第m号声音传声器位置向量；c为声速，j为虚数单位。
[0022]
作为优选，步骤3.1中声音传声器与相机l、相机r集成在同一硬件平台上。
[0023]
作为优选，步骤4具体为：
[0024]
建立泄漏速率模型：
[0025][0026]
上式中，κ为气体热容比；m为气体摩尔质量，单位为g/mol；g为单位体积的气体质量流量，单位为kg/(m2s)；g为管道直径，单位为m；le为泄漏点距离上游起点的距离，单位为m；f
′
为摩擦系数；r为气体常数；p1、t1、u1、ρ1分别表示上游起点的压力、温度、流速和气体密度；p2、t2、u2、ρ2分别表示与泄漏点同一直线的管内一点的压力、温度、流速和气体密度；p3、t3、u3、ρ3分别表示泄漏点的压力、温度、流速和气体密度；
[0027]
气体的状态方程为：
[0028]
pv＝znrt
[0029]
上式中，p为管道内压力；v为气体体积；r为气体常数；t为气体温度；z为气体的压缩因子；n为气体分子物质的量；
[0030]
气体的连续性方程为：
[0031]
u1a1ρ1＝u2a2ρ2[0032]
上式中，u1、a1、ρ1分别表示上游起点气体的流速、管道截面积和气体密度；u2、a2、ρ2分别表示与泄漏点同一直线的管内一点的气体的流速、管道截面积和气体密度；
[0033]
将气体的状态方程和连续性方程代入泄漏速率模型，求得泄漏孔处的气体泄漏速率q为：
[0034][0035]
上式中，q为泄漏速率，单位为kg/s；c为流量系数；a为泄漏孔的面积，单位为m2；pa为大气压力，单位为帕；流量系数c是一个和雷诺数re相关的系数，可以查表得到；雷诺数re是表征流体流动情况的无量纲数，通过流体的流速、密度、黏性系数计算得到，和流体本身性质相关。
[0036]
本发明的有益效果是：本发明通过一种基于双目摄像头的气体泄漏量测量方法，基于双目摄像头测距原理，错位布置双摄像头，计算出待测点距离，在复杂危险的工业环境中实现远距离非接触测量及定位；可实现远距离测量复杂环境中管道气体泄漏，解决工业管道测量环境复杂危险，甚至存在多管道遮挡、重叠等问题；保障测量人员安全。同时管道气体泄漏模型可计算管道气体泄漏速率。
附图说明
[0037]
图1为双目摄像头测距原理示意图；
[0038]
图2为气体管道孔隙泄漏的示意图。
具体实施方式
[0039]
下面结合实施例对本发明做进一步描述。下述实施例的说明只是用于帮助理解本发明。应当指出，对于本技术领域的普通人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干修饰，这些改进和修饰也落入本发明权利要求的保护范围内。
[0040]
实施例一
[0041]
本技术实施例一提供了一种如图1所示用于气体泄漏量测量的硬件平台，包括两台相同的相机l和相机r、气体泄漏探测器和主机；相机l和相机r用于远距离复杂环境中气体泄漏定位；气体泄漏探测器用于进行气体泄漏量测量；主机用于存储和显示相机l和相机r传输的图像数据。
[0042]
实施例二
[0043]
在实施例一的基础上，本技术实施例二提供了采用实施例1中硬件平台的基于双目摄像头的气体泄漏量测量方法：
[0044]
步骤1、搭建气体泄漏量测量硬件平台，硬件平台包括双目摄像头、气体泄漏探测器和主机；双目摄像头为两台相同的相机l和相机r；
[0045]
步骤2、基于双目摄像头测得测试点和管道待测泄漏点的距离；设待测泄漏点p在相机l和相机r上的成像点分别为pl和pr；相机l和相机r的焦距相同，均为f；相机l和相机r的之间的中心距为b；设测试点和待测泄漏点p的距离为z；
[0046]
根据相似三角形原理：
[0047]
[0048]
推导得：
[0049][0050]
上式中，x
l
为相机l的x坐标，xr为相机r的x坐标；x
l-xr为视差，视差是同一个空间点在相机l和相机r成像中对应的x坐标的差值；当相机l和相机r的相对位置处于水平时，将相机l和相机r得到的rgb图编码成灰度图，来反映测试点和待测泄漏点p的距离z的远近，离镜头越近的在灰度图中越亮；因此通过双目摄像原理可以测得泄漏点离测量点的实际距离；
[0051]
步骤3、基于阵列测量的波束形成信号处理技术测得待测泄漏点p的噪声信号强度；
[0052]
步骤3.1、根据各个传声器接收声音信号的时间差以及传声器的位置来获取高度方向性波束特征，识别声源；
[0053]
步骤3.2、对传声器阵列中各阵元的输出进行延迟，补偿各阵元上的传播延时；然后通过加运算和等运算，在所有传声器感受到指定声源的同一个瞬时波前，产生一个空间响应极大值(主瓣)，将各阵元的输出信号聚焦到声源的真实位置；具体公式为：
[0054][0055][0056]
上式中，是聚焦位置为信号角频率为ω时传声器阵列中阵元的输出；m是传声器个数；pm(ω)是第m号传声器接收到的声音信号；为聚焦位置处的假设点声源辐射声波传播至第m号传声器时相对参考点的延时量；为第m号传声器位置向量；c为声速，j为虚数单位；
[0057]
步骤4、根据已知管道内气体成分相关参数和气体泄漏探测器测得管道泄漏点处的气体泄漏量，根据能量守恒方程和动量守恒方程定量描述绝热流动过程，来建立泄漏速率模型，求得气体泄漏量；
[0058]
如图2所示，建立泄漏速率模型：
[0059][0060]
上式中，κ为气体热容比；m为气体摩尔质量，单位为g/mol；g为单位体积的气体质量流量，单位为kg/(m2s)；g为管道直径，单位为m；le为泄漏点距离上游起点的距离，单位为m；f
′
为摩擦系数；r为气体常数；p1、t1、u1、ρ1分别表示上游起点的压力、温度、流速和密度；p2、t2、u2、ρ2分别表示与泄漏点同一直线的管内一点的压力、温度、流速和密度；p3、t3、u3、ρ3分别表示泄漏点的压力、温度、流速和密度；
[0061]
气体的状态方程为：
[0062]
pv＝znrt
[0063]
上式中，p为管道内压力；v为气体体积；r为气体常数；t为气体温度；z为气体的压
缩因子；n为气体分子物质的量；
[0064]
气体的连续性方程为：
[0065]
u1a1ρ1＝u2a2ρ2[0066]
上式中，u1、a1、ρ1分别表示上游起点气体的流速、管道截面积和气体密度；u2、a2、ρ2分别表示与泄漏点同一直线的管内一点的气体的流速、管道截面积和气体密度；
[0067]
将气体的状态方程和连续性方程代入泄漏速率模型，求得泄漏孔处的气体泄漏速率q为：
[0068][0069]
上式中，q为泄漏速率，单位为kg/s；c为流量系数；a为泄漏孔的面积，单位为m2；pa为大气压力，单位为帕。