CN104654690B - 冷水机组控制方法及控制系统 - Google Patents
冷水机组控制方法及控制系统 Download PDFInfo
- Publication number
- CN104654690B CN104654690B CN201410658830.3A CN201410658830A CN104654690B CN 104654690 B CN104654690 B CN 104654690B CN 201410658830 A CN201410658830 A CN 201410658830A CN 104654690 B CN104654690 B CN 104654690B
- Authority
- CN
- China
- Prior art keywords
- state
- space model
- water chilling
- module
- space
- 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
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F25—REFRIGERATION OR COOLING; COMBINED HEATING AND REFRIGERATION SYSTEMS; HEAT PUMP SYSTEMS; MANUFACTURE OR STORAGE OF ICE; LIQUEFACTION SOLIDIFICATION OF GASES
- F25B—REFRIGERATION MACHINES, PLANTS OR SYSTEMS; COMBINED HEATING AND REFRIGERATION SYSTEMS; HEAT PUMP SYSTEMS
- F25B49/00—Arrangement or mounting of control or safety devices
- F25B49/02—Arrangement or mounting of control or safety devices for compression type machines, plants or systems
Abstract
本发明专利涉及一种冷水机组控制方法及控制系统,该控制方法包括:能量守恒和质量守恒定律建立冷水机组的状态空间模型,所述状态空间模型表征为状态空间方程;利用基于支持向量机辨识模块对所述状态空间模型中的状态空间方程中的未知参数进行辨识。通过上述实施方式,使用状态空间方法建立冷水机组状态空间模型,并利用基于支持向量机辨识模块对模型的未知参数进行辨识,其模型精确度高,利于对冷水机组的精确控制。
Description
技术领域
本发明涉及一种冷水机组控制方法及控制系统。
背景技术
传统的冷水机组控制系统中常采用经典的传递函数模型,传递函数模型一般只适用于单输入-单输出的定常系统动态仿真,而空调系统中大部分部件是多输入—多输出系统,另外,传递函数模型采用的是比较抽象的频域分析方法,不便于在时域中进行数值计算。显然,传递函数法并不是理想的空调系统动态仿真建模手段。状态空间方法在系统动态建模过程中越来越受到人们的青睐,非常适合多输入-多输出控制系统,但空调系统往往具有高度的非线性和参数不确定性,这无疑使建立在精确数学模型基础上的状态空间模型面临实际应用上的挑战。
发明内容
本发明实施例所要解决的技术问题在于,提供一种冷水机组控制方法及控制系统,使用状态空间方法建立冷水机组状态空间模型,并利用基于支持向量机辨识模块对模型的未知参数进行辨识,其模型精确度高,利于对冷水机组的精确控制。
为解决上述技术问题,本发明实施例提供一种冷水机组控制方法,包括如下步骤:能量守恒和质量守恒定律建立冷水机组的状态空间模型,所述状态空间模型表征为状态空间方程;利用基于支持向量机辨识模块对所述状态空间模型中的状态空间方程中的未知参数进行辨识。
进一步地,在所述能量守恒和质量守恒定律建立冷水机组的状态空间模型的步骤之中,包括:将所述状态空间模型中的冷冻水与蒸发器的对流换热系数、制冷剂与蒸发器的对流换热系数、冷却水与冷凝器的对流换热系数以及制冷剂与冷凝器的对流换热系数设定为各自相关变量的多项式函数,再将多项式函数代入状态空间方程;在所述利用基于支持向量机辨识模块对所述状态空间模型中的状态空间方程中的未知参数进行辨识的步骤之中,包括:根据所述基于支持向量机辨识模块对所述多项式函数中的未知参数进行辨识以获取所述未知参数。
进一步地,在所述根据所述基于支持向量机辨识模块对所述多项式函数中的未知参数进行辨识以获取所述未知参数的步骤之后,包括:采集冷水机组的运行数据并实时输入所述状态空间模型以计算出所述冷水机组的输出参数值;将所述输出参数值与预先设定的参数值进行比较,并根据比较结果控制所述冷水机组。
为解决上述技术问题,本发明实施例还提供一种冷水机组控制系统,包括:状态空间模型输入输出模块,用于能量守恒和质量守恒定律建立冷水机组的状态空间模型,所述状态空间模型表征为状态空间方程;基于支持向量机辨识模块,用于对所述状态空间模型中的状态空间方程中的未知参数进行辨识。
进一步地,所述状态空间模型输入输出模块,还用于将所述状态空间模型中的冷冻水与蒸发器的对流换热系数、制冷剂与蒸发器的对流换热系数、冷却水与冷凝器的对流换热系数以及制冷剂与冷凝器的对流换热系数设定为各自相关变量的多项式函数,再将多项式函数代入状态空间方程;所述基于支持向量机辨识模块,还用于对所述多项式函数中的未知参数进行辨识以获取所述未知参数。
进一步地,所述控制系统还包括数据采集模块、比较模块及控制器;所述数据采集模块,用于采集冷水机组的运行数据,并输入所述状态空间模型输入输出模块建立的状态空间模型;所述比较模块,用于将所述状态空间模型输入输出模块根据输入的运行数据计算出的输出参数值与预先设定的参数值进行比较;所述控制器,用于根据所述比较模块得到的比较结果控制所述冷水机组。
本发明实施例的有益效果是:通过提供基于支持向量机辨识模块对冷水机组状态空间模型未知参数进行辨识,从而实现空调系统的多输入-多输出控制系统,达到对冷水机组动态精确控制的目的,并能清晰呈现冷水机组的状态参数。
附图说明
图1是本发明冷水机组控制方法一实施例的流程图。
图2是本发明冷水机组控制方法一具体实施例的工作结构流程图。
图3是本发明冷水机组控制方法一具体实施例的流程图。
图4是本发明冷水机组控制系统的结构示意图。
图5是本发明冷水机组控制系统另一具体实施例的结构示意图。
具体实施方式
结合图1-图3参阅,本发明实施例的冷水机组控制方法,包括如下步骤:
步骤S1,能量守恒和质量守恒定律建立冷水机组的状态空间模型,状态空间模型表征为状态空间方程。
具体而言,在步骤S1中,需要将所述状态空间模型中的冷冻水与蒸发器的对流换热系数、制冷剂与蒸发器的对流换热系数、冷却水与冷凝器的对流换热系数以及制冷剂与冷凝器的对流换热系数设定为各自相关变量的多项式函数,再将多项式函数代入状态空间方程。其中,冷水机组各部件状态空间方程表征如下。
冷凝器内制冷剂能量守恒:
冷却水质量守恒方程:
Gcw,L=Gcw,E=Gcw (2)
冷凝器内冷却水能量守恒方程:
冷凝器传热能量守恒方程:
蒸发器内制冷剂能量守恒方程:
蒸发器内冷冻水质量守恒方程:
Gew,L=Gew,E=Gew (6)
蒸发器内冷冻水能量守恒方程:
蒸发器传热能量守恒方程:
压缩机功率模型:
膨胀阀能量方程:
hr,3=hr,4 (10)
冷水机组制冷量:
Qc=Grm(hr,1-hr,4)=Grm(fh(tk)-fh(tc))=GrmfQc(tk,tc) (11)
冷水机组COP(Coefficient of Performance,性能系数):
基本变量(θ)考虑由初始值与与其变化量之和:
θ=θo+Δθ (13)
集中参数(σ)使用一阶泰勒函数:
因此等式(1)-(12)可得:
ΔGcw,L=ΔGcw,E=ΔGcw (21)
ΔGew,L=ΔGew,E=ΔGew (22)
ΔNcom=XEw,1Δtc+XEw,2Δtk+XEw,3ΔGrm (23)
ΔQc=XQ,1Δtc+XQ,2Δtk+XQ,3ΔGrm (24)
ΔCOP=Xcop,1Δtc+Xcop,2Δtk (25)
其中,Tcr=(ccrMcr)o;Tcw=cwMcw2;Tcg=ccgMcg;
Ter=(cerMer)o;Tew=cwMew2;Teg=cegMeg;
根据等式(15)~(25),冷水机组动态模型可以表述如下:
Δychiller=Cchiller·Δxchiller+Dchiller·Δuchiller (27)
令:
ΔXchiller=Δxchiller+Achiller -1ξchiller,这样等式能转换成一标准状态方程:
Δychiller=CchillerΔXchiller+DchillerΔuchiller-CchillerAchiller -1ξchiller (29)
基础参数如下:
aco=C0+C1(tc-tcg)+C1(tc-tcg)2 (30)
aeo=D0+D1(teg-tk)+D2(teg-tk)2 (31)
aew=F0+F1Gew+F2Gew 2 (32)
acw=E0+E1Gcw+E2Gcw 2
上述描述中,C0,C1,C2,D0,D1,D2,F0,F1,F2,E0,E1,E2,为需要辨识的系数。其中,
Δxchiller=[Δtc,Δtcw,L,Δtcg,Δtk,Δtew,L,Δteg]T;
Δychiller=[Δtcw,L,ΔGcw,L,Δtew,L,ΔGew,L,ΔNcom,ΔQc,ΔCOP]T;
Δuchiller=[Δtcw,E,ΔGcw,E,Δtew,E,ΔGew,E,ΔGrm]T;
这里Achiller,Bchiller和Cchiller为状态输入和输出矩阵,Dchiller是直接耦合从输入到输出矩阵,Δxchiller,Δuchiller和Δychiller表示各自变量的量化量。
上述实施例中,各符号和下标具体示意如下表(一)和表(二)。
A | 面积,m2 | u | 速度,m/s |
a | 对流换热系数,W/(m2.℃) | N | 电量,W |
c | 质量比热容,J/(kg·℃) | Nu | Nusselt常数 |
d | 水力直径,m | Re | Reynolds Number |
G | 质量流量kg/s | Z | 压缩因子 |
h | 焓,J/kg | τ | 时间,s |
ho | 初始焓kJ/kg | ρ | 密度,kg/m3 |
M | 质量,kg | κs | 收缩因子 |
p | 压力,Pa | η | 效率 |
Q | 传热量,W | Δte,shr | 过热温度,℃ |
qr | 制冷剂潜热,kJ/kg | λ | 导热系数,W/(m·℃) |
Rr | 制冷剂气态比热容,J/(kg·K) | μ | 动力系数,kg/(m·s) |
t | 温度,℃ | ν | 粘滞系数,m3/kg |
T | 温度,K |
表(一):符号说明
atm | 标准大气压 | eg | 蒸发器壁面 |
bl | 饱和液体制冷剂 | ew | 蒸发器冷量 |
bg | 饱和气态制冷剂 | E | 进口 |
c | 冷凝或制冷热容 | k | 蒸发 |
cc | 临界状态 | L | 出口 |
cg | 冷凝器壁面 | lΔt | 制冷剂过冷 |
co | 冷凝器与蒸发器之间 | o | 初始 |
cr | 冷凝器中的制冷剂 | r | 制冷剂 |
cw | 冷却液 | rm | 循环制冷剂 |
com | 压缩机 | g | 气体状态 |
er | 蒸发器中的制冷剂 | f | 液体状态 |
eo | 蒸发器与制冷剂之间 | q | 蒸发过程 |
表(二):下标说明
步骤S2,利用基于支持向量机辨识模块对状态空间模型中的状态空间方程中的未知参数进行辨识。具体为该基于支持向量机辨识模块对多项式函数中的未知参数进行辨识以获取所述未知参数,参阅下文。
该基于支持向量机辨识模块算法拟采用ε损失函数的支持向量回归机(ε-SVR)的算法对空调部件模型进行辨识,其实现步骤为:
①已知样本的训练集T={(x1,y1),…(xl,yl)}∈(χ×γ)l,其中,
②选择线性核函数K(x,x’)和适当的惩罚函数C>0;
K(xi,xj)=(x·xT) (33)
③构造式(1)~(3),并求解最优问题:
约束条件
得到最优解
④计算辨识系数w:
上述步骤①-④中,符号说明:K(x,x’):基于支持向量机核函数;C:惩罚函数;α:对偶变量,为Lagrange乘子;αi:对偶变量α的第i个分量;l:训练点个数;w:辨识系数(权向量);Minimize:最小化。
参阅图4,本发明实施例还提供一种冷水机组控制系统。该控制系统具体包括如下模块。
状态空间模型输入输出模块10,用于能量守恒和质量守恒定律建立冷水机组的状态空间模型,状态空间模型表征为状态空间方程。
基于支持向量机辨识模块11,用于对状态空间模型中的状态空间方程中的未知参数进行辨识。基于支持向量机辨识模块11是结构风险最小化原理,而非经验风险最小化原理,克服了传统学习算法“对有限量的样本拟合精度越高,其模型的推广能力越差”的缺陷,从而保证了学习机器智能辨识的泛化能力,有利于提高状态空间模型的适应性。
具体而言,该状态空间模型输入输出模块10,还用于将状态空间模型中的冷冻水与蒸发器的对流换热系数、制冷剂与蒸发器的对流换热系数、冷却水与冷凝器的对流换热系数以及制冷剂与冷凝器的对流换热系数设定为各自相关变量的多项式函数,再将多项式函数代入状态空间方程;该基于支持向量机辨识模块11,还用于对多项式函数中的未知参数进行辨识以获取未知参数。
参阅图5,进一步地,控制系统还包括数据采集模块12和比较模块13;数据采集模块12,用于采集冷水机组的运行数据,并输入状态空间模型输入输出模块10建立的状态空间模型;比较模块13,用于将状态空间模型输入输出模块10根据输入的运行数据计算出的输出参数值与预先设定的参数值进行比较;控制器14,用于根据比较模块13得到的比较结果控制所述冷水机组。
当进行控制时,可选取状态空间模型输入输出模块10输出的一个或多个参数作为主控参数,其他输出参数通过反馈矩阵直接反馈给冷水机组的控制器14,达到状态反馈控制的目的。
综上所述,本发明实施例的冷水机组控制方法及控制系统,通过提供基于支持向量机辨识模块对冷水机组状态空间模型未知参数进行辨识,从而实现空调系统的多输入-多输出控制系统,达到对冷水机组动态精确控制的目的,并能清晰呈现冷水机组的状态参数。
Claims (4)
1.一种冷水机组控制方法,其特征在于,包括如下步骤:
能量守恒和质量守恒定律建立冷水机组的状态空间模型,所述状态空间模型表征为状态空间方程;
利用基于支持向量机辨识模块对所述状态空间模型中的状态空间方程中的未知参数进行辨识;
在所述能量守恒和质量守恒定律建立冷水机组的状态空间模型的步骤之中,包括:将所述状态空间模型中的冷冻水与蒸发器的对流换热系数、制冷剂与蒸发器的对流换热系数、冷却水与冷凝器的对流换热系数以及制冷剂与冷凝器的对流换热系数设定为各自相关变量的多项式函数,再将多项式函数代入状态空间方程;
在所述利用基于支持向量机辨识模块对所述状态空间模型中的状态空间方程中的未知参数进行辨识的步骤之中,包括:根据所述基于支持向量机辨识模块对所述多项式函数中的未知参数进行辨识以获取所述未知参数。
2.根据权利要求1所述的控制方法,其特征在于,在所述根据所述基于支持向量机辨识模块对所述多项式函数中的未知参数进行辨识以获取所述未知参数的步骤之后,包括:
采集冷水机组的运行数据并实时输入所述状态空间模型以计算出所述冷水机组的输出参数值;
将所述输出参数值与预先设定的参数值进行比较,并根据比较结果控制所述冷水机组。
3.一种冷水机组控制系统,其特征在于,包括:
状态空间模型输入输出模块,用于能量守恒和质量守恒定律建立冷水机组的状态空间模型,所述状态空间模型表征为状态空间方程;
基于支持向量机辨识模块,用于对所述状态空间模型中的状态空间方程中的未知参数进行辨识;
所述状态空间模型输入输出模块,还用于将所述状态空间模型中的冷冻水与蒸发器的对流换热系数、制冷剂与蒸发器的对流换热系数、冷却水与冷凝器的对流换热系数以及制冷剂与冷凝器的对流换热系数设定为各自相关变量的多项式函数,再将多项式函数代入状态空间方程;
所述基于支持向量机辨识模块,还用于对所述多项式函数中的未知参数进行辨识以获取所述未知参数。
4.根据权利要求3所述的控制系统,其特征在于:
所述控制系统还包括数据采集模块、比较模块及控制器;
所述数据采集模块,用于采集冷水机组的运行数据,并输入所述状态空间模型输入输出模块建立的状态空间模型;
所述比较模块,用于将所述状态空间模型输入输出模块根据输入的运行数据计算出的输出参数值与预先设定的参数值进行比较;
所述控制器,用于根据所述比较模块得到的比较结果控制所述冷水机组。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410658830.3A CN104654690B (zh) | 2014-11-18 | 2014-11-18 | 冷水机组控制方法及控制系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410658830.3A CN104654690B (zh) | 2014-11-18 | 2014-11-18 | 冷水机组控制方法及控制系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104654690A CN104654690A (zh) | 2015-05-27 |
CN104654690B true CN104654690B (zh) | 2017-01-11 |
Family
ID=53246111
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410658830.3A Active CN104654690B (zh) | 2014-11-18 | 2014-11-18 | 冷水机组控制方法及控制系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104654690B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110567101A (zh) * | 2019-08-14 | 2019-12-13 | 南京南瑞继保工程技术有限公司 | 一种基于支持向量机模型的冷水机组高能效控制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0883046B1 (en) * | 1997-06-05 | 2002-06-26 | Carrier Corporation | Prediction of chiller compressor motor overheating |
CN102283283A (zh) * | 2011-05-30 | 2011-12-21 | 广西大学 | 智能化浅层地能低温储粮系统 |
CN103869695A (zh) * | 2014-02-26 | 2014-06-18 | 杭州电子科技大学 | 印染车间空气温度控制方法 |
CN103912966A (zh) * | 2014-03-31 | 2014-07-09 | 武汉科技大学 | 一种地源热泵制冷系统优化控制方法 |
-
2014
- 2014-11-18 CN CN201410658830.3A patent/CN104654690B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0883046B1 (en) * | 1997-06-05 | 2002-06-26 | Carrier Corporation | Prediction of chiller compressor motor overheating |
CN1098512C (zh) * | 1997-06-05 | 2003-01-08 | 运载器有限公司 | 冷却器压缩机电机过热的预测方法 |
CN102283283A (zh) * | 2011-05-30 | 2011-12-21 | 广西大学 | 智能化浅层地能低温储粮系统 |
CN103869695A (zh) * | 2014-02-26 | 2014-06-18 | 杭州电子科技大学 | 印染车间空气温度控制方法 |
CN103912966A (zh) * | 2014-03-31 | 2014-07-09 | 武汉科技大学 | 一种地源热泵制冷系统优化控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104654690A (zh) | 2015-05-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Nunes et al. | Modeling, simulation and optimization of a vapor compression refrigeration system dynamic and steady state response | |
CN103718127A (zh) | 用于优化冷却水设备系统操作的装置和方法 | |
Oronel et al. | Heat and mass transfer in a bubble plate absorber with NH3/LiNO3 and NH3/(LiNO3+ H2O) mixtures | |
Wan et al. | A new method for prediction and analysis of heat and mass transfer in the counter-flow dew point evaporative cooler under diverse climatic, operating and geometric conditions | |
CN1987440A (zh) | 多项热工综合性能测试实验的组合方法 | |
Gazda et al. | The estimation of energy efficiency for hybrid refrigeration system | |
JP2009150640A (ja) | パッケージ型空調機による空調システムの冷房能力測定方法 | |
Zhang et al. | Experimental investigation on the geometric structure with perforated baffle for thermal stratification of the water tank | |
Saito et al. | Experimental and numerical performance evaluation of double-lift absorption heat transformer | |
Sutthivirode et al. | Experimental investigation of vapour compression chiller based on transient cooling performance influenced by expansion devices | |
CN104654690B (zh) | 冷水机组控制方法及控制系统 | |
Ohkura et al. | Numerical analysis on performance enhancement of a CO2 heat pump water heating system by extracting tepid water | |
CN109140628A (zh) | 一种直膨式水冷冷热风机组控制方法及其温度调节系统 | |
Albieri et al. | Advanced control systems for single compressor chiller units | |
CN104374591A (zh) | 一种电子膨胀阀的能力测量装置 | |
Kuznetsov et al. | Thermoeconomic optimization of supercritical refrigeration system with the refrigerant R744 (CO2) | |
CN203573558U (zh) | 一种水源热泵模拟试验台 | |
Saifi et al. | CFD parametric investigation for two phase flow of refrigerant 134a in an adiabatic capillary tube | |
Sisti et al. | Assessment of basic control strategies through dynamic simulations: A CO2-based chiller under extreme off-design conditions | |
CN204555205U (zh) | 一种肉食品加工车间水冷空调的恒温控制系统 | |
Sian et al. | A generalized log-linear poisson-modeled correlation to predict the optimal heat rejection pressure of transcritical CO2 systems | |
Chinsuwan et al. | An optimization of a pre-cooling system for batch processes: A case study of the tubular ice production process | |
CN206919327U (zh) | 一种具有热回收功能的空调器 | |
Hladky et al. | Modelling and control of thermal system | |
Hmood et al. | Heat transfer analysis of roll-bond evaporator in a small-scale household refrigerator operating with R600a |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20221026 Address after: 518000 3AB, floor 3, building 2, shuntong'an science and technology plant, No. 2, Zijing Road, Futian Free Trade Zone, Fubao community, Fubao street, Futian District, Shenzhen, Guangdong Province Patentee after: Shenzhen Zhongke Xingyuan Technology Co.,Ltd. Address before: 518000 Guangdong city of Shenzhen province Nanshan District Xili Lake Patentee before: SHENZHEN POLYTECHNIC |