CN106503388B - 空调系统特征识别方法 - Google Patents
空调系统特征识别方法 Download PDFInfo
- Publication number
- CN106503388B CN106503388B CN201610984289.4A CN201610984289A CN106503388B CN 106503388 B CN106503388 B CN 106503388B CN 201610984289 A CN201610984289 A CN 201610984289A CN 106503388 B CN106503388 B CN 106503388B
- Authority
- CN
- China
- Prior art keywords
- model
- water
- temperature
- air
- parameter
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Abstract
本发明提供一种空调系统特征识别方法,根据既有中央空调系统的实测运行数据,采用最小二乘法对空调系统中的冷水机组性能预测模型、表冷器性能预测模型、冷却塔性能预测模型、水泵模型以及流体输配管路阻力模型的模型参数进行特征识别,得到表征空调系统各部件结构特性的特征参数。该空调系统特征识别方法将理论分析与实际测试相结合,具有精度较高,适用性较广且所需实测参数在实际空调系统中能够方便获得的特点,可作为空调系统整体性能模拟预测的基础。
Description
技术领域
本发明属于中央空调系统仿真模拟领域,涉及为一种空调系统特征识别方法。
背景技术
常规的空调系统设备设计选型是以建筑的最大负荷为依据,而实测数据表明,一般中央空调系统80%以上的运行时间在60%以下的部分负荷下运行,因此需要通过运行调节以实现系统在部分负荷下的高效运行。对空调系统在不同负荷下进行运行参数优化调节及性能评估,借助仿真模拟手段是较易实现的方法。在对既有建筑中央空调系统进行建模仿真时,由于实际中央空调系统各部件的具体结构参数难以获取,导致对依据设备具体结构参数进行的精确模拟较难实现,使得常规的建模方法在实际应用中存在局限性。
针对实际建筑中央空调系统各部件具体结构参数缺乏的问题,国内外许多学者主要是通过海量的实际数据采集,并通过拟合的方法来建立仿真模型,其存在所需实测数据巨大,且数据需覆盖所有运行工况而引起实测所需时间长的不足,同时还存在精度较低,难以在不同空调系统中推广,适用性差等弊端。而在实际既有空调系统中,现实条件允许的测量实际运行参数的手段也存在限制。因此需要提出结构参数缺乏条件下,精度较高、适用性较广且所需实测参数较少并在实际空调系统中能够方便获得的建模方法。
发明内容
技术问题:本发明提供一种在既有中央空调系统各部件具体结构参数缺乏的条件下,具有精度较高,适用性较广且所需实测参数在实际空调系统中能够方便获得的空调系统特征识别方法。
技术方案:本发明的空调系统特征识别方法,包括以下步骤:
根据既有中央空调系统的实测运行数据,采用最小二乘法分别求解得到以下几种模型的模型参数:中央空调系统中冷水机组性能预测模型、表冷器性能预测模型、冷却塔性能预测模型、水泵模型、流体输配管路阻力模型,将所得到的模型参数作为表征系统各部件结构特性的特征参数;
所述冷水机组性能预测模型包括蒸发器模型、冷凝器模型、压缩机模型和节流阀模 型:
a.所述蒸发器模型在变水量工况下为:
Qe=mw,ecp,w(twi,e-two,e)=mr(heo-hei) (3)
所述蒸发器模型在变水温工况下为:
Qe=mw,ecp,w(twi,e-two,e)=mr(heo-hei) (8)
式中,Qe1、Δte1分别为两相区的换热量和换热温差;Qe2、Δte2分别为过热区的换热量和换热温差;mw,e为冷冻水流量;cp,w为水的比热;mr为制冷剂流量;A* 1,e为与蒸发器结构和污垢热阻有关的模型参数;x1为根据不同结构形式蒸发器的水侧换热系数经验公式得到的常系数;B* 1,e为与蒸发器结构和污垢热阻有关的模型参数;y1为根据不同结构形式蒸发器的制冷剂侧换热系数经验公式得到的常系数;C* 1,e为与蒸发器结构和污垢热阻有关的模型参数;A* 2,e为与蒸发器结构和污垢热阻有关的模型参数;B* 2,e为与蒸发器结构和污垢热阻有关的模型参数;C* 2,e为与蒸发器结构和污垢热阻有关的模型参数;y2为根据不同结构形式蒸发器的制冷剂侧换热系数经验公式得到的常系数; Qe为蒸发器总换热量;twi,e,two,e分别为冷冻水进口温度和出口温度;tw1,e为两相区冷冻水入口温度;hei,heo分别为蒸发器入口焓值和出口的焓值;te为蒸发温度;teo为压缩机吸气温度,即蒸发器出口的制冷剂温度;
b.所述冷凝器模型在变水量工况下为:
Qc=mw,ccp,w(two,c-twi,c)=mr(hci-hco) (11)
所述冷凝器模型在变水温工况下为:
Qc=mw,ccp,w(two,c-twi,c)=mr(hci-hco) (16)
式中,Qc为冷凝器换热量;mw,c为冷却水流量;twi,c,two,c分别为冷凝器的冷却水进口温度和出口温度;hci,hco分别为冷凝器入口焓值和出口焓值;Δtc1,Δtc2,Δtc3分别为冷凝器过冷区、两相区和过热区的换热温差;A* 1,c为需要确定的冷凝器模型参数;x2为根据不同结构形式冷凝器的水侧换热系数经验公式得到的常系数;B* 1,c为需要确定的 冷凝器模型参数;y3为根据不同结构形式冷凝器的制冷剂侧换热系数经验公式得到的常系数;C* 1,c为需要确定的冷凝器模型参数;A* 2,c为需要确定的冷凝器模型参数;B* 2,c为需要确定的冷凝器模型参数;y4为根据不同结构形式冷凝器的制冷剂侧换热系数经验公式得到的常系数;C* 2,c为需要确定的冷凝器模型参数;A* 3,c为需要确定的冷凝器模型参数;B* 3,c为需要确定的冷凝器模型参数;C* 3,c为需要确定的冷凝器模型参数;y5为根据不同结构形式冷凝器的制冷剂侧换热系数经验公式得到的常系数;tc为冷凝温度;tco为制冷剂在冷凝器出口的温度;tw1,c,tw2,c分别为两相区冷却水进口温度和出口温度;tci为压缩机排气温度;
c.所述压缩机模型为:
Vth=ψVth0 (22)
式中,mr为制冷剂的质量流量;λ为输气系数;Vth为压缩机理论输气量;v1为压缩机吸气比容;ψ为冷水机组负荷率;Vth0为压缩机额定工况下理论输气量;Teo为吸气温度;Tci为排气温度;pc为冷凝压力;pe为蒸发压力;k为压缩过程多变指数;Pth为压缩机理论功率;Pin为压缩机实际功率;ηe为压缩机的电能效率;
d.节流阀模型为:
hco=hei (26)
teo=te+△te (27)
tco=tc-△tc (28)
式中,hco为节流阀入口焓值;hei为节流阀出口焓值;Δte为过热度;Δtc为过冷度;
所述表冷器性能预测模型包括表冷器换热量模型、表冷器传热效能模型、表冷器接触系数模型、表冷器传热单元数模型、热容比模型、出风参数模型;
所述表冷器换热量模型为:
Qb=ma,b(hai,b-hao,b)=mw,bcp,w(two,b-twi,b) (29)
所述表冷器传热效能模型为:
所述表冷器接触系数模型为:
所述表冷器传热单元数模型为:
所述热容比模型为:
所述出风参数模型,在干工况时为:
tgo,b=tgi,b-ε1,b(tgi,b-twi,b) (34)
所述出风参数模型,在湿工况时为:
tgo,b=tgi,b-ε2,b(tgi,b-tb) (36)
tso,b=tgo,b-(1-ε2,b)(tgi,b-tsi,b) (37)
所述出风参数模型,在临界工况时,采用上述干工况时、湿工况时的出风参数模型均可;
式中,Qb为表冷器换热量;ma,b为空气质量流量;hai,b为空气入口焓值;hao,b为空气出口焓值;mw,b为水流量;twi,b,two,b分别为表冷器进口水温和出口水温;ε1,b为表冷器传热效能;tgi,tgo分别为空气进口干球温度和出口干球温度;γ为热容比;NTU为表冷器传热单元数;ε2,b为表冷器接触系数;tso,b,tsi,b分别为空气出口湿球温度和进口湿球温度;ao,b为空气侧换热系数;Fb为表冷器总换热面积;cp,a为空气比热;Kb为表冷器总传热系数;tLi为空气进口露点温度;
所述冷却塔性能预测模型,在冷却塔内空气在非饱和状态下时为:
所述冷却塔性能预测模型,在冷却塔内空气在饱和状态下时为:
式中,mw,t为冷却塔水质量流量;dz为冷却塔填料竖直方向微元长度;βt为冷却塔传质系数;a为填料比表面积;Fz为冷却塔填料横截面积;Xs,w为水温对应的饱和空气含湿量;X为冷却塔内空气水蒸气质量组分;Tw,t为水温;Le为刘易斯数;cp,v为水蒸气定压比热;r0为水的汽化潜热;Ta,t为冷却塔内空气温度;ma,t为冷却塔空气质量流量;Xs,a为饱和空气含湿量;(βtaFZ)j为换热过程中的等效值;Me为模型参数;mwi,t表示入口水流量;Ht表示填料高度;
所述水泵模型为:
Hp0=a0+a1V0+a2V0 2 (50)
Pin,p0=b0+b1V0+b2V0 2 (51)
式中,V0、Hp0、Pin,p0分别为泵在转速n0下的流量、扬程与功率;V1、Hp1、Pin,p1分别为泵在转速n1下的流量、扬程与功率;a0、a1、a2、b0、b1、b2为模型参数。
所述流体输配管路阻力模型为:
式中,Hf为设备或管路阻力;χi为第i段管路沿程阻力系数;li为第i段管路长度; d为管路水力直径;ζj为第j个管道局部阻力件阻力系数;V为管道体积流量;为模型参数。
进一步的,本发明方法中既有中央空调系统的实测运行数据具体为:求解冷水机组性能预测模型的模型参数,即冷凝器的模型参数、蒸发器的模型参数和压缩机的模型参数,需要的实测数据有压缩机吸气温度teo、压缩机吸气压力pe、压缩机排气温度tci、压缩机排气压力pc、冷水机组负荷率ψ、压缩机功率Pin、冷却水流量mw,c、冷凝器进口水温twi,c、冷凝器出口水温two,c、冷冻水流量mw,e、蒸发器进口水温twi,e、蒸发器出口水温two,e;
求解表冷器性能预测模型的模型参数,需要的实测数据为:冷冻水流量mw,e、表冷器进口水温twi,b、表冷器出口水温two,b、表冷器进口干球温度tgi,b、表冷器进口湿球温度tsi,b、表冷器出口干球温度tgo,b、表冷器出口湿球温度tso,b;
求解冷却塔性能预测模型的模型参数,需要的实测数据为:冷却水流量mw,c、冷却塔进口水温twi,t、冷却塔出口水温two,t、冷却塔进口干球温度tgi,t、冷却塔进口湿球温度tsi,t、冷却塔出口干球温度tgo,t、冷却塔出口湿球温度tso,t;
求解所述水泵模型的模型参数时,需要的实测数据为:水泵在转速n0下的流量、扬程和功率;
求解所述流体输配管路阻力模型的模型参数,需要实测的数据为:冷冻水流量mw,e、冷却水流量mw,c、水泵扬程H1。
进一步的,本发明方法中,采用最小二乘法分别求解中央空调系统中冷水机组性能预测模型、表冷器性能预测模型、冷却塔性能预测模型、水泵模型和流体输配管路阻力模型的模型参数时,使用克莱姆法则保证线性方程组解的存在性和唯一性,最终求得模型参数的唯一解。
进一步的,本发明方法中,冷水机组性能预测模型、表冷器性能预测模型、冷却塔性能预测模型、水泵模型以及流体输配管路阻力模型均是基于集总参数法,将既有空调系统中各部件的未知结构参数进行集总建立得到。
本发明将理论分析与实际测试相结合,在难以获得空调系统各部件具体结构参数的条件下,提出将空调系统未知结构参数集总的空调系统特征识别方法,建立空调系统各部件性能预测模型,并根据实测数据采用最小二乘法对未知模型参数进行特征识别,为研究空调系统部分负荷下的性能变化规律及其运行控制优化提供依据,突破性的解决现有建模方法存在的弊端。
有益效果:本发明与现有空调系统建模方法相比,具有以下优点:
现有结构参数缺乏条件下的中央空调系统设备模拟方法,主要是通过海量的实际数据采集,并通过拟合的方法来建立仿真模型,其存在所需实测数据巨大,实测数据需覆盖所有运行工况,实测时间长等不足,同时还存在仿真精度较低,难以在不同空调系统中推广,适用性差等局限。本发明将理论分析与实际测试相结合,在难以获得空调系统各部件具体结构参数的条件下,提出将空调系统各部件未知结构参数集总的特征识别方法,建立空调系统各部件性能预测模型,并根据实测数据采用最小二乘法对各部件由于结构参数未知所导致的未知模型参数进行特征识别,并由空调系统各部件性能预测模型组合形成完整的空调系统性能预测模型。本发明的空调系统特征识别方法具有所需实测数据少,实测时间短,仿真精度高,适用性较广且所需实测参数量在实际空调系统中能够方便获得的特点。
附图说明
图1为本发明为求解既有中央空调系统中各部件模型参数而设置的数据测点布置图。
图2为本发明的空调系统特征识别方法操作流程图。
具体实施方式
下面结合实施例和说明书附图对本发明作进一步的说明。
图1为求解既有中央空调系统中各部件模型参数而设置的数据测点布置图。根据基于参数集总的建模原理,求解冷水机组性能预测模型、表冷器性能预测模型、冷却塔性能预测模型、水泵模型以及流体输配管理组力模型的模型参数,需要结合有关实测数据,为此按图1所示布置测点:在压缩机进、出口处分别布置温度及压力传感器,实时检测压缩机吸气温度teo、压缩机吸气压力pe、压缩机排气温度tci、压缩机排气压力pc;在冷凝器进、出口处分别布置温度传感器,实时检测冷凝器进口水温twi,c和出口水温two,c;在蒸发器进、出口处分别布置温度传感器,实时检测蒸发器进口水温twi,e和出口水温two,e;在冷凝器、蒸发器的出水管上分别布置流量传感器,实时检测冷却水流量mw,c和冷冻水流量mw,e;在表冷器进、出口处分别布置温度传感器,实时检测表冷器进口水温twi,b和出口水温two,b;测量表冷器进口干球温度tgi,b、表冷器进口湿球温度tsi,b、表冷器出口干球温度tgo,b、表冷器出口湿球温度tso,b;在冷却塔进、出口处分别布置温度传感器,实 时检测冷却塔进口水温twi,t和出口水温two,t;测量冷却塔进口干球温度tgi,t、冷却塔进口湿球温度tsi,t、冷却塔出口干球温度tgo,t、冷却塔出口湿球温度tso,t;测量压缩机功率Pin、冷冻水泵功率Pin,dp以及冷却水泵功率Pin,lp;测量冷水机组负荷率ψ、冷冻水泵转速ndp、冷却水泵转速nlp。
图2为空调系统特征识别方法操作流程图,第一步基于集总参数法,将既有空调系统中各部件的未知结构参数进行集总,分别建立中央空调系统中冷水机组性能预测模型、表冷器性能预测模型、冷却塔的性能预测模型、水泵模型、流体输配管路阻力模型,主要包括:
(1)冷水机组性能预测模型:
a.蒸发器模型在变水量工况下为:
Qe=mw,ecp,w(twi,e-two,e)=mr(heo-hei) (3)
蒸发器模型在变水温工况下为:
Qe=mw,ecp,w(twi,e-two,e)=mr(heo-hei) (8)
式中,Qe1、Δte1分别为两相区的换热量和换热温差;Qe2、Δte2分别为过热区的换热量和换热温差;mw,e为冷冻水流量;cp,w为水的比热;mr为制冷剂流量;A* 1,e为与蒸发器结构和污垢热阻有关的模型参数;x1为根据不同结构形式蒸发器的水侧换热系数经验公式得到的常系数;B* 1,e为与蒸发器结构和污垢热阻有关的模型参数;y1为根据不同结构形式蒸发器的制冷剂侧换热系数经验公式得到的常系数;C* 1,e为与蒸发器结构和污垢热阻有关的模型参数;A* 2,e为与蒸发器结构和污垢热阻有关的模型参数;B* 2,e为与蒸发器结构和污垢热阻有关的模型参数;C* 2,e为与蒸发器结构和污垢热阻有关的模型参数;y2为根据不同结构形式蒸发器的制冷剂侧换热系数经验公式得到的常系数;Qe为蒸发器总换热量;twi,e,two,e分别为冷冻水进口温度和出口温度;tw1,e为两相区冷冻水入口温度;hei,heo分别为蒸发器入口焓值和出口的焓值;te为蒸发温度;teo为压缩机吸气温度,即蒸发器出口的制冷剂温度。
b.冷凝器模型在变水量工况下为:
Qc=mw,ccp,w(two,c-twi,c)=mr(hci-hco) (11)
冷凝器模型在变水温工况下为:
Qc=mw,ccp,w(two,c-twi,c)=mr(hci-hco) (16)
式中,Qc为冷凝器换热量;mw,c为冷却水流量;twi,c,two,c分别为冷凝器的冷却水进口温度和出口温度;hci,hco分别为冷凝器入口焓值和出口焓值;Δtc1,Δtc2,Δtc3分别为冷凝器过冷区、两相区和过热区的换热温差;A* 1,c为需要确定的冷凝器模型参数;x2为根据不同结构形式冷凝器的水侧换热系数经验公式得到的常系数;B* 1,c为需要确定的冷凝器模型参数;y3为根据不同结构形式冷凝器的制冷剂侧换热系数经验公式得到的常系数;C* 1,c为需要确定的冷凝器模型参数;A* 2,c为需要确定的冷凝器模型参数;B* 2,c为需要确定的冷凝器模型参数;y4为根据不同结构形式冷凝器的制冷剂侧换热系数经验公式得到的常系数;C* 2,c为需要确定的冷凝器模型参数;A* 3,c为需要确定的冷凝器模型参数;B* 3,c为需要确定的冷凝器模型参数;C* 3,c为需要确定的冷凝器模型参数;y5为根据不同结构形式冷凝器的制冷剂侧换热系数经验公式得到的常系数;tc为冷凝温度;tco为制冷剂在冷凝器出口的温度;tw1,c,tw2,c分别为两相区冷却水进口温度和出口温度;tci为压缩机排气温度。
c.压缩机模型为:
Vth=ψVth0 (22)
式中,mr为制冷剂的质量流量;λ为输气系数;Vth为压缩机理论输气量;v1为压缩机吸气比容;ψ为冷水机组负荷率,针对变转速运行压缩机,ψ即为压缩机实际运行转速与额定转速之比;针对变频压缩机,ψ即为压缩机实际运行频率与额定频率之比;Vth0 为压缩机额定工况下理论输气量;Te为蒸发温度;Teo为吸气温度;Tci为排气温度;pe为蒸发压力;k为压缩过程多变指数;Pth为压缩机理论功率;Pin为压缩机实际功率;ηe为压缩机的电能效率;
其中,对于高速、多杠压缩机(n.≥720r/min,C=3%~4%)输气系数λ可采用下面经验公式:
电能效率ηe表示为:
ηe=ηi·ηm·ηd·ηmo
式中ηi、ηm、ηd、ηmo分别是压缩机的指示效率、摩擦效率、联轴器传动效率和电机效率。通常取ηe=0.4~0.55。
d.节流阀模型为:
hco=hei (26)
teo=te+△te (27)
tco=tc-△tc (28)
式中,hco为节流阀入口焓值;hei为节流阀出口焓值;Δte为过热度;Δtc为过冷度;
(2)表冷器性能预测模型包括表冷器换热量模型、表冷器传热效能模型、表冷器接触系数模型、表冷器传热单元数模型、热容比模型、出风参数模型:
表冷器换热量模型:
Qb=ma,b(hai,b-hao,b)=mw,bcp,w(two,b-twi,b) (29)
表冷器传热效能模型:
表冷器接触系数模型:
表冷器传热单元数模型:
热容比模型:(33)
出风参数模型,在干工况时为:
tgo,b=tgi,b-ε1,b(tgi,b-twi,b) (34)
出风参数模型,在湿工况时为:
tgo,b=tgi,b-ε2,b(tgi,b-tb) (36)
tso,b=tgo,b-(1-ε2,b)(tgi,b-tsi,b) (37)
出风参数模型,在临界工况时,采用上述干工况时、湿工况时的出风参数模型均可;此处临界工况定义为,出风温度正好达到出风空气的露点温度。
式中,Qb为表冷器换热量;ma,b为空气质量流量;hai,b为空气入口焓值;hao,b为空气出口焓值;mw,b为水流量;twi,b,two,b分别为表冷器进口水温和出口水温;ε1,b为表冷器传热效能;tgi,tgo分别为空气进口干球温度和出口干球温度;γ为热容比;NTU为表冷器传热单元数;ε2,b为表冷器接触系数;tso,b,tsi,b分别为空气出口湿球温度和进口湿球温度;ao,b为空气侧换热系数;Fb为表冷器总换热面积;cp,a为空气比热;Kb为表冷器总传热系数;tLi为空气进口露点温度;
(3)冷却塔性能预测模型,在冷却塔内空气在非饱和状态下时为:
所述冷却塔性能预测模型,在冷却塔内空气在饱和状态下时为:
式中,mw,t为冷却塔水质量流量;dz为冷却塔填料竖直方向微元长度;βt为冷却塔传质系数;a为填料比表面积;Fz为冷却塔填料横截面积;Xs,w为水温对应的饱和空气含湿量;X为冷却塔内空气水蒸气质量组分;Tw,t为水温;Le为刘易斯数;cp,v为水蒸气定压比热;r0为水的汽化潜热;Ta,t为冷却塔内空气温度;ma,t为冷却塔空气质量流量;Xs,a为饱和空气含湿量;(βtaFZ)j为换热过程中的等效值;Me为模型参数;mwi,t表示入口水流量;Ht表示填料高度;
(4)水泵模型为:
Hp0=a0+a1V0+a2V0 2 (50)
Pin,p0=b0+b1V0+b2V0 2 (51)
式中,V0、Hp0、Pin,p0分别为泵在转速n0下的流量、扬程与功率;V1、Hp1、Pin,p1分别为泵在转速n1下的流量、扬程与功率;a0、a1、a2、b0、b1、b2为模型参数。
(5)流体输配管路阻力模型为:
式中,Hf为设备或管路阻力;χi为第i段管路沿程阻力系数;li为第i段管路长度;d为管路水力直径;ζj为第j个管道局部阻力件阻力系数;V为管道体积流量;为模型参数,可由空调系统实测某流量下水泵扬程减去设备总阻力获得其值。
第二步为根据实测数据,采用最小二乘法对各部件性能预测模型的模型参数进行特征识别。通过图1中空调系统运行时相关测点检测的数据,利用最小二乘法拟合未知模型参数,得到表征各部件结构特性的特征参数,并由空调系统各部件性能预测模型组合形成完整的空调系统性能预测模型。
上述实施例仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和等同替换,这些对本发明权利要求进行改进和等同替换后的技术方案,均落入本发明的保护范围。
Claims (4)
1.一种空调系统特征识别方法,其特征在于,该方法包括以下步骤:
根据既有中央空调系统的实测运行数据,采用最小二乘法分别求解得到以下几种模型的模型参数:中央空调系统中冷水机组性能预测模型、表冷器性能预测模型、冷却塔性能预测模型、水泵模型、流体输配管路阻力模型,将所求得的模型参数作为表征空调系统各部件结构特性的特征参数;
所述冷水机组性能预测模型包括蒸发器模型、冷凝器模型、压缩机模型和节流阀模型:
a.所述蒸发器模型在变水量工况下为:
Qe=mw,ecp,w(twi,e-two,e)=mr(heo-hei) (3)
所述蒸发器模型在变水温工况下为:
Qe=mw,ecp,w(twi,e-two,e)=mr(heo-hei) (8)
式中,Qe1、Δte1分别为两相区的换热量和换热温差;Qe2、Δte2分别为过热区的换热量和换热温差;mw,e为冷冻水流量;cp,w为水的比热;mr为制冷剂流量;A* 1,e为与蒸发器结构和污垢热阻有关的模型参数;x1为根据不同结构形式蒸发器的水侧换热系数经验公式得到的常系数;B* 1,e为与蒸发器结构和污垢热阻有关的模型参数;y1为根据不同结构形式蒸发器的制冷剂侧换热系数经验公式得到的常系数;C* 1,e为与蒸发器结构和污垢热阻有关的模型参数;A* 2,e为与蒸发器结构和污垢热阻有关的模型参数;B*2,e为与蒸发器结构和污垢热阻有关的模型参数;C* 2,e为与蒸发器结构和污垢热阻有关的模型参数;y2为根据不同结构形式蒸发器的制冷剂侧换热系数经验公式得到的常系数;Qe为蒸发器总换热量;twi,e,two,e分别为冷冻水进口温度和出口温度;tw1,e为两相区冷冻水入口温度;hei,heo分别为蒸发器入口焓值和出口的焓值;te为蒸发温度;teo为压缩机吸气温度,即蒸发器出口的制冷剂温度;
b.所述冷凝器模型在变水量工况下为:
Qc=mw,ccp,w(two,c-twi,c)=mr(hci-hco) (11)
所述冷凝器模型在变水温工况下为:
Qc=mw,ccp,w(two,c-twi,c)=mr(hci-hco) (16)
式中,Qc为冷凝器换热量;mw,c为冷却水流量;twi,c,two,c分别为冷凝器的冷却水进口温度和出口温度;hci,hco分别为冷凝器入口焓值和出口焓值;Δtc1,Δtc2,Δtc3分别为冷凝器过冷区、两相区和过热区的换热温差;A* 1,c为需要确定的冷凝器模型参数;x2为根据不同结构形式冷凝器的水侧换热系数经验公式得到的常系数;B* 1,c为需要确定的冷凝器模型参数;y3为根据不同结构形式冷凝器的制冷剂侧换热系数经验公式得到的常系数;C* 1,c为需要确定的冷凝器模型参数;A* 2,c为需要确定的冷凝器模型参数;B* 2,c为需要确定的冷凝器模型参数;y4为根据不同结构形式冷凝器的制冷剂侧换热系数经验公式得到的常系数;C* 2,c为需要确定的冷凝器模型参数;A* 3,c为需要确定的冷凝器模型参数;B*3,c为需要确定的冷凝器模型参数;C* 3,c为需要确定的冷凝器模型参数;y5为根据不同结构形式冷凝器的制冷剂侧换热系数经验公式得到的常系数;tc为冷凝温度;tco为制冷剂在冷凝器出口的温度;tw1,c,tw2,c分别为两相区冷却水进口温度和出口温度;tci为压缩机排气温度;
c.所述压缩机模型为:
Vth=ψVth0 (22)
式中,mr为制冷剂的质量流量;λ为输气系数;Vth为压缩机理论输气量;v1为压缩机吸气比容;ψ为冷水机组负荷率,针对结构形式一定的压缩机,额定工况下的理论输气量为定值,当进行能量调节时,将实际工况下的理论输气量与额定工况下的理论输气量比值定义为冷水机组负荷率,用ψ表示;Vth0为压缩机额定工况下理论输气量;Teo为吸气温度;Tci为排气温度;pc为冷凝压力;pe为蒸发压力;k为压缩过程多变指数;Pth为压缩机理论功率;Pin为压缩机实际功率;ηe为压缩机的电能效率;
d.节流阀模型为:
hco=hei (26)
teo=te+△te (27)
tco=tc-△tc (28)
式中,hco为节流阀入口焓值;hei为节流阀出口焓值;Δte为过热度;Δtc为过冷度;
所述表冷器性能预测模型包括表冷器换热量模型、表冷器传热效能模型、表冷器接触系数模型、表冷器传热单元数模型、热容比模型、出风参数模型;
所述表冷器换热量模型为:
Qb=ma,b(hai,b-hao,b)=mw,bcp,w(two,b-twi,b) (29)
所述表冷器传热效能模型为:
所述表冷器接触系数模型为:
所述表冷器传热单元数模型为:
所述热容比模型为:
所述出风参数模型,在干工况时为:
tgo,b=tgi,b-ε1,b(tgi,b-twi,b) (34)
所述出风参数模型,在湿工况时为:
tgo,b=tgi,b-ε2,b(tgi,b-tb) (36)
tso,b=tgo,b-(1-ε2,b)(tgi,b-tsi,b) (37)
所述出风参数模型,在临界工况时,采用上述干工况时、湿工况时的出风参数模型均可;
式中,Qb为表冷器换热量;ma,b为空气质量流量;hai,b为空气入口焓值;hao,b为空气出口焓值;mw,b为水流量;twi,b,two,b分别为表冷器进口水温和出口水温;ε1,b为表冷器传热效能;γ为热容比;NTU为表冷器传热单元数;ε2,b为表冷器接触系数;tso,b,tsi,b分别为空气出口湿球温度和进口湿球温度;ao,b为空气侧换热系数;Fb为表冷器总换热面积;cp,a为空气比热;Kb为表冷器总传热系数;tLi为空气进口露点温度,tgi,b为表冷器进口干球温度,tgo,b为表冷器出口干球温度;
所述冷却塔性能预测模型,在冷却塔内空气在非饱和状态下时为:
所述冷却塔性能预测模型,在冷却塔内空气在饱和状态下时为:
式中,mw,t为冷却塔水质量流量;dz为冷却塔填料竖直方向微元长度;βt为冷却塔传质系数;a为填料比表面积;Az为冷却塔填料横截面积;Xs,w为水温对应的饱和空气含湿量;X为冷却塔内空气水蒸气质量组分;Tw,t为水温;Le为刘易斯数;cp,v为水蒸气定压比热;r0为水的汽化潜热;Ta,t为冷却塔内空气温度;ma,t为冷却塔空气质量流量;Xs,a为饱和空气含湿量;(βtaFZ)j为换热过程中的等效值;Me为模型参数;mwi,t表示入口水流量;Ht表示填料高度;
所述水泵模型为:
Hp0=a0+a1V0+a2V0 2 (50)
Pin,p0=b0+b1V0+b2V0 2 (51)
式中,V0、Hp0、Pin,p0分别为泵在转速n0下的流量、扬程与功率;V1、Hp1、Pin,p1分别为泵在转速n1下的流量、扬程与功率;a0、a1、a2、b0、b1、b2为模型参数;
所述流体输配管路阻力模型为:
式中,Hf为设备或管路阻力;χi为第i段管路沿程阻力系数;li为第i段管路长度;d为管路水力直径;ζj为第j个管道局部阻力件阻力系数;V为管道体积流量;为模型参数。
2.根据权利要求1所述的一种空调系统特征识别方法,其特征在于,所述既有中央空调系统的实测运行数据具体为:求解所述冷水机组性能预测模型的模型参数,即求解冷凝器的模型参数、蒸发器的模型参数和压缩机的模型参数时,需要的实测数据有压缩机吸气温度teo、压缩机吸气压力pe、压缩机排气温度tci、压缩机排气压力pc、冷水机组负荷率ψ;压缩机功率Pin、冷却水流量mw,c、冷凝器进口水温twi,c、冷凝器出口水温two,c、冷冻水流量mw,e、蒸发器进口水温twi,e、蒸发器出口水温two,e;
求解所述表冷器性能预测模型的模型参数时,需要的实测数据为:冷冻水流量mw,e、表冷器进口水温twi,b、表冷器出口水温two,b、表冷器进口干球温度tgi,b、表冷器进口湿球温度tsi,b、表冷器出口干球温度tgo,b、表冷器出口湿球温度tso,b;
求解所述冷却塔性能预测模型的模型参数时,需要的实测数据为:冷却水流量mw,c、冷却塔进口水温twi,t、冷却塔出口水温two,t、冷却塔进口干球温度tgi,t、冷却塔进口湿球温度tsi,t、冷却塔出口干球温度tgo,t、冷却塔出口湿球温度tso,t;
求解所述水泵模型的模型参数时,需要的实测数据为:水泵在转速n0下的流量、扬程和功率;
求解所述流体输配管路阻力模型的模型参数,需要实测的数据为:冷冻水流量mw,e、冷却水流量mw,c、水泵扬程H1。
3.根据权利要求1所述的一种空调系统特征识别方法,其特征在于,所述采用最小二乘法分别求解中央空调系统中冷水机组性能预测模型、表冷器性能预测模型、冷却塔性能预测模型、水泵模型和流体输配管路阻力模型的模型参数时,使用克莱姆法则保证线性方程组解的存在性和唯一性,最终求得模型参数的唯一解。
4.根据权利要求1、2或3所述的一种空调系统特征识别方法,其特征在于,所述冷水机组性能预测模型、表冷器性能预测模型、冷却塔性能预测模型、水泵模型和流体输配管路阻力模型均是基于集总参数法,将既有空调系统中各部件的未知结构参数进行集总建立得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610984289.4A CN106503388B (zh) | 2016-11-09 | 2016-11-09 | 空调系统特征识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610984289.4A CN106503388B (zh) | 2016-11-09 | 2016-11-09 | 空调系统特征识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106503388A CN106503388A (zh) | 2017-03-15 |
CN106503388B true CN106503388B (zh) | 2019-09-20 |
Family
ID=58323878
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610984289.4A Active CN106503388B (zh) | 2016-11-09 | 2016-11-09 | 空调系统特征识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106503388B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107388524A (zh) * | 2017-08-31 | 2017-11-24 | 广东美的制冷设备有限公司 | 空调器及其能效计算方法 |
CN110083952A (zh) * | 2019-04-30 | 2019-08-02 | 蒋甫政 | 二氧化碳列车空调仿真方法 |
CN111125938B (zh) * | 2020-01-15 | 2021-07-16 | 华南理工大学 | 基于次优算法的大型中央空调冷冻水管网优化设计方法 |
CN112818527B (zh) * | 2021-01-21 | 2022-07-29 | 中国科学院合肥物质科学研究院 | 面向用户端的闭式冷却塔冷却水出口温度预测方法 |
CN113654255B (zh) * | 2021-08-02 | 2022-08-30 | 北京京仪自动化装备技术股份有限公司 | 制冷系统、压缩机频率控制方法、电子设备及存储介质 |
CN113757908B (zh) * | 2021-09-28 | 2023-07-07 | 广东电网有限责任公司 | 一种空调系统的热特性测量方法、系统、终端及存储介质 |
WO2024046434A1 (zh) * | 2022-08-31 | 2024-03-07 | 青岛海信日立空调系统有限公司 | 冷水机组空调系统及其控制方法 |
CN115796055B (zh) * | 2023-01-10 | 2023-04-07 | 北京云庐科技有限公司 | 一种基于完整空调系统仿真模型的优化运行调节方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7096925B2 (en) * | 2003-06-19 | 2006-08-29 | General Motors Corporation | Modular electric HVAC systems for vehicles |
CN101021374A (zh) * | 2007-03-23 | 2007-08-22 | 东南大学 | 地源热泵的热水装置 |
CN101051005A (zh) * | 2006-07-10 | 2007-10-10 | 东南大学 | 多功能制冷空调试验平台 |
CN202813887U (zh) * | 2012-07-24 | 2013-03-20 | 上海伯涵热能科技有限公司 | 单双级压缩顺序使用的热泵、热泵空调及热泵热水机组 |
CN203274053U (zh) * | 2013-04-01 | 2013-11-06 | 冯益安 | 一种利用大面积平板式换热器组成的高效能空调 |
CN204373109U (zh) * | 2014-12-26 | 2015-06-03 | 东南大学 | 一种单冷水机组诱导送风与辐射一体化空调装置 |
CN204665479U (zh) * | 2015-06-04 | 2015-09-23 | 袁玉平 | 压缩机吸气管路结构 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6925826B2 (en) * | 2003-05-05 | 2005-08-09 | Carrier Corporation | Modular bus air conditioning system |
-
2016
- 2016-11-09 CN CN201610984289.4A patent/CN106503388B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7096925B2 (en) * | 2003-06-19 | 2006-08-29 | General Motors Corporation | Modular electric HVAC systems for vehicles |
CN101051005A (zh) * | 2006-07-10 | 2007-10-10 | 东南大学 | 多功能制冷空调试验平台 |
CN101021374A (zh) * | 2007-03-23 | 2007-08-22 | 东南大学 | 地源热泵的热水装置 |
CN202813887U (zh) * | 2012-07-24 | 2013-03-20 | 上海伯涵热能科技有限公司 | 单双级压缩顺序使用的热泵、热泵空调及热泵热水机组 |
CN203274053U (zh) * | 2013-04-01 | 2013-11-06 | 冯益安 | 一种利用大面积平板式换热器组成的高效能空调 |
CN204373109U (zh) * | 2014-12-26 | 2015-06-03 | 东南大学 | 一种单冷水机组诱导送风与辐射一体化空调装置 |
CN204665479U (zh) * | 2015-06-04 | 2015-09-23 | 袁玉平 | 压缩机吸气管路结构 |
Non-Patent Citations (2)
Title |
---|
空调系统变冷冻水温性能实验与预测模型;江楚遥等;《化工学报》;20141231;第65卷(第S2期);第265-271页 * |
部分负荷下热湿分段处理空调系统性能实验研究;姜小伟等;《制冷学报》;20141031;第35卷(第5期);第25-29、48页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106503388A (zh) | 2017-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106529021B (zh) | 一种基于特征识别的空调系统仿真模拟方法 | |
CN106503388B (zh) | 空调系统特征识别方法 | |
CN201569203U (zh) | 蒸冷式调温除湿机组 | |
CN101603751B (zh) | 一种制冷系统的变频节能控制方法 | |
CN1987440A (zh) | 多项热工综合性能测试实验的组合方法 | |
CN106871391A (zh) | 基于有限测点的空调系统性能在线检测方法 | |
CN105716340B (zh) | 一种基于多区域结霜图谱的空气源热泵除霜控制方法 | |
CN106839340A (zh) | 一种空调制冷量测量方法、装置及空调器 | |
CN102779217A (zh) | 一种结霜工况下的制冷系统计算机仿真性能计算方法 | |
CN105910220B (zh) | Vrf分户计量方法 | |
CN102353403B (zh) | 中央空调主机冷冻水流量及冷却介质流量测量方法 | |
CN102966524B (zh) | 制冷压缩机低吸气过热度性能测试装置 | |
Cheng et al. | Study of heat-source-tower heat pump system efficiency | |
CN107063735A (zh) | 一种微通道蒸发器的性能测试实验台 | |
CN102269661B (zh) | 一种风冷压缩冷凝机组性能试验系统 | |
CN104298888B (zh) | 一种基于流量冷量关系模型的风机盘管冷量计量方法 | |
CN109975050A (zh) | 一种板翅式换热器流动及传热性能测试装置及其方法 | |
CN105181313A (zh) | 一种热力膨胀阀性能对比实验台 | |
CN109543221A (zh) | 一种用于空调冷水机组现场故障诊断方法 | |
CN109026120B (zh) | 一种沸腾冷凝式矿井热湿环境模拟实验平台 | |
Gayeski et al. | Empirical modeling of a rolling-piston compressor heat pump for predictive control in low-lift cooling | |
CN103257051B (zh) | 空调压缩机量热计测试台 | |
CN207976235U (zh) | 一种多联式空调系统中换热量测量仪和采集装置 | |
CN203324212U (zh) | 用于确定锅炉烟气工程酸露点的装置 | |
CN109297619B (zh) | 一种导热油能量计量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |