一种基于衰减指数函数模型的动态谐波相量测量方法技术

本发明专利技术公开了一种基于衰减指数函数模型的动态谐波相量测量方法。所述方法包括：选定模型参数值τ

【技术实现步骤摘要】
一种基于衰减指数函数模型的动态谐波相量测量方法


[0001]本专利技术涉及电力系统监测领域，尤其涉及动态谐波相量测量方法。

技术介绍

[0002]由于大量电力电子装置和调速电机的应用，电力系统中存在越来越多的谐波。测量谐波相量可以应用于谐波状态估计、谐波源定位、电能质量监测与控制、孤岛检测和电力系统高阻故障检测等。但由于系统的动态行为或负荷的动态特征，导致谐波可能存在动态特性，尤其当系统发生故障时，谐波幅值可能存在呈指数变化特征。当幅值指数变化的速率较快(即时间常数较小)时，这一动态特性将严重影响动态谐波相量的测量准确度，从而影响高阻故障检测和孤岛检测等高级应用的准确性。
[0003]当谐波幅值快速变化时，实现动态谐波相量高准确度测量的难点是反应谐波幅值变化的时间常数较小，导致建立的信号模型无法准确表征动态谐波向量，进而导致动态谐波向量的测量误差较大。已有的泰勒傅里叶变换方法(Platas
‑
Garza M A,Serna J A D L O.Dynamic Harmonic Analysis through Taylor
‑
Fourier Transform.IEEE Transactions on Instrumentation and Measurement,2011,60(3):804
‑
813.)和基于sinc插值函数的动态谐波相量测量方法(Chen L,Zhao W,Wang Q,et al.Dynamic harmonic synchrophasor estimator based on sinc interpolation functions.IEEE Transactions on Instrumentation&Measurement,2019,68(9):3054
‑
3065.)分别构建了基于泰勒信号模型和时域采样定理的动态谐波相量参数化模型，在谐波幅值和相位在窄带范围内变化情况下具有较高的准确度，但对于幅值呈指数形式快速变化的谐波相量而言，准确度仍然较低。

技术实现思路

[0004]为了解决已有的测量方法对谐波幅值快速变化情况下的动态谐波相量的测量准确度不高的问题，本专利技术的目的是提出一种新型动态谐波相量测量方法，该方法是对于幅值快速变化的动态谐波相量的高准确度测量方法，从而为高阻故障检测和孤岛事件检测提供实时信息。具体而言，本专利技术通过构建衰减指数函数模型实现了动态谐波相量的参数化建模，进一步实现了电网信号的近似表示，并基于此提出了动态谐波相量测量算法。该算法可保证谐波相量测量单元(harmonic phasor measurement unit,HPMU)在谐波幅值快速变化条件下，依然能准确测量动态谐波相量。
[0005]为实现实时目的，本专利技术采用以下技术方案：
[0006]一种基于衰减指数函数模型的动态谐波相量测量方法，其特征在于包括以下步骤：
[0007]1)针对谐波幅值呈指数变化情况而相应构建动态谐波相量参数化模型，而用有限项衰减指数函数对动态谐波相量进行表示；
[0008]2)根据所述动态谐波相量参数化模型实现电网信号的近似表示模型，从而可根据
这一表示结果，实现衰减指数函数特征相量的估计；
[0009]3)根据步骤2)构建的电网信号模型，基于最小二乘法实现衰减指数函数特征相量的估计。
[0010]对一个谐波幅值快速变化的电网信号，可以将其表示为
[0011][0012]其中，Re{
·
}表示取相量实部算子；H为信号中的最高谐波次数；f0是额定频率；a
h
(t)和θ
h
分别代表h次谐波分量的幅值和相位，这里假设幅值是时变的；为h次动态谐波相量。
[0013]针对幅值以指数形式快速变化的谐波相量而言，在一个很窄的时间窗内，可以用有限个衰减时间常数确定的衰减函数近似表达。如此，本专利技术通过构建以下衰减指数函数模型对p
h
(t)进行近似地表示，具体如下式所示：
[0014][0015]其中，τ
h,k
为模型中提前给定的第k个指数衰减函数的时间常数；p
k,h
表征第k个指数衰减函数的对应的相量值；K+1为模型中指数衰减函数的个数。如此，K可以被视为模型阶数。需要说明的是，这里的τ
h,k
是提前设定的参数。τ
h,k
给定后，p
k,h
可以根据实际信号序列计算得到。
[0016]以T
s
为时间间隔对s(t)进行采样。假设在整个时间窗
‑
T
w
/2≤t≤T
w
/2内，相应的采样点个数为N
w
。为了使t0＝0时刻位于时间窗的正中间，N
w
须为一个奇数。设t0＝0的两侧均各有N个采样点，则N
w
＝2N+1。如此，式(1)可以整理为：
[0017][0018]将式(3)表示成矩阵形式：
[0019][0020]其中，为一个包含信号s(t)的N
w
个采样点的列向量；和是两个矩阵，它们的每一列分别含有基函数和(分别见式(5)和式(6))的N
w
个采样点；和为两个包含p
k,h
及其共轭相量的列向量(分别见式(7)和式(8))；
*
是共轭算子。
[0021][0022][0023]p
K
＝[p
‑
K,1
,...,p
K,1
,...,p
‑
K,h
,...,p
K,h
,...,p
‑
K,H
,...,p
K,H
]T
ꢀꢀꢀꢀꢀꢀꢀ
(7)
[0024][0025]可以应用最小二乘法实现列向量p的最优估计，具体如下式所示：
[0026][0027]其中，H是埃尔米特算子；
^
表示对应的值为估计值。如此，h次谐波相量可以通过下式进行计算：
[0028][0029]本专利技术的有益效果是：
[0030](1)针对谐波幅值呈指数变化情况而相应构建的动态谐波相量参数化模型，可用有限项衰减指数函数对动态谐波相量进行表示。
[0031](2)基于提出的衰减指数函数模型，实现电网信号的近似表示，进而根据表示结果实现衰减指数函数特征相量的估计。
[0032](3)根据构建的电网信号模型，基于最小二乘法实现衰减指数函数特征相量的估计，进而可基于该特征相量计算动态谐波相量值。
附图说明
[0033]图1为本专利技术基于衰减指数函数模型的动态谐波相量测量方法流程图。
具体实施方式
[0034]参照附图，本专利技术提出了一种基于衰减指数函数模型的动态谐波相本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于衰减指数函数模型的动态谐波相量测量方法，其特征在于包括以下步骤：1)针对谐波幅值呈指数变化情况而相应构建动态谐波相量参数化模型，而用有限项衰减指数函数对动态谐波相量进行表示；2)根据所述动态谐波相量参数化模型实现电网信号的近似表示模型，从而可根据这一表示结果，实现衰减指数函数特征相量的估计；3)根据步骤2)构建的电网信号模型，基于最小二乘法实现衰减指数函数特征相量的估计。2.如权利要求1在所述的一种基于衰减指数函数模型的动态谐波相量测量方法，其特征在于步骤1)中，动态谐波相量参数化模型的构建为：选定模型参数值τ
k,h
，根据式(2)构建谐波相量的参数化模型，其中，τ
k,h
为模型中提前给定的第k个指数衰减函数的时间常数；p
k,h
表征第k个指数衰减函数的对应的相量值，h为谐波...

【专利技术属性】
技术研发人员：李景一，杨建军，谢小荣，陈垒，施朝晖，王克，傅春翔，郦洪柯，杨林刚，杨文斌，
申请(专利权)人：清华大学，
类型：发明
国别省市：