CN104332991A - 电网潮流阻塞调度方法和电网潮流裕度评估方法 - Google Patents

Abstract

本发明提供了一种电网潮流阻塞调度方法，包括：根据所有时刻实时交易量标定值、线路潮流标定值和断面潮流标定值之和获得目标函数；根据机组出力上下限约束和爬坡约束获得机组侧约束；根据各省负荷平衡约束、各省间实时交易约束、系统备用约束、线路传输极限约束和断面传输极限约束获得电网侧约束；以所述机组侧约束和电网侧约束为约束条件，根据所述约束条件和目标函数获得电网潮流阻塞模型；在约束条件限制下，求解所述电网潮流阻塞模型以获得目标函数的最小值。本发明给出了电网潮流阻塞调度方法，确定了大功率失却时的备用支援辅助决策，本发明还提出了电网潮流裕度评估方法，完善了电网备用智能调用控制技术手段。

Description

电网潮流阻塞调度方法和电网潮流裕度评估方法

技术领域

本发明电网安全领域，尤其是一种电网潮流阻塞调度方法和电网潮流评估方法。

背景技术

华东为典型的大受端电网。预计到2015年，华东电网区外直流总额定容量达到31800MW，考虑网损后落地电力最大可达29892MW，其中±800kV特高压直流占比近70％。区外直流满送时，来电量将达到华东电网预计最高负荷的14.5％，其中特高压直流来电占最高负荷负荷的10％。

随着华东电网的快速发展，电网负荷增长迅猛，虽然电网建设逐渐受到重视，发展迅速仍落后于电源建设，使得在部分时段或电网特殊方式下阻塞情况时有发生。消除电网阻塞的根本措施应该为合理的电网、电源统一规划和精细化的电网运行方式安排。

目前，对于我国当前调度模式，电网发生阻塞时调度员依据优先消除阻塞的原则，直接调节相应机组出力来解决电网阻塞问题，这与国外电力市场进行阻塞管理有着很大的不同。电力市场下阻塞管理往往通过分区电价或节点电价指导机组重新报价、出清，以系统最低购电成本为目标。

在我国“三公”调度下，当发生电网阻塞时，首要目标为确保电网安全、稳定运行，应立即控制各断面、线路潮流在限额范围内，对于省内潮流越限，优先调用本省机组消除阻塞，对于省内机组消除阻塞有困难的，则网调应组织各省进行跨省交易、备用支援。为此，本章立足于国情、以我国当前统一调度为背景、建立以消除电网阻塞为目标的省间备用支援辅助决策。

发明内容

本发明的目的在于针对目前传统消除电网阻塞的措施，以最优化理论为基础，建立电网各省市潮流阻塞模型，构造目标函数，确定约束条件，给出电网阻塞发生时的调度方法，并提出一种电网潮流裕度评估方法，以衡量省市间实时交易、电能置换对潮流产生的影响以及选取有利于减轻潮流的机组出力安排方式，完善电网备用智能调用控制技术手段。

为了达到上述目的，本发明提供了一种电网潮流阻塞调度方法，包括以下步骤：

根据至少一省在当前时刻向其他省份额外购电交易量和其他省份在当前时刻向该省额外购电交易量之和再乘以第一标定系数获得实时交易量标定值，根据线路在当前时刻的正向越限量和反向越限量之和再乘以第二标定系数获得线路潮流标定值，根据断面在当前时刻的正向越限量和反向越限量之和再乘以第三标定系数获得断面潮流标定值，根据所有时刻所述实时交易量标定值、线路潮流标定值和断面潮流标定值之和获得目标函数；

根据机组出力上下限约束和爬坡约束获得机组侧约束，根据机组有功出力的上限、下限及机组在当前时刻的有功出力获得所述机组出力上下限约束，根据机组一时间间隔内所能升降的最大出力的上限和下限获得所述机组爬坡约束；

根据各省负荷平衡约束、各省间实时交易约束、系统备用约束、线路传输极限约束和断面传输极限约束获得电网侧约束，根据机组在当前时刻的有功出力、该省在当前时刻向其他省份额外购电交易量、其他省份在当前时刻向该省额外购电交易量和该省在当前时刻的等效负荷预测值获得所述各省负荷平衡约束，根据所有省份在当前时刻向其他省份额外购电交易量之和以及其他省份在当前时刻向该省额外购电交易量之和获得所述各省间实时交易约束，根据该省要留取的正备用量、负备用量、机组有功出力的上限、下限及机组在当前时刻的有功出力获得所述系统备用约束，根据该线路在当前时刻潮流、机组对该线路的功率转移分布因子、机组在当前时刻的有功出力、电网中节点对该线路的功率转移分布因子、电网中节点在当前时刻的节点负荷预测值、该线路的正向最大越限量和该线路的反向最大越限量获得所述线路传输极限约束，根据断面在当前时刻的潮流、机组对断面的功率转移分布因子、机组在当前时刻的有功出力、电网中节点对断面的功率转移分布因子、所述电网中节点在当前时刻的节点负荷预测值、该断面的正向最大越限量和该断面的反向最大越限量获得所述断面传输极限约束；

以所述机组侧约束和电网侧约束为约束条件，根据所述约束条件和目标函数获得电网潮流阻塞模型；

在约束条件限制下，求解所述电网潮流阻塞模型以使所述实时交易量标定值、线路潮流标定值和断面潮流标定值均为零，从而获得目标函数的最小值。

进一步地，所述目标函数为：

$\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{M}_e(e_{a,t}^{+}+e_{a,t}^{-})+\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{M}_l(\mathrm{\Δ}{P}_{l,t}^{+}+\mathrm{\Δ}{P}_{l.t}^{-})+\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{M}_s\left(\mathrm{\Δ}{S}_{s,t}^{+}\mathrm{\Δ}{S}_{s,t}^{-}\right),s\∈S,a\∈A,l\∈L,$

优化目标为所述目标函数的最小值，即：

$\min \{\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{M}_e(e_{a,t}^{+}+e_{a,t}^{-})+\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{M}_l(\mathrm{\Δ}{P}_{l,t}^{+}+\mathrm{\Δ}{P}_{l,t}^{-})+\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{M}_s(\mathrm{\Δ}{S}_{s,t}^{+}+\mathrm{\Δ}{S}_{s,t}^{-})\},s\∈S,a\∈A,l\∈L$

其中，当前时刻为t时刻，为a省在t时刻向其他省份额外购电交易量，为其他省份在t时刻向a省额外购电交易量，和分别代表线路l在t时刻的正向和反向越限量，和分别代表断面s在t时刻的正向和反向越限量，t＝1,...,T，T为决策周期内的总时段数，A为省份的集合，S为所有断面集合，L为所有线路集合，M_l、M_e、M_s分别为第一标定系数、第二标定系数和第三标定系数。

进一步地，所述第一标定系数、第二标定系数和第三标定系数均为正数，且：

M_l＝M_s＞M_e。

进一步地，所述机组出力上下限约束为：

P_i,min≤P_i,t≤P_i,max，

所述机组爬坡约束为：

$\mathrm{\Δ}\stackrel{\‾}{P_i}\≤P_{i,t}-P_{i,t-1}\≤\mathrm{\Δ}\underset{\‾}{P_i},$

其中，i＝1,2,...,N，P_i,t和P_i,t-1分别为机组i在t时刻及t-1时刻的有功出力；P_i,min和P_i,max分别为机组i的出力下限和上限，和表示机组i在t时刻和t-1时刻这一个时间间隔内所能升降的最大出力范围。

进一步地，所述各省负荷平衡约束为：

$\mathrm{\Σ}{P}_{i,t}+e_{a,t}^{+}-e_{a,t}^{-}=D_{a,t},$

其中，i∈A_a，D_a,t为a省时段t的等效负荷预测值，A_a为a省机组集合，

所述各省间实时交易约束为：

$\underset{a}{\overset{A_a}{\mathrm{\Σ}}}{e}_{a,t}^{-}=\underset{a}{\overset{A_a}{\mathrm{\Σ}}}{e}_{a,t}^{+}$

$0\≤e_{a,t}^{-}\≤e_{a,\max }^{-}$

$0\≤e_{a.t}^{+}\≤e_{a,\max }^{+}$

其中，和为省间交易的最大量，

所述系统备用约束为：

$\underset{i=1}{\mathrm{\Σ}}(P_{i,\max }-P_{i,t})\≥r_{a,p},\underset{i=1}{\mathrm{\Σ}}(P_{i,t}-P_{i,\min })\≤r_{a,n},$

其中，r_a,p为规定a省要留取的正备用量，r_a,n为规定a省要留取的负备用量，

所述线路传输极限约束为：

$P_{l,t}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{l-i}{P}_{i,t}-\underset{k=1}{\overset{K}{\mathrm{\Σ}}}{G}_{l-k}{d}_{k,t}$

$P_{l,t}-\mathrm{\Δ}{P}_{l,t}^{+}\≤P_{l,\max }$

$P_{l,\min }\≤P_{l,t}+\mathrm{\Δ}{P}_{l,t}^{-}$

$0\≤P_{l,t}^{+}\≤{\mathrm{\ΔP}}_{l,\max }^{+}$

$0\≤P_{l,t}^{-}\≤{\mathrm{\ΔP}}_{l,\max }^{-}$

其中，l∈L,k＝1,2,...,K，P_l,t为线路l在时段t潮流，P_l,max和P_l,min分别为线路l的正传输功率极限和反传输功率极限，G_l-i为机组i对线路l的功率转移分布因子，G_l-k为电网中节点k对线路l的功率转移分布因子，d_k,t为电网中节点k在t时刻的节点负荷预测值，K为所有节点集合，L为所有线路集合；和分别为允许线路l反向的最大越限量和正向的最大越限量，

所述断面传输极限约束为：

$S_{s,t}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{s-i}{P}_{i,t}-\underset{k=1}{\overset{K}{\mathrm{\Σ}}}{G}_{s-k}{d}_{k,t}$

$S_{s,t}-\mathrm{\Δ}{S}_{s,t}^{+}\≤S_{s,\max }$

$S_{s,\min }\≤S_{s,t}+\mathrm{\Δ}{S}_{s,t}^{-}$

$0\≤\mathrm{\Δ}{S}_{s,t}^{+}\≤\mathrm{\Δ}{S}_{s,\max }^{+}$

$0\≤\mathrm{\Δ}{S}_{s,t}^{-}\≤\mathrm{\Δ}{S}_{s,\max }^{-}$

其中，i＝1,2,...,N，t＝1,...,T，s∈S，k＝1,2,...,K，G_s-i为机组i对断面s的功率转移分布因子，G_s-k为电网中节点k对断面s的功率转移分布因子，S_s,t为断面s在时段t的潮流，和分别为允许的最大越限量。

进一步地，在省内大机组故障时，根据所述电网侧约束和优化目标求解所述电网潮流阻塞模型以获得目标函数的最小值。

进一步地，在区外直流大功率失却时，即区外直流大功率为零时，所述各省负荷平衡约束为：

$\mathrm{\Σ}{P}_{i,t}+e_{a,t}^{+}-e_{a,t}^{-}=D_{a,t}+\mathrm{\Δ}{D}_{a,t},$

其中，ΔD_a,t为分摊到省市的口子变化量。

本发明还提供了一种电网潮流裕度评估方法，包括以下步骤：

根据断面在当前时刻的潮流和断面的正向最大越限量获得断面的目标系数，根据线路在当前时刻的潮流和线路的正向最大越限量获得线路的目标系数，根据断面的目标系数、机组对断面的功率转移分布因子、机组在当前时刻的有功出力、线路的目标系数和机组对该线路的功率转移分布因子获得潮流阻塞最严重情况下的目标函数以及潮流阻塞最佳情况下的目标函数；

根据机组出力上下限约束和爬坡约束获得机组侧约束，根据机组有功出力的上限、下限及机组在当前时刻的有功出力获得所述机组出力上下限约束，根据机组一时间间隔内所能升降的最大出力的上限和下限获得所述机组爬坡约束；

根据机组在当前时刻的有功出力、该省在当前时刻向其他省份额外购电交易量、其他省份在当前时刻向该省额外购电交易量和该省在当前时刻的等效负荷预测值获得各省负荷平衡约束；

以所述机组侧约束、各省负荷平衡约束为约束条件，根据所述约束条件和潮流阻塞最严重情况下的目标函数获得第一评估模型，根据所述约束条件和潮流阻塞最佳情况下的目标函数获得第二评估模型；

在约束条件限制下，求解所述第一评估模型以获得潮流阻塞最严重情况下的目标函数的最小值，求解所述第二评估模型以获得潮流阻塞最佳情况下的目标函数的最小值。

进一步地，所述潮流阻塞最严重情况下的目标函数为：

$-\underset{t=1}{\overset{T}{\mathrm{\Σ}}}(C_s\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{s-i}{P}_{i,t}+C_l^{/}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{l-i}{P}_{i,t}),s\∈S,l\∈L,$

第一优化目标为所述潮流阻塞最严重情况下的目标函数的最小值，即：

$\min \{-\underset{t=1}{\overset{T}{\mathrm{\Σ}}}(C_s\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{s-i}{P}_{i,t}+C_l^{/}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{l-i}{P}_{i,t})\},s\∈S,l\∈L$

所述潮流阻塞最佳情况下的目标函数为：

$\underset{t=1}{\overset{T}{\mathrm{\Σ}}}(C_s\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{s-i}{P}_{i,t}+C_l^{/}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{l-i}{P}_{i,t}),s\∈S,l\∈L,$

第二优化目标为所述潮流阻塞最佳情况下的目标函数的最小值，即：

$\min \left\{\underset{t=1}{\overset{T}{\mathrm{\Σ}}}(C_s\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{s-i}{P}_{i,t}+C_l^{/}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{l-i}{P}_{i,t})\right\},s\∈S,l\∈L$

其中，i＝1,2,...,N，i∈A_a,t＝1,...,T，l∈L,k＝1,2,...,K，A_a为a省机组集合，T为决策周期内的总时段数，C_s和分别为断面s和线路l的目标系数，S_s,t为所述断面在当前时刻的潮流，S_s,max为所述断面的正向最大越限量，当前时刻为t时刻，S为所有断面集合，L为所有线路集合，P_i,t为机组i在t时刻的有功出力，G_l-i为机组i对线路l的功率转移分布因子，G_s-i为机组i对断面s的功率转移分布因子。

进一步地，所述机组出力上下限约束为：

P_i,min≤P_i,t≤P_i,max，

所述机组爬坡约束为：

$\mathrm{\Δ}\stackrel{\‾}{P_i}\≤P_{i,t}-P_{i,t-1}\≤\mathrm{\Δ}\underset{\‾}{P_i},$

所述各省负荷平衡约束为：

$\mathrm{\Σ}{P}_{i,t}+e_{a,t}^{+}-e_{a,t}^{-}=D_{a,t},$

其中，i＝1,2,...,N，i∈A_a,i∈A_a,t＝1,...,T，P_i,t和P_i,t-1分别为机组i在t时刻及t-1时刻的有功出力；P_i,min和P_i,max分别为机组i的出力下限和上限，和表示机组i在t时刻和t-1时刻这一个时间间隔内所能升降的最大出力范围，D_a,t为a省时段t的等效负荷预测值。

与现有技术相比，本发明以最优化理论为基础，建立电网各省市潮流阻塞模型，构造目标函数，确定约束条件，给出了一种电网潮流阻塞调度方法。特别是针对华东电网大功率失却后潮流发生转移将会造成潮流越限情况，通过求解上述建立模型，确定了大功率失却时的备用支援辅助决策。同时，本发明还提出了一种电网潮流裕度评估方法，衡量省市间实时交易、电能置换对潮流产生的影响以及选取有利于减轻潮流的机组出力安排方式，完善了电网备用智能调用控制技术手段。

附图说明

图1为本发明实施例提供的电网潮流阻塞调度方法流程图；

图2为本发明实施例提供的断面的目标系数与断面潮流负载率的关系图；

图3为本发明实施例提供的线路的目标系数与线路潮流负载率的关系图。

具体实施方式

下面将结合示意图对本发明的具体实施方式进行更详细的描述。根据下列描述和权利要求书，本发明的优点和特征将更清楚。需说明的是，附图均采用非常简化的形式且均使用非精准的比例，仅用以方便、明晰地辅助说明本发明实施例的目的。

如图1所示，本发明提供了一种电网潮流阻塞调度方法，包括以下步骤：

根据至少一省在当前时刻向其他省份额外购电交易量和其他省份在当前时刻向该省额外购电交易量之和再乘以第一标定系数获得实时交易量标定值，根据线路在当前时刻的正向越限量和反向越限量之和再乘以第二标定系数获得线路潮流标定值，根据断面在当前时刻的正向越限量和反向越限量之和再乘以第三标定系数获得断面潮流标定值，根据所有时刻所述实时交易量标定值、线路潮流标定值和断面潮流标定值之和获得目标函数；

根据机组出力上下限约束和爬坡约束获得机组侧约束，根据机组有功出力的上限、下限及机组在当前时刻的有功出力获得所述机组出力上下限约束，根据机组一时间间隔内所能升降的最大出力的上限和下限获得所述机组爬坡约束；

根据各省负荷平衡约束、各省间实时交易约束、系统备用约束、线路传输极限约束和断面传输极限约束获得电网侧约束，根据机组在当前时刻的有功出力、该省在当前时刻向其他省份额外购电交易量、其他省份在当前时刻向该省额外购电交易量和该省在当前时刻的等效负荷预测值获得所述各省负荷平衡约束，根据所有省份在当前时刻向其他省份额外购电交易量之和以及其他省份在当前时刻向该省额外购电交易量之和获得所述各省间实时交易约束，根据该省要留取的正备用量、负备用量、机组有功出力的上限、下限及机组在当前时刻的有功出力获得所述系统备用约束，根据该线路在当前时刻潮流、机组对该线路的功率转移分布因子、机组在当前时刻的有功出力、电网中节点对该线路的功率转移分布因子、电网中节点在当前时刻的节点负荷预测值、该线路的正向最大越限量和该线路的反向最大越限量获得所述线路传输极限约束，根据断面在当前时刻的潮流、机组对断面的功率转移分布因子、机组在当前时刻的有功出力、电网中节点对断面的功率转移分布因子、所述电网中节点在当前时刻的节点负荷预测值、该断面的正向最大越限量和该断面的反向最大越限量获得所述断面传输极限约束；

以所述机组侧约束和电网侧约束为约束条件，根据所述约束条件和目标函数获得电网潮流阻塞模型；

在约束条件限制下，求解所述电网潮流阻塞模型以使所述实时交易量标定值、线路潮流标定值和断面潮流标定值均为零，从而获得目标函数的最小值。

所述目标函数为：

$\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{M}_e(e_{a,t}^{+}+e_{a,t}^{-})+\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{M}_l(\mathrm{\Δ}{P}_{l,t}^{+}+\mathrm{\Δ}{P}_{l.t}^{-})+\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{M}_s\left(\mathrm{\Δ}{S}_{s,t}^{+}\mathrm{\Δ}{S}_{s,t}^{-}\right),s\∈S,a\∈A,l\∈L,$

优化目标为所述目标函数的最小值，即：

$\min \{\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{M}_e(e_{a,t}^{+}+e_{a,t}^{-})+\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{M}_l(\mathrm{\Δ}{P}_{l,t}^{+}+\mathrm{\Δ}{P}_{l,t}^{-})+\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{M}_s(\mathrm{\Δ}{S}_{s,t}^{+}+\mathrm{\Δ}{S}_{s,t}^{-})\},s\∈S,a\∈A,l\∈L$

该优化目标的含义为通过各种手段使得各省市口子计划维持在给定值 $(e_{a,t}^{-}=0,e_{a,t}^{+}=0),$ 网内各线路潮流不越限 $({\mathrm{\ΔP}}_{l,t}^{+}=0,{\mathrm{\ΔP}}_{l,t}^{-}=0),$ 各断面潮流不越限额此时目标函数值为0，达到最优。

其中，当前时刻为t时刻，为a省在t时刻向其他省份额外购电交易量，为其他省份在t时刻向a省额外购电交易量，和分别代表线路l在t时刻的正向和反向越限量，和分别代表断面s在t时刻的正向和反向越限量，t＝1,...,T，T为决策周期内的总时段数，A为省份的集合，S为所有断面集合，L为所有线路集合，M_l、M_e、M_s分别为第一标定系数、第二标定系数和第三标定系数。

所述第一标定系数、第二标定系数和第三标定系数均为正数，均用于协调控制个各类目标的关系，电网应优先消除阻塞，则：

M_l＝M_s＞＞M_e。

所述机组出力上下限约束为：

P_i,min≤P_i,t≤P_i,max，

所述机组爬坡约束为：

$\mathrm{\Δ}\stackrel{\‾}{P_i}\≤P_{i,t}-P_{i,t-1}\≤\mathrm{\Δ}\underset{\‾}{P_i},$

其中，i＝1,2,...,N，t＝1,...,T，P_i,t和P_i,t-1分别为机组i在t时刻及t-1时刻的有功出力；P_i,min和P_i,max分别为机组i的出力下限和上限，和表示机组i在t时刻和t-1时刻这一个时间间隔内所能升降的最大出力范围。

t时刻和t-1时刻这一个时间间隔的长度是可设定的，比如考察t＝1和t＝2这两个时刻的时间间隔，将其设置为15分钟，且设置t＝1时刻为9时，那么t＝1时刻为9时15分，以此类推，在此不再赘述。

所述各省负荷平衡约束为：

$\mathrm{\Σ}{P}_{i,t}+e_{a,t}^{+}-e_{a,t}^{-}=D_{a,t},$

其中，i∈A_a,t＝1,...,T，D_a,t为a省时段t的等效负荷预测值(考虑口子计划，所述口子计划为某省和其他省交换的功率值)，Aa为a省机组集合，增加变量和是为了防止在求解过程中出现不可行解。

所述各省间实时交易约束为：

$\underset{a}{\overset{A_a}{\mathrm{\Σ}}}{e}_{a,t}^{-}=\underset{a}{\overset{A_a}{\mathrm{\Σ}}}{e}_{a,t}^{+}$

$0\≤e_{a,t}^{-}\≤e_{a,\max }^{-}$

$0\≤e_{a.t}^{+}\≤e_{a,\max }^{+}$

各省间实时交易约束公式的含义是各省市交易量之和应相等，其中，和为省间交易的最大量。

所述系统备用约束为：

$\underset{i=1}{\mathrm{\Σ}}(P_{i,\max }-P_{i,t})\≥r_{a,p},\underset{i=1}{\mathrm{\Σ}}(P_{i,t}-P_{i,\min })\≤r_{a,n},$

其中，i∈A_a,t＝1,...,T，r_a,p为规定a省要留取的正备用量，r_a,n为规定a省要留取的负备用量。

所述线路传输极限约束为：

$P_{l,t}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{l-i}{P}_{i,t}-\underset{k=1}{\overset{K}{\mathrm{\Σ}}}{G}_{l-k}{d}_{k,t}$

$P_{l,t}-\mathrm{\Δ}{P}_{l,t}^{+}\≤P_{l,\max }$

$P_{l,\min }\≤P_{l,t}+\mathrm{\Δ}{P}_{l,t}^{-}$

$0\≤P_{l,t}^{+}\≤{\mathrm{\ΔP}}_{l,\max }^{+}$

$0\≤P_{l,t}^{-}\≤{\mathrm{\ΔP}}_{l,\max }^{-}$

其中，i＝1,2,...,N，t＝1,...,T，l∈L,k＝1,2,...,K，P_l,t为线路l在时段t潮流，P_l,max和P_l,min分别为线路l的正传输功率极限和反传输功率极限，G_l-i为机组i对线路l的功率转移分布因子，G_l-k为电网中节点k对线路l的功率转移分布因子，d_k,t为电网中节点k在t时刻的节点负荷预测值，K为所有节点集合，L为所有线路集合；和分别为允许线路l反向的最大越限量和正向的最大越限量，

所述断面传输极限约束为：

$S_{s,t}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{s-i}{P}_{i,t}-\underset{k=1}{\overset{K}{\mathrm{\Σ}}}{G}_{s-k}{d}_{k,t}$

$S_{s,t}-\mathrm{\Δ}{S}_{s,t}^{+}\≤S_{s,\max }$

$S_{s,\min }\≤S_{s,t}+\mathrm{\Δ}{S}_{s,t}^{-}$

$0\≤\mathrm{\Δ}{S}_{s,t}^{+}\≤\mathrm{\Δ}{S}_{s,\max }^{+}$

$0\≤\mathrm{\Δ}{S}_{s,t}^{-}\≤\mathrm{\Δ}{S}_{s,\max }^{-}$

其中，i＝1,2,...,N，t＝1,...,T，s∈S，k＝1,2,...,K，G_s-i为机组i对断面s的功率转移分布因子，G_s-k为电网中节点k对断面s的功率转移分布因子，S_s,t为断面s在时段t的潮流，和分别为允许的最大越限量。

按照上述方法建立的模型为线性规划问题，该模型可以采取求解线性规划的一切算法进行计算，比如单纯型法、对偶单纯型法、内点法、可行方向法、惩罚函数法等方法。

因此，本发明提供了一种电网潮流阻塞调度方法，具体地说，通过构造目标函数、确定约束条件建立了电网各省市阻塞模型(即所述电网潮流阻塞模型)，给出了阻塞发生时的调度方法(根据该模型求解出的Pi,t、改变P_l,t和S_s,t，以使得线路潮流和断面潮流均在限额以内)。

而大机组跳闸、区外直流闭锁都会造成省市大功率失却，需要其他省市进行支援，同时大功率失却后潮流发生转移将会造成潮流越限等情况。这些情况都可以通过求解上述建立的模型，这里分大机组故障和区外直流闭锁两种情况进行讨论。

省内大机组发生故障后，首先应由本省省内机组加出力，若仍不能平衡负荷或消除阻塞则需要其他省市支援。对于传统的经济调度模型、阻塞管理模型未对负荷平衡约束进行松弛，若省内机组无法满足省内负荷，将无法求得可行解，这显然与实际情况不符，而采用本模型由于对负荷平衡约束和潮流约束进行了松弛，增加了和和和使得该问题存在可行解。这些变量不仅是求解优化问题的松弛变量而且有着十分明确的物理意义。

因此，在省内大机组故障时，根据所述电网侧约束和优化目标求解所述电网潮流阻塞模型以获得目标函数的最小值。

和代表了省市间需要进行的备用支援，如果和不为零，区域调度中心可以根据和进行相应的口子计划调整(即提供了相应的备用支援辅助决策)，例如和均为某给定量时，则区域调度中心通过a省在t时刻向其他省份额外购买该给定量的电以及其他省份在t时刻向a省额外出售该给定量的电实现口子计划调整。

和和代表潮流的越限量，若不为0，则说明即使通过调节全网的机组出力也不能进行消除潮流阻塞，此时需要调度员采取负荷控制、拍停机等其他措施(即提供了相应的备用支援辅助决策)，提前做好准备。

区外直流故障单级故障、双极闭锁等都会引起区外直流大功率失却。当直流大功率失却后将由动态ACE按照直流分摊比例，实时修改省市口子计划。此时，也应该将修改相应的各省负荷平衡约束，即区外直流大功率为零时，修改后的所述各省负荷平衡约束为：

$\mathrm{\Σ}{P}_{i,t}+e_{a,t}^{+}-e_{a,t}^{-}=D_{a,t}+\mathrm{\Δ}{D}_{a,t},$

其中，i∈A_a,t＝1,...,T，ΔD_a,t为分摊到省市的口子变化量。

区外直流大功率失却还容易造成潮流的大范围转移，因此应该重点关注潮流阻塞问题。上述优化目标中对潮流的越限量进行了相应的标定，若潮流越限量越大则标定的越多(标定的越多，则代表潮流越限越多，该省电网阻塞越大，电网性能越差)，因此自动实现了优先调节越限量大的线路或断面。

当省市发生负荷缺额、潮流阻塞时，根据备用共享原则，省市之间通过电能置换、实时交易进行备用支援确保电网安全、稳定、有序运行。省市之间进行的备用支援，有可能减轻线路、断面潮流，也可能加重潮流。如何衡量省市间实时交易、电能置换对潮流产生的影响以及选取有利于减轻潮流的机组出力安排方式是需要解决的问题。

设i省向j省在时段t支援备用，i省可以通过增加省内任意机组出力实现对j省的支援，j省则可以降低本省任意机组出力来实现口子调整，这样的组合有无穷多个，潮流分布影响也有无穷多个。我们将该备用支援最大程度加重线路、断面潮流甚至阻塞作为一种最严重情况，由此求得的机组出力向量P_b和潮流结果P_l,b，另一方面，将该备用支援最大程度减轻线路、断面潮流，使得各潮流分部均匀作为最好的结果，由此求得的机组出力向量P_g和潮流结果P_l,g，由此得到两种极端情况下的潮流结果，作为衡量该备用支援的潮流裕度，其他备用支援造成的结果都应该涵盖在这两种极端情况之间。

由此，本发明还提供了一种电网潮流的评估方法，包括以下步骤：

根据断面在当前时刻的潮流和断面的正向最大越限量获得断面的目标系数，根据线路在当前时刻的潮流和线路的正向最大越限量获得线路的目标系数，根据断面的目标系数、机组对断面的功率转移分布因子、机组在当前时刻的有功出力、线路的目标系数和机组对该线路的功率转移分布因子获得潮流阻塞最严重情况下的目标函数以及潮流阻塞最佳情况下的目标函数；

根据机组出力上下限约束和爬坡约束获得机组侧约束，根据机组有功出力的上限、下限及机组在当前时刻的有功出力获得所述机组出力上下限约束，根据机组一时间间隔内所能升降的最大出力的上限和下限获得所述机组爬坡约束；

根据机组在当前时刻的有功出力、该省在当前时刻向其他省份额外购电交易量、其他省份在当前时刻向该省额外购电交易量和该省在当前时刻的等效负荷预测值获得各省负荷平衡约束；

以所述机组侧约束、各省负荷平衡约束为约束条件，根据所述约束条件和潮流阻塞最严重情况下的目标函数获得第一评估模型，根据所述约束条件和潮流阻塞最佳情况下的目标函数获得第二评估模型；

在约束条件限制下，求解所述第一评估模型以获得潮流阻塞最严重情况下的目标函数的最小值，求解所述第二评估模型以获得潮流阻塞最佳情况下的目标函数的最小值。

为求得潮流阻塞最严重情况，在给定的省间实时交易情况下，最大化重载线路或断面的潮流，优化目标如下：

$\max \underset{t=1}{\overset{T}{\mathrm{\Σ}}}(C_s\mathrm{\Δ}{S}_{s,t}+C_l^{/}\mathrm{\Δ}{P}_{l,t})s\∈S,l\∈L$

其中ΔS_s,t＝S_s,t-S_s,max，ΔP_l,t＝P_l,t-P_l,max，C_s和分别为断面s和线路l的目标系数。

C_s和断面潮流负载率(S_s,t/S_s,max)关系如图2所示，和线路潮流负载率(P_s,t/P_s,max)关系如图3所示。此外，在优化过程中，优先考虑潮流负载率重的断面，让其先达到越限。

将上述优化目标中断面和线路潮流线性化，表示为以机组出力为决策变量，则：

$\mathrm{\Δ}{S}_{s,t}=S_{s,t}-S_{s,\max }=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{s-i}{P}_{i,t}-d_{s,t}-S_{s,\max }$

$\mathrm{\Δ}{L}_{l,t}=L_{l,t}-L_{l,\max }=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{l-i}{P}_{i,t}-d_{l,t}-S_{l,\max }$

可知，所述潮流阻塞最严重情况下的目标函数为：

$-\underset{t=1}{\overset{T}{\mathrm{\Σ}}}(C_s\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{s-i}{P}_{i,t}+C_l^{/}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{l-i}{P}_{i,t}),s\∈S,l\∈L,$

第一优化目标为所述潮流阻塞最严重情况下的目标函数的最小值，即：

$\min \{-\underset{t=1}{\overset{T}{\mathrm{\Σ}}}(C_s\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{s-i}{P}_{i,t}+C_l^{/}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{l-i}{P}_{i,t})\},s\∈S,l\∈L$

为求得潮流分布最佳情况，在给定的省间实时交易情况下，最小化重载断面和线路潮流，优化目标如下：

$\min \underset{t=1}{\overset{T}{\mathrm{\Σ}}}(C_s\mathrm{\Δ}{S}_{s,t}+C_l^{/}\mathrm{\Δ}{L}_{l,t})s\∈S,l\∈L$

将该优化目标表示成以发电机出力为决策变量，则所述潮流阻塞最佳情况下的目标函数为：

$\underset{t=1}{\overset{T}{\mathrm{\Σ}}}(C_s\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{s-i}{P}_{i,t}+C_l^{/}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{l-i}{P}_{i,t}),s\∈S,l\∈L,$

第二优化目标为所述潮流阻塞最佳情况下的目标函数的最小值，即：

$\min \left\{\underset{t=1}{\overset{T}{\mathrm{\Σ}}}(C_s\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{s-i}{P}_{i,t}+C_l^{/}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{l-i}{P}_{i,t})\right\},s\∈S,l\∈L$

其中，i＝1,2,...,N，t＝1,...,T，l∈L,k＝1,2,...,K，A_a为a省机组集合，T为决策周期内的总时段数，C_s和分别为断面s和线路l的目标系数，S_s,t为所述断面在当前时刻的潮流，S_s,max为所述断面的正向最大越限量，当前时刻为t时刻，S为所有断面集合，L为所有线路集合，P_i,t为机组i在t时刻的有功出力，G_l-i为机组i对线路l的功率转移分布因子，G_s-i为机组i对断面s的功率转移分布因子。

潮流最佳情况中约束条件和潮流阻塞最严重情况中约束条件相同，包含机组爬坡约束、机组出力上下限约束和负荷平衡约束，且：

所述机组出力上下限约束为：

P_i,min≤P_i,t≤P_i,max，

所述机组爬坡约束为：

$\mathrm{\Δ}\stackrel{\‾}{P_i}\≤P_{i,t}-P_{i,t-1}\≤\mathrm{\Δ}\underset{\‾}{P_i},$

所述各省负荷平衡约束为：

$\mathrm{\Σ}{P}_{i,t}+e_{a,t}^{+}-e_{a,t}^{-}=D_{a,t},$

其中，i＝1,2,...,N，i∈A_a,t＝1,...,T，P_i,t和P_i,t-1分别为机组i在t时刻及t-1时刻的有功出力；P_i,min和P_i,max分别为机组i的出力下限和上限，和表示机组i在t时刻和t-1时刻这一个时间间隔内所能升降的最大出力范围，D_a,t为a省时段t的等效负荷预测值。

通过求解上述潮流阻塞最严重情况和潮流分布最佳情况数学模型得到相应的机组出力和潮流结果。其中潮流分布最佳情况中的机组出力可以指导各个省市按照该结果调节本省机组出力，减轻重载潮流，而潮流阻塞最严重情况作为一个预警指标，让调度员对本次备用支援造成的结果有个提前的了解(即实现了采用所述电网潮流裕度评估方法实现了电网潮流裕度评估)。

综上，与现有技术相比，本发明以最优化理论为基础，建立电网各省市潮流阻塞模型，构造目标函数，确定约束条件，给出了一种电网潮流阻塞调度方法。特别是针对华东电网大功率失却后潮流发生转移将会造成潮流越限情况，通过求解上述建立模型，确定了大功率失却时的备用支援辅助决策。同时，本发明还提出了一种电网潮流裕度评估方法，衡量省市间实时交易、电能置换对潮流产生的影响以及选取有利于减轻潮流的机组出力安排方式，完善了电网备用智能调用控制技术手段。

上述仅为本发明的优选实施例而已，并不对本发明起到任何限制作用。任何所属技术领域的技术人员，在不脱离本发明的技术方案的范围内，对本发明揭露的技术方案和技术内容做任何形式的等同替换或修改等变动，均属未脱离本发明的技术方案的内容，仍属于本发明的保护范围之内。

Claims (10)

1.一种电网潮流阻塞调度方法，其特征在于，包括以下步骤：

根据至少一省在当前时刻向其他省份额外购电交易量和其他省份在当前时刻向该省额外购电交易量之和再乘以第一标定系数获得实时交易量标定值，根据线路在当前时刻的正向越限量和反向越限量之和再乘以第二标定系数获得线路潮流标定值，根据断面在当前时刻的正向越限量和反向越限量之和再乘以第三标定系数获得断面潮流标定值，根据所有时刻所述实时交易量标定值、线路潮流标定值和断面潮流标定值之和获得目标函数；

根据机组出力上下限约束和爬坡约束获得机组侧约束，根据机组有功出力的上限、下限及机组在当前时刻的有功出力获得所述机组出力上下限约束，根据机组一时间间隔内所能升降的最大出力的上限和下限获得所述机组爬坡约束；

根据各省负荷平衡约束、各省间实时交易约束、系统备用约束、线路传输极限约束和断面传输极限约束获得电网侧约束，根据机组在当前时刻的有功出力、该省在当前时刻向其他省份额外购电交易量、其他省份在当前时刻向该省额外购电交易量和该省在当前时刻的等效负荷预测值获得所述各省负荷平衡约束，根据所有省份在当前时刻向其他省份额外购电交易量之和以及其他省份在当前时刻向该省额外购电交易量之和获得所述各省间实时交易约束，根据该省要留取的正备用量、负备用量、机组有功出力的上限、下限及机组在当前时刻的有功出力获得所述系统备用约束，根据该线路在当前时刻潮流、机组对该线路的功率转移分布因子、机组在当前时刻的有功出力、电网中节点对该线路的功率转移分布因子、电网中节点在当前时刻的节点负荷预测值、该线路的正向最大越限量和该线路的反向最大越限量获得所述线路传输极限约束，根据断面在当前时刻的潮流、机组对断面的功率转移分布因子、机组在当前时刻的有功出力、电网中节点对断面的功率转移分布因子、所述电网中节点在当前时刻的节点负荷预测值、该断面的正向最大越限量和该断面的反向最大越限量获得所述断面传输极限约束；

以所述机组侧约束和电网侧约束为约束条件，根据所述约束条件和目标函数获得电网潮流阻塞模型；

在约束条件限制下，求解所述电网潮流阻塞模型以使所述实时交易量标定值、线路潮流标定值和断面潮流标定值均为零，从而获得目标函数的最小值。

2.如权利要求1所述的电网潮流阻塞调度方法，其特征在于，所述目标函数为：

$\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{M}_e(e_{a,t}^{+}+e_{a,t}^{-})+\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{M}_l({\mathrm{\ΔP}}_{l,t}^{+}+{\mathrm{\ΔP}}_{l,t}^{-})+\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{M}_s({\mathrm{\ΔS}}_{s,t}^{+}+{\mathrm{\ΔS}}_{s,t}^{-}),s\∈S,a\∈A,l\∈L,$

优化目标为所述目标函数的最小值，即：

$\min \{\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{M}_e(e_{a,t}^{+}+e_{a,t}^{-})+\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{M}_l({\mathrm{\ΔP}}_{l,t}^{+}+{\mathrm{\ΔP}}_{l,t}^{-})+\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{M}_s({\mathrm{\ΔS}}_{s,t}^{+}+{\mathrm{\ΔS}}_{s,t}^{-})\},s\∈S,a\∈A,l\∈L$

其中，当前时刻为t时刻，为a省在t时刻向其他省份额外购电交易量，为其他省份在t时刻向a省额外购电交易量，和分别代表线路l在t时刻的正向和反向越限量，和分别代表断面s在t时刻的正向和反向越限量，t＝1,...,T，T为决策周期内的总时段数，A为省份的集合，S为所有断面集合，L为所有线路集合，M_l、M_e、M_s分别为第一标定系数、第二标定系数和第三标定系数。

3.如权利要求2所述的电网潮流阻塞调度方法，其特征在于，所述第一标定系数、第二标定系数和第三标定系数均为正数，且：

M_l＝M_s＞M_e。

4.如权利要求3所述的电网潮流阻塞调度方法，其特征在于，所述机组出力上下限约束为：

P_i,min≤P_i,t≤P_i,max，

所述机组爬坡约束为：

$\mathrm{\Δ}{\stackrel{\‾}{P}}_i\≤P_{i,t}-P_{i,t-1}\≤\mathrm{\Δ}\underset{\‾}{P_i},$

其中，i＝1,2,...,N，P_i,t和P_i,t-1分别为机组i在t时刻及t-1时刻的有功出力；P_i,min和P_i,max分别为机组i的出力下限和上限，和表示机组i在t时刻和t-1时刻这一个时间间隔内所能升降的最大出力范围。

5.如权利要求3所述的电网潮流阻塞调度方法，其特征在于，所述各省负荷平衡约束为：

$\mathrm{\Σ}{P}_{i,t}+e_{a,t}^{+}-e_{a,t}^{-}=D_{a,t},$

其中，i∈A_a，D_a,t为a省时段t的等效负荷预测值，A_a为a省机组集合，

所述各省间实时交易约束为：

$\underset{a}{\overset{A_a}{\mathrm{\Σ}}}{e}_{a,t}^{-}=\underset{a}{\overset{A_a}{\mathrm{\Σ}}}{e}_{a,t}^{+}$

$0\≤e_{a,t}^{-}\≤w_{a,\max }^{-}$

$0\≤e_{a.t}^{+}\≤e_{a,\max }^{+}$

其中，和为省间交易的最大量，

所述系统备用约束为：

$\underset{i=1}{\mathrm{\Σ}}(P_{i,\max }-P_{i,t})\≥r_{a,p},\underset{i=1}{\mathrm{\Σ}}(P_{i,t}-P_{i,\min })\≤r_{a,n},$

其中，r_a,p为规定a省要留取的正备用量，r_a,n为规定a省要留取的负备用量，

所述线路传输极限约束为：

$P_{l,t}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{l-i}{P}_{i,t}-\underset{k=1}{\overset{K}{\mathrm{\Σ}}}{G}_{l-k}{d}_{k,t}$

$P_{l,t}-{\mathrm{\ΔP}}_{l,t}^{+}\≤P_{l,\max }$

$P_{l,\min }\≤P_{l,t}+{\mathrm{\ΔP}}_{l,t}^{-}$

$0\≤{\mathrm{\ΔP}}_{l,t}^{+}\≤{\mathrm{\ΔP}}_{l,\max }^{+}$

$0\≤{\mathrm{\ΔP}}_{l,t}^{-}\≤{\mathrm{\ΔP}}_{l,\max }^{-}$

其中，l∈L,k＝1,2,…,K，P_l,t为线路l在时段t潮流，P_l,max和P_l,min分别为线路l的正传输功率极限和反传输功率极限，G_l-i为机组i对线路l的功率转移分布因子，G_l-k为电网中节点k对线路l的功率转移分布因子，d_k,t为电网中节点k在t时刻的节点负荷预测值，K为所有节点集合，L为所有线路集合；和分别为允许线路l反向的最大越限量和正向的最大越限量，

所述断面传输极限约束为：

$S_{s,t}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{s-i}{P}_{i,t}-\underset{k=1}{\overset{K}{\mathrm{\Σ}}}{G}_{s-k}{d}_{k,t}$

$S_{s,t}-{\mathrm{\ΔS}}_{s,t}^{+}\≤S_{s,\max }$

$S_{s,\min }\≤S_{s,t}+{\mathrm{\ΔS}}_{s,t}^{-}$

$0\≤{\mathrm{\ΔS}}_{s,t}^{+}\≤{\mathrm{\ΔS}}_{s,\max }^{+}$

$0\≤{\mathrm{\ΔS}}_{s,t}^{-}\≤{\mathrm{\ΔS}}_{s,\max }^{-}$

其中，s∈S，k＝1,2,…,K，G_s-i为机组i对断面s的功率转移分布因子，G_s-k为电网中节点k对断面s的功率转移分布因子，S_s,t为断面s在时段t的潮流，和分别为允许的最大越限量。

6.如权利要求5所述的电网潮流阻塞调度方法，其特征在于，在省内大机组故障时，根据所述电网侧约束和优化目标求解所述电网潮流阻塞模型以获得目标函数的最小值。

7.如权利要求5所述的电网潮流阻塞调度方法，其特征在于，在区外直流大功率失却时，即区外直流大功率为零时，所述各省负荷平衡约束为：

$\mathrm{\Σ}{P}_{i,t}+e_{a,t}^{+}-e_{a,t}^{-}=D_{a,t}+{\mathrm{\ΔD}}_{a,t},$

其中，ΔD_a,t为分摊到省市的口子变化量。

8.一种电网潮流裕度评估方法，其特征在于，包括以下步骤：

根据断面在当前时刻的潮流和断面的正向最大越限量获得断面的目标系数，根据线路在当前时刻的潮流和线路的正向最大越限量获得线路的目标系数，根据断面的目标系数、机组对断面的功率转移分布因子、机组在当前时刻的有功出力、线路的目标系数和机组对该线路的功率转移分布因子获得潮流阻塞最严重情况下的目标函数以及潮流阻塞最佳情况下的目标函数；

根据机组出力上下限约束和爬坡约束获得机组侧约束，根据机组有功出力的上限、下限及机组在当前时刻的有功出力获得所述机组出力上下限约束，根据机组一时间间隔内所能升降的最大出力的上限和下限获得所述机组爬坡约束；

根据机组在当前时刻的有功出力、该省在当前时刻向其他省份额外购电交易量、其他省份在当前时刻向该省额外购电交易量和该省在当前时刻的等效负荷预测值获得各省负荷平衡约束；

以所述机组侧约束、各省负荷平衡约束为约束条件，根据所述约束条件和潮流阻塞最严重情况下的目标函数获得第一评估模型，根据所述约束条件和潮流阻塞最佳情况下的目标函数获得第二评估模型；

在约束条件限制下，求解所述第一评估模型以获得潮流阻塞最严重情况下的目标函数的最小值，求解所述第二评估模型以获得潮流阻塞最佳情况下的目标函数的最小值。

9.如权利要求8所述的电网潮流裕度评估方法，其特征在于，所述潮流阻塞最严重情况下的目标函数为：

$-\underset{t=1}{\overset{T}{\mathrm{\Σ}}}(C_s\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{s-i}{P}_{i,t}+C_l^{/}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{l-i}{P}_{i,t}),s\∈S,l\∈L,$

第一优化目标为所述潮流阻塞最严重情况下的目标函数的最小值，即：

$\min \{-\underset{t=1}{\overset{T}{\mathrm{\Σ}}}(C_s\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{s-i}{P}_{i,t}+C_l^{/}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{l-i}{P}_{i,t})\},s\∈S,l\∈L$

所述潮流阻塞最佳情况下的目标函数为：

$\underset{t=1}{\overset{T}{\mathrm{\Σ}}}(C_s\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{s-i}{P}_{i,t}+C_l^{/}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{l-i}{P}_{i,t}),s\∈S,l\∈L,$

第二优化目标为所述潮流阻塞最佳情况下的目标函数的最小值，即：

$\min \left\{\underset{t=1}{\overset{T}{\mathrm{\Σ}}}(C_s\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{s-i}{P}_{i,t}+C_l^{/}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{G}_{l-i}{P}_{i,t})\right\},s\∈S,l\∈L$

其中，i＝1,2,…,N，i∈A_a,t＝1,…,T，l∈L,k＝1,2,…,K，A_a为a省机组集合，T为决策周期内的总时段数，C_s和分别为断面s和线路l的目标系数，S_s,t为所述断面在当前时刻的潮流，S_s,max为所述断面的正向最大越限量，当前时刻为t时刻，S为所有断面集合，L为所有线路集合，P_i,t为机组i在t时刻的有功出力，G_l-i为机组i对线路l的功率转移分布因子，G_s-i为机组i对断面s的功率转移分布因子。

10.如权利要求9所述的电网潮流裕度评估方法，其特征在于，所述机组出力上下限约束为：

P_i,min≤P_i,t≤P_i,max，

所述机组爬坡约束为：

$\mathrm{\Δ}\stackrel{\‾}{P_i}\≤P_{i,t}-P_{i,t-1}\≤\mathrm{\Δ}\underset{\‾}{P_i},$

所述各省负荷平衡约束为：

$\mathrm{\Σ}{P}_{i,t}+e_{a,t}^{+}-e_{a,t}^{-}=D_{a,t},$

其中，i＝1,2,...,N,t＝1,...,T，P_i,t和P_i,t-1分别为机组i在t时刻及t-1时刻的有功出力；P_i,min和P_i,max分别为机组i的出力下限和上限，和表示机组i在t时刻和t-1时刻这一个时间间隔内所能升降的最大出力范围，D_a,t为a省时段t的等效负荷预测值。