CN110007141A - 一种基于电压电流谐波相似度的谐振点检测方法 - Google Patents
一种基于电压电流谐波相似度的谐振点检测方法 Download PDFInfo
- Publication number
- CN110007141A CN110007141A CN201910341957.5A CN201910341957A CN110007141A CN 110007141 A CN110007141 A CN 110007141A CN 201910341957 A CN201910341957 A CN 201910341957A CN 110007141 A CN110007141 A CN 110007141A
- Authority
- CN
- China
- Prior art keywords
- voltage
- similarity
- harmonic wave
- resonance point
- current 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
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/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
- G01R23/06—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage by converting frequency into an amplitude of current or voltage
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Current Or Voltage (AREA)
Abstract
本发明公开了一种基于电压电流谐波相似度的谐振点检测方法,包括以下步骤,步骤一,利用低频载波发送机向变压器低压侧注入幅值相等、频率不同的载波电流信号;步骤二,检测变压器低压侧的电压信号U(t)与电流信号I(t);步骤三,根据滑动DFT算法,提取出测量信号中的电压电流谐波;步骤四,根据余弦相似度计算公式,计算电压电流谐波波形之间的相似度;步骤五,根据计算得到的电压电流谐波之间的相似度,来确定谐振点。本发明此方法实现简单,只需测量变压器低压侧的电压电流信号,另外,该方法计算量小,测量数据无需上传到主站,可本地进行分析,可以实现对谐振点的实时在线检测。
Description
技术领域
本发明涉及配电网自动化***领域,涉及一种基于电压电流谐波相似度的谐振点检测方法。
背景技术
功率因数是电力***中一个很重要的性能指标,为了提高***的功率因数,普遍的做法是向***中投切无功补偿电容器组,并且一般根据***不同的负载大小,来投入不同容量的电容器组,提高电力***输电线路的电压,降低线路损耗,提高经济效益。但是由于电网***中普遍存在谐波源,当向***中投切补偿电容时,很可能导致***与补偿电容发生之间并联谐振,导致补偿电容产生很大的谐振电流,从而烧毁补偿电容器,甚至损坏其他的一些电力设备、威胁操作人员的生命安全,造成严重的电力事故。因此,提前准确检测出***的谐振点,对于抑制电网的谐振具有重要的意义。
目前已有公司和高校在谐振点检测方面取得了很多成果,并且已经制造出专门的检测设备,但是此类设备价格昂贵,检测过程操作复杂,速度较慢。当设备检测到所需的信号后,还要上传到主站来进行分析,检测过程具有滞后性,不能实时在线检测。因此,需要发展一种实现简单、能够实时在线检测谐振点的方法。
发明内容
本发明针对上述问题,克服现有技术的不足,提出一种基于电压电流谐波相似度的谐振点检测方法,该方法根据滑动DFT算法提取的电压电流谐波之间的相似性来确定谐振点。该方法实现简单,只需测量出变压器低压侧的电压电流信号,就能够实现谐振点的检测。另外,该方法计算量小,测量信号无需上传到主站,在本地设备即可进行分析,可实现实时在线检测,具有很好的工程实用性。
本发明通过计算电压电流谐波之间的相似度来检测谐振点。该方法首先利用低频载波发送机向变压器低压侧注入幅值相等、频率不同的载波电流信号,同时检测变压器低压侧的电压电流信号,然后根据滑动DFT算法提取出电压电流信号中的谐波波形,最后根据电压电流谐波之间的相似度来检测谐振点。
为了实现上述发明目的,本发明采取如下技术方案:
一种基于电压电流谐波相似度的谐振点检测方法,包括以下步骤,
步骤一,利用低频载波发送机向变压器低压侧注入幅值相等、频率不同的载波电流信号。
步骤二,检测变压器低压侧的电压信号U(t)与电流信号I(t)。
步骤三,根据滑动DFT算法,提取出测量信号中的电压电流谐波。
步骤四,根据余弦相似度计算公式,计算电压电流谐波波形之间的相似度。
步骤五,根据计算得到的电压电流谐波之间的相似度,来确定谐振点。
进一步地,步骤三中的滑动DFT算法的计算公式为,
其中,ak代表第k次谐波的实部,bk代表第k次谐波的虚部,N代表工频周期的数据点数,k代表谐波次数,T=20ms。
进一步地,步骤四中的余弦相似度计算公式为,
其中,Ik为提取的第k次电流谐波,Uk为提取的第k次电压谐波,N为用于计算相似度的数据点数,rk为计算得到的第k次电压电流谐波波形之间的相似度。
进一步地,步骤五中的根据计算得到的电压电流谐波之间的相似度来确定谐振点的原则为,rk越大,第k次谐波越靠近谐振点,当rk=1时,电压第k次谐波与电流第k次谐波之间的相似度达到最大,第k次谐波对应的频率为谐振点。
进一步地,步骤一中频率范围在1K~5KHz之间。
进一步地,步骤一中幅值范围在5~15V之间。
本发明的有益效果是:本发明一种基于电压电流谐波相似度的谐振点检测方法,根据计算得到的电压电流谐波之间的相似度来检测***的谐振点,此方法实现简单,只需测量变压器低压侧的电压电流信号,另外,该方法计算量小,测量数据无需上传到主站,可本地进行分析,可以实现对谐振点的实时在线检测。
附图说明
图1为本发明基于电压电流谐波相似度的谐振点检测方法总体流程图。
图2为本发明注入11次低频载波电流信号时测量的变压器低压侧电压电流信号波形。
图3为本发明注入22次低频载波电流信号时测量的变压器低压侧电压电流信号波形。
图4为本发明注入25次低频载波电流信号时测量的变压器低压侧电压电流信号波形。
图5为本发明注入59次低频载波电流信号时测量的变压器低压侧电压电流信号波形。
图6为本发明提取的电压电流测量信号中的11次谐波波形。
图7为本发明提取的电压电流测量信号中的22次谐波波形。
图8为本发明提取的电压电流测量信号中的25次谐波波形。
图9为本发明提取的电压电流测量信号中的59次谐波波形。
图10为本发明提取的注入11次低频载波电流信号时补偿电容产生的11次谐波电流波形。
图11为本发明提取的注入22次低频载波电流信号时补偿电容产生的22次谐波电流波形。
图12为本发明提取的注入25次低频载波电流信号时补偿电容产生的25次谐波电流波形。
图13为本发明提取的注入59次低频载波电流信号时补偿电容产生的59次谐波电流波形。
其中横坐标代表采样点,纵坐标代表电压电流的幅值,单位分别为V、A。
具体实施方式
下面结合附图1-13和实施例对本发明作进一步的说明,以具体阐述本发明的技术方案。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
结合附图1,本发明的一种基于电压电流谐波相似度的谐振点检测方法,包括以下步骤,
步骤一,利用低频载波发送机向变压器低压侧注入幅值为7A、频率分别为的11、22、25、59次的载波电流信号。
步骤二,检测变压器低压侧的电压信号U(t)与电流信号I(t),其波形如图2~图5所示。
步骤三,根据滑动DFT算法,分别提取注入11、22、25、59次载波电流信号时变压器低压侧电压电流信号中的11、22、25、59次谐波,其波形分别如图6~图9所示,其中滑动DFT算法的计算公式为,
其中,ak代表第k次谐波的实部,bk代表第k次谐波的虚部,N代表工频周期的采样点数,N=6250,k代表谐波次数,k=11,22,25,59,T=20ms。
步骤四,根据余弦相似度计算公式,分别计算电压电流11、22、25、59次谐波波形之间的相似度,余弦相似度的计算公式为,
其中,Ik为提取的第k次电流谐波,Uk为提取的第k次电压谐波,N为用于计算相似度的数据点数,rk为计算得到的第k次电压电流谐波波形之间的相似度,k=11,22,25,59。
步骤五,根据计算得到的电压电流谐波之间的相似度,来确定谐振点,其原则为rk越大,第k次谐波越靠近谐振点,当rk=1时,电压第k次谐波与电流第k次谐波之间的相似度达到最大,第k次谐波对应的频率为谐振点。经提取计算得到,第11次电压电流谐波之间的相似度为0.5209,第22次电压电流谐波之间的相似度为0.8917,第25次电压电流谐波之间的相似度为0.9421,第59次电压电流谐波之间的相似度为0.0729,经比较可见,第25次电压电流谐波之间的相似度最靠近1,可以判定第25谐波最靠近谐振点,因此,谐振点在第25次谐波频率附近。
本实施例中:利用实际的现场数据对本发明中的谐振点检测方法进行了测试验证。本发明根据电压电流谐波之间的相似度确定了谐振点在第25次谐波频率附近,同时,由图10~图13可知,注入25次载波电流信号时,补偿电容产生的谐波电流要比注入11、22、59次载波电流信号时补偿电容产生的电流大很多,说明在25次谐波频率附近发生了谐振,由此也进一步验证了本发明所提出方法的准确性。在实际的检测过程中,由于利用载波发送机向电网中注入载波电流信号,当检测到快接近谐振点的时候,会跳过谐振点所在的一段频率范围,避免引起***的谐振。因此,此方法最终的检测结果是检测出谐振点位于一定的频率范围。
本发明步骤一中频率范围在1K~5KHz之间,可以是1KHz、5KHz或1K~5KHz之间的任意数值,并不限于实施例所给出的数值。步骤一中幅值范围在5~15V之间,可以是5V、15V或5~15V之间的任意数值,并不限于实施例所给出的数值。
综上所述,本发明提出一种基于电压电流谐波相似度的谐振点检测方法,首先利用低频载波发送机向变压器低压侧注入幅值相等、频率不同的载波电流信号,同时检测变压器低压侧的电压电流信号,然后根据滑动DFT算法提取出电压电流信号中的电压电流谐波波形,最后根据电压电流谐波之间的相似度来检测谐振点。此方法只需要测量变压器低压侧的电压电流信号,实现简单。另外,此方法计算量小,测量数据不用上传到主站,在设备本地即可进行分析,能够对谐振点实现实时在线检测,该方法具有很好的工程实用性。
以上实施例是对本发明的具体实施方式的说明,而非对本发明的限制,有关技术领域的技术人员在不脱离本发明的精神和范围的情况下,还可以做出各种变换和变化而得到相对应的等同的技术方案,因此所有等同的技术方案均应该归入本发明的专利保护范围。
Claims (6)
1.一种基于电压电流谐波相似度的谐振点检测方法,其特征在于:包括以下步骤,
步骤一,利用低频载波发送机向变压器低压侧注入幅值相等、频率不同的载波电流信号;
步骤二,检测变压器低压侧的电压信号U(t)与电流信号I(t);
步骤三,根据滑动DFT算法,提取出测量信号中的电压电流谐波;
步骤四,根据余弦相似度计算公式,计算电压电流谐波波形之间的相似度;
步骤五,根据计算得到的电压电流谐波之间的相似度,来确定谐振点。
2.根据权利要求1所述的一种基于电压电流谐波相似度的谐振点检测方法,其特征在于:
步骤三中根据滑动DFT算法提取出测量信号中的电压电流谐波,其中滑动DFT算法的计算公式为,
其中,ak代表第k次谐波的实部,bk代表第k次谐波的虚部,N代表工频周波的数据点数,k代表谐波次数,T=20ms。
3.根据权利要求1或2所述的一种基于电压电流谐波相似度的谐振点检测方法,其特征在于:步骤四中根据余弦相似度计算公式计算电压电流谐波波形之间的相似度,其中余弦相似度计算公式为,
其中,Ik为提取的第k次电流谐波,Uk为提取的第k次电压谐波,N为用于计算相似度的数据点数,rk为计算得到的第k次电压电流谐波波形之间的相似度。
4.根据权利要求1或2所述的一种基于电压电流谐波相似度的谐振点检测方法,其特征在于:步骤五中根据计算得到的电压电流谐波之间的相似度来确定谐振点的原则为,rk越大,第k次谐波越靠近谐振点,当rk=1时,电压第k次谐波与电流第k次谐波之间的相似度达到最大,第k次谐波对应的频率为谐振点。
5.根据权利要求1或2所述的一种基于电压电流谐波相似度的谐振点检测方法,其特征在于:步骤一中频率范围在1K~5KHz之间。
6.根据权利要求1或2所述的一种基于电压电流谐波相似度的谐振点检测方法,其特征在于:步骤一中幅值范围在5~15V之间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910341957.5A CN110007141B (zh) | 2019-04-26 | 2019-04-26 | 一种基于电压电流谐波相似度的谐振点检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910341957.5A CN110007141B (zh) | 2019-04-26 | 2019-04-26 | 一种基于电压电流谐波相似度的谐振点检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110007141A true CN110007141A (zh) | 2019-07-12 |
CN110007141B CN110007141B (zh) | 2021-11-05 |
Family
ID=67174424
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910341957.5A Active CN110007141B (zh) | 2019-04-26 | 2019-04-26 | 一种基于电压电流谐波相似度的谐振点检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110007141B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110244123A (zh) * | 2019-06-24 | 2019-09-17 | 国网江西省电力有限公司鹰潭供电分公司 | 一种基于线变关系识别装置的降阶广义积分器谐振点检测方法 |
CN111913040A (zh) * | 2020-05-23 | 2020-11-10 | 青岛鼎信通讯股份有限公司 | 一种基于负载投切和相似度的谐振点的线变关系识别方法 |
CN112213578A (zh) * | 2020-09-23 | 2021-01-12 | 青岛鼎信通讯股份有限公司 | 一种中压线变关系识别方法 |
CN112379178A (zh) * | 2020-10-28 | 2021-02-19 | 国网安徽省电力有限公司合肥供电公司 | 存在延时的两种波形相似度判断方法、***及存储介质 |
CN112557822A (zh) * | 2020-12-05 | 2021-03-26 | 青岛鼎信通讯股份有限公司 | 一种中压台区拓扑关系定位识别方法 |
CN115589237A (zh) * | 2022-12-13 | 2023-01-10 | 青岛鼎信通讯股份有限公司 | 一种适用于电力领域的高频电流信号分支归属判决方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103066948A (zh) * | 2012-11-29 | 2013-04-24 | 中国科学院高能物理研究所 | 谐振励磁电源的电流给定波形的控制方法 |
CN104330627A (zh) * | 2014-10-24 | 2015-02-04 | 华中科技大学 | 一种基于有源注入电流的电网谐振自动检测方法 |
CN104655933A (zh) * | 2013-11-21 | 2015-05-27 | 西安丁子电子信息科技有限公司 | 一种基于lc谐振及相位检测的电感值测量方法及电路 |
CN104772450A (zh) * | 2015-01-15 | 2015-07-15 | 中南大学 | 一种铸造用超声电源***及谐振工作点自动跟踪方法 |
CN106645974A (zh) * | 2016-12-30 | 2017-05-10 | 北京金风科创风电设备有限公司 | 用于确定风机与电网谐振的方法及*** |
KR101735776B1 (ko) * | 2017-01-31 | 2017-05-15 | 주식회사 일산이엔지 | 신호의 혼합 및 공진회로를 포함하는 비접촉형 픽업 코일을 이용한 전력 전송 버스바 상태 모니터링을 위한 특정 고조파 주파수 선별 방법 및 장치 |
CN106771584A (zh) * | 2016-11-16 | 2017-05-31 | 合肥普望电子有限责任公司 | 一种应用于配电网的谐波检测方法 |
CN107769205A (zh) * | 2017-10-30 | 2018-03-06 | 许继电源有限公司 | 一种有源电力滤波器谐振抑制方法及装置 |
CN109100568A (zh) * | 2018-06-25 | 2018-12-28 | 南方电网科学研究院有限责任公司 | 电网谐波普测分析方法、装置、设备及介质 |
-
2019
- 2019-04-26 CN CN201910341957.5A patent/CN110007141B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103066948A (zh) * | 2012-11-29 | 2013-04-24 | 中国科学院高能物理研究所 | 谐振励磁电源的电流给定波形的控制方法 |
CN104655933A (zh) * | 2013-11-21 | 2015-05-27 | 西安丁子电子信息科技有限公司 | 一种基于lc谐振及相位检测的电感值测量方法及电路 |
CN104330627A (zh) * | 2014-10-24 | 2015-02-04 | 华中科技大学 | 一种基于有源注入电流的电网谐振自动检测方法 |
CN104772450A (zh) * | 2015-01-15 | 2015-07-15 | 中南大学 | 一种铸造用超声电源***及谐振工作点自动跟踪方法 |
CN106771584A (zh) * | 2016-11-16 | 2017-05-31 | 合肥普望电子有限责任公司 | 一种应用于配电网的谐波检测方法 |
CN106645974A (zh) * | 2016-12-30 | 2017-05-10 | 北京金风科创风电设备有限公司 | 用于确定风机与电网谐振的方法及*** |
KR101735776B1 (ko) * | 2017-01-31 | 2017-05-15 | 주식회사 일산이엔지 | 신호의 혼합 및 공진회로를 포함하는 비접촉형 픽업 코일을 이용한 전력 전송 버스바 상태 모니터링을 위한 특정 고조파 주파수 선별 방법 및 장치 |
CN107769205A (zh) * | 2017-10-30 | 2018-03-06 | 许继电源有限公司 | 一种有源电力滤波器谐振抑制方法及装置 |
CN109100568A (zh) * | 2018-06-25 | 2018-12-28 | 南方电网科学研究院有限责任公司 | 电网谐波普测分析方法、装置、设备及介质 |
Non-Patent Citations (3)
Title |
---|
AHMET S. YILMAZ: "《Applications of parametric spectral estimation methods on detection of power system harmonics》", 《ELECTRIC POWER SYSTEMS RESEARCH》 * |
成剑等: "《一种简化DFT的滑窗迭代算法在电力谐波检测中的应用》", 《计算机测量与控制》 * |
魏臣: "《基于三次谐波频率跟踪方法的超声波驱动电路设计》", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110244123A (zh) * | 2019-06-24 | 2019-09-17 | 国网江西省电力有限公司鹰潭供电分公司 | 一种基于线变关系识别装置的降阶广义积分器谐振点检测方法 |
CN110244123B (zh) * | 2019-06-24 | 2021-07-02 | 国网江西省电力有限公司鹰潭供电分公司 | 一种基于线变关系识别装置的降阶广义积分器谐振点检测方法 |
CN111913040A (zh) * | 2020-05-23 | 2020-11-10 | 青岛鼎信通讯股份有限公司 | 一种基于负载投切和相似度的谐振点的线变关系识别方法 |
CN112213578A (zh) * | 2020-09-23 | 2021-01-12 | 青岛鼎信通讯股份有限公司 | 一种中压线变关系识别方法 |
CN112379178A (zh) * | 2020-10-28 | 2021-02-19 | 国网安徽省电力有限公司合肥供电公司 | 存在延时的两种波形相似度判断方法、***及存储介质 |
CN112557822A (zh) * | 2020-12-05 | 2021-03-26 | 青岛鼎信通讯股份有限公司 | 一种中压台区拓扑关系定位识别方法 |
CN115589237A (zh) * | 2022-12-13 | 2023-01-10 | 青岛鼎信通讯股份有限公司 | 一种适用于电力领域的高频电流信号分支归属判决方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110007141B (zh) | 2021-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110007141A (zh) | 一种基于电压电流谐波相似度的谐振点检测方法 | |
CN101459334B (zh) | 电力***故障信息获取方法 | |
CN102253283B (zh) | 一种基于小波包能量谱的分布式微网并网孤岛检测方法 | |
CN100592441C (zh) | 特高压输电线路用电子式电流互感器及其电晕损失测量装置 | |
CN105911499B (zh) | 现场环境下超声波局部放电计量***及方法 | |
CN107942192A (zh) | 一种配电设备ct自动测试方法及*** | |
CN203025281U (zh) | 变压器异常振动实验装置 | |
CN105092997B (zh) | 一种特高压输电线路雷电绕击与反击的识别方法 | |
CN208818818U (zh) | 一种低功耗程控式电流互感器现场检定装置 | |
CN108548962A (zh) | 一种接地网交流、直流和冲击电阻综合测试***及方法 | |
CN201607479U (zh) | 一种电网过电压无接触在线监测与识别一体化装置 | |
CN104375007A (zh) | 一种基于相位识别的土壤电阻率测量方法及装置 | |
CN102628890A (zh) | 智能型接地安全参数监测方法及装置 | |
CN104931793B (zh) | 一种变电站接地网接地阻抗获取方法 | |
CN108362928A (zh) | 接触网电压非接触式测量设备及方法 | |
CN104655967B (zh) | 配电变压器绕组振动信号特征量提取方法 | |
CN101819224A (zh) | 基于微分环的雷电流测量传感器 | |
CN202975173U (zh) | 一种变电站配电***中的对地电容测量仪 | |
CN102749595B (zh) | 一种变压器铁芯剩余磁感应强度的测量方法 | |
CN110133371A (zh) | 一种基于恒压源变换器的电网谐振点检测方法 | |
CN105182116B (zh) | 一种基于加权梯度结构相似度的变压器绕组工作状态检测方法 | |
CN110007146A (zh) | 一种基于电压电流谐波相位的谐振点检测方法 | |
CN105866592A (zh) | 动态无功补偿响应波形采集***及采集方法 | |
CN203204090U (zh) | 变电站氧化锌避雷器在线监测装置 | |
CN104111438B (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 |