沉铁过程出口离子预测方法及系统

技术领域

本发明涉及湿法炼锌沉铁过程控制领域，尤其涉及一种沉铁过程出口离子预测方法及系统。

背景技术

锌是一种能应用在各个领域的重要金属。湿法冶炼由于高产出和低能耗成为最主要的锌冶炼方法，主要包括磨矿、浸出、净化、电解等工序。该方法中的硫酸锌溶液需要经过净化除杂后，方可进行电解得到锌单质。目前，针铁矿法常用来去除硫酸锌溶液中的主要杂质铁离子。沉铁过程的主要生产设备为连续搅拌釜式反应器(CSTR)，然而单个反应器不能将浸出液中过量铁离子直接降低到工艺要求范围内，因此需要在四个由高到低级联的反应器中进行除铁。前一反应器的出口离子浓度为下一个反应器的入口离子浓度，各个反应器出口铁离子浓度需按要求依次降低，保证硫酸锌溶液离开最后一个反应器后其中铁离子含量降低至工艺指标要求范围内。为了达到这个目的，需根据每个反应器入口离子浓度(即上一反应器的出口离子浓度)分别对其中添加的氧气和焙砂进行调节，然而在实际沉铁过程中，反应器自身的密闭性及检测装置的局限性使得反应器出口溶液中各离子的浓度只能通过人工定期抽样检测获得。导致对操作参数(氧气和焙砂)的调节具有很大的滞后性，降低沉铁效率。因此，建立沉铁过程模型，预测反应器出口离子浓度，对调整操作参数的盲目性、提高沉铁效率具有重大意义。

实际沉铁过程中，合理控制Fe

因此，现有的沉铁过程出口离子浓度预测模型中难以检测溶解氧浓度，造成预测的出口离子浓度精确度不高已成为本领域技术人员亟待解决的技术问题。

发明内容

本发明目的在于公开一种沉铁过程出口离子预测方法及系统，以解决现有的沉铁过程出口离子浓度预测模型中难以检测溶解氧浓度，造成预测的出口离子浓度精确度不高的技术问题。

为达上述目的，本发明公开一种沉铁过程出口离子预测方法，包括以下步骤：

S1：构建影响沉铁过程出口离子浓度效果的特征向量，获取有标签数据样本集，采用仿射聚类传播算法得到c个有代表性的中心样本点，分别围绕c个中心样本选取相似度高的L个时间连续样本构成c个子训练样本集；

S2：将对应不同工况的各子训练样本集分别进行局部邻域标准化处理，然后以处理后的各子训练样本集数据分别对c个动态概率隐变量子模型参数进行辨识，以及根据EM算法求得各动态概率隐变量子模型的参数，所述EM算法利用现有的子训练样本集数据信息并使用矩估计对初值进行选择；

S3：基于S1所构建的特征向量获取用于预测离线沉铁过程出口离子浓度特征的无标签的测试样本数据，对同一测试样本数据进行局部邻域标准化处理后输入到各动态概率隐变量子模型得到出口离子浓度预测值的概率分布；

S4:采用贝叶斯决策进行联合估计将c个动态概率隐变量子模型输出的预测值的概率分布结果进行融合，得到对应所述测试样本数据最终的出口离子浓度的预测值。

优选地，上述S1所构建的特征向量x

其中，从左到右分别为：二价铁离子浓度、三价铁离子浓度、铜离子浓度、锌离子浓度、液体流速、FeOOH的质量、氧化锌的质量和通入的氧气的流量；

预测的沉铁过程出口离子浓度向量y

其中，下标中的in和out分别代表输入和输出，

优选地，S1具体包括：

对数据集{X,Y}＝{(x

分别根据聚类中心v

其中，d

可选的，本发明局部邻域标准化处理的计算公式为：

其中，局部领域N

优选地，S2具体包括：

以H

H

x

y

其中，A

基于隐变量H

根据得到的第k个工况样本集X

采用EM算法求解似然函数J

其中，第k个动态概率隐变量子模型对应的参数为Θ

利用现有的子训练样本集数据信息{X

E

E

E

得到第k个动态概率隐变量子模型EM算法参数的初始值如下：

A

μπ

P

Σ

Σ

利用EM算法的M步求解出优化的模型参数为：

μ

为求解优化的参数，EM算法的E步需要求得隐变量的一阶和二阶条件数学期望，E(H

当t＝1时有p(H

当t≥2时有p(H

然后使用后向平滑算法推导出每个时刻隐变量的后验概率分布H

m

M

类似地，当t＝L时，p(H

最后，根据以上概率分布可得到M步所需的隐变量的一阶和二阶期望如下：

E(H

通过EM算法反复迭代E步和M步，直到参数矩阵Θ

优选地，S3具体包括：

首先需求出初始时刻，即t＝1时的测试样本隐变量后验分布

得到初始时刻测试样本隐变量概率函数的均值和方差分别为：

在已知输入变量序列x

类似地，得到参数：

根据贝叶斯法则可得到，隐变量关于测试样本输入序列

当t＝1时，

当t≥2时，

其中，下标m表示为测试样本的数量，上标q表示测试样本；最后，第t个测试样本序列

优选地，得到各动态概率隐变量子模型的预测输出为{y

对新的预测样本x

通过BCM算法得到的全局输出也是高斯分布的，其均值和方差如下：

得到的均值E(y

为达上述目的，本发明还公开一种沉铁过程出口离子预测系统，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其中，所述处理器执行所述计算机程序时实现上述相对应方法的步骤。

本发明具有以下有益效果：

对于复杂的沉铁工业过程，由于传感器采集到的过程数据与容器的出口离子浓度存在复杂的未知函数关系，不能单单采用机理模型来得到过程函数关系。但是过程能够产生很多包含过程信息的数据，本发明采用数据驱动的软测量建模方法来预测沉铁过程出口离子浓度，充分考虑了沉铁过程存在多工况、过程的非线性、高噪音以及过程存在动态特性等，基于贝叶斯决策混合动态概率隐变量模型的沉铁过程出口离子预测方法，通过工业过程数据得到沉铁过程出口离子预测模型，提高模型的可信度。相比现有技术而言，使用本发明中的沉铁过程出口离子预测模型预测出的出口离子浓度准确性更高，能达到高效除铁的作用。

下面将参照附图，对本发明作进一步详细的说明。

附图说明

构成本申请的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：

图1是本发明实施例的沉铁过程出口离子预测方法流程示意图。

图2是对沉铁过程进行机理分析采用的单个反应器的CSTR系统示意图。

图3是本发明实施例预测出口离子浓度的动态概率隐变量模型结构图。

图4是本发明实施例BCM算法结构示意图。

图5是本发明实施例方法预测值与真实值的对比示意图。

具体实施方式

以下结合附图对本发明的实施例进行详细说明，但是本发明可以由权利要求限定和覆盖的多种不同方式实施。

实施例1

一种沉铁过程出口离子预测方法，如图1所示，包括以下步骤：

S1：构建影响沉铁过程出口离子浓度效果的特征向量，获取有标签数据样本集，采用仿射聚类传播算法得到c个有代表性的中心样本点，分别围绕c个中心样本选取相似度高的L个时间连续样本构成c个子训练样本集。

以中国某锌冶炼企业中实际生产数据为例，沉铁过程中的测量变量及其采样频率如表1所示。

表1：

气体溶于液体的过程就是气体分子转移到液相主体的过程。影响气体在液体中的溶解度的因素通常有压强、温度、液体流速、溶液中离子浓度、溶液中固体含量、机械搅拌等。压强和温度变化时，液相主体间隙的变化使气体溶解度改变；通入氧气速率、溶液流速、机械搅拌通过影响与液面接触的氧气分子数目，改变溶解氧浓度。溶液中离子浓度引起溶液极性变化，氧气作为非极性分子的溶解度会发生改变。此外，沉铁过程中存在反应的动态平衡，其中针铁矿的存在对氧气溶解度也有影响。

如图2所示，以单个反应器为对象，分析针铁矿法沉铁过程影响氧气溶解过程的因素。设

氧化反应：4Fe

水解反应：

中和反应：2H

实际沉铁过程中，压强和温度的影响忽略不计，机械搅拌速度恒定。小试法中的溶解氧浓度并不都属于同一时刻，也并非每个时刻都能采到，但存在包含上一采样时刻溶解氧浓度的过程数据及同时包含当前采样时刻和上一采样时刻溶解氧浓度的过程数据。因此，为了考虑过程的动态特性，还利用其前采样时刻的输入变量来预测当前采样时刻的输出变量。因此，浸出液中溶解氧浓度的动态模型可以描述成：

y

在工业过程中采集连续n个有标签的数据集{X,Y}＝{(x

之后，用仿射聚类传播算法对数据集{X,Y}＝{(x

仿射传播聚类算法可以描述为：在初始时刻，每个样本点都被当成可能存在的聚类中心，然后通过相似度计算，不断进行聚类过程循环迭代，每个数据点去竞争当选聚类中心，直到找到最优的聚类中心为止。

基于上述采集得到的数据集{X}＝{(x

R(i,j)＝S(i,j)-max{A(i,k)+S(i,k)}(k＝1,...,n,k≠j)

A(i,j)＝min{0,R(j,J)+∑{max(0,R(k,j))}}(k＝1,...,n,k≠i,j)

其中，证据可信度R(i,j)用来描述x

R

A

仿射聚类传播算法通过以上公式进行循环更新，对各个证据进行搜索、传递和整理，λ来控制算法的迭代更新速度，对于数据点x

利用聚类中心在数据集{X,Y}中找到样本大小为L的连续子数据集{X

对每个聚类中心v

其中，d

在该步骤中，优选的局部邻域标准化处理的计算公式为：

其中，局部领域N

S2：将对应不同工况的各子训练样本集分别进行局部邻域标准化处理，然后以处理后的各子训练样本集数据分别对c个动态概率隐变量子模型参数进行辨识，以及根据EM算法求得各子动态概率隐变量子模型的参数，所述EM算法利用现有的子训练样本集数据信息并使用矩估计对初值进行选择。

优选地，参考图3所示的预测出口离子浓度的动态概率隐变量模型结构，以H

H

x

y

其中，A

基于隐变量H

根据得到的第k个工况样本集X

采用EM算法求解似然函数J

其中，第k个动态概率隐变量子模型对应的参数为Θ

利用现有的子训练样本集数据信息{X

E

E

E

得到第k个动态概率隐变量子模型EM算法参数的初始值如下：

A

μ

P

Σ

Σ

利用EM算法的M步求解出优化的模型参数为：

μ

为求解优化的参数，EM算法的E步需要求得隐变量的一阶和二阶条件数学期望，E(H

当t＝1时有p(H

当t≥2时有p(H

然后使用后向平滑算法推导出每个时刻隐变量的后验概率分布H

m

M

类似地，当t＝L时，p(H

最后，根据以上概率分布可得到M步所需的隐变量的一阶和二阶期望如下：

E(H

在该步骤中，由于传统的EM算法求解的不是全局最优解，故本实施例对其进行了改进，通过利用现有的子训练样本集数据信息并使用矩估计对初值进行选择，有效避免了陷入局部最优。

通过EM算法反复迭代E步和M步，直到参数矩阵Θ

S3：基于S1所构建的特征向量获取用于预测离线沉铁过程出口离子浓度特征的无标签的测试样本数据，对同一测试样本数据进行局部邻域标准化处理后输入到各动态概率隐变量子模型得到出口离子浓度预测值的概率分布。

优选地，该步骤具体包括：

首先需求出初始时刻，即t＝1时的测试样本隐变量后验分布

得到初始时刻测试样本隐变量概率函数的均值和方差分别为：

在已知输入变量序列x

类似地，得到参数：

根据贝叶斯法则可得到，隐变量关于测试样本输入序列

当t＝1时，

当t≥2时，

其中，下标m表示为测试样本的数量，上标q表示测试样本；最后，第t个测试样本序列

S4:采用贝叶斯决策进行联合估计将c个动态概率隐变量子模型输出的预测值的概率分布结果进行融合，得到对应所述测试样本数据最终的出口离子浓度的预测值。

得到各动态概率隐变量子模型的预测输出为{y

对新的预测样本x

通过BCM算法得到的全局输出也是高斯分布的，其均值和方差如下：

得到的均值E(y

基于本实施例贝叶斯决策混合动态概率隐变量模型的沉铁过程出口离子预测方法构造模型预测值与真实值对比图参照图5；由图可知，本实施例预测值与实测值基本吻合，准确度可靠。

实施例2

一种沉铁过程出口离子预测系统，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其中，所述处理器执行所述计算机程序时实现上述方法所对应的各步骤。

综上，本发明上述各实施例所分别公开的沉铁过程出口离子预测方法及系统，至少具有以下有益效果：

对于复杂的沉铁工业过程，由于传感器采集到的过程数据与容器的出口离子浓度存在复杂的未知函数关系，不能单单采用机理模型来得到过程函数关系。但是过程能够产生很多包含过程信息的数据，本发明采用数据驱动的软测量建模方法来预测沉铁过程出口离子浓度，充分考虑了沉铁过程存在多工况、过程的非线性、高噪音以及过程存在动态特性等，基于贝叶斯决策混合动态概率隐变量模型的沉铁过程出口离子预测方法，通过工业过程数据得到沉铁过程出口离子预测模型，提高模型的可信度。相比现有技术而言，使用本发明中的沉铁过程出口离子预测模型预测出的出口离子浓度准确性更高，能达到高效除铁的作用。

以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。