CN113866676A - 一种基于卡尔曼滤波的vsc微小故障诊断方法 - Google Patents

一种基于卡尔曼滤波的vsc微小故障诊断方法 Download PDF

Info

Publication number
CN113866676A
CN113866676A CN202111130636.4A CN202111130636A CN113866676A CN 113866676 A CN113866676 A CN 113866676A CN 202111130636 A CN202111130636 A CN 202111130636A CN 113866676 A CN113866676 A CN 113866676A
Authority
CN
China
Prior art keywords
model
vsc
current
fault
kalman filtering
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.)
Granted
Application number
CN202111130636.4A
Other languages
English (en)
Other versions
CN113866676B (zh
Inventor
郭楚佳
张莹
刘伟峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yunnan Aerospace Automation Technology Co ltd
Original Assignee
Shaanxi University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shaanxi University of Science and Technology filed Critical Shaanxi University of Science and Technology
Priority to CN202111130636.4A priority Critical patent/CN113866676B/zh
Publication of CN113866676A publication Critical patent/CN113866676A/zh
Application granted granted Critical
Publication of CN113866676B publication Critical patent/CN113866676B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/50Testing of electric apparatus, lines, cables or components for short-circuits, continuity, leakage current or incorrect line connections
    • G01R31/52Testing for short-circuits, leakage current or ground faults

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Testing Electric Properties And Detecting Electric Faults (AREA)

Abstract

本发明公开了一种基于卡尔曼滤波的VSC微小故障诊断方法,首先建立了VSC系统的离散预测模型,其次基于卡尔曼滤波算法建立系统模型的状态方程和量测方程,同时推导系统模型的卡尔曼滤波增益以及状态协方差,再根据对系统电流测量的预测值与系统实际电流的测量值比较,得到电流残差。对n个时间点的残差值进行加权平方和,计算得出故障指示信号,最后设定决策函数,将故障指示信号与所设定阈值进行比较判断系统模型是否故障。相对于常规的故障诊断方法,基于卡尔曼滤波算法对系统的内部结构有深层的认识,同时依赖于系统精确模型,将具有实时性,从而通过微小故障诊断增强系统的安全性与可靠性。

Description

一种基于卡尔曼滤波的VSC微小故障诊断方法
技术领域
本发明属于电力系统技术领域,具体涉及一种基于卡尔曼滤波的VSC微小故障诊断方法。
背景技术
近几年来,随着电力电子技术的成熟和电力电子器件的不断发展,两电平电压源型换流器(Voltage Source Converter,VSC)在电力系统中的应用越来越广泛。随着现代控制系统设备的复杂化和规模的大型化,未被检测出的微小故障可能会逐渐演变为对系统有较大危害的突变故障,一旦系统发生故障且没有及时发现和处理,故障很有可能发酵,造成更大的经济损失甚至人员伤亡。随着系统运行。因此,为了增强系统的安全性和可靠性,微小故障的检测和诊断尤为重要。
故障是指在一个过程中,观测变量或计算参数对可接受范围的偏离。而微小故障是一类更精细的诊断形式,针对故障诊断提出的在线近似法主要是通过在线自适应估计器对故障进行估计,估计结果可用于进行故障分离和故障辨识,但是这种方法故障估计的准确性现在还尚未能够在理论证明;人工智能定量诊断微小故障法是一种依赖于定量分析的智能诊断方法,其中包括的神经网络具有强的非线性拟合能力,自学习能力、自组织能力和容错能力,但其建模对偏差很敏感,任何微小的偏差或扰动都会导致拟合精度下降影响诊断结果,并且神经网络对诊断结果的可解释性相对较差;状态估计法对系统的内部结构有深层的认识,具有很好的实时性,但是其需要依赖于系统准确的数学模型。由于系统模型不可避免的存在不确定性和干扰,直接影响现有的电力电子器件微小故障诊断方法的准确性。
发明内容
为了克服上述现有技术存在的不足,本发明的目的在于提供一种基于卡尔曼滤波的VSC微小故障诊断方法,建立VSC的系统模型,得到系统的状态方程和量测方程,基于卡尔曼滤波算法对VSC系统的方程滤波,得到当前时刻的状态估计值,再与系统测量值进行比较得到残差,取n个时间点残差的加权平方和生成故障指示信号,设计决策函数,利用故障指示信号进行故障诊断。
为了实现上述目的,本发明采用的技术方案是:
一种基于卡尔曼滤波的VSC微小故障诊断方法,包括以下步骤;
Step1:建立主电路模型;
根据VSC的结构、电路结构及KVL定理,可得:
Figure BDA0003280351500000021
Step2:利用前向欧拉公式,得到电流的预测模型;
Figure BDA0003280351500000022
式中i是VSC模型交流测滤波电感电流、E是交流侧电压、u是直流侧电压、L是滤波电感,Ts是采样间隔;
Step3:基于卡尔曼滤波算法和模型方程建立VSC的状态方程和量测方程;
Figure BDA0003280351500000031
Figure BDA0003280351500000032
vf(k)和wf(k)是相互独立的模型噪声矩阵,式中H为单位阵;
Step4:计算系统模型当前时刻的电流残差;
由卡尔曼滤波算法得VSC系统模型电流的状态预测、量测预测方程:
Figure BDA0003280351500000033
由模型得测量残差:
Figure BDA0003280351500000034
由卡尔曼滤波算法得卡尔曼增益:
W(k+1)=P(k+1|k)R-1(k+1)
式中P(k+1|k)是VSC系统模型协方差:P(k+1|k)=P(k|k)+Q(k)
式中Q(k)和R(k)是VSC系统方程定义的误差统计量;
卡尔曼滤波算法得模型预测电流:
Figure BDA0003280351500000035
上式是在下一刻对下一刻模型电流的测量预测值,即在当前时刻对模型电流测量预测值:
Figure BDA0003280351500000041
VSC系统模型中已知实际的电流测量值i(k),卡尔曼滤波算法得到系统电流在当前时刻的测量预测值
Figure BDA0003280351500000042
分析实际测量值和预测测量值,其残差为:
Figure BDA0003280351500000043
Step5:设计决策函数;
计算一段时间内n个时间点残差的加权平方和:
Figure BDA0003280351500000044
式中Xm是m时刻的故障指示信号,ωm是m时刻的权重值,Rm是m时刻噪声矩阵;
设定故障指示信号为:
Figure BDA0003280351500000045
故障指示信号与决策函数比较,判断VSC系统模型是否出现故障;
所设置决策函数为
Figure BDA0003280351500000046
其中阈值e的选取需要在具体模型数据下,通过分析诊断电流在故障模式下和健康模式下的特征选取;
Step6:判断系统是否故障;
如故障指示信号大于所设定阈值e则系统出现故障,相反则系统无故障发生。
本发明的有益效果:
本发明针对卡尔曼滤波的VSC微小故障诊断,选用两电平电压源型换流器,得到系统模型状态方程和量测方程,基于卡尔曼滤波算法对模型计算,提高系统稳定性;由于故障元器件参数值突然发生变化,导致电路功能丧失而造成开路或短路等故障,其中有些微小故障不容判断但存在大隐患,采用卡尔曼率滤波算法得VSC模型电流的估计值,与实际值比较后的残差进行加权平方和,得到故障指示信号。将故障指示信号与所设定阈值比较,判断VSC系统微小故障,增强系统的安全性和可靠性。
附图说明
图1两电平电压源换流器原理图。
图2是本发明的算法框图。
图3是本发明的流程框图。
具体实施方式
下面结合附图对本发明作进一步详细说明。
本发明是针对VSC装置基于卡尔曼滤波算法实现微小故障诊断,以图1所示的结构及标号为例,针对三相两电平VSC微小故障诊断研究。详细说明本发明的实施方法。
具体检测步骤如下:
Step1:建立主电路模型;
本发明所涉及的VSC结构,以及接入系统时的原理图如图1所示,根据电路结构及KVL定理,可得:
Figure BDA0003280351500000051
Step2:利用前向欧拉公式,得到电流的预测模型;
Figure BDA0003280351500000061
Step3:基于卡尔曼滤波算法和模型方程建立VSC的状态方程和量测方程;
Figure BDA0003280351500000062
Figure BDA0003280351500000063
式中vf(k)和wf(k)是相互独立的模型的噪声矩阵,H为单位矩阵。
Step4:计算系统模型当前时刻的电流残差;
基于卡尔曼滤波算法,得到VSC系统模型电流的状态预测和量测预测:
Figure BDA0003280351500000064
由卡尔曼滤波算法得,在当前时刻模型的电流测量预测:
Figure BDA0003280351500000065
式中W(k)为卡尔曼增益,表示为:
W(k+1)=P(k+1|k)R-1(k+1)
式中P(k+1|k)是VSC系统模型的协方差:P(k+1|k)=P(k|k)+Q(k)
其中Q(k)和R(k)为VSC系统方程定义的误差统计量。
VSC系统模型已知实际的电流测量值i(k),卡尔曼滤波算法得到系统电流的测量预测值
Figure BDA0003280351500000066
即残差为实际测量值和预测测量值之差:
Figure BDA0003280351500000071
Step5:设计决策函数;
计算n个时间点残差的加权平方和:
Figure BDA0003280351500000072
式中Xm是m时刻的故障指示信号,ωm是m时刻的权重值,Rm是m时刻的噪声矩阵;
故障指示信号为:
Figure BDA0003280351500000073
设计决策函数为
Figure BDA0003280351500000074
其中e为所设定阈值。
Step6:判断系统是否故障;
根据Step5中所设计的决策函数,当故障指示信号大于所设定阈值e则系统出现故障,相反则系统无故障发生。
本发明针对电力系统中两电平电压源型换流器,由于元器件导致的开路、短路等情况进行故障诊断,采用卡尔曼滤波算法对VSC模型电流估计值与实际值之间的残差进行加权平方和,得到故障指示信号,并与所设定阈值进行比较,判断VSC系统微小故障,增强系统的安全性与可靠性。对于VSC微小故障面临的小偏差故障诊断效果显著。
如图1所示:两电平电压源型换流器电路图,图中直流侧(右侧)并联一个足够大的电容(C),保证在持续充/放电和元器件换相过程中保持电压不会发生打的变化。交流侧(左侧)通过电感(L)与电网(E)相连,串联电感(L)是在阻抗较小情况下防止直流侧电容(C)发生短路放电,损坏器件和装置。可关断器件T1-T6的导通或者关断是通过其控制端和一个主电路端子之间施加一定的控制信号控制,防止直流侧电压源短路,同一支路的上、下桥臂不能同时导通。如图2所示为卡尔曼滤波算法框图。如图3所示为微小故障诊断控制框图。

Claims (1)

1.一种基于卡尔曼滤波的VSC微小故障诊断方法,其特征在于,包括以下步骤;
Step1:建立主电路模型;
根据VSC的结构、电路结构及KVL定理,可得:
Figure FDA0003280351490000011
Step2:利用前向欧拉公式,得到电流的预测模型;
Figure FDA0003280351490000012
式中i是VSC模型交流测滤波电感电流、E是交流侧电压、u是直流侧电压、L是滤波电感,Ts是采样间隔;
Step3:基于卡尔曼滤波算法和模型方程建立VSC的状态方程和量测方程;
Figure FDA0003280351490000013
Figure FDA0003280351490000014
vf(k)和wf(k)是相互独立的模型噪声矩阵,式中H为单位阵;
Step4:计算系统模型当前时刻的电流残差;
由卡尔曼滤波算法得VSC系统模型电流的状态预测、量测预测方程:
Figure FDA0003280351490000021
由模型得测量残差:
Figure FDA0003280351490000022
由卡尔曼滤波算法得卡尔曼增益:
W(k+1)=P(k+1|k)R-1(k+1)
式中P(k+1|k)是VSC系统模型协方差:P(k+1|k)=P(k|k)+Q(k)
式中Q(k)和R(k)是VSC系统方程定义的误差统计量;
卡尔曼滤波算法得模型预测电流:
Figure FDA0003280351490000023
上式是在下一刻对下一刻模型电流的测量预测值,即在当前时刻对模型电流测量预测值:
Figure FDA0003280351490000024
VSC系统模型中已知实际的电流测量值i(k),卡尔曼滤波算法得到系统电流在当前时刻的测量预测值
Figure FDA0003280351490000025
分析实际测量值和预测测量值,其残差为:
Figure FDA0003280351490000026
Step5:设计决策函数;
计算一段时间内n个时间点残差的加权平方和:
Figure FDA0003280351490000027
式中Xm是m时刻的故障指示信号,ωm是m时刻的权重值,Rm是m时刻噪声矩阵;
设定故障指示信号为:
Figure FDA0003280351490000031
故障指示信号与决策函数比较,判断VSC系统模型是否出现故障;
所设置决策函数为
Figure FDA0003280351490000032
其中e为所设定阈值;
Step6:判断系统是否故障;
如故障指示信号大于所设定阈值e则系统出现故障,相反则系统无故障发生。
CN202111130636.4A 2021-09-26 2021-09-26 一种基于卡尔曼滤波的vsc微小故障诊断方法 Active CN113866676B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111130636.4A CN113866676B (zh) 2021-09-26 2021-09-26 一种基于卡尔曼滤波的vsc微小故障诊断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111130636.4A CN113866676B (zh) 2021-09-26 2021-09-26 一种基于卡尔曼滤波的vsc微小故障诊断方法

Publications (2)

Publication Number Publication Date
CN113866676A true CN113866676A (zh) 2021-12-31
CN113866676B CN113866676B (zh) 2024-10-01

Family

ID=78990748

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111130636.4A Active CN113866676B (zh) 2021-09-26 2021-09-26 一种基于卡尔曼滤波的vsc微小故障诊断方法

Country Status (1)

Country Link
CN (1) CN113866676B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4812995A (en) * 1987-05-19 1989-03-14 Girgis Adly A Adaptive Kalman Filtering in fault classification
CN108872790A (zh) * 2018-07-10 2018-11-23 武汉科技大学 一种基于卡尔曼滤波和支持向量机的mmc故障诊断方法
CN109857094A (zh) * 2019-03-14 2019-06-07 杭州电子科技大学 基于两阶段卡尔曼滤波算法的航空发动机故障诊断方法
CN110749835A (zh) * 2019-10-09 2020-02-04 三峡大学 基于卡尔曼滤波器的输电线路故障定位方法
KR20210075479A (ko) * 2019-12-13 2021-06-23 한밭대학교 산학협력단 컨버터를 위한 센서 고장 진단 시스템

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4812995A (en) * 1987-05-19 1989-03-14 Girgis Adly A Adaptive Kalman Filtering in fault classification
CN108872790A (zh) * 2018-07-10 2018-11-23 武汉科技大学 一种基于卡尔曼滤波和支持向量机的mmc故障诊断方法
CN109857094A (zh) * 2019-03-14 2019-06-07 杭州电子科技大学 基于两阶段卡尔曼滤波算法的航空发动机故障诊断方法
CN110749835A (zh) * 2019-10-09 2020-02-04 三峡大学 基于卡尔曼滤波器的输电线路故障定位方法
KR20210075479A (ko) * 2019-12-13 2021-06-23 한밭대학교 산학협력단 컨버터를 위한 센서 고장 진단 시스템

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈益广,李响,傅涛: "基于推广卡尔曼滤波算法的永磁同步电机无位置传感器控制方法", 微电机, no. 05 *

Also Published As

Publication number Publication date
CN113866676B (zh) 2024-10-01

Similar Documents

Publication Publication Date Title
CN109975656B (zh) 基于柔性直流配电网高阻接地故障检测方法
CN113219300B (zh) 一种基于相电流暂态稳态的配电网单相接地故障感知方法
CN112051481B (zh) 一种基于lstm的交直流混联电网故障区域诊断方法和系统
JP2021148778A (ja) 非接地配電システムのための過渡状態に基づいた故障位置特定方法
CN107016236B (zh) 基于非线性量测方程的电网假数据注入攻击检测方法
WO2019080367A1 (zh) 一种机械设备健康状态评估方法
CN106932686B (zh) 基于最大最小特征值法的电网异常快速检测方法
CN108535572B (zh) 基于基波零序特征的计量系统二次回路监测方法及装置
CN112187528B (zh) 基于sarima的工业控制系统通信流量在线监测方法
CN110516742A (zh) 一种基于组合神经网络模型的配电终端故障判别方法及系统
CN109739214A (zh) 工业过程间歇故障的检测方法
CN107329049A (zh) 一种基于卡尔曼滤波器的输电线路早期故障检测方法
CN114217164B (zh) 一种基于放电波形智能识别的电缆故障测距方法及系统
CN112034387B (zh) 一种基于预测序列的输电线路短路故障诊断方法及装置
CN111856209A (zh) 一种输电线路故障分类方法及装置
CN106301610A (zh) 一种超外差接收机的自适应故障检测和诊断方法及装置
CN109613372B (zh) 一种基于多元电网数据库的电网故障诊断方法
CN109406923A (zh) 一种列车供电系统接地漏电故障预测方法及装置
CN109635430B (zh) 电网输电线路暂态信号监测方法和系统
US20130158903A1 (en) Method and System for Detecting Transients in Power Grids
CN117849691A (zh) 一种电容式电压互感器多维协同运行监测预警系统和方法
US20220413032A1 (en) Machine learning based method and device for disturbance classification in a power transmission line
CN108092272A (zh) 一种基于渐消卡尔曼滤波的电压稳定在线监测方法
CN113866676B (zh) 一种基于卡尔曼滤波的vsc微小故障诊断方法
CN112290887A (zh) 一种基于l-m算法改进的bp神经网络电弧检测模型

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right

Effective date of registration: 20240409

Address after: 518000 1002, Building A, Zhiyun Industrial Park, No. 13, Huaxing Road, Henglang Community, Longhua District, Shenzhen, Guangdong Province

Applicant after: Shenzhen Wanzhida Technology Co.,Ltd.

Country or region after: China

Address before: 710021 Shaanxi province Xi'an Weiyang university campus of Shaanxi University of Science and Technology

Applicant before: SHAANXI University OF SCIENCE & TECHNOLOGY

Country or region before: China

TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20240904

Address after: No. 21, Building 15, Xiaopo Industrial and Mining Electromechanical Market, Zhuangjiawan Community, Longhua Street, Zhanyi District, Qujing City, Yunnan Province, 655000

Applicant after: Yunnan Aerospace Automation Technology Co.,Ltd.

Country or region after: China

Address before: 518000 1002, Building A, Zhiyun Industrial Park, No. 13, Huaxing Road, Henglang Community, Longhua District, Shenzhen, Guangdong Province

Applicant before: Shenzhen Wanzhida Technology Co.,Ltd.

Country or region before: China

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant