发明名称 模糊智能最优丙烯聚合生产过程最优软测量仪表和方法
摘要 本发明公开了一种模糊智能最优丙烯聚合生产过程最优软测量仪表,包括丙烯聚合生产过程、现场智能仪表、控制站、存放数据的DCS数据库、模糊智能最优软测量仪以及熔融指数软测量值显示仪。现场智能仪表及控制站与丙烯聚合生产过程相连,与DCS数据库相连;最优软测量仪与DCS数据库及软测量值显示仪相连。所述的模糊智能最优软测量仪包括模型更新模块、数据预处理模块、PCA主成分分析模块、神经网络模型模块以及模糊智能优化模块。以及提供了一种用软测量仪表实现的软测量方法。本发明实现在线测量、在线参数优化、软测量速度快、模型自动更新、抗干扰能力强、精度高。
申请公布号 CN103838958B 申请公布日期 2017.02.22
申请号 CN201310659448.X 申请日期 2013.12.09
申请人 浙江大学 发明人 刘兴高;李九宝
分类号 G06F19/00(2011.01)I;G06N3/02(2006.01)I 主分类号 G06F19/00(2011.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 周烽
主权项 一种模糊智能最优丙烯聚合生产过程最优软测量仪表,包括丙烯聚合生产过程、用于测量易测变量的现场智能仪表、用于测量操作变量的控制站、存放数据的DCS数据库、模糊智能最优软测量仪以及熔融指数软测量显示仪,所述现场智能仪表、控制站与丙烯聚合生产过程连接,所述现场智能仪表、控制站与DCS数据库连接,所述DCS数据库与模糊智能最优软测量模型的输入端连接,所述模糊智能最优软测量模型的输出端与熔融指数软测量显示仪连接,其特征在于:所述模糊智能最优软测量仪包括:(1)、数据预处理模块,用于将从DCS数据库输入的模型输入变量进行预处理,对输入变量中心化,即减去变量的平均值;再进行归一化处理,即除以变量值的变化区间;(2)、PCA主成分分析模块,用于将输入变量预白化处理及变量去相关,通过对输入变量施加一个线性变换实现,即主成分由C=MU得到,其中M为输入变量,C为主成分得分矩阵,U为载荷矩阵;若对原始数据进行重构,可由M=CU<sup>T</sup>计算,其中上标T表示矩阵的转置;当选取的主成分数目小于输入变量的变量个数时,M=CU<sup>T</sup>+E,其中E为残差矩阵;(3)、神经网络模型模块,采用一个四层的模糊神经网络、通过误差函数最小化来完成输入到输出的一种高度非线性映射,映射中保持拓扑不变性;(4)、模糊智能优化模块,用于采用模糊智能优化模块对神经网络进行优化,包括:(4.1)算法初始化,根据待优化的模糊神经网络结构参数构造出初始的粒子群体X=(x<sup>1</sup>,x<sup>2</sup>,…,x<sup>N</sup>),初始移动速度V=(v<sup>1</sup>,v<sup>2</sup>,…,v<sup>N</sup>),初始的各粒子历代最优值OP=(p<sup>1</sup>,p<sup>2</sup>,…,p<sup>N</sup>)和全局最优值p<sup>g</sup>;(4.2)通过下式执行粒子群优化算法,让粒子群收敛:<maths num="0001"><math><![CDATA[<mrow><msubsup><mi>x</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mi>i</mi></msubsup><mo>=</mo><msubsup><mi>x</mi><mi>k</mi><mi>i</mi></msubsup><mo>+</mo><msubsup><mi>v</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mi>i</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001111806600000011.GIF" wi="1110" he="71" /></maths><maths num="0002"><math><![CDATA[<mrow><msubsup><mi>v</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mi>i</mi></msubsup><mo>=</mo><msubsup><mi>wv</mi><mi>k</mi><mi>i</mi></msubsup><mo>+</mo><msub><mi>c</mi><mn>1</mn></msub><msub><mi>r</mi><mn>1</mn></msub><mrow><mo>(</mo><msubsup><mi>p</mi><mi>k</mi><mi>i</mi></msubsup><mo>-</mo><msubsup><mi>x</mi><mi>k</mi><mi>i</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><msub><mi>r</mi><mn>2</mn></msub><mrow><mo>(</mo><msubsup><mi>p</mi><mi>k</mi><mi>g</mi></msubsup><mo>-</mo><msubsup><mi>x</mi><mi>k</mi><mi>i</mi></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001111806600000012.GIF" wi="1390" he="71" /></maths>式中x为粒子的位置向量,i为粒子的序号,k为算法迭代代数,v为粒子速度,p表示初始的各粒子历代最优值OP=(p<sup>1</sup>,p<sup>2</sup>,…,p<sup>N</sup>)和全局最优值p<sup>g</sup>的最优值集合;v<sub>k</sub><sup>i</sup>为第k次迭代中第i个粒子的速度;x<sub>k</sub><sup>i</sup>为第k次迭代中第i个粒子的位置向量;p<sub>k</sub><sup>i</sup>为第k次迭代中第i个粒子的历代最优位置,p<sup>i</sup>为第i个粒子的历代最优解,w为速度权重系数,c<sub>1</sub>、c<sub>2</sub>分别为粒子历代最优解和群体最优解的吸引系数,r<sub>1</sub>、r<sub>2</sub>分别为随机数;(4.3)所述粒子群优化算法体现在速度权重值w是基于一种进化状态来进行估计的;粒子群的进化状态包括了探索期、开拓期、聚合期和跳出期;由f来定量表示;<maths num="0003"><math><![CDATA[<mrow><mi>f</mi><mo>=</mo><mfrac><mrow><msub><mi>d</mi><mi>g</mi></msub><mo>-</mo><msub><mi>d</mi><mi>min</mi></msub></mrow><mrow><msub><mi>d</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>-</mo><msub><mi>d</mi><mi>min</mi></msub></mrow></mfrac><mo>&Element;</mo><mo>&lsqb;</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001111806600000021.GIF" wi="1302" he="142" /></maths>其中,f表示粒子群的进化状态,d<sub>i</sub>是粒子i到粒子群中其他粒子的平均距离;d<sub>g</sub>是粒子群中最优解到粒子群中其他粒子的平均距离;d<sub>max</sub>和d<sub>min</sub>分别是{d<sub>i</sub>}中的最大、最小值;d<sub>i</sub>按照下式计算:<maths num="0004"><math><![CDATA[<mrow><msub><mi>d</mi><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mi>P</mi><mo>-</mo><mn>1</mn></mrow></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>&NotEqual;</mo><mi>i</mi></mrow><mi>P</mi></munderover><mo>|</mo><mo>|</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>x</mi><mi>j</mi></msub><mo>|</mo><msub><mo>|</mo><mn>2</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001111806600000022.GIF" wi="1150" he="135" /></maths>x<sub>i</sub>为第i个粒子的位置,x<sub>j</sub>为第j个粒子的位置,‖·‖为范数表达式;P为粒子群中粒子的数目;由f值来得到粒子群速度权重的计算式为:<maths num="0005"><math><![CDATA[<mrow><mi>w</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><mn>1.5</mn><msup><mi>e</mi><mrow><mo>-</mo><mn>2.6</mn><mi>f</mi></mrow></msup></mrow></mfrac><mo>&Element;</mo><mo>&lsqb;</mo><mn>0.4</mn><mo>,</mo><mn>0.9</mn><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001111806600000023.GIF" wi="1174" he="118" /></maths>所述模糊智能最优软测量仪还包括:模型更新模块,用于模型的在线更新,定期将离线化验数据输入到训练集中,更新神经网络模型;所述模糊神经网络输入层到第二层的传递函数使用的是高斯函数:<maths num="0006"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msubsup><mi>u</mi><mrow><mi>q</mi><mi>j</mi></mrow><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></msubsup><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msup><mrow><mo>&lsqb;</mo><msubsup><mi>u</mi><mi>q</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mo>-</mo><msub><mi>m</mi><mrow><mi>q</mi><mi>j</mi></mrow></msub><mo>&rsqb;</mo></mrow><mn>2</mn></msup><msubsup><mi>&sigma;</mi><mrow><mi>q</mi><mi>j</mi></mrow><mn>2</mn></msubsup></mfrac><mo>)</mo></mrow><mo>,</mo></mrow></mtd><mtd><mrow><mi>q</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>n</mi><mo>,</mo></mrow></mtd><mtd><mrow><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>N</mi></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001111806600000024.GIF" wi="1477" he="165" /></maths>式中u<sup>(1)</sup>、u<sup>(2)</sup>分别为第二层的输入和输出向量,m为中心向量,σ为宽度参数,q、j分别为维度和节点序号;使用高斯非线性函数能够增强该模糊神经网络的非线性拟合能力;所述PCA主成分分析模块中,PCA方法实现输入变量的预白化处理,能够简化神经网络模型的输入变量,进而提高模型的性能。
地址 310058 浙江省杭州市西湖区余杭塘路866号