CN102914970B - 工业过程控制系统中不能直接测量的性能参数的在线数据驱动估计方法 - Google Patents
工业过程控制系统中不能直接测量的性能参数的在线数据驱动估计方法 Download PDFInfo
- Publication number
- CN102914970B CN102914970B CN201210431022.4A CN201210431022A CN102914970B CN 102914970 B CN102914970 B CN 102914970B CN 201210431022 A CN201210431022 A CN 201210431022A CN 102914970 B CN102914970 B CN 102914970B
- Authority
- CN
- China
- Prior art keywords
- industrial process
- control system
- directly measure
- process control
- physical quantity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Abstract
工业过程控制系统中不能直接测量的性能参数的在线数据驱动估计方法,本发明涉及一种在工业过程控制系统中性能参数的数据驱动估计方法。本发明是要解决现有的工业过程控制系统中不能直接测量的性能参数依靠经验估计的方法可靠性低的技术问题。本方法:一、确定工业过程控制系统中可直接测量的物理量;二、在不能直接测量的性能参数y的一个退化周期内测量不同时刻不能直接测量的性能参数y及可直接测量的物理量的值并进行归一化处理,然后对构建的函数模型进行对数估计,再剔除冗余项,最后选择关键参数,得到工业过程控制系统中不能直接测量的性能参数y与运行时间的关系式。该方法适于对催化剂老化性能、配件老化、炉膛结渣的分析预测。
Description
技术领域
本发明涉及一种在工业过程控制系统中,对不可测的性能参数进行估计的方法。
背景技术
工业过程控制系统中,有一些性能参数不能直接、实时地测量出来,但是最终产品质量有着巨大的影响,如在石化加工等过程控制系统中,催化剂的性能对系统性能、最终产品质量有着很大的影响,而且催化剂存在退化现象,即随着使用时间的增加,催化剂性能逐渐下降,但催化剂的性能不能直接、实时地测量出来,同时,由于生产过程复杂,产品输出存在一定的滞后性,当最终产品出现质量问题时,催化剂退化问题已经对生产过程造成了巨大损失。又如金属冶炼过程中,锅炉炉膛结渣情况同样存在不可直接测量,又对最终产品的品质有重大影响。传统上,由于难于对此类反应过程机理进行建模和分析,人们只能根据实际生产经验,通过定期检测的方式来解决此类问题。这样即需要大量的实际生产经验,又无法提高生产效率。另外,由于生产过程的不确定性,使得依靠经验的方法存在较大的不可靠性。在实际生产中,类似的问题还有元器件的老化等。因此,建立一套不依赖于反应机理和生产经验的数据模型,从而有效、实时评估反应系统中催化剂性能和产品质量的方法,对于工业过程生产具有极大的意义。
发明内容
本发明是要解决现有的催化剂的性能依靠经验估计的方法可靠性低的技术问题,而提供工业过程控制系统中不能直接测量的性能参数的在线数据驱动估计方法。
本发明的工业过程控制系统中不能直接测量的性能参数的在线数据驱动估计方法按以下步骤进行:
一、确定工业过程控制系统中可直接测量的物理量,记为x1,x2,……,xk,……,xm;m为可直接测量的物理量的数目;
二、在不能直接测量的性能参数y的一个退化周期内,运行工业过程系统,按时间间隔为Δt分别离线测量不能直接测量的性能参数y在t=0,t1,t2,……,tk,……,tn时刻的值y0,y1,y2,……,yk,……,yn,同时分别记录t1,t2,……,tk,……,tn时刻各物理量x1,x2,……,xk,……,xm的值,得到
x1,0′, x1,1′, x1,2′,……, x1,k′,……, x1,n′;
x2,0′, x2,1′, x2,2′……, x2,k′,……, x2,n′;
……,
xm,0′, xm,1′, xm,2′,……, xm,k′,……, xm,n′;
三、将x1,0′,x1,1′,x1,2′,……,x1,k′,……,x1,n′;
x2,0′, x2,1′, x2,2′……, x2,k′, ……, x2,n′;
……,
xm,0′, xm,1′, xm,2′, ……, xm,k′, ……, xm,n′;进行归一化处理;得到
x1,0, x1,1, x1,2,……, x1,k,……, x1,n;
x2,0, x2,1, x2,2……, x1,k,……, x2,n;
……,
xm,0, xm,1, xm,2, ……, xm,k, ……, xm,n;
四、构建函数模型:其中f(t)=ceat+b,
其中,y是不能直接测量的性能参数;
f(t)表示催化剂性能的退化趋势,根据实践经验,确定为指数函数形式;
t为不能直接测量的性能参数y在一个退化周期内的反应时间,t=0时为反应开始时;
a、b、c是常量,且a<0;
表示了所测参数与不能直接测量的性能参数间的关系,x1,x2,……,xk,……,xm为可直接测量的物理量,m为可直接测量的物理量的数目;p1,p2,……,pm分别为可直接测量的物理量x1,x2,……,xk,……,xm的系数,该系数是与各循环无关的全局参数;
五、将步骤四中的f(t)=ceat+b进行离散化处理,得到fk+1(t)=(aΔt+1)fk(t)-abΔt,再将其改写成fk+1(t)=-adfk(t)+bd,其中,ad=-(aΔt+1),bd=-abΔt,Δt为步骤二中的时间间隔;
再将 改写成:
其中yk为t=tk时刻不能直接测量的性能参数y的值;
yk+1为t=tk+1时刻不能直接测量的性能参数y的值;
六、将步骤二得到的y0,y1,y2,……,yk,……,yn及
x1,0, x1,1, x1,2,……, x1,k,……, x1,n;
x2,0, x2,1, x2,2……, x1,k,……, x2,n;
……,
xm,0, xm,1, xm,2, ……, xm,k, ……, xm,n;
代入步骤五得到的离散型函数模型中,用基本最小二乘法或偏最小二乘回归法对ad、bd、p1,p2,……,pj,……,pm进行参数估计,得到ad、bd、p1、p2、……、pj、……、pm的具体值;
七、对可测变量x1,x2,……,xk,……,xm用留一交叉检验法剔除其中存在的冗余项,得到变量记为x1′,x2′,……, xn′,n≤m;与变量x1′,x2′,……, xn′对应的权重系数记为p1′、p2′、……、pj′、……、pn′;
八、将步骤七得到的x1′,x2′,……, xn′按照与其对应的权重系数p1′、p2′、……、pj′、……、pn′的绝对值由大到小顺序进行排列,然后将第一个和第二个可直接测量的物理量进行留N交叉检验(leave-N-out交叉检验),得到均方根误差Re1,然后逐次增加1个可直接测量的物理量进行留N交叉检验(leave-N-out交叉检验),得到均方根误差Re2,Re3,……,Ren-1,将Re1,Re2,Re3,……,Ren-1进行比较找到最小的一个,将最小的均方根误差所对应的可直接测量的物理量作为最终的关键变量,即为x1″,x2″,……,xr″;r≤n;与变量x1″,x2″,……,xr″对应的权重系数记为p1″、p2″、……、pj″、……、pr″;
九、将步骤六计算得出的ad、bd及步骤八计算得出的x1″,x2″,……, xr″及其对应的系数为p1″、p2″、……、pj″、……、pr″代入中;
十、按ad=-(aΔt+1),bd=-adΔt计算出a,b,再将t=0时的y0及x1″,x2″,……,xr″所代表的可直接测量的物理量在t=0时刻的值代入中计算得到c的值,再将a、b和c代入中,得到工业过程控制系统中不能直接测量的性能参数y与运行时间的关系式。
本发明是一种工业过程中不能直接测量的性能参数的在线估计方法,如催化剂性能、元器件的老化、锅炉炉膛结渣等,该方法通过建立合理的数学模型,将不能直接测量的性能参数、在线可测变量等因素纳入考虑范围,通过将未知的指数函数转化为递归形式,去除了模型中的非线性因素,然后,采用基于回归分析的参数估计方法实现对模型中各参数的估计,并使用变量选择技术提高估计方法的预测性能。得到的工业过程控制系统中不能直接测量的性能参数y与运行时间的关系式可以方便地计算出在一个周期内y的值,从而可以对y进行合理的预测估计,该方法适于对催化剂老化性能、配件老化、炉膛结渣等问题的分析预测。
具体实施方式
具体实施方式一:本实施方式的工业过程控制系统中不能直接测量的性能参数的在线数据驱动估计方法按以下步骤进行:
一、确定工业过程控制系统中可直接测量的物理量,记为x1,x2,……,xk,……,xm;m为可直接测量的物理量的数目;
二、在不能直接测量的性能参数y的一个退化周期内,运行工业过程系统,按时间间隔为Δt分别离线测量不能直接测量的性能参数y在t=0,t1,t2,……,tk,……,tn时刻的值y0,y1,y2,……,yk,……,yn,同时分别记录t1,t2,……,tk,……,tn时刻各物理量x1,x2,……,xk,……,xm的值,得到
x1,0′, x1,1′, x1,2′,……, x1,k′,……, x1,n′;
x2,0′, x2,1′, x2,2′……, x2,k′,……, x2,n′;
……,
xm,0′, xm,1′, xm,2′,……, xm,k′,……, xm,n′;
三、将x1,0′,x1,1′,x1,2′,……,x1,k′,……,x1,n′;
x2,0′, x2,1′, x2,2′……, x2,k′, ……, x2,n′;
……,
xm,0′, xm,1′, xm,2′, ……, xm,k′, ……, xm,n′;进行归一化处理;得到
x1,0, x1,1, x1,2,……, x1,k,……, x1,n;
x2,0, x2,1, x2,2……, x1,k,……, x2,n;
……,
xm,0, xm,1, xm,2, ……, xm,k, ……, xm,n;
四、构建函数模型:其中f(t)=ceat+b,
其中,y是不能直接测量的性能参数;
f(t)表示不能直接测量的性能参数性能的退化趋势,根据实践经验,确定为指数函数形式;
t为不能直接测量的性能参数y在一个退化周期内的反应时间,t=0时为反应开始时;
a、b、c是常量,且a<0;
表示了所测参数与不能直接测量的性能参数性能间的关系,x1,x2,……,xk,……,xm为可直接测量的物理量,m为可直接测量的物理量的数目;p1,p2,……,pm分别为可直接测量的物理量x1,x2,……,xk,……,xm的系数,该系数是与各循环无关的全局参数;
五、将步骤四中的f(t)=ceat+b进行离散化处理,得到fk+1(t)=(aΔt+1)fk(t)-abΔt,再将其改写成fk+1(t)=-adfk(t)+bd,其中,ad=-(aΔt+1),bd=-abΔt,Δt为步骤二中的时间间隔;
再将 改写成: 其中yk为t=tk时刻不能直接测量的性能参数y的值;
yk+1为t=tk+1时刻不能直接测量的性能参数y的值;
六、将步骤二得到的y0,y1,y2,……,yk,……,yn及
x1,0, x1,1, x1,2,……, x1,k,……, x1,n;
x2,0, x2,1, x2,2……, x1,k,……, x2,n;
……,
xm,0, xm,1, xm,2, ……, xm,k, ……, xm,n;
代入步骤五得到的离散型函数模型中,用基本最小二乘法或偏最小二乘回归法对ad、bd、p1,p2,……,pj,……,pm进行参数估计,得到ad、bd、p1、p2、……、pj、……、pm的具体值;
七、对可测变量x1,x2,……,xk,……,xm用留一交叉检验法剔除其中存在的冗余项,得到变量记为x1′,x2′,……, xn′,n≤m;与变量x1′,x2′,……, xn′对应的权重系数记为p1′、p2′、……、pj′、……、pn′;
八、将步骤七得到的x1′,x2′,……, xn′按照与其对应的权重系数p1′、p2′、……、pj′、……、pn′的绝对值由大到小顺序进行排列,然后将第一个和第二个可直接测量的物理量进行留N交叉检验(leave-N-out),得到均方根误差Re1,然后逐次增加1个可直接测量的物理量进行leave-N-out交叉检验,得到均方根误差Re2,Re3,……,Ren-1,将Re1,Re2,Re3,……,Ren-1进行比较找到最小的一个,将最小的均方根误差所对应的可直接测量的物理量作为最终的关键变量,即为x1″,x2″,……, xr″;r≤n;与变量x1″,x2″,……, xr″对应的权重系数记为p1″、p2″、……、pj″、……、pr″;
九、将步骤六计算得出的ad、bd及步骤八计算得出的x1″,x2″,……, xr″及其对应的系数为p1″、p2″、……、pj″、……、pr″代入中;
十、按ad=-(aΔt+1),bd=-abΔt计算出a,b,再将t=0时的y0及x1″,x2″,……,xr″所代表的可直接测量的物理量在t=0时刻的值代入中计算得到c的值,再将a、b和c代入中,得到工业过程控制系统中不能直接测量的性能参数y与运行时间的关系式。
本实施方式步骤一中的可直接测量的物理量的数目m的取值范围根据实际工业过程确定;
本实施方式步骤二中的采样时间间隔Δt根据实际工业过程的反应速度确定;
本实施方式的不能直接测量的性能参数的在线估计方法,如催化剂性能、元器件的老化、锅炉炉膛结渣等,通过建立合理的数学模型,将不能直接测量的性能参数、在线可测变量等因素纳入考虑范围,通过将未知的指数函数转化为递归形式,去除了模型中的非线性因素,然后,采用基于回归分析的参数估计方法实现对模型中各参数的估计,并使用变量选择技术提高估计方法的预测性能,得到的工业过程控制系统中不能直接测量的性能参数y与运行时间t的关系式,可以方便地估计出在一个周期内y在任意时刻的值,从而可以对y进行合理的预测估计。
具体实施方式二:本实施方式与具体实施方式一不同的是步骤七中对可测变量x1,x2,……,xk,……,xm用留一交叉检验法剔除其中交叉验证相关系数大于0.9的冗余项;其它与具体实施方式一相同。
本实施方式通过用留一交叉检验(Leave-one-out)法对可测变量x1,x2,……,xk,……,xm进行相关性分析,检验其各变量间的相关性,剔除其中存在较大相关的冗余项。
具体实施方式三:本实施方式与具体实施方式一不同的是工业过程控制系统中不可直接测量的物理量是催化剂的催化性能。其它与具体实施方式一相同。
具体实施方式四:本实施方式与具体实施方式一不同的是工业过程控制系统中不可直接测量的物理量是元器件的老化、锅炉炉膛结渣性能。其它与具体实施方式一相同。
用以下试验验证本发明的有益效果:
试验一:本试验以石油炼化过程中的催化裂化环节中使用的USY型分子筛催化剂为研究对象, USY型分子筛催化剂的催化性能的在线数据驱动估计方法按以下步骤进行:
一、确定石油炼化过程中的催化裂化环节中可直接测量的物理量为温度、压力、原油流量、原油的各组分含量,记为x1,x2,x4,x3,x5,……,xk,……,xm;其中m=50个;
二、在USY型分子筛催化剂性能y的一个退化周期内,运行工业过程系统,分别离线测量USY型分子筛催化剂的催化性能y在t=0,t1,t2,……,tk,……,tn时刻的值y0,y1,y2,……,yk,……,yn,同时分别记录t1,t2,……,tk,……,tn时刻各物理量x1,x2,……,xk,……,xm的值,得到
x1,0′, x1,1′, x1,2′,……, x1,k′,……, x1,n′;
x2,0′, x2,1′, x2,2′……, x2,k′,……, x2,n′;
……,
xm,0′, xm,1′, xm,2′,……, xm,k′,……, xm,n′;
三、将x1,0′,x1,1′,x1,2′,……,x1,k′,……,x1,n′;
x2,0′, x2,1′, x2,2′……, x2,k′, ……, x2,n′;
……,
xm,0′, xm,1′, xm,2′, ……, xm,k′, ……, xm,n′;进行归一化处理;得到
x1,0, x1,1, x1,2,……, x1,k,……, x1,n;
x2,0, x2,1, x2,2……, x1,k,……, x2,n;
……,
xm,0, xm,1, xm,2, ……, xm,k, ……, xm,n;
四、构建函数模型:
其中,y是催化剂的催化性能;
f(t)表示催化剂催化性能的退化趋势,确定为指数函数形式;
t为催化剂的催化性能y在一个退化周期内的反应时间,t=0时为反应开始时;
a,b,c是常量;
p1,...,pm为权重系数,权重系数是与各循环无关的全局参数;
x1,x2,……,xk,……,xm为可直接测量的物理量;
五、将步骤四所述的函数模型改成离散型:
其中yk为t=tk时催化剂的性能参数y的值;
yk+1为t=tk+1时催化剂的性能参数y的值;
六、将步骤二得到的y0,y1,y2,……,yk,……,yn及
x1,0, x1,1, x1,2,……, x1,k,……, x1,n;
x2,0, x2,1, x2,2……, x1,k,……, x2,n;
……,
xm,0, xm,1, xm,2, ……, xm,k, ……, xm,n;
代入步骤五得到的离散型函数模型中,用基本最小二乘法或偏最小二乘回归法对ad、bd、p1,p2,……,pj,……,pm进行参数估计,得到ad、bd、p1、p2、……、pj、……、pm的具体值;
七、对可测变量x1,x2,……,xk,……,xm用留一交叉检验法剔除交叉验证相关系数大于0.9的冗余项,得到变量记为x1′,x2′,……,xn′,n≤m;与变量x1′,x2′,……,xn′对应的权重系数记为p1′、p2′、……、pj′、……、pn′;
八、将步骤七得到的x1′,x2′,……, xn′按照与其对应的p1,p2,……,pj,……,pn的绝对值由大到小顺序进行排列,然后将第一个和第二个可直接测量的物理量进行leave-N-out交叉检验,得到均方根误差Re1,然后逐次增加1个可直接测量的物理量进行leave-N-out交叉检验,得到均方根误差Re2,Re3,……,Ren-1,将Re1,Re2,Re3,……,Ren-1进行比较找到最小的一个,将最小的均方根误差所对应的可直接测量的物理量作为最终的关键变量,即为x1″,x2″,……, xr″;r≤n;
九、将步骤六计算得出的ad、bd及步骤八计算得出的x1″,x2″,……, xr″及其对应的系数为p1,p2,……,pj,……,pr代入
十、按ad=-(aΔt+1),bd=-abΔt计算出a,b,再将t=0时的y0及x1″,x2″,……,xr″所代表的可直接测量的物理量在t=0时刻的值代入中计算得到c的值,再将a、b和c代入中,得到催化剂的催化性能y与运行时间的关系式。
本试验是催化剂的催化性能的在线估计方法,该方法通过建立合理的数学模型,将不能直接测量的催化剂的催化性能、在线可测变量的温度、压力、各组分的浓度等因素纳入考虑范围,通过将未知的指数函数转化为递归形式,去除了模型中的非线性因素,然后,采用基于回归分析的参数估计方法实现对模型中各参数的估计,并使用变量选择技术提高估计方法的预测性能,得到的催化剂的催化性能y与运行时间的关系式可以方便地计算出在一个周期内y的值,从而可以对y进行合理的预测估计。
Claims (3)
1.工业过程控制系统中不能直接测量的性能参数的在线数据驱动估计方法,其特征在于该方法按以下步骤进行:
一、确定工业过程控制系统中可直接测量的物理量,记为x1,x2,……,xk,……,xm;m为可直接测量的物理量的数目;
二、在不能直接测量的性能参数y的一个退化周期内,运行工业过程系统,按采样时间间隔为△t分别离线测量不能直接测量的性能参数y在t=0,t1,t2,……,tk,……,tn时刻的值y0,y1,y2,……,yk,……,yn,同时分别记录t1,t2,……,tk,……,tn时刻各物理量x1,x2,……,xk,……,xm的值,得到
x1,0′,x1,1′,x1,2′,……,x1,k′,……,x1,n′;
x2,0′,x2,1′,x2,2′……,x2,k′,……,x2,n′;
……,
xm,0′,xm,1′,xm,2′,……,xm,k′,……,xm,n′;
三、将x1,0',x1,1',x1,2',……,x1,k',……,x1,n';
x2,0′,x2,1′,x2,2′……,x2,k′,……,x2,n′;
……,
xm,0′,xm,1′,xm,2′,……,xm,k′,……,xm,n′;
进行归一化处理;得到
x1,0,x1,1,x1,2,……,x1,k,……,x1,n;
x2,0,x2,1,x2,2……,x1,k,……,x2,n;
……,
xm,0,xm,1,xm,2,……,xm,k,……,xm,n;
四、构建函数模型:其中f(t)=ceat+b,
其中,y是不能直接测量的性能参数;
t为不能直接测量的性能参数y在一个退化周期内的反应时间,t=0时为反应开始时;
a、b、c是常量,且a<0;
p1,p2,……,pm分别为可直接测量的物理量x1,x2,……,xk,……,xm的权重系数,该系数是与各循环无关的全局参数;
五、将步骤四中的f(t)=ceat+b进行离散化处理,得到fk+1(t)=(a△t+1)fk(t)-ab△t,再将其改写成fk+1(t)=-adfk(t)+bd,其中,ad=-(a△t+1),bd=-ab△t,△t为步骤二中的时间间隔;
再将 改写成:
其中yk为t=tk时刻不能直接测量的性能参数y的值;
yk+1为t=tk+1时刻不能直接测量的性能参数y的值;
六、将步骤二得到的y0,y1,y2,……,yk,……,yn及
x1,0,x1,1,x1,2,……,x1,k,……,x1,n;
x2,0,x2,1,x2,2……,x1,k,……,x2,n;
……,
xm,0,xm,1,xm,2,……,xm,k,……,xm,n;
代入步骤五得到的离散型函数模型中,用基本最小二乘法或偏最小二乘回归法对ad、bd、p1,p2,……,pj,……,pm进行参数估计,得到ad、bd、p1、p2、……、pj、……、pm的具体值;
七、对可测变量x1,x2,……,xk,……,xm用留一交叉检验法剔除其中存在的冗余项,得到变量记为x1',x2',……,xn',n≤m;与变量x1',x2',……,xn'对应的权重系数记为p1'、p2'、……、pj'、……、pn';对可测变量x1,x2,……,xk,……,xm用留一交叉检验法剔除其中交叉验证相关系数大于0.9的冗余项;
八、将步骤七得到的x1',x2',……,xn'按照与其对应的权重系数p1'、p2'、……、pj'、……、pn'的绝对值由大到小顺序进行排列,然后将第一个和第二个可直接测量的物理量进行留N交叉检验,得到均方根误差Re1,然后逐次增加1个可直接测量的物理量进行留N交叉检验,得到均方根误差Re2,Re3,……,Ren-1,将Re1,Re2,Re3,……,Ren-1进行比较找到最小的一个,将最小的均方根误差所对应的可直接测量的物理量作为最终的关键变量,即为x1″,x2″,……,xr″;r≤n;与变量x1″,x2″,……,xr″对应的权重系数记为p1″、p2″、……、pj″、……、pr″;
九、将步骤六计算得出的ad、bd及步骤八计算得出的x1″,x2″,……,xr″及其对应的系数为p1″、p2″、……、pj″、……、pr″代入中;
十、按ad=-(a△t+1),bd=-ab△t计算出a,b,再将t=0时的y0及x1″,x2″,……,xr″所代表的可直接测量的物理量在t=0时刻的值代入中计算得到c的值,再将a、b和c代入中,得到工业过程控制系统中不能直接测量的性能参数y与运行时间的关系式。
2.根据权利要求1所述的工业过程控制系统中不能直接测量的性能参数的在线数据驱动估计方法,其特征在于工业过程控制系统中不可直接测量的物理量是催化剂的催化性能。
3.根据权利要求1所述的工业过程控制系统中不能直接测量的性能参数的在线数据驱动估计方法,其特征在于工业过程控制系统中不可直接测量的物理量是元器件的老化或锅炉炉膛结渣性能。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210431022.4A CN102914970B (zh) | 2012-11-01 | 2012-11-01 | 工业过程控制系统中不能直接测量的性能参数的在线数据驱动估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210431022.4A CN102914970B (zh) | 2012-11-01 | 2012-11-01 | 工业过程控制系统中不能直接测量的性能参数的在线数据驱动估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102914970A CN102914970A (zh) | 2013-02-06 |
CN102914970B true CN102914970B (zh) | 2014-12-03 |
Family
ID=47613383
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210431022.4A Expired - Fee Related CN102914970B (zh) | 2012-11-01 | 2012-11-01 | 工业过程控制系统中不能直接测量的性能参数的在线数据驱动估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102914970B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9791852B2 (en) * | 2014-08-21 | 2017-10-17 | General Electric Technology Gmbh | Apparatus and method for controlling at least one operational parameter of a plant |
CN108638060B (zh) * | 2018-05-03 | 2021-09-28 | 大连理工大学 | 多自由度机器人参数标定中冗余参数分析剔除方法 |
CN111158351B (zh) * | 2020-01-19 | 2021-06-01 | 哈尔滨工业大学 | 数据驱动的故障诊断与最优控制系统一体化设计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1950766A (zh) * | 2004-05-04 | 2007-04-18 | 费舍-柔斯芒特系统股份有限公司 | 加工厂中过程模块和专家系统的集成 |
CN101477112A (zh) * | 2009-01-12 | 2009-07-08 | 浙江大学 | 工业流化床气相聚乙烯装置的产品质量在线软测量方法 |
CN101598737A (zh) * | 2009-06-26 | 2009-12-09 | 华东理工大学 | 精对苯二甲酸中对羧基苯甲醛含量的软测量方法 |
CN101825723A (zh) * | 2004-03-17 | 2010-09-08 | 维斯特恩格科地震控股有限公司 | 海上地震测量方法和系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7894989B2 (en) * | 2005-06-09 | 2011-02-22 | Exxonmobil Upstream Research Co. | Method for determining earth vertical electrical anisotropy in marine electromagnetic surveys |
-
2012
- 2012-11-01 CN CN201210431022.4A patent/CN102914970B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101825723A (zh) * | 2004-03-17 | 2010-09-08 | 维斯特恩格科地震控股有限公司 | 海上地震测量方法和系统 |
CN1950766A (zh) * | 2004-05-04 | 2007-04-18 | 费舍-柔斯芒特系统股份有限公司 | 加工厂中过程模块和专家系统的集成 |
CN101477112A (zh) * | 2009-01-12 | 2009-07-08 | 浙江大学 | 工业流化床气相聚乙烯装置的产品质量在线软测量方法 |
CN101598737A (zh) * | 2009-06-26 | 2009-12-09 | 华东理工大学 | 精对苯二甲酸中对羧基苯甲醛含量的软测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102914970A (zh) | 2013-02-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101899563B (zh) | 基于pca模型的连续退火机组炉内温度、张力监测及故障追溯方法 | |
CA2679128C (en) | Method and system of using inferential measurements for abnormal event detection in continuous industrial processes | |
CN109884892B (zh) | 基于交叉相关时滞灰色关联分析的流程工业系统预测模型 | |
CN107533560A (zh) | 用于推断原料成分的数据清理系统和方法 | |
CN107784178B (zh) | 一种基于多故障机理耦合的机械结构可靠性分析方法 | |
CN101211383A (zh) | 一种高炉铁水硅含量的特征分析预报方法 | |
CN104714537A (zh) | 一种基于联合相对变化分析和自回归模型的故障预测方法 | |
CN110569566B (zh) | 一种板带力学性能预测方法 | |
CN102914970B (zh) | 工业过程控制系统中不能直接测量的性能参数的在线数据驱动估计方法 | |
CN106156434A (zh) | 基于局部时滞重构的滑动窗时间差‑高斯过程回归建模方法 | |
CN106227910A (zh) | 一种基于灰色系统理论的加速退化试验可靠性评估方法 | |
CN110957011A (zh) | 连续搅拌反应器在未知时变测量噪声下的在线生产参数估计方法 | |
CN101833330A (zh) | 基于无激励闭环辨识的控制性能测试方法 | |
CN110942258B (zh) | 一种性能驱动的工业过程异常监测方法 | |
CN110222825B (zh) | 一种水泥成品比表面积预测方法及系统 | |
CN106326525A (zh) | 基于广义可加模型的热轧变形抗力预报方法 | |
CN111637925A (zh) | 一种桥梁伸缩缝状态的预警方法和预警系统 | |
CN103028614A (zh) | 一种带钢热连轧生产过程控制系统的优化方法 | |
Cong et al. | Multirate partial least squares for process monitoring | |
CN107679659A (zh) | 新型工业用电量预测方法 | |
CN114527731B (zh) | 基于监督概率慢特征分析的工业过程运行状态评价方法 | |
US10456704B2 (en) | Method of real-time prognosis of flooding phenomenon in packed column | |
Lee et al. | Adaptive multivariate regression modeling based on model performance assessment | |
JP2002108412A (ja) | モデル構築方法およびモデル構築システム | |
CN114819743B (zh) | 一种化工企业能耗诊断分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141203 Termination date: 20211101 |
|
CF01 | Termination of patent right due to non-payment of annual fee |