发明名称 最优混合丙烯聚合生产过程最优软测量仪表和方法
摘要 本发明公开了一种最优混合丙烯聚合生产过程最优软测量仪表,包括丙烯聚合生产过程、现场智能仪表、控制站、存放数据的DCS数据库、基于最优混合寻优的最优软测量仪以及熔融指数软测量值显示仪。现场智能仪表及控制站与丙烯聚合生产过程相连,与DCS数据库相连;最优软测量仪与DCS数据库及软测量值显示仪相连。所述的基于最优混合寻优的最优软测量仪包括模型更新模块、数据预处理模块、PCA主成分分析模块、神经网络模型模块以及最优混合优化模块。以及提供了一种用软测量仪表实现的软测量方法。本发明实现在线测量、在线参数优化、软测量速度快、模型自动更新、抗干扰能力强、精度高。
申请公布号 CN103838955A 申请公布日期 2014.06.04
申请号 CN201310658777.2 申请日期 2013.12.09
申请人 浙江大学 发明人 刘兴高;李九宝
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 周烽
主权项 1.一种最优混合丙烯聚合生产过程最优软测量仪表,包括丙烯聚合生产过程、用于测量易测变量的现场智能仪表、用于测量操作变量的控制站、存放数据的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)、神经网络模型模块,用于采用RBF神经网络、通过误差函数最小化来完成输入到输出的一种高度非线性映射,映射中保持拓扑不变性;(4)、最优混合优化模块,用于采用基于改进粒子群-模拟退火(IPSO-SA)最优混合算法的优化模块对神经网络进行优化,包括:(4.1)算法初始化,根据待优化的RBF神经网络结构构造出初始的粒子群体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)通过下式执行IPSO算法,让粒子群收敛:<maths num="0001"><![CDATA[<math><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>]]></maths><maths num="0002"><![CDATA[<math><mrow><msubsup><mi>v</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mi>i</mi></msubsup><mo>=</mo><msub><mi>w</mi><mi>k</mi></msub><msubsup><mi>v</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>]]></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>分别为随机数。这里的IPSO是改进的粒子群算法,其特征在于执行粒子群移动的公式中的w、c1和c2是自适应变化的;(4.3)待粒子群收敛到合适的程度时,开始对各个粒子的历代最优解组成的集合(best<sup>1</sup>,best<sup>2</sup>,…,best<sup>N</sup>)执行SA算法,得到的最优解作为算法的最后结果返回。所述基于最优混合寻优的最优软测量仪还包括:模型更新模块,用于模型的在线更新,定期将离线化验数据输入到训练集中,更新神经网络模型。所述的基于最优混合寻优的丙烯聚合生产过程最优软测量仪表,其特征在于:在所述最优混合优化模块中,粒子群执行迭代时所用到的w、c1和c2的自适应变化式为:<maths num="0003"><![CDATA[<math><mrow><mi>w</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>w</mi><mi>min</mi></msub><mo>+</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>iter</mi><mi>max</mi></msub><mo>-</mo><mi>k</mi></mrow><msub><mi>iter</mi><mi>max</mi></msub></mfrac><mo>)</mo></mrow><mi>&alpha;</mi></msup><mrow><mo>(</mo><msub><mi>w</mi><mi>max</mi></msub><mo>-</mo><msub><mi>w</mi><mi>min</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mi>c</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>c</mi><mrow><mn>1</mn><mi>min</mi></mrow></msub><mo>+</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>iter</mi><mi>max</mi></msub><mo>-</mo><mi>k</mi></mrow><msub><mi>iter</mi><mi>max</mi></msub></mfrac><mo>)</mo></mrow><mi>&beta;</mi></msup><mrow><mo>(</mo><msub><mi>c</mi><mrow><mn>1</mn><mi>max</mi></mrow></msub><mo>-</mo><msub><mi>c</mi><mrow><mn>1</mn><mi>min</mi></mrow></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>c<sub>2</sub>(k)=0.1                     (5)式中,k为算法迭代代数,w<sub>min</sub>、w<sub>max</sub>、c<sub>1min</sub>、c<sub>1max</sub>分别是常数,表示w和c<sub>1</sub>的最大、最小值,iter<sub>max</sub>为设定的算法最大迭代代数,α、β为常数。这样的自适应设计目的在于使得粒子群能够找到整个解空间内最多的局部最优解,为后续的SA算法利用这些局部最优解进一步发掘找到真正的全局最优解做好准备。所述的基于最优混合寻优的丙烯聚合生产过程最优软测量仪表,其特征在于:PCA主成分分析模块中,PCA方法实现输入变量的预白化处理,能够简化神经网络模型的输入变量,进而提高模型的性能。
地址 310058 浙江省杭州市西湖区余杭塘路866号