CN114944649A - 基于电量频谱的电网状态辨识方法 - Google Patents

Abstract

本发明公开了一种基于电量频谱的电网状态辨识方法，属于电能质量分析及信号分析领域。具体地，该方法是基于递归离散傅里叶变换的信号处理的电网状态在线辨识方法，该方法辨识所研究的系统包括逆变器主电路(直流侧电压源，三相桥路逆变电路，三相LCL滤波器)，控制单元，三相线路阻抗(三项线路电感，三相线路电容)，交流侧电网(三相电网，接地模块)和电网状态辨识单元(状态辨识单元和状态输出端)。该方法不需要向电网中注入扰动，仅通过对交流并网侧的电流进行递归离散傅里叶变换，根据其谐波频率的分布规律，判断电网的状态。

Description

基于电量频谱的电网状态辨识方法

技术领域

本发明涉及电能质量分析及信号分析领域，涉及一种基于电量频谱的电网状态辨识方法。

背景技术

随着新能源在电网中渗透率的增加，新能源电站的占有率甚至与传统的火力发电相持平，但是新能源的投入会降低电网的强度；同时由于天然资源的分配不均衡，长距离的输电也会造成电网强度变弱。此时为了提升远距离输电能力，会在输电线路中接入串补设备，来抵消弱电网。也因此线路会出现弱电网和串补状态

由于弱电网和串补电网是两种不同的电网状态，对于逆变器控制所采取的控制也是不同类型。当电网在两种状态中切换时，并网系统的稳定裕度也会有所变化，所以需要辨识电网的两种状态，来自适应的选择控制策略。

目前，对于纯粹的电网状态辨识的研究不多，更多的是对电网阻抗测量的研究，例如：

文献1、题为“Alves D K，Ribeiro R，Costa F B，et al.Wavelet-Based Monitorfor Grid Impedance Estimation of Three-Phase Networks”，IEEE Transactions onIndustrial Electronics，2020，PP(99)：1-1.(基于小波的三相电网阻抗估计，《IEEETransactions on Industrial Electronics》，2020年网络发表)的文章，通过采用小波系数能量分析来检测电网阻抗变化时，通过使用与信号注入方案相关联的实时平稳离散小波包变换来估计电网阻抗。虽然这种方法间歇性的注入谐波，但仍然影响了电网的电能质量。

文献2、题为“基于伪随机二进制序列的宽频带电网阻抗在线测量”在2017年中国高校电力电子与电力传动学术年会中发表的文章。也是通过注入谐波来测量电网阻抗，但是它的创新在于提出了一种新型的dq轴正交信号注入法，在测量精度不变的基础上能够一定程度上减少电网谐波幅值，从而有效提高谐波注入期间的电能质量。虽然这种方法能够一定程度上减少对电网电能质量的影响，但是本质上仍然属于主动测量方法

文献3、题为“基于同步挤压变换的电力系统谐波分析”《太阳能学报》，2021年第8期的文章。文章提出用于谐波分析的傅里叶同步挤压变换和同步挤压小波变换方法。这2种方法分别对电力谐波信号进行时频分析，并对谐波瞬时频率同步挤压锐化，得到更精细的谐波时频曲线，然后利用其可逆性将电力谐波分解为一组内蕴模态类函数分量，完成谐波各分量的提取。文章仅对谐波进行提取，并未进行后续的状态辨识分析。

综合以上文献，现有的技术存在以下不足：

1、现有的辨识方法主要研究重点放在阻抗测量上面，对于直接的电网状态辨识的研究少有研究；

2、现有通过阻抗来辨识电网状态的方法，通常是注入扰动实现的，会影响电网的电能质量，同时，对于弱电网和串补电网的辨识少有存在；

3、现有的状态辨识方法通常是加了硬件外设的，增加了成本。因此有必要研究基于被动测量的弱电网和串补电网辨识方法。

发明内容

本发明提出了一种基于电量频谱的电网状态辨识方法，无需注入扰动，无需增加额外设备，仅通过测量公共耦合点PCC点自身电流谐波来辨识弱电网与串补电网，进而进行下一步的逆变器自适应的控制，提高并网系统稳定性。

本发明的目的是这样实现的。本发明提供了一种基于电量频谱的电网状态辨识方法，应用该方法的系统包括逆变器主电路，控制单元，三相线路阻抗，交流侧电网和电网状态辨识单元；所述逆变器主电路包括依次串联的直流电源、三相全桥逆变电路和三相LCL滤波器，所述三相线路阻抗包括三相线路电感和三相线路电容，所述交流侧电包括三相电网和接地模块，所述电网状态辨识单元包括状态辨识单元和状态输出端；所述逆变器主电路的输出端与三相线路阻抗的一端相连，其连接点为公共耦合点PCC，三相线路阻抗的另一端接入三相电网，三相电网的另一侧与接地模块相连；所述控制单元的输入端与公共耦合点PCC相连、输出端与逆变器主电路中的三相全桥逆变电路相连，所述电网状态辨识单元的输入端与公共耦合点PCC相连、输出端与状态输出端相连；

所述辨识方法通过分析并网电流谐波频率分布规律，对电网状态进行辨识，具体的，本辨识方法的步骤如下：

步骤1，参数的设定

记电网基频为f_g、采样频率为f_s，在一个电网基波周期内按照等时采样间隔τ进行采样，且在初始采样中共采样N次，N为正整数，且N＝f_s/f_g；

选取包含基波在内的一个频段作为频率区间，以1Hz为间隔均分该频率区间，并将均分后的区间简化成N_Har个频率点，将N_Har个频率点中的任一个记为频率点Γ_k，其中k为N_Har个频率点按照频率值由小到大顺序排列的序号，k＝1，2，3，...，N_Har，N_Har为正整数；

引入数组X，X＝{x₁，x₂，...，x_i，...x_N}，其中，x_i为该数组中的任意一个元素，记为元素x_i，i＝1，2，...，N；设定数组X中的所有元素的初始值均为0；

步骤2，初始采样及计算

在频率点Γ_k处采样逆变器主电路输出的A相电流N次，得到N个初始A相电流采样值I₀₁，I₀₂，...，I_0N，将该N个初始A相电流采样值I₀₁，I₀₂，...，I_0N组成一个初始电流值队列E₀，E₀＝{I₀₁，I₀₂，...，I_0N}；

将初始电流值队列E₀中的N个初始A相电流采样值按照排列顺序填充进数组X，并输入到电网状态辨识单元，进行辨识及离散傅里叶变换，得到频率点Γ_k处的余弦分量和频率点Γ_k处的正弦分量，并记为初始余弦分量A_k0和初始正弦分量B_k0，所述离散傅里叶变换公式为：

步骤3，递归采样及递归计算

步骤3.1，给定仿真时长T_m，且在仿真时长中按照等时采样间隔τ共进行递归采样N_c次，N_c为正整数，N_c＝T_m/τ；

在步骤2的初始采样完成后，间隔一个采样间隔τ进行第1次递归采样及计算，具体的，在频率点Γ_k处采样逆变器主电路输出的A相电流1次得到第1次递归A相电流采样值I_A1，将该第1次递归A相电流采样值I_A1和初始电流值队列E₀中的第2个～第N个初始A相电流采样值合并，且该第1次递归A相电流采样值I_A1为队尾，组成一个由N个A相电流采样值组成的第1次递归电流值队列E₁，E₁＝{I₀₂，I₀₃，...，I_0N+1}，其中，I_0N+1＝I_A1；将该第1次递归电流值队列E₁中的N个A相电流采样值按照排列顺序填充进数组X，并根据递归公式得到第1次递归采样时刻频率点Γ_k处的余弦分量和第1次递归采样时刻频率点Γ_k处的正弦分量，并分别记为第1次递归采样时刻余弦分量A_k1和第1次递归采样时刻正弦分量B_k1；

在第1次递归采样完成后，间隔一个采样间隔τ进行第2次递归采样及计算，具体的，在频率点Γ_k处采样逆变器主电路输出的A相电流1次得到第2次递归A相电流采样值I_A2，将该第2次递归A相电流采样值I_A2和第1次递归电流值队列E₁中的第2个～第N个A相电流采样值合并，且该第2次递归A相电流采样值I_A2为队尾，组成一个由N个A相电流采样值组成的第2次递归电流值队列E₂，将该第2次递归电流值队列E₂中的N个A相电流采样值按照排列顺序填充进数组X，并根据递归公式得到第2次递归采样时刻频率点Γ_k处的余弦分量和第2次递归采样时刻频率点Γ_k处的正弦分量，并分别记为第2次递归采样时刻余弦分量A_k2和第2次递归采样时刻正弦分量B_k2；

依次类推，第n次递归采样得到一个第n次递归A相电流采样值I_An和一个由N个A相电流采样值组成的第n次递归电流值队列E_n，n为N_c次递归采样中的任意一次，n＝1，2，...，N_c；将该第n次递归电流值队列E_n中的N个A相电流采样值按照排列顺序填充进数组X，并根据递归公式得到第n次递归采样时刻频率点Γ_k处的余弦分量和第n次递归采样时刻频率点Γ_k处的正弦分量，并分别记为当前时刻余弦分量A_kn和当前时刻正弦分量B_kn；

所述递归公式的表达式如下：

式中，A′_k为当前时刻的上一采样时刻得到的余弦分量，B′_k为当前时刻的上一采样时刻得到的正弦分量，具体的：

当n＝1时，A′_k＝A_k0，B′_k＝B_k0；

当n＞1时，A′_k＝A_kn-1即第n-1次递归采样时刻得到的余弦分量，B′_k＝B_kn-1，即为第n-1次递归采样时刻得到的正弦分量；

步骤3.2，在仿真时长内，按照步骤3.1的方法，每隔一个采样间隔τ采样一次，并进行一次递归计算，最后得到第N_c次采样时刻的余弦分量和正弦分量，并分别记为仿真余弦分量A_k和仿真正弦分量B_k；

将A相电流在频率点Γ_k处的谐波幅值记为谐波幅值Amp_I_k，

步骤4，寻找低频谐波峰值g和高频谐波峰值h

将k＝1，2，3，...，N_Har分别代入步骤2、步骤3进行计算，得到N_Har个频率点对应的谐波幅值Amp_I_k，该N_Har个频率点对应的谐波幅值Amp_I_k组成一个谐波幅值序列Amp_I，Amp_I＝{Amp_I₁，Amp_I₂，...Amp_I_k，...，Amp_I_{N_Har}}；

将该谐波幅值序列Amp_I分成k＜50的谐波幅值区段和k≥50的谐波幅值区段，在k＜50的谐波幅值区段中，取谐波幅值Amp_I_k的最大值记为低频最大值Amp_I_g，其中，g为与低频最大值Amp_I_g对应的频率值，并记为低频谐波峰值g；在k≥50的谐波幅值区段中，取谐波幅值Amp_I_k的最大值记为高频最大值Amp_I_h，其中，h为与高频最大值Amp_I_h对应的频率值，记为高频谐波峰值h；

步骤5，电网状态的辨识

给定阈值ε，给定第一常数a、第二常数b、第三常数c、第四常数d，并设[a，b]为串补电网的特征频段，[c，d]为弱电网的特征频段；

判断低频谐波峰值g和高频谐波峰值h所处的频率是否在特征频段内，且谐波幅值是否达到给定的阈值ε，具体判断如下：

当a≤g≤b且Amp_I_g/Amp_I₅₀＞ε时，认定此时的电网状态为串补电网，输出状态标志位Flag＝2；其中，Amp_I₅₀为基波幅值Amp_I₅₀；

当c≤h≤d且Amp_I_h/Amp_I₅₀＞ε时，认定此时的电网状态为弱电网，输出状态标志位Flag＝1；

电网状态辨识单元中的状态输出端输出状态标志位F1ag，辨识结束。

相对于现有的技术，本发明的有益效果为：

1、本发明辨识方法，无需注入扰动，仅通过测量PCC点自身电流谐波来辨识弱电网与串补电网，进而进行下一步的逆变器自适应的控制，提高并网系统稳定性；

2、本发明辨识方法仅通过分析PCC点电流波形，无需分析电压的波形，判断电流谐波分布情况，就能够判断电网状态；

3、本发明辨识方法仅通过方法改进，无需增加外设硬件设备，无需增加硬件成本。

附图说明

图1为应用本发明辨识方法的系统结构图。

图2为本发明实施例中强电网-弱电网-串补电网的辨识波形。

图3为本发明实施例中的串补电网-弱电网-强电网的辨识波形。

具体实施方式

下面结合附图对本实施例进行具体的描述。

图1是应用本发明辨识方法的系统结构图，由图1可见，应用该方法的系统包括逆变器主电路10，控制单元20，三相线路阻抗30，交流侧电网40和电网状态辨识单元50。所述逆变器主电路10包括依次串联的直流电源、三相全桥逆变电路和三相LCL滤波器，所述三相线路阻抗30包括三相线路电感和三相线路电容，所述交流侧电网40包括三相电网和接地模块，所述电网状态辨识单元50包括状态辨识单元和状态输出端。所述逆变器主电路10的输出端与三相线路阻抗30的一端相连，其连接点为公共耦合点PCC，三相线路阻抗30的另一端接入三相电网，三相电网的另一侧与接地模块相连。所述控制单元20的输入端与公共耦合点PCC相连、输出端与逆变器主电路10中的三相全桥逆变电路相连，所述电网状态辨识单元50的输入端与公共耦合点PCC相连、输出端与状态输出端相连。

在图1中，L₁、L₂为三相LCL滤波器中的滤波电感、C_f为三相LCL滤波器中的滤波电容，L_g为三相线路阻抗中的线路电感，C_g为线路电容，Grid为三相电网，V_dc为直流电源处的直流侧电压。

本实施例中主电路参数为：V_dc＝770V，逆变器额定输出线电压为380V/50Hz，逆变器额定功率为20kW，滤波电感L₁的值为1.832mH，滤波电感L₂的值为0.7mH，滤波电容C_f的值为16.4uF/1.5Ω，强电网时线路电感L_g的值为4.6mH，弱电网时线路电感L_g的值为7mH，串补电网时线路电容C_g的值为16.4uF。

本发明提供了一种基于电量频谱的电网状态辨识方法，所述辨识方法通过分析并网电流谐波频率分布规律，对电网状态进行辨识，具体的，本辨识方法的步骤如下：

步骤1，参数的设定

记电网基频为f_g、采样频率为f_s，在一个电网基波周期内按照等时采样间隔τ进行采样，共采样N次，N为正整数，且N＝f_s/f_g；

选取包含基波在内的一个频段作为频率区间，以1Hz为间隔均分该频率区间，并将均分后的区间简化成N_Har个频率点，将N_Har个频率点中的任一个记为频率点Γ_k，其中k为N_Har个频率点按照频率值由小到大顺序排列的序号，k＝1，2，3，...，N_Har，N_Har为正整数；

引入数组X，X＝{x₁，x₂，...，x_i，...x_N}，其中，x_i为该数组中的任意一个元素，记为元素x_i，i＝1，2，...，N；设定数组X中的所有元素的初始值均为0。

在本实施例中，N＝150，N_Har＝50。

步骤2，初始采样及计算

在频率点Γ_k处采样逆变器主电路10输出的A相电流N次，得到N个初始A相电流采样值I₀₁，I₀₂，...，I_0N，将该N个初始A相电流采样值I₀₁，I₀₂，...，I_0N组成一个初始电流值队列E₀，E₀＝{I₀₁，I₀₂，...，I_0N}；

将初始电流值队列E₀中的N个初始A相电流采样值按照排列顺序填充进数组X，并输入到电网状态辨识单元50，进行辨识及离散傅里叶变换，得到频率点Γ_k处的余弦分量和频率点Γ_k处的正弦分量，并记为初始余弦分量A_k0和初始正弦分量B_k0，所述离散傅里叶变换公式为：

步骤3，递归采样及递归计算

步骤3.1，给定仿真时长T_m，且在仿真时长中按照等时采样间隔τ共进行递归采样N_c次，N_c为正整数，N_c＝T_m/τ；

在步骤2的初始采样完成后，间隔一个采样间隔τ进行第1次递归采样及计算，具体的，在频率点Γ_k处采样逆变器主电路10输出的A相电流1次得到第1次递归A相电流采样值I_A1，将该第1次递归A相电流采样值I_A1和初始电流值队列E₀中的第2个～第N个初始A相电流采样值合并，且该第1次递归A相电流采样值I_A1为队尾，组成一个由N个A相电流采样值组成的第1次递归电流值队列E₁，E₁＝{I₀₂，I₀₃，...，I_0N+1}，其中，I_0N+1＝I_A1；将该第1次递归电流值队列E₁中的N个A相电流采样值按照排列顺序填充进数组X，并根据递归公式得到第1次递归采样时刻频率点Γ_k处的余弦分量和第1次递归采样时刻频率点Γ_k处的正弦分量，并分别记为第1次递归采样时刻余弦分量A_k1和第1次递归采样时刻正弦分量B_k1；

在第1次递归采样完成后，间隔一个采样间隔τ进行第2次递归采样及计算，具体的，在频率点Γ_k处采样逆变器主电路10输出的A相电流1次得到第2次递归A相电流采样值I_A2，将该第2次递归A相电流采样值I_A2和第1次递归电流值队列E₁中的第2个～第N个A相电流采样值合并，且该第2次递归A相电流采样值I_A2为队尾，组成一个由N个A相电流采样值组成的第2次递归电流值队列E₂，将该第2次递归电流值队列E₂中的N个A相电流采样值按照排列顺序填充进数组X，并根据递归公式得到第2次递归采样时刻频率点Γ_k处的余弦分量和第2次递归采样时刻频率点Γ_k处的正弦分量，并分别记为第2次递归采样时刻余弦分量A_k2和第2次递归采样时刻正弦分量B_k2；

依次类推，第n次递归采样得到一个第n次递归A相电流采样值I_An和一个由N个A相电流采样值组成的第n次递归电流值队列E_n，n为N_c次递归采样中的任意一次，n＝1，2，...，N_c；将该第n次递归电流值队列E_n中的N个A相电流采样值按照排列顺序填充进数组X，并根据递归公式得到第n次递归采样时刻频率点Γ_k处的余弦分量和第n次递归采样时刻频率点Γ_k处的正弦分量，并分别记为当前时刻余弦分量A_kn和当前时刻正弦分量B_kn；

所述递归公式的表达式如下：

式中，A′_k为当前时刻的上一采样时刻得到的余弦分量，B′_k为当前时刻的上一采样时刻得到的正弦分量，具体的：

当n＝1时，A′_k＝A_k0，B′_k＝B_k0；

当n＞1时，A′_k＝A_kn-1即第n-1次递归采样时刻得到的余弦分量，B′_k＝B_kn-1，即为第n-1次递归采样时刻得到的正弦分量；

步骤3.2，在仿真时长内，按照步骤3.1的方法，每隔一个采样间隔τ采样一次，并进行一次递归计算，最后得到第N_c次采样时刻的余弦分量和正弦分量，并分别记为仿真余弦分量A_k和仿真正弦分量B_k；

将A相电流在频率点Γ_k处的谐波幅值记为谐波幅值Amp_I_k，

步骤4，寻找低频谐波峰值g和高频谐波峰值h

将k＝1，2，3，...，N_Har分别代入步骤2、步骤3进行计算，得到N_Har个频率点对应的谐波幅值Amp_I_k，该N_Har个频率点对应的谐波幅值Amp_I_k组成一个谐波幅值序列Amp_I，Amp_I＝{Amp_I₁，Amp_I₂，...Amp_I_k，...，Amp_I_{N_Har}}；

将该谐波幅值序列Amp_I分成k＜50的谐波幅值区段和k≥50的谐波幅值区段，在k＜50的谐波幅值区段中，取谐波幅值Amp_I_k的最大值记为低频最大值Amp_I_g，其中，g为与低频最大值Amp_I_g对应的频率值，并记为低频谐波峰值g；在k≥50的谐波幅值区段中，取谐波幅值Amp_I_k的最大值记为高频最大值Amp_I_h，其中，h为与高频最大值Amp_I_h对应的频率值，记为高频谐波峰值h；

步骤5，电网状态的辨识

给定阈值ε，给定第一常数a、第二常数b、第三常数c、第四常数d，并设[a，b]为串补电网的特征频段，[c，d]为弱电网的特征频段；

判断低频谐波峰值g和高频谐波峰值h所处的频率是否在特征频段内，且谐波幅值是否达到给定的阈值ε，具体判断如下：

当a≤g≤b且Amp_I_g/Amp_I₅₀＞ε时，认定此时的电网状态为串补电网，输出状态标志位Flag＝2；其中，Amp_I₅₀为基波幅值Amp_I₅₀；

当c≤h≤d且Amp_I_h/Amp_I₅₀＞ε时，认定此时的电网状态为弱电网，输出状态标志位Flag＝1；

电网状态辨识单元50中的状态输出端输出状态标志位Flag，辨识结束。

在本实施例中，a＝1，b＝15，c＝110，d＝140，ε＝0.00167。

图2是本发明实施例经强电网电变至弱电网，再变至串补电网的辨识结果波形。由图可见，系统在1.6s时缓起完成，此时是强电网状态；在3s时，电网变为弱电网，经过过渡时间0.2s判断在(110Hz，140Hz)特征频段有符合条件的幅值，标志位Flag＝1，辨识为弱电网；在5s时，电网由弱电网变成串补电网，经过过渡时间0.6s判断在(1Hz，15Hz)特征频段由符合条件的幅值，标志位Flag＝2，辨识为串补电网。

图3是本发明实施例经串补电网变至弱电网，再变至强电网的辨识结果波形。由图可见，在1.6s时缓起完成后，辨识出是串补电网，标志位Flag＝2；在5s时，电网变为弱电网，经过过渡时间0.6s判断在(110Hz，140Hz)特征频段有符合条件的幅值，标志位Flag＝1，辨识为弱电网；在8s时，电网由弱电网变成强电网，经过过渡时间0.9s判断在(1Hz，15Hz)、(110Hz，140Hz)特征频段均没有符合条件的幅值，标志位Flag＝0，辨识为强电网。

图2、3证明该发明可以很准确的辨识弱电网和串补电网。

Claims (1)

1.一种基于电量频谱的电网状态辨识方法，应用该方法的系统包括逆变器主电路(10)，控制单元(20)，三相线路阻抗(30)，交流侧电网(40)和电网状态辨识单元(50)；所述逆变器主电路(10)包括依次串联的直流电源、三相全桥逆变电路和三相LCL滤波器，所述三相线路阻抗(30)包括三相线路电感和三相线路电容，所述交流侧电网(40)包括三相电网和接地模块，所述电网状态辨识单元(50)包括状态辨识单元和状态输出端；所述逆变器主电路(10)的输出端与三相线路阻抗(30)的一端相连，其连接点为公共耦合点PCC，三相线路阻抗(30)的另一端接入三相电网，三相电网的另一侧与接地模块相连；所述控制单元(20)的输入端与公共耦合点PCC相连、输出端与逆变器主电路(10)中的三相全桥逆变电路相连，所述电网状态辨识单元(50)的输入端与公共耦合点PCC相连、输出端与状态输出端相连；

其特征在于，所述辨识方法通过分析并网电流谐波频率分布规律，对电网状态进行辨识，具体的，本辨识方法的步骤如下：

步骤1，参数的设定

记电网基频为f_g、采样频率为f_s，在一个电网基波周期内按照等时采样间隔τ进行采样，且在初始采样中共采样N次，N为正整数，且N＝f_s/f_g；

选取包含基波在内的一个频段作为频率区间，以1Hz为间隔均分该频率区间，并将均分后的区间简化成N_Har个频率点，将N_Har个频率点中的任一个记为频率点Γ_k，其中k为N_Har个频率点按照频率值由小到大顺序排列的序号，k＝1，2，3，...，N_Har，N_Har为正整数；

引入数组X，X＝{x₁，x₂，...，x_i，...x_N}，其中，x_i为该数组中的任意一个元素，记为元素x_i，i＝1，2，...，N；设定数组X中的所有元素的初始值均为0；

步骤2，初始采样及计算

在频率点Γ_k处采样逆变器主电路(10)输出的A相电流N次，得到N个初始A相电流采样值I₀₁，I₀₂，...，I_0N，将该N个初始A相电流采样值I₀₁，I₀₂，...，I_0N组成一个初始电流值队列E₀，E₀＝{I₀₁，I₀₂，...，I_0N}；

将初始电流值队列E₀中的N个初始A相电流采样值按照排列顺序填充进数组X，并输入到电网状态辨识单元(50)，进行辨识及离散傅里叶变换，得到频率点Γ_k处的余弦分量和频率点Γ_k处的正弦分量，并记为初始余弦分量A_k0和初始正弦分量B_k0，所述离散傅里叶变换公式为：

步骤3，递归采样及递归计算

步骤3.1，给定仿真时长T_m，且在仿真时长中按照等时采样间隔τ共进行递归采样N_c次，N_c为正整数，N_c＝T_m/τ；

在步骤2的初始采样完成后，间隔一个采样间隔τ进行第1次递归采样及计算，具体的，在频率点Γ_k处采样逆变器主电路(10)输出的A相电流1次得到第1次递归A相电流采样值I_A1，将该第1次递归A相电流采样值I_A1和初始电流值队列E₀中的第2个～第N个初始A相电流采样值合并，且该第1次递归A相电流采样值I_A1为队尾，组成一个由N个A相电流采样值组成的第1次递归电流值队列E₁，E₁＝{I₀₂，I₀₃，...，I_0N+1}，其中，I_0N+1＝I_A1；将该第1次递归电流值队列E₁中的N个A相电流采样值按照排列顺序填充进数组X，并根据递归公式得到第1次递归采样时刻频率点Γ_k处的余弦分量和第1次递归采样时刻频率点Γ_k处的正弦分量，并分别记为第1次递归采样时刻余弦分量A_k1和第1次递归采样时刻正弦分量B_k1；

在第1次递归采样完成后，间隔一个采样间隔τ进行第2次递归采样及计算，具体的，在频率点Γ_k处采样逆变器主电路(10)输出的A相电流1次得到第2次递归A相电流采样值I_A2，将该第2次递归A相电流采样值I_A2和第1次递归电流值队列E₁中的第2个～第N个A相电流采样值合并，且该第2次递归A相电流采样值I_A2为队尾，组成一个由N个A相电流采样值组成的第2次递归电流值队列E₂，将该第2次递归电流值队列E₂中的N个A相电流采样值按照排列顺序填充进数组X，并根据递归公式得到第2次递归采样时刻频率点Γ_k处的余弦分量和第2次递归采样时刻频率点Γ_k处的正弦分量，并分别记为第2次递归采样时刻余弦分量A_k2和第2次递归采样时刻正弦分量B_k2；

依次类推，第n次递归采样得到一个第n次递归A相电流采样值I_An和一个由N个A相电流采样值组成的第n次递归电流值队列E_n，n为N_c次递归采样中的任意一次，n＝1，2，...，N_c；将该第n次递归电流值队列E_n中的N个A相电流采样值按照排列顺序填充进数组X，并根据递归公式得到第n次递归采样时刻频率点Γ_k处的余弦分量和第n次递归采样时刻频率点Γ_k处的正弦分量，并分别记为当前时刻余弦分量A_kn和当前时刻正弦分量B_kn；

所述递归公式的表达式如下：

式中，A′_k为当前时刻的上一采样时刻得到的余弦分量，B′_k为当前时刻的上一采样时刻得到的正弦分量，具体的：

当n＝1时，A′_k＝A_k0，B′_k＝B_k0；

当n＞1时，A′_k＝A_kn-1即第n-1次递归采样时刻得到的余弦分量，B′_k＝B_kn-1，即为第n-1次递归采样时刻得到的正弦分量；

步骤3.2，在仿真时长内，按照步骤3.1的方法，每隔一个采样间隔τ采样一次，并进行一次递归计算，最后得到第N_c次采样时刻的余弦分量和正弦分量，并分别记为仿真余弦分量A_k和仿真正弦分量B_k；

将A相电流在频率点Γ_k处的谐波幅值记为谐波幅值Amp_I_k，

步骤4，寻找低频谐波峰值g和高频谐波峰值h

将k＝1，2，3，...，N_Har分别代入步骤2、步骤3进行计算，得到N_Har个频率点对应的谐波幅值Amp_I_k，该N_Har个频率点对应的谐波幅值Amp_I_k组成一个谐波幅值序列Amp_I，Amp_I＝{Amp_I₁，Amp_I₂，...Amp_I_k，...，Amp_I_{N_Har}}；

将该谐波幅值序列Amp_I分成k＜50的谐波幅值区段和k≥50的谐波幅值区段，在k＜50的谐波幅值区段中，取谐波幅值Amp_I_k的最大值记为低频最大值Amp_I_g，其中，g为与低频最大值Amp_I_g对应的频率值，并记为低频谐波峰值g；在k≥50的谐波幅值区段中，取谐波幅值Amp_I_k的最大值记为高频最大值Amp_I_h，其中，h为与高频最大值Amp_I_h对应的频率值，记为高频谐波峰值h；

步骤5，电网状态的辨识

给定阈值ε，给定第一常数a、第二常数b、第三常数c、第四常数d，并设[a，b]为串补电网的特征频段，[c，d]为弱电网的特征频段；

判断低频谐波峰值g和高频谐波峰值h所处的频率是否在特征频段内，且谐波幅值是否达到给定的阈值ε，具体判断如下：

当a≤g≤b且Amp_I_g/Amp_I₅₀＞ε时，认定此时的电网状态为串补电网，输出状态标志位Flag＝2；其中，Amp_I₅₀为基波幅值Amp_I₅₀；

当c≤h≤d且Amp_I_h/Amp_I₅₀＞ε时，认定此时的电网状态为弱电网，输出状态标志位Flag＝1；

电网状态辨识单元(50)中的状态输出端输出状态标志位Flag，辨识结束。

Images