CN106650997A - 一种选矿生产指标的优化方法 - Google Patents

一种选矿生产指标的优化方法 Download PDF

Info

Publication number
CN106650997A
CN106650997A CN201610948372.6A CN201610948372A CN106650997A CN 106650997 A CN106650997 A CN 106650997A CN 201610948372 A CN201610948372 A CN 201610948372A CN 106650997 A CN106650997 A CN 106650997A
Authority
CN
China
Prior art keywords
metal
amount
grade
mineral processing
month
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201610948372.6A
Other languages
English (en)
Inventor
余刚
周俊武
邹国斌
宋涛
杨佳伟
王庆凯
边维天
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing General Research Institute of Mining and Metallurgy
Original Assignee
Beijing General Research Institute of Mining and Metallurgy
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing General Research Institute of Mining and Metallurgy filed Critical Beijing General Research Institute of Mining and Metallurgy
Priority to CN201610948372.6A priority Critical patent/CN106650997A/zh
Publication of CN106650997A publication Critical patent/CN106650997A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Mining

Abstract

本发明公开了一种选矿生产指标的优化方法,首先设定选矿生产指标优化方法相关的参数、决策变量和辅助变量;建立选矿生产指标优化的目标模型;建立选矿生产的约束条件,所述约束条件包括综合指标目标范围约束,以及生产能力、工艺水平、生产材料和能源约束;采用进化算法对带约束的目标模型进行求解,得到选矿生产作业的原矿处理量、金精矿品位和金属量优化设定值,实现选矿生产指标的优化。利用该优化方法可以在选矿厂现有生产流程、设备、工艺水平及其他资源条件下,实现选矿厂的综合效益最大化,提高选矿厂的生产技术水平。

Description

一种选矿生产指标的优化方法
技术领域
本发明涉及选矿控制与优化技术领域,尤其涉及一种选矿生产指标的优化方法。
背景技术
目前,针对黄金的选矿生产过程是一个复杂的多阶段连续生产过程,而选矿生产指标的决策主要是在选矿厂现有生产流程、设备、工艺水平及其他资源条件下,确定合理的月、日生产指标包括作业量、精矿品位等,并根据实际生产指标情况进行调整,以实现矿山企业的经营目标,并力求提高生产经营效益。选矿生产指标的决策与调度调整是一个复杂动态过程,具有下列难度,一是跨越时间尺度长,需要在年月日时间尺度内,考虑合理的产能分配;二是存在多种约束,包括生产能力、工艺水平、生产材料、药剂、能源等多类约束;三是存在不确定干扰,包括原料性质波动、设备随机故障等,对稳定和提升生产指标带来困难。
现有技术方案一般采用人工决策法,该方法主要依靠经验,不仅耗费人力,而且不能及时应对变化,调整手段比较粗放,生产指标波动大,影响经济效益提升。
发明内容
本发明的目的是提供一种选矿生产指标的优化方法,利用该优化方法可以在选矿厂现有生产流程、设备、工艺水平及其他资源条件下,实现选矿厂的综合效益最大化,提高选矿厂的生产技术水平。
一种选矿生产指标的优化方法,所述方法包括:
设定选矿生产指标优化方法相关的参数、决策变量和辅助变量;
建立选矿生产指标优化的目标模型;
建立选矿生产的约束条件,所述约束条件包括综合指标目标范围约束,以及生产能力、工艺水平、生产材料和能源约束;
采用进化算法对带约束的目标模型进行求解,得到选矿生产作业的原矿处理量、金精矿品位和金属量优化设定值,实现选矿生产指标的优化。
由上述本发明提供的技术方案可以看出,利用该优化方法可以在选矿厂现有生产流程、设备、工艺水平及其他资源条件下,实现选矿厂的综合效益最大化,提高选矿厂的生产技术水平。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他附图。
图1为本发明实施例所提供选矿生产指标的优化方法流程示意图。
具体实施方式
下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明的保护范围。
下面将结合附图对本发明实施例作进一步地详细描述,如图1所示为本发明实施例所提供选矿生产指标的优化方法流程示意图,所述方法包括:
步骤11:设定选矿生产指标优化方法相关的参数、决策变量和辅助变量;
在该步骤中,首先设定相关变量,具体所设定的选矿生产指标优化方法相关的参数包括:
αm第m月的平均原矿品位;
w1原矿处理量对经济效益贡献的权重系数;
w2金属量对经济效益贡献的权重系数;
C1,1,basic核定工资单价,按单位原矿处理量提取的费用;
C2,1,basic核定工资单价,按单位金属量提取的费用;
C1,3原矿处理量的计划调整系数;
C2,3金属量的计划调整系数;
C4全月尾矿品位低于或超过基准值时的奖惩费用,元/(0.001g/t);
C5月生产过程中已经违反班次尾矿品位基准的惩罚费用,元;
C6月生产过程中已经违反车次精矿品位基准的惩罚费用,元;
C7月份金属平衡误差率超限惩罚费用,元;
C8峰平谷用电超限惩罚费用,元;
Cother其他惩罚费用,元;
QRW,basic月原矿处理量目标值,吨;
QRW,m,L,QRW,m,H第m月原矿处理量目标下限和上限;
QRW,davg,H日原矿处理量能力上限;
QRW,consumed当月累计已完成的原矿处理量,吨;
MCC,dct_basic月金属量目标值(返金率折算后),kg;
MCC,m,L,MCC,m,H第m月金属量(返金率折算后)目标下限和上限,kg;
MCC,dct_produced当月累计已生产的金属量(返金率折算量),kg;
MTL,produced当月累计已生产尾矿中所含的金属量(未折算),kg;
LH]精矿品位目标区间,g/t;
vbasic尾矿品位基准值g/t;
vL尾矿品位工艺水平下限g/t;
vH尾矿品位约束上限g/t;
εL金属回收率下限%;
qwater水量单耗t/t;
Qwater_day,H日供水能力上限t;
Tremain第m月中剩余的生产运行时间,单位天(如20.5天);
L1H1],[βL2H2],...[βLnHn]各精矿品位等级范围;
D1,D2,...Dn不同等级精矿品位的金属量折算率;
P1,P2生产完成率。
所设定的决策变量具体包括:
X=(X1,davg,X2,davg,X3,davg)T决策变量矢量;其中,X1,davg为日均原矿处理量,X2,davg为日均金属量,X3,davg为日均精矿品位。
所设定得辅助变量具体包括:
QRW,m(X)第m月的原矿总处理量;
QCC,m(X)第m月内的精矿产量;
MCC,m(X)第m月精矿产品中的金属量;
QTL,m(X)第m月内的尾矿产量;
MTL,m(X)第m月尾矿产品中的金属量;
QTL,davg(X)第m月的日均尾矿产量;
MTL,davg(X)第m月的日均尾矿金属量;
fdiscount(X3,davg)销售返金率,即销售金精矿中的金属量按品位等级折算的比例;
ftail_grade,davg(X)日均尾矿品位;
ftail_grade,m(X)第m月总尾矿品位;
C1,1(X)按原矿处理量的计划完成率核定的单位效益;
C1,2(X)原矿处理量的计划完成率;
C2,1(X)按金属量的计划完成率核定的单位效益;
C2,2(X)金属量的计划完成率。
步骤12:建立选矿生产指标优化的目标模型;
在该步骤中,所建立的目标模型具体为:
J1,m=A1(X)+A2(X)-C(X);
A1(X)表示按金矿原矿处理量完成率折算的相对经济效益;
A2(X)表示按黄金金属量完成率折算的相对经济效益;
C(X)表示品位超限、精矿品位超限、实际金属平衡误差率超限、实际峰平谷用电超限后的效益惩罚。
具体实现中,需要针对该目标模型实现对选矿厂的相对经济效益最大化maxJ1,m
举例来说,依据上述步骤11所创建的参数、决策变量和辅助变量来建立目标模型,一种典型的方法为:
选矿厂的相对经济效益最大化maxJ1,m (1)
其中:
QRW,m(X)=(X1,davgTremain+QRW,consumed) (3)
QTL,m(X)=QRW,m(X)-QCC,m(X) (5)
MCC,m(X)=(X2,davgTremainfdiscount(X3,davg)+MCC,dct_produced) (6)
MTL,m(X)=(0.001×αmX1,davg-X2,davg)Tremain+MTL,produced (7)
上述所建立的公式中,目标函数为最大化相对经济效益maxJ1,m,如式(1)所示;相对经济效益性能指标函数如式(2)所示。
指标函数包括三部分,其中w1QRW,m(X)×C1,1(X)×C1,2(X)×C1,3,表示按金矿原矿处理量计划完成率折算的相对经济效益部分,w2MCC,m(X)×C2,1(X)×C2,2(X)×C2,3表示按黄金金属量计划完成率折算的相对经济效益部分,1000×C4×(ftail_grade,m(X)-vbasic)表示全月尾矿品位低于或超过基准值时的奖惩费用,尾矿品位高于基准值时,金属流失浪费进行惩罚,反之进行奖励。
-C5-C6-C7-C8-C9表示对计划期内实际发生的班次尾矿品位超限、精矿品位超限、实际金属平衡误差率超限、实际峰平谷用电超限进行惩罚。
式(3)、(4)、(5)分别为金矿的原矿、精矿、尾矿处理量与决策变量的关系表达式。
式(6)、(7)分别为精矿、尾矿中所含金属量与决策变量关系表达式。
式(8)表示计划期总尾矿品位的关系表示式。
式(9)为按原矿处理量的完成率核定的效益系数计算方法,按基准值和完成率进行折算。
式(10)为原矿处理量的完成率与决策变量的关系表达式。
式(11)为销售返金率的表达式,即销售金精矿中的金属量按品位等级折算的金属量比例。
式(12)为按金属量的完成率核定的效益系数计算方法,按基准值和完成率进行折算。
式(13)为金属量的完成率与决策变量的关系表达式。
步骤13:建立选矿生产的约束条件;
在该步骤中,所述约束条件包括综合指标目标范围约束,以及生产能力、工艺水平、生产材料和能源约束,其中具体包括有:
1)综合指标目标范围约束表示为:
QRW,m,L≤QRW,m(X)≤QRW,m,H
MCC,m,L≤MCC,m(X)≤MCC,m,H
2)生产能力约束表示为:
X1,davg≤QRW,davg,H
3)物料与金属平衡关系约束
0.001×αmX1,davg=X2,davg+MTL,davg(X)
4)精矿品位目标约束:
βL≤X3,davg≤βH
5)尾矿品位目标范围及工艺水平约束:
ftail_grade,m(X)≤vH
vL≤ftail_grade,davg(X)≤vH
6)金属回收率约束:
7)供水能力约束:
X1,davgqwater≤Qwater_day,H
步骤14:采用进化算法对带约束的目标模型进行求解,得到选矿生产作业的原矿处理量、金精矿品位和金属量优化设定值,实现选矿生产指标的优化。
在该步骤中,一种典型的计算过程为:
1)初始化及产生种群:主要是设置初始化参数(如种群大小、进化代数、初始步长等),并产生初始化种群;
2)适应度分配:计算种群中个体的适应值,确定约束支配关系;
3)交配选择:采用二元竞争选择来确定交配个体;
4)个体更新:针对交配池中的个体,以一定的概率应用交叉(重组)、变异算子产生新的个体;
5)环境选择:根据适应度值和约束支配关系选择子代个体,并更新代数gen=gen+1;
6)终止准则:如果达到终止条件(最大代数或最佳个体的适应度值变化率阈值),则结束,否则根据新形成的种群回到第2)步开始下一代搜索。
综上所述,本发明实施例所述方法具有以下优点:
1、通过设定生产指标优的参数、决策变量、辅助变量,建立相对经济效益目标函数与约束函数,采用智能算法快速求解,在选矿厂综合生产目标范围和资源约束条件下,有利于实现选矿厂的综合经济效益最大化;
2、考虑了选矿厂的实际指标完成情况,以及矿石属性、设备能力等资源条件随生产周期变化情况,有利于生产指标的稳定和快速调整,为生产过程提供了控制目标和方向;
3、在同样满足生产约束条件下,本发明实施例所提供的生产指标优化方法比人工经验设定方法得到的相对经济效益超过20%以上。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书的保护范围为准。

Claims (5)

1.一种选矿生产指标的优化方法,其特征在于,所述方法包括:
设定选矿生产指标优化方法相关的参数、决策变量和辅助变量;
建立选矿生产指标优化的目标模型;
建立选矿生产的约束条件,所述约束条件包括综合指标目标范围约束,以及生产能力、工艺水平、生产材料和能源约束;
采用进化算法对带约束的目标模型进行求解,得到选矿生产作业的原矿处理量、金精矿品位和金属量优化设定值,实现选矿生产指标的优化。
2.根据权利要求1所述选矿生产指标的优化方法,其特征在于,所设定的选矿生产指标优化方法相关的参数具体包括:
αm第m月的平均原矿品位;
w1原矿处理量对经济效益贡献的权重系数;
w2金属量对经济效益贡献的权重系数;
C1,1,basic核定工资单价,按单位原矿处理量提取的费用;
C2,1,basic核定工资单价,按单位金属量提取的费用;
C1,3原矿处理量的计划调整系数;
C2,3金属量的计划调整系数;
C4全月尾矿品位低于或超过基准值时的奖惩费用;
C5月生产过程中违反班次尾矿品位基准的惩罚费用;
C6月生产过程中违反车次精矿品位基准的惩罚费用;
C7月份金属平衡误差率超限惩罚费用;
C8峰平谷用电超限惩罚费用;
Cother其他惩罚费用;
QRW,basic月原矿处理量目标值;
QRW,m,L,QRW,m,H第m月原矿处理量目标下限和上限;
QRW,davg,H日原矿处理量能力上限;
QRW,consumed当月累计已完成的原矿处理量;
MCC,dct_basic月金属量目标值;
MCC,m,L,MCC,m,H第m月金属量目标下限和上限;
MCC,dct_produced当月累计已生产的金属量;
MTL,produced当月累计已生产尾矿中所含的金属量;
LH]精矿品位目标区间;
vbasic尾矿品位基准值;
vL尾矿品位工艺水平下限;
vH尾矿品位约束上限;
εL金属回收率下限;
qwater水量单耗;
Qwater_day,H日供水能力上限;
Tremain第m月中剩余的生产运行时间;
L1H1],[βL2H2],...[βLnHn]各精矿品位等级范围;
D1,D2,...Dn不同等级精矿品位的金属量折算率;
P1,P2生产完成率。
3.根据权利要求1所述选矿生产指标的优化方法,其特征在于,所设定的决策变量具体包括:
X=(X1,davg,X2,davg,X3,davg)T决策变量矢量;其中,X1,davg为日均原矿处理量,X2,davg为日均金属量,X3,davg为日均精矿品位。
4.根据权利要求1所述选矿生产指标的优化方法,其特征在于,所设定得辅助变量具体包括:
QRW,m(X)第m月的原矿总处理量;
QCC,m(X)第m月内的精矿产量;
MCC,m(X)第m月精矿产品中的金属量;
QTL,m(X)第m月内的尾矿产量;
MTL,m(X)第m月尾矿产品中的金属量;
QTL,davg(X)第m月的日均尾矿产量;
MTL,davg(X)第m月的日均尾矿金属量;
fdiscount(X3,davg)销售返金率;
ftail_grade,davg(X)日均尾矿品位;
ftail_grade,m(X)第m月总尾矿品位;
C1,1(X)按原矿处理量的计划完成率核定的单位效益;
C1,2(X)原矿处理量的计划完成率;
C2,1(X)按金属量的计划完成率核定的单位效益;
C2,2(X)金属量的计划完成率。
5.根据权利要求1所述选矿生产指标的优化方法,其特征在于,所建立的目标模型具体为:
J1,m=A1(X)+A2(X)-C(X);
A1(X)表示按金矿原矿处理量完成率折算的相对经济效益;
A2(X)表示按黄金金属量完成率折算的相对经济效益;
C(X)表示品位超限、精矿品位超限、实际金属平衡误差率超限、实际峰平谷用电超限后的效益惩罚。
CN201610948372.6A 2016-10-26 2016-10-26 一种选矿生产指标的优化方法 Pending CN106650997A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610948372.6A CN106650997A (zh) 2016-10-26 2016-10-26 一种选矿生产指标的优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610948372.6A CN106650997A (zh) 2016-10-26 2016-10-26 一种选矿生产指标的优化方法

Publications (1)

Publication Number Publication Date
CN106650997A true CN106650997A (zh) 2017-05-10

Family

ID=58821470

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610948372.6A Pending CN106650997A (zh) 2016-10-26 2016-10-26 一种选矿生产指标的优化方法

Country Status (1)

Country Link
CN (1) CN106650997A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108090695A (zh) * 2018-01-04 2018-05-29 中国有色金属长沙勘察设计研究院有限公司 一种堆积型铝土矿多矿区的开采规划方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009100932A2 (de) * 2008-02-15 2009-08-20 Brandenburgische Technische Universität Cottbus Vorrichtung und verfahren zur aufladung von aerosolen
CN103617456A (zh) * 2013-12-04 2014-03-05 东北大学 一种选矿过程运行指标优化方法
CN103617470A (zh) * 2013-12-19 2014-03-05 东北大学 一种设备能力变化条件下的选矿综合生产指标优化方法
CN104267695A (zh) * 2014-09-23 2015-01-07 东北大学 一种选矿过程工艺指标优化决策系统及方法
CN104570739A (zh) * 2015-01-07 2015-04-29 东北大学 基于云和移动终端的选矿多生产指标优化决策系统及方法
CN105447585A (zh) * 2014-08-23 2016-03-30 沈阳东大自动化有限公司 一种选矿生产全流程综合生产指标优化决策系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009100932A2 (de) * 2008-02-15 2009-08-20 Brandenburgische Technische Universität Cottbus Vorrichtung und verfahren zur aufladung von aerosolen
CN103617456A (zh) * 2013-12-04 2014-03-05 东北大学 一种选矿过程运行指标优化方法
CN103617470A (zh) * 2013-12-19 2014-03-05 东北大学 一种设备能力变化条件下的选矿综合生产指标优化方法
CN105447585A (zh) * 2014-08-23 2016-03-30 沈阳东大自动化有限公司 一种选矿生产全流程综合生产指标优化决策系统
CN104267695A (zh) * 2014-09-23 2015-01-07 东北大学 一种选矿过程工艺指标优化决策系统及方法
CN104570739A (zh) * 2015-01-07 2015-04-29 东北大学 基于云和移动终端的选矿多生产指标优化决策系统及方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108090695A (zh) * 2018-01-04 2018-05-29 中国有色金属长沙勘察设计研究院有限公司 一种堆积型铝土矿多矿区的开采规划方法
CN108090695B (zh) * 2018-01-04 2021-09-14 中国有色金属长沙勘察设计研究院有限公司 一种堆积型铝土矿多矿区的开采规划方法

Similar Documents

Publication Publication Date Title
CN110378799B (zh) 基于多尺度深度卷积网络的氧化铝综合生产指标决策方法
CN102266927B (zh) 一种连铸机钢水浇铸顺序的控制方法
CN103617470A (zh) 一种设备能力变化条件下的选矿综合生产指标优化方法
McCarthy et al. Macro food policy planning: a general equilibrium model for Pakistan
CN104835074A (zh) 一种露天矿生产配矿方法
CN107704965B (zh) 一种兼顾生态保护需求的水库兴利调度图优化方法
CN102999019A (zh) 海水淡化工厂系统
CN104123587A (zh) Mto-mts管理模式下钢铁生产合同计划和多级库存匹配优化方法
CN106650999B (zh) 一种啤酒生产调度优化方法
CN102621941B (zh) 一种冷轧电镀锌机组钢卷排序方法及系统
CN105631759A (zh) 炼钢厂考虑铁水供应条件的多目标调度计划编制方法
CN106650997A (zh) 一种选矿生产指标的优化方法
CN103488085A (zh) 甲醇四塔精馏系统多目标优化控制方法
CN104299045B (zh) 一种选矿生产全流程精矿产量预报系统及方法
CN104281134A (zh) 基于人机交互的原矿选别过程多生产指标优化系统及方法
CN104318313B (zh) 基于案例推理的选矿全流程运行指标智能决策系统及方法
CN103246935B (zh) 梯级电站的发电后评估系统及方法
CN116342159A (zh) 一种多金属矿石可采经济品位模型的构建方法
CN114742422A (zh) 一种考虑低碳排放约束的原油调度优化方法及装置
Furen The Relationship Between Accumulation and Consumption in China's Economic Development
CN101604420A (zh) 煤气平衡模型的建立方法
CN101333669A (zh) 锌电解过程电能优化方法
CN107657352A (zh) 一种聚氯乙烯全流程生产规划的求解方法和装置
CN110619395A (zh) 一种回收产品重用的产品配置方法
Sen Investment Allocation in Development Planning

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20170510

RJ01 Rejection of invention patent application after publication