CN105092430B - 一种基于发散超声波衰减的颗粒粒度测量装置及方法 - Google Patents
一种基于发散超声波衰减的颗粒粒度测量装置及方法 Download PDFInfo
- Publication number
- CN105092430B CN105092430B CN201410215801.XA CN201410215801A CN105092430B CN 105092430 B CN105092430 B CN 105092430B CN 201410215801 A CN201410215801 A CN 201410215801A CN 105092430 B CN105092430 B CN 105092430B
- Authority
- CN
- China
- Prior art keywords
- ultrasonic wave
- ultrasonic
- signal
- measurement
- transducer
- 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.)
- Active
Links
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开一种基于发散声衰减测量两相体系总颗粒粒度的测量装置,包括超声波信号发生器、前置信号处理电路、超声波发射聚焦换能器、超声波接收换能器、信号调理电路、数据采集器和中心计算机,所述超声波发射聚焦换能器和所述超声波接收换能器之间为测量区。本发明测量***结构简单、价廉,具有能提高测量的准确度,能适用在不同的测量条件下。
Description
技术领域
本发明涉及超声检测技术,特别涉及一种利用发散超声波衰减对两相流中颗粒粒度进行测量的装置及方法。
背景技术
颗粒是指在一尺寸范围内处于离散状态、具有特定形状的几何体,通常以固体形式存在,也可以是液体或气泡。电厂锅炉中的煤粉、大气中的粉尘、河流泥沙、血液中的红血球、牛奶中的蛋白质均是以颗粒状存在。两相或多相体系中颗粒的粒度大小对现代工业生产如能源、环境、材料、生物、化工等领域有极其重要的意义。生产过程中具有合适粒度的颗粒不仅可以提高生产效率,保证产品质量,还可以节约能源,减少污染排放。
超声工业检测技术是利用超声波为信息载体来进行各种检测和测量的技术。超声工业检测通常分超声无损探伤和超声工业测量两个方面。超声工业测量,是利用待测的与介质特性和状态有关的非声学量来描述介质声学特性的超声量之间存在的关系,通过这些超声量的测定来分析介质特性、评价介质质量、并测出某些与工程有关的参量。超声波在两相系中的传播规律及衰减大小与颗粒物的粒度密切相关,所以可用作颗粒粒度的测量。
相比于采用其它原理的颗粒测量方法如电感应法、图像法、光散射法等测量方法,超声波具有强的穿透力,可在有色甚至不透明的物质中传播并具有测量速度快,容易实现测量和数据的自动化等优点,超声波传感器价格低且耐污损,测量***简单方便。目前常规的主动式超声波测量颗粒浓度方法是利用声衰减原理进行测量,在这种测量方法中,声波被假设为垂直入射的平面波,亦即入射声波抵达颗粒两相区域前没有任何声扩散效应,或者声扩散效应可以忽略,声衰减方法中接收探头在测量中必须要收到衰减后良好的信号,才能保证测量的准确度。这一假设使得以往测量方法中采用的经典声学模型得以应用,例如Epstein和Carhart,Allegra和Hawley先后发展起来的ECAH模型就可以较好地用来描述平面入射波的声波动。这种方法目前也已经得到较多的应用,但是在实际应用中还会面临不少问题,常用的超声波传感器发出声波一般采用活塞波来描述,由于传感器的尺寸有限,理论上不可能产生完全的平面波,在远场往往近似为球面声波。这样,就造成了实际波动和理论假设的差异,并且不可避免给测量带来误差。
发明内容
本发明所要解决的技术问题是针对上述背景技术的不足,提供了一种基于发散超声波信号测量颗粒粒度的测量装置及方法,入射超声波由于受到颗粒介质的散射和吸收,透射声波强度会衰减,为了正确描述颗粒体系中发散声的散射和衰减,通过Monte-carlo方法建立一个基于发散声波的理论模型,求得颗粒的粒度。
为解决上述技术问题,本发明一方面提供一种基于发散声衰减测量两相体系总颗粒粒度的测量装置,其中,包括超声波信号发生器、前置信号处理电路、超声波发射聚焦换能器、超声波接收换能器、信号调理电路、数据采集器和中心计算机,所述超声波发射聚焦换能器和所述超声波接收换能器之间为测量区;所述超声波信号发生器、所述前置信号处理电路和所述超声波发射聚焦换能器依次连接,所述超声波接收换能器、所述信号调理电路、所述数据采集器和所述中心计算机依次连接;
所述超声波信号发生器由操作者设定发出的超声波参数,所述前置信号处理电路对发出的信号进行放大稳定,经所述超声波发射聚焦换能器发出发散超声波,发散超声波通过所述测量区内具有的颗粒系后由所述超声波接收换能器接收,并经所述信号调理电路对信号进行放大和基本滤波,再由所述数据采集器传输给所述中心计算机处理。
上述测量装置,其中,所述前置信号处理电路设有前置功率放大电路和直流阻断器。
上述测量装置,其中,所述信号调理电路包括非反向放大器、对数放大器、低通滤波器、高通滤波器和积分电路。
另一方面,本发明还提供一种基于上述的测量装置的方法,其中,通过Monte-Carlo方法建模,描述和分析超声的波动,建立颗粒粒度和超声波衰减的关系,求得颗粒的粒度,具体建模步骤如下:
步骤1:采用Monte-carlo方法描述声波动,以声子的形式将声能量离散,通过每个声子的行为描述声波的前进、散射、吸收,当一个声子与颗粒发生散射时,通过随机数判断它是否被吸收或者是散射,通过归一化消声系数Kext和[0,1]区间服从均匀分布随机数ξ的大小比较进行判断,
消声系数Kext的计算通过下述计算得到:
其中,k为声波波数,σ为颗粒投影面积,Re()为取实部运算,An称为散射系数,通过Alex E.Hay和Douglas G.Mercer方法计算;
步骤2:声子遇到颗粒后的随机弹性散射方向散射声压分布的概率模型计算:
式中:θ是散射角;p(θ)是颗粒表面散射声压分布函数,可以由Faran理论计算:
式中:jn和nn分别是第一类球Bessel函数和第二类球Bessel函数,k为声波波数,r为接收点距离,取颗粒半径的1000倍;Pn(cosθ)是勒让德多项式,散射系数Bn由Faran理论公式计算;为确定散射方向,通过另一个[0,1]区间随机数ξ1与声压分布函数比较;
步骤3:声衰减计算,衰减系数A:
式中:Ndet是探测器接收的声子数目;Ntotal是声子样本容量;d是超声发射器与探测接收器之间的距离。
本发明采用上述技术方案,具有以下有益效果:测量***结构简单、廉价,可实现在线测量,可用于实验室科学研究,特别适用于工业现场的应用,测量精度较高。
附图说明
图1示出了本发明的声子传播示意图;
图2示出了本发明声衰减和颗粒粒度关系图;
图3示出了本发明测量装置***结构图。
具体实施方式
以下结合附图说明对本发明的实施例做进一步详细说明,但本实施例并不用于限制本发明,凡是采用本发明的相似方法及其相似变化,均应列入本发明的保护范围。
本发明实施例所提供的一种基于发散超声波方法测量颗粒粒度的装置,其特点是利用发散超声波声变化来对颗粒的粒度进行准确的测量。
发散超声波测量的基本原理是其利用入射发散超声波由于受到颗粒介质的散射和吸收,其透射声强会衰减,通过基于Monte-Carlo方法的建模,可以从理论上预测衰减值,并获得该衰减和颗粒粒度大小的关系,如此一来,通过这一关系,在浓度条件已知的情况下,即可根据超声波衰减求得颗粒度大小。
本发明方法具体实施如下:
如图1所示,方框内部表示介质,小圆圈代表颗粒物,左右小方框分别代表发射器和接收器,本发明由发散超声波探头激发出超声波,按照Monte-Carlo方法,将超声波将设为离散的声子,在传播过程中,每一个声子都会有明确的去向,声子的吸收、散射、前进均可以通过随机数方式来加以判断,在具体的实施中,为了达到统计效果,声子的数目应该不小于1百万个。
接收超声换能器捕获通过测量区的声子信号,并进行统计和计算可以获得声衰减信号,这样即可建立发散声的衰减模型,并且获得衰减与颗粒粒度之间的关系式。
对应的实验中,同样可以先后测出超声发射能量和超声波通过某一颗粒度大小未知的颗粒系能量来计算实验衰减值Amea。在浓度已知情况,根据步骤2得出的数值计算结果作衰减A与颗粒直径D或半径R的函数曲线,由经步骤3得到的实验值Amea在曲线上查到对应的D值。
本发明用于实现基于发散声方法测量颗粒粒度的装置,如图3所示,包括激励产生声波的超声波信号发生器1、对发生信号进行预处理的前置信号处理电路2(超声波发生电路)、把设定的电信号转换为声信号的超声波发射聚焦换能器3、把透过的声信号转换为电信号的超声波接收换能器4和进行信号调理的信号调理电路5,数据采集器6和进行结果处理和分析的中心计算机7。
超声波信号发生器1发出超声波经前置信号处理电路2对原始信号进行放大和调理,激发超声波发射聚焦换能器3组件内部的压电晶片产生指定频率的发散超声波,其波形特征取决于晶片几何形状,原始声波经测量区8与颗粒作用,由超声吸收、散射、前进等构成,超声波接收换能器4收到的声信号经信号调理电路5处理被数据采集器6采集后在中心计算机7上显示保存结果,并进行结果的处理。
超声波信号发生器1可设定产生的超声波参数,包括超声波的频率、超声波的幅值等,前置信号处理电路2是由前置功率放大电路和直流阻断器组成,超声波发射聚焦换能器3的作用是产生发散超声波,超声波接收换能器4是与超声波发射聚焦换能器3有同样频率的换能器,信号调理电路5包括非反向放大器、对数放大器、低通滤波器、高通滤波器、积分电路,非反向放大器用于从超声波信号中去除直流成分,并放大超声波信号,对数放大器用于将放大的声信号振幅保持稳定,积分电路用于从声信号中去除不必要的短信号,中心计算机7用于处理和显示测量结果,其安装有后续进行信号处理的软件。
从上述实施例可以看出,本发明的优势在于:
测量***结构简单、廉价,可实现在线测量,可用于实验室科学研究,特别适用于工业现场的应用,并且测量精度较高。
以上对本发明的具体实施例进行了详细描述,但本发明并不限制于以上描述的具体实施例,其只是作为范例。对于本领域技术人员而言,任何等同修改和替代也都在本发明的范畴之中。因此,在不脱离本发明的精神和范围下所作出的均等变换和修改,都应涵盖在本发明的范围内。
Claims (3)
1.一种基于发散声衰减测量两相体系总颗粒粒度的测量装置,其特征在于,包括超声波信号发生器、前置信号处理电路、超声波发射聚焦换能器、超声波接收换能器、信号调理电路、数据采集器和中心计算机,所述超声波发射聚焦换能器和所述超声波接收换能器之间为测量区;所述超声波信号发生器、所述前置信号处理电路和所述超声波发射聚焦换能器依次连接,所述超声波接收换能器、所述信号调理电路、所述数据采集器和所述中心计算机依次连接;
所述超声波信号发生器由操作者设定发出的超声波参数,所述前置信号处理电路对发出的信号进行放大稳定,经所述超声波发射聚焦换能器发出发散超声波,发散超声波通过所述测量区内具有的颗粒系后由所述超声波接收换能器接收,并经所述信号调理电路对信号进行放大和基本滤波,再由所述数据采集器传输给所述中心计算机处理。
2.根据权利要求1所述的测量装置,其特征在于,所述前置信号处理电路设有前置功率放大电路和直流阻断器。
3.根据权利要求1所述的测量装置,其特征在于,所述信号调理电路包括非反向放大器、对数放大器、低通滤波器、高通滤波器和积分电路。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410215801.XA CN105092430B (zh) | 2014-05-22 | 2014-05-22 | 一种基于发散超声波衰减的颗粒粒度测量装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410215801.XA CN105092430B (zh) | 2014-05-22 | 2014-05-22 | 一种基于发散超声波衰减的颗粒粒度测量装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105092430A CN105092430A (zh) | 2015-11-25 |
CN105092430B true CN105092430B (zh) | 2018-01-19 |
Family
ID=54573352
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410215801.XA Active CN105092430B (zh) | 2014-05-22 | 2014-05-22 | 一种基于发散超声波衰减的颗粒粒度测量装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105092430B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE202017102417U1 (de) * | 2016-04-25 | 2017-07-27 | Gwf Messsysteme Ag | Kompakter Weitwinkel-Schallwandler |
CN106769711A (zh) * | 2017-01-10 | 2017-05-31 | 中国计量大学 | 扫频式超声波聚焦式泥沙粒径分布在线测量仪 |
CN106769733B (zh) * | 2017-01-10 | 2020-09-18 | 中国计量大学 | 超声波聚焦式河流泥沙浓度在线测量仪 |
CN106644859B (zh) * | 2017-01-10 | 2019-05-31 | 中国计量大学 | 一种扫频式超声波聚焦式泥沙粒径分布在线测量方法 |
CN106769734B (zh) * | 2017-01-10 | 2019-03-29 | 中国计量大学 | 一种超声波聚焦式河流泥沙浓度在线测量方法 |
CN108535158A (zh) * | 2018-03-27 | 2018-09-14 | 北京理工大学 | 一种声-电复合的粉尘浓度检测***及检测方法 |
CN108663296B (zh) * | 2018-03-28 | 2020-06-12 | 北京理工大学 | 一种基于双频超声的粉尘浓度检测***及检测方法 |
EP3626174A1 (en) * | 2018-09-20 | 2020-03-25 | Koninklijke Philips N.V. | Switchable grating |
CN110426333B (zh) * | 2019-08-30 | 2021-09-17 | 河海大学常州校区 | 一种利用圆柱体散射声压检测悬浮液颗粒含量的方法 |
CN112525785B (zh) * | 2020-11-25 | 2023-03-24 | 长江水利委员会长江科学院 | 一种基于新型多频超声探头的悬移质含沙量测量方法 |
CN113504159B (zh) * | 2021-06-16 | 2023-01-03 | 中国农业大学 | 不透明散粒体检测分析的方法、装置及电子设备 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6082180A (en) * | 1998-10-21 | 2000-07-04 | Battelle Memorial Institute | Ultrasonic fluid densitometer for process control |
CN201096702Y (zh) * | 2007-09-27 | 2008-08-06 | 上海理工大学 | 一种颗粒粒度和浓度的测量装置 |
CA2761431A1 (en) * | 2009-05-11 | 2010-11-18 | The University Of Western Ontario | Ultrasonic method of monitoring particle size distribution of a medium |
CN102353621B (zh) * | 2011-06-28 | 2013-04-03 | 上海理工大学 | 一种光散射颗粒测量装置及方法 |
CN202661357U (zh) * | 2012-06-27 | 2013-01-09 | 上海理工大学 | 一种测量处于离散状态颗粒粒度分布的装置 |
CN102735595B (zh) * | 2012-07-31 | 2014-07-23 | 上海理工大学 | 基于连续波和猝发波测量离散状态颗粒粒度分布的方法 |
CN102788738A (zh) * | 2012-09-03 | 2012-11-21 | 北京理工大学 | 多相流体密度和浓度超声阵列检测方法 |
CN102830049A (zh) * | 2012-09-03 | 2012-12-19 | 北京理工大学 | 多相流体颗粒粒度超声阵列微纳检测方法 |
-
2014
- 2014-05-22 CN CN201410215801.XA patent/CN105092430B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105092430A (zh) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105092430B (zh) | 一种基于发散超声波衰减的颗粒粒度测量装置及方法 | |
CN101169363B (zh) | 颗粒粒度、浓度和密度测量方法及其装置 | |
Zhang et al. | Longitudinal wave scattering from rough crack-like defects | |
CN102735595B (zh) | 基于连续波和猝发波测量离散状态颗粒粒度分布的方法 | |
CN206192829U (zh) | 一种超低排放烟尘监测装置 | |
CN103983549A (zh) | 一种基于超声脉动原理测量颗粒粒径和浓度的方法 | |
CN108663296B (zh) | 一种基于双频超声的粉尘浓度检测***及检测方法 | |
CN103185554B (zh) | 一种板材厚度检测的装置 | |
CN105300856A (zh) | 基于超声阻抗谱对颗粒浓度和尺寸的测量方法 | |
CN104075968A (zh) | 一种超声粒径分析仪 | |
CN105891416A (zh) | 超声辅助流体中污染物高感度检测***及其工作方法 | |
Divós et al. | Resolution of stress wave based acoustic tomography | |
CN104569483A (zh) | 同时测量泥沙浓度和三维瞬时流速的超声传感器 | |
CN105651662A (zh) | 气溶胶质量浓度光学检测装置及其检测方法 | |
CN111351540A (zh) | 一种气力输送过程中颗粒质量流率的检测方法和*** | |
CN201096702Y (zh) | 一种颗粒粒度和浓度的测量装置 | |
CN110296913B (zh) | 一种可燃粉尘扩散动态浓度的检测***及其检测方法 | |
CN202661357U (zh) | 一种测量处于离散状态颗粒粒度分布的装置 | |
CN106769711A (zh) | 扫频式超声波聚焦式泥沙粒径分布在线测量仪 | |
CN112525785B (zh) | 一种基于新型多频超声探头的悬移质含沙量测量方法 | |
CN103411858B (zh) | 基于阵列式声传感器的颗粒粒度分布在线测量装置及方法 | |
CN209559714U (zh) | 一种超低排放烟尘浓度监测装置 | |
Ibrahim et al. | A review on ultrasonic process tomography system | |
CN105424557A (zh) | 锅炉煤粉颗粒参数的在线测量装置及方法 | |
CN112525778B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20190626 Address after: Room 213-02, Block A, Building 11, 128 Xiangyin Road, Yangpu District, Shanghai Patentee after: Shanghai surplus Precision Instrument Co., Ltd. Address before: 315336 Room 101, Unit 2, Building 15, Tongji Industrial Park, 77 Binhai Road, Ningbo City, Zhejiang Province Patentee before: NINGBO YINGNUO INSTRUMENT MANUFACTURING CO., LTD. |