留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

一种计及能量枢纽不同运行模式的综合能源系统混合能量流求解方法

潘益 王明深 叶昱媛 张宸宇 缪惠宇 郑建勇

潘益, 王明深, 叶昱媛, 张宸宇, 缪惠宇, 郑建勇. 一种计及能量枢纽不同运行模式的综合能源系统混合能量流求解方法[J]. 现代电力, 2021, 38(3): 277-285. doi: 10.19725/j.cnki.1007-2322.2021.0005
引用本文: 潘益, 王明深, 叶昱媛, 张宸宇, 缪惠宇, 郑建勇. 一种计及能量枢纽不同运行模式的综合能源系统混合能量流求解方法[J]. 现代电力, 2021, 38(3): 277-285. doi: 10.19725/j.cnki.1007-2322.2021.0005
Yi PAN, Mingshen WANG, Yuyuan YE, Chenyu ZHANG, Huiyu MIAO, Jianyong ZHENG. A Multi-energy Flow Calculation Method for Integrated Energy System Considering Different Operation Modes of Energy Hub[J]. Modern Electric Power, 2021, 38(3): 277-285. doi: 10.19725/j.cnki.1007-2322.2021.0005
Citation: Yi PAN, Mingshen WANG, Yuyuan YE, Chenyu ZHANG, Huiyu MIAO, Jianyong ZHENG. A Multi-energy Flow Calculation Method for Integrated Energy System Considering Different Operation Modes of Energy Hub[J]. Modern Electric Power, 2021, 38(3): 277-285. doi: 10.19725/j.cnki.1007-2322.2021.0005

一种计及能量枢纽不同运行模式的综合能源系统混合能量流求解方法

doi: 10.19725/j.cnki.1007-2322.2021.0005
详细信息
    作者简介:

    潘益(1993),男,博士,工程师,通信作者,研究方向:综合能源系统建模、优化,E-mail:panyi_seu@163.com

    王明深(1990),男,博士,工程师,研究方向:新能源、群体电动汽车充电规划、优化调控与运营机制,E-mail:wangmingshen@tju.edu.cn

    叶昱媛(1994),女,硕士,研究方向:故障诊断与评估算法、变电运维,E-mail:494666795@qq.com

    张宸宇(1989),男,博士,高级工程师,研究方向:交直流配用电技术,E-mail:njjsdky@163.com

    缪惠宇(1992),男,博士,工程师,研究方向:电能质量评估与控制,E-mail:2392341409@qq.com

    郑建勇(1966),男,教授,博士生导师,研究方向:智能电网与新能源发电,E-mail:zhengjy_seu@163.com

  • 中图分类号: TM73

A Multi-energy Flow Calculation Method for Integrated Energy System Considering Different Operation Modes of Energy Hub

  • 摘要: 综合能源系统的混合能量流求解是后续能量综合管理与调控的基础环节。为避免传统统一求解能量流时可能存在的模型维度高、收敛性差等问题,同时充分考虑能量枢纽(energy hub,EH)在不同运行模式下对系统能流分布的影响,提出了一种计及EH不同运行模式的综合能源系统混合能量流求解方法:基于EH建模以及建立电、气、热各能源子系统的能流独立求解模型,划分典型运行模式,利用EH单元输入输出的钳制约束对各能源子系统的能流独立计算结果进行迭代校正,实现EH在不同运行模式下的能流流程化求解。算例结果表明,此算例适用于不同运行模式的混合能量流计算,且具有较好的收敛性与合理性。
  • 图  1  第1类能量枢纽结构

    Figure  1.  Structure of the first type of energy hub

    图  2  第2类能源枢纽结构

    Figure  2.  Structure of the second type of energy hub

    图  3  综合能源系统多能流求解流程

    Figure  3.  Flowchart for multi-energy flow calculation of IES

    图  4  电-热-气综合能源系统算例结构

    Figure  4.  Case diagram of IES containing electricity, thermal and gas networks

    图  5  不同场景下的收敛曲线

    Figure  5.  Convergence curves in different cases

    B1  场景S1下关键元件N-1校验

    B1.  N-1 static security analysis of key components in scene S1

    B2  场景S2下关键元件N-1校验

    B2.  N-1 static security analysis of key components in scene S2

    表  1  电–热–气综合能源系统状态量信息

    Table  1.   Information of state variables for IES containing electricity, thermal and gas networks

    能量网络类型节点类型已知状态量未知状态量
    电力网络电力平衡节点节点电压幅值、相角节点注入有功无功功率
    PV节点节点注入有功功率、节点电压幅值节点注入无功功率、节点电压相角
    PQ节点节点注入有功无功功率节点电压幅值、相角
    热力网络热力平衡节点供水温度回水温度、供热功率、管道热媒流量
    热负荷节点回水温度、热功率供水温度、管道热媒流量
    其他热源节点供水温度、供热功率回水温度、管道热媒流量
    天然气网络天然气平衡节点节点压力节点注入气流量、天然气管道流量
    流量已知节点节点注入气流量节点压力、天然气管道流量
    下载: 导出CSV

    表  2  电力子系统潮流计算结果

    Table  2.   Calculation results for electricity flow

    支路编号支路有功功率/MW节点编号电压/pu
    S1S2S1S2
    d14.6153.479Bs11.0601.035
    d21.8071.128Bs21.0451.031
    d30.7630.508Bs31.0101.010
    d42.9112.478Bs41.0381.029
    d51.6561.207Bs51.0381.027
    d6−0.421−0.570Bs61.0701.070
    d7−5.319−5.360Bs71.0791.072
    d81.9021.560Bs81.0781.072
    d90.9450.746Bs91.0881.083
    d10−2.679−3.848Bs101.0841.080
    d111.3921.725Bs111.0771.075
    d120.7390.780Bs121.0701.070
    d132.0702.242Bs131.0711.070
    d140.8000.800Bs141.0791.076
    d151.1020.760
    d16−0.101−0.432
    d170.6740.463
    d18−1.009−1.340
    d190.1280.169
    d200.8431.056
    下载: 导出CSV

    表  3  电力子系统各发电单元出力

    Table  3.   Output of each electricity generation unit

    发电单元有功出力/MW
    S1S2
    EH15.98195.233
    EH26.48178
    G133
    G286.706
    下载: 导出CSV

    表  4  热力子潮流计算结果

    Table  4.   Calculation results for thermal flow

    管道编号管道流量/(kg/s)节点编号供水网络节点输出温度/℃
    S1S2S1S2
    l19.1739.174Hs189.97689.987
    l29.1749.174Hs289.90089.978
    l39.1939.178Hs389.88789.956
    l49.1929.177Hs489.88489.985
    l59.1949.184Hs589.99689.994
    l69.1939.182lh189.95489.952
    l77.5947.584lh289.94989.948
    l87.5947.584lh389.82489.925
    l97.5847.583lh489.82989.931
    l107.5847.583lh589.81689.885
    l114.86122.427lh689.82889.897
    l12−10.3277.259lh789.82789.905
    l13−28.714−11.107lh889.82889.906
    l14−47.100−29.462lh989.90389.914
    l15−65.447−47.810lh1089.90389.914
    l16−20.030−37.593
    下载: 导出CSV

    表  5  天然气子系统潮流计算结果

    Table  5.   Calculation results of natural gas flow of natural gas subsystem

    管道编号管道气流量/MMCFD节点编号节点压力/psia
    S1S2S1S2
    gl1201.448200.375Gs11000.0001000.000
    gl2196.281193.406Gs2643.607608.896
    gl349.42136.356Gs3627.115596.632
    gl4118.502110.748Gs4498.023476.758
    gl5145.009−190.645Gs5647.429619.786
    gl677.456−86.392Gs6384.860909.544
    gl7−109.080155.356Gs7538.8041273.362
    gl8117.976117.855Gs8485.013719.426
    gl9123.748123.653Gs91000.0001000.000
    gl10−33.631−33.627Gs10883.147864.692
    gl1141.60741.482Gs111100.0001100.000
    gl1250.11750.026Gs12691.885692.655
    Gs13723.185723.913
    Gs14645.180646.296
    下载: 导出CSV

    表  6  天然气子系统压缩机耗气量计算结果

    Table  6.   Gas consumption calculation results of the compressors of natural gas subsystem

    压缩机编号压缩机耗气量/MMCFD压缩机编号压缩机耗气量/MMCFD
    S1S2S1S2
    10.3410.31930.6570.315
    20.5330.70140.4660.405
    下载: 导出CSV
  • [1] CHEN Sheng, WEI Zhinong, SUN Guoqiang, et al. Adaptive robust day-ahead dispatch for urban energy systems[J]. IEEE Transactions on Industrial Electronics, 2019, 66(2): 1379−90. doi:  10.1109/TIE.2017.2787605
    [2] OMALLEY M, KROPOSKI B. Energy comes together: the integration of all systems[J]. IEEE Power & Energy Magazine, 2013, 11(5): 18−23.
    [3] 周晟锐, 刘继春, 张浩禹, 等. 基于机会约束的多能源枢纽电气互联综合能源系统日前经济调度[J]. 现代电力, 2020, 37(2): 197−204.

    ZHOU Shengrui, LIU Jichun, ZHANG Haoyu, et al. Day-ahead economic dispatch for electricity-gas integrated energy system with multiple energy hubs based on chance constraints[J]. Modern Electric Power, 2020, 37(2): 197−204(in Chinese).
    [4] BARATI F, SEIFI H, SEPASIAN M S, et al. Multi-period integrated framework of generation, transmission, and natural gas grid expansion planning for largescale systems[J]. IEEE Transactions on Power Systems, 2015, 30(5): 2527−2537. doi:  10.1109/TPWRS.2014.2365705
    [5] CLEGG S, MANCARELLA P. Storing renewables in the gas network: modelling of power-to-gas seasonal storage flexibility in low carbon power systems[J]. IET Generation, Transmission & Distribution, 2016, 10(3): 566−575.
    [6] 黎静华, 黄玉金, 张鹏. 综合能源系统多能流潮流计算模型与方法综述[J]. 电力建设, 2018, 39(3): 1−11. doi:  10.3969/j.issn.1000-7229.2018.03.001

    Li Jinghua, HUANG Yujin, ZHANG Peng. Review of multi-energy flow calculation model and method in integrated energy system[J]. Electric Power Construction, 2018, 39(3): 1−11(in Chinese). doi:  10.3969/j.issn.1000-7229.2018.03.001
    [7] 林威, 靳小龙, 穆云飞, 等. 区域综合能源系统多目标最优混合潮流算法[J]. 中国电机工程学报, 2017, 37(20): 5829−5839.

    LIN Wei, JIN Xiaolong, MU Yunfei, et al. Multi-objective optimal hybrid power flow algorithm for integrated local area energy system[J]. Proceeding of the CSEE, 2017, 37(20): 5829−5839(in Chinese).
    [8] MARTINE-MARES A, Fuerte-Esquivel CR. A unified gas and power flow analysis in natural gas and electricity coupled networks[J]. IEEE Transactions on Power Systems, 2012, 27(4): 2156−2166. doi:  10.1109/TPWRS.2012.2191984
    [9] CHEN Sheng, SUN Guoqiang, WEI Zhinong, et al. Multi-time combined gas and electric system optimal power flow incorporating wind power[C]//3rd International Conference on Power and Energy Systems Engineering. Kitakyushu, JAPAN, 2016: 111-116.
    [10] Liu X Z, Wu J Z, NICK J, et al. Combined analysis of electricity and heat networks[J]. Applied Energy, 2016(162): 1238−1250.
    [11] 王英瑞, 曾博, 郭经, 等. 电–热–气综合能源系统多能流计算方法[J]. 电网技术, 2016, 40(10): 2942−2951.

    WANG Yingrui, ZENG Bo, GUO Jing, et al. Multi-energy flow calculation method for integrated energy system containing electricity, heat and gas[J]. Power System Technology, 2016, 40(10): 2942−2951(in Chinese).
    [12] 徐宪东, 贾宏杰, 靳小龙, 等. 区域综合能源系统电/气/热混合潮流算法研究[J]. 中国电机工程学报, 2015, 35(14): 3634−3642.

    XU Xiandong, JIA Hongjie, JIN Xiaolong, et al. Study on hybrid heat-gas-power flow algorithm for integrated community energy system[J]. Proceeding of the CSEE, 2015, 35(14): 3634−3642(in Chinese).
    [13] 刘述欣, 戴赛, 胡林献, 等. 计及回水管网热损失的电热联合系统潮流模型及算法[J]. 电力系统自动化, 2018, 42(4): 77−81+134. doi:  10.7500/AEPS20170508013

    LIU Shuxin, DAI Sai, HU Linxian, et al. Power flow model and algorithm of combined power and heat system considering heat loss in return pipe network[J]. Automation of Electric Power Systems, 2018, 42(4): 77−81+134(in Chinese). doi:  10.7500/AEPS20170508013
    [14] 赵霞, 杨仑, 瞿小斌, 等. 电-气综合能源系统能流计算的改进方法[J]. 电工技术学报, 2018, 33(3): 467−477.

    ZHAO Xia, YANG Lun, QU Xiaobin, et al. An improved energy flow calculation method for integrated electricity and natural gas system[J]. Transactions of China Electrotechnical Society, 2018, 33(3): 467−477(in Chinese).
    [15] 王伟亮, 王丹, 贾宏杰, 等. 考虑天然气网状态的电力–天然气区域综合能源系统稳态分析[J]. 中国电机工程学报, 2017, 37(5): 1293−1305.

    WANG Weiliang, WANG Dan, JIA Hongjie, et al. Steady state analysis of electricity-gas regional integrated energy system with consideration of NGS network status[J]. Proceeding of the CSEE, 2017, 37(5): 1293−1305(in Chinese).
    [16] 张刚, 张峰, 张利, 等. 考虑多种耦合单元的电气热联合系统潮流分布式计算方法[J]. 中国电机工程学报, 2018, 38(22): 6594−6605.

    ZHANG Gang, ZHANG Feng, ZHANG Li, et al. Distributed algorithm for the power flow calculation of integrated electrical, gas, and heating network considering various coupling units[J]. Proceeding of the CSEE, 2018, 38(22): 6594−6605(in Chinese).
    [17] 潘益, 梅飞, 郑建勇, 等. 计及静态安全因素与热电最优潮流的综合能源系统联合运行优化模型[J]. 电网技术, 2019, 43(1): 50−59.

    PAN Yi, MEI Fei, ZHENG Jianyong, et al. Operation optimization model for multi-integrated energy systems considering static security and optimal energy flow[J]. Power system Technology, 2019, 43(1): 50−59(in Chinese).
    [18] 江茂泽. 输配气管网的模拟与分析[M]. 北京: 石油工业出版社, 1995.
    [19] 艾芊, 郝然. 多能互补、集成优化能源系统关键技术及挑战[J]. 电力系统自动化, 2018, 42(4): 2−10+46. doi:  10.7500/AEPS20170927008

    AI Qian, HAO Ran. Key technologies and challenges for multi-energy complementarity and optimization of integrated energy system[J]. Automation of Electric Power Systems, 2018, 42(4): 2−10+46(in Chinese). doi:  10.7500/AEPS20170927008
    [20] SHARIF A, ALMANSOONRI A, FOWLER M, et al. Design of an energy hub based on natural gas and renewable energy sources[J]. International Journal of Energy Research, 2014, 38(3): 363−373. doi:  10.1002/er.3050
    [21] 包铭磊, 杨阳, 丁一, 等. 考虑天然气系统影响的电力系统连锁故障评估[J]. 电网技术, 2019, 43(1): 32−40.

    BAO Minglei, YANG Yang, DING Yi, et al. Assessment of cascading failures in power system considering the effects of natural gas system[J]. Power System Technology, 2019, 43(1): 32−40(in Chinese).
  • [1] 李瑶虹, 黄伟, 孙贝贝, 张新鹤, 薛一鸣.  拒绝服务攻击下基于分布式事件触发一致性预测补偿的微电网能量优化管理 . 现代电力, 2021, 38(2): 178-186. doi: 10.19725/j.cnki.1007-2322.2020.0233
    [2] 周荔丹, 曹祖加, 姚钢, 刘东.  泛在电力物联网的发展分析 . 现代电力, 2021, 38(2): 119-128. doi: 10.19725/j.cnki.1007-2322.2020.0326
    [3] 王晨, 徐光福, 侯炜, 陈俊, 朱皓斌, 王健, 高姗姗.  弱电网下变电站级微网运行模式平滑切换控制策略 . 现代电力, 2021, 38(3): 258-267. doi: 10.19725/j.cnki.1007-2322.2020.0309
    [4] 孟明, 马辰南, 薛宛辰, 罗洋, 商聪.  基于多代理的综合能源系统分层分布式能量协调方法 . 现代电力, 2021, 38(2): 129-137. doi: 10.19725/j.cnki.1007-2322.2020.0256
    [5] 周晟锐, 刘继春, 张浩禹, 文杰, 张林, 冯麒铭.  基于机会约束的多能源枢纽电气互联综合能源系统日前经济调度 . 现代电力, 2020, 37(2): 197-204. doi: 10.19725/j.cnki.1007-2322.2019.0966
    [6] 杨秀媛, 刘凤鸣, 陈麒宇, 周博文, 卜思齐, 徐智蔷.  利益驱动的泛在电力物联网 . 现代电力, 2020, 37(1): 1-9. doi: 10.19725/j.cnki.1007-2322.2019.1063
    [7] 黄伟, 杨子力, 柳思岐.  基于物元可拓模型的特色小镇能源系统综合评价 . 现代电力, 2020, 37(5): 448-455. doi: 10.19725/j.cnki.1007-2322.2019.0620
    [8] 祁兵, 何承瑜, 李彬, 陈宋宋, 石坤, 薛溟枫.  基于蓄热电锅炉不同工作模式的区域综合能源系统优化调度 . 现代电力, 2019, 36(6): 45-51.
    [9] 李子林, 刘蓉晖.  考虑需求侧响应的含储能区域综合能源系统运行优化 . 现代电力, 2019, 36(6): 61-67.
    [10] 周钰童, 华亮亮, 黄伟, 刘明昌.  计及电热交易的区域综合能源多目标优化配置 . 现代电力, 2019, 36(4): 24-30.
    [11] 蔡佳铭, 王承民, 谢 宁, 彭 石.  基于能量函数模型的风电功率在线预测方法研究 . 现代电力, 2018, 35(2): 71-79.
    [12] 李德智, 董明宇, 尹忠东, 李晓强, 林挚.  一种基于模糊控制的智能家庭能量管理系统 . 现代电力, 2018, 35(6): 39-46.
    [13] 陈丽雪, 房方.  计及蓄电池寿命的风光储微网系统能量优化管理 . 现代电力, 2018, 35(3): 62-69.
    [14] 姬 源, 黄育松.  基于负荷聚集商的多家庭能量管理系统建模 . 现代电力, 2018, 35(2): 16-21.
    [15] 刘文霞, 李校莹, 刘座铭, 王 皑, 王佳伟.  基于主从博弈的能源服务商运营模式与能量管理 . 现代电力, 2018, 35(2): 8-15.
    [16] 林岩, 张建成, 王宁.  一种基于虚拟同步发电机的光储发电系统能量管理控制策略 . 现代电力, 2017, 34(5): 64-73.
    [17] 刘方, 杨秀.  含电蓄热解决高渗透率风电消纳问题的微网热电能量综合优化 . 现代电力, 2015, 32(4): 19-26.
    [18] 刘永前, 王一妹, 韩爽, 李莉.  基于CFD流场预计算的风电功率预测误差分布研究 . 现代电力, 2014, 31(5): 64-69.
    [19] 任思源, 何 青.  基于压电材料的振动能量收集试验研究 . 现代电力, 2010, 27(3): 70-74.
    [20] 隋丽颖, 黄孝彬, 谭 文, 田志强, 李 娜.  基于数据挖掘的火电机组运行模式建立方法研究 . 现代电力, 2010, 27(2): 74-77.
  • 加载中
图(7) / 表 (6)
计量
  • 文章访问数:  14
  • HTML全文浏览量:  1
  • PDF下载量:  3
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-01-06
  • 刊出日期:  2021-06-03

一种计及能量枢纽不同运行模式的综合能源系统混合能量流求解方法

doi: 10.19725/j.cnki.1007-2322.2021.0005
    作者简介:

    潘益(1993),男,博士,工程师,通信作者,研究方向:综合能源系统建模、优化,E-mail:panyi_seu@163.com

    王明深(1990),男,博士,工程师,研究方向:新能源、群体电动汽车充电规划、优化调控与运营机制,E-mail:wangmingshen@tju.edu.cn

    叶昱媛(1994),女,硕士,研究方向:故障诊断与评估算法、变电运维,E-mail:494666795@qq.com

    张宸宇(1989),男,博士,高级工程师,研究方向:交直流配用电技术,E-mail:njjsdky@163.com

    缪惠宇(1992),男,博士,工程师,研究方向:电能质量评估与控制,E-mail:2392341409@qq.com

    郑建勇(1966),男,教授,博士生导师,研究方向:智能电网与新能源发电,E-mail:zhengjy_seu@163.com

  • 中图分类号: TM73

摘要: 综合能源系统的混合能量流求解是后续能量综合管理与调控的基础环节。为避免传统统一求解能量流时可能存在的模型维度高、收敛性差等问题,同时充分考虑能量枢纽(energy hub,EH)在不同运行模式下对系统能流分布的影响,提出了一种计及EH不同运行模式的综合能源系统混合能量流求解方法:基于EH建模以及建立电、气、热各能源子系统的能流独立求解模型,划分典型运行模式,利用EH单元输入输出的钳制约束对各能源子系统的能流独立计算结果进行迭代校正,实现EH在不同运行模式下的能流流程化求解。算例结果表明,此算例适用于不同运行模式的混合能量流计算,且具有较好的收敛性与合理性。

English Abstract

潘益, 王明深, 叶昱媛, 张宸宇, 缪惠宇, 郑建勇. 一种计及能量枢纽不同运行模式的综合能源系统混合能量流求解方法[J]. 现代电力, 2021, 38(3): 277-285. doi: 10.19725/j.cnki.1007-2322.2021.0005
引用本文: 潘益, 王明深, 叶昱媛, 张宸宇, 缪惠宇, 郑建勇. 一种计及能量枢纽不同运行模式的综合能源系统混合能量流求解方法[J]. 现代电力, 2021, 38(3): 277-285. doi: 10.19725/j.cnki.1007-2322.2021.0005
Yi PAN, Mingshen WANG, Yuyuan YE, Chenyu ZHANG, Huiyu MIAO, Jianyong ZHENG. A Multi-energy Flow Calculation Method for Integrated Energy System Considering Different Operation Modes of Energy Hub[J]. Modern Electric Power, 2021, 38(3): 277-285. doi: 10.19725/j.cnki.1007-2322.2021.0005
Citation: Yi PAN, Mingshen WANG, Yuyuan YE, Chenyu ZHANG, Huiyu MIAO, Jianyong ZHENG. A Multi-energy Flow Calculation Method for Integrated Energy System Considering Different Operation Modes of Energy Hub[J]. Modern Electric Power, 2021, 38(3): 277-285. doi: 10.19725/j.cnki.1007-2322.2021.0005
    • 在能源系统转型与能源互联网的大背景下,综合能源系统(integrated energy system,IES)作为一种以电能为核心,协调优化多种异质能源形式、提供优质能源服务的物理载体,对提高能源综合利用率,促进可再生能源的消纳,实现多能互补以及能量的阶梯利用具有重要意义[1-3]

      综合能源系统多能流建模与能流求解是综合能源系统构建、运行管理的前提与基础。基于系统建模,开展混合潮流计算是后续优化调度与安全分析的核心。目前单一供能网络的能量流计算已有较为成熟的模型,但对于多种供能网络耦合的联合潮流计算,其研究还有待深入[4-5]。对于现有计算方法的研究,主要可以分为统一求解与解耦求解2大类[6-9]。统一求解法多通过建立多能源系统的矩阵方程并构造雅克比矩阵,结合牛顿–拉夫逊算法进行联合求解;解耦求解对各能源子系统进行独立求解,并通过能源耦合元件实现变量的迭代更新至误差在允许范围内。文献[10]针对热电耦合网络,同时利用这2类方法对热电潮流进行求解并进行了对比分析;文献[11]提出了一种适用于电/热/气综合能源系统的扩展牛–拉法多能流计算方法,实现多能流的统一计算。多能流的统一求解具有较高的普适性,较为直观,但模型较为复杂,对求解算法的要求较高,收敛性无法保证。文献[12]针对区域综合能源系统各能源子系统的不同耦合程度,提出了一种基于顺序求解思路的电力–天然气混合潮流算法,将天然气与电力网络分开求解,简化了计算的复杂度;文献[13]针对辐射型供热网,提出了热网潮流的前推回代计算方法并利用牛–拉法计算电网潮流,实现热电联合系统的解耦计算。

      目前,上述的多能流求解方法多针对特定的拓扑进行构建,而不同的能量枢纽(energy hub,EH)运行方式将对整个能源网络的拓扑参数产生影响,进而改变混合能量流的分布情况,因此这些多能流求解方法未必兼顾不同运行模式,使得普适性难以得到保障。

      本文针对电、气、热耦合型综合能源系统,构建不同能源子系统能量流的独立求解模型及方法,并提出以分布式顺序求解算法为核心的混合能量流计算流程,依托EH单元输入输出的钳制约束对各能源子系统的能流进行迭代校正,实现EH在不同运行模式下的能流流程化求解。

    • 对于电、热、气各能源子网络的独立建模,目前已有较为丰富的研究,具体建模方法可参考文献[14-16],本文不再赘述。

      综合能源系统中的电力子系统潮流采用传统牛–拉法进行求解。对于结构简单的辐射状热力网,可采用前推回代法[13]进行温度、流量的迭代计算。对于较复杂的多源环状热网,为提高求解效率,考虑以温度–流量模型[17]为基础,改写为如式(1)描述的热网状态矩阵方程:

      $$\Delta { F}{\rm{ = }}\Delta \left[ {\begin{array}{*{20}{c}} {\Delta {{H}}} \\ {\Delta { p}} \\ {\Delta {{ T}_{\rm{s}}}} \\ {\Delta {{ T}_{\rm{r}}}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{{ C}_{\rm{p}}}{ M}\left( {{{ T}_{\rm{s}}}-{{ T}_{\rm{r}}}} \right)-{{{H}}^{\rm{SP}}}} \\ {{{BKm}}\left| { m} \right|} \\ {{{ C}_{\rm{s}}}{{ T}_{\rm{s}}}-{{ b}_{\rm{s}}}} \\ {{{ C}_{\rm{r}}}{{ T}_{\rm{r}}}-{{ b}_{\rm{r}}}} \end{array}} \right] = 0$$ (1)

      式中:HSP为系统已知的热功率信息,当HSP对应EH与热网交互的热功率,M=mq=−Am;当HSP对应节点热功率,M=mmmq分别表示各管道支路流量以及流经热源或负荷节点的流量;K为管道阻力系数矩阵[10]B表示支路–回路关联矩阵;CsbsCrbr分别表示供热、回水网络拓扑、热媒流量与节点温度有关的矩阵,具体构造方法如文献[11]所述。TsTr分别包括热源EH侧供回水温度以及热媒汇聚节点供回水温度2部分。

      根据状态矩阵方程,利用牛–拉法,构造雅克比矩阵对式(1)进行迭代求解,其中待求状态量包括xh=[Ts,Tr,m]。

      天然气网络能流计算时同样考虑先将系统潮流方程组改写为式(2)所示的矩阵形式,进而构造天然气系统潮流矩阵方程:

      $$\left\{ \begin{aligned} & {{ A}_{_1}}{{ Q}_{\rm{l}}}{\rm{ - }}{ L}{\rm{ = }}{{ A}_{_1}}\left[ {{{\varPhi }^{'}}\left( {{\rm{ - }}{{ A}^{\rm{T}}}{{ P}_{\rm{f}}}} \right)} \right] - { L} = 0 \\ & {{ B}_1}{\varPhi }\left( {{{ Q}_{\rm{l}}}} \right){\rm{ = }}{{ B}_1}\Delta {{ P}_{\rm{l}}}{\rm{ = }}0 \end{aligned} \right.$$ (2)

      式中:A1B1分别表示为简约的节点-管路关联矩阵与管路–节点关联矩阵;AT为节点–管路关联矩阵的转置;Pf与ΔPl分别代表节点压力降矩阵与管道压力降矩阵;QlL分别代表管道流量矩阵以及节点气负荷矩阵;Φ(•)与Φ'(•)分别代表流量函数与压力降函数,两者互为反函数,天然气系统潮流矩阵方程的具体构造流程可详见参考文献[18]。基于天然气系统潮流矩阵方程,同样利用牛–拉法,构造雅克比矩阵对式(2)进行迭代求解,其中待求状态量包括xg=[Pf, Ql]。

    • 各供能网络间的耦合与多能形态的转化通常由IES的EH完成。目前通常利用能源集线器模型[19-20]对EH单元进行描述。

      EH内部包含较多的能量转换与存储环节,不同的能源元件配置与组合形式将对EH的外部输出特性产生影响,本节选择2种典型的结构进行分析说明。

      图1所示为第1类EH结构,主要分布在负荷侧,输入端连接外部多能流网络负荷节点,输出端直接向用户侧供能,外部能源网络输入的各种形式能量经EH单元转换为满足用户多元需求的电、热、冷、气负荷,此类EH的输入输出耦合关系可表示为

      $$\left[ {\begin{array}{*{20}{c}} {{L_{\rm{e}}}}\\ {}\\ {{L_{\rm{h}}}}\\ {}\\ {{L_{\rm{c}}}}\\ {}\\ {{L_{\rm{g}}}}\\ {} \end{array}} \right]{\rm{ = }}\left[ {\begin{array}{*{20}{c}} {{\eta _{\rm{T}}}(1 - {\omega _1} - {\omega _2})}&1&{\sigma \eta _{{\rm{chp}}}^{\rm{e}}(1 - {\omega _1} - {\omega _2})}&0\\ {{\eta _{\rm{T}}}{\omega _1}{\eta _{{\rm{EB}}}}\theta }&{{\omega _1}{\eta _{{\rm{EB}}}}}&{\begin{array}{*{20}{l}} {\theta \sigma (\eta _{{\rm{chp}}}^{\rm{h}} + }\\ {\eta _{{\rm{chp}}}^{\rm{e}}{\eta _{{\rm{EB}}}}{\omega _1})} \end{array}}&{\theta {\eta _{{\rm{ex}}}}}\\ {\begin{array}{*{20}{l}} {(1 - \theta ){\omega _1}{\eta _{{\rm{AC}}}}{\eta _{\rm{T}}}}\\ { + {\omega _2}{\eta _{\rm{T}}}{\eta _{{\rm{EC}}}}} \end{array}}&{\begin{array}{*{20}{l}} {(1 - \theta ){\omega _1}{\eta _{{\rm{EB}}}}{\eta _{{\rm{AC}}}}}\\ { + {\omega _2}{\eta _{{\rm{EC}}}}} \end{array}}&{\begin{array}{*{20}{l}} {(1 - \theta ){\eta _{\rm{AC}}}(\sigma \eta _{{\rm{chp}}}^h}\\ { + \sigma {\omega _1}{\eta _{{\rm{EB}}}}\eta _{{\rm{chp}}}^{\rm{e}}) + }\\ {\sigma {\omega _2}\eta _{{\rm{chp}}}^{\rm{e}}{\eta _{{\rm{EC}}}}} \end{array}}&{(1 - \theta ){\eta _{{\rm{ex}}}}{\eta _{{\rm{AC}}}}}\\ 0&0&{1 - \sigma }&0 \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {\begin{array}{*{20}{l}} {P_{\rm{e}}^{{\rm{net}}}}\\ {}\\ {{P_{{\rm{dg}}}}} \end{array}}\\ {}\\ {P_{\rm{g}}^{{\rm{net}}}}\\ {\begin{array}{*{20}{l}} {}\\ {P_{\rm{h}}^{{\rm{net}}}} \end{array}} \end{array}} \right]{\rm{ + }}\left[ {\begin{array}{*{20}{c}} {\begin{array}{*{20}{l}} {P_{\rm{e}}^s}\\ {}\\ {P_{\rm{h}}^{\rm{s}}} \end{array}}\\ {}\\ 0\\ {\begin{array}{*{20}{l}} {}\\ {P_{\rm{g}}^{\rm{s}}} \end{array}} \end{array}} \right]$$ (3)

      图2所示为第2类EH结构,主要分布在源测,为能源网络末端的多元负荷供能。输入端主要连接天然气站负荷节点以及各类分布式能源,输出端主要连接能源网络的源节点,输出的冷热电气能通过各类输能网络输送至用户侧,此类EH的输入输出耦合关系可表示为

      $$\left[ {\begin{array}{*{20}{c}} \begin{array}{l} {P_{\rm{o,e}}}\\ \\ {P_{\rm{o,h}}} \end{array}\\ {}\\ {{P_{\rm{o,c}}}}\\ \begin{array}{l} \\ {P_{\rm{o,g}}}\\ \end{array} \end{array}} \right]{\rm{ = }}\left[ {\begin{array}{*{20}{c}} {{\eta _{\rm{T}}}(1 - {\omega _1} - {\omega _2})}&{{\eta _{\rm{T}}}(1 - {\omega _1} - {\omega _2})}&{{\sigma _1}\eta _{{\rm{chp}}}^e{\eta _{\rm{T}}}(1 - {\omega _1} - {\omega _2})}\\ {{\omega _1}\theta \eta _{\rm{hp}}^h{\eta _{ex}}}&{{\omega _1}\theta \eta _{\rm{hp}}^h{\eta _{ex}}}&\begin{array}{l} ({\sigma _2}{\eta _{\rm{GB}}} + {\sigma _1}\eta _{\rm{chp}}^{\rm{h}} + \\ {\omega _1}{\sigma _1}\eta _{\rm{chp}}^{\rm{e}}\eta _{\rm{hp}}^{\rm{h}})\theta {\eta _{\rm{ex}}} \end{array}\\ {(1{\rm{ - }}\theta ){\omega _1}\eta _{\rm{hp}}^{\rm{h}}{\eta _{\rm{ex}}}}&{(1{\rm{ - }}\theta ){\omega _1}\eta _{\rm{hp}}^{\rm{h}}{\eta _{\rm{ex}}}}&\begin{array}{l} ({\sigma _2}{\eta _{\rm{GB}}} + {\sigma _1}\eta _{\rm{chp}}^{\rm{h}} + \\ {\omega _1}{\sigma _1}\eta _{\rm{chp}}^{\rm{e}}\eta _{\rm{hp}}^{\rm{h}}){\eta _{\rm{ex}}}(1 - \theta ) \end{array}\\ {{\omega _2}{\eta _{\rm{P2G}}}}&{{\omega _2}{\eta _{\rm{P2G}}}}&\begin{array}{l} (1 - {\sigma _1} - {\sigma _2}){\rm{ + }}\\ {\omega _2}{\sigma _1}\eta _{\rm{chp}}^{\rm{e}}{\eta _{\rm{P2G}}} \end{array} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} \begin{array}{l} P_{\rm{pv}}^{\rm{dg}}\\ \\ \end{array}\\ {P_{\rm{wt}}^{\rm{dg}}}\\ {}\\ \begin{array}{l} P_{\rm{g}}^{\rm{net}}\\ \end{array} \end{array}} \right] + \left[ {\begin{array}{*{20}{c}} {P_{\rm{e}}^{\rm{s}}}\\ {}\\ {{\eta _{\rm{ex}}}P_{\rm{h}}^{\rm{s}}}\\ {}\\ {P_{\rm{c}}^{\rm{s}}}\\ {}\\ {P_{\rm{g}}^{\rm{s}}} \end{array}} \right]$$ (4)

      式中:LeLhLcLg分别为电、热/冷、气负荷;ω1ω2为电能分配系数;θ为热能分配系数;σ为天然气分配系数;ηT$\eta _{{\rm{chp}}}^{\rm{e}}$$\eta _{{\rm{chp}}}^{\rm{h}}$ηEBηECηACηex分别表示变压器效率、热电联产(combined heat and power,CHP)机组的电效率、热效率、电锅炉制热效率、电制冷机与吸收式制冷机的制冷效率以及热交换器效率;$\eta _{{\rm{hp}}}^{\rm{h}}$ηGB分别表示热泵与燃气锅炉的制热效率;$P_{{\rm{pv}}}^{{\rm{dg}}}$$ P_{{\rm{wt}}}^{{\rm{dg}}} $分别代表分布式光伏与风机的发电功率;ηP2G为P2G设备的转换效率;$ P_{\rm{e}}^{\rm{s}} $$ P_{\rm{h}}^{\rm{s}} $$ P_{\rm{c}}^{\rm{s}} $以及$ P_{\rm{g}}^{\rm{s}}$分别代表电、热、冷、气储能单元功率。

    • 根据EH在能源传输网络中的作用,可从能源网络层面将其划分为4种运行模式。

      图  1  第1类能量枢纽结构

      Figure 1.  Structure of the first type of energy hub

      图  2  第2类能源枢纽结构

      Figure 2.  Structure of the second type of energy hub

      模式1:此模式下各EH节点对于电–热–气输能网络而言均为负荷节点,综合能源网络的平衡节点分别由外部电源(上级电网)、热源(上级换热站)及气源(上级燃气站)承担,各EH从各输能网络获得能量后经转换直接向所接负荷供能,其输入输出约束关系同式(3),在多能流稳态计算时各能源网络可独立进行求解。

      模式2:此模式下选择某一具有充足容量、稳定热能输出的EH作为热力网平衡节点,即热网侧的热力流平衡由此EH实现,其余EH视为PV/PQ节点、普通热源节点及天然气负荷节点,EH的输入输出约束关系同式(4)。此时,各EH所接入的电力网侧潮流由外部电网平衡,因此各EH单元的功率输出波动仅对平衡节点处的EH单元出力产生影响,各能源网络处于弱耦合状态。

      模式3:此模式下选择某一具有充足容量、稳定、动态响应快速的EH作为电力网络侧的平衡节点,其余EH节点均视为PV/PQ节点、普通热源节点及天然气负荷节点。此时,各能源网络的耦合状态与模式2类似。

      模式4:此模式下选择两个EH分别作为电力网与热力网的平衡节点,其余EH节点均视为PV/PQ节点、普通热源节点及天然气负荷节点。此时,各能源网络处于强耦合状态,电网侧及热网侧的潮流均由EH单元实现平衡与调节。

    • 在开展综合能源系统稳态能流计算前,需要明确各能源网络所包含的状态量。电力/热力/天然气供能网络通过EH进行耦合,各网络间根据耦合程度的不同将获得不同的运行状态量。

      对于电力网络,状态量主要包括节点电压幅值与相角、节点注入的有功与无功功率,节点类型主要有平衡节点、PV与PQ节点3类;对于热力网络,状态量主要包括管道热媒流量、供回水温度以及节点注入的热功率,节点类型主要有热源节点(又分为平衡热源节点与普通热源节点2类)、负荷节点(已知负荷热功率及节点回水温度)以及其他节点;对于天然气网络,状态量主要包括节点压力及管道气流量,节点类型主要包括压力已知节点和流量已知节点,其中压力已知节点多为气源节点且至少有一个气源节点设置为天然气平衡节点,流量已知节点则多为负荷节点。各能源网络具体的节点类型及待求状态量等信息如表1所示。

      表 1  电–热–气综合能源系统状态量信息

      Table 1.  Information of state variables for IES containing electricity, thermal and gas networks

      能量网络类型节点类型已知状态量未知状态量
      电力网络电力平衡节点节点电压幅值、相角节点注入有功无功功率
      PV节点节点注入有功功率、节点电压幅值节点注入无功功率、节点电压相角
      PQ节点节点注入有功无功功率节点电压幅值、相角
      热力网络热力平衡节点供水温度回水温度、供热功率、管道热媒流量
      热负荷节点回水温度、热功率供水温度、管道热媒流量
      其他热源节点供水温度、供热功率回水温度、管道热媒流量
      天然气网络天然气平衡节点节点压力节点注入气流量、天然气管道流量
      流量已知节点节点注入气流量节点压力、天然气管道流量
    • 综合能源系统多能流求解方法主要有统一求解与分布式顺序求解2类。统一求解法虽然能提供统一的求解模式,便于理解,但在求解时包含大量的待求解变量,雅克比矩阵维度巨大,求解复杂度较高,且模型收敛性也无法保证。而分布式顺序求解法能够充分利用电、气、热各网络潮流已有的计算求解模型,实现各能源网络的独立求解,大大提高了求解效率,缩减了求解复杂度。因此,本文主要依托分布式顺序求解法对前节所述的EH 4种典型运行模式开展多能流求解方法的研究。

      图3所示为本文综合能源系统多能流求解流程。模式1、模式2与模式3求解流程类似,各能源子系统利用各自现有成熟的求解模型依次进行独立求解;模式4由于各能源系统耦合程度较深,求解流程相对更为复杂。选择模式1与模式4两种具有代表性的运行场景下的能流计算过程进行分析说明。

      图  3  综合能源系统多能流求解流程

      Figure 3.  Flowchart for multi-energy flow calculation of IES

      1)对于模式1,其求解流程相对简单,归纳如下:

      步骤1:获取各EH输出端所接的气、电、热、冷负荷信息;LgLeLhLc,根据式(3)描述的EH单元内部耦合关系矩阵分别计算各EH输入功率$P_{\rm{e}}^{{\rm{net}}} $$ P_{\rm{g}}^{{\rm{net}}} $$ P_{\rm{h}}^{{\rm{net}}} $,各EH单元可视为电力、热力及天然气网络对应接入节点的等效负荷;

      步骤2:热力系统潮流的独立求解;基于热网侧潮流计算,热网中各循环水泵的耗电功率等效为电力网络侧的电负荷;

      步骤3:天然气系统潮流的独立求解。若天然气系统中含有电力驱动型压缩机,则还需额外计算压缩机的耗电功率,等效至电力网络侧的电负荷;

      步骤4:基于电力网络侧已知参数,利用牛–拉法求解电力网络潮流;

      步骤5:输出电/气/热网络多能流计算结果。

      2)对于模式4,由于同时选择两个EH单元分别作为电力与热力系统的平衡节点,因此各能源子系统间的耦合程度最高。与其他几种模式相比,模式4在各能源子系统能流独立求解的基础上,需要额外增加考虑耦合元件输出功率偏差约束的联合迭代求解过程。可归纳如下。

      步骤1:首先选择两个EH单元分别作为电力侧与热力侧的平衡节点,同时确定除平衡节点外的其他EH单元运行模式;

      步骤2:依据EH耦合关系及运行模式计算除平衡节点外的EH单元热电输出功率Po,hPo,e及输入侧天然气功率$P_{\rm{g}}^{{\rm{net}}} $

      步骤3:进入热电系统潮流联合迭代求解过程。

      ①首先初始化电力平衡节点处EH单元的输出电功率初始值$ P_{{\rm{o}},{\rm{e}}}^0 $

      ②其次,根据EH耦合关系及运行模式计算电力平衡节点的输出热功率$P_{{\rm{o}},{\rm{h}}}^{{\rm{hslack}}}$

      ③基于热力侧已知条件,进行热力系统的潮流求解,若热力侧潮流不收敛,则需修正EH单元的运行模式及相关运行参数,返回步骤2,直至热力潮流收敛;

      ④基于当前迭代过程中的热力潮流计算结果,获取热力平衡节点处的EH单元输出热电功率并计算水泵耗电功率;

      ⑤基于电力侧已知条件,进行电力子系统的潮流计算,若电力侧潮流不收敛,则同样需修正EH单元的运行模式及相关运行参数,返回步骤2,直至电力潮流收敛;

      ⑥获取电力平衡节点处EH单元的热电输出功率,计算相邻两次迭代过程的电力平衡节点处EH单元电输出功率偏差,若满足收敛条件$\left| {\Delta P_{{\rm{o,e}}}^t} \right| \leqslant {\varepsilon _2}$,则热电系统潮流联合迭代求解过程完成,输出热力/电力侧的潮流计算结果,反之,则需用当前第t次迭代时电力平衡节点处EH单元的电输出功率$P_{{\rm{o,e}}}^t$初始化第(t+1)次迭代过程的$ P_{{\rm{o}},{\rm{e}}}^0 $,重新依次执行步骤3,直至迭代过程完成;

      步骤4:进行天然气系统潮流的独立求解。同样需要处理不同压缩机类型的等效电力负荷计算问题,处理方式与模式1相同;

      步骤5:输出电/气/热网络多能流计算结果。

      此外,若考虑P2G技术的应用,则考虑采用合同电转气的商业模式[21],将P2G环节对应的天然气与电力节点分别等效为定流量节点与PV节点,进而参与多能流的计算。

    • 为验证本文所提多能流计算方法的有效性,以改进的IEEE14节点电力系统、15节点区域热力系统、14节点天然气系统以及两个EH(EH1和EH2)构成的电–热–气综合能源系统作为算例进行仿真分析。

      图4所示,电、气、热各能源子系统通过两个EH单元进行耦合。天然气子系统的平衡节点为Gs1,通过外部气源输入调节整个天然气网络的平衡,气源节点的压力已知为1000 psia,压缩机为燃气轮机驱动型。EH1的热网侧接入点为Hs6,电网侧接入点为Bs3,天然气网接入点Gs14视为天然气负荷节点;EH2的电网侧接入点为Bs1,热网侧接入点为Hs7,天然气接入点Gs3视为天然气负荷节点。其中,热力子系统中各热力负荷节点回水温度设置为30 ℃,热源节点供水温度为90 ℃。两个EH的拓扑结构如图2所示,但暂不考虑输冷环节,因此忽略吸收式制冷机及蓄冷等设备。EH单元的配置参数以及各能源子系统的基本参数见附表A1A6

      图  4  电-热-气综合能源系统算例结构

      Figure 4.  Case diagram of IES containing electricity, thermal and gas networks

    • 基于本文所提方法,利用MATLAB编写程序,针对2种不同场景对算例进行多能流求解与分析,计算结果如表26所示。

      表 2  电力子系统潮流计算结果

      Table 2.  Calculation results for electricity flow

      支路编号支路有功功率/MW节点编号电压/pu
      S1S2S1S2
      d14.6153.479Bs11.0601.035
      d21.8071.128Bs21.0451.031
      d30.7630.508Bs31.0101.010
      d42.9112.478Bs41.0381.029
      d51.6561.207Bs51.0381.027
      d6−0.421−0.570Bs61.0701.070
      d7−5.319−5.360Bs71.0791.072
      d81.9021.560Bs81.0781.072
      d90.9450.746Bs91.0881.083
      d10−2.679−3.848Bs101.0841.080
      d111.3921.725Bs111.0771.075
      d120.7390.780Bs121.0701.070
      d132.0702.242Bs131.0711.070
      d140.8000.800Bs141.0791.076
      d151.1020.760
      d16−0.101−0.432
      d170.6740.463
      d18−1.009−1.340
      d190.1280.169
      d200.8431.056

      表 3  电力子系统各发电单元出力

      Table 3.  Output of each electricity generation unit

      发电单元有功出力/MW
      S1S2
      EH15.98195.233
      EH26.48178
      G133
      G286.706

      表 4  热力子潮流计算结果

      Table 4.  Calculation results for thermal flow

      管道编号管道流量/(kg/s)节点编号供水网络节点输出温度/℃
      S1S2S1S2
      l19.1739.174Hs189.97689.987
      l29.1749.174Hs289.90089.978
      l39.1939.178Hs389.88789.956
      l49.1929.177Hs489.88489.985
      l59.1949.184Hs589.99689.994
      l69.1939.182lh189.95489.952
      l77.5947.584lh289.94989.948
      l87.5947.584lh389.82489.925
      l97.5847.583lh489.82989.931
      l107.5847.583lh589.81689.885
      l114.86122.427lh689.82889.897
      l12−10.3277.259lh789.82789.905
      l13−28.714−11.107lh889.82889.906
      l14−47.100−29.462lh989.90389.914
      l15−65.447−47.810lh1089.90389.914
      l16−20.030−37.593

      表 5  天然气子系统潮流计算结果

      Table 5.  Calculation results of natural gas flow of natural gas subsystem

      管道编号管道气流量/MMCFD节点编号节点压力/psia
      S1S2S1S2
      gl1201.448200.375Gs11000.0001000.000
      gl2196.281193.406Gs2643.607608.896
      gl349.42136.356Gs3627.115596.632
      gl4118.502110.748Gs4498.023476.758
      gl5145.009−190.645Gs5647.429619.786
      gl677.456−86.392Gs6384.860909.544
      gl7−109.080155.356Gs7538.8041273.362
      gl8117.976117.855Gs8485.013719.426
      gl9123.748123.653Gs91000.0001000.000
      gl10−33.631−33.627Gs10883.147864.692
      gl1141.60741.482Gs111100.0001100.000
      gl1250.11750.026Gs12691.885692.655
      Gs13723.185723.913
      Gs14645.180646.296

      表 6  天然气子系统压缩机耗气量计算结果

      Table 6.  Gas consumption calculation results of the compressors of natural gas subsystem

      压缩机编号压缩机耗气量/MMCFD压缩机编号压缩机耗气量/MMCFD
      S1S2S1S2
      10.3410.31930.6570.315
      20.5330.70140.4660.405

      场景1:运行于模式4。EH1的热网侧接入点Hs6视为热力平衡节点,电网侧接入点Bs3视为PV节点,天然气网接入点Gs14视为天然气负荷节点;EH2的电网侧接入点Bs1视为电力平衡节点,热网侧接入点Hs7视为普通热源节点,天然气接入点Gs3视为天然气负荷节点。

      场景2:运行于模式2。电力侧并网运行,并网点(电力平衡节点)为Bs6,外部电网等效为G2,即电力侧网络的平衡由G2调节。EH1的热网侧接入点Hs6视为热力平衡节点,电网侧接入点Bs3视为PV节点,天然气网接入点Gs14视为天然气负荷节点。而EH2的电网侧接入点Bs1视为PV节点,热网侧接入点Hs7视为普通热源节点,天然气接入点Gs3视为天然气负荷节点。

      根据混合潮流计算结果可知,对于电力子系统,2种不同场景下支路的有功潮流分布及节点电压分布产生一定差异,但整体数值较为接近。场景2中EH2的有功出力设置为8 MW,且各PV节点的无功储备充足,经计算平衡节点对应的G2机组有功出力为6.76 MW,2种场景下的节点电压幅值均在0.9~1.1 pu之间,未出现电压越限情况。此外2种场景仅运行模式不同,网络参数基本保持一致,因此2种场景下的潮流未发生明显转移,潮流计算结果较为接近。

      对于热力子系统,由计算结果可知,不同运行场景下的负荷节点管道流量变化较小,但主管道l11—l16对应的管道热水流量差异明显。两种场景下,EH1的热网接入点均为热力平衡节点,承担热网侧的潮流平衡调节。场景1下EH1与EH2的热功率输出分别为16.43 MW与5.03 MW,管道热水汇流节点为Hs2;场景2下EH1与EH2的热功率输出分别为12 MW与9.43 MW,管道热水汇流节点变为Hs3。可见,EH1热源输出热功率的降低使得EH1附近的管道l14、l15流量降低,而EH2热源附近管道l11与l16的流量随着EH2单元热功率出力的增加而增大。

      对于天然气子系统,部分节点等效气负荷的改变使得天然气管道流量与节点压力分布产生变化。场景2下管道gl5、gl6、gl7的气流量相较于场景1明显增加,节点Gs6、Gs7、Gs8压力值也存在明显提升,可能对系统的安全稳定运行产生影响。因此有必要考虑IES系统在不同运行模式下的能流分布情况,识别脆弱环节,通过合理配置EH单元及网络参数,保障系统在不同运行模式下均具有足够的安全裕度。

      此外,图5展示了本文所提多能流求解算法在不同场景下的收敛曲线。由图5可知,此算法有着较为理想的收敛速度,仅需经过少量迭代即可收敛到很小的误差范围内,从而验证了本节算法的有效性。同时,由于场景2下热/电/气网络的耦合程度相对较低,求解复杂度较低,因此与其他场景相比具有更快的收敛速度。

      图  5  不同场景下的收敛曲线

      Figure 5.  Convergence curves in different cases

      此外,基于本文所提能量流求解方法,进一步分析EH单元在不同运行模式下IES的N-1静态安全校验情况。同样以场景1与场景2为例进行分析,对应的关键元件N-1校验结果如附图B1附图B2所示。由图可见,对于场景1,EH2单元处于离网运行状态,各能源子系统耦合程度较高,在关键元件(发电机组G1、G2)N-1情况下,电力侧的潮流越限将通过耦合单元影响其他能源系统的运行状态,可能引起有故障的跨能源系统传播,增大了系统整体的安全运行风险。对于场景S2,EH2单元运行于并网状态,系统整体的耦合程度较低,电力侧安全运行的调节能力更强,使得系统具有比场景S1更为理想的静态安全性能。

      图  B1  场景S1下关键元件N-1校验

      Figure B1.  N-1 static security analysis of key components in scene S1

      图  B2  场景S2下关键元件N-1校验

      Figure B2.  N-1 static security analysis of key components in scene S2

      可见,EH单元的不同运行模式将对IES系统整体运行状态产生影响,强耦合下的IES系统在具有高灵活性与互补互济特性的同时往往更容易引起较低的安全运行裕度。同时也验证了本文所提能流计算方法在系统运行状态与安全性分析方面的有效性。

    • 1)EH作为各能源子网络的耦合互动桥梁,不同的运行模式将对不同能源子系统间的耦合深度以及整个综合能源能量流分布产生显著影响,体现了考虑能量枢纽运行模式的必要性。

      2)本文所提混合能量流求解算法具有较好的收敛性,避免了牛–拉法统一求解系统能量流时可能存在的模型维度大、收敛性无法保证等问题,具有一定的应用价值。

      此外,对于多时间尺度下能源网络的动态特性研究还有待深入,同时可考虑结合系统实际运行数据,对模型作进一步修正,提高其准确性与适应性。

      (本刊附录请见网络版,印刷版略)

参考文献 (21)

目录

    /

    返回文章
    返回