CN110007141B - 一种基于电压电流谐波相似度的谐振点检测方法 - Google Patents

一种基于电压电流谐波相似度的谐振点检测方法 Download PDF

Info

Publication number
CN110007141B
CN110007141B CN201910341957.5A CN201910341957A CN110007141B CN 110007141 B CN110007141 B CN 110007141B CN 201910341957 A CN201910341957 A CN 201910341957A CN 110007141 B CN110007141 B CN 110007141B
Authority
CN
China
Prior art keywords
voltage
current
harmonic
similarity
resonance point
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.)
Active
Application number
CN201910341957.5A
Other languages
English (en)
Other versions
CN110007141A (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.)
Qingdao Tuowei Technology Co.,Ltd.
Qingdao Zhidian New Energy Technology Co ltd
Qingdao Topscomm Communication Co Ltd
Original Assignee
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 Qingdao Topscomm Communication Co Ltd filed Critical Qingdao Topscomm Communication Co Ltd
Priority to CN201910341957.5A priority Critical patent/CN110007141B/zh
Publication of CN110007141A publication Critical patent/CN110007141A/zh
Application granted granted Critical
Publication of CN110007141B publication Critical patent/CN110007141B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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算法的计算公式为，
Figure BDA0002041005510000021
Figure BDA0002041005510000022
其中，ak代表第k次谐波的实部，bk代表第k次谐波的虚部，N代表工频周期的数据点数，k代表谐波次数，T＝20ms。
进一步地，步骤四中的余弦相似度计算公式为，
Figure BDA0002041005510000023
其中，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算法的计算公式为，
Figure BDA0002041005510000031
Figure BDA0002041005510000032
其中，ak代表第k次谐波的实部，bk代表第k次谐波的虚部，N代表工频周期的采样点数，N＝6250，k代表谐波次数，k＝11,22,25,59，T＝20ms。
步骤四，根据余弦相似度计算公式，分别计算电压电流11、22、25、59次谐波波形之间的相似度，余弦相似度的计算公式为，
Figure BDA0002041005510000041
其中，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 (3)

1.一种基于电压电流谐波相似度的谐振点检测方法，其特征在于：包括以下步骤，
步骤一，利用低频载波发送机向变压器低压侧注入幅值相等、频率不同的载波电流信号，频率范围在1K～5KHz之间，幅值范围在5～15V之间；
步骤二，检测变压器低压侧的电压信号U(t)与电流信号I(t)；
步骤三，根据滑动DFT算法，提取出测量信号中的电压电流谐波；
步骤四，根据余弦相似度计算公式，计算电压电流谐波波形之间的相似度；
步骤五，根据计算得到的电压电流谐波之间的相似度，来确定谐振点，该方法检测出的谐振点位于一定的频率范围；
步骤三中根据滑动DFT算法提取出测量信号中的电压电流谐波，其中滑动DFT算法的计算公式为，
Figure FDA0003245265850000011
Figure FDA0003245265850000012
其中，ak代表第k次谐波的实部，bk代表第k次谐波的虚部，N代表工频周波的数据点数，k代表谐波次数，T＝20ms。
2.根据权利要求1所述的一种基于电压电流谐波相似度的谐振点检测方法，其特征在于：步骤四中根据余弦相似度计算公式计算电压电流谐波波形之间的相似度，其中余弦相似度计算公式为，
Figure FDA0003245265850000013
其中，Ik为提取的第k次电流谐波，Uk为提取的第k次电压谐波，N为用于计算相似度的数据点数，rk为计算得到的第k次电压电流谐波波形之间的相似度。
3.根据权利要求1所述的一种基于电压电流谐波相似度的谐振点检测方法，其特征在于：步骤五中根据计算得到的电压电流谐波之间的相似度来确定谐振点的原则为，rk越大，第k次谐波越靠近谐振点，当rk＝1时，电压第k次谐波与电流第k次谐波之间的相似度达到最大，第k次谐波对应的频率为谐振点。
CN201910341957.5A 2019-04-26 2019-04-26 一种基于电压电流谐波相似度的谐振点检测方法 Active CN110007141B (zh)

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 CN110007141A (zh) 2019-07-12
CN110007141B true 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)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110244123B (zh) * 2019-06-24 2021-07-02 国网江西省电力有限公司鹰潭供电分公司 一种基于线变关系识别装置的降阶广义积分器谐振点检测方法
CN111913040A (zh) * 2020-05-23 2020-11-10 青岛鼎信通讯股份有限公司 一种基于负载投切和相似度的谐振点的线变关系识别方法
CN112213578A (zh) * 2020-09-23 2021-01-12 青岛鼎信通讯股份有限公司 一种中压线变关系识别方法
CN112379178B (zh) * 2020-10-28 2022-11-22 国网安徽省电力有限公司合肥供电公司 存在延时的两种波形相似度判断方法、系统及存储介质
CN112557822A (zh) * 2020-12-05 2021-03-26 青岛鼎信通讯股份有限公司 一种中压台区拓扑关系定位识别方法
CN115589237B (zh) * 2022-12-13 2023-03-17 青岛鼎信通讯股份有限公司 一种适用于电力领域的高频电流信号分支归属判决方法
CN119154206A (zh) * 2024-11-11 2024-12-17 湖南大学 基于多频段特征的逆变型微电网线路故障方向检测方法

Citations (7)

* Cited by examiner, † Cited by third party
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谐振及相位检测的电感值测量方法及电路
CN104772450A (zh) * 2015-01-15 2015-07-15 中南大学 一种铸造用超声电源系统及谐振工作点自动跟踪方法
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 南方电网科学研究院有限责任公司 电网谐波普测分析方法、装置、设备及介质

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104330627B (zh) * 2014-10-24 2017-04-26 华中科技大学 一种基于有源注入电流的电网谐振自动检测方法
CN106771584A (zh) * 2016-11-16 2017-05-31 合肥普望电子有限责任公司 一种应用于配电网的谐波检测方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
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谐振及相位检测的电感值测量方法及电路
CN104772450A (zh) * 2015-01-15 2015-07-15 中南大学 一种铸造用超声电源系统及谐振工作点自动跟踪方法
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 (2)

* Cited by examiner, † Cited by third party
Title
《一种简化DFT的滑窗迭代算法在电力谐波检测中的应用》;成剑等;《计算机测量与控制》;20050831;第771页-第774页 *
《基于三次谐波频率跟踪方法的超声波驱动电路设计》;魏臣;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20150515;第4页、第9页 *

Also Published As

Publication number Publication date
CN110007141A (zh) 2019-07-12

Similar Documents

Publication Publication Date Title
CN110007141B (zh) 一种基于电压电流谐波相似度的谐振点检测方法
CN109375060B (zh) 一种配电网故障波形相似度计算方法
CN101738551B (zh) 基于联网的暂态电能质量扰动智能分析方法
CN104569744B (zh) 一种适用于配电网线路的综合单端故障定位方法
CN104090214B (zh) 一种电缆故障检测及老化分析方法
CN105004939B (zh) 一种复合电能质量扰动信号定量分析方法
CN110687399B (zh) 一种配电网故障指示器波形故障开始时刻判断方法
Iksan et al. Appliances identification method of non-intrusive load monitoring based on load signature of VI trajectory
CN111308260B (zh) 一种基于小波神经网络的电能质量监测和电器故障分析系统及其工作方法
CN108535613A (zh) 一种基于组合窗函数的电压闪变参数检测方法
CN107942192A (zh) 一种配电设备ct自动测试方法及系统
CN104833935A (zh) 一种针对有功功率方向法定位谐波源局限性的检验方法
CN110007146B (zh) 一种基于电压电流谐波相位的谐振点检测方法
CN111157843A (zh) 一种基于时频域行波信息的配电网选线方法
CN104777356A (zh) 一种基于神经网络的实时高精度谐波检测方法
CN103245830B (zh) 一种结合ar谱估计与非线性优化的间谐波检测方法
CN110133371A (zh) 一种基于恒压源变换器的电网谐振点检测方法
CN111913040A (zh) 一种基于负载投切和相似度的谐振点的线变关系识别方法
CN110007145B (zh) 一种基于恒压源变换器和电压电流相位差的谐振点检测方法
CN110007143B (zh) 一种基于电压谐波幅值的谐振点检测方法
CN106596025A (zh) 基于冲激响应的公路隧道悬挂风机基础稳定性检测方法及系统
CN104730384A (zh) 基于不完全s变换的电能扰动识别与定位方法
Zaman et al. Multimode synchronous resonance detection in converter dominated power system using synchro-waveforms
CN104808643B (zh) 一种基于改进的倒双谱分析的控制回路非线性检测方法
CN106841931A (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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20241122

Address after: No. 858, Huaguan Road, high tech Zone, Qingdao, Shandong Province 266000

Patentee after: QINGDAO TOPSCOMM COMMUNICATION Co.,Ltd.

Country or region after: China

Patentee after: QINGDAO DINGJUN ELECTRIC Co.,Ltd.

Address before: 266024 area B, 5th floor, building 6, 288 Ningxia road, Qingdao City, Shandong Province

Patentee before: QINGDAO TOPSCOMM COMMUNICATION Co.,Ltd.

Country or region before: China

CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: No. 858, Huaguan Road, high tech Zone, Qingdao, Shandong Province 266000

Patentee after: QINGDAO TOPSCOMM COMMUNICATION Co.,Ltd.

Country or region after: China

Patentee after: Qingdao Zhidian New Energy Technology Co.,Ltd.

Patentee after: QINGDAO DINGJUN ELECTRIC Co.,Ltd.

Address before: No. 858, Huaguan Road, high tech Zone, Qingdao, Shandong Province 266000

Patentee before: QINGDAO TOPSCOMM COMMUNICATION Co.,Ltd.

Country or region before: China

Patentee before: QINGDAO DINGJUN ELECTRIC Co.,Ltd.

CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: No. 858, Huaguan Road, high tech Zone, Qingdao, Shandong Province 266000

Patentee after: QINGDAO TOPSCOMM COMMUNICATION Co.,Ltd.

Country or region after: China

Patentee after: Qingdao Zhidian New Energy Technology Co.,Ltd.

Patentee after: Qingdao Tuowei Technology Co.,Ltd.

Address before: No. 858, Huaguan Road, high tech Zone, Qingdao, Shandong Province 266000

Patentee before: QINGDAO TOPSCOMM COMMUNICATION Co.,Ltd.

Country or region before: China

Patentee before: Qingdao Zhidian New Energy Technology Co.,Ltd.

Patentee before: QINGDAO DINGJUN ELECTRIC Co.,Ltd.