CN116558049A - 一种基于中央空调负荷动态预测的***及优化控制方法 - Google Patents
一种基于中央空调负荷动态预测的***及优化控制方法 Download PDFInfo
- Publication number
- CN116558049A CN116558049A CN202310308357.5A CN202310308357A CN116558049A CN 116558049 A CN116558049 A CN 116558049A CN 202310308357 A CN202310308357 A CN 202310308357A CN 116558049 A CN116558049 A CN 116558049A
- Authority
- CN
- China
- Prior art keywords
- unit
- refrigerating
- chw
- water pump
- chilled water
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 182
- 238000005265 energy consumption Methods 0.000 claims abstract description 76
- 239000000498 cooling water Substances 0.000 claims abstract description 75
- 238000001816 cooling Methods 0.000 claims abstract description 54
- 238000005457 optimization Methods 0.000 claims abstract description 44
- 238000004378 air conditioning Methods 0.000 claims abstract description 33
- 238000005057 refrigeration Methods 0.000 claims abstract description 21
- 238000011217 control strategy Methods 0.000 claims abstract description 9
- 239000003570 air Substances 0.000 claims description 98
- 230000017525 heat dissipation Effects 0.000 claims description 15
- 238000010219 correlation analysis Methods 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 5
- 230000005855 radiation Effects 0.000 claims description 5
- 238000013528 artificial neural network Methods 0.000 claims description 4
- 239000012080 ambient air Substances 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000004140 cleaning Methods 0.000 claims description 3
- 238000007418 data mining Methods 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 claims description 3
- 238000000513 principal component analysis Methods 0.000 claims description 3
- 229920006395 saturated elastomer Polymers 0.000 claims description 3
- 239000008400 supply water Substances 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 2
- 238000010438 heat treatment Methods 0.000 abstract description 2
- 238000009423 ventilation Methods 0.000 abstract description 2
- 238000004422 calculation algorithm Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 238000012821 model calculation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F24—HEATING; RANGES; VENTILATING
- F24F—AIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
- F24F11/00—Control or safety arrangements
- F24F11/30—Control or safety arrangements for purposes related to the operation of the system, e.g. for safety or monitoring
- F24F11/46—Improving electric energy efficiency or saving
- F24F11/47—Responding to energy costs
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F24—HEATING; RANGES; VENTILATING
- F24F—AIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
- F24F11/00—Control or safety arrangements
- F24F11/62—Control or safety arrangements characterised by the type of control or by internal processing, e.g. using fuzzy logic, adaptive control or estimation of values
- F24F11/63—Electronic processing
- F24F11/64—Electronic processing using pre-stored data
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F24—HEATING; RANGES; VENTILATING
- F24F—AIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
- F24F11/00—Control or safety arrangements
- F24F11/70—Control systems characterised by their outputs; Constructional details thereof
- F24F11/80—Control systems characterised by their outputs; Constructional details thereof for controlling the temperature of the supplied air
- F24F11/83—Control systems characterised by their outputs; Constructional details thereof for controlling the temperature of the supplied air by controlling the supply of heat-exchange fluids to heat-exchangers
- F24F11/85—Control systems characterised by their outputs; Constructional details thereof for controlling the temperature of the supplied air by controlling the supply of heat-exchange fluids to heat-exchangers using variable-flow pumps
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B30/00—Energy efficient heating, ventilation or air conditioning [HVAC]
- Y02B30/70—Efficient control or regulation technologies, e.g. for control of refrigerant flow, motor or heating
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Fuzzy Systems (AREA)
- Mathematical Physics (AREA)
- Air Conditioning Control Device (AREA)
Abstract
本发明属于暖通空调及自动控制技术领域,尤其涉及一种基于中央空调负荷动态预测的***及优化控制方法,基于室内变量参数、室外变量参数和模型输入参数,建立冷负荷动态预测模型,用于获取下一个时刻的冷负荷预测值;基于中央空调制冷***涉及的能耗设备,建立***设备能耗模型;基于下一个时刻的冷负荷预测值作为***设备能耗模型的输入值,建立***综合能耗模型;基于***综合能效的参数约束关系和下一个时刻的冷负荷预测值,建立中央空调制冷***的控制优化模型,通过建立冷负荷动态预测模型、***设备能耗模型、控制优化模型,优化冷冻机房中冷水机组、冷冻水泵、冷却水泵、冷却塔的运行控制策略,能够实现综合能效最优。
Description
技术领域
本发明涉及暖通空调及自动控制技术领域,尤其涉及一种基于中央空调负荷动态预测的***及优化控制方法。
背景技术
大型公建空调能耗占比较高,控制优化目的是在满足末端负荷需求下,调整设备开启台数以及供水温度、温差、电机频率以及水流量等参数,以提升***整体能效。现有技术方案中,主要通过BA机房群控***实现远程开关机以及水泵风机变频等自动控制,缺少较精准的***能耗模型及优化控制方法,导致***在变工况运行状态下,各项设备联动控制反馈滞后,供需不匹配,仍需依赖人工操作,实际运行综合能效偏低。
综上所述目前亟需一种基于中央空调负荷动态预测的***及优化控制方法,能够协调自控***智能优化控制,实现制冷机房***供需精确匹配,高效节能运行。
发明内容
为了解决上述现有技术中存在的问题,本发明提供了一种基于中央空调负荷动态预测的***及优化控制方法,拟解决目前的大型公建空调能耗占比较高,各项设备联动控制反馈滞后,供需不匹配,仍需依赖人工操作,实际运行综合能效偏低等问题。
为解决上述技术问题,本发明采用的技术方案如下:
一种基于中央空调负荷动态预测的***及优化控制方法,基于室内变量参数、室外变量参数和模型输入参数,建立冷负荷动态预测模型,用于获取下一个时刻的冷负荷预测值;基于中央空调制冷***涉及的能耗设备,建立***设备能耗模型;将所述下一个时刻的冷负荷预测值作为***设备能耗模型的输入值,建立***综合能耗模型;基于***综合能效的参数约束关系和下一个时刻的冷负荷预测值,建立中央空调制冷***的控制优化模型,用于对中央空调的运行控制策略和参数进行优化。
优先的,所述室内变量参数包括室内人员、室内温度和照明及设备用电量;所述室外变量参数包括室外干球温度、湿度、室外自然风速和室外太阳辐射;所述模型输入参数包括数据挖掘清洗、基于室内外变量相关性分析确定影响因子和主成分分析确定输入参数。
优选的,采用相关性分析初步筛选室内变量参数、室外变量参数和模型输入参数,获得冷负荷动态预测模型实际输入参数,采用人工神经网络或支持向量回归两种预测方法,获得下一个时刻的冷负荷预测值。
优选的,所述***设备能耗模型包括制冷机组能耗模型、冷却塔能耗模型、冷冻水泵能耗模型和冷却水泵能耗模型。
优选的,制冷机组能耗模型中的参数包括:
蒸发器制冷量:Qe_chw=CP×mchw×(Tchw_t-Tchw_s),
式中,Qe_chw为单台制冷机组蒸发器制冷量,单位为kW;CP为冷冻水比热容;Tchw_r为冷冻水回水温度,单位℃;Tchw_s为冷冻水供水温度,单位℃;mchw为冷冻水循环流量,单位kg/s;冷凝器散热量:Qe_ch=CP×mcw×(Tcw_t-Tcw_s),
式中,Qe_ch为单台制冷机组冷凝器散热量,单位为kW;CP为冷却水比热容;Tcw_r为冷却水回水温度,单位为℃;Tcw_s为冷却水供水温度,单位为℃;mcw为冷却水循环流量,单位为kg/s;
根据热量平衡关系:Qe_ch=Qe_chw+Pch,
式中,Pch为单台制冷机组压缩机耗电量,单位为kW;
压缩机耗电量:Pch=Qe_chw/COPch,
式中,COPch为制冷机组制冷能效;
COPch函数关系式:
式中,为制冷机组荷载率;Tchw_s为蒸发器供水温度,单位为℃;Tcw_r为冷凝器回水温度,单位为℃;
COPch拟合曲线函数:式中,a0~a4为拟合函数回归系数;
多台制冷机组制冷量之和需满足末端冷负荷需求:∑Qe_chw=Qe0,
式中,Qe0为预测冷负荷,由冷负荷动态预测模型结果输入。
优选的,冷却塔散热量模型:根据经典的传热单元数(ε-NTU)方法进行计算,
式中,cs为平均饱和空气定压比热容kg/(kJ.℃);hs,w,i,hs,w,o分别为进出空气饱和焓值;tw,i,tw,o分别为进出空气干球温度;
式中,ma为风机质量风量;cpw为水的定压比热容;m为热容率;mw,i为进入冷却塔的水质量流量,单位为kg/s;
式中,εa为换热效率;ha,0,ha,i分别为进出空气焓值;hs,w,i为进空气饱和焓值;
NTU为传热单元数;m为热容比率;
循环水量和风机风量的关系式如下:
式中,c的范围为0.5~5;n的范围为-1.1~0.35;mcw为冷却水流量;NTU为传热单元数;
冷却塔风机能耗模型:Pfun=b0+b1ma+b2m2 a,
式中,Pfun为风机耗功率,kW;ma为风机质量风量,kg/s;b0~b2为拟合系数;
根据风机风量mda与频率f或转速n的关系转换得到:
式中,ma0为风机额定工况风量,单位为kg/s;ma1为风机运行工况风量,单位为kg/s;nfun0为风机额定工况转速,单位为r/min;nfun1为风机运行工况转速,单位为r/min;ffun0为风机额定工况频率,单位为HZ;ffun1为风机运行工况频率,单位为HZ;
式中,Pfun0为风机额定工况耗电功率,单位为kW;Pfun1为风机运行工况耗电功率,单位为kW;
优选的,冷冻水泵能耗模型:Pchwp=c0+c1mchwp+c2mchwp 2,式中,c0~c2为拟合系数;
式中,mchwp0为冷冻水泵额定工况水流量,单位为kg/s;mchwp1为冷冻水泵运行工况水流量,单位为kg/s;nchwp0为冷冻水泵额定工况转速,单位为r/min;nchwp1为冷冻水泵运行工况转速,单位为r/min;fchwp0为冷冻水泵额定工况频率,单位为HZ;ffun1为冷冻水泵运行工况频率,单位为HZ;
式中,Pchwp0为冷冻水泵额定工况耗电功率,单位为kW;Pchwp1为冷冻水泵运行工况耗电功率,单位为kW;
优选的,冷却水泵能耗模型:Pcwp=d0+d1mcwp+d2mcwp 2,式中,d0~d2为拟合系数,
Pcwp为冷却水泵耗电功率,单位为kW;mcwp为冷却水质量流量,单位为kg/s;
式中,mcwp0为冷冻水泵额定工况水流量,单位为kg/s;mcwp1为冷冻水泵运行工况水流量,单位为kg/s;ncwp0为冷冻水泵额定工况转速,单位为r/min;ncwp1为冷冻水泵运行工况转速,单位为r/min;fcwp0为冷冻水泵额定工况频率,单位为HZ;fcwp0为冷冻水泵运行工况频率,单位为HZ;
式中,Pcwp0为冷却水泵额定工况耗电功率,单位为kW;Pcwp1为冷却水泵运行工况耗电功率,单位为kW。
优选的,总能耗模型如下:
式中,Ptot为制冷机房***总能耗,单位为kW;Pch为制冷机组耗电量,单位为kW;Pfun为冷却塔风机耗电量,单位为kW;Pchwp为冷冻水泵耗电量,单位为kW;Pcwp为冷却水泵耗电量,单位为kW;N1~N4分别为制冷机组、冷却塔、冷冻水泵及冷却水泵的配置台数;Zi为设备开启状态,Zi=0或1,0为未开启,1为开启;
***综合能效模型为:COPss=Qe0/Ptot,式中,COPss为***综合能效;Qe0为中央空调的预测负荷,单位为kW;Ptot为制冷机房***总能耗,单位为kW。
优选的,中央空调制冷***的控制优化模型的优化目标设定:
其数学模型为:
式中,Ptot为制冷机房***总能耗,单位为kW;Pch为制冷机组耗电量,单位为kW;Pfun为冷却塔风机耗电量,单位为kW;Zi为设备开启状态,Zi=0或1,0为未开启,1为开启;Pchwp为冷冻水泵耗电量,单位为kW;Pcwp为冷却水泵耗电量,单位为kW;
maxCOPss=max(Qe0/Ptot);
式中,COPss为***综合能效,Qe0为中央空调的预测负荷,kW;Ptot为制冷机房***总能耗,kW;
约束条件设定:
中央空调***制冷量的约束:制冷量需满足冷负荷需求的约束:
式中,Qe_chw为单台制冷机组蒸发器制冷量,kW;Zi为设备开启状态,Zi=0或1,0为未开启,1为开启;Qe0为中央空调的预测负荷,kW;
中央空调***变频水泵及风机对应的调频范围需在设备运行规定范围内,约束为:
式中,ffun为风机频率,fchwp为冷冻水泵频率,fcwp为冷却水泵频率;
中央空调***的进出水温度约束:冷冻水水温需在控制器设定范围内,冷却水进水温度若太高会降低制冷机组能效,由于冷热塔与环境空气进行换热,所以要高于环境的湿球温度,一般按不小于3℃逼近度考虑,因此约束为:
式中,Tchw_s为制冷机组冷冻水供水温度,Tcw_s为制冷机组冷却水进水温度,Twb为冷却塔周边环境湿球温度,带有min是该水温设定的最小值,带有max是该水温设定的最大值。
中央空调***各设备之间的约束:冷负荷由末端用户产生,冷冻水在末端设备与室内空气换热产生制冷效果,冷冻水温度升高,通过水泵循环至制冷机组蒸发器侧换热,冷冻水温度降低,冷却水循环带走冷凝器散热,并通过冷却塔风机将热量散发至室外空气。
制冷机组蒸发器产生制冷量满足末端冷负荷需求:式中,Qe_chw为单台制冷机组蒸发器制冷量,单位为kW,Zi为设备开启状态,Zi=0或1,0为未开启,1为开启,Qe0为中央空调的预测负荷,kW;
蒸发器吸收冷冻水中的热量:Qe_chw=CP×mchw×(Tchw_r-Tchw_s),
式中,Qe_chw为单台制冷机组蒸发器制冷量,kW;CP为冷冻水比热容,取4.18kJ/kg.℃;Tchw_s为蒸发器供水温度,单位为℃;Tcw_r为冷凝器回水温度,单位为℃;
冷却水带走冷凝器散热量:Qe_chw+Pch=CP×mcw×(Tcw_r-Tcw_s),式中,Qe_chw为单台制冷机组蒸发器制冷量,单位为kW;CP为冷冻水比热容;mcw为冷却水质量流量,单位为kg/s;Tcw_s、Tcw_r为冷却水供回水温度,单位为℃;
冷却水泵与冷却塔之间的关系:Qe_chw+Pch=εama(hs,w,i-ha,i)=ma(ha,o-ha,i),式中,Qe_chw为单台制冷机组蒸发器制冷量,单位为kW;Pch为单台制冷机组压缩机耗电量,单位为kW;εa为换热效率;ma为风机质量风量,单位为kg/s;ha,o,ha,i分别为进出空气焓值;hs,w,i为进空气饱和焓值;
整体优化模型建立:结合优化目标和约束条件得到整体优化模型
maxCOPSs=max(Qe0/Ptot),
式中Ptot为制冷机房***总能耗,单位为kW;Pch为制冷机组耗电量,单位为kW;Pfun为冷却塔风机耗电量,单位为kW;Pchwp为冷冻水泵耗电量,单位为kW;Pcwp为冷却水泵耗电量,单位为kW;N1~N4分别为制冷机组、冷却塔、冷冻水泵及冷却水泵的配置台数;Zi为设备开启状态,Zi=0或1,0为未开,启,1为开启;
式中,Qe_chw为单台制冷机组蒸发器制冷量,kW;Zi为设备开启状态,Zi=0或1,0为未开启,1为开启;Qe0为中央空调的预测负荷,单位为kW;CP为冷冻水比热容;mchw为冷冻水质量流量,单位为kg/s;Tchw_s、Tchw_r为冷冻水供回水温度,单位为℃;Tcw_s、Tcw_r为冷却水供回水温度,单位为℃;Pch为单台制冷机组压缩机耗电量,单位为kW;mcw为冷却水质量流量,单位为kg/s;ha,o,ha,i分别为进出空气焓值;ffun为冷却塔风机频率,单位为HZ;fchwp为冷冻水泵频率,单位为HZ;fcwp为冷却水泵频率,单位为HZ。
本发明的有益效果包括:
本发明基于中央空调制冷***,通过建立冷负荷动态预测模型、***设备能耗模型、控制优化模型,协同楼宇自动控制***实际控制要求,优化冷冻机房中冷水机组、冷冻水泵、冷却水泵、冷却塔的运行控制策略,能够实现综合能效最优。
附图说明
图1为本发明原理示意图。
图2为本发明的冷负荷预测原理示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚,下面将结合本申请实施例中附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本申请一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本申请实施例的组件可以各种不同的配置来布置和设计。因此,以下对在附图中提供的本申请的实施例的详细描述并非旨在限制要求保护的本申请的范围,而是仅仅表示本申请的选定实施例。基于本申请的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本申请保护的范围。
下面结合附图1-附图2对本发明作进一步的详细说明:
一种基于中央空调负荷动态预测的***及优化控制方法,基于室内变量参数、室外变量参数和模型输入参数,建立冷负荷动态预测模型,用于获取下一个时刻的冷负荷预测值;基于中央空调制冷***涉及的能耗设备,建立***设备能耗模型;基于下一个时刻的冷负荷预测值作为***设备能耗模型的输入值,建立***综合能耗模型;基于***综合能效的参数约束关系和下一个时刻的冷负荷预测值,建立中央空调制冷***的控制优化模型,用于对中央空调的运行控制策略和参数进行优化,图1中,Tchw_s、Tchw_r为冷冻水供回水温度,m_chw为冷冻水流量;Tcw_s、Tcw_r为冷却水供回水温度,m_cw为冷却水流量;Tn_s、Tn_r为末端设备送回风温度,m_n为末端设备循环风量;h_ai、h_ao为冷却塔风机侧进出风焓值,m_a为冷却塔风机风量;Twb为室外空气湿球温度;Q_t,Q_t+1分别为t及t+1时刻的负荷预测值,对所完成的中央空调制冷***模型和优化模型利用遗传算法进行计算,最终得出最优控制策略,比较模型计算结果与实测数据在不同时间点的***总能耗及综合效率,验证中央空调的***模型和节能优化效果。基于***优化控制模型,可实现对***中冷冻水流量、冷却水流量、冷却塔风机及循环水泵电机的控制频率、冷冻水出水温度、冷却水进水温度,以及各类设备的启停进行优化控制。
所述室内变量参数包括室内人员变化、室内温度和照明及设备用电量;所述室外变量参数包括室外干球温度、湿度、室外自然风速和室外太阳辐射;所述模型输入参数包括数据挖掘清洗、基于室内外变量相关性分析确定影响因子和主成分分析确定输入参数。
采用相关性分析初步筛选室内变量参数、室外变量参数和模型输入参数,获得冷负荷动态预测模型实际输入参数,采用人工神经网络或支持向量回归两种预测方法,获得下一个时刻的冷负荷预测值,分析室内外各变量参数对冷、热负荷的影响,室外变量参数包括室外干球温度、风速、太阳辐射等,室内变量包括人数、逐时用电量等,同时考虑各类变量的历史值对冷、热负荷的影响,为节省信息采集成本,采用相关性分析初步筛选主要影响因子,获得模型实际输入参数,其次,采用人工神经网络或支持向量回归两种预测方法,获得未来时刻的冷负荷预测值,负荷预测结果作为***运行调节优化控制模型的输入值,以提升***综合能效。
所述***设备能耗模型包括制冷机组能耗模型、冷却塔能耗模型、冷冻水泵能耗模型和冷却水泵能耗模型。
模型2:建立***设备能耗模型
(1)制冷机组能耗模型发器制冷量:Qe_chw=CP×mchw×(Tchwr-Tchw_s) (1)
上式中,Qe_chw为单台制冷机组蒸发器制冷量,kW;CP为冷冻水比热容,取4.18kJ/kg.℃;Tchw_r为冷冻水回水温度,℃;Tchw_s为冷冻水供水温度,℃。
冷凝器散热量:Qe_ch=CP×mcw×(Tcw_r-Tcw_s) (2)
上式中,Qe_ch为单台制冷机组冷凝器散热量,kW;Cp为冷却水比热容,取4.18kJ/kg.℃;Tcw_r为冷却水回水温度,℃;Tcw_s为冷却水供水温度,℃。
根据热量平衡关系:Qe_ch=Qe_chw+Pch (3)
上式中,Pch为单台制冷机组压缩机耗电量,kW。
压缩机耗电量:Pch=Qe_chw/COPch (4)
上式中,COPch为制冷机组制冷能效,与设备结构及运行工况有关,具体可通过拟合COPch与设备荷载率、以及经过蒸发器、蒸发器前后的水温的函数关系式。
COPch函数关系式:COPch=f(φ,Tchw_s,Tcw_r) (5)
上式中,φ为制冷机组荷载率,Tchw_s为蒸发器供水温度,Tcw_r为冷凝器回水温度。
COPch拟合曲线函数:COPch=a0+a1φ+a2φ2+a3(Tchw_s-Tcw_r)+a4(Tchw_s-Tcw_r)2 (6)
上式中,a0~a4为拟合函数回归系数,具体数值可根据设备性能参数拟合得到。
多台制冷机组制冷量之和需满足末端冷负荷需求:∑Qe_chw=Qe0 (7)
上式中,Qe0为预测冷负荷,由动态负荷预测模型结果输入。
(2)冷却塔能耗模型
冷却塔散热量模型:冷热塔的散热特性与冷却塔结构、室外湿球温度、循环水量、冷却塔进水温度以及风机风量等多种因素有关,参考冷却塔“ε-NTU”计算模型,循环水量和风机风量的关系式如下:
上式中,c和n拟合参数,可通过厂商提供的冷却塔测试数据拟合得到,c的范围一般在0.5~5,而n的范围一般在-1.1~0.35,NTU为传热单元数,可根据经典的传热单元数(ε-NTU)方法进行计算,NTU根据式(8-1)~(8-4)计算得到。
上式中,εa为换热效率,m为热容率,cs为平均饱和空气定压比热容kg/(kJ.℃),cpw为水的定压比热容,取4.18kg/(kJ.℃),ha,o,ha,i分别为进出空气焓值,hs,w,i,hs,w,o分别为进出空气饱和焓值,tw,i,tw,o分别为进出空气干球温度,通过式8-1,计算出cs,将cs带入式8-2中,计算出m,将m带入式8-3中,计算出εa,将εa带入式8-4中,计算出NTU
冷却塔风机能耗模型:Pfun=b0+b1ma+b2m2 a (9)
上式中,Pfun为风机耗功率,kW;ma为风机质量风量,kg/s。b0~b2为拟合系数,可根据厂商提供的风量和风机能耗性能参数数据拟合得到。实际控制中是改变风机的供电频率来控制风机的转速,从而实现对风量的调节,调频范围一般为25~50Hz。由于实际风量难以测量,可根据风机风量ma与频率f或转速n的关系转换得到:
(3)冷冻水泵能耗模型
制冷***冷冻水泵一般采用离心泵,功率-流量特性关系可按二次多项式进行拟合:
冷冻水泵能耗模型:Pchwp=c0+c1mchw+c2mchw 2 (11)
上式中,c0~c2为拟合系数,根据厂家提供的设备性能参数可拟合得到。通过变频器改变冷冻泵电机的供电频率可控制其转速,以实现对水流量的调节,频率调节范围一般为30~50Hz。
(4)冷却水泵能耗模型
制冷***冷却水泵一般采用离心泵,功率-流量特性关系可按二次多项式进行拟合:
冷却水泵能耗模型:Pchp=d0+d1mcw+d2mcw 2 (13)
上式中,d0~d2为拟合系数,根据厂家提供的设备性能参数可拟合得到。通过改变冷却泵的供电频率可控制其转速,以实现对水流量的调节,频率调节范围为30~50Hz。
(5)***总的能耗模型
制冷机房***总能耗包括冷水机组、冷却塔风机、冷冻水泵及冷却水泵四类设备的总能耗,总能耗模型如下:
上式中,Ptot为制冷机房***总能耗,kW;Pch为制冷机组耗电量,kW;Pfun为冷却塔风机耗电量,kW;Pchwp为冷冻水泵耗电量,kW;Pcwp为冷却水泵耗电量,kW;N1~N4分别为制冷机组、冷却塔、冷冻水泵及冷却水泵的配置台数;Zi为设备开启状态,Zi=0或1,0为未开启,1为开启。
中央空调***的综合能效模型为:COPss=Qe0/Ptot (16)
上式中,COPss为***综合能效,Qe0为中央空调的预测负荷,kW;Ptot为制冷机房***总能耗,kW。
模型3:建立控制优化模型
中央空调***的总能耗是四类主要设备能耗的总和,不同设备的某个参数增大或减小,会带来***能耗的改变,通过这些参数的约束关系,建立中央空调制冷***的整体优化模型,再利用整体优化模型,对中央空调的运行控制策略和参数进行优化。根据优化目标、各参数的约束关系和设备之间的约束关系,建立优化模型:
优化目标设定:
中央空调***设备运行最优控制策略的目标是使得整个***在满足负荷需求和设备正常运行的前提下,总能耗最低,同时综合***效率最高,其数学模型为:
maxCOPss=max(Qe0/Ptot) (18)
(2)约束条件设定:
中央空调***制冷量的约束:制冷量需满足冷负荷需求的约束:
中央空调***变频水泵及风机对应的调频范围需在设备运行规定范围内,约束为:
中央空调***的进出水温度约束:冷冻水水温需在控制器设定范围内,冷却水进水温度若太高会降低制冷机组能效,由于冷热塔与环境空气进行换热,所以要高于环境的湿球温度,一般按不小于3℃逼近度考虑,因此约束为:
上式中,Tchw_s为制冷机组冷冻水供水温度,Tcw_s为制冷机组冷却水进水温度,Twb为冷却塔周边环境湿球温度,带有min是该水温设定的最小值,带有max是该水温设定的最大值。
中央空调***各设备之间的约束:冷负荷由末端用户产生,冷冻水在末端设备与室内空气换热产生制冷效果,冷冻水温度升高,通过水泵循环至制冷机组蒸发器侧换热,冷冻水温度降低,冷却水循环带走冷凝器散热,并通过冷却塔风机将热量散发至室外空气。
制冷机组蒸发器产生制冷量满足末端冷负荷需求:
蒸发器吸收冷冻水中的热量:Qe_chw=CP×mchw×(Tchw_r-Tchw_s) (23)
冷却水带走冷凝器散热量:Qe_chw+Pch=CP×mcw×(Tcw_r-Tcw_s) (24)
冷却水泵与冷却塔之间的关系:Qe_chw+Pch=εama(hs,w,i-ha,i)=ma(ha,o-ha,i) (25)
整体优化模型建立:结合优化目标和约束条件得到整体优化模型
maxCOPss=max(Qe0/Ptot), (27)
模型求解算法
对所完成的中央空调制冷***模型和优化模型利用遗传算法进行计算,最终得出最优控制策略,比较模型计算结果与实测数据在不同时间点的***总能耗及综合效率,验证中央空调的***模型和节能优化效果。基于***优化控制模型,可实现对***中冷冻水流量、冷却水流量、冷却塔风机及循环水泵电机的控制频率、冷冻水出水温度、冷却水进水温度,以及各类设备的启停进行优化控制。
以上所述实施例仅表达了本申请的具体实施方式,其描述较为具体和详细,但并不能因此而理解为对本申请保护范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请技术方案构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。
Claims (10)
1.一种基于中央空调负荷动态预测的***及优化控制方法,其特征在于,基于室内变量参数、室外变量参数和模型输入参数,建立冷负荷动态预测模型,用于获取下一个时刻的冷负荷预测值;基于中央空调制冷***涉及的能耗设备,建立***设备能耗模型;将所述下一个时刻的冷负荷预测值作为***设备能耗模型的输入值,建立***综合能耗模型;基于***综合能效的参数约束关系和下一个时刻的冷负荷预测值,建立中央空调制冷***的控制优化模型,用于对中央空调的运行控制策略和参数进行优化。
2.根据权利要求1所述的一种基于中央空调负荷动态预测的***及优化控制方法,其特征在于,所述室内变量参数包括室内人员、室内温度和照明及设备用电量;所述室外变量参数包括室外干球温度、湿度、室外自然风速和室外太阳辐射;所述模型输入参数包括数据挖掘清洗、基于室内外变量相关性分析确定影响因子和主成分分析确定输入参数。
3.根据权利要求2所述的一种基于中央空调负荷动态预测的***及优化控制方法,其特征在于,采用相关性分析初步筛选室内变量参数、室外变量参数和模型输入参数,获得冷负荷动态预测模型实际输入参数,采用人工神经网络或支持向量回归两种预测方法,获得下一个时刻的冷负荷预测值。
4.根据权利要求1所述的一种基于中央空调负荷动态预测的***及优化控制方法,其特征在于,所述***设备能耗模型包括制冷机组能耗模型、冷却塔能耗模型、冷冻水泵能耗模型和冷却水泵能耗模型。
5.根据权利要求1所述的一种基于中央空调负荷动态预测的***及优化控制方法,其特征在于,制冷机组能耗模型中的参数包括:
蒸发器制冷量:Qe_chw=CP×mchw×(Tchw_r-Tchw_s),
式中,Qe_chw为单台制冷机组蒸发器制冷量,单位为kW;CP为冷冻水比热容;Tchw_r为冷冻水回水温度,单位℃;Tchw_s为冷冻水供水温度,单位℃;mchw为冷冻水循环流量,单位kg/s;冷凝器散热量:Qe_ch=CP×mcw×(Tcw_r-Tcw_s),
式中,Qe_ch为单台制冷机组冷凝器散热量,单位为kW;CP为冷却水比热容;Tcw_r为冷却水回水温度,单位为℃;Tcw_s为冷却水供水温度,单位为℃;mcw为冷却水循环流量,单位为kg/s;
根据热量平衡关系:Qe_ch=Qe_chw+Pch,
式中,Pch为单台制冷机组压缩机耗电量,单位为kW;
压缩机耗电量:Pch=Qe_chw/COPch,
式中,COPch为制冷机组制冷能效;
COPch函数关系式:
式中,为制冷机组荷载率;Tchw_s为蒸发器供水温度,单位为℃;Tcw_r为冷凝器回水温度,单位为℃;
COPch拟合曲线函数:式中,a0~a4为拟合函数回归系数;
多台制冷机组制冷量之和需满足末端冷负荷需求:∑Qe_chw=Qe0,
式中,Qe0为预测冷负荷,由冷负荷动态预测模型结果输入。
6.根据权利要求1所述的一种基于中央空调负荷动态预测的***及优化控制方法,其特征在于,冷却塔散热量模型:根据经典的传热单元数(ε-NTU)方法进行计算,
式中,cs为平均饱和空气定压比热容kg/(kJ.℃);hs,w,i,hs,w,o分别为进出空气饱和焓值;tw,i,tw,o分别为进出空气干球温度;
式中,ma为风机质量风量;cpw为水的定压比热容;m为热容率;mw,i为进入冷却塔的水质量流量,单位为kg/s;
式中,εa为换热效率;ha,o,ha,i分别为进出空气焓值;hs,w,i为进空气饱和焓值;
NTU为传热单元数;m为热容比率;
循环水量和风机风量的关系式如下:
式中,c的范围为0.5~5;n的范围为-1.1~0.35;mcw为冷却水流量;NTU为传热单元数;
冷却塔风机能耗模型:Pfun=b0+b1ma+b2m2 a,
式中,Pfun为风机耗功率,kW;ma为风机质量风量,kg/s;b0~b2为拟合系数;
根据风机风量mda与频率f或转速n的关系转换得到:
式中,ma0为风机额定工况风量,单位为kg/s;ma1为风机运行工况风量,单位为kg/s;nfun0为风机额定工况转速,单位为r/min;nfun1为风机运行工况转速,单位为r/min;ffun0为风机额定工况频率,单位为HZ;ffun1为风机运行工况频率,单位为HZ;
式中,Pfun0为风机额定工况耗电功率,单位为kW;Pfun1为风机运行工况耗电功率,单位为kW。
7.根据权利要求6所述的一种基于中央空调负荷动态预测的***及优化控制方法,其特征在于,冷冻水泵能耗模型:Pchwp=c0+c1mchwp+c2mchwp 2,式中,c0~c2为拟合系数;
式中,mchwp0为冷冻水泵额定工况水流量,单位为kg/s;mchwp1为冷冻水泵运行工况水流量,单位为kg/s;nchwp0为冷冻水泵额定工况转速,单位为r/min;nchwp1为冷冻水泵运行工况转速,单位为r/min;fchwp0为冷冻水泵额定工况频率,单位为HZ;ffun1为冷冻水泵运行工况频率,单位为HZ;
式中,Pchwp0为冷冻水泵额定工况耗电功率,单位为kW;Pchwp1为冷冻水泵运行工况耗电功率,单位为kW。
8.根据权利要求6所述的一种基于中央空调负荷动态预测的***及优化控制方法,其特征在于,冷却水泵能耗模型:Pcwp=d0+d1mcwp+d2mcwp 2,式中,d0~d2为拟合系数,Pcwp为冷却水泵耗电功率,单位为kW;mcwp为冷却水质量流量,单位为kg/s;
式中,mcwp0为冷冻水泵额定工况水流量,单位为kg/s;mcwp1为冷冻水泵运行工况水流量,单位为kg/s;ncwp0为冷冻水泵额定工况转速,单位为r/min;ncwp1为冷冻水泵运行工况转速,单位为r/min;fcwp0为冷冻水泵额定工况频率,单位为HZ;fcwp0为冷冻水泵运行工况频率,单位为HZ;
式中,Pcwp0为冷却水泵额定工况耗电功率,单位为kW;Pcwp1为冷却水泵运行工况耗电功率,单位为kW。
9.根据权利要求1所述的一种基于中央空调负荷动态预测的***及优化控制方法,其特征在于,总能耗模型如下:
式中,Ptot为制冷机房***总能耗,单位为kW;Pch为制冷机组耗电量,单位为kW;Pfun为冷却塔风机耗电量,单位为kW;Pchwp为冷冻水泵耗电量,单位为kW;Pcwp为冷却水泵耗电量,单位为kW;N1~N4分别为制冷机组、冷却塔、冷冻水泵及冷却水泵的配置台数;Zi为设备开启状态,Zi=0或1,0为未开启,1为开启;
***综合能效模型为:COPss=Qe0/Ptot,式中,COPss为***综合能效;Qe0为中央空调的预测负荷,单位为kW;Ptot为制冷机房***总能耗,单位为kW。
10.根据权利要求9所述的一种基于中央空调负荷动态预测的***及优化控制方法,其特征在于,中央空调制冷***的控制优化模型的优化目标设定:
其数学模型为:
式中,Ptot为制冷机房***总能耗,单位为kW;Pch为制冷机组耗电量,单位为kW;Pfun为冷却塔风机耗电量,单位为kW;Zi为设备开启状态,Zi=0或1,0为未开启,1为开启;Pchwp为冷冻水泵耗电量,单位为kW;Pcwp为冷却水泵耗电量,单位为kW;
maxCOPss=max(Qe0/Ptot);
式中,COPss为***综合能效,Qe0为中央空调的预测负荷,kW;Ptot为制冷机房***总能耗,kW;
约束条件设定:
中央空调***制冷量的约束:制冷量需满足冷负荷需求的约束:
式中,Qe_chw为单台制冷机组蒸发器制冷量,kW;Zi为设备开启状态,Zi=0或1,0为未开启,1为开启;Qe0为中央空调的预测负荷,kW;
中央空调***变频水泵及风机对应的调频范围需在设备运行规定范围内,约束为:
式中,ffun为风机频率,fchwp为冷冻水泵频率,fcwp为冷却水泵频率;
中央空调***的进出水温度约束:冷冻水水温需在控制器设定范围内,冷却水进水温度若太高会降低制冷机组能效,由于冷热塔与环境空气进行换热,所以要高于环境的湿球温度,一般按不小于3℃逼近度考虑,因此约束为:
式中,Tchw_s为制冷机组冷冻水供水温度,Tcw_s为制冷机组冷却水进水温度,Twb为冷却塔周边环境湿球温度,带有min是该水温设定的最小值,带有max是该水温设定的最大值。
中央空调***各设备之间的约束:冷负荷由末端用户产生,冷冻水在末端设备与室内空气换热产生制冷效果,冷冻水温度升高,通过水泵循环至制冷机组蒸发器侧换热,冷冻水温度降低,冷却水循环带走冷凝器散热,并通过冷却塔风机将热量散发至室外空气。
制冷机组蒸发器产生制冷量满足末端冷负荷需求:式中,Qe_chw为单台制冷机组蒸发器制冷量,单位为kW,Zi为设备开启状态,Zi=0或1,0为未开启,1为开启,Qe0为中央空调的预测负荷,kW;
蒸发器吸收冷冻水中的热量:Qe_chw=CP×mchw×(Tchw_r-Tchw_s),
式中,Qe_chw为单台制冷机组蒸发器制冷量,kW;CP为冷冻水比热容,取4.18kJ/kg.℃;Tchw_s为蒸发器供水温度,单位为℃;Tcw_r为冷凝器回水温度,单位为℃;
冷却水带走冷凝器散热量:Qe_chw+Pch=CP×mcw×(Tcw_r-Tcw_s),式中,Qe_chw为单台制冷机组蒸发器制冷量,单位为kW;CP为冷冻水比热容;mcw为冷却水质量流量,单位为kg/s;Tcw_s、Tcw_r为冷却水供回水温度,单位为℃;
冷却水泵与冷却塔之间的关系:Qe_chw+Pch=εama(hs,w,i-ha,i)=ma(ha,o-ha,i),式中,Qe_chw为单台制冷机组蒸发器制冷量,单位为kW;Pch为单台制冷机组压缩机耗电量,单位为kW;εa为换热效率;ma为风机质量风量,单位为kg/s;ha,o,ha,i分别为进出空气焓值;hs,w,i为进空气饱和焓值;
整体优化模型建立:结合优化目标和约束条件得到整体优化模型
maxCOPss=max(Qe0/Ptot),
式中Ptot为制冷机房***总能耗,单位为kW;Pch为制冷机组耗电量,单位为kW;Pfun为冷却塔风机耗电量,单位为kW;Pchwp为冷冻水泵耗电量,单位为kW;Pcwp为冷却水泵耗电量,单位为kW;N1~N4分别为制冷机组、冷却塔、冷冻水泵及冷却水泵的配置台数;Zi为设备开启状态,Zi=0或1,0为未开,启,1为开启;
式中,Qe_chw为单台制冷机组蒸发器制冷量,kW;Zi为设备开启状态,Zi=0或1,0为未开启,1为开启;Qe0为中央空调的预测负荷,单位为kW;CP为冷冻水比热容;mchw为冷冻水质量流量,单位为kg/s;Tchw_s、Tchw_r为冷冻水供回水温度,单位为℃;Tcw_s、Tcw_r为冷却水供回水温度,单位为℃;Pch为单台制冷机组压缩机耗电量,单位为kW;mcw为冷却水质量流量,单位为kg/s;ha,o,ha,i分别为进出空气焓值;ffun为冷却塔风机频率,单位为HZ;fchwp为冷冻水泵频率,单位为HZ;fcwp为冷却水泵频率,单位为HZ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310308357.5A CN116558049A (zh) | 2023-03-27 | 2023-03-27 | 一种基于中央空调负荷动态预测的***及优化控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310308357.5A CN116558049A (zh) | 2023-03-27 | 2023-03-27 | 一种基于中央空调负荷动态预测的***及优化控制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116558049A true CN116558049A (zh) | 2023-08-08 |
Family
ID=87485121
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310308357.5A Pending CN116558049A (zh) | 2023-03-27 | 2023-03-27 | 一种基于中央空调负荷动态预测的***及优化控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116558049A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117062419A (zh) * | 2023-10-11 | 2023-11-14 | 北京科技大学 | 多末端供需匹配的数据中心冷源***参数优化方法及装置 |
CN117232097A (zh) * | 2023-11-09 | 2023-12-15 | 上海轻环能源科技有限公司 | 基于自学习融合模型中央空调制冷站优化控制方法及*** |
CN117515808A (zh) * | 2024-01-05 | 2024-02-06 | 深圳市博久科技工程有限公司 | 一种中央空调冷热源节能智控*** |
CN117557070A (zh) * | 2024-01-11 | 2024-02-13 | 江西南昌济生制药有限责任公司 | 能耗优化方法、装置及电子设备 |
CN117421885B (zh) * | 2023-10-18 | 2024-05-28 | 广东德尔智慧科技股份有限公司 | 基于事件驱动机制的冷冻水***优化控制方法 |
-
2023
- 2023-03-27 CN CN202310308357.5A patent/CN116558049A/zh active Pending
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117062419A (zh) * | 2023-10-11 | 2023-11-14 | 北京科技大学 | 多末端供需匹配的数据中心冷源***参数优化方法及装置 |
CN117062419B (zh) * | 2023-10-11 | 2023-12-19 | 北京科技大学 | 多末端供需匹配的数据中心冷源***参数优化方法及装置 |
CN117421885B (zh) * | 2023-10-18 | 2024-05-28 | 广东德尔智慧科技股份有限公司 | 基于事件驱动机制的冷冻水***优化控制方法 |
CN117232097A (zh) * | 2023-11-09 | 2023-12-15 | 上海轻环能源科技有限公司 | 基于自学习融合模型中央空调制冷站优化控制方法及*** |
CN117232097B (zh) * | 2023-11-09 | 2024-02-20 | 上海轻环能源科技有限公司 | 基于自学习融合模型中央空调制冷站优化控制方法及*** |
CN117515808A (zh) * | 2024-01-05 | 2024-02-06 | 深圳市博久科技工程有限公司 | 一种中央空调冷热源节能智控*** |
CN117515808B (zh) * | 2024-01-05 | 2024-06-04 | 深圳市博久科技工程有限公司 | 一种中央空调冷热源节能智控*** |
CN117557070A (zh) * | 2024-01-11 | 2024-02-13 | 江西南昌济生制药有限责任公司 | 能耗优化方法、装置及电子设备 |
CN117557070B (zh) * | 2024-01-11 | 2024-04-12 | 江西南昌济生制药有限责任公司 | 能耗优化方法、装置及电子设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116558049A (zh) | 一种基于中央空调负荷动态预测的***及优化控制方法 | |
EP1387988B1 (en) | Air-conditioning system | |
CN110748970B (zh) | 空调***以及空调***的控制方法 | |
CN110107989A (zh) | 基于冷冻水回水温度最佳设定点的小型定频冷水机组变水温控制方法 | |
CN110895016A (zh) | 一种基于模糊自适应的中央空调***节能群控方法 | |
AU2002310859A1 (en) | Air-conditioning system | |
CN106765860B (zh) | 一种核电站中央空调的控制***及方法 | |
CN111947275A (zh) | 地铁站风水联动节能优化方法 | |
CN108105919B (zh) | 一种干工况制冷的变频空调***及其控制方法 | |
CN105180538A (zh) | 变频压缩机功率分配装置及其运行方法 | |
CN112413823A (zh) | 一种需求响应模式下中央空调***的分布式能量优化管理方法 | |
Liang et al. | Energy-efficient air conditioning system using a three-fluid heat exchanger for simultaneous temperature and humidity control | |
CN212720195U (zh) | 一种基于***整体能效比cop最佳的冷却水***控制装置 | |
CN110131819B (zh) | 建筑节能空调***及其运行方法 | |
CN111706943A (zh) | 一种恒温层与空气热泵耦合的辐射空调*** | |
CN113218040A (zh) | 一种中央空调***能效提升控制方法 | |
CN110748965A (zh) | 空调***及空调***控制方法 | |
CN211146702U (zh) | 空调*** | |
CN114992731A (zh) | 一种仅利用地热能和空气能实现建筑制冷的***及控制方法 | |
CN204830613U (zh) | 变频压缩机功率分配装置 | |
CN111473482A (zh) | 一种水冷式中央空调冷却循环控制装置及其方法 | |
CN115796055B (zh) | 一种基于完整空调***仿真模型的优化运行调节方法 | |
CN211903360U (zh) | 一种采用高温热泵的冷热联供控制*** | |
CN217088476U (zh) | 一种用于严寒地区的vfd与ahp联合运行*** | |
CN220087791U (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 |