CN114814285A - 一种超声波测流方法 - Google Patents

一种超声波测流方法 Download PDF

Info

Publication number
CN114814285A
CN114814285A CN202210719013.9A CN202210719013A CN114814285A CN 114814285 A CN114814285 A CN 114814285A CN 202210719013 A CN202210719013 A CN 202210719013A CN 114814285 A CN114814285 A CN 114814285A
Authority
CN
China
Prior art keywords
ultrasonic
ultrasonic transducer
pipeline
liquid
flow rate
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
CN202210719013.9A
Other languages
English (en)
Other versions
CN114814285B (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.)
Shenyang Jiade Lianyi Energy Technology Co ltd
Original Assignee
Shenyang Jiade Lianyi Energy Technology 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 Shenyang Jiade Lianyi Energy Technology Co ltd filed Critical Shenyang Jiade Lianyi Energy Technology Co ltd
Priority to CN202210719013.9A priority Critical patent/CN114814285B/zh
Publication of CN114814285A publication Critical patent/CN114814285A/zh
Application granted granted Critical
Publication of CN114814285B publication Critical patent/CN114814285B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/24—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave
    • G01P5/245—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave by measuring transit time of acoustical waves

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Volume Flow (AREA)

Abstract

本发明涉及流体测量技术领域，尤其涉及一种超声波测流方法，所述方法包括：S1、获取设置在管道上的超声换能器组在预先设定时间段内所对应的第一传播时间信息，并基于该第一传播时间信息获取管道中液体的第一流速；S2、获取设置在管道上的超声波发射组当前所对应的第二传播时间信息；S3、基于所述第一流速、所述第二传播时间信息以及管道中液体的温度，确定管道中液体的最终流速。由于采用预先设定时间段内所对应的第一传播时间信息，并基于该第一传播时间信息获取管道中液体的第一流速，并进一步基于所述第一流速、所述第二传播时间信息以及管道中液体的温度，确定管道中液体的最终流速，使得管道中液体的最终流速更为精确。

Description

一种超声波测流方法
技术领域
本发明涉及流体测量技术领域，尤其涉及一种超声波测流方法。
背景技术
现有技术中超声波测流的方法有很多，比如时差法、相差法、多普勒法等，这些方法各有各的特点。
如图4所示，现有的时差法，是使用公式（A）来求出管道内液体的流速v，其中，公式 （A）为：
Figure 543534DEST_PATH_IMAGE003
，其中，已知L为超声波传输长度，
Figure 949107DEST_PATH_IMAGE004
为超声波传输路径 与管道夹角，T 12为从第一换能器1到第二换能器2的超声波的飞行时间，T 21为从第二换能器 2到第一换能器1的超声波的飞行时间，求管道内的液体流速v，由公式（A）可以看出管道内 液体的流速v是等于一个常量
Figure 591441DEST_PATH_IMAGE005
和一个关于时间的多项式
Figure 195598DEST_PATH_IMAGE006
的乘积，在这 个时间的多项式
Figure 897975DEST_PATH_IMAGE006
的分子是从第一换能器1到第二换能器2的超声波的飞行时 间T 12与从第二换能器2到第一换能器1的超声波的飞行时间T 21的差值，分母是从第一换能 器1到第二换能器2的超声波的飞行时间T 12和从第二换能器2到第一换能器1的超声波的飞 行时间T 21两个时间的乘积，在实际应用中，从第一换能器1到第二换能器2的超声波的飞行 时间T 12与从第二换能器2到第一换能器1的超声波的飞行时间T 21的差值的数值一般是在皮 秒的量级上（几十皮秒到几百皮秒），并且从第一换能器1到第二换能器2的超声波的飞行时 间T 12与从第二换能器2到第一换能器1的超声波的飞行时间T 21的差值是分子，所以该分子 的准确度直接影响到流速v的精确与否，而相对于时间的多项式
Figure 557888DEST_PATH_IMAGE006
的分母，实际 应用中从第一换能器1到第二换能器2的超声波的飞行时间T 12和从第二换能器2到第一换能 器1的超声波的飞行时间T 21一般都在几十微秒的量级上，并且从第一换能器1到第二换能器 2的超声波的飞行时间T 12和从第二换能器2到第一换能器1的超声波的飞行时间T 21之间还 是相乘的关系，那么分母的数值的量级由于相对于分子很大，所以，最终对于流速v的精确 度影响要小于分子的影响，但是，现有的超声波测流的时差法中都是根据一个从第一换能 器1到第二换能器2的超声波的飞行时间T 12和从第二换能器2到第一换能器1的超声波的飞 行时间T 21的时间差去计算流速v，在实际应用中难免会有误差，并不能保证这个从第一换能 器1到第二换能器2的超声波的飞行时间T 12和从第二换能器2到第一换能器1的超声波的飞 行时间T 21的时间差是准确，从而导致管道中液体流速v不能被准确得到。
发明内容
要解决的技术问题
鉴于现有技术的上述缺点、不足，本发明提供一种超声波测流方法，其解决了现有技术中管道中液体流速v测量不准确的技术问题。
技术方案
为了达到上述目的，本发明采用的主要技术方案包括：
本发明实施例提供一种超声波测流方法，包括：
S1、获取设置在管道上的超声换能器组在预先设定时间段内所对应的第一传播时间信息，并基于该第一传播时间信息获取管道中液体的第一流速；
所述超声换能器组包括交叉式配置的第一超声换能器对和第二超声换能器对，每个超声换能器对均包括第一超声波换能器和第二超声波换能器；
所述第一传播时间信息包括：在预先设定时间段内该超声换能器组中每个超声换能器对中的第一超声波换能器每次所发出的超声波传输到第二超声波换能器的时间以及第二超声波换能器每次所发出的超声波传输到第一超声波换能器的时间；
S2、获取设置在管道上的超声波发射组当前所对应的第二传播时间信息；
所述超声波发射组包括第一超声波传感器对和第二超声波传感器对，每个超声波传感器对均包括在管道上间隔第一距离设置的超声波发射器和超声波接收器；
所述第二传播时间信息包括：第一超声波传感器对中超声波发射器按管道液体顺流方向发射的超声波传输到对应的超声波接收器的时间、第二超声波传感器对中超声波发射器按管道液体逆流方向发射的超声波传输到对应的超声波接收器的时间；
S3、基于所述第一流速、所述第二传播时间信息以及管道中液体的温度，确定管道中液体的最终流速。
优选地，所述S1具体包括：
S11、基于管道上超声换能器组在预先设定时间段内所对应的第一传播时间信息，采用公式（1）获取该超声换能器组所对应的第一传播时间差值；
其中，所述公式（1）为：
Figure 331809DEST_PATH_IMAGE007
其中，
Figure 587210DEST_PATH_IMAGE008
为第一传播时间差值；
T nAB 为该超声换能器组中第一超声换能器对中第一超声波换能器在预先设定时间段内第n次发出的超声波传输到第一超声换能器对中对应的第二超声波换能器的时间；
T nBA 为该超声换能器组中第一超声换能器对中第二超声波换能器在预先设定时间段内第n次发出的超声波传输到第一超声换能器对中对应的第一超声波换能器的时间；
S nAB 为该超声换能器组中第二超声换能器对中第一超声波换能器在预先设定时间段内第n次发出的超声波传输到第二超声换能器对中对应的第二超声波换能器的时间；
S nBA 为该超声换能器组中第二超声换能器对中第二超声波换能器在预先设定时间段内第n次发出的超声波传输到第二超声换能器对中对应的第一超声波换能器的时间；
S12、基于管道上超声换能器组在预先设定时间段内所对应的第一传播时间信息，采用公式（2）获取超声换能器组所对应的第二数值；
所述公式（2）为：
Figure 460488DEST_PATH_IMAGE009
其中，K为第二数值；
S13、根据所述超声换能器组所对应的第一传播时间差值
Figure 342119DEST_PATH_IMAGE008
和所述第二数值K， 采用公式（3）获取管道中液体的第一流速；
所述公式（3）为：
Figure 450889DEST_PATH_IMAGE010
；
其中，V 1为管道中液体的第一流速；
Figure DEST_PATH_IMAGE012A
为管道中液体流向与第一直线的夹角；所述第一直线为第一超声换能器对中第 一超声波换能器和第二超声波换能器所在的直线；
L 1为第一超声换能器对或第二超声换能器对中第一超声波换能器到第二超声波换能器的距离。
优选地，所述S3包括：
S31、基于所述第一流速、所述第二传播时间信息以及管道中液体的温度，确定管道中液体的第二流速；
S32、基于所述第一流速和所述第二流速，确定管道中液体的最终流速。
优选地，所述S31具体包括：
S311、基于所述第二传播时间信息，采用公式（4）获取第二传播时间差；
所述公式（4）为：
Figure 124578DEST_PATH_IMAGE013
；
其中，
Figure 637599DEST_PATH_IMAGE014
为第二传播时间差；
t 1为第一超声波传感器对中超声波发射器按管道液体顺流方向发射的超声波传输到对应的超声波接收器的时间；
t 2为第二超声波传感器对中超声波发射器按管道液体逆流方向发射的超声波传输到对应的超声波接收器的时间；
S312、基于管道中液体的温度，确定超声波发射器发射的超声波在管道中液体的传播速度；
S313、基于所述第二传播时间差、超声波在管道中液体的传播速度、第一流速，采用公式（5）确定管道中液体的第二流速；
其中，所述公式（5）为：
Figure 567378DEST_PATH_IMAGE015
；
其中，V 2为管道中液体的第二流速；
c为超声波在管道中液体的传播速度；
L 2为第一距离。
优选地，所述S312具体包括：
基于所述管道中液体的温度，采用公式（6）确定超声波发射器发射的超声波在管道中液体的传播速度；
其中，所述公式（6）为：
c=331.5+0.607T；
其中，T为管道中液体的温度。
优选地，所述S32包括：
S321、基于所述第一流速和所述第二流速，采用公式（7）获取第三数值；
所述公式（7）为：
Figure 89626DEST_PATH_IMAGE016
；
其中，
Figure 522881DEST_PATH_IMAGE017
为第三数值；
S322、判断所述第三数值是否超过预先设定阈值，获取判断结果；
S323、基于所述判断结果、所述第一流速和/或所述第二流速，确定管道中液体的最终流速。
优选地，所述S323具体包括：
若所述判断结果为所述第三数值没有超过预先设定阈值，则根据所述第一流速和所述第二流速，采用公式（8）确定管道中液体的最终流速；
其中，所述公式（8）为：
Figure 65858DEST_PATH_IMAGE018
；
V 3为管道中液体的最终流速；
a为第一预先设定值，且a≠0；
b为第二预先设定值，且a+b=1。
优选地，所述S323具体包括：
若所述判断结果为所述第三数值超过了预先设定阈值，则将所述第一流速作为管道中液体的最终流速。
优选地，
所述超声换能器组中的第一超声换能器对和第二超声换能器对中的第一超声波换能器和第二超声波换能器均设置在管道的同一水平面中。
优选地，
所述超声波发射组与所述超声换能器组在管道上为接触设置。
有益效果
本发明的有益效果是：本发明的一种超声波测流方法，由于采用预先设定时间段内所对应的第一传播时间信息，并基于该第一传播时间信息获取管道中液体的第一流速，并进一步基于所述第一流速、所述第二传播时间信息以及管道中液体的温度，确定管道中液体的最终流速，相对于现有技术而言，在计算第一流速的过程中使用第一传播时间差值更为准确，使得管道中液体的最终流速确定地更为精确。
附图说明
图1为本发明的一种超声波测流方法流程图；
图2为本发明实施例中的管道上的超声换能器组结构示意图；
图3为本发明实施例中的超声波发射组结构示意图：
图4为现有技术中采用时差法测流时换能器示意图。
附图标记说明
A：第一超声波换能器；
B：第二超声波换能器；
C：超声波发射器；
D：超声波接收器；
1：第一换能器；
2：第二换能器。
具体实施方式
为了更好的解释本发明，以便于理解，下面结合附图，通过具体实施方式，对本发明作详细描述。
为了更好的理解上述技术方案，下面将参照附图更详细地描述本发明的示例性实施例。虽然附图中显示了本发明的示例性实施例，然而应当理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更清楚、透彻地理解本发明，并且能够将本发明的范围完整的传达给本领域的技术人员。
实施例一
参见图1，本实施例提供一种超声波测流方法，包括：
S1、获取设置在管道上的超声换能器组在预先设定时间段内所对应的第一传播时间信息，并基于该第一传播时间信息获取管道中液体的第一流速。
在具体应用中，换能器在管道的配置方式有3种，Z式是最常见的方式，即单声道，装置简单，适用于有足够长的直管段，流速分布为管道轴对称的场合；V式适用于流速不对称的流动流体的测量；当安装距离受到限制时，可采用X式。本实施例中，参见图2，所述超声换能器组包括交叉式（也就是X式）配置的第一超声换能器对和第二超声换能器对，每个超声换能器对均包括第一超声波换能器A和第二超声波换能器B；其中，每一对的超声换能器对中的第一超声波换能器A和第二超声波换能器B之间距离都为L 1。本实施例中所述超声换能器组中的第一超声换能器对和第二超声换能器对中的第一超声波换能器和第二超声波换能器均设置在管道的同一水平面中。
所述第一传播时间信息包括：在预先设定时间段内该超声换能器组中每个超声换能器对中的第一超声波换能器A每次所发出的超声波传输到第二超声波换能器B的时间以及第二超声波换能器B每次所发出的超声波传输到第一超声波换能器A的时间。
具体应用中，所述S1具体包括：
S11、基于管道上超声换能器组在预先设定时间段内所对应的第一传播时间信息，采用公式（1）获取该超声换能器组所对应的第一传播时间差值。
其中，所述公式（1）为：
Figure 922081DEST_PATH_IMAGE007
其中，
Figure 44758DEST_PATH_IMAGE008
为第一传播时间差值。
T nAB 为该超声换能器组中第一超声换能器对中第一超声波换能器A在预先设定时间段内第n次发出的超声波传输到第一超声换能器对中对应的第二超声波换能器B的时间。
T nBA 为该超声换能器组中第一超声换能器对中第二超声波换能器B在预先设定时间段内第n次发出的超声波传输到第一超声换能器对中对应的第一超声波换能器A的时间。
S nAB 为该超声换能器组中第二超声换能器对中第一超声波换能器A在预先设定时间段内第n次发出的超声波传输到第二超声换能器对中对应的第二超声波换能器B的时间。
S nBA 为该超声换能器组中第二超声换能器对中第二超声波换能器B在预先设定时间段内第n次发出的超声波传输到第二超声换能器对中对应的第一超声波换能器A的时间。
S12、基于管道上超声换能器组在预先设定时间段内所对应的第一传播时间信息，采用公式（2）获取超声换能器组所对应的第二数值。
所述公式（2）为：
Figure 598099DEST_PATH_IMAGE009
其中，K为第二数值。
S13、根据所述超声换能器组所对应的第一传播时间差值
Figure 46398DEST_PATH_IMAGE008
和所述第二数值K， 采用公式（3）获取管道中液体的第一流速。
所述公式（3）为：
Figure 560556DEST_PATH_IMAGE010
；
其中，V 1为管道中液体的第一流速。
Figure DEST_PATH_IMAGE012AA
为管道中液体流向与第一直线的夹角。
所述第一直线为第一超声换能器对中第一超声波换能器和第二超声波换能器所在的直线。
L 1为第一超声换能器对或第二超声换能器对中第一超声波换能器到第二超声波换能器的距离。
S2、获取设置在管道上的超声波发射组当前所对应的第二传播时间信息。
参见图3，本实施例中的所述超声波发射组包括第一超声波传感器对和第二超声波传感器对，每个超声波传感器对均包括在管道上间隔第一距离设置的超声波发射器C和超声波接收器D。
本实施例中所述超声波发射组与所述超声换能器组在管道上为接触设置。
所述第二传播时间信息包括：第一超声波传感器对中超声波发射器C按管道液体顺流方向发射的超声波传输到对应的超声波接收器D的时间、第二超声波传感器对中超声波发射器C按管道液体逆流方向发射的超声波传输到对应的超声波接收器D的时间。
具体的，假设超声波在静止的流体中的传播速度为
Figure DEST_PATH_IMAGE019
，流体的流速为u，则超声波 沿顺流和逆流的传播速度将不同。当第一超声波传感器对中超声波发射器C按顺流方向、第 二超声波传感器对中超声波发射器C按逆流方向发射超声波时，超声波到达超声波接收器D 所需要的时间t1和t2与流速之间的关系如下：
Figure 440918DEST_PATH_IMAGE020
；
Figure DEST_PATH_IMAGE021
。
S3、基于所述第一流速、所述第二传播时间信息以及管道中液体的温度，确定管道中液体的最终流速。
本实施例中所述S3包括：
S31、基于所述第一流速、所述第二传播时间信息以及管道中液体的温度，确定管道中液体的第二流速。
所述S31具体包括：
S311、基于所述第二传播时间信息，采用公式（4）获取第二传播时间差。
所述公式（4）为：
Figure 52028DEST_PATH_IMAGE013
；
其中，
Figure 608911DEST_PATH_IMAGE014
为第二传播时间差。
t 1为第一超声波传感器对中超声波发射器按管道液体顺流方向发射的超声波传输到对应的超声波接收器的时间。
t 2为第二超声波传感器对中超声波发射器按管道液体逆流方向发射的超声波传输到对应的超声波接收器的时间。
S312、基于管道中液体的温度，确定超声波发射器发射的超声波在管道中液体的传播速度，本实施例中管道的液体为水。
S313、基于所述第二传播时间差、超声波在管道中液体的传播速度、第一流速，采用公式（5）确定管道中液体的第二流速。
其中，所述公式（5）为：
Figure 970885DEST_PATH_IMAGE015
；
其中，V 2 为管道中液体的第二流速。
c为超声波在管道中液体的传播速度。本实施例中基于所述管道中液体的温度，采用公式（6）确定超声波发射器发射的超声波在管道中液体的传播速度c。
其中，所述公式（6）为：
c=331.5+0.607T；
其中，T为管道中液体的温度。
L 2为第一距离。
S32、基于所述第一流速和所述第二流速，确定管道中液体的最终流速。
所述S32包括：
S321、基于所述第一流速和所述第二流速，采用公式（7）获取第三数值。
所述公式（7）为：
Figure 638627DEST_PATH_IMAGE016
；
其中，
Figure 838664DEST_PATH_IMAGE017
为第三数值。
S322、判断所述第三数值是否超过预先设定阈值，获取判断结果。
S323、基于所述判断结果、所述第一流速和/或所述第二流速，确定管道中液体的最终流速。
本实施例中，所述S323具体包括：
若所述判断结果为所述第三数值没有超过预先设定阈值，则根据所述第一流速和所述第二流速，采用公式（8）确定管道中液体的最终流速。
其中，所述公式（8）为：
Figure 363186DEST_PATH_IMAGE018
；
V 3为管道中液体的最终流速。
a为第一预先设定值，且a≠0。
b为第二预先设定值，且a+b=1。
若所述判断结果为所述第三数值超过了预先设定阈值，则将所述第一流速作为管道中液体的最终流速。
实施例二
为更好的理解本发明实施例的方案，下面对本发明实施例的各步骤进行详细说明。
101、基于管道上超声换能器组在预先设定时间段内所对应的第一传播时间信息，采用公式（1）获取该超声换能器组所对应的第一传播时间差值。
其中，所述公式（1）为：
Figure 914253DEST_PATH_IMAGE007
其中，
Figure 244740DEST_PATH_IMAGE008
为第一传播时间差值。
T nAB 为该超声换能器组中第一超声换能器对中第一超声波换能器A在预先设定时间段内第n次发出的超声波传输到第一超声换能器对中对应的第二超声波换能器B的时间。
T nBA 为该超声换能器组中第一超声换能器对中第二超声波换能器B在预先设定时间段内第n次发出的超声波传输到第一超声换能器对中对应的第一超声波换能器A的时间。
S nAB 为该超声换能器组中第二超声换能器对中第一超声波换能器A在预先设定时间段内第n次发出的超声波传输到第二超声换能器对中对应的第二超声波换能器B的时间。
S nBA 为该超声换能器组中第二超声换能器对中第二超声波换能器B在预先设定时间段内第n次发出的超声波传输到第二超声换能器对中对应的第一超声波换能器A的时间。
在具体应用中，换能器在管道的配置方式有3种，Z式是最常见的方式，即单声道，装置简单，适用于有足够长的直管段，流速分布为管道轴对称的场合；V式适用于流速不对称的流动流体的测量；当安装距离受到限制时，可采用X式。本实施例中，参见图2，所述超声换能器组包括交叉式（也就是X式）配置的第一超声换能器对和第二超声换能器对，每个超声换能器对均包括第一超声波换能器A和第二超声波换能器B；其中，每一对的超声换能器对中的第一超声波换能器A和第二超声波换能器B之间距离都为L 1。
所述第一传播时间信息包括：在预先设定时间段内该超声换能器组中每个超声换能器对中的第一超声波换能器A每次所发出的超声波传输到第二超声波换能器B的时间以及第二超声波换能器B每次所发出的超声波传输到第一超声波换能器A的时间。
102、基于管道上超声换能器组在预先设定时间段内所对应的第一传播时间信息，采用公式（2）获取超声换能器组所对应的第二数值。
所述公式（2）为：
Figure 236967DEST_PATH_IMAGE009
其中，K为第二数值。
103、根据所述超声换能器组所对应的第一传播时间差值
Figure 932391DEST_PATH_IMAGE008
和所述第二数值K， 采用公式（3）获取管道中液体的第一流速。
所述公式（3）为：
Figure 194921DEST_PATH_IMAGE010
；
其中，V 1为管道中液体的第一流速。
Figure DEST_PATH_IMAGE012AAA
为管道中液体流向与第一直线的夹角。
所述第一直线为第一超声换能器对中第一超声波换能器和第二超声波换能器所在的直线。
L 1为第一超声换能器对或第二超声换能器对中第一超声波换能器到第二超声波换能器的距离。
本实施例中，参见图4，相对现有技术中仅根据一个从第一换能器1到第二换能器2的超声波的飞行时间T 12和从第二换能器2到第一换能器1的超声波的飞行时间T 21的时间差去计算流速v，在实际应用中难免会有误差，并不能保证这个从第一换能器1到第二换能器2的超声波的飞行时间T 12和从第二换能器2到第一换能器1的超声波的飞行时间T 21的时间差是准确，从而导致管道中液体流速v不能被准确得到。本实施例中的一种超声波测流方法，由于采用预先设定时间段内所对应的第一传播时间信息，并基于该第一传播时间信息获取管道中液体的第一流速，由于第一流速是根据公式（1）、公式（2）、公式（3）得到的，因此，本实施例中得到第一流速更为合理，因此进一步基于所述第一流速得到的管道中液体的最终流速更为精确。
104、获取设置在管道上的超声波发射组当前所对应的第二传播时间信息。
参见图3，本实施例中的所述超声波发射组包括第一超声波传感器对和第二超声波传感器对，每个超声波传感器对均包括在管道上间隔第一距离设置的超声波发射器C和超声波接收器D；本实施例中第一超声波传感器对和第二超声波传感器对是对齐设置的。
所述第二传播时间信息包括：第一超声波传感器对中超声波发射器C按管道液体顺流方向发射的超声波传输到对应的超声波接收器D的时间、第二超声波传感器对中超声波发射器C按管道液体逆流方向发射的超声波传输到对应的超声波接收器D的时间。
实际上，假设超声波在静止的流体中的传播速度为
Figure 391416DEST_PATH_IMAGE019
，流体的流速为u，则超声波 沿顺流和逆流的传播速度将不同。当第一超声波传感器对中超声波发射器C按顺流方向、第 二超声波传感器对中超声波发射器C按逆流方向发射超声波时，超声波到达超声波接收器D 所需要的时间t 1和t 2与流速之间的关系为：
Figure 503728DEST_PATH_IMAGE020
；
Figure 697949DEST_PATH_IMAGE021
。
105、基于所述第二传播时间信息，采用公式（4）获取第二传播时间差。
所述公式（4）为：
Figure 895712DEST_PATH_IMAGE013
；
其中，
Figure 69467DEST_PATH_IMAGE014
为第二传播时间差。
t 1为第一超声波传感器对中超声波发射器按管道液体顺流方向发射的超声波传输到对应的超声波接收器的时间。
t 2为第二超声波传感器对中超声波发射器按管道液体逆流方向发射的超声波传输到对应的超声波接收器的时间。
106、基于管道中液体的温度，确定超声波发射器发射的超声波在管道中液体的传播速度。
107、基于所述第二传播时间差、超声波在管道中液体的传播速度、第一流速，采用公式（5）确定管道中液体的第二流速。
其中，所述公式（5）为：
Figure 36286DEST_PATH_IMAGE015
；
其中，V 2为管道中液体的第二流速。
c为超声波在管道中液体的传播速度。
L 2为第一距离。
108、基于所述第一流速和所述第二流速，采用公式（7）获取第三数值。
所述公式（7）为：
Figure 135829DEST_PATH_IMAGE016
；
其中，
Figure 820888DEST_PATH_IMAGE017
为第三数值。
109、判断所述第三数值是否超过预先设定阈值，获取判断结果。
若所述判断结果为所述第三数值没有超过预先设定阈值，则根据所述第一流速和所述第二流速，采用公式（8）确定管道中液体的最终流速。
其中，所述公式（8）为：
Figure 296869DEST_PATH_IMAGE018
；
V 3为管道中液体的最终流速。
a为第一预先设定值，且a≠0。
b为第二预先设定值，且a+b=1。
若所述判断结果为所述第三数值超过了预先设定阈值，则将所述第一流速作为管道中液体的最终流速。
本实施例中的一种超声波测流方法，由于采用预先设定时间段内所对应的第一传播时间信息，并基于该第一传播时间信息获取管道中液体的第一流速，并进一步基于所述第一流速、所述第二传播时间信息以及管道中液体的温度，确定管道中液体的最终流速，相对于现有技术而言，在计算第一流速的过程中使用第一传播时间差值更为准确，同时，也考虑了液体温度因素使得管道中液体的最终流速更为精确。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例，或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质（包括但不限于磁盘存储器、CD-ROM、光学存储器等）上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备（系统）和计算机程序产品的流程图和／或方框图来描述的。应理解可由计算机程序指令实现流程图和／或方框图中的每一流程和／或方框、以及流程图和／或方框图中的流程和／或方框的结合。
应当注意的是，在权利要求中，不应将位于括号之间的任何附图标记理解成对权利要求的限制。词语“包含”不排除存在未列在权利要求中的部件或步骤。位于部件之前的词语“一”或“一个”不排除存在多个这样的部件。本发明可以借助于包括有若干不同部件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的权利要求中，这些装置中的若干个可以是通过同一个硬件来具体体现。词语第一、第二、第三等的使用，仅是为了表述方便，而不表示任何顺序。可将这些词语理解为部件名称的一部分。
此外，需要说明的是，在本说明书的描述中，术语“一个实施例”、“一些实施例”、“实施例”、“示例”、“具体示例”或“一些示例”等的描述，是指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管已描述了本发明的优选实施例，但本领域的技术人员在得知了基本创造性概念后，则可对这些实施例作出另外的变更和修改。所以，权利要求应该解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然，本领域的技术人员可以对本发明进行各种修改和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也应该包含这些修改和变型在内。

Claims (10)

1.一种超声波测流方法，其特征在于，包括：
S1、获取设置在管道上的超声换能器组在预先设定时间段内所对应的第一传播时间信息，并基于该第一传播时间信息获取管道中液体的第一流速；
所述超声换能器组包括交叉式配置的第一超声换能器对和第二超声换能器对，每个超声换能器对均包括第一超声波换能器和第二超声波换能器；
所述第一传播时间信息包括：在预先设定时间段内该超声换能器组中每个超声换能器对中的第一超声波换能器每次所发出的超声波传输到第二超声波换能器的时间以及第二超声波换能器每次所发出的超声波传输到第一超声波换能器的时间；
S2、获取设置在管道上的超声波发射组当前所对应的第二传播时间信息；
所述超声波发射组包括第一超声波传感器对和第二超声波传感器对，每个超声波传感器对均包括在管道上间隔第一距离设置的超声波发射器和超声波接收器；
所述第二传播时间信息包括：第一超声波传感器对中超声波发射器按管道液体顺流方向发射的超声波传输到对应的超声波接收器的时间、第二超声波传感器对中超声波发射器按管道液体逆流方向发射的超声波传输到对应的超声波接收器的时间；
S3、基于所述第一流速、所述第二传播时间信息以及管道中液体的温度，确定管道中液体的最终流速。
2.根据权利要求1所述的超声波测流方法，其特征在于，所述S1具体包括：
S11、基于管道上超声换能器组在预先设定时间段内所对应的第一传播时间信息，采用公式（1）获取该超声换能器组所对应的第一传播时间差值；
其中，所述公式（1）为：
Figure 651071DEST_PATH_IMAGE001
其中，
Figure 408812DEST_PATH_IMAGE002
为第一传播时间差值；
T nAB 为该超声换能器组中第一超声换能器对中第一超声波换能器在预先设定时间段内第n次发出的超声波传输到第一超声换能器对中对应的第二超声波换能器的时间；
T nBA 为该超声换能器组中第一超声换能器对中第二超声波换能器在预先设定时间段内第n次发出的超声波传输到第一超声换能器对中对应的第一超声波换能器的时间；
S nAB 为该超声换能器组中第二超声换能器对中第一超声波换能器在预先设定时间段内第n次发出的超声波传输到第二超声换能器对中对应的第二超声波换能器的时间；
S nBA 为该超声换能器组中第二超声换能器对中第二超声波换能器在预先设定时间段内第n次发出的超声波传输到第二超声换能器对中对应的第一超声波换能器的时间；
S12、基于管道上超声换能器组在预先设定时间段内所对应的第一传播时间信息，采用公式（2）获取超声换能器组所对应的第二数值；
所述公式（2）为：
Figure 384858DEST_PATH_IMAGE003
其中，K为第二数值；
S13、根据所述超声换能器组所对应的第一传播时间差值
Figure 987878DEST_PATH_IMAGE002
和所述第二数值K，采用 公式（3）获取管道中液体的第一流速；
所述公式（3）为：
Figure 365770DEST_PATH_IMAGE004
；
其中，V 1为管道中液体的第一流速；
Figure 530297DEST_PATH_IMAGE005
为管道中液体流向与第一直线的夹角；所述第一直线为第一超声换能器对中第一超 声波换能器和第二超声波换能器所在的直线；
L 1为第一超声换能器对或第二超声换能器对中第一超声波换能器到第二超声波换能器的距离。
3.根据权利要求2所述的超声波测流方法，其特征在于，所述S3包括：
S31、基于所述第一流速、所述第二传播时间信息以及管道中液体的温度，确定管道中液体的第二流速；
S32、基于所述第一流速和所述第二流速，确定管道中液体的最终流速。
4.根据权利要求3所述的超声波测流方法，其特征在于，所述S31具体包括：
S311、基于所述第二传播时间信息，采用公式（4）获取第二传播时间差；
所述公式（4）为：
Figure 728060DEST_PATH_IMAGE006
；
其中，
Figure 400350DEST_PATH_IMAGE007
为第二传播时间差；
t 1为第一超声波传感器对中超声波发射器按管道液体顺流方向发射的超声波传输到对应的超声波接收器的时间；
t 2为第二超声波传感器对中超声波发射器按管道液体逆流方向发射的超声波传输到对应的超声波接收器的时间；
S312、基于管道中液体的温度，确定超声波发射器发射的超声波在管道中液体的传播速度；
S313、基于所述第二传播时间差、超声波在管道中液体的传播速度、第一流速，采用公式（5）确定管道中液体的第二流速；
其中，所述公式（5）为：
Figure 101590DEST_PATH_IMAGE008
；
其中，V 2 为管道中液体的第二流速；
c为超声波在管道中液体的传播速度；
L 2 为第一距离。
5.根据权利要求4所述的超声波测流方法，其特征在于，所述S312具体包括：
基于所述管道中液体的温度，采用公式（6）确定超声波发射器发射的超声波在管道中液体的传播速度；
其中，所述公式（6）为：
c=331.5+0.607T；
其中，T为管道中液体的温度。
6.根据权利要求5所述的超声波测流方法，其特征在于，所述S32包括：
S321、基于所述第一流速和所述第二流速，采用公式（7）获取第三数值；
所述公式（7）为：
Figure 466712DEST_PATH_IMAGE009
；
其中，
Figure 886192DEST_PATH_IMAGE010
为第三数值；
S322、判断所述第三数值是否超过预先设定阈值，获取判断结果；
S323、基于所述判断结果、所述第一流速和/或所述第二流速，确定管道中液体的最终流速。
7.根据权利要求6所述的超声波测流方法，其特征在于，所述S323具体包括：
若所述判断结果为所述第三数值没有超过预先设定阈值，则根据所述第一流速和所述第二流速，采用公式（8）确定管道中液体的最终流速；
其中，所述公式（8）为：
Figure 362173DEST_PATH_IMAGE011
；
V 3为管道中液体的最终流速；
a为第一预先设定值，且a≠0；
b为第二预先设定值，且a+b=1。
8.根据权利要求6所述的超声波测流方法，其特征在于，所述S323具体包括：
若所述判断结果为所述第三数值超过了预先设定阈值，则将所述第一流速作为管道中液体的最终流速。
9.根据权利要求7或8所述的超声波测流方法，其特征在于，
所述超声换能器组中的第一超声换能器对和第二超声换能器对中的第一超声波换能器和第二超声波换能器均设置在管道的同一水平面中。
10.根据权利要求9所述的超声波测流方法，其特征在于，
所述超声波发射组与所述超声换能器组在管道上为接触设置。
CN202210719013.9A 2022-06-23 2022-06-23 一种超声波测流方法 Active CN114814285B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210719013.9A CN114814285B (zh) 2022-06-23 2022-06-23 一种超声波测流方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210719013.9A CN114814285B (zh) 2022-06-23 2022-06-23 一种超声波测流方法

Publications (2)

Publication Number Publication Date
CN114814285A true CN114814285A (zh) 2022-07-29
CN114814285B CN114814285B (zh) 2022-09-09

Family

ID=82521344

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210719013.9A Active CN114814285B (zh) 2022-06-23 2022-06-23 一种超声波测流方法

Country Status (1)

Country Link
CN (1) CN114814285B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117368518A (zh) * 2023-12-04 2024-01-09 湖北工业大学 一种管道内超声波飞行时间的测量系统及方法

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1248703A (zh) * 1998-09-03 2000-03-29 昌民技术株式会社 超声流速测量方法
EP1719980A1 (en) * 2004-02-26 2006-11-08 Fuji Electric Systems Co., Ltd. Ultrasonic flowmeter and ultrasonic flow rate measurement method
JP2010101767A (ja) * 2008-10-24 2010-05-06 Yokogawa Electric Corp 超音波流量計
WO2011010097A1 (en) * 2009-07-24 2011-01-27 Wayne Rudd A method for determining flow characteristics of a medium and associated apparatus
US20130238260A1 (en) * 2012-03-08 2013-09-12 Yokogawa Electric Corporation Ultrasonic measuring device
CN103868555A (zh) * 2012-12-11 2014-06-18 南京理工大学 一种用于超声波流量计的循环时差检测方法
CN105698884A (zh) * 2016-03-07 2016-06-22 上海电气自动化设计研究所有限公司 一种时差式超声波流量计改进的测量方法
CN111487437A (zh) * 2020-04-20 2020-08-04 东南大学 一种声学法测量烟道内烟气流速的装置及方法
US20200355532A1 (en) * 2019-03-29 2020-11-12 Krohne Ag Ultrasonic flowmeter, method for operating an ultrasonic flowmeter, measuring system and method for operating a measuring system
CN112665670A (zh) * 2021-01-21 2021-04-16 沈阳佳德联益能源科技股份有限公司 一种NB-IoT物联网超声波水表
CN112964317A (zh) * 2021-02-07 2021-06-15 山东锋士信息技术有限公司 一种基于温度补偿的流量测算方法
CN113532560A (zh) * 2021-07-27 2021-10-22 北京奥特美克科技股份有限公司 一种流量测量方法及装置

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1248703A (zh) * 1998-09-03 2000-03-29 昌民技术株式会社 超声流速测量方法
EP1719980A1 (en) * 2004-02-26 2006-11-08 Fuji Electric Systems Co., Ltd. Ultrasonic flowmeter and ultrasonic flow rate measurement method
JP2010101767A (ja) * 2008-10-24 2010-05-06 Yokogawa Electric Corp 超音波流量計
WO2011010097A1 (en) * 2009-07-24 2011-01-27 Wayne Rudd A method for determining flow characteristics of a medium and associated apparatus
US20130238260A1 (en) * 2012-03-08 2013-09-12 Yokogawa Electric Corporation Ultrasonic measuring device
CN103868555A (zh) * 2012-12-11 2014-06-18 南京理工大学 一种用于超声波流量计的循环时差检测方法
CN105698884A (zh) * 2016-03-07 2016-06-22 上海电气自动化设计研究所有限公司 一种时差式超声波流量计改进的测量方法
US20200355532A1 (en) * 2019-03-29 2020-11-12 Krohne Ag Ultrasonic flowmeter, method for operating an ultrasonic flowmeter, measuring system and method for operating a measuring system
CN111487437A (zh) * 2020-04-20 2020-08-04 东南大学 一种声学法测量烟道内烟气流速的装置及方法
CN112665670A (zh) * 2021-01-21 2021-04-16 沈阳佳德联益能源科技股份有限公司 一种NB-IoT物联网超声波水表
CN112964317A (zh) * 2021-02-07 2021-06-15 山东锋士信息技术有限公司 一种基于温度补偿的流量测算方法
CN113532560A (zh) * 2021-07-27 2021-10-22 北京奥特美克科技股份有限公司 一种流量测量方法及装置

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
MINGWEI L等: "《Research on improving the accuracy of the ultrasonic flow-meter with time difference method》", 《2010 INTERNATIONAL CONFERENCE ON ELECTRICAL AND CONTROL ENGINEERING.IEEE》 *
党瑞荣等: "《基于X型的井下流量测量误差校正系统及计算方法》", 《仪表技术》 *
罗永: "《提升时差法超声流量计计量精度关键技术研究》", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 *
黄军垒等: "《时差法超声波热能流量表的高精度补偿算法研究》", 《轻工学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117368518A (zh) * 2023-12-04 2024-01-09 湖北工业大学 一种管道内超声波飞行时间的测量系统及方法
CN117368518B (zh) * 2023-12-04 2024-02-23 湖北工业大学 一种管道内超声波飞行时间的测量系统及方法

Also Published As

Publication number Publication date
CN114814285B (zh) 2022-09-09

Similar Documents

Publication Publication Date Title
CN106643939B (zh) 用于超声波流量计计算超声波传播时间的方法
US5835884A (en) Method of determining a characteristic of a fluid
EP2687828A1 (en) Ultrasonic wedge and method for determining the speed of sound in same
RU2590318C2 (ru) Проверка температуры ультразвуковых расходомеров
CN114814285B (zh) 一种超声波测流方法
CN103196504B (zh) 一种多声道超声波流量测量方法及装置
CN105157771A (zh) 一种时差式超声波流量测量方法及装置
CN106768103A (zh) 一种超声波流量计自动校准时间偏差的方法
CN104457871A (zh) 一种流量计及流体测量方法
Tian et al. Energy peak fitting of echo based signal processing method for ultrasonic gas flow meter
CN104236648A (zh) 超声波流量计
US6816808B2 (en) Peak switch detector for transit time ultrasonic meters
US3788140A (en) Electroacoustical flow metering apparatus
WO2018086086A1 (zh) 一种提高检测氧气浓度准确性的方法
CN101922954A (zh) 超声波测量信号的包络线处理方法
CN206291930U (zh) 一种超声波质量流量计
Buchholz et al. Bubble column bioreactors: I. Methods for measuring the bubble size
CN102095889A (zh) 三通道超声时差流速测量方法
US20230243682A1 (en) Ultrasonic flow measurement
CN102023038A (zh) 一种管道流量的超声波测量方法
Waluś Mathematical modelling of an ultrasonic flowmeter primary device
Li et al. Peak ratio characteristic value sequence based signal processing method for transit-time ultrasonic gas flowmeter
CN106768104A (zh) 一种超声波质量流量计
JPS59173715A (ja) 相関型流速計
CN111337092B (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