CN103823964A - 一种全局最优丙烯聚合生产过程最优软测量仪表和方法 - Google Patents

Abstract

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

Description

一种全局最优丙烯聚合生产过程最优软测量仪表和方法

技术领域

本发明涉及一种最优软测量仪表及方法，具体是一种一种全局最优丙烯聚合生产过程最优软测量仪表及方法。

背景技术

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

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

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

发明内容

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

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

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

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

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

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

（4）、全局优化模块，用于采用全局优化模块对神经网络进行优化，包括：

（4.1）算法初始化，根据待优化的RBF神经网络结构构造出初始的解集S=(s₁,s₂,…,s_n)，n为初始解的个数，sn为第n个初始解,确定蚁群的大小m，设置蚁群寻优算法迭代次数的阈值MaxGen并初始化蚁群寻优的迭代次数序号gen=0；

（4.2）计算出解集S对应的适应度值F_i(i=1,2,…,n)，适应度值越大代表解越好；再根据下式确定解集中每个解被取到作为蚂蚁寻优的初始解的概率P_i(i=1,2,…,n)

$P_i\left(k\right)=\frac{F_i}{{\mathrm{\Σ}}_{i=1}^n{F}_i}(i=\mathrm{1,2},\·\·\·,n)---\left(1\right)$

n为初始解的个数，sn为第n个初始解，k为迭代次数。初始化执行寻优算法的蚂蚁编号a=0；

（4.3）蚂蚁a选取S中的一个解作为寻优的初始解，选取规则是根据P来做轮盘选；

（4.4）蚂蚁a在选取的初始解的基础上进行寻优，找到更好的解s_a′；

（4.5）如果a<m，则a=a+1，返回步骤4.3；否则继续向下执行步骤4.6；

（4.6）如果gen<MaxGen，则gen=gen+1，使用步骤4.4中所有蚂蚁得到的更好的解取代S中的对应解，返回步骤4.2；否则向下执行步骤4.7；

（4.7）计算出解集S对应的适应度值F_a(a=1,2,…,n)，选取适应度值最大的解作为算法的最优解，结束算法并返回。

作为优选的一种方案，所述全局最优软测量仪还包括：模型更新模块，用于模型的在线更新，将定期将离线化验数据输入到训练集中，更新神经网络模型。

作为优选的再一种方案：在所述的全局优化模块中，每一只蚂蚁在它选定的初始解的基础上寻优时会循环固定的次数，如果本次循环得到了更好的解，则在下次循环中会基于该解并保持搜索方向不变；否则在下次循环中仍基于原来的解但会调整搜索方向；

同时随着整个蚁群寻优代数的增加，蚂蚁搜索的步长会智能减小，以适合整个蚁群寻优的收敛：

del_k=Random·k^r     (2)

式中，del_k为蚂蚁第k代迭代的初始步长，k为迭代代数，Random为随机向量，r为负常实数。

对于解集S中长期不被蚂蚁选作寻优初始解的那些解，会采用遗传算法中的变异和交叉策略进行处理，避免群体早熟，提高智能连续空间蚁群算法算法的全局寻优性能。

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

一种全局最优聚丙烯生产过程最优软测量仪表实现的软测量方法，所述软测量方法具体实现步骤如下：

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

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

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

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

（5）全局优化模块，用于采用全局优化模块对神经网络进行优化，包括：

（5.5.1）算法初始化，根据待优化的RBF神经网络结构构造出初始的解集S=(s₁,s₂,…,s_n)，n为初始解的个数，sn为第n个初始解,确定蚁群的大小m，设置蚁群寻优算法迭代次数的阈值MaxGen并初始化蚁群寻优的迭代次数序号gen=0；

（5.5.2）计算出解集S对应的适应度值F_i(i=1,2,…,n)，适应度值越大代表解越好；再根据下式确定解集中每个解被取到作为蚂蚁寻优的初始解的概率P_i(i=1,2,…,n)

$P_i\left(k\right)=\frac{F_i}{{\mathrm{\&Sigma;}}_{i=1}^n{F}_i}(i=\mathrm{1,2},\&CenterDot;\&CenterDot;\&CenterDot;,n)---\left(1\right)$

n为初始解的个数，sn为第n个初始解，k为迭代次数。初始化执行寻优算法的蚂蚁编号a=0；

（5.5.3）蚂蚁a选取S中的一个解作为寻优的初始解，选取规则是根据P来做轮盘选；

（5.5.4）蚂蚁a在选取的初始解的基础上进行寻优，找到更好的解s_a′；

（5.5.5）如果a<m，则a=a+1，返回步骤4.3；否则继续向下执行步骤4.6；

（5.5.6）如果gen<MaxGen，则gen=gen+1，使用步骤4.4中所有蚂蚁得到的更好的解取代S中的对应解，返回步骤4.2；否则向下执行步骤4.7；

（5.5.7）计算出解集S对应的适应度值F_a(a=1,2,…,n)，选取适应度值最大的解作为算法的最优解，结束算法并返回。

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

作为优选的再一种方案：在所述的智能连续空间蚁群算法寻优步骤（5.5）中，每一只蚂蚁在它选定的初始解的基础上寻优时会循环固定的次数，如果本次循环得到了更好的解，则在下次循环中会基于该解并保持搜索方向不变；否则在下次循环中仍基于原来的解但会调整搜索方向；

同时随着整个蚁群寻优代数的增加，蚂蚁搜索的步长会智能减小，以适合整个蚁群寻优的收敛：

del_k=Random·k^r     (2)

式中，del_k为蚂蚁第k代迭代的初始步长，k为迭代代数，Random为随机向量，r为负常实数。

对于解集S中长期不被蚂蚁选作寻优初始解的那些解，会采用遗传算法中的变异和交叉策略进行处理，从而提高智能连续空间蚁群算法算法的全局寻优性能。

进一步，在所述的步骤（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=xU得到，其中x为输入变量，C为主成分得分矩阵，U为载荷矩阵。若对原始数据进行重构，可由x=CU^T计算，其中上标T表示矩阵的转置。当选取的主成分数目小于输入变量的变量个数时，x=CU^T+E，其中E为残差矩阵；

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

（4）、全局优化模块，用于采用全局优化模块对神经网络进行优化，包括：

（4.1）算法初始化，根据待优化的RBF神经网络结构构造出初始的解集S=(s₁,s₂,…,s_n)，n为初始解的个数，sn为第n个初始解,确定蚁群的大小m，设置蚁群寻优算法迭代次数的阈值MaxGen并初始化蚁群寻优的迭代次数序号gen=0；

（4.2）计算出解集S对应的适应度值F_i(i=1,2,…,n)，适应度值越大代表解越好；再根据下式确定解集中每个解被取到作为蚂蚁寻优的初始解的概率P_i(i=1,2,…,n)

$P_i\left(k\right)=\frac{F_i}{{\mathrm{\&Sigma;}}_{i=1}^n{F}_i}(i=\mathrm{1,2},\&CenterDot;\&CenterDot;\&CenterDot;,n)---\left(1\right)$

n为初始解的个数，sn为第n个初始解，k为迭代次数。初始化执行寻优算法的蚂蚁编号a=0；

（4.3）蚂蚁a选取S中的一个解作为寻优的初始解，选取规则是根据P来做轮盘选；

（4.4）蚂蚁a在选取的初始解的基础上进行寻优，找到更好的解s_a′；

（4.5）如果a<m，则a=a+1，返回步骤4.3；否则继续向下执行步骤4.6；

（4.6）如果gen<MaxGen，则gen=gen+1，使用步骤4.4中所有蚂蚁得到的更好的解取代S中的对应解，返回步骤4.2；否则向下执行步骤4.7；

（4.7）计算出解集S对应的适应度值F_a(a=1,2,…,n)，选取适应度值最大的解作为算法的最优解，结束算法并返回。

作为优选的一种方案，所述全局最优软测量仪还包括：模型更新模块，用于模型的在线更新，将定期将离线化验数据输入到训练集中，更新神经网络模型。

作为优选的再一种方案：在所述的全局优化模块中，每一只蚂蚁在它选定的初始解的基础上寻优时会循环固定的次数，如果本次循环得到了更好的解，则在下次循环中会基于该解并保持搜索方向不变；否则在下次循环中仍基于原来的解但会调整搜索方向；

同时随着整个蚁群寻优代数的增加，蚂蚁搜索的步长会智能减小，以适合整个蚁群寻优的收敛：

del_k=Random·k^r     (2)

式中，del_k为蚂蚁第k代迭代的初始步长，k为迭代代数，Random为随机向量，r为负常实数。

对于解集S中长期不被蚂蚁选作寻优初始解的那些解，会采用遗传算法中的变异和交叉策略进行处理，避免群体早熟，提高智能连续空间蚁群算法算法的全局寻优性能。

作为优选的再一种方案：在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=xU得到，其中x为输入变量，C为主成分得分矩阵，U为载荷矩阵。若对原始数据进行重构，可由x=CU^T计算，其中上标T表示矩阵的转置。当选取的主成分数目小于输入变量的变量个数时，x=CU^T+E，其中E为残差矩阵。

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

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

（4.1）算法初始化，根据待优化的RBF神经网络结构构造出初始的解集S=(s₁,s₂,…,s_n)，n为初始解的个数，sn为第n个初始解,确定蚁群的大小m，设置蚁群寻优算法迭代次数的阈值MaxGen并初始化蚁群寻优的迭代次数序号gen=0；

（4.2）计算出解集S对应的适应度值F_i(i=1,2,…,n)，适应度值越大代表解越好；再根据下式确定解集中每个解被取到作为蚂蚁寻优的初始解的概率P_i(i=1,2,…,n)

$P_i\left(k\right)=\frac{F_i}{{\mathrm{\&Sigma;}}_{i=1}^n{F}_i}(i=\mathrm{1,2},\&CenterDot;\&CenterDot;\&CenterDot;,n)---\left(1\right)$

n为初始解的个数，sn为第n个初始解，k为迭代次数。初始化执行寻优算法的蚂蚁编号a=0；

（4.3）蚂蚁a选取S中的一个解作为寻优的初始解，选取规则是根据P来做轮盘选；

（4.4）蚂蚁a在选取的初始解的基础上进行寻优，找到更好的解s_a′；

（4.5）如果a<m，则a=a+1，返回步骤4.3；否则继续向下执行步骤4.6；

（4.6）如果gen<MaxGen，则gen=gen+1，使用步骤4.4中所有蚂蚁得到的更好的解取代S中的对应解，返回步骤4.2；否则向下执行步骤4.7；

（4.7）计算出解集S对应的适应度值F_a(a=1,2,…,n)，选取适应度值最大的解作为算法的最优解，结束算法并返回。

在所述的全局优化模块中，每一只蚂蚁在它选定的初始解的基础上寻优时会循环固定的次数，如果本次循环得到了更好的解，则在下次循环中会基于该解并保持搜索方向不变；否则在下次循环中仍基于原来的解但会调整搜索方向；

同时随着整个蚁群寻优代数的增加，蚂蚁搜索的步长会智能减小，以适合整个蚁群寻优的收敛：

del_k=Random·k^r     (2)

式中，del_k为蚂蚁第k代迭代的初始步长，k为迭代代数，Random为随机向量，r为负常实数。

对于解集S中长期不被蚂蚁选作寻优初始解的那些解，会采用遗传算法中的变异和交叉策略进行处理，避免群体早熟，提高智能连续空间蚁群算法算法的全局寻优性能。

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

实施例2

1.参照图1、图2和图3，一种全局最优丙烯聚合生产过程最优软测量方法包括以下步骤：

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

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

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

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

（5）全局优化模块，用于采用全局优化模块对神经网络进行优化，包括：

（5.5.1）算法初始化，根据待优化的RBF神经网络结构构造出初始的解集S=(s₁,s₂,…,s_n)，n为初始解的个数，sn为第n个初始解,确定蚁群的大小m，设置蚁群寻优算法迭代次数的阈值MaxGen并初始化蚁群寻优的迭代次数序号gen=0；

（5.5.2）计算出解集S对应的适应度值F_i(i=1,2,…,n)，适应度值越大代表解越好；再根据下式确定解集中每个解被取到作为蚂蚁寻优的初始解的概率P_i(i=1,2,…,n)

$P_i\left(k\right)=\frac{F_i}{{\mathrm{\&Sigma;}}_{i=1}^n{F}_i}(i=\mathrm{1,2},\&CenterDot;\&CenterDot;\&CenterDot;,n)---\left(1\right)$

n为初始解的个数，sn为第n个初始解，k为迭代次数。初始化执行寻优算法的蚂蚁编号a=0；

（5.5.3）蚂蚁a选取S中的一个解作为寻优的初始解，选取规则是根据P来做轮盘选；

（5.5.4）蚂蚁a在选取的初始解的基础上进行寻优，找到更好的解s_a′；

（5.5.5）如果a<m，则a=a+1，返回步骤4.3；否则继续向下执行步骤4.6；

（5.5.6）如果gen<MaxGen，则gen=gen+1，使用步骤4.4中所有蚂蚁得到的更好的解取代S中的对应解，返回步骤4.2；否则向下执行步骤4.7；

（5.5.7）计算出解集S对应的适应度值F_a(a=1,2,…,n)，选取适应度值最大的解作为算法的最优解，结束算法并返回。

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

在所述的智能连续空间蚁群算法寻优步骤（5.5）中，每一只蚂蚁在它选定的初始解的基础上寻优时会循环固定的次数，如果本次循环得到了更好的解，则在下次循环中会基于该解并保持搜索方向不变；否则在下次循环中仍基于原来的解但会调整搜索方向；

同时随着整个蚁群寻优代数的增加，蚂蚁搜索的步长会智能减小，以适合整个蚁群寻优的收敛：

del_k=Random·k^r     (2)

式中，del_k为蚂蚁第k代迭代的初始步长，k为迭代代数，Random为随机向量，r为负常实数。

对于解集S中长期不被蚂蚁选作寻优初始解的那些解，会采用遗传算法中的变异和交叉策略进行处理，从而提高智能连续空间蚁群算法算法的全局寻优性能。

进一步，在所述的步骤（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=xU得到，其中x为输入变量，C为主成分得分矩阵，U为载荷矩阵。若对原始数据进行重构，可由x=CU^T计算，其中上标T表示矩阵的转置。当选取的主成分数目小于输入变量的变量个数时，x=CU^T+E，其中E为残差矩阵；

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

（4）、全局优化模块，用于采用全局优化模块对神经网络进行优化，包括：

（4.1）算法初始化，根据待优化的RBF神经网络结构构造出初始的解集S=(s₁,s₂,…,s_n)，n为初始解的个数，sn为第n个初始解,确定蚁群的大小m，设置蚁群寻优算法迭代次数的阈值MaxGen并初始化蚁群寻优的迭代次数序号gen=0；

（4.2）计算出解集S对应的适应度值F_i(i=1,2,…,n)，适应度值越大代表解越好；再根据下式确定解集中每个解被取到作为蚂蚁寻优的初始解的概率P_i(i=1,2,…,n)

$P_i\left(k\right)=\frac{F_i}{{\mathrm{\&Sigma;}}_{i=1}^n{F}_i}(i=\mathrm{1,2},\&CenterDot;\&CenterDot;\&CenterDot;,n)---\left(1\right)$

n为初始解的个数，sn为第n个初始解，k为迭代次数。初始化执行寻优算法的蚂蚁编号a=0；

（4.3）蚂蚁a选取S中的一个解作为寻优的初始解，选取规则是根据P来做轮盘选；

（4.4）蚂蚁a在选取的初始解的基础上进行寻优，找到更好的解s_a′；

（4.5）如果a<m，则a=a+1，返回步骤4.3；否则继续向下执行步骤4.6；

（4.6）如果gen<MaxGen，则gen=gen+1，使用步骤4.4中所有蚂蚁得到的更好的解取代S中的对应解，返回步骤4.2；否则向下执行步骤4.7；

（4.7）计算出解集S对应的适应度值F_a(a=1,2,…,n)，选取适应度值最大的解作为算法的最优解，结束算法并返回。

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

在所述全局优化模块中，每一只蚂蚁在它选定的初始解的基础上寻优时会循环固定的次数，如果本次循环得到了更好的解，则在下次循环中会基于该解并保持搜索方向不变；否则在下次循环中仍基于原来的解但会调整搜索方向；同时随着整个蚁群寻优代数的增加，蚂蚁搜索的步长会智能减小，以适合整个蚁群寻优的收敛：

del_k=Random·k^r     (2)

式中，del_k为蚂蚁第k代迭代的初始步长，k为迭代代数，Random为随机向量，r为负常实数。

对于解集S中长期不被蚂蚁选作寻优初始解的那些解，会采用遗传算法中的变异和交叉策略进行处理，避免群体早熟，提高智能连续空间蚁群算法算法的全局寻优性能。

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

2.一种用如权利要求1所述的基于智能连续空间蚁群算法寻优的聚丙烯生产过程最优软测量仪表实现的软测量方法，其特征在于，所述软测量方法具体实现步骤如下：

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

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

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

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

(5.5)采用全局优化模块对神经网络进行优化，包括：

（5.5.1）算法初始化，根据待优化的RBF神经网络结构构造出初始的解集S=(s₁,s₂,…,s_n)，n为初始解的个数，sn为第n个初始解,确定蚁群的大小m，设置蚁群寻优算法迭代次数的阈值MaxGen并初始化蚁群寻优的迭代次数序号gen=0；

（5.5.2）计算出解集S对应的适应度值F_i(i=1,2,…,n)，适应度值越大代表解越好；再根据下式确定解集中每个解被取到作为蚂蚁寻优的初始解的概率P_i(i=1,2,…,n)

$P_i\left(k\right)=\frac{F_i}{{\mathrm{\&Sigma;}}_{i=1}^n{F}_i}(i=\mathrm{1,2},\&CenterDot;\&CenterDot;\&CenterDot;,n)---\left(1\right)$

n为初始解的个数，sn为第n个初始解，k为迭代次数。初始化执行寻优算法的蚂蚁编号a=0；

（5.5.3）蚂蚁a选取S中的一个解作为寻优的初始解，选取规则是根据P来做轮盘选；

（5.5.4）蚂蚁a在选取的初始解的基础上进行寻优，找到更好的解s_a′；

（5.5.5）如果a<m，则a=a+1，返回步骤4.3；否则继续向下执行步骤4.6；

（5.5.6）如果gen<MaxGen，则gen=gen+1，使用步骤4.4中所有蚂蚁得到的更好的解取代S中的对应解，返回步骤4.2；否则向下执行步骤4.7；

（5.5.7）计算出解集S对应的适应度值F_a(a=1,2,…,n)，选取适应度值最大的解作为算法的最优解，结束算法并返回。

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

在所述的智能连续空间蚁群算法寻优步骤（5.5）中，每一只蚂蚁在它选定的初始解的基础上寻优时会循环固定的次数，如果本次循环得到了更好的解，则在下次循环中会基于该解并保持搜索方向不变；否则在下次循环中仍基于原来的解但会调整搜索方向；

同时随着整个蚁群寻优代数的增加，蚂蚁搜索的步长会智能减小，以适合整个蚁群寻优的收敛：

del_k=Random·k^r     (2)

式中，del_k为蚂蚁第k代迭代的初始步长，k为迭代代数，Random为随机向量，r为负常实数。

对于解集S中长期不被蚂蚁选作寻优初始解的那些解，会采用遗传算法中的变异和交叉策略进行处理，从而提高智能连续空间蚁群算法算法的全局寻优性能。

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

Images