发明内容
为了克服已有聚丙烯熔融指数在线预报方法的操作复杂、在线投运困难、实用性差的不足,本发明提供一种设计简单、易于在线投运、实用性强的聚丙烯装置熔融指数在线动态预报方法。
本发明解决其技术问题所采用的技术方案是:
一种聚丙烯熔融指数在线动态预报方法,所述预报方法包括如下步骤:
1)、建立聚丙烯熔融指数Q阶离散形式的参数模型结构,参见式(1):
其中,正整数Q是模型式(1)的阶次;b0(k)为未知参数,由本步骤第3)步辨识计算得到;y(k)为第k时刻聚丙烯熔融指数;u1(k-d)和u2(k-d)分别表示第k-d时刻聚合装置1#管的氢气输入量与丙烯输入量之比和聚合反应温度;u3(k-d)和u4(k-d)分别表示第k-d时刻聚合装置2#管的氢气输入量与丙烯输入量之比和聚合反应温度;d为聚合物平均反应停留时间;e(k)表示第k时刻均值为0的不相关随机噪声;系数A(q-1)和Bij...t(q-1)是m阶多项式,标号i=1,...,4;j=1,...,i;...;t=1,...,s;系数C(q-1)是n阶多项式。多项式A(q-1)、Bij...t(q-1)和C(q-1)参见式(2):
A(q-1)=a1(k)q-1+…+am(k)q-m
Bij...t(q-1)=bij...t,1(k)q-1+…+bij...t,m(k)q-m (2)
C(q-1)=1+c1(k)q-1+…+cn(k)q-n
式中,q-1表示时延算子符号;系数a1(k),...am(k),c1(k),...,cn(k),bij...t,1(k),...,bij...t,m(k),(i=1,...,4;j=1,...,i;...;t=1,...,s)为模型式(1)在第k时刻的未知参数,由本步骤第3)步在线辨识计算得到。
2)、对操作参数的数据作归一化处理,参见式(3):
其中,u′为归一化后的值;uh和ul分别表示该操作参数的变化上限和下限,是通过对操作数据统计分析得到。
3)、运用渐消记忆增广最小二乘算法递推辨识熔融指数模型式(1)在第k时刻的参数b0(k),a1(k),...am(k),c1(k),...,cn(k),bij...t,1(k),...,bij...t,m(k),(i=1,...,4;j=1,...,i;...;t=1,...,s),其算法参见式(4):
θ(k)=θ(k-1)+K(k)[y(k)-h(k)Tθ(k-1)],
其中,λ∈(0,1]为遗忘因子;θ(k)为第k时刻模型式(1)的辨识参数向量;K(k)表示第k时刻校正系数矩阵;P(k)表示第k时刻记忆系数矩阵;h(k)表示第k时刻过程参量的历史数据向量。θ(k)和h(k)表达式分别如下:
θ(k)T=[b0(k),b1,1(k),…,b1,m(k),…,bij...t,1(k),…,bij...t,m(k),
a1(k),…,am(k),c1(k),…,cn(k)],
其中,符号“T”表示向量的转置;
表示第s时刻的随机噪声的估计,s=k,k-1,...,k-n;变量u
ij...t(k-d-l)=u
i(k-d-l)u
j(k-d-l)…u
t(k-d-l),l=1,..,m。
4)、在线测量丙烯输入量、氢气输入量和聚合反应温度,并计算氢气输入量和丙烯输入量之比,再根据在线辨识的参数θ和熔融指数预报模型式(1),在线动态预报聚丙烯生产装置中产品的熔融指数值。
熔融指数在线预报器在每次实验室化验结束后自动启动模型式(1)参数更新程序,继而根据最新的模型参数,在每个采样时刻在线测量计算氢气输入量和丙烯输入量之比和聚合反应温度,利用模型式(1)在线预报聚丙烯产品中的熔融指数。周而复始,实现聚丙烯生产过程中熔融指数的在线动态预报。
本发明的技术构思为:以目前广泛采用的Spheripol丙烯聚合装置为对象(生产工艺流程如图1所示),建立聚丙烯熔融指数Q阶离散形式模型结构,并运用渐消记忆增广最小二乘辨识方法,结合工业现场数据在线更新熔融指数的模型参数,继而通过实时测量计算氢气输入量和丙烯输入量之比和聚合反应温度,实现聚丙烯熔融指数的在线动态预报操作。本发明设计方法的优点是容易理解、使用方便,可以在线动态预报诸如牌号切换等非平稳生产过程中的聚丙烯熔融指数值。本发明与何德峰,2010中聚丙烯熔融指数检测方法的区别在于,该聚丙烯熔融检测方法是针对特定的聚丙烯牌号建立熔融指数聚合机理模型的设计方法,并利用观测器实现预报,设计复杂,专业知识抽象;与陈果,2006;李文义,2008中聚丙烯熔融指数检测方法的区别在于,这些聚丙烯熔融检测方法是基于聚合动力学非线性机理模型,并应用序列二次规划算法实现对平稳生产过程中聚丙烯熔融指数的设计方法。
本发明主要执行部分在工业控制计算机上运行实施。本方法应用过程可以大致分为3个阶段:
1、预报器参数设置,即在组态界面中输入熔融指数模型的阶次Q,多项式A(q-1)、Bij...t(q-1)和C(q-1)的阶次m和n,遗忘因子λ,聚丙烯生产装置平均反应停留时间d,预报周期T1,更新周期T2。在界面中可以通过点击“+”和“-”完成参数的设置。设置参数确认后,点击“模型更新”开始熔融指数预报模型参数的更新操作,同时由工控机将数据送入动态数据库中保存。
2、熔融指数模型参数的辨识。启动工控机的CPU调用事先编制好的“模型参数辨识程序”软件包在线辨识熔融指数模型参数,包括自动调用动态数据库存储的生产数据和预报器的设置参数。模型参数更新完毕后,在组态界面的上框显示模型参数,下框显示本次参数辨识的效果。
点击组态界面中的“调试”按钮,进入模型参数辨识调试阶段。调整组态界面中的多项式A(q-1)、Bij...t(q-1)和C(q-1)的阶次m和n、遗忘因子λ等参数,并观察辨识效果,由此确定一组能良好反应实际熔融指数变化的模型参数。参数Q、m和n的调整规则:增大参数Q、m和n将提高辨识的效果,但增加预报模型的复杂度和辨识所需的操作数据量,从而增大辨识算法的在线计算量;相反,调小参数Q、m和n可以简化预报模型的复杂度,从而降低在线辨识的操作数据量,提高在线参数辨识的速度,但算法所具有的信息量少可恶化辨识效果。参数λ的调节规则:通常取0.95<λ<1;调小λ即对历史数据“遗忘”减少,历史数据比重增大;相反,增大λ即对历史数据“遗忘”增加,当前数据比重增大。因此,实际调试Q、m、n和λ参数时,应在辨识效果、模型复杂度和所需操作数据量等性能之间权衡。
3、点击组态界面“在线预报”按钮,启动工控机的CPU读取最佳熔融指数模型参数,并执行“熔融指数在线预报程序”预报当前时刻聚丙烯产品中熔融指数。此时组态界面上显示聚丙烯熔融指数动态预报曲线。在下一个采样预报时刻到来时,在线检测计算氢气输入量和丙烯输入量之比和聚合反应温度,之后重复整个执行过程。如此周而复始,实现聚丙烯生产装置熔融指数的在线动态预报。
全套聚丙烯熔融指数在线预报过程可以在工控机组态界面上完成,此过程可以参考本说明书下文中提供的工业实例应用。与传统聚丙烯熔融指数预报方法相比,本发明给出的聚丙烯熔融指数在线预报方法的最大特点是可以根据系统运行的实时信息实现熔融指数的在线动态预报。下文具体实施方法以聚丙烯生产过程为例说明本发明的实际效果,但本发明的应用范围并不以本实施例中的熔融指数在线预报为限。如前所述,本发明除了可以用于聚丙烯生产过程外,还可用于聚乙烯、聚氯乙烯等各类聚合物熔融指数的预报操作。
本发明的有益效果主要表现在:1、设计简单、容易理解、在线实施简便、实用性强;2、在相当宽的范围内实现聚丙烯熔融指数特别是非平稳生产过程的在线动态预报,可极大地提高聚丙烯生产的技术水平,提高企业生产经济效益和市场竞争力。
具体实施方式
下面结合附图对本发明作进一步描述。
一种聚丙烯熔融指数在线动态预报方法,所述预报方法包括如下步骤:
1)、建立聚丙烯熔融指数Q阶离散形式的参数模型结构,参见式(1):
其中,正整数Q是模型式(1)的阶次;b0(k)为未知参数,由本步骤第3)步辨识计算得到;y(k)为第k时刻聚丙烯熔融指数;u1(k-d)和u2(k-d)分别表示第k-d时刻聚合装置1#管的氢气输入量与丙烯输入量之比和聚合反应温度;u3(k-d)和u4(k-d)分别表示第k-d时刻聚合装置2#管的氢气输入量与丙烯输入量之比和聚合反应温度;d为聚合物平均反应停留时间;e(k)表示第k时刻均值为0的不相关随机噪声;系数A(q-1)和bij...t(q-1)是m阶多项式,标号i=1,...,4;j=1,...,i;...;t=1,...,s;系数C(q-1)是n阶多项式。多项式A(q-1)、Bij...t(q-1)和C(q-1)参见式(2):
A(q-1)=a1(k)q-1+…+am(k)q-m
Bij...t(q-1)=bij...t,1(k)q-1+…+bij...t,m(k)q-m (2)
C(q-1)=1+c1(k)q-1+…+cn(k)q-n
式中,q-1表示时延算子符号;系数a1(k),...am(k),c1(k),...,cn(k),bij...t,1(k),...,bij...t,m(k),(i=1,...,4;j=1,...,i;...;t=1,...,s)为模型式(1)在第k时刻的未知参数,由本步骤第3)步在线辨识计算得到。
2)、对操作参数的数据作归一化处理,参见式(3):
其中,u′为归一化后的值;uh和ul分别表示该操作参数的变化上限和下限,是通过对操作数据统计分析得到。
3)、运用渐消记忆增广最小二乘算法递推辨识熔融指数模型式(1)在第k时刻的参数b0(k),a1(k),...am(k),c1(k),...,cn(k),bij...t,1(k),...,bij...t,m(k),(i=1,...,4;j=1,...,i;...;t=1,...,s),其算法参见式(4):
θ(k)=θ(k-1)+K(k)[y(k)-h(k)Tθ(k-1)],
其中,λ∈(0,1]为遗忘因子;θ(k)为第k时刻模型式(1)的辨识参数向量;K(k)表示第k时刻校正系数矩阵;P(k)表示第k时刻记忆系数矩阵;h(k)表示第k时刻过程参量的历史数据向量。θ(k)和h(k)表达式分别如下:
θ(k)T=[b0(k),b1,1(k),…,b1,m(k),…,bij...t,1(k),…,bij...t,m(k),
a1(k),…,am(k),c1(k),…,cn(k)],
其中,符号“T”表示向量的转置;表示第s时刻的随机噪声的估计,s=k,k-1,...,k-n;变量uij...t(k-d-l)=ui(k-d-l)uj(k-d-l)…ut(k-d-l),l=1,..,m。
4)、在线测量丙烯输入量、氢气输入量和聚合反应温度,并计算氢气输入量和丙烯输入量之比,再利用在线辨识参数θ和熔融指数预报模型式(1),在线动态预报聚丙烯生产装置中产品的熔融指数值。
熔融指数在线预报器在每次实验室化验结束后自动启动模型式(1)参数更新程序,继而根据最新的模型参数,在每个采样时刻在线测量计算氢气输入量和丙烯输入量之比和聚合反应温度,利用模型式(1)在线预报聚丙烯产品中的熔融指数值。周而复始,实现聚丙烯生产过程中熔融指数值的在线动态预报。
本实施例为聚丙烯生产过程,具体操作过程:
1、在预报器参数设置界面中,通过点击“+”或“-”键完成预报器各参数值的设置,如下:阶次Q=2,m=2,n=1,遗忘因子λ=0.98,平均反应停留时间d=60分钟,预报周期T1=10分钟,更新周期T2=120分钟。
2、在组态界面上点击“模型更新”按钮进入下一组态界面,启动工控机的CPU调用事先编制好的“模型参数辨识程序”软件包在线辨识熔融指数模型参数。具体计算过程如下:
1)根据预报器参数设置,由公式(1)可得聚丙烯熔融指数2阶离散形式模型结构为
其中,b0(k)为未知参数,由本实施例步骤第3)步辨识计算得到;y(k)为第k时刻聚丙烯熔融指数;u1(k-d)和u2(k-d)分别表示第k-d时刻聚合装置1#管的氢气输入量与丙烯输入量之比和聚合反应温度;u3(k-d)和u4(k-d)分别表示第k-d时刻聚合装置2#管的氢气输入量与丙烯输入量之比和聚合反应温度;d为聚合物平均反应停留时间;e(k)表示第k时刻均值为0的不相关随机噪声;系数A(q-1)和Bij(q-1)是m阶多项式,标号i=1,...,4;j=1,...,i;系数C(q-1)是n阶多项式。多项式A(q-1)、Bij(q-1)和C(q-1)分别如下:
A(q-1)=a1(k)q-1+…+am(k)q-m
Bij(q-1)=bij,1(k)q-1+…+bij,m(k)q-m
C(q-1)=c1(k)q-1+…+cn(k)q-n
式中,q-1表示时延算子符号;系数a1(k),...am(k),c1(k),...,cn(k),bij,1(k),...,bij,m(k),(i=1,...,4;j=1,...,i)为本实施例熔融指数模型在第k时刻的未知参数,由本实施例步骤第3)步在线辨识计算得到。
2)对操作参数的数据做统计分析,得该变量的变化范围为uh和ul,据此应用公式(3)对操作数据作归一化处理,得
并将这些标定数据送入工控机动态数据库系统保存。
3)运用渐消记忆增广最小二乘算法递推辨识熔融指数模型在k时刻的参数b0(k),a1(k),...am(k),c1(k),...,cn(k),bij,1(k),...,bij,m(k),(i=1,...,4;j=1,...,i),得
θ(k)=θ(k-1)+K(k)[Δy(k)-h(k)Tθ(k-1)],
其中,遗忘因子λ=0.98;K(k)表示第k时刻校正系数矩阵;P(k)表示第k时刻记忆系数矩阵;第k时刻辨识参数向量θ(k)和过程参量历史数据向量h(k)分别如下:
θ(k)T=[b0(k),b1,1(k),…,b1,m(k),…,bij,1(k),…,bij,m(k),
a1(k),…,am(k),c1(k),…,cn(k)],
其中,符号“T”表示向量的转置;
表示第s时刻的随机噪声的估计,s=k,k-1,...,k-n;变量u
ij(k-d-l)=u
i(k-d-l)u
j(k-d-l),l=1,..,m。参数辨识算法在每次产品化验后利用最新熔融指数数据更新其模型参数θ,再利用参数更新后的熔融指数动态数学模型式(1),在线动态预报聚丙烯装置的熔融指数。
4)调试预报器的设计参数m、n和λ,得m=2、n=1和λ=0.998。
3、点击组态界面“在线预报”按钮,启动工控机的CPU读取最佳熔融指数模型参数,并执行“熔融指数在线预报程序”预报当前时刻聚丙烯产品中熔融指数。此时组态界面上显示聚丙烯熔融指数动态实时预报曲线。在下一个采样预报时刻到来时,在线测量计算氢气输入量与丙烯输入量之比和聚合反应温度,再重复整个执行过程,如此周而复始,实现聚丙烯熔融指数的在线动态预报。
以上阐述的是本发明给出的一个实施例所表现出的优良熔融指数在线动态预报效果。需要指出,上述实施例用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明做出的任何修改,都落入本发明的保护范围。