2

Course | SBM 工控课程汇报 PPT

 2 years ago
source link: https://zjcqn.github.io/posts/SBM/
Go to the source link to view the article. You can view the picture content, updated content and better typesetting reading experience. If the link is broken, please click the button below to view the snapshot at that time.

Course | SBM 工控课程汇报 PPT

Posted 1 year ago2020-10-10T12:33:00+08:00 by 陈钱牛
Updated 1 year ago2021-01-23T19:25:32+08:00

工业控制系统安全至关重要,在事故表现出人能感知的征兆之前对异常进行预测是减小事故损失的有效方法之一。

SBM(Simularity-based modeling,相似性建模)算法能够有效实现正常工况估计这一功能,并在浙能集团的生产中得到了大规模的应用。课程要求使用SBM来实现正常工况的预测和异常检测功能。

  • 数据集:训练集包括2461个时间点,时间点间隔为1小时,每个时间点包括10个测点数据,这10个测点经过专家筛选。测试集格式和训练集一样,包含1000个时间点,包含若干个异常点,每处时间长度50单位左右,均为单一参数故障。
  • 输出:每一个时间点XiXi的正常工况估计值XiestXesti,相对误差,以及相关曲线,异常点位置。
  • 评分标准:共三个得分项,
    • N处异常,识别出1处位置加5分,判定出异常参数加5分。
    • 测试数据中正常工况数据预估值与真实值相对误差绝对值的平均数,0.05-0.01得10分,0.01以下,每降低0.001加5分。
    • 状态矩阵向量数,100及以上得10分,100以下每减少10条加5分。

SBM基本算法

  1. 构建D矩阵
D=[X(1),X(2),⋯,X(m)]=⎡⎣⎢⎢⎢⎢X1(1)X2(1)⋯Xn(1)X1(2)X2(2)⋯Xn(2)⋯⋯⋯⋯X1(m)X2(m)⋯Xn(m)⎤⎦⎥⎥⎥⎥n×mD=[X(1),X(2),⋯,X(m)]=[X1(1)X1(2)⋯X1(m)X2(1)X2(2)⋯X2(m)⋯⋯⋯⋯Xn(1)Xn(2)⋯Xn(m)]n×m
  1. 输入真实值
Xin=[X1in,X2in,⋯,Xnin]TXin=[Xin1,Xin2,⋯,Xinn]T
  1. 产生估计值
aω0GωXest=DT⋅Xin=G−1⋅a=DT⊗D=ω0∑mi=0ωi0=D⋅ωa=DT⋅Xinω0=G−1⋅aG=DT⊗Dω=ω0∑i=0mω0iXest=D⋅ω

Ps: ⊗⊗ is a distance function

  1. 判断是否工作正常
⎧⎩⎨      ∣∣∣Xest−XX∣∣∣>σothers,alarm,no alarm{|Xest−XX|>σ,alarmothers,noalarm

pic1


pic2


pic3


pic4


pic5


pic6


pic7


pic8


pic9


pic10


pic11


This post is licensed under CC BY 4.0 by the author.
Share

About Joyk


Aggregate valuable and interesting links.
Joyk means Joy of geeK