蘭 川，董宏軍，劉佳鵬，徐敬元，金永新，田全寶，陳天浩

北方礦業(yè)有限責(zé)任公司 北京 100053

J KSimMet 是澳大利亞昆士蘭大學(xué) JKMRC 礦物研究中心基于 60 多年的研究基礎(chǔ)開(kāi)發(fā)的一款選礦廠碎磨回路穩(wěn)態(tài)模擬軟件，主要應(yīng)用于選礦流程計(jì)算、碎磨設(shè)備選型和碎磨流程優(yōu)化領(lǐng)域。選礦流程計(jì)算功能是一個(gè)獨(dú)立的軟件，整合到 JKSimMet 和 JKSimFloat中，用來(lái)調(diào)整選礦試驗(yàn)和選礦廠流程考察的數(shù)據(jù)。碎磨設(shè)備選型是 JKSimMet 應(yīng)用最多的地方，在設(shè)備選型過(guò)程中調(diào)用了 JKSimMet 的仿真引擎，通過(guò)對(duì)不同的流程結(jié)構(gòu)[1-2]或不同的磨機(jī)直徑和長(zhǎng)度進(jìn)行仿真運(yùn)算[3-4]，從而選擇最佳的流程結(jié)構(gòu)和設(shè)備規(guī)格。然而，新建項(xiàng)目在設(shè)備選型階段，礦石碎磨性質(zhì)參數(shù)、設(shè)備參數(shù)和工藝操作參數(shù)不詳盡，但為了設(shè)備選型，通常會(huì)采用附近礦山或性質(zhì)相似礦山的技術(shù)指標(biāo)，或者采用軟件的預(yù)設(shè)值。對(duì)于在產(chǎn)選礦廠的仿真優(yōu)化而言，這種處理無(wú)法通過(guò)有效的模型擬合來(lái)發(fā)揮仿真軟件的預(yù)測(cè)功能，不能對(duì)優(yōu)化措施提供量化指導(dǎo)。

筆者在對(duì)某 S AB 流程詳細(xì)考察和落重試驗(yàn)的基礎(chǔ)上[5]，用 JKSimMet 軟件對(duì)剛果(金)某大型銅鈷礦D6400×L3500 半自磨機(jī)—圓筒篩回路進(jìn)行了模型擬合；對(duì)裝球率、處理量、礦石硬度、頑石口尺寸和數(shù)量等因素進(jìn)行了仿真試驗(yàn)；并制作了一套頑石口孔隙為 60 mm 的格子板，逐步替換 30 mm 的格子板，驗(yàn)證了相關(guān)仿真試驗(yàn)的結(jié)論。

1 研究方法

JKSimMet 軟件以運(yùn)行中的選礦廠作為研究對(duì)象，通過(guò)流程考察全面了解現(xiàn)有選廠的運(yùn)行狀況，在此基礎(chǔ)上選擇合適的模型和參數(shù)，通過(guò)模型擬合使軟件再現(xiàn)生產(chǎn)廠實(shí)際運(yùn)行狀況，然后改變操作參數(shù)、設(shè)備參數(shù)，預(yù)測(cè)回路的運(yùn)行指標(biāo)。作為仿真工具使用時(shí)，按照畫(huà)流程圖、平衡數(shù)據(jù)、模型選擇和數(shù)據(jù)輸入、模型擬合和模擬試驗(yàn) 5 個(gè)過(guò)程開(kāi)展研究，具體工作流程如圖1 所示。

圖1 JKSimMet 軟件模擬仿真流程Fig.1 Simulation process of JKSimMet software

1.1 工藝流程

SAB 碎磨工藝流程為破碎—半自磨—球磨—旋流器分級(jí)。破碎合格的礦石經(jīng)礦倉(cāng)給入半自磨機(jī)，半自磨機(jī)排礦經(jīng)圓筒篩分級(jí)，篩上物由返料輸送帶返回半自磨機(jī)，篩下物由泵送至旋流器進(jìn)行分級(jí)，旋流器底流作為球磨機(jī)給料，溢流作為磨礦分級(jí)最終產(chǎn)品送至下游濕法工序。球磨機(jī)出料和半自磨機(jī)篩下物料共同匯入泵池，由泵給入旋流器進(jìn)行分級(jí)，形成磨礦和檢查分級(jí)閉路流程[5]。半自磨機(jī)—圓筒篩回路設(shè)計(jì)工藝流程如圖2 所示。主要設(shè)備為半自磨機(jī)和圓筒篩，輔助設(shè)備為礦倉(cāng)、給礦水設(shè)備、返料輸送帶，料流以帶箭頭的線條表示。

圖2 半自磨機(jī)—圓筒篩回路工藝流程Fig.2 Process flow of circuit of SAG mill with trommel

1.2 質(zhì)量平衡

質(zhì)量平衡是將流程考察中獲得的原始數(shù)據(jù)輸入軟件，根據(jù)質(zhì)量平衡原理對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行調(diào)整的過(guò)程。JKSimMet 軟件把料流數(shù)據(jù)分為試驗(yàn)數(shù)據(jù)(標(biāo)記為Exp)、平衡數(shù)據(jù)(標(biāo)記為 Bal)、仿真數(shù)據(jù)(標(biāo)記為 Sim)3 種類型。

1.3 模型選擇和數(shù)據(jù)輸入

數(shù)學(xué)模型是計(jì)算機(jī)仿真的基礎(chǔ)，不同模型需要輸入的參數(shù)和能夠預(yù)測(cè)的參數(shù)也不一樣。JKSimMet 6.3 版的半自磨機(jī)有 4 個(gè)模型，分別是 V ariable r ates AG/SAG Mill、Semi-Autogenous Mill(Leung)、Ball Mill(Perfect Mixing)和 Size Convertor。后 2 個(gè)模型分別是球磨機(jī)和碎磨設(shè)備的通用模型，但在半自磨機(jī)仿真試驗(yàn)中的作用有限，半自磨機(jī)選用 Variable rates AG/SAG Mill 模型(簡(jiǎn)稱變速模型)。該模型以 Semi-Autogenous Mill(Leung)為基礎(chǔ)發(fā)展形成，允許合適的料流返回半自磨機(jī)，并拓展了磨機(jī)直徑的放大范圍，突破了原來(lái) 8～9 m 的上限。篩上物由返料輸送帶返回半自磨機(jī)，滿足變速模型的要求。

輸入的數(shù)據(jù)除磨機(jī)內(nèi)部尺寸、轉(zhuǎn)速、料位(裝球率和總充填率)、礦石的可磨性、進(jìn)出料的粒度外，還包含排礦端格子板開(kāi)孔的大小和位置。變速模型的磨機(jī)參數(shù)輸入窗口如圖3 所示。從上到下數(shù)據(jù)依次為：直徑(襯板間)，m；筒體長(zhǎng)度(襯板間)，m；進(jìn)料端中空軸直徑，m；進(jìn)料端錐角(平頭=0)，(°)；出料端錐角(平頭=0)，(°)；格子孔尺寸-XG，mm；細(xì)粒尺寸 -XM，mm；格子板開(kāi)口面積分?jǐn)?shù)；頑石口占總開(kāi)口面積分?jǐn)?shù)；頑石口尺寸 -XP，mm；格子孔徑向平均相對(duì)位置。

圖3 半自磨機(jī)尺寸參數(shù)輸入窗口Fig.3 Size parameter input window for SAG mill

1.4 模型擬合

為半自磨機(jī)和圓筒篩選擇合適的模型，輸入流程考察得到的礦石參數(shù)、設(shè)備參數(shù)、料流參數(shù)，進(jìn)行模型擬合。模型擬合是應(yīng)用尋優(yōu)算法，尋找合適的參數(shù)值代替缺省值，使軟件預(yù)測(cè)的技術(shù)指標(biāo)最大程度地再現(xiàn)流程考察條件下的技術(shù)指標(biāo)。

半自磨機(jī)—圓筒篩回路的模型擬合結(jié)果如圖4所示。比較各料流的試驗(yàn)值和擬合值，半自磨機(jī)的總充填率為 20.96%，與流程考察實(shí)測(cè)值 21% 基本符合；物料的總密度為 3.32 t/m3，與流程考察實(shí)測(cè)值3.27 t/m3比較接近；綜合比較回路中各產(chǎn)品的礦石流量、P80和磨機(jī)工作指標(biāo)(功率、充填率)，模型預(yù)測(cè)指標(biāo)與流程考察指標(biāo)符合較好。結(jié)果表明，該模型擬合的參數(shù)較好地反映了設(shè)備工作狀況和礦石性質(zhì)，可以用于仿真試驗(yàn)。

模型擬合得到半自磨機(jī)變速模型的破碎速率如圖5 所示。破碎速率可以理解為單位時(shí)間內(nèi)某粒級(jí)礦石經(jīng)歷的破碎事件的次數(shù)。從圖5 可以看出：0.1 mm 以下的細(xì)顆粒破碎速率幾乎為零，說(shuō)明細(xì)粒級(jí)礦石幾乎不經(jīng)過(guò)破碎就直接排出；1～10 mm 之間的礦石破碎速率有一個(gè)峰值，說(shuō)明這個(gè)區(qū)間的礦石最適合半自磨機(jī)破碎；破碎速率有一個(gè)波谷，這是頑石容易出現(xiàn)的粒度區(qū)間；60 m m 以上的礦石破碎速率快速增大，說(shuō)明大塊礦石強(qiáng)度小，容易破碎。

圖5 半自磨機(jī)變速模型破碎速率Fig.5 Breakage rates of variable speed model of SAG mill

1.5 仿真模式

JKSimMet 軟件的過(guò)程仿真有兩種模式：一種是單一場(chǎng)景的仿真試驗(yàn)，即改變一個(gè)或多個(gè)操作條件，預(yù)測(cè)全回路的各項(xiàng)指標(biāo)；另一種是多因素-多水平的批量仿真試驗(yàn)，允許改變的操作參數(shù)較多，如磨機(jī)轉(zhuǎn)速、裝球率、磨機(jī)尺寸、排礦格子板開(kāi)孔的大小和位置、給礦量、給礦粒度、礦石硬度等。筆者分別進(jìn)行了半自磨機(jī)裝球率仿真試驗(yàn)，半自磨機(jī)裝球率、處理量?jī)梢蛩?多水平仿真試驗(yàn)和頑石口尺寸、數(shù)量與礦石硬度仿真試驗(yàn)。

2 仿真試驗(yàn)

2.1 半自磨機(jī)裝球率仿真試驗(yàn)

該試驗(yàn)為單因素-多水平仿真試驗(yàn)，給礦量設(shè)定為 200 t/h，考察半自磨機(jī)的裝球率 5%～15% 時(shí)，設(shè)備參數(shù)(半自磨機(jī)功率、總充填率)和料流參數(shù)(圓筒篩篩上物料、篩下物料)的情況。仿真試驗(yàn)結(jié)果如表1 所列。

表1 半自磨機(jī)裝球率仿真試驗(yàn)結(jié)果Tab.1 Simulation experiment results of ball loading rate of SAG mill

由表1 可知，增加鋼球充填率，磨機(jī)功率增加，總充填率逐步降低，產(chǎn)品粒度變細(xì)，返砂量減少，返砂粒度稍微減?。谎b球率為 5% 時(shí)，半自磨機(jī)總充填率達(dá) 29.3%，超過(guò) J KSimMet 軟件平均總充填率為25% 的限定條件[6]，可能會(huì)導(dǎo)致磨機(jī)“脹肚”；當(dāng)返砂量為 48.92 t/h，也同樣超過(guò)設(shè)計(jì)限值；增加鋼球充填率，產(chǎn)品粒度緩慢變細(xì)，但磨機(jī)功率增加較快，生產(chǎn)運(yùn)營(yíng)的經(jīng)濟(jì)性差。推薦的半自磨機(jī)鋼球充填率為7%～8%，總充填率在 23%～25% 之間，生產(chǎn)運(yùn)營(yíng)經(jīng)濟(jì)性相對(duì)較好。

2.2 半自磨機(jī)裝球率、處理量?jī)梢蛩?多水平仿真試驗(yàn)

為了探索通過(guò)增加半自磨機(jī)裝球率提高半自磨機(jī)處理量的可能性，進(jìn)行了 4 個(gè)裝球率和 5 個(gè)處理量組合的 20 個(gè)仿真試驗(yàn)。裝球率是鋼球體積占磨機(jī)容積的百分比；處理量百分比是以生產(chǎn)回路的處理量(200 t/h)為基點(diǎn)進(jìn)行變化的百分比；回路的實(shí)際處理量是圓筒篩篩下的礦量。試驗(yàn)考察了半自磨機(jī)的功率和總充填率、圓筒篩篩下物料的粒度和流量、圓筒篩篩上物料的粒度和流量等 6 個(gè)參數(shù)，結(jié)果如表2 所列。

設(shè)立適當(dāng)?shù)臉?biāo)準(zhǔn)，判斷仿真試驗(yàn)的處理量是否實(shí)現(xiàn)。首要限定條件為 25% 的總充填率；另一判定條件為返砂量，一般按新給礦量的 20%～30% 選擇返砂輸送帶和頑石破碎機(jī)。該試驗(yàn)中，返砂輸送帶是按新給礦量的 20% 設(shè)計(jì)，即返砂量需低于 40 t/h。由表2 可知：依據(jù) 2 個(gè)判定條件，在 8% 的裝球率下，合理的處理量為 220 t/h(仿真試驗(yàn) 12)；在 10% 的裝球率下，合理的處理量為 260 t/h(仿真試驗(yàn) 19)；在12% 的裝球率下，合理的處理量為 300 t/h(仿真試驗(yàn)26)；裝球率為 15% 時(shí)，處理量可以達(dá)到 360 t/h(仿真試驗(yàn) 31)，此時(shí)半自磨機(jī)的功率為 2 035 k W，接近電動(dòng)機(jī)的額定功率 2 100 kW。

2.3 頑石口尺寸、數(shù)量與礦石硬度仿真試驗(yàn)

格子板是半自磨機(jī)的排礦通道。格子板上開(kāi)孔的大小、數(shù)量和徑向位置會(huì)影響半自磨機(jī)的排礦粒度和排礦速度。格子板的開(kāi)孔量為開(kāi)孔面積與磨機(jī)橫截面積的比值，開(kāi)孔位置用格子孔至磨機(jī)軸心的距離與磨機(jī)半徑之比的加權(quán)平均值表示。磨機(jī)格子板開(kāi)孔量為9.38%，徑向平均位置 0.74，格子孔為長(zhǎng)方形，兩端呈半圓形，孔寬為 30 mm。

半自磨機(jī)內(nèi)的頑石量與礦石硬度、粒度、裝球率、球徑、襯板外形、磨機(jī)轉(zhuǎn)速等多種因素相關(guān)，在襯板的一個(gè)生命周期內(nèi)，變化最大的是襯板外形和礦石硬度。因此，實(shí)際的頑石量是不斷變化且難以控制的。該試驗(yàn)中，在礦山的整個(gè)生命周期內(nèi)，礦石硬度A×b為 43.0，屬于中偏硬礦石，有可能產(chǎn)生頑石；流程考察中，礦石硬度A×b為 61.2，屬于中軟礦石，返砂量實(shí)測(cè)為 8.34 t/h。對(duì)礦體不同深度礦石的可磨性測(cè)試表明，礦石硬度隨開(kāi)采深度逐漸變大。為了應(yīng)對(duì)未來(lái)可能的頑石量增多問(wèn)題，設(shè)計(jì)了頑石外排方案，利用輸送帶返回半自磨機(jī)，并預(yù)留了頑石破碎機(jī)的安裝空間。

JKSimMet 軟件的變速模型包含了頑石口的尺寸、位置和數(shù)量，可以模擬頑石口尺寸和位置對(duì)磨礦指標(biāo)的影響。試驗(yàn)分別模擬頑石口為 40、50 和 60 mm 3 種情況，替換量為 16 組格子板的 1/4，頑石口在外圈和內(nèi)圈交叉配置，均勻排出頑石，如圖6 所示。

圖6 1/4 的格子板開(kāi)設(shè)頑石口的布置方案Fig.6 Layout scheme for setting hard stone outlet on 1/4 grid plate

在裝球率為 8%、處理量為 200 t/h 條件下，3 種頑石口尺寸(40、50 和 60 mm)和 2 種礦石硬度(A×b=61.2，A×b=43.0)，共 6 種條件組合，對(duì)頑石量影響的試驗(yàn)結(jié)果如圖7 所示。由圖7 可以看出，頑石口階梯變化時(shí)，頑石口每增大 10 mm，頑石量增加 10%左右；同時(shí)礦石硬度變大后，相對(duì)于每個(gè)頑石口，頑石量平均增加 66%，由此可知，礦石硬度對(duì)頑石量的影響更為顯著。

圖7 礦石硬度和頑石口尺寸對(duì)頑石量的影響Fig.7 Influence of ore hardness and size of hard stone outlet on amount of hard stone

2.4 仿真結(jié)果的初步驗(yàn)證

JKSimMet 軟件模型是在大量實(shí)驗(yàn)室試驗(yàn)和工業(yè)生產(chǎn)指標(biāo)統(tǒng)計(jì)的基礎(chǔ)上逐步形成的[7]，仿真結(jié)果的試驗(yàn)驗(yàn)證是增強(qiáng)預(yù)測(cè)結(jié)果的可信度和軟件持續(xù)更新的基礎(chǔ)。為了驗(yàn)證引入頑石口對(duì)頑石量的預(yù)測(cè)結(jié)果，初步將 1 塊 30 mm 的外圈格子板替換成了 60 mm 的格子板(見(jiàn)圖8 中 1 點(diǎn)鐘位置)。測(cè)驗(yàn)結(jié)果表明，返砂量由8.34 t/h 增加至 10.00 t/h，返砂P80也由 26.34 mm 增加至 33.48 mm，該變化趨勢(shì)與仿真試驗(yàn)的預(yù)測(cè)一致。

圖8 裝有 1 塊 60 mm 格子板的 B 系列半自磨機(jī)排礦端Fig.8 Discharge end of B series SAG mill equipped with a 60 mm grid plate

3 結(jié)論

利用前期半自磨機(jī)—圓筒篩回路流程考察結(jié)果，對(duì)半自磨機(jī)—圓筒篩回路進(jìn)行模型擬合，得到了相應(yīng)的模型參數(shù)；運(yùn)用 JKSimMet 軟件，進(jìn)行了裝球率、處理量、頑石口和礦石硬度等參數(shù)的仿真試驗(yàn)，得出以下結(jié)論。

(1)半自磨機(jī)功率隨裝球率的增加而增加。半自磨機(jī)可通過(guò)增加鋼球充填率提高處理量，在處理量為200 t/h 的目標(biāo)下，鋼球充填率在 7%～8% 范圍內(nèi)，半自磨機(jī)可以穩(wěn)定運(yùn)行，控制性良好。

(2)裝球率、處理量?jī)梢蛩?多水平仿真試驗(yàn)表明，在總充填率保持 25% 的限制條件下，裝球率為12% 時(shí)，處理量可達(dá) 300 t/h。此時(shí)，裝球率、處理量均提高 50%，而磨機(jī)功率只提高了 24.1%。這為提高磨機(jī)電能利用效率指明了方向。

(3)頑石口仿真試驗(yàn)和工業(yè)實(shí)踐表明，頑石口的引入可以增加頑石排出量，但是礦石硬度對(duì)頑石產(chǎn)量的影響更大。隨著礦山采礦深度的增加，礦石硬度越來(lái)越大，當(dāng)這種趨勢(shì)影響到預(yù)期的產(chǎn)量時(shí)，磨機(jī)內(nèi)會(huì)產(chǎn)生更多的頑石。因此，安裝更多的頑石口乃至頑石破很有必要，屆時(shí)可為驗(yàn)證仿真結(jié)果提供更多的生產(chǎn)數(shù)據(jù)。