发明名称 颗粒粒度、浓度和密度测量方法及其装置
摘要 一种颗粒粒度、浓度和密度测量方法及其装置,涉及超声测量技术领域;所要解决的是提高超声测量通用性和准确性的技术问题;该测量装置包括依次连接的计算机、信号处理电路、脉冲波发射/接收电路、宽带换能器,宽带换能器设于缓冲板外侧;所述计算机处理程序:将信号作快速傅立叶变换获得超声幅值相位谱,换算得反射系数,声特性阻抗,声衰减系数和声速;由测量值能计算得混合物密度,由已知的颗粒物和连续介质密度计算得浓度值;再由直接反射波和透射回波作频谱分析得到的声衰减谱和理论声衰减谱的误差作为目标函数,按照优化方法进行优化,计算得颗粒粒度分布;本发明具有基于反射式超声信号测量和分析的,通用性强,测量结果准确的特点。
申请公布号 CN101169363B 申请公布日期 2010.12.22
申请号 CN200710046525.9 申请日期 2007.09.27
申请人 上海理工大学 发明人 苏明旭;蔡小舒;薛明华
分类号 G01N15/00(2006.01)I;G01N15/02(2006.01)I;G01N15/06(2006.01)I;G01N29/02(2006.01)I;G06F19/00(2006.01)I 主分类号 G01N15/00(2006.01)I
代理机构 上海申汇专利代理有限公司 31001 代理人 吴宝根
主权项 1.一种颗粒粒度、浓度和密度测量方法,其特征是,包括以下步骤:1)采集直接反射波和透射回波:由收/发换能器发射的单个窄脉冲超声波反射后所述收/发换能器获得直接反射波和透射回波信号;2)获得反射系数和声学特征阻抗:在数据处理的计算机中对时域信号作快速傅立叶变换得超声幅值和相位谱,由公式:<img file="FSB00000229414600011.GIF" wi="182" he="123" />和<img file="FSB00000229414600012.GIF" wi="323" he="143" />计算获得反射系数和声学特征阻抗;其中,R<sub>s</sub>和R<sub>c</sub>分别为对样品和标定物质进行测试时的反射系数,M<sub>s</sub>和M<sub>c</sub>为对应超声信号的幅度;Z<sub>s</sub>表示颗粒和连续介质的两相体混合物的声学特征阻抗,Z<sub>b</sub>表示壁面的声学特征阻抗;3)获得衰减系数和声速:在数据处理的计算机中由快速傅立叶变换得直接反射波和透射回波的幅值谱并修正了声反射,透射和扩散损失后的声衰减谱;由直接反射波和透射回波的时差可测声速<maths num="0001"><![CDATA[<math><mrow><msub><mi>&alpha;</mi><mi>meas</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>L</mi></mfrac><mi>ln</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>M</mi><mi>A</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>M</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><msub><mi>&alpha;</mi><mi>c</mi></msub><mo>;</mo></mrow></math>]]></maths>其中M<sub>A</sub>(f)为直接反射波幅值,M<sub>B</sub>(f)为透射回波幅值,L为声程,而α<sub>c</sub>为考虑声反射,扩散效果的超声衰减修正系数;4)计算颗粒和连续介质的两相混合物等效密度,在数据处理的计算机中由公式:ρ<sub>s</sub>=Z<sub>s</sub>/c能直接计算混合物等效密度;5)计算颗粒浓度,在数据处理的计算机中由公式:φ=[ρ<sub>p</sub>(ρ<sub>s</sub>-ρ<sub>l</sub>)]/[ρ<sub>s</sub>(ρ<sub>p</sub>-ρ<sub>l</sub>)]能直接计算颗粒浓度;其中ρ<sub>s</sub>为混合物等效密度,ρ<sub>p</sub>为颗粒相密度,ρ<sub>l</sub>为连续相密度;6)计算颗粒粒度:在数据处理的计算机中按照高浓度颗粒两相体系中复波数公式:<maths num="0002"><![CDATA[<math><mrow><msup><mi>&kappa;</mi><mn>2</mn></msup><mo>=</mo><msup><mi>&omega;</mi><mn>2</mn></msup><msubsup><mi>k</mi><mi>a</mi><mo>*</mo></msubsup><mo>&times;</mo><mfrac><mrow><mi>&rho;</mi><mo>[</mo><msup><mi>&rho;</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&phi;</mi><mo>+</mo><mi>&phi;S</mi><mo>)</mo></mrow><mo>+</mo><mi>&rho;S</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>]</mo></mrow><mrow><msup><mi>&rho;</mi><mo>&prime;</mo></msup><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&phi;</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>&rho;</mi><mo>[</mo><mi>S</mi><mo>+</mo><mi>&phi;</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>]</mo></mrow></mfrac></mrow></math>]]></maths>式中:φ为颗粒浓度,R为颗粒半径,<img file="FSB00000229414600021.GIF" wi="41" he="55" />为绝热压缩系数,S是浓度,粒度和物性的函数,ω为角频率;声衰减系数和声速可按下式计算,衰减系数:α=-Im(κ)声速:c=ω/Re(κ)为求解颗粒粒度分布,按照上面计算的理论声衰减系数(α<sub>theory</sub>(f))与实验测量信号换算得的声衰减系数α<sub>meas</sub>(f)=ln[M<sub>A</sub>(f)/M<sub>B</sub>(f)]/L-α<sub>c</sub>构造误差函数,<maths num="0003"><![CDATA[<math><mrow><mi>E</mi><mo>=</mo><msqrt><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>f</mi></msub></munderover><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>&alpha;</mi><mi>meas</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>&alpha;</mi><mi>theory</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow><msub><mi>&alpha;</mi><mrow><mi>theory</mi><mo>,</mo><mi>i</mi></mrow></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow><msub><mi>N</mi><mi>f</mi></msub></mfrac></msqrt><mo>;</mo></mrow></math>]]></maths>然后按照优化方法进行优化求解得颗粒粒度分布。
地址 200093 上海市杨浦区军工路516号