专利名称：基于超声导波技术的管道液体流量测量方法
技术领域：
本发明涉及基于超声导波技术的管道液体流量测量方法，属于测试计量 及无损检测技术领域。
背景技术：
在航空航天、核工业、电力、制冷、食品以及制药等行业，广泛使用小
管径(小于100mm，甚至小于10mm)管路用于传输液体。而流量作为一个重 要参数，是工艺流程中的必控参数之一。如何实现其实时、准确、方便的测 量，这是流量测量领域面临的难题之一。根据不同工况，采用不同类型的流 量计(诸如，差压式流量计、电磁式流量计、科氏质量流量计等等)进行流 量测量是现今工业界常常使用的方法，已经得到了较为广泛的应用。而基于 现场管路工况的复杂性、传输介质的多样性，以及管路安装的限制，对流量 计的技术指标的要求越来越高，适应范围越来越广，很难有一种流量计能满 足不同场合的需求。现有流量计的种类繁多，其工作原理多达10余种，类型 多达200多种，但仍不能完全满足各领域提出的新课题。近年来，人们仍然致 力于创新，不断推出新技术、新产品。
超声波流量计，与电磁流量计和科氏质量流量计等相同，是近来发展最 快，使用范围较广的流量计之一。但是，超声流量计，由于其原理上的原因， 也存在一些不足。在某些情况下，测量精度低，且很难有较大的提高；甚至 在一些特殊情况下，不能进行测量。如当管路的管径等于或小于超声波的波 长时，由于声程较短、声束旁瓣的干扰等因素，声波不可能以体波的形式传 播，因此，从原理上讲，已经不能采用常规超声体波流量计进行测量。而其他形式的小孔径流量计，如电磁流量计，由于其本身测量原理的原因，对于 输送非导磁介质的小管径管路，也不能完全代替超声流量计。因此，发展一 种新型的，适合管道尤其是小管径管道中液体流量测量方法显得十分重要。
目前，利用超声导波对管道中液体流量测量的研究已取得一定进展，证 实超声导波技术用于管道中液体流量测量的可行性和应用潜力。但目前取得
的研究成果较少且较为初步，仅见Haruminchi Sato等在2007年46巻第7B 期Japanese Journal of Applied Physics中发表的"Theoretical investigation of guide wave flowmeter",对超声导波技术用于充液管道液体流量测量进行了报 道。但该文中没有考虑超声导波的频散、能量分布等传播特性，对用于液体 流量检测的超声导波模态没有进行优化选取。文中选取的纵向模态对液体流 量变化不太敏感且频散严重，影响了超声导波技术对液体流量测量的准确性 和超声导波技术对液体流量测量的应用潜力，并且没有考虑超声导波纵向模 态的激励方式，文中采用的非轴对称激励方式无法有效地产生所选取的轴对 称纵向模态。
目前，对用于管道液体流量测量的超声导波纵向模态及其激励方式的优 化选取迄今国内外尚未见相关报道。
发明内容
本发明提出了一种基于超声导波技术的管道液体流量测量方法，该方法 可以对管道尤其是小管径管道中液体流量的大小进行快速、有效地测量。
该方法基于对超声导波纵向模态的理论分析，选取激励频率处的群速度 随频率变化率的绝对值低于0.002m，轴向位移在管中液体分布的平均值为在 管壁中分布的平均值的30%以上，且群速度随管中液体流速呈单调增加或减小的超声导波纵向模态用于管道液体流量的测量。
本发明所采用的装置包括充液管道、两个分别用于激励和接收超声导 波纵向模态的轴对称分布的传感器组(环)、功率放大器、转换装置、函数发 生器、示波器和计算机等。其中一个传感器组(环)安装于充液管道某一位 置，另一传感器组(环)安装于充液管道另一位置。两个传感器组(环)与 转换装置的两个端口连接，功率放大器、示波器分别与转换装置的另外两个 端口连接，计算机和示波器连接，函数发生器分别与功率放大器、示波器连 接。
本发明的基于超声导波技术的管道液体流量测量方法包括以下步骤-
1) 在管道l的位置A处安装轴对称分布的第一传感器组(环)2，在管 道1的另一位置B处安装轴对称分布的第二传感器组(环)3，设位置A和位 置B的轴向间距为丄；
2) 由函数发生器5产生一个具有一定中心频率的单音频信号，该中心频 率为所选取超声导波纵向模态的频率点；在该频率点处，超声导波纵向模态 的群速度随频率变化率的绝对值低于0.002m,轴向位移在管中液体分布的平 均值为在管壁中分布的平均值的30%以上，且群速度随管中液体流速呈单调 增加或减小；
3) 将转换装置8的C端和E端连接，D端和F端连接，即第一传感器组 (环)2与功率放大器4连接，用来激励超声导波纵向模态，将第二传感器组 (环)3与示波器6连接，用来接收超声导波纵向模态的信号；
4) 函数发生器5产生一激励信号，该激励信号经功率放大器4进行功率 放大；通过第一传感器组(环)2，在管道1中激励所选取的超声导波纵向模态；
5) 激励的超声导波纵向模态信号在管道1中从第一传感器组(环)2传 播至第二传感器组(环)3,即在管道中传播距离丄后，第二传感器组(环) 3接收信号a，在示波器6上显示，并通过数据端口存储到计算机7中；
6) 将转换装置8的C端和F端连接，D端和E端连接，即将第二传感器 组(环3与功率放大器4连接，用来激励超声导波纵向模态，将第一传感器 组(环)2与示波器6连接，用来接收超声导波纵向模态的信号；
7) 由函数发生器5产生与步骤4)相同的激励信号，经功率放大器4进 行功率放大；通过第二传感器组(环)3，在管道l中激励与步骤4)相同的 超声导波纵向模态；
8) 再次激励的超声导波纵向模态信号在管道1中从第二传感器组(环)3 传播至第一传感器组(环)2，即在管道中传播了距离丄后，第一传感器组(环) 2接收信号b,在示波器6上显示，并通过数据端口存储到计算机7中；
9) 确定所得到的超声导波纵向模态信号a和信号b传播距离丄所需的时 间分别为A和/2 ，从而得到超声导波纵向模态的传播时间差为= -/2;
10) 将管道中的液体流速v代入考虑液体流速的充液管道中超声导波纵 向模态的频散方程，通过数值求解得到不同液体流速v时的相速度cv频散曲
线，进而利用^=《/
l" 众
c。 一必—11 、p cto 乂
,其中必为圆频率，即得到超声导波纵向
模态的群速度^频散曲线，通过分析所激励的超声导波纵向模态在不同液体
流速v下的群速度频散曲线，可以得到某一频率点的一组液体流速v-超声导波 纵向模态群速度C。的数据点，通过线性拟合，得到超声导波纵向模态的群速度^随液体流速v的变化率5^」和液体流速为0即液体静止时该频率点的
Av
超声导波纵向模态的群速度^。；
11) 利用式p^^，即可计算得到管道中液体流速v。当计算得到的液 体流速v大于零时，表示液体流向为从A流向B，当液体流速v小于零时，表
示液体流向为从B流向A;
12) 利用式[7-7rw;2，即可求得充液管道中液体的体积流量C7，进而可以
确定充液管道中液体的质量流量# = /7" = 7t/t^2。其中，G为管道的内半径，
p为管道中液体密度。
与现有管道液体流量的测量方法相比，本发明具有以下优点1)可以对 管道尤其是小管径管道中液体流量的大小进行快速、有效地测量；2)不需改 变管道的结构，只需将传感器安装在管道外壁即可实行在线、无损测量；3)可 充分利用超声导波的多模态特性，优化选取群速度对液体流速敏感的超声导 波纵向模态用于管道中液体流量大小的测量，检测方便，测量可靠性高，效 率高，劳动强度低。
图1基于超声导波技术的管道液体流量测量装置原理图2水的流速v为0即水静止时充水钢管中超声导波纵向模态的群速度频散曲
线；
图3水的流速v不同时充水钢管中超声导波最低阶纵向模态^的群速度频散
曲线；其中，水的流速v变化范围为-100 100m/s，步长为10m/s;
图4充水钢管中频率50kHz的超声导波最低阶纵向模态a,的群速度与水的流速V之间的变化关系图。
图中1、充液管道，2、传感器组(环)，3、传感器组(环)，4、功率 放大器，5、函数发生器，6、示波器，7、计算机，8、转换装置，9、水的流 速v为0即水静止时充水钢管中超声导波最低阶纵向模态^的群速度频散曲 线，10、充水钢管中频率50kHz的超声导波最低阶纵向模态^的群速度与水 的流速v之间的线性拟合关系直线。
具体实施例方式
结合本发明方法的内容提供以下数值仿真得到的实施例
(1) 将沿充液管道的轴线对称分布的第一传感器组(环)2和第二传感组
(环)3安装在同一充水钢管i中，两者相距i:。本实施例中充水钢管的内半
径G = 4mm和外半径^ = 5mm ;
(2) 理论分析得到考虑管中水的流速v的充水钢管中超声导波纵向模态的
相速度c。频散方程lDh
化
A2A5
A,A2A5
A,A2
A,A4D55
=0，通过数值求解得到不同
水的流速v时充水钢管中超声导波纵向模态的相速度^频散曲线，进而利用
—，得到不同水的流速v时充水钢管中超声导波纵向模态的群速
血
c 一必'
d必
度c。频散曲线;
其中，D =//
L、Cp
^( r2)+2^^—2) 八D12 =//
Q 、
■、cp 乂 2辯
，"13 =
2，
p 匕
风,
> 一 v/ > —
一 2"《(or2)
"22 — , "23 一
Q 、
、cp 乂
V 、
一,
、cp 乂
2)，
"25 一 V ， "31 一 ， "32 — ，
、cp J 、cp J
Al,'W， KW， D43=^^， A4=^^，
/2 ^
、cp 乂
A4 =
2"
//02(^)+2-7^(c^)
，A3 =
2，
必2 <y2
风W-
AW
，A5 =々
必
1
々2=^_^，《)(z)和《)(z)为汉克尔函数(z为任意值)，
c， c.
且
(v = 0，l)， J。(z)和A(》分别为0阶和1阶第一
类贝塞尔函数，^z)和l^;)分别为O阶和1阶第二类贝塞尔函数，Cl =5960m/s 和q =3260111/8分别为钢的纵波波速和横波波速，// = /^22为钢的Lame常数，
/^7932kg/m3为钢的密度，^1^， ^:1500m/s为水的纵波波
c.,,lc_ 一 V)
11速，^-2为水的Lame常数，/ w = 1000kg/m3为水的密度；
(3) 图2给出了频带0-2.0MHz范围内，水的流速v为O即水静止时充水 钢管中超声导波纵向模态的群速度频散曲线，选取频率50kHz的超声导波最 低阶纵向模态^用于钢管中水的流量v测量，在该频率点处，纵向模态^的 群速度随频率变化率的绝对值低于0.002m，轴向位移在管中液体分布的平均 值为在管壁中分布的平均值的30%以上，且群速度随管中液体流速呈单调增 加或减小；
(4) 图3给出了充水钢管中超声导波最低阶纵向模态^在不同水的流速v
下的群速度频散曲线，频带为0-0.5MHz。图4给出了充水钢管中频率50kHz的
超声导波最低阶纵向模态"L的群速度与水的流速v之间的变化关系图。通过 线性拟合，可以确定频率50kHz的纵向模态^的群速度^随水的流速v的变化
率S = f = 0.995(m/sy(m/s)和水的流速v为O即水静止时该频率点的超声导波
△V
纵向模态的群速度^。 = 1413.08m/s ;
(5) 由函数发生器5产生一个具有一定中心频率的单音频信号。在本实施 中，该单音频信号的中心频率为50kHz。将转换装置8的C端和E端连接，D端 和F端连接，即第一传感器组(环)2与功率放大器4连接，用来激励超声导波 纵向模态，将第二传感器组(环)3与示波器6连接，用来接收超声导波纵向 模态的信号；
(6) 函数发生器5产生的单音频信号经功率放大器4放大，由第一传感器组 (环)2激励中心频率50kHz的超声导波纵向模态^。在该频率点处，纵向模
态"L的群速度随频率变化率的绝对值低于0.002m,轴向位移在管中液体分布的平均值为在管壁中分布的平均值的30%以上，且群速度随管中液体流速呈单 调增加或减小。经传播距离丄后，第二传感器组(环)3接收信号4，在示波 器6显示，并通过以太网端口存储到计算机7;
(7) 将转换装置8的C端和F端连接，D端和E端连接，即将第二传感 器组(环)3与功率传感器4连接，用来激励超声导波纵向模态，将第一传感 器组(环)2与示波器6连接，用来接收超声导波纵向模态的信号。由函数发 生器5产生同一激励信号，即中心频率50kHz的单音频信号。经功率放大器 4放大，由第二传感器组(环)3激励相同的中心频率50kHz的超声导波纵向 模态^，经传播距离丄后，第一传感器组(环)2接收信号4，在示波器6 显示，并通过以太网端口存储到计算机7;
(8) 数据分析得到接收的超声导波最低阶纵向模态a,的信号羊和4在传 播同一距离丄的时间^和/2，确定两者之间的传播时间差为Af-^-^
(9) 可以确定充水钢管中水的流速"^^-1.0034xl06x^(m/s);
(10) 钢管中水的体积流量￡/ =卩乂(4乂10-3)271==50.436><玍(1113/3),水的质
丄
量流量M = .= 50436x^(kg/s)，从而实现管道中液体流量的测量。1权利要求
1、基于超声导波技术的管道液体流量测量方法，其特征在于该方法包括以下步骤1)在管道(1)的位置A处安装轴对称分布的第一传感器组(2)，在管道(1)的另一位置B处安装轴对称分布的第二传感器组(3)，设位置A和位置B的轴向间距为L；2)由函数发生器(5)产生一个具有一定中心频率的单音频信号，该中心频率为所选取超声导波纵向模态的频率点，在该频率点处，超声导波纵向模态的群速度随频率变化率的绝对值低于0.002m，轴向位移在管中液体分布的平均值为在管壁中分布的平均值的30％以上，且群速度随管中液体流速呈单调增加或减小；3)调整转换装置，使第一传感器组(2)与功率放大器(4)连接，用来激励超声导波纵向模态，将第二传感器组(3)与示波器(6)连接，用来接收超声导波纵向模态的信号；4)函数发生器(5)产生一激励信号，该激励信号经功率放大器(4)进行功率放大，通过第一传感器组(2)在管道(1)中激励步骤2)所选取的超声导波纵向模态；5)激励的超声导波纵向模态信号在管道(1)中从第一传感器组(2)传播至第二传感器组(3)，即在管道中传播距离L后，第二传感器组(3)接收信号a，在示波器(6)上显示，并通过数据端口存储到计算机(7)中；6)调整转换装置，使第二传感器组(3)与功率放大器(4)连接，用来激励超声导波纵向模态，将第一传感器组(2)与示波器(6)连接，用来接收超声导波纵向模态的信号；7)由函数发生器(5)产生与步骤4)中相同的激励信号，该激励信号经功率放大器(4)进行功率放大，通过第二传感器组(3)在管道(1)中激励与步骤4)相同的超声导波纵向模态；8)再次激励的超声导波纵向模态信号在管道(1)中从第二传感器组(3)传播至第一传感器组(2)，即在管道中传播了距离L后，第一传感器组(2)接收信号b，在示波器(6)上显示，并通过数据端口存储到计算机(7)中；9)确定所得到的超声导波纵向模态信号a和信号b传播距离L所需的时间分别为t1和t2，从而得到超声导波纵向模态的传播时间差为Δt＝t1-t2；10)将管道中的液体流速v代入考虑液体流速的充液管道中超声导波纵向模态的频散方程，通过数值求解得到不同液体流速v时的相速度cp频散曲线，进而利用公式<maths id="math0001" num="0001" ><math><![CDATA[ <mrow><msub> <mi>c</mi> <mi>g</mi></msub><mo>=</mo><msubsup> <mi>c</mi> <mi>p</mi> <mn>2</mn></msubsup><mo>/</mo><mrow> <mo>(</mo> <msub><mi>c</mi><mi>p</mi> </msub> <mo>-</mo> <mi>ω</mi> <mo>·</mo> <mfrac><msub> <mi>dc</mi> <mi>p</mi></msub><mi>dω</mi> </mfrac> <mo>)</mo></mrow><mo>,</mo> </mrow>]]></math> id="icf0001" file="A2009100888420003C1.tif" wi="45" he="13" top= "137" left = "61" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>其中ω为圆频率，即得到超声导波纵向模态的群速度cg频散曲线，通过分析所激励的超声导波纵向模态在不同液体流速v下的群速度频散曲线，得到某一频率点的一组该超声导波纵向模态的液体流速v和超声导波纵向模态群速度cg的数据点，通过线性拟合，得到超声导波纵向模态的群速度cg随液体流速v的变化率<maths id="math0002" num="0002" ><math><![CDATA[ <mrow><mi>S</mi><mo>=</mo><mfrac> <mrow><mi>Δ</mi><msub> <mi>c</mi> <mi>g</mi></msub> </mrow> <mi>Δv</mi></mfrac> </mrow>]]></math> id="icf0002" file="A2009100888420003C2.tif" wi="16" he="11" top= "189" left = "136" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>和液体流速为0即液体静止时该频率点的超声导波纵向模态的群速度cg0；11)利用公式<maths id="math0003" num="0003" ><math><![CDATA[ <mrow><mi>v</mi><mo>=</mo><mfrac> <mrow><mi>Δt</mi><mo>·</mo><msubsup> <mi>c</mi> <mrow><mi>g</mi><mn>0</mn> </mrow> <mn>2</mn></msubsup> </mrow> <mrow><mn>2</mn><mi>S</mi><mo>·</mo><mi>L</mi> </mrow></mfrac> </mrow>]]></math> id="icf0003" file="A2009100888420003C3.tif" wi="20" he="11" top= "217" left = "60" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>计算得到管道中液体流速v；当计算得到的液体流速v大于零时，表示液体流向为从A流向B，当液体流速v小于零时，表示液体流向为从B流向A；12)利用式U＝πvr12，求得充液管道中液体的体积流量U，进而确定充液管道中液体的质量流量M＝ρU＝πρvr12，其中，r1为管道的内半径，ρ为管道中液体密度。
2、根据权利要求1所述的基于超声导波技术的管道液体流量测量方法，其特 征在于所述的第一传感器组和第二传感器组能够由传感器环代替。
全文摘要
本发明涉及基于超声导波技术的管道液体流量测量方法，属于测试计量及无损检测技术领域。本发明选取激励频率处的群速度随频率变化率的绝对值低于0.002m，轴向位移在管中液体分布的平均值为在管壁中分布的平均值的30％以上，且群速度随管中液体流速呈单调增加或减小的超声导波纵向模态用于管道液体流量的测量。本发明可以对管道尤其是小管径管道中液体流量的大小进行快速、有效地测量。
文档编号G01F1/66GK101614569SQ20091008884
公开日2009年12月30日 申请日期2009年7月20日 优先权日2009年7月20日
发明者何存富, 刘增华, 斌 吴, 张力新, 樊尚春 申请人:北京工业大学