共生礦品位指標(biāo)的聯(lián)合優(yōu)化

1、<p>　　共生礦品位指標(biāo)的聯(lián)合優(yōu)化</p><p>　　摘 要：B選礦廠以共生礦為原料，經(jīng)弱磁、強(qiáng)磁、反浮選等多個(gè)階段，產(chǎn)出鐵精礦和稀土精礦（聯(lián)產(chǎn)品），本文針對(duì)此工藝進(jìn)行了聯(lián)合優(yōu)化，并給出操作建議。首先，在參數(shù)校準(zhǔn)的基礎(chǔ)上，通過對(duì)品位等指標(biāo)的統(tǒng)計(jì)分析建立系統(tǒng)模型；然后，選擇弱磁選的精礦品位和產(chǎn)率為關(guān)鍵操作參數(shù)，并根據(jù)歷史數(shù)據(jù)統(tǒng)計(jì)出其可行組合及概率分布；最后，建立優(yōu)化模型，并利用LINGO編程與求解。

2、</p><p>　　關(guān)鍵詞：品位；校準(zhǔn)；聯(lián)產(chǎn)品；聯(lián)合優(yōu)化 </p><p>　　近年來，受產(chǎn)能過剩等因素的影響，鋼鐵企業(yè)的利潤空間進(jìn)一步收窄，上游工序降低成本成為緩解這一狀況的主要舉措[1]。品位和產(chǎn)率是選礦工藝中顯著影響工藝成本的指標(biāo)，其優(yōu)化模型及求解方法成為研究熱點(diǎn)。李萬亨等創(chuàng)立了品位指標(biāo)的綜合優(yōu)化法，通過數(shù)學(xué)方法進(jìn)行技術(shù)經(jīng)濟(jì)分析來選出最佳的品位指標(biāo)方案[2]；袁懷雨等創(chuàng)立了整體動(dòng)態(tài)

3、優(yōu)化法，其特點(diǎn)是通過采-選系統(tǒng)建模，采用模糊綜合評(píng)判方法對(duì)邊界品位進(jìn)行多目標(biāo)優(yōu)化[3]；諸克軍等引入神經(jīng)網(wǎng)絡(luò)、遺傳算法等智能優(yōu)化方法，對(duì)截止品位和入選品位進(jìn)行優(yōu)化[4]。 </p><p>　　本文針對(duì)共生礦通過多個(gè)選別階段產(chǎn)出聯(lián)產(chǎn)品的工藝，在上述文獻(xiàn)的基礎(chǔ)上，需應(yīng)用聯(lián)合優(yōu)化方法對(duì)品位進(jìn)行優(yōu)化。Lamothe等研究了一種混合的整數(shù)線性規(guī)劃模型，來優(yōu)化包含多種產(chǎn)品的合成供應(yīng)鏈的運(yùn)作成本[5]。Xiuli Chao等

4、針對(duì)供給能力隨機(jī)的定期更新庫存系統(tǒng)，研究了其補(bǔ)給決策和定價(jià)決策的聯(lián)合優(yōu)化問題[6]。 </p><p>　　1 選礦過程模型的建立 </p><p>　　B選礦廠的入廠原礦是包含赤鐵礦、氟碳鈰等礦物的共生礦，產(chǎn)品為鐵精礦和稀土精礦，選礦工藝為弱磁—強(qiáng)磁—反浮選聯(lián)合工藝流程[7]，如圖1所示。 </p><p>　　圖1 B選礦廠選礦原則流程 </p>

5、<p><b>　　1.1 校準(zhǔn)模型 </b></p><p>　　B選礦廠的品位化驗(yàn)誤差為，不能滿足優(yōu)化要求，所以首先應(yīng)對(duì)生產(chǎn)數(shù)據(jù)進(jìn)行校準(zhǔn)。校準(zhǔn)是在規(guī)定的條件下，用一個(gè)可參考的標(biāo)準(zhǔn)，對(duì)包括參考物質(zhì)在內(nèi)的測(cè)量器具的特性賦值，并確定其示值誤差[8]。本文以質(zhì)量守恒和鐵量守恒為標(biāo)準(zhǔn)[9]，對(duì)弱磁、強(qiáng)磁、浮選三個(gè)選別階段分別進(jìn)行了校準(zhǔn)。限于篇幅，僅介紹強(qiáng)磁部分的校準(zhǔn)模型，其原料為弱磁選尾

6、礦（rw），選別產(chǎn)物有中磁選精礦（zj）、強(qiáng)磁精選精礦（qj）、強(qiáng)磁精選尾礦（qz）、強(qiáng)磁粗選尾礦（qw）四種。 </p><p>　　目標(biāo)函數(shù)： （1） </p><p><b>　　S.t （2） </b></p><p><b>　　（3） </b></p><p><b>　?。?

7、） </b></p><p><b>　?。?） </b></p><p><b>　?。?） </b></p><p><b>　?。?） </b></p><p>　　式中， 分別為重量和校準(zhǔn)后重量， 分別為品位和校準(zhǔn)后品位， 為校正量， 為 的最大測(cè)量誤差，

8、為校正系數(shù)， 為最大校正系數(shù)。約束1為質(zhì)量平衡，約束2為鐵平衡。 </p><p>　　1.2 數(shù)理統(tǒng)計(jì)模型 </p><p>　　通過對(duì)校準(zhǔn)后的生產(chǎn)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析[10]，得出了下列一元數(shù)據(jù)關(guān)系模型，以下數(shù)學(xué)模型皆已通過檢驗(yàn)。本研究選取2012年為優(yōu)化期，文中數(shù)據(jù)皆為優(yōu)化期內(nèi)B選礦廠的生產(chǎn)數(shù)據(jù)。 </p><p>　?。?）強(qiáng)磁粗選給礦鐵品位 與弱磁選尾礦鐵品位

9、 的數(shù)學(xué)模型 </p><p><b>　?。?） </b></p><p>　?。?）強(qiáng)磁精選尾礦鐵品位 與強(qiáng)磁粗選給礦鐵品位 的數(shù)學(xué)模型 </p><p><b>　　（9） </b></p><p>　?。?）強(qiáng)磁精選精礦鐵品位 與強(qiáng)磁粗選給礦鐵品位 的數(shù)學(xué)模型 </p><

10、;p><b>　?。?0） </b></p><p>　?。?）中磁選精礦鐵品位 與強(qiáng)磁精選精礦鐵品位 </p><p><b>　?。?1） </b></p><p>　?。?）強(qiáng)磁粗選尾礦鐵品位 與弱磁選尾礦鐵品位 的數(shù)學(xué)模型 </p><p><b>　?。?2） </b

11、></p><p>　?。?）強(qiáng)磁粗選尾礦稀土品位 與強(qiáng)磁粗選尾礦鐵品位 的數(shù)學(xué)模型 </p><p><b>　?。?3） </b></p><p>　?。?）強(qiáng)磁精選尾礦鐵品位 與弱磁選尾礦鐵品位 的數(shù)學(xué)模型 </p><p><b>　?。?4） </b></p><

12、p>　?。?）反浮選精礦產(chǎn)率 與反浮選給礦鐵品位 的數(shù)學(xué)模型 </p><p><b>　　（15） </b></p><p>　?。?）反浮選富集比 與反浮選給礦鐵品位 的數(shù)學(xué)模型 </p><p><b>　?。?6） </b></p><p>　　2 聯(lián)合優(yōu)化與求解 </p>

13、<p>　　2.1 優(yōu)化方案的制定 </p><p>　　共生原礦經(jīng)過碎磨后進(jìn)入弱磁選階段，弱磁選精礦經(jīng)反浮選產(chǎn)出鐵精礦，弱磁選尾礦經(jīng)強(qiáng)磁選產(chǎn)出稀土浮選的原料。弱磁選作為聯(lián)產(chǎn)品的分離點(diǎn)，其相關(guān)技術(shù)指標(biāo)的變化會(huì)對(duì)最終產(chǎn)品產(chǎn)生較大影響，因此，本文把弱磁選精礦的品位和產(chǎn)率作為關(guān)鍵操作參數(shù)。 </p><p>　　根據(jù)B廠操作規(guī)程和相關(guān)設(shè)備的運(yùn)行狀況，弱磁選精礦品位、產(chǎn)率指標(biāo)的調(diào)整范

14、圍分別限制在57%～63.5%和27%～46%之間，將以上兩個(gè)指標(biāo)變動(dòng)范圍構(gòu)成的區(qū)域劃分為若干小區(qū)域，并對(duì)每一個(gè)區(qū)域的概率進(jìn)行統(tǒng)計(jì)，確立了49個(gè)（弱磁精礦品位，弱磁精礦產(chǎn)率）組合方案，見表1。 </p><p>　　表1 優(yōu)化方案概率表 </p><p>　　2.2 聯(lián)合優(yōu)化模型 </p><p>　　本模型為聯(lián)合優(yōu)化，目標(biāo)函數(shù)中包括鐵精礦和稀土精礦兩種產(chǎn)品的銷售收

15、入。 </p><p><b>　?。?7） </b></p><p>　　式中，G為選礦廠總利潤，P1、P2、P3為鐵精礦、稀土精礦、原礦的基準(zhǔn)價(jià)格， 為鐵精礦、稀土精礦、原礦的單位基準(zhǔn)浮動(dòng)價(jià)格（品位每浮動(dòng)一個(gè)百分點(diǎn)浮動(dòng)價(jià)格）， 為鐵精礦、稀土精礦、原礦基準(zhǔn)品位，C為制造費(fèi)用，其它同前。 </p><p>　　模型的約束條件包括：（1）質(zhì)量守

16、恒和鐵量平衡式，見公式（2）～（3）；（2）統(tǒng)計(jì)式，見公式（8）～（16）；（3）可行操作范圍，其中弱精品位和產(chǎn)率的操作范圍，見表1.。（4）變量類型約束，本模型的變量均為大于零的實(shí)數(shù)。 </p><p><b>　　3 算例及結(jié)論 </b></p><p>　　采用LINGO編程，對(duì)本模型進(jìn)行規(guī)劃求解，優(yōu)化期內(nèi)鐵精礦價(jià)格最高和稀土精礦價(jià)格最高時(shí)的計(jì)算結(jié)果分別見圖2、

17、圖3，圖中A點(diǎn)為常用操作點(diǎn)，箭頭所指方向?yàn)槔麧櫶荻确较?，k為斜率。 </p><p>　　圖2 鐵精礦價(jià)格較高時(shí)利潤圖 </p><p>　　圖3 稀土精礦價(jià)格較高時(shí)利潤圖 　　通過對(duì)算例結(jié)果分析，可得出如下結(jié)論： </p><p>　?。?）針對(duì)現(xiàn)場(chǎng)測(cè)量誤差較大的情形，在進(jìn)行統(tǒng)計(jì)分析建模之前，先對(duì)數(shù)據(jù)進(jìn)行校準(zhǔn)，可減少稱量誤差及化驗(yàn)誤差，確保數(shù)理統(tǒng)計(jì)模型預(yù)測(cè)的合理

18、性； </p><p>　　（2）對(duì)于共生礦生產(chǎn)兩種精礦的生產(chǎn)工藝，當(dāng)兩種產(chǎn)品的價(jià)格波動(dòng)時(shí)，最優(yōu)（品位，產(chǎn)率）操作點(diǎn)會(huì)發(fā)生變化。就文中算例，當(dāng)鐵精礦價(jià)格較高時(shí)，品位和產(chǎn)率優(yōu)化的梯度方向?yàn)椋?：76.9）；當(dāng)稀土精礦價(jià)格較高時(shí)，品位和產(chǎn)率的優(yōu)化梯度方向?yàn)椋?：38.5）。 </p><p><b>　　參考文獻(xiàn)： </b></p><p>　　[

19、1]康志偉.中國鋼鐵工業(yè)應(yīng)對(duì)鐵礦石上漲策略分析[J].現(xiàn)代商貿(mào)工業(yè)，2008（11）：166～167. </p><p>　　[2]何煥學(xué)，李萬亨.品位指標(biāo)多方案優(yōu)化技術(shù)研究[J].中國地質(zhì)礦產(chǎn)經(jīng)濟(jì)，1995（11）：17～20. </p><p>　　[3] 諸克軍，郭海湘，高艷卉，王德銀.金屬礦山組合品位的智能系統(tǒng)優(yōu)化[J].系統(tǒng)管理學(xué)報(bào)，2010（11）：509～513. </

20、p><p>　　[4]李克慶，牛京考，袁懷雨，劉保順.白云鄂博鐵礦磁鐵礦石品位指標(biāo)的優(yōu)化[J].北京科技大學(xué)學(xué)報(bào)，2007，29（3）：334～337. </p><p>　　[5]Jacques Lamothe，Khaled Hadj-Hamou，Michel Aldanondo.An optimization model for selecting a product family and

21、 designing its supply chain[J].European Journal of Operational Research，2006（169）：1030–1047. </p><p>　　[6]Xiuli Chao，Hong Chen，Shaohui Zheng.Joint replenishment and pricing decisions in inventory systems wit

22、h stochastically dependent supply capacity[J].European Journal of Operational Research，2008（191）：142–155. </p><p>　　[7]伍喜慶，王振平，陳泉源.包鋼氧化礦選礦技術(shù)過程優(yōu)化[J].礦冶工程，1996，16（2）：33～36. </p><p>　　[8]唐超，李世平.非線性