CN103825281B - 基于动态无功设备的电力系统中中枢母线电压的控制方法 - Google Patents

Abstract

本发明涉及一种基于动态无功设备的电力系统中中枢母线电压的控制方法，属于电力系统评估和控制技术领域。首先，监测系统当中中枢节点电压值，根据该电压值计算出参与电压控制的发电机与动态无功补偿设备的控制指令。将指令下发给相应的发电机与动态无功补偿设备。判断是否需要进行动态无功优化，如果不需要，则在时刻T₀+△T，返回第一步骤，若需要，则根据发电机与动态无功补偿装置的无功功率及上限和下限值，计算动态无功优化的控制指令并下发。发电机与动态无功补偿装置根据指令调整各自的无功功率。本发明方法能够有效的进行暂态电压稳定的预防控制，提高电力系统中动态无功储备水平和扰动后的系统电压的恢复效果。

Description

基于动态无功设备的电力系统中中枢母线电压的控制方法

技术领域

本发明涉及一种基于动态无功设备的电力系统中中枢母线电压的控制方法，属于电力系统评估和控制技术领域。

背景技术

近些年来，电力系统的电压稳定问题时有发生，并且有的已经造成了相当巨大的影响。电压控制方法显得尤为重要。但是由于传统的控制手段有限，导致控制响应速度很慢，本发明中，将动态无功补偿装置作为一种新的控制手段引入进来，提高了电压控制的响应速度。另外，在本发明中，还采用了动态无功优化的方法保证动态无功补偿设备保持在一个有足够的无功功率储备，以防在电力系统发生故障时，没有足够的无功功率容量响应控制的需求。

发明内容

本发明的目的是提出一种基于动态无功设备的电力系统中中枢母线电压的控制方法，该方法将动态无功补偿装置作为一种控制手段，纳入电力系统中枢母线电压控制中来，并且通过动态无功储备优化，实现了对电网发生大故障时的预防性控制。

本发明提出的基于动态无功设备的电力系统中中枢母线电压的控制方法，包括以下步骤：

（1）从电力系统的调度中心获取中枢母线电压的设定值当前采集时刻T0的电力系统中中枢母线的电压值V_p，通过下式计算得到中枢母线电压控制阶段的参与控制的发电机的电压调整量ΔV_g和参与控制的动态无功补偿装置的电压调整量ΔV_s：

$\begin{array}{c}\underset{\mathrm{\Δ}{V}_g,\mathrm{\Δ}{V}_s}{\min }\{W_p{\left|\right|a\·(V_p-V_p^{\mathrm{ref}})+C_{\mathrm{pg}}\mathrm{\Δ}{V}_g+C_{\mathrm{ps}}\mathrm{\Δ}{V}_s\left|\right|}^2+W_q({\left|\right|{\mathrm{\Θ}}_g\left|\right|}^2+{\left|\right|{\mathrm{\Θ}}_s\left|\right|}^2)\}\\ s.t.V_p^{\min }\≤V_p+C_{\mathrm{pg}}\mathrm{\Δ}{V}_g+C_{\mathrm{ps}}\mathrm{\Δ}{V}_s\≤V_p^{\max }\\ Q_g^{\min }\≤Q_g+C_g\mathrm{\Δ}{V}_g+C_{\mathrm{gs}}\mathrm{\Δ}{V}_s\≤Q_g^{\max }\end{array}$

$\begin{array}{c}{Q}_s^{\min }\≤Q_s+C_{\mathrm{sg}}\mathrm{\Δ}{V}_g+C_s\mathrm{\Δ}{V}_s\≤Q_s^{\max }\\ \left|\mathrm{\Δ}{V}_g\right|\≤\mathrm{\Δ}{V}_g^{\max }\\ \left|\mathrm{\Δ}{V}_s\right|\≤\mathrm{\Δ}{V}_s^{\max }\end{array}$

其中，Θ_g为发电机无功功率的裕度矢量，该裕度矢量Θ_g中的第i个分量为：

${\mathrm{\Θ}}_{g_i}=\frac{Q_{g_i}+\underset{j\∈\mathrm{Ng}}{\mathrm{\Σ}}{C}_{g_{\mathrm{ij}}}\mathrm{\Δ}{V}_{g_j}+\underset{k\∈\mathrm{Ns}}{\mathrm{\Σ}}{C}_{g{s}_{\mathrm{ik}}}\mathrm{\Δ}{V}_{s_k}-Q_{g_i}^{\min }}{Q_{g_i}^{\max }-Q_{g_i}^{\min }}$

Θ_s为动态无功补偿设备无功功率的裕度矢量，裕度矢量Θ_s的第i个分量为：

${\mathrm{\Θ}}_{s_i}=\frac{Q_{s_i}+\underset{j\∈\mathrm{Ns}}{\mathrm{\Σ}}{C}_{s_{\mathrm{ij}}}\mathrm{\Δ}{V}_{s_j}+\underset{k\∈\mathrm{Ng}}{\mathrm{\Σ}}{C}_{{\mathrm{sg}}_{\mathrm{ik}}}\mathrm{\Δ}{V}_{g_k}-Q_{s_i}^{\min }}{Q_{s_i}^{\max }-Q_{s_i}^{\min }}$

其中下标p表示中枢母线，下标g表示发电机，下标s表示动态无功补偿装置；C_pg是从调度中心获取的中枢母线电压偏差对参与控制的发电机的电压调整量ΔV_g的灵敏度矩阵；C_ps是从调度中心获取的中枢母线电压偏差对参与控制的动态无功补偿装置的电压调整量ΔV_s的灵敏度矩阵；W_p为中枢母线电压偏差的权重系数，W_p的取值范围为0.1-1，W_q为发电机和动态无功补偿装置的无功功率储备裕度的权重系数，W_q的取值范围为10-100，α是中枢母线电压的增益系数，α的取值范围为1-10，Q_g是参与控制的发电机当前发出的无功功率，分别为参与控制的发电机的无功功率下限与无功功率上限，Q_s是参与控制的动态无功补偿装置当前发出的无功功率，分别为参与控制的动态无功补偿装置发出的无功功率下限与无功功率的上限，C_g为从调度中心获取的参与控制的发电机的无功功率调整量对参与控制的发电机的电压调整量ΔV_g的灵敏度矩阵，C_s为从调度中心获取的参与控制的动态无功补偿装置的无功功率调整量对参与控制的动态无功补偿装置的电压调整量ΔV_s的灵敏度矩阵，C_gs为从调度中心获取的参与控制的发电机的无功功率调整量对参与控制的动态无功补偿设备的电压调整量ΔV_s的灵敏度矩阵；C_sg为从调度中心获取的参与控制的动态无功补偿装置的无功功率调整量对参与控制的发电机的电压调整量ΔV_g的灵敏度矩阵，为设定的中枢母线电压的下限，为设定的中枢母线电压上限，为从调度中心获取的参与控制的发电机允许的单步最大调整量，为从调度中心获取的参与控制的动态无功补偿装置允许的单步最大调整量，Ns为参与控制的动态无功补偿装置的编号的集合，Ng为参与控制的发电机的编号的集合；

将上述计算得到的参与控制的发电机的电压调整量ΔV_g和参与控制的动态无功补偿装置的电压调整量ΔV_s作为控制指令下达给相应的发电机和动态无功补偿装置，各发电机和动态无功补偿装置根据控制指令，完成中枢母线电压阶段控制；

（2）设定动态无功补偿装置无功功率的裕度阈值为将当前动态无功补偿装置的裕度值与设定的裕度阈值比较，若当前动态无功补偿装置的裕度值大于或等于设定的裕度阈值，则经过中枢母线电压两次控制之间的时间间隔ΔT后重复步骤（1），ΔT的取该数值为5至10分钟，若当前动态无功补偿装置的裕度值小于设定的裕度阈值，则进行步骤（3）；

（3）根据下式计算，得到动态无功补偿装置无功功率裕度优化阶段的参与控制的发电机和动态无功补偿装置的电压调整量ΔV_g和ΔV_s：

$\begin{array}{c}\underset{\mathrm{\Δ}{V}_g,\mathrm{\Δ}{V}_s}{\min }\{W_s{\left|\right|{\mathrm{\μ}}_s\left|\right|}^2+W_g{\left|\right|{\mathrm{\μ}}_g\left|\right|}^2\}\\ s.t.\\ V_p-V_{\mathrm{band}}\≤V_p+C_{\mathrm{pg}}\mathrm{\Δ}{V}_g+C_{\mathrm{ps}}\mathrm{\Δ}{V}_s\≤V_p+V_{\mathrm{band}}\\ Q_g^{\min }\≤Q_g+C_g\mathrm{\Δ}{V}_g+C_{\mathrm{gs}}\mathrm{\Δ}{V}_s\≤Q_s^{\max }\\ Q_s^{\min }\≤Q_s+C_{\mathrm{sg}}\mathrm{\Δ}{V}_g+C_s\mathrm{\Δ}{V}_s\≤Q_s^{\max }\\ \left|\mathrm{\Δ}{V}_g\right|\≤\mathrm{\Δ}{V}_g^{\max }\\ \left|\mathrm{\Δ}{V}_s\right|\≤\mathrm{\Δ}{V}_s^{\max }\end{array}$

上述目标函数中，μ_s为动态无功补偿装置的无功功率矢量，表示动态无功补偿装置经控制后无功功率偏离设定值的大小，μ_g为发电机的无功功率矢量，表示发电机经控制后的无功功率偏离所有发电机无功功率平均值的大小；

${\mathrm{\μ}}_{\mathrm{si}}=\frac{Q_{\mathrm{si}}+\underset{j\∈\mathrm{Ns}}{\mathrm{\Σ}}{C}_{s_{\mathrm{ij}}}\mathrm{\Δ}{V}_{s_j}+\underset{k\∈\mathrm{Ng}}{\mathrm{\Σ}}{C}_{s{g}_{\mathrm{ik}}}\mathrm{\Δ}{V}_{g_k}-Q_{s_i}^{\mathrm{reference}}}{Q_{\mathrm{si}}^{\max }-Q_{\mathrm{si}}^{\min }}$

${\mathrm{\μ}}_{\mathrm{gi}}=\frac{Q_{\mathrm{gi}}+\underset{j\∈\mathrm{Ng}}{\mathrm{\Σ}}{C}_{g_{\mathrm{ij}}}{\mathrm{\ΔV}}_{g_j}+\underset{k\∈\mathrm{Ns}}{\mathrm{\Σ}}{C}_{g{s}_{\mathrm{ik}}}\mathrm{\Δ}{V}_{s_k}-Q_g^{\mathrm{aver}}}{Q_g^{\mathrm{aum}}}$

为动态无功补偿装置的裕度设定值，为当前时刻所有发电机的无功功率平均值，为当前时刻发电机无功功率总和，W_s为动态无功补偿装置的无功功率矢量的权重系数，W_s的取值范围为0.1-1，W_g为发电机的无功功率矢量的权重系数，W_g的取值范围为10-100；

将计算得到的无功功率裕度优化阶段的电压调整量ΔV_g和ΔV_s作为控制指令下达给参与控制的发电机与动态无功储备装置，完成动态无功功率裕度优化阶段控制；

（4）在时刻T₀+△T，返回步骤（1）。

本发明提出的基于动态无功设备的电力系统中中枢母线电压的控制方法，其优点是：

1、本发明提出的基于动态无功设备的电力系统中中枢母线电压的控制方法，将动态无功补偿装置纳入传统的电压控制模型中，充分利用动态无功补偿装置的快速调节特性。

2、本发明方法中建立了动态无功储备优化模型，利用慢速动态无功补偿装置置换出动态无功补偿装置的无功功率，提供电力系统中动态无功储备水平。

3、本发明方法能够有效的进行暂态电压稳定的预防控制，提高扰动后的系统电压的恢复效果。

附图说明

图1是本发明方法的流程框图。

具体实施方式

本发明提出的基于动态无功设备的电力系统中中枢母线电压的控制方法，其流程框图如图1所示，包括以下步骤：

（1）从电力系统的调度中心获取中枢母线电压的设定值当前采集时刻T₀的电力系统中中枢母线的电压值V_p，通过下式计算得到中枢母线电压控制阶段的参与控制的发电机的电压调整量ΔV_g和参与控制的动态无功补偿装置的电压调整量ΔV_s：

$\begin{array}{c}\underset{\mathrm{\Δ}{V}_g,\mathrm{\Δ}{V}_s}{\min }\{W_p{\left|\right|a\·(V_p-V_p^{\mathrm{ref}})+C_{\mathrm{pg}}\mathrm{\Δ}{V}_g+C_{\mathrm{ps}}\mathrm{\Δ}{V}_s\left|\right|}^2+W_q({\left|\right|{\mathrm{\Θ}}_g\left|\right|}^2+{\left|\right|{\mathrm{\Θ}}_s\left|\right|}^2)\}\\ s.t.V_p^{\min }\≤V_p+C_{\mathrm{pg}}\mathrm{\Δ}{V}_g+C_{\mathrm{ps}}\mathrm{\Δ}{V}_s\≤V_p^{\max }\\ Q_g^{\min }\≤Q_g+C_g\mathrm{\Δ}{V}_g+C_{\mathrm{gs}}\mathrm{\Δ}{V}_s\≤Q_g^{\max }\\ Q_s^{\min }\≤Q_s+C_{\mathrm{sg}}\mathrm{\Δ}{V}_g+C_s\mathrm{\Δ}{V}_s\≤Q_s^{\max }\\ \left|\mathrm{\Δ}{V}_g\right|\≤\mathrm{\Δ}{V}_g^{\max }\\ \left|\mathrm{\Δ}{V}_s\right|\mathrm{\Δ}{V}_s^{\max }\end{array}$

其中，Θ_g为发电机无功功率的裕度矢量，该裕度矢量Θ_g中的第i个分量为：

${\mathrm{\Θ}}_{s_i}=\frac{Q_{s_i}+\underset{j\∈\mathrm{Ns}}{\mathrm{\Σ}}{C}_{s_{\mathrm{ij}}}\mathrm{\Δ}{V}_{s_j}+\underset{k\∈\mathrm{Ng}}{\mathrm{\Σ}}{C}_{{\mathrm{sg}}_{\mathrm{ik}}}\mathrm{\Δ}{V}_{g_k}-Q_{s_i}^{\min }}{Q_{s_i}^{\max }-Q_{s_i}^{\min }}$

Θ_s为动态无功补偿设备无功功率的裕度矢量，裕度矢量Θ_s的第i个分量为：

${\mathrm{\Θ}}_{s_i}=\frac{Q_{s_i}+\underset{j\∈\mathrm{Ns}}{\mathrm{\Σ}}{C}_{s_{\mathrm{ij}}}\mathrm{\Δ}{V}_{s_j}+\underset{k\∈\mathrm{Ng}}{\mathrm{\Σ}}{C}_{{\mathrm{sg}}_{\mathrm{ik}}}\mathrm{\Δ}{V}_{g_k}-Q_{s_i}^{\min }}{Q_{s_i}^{\max }-Q_{s_i}^{\min }}$

其中下标p表示中枢母线，下标g表示发电机，下标s表示动态无功补偿装置；C_pg是从调度中心获取的中枢母线电压偏差对参与控制的发电机的电压调整量ΔV_g的灵敏度矩阵；C_ps是从调度中心获取的中枢母线电压偏差对参与控制的动态无功补偿装置的电压调整量ΔV_s的灵敏度矩阵；W_p为中枢母线电压偏差的权重系数，W_p的取值范围为0.1-1，W_q为发电机和动态无功补偿装置的无功功率储备裕度的权重系数，W_q的取值范围为10-100，α是中枢母线电压的增益系数，α的取值范围为1-10，Q_g是参与控制的发电机当前发出的无功功率，分别为参与控制的发电机的无功功率下限与无功功率上限，Q_s是参与控制的动态无功补偿装置当前发出的无功功率，分别为参与控制的动态无功补偿装置发出的无功功率下限与无功功率的上限，C_g为从调度中心获取的参与控制的发电机的无功功率调整量对参与控制的发电机的电压调整量ΔV_g的灵敏度矩阵，C_s为从调度中心获取的参与控制的动态无功补偿装置的无功功率调整量对参与控制的动态无功补偿装置的电压调整量ΔV_s的灵敏度矩阵，C_gs为从调度中心获取的参与控制的发电机的无功功率调整量对参与控制的动态无功补偿设备的电压调整量ΔV_s的灵敏度矩阵；C_sg为从调度中心获取的参与控制的动态无功补偿装置的无功功率调整量对参与控制的发电机的电压调整量ΔV_g的灵敏度矩阵，为设定的中枢母线电压的下限，为设定的中枢母线电压上限，为从调度中心获取的参与控制的发电机允许的单步最大调整量，为从调度中心获取的参与控制的动态无功补偿装置允许的单步最大调整量，Ns为参与控制的动态无功补偿装置的编号的集合，Ng为参与控制的发电机的编号的集合；

将上述计算得到的参与控制的发电机的电压调整量ΔV_g和参与控制的动态无功补偿装置的电压调整量ΔV_s作为控制指令下达给相应的发电机和动态无功补偿装置，各发电机和动态无功补偿装置根据控制指令，完成中枢母线电压阶段控制；

（2）设定动态无功补偿装置无功功率的裕度阈值为对于裕度阈值的设定，一般情况下设定为动态无功补偿装置无功功率的中间值，但是也可以根据电网的实际情况略作调整，比如在可能大量缺少无功功率的地点，即可能发生严重故障的节点附近，设定的最优动态无功储备值可以取动态无功装置出力中间偏小的位置，使得发生严重故障的时候此台动态无功储备装置能够补偿更多的无功功率。将当前动态无功补偿装置的裕度值与设定的裕度阈值比较，若当前动态无功补偿装置的裕度值大于或等于设定的裕度阈值，则经过中枢母线电压两次控制之间的时间间隔ΔT后重复步骤（1），ΔT的取该数值为5至10分钟，若当前动态无功补偿装置的裕度值小于设定的裕度阈值，则进行步骤（3）；

（3）动态无功功率裕度优化控制是在保证中枢母线的电压合格、优质的前提下，将发电机与动态无功补偿装置的无功功率进行置换，使得当系统发生故障需要快速的调节手段时，动态无功补偿装置有足够的无功功率裕度来响应系统的控制需求，提高扰动后的系统电压恢复效果，实现暂态电压稳定的预防控制。

根据下式计算，得到动态无功补偿装置无功功率裕度优化阶段的参与控制的发电机和动态无功补偿装置的电压调整量ΔV_g和ΔV_s：

$\begin{array}{c}\underset{\mathrm{\Δ}{V}_g,\mathrm{\Δ}{V}_s}{\min }\{W_s{\left|\right|{\mathrm{\μ}}_s\left|\right|}^2+W_g{\left|\right|{\mathrm{\μ}}_g\left|\right|}^2\}\\ s.t.\\ V_p-V_{\mathrm{band}}\≤V_p+C_{\mathrm{pg}}\mathrm{\Δ}{V}_g+C_{\mathrm{ps}}\mathrm{\Δ}{V}_s\≤V_p+V_{\mathrm{band}}\end{array}$

$\begin{array}{c}{Q}_g^{\min }\≤Q_g+C_g\mathrm{\Δ}{V}_g+C_{\mathrm{gs}}\mathrm{\Δ}{V}_s\≤Q_g^{\max }\\ Q_s^{\min }\≤Q_s+C_{\mathrm{sg}}\mathrm{\Δ}{V}_g+C_s\mathrm{\Δ}{V}_s\≤Q_s^{\max }\\ \left|\mathrm{\Δ}{V}_g\right|\≤\mathrm{\Δ}{V}_g^{\max }\\ \left|\mathrm{\Δ}{V}_s\right|\≤V_s^{\max }\end{array}$

上述目标函数中，μ_s为动态无功补偿装置的无功功率矢量，表示动态无功补偿装置经控制后无功功率偏离设定值的大小，μ_g为发电机的无功功率矢量，表示发电机经控制后的无功功率偏离所有发电机无功功率平均值的大小；

${\mathrm{\μ}}_{\mathrm{si}}=\frac{Q_{\mathrm{si}}+\underset{j\∈\mathrm{Ns}}{\mathrm{\Σ}}{C}_{s_{\mathrm{ij}}}\mathrm{\Δ}{V}_{s_j}+\underset{k\∈\mathrm{Ng}}{\mathrm{\Σ}}{C}_{s{g}_{\mathrm{ik}}}\mathrm{\Δ}{V}_{g_k}-Q_{s_i}^{\mathrm{reference}}}{Q_{\mathrm{si}}^{\max }-Q_{\mathrm{si}}^{\min }}$

${\mathrm{\μ}}_{\mathrm{gi}}=\frac{Q_{\mathrm{gi}}+\underset{j\∈\mathrm{Ng}}{\mathrm{\Σ}}{C}_{g_{\mathrm{ij}}}{\mathrm{\ΔV}}_{g_j}+\underset{k\∈\mathrm{Ns}}{\mathrm{\Σ}}{C}_{g{s}_{\mathrm{ik}}}\mathrm{\Δ}{V}_{s_k}-Q_g^{\mathrm{aver}}}{Q_g^{\mathrm{aum}}}$

为动态无功补偿装置的裕度设定值，一般情况的取值为动态无功补偿装置无功功率的中间值，但是也可以根据电力系统中的实际情况略作调整，比如在可能大量缺少无功功率的地点，即可能发生严重故障的节点附近，设定的最优动态无功储备值可以取动态无功补偿装置出力中间偏小的位置，使得发生严重故障的时候此台动态无功补偿装置可以补偿更多的无功功率；为当前时刻所有发电机的无功功率平均值，为当前时刻发电机无功功率总和，W_s为动态无功补偿装置的无功功率矢量的权重系数，W_s的取值范围为0.1-1，W_g为发电机的无功功率矢量的权重系数，W_g的取值范围为10-100；

将计算得到的无功功率裕度优化阶段的电压调整量ΔV_g和ΔV_s作为控制指令下达给参与控制的发电机与动态无功储备装置，完成动态无功功率裕度优化阶段控制；

（4）在时刻T₀+△T，返回步骤（1）。时间T₀+ΔT（T₀为步骤（1）中所设置的时间起点，ΔT为本文所提出的中枢母线电压控制方法两次控制的时间间隔，一般取该数值为5至10分钟）。

Claims (1)

1.一种基于动态无功补偿设备的电力系统中中枢母线电压的控制方法，其特征在于该控制方法包括以下步骤：

(1)从电力系统的调度中心获取中枢母线电压的设定值当前采集时刻T₀的电力系统中中枢母线的电压值V_p，通过下式计算得到中枢母线电压控制阶段的参与控制的发电机的电压调整量ΔV_g和参与控制的动态无功补偿设备的电压调整量ΔV_s：

$\underset{{\mathrm{\ΔV}}_g,{\mathrm{\ΔV}}_s}{min}\left\{W_p\right||a\·(V_p-V_p^{ref})+C_{pg}{\mathrm{\ΔV}}_g+C_{ps}{\mathrm{\ΔV}}_s|{|}^2+W_q\left(\right|\left|{\mathrm{\Θ}}_g\right|{|}^2+\left|\right|{\mathrm{\Θ}}_s\left|{|}^2\right)\}$

受约束于： $V_p^{\min }\≤V_p+C_{pg}{\mathrm{\ΔV}}_g+C_{ps}{\mathrm{\ΔV}}_s\≤V_p^{max}$

$Q_g^{\min }\≤Q_g+C_g{\mathrm{\ΔV}}_g+C_{gs}{\mathrm{\ΔV}}_s\≤Q_g^{\max }$

$Q_s^{\min }\≤Q_s+C_{sg}{\mathrm{\ΔV}}_g+C_s{\mathrm{\ΔV}}_s\≤Q_s^{\max }$

$\left|{\mathrm{\ΔV}}_g\right|\≤{\mathrm{\ΔV}}_g^{max}$

$\left|{\mathrm{\ΔV}}_s\right|\≤{\mathrm{\ΔV}}_s^{\max }$

其中，Θ_g为发电机无功功率的裕度矢量，该裕度矢量Θ_g中的第i个分量为：

${\mathrm{\Θ}}_{g_i}=\frac{Q_{g_i}+\underset{j\∈Ng}{\Σ}{C}_{g_{ij}}{\mathrm{\ΔV}}_{g_j}+\underset{k\∈Ns}{\Σ}{C}_{{\mathrm{gs}}_{ik}}{\mathrm{\ΔV}}_{s_k}-Q_{g_i}^{\min }}{Q_{g_i}^{max}-Q_{g_i}^{\min }}$

Θ_s为动态无功补偿设备无功功率的裕度矢量，裕度矢量Θ_s的第i个分量为：

${\mathrm{\Θ}}_{s_i}=\frac{Q_{s_i}+\underset{j\∈Ns}{\Σ}{C}_{s_{ij}}{\mathrm{\ΔV}}_{s_j}+\underset{k\∈Ng}{\Σ}{C}_{{\mathrm{sg}}_{ik}}{\mathrm{\ΔV}}_{g_k}-Q_{s_i}^{\min }}{Q_{s_i}^{\max }-Q_{s_i}^{\min }}$

其中下标p表示中枢母线，下标g表示发电机，下标s表示动态无功补偿设备；C_pg是从调度中心获取的中枢母线电压偏差对参与控制的发电机的电压调整量ΔV_g的灵敏度矩阵；C_ps是从调度中心获取的中枢母线电压偏差对参与控制的动态无功补偿设备的电压调整量ΔV_s的灵敏度矩阵；W_p为中枢母线电压偏差的权重系数，W_p的取值范围为0.1-1，W_q为发电机和动态无功补偿设备的无功功率储备裕度的权重系数，W_q的取值范围为10-100，a是中枢母线电压的增益系数，a的取值范围为1-10，Q_g是参与控制的发电机当前发出的无功功率，分别为参与控制的发电机的无功功率下限与无功功率上限，Q_s是参与控制的动态无功补偿设备当前发出的无功功率，分别为参与控制的动态无功补偿设备发出的无功功率下限与无功功率的上限，C_g为从调度中心获取的参与控制的发电机的无功功率调整量对参与控制的发电机的电压调整量ΔV_g的灵敏度矩阵，C_s为从调度中心获取的参与控制的动态无功补偿设备的无功功率调整量对参与控制的动态无功补偿设备的电压调整量ΔV_s的灵敏度矩阵，C_gs为从调度中心获取的参与控制的发电机的无功功率调整量对参与控制的动态无功补偿设备的电压调整量ΔV_s的灵敏度矩阵，C_sg为从调度中心获取的参与控制的动态无功补偿设备的无功功率调整量对参与控制的发电机的电压调整量ΔV_g的灵敏度矩阵，为设定的中枢母线电压的下限，为设定的中枢母线电压上限，为从调度中心获取的参与控制的发电机允许的单步最大调整量，为从调度中心获取的参与控制的动态无功补偿设备允许的单步最大调整量，Ns为参与控制的动态无功补偿设备的编号的集合，Ng为参与控制的发电机的编号的集合；

将上述计算得到的参与控制的发电机的电压调整量ΔV_g和参与控制的动态无功补偿设备的电压调整量ΔV_s作为控制指令下达给相应的发电机和动态无功补偿设备，各发电机和动态无功补偿设备根据控制指令，完成中枢母线电压控制阶段；

(2)设定动态无功补偿设备无功功率的裕度阈值为将当前动态无功补偿设备的裕度值与设定的裕度阈值比较，若当前动态无功补偿设备的裕度值大于或等于设定的裕度阈值，则经过中枢母线电压两次控制之间的时间间隔ΔT后重复步骤(1)，ΔT的取值范围为5至10分钟，若当前动态无功补偿设备的裕度值小于设定的裕度阈值，则进行步骤(3)；

(3)根据下式计算，得到动态无功补偿设备动态无功功率裕度优化阶段的参与控制的发电机和动态无功补偿设备的电压调整量ΔV_g和ΔV_s：

$\underset{{\mathrm{\ΔV}}_g,{\mathrm{\ΔV}}_s}{min}\left\{W_s\right|\left|{\mathrm{\μ}}_s\right|{|}^2+W_g\left|\right|{\mathrm{\μ}}_g\left|{|}^2\right\}$

受约束于：V_p-V_band≤V_p+C_pgΔV_g+C_psΔV_s≤V_p+V_band

$Q_g^{\min }\≤Q_g+C_g{\mathrm{\ΔV}}_g+C_{gs}{\mathrm{\ΔV}}_s\≤Q_g^{\max }$

$Q_s^{\min }\≤Q_s+C_{sg}{\mathrm{\ΔV}}_g+C_s{\mathrm{\ΔV}}_s\≤Q_s^{\max }$

$\left|{\mathrm{\ΔV}}_g\right|\≤{\mathrm{\ΔV}}_g^{max}$

$\left|{\mathrm{\ΔV}}_s\right|\≤{\mathrm{\ΔV}}_s^{\max }$

上述目标函数中，μ_s为动态无功补偿设备的无功功率矢量，表示动态无功补偿设备经控制后无功功率偏离设定值的大小，μ_g为发电机的无功功率矢量，表示发电机经控制后的无功功率偏离所有发电机无功功率平均值的大小；

${\mathrm{\μ}}_{sr}=\frac{Q_{si}+\underset{j\∈Ns}{\Σ}{C}_{s_{ij}}{\mathrm{\ΔV}}_{s_j}+\underset{k\∈Ng}{\Σ}{C}_{{\mathrm{sg}}_{ik}}{\mathrm{\ΔV}}_{g_k}-Q_{s_i}^{reference}}{Q_{si}^{\max }-Q_{si}^{\min }}$

${\mathrm{\μ}}_{gi}=\frac{Q_{gi}+\underset{j\∈Ng}{\Σ}{C}_{g_{ij}}{\mathrm{\ΔV}}_{g_j}+\underset{k\∈Ns}{\Σ}{C}_{{\mathrm{gs}}_{ik}}{\mathrm{\ΔV}}_{s_k}-Q_g^{aver}}{Q_g^{sum}}$

为动态无功补偿设备的裕度设定值，为当前时刻所有发电机的无功功率平均值，为当前时刻发电机无功功率总和，W_s为动态无功补偿设备的无功功率矢量的权重系数，W_s的取值范围为0.1-1，W_g为发电机的无功功率矢量的权重系数，W_g的取值范围为10-100；

将计算得到的动态无功功率裕度优化阶段的电压调整量ΔV_g和ΔV_s作为控制指令下达给参与控制的发电机与动态无功补偿设备，完成动态无功功率裕度优化阶段控制；

(4)在时刻T₀+△T，返回步骤(1)。