CN106099996B - 含状态预测器的虚拟同步发电机一致性的控制方法 - Google Patents
含状态预测器的虚拟同步发电机一致性的控制方法 Download PDFInfo
- Publication number
- CN106099996B CN106099996B CN201610459615.XA CN201610459615A CN106099996B CN 106099996 B CN106099996 B CN 106099996B CN 201610459615 A CN201610459615 A CN 201610459615A CN 106099996 B CN106099996 B CN 106099996B
- Authority
- CN
- China
- Prior art keywords
- output
- difference delta
- angular frequency
- vsg
- consistency
- 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
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
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/48—Controlling the sharing of the in-phase component
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Eletrric Generators (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开一种含状态预测器的虚拟同步发电机一致性的控制方法,所述控制方法包括:根据第i台虚拟同步发电机VSG的输出功率Pi和额定功率Pi *确定第一输出差值Δyi;根据所述第一输出差值Δyi和第i台VSG的下垂特性曲线平移量Yi确定第二输出差值ΔYi,且所述第二输出差值ΔYi与一致性控制参数Qi调整所述下垂特性曲线平移量Yi;根据所述第二输出差值、第i台VSG的阻尼系数Dp,i、额定角频率ωb、实测的角频率ωi、虚拟转动惯量Ji确定输出角频率Δωi;根据实测的角频率ωi确定角度θi,并根据角度θi影响输出功率Pi及其一致收敛的速度。本发明含状态预测器的虚拟同步发电机一致性的控制方法通过第二输出差值ΔYi与一致性控制参数Qi调整下垂特性曲线平移量Yi,从而可提高收敛速度,以迅速到达收敛值。
Description
技术领域
本发明涉及虚拟同步发电机的技术领域,特别是涉及一种含状态预测器的虚拟同步发电机一致性的控制方法。
背景技术
随着分布式电源的规模不断增大,以逆变器为主要接口的电源越来越多接入微电网。传统逆变器几乎没有转动惯量,难以为电网提供惯性和阻尼,因此无法满足支撑频率和电压的需求。当大量分布式电源接入微电网后,会给微电网安全稳定运行带来巨大威胁。基于虚拟同步发电机(Virtual Synchronous Generator,VSG)的并网逆变器控制策略能够起到为微电网提供惯性和阻尼的作用。基于一致性的虚拟同步发电机技术能够通过相邻VSG间的通信,实现系统频率的一致恢复以及有功出力的合理分配,然而由于没有考虑系统一致收敛速度的问题,导致系统收敛速度较慢。
发明内容
本发明的目的是提供一种含状态预测器的虚拟同步发电机一致性的控制方法,可有效提高系统输出有功功率的一致性的收敛速度。
为实现上述目的,本发明提供了如下方案:
一种含状态预测器的虚拟同步发电机一致性的控制方法,所述控制方法包括:
根据第i台虚拟同步发电机VSG的输出功率Pi和额定功率Pi *确定第一输出差值Δyi;
根据所述第一输出差值Δyi和第i台VSG的下垂特性曲线平移量Yi确定第二输出差值ΔYi,且所述第二输出差值ΔYi与一致性控制参数Qi调整所述下垂特性曲线平移量Yi;
根据所述第二输出差值、第i台VSG的阻尼系数Dp,i、额定角频率ωb、实测的角频率ωi、虚拟转动惯量Ji调整输出角频率Δωi,且所述输出角频率Δωi与所述额定角频率ωb确定实测的角频率ωi;
根据实测的角频率ωi确定角度θi,并根据角度θi影响输出功率Pi及所述输出功率Pi的一致收敛的速度。
可选的,所述确定第一输出差值Δyi的方法包括:
根据输出功率Pi确定输出目标函数f(Pi);
根据额定功率Pi *确定额定目标函数f(Pi *);
根据以下公式以及输出目标函数f(Pi)和额定目标函数f(Pi *)确定第一输出差值Δyi:
其中a1,b1分别为选择系数。
可选的,所述输出目标函数f(Pi)的确定方法包括:
按照额定功率进行输出功率的分配,即满足下式时:
其中,m表示VSG的台数,i=1,...,m;
定目标函数f(Pi)选取如下:
f(Pi)=Pi *-1Pi;选择系数a1=1,b1=0。
可选的,所述输出目标函数f(Pi)的确定方法包括:
按照等微增率进行输出功率的分配,即满足下式时:
λ1(P1)=...λi(Pi)=...=λm(Pm);其中,m表示VSG的台数,i=1,...,m;λi(Pi)为第i台VSG的微增率函数;
根据以下公式确定目标函数f(Pi):
f(Pi)=λi(Pi);选择系数a1=0,b1=1。
可选的,所述第二输出差值ΔYi的确定方法包括:
根据以下公式,使用状态预测器Yi t确定一致性控制参数Qi:
其中,L表示拉普拉斯矩阵,lij表示所述L中的元素;ky,i表示过程参数,γ表示状态预测器的影响因子;
根据以下公式确定所述第二输出差值ΔYi:
其中,a2,b2分别选择系数,且a2取值与a1相同,b2取值与b1相同;ki表示恢复系数。
可选的,所述根据所述第二输出差值ΔYi、第i台VSG的阻尼系数Dp,i、额定角频率ωb、实测的角频率ωi、虚拟转动惯量Ji调整输出角频率Δωi的方法包括:
其中,Δωi=ωi-ωb。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明含状态预测器的虚拟同步发电机一致性的控制方法通过第二输出差值ΔYi与一致性控制参数Qi调整下垂特性曲线平移量Yi,从而可有效提高系统输出有功功率一致性的收敛速度,以迅速到达收敛值。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明含状态预测器的虚拟同步发电机一致性的控制方法的流程图;
图2为并联一致性的控制框图;
图3为三个VSG的通信拓扑图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种含状态预测器的虚拟同步发电机一致性的控制方法,通过第二输出差值ΔYi与一致性控制参数Qi调整下垂特性曲线平移量Yi,从而可有效提高系统输出有功功率一致性的收敛速度,以迅速到达收敛值。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
如图1和图2所示,本发明含状态预测器的虚拟同步发电机一致性的控制方法包括:
步骤100:根据第i台虚拟同步发电机VSG的输出功率Pi和额定功率Pi *确定第一输出差值Δyi;
步骤200:根据所述第一输出差值Δyi和第i台VSG的下垂特性曲线平移量Yi确定第二输出差值ΔYi,且所述第二输出差值ΔYi与一致性控制参数Qi调整所述下垂特性曲线平移量Yi;
步骤300:根据所述第二输出差值ΔYi、第i台VSG的阻尼系数Dp,i、额定角频率ωb、实测的角频率ωi、虚拟转动惯量Ji调整输出角频率Δωi,且所述输出角频率Δωi与所述额定角频率ωb确定实测的角频率ωi;
步骤400:根据实测的角频率ωi确定角度θi,并根据角度θi影响输出功率Pi及所述输出功率Pi的一致收敛的速度。
其中,在步骤100中,所述确定第一输出差值Δyi的方法包括:
步骤101:根据输出功率Pi确定输出目标函数f(Pi);
步骤102:根据额定功率Pi *确定额定目标函数f(Pi *);
步骤103:根据以下公式以及输出目标函数f(Pi)和额定目标函数f(Pi *)确定第一输出差值Δyi:
其中a1,b1分别为选择系数。
输出目标函数f(Pi)和额定目标函数f(Pi *)一般为一次函数,其选取方法较多,可根据具体的需求确定。例如:
(1)按照额定功率进行输出功率的分配,即满足以下公式,:
其中,m表示VSG的台数,i=1,...,m;
目标函数f(Pi)选取如下:
f(Pi)=Pi *-1Pi------公式(3);
选择系数a1=1,b1=0。
(2)按等微增率原则进行功率分配,即满足以下公式:
λ1(P1)=...λi(Pi)=...=λm(Pm)------公式(4);
其中,m表示VSG的台数,i=1,...,m;λi(Pi)为第i台VSG的微增率函数;
目标函数f(Pi):f(Pi)=λi(Pi)------公式(5);
选择系数a1=0,b1=1。
在步骤200中,所述第二输出差值ΔYi的确定方法包括:
步骤201:根据以下公式,使用状态预测器Yi t确定一致性控制参数Qi:
其中,L表示拉普拉斯矩阵,lij表示所述L中的元素;ky,i表示过程参数,γ表示状态预测器的影响因子。
步骤202:根据以下公式确定所述第二输出差值ΔYi:
其中,a2,b2分别选择系数,且a2取值与a1相同,b2取值与b1相同;ki表示恢复系数。
进一步地,在步骤300中,所述根据所述第二输出差值ΔYi、第i台VSG的阻尼系数Dp,i、额定角频率ωb、实测的角频率ωi、虚拟转动惯量Ji调整输出角频率Δωi的方法包括:
其中,Δωi=ωi-ωb。
其中,当多个VSG并联连接时,形成通信拓扑结构,则邻接矩阵A=[aij](两者相连为1,否则为0,对角线上元素为0),di=∑aij,度矩阵D=diag{di}(度矩阵中的非对角元素均为0,对角线上元素大小等于VSG连接的通信线路条数);而拉普拉斯矩阵L=D-A。
如图3所示,3个VSG并联连接时,形成通信拓扑结构,则邻接矩阵度矩阵则拉普拉斯矩阵
当系统中各VSG构成的通讯拓扑满足连通性时,即存在生成树是系统能够一致收敛的必要条件。对于无向图,拉普拉斯矩阵L为半正定矩阵,如果无向图是连通的,则它的最小非零特征值称为L的代数连通度。即:
其中,x表示系统状态变量。
若λ2(L)>0,则整个通讯拓扑是连通的,而且λ2(L)越大,则说明在系统中各VSG之间的直接通讯链路越多。若λ2(L)=0,则说明整个通讯拓扑不是连通的。且系统一致性演化速度和拉普拉斯矩阵的最小非零特征值λ2相关,λ2越大,VSG系统演化到平衡状态的速度也越快。
拉普拉斯矩阵L在无向图通信模型下,且满足生成树条件时,存在一个非奇异矩阵T,拉普拉斯矩阵可表示为
可以发现,系统的最小非零特征值λ2是由系统本身的通信拓扑所决定的,且无法改变。
而在本发明中引入状态预测器后,可将公式(5)和(6)转换为
公式(9)中等号右侧包括两部分,一部分为ΔYi=Δyi-ky,iYi,剩下的为第二部分通过对比,发现在时间尺度上,第一部分的收敛速度远远快于第二部分,所以考虑一致性收敛速度仅需考虑第二部分,即将公式(9)转化为:
公式(10)可写成矩阵形式为:
将Yi t=-LYi代入公式(11)可得
合并同类项可得
根据公式(8)可知,矩阵L+γL2结果为:
显然对于任意γ>0,则具有状态预测器的系统比没有状态预测器的系统有更大的最小非零特征值,因此能更快达到一致。且最小非零特征值越大,系统一致性收敛速度越快,从而能够迅速到达收敛值。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (3)
1.一种含状态预测器的虚拟同步发电机一致性的控制方法,其特征在于,所述控制方法包括:
根据第i台虚拟同步发电机VSG的输出功率Pi和额定功率Pi *确定第一输出差值Δyi;
根据所述第一输出差值Δyi和第i台VSG的下垂特性曲线平移量Yi确定第二输出差值ΔYi,且所述第二输出差值ΔYi与一致性控制参数Qi调整所述下垂特性曲线平移量Yi;
根据所述第二输出差值ΔYi、第i台VSG的阻尼系数Dp,i、额定角频率ωb、实测的角频率ωi、虚拟转动惯量Ji调整输出角频率Δωi,且所述输出角频率Δωi与所述额定角频率ωb确定实测的角频率ωi;
根据实测的角频率ωi确定角度θi,并根据角度θi影响输出功率Pi及所述输出功率Pi的一致收敛的速度;
其中,所述确定第一输出差值Δyi的方法包括:
根据输出功率Pi确定输出目标函数f(Pi);
根据额定功率Pi *确定额定目标函数f(Pi *);
根据以下公式以及输出目标函数f(Pi)和额定目标函数f(Pi *)确定第一输出差值Δyi:
其中a1,b1分别为选择系数;
所述第二输出差值ΔYi的确定方法包括:
根据以下公式,使用状态预测器确定一致性控制参数Qi:
其中,L表示拉普拉斯矩阵,lij表示所述L中的元素;ky,i表示过程参数,γ表示状态预测器的影响因子;
根据以下公式确定所述第二输出差值ΔYi:
其中,a2,b2分别选择系数,且a2取值与a1相同,b2取值与b1相同;ki表示恢复系数;
所述根据所述第二输出差值ΔYi、第i台VSG的阻尼系数Dp,i、额定角频率ωb、实测的角频率ωi、虚拟转动惯量Ji调整输出角频率Δωi的方法包括:
其中,Δωi=ωi-ωb。
2.根据权利要求1所述的含状态预测器的虚拟同步发电机一致性的控制方法,其特征在于,所述输出目标函数f(Pi)的确定方法包括:
按照额定功率进行输出功率的分配,即满足下式时:
其中,m表示VSG的台数,i=1,...,m;
目标函数f(Pi)选取如下:
选择系数a1=1,b1=0。
3.根据权利要求1所述的含状态预测器的虚拟同步发电机一致性的控制方法,其特征在于,所述输出目标函数f(Pi)的确定方法包括:
按照等微增率进行输出功率的分配,即满足下式时:
λ1(P1)=...λi(Pi)=...=λm(Pm);其中,m表示VSG的台数,i=1,...,m;λi(Pi)为第i台VSG的微增率函数;
根据以下公式确定目标函数f(Pi):
f(Pi)=λi(Pi);选择系数a1=0,b1=1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610459615.XA CN106099996B (zh) | 2016-06-22 | 2016-06-22 | 含状态预测器的虚拟同步发电机一致性的控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610459615.XA CN106099996B (zh) | 2016-06-22 | 2016-06-22 | 含状态预测器的虚拟同步发电机一致性的控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106099996A CN106099996A (zh) | 2016-11-09 |
CN106099996B true CN106099996B (zh) | 2018-10-12 |
Family
ID=57252097
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610459615.XA Active CN106099996B (zh) | 2016-06-22 | 2016-06-22 | 含状态预测器的虚拟同步发电机一致性的控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106099996B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108931721B (zh) * | 2017-05-27 | 2021-04-06 | 中国电力科学研究院 | 一种发电机阻尼性质判别方法及装置 |
CN109149658B (zh) * | 2018-08-23 | 2021-11-19 | 河海大学 | 基于一致性理论的独立微电网分布式动态经济调度方法 |
CN110071514B (zh) * | 2019-05-14 | 2022-11-29 | 电子科技大学 | 一种用于功率分配和电压频率恢复的一致性下垂控制方法 |
CN113162122B (zh) * | 2021-02-26 | 2022-08-02 | 华中科技大学 | 一种基于一致性算法的虚拟同步机控制方法和系统 |
CN113381442A (zh) * | 2021-05-21 | 2021-09-10 | 青海大学 | 机组分布式控制方法和系统、电子设备及可读存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104158181A (zh) * | 2014-08-14 | 2014-11-19 | 清华大学 | 一种基于虚拟同步发电机的可逆变流岸电控制方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013046503A (ja) * | 2011-08-25 | 2013-03-04 | Waseda Univ | 蓄電システム及びその制御方法 |
-
2016
- 2016-06-22 CN CN201610459615.XA patent/CN106099996B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104158181A (zh) * | 2014-08-14 | 2014-11-19 | 清华大学 | 一种基于虚拟同步发电机的可逆变流岸电控制方法 |
Non-Patent Citations (2)
Title |
---|
《具有状态预测期的多智能体系统一致性研究》;席裕庚等;《控制与决策》;20100531;第25卷(第5期);第769-772页 * |
基于虚拟同步发电机的柔性虚拟调速器模型;赵杨阳等;《电力系统自动化》;20160525;第40卷(第10期);第8-15页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106099996A (zh) | 2016-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106099996B (zh) | 含状态预测器的虚拟同步发电机一致性的控制方法 | |
Pulgar-Painemal et al. | On inertia distribution, inter-area oscillations and location of electronically-interfaced resources | |
Zhao et al. | Distributed plug-and-play optimal generator and load control for power system frequency regulation | |
CN106208102B (zh) | 一种基于辅助问题原理的主动配电网分布式无功优化方法 | |
Zhou et al. | Distributed secondary control for islanded microgrids with mobile emergency resources | |
Hasanien et al. | Water cycle algorithm‐based optimal control strategy for efficient operation of an autonomous microgrid | |
CN107579543A (zh) | 一种基于分层控制策略的孤岛微电网分布式协调控制方法 | |
Lai et al. | Distributed power control for DERs based on networked multiagent systems with communication delays | |
CN106374513B (zh) | 一种基于主从博弈的多微网联络线功率优化方法 | |
de Haan et al. | Stabilising system frequency using HVDC between the Continental European, Nordic, and Great Britain systems | |
CN106340907A (zh) | 一种电力系统安全稳定控制策略确定方法及装置 | |
WO2015092705A1 (en) | Methods and systems for power injection or extraction in a power network | |
CN109449981A (zh) | 一种考虑通信延时的分布式频率控制稳定性判断方法 | |
CN106300358B (zh) | 一种基于π型拓扑结构的交直流并联分散控制方法 | |
CN109193765A (zh) | 一种分布式光伏集群调控方法及装置 | |
CN108959671B (zh) | 半桥和全桥型模块化多电平换流器的实时仿真建模方法 | |
CN108493940A (zh) | 一种含柔性直流的交直流混合系统潮流交替迭代计算方法 | |
CN108448563A (zh) | 一种直流微电网分布式协同控制系统及直流微电网 | |
Ferreira et al. | Distributed control strategy for low-voltage three-phase four-wire microgrids: Consensus power-based control | |
Renner et al. | Potential of using DC voltage restoration reserve for HVDC grids | |
Cai et al. | A network model for the real-time communications of a smart grid prototype | |
CN108988320A (zh) | 电力系统动态元件响应特性对暂态电压稳定性影响分析方法 | |
Pant et al. | Periodic output feedback technique based design of STATCOM damping controller | |
CN109586275A (zh) | 交直流混联系统电磁暂态仿真的快速启动方法与装置 | |
CN106451421A (zh) | 一种含电磁环网交直流电网的建模与解环方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |