CN103838209A - 丙烯聚合生产过程自适应最优预报系统和方法 - Google Patents

Abstract

本发明公开了一种丙烯聚合生产过程自适应最优预报系统，包括丙烯聚合生产过程、现场智能仪表、控制站、存放数据的DCS数据库、自适应最优预报系统以及熔融指数预报值显示仪。现场智能仪表及控制站与丙烯聚合生产过程相连，与DCS数据库相连；最优预报系统与DCS数据库及预报值显示仪相连。所述的自适应最优预报系统包括模型更新模块、数据预处理模块、PCA主成分分析模块、神经网络模型模块以及自适应智能优化模块。以及提供了一种用预报系统实现的预报方法。本发明实现在线测量、在线参数优化、预报速度快、模型自动更新、抗干扰能力强、精度高。

Description

丙烯聚合生产过程自适应最优预报系统和方法

技术领域

本发明涉及一种最优预报系统及方法，具体是一种丙烯聚合生产过程自适应最优预报系统及方法。

背景技术

聚丙烯是由丙烯聚合而制得的一种热塑性树脂，丙烯最重要的下游产品，世界丙烯的50%，我国丙烯的65%都是用来制聚丙烯，是五大通用塑料之一，与我们的日常生活密切相关。聚丙烯是世界上增长最快的通用热塑性树脂，总量仅仅次于聚乙烯和聚氯乙烯。为使我国聚丙烯产品具有市场竞争力，开发刚性、韧性、流动性平衡好的抗冲共聚产品、无规共聚产品、BOPP和CPP薄膜料、纤维、无纺布料，及开发聚丙烯在汽车和家电领域的应用，都是今后重要的研究课题。

熔融指数是聚丙烯产品确定产品牌号的重要质量指标之一，它决定了产品的不同用途，对熔融指数的测量是聚丙烯生产中产品质量控制的一个重要环节，对生产和科研，都有非常重要的作用和指导意义。

然而，熔融指数的在线分析测量目前很难做到，一方面是在线熔融指数分析仪的缺乏，另一方面是现有的在线分析仪由于经常会堵塞而测量不准甚至无法正常使用所导致的使用上的困难。因此，目前工业生产中MI的测量，主要是通过人工取样、离线化验分析获得，而且一般每2-4小时只能分析一次，时间滞后大，给丙烯聚合生产的质量控制带来了困难，成为生产中急需解决的一个瓶颈问题。聚丙烯熔融指数的在线预报系统及方法研究，从而成为学术界和工业界的一个前沿和热点。

发明内容

为了克服目前已有的丙烯聚合生产过程的测量精度不高、易受人为因素的影响的不足，本发明的目的在于提供一种在线测量、在线参数优化、预报速度快、模型自动更新、抗干扰能力强、精度高的基于智能粒子群算法寻优的丙烯聚合生产过程熔融指数最优预报系统及方法。

本发明解决其技术问题所采用的技术方案是：

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

所述自适应最优预报系统包括：

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

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

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

（4）、自适应智能优化模块，用于采用基于智能粒子群算法的优化模块对神经网络进行优化，包括：

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

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

$x_{k+1}^i=x_k^i+v_{k+1}^i---\left(1\right)$

$v_{k+1}^i={\mathrm{wv}}_k^i+c_1{r}_1(p_k^i-x_k^i)+c_2{r}_2(p_k^g-x_k^i)---\left(2\right)$

式中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）这里的智能粒子群算法体现在速度权重值w是基于一种进化状态来进行估计的；粒子群的进化状态包括了探索期、开拓期、聚合器和跳出期；由f来定量表示；

$f=\frac{d_g-d_{\min }}{d_{\max }-d_{\min }}\∈\left[\mathrm{0,1}\right]---\left(3\right)$

d_i是粒子i到粒子群中其他粒子的平均距离；d_g是粒子群中最优解到粒子群中其他粒子的平均距离；d_max和d_min分别是{d_i}中的最大、最小值；d_i按照下式计算：

$d_i=\frac{1}{P-1}\underset{j=1,j\≠i}{\overset{P}{\mathrm{\Σ}}}{\left|\right|x_i-x_j\left|\right|}_2---\left(4\right)$

x_i为第i个粒子的位置，x_j为第j个粒子的位置，‖·‖为范数表达式，P为粒子群中粒子的数目；

由f值来得到粒子群速度权重的计算式为：

$w=\frac{1}{1+{1.5e}^{-2.6f}}\∈\left[\mathrm{0.4,0.9}\right]---\left(5\right)$

作为优选的一种方案，所述基于智能粒子群算法寻优的最优预报模型还包括：模型更新模块，用于模型的在线更新，将定期将离线化验数据输入到训练集中，更新神经网络模型。

作为优选的再一种方案：在所述自适应智能优化模块中，粒子群执行迭代时所用到的w是根据粒子群的进化状态来智能变化的；这样的智能设计目的在于使得粒子群能够在不同的阶段采用不同的速度权重，控制群体的收敛速度，提高整个群体的全局寻优能力。

作为优选的再一种方案：在PCA主成分分析模块中，PCA方法实现输入变量的预白化处理，能够简化神经网络模型的输入变量，进而提高模型的性能。

2、一种丙烯聚会生产过程自适应最优预报系统实现的预报方法，所述预报方法具体实现步骤如下：

（1）对丙烯聚合生产过程对象，根据工艺分析和操作分析，选择操作变量和易测变量作为模型的输入，操作变量和易测变量由DCS数据库获得；

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

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

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

（5）采用基于智能粒子群算法的优化模块对神经网络进行优化，包括：

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

（5.2）通过下式执行智能粒子群算法，让粒子群收敛：

$x_{k+1}^i=x_k^i+v_{k+1}^i---\left(1\right)$

$v_{k+1}^i={\mathrm{wv}}_k^i+c_1{r}_1(p_k^i-x_k^i)+c_2{r}_2(p_k^g-x_k^i)---\left(2\right)$

式中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₂分别为随机数。

（5.3）这里的智能粒子群算法体现在速度权重值w是基于一种进化状态来进行估计的；粒子群的进化状态包括了探索期、开拓期、聚合器和跳出期；由f来定量表示；

$f=\frac{d_g-d_{\min }}{d_{\max }-d_{\min }}\∈\left[\mathrm{0,1}\right]---\left(3\right)$

d_i是粒子i到粒子群中其他粒子的平均距离；d_g是粒子群中最优解到粒子群中其他粒子的平均距离；d_max和d_min分别是{d_i}中的最大、最小值；d_i按照下式计算：

$d_i=\frac{1}{P-1}\underset{j=1,j\≠i}{\overset{P}{\mathrm{\Σ}}}{\left|\right|x_i-x_j\left|\right|}_2---\left(4\right)$

x_i为第i个粒子的位置，x_j为第j个粒子的位置，‖·‖为范数表达式，P为粒子群中粒子的数目；

由f值来得到粒子群速度权重的计算式为：

$w=\frac{1}{1+{1.5e}^{-2.6f}}\∈\left[\mathrm{0.4,0.9}\right]---\left(5\right)$

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

本发明的技术构思为：对丙烯聚合生产过程的重要质量指标熔融指数进行在线最优预报，克服已有的聚丙烯熔融指数测量仪表测量精度不高、易受人为因素的影响的不足，引入自适应智能优化模块对神经网络参数和结构进行自动优化，不需要人为经验或多次测试来调整神经网络，从而得到具有最优的熔融指数预报功能的在线最优预报系统。

本发明的有益效果主要表现在：1、在线测量；2、在线参数自动优化；3、预报速度快；4、模型自动更新；5、抗干扰能力强；6、精度高。

附图说明

图1是丙烯聚合生产过程自适应最优预报系统及方法的基本结构示意图；

图2是自适应最优预报系统结构示意图；

图3是丙烯聚合生产过程Hypol工艺生产流程图。

具体实施方式

下面结合附图对本发明做进一步描述。本发明实施例用来解释说明本发明，而不是对本发明进行限制，在本发明的精神和权利要求的保护范围内，对本发明作出的任何修改和改变，都落入本发明的保护范围。

实施例1

1.参照图1、图2和图3，一种丙烯聚合生产过程自适应最优预报系统，包括丙烯聚合生产过程1、用于测量易测变量的现场智能仪表2、用于测量操作变量的控制站3、存放数据的DCS数据库4、基于智能蚁群算法寻优的最优预报系统5以及熔融指数预报值显示仪6，所述现场智能仪表2、控制站3与丙烯聚合生产过程1连接，所述现场智能仪表2、控制站3与DCS数据库4连接，所述DCS数据库4与自适应最优预报系统5的输入端连接，所述自适应最优预报系统5的输出端与熔融指数预报值显示仪6连接，所述自适应最优预报系统包括：

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

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

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

（4）、自适应智能优化模块，用于采用基于智能粒子群算法的优化模块对神经网络进行优化，包括：

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

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

$x_{k+1}^i=x_k^i+v_{k+1}^i---\left(1\right)$

$v_{k+1}^i={\mathrm{wv}}_k^i+c_1{r}_1(p_k^i-x_k^i)+c_2{r}_2(p_k^g-x_k^i)---\left(2\right)$

式中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）这里的智能粒子群算法体现在速度权重值w是基于一种进化状态来进行估计的；粒子群的进化状态包括了探索期、开拓期、聚合器和跳出期；由f来定量表示；

$f=\frac{d_g-d_{\min }}{d_{\max }-d_{\min }}\∈\left[\mathrm{0,1}\right]---\left(3\right)$

d_i是粒子i到粒子群中其他粒子的平均距离；d_g是粒子群中最优解到粒子群中其他粒子的平均距离；d_max和d_min分别是{d_i}中的最大、最小值；d_i按照下式计算：

$d_i=\frac{1}{P-1}\underset{j=1,j\≠i}{\overset{P}{\mathrm{\Σ}}}{\left|\right|x_i-x_j\left|\right|}_2---\left(4\right)$

x_i为第i个粒子的位置，x_j为第j个粒子的位置，‖·‖为范数表达式，P为粒子群中粒子的数目；

由f值来得到粒子群速度权重的计算式为：

$w=\frac{1}{1+{1.5e}^{-2.6f}}\∈\left[\mathrm{0.4,0.9}\right]---\left(5\right)$

在PCA主成分分析模块中，PCA方法实现输入变量的预白化处理，能够简化神经网络模型的输入变量，进而提高模型的性能。

2.丙烯聚合生产过程流程图如图3所示，根据反应机理以及流程工艺分析，考虑到聚丙烯生产过程中对熔融指数产生影响的各个因素，取实际生产过程中常用的九个操作变量和易测变量作为模型输入变量，有：三股丙烯进料流率，主催化剂流率，辅催化剂流率，釜内温度、压强、液位，釜内氢气体积浓度。

表1 自适应最优预报系统所需模型输入变量

表1列出了作为自适应最优预报系统5输入的9个模型输入变量，分别为釜内温度（T）、釜内压强（p）、釜内液位（L）、釜内氢气体积浓度（X_v）、3股丙烯进料流率（第一股丙烯进料流率f1，第二股丙烯进料流率f2，第三股丙烯进料流率f3）、2股催化剂进料流率（主催化剂流率f4，辅催化剂流率f5）。反应釜中的聚合反应是反应物料反复混合后参与反应的，因此模型输入变量涉及物料的过程变量采用前若干时刻的平均值。此例中数据采用前一小时的平均值。熔融指数离线化验值作为自适应最优预报系统5的输出变量。通过人工取样、离线化验分析获得，每4小时分析采集一次。

现场智能仪表2及控制站3与丙烯聚合生产过程1相连，与DCS数据库4相连；最优预报系统5与DCS数据库4及预报值显示仪6相连。现场智能仪表2测量丙烯聚合生产对象的易测变量，将易测变量传输到DCS数据库4；控制站3控制丙烯聚合生产对象的操作变量，将操作变量传输到DCS数据库4。DCS数据库4中记录的变量数据作为自适应最优预报系统5的输入，预报值显示仪6用于显示自适应最优预报系统5的输出，即预报值。

自适应最优预报系统5，包括：

（1）数据预处理模块7，用于对模型输入进行预处理，即中心化和归一化。对输入变量中心化，就是减去变量的平均值，使变量为零均值的变量，从而简化算法；对输入变量归一化，就是除以输入变量值的变化区间，是变量的值落到-0.5～0.5之内，进一步简化。

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

（3）神经网络模型模块9，采用RBF神经网络，多层前馈神经网络在网络结构上通常由输入层、隐含层和输出层组成。在网络特征上主要表现为既无层内神经元的互联，也无层间的反联络。这种网络实质上是一种静态网络，其输出只是现行输入的函数，而于过去和将来的输入或输出无关。RBF神经网络模型有一个输入层、一个输出层和一个隐藏层。理论上可以证明，RBF神经网络可以任意逼近非线性系统。RBF神经网络训练算法通过误差函数最小化来完成输入到输出的一种高度非线性映射，映射中保持拓扑不变性。

（4）自适应智能优化模块10：采用基于自适应智能优化模块对神经网络进行优化，通过智能粒子群算法寻优的强大全局寻优能力来优化神经网络的输入和隐藏层结构，并以此进行神经网络学习，从而建立丙烯聚合熔融指数的智能粒子群算法寻优优化的RBF神经网络最优预报模型。具体实现步骤如下：

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

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

$x_{k+1}^i=x_k^i+v_{k+1}^i---\left(1\right)$

$v_{k+1}^i={\mathrm{wv}}_k^i+c_1{r}_1(p_k^i-x_k^i)+c_2{r}_2(p_k^g-x_k^i)---\left(2\right)$

式中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）这里的智能粒子群算法体现在速度权重值w是基于一种进化状态来进行估计的；粒子群的进化状态包括了探索期、开拓期、聚合器和跳出期；由f来定量表示；

$f=\frac{d_g-d_{\min }}{d_{\max }-d_{\min }}\∈\left[\mathrm{0,1}\right]---\left(3\right)$

d_i是粒子i到粒子群中其他粒子的平均距离；d_g是粒子群中最优解到粒子群中其他粒子的平均距离；d_max和d_min分别是{d_i}中的最大、最小值；d_i按照下式计算：

$d_i=\frac{1}{P-1}\underset{j=1,j\≠i}{\overset{P}{\mathrm{\Σ}}}{\left|\right|x_i-x_j\left|\right|}_2---\left(4\right)$

x_i为第i个粒子的位置，x_j为第j个粒子的位置，‖·‖为范数表达式，P为粒子群中粒子的数目；

由f值来得到粒子群速度权重的计算式为：

$w=\frac{1}{1+{1.5e}^{-2.6f}}\∈\left[\mathrm{0.4,0.9}\right]---\left(5\right)$

（5）模型更新模块11，用于模型的在线更新，定期将离线化验数据输入到训练集中，更新神经网络模型。

实施例2

1.参照图1、图2和图3，一种基于智能粒子群算法寻优的丙烯聚合生产过程最优预报方法包括以下步骤：

（1）对丙烯聚合生产过程对象，根据工艺分析和操作分析，选择操作变量和易测变量作为模型的输入，操作变量和易测变量由DCS数据库获得；

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

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

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

（5）采用基于智能粒子群算法的优化模块对神经网络进行优化，包括：

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

（5.2）通过下式执行智能粒子群算法，让粒子群收敛：

$x_{k+1}^i=x_k^i+v_{k+1}^i---\left(1\right)$

$v_{k+1}^i={\mathrm{wv}}_k^i+c_1{r}_1(p_k^i-x_k^i)+c_2{r}_2(p_k^g-x_k^i)---\left(2\right)$

式中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₂分别为随机数。

（5.3）这里的智能粒子群算法体现在速度权重值w是基于一种进化状态来进行估计的；粒子群的进化状态包括了探索期、开拓期、聚合器和跳出期；由f来定量表示；

$f=\frac{d_g-d_{\min }}{d_{\max }-d_{\min }}\∈\left[\mathrm{0,1}\right]---\left(3\right)$

d_i是粒子i到粒子群中其他粒子的平均距离；d_g是粒子群中最优解到粒子群中其他粒子的平均距离；d_max和d_min分别是{d_i}中的最大、最小值；d_i按照下式计算：

$d_i=\frac{1}{P-1}\underset{j=1,j\≠i}{\overset{P}{\mathrm{\Σ}}}{\left|\right|x_i-x_j\left|\right|}_2---\left(4\right)$

x_i为第i个粒子的位置，x_j为第j个粒子的位置，‖·‖为范数表达式，P为粒子群中粒子的数目；

由f值来得到粒子群速度权重的计算式为：

$w=\frac{1}{1+{1.5e}^{-2.6f}}\∈\left[\mathrm{0.4,0.9}\right]---\left(5\right)$

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

2.本实施例的方法具体实施步骤如下：

步骤1：对丙烯聚合生产过程对象1，根据工艺分析和操作分析，选择操作变量和易测变量作为模型的输入。

步骤2：对样本数据进行预处理，由数据预处理模块7完成。

步骤3：对经过预处理的数据进行主成分分析，由PCA主成分分析模块8完成。

步骤4：基于模型输入、输出建立初始神经网络模型9。输入数据如步骤1所述获得，输出数据由离线化验获得。

步骤5：由自适应智能优化模块10优化初始神经网络9的输入和隐藏层结构。

步骤6：模型更新模块11定期将离线化验数据输入到训练集中，更新神经网络模型，自适应最优预报系统5建立完成。

步骤7：建立好的自适应最优预报系统5基于DCS数据库4传来的实时模型输入变量数据对丙烯聚合生产过程1的熔融指数进行基于智能粒子群算法寻优的最优预报。

步骤8：熔融指数预报显示仪6显示自适应最优预报系统5的输出，完成对丙烯聚合生产过程熔融指数的最优预报的显示。

Claims (2)

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）、自适应智能优化模块，用于采用基于智能粒子群算法的优化模块对神经网络进行优化，包括：

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

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

$x_{k+1}^i=x_k^i+v_{k+1}^i---\left(1\right)$

$v_{k+1}^i={\mathrm{wv}}_k^i+c_1{r}_1(p_k^i-x_k^i)+c_2{r}_2(p_k^g-x_k^i)---\left(2\right)$

式中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）这里的智能粒子群算法体现在速度权重值w是基于一种进化状态来进行估计的；粒子群的进化状态包括了探索期、开拓期、聚合器和跳出期；由f来定量表示；

$f=\frac{d_g-d_{\min }}{d_{\max }-d_{\min }}\∈\left[\mathrm{0,1}\right]---\left(3\right)$

d_i是粒子i到粒子群中其他粒子的平均距离；d_g是粒子群中最优解到粒子群中其他粒子的平均距离；d_max和d_min分别是{d_i}中的最大、最小值；d_i按照下式计算：

$d_i=\frac{1}{P-1}\underset{j=1,j\≠i}{\overset{P}{\mathrm{\Σ}}}{\left|\right|x_i-x_j\left|\right|}_2---\left(4\right)$

x_i为第i个粒子的位置，x_j为第j个粒子的位置，‖·‖为范数表达式，P为粒子群中粒子的数目；

由f值来得到粒子群速度权重的计算式为：

$w=\frac{1}{1+{1.5e}^{-2.6f}}\∈\left[\mathrm{0.4,0.9}\right]---\left(5\right)$

所述自适应最优预报系统还包括：模型更新模块，用于模型的在线更新，定期将离线化验数据输入到训练集中，更新神经网络模型。

所述的自适应最优预报系统，其特征在于：在所述自适应智能优化模块中，粒子群执行迭代时所用到的w是根据粒子群的进化状态来智能变化的；这样的智能设计目的在于使得粒子群能够在不同的阶段采用不同的速度权重，控制群体的收敛速度，提高整个群体的全局寻优能力。PCA主成分分析模块中，PCA方法实现输入变量的预白化处理，能够简化神经网络模型的输入变量，进而提高模型的性能。

2.一种用如权利要求1所述的丙烯聚会生产过程自适应最优预报系统实现的预报方法，其特征在于：所述预报方法具体实现步骤如下：

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

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

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

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

(5.5)采用基于智能粒子群算法的优化模块对神经网络进行优化，包括：

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

（5.5.2）通过下式执行智能粒子群算法，让粒子群收敛：

$x_{k+1}^i=x_k^i+v_{k+1}^i---\left(1\right)$

$v_{k+1}^i={\mathrm{wv}}_k^i+c_1{r}_1(p_k^i-x_k^i)+c_2{r}_2(p_k^g-x_k^i)---\left(2\right)$

式中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₂分别为随机数。

（5.5.3）这里的智能粒子群算法体现在速度权重值w是基于一种进化状态来进行估计的；粒子群的进化状态包括了探索期、开拓期、聚合器和跳出期；由f来定量表示；

$f=\frac{d_g-d_{\min }}{d_{\max }-d_{\min }}\∈\left[\mathrm{0,1}\right]---\left(3\right)$

d_i是粒子i到粒子群中其他粒子的平均距离；d_g是粒子群中最优解到粒子群中其他粒子的平均距离；d_max和d_min分别是{d_i}中的最大、最小值；d_i按照下式计算：

$d_i=\frac{1}{P-1}\underset{j=1,j\≠i}{\overset{P}{\mathrm{\Σ}}}{\left|\right|x_i-x_j\left|\right|}_2---\left(4\right)$

x_i为第i个粒子的位置，x_j为第j个粒子的位置，‖·‖为范数表达式，P为粒子群中粒子的数目；

由f值来得到粒子群速度权重的计算式为：

$w=\frac{1}{1+{1.5e}^{-2.6f}}\∈\left[\mathrm{0.4,0.9}\right]---\left(5\right)$

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

所述的预报方法，其特征在于：在所述自适应智能优化模块中，粒子群执行迭代时所用到的w是根据粒子群的进化状态来智能变化的；这样的智能设计目的在于使得粒子群能够在不同的阶段采用不同的速度权重，控制群体的收敛速度，提高整个群体的全局寻优能力。在所述的步骤（5.3）中采用PCA主成分分析方法实现输入变量的预白化处理，能够简化神经网络模型的输入变量，进而提高模型的性能。

Images