Distributed Energy 2018 , 3 (4): 1-8 https://doi.org/10.16513/j.cnki.10-1427/tk.2018.04.001

分布式能源系统优化运行

基于能量网络理论的区域综合能源系统运行优化研究

李明^,, 陈皓勇, 邱明, 黄龙, 陈思敏

华南理工大学电力学院，广东广州 510641

Operation Optimization of Regional Integrated Energy System Based on Energy Network Theory

LI Ming, CHEN Haoyong, QIU Ming, HUANG Long, CHEN Simin

School of Electric Power, South China University of Technology, Guangzhou 510641, Guangdong Province, China

中图分类号: TK 01

文献标识码: A

文章编号: 2096-2185(2018)04-0001-08

责任编辑: LI MingCHEN HaoyongQIU MingHUANG LongCHEN Simin

收稿日期: 2018-05-02

基金资助: 国家重点基础研究发展计划项目(973项目)(2016YFB0900102)

作者简介:

李　明(1996—)，男，硕士研究生，研究方向为能源互联网，1154954144@qq.com；

陈皓勇(1975—)，男，博士，教授，研究方向为电力系统调度、规划与控制，可再生能源发电与综合消纳技术等，eehychen@scut.edu.cn；

邱　明(1993—)，男，硕士研究生，研究方向为能源互联网；

黄　龙(1994—)，男，硕士研究生，研究方向为能源互联网；

陈思敏(1996—)，女，硕士研究生，研究方向为能源互联网。

展开

摘要

综合能源系统通过电、气、冷/热等多种形式能源的相互耦合、转化和资源分配，实现了能量梯级利用，提升能源综合利用率以及能源供应的灵活度。然而，不同类型的能源系统其能流性质以及传输特性各异，这阻碍了综合能源系统的统一建模和整体分析，难以实现和发挥综合能源系统的整体协同。从能量网络理论出发，基于能量公理化的思想，建立流体网络以及热力网络传输的等值电路模型；结合能量转换设备模型，完成多能流耦合系统的统一建模。随后，以系统总能量损失最小为目标，建立区域综合能源系统运行优化模型，通过粒子群算法对系统内部控制变量进行决策寻优。最后，对某岛屿32热力节点综合能源系统进行建模以及运行优化，验证所提模型以及优化方法的有效性。

关键词： 综合能源 ; 能量网络 ; 统一建模 ; 公理化 ; 运行优化

Abstract

Coupling various forms of energy, such as electricity, gas, cold/heat and other forms of energy, the integrated energy system has realized the cascade utilization of energy, enhanced the comprehensive utilization of energy and the flexibility of energy supply through the energy conversion and resource allocation between its sub-systems. However, the energy flows in different types of energy systems vary in their transmission characteristics, which hinders the unified modeling and overall analysis of the integrated energy system, so it is difficult to achieve and bring into play the overall coordination of the integrated energy system. In this paper, based on the theory of energy network and the axiomatic energy theory, the electrical equivalent model of the transmission line in thermal network and the fluid network are established. Combined with the energy conversion equipment model, the unified modeling of the multi-energy flow systems is completed. Then, taking the minimum total energy loss of the system as the goal, a regional integrated energy system operation optimization model is established, in which the particle swarm optimization algorithm is used to optimize the internal control variables of the system. Finally, an integrated energy system of an island with 32 thermal nodes is modeled and optimized, which verifies the validity of the proposed model and optimization method.

Keywords： integrated energy ; energy network ; unified modeling ; axiomatic ; operation optimization

PDF (1257KB) 元数据多维度评价相关文章收藏文章

本文引用格式导出 EndNote Ris Bibtex

李明, 陈皓勇, 邱明, 黄龙, 陈思敏. 基于能量网络理论的区域综合能源系统运行优化研究[J]. , 2018, 3(4): 1-8 https://doi.org/10.16513/j.cnki.10-1427/tk.2018.04.001

LI Ming, CHEN Haoyong, QIU Ming, HUANG Long, CHEN Simin. Operation Optimization of Regional Integrated Energy System Based on Energy Network Theory[J]. Distributed Energy, 2018, 3(4): 1-8 https://doi.org/10.16513/j.cnki.10-1427/tk.2018.04.001

0　引言

随着经济科技的迅速发展，人类社会对于能源的需求急剧上升。这直接导致了煤炭、石油等常规化石能源资源的迅速枯竭，并由此引发了全球范围内人类对于资源环境能否持续发展的担忧。世界各国开始探索一种能够提高能源的利用效率，适应可再生能源的开发利用并优化各类能源的统筹分配以实现可持续发展的新途径^[1,2]。综合能源系统(integrated energy system, IES)以其能量转化设备(如热电联产机组、热泵及循环泵等)为纽带，通过供热网络、电力网络及供气网络等能源子网间的能量流动及互补，可有效实现能流的合理分配，并提升能量的利用效率，无疑将成为未来能源领域的发展方向^[3]。

国内外对于综合能源系统的研究，目前主要集中于IES耦合多能源系统的稳态潮流求解、IES优化规划方法及IES市场经济理论等几个方面。文献[4]基于能源集线器模型，构建了热电联产系统模型，针对区域综合能源系统(regional integrated energy system，RIES)的完全耦合、部分耦合及完全解耦3种运行模式，形成了适用于含电/气/热等3种能流的RIES稳态潮流求解方法。文献[5]建立了考虑安全约束的电－气最优潮流，该模型在求解最优潮流时考虑了系统发生N－1故障的情况，通过灵敏度分析可快速得到稳定的故障应对方案，避免线路过负荷。文献[6]考虑风力发电并网带来的不确定性以及天然气管道气压安全约束问题，建立了一个电－气耦合IES的二阶段鲁棒发电日前调度模型，通过分段线性化方法处理了燃气管网气流约束的非线性，提高了算法的收敛性。文献[7]针对IES的设备选型、设备容量配置及系统运行策略的整体优化，设计了一种冷热电联供系统的三级协同整体优化方法，并提出一种“以电定热”和“以热定电”混合运行策略，实现了系统产能以及用能的良好匹配。而在市场经济理论部分，吴利兰等人提出了一个基于多主多从Stackelberg博弈的能源交易模型，通过分析IES中多个分布式能源站和用户之间的多种能源的交互方式，求解它们之间的均衡交互策略，通过竞争决定能源价格，从而实现各自收益的最大化^[8]。

在IES理论研究方面，大量学者开展了深入的研究，并取得了丰硕的成果。然而，当前大多数的研究通常只涉及2种能量的耦合，同时各网络模型均是基于分立的异质能流系统建模机制而开展的，不同的能流领域沿用本领域当前典型的能量潮流计算方法。不同领域的模型互换性差，数据表达不统一，在时间尺度统一及数据交换方面存在困难，影响IES的统一建模及分析。为此，本文从能量网络理论出发，基于能量公理化的思想，建立流体网络及热力网络传输的等值电路模型；随后结合能量转换设备模型，实现多能流耦合系统的统一建模，并以IES稳态能量潮流作为IES优化模型的等式约束进行优化；最后，以所搭建的模型对某岛屿RIES进行建模，同时以RIES能量损失最少为目标，通过粒子群寻优算法实现该岛屿RIES的运行决策优化。

1　IES网络等值模型

1.1　能量网络理论

作为能源互联网的关键组成部分，综合能源系统是典型的多能域耦合系统，是能量网络理论的研究对象。能量网络理论从能量角度出发，以能量作为不同物理过程之间相互耦合的桥梁，研究不同形式能量的传递、存储及转化的普适化规律，应用集总参数的分析方法，实现异质能流子网络的融合分析^[9]。

为实现不同能流系统的统一描述，能量网络理论将综合能源系统的异质能流子系统的状态变量统一使用广延量及强度量来表述：广延量指在状态变量中与物质的量密切联系的，且具有可叠加性的量，典型的广延量包括熵、体积等；而强度量则为状态变量中不随物体物质的量变化而变化的，且不具有加和性的量，典型的强度量包括温度、电压等。

考虑能量传递过程是1维的，且传输过程中不发生泄漏，以工程中常用的圆管传输为例，能量的空间平衡方程可简化为

(1) $\frac{d E}{d t} = (X_{A} H_{x A} - X_{E} H_{x E}) + \int_{A}^{E} H_{x} d X + \int_{A}^{E} X d x_{g}$

式中：E为管道中的能量；H_x为流经截面的广延量x的流率，为流量密度J_x和管(线)面积之积；X为强度量；下标A和E分别为传递线(管)路的首端和末端截面积；x_g为传递过程中自发产生的、与长度L相关的广延量的函数。

式(1)中，等式左边的表达式是该圆柱形传递线(管)路中的能量因广延量流动产生的变化率；右边的第1项和第2项是流经管道横截面的广延量所携带的能量之和，第3项是传递过程中因过程不可逆出现的能量耗散，第4项是管道中自发生成的广延量产生过程中释放的能量。此外，由分析力学的基本知识可知，能量传递转换过程中还会产生广义动能(惯性现象)和势能(弹性现象)及其相互转化^[10]。在稳态传输情况下，暂时忽略能量传递的弹性和惯性现象。

1.2　传输管线等值模型

虽然不同形式的能量传递过程从物质形态角度来看存在巨大差异，但从介质流的角度来看却呈现一定的规律性，如描述导电过程的欧姆定律、描述粘性流动的牛顿定律和描述热传导的傅里叶定律，其形式具有一致性^[11]。传递公理指出，在相应基本强度量差的推动下，基本广延量实现流动运输^[12]。广延量的传递规律可描述为

(2) $J_{x} = K F = - K \nabla X$

式中：F为促使基本广延量流动的驱动力；K为基本广延量的传递系数；∇X为共轭的基本强度量之差。

1.2.1　电能传递的能量网络等值模型

对于电力网络，能量流密度J_x实际上就是电流的密度，强度量X为电势。因此，当系统以直流形式平稳供电时，根据式(2)，可求得电流的通式：

(3) $H_{x} = J_{x} A = - K A \nabla X \Rightarrow H_{x} d L = - K A d X$

式中L为传输线(管)路的长度。再对式(3)两边进行积分，于是有

(4) $H_{x} L = - K A (X_{A} - X_{E}) \Rightarrow H_{x} = \frac{K A (X_{A} - X_{E})}{L}$

如此，式(4)的含义就十分明显了：H_x为电流，而(X_A－X_E)则为电线首末端的电压差；K为基本广延量的传递系数，在此表示电导率；A为导线的面积。这便是熟悉的欧姆定律。

由于电力系统在这方面的研究比较成熟，本文直接沿用交流电力系统，即动态情况下的能量网络的计算方法。其能量网络节点方程如下式所示，未知量为节点的强度量。

(5) $H = Y X$

式中：Y为能量网络节点导纳矩阵；H和X分别为支路电流向量和支路电压向量。

1.2.2　压能传递的能量网络等值模型

压能传递过程中，强度量是压强p，基本广延量流是体积流率J_v。由质量守恒定律可知，质量的源强度g_m＝0。结合比体积V与密度ρ的关系式V＝ρ^－1，推导可得随体形式的体积平衡方程^[13]：

(6) $ρ \frac{d V}{d t} = \nabla \cdot J_{v}$

同理，可对式(2)作如下变换：

(7) $\begin{matrix} H_{x} = - K A \frac{d X}{d L} \Rightarrow H_{x} \int \frac{1}{K A} d L = X_{A} - X_{E} \Rightarrow \\ H_{x} R = X_{A} - X_{E}, R = \int \frac{1}{K A} d L \end{matrix}$

特别地，当基本广延量传递系数K为定值时(如稳态层流状态下的压能传递)，线(管)路的传阻可进一步简化为 $R = \frac{L}{K A}$ ^[14]。

根据流体网络理论，有如下方程组：

(8) ${\begin{cases} A H = Φ \\ B_{f} Δ X = 0 \end{cases}$

式中：A为流体网络的关联矩阵；B_f为流体网络基本回路矩阵；H为支路流量；Φ为流经负荷或热源的流量；ΔX为支路两端的压强差。

1.2.3　热能传递的能量网络等值模型

热能传递过程中，以温度为强度量。同时由热力学第二定律可知，在一个与外界不存在相互作用的系统中，熵永不自动减少。在不可逆过程中，系统由不平衡状态过渡到平衡状态，熵单调增加，即熵的源强度g_s≠0^[15]。此时以熵作为基础广延量将带来诸多不便。

定义H_e＝XH_x为能流(在热能传递过程中H_e即热流量)，λ＝KX为热传输系数，并根据连续性原理，假设能量的耗散现象只存在于管线的径向，则可对式(2)作如下变换：

(9) $\begin{matrix} H_{x} = - \frac{K X A}{X} \frac{d X}{d r} \Rightarrow H_{e} \int_{r_{1}}^{r_{2}} \frac{1}{λ A} d r = X_{r_{1}} - X_{r_{2}} \Rightarrow \\ H_{e} R_{e} = X_{r_{1}} - X_{r_{2}}, R_{e} = \int_{r_{1}}^{r_{2}} \frac{1}{λ A} d r \end{matrix}$

式中：下标e为热耗散；下标r₁与r₂分别为热传输管道的内径与外径；R_e即为热流量传输的热阻。

取圆柱形管路始端温度X_in及末端温度X_out的平均值作为管道平均温度X_av，可得

(10) $X_{av} = \frac{X_{in} + X_{out}}{2}$

根据式(9)所描述的管道内部的热量平衡关系，有

(11) $c_{p} m (X_{in} - X_{out}) = \frac{X_{av} - X_{0}}{R_{e}}$

式中：c_p为水工质的比热容；m为质量流量；下标in和out分别代表流入与流出；X₀为环境温度。

联合式(10)、(11)，得

(12) $X_{out} = \frac{2 c_{p} m R_{e} - 1}{2 c_{p} m R_{e} + 1} X_{in} + \frac{2}{2 c_{p} m R_{e} + 1} X_{0}$

由于热能传递同压能，电能等传输过程相差甚远，在能流汇集节点，当来自不同管道热流载体混合时，热流载体的温度将迅速统一。其热量平衡关系可描述为

(13) $(\sum m_{out}) X_{out} = \sum (m_{in} X_{in})$

1.3　能量耦合元件模型

1) 热电联产机组是当前应用较为广泛的一类能量转换设备，它通过合理分配汽轮机内高压蒸汽及乏汽的去向，使得发电和供热能同时进行。热电联产机组一般分为燃气轮机式机组、往复式内燃机机组和抽汽式机组。其中，燃气轮机式、往复式内燃机等热电联产机组的功率表达式为

(14) $c_{m} = \frac{ϕ}{P}$

该类机组为定热电比，c_m为热电比率，ϕ为机组输出热量，P为机组输出电功率。

抽汽式热电联产机组的功率表达式由下式描述：

(15) $Z = \frac{ϕ}{P_{con} - P}$

该类机组为变热电比，p_con为全冷凝模式下的输出电功率，Z为常数。

2) 管道中热水的流动需要克服阻力，消耗压能。为维持一定的压强，需要循环水泵消耗的电功率P_p为

(16) $P_{p} = \frac{m_{p} g H_{p}}{10^{6} η_{p}}$

式中：m_p为送水量；η_p为泵效率；H_p为网络的泵压；g为重力加速度。

3) 热泵消耗少量的电能把环境中的低品位热能提升至高品位，热泵的能效系数表达式为

(17) $C_{OP} = \frac{ϕ_{rb}}{P_{rb}}$

式中：ϕ_rb为热泵的供热功率；P_rb为热泵的功率消耗值；C_OP为热泵的能效系数，其值大小随热源和热负荷之间温差的变化而变化。

2　IES优化模型与求解

2.1　优化模型

本文采用粒子群算法针对综合能源系统进行优化，模型的优化目标为系统的总能量损失最小，此外，综合能源系统中还应该对变量进行约束。

1) 目标函数。

本文基于单目标函数来优化综合能源系统，以系统的总能量损失最小为目标函数。

(18) ${\begin{cases} \min S_{loss} = P_{loss} + Q_{loss} + L_{loss} \\ P_{loss} = \sum P_{in} - \sum P_{out} \\ Q_{loss} = \sum Q_{in} - \sum Q_{out} \\ L_{loss} = 2 \sum_{i \in S_{hl}} R_{V i} H_{V i}^{2} \end{cases}$

式中：P_loss为电力系统的有功损失，其值等于输入功率之和与输出功率之和的差值；Q_loss为供热系统的热能损失，其值等于输入热能之和与输出热能之和的差值；L_loss为供热系统的压强损失(包含供回水管路)；S_hl为供热系统的管路集；H_Vi为流体管道内的广延量流；R_Vi为管道流阻。

2) 等式约束。

综合能源系统包含电力系统和供热系统，因此也存在着相应的等式约束，如：电力系统的潮流平衡方程、管道的热力平衡方程、管路压力方程、供热温度和回水温度方程。由式(1)可知，稳态情况下式(1)左端为0，同时右端第4项源项为0，则可得支路的能量平衡方程。于是，模型中的等式约束数学表述如下：

(19) ${\begin{cases} P_{s} - A_{e} H_{e} X_{e} = 0 \\ ρ_{w} C_{p} H_{V} A_{h} (X_{s}^{h} - X_{0}^{h}) - Q_{h} = 0 \\ B_{h} R_{V} H_{V} = 0 \\ E_{s} (X_{s, load}^{h} - X_{0}^{h}) - b_{s} = 0 \\ E_{r} (X_{r, load}^{h} - X_{0}^{h}) - b_{r} = 0 \end{cases}$

式中：P_s为注入节点的功率；C_p为热水的比热容；ρ_w为流体的密度；Q_h为负荷节点需求的热功率；A_e、A_h为供电、供热网络的关联矩阵；B_h为供热系统流体网络的基本回路矩阵；H_e为电网络的广延量流矩阵，表征电流；X_e为电网络的强度矩阵，表征电压；R_V为管道流阻；H_V为热流体网的广延量流，表示体积流量； $X_{s,load}^{h} 、 X_{r,load}^{h}$ 分别为负荷节点的供水温度、返回温度； $X_{s}^{h} 、 X_{0}^{h}$ 分别为供水温度、环境温度；E_s、b_s、E_r和b_r分别为供热系统负荷节点的供水温度系数矩阵和返回温度系数矩阵。

3) 不等式约束。

对于该综合能源系统的优化模型，本文考虑控制变量有：发电机的端电压大小X_e，G、变压器分接头投切位置T_ap以及供热系统中供热点的热水供应温度X_h，supply，不等式约束表述如下：

(20) ${\begin{array}{l} X_{e, G i}^{\min} \leq X_{e, G i} \leq X_{e, G i}^{\max}, & i \in S_{G} \\ T_{ap}_{i}^{\min} \leq T_{ap i} \leq T_{ap}_{i}^{\max}, & i \in S_{Tap} \\ X_{h, supply, i}^{\min} \leq X_{h, supply, i} \leq X_{h, supply, i}^{\max}, & i \in S_{source} \end{array}$

式中：上标max、min分别表示该控制变量调节范围的上下限；S_G、S_Tap、S_source分别为所有发电机组集合、变压器分接头集合、供热点集合。

系统中的状态变量随控制变量变化而改变，包含电力系统各负荷节点电压幅值X_e、发电机产生的无功功率Q_G、管路的流量大小H_V以及供热系统中负荷节点热水供应温度X_h，s。不等式约束表述如下：

(21) ${\begin{array}{l} X_{e, i}^{\min} \leq X_{e, i} \leq X_{e, i}^{\max}, & i \in S_{P Q} \\ Q_{G i}^{\min} \leq Q_{G i} \leq Q_{G i}^{\max}, & i \in S_{G} \\ H_{V i}^{\min} \leq H_{V i} \leq H_{V i}^{\max}, & i \in S_{hl} \\ X_{h, s i}^{\min} \leq X_{h, s i} \leq X_{h, s i}^{\max}, & i \in S_{hload} \end{array}$

式中，各状态变量受相应的集合约束，S_PQ、S_G、S_hl、S_hload分别为PQ节点集合、发电机组集合、供热系统管路集合、热负荷集合。

2.2　粒子群算法

本文通过粒子群算法来求解某综合能源系统的优化模型，该算法源于生物捕食行为的分析，动物间通过集体的分工协作及信息共享使得种群能够快速寻优。假设模型中待求的问题解都可表示为粒子，对于这些粒子，在一个N维的空间内搜索，赋予每个粒子速度和位置的特征，通过自身信息以及和其他粒子的对比获得最终的调整方向。

粒子群算法流程如图1所示。

显示原图| 下载原图ZIP| 生成PPT

图1 粒子群算法流程图

Fig.1 Flow chart of particle swarm optimization

使用粒子群算法时，粒子群算法包括粒子的速度、位置变量以及适应度函数。在该运算场景中，综合能源系统的适应度值函数设定为总能量损失的倒数(总能量损失最小时，适应度值最大)，粒子的位置变量对应系统中的控制变量(包括变压器分接头位置、机组机端电压和热水供应温度)，粒子的速度变量对应变量的变化值，粒子位置可表示为(x₁，x₂，…，x_m)，粒子速度表示为(v₁，v₂，…，v_m)，m表示该综合能源系统具有m维变量。迭代过程中，每代更新当前粒子对应的控制变量信息，通过目标函数值与适应度值反映当前优化参数的优异性。

3　算例分析

本文针对为某岛屿提供电力和热的综合能源系统进行分析，该综合能源系统如图2所示^[16]，其中电力网络节点负荷、供热网络节点负荷分别如表1、表2所示。电力网络和供热网络通过3台热电联产机组、循环水泵等互联实现能量的耦合。

表1 电力网络节点负荷

Table 1 Load buses in power network

电力网络节点编号	1	2	3	4	5	6
负荷/MW	0.2	0.0	0.5	0.5	0.2	0.2

新窗口打开

表2 供热网络节点负荷

Table 2 Load buses in thermal network

供热网络节点编号	负荷/MW	供热网络节点编号	负荷/MW	供热网络节点编号	负荷/MW	供热网络节点编号	负荷/MW
1	0.000 0	9	0.107 0	17	0.080 5	25	0.000 0
2	0.000 0	10	0.107 0	18	0.080 5	26	0.107 0
3	0.107 0	11	0.145 0	19	0.000 0	27	0.107 0
4	0.145 0	12	0.107 0	20	0.080 5	28	0.000 0
5	0.000 0	13	0.000 0	21	0.080 5	29	0.107 0
6	0.107 0	14	0.080 5	22	0.000 0	30	0.107 0
7	0.107 0	15	0.000 0	23	0.107 0	31	0.000 0
8	0.107 0	16	0.080 5	24	0.107 0	32	0.000 0

新窗口打开

显示原图| 下载原图ZIP| 生成PPT

图2 某岛屿综合能源系统

Fig.2 An island integrated energy system

在图2中，机组G1、G2、G3分别为燃气轮机机组、抽汽式汽轮机机组和往复式内燃机机组，电力网络包含9条母线和5处负荷，供热网络包含32个节点和32条管路。算例的状态变量和控制变量情况分别如表3、表4所示。

表3 算例的状态变量情况

Table 3 State variables in case

状态变量	个数
节点电压相角	8
节点电压幅值	6
管道质量流量	32
热水供应温度	29
热水返回温度	29
总数目	104

新窗口打开

表4 算例的控制变量情况

Table 4 Control variables in case

控制变量	个数
变压器分接头位置	1
机组1—3机端电压	3
热水供应温度	3
总数目	7

新窗口打开

本例中粒子群算法的维数为7，设定粒子群算法学习因子c₁和c₂为2，惯性因子w为1，算法的迭代次数为60，粒子数为40。

1) 迭代曲线。

算例中，所选择的适应度函数为系统能量损失的倒数，通过描绘群体中每代最佳适应度值可以得到算法的收敛特性。图3为粒子群算法迭代60次的迭代曲线，可得知系统优化后的能量损失降至0.087 9 MW，其中热能损失为0.079 5 MW，电能有功损失为0.005 8 MW，压能损失为0.002 6 MW。

显示原图| 下载原图ZIP| 生成PPT

图3 算法迭代曲线

Fig.3 Iterative curve of algorithm

由图3可知：应用粒子群算法进行计算，每代粒子的适应度值如黑点所示；将黑点连接起来便形成算法收敛特性曲线，其反映算法的收敛趋势；该算例的适应度曲线在迭代20多次后达到最优，因此其具备较强的收敛能力和较快的收敛速度。

2) 结果对比。

通过粒子群算法可获得目标函数最优值和控制变量的优化值，在该最优粒子结果下，该综合能源系统的电压相角、电压幅值、管道质量流量及供应(返回)温度在优化前后的变化分别如图4，图5，图6所示。

显示原图| 下载原图ZIP| 生成PPT

图4 电压相角与幅值对比

Fig.4 Comparison of voltage angle and amplitude

显示原图| 下载原图ZIP| 生成PPT

图5 质量流量对比

Fig.5 Comparison of mass flows

显示原图| 下载原图ZIP| 生成PPT

图6 节点温度对比

Fig.6 Comparison of nodes' temperature

由图4可看出，通过粒子群算法计算所得优化后电力系统中各节点的电压相角值有所改善(向0趋近)，而电力网络的节点电压幅值接近于1 pu，说明粒子群算法优化的方向能够有效降低系统的有功损失，同时也增加了电力系统运行的稳定性与可靠性。

图5中质量流量的正负代表管道流量的流向，可见利用粒子群算法优化后部分供热系统管道内的质量流量较优化前有所增加，同时不影响管路内流体的流向，说明采用粒子群算法寻优能有效改善原管道内的流量分布，从而减小管道内的热能损失。

图6中，优化前后供热系统节点的供应温度有所下降(受限于用户侧温度变化)，返回温度值基本不变，使得在用户侧热负荷满足需求的前提下，能较大程度地降低热水流经管路中所造成的热量损失，从而减小了能量的总损失。

采用综合能源系统优化模型后，该系统的安全裕度得到有效增加，增强了系统的安全可靠性；在系统参数变化不大的情况下，原小岛综合能源系统的能量损失为0.095 1 MW，经优化后能量损失为0.087 9 MW，系统能量损失的下降比例达到7.571%，有效地提升了综合能源系统的能量利用效率。

4　结论

随着能源技术的高速发展，传统的供能网络面临的挑战也日益加深，IES概念的提出，使其结构与功能更是日趋复杂与异构。未来的供能领域，多能流耦合已经成为了IES的显著特征，为了优化IES的规划调度，得到其整体性能的准确运行模拟，必然要进行多能流系统的统一建模以及运行模拟。

本文从能量网络理论出发，基于能量公理化模型完成了某32节点综合能源系统的最优能量潮流计算。从结果来看，所建模型在减少IES运行过程的能量损耗方面具有显著的效果。然而，该模型在改善热力网络能量潮流分布的过程中，也一定程度加大了压能(即电能)的损耗。就能量品质及价值而言，电能远大于低品质的热能。这无疑有悖于IES的“能质匹配”原则。

后续将从能量品质()的角度出发，开展基于“能质匹配”原则的IES运行优化研究。

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]	刘振亚. 全球能源互联网[M]. 北京：中国电力出版社，2015. [本文引用: 1]
[2]	BHATTACHARYA M, PARAMATI S R, OZTURK I, et al. The effect of renewable energy consumption on economic growth: evidence from top 38 countries [J]. Applied Energy, 2016, 162: 733-741. [本文引用: 1]
[3]	孙宏斌，郭庆来，潘昭光. 能源互联网：理念、架构与前沿展望 [J]. 电力系统自动化，2015, 39(19): 1-8. [本文引用: 1] SUN Hongbin, GUO Qinglai, PAN Zhaoguang. Energy internet: concept, architecture and frontier outlook [J]. Automation of Electric Power Systems. 2015, 39(19): 1-8. [本文引用: 1]
[4]	徐宪东，贾宏杰，靳小龙，等. 区域综合能源系统电/气/热混合潮流算法研究 [J]. 中国电机工程学报，2015, 35(14): 3634-3642. [本文引用: 1] XU Xiandong, JIA Hongjie, JIN Xiaolong, et al. Study on hybrid heat-gas-power flow algorithm for integrated community energy system [J]. Proceedings of the CSEE, 2015, 35(14): 3634-3642. [本文引用: 1]
[5]	CORREA-POSADA C M, SANCHEZ-MARTIN P. Security-constrained optimal power and natural-gas flow [J]. IEEE Transactions on Power Systems, 2014, 29(4): 1780-1787. [本文引用: 1]
[6]	YANG J, ZHANG N, KANG C, et al. Effect of natural gas flow dynamics in robust generation scheduling under wind uncertainty [J]. IEEE Transactions on Power Systems, 2017, 33(2): 2087-2097. [本文引用: 1]
[7]	赵峰，张承慧，孙波，等. 冷热电联供系统的三级协同整体优化设计方法 [J]. 中国电机工程学报，2015, 35(15): 3785-3793. [本文引用: 1] ZHAO Feng, ZHANG Chenghui, SUN Bo, et al. Three-stage collaborative global optimization design method of combined cooling heating and power [J]. Proceedings of the CSEE. 2015, 35(15): 3785-3793. [本文引用: 1]
[8]	吴利兰，荆朝霞，吴青华，等. 基于Stackelberg博弈模型的综合能源系统均衡交互策略 [J]. 电力系统自动化，2018, 42(4): 142-150. [本文引用: 1] WU Lilan, JING Zhaoxia, WU Qinghua, et al. Equilibrium strategies for integrated energy systems based on stackelberg game model [J]. Automation of Electric Power Systems. 2018, 42(4): 142-150. [本文引用: 1]
[9]	陈皓勇，文俊中，王增煜，等. 能量网络的传递规律与网络方程 [J]. 西安交通大学学报，2014, 48(10): 66-76. [本文引用: 1] CHEN Haoyong, WEN Junzhong, WANG Zengyi, et al. Transfer laws and equations of energy networks [J]. Journal of Xi'an Jiaotong University. 2014, 48(10): 66-76. [本文引用: 1]
[10]	梅凤翔，刘桂林. 分析力学基础[M]. 西安：西安交通大学出版社，1987. [本文引用: 1]
[11]	陈群，郝俊红，陈磊，等. 电－热综合能源系统中能量的整体输运模型 [J]. 电力系统自动化，2017, 41(13): 7-13. [本文引用: 1] CHEN Qun, HAO Junhong, CHEN Lei, et al. Integral transport model for energy of electric-thermal integrated energy system [J]. Automation of Electric Power Systems. 2017, 41(13): 7-13. [本文引用: 1]
[12]	韩光泽，华贲，陈清林，等. 传递过程的统一描述 [J]. 工程热物理学报，2003, 24(6): 1025-1027. [本文引用: 1] HAN Guangze, HUA Ben, CHEN Qinglin, et al. A unified description of the transfer process [J]. Journal of Engineering Thermophysics. 2003, 24(6): 1025-1027. [本文引用: 1]
[13]	韩光泽. 寂态热动力学基础理论研究 [D]. 广州：华南理工大学，2002. [本文引用: 1] HAN Guangze. Research on the basic theory of dead state thermodynamics [D]. Guangzhou: South China University of Technology, 2002. [本文引用: 1]
[14]	罗志昌. 流体网络理论[M]. 北京：机械工业出版社，1988. [本文引用: 1]
[15]	德格鲁脱，梅休尔. 非平衡态热力学[M]. 上海：上海科学技术出版社，1981. [本文引用: 1]
[16]	LIU X, JENKINS N, WU J, et al. Combined analysis of electricity and heat networks [J]. Energy Procedia, 2014, 61: 155-159. [本文引用: 1]

2015

... 随着经济科技的迅速发展，人类社会对于能源的需求急剧上升.这直接导致了煤炭、石油等常规化石能源资源的迅速枯竭，并由此引发了全球范围内人类对于资源环境能否持续发展的担忧.世界各国开始探索一种能够提高能源的利用效率，适应可再生能源的开发利用并优化各类能源的统筹分配以实现可持续发展的新途径^[1,2].综合能源系统(integrated energy system, IES)以其能量转化设备(如热电联产机组、热泵及循环泵等)为纽带，通过供热网络、电力网络及供气网络等能源子网间的能量流动及互补，可有效实现能流的合理分配，并提升能量的利用效率，无疑将成为未来能源领域的发展方向^[3]. ...

The effect of renewable energy consumption on economic growth: evidence from top 38 countries

2016

能源互联网：理念、架构与前沿展望

2015

能源互联网：理念、架构与前沿展望

2015

区域综合能源系统电/气/热混合潮流算法研究

2015

... 国内外对于综合能源系统的研究，目前主要集中于IES耦合多能源系统的稳态潮流求解、IES优化规划方法及IES市场经济理论等几个方面.文献[4]基于能源集线器模型，构建了热电联产系统模型，针对区域综合能源系统(regional integrated energy system，RIES)的完全耦合、部分耦合及完全解耦3种运行模式，形成了适用于含电/气/热等3种能流的RIES稳态潮流求解方法.文献[5]建立了考虑安全约束的电－气最优潮流，该模型在求解最优潮流时考虑了系统发生N－1故障的情况，通过灵敏度分析可快速得到稳定的故障应对方案，避免线路过负荷.文献[6]考虑风力发电并网带来的不确定性以及天然气管道气压安全约束问题，建立了一个电－气耦合IES的二阶段鲁棒发电日前调度模型，通过分段线性化方法处理了燃气管网气流约束的非线性，提高了算法的收敛性.文献[7]针对IES的设备选型、设备容量配置及系统运行策略的整体优化，设计了一种冷热电联供系统的三级协同整体优化方法，并提出一种“以电定热”和“以热定电”混合运行策略，实现了系统产能以及用能的良好匹配.而在市场经济理论部分，吴利兰等人提出了一个基于多主多从Stackelberg博弈的能源交易模型，通过分析IES中多个分布式能源站和用户之间的多种能源的交互方式，求解它们之间的均衡交互策略，通过竞争决定能源价格，从而实现各自收益的最大化^[8]. ...

区域综合能源系统电/气/热混合潮流算法研究

2015

Security-constrained optimal power and natural-gas flow

2014

Effect of natural gas flow dynamics in robust generation scheduling under wind uncertainty

2017

冷热电联供系统的三级协同整体优化设计方法

2015

冷热电联供系统的三级协同整体优化设计方法

2015

基于Stackelberg博弈模型的综合能源系统均衡交互策略

2018

基于Stackelberg博弈模型的综合能源系统均衡交互策略

2018

能量网络的传递规律与网络方程

2014

... 作为能源互联网的关键组成部分，综合能源系统是典型的多能域耦合系统，是能量网络理论的研究对象.能量网络理论从能量角度出发，以能量作为不同物理过程之间相互耦合的桥梁，研究不同形式能量的传递、存储及转化的普适化规律，应用集总参数的分析方法，实现异质能流子网络的融合分析^[9]. ...

能量网络的传递规律与网络方程

2014

1987

... 式(1)中，等式左边的表达式是该圆柱形传递线(管)路中的能量因广延量流动产生的变化率；右边的第1项和第2项是流经管道横截面的广延量所携带的能量之和，第3项是传递过程中因过程不可逆出现的能量耗散，第4项是管道中自发生成的广延量产生过程中释放的能量.此外，由分析力学的基本知识可知，能量传递转换过程中还会产生广义动能(惯性现象)和势能(弹性现象)及其相互转化^[10].在稳态传输情况下，暂时忽略能量传递的弹性和惯性现象. ...

电－热综合能源系统中能量的整体输运模型

2017

... 虽然不同形式的能量传递过程从物质形态角度来看存在巨大差异，但从介质流的角度来看却呈现一定的规律性，如描述导电过程的欧姆定律、描述粘性流动的牛顿定律和描述热传导的傅里叶定律，其形式具有一致性^[11].传递公理指出，在相应基本强度量差的推动下，基本广延量实现流动运输^[12].广延量的传递规律可描述为 ...

电－热综合能源系统中能量的整体输运模型

2017

传递过程的统一描述

2003

传递过程的统一描述

2003

寂态热动力学基础理论研究

2002

... 压能传递过程中，强度量是压强p，基本广延量流是体积流率J_v.由质量守恒定律可知，质量的源强度g_m＝0.结合比体积V与密度ρ的关系式V＝ρ^－1，推导可得随体形式的体积平衡方程^[13]： ...

寂态热动力学基础理论研究

2002

1988

... 特别地，当基本广延量传递系数K为定值时(如稳态层流状态下的压能传递)，线(管)路的传阻可进一步简化为

R = \frac{L}{K A}

^[14]. ...

1981

... 热能传递过程中，以温度为强度量.同时由热力学第二定律可知，在一个与外界不存在相互作用的系统中，熵永不自动减少.在不可逆过程中，系统由不平衡状态过渡到平衡状态，熵单调增加，即熵的源强度g_s≠0^[15].此时以熵作为基础广延量将带来诸多不便. ...

Combined analysis of electricity and heat networks

2014

... 本文针对为某岛屿提供电力和热的综合能源系统进行分析，该综合能源系统如图2所示^[16]，其中电力网络节点负荷、供热网络节点负荷分别如表1、表2所示.电力网络和供热网络通过3台热电联产机组、循环水泵等互联实现能量的耦合. ...

〈

〉

基于能量网络理论的区域综合能源系统运行优化研究

Operation Optimization of Regional Integrated Energy System Based on Energy Network Theory

0 引言

1 IES网络等值模型

1.1 能量网络理论

1.2 传输管线等值模型

1.2.1 电能传递的能量网络等值模型

1.2.2 压能传递的能量网络等值模型

1.2.3 热能传递的能量网络等值模型

1.3 能量耦合元件模型

2 IES优化模型与求解

2.1 优化模型

2.2 粒子群算法

3 算例分析

4 结论

参考文献 文献选项 原文顺序 文献年度倒序 文中引用次数倒序 被引期刊影响因子

0　引言

1　IES网络等值模型

1.1　能量网络理论

1.2　传输管线等值模型

1.2.1　电能传递的能量网络等值模型

1.2.2　压能传递的能量网络等值模型

1.2.3　热能传递的能量网络等值模型

1.3　能量耦合元件模型

2　IES优化模型与求解

2.1　优化模型

2.2　粒子群算法

3　算例分析

4　结论

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子