CN107764369A - 超声波液位测量方法 - Google Patents

超声波液位测量方法 Download PDF

Info

Publication number
CN107764369A
CN107764369A CN201610705713.7A CN201610705713A CN107764369A CN 107764369 A CN107764369 A CN 107764369A CN 201610705713 A CN201610705713 A CN 201610705713A CN 107764369 A CN107764369 A CN 107764369A
Authority
CN
China
Prior art keywords
ultrasonic
signal
ultrasonic wave
flag
energy
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.)
Granted
Application number
CN201610705713.7A
Other languages
English (en)
Other versions
CN107764369B (zh
Inventor
郭宇
赵冬泉
李磊
郭实敏
唐兰贵
李雪森
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Qing Kong Human Settlement Research Institute Co Ltd
Original Assignee
Beijing Qing Kong Human Settlement Research Institute Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Qing Kong Human Settlement Research Institute Co Ltd filed Critical Beijing Qing Kong Human Settlement Research Institute Co Ltd
Priority to CN201610705713.7A priority Critical patent/CN107764369B/zh
Publication of CN107764369A publication Critical patent/CN107764369A/zh
Application granted granted Critical
Publication of CN107764369B publication Critical patent/CN107764369B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F23/00Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm
    • G01F23/22Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measuring physical variables, other than linear dimensions, pressure or weight, dependent on the level to be measured, e.g. by difference of heat transfer of steam or water
    • G01F23/28Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measuring physical variables, other than linear dimensions, pressure or weight, dependent on the level to be measured, e.g. by difference of heat transfer of steam or water by measuring the variations of parameters of electromagnetic or acoustic waves applied directly to the liquid or fluent solid material
    • G01F23/296Acoustic waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F23/00Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm
    • G01F23/22Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measuring physical variables, other than linear dimensions, pressure or weight, dependent on the level to be measured, e.g. by difference of heat transfer of steam or water
    • G01F23/28Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measuring physical variables, other than linear dimensions, pressure or weight, dependent on the level to be measured, e.g. by difference of heat transfer of steam or water by measuring the variations of parameters of electromagnetic or acoustic waves applied directly to the liquid or fluent solid material
    • G01F23/296Acoustic waves
    • G01F23/2965Measuring attenuation of transmitted waves

Landscapes

  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Electromagnetism (AREA)
  • Thermal Sciences (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Levels Of Liquids Or Fluent Solid Materials (AREA)

Abstract

本发明提供了一种超声波液位测量方法,包括下列步骤:控制装置控制超声波传感器向液面发送固定频率的超声波探测信号,超声波探测信号在接触液面后,向超声波传感器反弹回超声波回传信号;控制装置实时采集接收的超声波回传信号,同时在控制装置每采集一个超声波回传信号数据之后,实时动态调节可变增益放大器的输入电压来对超声波回传信号进行距离衰减动态补偿,得到超声波补偿信号;控制装置对超声波补偿信号进行加窗处理,以特征窗口快速从超声波补偿信号中提取目标频率信号。控制装置在超声波补偿信号中选取目标频率最大能量值对应的信息,并在距离能量最大值位置特征范围外,选取第二大能量值对应的信息。本发明能够给出精确的液位结果。

Description

超声波液位测量方法
技术领域
本发明涉及测量技术和传感器数据处理技术,具体地说,本发明涉及一种超声波液位测量方法。
背景技术
随着科学技术的发展,人们对雨水、污水、中水、河道、湖泊、道路等市政与水利领域的液位监测需求越来越大。近些年,我国多个城市在每年雨季都会发生内涝事件,呈现了发生范围广、积水深度大、积水时间长的特点。内涝事件的发生,不仅给人民群众在雨天出行带来诸多不便,而且导致了个人和公共财产的巨大损失,甚至威胁了人民群众的生命安全,引起了社会各方面的广泛关注。
目前在超声波液位监测方面,其它设备的功能性、稳定性和硬件集成性不能很好的满足在恶劣环境中长期监测的需求。而且价格过高、功能不足、设备零配件供货周期长等问题,严重影响了液位监测等效果。主要表现如下:
(1)现有超声波液位计没有对回传的超声波能量做距离衰减动态补偿设计,可以满足环境简单的测距需求,但是不能满足测试环境复杂的应用场合,会造成测距的不准确和跳变。
(2)现有超声波液位计采用复杂的检波方法对回波中的信号进行提取工作,造成硬件资源的浪费以及计算的复杂,从而在硬件成本比较浪费,开发周期时间长,设计工作上比较复杂。
(3)现有超声波液位计对回波中的信号进行提取单一的能量最大值运算,可以满足超声波传感器与目标物体之间的障碍物比较少,而且对超声波液位计的安装要求比较苛刻。使得超声波液位计的使用场合要求比较严格,超声波液位计的测距效果上不是很理想。
因此,亟待出现一种能够适应复杂环境的液位测量要求和简单有效的液位测量方法。
发明内容
为解决上述技术问题,本发明提供了一种能够适应复杂环境的液位测量要求和简单有效的液位测量方法。
本发明采用的技术方案为:
本发明的实施例提供一种超声波液位测量方法,包括:控制装置控制超声波传感器向被测液面发送固定频率的超声波探测信号,所述超声波探测信号在接触所述液面后,向所述超声波传感器反弹回超声波回传信号;所述控制装置实时采集接收的超声波回传信号,同时在控制装置每采集一个超声波回传信号之后,实时动态调节可变增益放大器的输入电压来对所述超声波回传信号进行距离衰减动态补偿,得到超声波补偿信号;所述控制装置对所述超声波补偿信号进行加窗处理,以特征窗口快速从所述超声波补偿信号中提取目标频率信号;所述控制装置在所述超声波补偿信号中选取目标频率最大能量值对应的信息,并在距离能量最大值位置特征范围外,选取第二大能量值对应的信息;所述控制装置基于所选取目标频率的最大能量值对应的信息和选取目标频率的第二大能量值对应的信息,得到有效的超声波反射所述液面所对应的信息,计算所述超声波传感器距离所述液面的液位高度。
优选地,所述控制装置驱动第一数字模拟转换器发送固定频率的超声波探测信号,发送的探测信号经经功率放大器放大后,通过所述超声波传感器发送;所述控制装置控制模拟数字转换器以大于所述超声波回传信号中的最高频率的2倍的采样频率来采集所述超声波回传信号;所述可变增益放大器的输入电压与对所述超声波回传信号进行距离衰减动态补偿的增益能量值成正比。
优选地,所述控制装置每采集一个超声波回传信号之后,控制装置通过控制第二数字模拟转换器输出连续的变化的电压来动态调节所述可变增益放大器的输入电压,进而动态调节所述增益能量值,以对所述超声波回传信号的下一个采集点进行距离衰减补偿。
优选地,控制装置控制超声波传感器发送的超声波探测信号中包括固定Nperiod个周期的正弦波信号,所述控制装置对所述超声波补偿信号进行分析,提取所述超声波补偿信号中Nperiod个周期正弦波信号的位置信息,并计算超声波传感器到被测液面的距离。
优选地,控制装置对所述超声波补偿信号进行加窗处理,包括利用特征窗口大小对所述超声波补偿信号分数据段进行分析,其中特征窗口大小通过如下公式获得:SIZEwindow=Nperiod*fsample/fsend,其中,SIZEwindow为信号处理窗口大小,Nperiod为超声波传感器发送波形周期数,fsample为控制装置采集信号频率,fsend为超声波探头发送信号频率。
优选地,控制装置采用1和-1两个特征值做信号提取特征点来从所述超声波补偿信号中快速提取目标频率信号的位置信息。
优选地,所述最大能量值对应的信息和所述第二大能力最大值对应的信息包括:超声波回传信号中能量最大的点数Pfmax,超声波回传信号中能量最大的点数的能量Efmax,超声波回传信号中距离能量最大值位置特征值范围外能量次大的点数Psmax,超声波回传信号中能量次大的点数的能量Esmax
优选地,还包括:控制装置对超声波回传信号中能量最大点的信息和超声波回传信号中能量次大点的信息进行过滤处理,过滤处理过程中涉及的参数包括:超声波回传信号中能量一级阀值Evalid1,超声波回传信号中能量二级阀值Evalid2,最近一次有效的超声波回传信号中能量最大的点数Plastmax,最近一次有效的超声回传信号中能量最大的点数的能量Elastmax,超声波反射物体回传信号中最大点标志位FlagEfmax,超声波回传信号中能量次大点标志位FlagEsmax,超声波回传信号中能量最大的点数与上一次测量有效值允许变化阀值Pchange,超声波回传信号中中能量存在有效值标志位Flagsuccess
优选地,控制装置基于以下过滤处理原则对超声波回传信号中能量最大点的信息和超声波回传信号中能量次大点的信息进行过滤处理:
1)
当Efmax≥Evalid1时,如果Efmax≥Evalid2,则FlagEfmax=2
当Efmax≥Evalid1时,如果|Pfmax-Plastmax|≤Pchange,则FlagEfmax=1
当Efmax≥Evalid1时,如果|Pfmax-Plastmax|>Pchange,则FlagEfmax=0
当Efmax<Evalid1时,FlagEfmax=0
2)
当Esmax≥Evalid1时,如果Esmax≥Evalid2,FlagEsmax=2
当Esmax≥Evalid1时,如果|Psmax-Plastmax|≤Pchange,FlagEsmax=1
当Esmax≥Evalid1时,如果|Psmax-Plastmax|>Pchange,FlagEsmax=0
当Esmax<Evalid1时,FlagEsmax=0
3)
当FlagEfmax=FlagEsmax时,如果FlagEfmax=0,则Flagsuccess=0;
当FlagEfmax=FlagEsmax时,如果FlagEfmax≠0且
|Pfmax-Plastmax|≤|Psmax-Plastmax|,则Flagsuccess=1,Plastmax=Pfmax,Elastmax=Efmax
当FlagEfmax=FlagEsmax时,如果FlagEfmax≠0且
|Pfmax-Plastmax|>|Psmax-Plastmax|,则Flagsuccess=1,Plastmax=Psmax,Elastmax=Esmax
当FlagEfmax≠FlagEsmax时,Flagsuccess=1,Plastmax=Pfmax,Elastmax=Efmax
优选地,控制装置根据所述过滤处理原则得到Plastmax,Elastmax和Flagsuccess的具体数值,并根据下述公式来计算超声波传感器距离液面的距离:
H2=Plastmax*Vspeed/2/fsample其中:H2为超声波传感器距离液面的距离;Plastmax为最近一次有效的超声波回传信号中能量最大的点数;Vspeed为声音在空气中的传播速度;fsample为模拟数字转换器的采样率。
本发明的超声波液位测量方法对回传的超声波能量做距离衰减动态补偿设计,可以满足测试环境复杂的应用场合,能够获得精确稳定的测距;同时采用了简单的1和-1两个特征值做信号提取特征点来从所述超声波补偿信号中快速提取目标频率信号的位置信息,使得计算过程简单化,节约了处理器的资源。此外,提取了最大能量值和第二大能量最大值两个值,同时利用过滤原则对最大能量值对应的信息和第二大能量最大值对应的信息进行了过滤,从而可以提高测距的抗扰能力。
附图说明
图1示出了本发明一个实施例的超声波液位测量方法的流程图。
图2示出了本发明一个实施例的测距装置的硬件模块示意图。
图3示出了本发明一个实施例的测距设备的现场使用安装图。
图4示出了超声波回波信号没有进行距离动态补偿的原始信号波形。
图5示出了超声波回波信号进行距离动态补偿的原始信号波形。
图6示出了超声波回波信号没有进行距离动态补偿提取能量值后的波形。
图7示出了超声波回波信号进行动态能量补偿提取能量值后的波形。
图8示出了超声波回波信号进行动态能量补偿的原始信号波形。
图9示出了从超声波回波信号中快速提取能量值后的波形。
具体实施方式
以下结合附图对本发明的具体实施方式进行描述。
图1示出了本发明一个实施例的超声波测距方法的流程图,图2示出了本发明一个实施例的测距装置的硬件模块示意图,图3示出了本发明一个实施例的测距设备的现场使用安装图。
首先,参考图2和图3对本发明中应用到的测距装置进行描述。如图3所示,作为一示意性的实施例,可将具有处理功能的控制装置和超声波传感器安装在井的上部,以对井的液位进行监测,当然,本发明并不局限于对距离液面的液位进行测量,还可以对距离其他目标物体的高度进行测量。在图3中,各标记的含义如下:
1)H1:超声传感器底部距井口距离=传感器支架距井口距离+传感器支架距超声传感器底部距离。
2)H2:超声传感器底部距水面距离。
3)H3:水深高度。
4)H4:井深高度。
5)H5:预警液位。
6)H6:报警液位。
在本发明的一实施例中,H1=0.78米,H4=8.12米,H5=2.00米,H6=5.00米,液面高度H3=1.51米。本发明的计算目标为H2,随后将进行介绍。
图2示出了本发明一个实施例的测距装置硬件模块示意图,应用硬件模块主要包括:处理器(MCU),可编程逻辑器件(CPLD),数字模拟转换器(DA1),功率放大器(PA)、数字模拟转换器(DA2),模拟数字转换器(AD),可变增益放大器(VGA),收发一体超声波传感器,其中处理器和可编程逻辑器件可构成测距装置的控制装置,对其余模块进行控制,各模块可按照图2所示出的关系进行通讯连接,本发明也可使用两个超声波传感器来发送和接收超声波信号,而不仅仅使用收发一体超声波传感器。具体地,在本发明中,处理器主要起控制作用,处理器与可编程逻辑器件之间通过串行外设接口(SPI)实现控制命令的下发与数据的传输。可编程逻辑器件CPLD接收处理器的发波命令,控制数字模拟转换器DA1发送信号,经过功率放大器驱动超声波传感器向目标物体(本发明中为液面)发送超声波,发送的超声波在碰到液面后反弹回超声波传感器;可编程逻辑器件CPLD接收处理器发送的补偿命令,控制可变增益放大器VGA对回传的超声波信息进行距离衰减动态补偿,以尽可能的消除液面周边环境对回传的超声波信号产生的干扰;并且,可编程逻辑器件接收处理器发送的信号采集命令,控制模拟数字转换器采集超声波回传信号,可编程逻辑器件对采集的超声波回传信号进行分析处理,得到最大能量值对应的信息和第二大能量值对应的信息,可编程逻辑器件将该最大能量值对应的信息和第二大能量值对应的信息回传给处理器,处理器最后得到超声波传感器与目标物体的距离数值。
接下来,结合图1、图4至图9(参考图2和图3)和优选实施例对本发明的超声波液位测量方法进行介绍。
如图1所示,本发明的超声波液位测量方法包括以下步骤:
步骤101:测量***初始化,***中的各个模块上电,开始测距过程。具体地,处理器通过串行外设接口(SPI)与可编程逻辑器件(CPLD)通信,可编程逻辑器件(CPLD)初始化数字模拟转换器(DA1),将数字模拟转换器(DA1)的输出电压调到中心位置,功率放大器(PA)的输出调到中心位置,处理器通过串行外设接口(SPI)与可编程逻辑器件通信,可编程逻辑器件(CPLD)初始化数字模拟转换器(DA2),将数字模拟转换器(DA2)的输出电压调到零点位置,数字模拟转换器(DA2)调节可变增益放大器(VGA)对超声波传感器回传回传信号进行零补偿处理。同时,处理器通过串行外设接口(SPI)与可编程逻辑器件通信,可编程逻辑器件(CPLD)初始化模拟数字转换器(AD),将模拟数字转换器(AD)工作在非低功耗模式,使得模拟数字转换器(AD)开始工作。上述所有模块准备好,等待处理器向可编程逻辑器件发送测距命令。
步骤102:可编程逻辑器件中预存了16个正弦波特征点,这些正弦波特征点按顺序画成曲线,可以绘制成一个周期的正弦波曲线。处理器通过SPI接口向可编程逻辑器件发送测距命令,可编程逻辑器件驱动数字模拟转换器(DA1)发送16个周期频率为40Khz的正弦波,数字模拟转换器(DA1)将该波形输出给功率放大器(PA),功率放大器PA将接收的正弦波信号的功率进行放大,驱动超声波传感器向被测物体(液面)发送超声波探测信号。
步骤103:在超声波传感器向被测液面发送超声波探测信号后,处理器通过SPI接口向可编程逻辑器件发送信号补偿命令,可编程逻辑器调节数字模拟转换器(DA2)的输入数字量,驱动VGA动态补偿超声波能量的距离衰减,本发明的数字模拟转换器DA2输入的初始数字量为15000,后一个数字量以步长大小为5逐个递增,数字模拟转换器DA2输入的数字量总个数为8000,但并不局限于此,例如,也可以以16个点为一个周期,前9个点输出递增5,后7个点输出递增6;同时,可编程逻辑器件通过模拟数字转换器AD采集补偿的超声波回传信号,采样频率应符合奈奎斯特定理,因为在进行模拟/数字信号的转换过程中,当采样频率fsmax大于回传信号中最高频率fmax的2倍时(fsmax>2fmax),采样之后的数字信号完整地保留了原始信号中的信息。本发明中的模拟数字转换器AD在进行模拟数字信号的转换过程中,一实施例示出的采样频率为最高回传信号频率的4倍,即采样频率为160Khz。
在本发明中,超声波传感器发送固定Nperiod个周期超声波信号,这些信号接触到目标物体后,会反弹回超声波,反弹回的超声波信号回传给超声波传感器。超声波在空气中回传的过程中,随着传输距离的增加,超声波的总能量值会衰减,并且超声波传感器接收到的能量值与超声波传输距离成反比,通过公式表示为:Esound=Et–K1*Sdistance。其中:Esound为超声波传感器发送的超声波,与前方物体接触,反弹回来的超声波能量值;K1为超声波在空气之中传播,随着距离的增加,超声波能量值的衰减斜率;Sdistance为超声波传感器与前方物体的距离,为超声波在空气中经过总路径的一半;Et为超声波传感器发射超声波总能量。
对回传的超声波信号的补偿具体为:CPLD调节可变增益放大器(VGA)的输入电压,实时动态的补偿超声波距离衰减的能量值,可变增益放大器的动态补偿能量值与可变增益放大器的输入电压成正比,通过公式表示为:Evga=K2*VDA2。其中:Evga为可变增益放大器增益值;K2为可变增益放大器增益斜率;VDA2为数字模拟转换器DA2输出电压值。
其中,可变增益放大器VGA的输入电压,与数字模拟转换器DA2的输入数字值成正比,通过公式表示为:VDA2=VDA2ref*DDAin/Dnum,其中,VDA2ref为数字模拟转换器DA2基准电压;DDain为数字模拟转换器DA2的输出数字量数值;VDA2为可变增益放大器VGA的输入电压;Dnum为数字模拟转换器DA2输出数字量总值。
可通过调节可变增益器(VGA)的输入电压的大小来动态的补偿距离衰减的能量值。距离补偿因子Vm=K1/K2,其中,K1为超声波能量值的距离衰减斜率,K2为为可变增益放大器增益斜率,Vm为VGA距离补偿因子。
可通过调节数字模拟转换器DA2的输入数字量大小,来调节可变增益放大器VGA的输入电压。数字补偿因子Dm=Vm*Dnum/VDA2ref,其中:Vm为距离补偿因子,Dm为数字模拟转换器DA2的数字补偿因子。
在本发明中,通过控制数字模拟转换器DA2的电压来对超声波距离衰减动态补偿方法为Dvga=Dt/f,其中:Dvga为DA2每个点输出的数字增量,Dt为DA2每秒输出的数字增量,f为DA2输出频率,在本发明中,VGA距离衰减动态补偿的补偿频率为40Khz。
即,可编程逻辑器件CPLD通过调节向数字模拟转换器DA2输入的数字量大小来调节DA2的输出电压,进而调节可变增益放大器VGA的输入电压,最终对超声波回传信号的能量值进行动态补偿。
图4中示出了超声波回波信号没有进行距离动态补偿的原始信号图形。图5示出了超声波回波信号进行动态能量补偿的原始信号图形。图4和图5所示超声波原始回波信号都是在800点(干扰信号)和5450点(液位信号)附近波形幅值比较大,说明在这两处有反射物体回波信号。图4和图5的区别在于超声波原始回波信号进行动态能量补偿后,超声波原始回波信号可以弥补随着距离造成的能量衰减。将在步骤104中描述图4和图5所示出的信号的快速提取能量值的结果。
步骤104:可编程逻辑器件实时的处理模拟数字转换器AD采集回来的超声波补偿回传信号,对补偿信号进行加窗处理,采用的特征窗口大小的计算公式为:SIZEwindow=Nperiod*fsample/fsend,其中:SIZEwindow为信号处理窗口大小,Nperiod为超声波传感器发送的波形周期数,fsample为处理器采集信号的频率,fsend为超声波探头发送的信号频率,在本发明中,由于Nperiod为16,fsample为160Khz,fsend为40Khz,这样得到信号处理窗口的大小为64。同时可编程逻辑器件快速提取40Khz的目标频率信号的能量值,也就是将补偿信号的每个点与对应正弦波特征点进行乘法运算,然后把这个信号点的乘法结果和后面信号点的乘法结果进行求和运算,求和运算的次数为窗口大小,最后得到这个信号点的能量值。本发明采用1和-1两个特征值做检波信号的正弦波特征点,每个周期取4个特征点,也就是1,1,-1,-1;信号中每个数据乘以1的计算步骤可以简化为取信号本身,信号中每个数据乘以-1的计算步骤可以简化为取信号的相反数,然后,将特征窗口内的64个信号点的运算结果进行求和运算,得到该补偿信号点的能量值,如此类推,可得到每个信号点的能量值,每个特征窗口内的能量值所对应的位置信息为该窗口开始的第一个信号点。本发明采用1,1,-1,-1作为特征值,既简化了计算过程,又节约了处理器的资源。
图6示出了超声波回波信号没有进行距离动态补偿提取能量值后的波形。图7示出了超声波回波信号进行动态能量补偿提取能量值后的波形。图6和图7所示的超声波回波信号提取能量值后的波形都是在800(干扰信号)点和5450(液位信号)点附近出现峰值,说明超声波信号在这两处有反射。图6和图7的区别在于,图6所示离超声波传感器近的干扰物会对测距造成错误影响,图7所示超声波传感器近的干扰物会对测距造微小影响,图6和图7对比说明超声波回波信号进行动态能量补偿后可以提高测距的抗干扰能力。
图8示出了超声波回波信号进行动态能量补偿的原始信号图形。超声波原始回波信号在3750点(流槽干扰信号)和5450点(液位信号)附近波形幅值比较大,说明在这两处有反射物体回波信号。图9示出了图8超声波回波信号中快速提取能量值后的波形,步骤105中描述了图9波形取双峰,过滤算法后的结果。
步骤105:可编程逻辑器件在补偿信号的能量值中找出最大能量值和第二大能量最大值,并提取他们的有效位置信息,包括:超声波反射物体回波(回传信号)中能量最大的点数Pfmax,超声波反射物体回波中能量最大的点数的能量Efmax,超声波反射物体回波中距离能量最大值位置特征值范围外能量次大的点数Psmax,超声波反射物体回波中能量次大的点数的能量Esmax。从图9示出超声波回波信号能量值的波形中提取双峰,可以得到Pfmax=3777,Efmax=46990000,Psmax=5493,Esmax=39150000。
可编程逻辑器件将提取的Pfmax,Efmax,Psmax,Esmax数值通过SPI接口传给处理器,接着,处理器对Pfmax,Efmax,Psmax,Esmax进行过滤处理,过滤处理过程中涉及的参数包括:超声波反射物体回波能量一级阀值Evalid1,超声波反射物体回波能量二级阀值Evalid2,最近一次有效的超声波反射物体回波中能量最大的点数Plastmax,最近一次有效的超声波反射物体回波中能量最大的点数的能量Elastmax;超声波反射物体回波中能量最大点标志位FlagEfmax,超声波反射物体回波中能量次大点标志位FlagEsmax;超声波反射物体回波中能量最大的点数与上一次测量有效值允许变化阀值Pchange;超声波反射物体回波中能量存在有效值标志位Flagsuccess。其中,Evalid1、Evalid2、Elastmax、Pchange的值存储在处理器中,为定值。超声波反射物体回波中能量最大点标志位FlagEfmax和超声波反射物体回波中能量次大点标志位FlagEsmax反应的是能量最大值和能量次大值的能量等级,在本发明中,采用0、1和2作为标志位的能量等级,以确定选取能量最大值和能量次大值中的那个作为有效值,其中0表示最低级,2表示最高级。超声波反射物体回波中能量存在有效值标志位Flagsuccess反应的是本次测距的数据有效性。0表示本次测距无效,1表示本次测距有效。
步骤106:处理器对超声波反射物体回波中能量最大点的信息和超声波反射物体回波中能量次大点的信息进行过滤处理,从超声波反射物体回波中能量最大点和超声波反射物体回波中能量次大点中找出有效的超声波反射物体回波中的点的信息,过滤处理按照如下原则:
(1)
当Efmax≥Evalid1时,如果Efmax≥Evalid2,则FlagEfmax=2
当Efmax≥Evalid1时,如果|Pfmax-Plastmax|≤Pchange,则FlagEfmax=1
当Efmax≥Evalid1时,如果|Pfmax-Plastmax|>Pchange,则FlagEfmax=0
当Efmax<Evalid1时,FlagEfmax=0
(2)
当Esmax≥Evalid1时,如果Esmax≥Evalid2,FlagEsmax=2
当Esmax≥Evalid1时,如果|Psmax-Plastmax|≤Pchange,FlagEsmax=1
当Esmax≥Evalid1时,如果|Psmax-Plastmax|>Pchange,FlagEsmax=0
当Esmax<Evalid1时,FlagEsmax=0
(3)
当FlagEfmax=FlagEsmax时,如果FlagEfmax=0,则Flagsuccess=0;
当FlagEfmax=FlagEsmax时,如果FlagEfmax≠0且
|Pfmax-Plastmax|≤|Psmax-Plastmax|,则Flagsuccess=1,Plastmax=Pfmax,Elastmax=Efmax
当FlagEfmax=FlagEsmax时,如果FlagEfmax≠0且
|Pfmax-Plastmax|>|Psmax-Plastmax|,则Flagsuccess=1,Plastmax=Psmax,Elastmax=Esmax
当FlagEfmax≠FlagEsmax时,Flagsuccess=1,Plastmax=Pfmax,Elastmax=Efmax
上述过滤原则(1)用于确定能量最大点标志位FlagEfmax,即确定能量最大值的等级,上述过滤原则(2)用于确定能量次大点标志位FlagEsmax,即确定能量次大值的等级,然后根据过滤原则(1)和(2)的结果,来选取有效标志位,并决定将能量最大点数及其能量值和能量次大值点数及其能量值中的那个作为最近一次有效的超声波反射物体回波点数和能量值。Flagsuccess=1表示测距成功,Flagsuccess=0表示测距不成功,测距成功能够得到有效的液位信息,而测试不成功则不能得到有效的液位信息。
在本发明的一实施例中,Evalid1=10000000,Evalid2=30000000,Plastmax=5470,Elastmax=32230000,Pchange=50。这样,结合图8和图9所示出的波形,单纯的提取最大能量值位置信息作计算,得到的结果为:Flagsuccess=1,Plastmax=3777,Elastmax=46990000;单纯的提取最大能量值位置信息作计算,没有提取两个峰值做过滤处理,得到流槽干扰信号的位置信息,为错误液位信息。取双峰经过滤处理后的结果为:Flagsuccess=1,Plastmax=5493,Elastmax=39150000,为正确液位信息。根据上述结果可知,选取了最大点能量值的位置信息和能量次大点的位置信息并对它们进行过滤处理后,得到的Plastmax和Elastmax与最近一次有效的Plastmax和Elastmax较相近,由此可知,超声波反射物体回波中能量最大点的信息和超声波反射物体回波中能量次大点的信息经过上述过滤原则过滤处理后,能够提高测距的抗扰能力。如图9所示,提取两个两个峰值做过滤处理,得到正确液位,液位的具体算法参见下述步骤107。
步骤107:计算液位H2
处理器根据过滤处理得到有效的超声波反射物体回波中的点的信息Plastmax,采用公式H2=Plastmax*Vspeed/2/fsample来计算液位H2。其中:Plastmax为最近一次有效的超声波反射物体回波中能量最大的点数,Vspeed为声音在空气中的传播速度,fsample为模拟数字转换器AD的采样频率。
在本发明的一实施例中,根据步骤106计算的Plastmax、处理器采集信号的频率fsample以及Vspeed来计算H2,最终得到H2=5.83米。本次液位测量完成,各模块关电。
以上所描述的装置实施例仅仅是示意性的,其中作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分的方法。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (10)

1.一种超声波液位测量方法,其特征在于,包括:
控制装置控制超声波传感器向被测液面发送固定频率的超声波探测信号,所述超声波探测信号在接触所述液面后,向所述超声波传感器反弹回超声波回传信号;
所述控制装置实时采集接收的超声波回传信号,同时在控制装置每采集一个超声波回传信号之后,实时动态调节可变增益放大器的输入电压来对所述超声波回传信号进行距离衰减动态补偿,得到超声波补偿信号;
所述控制装置对所述超声波补偿信号进行加窗处理,以特征窗口快速从所述超声波补偿信号中提取目标频率信号;
所述控制装置在所述超声波补偿信号中选取目标频率最大能量值对应的信息,并在距离能量最大值位置特征范围外,选取第二大能量值对应的信息;
所述控制装置基于所选取目标频率的最大能量值对应的信息和选取目标频率的第二大能量值对应的信息,得到有效的超声波反射所述液面所对应的信息,计算所述超声波传感器距离所述液面的液位高度。
2.根据权利要求1所述的超声波液位测量方法,其特征在于,所述控制装置驱动第一数字模拟转换器发送固定频率的超声波探测信号,发送的探测信号经经功率放大器放大后,通过所述超声波传感器发送;
所述控制装置控制模拟数字转换器以大于所述超声波回传信号中的最高频率的2倍的采样频率来采集所述超声波回传信号;
所述可变增益放大器的输入电压与对所述超声波回传信号进行距离衰减动态补偿的增益能量值成正比。
3.根据权利要求2所述的超声波液位测量方法,其特征在于,所述控制装置每采集一个超声波回传信号之后,控制装置通过控制第二数字模拟转换器输出连续的变化的电压来动态调节所述可变增益放大器的输入电压,进而动态调节所述增益能量值,以对所述超声波回传信号的下一个采集点进行距离衰减补偿。
4.根据权利要求1所述的超声波液位测量方法,其特征在于,控制装置控制超声波传感器发送的超声波探测信号中包括固定Nperiod个周期的正弦波信号,所述控制装置对所述超声波补偿信号进行分析,提取所述超声波补偿信号中Nperiod个周期正弦波信号的位置信息,并计算超声波传感器到被测液面的距离。
5.根据权利要求4所述的超声波液位测量方法,其特征在于:控制装置对所述超声波补偿信号进行加窗处理,包括利用特征窗口大小对所述超声波补偿信号分数据段进行分析,其中特征窗口大小通过如下公式获得:SIZEwindow=Nperiod*fsample/fsend,其中,SIZEwindow为信号处理窗口大小,Nperiod为超声波传感器发送波形周期数,fsample为控制装置采集信号频率,fsend为超声波探头发送信号频率。
6.根据权利要求4所述的超声波液位测量方法,其特征在于,控制装置采用1和-1两个特征值做信号提取特征点来从所述超声波补偿信号中快速提取目标频率信号的位置信息。
7.根据权利要求1所述的超声波液位测量方法,其特征在于,所述最大能量值对应的信息和所述第二大能力最大值对应的信息包括:超声波回传信号中能量最大的点数Pfmax,超声波回传信号中能量最大的点数的能量Efmax,超声波回传信号中距离能量最大值位置特征值范围外能量次大的点数Psmax,超声波回传信号中能量次大的点数的能量Esmax
8.根据权利要求7所述的超声波液位测量方法,其特征在于,还包括:控制装置对超声波回传信号中能量最大点的信息和超声波回传信号中能量次大点的信息进行过滤处理,过滤处理过程中涉及的参数包括:超声波回传信号中能量一级阀值Evalid1,超声波回传信号中能量二级阀值Evalid2,最近一次有效的超声波回传信号中能量最大的点数Plastmax,最近一次有效的超声回传信号中能量最大的点数的能量Elastmax,超声波反射物体回传信号中最大点标志位FlagEfmax,超声波回传信号中能量次大点标志位FlagEsmax,超声波回传信号中能量最大的点数与上一次测量有效值允许变化阀值Pchange,超声波回传信号中中能量存在有效值标志位Flagsuccess
9.根据权利要求8所述的超声波液位测量方法,其特征在于,控制装置基于以下过滤处理原则对超声波回传信号中能量最大点的信息和超声波回传信号中能量次大点的信息进行过滤处理:
1)
当Efmax≥Evalid1时,如果Efmax≥Evalid2,则FlagEfmax=2
当Efmax≥Evalid1时,如果|Pfmax-Plastmax|≤Pchange,则FlagEfmax=1
当Efmax≥Evalid1时,如果|Pfmax-Plastmax|>Pchange,则FlagEfmax=0
当Efmax<Evalid1时,FlagEfmax=0
2)
当Esmax≥Evalid1时,如果Esmax≥Evalid2,FlagEsmax=2
当Esmax≥Evalid1时,如果|Psmax-Plastmax|≤Pchange,FlagEsmax=1
当Esmax≥Evalid1时,如果|Psmax-Plastmax|>Pchange,FlagEsmax=0
当Esmax<Evalid1时,FlagEsmax=0
3)
当FlagEfmax=FlagEsmax时,如果FlagEfmax=0,则Flagsuccess=0;
当FlagEfmax=FlagEsmax时,如果FlagEfmax≠0且
|Pfmax-Plastmax|≤|Psmax-Plastmax|,则Flagsuccess=1,Plastmax=Pfmax,Elastmax=Efmax
当FlagEfmax=FlagEsmax时,如果FlagEfmax≠0且
|Pfmax-Plastmax|>|Psmax-Plastmax|,则Flagsuccess=1,Plastmax=Psmax,Elastmax=Esmax
当FlagEfmax≠FlagEsmax时,Flagsuccess=1,Plastmax=Pfmax,Elastmax=Efmax
10.根据权利要求9所述的超声波液位测量方法,其特征在于,控制装置根据所述过滤处理原则得到Plastmax,Elastmax和Flagsuccess的具体数值,并根据下述公式来计算超声波传感器距离液面的距离:
H2=Plastmax*Vspeed/2/fsample其中:H2为超声波传感器距离液面的距离;Plastmax为最近一次有效的超声波回传信号中能量最大的点数;Vspeed为声音在空气中的传播速度;fsample为模拟数字转换器的采样率。
CN201610705713.7A 2016-08-23 2016-08-23 超声波液位测量方法 Active CN107764369B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610705713.7A CN107764369B (zh) 2016-08-23 2016-08-23 超声波液位测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610705713.7A CN107764369B (zh) 2016-08-23 2016-08-23 超声波液位测量方法

Publications (2)

Publication Number Publication Date
CN107764369A true CN107764369A (zh) 2018-03-06
CN107764369B CN107764369B (zh) 2019-11-19

Family

ID=61264526

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610705713.7A Active CN107764369B (zh) 2016-08-23 2016-08-23 超声波液位测量方法

Country Status (1)

Country Link
CN (1) CN107764369B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109374921A (zh) * 2018-11-24 2019-02-22 南华机电(太仓)有限公司 一种超声波风速风向测量装置及方法
CN110686386A (zh) * 2019-10-12 2020-01-14 珠海格力电器股份有限公司 区域控制方法、装置、***及空调***
CN111076793A (zh) * 2019-12-30 2020-04-28 江河瑞通(北京)技术有限公司 窨井专用超声波液位测量装置及方法
CN111308472A (zh) * 2020-03-27 2020-06-19 浙江清环智慧科技有限公司 超声波测距方法、装置、***、电子设备及存储介质
CN111650397A (zh) * 2020-04-02 2020-09-11 浙江清环智慧科技有限公司 超声波回波信号的处理方法、装置、流量计及存储介质
CN112630781A (zh) * 2020-11-27 2021-04-09 深圳数联天下智能科技有限公司 超声波测距装置及超声波测距方法
US11268845B2 (en) * 2019-08-30 2022-03-08 National Applied Research Laboratories Liquid level monitoring system and method

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5428984A (en) * 1993-08-30 1995-07-04 Kay-Ray/Sensall, Inc. Self test apparatus for ultrasonic sensor
CN1350164A (zh) * 2000-10-25 2002-05-22 北京特力声技术有限责任公司 超声波测量大型油罐液位的方法
CN101464171A (zh) * 2007-12-18 2009-06-24 深圳职业技术学院 一种超声波流量检测***及检测方法
CN102636252A (zh) * 2012-04-10 2012-08-15 吉林大学 一种超声波到达精确时刻检测的方法及装置
CN102944288A (zh) * 2012-11-05 2013-02-27 中环天仪股份有限公司 具有多量程自适应功能的超声波液位计
CN104132713A (zh) * 2014-05-22 2014-11-05 深圳万讯自控股份有限公司 一种导波雷达液位计回波信号处理方法及装置
CN105241459A (zh) * 2015-10-19 2016-01-13 东南大学 一种用于室内水下目标定位的时延估计方法及装置
CN105758499A (zh) * 2016-04-18 2016-07-13 南京理工大学 基于超声脉冲回波法的非接触式液体液位检测***及方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5428984A (en) * 1993-08-30 1995-07-04 Kay-Ray/Sensall, Inc. Self test apparatus for ultrasonic sensor
CN1350164A (zh) * 2000-10-25 2002-05-22 北京特力声技术有限责任公司 超声波测量大型油罐液位的方法
CN101464171A (zh) * 2007-12-18 2009-06-24 深圳职业技术学院 一种超声波流量检测***及检测方法
CN102636252A (zh) * 2012-04-10 2012-08-15 吉林大学 一种超声波到达精确时刻检测的方法及装置
CN102944288A (zh) * 2012-11-05 2013-02-27 中环天仪股份有限公司 具有多量程自适应功能的超声波液位计
CN104132713A (zh) * 2014-05-22 2014-11-05 深圳万讯自控股份有限公司 一种导波雷达液位计回波信号处理方法及装置
CN105241459A (zh) * 2015-10-19 2016-01-13 东南大学 一种用于室内水下目标定位的时延估计方法及装置
CN105758499A (zh) * 2016-04-18 2016-07-13 南京理工大学 基于超声脉冲回波法的非接触式液体液位检测***及方法

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109374921A (zh) * 2018-11-24 2019-02-22 南华机电(太仓)有限公司 一种超声波风速风向测量装置及方法
US11268845B2 (en) * 2019-08-30 2022-03-08 National Applied Research Laboratories Liquid level monitoring system and method
CN110686386A (zh) * 2019-10-12 2020-01-14 珠海格力电器股份有限公司 区域控制方法、装置、***及空调***
CN110686386B (zh) * 2019-10-12 2021-09-14 珠海格力电器股份有限公司 区域控制方法、装置、***及空调***
CN111076793A (zh) * 2019-12-30 2020-04-28 江河瑞通(北京)技术有限公司 窨井专用超声波液位测量装置及方法
CN111076793B (zh) * 2019-12-30 2021-12-21 江河瑞通(北京)技术有限公司 窨井专用超声波液位测量装置及方法
CN111308472A (zh) * 2020-03-27 2020-06-19 浙江清环智慧科技有限公司 超声波测距方法、装置、***、电子设备及存储介质
CN111650397A (zh) * 2020-04-02 2020-09-11 浙江清环智慧科技有限公司 超声波回波信号的处理方法、装置、流量计及存储介质
CN112630781A (zh) * 2020-11-27 2021-04-09 深圳数联天下智能科技有限公司 超声波测距装置及超声波测距方法

Also Published As

Publication number Publication date
CN107764369B (zh) 2019-11-19

Similar Documents

Publication Publication Date Title
CN107764369B (zh) 超声波液位测量方法
CN106886015B (zh) 一种多波束声呐主要声学指标的检测装置及检测方法
CN102095485B (zh) 一种便携式超声声场参数测量***及其测量方法
CN101595396A (zh) 高分辨率超声位移测量的装置与方法
WO2007010460A3 (en) Apparatus for the detection of heart activity
US6530890B2 (en) Ultrasound diagnostic apparatus and method for measuring blood flow velocity using doppler effect
CN103444613A (zh) 一种用于鱼类养殖的投喂控制***及方法
CN102590345A (zh) 沉积物辨别装置、超声波探测机、参数设定方法及其程序
CN109029602B (zh) 基于超声波的流量测量方法及流量计
CN103499374A (zh) 一种基于神经网络的超声波动态液位检测方法和***
CN106568988A (zh) 一种排水管道水体速度的测量***
CN103499954B (zh) 鱼类养殖投饲控制***及方法
CN110296913B (zh) 一种可燃粉尘扩散动态浓度的检测***及其检测方法
CN201884023U (zh) 钻孔沉渣厚度超声测量***
CN102928513B (zh) 西瓜成熟度超声无损检测装置
CN107064304A (zh) 一种水果结构无损检测装置及方法
CN108614036B (zh) 一种锅炉管道内氧化皮的分级检测方法和***
CN112525785B (zh) 一种基于新型多频超声探头的悬移质含沙量测量方法
CN106093181A (zh) 磁特性参数、磁巴克豪森、磁声发射检测方法和装置
CN202075264U (zh) 用于桩基超声波检测的超声波传感器
CN109541026A (zh) 一种反射式金属蠕变性能的非线性超声检测***和检测方法
CN108303357A (zh) 基于多频水声信号的悬移质测量***及信号处理方法
CN112799072A (zh) 一种超声波测距传感器和测距方法
CN100485413C (zh) 扫频信号的时频互相关检测和时延估计方法
CN110006996B (zh) 金属复合材料超声波测厚方法及层间界面在超声波测厚中的应用

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant