专利名称：超声波测量信号的包络线处理方法
超声波测量信号的包络线处理方法本发明涉及气体或液体超声波流量计的一种超声波测量信号包络线处理方法。 [背景技术]漩涡流量计、电磁流量计、转子流量计、孔板流量计和涡轮流量计等传统流量计在 液体计量和气体计量市场上占着主导地位，显示了他们在计量方面的优越性。但是，它们自 身存在着无法克服的不足，比如说测量管段中存在可动部件、量程比窄、对测量介质要求过 高等等。面对不断提高的计量要求，超声波流量计应运而生。超声波流量计发展短短不过几十年，特别是90年代以后由于电子和微处理技术 的飞速发展为超声波流量计的发展提供了更高的平台，超声波流量计开始在各个行业计量 管理中得到广泛的应用，相对与传统流量计的优势也已渐渐显现出来，比如特别适合大口 径测量、量程比特别宽、无可动部件、拆卸方便故障率低、对测量介质要求少等等。但是现在的超声波流量计也存在技术难题，例如不同介质的声阻抗不同，在两种 介质的交接面上透过率会有不同，特别在气体超声波流量计中换能器和气体的声阻抗相差 很大，声波在两者交接面上大部分被反射，透过的只是一小部分，再加上气体对声波的衰减 比较利害，接收探头接收到的回波已经很微弱了。困难不仅仅在超声波信号的发射与接收上面，在接收到的超声波测量信号处理上 面也有很多困难，测量信号处理的好坏直接影响到流量计的精度，目前应用的比较多的是 多普勒法、互相关法等测量信号处理方法，分辨率及精度不高。本发明提出一种超声波测量信号包络线处理方法，用以提高超声波流量计特别是 气体超声波流量计的分辨率及精度。本发明的方法按以下步骤进行1)在测量管道内轴线上安装一对超声波测量超声A探头和超声B探头，测得超声 波信号，设由A探头发射超声波至B探头接收到超声波所经历的传播时间为tl，由B探头发 射超声波至A探头接收到超声波所经历的传播时间为t2 ；2)测量得到管段直径D、超声探头与轴线夹角θ和由电路系统及机械系统共同带 来的延时τ的数值；3)对所述超声波信号进行放大、滤波、整形，得到整形信号；4)对所述整形信号进行高速采样，再经过软件归一化处理，得到信号、幅值基本一 致的处理信号；5)用最小二乘法对所述处理信号的幅值拟合出一条顶点包络线；6)用该包络线与设定的门限电压比较，得到所述的传播时间tl和t2，7)将上述测得的数据代入下列公式，计算出介质流速ν和介质的声速c
介质的流速为ν 介质的声速为 实验证明本发明所述方法能较大程度提高超声波信号的分辨率，提高测量精度， 切实有效可行。
图1为本发明测量超声探头位置示意图。图2为测得超声波信号示意图。图3为整形后的处理信号示意图。图4为整形后处理信号与门限电压比较示意图。下面结合和附图对本发明作进一步详细说明。如图1所示，安装在轴线上的一对超声波测量探头被用来测量管道内的流速。设由A探头发射声波至B探头接收到声波所经历的时间为、；由B探头发射声波 至A探头接收到声波所经历的时间为t2;D为测量管段直径，τ为电路系统及机械系统共 同带来的延时，θ为超声探头与轴线的夹角，D、τ、θ的数值可由常用仪器测量得到；ν为 介质流速，c为介质中的声速。
2D测量管道内的声程为 当测量管道中待测介质以速度ν向前流动时超声波脉冲在介质中的顺、逆流传播 时间、和t2 由上式求得重复频率分别为 重复频率相减得 这样求出管道中介质流速ν 另外，还可以精确地测出待测介质的声速c 在一般超声波流量计中，通常是采用固定门限电压加自动增益的方式来确定、和 t2，如图2所示。由于测量过程中信号幅度的随流速不断的变化以及信号幅度变化的不确定性，例 如，假定门限电压设定在α+处，实际获得可能是adi处或者是adk处，这两处都会给测量带 来比较大的误差。如何消除因为信号幅度变化带来的误差就是本发明所具有的特点。具体实施步骤 是对接收到的超声波信号进行放大、滤波、整形，得到如图3所示的整形信号，再由 微处理器对整形信号进行高速采样，再经过软件归一化处理，得到信号、幅值基本一致的处 理信号，如图3所示。最后用最小二乘法对处理信号拟合出一条顶点包络线，如图4所示。拟合过程如下设{(Xk，yk)}nκ = 1有n个点，横坐标是确定的.其对应的抛物线表示y = f(x)= Ax2+bx+c。由最小二乘法的基本原理可知对于给定的一组数据点坐标的分布，数据(Xi，yi) (i = 1，2，. . .m)可构造一个 η (n < m)次的多项式y = f(x) = ax2+bx+c，使得 得到最小值。令偏导数况/拟，况/湖,况为零，可得ο = dE{A^C) = 2± (Αχ] + BXk + C 求解a，b, c的线性方程组为 令
JV 4NN 2NN2NN
S,=Σχ1'53=ΣχρΞ2=ΣχρΞι=ΣχρΞ/=Σ);^Ιι^γ,=Σ)^ JlXeSy=Tyk
k=\Jt=Ik=lk=\k=\Kk=lk=l则方程组可改为S4A + S3B + S2C = Syl---(1)S,A + S2B + SxC^Syl— (2)S^+S^+NC = Sy—(3)分别将⑴、⑵与(3)通分变换可得j _ [ (^4-S22XNS2-S^-jNS,-S2S1)21 X 1 ^ 1 3 1 2)\NSyl-S2S ^(NS2-S12)-(NS,-S1S^(NS3-S1S2)
~NS2-S12将A、B分别带入方程组任何一个式中均可求出C。从而我们可以根据采样数据的时间与峰值的位置及幅度得到顶点包络线，如图4 所示。再用软件设定的门限电压与得到的包络线进行比较，可分别得到传播时间、和 t2,代入前述公式
D ( 1 1 )介质的流速为V =----
sm6'cos6' ^i1 - τ t2-r )
D ( 1 1 ^介质的声速为C = ^7 -+-
Sinef^i1 -τ t2-r J计算出介质流速ν和介质的声速C。
权利要求
一种超声波测量信号包络线处理方法，按以下步骤进行1)在测量管道内轴线上安装一对超声波测量超声A探头和超声B探头，测得超声波信号，设由A探头发射超声波至B探头接收到超声波所经历的传播时间为t1，由B探头发射超声波至A探头接收到超声波所经历的传播时间为t2；2)测量得到管段直径D、超声探头与轴线夹角θ和由电路系统及机械系统共同带来的延时τ的数值；3)对所述超声波信号进行放大、滤波、整形，得到整形信号；4)对所述整形信号进行高速采样，再经过软件归一化处理，得到信号、幅值基本一致的处理信号；5)用最小二乘法对所述处理信号的幅值拟合出一条顶点包络线；6)用该包络线与设定的门限电压比较，得到所述的传播时间t1和t2，7)将上述测得的数据代入下列公式，计算出介质流速v和介质的声速c介质的流速为 <mrow><mi>v</mi><mo>=</mo><mfrac> <mi>D</mi> <mrow><mi>sin</mi><mi></mi><mi>θ</mi><mi>cos</mi><mi>θ</mi> </mrow></mfrac><mrow> <mo>(</mo> <mfrac><mn>1</mn><mrow> <msub><mi>t</mi><mn>1</mn> </msub> <mo>-</mo> <mi>τ</mi></mrow> </mfrac> <mo>-</mo> <mfrac><mn>1</mn><mrow> <msub><mi>t</mi><mn>2</mn> </msub> <mo>-</mo> <mi>τ</mi></mrow> </mfrac> <mo>)</mo></mrow><mo>;</mo> </mrow>介质的声速为 <mrow><mi>c</mi><mo>=</mo><mfrac> <mi>D</mi> <mrow><mi>sin</mi><mi>θ</mi> </mrow></mfrac><mrow> <mo>(</mo> <mfrac><mn>1</mn><mrow> <msub><mi>t</mi><mn>1</mn> </msub> <mo>-</mo> <mi>τ</mi></mrow> </mfrac> <mo>+</mo> <mfrac><mn>1</mn><mrow> <msub><mi>t</mi><mn>2</mn> </msub> <mo>-</mo> <mi>τ</mi></mrow> </mfrac> <mo>)</mo></mrow><mo>.</mo> </mrow>
全文摘要
一种超声波测量信号的包络线处理方法在测量管道内轴线上安装一对超声波测量超声A探头和超声B探头，测得超声波信号，设由A探头发射超声波至B探头接收到超声波所经历的传播时间为t1，由B探头发射超声波至A探头接收到超声波所经历的传播时间为t2；测量得到管段直径D、超声探头与轴线夹角θ和延时数值τ；对所述超声波信号进行放大、滤波、整形，得到整形信号；对所述整形信号进行高速采样，再经过软件归一化处理，得到信号、幅值基本一致的处理信号；拟合出一条顶点包络线与设定的门限电压比较，计算出介质流速v和介质的声速c介质的流速为介质的声速为实验证明本发明能较大程度提高超声波信号的分辨率和测量精度。
文档编号G01P5/24GK101922954SQ20091005326
公开日2010年12月22日 申请日期2009年6月17日 优先权日2009年6月17日
发明者傅金林, 叶荣辉, 吴叶军, 黄根涛 申请人:上海一诺仪表有限公司