CN109193628A - 一种基于一致性的综合能源系统能量管理方法 - Google Patents

Abstract

本发明提出一种基于一致性的综合能源系统能量管理方法，热力网、电力网，能量供应侧和负荷需求侧的参与主体，都参与到能量市场中来追求自身的利益最大化，构成全社会福利最大化模型。CHP机组运行在改进的运行策略中，电负荷和热负荷都参与到需求响应中。原始的优化问题采用拉格朗日乘子法解耦全局功率平衡后，被分解为4个完全由就地的约束和目标函数组成的子问题，因此适合采用完全分布式一致性算法求解。该方法充分考虑能量供应侧的热电耦合和负荷侧的热电耦合，完全的分布式求解在保证用户隐私的基础上具有很高的效率，对于未来综合能源系统在配网层面广泛应用下的分布式能量管理具有很好的指导意义。

Description

一种基于一致性的综合能源系统能量管理方法

技术领域

本发明属于电力系统能量管理技术领域，尤其涉及一种基于一致性的综合能源系统能量管理方法。

背景技术

能源和环境问题影响着社会和经济的发展，综合能源系统在提高能源利用效率、减少环境污染等方面具有很大的优势，近年来综合能源系统得到广泛的研究和应用。相比于综合能源系统的规划、协调控制等，能量管理在减少运行成本方面备受研究者关注，传统的能量管理方法有一个共同的特征：都需要集中式优化，也就是说需要一个控制中心来收集并处理各个参与主体的信息，然而在配网层面，不同的主体属于不同的运营单位，因此考虑到信息的隐私性和参与主体归属的不确定性，集中式优化不再适用于我们的场景，相反的，分布式优化展现出更多的优势。综合能源系统的能量管理问题，有别于单纯的电力网络，网络中是由电网络和热网络组成的，且能量供应侧和负荷消耗侧都存在热和电的耦合，对于分布式的综合能源系统能量管理方法还有待进一步研究。

发明内容

针对上述问题，本发明提出一种基于一致性的综合能源系统能量管理方法，其特征在于，所述方法包括以下步骤：

步骤1：构建关于综合能源系统的社会福利最大化模型；

步骤2：对社会福利最大化模型进行解耦分析；

步骤3：采用完全分布式的一致性算法对解耦后的问题进行求解。

进一步，综合能源系统包括发电侧和需求侧两部分，其中发电侧参与主体包括变电站、分布式电源(光伏和风机)和热电联产(CHP)单元，负荷侧参与主体包括热负荷和电负荷。

进一步，社会福利最大化模型为：

目标函数：以任一时刻全社会各个参与主体总利润最大为目标，该目标函数的数学模型表示为：

max∑[P_SUB(P_i,p_e)+P_PV(P_i,p_e)+P_WT(P_i,p_e)+P_CHP(P_i,H_i,p_e,p_h)+P_d(P_i,H_i,p_e,p_h)]

其中，P_SUB(P_i,p_e)为变电站的利润，P_PV(P_i,p_e)为光伏的利润，P_WT(P_i,p_e)为风机的利润，P_CHP(P_i,H_i,p_e,p_h)为CHP机组的利润，P_d(P_i,H_i,p_e,p_h)为热和电负荷的利润；

所述总利润由总收益减去总成本，其中变电站的利润为卖电给配网内部用户的收益扣除向公共电网购电产生的成本，公共电网的电由火力发电机组提供，

C_i(P_i)＝a_iP_i ²+b_iP_i+c_i

其中，V_SUB为变电站组成的集合，p_e为电价，P_i为电功率，C_i(P_i)为火力发电机组的发电成本，a_i，b_i，c_i为相应的费用系数；

可再生能源发电无需发电成本，利润仅由卖电收益构成，

其中，V_PV为光伏组成的集合，V_WT为风机组成的集合；

CHP机组既可以产生热又可以产生电，利润如下，

其中，V_CHP为CHP机组组成的集合，p_h为热价，H_i为热功率，F_i(P_i,H_i)为CHP机组产生对应热和电的生产成本；具体的CHP生产成本为：

其中，φ，β，γ，δ，θ，为CHP机组相应的生产成本系数；

负荷侧既有电负荷又有热负荷，消耗热和电产生的效益采用用能的满意度来衡量，也称效用函数，

U_j(P_j,H_j)＝U_j(P_j)+U_j(H_j)

其中，V_d为负荷组成的集合，U_j(P_j,H_j)为负荷用能的满意度，U_j(P_j)为用电的满意度，U_j(H_j)为用热的满意度，ω_j,α_j为相应的系数；

决策变量：包括各个参与主体的电出力P_j和热出力H_j，整个综合能源系统的电价p_e和热价p_h；

约束条件：社会福利最大化模型中包含全局功率平衡约束和就地功率范围约束。

进一步，全局功率平衡约束包括电功率平衡约束和热功率平衡约束，就地功率范围约束包括变电站功率范围约束、风机出力范围约束、光伏出力范围约束、CHP机组出力范围约束和负荷需求范围约束。

进一步，对社会福利最大化模型进行解耦分析为：

采用拉格朗日乘子法将全局功率平衡约束加到原目标函数中，重新构造函数，从而达到全局功率平衡解耦的目的，其中，新的目标函数如下：

其中，V_el为电负荷组成的集合，V_tl为热负荷组成的集合，λ为电功率平衡对应的拉格朗日乘子，μ为热功率平衡对应的拉格朗日乘子；

约束条件：只包括就地功率范围约束；

决策变量：包括各个参与主体的电出力P_j和热出力H_j，整个综合能源系统的电价p_e和热价p_h。

进一步，重新定义此处的拉格朗日乘子为增量费用或者增量效用，当参与主体为能量供应源侧则称为增量费用，当参与主体为负荷能量消耗侧，则称为增量效用，此处的增量费用或者增量效用可通过对新的目标函数求偏导获得。

进一步，当CHP机组生产成本只由热功率表示时：

P_CHP(H_i,p_e,p_h)＝p_eP_i+p_hH_i-F_i(H_i)

负荷侧的耦合受人们的生产和生活的影响，满足一定的线性关系：

P_j＝k_jH_j

其中，k_j为相关系数；

热电耦合模型如下：

其中，η_e为电效率，η_r为热回收系数；

因此，原来的优化问题可以被分解为4个子问题，这4个子问题完全由就地的目标函数和约束条件组成，其中，4个子问题为：

1)CHP子问题：min F_i(P_i,H_i)-λ_iP_i-μ_iH_i

H_i ^min≤H_i≤H_i ^max

2)负荷子问题：min μ_jH_j+λ_j(P_j)-U_j(P_j,H_j)

s.t.P_j＝k_jH_j

H_j ^min≤H_j≤H_j ^max

P_j ^min≤P_j≤P_j ^max

3)变电站子问题：min C_i(P_i)-λ_iP_i

s.t.P_i ^min≤P_i≤P_i ^max

4)分布式电源子问题，包括光伏子问题和风机子问题，其中光伏子问题为：

min λ_iP_i

s.t.P_i ^min≤P_i≤P_i ^max

风机子问题为：

min λ_iP_i

s.t.P_i ^min≤P_i≤P_i ^max。

进一步，完全分布式求解算法为：ξ_i ^E(k)为电功率在第k次迭代中的就地偏差量，ξ_i ^H(k)为热功率在第k次迭代中的就地偏差量，一致性算法迭代过程如下：

1)初始化：设置整个综合能源系统中电力网和热力网的初始就地功率偏差ξ_i ^E(0)、ξ_i ^H(0)，初始增量费用和效用λ_i(0)、μ_i(0)，以及初始能量供应源的功率出力和负荷功率消耗P_i(0)、H_i(0)；

2)循环：更新增量费用和效用、更新能量供应源的功率出力和负荷的功率消耗、更新就地功率偏差、判断是否满足收敛条件；

3)输出：整个综合能源系统中各个参与主体的增量费用和效用λ_i(k)、μ_i(k)和功率出力和消耗P_i(k)、H_i(k)。

进一步，一致性算法迭代过程，所述更新增量费用和效用分为热力网和电力网两个，其中热力网为：

其中，V_t为热力网所有参与主体组成的集合，x_ij为热力网中通信拓扑对应的行随机矩阵的元素，ρ_t为(0，1)区间内的正反馈系数；

电力网中的更新增量费用和效用为：

其中，V_e为电力网所有参与主体组成的集合，r_ij为电力网中通信拓扑对应的行随机矩阵的元素，ρ_e为(0，1)区间内的正反馈系数；

进一步，能量供应源的功率出力和负荷的功率消耗分为热力网和电力网两个，其中热力网的能量供应源的功率出力为：

热力网中负荷的功率消耗为：

H_i(k+1)＝argmin[λ_i(k+1)k_iH_i(k)+μ_i(k+1)H_i(k)-U_i(H_i(k))],H_i(k)∈[H_i ^min,H_i ^max]

电力网中的能量供应源的功率出力由变电站、CHP、光伏和风机组成，其中变电站功率出力为：

P_i(k+1)＝argmin[C_i(P_i(k))-λ_i(k+1)P_i(k)],P_i(k)∈[P_i ^min,P_i ^max]

CHP功率出力由热力网功率平衡决定：

分布式电源，包括光伏和风机，都采用优先消纳的策略：

P_i(k+1)＝P_i ^max

电力网中负荷的功率消耗与热力网中负荷消耗有关：

P_i(k+1)＝k_iH_i(k+1),P_i(k)∈[P_i ^min,P_i ^max]。

进一步，更新就地功率偏差分为热力网和电力网两部分，其中热力网源侧新的就地功率偏差为：

热力网负荷侧新的就地功率偏差为：

其中，y_ij为热力网中通信拓扑对应的列随机矩阵的元素；

其中电力网源侧新的就地功率偏差为：

电力网负荷侧新的就地功率偏差为：

进一步，是否满足收敛条件，收敛条件为所有就地功率偏差的绝对值小于0.01，所有增量费用和效用临近的两次迭代结果之差的绝对值小于0.01，即，其中ε_μ，ε_ξE，ε_μ，ε_λ都取0.01。

本发明的有益效果在于：

本发明提出的基于一致性的综合能源系统能量管理方法，使综合能源系统里面的各个参与主体都达到利润最大化的状态。采用完全分布式的一致性求解算法，既保证了网络中各个参与主体的隐私，又提高了运行的效率。在综合能源系统中，能量供应侧和负荷需求侧都存在热电耦合关系，网络中的CHP机组采用改进的运行策略和热电负荷假定满足线性关系很好的解决了该双侧的热电耦合问题。

附图说明

附图1为本发明的一种基于一致性的综合能源系统能量管理方法流程图；

附图2为本发明的实施例中综合能源系统算例系统结构图；

附图3a-3d为本发明实施例中热力网迭代收敛结果；

附图4a-4d为本发明实施例中电力网迭代收敛结果；

附图5为本发明实施例中全社会福利结果迭代收敛曲线。

具体实施方式

下面结合附图和实施例对本发明进行详细说明。

本发明提出一种基于一致性的综合能源系统能量管理方法，将能量供应侧的热电耦合和负荷侧的热电耦合结合起来，并考虑了负荷的需求响应，最终达到全社会福利最大化，采用完全分布式的一致性算法求解，既保证了参与主体的隐私，也提高了计算的效率。具体方法如附图1所示。附图1为本发明的一种基于一致性的综合能源系统能量管理方法流程图，包括以下步骤：

步骤1：构建关于综合能源系统的社会福利最大化模型；

步骤2：对社会福利最大化模型进行解耦分析；

步骤3：采用完全分布式的一致性算法对解耦后的问题进行求解。

综合能源系统包括发电侧和需求侧两部分，其中发电侧参与主体包括变电站、分布式电源(光伏和风机)和热电联产(CHP)单元，负荷侧参与主体包括热负荷和电负荷。发电测和需求侧都参与到能量市场中以追求自己最大的利益。

社会福利最大化模型为：

目标函数：以任一时刻全社会各个参与主体总利润最大为目标，该目标函数的数学模型表示为：

max∑[P_SUB(P_i,p_e)+P_PV(P_i,p_e)+P_WT(P_i,p_e)+P_CHP(P_i,H_i,p_e,p_h)+P_d(P_i,H_i,p_e,p_h)]

其中，P_SUB(P_i,p_e)为变电站的利润，P_PV(P_i,p_e)为光伏的利润，P_WT(P_i,p_e)为风机的利润，P_CHP(P_i,H_i,p_e,p_h)为CHP机组的利润，P_d(P_i,H_i,p_e,p_h)为热和电负荷的利润。

所述总利润由总收益减去总成本，其中变电站的利润为卖电给配网内部用户的收益扣除向公共电网购电产生的成本，公共电网的电由火力发电机组提供，

C_i(P_i)＝a_iP_i ²+b_iP_i+c_i

其中，V_SUB为变电站组成的集合，p_e为电价，P_i为电功率，C_i(P_i)为火力发电机组的发电成本，a_i，b_i，c_i为相应的费用系数。

可再生能源发电无需发电成本，利润仅由卖电收益构成，

其中，V_PV为光伏组成的集合，V_WT为风机组成的集合。

CHP机组既可以产生热又可以产生电，利润如下，

其中，V_CHP为CHP机组组成的集合，p_h为热价，H_i为热功率，F_i(P_i,H_i)为CHP机组产生对应热和电的生产成本。具体的CHP生产成本为：

其中，φ，β，γ，δ，θ，为CHP机组相应的生产成本系数。

负荷侧既有电负荷又有热负荷，消耗热和电产生的效益采用用能的满意度来衡量，也称效用函数，

U_j(P_j,H_j)＝U_j(P_j)+U_j(H_j)

其中，V_d为负荷组成的集合，U_j(P_j,H_j)为负荷用能的满意度，U_j(P_j)为用电的满意度，U_j(H_j)为用热的满意度，ω_j,α_j为相应的系数。

决策变量：包括各个参与主体的电出力P_j和热出力H_j，整个综合能源系统的电价p_e和热价p_h；

约束条件：包括电功率平衡约束、热功率平衡约束、变电站功率范围约束、风机出力范围约束、光伏出力范围约束、CHP机组出力范围约束、负荷需求范围约束。

因此，模型可以写成一个优化问题：

P_i ^min≤P_i≤P_i ^max,i∈V_SUB

P_i ^min≤P_i≤P_i ^max,i∈V_CHP

P_j ^min≤P_j≤P_j ^max,j∈V_el

H_j ^min≤H_j≤H_j ^max,j∈V_tl

P_i ^min≤P_i≤P_i ^max,i∈V_PV

P_i ^min≤P_i≤P_i ^max,i∈V_WT

其中，P^min，P^max分别为电功率下限和上限，H^min，H^max分别为热功率下限和上限，下标i代表能量供应侧，下标j代表能量消耗侧，V_tl为热负荷组成的集合，V_el为电负荷组成的集合。

将电功率平衡和热功率平衡带入目标函数，化简得到：

因此得到如下凸优化问题：

P_i ^min≤P_i≤P_i ^max,i∈V_SUB

P_i ^min≤P_i≤P_i ^max,i∈V_CHP

P_j ^min≤P_j≤P_j ^max,j∈V_el

H_j ^min≤H_j≤H_j ^max,j∈V_tl

P_i ^min≤P_i≤P_i ^max,i∈V_PV

P_i ^min≤P_i≤P_i ^max,i∈V_WT

对社会福利最大化模型进行解耦分析为：

社会福利最大化模型中包含全局功率平衡约束(包括电功率平衡约束和热功率平衡约束)和就地功率范围约束(包括变电站功率范围约束、风机出力范围约束、光伏出力范围约束、CHP机组出力范围约束和负荷需求范围约束)，采用拉格朗日乘子法将全局功率平衡约束加到原目标函数中，重新构造函数，从而达到全局功率平衡解耦的目的。其中，新的目标函数如下：

其中，V_el为电负荷组成的集合，V_tl为热负荷组成的集合，λ为电功率平衡对应的拉格朗日乘子，μ为热功率平衡对应的拉格朗日乘子；

约束条件：只包括就地功率范围约束，即变电站功率范围约束、风机出力范围约束、光伏出力范围约束、CHP机组出力范围约束和负荷需求范围约束；

决策变量：包括各个参与主体的电出力P_j和热出力H_j，整个综合能源系统的电价p_e和热价p_h。

该综合能源系统中每个参与主体都对应着一个拉格朗日乘子，因此网络被分解成热力网和电力网来求解，由于该处的拉格朗日乘子很巧妙的代替了原来的电价和热价，因此重新定义此处的拉格朗日乘子为增量费用或者增量效用，当参与主体为能量供应源则称为增量费用，当参与主体为负荷能量消耗侧，则称为增量效用，此处的增量费用或者增量效用可通过对新的目标函数求偏导获得；通过CHP机组改进的运行策略和假设负荷侧的热负荷和电负荷线性相关，我们很好的解决了能量源侧CHP机组的热电耦合和负荷侧热电负荷的耦合。其中传统的CHP机组运行在以热定电或者以电定热的运行策略上，针对每个CHP机组，其热电耦合模型如下：

其中，η_e为电效率，η_r为热回收系数。

综合能源系统在配网层面上，各个CHP机组只需满足全网的热力和电力的能量平衡，机组热和电的出力不受以热定电或以电定热运行策略的影响，运行在改进的运行策略上：由于CHP机组是综合能源系统中唯一的热源，因此，每个CHP机组在满足全局热功率平衡的基础上产生自己的热功率来最大化自己的利益。当CHP机组生产成本只由热功率表示时：

P_CHP(H_i,p_e,p_h)＝p_eP_i+p_hH_i-F_i(H_i)

负荷侧的耦合受人们的生产和生活的影响，满足一定的线性关系：

P_j＝k_jH_j

其中，k_j为相关系数。

因此，当负荷的效用模型只由电功率表示时：

U_j ^*(H_j)＝U_j(k_jH_j)+U_j(H_j)

当负荷的效用模型只由电功率表示时：

所述的增量费用(效用)分为热力网的和电力网的，定义热力网的增量费用(效用)为：

定义电力网的增量费用(效用)为：

因此，原来的优化问题可以被分解为4个子问题，这4个子问题完全由就地的目标函数和约束条件组成，其中，4个子问题为：

1)CHP子问题：min F_i(P_i,H_i)-λ_iP_i-μ_iH_i

H_i ^min≤H_i≤H_i ^max

2)负荷子问题：min μ_jH_j+λ_j(P_j)-U_j(P_j,H_j)

s.t.P_j＝k_jH_j

H_j ^min≤H_j≤H_j ^max

P_j ^min≤P_j≤P_j ^max

3)变电站子问题：min C_i(P_i)-λ_iP_i

s.t.P_i ^min≤P_i≤P_i ^max

4)分布式电源子问题，包括光伏子问题和风机子问题，其中光伏子问题为：

min λ_iP_i

s.t.P_i ^min≤P_i≤P_i ^max

风机子问题为：

min λ_iP_i

s.t.P_i ^min≤P_i≤P_i ^max

采用完全分布式的一致性算法对全局功率平衡解耦后的问题进行求解为，采用一致性算法进行求解，该算法是完全分布式的，没有任何求解中心；且该算法在强连接的有向图中具有很好的收敛性。其中，完全分布式求解算法为，

为保证全局功率平衡的收敛性，需要引入新的变量来代表就地功率偏差量，其中ξ_i ^E(k)为电功率在第k次迭代中的就地偏差量，ξ_i ^H(k)为热功率在第k次迭代中的就地偏差量。一致性算法迭代过程如下：

1)初始化：设置整个综合能源系统中电力网和热力网的初始就地功率偏差ξ_i ^E(0)，ξ_i ^H(0)、初始增量费用(效用)λ_i(0)，μ_i(0)和初始能量供应源的功率出力和负荷功率消耗P_i(0)，H_i(0)；

2)循环：更新增量费用(效用)—>更新能量供应源的功率出力和负荷的功率消耗—>更新就地功率偏差—>判断是否满足收敛条件；

3)输出：整个综合能源系统中各个参与主体的增量费用(效用)λ_i(k)，μ_i(k)和功率出力(消耗)P_i(k)，H_i(k)。

设置整个综合能源系统中电力网和热力网的初值，热力网中满足：

H_i(0)＝0,i∈V_t

ξ^H _i(0)＝0,i∈V_t

电力网中满足：

P_i(0)＝0,i∈V_e

ξ_i ^E(0)＝0,i∈V_e

其中，V_t为热力网所有参与主体组成的集合，V_e为电力网所有参与主体组成的集合。

更新增量费用(效用)分为热力网和电力网两个，其中热力网为：

其中，x_ij为热力网中通信拓扑对应的行随机矩阵的元素，ρ_t为(0，1)区间内的正反馈系数。

电力网中的更新增量费用(效用)为：

其中，r_ij为电力网中通信拓扑对应的行随机矩阵的元素，ρ_e为(0，1)区间内的正反馈系数。

能量供应源的功率出力和负荷的功率消耗分为热力网和电力网两个，其中热力网的能量供应源的功率出力为：

热力网中负荷的功率消耗为：

H_i(k+1)＝argmin[λ_i(k+1)k_iH_i(k)+μ_i(k+1)H_i(k)-U_i(H_i(k))],H_i(k)∈[H_i ^min,H_i ^max]

电力网中的能量供应源的功率出力由变电站、CHP、光伏和风机组成，其中变电站功率出力为：

P_i(k+1)＝argmin[C_i(P_i(k))-λ_i(k+1)P_i(k)],P_i(k)∈[P_i ^min,P_i ^max]

CHP功率出力由热力网功率平衡决定：

分布式电源，包括光伏和风机，都采用优先消纳的策略：

P_i(k+1)＝P_i ^max

电力网中负荷的功率消耗与热力网中负荷消耗有关：

P_i(k+1)＝k_iH_i(k+1),P_i(k)∈[P_i ^min,P_i ^max]

更新就地功率偏差分为热力网和电力网两部分，其中热力网源侧新的就地功率偏差为：

热力网负荷侧新的就地功率偏差为：

其中，y_ij为热力网中通信拓扑对应的列随机矩阵的元素。

其中电力网源侧新的就地功率偏差为：

电力网负荷侧新的就地功率偏差为：

是否满足收敛条件，收敛条件为所有就地功率偏差的绝对值小于0.01，所有增量费用(效用)临近的两次迭代结果之差的绝对值小于0.01，即，其中ε_μ，ε_μ，ε_λ都取0.01。

下面以一个实施例来说明本发明的方法。

本实施例采用改进的IEEE-33节点综合能源系统进行详细说明，如图2所示的综合能源系统算例系统结构图，图中包括32个负荷节点、2个光伏、2个风机和3个CHP机组，通信拓扑沿实际综合能源系统物理线路配置，整个网络由热力网和电力网组成。本算例的能量源侧的相关参数如表1所示,负荷侧的相关参数如表2所示。

表1 能量源侧的相关参数

表2 负荷侧的相关参数(功率单位为kW)

附图3为热力网迭代收敛结果，其中，附图3a-3c分别表示热力网的增量费用(效用)、就地功率出力(消耗)和就地功率偏差。图3(a)可以看出增量费用(效用)随着迭代最终收敛到6.356美分/kWh，与此同时，图3(b)可以看出就地功率偏差趋向于0；就地功率出力在图3(c)中表示为负的，就地功率消耗在图3(c)中表示为正的；随着迭代次数增加，可以看出热力网的源侧和负荷侧的总的功率偏差接近于0，这个结论可以进一步在图3(d)中得到验证。更重要的，所有的增量费用(效用)都趋向于同一个值，在这个值下，所有的参与主体都达到了最优，因此热力网问题求解完毕。

附图4为电力网迭代收敛结果。从图4(a)我们可以看出电力网增量费用(效用)的最终收敛值在第850次迭代获得，为13.41美分/kWh，从图4(b)看出，就地的功率偏差随着迭代次数的增加也逐渐趋向于0。在热力网的迭代结果基础上，我们很容易获得图4(c)所示的CHP机组电出力结果，也用负的表示以区别于正的功率消耗；变电站节点一开始为0，逐渐增加出力达到最终的575kW。与此同时，电力网通过优先消纳光伏和风机产生的电来充分利用可再生能源。图4(d)中所示的总电力需求和总电力消耗接近于相等，表明电功率平衡得到满足。

附图5为全社会福利结果迭代收敛曲线，从曲线上我们可以看出，随着迭代次数的增加，在热力网和电力网达到稳定的同时，全社会的福利也逐渐增大并趋向于最终的2.007*10⁴美元。因此，整个综合能源系统的能量管理问题得到了解决。

此实施例仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。

Claims (12)

1.一种基于一致性的综合能源系统能量管理方法，其特征在于，所述方法包括以下步骤：

步骤1：构建关于综合能源系统的社会福利最大化模型；

步骤2：对社会福利最大化模型进行解耦分析；

步骤3：采用完全分布式的一致性算法对解耦后的问题进行求解。

2.根据权利要求1所述的一种基于一致性的综合能源系统能量管理方法，其特征在于，所述综合能源系统包括发电侧和需求侧两部分，其中发电侧参与主体包括变电站、分布式电源(光伏和风机)和热电联产(CHP)单元，负荷侧参与主体包括热负荷和电负荷。

3.根据权利要求2所述的一种基于一致性的综合能源系统能量管理方法，其特征在于，所述社会福利最大化模型为：

目标函数：以任一时刻全社会各个参与主体总利润最大为目标，该目标函数的数学模型表示为：

max∑[P_SUB(P_i,p_e)+P_PV(P_i,p_e)+P_WT(P_i,p_e)+P_CHP(P_i,H_i,p_e,p_h)+P_d(P_i,H_i,p_e,p_h)]

其中，P_SUB(P_i,p_e)为变电站的利润，P_PV(P_i,p_e)为光伏的利润，P_WT(P_i,p_e)为风机的利润，P_CHP(P_i,H_i,p_e,p_h)为CHP机组的利润，P_d(P_i,H_i,p_e,p_h)为热和电负荷的利润；

所述总利润由总收益减去总成本，其中变电站的利润为卖电给配网内部用户的收益扣除向公共电网购电产生的成本，公共电网的电由火力发电机组提供，

C_i(P_i)＝a_iP_i ²+b_iP_i+c_i

其中，V_SUB为变电站组成的集合，p_e为电价，P_i为电功率，C_i(P_i)为火力发电机组的发电成本，a_i，b_i，c_i为相应的费用系数；

可再生能源发电无需发电成本，利润仅由卖电收益构成，

其中，V_PV为光伏组成的集合，V_WT为风机组成的集合；

CHP机组既可以产生热又可以产生电，利润如下，

其中，V_CHP为CHP机组组成的集合，p_h为热价，H_i为热功率，F_i(P_i,H_i)为CHP机组产生对应热和电的生产成本；具体的CHP生产成本为：

其中，φ，β，γ，δ，θ，为CHP机组相应的生产成本系数；

负荷侧既有电负荷又有热负荷，消耗热和电产生的效益采用用能的满意度来衡量，也称效用函数，

U_j(P_j,H_j)＝U_j(P_j)+U_j(H_j)

其中，V_d为负荷组成的集合，U_j(P_j,H_j)为负荷用能的满意度，U_j(P_j)为用电的满意度，U_j(H_j)为用热的满意度，ω_j,α_j为相应的系数；

决策变量：包括各个参与主体的电出力P_j和热出力H_j，整个综合能源系统的电价p_e和热价p_h；

约束条件：社会福利最大化模型中包含全局功率平衡约束和就地功率范围约束。

4.根据权利要求3所述的一种基于一致性的综合能源系统能量管理方法，其特征在于，全局功率平衡约束包括电功率平衡约束和热功率平衡约束，就地功率范围约束包括变电站功率范围约束、风机出力范围约束、光伏出力范围约束、CHP机组出力范围约束和负荷需求范围约束。

5.根据权利要求3所述的一种基于一致性的综合能源系统能量管理方法，其特征在于，所述对社会福利最大化模型进行解耦分析为：

采用拉格朗日乘子法将全局功率平衡约束加到原目标函数中，重新构造函数，从而达到全局功率平衡解耦的目的，其中，新的目标函数如下：

其中，V_el为电负荷组成的集合，V_tl为热负荷组成的集合，λ为电功率平衡对应的拉格朗日乘子，μ为热功率平衡对应的拉格朗日乘子；

约束条件：只包括就地功率范围约束；

决策变量：包括各个参与主体的电出力P_j和热出力H_j，整个综合能源系统的电价p_e和热价p_h。

6.根据权利要求5所述的一种基于一致性的综合能源系统能量管理方法，其特征在于，重新定义此处的拉格朗日乘子为增量费用或者增量效用，当参与主体为能量供应源侧则称为增量费用，当参与主体为负荷能量消耗侧，则称为增量效用，此处的增量费用或者增量效用可通过对新的目标函数求偏导获得。

7.根据权利要求5所述的一种基于一致性的综合能源系统能量管理方法，其特征在于，当CHP机组生产成本只由热功率表示时：

P_CHP(H_i,p_e,p_h)＝p_eP_i+p_hH_i-F_i(H_i)

负荷侧的耦合受人们的生产和生活的影响，满足一定的线性关系：

P_j＝k_jH_j

其中，k_j为相关系数；

热电耦合模型如下：

其中，η_e为电效率，η_r为热回收系数；

因此，原来的优化问题可以被分解为4个子问题，这4个子问题完全由就地的目标函数和约束条件组成，其中，4个子问题为：

1)CHP子问题：min F_i(P_i,H_i)-λ_iP_i-μ_iH_i

H_i ^min≤H_i≤H_i ^max

2)负荷子问题：minμ_jH_j+λ_j(P_j)-U_j(P_j,H_j)

s.t.P_j＝k_jH_j

H_j ^min≤H_j≤H_j ^max

P_j ^min≤P_j≤P_j ^max

3)变电站子问题：min C_i(P_i)-λ_iP_i

s.t.P_i ^min≤P_i≤P_i ^max

4)分布式电源子问题，包括光伏子问题和风机子问题，其中光伏子问题为：

min λ_iP_i

s.t.P_i ^min≤P_i≤P_i ^max

风机子问题为：

minλ_iP_i

s.t.P_i ^min≤P_i≤P_i ^max。

8.根据权利要求7所述的一种基于一致性的综合能源系统能量管理方法，其特征在于，完全分布式求解算法为：ξ_i ^E(k)为电功率在第k次迭代中的就地偏差量，为热功率在第k次迭代中的就地偏差量，一致性算法迭代过程如下：

1)初始化：设置整个综合能源系统中电力网和热力网的初始就地功率偏差ξ_i ^E(0)、初始增量费用和效用λ_i(0)、μ_i(0)，以及初始能量供应源的功率出力和负荷功率消耗P_i(0)、H_i(0)；

2)循环：更新增量费用和效用、更新能量供应源的功率出力和负荷的功率消耗、更新就地功率偏差、判断是否满足收敛条件；

3)输出：整个综合能源系统中各个参与主体的增量费用和效用λ_i(k)、μ_i(k)和功率出力和消耗P_i(k)、H_i(k)。

9.据权利要求8所述的一种基于一致性的综合能源系统能量管理方法，其特征在于，所述的一致性算法迭代过程，所述更新增量费用和效用分为热力网和电力网两个，其中热力网为：

其中，V_t为热力网所有参与主体组成的集合，x_ij为热力网中通信拓扑对应的行随机矩阵的元素，ρ_t为(0，1)区间内的正反馈系数；

电力网中的更新增量费用和效用为：

其中，V_e为电力网所有参与主体组成的集合，r_ij为电力网中通信拓扑对应的行随机矩阵的元素，ρ_e为(0，1)区间内的正反馈系数。

10.据权利要求8所述的一种基于一致性的综合能源系统能量管理方法，其特征在于，所述能量供应源的功率出力和负荷的功率消耗分为热力网和电力网两个，其中热力网的能量供应源的功率出力为：

热力网中负荷的功率消耗为：

H_i(k+1)＝argmin[λ_i(k+1)k_iH_i(k)+μ_i(k+1)H_i(k)-U_i(H_i(k))],H_i(k)∈[H_i ^min,H_i ^max]

电力网中的能量供应源的功率出力由变电站、CHP、光伏和风机组成，其中变电站功率出力为：

P_i(k+1)＝argmin[C_i(P_i(k))-λ_i(k+1)P_i(k)],P_i(k)∈[P_i ^min,P_i ^max]

CHP功率出力由热力网功率平衡决定：

分布式电源，包括光伏和风机，都采用优先消纳的策略：

P_i(k+1)＝P_i ^max

电力网中负荷的功率消耗与热力网中负荷消耗有关：

P_i(k+1)＝k_iH_i(k+1),P_i(k)∈[P_i ^min,P_i ^max]。

11.据权利要求8所述的一种基于一致性的综合能源系统能量管理方法，其特征在于，所述更新就地功率偏差分为热力网和电力网两部分，其中热力网源侧新的就地功率偏差为：

热力网负荷侧新的就地功率偏差为：

其中，y_ij为热力网中通信拓扑对应的列随机矩阵的元素；

其中电力网源侧新的就地功率偏差为：

电力网负荷侧新的就地功率偏差为：

12.据权利要求8所述的一种基于一致性的综合能源系统能量管理方法，其特征在于，所述是否满足收敛条件，收敛条件为所有就地功率偏差的绝对值小于0.01，所有增量费用和效用临近的两次迭代结果之差的绝对值小于0.01，即，其中ε_μ，ε_ξE，ε_μ，ε_λ都取0.01。