CN115372698A - 一种用于电力系统高次谐波抑制的测控装置及方法 - Google Patents

一种用于电力系统高次谐波抑制的测控装置及方法 Download PDF

Info

Publication number
CN115372698A
CN115372698A CN202211314350.6A CN202211314350A CN115372698A CN 115372698 A CN115372698 A CN 115372698A CN 202211314350 A CN202211314350 A CN 202211314350A CN 115372698 A CN115372698 A CN 115372698A
Authority
CN
China
Prior art keywords
harmonic
amplitude
measurement
sampling
calculation
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
Application number
CN202211314350.6A
Other languages
English (en)
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.)
Dongfang Electronics Co Ltd
Original Assignee
Dongfang Electronics 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 Dongfang Electronics Co Ltd filed Critical Dongfang Electronics Co Ltd
Priority to CN202211314350.6A priority Critical patent/CN115372698A/zh
Publication of CN115372698A publication Critical patent/CN115372698A/zh
Pending legal-status Critical Current

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/16—Spectrum analysis; Fourier analysis
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/141—Discrete Fourier transforms
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/01—Arrangements for reducing harmonics or ripples
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/40—Arrangements for reducing harmonics

Abstract

本申请公开了一种用于电力系统高次谐波抑制的测控装置及方法，属于自动化测控的技术领域，包括：根据设定采样频率对模拟信号进行采样，得到采样信号；根据信号采样模块输出的采样数据，进行重采样存入计算缓冲区，为谐波检测模块提供检测数据，并为测量值计算模块提供计算数据；再根据采样数据转换模块输出的重采样数据作为输入，检测各次谐波分量大小，从而确定是否存在各次高次谐波；最后依据谐波检测模块的检测结果，按照测量幅值计算方案，并基于傅里叶算法，精确计算测量值。本申请具有能够对高次谐波进行有效的检测和抑制，极大程度地提高谐波频率分量的测量精度，进而维持电力系统的稳定运行的效果。

Description

一种用于电力系统高次谐波抑制的测控装置及方法
技术领域
本申请涉及自动化测控的技术领域，尤其涉及一种用于电力系统高次谐波抑制的测控装置及方法。
背景技术
目前，随着非线性元器件的大量使用，电力系统面临着严重的谐波污染。在电力系统中，发电设备（电源端）、输配电设备和负荷的非线性等容易造成波形畸变，产生不同分量的谐波和间谐波。谐波会引起变压器、旋转设备和线路的附加损耗和发热，而且因长期振动导致金属疲劳和机械损坏，损坏系统设备，还会引起系统谐振，危机电力系统的安全运行。另一方面，电网的运行监控和调度依赖于变电站测控装置对电参量的精确测量，而高次谐波的注入直接影响工频状态下测控装置的稳态测量精度，引起较大的测量误差。
现有公告号为CN108377154B的中国专利公开了一种谐波抑制装置、谐波抑制方法及电子装置，谐波抑制装置包括：第一获取单元，用于获取电子装置所使用的运营商的信息，以得到电子装置中进行通信时使用的电磁波信号的第一频段；第一反相单元，用于根据第一频段，产生第一反相谐波信号；第一输出单元，用于将第一反相谐波信号输出，第一反相谐波信号用于消除电子装置使用当前运营商时产生的至少部分辐射杂散。
针对上述相关技术，采用有源滤波电路进行谐波抑制操作，由于上述操作过程无法对谐波分量进行精确区分，对谐波频率分量的量测精度偏低，因此，难以对谐波频率分量进行较为精确的测量操作。
发明内容
本申请针对现有技术存在的不足，提供一种用于电力系统高次谐波抑制的测控装置及方法，能够对高次谐波进行有效的检测，极大程度地提高谐波频率分量的测量精度，进而维持电力系统的稳定运行。
本申请解决上述技术问题的技术方案如下：
一种用于电力系统高次谐波抑制的测控方法，包括：
根据离散傅里叶变换算法，得出频率系数
Figure 795848DEST_PATH_IMAGE001
，公式如下：
Figure 598719DEST_PATH_IMAGE002
（7）
公式（7）中，u(n)表示已知的电压输入信号，N为周期采样点数，其中，0≤k≤（N-1）；
利用欧拉公式展开所述频率系数
Figure 268735DEST_PATH_IMAGE001
，得到所述频率系数
Figure 524267DEST_PATH_IMAGE001
的实部
Figure 180507DEST_PATH_IMAGE003
和虚部
Figure 696939DEST_PATH_IMAGE004
，公式如（8）和（9）所示：
Figure 678802DEST_PATH_IMAGE005
（8）
Figure 280684DEST_PATH_IMAGE006
（9）
其中，N为周期采样点数，k为谐波次数，k=1，2……13；
计算得出k个频率分量的幅值
Figure 737686DEST_PATH_IMAGE007
，即电压的基波分量和各次谐波分量的幅值，公 式如下：
Figure 249570DEST_PATH_IMAGE008
（10）
其中，k为谐波次数，k=1，2……13；
根据所述幅值
Figure 261388DEST_PATH_IMAGE007
，计算得到基波和2-13次谐波的幅值大小；
当2-13次中的某次谐波幅值超过设定阈值后，认定存在该次谐波分量，根据所述 幅值
Figure 960354DEST_PATH_IMAGE007
，计算得到电压的有效值
Figure 348610DEST_PATH_IMAGE009
，公式如下：
Figure 449421DEST_PATH_IMAGE010
（11）
根据所述电压的有效值
Figure 897720DEST_PATH_IMAGE009
，计算得到电压和电流的有效值；
对于13次以上的谐波分量，不进行检测，不加入有效值计算，实现高次谐波抑制作用；
同理，根据离散傅里叶变换算法以及欧拉公式，计算得到电流频率系数的实部
Figure 349561DEST_PATH_IMAGE011
和虚部
Figure 275929DEST_PATH_IMAGE012
，得出电流的基波和各次谐波幅值
Figure 496826DEST_PATH_IMAGE013
，其中，计算电流有效值
Figure 991392DEST_PATH_IMAGE014
，公式如 下：
Figure 789584DEST_PATH_IMAGE015
（12）
根据有功功率P和无功功率Q的计算公式：
Figure 395009DEST_PATH_IMAGE016
（13）
Figure 595046DEST_PATH_IMAGE017
（14）
分别得出电力信号有功功率P和无功功率Q的推导公式：
Figure 251724DEST_PATH_IMAGE018
（15）
Figure 412578DEST_PATH_IMAGE019
（16）。
通过采用上述技术方案，通过离散傅里叶变换算法可将电压和电流信号按照基波频率的整数倍进行分解，分解出基波和各次谐波分量。能够对谐波频率分量进行较为精确的测量操作，能够对高次谐波进行有效的检测，极大程度地提高谐波频率分量的测量精度，进而有助于维持电力系统的稳定运行。
进一步的，当k取值为2-13时，基波幅值和谐波幅值采用40点采样数据进行计算，电压、电流有效值等于基波幅值和谐波幅值之和。
进一步的，当k取值为14-38或40-50时，谐波幅值不参与有效值计算，基波幅值采用40点采样数据进行计算，电压、电流有效值等于基波幅值。
进一步的，当k取值大于等于39时，谐波幅值不参与有效值计算，基波幅值采用80点采样数据进行计算，电压、电流有效值等于基波幅值。
进一步的，谐波分量全部采用80点采样数据进行检测，电压和电流幅值同样采用80点采样数据进行计算，根据谐波检测算法：
Figure 946328DEST_PATH_IMAGE020
（8）
Figure 876238DEST_PATH_IMAGE021
（9）
其中，k为谐波次数，k=1，2……13；
分解得到电压和电流信号的基波和2-13次谐波分量，再根据
Figure 571661DEST_PATH_IMAGE022
（10）
计算得到基波和2-13次谐波的幅值大小。
进一步的，计算缓冲区采样数据填充完成之前，谐波检测任务正常执行，测量任务也正常运行，只是不进行测量值计算，测量幅值结果直接写入“0”值；当计算缓冲区填充一个完整周波采样数据时，测量幅值正常计算。
一种用于电力系统高次谐波抑制的测控装置，包括：信号采样模块、频率计算模块、采样数据转换模块、谐波检测模块、测量值计算模块；
所述信号采样模块以设定采样频率对模拟信号进行采样，得到采样信号，每周波采样点数为40点；
所述频率计算模块是基于信号采样模块输出的采样数据，采用过零点算法计算系统工频，为测量值计算模块提供判断依据；
所述采样数据转换模块是基于信号采样模块输出的采样数据，进行重采样存入计算缓冲区，此时每周波采样点数大于等于80点，为谐波检测模块提供检测数据，并为测量值计算模块提供计算数据；
所述谐波检测模块以采样数据转换模块输出的重采样数据作为输入，经过80点谐波检测算法，检测各次谐波分量大小，从而确定是否存在各次高次谐波，检测结果用于计算电压、电流和功率等测量值；
所述测量值计算模块依据谐波检测模块的检测结果，按照测量幅值计算方案，采用40点和80点采样数据结合的方式，并基于傅里叶算法，精确计算测量值。
综上所述，与现有技术相比，上述技术方案的有益效果是：
通过离散傅里叶变换算法可将电压和电流信号按照基波频率的整数倍进行分解，分解出基波和各次谐波分量。能够对谐波频率分量进行较为精确的测量操作，能够对高次谐波进行有效的检测，极大程度地提高谐波频率分量的测量精度，进而有助于维持电力系统的稳定运行。
附图说明
图1为本申请实施例的高次谐波抑制测控装置流程图；
图2为本申请实施例的检测十三次谐波正、余弦系数波形图；
图3为本申请实施例的电压信号叠加谐波仿真波形图；
图4为本申请实施例的谐波检测各次谐波分量分布图。
具体实施方式
以下结合附图对本申请的原理和特征进行描述，所举实例只用于解释本申请，并非用于限定本申请的范围。
本申请实施例公开一种用于电力系统高次谐波抑制的测控装置及方法。
本申请主要基于离散傅里叶级数理论检测和抑制谐波分量，针对于离散傅里叶级 数，首先，一个周期的离散时间信号
Figure 219811DEST_PATH_IMAGE023
须满足：
Figure 291673DEST_PATH_IMAGE024
（1）
在式（1）中，N为正整数，表示
Figure 341668DEST_PATH_IMAGE025
的周期。例如，复指数序列
Figure 473572DEST_PATH_IMAGE026
是一个周期序 列，周期为N，基波频率
Figure 343439DEST_PATH_IMAGE027
为：
Figure 218992DEST_PATH_IMAGE028
（2）
由该复指数序列可知，呈谐波关系的复指数序列集为：
Figure 389073DEST_PATH_IMAGE029
（3）
其中，k=0，±1，±2……
每个分量
Figure 301665DEST_PATH_IMAGE030
的频率是基波频率
Figure 49041DEST_PATH_IMAGE027
的整数倍。
同时，
Figure 335142DEST_PATH_IMAGE031
也是周期序列，公式如下：
Figure 484363DEST_PATH_IMAGE032
（4）
其中，r为整数。
即当k变化了一个N的整数倍时，能够得到一个完全相同的序列。所以，基波周期为 N的周期序列
Figure 567857DEST_PATH_IMAGE033
，可用N个呈谐波关系的复指数序列的加权和进行表示，公式如下：
Figure 536950DEST_PATH_IMAGE034
（5）
上式（5）求和限
Figure 629671DEST_PATH_IMAGE035
表示求和仅需要包含N项，k既可以取0，1，2……N-1，也 可以取2，3，4……N+1，以此类推。无论k如何取值，由于公式（4）的周期性，公式（5）右边求和 结果都是相同的。
将周期序列按照公式（5）的形式进行表示，即一组呈谐波关系的复指数序列的加 权和，称为离散傅里叶级数。其中，系数
Figure 508765DEST_PATH_IMAGE036
称为离散傅里叶系数，公式如下：
Figure 153373DEST_PATH_IMAGE037
（6）
上式（5）、（6）是周期离散时间信号
Figure 485128DEST_PATH_IMAGE038
和其傅里叶系数
Figure 240595DEST_PATH_IMAGE036
相互转换的公式。其 中，公式（5）可以根据计算出的
Figure 239775DEST_PATH_IMAGE036
值还原出原始信号
Figure 55284DEST_PATH_IMAGE038
，而公式（6）可以根据已知信号
Figure 608756DEST_PATH_IMAGE038
分析提取出信号中所含的频谱。
另外，傅里叶系数
Figure 167914DEST_PATH_IMAGE036
称为信号
Figure 756021DEST_PATH_IMAGE038
的频谱系数。这些频谱系数说明信号
Figure 742432DEST_PATH_IMAGE039
可被 分解成N个呈谐波关系的复指数序列的加权和，本发明谐波分量的检测和测量值计算的原 理主要来自于此。
本申请实施例公开一种用于电力系统高次谐波抑制的测控方法。
实施例1
一种用于电力系统高次谐波抑制的测控装置及方法，包括：
通过离散傅里叶变换算法可将电压和电流信号按照基波频率的整数倍进行分解， 分解出基波和各次谐波分量。根据上述傅里叶级数理论，利用公式（6）可推导出，是基波频 率k倍的频率系数
Figure 317288DEST_PATH_IMAGE036
。
根据离散傅里叶变换算法，得出频率系数
Figure 289923DEST_PATH_IMAGE001
，公式如下：
Figure 857171DEST_PATH_IMAGE002
（7）
公式（7）中，u(n)表示已知的电压输入信号，N为周期采样点数，其中，0≤k≤（N-1）；
利用欧拉公式展开所述频率系数
Figure 889849DEST_PATH_IMAGE001
，得到所述频率系数
Figure 808126DEST_PATH_IMAGE001
的实部
Figure 584452DEST_PATH_IMAGE003
和虚部
Figure 412731DEST_PATH_IMAGE004
，公式如（8）和（9）所示：
Figure 740944DEST_PATH_IMAGE005
（8）
Figure 756305DEST_PATH_IMAGE006
（9）
其中，N为周期采样点数，k为谐波次数，k=1，2……13；
计算得出k个频率分量的幅值
Figure 460956DEST_PATH_IMAGE007
，即电压的基波分量和各次谐波分量的幅值，公 式如下：
Figure 878162DEST_PATH_IMAGE008
（10）
其中，k为谐波次数，k=1，2……13；
根据所述幅值
Figure 377276DEST_PATH_IMAGE007
，计算得到基波和2-13次谐波的幅值大小；
当2-13次中的某次谐波幅值超过设定阈值后，认定存在该次谐波分量，根据所述 幅值
Figure 145512DEST_PATH_IMAGE007
，计算得到电压的有效值
Figure 526290DEST_PATH_IMAGE009
，公式如下：
Figure 188216DEST_PATH_IMAGE010
（11）
根据所述电压的有效值
Figure 733598DEST_PATH_IMAGE009
，计算得到电压和电流的有效值；
对于13次以上的谐波分量，不进行检测，不加入有效值计算，实现高次谐波抑制作用；
同理，根据离散傅里叶变换算法以及欧拉公式，计算得到电流频率系数的实部
Figure 582605DEST_PATH_IMAGE011
和虚部
Figure 504425DEST_PATH_IMAGE012
，得出电流的基波和各次谐波幅值
Figure 161802DEST_PATH_IMAGE013
，其中，计算电流有效值
Figure 2719DEST_PATH_IMAGE014
，公式如 下：
Figure 214389DEST_PATH_IMAGE015
（12）
根据有功功率P和无功功率Q的计算公式：
Figure 798954DEST_PATH_IMAGE016
（13）
Figure 310838DEST_PATH_IMAGE017
（14）
分别得出电力信号有功功率P和无功功率Q的推导公式：
Figure 463601DEST_PATH_IMAGE018
（15）
Figure 287201DEST_PATH_IMAGE019
（16）。
实现了电力系统在电网电力电子化引发的高次谐波严重影响下对电压、电流和功率等遥测值的精准测量，填补了测控装置谐波抑制的算法空缺。基于傅氏算法，算法简单易实现，并且有效地提高了高次谐波影响下工频遥测值的测量精度。
实施例1的实施原理为：通过离散傅里叶变换算法可将电压和电流信号按照基波频率的整数倍进行分解，分解出基波和各次谐波分量。通过利用谐波检测算法，能够对谐波频率分量进行较为精确的测量操作，能够对高次谐波进行有效的检测，极大程度地提高谐波频率分量的测量精度，进而有助于维持电力系统的稳定运行。
实施例2
本实施例与实施例1的不同之处在于，数据采样形式不同。
根据奈奎斯特采样定理：fs≥2fc（19）
可知，当采样频率fs小于信号最高频率的2倍时，就会产生频率混叠现象。就目前 叠加谐波次数最大可到40次，而采样数据只有40点（即采样频率为基波频率
Figure 550823DEST_PATH_IMAGE040
的40倍，信 号最大频率为基波频率
Figure 41847DEST_PATH_IMAGE040
的40倍），未大于或等于
Figure 380161DEST_PATH_IMAGE040
的80倍，显然不满足奈奎斯特定理，必 将产生频率混叠。因此，至少要选择80点采样数据，才能避免频率混叠，防止错误检测和抑 制20次以上的高次谐波。
考虑CPU负载率和频率混叠的影响，本申请实施例提出40点和80点采样数据结合进行谐波检测和测量值计算的方案。其中，当k取值为2-13时，基波幅值和谐波幅值采用40点采样数据进行计算，电压、电流有效值等于基波幅值和谐波幅值之和；当k取值为14-38或40-50时，谐波幅值不参与有效值计算，基波幅值采用40点采样数据进行计算，电压、电流有效值等于基波幅值；当k取值大于等于39时，谐波幅值不参与有效值计算，基波幅值采用80点采样数据进行计算，电压、电流有效值等于基波幅值。
为了验证该方案对谐波干扰信号的检测和抑制效果，使用昂立测试仪模拟电力系统的电压、电流信号，叠加上2-50次谐波信号进行功能测试。测试仪施加的基波电压幅值为57.735V，谐波电压幅值是基波电压幅值的20%，设为11.547V，基波电流幅值设为5A，谐波电流为1A。理论上，施加2-13次谐波时测量有效值为基波与谐波的幅值叠加，施加13次以上谐波时被有效抑制同时测量有效值为基波幅值，电压、电流和功率的理论值下表已给出。表1为施加2-13次谐波的最大误差对比情况，表2为施加14-50次高次谐波的最大误差对比情况。
表1为40点和80点采样数据结合施加2~13次谐波最大误差对比表
Figure 691057DEST_PATH_IMAGE041
表2为40点和80点采样数据结合施加14~50次高次谐波最大误差对比表
Figure 492790DEST_PATH_IMAGE042
测控装置技术规范要求，电压和电流的测量误差应小于0.2%，功率测量误差应小于0.5%，叠加谐波后允许的误差改变量不应超过200％，即叠加谐波后电压和电流测量误差应小于0.4%，功率测量误差应小于1%。从测试数据可以看出，叠加2-13次谐波时，电压、电流测量误差小于0.4%，功率测量误差小于1%；叠加14-50次谐波时，电压、电流测量误差小于0.4%，功率测量误差小于1%。以上不同测试情况下的测量误差均满足规范误差要求。
实施例2的实施原理为：采用40点和80点采样相结合的检测和计算方式，减少硬件配置的限制，可使计算量大大降低。不仅良好适应装置的硬件性能，同时防止CPU负载率过高，实现多功能测控装置计算量的精准度，当装置采用40点采样数据进行测量值计算时，能够保证装置所有功能正常运行。
实施例3
本实施例与实施例1的不同之处在于，数据采样形式不同，具体方案如下：谐波分量全部采用80点采样数据进行检测，电压和电流幅值同样采用80点采样数据进行计算。谐波分量全部采用80点采样数据进行检测，电压和电流幅值同样采用80点采样数据进行计算，根据谐波检测算法：
Figure 838321DEST_PATH_IMAGE020
（8）
Figure 67308DEST_PATH_IMAGE021
（9）
其中，k为谐波次数，k=1，2……13；
分解得到电压和电流信号的基波和2-13次谐波分量，再根据
Figure 6445DEST_PATH_IMAGE022
（10）
计算得到基波和2-13次谐波的幅值大小；
当2-13次中的某次谐波幅值超过一定阈值后，认定存在该次谐波分量，并通过公式
Figure 2083DEST_PATH_IMAGE043
（11）
计算得到电压和电流的有效值；
对于13次以上的谐波分量，不进行检测，不加入有效值计算，实现高次谐波抑制作用。
当装置硬件性能提高，CPU负载率不受限制时，本申请实施例提出80点谐波检测、抑制以及80点测量值计算的方案。采用昂立测试仪向测控装置模拟输入电压和电流信号，并叠加2-50次谐波进行功能测试。测试仪施加的电压、电流的基波和谐波信号同上述实施例2一致。表3为施加2-13次谐波的最大误差对比情况，表4为施加14-50次高次谐波的最大误差对比情况。
表3为80点谐波检测时施加2-13次谐波的最大误差对比表
Figure 77487DEST_PATH_IMAGE044
表4为80点谐波检测时施加14-50次高次谐波最大误差对比表
Figure 602009DEST_PATH_IMAGE045
测控装置技术规范对电压、电流和功率的测量误差要求同上述实施例2一致。从测试结果可以看出，当叠加2-13次谐波时，电压、电流测量误差小于0.02%，不大于0.2%，功率测量误差小于0.03%，不大于1%；叠加14-50次谐波时，电压、电流测量误差小于0.01%，不大于0.2%，功率测量误差小于0.03%，不大于1%。以上不同测试情况下的测量误差均在规范误差要求范围内。
实施例3的实施原理为：当装置的硬件配置较高，且CPU负载率较低时，可支持测控 装置的高计算量，能够保证装置所有功能正常运行。因此，装置可以采用80点采样数据进行 测量值计算。由奈奎斯特定理可知，当采样频率为基波频率
Figure 762863DEST_PATH_IMAGE046
的80倍时，能够较为准确地分 解出2-40次的谐波分量，而不发生频率混叠，满足技术要求。
实施例4
本实施例与实施例1的不同之处在于，幅值跳变解决方案不同，当测试仪向测控装置初始加入遥测量的瞬间，由于计算缓冲区在谐波检测任务执行时刻未及时填入完整的一周期采样数据，导致谐波检测算法发生误判，不能准确检测出相应的谐波分量，进而产生幅值跳变。如果不加以处理，会导致测量值计算异常，影响电力系统的安全运行。针对这一现象，本申请实施例提出一种抑制幅值跳变产生计算误差的解决方案。
导致上述异常问题的主要原因在于，谐波检测时计算缓冲区不是一个完整的周波数据，算法不能准确判断。因此，可在计算缓冲区采样数据填充完成之前，让谐波检测任务正常执行，测量任务也正常运行，只是不进行测量值计算，测量幅值结果直接写入“0”值；在计算缓冲区填充好一个完整周波采样数据后，测量幅值才正常计算。
实施例4的实施原理为：消除幅值跳变产生的测量误差影响可使电力系统安全、稳定运行，提高了系统的安全性和稳定性。
参照图1，一种用于电力系统高次谐波抑制的测控装置，包括：信号采样模块、频率计算模块、采样数据转换模块、谐波检测模块、测量值计算模块；
所述信号采样模块由采样脉冲对模拟信号进行采样，得到用于频率计算和采样数据转换的采样信号，每周波采样点数为40点；
所述频率计算模块是基于信号采样模块输出的采样数据，采用过零点算法计算系统侧频率值，为测量值计算模块提供判断依据；
所述采样数据转换模块是基于信号采样模块输出的采样数据，进行重采样存入计算缓冲区，此时每周波采样点数为80点，为谐波检测模块提供检测数据，并为测量值计算模块提供计算数据；
所述谐波检测模块以采样数据转换模块输出的80点采样数据作为输入，经过80点谐波检测算法，检测各次谐波分量大小，从而确定是否存在高次谐波，检测结果用于计算电压、电流和功率等测量值；
所述测量值计算模块依据谐波检测模块的检测结果，按照测量幅值计算方案，采用40点和80点采样数据结合的方式，并基于离散傅里叶算法，精确计算测量值。
在谐波检测过程中，谐波检测算法的正、余弦系数事先通过MATLAB工具生成。根据公式（8）和（9），余弦系数和正弦系数的MATLAB生成代码分别为：
cos(n, k) = sqrt(2)/N * cos(n*k*2*pi/N)（17）
sin(n, k) = sqrt(2)/N * sin(n*k*2*pi/N)（18）
参照图2，其中，N是周期采样点个数，等于80；n是采样点序号，取值为1-N；k是谐波次数，取值为1-13。谐波检测算法生成的正、余弦系数波形图如附图2所示。然后，将生成的正、余弦系数cos(n, k)和sin(n, k)经过一定倍数的放大，作为实际计算系数与测控装置采样数据进行对应采样点的乘积、求和运算。
于是，可得到基波和2至13次谐波分量的实部、虚部。然后，经公式（10）可计算获得基波和2至13次谐波分量的幅值大小，再根据公式（11）和（12）即可计算出电压、电流的有效值。最后，由公式（14）、（15）得到系统的有功功率和无功功率。
参照图3和图4，在仿真时，周期采样点数N设置为80点，采样频率fs设为50Hz，电压幅值Ua为57.735V，2次谐波和13次谐波幅值设置为基波幅值的20%，均为11.547V。附图3和附图4分别给出了电压信号叠加谐波仿真波形图及谐波检测各次谐波分量分布图。其中，附图3模拟正弦电压信号受到2次和13次谐波干扰后的波形情况，附图4是经过谐波检测算法分解出的基波和谐波分量。
可见，系统能够分解得出了基波和2至13次谐波分量，且各次谐波幅值Uk=[57.7350000000000,11.5470000000000,1.70986060406632e-14,1.55510504872441e-14,1.48045844976432e-14,5.55111512312578e-15,1.15484543030194e-14,3.22462457067851e-14,3.33364684119621e-14,2.78468253768464e-14,3.88226593384640e-14,1.67467124462598e-14,11.5470000000000]，其中，k为谐波次数，取值1-13。
以上所述仅为本申请的较佳实施例，并不用以限制本申请，凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。

Claims (7)

1.一种用于电力系统高次谐波抑制的测控方法，其特征在于，包括：
根据离散傅里叶变换算法，得出频率系数
Figure 163604DEST_PATH_IMAGE001
，公式如下：
Figure 778257DEST_PATH_IMAGE002
（7）
公式（7）中，u(n)表示已知的电压输入信号，N为周期采样点数，其中，0≤k≤（N-1）；
利用欧拉公式展开所述频率系数
Figure 747350DEST_PATH_IMAGE001
，得到所述频率系数
Figure 574491DEST_PATH_IMAGE001
的实部
Figure 578219DEST_PATH_IMAGE003
和虚部
Figure 95264DEST_PATH_IMAGE004
，公 式如（8）和（9）所示：
Figure 551653DEST_PATH_IMAGE005
（8）
Figure 182486DEST_PATH_IMAGE006
（9）
其中，N为周期采样点数，k为谐波次数，k=1，2……13；
计算得出k个频率分量的幅值
Figure 40720DEST_PATH_IMAGE007
，即电压的基波分量和各次谐波分量的幅值，公式如 下：
Figure 997175DEST_PATH_IMAGE008
（10）
其中，k为谐波次数，k=1，2……13；
根据所述幅值
Figure 81806DEST_PATH_IMAGE007
，计算得到基波和2-13次谐波的幅值大小；
当2-13次中的某次谐波幅值超过设定阈值后，认定存在该次谐波分量，根据所述幅值
Figure 375384DEST_PATH_IMAGE007
，计算得到电压的有效值
Figure 229070DEST_PATH_IMAGE009
，公式如下：
Figure 90847DEST_PATH_IMAGE010
（11）
根据所述电压的有效值
Figure 787407DEST_PATH_IMAGE009
，计算得到电压和电流的有效值；
对于13次以上的谐波分量，不进行检测，不加入有效值计算，实现高次谐波抑制作用；
同理，根据离散傅里叶变换算法以及欧拉公式，计算得到电流频率系数的实部
Figure 760043DEST_PATH_IMAGE011
和 虚部
Figure 468236DEST_PATH_IMAGE012
，得出电流的基波和各次谐波幅值
Figure 625548DEST_PATH_IMAGE013
，其中，计算电流有效值
Figure 675981DEST_PATH_IMAGE014
，公式如下：
Figure 311362DEST_PATH_IMAGE015
（12）
根据有功功率P和无功功率Q的计算公式：
Figure 874062DEST_PATH_IMAGE016
（13）
Figure 343220DEST_PATH_IMAGE017
（14）
分别得出电力信号有功功率P和无功功率Q的推导公式：
Figure 217635DEST_PATH_IMAGE018
（15）
Figure 63231DEST_PATH_IMAGE019
（16）。
2.根据权利要求1所述的一种用于电力系统高次谐波抑制的测控方法，其特征在于：当k取值为2-13时，基波幅值和谐波幅值采用40点采样数据进行计算，电压、电流有效值等于基波幅值和谐波幅值之和。
3.根据权利要求1所述的一种用于电力系统高次谐波抑制的测控方法，其特征在于：当k取值为14-38或40-50时，谐波幅值不参与有效值计算，基波幅值采用40点采样数据进行计算，电压、电流有效值等于基波幅值。
4.根据权利要求1所述的一种用于电力系统高次谐波抑制的测控方法，其特征在于：当k取值大于等于39时，谐波幅值不参与有效值计算，基波幅值采用80点采样数据进行计算，电压、电流有效值等于基波幅值。
5.根据权利要求1所述的一种用于电力系统高次谐波抑制的测控方法，其特征在于：谐波分量全部采用80点采样数据进行检测，电压和电流幅值同样采用80点采样数据进行计算，根据谐波检测算法：
Figure 870650DEST_PATH_IMAGE020
（8）
Figure 245131DEST_PATH_IMAGE021
（9）
其中，k为谐波次数，k=1，2……13；
分解得到电压和电流信号的基波和2-13次谐波分量，再根据
Figure 13367DEST_PATH_IMAGE022
（10）
计算得到基波和2-13次谐波的幅值大小。
6.根据权利要求1所述的一种用于电力系统高次谐波抑制的测控方法，其特征在于：计算缓冲区采样数据填充完成之前，谐波检测任务正常执行，测量任务也正常运行，只是不进行测量值计算，测量幅值结果直接写入“0”值；当计算缓冲区填充一个完整周波采样数据时，测量幅值正常计算。
7.一种用于电力系统高次谐波抑制的测控装置，其特征在于，包括：信号采样模块、频率计算模块、采样数据转换模块、谐波检测模块、测量值计算模块；
所述信号采样模块以设定采样频率对模拟信号进行采样，得到采样信号，每周波采样点数为40点；
所述频率计算模块是基于信号采样模块输出的采样数据，采用过零点算法计算系统工频，为测量值计算模块提供判断依据；
所述采样数据转换模块是基于信号采样模块输出的采样数据，进行重采样存入计算缓冲区，此时每周波采样点数大于等于80点，为谐波检测模块提供检测数据，并为测量值计算模块提供计算数据；
所述谐波检测模块以采样数据转换模块输出的重采样数据作为输入，经过80点谐波检测算法，检测各次谐波分量大小，从而确定是否存在各次高次谐波，检测结果用于计算电压、电流和功率等测量值；
所述测量值计算模块依据谐波检测模块的检测结果，按照测量幅值计算方案，采用40点和80点采样数据结合的方式，并基于傅里叶算法，精确计算测量值。
CN202211314350.6A 2022-10-26 2022-10-26 一种用于电力系统高次谐波抑制的测控装置及方法 Pending CN115372698A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211314350.6A CN115372698A (zh) 2022-10-26 2022-10-26 一种用于电力系统高次谐波抑制的测控装置及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211314350.6A CN115372698A (zh) 2022-10-26 2022-10-26 一种用于电力系统高次谐波抑制的测控装置及方法

Publications (1)

Publication Number Publication Date
CN115372698A true CN115372698A (zh) 2022-11-22

Family

ID=84074358

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211314350.6A Pending CN115372698A (zh) 2022-10-26 2022-10-26 一种用于电力系统高次谐波抑制的测控装置及方法

Country Status (1)

Country Link
CN (1) CN115372698A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115754433A (zh) * 2023-01-09 2023-03-07 北京智芯微电子科技有限公司 特征电流编码识别方法、装置及网络拓扑识别方法
CN116223905A (zh) * 2022-12-14 2023-06-06 广州汇电云联互联网科技有限公司 谐波频率检测系统、方法、装置及存储介质

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101581752A (zh) * 2009-06-09 2009-11-18 深圳市双合电脑系统股份有限公司 一种基于故障录波的输电线路参数在线测量方法
CN202837498U (zh) * 2012-10-19 2013-03-27 重庆飞尔达机械有限责任公司 煤矿用馈电开关检测装置
CN105067882A (zh) * 2015-07-24 2015-11-18 许继电气股份有限公司 一种电量幅值的测定方法
CN109358223A (zh) * 2018-09-29 2019-02-19 海特尔机电工程技术（马鞍山）有限公司 一种滑动窗dft谐波电流检测方法及装置、存储介质
CN110602247A (zh) * 2019-09-27 2019-12-20 国网河南省电力公司博爱县供电公司 一种输电线路异物和缺失智能视频在线监测平台
CN110988470A (zh) * 2019-12-17 2020-04-10 国网江苏省电力有限公司检修分公司 一种基于自适应电网次谐波提取与控制方法
CN111190042A (zh) * 2019-12-30 2020-05-22 中国电力科学研究院有限公司 一种智能传感终端及对电网信号全频段测量的方法
CN114337204A (zh) * 2021-12-02 2022-04-12 北京科技大学顺德研究生院 一种具有低开关频率特性的预测控制指定谐波抑制开关策略

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101581752A (zh) * 2009-06-09 2009-11-18 深圳市双合电脑系统股份有限公司 一种基于故障录波的输电线路参数在线测量方法
CN202837498U (zh) * 2012-10-19 2013-03-27 重庆飞尔达机械有限责任公司 煤矿用馈电开关检测装置
CN105067882A (zh) * 2015-07-24 2015-11-18 许继电气股份有限公司 一种电量幅值的测定方法
CN109358223A (zh) * 2018-09-29 2019-02-19 海特尔机电工程技术（马鞍山）有限公司 一种滑动窗dft谐波电流检测方法及装置、存储介质
CN110602247A (zh) * 2019-09-27 2019-12-20 国网河南省电力公司博爱县供电公司 一种输电线路异物和缺失智能视频在线监测平台
CN110988470A (zh) * 2019-12-17 2020-04-10 国网江苏省电力有限公司检修分公司 一种基于自适应电网次谐波提取与控制方法
CN111190042A (zh) * 2019-12-30 2020-05-22 中国电力科学研究院有限公司 一种智能传感终端及对电网信号全频段测量的方法
CN114337204A (zh) * 2021-12-02 2022-04-12 北京科技大学顺德研究生院 一种具有低开关频率特性的预测控制指定谐波抑制开关策略

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘 毅: "差分与傅氏变换相结合算法在微机电流保护中的应用" *
张晔: "《信号与系统》", 31 December 2020 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116223905A (zh) * 2022-12-14 2023-06-06 广州汇电云联互联网科技有限公司 谐波频率检测系统、方法、装置及存储介质
CN116223905B (zh) * 2022-12-14 2023-12-08 广州市均能科技有限公司 谐波频率检测系统、方法、装置及存储介质
CN115754433A (zh) * 2023-01-09 2023-03-07 北京智芯微电子科技有限公司 特征电流编码识别方法、装置及网络拓扑识别方法

Similar Documents

Publication Publication Date Title
CN115372698A (zh) 一种用于电力系统高次谐波抑制的测控装置及方法
Asiminoaei et al. A digital controlled PV-inverter with grid impedance estimation for ENS detection
Wen et al. Simple interpolated FFT algorithm based on minimize sidelobe windows for power-harmonic analysis
CN104897960B (zh) 基于加窗四谱线插值fft的谐波快速分析方法及系统
CN113064021B (zh) 用于实现电力电子化电网高次谐波抑制的测控装置和方法
Tarasiuk Comparative study of various methods of DFT calculation in the wake of IEC Standard 61000-4-7
Roscoe et al. Frequency and fundamental signal measurement algorithms for distributed control and protection applications
Kocewiak et al. The impact of harmonics calculation methods on power quality assessment in wind farms
Kusljevic et al. Simultaneous reactive-power and frequency estimations using simple recursive WLS algorithm and adaptive filtering
CN104007408A (zh) 一种pmu动态性能的在线检测方法及装置
CN109444525B (zh) 变电站测控装置交流采样数据有效值及功率因数计算方法
CN113358922B (zh) 非稳态电网信号的基波、谐波、间谐波电能功率测量方法
Radil et al. Methods for estimation of voltage harmonic components
CN103592514A (zh) 一种新型谐波高精度检测方法
Redondo et al. A strategy for improving the accuracy of flicker emission measurement from wind turbines
Nguyen et al. Application of a z-transform signal model and median filtering for power system frequency and phasor measurements
CN112180161A (zh) 一种非同步高采样率采样条件下谐波间谐波群测量方法
Jia et al. A novel adaptive linear prediction-based parameter estimation method for monitoring sub-/inter-harmonics during SSI events
Luo et al. Research and design of data acquisition and processing algorithm based on improved FFT
Rabehi et al. Parameters Selection Criteria of Prony Method for Accurate Harmonics and Inter-Harmonic Component Identification.
Petrović et al. Digital processing of Synchronously Sampled AC Signals in Presence of Interharmonics and Subharmonics
Zhang et al. Study of harmonic analysis based on improved discrete Fourier transform
Xiong et al. A Multiple-Frequency-Taylor-Model based estimator for dynamic synchrophasor considering decaying DC component
Yang Method of measuring power quality and development of monitor device
Zhang et al. Research on high-speed ADC test and evaluation system

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20221122

RJ01 Rejection of invention patent application after publication