CN108270241B - 风机并网逆变器虚拟同步发电机的控制方法 - Google Patents
风机并网逆变器虚拟同步发电机的控制方法 Download PDFInfo
- Publication number
- CN108270241B CN108270241B CN201810116498.6A CN201810116498A CN108270241B CN 108270241 B CN108270241 B CN 108270241B CN 201810116498 A CN201810116498 A CN 201810116498A CN 108270241 B CN108270241 B CN 108270241B
- Authority
- CN
- China
- Prior art keywords
- virtual
- inverter
- voltage
- virtual synchronous
- power
- 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
Links
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 125
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000005284 excitation Effects 0.000 claims abstract description 60
- 238000004804 winding Methods 0.000 claims abstract description 22
- 230000001052 transient effect Effects 0.000 claims abstract description 19
- 238000013016 damping Methods 0.000 claims abstract description 16
- 238000013178 mathematical model Methods 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000009795 derivation Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 5
- 238000011217 control strategy Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 238000001914 filtration Methods 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- H02J3/386—
-
- 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/48—Controlling the sharing of the in-phase component
-
- 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
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Abstract
一种风机并网逆变器虚拟同步发电机的控制方法,通过建立励磁绕组电磁动态方程和同步发电机的转子运动方程,得到适用于风机接入电网的逆变器的虚拟同步发电机控制数学模型,采用递归设计方法设计励磁控制器,通过构造系统的李雅普诺夫函数,求解得到系统平衡点,使得系统在直接补偿虚拟同步发电机的阻尼系数的同时保证在平衡点附近局部稳定,由此得到风机逆变器PWM调制信号的幅值E;采用最大功率点跟踪与同步发电机的转子运动方程结合来求得PWM调制信号的相角δ。由调制波相角δ、电压幅值E构成逆变器调制信号,实现风机并网的虚拟同步发电机控制。该方法具有物理意义明显,性能优越,接入电力系统后暂态稳定性和动态稳定性得到提升。
Description
技术领域
本发明涉及电力电子变流器控制,特别是一种风机并网逆变器虚拟同步发电机的控制方法。
背景技术
人类社会可持续发展面临着能源安全、环境污染、气候变化等诸多挑战。化石能源大量、广泛的使用,在创造了工业文明的同时也造成环境污染,气候变暖,生态恶化,对人类的生存与发展构成了严重威胁,应对这些挑战的关键是加快能源变革,开发清洁能源,优化能源结构,提高能源效率,具体来讲就是要大力发展风能、太阳能、地热能、潮汐能等可再生能源。我国由于独特的地理优势,陆上可利用风力资源达2.53亿千瓦,海上可利用风力资源达7.5亿千瓦,合计共达10亿千瓦,其分布广泛,可利用风力资源相当丰富,2013年总装机容量一度跃居世界第一。我国风力发电事业潜力巨大,发展迅速,前景广阔。
随着分布式的风力发电机大量接入电网,由于其具有间歇性、随机性、不稳定性,会造成发电功率与负荷需求不匹配的问题,从而影响电力系统的安全运行和电能质量。分布式电源一般通过电力电子器件比如整流器、逆变器将电能并网或者输送到负荷,相比于运行成熟的传统发电设备(如同步发电机),电力电子变换器具有响应迅速,控制灵活等优势,但是也有缺乏惯性和阻尼等缺点,从而为电力系统的稳定性带来了一系列的挑战。
在传统的电网运行过程中,由于同步发电机的转子固有的惯性,当电源侧和负载侧有功功率不平衡时,同步机转子转速变化从而补偿系统的功率缺额,以平缓电网频率波动。从这个角度看,同步发电机具有充足的旋转备用容量和转动惯量,能够在负荷变化时为配电网提供必要的电压和频率支撑,即“电网友好型”的电源单元。为了得到“电网友好型”的风电单元,可以模拟并网逆变器具有同步发电机的外特性,这样就可以方便分布式风机并网,而不至于降低电力系统稳定性。这种模拟同步发电机外特性的控制方法被称为虚拟同步发电机(virtual synchronous generator,VSG)技术。该技术主要通过模拟同步发电机的机械方程和电磁方程来控制并网逆变器,从而使得其外特性同步机别无二致,与此同时还可以将传统电网的运行控制策略移植到含大量分布式风机作为电源的电网中。
在传统电力系统中,为改善电力系统稳定性,通常采用对发电机组的励磁控制,可以有效地维持发电机端电压的调节精度在一个给定的范围内,传统同步机组的励磁调节一般应用自动电压调节器AVR(automatic voltage regulator),其原理是利用发电机端电压偏差进行PID(proportion integral derivative regulating)调节。基于虚拟同步机控制策略的风机并网逆变器也大多采用PID调节手段,这种调节方式优点是保证了端电压的调节精度,也存在着降低电力系统静态与暂态稳定性的缺点。在随后的励磁调节器设计方面,新的控制律利用微分几何方法,对非线性系统采用精确的反馈线性化方法,通过非线性反馈将非线性电力系统转化为线性系统,然后按照线性系统的设计理论进行励磁调节的设计,但是对模型的精确度要求较高,在面对系统参数或结构不确定等情况下,控制效果不好。
现有的并网型风力发电机组大多基于有功、无功与同步坐标轴系下d轴和q轴电流的关系,实现对机组的有功、无功控制。这种功率控制方式使得风力机组对电网表现出反负荷的接口特性,使得风电场内风电机组集群只能做电流源汇集,对电力系统的机械惯量无贡献,使电力系统的电压和频率的稳定性降低。
发明内容
为应对上述传统方法的不足,本发明提供一种风机并网逆变器虚拟同步发电机的控制方法,能够直接补偿虚拟同步机阻尼系数的控制策略,以通过风机并网逆变器改善电力系统暂态稳定性和动态品质。该控制器构造的李雅普诺夫函数采用了递归的设计方法,并充分利用了虚拟同步机本身的暂态能量函数。
本发明的技术解决方案如下:
一种风机并网逆变器虚拟同步发电机的控制方法,其特点在于,包括以下步骤:
1)获取逆变器输入输出端口的电压、电流、频率特征量,计算得到虚拟同步发电机输入及输出功率;
2)结合最大功率点跟踪建立背靠背风机并网系统的逆变器的虚拟同步发电机的转子运动方程,得到逆变器PWM调制波相角;
3)建立虚拟同步机的虚拟励磁绕组电磁动态方程;
4)建立具有励磁控制的虚拟同步发电机单机无穷大系统数学模型并定义其平衡点;
5)递归设计直接补偿阻尼系数的虚拟励磁控制器;
6)结合虚拟同步机转子运动方程得到的虚拟功角δ、逆变器调制波电压幅值E得出逆变器调制信号,实现风机并网的虚拟同步发电机控制。
所述获取所述虚拟同步发电机输出到电网端的电压电流以及电网频率,根据所测电压电流计算虚拟同步机输出功率,包括:
1)获取虚拟同步机输出电压电流及电网频率:
通过电压电流互感器测量风机并网逆变器在三相abc坐标系下的输出侧电压Uabc及输出电流Iabc;通过锁相环模块测量电网频率ωg;
2)计算虚拟同步机输出功率:
利用前述所测输出电压Uabc、输出电流Iabc,按式(1)计算虚拟同步机在abc坐标系输出功率,或将电压电流转换到dq坐标系再进行计算:
Pout=UaIa+UbIb+UcIc=UdId+UqIq (1)
式中,Ud、Uq,Id、Iq分别为母线电压和母线电流在d轴和q轴分量。
所述根据最大功率点跟踪(MPPT)建立背靠背风机并网系统的逆变器的虚拟同步发电机的转子运动方程,得到逆变器PWM调制波相角,包括:
1)根据当前风速和风机的转子转速计算获得的MPPT指令值,结合有功频率下垂控制器的下垂系数的调节信号,通过对下垂系数的改变,使得虚拟拖动转矩的指令值等于MPPT指令值,然后根据虚拟转轴的机械运动特性实现调速器功能:
对于有功-频率下垂控制环节,即
式中,Pset为虚拟同步控制下逆变器的参考有功输入或上层调度指令;Dp为下垂控制系数(由MPPT指令值决定);ωref为预先设定角速度参考值;
2)建立虚拟同步机的转子运动特性方程为:
式中,δ[rad]为逆变器虚拟功角,表示无穷大母线电压与逆变器出口电压之间的夹角;Kd为阻尼系数;ωr[p.u.]为虚拟转子角速度与同步角速度的偏差;ωs=2πf为同步角速度;Pset[p.u.]为虚拟同步控制下逆变器的参考功率指令或上层调度功率指令,对应于传统同步机的机械功率Pm;Pe[p.u.]为逆变器输出电磁功率;
3)基于上式(2)、(3),建立有功频率控制模型,得到虚拟转子角速度及调制波信号相角为:
式中,δ[rad]为逆变器虚拟功角,也即控制逆变器的PWM调制波信号相角。
所述虚拟同步机的虚拟励磁绕组电磁动态方程,包括:
1)采用与同步机励磁绕组电磁动态方程中参数相对应的方式,建立虚拟同步发电机励磁绕组电磁动态方程为:
其中,K′d是励磁绕组的时间常数(s);Vset[p.u.]为对应于系统稳态运行时的给定励磁电压;uf为对应于励磁电压的调节量;Vq是逆变器出口电压,对应于发电机空载感应电动势;是暂态电势,Vq与V′q存在关系为:
Vq=V′q+(xvir-x′vir)Id (7)
其中,xvir表示虚拟定子电抗,可用虚拟阻抗方法进行模拟;x′vir表示虚拟暂态同步电抗,同样可用虚拟阻抗方法模拟;Id为母线电流d轴分量,可表示为:
其中,U为无限大系统的母线电压,可看作常数,x′d∑=x′vir+xl表示为虚拟定子暂态电抗与线路电抗之和。
所述的建立具有励磁控制的虚拟同步发电机单机无穷大系统数学模型并定义其平衡点,包括:
1)根据虚拟同步机虚拟励磁绕组电磁动态方程,将(7)代入(8),得到Vq与V′q之间有如下关系成立:
2)已知电动势、电压和电流的关系为:
将(10)代入(1)得到虚拟同步机送入无穷大系统的有功功率表达式为:
3)根据虚拟同步机转子运动方程(1)和虚拟励磁绕组电磁动态方程(2),最终建立具有励磁控制的虚拟同步发电机单机无穷大系统数学模型为:
5)根据控制理论中对系统平衡点的定义,即系统的平衡点是指满足f(x)≡0的点。定义系统(12)的状态变量为δ,ωr,V′q,并且其平衡点用(δs,ωrs,V′qs)表示,根据该定义,系统(12)的平衡点应满足如下条件:
当给定逆变器预设功率Pset和虚拟励磁电压Vset时,求解(13)非线性代数方程得到系统平衡点。
所述的直接补偿阻尼系数的虚拟励磁控制器的递归设计方案,具体包括:
1)考虑化简后的系统模型(12)中的子系统(δ,ωr),并将暂态电势V′q表示为V′q=V′qs+ΔV′q,则定义ΔV′q=V′q-V′qs,子系统(δ,ωr)可写为:
其中,Pe0=a1V′qssinδ-a2sin2δ、M=2H。
2)将ΔV′q看作该子系统的虚拟控制律,通过调节逆变器虚拟励磁使q轴暂态电势的变化ΔV′q满足a1ΔV′qsinδ=M·kDωr,由于存在调节误差,定义调节误差为:
e=M·kDωr-a1ΔV′qsinδ (15)
其中,kD为自定义变量且kD>0,MkD理想补偿逆变器虚拟阻尼系数。由此可得:
对上式求导得到:
选取控制律r为:
r=-ωsωr-kee (19)
由此,得到直接改善虚拟同步机阻尼系数的励磁控制律u为:
所述的步骤6)是根据得到的控制风机并网逆变器的PWM调制波信号相角δ;再将控制律u(21)代入式(9)(10),通过递归计算得逆变器出口电压参考值Vq,逆变器调制波电压幅值E可由(22)求出:
Vq=Vdc·E (22)
从而得到逆变器调制波电压幅值E。即可得到PWM调制波信号E∠δ,作用于风机并网逆变器,实现风机并网的虚拟同步发电机控制。
本发明方法获得的调制波与载波作比较生成的脉冲信号用于控制风机并网的三相逆变器,实现风机并网的虚拟同步发电机控制。
附图说明
图1是本发明适用的风机并网单元的结构图。
图2是虚拟同步控制下逆变器对应传统同步机示意图。
图3是本发明控制方法流程图。
图4整体控制流程示意图。
具体实施方式
下面结合附图对本发明作进一步说明,但不应以此限制本发明的保护范围。
图1为适用于该方法的直驱风机并网单元的结构图。图中直驱风机输出的交流电先经三相整流器整流,再经三相逆变器逆变输出并经过滤波环节后传输到电网。控制逆变器内电力电子器件行为的信号为本发明的虚拟同步发电机控制方法得到的调制波信号E∠δ。
图2是本发明提出的控制方法流程图。本发明风机并网逆变器虚拟同步发电机的控制方法,包括以下步骤:
1)获取逆变器输入输出端口的电压、电流、频率特征量,计算得到虚拟同步发电机输入及输出功率;
2)结合最大功率点跟踪建立背靠背风机并网系统的逆变器的虚拟同步发电机的转子运动方程,得到逆变器PWM调制波相角;
3)建立虚拟同步机的虚拟励磁绕组电磁动态方程;
4)建立具有励磁控制的虚拟同步发电机单机无穷大系统数学模型并定义其平衡点;
5)递归设计直接补偿阻尼系数的虚拟励磁控制器;
6)结合虚拟同步机转子运动方程得到的虚拟功角δ、逆变器调制波电压幅值E得出逆变器调制信号,实现风机并网的虚拟同步发电机控制。
所述获取所述虚拟同步发电机输出到电网端的电压电流以及电网频率,根据所测电压电流计算虚拟同步机输出功率,包括:
1)获取虚拟同步机输出电压电流及电网频率:
通过电压电流互感器测量风机并网逆变器在三相abc坐标系下的输出侧电压Uabc及输出电流Iabc;通过锁相环模块测量电网频率ωg;
(2)计算虚拟同步机输出功率:
利用前述所测输出电压Uabc、输出电流Iabc,按式(1)计算虚拟同步机在abc坐标系输出功率,或将电压电流转换到dq坐标系再进行计算:
Pout=UaIa+UbIb+UcIc=UdId+UqIq (1)
式中,Ud、Uq,Id、Iq分别为母线电压和母线电流在d轴和q轴分量。
所述根据最大功率点跟踪(MPPT)建立背靠背风机并网系统的逆变器的虚拟同步发电机的转子运动方程,得到逆变器PWM调制波相角,包括:
1)根据当前风速和风机的转子转速计算获得的MPPT指令值,结合有功频率下垂控制器的下垂系数的调节信号,通过对下垂系数的改变,使得虚拟拖动转矩的指令值等于MPPT指令值,然后根据虚拟转轴的机械运动特性实现调速器功能:
对于有功-频率下垂控制环节,即
式中,Pset为虚拟同步控制下逆变器的参考有功输入或上层调度指令;Dp为下垂控制系数(由MPPT指令值决定);ωref为预先设定角速度参考值;
2)建立虚拟同步机的转子运动特性方程为:
式中,δ[rad]为逆变器虚拟功角,表示无穷大母线电压与逆变器出口电压之间的夹角;Kd为阻尼系数;ωr[p.u.]为虚拟转子角速度与同步角速度的偏差;ωs=2πf为同步角速度;Pset[p.u.]为虚拟同步控制下逆变器的参考功率指令或上层调度功率指令,对应于传统同步机的机械功率Pm;Pe[p.u.]为逆变器输出电磁功率;
3)基于上式(2)、(3),建立有功频率控制模型,得到虚拟转子角速度及调制波信号相角为:
式中,δ[rad]为逆变器虚拟功角,也即控制逆变器的PWM调制波信号相角。
所述虚拟同步机的虚拟励磁绕组电磁动态方程,包括:
1)采用与同步机励磁绕组电磁动态方程中参数相对应的方式,建立虚拟同步发电机励磁绕组电磁动态方程为:
其中,K′d是励磁绕组的时间常数(s);Vset[p.u.]为对应于系统稳态运行时的给定励磁电压;uf为对应于励磁电压的调节量;Vq是逆变器出口电压,对应于发电机空载感应电动势;是暂态电势,Vq与V′q存在关系为:
Vq=V′q+(xvir-x′vir)Id (7)
其中,xvir表示虚拟定子电抗,可用虚拟阻抗方法进行模拟;x′vir表示虚拟暂态同步电抗,同样可用虚拟阻抗方法模拟;Id为母线电流d轴分量,可表示为:
其中,U为无限大系统的母线电压,可看作常数,x′d∑=x′vir+xl表示为虚拟定子暂态电抗与线路电抗之和。
所述的建立具有励磁控制的虚拟同步发电机单机无穷大系统数学模型并定义其平衡点,包括:
1)根据虚拟同步机虚拟励磁绕组电磁动态方程,将(7)代入(8),得到Vq与V′q之间有如下关系成立:
2)已知电动势、电压和电流的关系为:
将(10)代入(1)得到虚拟同步机送入无穷大系统的有功功率表达式为:
3)根据虚拟同步机转子运动方程(1)和虚拟励磁绕组电磁动态方程(2),最终建立具有励磁控制的虚拟同步发电机单机无穷大系统数学模型为:
4)根据控制理论中对系统平衡点的定义,即系统的平衡点是指满足f(x)≡0的点。定义系统(12)的状态变量为δ,ωr,V′q,并且其平衡点用(δs,ωrs,V′qs)表示,根据该定义,系统(12)的平衡点应满足如下条件:
当给定逆变器预设功率Pset和虚拟励磁电压Vset时,求解(13)非线性代数方程得到系统平衡点。
所述的直接补偿阻尼系数的虚拟励磁控制器的递归设计方案,具体包括:
1)考虑化简后的系统模型(12)中的子系统(δ,ωr),并将暂态电势V′q表示为V′q=V′qs+ΔV′q,则定义ΔV′q=V′q-V′qs,子系统(δ,ωr)可写为:
其中,Pe0=a1V′qssinδ-a2sin2δ、M=2H。
2)将ΔV′q看作该子系统的虚拟控制律,希望通过调节逆变器虚拟励磁使q轴暂态电势的变化ΔV′q满足a1ΔV′qsinδ=M·kDωr,由于存在调节误差,定义调节误差为:
e=M·kDωr-a1ΔV′qsinδ (15)
其中,kD为自定义变量且kD>0,MkD理想补偿逆变器虚拟阻尼系数。由此可得:
对上式求导得到:
选取控制律r为:
r=-ωsωr-kee (19)
由此,得到直接改善虚拟同步机阻尼系数的励磁控制律u为:
所述的步骤6)是根据得到的控制风机并网逆变器的PWM调制波信号相角δ;再将控制律u(21)代入式(9)(10),通过递归计算得逆变器出口电压参考值Vq,逆变器调制波电压幅值E可由(22)求出:
Vq=Vdc·E (22)
从而得到逆变器调制波电压幅值E。即可得到PWM调制波信号E∠δ,作用于风机并网逆变器,实现风机并网的虚拟同步发电机控制。
图3是本发明控制方法流程图,其步骤包括:
S10,测量逆变器输入输出端口的电压电流频率特征量,计算得到虚拟同步发电机输入及输出功率;图1为适用于该方法的风机并网系统结构图。图中风机输出的频率幅值变化的交流电通过三相整流器得到幅值恒定的直流电,类似地,该三相整流器也可以采用虚拟同步电动机的控制策略;接着再经DC/AC逆变器输出,通过LCL滤波环节后输送到电网。控制逆变器内电力电子器件导通和关断的PWM脉冲由本发明控制方法得到的调制波信号E∠δ产生,以实现风机并网的VSG控制。
S20,结合最大功率点跟踪建立背靠背风机并网系统的逆变器的虚拟同步发电机的转子运动方程,如图4的上半部分所示,其中MPPT控制算法是通过结合有功/频率下垂控制器的下垂系数的调节信号,通过对下垂系数的改变,使得虚拟拖动转矩的指令值等于MPPT指令值,然后根据虚拟转轴的机械运动特性实现调速器功能。再由虚拟同步发电机的转子运动方程得到逆变器PWM调制波相角;
S30建立虚拟同步机的虚拟励磁绕组电磁动态方程;
S40建立具有励磁控制的虚拟同步发电机单机无穷大系统数学模型并定义其平衡点:
具体地,该步骤包括:参照图2,结合步骤二、步骤三得到的转子运动方程、电磁动态方程,推导出虚拟同步发电机控制下的逆变器系统数学模型表达式。根据虚拟同步发电机控制下的逆变器系统数学模型表达式,定义具有励磁控制的虚拟同步发电机单机无穷大系统平衡点,并通过解非线性代数方程得到系统平衡点,用于虚拟励磁控制器的设计。
S50,结合步骤四所得模型及其平衡点,如图4所示,结合前面的有功频率控制模型得出的ωr、Pset、δ等参数,通过构造李雅普诺夫函数,递归设计直接补偿阻尼系数的虚拟励磁控制器,即通过图4中递归计算得到用于控制风机并网逆变器的PWM调制波信号幅值。
S60,根据S50所设计的控制律,得到用于逆变器控制的调制波信号幅值E,如图4所示。结合S20得到的虚拟功角δ,即可生成调制波信号E∠δ。该信号与载波作比较,生成的脉冲信号用于控制风机逆变器开关管的导通和关断。
Claims (6)
1.一种风机并网逆变器虚拟同步发电机的控制方法,其特征在于,包括以下步骤:
1)获取逆变器输入输出端口的电压、电流、频率特征量,计算得到虚拟同步发电机输入及输出功率;
2)结合最大功率点跟踪建立背靠背风机并网系统的逆变器的虚拟同步发电机的转子运动方程,得到逆变器PWM调制波相角,具体包括:
2.1)根据当前风速和风机的转子转速计算获得的MPPT指令值,结合有功频率下垂控制器的下垂系数的调节信号,通过对下垂系数的改变,使得虚拟拖动转矩的指令值等于MPPT指令值,然后根据虚拟转轴的机械运动特性实现调速器功能:
对于有功-频率下垂控制环节,即
式中,Pset为虚拟同步控制下逆变器的参考有功输入或上层调度指令;Dp为下垂控制系数,由MPPT指令值决定;ωref为预先设定角速度参考值;
2.2)建立虚拟同步机的转子运动特性方程为:
式中,δ[rad]为逆变器虚拟功角,表示无穷大母线电压与逆变器出口电压之间的夹角;Kd为阻尼系数;ωr[p.u.]为虚拟转子角速度与同步角速度的偏差;ωs=2πf为同步角速度;Pset[p.u.]为虚拟同步控制下逆变器的参考功率指令或上层调度功率指令,对应于传统同步机的机械功率Pm;Pe[p.u.]为逆变器输出电磁功率;
2.3)基于上式(2)、(3),建立有功频率控制模型,得到虚拟转子角速度及调制波信号相角为:
式中,δ[rad]为逆变器虚拟功角,也即控制逆变器的PWM调制波信号相角;
3)建立虚拟同步机的虚拟励磁绕组电磁动态方程;
4)建立具有励磁控制的虚拟同步发电机单机无穷大系统数学模型并定义其平衡点;
5)递归设计直接补偿阻尼系数的虚拟励磁控制器;
6)结合虚拟同步机转子运动方程得到的虚拟功角δ、逆变器调制波电压幅值E得出逆变器调制信号,实现风机并网的虚拟同步发电机控制。
2.根据权利要求1所述的方法,其特征在于,所述获取所述虚拟同步发电机输出到电网端的电压电流以及电网频率,根据所测电压电流计算虚拟同步机输出功率,包括:
1)获取虚拟同步机输出电压电流及电网频率:
通过电压电流互感器测量风机并网逆变器在三相abc坐标系下的输出电压Uabc及输出电流Iabc;通过锁相环模块测量电网频率ωg;
(2)计算虚拟同步机输出功率:
利用所述的输出电压Uabc、输出电流Iabc,按式(1)计算虚拟同步机在abc坐标系输出功率,或将电压电流转换到dq坐标系再进行计算:
Pout=UaIa+UbIb+UcIc=UdId+UqIq (1)
式中,Ud、Uq,Id、Iq分别为母线电压和母线电流在d轴和q轴分量。
3.根据权利要求2所述的方法,其特征在于,所述虚拟同步机的虚拟励磁绕组电磁动态方程,包括:
1)采用与同步机励磁绕组电磁动态方程中参数相对应的方式,建立虚拟同步发电机励磁绕组电磁动态方程为:
其中,K′d是励磁绕组的时间常数(s);Vset[p.u.]为对应于系统稳态运行时的给定励磁电压;uf为对应于励磁电压的调节量;Vq是逆变器出口电压,对应于发电机空载感应电动势;是暂态电势,Vq与V′q存在关系为:
Vq=V′q+(xvir-x′vir)Id (7)
其中,xvir表示虚拟定子电抗,用虚拟阻抗方法进行模拟;x′vir表示虚拟暂态同步电抗,同样用虚拟阻抗方法模拟;Id为母线电流d轴分量,表示为:
其中,U为无限大系统的母线电压,看作常数,x′d∑=x′vir+xl表示为虚拟定子暂态电抗与线路电抗之和。
4.根据权利要求3所述的方法,其特征在于,建立具有励磁控制的虚拟同步发电机单机无穷大系统数学模型并定义其平衡点,包括:
1)根据虚拟同步机虚拟励磁绕组电磁动态方程,将(7)代入(8),得到Vq与V′q之间有如下关系成立:
2)已知电动势、电压和电流的关系为:
将(10)代入(1)得到虚拟同步机送入无穷大系统的有功功率表达式为:
3)根据式(1)和式(2),最终建立具有励磁控制的虚拟同步发电机单机无穷大系统模型(12)如下:
4)根据控制理论中对系统平衡点的定义,即系统的平衡点是指满足f(x)≡0的点,定义系统模型(12)的状态变量为δ,ωr,V′q,并且其平衡点用(δs,ωrs,V′qs)表示,根据该定义,系统(12)的平衡点应满足如下条件:
当给定逆变器预设功率Pset和虚拟励磁电压Vset时,求解(13)非线性代数方程得到系统平衡点。
5.根据权利要求4所述的方法,其特征在于,所述直接补偿阻尼系数的虚拟励磁控制器的递归设计方案,具体包括:
1)考虑化简后的系统模型(12)中的子系统(δ,ωr),并将暂态电势V′q表示为V′q=V′qs+ΔV′q,则定义ΔV′q=V′q-V′qs,子系统(δ,ωr)写为:
其中,Pe0=a1V′qssinδ-a2sin2δ、M=2H;
2)将ΔV′q看作该子系统的虚拟控制律,通过调节逆变器虚拟励磁使q轴暂态电势的变化ΔV′q满足a1ΔV′qsinδ=M·kDωr,由于存在调节误差,定义调节误差为:
e=M·kDωr-a1ΔV′qsinδ (15)
其中,kD为自定义变量且kD>0,MkD理想补偿逆变器虚拟阻尼系数,由此得到:
对上式求导得到:
选取控制律r为:
r=-ωsωr-kee (19)
由此,得到直接改善虚拟同步机阻尼系数的励磁控制律u为:
6.根据权利要求5所述的方法,其特征在于,所述的步骤6)是根据得到的控制风机并网逆变器的PWM调制波信号相角δ;再将式(21)代入式(9)(10),通过递归计算得逆变器出口电压参考值Vq,逆变器调制波电压幅值E可由(22)求出:
Vq=Vdc·E (22)
从而得到逆变器调制波电压幅值E,即得到PWM调制波信号E∠δ,作用于风机并网逆变器,实现风机并网的虚拟同步发电机控制。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810116498.6A CN108270241B (zh) | 2018-02-06 | 2018-02-06 | 风机并网逆变器虚拟同步发电机的控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810116498.6A CN108270241B (zh) | 2018-02-06 | 2018-02-06 | 风机并网逆变器虚拟同步发电机的控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108270241A CN108270241A (zh) | 2018-07-10 |
CN108270241B true CN108270241B (zh) | 2021-06-01 |
Family
ID=62773552
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810116498.6A Active CN108270241B (zh) | 2018-02-06 | 2018-02-06 | 风机并网逆变器虚拟同步发电机的控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108270241B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109149604B (zh) * | 2018-09-19 | 2021-11-09 | 湖南工业大学 | 一种保证电力系统稳定性的采样周期优化方法 |
CN110611326A (zh) * | 2019-08-23 | 2019-12-24 | 中国电力科学研究院有限公司 | 一种并网逆变器的下垂控制方法及系统 |
CN110429655B (zh) * | 2019-09-06 | 2023-06-27 | 国网辽宁省电力有限公司 | 基于同步机三阶模型的储能单元主动支撑控制方法及系统 |
CN111313401B (zh) * | 2019-12-06 | 2023-08-15 | 国网天津市电力公司电力科学研究院 | 一种基于哈密顿控制原理的非车载充放电控制系统及方法 |
CN111564850B (zh) * | 2020-06-11 | 2023-11-10 | 王业勤 | 一种基于有界性pid控制的虚拟同步发电机型逆变器 |
CN113328457A (zh) * | 2020-12-10 | 2021-08-31 | 国网山东省电力公司德州供电公司 | 一种提高并网友好性的风电机组虚拟控制方法 |
CN112751346B (zh) * | 2020-12-30 | 2023-02-28 | 郑州轻工业大学 | 一种基于虚拟阻抗的dfig-pss控制器设计方法 |
CN113193793A (zh) * | 2021-03-18 | 2021-07-30 | 武汉科力源电气有限公司 | 基于vsg技术的同步电机背靠背调速方法及相关设备 |
CN114285100A (zh) * | 2021-12-29 | 2022-04-05 | 华北电力大学 | 提高组网型并网逆变器暂态稳定性的控制方法及装置 |
CN115882468B (zh) * | 2022-10-24 | 2024-01-23 | 国网湖北省电力有限公司电力科学研究院 | 一种基于端口能量重塑的虚拟同步控制方法 |
CN117200360B (zh) * | 2023-09-13 | 2024-03-12 | 山东大学 | 自同步电压源型双馈风电机组高低电压穿越控制方法及系统 |
CN117394432A (zh) * | 2023-10-24 | 2024-01-12 | 兰州理工大学 | 一种虚拟同步发电机的多参数自适应协同控制方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9166495B2 (en) * | 2011-06-24 | 2015-10-20 | Abb Technology Ltd. | Wind power converter |
CN107465212A (zh) * | 2017-09-25 | 2017-12-12 | 江苏大学 | 基于虚拟同步发电机技术的风电机组微网运行虚拟惯性控制系统及方法 |
-
2018
- 2018-02-06 CN CN201810116498.6A patent/CN108270241B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9166495B2 (en) * | 2011-06-24 | 2015-10-20 | Abb Technology Ltd. | Wind power converter |
CN107465212A (zh) * | 2017-09-25 | 2017-12-12 | 江苏大学 | 基于虚拟同步发电机技术的风电机组微网运行虚拟惯性控制系统及方法 |
Non-Patent Citations (2)
Title |
---|
"基于虚拟同步电机的微网逆变器控制方法研究";李国进等;《自动化与仪器仪表》;20171125(第11期);第138-140页 * |
"直接补偿发电机机械阻尼系数的新型非线性励磁控制";孙元章等;《第十一届全国电工数学学术年会论文集》;20080828;第323-336页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108270241A (zh) | 2018-07-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108270241B (zh) | 风机并网逆变器虚拟同步发电机的控制方法 | |
CN108493984B (zh) | 适用于光伏并网系统的虚拟同步发电机控制方法 | |
Guan et al. | Frequency stability of hierarchically controlled hybrid photovoltaic-battery-hydropower microgrids | |
Wang et al. | Modelling of a PMSG wind turbine with autonomous control | |
Abolhassani et al. | Integrated doubly fed electric alternator/active filter (IDEA), a viable power quality solution, for wind energy conversion systems | |
CN102074967B (zh) | 一种具有并网特性的储能型风电场控制方法 | |
Kassem et al. | Robust control of an isolated hybrid wind–diesel power system using linear quadratic Gaussian approach | |
Zheng et al. | Integral-type terminal sliding-mode control for grid-side converter in wind energy conversion systems | |
CN108199396B (zh) | 储能逆变器虚拟励磁闭环控制系统及其设计方法 | |
Mensou et al. | Performance of a vector control for DFIG driven by wind turbine: real time simulation using DS1104 controller board | |
CN107565599B (zh) | 一种基于vsg的风电并网半实物仿真系统 | |
Bharti et al. | Controller design for DFIG driven by variable speed wind turbine using static output feedback technique | |
CN105896614A (zh) | 一种光伏逆变器稳态电压平衡控制方法及系统 | |
WO2021110171A1 (zh) | 一种基于p-u下垂特性的虚拟直流电机控制方法 | |
Bouchiba et al. | Implementation and comparative study of control strategies for an isolated DFIG based WECS | |
CN111049178A (zh) | 一种直驱永磁风电机组经vsc-hvdc并网稳定控制分析方法 | |
Azongha et al. | Grid integration studies of a switched reluctance generator for future hardware-in-the-loop experiments | |
Teshager et al. | Direct power control strategy based variable speed pumped storage system for the reduction of the wind power fluctuation impact on the grid stability | |
Kenne et al. | Adaptive PI control strategy for a self-excited induction generator driven by a variable speed wind turbine | |
CN114123244A (zh) | 一种考虑风-储-直联合调频的电网频率特征计算方法 | |
Ntuli et al. | A comparative study of DFIG-based wind energy conversion systems with Fuzzy-PID and Fuzzy-logic controllers | |
Cao et al. | Impedance-besed analysis on subsynchronous oscillation mechanism caused of VSC-HVDC with double-fed induction generator | |
Raina et al. | Variable speed wind energy conversion using synchronous machine | |
Charles et al. | Field oriented control of Doubly Fed Induction Generator in wind power system | |
CN113328457A (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 |