鄢發(fā)明， 艾光華， 吳彩斌， 李曉東， 石志中， 毛文明

（1.江西理工大學(xué)江西省礦業(yè)工程重點實驗室，江西 贛州341000；2.湖南柿竹園有色金屬有限責(zé)任公司，湖南 郴州423037）

磨礦作業(yè)在選礦廠中占有非常重要的地位，其能耗約占整個選礦廠能耗的50%左右[1].磨礦是一個復(fù)雜的高能耗、低效率的作業(yè)，其過程還不是完全清楚，為了更好地進行磨礦作業(yè)，研究磨礦動力學(xué)很有必要[2-4].目前，對礦石磨碎的動力學(xué)模型有許多種，主要分為2類，一類是粒度及其相關(guān)量與時間及其相關(guān)量的動力學(xué)模型，另一類是總體平衡動力學(xué)模型[5].本文選用E.W.戴維斯提出的動力學(xué)模型來分析討論鎢礦石在磨礦過程中的行為，此模型研究的是被磨物料磨碎速率與其磨礦時間的關(guān)系，而物料的磨碎速率遵循磨礦動力學(xué)方程，磨礦動力學(xué)方程的參數(shù)k和n在一定程度上可以反映物料被磨碎的快慢以及為選擇性磨礦和現(xiàn)場生產(chǎn)指標評價提供依據(jù)[6-15].磨礦動力學(xué)方程式還可以用來預(yù)測不同磨礦時間磨礦產(chǎn)品的粒度分布[16-18].在很多情況n階磨礦動力學(xué)較一階磨礦動力學(xué)更符合實際，因此大多情況下，此動力學(xué)模型均是用n階磨礦動力學(xué)來分析礦石在磨礦過程中的行為[19-21]，而本文針對湖南柿竹園鎢礦石特性，將其磨礦產(chǎn)品分為不大于0.1 mm粒級的細粒級產(chǎn)品和大于0.1 mm粒級的粗粒級產(chǎn)品，分別探討其磨礦動力學(xué)規(guī)律并建立數(shù)學(xué)模型，為選礦廠磨礦實踐提供指導(dǎo).

1 試 驗

1.1 礦石性質(zhì)

試驗所用鎢礦石取自湖南柿竹園多金屬選廠，主要礦物為白鎢礦、黑鎢礦、假象及半假象黑鎢礦，其中白鎢礦與黑鎢礦的比例為7∶3.白鎢礦密度為6.1 t/m3，解離度平行中等，普氏硬度系數(shù)f為4.5～5，屬中等硬度；黑鎢礦密度為7.30 t/m3，解離度平行完全，普氏硬度系數(shù)f為4～5.5，屬中等硬度.鎢礦石不規(guī)則抗壓強度為400 kg/cm2左右，即換算成普氏硬度系數(shù)f為4左右.

1.2 試驗條件

在XMB三輥四筒磨機上分別進行了鎢礦石不同磨礦時間（2min、3min、4min、5min、6min、7min）的分批磨礦試驗.所用磨筒尺寸為φ135mm×146mm，磨筒轉(zhuǎn)速為85 r/min，磨礦介質(zhì)為鋼球，磨礦濃度為70%，介質(zhì)充填率為45%，磨礦量為200 g.

1.3 粒度測定

采用篩孔尺寸為 0.074 mm、0.1 mm、0.15 mm、0.3 mm、0.5 mm、1.0 mm、2.0 mm 的套篩測定物料0.074～2.00mm間的各粒級分布；采用激光粒度分析儀測定物料≤0.074 mm的各粒級分布.

1.4 試驗方法

對不同磨礦時間下的磨礦產(chǎn)品粒度分布，用一階磨礦動力學(xué)進行驗證，并采用Origin軟件進行擬合，求解不同粒級下磨礦動力學(xué)參數(shù)k值，觀察并比較求解不同粒級的磨礦動力學(xué)參數(shù)k值時的擬合程度R2，確定磨礦產(chǎn)品符合一階磨礦動力學(xué)的粒級范圍，并建立一階磨礦動力學(xué)方程式[22-24].對于磨礦產(chǎn)品不在一階磨礦動力學(xué)粒級范圍內(nèi)的，同樣采用Origin軟件進行求解不同粒級的磨礦動力學(xué)參數(shù)k和n值，并在此粒級范圍內(nèi)建立n階磨礦動力學(xué)方程式[25].最后分別運用所求得的各自粒級范圍內(nèi)的一階或n階磨礦動力學(xué)方程式進行理論計算，將計算結(jié)果與實際磨礦結(jié)果進行比較，來驗證本研究所求解的磨礦動力學(xué)方程式的準確性，為選礦廠磨礦提供理論依據(jù)[26-28].

2 磨礦動力學(xué)原理

磨礦動力學(xué)是指被磨物料的磨碎速率與磨礦時間的關(guān)系，其n階磨礦動力學(xué)方程式為：

式（1）中：R為磨礦為t時刻磨礦產(chǎn)品中粗粒級產(chǎn)率；R0為原礦中該粗粒級比率；t為磨礦時間；k為動力學(xué)參數(shù)；n為動力學(xué)參數(shù)，當n=0時稱為零階動力學(xué)，n=1時稱為一階動力學(xué)[12，29].

對式（1）進行變換得：

式（1）中，動力學(xué)參數(shù)n和k分別是粒度d的函數(shù)，表達式如式（4）和式（5）所示[30]：

式（4）、式（5）中 c0，c1，a0，a1及 x1，x2分別為待定參數(shù).

3 結(jié)果與討論

3.1 一階磨礦動力學(xué)

3.1.1 磨礦動力學(xué)參數(shù)k的求解

在磨礦條件為磨筒轉(zhuǎn)速170 r/min，礦漿濃度為70%，介質(zhì)充填率為45%下，依次對物料進行不同時長（2 min、3 min、4 min、5 min、6 min、7 min）的磨礦，結(jié)果如表1所示.

由表1可知，0.1 mm及其以下粒級篩上累積產(chǎn)率隨著時間的延長逐漸減少，且趨勢較平緩；0.1mm以上粒級篩上累積產(chǎn)率在5 min之前減少的較快，而5min之后減少的較慢.這表明粗粒級礦粒相對于較細粒級礦粒更好磨，隨著磨礦時間的延長，各類粒級礦石裂縫、缺陷越少，磨細越困難，因而使得礦石的可磨性逐漸變差[8].

表1 不同磨礦時間下磨礦產(chǎn)品篩上累積粒度分布/%Table 1 Cum ulative oversize distribution of grinding products under different grinding tim e/%

由表1可計算出磨礦產(chǎn)品不同時間不同粒級的ln（R0/R）值，結(jié)果如表 2 所示.

表2 不同磨礦時間下不同粒級的ln（R0/R）值Tab le 2 ln（R0/R） of different particle sizes under different grinding tim e

根據(jù)表2繪制不同粒級的ln（R0/R）與時間t的關(guān)系曲線，如圖1所示.并根據(jù)式（1）對不同粒級的ln（R0/R）與時間t的關(guān)系曲線進行線性擬合，得到鎢礦石不同粒級下的磨礦動力學(xué)參數(shù)k值，結(jié)果如表3所示.

圖1 不同粒級下ln（R0/R）與時間t的關(guān)系曲線Fig.1 Relationship curve between the ln（R0/R） values of different particle sizes and grinding tim e t

表3 不同粒級的磨礦動力學(xué)參數(shù)kTable 3 G rinding kinetics parameters k of different particle sizes

由表3可知：隨著粒級的減小，k值逐漸變小，礦粒被磨碎速率越小[5]，這是因為隨著礦石粒徑的減少，剩余礦粒表面裂縫相對較少，強度較高.礦粒在0.1mm及其以下粒級的線性相關(guān)系數(shù)R2均在0.99以上，而在0.1mm以上粒級的線性相關(guān)系數(shù)R2均在0.96以上，擬合程度不高.因此礦粒在0.1mm及其以下粒級符合一階磨礦動力學(xué)，而在0.1mm以上粒級則符合n階磨礦動力學(xué).

3.1.2 一階磨礦動力學(xué)方程的建立

在上述基礎(chǔ)上，將表3中不同粒級的磨礦動力學(xué)參數(shù)k在0.1mm及其以下粒級進行曲線擬合，求出鎢礦石的磨礦動力學(xué)參數(shù)k的表達式，結(jié)果如圖2所示.

圖2 磨礦動力學(xué)參數(shù)k與d的關(guān)系曲線Fig.2 Relationship curve between the grinding kinetics param eters k and d

由圖2可知，參數(shù)Ki呈冪函數(shù)的形式，可表示為：

利用Origin軟件對式（6）進行曲線擬合，求得參數(shù)a0、a1和x，從而求得鎢礦石的一階磨礦動力學(xué)參數(shù)k的表達式為：

將式（7）代入式（1），可建立鎢礦石在 0.1 mm 及其以下粒級的一階磨礦動力學(xué)方程為：

3.2 n階磨礦動力學(xué)

3.2.1 磨礦動力學(xué)參數(shù)n和k的求解

由上述分析可知，鎢礦石在0.1 mm以上粒級時不符合一階磨礦動力學(xué).為驗證其是否符合n階磨礦動力學(xué)，由表1中0.1 mm以上粒級數(shù)據(jù)通過式（3）進行計算，得到不同粒級的lg（lg R0/R）與t的關(guān)系，結(jié)果如表4所示.

表4 不同磨礦時間下不同粒級lg（lg R0/R）值Tab le 4 lg（lg R0/R） of different particle sizes under different grinding tim e

根據(jù)表4繪制不同粒級的 lg（lg R0/R）與 lg t的關(guān)系曲線，如圖 3所示.根據(jù)式（2）對不同粒級的lg（lg R0/R）與lg t的關(guān)系曲線進行線性擬合，得到不同粒級的動力學(xué)參數(shù)n值，進一步利用式（3）計算動力學(xué)參數(shù)k值，結(jié)果如表5所示.

圖3 不同粒級的lg（lg R0/R）與lg t的關(guān)系曲線圖Fig.3 Relationship curve between the lg（lg R0/R）values of different particle sizes and lgt

表5 不同粒級的磨礦動力學(xué)參數(shù)n和kTable 5 G rinding kinetics param eters n and k of different particle sizes

由表5可以看出，隨著粒級的減小，礦粒的動力學(xué)參數(shù)n值逐漸增大，k值逐漸減小.這是因為隨著礦石粒徑的減少，礦粒表面裂縫相對越少，礦粒強度越高，因而表現(xiàn)為隨著礦石粒徑逐漸減小，礦石的可磨性逐漸變差直至變化不大；隨著礦石粒徑的減少，礦石逐漸變脆，礦石的均勻性逐漸變差[16].

3.2.2 n階磨礦動力學(xué)方程的建立

根據(jù)表5中數(shù)據(jù)繪制動力學(xué)參數(shù)n和k與粒徑d的關(guān)系曲線，并進行擬合，如圖4所示.

圖4 動力學(xué)參數(shù)n和k與粒徑d的關(guān)系曲線圖Fig.4 Relationship curve between grinding kinetics param eters n and k and diam eter d

由圖4可知，曲線k=f（d）呈冪函數(shù)形式，而曲線n=f（d）更接近于對數(shù)形式， 如式（9）和式（10）所示.

根據(jù)表5數(shù)據(jù)，利用Origin軟件分別對式（9）和式（10）進行曲線擬合，求得未知參數(shù) c0，c1，a0，a1和x，進而求得鎢礦石的磨礦動力學(xué)參數(shù)k和n的表達式分別為：

將式（11）和式（12）代入式（1）得在 0.1 mm 以上粒級的n階磨礦動力學(xué)方程為：

3.3 磨礦動力學(xué)方程的驗證

為了驗證所建立的磨礦動力學(xué)方程對柿竹園鎢礦石磨礦的適用性，將式（8）和式（13）建立的0.1 mm及其以下粒級與0.1 mm以上粒級的磨礦動力學(xué)方程進行計算，得到各自粒級范圍內(nèi)不同磨礦時間下的磨礦產(chǎn)品篩上累積產(chǎn)率理論值，結(jié)果如表6所示.

表6 不同磨礦時間下的磨礦產(chǎn)品篩上累積產(chǎn)率理論值/%Tab le 6 Theoreticalvaluesof cum ulative oversize distribution of grinding p roducts under different grinding time/%

將表6數(shù)據(jù)與表1數(shù)據(jù)進行對比，得到利用磨礦動力學(xué)方程計算的不同磨礦時間下磨礦產(chǎn)品不同粒級的篩上累積產(chǎn)率理論值與實際磨礦試驗值的偏差值，結(jié)果如表7所示.

表7 篩上累積產(chǎn)率理論值與實際磨礦試驗值的偏差值△R/%Tab le 7 Deviation values△R of theoretical values and actual grinding test values of cum u lative oversize/%

從表7可以看出，所建立的磨礦動力學(xué)方程對柿竹園鎢礦石實際磨碎規(guī)律的吻合程度較高，與其它磨礦動力學(xué)模型相比，該磨礦動力學(xué)方程相對結(jié)構(gòu)簡單.因此，利用所建立的磨礦動力學(xué)方程為選礦廠預(yù)測磨礦產(chǎn)品粒度分布及進一步為通過調(diào)節(jié)磨機給礦量和分級機返砂量來控制磨礦產(chǎn)品粒度分布情況提供了理論依據(jù).

4 結(jié) 論

1）對柿竹園鎢礦石而言，礦粒在0.1 mm及其以下粒級符合一階磨礦動力學(xué)，在0.1mm以上粒級則符合n階磨礦動力學(xué).

2）一階磨礦動力學(xué)參數(shù)k與粒徑d呈冪函數(shù)關(guān)系，其一階磨礦動力學(xué)方程為R=R0exp[-（0.014 42+9.538 65d1.76128）t].

3）n階磨礦動力學(xué)參數(shù)中，k與粒徑d呈冪函數(shù)關(guān)系，n與粒徑d呈對數(shù)函數(shù)關(guān)系，其n階磨礦動力學(xué)方程為 R=R0exp[-（-0.134 1+1.151 42d0.68169）t0.71818-0.38656lnd].

4）所建立的磨礦動力學(xué)方程相對結(jié)構(gòu)簡單，而且所得到的磨礦產(chǎn)品篩上累積產(chǎn)率理論值與磨礦試驗值吻合程度較高，可以用來指導(dǎo)選礦廠磨礦實踐.

[1]馬天雨，桂衛(wèi)華，王雅琳，等.基于分批試驗的工業(yè)球磨機粒級分布預(yù)測模型[J].北京工業(yè)大學(xué)學(xué)報， 2012，38 （9）：1281-1286.

[2]T?nshoff H K， Friemuth T， Becker J C.Process monitoring in grinding[J].CIRP Annals-Manufacturing Technology，2002，51（2）：551-571.

[3]Mizonov V E，Berthiaux H，Zhukov V P，et al.Application of multi-dimensional markov chains to model kinetics of grinding with internal classification[J].International Journal of Mineral Processing， 2004，74 （增刊 1）：307-315.

[4]劉全軍，姜美光.碎礦與磨礦技術(shù)發(fā)展及現(xiàn)狀[J].云南冶金，2012，41（5）：21-28.

[5]王志輝.建陶原料粉磨特性的試驗研究[J].福建建材，2004（1）：17-19.

[6]段希祥.磨礦動力學(xué)參數(shù)與磨礦時間的關(guān)系研究[J].昆明工學(xué)院學(xué)報，1988，13（5）：23-33.

[7]侯英，丁亞卓，印萬忠，等.磨礦動力學(xué)參數(shù)對磨礦速度的影響[J].東北大學(xué)學(xué)報，2013，34（5）： 708-711.

[8]李茂林，崔瑞，王非，等.階段磨選磨礦粒度劃分的理論分析與計算[J].礦冶工程，2011，31（1）：27-32.

[9]劉全軍，唐榮，文書明.石英-磁鐵礦體系選擇性磨細研究[J].昆明工學(xué)院學(xué)報，1994，19（3）：85-89.

[10]田金星.石墨及其混合物料的磨礦動力學(xué)行為[J].中國有色金屬學(xué)報，1996，6（4）：47-50.

[11]田金星，楊華，尚保平.石墨-方解石混合物料方解石組分的碎裂參數(shù)特性[J].金屬礦山，2000（2）：41-43.

[12]葉賢東，文書明，張文彬.超臨速磨機磨礦動力學(xué)研究[J].有色金屬，2004，56（1）： 83-85.

[13]吳明珠，張成名.磨礦動力學(xué)在磨礦操作中的應(yīng)用[J].有色金屬，1989（5）：10-15.

[14]侯英，印萬忠，丁亞卓，等.不同破碎方式下產(chǎn)品磨礦特性的對比研究[J].有色金屬，2014（1）：5-8.

[15]印萬忠，吳凱，王紀鎮(zhèn)，等.破碎方式對紫金山銅金礦石可磨性及浮選的影響[J].福州大學(xué)學(xué)報，2014，42（2）： 321-326.

[16]王力，張常法，張軍，等.助磨劑對煤瀝青磨礦動力學(xué)的影響[J].山東科技大學(xué)學(xué)報，2008，27（5）：23-26.

[17]劉曉，徐超，李波.煤瀝青水漿干法制備的磨礦試驗研究[J].廣東化工，2011，39（23）：82-84.

[18]王舒婭，龍光明，祁米，等.磨礦時間對天青石粒度和鍶質(zhì)量分數(shù)分布的影響[J].化工礦物與加工，2011（7）： 8-10.

[19]王敏，丁浩.攪拌磨濕法超細磨礦過程動力學(xué)的研究[J].中國非金屬礦工業(yè)導(dǎo)刊，1999（2）：21-23.

[20]陳炳辰，周凌嘉，仲崇波.邦德功指數(shù)模擬計算方法的研究[J].金屬礦山，1990（8）：36-39.

[21]仲崇波，王成功，陳炳辰.邦德棒磨功指數(shù)模擬計算方法的研究[J].中國礦山，1997，6（2）：48-52.

[22]Ahmadi R，Hashemzadehfini M，Amiri Parian M.Rapid determination of bond rod-mill work index bymodeling the grinding kinetics[J].Advanced Powder Technology，2013（24） ：441-445.

[23]Bozkurt V，?zgür I.Dry grinding kineticsof colemanite[J].Powder Technology， 2007，176：88-92.

[24]Ipek H.The effects of grinding media shape on breakage rate[J].Minerals Engineering， 2006（19）：91-93.

[25]Ian Bailon-Poujol， Jean-Paul Bailon， Gilles L’Espérance.Ballmill grinding kinetics ofmaster alloys for steel powdermetallurgy applications[J].Powder Technology， 2011，210：267-272.

[26]Heekyu Choi， Woong Lee， Seongsoo Kim.Application of grinding kinetics analysis of inorganic powders by a stirred ball mill[J].Korean Journal of Chemical Engineering， 2009，26 （6）：1806-1812.

[27]Fuerstenau DW，De A，Kapur P C.Linear and nonlinear particle breakage processes in comminution[J].International Journal of Mineral Processing， 2004，74：317-327.

[28]Teke E，Yekeler M，Ulusoy U，et al.Kinetics of dry grinding of industrialminerals：calcite and barite[J].International Journal of Mineral Processing， 2002，67（1）：29-42.

[29]段希祥.碎礦與磨礦[M].北京：冶金工業(yè)出版社，2012.

[30]陳丙辰.磨礦原理[M].北京：冶金工業(yè)出版社，1989.