CN102103147B - 超声自相关横向流速测量方法 - Google Patents
超声自相关横向流速测量方法 Download PDFInfo
- Publication number
- CN102103147B CN102103147B CN2011100040999A CN201110004099A CN102103147B CN 102103147 B CN102103147 B CN 102103147B CN 2011100040999 A CN2011100040999 A CN 2011100040999A CN 201110004099 A CN201110004099 A CN 201110004099A CN 102103147 B CN102103147 B CN 102103147B
- Authority
- CN
- China
- Prior art keywords
- ultrasonic
- particle
- focus
- flow velocity
- transverse flow
- 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
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种超声自相关横向流速测量方法,它利用单探测声束直接准确地测量与探测超声发射方向的垂直方向上的流速分量,即横向流速。这种利用超声直接测量横向流速的方法,不用知道超声束和粒子移动方向的夹角,当有粒子穿过超声聚焦区域时,产生较强的超声反射信号,由超声反射信号的归一化自相关函数,得到粒子穿过超声聚焦区域的渡越时间,用超声聚焦区域的横向宽度除以渡越时间,即是横向流速。
Description
技术领域
本发明涉及超声检测流速技术领域,特别涉及一种超声自相关横向流速测量方法,此处流速指的是含固态粒子的流体的流速。
背景技术
流速测量在医学和工业领域有广泛的应用,如水流监测,血流测量等。目前常用的测量方法是多普勒超声技术,这种方法是根据多普勒效应,通过计算流动粒子对入射超声的多普勒频移,而计算出粒子的移动速度。超声多普勒频移正比于超声束和粒子流速方向的夹角的余弦值,因此,多普勒超声技术只能给出流速在探测声束方向的流速分量,即纵向流速,而无法给出和探测声束垂直的流速分量,即横向流速。为了解决此问题,余建国等提出双超声束多普勒血流速度测量方法(中国专利公开号:CN1257695),叶国帆提出双向定量测量血流绝对速度的超声多普勒方法与仪器(中国专利公开号:CN85100528),这两种方法都是利用多普勒效应及双探测超声束,进行绝对速度测量。但是当流速方向和探测声束方向接近垂直时,多普勒角接近九十度,因此引起的多普勒频移很小,会导致较大的检测误差。
发明内容
本发明目的在于提出一种超声自相关横向流速测量方法,它利用单探测声束直接准确地测量与探测超声发射方向的垂直方向上的流速分量。
本发明所提出的技术解决方案是这样的:一种超声自相关横向流速测量方法,这里的流速指的是含固态粒子的流体(6)的流速,包括有:第1超声聚焦换能器(1)、第2超声聚焦换能器(2)、信号采集及处理系统(5)、归一化自相关运算,第1超声聚焦换能器(1)和第2超声聚焦换能器器(2)聚焦位置相同,第1超声聚焦换能器(1)发射聚焦超声(3),当有粒子穿过超声聚焦区域时,产生较强的超声反射信号(4),由第2超声聚焦换能器(2)接受,由信号采集及处理系统(5)进行信号采集并计算超声反射信号f(t)的归一化自相关函数,得到沿探测超声发射方向的垂直方向上的流速分量,即横向流速,由粒子反射的超声反射信号f(t)得到粒子穿过超声聚焦区域的渡越时间τ0,则横向流速为V=L/τ0,这里L是超声聚焦区域的横向宽度,参见图1所示。
粒子穿过超声聚焦区域的渡越时间τ0通过下列方式获得:由粒子反射的超声反射信号f(t)的归一化自相关函数计算粒子穿过超声聚焦区域的渡越时间τ0,粒子超声反射信号f(t)的归一化自相关函数的斜率的倒数等于粒子穿过超声聚焦区域的渡越时间τ0,如公式(1)所示,
在公式(1),f(t)为测量到的超声反射信号,τ0为粒子穿过超声聚焦区域的渡越时间,τ为自相关计算中的时间延迟。则横向流速为V=L/τ0,这里L是超声聚焦区域的横向宽度。
本发明的原理如图2所示,在图2中,弧线(7)和弧线(8)之间的区域表示放大的第1超声聚焦换能器(1)和第2超声聚焦换能器(2)两者聚焦区域重合之处,即探测区域(9),探测区域(9)的腰宽L大于粒子尺寸。当没有粒子流过探测区域(9)时,由探测区域(9)反射回去的超声很弱,当有粒子穿过探测区域(9)时,会有较强的超声反射,因此第2超声聚焦换能器(2)会接受到一个较强的脉冲信号,当粒子移动速度较慢时,粒子穿过探测区域(9)的时间较长,产生的脉冲信号的脉宽较大,当粒子移动速度较快时,粒子穿过探测区域(9)的时间较短,产生的脉冲信号的脉宽较窄,脉冲信号的脉宽等于粒子穿过探测区域(9)的时间,即渡越时间τ0。超声聚焦换能器器(2)也可以接收到在探测区域(9)之外的粒子反射的超声信号,但是由于处于探测区域(9)之外的粒子反射的超声信号远小于探测区域(9)内部粒子反射的超声信号,因此,探测区域(9)之外的粒子反射的超声信号可以忽略,所以,脉冲信号的宽度就等于粒子经过探测区域(9)的渡越时间τ0。如公式(1)所示,超声反射信号的归一化自相关函数的斜率的倒数等于脉冲信号的脉宽。因此如果知道探测区域(9)的宽度L,如图(2)所示,那么就可以知道粒子穿过探测区域(9)的速度,即流速V=L/τ0。
当粒子流动方向和探测区域(9)垂直时,如图3(a)所示,箭头表示粒子的流动方向,L是探测区域(9)的横向宽度,则这时测量的速度V=L/τ0,和探测声束垂直,即测量的是和探测声束垂直的流速,即横向流速。当粒子流动方向和探测声束不垂直时,如图3(b)所示,箭头表示粒子的流动方向,粒子在探测区域(9)中实际经过的距离为L1,因此粒子的绝对速度为V=L1/τ0,那么粒子的横向速度为cos(θ)L1/τ0,而cos(θ)L1/τ0等于L/τ0。因此,不论粒子的运动方向是否和探测声束垂直,由公示V=L/τ0得到的就是粒子的横向速度,L是超声探测区域(9)的横向宽度,这就是本方法的基本原理。
与现有技术相比,本发明具有如下显著效果:
由于本发明不是基于多普勒效应,因此避免了目前的基于多普勒效应的流速测量方法的缺陷,可以有效的测量横向流速。基于多普勒效应的流速测量方法适合于测量纵向流速,而本发明的方法适合于测量横向流速,二者结合,就可以给出粒子的绝对流速,而不用知道探测声束和粒子移动方向的夹角。
本发明在医学和工业领域有广泛的应用,可用于进行含固态粒子的流体的流速测量,如水流监测、烟尘气流速度检测、血流测量等。
附图说明
图1是本发明超声自相关横向流速测量方法示意图。
图2是图1所示测量方法的原理图,图中表示放大的探测区域,显示了粒子流过探测区域时的情况。
图3是本发明的原理图,在放大的探测区域中,显示了粒子垂直流过探测区域图3(a)和非垂直流过探测区域图3(b)时的情况。
图4是本发明一个实施例的示意图。
具体实施方式
通过下面实施例对本发明作进一步详细阐述。
参见图4所示,待测样品(12)为待测的含固态粒子的水流,计算机(11)通过控制超声控制器(10)来控制第1超声聚焦换能器(1)发射聚焦的脉冲超声(3),粒子反射的超声(4)由第2超声聚焦换能器(2)接收,并通过信号采集及处理系统(5)传给计算机(11),设采集到的信号为f(t),用公式(1)计算信号f(t)的归一化自相关函数,f(t)的归一化自相关函数的斜率的倒数等于粒子穿过超声探测区域的渡越时间τ0,则横向流速为V=L/τ0,这里L是超声探测区域(9)宽度L。
探测区域(9)宽度L由以下方法测量:用一个已知流速为V1的粒子流作为样品,首先,调节探测声束和粒子流动方向接近垂直,这个由多普勒频移进行判断,调节探测声束的方向,当多普勒频移接近零时,这时粒子的运动方向和探测声束接近垂直。则按照上述同样的方法,测量得到粒子穿过探测区域(9)的时间τ1,则可以得到探测区域(9)的宽度L=v1τ1。
因此,待测样品(12)的横向流速为v1τ1/τ0。
这里请注意,探测区域(9)宽度L是系统参数,用一个流速已知的样品进行一次实验即可确定。
Claims (2)
1.一种超声自相关横向流速测量方法,用第1超声聚焦换能器(1)发射聚焦超声(3),当有粒子穿过超声聚焦区域时,产生较强的超声反射信号(4),由第2超声聚焦换能器(2)接收,第1超声聚焦换能器(1)和第2超声聚焦换能器(2)聚焦位置相同,由信号采集及处理系统(5)进行信号采集并计算超声反射信号f(t)的归一化自相关函数,得到沿探测超声发射方向的垂直方向上的流速分量,即横向流速,其特征在于:由粒子反射的超声反射信号f(t)得到粒子穿过超声聚焦区域的渡越时间τ0,则横向流速为V=L/τ0,这里L是超声聚焦区域的横向宽度。
2.根据权利要求1所述的超声自相关横向流速测量方法,其特征在于:所述粒子穿过超声聚焦区域的渡越时间τ0通过下列方式获得:由粒子反射的超声反射信号f(t)的归一化自相关函数计算粒子穿过超声聚焦区域的渡越时间τ0,粒子超声反射信号f(t)的归一化自相关函数的斜率的倒数等于粒子穿过超声聚焦区域的渡越时间τ0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100040999A CN102103147B (zh) | 2011-01-05 | 2011-01-05 | 超声自相关横向流速测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100040999A CN102103147B (zh) | 2011-01-05 | 2011-01-05 | 超声自相关横向流速测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102103147A CN102103147A (zh) | 2011-06-22 |
CN102103147B true CN102103147B (zh) | 2012-04-18 |
Family
ID=44156076
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011100040999A Expired - Fee Related CN102103147B (zh) | 2011-01-05 | 2011-01-05 | 超声自相关横向流速测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102103147B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102507973A (zh) * | 2011-09-29 | 2012-06-20 | 河南安宏信息科技有限公司 | 一种检测水流的方法及装置 |
US10646197B2 (en) * | 2016-07-06 | 2020-05-12 | Biosense Webster (Israel) Ltd. | Ascertaining tissue thickness |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE2813754A1 (de) * | 1977-03-30 | 1978-10-05 | Nat Res Dev | Verfahren und vorrichtung zum erkennen der stroemung eines materials durch wellenenergiestrahlung |
CN1248703A (zh) * | 1998-09-03 | 2000-03-29 | 昌民技术株式会社 | 超声流速测量方法 |
CN1320809A (zh) * | 2000-04-24 | 2001-11-07 | 昌民技术有限公司 | 超声波流速测量仪 |
CN101162164A (zh) * | 2007-11-16 | 2008-04-16 | 浙江理工大学 | 时差法超声波流量计的频率调制波形标记方法 |
-
2011
- 2011-01-05 CN CN2011100040999A patent/CN102103147B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE2813754A1 (de) * | 1977-03-30 | 1978-10-05 | Nat Res Dev | Verfahren und vorrichtung zum erkennen der stroemung eines materials durch wellenenergiestrahlung |
CN1248703A (zh) * | 1998-09-03 | 2000-03-29 | 昌民技术株式会社 | 超声流速测量方法 |
CN1320809A (zh) * | 2000-04-24 | 2001-11-07 | 昌民技术有限公司 | 超声波流速测量仪 |
CN101162164A (zh) * | 2007-11-16 | 2008-04-16 | 浙江理工大学 | 时差法超声波流量计的频率调制波形标记方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102103147A (zh) | 2011-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101266228B (zh) | 一种材料声速测量方法 | |
CN101004353A (zh) | 一种用于超声波流量计时差交叉检测方法 | |
JP2019502119A (ja) | 改良型ビーム整形音響信号伝搬時間差式流量計 | |
CN1997911A (zh) | 使用超声波的距离测量方法与装置 | |
CN102261937A (zh) | 一种高精度时差式超声波流量计及其流量测量方法 | |
CN105699486B (zh) | 一种斜角表面裂纹倾斜角度和深度的检测方法 | |
Andruszkiewicz et al. | Gas bubble detection in liquid metals by means of the ultrasound transit-time-technique | |
CN100380101C (zh) | 多普勒型超声波流量计 | |
NGUYEN et al. | Development of multiwave method using ultrasonic pulse Doppler method for measuring two-phase flow | |
CN101828929A (zh) | 利用视在位移的多普勒血流速度矢量测量方法 | |
CN106153132A (zh) | 基于Lamb波的非接触式流体流量测量系统及方法 | |
CN102103147B (zh) | 超声自相关横向流速测量方法 | |
Han et al. | Studies on the transducers of clamp-on transit-time ultrasonic flow meter | |
CN102508249B (zh) | 一种基于dsp的高精度超声波测距系统及其测距方法 | |
CN102829830B (zh) | 用于超声波流量检测中检测超声波传播速度的方法及装置 | |
CN105467395B (zh) | 一种超远距离超声波测量仪 | |
CN103932737A (zh) | 一种心血管血液流速传感器 | |
US5325726A (en) | Method and device for measuring volume flows in liquids and gases | |
US8714029B2 (en) | Flow measuring device and method including both a doppler frequency shift measurement method and travel time measurement method | |
CN116295149A (zh) | 一种基于时差式超声波流量计的管道气泡大小测量系统 | |
CN103217196A (zh) | 超声波流量传感器 | |
Kou et al. | Research of long range accurate ranging technology based on ultrasonic sensor measurement | |
CN105300654A (zh) | 一种低速风洞内的超声精确测距系统 | |
Mahadeva et al. | Studies of the accuracy of clamp-on transit time ultrasonic flowmeters | |
CN203177907U (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: 20120418 Termination date: 20160105 |
|
EXPY | Termination of patent right or utility model |