首页 / 专利库 / 信号处理 / 采样频率 / 利用频域信息滤除衰减直流的电力信号同步相量测量方法

利用频域信息滤除衰减直流的电信号同步相量测量方法

阅读:1发布:2022-02-11

专利汇可以提供利用频域信息滤除衰减直流的电信号同步相量测量方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了利用频域信息滤除衰减直流的电 力 信号 同步相量测量方法,其步骤是,A、 数据采集 与预估计处理:采集电力信号x(t)得到 电网 的电力信号离散序列x(n),利用 加窗 傅里叶变换对电力信号离散序列x(n)滤波,得到电力信号离散序列x(n)在5个 角 速度 ω处的相量测量预估计值X(ω);B、修正相量预估计值:对得到的5个相量测量预估计值X(ω),用考虑衰减 直流分量 的动态信号模型,求得电力信号基波分量的时变特性的各阶导数值,并进行 相移 运算得到当前时刻trep电网的电力信号的基波相量值该方法能够消除或减弱信号基波分量的时变性以及衰减直流分量对相量测量的影响,并大大提高信号的相量测量 精度 。,下面是利用频域信息滤除衰减直流的电信号同步相量测量方法专利的具体信息内容。

1.利用频域信息滤除衰减直流的电信号同步相量测量方法,其步骤为:
A、数据采集与预估计处理
电网中的电流互感器或电压互感器测量电网的电流信号或电压信号得到电网的电力信号x(t),采集电力信号x(t)得到电网的电力信号离散序列x(n),其中t为时间,n为采样时刻点;数字信号处理器利用加窗傅里叶变换对电力信号离散序列x(n)滤波,得到电力信号离散序列x(n)在5个速度ω处的相量测量预估计值X(ω);其中ω=2πf,f为角速度ω对应的频率,f=0,1,49,50,51Hz;
B、修正相量预估计值
数字信号处理器对得到的5个相量测量预估计值X(ω),用考虑衰减直流分量的动态信号模型,求得电力信号基波分量的时变特性的各阶导数值,并进行相移运算得到当前时刻trep电网的电力信号的基波相量值
2.如权利1所述的利用频域信息滤除衰减直流的电力信号同步相量测量方法,其特征在于:所述的步骤B中考虑衰减直流分量的动态信号模型用以下方法建立:
(B1)数字信号处理器利用低频带限相量a(t0)的泰勒级数近似表示电力信号基波分量的时变特性xg(t0), 其中,Δt为当
前时刻trep与基准时刻t0的差,基准时刻t0为加窗傅里叶变换的数据窗中心时刻,k为泰(k)
勒级数的阶次,k=0,1,2;K为泰勒级数的最高阶次,取为2;a 为电力信号基波分量的时变特性xg(t0)也即低频带限相量a(t0)的k阶导数,Δa为用泰勒级数表示电力信号基波分量的时变特性xg(t0)产生的误差;
(B2)数字信号处理器利用衰减指数函数 的泰勒级数近似表示电力信号x(t)中的衰减直流分量xd(t0), 其中,e为
(k)
自然对数的底;b、τ分别为衰减指数函数的幅值初始值和衰减时间常数;d 为衰减直流分量xd(t0)的k阶导数,Δd为用衰减指数函数 的泰勒级数表示电力信号x(t)中的衰减直流分量xd(t0)所产生的误差;
(B3)将B1步的电力信号基波分量的时变特性xg(t0)与常量旋转相量 相乘,得到电力信号基波相量的模型X′(t0), 将电力信号基波相量
的模型X′(t0)的实部加上B2步的衰减直流分量xd(t0),即得到考虑衰减直流分量的电力信号x(t)的动态模型x′(t0), 其中,j为虚
数单位,f0=50Hz为信号的基波频率,Re表示取相量实部;
然后利用欧拉公式,对A 步中的电力信号 x(t)进行离散化处理,
得 到 电 力 信 号 离 散 序 列x(n) 的 泰勒 级 数 形 式 的 数 学 模 型 x ′ (n):
其 中,ω0=2πf0/fs,ω0为
采样角速度,fs为采样频率;α(k)为电力信号基波分量的时变特性xg(t0)的k阶导数(k) (k) (k)
a 的离散化值, β 为衰减直流分量xd(t0)的k阶导数d 的离散化值,
(B4)将电力信号离散序列x(n)的数学模型x′(n)进行加窗傅里叶变换,得到电力信号离散序列x(n)在5个角速度ω处的相量测量预估计值X(ω)的数学模型X′(ω):
其中,N为加窗傅里叶变换的数据窗中的数据个数;h(n)为窗函数序列。
3.如权利要求1所述的利用频域信息滤除衰减直流的电力信号同步相量测量方法,其特征在于:所述的步骤B中求得电力信号基波分量时变特性的各阶导数值的具体方法为:
(B5)数字信号处理器将步骤A中得到的5个相量估计值X(ω)分别代入相量测量预估(k) (k)
计值X(ω)的数学模型X′(ω)中,得到5个含α 以及β 的方程;将这些方程联立组(k) (k)
成方程组,求出电力信号基波分量的时变特性xg(t0)的k阶导数a 的离散化值α ,再根(k)
据 算出电力信号基波分量的时变特性xg(t0)的各阶导数值a (t0)。
4.如权利要求1所述的利用频域信息滤除衰减直流的电力信号同步相量测量方法,其特征在于:所述的步骤B中相移运算得到当前时刻trep电网的电力信号的基波相量值的具体方法是:
(B6)当前时刻trep由数字信号处理器从GPS系统中获得,将当前时刻trep与基准时刻t0的差Δt代入B3步中电力信号基波相量的模型X′(t0)中,同时将B5步得到的电力信(k)
号基波分量的时变特性xg(t0)的各阶导数a (t0),代入B3步中电力信号基波相量的模型X′(t0)中,即得到当前时刻trep电网的电力信号的基波相量值

说明书全文

利用频域信息滤除衰减直流的电信号同步相量测量方法

技术领域

[0001] 本发明涉及电力系统中的电力信号同步相量测量方法。

背景技术

[0002] 随着全球电力市场化和电网区域互联的发展,电网的运行环境日益复杂,其安全稳定运行问题日益突出,迫切需要提高电网的动态安全监控能力。近年来,广域测量系统(Wide Area Measurement System,WAMS)作为一种新的电网动态监控系统,为电网动态安全监控提供了新的技术手段。WAMS中的数据采集和处理方式要求对采集数据的处理结果是带有精确时标的相量数据。因此,能同步采集电力信号的时间、相位和幅值的同步相量测量技术是实现WAMS的基础和核心。而同步相量测量技术的核心是相量估计方法的设计,估计方法的估计精度将直接影响到WAMS的应用效果。
[0003] 离散傅里叶变换(Discrete Fourier Transform,DFT)算法具有良好的谐波抑制特性和快速运算特性,在静态条件下具有较好应用价值,已被广泛应用到同步相量测量中。
[0004] 现有的DFT算法,在非同步采样的情况下存在幅值和频率误差,且随着其非同步性的增强,误差急剧增大,往往达不到实际应用的要求。当系统频率是随着系统配置参数以及状态变化而变化时,不能保证采样系统对实测信号一直保持同步采样。对于频率为f的正弦序列,它的频谱应该只是在f处有离散谱。但是,在利用DFT求它的频谱做了截短,结果使信号的频谱不只是在f处有离散谱,而是在以f为中心的频带范围内都有谱线出现,它们可以理解为是从f频率上泄露出去的;同时,对一函数实行采样,即是抽取采样点上的对应函数值。其效果如同透过栅栏的缝隙观看外景一样,只有落在缝隙前的少数景象被看到,其余景象均被栅栏挡住而视为零,这种现象视为栅栏效应。
[0005] 因此应用DFT进行相量测量时所引起的频谱泄露和栅栏现象可能会使得DFT算法产生较大误差,甚至得到一个不可用的结果。并且当电网处于动态过程中,电力信号的频率及幅值均随着时间变化,算法对信号动态特性的不完全表示也可能会加大DFT算法产生的误差。再者,当电力系统发生接地故障时,故障电流信号中可能包含衰减直流分量,这也会加大DFT算法产生的误差。

发明内容

[0006] 本发明的目的是提供利用频域信息滤除衰减直流的电力信号同步相量测量方法,该方法能在电网处于频率偏移等动态过程中时,更精确的实现电力信号的同步相量测量,并且能够减弱或消除信号中可能存在的衰减直流分量对相量测量精度的影响。
[0007] 本发明为实现其发明目的所采用的技术方案为:利用频域信息滤除衰减直流的电力信号同步相量测量方法,其步骤为:
[0008] A、数据采集与预估计处理
[0009] 电网中的电流互感器或电压互感器测量电网的电流信号或电压信号得到电网的电力信号x(t),采集电力信号x(t)得到电网的电力信号离散序列x(n),其中t为时间,n为采样时刻点;数字信号处理器利用加窗傅里叶变换对电力信号离散序列x(n)滤波,得到电力信号离散序列x(n)在5个速度ω处的相量测量预估计值X(ω);其中ω=2πf,f为角速度ω对应的频率,f=0,1,49,50,51Hz;
[0010] B、修正相量预估计值
[0011] 数字信号处理器对得到的5个相量测量预估计值X(ω),用考虑衰减直流分量的动态信号模型,求得电力信号基波分量的时变特性的各阶导数值,并进行相移运算得到当前时刻trep电网的电力信号的基波相量值
[0012] 与现有技术相比,本发明的有益效果是:
[0013] 本发明首先对电力信号离散序列进行加窗傅里叶变换,利用0,1,49,50,51Hz共5个频点(角速度)处的相量分量获得5个相量测量预估计值,其中0,1Hz频点处的相量测量预估计值含有丰富的衰减直流分量信息,而49,50,51Hz频点处的相量测量预估计值含有丰富的频率偏移等基波时变特性信息。再利用低频带限相量及泰勒级数合理表示信号基波分量的时变特性,利用衰减指数函数及泰勒级数合理表示信号的衰减直流分量,从而得到计算简单、模拟足够准确的考虑衰减直流分量的电力信号动态模型;进而利用该模型有效计算出在考虑衰减直流分量情况下的电力信号基波分量时变特性的各阶导数值;最后进行相移运算,得到修正后的当前时刻的精确相量值。求出的相量值既消除了衰减直流分量对相量计算的影响,又能准确地测出频率偏移等动态过程形成的电力信号基波分量的时变特性对相量值产生的偏移,其对电力信号基波相量值的测量更加精确、可靠。
[0014] 进一步,所述的步骤B中考虑衰减直流分量的动态信号模型用以下方法建立:
[0015] (B1)数字信号处理器利用低频带限相量a(t0)的泰勒级数近似表示电力信号基波分量的时变特性xg(t0), 其中,Δt为当前时刻trep与基准时刻t0的差,基准时刻t0为加窗傅里叶变换的数据窗中心时刻,k为(k)
泰勒级数的阶次,k=0,1,2;K为泰勒级数的最高阶次,取为2;a 为电力信号基波分量的时变特性xg(t0)也即低频带限相量a(t0)的k阶导数,Δa为用泰勒级数表示电力信号基波分量的时变特性xg(t0)产生的误差;
[0016] (B2)数字信号处理器利用衰减指数函数 的泰勒级数近似表示电力信号x(t)中的衰减直流分量xd(t0), 其中,e为自然对数的底;b、τ分别为衰减指数函数的幅值初始值和衰减时间常数;d(k)为衰减直流分量xd(t0)的k阶导数,Δd为用衰减指数函数 的泰勒级数表示电力信号x(t)中的衰减直流分量xd(t0)所产生的误差;
[0017] (B3)将B1步的电力信号基波分量的时变特性xg(t0)与常量旋转相量相乘,得到电力信号基波相量的模型X′(t0), 将电力信号基波相量的模型X′(t0)的实部加上B2步的衰减直流分量xd(t0),即得到考虑衰减直流分量的电力信号x(t)的动态模型x′(t0), 其中,
j为虚数单位,f0=50Hz为信号的基波频率,Re表示取相量实部;
[0018] 然后利用欧拉公式,对A步中的电力信号x(t)进行离散化处理,得到电力信号离散序列x(n) 的泰勒级数形式的数学模型x′ (n):
其中,ω0=2πf0/fs,ω0为采
(k) (k)
样角速度,fs为采样频率;α 为电力信号基波分量的时变特性xg(t0)的k阶导数a 的离(k) (k) k (k) (k) (k)
散化值,α =a /fs;β 为衰减直流分量xd(t0)的k阶导数d 的离散化值,β =(k) k
d /fs;
[0019] (B4)将电力信号离散序列x(n)的数学模型x′(n)进行加窗傅里叶变换,得到电力信号离散序列x(n)在5个角速度ω处的相量测量预估计值X(ω)的数学模型X′(ω):
[0020]
[0021] 其中,N为加窗傅里叶变换的数据窗中的数据个数;h(n)为窗函数序列。
[0022] 进一步,本发明的所述的步骤B中求得电力信号基波分量时变特性的各阶导数值的具体方法为:
[0023] (B5)数字信号处理器将步骤A中得到的5个相量估计值X(ω)分别代入相量测量(k) (k)预估计值X(ω)的数学模型X′(ω)中,得到5个含α 以及β 的方程;将这些方程联(k) (k)
立组成方程组,求出电力信号基波分量的时变特性xg(t0)的k阶导数a 的离散化值α ,(k) (k) k (k)
再根据α =a /fs,算出电力信号基波分量的时变特性xg(t0)的各阶导数值a (t0)。
[0024] 进一步,本发明所述的步骤B中相移运算得到当前时刻trep电网的电力信号的基波相量值 的具体方法是:
[0025] (B6)当前时刻trep由数字信号处理器从GPS系统中获得,将当前时刻trep与基准时刻t0的差Δt代入B3步中电力信号基波相量的模型X′(t0)中,同时将B5步得到的电力信号基波分量的时变特性xg(t0)的各阶导数a(k)(t0),代入B3步中电力信号基波相量的模型X′(t0)中,即得到当前时刻trep电网的电力信号的基波相量值
[0026]

具体实施方式

[0027] 下面结合具体实施方式对本发明做进一步的详细说明。
[0028] 实施例
[0029] 本发明的一种具体实施方式是:利用频域信息滤除衰减直流的电力信号同步相量测量方法,其步骤为:
[0030] A、数据采集与预估计处理
[0031] 电网中的电流互感器或电压互感器测量电网的电流信号或电压信号得到电网的电力信号x(t),采集电力信号x(t)得到电网的电力信号离散序列x(n),其中t为时间,n为采样时刻点;数字信号处理器利用加窗傅里叶变换对电力信号离散序列x(n)滤波,得到电力信号离散序列x(n)在5个角速度ω处的相量测量预估计值X(ω);其中ω=2πf,f为角速度ω对应的频率,f=0,1,49,50,51Hz;
[0032] B、修正相量预估计值
[0033] 数字信号处理器对得到的5个相量测量预估计值X(ω),用考虑衰减直流分量的动态信号模型,求得电力信号基波分量的时变特性的各阶导数值,并进行相移运算得到当前时刻trep电网的电力信号的基波相量值
[0034] 本例中所述的步骤B中考虑衰减直流分量的动态信号模型用以下方法建立:
[0035] (B1)数字信号处理器利用低频带限相量a(t0)的泰勒级数近似表示电力信号基波分量的时变特性xg(t0), 其中,Δt为当前时刻trep与基准时刻t0的差,基准时刻t0为加窗傅里叶变换的数据窗中心时刻,k为泰勒级数的阶次,k=0,1,2;K为泰勒级数的最高阶次,取为2;a(k)为电力信号基波分量的时变特性xg(t0)也即低频带限相量a(t0)的k阶导数,Δa为用泰勒级数表示电力信号基波分量的时变特性xg(t0)产生的误差;
[0036] (B2)数字信号处理器利用衰减指数函数 的泰勒级数近似表示电力信号x(t)中的衰减直流分量xd(t0), 其中,e为自然对数的底;b、τ分别为衰减指数函数的幅值初始值和衰减时间常数;d(k)为衰减直流分量xd(t0)的k阶导数,Δd为用衰减指数函数 的泰勒级数表示电力信号x(t)中的衰减直流分量xd(t0)所产生的误差;
[0037] (B3)将B1步的电力信号基波分量的时变特性xg(t0)与常量旋转相量相乘,得到电力信号基波相量的模型X′(t0), 将电力信号基波相量的模型X′(t0)的实部加上B2步的衰减直流分量xd(t0),即得到考虑衰减直流分量的电力信号x(t)的动态模型x′(t0), 其中,
j为虚数单位,f0=50Hz为信号的基波频率,Re表示取相量实部;
[0038] 然后利用欧拉公式,对A步中的电力信号x(t)进行离散化处理,得到电力信号离散序列x(n) 的泰勒级数形式的数学模型x′ (n):
其中,ω0=2πf0/fs,ω0为采
(k) (k)
样角速度,fs为采样频率;α 为电力信号基波分量的时变特性xg(t0)的k阶导数a 的离(k) (k) k (k) (k) (k)
散化值,α =a /fs;β 为衰减直流分量xd(t0)的k阶导数d 的离散化值,β =(k) k
d /fs;
[0039] (B4)将电力信号离散序列x(n)的数学模型x′(n)进行加窗傅里叶变换,得到电力信号离散序列x(n)在5个角速度ω处的相量测量预估计值X(ω)的数学模型X′(ω):
[0040]
[0041] 其中,N为加窗傅里叶变换的数据窗中的数据个数;h(n)为窗函数序列。
[0042] 本例中所述的步骤B中求得电力信号基波分量时变特性的各阶导数值的具体方法为:
[0043] (B5)数字信号处理器将步骤A中得到的5个相量估计值X(ω)分别代入相量测量(k) (k)预估计值X(ω)的数学模型X′(ω)中,得到5个含α 以及β 的方程;将这些方程联(k) (k)
立组成方程组,求出电力信号基波分量的时变特性xg(t0)的k阶导数a 的离散化值α ,(k) (k) k (k)
再根据α =a /fs,算出电力信号基波分量的时变特性xg(t0)的各阶导数值a (t0)。
[0044] 本例中所述的步骤B中相移运算得到当前时刻trep电网的电力信号的基波相量值的具体方法是:
[0045] (B6)当前时刻trep由数字信号处理器从GPS系统中获得,将当前时刻trep与基准时刻t0的差Δt代入B3步中电力信号基波相量的模型X′(t0)中,同时将B5步得到的电力信号基波分量的时变特性xg(t0)的各阶导数a(k)(t0),代入B3步中电力信号基波相量的模型X′(t0)中,即得到当前时刻trep电网的电力信号的基波相量值
[0046]
[0047] 仿真实验
[0048] 为验证本发明对处于动态条件下并且含有衰减直流分量的电力信号的相量测量精度,建立频率偏移情况下含有衰减直流分量的理想动态信号模型,分别用本发明方法和傅里叶算法对该信号进行测量,从而对比算法性能。仿真中,旋转调制相量的频率f0取为50Hz,以采样频率fs=2.4kHz对信号进行采样,采用矩形数据窗且其数据个数取为N=48。
依据本发明,对一个数据窗使用5个角速度ω=2πf,f=0,1,49,50,51Hz处的数据。
[0049] 仿真中以幅值误差和相角误差以及总相量误差(Total Vector Error,TVE)作为衡量算法性能的指标。
[0050] 仿真中利用一频率为非额定频率且含有衰减直流分量的信号来表示动态情况下
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈