[发明专利]一种风力发电机组电压闪变的测量方法在审

专利信息
申请号: 201711329176.1 申请日: 2017-12-13
公开(公告)号: CN108169540A 公开(公告)日: 2018-06-15
发明(设计)人: 孟克其劳;王旭;张占强 申请(专利权)人: 内蒙古工业大学
主分类号: G01R19/04 分类号: G01R19/04;G01R23/02;G01R23/165
代理公司: 西安铭泽知识产权代理事务所(普通合伙) 61223 代理人: 韩晓娟
地址: 010080 内蒙古自治区呼*** 国省代码: 内蒙古;15
权利要求书: 查看更多 说明书: 查看更多
摘要: 发明涉及风力发电机组电能质量测量技术领域,公开了一种风力发电机组电压闪变的测量方法,包括以下步骤:S1、测量风力发电机组的输出端,得到电压和电流的时间序列;S2、在虚拟电网中模拟电压波动,S3、对虚拟电网的输出信号进行平方;S4、通过Nuttall窗函数对平方后的输出信号进行截断;S5、利用FFT进行滤波;S6、对只包含直流分量和低频闪变包络的信号通过经验模态分解筛分出频率最高的闪变包络信号;S7、对每个IMF波动信号进行HHT变换;S8、对多个闪变电压进行频谱分析,获得包含多个频率成分的短时闪变严重程度,这种测量方法,能够准确测量出含有多个频率成分的闪变电压的短时闪变严重程度。 1
搜索关键词: 闪变 风力发电机组 测量 电压闪变 输出信号 虚拟电网 经验模态分解 包络信号 波动信号 模拟电压 频谱分析 时间序列 信号通过 直流分量 质量测量 准确测量 输出端 包络 截断 滤波
【主权项】:
1.一种风力发电机组电压闪变的测量方法,其特征在于,包括以下步骤:

S1、测量风力发电机组的输出端,得到电压和电流的时间序列;

S2、对电压和电流的时间序列,在没有其他波动源的虚拟电网中模拟电压波动;

S3、通过平方器对虚拟电网的输出信号进行平方;

S4、通过Nuttall窗函数对平方后的输出信号进行截断;

S5、利用FFT对Nuttall窗函数截断后的信号中高于35Hz的频率分量进行滤波,获得只包含直流分量和低频闪变包络的信号;

S6、对只包含直流分量和低频闪变包络的信号通过经验模态分解EMD减去过零点波动的方法筛分出频率最高的闪变包络信号,获得频率由高到底依次排列的多个IMF波动信号;

S7、对每个IMF波动信号进行HHT变换,得到多个闪变电压;

S8、对多个闪变电压进行频谱分析,获得包含多个频率成分的闪变电压的短时闪变严重程度。

2.如权利要求1所述的风力发电机组电压闪变的测量方法,其特征在于,所述步骤S2中的虚拟电网由理想电压源u0(t)、电流源、电感和电阻串联组成,风力发电机组用电流源表示,im(t)为线电流瞬时值,其中模拟电压的瞬时值ufic为:

式(1)中Rfic为虚拟电网的电阻,Lfic为虚拟电网的电感;

理想电压源u0(t)应满足理想电压的如下特性:理想电压不应有任何波动,且电压闪变为零;u0(t)应与测量电压的基波具有相同的电气角αm(t);

为满足以上条件u0(t)定义为:

式(2)中,un为电网额定电压的有效值,测量电压的基波的电气角αm(t)为:

式(3)中,f(t)为电网频率,随时间变化;t为开始记录的起始时间;α0为t=0时的电气角;

机组电压和电流时间序列通过虚拟电网,得到电压波动完全来自于风力发电机组的模拟电压的瞬时值ufic

3.如权利要求1所述的风力发电机组电压闪变的测量方法,其特征在于,所述步骤S4中Nuttall窗函数的时域表达式为:

式(4)中M为四项最小旁瓣Nuttall窗的项数,bm需满足式(5)

使用四项最小旁瓣Nuttall窗对u2(t)进行截断,即时间域相乘,频率域卷积,得到z(t)

z(t)=ufic2WN(n)                (6)

对四项最小旁瓣Nuttall窗截断后的信号z(t)进行傅里叶变换,得到其频谱F(ω)

式(7)中在频率域去掉F(ω)的大于35Hz高频分量,得到F1(ω),对F1(ω)求傅里叶逆变换得到只包含直流和低频闪变包络的信号z1(t)

式(8)中z1(t)仅包含直流分量和闪变包络。

4.如权利要求1所述的风力发电机组电压闪变的测量方法,其特征在于,所述步骤S6中所述经验模态分解EMD通过减去过零点波动的方法筛分出频率最高的成分,依照此方法,IMF根据频率从高到低依次排列,每个IMF代表了信号包含的一个波动,EMD过程如下:

1)令只包含直流和低频闪变包络的信号z1(t)=x(t),找到信号x(t)的极值点,使用三次样条差值的方法得到信号的上包络v1(t)和下包络v(t);

2)均值包络m(t)可以按照式(9)得到:

m(t)=[v1(t)+v(t)]/2                      (9)

3)由式(10)得到h1(t):

h1(t)=x(t)‑m(t)                         (10)

4)如果由式(10)得到的h1(t)不满足分量终止条件,则将h1(t)迭代为新的x(t),重复1)、2)和3)步,直到h1(t)满足分解终止条件时,迭代停止,此时得到的一个IMF记为c1(t),均值包络m(t)和减去均值包络m(t)之后的信号h(t)继续按照式(11)迭代:

下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于内蒙古工业大学,未经内蒙古工业大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201711329176.1/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

tel code back_top