CN107728664B - 一种冶金过程多反应器级联pH值优化控制方法 - Google Patents

一种冶金过程多反应器级联pH值优化控制方法 Download PDF

Info

Publication number
CN107728664B
CN107728664B CN201710378120.9A CN201710378120A CN107728664B CN 107728664 B CN107728664 B CN 107728664B CN 201710378120 A CN201710378120 A CN 201710378120A CN 107728664 B CN107728664 B CN 107728664B
Authority
CN
China
Prior art keywords
reactor
value
calcining
concentration
orp
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
Application number
CN201710378120.9A
Other languages
English (en)
Other versions
CN107728664A (zh
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.)
Central South University
Original Assignee
Central South University
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 Central South University filed Critical Central South University
Priority to CN201710378120.9A priority Critical patent/CN107728664B/zh
Publication of CN107728664A publication Critical patent/CN107728664A/zh
Application granted granted Critical
Publication of CN107728664B publication Critical patent/CN107728664B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D21/00Control of chemical or physico-chemical variables, e.g. pH value
    • CCHEMISTRY; METALLURGY
    • C22METALLURGY; FERROUS OR NON-FERROUS ALLOYS; TREATMENT OF ALLOYS OR NON-FERROUS METALS
    • C22BPRODUCTION AND REFINING OF METALS; PRETREATMENT OF RAW MATERIALS
    • C22B19/00Obtaining zinc or zinc oxide
    • C22B19/20Obtaining zinc otherwise than by distilling
    • C22B19/26Refining solutions containing zinc values, e.g. obtained by leaching zinc ores
    • CCHEMISTRY; METALLURGY
    • C25ELECTROLYTIC OR ELECTROPHORETIC PROCESSES; APPARATUS THEREFOR
    • C25CPROCESSES FOR THE ELECTROLYTIC PRODUCTION, RECOVERY OR REFINING OF METALS; APPARATUS THEREFOR
    • C25C1/00Electrolytic production, recovery or refining of metals by electrolysis of solutions
    • C25C1/16Electrolytic production, recovery or refining of metals by electrolysis of solutions of zinc, cadmium or mercury
    • CCHEMISTRY; METALLURGY
    • C25ELECTROLYTIC OR ELECTROPHORETIC PROCESSES; APPARATUS THEREFOR
    • C25CPROCESSES FOR THE ELECTROLYTIC PRODUCTION, RECOVERY OR REFINING OF METALS; APPARATUS THEREFOR
    • C25C7/00Constructional parts, or assemblies thereof, of cells; Servicing or operating of cells
    • C25C7/06Operating or servicing
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P10/00Technologies related to metal processing
    • Y02P10/20Recycling

Landscapes

  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Metallurgy (AREA)
  • Materials Engineering (AREA)
  • Organic Chemistry (AREA)
  • Electrochemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Automation & Control Theory (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Manufacturing & Machinery (AREA)
  • Mechanical Engineering (AREA)
  • Compounds Of Iron (AREA)

Abstract

本发明提供一种冶金过程多反应器级联pH值优化控制方法,针对冶金过程存在的焙砂添加量影响因素多、离子化验滞后性及工况变化频繁导致的焙砂添加量控制困难等问题,通过入口工况划分,利用在线检测获得的ORP和pH值,结合入口溶液流量、Fe2+浓度、Cu2+浓度等信息,根据出口溶液离子浓度建立焙砂添加量优化控制模型,得到pH设定值,并基于模糊控制规则建立焙砂调整量修正模型,对当前焙砂添加量进行修正,使pH达到最优设定值,从而保证针铁矿沉铁过程的稳定运行。

Description

一种冶金过程多反应器级联pH值优化控制方法
技术领域
本发明涉及湿法冶炼技术领域,更具体地,涉及一种冶金过程多反应器级联pH值优化控制方法。
背景技术
湿法炼锌是一个工序流程复杂的冶金过程,包括锌精矿的焙烧、浸出、净化、电解、阴极锌的熔铸等过程。锌精矿焙烧后用电解废液进行中性浸出,使大部分氧化锌溶解,得到的矿浆分离出上清液和底流矿浆。上清液净化后电积产出金属锌,熔铸成锭。底流矿浆进行酸性浸出以溶解残余的氧化锌,酸性浸出液返回到中性浸出;含锌约20%的酸性浸出渣,须进一步处理,传统方法采用回转窑挥发,回收其中的锌、铅和部分稀散金属。
目前在湿法炼锌过程中,硫酸锌溶液中铁离子浓度过高会影响后续净化工序,降低电解过程的电流效率。而针铁矿法沉铁过程中,通过添加氧气及中和剂(如焙砂),使溶液中的Fe2+氧化为Fe3+,同时调节溶液中的pH值,以保证溶液处于合适的pH值范围内,进而使Fe3+水解并以针铁矿的形式沉淀,以控制出口溶液中的铁离子含量处于目标范围内。因此沉铁过程中pH值的控制显得尤为重要,pH值过高或过低都不利于亚铁离子氧化及高铁离子水解。
在实际生产过程中,往往存在一定程度的变化因素,如矿源波动大、入口溶液成分变化大、沉铁反应机理复杂及工况多变等,容易使得溶液pH值波动较大,且由于离子浓度检测信息具有一定的滞后性,从而无法实时且准确地获取溶液中离子浓度的变化情况。
发明内容
本发明的目的是提供一种冶金过程多反应器级联pH值优化控制方法,以解决现有冶金过程中溶液离子浓度的变化情况不能准确实时调控的问题。
根据本发明的技术方案,提供一种冶金过程多反应器级联pH值优化控制方法,该方法包括:
步骤1,根据反应器的入口工况参数,设定多反应器级联pH值;
步骤2,根据pH值和工况参数,建立焙砂添加量数学模型;
步骤3,基于模糊控制规则对焙砂添加量数学模型进行修正,得到pH最优设定值。
作为上述技术方案的进一步改进方案,所述步骤1中,包括:
步骤11,将反应器的入口参数模糊化处理,对工况进行分类;
步骤12,根据模糊化的入口参数,判断当前状态的工况,设定多反应器级联pH值。
所述步骤11中,入口参数包括入口流量、亚铁离子浓度和铜离子浓度信息,对入口流量、亚铁离子浓度、铜离子浓度进行分类,并对ΔO2、ORP、pH值、ΔORP和ΔpH值进行模糊化处理。ΔO2为氧气变化量,ΔORP为ORP变化量,ΔpH为pH变化量。
所述步骤1中,将3个入口工况参数模糊划分区域如下:
入口流量(F)的范围分别为[25,70],[70,125],[125,175](单位为m3/h);
溶液中亚铁离子浓度(cFe2+)的取值范围分别为[0,7.5],[7.5,10],[10,13.5](单位为g/L);
铜离子浓度(cCu2+)的取值范围分别为[0,0.8],[0.8,1.5],[1.5,3](单位为g/L);
上述3个入口工况参数的模糊语言表述如下:
入口流量(F)={Smal l,Medium,Large};
溶液中亚铁离子浓度(cFe2+)={Smal l,Medium,Large};
铜离子浓度(cCu2+)={Small,Medium,Large};
上述5个工况参数的模糊语言变量如下:
ΔO2={H,S,L};
ORP={NB,NS,ZO,PS,PB};
pH={H,S,L};
ΔORP={NB,NS,ZO,PS,PB};
ΔpH={NB,NS,ZO,PS,PB};
其中,NB、NS、ZO、PS和PB的含义依次为“负大”、“负小”、“零”、“正小”和“正大”,H、S和L的含义依次为“偏高区域”、“稳定区域”和“偏低区域”。
步骤11中,通过采用较为平滑的钟形隶属度函数对工况参数模糊处理,钟形隶属度函数为:
所述步骤12中,设定多反应器级联pH值的方法包括:定义当前工况,并初步设定各反应器的pH值决策集;根据各反应器在规定时间内的氧气变化情况,得到反应器的pH值最终决策集。具体为:
Step 1:定义当前第m(m=1,2,…,27)个运行工况描述为:
Um=(C1x,C2y,C3z),其中,C1x表示第x(x=1,2,3)个入口流量条件,C2y表示第y(y=1,2,3)个亚铁离子浓度条件,C3z表示第z(z=1,2,3)个铜离子浓度条件;
Step 2:划分工况后,根据工况信息Um及第i个(i=1,2,3)反应器的ORP初步设定决策集:V=(v1o,v2p,v3q),其中,v1o表示1#反应器的pH设定值(o代表H,S,L),v2p表示2#反应器的pH设定值(p代表H,S,L),v3q表示3#反应器的pH设定值(q代表H,S,L);
Step 3:执行i=i+1;
Step 4:在时间内,判断(i-1)#反应器氧气的变化情况(ΔOi-1);若ΔOi-1=S(t≥30min),跳转Step 6;否则,跳转Step 5;
其中,F为溶液流量;F′为晶种返流流量;ΔOi-1为(i-1)#反应器氧气变化情况,根据实际工业现场数据可知,氧气变化量持续保持30min以上时ORP将会发生显著改变;
Step 5:若ΔOi-1不处于稳定区域,将会导致(i-1)#反应器出口Fe2+发生变化,使得pH值初始决策集中i#、(i+1)#反应器初始设定值发生偏移:
vip′=vip+σ,且令vip=vip′,
viq′=viq+σ,且令viq=viq′;
Step 6:判断i是否大于3,如果i大于3,则结束并输出最终决策集V′=(v1o,v2p,v3q),否则,跳转Step3。
其中,ΔO2为当前反应器氧气与工况预设氧气的差值,O2预设工艺参数如表1所示。
表1 O2预设工艺参数
步骤12中,沉铁过程会产生大量的H+导致溶液中pH值降低,为了满足后续净化和电解工序的安全稳定地运行,需要保证沉铁过程各个反应器的pH值在现场工艺要求范围内,如表2所示。
表2沉铁过程pH值工艺现状
所述步骤2中,针铁矿法沉铁过程中影响焙砂添加量的主要因素为pH值和溶液中铁离子浓度,基于物料平衡和CSTR建立的焙砂添加量数学模型为:
上式代表1#反应器的焙砂添加量数学模型。1#反应器入口流量由两部分组成,即沉铟后液流量F和晶种返流流量F',且各离子均含于沉铟后液中,而1#出口流量则为沉铟流量与入口晶种返流流量之和F+F'。
上式代表2#、3#反应器的焙砂添加量数学模型。2#、3#反应器出入口流量皆为沉铟流量与入口晶种返流流量之和F+F'。
式中各参数含义为:F为溶液流量;F'为晶种返流流量;Qk为第k个反应焙砂流量;MZnO为ZnO原子质量;V为反应器体积;ηA为焙砂中ZnO含量;ηB为ZnO利用率;α、β为反应级数。为第k个反应器H+反应速率;为第k个反应器出口H+浓度;为第k个反应器出口Fe2+浓度,k=1,2,3。
所述步骤3中,由于溶液中的亚铁离子浓度、铁离子浓度等工艺参数化验的大滞后性,使得焙砂添加数学模型的计算值会出现偏差,为了使pH达到最优设定值,故需基于模糊规则对当前焙砂添加量数学模型进行修正,具体包括:根据分类后的工况,通过当前反应器与前一反应器的ORP变化值ΔORP和当前反应器pH实际化验值与预设值的偏差值ΔpH,基于模糊控制规则得到焙砂修正量,通过焙砂对氢离子的中和作用,从而得到pH最优设定值。
以ΔORP与ΔpH值作为输入变量,以焙砂修正量作为输出变量;其中,时序对应后,ΔORP为ORPi与上个反应器的ORPi-1的变化值,ΔpHi值为当前反应器实际化验值(pHi实)与预设值(pHi预)的偏差值,用ΔZnO表示焙砂修正量。
其中,以2#反应器的入口工况为例,其模糊规则表如下所示:
表3模糊规则表
基于上述技术方案,本发明提出一种冶金过程多反应器级联pH值优化控制方法,针对冶金过程存在的焙砂添加量影响因素多、离子化验滞后性及工况变化频繁导致的焙砂添加量控制困难等问题,通过入口工况划分,利用在线检测获得的ORP和pH值,结合入口溶液流量、Fe2+浓度、Cu2+浓度等信息,根据出口溶液离子浓度建立焙砂添加量优化控制模型,得到pH设定值,并基于模糊控制规则建立焙砂调整量修正模型,对当前焙砂添加量进行修正,使pH达到最优设定值,从而保证针铁矿沉铁过程的稳定运行。本发明的控制方法具有以下特点:
克服了目前现场离子浓度检测对人工的依赖性,本发明消除了检测的滞后性,可达到实时调控的目的;
克服了现有的pH值控制方法只将pH值控制在一定范围之内的缺陷,本发明可提供pH的最优设定值;
相比于现有pH控制方法中只从当前反应器出入口离子浓度角度出发进行模型的建立及控制,本发明考虑了入口工况波动对后续反应的影响,提高了pH值优化控制的准确性。
附图说明
图1为根据本发明实施例的沉铁过程的工艺流程;
图2为根据本发明实施例的沉铁过程的系统图;
图3根据本发明实施例的pH值控制规则流程图;
图4为根据本发明实施例的ΔO2与σ关系图;
图5为根据本发明实施例的1#反应器焙砂添加量仿真结果图;
图6为根据本发明实施例的2#反应器焙砂添加量仿真结果图;
图7为根据本发明实施例的3#反应器焙砂添加量仿真结果图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例及其之间任意组合,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明中,氧化还原电位,Oxidation-Reduction Potential,简称ORP,表征介质氧化性或还原性的相对程度。连续搅拌反应器系统,或称连续搅拌槽,ContinuousStirred Tank Reactor,简称CSTR,连续搅拌槽反应器是指带有搅拌浆的槽式反应器,通过搅拌使物料体系达到均匀状态,以有利于反应的均匀和传热。反应过程包括体系中物料的物理和化学的变化,表征其体系特性的参数包括温度、压力、液位及体系组分等。
在根据本申请的一个实施例中,参考图1和图2,多反应器系统包括1#~3#三个相同且相互串联的非封闭连续搅拌槽式反应器,各个反应器之间以溢流槽相连通,呈高到低梯度排列,由沉铟浓密机产生的沉铟后液经溢流槽流入1#反应器中,经1#反应器沉铁反应并溢出后,经溢流槽依次流入2#和3#反应器,以完成整个沉铁过程。整个沉铁过程中3个反应器的温度都保持在80℃左右,故可假设为恒温过程。
在沉铁过程中,分别向三个反应器的底部通入氧气,使溶液中的Fe2+氧化为Fe3+;同时向1#~3#反应器中添加中和剂(焙砂),来调节溶液中的pH值,以保证溶液中的酸性环境,使Fe3+水解并以针铁矿的形式沉淀,从而达到除铁的目的。
在根据本申请的一个实施例中,提供一种冶金过程多反应器级联pH值优化控制方法,该方法包括:
步骤1,根据反应器的入口工况条件的工况参数,设定多反应器级联pH值;
步骤2,根据pH值和工况参数,基于物料平衡和CSTR建立焙砂添加量数学模型;
步骤3,基于模糊控制规则对焙砂添加量数学模型进行修正,得到pH最优设定值。
所述步骤1中,进一步包括:
步骤11,将反应器的入口参数模糊化处理,对工况进行分类;
步骤12,根据模糊化的入口参数,判断当前状态的工况,基于建立的专家规则,设定多反应器级联pH值。
所述步骤11中,由于入口工况实时处于变化波动之中,且入口参数的变化决定了后续反应过程控制参数的调整,故采集入口工况(入口流量、Fe2+浓度、Cu2+浓度)以及各反应器O2、ORP、pH值的信息,对入口流量、Fe2+浓度、Cu2+浓度进行分类,并对ΔO2、ORP、pH值、ΔORP和ΔpH值进行模糊化处理。ΔO2为氧气变化量,ΔORP为ORP变化量,ΔpH为pH变化量。
所述步骤1中,将3个入口工况参数模糊划分区域如下:
入口流量(F)的取值范围分别为[25,70],[70,125],[125,175](单位为m3/h);
溶液中亚铁离子浓度(cFe2+))的取值范围分别为[0,7.5],[7.5,10],[10,13.5](单位为g/L);
铜离子浓度(cCu2+)的取值范围分别为[0,0.8],[0.8,1.5],[1.5,3](单位为g/L);
上述3个入口工况参数的模糊语言表述如下:
入口流量(F)={Smal l,Medium,Large};
溶液中亚铁离子浓度(cFe2+)={Smal l,Medium,Large};
铜离子浓度(cCu2+)={Smal l,Medium,Large};
上述5个工况参数的模糊语言变量如下:
ΔO2={H,S,L};
ORP={NB,NS,ZO,PS,PB};
pH={H,S,L};
ΔORP={NB,NS,ZO,PS,PB};
ΔpH={NB,NS,ZO,PS,PB};
其中,NB、NS、ZO、PS和PB的含义依次为“负大”、“负小”、“零”、“正小”和“正大”,H、S和L的含义依次为“偏高区域”、“稳定区域”和“偏低区域”。
步骤11中,通过采用较为平滑的钟形隶属度函数对工况参数模糊处理,钟形隶属度函数为:
所述步骤12中,基于工况划分设定多反应器级联pH值的方法包括:定义当前工况,并初步设定各反应器的pH值决策集;根据各反应器在规定时间内的氧气变化情况,得到最终决策集。
进一步,参考图3,pH值预设定方法具体包括以下步骤:
Step 1:定义当前第m(m=1,2,…,27)个运行工况描述为:
Um=(C1x,C2y,C3z),
其中,C1x表示第x(x=1,2,3)个入口流量条件,C2y表示第y(y=1,2,3)个亚铁离子浓度条件,C3z表示第z(z=1,2,3)个铜离子浓度条件;
Step 2:划分工况后,根据工况信息Um及第i个(i=1,2,3)反应器的ORP初步设定决策集:
V=(v1o,v2p,v3q),
其中,v1o表示1#反应器的pH设定值(o代表H,S,L),即V1H表示1#反应器的pH设定值为H,V1S表示1#反应器的pH设定值为S,V1L表示1#反应器的pH设定值为L,下同理;v2p表示2#反应器的pH设定值(p代表H,S,L);v3q表示3#反应器的pH设定值(q代表H,S,L);
Step 3:执行i=i+1;
Step 4:在时间内,判断(i-1)#(i=1,2,3)反应器氧气的变化情况(ΔOi-1);若Δoi-1=S(t≥30min),跳转Step 6;否则,跳转Step 5;
其中,F为溶液流量;F′为晶种返流流量;ΔOi-1为(i-1)#反应器氧气变化情况,根据实际工业现场数据可知,氧气变化量持续保持30min以上时ORP将会发生显著改变;
Step 5:若Δoi-1不处于稳定区域,将会导致(i-1)#反应器出口Fe2+发生变化,使得pH值初始决策集中i#、(i+1)#反应器初始设定值发生偏移:
vip′=vip+σ,且令vip=vip
viq′=viq+σ,且令viq=viq
其中,ΔO2与σ关系如图4所示;
Step 6:判断i是否大于3,如果i大于3,则结束并输出最终决策集V′=(v1o,v2p,v3q),否则,跳转Step3。
其中,ΔO2为当前反应器氧气与工况预设氧气的差值,O2预设工艺参数如表1所示。
表1 O2预设工艺参数
步骤12中,沉铁过程会产生大量的H+导致溶液中pH值降低,为了满足后续净化和电解工序的安全稳定地运行,需要保证沉铁过程各个反应器的pH值在现场工艺要求范围内,如表2所示。
表2沉铁过程pH值工艺现状
焙砂作为中和剂可中和沉铁料液中的H+离子,即增加pH值。由勒夏特列原理可得,当pH值变大时,将会促进亚铁氧化反应速率,进一步生成针铁矿。据此建立pH值控制规则如下:
1)IF(U=(C13,C23,C31)and ORP=NS)
THENV′=(v1H,v2H,v3H);
(入口工况为高铁、低铜、大流量,此时各反应器反应时间较短,且入口Fe2+浓度高,极易导致出口铁离子过高不达标;若当前反应器ORP较低时,代表当前反应器氧化程度较低,故应及时提高pH预设值。)
2)IF(U=(C12,C22,C32)and ORP=ZO)THENV′=(v1S,v2S,v3S);
3)IF(U=(C11,C21,C33)and ORP=PB)THENV′=(v1L,V2L,v3L);(入口工况为低铁、高铜、小流量,此时各反应器反应时间较长,且入口Fe2+浓度低,极易导致出口铁离子过低不达标。若当前反应器ORP高时,代表当前反应器氧化程度大,故应及时降低pH预设值。)
在针铁矿法沉铁过程中,反应器中主要存在着氧化反应、水解反应和中和反应,可用以下公式简化表示:
氧化反应:4Fe2++4H++O2→4Fe3++2H2O (1)
水解反应:Fe3++2H2O→FeOOH↓+3H+ (2)
中和反应:2H++ZnO→Zn2++H2O (3)
由式(1)可知亚铁离子氧化消耗少量的氢离子,但由式(2)可知三价铁离子水解会产生大量的氢离子,其总反应式如下式所示:
4Fe2++O2+6H2O→4FeOOH↓+8H+ (4)
由式(3)、(4)可知,影响焙砂添加量的因素为pH值及溶液中铁离子浓度。在针铁矿法沉铁过程中,需要加入中和剂(焙砂)来调节溶液中的pH值,因此需建立焙砂添加量数学模型,如下所示:
上式为公式(5),代表1#反应器的焙砂添加量数学模型。1#反应器入口流量由两部分组成,即沉铟后液流量F和晶种返流流量F',且各离子均含于沉铟后液中,而1#出口流量则为沉铟流量与入口晶种返流流量之和F+F'。
上式为公式(6),代表2#、3#反应器的焙砂添加量数学模型。2#、3#反应器出入口流量皆为沉铟流量与入口晶种返流流量之和F+F'。
式(5)、(6)中各参数含义为:F为溶液流量;F'为晶种返流流量;Qk为第k个反应焙砂流量;MZnO为ZnO原子质量;V为反应器体积,V=250m3;ηA为焙砂中ZnO含量;ηB为ZnO利用率;α、β为反应级数。为第k个反应器H+反应速率;为第k个反应器出口H+浓度;为第k个反应器出口Fe2+浓度,k=1,2,3。
由于溶液中的亚铁离子浓度、铁离子浓度等工艺参数化验的大滞后(每8h化验2次),影响了焙砂的准确添加,故焙砂添加数学模型的计算值会出现偏差,为了使pH达到最优设定值,需要对焙砂添加数学模型计算的焙砂添加量进行修正。模糊控制能够高效率、准确和迅速地解决实际问题。因此,建立数学模型与模糊专家规则修正模型相结合的焙砂添加控制方法。
对焙砂添加量数学模型进行修正,具体包括:根据分类后的工况,通过当前反应器与前一反应器的ORP变化值ΔORP和当前反应器pH实际化验值与预设值的偏差值ΔpH,基于模糊控制规则得到焙砂修正量。
以ΔORP与ΔpH值作为输入变量,以焙砂修正量作为输出变量;其中,时序对应后,ΔORP为ORPi与上个反应器的ORPi-1的变化值,ΔpHi值为当前反应器实际化验值(pHi实)与预设值(pHi预)的偏差值,用ΔZnO表示焙砂修正量。
以2#反应器入口工况为例,输入变量选择ΔORP与ΔpH值,对输出变量焙砂调整量(ΔZnO)建立如下的模糊规则表:
表3模糊规则表
其中:
1)ifΔORP=PB,ΔpH=PB,ThenΔZnO=NB。(当前反应器ORP值远大于上个反应器ORP值,表明反应器氧化状态过强;且pH值的实际化验值远大于预设值,故表明上个时刻焙砂添加量较多,应大大减少焙砂添加量。)
2)ifΔORP=PB,ΔpH=NB,ThenΔZnO=PS。(当前反应器ORP值远大于上个反应器ORP值,表明ORP处于最优范围内;pH值的实际化验值远小于预设值,故表明上个时刻焙砂添加量较少,应适量增加焙砂添加量。)
3)ifΔORP=ZO,ΔpH=ZO,ThenΔZnO=PS。(当前反应器ORP值与上个反应器ORP值相同,表明反应器氧化状态没有发生变化,ORP变化值较小;pH值的实际化验值与预设值相同,故表明上个时刻焙砂添加量适中;但由于氧化状态较弱,故应适当增加焙砂添加量。)
4)ifΔORP=NB,ΔpH=PB,ThenΔZnO=ZO。(当前反应器ORP值远小于上个反应器ORP值,表明当前反应器氧化状态很弱;pH值的实际化验值远大于预设值,故表明上个时刻焙砂添加量较多;但由于氧化状态很弱,故应保持焙砂添加量不变。)
5)ifΔORP=NB,ΔpH=NB,ThenΔZnO=PB。(当前反应器ORP值远小于上个反应器ORP值,表明当前反应器氧化状态很弱;pH值的实际化验值远小于预设值,故表明上个时刻焙砂添加量较少;应大大增加焙砂添加量。)
利用这些现场数据对模型进行仿真,得出如下结论:本发明实施例所提出的综合考虑入口工况波动等因素预设定pH值,根据沉铁过程化学反应机理和物料平衡建立焙砂添加量模型,由于沉铁过程化学机理复杂,单纯机理模型可能会对其产生偏差,故提出了基于模糊规则的焙砂添加量修正模型。图5、6、7分别为1#、2#、3#反应器焙砂添加量的仿真结果,图中左侧为人工操作值,右侧为模型输出值,该方法有效地减小了pH值波动范围,为实现针铁矿沉铁过程的稳定运行提供了条件。
本发明实施例提供了一种冶金过程多反应器级联pH值优化控制方法,通过入口工况划分,利用在线检测获得的ORP和pH值,结合入口溶液流量、Fe2+浓度、Cu2+浓度等信息,根据出口溶液离子浓度建立焙砂添加量优化控制模型,得到pH设定值,并基于模糊控制规则建立焙砂调整量修正模型,对当前焙砂添加量进行修正,使pH达到最优设定值,从而保证针铁矿沉铁过程的稳定运行。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (7)

1.一种冶金过程多反应器级联pH值优化控制方法,该方法包括:
步骤1,根据反应器的入口工况参数,设定多反应器级联pH值;
步骤2,根据pH值和工况参数,建立焙砂添加量数学模型;
步骤3,基于模糊控制规则对焙砂添加量数学模型进行修正,得到pH最优设定值;
所述步骤1中,包括:
步骤11,将反应器的入口参数模糊化处理,对工况进行分类;
步骤12,根据模糊化的入口参数,判断当前状态的工况,设定多反应器级联pH值;
所述步骤12中,多反应器级联pH值的设定方法包括:定义当前工况,并初步设定各反应器的pH值决策集,根据各反应器在规定时间内的氧气变化情况,得到反应器的pH值最终决策集,具体为:
Step 1:定义当前第m,m=1,2,…,27个运行工况描述为:
Um=(C1x,C2y,C3z),其中,C1x表示第x,x=1,2,3个入口流量条件,C2y表示第y,y=1,2,3个亚铁离子浓度条件,C3z表示第z,z=1,2,3个铜离子浓度条件;
Step 2:划分工况后,根据工况信息Um及第i个i=1,2,3反应器的ORP初步设定决策集:V=(v1o,v2p,v3q),其中,v1o表示1#反应器的pH设定值,o代表H,S,L,v2p表示2#反应器的pH设定值,p代表H,S,L,v3q表示3#反应器的pH设定值,q代表H,S,L,H、S和L的含义依次为“偏高区域”、“稳定区域”和“偏低区域”;
Step 3:当i=2时,在时间内,判断1#反应器氧气的变化情况ΔO1;若ΔO1=S,t≥30min,跳转Step 5;否则,跳转Step 4;
其中,F为溶液流量;F′为晶种返流流量;V为反应器体积;ΔO1为1#反应器氧气变化情况;
Step 4:若ΔO1未处于稳定区域,将导致1#反应器出口Fe2+发生变化,使得pH值初始决策集中2#、3#反应器初始设定值发生偏移:
v2p′=v2p+σ,且令v2p=v2p′,
v3q′=v3q+σ,且令v3q=v3q′;
其中,v2p′表示2#反应器发生偏移后的pH设定值,v3q′表示3#反应器发生偏移后的pH设定值,σ为氧气偏移量;
Step 5:当i=3时,在时间内,判断2#反应器氧气的变化情况ΔO2;若ΔO2=S,t≥30min,跳转Step 7;否则,跳转Step 6;
其中,F为溶液流量;F′为晶种返流流量;V为反应器体积;ΔO2为2#反应器氧气变化情况;
Step 6:若ΔO2未处于稳定区域,将导致2#反应器出口Fe2+发生变化,使得pH值初始决策集中3#反应器初始设定值发生偏移:
v3q′=v3q+σ,且令v3q=v3q′;
其中,v3q′表示3#反应器发生偏移后的pH设定值,σ为氧气偏移量;
Step 7:结束并输出最终决策集V′=(v1o,v2p,v3q)。
2.根据权利要求1所述的pH值优化控制方法,其特征在于,所述步骤11中,入口参数包括入口流量、亚铁离子浓度和铜离子浓度信息,对入口流量、亚铁离子浓度、铜离子浓度进行分类,并对ΔO2、ORP、pH值、ΔORP和ΔpH值进行模糊化处理。
3.根据权利要求2所述的pH值优化控制方法,其特征在于,工况参数的模糊语言变量如下:
入口流量(F)={Small,Medium,Large},Small、Medium和Large分别表示入口流量(F)取值为[25,70]、[70,125]和[125,175],单位为m3/h;
溶液中亚铁离子浓度(cFe2+)={Small,Medium,Large},Small、Medium和Large分别表示溶液中亚铁离子浓度(cFe2+)取值为[0,7.5]、[7.5,10]和[10,13.5],单位为g/L;
铜离子浓度(cCu2+)={Small,Medium,Large},Small、Medium和Large分别表示铜离子浓度(cCu2+)取值为[0,0.8]、[0.8,1.5]和[1.5,3],单位为g/L;
ΔO2={H,S,L};
ORP={NB,NS,ZO,PS,PB};
pH={H,S,L};
ΔORP={NB,NS,ZO,PS,PB};
ΔpH={NB,NS,ZO,PS,PB};
其中,NB、NS、ZO、PS和PB的含义依次为“负大”、“负小”、“零”、“正小”和“正大”,H、S和L的含义依次为“偏高区域”、“稳定区域”和“偏低区域”。
4.根据权利要求3所述的pH值优化控制方法,其特征在于,步骤11中,通过采用钟形隶属度函数对工况参数模糊处理。
5.根据权利要求1所述的pH值优化控制方法,其特征在于,所述步骤2中,基于物料平衡和CSTR,建立反应器的焙砂添加量数学模型。
6.根据权利要求5所述的pH值优化控制方法,其特征在于,所述焙砂添加量数学模型为:
式中,F为溶液流量;F'为晶种返流流量;Qk为第k个反应焙砂流量;MZnO为ZnO原子质量;V为反应器体积;ηA为焙砂中ZnO含量;ηB为ZnO利用率;α、β为反应级数;rH+,1为第一个反应器H+反应速率;cH+,0为第一个反应器入口H+浓度;cFe2+,0为第一个反应器入口Fe2+浓度;cFe2+,1为第一个反应器出口Fe2+浓度;cH+,1为第一个反应器出口H+浓度;为第k个反应器H+反应速率;为第k个反应器出口H+浓度;为第k个反应器出口Fe2+浓度,k=2,3。
7.根据权利要求1所述的pH值优化控制方法,其特征在于,所述步骤3中,根据分类后的工况,通过当前反应器与前一反应器的ORP变化值ΔORP和当前反应器pH实际化验值与预设值的偏差值ΔpH,基于模糊控制规则得到焙砂修正量,从而得到pH最优设定值。
CN201710378120.9A 2017-05-25 2017-05-25 一种冶金过程多反应器级联pH值优化控制方法 Active CN107728664B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710378120.9A CN107728664B (zh) 2017-05-25 2017-05-25 一种冶金过程多反应器级联pH值优化控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710378120.9A CN107728664B (zh) 2017-05-25 2017-05-25 一种冶金过程多反应器级联pH值优化控制方法

Publications (2)

Publication Number Publication Date
CN107728664A CN107728664A (zh) 2018-02-23
CN107728664B true CN107728664B (zh) 2019-10-29

Family

ID=61201235

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710378120.9A Active CN107728664B (zh) 2017-05-25 2017-05-25 一种冶金过程多反应器级联pH值优化控制方法

Country Status (1)

Country Link
CN (1) CN107728664B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109813782B (zh) * 2019-02-01 2020-12-29 杭州原创软件有限公司 一种自诊断ph在线监测系统及方法
CN110672792B (zh) * 2019-10-09 2020-08-28 中南大学 湿法炼锌中性浸出过程pH值软测量方法和系统
CN115081305B (zh) * 2022-04-18 2024-08-23 中南大学 多工况级联冶金反应过程氧化还原电位优化设定方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1993015236A1 (fr) * 1992-01-28 1993-08-05 N.V. Union Miniere S.A. Procede pour separer du fer d'une solution aqueuse
CN1122258A (zh) * 1994-07-11 1996-05-15 三菱重工业株式会社 控制亚硫酸盐氧化的方法
CN102010994A (zh) * 2010-12-29 2011-04-13 株洲冶炼集团股份有限公司 一种湿法炼锌过程中高酸高铁溶液针铁矿沉铁方法
CN105525093A (zh) * 2015-12-25 2016-04-27 云南云铜锌业股份有限公司 一种简化湿法炼锌工艺的方法
CN105950866A (zh) * 2016-04-29 2016-09-21 中南大学 一种基于氧化还原电位的锌液净化除铜过程优化控制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1993015236A1 (fr) * 1992-01-28 1993-08-05 N.V. Union Miniere S.A. Procede pour separer du fer d'une solution aqueuse
CN1122258A (zh) * 1994-07-11 1996-05-15 三菱重工业株式会社 控制亚硫酸盐氧化的方法
CN102010994A (zh) * 2010-12-29 2011-04-13 株洲冶炼集团股份有限公司 一种湿法炼锌过程中高酸高铁溶液针铁矿沉铁方法
CN105525093A (zh) * 2015-12-25 2016-04-27 云南云铜锌业股份有限公司 一种简化湿法炼锌工艺的方法
CN105950866A (zh) * 2016-04-29 2016-09-21 中南大学 一种基于氧化还原电位的锌液净化除铜过程优化控制方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Optimal control for zinc solution purification based on interacting CSTR models;Yong Gang Li etc;《Journal of Process Control》;20120910;全文 *
杨根.针铁矿法沉铁过程pH值的稳定控制及应用.《中国优秀硕士学位论文全文数据库•信息科技辑》.2014,(第5期), *
湿法炼锌过程中pH值的预测控制;何献忠;《电气自动化》;20090228;第31卷(第1期);全文 *

Also Published As

Publication number Publication date
CN107728664A (zh) 2018-02-23

Similar Documents

Publication Publication Date Title
CN107728664B (zh) 一种冶金过程多反应器级联pH值优化控制方法
CN103276205B (zh) 一种钒铬浸出液中分离提取钒铬的方法
Han et al. Dynamic multi-objective optimization arising in iron precipitation of zinc hydrometallurgy
Sun et al. An integrated prediction model of cobalt ion concentration based on oxidation–reduction potential
Sun et al. Modeling, optimization, and control of solution purification process in zinc hydrometallurgy
Xie et al. Hybrid fuzzy control for the goethite process in zinc production plant combining type-1 and type-2 fuzzy logics
CN102732723B (zh) 一种锌湿法冶炼砷盐净化除钴过程废酸添加控制方法
CN112899496A (zh) 湿法炼锌净化除钴过程入口铜离子浓度估计方法及系统
Yang et al. Modeling, optimization, and control of zinc hydrometallurgical purification process
CN110675918B (zh) 一种沉铁过程出口离子预测方法及其系统
CN114420220A (zh) 一种湿法炼锌中性浸出过程pH值的预测方法
Niu et al. Optimization control of hydrometallurgical leaching process based on IDE-MPC
Zhu et al. A fuzzy control method based on rule extraction for zinc leaching process of zinc hydrometallurgy
CN115202190B (zh) 一种基于工艺指标调节潜质的工况划分和识别方法
Zhang et al. Optimal setting strategy of electrocoagulation process in heavy metal wastewater treatment plant
CN105950866B (zh) 一种基于氧化还原电位的锌液净化除铜过程优化控制方法
Gu et al. Highly oriented (002) crystals zinc powder from recovery of zinc smelting slag
CN107657143B (zh) 液体中氧气通入量的设置方法
Zhang et al. Control strategy for hydrometallurgical removal process based on modelling and evaluation
Min et al. Intelligent optimal control for lead-zinc sintering process state
CN114496110A (zh) 湿法炼锌净化除钴反应器出口钴离子浓度预测方法和系统
CN112859597A (zh) 一种基于数据驱动的沉铁过程优化控制方法及其系统
CN105740605A (zh) 一种湿法炼锌除铜过程生产工况评估方法
CN115081305B (zh) 多工况级联冶金反应过程氧化还原电位优化设定方法
AU2021368016B2 (en) Cobalt extraction

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