CN111598725A - 一种基于网损的节点边际电价计算方法及装置 - Google Patents

Abstract

本申请公开了一种基于网损的节点边际电价计算方法及装置，方法包括：通过预置权重优化预置经济调度模型；构建优化后的预置经济调度模型的第一约束条件，第一约束条件包括：节点功率平衡约束、发电机组出力约束、直流潮流关系约束和网络传输能力约束；建立输电线路潮流与网损之间的第二约束条件，第二约束条件满足二次项约束；将优化后的预置经济调度模型、第一约束条件和第二约束条件转化为模型拉格朗日函数；根据模型拉格朗日函数和原始经济调度潮流计算节点边际电价。本申请解决了现有节点边际电价机制未考虑系统网损的影响，使得降低系统购电成本的同时，系统网损增大的技术问题。

Description

一种基于网损的节点边际电价计算方法及装置

技术领域

本申请涉及电力市场和调度运行技术领域，尤其涉及一种基于网损的节点边际电价计算方法及装置。

背景技术

节点边际电价是电力现货市场中重要的价格机制。与统一出清模式相比，节点边际电价充分考虑了电网传输能力导致的在空间上不同节点的电价差异，能够反映出当节点负荷波动时对全网购电成本造成的影响，因此，节点边际电价在电力市场得到了广泛的应用。但是，传统的节点边际电价并未考虑到系统的网损差异，容易造成降低系统购电成本的同时，系统网损增大的问题。

发明内容

本申请提供了一种基于网损的节点边际电价计算方法及装置，用于解决现有节点边际电价机制未考虑系统网损的影响，使得降低系统购电成本的同时，系统网损增大的技术问题。

有鉴于此，本申请第一方面提供了一种基于网损的节点边际电价计算方法，包括：

通过预置权重优化预置经济调度模型；

构建优化后的所述预置经济调度模型的第一约束条件，所述第一约束条件包括：节点功率平衡约束、发电机组出力约束、直流潮流关系约束和网络传输能力约束；

建立输电线路潮流与网损之间的第二约束条件，所述第二约束条件满足二次项约束；

将优化后的所述预置经济调度模型、所述第一约束条件和所述第二约束条件转化为模型拉格朗日函数，所述模型拉格朗日函数满足极值最优条件；

根据所述模型拉格朗日函数和原始经济调度潮流计算节点边际电价。

优选地，所述通过预置权重优化预置经济调度模型，包括:

根据预置优化公式优化购电成本和所述系统网损对所述预置经济调度模型的影响，所述预置优化公式为：

其中，NG、NL分别为系统中发电机组台数和输电线路条数，

为发电机组g的电能报价，

为发电机组g的中标电力，P_l ^Loss为输电线路l的传输功率损耗。α₁、α₂分别为购电成本和系统网损的所述预置权重。

优选地，所述节点功率平衡约束为：任一节点的注入功率和流出功率相等，所述注入功率包括发电机组有功功率和输电线路有功潮流，所述流出功率包括负荷功率。

优选地，所述发电机组出力约束为：每个发电机组的出力在预置出力上限值和预置出力下限值的范围内。

优选地，所述直流潮流关系约束为：线路有功潮流和节点相角之间满足约束公式，所述约束公式为：

其中，

分别为线路从节点i流向节点j的有功潮流和节点j流向节点i的有功潮流，θ_i、θ_j分别为线路两端节点i和节点j的电压相角，x_l(i,j)为输电线路l的阻抗值。

优选地，所述网络传输能力约束为：网络传输功率在输电线路的预置功率上限值和预置功率下限值范围内。

优选地，所述建立输电线路潮流与网损之间的第二约束条件，包括：

通过网损公式描述根据所述输电线路潮流与网损建立的第二约束条件，所述网损公式为：

相应地，二次项约束为：

P_l ^PFloss≥g_l(i,j)(θ_i-θ_j)²；

其中，P_l ^PFloss为输电线路l的损耗，g_l(i,j)为该线路的电阻。

优选地，所述将优化后的所述预置经济调度模型、所述第一约束条件和所述第二约束条件转化为模型拉格朗日函数，包括：

将优化后的所述预置经济调度模型、所述第一约束条件和所述第二约束条件抽象为数学模型；

将所述数学模型转化为所述模型拉格朗日函数。

本申请第二方面提供了一种基于网损的节点边际电价计算装置，包括：

第一构建模块，用于构建优化后的所述预置经济调度模型的第一约束条件，所述第一约束条件包括：节点功率平衡约束、发电机组出力约束、直流潮流关系约束和网络传输能力约束；

第二构建模块，用于建立输电线路潮流与网损之间的第二约束条件，所述第二约束条件满足二次项约束；

转化模块，用于将优化后的所述预置经济调度模型、所述第一约束条件和所述第二约束条件转化为模型拉格朗日函数，所述模型拉格朗日函数满足极值最优条件；

计算模块，用于根据所述模型拉格朗日函数和原始经济调度潮流计算节点边际电价。

所述第二构建模块具体用于：

通过网损公式描述根据所述输电线路潮流与网损建立的第二约束条件，所述网损公式为：

相应地，二次项约束为：

P_l ^PFloss≥g_l(i,j)(θ_i-θ_j)²；

其中，P_l ^PFloss为输电线路l的损耗，g_l(i,j)为该线路的电阻。

从以上技术方案可以看出，本申请实施例具有以下优点：

本申请中，提供了一种基于网损的节点边际电价计算方法，包括：通过预置权重优化预置经济调度模型；构建优化后的预置经济调度模型的第一约束条件，第一约束条件包括：节点功率平衡约束、发电机组出力约束、直流潮流关系约束和网络传输能力约束；建立输电线路潮流与网损之间的第二约束条件，第二约束条件满足二次项约束；将优化后的预置经济调度模型、第一约束条件和第二约束条件转化为模型拉格朗日函数，模型拉格朗日函数满足极值最优条件；根据模型拉格朗日函数和原始经济调度潮流计算节点边际电价。

本申请提供的基于网损的节点边际电价计算方法中，在经济调度模型中加入节点功率平衡约束、发电机组出力约束、直流潮流关系约束和网络传输能力约束，能够得到更加符合实际的经济调度模型，在预置的经济调度模型中加入网损的约束，再通过拉格朗日函数辅助进行模型优化，便于优化计算，通过优化计算得到的节点边际电价更能反映网损对其影响力的大小，便于操作人员根据实际的输电线路网损调整节点边际电价机制。因此，本申请提供的基于网损的几点边际电价计算方法解决了现有节点边际电价机制未考虑输电线路网损的影响，使得降低系统购电成本的同时，系统网损增大的技术问题。

附图说明

图1为本申请实施例提供的一种基于网损的节点边际电价计算方法的一个流程示意图；

图2为本申请实施例提供的一种基于网损的节点边际电价计算方法的另一个流程示意图；

图3为本申请实施例提供的一种基于网损的节点边际电价计算装置的一个结构示意图。

具体实施方式

为了使本技术领域的人员更好地理解本申请方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。

为了便于理解，请参阅图1，本申请提供的一种基于网损的节点边际电价计算方法的实施例一，包括：

步骤101、通过预置权重优化预置经济调度模型。

需要说明的是，电力系统经济调度是指在满足安全和电能质量的前提下，合理利用能源和设备，以最低的发电成本或燃料费用保证对用户可靠地供电的一种调度方法，而经济调度模型就是考量多方因素，建立合理的调度模型，提高对用户供电的高效性的同时也能保证供电商的效益。预置权重是人为设置的，主要用于优化经济调度模型中的购电成本和系统网损，一般购电成本和系统网损都是越低越好，但是不可能同时都达到最低，因此，需要通过人为设置的权重平衡二者的影响。

步骤102、构建优化后的预置经济调度模型的第一约束条件，第一约束条件包括：节点功率平衡约束、发电机组出力约束、直流潮流关系约束和网络传输能力约束。

需要说明的是，影响经济调度模型的优化的因素有很多，准确的选取这些影响因素用于约束经济调度模型的优化更加有利于电力调度市场的发展，更加迎合多方的个性化需求。经济调度模型的约束首先应当考虑的是电力系统的内部约束，例如节点功率平衡约束、发电机组处理约束和直流潮流关系约束等，这些属于根源处的约束，是衡量一个电力系统的供电效能的重要因素。而供电的传输过程也是关键环节，传输硬件自身的传输性能决定了电力传输的上限和下限，只有在这个范围内，才能完成正常的电力传输。

步骤103、建立输电线路潮流与网损之间的第二约束条件，第二约束条件满足二次项约束。

需要说明的是，现有的经济调度模型没有考虑到输电线路网损的差异对经济调度的影响，容易在降低购电费用的同时造成系统网损增大；系统的网损主要与输电线路潮流有关，因此可以根据输电线路潮流与网损之间的关系建立第二约束条件，与此同时，在这个约束条件的基础上增加二次项约束，能更加准确的表述系统网损的约束对经济调度模型的影响。

步骤104、将优化后的预置经济调度模型、第一约束条件和第二约束条件转化为模型拉格朗日函数，模型拉格朗日函数满足极值最优条件。

需要说明的是，获得了优化的预置经济调度模型，建立了第一约束条件和第二约束条件，将模型和约束条件联合起来就可以计算得到最后的约束结果，联合后的抽象数学模型计算较为复杂，因此，需要将联合后数学模型转化为模型拉格朗日函数，模型拉格朗日函数是数学优化问题计算分析中较为经典的辅助函数，可以将含有约束条件的优化问题转化为无约束条件的规划问题，所以，这个转化过程是一种优化计算的过程。

步骤105、根据模型拉格朗日函数和原始经济调度潮流计算节点边际电价。

需要说明的是，原始经济调度潮流可以根据本实施例中的预置经济调度模型计算得到，只需要将建立的第二约束条件去除，联立预置经济调度模型和第一约束条件就可以得到经典的经济调度模型，优化计算就能得到原始经济调度潮流。节点边际电价定义为当节点的负荷发生单位变化时，系统购电成本随之而变。

本实施例提供的基于网损的节点边际电价计算方法中，在经济调度模型中加入节点功率平衡约束、发电机组出力约束、直流潮流关系约束和网络传输能力约束，能够得到更加符合实际的经济调度模型，在预置的经济调度模型中加入网损的约束，再通过拉格朗日函数辅助进行模型优化，便于优化计算，通过优化计算得到的节点边际电价更能反映输电线路网损对其影响力的大小，便于操作人员根据实际的输电线路网损调整节点边际电价机制。因此，本实施例提供的基于网损的几点边际电价计算方法解决了现有节点边际电价机制未考虑输电线路网损的影响，使得降低系统购电成本时造成系统网损增大的技术问题。

为了便于理解，请参阅图2，本申请实施例中提供了一种基于网损的节点边际电价计算方法的实施例二，包括：

步骤201、根据预置优化公式优化购电成本和系统网损对预置经济调度模型的影响。

需要说明的是，电力系统经济调度是指在满足安全和电能质量的前提下，合理利用能源和设备，以最低的发电成本或燃料费用保证对用户可靠地供电的一种调度方法，而经济调度模型就是考量多方因素，建立合理的调度模型，提高对用户供电的高效性的同时也能保证供电商的效益。预置权重是人为设置的，主要用于优化经济调度模型中的购电成本和系统网损，一般购电成本和系统网损都是越低越好，但是系统中的不同输电线路的网损是存在差异的，倘若统一购电成本，显然是不能满足供电系统差异化需求的，也不是最合理的机制，因此，需要通过人为设置的权重平衡二者的影响。

其中，预置优化公式为：

其中，NG、NL分别为系统中发电机组台数和输电线路条数，

为发电机组g的电能报价，

为发电机组g的中标电力，P_l ^Loss为输电线路l的传输功率损耗。α₁、α₂分别为购电成本和系统网损的预置权重。

步骤202、构建优化后的预置经济调度模型的第一约束条件。

其中，第一约束条件包括：节点功率平衡约束、发电机组出力约束、直流潮流关系约束和网络传输能力约束。

需要说明的是，节点功率平衡约束为：任一节点的注入功率和流出功率相等，注入功率包括发电机组有功功率和输电线路有功潮流，流出功率包括负荷功率。具体的可以用节点功率平衡公式描述，节点功率平衡公式如下：

其中，g∈ψ₁(i)是与节点i相连接的发电机组的集合，j∈ψ₂(i)是以节点i为终点的线路集合，

为发电机组g的有功功率，

为输电线路l有功潮流，

为节点i的负荷功率。

需要说明的是，发电机组出力约束为：每个发电机组的出力在预置出力上限值和预置出力下限值的范围内。具体可以通过出力功率约束如下：

其中，

和

分别为发电机组g的出力下限功率和出力上限功率。

需要说明的是，直流潮流关系约束为：线路有功潮流和节点相角之间满足约束公式，约束公式为：

其中，

分别为线路从节点i流向节点j的有功潮流和节点j流向节点i的有功潮流，θ_i、θ_j分别为线路两端节点i和节点j的电压相角，x_l(i,j)为输电线路l的阻抗值。

需要说明的是，网络传输能力约束为：网络传输功率在输电线路的预置功率上限值和预置功率下限值范围内。具体的网络传输功率的约束表述如下：

其中，P_l ^PFow,min、P_l ^PFow,max分别为输电线路l的预置功率上限值和预置功率下限值。

需要说明的是，需要说明的是，影响经济调度模型的优化的因素有很多，准确的选取这些影响因素用于约束经济调度模型的优化更加有利于电力调度市场的发展，更加迎合多方的个性化需求。经济调度模型的约束首先应当考虑的是电力系统的内部约束，例如节点功率平衡约束、发电机组处理约束和直流潮流关系约束等，这些属于根源处的约束，是衡量一个电力系统的供电效能的重要因素。而供电的传输过程也是关键环节，传输硬件自身的传输性能决定了电力传输的上限和下限，只有在这个范围内，才能完成正常的电力传输。

步骤203、通过网损公式描述根据输电线路潮流与网损建立的第二约束条件。

需要说明的是，现有的经济调度模型没有考虑到输电线路网损的差异对经济调度的影响，容易在降低购电费用的同时，系统网损增大；系统的网损主要与输电线路潮流有关，因此可以根据输电线路潮流与网损之间的关系建立第二约束条件，与此同时，在这个约束条件的基础上增加二次项约束，能更加准确的表述输电线路潮流网损的约束对经济调度模型的影响。网损公式为：

相应地，二次项约束为：

P_l ^PFloss≥g_l(i,j)(θ_i-θ_j)²；

其中，P_l ^PFloss为输电线路l的损耗，g_l(i,j)为该线路的电阻。

步骤204、将优化后的预置经济调度模型、第一约束条件和第二约束条件抽象为数学模型。

需要说明的是，将上述优化的模型和所有的约束条件联立得到基于网损的经济调度模型，然后将其抽象为数学模型的形式，具体如下：

min F＝a^Gx；

其中，x为决策变量向量，a^G为目标函数一次项系数，

分别为一次等式约束中的一次项和常数项，

分别为一次不等式约束中的一次项和常数项，

分别为二次不等式约束中的二次项、一次项和常数项。则预置优化公式

对应上述抽象数学模型中的目标项minf＝a^Gx；线性等式约束

对应节点功率平衡约束

直流潮流关系约束

网损公式

线性不等式约束

对应发电机组出力约束

网络传输能力约束

二次不等式约束

对应网损二次项约束P_l ^PFloss≥g_l(i,j)(θ_i-θ_j)²。

需要说明的是，这里的目标函数即为需要优化的经济化调度模型，决策变量向量即为约束条件。决策变量向量包括发电机组出力

线路有功潮流

线路损耗P_l ^PFloss和节点电压相角θ_i。优化的数学模型中虽然包括二次项，但是在本实施例中，该二次项可转化为凸二次项约束。

步骤205、将数学模型转化为模型拉格朗日函数。

其中，模型拉格朗日函数满足极值最优条件。

需要说明的是，联合后的抽象数学模型计算较为复杂，因此，需要将联合后数学模型转化为模型拉格朗日函数，模型拉格朗日函数是数学优化问题计算分析中较为经典的辅助函数，可以将含有约束条件的优化问题转化为无约束条件的规划问题，所以，这个转化过程是一种优化计算的过程。具体的可以得到模型的拉格朗日函数如下：

其中，β₁、β₂、β₃分别是线性等式约束、线性不等式约束、凸二次不等式约束对应的拉格朗日乘子。

模型拉格朗日函数满足机制最优条件，是指根据数学优化理论，上述凸二次约束的经济调度问题的拉格朗日函数满足最优性条件：

也就是说，该凸二次约束存在最大极值点。

步骤206、根据模型拉格朗日函数和原始经济调度潮流计算节点边际电价。

需要说明的是，原始经济调度潮流可以根据本实施例中的预置经济调度模型计算得到，只需要将建立的第二约束条件去除，联立预置经济调度模型和第一约束条件就可以得到经典的经济调度模型，优化计算就能得到原始经济调度潮流。节点边际电价定义为当节点的负荷发生单位变化时，系统购电成本随之而变。具体的节点边际电价公式表示为：

其中，Pri_i为节点i的节点边际电价，由于考虑网损影响，因此相比于未考虑网损影响时增加了凸二次不等式约束对应的乘子β₃对应部分，其显性表达式则需要根据具体的计算案例表示。

为了便于理解，请参阅图3，本申请中还提供了一种基于网损的节点边际电价计算装置的实施例，包括：

优化模块301，用于通过预置权重优化预置经济调度模型；

第一构建模块302，用于构建优化后的预置经济调度模型的第一约束条件，第一约束条件包括：节点功率平衡约束、发电机组出力约束、直流潮流关系约束和网络传输能力约束；

第二构建模块303，用于建立输电线路潮流与网损之间的第二约束条件，第二约束条件满足二次项约束；

转化模块304，用于将优化后的预置经济调度模型、第一约束条件和第二约束条件转化为模型拉格朗日函数，模型拉格朗日函数满足极值最优条件；

计算模块305，用于根据模型拉格朗日函数和原始经济调度潮流计算节点边际电价。

其中，第二构建模块具体用于：

通过网损公式描述根据输电线路潮流与网损建立的第二约束条件，网损公式为：

相应地，二次项约束为：

P_l ^PFloss≥g_l(i,j)(θ_i-θ_j)²；

其中，P_l ^PFloss为输电线路l的损耗，g_l(i,j)为该线路的电阻。

在本申请所提供的几个实施例中，应该理解到，所揭露的装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。

所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。

另外，在本申请各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。

所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以通过一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(英文全称：Read-OnlyMemory，英文缩写：ROM)、随机存取存储器(英文全称：Random Access Memory，英文缩写：RAM)、磁碟或者光盘等各种可以存储程序代码的介质。

以上所述，以上实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围。

Claims (10)

1.一种基于网损的节点边际电价计算方法，其特征在于，包括：

通过预置权重优化预置经济调度模型；

构建优化后的所述预置经济调度模型的第一约束条件，所述第一约束条件包括：节点功率平衡约束、发电机组出力约束、直流潮流关系约束和网络传输能力约束；

建立输电线路潮流与网损之间的第二约束条件，所述第二约束条件满足二次项约束；

将优化后的所述预置经济调度模型、所述第一约束条件和所述第二约束条件转化为模型拉格朗日函数，所述模型拉格朗日函数满足极值最优条件；

根据所述模型拉格朗日函数和原始经济调度潮流计算节点边际电价。

2.根据权利要求1所述的基于网损的节点边际电价计算方法，其特征在于，所述通过预置权重优化预置经济调度模型，包括:

根据预置优化公式优化购电成本和系统网损对所述预置经济调度模型的影响，所述预置优化公式为：

其中，NG、NL分别为系统中发电机组台数和输电线路条数，

为发电机组g的电能报价，

为发电机组g的中标电力，P_l ^Loss为输电线路l的传输功率损耗，α₁、α₂分别为购电成本和系统网损的所述预置权重。

3.根据权利要求1所述的基于网损的节点边际电价计算方法，其特征在于，所述节点功率平衡约束为：任一节点的注入功率和流出功率相等，所述注入功率包括发电机组有功功率和输电线路有功潮流，所述流出功率包括负荷功率。

4.根据权利要求1所述的基于网损的节点边际电价计算方法，其特征在于，所述发电机组出力约束为：每个发电机组的出力在预置出力上限值和预置出力下限值的范围内。

5.根据权利要求1所述的基于网损的节点边际电价计算方法，其特征在于，所述直流潮流关系约束为：线路有功潮流和节点相角之间满足约束公式，所述约束公式为：

其中，

分别为线路从节点i流向节点j的有功潮流和节点j流向节点i的有功潮流，θ_i、θ_j分别为线路两端节点i和节点j的电压相角，x_l(i,j)为输电线路l的阻抗值。

6.根据权利要求1所述的基于网损的节点边际电价计算方法，其特征在于，所述网络传输能力约束为：网络传输功率在输电线路的预置功率上限值和预置功率下限值范围内。

7.根据权利要求5所述的基于网损的节点边际电价计算方法，其特征在于，所述建立输电线路潮流与网损之间的第二约束条件，包括：

通过网损公式描述根据所述输电线路潮流与网损建立的第二约束条件，所述网损公式为：

相应地，二次项约束为：

P_l ^PFloss≥g_l(i,j)(θ_i-θ_j)²；

其中，P_l ^PFloss为输电线路l的损耗，g_l(i,j)为该线路的电阻。

8.根据权利要求1中所述的基于网损的节点边际电价计算方法，其特征在于，所述将优化后的所述预置经济调度模型、所述第一约束条件和所述第二约束条件转化为模型拉格朗日函数，包括：

将优化后的所述预置经济调度模型、所述第一约束条件和所述第二约束条件抽象为数学模型；

将所述数学模型转化为所述模型拉格朗日函数。

9.一种基于网损的节点边际电价计算装置，其特征在于，包括：

优化模块，用于通过预置权重优化预置经济调度模型；

第一构建模块，用于构建优化后的所述预置经济调度模型的第一约束条件，所述第一约束条件包括：节点功率平衡约束、发电机组出力约束、直流潮流关系约束和网络传输能力约束；

第二构建模块，用于建立输电线路潮流与网损之间的第二约束条件，所述第二约束条件满足二次项约束；

转化模块，用于将优化后的所述预置经济调度模型、所述第一约束条件和所述第二约束条件转化为模型拉格朗日函数，所述模型拉格朗日函数满足极值最优条件；

计算模块，用于根据所述模型拉格朗日函数和原始经济调度潮流计算节点边际电价。

10.根据权利要求9中所述的基于网损的节点边际电价计算装置，其特征在于，所述第二构建模块具体用于：

通过网损公式描述根据所述输电线路潮流与网损建立的第二约束条件，所述网损公式为：

相应地，二次项约束为：

P_l ^PFloss≥g_l(i,j)(θ_i-θ_j)²；

其中，P_l ^PFloss为输电线路l的损耗，g_l(i,j)为该线路的电阻。

Images