基于热网偏微分方程约束的热电联产系统日前调度方法技术方案

本发明专利技术公开了一种基于热网偏微分方程约束的热电联产系统日前调度方法。本发明专利技术采用的技术方案为：首先基于热力学第一定律，建立能够反映热网暂态传热特性的偏微分方程模型；其次将上述热网偏微分方程进行差分处理，并将差分得到的代数方程组作为部分约束条件，以日运行成本最小为目标函数，建立热电联产系统日前优化调度模型；最后根据基于热网偏微分方程约束的热电联产系统日前调度模型，确定系统中各个设备的调度计划，联络线的计划功率，以及热网的运行状态。本发明专利技术所提出的优化调度方法能够反映热网的真实运行状态，并提升热电联产系统运行经济性。

Day-ahead scheduling method for cogeneration system based on partial differential equation constraints of heat network

【技术实现步骤摘要】
基于热网偏微分方程约束的热电联产系统日前调度方法
本专利技术属于能源系统运行优化
，涉及一种基于热网偏微分方程约束的热电联产系统日前调度方法。
技术介绍
热电联产系统是综合能源系统的一种表现形式，以热电联产机组为核心，通过对发电后余热进行回收利用实现电能和热能的联合供应。热电联产系统基于热能的综合梯级利用，一次能源利用率可超过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：该微元的内能增量为：式中：c为管道内热水比热容；为热水的质量流量；v为热水的流速；T为热水温度。步骤1.3：所述的和可由式(4)计算：式中：qx和qx+dx分别为管道x处和x+dx处的热流密度；A为热水横截面积。由热传导定律可知，热流密度与温度梯度成正比，但方向相反，故qx和qx+dx可由式(5)计算：式中：λw为管道内热水的导热系数。步骤1.4：所述的和可由式(6)计算：式中：Tx和Tx+dx分别为管道x处和x+dx处的热水温度，满足式(7)：步骤1.5：所述的可由式(8)计算：式中：Ta为管壁外环境温度；R为管壁单位长度热阻。步骤1.6：将各物理量代入式(2)化简可得描述管道暂态传热特性的控制方程：式中：ρ为热水的密度。为求解该控制方程，设置如下初边值条件：式中：为初始时刻管道内各处的水温分布函数；ψ(t)为管道热水流入端的水温变化与时间的函数关系。步骤2中所述的热电联产系统日前优化调度模型的具体建立步骤如下：步骤2.1：将上述热网偏微分方程进行差分处理；将区域Γ＝{(x,t)|0≤x≤L,0≤t≤P}用两簇平行直线分割成矩形网络，其中L和P分别为管道长度和仿真时长；h和τ分别为空间步长和时间步长，满足：在节点(xi,tk)处考虑所述偏微分方程，有：对于其中的偏导数，由下式差商代替：略去小项：并用Tik代替T(xi,tk)，得到如下差分格式：步骤2.2：分析所选差分格式的稳定性；设为式(17)所示方程组的解，则误差满足式(18)所示方程：若式(19)成立，则误差满足式(20)：式中：||·||∞为L∞范数。又注意到且有||εk+1||∞≤||εk||∞,0≤k≤N-1(21)递推可得：||εk+1||∞≤||εk||∞,0≤k≤N-1(22)即式(19)成立时差分格式关于L∞范数是稳定的。步骤2.3：分析所选差分格式的收敛性；由所选差分格式可知误差满足式(23)：由式(15)可知：式中：若式(19)成立，则误差满足式(26)再注意到且有||εk+1||∞≤||εk||∞+C1τ(τ+h),0≤k≤N-1(27)递推可得：||εk||∞≤||ε0||∞+C1kτ(τ+h)≤C1P(τ+h),1≤k≤N(28)即式(19)成立时差分格本文档来自技高网...

【技术保护点】
1.基于热网偏微分方程约束的热电联产系统日前调度方法，其特征在于，包括如下步骤：步骤1，基于热力学第一定律，建立能够反映热网暂态传热特性的偏微分方程模型；步骤2，将上述热网偏微分方程进行差分处理，并将差分得到的代数方程组作为部分约束条件，以日运行成本最小为目标函数，建立热电联产系统日前优化调度模型；步骤3，根据基于热网偏微分方程约束的热电联产系统日前调度模型，确定系统中各个设备的调度计划，联络线的计划功率，以及热网的运行状态，所述各个设备包括抽汽冷凝式热电机组、发电机、电锅炉、燃气锅炉、热泵和储能电站。

【技术特征摘要】
1.基于热网偏微分方程约束的热电联产系统日前调度方法，其特征在于，包括如下步骤：步骤1，基于热力学第一定律，建立能够反映热网暂态传热特性的偏微分方程模型；步骤2，将上述热网偏微分方程进行差分处理，并将差分得到的代数方程组作为部分约束条件，以日运行成本最小为目标函数，建立热电联产系统日前优化调度模型；步骤3，根据基于热网偏微分方程约束的热电联产系统日前调度模型，确定系统中各个设备的调度计划，联络线的计划功率，以及热网的运行状态，所述各个设备包括抽汽冷凝式热电机组、发电机、电锅炉、燃气锅炉、热泵和储能电站。2.根据权利要求1所述的基于热网偏微分方程约束的热电联产系统日前调度方法，其特征在于，所述步骤1的具体内容包括：所述热网由多根相互连接的长直管道构成，取其中一根长直热水管道作为研究对象，将热源端设为原点位置，沿着水流方向设为x轴正方向，建立一维坐标系；步骤1.1：假设水力工况稳定，对位于x处长为dx的管道微元应用热力学第一定律得：dU＝dQ+dW(1)式中：dU为该微元的内能增量；dQ为外界向该微元传递的热量；dW为外界对该微元做的功；忽略管壁的摩擦热影响，dW取零，式(1)写为：式中：为前一个微元通过热传导传递给该微元的热量；为该微元通过热传导传递给下一个微元的热量；为前一个微元内热水流入该微元时带来的热量；为该微元内热水流向下一个微元时带走的热量；为该微元通过管壁向外界散失的热量；步骤1.2：该微元的内能增量为：式中：c为管道内热水比热容；为热水的质量流量；v为热水的流速；T为热水温度；步骤1.3：所述的和由式(4)计算：式中：qx和qx+dx分别为管道x处和x+dx处的热流密度；A为热水横截面积；由热传导定律知，热流密度与温度梯度成正比，但方向相反，故qx和qx+dx由式(5)计算：式中：λw为管道内热水的导热系数；步骤1.4：所述的和由式(6)计算：式中：Tx和Tx+dx分别为管道x处和x+dx处的热水温度，满足式(7)：步骤1.5：所述的由式(8)计算：式中：Ta为管壁外环境温度；R为管壁单位长度热阻；步骤1.6：将各物理量代入式(2)化简得描述管道暂态传热特性...

【专利技术属性】
技术研发人员：李志浩，汪湘晋，张雪松，倪筹帷，赵波，顾伟，姚帅，陆帅，潘光胜，吴晨雨，
申请(专利权)人：国网浙江省电力有限公司电力科学研究院，国家电网有限公司，
类型：发明
国别省市：浙江,33