一种计及延迟的区域蒸汽热网拓扑结构优化方法及系统技术方案

本发明专利技术申请公开了一种计及延迟的区域蒸汽热网拓扑结构优化方法及系统，本模型预测控制方法包括如下步骤：步骤S1，建立基于Modelica语言的蒸汽管道仿真模型；步骤S2，计算管道末端蒸汽温度随源端流量、温度变化的响应特性；步骤S3，设定管网最低延迟目标函数，建立区域源荷空间分布约束，采用粒子群优化算法进行目标函数的求解；步骤S4：将基于Modelica语言的模型搭建、管道末端温度对源端参数变化的响应特性计算与基于粒子群优化算法的热网拓扑规划分别进行模块封装并建立数据传输通道。

A topology optimization method and system of regional steam heating network considering delay

【技术实现步骤摘要】
一种计及延迟的区域蒸汽热网拓扑结构优化方法及系统
本专利技术涉及智能供热领域，特别涉及一种计及延迟的区域蒸汽热网拓扑结构优化方法及系统。
技术介绍
蒸汽管道是用于蒸汽供热系统对蒸汽进行运输的载体。在一个复杂多源环状的蒸汽管网中，往往具有多种不同的蒸汽供应方案，而且蒸汽在输送过程中状态参数变化大，甚至伴随相变发生，所以蒸汽管道模型在热力学、传热学和流体力学等方程的求解方面存在难题。目前,在对蒸汽管网进行设计和运行调节计算时，水力计算和热力计算是各自分开的。一般来说，这种简化方法将蒸汽处理为单相可压缩流体，认为密度只是压力的函数，这种人为地将压力和焓值分开求解的离式算法给蒸汽管网的计算带来了较大误差。此外，对蒸汽管道流动的模拟计算多是基于稳态过程的仿真模型，无法准确预测和分析蒸汽管道流动的瞬变特性，即蒸汽在管道的流动过程中温度、压力和密度随时间的变化特性。考虑到蒸汽管网工质流动过程的蒸汽参数波动、可压缩性、状态变化、摩擦和传热等多种因素的作用，如何从机理角度出发，建立复杂多源环状蒸汽热网的数学模型，并基于该模型建立区域蒸汽热网的拓扑结构优化方法是智慧供热背景下高效利用区域现有资源禀赋的关键技术问题。本专利技术针对上述问题，提出了考虑末端响应的热网机理模型，建立了基于粒子群优化算法的区域蒸汽热网拓扑结构优化方法，并将二者按照统一的系统框架进行集成。
技术实现思路
本专利技术的目的是提供一种计及延迟的区域蒸汽热网拓扑结构优化方法及系统，以基于Modelica语言建模的仿真模型和基于粒子群算法的优化手段，对区域内给定热源、热用户分布约束条件下的热网规划方案进行寻优，从而以较低的人工成本获取管网输运延迟最低的热网拓扑规划。一种计及延迟的区域蒸汽热网拓扑结构优化方法及系统，包含如下步骤：步骤S1，建立基于Modelica语言的蒸汽管道仿真模型；步骤S2，计算管道末端蒸汽温度随源端流量、温度变化的响应特性；步骤S3，设定管网最低延迟目标函数，建立区域源荷空间分布约束，采用粒子群优化算法进行目标函数的求解；步骤S4：将基于Modelica语言的模型搭建、管道末端温度对源端参数变化的响应特性计算与基于粒子群优化算法的热网拓扑规划分别置于三个相应模块，并建立三者间数据通信，以便优化过程调取模型计算数据。进一步地，建立基于Modelica语言的蒸汽管道仿真模型，即通过Modelica语言的管道接口、流动模型和传热模型组成具有分布质量、能量和动量平衡的直管模型；所述蒸汽管道仿真模型含括的质量连续方程、能量平衡方程、管道摩擦阻力方程和动量守恒方程的偏微分方程式表示如下：所述质量连续方程：式中：x为管道的沿程距离(流量沿坐标x)，m；t为时间；ρ为ρ(x,t)表示蒸汽密度，kg/m3；v为v(x,t)表示蒸汽流速，m/s；A、A(x)为蒸汽流通面积，m2；所述能量平衡方程：式中：T为T(x,t)表示蒸汽温度，K；u为u(x,t)表示比内能，J；z为z(x)表示地面高度，m；g为重力m/s2；所述动量守恒方程：式中：FF为阻力，N；所述管道摩擦阻力方程：式中：f为摩擦系数；S为周长，m；节点与外界环境的换热方程式：q＝α(Ts-Ta)(5)式中：q表示换热量，W/m2；Ts为蒸汽温度，K；Ta为蒸汽温度，K；α为传热系数，W/m2·K；所述质量连续方程、能量平衡方程、管道摩擦阻力方程和动量守恒方程均通过有限体积法和动量平衡的交错网格方案处理计算。进一步地，计算管道末端蒸汽温度随源端流量、温度变化的响应特性，具体流程如下：步骤S201，向Modelica源测模型输入供给蒸汽的温度和流量，向末端模型输入用户的压力需求，并对保温层、管道厚度的影响因素作相应的传热系数转换，以获取管道内部的温度、压力、和流量分布；步骤S202，在Modelica模型中保持源侧蒸汽参数不变，让管道质量流量以阶跃的方式增大，待末端蒸汽温度稳定后，又以阶跃的方式减小到原来的流量水平，从而获得该工况下温度响应特性；步骤S203，在Modelica模型中保持源侧蒸汽压力与流量不变，令源端供应蒸汽温度以阶跃的方式增大，待末端蒸汽温度稳定后，又以阶跃的方式减小到原来的供应温度水平，从而计算该工况下的温度响应完成时间。进一步地，设定管网最低延迟的目标函数，建立区域源荷空间分布约束，采用粒子群优化算法进行目标函数的求解，求解过程中每次计算都需要获取Modelica模型所求解的延迟特性数据；具体流程包括：步骤S301，以热网内n个热用户归一化后的热力延迟时间为目标函数：式中：DTi和Dqi分别为第i个用户末端温度随源端温度、流量变化的响应时间，s；DTmax和Dqmax分别为所有用户中温度、流量响应的最大值，s；DTmin和Dqmin分别为所有用户中温度、流量响应的最小值，s；k为流量响应相对于温度响应的重要性系数，一般取1；步骤S302，建立源端和末端用户空间分布约束；考虑到实际中热源和用户一般位于公路旁边，而蒸汽管道也往往沿公路建设，将源荷间的蒸汽管线连接方案集定义为以公路拓扑为网格，以热源、热用户为节点的网格路径搜索集，通过递归的方式搜寻出区域内源荷所有的连接方案，并以此方案集作为优化问题的可行域；步骤S303，采用粒子群优化算法进行目标函数的求解，其具体求解过程为：步骤S303.1，将可行域内的具体连接方案组合为d维粒子，并初始化设置惯性因子、加速因子等参数；步骤S303.2，初始化管网连接方案的速度向量，将每个方案的个体历史最优位置设置为当前方案的位置，并计算所有连接规划方案的最优位置；步骤S303.3，更新每个调度方案的速度和位置，更新公式如下：vij(t+1)＝vij(t)+c1r1(pbestij(t)-xij(t))+c2r2(gbestj(t)-xij(t))(7)xij(t+1)＝xij(t)+vij(t+1)(8)其中，i表示第i个粒子；j表示粒子的第j维；vij(t)表示粒子i在进化到t代时的第j维飞行速度分量；xij(t)表示粒子i在进化到t代时的第j维位置分量；pbestij(t)表示粒子i在进化到t代时的第j维个体最优位置pbesti分量；gbestj(t)表示进化到t代时整个粒子群的最优位置gbest的第j维分量；c1，c2为加速因子或称学习因子；r1，r2为[0,1]的随机数；步骤S303.4，计算每个方案位置的目标函数值，更新每个方案的个体历史最优位置与所有调度方案的最优位置；步骤S303.5，若达到最大迭代次数则停止计算，输出结果；否则返回步骤S303.3继续搜索计算。进一步地，将Modelica的热力水力建模、末端用户热力响应时间计算与基于粒子群优化算法的热网拓扑规划分别进行模块封装，建立数据传输通道，以便优化规划模块每次迭代都能提取经建本文档来自技高网...

【技术保护点】
1.一种计及延迟的区域蒸汽热网拓扑结构优化方法，其特征在于，包含如下步骤：/n步骤S1，建立基于Modelica语言的蒸汽管道仿真模型；/n步骤S2，计算管道末端蒸汽温度随源端流量、温度变化的响应特性；/n步骤S3，设定管网最低延迟目标函数，建立区域源荷空间分布约束，采用粒子群优化算法进行目标函数的求解；/n步骤S4：将基于Modelica语言的模型搭建、管道末端温度对源端参数变化的响应特性计算与基于粒子群优化算法的热网拓扑规划分别置于三个相应模块，并建立三者间数据通信，以便优化过程调取模型计算数据。/n

【技术特征摘要】
1.一种计及延迟的区域蒸汽热网拓扑结构优化方法，其特征在于，包含如下步骤：
步骤S1，建立基于Modelica语言的蒸汽管道仿真模型；
步骤S2，计算管道末端蒸汽温度随源端流量、温度变化的响应特性；
步骤S3，设定管网最低延迟目标函数，建立区域源荷空间分布约束，采用粒子群优化算法进行目标函数的求解；
步骤S4：将基于Modelica语言的模型搭建、管道末端温度对源端参数变化的响应特性计算与基于粒子群优化算法的热网拓扑规划分别置于三个相应模块，并建立三者间数据通信，以便优化过程调取模型计算数据。


2.根据权利要求1所述的一种计及延迟的区域蒸汽热网拓扑结构优化方法，其特征在于，所述步骤S1，建立基于Modelica语言的蒸汽管道仿真模型，即通过Modelica语言的管道接口、流动模型和传热模型组成具有分布质量、能量和动量平衡的直管模型；所述蒸汽管道仿真模型含括的质量连续方程、能量平衡方程、管道摩擦阻力方程和动量守恒方程的偏微分方程式表示如下：
所述质量连续方程：



式中：x为管道的沿程距离(流量沿坐标x)，m；t为时间；ρ为ρ(x,t)表示蒸汽密度，kg/m3；v为v(x,t)表示蒸汽流速，m/s；A、A(x)为蒸汽流通面积，m2；
所述能量平衡方程：



式中：T为T(x,t)表示蒸汽温度，K；u为u(x,t)表示比内能，J；z为z(x)表示地面高度，m；g为重力m/s2；
所述动量守恒方程：



式中：FF为阻力，N；
所述管道摩擦阻力方程：



式中：f为摩擦系数；S为周长，m；
节点与外界环境的换热方程式：
q＝α(Ts-Ta)(5)
式中：q表示换热量，W/m2；Ts为蒸汽温度，K；Ta为蒸汽温度，K；α为传热系数，W/m2·K；
所述质量连续方程、能量平衡方程、管道摩擦阻力方程和动量守恒方程均通过有限体积法和动量平衡的交错网格方案处理计算。


3.根据权利要求1所述的一种计及延迟的区域蒸汽热网拓扑结构优化方法，其特征在于，所述步骤S2，计算管道末端蒸汽温度随源端流量、温度变化的响应特性，具体流程如下：
步骤S201，向Modelica源测模型输入供给蒸汽的温度和流量，向末端模型输入用户的压力需求，并对保温层、管道厚度的影响因素作相应的传热系数转换，以获取管道内部的温度、压力、和流量分布；
步骤S202，在Modelica模型中保持源侧蒸汽参数不变，让管道质量流量以阶跃的方式增大，待末端蒸汽温度稳定后，又以阶跃的方式减小到原来的流量水平，从而获得该工况下温度响应特性；
步骤S203，在Modelica模型中保持源侧蒸汽压力与流量不变，令源端供应蒸汽温度以阶跃的方式增大，待末端蒸汽温度稳定后，又以阶跃的方式减小到原来的供应温度水平，从而计算该工况下的温度响应完成时间。


4.根据权利要求1所述的一种计及延迟的区域蒸汽热网拓扑结构优化方法，其特征在于，所述步骤S3，设定管网最低延迟的目标函数，建立区域源荷空间分布约束，采用粒子群优化算法进行目标函数的求解，求解过程中每次计算都需要获取Modelica模型所求解的延迟特性数据；具体流程包括：
步骤S301，以热网内n个热用户归一化后的热力延迟时间为目标函数：



式中：DTi和Dqi分别为第i个用户末端温度随源端温度、流量变化的响应时间，s；DTmax和Dqmax分别为所有用户中温度、流量响应的最大值，s；DTmin和Dqmin分别为所有用户中温度、流量响应的最小值，s；k为流量响应相对于温度响应的重要性系数，一般取1；
步骤S302，建立源端和末端用户空间分布约束；考虑到实际中热源和用户一般位于公路旁边，而蒸汽管道也往往沿...

【专利技术属性】
技术研发人员：陆海，严玉廷，苏适，李耀华，杨洋，李翔，张旭东，栾思平，
申请(专利权)人：云南电网有限责任公司电力科学研究院，云南电网有限责任公司大理供电局，
类型：发明
国别省市：云南;53