CN106950946B - 一种基于优化原则的湿法冶金异常控制方法 - Google Patents
一种基于优化原则的湿法冶金异常控制方法 Download PDFInfo
- Publication number
- CN106950946B CN106950946B CN201710343173.7A CN201710343173A CN106950946B CN 106950946 B CN106950946 B CN 106950946B CN 201710343173 A CN201710343173 A CN 201710343173A CN 106950946 B CN106950946 B CN 106950946B
- Authority
- CN
- China
- Prior art keywords
- abnormal
- hydrometallurgy
- data
- optimization problem
- working condition
- 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.)
- Active
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 133
- 238000000034 method Methods 0.000 title claims abstract description 113
- 238000009854 hydrometallurgy Methods 0.000 title claims abstract description 80
- 230000002159 abnormal effect Effects 0.000 claims abstract description 147
- 230000008569 process Effects 0.000 claims abstract description 73
- 239000002562 thickening agent Substances 0.000 claims description 34
- 238000012545 processing Methods 0.000 claims description 31
- 239000002002 slurry Substances 0.000 claims description 29
- 239000007788 liquid Substances 0.000 claims description 27
- 238000012544 monitoring process Methods 0.000 claims description 15
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 10
- 230000007246 mechanism Effects 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 6
- 238000004062 sedimentation Methods 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 claims description 4
- 238000009792 diffusion process Methods 0.000 claims description 3
- 238000012847 principal component analysis method Methods 0.000 claims description 3
- 230000002547 anomalous effect Effects 0.000 claims 2
- 238000004519 manufacturing process Methods 0.000 description 18
- 238000001514 detection method Methods 0.000 description 15
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 9
- 239000010931 gold Substances 0.000 description 9
- 229910052737 gold Inorganic materials 0.000 description 9
- 238000003745 diagnosis Methods 0.000 description 7
- 238000003825 pressing Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 6
- 229910052802 copper Inorganic materials 0.000 description 6
- 239000010949 copper Substances 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000002386 leaching Methods 0.000 description 5
- 229910052717 sulfur Inorganic materials 0.000 description 5
- 239000011593 sulfur Substances 0.000 description 5
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 4
- 229910052500 inorganic mineral Inorganic materials 0.000 description 4
- 239000011707 mineral Substances 0.000 description 4
- 230000033228 biological regulation Effects 0.000 description 3
- 230000005856 abnormality Effects 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005265 energy consumption Methods 0.000 description 2
- 238000002156 mixing Methods 0.000 description 2
- 239000000047 product Substances 0.000 description 2
- 238000004611 spectroscopical analysis Methods 0.000 description 2
- 238000005406 washing Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000018044 dehydration Effects 0.000 description 1
- 238000006297 dehydration reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000012065 filter cake Substances 0.000 description 1
- 238000005188 flotation Methods 0.000 description 1
- 238000011010 flushing procedure Methods 0.000 description 1
- 238000000227 grinding Methods 0.000 description 1
- 231100000086 high toxicity Toxicity 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 238000005272 metallurgy Methods 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
- 238000009853 pyrometallurgy Methods 0.000 description 1
- 238000012797 qualification Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000010802 sludge Substances 0.000 description 1
- 239000002910 solid waste Substances 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
- 230000008719 thickening Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
- G05B23/0205—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
- G05B23/0218—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
- G05B23/0243—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults model based detection method, e.g. first-principles knowledge model
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/20—Pc systems
- G05B2219/24—Pc safety
- G05B2219/24065—Real time diagnostics
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Paper (AREA)
Abstract
本发明提供一种基于优化原则的湿法冶金异常控制方法,包括:获取湿法冶金过程的在线数据;判断在线数据中是否有异常数据;若有,则辨识异常数据,并获取异常数据的异常工况;根据所述异常数据、异常工况确定湿法冶金异常控制的优化原则;根据优化原则,将湿法冶金异常控制问题转换为单目标优化问题或多目标优化问题;根据单目标优化问题或多目标优化问题的约束条件,求解所述单目标优化问题或多目标优化问题,获得安全处理策略;将安全处理策略发送安全控制系统,以使安全控制系统执行安全处理策略。上述方法能够识别异常工况,并针对异常工况制定有效的安全处理策略,可以降低故障的发生概率。
Description
技术领域
本发明属于湿法冶金技术,尤其涉及一种基于优化原则的湿法冶金异常控制方法。
背景技术
湿法冶金工艺逐渐成熟,其工业自动化水平急需提高,与传统的火法冶金相比,湿法冶金技术具有很多优势,例如:清洁、高效、低品位复杂金属矿产更易于回收等。我国矿产资源贫矿多,杂质含量高,湿法冶金更加适用于我国矿产资源的现状,其工业化水平的提高,对于提高矿产资源的综合利用率,降低固体废弃物产量,减少环境污染,都有重大的意义。
湿法冶金全流程采用了磨矿、浮选、脱水调浆、氰化浸出、压滤洗涤、置换等湿法冶炼的工艺流程。通过选矿工序,根据调浆后的矿浆中金品位及硫含量的高低,3种矿浆(低硫尾矿、高硫矿和高铜矿)分别进入不同的处理流程,低硫尾矿进入低品位碳吸附氰化浸出提金流程,高硫矿进入低品位氰化浸出置换流程,高铜矿进入高品位氰化浸出置换流程。金具有很高的经济价值,因此为了获得更高的经济效益,有必要在金的生产中保持良好的操作性能及安全生产。
近几年湿法冶金工艺、设备研究进展迅速。但是湿法冶金工艺流程复杂,设备类型多样,整个流程具有多变量、变量之间强耦合等综合复杂性、其特性随生产条件变化而变化、受到原料成分、工况、设备状态等多种不确定性因素干扰等特点,如有操作不当、疏于检测或因不可抗拒的自然因素引起设备故障而导致生产中断等事故,将会给生产带来巨大的浪费和损失。同时,湿法冶金高腐蚀、高温甚至高毒的工艺条件更对生产的安全性提出了更高的要求。此外,湿法冶金工艺流程由多个彼此密切关联的子过程构成,一旦出现一个异常工况没有被及时调整,则会产生更多的连锁异常甚至无法挽回的故障,对设备及操作人员的安全造成严重的威胁。
目前,针对金湿法冶金生产中出现的异常工况,工厂所采用的方法是对与生产过程息息相关的过程变量信息进行人为或者仪表的记录,通过人为查看报表来获取过程运行状态信息,当发生异常工况时,操作人员基于机理知识和工艺调整经验给出相应的异常处理决策。但是现场工作人员的人为调整方法完全依赖于操作者各自的经验,很难保证决策的及时性及准确性。
发明内容
针对现有存在的技术问题,本发明提供一种基于优化原则的湿法冶金异常控制方法,该方法可以确定过程变量的正常状态信息或异常状态信息,在发生异常状态时,可及时制定安全策略,提高生产安全性、生产效率、降低事故发生率。
本发明的基于优化原则的湿法冶金异常控制方法,包括:
获取湿法冶金过程的在线数据;
判断在线数据中是否有异常数据;
若有,则辨识异常数据,并获取异常数据的异常工况;
根据所述异常数据、异常工况确定湿法冶金异常控制的优化原则;
根据优化原则,将湿法冶金异常控制问题转换为单目标优化问题或多目标优化问题;
根据单目标优化问题或多目标优化问题的约束条件,求解所述单目标优化问题或多目标优化问题,获得安全处理策略;
将安全处理策略发送安全控制系统,以使安全控制系统执行安全处理策略。
可选地,所述方法还包括:
在湿法冶金过程的异常工况恢复正常之后,重复获取湿法冶金过程的在线数据,以及判断在线数据中异常数据的步骤;
或者,
在安全处理策略发送安全控制系统的步骤之后,间隔预设时间段,重复获取湿法冶金过程的在线数据,以及判断在线数据中异常数据的步骤。
可选地,判断在线数据中是否有异常数据,包括:
查看在线数据中每一数据的数值是否在正常范围内,若否,则为异常数据;
相应地,获取异常数据的异常工况,包括:
确定异常数据所处的工况环节,将该工况环节作为异常工况。数据间的相关关系发生变化并超过一定范围也作为异常工况。
可选地,求解所述单目标优化问题或多目标优化问题,获得安全处理策略,包括:
采用高斯伪谱法求解所述单目标优化问题或多目标优化问题,获得安全处理策略。
可选地,判断在线数据中是否有异常数据,包括:
采用主成分分析法的监测方式实时监测湿法冶金过程中各个工况的状态。
可选地,若异常工况为湿法冶金中压滤机前缓冲槽的矿浆液位工况异常,则
单目标优化问题的约束条件包括:
其中,qe(t)=Qe(t)/A是单位面积的溢流流量;qu(t)=Qu(t)/A是单位面积的底流流量;qf(t)=Qf(t)/A是单位面积的进料流量;A是浓密机的横截面积;Qe(t)是溢流流量;Qu(t)是底流流量;Qf(t)是进料流量;Cf(t)是进料浓度;D是扩散系数;n是层数;z是浓密机高度;Δz=z/n是浓密机中每层的高度;νs(t)是沉降速度;C(t)是矿浆浓度;Cn(t)是浓密机底流浓度;溢流流量Qe(t)为固定值;
Qr(t)是压滤机的返回水流量;h1(t)是浓密机中的矿浆液位;
其中,f1(t)是第一渣浆泵的工作频率;f2(t)是第二渣浆泵的工作频率;h2(t)是第一缓冲槽的矿浆液位;h3(t)是第二缓冲槽的矿浆液位;r1是第一缓冲槽底部半径;r2是第二缓冲槽底部半径;k是频率和流量在数值上的比例系数;
状态限制条件包括:
h1min≤h1(t)≤h1max (5)
h2min≤h2(t)≤h2max (6)
h3min≤h3(t)<h3max (7)
Cmin≤Cn(t)≤Cmax (8)
其中,h1min和h1max分别是状态变量h1(t)的下限和上限;下标min表示对应状态变量的下限,下标max表示对应状态变量的上限。
控制限制条件包括:
f1min≤f1(t)≤f1max (9)
f2min≤f2(t)≤f2max (10)
其中,f1min和f1max分别是f1(t)的低限和高限,下标min表示对应的低限,下标max表示的对应的高限;q1,q2和q3分别是Qf(t),Qu(t)和Qr(t)变化率的上限。
可选地,若异常工况为湿法冶金中压滤机前缓冲槽的矿浆液位工况异常,则湿法冶金异常控制的优化原则是时间最短;
相应地,根据优化原则,将湿法冶金异常控制问题转换为单目标优化问题,包括:
单目标优化的目标函数为从异常状态恢复至正常状态的时间tf最小,tf代表第二缓冲槽中的矿浆异常液位调整回正常的时间,
min tf s.t.(1)-(13) (14)。
可选地,单目标优化问题或多目标优化问题的约束条件包括:湿法冶金过程中工艺机理模型限制、湿法冶金过程各相关变量的工艺限制;
预先建立湿法冶金过程中每一工况的约束条件,以及每一工况中每一变量的约束条件/限定条件;
每一变量为在线数据对应的变量数据。
第二方面,本发明提供一种基于优化原则的湿法冶金异常控制装置,包括:
异常数据分析单元,用于获取湿法冶金过程的在线数据;判断在线数据中是否有异常数据;若有,则辨识异常数据,并获取异常数据的异常工况;
优化单元,用于根据所述异常数据、异常工况确定湿法冶金异常控制的优化原则;根据优化原则,将湿法冶金异常控制问题转换为单目标优化问题或多目标优化问题;根据单目标优化问题或多目标优化问题的约束条件,求解所述单目标优化问题或多目标优化问题,获得安全处理策略;
发送单元,用于将安全处理策略发送安全控制系统,以使安全控制系统执行安全处理策略。
可选地,所述优化单元,具体用于
根据所述异常数据、异常工况确定湿法冶金异常控制的优化原则;根据优化原则,将湿法冶金异常控制问题转换为单目标优化问题或多目标优化问题;根据单目标优化问题或多目标优化问题的约束条件,采用高斯伪谱法求解所述单目标优化问题或多目标优化问题,获得安全处理策略。
本发明的基于优化原则的湿法冶金异常控制方法,能够识别在线数据的异常数据和异常数据对应的异常工况,进而将异常控制问题转换为单目标优化问题或多目标优化问题,根据预先建立的约束条件,求解上述单目标优化问题或多目标优化问题,获得安全处理策略,以使安全控制系统根据安全处理策略进行处理,由此可有效降低故障发生概率,同时提高生产安全性、生产效率。
附图说明
图1为本发明基于优化原则的异常控制过程示意图;
图2为高铜矿浓密压滤过程工艺流程示意图;
图3为缓冲槽2中矿浆液位h3(t)示意图;
图4为浓密机的液位h1(t)示意图;
图5为缓冲槽1中矿浆液位h2(t)示意图;
图6为浓密机底流浓度Cn(t)示意图;
图7为渣浆泵1的工作频率f1(t)示意图;
图8为渣浆泵2的工作频率f2(t)示意图;
图9为浓密机进料流量Qf(t)示意图;
图10为浓密机底流流量Qu(t)示意图;
图11为压滤机返回水流量Qr(t)示意图。
具体实施方式
为了更好的解释本发明,以便于理解,下面结合附图,通过具体实施方式,对本发明作详细描述。
针对湿法冶金的采集和处理装置可包括:湿法冶金过程监测与故障诊断系统、湿法冶金安全控制系统、上位机、PLC、现场传感变送部分。其中现场传感变送部分包括浓度、压力、流量、液位等检测仪表。在湿法冶金过程现场安装检测仪表,检测仪表将采集的信号通过Prof i bus-DP总线送到PLC,PLC通过以太网定时将采集信号传送到上位机,上位机把接受的数据传到湿法冶金过程监测与故障诊断系统,进行生产过程实时监测与诊断,湿法冶金安全控制系统根据监测结果制定调整方案,自动调整相关变量使异常恢复为正常。
上述装置的各部分功能:
①现场传感变送部分:包括浓度、压力、流量、液位等检测仪表由传感器组成,负责过程数据的采集与传送;
②PLC:负责把采集的信号A/D转换,并通过以太网把信号传送给上位机;PLC控制器采用S imens 400系列的CPU 414-2,具有Prof i bus DP口连接分布式I O;为PLC配备以太网通讯模块,用于上位机访问PLC数据;PLC控制器和以太网通讯模块放置在中央控制室中的PLC柜中;
③上位机:收集本地PLC数据,传送给湿法冶金过程监测与故障诊断系统,上位机选用i 7联想计算机,采用W I NDOW XP操作系统;
④全流程的过程监测及故障诊断系统对湿法冶金的运行状态进行实时监测和诊断;
⑤湿法冶金安全控制系统根据监测结果制定调整方案,将结果传送给上位机,控制相关变量从异常恢复到正常。
现有的高铜矿的浓密压滤工艺流程可包括浓密和压滤洗涤工序;过程检测系统主要由浓度检测、压力检测、流量检测、液位检测等构成。
在湿法冶金过程现场安装检测仪表,检测仪表将采集的信号通过Prof i bus-DP传送到PLC中,PLC定时将采集信号通过以太网传送给上位机,上位机把接收的数据传给湿法冶金过程监测与故障诊断系统,若异常工况发生,将异常信号传送给湿法冶金安全控制系统,该系统将制定安全控制决策即安全处理策略用于排除异常。
如图1所示,图1示出了本发明一实施例的湿法冶金异常控制方法,该方法包括:
101、获取湿法冶金过程的在线数据;
102、判断在线数据中是否有异常数据,若有,执行步骤103,否则,结束或重复上述步骤101的步骤。
举例来说,可以查看在线数据中每一数据的数值是否在正常范围内,若否,则为异常数据;以及查看异常数据所处的工况环节,将该工况环节作为异常工况。此外,数据间的相关关系发生变化并超过一定范围也作为异常工况。
或者,采用主成分分析法的监测方式实时监测湿法冶金过程中各个工况的状态。
也就是说,在工业现场,操作人员经常使用简单的限制报警去判定关注的物理量是否超过了设定的阈值。若超过了设定阈值,则认为异常工况发生。也可使用其它的数据监测方法监测系统的状态,例如主成分分析法。
103、若有,则辨识异常数据,并获取异常数据的异常工况。
本实施例中,异常的定义如下:当系统的性能指标不能满足过程需求时,操作者能够通过自动或者手动的调节方式将其返回到满足性能指标的程度,这种工况叫做异常工况。
通常,异常工况体现在下面几个方面:产品的合格率不达标;经济效益达不到预期计划;系统存在安全隐患;能量消耗超过阈值等等。这些方面并不是完全独立的,它们相互影响,相互转化。最直观的异常体现在某些被关注的物理量超出限定范围。
104、根据所述异常数据、异常工况确定湿法冶金异常控制的优化原则。
105、根据优化原则,将湿法冶金异常控制问题转换为单目标优化问题或多目标优化问题。
举例来说,若异常工况为湿法冶金中压滤机前缓冲槽的矿浆液位工况异常,则湿法冶金异常控制的优化原则是时间最短;
相应地,步骤105则具体为,将湿法冶金异常控制问题转换为单目标优化问题,即时间最短目标的优化问题,单目标优化的目标函数为从异常状态恢复至正常状态的时间最小。
在具体应用中,异常控制的优化原则与时间是息息相关的。通常,将异常控制问题转化为时间优化控制问题。即优化的目标函数为从异常工况调节回正常工况的时间最短。
然而,在实践中,还需考虑质量、经济、安全和能量消耗等其它目标,此时,异常控制问题可转换为多目标优化问题。
应该指出的是,在多目标优化问题中,时间指标可给予最高的优先级。
106、根据单目标优化问题或多目标优化问题的约束条件,求解所述单目标优化问题或多目标优化问题,获得安全处理策略。
在本实施例中,约束条件可为预先根据工艺过程或其他模型建立的限制条件。
通常,优化问题的限定条件可包括工艺机理模型限制和过程各相关变量的工艺限制。
另外,在具体应用中,可采用高斯伪谱法求解所述单目标优化问题或多目标优化问题,获得安全处理策略。
目前,高斯伪谱法是一种优越的动态优化问题求解方法,它将最优控制问题转化为非线性规划问题,然后采用适当的约束非线性优化算法可求得数值最优解。在近似精度、收敛速度和计算效率等方面,高斯伪谱法具有更高优势。
107、将安全处理策略发送安全控制系统,以使安全控制系统执行安全处理策略。
上述异常控制方法的目的是制定有效的措施使异常工况恢复到正常工况,以保证系统运行在设定的范围内。在收到异常的信号之后,操作人员应该制定安全措施即安全处理策略去尽快排除异常。
在实际应用中,执行上述步骤101至步骤107之后,还可在湿法冶金过程的异常工况恢复正常之后,重复上述步骤101和102过程。
或者,在一种可选的实现方式中,在安全处理策略发送安全控制系统的步骤之后,间隔预设时间段,重复上述步骤101和102过程。
本实施例的方法,能够识别在线数据的异常数据和异常数据对应的异常工况,进而将异常控制问题转换为单目标优化问题或多目标优化问题,根据预先建立的约束条件,求解上述单目标优化问题或多目标优化问题,获得安全处理策略,以使安全控制系统根据安全处理策略进行处理,由此可有效降低故障发生概率,同时提高生产安全性、生产效率。
结合图2详细来说,高铜矿浓密压滤工艺流程如图2所示,此过程包括两个主要设备浓密机和压滤机,它们被安放在不同的车间。矿浆经过浓密机的沉淀作用以后,进入缓冲槽1(即第一缓冲槽)中,通过渣浆泵1(即第一渣浆泵),矿浆泵进入缓冲槽2(即第二缓冲槽),通过渣浆泵2(即第二渣浆泵),矿浆进入压滤机中,压滤过程得到的滤饼进入后续的浸出工序,压滤机滤出的水则进入浓密机,一部分水进入浓密机底部,用于浓密机底部阻塞时冲开底部淤泥;另一部分水直接进入浓密机,参与浓密过程。湿法冶金过程检测系统主要由浓度检测、压力检测、流量检测、液位检测等构成。
本发明以高铜矿浓密压滤工艺中压滤机前缓冲槽矿浆液位异常这一工况为例,基于优化原则进行异常控制,包括以下步骤:
步骤1:采集数据,分析异常工况,确定优化问题的目标函数及限定条件。
因为金高昂的经济价值,若含金的矿浆冒槽,则会造成严重的经济损失,所以考虑压滤机前缓冲槽矿浆液位异常这一工况具有重要意义。
基于机理分析,若压滤机前缓冲槽2矿浆液位异常,渣浆泵1应该停止,渣浆泵2应该设为最大功率,以使缓冲槽2中的矿浆尽快排出。
然而,在有限时间内,缓冲槽1的容积以及压滤机的处理能力都是有限的。若渣浆泵1被关闭,同时阀3的开度不变,则缓冲槽1中的矿浆液位将可能异常。若渣浆泵2设为最大功率的时间过长,则压滤机可能出现异常。
因此,针对压滤机前缓冲槽2矿浆液位异常这一工况,必须满足一定的限定条件,以便在处理这个异常工况的同时不产生其他异常工况,这其中涉及到两个车间的配合。操作人员针对异常工况制定安全策略时,总是在时间上追求尽快将异常工况恢复为正常工况。因此基于优化原则可将异常控制问题转化为单目标优化问题即时间优化问题。调控时间最短为优化目标,机理分析和工艺限制分别构成优化问题的等式约束和不等式约束。
机理分析构成优化问题的等式约束:
其中qe(t)=Qe(t)/A是单位面积的溢流流量;qu(t)=Qu(t)/A是单位面积的底流流量;qf(t)=Qf(t)/A是单位面积的进料流量;A是浓密机的横截面积;Qe(t)是溢流流量;Qu(t)是底流流量;Qf(t)是进料流量;Cf(t)是进料浓度;D是扩散系数;n是层数;z是浓密机高度;Δz=z/n是浓密机中每层的高度;νs(t)是沉降速度;C(t)是矿浆浓度;Cn(t)是浓密机底流浓度。为了简化分析过程,不考虑溢流的影响。在接下来的分析中,溢流流量Qe(t)被设为固定值。在浓密机中的矿浆液位满足下面的等式
其中Qr(t)是压滤机的返回水流量;h1(t)是浓密机中的矿浆液位。
缓冲槽的工作原理类似于水箱,因此,针对缓冲槽1和缓冲槽2的机理可以用下面的等式来表示
其中f1(t)是渣浆泵1的工作频率;f2(t)是渣浆泵2的工作频率;h2(t)是缓冲槽1的矿浆液位;h3(t)是缓冲槽2的矿浆液位;r1是缓冲槽1底部半径;r2是缓冲槽2底部半径;k是频率和流量在数值上的比例系数。渣浆泵工作频率越高,进出缓冲槽的流量越大。
工艺限制构成优化问题的不等式约束:状态限制和控制限制。状态限制用下面的不等式表示
h1min≤h1(t)≤h1max (5)
h2min≤h2(t)≤h2max (6)
h3min≤h3(t)<h3max (7)
Cmin≤Cn(t)≤Cmax (8)
其中h1min和h1max分别是状态变量h1(t)的下限和上限。针对式(6)-(8),类似的解释可以被获得。即,下标min表示对应状态变量的下限,下标max表示对应状态变量的上限,如h2min、h2max分别是h2(t)的下限、上限;h3min、h3max分别是h3(t)的下限、上限;Cmin、Cmax分别是Cn(t)的下限、上限。
控制限制用下面的不等式表示
f1min≤f1(t)≤f1max (9)
f2min≤f2(t)≤f2max (10)
其中f1min和f1max分别是f1(t)的低限和高限。针对式(9)-(13),类似的解释可以被获得,即,下标min,如Qfmin、Qumin、Qrmin表示Qf(t)、Qu(t)、Qr(t)各自对应的低限,下标maxQfmax、Qumax、Qrmax表示Qf(t)、Qu(t)、Qr(t)各自对应的高限。q1,q2和q3分别是Qf(t),Qu(t)和Qr(t)变化率的上限。
目标函数为时间tf最小。tf代表缓冲槽2中的矿浆异常液位调整回正常的时间。因此形成的时间优化问题为
min tfs.t.(1)-(13) (14)
步骤2:辨识异常
压滤机前缓冲槽矿浆液位异常即h3(t)异常。对于异常工况的辨识,简单的限制报警装置用于确定是否缓冲槽2中的液位超出界限。这个界限通过操作经验获得。根据式(1)-(14),则优化问题形成,具体的参数如下
3m≤h1(t0)≤3.6m,3m≤h1(t)≤4m,3m≤h1(tf)≤3.6m
1m≤h2(t0)≤2.1m,1m≤h2(t)≤2.5m,1m≤h2(tf)≤2.1m
2m≤h3(t0)≤4m,2m≤h3(t)≤4.5m,2m≤h3(tf)≤4m
8kg/m3≤Cn(t0)≤11kg/m3,8kg/m3≤Cn(t)≤12.65kg/m3,8kg/m3≤Cn(tf)≤11kg/m3
0≤f1(t)≤49.8Hz
0≤f2(t)≤49.8Hz
其它参数为D=0.542,k=2.34,r1=10m,r2=20m,A=500m2,Qe(t)=250m3/h和Cf(t)=2.3kg/m3。
上述参数是在利用计算机中的mat lab软件包求解的过程中,需要用到的参数。本实施例仅为举例说明。
步骤3:求解时间优化问题,获得异常条件下的调节决策。
若矿浆液位h3(t)超过4米,则视为矿浆液位异常。为验证提出算法的有效性,考虑异常的液位状态h3(t)为4.3米。使用高斯伪谱法求解上面形成的时间优化问题,设置高斯伪谱法的插入点为5。
步骤4:将获得的调节策略实施,若异常工况排除,则结束,否则,重新执行上述步骤1,重新获取决策方案。
状态变量的变化过程如图3-图6所示。图3为缓冲槽2中矿浆液位h3(t)的示意图;图4为浓密机的液位h1(t)的示意图;图5为缓冲槽1中矿浆液位h2(t)的示意图;图6为浓密机底流浓度Cn(t)的示意图。由图可知:缓冲槽2中矿浆液位h3(t)有下降的趋势;浓密机的液位h1(t),缓冲槽1中矿浆液位h2(t)和浓密机底流浓度Cn(t)均上升。异常状态变量h3(t)恢复到正常的同时,其它状态变量均保持在正常范围内。
在图3中,矿浆液位h3超过4m,就是异常液位,设定现在异常液位为4.3m,图3中h3液位从4.3开始变化,经过0.3168h,异常恢复为4m。与此同时,图4-图6表明在控制此异常状态量的同时,其它状态量也是正常的,即没有出现新的异常状态。
图7-图11为在调节异常状态时应该设定的控制变量的值,即异常的控制决策。一共有5个控制变量,所以有5个图。
控制变量的变化过程如图7-图11所示。图7为渣浆泵1的工作频率f1(t)示意图;图8为渣浆泵2的工作频率f2(t)示意图;图9为浓密机进料流量Qf(t)示意图;图10为浓密机底流流量Qu(t)示意图;图11为压滤机返回水流量Qr(t)示意图。由图可知:渣浆泵1的工作频率f1(t)有下降的趋势,渣浆泵2的工作频率f2(t)有上升的趋势,可知这符合机理分析的调节措施。浓密机进料流量Qf(t),浓密机底流流量Qu(t)和压滤机返回水流量Qr(t)相互配合进行调节,使其它状态变量保持在正常范围内。
通过上面的实例,基于优化原则的金湿法冶金异常控制方法能够根据识别的异常工况制定有效的安全措施,措施实施后能够尽快使异常工况恢复为正常。
本实施例能有效根据金湿法冶金生产过程出现的异常工况,基于优化原则将异常控制问题转化为优化问题求解,制定有效的安全控制决策。
上述方法了解过程知识,明确过程变量的正常和异常状态信息,在发生异常状态时,及时制定安全策略,对于提高生产安全性、提高生产效率、降低事故发生率、提高企业竞争力具有重要意义。
另一方面,本发明还提供一种基于优化原则的湿法冶金异常控制装置,包括:
异常数据分析单元,用于获取湿法冶金过程的在线数据;判断在线数据中是否有异常数据;若有,则辨识异常数据,并获取异常数据的异常工况;
优化单元,用于根据所述异常数据、异常工况确定湿法冶金异常控制的优化原则;根据优化原则,将湿法冶金异常控制问题转换为单目标优化问题或多目标优化问题;根据单目标优化问题或多目标优化问题的约束条件,求解所述单目标优化问题或多目标优化问题,获得安全处理策略;
发送单元,用于将安全处理策略发送安全控制系统,以使安全控制系统执行安全处理策略。
可选地,所述优化单元,具体用于
根据所述异常数据、异常工况确定湿法冶金异常控制的优化原则;根据优化原则,将湿法冶金异常控制问题转换为单目标优化问题或多目标优化问题;根据单目标优化问题或多目标优化问题的约束条件,采用高斯伪谱法求解所述单目标优化问题或多目标优化问题,获得安全处理策略。
以上结合具体实施例描述了本发明的技术原理,这些描述只是为了解释本发明的原理,不能以任何方式解释为对本发明保护范围的限制。基于此处解释,本领域的技术人员不需要付出创造性的劳动即可联想到本发明的其它具体实施方式,这些方式都将落入本实用新型的保护范围之内。
Claims (7)
1.一种基于优化原则的湿法冶金异常控制方法,其特征在于,包括:
获取湿法冶金过程的在线数据;
判断在线数据中是否有异常数据;
若有,则辨识异常数据,并获取异常数据的异常工况;
根据所述异常数据、异常工况确定湿法冶金异常控制的优化原则;
根据优化原则,将湿法冶金异常控制问题转换为单目标优化问题或多目标优化问题;
根据单目标优化问题或多目标优化问题的约束条件,求解所述单目标优化问题或多目标优化问题,获得安全处理策略;
将安全处理策略发送安全控制系统,以使安全控制系统执行安全处理策略;
求解所述单目标优化问题或多目标优化问题,获得安全处理策略,包括:
采用高斯伪谱法求解所述单目标优化问题或多目标优化问题,获得安全处理策略。
2.根据权利要求1所述的方法,其特征在于,所述方法还包括:
在湿法冶金过程的异常工况恢复正常之后,重复获取湿法冶金过程的在线数据,以及判断在线数据中异常数据的步骤;
或者,
在安全处理策略发送安全控制系统的步骤之后,间隔预设时间段,重复获取湿法冶金过程的在线数据,以及判断在线数据中异常数据的步骤。
3.根据权利要求1所述的方法,其特征在于,判断在线数据中是否有异常数据,包括:
查看在线数据中每一数据的数值是否在正常范围内,若否,则为异常数据;
相应地,获取异常数据的异常工况,包括:
确定异常数据所处的工况环节,将该工况环节作为异常工况,以及数据间的相关关系发生变化并超过一定范围也作为异常工况。
4.根据权利要求1所述的方法,其特征在于,判断在线数据中是否有异常数据,包括:
采用主成分分析法的监测方式实时监测湿法冶金过程中各个工况的状态。
5.根据权利要求1所述的方法,其特征在于,若异常工况为湿法冶金中压滤机前缓冲槽的矿浆液位工况异常,则
单目标优化问题的约束条件包括:
其中,qe(t)=Qe(t)/A是单位面积的溢流流量;qu(t)=Qu(t)/A是单位面积的底流流量;qf(t)=Qf(t)/A是单位面积的进料流量;A是浓密机的横截面积;Qe(t)是溢流流量;Qu(t)是底流流量;Qf(t)是进料流量;Cf(t)是进料浓度;D是扩散系数;n是层数;z是浓密机高度;Δz=z/n是浓密机中每层的高度;νs(t)是沉降速度;C(t)是矿浆浓度;Cn(t)是浓密机底流浓度;溢流流量Qe(t)为固定值;
Qr(t)是压滤机的返回水流量;h1(t)是浓密机中的矿浆液位;
其中,f1(t)是第一渣浆泵的工作频率;f2(t)是第二渣浆泵的工作频率;h2(t)是第一缓冲槽的矿浆液位;h3(t)是第二缓冲槽的矿浆液位;r1是第一缓冲槽底部半径;r2是第二缓冲槽底部半径;k是频率和流量在数值上的比例系数;
状态限制条件包括:
h1min≤h1(t)≤h1max (5)
h2min≤h2(t)≤h2max (6)
h3min≤h3(t)<h3max (7)
Cmin≤Cn(t)≤Cmax (8)
其中,h1min和h1max分别是状态变量h1(t)的下限和上限;下标min表示对应状态变量的下限,下标max表示对应状态变量的上限;
控制限制条件包括:
f1min≤f1(t)≤f1max (9)
f2min≤f2(t)≤f2max (10)
其中,f1min和f1max分别是f1(t)的低限和高限,下标min表示对应的低限,下标max表示的对应的高限;q1,q2和q3分别是Qf(t),Qu(t)和Qr(t)变化率的上限。
6.根据权利要求5所述的方法,其特征在于,若异常工况为湿法冶金中压滤机前缓冲槽的矿浆液位工况异常,则湿法冶金异常控制的优化原则是时间最短;
相应地,根据优化原则,将湿法冶金异常控制问题转换为单目标优化问题,包括:
单目标优化的目标函数为从异常状态恢复至正常状态的时间tf最小,tf代表第二缓冲槽中的矿浆异常液位调整回正常的时间,
min tf s.t.(1)-(13) (14)。
7.根据权利要求1至6任一所述的方法,其特征在于,单目标优化问题或多目标优化问题的约束条件包括:湿法冶金过程中工艺机理模型限制、湿法冶金过程各相关变量的工艺限制;
预先建立湿法冶金过程中每一工况的约束条件,以及每一工况中每一变量的约束条件/限定条件;
每一变量为在线数据对应的变量数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710343173.7A CN106950946B (zh) | 2017-05-16 | 2017-05-16 | 一种基于优化原则的湿法冶金异常控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710343173.7A CN106950946B (zh) | 2017-05-16 | 2017-05-16 | 一种基于优化原则的湿法冶金异常控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106950946A CN106950946A (zh) | 2017-07-14 |
CN106950946B true CN106950946B (zh) | 2019-01-22 |
Family
ID=59478771
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710343173.7A Active CN106950946B (zh) | 2017-05-16 | 2017-05-16 | 一种基于优化原则的湿法冶金异常控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106950946B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109276945B (zh) * | 2018-10-24 | 2020-10-23 | 东北大学 | 一种浓密脱水过程入料异常工况的自愈控制方法 |
CN109885012B (zh) * | 2019-04-12 | 2020-12-25 | 东北大学 | 一种金湿法冶金全流程实时优化补偿方法 |
CN110947237B (zh) * | 2019-11-29 | 2021-08-06 | 北京矿冶科技集团有限公司 | 一种基于浓密机-压滤机生产平衡的协同优化控制方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104199389A (zh) * | 2014-06-12 | 2014-12-10 | 东北大学 | 金湿法冶金置换过程的优化方法 |
CN104298187A (zh) * | 2014-06-12 | 2015-01-21 | 东北大学 | 金湿法冶金全流程三层结构过程监测方法 |
-
2017
- 2017-05-16 CN CN201710343173.7A patent/CN106950946B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104199389A (zh) * | 2014-06-12 | 2014-12-10 | 东北大学 | 金湿法冶金置换过程的优化方法 |
CN104298187A (zh) * | 2014-06-12 | 2015-01-21 | 东北大学 | 金湿法冶金全流程三层结构过程监测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106950946A (zh) | 2017-07-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104062953B (zh) | 金湿法冶金过程运行状态在线评价方法 | |
CN106950946B (zh) | 一种基于优化原则的湿法冶金异常控制方法 | |
CN107065834B (zh) | 湿法冶金过程中浓密机的故障诊断方法 | |
CN104298187B (zh) | 金湿法冶金全流程三层结构过程监测方法 | |
CN106959662B (zh) | 一种电熔镁炉异常工况识别及控制方法 | |
CN116823175B (zh) | 一种石油化工生产信息化系统智能运维方法及系统 | |
CN116991130A (zh) | 一种石化生产智能化自动化控制系统和方法 | |
CN115470831B (zh) | 一种基于频繁项集推理的水电信号异常判断方法 | |
CN116187725B (zh) | 一种用于锻造自动线的锻造设备管理系统 | |
CN115358626A (zh) | 一种选矿生产设备安全检查预警系统 | |
CN106296440A (zh) | 基于ann和es集成的变电站告警信息分析决策系统和方法 | |
CN116938676A (zh) | 一种基于数据源分析的通信风险联合预警方法 | |
CN117071046B (zh) | 一种用于镀锌自动滚镀生产线的加工智能管理系统 | |
CN108037747A (zh) | 金湿法冶金全流程分层分块运行状态在线评价方法 | |
CN110532699B (zh) | 基于模糊dcd的湿法冶金浓密洗涤过程故障诊断方法 | |
CN201737784U (zh) | 油田注入水处理及监测装置 | |
CN106886154B (zh) | 基于区间数的湿法冶金全流程建模方法与优化方法 | |
CN108709426B (zh) | 基于频谱特征双边检测法烧结机漏风故障在线诊断方法 | |
CN107121490A (zh) | 一种湿法冶金过程中运行状态的评价方法 | |
CN113882849B (zh) | 一种多相混输管道危险流型预警及消除系统及其工作方法 | |
CN217724695U (zh) | 一种尾矿车间传统浓密机自动控制系统 | |
CN211420084U (zh) | 一种焦化生产剩余氨水内焦油的深处理装置 | |
CN115526196A (zh) | 基于实时数据及模型的潜油电泵工况在线动态诊断方法 | |
CN114636654A (zh) | 一种通过噪声数据特征实现纳滤膜污染预测的方法 | |
CN109840722B (zh) | 一种金湿法冶金全流程运行状态在线评价方法及系统 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |