发明名称 一种聚丙烯熔融指数在线动态预报方法
摘要 本发明公开了一种聚丙烯熔融指数的动态在线预报方法。通过点击工控机组态界面上的按钮启动事先编写好的熔融指数在线预报的程序,以聚丙烯全流程工艺熔融指数动态预报数学模型为基础,结合渐消记忆增广最小二乘法实现聚丙烯熔融指数的在线预报。本发明的最重要的特征是在聚合装置连续生产的状态下,仅通过在线测量计算氢气与丙烯输入流量比值和聚合反应温度,即可在相当宽的范围内在线预报聚丙烯树脂的熔融指数,具有简便有效、预报方法实施方便和低成本的优势,可极大地提高聚丙烯关键质量的检测频率和效率,降低不合格产品的排放,提高生产经济效益和市场竞争力。
申请公布号 CN102663220A 申请公布日期 2012.09.12
申请号 CN201210051707.6 申请日期 2012.03.01
申请人 浙江工业大学 发明人 何德峰;俞立;宋秀兰
分类号 G06F19/00(2006.01)I 主分类号 G06F19/00(2006.01)I
代理机构 杭州天正专利事务所有限公司 33201 代理人 王兵;黄美娟
主权项 1.一种聚丙烯熔融指数在线动态预报方法,所述预报方法包括如下步骤:1)、建立聚丙烯熔融指数Q阶离散形式的参数模型结构,参见式(1):<img file="FDA0000139937750000011.GIF" wi="760" he="74" /><img file="FDA0000139937750000012.GIF" wi="847" he="79" /><img file="FDA0000139937750000013.GIF" wi="1408" he="156" /><img file="FDA0000139937750000014.GIF" wi="490" he="54" />其中,正整数(Q是模型式(1)的阶次;b<sub>0</sub>(k)为未知参数,由本步骤第3)步辨识计算得到;y(k)为第k时刻聚丙烯熔融指数;u<sub>1</sub>(k-d)和u<sub>2</sub>(k-d)分别表示第k-d时刻聚合装置1#管的氢气输入量与丙烯输入量之比和聚合反应温度;u<sub>3</sub>(k-d)和u<sub>4</sub>(k-d)分别表示第k-d时刻聚合装置2#管的氢气输入量与丙烯输入量之比和聚合反应温度;d为聚合物平均反应停留时间;e(k)表示第k时刻均值为0的不相关随机噪声;系数A(q<sup>-1</sup>)和B<sub>ij...jt</sub>(q<sup>-1</sup>)是m阶多项式,标号i=1,...,4;j=1,...,i;...;t=1,...,s;系数C(q<sup>-1</sup>)是n阶多项式。多项式A(q<sup>-1</sup>)、B<sub>ij...t</sub>(q<sup>-1</sup>)和C(q<sup>-1</sup>)参见式(2):A(q<sup>-1</sup>)=a<sub>1</sub>(k)q<sup>-1</sup>+…+a<sub>m</sub>(k)q<sup>-m</sup>B<sub>ij...t</sub>(q<sup>-1</sup>)=b<sub>ij...t,1</sub>(k)q<sup>-1</sup>+…+b<sub>ij...t,m</sub>(k)q<sup>-m</sup>(2)C(q<sup>-1</sup>)=1+c<sub>1</sub>(k)q<sup>-1</sup>+…+c<sub>n</sub>(k)q<sup>-n</sup>式中,q<sup>-1</sup>表示时延算子符号;系数a<sub>1</sub>(k),...a<sub>m</sub>(k),c<sub>1</sub>(k),...,c<sub>n</sub>(k),b<sub>ij...t,1</sub>(k),...,b<sub>ij...t,m</sub>(k),(i=1,...,4;j=1,...,i;...;t=1,...,s)为模型式(1)在第k时刻的未知参数,由本步骤第3)步在线辨识计算得到。2)、对操作参数的数据作归一化处理,参见式(3):<maths num="0001"><![CDATA[<math><mrow><msup><mi>u</mi><mo>&prime;</mo></msup><mo>=</mo><mfrac><mrow><mi>u</mi><mo>-</mo><msub><mi>u</mi><mi>l</mi></msub></mrow><mrow><msub><mi>u</mi><mi>h</mi></msub><mo>-</mo><msub><mi>u</mi><mi>l</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,u′为归一化后的值;u<sub>h</sub>和u<sub>l</sub>分别表示该操作参数的变化上限和下限,是通过对操作数据统计分析得到。3)、运用渐消记忆增广最小二乘算法递推辨识熔融指数模型式(1)在第k时刻的参数b<sub>0</sub>(k),a<sub>1</sub>(k),...a<sub>m</sub>(k),c<sub>1</sub>(k),...,c<sub>n</sub>(k),b<sub>ij...t,1</sub>(k),...,b<sub>ij...t,m</sub>(k),(i=1,...,4;j=1,...,i;..;t=1,...,s),其算法参见式(4):θ(k)=θ(k-1)+K(k)[y(k)-h(k)<sup>T</sup>θ(k-1)],<maths num="0002"><![CDATA[<math><mrow><mi>K</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>P</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>h</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><mi>h</mi><msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mi>T</mi></msup><mi>P</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>h</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>&lambda;</mi></mrow></mfrac><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><mi>P</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>&lambda;</mi></mfrac><mo>[</mo><mi>I</mi><mo>-</mo><mi>K</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mi>h</mi><msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mi>T</mi></msup><mo>]</mo><mi>P</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>.</mo></mrow></math>]]></maths>其中,λ∈(0,1]为遗忘因子;θ(k)为第k时刻模型式(1)的辨识参数向量;K(k)表示第k时刻校正系数矩阵;P(k)表示第k时刻记忆系数矩阵;h(k)表示第k时刻过程参量的历史数据向量。θ(k)和h(k)表达式分别如下:θ(k)<sup>T</sup>=[b<sub>0</sub>(k),b<sub>1,1</sub>(k),…,b<sub>1,m</sub>(k),…,b<sub>ij...t,1</sub>(k),…,b<sub>ij...t,m</sub>(k),a<sub>1</sub>(k),…,a<sub>m</sub>(k),c<sub>1</sub>(k),…,c<sub>n</sub>(k)],<maths num="0004"><![CDATA[<math><mrow><mi>h</mi><msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mi>T</mi></msup><mo>=</mo><mo>[</mo><mn>1</mn><mo>,</mo><msub><mi>u</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mi>d</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>u</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mi>d</mi><mo>-</mo><mi>m</mi><mo>)</mo></mrow><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mi>u</mi><mrow><mi>ij</mi><mo>.</mo><mo>.</mo><mo>.</mo><mi>t</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mi>d</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>u</mi><mrow><mi>ij</mi><mo>.</mo><mo>.</mo><mo>.</mo><mi>t</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mi>d</mi><mo>-</mo><mi>m</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mi>m</mi><mo>)</mo></mrow><mo>,</mo><mover><mi>e</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mover><mi>e</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mi>n</mi><mo>)</mo></mrow><mo>]</mo><mo>,</mo></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><mover><mi>e</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mi>h</mi><msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mi>T</mi></msup><mi>&theta;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></math>]]></maths>其中,符号“T”表示向量的转置;<img file="FDA0000139937750000025.GIF" wi="403" he="58" />表示第s时刻的随机噪声的估计,s=k,k-1,...,k-n;变量u<sub>ij...t</sub>(k-d-l)=u<sub>i</sub>(k-d-l)u<sub>j</sub>(k-d-l)…u<sub>t</sub>(k-d-l),l=1,..,m。4)、在线测量丙烯输入量、氢气输入量和聚合反应温度,并计算氢气输入量和丙烯输入量之比,再根据在线辨识的参数θ和熔融指数预报模型式(1),在线动态预报聚丙烯生产装置中产品的熔融指数值。熔融指数在线预报器在每次实验室化验结束后自动启动模型式(1)参数更新程序,继而根据最新的模型参数,在每个采样时刻在线测量计算氢气输入量和丙烯输入量之比和聚合反应温度,利用模型式(1)在线预报聚丙烯产品中的熔融指数。周而复始,实现聚丙烯生产过程中熔融指数的在线动态预报。
地址 310014 浙江省杭州市下城区潮王路18号