CN109726906B - 基于热网偏微分方程约束的热电联产系统日前调度方法 - Google Patents
基于热网偏微分方程约束的热电联产系统日前调度方法 Download PDFInfo
- Publication number
- CN109726906B CN109726906B CN201811568500.XA CN201811568500A CN109726906B CN 109726906 B CN109726906 B CN 109726906B CN 201811568500 A CN201811568500 A CN 201811568500A CN 109726906 B CN109726906 B CN 109726906B
- Authority
- CN
- China
- Prior art keywords
- heat
- equation
- formula
- supply network
- heat supply
- 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
- 238000005457 optimization Methods 0.000 claims abstract description 12
- 238000012546 transfer Methods 0.000 claims abstract description 10
- 230000001052 transient effect Effects 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims abstract description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 36
- 238000004146 energy storage Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 6
- 239000000446 fuel Substances 0.000 claims description 6
- 241000764238 Isis Species 0.000 claims description 4
- 238000009833 condensation Methods 0.000 claims description 4
- 230000005494 condensation Effects 0.000 claims description 4
- 238000011160 research Methods 0.000 claims description 4
- 238000004088 simulation Methods 0.000 claims description 3
- 238000005315 distribution function Methods 0.000 claims description 2
- 230000004907 flux Effects 0.000 claims description 2
- 230000003993 interaction Effects 0.000 claims description 2
- 238000012423 maintenance Methods 0.000 claims description 2
- 239000007789 gas Substances 0.000 description 5
- 238000010438 heat treatment Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 239000002918 waste heat Substances 0.000 description 1
Images
Abstract
本发明公开了一种基于热网偏微分方程约束的热电联产系统日前调度方法。本发明采用的技术方案为:首先基于热力学第一定律,建立能够反映热网暂态传热特性的偏微分方程模型;其次将上述热网偏微分方程进行差分处理,并将差分得到的代数方程组作为部分约束条件,以日运行成本最小为目标函数,建立热电联产系统日前优化调度模型;最后根据基于热网偏微分方程约束的热电联产系统日前调度模型,确定系统中各个设备的调度计划,联络线的计划功率,以及热网的运行状态。本发明所提出的优化调度方法能够反映热网的真实运行状态,并提升热电联产系统运行经济性。
Description
技术领域
本发明属于能源系统运行优化技术领域,涉及一种基于热网偏微分方程约束的热电联产系统日前调度方法。
背景技术
热电联产系统是综合能源系统的一种表现形式,以热电联产机组为核心,通过对发电后余热进行回收利用实现电能和热能的联合供应。热电联产系统基于热能的综合梯级利用,一次能源利用率可超过60%,远高于燃煤电厂。目前,北欧部分国家的热电联产机组已经占据了火电总机组数目的60%以上;我国80%的工业供热和30%的民用供热也均来自热电联产机组。
热电联产是解决我国城市供热热源结构不合理、热电供需矛盾突出、供热热源能效低污染重等问题的主要途径之一。国家发改委、国家能源局等五部门联合下发了《关于印发<热电联产管理办法>的通知》(发改能源〔2016〕617号),对我国热电联产行业的各个环节进行了规范,并提出力争实现北方大中型以上城市热电联产集中供热率达到60%以上的目标。积极发展热电联产有利于提高能源综合利用效率、防治大气污染问题、缓和区域性用电用热矛盾等,是构建低碳、高效、节能、环保的能源供应体系的重要手段之一。
区域热网是连接热电厂与热用户之间的媒介,通常呈辐射状。我国城市热网集中供热半径可达数十公里,供热规模大,传输距离远,传热特性复杂,难以精细化调度。目前已有诸多文献研究区域热网的传输时延与传热特性,但多采用数值仿真或者实验测试方法,模型复杂且难以应用于系统的优化调度;部分关于热电联产系统优化调度的研究中所用热网模型则过于简化,难以反映热网的暂态传热特性,优化结果与实际运行状态差距大,缺乏实用性。
发明内容
为真实反映热网运行状态,并最大限度地减少热电联产系统的运行成本,提升热电联产系统经济性,本发明提供一种基于热网偏微分方程约束的热电联产系统日前调度方法。
本发明采用的技术方案如下:基于热网偏微分方程约束的热电联产系统日前调度方法,其包括:
步骤1,基于热力学第一定律,建立能够反映热网暂态传热特性的偏微分方程模型;
步骤2,将上述热网偏微分方程进行差分处理,并将差分得到的代数方程组作为部分约束条件,以日运行成本最小为目标函数,建立热电联产系统日前优化调度模型;
步骤3,根据基于热网偏微分方程约束的热电联产系统日前调度模型,确定系统中各个设备的调度计划,联络线的计划功率,以及热网的运行状态,所述各个设备包括抽汽冷凝式热电机组、发电机、电锅炉、燃气锅炉、热泵和储能电站。
本发明具有的有益效果如下:本发明所提出的基于热网偏微分方程约束的优化调度方法,充分考虑了电、热子系统的传输网络特性,能够反映热网的真实运行状态,最大限度提升热电联产系统运行的经济性。
附图说明
图1是本发明的系统结构图;
图2是本发明的流程图;
图3是本发明的电、热功率调度结果图;
图4是本发明的供水管网节点温度运行状态图;
图5是本发明的回水管网节点温度运行状态图。
具体实施方式
构成本发明的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。
如图1所示,本实施例的热电联产系统由联络线、抽汽冷凝式热电机组、发电机、电锅炉、燃气锅炉、热泵、风机和储能电站组成。该系统中,热电联产机组承担供热基荷,电功率缺额由发电机、电网、风机和储能电站补偿,热功率缺额由电锅炉、燃气锅炉和热泵补偿。系统通过一条联络线与其他电网相连,可从电网购电也可向电网售电,系统同时向电、热两种负荷提供能量需求,调度中心通过输入的预测信息、价格信息、室温约束和系统运行约束进行统一优化,然后将优化运行指令下发给各设备执行。
如图2所示,本实施的一种基于热网偏微分方程约束的热电联产系统日前调度方法,具体包括如下步骤:
步骤1:基于热力学第一定律,建立能够反映热网暂态传热特性的偏微分方程模型;
步骤2:将上述热网偏微分方程进行差分处理,并将差分得到的代数方程组作为部分约束条件,以日运行成本最小为目标函数,建立热电联产系统日前优化调度模型;
步骤3:根据基于热网偏微分方程约束的热电联产系统日前调度模型,确定系统中各个设备的调度计划,联络线的计划功率,以及热网的运行状态,所述各个设备包括:抽汽冷凝式热电机组、电锅炉和发电机。
步骤1中所述的反映热网暂态传热特性的偏微分方程模型为:
所述热网由多根相互连接的长直管道构成,取其中一根长直热水管道作为研究对象,将热源端设为原点位置,沿着水流方向设为x轴正方向,建立一维坐标系。
步骤1.1:假设水力工况稳定,对位于x处长为dx的管道微元应用热力学第一定律可得:
dU=dQ+dW (1)
式中:dU为该微元的内能增量;dQ为外界向该微元传递的热量;dW为外界对该微元做的功。忽略管壁的摩擦热影响,dW取零,式(1)可以写为:
式中:为前一个微元通过热传导传递给该微元的热量;为该微元通过热传导传递给下一个微元的热量;为前一个微元内热水流入该微元时带来的热量;为该微元内热水流向下一个微元时带走的热量;为该微元通过管壁向外界散失的热量。
步骤1.2:该微元的内能增量为:
式中:qx和qx+dx分别为管道x处和x+dx处的热流密度;A为热水横截面积。由热传导定律可知,热流密度与温度梯度成正比,但方向相反,故qx和qx+dx可由式(5)计算:
式中:λw为管道内热水的导热系数。
式中:Tx和Tx+dx分别为管道x处和x+dx处的热水温度,满足式(7):
式中:Ta为管壁外环境温度;R为管壁单位长度热阻。
步骤1.6:将各物理量代入式(2)化简可得描述管道暂态传热特性的控制方程:
式中:ρ为热水的密度。
为求解该控制方程,设置如下初边值条件:
步骤2中所述的热电联产系统日前优化调度模型的具体建立步骤如下:
步骤2.1:将上述热网偏微分方程进行差分处理;
将区域Γ={(x,t)|0≤x≤L,0≤t≤P}用两簇平行直线
分割成矩形网络,其中L和P分别为管道长度和仿真时长;h和τ分别为空间步长和时间步长,满足:
在节点(xi,tk)处考虑所述偏微分方程,有:
对于其中的偏导数,由下式差商代替:
略去小项:
并用Ti k代替T(xi,tk),得到如下差分格式:
步骤2.2:分析所选差分格式的稳定性;
若式(19)成立,
则误差满足式(20):
式中:||·||∞为L∞范数。
||εk+1||∞≤||εk||∞,0≤k≤N-1 (21)
递推可得:
||εk+1||∞≤||εk||∞,0≤k≤N-1 (22)
即式(19)成立时差分格式关于L∞范数是稳定的。
步骤2.3:分析所选差分格式的收敛性;
由式(15)可知:
式中:
若式(19)成立,则误差满足式(26)
||εk+1||∞≤||εk||∞+C1τ(τ+h),0≤k≤N-1 (27)
递推可得:
||εk||∞≤||ε0||∞+C1kτ(τ+h)≤C1P(τ+h),1≤k≤N (28)
即式(19)成立时差分格式在L∞范数下关于时间步长是1阶、空间步长是1阶收敛的。
步骤2.4:建立目标函数,以日运行费用的最小化为优化目标,具体目标函数为:
步骤3根据日前优化调度的结果,系统直接获得抽汽冷凝式热电机组、发电机、电锅炉、燃气锅炉、热泵和储能电站的调度计划,联络线的计划功率,以及热网的运行状态。
如图3所示,展示了系统电功率和热功率的日前调度计划,可以看出,由于考虑了供热系统的热惯性,热功率出力与热负荷不必实时匹配,而是围绕热负荷小范围波动。
如图4、图5所示,展示了热网供、回水管道中三个节点的温度变化情况,可以看出,供水管网节点温度在90~100℃,回水管网温度在65~75℃,符合热网运行约束。
另外需要说明的是,在上述具体实施方式中所描述的各个具体技术特征,在不矛盾的情况下,可以通过任何合适的方式进行组合。为了避免不必要的重复,本发明对各种可能的组合方式不再另行说明。
Claims (3)
1.基于热网偏微分方程约束的热电联产系统日前调度方法,其特征在于,包括如下步骤:
步骤1,基于热力学第一定律,建立能够反映热网暂态传热特性的偏微分方程模型;
步骤2,将上述热网偏微分方程进行差分处理,并将差分得到的代数方程组作为部分约束条件,以日运行成本最小为目标函数,建立热电联产系统日前优化调度模型;
步骤3,根据基于热网偏微分方程约束的热电联产系统日前调度模型,确定系统中各个设备的调度计划,联络线的计划功率,以及热网的运行状态,所述各个设备包括抽汽冷凝式热电机组、发电机、电锅炉、燃气锅炉、热泵和储能电站;
所述步骤1的具体内容包括:
所述热网由多根相互连接的长直管道构成,取其中一根长直热水管道作为研究对象,将热源端设为原点位置,沿着水流方向设为x轴正方向,建立一维坐标系;
步骤1.1:假设水力工况稳定,对位于x处长为dx的管道微元应用热力学第一定律得:
dU=dQ+dW (1)
式中:dU为该微元的内能增量;dQ为外界向该微元传递的热量;dW为外界对该微元做的功;忽略管壁的摩擦热影响,dW取零,式(1)写为:
式中:为前一个微元通过热传导传递给该微元的热量;为该微元通过热传导传递给下一个微元的热量;为前一个微元内热水流入该微元时带来的热量;为该微元内热水流向下一个微元时带走的热量;为该微元通过管壁向外界散失的热量;
步骤1.2:该微元的内能增量为:
式中:qx和qx+dx分别为管道x处和x+dx处的热流密度;A为热水横截面积;由热传导定律知,热流密度与温度梯度成正比,但方向相反,故qx和qx+dx由式(5)计算:
式中:λw为管道内热水的导热系数;
式中:Tx和Tx+dx分别为管道x处和x+dx处的热水温度,满足式(7):
式中:Ta为管壁外环境温度;R为管壁单位长度热阻;
步骤1.6:将各物理量代入式(2)化简得描述管道暂态传热特性的控制方程:
式中:ρ为热水的密度;
为求解该控制方程,设置如下初边值条件:
2.根据权利要求1所述的基于热网偏微分方程约束的热电联产系统日前调度方法,其特征在于,所述热电联产系统日前优化调度模型的具体建立步骤如下:
步骤2.1:将上述热网偏微分方程进行差分处理;
将区域Γ={(x,t)|0≤x≤L,0≤t≤P}用两簇平行直线
分割成矩形网络,其中L和P分别为管道长度和仿真时长;h和τ分别为空间步长和时间步长,满足:
在节点(xi,tk)处考虑所述偏微分方程,有:
对于其中的偏导数,由下式差商代替:
略去小项:
并用Ti k代替T(xi,tk),得到如下差分格式:
步骤2.2:分析所选差分格式的稳定性;
若式(19)成立,
则误差满足式(20):
式中:||·||∞为L∞范数;
||εk+1||∞≤||εk||∞,0≤k≤N-1 (21)
递推得:
||εk+1||∞≤||εk||∞,0≤k≤N-1 (22)
即式(19)成立时差分格式关于L∞范数稳定;
步骤2.3:分析所选差分格式的收敛性;
由式(15)知:
式中:
若式(19)成立,则误差满足式(26)
||εk+1||∞≤||εk||∞+C1τ(τ+h),0≤k≤N-1 (27)
递推得:
||εk||∞≤||ε0||∞+C1kτ(τ+h)≤C1P(τ+h),1≤k≤N (28)
即式(19)成立时差分格式在L∞范数下关于时间步长是1阶、空间步长是1阶收敛;
步骤2.4:建立目标函数,以日运行费用的最小化为优化目标,具体目标函数为:
3.根据权利要求1所述的基于热网偏微分方程约束的热电联产系统日前调度方法,其特征在于:根据日前优化调度的结果,热电联产系统直接获得抽汽冷凝式热电机组、发电机、电锅炉、燃气锅炉、热泵和储能电站的调度计划,联络线的计划功率以及热网运行状态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811568500.XA CN109726906B (zh) | 2018-12-21 | 2018-12-21 | 基于热网偏微分方程约束的热电联产系统日前调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811568500.XA CN109726906B (zh) | 2018-12-21 | 2018-12-21 | 基于热网偏微分方程约束的热电联产系统日前调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109726906A CN109726906A (zh) | 2019-05-07 |
CN109726906B true CN109726906B (zh) | 2021-08-27 |
Family
ID=66297803
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811568500.XA Active CN109726906B (zh) | 2018-12-21 | 2018-12-21 | 基于热网偏微分方程约束的热电联产系统日前调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109726906B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110661266B (zh) * | 2019-09-24 | 2022-07-15 | 东南大学 | 用于热电联产系统动态最优能流计算的差分步长优选方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108491973A (zh) * | 2018-03-21 | 2018-09-04 | 浙江大学 | 一种考虑需求侧响应的电-热联合系统优化调度方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104616208B (zh) * | 2015-02-04 | 2017-10-13 | 东南大学 | 一种基于模型预测控制的冷热电联供型微电网运行方法 |
CN105279709A (zh) * | 2015-11-17 | 2016-01-27 | 东南大学 | 一种基于热网热惯性的电网日前优化调度方法 |
CN106505634B (zh) * | 2016-12-14 | 2019-03-19 | 东南大学 | 基于两阶段协调优化与控制的冷热电联供型微网运行方法 |
CN106786793B (zh) * | 2016-12-14 | 2019-04-09 | 东南大学 | 一种基于鲁棒优化的冷热电联供型微网运行方法 |
CN108599137B (zh) * | 2017-12-28 | 2020-03-31 | 东南大学 | 考虑区域热网暂态传热特性的多能流系统优化运行方法 |
-
2018
- 2018-12-21 CN CN201811568500.XA patent/CN109726906B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108491973A (zh) * | 2018-03-21 | 2018-09-04 | 浙江大学 | 一种考虑需求侧响应的电-热联合系统优化调度方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109726906A (zh) | 2019-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Huang et al. | Heat and power load dispatching considering energy storage of district heating system and electric boilers | |
Zhou et al. | An engineering approach to the optimal design of distributed energy systems in China | |
Sarbu et al. | A review of modelling and optimisation techniques for district heating systems | |
Tian et al. | System modeling and optimal dispatching of multi-energy microgrid with energy storage | |
Chen et al. | Optimal operation of integrated energy system based on exergy analysis and adaptive genetic algorithm | |
Lin et al. | A multi-agent based optimization architecture for energy hub operation | |
Jayasuriya et al. | Energy hub modelling for multi-scale and multi-energy supply systems | |
Tan et al. | Modeling and daily operation optimization of a distributed energy system considering economic and energy aspects | |
Li et al. | Economic dispatch of an integrated heat-power energy distribution system with a concentrating solar power energy hub | |
CN109726906B (zh) | 基于热网偏微分方程约束的热电联产系统日前调度方法 | |
CN113886761A (zh) | 一种综合能源系统能效分析评估方法 | |
Werner et al. | District heating research in China | |
Wang et al. | Investigating the potential for district heating networks with locally integrated solar thermal energy supply | |
Kapitonov et al. | Principles regulation of electricity tariffs for the integrated generation of traditional and alternative energy sources | |
Zhou et al. | Digital Twin Modeling and Simulation of Distributed and Centralized Integrated Energy System on CloudPSS-IESLab | |
Rudra et al. | Techno-economic analysis of a local district heating plant under fuel flexibility and performance | |
Park et al. | Short-term forecasting of city gas daily demand | |
Zhao et al. | Optimal scheduling method for electrical-thermal integrated energy system considering heat storage characteristics of heating network | |
Chen et al. | A multi-energy microgrid modelling and optimization method based on exergy theory | |
Liu et al. | A method to calculate energy station’s output in community energy planning considering the attenuation and delay of pipe network | |
Long et al. | Optimization strategy of CCHP integrated energy system based on source-load coordination | |
Zhao et al. | A multi-objective power flow optimization model of electric and thermal distribution network using the power flow method | |
CN107451929B (zh) | 一种微型燃气轮机冷热电联供微电网的经济性分析系统 | |
Yao et al. | A transient thermodynamic model of district heating network for operational optimization of the energy integration system | |
Wrzalik | Smart district heating networks in the era of energy transformation |
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 |