CN110781598B - 一种区域综合能源***运行状态计算方法 - Google Patents
一种区域综合能源***运行状态计算方法 Download PDFInfo
- Publication number
- CN110781598B CN110781598B CN201911040182.4A CN201911040182A CN110781598B CN 110781598 B CN110781598 B CN 110781598B CN 201911040182 A CN201911040182 A CN 201911040182A CN 110781598 B CN110781598 B CN 110781598B
- Authority
- CN
- China
- Prior art keywords
- representing
- model
- node
- temperature
- heat
- 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 32
- 238000004364 calculation method Methods 0.000 claims abstract description 39
- 230000008878 coupling Effects 0.000 claims abstract description 38
- 238000010168 coupling process Methods 0.000 claims abstract description 38
- 238000005859 coupling reaction Methods 0.000 claims abstract description 38
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 127
- 238000012546 transfer Methods 0.000 claims description 38
- 238000005265 energy consumption Methods 0.000 claims description 33
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 18
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 17
- 239000003546 flue gas Substances 0.000 claims description 17
- 239000002918 waste heat Substances 0.000 claims description 17
- 239000000498 cooling water Substances 0.000 claims description 15
- 239000012530 fluid Substances 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 15
- 239000003507 refrigerant Substances 0.000 claims description 12
- 239000007789 gas Substances 0.000 claims description 9
- 239000003345 natural gas Substances 0.000 claims description 9
- 238000010438 heat treatment Methods 0.000 claims description 8
- 230000006835 compression Effects 0.000 claims description 6
- 238000007906 compression Methods 0.000 claims description 6
- 239000011521 glass Substances 0.000 claims description 6
- 238000009413 insulation Methods 0.000 claims description 6
- 229910052751 metal Inorganic materials 0.000 claims description 6
- 239000002184 metal Substances 0.000 claims description 6
- 239000002689 soil Substances 0.000 claims description 6
- 238000010521 absorption reaction Methods 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000001816 cooling Methods 0.000 claims description 3
- 238000001704 evaporation Methods 0.000 claims description 3
- 230000008020 evaporation Effects 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 3
- 239000007924 injection Substances 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 238000012887 quadratic function Methods 0.000 claims description 3
- 238000009423 ventilation Methods 0.000 claims description 3
- 239000003517 fume Substances 0.000 claims 1
- 150000002739 metals Chemical class 0.000 claims 1
- 238000011217 control strategy Methods 0.000 abstract description 2
- 238000010248 power generation Methods 0.000 description 10
- 239000000243 solution Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 238000004321 preservation Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005485 electric heating Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Landscapes
- Air Conditioning Control Device (AREA)
Abstract
本发明公开了一种区域综合能源***运行状态计算方法,所述方法包括,建立以下计算模型:建筑热负荷模型、热力管道模型、耦合设备模型、电力***模型、热电联产机组模型、地源热泵机组模型;基于上述模型,依次获取负荷节点热功率、热力***的运行状态结果、热电联产机组模型运行状态结果、地源热泵机组模型的运行状态结果及节点的电功率,最终获取区域综合能源***的运行状态。本发明的计算方法可以用于分析区域综合能源***的运行状态,并可同时计算各种并网能源设备的内部运行状态,为研究***和设备特性、制定相应的运行控制策略提供数据基础。
Description
技术领域
本发明属于综合能源调控领域,特别涉及一种区域综合能源***运行状态计算方法。
背景技术
目前,区域综合能源***日益成为能源形态发展的新趋势,区域综合能源***提供的能源形式包括冷、热、电,对这些混合能量流的状态分析是了解区域综合能源***运行特性和实时运行状态、制定***运行控制策略的重要基础。但传统的状态估计方法主要面向单种能源形式。近些年来,多能流耦合***的状态估计已经成为一个研究热点,众多研究人员开展了相关工作。现有技术的不足之处主要是多数研究工作侧重考虑能流平衡,仅仅考虑了设备的外特性,利用大量的经验数据建模,导致多能流综合能源***的状态估计结果误差较大。同时,一些重要的设备参数不能通过基于简化模型的状态估计方法获得。
随着现在电力物联网概念的提出和试点工程的不断推进,区域综合能源***将作为智慧能源终端接入泛在电力物联网,并通过综合能源管理实现各种能流的优化控制、发掘能源***潜力和保证能源***安全经济运行。同时,区域综合能源***能够响应泛在电力物联网的调度指令,向上级电网提供辅助服务,并能够灵活地参与电力市场交易。由此可见,推进区域综合能源***建设意义重大,如何实现对区域综合能源***的运行状态进行分析评估进行越来越成为亟待解决的问题。
发明内容
针对上述问题,本发明提供了一种区域综合能源***运行状态计算方法,通过该方法能够高效地获取所述区域综合能源***的运行状态。
本发明的目的在于提供一种区域综合能源***运行状态计算方法,所述方法包括,
建立以下计算模型:
建筑热负荷模型、热力管道模型、耦合设备模型、电力***模型、热电联产机组模型、地源热泵机组模型;
基于所述建筑热负荷模型,获取热力***中负荷节点热功率;
基于所述热力管道模型,利用牛顿拉夫逊法计算热力***的运行状态结果;
根据所述热力***的运行状态结果,利用牛顿拉夫逊法迭代计算热电联产机组模型与地源热泵机组模型的运行状态结果;
根据热电联产机组模型运行状态结果、地源热泵机组模型的运行状态结果、热力***的运行状态结果以及负荷节点热功率,在所述耦合设备模型下,计算热力***中节点的电功率;
基于电力***模型与所述节点的电功率,获取区域综合能源***的运行状态。
进一步,所述热电联产机组模型包括压气机模型、中冷器模型、涡轮模型、涡轮增压器模型、气缸模型、余热换热器模型、蓄水箱模型及余热锅炉模型,其中,
所述压气机模型为:
p1=p0πc
其中,T0和T1分别表示压气机进口温度和出口温度,单位为℃;ηc表示压气机效率;πc表示压气机压比;k为空气绝热指数;p0和p1分别表示压气机进口压力和出口压力,单位为MPa;Pqtc表示压气机消耗的功率,单位为W;Da表示压气机空气流量,单位为kg/s;R表示气体常数;
所述中冷器模型为:
T2=T1(1-φ1)+φ1Tw
Tc0=K1Tw+K2T1
Dc0=K3Da
其中,T1表示压气机的出口温度,T2表示中冷器出口空气温度,单位为℃;Tw表示中冷器中冷却水进口温度,单位为℃;φ1为中冷器冷却效率;Tc0表示中冷器中冷却水出口温度,单位为℃;Dc0表示中冷器中冷却水流量,单位为kg/s;Da表示压气机空气流量,K1为中冷器中与冷却水进口温度相应的换热系数,K2均为中冷器中与压气机出口温度相应的换热系数,K3为中冷器的流量比例;
所述涡轮模型为:
其中,Tt和Tt1分别表示涡轮入口烟气温度和出口烟气温度,单位为℃;ηt表示涡轮效率;πt表示涡轮膨胀比;k为空气绝热指数;R表示气体常数;Pttc表示涡轮输出扭矩,单位为W;Da表示压气机空气流量;Df表示压天然气流量,单位为kg/s;
所述涡轮增压器模型为:
Pttcηtc-Pqtc-103×Pm=0
其中,ηtc表示涡轮增压器机械效率;Pm表示热电联产机组对外输出的功率,单位为kW;Pttc表示涡轮输出扭矩;Pqtc表示压气机消耗的功率;
所述气缸模型为:
DfHu×103+DacpaT2-Pm×103-(Da+Df)cpgTt-Qw×103=0
Qw-τwDfHu=0
Qw×103-Dmgcw(Tg1-Tg0)=0
其中,Hu表示天然气低位热值,单位为kJ/kg;cpa和cpg分别表示空气和烟气的定压比热容,单位为J/(kg·℃);Qw表示单位时间气缸中缸套水带走的热量,单位为kW;τw表示缸套水吸热效率;Dmg表示缸套水流量,单位为kg/s;cw表示缸套水的比热容,单位为J/(kg·℃);Tg0和Tg1分别表示缸套水进口和出口温度,单位为℃;kmg表示气缸传热系数,单位为W/(m2·℃);Amg表示缸套水换热壁面积,单位为m2;Tmg表示缸套水换热壁温度,单位为℃;Df表示压天然气流量;Da表示压气机空气流量;Pm表示热电联产机组对外输出的功率;Tt表示涡轮入口烟气温度;
所述余热换热器模型为:
Dmgcp1(Tg1-Tg0)-Dc1cp2(Th0-Tc1)=0
其中,cp1和cp2分别表示热水和冷水的定压比热容,单位为J/(kg·℃);Ke1和Ke2分别表示热水侧和冷水侧的传热系数,单位为W/(m2·℃);A1和A2分别表示壳程和管程的换热面积,单位为m2;Thp表示换热壁的平均温度,单位为℃;Dc1表示冷水箱出口水流量,单位为kg/s;Th0和Tc1分别表示冷水箱出口和进口温度,单位为℃;Tg0和Tg1分别表示缸套水进口和出口温度;Dmg表示缸套水流量;
所述蓄水箱模型为:
其中,Qch2表示蓄水箱的热损失,单位为W;λc表示蓄水箱的导热系数,单位为W/(m·℃);hc表示蓄水箱的高度,单位为m;dc1和dc2分别表示蓄水箱的内、外直径,单位为m;Tc2表示蓄水箱出口水温度,单位为℃;δc表示蓄水箱保温层厚度,单位为m;Dout表示外部供水流量,单位为kg/s;Tout为外部供水温度,单位为℃;T0表示压气机进口温度;Th0表示冷水箱出口温度;Dc1表示冷水箱出口水流量;cw表示缸套水的比热容;
所述余热锅炉模型为:
φ2(Da+Df)cpg(Tt1-Tt2)-Q1×103=0
Q1×103+Dc1cp2Tc2-Dc1cp2Tp=0
其中,Tt2表示余热锅炉出口烟气温度,单位为℃;φ2为余热锅炉热效率;Q1为烟气与金属之间的总换热量,单位为kW;Ky1和Ky2表示传热系数,单位为W/(m2·℃);Ah表示换热面积,单位为m2;Tm表示金属温度,单位为℃;Tp表示余热锅炉蒸汽出口温度,单位为℃。
进一步,所述地源热泵机组模型包括埋地换热器、蒸发器模型、冷凝器模型及压缩机模型,其中,
所述埋地换热器为:
GecwTew1-GecwTew2-Grh1+Grh4=0
其中,Tew1和Tew2分别表示蒸发器进口与出口冷冻水温度,单位为℃;Kd表示土壤的传热系数,单位为kW/(m2·℃);Fd表示土壤的换热面积,单位为m2;Ge表示蒸发器侧冷冻水的质量流量,单位为kg/s;Tw表示埋管井壁温度,单位为℃;Gr表示蒸发器内制冷剂质量流量,单位为kg/s;h1和h4分别表示蒸发器出口和进口制冷剂焓,单位为kJ/kg;Ke表示蒸发器传热系数,单位为kW/(m2·℃);Fe表示蒸发器传热面积,单位为m2;Te表示蒸发器的蒸发温度,单位为℃;
所述蒸发器模型为:
h1+6.5323×10-5Te 3-0.0011Te 2-0.3897Te-621.9697=0
vd-e-1.2013×106Te 3+1.3666×10-4Te 2-0.0313Te-3.0703=0
其中,vd表示制冷剂的比容,单位为m3/kg;pe为蒸发器中制冷剂温度为Te时的饱和压力,单位为Pa;
所述冷凝器模型为:
GccwTcw2-Grh2+Grh4-GccwTcw1=0
h4-1.3587Tc-418.2381=0
其中,Gc表示冷凝器中冷却水的质量流量,单位为kg/s;Tcw1和Tcw2分别表示冷凝器进口、出口冷冻水温度,单位为℃;Kc表示冷凝器传热系数,单位为kW/(m2·℃);Fc表示冷凝器传热面积,单位为m2;Tc表示冷凝器中的冷凝温度,单位为℃;
所述压缩机模型为:
ηe-0.0065ε3+0.0901ε2-0.384ε-0.321=0
ηv-0.0039ε2+0.058ε-0.9662=0
其中,pc表示冷凝器内的压力,单位为Pa;ε表示压缩机压缩比;ηe表示压缩机效率;ηv表示压缩机的容积效率;Vc表示压缩机的排气量,单位为m3;h2表示冷凝器进口焓;Z1和Z2分别表示压缩机进口和出口气体的压缩系数。
进一步,所述建筑热负荷模型为:
Kis-1,wTis,w+K1-os,wTos,w-(Kis-1,w+K1-os,w)T1,w=0
(h0,w+K1-os,w)Tos,w-K1-os,wT1,w-h0,wT0-qsolar=0
Kis-1,wT1,w-(hi+Kis-1,w)Tis,w+hiTi=0
Kis-1,rTis,r+K1-os,rTos,r-(Kis-1,r+K1-os,r)T1,r=0
(h0,r+K1-os,r)Tos,r-K1-os,rT1,r-h0,rT0-qsolar=0
Kis-1,rT1,r-(hi+Kis-1,r)Tis,r+hiTi=0
Kis-1,fTis,f+K1-os,fTos,f-(Kis-1,f+K1-os,f)T1,f=0
(h0,f+K1-os,f)Tos,f-K1-os,fT1,f-h0,fT0-qsolar=0
Kis-1,fT1,f-(hi+Kis-1,f)Tis,f+hiTi=0
-hiTi(Aw+Ar+Af)=0
qh-qz+Agqsolar=0
其中,Kis-1,w、Kis-1,r和Kis-1,f分别表示墙体、屋顶和地板与屋内空气的传热系数,单位为kW/(m2·℃);K1-os,w、K1-os,r和K1-os,f表示墙体、屋顶和地板与外表面传热系数,单位为kW/(m2·℃);Tis,w、Tis,r和Tis,f分别表示墙体、屋顶和地板内表面温度,单位为℃;Tos,w、Tos,r和Tos,f分别表示墙体、屋顶和地板外表面温度,单位为℃;T1,w、T1,r和T1,f分别表示墙体、屋顶和地板内部温度,单位为℃;h0,w、h0,r和h0,f分别表示墙体、屋顶和地板外表面与室外空气对流换热系数,单位为kW/(m2·℃);qsolar表示日照强度,单位为kW/m2;Ti表示室内温度;hi表示墙体、屋顶和地板内表面与室内空气对流换热系数,单位为kW/(m2·℃);qz表示建筑耗热量,单位为kW;Aw、Ar、Af和Ag分别表示墙体、屋顶、地板和玻璃面积,单位为m2;nk表示建筑换气次数,单位为次/h;Vi表示建筑体积,单位为m3;ρi表示室内空气密度,单位为kg/m3;cpi表示室内空气比热容,单位为kJ/(kg·℃);Kg表示玻璃材料传热系数,单位为kW/(m2·℃);qh表示供热***末端装置供热量,单位为kW。
进一步,所述热力管道模型包括水力模型与热力模型,其中,
所述水力模型为:
AM-Q=0
BfΔH=0
其中,A表示热网的关联矩阵;Bf表示热网的回路矩阵;M表示热网中的流量列向量;Q表示节点注入流量列向量;ΔH表示压差列向量;Z表示节点高度差列向量;HP表示水泵压头列向量;Hv表示阀门处压力损失列向量;F表示管道阻力系数列向量;L表示管道长度列向量;D表示管道直径列向量;S表示管道横截面积列向量;表示流体的平均密度;
水力模型中阀门的压差满足:
其中,pin和pout分别表示阀门入口和阀门出口的压力,单位为MPa;ξ表示阀门开度;d表示管道直径,单位为m;ρ表示工作流体的密度,单位为kg/m3;m表示工作流体的流量,单位为kg/s;
水力模型中泵扬程满足:
hp=pout.p-pin.p=a2m2+a1m+a0
其中,pin.p和pout.p分别表示泵入口和出口压力,单位为MPa;n表示泵实际转速,单位为r/min;n0表示泵额定转速,单位为r/min;ai表示泵多项式拟合系数,ai0表示额定工况下多项式拟合系数,i为整数;m表示工作流体的流量。
所述热力模型为:
Te-Ta-E(Ad TTn-Ta)=0
CpAudiag(M)Te-CpAddiag(M)Ad TTn-QJ=0
其中,Au和Ad分别表示上、下三角关联矩阵;Te表示管道末端处的温度向量;Tn表示节点温度列向量;Ta表示环境温度列向量;E表示温度衰减系数对角矩阵;M表示热网中的流量列向量;Cp表示热力管网节点处流体的比热容列向量;QJ表示节点处热负荷列向量。
进一步,所述方法还包括对所述热力模型进行求解时,指定出热力***中节点温度的值或指定出热力***中单个节点的热功率的值,具体包括:
当所述热力***节点为非源非负荷节点时:单个节点的热功率满足QJ=0;
当所述热力***节点为非源负荷节点时:QJ=Qh;
当所述热力***节点为定温节点时,所述节点温度为指定值。
进一步,所述耦合设备模型包括:热电联产机组耦合接口模型、地源热泵机组能耗模型、空调风机盘管、水泵;其中,
所述热电联产机组耦合接口模型对应电力***中的电源节点;
所述地源热泵机组能耗模型、空调风机盘管、水泵对应电力***中的电负荷节点。
进一步,所述计算热力***中节点的电功率以及基于电力***模型与所述节点的电功率,获取区域综合能源***的运行状态具体包括:
对所述热电联产机组耦合接口模型进行热力***与电力***联合计算,其中,
根据热力潮流计算,获取所述热电联产机组耦合接口模型的对外输出功率Pm;
将所述对外输出功率Pm与第一比例系数相乘,计算得出热电联产机组耦合接口在电力***中对应电源节点的电功率Pe;
对所述地源热泵机组能耗模型进行热力***与电力***联合计算,其中,
根据热力潮流计算,获取地源热泵机组能耗模型中所述压缩机的输出功率;
将所述压缩机的输出功率与第二比例系数相乘,计算得出地源热泵机组能耗模型在电力***中对应负荷节点的电功率Pl;
对所述水泵进行热力***与电力***联合计算,其中,
获取所述水泵的能耗P,其中,
其中,P表示水泵能耗,单位为kW;G表示质量流量,单位为kg/s;g表示重力加速度,数值为9.8N/kg;H表示扬程,单位为Pa;η表示水泵效率,是流量G的二次函数,记为η=k2G2+k1G+k0,k2、k1、k0为拟合系数;
基于所述水泵能耗P,计算得出水泵的电机功率为Pe(kW),其中,
对所述空调风机盘管进行热力***与电力***联合计算,其中,
根据空调风机盘管的换热量qh乘以第三比例系数,计算得出空调风机盘的电负荷量PL。
进一步,所述电力***模型为:
i,j=1,2,…,n
其中,Pi和Qi分别表示电力线路节点i处的有功功率和无功功率;Gij和Bij分别表示节点i与节点j的电导矩阵和电纳矩阵;θij表示节点i与节点j的相角差;Vi和Vj分别表示节点i和节点j处的电压。
本发明提出的计算方法能够促进区域综合能源***安全高效运行,具体表现为:基于区域综合能源***中的一个或对多个模型对运行状态计算显著提高了计算结果的准确性,也能够及时获取综合能源***重要运行参数;借助能源***运行状态计算,能够仿真不同运行策略的实施效果,为区域综合能源***优化运行策略的制定提供数据支撑。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1示出了本发明实施例中一种区域综合能源***运行状态计算方法流程示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地说明,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明实施例中一种区域综合能源***运行状态计算方法,所述方法包括,首先,建立以下计算模型:建筑热负荷模型、热力管道模型、耦合设备模型、电力***模型、热电联产机组模型、地源热泵机组模型;其次,基于所述建筑热负荷模型,获取热力***中负荷节点热功率;然后,基于所述热力管道模型,利用牛顿拉夫逊法计算热力***的运行状态结果;然后,根据所述热力***的运行状态结果,利用牛顿拉夫逊法迭代计算热电联产机组模型与地源热泵机组模型的运行状态结果;然后,根据热电联产机组模型运行状态结果、地源热泵机组模型的运行状态结果、热力***的运行状态结果以及负荷节点热功率,在所述耦合设备模型下,计算热力***中节点的电功率;最后,基于电力***模型与所述节点的电功率,获取区域综合能源***的运行状态。采用上述计算方法显著提高了区域综合能源***运行状态计算结果的准确性,也能够及时获取综合能源***重要运行参数;借助能源***运行状态计算,能够进一步仿真不同运行策略的实施效果,为区域综合能源***优化运行策略的制定提供了数据支撑。
具体的,(1)所述热电联产机组模型具体包括:
(1-1)压气机模型:
p1=p0πc
其中,T0和T1分别表示压气机进口温度和出口温度,单位为℃;ηc表示压气机效率;πc表示压气机压比;k为空气绝热指数;p0和p1分别表示压气机进口压力和出口压力,单位为MPa;Pqtc表示压气机消耗的功率,单位为W;Da表示压气机空气流量,单位为kg/s;R表示气体常数。
(1-2)中冷器模型:
T2=T1(1-φ1)+φ1Tw
Tc0=K1Tw+K2T1
Dc0=K3Da
其中,T1表示压气机的出口温度,T2表示中冷器出口空气温度,单位为℃;Tw表示中冷器中冷却水进口温度,单位为℃;φ1为中冷器冷却效率;Tc0表示中冷器中冷却水出口温度,单位为℃;Dc0表示中冷器中冷却水流量,单位为kg/s;Da表示压气机空气流量,K1为中冷器中与冷却水进口温度相应的换热系数,K2均为中冷器中与压气机出口温度相应的换热系数,K3为中冷器的流量比例。
(1-3)涡轮模型:
其中,Tt和Tt1分别表示涡轮入口烟气温度和出口烟气温度,单位为℃;ηt表示涡轮效率;πt表示涡轮膨胀比;k为空气绝热指数;R表示气体常数;Pttc表示涡轮输出扭矩,单位为W;Da表示压气机空气流量;Df表示压天然气流量,单位为kg/s。
(1-4)涡轮增压器模型:
Pttcηtc-Pqtc-103×Pm=0
其中,ηtc表示涡轮增压器机械效率;Pm表示热电联产机组对外输出的功率,单位为kW;Pttc表示涡轮输出扭矩;Pqtc表示压气机消耗的功率。
(1-5)气缸模型:
DfHu×103+DacpaT2-Pm×103-(Da+Df)cpgTt-Qw×103=0
Qw-τwDfHu=0
Qw×103-Dmgcw(Tg1-Tg0)=0
其中,Hu表示天然气低位热值,单位为kJ/kg;cpa和cpg分别表示空气和烟气的定压比热容,单位为J/(kg·℃);Qw表示单位时间气缸中缸套水带走的热量,单位为kW;τw表示缸套水吸热效率;Dmg表示缸套水流量,单位为kg/s;cw表示缸套水的比热容,单位为J/(kg·℃);Tg0和Tg1分别表示缸套水进口和出口温度,单位为℃;kmg表示气缸传热系数,单位为W/(m2·℃);Amg表示缸套水换热壁面积,单位为m2;Tmg表示缸套水换热壁温度,单位为℃;Df表示压天然气流量;Da表示压气机空气流量;Pm表示热电联产机组对外输出的功率;Tt表示涡轮入口烟气温度。
(1-6)余热换热器模型:
Dmgcp1(Tg1-Tg0)-Dc1cp2(Th0-Tc1)=0
其中,cp1和cp2分别表示热水和冷水的定压比热容,单位为J/(kg·℃);Ke1和Ke2分别表示热水侧和冷水侧的传热系数,单位为W/(m2·℃);A1和A2分别表示壳程和管程的换热面积,单位为m2;Thp表示换热壁的平均温度,单位为℃;Dc1表示冷水箱出口水流量,单位为kg/s;Th0和Tc1分别表示冷水箱出口和进口温度,单位为℃;Tg0和Tg1分别表示缸套水进口和出口温度;Dmg表示缸套水流量。
(1-7)蓄水箱模型:
其中,Qch2表示蓄水箱的热损失,单位为W;λc表示蓄水箱的导热系数,单位为W/(m·℃);hc表示蓄水箱的高度,单位为m;dc1和dc2分别表示蓄水箱的内、外直径,单位为m;Tc2表示蓄水箱出口水温度,单位为℃;δc表示蓄水箱保温层厚度,单位为m;Dout表示外部供水流量,单位为kg/s;Tout为外部供水温度,单位为℃;T0表示压气机进口温度;Th0表示冷水箱出口温度;Dc1表示冷水箱出口水流量;cw表示缸套水的比热容。
(1-8)余热锅炉模型:
φ2(Da+Df)cpg(Tt1-Tt2)-Q1×103=0
Q1×103+Dc1cp2Tc2-Dc1cp2Tp=0
其中,Tt2表示余热锅炉出口烟气温度,单位为℃;φ2为余热锅炉热效率;Q1为烟气与金属之间的总换热量,单位为kW;Ky1和Ky2表示传热系数,单位为W/(m2·℃);Ah表示换热面积,单位为m2;Tm表示金属温度,单位为℃;Tp表示余热锅炉蒸汽出口温度,单位为℃。
本实施例中,(1-1)-(1-8)所涉及的模型方程中,待定的状态参数为T1、T2、Pqtc、Dc0、Tc0、Tt、Tt1、Pttc、Pm、Df、Qw、Tg1、Tg0、Tmg、Thp、Th0、Tc2、Qch2、Q1、Tt2和Tm。
(2)所述地源热泵机组模型具体包括:
(2-1)埋地换热器:
GecwTew1-GecwTew2-Grh1+Grh4=0
其中,Tew1和Tew2分别表示蒸发器进口与出口冷冻水温度,单位为℃;Kd表示土壤的传热系数,单位为kW/(m2·℃);Fd表示土壤的换热面积,单位为m2;Ge表示蒸发器侧冷冻水的质量流量,单位为kg/s;Tw表示埋管井壁温度,单位为℃;Gr表示蒸发器内制冷剂质量流量,单位为kg/s;h1和h4分别表示蒸发器出口和进口制冷剂焓,单位为kJ/kg;Ke表示蒸发器传热系数,单位为kW/(m2·℃);Fe表示蒸发器传热面积,单位为m2;Te表示蒸发器的蒸发温度,单位为℃。
(2-2)蒸发器模型:
h1+6.5323×10-5Te 3-0.0011Te 2-0.3897Te-621.9697=0
vd-e-1.2013×106Te 3+1.3666×10-4Te 2-0.0313Te-3.0703=0
其中,vd表示制冷剂的比容,单位为m3/kg;pe为蒸发器中制冷剂温度为Te时的饱和压力,单位为Pa。
(2-3)冷凝器模型:
GccwTcw2-Grh2+Grh4-GccwTcw1=0
h4-1.3587Tc-418.2381=0
其中,Gc表示冷凝器中冷却水的质量流量,单位为kg/s;Tcw1和Tcw2分别表示冷凝器进口、出口冷冻水温度,单位为℃;Kc表示冷凝器传热系数,单位为kW/(m2·℃);Fc表示冷凝器传热面积,单位为m2;Tc表示冷凝器中的冷凝温度,单位为℃。
(2-4)压缩机模型:
ηe-0.0065ε3+0.0901ε2-0.384ε-0.321=0
ηv-0.0039ε2+0.058ε-0.9662=0
其中,pc表示冷凝器内的压力,单位为Pa;ε表示压缩机压缩比;ηe表示压缩机效率;ηv表示压缩机的容积效率;Vc表示压缩机的排气量,单位为m3;h2表示冷凝器进口焓;Z1和Z2分别表示压缩机进口和出口气体的压缩系数。
本实施例中,(2-1)-(2-4)所涉及的模型方程中,待定的状态参数为Tew1、Tew2、h1、h4、Gr、Te、νd、pe、Ge、h2、Tc、pc、ε、ηe和ηv。
本实施例中,所述电联产机组模型运行状态结果与地源热泵机组模型的运行状态结果至少包括模型中任一节点的温度和流经节点的流量。
(3)所述建筑热负荷模型:
Kis-1,wTis,w+K1-os,wTos,w-(Kis-1,w+K1-os,w)T1,w=0
(h0,w+K1-os,w)Tos,w-K1-os,wT1,w-h0,wT0-qsolar=0
Kis-1,wT1,w-(hi+Kis-1,w)Tis,w+hiTi=0
Kis-1,rTis,r+K1-os,rTos,r-(Kis-1,r+K1-os,r)T1,r=0
(h0,r+K1-os,r)Tos,r-K1-os,rT1,r-h0,rT0-qsolar=0
Kis-1,rT1,r-(hi+Kis-1,)Tis,r+hiTi=0
Kis-1,fTis,f+K1-os,fTos,f-(Kis-1,f+K1-os,f)T1,f=0
(h0,f+K1-os,f)Tos,f-K1-os,fT1,f-h0,fT0-qsolar=0
Kis-1,fT1,f-(hi+Kis-1,f)Tis,f+hiTi=0
-hiTi(Aw+Ar+Af)=0
qh-qz+Agqsolar=0
其中,Kis-1,w、Kis-1,r和Kis-1,f分别表示墙体、屋顶和地板与屋内空气的传热系数,单位为kW/(m2·℃);K1-os,w、K1-os,r和K1-os,f表示墙体、屋顶和地板与外表面传热系数,单位为kW/(m2·℃);Tis,w、Tis,r和Tis,f分别表示墙体、屋顶和地板内表面温度,单位为℃;Tos,w、Tos,r和Tos,f分别表示墙体、屋顶和地板外表面温度,单位为℃;T1,w、T1,r和T1,f分别表示墙体、屋顶和地板内部温度,单位为℃;h0,w、h0,r和h0,f分别表示墙体、屋顶和地板外表面与室外空气对流换热系数,单位为kW/(m2·℃);qsolar表示日照强度,单位为kW/m2;Ti表示室内温度;hi表示墙体、屋顶和地板内表面与室内空气对流换热系数,单位为kW/(m2·℃);qz表示建筑耗热量,单位为kW;Aw、Ar、Af和Ag分别表示墙体、屋顶、地板和玻璃面积,单位为m2;nk表示建筑换气次数,单位为次/h;Vi表示建筑体积,单位为m3;ρi表示室内空气密度,单位为kg/m3;cpi表示室内空气比热容,单位为kJ/(kg·℃);Kg表示玻璃材料传热系数,单位为kW/(m2·℃);qh表示供热***末端装置供热量(即热力***中负荷节点热功率),单位为kW。
本实施例中,建筑热负荷模型中的待定的状态参数为T1,w、Tos,w、Tis,w、T1,r、Tos,r、Tis,r、T1,f、Tos,f、Tis,f、qz和qh。
需要说明的是:所述建筑热负荷模型中的i起标记作用,表示内侧或室内,is表示内表面,os表示外表面。
(4)所述热力管道模型具体包括:
(4-1)水力模型:
AM-Q=0
BfΔH=0
其中,A表示热网的关联矩阵;Bf表示热网的回路矩阵;M表示热网中的流量列向量;Q表示节点注入流量列向量;ΔH表示压差列向量;Z表示节点高度差列向量;HP表示水泵压头列向量;Hv表示阀门处压力损失列向量;F表示管道阻力系数列向量;L表示管道长度列向量;D表示管道直径列向量;S表示管道横截面积列向量;表示流体的平均密度。
进一步,所述水力模型中阀门压差公式为:
其中,pin和pout分别表示阀门入口和阀门出口的压力,单位为MPa;ξ表示阀门开度;d表示管道直径,单位为m;ρ表示工作流体的密度,单位为kg/m3;m表示工作流体的流量,单位为kg/s;
进一步,所述水力模型中泵扬程公式为:
hp=pout.p-pin.p=a2m2+a1m+a0
其中,pin.p和pout.p分别表示泵入口和出口压力,单位为MPa;n表示泵实际转速,单位为r/min;n0表示泵额定转速,单位为r/min;ai表示泵多项式拟合系数,ai0表示额定工况下多项式拟合系数,i为整数;m表示工作流体的流量。
(4-2)热力模型为:
Te-Ta-E(Ad TTn-Ta)=0
CpAudiag(M)Te-CpAddiag(M)Ad TTn-QJ=0
其中,Au和Ad分别表示上、下三角关联矩阵,即热力管网节点-支路的关联矩阵;Te表示管道末端处的温度向量;Tn表示节点温度列向量;Ta表示环境温度列向量;E表示温度衰减系数对角矩阵;M表示热网中的流量列向量;Cp表示热力***节点处流体的比热容列向量;QJ表示节点处热负荷列向量;
本实施例中,在求解所述热力模型时,上述热力模型方程个数是B+N个(B为管道数,N为节点数),而未知状态量Te、Tn、Qj是B+2N个,从而对于热力管网中节点温度和热功率需要指定一个值才能够对上述热力模型进行求解,其中在指定时,充分考虑了热力管网节点的属性,具体的,
当所述热力管网节点为非源非负荷节点时:单个节点的热功率满足QJ=0;
当所述热力管网节点为非源负荷节点时:QJ=Qh,即节点处的热负荷值;
当所述热力管网节点为源节点时:通过地源热泵机组模型计算得到冷凝器出口冷冻水温度Tcw2,所述Tcw2此时为单个节点的温度。更具体地,本发明实施例中,热力管网为热力管道,热力管网在计算Tcw2时,首先根据热力管道模型中的水力模型计算得到的地源热泵机组模型中的Gr和指定的回水温度Tcw1(也作冷凝器进口冷冻水温度),再进一步通过地源热泵机组模型计算得到Tcw2(供热管道首端温度),作为定温节点参与计算。进一步,需要说明的是,热力模型针对的是热力管道(管网)的精细模型,从而在整个***中,所述冷凝器出口与供热管道首端对应连接,从而冷凝器出口温度等于供热管道首端温度。
本实施例中,所述热力***的运行状态结果包括热力***中各节点的热功率/温度以及各管道的流量、首末端温度和压差等。
(5)所述耦合设备模型具体的包括热电联产机组耦合接口模型、地源热泵机组能耗模型、空调风机盘管、水泵。
在耦合热备模型中,所述热电联产机组耦合接口模型对应电力***中的电源节点;所述地源热泵机组能耗模型、空调风机盘管、水泵对应电力***中的电负荷节点。本实施例中,通过建立耦合设备模型将热力***与电力***联合计算,且热电耦合计算一般遵循以热定电的原则,即先根据热网潮流计算耦合设备的电功率,再对应到电力***中计算出对应节点的有功功率值,之后完成电力***潮流计算。具体的:
(5-1)热电联产机组耦合接口模型:
通过热力潮流计算可以得到热电联产机组耦合接口模型的对外输出功率Pm,即原动机机械功率,乘以第一比例系数得到电力***中对应电源点的电功率Pe值。其中,所述第一比例系数为热电联产机组耦合接口设备在热力***中的输出功率与在电力***中的电功率的比值,该比值为已知值。进一步需要说明的是:热电联产机组耦合接口模型与热电联产机组模型都是以热电联产机组为基础,其中,热电联产机组模型为热电联机组的热力细化模型,电联产机组耦合接口模型为热电联机组在耦合设备中耦合接口部分的模型,从而电联产机组耦合接口模型与热电联产机组模型的输出功率一致。
(5-2)地源热泵机组能耗模型
地源热泵机组能耗模型的耗能主要由热泵机组(即压缩机能耗)、循环泵***、空调、末端设备(风机盘管)、换热器传热能耗构成。其中,热泵机组能耗占40%,空调、循环泵能耗20%,换热器传热能耗20%。
本实施例中,以压缩机的能耗作为地源热泵机组能耗模型的能耗,具体的,根据热力潮流计算得到压缩机的输出功率,即压缩机能耗,乘以第二比例系数得到电力***对应负荷点的电功率Pl值。其中,所述第二比例系数为地源热泵机组能耗模型在热力***中的输出功率与在电力***中的电功率的比值,该比值为已知值。
(5-3)水泵
上式中,P表示水泵能耗,单位为kW;G表示质量流量,单位为kg/s;g表示重力加速度,数值为9.8N/kg;H表示扬程,单位为Pa;η表示水泵效率,是流量G的二次函数,记为η=k2G2+k1G+k0,k2、k1、k0为拟合系数,通过实验手段得到。水泵电机功率为Pe(kW),有:
(5-4)空调风机盘管
空调、风机盘管装设于用户处,对于变频空调,可以根据换热量qh乘以第三比例系数计算电负荷量(电动率)PL。本实施例中,换热器传热损耗不计入电损耗中。所述第三比例系数为空调风机盘管在热力***中的换热量与在电力***中电负荷量的比值,该比值为已知值。
(6)电力***模型:
i,j=1,2,…,n
上式中,Pi和Qi分别表示节点i处的有功功率和无功功率;Gij和Bij分别表示节点i与节点j的电导矩阵和电纳矩阵;θij表示节点i与节点j的相角差;Vi和Vj分别表示节点i和节点j处的电压。
(7)将以上所有模型方程联立,基于牛顿拉夫逊法迭代进行求解,即可求取区域综合能源***的运行状态。
优选地,主要是针对热力管道模型、电力***模型、耦合设备模型利用牛顿拉夫逊法进行联合迭代计算,计算结果中的一个或多个结果作为热电联产机组、地源热泵机组模型计算迭代的输入值,热电联产机组、地源热泵小模型再次利用牛顿拉夫逊法单独迭代计算,最终得到结果。建筑热负荷模型不需要迭代,是由输入条件根据热负荷模型计算出热负荷功率作为电热联合计算的输入参数之一。进一步具体的:1、根据已知热负荷参数利用热负荷模型计算热负荷节点的热负荷功率值;2、根据热力管道模型,利用牛顿拉夫逊法得到各节点的热功率/温度以及各管道的流量、首末端温度和压差;3、根据热力管道计算结果,将对应热电联产机组模型、地源热泵机组模型节点的相关结果(节点的温度和流经节点的流量)作为输入参数与已知的设备参数一起,利用牛顿拉夫逊法计算热电联产机组、地源热泵机组模型的运行状态;4、根据热网的计算结果,在耦合设备处根据热电耦合关系计算该节点电功率(热电联产机组模型是发电出力,其他是电负荷);5、根据计算出的节点电功率,利用电力***模型,根据牛顿拉夫逊法或其他电力***潮流算法计算电力***潮流。
尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (8)
1.一种区域综合能源***运行状态计算方法,其特征在于,所述方法包括,
建立以下计算模型:
建筑热负荷模型、热力管道模型、耦合设备模型、电力***模型、热电联产机组模型、地源热泵机组模型;
其中,所述耦合设备模型包括:热电联产机组耦合接口模型、地源热泵机组能耗模型、空调风机盘管、水泵;
基于所述建筑热负荷模型,获取热力***中负荷节点热功率;
基于所述热力管道模型,利用牛顿拉夫逊法计算热力***的运行状态结果;
根据所述热力***的运行状态结果,利用牛顿拉夫逊法迭代计算热电联产机组模型与地源热泵机组模型的运行状态结果;
根据热电联产机组模型运行状态结果、地源热泵机组模型的运行状态结果、热力***的运行状态结果以及负荷节点热功率,在所述耦合设备模型下,计算热力***中节点的电功率;
基于电力***模型与所述节点的电功率,获取区域综合能源***的运行状态;
其中,所述计算热力***中节点的电功率以及基于电力***模型与所述节点的电功率,获取区域综合能源***的运行状态具体包括:
对所述热电联产机组耦合接口模型进行热力***与电力***联合计算;
其中,根据热力潮流计算,获取所述热电联产机组耦合接口模型的对外输出功率Pm;
将所述对外输出功率Pm与第一比例系数相乘,计算得出热电联产机组耦合接口在电力***中对应电源节点的电功率Pe;
对所述地源热泵机组能耗模型进行热力***与电力***联合计算;
其中,根据热力潮流计算,获取地源热泵机组能耗模型中压缩机的输出功率;
将所述压缩机的输出功率与第二比例系数相乘,计算得出地源热泵机组能耗模型在电力***中对应负荷节点的电功率Pl;
对所述水泵进行热力***与电力***联合计算;
其中,获取所述水泵的能耗P,
P表示水泵能耗;G表示质量流量;g表示重力加速度;H表示扬程;η表示水泵效率,是流量G的二次函数,记为η=k2G2+k1G+k0,k2、k1、k0为拟合系数;
基于所述水泵能耗P,计算得出水泵的电机功率为Pe(kW),
对所述空调风机盘管进行热力***与电力***联合计算;
其中,根据空调风机盘管的换热量qh乘以第三比例系数,计算得出空调风机盘的电负荷量PL。
2.根据权利要求1所述的计算方法,其特征在于,所述热电联产机组模型包括压气机模型、中冷器模型、涡轮模型、涡轮增压器模型、气缸模型、余热换热器模型、蓄水箱模型及余热锅炉模型,其中,
所述压气机模型为:
p1=p0πc
其中,T0和T1分别表示压气机进口温度和出口温度;ηc表示压气机效率;πc表示压气机压比;k为空气绝热指数;p0和p1分别表示压气机进口压力和出口压力;Pqtc表示压气机消耗的功率;Da表示压气机空气流量;R表示气体常数;
所述中冷器模型为:
T2=T1(1-φ1)+φ1Tw
Tc0=K1Tw+K2T1
Dc0=K3Da
其中,T1表示压气机的出口温度,T2表示中冷器出口空气温度;Tw表示中冷器中冷却水进口温度;φ1为中冷器冷却效率;Tc0表示中冷器中冷却水出口温度,单位为℃;Dc0表示中冷器中冷却水流量;Da表示压气机空气流量,K1为中冷器中与冷却水进口温度相应的换热系数,K2均为中冷器中与压气机出口温度相应的换热系数,K3为中冷器的流量比例;
所述涡轮模型为:
其中,Tt和Tt1分别表示涡轮入口烟气温度和出口烟气温度;ηt表示涡轮效率;πt表示涡轮膨胀比;k为空气绝热指数;R表示气体常数;Pttc表示涡轮输出扭矩;Da表示压气机空气流量;Df表示压天然气流量;
所述涡轮增压器模型为:
Pttcηtc-Pqtc-103×Pm=0
其中,ηtc表示涡轮增压器机械效率;Pm表示热电联产机组对外输出的功率;Pttc表示涡轮输出扭矩;Pqtc表示压气机消耗的功率;
所述气缸模型为:
DfHu×103+DacpaT2-Pm×103-(Da+Df)cpgTt-Qw×103=0
Qw-τwDfHu=0
Qw×103-Dmgcw(Tg1-Tg0)=0
其中,Hu表示天然气低位热值;cpa和cpg分别表示空气和烟气的定压比热容;Qw表示单位时间气缸中缸套水带走的热量;τw表示缸套水吸热效率;Dmg表示缸套水流量;cw表示缸套水的比热容;Tg0和Tg1分别表示缸套水进口和出口温度;kmg表示气缸传热系数;Amg表示缸套水换热壁面积;Tmg表示缸套水换热壁温度;Df表示压天然气流量;Da表示压气机空气流量;Pm表示热电联产机组对外输出的功率;Tt表示涡轮入口烟气温度;
所述余热换热器模型为:
Dmgcp1(Tg1-Tg0)-Dc1cp2(Th0-Tc1)=0
其中,cp1和cp2分别表示热水和冷水的定压比热容;Ke1和Ke2分别表示热水侧和冷水侧的传热系数;A1和A2分别表示壳程和管程的换热面积;Thp表示换热壁的平均温度;Dc1表示冷水箱出口水流量;Th0和Tc1分别表示冷水箱出口和进口温度;Tg0和Tg1分别表示缸套水进口和出口温度;Dmg表示缸套水流量;
所述蓄水箱模型为:
其中,Qch2表示蓄水箱的热损失;λc表示蓄水箱的导热系数;hc表示蓄水箱的高度;dc1和dc2分别表示蓄水箱的内、外直径;Tc2表示蓄水箱出口水温度;δc表示蓄水箱保温层厚度;Dout表示外部供水流量;Tout为外部供水温度;T0表示压气机进口温度;Th0表示冷水箱出口温度;Dc1表示冷水箱出口水流量;cw表示缸套水的比热容;
所述余热锅炉模型为:
φ2(Da+Df)cpg(Tt1-Tt2)-Q1×103=0
Q1×103+Dc1cp2Tc2-Dc1cp2Tp=0
其中,Tt2表示余热锅炉出口烟气温度;φ2为余热锅炉热效率;Q1为烟气与金属之间的总换热量;Ky1和Ky2表示传热系数;Ah表示换热面积;Tm表示金属温度;Tp表示余热锅炉蒸汽出口温度。
3.根据权利要求1所述的计算方法,其特征在于,所述地源热泵机组模型包括埋地换热器、蒸发器模型、冷凝器模型及压缩机模型,其中,
所述埋地换热器为:
GecwTew1-GecwTew2-Grh1+Grh4=0
其中,Tew1和Tew2分别表示蒸发器进口与出口冷冻水温度;Kd表示土壤的传热系数;Fd表示土壤的换热面积;Ge表示蒸发器侧冷冻水的质量流量;Tw表示埋管井壁温度;Gr表示蒸发器内制冷剂质量流量;h1和h4分别表示蒸发器出口和进口制冷剂焓;Ke表示蒸发器传热系数;Fe表示蒸发器传热面积;Te表示蒸发器的蒸发温度;
所述蒸发器模型为:
h1+6.5323×10-5Te 3-0.0011Te 2-0.3897Te-621.9697=0
其中,vd表示制冷剂的比容;pe为蒸发器中制冷剂温度为Te时的饱和压力;
所述冷凝器模型为:
GccwTcw2-Grh2+Grh4-GccwTcw1=0
h4-1.3587Tc-418.2381=0
其中,Gc表示冷凝器中冷却水的质量流量;Tcw1和Tcw2分别表示冷凝器进口、出口冷冻水温度;Kc表示冷凝器传热系数;Fc表示冷凝器传热面积;Tc表示冷凝器中的冷凝温度;
所述压缩机模型为:
ηe-0.0065ε3+0.0901ε2-0.384ε-0.321=0
ηv-0.0039ε2+0.058ε-0.9662=0
其中,pc表示冷凝器内的压力;ε表示压缩机压缩比;ηe表示压缩机效率;ηv表示压缩机的容积效率;Vc表示压缩机的排气量;h2表示冷凝器进口焓;Z1和Z2分别表示压缩机进口和出口气体的压缩系数。
4.根据权利要求1所述的计算方法,其特征在于,所述建筑热负荷模型为:
Kis-1,wTis,w+K1-os,wTos,w-(Kis-1,w+K1-os,w)T1,w=0
(h0,w+K1-os,w)Tos,w-K1-os,wT1,w-h0,wT0-qsolar=0
Kis-1,wT1,w-(hi+Kis-1,w)Tis,w+hiTi=0
Kis-1,rTis,r+K1-os,rTos,r-(Kis-1,r+K1-os,r)T1,r=0
(h0,r+K1-os,r)Tos,r-K1-os,rT1,r-h0,rT0-qsolar=0
Kis-1,rT1,r-(hi+Kis-1,r)Tis,r+hiTi=0
Kis-1,fTis,f+K1-os,fTos,f-(Kis-1,f+K1-os,f)T1,f=0
(h0,f+K1-os,f)Tos,f-K1-os,fT1,f-h0,fT0-qsolar=0
Kis-1,fT1,f-(hi+Kis-1,f)Tis,f+hiTi=0
qh-qz+Agqsolar=0
其中,Kis-1,w、Kis-1,r和Kis-1,f分别表示墙体、屋顶和地板与屋内空气的传热系数;K1-os,w、K1-os,r和K1-os,f表示墙体、屋顶和地板与外表面传热系数;Tis,w、Tis,r和Tis,f分别表示墙体、屋顶和地板内表面温度;Tos,w、Tos,r和Tos,f分别表示墙体、屋顶和地板外表面温度;T1,w、T1,r和T1,f分别表示墙体、屋顶和地板内部温度;h0,w、h0,r和h0,f分别表示墙体、屋顶和地板外表面与室外空气对流换热系数;qsolar表示日照强度;Ti表示室内温度;hi表示墙体、屋顶和地板内表面与室内空气对流换热系数;qz表示建筑耗热量;Aw、Ar、Af和Ag分别表示墙体、屋顶、地板和玻璃面积;nk表示建筑换气次数;Vi表示建筑体积;ρi表示室内空气密度;cpi表示室内空气比热容;Kg表示玻璃材料传热系数;qh表示供热***末端装置供热量。
5.根据权利要求1所述的计算方法,其特征在于,所述热力管道模型包括水力模型与热力模型,其中,
所述水力模型为:
AM-Q=0
BfΔH=0
其中,A表示热网的关联矩阵;Bf表示热网的回路矩阵;M表示热网中的流量列向量;Q表示节点注入流量列向量;ΔH表示压差列向量;Z表示节点高度差列向量;HP表示水泵压头列向量;Hv表示阀门处压力损失列向量;F表示管道阻力系数列向量;L表示管道长度列向量;D表示管道直径列向量;S表示管道横截面积列向量;表示流体的平均密度;
水力模型中阀门的压差满足:
其中,pin和pout分别表示阀门入口和阀门出口的压力;ξ表示阀门开度;d表示管道直径;ρ表示工作流体的密度;m表示工作流体的流量;
水力模型中泵扬程满足:
hp=pout.p-pin.p=a2m2+a1m+a0
其中,pin.p和pout.p分别表示泵入口和出口压力;n表示泵实际转速;n0表示泵额定转速;ai表示泵多项式拟合系数,ai0表示额定工况下多项式拟合系数,i为整数;m表示工作流体的流量;
所述热力模型为:
Te-Ta-E(Ad TTn-Ta)=0
CpAudiag(M)Te-CpAddiag(M)Ad TTn-QJ=0
其中,Au和Ad分别表示上、下三角关联矩阵;Te表示管道末端处的温度向量;Tn表示节点温度列向量;Ta表示环境温度列向量;E表示温度衰减系数对角矩阵;M表示热网中的流量列向量;Cp表示热力管网节点处流体的比热容列向量;QJ表示节点处热负荷列向量。
6.根据权利要求3-5任一所述的计算方法,其特征在于,所述方法还包括对热力模型进行求解时,指定出热力***中节点温度的值或指定出热力***中单个节点的热功率的值,具体包括:
当所述热力***节点为非源非负荷节点时:单个节点的热功率满足QJ=0;
当所述热力***节点为非源负荷节点时:QJ=Qh;
当所述热力***节点为定温节点时,所述节点温度为指定值。
7.根据权利要求1所述的计算方法,其特征在于,
所述热电联产机组耦合接口模型对应电力***中的电源节点;
所述地源热泵机组能耗模型、空调风机盘管、水泵对应电力***中的电负荷节点。
8.根据权利要求1所述的计算方法,其特征在于,所述电力***模型为:
i,j=1,2,...,n
其中,Pi和Qi分别表示电力线路节点i处的有功功率和无功功率;Gij和Bij分别表示节点i与节点j的电导矩阵和电纳矩阵;θij表示节点i与节点j的相角差;Vi和Vj分别表示节点i和节点j处的电压。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911040182.4A CN110781598B (zh) | 2019-10-29 | 2019-10-29 | 一种区域综合能源***运行状态计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911040182.4A CN110781598B (zh) | 2019-10-29 | 2019-10-29 | 一种区域综合能源***运行状态计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110781598A CN110781598A (zh) | 2020-02-11 |
CN110781598B true CN110781598B (zh) | 2024-02-09 |
Family
ID=69387486
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911040182.4A Active CN110781598B (zh) | 2019-10-29 | 2019-10-29 | 一种区域综合能源***运行状态计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110781598B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112303705A (zh) * | 2020-09-24 | 2021-02-02 | 株洲麦格米特电气有限责任公司 | 采暖热泵运行控制方法、装置、控制器及计算机可读存储介质 |
CN113156835B (zh) * | 2021-04-26 | 2022-12-06 | 杭州鸿晟电力设计咨询有限公司 | 一种用于电-热综合能源***运行控制的建模仿真方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108960503A (zh) * | 2018-07-02 | 2018-12-07 | 河海大学 | 基于内点法的综合能源***多场景优化分析方法 |
CN109255194A (zh) * | 2018-09-27 | 2019-01-22 | 河海大学 | 一种含p2h的电-热互联综合能源***光伏消纳分析方法 |
CN109347107A (zh) * | 2018-09-29 | 2019-02-15 | 河海大学 | 一种并网运行电热互联综合能源***潮流计算方法 |
CN109830957A (zh) * | 2019-02-22 | 2019-05-31 | 南方电网科学研究院有限责任公司 | 一种面向区域综合能源***的节点运算负荷计算方法 |
CN109978625A (zh) * | 2019-03-28 | 2019-07-05 | 河海大学 | 一种计及电热气网络的综合能源***多目标运行优化方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105226643B (zh) * | 2015-09-23 | 2017-11-28 | 清华大学 | 安全约束下电力***运行模拟模型快速生成与求解方法 |
-
2019
- 2019-10-29 CN CN201911040182.4A patent/CN110781598B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108960503A (zh) * | 2018-07-02 | 2018-12-07 | 河海大学 | 基于内点法的综合能源***多场景优化分析方法 |
CN109255194A (zh) * | 2018-09-27 | 2019-01-22 | 河海大学 | 一种含p2h的电-热互联综合能源***光伏消纳分析方法 |
CN109347107A (zh) * | 2018-09-29 | 2019-02-15 | 河海大学 | 一种并网运行电热互联综合能源***潮流计算方法 |
CN109830957A (zh) * | 2019-02-22 | 2019-05-31 | 南方电网科学研究院有限责任公司 | 一种面向区域综合能源***的节点运算负荷计算方法 |
CN109978625A (zh) * | 2019-03-28 | 2019-07-05 | 河海大学 | 一种计及电热气网络的综合能源***多目标运行优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110781598A (zh) | 2020-02-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | Performance and feasibility study of hybrid ground source heat pump system assisted with cooling tower for one office building based on one Shanghai case | |
CN103020481B (zh) | 一种基于节能的确定空气源热泵地暖最佳运行工况的方法 | |
CN106567748B (zh) | 非绝热气体膨胀的压缩空气储能*** | |
CN107514667A (zh) | 采用电动热泵实现热电厂跨季节蓄热放热的集中供热*** | |
CN107741070A (zh) | 一种空气源热泵热水器高密度储热一体机 | |
CN110781598B (zh) | 一种区域综合能源***运行状态计算方法 | |
CN111829059B (zh) | 一种供热***动态建模方法、模型及调控*** | |
Cai et al. | Theoretical and experimental study of combined heat and power (CHP) system integrated with ground source heat pump (GSHP) | |
CN106357160A (zh) | 一种利用热管温差发电的***及方法 | |
CN110502791A (zh) | 基于能源集线器的综合能源***稳态建模方法 | |
CN104456699B (zh) | 一种面向被动房的空气式太阳能热水新风*** | |
Kong et al. | Influence of refrigerant charge and condenser area on direct-expansion solar-assisted heat pump system for radiant floor heating | |
Xu et al. | Performance evaluation of a pipe-embedded phase change material (PE-PCM) roof integrated with solar collector | |
CN106096248A (zh) | 不同环境下的最佳供暖方式的甄选方法 | |
CN206517335U (zh) | 一种利用太阳能和温差发电的*** | |
Wang et al. | Heat transfer and thermal characteristics analysis of direct air-cooled combined heat and power plants under off-design conditions | |
CN106482266A (zh) | 一种高效节能减排的地源恒温*** | |
Al-Dabbas | Achievement of geothermal energy using ground heat exchanger in Ma’en | |
CN208205498U (zh) | 一种烟气余热驱动的双级吸收式土壤源热泵*** | |
CN103266998B (zh) | 一种循环热管式高层建筑底层温差通风发电*** | |
Tan et al. | Modeling and simulation of layered water tank based on MATLAB/SIMULINK | |
CN209214149U (zh) | 一种无冷却塔的楼宇式分布式能源*** | |
WO2007030995A1 (fr) | Conditionneur d'air a echange geothermique | |
CN106679016A (zh) | 一种高效节能减排的地源恒温装置 | |
CN206517319U (zh) | 一种利用热管温差发电的*** |
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 |