CN102866261B - 检测超声波在测流速中飞行时间的方法 - Google Patents
检测超声波在测流速中飞行时间的方法 Download PDFInfo
- Publication number
- CN102866261B CN102866261B CN201210344506.5A CN201210344506A CN102866261B CN 102866261 B CN102866261 B CN 102866261B CN 201210344506 A CN201210344506 A CN 201210344506A CN 102866261 B CN102866261 B CN 102866261B
- Authority
- CN
- China
- Prior art keywords
- frequency
- ultrasonic
- phase
- time
- flight time
- 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.)
- Expired - Fee Related
Links
Abstract
本发明公开了一种检测超声波在测流速中飞行时间的方法,属于超声波仪器的检测技术领域,其技术方案采用了双频超声波相位检查的方法计算超声波飞行时间,有效增加了仪器测量的量程。并且不同于其他多频超声波测量方法,在对传感器存在相位漂移的情况给出了解决办法,使其在超声波流速或者流量的检查中更有实用性。
Description
技术领域
本发明涉及一种超声波仪器的检测方法,具体涉及一种超声波飞行时间的检测方法。
背景技术
超声波在媒质中传播时,其传播的时间会受到流体流动速度的影响。当声波和流体传播方向相同时,声波飞行时间会减小;当声波和流体传播反向时,声波的飞行时间会增加。声速,流体速度,和飞行时间的关系如下:
其中t为声波飞行时间,c为声速,v为流体流速,L为距离。只要测得声波飞行时间,则可求出流体流速。进而也可以计算流体流量。因此超声波飞行时间的检测是超声波风速仪,超声波流量计等仪器的重要测量物理量。
现有技术一(阈值比较法):
目前,超声波飞行时间的检测最广泛使用的方法是声波阈值检测法。其结构基本如图1所示。
处理器首先向超声波驱动电路发送方波,此时开始启动处理器内部的计数器进行计数。驱动电路将方波放大并驱动超声波发射器,向流体媒质中发射超声波。超声波在流体中飞行一段时间后进入接收传感器,接收到的波形如图2。当接收的波形经过硬件的电压阈值比较后,转换成方波。当处理器检测到方波的上升沿时,处理器内部的计数器停止计数。此时计数器所计的时间即为超声波传播的时间。
现有技术二(相位法):
除了技术一以外,目前在超声波流速测量中使用比较多的还有相差法。当媒介中的流体静止时,超声波探头接收到的波形是一条正弦曲线或者余弦曲线。当流体发生流动时,超声波的飞行时间变化,接收到的波形和原有波形发生一定相差。通过模拟或者数字鉴相器测量这个相位差来计算超声波的飞行时间,如图3所示。
现有技术一是目前超声波流速或者流量测量里最广泛的方法。但为了保证测量精度,就需要使用中心频率较高的高频超声波换能器,增加了产品成本。且由于超声波的接收波形幅值会由于流速的变化而改变,因此在电压阈值比较的时候可能会出现声波整周期的计算误差。
现有技术二在实际产品中也有一定应用。如果使用单频的相位差来计算飞行时间,则所测飞行时间的变化范围只能在一个超声波周期内,测量范围较小,很难满足实际应用的要求。也有使用多频法测量的方法,但都停留在理论阶段,很难用于实际应用。
发明内容
本发明克服了现有技术的不足,提供效果更好的检查超声波在测流速中飞行时间的方法。
考虑到现有技术的上述问题,根据本发明的一个方面,为解决上述的技术问题,本发明采用以下技术方案:
一种检测超声波在测流速中飞行时间的方法,包括以下步骤:
步骤1,处理器发送为频率f1的超声波信号,经过驱动电路驱动放大以后,再由超声波发射器发射出;同时,计数器开始计数;
步骤2,计数器的计数一段时间为t0后停止计数;
步骤3,,接收传感器将接收到的为频率f1的超声波信号经模拟电路放大和滤波处理后传送给模数转换器;
步骤4,模数转换器将转化后的超声波数字信号发送给处理器,在所述处理器中完成频率f1的相位差的计算;本步骤中,相位差,其中当超声波在媒质中传播一段时间后,接收到的波形相位差在媒质静止情况下为,当媒质为流体时,接收到的波形相位差为;
步骤5,由步骤1中所述的处理器发出频率f2的超声波信号,所述频率f2的超声波经过驱动电路驱动放大以后,再由超声波发射器发射出;同时,计数器开始计数;
步骤6,所述计数器的计数一段时间为t0后停止计数;
步骤7,在所述计数器停止计数后,接收传感器将接收到的为频率f2的超声波信号经模拟电路放大和滤波处理后传送给模数转换器;
步骤8,所述模数转换器再将转化后的超声波数字信号发送给处理器,在所述处理器中完成频率f2的相位差计算;本步骤中,相位差,其中其中当超声波在媒质中传播一段时间后,接收到的波形相位差在媒质静止情况下为,当媒质为流体时,接收到的波形相位差为;
步骤9,根据两个频率f1和f2的相频关系推导计算没有相位偏移量时的传播时间
其中k和k’为…-2,-1,0,1,2…
或在实际测量中相位会因为存在偏移误差,设此偏移误差全部产生在上且偏移量为,则传播时间为
或根据两个频率f1和f2的和频计算传播时间为,
或根据两个频率f1和f2的差频计算传播时间为,
在所述的式(1)、式(2)、式(3)、式(4)中,通过计算的大小判断所测量的与是否在同一周期内的情况,以确定k与k’的关系,计算出Δt,得出超声波的飞行时间为t=t0+Δt。
作为优选,所述处理器为FPGA。
作为优选,所述处理器为CPLD。
本发明还可以是:
作为优选,所述频率f1的相位计算采用DFT运算方法。
作为优选,所述频率f2的相位计算采用DFT运算方法。
与现有技术相比,本发明的有益效果之一是:
本发明的技术方案中,使用双频超声波相位检查的方法计算超声波飞行时间。此种方法可以有效增大测试的量程,并降低由于温度等原因造成的传感器相位漂移。同时,本方法有效增加了仪器测量的量程;并且不同于其他多频超声波测量方法,在对传感器存在相位漂移的情况给出了解决办法,使其在超声波流速或者流量的检查中更有实用性。
附图说明
为了更清楚的说明本申请文件实施例或现有技术中的技术方案,下面将对实施例或现有技术的描述中所需要使用的附图作简单的介绍,显而易见地,下面描述中的附图仅是对本申请文件中一些实施例的参考,对于本领域技术人员来讲,在不付出创造性劳动的情况下,还可以根据这些附图得到其它的附图。
图1、图2和图3为背景技术中的参考示意图;
图4示出了根据本发明一个实施例所使用的硬件结构框图。
具体实施方式
下面结合实施例对本发明作进一步地详细说明,但本发明的实施方式不限于此。
本方案中测量超声波飞行时间采用发射双频超声波的方法,方案的具体方法如下:
超声波发送换能器发射两个频率f1和f2的超声波。当超声波在媒质中传播一段时间后,接收到的波形在媒质静止情况下为与。当流体流动时候,接收的相位变为与。设两个频率的相位差为和,则由三角函数公式我们可以得出:
首先根据相位差的大小范围可以判断出k与k’的关系,当与在同一周期里时,则k=k’。在实际测量中相位会因为环境和器件等因素有一定的偏移误差,设此偏移误差全部产生在上且偏移量为,则上式改为
根据测量仪器的量程要求,可以确定k为几个固定值,如k=-1,0,1。当k取值不同时,对仪器的偏移误差的要求也不同。当k=-1,0,1时,只要满足远小于即可以解出k的唯一解。当我们求得k值以后,根据式2就可以求得飞行时间t。
当与不在同一周期里时,首先的大小可以判断出k与k’的关系,例如k’=k+1,再代入上式计算出k。
此方法不仅可以由公式1中使用两个单频f1和f2的相频关系推导,也可使用两个频率的和频或者差频来推导,例如使用公式或,推导过程和上面相同。也为本专利的保护范围。
图4示出了根据本发明一个实施例的结构框图。
由FPGA或者CPLD组成的处理器首先发送频率f1的信号,经过驱动电路驱动放大以后,经过超声波发射器发射出超声波。在发送信号的同时开启FPGA内部计数器开始计数。当计数器计数到一定数值时候,停止计数,并开始启动ADC采集信号。此时计数器计数时间为t0,且超声波在媒质中传播了一段时间并进入到接收传感器。接收到的信号经过模拟电路的放大,滤波处理后传送给模数转换器ADC。ADC再将转化的数字信号发送给FPGA,在FPGA中完成频率f1的相位计算。此计算相位的方法可以使用DFT运算,但不限于DFT方法。计算完f1的相位后再由FPGA发送频率f2的信号,与f1的操作过程相同,完成f2的相位计算。在流速或者流量计设计之初,我们首先会确定仪器测量流速的量程,根据此量程可确定k为几个固定整数值。在这个量程内,通过计算的大小判断所测量的与是否在同一周期内,即得到k与k’的关系。将所测量的与代入公式2中,并使用k与k’的约束条件,在FPGA中解出k值。根据上述式2计算出Δt,则最终的飞行时间即为t=t0+Δt。
本发明中使用双频超声波相位检查的方法计算超声波飞行时间。此种方法可以有效增大测试的量程,并降低由于温度等原因造成的传感器相位漂移。本发明的保护点即为此双频飞行时间测量方法,以及专利中提到同理的扩展方法。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其它实施例的不同之处,各个实施例之间相同相似部分相互参见即可。
在本说明书中所谈到的“一个实施例”、“另一个实施例”、 “实施例”、等,指的是结合该实施例描述的具体特征、结构或者特点包括在本申请概括性描述的至少一个实施例中。在说明书中多个地方出现同种表述不是一定指的是同一个实施例。进一步来说,结合任一实施例描述一个具体特征、结构或者特点时,所要主张的是结合其他实施例来实现这种特征、结构或者特点也落在本发明的范围内。
尽管这里参照本发明的多个解释性实施例对本发明进行了描述,但是,应该理解,本领域技术人员可以设计出很多其他的修改和实施方式,这些修改和实施方式将落在本申请公开的原则范围和精神之内。更具体地说,在本申请公开、附图和权利要求的范围内,可以对主题组合布局的组成部件和/或布局进行多种变型和改进。除了对组成部件和/或布局进行的变型和改进外,对于本领域技术人员来说,其他的用途也将是明显的。
Claims (5)
1.一种检测超声波在测流速中飞行时间的方法,其特征在于,包括以下步骤:
步骤1,处理器发送频率为f 1的超声波信号,经过驱动电路驱动放大以后,再由超声波发射器发射出去;同时,计数器开始计数;
步骤2,计数器计数一段时间t 0后停止计数;
步骤3,接收传感器将接收到的频率为f 1的超声波信号经模拟电路放大和滤波处理后传送给模数转换器;
步骤4,模数转换器将转化后的超声波数字信号发送给处理器,在所述处理器中完成频率为f 1的超声波数字信号的相位差 的计算;本步骤中,相位差,其中当超声波在媒质中传播一段时间后,接收到的波形相位在媒质静止情况下为,当媒质为流体时,接收到的波形相位为;
步骤5,由步骤1中所述的处理器发出频率为f 2的超声波信号,所述频率为f 2的超声波信号经过驱动电路驱动放大以后,再由超声波发射器发射出去;同时,计数器开始计数;
步骤6,所述计数器计数一段时间t 0后停止计数;
步骤7,在所述计数器停止计数后,接收传感器将接收到的频率为f 2的超声波信号经模拟电路放大和滤波处理后传送给模数转换器;
步骤8,所述模数转换器再将转化后的超声波数字信号发送给处理器,在所述处理器中完成频率为f 2的超声波数字信号的相位差计算;本步骤中,相位差,其中当超声波在媒质中传播一段时间后,接收到的波形相位在媒质静止情况下为,当媒质为流体时,接收到的波形相位为;
步骤9,根据两个频率和的相频关系推导计算没有相位偏移量时的传播时间:
……(1)
其中k和k’为…-2,-1,0,1,2…,
但是实际测量中相位会存在相位偏移量,设此相位偏移量全部产生在上且偏移量为,则传播时间为,……(2)
或根据两个频率f 1和 f 2的和频计算传播时间为,……(3)
或根据两个频率和的差频计算传播时间为,……(4)
在所述的式(1)、式(2)、式(3)、式(4)中,通过计算-的大小判断所测量的与是否在同一周期内的情况,以确定k与k’的关系,计算出Δt,得出超声波的飞行时间为t= t 0 +Δt。
2.根据权利要求1所述的检测超声波在测流速中飞行时间的方法,其特征在于,步骤1中所述的处理器为FPGA。
3.根据权利要求1所述的检测超声波在测流速中飞行时间的方法,其特征在于,步骤1中所述的处理器为CPLD。
4.根据权利要求1所述的检测超声波在测流速中飞行时间的方法,其特征在于,所述频率为f 1的超声波数字信号的相位差计算采用DFT运算方法。
5.根据权利要求1所述的检测超声波在测流速中飞行时间的方法,其特征在于,所述频率为f 2的超声波数字信号的相位差计算采用DFT运算方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210344506.5A CN102866261B (zh) | 2012-09-18 | 2012-09-18 | 检测超声波在测流速中飞行时间的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210344506.5A CN102866261B (zh) | 2012-09-18 | 2012-09-18 | 检测超声波在测流速中飞行时间的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102866261A CN102866261A (zh) | 2013-01-09 |
CN102866261B true CN102866261B (zh) | 2014-11-05 |
Family
ID=47445252
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210344506.5A Expired - Fee Related CN102866261B (zh) | 2012-09-18 | 2012-09-18 | 检测超声波在测流速中飞行时间的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102866261B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103471668B (zh) * | 2013-09-23 | 2016-04-06 | 浙江理工大学 | 一种逐级升高电压激励测量超声波渡越时间的电路 |
CN107179676B (zh) * | 2017-05-24 | 2019-04-23 | 成都理工大学 | 基于数字解调的超声波飞行时间测量方法 |
JP7012497B2 (ja) | 2017-10-04 | 2022-01-28 | 上田日本無線株式会社 | 伝搬時間測定装置、気体濃度測定装置および伝搬時間測定プログラム |
CN109782018B (zh) * | 2019-01-18 | 2021-02-19 | 济南大学 | 一种超声波风速风向测量实验仪及测量方法 |
EP3779374B1 (en) | 2019-03-20 | 2022-08-17 | Shenzhen Goodix Technology Co., Ltd. | Time of flight generation circuit, and relevant chip, flow meter and method |
CN113030248B (zh) * | 2021-03-08 | 2024-02-09 | 哈尔滨理工大学 | 基于超声波双频相位差的氢气测量系统及方法 |
CN114088972B (zh) * | 2021-11-21 | 2024-04-05 | 吉林大学 | 一种基于相位双频法的超声波风速风向测量系统及方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4538163B2 (ja) * | 2001-03-30 | 2010-09-08 | サーパス工業株式会社 | 流速測定方法及び流速測定装置並びに流量測定方法及び流量測定装置 |
CN100403055C (zh) * | 2006-03-24 | 2008-07-16 | 成都中科动态仪器有限公司 | 双频信号飞行时间测量方法及其测量装置 |
CN202083707U (zh) * | 2011-04-22 | 2011-12-21 | 哈尔滨理工大学 | 基于相位差法测量超声波风速风向的装置 |
-
2012
- 2012-09-18 CN CN201210344506.5A patent/CN102866261B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102866261A (zh) | 2013-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102866261B (zh) | 检测超声波在测流速中飞行时间的方法 | |
CN102589627B (zh) | 一种用于超声波流量计的绝对传播时间测量方法 | |
CN105157771A (zh) | 一种时差式超声波流量测量方法及装置 | |
WO2014205884A1 (zh) | 一种基于零点分析的时差式超声波流量计测量方法 | |
CN104501889B (zh) | 基于互相关时差法超声波流量的检测方法 | |
EP2653840A1 (en) | Flow volume measuring device | |
CN101598580B (zh) | 一种提高时差式超声波流量计精度的方法 | |
CN104457871A (zh) | 一种流量计及流体测量方法 | |
CN103162752A (zh) | 用于超声波流量计的相位编码同步时差检测装置及方法 | |
Tian et al. | Energy peak fitting of echo based signal processing method for ultrasonic gas flow meter | |
CN109029598A (zh) | 一种互相关法气体超声波流量计中基于主次包络结合处的参考波形选取方法 | |
CN111157065A (zh) | 气体超声流量计超声波信号传输回路中声延时测量方法 | |
CN106643937A (zh) | 一种基于超声波流量计的流量测量方法及装置 | |
CN203100824U (zh) | 一种基于声驻波法和声卡的液位测量装置 | |
Mu et al. | Echo signal envelope fitting based signal processing methods for ultrasonic gas flow-meter | |
Han et al. | Studies on the transducers of clamp-on transit-time ultrasonic flow meter | |
EP3164680B1 (en) | Method of measuring time of flight of an ultrasound pulse | |
CN203069223U (zh) | 用于超声波流量计的相位编码同步时差检测装置 | |
CN102967334A (zh) | 利用对信号包络线处理测量流体流量的系统及方法 | |
CN102360024A (zh) | 一种纸浆流速及流量的测量方法 | |
CN204165597U (zh) | 用于流量检测系统的回波信号调理电路 | |
JP2014182024A (ja) | 超音波測定器 | |
CN202974355U (zh) | 一种超声波流量测量装置 | |
RU118743U1 (ru) | Ультразвуковой расходомер | |
CN103063868B (zh) | 流体流速测量装置及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141105 Termination date: 20150918 |
|
EXPY | Termination of patent right or utility model |