CN106099996A - 含状态预测器的虚拟同步发电机一致性的控制方法 - Google Patents
含状态预测器的虚拟同步发电机一致性的控制方法 Download PDFInfo
- Publication number
- CN106099996A CN106099996A CN201610459615.XA CN201610459615A CN106099996A CN 106099996 A CN106099996 A CN 106099996A CN 201610459615 A CN201610459615 A CN 201610459615A CN 106099996 A CN106099996 A CN 106099996A
- Authority
- CN
- China
- Prior art keywords
- output
- difference delta
- vsg
- output difference
- virtual synchronous
- 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.)
- Granted
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 (6)
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的一致收敛的速度。
2.根据权利要求1所述的含状态预测器的虚拟同步发电机一致性的控制方法,其特征在于,所述确定第一输出差值Δyi的方法包括:
根据输出功率Pi确定输出目标函数f(Pi);
根据额定功率Pi *确定额定目标函数f(Pi *);
根据以下公式以及输出目标函数f(Pi)和额定目标函数f(Pi *)确定第一输出差值Δyi:
其中a1,b1分别为选择系数。
3.根据权利要求2所述的含状态预测器的虚拟同步发电机一致性的控制方法,其特征在于,所述输出目标函数f(Pi)的确定方法包括:
按照额定功率进行输出功率的分配,即满足下式时:
其中,m表示VSG的台数,i=1,...,m;
目标函数f(Pi)选取如下:
f(Pi)=Pi -1Pi;选择系数a1=1,b1=0。
4.根据权利要求2所述的含状态预测器的虚拟同步发电机一致性的控制方法,其特征在于,所述输出目标函数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。
5.根据权利要求2所述的含状态预测器的虚拟同步发电机一致性的控制方法,其特征在于,所述第二输出差值ΔYi的确定方法包括:
根据以下公式,使用状态预测器确定一致性控制参数Qi:
其中,L表示拉普拉斯矩阵,lij表示所述L中的元素;ky,i表示过程参数,γ表示状态预测器的影响因子;
根据以下公式确定所述第二输出差值ΔYi:
其中,a2,b2分别选择系数,且a2取值与a1相同,b2取值与b1相同;ki表示恢复系数。
6.根据权利要求1-5中任一项所述的含状态预测器的虚拟同步发电机一致性的控制方法,其特征在于,所述根据所述第二输出差值ΔYi、第i台VSG的阻尼系数Dp,i、额定角频率ωb、实测的角频率ωi、虚拟转动惯量Ji调整输出角频率Δωi的方法包括:
其中,Δωi=ωi-ωb。
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 true CN106099996A (zh) | 2016-11-09 |
CN106099996B 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) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108931721A (zh) * | 2017-05-27 | 2018-12-04 | 中国电力科学研究院 | 一种发电机阻尼性质判别方法及装置 |
CN109149658A (zh) * | 2018-08-23 | 2019-01-04 | 河海大学 | 基于一致性理论的独立微电网分布式动态经济调度方法 |
CN110071514A (zh) * | 2019-05-14 | 2019-07-30 | 电子科技大学 | 一种用于功率分配和电压频率恢复的一致性下垂控制方法 |
CN113162122A (zh) * | 2021-02-26 | 2021-07-23 | 华中科技大学 | 一种基于一致性算法的虚拟同步机控制方法和系统 |
CN113381442A (zh) * | 2021-05-21 | 2021-09-10 | 青海大学 | 机组分布式控制方法和系统、电子设备及可读存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013046503A (ja) * | 2011-08-25 | 2013-03-04 | Waseda Univ | 蓄電システム及びその制御方法 |
CN104158181A (zh) * | 2014-08-14 | 2014-11-19 | 清华大学 | 一种基于虚拟同步发电机的可逆变流岸电控制方法 |
-
2016
- 2016-06-22 CN CN201610459615.XA patent/CN106099996B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013046503A (ja) * | 2011-08-25 | 2013-03-04 | Waseda Univ | 蓄電システム及びその制御方法 |
CN104158181A (zh) * | 2014-08-14 | 2014-11-19 | 清华大学 | 一种基于虚拟同步发电机的可逆变流岸电控制方法 |
Non-Patent Citations (2)
Title |
---|
席裕庚等: "《具有状态预测期的多智能体系统一致性研究》", 《控制与决策》 * |
赵杨阳等: "基于虚拟同步发电机的柔性虚拟调速器模型", 《电力系统自动化》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108931721A (zh) * | 2017-05-27 | 2018-12-04 | 中国电力科学研究院 | 一种发电机阻尼性质判别方法及装置 |
CN108931721B (zh) * | 2017-05-27 | 2021-04-06 | 中国电力科学研究院 | 一种发电机阻尼性质判别方法及装置 |
CN109149658A (zh) * | 2018-08-23 | 2019-01-04 | 河海大学 | 基于一致性理论的独立微电网分布式动态经济调度方法 |
CN109149658B (zh) * | 2018-08-23 | 2021-11-19 | 河海大学 | 基于一致性理论的独立微电网分布式动态经济调度方法 |
CN110071514A (zh) * | 2019-05-14 | 2019-07-30 | 电子科技大学 | 一种用于功率分配和电压频率恢复的一致性下垂控制方法 |
CN110071514B (zh) * | 2019-05-14 | 2022-11-29 | 电子科技大学 | 一种用于功率分配和电压频率恢复的一致性下垂控制方法 |
CN113162122A (zh) * | 2021-02-26 | 2021-07-23 | 华中科技大学 | 一种基于一致性算法的虚拟同步机控制方法和系统 |
CN113162122B (zh) * | 2021-02-26 | 2022-08-02 | 华中科技大学 | 一种基于一致性算法的虚拟同步机控制方法和系统 |
CN113381442A (zh) * | 2021-05-21 | 2021-09-10 | 青海大学 | 机组分布式控制方法和系统、电子设备及可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN106099996B (zh) | 2018-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106099996A (zh) | 含状态预测器的虚拟同步发电机一致性的控制方法 | |
Robbins et al. | Optimal tap setting of voltage regulation transformers in unbalanced distribution systems | |
CN107147151B (zh) | 一种同步发电机逆变器的并联功率分配控制方法 | |
CN105140907B (zh) | 直流微网多智能体自适应下垂一致性协调控制方法及装置 | |
CN106777827A (zh) | 一种机电‑电磁混合仿真方法及系统 | |
CN103995734A (zh) | 基于rtds的电力系统混合实时仿真系统及仿真方法 | |
CN107579543A (zh) | 一种基于分层控制策略的孤岛微电网分布式协调控制方法 | |
CN107342586A (zh) | 一种用于确定配电网静态等值的方法及装置 | |
CN108959671B (zh) | 半桥和全桥型模块化多电平换流器的实时仿真建模方法 | |
Tu et al. | Transient stability constrained optimal power flow using independent dynamic simulation | |
CN103678798B (zh) | 一种用于含分布式电源配电网的电磁暂态实时仿真方法 | |
CN109149665A (zh) | 用于柔性直流交流电网联合仿真的多速率仿真方法及系统 | |
Wang et al. | Phase balancing in power distribution network with data center | |
CN103715687B (zh) | 一种主动配电网支路功率拥塞实时控制方法 | |
CN108039720A (zh) | 一种交直流混联系统最大输电能力的确定方法和装置 | |
CN107194039A (zh) | 一种基于改进高斯伪谱法的空间柔性系统展开控制方法 | |
CN108347057B (zh) | 一种lcc-mmc混合直流电网交替迭代潮流计算方法 | |
CN108614416B (zh) | 光伏虚拟同步发电机控制参数的整定方法及系统 | |
CN104578055A (zh) | 一种基于网络节点编号优化的电力系统暂态稳定仿真方法 | |
CN106877309A (zh) | 主从控制模式下直流微电网的稳定性判定方法 | |
CN108427821A (zh) | 一种交直流混联电网运行特性的数模混合仿真方法及系统 | |
CN107196337A (zh) | 一种孤岛静态工作点校验方法及系统 | |
CN107370149A (zh) | 确定电力系统负荷削减量的方法和系统 | |
CN108494017B (zh) | 一种基于逆变器的自治型微电网系统分布式协调控制方法 | |
CN104184165B (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 |