基于超聲波阻抗譜的顆粒粒徑表征方法

紀(jì)曉明，蘇明旭，汪雪，蔡小舒(上海理工大學(xué)顆粒與兩相流測(cè)量研究所，上海市動(dòng)力工程多相流動(dòng)與傳熱重點(diǎn)實(shí)驗(yàn)室，上海 200093)

?

基于超聲波阻抗譜的顆粒粒徑表征方法

紀(jì)曉明，蘇明旭，汪雪，蔡小舒
(上海理工大學(xué)顆粒與兩相流測(cè)量研究所，上海市動(dòng)力工程多相流動(dòng)與傳熱重點(diǎn)實(shí)驗(yàn)室，上海 200093)

摘要：研究超聲反射波譜與顆粒粒徑及濃度之間的關(guān)系并發(fā)展一種基于超聲波阻抗譜的顆粒粒徑表征方法。對(duì)超聲波波動(dòng)理論模型適當(dāng)變形，建立超聲阻抗譜與顆粒粒徑及濃度之間的關(guān)系；進(jìn)一步通過(guò)數(shù)值模擬分析超聲阻抗譜對(duì)顆粒濃度及粒徑變化的敏感性；實(shí)驗(yàn)中使用中心頻率10～100 MHz超聲波換能器，利用自發(fā)自收模式對(duì)超聲波在緩沖層與介質(zhì)界面上的反射波信號(hào)進(jìn)行測(cè)量并分析，對(duì)體積中位徑分別為7.69、21.58、66.64 μm的聚苯乙烯懸浮液進(jìn)行實(shí)驗(yàn)，獲得阻抗譜并與數(shù)值模擬結(jié)果對(duì)比，根據(jù)實(shí)驗(yàn)阻抗譜進(jìn)行數(shù)據(jù)反演獲得樣品顆粒粒度分布，并與圖像分析結(jié)果進(jìn)行了對(duì)比，結(jié)果顯示，本方法可有效分辨3種顆粒樣品的粒徑。

關(guān)鍵詞：超聲阻抗譜；聚苯乙烯懸浮液；粒度分布；濃度；測(cè)量；數(shù)值模擬；反演

2015-12-23收到初稿，2016-03-16收到修改稿。

聯(lián)系人：蘇明旭。 第一作者：紀(jì)曉明（1991—），男，碩士研究生。

Received date: 2015-12-23.

Foundation item: supported by the National Natural Science Foundation of China (51176128,5120611).

引　言

顆粒存在于工業(yè)中的各個(gè)領(lǐng)域，其粒徑大小及濃度高低會(huì)影響相關(guān)醫(yī)藥產(chǎn)品質(zhì)量、食品保質(zhì)期、能源利用效率以及設(shè)備使用壽命，所以對(duì)顆粒特性的表征得到了廣泛的關(guān)注和研究。目前已有較多方法得到發(fā)展和應(yīng)用，如動(dòng)態(tài)光散射[1]、多波長(zhǎng)消光法[2]和圖像法[3-6]等，但是上述方法都存在著不足，即對(duì)于高濃度顆粒懸浮液不宜直接測(cè)量而需要預(yù)先稀釋?zhuān)芏嗲闆r下不恰當(dāng)?shù)南♂屵^(guò)程會(huì)額外引入誤差或使得待測(cè)樣品缺乏代表性。由于超聲波穿透性較強(qiáng)，利用超聲法對(duì)濃度相對(duì)較高[7]（體積分?jǐn)?shù)大于20%）的顆粒懸浮液進(jìn)行測(cè)量較為合適。

目前，常見(jiàn)的懸浮液顆粒尺寸和濃度超聲法測(cè)量原理是基于超聲衰減譜和相速度譜與顆粒粒徑及濃度之間的關(guān)系，這一關(guān)系在數(shù)學(xué)上形成超聲顆粒粒徑表征的理論模型[8-9]。實(shí)驗(yàn)中，超聲波通過(guò)測(cè)量區(qū)后由換能器接收，通過(guò)分析接收信號(hào)獲得衰減譜和相速度譜，結(jié)合模型反演出顆粒粒徑及濃度[10-11]。然而，當(dāng)測(cè)量對(duì)象濃度更高或介質(zhì)本身對(duì)超聲有較強(qiáng)衰減（如甘油）時(shí)，可能會(huì)由于超聲信號(hào)衰減太大導(dǎo)致接收換能器無(wú)法獲得信號(hào)，或所獲信號(hào)信噪比太差，從而影響測(cè)量結(jié)果的精度甚至不能測(cè)量。一種處理方法是盡量減小發(fā)射和接收換能器間距（聲程），但這一方面可能使得測(cè)量區(qū)處于近場(chǎng)區(qū)，另一方面在很多現(xiàn)場(chǎng)應(yīng)用中聲程不可調(diào)。

與透射式測(cè)量不同，對(duì)超聲在緩沖層與介質(zhì)界面上的反射波信號(hào)進(jìn)行測(cè)量時(shí)，只要選取合適的緩沖層材料（即恰當(dāng)?shù)淖杩共睿┚涂梢缘玫叫旁氡容^好的信號(hào)，且由于無(wú)須考慮超聲波在介質(zhì)中傳播過(guò)程，所以信號(hào)強(qiáng)度對(duì)其衰減特性也并不敏感。此外，超聲特性阻抗與顆粒介質(zhì)本身特性緊密相關(guān)，原理上利用超聲反射波信號(hào)表征顆粒特性是可行的。目前已有一些研究從實(shí)驗(yàn)角度研究反射信號(hào)與被測(cè)介質(zhì)特性，如含氣食品中氣泡性質(zhì)的表征[12]、蜂蜜的動(dòng)態(tài)剪切流變性質(zhì)[13]、靜態(tài)軸承內(nèi)潤(rùn)滑油的黏度分布[14]，但是對(duì)于顆粒懸浮液中顆粒度進(jìn)行表征還未見(jiàn)報(bào)道。為此，本文研究一種超聲阻抗實(shí)部譜表征顆粒粒徑的方法，并將其與顯微鏡圖像法進(jìn)行對(duì)比。

1　超聲阻抗譜顆粒表征原理

1.1超聲阻抗原理

由超聲顆粒表征中ECAH[8-9]（Epstein-Carhart-Allegra-Hawley）理論模型出發(fā)，建立懸浮液中顆粒粒徑及濃度與超聲阻抗之間的關(guān)系。復(fù)波數(shù)κ表達(dá)式如下

式中，κc=2π/λ為入射壓縮波波數(shù)，φ為介質(zhì)濃度，R為顆粒半徑，n為顆粒粒度分散等級(jí)，An為散射系數(shù)，由壓縮波函數(shù)結(jié)合邊界條件得到。

按照復(fù)波數(shù)定義

可見(jiàn)，復(fù)波數(shù)的實(shí)部與相速度[c(ω)]有關(guān)，而其虛部則等于衰減系數(shù)[α(ω)]，其中ω為超聲波角頻率。超聲復(fù)阻抗（Z）定義由式（3）[15]給出。

式中，ρ為介質(zhì)等效密度，λ為超聲波波長(zhǎng)。聯(lián)立式(1)～式(3)，有

式中，Re為取實(shí)部運(yùn)算。式（4）建立了超聲復(fù)阻抗與顆粒特性（如粒徑、濃度及密度等）之間的關(guān)系。

1.2數(shù)值計(jì)算

利用表1的物性參數(shù)[16]，結(jié)合式（4），進(jìn)行了一系列的數(shù)值計(jì)算，討論了聚苯乙烯顆粒-水懸浮液中顆粒的粒徑及濃度的變化對(duì)超聲阻抗實(shí)部的影響。模擬中聚苯乙烯顆粒粒徑大小為10、30、50、100 μm，聚苯乙烯-水懸浮液體積分?jǐn)?shù)分別為10%、30%和50%。

表1　水、聚苯乙烯物性參數(shù)(20℃)Table 1　Physical parameters of water and polystyrene

圖1為聚苯乙烯懸浮液中聲阻抗實(shí)部在不同粒徑和濃度條件下隨超聲頻率的變化曲線。每條阻抗譜曲線均對(duì)應(yīng)一特征頻率，該特征頻率與顆粒對(duì)超聲波的共振散射有關(guān)[17]。低于此頻率，阻抗穩(wěn)定在某一水平值，直至到達(dá)該頻率時(shí)出現(xiàn)明顯振蕩，信號(hào)迅速降至最小值（谷），又迅速升至最大值（峰），之后振蕩逐次衰減直至重新趨于穩(wěn)定。在圖1(a)中，該特征頻率隨著顆粒粒徑增大而變小，但曲線之間變化趨勢(shì)相似。在圖1(b)中，改變懸浮液顆粒濃度，振蕩幅度(峰值和谷值差)隨濃度增加而遞增，但其特征頻率并無(wú)明顯變化。簡(jiǎn)單地說(shuō)，阻抗譜特征頻率和顆粒粒徑具有直接關(guān)聯(lián)而其振蕩幅度與濃度直接關(guān)聯(lián)。在超聲衰減和速度譜方法中，顆粒粒徑和濃度會(huì)同時(shí)影響衰減或速度譜，這使得粒徑和濃度被嚴(yán)重地耦合在一起，而阻抗譜中特征頻率和幅度的分離使得上述兩個(gè)參數(shù)的檢測(cè)相對(duì)獨(dú)立，對(duì)于降低問(wèn)題復(fù)雜性非常有利。理論上，阻抗譜虛部同樣可以用于顆粒表征，但前期分析發(fā)現(xiàn)其敏感度較實(shí)部低，所以本文只討論利用阻抗譜實(shí)部對(duì)顆粒特性進(jìn)行表征。

圖1　聚苯乙烯-水懸浮液中聲阻抗譜數(shù)值模擬結(jié)果（實(shí)部）Fig.1　Numerical results of real part of acoustic impedance spectra for polystyrene suspensions

2　實(shí)驗(yàn)方法和過(guò)程

2.1材料

實(shí)驗(yàn)所用的聚苯乙烯標(biāo)準(zhǔn)顆粒由北京海岸鴻蒙標(biāo)準(zhǔn)物質(zhì)技術(shù)有限公司提供，共采用了3種樣品，利用光學(xué)顯微鏡分析，其顆粒圖像和粒度分布如圖2所示。測(cè)量結(jié)果顯示樣品1、樣品2和樣品3的體積中位徑分別為7.69、21.58、66.64 μm，同時(shí)可見(jiàn)粒度分布范圍較窄并且具有較好球形度。

圖2 顯微鏡法實(shí)驗(yàn)樣品圖像和顆粒粒度分布Fig.2　Image and particle size distribution of samples obtained by optical microscopy-based technique

2.2實(shí)驗(yàn)裝置和方法

圖3為實(shí)驗(yàn)中測(cè)量反射波譜裝置。裝置基本工作原理和其他超聲波譜儀類(lèi)似[18-19]。其中硬件包括超聲脈沖波發(fā)射接收儀、超聲換能器、懸浮液樣品池、雙通道高速數(shù)字化儀以及用于信號(hào)采集和分析的計(jì)算機(jī)。

圖3　懸浮液顆粒尺寸分布的超聲測(cè)量裝置Fig.3　Sketch of ultrasonic spectrometer used to measure particle size distribution in suspension

由前述數(shù)值模擬結(jié)果可知，不同顆粒粒徑具有對(duì)應(yīng)的特征頻率，實(shí)驗(yàn)中應(yīng)選用對(duì)應(yīng)中心頻率的超聲換能器。盡管商業(yè)上已經(jīng)出現(xiàn)頻率1～200 MHz的超寬頻掃描式超聲波發(fā)射和測(cè)量設(shè)備，但限于目前實(shí)驗(yàn)條件，本文針對(duì)樣品1、樣品2和樣品3分別選擇中心頻率為100 MHz (Model V2012 (BC), Panametrics)、30 MHz( Model V356-SU, Panametrics) 和10 MHz(Model V312-SU, Panametrics) 超聲換能器。

實(shí)驗(yàn)過(guò)程如下：在樣品池中置入一定量的聚苯乙烯懸浮液（20℃），混合均勻。由超聲換能器發(fā)出一束超聲波A，穿過(guò)緩沖層后，在緩沖層與樣品的界面處反射形成反射波B，再次經(jīng)過(guò)緩沖層后反射波被同一個(gè)超聲換能器接收，經(jīng)分析可獲得反射波B對(duì)應(yīng)的幅值（M）和相位(θ)。為計(jì)算反射系數(shù)，測(cè)出聲學(xué)特性已知的水的超聲反射信號(hào)作為背景，由式（5）計(jì)算反射系數(shù)Rs，該法可抵消換能器與緩沖層耦合處產(chǎn)生的聲能損失和聲場(chǎng)中聲發(fā)散造成的擴(kuò)散損失。

式（5）中下角標(biāo)c和s分別代表標(biāo)定物和測(cè)量樣品。利用已知的水和緩沖層的超聲特性，在溫度為20℃時(shí)，以石英玻璃緩沖層的聲阻抗為Zdl=1.2× 107kg·m?2·s?1，相對(duì)于水的反射系數(shù)為Rc=0.7778。將式(5)算得Rs代入式（6）中即可得到復(fù)阻抗。

Zs即為實(shí)驗(yàn)超聲復(fù)阻抗。并可按Zs=Zr+iZi分為實(shí)部Zr和虛部Zi。

3　實(shí)驗(yàn)現(xiàn)象和結(jié)果

3.1實(shí)驗(yàn)現(xiàn)象與分析

為了分析顆粒對(duì)于超聲波反射信號(hào)的影響，首先對(duì)分別安裝換能器于樣品池側(cè)面和底部(圖3)的兩類(lèi)樣品池布置方式對(duì)比測(cè)試，其主要差異在于超聲換能器安裝位置和顆粒壁面接觸形式。實(shí)驗(yàn)過(guò)程中，由于聚苯乙烯的密度大于水的密度，樣品池中顆粒會(huì)逐漸沉降并在底部累積，逐漸形成樣品池上部濃度逐漸降低，下部則增高情形。圖4以樣品池1為例示意了兩個(gè)不同時(shí)刻的顆粒分布，應(yīng)注意一部分聚苯乙烯顆粒會(huì)始終附著于壁面而并未下沉。

圖4　樣品池1內(nèi)樣品的沉降變化Fig. 4　Sedimentation variation of sample in measurement cell 1

對(duì)于樣品池1，從某一時(shí)刻起，按時(shí)間間隔60 s，記錄顆粒沉降過(guò)程中t1～t5時(shí)刻接收的超聲時(shí)域信號(hào)，由圖5(a)給出?？梢钥闯?，不同時(shí)刻信號(hào)彼此幾乎重合，但與背景信號(hào)存在一定差別。這說(shuō)明即便樣品池中顆粒濃度變化，所獲得的超聲時(shí)域信號(hào)并沒(méi)有明顯差異。而將附著于緩沖層與樣品界面處的顆粒去除之后T時(shí)刻，再次測(cè)量接收到的超聲信號(hào)BT[圖5(b)]，可見(jiàn)其與背景信號(hào)完全重合。由此可推斷測(cè)得聲阻抗信號(hào)與附著反射界面上的局部顆粒物有關(guān)，而與樣品區(qū)整體的顆粒粒徑和濃度并無(wú)關(guān)系。故后續(xù)顆粒表征實(shí)驗(yàn)中，均直接以附著在緩沖層表面的顆粒為測(cè)量對(duì)象，為此采用樣品池2進(jìn)行實(shí)驗(yàn)。同時(shí)由于測(cè)量區(qū)域只局限于緩沖層表層，與樣品池中其他區(qū)域的樣品無(wú)關(guān)，所以本文不討論峰的展寬；② 理論計(jì)算時(shí)事先設(shè)定了10%體積分?jǐn)?shù)，而實(shí)驗(yàn)信號(hào)由樣品池2獲得，其時(shí)樣品堆積于樣品的濃度測(cè)量。

圖5　樣品信號(hào)與背景信號(hào)的對(duì)比Fig.5　Comparison of signal of samples and background

3.2實(shí)驗(yàn)結(jié)果及分析

3種樣品的實(shí)驗(yàn)阻抗譜實(shí)部與數(shù)值模擬結(jié)果對(duì)比如圖6所示。實(shí)驗(yàn)結(jié)果與理論曲線趨勢(shì)基本一致，3種樣品的實(shí)驗(yàn)阻抗譜在理論特征頻率附近都出現(xiàn)了對(duì)應(yīng)的振蕩（最大值和最小值）。

圖6　3種樣品的實(shí)驗(yàn)阻抗譜與理論阻抗譜（實(shí)部）Fig. 6　Real parts of experimental and theoretical ultrasonic impedance

對(duì)比圖6(a)、(b)，雖然整體趨勢(shì)的一致性很好，但還是存在著略微的差別：實(shí)驗(yàn)的振蕩頻率區(qū)間（最大值與最小值對(duì)應(yīng)頻率差）較理論結(jié)果大；同時(shí)實(shí)驗(yàn)特征頻率均比理論值略小；實(shí)驗(yàn)結(jié)果的幅值明顯大于理論結(jié)果的幅值。上述差異的存在主要有如下幾個(gè)原因：① 實(shí)驗(yàn)用聚苯乙烯樣品顆粒粒徑具有一定分布，而理論值僅按其平均粒徑計(jì)算，造成特征樣品池底部致使此處顆粒濃度較高，使得振蕩幅度加大；③ 理論模型中并沒(méi)有考慮聲波在顆粒間的復(fù)散射，而實(shí)驗(yàn)中樣品局部濃度會(huì)很高，超聲與顆粒間的作用也更為復(fù)雜。

為了通過(guò)超聲波阻抗譜分析顆粒粒度分布，采用了一種求解病態(tài)方程的ORT算法[2]。該法是在Twomey方法的基礎(chǔ)上，利用廣義交叉校驗(yàn)（GCV）、L-Curve準(zhǔn)則[20]等方式對(duì)正則化因子優(yōu)化的一種獨(dú)立模式算法。為了用于阻抗譜的反演，在沒(méi)有改變核心求解過(guò)程的同時(shí)對(duì)系數(shù)矩陣計(jì)算和導(dǎo)入譜信息進(jìn)行了改動(dòng)。對(duì)于多分散顆粒系，參照衰減譜方法[16]將阻抗方程離散，并在多頻條件下矩陣方程AF=G。

式中，A為系數(shù)矩陣；F為離散化的顆粒尺寸頻度分布，是待求值；G為實(shí)際測(cè)得不同頻率下的超聲阻抗譜實(shí)部構(gòu)成的向量；下角標(biāo)i、j分別代表粒徑和阻抗譜的離散分檔數(shù)。

圖7為實(shí)驗(yàn)阻抗譜反演粒度分布曲線與圖像法結(jié)果的對(duì)比。由圖可知，超聲阻抗譜法很好區(qū)分3種不同粒徑區(qū)間的顆粒。粒度分布范圍和圖像分析結(jié)果也吻合較好，對(duì)應(yīng)表2給出的體積中位徑，樣品1的結(jié)果相對(duì)偏差最小。但隨著粒徑增大，阻抗法與圖像法的粒度分布寬度偏差加大，其體積中位徑與圖像分析結(jié)果也存在一定的偏差，但偏差小于11%。對(duì)于粒徑更大、粒徑范圍分布更廣的樣品，大顆粒的沉降較小顆?？?，所以其首先附著緩沖層界面，可能導(dǎo)致測(cè)量平均值偏大。

圖7　3種樣品的超聲阻抗譜法和圖像法測(cè)量結(jié)果的對(duì)比Fig. 7　Comparison between impedance and image measurement results for three samples

表2　超聲譜法與圖像法測(cè)量樣品體積中位徑Table 2　Comparison of sample volume median diameter measured by ultrasonic spectrum and image approach

4　結(jié)　論

通過(guò)阻抗定義對(duì)ECAH模型變形，進(jìn)行了顆粒系超聲波阻抗譜數(shù)值模擬，發(fā)現(xiàn)超聲阻抗譜實(shí)部對(duì)顆粒系特征變化敏感，特征頻率隨顆粒粒徑增大而遞減；特征幅度則隨濃度增加而遞增。設(shè)計(jì)一種底部入射式樣品池裝置測(cè)量并分析緩沖層與樣品界面反射波，通過(guò)獲得阻抗譜驗(yàn)證了數(shù)值結(jié)果。根據(jù)實(shí)驗(yàn)現(xiàn)象初步斷定：樣品的聲阻抗只與附著在緩沖層表面的顆粒特性有關(guān)，而與其他區(qū)域顆粒無(wú)關(guān)。利用ORT反演算法求得顆粒粒度分布，3種樣品粒度分布區(qū)間被很好地區(qū)分，與圖像法測(cè)量結(jié)果較為吻合，中位徑偏差小于11%。

為了發(fā)展超聲阻抗法并實(shí)現(xiàn)顆粒系的無(wú)損、非接觸在線測(cè)量，理論模型需要進(jìn)一步考慮聲復(fù)散射等高濃度效應(yīng)，實(shí)驗(yàn)中應(yīng)對(duì)顆粒運(yùn)動(dòng)和顆粒間相互作用對(duì)聲信號(hào)的影響加以研究。

符號(hào)說(shuō)明

A——入射波

An——壓縮波散射系數(shù)

B——反射波

c(ω)——超聲相速度，m·s?1

M——頻域信號(hào)幅值，V

n——顆粒分散等級(jí)

R——顆粒粒徑，m

Zs,Zdl——分別為介質(zhì)阻抗和緩沖層阻抗，kg·m?2·s?1

α(ω)——超聲衰減系數(shù)，Np·m?1

θ——相位

κ,κc——分別為超聲復(fù)波數(shù)、入射波壓縮波波數(shù)，m?1

λ——超聲波長(zhǎng)，m

ρ——懸浮液密度，kg·m?3

φ——懸浮液濃度

ω——超聲波角頻率，s?1

下角標(biāo)

c——標(biāo)定物

s——測(cè)量樣品

References

[1]WANG X Z, LIU L D, LI R F, et al. Online characterisation of nanoparticle suspensions using dynamic light scattering, ultrasound spectroscopy and process tomography [J]. Chemical Engineering Research and Design, 2009, 87(6): 874-884.

[2]SU M X, XU F, C X S, et al. Optimization of regularization parameter of inversion in particle sizing using light extinction method [J]. Particuology, 2007, 5(4): 295-299.

[3]BLAISOT J B, YON J. Droplet size and morphology characterization for dense sprays by image processing: application to the Diesel spray [J]. Exp. Fluids, 2005, 39(6): 977-994.

[4]CHEN X Z, ZHOU W, CAI X S, et al. In-line imaging measurementsof particle size, velocity and concentration in a particulate two-phase flow [J]. Particuology, 2014, 13: 106-113.

[5]成林虎, 蔡小舒, 周騖. 基于顆粒離焦模糊影像梯度算法的顆粒粒徑測(cè)量方法[J]. 化工學(xué)報(bào), 2012, 63(12): 3832-3838. CHENG L H, CAI X S, ZHOU W . Particle sizing from defocus image of spherical particles by image transition region gradient method [J]. CIESC Journal, 2012, 63(12): 3832-3838.

[6]劉海龍, 陳孝震, 蔡小舒, 等. 基于軌跡圖像的氣液旋風(fēng)分離器液滴粒度、濃度、速度的在線測(cè)量[J].化工學(xué)報(bào), 2012, 63(6): 1729-1734. LIU H L, CHEN X Z, CAI X S, et al. In-line measurement of size, concentration and velocity of drops from gas-liquid cyclone separator based on trajectory image processing[J]. CIESC Journal, 2012, 63(6): 1729-1734.

[7]Measurement and characterization of particles by acoustic methods(Ⅰ): Concepts and procedures in ultrasonic attenuation spectroscopy: ISO 20998-1—2006 [S].

[8]EPSTEIN P S, CARHART R R. The absorption of sound in suspensions and emulsions(Ⅰ): Water fog in air [J]. The Journal of the Acoustical Society of America, 1953, 25(3): 553-565.

[9]ALLEGRA J R, HAWLEY S A. Attenuation of sound in suspensions and emulsions: theory and experiments [J]. The Journal of the Acoustical Society of America, 1972, 51(5B): 1545-1564.

[10]SU M X, CAI X S, XUE M H, et al. Particle sizing in dense two-phase droplet systems by ultrasonic attenuation and velocity spectra [J]. Sci. China Ser. E-Technol. Sci., 2009, 52(6): 1502-1510.

[11]呼劍, 蘇明旭, 蔡小舒, 等. 高頻寬帶超聲衰減譜表征納米顆粒粒度的方法[J]. 化工學(xué)報(bào), 2010, 61(11): 2985-2991. HU J, SU M X, CAI X S, et al. Broad-band high-frequency ultrasonic attenuation spectrum method for measuring nanoparticle size distribution[J]. CIESC Journal, 2010, 61(11): 2985-2991.

[12]KULMYRZAEV A, CANCELLIERE C, MCCLEMENTS D J. Characterization of aerated foods using ultrasonic reflectance spectroscopy [J]. Journal of Food Engineering, 2000, 46(4): 235-241.

[13]KULMYRZAEV A, MCCLEMENTS D J. High frequency dynamic shear rheology of honey [J]. Journal of Food Engineering, 2000, 45(4): 219-224.

[14]KASOLANG S, AHMAD M A, DWYER J R S. Measurement of circumferential viscosity profile in stationary journal bearing by shear ultrasonic reflection [J]. Tribology International, 2011, 44(11): 1264-1270.

[15]DUKHIN A S, GOETZ P J. Fundamentals of interface and colloid science [M]//Studies in Interface Science. Amsterdam: Elsevier, 2002: 17-73.

[16]章維, 蘇明旭, 蔡小舒. 基于超聲衰減譜和相速度的顆粒粒徑測(cè)量[J]. 化工學(xué)報(bào), 2014, 65(3): 898-904. ZHANG W, SU M X, CAI X S. Particle size distribution measurement based on ultrasonic attenuation and phase velocity spectrum [J]. CIESC Journal, 2014, 65(3): 898-904.

[17]HAY A E, SCHAAFSMA A S. Resonance scattering in suspensions [J]. The Journal of the Acoustical Society of America, 1989, 85(3): 1124-1138.

[18]MCCLEMENTS D J, FAIRLEY P. Ultrasonic pulse echo reflectometer [J]. Ultrasonics, 1991, 29(1): 58-62.

[19]MCCLEMENTS D J, FAIRLEY P. Frequency scanning ultrasonic pulse echo reflectometer [J]. Ultrasonics, 1992, 30(6): 403-405.

[20]ELI?ABE G E, GARCíA R, LUIS H. Latex particle size distribution from turbidimetry using inversion techniques [J]. Journal of Colloid and Interface Science, 1989, 129(1): 192-200.

Particle size characterization based on ultrasonic impedance spectrum

JI Xiaoming, SU Mingxu, WANG Xue, CAI Xiaoshu
(Institute of Particle and Two-phase Flow Measurement, Shanghai Key Laboratory of Multiphase Flow and Heat Transfer in Power Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China)

Abstract:The purpose of this paper is to explore the relationship between the characteristics of particles and ultrasonic impedance spectra. Although the ultrasonic attenuation and velocity spectra methods are becoming the sophisticated technology for particle size characterization, they are not available in some actual projects, either because the particulate matters to be investigated are too dense or because there are some bubbles inside. However, ultrasonic impedance spectra technology emerges as a potential technique to provide a non-destructively detecting way in the aforementioned occasions. By modifying the ultrasonic attenuation spectral model, the theories related the ultrasonic impedance spectra to particle concentration and size are established to implement quite a few numerical simulations, showing that the acoustic impedances of particles is sensitive to the changes in particle concentration and size. After that, three kinds of polystyrene suspensions with different volume median diameters (7.69 μm, 21.58 μm and 66.64 μm) are investigated experimentally using the transducers with center frequencies varied from 10 MHz to 100 MHz. The results are quantitatively validated by the theory whereby a conclusion can be yielded that it is possible to distinguish the particles with different sizes by using ultrasonic impedance spectra. The comparison of the inversion results and image analysis indicates that the impedance technique is available to particle size characterization.

Key words:ultrasonic impedance spectra; polystyrene suspensions; particle size distribution; concentration; measurement; numerical simulation; inversion

中圖分類(lèi)號(hào)：TB 52；O 426.9

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：0438—1157（2016）06—2284—07

DOI：10.11949/j.issn.0438-1157.20151954

基金項(xiàng)目：國(guó)家自然科學(xué)基金項(xiàng)目（51176128,5120611）。

Corresponding author:Prof. SU Mingxu, sumx@usst.edu.cn