CN104092213A - 一种基于优化方法的不确定潮流支路功率分析方法 - Google Patents

Abstract

本发明提供一种基于优化方法的不确定潮流支路功率分析方法，其步骤包括：遍历整个网络，统计网络中节点数量与类型，并为节点编号；将每个节点的负荷波动表示为包含上限与下限信息的区间形式；分别建立基于潮流方程、节点参数及电力系统运行限制的约束条件；分别以每条支路上流过的有功功率及无功功率作为目标函数，结合已构建的约束条件，建立支路有功功率优化模型和支路无功功率优化模型；使用优化算法求解优化模型，得到各条支路上流过功率的波动范围。根据本发明所述的方法，可以完全规避保守性问题对支路功率波动范围的影响。此外，本发明所述的方法中，不同的优化模型具有相同的约束条件，可以使用并行计算技术，提高计算效率。

Description

一种基于优化方法的不确定潮流支路功率分析方法

技术领域

本发明属于电力系统稳态分析领域，涉及一种电力系统不确定潮流支路功率分析方法，更具体地，涉及一种基于优化方法的不确定潮流支路功率分析方法。

背景技术

不确定潮流是针对电力系统内负荷和发电机的出力不能精确知道，但知道其一定包含在某个给定区间范围内的情况，通过分析得到系统中各状态量及潮流分布的区间解。支路功率分析是潮流分析过程中的一个重要组成部分，是对系统中各条支路上流过的有功功率及无功功率的求解与分析，其是研究电力系统运行、规划以及安全性、可靠性的基础，属于电力系统稳态分析。

针对不确定潮流问题，前人提出了基于电流注入方程的区间潮流分析方法，该方法使用Krawczyk迭代算子求解区间非线性方程组，从而达到求解不确定潮流分布的目的。然而，受限于区间运算的保守性，使用传统的支路功率计算公式，应用Krawczyk迭代算子分析所得的节点电压幅值与相角的区间结果计算支路功率，结果往往过于保守，包含很多实际上不可能出现的电力系统运行状态，以至于所得支路功率的波动范围毫无参考价值。此外，Krawczyk算子迭代过程将花费大量的时间成本，降低不确定潮流支路功率的分析效率。

发明内容

技术问题：本发明的目的是提供一种能够有效求解不确定潮流支路功率波动范围的基于优化方法的不确定潮流支路功率分析方法。

技术方案：本发明的基于优化方法的不确定潮流支路功率分析方法，包括以下步骤：

1)遍历整个电力网络，统计PQ节点、PV节点、平衡节点3种节点，确定节点个数N并为节点编号；

2)将每个节点的负荷波动表示为包含上限与下限信息的区间形式；

3)分别建立基于潮流方程的等式约束、基于PQ节点给定参数的有功功率不等式约束和无功功率不等式约束、基于PV节点给定参数的有功功率不等式约束和电压幅值等式约束、基于平衡节点给定参数的电压实部等式约束和电压虚部等式约束，以及在电力网络包含系统状态量阈值时，还需建立基于电力系统运行限制的约束；

4)以每条支路上流过的有功功率表达式作为目标函数，结合所述步骤3)中建立的全部约束条件，构建得到支路有功功率优化模型；同时以每条支路上流过的无功功率表达式作为目标函数，结合所述步骤3)中建立的全部约束条件，构建得到支路无功功率优化模型；

5)求解所述步骤4)中得到的两种优化模型，得到各条支路上流过的功率的波动范围。

本发明的优选方案中，步骤1)中，节点编号由自然数1开始，直至节点个数N，在编号过程中无需考虑节点类型的影响。

本发明的步骤2)中，分别设置PQ节点有功功率的最大值和最小值、PQ节点无功功率的最大值和最小值、PV节点有功功率的最大值和最小值，从而将节点的负荷波动表征为包含上限与下限信息的区间形式。

本发明的优选方案中，步骤3)中：

所述基于潮流方程的等式约束为：

$P_i-e_i\underset{j=1}{\overset{N}{\mathrm{\Σ}}}(G_{\mathrm{ij}}{e}_j-B_{\mathrm{ij}}{f}_j)-f_i\underset{j=1}{\overset{N}{\mathrm{\Σ}}}(G_{\mathrm{ij}}{f}_j+B_{\mathrm{ij}}{e}_j)=0$

$Q_i-f_i\underset{j=1}{\overset{N}{\mathrm{\Σ}}}(G_{\mathrm{ij}}{e}_j-B_{\mathrm{ij}}{f}_i)+e_i\underset{j=1}{\overset{N}{\mathrm{\Σ}}}(G_{\mathrm{ij}}{f}_j+B_{\mathrm{ij}}{e}_j)=0$

其中，P_i、Q_i分别表示节点i处注入的有功功率与无功功率，e_i、f_i分别表示节点i处电压的实部与虚部，e_j、f_j分别表示节点j处电压的实部与虚部，G_ij、B_ij分别表示连接节点i与节点j的支路的电导与电纳，N表示系统中的节点个数；

所述基于PQ节点给定参数的有功功率不等式约束为：

$\underset{\‾}{P_{\mathrm{is}}}\≤P_i\≤\stackrel{\‾}{P_{\mathrm{is}}}$

所述基于PQ节点给定参数的无功功率不等式约束为：

$\underset{\‾}{Q_{\mathrm{is}}}\≤Q_i\≤\stackrel{\‾}{Q_{\mathrm{is}}}$

其中，P_i表示节点i处注入的有功功率，Q_i表示节点i处注入的无功功率，P _is 表示步骤2)中节点i处的有功功率的区间下限，表示步骤2)中节点i处的有功功率的区间上限，Q _is 表示步骤2)中节点i处的无功功率的区间下限，表示步骤2)中节点i处的无功功率的区间上限；

所述基于PV节点给定参数的有功功率不等式约束为：

$\underset{\‾}{P_{\mathrm{is}}}\≤P_i\≤\stackrel{\‾}{P_{\mathrm{is}}}$

所述基于PV节点给定参数的电压幅值等式约束为：

U_is ²＝e_i ²+f_i ²

其中，P_i表示节点i处注入的有功功率，P _is 表示步骤2)中节点i处的有功功率的区间下限，表示步骤2)中节点i处的有功功率的区间上限，U_is表示节点i处给定的电压幅值，e_i、f_i分别表示节点i处电压的实部与虚部；

所述基于平衡节点给定参数的电压实部等式约束为：

e_refs＝e_ref

所述基于平衡节点给定参数的电压虚部等式约束为：

f_refs＝f_ref

其中，e_ref、f_ref分别表示计算所得平衡节点处的电压的实部与虚部，e_refs、f_refs分别表示给定的平衡节点处的电压的实部与虚部。

所述基于电力系统运行限制的约束为：

f(x)≤0

其中，x为系统中的任一具有限制条件的状态量，f(x)是关于该状态量的函数。

本发明的优选方案中，步骤4)中构建得到的支路有功功率优化模型如下：

obj.   min(max)  P_ij

s.t.   h_k(x)＝0  k＝1,2,…,m

       g_t(x)≤0  t＝1,2,…,n

步骤4)中构建得到的支路无功功率优化模型如下：

obj.   min(max)  Q_ij

s.t.   h_k(x)＝0  k＝1,2,…,m

       g_t(x)≤0  t＝1,2,…,n

其中，obj.是优化模型中目标函数的标志符，s.t.是优化模型中约束条件的标志符，P_ij表示连接节点i与节点j的支路上流过的有功功率， $P_{\mathrm{ij}}=[R_{\mathrm{ij}}\·e_i\·(e_i-e_j)+X_{\mathrm{ij}}\·e_i\·(f_i-f_j)+R_{\mathrm{ij}}\·f_i\·(f_i-f_j)-X_{\mathrm{ij}}\·f_i\·(e_i-e_j)]/({R_{\mathrm{ij}}}^2+X_{\mathrm{ij}}^2),$ Q_ij表示连接节点i与节点j的支路上流过的无功功率， $Q_{\mathrm{ij}}=[R_{\mathrm{ij}}\·f_i\·(e_i-e_j)+X_{\mathrm{ij}}\·f_i\·(f_i-f_j)-R_{\mathrm{ij}}\·e_i\·(f_i-f_j)+X_{\mathrm{ij}}\·e_i\·(e_i-e_j)]/({R_{\mathrm{ij}}}^2+X_{\mathrm{ij}}^2)$ R_ij、X_ij分别表示连接节点i与节点j的支路上的串联电阻与电抗，B_ij表示连接节点i与节点j的支路的电纳，e_i、f_i分别表示节点i处电压的实部与虚部，e_j、f_j分别表示节点j处电压的实部与虚部，min表示以目标函数表达式的最小值作为优化目标，max表示以目标函数表达式的最大值作为优化目标，h_k(x)，k＝1,2,…,m，是定义在实数域内的关于系统状态量的函数，k为该函数的序号，m表示约束条件中等式约束的个数，g_t(x)，t＝1,2,…,n，是定义在实数域内的关于系统状态量的函数，t为该函数的序号，n表示约束条件中不等式约束的个数。

本发明的优选方案中，步骤5)的具体方法为：求解支路有功功率优化模型，将得到的各目标函数的最大值和最小值分别作为对应支路有功功率波动范围的上限与下限；求解支路无功功率优化模型，将得到的各目标函数的最大值和最小值分别作为对应支路无功功率波动范围的上限与下限。

有益效果：本发明与现有技术相比，具有以下优点：

前人提出的使用Krawczyk迭代算子的基于电流注入方程的区间潮流分析方法，所得的潮流结果往往过于保守，应用该方法的节点电压幅值及相角结果计算支路功率，所得结果的保守性更为严重，包含很多实际上不可能出现的电力系统运行状态，范围巨大以至于毫无参考价值。本发明提供的基于优化方法的不确定潮流支路功率分析方法，将传统的不确定潮流支路功率分析方法中的迭代求解过程转化为优化问题，将支路上流过的有功功率及无功功率作为优化模型的目标函数，通过求解优化模型，得到每条支路上流过功率的上限与下限，进而得到支路功率的波动范围。在优化问题中，目标函数的极大值与极小值都是客观存在的。因此，使用本发明所述方法，所得每条支路上的功率的上限与下限都是客观存在的，即，该方法完全规避了保守性问题对支路功率波动范围的影响。此外，前人提出的迭代方法在迭代过程中将花费大量的时间成本，降低了不确定潮流支路功率的分析效率。在本发明提供的基于优化方法的不确定潮流支路功率分析方法中，由于求解不同支路上流过的功率时所用优化模型的约束条件是相同的，因此可以使用并行计算技术，同时进行多条支路上功率的求解，从而达到提高分析效率的目的。

附图说明

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

图2是编号后的IEEE-14节点系统结构图。

具体实施方式

下面结合实施例和说明书附图，对本发明的技术方案进行具体介绍。

图1是本发明的方法流程示意图，介绍了本发明方法的基本步骤。图2是编号后的IEEE-14节点系统，给定了网络中平衡节点的电压幅值及相角大小、PQ节点的负荷波动范围(在本例中设定为有功功率及无功功率均关于确定值波动±10％)、PV节点的电压幅值大小及有功功率波动范围(在本例中设定为有功功率关于确定值波动±10％)，以下以该系统为例说明本发明方法的具体实现。

1)遍历整个网络，统计PQ节点、PV节点、平衡节点3种节点，确定网络节点个数N并为节点编号。电力系统中存在3种类型的节点，其中，有功功率和无功功率已知的节点为PQ节点，有功功率和电压幅值已知的节点为PV节点，电压幅值和电压相角已知的节点为平衡节点。以图2所示系统为例，该网络共有14个节点且已从1至14编号，其中编号1的节点为平衡节点，编号4、5、7、9、10、11、12、13、14的节点为PQ节点，编号2、3、6、8的节点为PV节点。

2)将步骤1)所述网络中的负荷波动表示为包含上限与下限信息的区间形式。以图2所示系统中的节点2和节点9为例，其中节点2为PV节点，负荷额定有功功率标幺化后为0.217；节点9为PQ节点，负荷额定有功功率标幺化后为0.295，额定无功功率标幺化后为0.166。考虑±10％的负荷波动性，则节点2的负荷有功功率可以表征为区间[0.1953,0.2387]，节点9的负荷有功功率和无功功率可以分别表征为区间[0.2655,0.3245]和[0.1494,0.1826]。

3)针对已有系统网络，分别建立基于潮流方程的等式约束、基于PQ节点给定参数的有功功率不等式约束和无功功率不等式约束、基于PV节点给定参数的有功功率不等式约束和电压幅值等式约束、基于平衡节点给定参数的电压实部等式约束和电压虚部等式约束，以及在电力网络包含系统状态量阈值时，还需建立基于电力系统运行限制的约束。

基于潮流方程的等式约束为：

$P_i-e_i\underset{j=1}{\overset{N}{\mathrm{\Σ}}}(G_{\mathrm{ij}}{e}_j-B_{\mathrm{ij}}{f}_j)-f_i\underset{j=1}{\overset{N}{\mathrm{\Σ}}}(G_{\mathrm{ij}}{f}_j+B_{\mathrm{ij}}{e}_j)=0$

$Q_i-f_i\underset{j=1}{\overset{N}{\mathrm{\Σ}}}(G_{\mathrm{ij}}{e}_j-B_{\mathrm{ij}}{f}_i)+e_i\underset{j=1}{\overset{N}{\mathrm{\Σ}}}(G_{\mathrm{ij}}{f}_j+B_{\mathrm{ij}}{e}_j)=0$

其中，P_i、Q_i分别表示节点i处注入的有功功率与无功功率，e_i、f_i分别表示节点i处电压的实部与虚部，e_j、f_j分别表示节点j处电压的实部与虚部，G_ij、B_ij分别表示连接节点i与节点j的支路的电导与电纳，N表示系统中的节点个数。每个节点均存在有功功率潮流方程约束与无功功率潮流方程约束，因此，该类约束条件共有2N个。图2所示的系统中共有14个节点，该类约束条件共有28个。以节点3为例，与其相连的有节点2与节点4，关于节点3的潮流方程约束可以表示为：

P₃-e₃[(G₂₃e₂-B₂₃f₂)+(G₃₄e₄-B₃₄f₄)]-f₃[(G₂₃f₂+B₂₃e₂)+(G₃₄f₄+B₃₄e₄)]＝0

Q₃-f₃[(G₂₃e₂-B₂₃f₂)+(G₃₄e₄-B₃₄f₄)]+e₃[(G₂₃f₂+B₂₃e₂)+(G₃₄f₄+B₃₄e₄)]＝0

基于PQ节点给定参数的有功功率不等式约束为：

$\underset{\‾}{P_{\mathrm{is}}}\≤P_i\≤\stackrel{\‾}{P_{\mathrm{is}}}$

基于PQ节点给定参数的无功功率不等式约束为：

$\underset{\‾}{Q_{\mathrm{is}}}\≤Q_i\≤\stackrel{\‾}{Q_{\mathrm{is}}}$

其中，P_i表示节点i处注入的有功功率，Q_i表示节点i处注入的无功功率，P _is 表示步骤2)中节点i处的有功功率的区间下限，表示步骤2)中节点i处的有功功率的区间上限，Q _is 表示步骤2)中节点i处的无功功率的区间下限，表示步骤2)中节点i处的无功功率的区间上限。将系统中PQ节点的个数记为N_PQ，则基于PQ节点给定参数的约束条件共有2N_PQ个。图2所示的系统中共有9个PQ节点，该类约束条件共有18个。以节点9为例，其基于PQ节点给定参数的约束可以表示为：

$\underset{\‾}{P_{9s}}\≤P_9\≤\stackrel{\‾}{P_{9s}}$

$\underset{\‾}{Q_{9s}}\≤Q_9\≤\stackrel{\‾}{Q_{9s}}$

其中，节点9处的负荷有功功率与无功功率的下限和上限如步骤2)中所述。

基于PV节点给定参数的有功功率不等式约束为：

$\underset{\‾}{P_{\mathrm{is}}}\≤P_i\≤\stackrel{\‾}{P_{\mathrm{is}}}$

基于PV节点给定参数的电压幅值等式约束为：

U_is ²＝e_i ²+f_i ²

其中，P_i表示节点i处注入的有功功率，P _is 表示步骤2)中节点i处的有功功率的区间下限，表示步骤2)中节点i处的有功功率的区间上限，U_is表示节点i处给定的电压幅值，e_i、f_i分别表示节点i处电压的实部与虚部。将系统中PV节点的个数记为N_PV，则基于PV节点给定参数的约束条件共有2N_PV个。图2所示的系统中共有4个PV节点，该类约束条件共有8个。以节点2为例，其基于PV节点给定参数的约束可以表示为：

$\underset{\‾}{P_{2s}}\≤P_2\≤\stackrel{\‾}{P_{2s}}$

U_2s ²＝e₂ ²+f₂ ²

其中，节点2处的电压幅值为给定的确定量，负荷有功功率下限和上限如步骤2)中所述。

基于平衡节点给定参数的电压实部等式约束为：

e_refs＝e_ref

基于平衡节点给定参数的电压虚部等式约束为：

f_refs＝f_ref

其中，e_ref、f_ref分别表示计算所得平衡节点处的电压的实部与虚部，e_refs、f_refs分别表示给定的平衡节点处的电压的实部与虚部。通常的电力系统中只存在一个平衡节点，因此，基于平衡节点给定参数的约束条件只有2个。图2所示的系统中，节点

1为平衡节点，则基于平衡节点给定参数的约束为：

e_1s＝e₁

f_1s＝f₁

基于电力系统运行限制的约束为：

f(x)≤0

其中，x为系统中任一具有限制条件的状态量，f(x)是关于该状态量的函数。该类约束条件在设置时较为随意，不同的电力系统具有不同的运行限制，随着电力系统运行限制的改变，该类约束条件也会随之改变。假设在图2所示的系统中，电力系统运行要求节点电压满足：

0.8≤e_i≤1.2

-0.3≤f_i≤0.3

其中，e_i、f_i分别表示节点i处电压的实部与虚部。图2所示的系统中共有14个节点，该类基于电力系统运行限制的约束共有28个。

4)在步骤3)的基础上，构建得到支路有功功率优化模型和支路无功功率优化模型：

根据支路功率的计算公式

$S_{\mathrm{ij}}={\stackrel{\·}{U}}_i\·{(\frac{{\stackrel{\·}{U}}_i-{\stackrel{\·}{U}}_j}{R_{\mathrm{ij}}+{\mathrm{jX}}_{\mathrm{ij}}}+{\stackrel{\·}{U}}_i\·j\frac{1}{2}{B}_{\mathrm{ij}})}^{*}$

可以推导出支路有功功率的表达式

$P_{\mathrm{ij}}=[R_{\mathrm{ij}}\·e_i\·(e_i-e_j)+X_{\mathrm{ij}}\·e_i\·(f_i-f_j)+R_{\mathrm{ij}}\·f_i\·(f_i-f_j)-X_{\mathrm{ij}}\·f_i\·(e_i-e_j)]/({R_{\mathrm{ij}}}^2+X_{\mathrm{ij}}^2)$

以及支路无功功率的表达式

$\begin{array}{c}{Q}_{\mathrm{ij}}=[R_{\mathrm{ij}}\·f_i\·(e_i-e_j)+X_{\mathrm{ij}}\·f_i\·(f_i-f_j)-R_{\mathrm{ij}}\·e_i\·(f_i-f_j)+X_{\mathrm{ij}}\·e_i\·(e_i-e_j)]/({R_{\mathrm{ij}}}^2+X_{\mathrm{ij}}^2)\\ -\frac{1}{2}\·{f_i}^2\·B_{\mathrm{ij}}-\frac{1}{2}\·{e_i}^2\·B_{\mathrm{ij}}\end{array}$

其中，S_ij表示连接节点i与节点j的支路上流过的复功率，分别表示节点i处和节点j处的电压相量，R_ij、X_ij分别表示连接节点i与节点j的支路上的串联电阻与电抗，B_ij表示连接节点i与节点j的支路上的电纳，P_ij表示连接节点i与节点j的支路上流过的有功功率，Q_ij表示连接节点i与节点j的支路上流过的无功功率，e_i、f_i分别表示节点i处电压的实部与虚部，e_j、f_j分别表示节点j处电压的实部与虚部。

以每条支路上流过的有功功率表达式作为目标函数，结合步骤3)中所述的全部约束条件，构建得到的支路有功功率优化模型如下：

obj.   min(max)  P_ij

s.t.   h_k(x)＝0  k＝1,2,…,m

       g_t(x)≤0  t＝1,2,…,n

以每条支路上流过的无功功率表达式作为目标函数，结合步骤3)中所述的全部约束条件，构建得到的支路无功功率优化模型如下：

obj.   min(max)  Q_ij

s.t.   h_k(x)＝0  k＝1,2,…,m

       g_t(x)≤0  t＝1,2,…,n

其中，obj.是优化模型中目标函数的标志符，s.t.是优化模型中约束条件的标志符，P_ij表示连接节点i与节点j的支路上流过的有功功率，Q_ij表示连接节点i与节点j的支路上流过的无功功率，R_ij、X_ij分别表示连接节点i与节点j的支路上的串联电阻与电抗，B_ij表示连接节点i与节点j的支路的电纳，e_i、f_i分别表示节点i处电压的实部与虚部，e_j、f_j分别表示节点j处电压的实部与虚部，min表示以目标函数表达式的最小值作为优化目标，max表示以目标函数表达式的最大值作为优化目标，h_k(x)，k＝1,2,…,m，是定义在实数域内的关于系统状态量的函数，k为该函数的序号，m表示约束条件中等式约束的个数，g_t(x)，t＝1,2,…,n，是定义在实数域内的关于系统状态量的函数，t为该函数的序号，n表示约束条件中不等式约束的个数。图2所示系统中共有9个PQ节点、4个PV节点、1个平衡节点，且有28个基于电力系统运行限制的约束，因此，在支路有功功率优化模型和支路无功功率优化模型中m＝34，n＝50。以连接节点2与节点3的支路为例，该支路上的有功功率优化模型可以表示为：

obj.   min(max)  P₂₃

s.t.   h_k(x)＝0  k＝1,2,…,34

       g_t(x)≤0  t＝1,2,…,50

以连接节点2与节点3的支路为例，该支路上的无功功率优化模型可以表示为：

obj.   min(max)  Q₂₃

s.t.   h_k(x)＝0  k＝1,2,…,34

       g_t(x)≤0  t＝1,2,…,50

5)使用优化算法求解步骤4)中所得的优化模型，可以得到每条支路上流过的有功功率及无功功率波动范围，具体方法为：求解支路有功功率优化模型，将得到的各目标函数的最大值和最小值分别作为对应支路有功功率波动范围的上限与下限；求解支路无功功率优化模型，将得到的各目标函数的最大值和最小值分别作为对应支路无功功率波动范围的上限与下限。求解非线性优化模型的方法多种多样，也有诸多功能强大的优化软件可供使用，本发明不再赘述。仍以图2所示系统中连接节点2与节点3的支路为例，通过求解该支路的有功功率优化模型，可以得到该支路上从节点2流向节点3的有功功率的波动范围；通过求解该支路的无功功率优化模型，可以得到该支路上从节点2流向节点3的无功功率的波动范围。

以上仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干可以预期的改进和等同替换，这些对本发明权利要求进行改进和等同替换后的技术方案，均落入本发明的保护范围。

Claims (6)

1.一种基于优化方法的不确定潮流支路功率分析方法，其特征在于，该方法包括如下步骤：

1)遍历整个电力网络，统计PQ节点、PV节点、平衡节点3种节点，确定节点个数N并为节点编号；

2)将每个节点的负荷波动表示为包含上限与下限信息的区间形式；

3)分别建立基于潮流方程的等式约束、基于PQ节点给定参数的有功功率不等式约束和无功功率不等式约束、基于PV节点给定参数的有功功率不等式约束和电压幅值等式约束、基于平衡节点给定参数的电压实部等式约束和电压虚部等式约束，以及在电力网络包含系统状态量阈值时，还需建立基于电力系统运行限制的约束；

4)以每条支路上流过的有功功率表达式作为目标函数，结合所述步骤3)中建立的全部约束条件，构建得到支路有功功率优化模型；同时以每条支路上流过的无功功率表达式作为目标函数，结合所述步骤3)中建立的全部约束条件，构建得到支路无功功率优化模型；

5)求解所述步骤4)中得到的两种优化模型，得到各条支路上流过的功率的波动范围。

2.根据权利要求1所述的基于优化方法的不确定潮流支路功率分析方法，其特征在于，所述步骤1)中，节点编号由自然数1开始，直至节点个数N，在编号过程中无需考虑节点类型的影响。

3.根据权利要求1所述的基于优化方法的不确定潮流支路功率分析方法，其特征在于，所述步骤2)中，分别设置PQ节点有功功率的最大值和最小值、PQ节点无功功率的最大值和最小值、PV节点有功功率的最大值和最小值，从而将节点的负荷波动表征为包含上限与下限信息的区间形式。

4.根据权利要求1、2或3所述的基于优化方法的不确定潮流支路功率分析方法，其特征在于，所述步骤3)中：

所述基于潮流方程的等式约束为：

$P_i-e_i\underset{j=1}{\overset{N}{\mathrm{\Σ}}}(G_{\mathrm{ij}}{e}_j-B_{\mathrm{ij}}{f}_j)-f_i\underset{j=1}{\overset{N}{\mathrm{\Σ}}}(G_{\mathrm{ij}}{f}_j+B_{\mathrm{ij}}{e}_j)=0$

$Q_i-f_i\underset{j=1}{\overset{N}{\mathrm{\Σ}}}(G_{\mathrm{ij}}{e}_j-B_{\mathrm{ij}}{f}_i)+e_i\underset{j=1}{\overset{N}{\mathrm{\Σ}}}(G_{\mathrm{ij}}{f}_j+B_{\mathrm{ij}}{e}_j)=0$

其中，P_i、Q_i分别表示节点i处注入的有功功率与无功功率，e_i、f_i分别表示节点i处电压的实部与虚部，e_j、f_j分别表示节点j处电压的实部与虚部，G_ij、B_ij分别表示连接节点i与节点j的支路的电导与电纳，N表示系统中的节点个数；

所述基于PQ节点给定参数的有功功率不等式约束为：

$\underset{\‾}{P_{\mathrm{is}}}\≤P_i\≤\stackrel{\‾}{P_{\mathrm{is}}}$

所述基于PQ节点给定参数的无功功率不等式约束为：

$\underset{\‾}{Q_{\mathrm{is}}}\≤Q_i\≤\stackrel{\‾}{Q_{\mathrm{is}}}$

其中，P_i表示节点i处注入的有功功率，Q_i表示节点i处注入的无功功率，P _is 表示步骤2)中节点i处的有功功率的区间下限，表示步骤2)中节点i处的有功功率的区间上限，Q _is 表示步骤2)中节点i处的无功功率的区间下限，表示步骤2)中节点i处的无功功率的区间上限；

所述基于PV节点给定参数的有功功率不等式约束为：

$\underset{\‾}{P_{\mathrm{is}}}\≤P_i\≤\stackrel{\‾}{P_{\mathrm{is}}}$

所述基于PV节点给定参数的电压幅值等式约束为：

U_is ²＝e_i ²+f_i ²

其中，P_i表示节点i处注入的有功功率，P _is 表示步骤2)中节点i处的有功功率的区间下限，表示步骤2)中节点i处的有功功率的区间上限，U_is表示节点i处给定的电压幅值，e_i、f_i分别表示节点i处电压的实部与虚部；

所述基于平衡节点给定参数的电压实部等式约束为：

e_refs＝e_ref

所述基于平衡节点给定参数的电压虚部等式约束为：

f_refs＝f_ref

其中，e_ref、f_ref分别表示计算所得平衡节点处的电压的实部与虚部，e_refs、f_refs分别表示给定的平衡节点处的电压的实部与虚部。

所述基于电力系统运行限制的约束为：

f(x)≤0

其中，x为系统中的任一具有限制条件的状态量，f(x)是关于该状态量的函数。

5.根据权利要求1、2或3所述的基于优化方法的不确定潮流支路功率分析方法，其特征在于，所述步骤4)中构建得到的支路有功功率优化模型如下：

obj.   min(max)  P_ij

s.t.   h_k(x)＝0  k＝1,2,…,m

       g_t(x)≤0  t＝1,2,…,n

步骤4)中构建得到的支路无功功率优化模型如下：

obj.   min(max)  Q_ij

s.t.   h_k(x)＝0  k＝1,2,…,m

       g_t(x)≤0  t＝1,2,…,n

其中，obj.是优化模型中目标函数的标志符，s.t.是优化模型中约束条件的标志符，P_ij表示连接节点i与节点j的支路上流过的有功功率， $P_{\mathrm{ij}}=[R_{\mathrm{ij}}\·e_i\·(e_i-e_j)+X_{\mathrm{ij}}\·e_i\·(f_i-f_j)+R_{\mathrm{ij}}\·f_i\·(f_i-f_j)-X_{\mathrm{ij}}\·f_i\·(e_i-e_j)]/({R_{\mathrm{ij}}}^2+X_{\mathrm{ij}}^2),$ Q_ij表示连接节点i与节点j的支路上流过的无功功率， $Q_{\mathrm{ij}}=[R_{\mathrm{ij}}\·f_i\·(e_i-e_j)+X_{\mathrm{ij}}\·f_i\·(f_i-f_j)-R_{\mathrm{ij}}\·e_i\·(f_i-f_j)+X_{\mathrm{ij}}\·e_i\·(e_i-e_j)]/({R_{\mathrm{ij}}}^2+X_{\mathrm{ij}}^2)$ R_ij、X_ij分别表示连接节点i与节点j的支路上的串联电阻与电抗，B_ij表示连接节点i与节点j的支路的电纳，e_i、f_i分别表示节点i处电压的实部与虚部，e_j、f_j分别表示节点j处电压的实部与虚部，min表示以目标函数表达式的最小值作为优化目标，max表示以目标函数表达式的最大值作为优化目标，h_k(x)，k＝1,2,…,m，是定义在实数域内的关于系统状态量的函数，k为该函数的序号，m表示约束条件中等式约束的个数，g_t(x)，t＝1,2,…,n，是定义在实数域内的关于系统状态量的函数，t为该函数的序号，n表示约束条件中不等式约束的个数。

6.根据权利要求1、2或3所述的基于优化方法的不确定潮流支路功率分析方法，其特征在于，所述步骤5)的具体方法为：求解支路有功功率优化模型，将得到的各目标函数的最大值和最小值分别作为对应支路有功功率波动范围的上限与下限；求解支路无功功率优化模型，将得到的各目标函数的最大值和最小值分别作为对应支路无功功率波动范围的上限与下限。