发明名称 一种基于超声波的管道流速成像和流量测量方法
摘要 一种基于超声波的管道流速成像和流量测量方法,加工制作具有三个压电单元组成的超声波探头,用于发射和接收超声波信号,将超声波探头分别布置在与管道轴线成45°的两个倾斜截面上,测量各声道上超声波顺流和逆流传播的飞行时间,计算管道流体在各声道上的平均流速:将每组声道等间隔内插,获取内插后各声道上的平均流速;基于平行束投影的滤波反投影方法从声道的平均流速数据出发重建流速分布剖面,基于重建得到的流速分布剖面图等间隔选取流速点,由各点流速简单平均并与管道截面积相乘即为管道流量数值,本发明可以实现在复杂流态下对管道内部流场分布的监控,减小超声波流量计的模型误差,最终获取精确的流量数值。
申请公布号 CN102607653B 申请公布日期 2013.09.04
申请号 CN201210100169.5 申请日期 2012.04.06
申请人 清华大学 发明人 王伯雄;崔园园;柳建楠;罗秀芝
分类号 G01F1/66(2006.01)I;G01P5/24(2006.01)I 主分类号 G01F1/66(2006.01)I
代理机构 西安智大知识产权代理事务所 61215 代理人 贾玉健
主权项 1.一种基于超声波的管道流速成像和流量测量方法,其特征在于,包括以下步骤:步骤一:加工制作具有三个压电单元组成的超声波探头,用于发射和接收超声波信号,每个压电单元均可完成发射和接收超声波任务,每个压电单元均有一定的声束角,声束角满足这样的条件:使得边界压电单元的声束角分别涉及相邻的四个探头,中心单元的声束角涉及正对面的三个探头;步骤二:将超声波探头分别布置在与管道轴线成45°的两个倾斜截面上,每个测量平面上均匀放置12个兼具发射和接收任务的超声波探头,探头在每个测量平面上形成66条声道;步骤三:分别测量各声道上超声波顺流和逆流传播的飞行时间,由下式计算管道流体在各声道上的平均流速:<maths num="0001"><![CDATA[<math><mrow><msub><mover><mi>V</mi><mo>&OverBar;</mo></mover><mi>path</mi></msub><mo>=</mo><mfrac><mi>L</mi><mrow><mn>2</mn><mi>cos</mi><mi>&phi;</mi></mrow></mfrac><mrow><mo>(</mo><mfrac><mrow><msub><mi>t</mi><mrow><mi>B</mi><mo>&RightArrow;</mo><mi>A</mi></mrow></msub><mo>-</mo><msub><mi>t</mi><mrow><mi>A</mi><mo>&RightArrow;</mo><mi>B</mi></mrow></msub></mrow><mrow><msub><mi>t</mi><mrow><mi>A</mi><mo>&RightArrow;</mo><mi>B</mi></mrow></msub><msub><mi>t</mi><mrow><mi>B</mi><mo>&RightArrow;</mo><mi>A</mi></mrow></msub></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,t<sub>B→A</sub>和t<sub>A→B</sub>分别为超声波顺流和逆流的飞行时间,<img file="FDA00003261623100018.GIF" wi="39" he="52" />为测量截面的倾角,L为声传播路径长度;步骤四:将每个测量平面上的声道分为12组,每组中的声道相互平行,将每组声道等间隔内插,并按步骤三的公式获取内插后各声道上的平均流速;步骤五:基于平行束投影的滤波反投影方法从声道的平均流速数据出发重建流速分布剖面,公式如下:<img file="FDA00003261623100012.GIF" wi="1396" he="109" />式中,<img file="FDA00003261623100013.GIF" wi="228" he="71" />为任意一点,<img file="FDA00003261623100014.GIF" wi="170" he="73" />代表过<img file="FDA00003261623100015.GIF" wi="42" he="56" />的任意直线L,r为原点到直线L的距离,<img file="FDA00003261623100016.GIF" wi="355" he="79" />为直线L的法向,<img file="FDA00003261623100017.GIF" wi="169" he="71" />为平行束投影采样<img file="FDA00003261623100021.GIF" wi="144" he="75" />关于第一个变量的傅里叶变换;该公式的具体实施过程为:(1)对平行束投影采样<img file="FDA00003261623100022.GIF" wi="148" he="76" />关于r做傅里叶变换,得到<img file="FDA00003261623100023.GIF" wi="194" he="74" />(2)进行频域滤波操作,滤波器传递函数为|σ|;(3)对滤波后的数据关于σ做反傅里叶变换,得到<img file="FDA00003261623100024.GIF" wi="187" he="80" />(4)对所有经过点<img file="FDA00003261623100025.GIF" wi="46" he="56" />的直线数据<img file="FDA00003261623100026.GIF" wi="154" he="72" />积分求和,即为该点的速度值<img file="FDA00003261623100027.GIF" wi="147" he="64" />步骤六:基于重建得到的流速分布剖面图等间隔选取流速点,由各点流速简单平均并与管道截面积相乘即为管道流量数值。
地址 100084 北京市海淀区100084信箱82分箱清华大学专利办公室