【技术实现步骤摘要】
一种考虑区域热网暂态传热特性的多能流系统优化运行方法
本专利技术涉及一种多能流系统优化运行方法，具体属于能源系统运行优化
技术介绍
随着世界能源危机和环境污染问题的日益加重，能源互联网、综合能源系统、“互联网+”智慧能源等创新想法不断涌现，转变能源供给模式，建设清洁高效的现代能源体系具有重要意义。多能流系统依托能源传输技术、能量转换设备和可再生能源等相关技术的不断革新，耦合电网、热网和天然气网等多种能源网络，有利于实现多能协同供应和能源综合梯级利用，逐渐成为能源领域的重要发展趋势之一。多种能源相互耦合是多能流系统区别于传统能源系统的重要标志之一，各类能源的特性差异对系统的规划运行提出了新的挑战。区域热网作为传输热能的媒介，由于在传输损耗和时间尺度方面与电网具有明显差异，在构建多能协同供应的多能流系统中吸引了广泛的研究关注。目前对于热网水力工况的研究已趋于成熟，但对热力工况的研究尚处于探索阶段。部分学者仅关注热网的传热特性，使用数值模拟、实验建模等方法，虽然能够获得较高精度的热力数据，但是所用模型较为复杂，求解耗时长，并且需要给定用户热负荷数据和热源供水温...
【技术保护点】
1.一种考虑区域热网暂态传热特性的多能流系统优化运行方法，其特征在于，该运行方法包括以下步骤：步骤1)：基于图论建立考虑回水管网拓扑结构、管段温度变化动态过程和管段传输时延的区域热网暂态传热特性模型，该模型包括六个方面内容：热源特性、热负荷特性、节点流量平衡、节点功率融合特性、管段传热特性和管段暂/稳态传热特性判据；步骤2)：建立多能流系统运行优化模型，以系统一个滚动调度周期内的运行成本最低为目标函数，结合系统应该满足的运行约束条件，确定多能流系统中各设备的实时出力、热网各热源的实时供回水温度，实现系统的统一调度和优化运行。
【技术特征摘要】
1.一种考虑区域热网暂态传热特性的多能流系统优化运行方法，其特征在于，该运行方法包括以下步骤：步骤1)：基于图论建立考虑回水管网拓扑结构、管段温度变化动态过程和管段传输时延的区域热网暂态传热特性模型，该模型包括六个方面内容：热源特性、热负荷特性、节点流量平衡、节点功率融合特性、管段传热特性和管段暂/稳态传热特性判据；步骤2)：建立多能流系统运行优化模型，以系统一个滚动调度周期内的运行成本最低为目标函数，结合系统应该满足的运行约束条件，确定多能流系统中各设备的实时出力、热网各热源的实时供回水温度，实现系统的统一调度和优化运行。2.按照权利要求1所述的考虑区域热网暂态传热特性的多能流系统优化运行方法，其特征在于：所述的步骤1)中，将热水管网看作流体网络，每根管道为一条支路，热源、热负荷和管道的连接点为节点；建立区域热网暂态传热特性模型的过程为：步骤101)考虑热源特性：设TSO和TSI分别为热源的供、回水温度，为了保证热源的正常供热，需要对热源的供回水温度进行限制，如式(1)所示：式中，TSO,max和TSO,min分别表示热源供水温度的上、下限；TSI,max和TSI,min分别表示热源回水温度的上、下限；另一方面，热源均由物理设备构成，应满足功率爬坡约束，如式(2)所示：式中，Qs(t)表示当前调度时刻点热源的输出热功率；Qs(t+Δt)表示下一个调度时刻点热源的输出热功率；Δt表示调度时间间隔；和分别表示热源爬坡热功率的上、下限；步骤102)考虑热负荷特性：对于包含热负荷的热网支路，供、回水管道的温度与热负荷的取用功率满足式(3)：式中：TLI和TLO分别为所有热负荷处的供水温度和回水温度构成的列向量；c为热媒比热容，取4.2kJ/(kg·℃)；ρ为热媒密度，取934.667kg/m3；ηHX为热网负荷端各换热设备的平均效率；QL/F为所有热负荷的取用功率与管道流量组成的列向量，如式(4)所示：QL/F＝[QL1/fL1,QL2/fL2,…QLk/fLk]T式(3)式中：QLi和fLi分别表示第i个负荷取用的功率和其所在管道的运行流量，i＝1,2,…,k，k为负荷总数；为了保证热负荷的供热质量，需要对热负荷处的供回水温度进行限制，如式(5)所示：式中，TLI和TLO分别为热负荷处的供、回水温度；TLI,max和TLI,min分别表示热负荷供水温度的上、下限；TLO,max和TLO,min分别表示热负荷回水温度的上、下限；步骤103)考虑节点流量平衡：对于热网中任一节点，流入的热媒流量之和等于流出的热媒流量之和，即:式中，和分别表示热网的上、下关联矩阵；F为流量列向量，上、下关联矩阵的定义如式(7)所示：热网流量列向量的定义如式(8)所示：F＝[f1,f2,…,fb]T式(7)式中：fi表示第i条支路的热媒流量，i＝1,2,…,b，b为支路总数，n为节点总数；步骤104)考虑节点功率融合特性：由能量守恒定律可知，对于热网中的任一节点，流入的功率之和等于流出的功率之和，即：式中：TSF和TEF是由管段始、末温度和流量组成的列向量，定义如式(10)：式中：Tsi和Tei分别表示第i根管道的起始温度和末端温度；节点处发生功率融合后，该节点的温度与流出该节点的热媒温度相等，即：Tni＝Ts1,i＝Ts2,i＝…＝Tsk,i式(10)式中：Tni为第i个节点的温度；Ts1,i,Ts2,i,…,Tsk,i为所有与第i个节点直接相连，并且热媒流出该节点的管道的起始温度；步骤105)考虑管段传热特性，当管段的传输时延远小于系统调度策略的时间尺度时，在每个调度时刻点处，热网已经达到稳态，使用稳态传热公式描述热网的传热特性；而当管段的传输时延与系统调度策略的时间尺度相当时，使用暂态传热公式描述热网的传热特性，以保证调度策略的精确性；具体包括步骤1051)～1052)：步骤1051)管段稳态传热特性：对于一维管道，由傅里叶定律知：传导的热功率与温度梯度和传热面积成正比，即：式中：Qloss代表管道的热功率损失，λ表示管道的热导率；A表示传热的正对面积；T表示温度；δ表示厚度；等号右边的负号表示热流方向总是与温度梯度方向相反；流体在管道内的热功率损失采用多层圆筒壁热传导模型计算：式中：Ta为外界环境温度，r为圆筒的半径，ri为第i层圆筒的半径，λi为第i层材料的热导率，dl为管道微元的长度；记管道单位长度热阻R为：则单位长度管道内流体损失的热功率为：对于一根管段，其首端的热功率为：Q0＝cρf(Ts-Tr)式(15)式中：c表示热媒的比热容；ρ表示热媒的密度；f表示热媒的运行流量；Ts表示热源的供水温度；Tr表示回水温度；距离首端x处的热网传输热功率等于热源处热功率与热网传输损失功率之差，即：解得管道内流体温度随管长分布为：考虑到式(18)中指数项x的系数接近于零，将其泰勒展开并略去二次项及更高次项，得到线性化后的形式：将式(19)改成矩阵形式：TE'＝ηTs'式(19)式中：TE'和Ts'分别是由管段首端温度、x处温度和外界温度构成的等效管段起始温度列向量和等效管段末端温度列向量；η是由管段参数构成的矩阵：TS'＝[Ta-Ts1,Ta-Ts2,…,Ta-Tsb]TTE'＝[Te1-Ts1,Te2-Ts2,…,Teb-Tsb]T步骤1052)管段暂态传热特性：取管段内长度为dl的流体微元进行分析，不计流体微元之间的相互作用，仅考虑流体微元通过管壁对外界散失热量Hloss，由傅里叶定律可知：式中：T(t)表示该流体微元在t时刻的温度；初始时刻，该流体微元刚从热源处发出，具有的能量为：式中：d为管段的内径；t时刻该流体微元具有的能量等于初始能量减去散失的热量，即：联立式(23)和(24)可以解得管段内该流体微元的温度随时间的变化关系，如式(25)所示：若该流体微元经过时...
【专利技术属性】
技术研发人员：姚帅，顾伟，陆帅，吴晨雨，
申请(专利权)人：东南大学，
类型：发明
国别省市：江苏,32
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。