CN105137198A - 一种基于Nuttall窗-五点变换FFT的介损测量新方法 - Google Patents
一种基于Nuttall窗-五点变换FFT的介损测量新方法 Download PDFInfo
- Publication number
- CN105137198A CN105137198A CN201510583456.XA CN201510583456A CN105137198A CN 105137198 A CN105137198 A CN 105137198A CN 201510583456 A CN201510583456 A CN 201510583456A CN 105137198 A CN105137198 A CN 105137198A
- Authority
- CN
- China
- Prior art keywords
- voltage
- dielectric loss
- fft
- sequence
- nuttall window
- 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.)
- Pending
Links
Landscapes
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明涉及一种基于Nuttall窗-五点变换FFT算法的介质损耗角高精度测量方法,该方法采用Nuttall窗对输入信号进行截断,以减小截断误差。然后利用FFT计算出截断后的信号相量序列。再利用五点变换对FFT输出序列进行加权平均以加大主频点外的谱线衰减速度,从而减小频谱泄漏效应。最后通过对加权后的新信号相量序列进行分析计算出介质损耗角。本发明能有效克服直流分量、基波频率、谐波、随机噪声对介质损耗角测量的影响,在随机噪声干扰下测量精度能达到0.5×10-5,在上述其它干扰下测量精度能达到0.5×10-6。
Description
技术领域
本发明涉及一种基于Nuttall窗-五点变换的介质损耗角测量方法。
背景技术
容性设备是电力***设备的重要组成部分。若其发生设备故障,则在很大程度上能影响到电力***的正常运行。而容性设备的设备故障绝大部分是由于绝缘问题引发的,因此有必要对其绝缘状态进行监测。介质损耗角δ是表征设备绝缘状态的重要参量之一,通过对设备的介质损耗角的监测,以及对其绝缘状态变化趋势的分析能对设备绝缘状态进行有效的判断。
对于介质损耗角监测的软件算法中主要以谐波分析法为代表。谐波分析法中又以FFT方法应用最为广泛,具有处理速度快,不受谐波干扰等优点。
FFT算法需要在同步采样的条件下,当采样频率与信号频率之间不为整数倍关系时由传统FFT算法得到的基波幅值相角就和其真实值存在较大的误差。实际电网频率总是在50Hz附近波动,传统的消除非同步采样误差的方法有以下几种:一是在***中加入锁相环实现频率跟踪从而达到同步采样的条件。锁相环实现简单,但是由于其响应较慢不能真正意义上实现完全同步采样。另一种则是通过对传统的FFT算法进行优化,以减小非同步采样对其计算精度的影响。
发明内容
本发明的目的提供一种基于Nuttall窗-五点变换FFT算法的介质损耗角高精度测量方法。
本发明采用以下技术方案实现:一种基于Nuttall窗-五点变换FFT的介损测量新方法,其特征在于:先在时域内采用Nuttall窗对电压、电流信号进行截断;通过快速傅里叶算法求出电压、电流的相量序列;对求得的电压、电流的相量序列采用五点加权变换,再利用变换得到的新电压、电流的相量序列计算出基波电压、电流的相位角,从而得到介质损耗角。
进一步的,包括以下具体计算步骤:步骤S1:将含有谐波、噪声以及直流分量的电压信号u(t)和电流信号i(t)以采样频率fs进行采样,得到离散时间信号:
其中random(n)为随机干扰,ku和ki为直流分量,H为谐波总次数,f为电压电流信号频率,δhu和δhi为电压电流信号初相角,Uh和Ih分别为第h次谐波的电压、电流幅值;
步骤S2:用Nuttall窗对公式(1)和(2)中的电压电流信号u(n)和i(n)进行截断得到新的信号序列:
uw(n)=w(n)u(n)(3)
iw(n)=w(n)i(n)(4)
其中,Nuttall窗w(n)的时域表达式为:
式中,a0=0.3635819,a1=0.4891775,a2=0.1365995,a3=0.0106411;
步骤S3:设Uw(n)、Iw(n)分别为电压电流信号uw(n)和iw(n)的FFT输出序列,则通过公式(6)对序列Uw(n)、Iw(n)进行五点加权变换得到新序列Uh5(n)、Ih5(n):
Xh5(n)=Xw(n)/4-[Xw(n-1)+Xw(n+1)]/6+[Xw(n-2)+Xw(n+2)]/24(6)
其中,n=1,2,…,N-1,N;
步骤S4:由步骤S3可以得到五点变换后的新的电压电流信号FFT输出序列Uh5(n)、Ih5(n),则介质损耗角计算公式为;
在本发明一实施例中,步骤S3中当n=1,2,N-1,N时,公式(6)中不存在的项用零代替。
与现有技术相比,本发明的有益效果是利用Nuttall窗减小非同步采样的误差,利用五点变换加快主频点以外的谱线的衰减速度,从而抑制频谱泄漏效应。提高了非同步采样下基于FFT的介质损耗角测量误差,测量精度高,硬件要求小。
附图说明
图1是介质损耗角测量算法流程图。
图2是Nuttall窗时域特性图。
图3是Nuttall窗频域特性图。
具体实施方式
一种基于Nuttall窗-五点变换FFT的介损测量新方法,先在时域内采用Nuttall窗对电压、电流信号进行截断;通过快速傅里叶算法求出电压、电流的相量序列;对求得的电压、电流的相量序列采用五点加权变换,再利用变换得到的新电压、电流的相量序列计算出基波电压、电流的相位角,从而得到介质损耗角。具体计算流程步骤参见图1。
主要计算步骤如下:
步骤S1:将含有谐波,噪声以及直流分量的电压信号u(t)和电流信号i(t)以采样频率fs进行采样,得到离散时间信号:
其中random(n)为随机干扰,ku和ki为直流分量,H为谐波总次数,f为电压电流信号频率,δhu和δhi为电压电流信号初相角,Uh和Ih为第h次谐波的幅值。。
步骤S2:步骤S2:用Nuttall窗对公式(1)和(2)中的电压电流信号u(n)和i(n)进行截断得到新的信号序列:
uw(n)=w(n)u(n)(10)
iw(n)=w(n)i(n)(11)
其中,Nuttall窗w(n)的时域表达式为:
式中,a0=0.3635819,a1=0.4891775,a2=0.1365995,a3=0.0106411;
步骤S3:设Uw(n)、Iw(n)分别为电压电流信号uw(n)和iw(n)的FFT输出序列,则
通过公式(6)对序列Uw(n)、Iw(n)进行五点加权变换得到新序列Uh5(n)、Ih5(n):
Xh5(n)=Xw(n)/4-[Xw(n-1)+Xw(n+1)]/6+[Xw(n-2)+Xw(n+2)]/24(13)
其中,n=1,2,…,N-1,N;特别地,当n=1,2,N-1,N时,公式(6)中不存在的项用零代替;
步骤S4:由步骤S3可以得到五点变换后的新的电压电流信号FFT输出序列Uh5(n)、Ih5(n),则介质损耗角计算公式为;
本发明一具体实施例中的容性绝缘设备的等效电路模型采用电阻电容RC并联等效电路模型,其中电容C为100uF,电阻R为5000Ω,介损真值计算公式为
(1)利用离散采样与AD转换
利用DSP的AD功能将所测得的电压电流模拟信号转换成数字信号,采样频率fs=500Hz。
(2)利用Nuttall窗对信号进行截断
选用合适长度的Nuttall窗对电压电流信号进行截断,其时域特性与频域特性如图2与图3所示。
(3)对截断后的数据进行傅里叶变换
(4)对傅里叶变换后得到的相量序列进行五点变换。
利用变换公式(6)对傅里叶变换后的结果进行五点加权变换。
(5)计算介质损耗角
根据公式(7),利用上一步得到的相量序列计算出介质损耗角。
在本实例中,当基波频率从49.5Hz变化到50.5Hz时,介质损耗角测量结果如表1所示。
表1频率波动时介质损耗角测量结果
在本实例中,当基波频率为50Hz,三次谐波含量变化时,介质损耗角测量结果如表2所示。
表2三次谐波含量变化时介质损耗角测量结果
在本实例中,当基波频率为50Hz,随机噪声含量变化时,介质损耗角测量结果如表3所示。
表3随机噪声含量变化时介质损耗角测量结果
采用本发明的技术方案从而抑制频谱泄漏效应。提高了非同步采样下基于FFT的介质损耗角测量误差,测量精度高,硬件要求小。
以上所述仅为本发明的较佳实施例,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本发明的涵盖范围。因此本发明的保护范围应该以权利要求的保护范围为准。
Claims (3)
1.一种基于Nuttall窗-五点变换FFT的介损测量新方法,其特征在于:先在时域内采用Nuttall窗对电压、电流信号进行截断;通过快速傅里叶算法求出电压、电流的相量序列;对求得的电压、电流的相量序列采用五点加权变换,再利用变换得到的新电压、电流的相量序列计算出基波电压、电流的相位角,从而得到介质损耗角。
2.根据权利要求1所述的基于Nuttall窗-五点变换FFT的介损测量新方法,其特征在于:包括以下具体计算步骤:
步骤S1:将含有谐波、噪声以及直流分量的电压信号u(t)和电流信号i(t)以采样频率fs进行采样,得到离散时间信号:
其中random(n)为随机干扰,ku和ki为直流分量,H为谐波总次数,f为电压电流信号频率,δhu和δhi为电压电流信号初相角,Uh和Ih分别为第h次谐波的电压、电流幅值;
步骤S2:用Nuttall窗对公式(1)和(2)中的电压电流信号u(n)和i(n)进行截断得到新的信号序列:
uw(n)=w(n)u(n)(3)
iw(n)=w(n)i(n)(4)
其中,Nuttall窗w(n)的时域表达式为:
式中,a0=0.3635819,a1=0.4891775,a2=0.1365995,a3=0.0106411;
步骤S3:设Uw(n)、Iw(n)分别为电压电流信号uw(n)和iw(n)的FFT输出序列,则通过公式(6)对序列Uw(n)、Iw(n)进行五点加权变换得到新序列Uh5(n)、Ih5(n):
Xh5(n)=Xw(n)/4-[Xw(n-1)+Xw(n+1)]/6+[Xw(n-2)+Xw(n+2)]/24
(6)
其中,n=1,2,…,N-1,N;
步骤S4:由步骤S3可以得到五点变换后的新的电压电流信号FFT输出序列Uh5(n)、Ih5(n),则介质损耗角计算公式为;
3.根据权利要求2所述的基于Nuttall窗-五点变换FFT的介损测量新方法,其特征在于:步骤S3中当n=1,2,N-1,N时,公式(6)中不存在的项用零代替。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510583456.XA CN105137198A (zh) | 2015-09-15 | 2015-09-15 | 一种基于Nuttall窗-五点变换FFT的介损测量新方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510583456.XA CN105137198A (zh) | 2015-09-15 | 2015-09-15 | 一种基于Nuttall窗-五点变换FFT的介损测量新方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105137198A true CN105137198A (zh) | 2015-12-09 |
Family
ID=54722614
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510583456.XA Pending CN105137198A (zh) | 2015-09-15 | 2015-09-15 | 一种基于Nuttall窗-五点变换FFT的介损测量新方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105137198A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104360170A (zh) * | 2014-10-31 | 2015-02-18 | 国家电网公司 | 一种基于迭代匹配追踪容性设备介损角计算方法 |
CN105588985A (zh) * | 2016-03-18 | 2016-05-18 | 国网上海市电力公司 | 一种基于相位因子变换插值的介损在线精确测量方法 |
CN106019101A (zh) * | 2016-06-24 | 2016-10-12 | 兴义供电局 | 一种电力设备绝缘状态评估方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101113995A (zh) * | 2007-08-29 | 2008-01-30 | 湖南大学 | 基于Nuttall窗双峰插值FFT的基波与谐波检测方法 |
CN102253282A (zh) * | 2011-06-10 | 2011-11-23 | 中国矿业大学 | 纳托尔窗函数连续频谱内插电力谐波参数获取方法 |
CN104391178A (zh) * | 2014-12-05 | 2015-03-04 | 国家电网公司 | 一种基于Nuttall窗的时移相位差稳态谐波信号校正方法 |
-
2015
- 2015-09-15 CN CN201510583456.XA patent/CN105137198A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101113995A (zh) * | 2007-08-29 | 2008-01-30 | 湖南大学 | 基于Nuttall窗双峰插值FFT的基波与谐波检测方法 |
CN102253282A (zh) * | 2011-06-10 | 2011-11-23 | 中国矿业大学 | 纳托尔窗函数连续频谱内插电力谐波参数获取方法 |
CN104391178A (zh) * | 2014-12-05 | 2015-03-04 | 国家电网公司 | 一种基于Nuttall窗的时移相位差稳态谐波信号校正方法 |
Non-Patent Citations (3)
Title |
---|
万全等: "基于Nuttall 窗频谱校正的介质损耗因数测量", 《仪器仪表学报》 * |
孙鹏等: "五点加权FFT介质损耗角测量算法的研究", 《高压电器》 * |
陈国志等: "基于Nuttall窗插值FFT的谐波分析方法", 《电力自动化设备》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104360170A (zh) * | 2014-10-31 | 2015-02-18 | 国家电网公司 | 一种基于迭代匹配追踪容性设备介损角计算方法 |
CN104360170B (zh) * | 2014-10-31 | 2017-11-03 | 国家电网公司 | 一种基于迭代匹配追踪容性设备介损角计算方法 |
CN105588985A (zh) * | 2016-03-18 | 2016-05-18 | 国网上海市电力公司 | 一种基于相位因子变换插值的介损在线精确测量方法 |
CN106019101A (zh) * | 2016-06-24 | 2016-10-12 | 兴义供电局 | 一种电力设备绝缘状态评估方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103576002B (zh) | 一种容性绝缘设备介质损耗角的计算方法 | |
CN103487652B (zh) | 一种频率自适应实时分次谐波检测方法 | |
CN102435844B (zh) | 一种频率无关的正弦信号相量计算方法 | |
CN103454497B (zh) | 基于改进加窗离散傅立叶变换的相位差测量方法 | |
CN103869162B (zh) | 一种基于时域准同步的动态信号相量测量方法 | |
CN103308766A (zh) | 一种基于凯撒自卷积窗双谱线插值fft谐波分析方法及其装置 | |
WO2016138764A1 (zh) | 一种提升pmu同步相量测量精度的方法 | |
CN103983849B (zh) | 一种实时高精度的电力谐波分析方法 | |
Wen et al. | Hanning self-convolution window and its application to harmonic analysis | |
CN101441233A (zh) | 基于Kaiser窗双谱线插值FFT的基波与谐波检测方法 | |
CN105487034A (zh) | 一种0.05级电子式互感器校验方法及*** | |
CN104407198A (zh) | 一种用于dvr装置中检测sag信号的方法和*** | |
CN105137198A (zh) | 一种基于Nuttall窗-五点变换FFT的介损测量新方法 | |
CN113253052A (zh) | 一种基于改进smmg的高压直流输电线路故障测距方法 | |
Kumar et al. | A novel sub-cycle-based method for estimation of decaying DC component and fundamental phasor | |
CN110967658A (zh) | 一种基于数字微差法的模拟量输入合并单元校验仪溯源的方法 | |
CN103543331B (zh) | 一种计算电信号谐波和间谐波的方法 | |
Ren et al. | Use of recursive wavelet transform for estimating power system frequency and phasors | |
Berdin et al. | Estimating the instantaneous values of the state parameters during electromechanical transients | |
CN103487665B (zh) | 测量高压断路器均压电容试验方法 | |
CN109444539B (zh) | 一种基于克拉克变换的同步相量测量方法 | |
CN110568309A (zh) | 一种滤波器、同步相量测量***及方法 | |
Zhang et al. | Frequency shifting and filtering algorithm for power system harmonic estimation | |
CN105467209B (zh) | 一种新的金属氧化物避雷器泄漏电流分析方法 | |
CN106053936B (zh) | 一种获取电学信号瞬时频率的方法及*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20151209 |
|
RJ01 | Rejection of invention patent application after publication |