具体实施方式

下面结合附图对本发明的具体实施方式作进一步的详细说明。

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。

实施例1

本发明提出了一种适用于电力电子装备组网控制的频率检测方法，如图1所示，所述方法包括：

步骤101：基于预设固定频率的虚拟坐标系，对采集的电力电子并网装备端三相电压瞬时值进行d/q轴电压计算；

步骤102：将计算得到的d/q轴电压与经过一阶广义积分器选频特性后的电压进行外积计算，得到虚拟相角与检测相位角差值的正弦；

步骤103：将所述差值的正弦经过积分累加直到所述差值的正弦为零时，将预设固定频率和所述实际电网频率与所述预设固定频率的偏差相加，检测实际电网频率。

步骤101，具体包括：

通过传感器，采集电力电子并网装备端三相电压瞬时值，u_sa、u_sb、u_sc分别表示采集得到的三相电压的瞬时值。

利用Clarke变换对三相电压u_sa、u_sb、u_sc进行Clarke变换得到三相电压的α轴分量u_sα和β轴分量u_sβ；

Clarke变换的表达式如下：

利用Parker变换，以50Hz固定频率作为虚拟基准频率，即由周期为20ms幅值为2π的锯齿波产生的相位角θ分别对u_sα～u_sβ进行Parker变换得到在固定频率虚拟坐标系(旋转角频率ω₁＝100πrad/s)下三相电压的d轴分量u_sd和q轴分量u_sq；

Parker变换的表达式如下：

其中，相位角θ是一周期为20ms幅值为2π的锯齿波状的相角信号。

步骤102，具体包括：

将所获得的u_sd和u_sq分别经过两个参数相同的一阶广义积分器(也称将降阶广义积分器)(reduced-order generalized integrator,ROGI)，而后可获得所检测到的频率处的d轴电压和q轴电压u′_sd和u′_sq，在频域内的表达式为，

式中，Δω＝2π·Δf为谐振频率，Δω是旋转角频率偏差(虚拟角频率与实际角频率的差值)，Δf为频率偏差(虚拟频率与实际频率的差值)，如附图1所示。u′_sd为电力电子并网装备端三相电压瞬时值的d轴分量对应的在所检测到的频率处的电压幅值，u′_sq为电力电子并网装备端三相电压瞬时值的q轴分量对应的在所检测到的频率处的电压幅值，G(s)为一阶广义积分器，s为拉普拉斯算子，ω_c1为截止频率，一般设置在5-25rad/s间，j为复向量的复数因子。

为利用虚拟相角与检测相角之间的相角差Δθ进行频率估算，考虑到Δθ在零度附近变化，可引入矢量外积。构造虚拟相位角下的电压矢量为：

U_sdq＝u_sd+j·u_sq

经过一阶广义积分器后的电压矢量为：

U′_sdq＝u′_sd+j·u′_sq

U′_sdq＝G(s)·U_sdq

式中，其中，u_sd、u_sq为电压在固定频率下虚拟相位角下的d、q轴分量，u′_sd、u′_sq为检测频率下经过一阶广义积分器后的电压d、q轴分量，j复向量的复数因子。

二者间的相角(U_sdq相角减去U′_sdq相角)的正弦可表示为，

式中，sinΔθ为虚拟相角与检测相位角之间的差值的正弦，为三相电压瞬时值的d轴分量的平方值，为三相电压瞬时值的q轴分量的平方值，为经过一阶积分器的三相电压瞬时值的d轴分量对应的在所检测到的频率处的电压幅值的平方值，为经过一阶积分器的三相电压瞬时值的q轴分量对应的在所检测到的频率处的电压幅值的平方值。

步骤103，具体包括：

将虚拟相角与检测相位角之间的差值的正弦带入积分器或准积分器，积分器或准积分器输出检测频率与实际电网频率的偏差的估测值，将检测频率与实际电网频率的偏差的估测值带入一阶广义积分器进行积分累加，一阶广义积分器输出的值与所检测到的频率处的d轴电压和q轴电压u′_sd和u′_sq计算新的虚拟相角与检测相位角之间的差值的正弦，再将新的虚拟相角与检测相位角之间的差值的正弦带入积分器或准积分器，进行迭代，直至虚拟相角与检测相位角之间的差值的正弦为零，得到检测频率与实际电网频率的偏差，最终，将检测频率与实际电网频率的偏差与固定频率相加，得到实际电网频率。

在稳态情况下，检测频率与实际电网频率相等，则频率偏差为零，G(s)在谐振频率处相角为零度；当出现频率扰动导致实际电网频率下降时，频率偏差为正数，G(s)在谐振频率处相角为负值，通过积分累加，使得检测频率逐步降低，并趋近于实际频率，当检测频率与实际电网频率相等，G(s)相角为零，进入新的稳态；当出现频率扰动导致实际电网频率上升时，频率偏差为负数，G(s)在谐振频率处相角为正值，通过积分累加，使得检测频率逐步上升，并趋近于实际频率，当检测频率与实际电网频率相等，G(s)相角为零，进入新的稳态。

在实际电网频率在49-51H_Z的区间内波动时，由于电网频率偏差角度相对较小，可近似认为，

sinΔθ≈Δθ

这也反应了电网频率f_g的调整趋势，当谐振频率收敛至实际电网频率与虚拟基准频率的差值时，U_sdq和U′_sdq的外积(U_sdq×U′_sdq)为零。引入积分器或者准积分器(表达式如下)，并将U_sdq和U′_sdq的外积(U_sdq×U′_sdq)作为频率估测的误差输入，并与降阶广义积分器构成频率估测闭环控制，得到实现对电网电压频率的无静差估测。

积分器：

准积分器：

式中，H₁(s)为积分器，H₂(s)为准积分器，s为拉普拉斯算子，ω_c2为准积分器的截止频率。

同时，积分节器(准积分调节器)对电网电压频率的估测能起一定的滤波作用，同时本发明采用的降级广义积分器的谐振频率在基频附近，并与频率偏差量相同，对负序、谐波电压扰动能够起到一定的衰减作用，增强了该频率估测方法的抗干扰性。

实施例2

为验证提出的适用于电力电子装备并网控制的频率检测方法理论的正确性，在MATLAB/Simulink仿真软件中搭建了主电路图如图2所示的电网频率估算控制仿真模型。仿真中，电网实际频率在1.0s发生波动。此时将根据控制策略所得到的频率检测值与实际频率相对比，可得到仿真结果如图3所示。

由仿真结果可知，本专利所提出的频率检测方法能在电网频率变化情况下仍保持检测频率对实际频率的准确跟踪，动态响应时间短，能够对电网电频率中的估量值起到一定的滤波作用，抗干扰能力强。

实施例3

本发明提供了一种适用于电力电子装备组网控制的频率检测装置，如图4所示，所述装置包括：

第一计算模块，用于基于预设固定频率的虚拟坐标系，对采集的电力电子并网装备端三相电压瞬时值进行d/q轴电压计算；

第二计算模块，用于将计算得到的d/q轴电压与经过一阶广义积分器选频特性后的电压进行外积计算，得到虚拟相角与检测相位角差值的正弦；

累加和检测模块，用于将所述差值的正弦经过积分累加直到所述差值的正弦为零时，将预设固定频率和所述实际电网频率与所述预设固定频率的偏差相加，检测实际电网频率。

优选的，所述第一计算模块具体用于：

对采集的电力电子并网装备端三相电压瞬时值进行Clarke变换后的结果在预设固定频率的虚拟坐标系下进行Parker变换，得到三相电压瞬时值的d/q轴电压；

所述经过一阶广义积分器选频特性后的d轴电压的计算式如下：

其中，u′_sd为三相电压瞬时值的d轴电压经过一阶广义积分器选频特性后的电压值，s为拉普拉斯算子，ω_c1为一阶广义积分器的截止频率，j为复向量的复数因子，Δω为虚拟角频率与实际角频率的差值，u_sd为三相电压瞬时值在预设固定频率的虚拟坐标系下的d轴电压。

所述经过一阶广义积分器选频特性后的q轴电压的计算式如下：

其中，u′_sq为三相电压瞬时值的q轴电压经过一阶广义积分器选频特性后的电压值，u_sq为三相电压瞬时值在预设固定频率的虚拟坐标系下的q轴电压，s为拉普拉斯算子，ω_c1为一阶广义积分器的截止频率，j为复向量的复数因子，Δω为虚拟角频率与实际角频率的差值。

优选的，所述第二计算模块具体用于：

所述虚拟相角与检测相位角差值的正弦的计算式如下：

其中，sinΔθ为虚拟相角与检测相位角差值的正弦，u_sd为三相电压瞬时值在预设固定频率的虚拟坐标系下的d轴电压，u_sq为三相电压瞬时值在预设固定频率的虚拟坐标系下的q轴电压，u′_sd为三相电压瞬时值的d轴电压经过一阶广义积分器选频特性后的电压值，u′_sq为三相电压瞬时值的q轴电压经过一阶广义积分器选频特性后的电压值。

优选的，所述累加和检测模块具体用于：

步骤1：将所述相位角差值输入至积分器或准积分器，进行积分累加得到实际电网频率与所述预设固定频率的偏差；

步骤2：将所述实际电网频率与所述预设固定频率的偏差分别乘以2π得到实际电网角频率与所述预设固定角频率的偏差；

步骤3：基于一阶广义积分器的选频特性，利用实际电网角频率与所述预设固定角频率的偏差对相应的三相电压瞬时值的d/q轴电压进行低通滤波，得到选频特性后的三相电压瞬时值的d/q轴电压；

步骤4：基于选频特性后的三相电压瞬时值的d/q轴电压和对应的三相电压瞬时值的d/q轴电压计算虚拟相角与检测相位角差值的正弦；

步骤5：若所述差值的正弦不为零，则返回步骤1，否则，将预设固定频率和所述实际电网频率与所述预设固定频率的偏差相加，检测实际电网频率；

所述一阶广义积分器的表达式如下：

其中，G(s)为一阶广义积分器，s为拉普拉斯算子，ω_c1为一阶广义积分器的截止频率，j为复向量的复数因子，Δω为虚拟角频率与实际角频率的差值的估测值；

所述积分器的表达式如下：

其中，H₁(s)为积分器，s为拉普拉斯算子。

所述准积分器的表达式如下：

其中，H₂(s)为准积分器，s为拉普拉斯算子，ω_c2为准积分器的截止频率。

本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。