李少華，朱明亮，張立棟，王 擎，劉朝青，于丁一

（1 中國大唐科學技術(shù)研究院，北京 100032；2 東北電力大學油頁巖綜合利用教育部工程研究中心，吉林 吉林 132012）

混合是物料在外力（重力或機械力等）作用下發(fā)生運動速度和運動方向的改變，使各組分顆粒均勻分布的操作過程；其目的在于通過混合過程獲得組成和性質(zhì)均勻的混合物，以保證所得產(chǎn)物組成、結(jié)構(gòu)和性能的均勻一致性[1-2]。

國內(nèi)外對顆粒混合的研究從20 世紀中葉開始至今一直沒有停息過。1966年Hogg 等[3]對水平攪拌混合裝置進行實驗研究并提出了混合動力方程；李水清等[4-5]利用分子動力學在中等轉(zhuǎn)速情況下經(jīng)受力分析模擬顆粒在水平回轉(zhuǎn)設(shè)備中的運動軌跡，并用PIV 試驗進行驗證。Ottino 小組[6-8]對顆粒物料的混合與分離進行了基本闡述，并提出建立顆粒物料運動模型所需考慮的基本問題與顆?；旌嫌嬎惴椒ǎㄟ^實驗對EDEM 等模擬軟件的模擬結(jié)果進行論證。Jha 等[9]研究了不同應(yīng)變速率物料對混合的影響，并得到當應(yīng)變速率為0.25 與0.5 時二元物料混合最好。Cleary 等[10]采用DEM 方法模擬了滾筒中的顆?；旌希⒉捎昧藘煞N不同方法估算混合度。Simo Siiri? 等[11]采用DEM 方法研究混合參數(shù)的選擇對混合效果的影響，得出選擇合適的混合參數(shù)，不但提高了混合的最終效果，而且在減少混合器耗損上也有著良好的效果。趙永志等[12-14]對薄滾筒內(nèi)二元S 型顆粒體系采用DEM 方法對顆粒流動進行數(shù)值模擬，討論顆粒在薄滾筒內(nèi)混合、運動及分布情況。

在這些研究中，顆?；旌隙鹊脑u價一直沒有一個定性的分析，雖然各學者均采用與自己研究相適合的混合度評價方式，但并無具體分析。本文作者以回轉(zhuǎn)裝置內(nèi)三組元顆粒徑向混合實驗為研究背景，對比較常用的混合度評價方法進行簡單介紹與 對比，并針對Lacey 指數(shù)取樣及計算相關(guān)問題進行探討。

1 混合度評價方法簡述

混合作為一個顆粒與顆粒間的物理行為，其混合過程極其復雜，對其混合程度的描述也難以詳盡。通常以定性描述的方式評價混合度，比較其混合優(yōu)劣。圖1 是經(jīng)常被用來描述顆粒的3 種混合狀態(tài)。圖1(a)為完全分離狀態(tài)，又稱為顆粒混合初始狀態(tài)；圖1(b)為部分混合狀態(tài)，又稱為隨機混合狀態(tài)，一般實際混合狀態(tài)就是這一狀態(tài)，兩種顆粒呈無序不規(guī)則排列方式；圖1(c)為完全混合狀態(tài)，又稱為理想混合狀態(tài)，兩種顆粒完全均勻地分布在一個顆粒群里，且互相之間接觸最大而同種顆粒間無接觸，但這種狀態(tài)一般在生產(chǎn)中無法達到。

為了能夠定量的分析顆粒間的混合程度，以便直觀地了解其混合程度的高低，通常采用標準差、 變異系數(shù)、接觸數(shù)、混合指數(shù)等來評價混合程度。

圖1 顆?；旌夏Ｐ?/p>

1.1 標準差評價法

該方法主要利用多組樣本數(shù)據(jù)的標準偏差σ 來評價顆?；旌铣潭?，是一種較為簡單的評價方法，但標準偏差的大小與測定次數(shù)有關(guān)。另外，標準偏差值與各測定值相對于平均值的殘差有關(guān)，而與各測定值本身的大小無關(guān)。當混合物料中的組分含量相差懸殊時，用標準偏差很難說明混合程度。因此僅用標準偏差還不足以充分說明混合程度。

1.2 變異系數(shù)評價法

變異系數(shù)（Cv）又稱離散度，為了客觀地反映混合程度，在上述標準偏差方法的基礎(chǔ)上對平均值這一參數(shù)突顯出來，用以彌補標準偏差法的不足。其算法為在上述標準偏差的基礎(chǔ)上再利用式（1）進行計算。

變異系數(shù)能夠較客觀地反映出顆?；旌锨闆r，曾被應(yīng)用于測量回轉(zhuǎn)裝置內(nèi)顆粒軸向混合度[15]。但對于徑向混合來說，如果對混合各時刻圖像進行網(wǎng)格取樣，計算所得混合受局部影響很大，如本實驗混合穩(wěn)定后利用變異系數(shù)法計算得52 s 時Cv=1.6，而在54 s 時Cv=0.874，其值變化極大，不能夠正確地反映出徑向混合變化規(guī)律。

1.3 接觸數(shù)評價法

Van Puyvelde 等[16]通過對顆?；旌蠄D像的研究提出了通過計算顆粒接觸數(shù)來衡量顆粒混合程度的方法，這種方法能計算出不同區(qū)域內(nèi)的混合程度。由圖1 可以得出，(a)接觸數(shù)為5，(b)接觸數(shù)為23，(c)接觸數(shù)為31，通過這組接觸數(shù)可以很好地反映出不同時刻特定區(qū)域內(nèi)顆粒的混合狀態(tài)，在區(qū)域(a)、(b)、(c)內(nèi)接觸數(shù)有很大的區(qū)別，顆粒接觸數(shù)呈現(xiàn)出遞增的趨勢，混合度也呈現(xiàn)出遞增的趨勢，采用接觸數(shù)評價混合度能很好地反應(yīng)混合程度。

張立棟等[17]采用大顆粒和小顆粒之間的接觸數(shù)與總的接觸數(shù)（大顆粒與大顆粒，大顆粒與小顆粒，小顆粒與小顆粒）的比值來衡量混合程度q，其數(shù)學表達式如式（2）所示。

式中，Csl為大顆粒與小顆粒的接觸數(shù)；Ctotal為總的接觸數(shù)。q 值越小表明混合效果越差，越大表明混合效果越好。

但這種方法對顆粒統(tǒng)計要求極高，很難在實驗分析中實現(xiàn)，因此多數(shù)被用于數(shù)值模擬中的顆?；旌显u價。

1.4 Lacey 指數(shù)評價法

該方法也可稱為混合指數(shù)法，是Lacey[18]在1954年利用統(tǒng)計學方法提出的一種針對顆?；旌隙鹊挠嬎惴椒?。他通過對顆粒混合標準差評價法進行分析，結(jié)合顆?；旌蠈嶋H情況，提出混合指數(shù)計算式（3）。

式中，S2為兩種顆粒實際混合方差；S02為完全分離時的混合方差；Sr2為兩種顆粒完全隨機混合方差。各量計算公式如式（4）。

式中，Ns為樣本總數(shù)；ai為兩類顆粒任意一類在樣本i 中所占的比例；為相應(yīng)顆粒在滾筒內(nèi)所占顆粒比例；ki為樣本i 的權(quán)重，可表示為式（5）；k 表示為式（6）。

式中，Ni為樣本i 內(nèi)的顆粒數(shù)；Nt為顆粒數(shù) 總和。

完全分離時的混合方差S02計算公式為式（7）。

式中，P 為一種顆粒所占比例；（1-P）為另一種顆粒所占比例。

兩種顆粒完全隨機混合方差Sr2計算公式為 式（8）。

式中，n 為每一個樣本內(nèi)平均顆粒數(shù)。

該方法與變異系數(shù)法有著共同的優(yōu)點，但在取樣上差別比較明顯，其取樣采用對徑向圖像網(wǎng)格劃分取樣法，可利用攝像等設(shè)備對實驗進行連續(xù)取樣。

1.5 其它評價法

在混合度研究上國內(nèi)外學者一直沒有停止其前進的腳步，隨著對顆?；旌险J識的深入，對顆?；旌系乃惴ㄒ苍谥饾u發(fā)展著。Finnie 等[19]用統(tǒng)計焓來衡量混合度，這種方法注重對某一區(qū)域內(nèi)物料比率的計算。Van Puyvelde[20]研究表明，用接觸數(shù)和統(tǒng)計焓評價混合度分別適用于不同的領(lǐng)域。采用統(tǒng)計焓這種方法重點衡量某一區(qū)域內(nèi)的總體混合，而接觸數(shù)則用于發(fā)展互相混合的顆粒間的換熱模型。

通過上文對比發(fā)現(xiàn)，混合程度評價方法較合理且應(yīng)用較廣泛的有變異系數(shù)法、接觸數(shù)法及Lacey指數(shù)算法，其在混合度評價上各有優(yōu)劣，且適用范圍不同，其中，Lacey 指數(shù)算法比較適合評價徑向混合；變異系數(shù)比較適合評價軸向混合，也可應(yīng)用于易取料的工業(yè)混合測量；變異系數(shù)法主要適用于計算機數(shù)值模擬中涉及的混合。因此本文作者應(yīng)用Lacey 指數(shù)算法對顆?；旌线M行評價。

2 Lacey 指數(shù)算法探討

為了適應(yīng)不同粒徑顆粒群混合度的計算，國內(nèi)外學者普遍采用等效顆粒數(shù)算法對獲得的大小顆粒數(shù)進行處理，主要方式就是設(shè)一標準顆粒，之后將其它顆粒利用等效關(guān)系轉(zhuǎn)換成標準顆粒數(shù)，如本文中，以1 mm 直徑顆粒為標準顆粒，由于只研究徑向即二維顆?；旌?，故可以將1 個3 mm 直徑顆?；? mm 直徑顆粒按面積等效法轉(zhuǎn)換成9 個標準顆粒（1 mm 顆粒）或25 個標準顆粒，再利用標準顆粒總數(shù)代入式（3）～式（8）中進行計算。

Lacey 指數(shù)算法在應(yīng)用過程中可分為取樣、提取數(shù)據(jù)、計算混合度三部分，其中提取數(shù)據(jù)與計算混合度對于徑向混合度并無太大影響，只要能夠保證其準確性即可，但在取樣上對于顆?；旌嫌绊懕容^大。Lacey 指數(shù)算法取樣如圖2 所示，其取樣樣本多少將直接影響到樣本內(nèi)顆粒分布，進而影響總體混合度。因此一種合適的取樣將影響顆粒混合評價正確與否。本文作者以33.3%填充率無傾角3.4 r/min 工況為例，分別采用4×4、6×6、8×8、 10×10 共4 種取樣方式，分析不同取樣對顆?；旌系挠绊懬闆r。為了全面地研究取樣對顆?；旌隙鹊挠绊?，取混合初始時0～8 s 及混合穩(wěn)定后的52～60 s 兩個階段進行分析。

圖2 Lacey 指數(shù)取樣示例

圖3 混合初始階段各取樣混合度

2.1 混合數(shù)值變化分析

初始0～8 s 各個取樣方法混合度如圖3 所示，隨著樣本的增加，混合度逐漸降低，但混合度變化趨于平緩。其中4×4 取樣方式下，顆?；旌隙鹊淖兓顬閯×遥m然總體變化趨勢還是上升的，但上下波動幅度最大；而10×10 取樣方式下顆?；旌隙茸兓顬槠骄?，整體呈現(xiàn)出先略降后緩慢上升的趨勢；6×6 取樣方式及8×8 取樣方式所獲得的混合度變化介于4×4 和10×10 兩者之間。

為了獲得引起這一變化的原因，將取樣原圖及Lacey 指數(shù)算法與混合度曲線進行對比分析。

首先按Lacey 指數(shù)算法分析。Lacey 指數(shù)計算混合度M 直接涉及的量有S2、S02、Sr2三個，因此針對這3 個量進行分析：其中，根據(jù)式（7）可以得出S02只與大小顆?？倲?shù)有關(guān)，與取樣方式無關(guān)；Sr2可看作S02/n，其中n值為顆?？倲?shù)除以樣本個數(shù)，在本文中n 值均較大，求出的Sr2均為10-3數(shù)量級，與S02（10-1數(shù)量級）相差很大，對M 值影響極小，故假設(shè)引起圖3 現(xiàn)象與Sr2無直接關(guān)系；根據(jù)式（8）得出S2與單個樣本內(nèi)顆粒總是有直接關(guān)系，因此假設(shè)S2是引起圖3 現(xiàn)象的主要原因。

為了驗證上述兩點假設(shè)，對實驗數(shù)據(jù)作以下分析：用顆粒總數(shù)代替n，這樣n 將不再受取樣方式影響，計算發(fā)現(xiàn)其規(guī)律與圖3 所示規(guī)律一致且數(shù)據(jù)基本無變化，這表明假設(shè)一正確；如圖4 所示，對比 4×4 取樣與8×8 取樣0 s 時部分圖片，發(fā)現(xiàn)其取樣空間相同，切割外框線相同，僅是將4×4 中每一樣本再次分成4 個樣本，但其整體混合度值相差很大，其4×4 時M 值為0.299，而8×8 時M 值為0.1496，這表明樣本內(nèi)顆粒數(shù)是引起混合度變化的關(guān)鍵。由于Sr2假設(shè)成立，所以S2假設(shè)也是成立的。

圖4 4×4 取樣與8×8 取樣對比

在S2計算公式中，k 在本文計算中為恒定值1，為小顆粒在滾筒內(nèi)所占顆粒比例，即q，并不會因樣本取法不同而改變，故這兩個量不予分析。因此引起S2變化的只可能與ki、ai有關(guān)。

ki與ai在計算過程中都需要樣本內(nèi)等效顆粒數(shù)及樣本內(nèi)大小顆粒數(shù)，通過圖4 能夠看出，4×4 取樣與8×8 取樣其每一樣本內(nèi)顆粒數(shù)均有極大變化，特別是8×8 中47、54 與55 三個樣本內(nèi)部包含了1 mm、 3 mm、5 mm 顆粒，這將直接影響ki與ai兩個值，從而影響到S2，這也是不同取樣所得混合度不同的根本原因，即單一樣本內(nèi)大小顆粒數(shù)比例及個數(shù)改變引起混合度數(shù)值的改變。

為了證明同一樣本內(nèi)顆粒分布是影響混合度數(shù)值變化的根本原因，本文對2 s 時刻圖像進行如下處理：在3.4 r/min 轉(zhuǎn)速下2 s 時間內(nèi)可使?jié)L筒旋轉(zhuǎn)40.8°，但這一角度并未引起顆粒間相對位置較大改變，所以2 s 時刻與0 s 時刻混合度理論值應(yīng)該基本相同，因此將2 s 時刻圖像順時針（滾筒為逆時針旋轉(zhuǎn)）旋轉(zhuǎn)40.8°，獲得圖像與2 s 時刻初始位置及0 s 時刻圖像對比如圖5 所示，圖5(b)與圖5(c)顆粒分布基本無變化，按4×4 取樣計算圖5(b)混合度為0.3026，這與圖5(c)即0 s 時刻混合度0.2992基本相同，但與圖5(a)混合度0.4878 相差甚遠。

對比其4×4 取樣圖（圖6）發(fā)現(xiàn)，在大小顆粒交匯位置樣本內(nèi)顆粒數(shù)有很大不同，圖6(a)中這一交匯位置單一樣本內(nèi)大小顆粒數(shù)相差不多，而圖6(b)與圖6(c)這一位置單一樣本內(nèi)大小顆粒數(shù)相差極大，這就導致這一區(qū)域單一樣本內(nèi)混合度前者要遠遠優(yōu)于后者，從而引起整體混合度之間的巨大 差異。

圖5 滾筒旋轉(zhuǎn)對比圖

圖6 滾筒旋轉(zhuǎn)取樣對比圖

通過上述討論可得，Lacey 指數(shù)計算中取樣方式對顆?；旌隙葦?shù)值影響主要體現(xiàn)在單一樣本內(nèi)顆粒分布上，單一樣本內(nèi)顆粒個數(shù)相差越小，則整體混合值越高。

取上述工況混合穩(wěn)定階段52～60 s 圖像進行計算得混合度如圖7 所示。隨著取樣數(shù)量的增加，其混合度值降低，但其總體趨勢變化不大。表明取樣不同對混合趨勢基本沒有影響，但對混合度數(shù)值影響較大，這與初始階段所呈現(xiàn)的規(guī)律相同。這表明上述分析結(jié)論適用于混合各個階段。

2.2 取樣標準分析

從上文分析結(jié)果可以看出，雖然取樣方式對顆?；旌锨€趨勢影響較小，但對混合度數(shù)值影響極大。通過相關(guān)性計算得：取樣網(wǎng)格大小與顆?；旌隙缺憩F(xiàn)為顯著性正相關(guān)，其相關(guān)系數(shù)為0.9。因此如何選取取樣方式將影響實驗數(shù)據(jù)處理的準確性。

通過圖6 及相關(guān)結(jié)論可以得出，取樣空間過大并不利于顆粒混合度的準確分析，由于取樣空間過大，其內(nèi)部易產(chǎn)生大小顆粒群間混合，被誤認為顆粒間混合進行計算，這使得計算結(jié)果產(chǎn)生較大誤差。

圖7 混合穩(wěn)定階段各取樣混合度

圖8 取樣引起混合變低示意圖

但取樣空間過小時，會產(chǎn)生大小顆粒被人為的分離現(xiàn)象，如圖8 暗色標記顆粒附近，其上部大顆粒與小顆粒完全混合在一起，但由于網(wǎng)格過小無法將大顆粒附近的全部信息包含在一起，使得其混合度被人為的降低。針對本實驗工況計算得：所用顆粒為1 mm、3 mm、5 mm，其大小顆粒理想化均勻橫向混合為1 mm 顆粒、3 mm 顆粒，1 mm 顆粒、5 mm 顆粒、1 mm 顆粒分布模式，其橫向尺寸為11 mm，則單一樣本橫向尺寸最好在11 mm 左右，即可將其信息全部包含在內(nèi)。本研究所取樣本4×4、6×6、8×8、10×10 計算成單一樣本單邊尺寸分別為21 mm、14 mm、10.5 mm 及8.4 mm，則14 mm樣本與10.5 mm 樣本接近最佳樣本，但10.5 mm 樣本容易引起圖8 所示問題，故14 mm 取樣方式即 6×6 取樣方式為最佳。通過對比圖3 與圖7 發(fā)現(xiàn)，相對于其它兩種取樣方式而言，這兩種取樣方式各項數(shù)據(jù)變化趨勢與其它相同而數(shù)值居于中間位置，初始階段兩種取樣方式下混合度變化基本相同，但穩(wěn)定階段8×8 取樣方式下顆粒混合度上下變化較大，相對而言6×6 取樣方式優(yōu)于8×8 取樣方式。其結(jié)果驗證了理想混合度計算選取的取樣方式的正確性。

為了進一步驗證此結(jié)論對三組元顆粒徑向混合評價上的廣泛性，對6.8 r/min、10 r/min 及16.7%填充率3.4 r/min 工況穩(wěn)定階段進行計算得如圖9 所示結(jié)果。結(jié)果表明，不同取樣方式所引起的混合度變化基本相同。這一結(jié)果表明上述取樣方式普遍適用于三組元顆粒混合實驗。

3 結(jié) 論

對比分析多種混合度評價方法，較常用的混合度評價方法為變異系數(shù)、接觸數(shù)及Lacey 指數(shù)，其中變異系數(shù)比較適合評價軸向混合，也可應(yīng)用于易取料的工業(yè)混合測量；接觸數(shù)法主要適用于計算機數(shù)值模擬中涉及的混合；Lacey 指數(shù)算法比較適合評價徑向混合。

圖9 混合度綜合取樣

通過對Lacey 指數(shù)算法取樣方式對比分析可知，在同一工況下不同取樣方式所得混合度有很大差異，其原因在于不同取樣方式引起單一樣本內(nèi)大小顆粒比例不同，從而引起S2值產(chǎn)生較大變化，最終引起M 值變化；在實驗處理上最佳取樣尺寸應(yīng)該略大于顆粒最佳混合時橫向尺寸；本實驗所得最 佳取樣尺寸為單一樣本單邊長14 mm，即6×6 取樣方式。

[1] 盧壽慈. 粉體技術(shù)手冊[M]. 北京：化學工業(yè)出版社，2004：590-603.

[2] Williams J C. The segregation of Powders and granular materials[J]. Fuel Soeiety Journal，1963，14：29-35.

[3] Hogg R，Cahn D S，Healy T W，et al. Diffusional mixing in an deal system[J]. Chemical Engineering Science，1966，21（11）：1025-1038.

[4] Li Shuiqing，Yao Qiang，Chen Bing，et al. Molecular dynamics simulation and continuum modeling of granular surface flow in rotating drums[J]. Chinese Science Bulletin，2007，52（5）：692-700.

[5] 李水清，姚強，陳冰，等. 回轉(zhuǎn)筒顆粒表面流的分子動力學模擬和連續(xù)理論[J]. 科學通報，2006，51（19）：2319-2326.

[6] Ottino J M，Khakhar D V. Mixing and segregation of granular materials[J]. Annual Review of Fluid Mechanics，2000，32：55-91.

[7] Ottino J M，Lueptow R M. On mixing and demixing[J]. Science，2008，319（5865）：912-913.

[8] Pérez-Alonso C，Delgadillo J A. Experimental validation of 2D DEM code by digital image analysis in tumbling mills[J]. Minerals Engineering，2012，25（1）：20-27.

[9] Jha A K，Puri V M. Percolation segregation of multi-size and multi-component particulate materials[J]. Powder Technology，2010，197（3）：274-282.

[10] Cleary P W，Metcalfe G，Liffman K. How well do discrete element granular flow models capture the essentials of mixing processes?[J]. Applied Mathematical Modelling，1998，22（12）：995-1008.

[11] Siiri? S，Yliruusi J. DEM simulation of influence of parameters on energy and degree of mixing[J]. Particuology，2011，9（4）：406-413.

[12] 趙永志，程易. 水平滾筒內(nèi)二元顆粒體系徑向分離模式的數(shù)值模擬研究[J]. 物理學報，2008，57（1）：321-328.

[13] 趙永志，張憲旗，劉延雷，等. 滾筒內(nèi)非等粒徑二元顆粒體系增混機理研究[J]. 物理學報，2009，58（12）：8386-8393.

[14] Jiang Maoqiang，Zhao Yongzhi，Liu Gesi，et al. Enhancing mixing of particles by baffles in a rotating drum mixer[J]. Particuology，2011，9（3）：270-278.

[15] 張立棟，李少華，余侃勝，等. 油頁巖與固體熱載體在回轉(zhuǎn)干餾爐內(nèi)混合特性的冷態(tài)實驗[J]. 化工進展，2011，30（3）：492-497.

[16] Van Puyvelde D R，Young B R，Wilson M A，et al. Experimental determination of transverse mixing kinetics in a rolling drum by image analysis[J]. Powder Technology，1999，106（3）：183-191.

[17] 李少華，張立棟，張軒，等. 回轉(zhuǎn)干餾爐內(nèi)影響顆?；旌线\動因素的數(shù)值分析[J]. 中國電機工程學報，2011，31（2）：32-38.

[18] Lacey B P M C. Developments in the theory of particle mixing[J]. Journal of Applied Chemistry，1954，4：257-268.

[19] Schutyser M A I，Padding J T，Weber F J，et al. Discrete particle simulations predicting mixing behavior of solid substrate particles in a rotating drum fermenter[J]. Biotechnology and Bioengineering，2001，75（6）：666-675.

[20] Van Puyvelde D R. Comparison of discrete elemental modeling to experimental data regarding mixing of solids in the transverse direction of a rotating kiln[J]. Chemical Engineering Science，2006，61（13）：4462-4465.