[发明专利]一种配电网故障行波波头的辨识方法有效
申请号: | 201611070835.X | 申请日: | 2016-11-29 |
公开(公告)号: | CN106646121B | 公开(公告)日: | 2019-01-22 |
发明(设计)人: | 黄哲洙;王洋;冯喜强;李悦悦;刘文娟;郎斌;钟立波;马煜;张家豪;孟繁军;王彦静;刘文祥 | 申请(专利权)人: | 国网辽宁省电力有限公司沈阳供电公司;国家电网公司;北京丹华昊博电力科技有限公司 |
主分类号: | G01R31/08 | 分类号: | G01R31/08 |
代理公司: | 沈阳亚泰专利商标代理有限公司 21107 | 代理人: | 许宇来 |
地址: | 110003 辽*** | 国省代码: | 辽宁;21 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明公开了一种配电网故障行波波头的辨识方法。对故障电压行波进行凯伦贝尔变换得到故障电压线模分量,接着构建行波波头数学模型,根据该模型构造行波信号状态方程和观测方程,然后运用无迹卡尔曼滤波对故障电压行波去噪,并根据收敛的状态变量得到故障行波到达时刻。本发明能够在信噪比低的情况下依旧准确辨识行波波头,具有很好的工程应用前景。 | ||
搜索关键词: | 波头 行波 辨识 故障电压行波 配电网故障 无迹卡尔曼滤波 工程应用 故障电压 故障行波 观测方程 模型构造 数学模型 线模分量 行波信号 状态变量 状态方程 信噪比 构建 去噪 收敛 | ||
【主权项】:
1.一种配电网故障行波波头的辨识方法,其特征在于:所述方法包括以下步骤:步骤1:对故障采集器得到的三相故障电压行波进行凯伦贝尔变换,得到线模电压分量,其中,凯伦贝尔变换为:
式中u1,u2为线模电压,u0为零模电压,ua,ub,uc为各相电压;步骤2:构建行波波头数学表达式,其表达式可以由以下函数表征:
其中,u(k)为k时刻故障线模电压信号的值,T是行波到达故障检测装置的时刻,当k<T时,故障行波还没有到达,系统正常,As表示基波信号的幅值,ωc为基波信号的角频率,
为基波信号的初始相角;当k>T时,信号发生突变,在基波信号的基础上,叠加了暂态行波,其中Ae表示初始暂态行波波头的幅值,Ts为行波衰减时间常数;步骤3:根据步骤2构建的行波信号表达式,选取状态变量
构建离散化的行波信号状态方程和观测方程如下:x(k)=x(k‑1) k=1,2...u(k)=H[x(k‑1),k]+V(k) k=1,2...式中,x(K)为k时刻系统状态变量,x(k‑1)为k‑1时刻系统状态变量,H[]表示故障电压信号u(k)与状态变量的函数关系,V(K)表示观测噪声;步骤4:运用无迹卡尔曼滤波对故障行波电压进行处理,得到去噪的行波和状态变量T的收敛值;4‑1、初始化状态变量均值:![]()
其中,x0为初始值,
为初始向量均值,P0是初始协方差矩阵;4‑2、计算sigma采样点:
其中L为状态空间的维数,为6,则sigma点的个数为2L+1=13;
为第k‑1时刻的状态变量均值;Xi(k‑1)为第k‑1时刻的sigma点集;Sigma点对应的权重系数为:
其中,λ为比例系数,改变λ的大小可以调节sigma点与均值点之间的距离;4‑3、计算状态变量均值的一步估计,将k‑1时刻的sigma点分别代入状态方程,得到k时刻的sigma点集
以及其对应的均值![]()
![]()
4‑4、计算协方差矩阵的一步估计:
4‑5、获得卡尔曼增益:![]()
![]()
![]()
![]()
其中,
是观测点集,
是观测点集的均值,PU(k),U(k)是观测向量的自相关矩阵,PX(k),U(k)是观测向量与状态向量的互相关矩阵,Kk为卡尔曼增益;4‑6、更新状态估计值和协方差矩阵:![]()
式中u(k)为步骤1变换得到的故障行波线模分量在k时刻的值;完成k时刻的状态变量和协方差矩阵的更新后,返回步骤4‑2,进行k+1时刻的估算;当k=N,N为故障信号的采样点数,状态变量和协方差矩阵停止更新,此时输出变换得到的电压信号U(k)和状态变量T的值,此时的U(k)为去噪后的行波信号,T为故障行波波头到达检测装置的时刻。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于国网辽宁省电力有限公司沈阳供电公司;国家电网公司;北京丹华昊博电力科技有限公司,未经国网辽宁省电力有限公司沈阳供电公司;国家电网公司;北京丹华昊博电力科技有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201611070835.X/,转载请声明来源钻瓜专利网。