第63卷第7期 化工 学报 Vo1.63 NO.7 2012年7月 CIESC Journal July 2O12 基于D—FNN的聚合过程转化速率软测量 建模及重构 王介生,郭秋平 (辽宁科技大学电子与信息工程学院,辽宁鞍山114044) 关键词:聚合过程;动态模糊神经网络;核主元分析;软测量;模型迁移 DOI:10.3969/j.issn.0438—1157.2012.07.024 中图分类号:TP 273 文献标志码:A 文章编号:O438—1157(2012)07—2163~O7 Soft-sensor modeling and recOnfigurati0n of conversion velocity in PVC polymerization process based on D—FNN WANG Jiesheng。GUO Qiuping (School of Electronic and Information Engineering,University of Science and Technology Liaoning, Anshan 114044,Liaoning,China) Abstract:For forecasting the key technology indicator conversion velocity of vinyl chloride monomer (VCM) in the polyvinylchloride(PVC)polymerization process,a soft—sensor modeling method based on dynamic fuzzy neural network(D—FNN)was proposed.Firstly,kernel principal component analysis (KPCA)method was adopted to select the auxiliary variables of soft—sensing model in order to reduce the model dimensionality.Then the learning algorithm of D-FNN included the rule extraction principles,the classification learning strategy,the precedent parameters arrangements,the rule trimming technology based on error descendent ratio and the consequent parameters decision based on extended Kalman filter (EKF).The proposed soft—sensor model could automatically decide the fuzzy rules SO as to realize the nonlinear mapping between input and output variables of the discussed soft—sensor mode1.Model migration method was adopted to realize the on—line adaptive revision and reconfiguration of soft—sensor mode1.In the end,simulation results showed that the proposed model could significantly enhance the prediction accuracy and robustness of the technico。。economic indices and satisfy the real・-time control requirements of PVC polymerization process. Key words:polymerize process;dynamic fuzzy neural network;kernel principal component analysis; soft—sensor;model migration 士 丘 聚合工艺生产聚氯乙烯(PVC)树脂是一种典型 的间歇式化工生产过程。VCM的转化率对PVC 树脂产品质量有很大影响,不同转化率时对PVC 以氯乙烯单体(VCM)为原料,采用悬浮法 2011—12一O5收到初稿,2011—04—05收到修改稿。 联系人及第一作者:王介生(1977一),男,博士,副教授。 基金项目:中国博士后科学基金面上项目(20110491510); 辽宁省优秀人才计划项目(LJQ2011027)。 Received date:2011—12—05. Corresponding author:WANG Jiesheng,wangjiesheng@126. com Foundation item:supported by the China Postdoctoral Science Foundation(20110491510)and the Program for Liaoning Excellent Talents in University(LJQ2011027). 化工学报 第63卷 树脂的分子量、热稳定性、孔隙率、VCM的残留 软测量模型的动态调整难题。本文则提出一种数据 驱动的模型迁移方法,对所建立的聚氯乙烯生产过 程VCM转化速率DFNN软测量模型进行动态重 构,仿真结果验证了所提方法的有效性。 量、增塑剂的吸收率、加工流动性等都有一定的影 响。由于现场条件限制和缺乏成熟的检测装置,在 实际生产过程中氯乙烯的转化速率很难实时获取, 难以实现直接质量闭环控制_1]。 动态模糊神经网络(dynamic fuzzy neural net— work,D—FNN)是模糊理论和神经网络相互混合 1 聚氯乙烯(PVC)聚合过程 在聚氯乙烯(PVC)聚合过程,将各种原料 与助剂加入到反应釜内,在搅拌的作用下,充分均 匀分散,然后加入适量的引发剂开始反应,并不断 的一种神经网络模型l_2 ]。文献Es]提出了一种改 进的模糊神经网络软测量建模方法,采用LM算 法和一阶梯度下降法相结合来训练网络参数,最后 用该建模方法建立了聚合反应中熔融指数的软测量 模型,结果表明该模型具有对初始值的选择不敏 感、收敛性快、预测精度高等优点,很适合工程应 用。文献[6]采用粒子群算法优化模糊神经网络 的参数,建立丙烯腈收率软测量模型。针对软测量 模型的动态校正问题,文献I-7]利用冲激响应模 型给出了一个基于最小二乘方法的动态软测量建模 方法;文献[8]则提出一种基于在线聚类和V一支 地向反应釜的夹套和挡板通入冷却水,达到移出反 应热的目的。当氯乙烯(VCM)的转化率达到一 定值时,出现一个适当的压降,即终止反应出料。 反应完成后的浆料经汽提脱析出内含VCM后送到 干燥工序脱水干燥。典型PVC聚合釜工艺流程如 图1所示_】]。 根据聚合反应工艺特点,确定了与VCM转化 率和转化速率相关的l0个过程变量作为软测量建 模的辅助变量,分别是釜内温度(TIC—PIO1)、釜 内压力(PIC—P102)、挡板水流量(FIC—PIOI)、夹 持向量回归机的多模型软测量建模方法,试图解决 FV-PX05 VS _PX19 图1聚合釜工艺流程 Fig.1 Technique flowchart of polymerization apparatus 第7期 王介生等:基于D-FNN的聚合过程转化速率软测量建模及重构 ‘2165。 套水流量(FIC—P102)、注入水流量(FIC-P104)、 密封水流量(FI&P105)、冷却水进El温度(为夹套 水和挡板水共用,TI-PlO7)、夹套水出口温度(TI- P109)、挡板水出口温度(TI-PllO)、注入水和密封 水进口温度(即冷水槽出口温度,TIC-WA01)。 2 动态模糊神经网络 2.1 网络结构 动态模糊神经网络的结构基于扩展径向基神经 网络,不考虑反馈网络,D-FNN的结构 如图2 所示。 M} Rw N 图2动态模糊神经网络结构 Fig.2 Structure of dynamical fuzzy neural networks D-FNN本质是一个基于TSK模型的模糊系 统_9 。动态模糊神经网络共5层,具体描述 如下。 (1)输人层 , ,…,z 代表输入变量, r是输入变量个数。 (2)隶属函数层MF 代表第i个输入变量的 第J个隶属函数,本文采用高斯函数作为动态模糊 神经网络的隶属函数 (xi)一 p『L 一 0』 ].J (1) 式中 是 的第J个隶属函数,U是系统总的规 则数,c 是z 的第J个隶属函数的中心, 是 的 第J个隶属函数的宽度。 { (3)T范数层RJ代表第J条模糊规则,该层 节点代表模糊规则中的IF部分,第J个规则R 的 输出为 , 唧[一 ]=e)【p[一 ](2) 式中 一(Xl, 2,…,z )∈R ,cJ一 一 L ( 一1,2,…,“) (3) … 一—c — +.. (一 一)] ∑e (一 ) ・2 66・ 化工学报 第63卷 式中 …是预先设定的最大误差,e i 是最小误 差,也是期望的D—FNN精度, (0<fl<1)是收敛 常数,d…是输入空间的最大长度,d 是最小长 度,),(0<y<1)是衰减常数。 2.2.3 前提参数分配新模糊规则产生后,初始 参数分配方式为 式中k(是>1)是重叠因子。 需要注意的是:当第一个观测数据 (X ,t )得到后,此时D—FNN还没有建立起来, 因此这个数据被选为第一条模糊规则:C1一x , 一-O"。,其中 。是预先设定的常数。还需要指出的 是只有当IIe II>忌 ,d…>忌 时,才需增加一条模 糊规则。 2.2.4 结果参数确定 假设观测数据为 个,产 生的模糊规则数为 ,根据规则的产生准则和式 (6)可以得到N节点的输出为 一 :] , 于是,对于任意的输人x,(z z。 ,…, z,,),由式(4)计算系统的输出y 。 模型为TSK模型时,w和If,的计算公式为 11 … 1n ; 丸1 … ‰ l1‘ 11 … ≠l ・-z1 ; , 。、 丸1・zl1 … 丸 ・z1 ; 11・ ,1 … pl ‘zm j l・上 l … ・zm 线性最小二乘法(1inear least square,LLS) 和卡尔曼滤波(Kalman filter,KF)方法可以用 于结果参数调整。本文采用扩展卡尔曼滤波(ex— tended Kalman filter,EKF)方法来调整网络参 数。扩展的卡尔曼滤波是基于梯度的在线学习算 法口 ,它的基本思想是可以把全局算法划分为线 性和非线性部分,图3为EKF用于D-FNN参数 调整的学习流程,KF算法调整结果参数,EKF调 图3 EKF学习流程图 Fig.3 Learning flow chat of EKF 整前提参数。 用EKF算法对于高斯宽度优化 一 (19) S 一 H一 F sH( 一1,2,…, ) (2o) 五一 + (r7一 W )(21) 初始条件为 >O、S。一 ,其中,p)0并且 数值不大, 是第i次观测的增益矩阵,Js 是第i 次观测误差协方差矩阵, 一( , 。,…, )为 经过i次迭代后的高斯宽度向量,F 一( ,…, 3zi,…, )是第i次观测宽度梯度向量。 2.2.5修剪技术对于一个动态时变非线性系统 来说,本文采用误差下降率(ERR)方法作为模 糊规则的剔除策略以使D—FNN结构更为紧凑。如 果 反映第i个规则的重要性, 值越大说明第i 个规则越重要。若式(22)成立,则第i个规则就 可以剔除。 协<k (22) 式中k 是预先设定的阈值。 3 D—FNN软测量建模及重构 3.1 基于核主元分析法的数据降维 对于神经网络软测量模型来说,输入向量维数 过大,会使得网络拓扑结构很大和训练变得很复 杂,所以本文采用核主元分析方法(KPCA)对高 维输入向量进行降维。KPCA是将核函数的概念 引入PCA(主元分析)的一种非线性推广,它比 PCA具有更好的处理非线性问题的能力 。对 由软测量模型的输入组成的历史数据进行核主元分 析,分析结果如表1所示,前5个主元的方差贡献 率就已经达到95 以上,因此原始数据经KPCA 处理后的5个主元作为神经网络模型的输入,既 保留了原始变量的特征信息,又简化了神经网络 结构。 第7期 王介生等:基于D-FNN的聚合过程转化速率软测量建模及重构 表1不同主元个数的累计贡献率 Table 1 Cumulative contribution rate of different principle number Percent Cumulative contribution Principl “ mb Var.1ance/, 0 ra1;e/…Z0 3.2 D—FNN软测量建模 本文采用图2的动态模糊神经网络结构,建立 聚氯乙烯过程VCM转化速率的软测量模型。由前 述可知,D-FNN的5层结构分别为输入层、隶属 函数层、T范数层、归一化层和输出层。 输入层节点个数由系统的输入数决定,有几个 输入数据就产生几个节点。本文是以KPCA处理 后的5个主元变量作为输入,故输入层有5个节 点。T范数层的每一个节点与一条模糊规则的IF 部分相对应,因此节点数多少由模糊规则自动生成 策略在D-FNN训练过程中动态调整。归一化层节 点个数与T范数层节点个数相同,归一化层负责 将第二层模糊化后的数据进行去模糊化处理。输出 层显示网络输出,由具体问题决定,本文以VCM 转化速率作为输出,故输出层节点个数为1。确定 了【)_FNN软测量模型的输入和输出之后,采用基 于EKF参数调整方法建立D—FNN软测量模型。 3.3基于模型迁移的软测量模型重构 PVC聚合过程是典型的间歇式生产过程,每 釜之间反应过程本质是一致的,但由于工况不同而 有所区别。因此,本文采用Lu等口 提出的基于过 程相似性的模型迁移建模方法对软测量模型进行校 正,即在已经建立好的基模型上利用较少的新数据 开发与其相似过程的新模型,其基本原理如图4 所示。 本文利用输入输出偏差修正(IOSBC)方法进 行迁移建模 ,其基本原理如图5所示。 假设原始的软测量模型为 Yb 。 一,(xh 。 ) (23) 式中x 和y 是原模型的输入和输出。 单独经过输出空问迁移和规划获得新模型 Y…一SoYb +Bo—Sof(Xb一)+Bo (24) 式中 s。和B。代表原模型输出空间上的尺度因子 图4基于过程相似性的模型迁移建模方法 Fig.4 Model migration modeling method based on process similarity 图5 IOSBC迁移建模原理 Fig.5 IOSBC migration modeling principle 与偏移参数。然后,对输入空间进行偏移与修正, 新模型的输入x…可以由原模型输入x 通过输入 偏移修正(ISBC)获得。 X =SI X +BI (25) 式中s 和BI表示输入空间的尺度因子与偏移参 数。所以,新模型由输入输出偏移修正(IOSBC) 原模型获得。 Y…一Sof(SIXn +BI)+B0 (26) 新模型的数据可以用来训练修正参数 MinJ(S0,B0,SI,BI)一ee s.t. e 一Y 一Es0f(sI Xn… +BI)+Bo] (27) 式中Y 是新过程的第i个观察数据,e代表新过 程的测量值与新模型之间的预测误差。应该表明的 是,所有识别的修正和规划参数的维数大小是由输 入输出空间维数确定的。 4 实验结果及分析 本文以某化工厂4万吨/年PVC生产装置的聚 合工序为背景,它是由美国古德里奇(BFG)公司 引进的技术,以VCM为原料,采用悬浮法聚合工 艺生产PVC树脂。基于D-FNN算法建立了PVC 生产过程中VCM转化速率的软测量模型。收集 PVC聚合过程部分生产历史数据,选出具有均匀 性和代表性的2釜共640组历史数据。然后把数据 预处理后的数据分为两部分,前500组数据作为训 练数据,后140组数据作为本实验的测试数据。 ・2j68・ 化工学报 第63卷 input samples 图6模糊规则 Fig.6 Fuzzy rules 8 旦 舌 舀 占 8 o o o o input samples 4 3 2 1 O加 ●2 3 4 5 fa1 D’FNN predictive output based on EKF 南 g 旦 善 琶 暑 input samples (b)D—FNN predictive output based on LLS mput samples (c)D—FNN predictive output based on KF 图7 VCM转化速率的预测输出 Fig.7 Predictive outputs of VCM conversion velocity 140组测试数据中前100组来自同一个聚合釜聚合 过程,后40组来自新的聚合过程,用来对模型迁 6 号 善 暑 旨 input samples 图8 VCM转化速率的预测误差 Fig.8 Predictive error of VCM convention velocity 移方法进行验证。 图6是3种方法产生的模糊规则数, 图7(a)~(c)分别是基于LLS、KF和EKF参数 调节方法的D-FNN软测量模型的预测输出与实际 输出的比较,图8是3种模型的预测输出误差。 从以上仿真结果可以看出,3种参数调节方法 下的D-FNN软测量模型的预测精度都非常高,其 中基于EKF的D-FNN软测量模型预测误差最小, 最适合用于聚合过程VCM转化速率的预测。从图 8可以明显看出聚合过程发生变化时,即在新的一 釜聚合数据预测时,预测误差比之前的要大一些, 但是在基于模型迁移思想实现软测量模型的在线自 适应校正的情况下,也得到更为准确的预测结果。 总体来说,动态模糊神经网络的学习速度较快,泛 化能力和鲁棒性都较好,具有广泛的应用价值。 5 结 论 采用D—FNN建立了PVC聚合生产过程中的 氯乙烯转化速率的软测量模型,并基于模型迁移思 想实现了PVC聚合生产过程这种典型间歇式复杂 工业过程关键技术指标软测量模型的在线自适应校 正及软测量模型的重构。仿真结果表明该模型能够 显著提高PVC聚合过程中经济技术指标的预测精 度,可以满足聚合釜生产过程的实时控制要求。 RefeFences [1] Gao Shuzhi(高淑芝),Gao Xianwen(高宪文),Wang Jiesheng(王介生).Soft—sensor modeling of polymerizing process based on multi T—S fuzzy neural network optimized with harmony search algorithm[J].Chinese Journal of Scientific Instrument(仪器仪表学报),2011,32(3):495— 500 [2] Wu Shiqian(伍世虔),Xu Jun(徐军).Dynamic Fuzzy Neural Network--Design and Application(动态模糊神经网 络——设计与应用)[M].BeOing:Tsinghua University 第7期 王介生等:基于D-FNN的聚合过程转化速率软测量建模及重构 ・2j69・ Press,2008:27—35 -I3] Wu S Q,Meng Joo Er.Dynamic fuzzy neural networks--a novel approach to function approximation[J].IEEE Transactions on Systems,Man,and Cybernetics,Part B: Cybernetics,2000,30(2):358 364 E4] Faa—Jeng Lin,Chih—Hong Lin,Po—Hung Shen.Self- constructing fuzzy neural network speed controller for permanent—magnet synchronous motor drive[J].Fuzzy Systems,IEEE Transactions on Fuzzy Systems,2001,9 (5):751-759 Es] Liu Ruilan(刘瑞兰),Su Hongye(苏宏业),Chu Jian (褚健).A soft sensor modeling algorithm based on modified fuzzy neural network[J].hformation and Control(信息与控制),2003,32(4):367—370 [6] Chen Guochu(陈国初),Xu Yufa(徐余法),Yu Jinshou (俞金寿).Soft sensor modelling of acrylonitrile yield based on particle swarm optimization fuzzy neural networks[J]. Journal of System Simulation(系统仿真学报),2007,19 (23):5370—5372 [7] Ma Yong(马勇),Huang Dexian(黄德先),Jin Yihui(金 以慧).Discussion about dynamic soft—sensing modeling [J].Journal of Chemical Industry and Engineering (China)(化工学报),2005,56(8):1516 1519O E83 IA Xiuliang(李修亮),Su Hongye(苏宏业),Chu Jian(褚 健).Multiple models soft—sensing technique based on online clustering arithmetic[J].Journal of Chemical Industry and Engineering(China)(化工学报),2007,s8(u): Z834一Z839 Eg] Kasabov N K, Qun Song.DENFIS:dynamic evolving neural—fuzzy inference system and its application for time series prediction[J].IEEE Transactions on Fuzzy Systems,2002,10(2):144—154 [10] Chia—Feng Juang.A TSK—type recurrent fuzzy network for dynamic systems processing by neural network and genetic algorithms[J].IEEE Transactions on Fuzzy Systems, 2002,10(2):155 17O [11] Wu S Q,Meng Joo Er.Gao Y.A fast approach for automatic generation of fuzzy rules by generalized dynamic fuzzy neural networks[J].IEEE Transactions on Fuzzy Systems,2001,9(4):578—594 [12] Gao L J,Chua K S,Chong W K,et a1.A comparison of PCA,KPCA and ICA for dimensionality reduction in support vector machine[J].Neurocomputing,2003,55 (1/2):32l一336 [13] Jong—Min Lee,ChangKyoo Yoo,Sang Wook Choi,et a1. Nonlinear process monitoring using kernel principal component analysis[J].Chemical Engineering Science, 2004,59(1):223-234 [14] Lu J D,Gao F R.Process modeling based on process similarity[J].Industrial&Engineering Chemistry Research,2008,47(63):1967—1974 [15] Lu J D,Yao K,Gao F R.Process similarity and developing new process models through migration[J].AIChE Journal,2009,S5(9):2318-2328