一种基于制冷性能曲线的智能控制方法与流程

本发明属于空调制冷
技术领域：
，尤其涉及一种基于制冷性能曲线的智能控制方法。
背景技术：
：目前应用于制冷空调系统冷水机组的控制方法主要有电流百分比的控制方式、供回水温差控制方式、压差控制和流量控制等。当前，国内外制冷空调群控技术主要应用在冷水机组群控方面，冷水机组的加减机控制逻辑有多种选择，目前的节能技术中以压缩机运行电流RLA(RatedLoadAmps)与额定电流的比值为依据的方式相对可靠节能。主要控制思路为，在系统负荷增大时，控制系统会自动对当前系统实际的冷冻水总供水温度与冷冻水供水温度设定值进行比较，并且自控系统会根据事先设定好的加载参数进行判断，如果满足加载条件，则自动控制系统会自动启动下一台机组以满足系统的需要。当冷冻机组的电流百分比低于设定的下限值，则减载该台冷冻机组。COP(CoefficientOfPerformance)即能量与热量之间的转换比率，简称能效比，是每台主机自身的属性。现有的空调系统群控方式中往往会忽视冷水机组COP值在不同工况下的变化，从而使得机组在低COP值下运行，普遍达不到很好的节能效果。公开号为CN104913559的中国专利申请提出一种基于主机COP值的制冷机组群控方法，重点强调一种在满足冷量需求前提下，可以充分节能的制冷机组群控方法。其主要根据冷冻水总管的供水或回水温度来判断加减载的，并采用了模糊PID控制器(Proportion,Integral,Derivative)复合控制温度的方式，使整个系统具有一定的稳定性，控制冷水机组在最高COP值点左右运行。该专利本具有一定的可行性并取得了一定的效果；但是从负荷分配的角度来看，要想让多台冷水机组都在最低的负载率下达到最高的COP值，使得能耗进一步降低，该专利还缺乏整体上的协调性。公开号为CN204853826的中国专利申请提出一种基于冷水主机COP能效曲线优化控制系统，该控制系统使用“COP优化控制系统”进行信号采集及冷水主机在不同运行环境下的COP能效曲线值进行比较，找出冷水主机最佳COP能效曲线及其最佳运行环境。优化冷水主机运行效率，降低能耗。但是该专利并未提供具体的实施方案与内容，而是否能够达到节能也没有进一步的验证，还有待更深入的研究与分析。技术实现要素：技术问题：本发明根据现有的控制方式的不足之处，在保证安全可靠以及制冷量的要求的情况下，提出了一种基于制冷性能曲线的智能控制方法，即根据冷水机组的性能曲线实时合理控制机组的启停及负荷分配。本发明包括如下步骤：步骤1，测量当前制冷空调系统冷水机组性能系数COP(CoefficientOfPerformance)值和负载率LR(LoadRate)值，进行曲线拟合得到COP值和LR值的关系模型；步骤2，验证关系模型的准确性，如果关系模型满足条件，执行步骤3，否则返回步骤1；步骤3，获取当前制冷空调系统总负荷，检测当前制冷空调系统总负荷在距上一次测量是否发生变化，如果是，执行步骤4，否则执行步骤6；步骤4，计算冷水机组的总能耗，并得到当前空调系统约束条件；步骤5，优化目标能耗方程，计算适应度值，判定是否满足终止条件，即适应度值是否达到全局最低，如果是，输出方程最优解，根据最优解中的负荷分配结果制定负荷分配控制方案，否则继续优化；步骤6，根据负荷分配控制方案对冷水机组进行控制。步骤7，根据供回水温度检测空调系统的稳定性，如果空调系统稳定，空调系统按得到的分配方案持续运行，否则返回执行步骤1。步骤1包括：测量当前冷水机组COP值和LR值，进行曲线拟合得到如下COP值和LR值的关系模型：COPi＝ai+biLRi+ciLRi2，其中，ai,bi,ci是固定不变的的参数，COPi表示第i台机组的COP值，LRi表示第i台机组的LR值。步骤2中，采用如下公式验证关系模型的准确性：其中，K为辨识结果，h为观测的数据向量，z为随机干扰，d为种群数量，P为估计误差协方差阵，θ用于存储参数估计结果，γ为遗忘因子，作用是加强新的数据提供的信息量，逐渐削弱老的数据,防止数据饱和，Kd表示d粒子的辨识结果，Kd+1表示更新的下一代辨识结果，θd+1表示更新的下一代参数估计结果，Pd+1表示示更新的下一代估计误差协方差阵。实时获取制冷空调系统冷水机组COP值和LR值，利用该公式进行系数ai,bi,ci的回归，将结果代入步骤1的关系模型中并绘制成以COP值间接表示的最小能耗曲线，将系数回归得到的模型与实际测量并计算的到的模型进行误差比较，如果COP值与实际测量的误差在阈值T1范围内(T1取值一般为3％)，执行步骤3，否则返回步骤1。步骤4包括：步骤4-1，计算冷水机组的总能耗J：其中，A代表系统总负荷，Xi代表第i台冷水机组的制冷量占总负荷的百分比，n表示运行的冷水机组的台数；步骤4-2，约束条件为：0≤Xi≤1，步骤5包括：步骤5-1，初始化一个规模为m的粒子群，初始化过程如下：设定群体规模m；对任意的粒子i及其维度s，在[-xmax,xmax]内服从均匀分布产生xis；对任意的粒子i及其维度s，在[-vmax,vmax]内服从均匀分布产生vis；其中vmax表示最大速度，xmax表示搜索空间的最大值，xis表示对任意的粒子i及其维度s的搜索空间的值，vis表示对任意的粒子i及其维度s的速度；步骤5-2，计算每个粒子的适应度值：目标函数fun2为：y＝-(-20.34*x(1)^2+19.41*x(1)+0.71)-(-52.96*x(2)^2+29.51*x(2)+1.04)-(-185.75*x(3)^2+57.92*x(3)+0.69)；此处x表示一个函数的自变量，y表示应变量；适应度值计算公式为：1.0*pop(j,1)<0.683，表示根据计算所求的对象1的负载率小于0.863；1.0*pop(j,2)<0.367，表示根据计算所求的对象2的负载率小于0.367；1.0*pop(j,3)<0.216，表示根据计算所求的对象3的负载率小于0.216；1.0*pop(j,1)+1.0*pop(j,2)+1.0*pop(j,3)>＝1，表示所求的对象综合大于等于1；适应度值fitness(j)＝fun2(pop(j,:))，表示根据函数fun2计算适应度值fitness；Pop()表示程序中对象的位置，以便识别粒子。步骤5-3，对每个粒子将其适应度值和其经历过的最好的位置的适应度值进行比较，若该粒子适应度值相比于其它粒子较低，表示负荷分配方案的能耗较小，则将该粒子作为当前最好位置；步骤5-4，对每个粒子将其适应度值和全局经历过的最好的位置的适应度值进行比较，若该粒子适应度值优于全局经历过的最好的位置的适应度值，则将其作为当前的全局最好位置，否则舍弃该粒子；步骤5-5，根据计算公式：vis(t+1)＝vis(t)+c1r1s(t)(pis(t)-xis(t))+c2r2s(t)(pgs(t)-xis(t))和xis(t+1)＝xis(t)+vis(t+1)分别对粒子的速度和位置进行更新；下标is表示i粒子的s维度，p表示当前的粒子，r是介于(0,1)之间的随机数，学习因子C1、C2是粒子群算法中的两个初始参数；vis(t)表示第t次迭代时i粒子的s维度的速度，xis(t)表示第t次迭代时i粒子的s维度的的位置；pis(t)表示第t次迭代时当前粒子在s维度经历的最好位置；pgs(t)表示第t次迭代时当前种群中在s维度的的经历的最好位置。步骤5-6，如果满足在所有迭代次数中的适应度值最小，则输出解；否则返回到步骤5-2。步骤7中，采用如下方法检测空调系统的稳定性：根据如下公式计算供水温度指标α：其中，Mgi为当前空调系统制冷机组设定的供水温度，Ngi为当前空调系统制冷机组测量的供水温度，Mgi+1为下一时刻空调系统制冷机组设定的供水温度，Ngi+1为下一时刻空调系统制冷机组测量的供水温度；根据如下公式计算回水温度指标β：其中，为当前空调系统制冷机组设定的回水温度，为当前空调系统制冷机组测量的回水温度，为下一时刻空调系统制冷机组设定的回水温度，为下一时刻空调系统制冷机组测量的回水温度，供水温度指标α和回水温度指标β的值均为1，空调系统的供回水温度达到稳定。有益效果：本发明与现有技术相比，具有以下优点：现有的控制方式中主要根据冷冻水总管的供水或回水温度来判断加减载的，本专利在同样保障安全的前提下以现有的数据测量，根据用户侧的需求通过数学模型寻优计算来控制空调系统冷水机组的启停和负荷分配，利用了寻优算法即能保证用户侧的安全可靠运行，又能合理的计算负荷，保证了冷水机组的总负载率在安全的情况下降至最低，达到节能的目的。本发明利用了带遗忘因子的最小二乘法验证了模型的准确性，考虑了时域不同情况下对机组的影响，并通过曲线拟合技术，将冷水机组的性能转化为数学模型，再根据寻优计算的结果来对机组控制，PSO寻优算法方便简单易操作，使得工作下的冷水机组整体的COP值达到最高，最终能从根本上达到节能的目的。附图说明下面结合附图和具体实施方式对本发明做更进一步的具体说明，本发明的上述或其他方面的优点将会变得更加清楚。图1是冷水系统连接示意图。图2是负荷分配适应度曲线。图3是总能耗对比图。图4是本发明流程图。具体实施方式下面结合附图及实施例对本发明做进一步说明。如图4所示，本发明包括如下步骤：步骤1，测量当前制冷空调系统冷水机组COP值和LR值，进行曲线拟合得到COP值和LR值的关系模型；步骤2，验证关系模型的准确性，如果关系模型满足条件，执行步骤3，否则返回步骤1；步骤3，获取当前制冷空调系统总负荷，检测当前制冷空调系统总负荷在距上一次测量是否发生变化，如果是，执行步骤4，否则执行步骤6；步骤4，计算冷水机组的总能耗，并得到当前空调系统约束条件；步骤5，优化目标能耗方程，计算适应度值，判定是否满足终止条件，即适应度值是否达到全局最低，如果是，输出方程最优解，根据最优解中的负荷分配结果制定负荷分配控制方案(负荷分配控制方案具体指不同冷水机组分配的负载)，否则继续优化；步骤6，根据负荷分配控制方案对冷水机组进行控制。步骤7，根据供回水温度检测空调系统的稳定性，如果空调系统稳定，空调系统按得到的分配方案持续运行，否则返回执行步骤1。步骤1包括：一般情况下，制冷的空调系统有较多的台冷水机组组成，其规格也是不尽相同，其主要采用的是冷媒水供水温度的控制方法，假如系统中所有的冷水机组额定制冷量均相同，各台机组均匀分担的冷量；若机组额定制冷量不完全相同，则每一台机组按照其额定制冷量占运行机组总制冷量的比例来提供冷量。本发明提出的冷水机组负荷分配优化策略，是根据COP值和LR之间关系，建立冷水机组总能耗目标函数，然后求解得到各冷水机组所承担的负荷比率。通过了解数据机房需求的制冷量，对冷水机组的负荷进行分配。COP是冷水机组自己自身的属性，其变化为物理过程，其状态方程中的参数属于慢时变参数，随时间的变化较为缓慢。本文采用曲线拟合法进行计算，来得到COP与LR的关系。因为COP与LR之间不存在线性关系，所以通过观测的数据拟合出它们之间的关系曲线，对3台不同额定制冷量的冷水机组进行了实验，测量当前冷水机组COP值和LR值，进行曲线拟合得到如下COP值和LR值的关系模型：其中，ai,bi,ci是固定不变的的参数，COPi表示第i台机组的COP值，LRi表示第i台机组的LR值。步骤2包括：本发明采用的曲线拟合获取COP与LR关系式，拟合得到的系数ai,bi,ci是固定不变的。机组性能系数COP不仅与LR的有关，还与机组其它属性有关，但都是属于慢时变参数，为了使模型更为准确，采用了带遗忘因子的最小二乘法来验证模型的准确性，采用如下公式验证关系模型的准确性：其中，K为辨识结果，h为观测的数据向量，z为随机干扰，d为种群数量，P为估计误差协方差阵，θ用于存储参数估计结果，γ为遗忘因子，作用是加强新的数据提供的信息量,逐渐削弱老的数据,防止数据饱和.实时获取制冷空调系统冷水机组COP值和LR值，利用该公式进行系数ai,bi,ci的回归，将结果代入步骤1的关系模型中并绘制成以COP值间接表示的最小能耗曲线，将系数回归得到的模型与实际测量并计算的到的模型进行误差比较。通过实时获取制冷空调系统冷水机组COP值和LR值，利用上述的带遗忘因子的最小二乘法进行系数回归，取额定功率为4100kW的冷水机组为例，其回归系数结果如下表1所示：表1将取得的结果带入原方程并绘制成曲线，与实际的参数模型进行比较发现系统的误差很小。在一个稳定的空调系统下，采用慢时变系数回归获取COP值与实际相差2.37％，在误差允许的范围之内；但在系统不稳定的情况下，某些时段负荷的不规则震荡会使得数据偏离，所以出现数据异常时可采用多次多时段测量以减小错误的出现。步骤4包括：步骤4-1，计算冷水机组的总能耗J：其中，A代表系统总负荷，Xi代表第i台冷水机组的制冷量占总负荷的百分比，n表示运行的冷水机组的台数；步骤4-2，约束条件为：0≤Xi≤1，步骤5包括：步骤5-1，初始化一个规模为m的粒子群，初始化过程如下：设定群体规模m；对任意的粒子i及其维度s，在[-xmax,xmax]内服从均匀分布产生xis；对任意的粒子i及其维度s，在[-vmax,vmax]内服从均匀分布产生vis；其中vmax表示最大速度，xmax表示搜索空间的最大值；步骤5-2，计算每个粒子的适应度值：目标函数fun2为：y＝-(-20.34*x(1)2+19.41*x(1)+0.71)-(-52.96*x(2)2+29.51*x(2)+1.04)-(-185.75*x(3)^2+57.92*x(3)+0.69)；适应度值计算公式为：1.0*pop(j,1)<0.683，1.0*pop(j,2)<0.367，1.0*pop(j,3)<0.216，1.0*pop(j,1)+1.0*pop(j,2)+1.0*pop(j,3)>＝1，适应度值fitness(j)＝fun2(pop(j,:))；步骤5-3，对每个粒子将其适应度值和其经历过的最好的位置的适应度值进行比较，若该粒子适应度值相比于其它粒子较低，表示负荷分配方案的能耗较小，则将该粒子作为当前最好位置；步骤5-4，对每个粒子将其适应度值和全局经历过的最好的位置的适应度值进行比较，若该粒子适应度值优于全局经历过的最好的位置的适应度值，则将其作为当前的全局最好位置，否则舍弃该粒子；步骤5-5，根据计算公式：vis(t+1)＝vis(t)+c1r1s(t)(pis(t)-xis(t))+c2r2s(t)(pgs(t)-xis(t))和xis(t+1)＝xis(t)+vis(t+1)分别对粒子的速度和位置进行更新；下标is表示i粒子的s维度，p表示当前的粒子，r是介于(0,1)之间的随机数，学习因子C1、C2是粒子群算法中的两个初始参数；vis(t)表示第t次迭代时i粒子的s维度的速度，xis(t)表示第t次迭代时i粒子的s维度的的位置；pis(t)表示第t次迭代时当前粒子在s维度经历的最好位置；pgs(t)表示第t次迭代时当前种群中在s维度的的经历的最好位置。步骤5-6，如果满足在所有迭代次数中的适应度值最小，则输出解；否则返回到步骤5-2。步骤7中，采用如下方法检测空调系统的稳定性：根据如下公式计算供水温度指标α：其中，Mgi为当前空调系统制冷机组设定的供水温度，Ngi为当前空调系统制冷机组测量的供水温度，Mgi+1为下一时刻空调系统制冷机组设定的供水温度，Ngi+1为下一时刻空调系统制冷机组测量的供水温度；根据如下公式计算回水温度指标β：其中，为当前空调系统制冷机组设定的回水温度，为当前空调系统制冷机组测量的回水温度，为下一时刻空调系统制冷机组设定的回水温度，为下一时刻空调系统制冷机组测量的回水温度，供水温度指标α和回水温度指标β的值均为1，空调系统的供回水温度达到稳定。实施例以北京亦庄制冷工程项目空调系统为例，本工程装机容量较大，耗能较高，从运营效率、运营期成本以及节能减排等的方面考虑，本工程采用了水冷冷水机组作为集中空调的冷源，按照分期分批建设的规模和进度，现配置三台额定制冷量分别为4100kW、2200kW和1300kW的离心式冷水机组。冷水系统的连接示意图如图1：图中有：1-机房；2-机密空调等设备；3-冷冻回水；4-冷冻出水；5-冷冻站；6-冷水机组；7-冷却出水；8-冷却回水；9-冷却塔；10-蓄冷罐：取环境温度5℃、相对湿度为50％的情况下空调的复合负载率LR和COP值的关系为例，本次北京亦庄电信具体情况如下表：表2不同负载率下的COP值5℃4100kW2200kW1300kW负载率LRCOPCOPCOP0.152.22.42.20.34.33.94.00.454.44.24.50.65.14.85.00.755.35.15.10.954.94.64.8现北京亦庄制冷工程项目用户侧总需求量约为6000kW，如按正常开启一台4100kW和一台2200kW的冷水机组即可以满足要求。通过本文的分析采用PSO寻优计算获取更加优化的方案，具体如下：(1)经过计算与转化得到测试函数的表达式：其约束条件通过PSO算法仿真的具体步骤如下：步骤一：初始化一个规模为200的粒子群并设置好初始位置和速度；步骤二：根据计算求出每一个粒子的适应度值；步骤三：将每一个新测量的粒子的适应度与之前最好的适应度值进行对比，如果较好则替换；步骤四：更新粒子速度与位置；步骤五：如果计算的负载率符合之前的要求输出该解，否则返回步骤二。得出适应度曲线如图2：从图中可以看出PSO算法在带约束的函数极值寻优方面体现了较好的寻优能力，收敛速度也较快，且比较简单容易操作。运行算法取出三组数据和最初方案四的分配如下表：表3优化后的负荷分配表序号x1x2x310.53390.31570.173020.53730.30460.165030.52480.31890.174140.66660.33330在通过仿真得到的结果代入到最初始的方程进行计算COP值：表4冷水机组COP值从表中可以看出前三种放案的负荷分配使得冷水机组都在高COP值下运行，而方案四的COP值不太理想。再通过计算对比前后计算能源消耗可以得出，本次北京亦庄制冷冷水机组运行方案相对于最初的方案可以节能7.6％左右，对于空调系统长期运行的制冷来说，将会是非常可观的能源，计算结果如图3所示。本发明提供了一种基于制冷冷水机组性能曲线的群控方法，具体实现该技术方案的方法和途径很多，以上所述仅是本发明的优选实施方式，应当指出，对于本
技术领域：
的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。当前第1页1 2 3