CN110348602B - 计及天然气管网和热力管网特性的综合能源系统优化方法 - Google Patents
计及天然气管网和热力管网特性的综合能源系统优化方法 Download PDFInfo
- Publication number
- CN110348602B CN110348602B CN201910491982.1A CN201910491982A CN110348602B CN 110348602 B CN110348602 B CN 110348602B CN 201910491982 A CN201910491982 A CN 201910491982A CN 110348602 B CN110348602 B CN 110348602B
- Authority
- CN
- China
- Prior art keywords
- energy
- power
- network
- pipeline
- natural gas
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 title claims abstract description 190
- 239000003345 natural gas Substances 0.000 title claims abstract description 97
- 238000005457 optimization Methods 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000004146 energy storage Methods 0.000 claims abstract description 36
- 230000008878 coupling Effects 0.000 claims abstract description 28
- 238000010168 coupling process Methods 0.000 claims abstract description 28
- 238000005859 coupling reaction Methods 0.000 claims abstract description 28
- 238000009826 distribution Methods 0.000 claims abstract description 27
- 238000010438 heat treatment Methods 0.000 claims abstract description 18
- 238000010276 construction Methods 0.000 claims abstract description 14
- 239000007789 gas Substances 0.000 claims description 75
- 230000005611 electricity Effects 0.000 claims description 18
- 238000003860 storage Methods 0.000 claims description 18
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 17
- 230000005540 biological transmission Effects 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 10
- 239000000126 substance Substances 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 7
- 238000007599 discharging Methods 0.000 claims description 6
- 150000001875 compounds Chemical class 0.000 claims description 4
- 230000006835 compression Effects 0.000 claims description 4
- 238000007906 compression Methods 0.000 claims description 4
- 238000005338 heat storage Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 3
- 238000009434 installation Methods 0.000 claims description 3
- 230000014759 maintenance of location Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 238000005303 weighing Methods 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 2
- 239000007924 injection Substances 0.000 claims description 2
- 230000002708 enhancing effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 7
- 238000005485 electric heating Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000006467 substitution reaction Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009194 climbing Effects 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000001737 promoting effect Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000003139 buffering effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000010206 sensitivity analysis Methods 0.000 description 1
- 230000002195 synergetic effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种计及天然气管网和热力管网特性的综合能源系统优化方法。本发明的综合能源系统优化方法,其包括步骤:1)构建含多能流耦合设备和储能设备的能源中心设备模型;2)构建含电力网络、天然气管网和热力管网的能源网络模型;3)以综合能源系统运行周期内总成本最小为优化目标,考虑综合能源系统建设约束和运行约束,建立计及天然气管网和热力管网特性的综合能源系统优化模型。本发明所提方法能够提供综合能源系统协同规划的技术方案,在增强系统灵活性的同时,提高综合能源利用率。
Description
技术领域
本发明涉及电力系统优化领域,具体地说是一种计及天然气管网和热力管网特性的综合能源系统优化方法。
背景技术
能源可持续发展形势的日益严峻,促使各国打破各能源系统单独规划、独立运行的既有模式,开展多能流综合利用的研究。多种能源系统在规划、设计、建设和运行阶段协调配合,能够推动多能流互补共济,促进可再生能源消纳,提升能源整体利用效率,增强能源系统灵活性。能源中心将综合能源系统中的多能流耦合设备和储能设备抽象为一个输入-输出双端口网络模型,模型中多种能流分别从两个端口输入与输出,简化了综合能源系统中复杂的多能流耦合关系。在此基础上,综合能源系统规划问题可以分为能源中心规划与能源网络规划两部分内容。目前,关于能源中心的规划问题已经开展了较为充分的研究。
能源中心的规划大多建立在优化运行的基础上,主要关注能源中心的多能流耦合设备和储能设备的选址定容,但忽略了能源网络特性的影响。然而,综合能源系统中的能源中心往往不是独立运行的,对于综合能源系统规划问题,除考虑能源网络规划之外,还需要考虑能源网络特性对能源中心运行的影响。目前,关于能源网络特性对能源中心运行影响的研究主要集中在气网管存、热网损耗和热网延时方面。针对含多个能源中心的综合能源系统规划问题,部分研究考虑了天然气管网或热力管网,还未有兼顾天然气与热力管网特性的综合能源系统优化方法研究。
发明内容
本发明所要解决的技术问题是克服上述现有技术存在的缺陷,提供一种计及天然气管网和热力管网特性的综合能源系统优化方法,其在满足多区域电、气、热负荷需求以及综合能源系统优化运行的基础上,得到多能量枢纽与能源网络的协同规划方案,增强能源系统的灵活性,提高综合能源利用效率。
为此,本发明采用如下的技术方案:计及天然气管网和热力管网特性的综合能源系统优化方法,其包括步骤:
1)构建含多能流耦合设备和储能设备的能源中心设备模型;
2)构建含电力网络、天然气管网和热力管网的能源网络模型;
3)以综合能源系统运行周期内总成本最小为优化目标,考虑综合能源系统建设约束和运行约束,建立计及天然气管网和热力管网特性的综合能源系统优化模型。
本发明构造了计及管网特性的综合能源系统优化数学模型,通过MATLAB环境下的YALMIP/GUROBI求解便可得到综合能源系统中能量枢纽设备与能源网络的最优规划技术方案。
进一步地,步骤1)中,所述能源中心设备模型抽象为一个输入-输出双端口网络模型,多种能流分别从两个端口输入与输出,多能流耦合设备和储能设备的输入和输出端按照能源形式,分别汇集到同一端点。
进一步地,步骤1)中,所述的多能流耦合设备包括电热锅炉、燃气锅炉、燃气轮机和热电联产机组,其能源转移效率统一表示为:
所述储能设备包括储电、储气和储热设备,储能设备的运行约束统一表示为:
所述能源中心设备模型两个端口的输入与输出功率需满足:
进一步地,步骤2)中,采用直流潮流模型描述电力网络:
电力网络节点能量平衡约束表示为:
进一步地,步骤2)中,所述的天然气管网中,天然气管道约束有:
根据气体状态方程与波义耳定律,管存相关的计算式如下:
且其满足质量守恒定律,如下式所示:
其中,
式中:Vij,t为天然气管道ij内的管存量;pi,t和pj,t分别为管道ij首末端的气压;和分别为管道ij的出、入口流量;和分别为管道ij内径和长度;Rgas为通用气体常数;为管道ij的管存系数;Mgas为天然气分子量;Tg、ψ和ρg分别为天然气温度、压缩因子和相对空气的密度;Δt为单位时间段的时长;
此外,天然气管道传输的气流量与首末端气压有关,实际运行中大多数输气管道以高雷诺数的流速运行即处于湍流状态,满足管道气流方程,如下式所示,参数均折算到标准状况下:
pi,min≤pi,t≤pi,max,
其中
所述的天然气管网中,加压站约束表示为:
pi,t≤ξcompj,t,
天然气网络节点能量平衡约束表示为:
式中:为天然气管网中与节点k相连的节点集合;和分别为管道jk出口和入口端的气功率;为外部气源注入综合能源系统的气功率;为注入能源中心的气功率;为天然气热值;和分别为管道ik的出、入口流量;ξcom表示加压站最大加压系数。
进一步地,步骤2)中,所述的热力管网中,换热站约束有:
供水管和回水管的出入口温度约束表示如下:
热负荷和能源中心与换热站热交换约束表示如下:
热力管网节点热量守恒约束表示如下:
式中:和和分别为第k个能源中心/第f个热负荷的供水管与回水管的出入口温度;和分别为第k个能源中心和第f个热负荷与其换热站的热交换功率;cw为水的比热容;和分别为单位时间内流过换热站的工质质量;NZ为流入汇集点z的管道集合;Tz,t和分别为汇集点z和管道b出口的工质温度;为单位时间内管道b流出的工质质量;
所述的热力管网中,热网延时效应约束有:
其中,
式中:和分别为热力传输延时时长的上、下限;和分别为未计及温度损失时的管道出、入口温度;ρw为热力管网工质的密度;和分别为t-γb,t和t-φb,t+1时刻到t时刻注入管道的工质质量;Ν为正整数集合,n表示其中的元素;和分别为δt时间内流入和流出管道b的工质质量;和分别表示t-φb,t和t-γb,t时刻注入管道的工质温度;Ab和分别表示管道的横截面面积和长度。
所述的热力管网中,热网损耗约束有:
由于工质在传输过程中不可避免地与管道进行热交换而产生热量损耗,故管道出口温度根据苏霍夫温降公式进行修正:
其中,
所述的热力管网中,热网节点能量平衡约束有:
能源中心和热负荷的热能平衡约束:
进一步地,步骤3)中,综合能源系统优化模型中,目标函数表示为:
其中,
式中:用下标s表示第s个场景;Cinv、和Ctotal分别表示计及设备残值的投资成本、第τ年的外部能源购买成本和系统运行周期内总成本;r为贴现率;Hor为规划年限;D为一年的天数;NS为一年内的场景集合;Νeh和Nbr分别为综合能源系统拓扑结构中的节点集合与支路集合;ΝX和Nnet分别为综合能源系统中的能源中心设备种类集合和能源网络种类集合;和为第k个能源中心中候选X类设备的集合和能源网络κ中节点i与j之间的候选线路或管道集合;ωs为场景s发生的概率;Φ为一个典型日的时间分段数;和分别为从外部购电和购气功率;和分别为电能和天然气的单位购买成本;假定投运都发生在年初,Rx、cx、βx和Sx分别为x的规划期末残值率、单位容量投资成本、候选设备投运状态和单台/条/回容量;Δt为单位时间段的时长;
假设能源中心设备与能源网络折旧程度与投运时间呈线性关系,x的残值率统一描述为:
进一步地,步骤3)中,综合能源系统优化模型中,建设约束为:
综合能源系统的投资成本包括多能流耦合设备、储能设备以及电力网络、天然气管网和热力管网的建设成本,投资成本存在上限,如下式所示:
对于能源中心设备和能源网络,设备安装台数和线路或管道的建设条/回数需要满足下面的约束:
进一步地,步骤3)中,综合能源系统优化模型中,运行约束为:
能源中心中的设备输入功率和爬/滑坡速度约束条件统一表示为:
在能源网络中,两节点之间建设多条并行线路,由于能源网络的非线性,需要分别计算每一条线路的运行状态,能源网络线路功率约束统一表示为:
式中:和分别为节点i与j之间第l条电力网线路的传输功率和天然气管道的出入口功率;为向热负荷f供热的第l回管道的传输功率;ζe,tran、ζg,tran和ζhex为候选电力线路、天然气管道和热力管道的容量裕度;0-1变量和为候选电力线路、天然气管道和热力管道的投运状态;和为候选电力线路、天然气管道和热力管道的容量;
从外部注入的电功率和气功率需要满足下面的约束:
进一步地,利用增量法对非线性约束进行线性化处理;
对于非线性函数h(y),线性化方法简述如下:权衡计算精度与计算量,将自变量的取值范围分为υ个区间;计算区间各个分段点Yi处的函数值;函数即表示为下式:
本发明构建了含电力网络、天然气管网和热力管网的综合能源系统模型,提出了计及管网特性影响的综合能源系统优化规划方法。所提出的规划方法在满足多区域电、气、热负荷需求以及综合能源系统优化运行的基础上,得到多能量枢纽与能源网络的协同规划方案。通过算例结果的分析,验证了在综合能源系统规划中考虑能源网络的必要性与可行性。
本发明所提方案表现出储气设备协同CHP机组与燃气轮机进行配置,以负荷定天然气管网配置,以及储热设备协同电热锅炉进行配置等特点,增强了能源系统的灵活性,提高了综合能源利用效率;天然气管网特性影响主要体现在储气设备规划方面;热力管网特性的影响主要体现在能量耦合设备(电热锅炉、CHP机组等)以及热力管网的选址定容方面。本发明所提方案还可以在引导电能替代、推动“电-气-热”多能流互补共济以及促进“源-网-荷-储”协同发展等方面发挥作用。
附图说明
图1为本发明实施例中综合能源系统能源中心架构图;
图2为本发明实施例中环式热力管网典型结构图;
图3为本发明实施例中热网延时效应示意图;
图4为本发明应用例中6节点综合能源系统结构框架图;
图5为本发明应用例中网络特性对综合能源系统优化规划的影响图;
图6为本发明应用例中不同负荷规模的综合能源系统优化规划图;
图7为本发明应用例中不同负荷热电比的综合能源系统优化规划图;
图8为本发明实施例中综合能源系统规划方法的流程图。
具体实施方式
下面结合附图并通过具体实施方式来进一步说明本发明的技术方案。本领域的技术人员应该明了,所述的实施例仅是帮助理解本发明,不应视为对本发明的具体限制。
实施例
本实施例为一种计及天然气管网和热力管网特性的综合能源系统规划方法,其包括以下步骤:
步骤1、构建能源中心设备模型
构建能源中心的架构如图1所示,包括多能流耦合设备和储能设备。能源中心抽象为一个输入-输出双端口网络模型,模型中多种能流分别从两个端口输入与输出,多能流耦合设备和储能设备的输入和输出端可以看作分别按照能流种类汇集到同一点。
1)多能流耦合设备
能源中心内部的多能流耦合设备发挥能源转换器的作用,通过其内部的电、气、热多能流互补共济,可以满足多种负荷的用能需求。多能流耦合设备包括电热锅炉、燃气锅炉、燃气轮机和热电联产机组等。多能流耦合设备可以统一表示为:
2)储能设备
储能设备是能源中心中的重要设备,启停时间短,功率爬坡速度快,能够在短时间内响应供用能侧的功率变化。储能设备包括储电、储气和储热设备等。储能设备的运行约束可以统一表示为:
3)能源中心端口
能源中心两个端口的输入与输出功率需满足:
步骤2、能源网络模型构建
1)电力网络
采用直流潮流模型描述电力网络:
电力网络节点能量平衡约束表示为:
2)天然气管网
综合能源系统中的天然气系统通常由气源、管道、压缩机和气负荷等组成。
A.天然气管道约束
天然气的传输速度远不及电力且具有压缩性,因此管道输入与输出流量不必时刻相等,管存表现出一定的缓冲作用。根据气体状态方程与波义耳定律,管存相关的计算式如式(8)所示,且其满足质量守恒定律,如式(9)所示。
其中
式中:Vij,t为天然气管道ij内的管存量;pi,t和pj,t分别为管道ij首末端的气压;和分别为管道ij的出入口流量;和分别为管道ij内径和长度;Rgas为通用气体常数;为管道ij的管存系数;Mgas为天然气分子量;Tg、ψ和ρg分别为天然气温度、压缩因子和相对空气的密度。
此外,天然气管道传输的气流量与首末端气压有关。实际运行中大多数输气管道以高雷诺数的流速运行即处于湍流状态,满足管道气流方程,如式(11)-(12)所示;式(13)则表示气网节点气压的上下限约束。本实施例中参数均折算到标准状况下。
pi,min≤pi,t≤pi,max (13)
其中
B.加压站约束
由于天然气管网内部存在摩擦力,气压会逐步衰减,因此天然气管网中一般安装有加压站,用于提升天然气管道内的气压。加压站模型可以简单表示为:
pi,t≤ξcompj,t (15)
C.气网节点能量平衡约束
天然气网络节点能量平衡约束表示为:
3)热力管网
热力系统通常由热源、环式管网、换热站和热负荷等组成。环式热力管网典型结构如图2所示。
A.换热站约束
综合能源系统中的热能由热力管网的工质进行传递,并在换热站进行热量交换,传递与交换的热功率大小与各节点温度有关。下列约束中,式(19)-(22)表示供水管和回水管的出入口温度约束;式(23)和式(24)分别表示热负荷和能源中心与换热站热交换约束;式(25)表示热力管网节点热量守恒约束。
式中:和和为第k个能源中心/第f个热负荷的供水管与回水管的出入口温度;和分别为第k个能源中心和第f个热负荷与其换热站的热交换功率;cw为水的比热容;和分别为单位时间内流过换热站的工质质量;NZ为流入汇集点z的管道集合;Tz,t和分别为汇集点z和管道b出口的工质温度;为单位时间内管道b流出的工质质量。
B.热网延时效应约束
热力管网的工质在管网中流动需要足够时间且存在一定损耗。热力传播速度近似等于载体流动速度,故热力管网延时特性可以用加权平均的方法进行描述。图3为热力管网纵截面,右侧阴影部分为t时段流出管道的工质,和分别为δt时间内流入和流出管道b的工质质量。如式(26)所示,流出的工质温度可以由三部分温度的加权平均数表示。
其中
C.热网损耗约束
由于工质在传输过程中不可避免地与管道进行热交换而产生热量损耗,故管道出口温度可以根据苏霍夫温降公式进行修正:
其中
D.热网节点能量平衡约束
能源中心与热负荷均满足热能平衡。式(33)-(34)分别为能源中心和热负荷的热能平衡约束。
步骤3、优化规划模型
1)目标函数
采用场景分析的方法,将一年的负荷情况削减成s个场景。以综合能源系统规划期内能源投资运行总成本最小为优化目标,决策变量为候选的多能流耦合设备、储能设备、电力线路、天然气管道和热力管道的投运状态;能源中心中可以投运多台多种型号的设备,能源网络两节点之间可以投运多条并行线路或管道。此外,规划过程中不改变综合能源系统的拓扑结构。本实施例中用天然气管道功率上限描述天然气系统容量,用热力管道功率上限描述热力系统容量。目标函数可以表示为:
其中,
式中:s表示场景编号;Cinv、和Ctotal分别表示计及设备残值的投资成本、第τ年的外部能源购买成本和系统投资运行总成本;r为贴现率;Hor为规划年限;D为一年的天数;NS为一年内的场景集合;Νeh和Nbr分别为综合能源系统拓扑结构中的节点集合与支路集合;ΝX和Nnet分别为综合能源系统中的能源中心设备种类集合和能源网络种类集合;和为第k个能源中心中候选X类设备的集合和能源网络κ中节点i与j之间的候选线路或管道集合;ωs为场景s发生的概率;Φ为一个典型日的时间分段数,本文设定Φ为24小时;和分别为从外部购电和购气功率;和分别为电能和天然气的单位购买成本;假定投运都发生在年初,Rx、cx、βx和Sx分别为x的规划期末残值率、单位容量投资成本、候选设备(线路或管道)投运状态和单台(条或回)容量。
假设能源中心设备与能源网络折旧程度与投运时间呈线性关系,x的残值率可以统一描述为:
2)约束条件
A.建设约束
综合能源系统的投资成本包括多能流耦合设备、储能设备以及电力网络、天然气管网和热力管网的建设成本,投资成本通常存在上限,如式(39)所示:
式中为综合能源系统投资成本上限。
对于能源中心设备和能源网络,设备安装台数和线路或管道的建设条(回)数需要满足式(40)和(41)的约束:
B.运行约束
能源中心中的设备输入功率和爬(滑)坡速度约束条件统一表示为:
在能源网络中,两节点之间可以建设多条并行线路。由于能源网络的非线性,需要分别计算每一条线路的运行状态。节点与之间第条(回)线路的功率约束统一表示为:
式中:和分别为节点i与j之间第l条电力网线路的传输功率和天然气管道的出入口功率;为向热负荷f供热的第l回管道的传输功率;ζe,tran、ζg,tran和ζhex为候选电力线路、天然气管道和热力管道的容量裕度;0-1变量和为候选电力线路、天然气管道和热力管道的投运状态;和为候选电力线路、天然气管道和热力管道的容量。
从外部注入的电功率和气功率需要满足式(47)和(48)的约束:
利用增量法对非线性约束进行线性化处理。对于非线性函数h(y),线性化方法简述如下:权衡计算精度与计算量,将自变量的取值范围分为υ个区间;计算区间各个分段点Yi处的函数值;函数即可以表示为式(49)。其中,μi为连续变量,代表每个分段上的占比;为0-1变量,用于确保增量法可以表示可行域内的所有函数值。对于天然气管网非线性约束,依次线性化式(11)中的三个平方项,然后进行线性叠加,即完成线性化。
对于建立的混合整数线性优化模型,采用YALMIP/GUROBI求解器进行求解,便可得到综合能源系统中能量枢纽与能源网络的协同规划结果。
应用例
参数设置:以包含一个6节点综合能源系统为例进行说明,综合能源系统如图4所示。能源中心1、能源中心2和能源中心3均承载电、气、热三类负荷;其他能源中心只承载电、气两类负荷。外部电网通过节点1、2和6向综合能源系统供电,外部气源通过节点3和6向综合能源系统供气,能源中心1、能源中心2和能源中心3均通过环式网络向各自区域内两个热负荷供热;能源网络线路编号见表1。将电、气、热的日负荷曲线分为夏季、过渡季和冬季三个典型场景;在规划期内,从外部购买天然气和电能的单位成本按照折现率增长,初年电价选用浙江省峰谷电价,初年天然气价格设为3.25元/m3。其他参数见表1至表5。
表1综合能源系统候选网络参数
表2综合能源系统投资参数
表3综合能源系统候选设备参数
表4天然气管网参数
表5热力管网参数
采用YALMIP/GUROBI求解器进行求解,计及管网特性的综合能源系统优化规划方案如表6和表7所示。
表6多能流耦合设备与储能设备优化规划方案
表7综合能源系统能源网络优化规划结果
计及天然气管网和热力管网特性的综合能源系统规划与网络特性密切相关。以表1中的规划方案为参照情景,对下述三个场景进行对比分析:
情景1:忽略天然气管存效应影响,假设天然气管道出入口流量时刻保持一致。
情景2:忽略热力管网延时效应,假设热力管网出入口温度变化趋势时刻保持一致。
情景3:忽略热力管网热量损耗,假设热力管网不发生热量损耗。
在情景1中,优化规划方案中储气设备总容量由88MW增加到为104MW,储气设备容量如图5所示;这是由于在参照情景中,天然气管道出入口流量不必时刻相等,管网表现出储能设备的特性,即可以在一定范围内通过调节气压控制天然气流量,对储能设备起到替代作用。
在情景2中,优化规划方案中电热锅炉数目减少,供热总容量如图5所示。结合热力管网1的供水管运行数据可以发现,在参照场景中,t时段流出管道的工质温度等于t-1和t-2两个时段流入工质温度的加权平均数,该特性会增大供水管中流入工质的温度调节幅度,即在忽略延时效应的情形下,供热负荷(热负荷与管网热量损耗之和)的波动性小于参照场景,故对电热锅炉的容量需求小于参照场景。
在情景3中,优化规划方案用一台热功率较低的CHP机组替换了原有的两台电热锅炉,供热总容量如图5所示。这是由于忽略热量损耗后,能源中心供热负荷的计算值小于实际值,而该规划方案在实际运行中无法完全达到供需平衡。
计及天然气管网和热力管网特性的综合能源系统规划与负荷规模和负荷热电比密切相关。对负荷规模和负荷热电比两个因素进行敏感性分析:
将负荷规模从减少30%逐次调整到增加30%,分析其对综合能源系统优化规划的影响。如图6所示,随着电、气、热负荷需求的增长,优化规划方案主要变化趋势如下:其一是外部电能对燃气轮机的替代;在规划方案中能源中心4、能源中心5之间天然气管道容量需求减小,能源中心4处不再配置燃气轮机,同时也伴随储气设备减少。其二是电热锅炉对CHP机组的替代;能源中心中CHP机组总数减少,并伴随电热锅炉增多。这是因为天然气气源每日提供的天然气容量有限,需要优先满足各能源中心的气负荷;在天然气较为充裕时,CHP机组可表现出高效利用能源的优势。
原始数据的夏季、过渡季和冬季负荷热电比依次为0.149、0.260和0.962,逐次调整热负荷占比并保持总负荷不变。如图7所示,随着热负荷占比的增加,先后出现电热锅炉增加、输电网容量减少、热力管网扩容和CHP机组替代燃气轮机。也即在本文算例中,当热负荷占比增加、电负荷占比减小时,优化规划方案首先考虑投建电热锅炉,将电能转化为热能,实现多能互补共济;然后减少输电网容量、扩建热力管网;最后考虑增加CHP机组。
上述实施方式已经对本发明的一些细节进行了描述,但是不能理解为对本发明的限制,本领域的技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对其进行变化、修改、替换和变型。
Claims (3)
1.计及天然气管网和热力管网特性的综合能源系统优化方法,其特征在于,包括步骤:
1)构建含多能流耦合设备和储能设备的能源中心设备模型;
2)构建含电力网络、天然气管网和热力管网的能源网络模型;
3)以综合能源系统运行周期内总成本最小为优化目标,考虑综合能源系统建设约束和运行约束,建立计及天然气管网和热力管网特性的综合能源系统优化模型;
步骤1)中,所述能源中心设备模型抽象为一个输入-输出双端口网络模型,多种能流分别从两个端口输入与输出,多能流耦合设备和储能设备的输入和输出端按照能源形式,分别汇集到同一端点;
步骤1)中,所述的多能流耦合设备包括电热锅炉、燃气锅炉、燃气轮机和热电联产机组,其能源转移效率统一表示为:
所述储能设备包括储电、储气和储热设备,储能设备的运行约束统一表示为:
式中:下标t表示t时刻,为储能设备x的储能量;Pt κ,xi和Pt κ,xo分别为储能设备x的充、放能速率;ηκ,xi和ηκ,xo分别为储能设备x的充、放能效率;Δt为单位时间段的时长;和分别为储能设备x存储能量的上、下限;
所述能源中心设备模型两个端口的输入与输出功率需满足:
步骤2)中,采用直流潮流模型描述电力网络:
电力网络节点能量平衡约束表示为:
步骤2)中,所述的天然气管网中,天然气管道约束有:
根据气体状态方程与波义耳定律,管存相关的计算式如下:
且其满足质量守恒定律,如下式所示:
其中,
式中:Vij,t为天然气管道ij内的管存量;pi,t和pj,t分别为管道ij首末端的气压;和分别为管道ij的出、入口流量;和分别为管道ij内径和长度;Rgas为通用气体常数;为管道ij的管存系数;Mgas为天然气分子量;Tg、ψ和ρg分别为天然气温度、压缩因子和相对空气的密度;Δt为单位时间段的时长;
此外,天然气管道传输的气流量与首末端气压有关,实际运行中大多数输气管道以高雷诺数的流速运行即处于湍流状态,满足管道气流方程,如下式所示,参数均折算到标准状况下:
pi,min≤pi,t≤pi,max,
其中
所述的天然气管网中,加压站约束表示为:
pi,t≤ξcompj,t,
天然气网络节点能量平衡约束表示为:
式中:为天然气管网中与节点k相连的节点集合;和分别为管道jk出口和入口端的气功率;为外部气源注入综合能源系统的气功率;为注入能源中心的气功率;为天然气热值;和分别为管道ik的出、入口流量;ξcom表示加压站最大加压系数;
步骤2)中,所述的热力管网中,换热站约束有:
供水管和回水管的出入口温度约束表示如下:
热负荷和能源中心与换热站热交换约束表示如下:
热力管网节点热量守恒约束表示如下:
式中:和和分别为第k个能源中心/第f个热负荷的供水管与回水管的出入口温度;和分别为第k个能源中心和第f个热负荷与其换热站的热交换功率;cw为水的比热容;和分别为单位时间内流过换热站的工质质量;NZ为流入汇集点z的管道集合;Tz,t和分别为汇集点z和管道b出口的工质温度;为单位时间内管道b流出的工质质量;
所述的热力管网中,热网延时效应约束有:
其中,
式中:γb,t和φb,t分别为热力传输延时时长的上、下限;和分别为未计及温度损失时的管道出、入口温度;ρw为热力管网工质的密度;和分别为t-γb,t和t-φb,t+1时刻到t时刻注入管道的工质质量;Ν为正整数集合,n表示其中的元素;和分别为δt时间内流入和流出管道b的工质质量;和分别表示t-φb,t和t-γb,t时刻注入管道的工质温度;Ab和分别表示管道的横截面面积和长度;
所述的热力管网中,热网损耗约束有:
由于工质在传输过程中不可避免地与管道进行热交换而产生热量损耗,故管道出口温度根据苏霍夫温降公式进行修正:
其中,
所述的热力管网中,热网节点能量平衡约束有:
能源中心和热负荷的热能平衡约束:
步骤3)中,综合能源系统优化模型中,目标函数表示为:
其中,
式中:用下标s表示第s个场景;Cinv、和Ctotal分别表示计及设备残值的投资成本、第τ年的外部能源购买成本和系统运行周期内总成本;r为贴现率;Hor为规划年限;D为一年的天数;NS为一年内的场景集合;Νeh和Nbr分别为综合能源系统拓扑结构中的节点集合与支路集合;ΝX和Nnet分别为综合能源系统中的能源中心设备种类集合和能源网络种类集合;和为第k个能源中心中候选X类设备的集合和能源网络κ中节点i与j之间的候选线路或管道集合;ωs为场景s发生的概率;Φ为一个典型日的时间分段数;和分别为从外部购电和购气功率;和分别为电能和天然气的单位购买成本;假定投运都发生在年初,Rx、cx、βx和Sx分别为x的规划期末残值率、单位容量投资成本、候选设备投运状态和单台/条/回容量;Δt为单位时间段的时长;
假设能源中心设备与能源网络折旧程度与投运时间呈线性关系,x的残值率统一描述为:
步骤3)中,综合能源系统优化模型中,建设约束为:
综合能源系统的投资成本包括多能流耦合设备、储能设备以及电力网络、天然气管网和热力管网的建设成本,投资成本存在上限,如下式所示:
对于能源中心设备和能源网络,设备安装台数和线路或管道的建设条/回数需要满足下面的约束:
步骤3)中,综合能源系统优化模型中,运行约束为:
能源中心中的设备输入功率和爬/滑坡速度约束条件统一表示为:
在能源网络中,两节点之间建设多条并行线路,由于能源网络的非线性,需要分别计算每一条线路的运行状态,能源网络线路功率约束统一表示为:
式中:和分别为节点i与j之间第l条电力网线路的传输功率和天然气管道的出入口功率;为向热负荷f供热的第l回管道的传输功率;ζe,tran、ζg,tran和ζhex为候选电力线路、天然气管道和热力管道的容量裕度;0-1变量和为候选电力线路、天然气管道和热力管道的投运状态;和为候选电力线路、天然气管道和热力管道的容量;
从外部注入的电功率和气功率需要满足下面的约束:
2.根据权利要求1所述的计及天然气管网和热力管网特性的综合能源系统优化方法,其特征在于,利用增量法对非线性约束进行线性化处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910491982.1A CN110348602B (zh) | 2019-06-06 | 2019-06-06 | 计及天然气管网和热力管网特性的综合能源系统优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910491982.1A CN110348602B (zh) | 2019-06-06 | 2019-06-06 | 计及天然气管网和热力管网特性的综合能源系统优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110348602A CN110348602A (zh) | 2019-10-18 |
CN110348602B true CN110348602B (zh) | 2021-09-21 |
Family
ID=68181568
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910491982.1A Active CN110348602B (zh) | 2019-06-06 | 2019-06-06 | 计及天然气管网和热力管网特性的综合能源系统优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110348602B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111353128B (zh) * | 2020-01-17 | 2023-07-25 | 浙江工业大学 | 一种基于非合作博弈的多能量枢纽优化运行方法 |
CN111611690B (zh) * | 2020-04-17 | 2021-06-22 | 清华大学 | 一种用于综合能源网中热力管网运行参数的动态计算方法 |
CN111815111B (zh) * | 2020-06-02 | 2022-05-13 | 天津大学 | 一种考虑管线风险水平的区域综合能源扩展规划方法 |
CN111928294B (zh) * | 2020-08-06 | 2023-03-24 | 华能太原东山燃机热电有限责任公司 | 一种燃气-蒸汽联合循环机组热电成本分摊的方法 |
CN111882137B (zh) * | 2020-08-07 | 2022-10-11 | 西南石油大学 | 计及天然气管网压力能消纳的充电设施优化规划方法 |
CN112366697B (zh) * | 2020-10-30 | 2022-06-17 | 杭州意能电力技术有限公司 | 一种多能流配电网日前能量管理模型的管理方法 |
CN112580994A (zh) * | 2020-12-23 | 2021-03-30 | 华北电力大学 | 含分布式能源接入的园区综合能源系统规划方法 |
CN112989612B (zh) * | 2021-03-18 | 2023-09-19 | 贵州电网有限责任公司 | 一种基于麦考密克包络的电热综合能源系统线性建模方法 |
CN116843070B (zh) * | 2023-07-03 | 2024-01-26 | 上海轻环能源科技有限公司 | 电力现货市场天然气长输管网运行调度最优化方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106447529A (zh) * | 2016-08-30 | 2017-02-22 | 上海交通大学 | 考虑热水管网的分布式能源系统建模和运行优化方法 |
CN107808218A (zh) * | 2017-10-25 | 2018-03-16 | 国网天津市电力公司 | 基于热电比调节的城市能源互联网潮流计算方法 |
CN108258679A (zh) * | 2017-12-25 | 2018-07-06 | 国网浙江省电力有限公司经济技术研究院 | 考虑供热网储热特性的电-热综合能源系统优化调度方法 |
CN108596453A (zh) * | 2018-04-10 | 2018-09-28 | 山东大学 | 考虑网络动态特性的综合能源系统日前优化调度方法及系统 |
CN109255471A (zh) * | 2018-08-17 | 2019-01-22 | 国网山东省电力公司电力科学研究院 | 一种含风电的电-气-热综合能源系统扩展规划优化方法 |
CN109524957A (zh) * | 2018-11-07 | 2019-03-26 | 国网浙江省电力有限公司经济技术研究院 | 考虑碳交易机制和柔性负荷的综合能源系统优化调度方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040054564A1 (en) * | 2002-09-17 | 2004-03-18 | Fonseca Adolfo M. | Systems and methods for the optimization of resources in energy markets |
-
2019
- 2019-06-06 CN CN201910491982.1A patent/CN110348602B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106447529A (zh) * | 2016-08-30 | 2017-02-22 | 上海交通大学 | 考虑热水管网的分布式能源系统建模和运行优化方法 |
CN107808218A (zh) * | 2017-10-25 | 2018-03-16 | 国网天津市电力公司 | 基于热电比调节的城市能源互联网潮流计算方法 |
CN108258679A (zh) * | 2017-12-25 | 2018-07-06 | 国网浙江省电力有限公司经济技术研究院 | 考虑供热网储热特性的电-热综合能源系统优化调度方法 |
CN108596453A (zh) * | 2018-04-10 | 2018-09-28 | 山东大学 | 考虑网络动态特性的综合能源系统日前优化调度方法及系统 |
CN109255471A (zh) * | 2018-08-17 | 2019-01-22 | 国网山东省电力公司电力科学研究院 | 一种含风电的电-气-热综合能源系统扩展规划优化方法 |
CN109524957A (zh) * | 2018-11-07 | 2019-03-26 | 国网浙江省电力有限公司经济技术研究院 | 考虑碳交易机制和柔性负荷的综合能源系统优化调度方法 |
Non-Patent Citations (1)
Title |
---|
计及风电不确定性的电-气-热综合能源系统扩展规划方法;李哲,王成福,梁军,赵鹏辉,张哲;《电网技术》;20181130;第3477-3485页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110348602A (zh) | 2019-10-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110348602B (zh) | 计及天然气管网和热力管网特性的综合能源系统优化方法 | |
CN108596453B (zh) | 考虑网络动态特性的综合能源系统日前优化调度方法及系统 | |
Zhang et al. | Day-ahead optimal dispatch for integrated energy system considering power-to-gas and dynamic pipeline networks | |
CN110263966B (zh) | 考虑动态传热过程的电-热综合能源系统优化调度方法 | |
CN111815042B (zh) | 计及精细化热网模型的电热综合能源系统优化调度方法 | |
CN109389248B (zh) | 基于区域热网的综合能源系统热电功率协调调度方法 | |
CN112347607B (zh) | 一种基于凸松弛的热电联合调度方法 | |
CN111222257B (zh) | 基于凸优化的电热水多能流协同调度方法 | |
CN112668188B (zh) | 多园区综合能源系统的分布式鲁棒协同优化调度方法 | |
CN109359839A (zh) | 一种综合能源系统中区域热网的节点热价计算方法 | |
CN109447323A (zh) | 一种计及节点热价的综合能源系统两阶段容量配置方法 | |
Li et al. | Gradient descent iterative method for energy flow of integrated energy system considering multiple modes of compressors | |
Chen et al. | Optimal low‐carbon scheduling of integrated local energy system considering oxygen‐enriched combustion plant and generalized energy storages | |
CN116822683A (zh) | 一种基于碳交易计及不确定性的综合能源系统优化运行方法 | |
CN114139958A (zh) | 一种计及管存和热惯性的综合能源系统运行优化方法 | |
CN111724026B (zh) | 一种多能源网络与配水网耦合运行的优化方法 | |
CN116341881B (zh) | 计及热网灵活性的电-热系统鲁棒超前调度方法及系统 | |
CN112926835A (zh) | 一种考虑热网动态特性的综合能源系统优化调度方法 | |
CN110020506B (zh) | 基于电热型综合能源系统运行优化的差分格式选择方法 | |
CN116542447A (zh) | 一种电热系统优化调度方法 | |
CN113077173B (zh) | 一种基于综合能源系统节点能价策略的综合需求响应方法 | |
CN110728441B (zh) | 基于序列线性规划的电-气联合市场集中式出清方法 | |
CN112508730A (zh) | 一种综合能源虚拟电厂的风电消纳策略 | |
CN114066056A (zh) | 考虑热电协同综合能源系统灵活性的优化调度方法及系统 | |
CN115577479B (zh) | 一种区域冷、热、气碳素流计算模型的构建方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |