一种气动调节阀粘滞特性参数的估计方法与流程

本发明属于工业过程控制领域，具体涉及一种气动调节阀粘滞特性参数的在线估计方法。

背景技术：

典型的工业过程控制领域的单回路控制系统如图1所示。气动调节阀是工业过程系统中常见的执行设备之一，其性能的好坏直接决定了工业过程控制系统的性能，对产品质量、生产装置运行的经济性和稳定性有着直接影响。由于受到支撑件损坏、阀芯及阀座受腐蚀、弹簧膜片破损、密封性下降、阀杆封装过紧、金属过热膨胀、润滑性能下降等因素的影响，气动调节阀在实际生产过程中呈现严重的非线性特征。这些非线性特征是导致控制系统性能恶化和回路振荡的主要原因之一。某一回路的振荡往往会波及其它控制回路，导致生产装置整体性能下降、原材料消耗过大、能耗增加、执行器磨损加速，甚至会导致系统不稳定。因此，通常需要对气动调节阀的输入输出特性进行检测，确定其非线性特征，然后设计相应的补偿措施，使其输入输出特征更为线性化。这样一方面有利于系统控制性能的提升，另一方面也有利于生产过程运行的稳定性和经济性。

气动调节阀中的非线性特征一般可以用死区、滞环等非线性环节描述，如图2所示。其突出的特点是调节阀上升行程和下降行程不重合，上下行程存在变差，导致控制器发出的控制指令不能得到有效执行，进而影响控制系统性能。

常见的气动执行器非线性特征的检测方法有交叉互相关方法(horch,2000)、定性hammerstein模型方法(srinivasan,2005)、面积比率法(singhal,2005)、基于继电器方法(rossi,2005)、模式识别方法(yamashita,2005)、双相关方法(choudhury,2006)、hammerstein模型方法(choudhury,2008)、全局收索方法(jelali,2008)等。上述方法利用控制输入信号和被控过程输出信号对阀门非线性特征进行检测，通常需要在控制回路出现振荡的情况下才有效，在回路没有出现振荡时其有效性将大大降低或者根本无法使用。

技术实现要素：

本发明针对现有技术的不足，给出一种气动调节阀非线性特性的在线检测方法，采用估计器对调节阀非线性特征进行实时估计，为故障诊断以及控制系统的非线性补偿提供有效信息。

一种气动调节阀粘滞特性参数的检测方法，具体为：将气动调节阀的粘滞非线性特征通过线性化的方式转化为已知输入和未知输入两个部分，借助于未知输入估计器利用现场数据对未知输入进行在线估计；再在现有控制信号的基础上叠加测试信号以获取估计信息，进而利用估计信息对调节阀的非线性特征进行拟合。

作为优选：将气动阀的粘滞非线性特征描述为未知非线性函数；

作为优选：在现有控制信号的基础上叠加的测试信号为方波、prbs或gbn信号。

作为优选：利用估计器估计的未知输入值和原有的已知输入值绘制散点图，从图中直接读取粘滞特性的基本参数。

作为优选：一种气动调节阀粘滞特性参数的检测方法，该方法具体包括以下步骤：

步骤1：建立被控过程的数学模型。

被控过程的数学模型可以用如下形式的线性离散状态空间模型来描述：

其中x为系统状态，n为气动调节阀开度，y为被控过程中的被控工艺参数的测量值，a、b、c为系统矩阵，w和v为过程噪声和测量噪声，一般假设w和v为零均值的白噪声，且方差为和上述符号中的下标k表示采样时刻。

步骤2:将气动调节阀的非线性特征转化为输入不确定性。

实际系统中气动调节阀由于各种原因导致其输出的开度n同控制信号u不一致，二者之间存在一定的非线性关系。为描述这一非线性特征，可以在线性模型(1)中用非线性函数n＝f(u)，这样模型(1)可以表示为：

其中f(u)为一非线性函数，用以描述气动调节阀的非线性特征，其余参数意义同模型(1)相同。由于描述气动阀非线性函数f(u)通常是未知的，而且在生产装置实际运行中其特性还有可能发生变化，因此无法直接确定其非线性函数的具体表达式。为解决这一问题，在k时刻，将非线性函数f(uk)在任一工作点u0处进行taylor展开，可得

其中是非线性环节的展开的一阶项，εk为其余高阶项之和，引入未知输入项dk＝(α-i)uk+ok，由式(3)可得：

f(uk)＝uk+dk(4)

这样，过程模型(2)可以表示为：

通过上述步骤将气动调节阀的非线性特征f(u)转化为未知输入d来描述。

步骤3：设计未知输入估计器对输入不确定性进行估计。

在将非线性特征f(u)转换为输入不确定d后，需要对d的值进行在线估计。接下来利用未知输入估计器对模型中的状态x和未知输入d进行在线估计，其设计方法如下：

其中为利用k-1时刻对k时刻的状态进行预测的估计值；为k-1时刻对d的估计值；为状态估计过程中的中间变量；为k时刻对状态x的估计值；mk和kk为估计器在k时刻的增益矩阵。mk和kk的计算由下式决定：

其中κk为任选矩阵，但必须保证满秩。

步骤4：在控制器输出端叠加测试信号ud。

由于控制器的控制量u在通常情况下变化信息不够丰富，并且在闭环控制作用下由于控制器的调节作用，其变动形式和幅值都非常有限，利用这些信息并不一定能准确分析出f(u)具体的非线性形式。为了进一步充分激发非线性特征f(u)中的动态信息，保证估计结果的全面性和合理性，本发明考虑在控制量输入端叠加测试信号ud，如方波、prbs或gbn信号。这里，信号ud的大小和形式都是可以离线确定的，即ud为已知信号。

步骤5：气动调节阀非线性特性的估计。

通过前三个步骤的分析可以看出，气动调节阀的输入为u，输出为f(u)，而通过式(5)可以看出，f(u)的值可以通过当前的控制作用u和未知输入d的估计值进行近似：

其中ud,k为测试信号ud在k时刻的值。这样气动调节阀的非线性特征f(u)可以利用u+ud和进行分析和确定。

本发明的主要优点在于不依赖于被控过程是否处于振荡状态以及气动调节阀是否存在阀位反馈信号n，在一般工况下仍然能准确获取调节阀的非线性特征。本发明可以有效解决气动阀非线性特征的检测问题，可以为调节阀非线性补偿提供所需信息，检测结果可以用于系统的性能分析和故障诊断，对提升系统运行的安全性和准确性有显著作用。

附图说明

图1典型工业过程控制回路示意图；

图2气动调节阀非线性特征描述示意图；

图3液位控制系统工艺流程图；

图4液位控制系统方框图；

图5测试信号叠加示意图；

图6粘滞参数检测结果图。

具体实施方式

为了验证本发明的有效性，采用某生产装置的液位系统实例证明。本液位控制系统的工艺流程如图3所示，工艺液体通过气动调节阀1进入储液罐，然后经气动调节阀2流向下一生产工序。为保证后续生产工艺的正常运转，需要保证储液罐的液位稳定在一定的设定值。气动调节阀1用于实现对液位的控制，调节阀2保持开度不变，保证下一工序环节的工艺液体供应，整个控制系统的框图如图4所示。以下按照本发明的实时步骤给出简要介绍。

步骤1：系统数学模型的建立。

假设通过调节阀1的瞬时流量增量为δq1，通过调节阀2的瞬时流量增量为δq2，储液罐的截面积为a，调节阀1的流量系数为ku，调节阀2的阻力系数为rs，液位控制器的调节增量为δu。

在工作点(u0,h0)附近，根据物料平衡原理有：

根据调节阀1的开度与流量关系可得：

δq1＝kuδu(9)

根据流体力学基本原理有：

由式(8)、(9)和(10)可得：

对式(11)进行拉普拉斯变换可得：

其中s为拉普拉斯算子，k＝kurs，t＝ars。模型中的参数a、ku和rs可以参考储液罐规格、调节阀的相关特性参数确定。

步骤2：将模型(12)转换为带未知输入的线性状态空间方程。

式(12)为传递函数模型，可以转换为状态空间形式，再经过离散化可得如下离散线性状态空间模型：

其中噪声w和v的统计特征参数q和r一般可以根据实际系统情况进行设定，传感器获得的测量值yk＝y0+δyk，调节阀1的实际输出值为uk＝u0+δuk。在当前测试系统中a＝0.8187，b＝0.4532，c＝0.48，r＝10^-5，q＝10^-5。

步骤3：在控制量δu端叠加方波信号、prbs或者gbn信号ud。

在控制量δu端叠加方波、prbs或者gbn信号ud并记如图5所示，这样模型(14)变为：

步骤4：根据未知输入估计器对d进行在线估计。

根据当前测量值hk计算当前液位增量δhk＝hk-h0。根据δhk和利用式(6)和(7)对d的值进行估计，可得依次迭代，并将和记录并保存。为使用方便，这里用un和dn表示和的n组历史数据：

步骤5：绘制un和dn的散点图，如图6所示，从而可以确定该气动调节阀的变差约为0.2。