CN101922954A - 超声波测量信号的包络线处理方法 - Google Patents
超声波测量信号的包络线处理方法 Download PDFInfo
- Publication number
- CN101922954A CN101922954A CN2009100532642A CN200910053264A CN101922954A CN 101922954 A CN101922954 A CN 101922954A CN 2009100532642 A CN2009100532642 A CN 2009100532642A CN 200910053264 A CN200910053264 A CN 200910053264A CN 101922954 A CN101922954 A CN 101922954A
- Authority
- CN
- China
- Prior art keywords
- ultrasonic
- probe
- signal
- medium
- tau
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
Description
[技术领域]
本发明涉及气体或液体超声波流量计的一种超声波测量信号包络线处理方法。
[背景技术]
漩涡流量计、电磁流量计、转子流量计、孔板流量计和涡轮流量计等传统流量计在液体计量和气体计量市场上占着主导地位,显示了他们在计量方面的优越性。但是,它们自身存在着无法克服的不足,比如说测量管段中存在可动部件、量程比窄、对测量介质要求过高等等。面对不断提高的计量要求,超声波流量计应运而生。
超声波流量计发展短短不过几十年,特别是90年代以后由于电子和微处理技术的飞速发展为超声波流量计的发展提供了更高的平台,超声波流量计开始在各个行业计量管理中得到广泛的应用,相对与传统流量计的优势也已渐渐显现出来,比如特别适合大口径测量、量程比特别宽、无可动部件、拆卸方便故障率低、对测量介质要求少等等。
但是现在的超声波流量计也存在技术难题,例如:不同介质的声阻抗不同,在两种介质的交接面上透过率会有不同,特别在气体超声波流量计中换能器和气体的声阻抗相差很大,声波在两者交接面上大部分被反射,透过的只是一小部分,再加上气体对声波的衰减比较利害,接收探头接收到的回波已经很微弱了。
困难不仅仅在超声波信号的发射与接收上面,在接收到的超声波测量信号处理上面也有很多困难,测量信号处理的好坏直接影响到流量计的精度,目前应用的比较多的是多普勒法、互相关法等测量信号处理方法,分辨率及精度不高。
[发明内容]
本发明提出一种超声波测量信号包络线处理方法,用以提高超声波流量计特别是气体超声波流量计的分辨率及精度。
本发明的方法按以下步骤进行:
1)在测量管道内轴线上安装一对超声波测量超声A探头和超声B探头,测得超声波信号,设由A探头发射超声波至B探头接收到超声波所经历的传播时间为t1,由B探头发射超声波至A探头接收到超声波所经历的传播时间为t2;
2)测量得到管段直径D、超声探头与轴线夹角θ和由电路系统及机械系统共同带来的延时τ的数值;
3)对所述超声波信号进行放大、滤波、整形,得到整形信号;
4)对所述整形信号进行高速采样,再经过软件归一化处理,得到信号、幅值基本一致的处理信号;
5)用最小二乘法对所述处理信号的幅值拟合出一条顶点包络线;
6)用该包络线与设定的门限电压比较,得到所述的传播时间t1和t2,
7)将上述测得的数据代入下列公式,计算出介质流速v和介质的声速c:
介质的流速为
介质的声速为
实验证明本发明所述方法能较大程度提高超声波信号的分辨率,提高测量精度,切实有效可行。
[附图说明]
图1为本发明测量超声探头位置示意图。
图2为测得超声波信号示意图。
图3为整形后的处理信号示意图。
图4为整形后处理信号与门限电压比较示意图。
[具体实施方式]
下面结合和附图对本发明作进一步详细说明。
如图1所示,安装在轴线上的一对超声波测量探头被用来测量管道内的流速。
设由A探头发射声波至B探头接收到声波所经历的时间为t1;由B探头发射声波至A探头接收到声波所经历的时间为t2;D为测量管段直径,τ为电路系统及机械系统共同带来的延时,θ为超声探头与轴线的夹角,D、τ、θ的数值可由常用仪器测量得到;v为介质流速,c为介质中的声速。
测量管道内的声程为:
当测量管道中待测介质以速度v向前流动时超声波脉冲在介质中的顺、逆流传播时间t1和t2:
由上式求得重复频率分别为:
重复频率相减得:
这样求出管道中介质流速v:
另外,还可以精确地测出待测介质的声速c:
∴
在一般超声波流量计中,通常是采用固定门限电压加自动增益的方式来确定t1和t2,如图2所示。
由于测量过程中信号幅度的随流速不断的变化以及信号幅度变化的不确定性,例如,假定门限电压设定在Adj处,实际获得可能是Adi处或者是Adk处,这两处都会给测量带来比较大的误差。
如何消除因为信号幅度变化带来的误差就是本发明所具有的特点。具体实施步骤是:
对接收到的超声波信号进行放大、滤波、整形,得到如图3所示的整形信号,再由微处理器对整形信号进行高速采样,再经过软件归一化处理,得到信号、幅值基本一致的处理信号,如图3所示。最后用最小二乘法对处理信号拟合出一条顶点包络线,如图4所示。
拟合过程如下:
设{(xk,yk)}N K=1有N个点,横坐标是确定的.其对应的抛物线表示:y=f(x)=Ax2+Bx+C。
由最小二乘法的基本原理可知:对于给定的一组数据点坐标的分布,数据(xi,yi)(i=1,2,...m)可构造一个n(n<m)次的多项式:y=f(x)=Ax2+Bx+C,使得
得到最小值。
求解A,B,C的线性方程组为:
令
则方程组可改为:
S2A+S1B+NC=Sy---(3)
分别将(1)、(2)与(3)通分变换可得:
将A、B分别带入方程组任何一个式中均可求出C。
从而我们可以根据采样数据的时间与峰值的位置及幅度得到顶点包络线,如图4所示。
再用软件设定的门限电压与得到的包络线进行比较,可分别得到传播时间t1和t2,代入前述公式:
介质的流速为
介质的声速为
计算出介质流速v和介质的声速c。
Claims (1)
1.一种超声波测量信号包络线处理方法,按以下步骤进行:
1)在测量管道内轴线上安装一对超声波测量超声A探头和超声B探头,测得超声波信号,设由A探头发射超声波至B探头接收到超声波所经历的传播时间为t1,由B探头发射超声波至A探头接收到超声波所经历的传播时间为t2;
2)测量得到管段直径D、超声探头与轴线夹角θ和由电路系统及机械系统共同带来的延时τ的数值;
3)对所述超声波信号进行放大、滤波、整形,得到整形信号;
4)对所述整形信号进行高速采样,再经过软件归一化处理,得到信号、幅值基本一致的处理信号;
5)用最小二乘法对所述处理信号的幅值拟合出一条顶点包络线;
6)用该包络线与设定的门限电压比较,得到所述的传播时间t1和t2,
7)将上述测得的数据代入下列公式,计算出介质流速v和介质的声速c:
介质的流速为
介质的声速为
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100532642A CN101922954A (zh) | 2009-06-17 | 2009-06-17 | 超声波测量信号的包络线处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100532642A CN101922954A (zh) | 2009-06-17 | 2009-06-17 | 超声波测量信号的包络线处理方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101922954A true CN101922954A (zh) | 2010-12-22 |
Family
ID=43338020
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009100532642A Pending CN101922954A (zh) | 2009-06-17 | 2009-06-17 | 超声波测量信号的包络线处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101922954A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102967334A (zh) * | 2012-09-26 | 2013-03-13 | 朱作行 | 利用对信号包络线处理测量流体流量的系统及方法 |
CN103063868A (zh) * | 2013-01-05 | 2013-04-24 | 浙江大学 | 流体流速测量装置及方法 |
US20140303910A1 (en) * | 2013-04-04 | 2014-10-09 | Texas Instruments Incorporated | Extended range adc flow meter |
WO2015191785A1 (en) * | 2014-06-10 | 2015-12-17 | Texas Instruments Incorporated | Extended range adc flow meter |
CN106768109A (zh) * | 2017-02-21 | 2017-05-31 | 合肥工业大学 | 基于回波上升段峰值拟合和基于回波能量点定位的气体超声波流量计信号处理方法 |
CN107843305A (zh) * | 2017-10-31 | 2018-03-27 | 合肥工业大学 | 一类基于回波信号包络拟合的气体超声波流量计信号处理方法 |
CN113418573A (zh) * | 2021-06-18 | 2021-09-21 | 重庆市山城燃气设备有限公司 | 层流式气体测量流道、流量计、流量测量方法和系统 |
-
2009
- 2009-06-17 CN CN2009100532642A patent/CN101922954A/zh active Pending
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102967334A (zh) * | 2012-09-26 | 2013-03-13 | 朱作行 | 利用对信号包络线处理测量流体流量的系统及方法 |
CN103063868A (zh) * | 2013-01-05 | 2013-04-24 | 浙江大学 | 流体流速测量装置及方法 |
CN103063868B (zh) * | 2013-01-05 | 2014-08-20 | 浙江大学 | 流体流速测量装置及方法 |
US10955273B2 (en) * | 2013-04-04 | 2021-03-23 | Texas Instruments Incorporated | Extended range ADC flow meter |
US20140303910A1 (en) * | 2013-04-04 | 2014-10-09 | Texas Instruments Incorporated | Extended range adc flow meter |
WO2015191785A1 (en) * | 2014-06-10 | 2015-12-17 | Texas Instruments Incorporated | Extended range adc flow meter |
EP3155382A4 (en) * | 2014-06-10 | 2018-02-28 | Texas Instruments Incorporated | Extended range adc flow meter |
CN106461437B (zh) * | 2014-06-10 | 2020-07-17 | 德克萨斯仪器股份有限公司 | 扩展范围adc流量计 |
US10801868B2 (en) | 2014-06-10 | 2020-10-13 | Texas Instruments Incorporated | Extended range ADC flow meter |
CN106461437A (zh) * | 2014-06-10 | 2017-02-22 | 德克萨斯仪器股份有限公司 | 扩展范围adc流量计 |
US11255708B2 (en) | 2014-06-10 | 2022-02-22 | Texas Instruments Incorporated | Extended range ADC flow meter |
US11747181B2 (en) | 2014-06-10 | 2023-09-05 | Texas Instruments Incorporated | Extended range ADC flow meter |
CN106768109A (zh) * | 2017-02-21 | 2017-05-31 | 合肥工业大学 | 基于回波上升段峰值拟合和基于回波能量点定位的气体超声波流量计信号处理方法 |
CN107843305A (zh) * | 2017-10-31 | 2018-03-27 | 合肥工业大学 | 一类基于回波信号包络拟合的气体超声波流量计信号处理方法 |
CN107843305B (zh) * | 2017-10-31 | 2019-11-12 | 合肥工业大学 | 一种基于回波信号包络拟合的气体超声波流量计信号处理方法 |
CN113418573A (zh) * | 2021-06-18 | 2021-09-21 | 重庆市山城燃气设备有限公司 | 层流式气体测量流道、流量计、流量测量方法和系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101922954A (zh) | 超声波测量信号的包络线处理方法 | |
US6067861A (en) | Method and apparatus for ultrasonic doppler velocimetry using speed of sound and reflection mode pulsed wideband doppler | |
Chen et al. | Realization of a multipath ultrasonic gas flowmeter based on transit-time technique | |
CN106768103B (zh) | 一种超声波流量计自动校准时间偏差的方法 | |
CN101004353A (zh) | 一种用于超声波流量计时差交叉检测方法 | |
CN102261937A (zh) | 一种高精度时差式超声波流量计及其流量测量方法 | |
EP3155382B1 (en) | Extended range adc flow meter | |
CN105300508A (zh) | 一种超声波信号传播时间的双阈值检测电路及检测方法 | |
CN111157065A (zh) | 气体超声流量计超声波信号传输回路中声延时测量方法 | |
JPS5824816A (ja) | ドツプラ−型超音波流量計 | |
CN114088151B (zh) | 外夹式多声道超声波流量检测装置及检测方法 | |
CN100380101C (zh) | 多普勒型超声波流量计 | |
CN102866261B (zh) | 检测超声波在测流速中飞行时间的方法 | |
CN103162752A (zh) | 用于超声波流量计的相位编码同步时差检测装置及方法 | |
CN104457871A (zh) | 一种流量计及流体测量方法 | |
CN102914333B (zh) | 利用超声波检测流量的检测方法 | |
US20030125902A1 (en) | Peak switch detector for transit time ultrasonic meters | |
Han et al. | Studies on the transducers of clamp-on transit-time ultrasonic flow meter | |
CN102980539A (zh) | 锅炉受热面管管壁金属层和氧化层厚度的测量方法 | |
CN102829829A (zh) | 一种时差法超声波流量检测方法及装置 | |
CN102095889B (zh) | 三通道超声时差流速测量方法 | |
CN203069223U (zh) | 用于超声波流量计的相位编码同步时差检测装置 | |
CN102829830A (zh) | 用于超声波流量检测中检测超声波传播速度的方法及装置 | |
CN112903043A (zh) | 一种多声道超声波流量计系统 | |
KR101630301B1 (ko) | 자기왜곡방식의 거리측정을 이용한 테이퍼관형 면적식 유량계 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20101222 |