发明名称 低含气率的气液两相流流型识别方法
摘要 本发明属流体测量技术领域,涉及公开了一种低含气率的气液两相流流型识别方法。该方法首先应用科氏质量流量计测量气液两相流的瞬时流量波动数据,然后,采用经验模态分解方法提取气液两相流瞬时流量波动的IMF分量,并计算各IMF分量的能量占总能量的比例,最后,根据各IMF分量的能量比例识别低含气率气液两相流的流型。该流型识别方法简单、易于实现、实时性好。适用于低含气率的气液两相流流型识别。
申请公布号 CN103674478A 申请公布日期 2014.03.26
申请号 CN201310671989.4 申请日期 2013.12.11
申请人 中国石油大学(华东) 发明人 王微微;陈宇;梁霄;白明雷;陈静静
分类号 G01M10/00(2006.01)I 主分类号 G01M10/00(2006.01)I
代理机构 代理人
主权项 1.一种低含气率的气液两相流流型识别方法,其特征在于包括如下具体步骤: 步骤A:应用经验模态分解方法对流量波动信号进行分解,(1)找出信号x(t)中的所有局部极大值点,用3次样条函数连接成上包络,再用3次样条函数连接所有局部极小值点构成下包络,求上下包络线的均值m<sub>1</sub>(t),其中,x(t)为科氏质量流量计测得的质量流量瞬时波动数据;(2)求x(t)与m<sub>1</sub>(t)的差值d<sub>1</sub>(t),并计算d<sub>11</sub>(t)=d<sub>1</sub>(t)-m<sub>11</sub>(t),其中,m<sub>11</sub>(t)为x(t)的局部均值,该过程重复k次,得到d<sub>11</sub>(t),d<sub>12</sub>(t),…,d<sub>1(k-1)</sub>(t),d<sub>1k</sub>(t),如果满足<img file="FDA0000434763940000011.GIF" wi="723" he="177" />则c<sub>1</sub>(t)=d<sub>1k</sub>(t)就是由信号x(t)中得到的第一个IMF分量,其中,T为数据点数,ε为设置值;(3)计算剩余分量r<sub>1</sub>(t)=x(t)-c<sub>1</sub>(t),并将r<sub>1</sub>(t)作为新的波动信号x(t),重复步骤A中(2)的过程,依次得到c<sub>2</sub>(t),c<sub>3</sub>(t),…,c<sub>n</sub>(t)和<img file="FDA0000434763940000012.GIF" wi="600" he="317" />以上过程直到满足c<sub>n</sub>(t)≤δ或r<sub>n</sub>(t)≤δ或剩余分量r<sub>n</sub>(t)成为单调函数时结束,其中,δ为设置值;步骤B:根据<img file="FDA0000434763940000013.GIF" wi="842" he="246" />计算每个IMF分量所占的能量比例,式中,<img file="FDA0000434763940000016.GIF" wi="68" he="68" />表示第i个IMF分量c<sub>i</sub>的能量占总能量的比例;步骤C:根据<img file="FDA0000434763940000017.GIF" wi="65" he="68" />识别低含气率的气液两相流流型,如果满足<img file="FDA0000434763940000014.GIF" wi="281" he="146" />且<img file="FDA0000434763940000018.GIF" wi="367" he="79" />i=1,2,3,那么该流型为泡状流型;如果满足<img file="FDA0000434763940000019.GIF" wi="254" he="74" />那么该流型为过渡流型;如果满足<img file="FDA0000434763940000015.GIF" wi="284" he="149" />且<img file="FDA00004347639400000110.GIF" wi="239" he="77" />那么该流型为段塞流型。
地址 266580 山东省青岛市经济技术开发区长江西路66号