CN117394432B - 一种虚拟同步发电机的多参数自适应协同控制方法 - Google Patents
一种虚拟同步发电机的多参数自适应协同控制方法 Download PDFInfo
- Publication number
- CN117394432B CN117394432B CN202311380654.7A CN202311380654A CN117394432B CN 117394432 B CN117394432 B CN 117394432B CN 202311380654 A CN202311380654 A CN 202311380654A CN 117394432 B CN117394432 B CN 117394432B
- Authority
- CN
- China
- Prior art keywords
- virtual
- inertia
- angular frequency
- obtaining
- damping coefficient
- 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
- 238000000034 method Methods 0.000 title claims abstract description 46
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 31
- 238000013016 damping Methods 0.000 claims abstract description 84
- 230000008859 change Effects 0.000 claims abstract description 52
- 230000003044 adaptive effect Effects 0.000 claims abstract description 30
- 230000008569 process Effects 0.000 claims abstract description 21
- 230000001052 transient effect Effects 0.000 claims abstract description 16
- 238000000638 solvent extraction Methods 0.000 claims abstract description 5
- 230000001965 increasing effect Effects 0.000 claims description 29
- 230000001133 acceleration Effects 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000012546 transfer Methods 0.000 claims description 7
- 230000004044 response Effects 0.000 description 11
- 230000007423 decrease Effects 0.000 description 6
- 238000005192 partition Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000010355 oscillation Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004134 energy conservation Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000003912 environmental pollution Methods 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 230000000087 stabilizing effect Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000013468 resource allocation Methods 0.000 description 1
Classifications
-
- 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
-
- 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/002—Flicker reduction, e.g. compensation of flicker introduced by non-linear load
-
- 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/381—Dispersed generators
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
- H02M1/00—Details of apparatus for conversion
- H02M1/0003—Details of control, feedback or regulation circuits
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
- H02M7/00—Conversion of ac power input into dc power output; Conversion of dc power input into ac power output
- H02M7/42—Conversion of dc power input into ac power output without possibility of reversal
- H02M7/44—Conversion of dc power input into ac power output without possibility of reversal by static converters
- H02M7/48—Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode
-
- 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/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- Nonlinear Science (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明公开了一种虚拟同步发电机的多参数自适应协同控制方法,包括:基于功率偏差和角频率变化率对虚拟同步发电机的暂态过程进行分区得到若干区域;基于有功功率方程得到虚拟阻抗对角频率变化率的第一影响规则,基于转子运动方程获得阻尼系数对角频率偏差的第二影响规则和转动惯量对角频率变化率的第三影响规则;基于第一影响规则、第二影响规则和第三影响规则得到各区域对应的调整规则;基于各区域对应的调整规则确定虚拟阻抗、转动惯量和阻尼系数的自适应函数;基于虚拟阻抗、转动惯量和阻尼系数的自适应函数得到各自适应函数对应的参数值。本发明利用虚拟阻抗、转动惯量和阻尼系数的灵活可调节性有效提高VSG系统的动态性能。
Description
技术领域
本发明涉及逆变器控制技术领域,特别涉及一种虚拟同步发电机的多参数自适应协同控制方法。
背景技术
经济的飞速发展,同时也带来了化石能源枯竭与环境污染的两大问题,随着双碳目标的提出,环保节能已经成为当下的热点话题,鉴于此,能源格局深化改革,分布式发电技术迅速发展,并逐渐在传统电力系统中占据一席之地。然而,分布式能源在电力系统中所占比重不断増加的同时,也给电力系统的安全稳定运行带来新的影响与挑战。因为基于电力电子逆变器接口的分布式电源不具备有利于保持系统稳定运行的旋转惯性和阻尼分量,而这些正是传统电力系统中占支配地位的同步发电机所固有的特质。故虚拟同步发电机技术(VSG)借鉴同步发电机的运行特性,使以电力电子器件为主的微电网具备了惯量和阻尼特性,不仅提高了系统稳定性,还可以更好地实现与传统电网的完美兼容。
近年来,恒功率控制、恒压/恒频控制和下垂控制等多种接口逆变器控制策略被应用到微电网运行控制中,但是往往运行模式比较单一,传统VSG控制方法其虚拟阻抗、转动惯量和阻尼系数往往是固定的,不能实时根据电力系统要求变化,彼此间配合和参数调控比较困难,不能解决高比例新能源所带来的问题。
发明内容
针对现有技术的不足之处,提供一种虚拟同步发电机的多参数自适应协同控制方法,本发明首先通过VSG暂态过程响应曲线,对VSG的暂态过程中的一个振荡周期进行分区,得到不同分区内功率偏差和角频率变化率的变化情况,然后通过对有功功率方程和转子运动方程进行理论分析得到虚拟阻抗和转动惯量、阻尼系数对角频率偏差和角频率变化率的影响规律,进一步与暂态过程分区相结合得到虚拟阻抗、转动惯量和阻尼系数的实时调整规则,依据实时调整规则设计虚拟阻抗、转动惯量和阻尼系数的自适应函数,最后确定VSG的自适应函数的参数。
本发明提供的一种虚拟同步发电机的多参数自适应协同控制方法,包括:
基于功率偏差和角频率变化率对虚拟同步发电机的暂态过程进行分区得到若干区域;
基于有功功率方程得到虚拟阻抗对所述角频率变化率的第一影响规则,基于转子运动方程获得阻尼系数对角频率偏差的第二影响规则和转动惯量对所述角频率变化率的第三影响规则;
基于所述第一影响规则、第二影响规则和第三影响规则得到若干所述区域对应的调整规则;
基于若干所述区域对应的调整规则确定虚拟阻抗、转动惯量和阻尼系数的自适应函数;
基于所述虚拟阻抗、转动惯量和阻尼系数的自适应函数得到各自适应函数对应的参数值。
可选地,基于有功功率方程得到虚拟阻抗对所述角频率变化率的第一影响规则的计算公式为:
式中,Ev为虚拟阻抗上的电压,E为逆变器的输出电压,Ev=IgZv,Zv=Rv+jωLv,Rv表示引入的虚拟电阻、Lv表示引入的虚拟电感,Zv为引入的虚拟阻抗,ω为VSG的同步转子角频率,Ig为逆变器的输出电流,Ug为逆变器的输出电压幅值,X为逆变器输出实际电抗,δ为功角。
可选地,基于转子运动方程获得阻尼系数对角频率偏差的第二影响规则的计算公式为:
式中,Δω表示角频率偏差,Tm表示机械转矩,Pm表示机械功率,Te表示电磁转矩,/>Pe为电磁功率,ω0为额定转子角频率,J为转动惯量,D为阻尼系数,ω为VSG的同步转子角频率,t表示时间。
可选地,对所述第二影响规则的计算公式进行阻尼项去除且有功功率偏差一定时得到转动惯量对所述角频率变化率的第三影响规则,所述第三影响规则的计算公式为:
式中,表示角频率变化率。
可选地,若干所述区域包括第一加速区域、第一减速区域、第二加速区域和第二减速区域;其中,各区域对应的调整规则包括:
当处于所述第一加速区域和第二加速区域时,采取减小虚拟阻抗,增大转动惯量和增加阻尼系数操作;
当处于所述第一减速区域和第二减速区域时,采取增大虚拟阻抗,减小转动惯量和增加阻尼系数操作。
可选地,构建的所述虚拟阻抗的自适应函数为:
式中,Lv0为初始虚拟电感,Δθv为虚拟阻抗角变化量,Rv0为初始虚拟电阻,Xv0为初始虚拟电抗,Zv0为初始虚拟阻抗,H表示角频率变化的阈值,ΔP表示额定功率与电磁功率的差值,Δω表示角频率偏差。
可选地,构建的所述转动惯量的自适应函数为:
式中,J0为初始虚拟惯量,k1为虚拟惯量调节系数,始终为正;m为虚拟惯量调节指数。
可选地,构建的所述阻尼系数的自适应函数为:
式中,D0为初始阻尼系数,kd为阻尼系数的调节系数。
可选地,基于所述虚拟阻抗的自适应函数得到初始虚拟电感值和虚拟电感变化量;
其中,获得所述初始虚拟电感值的过程包括:
基于VSG的二阶电压方程得到初始虚拟电感值;
获得所述虚拟电感变化量包括:
基于所述角频率偏差计算得到虚拟电感角补偿值的调整系数;
基于所述虚拟电感角补偿值的调整系数和所述角频率偏差得到虚拟阻抗角变化量;
通过改变所述虚拟阻抗角变化量得到虚拟电感的变化量。
可选地,基于所述转动惯量的自适应函数确定初始转动惯量,基于所述阻尼系数的自适应函数确定初始阻尼系数;
其中,确定所述初始转动惯量和所述初始阻尼系数的过程包括:
基于小信号模型得到有功闭环传递函数;
基于所述有功闭环传递函数得到闭环特征方程;
基于所述闭环特征方程得到转动惯量的根轨迹和阻尼系数的根轨迹;
基于所述转动惯量的根轨迹得到所述初始转动惯量,基于所述阻尼系数的根轨迹得到初始阻尼系数。
本发明具有如下技术效果:
本发明首先对VSG的暂态过程中的一个振荡周期进行分区,得到不同分区内功率偏差和角频率变化率的变化情况,然后通过对有功功率方程和转子运动方程进行理论分析得到虚拟阻抗和转动惯量、阻尼系数对角频率偏差和角频率变化率的影响规律,进一步与暂态过程分区相结合得到虚拟阻抗、转动惯量和阻尼系数的实时调节规则,设计相应的自适应函数,为了提升系统的稳定性,设计出虚拟阻抗的自适应函数,考虑到控制目的以及指数函数曲线的变化特点,本发明以指数函数为基础建立关于惯量的自适应函数,为使系统的综合动态性能达到最佳,进而建立阻尼系数的自适应函数,最后确定VSG的自适应函数的参数。本发明利用虚拟阻抗、转动惯量和阻尼系数的灵活可调节性可以有效提高VSG系统的动态性能,提高系统的频率稳定性,并且可以有效提高能源利用效率,降低碳排放,优化资源配置,满足绿色发展需求,有效改善环境污染问题,在环保节能方面有着广泛的应用前景和积极促进作用。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例中的实施流程图;
图2为本发明实施例中VSG的暂态过程响应曲线;
图3为本发明实施例中转动惯量变化下的根轨迹图;
图4为本发明实施例中阻尼系数变化下的根轨迹图;
图5为本发明实施例中虚拟惯量调节指数m变化下,J随k1和dωdt变化的三维曲线簇;
图6为本发明实施例中不同情形下有功功率变化对比图;
图7为本发明实施例中同情形下频率变化对比图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例一
如图1所示,本发明提出一种虚拟同步发电机的多参数自适应协同控制方法,步骤如下:
S1,根据功率偏差和角频率变化率对VSG的暂态过程进行分区,具体实施过程包括:
图2为VSG的暂态过程响应曲线,将其分为4个区间。
第一阶段t0-t1,又称其为加速区间,Pref-Pe>0,dω/d(t)>0,角频率ω逐渐减大,到t1时达到最大幅值,此区间内Δω>0,功角δ逐渐增大,电磁功率Pe也逐渐增大。
第二阶段t1-t2,又称其为减速区间,Pref-Pe<0,dω/d(t)<0,角频率ω逐渐减小,到t2时达到额定角频率,此区间内Δω>0,功角δ逐渐增大,电磁功率Pe也逐渐增大。
第三阶段t2-t3,又称其为加速区间,Pref-Pe<0,dω/d(t)<0,角频率ω逐渐增大,到t3时|ω|达到最大幅值,此区间内Δω<0,此过程中功角δ逐渐增大,电磁功率Pe在逐渐减小。
第四阶段t3-t4,又称其为减速区间,Pref-Pe>0,dω/d(t)>0,角频率ω逐渐减小,到t4时达到额定角频率,此区间内Δω<0,此过程中功角δ逐渐减小,电磁功率的绝对值|Pe|在逐渐增大。
S2,根据有功功率方程和转子运动方程获得虚拟阻抗和转动惯量、阻尼系数对角频率偏差和角频率变化率的影响规律,具体包括:基于有功功率方程得到虚拟电阻对角频率变化率的第一影响规则,基于转子运动方程获得阻尼系数对角频率偏差的第二影响规则和转动惯量对角频率变化率的第三影响规则。
S21,根据有功功率方程获得虚拟阻抗对角频率变化率的第一影响规律;
在VSG中线路大多是电阻性的,加入虚拟阻抗技术,使逆变器的输出阻抗趋于感性,实现功率解耦,极大地改善了无功功率的分配,减小了逆变器之间的环流。
式中:E*为参考电压,Ev为虚拟电阻上的电压,E为逆变器的输出电压,Rv、Lv分别为引入的虚拟电阻和虚拟电感,Zv为引入的虚拟阻抗,Ig为逆变器的输出电流。
VSG的输出有功功率方程为:
式中:Ug为逆变器的输出电压幅值,X为逆变器输出实际电抗,δ为功角。
将式(1)和式(2)结合,可得输出功率和虚拟阻抗电压的关系式:
若适当改变虚拟阻抗的值,则虚拟阻抗上的电压会发生相应变化,从而改变了VSG的输出有功功率,进一步会影响系统频率的波动。进一步当给定的有功功率突然下降时,适当增大虚拟阻抗的值,会使Pe加速追踪新的P′ref设定,减小了系统频率波动,缩短系统达到新的稳态的时间。
S22,基于转子运动方程获得阻尼系数对角频率偏差的第二影响规则和转动惯量对角频率变化率的第三影响规则;
令VSG的极对数为1,则VSG的转子运动方程可以表示为:
式中:J为转动惯量,D为阻尼系数,ω为VSG的同步转子角频率,ω0为额定转子角频率,Pm为机械功率,Pe为电磁功率,δ为功角。
根据式(4)和式(5)可得角频率偏差Δω与阻尼系数D之间的第二影响规则如式(6)所示:
假设Tm-Te-J(dω/dt)恒定时,当阻尼系数D增大,角频率偏差Δω会减小,理论上,增大阻尼系数D可对角频率偏差起抑制作用,但是阻尼系数D过大,会使得系统调节时的阻力过大,进而系统响应速度变慢。
忽略式(6)中的阻尼项,可得转动惯量J和角频率变化率dω/dt之间的关系:
由式(7)可知,当有功功率偏差一定时,转动惯量J与角频率变化率dω/dt成反比,进而可得到第三影响规则:
由式(8)可知,当系统频率受到扰动时,可以通过调整转动惯量J,来抑制系统频率的波动,理论上,在系统受到扰动时,转动惯量J设置越大,则角频率的变化率越小,对系统越有利,但是过大的转动惯量J会使系统的动态特性变差,不利于维持系统的稳定性。
S3,根据虚拟阻抗、转动惯量和阻尼系数对角频率偏差和角频率变化率的影响规律确定VSG在暂态过程中的实时调整规则,具体包括:根据第一影响规则、第二影响规则和第三影响规则确定VSG在暂态过程中的的实时调整规则。
现在电力电子逆变器接口的分布式电源存在惯性和阻尼不足的缺点,不利于电网稳定,而同步发电机具有天然的并网友好性,但其转动惯量、阻尼系数和同步电抗都是固定的,因此本专利提出一种VSG的多参数自适应协同控制方法,使得并网逆变器具有与同步发电机相媲美的外特性,同时其参数还是可调的,增强了电网的抗干扰能力,结合步骤S1所述暂态响应过程分区和步骤S2所述影响规律,可得:
第一阶段t0-t1,又称其为加速区间,在此阶段中若适当减小虚拟阻抗的值,则可以使电磁功率Pe加速追踪额定功率Pref,同时减小角频率偏差ω的波动;若增大转动惯量,则可以增大平抑有功功率偏差的动能,角频率变化率就会减小,相当于阻止了系统频率的快速变化;同时增加阻尼系数也能防止出现频率过大的现象。
第二阶段t1-t2,又称其为减速区间,在此阶段中若适当增大虚拟阻抗的值,则可以使电磁功率Pe加速追踪额定功率Pref,同时会减小频率ω的波动;若减小转动惯量值,则角频率变化率会增大,则可以加快系统频率恢复至稳态值;同时增加阻尼系数也能在角频率偏差Δω较大时,进一步抑制角频率偏差。
第三阶段t2-t3,又称其为加速区间,此阶段中若适当减小虚拟阻抗的值,可以加速Pe的减小,使其快速追踪额定功率Pref,增大转动惯量和阻尼系数的值,则会增大平抑有功功率偏差的动能,角频率变化率会减小,抑制了频率的波动。
第四阶段t3-t4,又称其为减速区间,此阶段中若适当增大虚拟阻抗的值,可以使电磁功率Pe加速追踪额定功率Pref,减小转动惯量,加快频率的变化,使得系统频率快速恢复到额定频率,适当增大阻尼系数,防止出现频率过大的现象。
综上所述,得到虚拟阻抗、转动惯量和阻尼系数在暂态过程中的实时调整规则,如下表1所示:
表1
注:+表示增加,-表示减少。
S4,根据所述VSG的实时调整规则确定虚拟阻抗、转动惯量和阻尼系数的自适应函数。
逆变器是新能源和分布式发电接入智能电网的通用关键设备,引入VSG技术使得逆变器具有与同步发电机相媲美的性能,但是传统的VSG技术,因为其参数固定,往往难以满足电网的运行要求,因此将VSG与自适应相结合,可以提高VSG的灵活性和可靠性。
为了使VSG在分布式能源投切或者系统受到扰动,能够不断调节转动惯量、阻尼系数和虚拟阻抗来衰减振荡幅度,抑制频率变化,在频率调节中获得较好的动态响应,快速渡过不稳定的暂态过程,为此在保证电力系统安全稳定的前提下,结合新能源发电的要求以及响应特性,设计出转动惯量、阻尼系数和虚拟阻抗的自适应函数表达式。
S41,设计虚拟阻抗的自适应函数;
为了提升系统的稳定性,从虚拟阻抗方面出发,设计出的自适应虚拟阻抗表达式为:
式中:Lv0为初始虚拟电感,Δθv为虚拟阻抗角变化量,Rv0为初始虚拟电阻,Xv0为初始虚拟电抗,Zv0为初始虚拟阻抗。
S42,设计转动惯量和阻尼系数的自适应函数;
考虑到控制目的以及指数函数曲线的变化特点,本发明以指数函数为基础建立关于惯量的自适应函数式如下:
式中:J0为初始虚拟惯量,k1为虚拟惯量调节系数,始终为正;m为虚拟惯量调节指数,H为角频率变化的阈值。
为使系统的综合动态性能达到最佳,进一步得阻尼系数的函数表达式为:
式中:D0为初始阻尼系数,kd为阻尼系数的调节系数。
J、D、Lv由微电网功率偏差和角频率偏差共同决定,充分发挥了逆变器参数灵活可调的优势,解决了动态调节性能和转动惯量J之间的矛盾,同时克服了传统转动惯量J算法中角频率直接微分引入的系统噪声的问题。
S5,根据虚拟阻抗角变化量、频率变化量和最佳阻尼比确定所述虚拟阻抗、转动惯量和阻尼系数自适应函数的参数。
根据VSG的二阶电压方程确定初始虚拟电感Lvo,由此兼顾系统的响应速度和调节时间。在VSG系统中加入虚拟阻抗后,因为ωLv>>Rv,则改变Lv即可,进一步的通过对虚拟电感变化量ΔLv和虚拟阻抗角变化量Δθv的描述,在实现了功率解耦的情况下,提高VSG的动态响应能力。在自适应虚拟阻抗表达式中,比较重要的参数是虚拟阻抗角变化量Δθv,它的选取将直接影响虚拟电感的变化量ΔLv,根据VSG的二阶电压方程获得初始虚拟电感值,如式(12)所示,
式中:udref为VSG的虚拟内电势d轴分量,uqref为VSG的虚拟内电势q轴分量,uf为VSG的虚拟内电势幅值,id为VSG的输出电流d轴分量,iq为VSG的输出电流q轴分量。
式中:Δθv为虚拟阻抗角变化量。
式中:β为虚拟电感角补偿值得调整系数,Δθvmax为最大虚拟阻抗角变化量。
式中:kβ为常数,k2为有功功率变比。
由VSG的小信号模型,得到有功功率的闭环传递函数,由闭环传递函数得到闭环特征方程,画出不同转动惯量和阻尼系数的根轨迹图,进一步得到转动惯量和阻尼系数的变化范围,同时通过给出了四组不同调节指数m大小下,转动惯量J随调节系数k1和角频率偏差频率变化率dω/dt变化的三维曲线簇,进而得到调节指数m和调节系数k1的变化规则。
根据VSG的小信号模型,得到系统有功闭环传递函数如下:
可得到闭环特征方程为:
由式(15)得到其对应的自然振荡角频率ωn和阻尼比ζ为:
式中:E为VSG电势,U为VSG端电压,Kω为调差系数。
由式(16)画出不同转动惯量J和阻尼系数D下的根轨迹图如图3和图4所示。
由图3可以看出存在一对共轭复根,当J逐渐增大时,共轭复根慢慢向虚轴靠近,此时系统的超调量逐渐增大。当J增大到一定程度时,两个极点逐渐向零点靠拢,VSG的并网稳定性将受到严重的威胁。因此,要在一定范围内选择J的大小,为了充分利用逆变器的输出功率能力,初始转动惯量J0可根据式(18)选取:
由图4可以看出存在一对共轭复根,当D增大时,共轭复根的虚部不断减小,逐渐向实轴靠近,此时系统的阻尼随着增加,而系统的超调将逐渐减小;但当D增大到一定程度时,两个根在实轴上相反而行。此时系统由欠阻尼状态过渡到过阻尼状态,故D不能太大,一般情况下系统工作在欠阻尼状态,当ζ<0.4时,超调量较大;当ζ>0.8时,又会使得响应速度迟缓;所以,取ζ=0.707为“最佳阻尼比”,则初始阻尼系数D0可根据式(19)选取:
式中,Kp为有功频率下垂系数。
由图5可看出,当调节指数m由小变大时,三维曲面由凸变凹,即转动惯量在dω/dt较小时的数值范围逐渐减小,当dω/dt增加时,转动惯量的增长速度逐渐加快。故当调节指数m较大时,在负荷突变的情况下,随着dω/dt的增大,转动惯量J增长迅速,可减小系统频率的变化幅度。考虑微电网接口逆变器的额定功率和系统允许的最大角频率偏差频率偏差,配置调节系数k1,其表达式如下:
结果分析:
分别对以下几种情形下的VSG控制方法进行对比。如图6不同情形下有功功率变化对比和图7不同情形下频率变化对比所示,情形1:采用固定参数的VSG控制;情形2:虚拟惯量自适应的VSG控制;情形3:虚拟惯量和虚拟阻抗协同自适应的VSG控制;情形4:转动惯量、阻尼系数和虚拟阻抗协同自适应的VSG控制。从图中可以看出采用情形1、情形2、情形3与情形4的超调量σ%分别为20%、13%、4%、2%;调节时间ts分别约为0.48s、0.44s、0.37s、0.41s。说明本专利所提出的一种虚拟同步发电机的多参数自适应协同控制方法能使VSG的动态性能最好,验证了该方法的有效性。
本发明利用虚拟阻抗、转动惯量和阻尼系数的灵活可调节性有效提高了VSG系统的动态响应性能,提高了系统的频率稳定性。
以上显示和描述了本发明的基本原理、主要特征和优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (8)
1.一种虚拟同步发电机的多参数自适应协同控制方法,其特征在于,包括:
基于功率偏差和角频率变化率对虚拟同步发电机的暂态过程进行分区得到若干区域;
基于有功功率方程得到虚拟阻抗对所述角频率变化率的第一影响规则,基于转子运动方程获得阻尼系数对角频率偏差的第二影响规则和转动惯量对所述角频率变化率的第三影响规则;
基于所述第一影响规则、第二影响规则和第三影响规则得到若干所述区域对应的调整规则;
基于若干所述区域对应的调整规则确定虚拟阻抗、转动惯量和阻尼系数的自适应函数;
基于所述虚拟阻抗、转动惯量和阻尼系数的自适应函数得到各自适应函数对应的参数值;
基于有功功率方程得到虚拟阻抗对所述角频率变化率的第一影响规则的计算公式为:
式中,Ev为虚拟阻抗上的电压,E为逆变器的输出电压,Ev=IgZv,Zv=Rv+jωLv,Rv表示引入的虚拟电阻、Lv表示引入的虚拟电感,Zv为引入的虚拟阻抗,ω为VSG的同步转子角频率,Ig为逆变器的输出电流,Ug为逆变器的输出电压幅值,X为逆变器输出实际电抗,δ为功角;
构建的所述虚拟阻抗的自适应函数为:
式中,Lv0为初始虚拟电感,Δθv为虚拟阻抗角变化量,Rv0为初始虚拟电阻,Xv0为初始虚拟电抗,Zv0为初始虚拟阻抗,H表示角频率变化的阈值,ΔP表示额定功率与电磁功率的差值,Δω表示角频率偏差。
2.根据权利要求1所述的虚拟同步发电机的多参数自适应协同控制方法,其特征在于,基于转子运动方程获得阻尼系数对角频率偏差的第二影响规则的计算公式为:
式中,Δω表示角频率偏差,Tm表示机械转矩,Pm表示机械功率,Te表示电磁转矩,/>Pe为电磁功率,ω0为额定转子角频率,J为转动惯量,D为阻尼系数,ω为VSG的同步转子角频率,t表示时间。
3.根据权利要求2所述的虚拟同步发电机的多参数自适应协同控制方法,其特征在于,对所述第二影响规则的计算公式进行阻尼项去除且功率偏差一定时得到转动惯量对所述角频率变化率的第三影响规则,所述第三影响规则的计算公式为:
式中,表示角频率变化率。
4.根据权利要求1所述的虚拟同步发电机的多参数自适应协同控制方法,其特征在于,若干所述区域包括第一加速区域、第一减速区域、第二加速区域和第二减速区域;其中,各区域对应的调整规则包括:
当处于所述第一加速区域和第二加速区域时,采取减小虚拟阻抗,增大转动惯量和增加阻尼系数操作;
当处于所述第一减速区域和第二减速区域时,采取增大虚拟阻抗,减小转动惯量和增加阻尼系数操作。
5.根据权利要求1所述的虚拟同步发电机的多参数自适应协同控制方法,其特征在于,构建的所述转动惯量的自适应函数为:
式中,J0为初始虚拟惯量,k1为虚拟惯量调节系数,始终为正;m为虚拟惯量调节指数。
6.根据权利要求5所述的虚拟同步发电机的多参数自适应协同控制方法,其特征在于,构建的所述阻尼系数的自适应函数为:
式中,D0为初始阻尼系数,kd为阻尼系数的调节系数。
7.根据权利要求1所述的虚拟同步发电机的多参数自适应协同控制方法,其特征在于,基于所述虚拟阻抗的自适应函数得到初始虚拟电感值和虚拟电感变化量;
其中,获得所述初始虚拟电感值的过程包括:
基于VSG的二阶电压方程得到初始虚拟电感值;
获得所述虚拟电感变化量包括:
基于所述角频率偏差计算得到虚拟电感角补偿值的调整系数;
基于所述虚拟电感角补偿值的调整系数和所述角频率偏差得到虚拟阻抗角变化量;
通过改变所述虚拟阻抗角变化量得到虚拟电感的变化量。
8.根据权利要求7所述的虚拟同步发电机的多参数自适应协同控制方法,其特征在于,基于所述转动惯量的自适应函数确定初始转动惯量,基于所述阻尼系数的自适应函数确定初始阻尼系数;
其中,确定所述初始转动惯量和所述初始阻尼系数的过程包括:
基于小信号模型得到有功闭环传递函数;
基于所述有功闭环传递函数得到闭环特征方程;
基于所述闭环特征方程得到转动惯量的根轨迹和阻尼系数的根轨迹;
基于所述转动惯量的根轨迹得到所述初始转动惯量,基于所述阻尼系数的根轨迹得到初始阻尼系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311380654.7A CN117394432B (zh) | 2023-10-24 | 2023-10-24 | 一种虚拟同步发电机的多参数自适应协同控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311380654.7A CN117394432B (zh) | 2023-10-24 | 2023-10-24 | 一种虚拟同步发电机的多参数自适应协同控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117394432A CN117394432A (zh) | 2024-01-12 |
CN117394432B true CN117394432B (zh) | 2024-05-28 |
Family
ID=89469794
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311380654.7A Active CN117394432B (zh) | 2023-10-24 | 2023-10-24 | 一种虚拟同步发电机的多参数自适应协同控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117394432B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108270241A (zh) * | 2018-02-06 | 2018-07-10 | 国网四川省电力公司电力科学研究院 | 风机并网逆变器虚拟同步发电机的控制方法 |
CN111446741A (zh) * | 2020-04-03 | 2020-07-24 | 太原理工大学 | 一种用于直驱风力发电的虚拟同步发电机参数自适应方法 |
CN111525617A (zh) * | 2020-05-11 | 2020-08-11 | 广西大学 | 一种基于励磁调节的vsg功率解耦控制方法及系统 |
CN113098058A (zh) * | 2021-04-06 | 2021-07-09 | 广东电网有限责任公司电力科学研究院 | 一种转动惯量自适应优化控制方法、装置、设备及介质 |
CN113890117A (zh) * | 2021-10-20 | 2022-01-04 | 东南大学 | 一种基于自适应虚拟阻抗的功率解耦改进控制方法 |
CN114243761A (zh) * | 2021-11-25 | 2022-03-25 | 国网山东省电力公司烟台供电公司 | 一种微电网车并离网运行模式切换的控制方法及系统 |
CN115764989A (zh) * | 2022-11-04 | 2023-03-07 | 昆明理工大学 | 一种基于构网型换流器的虚拟同步发电机系统 |
CN115912405A (zh) * | 2022-12-21 | 2023-04-04 | 华能国际电力股份有限公司大连电厂 | 一种用于复杂振荡环境下虚拟同步发电机自适应控制策略 |
CN116388224A (zh) * | 2023-04-11 | 2023-07-04 | 华北电力大学 | 一种基于自适应虚拟同步阻抗的动态解耦控制方法 |
CN116454978A (zh) * | 2023-05-09 | 2023-07-18 | 厦门理工学院 | 一种虚拟同步发电机控制方法 |
WO2023178873A1 (zh) * | 2022-03-22 | 2023-09-28 | 西安热工研究院有限公司 | 火电熔融盐储能换流装置的并网点电压调节方法及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102017106213A1 (de) * | 2017-03-22 | 2018-09-27 | Wobben Properties Gmbh | Verfahren zum Einspeisen elektrischer Leistung in ein elektrisches Versorgungsnetz |
-
2023
- 2023-10-24 CN CN202311380654.7A patent/CN117394432B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108270241A (zh) * | 2018-02-06 | 2018-07-10 | 国网四川省电力公司电力科学研究院 | 风机并网逆变器虚拟同步发电机的控制方法 |
CN111446741A (zh) * | 2020-04-03 | 2020-07-24 | 太原理工大学 | 一种用于直驱风力发电的虚拟同步发电机参数自适应方法 |
CN111525617A (zh) * | 2020-05-11 | 2020-08-11 | 广西大学 | 一种基于励磁调节的vsg功率解耦控制方法及系统 |
CN113098058A (zh) * | 2021-04-06 | 2021-07-09 | 广东电网有限责任公司电力科学研究院 | 一种转动惯量自适应优化控制方法、装置、设备及介质 |
CN113890117A (zh) * | 2021-10-20 | 2022-01-04 | 东南大学 | 一种基于自适应虚拟阻抗的功率解耦改进控制方法 |
CN114243761A (zh) * | 2021-11-25 | 2022-03-25 | 国网山东省电力公司烟台供电公司 | 一种微电网车并离网运行模式切换的控制方法及系统 |
WO2023178873A1 (zh) * | 2022-03-22 | 2023-09-28 | 西安热工研究院有限公司 | 火电熔融盐储能换流装置的并网点电压调节方法及系统 |
CN115764989A (zh) * | 2022-11-04 | 2023-03-07 | 昆明理工大学 | 一种基于构网型换流器的虚拟同步发电机系统 |
CN115912405A (zh) * | 2022-12-21 | 2023-04-04 | 华能国际电力股份有限公司大连电厂 | 一种用于复杂振荡环境下虚拟同步发电机自适应控制策略 |
CN116388224A (zh) * | 2023-04-11 | 2023-07-04 | 华北电力大学 | 一种基于自适应虚拟同步阻抗的动态解耦控制方法 |
CN116454978A (zh) * | 2023-05-09 | 2023-07-18 | 厦门理工学院 | 一种虚拟同步发电机控制方法 |
Non-Patent Citations (5)
Title |
---|
Adaptive Backstepping Design for Stabilizing Synchronverter Control Topology in AC Microgrids;Mobin Naderi 等;《2019 Smart Grid Conference (SGC)》;20191231;1-6 * |
Adaptive Virtual Inertia-Damping System Based on Model Predictive Control for Low-Inertia Microgrids;ASMAA FAWZY 等;《IEEE Access》;20210831;109718-109731 * |
光伏微电网故障电流控制研究;陈天胜 等;《电力电子技术》;20200229;第54卷(第2期);85-88 * |
孤岛微电网的电压不平衡二次控制策略;吴丽珍 等;《兰州理工大学学报》;20161031;第42卷(第5期);84-89 * |
虚拟同步发电机接口变换器并联运行 虚拟阻抗自适应控制;温春雪 等;《电工技术学报》;20201231;第35卷;494-502 * |
Also Published As
Publication number | Publication date |
---|---|
CN117394432A (zh) | 2024-01-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Power-flow control and stability enhancement of four parallel-operated offshore wind farms using a line-commutated HVDC link | |
CN110957714B (zh) | 平抑充电站直流微网电压波动的虚拟直流电机控制方法 | |
Salama et al. | Frequency and voltage control of microgrid with high WECS penetration during wind gusts using superconducting magnetic energy storage | |
CN110011364A (zh) | 一种降低有功功率载荷波动对系统稳定影响的控制方法 | |
CN110518600A (zh) | 一种基于vsg的改进型多机并联的pmsg并网主动支撑控制结构 | |
Li et al. | An improved VSG control strategy based on the amplitude-frequency characteristics of virtual power | |
CN111478365B (zh) | 一种直驱风电机组虚拟同步机控制参数的优化方法及系统 | |
CN112993991A (zh) | 一种风机双通道阻尼控制低频振荡广域阻尼控制方法 | |
CN115714435A (zh) | 基于虚拟同步发电机的光伏混合储能系统功率分配及虚拟惯性控制方法 | |
CN115764989A (zh) | 一种基于构网型换流器的虚拟同步发电机系统 | |
CN111327074A (zh) | 一种基于虚拟同步机技术的直流侧混合储能协调控制方法 | |
CN108879893A (zh) | 一种基于准pr控制的电动汽车虚拟同步电机快充方法 | |
CN116845924B (zh) | 基于相角自生成策略的构网型风电场电压源控制方法 | |
CN112787325A (zh) | 基于李雅普诺夫直接法的永磁同步风力发电机暂态同步稳定性的定量评估方法 | |
CN117394432B (zh) | 一种虚拟同步发电机的多参数自适应协同控制方法 | |
CN116937546A (zh) | 一种考虑风储并网的电网低频振荡抑制方法及系统 | |
Li et al. | VSG virtual inertial control strategy based on lead-lag link and fuzzy logic control | |
CN116683491A (zh) | 一种新能源微电网的惯量控制方法 | |
CN108736517B (zh) | 一种基于vsg的逆变型分布式电源自适应阻尼控制方法 | |
CN113937789B (zh) | 基于分数阶滤波的电压源型双馈风机前馈阻尼控制方法 | |
CN115425636A (zh) | 含飞轮储能的直流微电网虚拟惯量自适应控制方法 | |
CN115276084A (zh) | 一种应用于水电站柴油发电机黑启动的模糊pid与混合储能协同控制方法 | |
CN113839577A (zh) | 基于模糊控制的三电平npc逆变器结构虚拟同步发电机 | |
Shao et al. | Adaptive damping coefficient control of virtual synchronous generator of microgrid inverter | |
Li et al. | The Small‐Signal Stability of Offshore Wind Power Transmission Inspired by Particle Swarm Optimization |
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 |