发明名称 用于选矿生产过程的配矿方法
摘要 一种用于选矿生产过程的配矿方法,包括采集相关参数;确定每时段筛选期望的铁精矿品位;建立模型,配置参数;评价每个微粒的适应度;确定全局最好的适应度;输出原矿配矿比例等六个步骤,本发明为选矿生产提供了一种有效的配矿方法,改变了由人工计算造成的误差和凭人工经验进行金属平衡分析造成的主观失误,解决了选矿厂生产过程的科学配矿问题。
申请公布号 CN1680945B 申请公布日期 2011.05.11
申请号 CN200410021249.7 申请日期 2004.04.09
申请人 东北大学 发明人 柴天佑;李慧莹;黄肖玲
分类号 G06F11/00(2006.01)I;C22B1/00(2006.01)I 主分类号 G06F11/00(2006.01)I
代理机构 沈阳东大专利代理有限公司 21109 代理人 梁焱
主权项 1.一种用于选矿生产过程的配矿方法,其特征在于包括以下六个步骤:1)采集相关参数;采集的参数包括一周原矿种类及库存、设备生产能力及精矿需求量;2)确定每时段筛选期望的铁精矿品位;制定选矿过程周计划时决策者期望每时段的筛选的铁精矿的金属品位尽可能接近P<sub>0</sub>(t);3)建立模型,配置参数;所述建立的模型为:<img file="DEST_PATH_FSB00000304112600011.GIF" wi="1590" he="181" />s.t.<img file="DEST_PATH_FSB00000304112600012.GIF" wi="1643" he="151" /><img file="DEST_PATH_FSB00000304112600013.GIF" wi="1601" he="155" /><img file="DEST_PATH_FSB00000304112600014.GIF" wi="1595" he="181" /><img file="DEST_PATH_FSB00000304112600015.GIF" wi="1612" he="323" />x<sub>i</sub>(t)≥0,i=1,...,N,t=1,...,T;                            (6)其中,N为进行磁选的原矿种类,G<sub>i</sub>为第i种矿石的总量,μ<sub>il</sub>为粉矿和块矿的比例,ρ<sub>il</sub>为铁元素含量,η<sub>il</sub>为金属回收率,σ<sub>1</sub>为精矿产率;<img file="DEST_PATH_FSB00000304112600016.GIF" wi="94" he="72" />为在计划期[1,T]范围内t时段用于筛选粉矿和块矿的设备额定生产能力;S<sub>0</sub>为库存中铁精矿的原始储量,M(t)为各时段烧结厂所需的精矿量;X<sub>i</sub>(t)代表第i种原矿在时间t的投入量;P(t)代表实际综合精矿品位,P<sub>0</sub>(t)为金属平衡品位;其约束条件为: (1)是使得各时段金属平衡偏差最小;(2)是指各类原矿总量限制;(3)是指中间产品生产的设备能力限制;(4)保证铁精矿需求量限制;(5)是指综合铁精矿品位;(6)是指第i种原矿在时间t的投入量大于0;所述配置参数包括:微粒数量m、惯性权重、加速常数、最大速度/最大运行代数;4)评价每个微粒的适应度;定义ΔG<sub>i</sub>、ΔCE<sub>1</sub>和ΔM(t)代表各约束的违反量,其中:<img file="DEST_PATH_FSB00000304112600021.GIF" wi="547" he="149" /><img file="DEST_PATH_FSB00000304112600022.GIF" wi="781" he="151" /><img file="DEST_PATH_FSB00000304112600023.GIF" wi="990" he="152" />评价函数选择如下的乘法形式:eval(x)=P<sub>0</sub>-f(x)×P(x),这里:<img file="DEST_PATH_FSB00000304112600024.GIF" wi="1321" he="151" />其中:<img file="DEST_PATH_FSB00000304112600025.GIF" wi="149" he="69" /><img file="DEST_PATH_FSB00000304112600026.GIF" wi="197" he="63" />和ΔM<sup>max</sup>(t)代表当前种群中各约束的最大违反量;<img file="DEST_PATH_FSB00000304112600027.GIF" wi="413" he="69" /><img file="DEST_PATH_FSB00000304112600028.GIF" wi="549" he="64" />ΔM<sup>max</sup>(t)=max{ε,ΔM(t)};对于每个微粒,将其适应度与经历过的最好位置做比较;如果较好,则将其作为当前的最好位置;5)确定全局最好的适应度;对于每个微粒,将其适应度与全局所经历过的最好位置做比较;如果较好,则重新设置索引号;根据公式调整微粒的速度和位置;粒子速度和位置调整公式如下:<img file="DEST_PATH_FSB00000304112600029.GIF" wi="1155" he="111" />x<sub>i</sub>(t)=x<sub>i</sub>(t)+υ<sub>it</sub>;6)输出原矿配矿比例;根据最大代数判断是否达到算法终止条件,如达到则输出原矿配矿比例值,如未满足则 返回步骤4)。 
地址 110004 辽宁省沈阳市和平区文华路3号巷11号