工业园区综合能源系统蒸汽混合建模仿真研究.pdf
《工业园区综合能源系统蒸汽混合建模仿真研究.pdf》由会员分享,可在线阅读,更多相关《工业园区综合能源系统蒸汽混合建模仿真研究.pdf(8页珍藏版)》请在咨信网上搜索。
1、第 52 卷 第 8 期 Vol.52 No.8 2023 年 8 月 THERMAL POWER GENERATION Aug.2023 修 回 日 期:2022-12-07 基 金 项 目:国家自然科学基金项目(51806190);国家重点研发计划项目(2019YFE0126000)Supported by:National Natural Science Foundation of China(51806190);National Key Research and Development Program(2019YFE0126000)第一作者简介:周懿(1983),女,博士,助理研究员,
2、主要研究方向为供热管网稳动态分析与运行调控,。通信作者简介:林小杰(1992),男,博士,副研究员,主要研究方向为能源系统调度优化,。DOI:10.19666/j.rlfd.202212199 工业园区综合能源系统蒸汽 混合建模仿真研究 周 懿1,2,李建训3,王家乐4,林小杰2,4,钟 崴2,4(1.常州英集动力科技有限公司,江苏 常州 213022;2.浙江大学常州工业技术研究院,江苏 常州 213000;3.招远玲珑热电有限公司,山东 招远 265400;4.浙江大学能源工程学院,浙江 杭州 310027)摘要针对目前工业园区综合能源系统中蒸汽建模缺乏动态建模与分析以及求解时采用数值法需
3、要在精度和计算效率间平衡的问题,基于蒸汽流动基本方程,建立蒸汽流动过程的水力热力过程数学模型。通过解析与数值混合的方法进行水力热力建模,并采用有限差分法求解热力模型。将该方法与商业仿真软件仿真数据相比:求解误差在 0.43%左右,同时探究了不同求解条件对于热力过程的影响;发现空间步长增大会使温度曲线的波动性增加;随着管道内径增大,温度传递变慢,模型响应速度下降;初始压力对模型的求解过程和精度并无太大影响。关键词蒸汽系统;动态模型;热力模型;有限差分法;向后差分 引用本文格式周懿,李建训,王家乐,等.工业园区综合能源系统蒸汽混合建模仿真研究J.热力发电,2023,52(8):121-128.ZH
4、OU Yi,LI Jianxun,WANG Jiale,et al.Simulation study on steam hybrid modeling of integrated energy system in industrial parksJ.Thermal Power Generation,2023,52(8):121-128.Simulation study on steam hybrid modeling of integrated energy system in industrial parks ZHOU Yi1,2,LI Jianxun3,WANG Jiale4,LIN Xi
5、aojie2,4,ZHONG Wei2,4(1.Changzhou Engipower Technology Co.,Ltd.,Changzhou 213022,China;2.Changzhou Industrial Technology Research Institute of Zhejiang University,Changzhou 213000,China;3.Linglong Group ThermoElectron Corporation,Zhaoyuan 265400,China;4.College of Energy Engineering,Zhejiang Univers
6、ity,Hangzhou 310027,China)Abstract:In view of the current lack of dynamic modeling and analysis of steam modeling in the integrated energy system of industrial parks,the numerical method needs to balance between accuracy and computational efficiency.Based on the basic equation of steam flow,this pap
7、er establishes a mathematical model of the hydrothermal process of steam flow.Hydrothermal model is carried out by using a mixture of analytical and numerical methods,and the thermal model is solved by finite difference method for simulation.Compared with the simulation data of commercial simulation
8、 software,the solution error of this method is about 0.43%,and the influence of different solution conditions on the thermal process is explored.It was found that an increase in spatial step size would increase the volatility of the temperature curve.As the inner diameter of the pipeline increases,t
9、he temperature transfer slows down and the response speed of the model decreases;The initial pressure does not have a significant impact on the solving process and accuracy of the model.Key words:steam system;dynamic model;thermal model;finite difference method;backward differential 目前综合能源研究以工业园区综合能
10、源为主,电和热是综合能源研究主要方向,其中热网蒸汽规划设计及优化调度已成为研究的热点1-5。热网中蒸汽的输送十分关键,工业用供热系统的 2/3 来自蒸122 2023 年 http:/ 汽,但是蒸汽系统建模往往被忽略,目前石化领域考虑了蒸汽建模,但其建模方法较为独特难以整合进入综合能源系统。所以本文对热网中蒸汽的输送过程开展机理建模,为热网管网的有效规划及综合成本最小化打下理论基础。田兆斐等6基于质量守恒、能量守恒定律建立了两相箱体模型和流体网络模型。沈智鹏等7采用混合建模法建立了大型集装箱船舶辅助锅炉的炉膛、汽包等数学模型。刘笑驰等8利用半隐式差分法实现蒸汽管网水力、热力耦合的稳态求解。孙玉
11、宝等9基于节点方程法求解节点温度与节点压力,建立了蒸汽稳定流动数学模型。张增刚10基于节点方程法,反复迭代计算求解压力和温度。现有热力管网模型多为稳态模型,谢冬雷等11基于管道压力变化公式、节点流量方程和节点压力连续方程,建立了管网水力计算模型。章金华12建立了稳态水力计算模型,给出了压力、摩擦阻力、密度计算公式和流程图。郭永强13优化了管道压力变化公式,并提升了模型精度。李祥立等14基于分段积分法对蒸汽管网水力计算的逐段法作了改进,考虑了蒸汽密度变化对比摩阻的影响,建立了稳态水力模型。但是蒸汽在管网中的流动往往处于不稳定工作状态,其参数持续随时间发生变化,所以建立蒸汽管网的动态模型尤为必要1
12、5。刘斯斌等16利用 APROS软件研究了长输蒸汽供热管网动态延迟性及其影响因素。陆海17和孙鑫南18等利用 MODELICA 软件研究了单根蒸汽管道动态仿真模型及其延迟评估。高建强等19借助 Gensystem 模型开发平台,构建了管网仿真动态模型。周守军20基于节点和支路理论,建立了动态水力模型,通过改进平方根法求解。Chen等人21建立了动态水力模型,采用欧拉有限差分技术,引入虚拟节点求解。宋扬等22采用稀疏矩阵法求解蒸汽供热管网的动态水力热力模型,减少了求解时间,加快了求解速度。胡宪法等15建立了动态热力模型,准确定量刻画管网的蓄热、传热、耗散等关键过程。殷戈等23以单管热力计算的迭代
13、算法为基础,采用二叉树的结构,编制了支状管网的热力设计计算软件。王东24利用相关的热力学公式建立了一套数学模型,由此简化了水力计算方法。本文针对单管中蒸汽的输送过程,拟从蒸汽动态建模出发,将水力、热力模型分步求解,采用有限差分法中的显式差分格式对热力模型进行求解,重点分析其输送过程的动态热力特性,以及离散步长、管道参数、初始压力对热力模型的影响。1 单管动态物理模型 1.1 模型假设 本文针对单管蒸汽输送过程建立动态水力热力模型,实现对模型的求解和分析,从而得到管道参数的变化特性,以便蒸汽的合理调度和利用。针对单管中过热蒸汽进行动态建模,首先,对蒸汽输送过程进行机理建模,假设如下:1)由于管道
14、约束的存在,将蒸汽流动简化为沿管道的一维方向流动;2)假定蒸汽为理想气体,在压力较低、密度较低的情况下(分子间力可以忽略不计),气体可视为理想气体;3)蒸汽在管道内不发生相变,为单相可压缩气体流动过程;4)管道为水平管段,不考虑管道倾斜和垂直 管段;5)不考虑管道热损失。1.2 蒸汽流动水力热力模型 根据上述假设,基于蒸汽在单管中流动的质量守恒方程、流动方程和能量守恒方程,并加入气体状态方程和气体焓方程,建立单管中蒸汽流动的动态数学模型:22200202,vtxvvPvtxxDvhvuxtp Thh p T (1)式中:t 为时间,s;x 为空间距离,m;为蒸汽密度,kg/m3;v为蒸汽的流速
15、,m/s;p 为蒸汽压力,Pa;为管道的沿程摩阻系数;D 为管道的内径,m;h 为蒸汽比焓,J/kg;u 为蒸汽比内能,J/kg;K 为传热系数,W/(m2 K);T 为蒸汽温度,K。2 单管动态物理模型求解 本文中水力模型和热力模型之间存在耦合关系,动态水力模型和热力模型的耦合会大大增加模型求解过程中的计算难度和计算量,由于工业园区第 8 期 周 懿 等 工业园区综合能源系统蒸汽混合建模仿真研究 123 http:/ 综合能源系统热力学温度变化(5%10%)较小25,且目前有关动态建模的文献也提出了水力与热力分离求解的思路15,22,所以本文采用平均温度的方法解耦,最终将水力模型与热力模型分
16、离求解。即采用解析与数值结合的方法,利用平均温度求解水力模型,使用水力模型的求解结果利用数值法对热力模型进行求解。2.1 水力模型求解 蒸汽在管道中的一维流动的动态水力模型为:22002vtxvvpvtxxD (2)其中,当蒸汽流速远小于声速时,忽略对流项,即:20vx (3)并对摩擦力项中的速度平方项进行线性化处理:a2a2222vvvvDD (4)式中:va为管道中蒸汽流速的基值。根据假设,蒸汽符合理想气体状态方程:pRT (5)式中:R 为气体常数,J/(kg K)。根据假设,蒸汽为单相流体,不发生相变,可得:GvA (6)式中:G 为蒸汽的质量流量,kg/s;A 为管道面积,m2。将方
17、程转化到频域求解,得蒸汽管道末端的流量和压力值分别为:gg ggg22g gggglgg ggg4000222400g0g222g2e2242e224KZ YKllKZ YKllK GY pGGKZ YK GY pGKZ Y (7)gg ggg2g g2g4 00022gglgg gggg24000222g g2e2242e224KZ YKllKZ YKllK pZ GppKZ YKZ GpKZ Yp (8)式中:G1为管道出口流量,kg/s;G0为管道入口流量,kg/s;p1为管道出口压力,Pa;l 为管道长度,m;p0为管道入口压力,Pa;zg=jwLg,Rg=vb/AD,Rg为气阻,(m
18、2 s)1;Lg为气感,Lg=1/A,m1;Yg=jwCg;Cg为气容,Cg=A/RT,(kg m2)/J;Kg为受控气压源,gb2/2KvRTD,(kg m)/(J s2)。2.2 热力模型求解 针对蒸汽在管道中流动的热力模型,基于水力模型中压力和流量的计算结果,采用有限差分法进行求解。将能量方程以流量 G 和压力 p 的形式表示,根据假设,蒸汽视为理想气体,蒸汽的焓与内能均满足理想气体焓和内能公式:ddphcT (9)式中:cp为气体定压比热容,J/(kg K),依据 IAPWS-IF97 计算公式查表得到。vdducT (10)式中:cv为气体定容比热容,J/(kg K),依据 IAPW
19、S-IF97 计算公式查表得到。能量方程表示为:323332233323322vvpvvvhu RTTR G TGtA cxAcc TRGGR TTR G TxA cxAcpppppA Dp c (11)对式(11)进行向后差分,由于流量、压力和温度值变化不大,且与温度有关参数的变化也不大,所以将温度有关的时间步长和空间步长项的系数看作常数 a:nnnnp jjjjnv jnjc T R GaAp c (12)式中:在有限差分法中,设置 n 为时间项数,j 为空间项数。将其余项当作温度的函数:3231332233323132 2nnnnnnnnjjjjjjjjnjnnnnjv jjv jnnn
20、nnnnnjjjjjjjjnnnnjv jjv jnjhuR TR G TGGf TxAcAp cG R TppR G TxApcA Dppc(13)方程经过简化为:+11 0nnnnjjjjnjTTTTaf Ttx (14)当x0,t0 时,有限差分格式的截断误差为 M(x,y):222211,22TTM x txtatxtx (15)式中:为时间,s;为空间,m。124 2023 年 http:/ M(x,t)1 证明了向后差分格式与微分方程相容,满足相容性条件。假设方程差分格式的解为njT,微分方程的解为T(xj,tn)(表示微分方程在 xj位置、tn时刻的温度解),差分格式解与微分方程
21、解的误差为:+1111,1,njjjnnjjnnnnhTT x taha hM x tt (16)式中:=t/x,称为网格比。当满足 1a0 条件时,若x0,t0 时,有1jnh0,则11(,)jnnjTT x t,向后差分格式收敛。对差分格式稳定性进行分析,当时间步长比较 小时,经数据分析可得,()0njf Tt,将其扩展到整个实轴,并将两侧用傅里叶积分表示,化简后可得:111enikntTaT (17)式中 k 由傅里叶变换得到。可得到增长因子项为:22,141sin2k tG t kaa (18)当 1a0 时,|G(t,k)|211+M(x,t)t,满足 von Neumann 条件,
22、故差分格式稳定。总而言之,向后差分格式满足相容性条件,并且在 1a0 时,满足收敛性和稳定性条件,可用于微分方程的求解。3 模型验证与方法对比 3.1 研究对象 本文拟对长度为 6.0 km 的单管进行研究和分析,以蒸汽入口压力、入口温度和出口流量作为条件输入。管道参数以及变量范围见表 1,入口温度和出口流量随时间变化曲线如图 1 和图 2 所示,以入口流量、出口压力以及管道 1.5、3.0、4.5、6.0 km处的温度作为输出,设置仿真时间为 6 h。表 1 管道参数和输入变量范围 Tab.1 Scope of pipeline parameters and input variables
23、管道参数 数值 输入变量 数值 长度/km 6.0 入口压力/MPa 4.04.5 内径/mm 400500 入口温度/K 578.15588.15 沿程摩阻系数 0.03 出口流量/(kg.s1)7.98.2 3.2 模型求解及验证 3.2.1 模型求解条件 1)不同空间步长 在其他条件相同的情况下,设置管道入口压力为4.5 MPa,管道内径为400 mm,求解时间步长为 10 s,空间步长分别为 6.0、10.0、12.0、20.0、30.0 m,对热力模型求解结果分析。图 1 入口温度变化曲线 Fig.1 Inlet temperature change curve 图 2 出口流量变化
24、曲线 Fig.2 Egress flow variation curve 2)不同管道内径 在其他条件相同情况下,设置管道入口压力为 4.5 MPa,求解时间步长 10 s,空间步长为 6.0 m,管道内径分别为 420、440、460、480、500 mm,对热力模型求解结果分析。3)不同入口压力 在其他条件相同的情况下,设置管道内径为 400 mm,蒸汽管道时间步长 10 s,空间步长为6.0 m,设置蒸汽管道入口压力分别为4.0、4.1、4.2、4.3、4.4 MPa,对热力模型求解结果分析。3.2.2 模型验证 基于水力模型中压力和流量的计算结果,对温度进行求解,最终与 APROS 软
25、件仿真数据进行对比,从而得到不同条件下温度的求解误差,从而验证模型的准确性。在设置的不同离散步长、管道内径和初始压力条件下,运用仿真软件获得 1.5、3.0、4.5、6.0 km 处温度随时间的变化关系,并与有限差分法在同样条件下相同位置处温度随时间的变化关系对比,得到温度误差为:第 8 期 周 懿 等 工业园区综合能源系统蒸汽混合建模仿真研究 125 http:/ 121100%TTeT (19)式中:e为有限差分法求解误差。将误差求和汇总,并除以总数,得到某一条件下有限差分法的平均误差:0eEN (20)式中:E 为有限差分法求解平均误差;N0为参与误差计算的温度个数。3.3 软件仿真 A
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 工业园区 综合 能源 系统 蒸汽 混合 建模 仿真 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。