CN114460366A - Power grid harmonic phase analysis method, system and device - Google Patents

Power grid harmonic phase analysis method, system and device Download PDF

Info

Publication number
CN114460366A
CN114460366A CN202210093206.8A CN202210093206A CN114460366A CN 114460366 A CN114460366 A CN 114460366A CN 202210093206 A CN202210093206 A CN 202210093206A CN 114460366 A CN114460366 A CN 114460366A
Authority
CN
China
Prior art keywords
phase
harmonic
voltage
current
signal
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
CN202210093206.8A
Other languages
Chinese (zh)
Other versions
CN114460366B (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.)
Shenzhen Power Supply Bureau Co Ltd
Original Assignee
Shenzhen Power Supply Bureau 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 Shenzhen Power Supply Bureau Co Ltd filed Critical Shenzhen Power Supply Bureau Co Ltd
Priority to CN202210093206.8A priority Critical patent/CN114460366B/en
Publication of CN114460366A publication Critical patent/CN114460366A/en
Application granted granted Critical
Publication of CN114460366B publication Critical patent/CN114460366B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/40Arrangements for reducing harmonics

Landscapes

  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Remote Monitoring And Control Of Power-Distribution Networks (AREA)

Abstract

本发明涉及一种电网谐波相位分析方法、系统和装置。所述方法包括:获取关注母线的电压波形信号和电力用户馈线上的电流波形信号;将统计时间划分为多个计算窗口;且在每个计算窗口内,对所述电压波形信号和电流波形信号截取预设数量个周波;对所述预设数量个周波进行快速傅里叶变换,得到每相电压信号的谐波相量和每相电流信号的谐波相量;根据所述每相电压信号的谐波相量获取每相基波电压上升沿的过零点时刻,并根据所述每相电流信号的谐波相量获取每相电流信号各次谐波在各相所述过零点时刻的相位;基于在统计时间内获得的所有相位数据进行分析。上述电网谐波相位分析方法、系统和装置能够更准确地反映谐波电流相位分布,且兼具有相序容错能力。

Figure 202210093206

The invention relates to a power grid harmonic phase analysis method, system and device. The method includes: acquiring the voltage waveform signal of the bus of interest and the current waveform signal on the power user feeder; dividing the statistical time into a plurality of calculation windows; and in each calculation window, analyzing the voltage waveform signal and the current waveform signal Intercept a preset number of cycles; perform fast Fourier transform on the preset number of cycles to obtain a harmonic phasor of each phase voltage signal and a harmonic phasor of each phase current signal; according to the each phase voltage signal The harmonic phasor of each phase obtains the zero-crossing moment of the rising edge of the fundamental voltage of each phase, and obtains the phase of each harmonic of each phase current signal at the zero-crossing moment of each phase according to the harmonic phasor of the current signal of each phase ; The analysis is based on all phase data acquired within the statistical time. The above-mentioned power grid harmonic phase analysis method, system and device can more accurately reflect the harmonic current phase distribution, and also have the capability of phase sequence fault tolerance.

Figure 202210093206

Description

电网谐波相位分析方法、系统和装置Power grid harmonic phase analysis method, system and device

技术领域technical field

本发明涉及电网的电能质量分析技术领域,特别是涉及一种电网谐波相位分析方法、系统和装置。The invention relates to the technical field of power quality analysis of power grids, in particular to a method, system and device for analyzing harmonic phases of power grids.

背景技术Background technique

现代电力供应要求具有高可靠性、控制灵活、应用方便,但电力系统污染导致电能质量不佳的问题一直存在:一方面,随着各种新型用户负荷,尤其是具有非线性和会产生谐波的各种电力电子装置的大规模应用,给电网电能质量带来不良影响;另一方面,用户对供电可靠性的要求越来越高,特别是众多精密电子设备在电力系统大量使用,对供电质量的敏感程度越来越高。因此,亟需加强配电网谐波管理。Modern power supply requires high reliability, flexible control and convenient application, but the problem of poor power quality caused by power system pollution has always existed: on the one hand, with various new user loads, especially nonlinear and harmonic The large-scale application of various power electronic devices has brought adverse effects on the power quality of the power grid; on the other hand, users have higher and higher requirements for the reliability of power supply, especially the large-scale use of many precision electronic equipment in the power system. Quality is becoming more and more sensitive. Therefore, it is urgent to strengthen the harmonic management of distribution network.

配电网的谐波管理是现代电力系统的重要基础工作,是电力系统安全、经济、稳定运行的保证。谐波管理的主要内容包括负荷调查、现场测试、分析并提出治理建议和制订并实施谐波管理办法。其主要包括:1)通过实地调查和测试得到相应的谐波数据;2)利用相应的技术手段进行谐波分析;3)根据分析结果提出治理措施等。Harmonic management of distribution network is an important basic work of modern power system, and it is the guarantee of safe, economical and stable operation of power system. The main contents of harmonic management include load investigation, on-site test, analysis and put forward management suggestions, and formulate and implement harmonic management methods. It mainly includes: 1) Obtaining corresponding harmonic data through on-site investigation and testing; 2) Using corresponding technical means to carry out harmonic analysis; 3) Proposing control measures according to the analysis results.

而谐波分析中较为重要的一环为谐波相位分布特性的分析。目前利用谐波阻抗角表征谐波电流的相位特征,但其存在以下问题:The most important part of harmonic analysis is the analysis of harmonic phase distribution characteristics. At present, the harmonic impedance angle is used to characterize the phase characteristics of the harmonic current, but it has the following problems:

(1)由于关注点的谐波电压受背景谐波电压和负荷谐波电流的共同影响,采用谐波电压作为典型负荷谐波电流相位分布特性的参考基准会使得谐波电流的相位分布出现较大的偏差,因此谐波阻抗角并不能准确地表征谐波电流的相位特征。(1) Since the harmonic voltage of the point of interest is jointly affected by the background harmonic voltage and the load harmonic current, using the harmonic voltage as the reference for the phase distribution characteristics of the typical load harmonic current will make the phase distribution of the harmonic current appear relatively high. Therefore, the harmonic impedance angle cannot accurately characterize the phase characteristics of the harmonic current.

(2)谐波阻抗角不能具备相序容错能力。不同低压用电负荷中的ABC三相未必是完全对应的,即可能出现一台配变的ABC三相与另外一台配变的BCA三相是对应的,若谐波电流相位分布不具备相序容错能力,若错将B相数据当作A相数据,可能会得出错误的结果。(2) The harmonic impedance angle cannot have phase sequence fault tolerance. The ABC three-phase in different low-voltage electrical loads may not be completely corresponding, that is, the ABC three-phase of one distribution transformer may correspond to the BCA three-phase of another distribution transformer. Sequence fault tolerance, if the B-phase data is mistakenly regarded as the A-phase data, wrong results may be obtained.

发明内容SUMMARY OF THE INVENTION

基于此,针对采用谐波阻抗角表征相位谐波电流的相位特征存在的问题,提供一种电网谐波相位分析方法、系统和装置。Based on this, in order to solve the problem of using the harmonic impedance angle to characterize the phase characteristics of the phase harmonic current, a method, system and device for analyzing the harmonic phase of a power grid are provided.

一种电网谐波相位分析方法,所述的方法包括:A power grid harmonic phase analysis method, the method includes:

获取关注母线的电压波形信号和电力用户馈线上的电流波形信号;Obtain the voltage waveform signal of the concerned bus and the current waveform signal of the power user feeder;

将统计时间划分为多个计算窗口;且在每个计算窗口内,对所述电压波形信号和电流波形信号截取预设数量个周波;dividing the statistical time into a plurality of calculation windows; and in each calculation window, intercepting a preset number of cycles from the voltage waveform signal and the current waveform signal;

对所述预设数量个周波进行快速傅里叶变换,得到每相电压信号的谐波相量和每相电流信号的谐波相量;performing fast Fourier transform on the preset number of cycles to obtain the harmonic phasor of each phase voltage signal and the harmonic phasor of each phase current signal;

根据所述每相电压信号的谐波相量获取每相基波电压上升沿的过零点时刻,并以过零点时刻的相位作为参考基准,根据所述每相电流信号的谐波相量获取每相电流信号各次谐波在各相所述过零点时刻的相位;Obtain the zero-crossing time of the rising edge of the fundamental voltage of each phase according to the harmonic phasor of the voltage signal of each phase, and take the phase at the zero-crossing time as a reference, and obtain each phase according to the harmonic phasor of the current signal of each phase. The phase of each harmonic of the phase current signal at the zero-crossing point of each phase;

基于在统计时间内获得的所有相位数据进行分析。Analysis is based on all phase data acquired within statistical time.

在其中一个实施例中,所述计算窗口在统计时间内按固定时长等分。In one of the embodiments, the calculation window is divided equally into a fixed duration within the statistical time.

在其中一个实施例中,所述获取关注母线的电压波形信号和电力用户馈线上的电流波形信号,包括:In one embodiment, the acquiring the voltage waveform signal of the concerned bus and the current waveform signal of the power user feeder includes:

利用电压互感器、电流互感器配合数据采集与监视控制系统获取所述电压波形信号和电流波形信号。The voltage waveform signal and the current waveform signal are obtained by using the voltage transformer and the current transformer in cooperation with the data acquisition and monitoring control system.

在其中一个实施例中,所述电压波形信号和电流波形信号每个通道的采样频率不小于12.8kHz。In one embodiment, the sampling frequency of each channel of the voltage waveform signal and the current waveform signal is not less than 12.8 kHz.

在其中一个实施例中,所述预设数量为10。In one of the embodiments, the preset number is 10.

在其中一个实施例中,所述根据所述每相电压信号的谐波相量获取每相基波电压上升沿的过零点时刻,包括:In one embodiment, obtaining the zero-crossing time of the rising edge of the fundamental wave voltage of each phase according to the harmonic phasor of the voltage signal of each phase includes:

在每一当前时刻,判断每相基波电压相角是否符合以下条件:上一时刻的基波电压相角小于零,下一时刻的基波电压相角大于零;At each current moment, determine whether the phase angle of the fundamental wave voltage of each phase meets the following conditions: the phase angle of the fundamental wave voltage at the previous moment is less than zero, and the phase angle of the fundamental wave voltage at the next moment is greater than zero;

若满足,则当前时刻为过零点时刻。If satisfied, the current time is the zero-crossing time.

在其中一个实施例中,所述固定时长为5秒。In one of the embodiments, the fixed duration is 5 seconds.

在其中一个实施例中,所述各次谐波为3次谐波、5次谐波以及7次谐波。In one of the embodiments, the harmonics are the 3rd harmonic, the 5th harmonic and the 7th harmonic.

一种电网谐波相位分析系统,包括:A power grid harmonic phase analysis system, comprising:

信号获取模块,用于获取关注母线的电压波形信号和电力用户馈线上的电流波形信号;The signal acquisition module is used to acquire the voltage waveform signal of the concerned bus and the current waveform signal of the power user feeder;

波形截取模块,用于将统计时间划分为多个计算窗口;且在每个计算窗口内,对所述电压波形信号和电流波形信号截取预设数量个周波;a waveform interception module for dividing the statistical time into a plurality of calculation windows; and in each calculation window, intercepts a preset number of cycles from the voltage waveform signal and the current waveform signal;

谐波相量计算模块,用于对所述预设数量个周波进行快速傅里叶变换,得到每相电压信号的谐波相量和每相电流信号的谐波相量;a harmonic phasor calculation module, configured to perform fast Fourier transform on the preset number of cycles to obtain a harmonic phasor of each phase voltage signal and a harmonic phasor of each phase current signal;

谐波相位计算模块,用于根据所述每相电压信号的谐波相量获取每相基波电压上升沿的过零点时刻,并过零点时刻的相位作为参考基准,根据所述每相电流信号的谐波相量获取每相电流信号各次谐波在各相所述过零点时刻的相位;The harmonic phase calculation module is used to obtain the zero-crossing point moment of the rising edge of the fundamental wave voltage of each phase according to the harmonic phasor of the voltage signal of each phase, and use the phase of the zero-crossing point as a reference, according to the current signal of each phase The harmonic phasor of obtains the phase of each harmonic of each phase current signal at the zero-crossing time of each phase;

分析模块,用于基于在统计时间内获得的所有相位数据进行分析。Analysis module for analysis based on all phase data acquired in statistical time.

一种电网谐波相位分析装置,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现上述的方法的步骤。A power grid harmonic phase analysis device includes a memory and a processor, wherein the memory stores a computer program, and the processor implements the steps of the above method when executing the computer program.

上述电网谐波相位分析方法、系统和装置,通过以基波电压上升沿过零点为参考基准来得到各相各次谐波的相位,能够更准确地反映谐波电流相位分布。各相谐波电流相位都是对应于相应各相基波电压得到,例如根据基波电压A相上升沿过零点得到的是谐波电流A相的各次谐波相位,不会因为相序认定错误而受到影响,因此兼具有相序容错能力。The above power grid harmonic phase analysis method, system and device can more accurately reflect the harmonic current phase distribution by taking the zero-crossing point of the rising edge of the fundamental wave voltage as a reference to obtain the phase of each phase and each harmonic. The harmonic current phase of each phase is obtained corresponding to the fundamental wave voltage of each phase. For example, according to the zero-crossing point of the rising edge of phase A of the fundamental wave voltage, the phase of each harmonic current of phase A of the harmonic current is obtained, which will not be determined because of the phase sequence. It is affected by the error, so it has the capability of phase sequence fault tolerance.

附图说明Description of drawings

图1a为一实施例的电网谐波相位分析方法流程图;1a is a flowchart of a method for analyzing harmonic phase of a power grid according to an embodiment;

图1b为统计时间、计算窗口和截取的周波之间的关系;Figure 1b shows the relationship between statistical time, calculation window and intercepted cycles;

图2a为3次谐波电流相位分布图;Figure 2a is a phase distribution diagram of the third harmonic current;

图2b为5次谐波电流相位分布图;Figure 2b is the phase distribution diagram of the 5th harmonic current;

图2c为7次谐波电流相位分布图;Figure 2c is the phase distribution diagram of the 7th harmonic current;

图3为一实施例的电网谐波相位分析系统模块图。FIG. 3 is a block diagram of a power grid harmonic phase analysis system according to an embodiment.

具体实施方式Detailed ways

为了便于理解本发明,下面将参照相关附图对本发明进行更全面的描述。附图中给出了本发明的首选实施例。但是,本发明可以以许多不同的形式来实现,并不限于本文所描述的实施例。相反地,提供这些实施例的目的是使对本发明的公开内容更加透彻全面。In order to facilitate understanding of the present invention, the present invention will be described more fully hereinafter with reference to the related drawings. Preferred embodiments of the invention are shown in the accompanying drawings. However, the present invention may be embodied in many different forms and is not limited to the embodiments described herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete.

除非另有定义,本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在本发明的说明书中所使用的术语只是为了描述具体的实施例的目的,不是旨在于限制本发明。本文所使用的术语“及/或”包括一个或多个相关的所列项目的任意的和所有的组合。Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. The terms used herein in the description of the present invention are for the purpose of describing specific embodiments only, and are not intended to limit the present invention. As used herein, the term "and/or" includes any and all combinations of one or more of the associated listed items.

图1a为一实施例的电网谐波相位分析方法流程图。如图1a所示,所述的方法包括:FIG. 1a is a flowchart of a method for analyzing harmonic phase of a power grid according to an embodiment. As shown in Figure 1a, the method includes:

步骤S102:获取关注母线的电压波形信号和电力用户馈线上的电流波形信号。Step S102: Acquire the voltage waveform signal of the bus of interest and the current waveform signal on the power user feeder.

其中,关注母线是待分析谐波相位特性的母线,电力用户馈线为从关注母线分配出去的配电线路,与用户侧负荷相连。Among them, the busbar of interest is the busbar whose harmonic phase characteristics are to be analyzed, and the power user feeder is the distribution line allocated from the busbar of interest, and is connected to the load on the user side.

具体地,可以利用电压互感器(Potential Transformer,PT)、电流互感器(Current Transformer,CT)配合数据采集与监视控制系统(Supervisory Control AndData Acquisition System,SCADA)获取所述电压波形信号和电流波形信号。Specifically, the voltage waveform signal and the current waveform signal can be obtained by using a potential transformer (PT) and a current transformer (CT) in conjunction with a data acquisition and supervisory control system (Supervisory Control And Data Acquisition System, SCADA). .

SCADA是对分布距离远,生产单位分散的生产系统的一种数据采集、监视和控制系统。SCADA在远动系统中占重要地位,可以对现场的运行设备进行监视和控制,以实现数据采集、设备控制、测量、参数调节以及各类信号报警等各项功能。SCADA的重要组成部分包括远程终端单元(Remote Terminal Unit,RTU),和馈线终端单元(Feedback Terminal Unit,FTU)。SCADA is a data acquisition, monitoring and control system for the production system with long distribution distance and scattered production units. SCADA occupies an important position in the telecontrol system, which can monitor and control the running equipment in the field to realize various functions such as data acquisition, equipment control, measurement, parameter adjustment and various signal alarms. An important part of SCADA includes remote terminal unit (Remote Terminal Unit, RTU), and feeder terminal unit (Feedback Terminal Unit, FTU).

电压互感器和电流互感器是利用互感原理制作的电压和电流测量仪器。为方便数字化处理,采集的电压波形信号和电流波形信号为高频采样信号。在一个实施例中,采集关注母线的三相电压波形信号和电力用户馈线上的三相电流波形信号,所述电压波形信号和电流波形信号每相的采样频率不小于12.8kHz。Voltage transformers and current transformers are voltage and current measuring instruments made using the principle of mutual inductance. In order to facilitate digital processing, the collected voltage waveform signals and current waveform signals are high-frequency sampling signals. In one embodiment, the three-phase voltage waveform signal of the concerned bus and the three-phase current waveform signal on the power user feeder are collected, and the sampling frequency of each phase of the voltage waveform signal and the current waveform signal is not less than 12.8 kHz.

本实施例中,将电压波形信号记为uk(t),电流信号波形记为ik(t)。其中k代表配电网中的A、B、C三相,t代表时间。In this embodiment, the voltage waveform signal is denoted as uk (t), and the current signal waveform is denoted as ik (t). Where k represents the three phases of A, B, and C in the distribution network, and t represents time.

步骤S104:将统计时间划分为多个计算窗口;且在每个计算窗口内,对所述电压波形信号和电流波形信号截取预设数量个周波。Step S104: Divide the statistical time into a plurality of calculation windows; and in each calculation window, intercept a preset number of cycles from the voltage waveform signal and the current waveform signal.

为分析谐波相位特性,需要在一个较长的统计时间内采集信号波形进行分析和计算。在其中一个实施例中,所述计算窗口在统计时间内按固定时长等分。即按照固定时长进行一次截取和计算。例如统计时间为12个小时,若每隔5秒进行一次截取和计算,则需要进行8640次截取和计算。In order to analyze the harmonic phase characteristics, it is necessary to collect the signal waveform for analysis and calculation in a long statistical time. In one of the embodiments, the calculation window is divided equally into a fixed duration within the statistical time. That is, one interception and calculation is performed according to a fixed duration. For example, the statistical time is 12 hours. If the interception and calculation are performed every 5 seconds, 8640 interceptions and calculations are required.

计算窗口即计算时间段。本步骤在至少一个计算时间段内,对所述电压波形信号和电流波形信号截取预设数量个周波。周波是指一个信号周期的波形。例如对于工频电压而言,电压频率为50Hz,则一个信号周期为0.02秒。预设数量根据情况而定,例如根据计算能力、计算精度要求等确定截取的周波的数量。在其中一个实施例中,可以分别从所述电压波形信号和电流波形信号截取10个周波,即200毫秒时长的波形信号。统计时间、计算窗口和截取的周波之间的关系如图1b所示。The calculation window is the calculation time period. In this step, a preset number of cycles are intercepted from the voltage waveform signal and the current waveform signal during at least one calculation period. A cycle is the waveform of one signal period. For example, for the power frequency voltage, if the voltage frequency is 50Hz, one signal period is 0.02 seconds. The preset number is determined according to the situation, for example, the number of intercepted cycles is determined according to the calculation capability, calculation accuracy requirements, and the like. In one of the embodiments, 10 cycles, ie, waveform signals with a duration of 200 milliseconds, may be intercepted from the voltage waveform signal and the current waveform signal, respectively. The relationship between statistical time, calculation window and intercepted cycles is shown in Figure 1b.

步骤S106:对所述预设数量个周波进行快速傅里叶变换,得到每相电压信号的谐波相量和每相电流信号的谐波相量。Step S106: Perform fast Fourier transform on the preset number of cycles to obtain a harmonic phasor of each phase voltage signal and a harmonic phasor of each phase current signal.

以10个周波为例,在每个计算窗口内,将对10个周波信号进行快速傅里叶变换(Fast Fourier Transform,FFT)。即可分别得到一组电压谐波相量和一组电流谐波相量。Taking 10 cycles as an example, in each calculation window, fast Fourier transform (Fast Fourier Transform, FFT) will be performed on the 10 cycle signals. A set of voltage harmonic phasors and a set of current harmonic phasors can be obtained respectively.

在统计时间内,每次计算都能得到一组电压谐波相量和一组电流谐波相量。因此总共能得到N组电压谐波相量和N组电流谐波相量。本实施例中,将第n组电压谐波相量记为

Figure BDA0003489886950000061
将第n组电流谐波相量记为
Figure BDA0003489886950000062
其中n=1,2……N。k代表配电网中的A、B、C三相,h代表谐波次数,可以是3、5、7、……等。可以理解,
Figure BDA0003489886950000063
Figure BDA0003489886950000064
都是关于时间t的离散时间序列。其幅值相角的形式分别可表示为
Figure BDA0003489886950000071
Figure BDA0003489886950000072
In the statistical time, a set of voltage harmonic phasors and a set of current harmonic phasors can be obtained for each calculation. Therefore, a total of N groups of voltage harmonic phasors and N groups of current harmonic phasors can be obtained. In this embodiment, the nth group of voltage harmonic phasors are denoted as
Figure BDA0003489886950000061
Denote the nth group of current harmonic phasors as
Figure BDA0003489886950000062
where n=1,2...N. k represents the three phases of A, B, and C in the distribution network, and h represents the harmonic order, which can be 3, 5, 7, ... etc. understandably,
Figure BDA0003489886950000063
and
Figure BDA0003489886950000064
Both are discrete time series with respect to time t. The form of its amplitude and phase angle can be expressed as
Figure BDA0003489886950000071
and
Figure BDA0003489886950000072

步骤S108:根据所述每相电压信号的谐波相量获取每相基波电压上升沿的过零点时刻,并以过零点时刻的相位作为参考基准,根据所述每相电流信号的谐波相量获取每相电流信号各次谐波在各相所述过零点时刻的相位。Step S108: Obtain the zero-crossing time of the rising edge of the fundamental wave voltage of each phase according to the harmonic phasor of the voltage signal of each phase, and use the phase at the zero-crossing time as a reference, according to the harmonic phase of the current signal of each phase. Obtain the phase of each harmonic of each phase current signal at the zero-crossing point of each phase.

基波即频率与信号频率相同的波,为表示方便,用h=1表示基波。则每相基波的过零点时刻即各相基波电压相角为0的时刻,例如相基波电压相角为0可以表示为

Figure BDA0003489886950000073
The fundamental wave is the wave with the same frequency as the signal frequency. For convenience, h=1 is used to represent the fundamental wave. Then the zero-crossing moment of the fundamental wave of each phase is the moment when the phase angle of the fundamental wave voltage of each phase is 0. For example, when the phase angle of the fundamental wave voltage of each phase is 0, it can be expressed as
Figure BDA0003489886950000073

在其中一个实施例中,所述根据所述每相电压信号的谐波相量获取每相基波电压上升沿的过零点时刻,包括:In one embodiment, obtaining the zero-crossing time of the rising edge of the fundamental wave voltage of each phase according to the harmonic phasor of the voltage signal of each phase includes:

在每一当前时刻,判断每相基波电压相角是否符合以下条件:上一时刻的基波电压相角小于零,下一时刻的基波电压相角大于零;若满足,则当前时刻为过零点时刻。At each current moment, it is judged whether the phase angle of the fundamental wave voltage of each phase meets the following conditions: the phase angle of the fundamental wave voltage at the previous moment is less than zero, and the phase angle of the fundamental wave voltage at the next moment is greater than zero; if so, the current moment is zero o'clock time.

以A相为例,采用公式表示为:Taking phase A as an example, the formula is expressed as:

t0-1:

Figure BDA0003489886950000074
t 0 -1:
Figure BDA0003489886950000074

t0+1:

Figure BDA0003489886950000075
t 0 +1:
Figure BDA0003489886950000075

找到每相基波电压上升沿的过零点时刻t0之后,再找到同一组电流信号的谐波相量同相的各h次谐波相位,该各h次谐波相位为本申请所要求取的谐波相位。After finding the zero-crossing time t 0 of the rising edge of the fundamental wave voltage of each phase, find each h-order harmonic phase of the same phase of the harmonic phasors of the same group of current signals, and the h-order harmonic phase is the one required by the application. harmonic phase.

即:通过

Figure BDA0003489886950000076
获取A相基波电压相角
Figure BDA0003489886950000077
的时刻t0后,找到A相h次谐波电流
Figure BDA0003489886950000078
中时刻t0对应的相位
Figure BDA0003489886950000079
即为第n组相量数组中A相h次谐波电流以A相基波电压上升沿的过零点为参考的相位。That is: through
Figure BDA0003489886950000076
Get A-phase fundamental voltage phase angle
Figure BDA0003489886950000077
After the time t 0 of , find the h-th harmonic current of phase A
Figure BDA0003489886950000078
The phase corresponding to the middle time t 0
Figure BDA0003489886950000079
That is, the phase of the h-th harmonic current of phase A in the nth group of phasor arrays with the zero-crossing point of the rising edge of the fundamental voltage of phase A as the reference.

对B相和C相重复上述步骤,即可得到第n组相量数组中B、C相h次谐波电流相位

Figure BDA00034898869500000710
在其中一个实施例中,所述各次谐波为3次谐波、5次谐波以及7次谐波。Repeat the above steps for B-phase and C-phase to obtain the phase of the h-th harmonic current of B and C phases in the nth group of phasor arrays
Figure BDA00034898869500000710
In one of the embodiments, the harmonics are the 3rd harmonic, the 5th harmonic and the 7th harmonic.

步骤S110:基于在统计时间内获得的所有相位数据进行分析。Step S110: Perform analysis based on all phase data obtained within the statistical time.

对n=1,2……N中每组相量数组处理后,即可得到所有计算窗口内的谐波相位数据。利用所得到的所有谐波相位数据可进行电网谐波相位特性分析,从而为解决谐波污染提供依据,提升电能质量。在一个实施例中,分析谐波相位的方法是绘制谐波电流相位分布图。例如,对于3次谐波电流,将在每个计算窗口内获得的3次谐波电流相位数据表示在相位分布图中;对于5次谐波电流,将在每个计算窗口内获得的5次谐波电流相位数据表示在相位分布图中;等等以此类推。After processing each group of phasor arrays in n=1, 2...N, the harmonic phase data in all calculation windows can be obtained. All the obtained harmonic phase data can be used to analyze the harmonic phase characteristics of the power grid, thereby providing a basis for solving harmonic pollution and improving power quality. In one embodiment, the method of analyzing the harmonic phase is to plot the harmonic current phase distribution. For example, for the 3rd harmonic current, the phase data of the 3rd harmonic current obtained in each calculation window will be represented in the phase distribution diagram; for the 5th harmonic current, the 5th harmonic current obtained in each calculation window will be Harmonic current phase data is represented in a phase profile; and so on.

为了更好的说明,以下通过某公变居民用电负荷为例进行详细说明:For a better explanation, the following is a detailed description of the residential electricity load of a public substation as an example:

(1)利用电压互感器、电流互感器与SCADA系统获取该公变母线的电压波形信号uk(t)和电力用户馈线的电流波形信号ik(t),每个通道波形信号的采样频率为12.8kHz。(1) Use the voltage transformer, current transformer and SCADA system to obtain the voltage waveform signal u k (t) of the public transformer bus and the current waveform signal ik (t) of the power user feeder, and the sampling frequency of the waveform signal of each channel is 12.8kHz.

(2)总统计时间为12小时,快速傅里叶变换(FFT)计算间隔取5秒,则在统计时间内N为8640,即总共采集8640组电压波形信号和每相电流波形信号。在每组电压波形信号和每相电流波形信号中,对采集到的A相电压波形信号uA(t)和每相电流波形信号ik(t)分别截取10个周波,利用快速傅里叶变换(FFT)进行变换。得到的该公变母线上A相基波电压相量

Figure BDA0003489886950000081
以及每相电力用户馈线的3次、5次、7次谐波电流相量
Figure BDA0003489886950000082
其中k=A,B,C以及h=3,5,7。(2) The total statistical time is 12 hours, and the fast Fourier transform (FFT) calculation interval is 5 seconds, then N is 8640 in the statistical time, that is, a total of 8640 groups of voltage waveform signals and each phase current waveform signal are collected. In each group of voltage waveform signals and each phase current waveform signal, 10 cycles are intercepted from the collected A-phase voltage waveform signal u A (t) and each phase current waveform signal i k (t), respectively, using fast Fourier transform Transform (FFT) to transform. The obtained A-phase fundamental voltage phasor on the public transformer bus
Figure BDA0003489886950000081
And the 3rd, 5th and 7th harmonic current phasors of each phase power user feeder
Figure BDA0003489886950000082
where k=A,B,C and h=3,5,7.

(3)对于第n个相量数组,找到A相基波电压相角

Figure BDA0003489886950000083
的时刻t0作为A相基波电压上升沿的过零点,并以此为参考基准。(3) For the nth phasor array, find the phase angle of the fundamental voltage of phase A
Figure BDA0003489886950000083
The time t 0 is taken as the zero-crossing point of the rising edge of the fundamental voltage of phase A, and takes this as the reference.

(4)设k=A,在第n个相量数组中,找到A相h次谐波电流

Figure BDA0003489886950000084
中时刻t0对应的
Figure BDA0003489886950000085
即为第n组相量数组中A相h次谐波电流以A相基波电压上升沿的过零点为参考的相位。(4) Set k=A, in the nth phasor array, find the h-th harmonic current of phase A
Figure BDA0003489886950000084
Corresponding to the middle time t 0
Figure BDA0003489886950000085
That is, the phase of the h-th harmonic current of phase A in the nth group of phasor arrays with the zero-crossing point of the rising edge of the fundamental voltage of phase A as the reference.

(5)分别设k=B和k=C,重复步骤(4),得到第n个数组中B、C相h次谐波电流相位

Figure BDA0003489886950000091
(5) Set k=B and k=C respectively, and repeat step (4) to obtain the phase of the h-th harmonic current of B and C phases in the nth array
Figure BDA0003489886950000091

(6)重复步骤(3)、(4)、(5),直至n=8640。根据所得结果作出以A相基波电压过零点为基准,得到的3次、5次、7次谐波电流相位分布图。如图2a~图2c所示。(6) Repeat steps (3), (4), (5) until n=8640. According to the obtained results, the phase distribution diagrams of the 3rd, 5th, and 7th harmonic currents obtained based on the zero-crossing point of the A-phase fundamental wave voltage are made. As shown in Figure 2a ~ Figure 2c.

基于相同发明构思,还提供一种电网谐波相位分析系统。如图3所示,一实施例的电网谐波相位分析系统300包括:Based on the same inventive concept, a power grid harmonic phase analysis system is also provided. As shown in FIG. 3 , a power grid harmonic phase analysis system 300 according to an embodiment includes:

信号获取模块302,用于获取关注母线的电压波形信号和电力用户馈线上的电流波形信号。The signal acquisition module 302 is configured to acquire the voltage waveform signal of the bus of interest and the current waveform signal of the power user feeder.

其中,关注母线是待分析谐波相位特性的母线,电力用户馈线为从关注母线分配出去的配电线路,与用户侧负荷相连。Among them, the busbar of interest is the busbar whose harmonic phase characteristics are to be analyzed, and the power user feeder is the distribution line allocated from the busbar of interest, and is connected to the load on the user side.

具体地,可以利用电压互感器(Potential Transformer,PT)、电流互感器(Current Transformer,CT)配合数据采集与监视控制系统(Supervisory Control AndData Acquisition System,SCADA)获取所述电压波形信号和电流波形信号。Specifically, the voltage waveform signal and the current waveform signal can be obtained by using a potential transformer (PT) and a current transformer (CT) in conjunction with a data acquisition and supervisory control system (Supervisory Control And Data Acquisition System, SCADA). .

SCADA是对分布距离远,生产单位分散的生产系统的一种数据采集、监视和控制系统。SCADA在远动系统中占重要地位,可以对现场的运行设备进行监视和控制,以实现数据采集、设备控制、测量、参数调节以及各类信号报警等各项功能。SCADA的重要组成部分包括远程终端单元(Remote Terminal Unit,RTU),和馈线终端单元(Feedback Terminal Unit,FTU)。SCADA is a data acquisition, monitoring and control system for the production system with long distribution distance and scattered production units. SCADA occupies an important position in the telecontrol system, which can monitor and control the running equipment in the field to realize various functions such as data acquisition, equipment control, measurement, parameter adjustment and various signal alarms. An important part of SCADA includes remote terminal unit (Remote Terminal Unit, RTU), and feeder terminal unit (Feedback Terminal Unit, FTU).

电压互感器和电流互感器是利用互感原理制作的电压和电流测量仪器。为方便数字化处理,采集的电压波形信号和电流波形信号为高频采样信号。在一个实施例中,采集关注母线的三相电压波形信号和电力用户馈线上的三相电流波形信号,所述电压波形信号和电流波形信号每相的采样频率不小于12.8kHz。Voltage transformers and current transformers are voltage and current measuring instruments made using the principle of mutual inductance. In order to facilitate digital processing, the collected voltage waveform signals and current waveform signals are high-frequency sampling signals. In one embodiment, the three-phase voltage waveform signal of the concerned bus and the three-phase current waveform signal on the power user feeder are collected, and the sampling frequency of each phase of the voltage waveform signal and the current waveform signal is not less than 12.8 kHz.

本实施例中,将电压波形信号记为uk(t),电流信号波形记为ik(t)。其中k代表配电网中的A、B、C三相,t代表时间。In this embodiment, the voltage waveform signal is denoted as uk (t), and the current signal waveform is denoted as ik (t). Where k represents the three phases of A, B, and C in the distribution network, and t represents time.

波形截取模块304,用于将统计时间划分为多个计算窗口;且在每个计算窗口内,对所述电压波形信号和电流波形信号截取预设数量个周波。The waveform interception module 304 is configured to divide the statistical time into a plurality of calculation windows; and in each calculation window, intercept a preset number of cycles from the voltage waveform signal and the current waveform signal.

为分析谐波相位特性,需要在一个较长的统计时间内采集信号波形进行分析和计算。在其中一个实施例中,所述计算窗口在统计时间内按固定时长等分。即按照固定时长进行一次截取和计算。例如统计时间为12个小时,若每隔5秒进行一次截取和计算,则需要进行8640次截取和计算。In order to analyze the harmonic phase characteristics, it is necessary to collect the signal waveform for analysis and calculation in a long statistical time. In one of the embodiments, the calculation window is divided equally into a fixed duration within the statistical time. That is, one interception and calculation is performed according to a fixed duration. For example, the statistical time is 12 hours. If the interception and calculation are performed every 5 seconds, 8640 interceptions and calculations are required.

计算窗口即计算时间段。本步骤在至少一个计算时间段内,对所述电压波形信号和电流波形信号截取预设数量个周波。周波是指一个信号周期的波形。例如对于工频电压而言,电压频率为50Hz,则一个信号周期为0.02秒。预设数量根据情况而定,例如根据计算能力、计算精度要求等确定截取的周波的数量。在其中一个实施例中,可以分别从所述电压波形信号和电流波形信号截取10个周波,即200毫秒时长的波形信号。The calculation window is the calculation time period. In this step, a preset number of cycles are intercepted from the voltage waveform signal and the current waveform signal during at least one calculation period. A cycle is the waveform of one signal period. For example, for the power frequency voltage, if the voltage frequency is 50Hz, one signal period is 0.02 seconds. The preset number is determined according to the situation, for example, the number of intercepted cycles is determined according to the calculation capability, calculation accuracy requirements, and the like. In one of the embodiments, 10 cycles, ie, waveform signals with a duration of 200 milliseconds, may be intercepted from the voltage waveform signal and the current waveform signal, respectively.

谐波相量计算模块306,用于对所述预设数量个周波进行快速傅里叶变换,得到每相电压信号的谐波相量和每相电流信号的谐波相量。The harmonic phasor calculation module 306 is configured to perform fast Fourier transform on the preset number of cycles to obtain the harmonic phasor of each phase voltage signal and the harmonic phasor of each phase current signal.

以10个周波为例,在每个计算窗口内,将对10个周波信号进行快速傅里叶变换(Fast Fourier Transform,FFT)。即可分别得到一组电压谐波相量和一组电流谐波相量。Taking 10 cycles as an example, in each calculation window, fast Fourier transform (Fast Fourier Transform, FFT) will be performed on the 10 cycle signals. A set of voltage harmonic phasors and a set of current harmonic phasors can be obtained respectively.

在统计时间内,每次计算都能得到一组电压谐波相量和一组电流谐波相量。因此总共能得到N组电压谐波相量和N组电流谐波相量。本实施例中,将第n组电压谐波相量记为

Figure BDA0003489886950000101
将第n组电流谐波相量记为
Figure BDA0003489886950000102
其中n=1,2……N。k代表配电网中的A、B、C三相,h代表谐波次数,可以是3、5、7、……等。可以理解,
Figure BDA0003489886950000111
Figure BDA0003489886950000112
都是关于时间t的离散时间序列。其幅值相角的形式分别可表示为
Figure BDA0003489886950000113
Figure BDA0003489886950000114
In the statistical time, a set of voltage harmonic phasors and a set of current harmonic phasors can be obtained for each calculation. Therefore, N groups of voltage harmonic phasors and N groups of current harmonic phasors can be obtained in total. In this embodiment, the nth group of voltage harmonic phasors are denoted as
Figure BDA0003489886950000101
Denote the nth group of current harmonic phasors as
Figure BDA0003489886950000102
where n=1,2...N. k represents the three phases of A, B, and C in the distribution network, and h represents the harmonic order, which can be 3, 5, 7, ... etc. understandably,
Figure BDA0003489886950000111
and
Figure BDA0003489886950000112
are all discrete time series with respect to time t. The form of its amplitude and phase angle can be expressed as
Figure BDA0003489886950000113
and
Figure BDA0003489886950000114

谐波相位计算模块308,用于根据所述每相电压信号的谐波相量获取每相基波电压上升沿的过零点时刻,并以过零点时刻的相位作为参考基准,根据所述每相电流信号的谐波相量获取每相电流信号各次谐波在各相所述过零点时刻的相位。The harmonic phase calculation module 308 is configured to obtain the zero-crossing point moment of the rising edge of the fundamental wave voltage of each phase according to the harmonic phasor of the voltage signal of each phase, and use the phase of the zero-crossing point as a reference. The harmonic phasor of the current signal obtains the phase of each harmonic of each phase current signal at the zero-crossing point of each phase.

基波即频率与信号频率相同的波,为表示方便,用h=1表示基波。则每相基波的过零点时刻即各相基波电压相角为0的时刻,例如相基波电压相角为0可以表示为

Figure BDA0003489886950000115
The fundamental wave is the wave with the same frequency as the signal frequency. For convenience, h=1 is used to represent the fundamental wave. Then the zero-crossing moment of the fundamental wave of each phase is the moment when the phase angle of the fundamental wave voltage of each phase is 0. For example, when the phase angle of the fundamental wave voltage of each phase is 0, it can be expressed as
Figure BDA0003489886950000115

在其中一个实施例中,所述根据所述每相电压信号的谐波相量获取每相基波电压上升沿的过零点时刻,包括:In one embodiment, obtaining the zero-crossing time of the rising edge of the fundamental wave voltage of each phase according to the harmonic phasor of the voltage signal of each phase includes:

在每一当前时刻,判断每相基波电压相角是否符合以下条件:上一时刻的基波电压相角小于零,下一时刻的基波电压相角大于零;若满足,则当前时刻为过零点时刻。At each current moment, it is judged whether the phase angle of the fundamental wave voltage of each phase meets the following conditions: the phase angle of the fundamental wave voltage at the previous moment is less than zero, and the phase angle of the fundamental wave voltage at the next moment is greater than zero; if so, the current moment is zero o'clock time.

以A相为例,采用公式表示为:Taking phase A as an example, the formula is expressed as:

t0-1:

Figure BDA0003489886950000116
t 0 -1:
Figure BDA0003489886950000116

t0+1:

Figure BDA0003489886950000117
t 0 +1:
Figure BDA0003489886950000117

找到每相基波电压上升沿的过零点时刻t0之后,再找到同一组电流信号的谐波相量同相的各h次谐波相位,该各h次谐波相位为本申请所要求取的谐波相位。After finding the zero-crossing time t 0 of the rising edge of the fundamental wave voltage of each phase, find each h-order harmonic phase of the same phase of the harmonic phasors of the same group of current signals, and the h-order harmonic phase is the one required by the application. harmonic phase.

即:通过

Figure BDA0003489886950000118
获取A相基波电压相角
Figure BDA0003489886950000119
的时刻t0后,找到A相h次谐波电流
Figure BDA00034898869500001110
中时刻t0对应的相位
Figure BDA00034898869500001111
即为第n组相量数组中A相h次谐波电流以A相基波电压上升沿的过零点为参考的相位。That is: through
Figure BDA0003489886950000118
Get A-phase fundamental voltage phase angle
Figure BDA0003489886950000119
After the time t 0 of , find the h-th harmonic current of phase A
Figure BDA00034898869500001110
The phase corresponding to the middle time t 0
Figure BDA00034898869500001111
That is, the phase of the h-th harmonic current of phase A in the nth group of phasor arrays with the zero-crossing point of the rising edge of the fundamental voltage of phase A as the reference.

对B相和C相重复上述步骤,即可得到第n组相量数组中B、C相h次谐波电流相位

Figure BDA0003489886950000121
在其中一个实施例中,所述各次谐波为3次谐波、5次谐波以及7次谐波。Repeat the above steps for B-phase and C-phase to obtain the phase of the h-th harmonic current of B and C phases in the nth group of phasor arrays
Figure BDA0003489886950000121
In one of the embodiments, the harmonics are the 3rd harmonic, the 5th harmonic and the 7th harmonic.

分析模块310,用于基于在统计时间内获得的所有相位数据进行分析。An analysis module 310 for performing analysis based on all phase data obtained within the statistical time.

对n=1,2……N中每组相量数组处理后,即可得到所有计算窗口内的谐波相位数据。利用所得到的所有谐波相位数据可进行电网谐波相位特性分析,从而为解决谐波污染提供依据,提升电能质量。在一个实施例中,分析谐波相位的方法是绘制谐波电流相位分布图。例如,对于3次谐波电流,将在每个计算窗口内获得的3次谐波电流相位数据表示在相位分布图中;对于5次谐波电流,将在每个计算窗口内获得的5次谐波电流相位数据表示在相位分布图中;等等以此类推。After processing each group of phasor arrays in n=1, 2...N, the harmonic phase data in all calculation windows can be obtained. All the obtained harmonic phase data can be used to analyze the harmonic phase characteristics of the power grid, thereby providing a basis for solving harmonic pollution and improving power quality. In one embodiment, the method of analyzing the harmonic phase is to plot the harmonic current phase distribution. For example, for the 3rd harmonic current, the phase data of the 3rd harmonic current obtained in each calculation window will be represented in the phase distribution diagram; for the 5th harmonic current, the 5th harmonic current obtained in each calculation window will be Harmonic current phase data is represented in a phase profile; and so on.

基于相同发明构思,还提供一种电网谐波相位分析装置。该电网谐波相位分析装置包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现上述的方法的步骤。Based on the same inventive concept, a power grid harmonic phase analysis device is also provided. The power grid harmonic phase analysis device includes a memory and a processor, wherein the memory stores a computer program, and the processor implements the steps of the above method when executing the computer program.

上述电网谐波相位分析方法、系统和装置,通过以基波电压上升沿过零点为参考基准来得到各相各次谐波的相位,能够更准确地反映谐波电流相位分布。各相谐波电流相位都是对应于相应各相基波电压得到,例如根据基波电压A相上升沿过零点得到的是谐波电流A相的各次谐波相位,不会因为相序认定错误而受到影响,因此兼具有相序容错能力。The above power grid harmonic phase analysis method, system and device can more accurately reflect the harmonic current phase distribution by taking the zero-crossing point of the rising edge of the fundamental wave voltage as a reference to obtain the phase of each phase and each harmonic. The harmonic current phase of each phase is obtained corresponding to the fundamental wave voltage of each phase. For example, according to the zero-crossing point of the rising edge of phase A of the fundamental wave voltage, the phase of each harmonic current of phase A of the harmonic current is obtained, which will not be determined because of the phase sequence. It is affected by the error, so it has the capability of phase sequence fault tolerance.

以上所述实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。The technical features of the above-described embodiments can be combined arbitrarily. For the sake of brevity, all possible combinations of the technical features in the above-described embodiments are not described. However, as long as there is no contradiction between the combinations of these technical features, All should be regarded as the scope described in this specification.

以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。The above-mentioned embodiments only represent several embodiments of the present invention, and the descriptions thereof are specific and detailed, but should not be construed as a limitation on the scope of the invention patent. It should be pointed out that for those of ordinary skill in the art, without departing from the concept of the present invention, several modifications and improvements can also be made, which all belong to the protection scope of the present invention. Therefore, the protection scope of the patent of the present invention should be subject to the appended claims.

Claims (10)

1.一种电网谐波相位分析方法,其特征在于,所述的方法包括:1. a power grid harmonic phase analysis method, is characterized in that, described method comprises: 获取关注母线的电压波形信号和电力用户馈线上的电流波形信号;Obtain the voltage waveform signal of the concerned bus and the current waveform signal of the power user feeder; 将统计时间划分为多个计算窗口;且在每个计算窗口内,对所述电压波形信号和电流波形信号截取预设数量个周波;dividing the statistical time into a plurality of calculation windows; and in each calculation window, intercepting a preset number of cycles from the voltage waveform signal and the current waveform signal; 对所述预设数量个周波进行快速傅里叶变换,得到每相电压信号的谐波相量和每相电流信号的谐波相量;performing fast Fourier transform on the preset number of cycles to obtain the harmonic phasor of each phase voltage signal and the harmonic phasor of each phase current signal; 根据所述每相电压信号的谐波相量获取每相基波电压上升沿的过零点时刻,并以过零点时刻的相位作为参考基准,根据所述每相电流信号的谐波相量获取每相电流信号各次谐波在各相所述过零点时刻的相位;Obtain the zero-crossing time of the rising edge of the fundamental voltage of each phase according to the harmonic phasor of the voltage signal of each phase, and take the phase at the zero-crossing time as a reference, and obtain each phase according to the harmonic phasor of the current signal of each phase. The phase of each harmonic of the phase current signal at the zero-crossing time of each phase; 基于在统计时间内获得的所有相位数据进行分析。Analysis is based on all phase data acquired within statistical time. 2.根据权利要求1所述的电网谐波相位分析方法,其特征在于,所述计算窗口在统计时间内按固定时长等分。2 . The power grid harmonic phase analysis method according to claim 1 , wherein the calculation window is divided into equal parts by a fixed duration in the statistical time. 3 . 3.根据权利要求1所述的电网谐波相位分析方法,其特征在于,所述获取关注母线的电压波形信号和电力用户馈线上的电流波形信号,包括:3. The power grid harmonic phase analysis method according to claim 1, wherein the acquiring the voltage waveform signal of the concerned bus and the current waveform signal on the power user feeder comprises: 利用电压互感器、电流互感器配合数据采集与监视控制系统获取所述电压波形信号和电流波形信号。The voltage waveform signal and the current waveform signal are obtained by using the voltage transformer and the current transformer in cooperation with the data acquisition and monitoring control system. 4.根据权利要求1所述的电网谐波相位分析方法,其特征在于,所述电压波形信号和电流波形信号每个通道的采样频率不小于12.8kHz。4 . The power grid harmonic phase analysis method according to claim 1 , wherein the sampling frequency of each channel of the voltage waveform signal and the current waveform signal is not less than 12.8 kHz. 5 . 5.根据权利要求1所述的电网谐波相位分析方法,其特征在于,所述预设数量为10。5 . The power grid harmonic phase analysis method according to claim 1 , wherein the preset number is 10. 6 . 6.根据权利要求1所述的电网谐波相位分析方法,其特征在于,所述根据所述每相电压信号的谐波相量获取每相基波电压上升沿的过零点时刻,包括:6 . The power grid harmonic phase analysis method according to claim 1 , wherein, obtaining the zero-crossing time of the rising edge of the fundamental voltage of each phase according to the harmonic phasor of the voltage signal of each phase, comprising: 6 . 在每一当前时刻,判断每相基波电压相角是否符合以下条件:上一时刻的基波电压相角小于零,下一时刻的基波电压相角大于零;At each current moment, it is judged whether the phase angle of the fundamental wave voltage of each phase meets the following conditions: the phase angle of the fundamental wave voltage at the previous moment is less than zero, and the phase angle of the fundamental wave voltage at the next moment is greater than zero; 若满足,则当前时刻为过零点时刻。If satisfied, the current time is the zero-crossing time. 7.根据权利要求2所述的电网谐波相位分析方法,其特征在于,所述固定时长为5秒。7 . The power grid harmonic phase analysis method according to claim 2 , wherein the fixed duration is 5 seconds. 8 . 8.根据权利要求1所述的电网谐波相位分析方法,其特征在于,所述各次谐波为3次谐波、5次谐波以及7次谐波。8 . The power grid harmonic phase analysis method according to claim 1 , wherein the harmonics are the 3rd harmonic, the 5th harmonic and the 7th harmonic. 9 . 9.一种电网谐波相位分析系统,包括:9. A power grid harmonic phase analysis system, comprising: 信号获取模块,用于获取关注母线的电压波形信号和电力用户馈线上的电流波形信号;The signal acquisition module is used to acquire the voltage waveform signal of the concerned bus and the current waveform signal of the power user feeder; 波形截取模块,用于将统计时间划分为多个计算窗口;且在每个计算窗口内,对所述电压波形信号和电流波形信号截取预设数量个周波;a waveform interception module, configured to divide the statistical time into a plurality of calculation windows; and in each calculation window, intercept a preset number of cycles from the voltage waveform signal and the current waveform signal; 谐波相量计算模块,用于对所述预设数量个周波进行快速傅里叶变换,得到每相电压信号的谐波相量和每相电流信号的谐波相量;a harmonic phasor calculation module, configured to perform fast Fourier transform on the preset number of cycles to obtain the harmonic phasor of each phase voltage signal and the harmonic phasor of each phase current signal; 谐波相位计算模块,用于根据所述每相电压信号的谐波相量获取每相基波电压上升沿的过零点时刻,并以过零点时刻的相位作为参考基准,根据所述每相电流信号的谐波相量获取每相电流信号各次谐波在各相所述过零点时刻的相位;The harmonic phase calculation module is used to obtain the zero-crossing moment of the rising edge of the fundamental voltage of each phase according to the harmonic phasor of the voltage signal of each phase, and use the phase at the zero-crossing moment as a reference, according to the current of each phase The harmonic phasor of the signal obtains the phase of each harmonic of each phase current signal at the zero-crossing time of each phase; 分析模块,用于基于在统计时间内获得的所有相位数据进行分析。Analysis module for analysis based on all phase data acquired in statistical time. 10.一种电网谐波相位分析装置,其特征在于,包括存储器和处理器,所述存储器存储有计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至8中任一项所述的方法的步骤。10. A power grid harmonic phase analysis device, characterized in that it comprises a memory and a processor, wherein the memory stores a computer program, characterized in that, when the processor executes the computer program, the implementation of claims 1 to 8 The steps of any one of the methods.
CN202210093206.8A 2022-01-26 2022-01-26 Power grid harmonic phase analysis method, system and device Active CN114460366B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210093206.8A CN114460366B (en) 2022-01-26 2022-01-26 Power grid harmonic phase analysis method, system and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210093206.8A CN114460366B (en) 2022-01-26 2022-01-26 Power grid harmonic phase analysis method, system and device

Publications (2)

Publication Number Publication Date
CN114460366A true CN114460366A (en) 2022-05-10
CN114460366B CN114460366B (en) 2025-05-23

Family

ID=81411186

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210093206.8A Active CN114460366B (en) 2022-01-26 2022-01-26 Power grid harmonic phase analysis method, system and device

Country Status (1)

Country Link
CN (1) CN114460366B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114935687A (en) * 2022-05-11 2022-08-23 深圳供电局有限公司 Load harmonic emission characteristic extraction method, device and computer equipment
CN115469150A (en) * 2022-08-19 2022-12-13 国网浙江省电力有限公司电力科学研究院 Harmonic wave tracing method based on fault recording
CN117148260A (en) * 2023-03-02 2023-12-01 国网河北省电力有限公司营销服务中心 Measurement deviation positioning method, device, terminal and storage medium
CN118884135A (en) * 2024-09-29 2024-11-01 国网浙江省电力有限公司温州供电公司 A novel method, system, device and medium for real-time positioning of harmonics in power system

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1193102A (en) * 1966-06-14 1970-05-28 Hoffman Electronics Corp Phase Angle Measurement System.
US4667152A (en) * 1985-10-17 1987-05-19 Hayes Raymond M Method of and system for determining locations of sources of harmonics in a power distribution network
JP2005061939A (en) * 2003-08-11 2005-03-10 Yokogawa Electric Corp Harmonic analyzer
KR20060038527A (en) * 2004-10-30 2006-05-04 한국전력공사 Distribution automation terminal device and method with real-time monitoring of electrical quality of distribution system
JP2011229361A (en) * 2010-03-31 2011-11-10 Daihen Corp Phase detection device
KR20160064695A (en) * 2014-11-28 2016-06-08 서울과학기술대학교 산학협력단 Apparatus and method for detecting phase and harmonic using fast fourier transform
CN107543962A (en) * 2017-08-11 2018-01-05 安徽大学 The computational methods of leading m-Acetyl chlorophosphonazo spectrum distribution
CN111400649A (en) * 2020-03-27 2020-07-10 深圳供电局有限公司 Harmonic responsibility quantification method and device, computer equipment and storage medium
CN115144653A (en) * 2022-05-26 2022-10-04 国网北京市电力公司 Harmonic phase angle calculation method, device, storage medium and equipment

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1193102A (en) * 1966-06-14 1970-05-28 Hoffman Electronics Corp Phase Angle Measurement System.
US4667152A (en) * 1985-10-17 1987-05-19 Hayes Raymond M Method of and system for determining locations of sources of harmonics in a power distribution network
JP2005061939A (en) * 2003-08-11 2005-03-10 Yokogawa Electric Corp Harmonic analyzer
KR20060038527A (en) * 2004-10-30 2006-05-04 한국전력공사 Distribution automation terminal device and method with real-time monitoring of electrical quality of distribution system
JP2011229361A (en) * 2010-03-31 2011-11-10 Daihen Corp Phase detection device
KR20160064695A (en) * 2014-11-28 2016-06-08 서울과학기술대학교 산학협력단 Apparatus and method for detecting phase and harmonic using fast fourier transform
CN107543962A (en) * 2017-08-11 2018-01-05 安徽大学 The computational methods of leading m-Acetyl chlorophosphonazo spectrum distribution
CN111400649A (en) * 2020-03-27 2020-07-10 深圳供电局有限公司 Harmonic responsibility quantification method and device, computer equipment and storage medium
CN115144653A (en) * 2022-05-26 2022-10-04 国网北京市电力公司 Harmonic phase angle calculation method, device, storage medium and equipment

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
沈杰;孟晓芳;王英男;侯森;张帆;赵夕彤;: "配电网谐波分析的傅里叶变换方法", 东北电力技术, no. 10, 20 October 2015 (2015-10-20) *
罗楚军;袁召;尹小根;何俊佳;陈轩恕;潘垣;: "基于改进半波傅氏算法的故障电流相控开断零点预测研究", 高压电器, no. 06, 16 June 2012 (2012-06-16) *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114935687A (en) * 2022-05-11 2022-08-23 深圳供电局有限公司 Load harmonic emission characteristic extraction method, device and computer equipment
CN115469150A (en) * 2022-08-19 2022-12-13 国网浙江省电力有限公司电力科学研究院 Harmonic wave tracing method based on fault recording
CN117148260A (en) * 2023-03-02 2023-12-01 国网河北省电力有限公司营销服务中心 Measurement deviation positioning method, device, terminal and storage medium
CN117148260B (en) * 2023-03-02 2024-05-28 国网河北省电力有限公司营销服务中心 Method, device, terminal and storage medium for locating measurement deviation
CN118884135A (en) * 2024-09-29 2024-11-01 国网浙江省电力有限公司温州供电公司 A novel method, system, device and medium for real-time positioning of harmonics in power system
CN118884135B (en) * 2024-09-29 2024-12-13 国网浙江省电力有限公司温州供电公司 Novel power system harmonic wave real-time positioning method, system, equipment and medium

Also Published As

Publication number Publication date
CN114460366B (en) 2025-05-23

Similar Documents

Publication Publication Date Title
CN114460366A (en) Power grid harmonic phase analysis method, system and device
CN102967800B (en) Method and device for positioning single-phase ground fault section of power distribution network based on transient signal prony algorithm
CN104142422A (en) Transformation substation leakage current and capacitive current online monitoring management system and working method thereof
CN107219391A (en) Quality of voltage Monitoring Indexes device
Guo et al. Global-sensitivity-based theoretical analysis and fast prediction of traveling waves with respect to fault resistance on HVDC transmission lines
CN102183685A (en) Flexible wide area power grid phasor measurement method
Oleinikova et al. PMU measurements application for transmission line temperature and sag estimation algorithm development
Oliveira et al. A contribution for the measuring process of harmonics and interharmonics in electrical power systems with photovoltaic sources
Calin et al. Measurement data aggregation for active distribution networks
CN203204090U (en) Substation zinc oxide arrester on-line monitoring device
Zhurkin et al. Application of the short time Fourier transform for calculating electric energy losses under the conditions of a non-sinusoidal voltage
Seger et al. Power system monitoring through low-voltage distribution network using freePMU
Chen et al. Locating sub-cycle faults in distribution network applying half-cycle DFT method
CN108982954B (en) Method and system for calculating phase voltage amplitude and phase suitable for feeder line terminal
CN105637730A (en) Power system control
CN102004854A (en) Method for acquiring nonlinear load equivalent impedance parameters of distribution network in real time
Xu et al. Operating cycle performance, lost periodicity, and waveform distortion of switch-mode power supplies
CN117665691A (en) Remote online verification device for electric energy meters based on Internet of Things technology
Farantatos et al. Implementation of a 3-phase state estimation tool suitable for advanced distribution management systems
An et al. Supraharmonics Measurement Based on Hybrid Online Measurement and Offline Analysis
CN106771678A (en) A kind of phase detecting method and device based on Hilbert-Huang transform and expert system
CN114460359A (en) Method, device, equipment and storage medium for eliminating harmonic voltage measurement error of power grid
CN112713589A (en) Low-voltage distribution network line impedance estimation method based on measurement of power quality analyzer
Balabin et al. Evaluation of PMU performance during transients
Vats et al. Application of phasor & frequency estimation techniques in phasor measurement unit

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