[发明专利]一种面向智能电厂汽轮机轴承的闭环控制系统故障退化状态预测方法

权利要求书

1.一种面向智能电厂汽轮机轴承的闭环控制系统故障退化状态预测方法，其特征在于，包括以下步骤：

(1)数据采集：获取整个故障退化过程的数据，设汽轮机轴承的闭环控制系统中包含J个过程变量，每次采样可以得到一个1×J的向量，采样M次后得到一个故障退化过程下的二维矩阵X(M×J)；

(2)提取特征：用CVA对汽轮机轴承的闭环控制系统的故障数据进行分析，提取时序相关的部分Z，然后使用SFA对时序相关的部分Z^T进行分析，分别提取稳态特征S和时变特征具体包括如下子步骤：

(2.1)使用CVA提取时序相关的部分Z；CVA具体包括以下子步骤：

(2.1.1)为了考虑数据的时序相关性，在每个时间点t，按如下方式构造过去向量x_p(t)和未来向量x_f(t)：

x_p(t)＝[x_t-1^T,x_t-2^T,...,x_t-h^T]^T∈R^Jh (1)

x_f(t)＝[x_t^T,x_t+1^T,...,x_t+l-1^T]^T∈R^Jl (2)

其中，h、l代表滞后数，x_t表示t时刻所采集的过程变量向量；

(2.1.2)对过去向量x_p(t)和未来向量x_f(t)进行标准化，然后将不同时刻的过去向量和未来向量按如下形式排列成过去矩阵X_p和未来矩阵X_f：

X_p＝[x_p(t)^T,x_p(t+1)^T,...,x_p(t+M)^T]∈R^Jh×N (3)

X_f＝[x_f(t)^T,x_f(t+1)^T,...,x_f(t+M)^T]∈R^Jl×N (4)

其中，N＝M-h-l+1；

(2.1.3)CVA的目标是找到过程变量的线性组合，以最大化过去矩阵X_p和未来矩阵X_f之间的相关性；此优化问题可以通过下面的奇异值分解得到解决：

其中，Σ_pp，Σ_ff和Σ_pf分别是矩阵X_p和X_f的协方差矩阵和互协方差矩阵；对角阵Λ＝diag(γ₁,γ₂,...,γ_r,0,...,0)，γ₁≥...≥γ_r，γ_i表示典型相关值；

(2.1.4)求出状态空间和残差空间的投影矩阵：

其中矩阵U_k由矩阵U的前k列组成，I为单位阵；

(2.1.5)将过去矩阵X_p投影到状态空间和残差空间得到状态矩阵Z和残差矩阵E：

Z＝JX_p∈R^k×N (8)

E＝LX_p∈R^Jh×N (9)

分别代表时序相关的部分和时序不相关的部分；

(2.2)使用SFA对时序相关的部分Z^T进行分析，分别提取稳态特征S和时变特征SFA的具体步骤如下：

(2.2.1)对输入Z^T的协方差矩阵进行奇异值分解：

Z^TZ_t＝OΓO^T (10)

其中_t表示在时间上求均值；

(2.2.2)对输入Z^T进行白化处理，从而使Z^T中的变量不相关，同时为单位方差；白化后的数据为：

其中白化矩阵

(2.2.3)对A的一阶差分的协方差矩阵进行奇异值分解：

(2.2.4)求出投影矩阵W和慢特征S:

慢特征S代表稳态特征，其一阶差分代表时变特征；

(3)根据闭环控制系统的故障退化状态所经历的时间，对稳态特征S和时变特征进行划分，获得各个故障状态下的稳态特征和时变特征；

(4)建立模型：使用各个故障状态下的稳态特征和时变特征来训练对应的故障退化状态的CHMM模型，建立包含各个故障退化状态CHMM模型的故障退化状态库；

(5)在线预测：对于实时采集的数据，采用步骤(2)所述的方式提取特征，然后将提取的特征分别输入故障退化状态库中的CHMM模型中，比较其输出的概率，输出概率最大的CHMM模型所对应的故障退化状态即为当前闭环控制系统所处的故障退化状态。