本发明涉及新能源大规模接入系统中无功电压调节技术领域，具体地，涉及一种基于灵敏度的风电场无功调节优先级选择方法。
背景技术：
我国新能源大规模集群式发展，对于电网来说带来了很多挑战。近几年来随着风电并网容量直线上升，随之而来的无功电压运行和控制问题也更加复杂。通过对不同风电机组和各种无功补偿设备的无功电压特性研究可以看出，对于大规模风电接入的电网的无功电压控制需要风电场互相协调配合。
在实现本发明的过程中，发明人发现现有技术中至少存在以下缺陷:首先对于大规模风电并入的系统，在无功电压就地控制时无功补偿设备往往不能满足电压稳定运行的要求；其次除了鼠笼式风电机组以外，其他风电机组均具有良好的无功补偿能力，但无功补偿容量较少；大规模风电场集群并入电网的系统，多风电机组投入无功电压控制时亟需一套协调控制策略以期实现风电场间的互相协调配合。
技术实现要素：
本发明的目的在于，针对上述问题，提出一种基于灵敏度的风电场无功调节优先级选择方法，以期获得某风电场并网点母线电压与其相的风电场输出无功功率的灵敏度系数，以实现在无功电压就地控制时风电场之间的互相协调配合。
为实现上述目的，本发明采用的技术方案是：一种基于灵敏度的风电场无功调节优先级选择方法，主要包括：
步骤1：计算潮流列出节点功率平衡方程；
步骤2：通过对节点功率方程中变量划分为扰动变量，状态变量，控制变量，获得灵敏度矩阵Sxu，Sxp；
步骤3：通过灵敏度矩阵和节点功率平衡方程求得风电场i的电压关于风电场j的输出无功功率变化的灵敏度S1和风电场j的电压关于风电场j的输出无功功率变化的灵敏度S2；通过获得风电场并网点电压与其相邻风电场之间的灵敏度，最终获得大规模风电场无功电压就地调节，互相补偿的风电场动作的优先级顺序。
进一步地，步骤1中，所述节点功率平衡方程的计算步骤包括：
设电力系统中节点i的发电机功率为PGi+jQGi，负荷功率为PLDi+jQLDi，以极坐标表示节点电压时，节点的功率平衡方程为：
ΔPi＝PGi-PLDi-Ui∑Uj(Gijcosδij+Bijsinδij)＝0
ΔQi＝QGi-QLDi-Ui∑Uj(Gijsinδij-Bijcosδij)＝0 (1)
式(1)中，ΔPi,ΔQi为注入节点i的有功功率和无功功率；
PGi,QGi为注入节点i的发电机有功功率和无功功率；
PLDi,QLDi为注入节点i的负荷有功功率和无功功率；
Ui,Uj为节点i的电压；
Gij+jBij为节点i、j之间的导纳；
δij为节点i、j之间的电压相角差，δij＝δi-δj。
进一步地，步骤2中，获得灵敏度矩阵的Sxu，Sxp的计算步骤包括：
若以p、u、x分别表示扰动变量、控制变量和状态变量，电力系统节点功率平衡方程可表示为：
f(x,u,p)＝0 (2)
式(2)中p——扰动变量；
u——控制变量；
x——状态变量；
设包含风电场区域运行在状态(x(0),u(0),p(0))下，如某一风电场有功出力波动Δp，常规发电机做了相应的调整Δu，系统相应发生了状态偏离Δx，则系统的功率平衡方程应为：
f(x(0)+Δx,u(0)+Δu,p(0)+Δp)＝0 (3)
式中x(0),u(0),p(0)——系统运行状态，母线运行电压，风电场有功出力；
Δx——系统状态变化量；
Δu——电压变化量；
Δp——风电场有功出力变化量；
假设稳态情况下各变量的变化很小，以至可忽略其二阶及以上的各项，则式(3)可以展开化简为：
式(4)中均为雅可比矩阵；
将基于状态(x(0),u(0),p(0))的矩阵、矩阵和矩阵分别简记为Jx、Ju和Jp，则有
JxΔx+JuΔu+JpΔp＝0 (5)
式中Jx、Ju、Jp均为雅可比矩阵；
上式(5)称为灵敏度方程，由该式可解出状态变量的偏离值：
式中为灵敏度矩阵；
为灵敏度矩阵。
进一步地，步骤3中，获得风电场并网点电压与其相邻风电场之间的灵敏度计算步骤包括：
设J为系统极坐标形式的全阶雅可比矩阵，ΔY和ΔX分别为系统注入变化量和状态变化量，将ΔY分成两部分，即：ΔY＝[ΔYIΔYN]T，对应的ΔX则变为ΔX＝[ΔXIΔXN]T，相应地将J划分为4部分，即：
式(7)中JII、JIN、JNI、JNN——部分全阶雅克比矩阵，
注入量与状态变化量之间的关系为
式(8)中，
ΔYI——被控风电场接入地区节点的注入变化量；
ΔXI——被控风电场接入地区节点的状态变化量；
ΔYN——其余节点的注入变化量；
ΔXN——其余节点的状态变化量，
设ΔYN＝0，即实际注入量只包含ΔYI，则由式(8)可推出：
由于JNN非奇异从而可得：
对JSI求逆得：
则被控风电场接入点地区节点之间的关系为
式中 ΔUw——风电场出口节点的电压；
Δθw——风电场出口节点的相角变化向量；
ΔUc——接入点的电压；
Δθc——接入点的相角变化向量；
ΔPW——风电场出口节点注入有功变化向量；
ΔQW——风电场出口节点注入无功变化向量；
ΔPC——接入点注入有功变化向量；
ΔQC——接入点注入无功变化向量，
无功功率与电压的灵敏度关系，则由式(12)和式(13)可得风电场i的电压关于风电场j的输出无功功率变化的灵敏度关系为
S1＝A2m+n+i,m+n+j (14)
式(14)中S1——风电场i的电压关于风电场j的输出无功功率变化的灵敏度。
风电场j的电压关于风电场j的输出无功功率变化的灵敏度关系为
S2＝Am+j,m+n+j (15)
式(15)中S2——风电场j的电压关于风电场j的输出无功功率变化的灵敏度。
本发明各实施例的一种基于灵敏度的风电场无功调节优先级选择方法，由于主要包括：通过无功-电压灵敏度矩阵和节点功率平衡方程求得风电场并网点电压与其相邻风电场之间的灵敏度，即可获得大规模风电场无功电压就地调节，互相补偿的风电场动作的优先级顺序；从而可以克服现有技术中风电场群无功就地控制无序、滞后等问题。
本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。
下面通过实施例，对本发明的技术方案做进一步的详细描述。
具体实施方式
以下结合对本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。
具体地，基于灵敏度的风电场优先级选择法是由于不同执行站发出的无功对并网点电压响应特性不同。因此，为了让执行站能够尽可能准确快速地参与系统电压控制的核心就是能够估算出风电场输出无功功率的变化引起的相关节点电压变化，即两者之间的灵敏度关系。计算出灵敏度系数后，选取对并网点母线电压灵敏度最高的执行站。主要包括：
第一步：节点功率平衡方程；设电力系统中节点i的发电机功率为PGi+jQGi，负荷功率为PLDi+jQLDi，以极坐标表示节点电压时，节点的功率平衡方程为
ΔPi＝PGi-PLDi-Ui∑Uj(Gijcosδij+Bijsinδij)＝0
ΔQi＝QGi-QLDi-Ui∑Uj(Gijsinδij-Bijcosδij)＝0 (1)
式中：
ΔPi,ΔQi——注入节点i的有功功率和无功功率；
PGi,QGi——注入节点i的发电机有功功率和无功功率；
PLDi,QLDi——注入节点i的负荷有功功率和无功功率；
Ui,Uj——节点i的电压；
Gij+jBij——节点i、j之间的导纳；
δij——节点i、j之间的电压相角差，δij＝δi-δj。
第二步：灵敏度方程，获得灵敏度矩阵；若以p、u、x分别表示扰动变量、控制变量和状态变量。电力系统节点功率平衡方程可表示为：
f(x,u,p)＝0 (2)
式中 p——扰动变量；
u——控制变量；
x——状态变量。
设包含风电场区域运行在状态(x(0),u(0),p(0))下，如某一风电场有功出力波动Δp，常规发电机做了相应的调整Δu，系统相应发生了状态偏离Δx，则系统的功率平衡方程应为：
f(x(0)+Δx,u(0)+Δu,p(0)+Δp)＝0 (3)
式中x(0),u(0),p(0)——系统运行状态，母线运行电压，风电场有功出力；
Δx——系统状态变化量；
Δu——电压变化量；
Δp——风电场有功出力变化量。
假设稳态情况下各变量的变化很小，以至可忽略其二阶及以上的各项，则式3-10可以展开化简为：
式中——雅可比矩阵。
如将基于状态(x(0),u(0),p(0))的矩阵分别简记为Jx、Ju、Jp，则有
JxΔx+JuΔu+JpΔp＝0 (5)
式中Jx、Ju、Jp——雅可比矩阵。
上式称为灵敏度方程，由该式可解出状态变量的偏离值：
式中——灵敏度矩阵；
——灵敏度矩阵。
第三步：并网点电压与无功源输出无功功率的灵敏度
对于任意给定变化量的控制变量或扰动变量，采用上式均可求得各状态变量的变化量。
将上述方法应用于电力系统，求得无功-电压灵敏度矩阵。本文只关注部分节点，将电力系统中的无功补偿安装点，负荷节点和发电机节点作为注入节点，将全阶雅克比矩阵进行降阶，具体的求取方法如下：
设J为系统极坐标形式的全阶雅可比矩阵，ΔY和ΔX分别为系统注入变化量和状态变化量，将ΔY分成两部分，即：ΔY＝[ΔYIΔYN]T，对应的ΔX则变为ΔX＝[ΔXIΔXN]T。相应地将J划分为4部分，即：
式中JII、JIN、JNI、JNN——部分全阶雅克比矩阵。
注入量与状态变化量之间的关系为
式中 ΔYI——被控风电场接入地区节点的注入变化量；
ΔXI——被控风电场接入地区节点的状态变化量；
ΔYN——其余节点的注入变化量；
ΔXN——其余节点的状态变化量。
设ΔYN＝0，即实际注入量只包含ΔYI，则由式3-15可推出：
由于JNN非奇异从而可得：
对JSI求逆得：
则被控风电场接入点地区节点之间的关系为
式中 ΔUw——风电场出口节点的电压；
Δθw——风电场出口节点的相角变化向量；
ΔUc——接入点的电压；
Δθc——接入点的相角变化向量；
ΔPW——风电场出口节点注入有功变化向量；
ΔQW——风电场出口节点注入无功变化向量；
ΔPC——接入点注入有功变化向量；
ΔQC——接入点注入无功变化向量。
本文只研究无功功率与电压的灵敏度关系，则由式(12)和式(13)可得风电场i的电压关于风电场j的输出无功功率变化的灵敏度关系为
S1＝A2m+n+i,m+n+j (14)
式中S1——风电场i的电压关于风电场j的输出无功功率变化的灵敏度。
风电场j的电压关于风电场j的输出无功功率变化的灵敏度关系为
S2＝Am+j,m+n+j (15)
式中S2——风电场j的电压关于风电场j的输出无功功率变化的灵敏度。
至少可以达到以下有益效果：首先当某风电场无功补偿设备用尽时，还可以利用风电机组进行无功补偿，作为无功电压控制的最后一道防线；其次通过基于灵敏度的风电场无功调节优先级选择方法，对大规模风电场群提出了一套风电场之间无功电压协调控制方法，以期达到风电场无功电压就地控制的高效性和快速性。
最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。