CN108493984B - 适用于光伏并网系统的虚拟同步发电机控制方法 - Google Patents
适用于光伏并网系统的虚拟同步发电机控制方法 Download PDFInfo
- Publication number
- CN108493984B CN108493984B CN201810116489.7A CN201810116489A CN108493984B CN 108493984 B CN108493984 B CN 108493984B CN 201810116489 A CN201810116489 A CN 201810116489A CN 108493984 B CN108493984 B CN 108493984B
- Authority
- CN
- China
- Prior art keywords
- virtual
- voltage
- inverter
- synchronous generator
- control
- 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 111
- 238000000034 method Methods 0.000 title claims abstract description 42
- 230000001052 transient effect Effects 0.000 claims abstract description 20
- 230000005284 excitation Effects 0.000 claims description 54
- 238000013178 mathematical model Methods 0.000 claims description 18
- 238000004804 winding Methods 0.000 claims description 16
- 238000013016 damping Methods 0.000 claims description 15
- 230000014509 gene expression Effects 0.000 claims description 14
- 238000013461 design Methods 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 5
- 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 description 3
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 10
- 238000004088 simulation Methods 0.000 abstract description 5
- 238000010248 power generation Methods 0.000 description 10
- 238000011217 control strategy Methods 0.000 description 5
- 230000001276 controlling effect Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- H02J3/383—
-
- 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/24—Arrangements for preventing or reducing oscillations of power in networks
-
- 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/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
一种适用于光伏并网系统的虚拟同步发电机控制方法,该方法是控制生成PWM脉冲所需的调制波,调制波与载波作比较生成PWM脉冲信号,用于控制三相桥DC/AC逆变器开关管的导通和关断,实现光伏并网的VSG控制。实验仿真表明,本发明方法在光伏并网中能对DC/AC逆变器施加有效控制,对于维持系统功率平衡、电压和频率稳定具有较好的效果。改善了光伏并网系统暂态稳定性和动态品质。
Description
技术领域
本发明涉及分布式发电并网控制,特别是一种适用于光伏并网系统的虚拟同步发电机(简称为VSG)控制方法。
背景技术
能源是维系人类生存和社会发展的重要因素,在发挥不可或缺的积极作用的同时,也引发了能源短缺和化石能源污染等问题。为缓解这些问题,近年来针对太阳能、风能等可再生能源的利用形式逐渐得到各国重视而得到大力发展,其中太阳能光伏发电占了相当大的比重。截至2016年,全球累计的光伏发电装机容量已经超过300GW,而我国自2013年起已连续四年居于光伏装机总量首位——截至2017 年上半年,我国并网光伏装机总量已达到101.82GW,其中分布式光伏发电17.43GW,占比17%,未来这一比例仍将持续增大。
分布式发电具有诸多优点,如投资小、无污染、可靠性高以及发电方式灵活等,同时也存在不易调控、具有随机性和波动性等问题。分布式电源在并网过程中往往需要通过电力电子变换器(并网DC/AC逆变器)来输送电能,相比于运行成熟的传统发电设备(如同步发电机),电力电子变换器具有响应迅速,控制灵活等优势,但是自身不具备同步发电机所固有的旋转惯性和阻尼分量,当系统中接入的分布式电源达到一定规模时,必将对电力系统的安全稳定运行造成威胁。在传统的发电模式下,同步发电机凭借其所固有的旋转惯性和阻尼特性,能够为电网提供必要的电压和频率支撑,表现出良好的接入特性。受此启发,有学者借鉴同步发电机的机械方程和电磁方程提出虚拟同步发电机(virtualsynchronous generator,VSG)控制策略,使电力电子并网DC/AC逆变器能模拟出同步发电机的外特性,可以实现分布式电源的“友好型”接入,提高电力系统稳定性和优化电能质量,并能方便地将一些传统电网的运行控制策略移植到含分布式电源的电网中。
在传统电力系统中,针对发电机组的励磁控制是提高电力系统稳定性的一大有效手段,可以将发电机端电压的调节精度维持在一个给定的范围内,一般采用基于发电机端电压偏差进行的比例积分微分调节,简称PID调节方式。受此启发,基于VSG控制的DC/AC逆变器也多采用PID调节手段,虽然能够保证端电压的调节精度,但是难以有效改善电力系统的稳定性和故障后系统的动态品质。在随后的励磁调节器设计方面,新的控制律利用微分几何方法,对非线性系统采用精确的反馈线性化方法,通过非线性反馈将非线性电力系统转化为线性系统,然后按照线性系统的设计理论进行励磁调节的设计。这种方法对模型的精确度要求较高,在系统结构或参数不确定等情况下,难以保证控制效果。
针对分布式逆变电源,传统的控制策略是基于旋转坐标系解耦的电流型控制策略,在稳态并网条件下能够实现有功功率和无功功率的解耦控制,但在系统故障等暂态情况下,采用此种控制算法时,DC/AC逆变器的动态特性较差。此外,传统的解耦控制一般将系统近似为线性来进行控制,而实际上分布式DC/AC逆变器的数学模型呈现非线性特质,这样的近似处理无疑会对控制效果造成影响。
发明内容
为应对上述传统方法的不足,本发明提供一种适用于光伏并网系统的虚拟同步发电机控制方法,该方法在光伏并网中能对DC/AC逆变器施加有效控制,对于维持系统功率平衡、电压和频率稳定具有较好的效果。改善了光伏并网系统暂态稳定性和动态品质。
本发明的技术解决方案如下:
一种适用于光伏并网系统的虚拟同步发电机控制方法,其特点在于,包括下列步骤:
1)获取所述的虚拟同步发电机输出到电网端的电压、电流以及电网频率,根据所测电压、电流计算虚拟同步发电机的输出功率;
2)根据虚拟同步电机机的转子运动特性方程,结合下垂控制,建立有功频率控制模型,得到虚拟转子角速度及调制波信号相角;
3)建立采用励磁控制的虚拟同步发电机单机无穷大系统数学模型,并定义该系统运行平衡点;
4)基于上述数学模型和平衡点,采用递归设计实现直接补偿阻尼系数的虚拟励磁控制,获取调制波信号电压幅值;
5)结合所述的调制波信号电压幅值和调制波信号相角构建调制波信号,生成 PWM脉冲作用于系统中DC/AC逆变器以实现光伏并网的VSG控制。
所述的获取所述虚拟同步发电机输出到电网端的电压电流以及电网频率,根据所测电压电流计算虚拟同步发电机输出功率,具体步骤包括:
1)获取虚拟同步发电机输出电压、电流及电网频率:
通过互感器采集虚拟同步发电机在三相abc坐标系下的输出侧的电压Uabc及输出电流Iabc;电网频率通过锁相环(PLL)模块测量电网侧公共母线角速度ωg;
2)计算虚拟同步发电机输出功率:
利用所测输出侧的电压Uabc及输出电流Iabc,按式(1)计算虚拟同步发电机输出功率:
Pout=UaIa+UbIb+UcIc (1)
或先将电压电流转换到dq坐标系下再按式(2)求功率:
式中,Ud、Uq,Id、Iq分别为母线电压和母线电流在d轴和q轴分量。
所述的根据虚拟同步发电机的转子运动特性方程,结合下垂控制,建立有功频率控制模型,得到虚拟转子角速度及调制波信号相角,具体包括下列步骤:
1)模拟同步机,建立虚拟同步发电机的转子运动特性方程为:
式中,H为虚拟惯性时间常数,与同步机转动惯量J相对应;Pin、Pout为DC/AC逆变器的输入功率和输出功率,类比传统同步机的机械功率和电磁功率;ω、ωg为 DC/AC逆变器的虚拟转子角速度和电网侧公共母线角速度;Kd为阻尼系数;
再考虑有功-频率下垂控制环节,即
式中,Pset为虚拟同步控制下DC/AC逆变器的参考有功输入(可为上层调度指令);Dp为下垂控制系数;ωref为角速度参考值(一般即为50Hz,标幺值为1);
2)基于上式(3)(4),建立有功频率控制模型,得到DC/AC逆变器的虚拟转子角速度及调制波信号相角为:
所述的建立采用励磁控制的虚拟同步发电机单机无穷大系统数学模型,并定义该系统运行平衡点,具体包括:
1)基于上述转子运动方程,采用标幺值表示,并做部分简化,可得:
式中,ωr[p.u.]为虚拟转子角速度与电网侧同步角速度的偏差;ωs=2πf0为基准角速度,f0=50Hz;
2)采用与同步机励磁绕组电磁动态方程中参数相对应的方式,建立虚拟同步发电机励磁绕组电磁动态方程为:
式中,K′d是励磁绕组的时间常数(s);Vset[p.u.]为对应于系统稳态运行时给定励磁电压;uf为对应于励磁电压的调节量;Vq是所述的DC/AC逆变器出口电压,对应于发电机空载感应电动势;是暂态电势;
根据定义Vq与V′q存在关系为:
Vq=V′q+(xvir-x′vir)Id (9)
式中,xvir表示虚拟定子电抗,x′vir表示虚拟暂态同步电抗,均可用虚拟阻抗方法模拟;Id为母线电流d轴分量,可表示为:
其中,U为无限大系统的母线电压,可看作常数;x′d∑=x′vir+xl为虚拟定子暂态电抗与线路电抗之和;
3)根据虚拟同步发电机虚拟励磁绕组电磁动态方程,将式(9)代入(10)得到Vq与V′q之间有如下关系成立:
4)根据虚拟同步发电机电磁功率方程,可得到虚拟同步发电机有功功率表达式:
已知电动势、电压和电流的关系为
将(12)代入(2)得到虚拟同步发电机输出到电网的有功功率表达式为:
5)根据虚拟同步发电机转子运动方程(7)和虚拟励磁绕组电磁动态方程(8),最终建立具有励磁控制的虚拟同步发电机单机无穷大系统数学模型为:
6)定义系统运行平衡点,将上述模型简写为:
根据控制理论中对系统平衡点的定义,即系统的平衡点是指满足 f(x)≡0的点,定义系统(15)的状态变量为δ,ωr,V′q,并且其平衡点用(δs,ωrs,V′qs)表示,根据该定义,系统(15)的平衡点应满足如下条件:
当给定DC/AC逆变器预设功率Pset和虚拟励磁电压Vset时,可通过求解(16)非线性代数方程得到系统运行平衡点。
所述的基于上述数学模型和平衡点,采用递归设计实现直接补偿阻尼系数的虚拟励磁控制,获取调制波信号电压幅值,具体包括:
1)考虑化简后的系统模型(15)中的子系统(δ,ωr),并将暂态电势V′q表示为 V′q=V′qs+ΔV′q,则定义ΔV′q=V′q-V′qs,子系统(δ,ωr)可写为:
其中,Pout0=a1V′qssinδ-a2sin 2δ,记M=2H;
2)将ΔV′q看作该子系统的虚拟控制律,希望通过调节DC/AC逆变器虚拟励磁使 q轴暂态电势的变化ΔV′q满足a1ΔV′qsinδ=M·kDωr,由于存在调节误差,定义调节误差为:
e=M·kDωr-a1ΔV′qsinδ (18)
其中,kD为自定义变量且kD>0,MkD为理想补偿DC/AC逆变器虚拟阻尼系数,由此可得:
对上式求导得到:
选取控制律r为:
r=-ωsωr-kee (22)
由此,可得到直接改善虚拟同步发电机阻尼系数的励磁控制律u为:
3)已知DC/AC逆变器PWM调制波电压幅值E与DC/AC逆变器出口电压幅值 Vq存在下列关系:
Vq=Vdc·E (25)
其中,Vdc为直流侧电压估计值。
将上述控制律(24)代入式(8)、(9),通过递归计算可得DC/AC逆变器出口电压参考值Vq,再由式(25)即得DC/AC逆变器调制波信号电压幅值E。
所述的结合上述所得调制波信号相角δ、调制波信号电压幅值E构建调制波信号,再生成PWM脉冲信号作用于所述的DC/AC逆变器,以实现光伏并网的VSG 控制,这里PWM三相调制波信号的表达式为:
考虑将上述幅值E和相角δ输入到三相程控电压源以获得上述三相调制波信号,再作用于PWM发生器产生DC/AC逆变器控制脉冲。
本发明技术效果:
经实验仿真测试,本发明生成的三相调制波信号稳定且精度较高,进一步产生的PWM控制脉冲能有效控制DC/AC逆变器,保证了光伏并网的VSG控制效果。改善了光伏并网系统暂态稳定性和动态品质。
附图说明
图1是适用于本发明方法的光伏并网系统结构图;
图2是本发明适用于光伏并网系统的虚拟同步发电机控制方法流程图;
图3是虚拟同步控制下DC/AC逆变器对应传统同步机示意图;
图4是本发明控制方法示意图
图5实验仿真频率控制效果图
图6实验仿真有功控制效果图
具体实施方式
下面结合附图对本发明作进一步说明,但不应以此限制本发明的保护范围。
本发明提供一种适用于光伏并网系统的虚拟同步发电机控制方法,图1为适用于本发明方法的光伏并网系统结构图。图中光伏电池板输出的直流电先经DC/DC 变换器调压,再经DC/AC逆变器输出,通过LC滤波环节后输送到电网。控制DC/AC 逆变器内电力电子器件导通和关断的PWM脉冲由本发明控制方法得到的调制波信号E∠δ产生,以实现光伏并网的VSG控制。
图2是本发明适用于光伏并网系统的虚拟同步发电机控制方法流程图,本发明适用于光伏并网系统的虚拟同步发电机控制方法,包括下列步骤:
1)获取所述的虚拟同步发电机输出到电网端的电压、电流以及电网频率,根据所测电压、电流计算虚拟同步发电机的输出功率;
2)根据虚拟同步电机机的转子运动特性方程,结合下垂控制,建立有功频率控制模型,得到虚拟转子角速度及调制波信号相角;
3)建立采用励磁控制的虚拟同步发电机单机无穷大系统数学模型,并定义该系统运行平衡点;
4)基于上述数学模型和平衡点,采用递归设计实现直接补偿阻尼系数的虚拟励磁控制,获取调制波信号电压幅值;
5)结合所述的调制波信号电压幅值和调制波信号相角构建调制波信号,生成 PWM脉冲作用于系统中DC/AC逆变器以实现光伏并网的VSG控制。
所述的获取所述虚拟同步发电机输出到电网端的电压电流以及电网频率,根据所测电压电流计算虚拟同步发电机输出功率,具体步骤包括:
1)获取虚拟同步发电机输出电压、电流及电网频率:
通过互感器采集虚拟同步发电机在三相abc坐标系下的输出侧的电压Uabc及输出电流Iabc;电网频率通过锁相环(PLL)模块测量电网侧公共母线角速度ωg;
2)计算虚拟同步发电机输出功率:
利用所测输出侧的电压Uabc及输出电流Iabc,按式(1)计算虚拟同步发电机输出功率:
Pout=UaIa+UbIb+UcIc (1)
或先将电压电流转换到dq坐标系下再按式(2)求功率:
式中,Ud、Uq,Id、Iq分别为母线电压和母线电流在d轴和q轴分量。
所述的根据虚拟同步发电机的转子运动特性方程,结合下垂控制,建立有功频率控制模型,得到虚拟转子角速度及调制波信号相角,具体包括下列步骤:
1)模拟同步机,建立虚拟同步发电机的转子运动特性方程为:
式中,H为虚拟惯性时间常数,与同步机转动惯量J相对应;Pin、Pout为DC/AC逆变器的输入功率和输出功率,类比传统同步机的机械功率和电磁功率;ω、ωg为 DC/AC逆变器的虚拟转子角速度和电网侧公共母线角速度;Kd为阻尼系数;
再考虑有功-频率下垂控制环节,即
式中,Pset为虚拟同步控制下DC/AC逆变器的参考有功输入(可为上层调度指令);Dp为下垂控制系数;ωref为角速度参考值(一般即为50Hz,标幺值为1);
2)基于上式(3)(4),建立有功频率控制模型,得到DC/AC逆变器的虚拟转子角速度及调制波信号相角为:
所述的建立采用励磁控制的虚拟同步发电机单机无穷大系统数学模型,并定义该系统运行平衡点,具体包括:
1)基于上述转子运动方程,采用标幺值表示,并做部分简化,可得:
式中,ωr[p.u.]为虚拟转子角速度与电网侧同步角速度的偏差;ωs=2πf0为基准角速度,f0=50Hz;
2)采用与同步机励磁绕组电磁动态方程中参数相对应的方式,建立虚拟同步发电机励磁绕组电磁动态方程为:
式中,K′d是励磁绕组的时间常数(s);Vset[p.u.]为对应于系统稳态运行时给定励磁电压;uf为对应于励磁电压的调节量;Vq是所述的DC/AC逆变器出口电压,对应于发电机空载感应电动势;是暂态电势;
根据定义Vq与V′q存在关系为:
Vq=V′q+(xvir-x′vir)Id (9)
式中,xvir表示虚拟定子电抗,x′vir表示虚拟暂态同步电抗,均可用虚拟阻抗方法模拟;Id为母线电流d轴分量,可表示为:
其中,U为无限大系统的母线电压,可看作常数;x′d∑=x′vir+xl为虚拟定子暂态电抗与线路电抗之和;
3)根据虚拟同步发电机虚拟励磁绕组电磁动态方程,将式(9)代入(10)得到Vq与V′q之间有如下关系成立:
4)根据虚拟同步发电机电磁功率方程,可得到虚拟同步发电机有功功率表达式:
已知电动势、电压和电流的关系为
将(12)代入(2)得到虚拟同步发电机输出到电网的有功功率表达式为:
5)根据虚拟同步发电机转子运动方程(7)和虚拟励磁绕组电磁动态方程(8),最终建立具有励磁控制的虚拟同步发电机单机无穷大系统数学模型为:
6)定义系统运行平衡点,将上述模型简写为:
根据控制理论中对系统平衡点的定义,即系统的平衡点是指满足 f(x)≡0的点,定义系统(15)的状态变量为δ,ωr,V′q,并且其平衡点用(δs,ωrs, V′qs)表示,根据该定义,系统(15)的平衡点应满足如下条件:
当给定DC/AC逆变器预设功率Pset和虚拟励磁电压Vset时,通过求解(16)非线性代数方程得到系统运行平衡点。
基于上述数学模型和平衡点,采用递归设计实现直接补偿阻尼系数的虚拟励磁控制,获取调制波信号电压幅值,具体包括:
1)考虑化简后的系统模型(15)中的子系统(δ,ωr),并将暂态电势V′q表示为 V′q=V′qs+ΔV′q,则定义ΔV′q=V′q-V′qs,子系统(δ,ωr)可写为:
其中,Pout0=a1V′qssinδ-a2sin 2δ,记M=2H;
2)将ΔV′q看作该子系统的虚拟控制律,希望通过调节DC/AC逆变器虚拟励磁使 q轴暂态电势的变化ΔV′q满足a1ΔV′qsinδ=M·kDωr,由于存在调节误差,定义调节误差为:
e=M·kDωr-a1ΔV′qsinδ (18)
其中,kD为自定义变量且kD>0,MkD为理想补偿DC/AC逆变器虚拟阻尼系数,由此可得:
对上式求导得到:
选取控制律r为:
r=-ωsωr-kee (22)
由此,可得到直接改善虚拟同步发电机阻尼系数的励磁控制律u为:
3)已知DC/AC逆变器PWM调制波电压幅值E与DC/AC逆变器出口电压幅值 Vq存在下列关系:
Vq=Vdc·E (25)
其中,Vdc为直流侧电压估计值,将上述控制律(24)代入式(8)、(9),通过递归计算可得DC/AC逆变器出口电压参考值Vq,再由式(25)即得DC/AC逆变器调制波信号电压幅值E。
所述的结合上述所得调制波信号相角δ、调制波信号电压幅值E构建调制波信号,再生成PWM脉冲信号作用于所述的DC/AC逆变器,以实现光伏并网的VSG 控制,这里PWM三相调制波信号的表达式为:
考虑将上述幅值E和相角δ输入到三相程控电压源以获得上述三相调制波信号,再作用于PWM发生器产生DC/AC逆变器控制脉冲。
简略地说,本发明方法分为两大部分,分别求取调制波信号的相角和幅值,第一步获取DC/AC逆变器输出电压、电流以及电网频率并计算虚拟同步发电机输出功率,在此基础上第二步建立有功频率控制模型,推导虚拟转子角速度计算公式,进而求得调制波信号相角;第三步建立结合励磁控制的VSG系统数学模型,定义并求取系统运行平衡点,在此基础上第四步采用递归设计实现虚拟励磁控制,进而获取调制波信号电压幅值;最后一步结合所得幅值E和相角δ构建调制波信号E∠δ,生成PWM脉冲作用于光伏并网DC/AC逆变器。
具体地,步骤二中,建立有功频率控制模型,除了依照虚拟同步发电机的转子运动特性方程外,还考虑加入有功频率下垂控制环节,可以使分布式逆变电源具备一定的调频能力,从而在系统负荷变化的情况下可以给系统提供一定的频率支持,以提高系统的稳定性能;
步骤三中,参照图3,采用与同步发电机参数对应的方式,建立虚拟同步发电机控制下的DC/AC逆变器转子运动方程,励磁绕组电磁动态方程以及虚拟同步发电机电磁动态方程,最终得到虚拟同步发电机控制下的DC/AC逆变器系统数学模型表达式;
步骤四中,通过构造李雅普诺夫函数,采用递归设计实现直接补偿阻尼系数的虚拟励磁控制,充分利用了虚拟同步发电机本身的暂态能量函数;最后一步,构建的调制波与载波作比较生成PWM脉冲信号,用于控制三相桥DC/AC逆变器开关管的导通和关断,实现光伏并网的VSG控制。
图4为本发明控制方法示意图,其中上半部分对应图2中的第一步、第二步, PLL锁相环模块获取电网频率ωg,Pset、Pout作为参数输入,经过由下垂控制环节和转子运动特性方程组合而成的有功频率控制模型,得到虚拟转子角速度ω及调制波信号相角δ。而下半部分对应虚拟励磁控制模型,基于虚拟励磁控制器得到的控制律u,再通过一系列递归计算可得到用于DC/AC逆变器控制的调制波信号幅值E。最后结合上下两部分所得到的相角δ和幅值E,构建调制波信号,输入到PWM发生器中生成PWM脉冲。上下两部分相互沟通联系,构成一个整体控制策略,也即本发明的控制方法。
图5、图6为硬件在环(HIL)实验仿真结果。其中,图5所示为DC/AC逆变器虚拟转子角速度ω(p.u.),可见在初始化完成后,其一直稳定在标幺值1附近。当有功需求变动时,将自身波动限制在0.5%以内,间接反映了本控制方法对系统频率的良好控制。图6所示,上下分别为虚拟同步控制下DC/AC逆变器的有功输入参考值Pset(p.u.)和有功输出Pout(p.u.)。由图可见,当有功需求(即参考值)变动时,DC/AC逆变器实际有功输出也能快速反应跟进,反映了较好的功率跟随效果。通过实验仿真,说明本发明提出的虚拟同步控制方法在光伏并网中能对DC/AC逆变器施加有效控制,对于维持系统功率平衡、电压和频率稳定具有较好的效果。
Claims (5)
1.一种适用于光伏并网系统的虚拟同步发电机控制方法,其特征在于,包括下列步骤:
1)获取所述的虚拟同步发电机输出到电网端的电压、电流以及电网频率,根据所测电压、电流计算虚拟同步发电机的输出功率;
2)根据虚拟同步电机机的转子运动特性方程,结合下垂控制,建立有功频率控制模型,得到虚拟转子角速度及调制波信号相角,具体包括下列步骤:
2.1)建立虚拟同步发电机的转子运动特性方程为:
式中,H为虚拟惯性时间常数,与同步机转动惯量J相对应;Pin、Pout为DC/AC逆变器的输入功率和输出功率,类比传统同步机的机械功率和电磁功率;ω、ωg为DC/AC逆变器的虚拟转子角速度和电网侧公共母线角速度;Kd为阻尼系数;
再考虑有功-频率下垂控制环节,即
式中,Pset为虚拟同步控制下DC/AC逆变器的参考有功输入或上层调度指令;Dp为下垂控制系数;ωref为角速度参考值,即为50Hz,标幺值为1;
2.2)基于上式(3)、(4),建立有功频率控制模型,得到DC/AC逆变器的虚拟转子角速度及调制波信号相角为:
2.3)建立采用励磁控制的虚拟同步发电机单机无穷大系统数学模型,并定义该系统运行平衡点;
2.4)基于上述数学模型和平衡点,采用递归设计实现直接补偿阻尼系数的虚拟励磁控制,获取调制波信号的电压幅值;
2.5)结合所述的调制波信号电压幅值和调制波信号相角构建调制波信号,生成PWM脉冲作用于系统中DC/AC逆变器,以实现光伏并网虚拟同步发电机控制。
2.根据权利要求1所述的方法,其特征在于,所述的获取所述虚拟同步发电机输出到电网端的电压电流以及电网频率,根据所测电压电流计算虚拟同步发电机输出功率,包括:
1)获取虚拟同步发电机输出电压、电流及电网频率:
通过互感器采集虚拟同步发电机在三相abc坐标系下的输出侧的电压Uabc及输出电流Iabc;电网频率通过锁相环(PLL)模块测量电网侧公共母线角速度ωg;
2)计算虚拟同步发电机输出功率:
利用所述的电压Uabc及输出电流Iabc,按式(1)计算虚拟同步发电机输出功率:
Pout=UaIa+UbIb+UcIc (1)
或先将电压电流转换到dq坐标系下再按式(2)求功率:
式中,Ud、Uq,Id、Iq分别为母线电压和母线电流在d轴和q轴分量。
3.根据权利要求2所述的方法,其特征在于,所述的建立采用励磁控制的虚拟同步发电机单机无穷大系统数学模型,并定义该系统运行平衡点,具体包括:
1)基于上述转子运动特性方程,采用标幺值表示,并做部分简化,得:
式中,ωr[p.u.]为虚拟转子角速度与电网侧同步角速度的偏差;ωs=2πf0为基准角速度,f0=50Hz;
2)采用与同步机励磁绕组电磁动态方程中参数相对应的方式,建立虚拟同步发电机励磁绕组电磁动态方程为:
式中,K′d是励磁绕组的时间常数(s);Vset[p.u.]为对应于系统稳态运行时给定励磁电压;uf为对应于励磁电压的调节量;Vq是所述的DC/AC逆变器出口电压,对应于发电机空载感应电动势;是暂态电势;
根据定义Vq与Vq′存在关系为:
Vq=V′q+(xvir-x′vir)Id (9)
式中,xvir表示虚拟定子电抗,x′vir表示虚拟暂态同步电抗,均用虚拟阻抗方法模拟;Id为母线电流d轴分量,表示为:
其中,U为无限大系统的母线电压,看作常数;x′d∑=x′vir+xl为虚拟定子暂态电抗与线路电抗之和;
3)根据虚拟同步发电机虚拟励磁绕组电磁动态方程,将式(9)代入(10)得到Vq与Vq′之间有如下关系成立:
4)根据虚拟同步发电机电磁功率方程,得到虚拟同步发电机有功功率表达式:
已知电动势、电压和电流的关系为:
将(12)代入(2)得到虚拟同步发电机输出到电网的有功功率表达式为:
5)根据转子运动特性方程(7)和虚拟励磁绕组电磁动态方程(8),最终建立具有励磁控制的虚拟同步发电机单机无穷大系统数学模型为:
6)定义系统运行平衡点,将上述模型简写为:
根据控制理论中对系统平衡点的定义,即系统的平衡点是指满足f(x)≡0的点,定义系统(15)的状态变量为δ,ωr,V′q,并且其平衡点用(δs,ωrs,V′qs)表示,根据该定义,系统(15)的平衡点应满足如下条件:
当给定DC/AC逆变器预设功率Pset和虚拟励磁电压Vset时,通过求解(16)非线性代数方程得到系统运行平衡点。
4.根据权利要求3所述的方法,其特征在于,所述的基于上述数学模型和平衡点,采用递归设计实现直接补偿阻尼系数的虚拟励磁控制,获取调制波信号电压幅值,具体包括:
1)考虑化简后的系统模型(15)中的子系统(δ,ωr),并将暂态电势V′q表示为V′q=V′qs+ΔV′q,则定义ΔV′q=V′q-V′qs,子系统(δ,ωr)写为:
其中,Pout0=a1V′qssinδ-a2sin2δ,记M=2H;
2)将△V′q看作该子系统的虚拟控制律,通过调节DC/AC逆变器虚拟励磁使q轴暂态电势的变化ΔV′q满足a1ΔV′qsinδ=M·kDωr,由于存在调节误差,定义调节误差为:
e=M·kDωr-a1ΔV′qsinδ (18)
其中,kD为自定义变量且kD>0,MkD为理想补偿DC/AC逆变器虚拟阻尼系数,由此得:
对上式求导得到:
选取控制律r为:
r=-ωsωr-kee (22)
由此,得到直接改善虚拟同步发电机阻尼系数的励磁控制律u为:
3)已知DC/AC逆变器PWM调制波电压幅值E与DC/AC逆变器出口电压幅值Vq存在下列关系:
Vq=Vdc·E (25)
其中,Vdc为直流侧电压估计值,将上述控制律(24)代入式(8)、(9),通过递归计算得DC/AC逆变器出口电压参考值Vq,再由式(25)即得DC/AC逆变器调制波信号电压幅值E。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810116489.7A CN108493984B (zh) | 2018-02-06 | 2018-02-06 | 适用于光伏并网系统的虚拟同步发电机控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810116489.7A CN108493984B (zh) | 2018-02-06 | 2018-02-06 | 适用于光伏并网系统的虚拟同步发电机控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108493984A CN108493984A (zh) | 2018-09-04 |
CN108493984B true CN108493984B (zh) | 2021-10-08 |
Family
ID=63344279
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810116489.7A Active CN108493984B (zh) | 2018-02-06 | 2018-02-06 | 适用于光伏并网系统的虚拟同步发电机控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108493984B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109301858A (zh) * | 2018-09-07 | 2019-02-01 | 南昌大学 | 一种基于自适应动态规划优化虚拟同步机光伏并网逆变器的控制方法 |
CN109449979B (zh) * | 2018-11-23 | 2020-10-09 | 燕山大学 | 一种基于虚拟同步发电机的光伏振荡平抑方法及系统 |
CN110518626A (zh) * | 2019-05-22 | 2019-11-29 | 国网辽宁省电力有限公司电力科学研究院 | 一种基于同步发电机标准三阶模型的光伏并网逆变器主动支撑控制方法 |
CN110266062B (zh) * | 2019-06-20 | 2023-04-07 | 上海交通大学 | 逆变型分布式电源双层自适应惯量控制方法及装置 |
CN110504706B (zh) * | 2019-07-22 | 2021-05-04 | 中国农业大学 | 虚拟同步发电机并网的暂态稳定性预测方法及装置 |
JP6981528B2 (ja) * | 2019-08-02 | 2021-12-15 | 東芝三菱電機産業システム株式会社 | 発電システムおよびプラント制御装置 |
EP3832829A1 (en) * | 2019-12-02 | 2021-06-09 | ABB Schweiz AG | Grid forming vector current control |
US11735928B2 (en) | 2020-04-17 | 2023-08-22 | Hitachi Energy Switzerland Ag | Multi-port grid forming control for grid interties |
CN111564850B (zh) * | 2020-06-11 | 2023-11-10 | 王业勤 | 一种基于有界性pid控制的虚拟同步发电机型逆变器 |
CN113328457A (zh) * | 2020-12-10 | 2021-08-31 | 国网山东省电力公司德州供电公司 | 一种提高并网友好性的风电机组虚拟控制方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103066876A (zh) * | 2013-01-17 | 2013-04-24 | 电子科技大学 | 一种具有同步发电机特性的逆变控制器 |
CN107147151A (zh) * | 2017-04-20 | 2017-09-08 | 华北电力大学(保定) | 一种同步发电机逆变器的并联功率分配控制方法 |
-
2018
- 2018-02-06 CN CN201810116489.7A patent/CN108493984B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103066876A (zh) * | 2013-01-17 | 2013-04-24 | 电子科技大学 | 一种具有同步发电机特性的逆变控制器 |
CN107147151A (zh) * | 2017-04-20 | 2017-09-08 | 华北电力大学(保定) | 一种同步发电机逆变器的并联功率分配控制方法 |
Non-Patent Citations (2)
Title |
---|
基于虚拟同步电机的微网逆变器控制方法研究;李国进等;《自动化与仪器仪表》;20171125(第11期);第323-336页 * |
直接补偿发电机机械阻尼系数的新型非线性励磁控制;孙元章等;《第十一届全国电工数学学术年会论文集》;20080828;第138-140页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108493984A (zh) | 2018-09-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108493984B (zh) | 适用于光伏并网系统的虚拟同步发电机控制方法 | |
CN108270241B (zh) | 风机并网逆变器虚拟同步发电机的控制方法 | |
CN110277803B (zh) | 一种储能变流器的虚拟同步发电机控制方法及控制装置 | |
CN105914778B (zh) | 基于虚拟同步发电机的微网逆变器多环路控制方法 | |
CN108199396B (zh) | 储能逆变器虚拟励磁闭环控制系统及其设计方法 | |
WO2022077847A1 (zh) | 一种用于混合微电网mmc互联变换器的虚拟同步机控制方法 | |
CN102074967B (zh) | 一种具有并网特性的储能型风电场控制方法 | |
CN106877710B (zh) | 基于虚拟同步电动机的三相pwm整流器多环路控制电路及控制方法 | |
CN110429655B (zh) | 基于同步机三阶模型的储能单元主动支撑控制方法及系统 | |
Kassem et al. | Robust control of an isolated hybrid wind–diesel power system using linear quadratic Gaussian approach | |
CN103999315B (zh) | 用于提高发电厂有功功率生产量的发电机和svc的经协调的控制方法及其控制器 | |
CN108964040B (zh) | 电网不平衡下虚拟同步发电机功率-电流协调控制方法 | |
CN110994680B (zh) | 一种自适应惯性系数的光伏虚拟同步机稳定控制方法 | |
CN104638679A (zh) | 一种采用自适应调节的虚拟同步发电机频率控制方法 | |
CN115313524B (zh) | 一种基于构网型变流器的光伏发电并网控制方法及系统 | |
CN105811825A (zh) | 基于电流补偿的虚拟同步发电机功率解耦方法 | |
WO2021110171A1 (zh) | 一种基于p-u下垂特性的虚拟直流电机控制方法 | |
CN107565599B (zh) | 一种基于vsg的风电并网半实物仿真系统 | |
Erazo-Damian et al. | Stand-alone DFIG FOC sensitivity and stability under mismatched inductances | |
CN104064083B (zh) | 一种电网模拟控制系统 | |
CN106533289B (zh) | 一种非线性电压控制方法及系统 | |
CN107887919A (zh) | 一种分布式储能逆变器协调控制方法及系统 | |
CN109066793A (zh) | 一种微网逆变器的虚拟柴油发电机组控制方法 | |
CN116316848A (zh) | 一种基于虚拟同步发电机、虚拟同步电动机的微电网中直驱风机协同控制系统及方法 | |
CN112000018B (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 |