CN201096702Y - 一种颗粒粒度和浓度的测量装置 - Google Patents
一种颗粒粒度和浓度的测量装置 Download PDFInfo
- Publication number
- CN201096702Y CN201096702Y CNU2007200750624U CN200720075062U CN201096702Y CN 201096702 Y CN201096702 Y CN 201096702Y CN U2007200750624 U CNU2007200750624 U CN U2007200750624U CN 200720075062 U CN200720075062 U CN 200720075062U CN 201096702 Y CN201096702 Y CN 201096702Y
- Authority
- CN
- China
- Prior art keywords
- concentration
- measurement
- receiving transducer
- ultrasonic
- grain graininess
- 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
Abstract
一种颗粒粒度和浓度测量装置,涉及超声测量技术领域;所要解决的是颗粒粒度和浓度测量技术问题;该测量装置包括进行数据处理的计算机,连接计算机的信号处理电路,连接信号处理电路的脉冲波发射/接收电路,连接脉冲波发射和接收电路的发射/接收换能器组件;数据处理程序将测量信号与所设标准物质进行对比;实测的声衰减谱与按照声散射公式计算理论声衰减谱一起,构建矩阵和线性方程组并求解;用双频率衰减比值法实现对平均粒度的计算,以优化方法计算粒度分布,并根据颗粒粒度反算出颗粒浓度。本实用新型具有能实现完全的非接触测量的,基于超声散射衰减谱分析的,能降低测量误差的,能在线测量高浓度和高衰减的特点。
Description
技术领域
本实用新型涉及超声测量技术,特别是涉及一种用于两相悬浊液或乳浊液中离散状态颗粒的粒度分布和浓度的测量技术。
背景技术
两相流颗粒测量技术是指对其中分散状颗粒的粒度分布和浓度进行测量,在许多涉及到两相流的动力、化工、医药、环保、材料、水利等领域中具有广泛应用背景。该应用的共同特点是颗粒物质需要在一种在线,非接触,无需稀释(较高浓度)的状态下进行测量。现有的颗粒测量方法例如筛分法、显微镜法、全息照相法、电感应法、沉降法和光散射法等,都很难实现上述条件的测量。
由于超声波具有很宽的频带范围,可以穿透几乎任何物质,可在有色甚至不透明的物质中传播并具有测量速度快,容易实现测量和数据处理自动化等优点,并且超声在颗粒物质中的传播机理与颗粒粒度和浓度密切相关,因此可被用作颗粒物信息载体实现颗粒物质的粒度和浓度测量。
现有的超声颗粒浓度测量方法中,采用简单经验公式或事先用样品标定,没有考虑到处于多分散分布状态的颗粒粒度不同对超声衰减量(或声速)测量带来的影响,事实上颗粒粒度及分布的不同对于超声衰减(或声速)影响非常大,忽略该影响导致所测浓度误差大,不能得到准确的测量结果,更不可能获得颗粒的粒度。
实用新型内容
针对上述现有技术中存在的缺陷,所要解决的技术问题是提供一种能实现完全的非接触测量的,基于超声散射衰减谱分析的,能降低测量误差的,能在线测量高浓度和高衰减的颗粒粒度和浓度测量装置。
为了解决上述技术问题,本实用新型所提供的一种颗粒粒度浓度测量装置,其特征在于,包括进行数据处理的计算机,连接计算机的信号处理电路,连接信号处理电路的脉冲波发射/接收电路,连接脉冲波发射和接收电路的至少一个发射/接收换能器组件;所述计算机为能进行将测量信号根据设定公式计算得平均粒度和根据颗粒粒度反算出颗粒浓度的计算机。
进一步的,所述信号处理电路设有高速A/D转换单元和信号放大单元。
进一步的,所述一对以上的超声发射/接收换能器在样品池管道壁二侧平行排列。
进一步的,所述发射/接收换能器组件为一发一收结构,发射换能器设于样品池管道壁的一外侧,接收换能器设于样品池管道壁的对面的另一外侧。
进一步的,所述发射/接收换能器组件为自发自收结构,发射/接收换能器设于样品池管道壁的一外侧,反射板设于样品池管道壁的对面的另一外侧。
利用本实用新型提供的颗粒粒度浓度测量装置,由于采用基于超声散射衰减谱分析的超声换能器的非接触测量,降低对测量对象的干扰,能降低测量误差,能实现在线测量高浓度和高衰减的颗粒粒度和浓度。由于采用用双频率衰减比值方法实现对平均粒度的计算,以优化方法计算粒度分布,并根据粒度和分布反算出颗粒的浓度,进而降低测量误差,提高测量精度。
附图说明
图1是本实用新型实施例的颗粒粒度浓度测量装置的一发一收装置的结构示意图;
图2是本实用新型实施例的颗粒粒度浓度测量装置的自发自收装置的结构示意图;
图3是本实用新型实施例的颗粒粒度浓度测量装置的在线测量探针的结构示意图。
具体实施方式
以下结合附图说明对本实用新型的实施例作进一步详细描述,但本实施例并不用于限制本实用新型,凡是采用本实用新型的相似结构及其相似变化,均应列入本实用新型的保护范围。
本实用新型考虑到声衰减或速度的频谱(指其随声频率的频域曲线)与颗粒粒径密切相关,这一关系完全可以被应用到离散状态的颗粒粒度测量中;在本实用新型实施例中,考虑超声在颗粒介质中传播时的散射效应,计算超声波通过颗粒两相体系后的理论预测衰减谱;
式中:α为与超声频率f有关的超声衰减系数,φ为颗粒浓度,R为颗粒半径,k为连续介质波数,An为n阶声散射系数,该级数的最大阶数取值视颗粒,超声频率即精度要求而定。
将测量数据经快速傅立叶处理后计算得超声信号幅度谱和相位谱,与标准物质对比即可同时由测量超声信号计算出超声在颗粒两相体系中的衰减谱:
Mc和Ms分别为标定物质和样品信号作傅立叶变换后幅值,L为声程,而αc为考虑声反射,扩散效果的超声衰减修正系数。
计算颗粒粒度分布:由公式:
计算出不同颗粒粒度(R),不同频率(f)下超声衰减系数,并且据此构造矩阵T,
利用实验手段获得的数个频率f1至fm超声衰减系数构造列向量,
gmeas=[ameas(f1)…ameas(fm)]T,
可以写成线性方程组
Tx=gmeas;
其中,待求解向量x表征了从R1至Rn的不同颗粒粒度的离散相对频度分布。由于此线性方程组求解涉及病态问题且解向量必须为非负,通常不采用高斯消元法等求解。最优正则化方法具体为:重新构造一个线性方程组
(TTT+γH)x=TTgmeas
矩阵H具有如下形式
而光顺因子γ则采用对下式作最优化获得(即寻找使得该式取值最小的γmin)
式中,符号tr表示求迹运算,上标T为矩阵转置运算,上标-1指矩阵求逆。在H和γ获得后,可以用非负最小二乘法求解线性方程组,得到颗粒粒度分布;本实用新型除了采用最优正则化方法之外,还可用DFP(Davidon-Fletcher-Powell)方法、最速下降法(Steepest Descent method)和模拟退火法(Simulated Annealingmethod)等优化方法;
根据粒度分布反算出颗粒的浓度:按公式:
求得平均粒度(N为对应粒度R的颗粒粒度分布),将R32回代至公式
同时将测得超声衰减系数带入公式等号左侧,即可求出唯一的未知变量φ(即颗粒体积浓度)。
综上所述,本实用新型实施例的颗粒粒度和浓度测量原理,包括以下步骤:
1)反射法(或透射)测量样品和标定物质:(超)高频窄脉冲电信号触发收/发换能器发射的单个窄脉冲信号多次反射后所述收/发换能器采集多次超声反射信号,经信号处理电路对采集的多次超声反射信号进行放大、模数转换获得测量数据;
2)计算理论和实验超声衰减谱:在数据处理的计算机中利用公式计算超声波通过颗粒两相体系后的理论预测衰减谱:将测量数据经快速傅立叶处理后计算得超声信号幅度谱和相位谱,与标准物质对比即可同时由测量超声信号计算出超声在颗粒两相体系中的衰减谱:
φ为颗粒浓度,R为颗粒半径,k为连续介质波数,An为n阶声散射系数。该式子建立了超声衰减谱和颗粒粒度,浓度,以及超声频率之间联系。
将测量数据经快速傅立叶变换处理后计算得超声信号幅度谱和相位谱,与标准物质对比即可(由测得超声信号)换算出超声在颗粒两相体系中的衰减谱:
其中Mc和Ms分别为标定物质和样品信号作傅立叶变换后幅值,L为声程,而αc为考虑声反射,扩散效果的超声衰减修正系数;
3)采用双频(二个频率)超声衰减比值法在数据处理的计算机中计算平均粒度;即由步骤2)中理论公式计算多个不同颗粒半径R对应于选定的二个不同频率下超声衰减系数比值αtheory(f1,R)/αtheory(f2,R)绘制成一条曲线(横坐标为R,纵坐标为超声衰减系数比值αtheory(f1,R)/αtheory(f2,R)),同时将测得频率f1和f2下的超声衰减系数,换算成比值αmeas(f1)/αmeas(f2),过该比值作一与横坐标之平行线后,其与前面绘制曲线交点处横坐标值即对应了平均粒度(以半径表示)的求解结果。
4)计算颗粒粒度分布:在数据处理的计算机中由公式:
计算出不同颗粒粒度(R),不同频率(f)下超声衰减系数,并且据此构造矩阵T,
利用实验手段获得的数个频率f1至fm超声衰减系数构造列向量,
gmeas=[ameas(f1)…ameas(fm)]T,
可以写成线性方程组
Tx=gmeas;
其中,待求解向量x表征了从R1至Rn的不同颗粒粒度的离散相对频度分布。由于此线性方程组求解涉及病态问题且解向量必须为非负,通常不采用高斯消元法等求解。我们的最优正则化方法具体为:重新构造一个线性方程组
(TTT+γH)x=TTgmeas
矩阵H具有如下形式
而光顺因子γ则采用对下式作最优化获得(即寻找使得该式取值最小的γmin)
式中,符号tr表示求迹运算,上标T为矩阵转置运算,上标-1指矩阵求逆。在H和γ获得后,之后确定光顺参数,通过非负最小二乘方法计算得出颗粒粒度分布;本实用新型除了采用最优正则化方法之外,还可用DFP(Davidon-Fletcher-Powell)方法、最速下降法(Steepest Descent method)和模拟退火法(Simulated Annealingmethod)等优化方法;
5)根据粒度分布反算出颗粒的浓度:在数据处理的计算机中,按公式:
求得平均粒度,将其回代至公式
即可求出唯一的未知变量φ(即颗粒体积浓度);
本实用新型实施例所提供的一种颗粒粒度浓度测量装置,包括进行数据处理的计算机,连接计算机的设有高速A/D转换单元和信号放大单元的信号处理电路,连接信号处理电路的脉冲波发射/接收电路,分别连接脉冲波发射和接收电路的发射换能器和接收换能器;如图1所示,发射换能器1设于样品池2管道壁的一外侧,接收换能器3设于样品池2管道壁的对面的另一外侧;也可采用多对超声发射/接收换能器在样品池管道壁二侧平行排列;本测量装置采用一个换能器发射宽带超声波,由脉冲信号触发的超声波通过超声换能器间的测量区,并与颗粒两相介质作用;接收换能器的电信号经放大、数模转换传输入进行数据处理的计算机,经快速傅立叶处理后计算得超声信号幅度谱和相位谱,与标准物质对比计算出超声衰减谱和速度谱。部分声能由于散射作用不能到达对面的接收换能器,为求解超声波的衰减程度,将测量声信号与某一标准物质(如纯水)进行对比。实测的超声衰减,速度谱与按照声散射公式计算理论声衰减和速度谱对比,构建矩阵和线性方程组,并求解结果。
如图2所示,本实用新型实施例颗粒粒度浓度测量装置的超声发射/接收换能器组件也可采用自发自收的结构,即发射/接收换能器4设于样品池2’管道壁的一外侧,反射板5设于样品池2’管道壁的对面的另一外侧;采用一个发射/接收换能器4发射宽带超声波,通过待测颗粒两相介质,在对面反射板5反射后又由同样的发射/接收换能器4接收,接收换能器的电信号经放大、数模转换传输入进行数据处理的计算机,经快速傅立叶处理后计算得超声信号幅度谱和相位谱,与标准物质对比计算出超声衰减谱和速度谱。为求解超声波的衰减程度,将测量超声信号与某一标准物质(如纯水)进行对比。实测的声衰减,速度谱与按照超声散射公式计算理论超声衰减和速度谱对比,构建矩阵和线性方程组,并求解结果。
如图3所示,本实用新型另一实施例采用宽带探头,将测量部件做成探头形式,在探头凹槽的一侧设保护缓冲板8,保护缓冲板8外侧设宽带换能器7,并由电缆线6连接脉冲波发射/接收电路,在探头凹槽对面的另一侧设反射板9;宽带探头以法兰形式和管道连接或插入管道中作在线检测,对时域信号作快速傅立叶变换得超声幅值相位谱,为反演过程提供了丰富的频谱信息。
本实用新型实施例用本装置能对样品池或在线管道内颗粒两相介质的粒度和分布、浓度进行测量。在线管道测量时,颗粒粒度浓度测量装置的超声发射/接收换能器组件通过法兰连接,亦可直接附着在管道外壁上。测量结果以曲线和数据的形式直接显示在计算机屏幕上,同时保存在硬盘上供以后分析使用。
本实用新型所述计算机包括单板机、单片机、可编程微处理器和DSP(数字处理芯片等。
Claims (5)
1、一种颗粒粒度和浓度的测量装置,其特征在于,包括进行数据处理的计算机,连接计算机的信号处理电路,连接信号处理电路的脉冲波发射/接收电路,连接脉冲波发射和接收电路的至少一个发射/接收换能器组件;所述计算机为能进行将测量信号根据设定公式计算得平均粒度和根据颗粒粒度反算出颗粒浓度的计算机。
2、根据权利要求1所述的颗粒粒度和浓度的测量装置,其特征在于,所述信号处理电路设有高速A/D转换单元和信号放大单元。
3、根据权利要求1所述的颗粒粒度和浓度的测量装置,其特征在于,所述一对以上的超声发射/接收换能器在样品池管道壁二侧平行排列。
4、根据权利要求1所述的颗粒粒度和浓度的测量装置,其特征在于,所述发射/接收换能器组件为一发一收结构,发射换能器设于样品池管道壁的一外侧,接收换能器设于样品池管道壁的对面的另一外侧。
5、根据权利要求1所述的颗粒粒度和浓度的测量装置,其特征在于,所述发射/接收换能器组件为自发自收结构,发射/接收换能器设于样品池管道壁的一外侧,反射板设于样品池管道壁的对面的另一外侧。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNU2007200750624U CN201096702Y (zh) | 2007-09-27 | 2007-09-27 | 一种颗粒粒度和浓度的测量装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNU2007200750624U CN201096702Y (zh) | 2007-09-27 | 2007-09-27 | 一种颗粒粒度和浓度的测量装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN201096702Y true CN201096702Y (zh) | 2008-08-06 |
Family
ID=39923755
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNU2007200750624U Expired - Fee Related CN201096702Y (zh) | 2007-09-27 | 2007-09-27 | 一种颗粒粒度和浓度的测量装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN201096702Y (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102985816A (zh) * | 2010-07-08 | 2013-03-20 | 恩德斯+豪斯流量技术股份有限公司 | 超声颗粒测量系统 |
CN103969327A (zh) * | 2014-06-04 | 2014-08-06 | 江南大学 | 一种超声波油烟浓度传感器 |
CN105092430A (zh) * | 2014-05-22 | 2015-11-25 | 宁波盈诺仪器制造有限公司 | 一种基于发散超声波衰减的颗粒粒度测量装置及方法 |
CN105842130B (zh) * | 2016-05-18 | 2019-04-05 | 陕西煤业化工技术研究院有限责任公司 | 一种非牛顿基液流体纳米颗粒团聚超声监测装置及方法 |
CN111735741A (zh) * | 2020-07-03 | 2020-10-02 | 中国人民解放军国防科技大学 | 一种不规则生物颗粒毫米波透过率的测试装置 |
CN114252377A (zh) * | 2021-12-28 | 2022-03-29 | 聪明猪检测技术(成都)有限公司 | 基于多传感器管道介质中颗粒含量检测方法及系统 |
CN114252377B (zh) * | 2021-12-28 | 2024-05-03 | 聪明猪检测技术(成都)有限公司 | 基于多传感器管道介质中颗粒含量检测方法及系统 |
-
2007
- 2007-09-27 CN CNU2007200750624U patent/CN201096702Y/zh not_active Expired - Fee Related
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102985816A (zh) * | 2010-07-08 | 2013-03-20 | 恩德斯+豪斯流量技术股份有限公司 | 超声颗粒测量系统 |
CN105092430A (zh) * | 2014-05-22 | 2015-11-25 | 宁波盈诺仪器制造有限公司 | 一种基于发散超声波衰减的颗粒粒度测量装置及方法 |
CN103969327A (zh) * | 2014-06-04 | 2014-08-06 | 江南大学 | 一种超声波油烟浓度传感器 |
CN103969327B (zh) * | 2014-06-04 | 2016-07-06 | 江南大学 | 一种超声波油烟浓度传感器 |
CN105842130B (zh) * | 2016-05-18 | 2019-04-05 | 陕西煤业化工技术研究院有限责任公司 | 一种非牛顿基液流体纳米颗粒团聚超声监测装置及方法 |
CN111735741A (zh) * | 2020-07-03 | 2020-10-02 | 中国人民解放军国防科技大学 | 一种不规则生物颗粒毫米波透过率的测试装置 |
CN114252377A (zh) * | 2021-12-28 | 2022-03-29 | 聪明猪检测技术(成都)有限公司 | 基于多传感器管道介质中颗粒含量检测方法及系统 |
CN114252377B (zh) * | 2021-12-28 | 2024-05-03 | 聪明猪检测技术(成都)有限公司 | 基于多传感器管道介质中颗粒含量检测方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101135626A (zh) | 一种颗粒粒度和浓度测量方法及其装置 | |
CN101169363B (zh) | 颗粒粒度、浓度和密度测量方法及其装置 | |
CN101169364B (zh) | 对离散状态颗粒粒度分布测量的方法及其装置 | |
Calmon et al. | CIVA: An expertise platform for simulation and processing NDT data | |
Malyarenko et al. | Ultrasonic Lamb wave diffraction tomography | |
CN102735595B (zh) | 基于连续波和猝发波测量离散状态颗粒粒度分布的方法 | |
Howard et al. | On the probability of detecting wall thinning defects with dispersive circumferential guided waves | |
CN102980945B (zh) | 基于双谱分析的结构微裂纹混频非线性超声检测方法 | |
CN201096702Y (zh) | 一种颗粒粒度和浓度的测量装置 | |
Ghoshal et al. | Wigner distribution of a transducer beam pattern within a multiple scattering formalism for heterogeneous solids | |
Benstock et al. | The influence of surface roughness on ultrasonic thickness measurements | |
CN106226392A (zh) | 基于超声衰减机理模型的油水两相流相含率测量方法 | |
CN102830049A (zh) | 多相流体颗粒粒度超声阵列微纳检测方法 | |
CN105300856B (zh) | 基于超声阻抗谱对颗粒浓度和尺寸的测量方法 | |
Fan et al. | Non-contact ultrasonic gas flow metering using air-coupled leaky Lamb waves | |
CN105092430A (zh) | 一种基于发散超声波衰减的颗粒粒度测量装置及方法 | |
CN201096703Y (zh) | 颗粒粒度、浓度和密度测量装置 | |
US4301684A (en) | Ultrasonic non-destructive evaluation technique for structures of complex geometry | |
Tian et al. | Energy peak fitting of echo based signal processing method for ultrasonic gas flow meter | |
CN102607479B (zh) | 基于声压反射系数功率谱测量超声在薄层介质中往返时间的方法 | |
CN104075968A (zh) | 一种超声粒径分析仪 | |
Bloxham et al. | Combining simulated and experimental data to simulate ultrasonic array data from defects in materials with high structural noise | |
CN114994175A (zh) | 模态分解双谱分析的空耦超声应力检测装置及方法 | |
Kalashnikov et al. | Errors and uncertainties in the measurement of ultrasonic wave attenuation and phase velocity | |
Lohne et al. | Measurements and 3D simulations of ultrasonic directive beam transmission through a water-immersed steel plate |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20080806 Termination date: 20100927 |