一种分布式新能源互联微电网集群振荡控制方法

1.本发明涉及多源互联微电网系统稳定控制技术领域，特别是涉及一种分布式新能源互联微电网集群振荡控制方法。


背景技术：

2.随着双碳目标的提出与能源发展要求的不断提高，优先发展可再生能源是推动我国能源供给革命，实现高质量可持续发展的必由之路。分布式发电技术是高效利用可再生能源的供能手段，当前主要包含风力发电、光伏发电、储能和微型燃气轮机等形式。这些分布式电源大量采用具有弱阻尼惯量特征的电力电子变换器接入系统，在一定程度上提高了响应速度，但却降低了系统阻尼，尤其是当其自主运行时，由于缺乏上级大电网的支撑，需自我维持频率电压稳定，再加上其小容量、低惯性的特征，此时一个微小的故障也可能导致系统性能受到巨大影响。同时，以风电、光伏为代表的可再生能源，由于其随机性和间接性的特征将会引发功率波动造成系统的频率偏移，严重影响系统的供电质量，对系统的安全稳定运行带来了严峻挑战。
3.互联微电网系统中包含不同类型的微电源和负荷，极易产生微电源间的交互作用、动态负荷间的交互作用、微电源与负荷间的交互作用等，在遭遇扰动时容易产生振荡甚至失稳，多种控制主体可能同时参与系统振荡。因此，针对新能源发电机组高渗透率的多源互联微电网系统，研究其在不同扰动场景下的动态响应特性以及设计一种抑制振荡的方法显得尤为重要。
4.目前在实际电力系统中普遍以同步发电机的励磁调节器加装电力系统稳定器(power system stabilization，pss)为主要阻尼控制器，但pss对于低频振荡特别是区间振荡并不总是有效，加之随着以风电和光伏等新能源发电渗透率越来越高，传统阻尼控制手段将难以独立应对电力系统运行特性出现的新变化。同时，在以电力电子变换器为主导的多源互联微电网系统的振荡抑制策略中，主要分为附加装置或控制回路附加额外阻尼以及优化控制参数设计，这些控制方法均存在面临在工况复杂多变时适应性不强的问题。而且目前已有的新能源附加阻尼控制主要针对一种风电或光伏发电机组展开，针对多源互联微电网系统的新能源机组特性的协同阻尼控制研究并不多。为充分发掘新能源发电在参与抑制多源互联微电网系统振荡方面的潜力，需要在考虑新能源机组特性基础上设计更具工况适应性的控制策略，以适应未来新型电力系统以分布式风电、光伏发电单元为主导的发展趋势。
5.现有技术公开了一种孤岛微电网系统及其交互振荡抑制方法、系统，源侧虚拟同步机和负载侧虚拟同步机结合起来提供虚拟惯性和阻尼。该发明提出了一种d轴电感电流前馈控制和d轴电压反馈控制的方法来重塑vsg的阻抗，以减小vsg的阻抗幅值。因此，vsg和lvsm之间的低频交互作用可以得到抑制。该现有技术中针对的是单个微电网系统包含多个相同的电压源型换流器，且并未考虑新能源机组的特性，将直流侧用恒定电压源替代，此外，该方案采用附加控制方法中的虚拟同步机方法为系统提供虚拟惯性和阻尼，但却忽视
了新能源机组的随机性可能造成系统阻尼特性随着复杂运行工况的变化而变化，故该控制方案与多变运行方式之间存在不适应性。


技术实现要素：

6.本发明的目的是：提供一种分布式新能源互联微电网集群振荡控制方法，以解决现有技术中存在的新能源机组的随机性与间歇性对于系统存在影响，基于附加控制方法在扰动较大或存在不确定性导致运行工况发生改变时，其设计的振荡抑制控制器难以保证抑制效果，缺乏工况适应性的问题。
7.为了实现上述目的，本发明提供了一种分布式新能源互联微电网集群振荡控制方法，包括：
8.s1、建立多源互联微电网系统的数学模型，基于多源互联微电网系统的稳态运行点，对多源互联微电网系统进行线性化处理，获得多源互联微电网系统的线性化小信号状态空间模型，所述多源互联微电网系统的线性化小信号状态空间模型包括系统状态矩阵；
9.s2、计算系统状态矩阵的特征值及左、右特征向量，获得右特征向量模态矩阵和左特征向量模态矩阵；
10.s3、选取待选输入信号和待选反馈输出信号，构建状态空间描述的多源互联微电网系统线性化模型，选择特定振荡模式，利用几何测度理论结合右特征向量模态矩阵和左特征向量模态矩阵进行能观可控量化分析，获得可控性和能观性几何测度指标；
11.s4、依据可控性和能观性几何测度指标，构建联合能观可控量化指标，根据联合能观可控量化指标设计抑制振荡模式的最优控制回路，所述最优控制回路包括最优控制输入信号和最优反馈输出信号；
12.s5、根据最优控制输入信号建立机组预测模型，并进行离散化处理，获得离散状态空间模型；
13.s6、依据离散状态空间模型计算多步的状态信息与输出信息的偏差，将最优反馈输出信号作为成本函数的预测追踪目标，结合状态信息与输出信息的偏差和最优控制输入信号构造二次规划形式的成本函数；
14.s7、求解二次规划形式的成本函数，如果无解，将控制量置零，并重复步骤s6，如果有解，获得控制输入增量，所述控制输入增量为最优控制序列。
15.优选的，在步骤s1中，所述多源互联微电网系统的数学模型为：
[0016][0017]
其中，x＝[x
sg1
,x
wt1
,x
pv1
,x
ess1
…
x
sgi
,x
wti
,x
pvi
,x
essi
,x
line
,x
load
]
t
为互联微电网系统中各微电源、线路及负荷的状态变量，n表示为系统状态量的个数，u为系统的输入变量，f表示为描述系统状态量的微分量与系统状态变量和输入变量关系的函数。
[0018]
优选的，所述多源互联微电网系统的数学模型包括同步机的动态模型、新能源机组的一次能源动态模型和逆变器的三环控制模型；
[0019]
所述同步机的动态模型包括同步机机组主电路模型和同步机机组控制电路模型，所述同步机机组主电路模型为：
[0020][0021][0022]
其中，δ和ωr为sbdg的功角和转子角速度；m和d为转子惯量系数及阻尼系数；pe和qe分别为sbdg输出的有功功率和无功功率；pm为同步发电机机械功率；ω
com
为互联微电网系统公共坐标轴参考角速度；vq和vq'为交轴空载电动势和暂态电动势；v
fd
和vf分别为励磁电压和下垂励磁环节输出电压；xd和x'd分别为额定频率ωn下发电机的直轴电抗和直轴暂态电抗；td'o为转子时间常数；v
t
＝[v
td v
tq
]
t
表示为sbdg端口电压；i
t
＝[i
td i
tq
]
t
表示为输出电流；
[0023]
所述同步机机组控制电路模型为：
[0024][0025][0026]
其中，v
exc
、k、τ分别为延迟环节输出电压、增益系数和延时系数；vn、ωn为系统额度电压和角速度；v
ref
、ω
ref
为电压和角速度参考值；tg、ψ1/ψ2/ψ3/ψ4均为速度中间控制器状态量；k1和t1/t2为调节器增益及其时间常数；t3/t4/t5为制动器时间常数；td为同步发电机时间延迟；
[0027]
针对基于逆变器类分布式电源，新能源机组的一次能源动态模型包括风力发电机组的风机及其传动系统模型和光伏发电机组的一次能源动态与直流电容电压动态模型，所述风机及其传动系统模型为：
[0028][0029][0030]
其中，ω
wt
、ωr、h
wt
、hg分别为风力机和异步电机转速和惯量系数；kg、f、k
tg
、d
tg
、t
sh
分别表示为齿轮箱传动比、摩擦系数、机械刚度、阻尼系数和机械轴转矩；t
wt
、te分别为风力机转矩和电磁转矩；ρ、r、v、p、c
p
分别表示为空气密度、叶片长度、风速、极对数和风能利用系数；表示为定子磁通；is＝[i
sd i
sq
]
t
表示为定子电流；
[0031]
所述一次能源动态与直流电容电压动态模型为：
[0032][0033]
其中，v
dc
、i
pv
、c
dc
分别为光伏直流侧电压、电流和电容；n
p
、ns、nc分别是光伏系统串联数量、每串串联的模块数量和每个光伏组件的串联电池数量；i
rs
和i
ph
分别为p-n结反向饱和电流和每串pv板的短路电流；q、k、t
t
、a
if
分别为单位电荷、玻尔兹曼常数、温度及理想因子；
[0034]
所述逆变器的三环控制模型为：
[0035][0036][0037][0038][0039]
其中，v
od
、v
oq
、i
od
、i
oq
分别为dq参考系下的逆变器输出电压和电流；m
p
、nq为下垂环输出参考电压、有功下垂系数和无功下垂系数；ψd、ψq为电压控制器积分项的状态变量；γd、γq为电流控制器积分项的状态变量；为逆变器出口电流参考值；i
ldq
为逆变器出口实际电流。
[0040]
优选的，在步骤s1中，所述多源互联微电网系统的线性化小信号状态空间模型为：
[0041][0042]
其中，a
nmg
是系统状态矩阵。
[0043]
优选的，在步骤s2中，所述计算系统状态矩阵的特征值及左右特征向量的公式为：
[0044][0045]
其中，分别是对应λi的右特征向量和左特征向量，对于同一特征值归一化处理后有而对于不同特征值的左右特征向量则满足关系式不同特征值λ1，λ2，...λn对应的右特征向量模态矩阵φ和左特征向量模态矩阵表示为：
[0046][0047]
优选的，在步骤s3中，所述选取待选输入信号和待选反馈输出信号的方法包括工程实际选取或参与因子分析结果进行选取，所述参与因子为状态量对于某一振荡模态的影响力或者贡献度的大小。
[0048]
优选的，在步骤s3中，所述状态空间描述的多源互联微电网系统线性化模型为：
[0049][0050]
优选的，在步骤s3中，所述获得可控性和能观性几何测度指标的公式为：
[0051][0052]
其中，bk为输入矩阵b
nmg
的第k列，c
l
为输出矩阵c
nmg
的第l行，为输入向量bk与左特征向量之间的夹角，为输出向量c
l
与右特征向量φi之间的夹角；
[0053]
在步骤s4中，所述联合能观可控量化指标为：
[0054]gcoi
(k,l)＝g
ci
(k)
·goi
(l)。
[0055]
优选的，在步骤s5中，所述离散状态空间模型为：
[0056][0057]
其中，ad、bd均为离散化处理后所获得的离散状态矩阵及离散输入矩阵；δu
*
(k)为选取的针对该特定振荡模式且不影响其它振荡抑制作用的最优控制输入信号。
[0058]
优选的，在步骤s6中，二次规划形式的成本函数为：
[0059]
[0060]
subject to aiδui≤b
i a
eq,i
δui＝b
eq,i
[0061]
与现有技术相比，本发明的有益效果在于：对于多源互联微电网系统的振荡问题，采用了几何测度理论，分析了待选输入信号及待选反馈输出信号的综合可控性以及能观性，获取抑制特定振荡模式的最优控制回路，能够避免控制器之间的作用冲突；考虑到系统多种约束情况下构成二次规划形式的成本函数，形成抑制多源互联微电网系统的振荡的稳定控制策略，能够在线实现振荡抑制作用，更具工况适应性；
[0062]
进一步地，考虑了基于同步机类分布式电源和基于逆变器类分布式电源这两类微电源，同时对于新能源发电机组考虑到供能部分动态和机侧换流器控制特性，进而能够减少新能源机组的随机性与间歇性对于系统的影响。
附图说明
[0063]
图1是本发明实施例的分布式新能源互联微电网集群振荡控制方法的流程图；
[0064]
图2是本发明实施例的多源互联微电网系统的结构示意图；
[0065]
图3是本发明实施例的基于逆变器类分布式电源振荡抑制控制器结构示意图；
[0066]
图4是本发明实施例的基于同步机类分布式电源振荡抑制控制器结构示意图。
具体实施方式
[0067]
下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
[0068]
在本发明的描述中，需要说明的是，术语“中心”、“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”、“第三”仅用于描述目的，而不能理解为指示或暗示相对重要性。
[0069]
在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
[0070]
此外，在本发明的描述中，除非另有说明，“多个”的含义是两个或两个以上。
[0071]
实施例一
[0072]
如图1-2所示，本发明优选实施例的一种分布式新能源互联微电网集群振荡控制方法，包括：
[0073]
s1、建立多源互联微电网系统的数学模型，基于多源互联微电网系统的稳态运行点，对多源互联微电网系统进行线性化处理，获得多源互联微电网系统的线性化小信号状态空间模型，所述多源互联微电网系统的线性化小信号状态空间模型包括系统状态矩阵；
[0074]
s2、计算系统状态矩阵的特征值及左、右特征向量，获得右特征向量模态矩阵和左特征向量模态矩阵；
[0075]
s3、选取待选输入信号和待选反馈输出信号，构建状态空间描述的多源互联微电网系统线性化模型，选择特定振荡模式，利用几何测度理论结合右特征向量模态矩阵和左特征向量模态矩阵进行能观可控量化分析，获得可控性和能观性几何测度指标；
[0076]
s4、依据可控性和能观性几何测度指标，构建联合能观可控量化指标，根据联合能观可控量化指标设计抑制振荡模式的最优控制回路，所述最优控制回路包括最优控制输入信号和最优反馈输出信号；
[0077]
s5、根据最优控制输入信号建立机组预测模型，并进行离散化处理，获得离散状态空间模型；
[0078]
s6、依据离散状态空间模型计算多步的状态信息与输出信息的偏差，将最优反馈输出信号作为成本函数的预测追踪目标，结合状态信息与输出信息的偏差和最优控制输入信号构造二次规划形式的成本函数；
[0079]
s7、求解二次规划形式的成本函数，如果无解，将控制量置零，并重复步骤s6，如果有解，获得控制输入增量，所述控制输入增量为最优控制序列。
[0080]
实施例二
[0081]
多源互联微电网系统主要由微电源、负载、电力网络和网间换流器等元件组成，其中的微电源大致可分为基于同步机类分布式电源(synchronous-based distributed generator,sbdg)和基于逆变器类分布式电源(inverter-based distributed generator,ibdg)。
[0082]
如图3-4所示，在步骤s1中，储能微电源作为多源互联微电网系统的重要组成部分，可以通过控制自身充放电来调节微电网系统中的功率，平抑系统的功率波动，为系统提供频率与电压支撑。基于以上考虑机组特性的微型燃气轮机、风力发电机组和光伏机组的数学模型，各分布式发电单元主要通过电力网络连接形成互联微电网系统，一般可以用n个一阶非线性常微分方程来刻画多源互联微电网系统的运行特性，多源互联微电网系统的数学模型为：
[0083][0084]
其中，x＝[x
sg1
,x
wt1
,x
pv1
,x
ess1
…
x
sgi
,x
wti
,x
pvi
,x
essi
,x
line
,x
load
]
t
为互联微电网系统中各微电源、线路及负荷的状态变量，n表示为系统状态量的个数，u为系统的输入变量，f表示为描述系统状态量的微分量与系统状态变量和输入变量关系的函数。
[0085]
多源互联微电网系统的数学模型包括同步机的动态模型、新能源机组的一次能源动态模型和逆变器的三环控制模型；
[0086]
对于常见的微型燃气轮机发电机组可视为受下垂励磁电压控制环节和速度控制环节的同步机，根据四阶同步机转子运动方程和绕组电压方程，同步机的动态模型包括同步机机组主电路模型和同步机机组控制电路模型，同步机机组主电路模型为：
[0087]
[0088][0089]
其中，δ和ωr为sbdg的功角和转子角速度；m和d为转子惯量系数及阻尼系数；pe和qe分别为sbdg输出的有功功率和无功功率；pm为同步发电机机械功率；ω
com
为互联微电网系统公共坐标轴参考角速度；vq和vq'为交轴空载电动势和暂态电动势；v
fd
和vf分别为励磁电压和下垂励磁环节输出电压；xd和x'd分别为额定频率ωn下发电机的直轴电抗和直轴暂态电抗；td'o为转子时间常数；v
t
＝[v
td v
tq
]
t
表示为sbdg端口电压；i
t
＝[i
td i
tq
]
t
表示为输出电流；
[0090]
对于基于同步机类分布式电源，其励磁环节无功/电压下垂控制环节和速度控制环节的有功/频率下垂控制，同步机机组控制电路模型为：
[0091][0092][0093]
其中，v
exc
、k、τ分别为延迟环节输出电压、增益系数和延时系数；vn、ωn为系统额度电压和角速度；v
ref
、ω
ref
为电压和角速度参考值；tg、ψ1/ψ2/ψ3/ψ4均为速度中间控制器状态量；k1和t1/t2为调节器增益及其时间常数；t3/t4/t5为制动器时间常数；td为同步发电机时间延迟；
[0094]
针对基于逆变器类分布式电源，通常在忽视供能部分动态或机侧换流器控制特性时，通常以理想储能发电单元为载体，将供能部分处理为恒定电压源或电流源，但是风光新能源机组的随机性与间歇性会严重影响系统的动态响应，因此需要保证建模的完整性，考虑这些机组的一次能源动态与直流环节动态的建模。
[0095]
新能源机组的一次能源动态模型包括风力发电机组的风机及其传动系统模型和光伏发电机组的一次能源动态与直流电容电压动态模型，风机及其传动系统模型为：
[0096][0097][0098]
其中，ω
wt
、ωr、h
wt
、hg分别为风力机和异步电机转速和惯量系数；kg、f、k
tg
、d
tg
、t
sh
分别表示为齿轮箱传动比、摩擦系数、机械刚度、阻尼系数和机械轴转矩；t
wt
、te分别为风力机转矩和电磁转矩；ρ、r、v、p、c
p
分别表示为空气密度、叶片长度、风速、极对数和风能利用系数；表示为定子磁通；is＝[i
sd i
sq
]
t
表示为定子电流；
[0099]
对于光伏发电机组，在忽略换流器损耗时，一次能源动态与直流电容电压动态模型为：
[0100][0101]
其中，v
dc
、i
pv
、c
dc
分别为光伏直流侧电压、电流和电容；n
p
、ns、nc分别是光伏系统串联数量、每串串联的模块数量和每个光伏组件的串联电池数量；i
rs
和i
ph
分别为p-n结反向饱和电流和每串pv板的短路电流；q、k、t
t
、a
if
分别为单位电荷、玻尔兹曼常数、温度及理想因子；
[0102]
针对基于逆变器类分布式电源(ibdg)的控制回路部分，主要包括瞬时功率计算、下垂控制环路和内部电压电流控制环路，逆变器的三环控制模型为：
[0103][0104][0105][0106][0107]
其中，v
od
、v
oq
、i
od
、i
oq
分别为dq参考系下的逆变器输出电压和电流；m
p
、nq为下垂环输出参考电压、有功下垂系数和无功下垂系数；ψd、ψq为电压控制器积分项的状态变量；
γd、γq为电流控制器积分项的状态变量；为逆变器出口电流参考值；i
ldq
为逆变器出口实际电流。
[0108]
在步骤s1中，基于稳态运行点x0，可对系统进行线性化，并进行泰勒展开保留其展开各项中的常数项和一阶项，多源互联微电网系统的线性化小信号状态空间模型为：
[0109][0110]
其中，a
nmg
是系统状态矩阵。
[0111]
在步骤s2中，在获取多源互联微电网系统的动态模型后，为了研究多源互联微电网系统的振荡稳定性，还需要进一步详细分析系统状态矩阵的特征值以获得系统主导振荡模态、振荡频率与阻尼比、引起振荡的变量参与度分析等。基于系统线性化可得到系统状态矩阵a
nmg
的具体数值表示，计算系统状态矩阵的特征值及左右特征向量的公式为：
[0112][0113]
其中，分别是对应λi的右特征向量和左特征向量，对于同一特征值归一化处理后有而对于不同特征值的左右特征向量则满足关系式不同特征值λ1，λ2，...λn对应的右特征向量模态矩阵φ和左特征向量模态矩阵表示为：
[0114][0115]
实施例三
[0116]
由于产生振荡的位置和最佳阻尼控制位置常常不在同一地点，因此在进行具体的控制器设计之前，首先需要解决控制器的安装位置、控制输入-输出的选择和配对等控制结构问题，利用几何测度法对比不同类型信号的可控性和能观性优劣，用于设计抑制不同振荡模式下的最优控制回路。在进行几何测度分析之前，事先确定待选输入信号u(t)以及待选反馈输出信号y(t)，构成含有控制输入及系统输出的状态空间描述的多源互联微电网系统线性化模型。
[0117]
在步骤s3中，选取待选输入信号和待选反馈输出信号的方法包括工程实际选取或参与因子分析结果进行选取，参与因子为状态量对于某一振荡模态的影响力或者贡献度的大小。
[0118]
在步骤s3中，状态空间描述的多源互联微电网系统线性化模型为：
[0119][0120]
在步骤s3中，获得可控性和能观性几何测度指标的公式为：
[0121][0122]
其中，bk为输入矩阵b
nmg
的第k列，c
l
为输出矩阵c
nmg
的第l行，为输入向量bk与左特征向量之间的夹角，为输出向量c
l
与右特征向量φi之间的夹角；
[0123]
在步骤s4中，依据各振荡模式的可控性和能观性几何测度指标g
ci
、g
oi
，进而可以定义振荡模式λi的联合能观可控量化指标为：
[0124]gcoi
(k,l)＝g
ci
(k)
·goi
(l)。
[0125]
抑制振荡控制回路设计依据为首要关注弱阻尼振荡模式有最大联合能观可控量化指标g
coi
(k,l)，同时需要保证待选信号对各其他振荡模式的影响作用最小，从而可以获取抑制该振荡模式最优控制回路，同时还可以避免控制器之间的交互作用。
[0126]
在步骤s5中，根据选取的抑制特定振荡模式的最优控制输入位置，建立仅含该控制输入的机组预测模型，以采样周期ts将机组线性化小信号状态空间模型进行离散化处理，离散状态空间模型为：
[0127][0128]
其中，ad、bd均为离散化处理后所获得的离散状态矩阵及离散输入矩阵；δu
*
(k)为选取的针对该特定振荡模式且不影响其它振荡抑制作用的最优控制输入信号。
[0129]
在步骤s6中，根据选择的抑制特定振荡模式的最优控制输入信号，将其作为成本函数的预测追踪目标，并在考虑各机组等式及不等式约束的条件下，二次规划形式的成本函数为：
[0130][0131]
subject to aiδui≤b
i a
eq,i
δui＝b
eq,i
[0132]
其中，求取所得到成本函数为最小值时对应的控制输入增量δui即为抑制该振荡模式下的最优控制序列。
[0133]
本发明的工作过程为：
[0134]
s1、建立多源互联微电网系统的数学模型，基于多源互联微电网系统的稳态运行点，对多源互联微电网系统进行线性化处理，获得多源互联微电网系统的线性化小信号状态空间模型，所述多源互联微电网系统的线性化小信号状态空间模型包括系统状态矩阵；
[0135]
s2、计算系统状态矩阵的特征值及左、右特征向量，获得右特征向量模态矩阵和左特征向量模态矩阵；
[0136]
s3、选取待选输入信号和待选反馈输出信号，构建状态空间描述的多源互联微电网系统线性化模型，选择特定振荡模式，利用几何测度理论结合右特征向量模态矩阵和左特征向量模态矩阵进行能观可控量化分析，获得可控性和能观性几何测度指标；
[0137]
s4、依据可控性和能观性几何测度指标，构建联合能观可控量化指标，根据联合能
观可控量化指标设计抑制振荡模式的最优控制回路，所述最优控制回路包括最优控制输入信号和最优反馈输出信号；
[0138]
s5、根据最优控制输入信号建立机组预测模型，并进行离散化处理，获得离散状态空间模型；
[0139]
s6、依据离散状态空间模型计算多步的状态信息与输出信息的偏差，将最优反馈输出信号作为成本函数的预测追踪目标，结合状态信息与输出信息的偏差和最优控制输入信号构造二次规划形式的成本函数；
[0140]
s7、求解二次规划形式的成本函数，如果无解，将控制量置零，并重复步骤s6，如果有解，获得控制输入增量，所述控制输入增量为最优控制序列。
[0141]
综上，本发明实施例提供一种分布式新能源互联微电网集群振荡控制方法，其对于多源互联微电网系统的振荡问题，采用了几何测度理论，分析了待选输入信号及待选反馈输出信号的综合可控性以及能观性，获取抑制特定振荡模式的最优控制回路，能够避免控制器之间的作用冲突；考虑到系统多种约束情况下构成二次规划形式的成本函数，形成抑制多源互联微电网系统的振荡的稳定控制策略，能够在线实现振荡抑制作用，更具工况适应性；考虑了基于同步机类分布式电源和基于逆变器类分布式电源这两类微电源，同时对于新能源发电机组考虑到供能部分动态和机侧换流器控制特性，能够减少新能源机组的随机性与间歇性对于系统的影响。
[0142]
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和替换，这些改进和替换也应视为本发明的保护范围。

技术特征：
1.一种分布式新能源互联微电网集群振荡控制方法，其特征在于，包括：s1、建立多源互联微电网系统的数学模型，基于多源互联微电网系统的稳态运行点，对多源互联微电网系统进行线性化处理，获得多源互联微电网系统的线性化小信号状态空间模型，所述多源互联微电网系统的线性化小信号状态空间模型包括系统状态矩阵；s2、计算系统状态矩阵的特征值及左、右特征向量，获得右特征向量模态矩阵和左特征向量模态矩阵；s3、选取待选输入信号和待选反馈输出信号，构建状态空间描述的多源互联微电网系统线性化模型，选择特定振荡模式，利用几何测度理论结合右特征向量模态矩阵和左特征向量模态矩阵进行能观可控量化分析，获得可控性和能观性几何测度指标；s4、依据可控性和能观性几何测度指标，构建联合能观可控量化指标，根据联合能观可控量化指标设计抑制振荡模式的最优控制回路，所述最优控制回路包括最优控制输入信号和最优反馈输出信号；s5、根据最优控制输入信号建立机组预测模型，并进行离散化处理，获得离散状态空间模型；s6、依据离散状态空间模型计算多步的状态信息与输出信息的偏差，将最优反馈输出信号作为成本函数的预测追踪目标，结合状态信息与输出信息的偏差和最优控制输入信号构造二次规划形式的成本函数；s7、求解二次规划形式的成本函数，如果无解，将控制量置零，并重复步骤s6，如果有解，获得控制输入增量，所述控制输入增量为最优控制序列。2.根据权利要求1所述的一种分布式新能源互联微电网集群振荡控制方法，其特征在于：在步骤s1中，所述多源互联微电网系统的数学模型为：其中，x＝[x
sg1
,x
wt1
,x
pv1
,x
ess1
…
x
sgi
,x
wti
,x
pvi
,x
essi
,x
line
,x
load
]
t
为互联微电网系统中各微电源、线路及负荷的状态变量，n表示为系统状态量的个数，u为系统的输入变量，f表示为描述系统状态量的微分量与系统状态变量和输入变量关系的函数。3.根据权利要求2所述的一种分布式新能源互联微电网集群振荡控制方法，其特征在于：所述多源互联微电网系统的数学模型包括同步机的动态模型、新能源机组的一次能源动态模型和逆变器的三环控制模型；所述同步机的动态模型包括同步机机组主电路模型和同步机机组控制电路模型，所述同步机机组主电路模型为：
其中，δ和ω
r
为sbdg的功角和转子角速度；m和d为转子惯量系数及阻尼系数；p
e
和q
e
分别为sbdg输出的有功功率和无功功率；p
m
为同步发电机机械功率；ω
com
为互联微电网系统公共坐标轴参考角速度；v
q
和v
′
q
为交轴空载电动势和暂态电动势；v
fd
和v
f
分别为励磁电压和下垂励磁环节输出电压；x
d
和x
′
d
分别为额定频率ω
n
下发电机的直轴电抗和直轴暂态电抗；为转子时间常数；v
t
＝[v
td v
tq
]
t
表示为sbdg端口电压；i
t
＝[i
td i
tq
]
t
表示为输出电流；所述同步机机组控制电路模型为：所述同步机机组控制电路模型为：其中，v
exc
、k、τ分别为延迟环节输出电压、增益系数和延时系数；v
n
、ω
n
为系统额度电压和角速度；v
ref
、ω
ref
为电压和角速度参考值；t
g
、ψ1/ψ2/ψ3/ψ4均为速度中间控制器状态量；k1和t1/t2为调节器增益及其时间常数；t3/t4/t5为制动器时间常数；t
d
为同步发电机时间延迟；针对基于逆变器类分布式电源，新能源机组的一次能源动态模型包括风力发电机组的风机及其传动系统模型和光伏发电机组的一次能源动态与直流电容电压动态模型，所述风机及其传动系统模型为：
其中，ω
wt
、ω
r
、h
wt
、h
g
分别为风力机和异步电机转速和惯量系数；k
g
、f、k
tg
、d
tg
、t
sh
分别表示为齿轮箱传动比、摩擦系数、机械刚度、阻尼系数和机械轴转矩；t
wt
、t
e
分别为风力机转矩和电磁转矩；ρ、r、v、p、c
p
分别表示为空气密度、叶片长度、风速、极对数和风能利用系数；表示为定子磁通；i
s
＝[i
sd i
sq
]
t
表示为定子电流；所述一次能源动态与直流电容电压动态模型为：其中，v
dc
、i
pv
、c
dc
分别为光伏直流侧电压、电流和电容；n
p
、n
s
、n
c
分别是光伏系统串联数量、每串串联的模块数量和每个光伏组件的串联电池数量；i
rs
和i
ph
分别为p-n结反向饱和电流和每串pv板的短路电流；q、k、t
t
、a
if
分别为单位电荷、玻尔兹曼常数、温度及理想因子；所述逆变器的三环控制模型为：所述逆变器的三环控制模型为：所述逆变器的三环控制模型为：所述逆变器的三环控制模型为：其中，v
od
、v
oq
、i
od
、i
oq
分别为dq参考系下的逆变器输出电压和电流；m
p
、n
q
为下垂环输出参考电压、有功下垂系数和无功下垂系数；ψ
d
、ψ
q
为电压控制器积分项的状态变量；γ
d
、
γ
q
为电流控制器积分项的状态变量；为逆变器出口电流参考值；i
ldq
为逆变器出口实际电流。4.根据权利要求1所述的一种分布式新能源互联微电网集群振荡控制方法，其特征在于：在步骤s1中，所述多源互联微电网系统的线性化小信号状态空间模型为：其中，a
nmg
是系统状态矩阵。5.根据权利要求1所述的一种分布式新能源互联微电网集群振荡控制方法，其特征在于：在步骤s2中，所述计算系统状态矩阵的特征值及左右特征向量的公式为：其中，φ
i
,分别是对应λ
i
的右特征向量和左特征向量，对于同一特征值归一化处理后有而对于不同特征值的左右特征向量则满足关系式不同特征值λ1，λ2，...λ
n
对应的右特征向量模态矩阵φ和左特征向量模态矩阵表示为：6.根据权利要求1所述的一种分布式新能源互联微电网集群振荡控制方法，其特征在于：在步骤s3中，所述选取待选输入信号和待选反馈输出信号的方法包括工程实际选取或参与因子分析结果进行选取，所述参与因子为状态量对于某一振荡模态的影响力或者贡献度的大小。7.根据权利要求1所述的一种分布式新能源互联微电网集群振荡控制方法，其特征在于：在步骤s3中，所述状态空间描述的多源互联微电网系统线性化模型为：8.根据权利要求1所述的一种分布式新能源互联微电网集群振荡控制方法，其特征在于：在步骤s3中，所述获得可控性和能观性几何测度指标的公式为：其中，b
k
为输入矩阵b
nmg
的第k列，c
l
为输出矩阵c
nmg
的第l行，为输入向量b
k
与左
特征向量之间的夹角，为输出向量c
l
与右特征向量φ
i
之间的夹角；在步骤s4中，所述联合能观可控量化指标为：g
coi
(k,l)＝g
ci
(k)
·
g
oi
(l)。9.根据权利要求1所述的一种分布式新能源互联微电网集群振荡控制方法，其特征在于：在步骤s5中，所述离散状态空间模型为：其中，a
d
、b
d
均为离散化处理后所获得的离散状态矩阵及离散输入矩阵；δu
*
(k)为选取的针对该特定振荡模式且不影响其它振荡抑制作用的最优控制输入信号。10.根据权利要求1所述的一种分布式新能源互联微电网集群振荡控制方法，其特征在于：在步骤s6中，二次规划形式的成本函数为：minmizesubject to a
i
δu
i
≤b
i a
eq,i
δu
i
＝b
eq,i
。

技术总结
本发明涉及多源互联微电网系统稳定控制技术领域，公开了一种分布式新能源互联微电网集群振荡控制方法，对于多源互联微电网系统的振荡问题，采用了几何测度理论，分析了待选输入信号及待选反馈输出信号的综合可控性以及能观性，获取抑制特定振荡模式的最优控制回路，能够避免控制器之间的作用冲突；考虑到系统多种约束情况下构成二次规划形式的成本函数，形成抑制多源互联微电网系统的振荡的稳定控制策略，能够在线实现振荡抑制作用，更具工况适应性；考虑了基于同步机类分布式电源和基于逆变器类分布式电源这两类微电源，同时对于新能源发电机组考虑到供能部分动态和机侧换流器控制特性，能够减少新能源机组的随机性与间歇性对于系统的影响。间歇性对于系统的影响。间歇性对于系统的影响。


技术研发人员：赵卓立 杨庆刚 许志荣 徐家文 卢健钊 孟安波
受保护的技术使用者：广东工业大学
技术研发日：2023.06.13
技术公布日：2023/9/22