CN110244123A - 一种基于线变关系识别装置的降阶广义积分器谐振点检测方法 - Google Patents

一种基于线变关系识别装置的降阶广义积分器谐振点检测方法 Download PDF

Info

Publication number
CN110244123A
CN110244123A CN201910549603.XA CN201910549603A CN110244123A CN 110244123 A CN110244123 A CN 110244123A CN 201910549603 A CN201910549603 A CN 201910549603A CN 110244123 A CN110244123 A CN 110244123A
Authority
CN
China
Prior art keywords
frequency
current
voltage
depression
harmonic
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
Application number
CN201910549603.XA
Other languages
English (en)
Other versions
CN110244123B (zh
Inventor
喻旻
余仁山
尹怀强
曹乾磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yingtan Power Supply Branch Jiangxi Electric Power Co Ltd
Qingdao Topscomm Communication Co Ltd
Original Assignee
Yingtan Power Supply Branch Jiangxi Electric Power Co Ltd
Qingdao Topscomm Communication Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Yingtan Power Supply Branch Jiangxi Electric Power Co Ltd, Qingdao Topscomm Communication Co Ltd filed Critical Yingtan Power Supply Branch Jiangxi Electric Power Co Ltd
Priority to CN201910549603.XA priority Critical patent/CN110244123B/zh
Publication of CN110244123A publication Critical patent/CN110244123A/zh
Application granted granted Critical
Publication of CN110244123B publication Critical patent/CN110244123B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16Spectrum analysis; Fourier analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/088Aspects of digital computing

Landscapes

  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Measurement Of Resistance Or Impedance (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种基于线变关系识别装置的降阶广义积分器谐振点检测方法,该方法通过线变关系识别设备主动向电网发送指定频次的载波信号,首先,采集线变关系识别设备接入点处的电网电压和电流信号;然后,根据提取到的电压和电流信号数据,利用降阶广义积分器计算公式,提取注入频次的谐波电压和谐波电流信号;进而,利用余弦相似度公式,计算所提取的注入频次谐波信号的余弦相似系数;最后,通过计算的相似系数判断谐振频率范围,并通过更改注入频率,重复上述步骤进一步缩小范围,该方案计算量小,能够实时检测电网谐振频率范围。

Description

一种基于线变关系识别装置的降阶广义积分器谐振点检测 方法
技术领域
本发明涉及公开了一种线变关系识别装置的谐振点检测方案,具体是一种基于降阶广义积分器的谐振频率检测方法。
背景技术
随着电力电子技术的发展以及大量非线性装置的普及应用,电网电压谐波问题越来越严重,这将引发电网各种故障问题,因此谐波的治理问题受到广泛的关注。目前治理谐波的装置主要有无源滤波装置和有源电力滤波装置两种,其中无源滤波器结构简单,但是灵活性差,易与电网产生谐振,效果不理想。而有源电力滤波装置工作与电网参数无关,并且能够同时补偿多个谐波频率,因此应用更加的广泛。
有源电力滤波装置在补偿容性负载时,若补偿频率在电网谐振频率附近,这将导致电网电压产生畸变,过大的谐振电流会引起设备的损害。这类设备在工作前需要事先检测当前电网的谐振点,因此,对于电网谐振点的实时在线检测方案研究具有实际意义。
发明内容
本发明的目的是提供一种解决电网谐振频率实时在线检测的方案,该方法基于线变识别设备进行指定频次载波电流的注入,采集电网的电压和电流,利用降阶广义积分器进行该注入频次信号的提取,并根据余弦相似度公式计算相似系数,最后利用相似系数的幅值进行谐振频率范围的判断。
为实现上述目的,本发明将采用的技术方案如下所示:
一种基于线变关系识别装置的降阶广义积分器谐振点检测方法,包括以下步骤,
步骤一,将线变关系识别设备接入电网低压侧,向电网发送指定频次和幅值的载波电流信号;
步骤二,利用线变关系识别设备采集接入点的电压和电流信号u,i,利用降阶广义积分器提取采集到电压电流信号的注入频次谐波电压、电流;
步骤三,根据提取出的注入频次谐波电压、电流信号,利用余弦相似度公式,计算该频次的谐波电压和谐波电流的相似系数;
步骤四,根据计算的该频次谐波电压和谐波电流的相似系数,判断电网谐振频率范围,更改线变关系识别装置的发送频率,重复上述步骤。
所述的指定频次和幅值的载波电流信号,幅值范围为5A~15A,频率范围为 500Hz~5KHz;
所述的降阶广义积分器提取指定频次谐波电压、电流原理如下所示:
二阶广义积分器传递函数可以表示为:
其中K表示为传递函数的增益,wn表示为第n次谐波角频率。
根据上式进而可得降阶广义积分器的传递函数如下:
那么,基于降阶广义积分器的谐波提取公式可以表示为:
In(s)=Gr(s)×I(s)
Un(s)=Gr(s)×U(s)
其中,I(s)采集电流的傅里叶变换,U(s)为采集电压的傅里叶变换;In(s)为第n次谐波电流的傅里叶变换,Un(s)为第n次谐波电压的傅里叶变换。
所述的计算提取谐波电压和电流余弦相似系数的原理如下所示:
余弦相关系数公式如下:
其中,N为采集的电压和电流数据点数,ρn为第n次谐波电压和电流的相似系数。
所述的电网谐振频率范围判断方法,其原理如下所示:
根据计算的注入频次电压和电流余弦相关系数ρn判断谐振频率范围,其中0≤ρn≤1,理论上当ρn=1时,所注入的载波频率即为当前时刻电网谐振频率,并且ρn越小,所注入的载波频率越远离电网谐振频率。
本发明的有益效果:提供了一种基于线变关系识别设备的降阶广义积分器谐振频率检测方案,本发明是一种实时在线的谐波频率检测方案,通过线变关系识别设备注入指定频次电流信号,在线采集和提取所注入频次的谐波电压和电流信号,所用的降阶广义积分提取方案与传统的二阶广义积分方案相比,阶次更低,计算量小,解决了谐波信号的在线提取和谐振点的实时检测问题。
附图说明
图1为本发明基于线变关系识别设备的谐振频率检测原理框图;
图2为本发明线变关系装置注入11次谐波时采集的低压侧A相电压、电流波形;
图3为本发明线变关系装置注入22次谐波时采集的低压侧A相电压、电流波形;
图4为本发明线变关系装置注入25次谐波时采集的低压侧A相电压、电流波形;
图5为本发明线变关系装置注入59次谐波时采集的低压侧A相电压、电流波形;
图6为本发明提取的低压侧A相电压、电流的11次谐波波形;
图7为本发明提取的低压侧A相电压、电流的22次谐波波形;
图8为本发明提取的低压侧A相电压、电流的25次谐波波形;
图9为本发明提取的低压侧A相电压、电流的59次谐波波形;
具体实施方式
下面结合说明书附图1-9,对本发明具体实施方式进行完整、详细的描述。
如图1所示为基于线变关系识别设备谐振频率检测原理框图,该方法通过采集电网的电压和电流,利用降阶广义积分器进行该注入频次信号的提取,并根据余弦相似度公式计算相似系数,最后利用相似系数的幅值进行谐振频率范围的判断,实现步骤如下所示:
步骤一,将线变关系识别设备接入电网低压侧,向电网发送指定频次为11次,22次, 25次,59次,幅值为10A的载波电流信号;
步骤二,利用线变关系识别设备采集接入点的电压和电流信号u,i,利用降阶广义积分器提取采集到电压、电流信号的注入频次谐波电压电流;
降阶广义积分器提取指定频次谐波电压电流原理如下所示:
二阶广义积分器传递函数可以表示为:
其中K表示为传递函数的增益,wn表示为第11次,22次,25次,59次谐波角频率。根据上式进而可得降阶后广义积分器的传递函数如下:
那么,基于降阶广义积分器的谐波提取公式可以表示为:
In(s)=Gr(s)×I(s)
Un(s)=Gr(s)×U(s)
其中,I(s)采集电流的傅里叶变换,U(s)为采集电压的傅里叶变换;In(s)为第11次, 22次,25次,59次谐波电流的傅里叶变换,Un(s)为第11次,22次,25次,59次谐波电压的傅里叶变换。
步骤三,根据提取出的注入频次谐波电压电流信号,利用余弦相似度公式,计算该频次的谐波电压和谐波电流的相似系数,
余弦相关系数公式如下:
其中,N为采集的电压和电流数据点数,ρn为第11次,22次,25次,59次的谐波电压和谐波电流相似系数,分别为0.4868,0.8998,0.8617,0.0649。
步骤四,根据计算的该频次谐波电压和谐波电流的相似系数,判断电网谐振频率范围,从步骤三可以看出第22次谐波附近电压和电流相似系数最高,远离22次时,相似系数逐渐减小,其中第59次谐波时电压和电流相似系数最小,因此可以判断电网谐振频率在1100Hz附近。

Claims (4)

1.一种基于线变关系识别装置的降阶广义积分器谐振点检测方法,其特征在于:利用线变识别设备进行指定频次载波的注入,采集电网的电压和电流,利用降阶广义积分器进行该注入频次信号的提取,并根据余弦相似度公式计算相似系数,最后根据相似系数的幅值进行谐振频率范围的判断;具体步骤如下:
步骤一,将线变关系识别设备接入电网低压侧,测试发送有效幅值范围为5A~15A,频率范围为500Hz~5KHz的载波电流信号;
步骤二,利用线变关系识别设备采集接入点的电压和电流信号u,i,利用降阶广义积分器提取采集到电压、电流信号的注入频次谐波电压电流;
步骤三,根据提取出的注入频次谐波电压电流信号,利用余弦相似度公式,计算该频次谐波电压和谐波电流的相似系数;
步骤四,根据计算的该频次谐波电压和谐波电流相似系数,判断电网谐振频率范围,更改线变关系识别装置的发送电流频率,重复上述步骤。
2.根据权利要求1所述的基于线变关系识别装置的降阶广义积分器谐振点检测方法,其特征在于:
降阶广义积分器原理如下:
二阶广义积分器传递函数可以表示为:
其中K表示为传递函数的增益,wn表示为第n次谐波角频率。
根据上式进而可得降阶广义积分器的传递函数如下:
那么,基于降阶广义积分器的谐波提取公式可以表示为:
In(s)=Gr(s)×I(s)
Un(s)=Gr(s)×U(s)
其中,I(s)采集电流的傅里叶变换,U(s)为采集电压的傅里叶变换;In(s)为第n次谐波电流的傅里叶变换,Un(s)为第n次谐波电压的傅里叶变换。
3.根据权利要求1所述的提取的注入频次谐波电压电流信号余弦相似系数,其特征在于:余弦相关系数公式如下:
其中,N为采集的电压和电流数据点数,ρn为第n次谐波电压和电流的相似系数。
4.根据权利要求1所述的电网谐振频率范围判断方法,其特征在于:根据步骤三计算的注入频次电压和电流余弦相关系数ρn判断谐振频率范围,其中0≤ρn≤1,理论上当ρn=1时,当前时刻所注入的载波频率即为电网谐振频率,并且ρn越小,所注入的载波频率离电网谐振频率越远。
CN201910549603.XA 2019-06-24 2019-06-24 一种基于线变关系识别装置的降阶广义积分器谐振点检测方法 Active CN110244123B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910549603.XA CN110244123B (zh) 2019-06-24 2019-06-24 一种基于线变关系识别装置的降阶广义积分器谐振点检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910549603.XA CN110244123B (zh) 2019-06-24 2019-06-24 一种基于线变关系识别装置的降阶广义积分器谐振点检测方法

Publications (2)

Publication Number Publication Date
CN110244123A true CN110244123A (zh) 2019-09-17
CN110244123B CN110244123B (zh) 2021-07-02

Family

ID=67889050

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910549603.XA Active CN110244123B (zh) 2019-06-24 2019-06-24 一种基于线变关系识别装置的降阶广义积分器谐振点检测方法

Country Status (1)

Country Link
CN (1) CN110244123B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110994607A (zh) * 2019-12-13 2020-04-10 国网吉林省电力有限公司松原供电公司 线变关系识别系统
CN111693771A (zh) * 2020-06-22 2020-09-22 南方电网科学研究院有限责任公司 一种电网谐振点检测装置及方法
CN111834952A (zh) * 2019-09-21 2020-10-27 青岛鼎信通讯股份有限公司 一种故障指示器安装拆卸工具
CN112217213A (zh) * 2020-10-24 2021-01-12 青岛鼎信通讯股份有限公司 一种基于apf装置的谐振点实时计算方法
CN112595896A (zh) * 2020-11-28 2021-04-02 国网河南省电力公司电力科学研究院 一种对地电容参数检测方法
CN113612230A (zh) * 2021-01-05 2021-11-05 青岛鼎信通讯股份有限公司 一种基于线变关系识别设备的低压电网谐波电流补偿方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040133096A1 (en) * 2003-01-07 2004-07-08 Yo Taniguchi Inspection apparatus using nuclear magnetic resonance
CN106019026A (zh) * 2016-07-20 2016-10-12 东南大学 一种基于波形匹配的谐波电压责任的划分方法
CN106786780A (zh) * 2017-03-02 2017-05-31 江苏大学 一种基于虚拟同步发电机的并网控制方法及系统
CN107666151A (zh) * 2016-07-22 2018-02-06 郭琳 一种基于静止坐标系的并网逆变器直接功率控制
CN108471124A (zh) * 2018-04-24 2018-08-31 华北电力大学(保定) 光伏逆变器接入交直流混合微电网的谐振抑制方法
CN109067193A (zh) * 2018-08-17 2018-12-21 燕山大学 一种级联型电力电子变压器及其不平衡补偿控制方法
CN109308571A (zh) * 2018-08-29 2019-02-05 华北电力科学研究院有限责任公司 配电线路线变关系检测方法
CN110007141A (zh) * 2019-04-26 2019-07-12 青岛鼎信通讯股份有限公司 一种基于电压电流谐波相似度的谐振点检测方法
CN111693771A (zh) * 2020-06-22 2020-09-22 南方电网科学研究院有限责任公司 一种电网谐振点检测装置及方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040133096A1 (en) * 2003-01-07 2004-07-08 Yo Taniguchi Inspection apparatus using nuclear magnetic resonance
CN106019026A (zh) * 2016-07-20 2016-10-12 东南大学 一种基于波形匹配的谐波电压责任的划分方法
CN107666151A (zh) * 2016-07-22 2018-02-06 郭琳 一种基于静止坐标系的并网逆变器直接功率控制
CN106786780A (zh) * 2017-03-02 2017-05-31 江苏大学 一种基于虚拟同步发电机的并网控制方法及系统
CN108471124A (zh) * 2018-04-24 2018-08-31 华北电力大学(保定) 光伏逆变器接入交直流混合微电网的谐振抑制方法
CN109067193A (zh) * 2018-08-17 2018-12-21 燕山大学 一种级联型电力电子变压器及其不平衡补偿控制方法
CN109308571A (zh) * 2018-08-29 2019-02-05 华北电力科学研究院有限责任公司 配电线路线变关系检测方法
CN110007141A (zh) * 2019-04-26 2019-07-12 青岛鼎信通讯股份有限公司 一种基于电压电流谐波相似度的谐振点检测方法
CN111693771A (zh) * 2020-06-22 2020-09-22 南方电网科学研究院有限责任公司 一种电网谐振点检测装置及方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
钟建伟,等: "基于降阶广义积分器的三相并网逆变器锁频环技术", 《农村电气化》 *
顾伟,等: "基于波形匹配的谐波责任划分方法", 《电力系统自动化》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111834952A (zh) * 2019-09-21 2020-10-27 青岛鼎信通讯股份有限公司 一种故障指示器安装拆卸工具
CN111834952B (zh) * 2019-09-21 2021-10-01 青岛鼎信通讯股份有限公司 一种故障指示器安装拆卸工具
CN110994607A (zh) * 2019-12-13 2020-04-10 国网吉林省电力有限公司松原供电公司 线变关系识别系统
CN111693771A (zh) * 2020-06-22 2020-09-22 南方电网科学研究院有限责任公司 一种电网谐振点检测装置及方法
CN112217213A (zh) * 2020-10-24 2021-01-12 青岛鼎信通讯股份有限公司 一种基于apf装置的谐振点实时计算方法
CN112595896A (zh) * 2020-11-28 2021-04-02 国网河南省电力公司电力科学研究院 一种对地电容参数检测方法
CN113612230A (zh) * 2021-01-05 2021-11-05 青岛鼎信通讯股份有限公司 一种基于线变关系识别设备的低压电网谐波电流补偿方法

Also Published As

Publication number Publication date
CN110244123B (zh) 2021-07-02

Similar Documents

Publication Publication Date Title
CN110244123A (zh) 一种基于线变关系识别装置的降阶广义积分器谐振点检测方法
CN110426604B (zh) 一种谐振接地系统单相接地故障选线方法
CN111077375B (zh) 基于频域减的背景谐波影响下电网阻抗在线辨识方法
CN111693771A (zh) 一种电网谐振点检测装置及方法
CN109783768A (zh) 基于波形松弛迭代的时域多导体传输线电磁脉冲响应快速建模方法
CN112067948B (zh) 一种配电网单相接地故障的故障选线方法及系统
CN109917478A (zh) 一种提高采集信号幅度的磁共振探水系统及方法
CN111551821A (zh) 一种配电网接地故障辨识方法、装置及设备
CN108268837A (zh) 基于小波熵和混沌特性的辐射源指纹特征提取方法
CN110112776A (zh) 考虑电网背景谐波的并网逆变器电网阻抗辨识方法
CN103823158B (zh) 采用不变矩的谐振接地系统故障选线方法
Ye et al. Single pole‐to‐ground fault location method for mmc‐hvdc system using wavelet decomposition and dbn
CN118311376A (zh) 配电网故障选线方法及装置、介质、设备
CN103913632A (zh) 一种基于振荡分量功率方向的低频同步振荡源识别方法
CN112217213B (zh) 一种基于apf装置的谐振点实时计算方法
CN213957528U (zh) 一种智能网络型故障定位系统
CN114034973A (zh) 一种配电线路接地故障的故障区域识别方法及其装置和系统
CN110556799B (zh) 适用于逆变型电源场站送出线路的方向元件设计方法
CN115616330B (zh) 一种基于波形相似度的输电线路多重雷识别方法及系统
CN117117965A (zh) 基于磁敏电流传感器的配电网拓扑识别系统、方法及装置
Sumner et al. Improved power quality control and intelligent protection for grid connected power electronic converters, using real time parameter estimation
CN102983572A (zh) 谐波和无功电流提取方法
CN110007146A (zh) 一种基于电压电流谐波相位的谐振点检测方法
CN209055625U (zh) 一种电力电缆局部放电定位检测装置
CN107167663A (zh) 一种用于抑制工频感应电压的组合测量电路及测量方法

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