CN109444537A - 一种计及带外干扰的自适应同步相量测量方法 - Google Patents
一种计及带外干扰的自适应同步相量测量方法 Download PDFInfo
- Publication number
- CN109444537A CN109444537A CN201811215147.7A CN201811215147A CN109444537A CN 109444537 A CN109444537 A CN 109444537A CN 201811215147 A CN201811215147 A CN 201811215147A CN 109444537 A CN109444537 A CN 109444537A
- Authority
- CN
- China
- Prior art keywords
- band interference
- frequency
- phasor
- value
- fundamental wave
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
- G01R23/165—Spectrum analysis; Fourier analysis using filters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Noise Elimination (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明公开了一种计及带外干扰的自适应同步相量测量方法,涉及电力***同步相量测量领域,用以解决现有技术下用于测量电力***同步相量的算法并未在对基波相量进行泰勒建模时同时考虑带外干扰的泰勒模型,从而导致在带外干扰下相量测量误差大,TVE值难以满足标准。本发明新增了带外信号的泰勒相量模型以及对带外信号频率、幅值的提取和判断过程,通过判断是否存在带外干扰建立不同的相量模型,调用不同基波相量修正矩阵,在增加有限运算量的基础上补充了现有基于泰勒模型的动态同步相量测量算法的不足。能够有效判断带外干扰是否存在于电力信号中,当带外干扰和频率偏移同时存在时,提高了电力信号基波同步相量的测量精度,满足测量标准要求。
Description
技术领域
本发明涉及一种计及带外干扰的自适应同步相量测量方法,属于电力***同步相量测 量领域。
背景技术
随着电力***电能传输的发展,高压直流输电技术逐渐成为研究热点。但负载端负荷 动态变化导致的基波频率偏移和换流操作引入的带外干扰信号也给同步相量测量带来一定 负面影响。这两种现象的存在使得同步相量测量算法难以对基波信号进行提取,测量精度 急剧下降,进而限制了广域测量***的高级应用。
现有的同步相量测量算法大致分为两类,一类是以短时傅里叶变换和插值法为代表的 静态算法,另一类是泰勒模型法和其扩展算法的动态算法。这些算法在不同工况下都能应 对频率偏移带来的影响,但信号同时出现频率偏移和带外干扰时,短时傅里叶变换法、插 值法和泰勒模型法需要依赖较高的采样频率或较多的采样点数,增加了运算负担。
虽然通过级联滤波器及其可调整的基波滤波频率实现了对带外干扰以及频率偏移引起 的频谱泄露的抑制,但其滤波器参数设计较复杂且难以保证算法的响应时间满足测量标准。 而基于泰勒模型的动态同步相量测量算法,其改进的泰勒模型能够处理频率偏移较大的各 个动态工况,在运算量、响应时间、抗噪能力、测量精度等方面具有较好的性能体现,但 是其在相量模型建立时并未考虑加入带外干扰的相量模型,使得单纯利用单一频点对信号 进行滤波得到的相量值含有较多的带外泄露值,从而导致基于泰勒模型的动态同步相量测 量算法在带外干扰的工况下的同步相量测量误差较大,难以满足测量标准对同步相量测量 误差的要求。
因此有必要在基波相量模型的基础上增加带外干扰相量模型,改进基于泰勒模型的动 态同步相量测量算法,设计一种计及带外干扰的自适应同步相量测量方法。
发明内容
本发明的目的是:为解决当带外干扰和频率偏移同时存在时,现有的同步相量测算方法 测量精度急剧下降,难以满足同步相量测量要求,本发明提供一种计及带外干扰的自适应同 步相量测量方法。
本发明为实现其发明目的,所采用的技术方案如下:
S1:对电力信号采样得到电压/电流离散序列x(n),利用滤波频率ω0=2πf0/fs和ωb=2πfb/fs在一个数据窗内获得基波初估相量值和带外干扰初估相量值其中,其中n表示采样时刻,f0表示基波频率,fb表示定义带外干扰范围的中心 位置,fs表示对电力信号的采样频率,ω0表示额定基波角频率,ωb表示带外范围中心点角 频率,l表示数据窗序号;
S2:利用S1中得到的基波初估相量值求得基波粗估频率根据基波和带外 干扰间的泄露关系,再利用S1中得到的初估相量值以及基波粗估频率扣除基波对带外信号的泄露值,获得较为纯净的带外相量值Xb,以此获得带外粗估频率
S3:判断带外粗估频率是否在带外干扰定义区间内,
若在带外干扰定义区间内,再继续判断扣除泄露值后带外干扰相量幅值的大小,此时 的幅值大小分为三种情况:(1)大于δ2;(2)小于δ1;(3)在[δ1,δ2]区间内;δ1和δ2为根据标准定义并结合实际情况设置的幅值上限和幅值下限值;
若不在带外干扰定义区间内,无需进一步判断其相量幅值大小;
S4:如果幅值大于δ2,则它被视作间谐波,用相关滤波算法来处理;
如果幅值小于δ1,则直接用MDSEA算法处理;
如果幅值在[δ1,δ2]区间内,则根据S2中求得的基波粗估频率带外粗估频率建 立基波相量泰勒模型和带外相量泰勒模型,再用S1中得到的初估相量值求得其具体参数,得到基波相量测量值X0(l);
如果不含带外干扰信号,则直接用MDSEA算法处理。
具体地,所述S1具体为:用矩形截断电压电压/电流离散序列x(n),并选用滤波频率 ω0=2πf0/fs和ωb=2πfb/fs分别进行离散傅里叶变换得到基波和带外初估相量值 其中,n=tfs,t表示时间,fs表示采样频率;f0为50Hz,fb取带外干扰定义范围的中点17.5Hz、87.5Hz二者之一;离散傅里叶变换操作如下:
其中,x(n)表示离散采样序列,h(n)表示矩形窗序列,n表示采样时刻,N表示一个周波内的采样点数,l表示第n个采样点距数据窗中心相差的的采样点数。
具体地,所述S2具体为,先利用得到相位差,再根据频率定义公式、基波频 率和相位差与时间差的比值可得基波粗估频率相位差由S1中得到的历史相量值进行求 解;
S2更具体的操作如下:
设和为不同数据窗下傅里叶变换结果的历史求和, 则相位差根据频率定义,则长度为2M+1的数据窗的中点处的基波粗 估频率其中,angle()表示取角度函数,()*表示取共轭,∑()表 示求和函数,round()表示取整函数。
进一步地,所述S2中,基波和带外干扰之间的泄露关系满足以下公式:
其中, ωb=2πfb/fs,是基波在带外频率点处的泄露值,需要扣除这部分泄露值,由此获得较为纯净的带外相量值Xb;
S22:由相角的定义根据带外相量值Xb得到带外干扰的相角信息利用带 外相角高阶信息将表示成多项式形式:其中,表示带 外相角,表示带外相角变化率,由lF个历史带外干扰的相角信息求得;表示带 外相角变化的高阶信息,i表示当前采样点到之间的采样点数,i=1,2…lF,ΔT=1/fs, ΔT表示每个之间的时间间隔;
S23:再次根据频率定义求得带外干扰粗估频率
具体地,所述S3中,所述带外干扰定义频率区间是 [10Hz,f0-fs/2]∪[75Hz,2thharmonic],其中fs是对电力信号的采样频率即相量报告频率, 取50Hz;2th harmonic是100Hz,则带外干扰频率区间[10Hz,25Hz]∪[75Hz,100Hz];判断 S2中求得的带外粗估频率是否属于该范围,如果属于该范围,由S2中得到的带外相量 值Xb求得带外幅值|Xb|=abs(Xb),判断其是否大于0.01,如果是则说明考虑的电力信号 中包含带外干扰,反之则不含带外干扰。
具体地,所述S4中,所述MDSEA算法处理的过程为:将S2中求得的基波粗估频率建立基波相量泰勒模型,再用S1中得到的初估相量值求得其具体参数,得到基波相量测量值X0(l)。
优选地,所述S4中,所述δ1的取值为0.01,δ2的取值为0.1。
具体地,所述S4具体包括:
如果所述S3判断出所分析的信号中含有带外干扰成分,则:
S41:
使电力信号相量模型X(t)为:
式中,泰勒模型分别表示为和t表示时间,a(k)和b(k)分 别为a(t)和b(t)的k阶导数,a(t)表示基波信号的动态变化,b(t)表示带外信号的动态变化,
S42:对电力信号相量模型求取实部Re[X(t)]再进行离散化操作得到:
式中,α(k)(n)和β(k)(n)分别是离散后的a(k)(t)和b(k)(t)值;
S43:利用所述S2得到的基波粗估频率和带外粗估频率分别对S42中的公式进行 滤波得到:
式中,基波泰勒级数矩阵为A=[α(0),…,α(K)],带外干扰泰勒级数矩阵为 B=[β(0),…,β(K)],矩阵C,D,E,F的表达式分别如下:
矩阵中的H表达式为J表达式为
S44:再对S3中的滤波公式进行虚实部分离得到下式:
式中,[]R表示矩阵中每个元素的实部,[]I表示矩阵中每个元素的虚部,矩阵令G=[G1 G2], [AB]T=[AR AI BR BI]T,则获得一个数据窗内新的傅里叶变换值与待求解参 数之间的矩阵关系G=[G1 G2];
S45:利用最小二乘求得基波同步相量测量值,具体操作如下:
基于此式,又考虑到参考时刻与GPS报告时刻之间的时间差Δτ,进行如下相移操作 得到报告时刻的基波相量值:
其中,trep表示GPS的报告时间;
如果所述S3判断出所分析的信号中不含带外干扰成分,则:
电力信号相量模型X(t)中只含有基波信息a(t)表示基波 信号的动态变化,用S2中得到的基波粗估频率对数据窗内的离散序列进行滤 波,此时S43中的简化为S44 中虚实部分离后的G=[G1 G2]简化为G=[G1],再最小二乘法求得 从而得到基波泰勒级数矩阵A=[α(0),…,α(K)]中的参数, 最后同样考虑参考时刻与GPS报告时刻之间的时间差Δτ对上述参数进行相移操 作得到报告时刻基波同步相量值
综上所述,与现有技术相比,本发明的有益效果是:
1.本发明采用的同步相量测量算法能够有效判断信号中是否含有带外干扰,根据前期 判断结果选择建立不同的信号相量模型,再调用不同的修正矩阵G完成对基波相量参数的 求解,对具体信号有一定的适应性。
2.本发明在频率偏移和带外干扰同时存在的工况下,在原有基波相量模型基础上新建 带外相量模型,模型上更加精确。面对理想信号的带外干扰工况,能够满足IEEE StdC37.118.1TM-2011的测量标准对1.3%TVE要求,相比MDSEA算法,基波幅值和相角测量误差更小。
3.本发明噪声存在的实际工况下,算法具有一定的稳定性和抗噪性。相比MDSEA算法, 基波相量的TVE值、基波幅值和相角测量误差更小。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要 使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例, 对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得 其他的附图。通过附图所示,本发明的上述及其它目的、特征和优势将更加清晰。在全部 附图中相同的附图标记指示相同的部分。并未刻意按实际尺寸等比例缩放绘制附图,重点 在于示出本发明的主旨。
图1为本发明方法流程图;
图2为存在带外干扰时的相量模型图;
图3为仿真实验中fb范围是10Hz-25Hz时,MDSEA算法和本发明算法得到的TVE最大值 对比图;
图4为仿真实验中fb范围是75Hz-100Hz时,MDSEA算法和本发明算法得到的TVE最大 值;
图5为仿真实验中fb范围是10Hz-25Hz、信号含有40dB白噪声时,MDSEA算法和本发明算法得到的TVE最大值对比图;
图6为仿真实验中fb范围是75Hz-100Hz、信号含有40dB白噪声时,MDSEA算法和本发 明算法得到的TVE最大值。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的 附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本 发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员 在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
下面结合附图和具体实施方式对本发明做进一步的详细说明。
实施例
一种计及带外干扰的自适应同步相量测量方法,如图1所示,包括以下步骤:
S1:对电力信号采样得到电压/电流离散序列x(n),利用滤波频率ω0=2πf0/fs和ωb=2πfb/fs在一个数据窗内获得基波初估相量值和带外干扰初估相量值其中,其中n表示采样时刻,f0表示基波频率,fb表示定义带外干扰范围的中心 位置,fs表示对电力信号的采样频率,ω0表示额定基波角频率,ωb表示带外范围中心点角 频率,l表示数据窗序号。
S2:利用S1中得到的基波初估相量值求得基波粗估频率根据基波和带外 干扰间的泄露关系,再利用S1中得到的初估相量值以及基波粗估频率扣除基波对带外信号的泄露值,获得较为纯净的带外相量值Xb,以此获得带外粗估频率
S3:判断带外粗估频率是否在带外干扰定义区间内,则说明电力信号中包含带外信 号,再继续判断扣除泄露值后带外干扰相量幅值的大小,此时的幅值大小分为三种情况:(1) 大于δ2;(2)小于δ1;(3)在[δ1,δ2]区间内;δ1和δ2为根据标准定义并结合实际情况设置 的幅值上限和幅值下限值,所述δ1的取值为0.01,δ2的取值为0.1。
若不在带外干扰定义区间内,无需进一步判断其相量幅值大小;
S4:如果幅值大于δ2,则它被视作间谐波,用相关滤波算法来处理;
如果幅值小于δ1,则认为它对基波相量测量的影响可忽略,直接用MDSEA算法处理: 将S2中求得的基波粗估频率建立基波相量泰勒模型,再用S1中得到的初估相量值求得其具体参数,得到基波相量测量值X0(l);
如果幅值在[δ1,δ2]区间内,则表明所含带外干扰对基波相量测量有不良影响,则根据 S2中求得的基波粗估频率带外粗估频率建立基波相量泰勒模型和带外相量泰勒模 型,再用S1中得到的初估相量值求得其具体参数,得到基波相量测量 值X0(l);
如果不含带外干扰信号,则将S2中求得的基波粗估频率建立基波相量泰勒模型,再 用S1中得到的初估相量值求得其具体参数,得到基波相量测量值X0(l)。
作为一种优选的实施方式,所述S1具体为:用矩形截断电压电压/电流离散序列x(n), 并选用滤波频率ω0=2πf0/fs和ωb=2πfb/fs分别进行离散傅里叶变换得到基波和带外初估 相量值其中,n=tfs,t表示时间,fs表示采样频率;f0为50Hz,fb取带外干扰定义范围的中点17.5Hz、87.5Hz二者之一;离散傅里叶变换操作如下:
其中,x(n)表示离散采样序列,h(n)表示矩形窗序列,n表示采样时刻,N表示一个周波内的采样点数,l表示第n个采样点距数据窗中心相差的的采样点数。
作为一种优选的实施方式,所述S2具体为,先利用得到相位差,再根据频率 定义公式、基波频率和相位差与时间差的比值可得基波粗估频率相位差由S1中得到的 历史相量值进行求解;
S2更具体的操作如下:
设和为不同数据窗下傅里叶变换结果的历史求和, 则相位差根据频率定义,则长度为2M+1的数据窗的中点处的基波粗 估频率其中,angle()表示取角度函数,()*表示取共轭,∑()表 示求和函数,round()表示取整函数。
作为一种优选的实施方式,所述S2中,基波和带外干扰之间的泄露关系满足以下公式:
其中, ωb=2πfb/fs,是基波在带外频率点处的泄露值,需要扣除这部分泄露值,由此获得较为纯净的带外相量值Xb;
S22:由相角的定义根据带外相量值Xb得到带外干扰的相角信息利用带 外相角高阶信息将表示成多项式形式:其中,表示带 外相角,表示带外相角变化率,由lF个历史带外干扰的相角信息求得;表示带 外相角变化的高阶信息,i表示当前采样点到之间的采样点数,i=1,2…lF,ΔT=1/fs, ΔT表示每个之间的时间间隔;
S23:再次根据频率定义求得带外干扰粗估频率
作为一种优选的实施方式,所述S3中,所述带外干扰定义频率区间是[10Hz,f0-fs/2]∪[75Hz,2th harmonic],其中fs是对电力信号的采样频率即相量报告频率, 取50Hz;2th harmonic是100Hz,则带外干扰频率区间[10Hz,25Hz]∪[75Hz,100Hz];判断 S2中求得的带外粗估频率是否属于该范围,如果属于该范围,由S2中得到的带外相量 值Xb求得带外幅值|Xb|=abs(Xb),判断其是否大于0.01,如果是则说明考虑的电力信号 中包含带外干扰,反之则不含带外干扰。
作为一种优选的实施方式,如果S3判断出所分析的信号中含有带外干扰成分,则建立 如图2所示的相量模型,如果所述S3判断出所分析的信号中含有带外干扰成分,则:
S41:
使电力信号相量模型X(t)为:
式中,泰勒模型分别表示为和t表示时间,a(k)和b(k)分 别为a(t)和b(t)的k阶导数,a(t)表示基波信号的动态变化,b(t)表示带外信号的动态变化,
S42:对电力信号相量模型求取实部Re[X(t)]再进行离散化操作得到:
式中,α(k)(n)和β(k)(n)分别是离散后的a(k)(t)和b(k)(t)值;
S43:利用所述S2得到的基波粗估频率和带外粗估频率分别对S42中的公式进行 滤波得到:
式中,基波泰勒级数矩阵为A=[α(0),…,α(K)],带外干扰泰勒级数矩阵为 B=[β(0),…,β(K)],矩阵C,D,E,F的表达式分别如下:
矩阵中的H表达式为J表达式为
S44:再对S3中的滤波公式进行虚实部分离得到下式:
式中,[]R表示矩阵中每个元素的实部,[]I表示矩阵中每个元素的虚部,矩阵令G=[G1 G2], [AB]T=[AR AI BR BI]T,则获得一个数据窗内新的傅里叶变换值与待求解参 数之间的矩阵关系G=[G1 G2];
S45:利用最小二乘求得基波同步相量测量值,具体操作如下:
基于此式,又考虑到参考时刻与GPS报告时刻之间的时间差Δτ,进行如下相移操作 得到报告时刻的基波相量值:
其中,trep表示GPS的报告时间;
如果所述S3判断出所分析的信号中不含带外干扰成分,则:
电力信号相量模型X(t)中只含有基波信息a(t)表示基波 信号的动态变化,用S2中得到的基波粗估频率对数据窗内的离散序列进行滤 波,此时S43中的简化为S44 中虚实部分离后的G=[G1 G2]简化为G=[G1],再最小二乘法求得 从而得到基波泰勒级数矩阵A=[α(0),…,α(K)]中的参数, 最后同样考虑参考时刻与GPS报告时刻之间的时间差Δτ对上述参数进行相移操 作得到报告时刻基波同步相量值。
本发明新增了带外信号的泰勒相量模型以及对带外信号频率、幅值的提取和判断过程, 通过判断是否存在带外干扰建立不同的相量模型,调用不同基波相量修正矩阵,在增加有 限运算量的基础上补充了现有基于泰勒模型的动态同步相量测量算法的不足。能够有效判 断带外干扰是否存在于电力信号中,当带外干扰和频率偏移同时存在时,提高了电力信号 基波同步相量的测量精度,满足测量标准要求。
仿真实验:
为验证本发明方法在带外干扰下基波同步相量测量精度,分别用本发明算法 和已发表的考虑频率偏移的泰勒模型同步相量测量算法(A Modified DynamicSynchrophasor Estimation Algorithm Considering Frequency Deviation, MDSEA),根据IEEE Std C37.118.1TM-2011测量标准的内容,对频率偏移和带外 干扰同时存在以及噪声影响情况下的两种信号进行仿真测试。在实验中,基波额 定频率f0=50Hz,采样频率fs=N×f0=48×50=2400Hz,带外干扰频率范围 为[10Hz,25Hz]∪[75Hz,100Hz],S1中的带外滤波频率设置为fb=87.5Hz,基波
1.频率偏移与带外干扰工况下的仿真
根据IEEE Std C37.118.1TM-2011测量标准,设置被测信号表达式为:
x(t)=cos[2πt(f0+Δf)]+0.1cos(2πfbt),
上式在基波信号的基础上叠加了一个10%的带外干扰成分,fb测试范围为(10Hz-25Hz) 和(75Hz-100Hz),其步进为1Hz,每个频率测试时间为1s;基波频率偏移为Δf=-2.5-2.5Hz。 对其基波TVE最大值进行统计,绘制出图3和图4。由图3和图4可知,MDSEA受带外信号 的影响更为严重,部分测量结果超出IEEE标准对TVE为1.3%的测量要求。相比之下,本发 明因为能准确判断带外信号的存在,并在模型中考虑了对带外干扰信号的处理,分别提取 基波和带外信息,所以具有更高的同步相量测量精度。
2.频率偏移与带外干扰以及噪声工况下的仿真
在实际情况下,噪声和带外可能是同时存在的。为了验证算法的稳定性和抗噪性,设 置信号模型为:
x(t)=cos[2π(f0+Δf)+t]+0.1cos(2πfbt)+N(t)
上式中参数fb、Δf与仿真1设置相同,N(t)为高斯白噪声,信噪比设置为40dB;同样 对其基波TVE最大值进行统计,绘制出图5和图6。从图中可以看出,相同参数设置下,本发明算法在噪声工况下对也能对带外干扰信号进行提取,本发明算法得到的基波TVE最大值明显小于MDSEA算法,同步相量测量精度更高。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何属 于本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵 盖在本发明的保护范围之内。
Claims (8)
1.一种计及带外干扰的自适应同步相量测量方法,其特征在于,包括以下步骤:
S1:对电力信号采样得到电压/电流离散序列x(n),利用滤波频率ω0=2πf0/fs和ωb=2πfb/fs在一个数据窗内获得基波初估相量值和带外干扰初估相量值其中,其中n表示采样时刻,f0表示基波频率,fb表示定义带外干扰范围的中心位置,fs表示对电力信号的采样频率,ω0表示额定基波角频率,ωb表示带外范围中心点角频率,l表示数据窗序号;
S2:利用S1中得到的基波初估相量值求得基波粗估频率根据基波和带外干扰间的泄露关系,再利用S1中得到的初估相量值以及基波粗估频率扣除基波对带外信号的泄露值,获得较为纯净的带外相量值Xb,以此获得带外粗估频率
S3:判断带外粗估频率是否在带外干扰定义区间内,
若在带外干扰定义区间内,再继续判断扣除泄露值后带外干扰相量幅值的大小,此时的幅值大小分为三种情况:(1)大于δ2;(2)小于δ1;(3)在[δ1,δ2]区间内;δ1和δ2为根据标准定义并结合实际情况设置的幅值上限和幅值下限值;
若不在带外干扰定义区间内,无需进一步判断其相量幅值大小;
S4:如果幅值大于δ2,则它被视作间谐波,用相关滤波算法来处理;
如果幅值小于δ1,则直接用MDSEA算法处理;
如果幅值在[δ1,δ2]区间内,则根据S2中求得的基波粗估频率带外粗估频率建立基波相量泰勒模型和带外相量泰勒模型,再用S1中得到的初估相量值求得其具体参数,得到基波相量测量值X0(l);
如果不含带外干扰信号,则直接用MDSEA算法处理。
2.根据权利要求1所述的一种计及带外干扰的自适应同步相量测量方法,其特征在于,所述S1具体为:用矩形截断电压电压/电流离散序列x(n),并选用滤波频率ω0=2πf0/fs和ωb=2πfb/fs分别进行离散傅里叶变换得到基波和带外初估相量值其中,n=tfs,t表示时间,fs表示采样频率;f0为50Hz,fb取带外干扰定义范围的中点17.5Hz、87.5Hz二者之一;离散傅里叶变换操作如下:
其中,x(n)表示离散采样序列,h(n)表示矩形窗序列,n表示采样时刻,N表示一个周波内的采样点数,l表示第n个采样点距数据窗中心相差的的采样点数。
3.根据权利要求1所述的一种计及带外干扰的自适应同步相量测量方法,其特征在于,所述S2具体为,先利用得到相位差,再根据频率定义公式、基波频率和相位差与时间差的比值可得基波粗估频率相位差由S1中得到的历史相量值进行求解;
S2更具体的操作如下:
设和为不同数据窗下傅里叶变换结果的历史求和,则相位差根据频率定义,则长度为2M+1的数据窗的中点处的基波粗估频率其中,angle()表示取角度函数,()*表示取共轭,∑()表示求和函数,round()表示取整函数。
4.根据权利要求1所述的一种计及带外干扰的自适应同步相量测量方法,其特征在于,所述S2中,基波和带外干扰之间的泄露关系满足以下公式:
其中, ωb=2πfb/fs,是基波在带外频率点处的泄露值,需要扣除这部分泄露值,由此获得较为纯净的带外相量值Xb;
S22:由相角的定义根据带外相量值Xb得到带外干扰的相角信息利用带外相角高阶信息将表示成多项式形式:其中,表示带外相角,表示带外相角变化率,由lF个历史带外干扰的相角信息求得;表示带外相角变化的高阶信息,i表示当前采样点到之间的采样点数,i=1,2…lF,ΔT=1/fs,ΔT表示每个之间的时间间隔;
S23:再次根据频率定义求得带外干扰粗估频率
5.根据权利要求1所述的一种计及带外干扰的自适应同步相量测量方法,其特征在于,所述S3中,所述带外干扰定义频率区间是[10Hz,f0-fs/2]∪[75Hz,2th harmonic],其中fs是对电力信号的采样频率即相量报告频率,取50Hz;2th harmonic是100Hz,则带外干扰频率区间[10Hz,25Hz]∪[75Hz,100Hz];判断S2中求得的带外粗估频率是否属于该范围,如果属于该范围,由S2中得到的带外相量值Xb求得带外幅值|Xb|=abs(Xb),判断其是否大于0.01,如果是则说明考虑的电力信号中包含带外干扰,反之则不含带外干扰。
6.根据权利要求1所述的一种计及带外干扰的自适应同步相量测量方法,其特征在于,所述S4中,所述MDSEA算法处理的过程为:将S2中求得的基波粗估频率建立基波相量泰勒模型,再用S1中得到的初估相量值求得其具体参数,得到基波相量测量值X0(l)。
7.根据权利要求1所述的一种计及带外干扰的自适应同步相量测量方法,其特征在于,所述S4中,所述δ1的取值为0.01,δ2的取值为0.1。
8.根据权利要求1所述的一种计及带外干扰的自适应同步相量测量方法,其特征在于,所述S4具体包括:
如果所述S3判断出所分析的信号中含有带外干扰成分,则:
S41:
使电力信号相量模型X(t)为:
式中,泰勒模型分别表示为和t表示时间,a(k)和b(k)分别为a(t)和b(t)的k阶导数,a(t)表示基波信号的动态变化,b(t)表示带外信号的动态变化,
S42:对电力信号相量模型求取实部Re[X(t)]再进行离散化操作得到:
式中,α(k)(n)和β(k)(n)分别是离散后的a(k)(t)和b(k)(t)值;
S43:利用所述S2得到的基波粗估频率和带外粗估频率分别对S42中的公式进行滤波得到:
式中,基波泰勒级数矩阵为A=[α(0),…,α(K)],带外干扰泰勒级数矩阵为B=[β(0),…,β(K)],矩阵C,D,E,F的表达式分别如下:
矩阵中的H表达式为J表达式为
S44:再对S3中的滤波公式进行虚实部分离得到下式:
式中,[]R表示矩阵中每个元素的实部,[]I表示矩阵中每个元素的虚部,矩阵令G=[G1 G2],[A B]T=[AR AI BR BI]T,则获得一个数据窗内新的傅里叶变换值与待求解参数之间的矩阵关系G=[G1 G2];
S45:利用最小二乘求得基波同步相量测量值,具体操作如下:
基于此式,又考虑到参考时刻与GPS报告时刻之间的时间差Δτ,进行如下相移操作得到报告时刻的基波相量值:
其中,trep表示GPS的报告时间;
如果所述S3判断出所分析的信号中不含带外干扰成分,则:
电力信号相量模型X(t)中只含有基波信息a(t)表示基波信号的动态变化,用S2中得到的基波粗估频率对数据窗内的离散序列进行滤波,此时S43中的简化为S44中虚实部分离后的G=[G1 G2]简化为G=[G1],再最小二乘法求得从而得到基波泰勒级数矩阵A=[α(0),…,α(K)]中的参数,最后同样考虑参考时刻与GPS报告时刻之间的时间差Δτ对上述参数进行相移操作得到报告时刻基波同步相量值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811215147.7A CN109444537B (zh) | 2018-10-18 | 2018-10-18 | 一种计及带外干扰的自适应同步相量测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811215147.7A CN109444537B (zh) | 2018-10-18 | 2018-10-18 | 一种计及带外干扰的自适应同步相量测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109444537A true CN109444537A (zh) | 2019-03-08 |
CN109444537B CN109444537B (zh) | 2020-06-19 |
Family
ID=65546795
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811215147.7A Active CN109444537B (zh) | 2018-10-18 | 2018-10-18 | 一种计及带外干扰的自适应同步相量测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109444537B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110376435A (zh) * | 2019-08-05 | 2019-10-25 | 西南交通大学 | 一种基于带外-多频率模型的同步相量测量方法 |
CN112485522A (zh) * | 2020-12-09 | 2021-03-12 | 国网四川省电力公司电力科学研究院 | 基于电能数据感知的平顶窗函数同步相量测量方法及装置 |
CN112485524A (zh) * | 2020-11-10 | 2021-03-12 | 广东电网有限责任公司广州供电局 | 一种用于pmu测试的高精度校准器相量计算方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1477401A (zh) * | 2003-07-18 | 2004-02-25 | 清华大学 | 高精度的同步相量测量方法 |
CN104020352A (zh) * | 2014-06-09 | 2014-09-03 | 华北电力大学 | 一种适用于m类pmu单元的同步相量测量方法 |
CN104502703A (zh) * | 2014-12-13 | 2015-04-08 | 西南交通大学 | 基于频域动态模型的电力信号同步相量测量方法 |
CN105759118A (zh) * | 2016-03-02 | 2016-07-13 | 中国电力科学研究院 | 一种同步相量的量测方法 |
CN107345984A (zh) * | 2017-06-23 | 2017-11-14 | 华北电力大学 | 一种基于信号识别的自适应同步相量测量方法 |
CN107589299A (zh) * | 2017-08-03 | 2018-01-16 | 西南交通大学 | 基于多频率相量模型的电力信号同步相量测量方法 |
CN107656151A (zh) * | 2017-08-10 | 2018-02-02 | 国家电网公司 | 实现全时序电网量测数据同步采样及处理的测控方法 |
CN108614155A (zh) * | 2018-05-31 | 2018-10-02 | 许继集团有限公司 | 一种加入汉明窗的同步相量测量方法及*** |
-
2018
- 2018-10-18 CN CN201811215147.7A patent/CN109444537B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1477401A (zh) * | 2003-07-18 | 2004-02-25 | 清华大学 | 高精度的同步相量测量方法 |
CN104020352A (zh) * | 2014-06-09 | 2014-09-03 | 华北电力大学 | 一种适用于m类pmu单元的同步相量测量方法 |
CN104502703A (zh) * | 2014-12-13 | 2015-04-08 | 西南交通大学 | 基于频域动态模型的电力信号同步相量测量方法 |
CN105759118A (zh) * | 2016-03-02 | 2016-07-13 | 中国电力科学研究院 | 一种同步相量的量测方法 |
CN107345984A (zh) * | 2017-06-23 | 2017-11-14 | 华北电力大学 | 一种基于信号识别的自适应同步相量测量方法 |
CN107589299A (zh) * | 2017-08-03 | 2018-01-16 | 西南交通大学 | 基于多频率相量模型的电力信号同步相量测量方法 |
CN107656151A (zh) * | 2017-08-10 | 2018-02-02 | 国家电网公司 | 实现全时序电网量测数据同步采样及处理的测控方法 |
CN108614155A (zh) * | 2018-05-31 | 2018-10-02 | 许继集团有限公司 | 一种加入汉明窗的同步相量测量方法及*** |
Non-Patent Citations (10)
Title |
---|
LING FU: "A Modified Dynamic Synchrophasor Estimation Algorithm Considering Frequency Deviation", 《IEEE TRANSACTIONS ON SMART GRID 》 * |
常乃超: "动态条件下同步相量测量装置的数字滤波器及计算优化", 《电力***自动化》 * |
张佳怡: "计及频率偏移的动态谐波相量测量算法", 《中国电机工程学报》 * |
王科: "基于扩展卡尔曼滤波频率跟踪的DFT同步相量测量算法", 《电网技术》 * |
白莎: "基于多频率相量模型的动态同步相量测量算法", 《中国电机工程学报》 * |
符玲: "基于时频信息的动态同步相量测量算法", 《中国电机工程学报》 * |
符玲: "基于频域动态模型的同步相量测量算法", 《中国电机工程学报》 * |
麦瑞坤: "动态条件下的同步相量测量算法的研究", 《中国电机工程学报》 * |
麦瑞坤: "基于泰勒展开模型的同步相量估计新算法", 《电力***自动化》 * |
麦瑞坤: "滤除衰减直流分量的动态同步相量测量算法", 《中国电机工程学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110376435A (zh) * | 2019-08-05 | 2019-10-25 | 西南交通大学 | 一种基于带外-多频率模型的同步相量测量方法 |
CN112485524A (zh) * | 2020-11-10 | 2021-03-12 | 广东电网有限责任公司广州供电局 | 一种用于pmu测试的高精度校准器相量计算方法 |
CN112485524B (zh) * | 2020-11-10 | 2024-01-19 | 广东电网有限责任公司广州供电局 | 一种用于pmu测试的高精度校准器相量计算方法 |
CN112485522A (zh) * | 2020-12-09 | 2021-03-12 | 国网四川省电力公司电力科学研究院 | 基于电能数据感知的平顶窗函数同步相量测量方法及装置 |
CN112485522B (zh) * | 2020-12-09 | 2023-05-16 | 国网四川省电力公司电力科学研究院 | 基于电能数据感知的平顶窗函数同步相量测量方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN109444537B (zh) | 2020-06-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105223418B (zh) | 次同步和超同步谐波相量的测量方法及测量装置 | |
CN102539915B (zh) | 时延傅立叶变换测频法精确计算电力谐波参数方法 | |
CN103353558B (zh) | 一种电能质量监测方法 | |
CN109444537A (zh) | 一种计及带外干扰的自适应同步相量测量方法 | |
CN101363896B (zh) | 一种高精度电压暂态事件检测与录波方法 | |
CN102338827B (zh) | 用于电网波形畸变分析和电力谐波参数自动监测方法 | |
EP2690451A1 (en) | System frequency measurement method, synchrophasor measurement method and device thereof | |
CN102253282B (zh) | 纳托尔窗函数连续频谱内插电力谐波参数获取方法 | |
CN102331526B (zh) | 汉宁窗函数连续频谱内插电力谐波参数获取方法 | |
CN106405230B (zh) | 频率测量方法和装置 | |
CN108614155B (zh) | 一种加入汉明窗的同步相量测量方法及*** | |
CN104316768B (zh) | 一种三相不平衡扰动源定位的负序阻抗参数估算方法 | |
JP6416072B2 (ja) | 同期フェーザ測定装置およびパルス生成装置 | |
CN107247182A (zh) | 一种基于量测相量数据的间谐波分量还原方法 | |
Romano et al. | An enhanced interpolated-modulated sliding DFT for high reporting rate PMUs | |
CN106841778A (zh) | 基于pmu实现的次同步和超同步谐波参数的处理方法 | |
CN111122952B (zh) | 一种快速检测三相电压暂降的方法 | |
CN109298362B (zh) | 一种配电线路故障指示器三相同步精度测试方法、装置和*** | |
CN108761222A (zh) | 一种不平衡工况下的电网电压序分量快速提取***及方法 | |
CN109669072A (zh) | 一种配电网的自适应同步相量量测方法 | |
CN106154037A (zh) | 一种基于校验的同步相量自适应计算方法 | |
CN112748284A (zh) | 电力***同步波形与广谱相量测量方法及装置 | |
CN103983847A (zh) | 一种同步相量测量中基于rls的自适应频率跟踪测量方法 | |
CN103605904B (zh) | 基于误差估算的自补偿电力***幅值算法 | |
CN111625769A (zh) | 一种基于拉格朗日插值和三次指数平滑的pmu-scada数据对时与融合方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |