发明名称 丙烯聚合生产过程自适应最优预报系统和方法
摘要 本发明公开了一种丙烯聚合生产过程自适应最优预报系统,包括丙烯聚合生产过程、现场智能仪表、控制站、存放数据的DCS数据库、自适应最优预报系统以及熔融指数预报值显示仪。现场智能仪表及控制站与丙烯聚合生产过程相连,与DCS数据库相连;最优预报系统与DCS数据库及预报值显示仪相连。所述的自适应最优预报系统包括模型更新模块、数据预处理模块、PCA主成分分析模块、神经网络模型模块以及自适应智能优化模块。以及提供了一种用预报系统实现的预报方法。本发明实现在线测量、在线参数优化、预报速度快、模型自动更新、抗干扰能力强、精度高。
申请公布号 CN103838209A 申请公布日期 2014.06.04
申请号 CN201310659482.7 申请日期 2013.12.09
申请人 浙江大学 发明人 刘兴高;李九宝
分类号 G05B19/418(2006.01)I;G06N3/02(2006.01)I 主分类号 G05B19/418(2006.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)、自适应智能优化模块,用于采用基于智能粒子群算法的优化模块对神经网络进行优化,包括:(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)通过下式执行智能粒子群算法,让粒子群收敛:<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><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>]]></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"><![CDATA[<math><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><mi>max</mi></msub><mo>-</mo><msub><mi>d</mi><mi>min</mi></msub></mrow></mfrac><mo>&Element;</mo><mo>[</mo><mn>0,1</mn><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>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"><![CDATA[<math><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><mi>&Sigma;</mi><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><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>x</mi><mi>j</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>x<sub>i</sub>为第i个粒子的位置,x<sub>j</sub>为第j个粒子的位置,‖·‖为范数表达式,P为粒子群中粒子的数目;由f值来得到粒子群速度权重的计算式为:<maths num="0005"><![CDATA[<math><mrow><mi>w</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><msup><mrow><mn>1.5</mn><mi>e</mi></mrow><mrow><mo>-</mo><mn>2.6</mn><mi>f</mi></mrow></msup></mrow></mfrac><mo>&Element;</mo><mo>[</mo><mn>0.4,0.9</mn><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>所述自适应最优预报系统还包括:模型更新模块,用于模型的在线更新,定期将离线化验数据输入到训练集中,更新神经网络模型。所述的自适应最优预报系统,其特征在于:在所述自适应智能优化模块中,粒子群执行迭代时所用到的w是根据粒子群的进化状态来智能变化的;这样的智能设计目的在于使得粒子群能够在不同的阶段采用不同的速度权重,控制群体的收敛速度,提高整个群体的全局寻优能力。PCA主成分分析模块中,PCA方法实现输入变量的预白化处理,能够简化神经网络模型的输入变量,进而提高模型的性能。
地址 310058 浙江省杭州市西湖区余杭塘路866号