一种新能源发电振荡抑制方法与系统与流程

1.本发明涉及发电振荡抑制领域，具体涉及一种新能源发电振荡抑制方法与系统。


背景技术：

2.在中国的北部和西北部，大多数大型新能源发电都在偏远地区开发，并通过长距离输电线路接入主电网，新能源接入点电网强度弱。新能源发电控制和弱电网之间的复杂相互作用可能导致不稳定和振荡问题。例如，在2015年至2017年间，新疆哈密风电基地频繁出现次/超同步振荡问题，次同步振荡分量与超同步振荡分量关于基频50hz对称，如图1所示。此宽频段范围内的振荡分量易与同步发电机组轴系扭振相互作用，导致轴系永久性损伤。
3.近几年，上述问题成为国内外研究的热点。研究表明，新能源并网振荡产生的主要因素为新能源发电阻抗与电网阻抗形成了串联振荡回路，振荡回路在振荡频率点的阻尼不足，易产生振荡问题。
4.目前，实际系统对新能源并网振荡的防控措施主要有两方面，一是快速切除新能源场站，当振荡产生时，新能源场站的控制保护装置动作，直接跳开新能源场站送出线路，相当于破坏了系统串联振荡的条件，进而消除振荡；二是改进新能源发电机组控制，通过优化新能源发电机组控制参数，改变新能源发电的阻抗特性，例如，提升振荡频率点的阻尼特性等，提升系统的稳定裕度，离线刷入实际机组控制器后，进行现场试验验证。
5.上述直接切除场站和离线优化控制两种防控手段均存在不足。直接切除新能源场站使得场站整体断电，不利于新能源的送出和消纳，同时，可能对新能源机组造成损坏。离线优化新能源发电机组控制的方法在并网条件和运行工况发生变化情况下适应性不足，应对新的振荡问题工程量大。


技术实现要素：

6.为了解决现有技术中所存在的问题，本发明提供一种新能源发电振荡抑制方法，包括：
7.对于采集的机组端口电压、电流瞬时量，利用在线fft分析确定各频段的振荡分量；
8.对各频段的振荡分量进行振荡识别；
9.基于振荡识别结果按振荡频段和所述频段对应的触发条件启动控制参数的在线调整，实现振荡抑制。
10.优选的，所述对各频段的振荡分量进行振荡识别，包括：
11.基于不同频段分别按照振荡电流分量幅值排序，识别出各频段的主导振荡分量。
12.优选的，所述触发条件包括：基于振荡电流幅值触发和基于振荡电流幅值结合振荡持续时间触发。
13.优选的，所述基于振荡电流幅值触发，包括：
14.当任意频率振荡电流幅值达到第一阈值时，直接触发抑制控制。
15.优选的，所述基于振荡电流幅值结合振荡持续时间触发，包括：
16.振荡电流幅值未达到第一阈值，但达到第二阈值时，触发计时器，当振荡持续时间达到时间阈值时，触发振荡抑制控制；
17.其中，第二阈值小于第一阈值。
18.优选的，所述基于振荡识别结果按振荡频段和所述频段对应的触发条件启动控制参数的在线调整，包括：
19.按照不同短路比的适应性，针对各振荡频段分别按顺序设置至少三组控制参数；
20.当收到控制触发信号时，按照触发振荡告警信号基于所述顺序依次递增在线调整控制参数直到振荡抑制成功或触发停机指令；
21.所述各振荡频段的控制参数的设置顺序不同。
22.优选的，所述振荡频段包括：次/超同步频段、中频段和高频段；
23.其中，所述次/超同步频段为1～100hz范围内；所述中频段为100hz～f
mb
范围内；所述高频段为大于f
mb
；所述f
mb
为电流环带宽频率。
24.优选的，在采集机组端口电压、电流瞬时量之前，还包括：
25.为所述次/超同步频段设定不小于1khz的采样频率；
26.为所述高频段设定大于5khz的采样频率；
27.为所述中频段设定1khz到5khz之间的任意一个采样频率；
28.设定频谱分析时间窗长。
29.基于同一种发明构思，本发明还提供一种新能源发电振荡抑制系统，包括：
30.采集模块，用于对于采集的机组端口电压、电流瞬时量，利用在线fft分析确定各频段的振荡分量；
31.识别模块，用于对各频段的振荡分量进行振荡识别；
32.抑制模块，用于基于振荡识别结果按振荡频段和所述频段对应的触发条件启动控制参数的在线调整，实现振荡抑制。
33.优选的，所述识别模块具体用于：基于不同频段分别按照振荡电流分量幅值排序，识别出各频段的主导振荡分量。
34.优选的，所述抑制模块具体用于：按照不同短路比的适应性，针对各振荡频段分别按顺序设置至少三组控制参数；当收到控制触发信号时，按照触发振荡告警信号基于所述顺序依次递增在线调整控制参数直到振荡抑制成功或触发停机指令；所述各振荡频段的控制参数的设置顺序不同。
35.与现有技术相比，本发明的有益效果为：
36.本发明提出一种新能源发电振荡抑制方法与系统，包括对于采集的机组端口电压、电流瞬时量，利用在线fft分析确定各频段的振荡分量；对各频段的振荡分量进行振荡识别；基于振荡识别结果按振荡频段和所述频段对应的触发条件启动控制参数的在线调整，实现实时进行控制参数的调整优化和风电机组在弱电网下宽频带振荡问题的有效抑制；
37.本发明提供的技术手段对不同类型的新能源发电装置都适用，具有通用性，有利于解决实际各种类型新能源发电并网振荡难题；
38.本发明提供的技术手段是解决多频段振荡问题的通法，但针对不同频段的振荡问题应采用不同的参数和策略；
39.本发明设置具有两种触发阈值以及振荡幅值触发和振荡时间触发两种方式的控制触发机制；
40.本发明控制参数在线调整方法中设置多组控制参数，以便于针对不同振荡问题进行多样化的策略调整。
附图说明
41.图1为哈密某风电机组端口电流振荡分量傅里叶分析图；
42.图2为本发明的一种新能源发电振荡抑制方法流程图；
43.图3为本发明的振荡分量实时监测技术框图；
44.图4为本发明的次/超同步振荡抑制触发控制框图；
45.图5为本发明的次/超同步振荡控制参数在线调整框图；
46.图6为本发明的基于在线阻抗重塑的风电机组宽频带振荡抑制技术框图；
47.图7为本发明的基于在线阻抗重塑的直驱风电机组宽频带振荡抑制仿真验证。
具体实施方式
48.本发明的一种新能源发电振荡抑制方法，实时采集机组端口电压、电流瞬时量，通过在线fft分析，得出次/超同步、中频段、高频段等不同频段的振荡分量，通过振荡识别和判定触发告警信号，启动控制参数的在线调整，对风电机组锁相环、电流环、直流电压环等多回路控制参数进行调整，实现振荡抑制。为了更好地理解本发明，下面结合说明书附图和实例对本发明的内容做进一步的说明。
49.实施例1：
50.本发明提出一种新能源发电振荡抑制方法，如图2所示，包括：
51.s1、对于采集的机组端口电压、电流瞬时量，利用在线fft分析确定各频段的振荡分量；
52.s2、对各频段的振荡分量进行振荡识别；
53.s3、基于振荡识别结果按振荡频段和所述频段对应的触发条件启动控制参数的在线调整，实现振荡抑制。
54.下面，详述本发明一种新能源发电振荡抑制方法，其过程如图3所示。
55.步骤s1、对于采集的机组端口电压、电流瞬时量，利用在线fft分析确定各频段的振荡分量；具体包括：
56.为了进行宽频带振荡的实时监测，进行电压电流采样频率及频谱分析时间窗长的初始化设置。电压电流采样频率的设置与监测频率带宽范围有关，对于次/超同步振荡，1khz的采样频率能够满足振荡监测与控制的精度要求，而对于1khz以上的高频振荡，电压电流采样频率应达到5khz以上。频谱分析时间窗长的选择涉及频率分辨率和算法速度的要求，较长的数据窗能够达到更高精度的频率分辨率，同时带来的是计算时间的增加。对于风电机组的振荡抑制控制，更为重要的是分析和控制算法的执行速度，因此，频谱分析时间窗长设置为1秒。
57.步骤s2、对各频段的振荡分量进行振荡识别，包括：
58.对不同频段按照振荡电流分量幅值排序。由于风电机组并网系统不同频段振荡问题涉及的控制环节和参数不同，因此，在振荡分量实时监测环节，按照次/超同步频段、中频段和高频段不同频段分别按照振荡电流分量幅值排序，以便快速识别不同频段的主导振荡分量。次/超同步频段一般设为1～100hz范围，中频段一般设为100hz到风电机组电流环带宽频率左右，设为f
mb
，高频段为电流环带宽频率以上频率范围。
59.步骤s3、基于振荡识别结果按振荡频段和所述频段对应的触发条件启动控制参数的在线调整，包括：
60.当监测到不同频段振荡分量时，触发基于控制参数在线调整的振荡抑制控制。控制触发机制的核心是设置触发方式和条件，设置两种触发条件，如图4所示，一是当任意频率振荡电流幅值达到阈值1即i
m1,ssb
时，直接触发抑制控制；二是，振荡电流幅值未达到i
m1,ssb
，但达到阈值2即i
m2,ssb
时，触发计时器，当振荡持续时间达到t
m,ssb
时，触发振荡抑制控制，i
m2,ssb
＜i
m1,ssb
。通过设置上述振荡幅值触发和振荡时间触发两种方式，一方面可以减少运行过程中的扰动造成的误动作，另一方面，对于发散速度快的振荡，可实现快速响应和控制。由于不同的频段的振荡问题主导的控制环节和参数不同，故应分频段单独设定阈值等参数，单独触发。
61.依据前述风电机组控制参数优化方法，以及宽频带振荡的告警触发机制，提出风电机组次/超同步振荡控制参数在线调整方法。按照不同短路比的适应性，设置三组或多组控制参数，如图5所示，例如，和和和表示三组锁相环控制参数，和和和表示三组直流电压环控制参数，和和和表示三组电流环控制参数。当收到控制触发信号时，按照触发振荡告警信号依次递增在线调整控制参数。当第1次触发振荡告警信号时，将锁相环控制参数从和调整到和将电流环控制参数从和调整到和当第2次振荡告警信号触发时，将直流电压环控制参数从和调整到和当第3次振荡告警信号触发时，将锁相环控制参数从和调整到和将电流环控制参数从和调整到和当第4次振荡信号触发时，将直流电压环控制参数从和调整到和若振荡仍无法平息，则触发停机指令。由于不同频段的振荡问题主导的控制环节和参数不同，因此，不同频段的振荡控制可设置不同顺序的参数调整策略。对于次/超同步振荡，主导控制环节为锁相环和直流电压环，因此，针对次/超同步振荡，控制参数的调整策略以锁相环和直流电压环参数为主。
62.上述方法仅考虑次/超同步振荡下基于振荡分量识别及告警触发机制的振荡控制参数在线调整方法步骤。在实际应用中，由于工况的变化，振荡频率分布在各个频段内，不同频段的振荡问题主导因素也各不相同，故应在线进行振荡分量识别并针对振荡分量分布频段分别分析，设置不同顺序的参数调整策略，快速响应抑制振荡。
63.本实施例中控制参数的设置可以按照常规方式设置也可以按照实施例4的过程基于阻尼比灵敏度进行确定。
64.实施例2：
65.基于同一种发明构思，本发明还提供了一种新能源发电振荡抑制系统，包括：
66.采集模块，用于对于采集的机组端口电压、电流瞬时量，利用在线fft分析确定各频段的振荡分量；
67.识别模块，用于对各频段的振荡分量进行振荡识别；
68.抑制模块，用于基于振荡识别结果按振荡频段和所述频段对应的触发条件启动控制参数的在线调整，实现振荡抑制。
69.其中，识别模块具体用于：基于不同频段分别按照振荡电流分量幅值排序，识别出各频段的主导振荡分量。
70.抑制模块具体用于：按照不同短路比的适应性，针对各振荡频段分别按顺序设置至少三组控制参数；当收到控制触发信号时，按照触发振荡告警信号基于所述顺序依次递增在线调整控制参数直到振荡抑制成功或触发停机指令；所述各振荡频段的控制参数的设置顺序不同。
71.本实施例的一种新能源发电振荡抑制系统是为了实现本发明提供的一种新能源发电振荡抑制方法而设计的软件系统，因此本实施例中各模块实现的功能与实施例1一致，这里不再赘述。
72.实施例3
73.基于本发明的一种新能源发电振荡抑制方法，在实时仿真平台中搭建上图6所示系统的仿真模型，振荡监测与抑制控制模块采用c语言编程实现，直驱风电机组采用并网逆变器模型。
74.振荡抑制控制仿真结果如图7所示。在5s时，增加风电机组有功功率，振荡电流幅值达到i
m2,ssb
，触发计时器，达到设置的时间阈值t
m,ssb
＝1s时，控制触发，控制参数由初始参数组1调整到参数组2。随着功率的不断增大，11.3s时，再次达到振荡电流阈值，12.3s时控制触发，调整至参数组3，振荡逐步衰减，但是振荡累计时间再次达到t
m,ssb
，控制再次触发，调整至参数组4，振荡逐渐消除。
75.实施例4
76.基于阻尼比灵敏度确定风电机组控制参数之前需要，先构建新能源发电并网逆变器的阻抗模型，然后对新能源并网系统宽频带动态建模与稳定性分析，然后基于上述理论确定风电机组控制参数。
77.首先构建新能源发电并网逆变器的阻抗模型：
78.典型新能源并网逆变器的电路和控制结构主要包括直流母线电容，三相h桥开关电路，交流滤波电感等电路环节，以及锁相环、电流环、直流母线电压环、pmw调制等控制环节。
79.新能源发电侧等效为恒定的电流源id，直流母线电压控制的功能是控制直流母线电压保持在额定值v
dc
，并给出d轴电流控制的参考值i
dref
。锁相控制的功能是跟踪电网电压，得到相角θ
pll
，用于三相静止坐标系与dq旋转坐标系的变换。电流控制的功能是控制并网逆变器输出的电流跟踪指令值i
dref
和i
qref
。
80.直流母线电容动态、锁相环控制等电路和控制的非线性因素会引起逆变器的频率耦合效应，在三相静止坐标系下的具体表现为：在频率为f
p
的正序小信号电压扰动下，逆变
器除了会产生对应频率的正序电流响应外，还会产生频率为f
p-2f1的负序电流响应，反之同理。对于三相交流系统，当f
p
＜2f1时，负频率f
p-2f1下的负序分量在数学上与正频率2f
1-f
p
的正序分量的共轭等价。逆变器存在频率耦合效应，为了描述其扰动、响应特性，并网逆变器的序阻抗模型被推广为正、负序不同频率耦合的2
×
2矩阵模型，如下，
[0081][0082][0083]
其中，和分别为正、负序下频率f
p
和频率f
p-2f1的小信号电压扰动，和分别为对应的电流响应，y
p
(s)为逆变器在频率f
p
下的正序导纳，yn(s-j2ω1)为逆变器在频率f
p-2f1下的负序导纳，yc(s)和yr(s-j2ω1)为逆变器在频率f
p
下正序和频率f
p-2f1下负序之间的耦合导纳。
[0084]yp
(s)为逆变器在频率f
p
下的正序电压扰动下，产生的正序电流响应与正序电压扰动的比值，yc(s)为产生的频率为f
p-2f1的负序电流响应与正序电压扰动的比值。同时，由负频率的负序分量与正频率正序分量的共轭关系可得，
[0085][0086][0087]
其中，“*”表示复数共轭，因此，可在只计算y
p
(s)和yc(s)的基础上，通过上述关系得出yn(s-j2ω1)和yr(s-j2ω1)。
[0088]
根据已有的研究成果，可得考虑直流母线电容动态和电压控制的并网逆变器阻抗/导纳模型为，
[0089][0090][0091]
其中，
[0092][0093][0094]
[0095][0096]y00
(s)＝sc
dc
[0097][0098][0099][0100]
ω1为基波角频率。v1逆变器端口电压相量，i1为逆变器输出电流相量，ps为逆变器输出的有功功率，qs为逆变器输出的无功功率，v1、i1、ps和qs表示了逆变器运行的稳态工作点。l为逆变器的滤波电感，c
dc
为逆变器的直流母线电容，v
dc
为直流母线的额定电压，也是稳态运行点下的直流母线电压，kd为电流控制的解耦系数，一般等于ω1l。
[0101]
由以上式(1-5)和式(1-6)可以看出，阻抗/导纳模型与逆变器的电路参数、控制参数，以及运行工况均密切相关，故阻抗特性重塑的方法涉及控制参数优化和运行工况。可以通过基于多组运行工况对控制参数进行优化，实现阻抗特性重塑，从而改善和抑制新能源并网系统振荡问题。
[0102]
双馈风电机组的电路和控制结构包括：感应异步发电机，背靠背变流器(网侧变流器和机侧变流器)，以及相应的锁相环、电流环等控制回路。va、vb、vc为电网侧(即定子侧)电压；i
ca
、i
cb
、i
cc
为网侧变流器输出三相电流；i
sa
、i
sb
、i
sc
为感应发电机定子侧输出电流；i
ra
、i
rb
、i
rc
为机侧变流器输出电流；θm为转子位置角；lf为网侧变流器交流滤波电感；θ
pll
为锁相环输出的电网电压相位角；h
pll
(s)为锁相环控制器传递函数；i
cd
、i
cq
为网侧变流器dq轴电流指令；i
rd
、i
rq
为机侧变流器dq轴电流指令；h
si
(s)、h
ri
(s)分别为网侧和机侧变流器电流控制器传递函数，k
sd
、k
rd
分别为网侧和机侧变流器电流控制解耦系数；m
sa
、m
sb
、m
sc
为网侧变流器调制信号；m
ra
、m
rb
、m
rc
为机侧变流器调制信号。
[0103]
可以看出，双馈风电机组的感应发电机定子侧和网侧变流器均接入电网，因此，其阻抗相当于机侧阻抗(包括定子、转子及机侧变流器)和网侧变流器阻抗两部分的并联。由于网侧变流器的电路和控制与上述典型并网逆变器相同，这里不再详述其表达式，参考并网逆变器的频率耦合序阻抗模型，将其定义为，
[0104][0105]
双馈风电机组机侧阻抗与感应发电机的电路参数、控制环节、运行工作点均有关。
基本控制环节和网侧逆变器相同，均考虑锁相环和电流环，而机侧阻抗和网侧变流器阻抗的主要不同体现在基本的电路结构，感应电机内部定转子分量之间由定转子电阻、定转子漏电感、互感组成，且受电机转速的影响。因此，机侧阻抗表达式相对较为复杂，同理将其形式定义为，
[0106][0107]
其中，
[0108][0109][0110]
式中，s'＝s-j2ω1；i
r1
为稳定运行工作点的转子侧变流器的输出电流基波相量；m
r1
为稳定运行工作点的转子侧变流器的电压调制信号相量；ke为定、转子匝数比；rr和rs为折算到定子侧的定、转子电阻；lr和ls为折算到定子侧的定、转子电感；lm为折算到定子侧的定、转子间的等效互感；σ
p
(s)＝(s-j2πfm)/s为转差系数。
[0111]
根据式(1-3)和(1-4)的共轭关系，可同理得到y
n,r
(s
′
)和y
r,r
(s
′
)的表达式。最终，可得出双馈风电机组的阻抗模型为，
[0112]ydfig
(s)＝y
dfig,gsc
(s) y
dfig,r
(s)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1-11)
[0113]
由以上式(1-7)至式(1-11)可以看出，阻抗模型与逆变器的电路参数、控制参数，以及运行工况均密切相关，故可以得到类似的结论，即其阻抗特性重塑的方法涉及控制参数优化和运行工况。
[0114]
最后对新能源并网系统宽频带动态建模与稳定性分析：
[0115]
基于建立的新能源发电装置序阻抗模型，可为单机并网系统建模构建等效电路模型。
[0116]
设电网电压注入正序小扰动信号风电机组端口电压扰动分量在耦合频率下产生的负序电流响应由受控电流源表示，该电流通过电网和风电机组阻抗将在机组端口产生负序电压该电压又将在原扰动频率下产生正序电流响应和表达式如下，
[0117][0118][0119]
上述等效电路模型可描述为以正、负序电压扰动为输入，以正、负序响应为输出的两输入两输出系统，其开环传递函数矩阵如式(1-14)所示，系统的稳定性判别应采用推广的nyquist判据，通过分析式(1-7)两个特征根的轨迹对(-1，0j)点的环绕判定系统稳定性。
[0120][0121]
同时，经过整理，等效电路模型和上述两输入两输出系统，可变换为单输入单输出系统，单输入单输出系统的闭环传递函数为：
[0122][0123]
其中，yq(s)表征了逆变器与电网之间的频率耦合效应，可视为逆变器并联了一个额外的导纳/阻抗，
[0124][0125]
于是，系统的稳定性分析仍可以基于阻抗比(y
p
(s) yq(s))/yg(s)采用单输入单输出的nyquist稳定性判据进行分析。
[0126]
另一方面，系统的稳定性可以通过分析闭环传递函数的极点分布实现，由式(1-15)可知，求解系统的闭环极点，
[0127]y1wt
(s)＝y
p
(s) yq(s) y
gp
(s)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1-17)
[0128]
设风电机组接入给定的电网系统存在ns个振荡模式，即式(1-17)存在ns个零点即闭环极点，
[0129]
si＝-σi jωi，i＝1,2,...,ns[0130]
当所有闭环极点均处于复平面的左半平面时，即对任意si，均有σi＞0，则系统稳定。事实上，在基于阻抗的系统模型中，系统的闭环极点等价于系统的振荡/谐振模式。σi为振荡模式的衰减因子，ωi为振荡频率，若σi＞0，即衰减因子为正，则振荡模式i稳定，反之，则不稳定。
[0131]
针对每个振荡模式，可计算阻尼比，
[0132][0133]
基于上述内容，本实施例提出一种基于阻尼比灵敏度确定风电机组的控制参数的方法，包括：
[0134]
s1、基于风电机组发电装置结构确定阻抗模型；
[0135]
s2、基于所述阻抗模型，选取多组运行工况和多组控制参数；
[0136]
s3、基于各工况，依次对各组控制参数在各振荡模式下进行稳定性分析，确定所述工况下的优化参数组；
[0137]
s4、选取满足所有工况稳定性要求的至少一个优化参数组作为风电机组的控制参数；
[0138]
其中，所述阻抗模型至少包括：多个控制参数和多个运行工况；所述在各振荡模式下进行稳定性分析包括基于各工况利用各参数对所有振荡模式下的阻尼比的灵敏度确定所述工况下的优化参数组。
[0139]
在执行本发明的步骤s1之前，先确定控制参数与各振荡模式下阻尼比的灵敏度的关系：
[0140]
设新能源发电装置可优化的控制参数个数为nk，对每个参数取不同的值，组成参数组k，
[0141]
k＝{kj}，j＝1,2,...,nkꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2-1)
[0142]
定义控制参数对振荡模式阻尼比的灵敏度为，
[0143][0144]
若ρ
ij
＞0，则表明增大控制参数kj，振荡模式i的阻尼比ξi增大，该振荡模式趋于稳定，即控制参数kj与振荡模式i的稳定性正相关；反之，则说明控制参数kj与振荡模式i的稳定性负相关。同时，ρ
ij
绝对值越大，则表明控制参数kj对振荡模式i的稳定性作用越大。可通过控制参数对多个振荡模式阻尼比灵敏度的加权求和，以获取该控制参数对系统宽频带振荡稳定性提升最大的优化方向。
[0145]
下面，对本发明提出的一种基于阻尼比灵敏度的风电机组阻抗重塑方法进行详细介绍。
[0146]
步骤s1、基于风电机组发电装置结构确定阻抗模型，具体包括：
[0147]
1，确定电网阻抗(或导纳)，对于感性弱电网，给定短路比以及等效电网阻抗的阻感比等参数，对于串补电网，给定等效电网阻抗的电阻、电抗和串联电容等参数，写出电网阻抗形式如下，
[0148][0149][0150]
其中，rg为电网阻抗的等值电阻，lg为电网阻抗的等值电感，cg为电网阻抗的等值电容，z
gp
(s)和y
gp
(s)为扰动频率下的电网正序阻抗和导纳，z
gn
(s
′
)和y
gn
(s
′
)为耦合频率下的电网负序阻抗和导纳。
[0151]
需要注意的是，对于串补电网，等值阻抗表达式可能与串补电容接入的位置有关，这里仅给出标准形式以说明方法。
[0152]
2，确定风电机组的电路与控制结构，得出其阻抗模型，为描述控制参数和运行工况对阻抗特性的影响，将阻抗模型重写为，
[0153]yp
(s,k,o)，yc(s,k,o)，yr(s
′
,k,o)和yn(s
′
,k,o)
ꢀꢀꢀꢀꢀꢀꢀ
(2-5)
[0154]
其中，o＝{p,q,v,i}为风电机组的运行工况参数组，表明了机组输出有功、无功功率，端口电压和输出电流的基波相量；y
p
(s,k,o)为风电机组在扰动频率下的正序导纳，yc(s,k,o)为风电机组扰动频率到耦合频率的耦合导纳，yr(s
′
,k,o)为风电机组耦合频率到扰动频率的耦合导纳，yn(s
′
,k,o)为风电机组在耦合频率下的负序导纳。
[0155]
步骤s2、基于所述阻抗模型，选取多组运行工况和多组控制参数，具体包括：
[0156]
1，选取no组运行工况，om，m＝1,2,...,no，no一般取10～20，下面针对no组运行工况分别进行多回路参数优化。
[0157]
步骤s3、基于各工况，依次对各组控制参数在各振荡模式下进行稳定性分析，确定所述工况下的优化参数组，具体包括：
[0158]
步骤1，从多组运行工况中依次选取工况；
[0159]
步骤2，基于选取的工况，对各控制参数组依次计算各振荡模式下的阻尼比；
[0160]
步骤3，当所述控制参数组对应的各振荡模式下的阻尼比均满足稳定性要求时，所述参数组为所述工况对应的优化参数组，继续执行步骤1，直到计算完所有工况后退出；否则，执行步骤4；
[0161]
步骤4，依次对各参数组中各参数叠加扰动继续计算各参数对所有振荡模式下的阻尼比的灵敏度；并计算所述参数所有振荡模式下的阻尼比的灵敏度之和；
[0162]
步骤5，基于灵敏度之和对各参数基于加权法进行优化调整，执行步骤2。
[0163]
其过程具体为：
[0164]
将参数组k(首轮优化代入初始参数ko)和工况om代入，根据公式如下，
[0165]y1wt
(s)＝y
p
(s) yq(s) y
gp
(s)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2-6)
[0166]
根据式(2-4)获取系统闭环传递函数y
1wt
(s,k,om)，求取非零根，得出系统的ns个振荡模式，si＝-σi jωi，以及阻尼比ξi，i＝1,2,...,ns，
[0167]y1wt
(s,k,om)＝0
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2-7)
[0168]
对任意振荡模式si，均有阻尼比ξi≥ξ
ε
，ξ
ε
＞0，参数组k能够保证该工况下的稳定，定义为获取工况om下的优化参数组km，回到第三步，继续对下一组运行工况进行参数优化；反之，则获取所有阻尼比ξi＜ξ
ε
的风险振荡模式，设存在n
sr
个风险振荡模式，si＝-σi jωi，i＝1,2,...,n
sr
。
[0169]
依次对所有可优化参数kj叠加扰动δkj，j＝1,2,...,nk，δkj是以初始值为基准的标幺值，
[0170]
0.01≤δkj≤0.05
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2-8)
[0171]
由式(2-4)计算各个参数kj对各个振荡模式阻尼比的灵敏度，ρ
ij
，i＝1,2,...,n
sr
，j＝1,2,...,nk，接着，计算各个参数kj对所有振荡模式阻尼比的灵敏度之和，
[0172][0173]
值得注意的是，控制参数kj对不同振荡模式阻尼比变化的作用方向可能相反，即灵敏度的正负不同，通过阻尼比灵敏度求和，可综合得出该参数对所有振荡模式的整体作用，以确定提升整体宽频带特性的优化方向。同时，对于实际振荡问题中突出的不稳定振荡模式，可采用加权的方式加快对该振荡模式稳定性的提升效果和参数优化效率。
[0174]
对可优化参数kj，j＝1,2,...,nk，进行优化调整，
[0175][0176]
返回第四步进行下一轮迭代优化。
[0177]
最终获取no组运行工况下的优化参数组km，m＝1,2,...,no。
[0178]
上述方法通过综合考虑不同控制回路参数对风电机组并网系统不同频段振荡稳定性的影响及量化指标。在实际应用中，既可根据运行工况，动态调整控制参数，以满足不同工况下的稳定要求；也可选取最恶劣运行工况下的优化参数组，并对不同工况下的稳定性进行校验，选取一组或多组能够满足所有运行工况下稳定性的参数。
[0179]
本实施例基于阻尼比灵敏度确定风电机组控制参数的方法，基于频域阻抗方法，建立风电机组多回路控制参数与其宽频带动态特性的内在联系，定义了多回路控制参数对新能源并网振荡模式的阻尼比灵敏度指标，基于该阻尼比灵敏度指标，提出了多回路参数协同优化的新能源发电阻抗重塑方法，实现宽频带振荡的抑制。
[0180]
新能源发电装置的阻抗特性受电路参数、控制参数和运行工况影响，考虑到其中电路参数在装置设计时已经确定，故仅可通过控制参数的优化的方式重塑阻抗模型，使其在更大的工况范围内均适用。
[0181]
显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
[0182]
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0183]
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0184]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0185]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0186]
以上仅为本发明的实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均包含在申请待批的本发明的权利要求范围之内。