CN110112776B - 考虑电网背景谐波的并网逆变器电网阻抗辨识方法 - Google Patents

考虑电网背景谐波的并网逆变器电网阻抗辨识方法 Download PDF

Info

Publication number
CN110112776B
CN110112776B CN201910376070.XA CN201910376070A CN110112776B CN 110112776 B CN110112776 B CN 110112776B CN 201910376070 A CN201910376070 A CN 201910376070A CN 110112776 B CN110112776 B CN 110112776B
Authority
CN
China
Prior art keywords
current
voltage
grid
signal
frequency
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.)
Active
Application number
CN201910376070.XA
Other languages
English (en)
Other versions
CN110112776A (zh
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.)
China Electric Power Research Institute Co Ltd CEPRI
Zhengzhou University of Light Industry
Original Assignee
China Electric Power Research Institute Co Ltd CEPRI
Zhengzhou University of Light Industry
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 China Electric Power Research Institute Co Ltd CEPRI, Zhengzhou University of Light Industry filed Critical China Electric Power Research Institute Co Ltd CEPRI
Priority to CN201910376070.XA priority Critical patent/CN110112776B/zh
Publication of CN110112776A publication Critical patent/CN110112776A/zh
Application granted granted Critical
Publication of CN110112776B publication Critical patent/CN110112776B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
    • G01R27/02—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/01—Arrangements for reducing harmonics or ripples
    • H02J3/382—
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/40—Arrangements for reducing harmonics

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Inverter Devices (AREA)

Abstract

本发明提出一种考虑电网背景谐波的并网逆变器电网阻抗辨识方法，首先，利用电压传感器和电流传感器分别测得并网逆变器的电压和桥臂侧电流，通过比例积分控制器对电压和桥臂侧电流进行变换获得调制电压信号；然后，将高频电压信号注入调制电压信号中后再通过控制系统将其输入到并网逆变器中，更新并网逆变器的电压，并利用电流传感器测得并网逆变器的电网侧电流，将更新后的电压和电网侧电流分别代入电压提取模块和电流提取模块；最后，利用改进的复数滤波器提取并网逆变器的电网的高频电压信号和高频电流信号，并对其进行计算得到电网的阻抗值。本发明中改进的复数滤波器考虑了电网的5次和7次高频背景谐波的影响，提高了电网阻抗的辨识精度。

Description

考虑电网背景谐波的并网逆变器电网阻抗辨识方法
技术领域
本发明涉及电力电子技术领域，特别是指一种考虑电网背景谐波的并网逆变器电网阻抗辨识方法。
背景技术
近年来，随着新能源并网逆变器安装规模的快速扩大，电网越来越多的呈现出“弱电网”的特性，其阻抗也越来越大，并对并网逆变器的稳定运行产生了巨大的影响。针对弱电网下并网逆变器的控制，采用基于阻抗的稳定性判据是研究并网逆变器和电网之间交互作用的一项重要方法，而该判据需要获取准确的电网阻抗信息，因此，研究电网阻抗辨识方法具有重要意义。常用的电网阻抗辨识方法主要包括被动法和主动法两类。被动法通过检测电网固有的电压及电流谐波来计算电网阻抗，其优点是不会给电网增加谐波扰动，但是由于信噪比较低，导致被动法的阻抗辨识精度较低。主动法通过向电网中注入特性频率的电压谐波，并提取电网的谐波电流，从而实现电网阻抗辨识。该方法通过注入高频信号可以提高信噪比，从而提高电网阻抗辨识精度，因此其应用更广。
目前已经有许多电网阻抗辨识方法申请了专利，比如申请号为201710113861.4，发明名称为电网阻抗辨识的验证方法及实验装置，申请号为201710361584.9，发明名称为一种基于PRBS扰动注入的电网阻抗在线辨识方法及装置，申请号为201820339286.X，实用新型名称为基于在线阻抗辨识的辨识电路，这些方法都需要在电流指令中叠加高频信号，并通过电流环比例积分控制器将高频电流注入到电网中。由于电流环比例积分控制器只能对直流信号实现无静差跟踪，而无法实现对所注入的高频信号实现无静差跟踪，导致实际注入的高频信号效果较差。此外，由于电网中常含有幅值较大的5次和7次背景谐波，这些谐波的存在也会对电网阻抗的辨识精度产生影响，而上述专利均未考虑电网背景谐波的抑制问题。
文献[谢少军,季林,许津铭.并网逆变器电网阻抗检测技术综述[J].电网技术,2015,39(2):320-326.]总结了现有常用的电网阻抗辨识方法，主要包括主动法、被动法、准被动法等，然而，这些方法也未考虑电网5次和7次背景谐波的影响。
发明内容
针对现有的电网阻抗辨识方法均未考虑电网背景谐波，造成电网阻抗的辨识精度较低的技术问题，本发明提出了一种考虑电网背景谐波的并网逆变器电网阻抗辨识方法，采用改进的复数滤波器分别对电压提取模块和电流提取模块进行提取高频电压信号和高频电流信号，获得电网的阻抗值，本发明采用的复数滤波器中增加了电网5次和7次背景谐波抑制模块，能够提高电网阻抗辨识精度，从而减小了电网背景谐波对电网阻抗辨识的影响。
本发明的技术方案是这样实现的：
一种考虑电网背景谐波的并网逆变器电网阻抗辨识方法，其步骤如下：
S1、利用电压传感器对并网逆变器的电网进行采样得到并网逆变器的线电压ugɑb、线电压ugbc和线电压ugca，计算得到三相电网的相电压uga、相电压ugb和相电压ugc，将相电压uga、相电压ugb和相电压ugc变换到两相静止DQ坐标系中得到电压ugD和电压ugQ，计算得到电网电压角度θ0；
S2、利用电流传感器对并网逆变器的桥臂侧电流进行采样得到三相电流ia、三相电流ib和三相电流ic，并将三相电流ia、三相电流ib和三相电流ic变换到两相静止DQ坐标系中得到两个电流分量分别为电流iD和电流iQ，再利用电网电压角度θ0将电流iD和电流iQ通过坐标变换得到同步旋转dq坐标系上的两个电流分量分别为电流id和电流iq；
S3、设定电流参考值为电流idref和电流iqref，将电流idref、电流iqref和步骤S2得到的电流id、电流iq经过比例积分控制器得到同步旋转dq坐标系中的调制电压信号udref和调制电压信号uqref，再将调制电压信号udref和调制电压信号uqref变换到两相静止DQ坐标系中得到调制电压信号uDref和调制电压信号uQref；
S4、将高频电压信号uDh0和高频电压信号uQh0分别注入步骤S3得到的调制电压信号uDref和调制电压信号uQref中得到两个调制电压信号分别为高频调制电压信号uDhref和高频调制电压信号uQhref，再将高频调制电压信号uDhref和高频调制电压信号uQhref输入空间矢量调制单元，输出6路PWM信号；
S5、通过并网逆变器中控制系统将PWM信号输入到并网逆变器中，更新步骤S1中的电压ugD和电压ugQ；
S6、将步骤S5得到的电压ugD和电压ugQ分别代入uDh提取模块和uQh提取模块，利用改进的复数滤波器分别对uDh提取模块和uQh提取模块进行提取操作，获得高频电压信号uDh和高频电压信号uQh；
S7、利用电流传感器对并网逆变器的电网侧电流进行采样得到三相电流iga、三相电流igb和三相电流igc，并将三相电流iga、三相电流igb和三相电流igc变换到两相静止DQ坐标系中得到两个电流分量分别为电流igD和电流igQ；
S8、将步骤S7得到的电流igD和电流igQ分别代入iDh提取模块和iQh提取模块，通过改进的复数滤波器分别对iDh提取模块和iQh提取模块进行提取操作，获得高频电流信号iDh和高频电流信号iQh；
S9、根据步骤S6得到的高频电压信号uDh、高频电压信号uQh和步骤S8得到的高频电流信号iDh、高频电流信号iQh计算并网逆变器的电网的电阻值
Figure BDA0002051709430000031
和电感值
Figure BDA0002051709430000032
进而获得电网的阻抗值。
优选地，所述步骤S1中的电压ugD和电压ugQ为：
Figure BDA0002051709430000033
其中，
Figure BDA0002051709430000034
则利用电压ugD和电压ugQ计算得到电网电压角度θ0为：
Figure BDA0002051709430000035
优选地，所述步骤S2中的电流id和电流iq为：
Figure BDA0002051709430000036
优选地，所述步骤S3中的调制电压信号uDref和调制电压信号uQref为：
Figure BDA0002051709430000037
其中，
Figure BDA0002051709430000038
k1为比例积分调节器的比例系数，k2为比例积分调节器的积分系数，s为拉普拉斯算子。
优选地，所述步骤S4中的高频调制电压信号uDhref和高频调制电压信号uQhref为：
Figure BDA0002051709430000039
其中，
Figure BDA00020517094300000310
Uh为所注入的高频信号的幅值，t代表时间。
优选地，所述步骤S6中的高频电压信号uDh和高频电压信号uQh的提取方法为：
S61、利用步骤S5得到的电压ugD和电压ugQ分别计算误差电压信号ugDerr1和误差电压信号ugQerr1：
Figure BDA00020517094300000311
其中，ugD和ugQ分别为两相静止DQ坐标系上的电压，uDh和uQh均为待提取的高频电压信号，
Figure BDA0002051709430000041
和
Figure BDA0002051709430000042
均为电网电压正序分量，
Figure BDA0002051709430000043
和
Figure BDA0002051709430000044
均为电网电压谐波分量；
S62、根据步骤S61得到的误差电压信号ugDerr1和误差电压信号ugQerr1计算高频电压信号uDh、高频电压信号uQh、电网电压正序分量
Figure BDA0002051709430000045
电网电压正序分量
Figure BDA0002051709430000046
Figure BDA0002051709430000047
其中，ωhc,u为高频电压信号uDh提取单元和高频电压信号uQh提取单元的截止频率，ωc,u为电网电压正序电压
Figure BDA0002051709430000048
提取单元和电网电压正序电压
Figure BDA0002051709430000049
提取单元的截止频率，ω0为电网电压的同步角频率，
Figure BDA00020517094300000410
θ0为电网电压的角度，
Figure BDA00020517094300000411
j代表虚数；
S63、将步骤S61得到的误差电压信号ugDerr1和误差电压信号ugQerr1通过坐标变换得到同步旋转dq坐标系上的误差电压信号ugderr1和误差电压信号ugqerr1：
Figure BDA00020517094300000412
S64、根据步骤S63得到的误差电压信号ugderr1和误差电压信号ugqerr1，计算电网电压谐波分量
Figure BDA00020517094300000413
和电网电压谐波分量
Figure BDA00020517094300000414
Figure BDA00020517094300000415
其中，ωc6,u为电网电压谐波分量
Figure BDA00020517094300000416
提取单元和电网电压谐波分量
Figure BDA00020517094300000417
提取单元的截止频率；
S65、将步骤S64得到的电网电压谐波分量
Figure BDA00020517094300000418
和电网电压谐波分量
Figure BDA00020517094300000419
变换到两相静止DQ坐标系中，得到两相静止DQ坐标系下的电网电压谐波分量
Figure BDA00020517094300000420
和电网电压谐波分量
Figure BDA00020517094300000421
Figure BDA00020517094300000422
S66、将步骤S62得到的高频电压信号uDh、高频电压信号uQh、电网电压正序分量
Figure BDA00020517094300000423
电网电压正序分量
Figure BDA0002051709430000051
和步骤S65得到的电网电压谐波分量
Figure BDA0002051709430000052
电网电压谐波分量
Figure BDA0002051709430000053
代入步骤S61中，更新误差电压信号ugDerr1和误差电压信号ugQerr1；
S67、重复执行步骤S61至步骤S66，直到达到设定的命令信号时停止运行，输出提取后的高频电压信号uDh和高频电压信号uQh。
优选地，所述步骤S7中的电流igD和电流igQ为：
Figure BDA0002051709430000054
优选地，所述步骤S8中的高频电流信号iDh和高频电流信号iQh的提取方法为：
S81、利用步骤S7得到的电流igD和电流igQ分别计算误差电流信号igDerr1和误差电流信号igQerr1：
Figure BDA0002051709430000055
其中，igD和igQ分别为两相静止DQ坐标系下的电流，iDh和iQh均为待提取的高频电流信号，
Figure BDA0002051709430000056
和
Figure BDA0002051709430000057
均为电网电流正序分量，
Figure BDA0002051709430000058
和
Figure BDA0002051709430000059
均为电网电流谐波分量；
S82、根据步骤S81得到的误差电流信号igDerr1和误差电流信号igQerr1计算高频电流信号iDh、高频电流信号iQh、电网电流正序分量
Figure BDA00020517094300000510
电网电流正序分量
Figure BDA00020517094300000511
Figure BDA00020517094300000512
其中，ωhc,i为高频电流信号iDh和高频电流信号iQh提取单元的截止频率，且ωhc,i＝ωhc,u，ωc,i为电网电流正序分量
Figure BDA00020517094300000513
和电网电流正序分量
Figure BDA00020517094300000514
提取单元的截止频率，ωc,i＝ωc,u；
S83、将步骤S81得到的误差电流信号igDerr1和误差电流信号igQerr1通过坐标变换得到同步旋转dq坐标系上的误差电流信号igderr1和误差电流信号igqerr1：
Figure BDA00020517094300000515
S84、根据步骤S83得到的误差电流信号igderr1和误差电流信号igqerr1，计算电网电流谐波分量
Figure BDA0002051709430000061
和电网电流谐波分量
Figure BDA0002051709430000063
其中，ωc6,i为电网电流谐波分量
Figure BDA0002051709430000064
提取单元和电网电流谐波分量
Figure BDA0002051709430000065
提取单元的截止频率，且ωc6,i＝ωc6,u；
S85、将步骤S84得到的电网电流谐波分量
Figure BDA0002051709430000066
和电网电流谐波分量
Figure BDA0002051709430000067
变换到两相静止DQ坐标系中，得到两相静止DQ坐标系下的电网电压谐波分量
Figure BDA0002051709430000068
和电网电压谐波分量
Figure BDA0002051709430000069
Figure BDA00020517094300000610
S86、将步骤S82得到的高频电流信号iDh、高频电流信号iQh、电网电流正序分量
Figure BDA00020517094300000611
电网电流正序分量
Figure BDA00020517094300000612
和步骤S85得到的电网电流谐波分量
Figure BDA00020517094300000613
电网电流谐波分量
Figure BDA00020517094300000614
带入步骤S81更新误差电流信号igDerr1和误差电流信号igQerr1；
S87、重复执行步骤S81至步骤S86，直到达到设定的命令信号时停止运行，输出提取后的高频电流信号iDh和高频电流信号iQh。
优选地，所述电网的电阻值
Figure BDA00020517094300000615
和电感值
Figure BDA00020517094300000616
为：
Figure BDA00020517094300000617
本技术方案能产生的有益效果：与常规方案相比，本发明在复数滤波器的基础上增加了5次和7次谐波抑制模块，可以消除电网背景谐波里所包含的5次和7次谐波对阻抗辨识的影响，同时，由于电网中的5次和7次谐波等同于同步旋转dq坐标系上的6次谐波，因此，本发明在同步旋转dq坐标系上设计了一种6次谐波抑制模块，从而可对电网背景谐波中的5、7次谐波进行抑制，提高了电网阻抗辨识精度，降低了计算量。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明的电网阻抗辨识模块整体结构示意图。
图2为图1中的高频电压提取模块结构示意图。
图3为图1中的高频电流提取模块结构示意图。
图4为本发明的整体结构示意图。
图5为常规方法的阻抗辨识仿真结果图。
图6为图5中区域A的局部结果图。
图7为本发明的阻抗辨识仿真结果图。
图8为图7中区域B的局部结果图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1和图4所示，一种考虑电网背景谐波的并网逆变器电网阻抗辨识方法，首先，利用电压传感器和电流传感器分别测得并网逆变器的电压和桥臂侧电流，通过比例积分控制器对电压和桥臂侧电流进行变换获得调制电压信号；然后，将高频电压信号注入调制电压信号中后再通过控制系统将其输入到并网逆变器中，更新并网逆变器的电压，并利用电流传感器测得并网逆变器的电网侧电流，将更新后的电压和电网侧电流分别代入电压提取模块和电流提取模块；最后，利用改进的复数滤波器提取并网逆变器的电网的高频电压信号和高频电流信号，并对其进行计算得到电网的阻抗值。其具体步骤如下：
S1、利用电压传感器对并网逆变器的电网进行采样得到并网逆变器电网的线电压ugɑb、线电压ugbc和线电压ugca，并通过公式(1)对线电压ugɑb、线电压ugbc和线电压ugca进行计算得到三相电网的相电压uga、相电压ugb和相电压ugc：
Figure BDA0002051709430000071
再根据公式(2)将相电压uga、相电压ugb和相电压ugc变换到两相静止DQ坐标系中，得到电压ugD和电压ugQ：
Figure BDA0002051709430000081
因此，电网电压角度θ0为：
Figure BDA0002051709430000082
S2、利用电流传感器对并网逆变器的桥臂侧电流进行采样得到三相电流ia、三相电流ib和三相电流ic，并根据公式(3)将三相电流ia、三相电流ib和三相电流ic变换到两相静止DQ坐标系中得到两个电流分量分别为电流iD和电流iQ：
Figure BDA0002051709430000083
再根据公式(4)将电流iD和电流iQ通过坐标变换得到同步旋转dq坐标系上的两个电流分量分别为电流id和电流iq：
Figure BDA0002051709430000084
其中，θ0为电网电压角度。
S3、设定电流参考值为电流idref和电流iqref，并根据公式(5)将电流idref、电流iqref和步骤S2得到的电流id、电流iq经过比例积分控制器变换得到同步旋转dq坐标系中的两个调制电压信号分别为调制电压信号udref和调制电压信号uqref：
Figure BDA0002051709430000085
再根据公式(6)将调制电压信号udref和调制电压信号uqref变换到两相静止DQ坐标系中得到两个调制电压信号分别为调制电压信号uDref和调制电压信号uQref：
Figure BDA0002051709430000086
其中，k1为比例积分调节器的比例系数，k2为比例积分调节器的积分系数，s为拉普拉斯算子。
S4、将高频电压信号uDh0和高频电压信号uQh0分别注入步骤S3得到的调制电压信号uDref和调制电压信号uQref中得到两个调制电压信号分别为高频调制电压信号uDhref和高频调制电压信号uQhref：
Figure BDA0002051709430000091
其中，
Figure BDA0002051709430000092
Uh为所注入的高频信号的幅值，t代表时间；
再将高频调制电压信号uDhref和高频调制电压信号uQhref输入空间矢量调制单元，输出6路PWM信号以控制并网逆变器的运行。
S5、通过并网逆变器中控制系统将PWM信号输入到并网逆变器中，再次执行步骤S1对电压ugD和电压ugQ进行更新。
S6、如图2所示，将步骤S5得到的电压ugD和电压ugQ分别代入uDh提取模块和uQh提取模块，利用改进的复数滤波器分别对uDh提取模块和uQh提取模块进行提取操作，获得高频电压信号uDh和高频电压信号uQh，其具体步骤如下：
S61、利用步骤S5得到的电压ugD和电压ugQ分别计算误差电压信号ugDerr1和误差电压信号ugQerr1：
Figure BDA0002051709430000093
其中，ugD和ugQ分别为两相静止DQ坐标系上的电压，uDh和uQh均为待提取的高频电压信号，
Figure BDA0002051709430000094
和
Figure BDA0002051709430000095
均为电网电压正序分量，
Figure BDA0002051709430000096
和
Figure BDA0002051709430000097
均为电网电压谐波分量；初始时，高频电压信号uDh、高频电压信号uQh、电网电压正序分量
Figure BDA0002051709430000098
电网电压正序分量
Figure BDA0002051709430000099
电网电压谐波分量
Figure BDA00020517094300000910
和电网电压谐波分量
Figure BDA00020517094300000911
的值均设为零。
S62、根据步骤S61中的误差电压信号ugDerr1和误差电压信号ugQerr1计算高频电压信号uDh、高频电压信号uQh、电网电压正序分量
Figure BDA00020517094300000912
电网电压正序分量
Figure BDA00020517094300000913
Figure BDA0002051709430000101
其中，ωhc,u为高频电压信号uDh提取单元和高频电压信号uQh提取单元的截止频率，ωc,u为电网电压正序电压
Figure BDA0002051709430000102
提取单元和电网电压正序电压
Figure BDA0002051709430000103
提取单元的截止频率，ω0为电网电压的同步角频率，
Figure BDA0002051709430000104
其中，θ0为电网电压的角度，
Figure BDA0002051709430000105
j代表虚数。
S63、将步骤S61得到的误差电压信号ugDerr1和误差电压信号ugQerr1通过坐标变换得到同步旋转dq坐标系上的误差电压信号ugderr1和误差电压信号ugqerr1：
Figure BDA0002051709430000106
S64、根据步骤S63得到的误差电压信号ugderr1和误差电压信号ugqerr1，计算电网电压谐波分量
Figure BDA0002051709430000107
和电网电压谐波分量
Figure BDA0002051709430000108
Figure BDA0002051709430000109
其中，ωc6,u为电网电压谐波分量
Figure BDA00020517094300001010
和电网电压谐波分量
Figure BDA00020517094300001011
提取单元的截止频率。
S65、将步骤S64得到的电网电压谐波分量
Figure BDA00020517094300001012
和电网电压谐波分量
Figure BDA00020517094300001013
变换到两相静止DQ坐标系中，得到两相静止DQ坐标系下的电网电压谐波分量
Figure BDA00020517094300001014
和电网电压谐波分量
Figure BDA00020517094300001015
Figure BDA00020517094300001016
S66、将步骤S62得到的高频电压信号uDh、高频电压信号uQh、电网电压正序分量
Figure BDA00020517094300001017
电网电压正序分量
Figure BDA00020517094300001018
和步骤S65得到的电网电压谐波分量
Figure BDA00020517094300001019
电网电压谐波分量
Figure BDA00020517094300001020
代入步骤S61中，更新误差电压信号ugDerr1和误差电压信号ugQerr1。
S67、重复执行步骤S61至步骤S66，直到达到设定的命令信号时停止运行，输出提取后的高频电压信号uDh和高频电压信号uQh。
S7、利用电流传感器对并网逆变器的电网侧电流进行采样得到三相电流iga、三相电流igb和三相电流igc，并根据公式(14)将三相电流iga、三相电流igb和三相电流igc变换到两相静止DQ坐标系中，得到两个电流分量分别为电流igD和电流igQ：
Figure BDA0002051709430000111
S8、如图3所示，将步骤S7得到的电流igD和电流igQ分别代入iDh提取模块和iQh提取模块，通过改进的复数滤波器分别对iDh提取模块和iQh提取模块进行提取操作，获得高频电流信号iDh和高频电流信号iQh，其具体步骤如下：
S81、利用步骤S7得到的电流igD和电流igQ分别计算误差电流信号igDerr1和误差电流信号igQerr1：
Figure BDA0002051709430000112
其中，igD和igQ分别为两相静止DQ坐标系下的电流，iDh和iQh均为待提取的高频电流信号，
Figure BDA0002051709430000113
和
Figure BDA0002051709430000114
均为电网电流正序分量，
Figure BDA0002051709430000115
和
Figure BDA0002051709430000116
均为电网电流谐波分量。
S82、根据步骤S81得到的误差电流信号igDerr1和误差电流信号igQerr1计算高频电流信号iDh、高频电流信号iQh、电网电流正序分量
Figure BDA0002051709430000117
电网电流正序分量
Figure BDA0002051709430000118
Figure BDA0002051709430000119
其中，ωhc,i为高频电流信号iDh提取单元和高频电流信号iQh提取单元的截止频率，且ωhc,i＝ωhc,u，ωc,i为电网电流正序分量
Figure BDA00020517094300001110
提取单元和电网电流正序分量
Figure BDA00020517094300001111
提取单元的截止频率，ωc,i＝ωc,u。
S83、将步骤S81得到的误差电流信号igDerr1和误差电流信号igQerr1通过坐标变换得到同步旋转dq坐标系上的误差电流信号igderr1和误差电流信号igqerr1：
Figure BDA0002051709430000121
S84、根据步骤S83得到的误差电流信号igderr1和误差电流信号igqerr1，计算电网电流谐波分量
Figure BDA0002051709430000122
和电网电流谐波分量
Figure BDA0002051709430000123
Figure BDA0002051709430000124
其中，ωc6,i为电网电流谐波分量
Figure BDA0002051709430000125
提取单元和电网电流谐波分量
Figure BDA0002051709430000126
提取单元的截止频率，且ωc6,i＝ωc6,u。
S85、将步骤S84得到的电网电流谐波分量
Figure BDA0002051709430000127
和电网电流谐波分量
Figure BDA0002051709430000128
变换到两相静止DQ坐标系中，得到两相静止DQ坐标系下的电网电压谐波分量
Figure BDA0002051709430000129
和电网电压谐波分量
Figure BDA00020517094300001210
Figure BDA00020517094300001211
S86、将步骤S82得到的高频电流信号iDh、高频电流信号iQh、电网电流正序分量
Figure BDA00020517094300001212
电网电流正序分量
Figure BDA00020517094300001213
和步骤S85得到的电网电流谐波分量
Figure BDA00020517094300001214
电网电流谐波分量
Figure BDA00020517094300001215
代入步骤S81更新误差电流信号igDerr1和误差电流信号igQerr1。
S87、重复执行步骤S81至步骤S86，直到达到设定的命令信号时停止运行，输出提取后的高频电流信号iDh和高频电流信号iQh。
S9、根据步骤S6得到的高频电压信号uDh、高频电压信号uQh和步骤S8得到的高频电流信号iDh、高频电流信号iQh计算并网逆变器的电网的电阻值
Figure BDA00020517094300001216
和电感值
Figure BDA00020517094300001217
进而获得电网的阻抗值。其中，电网的电阻值
Figure BDA00020517094300001218
和电感值
Figure BDA00020517094300001219
的计算法方法如公式(20)所示：
Figure BDA00020517094300001220
为了验证本发明的有效性，进行了仿真验证。仿真采用并网逆变器的直流侧电压udc为700V，并网逆变器侧输出电感Li为5mH，滤波电容C为15.6μF，阻尼电阻Rd为2Ω，电网角频率ω0为314rad/s，电网相电压幅值为311V，注入的高频信号幅值Uh为121V，注入的高频信号频率为3424rad/s，截止频率ωhc,u和截止频率ωhc,i为400rad/s，截止频率ωc,u和截止频率ωc,i为221rad/s，截止频率ωc6,u和截止频率ωc6,i为221rad/s。仿真时，设定电流idref和电流iqref的参考值分别为40A和0A，电网电阻Rg为1Ω，电网电感Lg为0.6mH，在电网电压中注入了幅值为2V的7次谐波分量。为了验证本发明的有效性，与常规不加5次谐波分量和7次谐波分量的方法进行了对比研究。图5和图6给出了常规方案的仿真结果，图7和图8给出了本发明方案的仿真结果。如图5和图6所示，由于常规方案没有考虑电网中5次背景谐波和7次背景谐波的影响，且由于所注入的高频信号频率与5次谐波分量和7次谐波分量较接近，导致估算的电网阻抗值存在较大的高频波动。如图7和图8所示，本发明通过增加一个5次谐波和7次谐波抑制模块，消除了电网中5次谐波和7次谐波对阻抗辨识的影响，得到的电网阻抗波动较小。
在高频电压信号注入方面，本发明没有选择在电流指令上叠加高频信号，而是直接在两相静止DQ坐标系上的电压调制信号上叠加了两个高频电压信号，从而确保了高频信号的有效注入，且不需要重新修改调制电流环的比例积分调节器。在高频电压信号和高频电流信号提取方式上，本发明在复数滤波器的基础上增加了一个5次和7次谐波抑制模块，从而可以消除电网背景谐波里所包含的5次和7次谐波对阻抗辨识的影响。为了简化计算量，本发明在同步旋转dq坐标系上进行谐波抑制。由于电网中的5次和7次谐波等同于同步旋转dq坐标系上的6次谐波，因此，本发明在同步旋转dq坐标系上设计了一种6次谐波抑制模块，从而可对电网背景谐波中的5次和7次谐波进行抑制，并降低计算量。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (6)

1.一种考虑电网背景谐波的并网逆变器电网阻抗辨识方法，其特征在于，其步骤如下：
S1、利用电压传感器对并网逆变器的电网进行采样得到并网逆变器的线电压ugɑb、线电压ugbc和线电压ugca，计算得到三相电网的相电压uga、相电压ugb和相电压ugc，将相电压uga、相电压ugb和相电压ugc变换到两相静止DQ坐标系中得到电压ugD和电压ugQ，计算得到电网电压角度θ0；
S2、利用电流传感器对并网逆变器的桥臂侧电流进行采样得到三相电流ia、三相电流ib和三相电流ic，并将三相电流ia、三相电流ib和三相电流ic变换到两相静止DQ坐标系中得到两个电流分量分别为电流iD和电流iQ，再利用电网电压角度θ0将电流iD和电流iQ通过坐标变换得到同步旋转dq坐标系上的两个电流分量分别为电流id和电流iq；
S3、设定电流参考值为电流idref和电流iqref，将电流idref、电流iqref和步骤S2得到的电流id、电流iq经过比例积分控制器得到同步旋转dq坐标系中的调制电压信号udref和调制电压信号uqref，再将调制电压信号udref和调制电压信号uqref变换到两相静止DQ坐标系中得到调制电压信号uDref和调制电压信号uQref；
S4、将高频电压信号uDh0和高频电压信号uQh0分别注入步骤S3得到的调制电压信号uDref和调制电压信号uQref中得到两个调制电压信号分别为高频调制电压信号uDhref和高频调制电压信号uQhref，再将高频调制电压信号uDhref和高频调制电压信号uQhref输入空间矢量调制单元，输出6路PWM信号；
S5、通过并网逆变器中控制系统将PWM信号输入到并网逆变器中，更新步骤S1中的电压ugD和电压ugQ；
S6、将步骤S5得到的电压ugD和电压ugQ分别代入uDh提取模块和uQh提取模块，利用改进的复数滤波器分别对uDh提取模块和uQh提取模块进行提取操作，获得高频电压信号uDh和高频电压信号uQh，具体方法为：
S61、利用步骤S5得到的电压ugD和电压ugQ分别计算误差电压信号ugDerr1和误差电压信号ugQerr1：
Figure FDA0002506056750000011
其中，ugD和ugQ分别为两相静止DQ坐标系上的电压，uDh和uQh均为待提取的高频电压信号，
Figure FDA0002506056750000012
和
Figure FDA0002506056750000013
均为电网电压正序分量，
Figure FDA0002506056750000014
和
Figure FDA0002506056750000015
均为电网电压谐波分量；
S62、根据步骤S61得到的误差电压信号ugDerr1和误差电压信号ugQerr1计算高频电压信号uDh、高频电压信号uQh、电网电压正序分量
Figure FDA0002506056750000021
电网电压正序分量
Figure FDA0002506056750000022
Figure FDA0002506056750000023
其中，ωhc,u为高频电压信号uDh提取单元和高频电压信号uQh提取单元的截止频率，ωc,u为电网电压正序电压
Figure FDA0002506056750000024
提取单元和电网电压正序电压
Figure FDA0002506056750000025
提取单元的截止频率，ω0为电网电压的同步角频率，
Figure FDA0002506056750000026
其中，θ0为电网电压的角度，
Figure FDA0002506056750000027
j代表虚数；
S63、将步骤S61得到的误差电压信号ugDerr1和误差电压信号ugQerr1通过坐标变换得到同步旋转dq坐标系上的误差电压信号ugderr1和误差电压信号ugqerr1：
Figure FDA0002506056750000028
S64、根据步骤S63得到的误差电压信号ugderr1和误差电压信号ugqerr1，计算电网电压谐波分量
Figure FDA0002506056750000029
和电网电压谐波分量
Figure FDA00025060567500000210
其中，ωc6,u为电网电压谐波分量
Figure FDA00025060567500000211
提取单元和电网电压谐波分量
Figure FDA00025060567500000212
提取单元的截止频率；
S65、将步骤S64得到的电网电压谐波分量
Figure FDA00025060567500000213
和电网电压谐波分量
Figure FDA00025060567500000214
变换到两相静止DQ坐标系中，得到两相静止DQ坐标系下的电网电压谐波分量
Figure FDA00025060567500000215
和电网电压谐波分量
Figure FDA00025060567500000216
Figure FDA00025060567500000217
S66、将步骤S62得到的高频电压信号uDh、高频电压信号uQh、电网电压正序分量
Figure FDA00025060567500000218
电网电压正序分量
Figure FDA00025060567500000219
和步骤S65得到的电网电压谐波分量
Figure FDA00025060567500000220
电网电压谐波分量
Figure FDA00025060567500000221
代入步骤S61中，更新误差电压信号ugDerr1和误差电压信号ugQerr1；
S67、重复执行步骤S61至步骤S66，直到达到设定的命令信号时停止运行，输出提取后的高频电压信号uDh和高频电压信号uQh；
S7、利用电流传感器对并网逆变器的电网侧电流进行采样得到三相电流iga、三相电流igb和三相电流igc，并将三相电流iga、三相电流igb和三相电流igc变换到两相静止DQ坐标系中得到两个电流分量分别为电流igD和电流igQ：
Figure FDA0002506056750000031
S8、将步骤S7得到的电流igD和电流igQ分别代入iDh提取模块和iQh提取模块，通过改进的复数滤波器分别对iDh提取模块和iQh提取模块进行提取操作，获得高频电流信号iDh和高频电流信号iQh，具体方法为：
S81、利用步骤S7得到的电流igD和电流igQ分别计算误差电流信号igDerr1和误差电流信号igQerr1：
Figure FDA0002506056750000032
其中，igD和igQ分别为两相静止DQ坐标系下的电流，iDh和iQh均为待提取的高频电流信号，
Figure FDA0002506056750000033
和
Figure FDA0002506056750000034
均为电网电流正序分量，
Figure FDA0002506056750000035
和
Figure FDA0002506056750000036
均为电网电流谐波分量；
S82、根据步骤S81得到的误差电流信号igDerr1和误差电流信号igQerr1计算高频电流信号iDh、高频电流信号iQh、电网电流正序分量
Figure FDA0002506056750000037
电网电流正序分量
Figure FDA0002506056750000038
Figure FDA0002506056750000039
其中，ωhc,i为高频电流信号iDh和高频电流信号iQh提取单元的截止频率，且ωhc,i＝ωhc,u，ωc,i为电网电流正序分量
Figure FDA00025060567500000310
和电网电流正序分量
Figure FDA00025060567500000311
提取单元的截止频率，ωc,i＝ωc,u；
S83、将步骤S81得到的误差电流信号igDerr1和误差电流信号igQerr1通过坐标变换得到同步旋转dq坐标系上的误差电流信号igderr1和误差电流信号igqerr1：
Figure FDA0002506056750000041
S84、根据步骤S83得到的误差电流信号igderr1和误差电流信号igqerr1，计算电网电流谐波分量
Figure FDA0002506056750000042
和电网电流谐波分量
Figure FDA0002506056750000043
其中，ωc6,i为电网电流谐波分量
Figure FDA0002506056750000044
提取单元和电网电流谐波分量
Figure FDA0002506056750000045
提取单元的截止频率，且ωc6,i＝ωc6,u；
S85、将步骤S84得到的电网电流谐波分量
Figure FDA0002506056750000046
和电网电流谐波分量
Figure FDA0002506056750000047
变换到两相静止DQ坐标系中，得到两相静止DQ坐标系下的电网电压谐波分量
Figure FDA0002506056750000048
和电网电压谐波分量
Figure FDA0002506056750000049
Figure FDA00025060567500000410
S86、将步骤S82得到的高频电流信号iDh、高频电流信号iQh、电网电流正序分量
Figure FDA00025060567500000411
电网电流正序分量
Figure FDA00025060567500000412
和步骤S85得到的电网电流谐波分量
Figure FDA00025060567500000413
电网电流谐波分量
Figure FDA00025060567500000414
带入步骤S81更新误差电流信号igDerr1和误差电流信号igQerr1；
S87、重复执行步骤S81至步骤S86，直到达到设定的命令信号时停止运行，输出提取后的高频电流信号iDh和高频电流信号iQh；
S9、根据步骤S6得到的高频电压信号uDh、高频电压信号uQh和步骤S8得到的高频电流信号iDh、高频电流信号iQh计算并网逆变器的电网的电阻值
Figure FDA00025060567500000415
和电感值
Figure FDA00025060567500000416
进而获得电网的阻抗值。
2.根据权利要求1所述的考虑电网背景谐波的并网逆变器电网阻抗辨识方法，其特征在于，所述步骤S1中的电压ugD和电压ugQ为：
Figure FDA00025060567500000417
其中，
Figure FDA00025060567500000418
则利用电压ugD和电压ugQ计算得到电网电压角度θ0为：
Figure FDA00025060567500000419
3.根据权利要求1或2所述的考虑电网背景谐波的并网逆变器电网阻抗辨识方法，其特征在于，所述步骤S2中的电流id和电流iq为：
Figure FDA0002506056750000051
其中，
Figure FDA0002506056750000052
4.根据权利要求3所述的考虑电网背景谐波的并网逆变器电网阻抗辨识方法，其特征在于，所述步骤S3中的调制电压信号uDref和调制电压信号uQref为：
Figure FDA0002506056750000053
其中，
Figure FDA0002506056750000054
k1为比例积分调节器的比例系数，k2为比例积分调节器的积分系数，s为拉普拉斯算子。
5.根据权利要求1所述的考虑电网背景谐波的并网逆变器电网阻抗辨识方法，其特征在于，所述步骤S4中的高频调制电压信号uDhref和高频调制电压信号uQhref分别为：
Figure FDA0002506056750000055
其中，
Figure FDA0002506056750000056
Uh为所注入的高频信号的幅值，t代表时间。
6.根据权利要求1所述的考虑电网背景谐波的并网逆变器电网阻抗辨识方法，其特征在于，所述电网的电阻值
Figure FDA0002506056750000057
和电感值
Figure FDA0002506056750000058
为：
Figure FDA0002506056750000059
CN201910376070.XA 2019-05-07 2019-05-07 考虑电网背景谐波的并网逆变器电网阻抗辨识方法 Active CN110112776B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910376070.XA CN110112776B (zh) 2019-05-07 2019-05-07 考虑电网背景谐波的并网逆变器电网阻抗辨识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910376070.XA CN110112776B (zh) 2019-05-07 2019-05-07 考虑电网背景谐波的并网逆变器电网阻抗辨识方法

Publications (2)

Publication Number Publication Date
CN110112776A CN110112776A (zh) 2019-08-09
CN110112776B true CN110112776B (zh) 2020-09-25

Family

ID=67488644

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910376070.XA Active CN110112776B (zh) 2019-05-07 2019-05-07 考虑电网背景谐波的并网逆变器电网阻抗辨识方法

Country Status (1)

Country Link
CN (1) CN110112776B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113702706B (zh) * 2020-12-23 2023-10-31 中国科学院广州能源研究所 一种基于电力电子变换器的电网阻抗测量方法
CN114188983B (zh) * 2021-12-20 2022-06-21 国网黑龙江省电力有限公司佳木斯供电公司 并网逆变器的控制方法
CN115236404B (zh) * 2022-07-07 2024-07-30 新疆海纳同创智能科技有限公司 一种并网逆变器端口阻抗自测量方法
CN115684719A (zh) * 2023-01-03 2023-02-03 国网江西省电力有限公司电力科学研究院 一种新能源并网系统的宽频带耦合谐波分析方法及系统

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012166933A1 (en) * 2011-06-01 2012-12-06 Enphase Energy, Inc. Method and apparatus for grid impedance detection
CN102611143B (zh) * 2012-03-14 2014-03-19 电子科技大学 一种三相并网逆变器并网电流的控制方法
CN106505840B (zh) * 2016-11-30 2019-02-15 云南电网有限责任公司 一种并网光伏逆变器谐波治理方法
CN107026477A (zh) * 2017-05-26 2017-08-08 合肥工业大学 弱电网下具有电网电压前馈滞后补偿的并网逆变器控制方法
CN108306283A (zh) * 2017-12-29 2018-07-20 国网北京市电力公司 微电网系统的控制方法和装置

Also Published As

Publication number Publication date
CN110112776A (zh) 2019-08-09

Similar Documents

Publication Publication Date Title
CN110112776B (zh) 考虑电网背景谐波的并网逆变器电网阻抗辨识方法
CN110429886B (zh) 一种永磁同步电机低速域转子位置辨识方法
CN108923463B (zh) 考虑锁相环的单相lcl型并网逆变器的频率耦合建模方法
CN109742963B (zh) 一种单相脉冲整流器电网电压估算方法
CN106501574A (zh) 一种新型有源电力滤波器谐波电流检测方法
CN106505840A (zh) 一种并网光伏逆变器谐波治理方法
CN102628894B (zh) 一种选择性谐波提取方法
CN110224431B (zh) 并网逆变器系统中减小锁相环影响的控制方法
CN105785788A (zh) 一种快速三相电压锁相环方法及其动态响应性能分析方法
CN103094924A (zh) 基于无锁相环实现电网同步基准正弦的方法及其装置
CN109802385B (zh) 电压源型逆变器的阻抗建模方法
CN110045187B (zh) 基于高频信号注入的并网逆变器电网阻抗辨识方法
CN103546149A (zh) 一种三相电力系统的锁相方法
CN106208063A (zh) 一种有源电力滤波器的自抗扰控制方法和有源电力滤波器
CN106130351A (zh) 一种电动汽车直流充电器输出电压波纹抑制系统及方法
CN107294139A (zh) 一种电流源并网变换器控制方法
CN111435141B (zh) 一种电网阻抗测量装置及方法
CN109617076B (zh) 用于改善pcc处电能质量的电压检测型有源滤波器
CN106655270A (zh) 抑制三相lcl型并网变换器整流启动冲击电流的控制方法
CN114123206B (zh) 面向电动汽车充电站的谐波治理方法
CN115378025A (zh) 一种无锁相环的风电变流器并网控制方法、设备及介质
CN114512991A (zh) 一种有源谐波和无功发生系统及方法
CN110350533B (zh) 计及新能源发电控制策略的改进仿射谐波潮流计算方法
CN115065092A (zh) 单相并网变换器频率耦合调节控制方法
CN114928076A (zh) 无交流电压传感器的虚拟同步机双闭环控制方法

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