本发明涉及一种适用于电网中风电系统的动态无功设备无功功率优化方法，属于电力系统安全和控制技术领域。
背景技术：
我国张北地区千万千瓦级大规模风电馈入弱送端电力系统，出力随机波动强，运行特性非常复杂。而部分风电汇集区域动态无功储备不够，因此电压安全问题突出，在严重故障下的安全问题尤为显著。而动态无功储备能够在电网发生扰动后提供快速无功电压支撑，保障电网安全。由此可见，如何在线优化大规模风电系统的动态无功储备，保障电网安全运行，是十分必要的。
技术实现要素：
本发明的目的是提出一种适用于电网中风电系统的动态无功设备无功功率优化方法。基于已有技术的背景，解决针对大规模风电系统的动态无功储备在线优化问题，以确保电网发生扰动后，仍能维持稳定运行，显著提高用电可靠性。
本发明提出的适用于电网中风电系统的动态无功设备无功功率优化方法，包括以下步骤：
(1)计算电网中风电系统的各风电场内所有风电机组的无功功率调节范围和所有动态无功设备的无功功率的调节范围，过程如下：
分别建立计算电网中风电系统侧第k个风电场的所有风电机组的无功功率调节量上限值
其中，ng,k与nd,k分别为风电系统侧第k个风电场中风电机组和动态无功设备的数量，δqg,k,i为第k个风电场中第i个风电机组的无功功率变化量，δqd,k,j表示第k个风电场中第j个动态无功设备的无功功率变化量；
上述目标函数的约束条件如下：
当前状态下风电场母线电压的约束条件为：
其中，sg,k,m和sd,k,n分别表示第k个风电场中第m个风电机组和第n个动态无功设备的无功功率对该风电场母线电压的灵敏度系数，vk表示当前状态下第k个风电场母线电压值，vk和
当前状态下风电机组机端电压的约束条件为：
其中，cg,k,m,i和cd,k,n,i分别表示第k个风电场中第m个风电机组和第n个动态无功设备的无功功率对第i个风电机组的电压的灵敏度系数，vg,k,i表示第k个风电场中第i个风电机组的机端电压值，vg,k,i和
第k个风电场中第i个风电机组无功功率出力约束条件为：
其中，δqg,k,i和
第k个风电场中第j个动态无功设备无功功率出力约束条件为：
其中，δqd,k,j和
(2)根据电网中风电系统中风电汇集站记录的当前各风电场的运行状态，生成n个模拟场景，n为风电场汇集站区域内风电场的个数，并计算n个模拟场景下的安全裕度，包括以下步骤：
(2-1)根据当前各风电场的运行状态，风电汇集站生成n个模拟场景，利用下式，计算其中任意第t个模拟场景的无功功率，t＝1,2,...,n：
其中，
(2-2)利用牛顿-拉夫逊潮流计算方法，根据上述步骤(2-1)的n个模拟场景，计算每个风电场的潮流分布，得到第t个模拟场景中第k个风电场的母线电压值
(2-3)重复步骤(2-2)，分别得到与各模拟场景相对应的风电场汇集区域的安全裕度；
(3)根据上述与各模拟场景相对应的风电场汇集区域的安全裕度，拟合得到风电场的所有动态无功设备的无功功率调节量总和与风电场汇集区域安全裕度的关系如下：
其中，r表示风电场汇集的区域安全裕度，ak和bk为拟合得到的回归系数，ak和bk通过下式计算得到：
(4)根据上述安全裕度的关系，计算当前状态下风电场中各动态无功设备的无功功率设定值，包括以下步骤：
(4-1)分别建立各风电场中所有动态无功设备的最优无功功率调节量总和、所有风电机组最优无功功率调节量总和的目标函数如下：
上述目标函数的约束条件如下：
当前状态下各风电场的母线电压的约束条件为：
其中，hg,k,l表示风电场汇集站区域内第k个风电场的风电机组无功功率对第l个风电场母线电压的灵敏度系数，hd,k,l表示风电场汇集站区域内第k个风电场的动态无功设备无功功率对第l个风电场母线电压的灵敏度系数，vl表示当前状态下第l个风电场的母线电压值，vl和
第k个风电场的所有动态无功设备的无功功率总和的约束条件为：
第k个风电场的所有风电机组的无功功率总和的约束条件为：
根据约束条件，求解上述目标函数，得到风电场中所有动态无功设备的最优无功功率调节量总和
(4-2)分别建立计算各风电场内每个动态无功设备的无功功率设定值的目标函数如下：
上述目标函数的约束条件如下：
当前状态下风电场母线电压的约束条件为：
其中，sg,k,m和sd,k,n分别表示第k个风电场中第m个风电机组和第n个动态无功设备的无功功率对该风电场母线电压的灵敏度系数，vk表示当前状态下第k个风电场母线电压值，vk和
当前状态下风电机组机端电压的约束条件为：
其中，cg,k,m,i和cd,k,n,i分别表示第k个风电场中第m个风电机组和第n个动态无功设备的无功功率对第i个风电机组的电压的灵敏度系数，vg,k,i表示第k个风电场中第i个风电机组的机端电压值，vg,k,i和
第k个风电场中第i个风电机组无功功率出力约束条件为：
其中，δqg,k,i和
第k个风电场中第j个动态无功设备无功功率出力约束条件为：
其中，δqd,k,j和
第k个风电场所有风电机组的无功功率值总和约束条件为：
第k个风电场所有动态无功设备的无功功率值总和约束条件为：
根据上述步骤(4-1)的计算结果和上述约束条件，求解上述目标函数，分别得到风电场中各动态无功设备的最优无功功率调节量的
本发明提出的适用于大规模风电系统的动态无功储备优化方法，其优点是：
本发明方法解决了大规模风电系统的动态无功储备在线优化问题，确保电网发生扰动后，仍能维持稳定运行，因此显著提高用电可靠性。而且本发明方法可以在线优化大规模风电系统的动态无功储备配置，计算速度快，可以很大程度地提高电网运行的稳定性与可靠性。
具体实施方式
本发明提出的适用于电网中风电系统的动态无功设备无功功率优化方法，包括以下步骤：
(1)计算电网中风电系统的各风电场内所有风电机组的无功功率调节范围和所有动态无功设备(例如静态无功功率补偿器、静态无功功率发生器等等)的无功功率的调节范围，过程如下：
分别建立计算电网中风电系统侧第k个风电场的所有风电机组的无功功率调节量上限值
其中，ng,k与nd,k分别为风电系统侧第k个风电场中风电机组和动态无功设备的数量，δqg,k,i为第k个风电场中第i个风电机组的无功功率变化量，δqd,k,j表示第k个风电场中第j个动态无功设备的无功功率变化量；
上述目标函数的约束条件如下：
当前状态下风电场母线电压的约束条件为：
其中，sg,k,m和sd,k,n分别表示第k个风电场中第m个风电机组和第n个动态无功设备的无功功率对该风电场母线电压的灵敏度系数，vk表示当前状态下第k个风电场母线电压值，vk和
当前状态下风电机组机端电压的约束条件为：
其中，cg,k,m,i和cd,k,n,i分别表示第k个风电场中第m个风电机组和第n个动态无功设备的无功功率对第i个风电机组的电压的灵敏度系数，vg,k,i表示第k个风电场中第i个风电机组的机端电压值，vg,k,i和
第k个风电场中第i个风电机组无功功率出力约束条件为：
其中，δqg,k,i和
第k个风电场中第j个动态无功设备无功功率出力约束条件为：
其中，δqd,k,j和
(2)根据电网中风电系统中风电汇集站记录的当前各风电场的运行状态，生成n个模拟场景，n为风电场汇集站区域内风电场的个数，并计算n个模拟场景下的安全裕度，包括以下步骤：
(2-1)根据当前各风电场的运行状态，风电汇集站生成n个模拟场景，利用下式，计算其中任意第t个模拟场景的无功功率，t＝1,2,...,n：
其中，
(2-2)利用牛顿-拉夫逊潮流计算方法，根据上述步骤(2-1)的n个模拟场景，计算每个风电场的潮流分布，得到第t个模拟场景中第k个风电场的母线电压值
(2-3)重复步骤(2-2)，分别得到与各模拟场景相对应的风电场汇集区域的安全裕度；
(3)根据上述与各模拟场景相对应的风电场汇集区域的安全裕度，拟合得到风电场的所有动态无功设备的无功功率调节量总和与风电场汇集区域安全裕度的关系如下：
其中，r表示风电场汇集的区域安全裕度，ak和bk为拟合得到的回归系数，ak和bk通过下式计算得到：
(4)根据上述安全裕度的关系，计算当前状态下风电场中各动态无功设备的无功功率设定值，包括以下步骤：
(4-1)分别建立各风电场中所有动态无功设备的最优无功功率调节量总和、所有风电机组最优无功功率调节量总和的目标函数如下：
上述目标函数的约束条件如下：
当前状态下各风电场的母线电压的约束条件为：
其中，hg,k,l表示风电场汇集站区域内第k个风电场的风电机组无功功率对第l个风电场母线电压的灵敏度系数，hd,k,l表示风电场汇集站区域内第k个风电场的动态无功设备无功功率对第l个风电场母线电压的灵敏度系数，vl表示当前状态下第l个风电场的母线电压值，vl和
第k个风电场的所有动态无功设备的无功功率总和的约束条件为：
第k个风电场的所有风电机组的无功功率总和的约束条件为：
根据约束条件，求解上述目标函数，得到风电场中所有动态无功设备的最优无功功率调节量总和
(4-2)分别建立计算各风电场内每个动态无功设备的无功功率设定值的目标函数如下：
上述目标函数的约束条件如下：
当前状态下风电场母线电压的约束条件为：
其中，sg,k,m和sd,k,n分别表示第k个风电场中第m个风电机组和第n个动态无功设备的无功功率对该风电场母线电压的灵敏度系数，vk表示当前状态下第k个风电场母线电压值，vk和
当前状态下风电机组机端电压的约束条件为：
其中，cg,k,m,i和cd,k,n,i分别表示第k个风电场中第m个风电机组和第n个动态无功设备的无功功率对第i个风电机组的电压的灵敏度系数，vg,k,i表示第k个风电场中第i个风电机组的机端电压值，vg,k,i和
第k个风电场中第i个风电机组无功功率出力约束条件为：
其中，δqg,k,i和
第k个风电场中第j个动态无功设备无功功率出力约束条件为：
其中，δqd,k,j和
第k个风电场所有风电机组的无功功率值总和约束条件为：
第k个风电场所有动态无功设备的无功功率值总和约束条件为：
根据上述步骤(4-1)的计算结果和上述约束条件，求解上述目标函数，分别得到风电场中各动态无功设备的最优无功功率调节量的