阅读说明：本技术 一种滤波器、同步相量测量系统及方法 (Filter, synchronous phasor measurement system and method ) 是由 薛蕙 阮梦洁 程一帆 于 2019-08-08 设计创作，主要内容包括：本发明实施例提供一种滤波器、同步相量测量系统及方法。所述滤波器为固定窗口长度的自适应滑动平均滤波器；所述系统包括滤波器和锁相环控制环路,所述锁相环控制环路的输出端与所述滤波器的输入端进行级联；所述方法包括获取智能电网输出的电压信号；将所述电压信号输入至所述同步相量测量系统,获取所述同步相量测量系统输出的相量测量结果和频率测量结果。本发明实施例通过基于锁相环的同步相量测量算法,利用固定长度自适应滑动平均滤波器对1型锁相环输出的波纹进行滤波,并给出估计参数的时标,最后的测量结果表明基于同步相量测量系统的方法能够在电网畸变的情况下准确地跟踪电网频率和相量,其结果的响应时间和测量精度能满足IEEE C37.118的要求。(The embodiment of the invention provides a filter, a synchronous phasor measurement system and a synchronous phasor measurement method. The filter is an adaptive sliding average filter with a fixed window length; the system comprises a filter and a phase-locked loop control loop, wherein the output end of the phase-locked loop control loop is cascaded with the input end of the filter; the method comprises the steps of obtaining a voltage signal output by the smart grid; and inputting the voltage signal into the synchronous phasor measurement system, and acquiring a phasor measurement result and a frequency measurement result output by the synchronous phasor measurement system. According to the embodiment of the invention, through a synchronous phasor measurement algorithm based on the phase-locked loop, the fixed-length self-adaptive sliding average filter is used for filtering the ripple output by the type 1 phase-locked loop, and the time scale of an estimated parameter is given, and the final measurement result shows that the method based on the synchronous phasor measurement system can accurately track the frequency and the phasor of a power grid under the condition of power grid distortion, and the response time and the measurement precision of the result can meet the requirements of IEEE C37.118.)

1. A filter is characterized in that the filter is an adaptive sliding average filter with a fixed window length, and the value of the fixed window length is obtained based on the number of sampling points and the sampling frequency.

2. A synchronized phasor measurement system comprising a filter according to claim 1, and further comprising a phase locked loop control loop, the output of which is cascaded with the input of the filter.

3. a synchronized phasor measurement system according to claim 2, wherein said phase locked loop control loop includes a filter combination consisting of a number of non-adaptive sliding average filters connected in series, a first park transformation system connected in cascade to an input of said filter combination, a second park transformation system connected in cascade to an output of said filter combination, and a type 1 phase locked loop connected in cascade to an output of said second park transformation system;

wherein the filter combination is used for pre-filtering.

4. a synchronized phasor measurement system according to claim 2, wherein said output of said phase locked loop control loop includes a phase angle component output, a magnitude component output and a frequency component output;

Wherein the phase angle component output is cascaded with two of the filters in series;

The amplitude component output terminal and the frequency component output terminal are respectively cascaded with one of the filters.

5. A smart grid synchrophasor measurement method based on the synchrophasor measurement system according to any one of claims 2 to 4, comprising:

Acquiring a voltage signal output by the smart grid;

And inputting the voltage signal into the synchronous phasor measurement system, and acquiring a phasor measurement result and a frequency measurement result output by the synchronous phasor measurement system.

6. the smart grid synchronous phasor measurement method according to claim 5, wherein the smart grid is a three-phase system or a single-phase system.

7. The smart grid synchronized phasor measurement method according to claim 5, wherein the inputting the voltage signal into the synchronized phasor measurement system, obtaining phasor measurement results and frequency measurement results output by the synchronized phasor measurement system, further comprises:

Setting time scales of the phasor measurement result and the frequency measurement result;

The time stamp is obtained by measuring delays of the phasor measurement and the frequency measurement.

8. The smart grid synchronous phasor measurement method according to claim 7, wherein the measured delay is a measured delay obtained based on phase angle, amplitude, frequency and frequency change rate.

9. The smart grid synchronized phasor measurement method according to any one of claims 5 to 8, further comprising:

Resetting a rotation angle of a second park transformation system in the synchrophasor measurement system after a preset time.

10. The smart grid synchrophasor measurement method according to claim 5, further comprising: and respectively carrying out performance evaluation on the response time and the measurement precision of the phasor measurement result and the frequency measurement result.

Technical Field

The invention relates to the technical field of smart power grids, in particular to a filter, a synchronous phasor measurement system and a synchronous phasor measurement method.

Background

A synchronous Phasor Measurement Unit (PMU) plays an important role in monitoring power grid dynamics in a wide-area and real-time manner. Wherein, the synchrophasor measurement algorithm is the core of PMU. Therefore, an accurate and fast synchronous phasor measurement algorithm is a key to the wide application of the PMU in the full-system monitoring and reliable operation of the smart grid.

Disclosure of Invention

The embodiment of the invention provides a filter, a synchronized phasor measurement system and a synchronized phasor measurement method, which are used for solving the problems of low accuracy of a measurement result and slow response time caused by the fact that the frequency deviation is greatly influenced by power grid distortion when the state of an intelligent power grid is monitored based on the synchronized phasor measurement in the prior art.

In a first aspect, an embodiment of the present invention provides a filter, including:

The filter is an adaptive sliding average filter with a fixed window length, and the value of the fixed window length is obtained based on the number of sampling points and the sampling frequency.

In a second aspect, an embodiment of the present invention provides a synchronized phasor measurement system, including:

the filter also comprises a phase-locked loop control loop, and the output end of the phase-locked loop control loop is cascaded with the input end of the filter.

Preferably, the phase-locked loop control loop comprises a filter combination formed by connecting a plurality of non-adaptive sliding average filters in series, a first park transformation system cascaded with an input end of the filter combination, a second park transformation system cascaded with an output end of the filter combination, and a type-1 phase-locked loop cascaded with an output end of the second park transformation system;

Wherein the filter combination is used for pre-filtering.

Preferably, the output end of the phase-locked loop control loop comprises a phase angle component output end, an amplitude component output end and a frequency component output end;

Wherein the phase angle component output is cascaded with two of the filters in series;

The amplitude component output terminal and the frequency component output terminal are respectively cascaded with one of the filters.

In a third aspect, an embodiment of the present invention provides a smart grid synchronous phasor measurement method, including:

Acquiring a voltage signal output by the smart grid;

And inputting the voltage signal into the synchronous phasor measurement system, and acquiring a phasor measurement result and a frequency measurement result output by the synchronous phasor measurement system.

Preferably, the smart grid is a three-phase system or a single-phase system.

Preferably, the inputting the voltage signal into the synchronized phasor measurement system to obtain the phasor measurement result and the frequency measurement result output by the synchronized phasor measurement system further includes:

Setting time scales of the phasor measurement result and the frequency measurement result;

The time stamp is obtained by measuring delays of the phasor measurement and the frequency measurement.

Preferably, the measured delay is a measured delay based on phase angle, amplitude, frequency and frequency rate of change.

Preferably, the method further comprises:

Resetting a rotation angle of a second park transformation system in the synchrophasor measurement system after a preset time.

Preferably, the method further comprises: and respectively carrying out performance evaluation on the response time and the measurement precision of the phasor measurement result and the frequency measurement result.

According to the filter, the synchronous phasor measurement system and the method provided by the embodiment of the invention, the fixed-length self-adaptive sliding average filter is used for filtering the ripple output by the type 1 phase-locked loop through the synchronous phasor measurement algorithm based on the phase-locked loop, the time scale of an estimated parameter is given, and the final result verifies that the frequency and the phasor of a power grid can be accurately tracked under the condition of power grid distortion, and the response time and the measurement precision of the filter can meet the requirements of IEEE C37.118.

Drawings

In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.

FIG. 1 is a graph comparing frequency responses of various filters provided by embodiments of the present invention;

FIG. 2 is a graph comparing the filtering results of various filters provided by embodiments of the present invention;

Fig. 3 is a structural diagram of a synchronized phasor measurement system according to an embodiment of the present invention;

Fig. 4 is an internal structural diagram of a phase-locked loop control loop according to an embodiment of the present invention;

FIG. 5 is a block diagram of a cascade of a plurality of non-adaptive sliding average filters according to an embodiment of the present invention;

fig. 6 is a flowchart of a synchronized phasor measurement method for a smart grid according to an embodiment of the present invention;

FIG. 7 is a comparison graph of the total error of vector estimation under step change of amplitude, phase angle and frequency according to an embodiment of the present invention;

FIG. 8 is a graph illustrating frequency error versus amplitude, phase angle, and frequency step change provided by an embodiment of the present invention;

Fig. 9 is a comparison graph of frequency change rate error under amplitude, phase angle, and frequency step change provided by the embodiment of the present invention.

Detailed Description

In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be obtained by a person skilled in the art without any inventive step based on the embodiments of the present invention, are within the scope of the present invention.

In order to solve the problems in the prior art, an embodiment of the present invention provides a method for measuring a synchronous phasor based on a phase-locked loop, where the method is based on a synchronous phasor measurement system, and the system includes a phase-locked loop and an Adaptive sliding Average Filter (a Fixed length Adaptive Moving Average Filter, hereinafter referred to as FAMAF) with a Fixed window length.

An embodiment of the present invention provides a filter, including:

The filter is an adaptive sliding average filter with a fixed window length, and the value of the fixed window length is obtained based on the number of sampling points and the sampling frequency.

In particular, the MAF in the embodiment of the present invention is obtained based on a Moving Average Filter (MAF), which is widely used due to its simple, effective, and easy implementation characteristics. The course of a MAF can be expressed as:

Wherein N is_w＝T_w/T_s，T_wis the window length, T, of the MAF_s＝1/f_s，f_sIs the sampling frequency, h_MAF(m) is the impulse response function of MAF, h_MAF(m)＝[1,1,...,1]，m＝0,1,...,N_w-1. Equation (1) is equivalent to a pulse response of h_MAF(m) a Finite Impulse Response (FIR) filtering process. The frequency response of a MAF can be described as:

Wherein f is_w＝1/T_w. Can prove H_MAF(kf_w) 0, it means that MAF can be filtered out at f_wAt integer multiples of (a).

If T in the formulas (1) and (2)_wnot adjusted with the input signal, thenMAF is a non-adaptive MAF (NMAF). The NMAF has fixed filtering capability, and the filtering capability is sensitive to the frequency change of a power grid. For example, the window length of the NMAF is set as the nominal period of the power grid, so that harmonics can be better filtered; however, the filtering capability of the NMAF decreases when the frequency shifts.

on the contrary, if T_wThis MAF is an Adaptive MAF (AMAF) that can be adjusted with the input signal. AMAFs have a variable frequency characteristic from which interference at frequency offsets can be filtered out. However, the window length of the AMAF is variable, which makes it difficult to accurately set the time scale of estimating parameters in synchrophasor measurement. To address this problem, a fixed-length AMAF (FAMAF) is proposed that has adjustable filtering capabilities for a fixed data window length.

FAMAF is defined as:

wherein N is_w＝T_w/T_s，T_wIs the window length of the FAMAF, and r is a coefficient that adapts the FAMAF frequency. Equation (3) is equivalent to an impulse response of h_FAMAF(m)＝[1-r,1,......,1-r]Wherein m is 0,1_w-1。h_FAMAFthe transfer function and frequency response of (m) can be described as:

It can be demonstrated that if the value of r in equation (5) is set to:

Then it is possible to obtain:

|H_FAMAF[k(f_w+Δf_w)]|≈0 (7)

Equation (7) shows that FAMAF can be achieved with a fixed data window length T_w＝1/f_wFiltering out f_w+Δf_wIntegral multiples of the interference component, the frequency of the FAMAF may be set accordingly based on the adjustment of the value of r.

the embodiment of the invention explains the effectiveness of the proposed FAMAF by a specific example, wherein a test signal x (t) is 1+ sin (2 pi 65t), the frequency of an oscillation component is 65Hz, and the sampling frequency is 1.44 kHz. Three filters are used to filter the oscillation component. The filter includes NMAF, FAMAF shown in equation (3), and AMAF based on linear interpolation (LIP-AMAF). The data window lengths for NMAF and FAMAF are fixed within a nominal period, i.e., 16.7 ms. The length of the data window of the LIP-AMAF is adjusted to 16.0ms according to the grid frequency. The r value of FAMAF is 0.9223, see equation (6), the a value of LIP-AMAF is 0.1538, see references "L.Wang, Q.Jiang, L.Long, C.Zhang, and Y.Wei," A Novel Phase-locked Loop Based on Frequency Detector and Initial Phase angle Detector, "IEEE trans. Power Electron., vol.28, pp.4538-4549, Oct.2013", where a is a coefficient for adapting the Frequency of LIP-AMAF.

Fig. 1 is a graph comparing frequency responses of various filters provided by an embodiment of the present invention, as shown in fig. 1, and the frequency responses of NMAF, LIP-AMAF, and FAMAF are respectively shown by a dotted line, a solid line, and a dashed line, wherein the frequency responses of three filters are shown in (a) of fig. 1, and a local enlargement of (a) of fig. 1 is shown in (b) of fig. 1. Fig. 2 is a graph comparing filtering results of various filters provided in an embodiment of the present invention, as shown in fig. 2, and dotted lines, solid lines, and chain dotted lines are filtering results of NMAF, LIP-AMAF, and FAMAF, respectively, as shown in (a) of fig. 2, wherein detailed comparison of LIP-AMAF and FAMAF is shown in (b) of fig. 2. As can be seen from fig. 1(a), the three filters have similar magnitude responses. As can be seen from fig. 1(b), the proposed FAMAF attenuates the largest amplitude for the oscillating component at 65 Hz. This means that the FAMAF has the best filtering capability for this oscillating component. The filtering results as shown in fig. 2 verify this. In FIG. 2(a), the waveforms of LIP-AMAF and FAMAF almost coincide. This means that they have good filtering capability for interference, and a detailed comparison of the two filters is given in fig. 2 (b). As can be seen from FIG. 2(b), FAMAF has a stronger filtering capability than LIP-AMAF. Furthermore, the FAMAF has a fixed data window length and therefore a fixed response delay compared to LIP-AMAF.

the embodiment of the invention has the adjustable filtering capability of the fixed-length window by improving on the basis of the moving average filter, and effectively overcomes the defects that the NMAF is weak in anti-interference capability when influenced by frequency offset and the AMAF window is variable in length, so that the time scale of the estimated parameter is difficult to accurately set.

Therefore, in a synchrophasor measurement application, the time scale can be set for the estimated parameters efficiently.

on the basis of the foregoing embodiment, fig. 3 is a structural diagram of a synchronized phasor measurement system according to an embodiment of the present invention, as shown in fig. 3, including:

The filter also comprises a phase-locked loop control loop, and the output end of the phase-locked loop control loop is cascaded with the input end of the filter.

specifically, the synchronous phasor measurement system is configured from an input end as a phase-locked loop control loop, the control loop is configured by a phase-locked loop and a plurality of components arranged in front of the input end of the phase-locked loop, and an output end of the phase-locked loop control loop is cascaded with the filter obtained in the above embodiment to form the whole system.

According to the embodiment of the invention, FAMAF is cascaded behind the phase-locked loop control loop, so that the anti-interference capability of the system under frequency deviation is enhanced.

On the basis of the foregoing embodiment, fig. 4 is an internal structure diagram of a phase-locked loop control loop provided in an embodiment of the present invention, as shown in fig. 4, the phase-locked loop control loop includes a filter combination formed by connecting a plurality of non-adaptive sliding average filters in series, a first park transformation system cascaded with an input end of the filter combination, a second park transformation system cascaded with an output end of the filter combination, and a type 1 phase-locked loop cascaded with an output end of the second park transformation system;

Wherein the filter combination is used for pre-filtering.

In particular, the main body of the architecture of fig. 4 is a type 1 Phase Locked Loop (PLL), which is the simplest PLL, with only one integrator in its control loop, and thus has a high margin of stability and fast dynamic response. However, the interference rejection of type 1 PLLs is limited. In order to solve this problem, a pre-filtering element is designed in the structure of fig. 4, and a filter is added in front of a control loop of the PLL.

Before the PLL control loop of fig. 4, there is a cascade of two stages of Park transforms, where the rotation frequency of the two stages of Park transforms is different and the corresponding rotation angle is also different, i.e. both the rotation frequency and the rotation angle can be used to describe the characteristics of the Park transform rotation coordinate system. The rotating frequency (angle) of the first Park conversion is the nominal frequency f of the power grid₀(Angle θ)_rf＝∫2πf₀dt), the first Park transformation is equivalent to the input voltage vector at frequency (angle) f₀(θ_rf) The rotation is reversed such that after the first Park transformation, the fundamental positive sequence component (FFPS) is at frequency (phase angle) Δ f₁-f₀(Δθ＝θ₁-θ_rf) Converted to quasi-dc terms. The rotation frequency (angle) of the second Park transformation is approximately Δ f (Δ θ) due to the influence of the phase locked loop control loop.

based on the above embodiments, a plurality of NMAF cascades are used as a pre-filtering link of a PLL, fig. 5 is a cascade structure diagram of a plurality of non-adaptive sliding average filters provided by the embodiments of the present invention, and as shown in fig. 5, a cascade structure of 4 NMAFs is used, and each NMAF is connected in series to form a filter combination.

In an embodiment of the invention, the window length of the NMAF is set to T₀/2, wherein T₀is the nominal period of the grid.

The use of NMAF in a PLL as shown in fig. 4, at frequency offset, causes amplitude attenuation. Multiple NMAF cascade structures are applied after the first Park transformation, with Δ f ═ f₁-f₀the FFPS is converted to a quasi-dc term. Therefore, the amplitude attenuation caused by NMAF can be obtained from equation (2)：

wherein N is₀Is the number of sample points for a nominal period. Since this link is made up of four NMAF cascades, the amplitude decays to (K)_NMAF)⁴In FIG. 4, K is used₁Amplitude attenuation can be compensated for, where K₁＝(1/K_NMAF)⁴。

The NMAF cascade structure shown in fig. 4 can block even harmonics after the first Park transformation. Note that the first Park transformation is equivalent to a voltage signal in f₀so that the first even harmonic after Park transformation corresponds to the fundamental negative sequence component (FFNS) and odd harmonics in the electrical signal. However, in the presence of frequency offsets, the filtering capabilities of the NMAF cascade structure may be degraded. Furthermore, the NMAF cascade structure cannot block even harmonics in the electrical signal. The NMAF cascade structure filtering performance is not ideal, and the phase-locked loop output is fluctuated.

To further filter the ripple in the PLL output, a post value filtering stage based on FAMAF is provided in fig. 3. To achieve a more accurate phase angle measurement, a pair of FAMAFs is cascaded after the phase angle of the PLL output. A total of four FAMAFs are used throughout the post-filtering stage. In an embodiment of the invention, the window length of FAMAF is set to T_w＝(N₀+1)/f_s. At a fixed data window length, the FAMAF can effectively filter out the interference component in the PLL output at frequency offsets.

The calculation of r in FAMAF is given below. First, the frequency of the PLL output is passed through a FAMAF to obtain a preliminary frequency estimateIs used in the frequency adaptation process for FAMAF, as shown in fig. 3. The value of r in FAMAF is given by equation (6), where f_w＝f_s/(N₀+1)。

To pairThe response time for the estimation is fast, however, its measurement accuracy may be degraded in a dynamic environment. To solve this problem, fig. 3 gives a more accurate estimate f of the frequency₁。f₁Calculated using the Least Squares Fitting Method (LSFM).

Where Δ θ' (k) is the estimated phase angle error, as shown in fig. 3.Is the step size between two phase angles, L_fIs the number of phase angles used in the LSFM method. F_mIs a constant coefficient and can be calculated. Suppose F_s＝1.44 kHz。L_fF can be obtained when 9₁＝-16，F₂＝-12，F₃＝-8，F₄＝-4，F₅＝0，F₆＝4，F₇＝8，F₈＝12，F₉＝16。

The rate of frequency change (ROCOF) can also be calculated using LSFM:

Wherein f is₁(k) Is the estimated frequency in equation (9),Is the step size between two frequencies, L_rfIs the number of frequencies used in the LSFM.

L_f，L_rfa larger value of (a) indicates a longer data window length for estimating frequency and ROCOF, corresponding to a higher measurement accuracy and a slower dynamic response, and vice versa.

in the configuration shown in fig. 3, the FAMAF is applied after the PLL control loop, so that the frequency feedback branch for the FAMAF does not affect the stability margin of the PLL. Therefore, a wide bandwidth can be set for the control loop of the PLL to achieve a fast dynamic response. K in the examples of the present invention_pIs set to 0.5f_s。

Fig. 6 is a flowchart of a synchronized phasor measurement method for a smart grid according to an embodiment of the present invention, as shown in fig. 6, including:

S1, acquiring a voltage signal output by the smart grid;

And S2, inputting the voltage signal to the synchronous phasor measurement system, and obtaining a phasor measurement result and a frequency measurement result output by the synchronous phasor measurement system.

Specifically, the smart grid synchronized phasor Measurement method provided by the embodiment of the invention is a PLL-based synchronized phasor Measurement Algorithm (a Phase-locked loop based synchronized phased sensor Measurement Algorithm, referred to as PSMA hereinafter) that inputs a voltage signal output by a smart grid to be estimated and measured to a Phase-locked loop control loop, an output end of the Phase-locked loop control loop is cascaded with the FAMAF, and a phasor Measurement result and a frequency Measurement result are output from an output end of the FAMAF.

The embodiment of the invention has adjustable filtering capability based on the FAMAF through the proposed PSMA algorithm, has fixed data window length, enhances the anti-interference capability under frequency deviation by adopting the FAMAF after a phase-locked loop control loop, and effectively measures and estimates the phasor measurement result and the frequency measurement result of the input voltage signal.

On the basis of the above embodiment, the smart grid is a three-phase system or a single-phase system.

In particular, the PSMA method is designed for three-phase systems, pre-filtering of multiple NMAF cascadesThe loop section has high interference resistance to FFNS component. Furthermore, the PSMA can be applied equally to a single-phase system, which can be equivalent to an unbalanced three-phase system with very high FFNS components, as long as the three inputs are in terms of three components v in the three-phase system_a、v_bAnd v_cAccording to v, as shown in FIG. 3_a(t)＝2v(t)，v_b(t)＝v_csetting (t) ═ 0 enables application of PSMA in a single-phase system.

on the basis of the foregoing embodiment, the inputting the voltage signal to the synchronous phasor measurement system to obtain a phasor measurement result and a frequency measurement result output by the synchronous phasor measurement system further includes:

Setting time scales of the phasor measurement result and the frequency measurement result;

The time stamp is obtained by measuring delays of the phasor measurement and the frequency measurement.

The measured delay is a measured delay based on phase angle, amplitude, frequency, and frequency rate of change.

Specifically, the PSMA method is a recursive method, and the synchronous phasor measurement technique requires analyzing and calculating synchronous phasors such as frequency, phase angle, etc. at the time under the same time tag, so the synchronous phasor measurement needs to give not only the measured data, but also a time scale corresponding to the data to make sure that the data is the measured data at a certain time. However, NMAF, FAMAF, PLL control loop and LSFM all generate time delay, so that the data measured at a certain time t is not the data at the real time t, but the data at the time t minus the measurement delay time. the measurement delay time is subtracted from t to set the time scale of the data measured at the time t, so that the time error caused by the measurement delay can be avoided, and the measurement result is more accurate.

in order to accurately set the time scale of the estimated parameter, the measurement delay of the estimated parameter needs to be known. The delay in the measurement of the estimated parameters is mainly caused by the implementation of NMAF, FAMAF, PLL control loop and LSFM. The measured delays produced by the estimation of phase angle, amplitude, frequency and ROCOF are analyzed below.

The measured delay of the phase angle measurement is mainly caused by the NMAF, FAMAF and the control loop of the PLL. The delay of an NMAF or FAMAF generation is (N)_w-1)/2, a total of four NMAF and two FAMAF are used in the phase angle measurement, which results in a delay time of (2N)₀-2)T_s. The control loop is a one-step system with a time constant of 1/k_pDynamic response time and 1/k_pOn, if k_pis set to 0.5f_sThen the measurement delay caused by the control loop is 2T_s. Thus, the total delay time of the phase angle measurement is 2N₀T_s， t-2N₀T_sFor setting the time scale of the phase angle measurement at time t.

the measurement delay of the amplitude measurement is caused by NMAF and FAMAF. The amplitude measurement uses a total of four NMAF and one FAMAF, which generate a time delay of (1.5N)₀-2)T_sThus t- (1.5N)₀-2)T_sFor setting the time scale of the amplitude measurement at time t.

The measurement delay of the frequency measurement is mainly caused by NMAF, FAMAF, PLL control loop and LSFM in equation (9). A total of four NMAFs and two FAMAFs are used in the frequency measurement. Delay time caused by NMAF, FAMAF and PLL is 2N₀T_sThe delay caused by LSFM isThus, the total delay time iswhen in useL_FWhen 9, the total delay time measured is 3N₀T_sThus t-3N₀T_sFor setting the time scale of the frequency measurement at time t.

The measurement delay of ROOF is mainly caused by frequency measurement and LSFM, and the total delay time iswhen in useL_F＝9、L_RFWhen 9, the total delay time measured is 3.5N₀T_s，t-3.5N₀T_sFor setting the time scale for the ROCOF measurement at time t.

On the basis of the above embodiment, the method further includes:

Resetting a rotation angle of a second park transformation system in the synchrophasor measurement system after a preset time.

In particular, the rotation angle of the second Park transformation (Δ θ in fig. 4) is an increasing ramp signal, the value of which must be reset after a certain time to avoid data overflow. For example, when Δ θ (k) > 2M pi (M is an integer and may be preset), Δ θ (k) may be reset to Δ θ (k) ═ Δ θ (k) -2M pi to avoid data overflow. The reset process generates a phase jump in Δ θ (k), which leads to the effect of the FAMAF on Δ θ' (i) (i ═ k, …, k + N in fig. 3₀) Is estimated. To solve this problem, when Δ θ (k) is reset, Δ θ (i) (i ═ k-N) is reset₀…, k-1) front N₀The value is reset to Δ θ (i) -2M pi, which is good for avoiding erroneous phase estimation due to the reset process.

On the basis of the above embodiment, the method further includes: and respectively carrying out performance evaluation on the response time and the measurement precision of the phasor measurement result and the frequency measurement result.

Specifically, after a voltage signal of the smart grid is input to the synchronous phasor measurement system to obtain a phasor measurement result and a frequency measurement result, performance evaluation is performed on the obtained results, and evaluation is performed from two dimensions of response time and measurement accuracy.

(1) Response time

The response time of the PSMA at the step change is evaluated. The change in amplitude, phase angle and frequency of the step change can be modeled as:

x(t)＝X_m[1+k_xf₁(t)]×cos{2π[f₀+k_bf₁(t)]t+k_af₁(t)} (11)

wherein f is₁(t) is a unit step function, k_x，k_a，k_bThe amplitude, phase angle and frequency in the step function are respectively, the nominal frequency is 60Hz, and the sampling frequency is 1.44 kHz.

k_x，k_a，k_bThe value reference of (d) is set to 0.1,. pi./18 and 1. Step changes in frequencies not specified in the list are also considered for evaluating the performance of the PSMA method. Parameter setting of PSMAL_F＝9、L_RF9. The total error of vector estimation (TVE), Frequency Error (FE), ROCOF error (RFE) given in equation (12) are shown in fig. 7 to 9, respectively.

Wherein X_r(t) and X_i(t) is the real and imaginary parts of the theoretical value of the input signal at time t, x_r(t) and x_i(t) is the real and imaginary parts of the measurement.

as shown in Table 1, the results of the alignment are listed according to the standard, showing the response time of TVE, FE and RFE, and the response time of class P/M required by IEEE Standard C37.118 is shown in Table 1. Fig. 7 to 9 and table 1 show that the response time of PSMA satisfies both the P-class and M-class requirements specified by IEEE C37.118.

TABLE 1

The experimental results show that the parametersL_f、L_rfThere is no effect on the response time of the TVE. Therefore, the response time of the TVE always meets the P/M requirement. However, the response time of frequency and ROOF is parametrizedL_f、L_rfThe larger these parameter values, the slower the response time, the higher the frequency and ROCOF measurement accuracy, and vice versa.

(2) Measurement accuracy

In the measurement accuracy evaluation, the proposed PSMA method was compared with the prior art. The performance of the method was evaluated using maximum FE, maximum RFE, maximum Magnitude Error (ME), maximum phase Angle Error (AE), and TVE. In all comparative experiments, the PSMA parameters were setL_F＝9、L_RF9. Therefore, the response time of the PSMA meets the requirements of the P/M class in the IEEE standard. The performance of the PSMA was evaluated at different sampling frequencies under a 50Hz (or 60Hz) system.

1) Comparison with quasi-positive sequence discrete Fourier transform (Qps-DFT) method

The PSMA algorithm was compared to the Qps-DFT algorithm, as described in references "T.Xia and Y.Liu," Single-Phase Angle Measurement in Electric Power Systems ", IEEETrans. Power Syst," vol.24, No.2, pp.844-852, May 2010 ", and" L.Zhan, Y.Liu, J.Culliss, J.ZHao, and Y.Liu, "Dynamic Single-Phase Synchronized Phase and frequency Estimation of the Distribution Level," IEEE Trans.Smart Grid, vol.6, No.4, pp.2013-2022, July 2015 ". For comparison purposes, the nominal frequency was set to 60Hz and an approximate TVE (aTVE) was given for comparison. The aTVE means that in the TVE calculation, the amplitude error is assumed to be zero, and only the angle measurement error is converted into the TVE used in the PMU standard. The test conditions and the measurement accuracy are shown in table 2. For PSMA, the actual TVE is shown in table 2, which is calculated using the estimated magnitude and phase angle in PSMA.

table 2 shows that PSMA outperforms Qps-DFT in most parameters. The accuracy of PSMA is lower than Qps-DFT only at four parameters.

In case 2, the measurement error of PSMA in FE is higher than Qps-DFT. The main reason is that PSMA has faster dynamics than Qps-DFT and is therefore more sensitive to noise. The measurement accuracy of the PSMA can be improved by adding an additional noise filter or increasing the data window length of the NMAF and LSFM in the PSMA structure. For example, whenIs increased toThe maximum FE will be reduced to 0.5 mHz.

in case 3, the error of PSMA in AE, FE and TVE is higher than that of Qps-DFT. But the error of measurement of PSMA (4e-8mHz and 6.9 e-3% TVE) is much smaller than PMU standards (25mHz and 1% TVE). In addition, in practical applications, the fundamental frequency of the power grid may have frequency deviation. The advantages of PSMA are apparent at non-nominal frequencies. For example, the measurement accuracy of PSMA is higher than the Qps-DFT when the fundamental frequency deviates to 59.5/60.5 Hz.

TABLE 2

Remarking: PMUstd.: a synchronous phasor measurement standard; reference signal: a reference signal; AE: phase angle error; FE: frequency error, aTVE; approximate total vector error; TVE; a total vector error; steady state accuracuracy: steady state accuracy; dynamic state accuracycacy: dynamic precision; frequency range: a frequency range; noise: noise; harmonic disorders: harmonic distortion; frequency ramp: a frequency ramp; magnitude modulation: amplitude modulation; phase modulation: phase angle modulation.

2) Comparison with distributed-level phasor measurement unit (D-PMU) method and interpolated discrete Fourier transform (IpDFT) method

PSMA was compared with D-PMU and IpDFT, see references "K.Chauhan, M.V.Reddy, and R.Sodhi," A Novel Distribution-Level Phasor Estimation Algorithm Using empirical wave Transform, "IEEE Trans.Ind.Electron., vol.65, No.10, pp.7984-7995, Oct.2018" and "P.Romano and M.Paolone," Enhanced Interpolated-DFT for Synchrophasor Estimation FPs: Theory, amplification and variation of aPMU Prototype ", IEEE ns.Instrum.63, vol.no. 12, 2824-2836, Dec,2014. For objective comparison, the nominal frequency was 50Hz and the sampling frequency was 6.4 kHz. A total of five experiments were performed, with test conditions referenced:

Case 1 (frequency range test): the test signal is x (t) ═ cos (2 pi ft), and f is more than or equal to 45 and less than or equal to 55 Hz.

Case 2 (frequency ramp test): the test signal is x (t) ═ cos (2 pi ft + pi R)_ft²) Wherein R is_fIs the ramp rate (in Hz/s) with a frequency varying from 49.2Hz to 50.2Hz for 1 second in analog time, R_fIs 1 Hz/s.

case 3 (harmonic distortion test): the test signal is x (t) ═ cos (2 pi f)₀t)+0.1cos(2πhf₀t), h is the harmonic order varying from 2 to 50.

Case 4 (phase modulation test): the test signal is x (t) ═ 1+ k_xcos(2πf_mt)]cos[2πf₀t+k_acos(2πf_mt)]Wherein k is_x＝0，k_a0.1, modulation frequency f_mVarying from 0.1 to 5 Hz.

case 5 (amplitude modulation test): the test signal is x (t) ═ 1+ k_xcos(2πf_mt)]cos[2πf₀t+k_acos(2πf_mt)]Wherein k is_x＝0.1，k_a0, modulation frequency f_mvarying from 0.1 to 5 Hz.

The maximum measurement error is shown in table 3. Table 3 shows that all methods meet PMU criteria. However, PSMA outperforms D-PMU and IpDFT in most parameters.

In table 3, only in the harmonic distortion test of case 3, the measurement accuracy of PSMA is not the highest. However, the measurement error of PSMA is much smaller than the PMU standard. In addition, in practical applications, a frequency deviation may occur in the fundamental frequency of the power grid. The advantages of PSMA can be seen in the non-nominal frequency conditions described in case 3 in table 2.

TABLE 3

3) Comparison with DPMA, CT-DFT, and MDSEA methods

The PSMA method and Dynamic synchronous Phasor measurement Algorithm (DPMA), Clarke transform-Based discrete Fourier transform (CT-DFT), Modified Dynamic synchronous Phasor measurement Algorithm (MDSEA) were tested under power oscillation conditions, as detailed in references "R.K.Mai, Z.Y.He, L.Fu, B.Kirby, and Z.Q.Bo," A Dynamic synchronous Phasor Estimation Algorithm for on Application, "IEEE Transmission. Power Del., vol.25, Pha.2, pp.570-578," Apr. 2010. "," L.Zhan, Y.Liu, and Esy.Liu, "A Clarke Transformation-Based sensor and Frequency Algorithm for on Application," IEEE Transmission Algorithm for free X.640, "A.20111. macro.67, E.2018, Mark, DE.67, E.19, DE.8, E.8, E.J.8, Mark.7, E.8. J.7, E.J.J.. For objective comparison, the nominal frequency was 50Hz and the sampling frequency was 2.4 kHz. Two comparative experiments will be performed in this section.

The test signal is x (t) ═ 1+ k_xcos(2πf_mt)]×cos[2π(f₀+Δf)t+k_αcos(2πf_mt-π)]Wherein f is_mis the modulation frequency, and varies from 0.1Hz to 5Hz, k_xIs the amplitude modulation factor, k_ais the phase angle modulation factor.

Case 1 (amplitude shaking)oscillation test): k is a radical of_x＝0.1，k_a＝0。

Case 2 (phase angle oscillation test): k is a radical of_x＝0，k_a＝0.1。

The test conditions and the measurement accuracy are shown in table 1. Table 4 shows that PSMA provides the most accurate measurement for most parameters.

Only in the first case, PSMA does not have the highest accuracy in FE when the fundamental frequency deviates to 47.5Hz and 45 Hz. But the measurement errors (1.4e-2Hz and 7e-3Hz) are much smaller than the PMU standard (0.3 Hz). Furthermore, the PSMA measures TVE and REF with the highest accuracy under these conditions.

Comparison of measurement accuracy verifies the high accuracy of PSMA. And demonstrates that PSMA meets PMU requirements under all test conditions.

TABLE 4

According to the embodiment of the invention, by comparing the measurement result with various measurement algorithm results in the prior art, the response time and the measurement precision of the PSMA measurement result are integrally superior to those of the existing measurement algorithm results, and the PSMA algorithm can accurately track the frequency and phasor of the power grid under the condition of power grid distortion. The response time and the measurement accuracy meet the requirements of IEEEC 37.118.

the above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.

Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.