《武汉工程大学学报》 2012年10期
11-14,18
出版日期：2012-11-06
ISSN:1674-2869
CN:42-1779/TQ
磷矿堆场多目标优化配矿模型
0引言配矿是矿物资源质与量的相互协调与均衡，是一项多阶段的工艺过程，它贯穿于从矿山规划、计划到矿山生产整个过程和各个生产环节.堆场配矿是这一整体系统的重要环节之一，特别是对像磷矿矿石特征各异，原矿品质差异大，脉石对选矿、磷化工影响严重的一类矿产资源，堆场配矿具有举足轻重的地位\[1\].对配矿计划数学模型的分析和研究，国内外有很多研究成果可以借鉴.其解决方式也多种多样，比较成熟的技术如线性规划\[2-9\]、0-1整数规划法\[10\]、UM模型\[11\]、非线性回归分析等等.其中胡清淮\[6\]在建立线性规划模型基础上，对模型退化性进行了讨论，并借助计算机用逆矩阵单纯形法完成了求解计算；戚国安和舒航\[7\]对应用线性规划模型配矿中的随机误差进行了分析，通过计算期望产值，确定理想配矿品位，并指出配矿质量控制对矿山经济效益有着重要影响.张立溥和胡清淮\[8\]分析了配矿模型中的目标函数系数灵敏度的重要性.刘文生应用模糊线性规划的方法\[9\]，分析了矿山企业生产过程中包含许多非确定性的因素和模糊特征的主体网络系统.以上线性规划等配矿构模一般以经济效益为目标，对产能、有益有害组分、非负性等进行约束，建立模型求解.而在矿山实际配矿过程中，配矿目标往往不只一项，且很多情况下，各子目标可能是相互冲突的，非多目标规划模型就显得无能为力.众多学者尝试了多目标规划模型\[12-16\]通过引入模糊方法等，对目标重要性处理开展了大量研究\[17-19\]，取得一定成果.然而，常规的多目标解算方法不能有效地解决诸多配矿实际问题及要求\[20\]，模糊多目标等算法\[17-19\]较为复杂，现实生产中，多目标规划仍未能成为配矿领域内实用性及通用性更强的方法.目前仍以线性构模为主要手段，存在着不足之处.笔者基于磷矿堆场实际生产工艺过程，建立堆场多目标配矿模型，较完整描述该工艺阶段需要达到的配矿目标基础上，应用Matlab优化工具箱中多目标规划工具，实现了堆场配矿多目标规划模型简便、有效求解，实际操作性强.1多目标优化配矿模型1.1多目标优化原理多目标优化问题(multiobjective optimization problem, MOP)又称为多标准优化问题.不失一般性，一个具有n个决策变量，m个目标变量的多目标优化问题\[21\]可表述为
miny=F(x)=(f1(x),f2(x),…,fm(x))T
s,t.gi(x)≤0, i=1,2,…,q
hj(x)=0, j=1,2,…,p(1)其中，x=(x1,…,xn) ∈XRn为n维的决策矢量，X为n维的决策空间，y=(yl,…,ym) ∈YRm为m维的目标矢量，Y为m维的目标空间.目标函数F(x)定义了m个由决策空间向目标空间的映射函数；gi(x)≤0(i=1,2,…,q)定义了q个不等式约束；hj(x)=0(j=1,2,…,p)定义了p个等式约束.对于多目标优化问题有如下重要定义\[2223\]：定义1(可行解)：对于某个x∈X，如果x满足(1)中的约束条件gi(x)≤0 (i=1,2,…,q)和hj(x)=0 (j=1,2,…,p)，则称x为可行解.定义2(可行解集合)：由X中的所有的可行解组成的集合称为可行解集合，记为Xf，且XfX.定义3(Pareto占优)：假设xA,xB∈Xf是式(1)所示多目标优化问题的两个可行解，则称与xB相比xA是Pareto占优的，当且仅当i=1,2,…,m,fi(xA)≤fi(xB)Λj=1,2,…,m,fj(xA)<fj(xB)(2)记作xAxB，也称为xA支配xB.定义4(Pareto最优解)：一个解x*∈Xf被称为Pareto最优解(或非支配解)，当且仅当满足如下条件:x∈Xf∶xx*(3)定义5(Pareto最优解集)：Pareto最优解集是所有Pareto最优解的集合，定义如下:P*△={x*|x∈Xf∶xx*}(4)多目标优化问题各子目标之间具有较强的耦合性及冲突性，多个子目标不可能同时达到最优，只能在他们中间进行协调和折中处理，使各子目标函数尽可能达到最优.Pareto最优解即是一个可以接受的“不坏”的解，通常一个多目标问题会具有多个Pareto 最优解.在实际应用问题中，必须根据对问题的了解程度和决策人员的个人偏好，从 Pareto 最优解集中挑选一个或一些解作为多目标优化问题的最优解.第10期李志国，等：磷矿堆场多目标优化配矿模型
武汉工程大学学报第34卷
1.2磷矿堆场多目标优化配矿模型随着磷矿逐年开采，富矿越来越少，磷矿石品质逐渐变差，越来越多的矿石采出后不易直接进入浮选流程.通常运至堆场分堆堆放，根据入选原矿要求进行矿石混配后再进入破碎—磨浮流程.入选原矿一般要求满足各组分指标的同时，最大限度利用矿石.假设有n堆原矿参与配矿，各堆矿量为Mi吨（i=1,…,n），主要品质指标为P2O5、MgO及倍半氧化物，目前考察矿山倍半氧化物均满足矿石入选要求，该配矿模型暂不考虑倍半指标要求，故考虑品质参数如下：pi为第i堆原矿石P2O5的含量，%，（i=1,…,n）；mgi为第i堆原矿石MgO的含量，%，（i=1,…,n）；参照多目标规划一般形式（1）可对配矿模型表示入下：
max∑ni=1mi
maxpimi∑ni=1mi
minmgimi∑ni=1mi
s.t.mi≤Mi
mi≥0i=1,…,n(5)式(5)中：mi为决策变量，即第i堆矿石配矿量，吨.其中，目标函数具有理想目标值.通常提高矿石P2O5的含量，MgO等杂质的含量将会升高，两目标相互冲突，即某一子目标的改善将引起其他子目标性能的降低.故模型为典型多目标规划问题.2多目标优化模型解算处理多目标优化的方法比较多，如约束法、评价函数法、功效系数法等\[24\].随着计算机应用技术的发展，各种用于多目标优化解算的软件层出不穷.Matlab作为三大数学软件之一，其强大的科学计算功能和覆盖面广、专业性强的工具箱发展成为适合多学科、跨平台的大型实用科学计算软件，也为最优化问题的解决提供了有力工具.2.1基于Matlab多目标优化模型解算技术Matlab优化工具箱功能涵盖：无约束条件非线性极小值；约束条件下非线性极小值，包括目标逼近问题、极大-极小值问题和半无限极小值问题；二次规划和线性规划问题；非线性最小二乘逼近和曲线拟合；非线性系统的方程求解；约束条件下的线性最小二乘优化；复杂结构的大规模优化问题等.其中对于多目标规划问题，Matlab可以利用不同函数进行求解.例如：在评价函数法中我们所得最后评价函数为一线性函数，且约束条件也为线性函数，则可直接利用工具箱提供的线性规划函数linprog进行求解；若得到评价函数为非线性函数，则可利用fmincon函数求解；如果采用极大极小法求解，则可直接利用工具箱中fminimax函数求解.对于式(5)表示的配矿模型，原矿石P2O5及MgO的含量均有一确定目标值，可直接调用fgoalattain函数求解.其一般调用格式如下：\[x,fval\]=fgoalattain(fun,x0,goal,weight,A,b,Aeq,beq,lb,ub,nonlcon,options)其中，x、fval为函数返回的决策变量及目标函数向量值；fun为需要最小化的目标函数；x0为求解的初始值；goal为目标函数期望值；weight为目标权重，其他相关信息详见Matlab优化工具箱手册.选厂对入选原矿不同组分的品质要求，可以通过对goal与weight参数的控制设置来实现.输入参数中的weight变量为权重向量，可控制低于或超过函数指定目标的相对程度；设置硬约束，以表示某一组分含量严格达到目标值.2.2解算实例为了验证该算法的有效性，计算文献\[25\]实例如下:目标规划模型为：minf1(x)=3x2+2x4
maxf2(x)=10x1+9x2+8x3+7x4
s.t.x1+x2≥30
x3+x4≥30
x1+2x3≤120
x1,x2,x3,x4≥0（6）f1(x)和f2(x)在同一优先级下，f1(x)的上限值为20，f2(x)的下限值为800，采用本文解法及文献\[25\]解法结果比较如表1.表1多目标规划不同解算方法结果对比
Table 1Results comparison during different
solving method for multiobjective optimization
参数x1x2x3x4f1f2单纯型解030607.143104.286800文献\[25\]
算法解69.909
64.802
58.5390.794
0.790
0.20824.681
27.056
29.8848.688
3.579
1.02319.758
9.528
2.670964.501
896.635
833.493本文算法58.09303000820.927从求解结果对比可以发现，单纯型法得到的解不是决策者满意解，而采用Matlab优化工具箱fgoalattain函数得到的解及文献\[25\]遗传算法得到的三个解均为决策者满意解.并且，针对该类具有目标值约束或确定期望值的特殊小型多目标规划问题，fgoalattain函数具有参数设置方便，求解简单的优势.3磷矿堆场多目标优化配矿实现针对某磷矿堆场多目标优化配矿模型，应用Matlab优化工具箱中多目标优化函数求解实验如下：某矿浮选厂堆场一月份共堆存来自4处不同矿源的原矿待处理，各堆参数如表2所示：表2某矿堆场入选原矿库存情况表
Table 2Inventory of raw ore in a stockyard
矿堆参数ABCD矿量/t18 174.19 908.846 896.6825 634.29w(P2O5)/%23.1622.6921.6823.47w(MgO)/%4.094.263.664.7选厂入选原矿石主要指标要求如表3所示：表3某浮选厂入选原矿品质要求
Table 3Ore quality requirements for mineral processing
Cw(P2O5)/%w(MgO)/%指标参数23.054.2假设A、B、C、D矿堆配矿量分别为：m1、m2、m3、m4吨，根据式（5）构建优化模型如下：
max23.16m1+22.69m2+21.68m3+23.47m4m1+m2+m3+m4
min4.09m1+4.26m2+3.66m3+4.7m4m1+m2+m3+m4
maxm1+m2+m3+m4
s.t.0≤m1≤18 174.1
0≤m2≤9 908.84
0≤m3≤6 896.68
0≤m4≤25 634.29(7)设置目标函数期望值为60 613.91，23.05,4.2；严格控制P2O5含量不低于23.05%，MgO不超过4.2%.应用fgoalattain函数解算结果如下：m1=18 174.1，m2=0，m3=4 014.0，m4=8 333.4即A∶C∶D≈9∶2∶4对应目标函数结果为：30 521.5，23.05，4.2；计算结果显示按照A∶C∶D=9∶2∶4混配后，可以得到P2O5含量23.05%，MgO含量4.2%的矿石30 521.5 t.实际配矿操作后，混配矿石30 290 t，实测P2O5含量23.05%，MgO含量4.19%.矿石主要指标均满足入选矿石要求，充分利用了低磷（C）及高镁（D）磷矿石，配矿效果明显.4结语a. 通过分析磷矿堆场配矿过程，选取P2O5、MgO质量指标及磷矿石最大限度利用为目标，以堆场资源量及决策变量非负为约束，构建了磷矿堆场多目标配矿模型.b. 磷矿堆场多目标配矿模型简单、直观、较完整的描述了堆场配矿工艺过程及目的，很好体现了矿石充分利用的宗旨，符合资源节约型社会建设要求.c. 应用Matlab优化工具箱多目标规划相关函数，实现了磷矿堆场多目标配矿模型简便、有效求解，证明了Matlab优化工具箱对类似堆场多目标配矿模型求解的有效性.