CN109241676B - 综合能源系统中热网及建筑物时间分辨率的选择方法 - Google Patents
综合能源系统中热网及建筑物时间分辨率的选择方法 Download PDFInfo
- Publication number
- CN109241676B CN109241676B CN201811213952.6A CN201811213952A CN109241676B CN 109241676 B CN109241676 B CN 109241676B CN 201811213952 A CN201811213952 A CN 201811213952A CN 109241676 B CN109241676 B CN 109241676B
- Authority
- CN
- China
- Prior art keywords
- supply network
- heat
- heat supply
- time resolution
- simulation model
- 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
- 238000000034 method Methods 0.000 title claims abstract description 13
- 238000004088 simulation Methods 0.000 claims abstract description 28
- 238000005457 optimization Methods 0.000 claims abstract description 27
- 238000002347 injection Methods 0.000 claims abstract description 7
- 239000007924 injection Substances 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 4
- 239000007789 gas Substances 0.000 description 19
- 230000005611 electricity Effects 0.000 description 9
- 238000003860 storage Methods 0.000 description 9
- 238000010438 heat treatment Methods 0.000 description 7
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- 230000006872 improvement Effects 0.000 description 4
- 238000005338 heat storage Methods 0.000 description 3
- 239000003345 natural gas Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000007599 discharging Methods 0.000 description 2
- 238000004146 energy storage Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005485 electric heating Methods 0.000 description 1
- 230000017525 heat dissipation Effects 0.000 description 1
- 238000004519 manufacturing process 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
- 230000005855 radiation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明专利公开了一种综合能源系统中热网及建筑物时间分辨率的选择方法,建立综合能源系统运行优化模型及仿真模型,以运行优化结果中热网注入功率作为热网及热负荷状态仿真模型的输入,得到仿真模型结果,再建立热网及建筑物状态误差指标和时间分辨率选择判据,通过比较运行优化结果与仿真模型结果,判断误差指标值是否满足初设的判据,当不满足判据时则调整热网及热负荷时间分辨率并重复以上步骤,直到得到合适的热网及热负荷时间分辨率,本方法充分考虑不同时间分辨率对热网及热负荷运行状态的影响,选择出最佳的热网及热负荷时间分辨率,保证综合能源系统的热功率平衡与供热质量,保证综合能源系统的安全可靠运行。
Description
所属领域
本发明属于综合能源系统领域,具体涉及一种综合能源系统中热网及建筑物时间分辨率的选择方法。
背景技术
能源与环境问题是当前人类社会所面临的两大关键问题。能源系统是维持人类生产与生活的基础,构建一个绿色、可持续的能源系统是解决能源与环境问题的关键。与传统独立能源系统相比,综合能源系统融合了电、气、冷、热等多种能源,对多种能源进行协同优化,因此具有高效的能源利用效率与良好的排放特性,对实现人类社会的可持续发展具有重要意义。
运行优化是保证综合能源系统高效运行的前提。由于综合能源系统中热网及热负荷存在热惯性,因此需要使用(偏)微分方程进行描述,在运行优化中如何选择恰当的热网及热负荷时间分辨率是一个关键问题。时间分辨率过小则导致优化问题规划太大,难以求解;时间分辨率过大则无法精确描述热网及热负荷状态,导致热功率失衡及室内温度越限等问题,从而威胁综合能源系统的安全可靠运行。
发明内容
本发明正是针对现有技术中的问题,提供了一种综合能源系统中热网及建筑物时间分辨率的选择方法,建立综合能源系统运行优化模型及仿真模型,以运行优化结果中热网注入功率作为热网及热负荷状态仿真模型的输入,得到仿真模型结果,再建立热网及建筑物状态误差指标和时间分辨率选择判据,通过比较运行优化结果与仿真模型结果,判断误差指标值是否满足初设的判据,当不满足判据时则调整热网及热负荷时间分辨率并重复以上步骤,直到得到合适的热网及热负荷时间分辨率,本方法充分考虑不同时间分辨率对热网及热负荷运行状态的影响,选择出最佳的热网及热负荷时间分辨率,保证综合能源系统的热功率平衡与供热质量,保证综合能源系统的安全可靠运行。
为了实现上述目的,本发明采用的技术方案是:综合能源系统中热网及建筑物时间分辨率的选择方法,包括以下步骤:
S1,建立综合能源系统运行优化模型,并求解运行优化结果;
S2,建立热网及热负荷状态仿真模型,将步骤S1运行优化结果中的热电联产系统向热网注入热媒的温度作为输入,求解仿真模型;
S3,建立热网及热负荷状态误差指标,并计算误差指标值;
S4,建立热网及热负荷时间分辨率选择判据;
S5,比较步骤S1中运行优化结果与步骤S2中仿真模型结果,判断误差指标值是否满足步骤S4中的判据。
作为本发明的一种改进,所述步骤S5中,选定热网及热负荷初始时间分辨率,判断误差指标值是否满足步骤S4中的判据,若不满足,则调整时间分辨率后重复步骤S1-S5。
作为本发明的一种改进,所述步骤S3中热网及热负荷状态误差指标包括热功率误差百分比和最大室温偏差。
作为本发明的又一种改进,所述步骤S3中热功率误差百分比△Q%为:
其中,td为调度时段;th,s为及热网及热负荷仿真模型时间段;Δtd为综合能源系统调度指令周期;△th,s为热网及热负荷仿真模型的时间步长;为运行优化模型和热网及热负荷仿真模型中热电联产系统向热网注入的热功率结果;
作为本发明的另一种改进,所述步骤S4建立的热网及热负荷时间分辨率的判据为:
其中δ1和δ2分别为所设置的最大容许热功率误差与最大容许室温偏差。
与现有技术相比,本发明专利提出了一种综合能源系统运行优化中热网及热负荷时间分辨率选择方法,针对综合能源系统中热网及热负荷的热惯性特征,充分考虑不同时间分辨率对热网及热负荷运行状态的影响,选择出合适的热网及热负荷时间分辨率,保证综合能源系统的热功率平衡与供热质量,从而保证综合能源系统安全可靠运行。。
附图说明
图1是本发明热网管段建模示意图;
图2是本发明建筑物热力学模型;
图3是本发明综合能源系统热网及热负荷时间分辨率选择流程图;
图4是本发明实施例2中综合能源系统结构图。
具体实施方式
以下将结合附图和实施例,对本发明进行较为详细的说明。
实施例1
一种综合能源系统中热网及建筑物时间分辨率的选择方法,如图3所示,包括以下步骤:
S1,建立综合能源系统运行优化模型,所述步骤进一步包括:
S11,建立综合能源系统运行优化目标函数:
所述目标函数包括综合能源系统购买天然气成本、综合能源系统从电网购电成本和向电网售电收益:
其中,Δtd为综合能源系统调度指令周期;td为调度时段;Td为td集合;cgas、cgrid,s分别为天然气热值价格、向电网购电价格及向电网售电价格; 分别为第td时段燃气轮机输出电功率、燃气锅炉输出热功率、向电网购电功率及向电网售电功率;ηgt、ηgb分别为燃气轮机及燃气锅炉效率。
S12,建立综合能源系统运行优化约束条件,所述步骤进一步包括:
S121,建立热电联产系统运行约束条件,至少包括能量平衡约束、设备上下线约束及储能装置约束:
设备上下线约束包括燃气轮机出力约束、电网交互功率约束、可再生能源出力约束、以及燃气锅炉热出力约束。燃气轮机出力约束为其中 P gt分别为燃气轮机电功率上下限;分别为燃气轮机的状态及电出力;电网交互功率约束为
储能装置约束至少包括蓄电池运行约束、储热罐运行约束及储冷罐运行约束,所述蓄电池运行约束为
S122,建立热网及热负荷约束条件:所述热网约束条件包括管段延时与热损约束、节点能量平衡约束及源节点及负荷节点交换热功率:
热网建模方法如图1所示,首先按照下式计算热网中的管段相关参数。其中,γj为j中所包含的最大“能量块”数量;为取整函数;为管段j的热媒质量流量,kg/s;Δth为模型分辨率,h;ρw为密度,kg/m3;为j的截面积,m2;lj为j的长度,m;Rj为j所包含所有“能量块”的热媒质量,kg;,Φpipe为集合。
管段延时与热损约束分别为:
其中,分别为第th时段供水管段j和回水管段j的进口温度;kj为第th-γj时段“能量块”占出口处“能量块”的比例系数;分别为未考虑热损情况下第th时段出口处热媒温度;为第th时段管网环境温度;分别为第th时段出口处热媒温度;βj为管段j的热损系数;kj和βj计算公式为
其中,λj为j的传热系数,kW/(m·℃);cw为的比热容,kJ/(kg·℃)。
节点能量平衡约束为
源节点及负荷节点交换热功率计算公式为
热负荷考虑为建筑物供暖负荷,其热力学模型如图2所示,建筑物一阶热力学模型如下式所示:
其中,τin、τout分别为建筑物室内温度和室外温度;Rs为建筑物热阻;Cair为建筑物热容;Ql为向建筑物注入的热功率。对上式差分可得:
S123,建立热电联产系统与热网耦合约束:
热电联产系统与热网耦合约束为:
S2,建立热网及热负荷状态仿真模型:
其中,分别为第th时段供水管段j和回水管段j的进口温度;kj为第th-γj时段“能量块”占出口处“能量块”的比例系数;分别为未考虑热损情况下第th时段出口媒温度。βj为管段j的热损系数,λj为管段j的传热系数;cw为热媒的比热容;分别为第th时段节点n处的注入热功率和热负荷热功率;Φns、Φnl分别为源节点集合及负荷节点集合;分别为第th时段建筑物室内温度和室外温度;Rs为建筑物热阻;Cair为建筑物热容;为第th时段向建筑物注入的热功率;γ为散热片辐射系数;F为散热器总散热面积。
S3,建立热网及热负荷状态误差指标:
热负荷状态误差指标包括热功率误差百分比和最大室温偏差,计算公式分别为:
其中,△Q%为热功率误差百分比;td为调度时段;th,s为及热网及热负荷仿真模型时间段;Δtd为综合能源系统调度指令周期;△th,s为热网及热负荷仿真模型的时间步长;为运行优化模型和热网及热负荷仿真模型中热电联产系统向热网注入的热功率结果;为最大室温偏差;为th时段室温偏差;为热网及热负荷仿真结果中室内温度结果; τ in分别为允许的室内温度最大值与最小值。
S4,建立热网及热负荷时间分辨率选择判据:
热网及热负荷时间分辨率选择判据为:
其中δ1和δ2分别为所设置的最大容许热功率误差与最大容许室温偏差。
S5,建立热网及热负荷时间分辨率选择流程,选定热网及热负荷初始时间分辨率,设定最大容许热功率误差δ1和最大容许室温偏差δ2,比较运行优化模型结果与状态仿真模型结果,判断是否满足S4中的判据,若不满足,则调整时间分辨率后重复步骤S1-S5,若满足,则步骤结束,获得最佳时间分辨率。
实施例2
以图4所示的综合能源系统为例,该系统包含1台5MW的燃气轮机,1台5MW燃气锅炉,容量为1.5MW的风机,1MWh的蓄电池和5MWh的储热罐。热网包含6个节点,其中节点1连接CHP系统,节点4、5和6分别连接热负荷。运行优化周期为24h。最大容许热功率误差百分比设置为1.5%,最大容许室温偏差设置为0.5℃。初始热网及建筑物时间分辨率Δth设置为60min,Δth取值集合设置为{60min,30min,20min,15min,10min,1min}。按照图3所示的流程图执行。不同时间分辨率下热功率误差、室外偏差、以及运行成本如表1所示。可见随着时间分辨率的减小,热功率误差百分比和室温偏差均减小。当模型分辨率取60min和30min时,热功率误差与室温偏差均不满足要求,而当模型分辨率取20min时,热功率误差及室温偏差均满足要求。因此,适应于本例的热网及热负荷时间分辨率为20min。此外,不同时间分辨率下系统运行成本并无明显差异,因此可认为时间分辨率基本不影响系统运行经济性。
表1不同时间分辨率下的指标与运行成本
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实例的限制,上述实例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等同物界定。
Claims (3)
1.综合能源系统中热网及建筑物时间分辨率的选择方法,其特征在于,包括以下步骤:
S1,建立综合能源系统运行优化模型,并求解运行优化结果;
S2,建立热网及热负荷状态仿真模型,将步骤S1运行优化结果中的热电联产系统向热网注入热媒的温度作为输入,求解仿真模型;
S3,建立热网及热负荷状态误差指标,并计算误差指标值,热网及热负荷状态误差指标包括热功率误差百分比和最大室温偏差,其中,热功率误差百分比△Q%为:
其中,td为调度时段;th,s为及热网及热负荷仿真模型时间段;Δtd为综合能源系统调度指令周期;△th,s为热网及热负荷仿真模型的时间步长;为运行优化模型和热网及热负荷仿真模型中热电联产系统向热网注入的热功率结果;
S4,建立热网及热负荷时间分辨率选择判据;
S5,比较步骤S1中运行优化结果与步骤S2中仿真模型结果,判断误差指标值是否满足步骤S4中的判据。
2.如权利要求1所述的综合能源系统中热网及建筑物时间分辨率的选择方法,其特征在于所述步骤S5中,选定热网及热负荷初始时间分辨率,判断误差指标值是否满足步骤S4中的判据,若不满足,则调整时间分辨率后重复步骤S1-S5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811213952.6A CN109241676B (zh) | 2018-10-18 | 2018-10-18 | 综合能源系统中热网及建筑物时间分辨率的选择方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811213952.6A CN109241676B (zh) | 2018-10-18 | 2018-10-18 | 综合能源系统中热网及建筑物时间分辨率的选择方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109241676A CN109241676A (zh) | 2019-01-18 |
CN109241676B true CN109241676B (zh) | 2023-04-07 |
Family
ID=65053799
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811213952.6A Active CN109241676B (zh) | 2018-10-18 | 2018-10-18 | 综合能源系统中热网及建筑物时间分辨率的选择方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109241676B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110020506B (zh) * | 2019-03-08 | 2023-04-07 | 东南大学 | 基于电热型综合能源系统运行优化的差分格式选择方法 |
CN112364493A (zh) * | 2020-10-30 | 2021-02-12 | 国网福建省电力有限公司厦门供电公司 | 城市综合能源系统智能调度方法及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108599137B (zh) * | 2017-12-28 | 2020-03-31 | 东南大学 | 考虑区域热网暂态传热特性的多能流系统优化运行方法 |
CN108288855B (zh) * | 2018-01-23 | 2020-11-24 | 国电南京自动化股份有限公司 | 一种基于微电网的动态策略仿真与优化方法 |
-
2018
- 2018-10-18 CN CN201811213952.6A patent/CN109241676B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109241676A (zh) | 2019-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109190785B (zh) | 一种电热耦合综合能源系统运行优化方法 | |
CN111144668B (zh) | 考虑场景模拟的综合能源系统随机优化模型的建立方法 | |
CN107808216B (zh) | 电-气-热互联系统弃风弃光和电气热负荷削减综合最小优化模型构建方法 | |
CN110188492B (zh) | 一种考虑热网特性的冷热电联供微网优化调度方法 | |
CN111950807B (zh) | 计及不确定性与需求响应的综合能源系统优化运行方法 | |
CN111428351B (zh) | 基于前推回代法的电-热综合能源系统潮流计算方法 | |
CN109034508B (zh) | 考虑电热双重不确定性的综合能源系统鲁棒优化调度方法 | |
CN111030094B (zh) | 一种多能互补园区供能可靠性评估方法 | |
CN109389248B (zh) | 基于区域热网的综合能源系统热电功率协调调度方法 | |
CN112417651B (zh) | 一种基于后悔规避的用户级综合能源系统优化方法 | |
CN111191353B (zh) | 考虑热惯性的电-热综合能源系统协调运行方法 | |
CN109241676B (zh) | 综合能源系统中热网及建筑物时间分辨率的选择方法 | |
CN112036652A (zh) | 基于机会约束规划的光伏-储能综合能源系统规划方法 | |
CN115186902A (zh) | 温室综合能源系统的调控方法、装置、终端及存储介质 | |
CN110729720A (zh) | 一种蓄热罐容量规划方法和电转气参数敏感性分析方法 | |
CN111476394B (zh) | 一种适用于电热气等多能源系统的鲁棒运行优化方法 | |
CN113255198A (zh) | 一种含虚拟储能的冷热电联供微电网多目标优化方法 | |
CN112949940A (zh) | 一种综合能源虚拟电厂多设备选址方法 | |
CN116432824A (zh) | 基于多目标粒子群的综合能源系统优化方法及系统 | |
CN113011655B (zh) | 一种基于两阶段随机规划的社区多能源系统的混合储能系统规划方法 | |
CN107832873A (zh) | 基于双层母线式结构的综合能源系统优化规划方法及装置 | |
Li et al. | Planning model of integrated energy system considering P2G and energy storage | |
CN112926835A (zh) | 一种考虑热网动态特性的综合能源系统优化调度方法 | |
CN112001639A (zh) | 综合能源系统能源需求的可调能力评估方法及存储介质 | |
Yu et al. | Wind-CHP generation aggregation with storage capability of district heating network |
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 |