[发明专利]基于混合寻优的丙烯聚合生产过程最优软测量系统和方法

权利要求书

1.一种基于混合寻优的丙烯聚合生产过程最优软测量系统，包括丙烯聚合生产过程、用于测量易测变量的现场智能仪表、用于测量操作变量的控制站、存放数据的DCS数据库、基于混合寻优的最优软测量系统以及熔融指数软测量显示仪，所述现场智能仪表、控制站与丙烯聚合生产过程连接，所述现场智能仪表、控制站与DCS数据库连接，所述DCS数据库与基于混合寻优的的最优软测量模型的输入端连接，所述基于混合寻优的最优软测量系统的输出端与熔融指数软测量显示仪连接，其特征在于：所述基于混合寻优的最优软测量系统包括：

（1）、数据预处理模块，用于将从DCS数据库输入的模型输入变量进行预处理，对输入变量中心化，即减去变量的平均值；再进行归一化处理，即除以变量值的变化区间；

（2）、PCA主成分分析模块，用于将输入变量预白化处理及变量去相关，通过对输入变量施加一个线性变换实现，即主成分由C=MU得到，其中M为输入变量，C为主成分得分矩阵，U为载荷矩阵。若对原始数据进行重构，可由M=CU^T计算，其中上标T表示矩阵的转置。当选取的主成分数目小于输入变量的变量个数时，M=CU^T+E，其中E为残差矩阵；

（3）、神经网络模型模块，用于采用RBF神经网络、通过误差函数最小化来完成输入到输出的一种高度非线性映射，映射中保持拓扑不变性；

（4）、混合寻优优化模块，用于采用基于粒子群-模拟退火（HPSO-SA）混合算法的优化模块对神经网络进行优化，包括：

（4.1）算法初始化，根据待优化的RBF神经网络结构构造出初始的粒子群体X=(x¹,x²,…,x^N)，初始移动速度V=(v¹,v²,…,v^N)，初始的各粒子历代最优值OP=(p¹,p²,…,p^N)和全局最优值p^g；

（4.2）通过下式执行HPSO算法，让粒子群收敛：

xk+1i=xki+vk+1i---(1)]]>

vk+1i=wkvki+c1r1(pki-xki)+c2r2(plg-xki)---(2)]]>

式中x为粒子的位置向量，i为粒子的序号，k为算法迭代代数，v为粒子速度，p表示初始的各粒子历代最优值OP=(p¹,p²,…,p^N)和全局最优值p^g的最优值集合。v_kⁱ为第k次迭代代数中第i个粒子的速度；x_kⁱ第k次迭代代数中第i个粒子的位置向量；p_kⁱ为第k次迭代代数中第i个粒子的历代最优位置，pⁱ为第i个粒子的历代最优解，w为速度权重系数，c₁、c₂分别为粒子历代最优解和群体最优解的吸引系数，r₁、r₂分别为随机数。

（4.3）待粒子群收敛到合适的程度时，开始对各个粒子的历代最优解组成的集合(best¹,best²,…,best^N)执行SA算法，得到的最优解作为算法的最后结果返回。

所述基于混合寻优的最优软测量系统还包括：模型更新模块，用于模型的在线更新，定期将离线化验数据输入到训练集中，更新神经网络模型。

所述的基于混合寻优的的丙烯聚合生产过程最优软测量系统，其特征在于：PCA主成分分析模块中，PCA方法实现输入变量的预白化处理，能够简化神经网络模型的输入变量，进而提高模型的性能。

2.一种用如权利要求1所述的基于混合寻优的的聚丙烯生产过程最优软测量系统实现的软测量方法，其特征在于：所述软测量方法具体实现步骤如下：

(4.1)对丙烯聚合生产过程对象，根据工艺分析和操作分析，选择操作变量和易测变量作为模型的输入，一般操作变量和易测变量取温度、压力、液位、氢气气相百分数、3股丙烯进料流速和2股催化剂进料流速这些变量，由DCS数据库获得；

(4.2)对样本数据进行预处理，对输入变量中心化，即减去变量的平均值；再进行归一化处理，即除以变量值的变化区间；

(4.3)PCA主成分分析模块，用于将输入变量预白化处理及变量去相关，通过对输入变量施加一个线性变换实现，即主成分由C=MU得到，其中M为输入变量，C为主成分得分矩阵，U为载荷矩阵。若对原始数据进行重构，可由M=CU^T计算，其中上标T表示矩阵的转置。当选取的主成分数目小于输入变量的变量个数时，M=CU^T+E，其中E为残差矩阵；

(4.4)基于模型输入、输出数据建立初始神经网络模型，采用RBF神经网络，通过误差最小化来完成输入到输出的一种高度非线性映射，映射中保持拓扑不变性；

(4.5)采用基于HPSO-SA算法的优化模块对神经网络进行优化，包括：

（4.5.1）算法初始化，根据待优化的RBF神经网络结构构造出初始的粒子群体X=(x¹,x²,…,x^N)，初始移动速度V=(v¹,v²,…,v^N)，初始的各粒子历代最优值OP=(p¹,p²,…,p^N)和全局最优值p^g；

（4.5.2）通过下式执行HPSO算法，让粒子群收敛：

xk+1i=xki+vk+1i---(1)]]>

vk+1i=wkvki+c1r1(pki-xki)+c2r2(plg-xki)---(2)]]>

式中x为粒子的位置向量，i为粒子的序号，k为算法迭代代数，v为粒子速度，p表示初始的各粒子历代最优值OP=(p¹,p²,…,p^N)和全局最优值p^g的最优值集合。v_kⁱ为第k次迭代代数中第i个粒子的速度；x_kⁱ第k次迭代代数中第i个粒子的位置向量；p_kⁱ为第k次迭代代数中第i个粒子的历代最优位置，pⁱ为第i个粒子的历代最优解，w为速度权重系数，c₁、c₂分别为粒子历代最优解和群体最优解的吸引系数，r₁、r₂分别为随机数。

（4.5.3）待粒子群收敛到合适的程度时，开始对各个粒子的历代最优解组成的集合(best¹,best²,…,best^N)执行SA算法，得到的最优解作为算法的最后结果返回。

所述软测量方法还包括：定期将离线化验数据输入到训练集中，更新神经网络模型。

所述的软测量方法，其特征在于：在所述的步骤（4.3）中采用PCA主成分分析方法实现输入变量的预白化处理，能够简化神经网络模型的输入变量，进而提高模型的性能。