联合供冷系统的递阶优化方法
联合供冷系统的递阶优化方法
【技术领域】
[0001] 本发明涉及能源网调度领域，特别是一种由常规电制冷机和冰蓄冷系统构成的联 合供冷系统的递阶优化方法。
【背景技术】
[0002] 由于大型城市普遍存在用地紧张的问题，为了利用有限的土地面积尽可能增大城 市容量，高层建筑成为城市设计者们理想的选择。高层建筑往往具有独立的供冷网络，提供 整栋楼各种各样的用冷需求。随着各种新兴技术的出现以及分时电价政策的推广，仅仅使 用常规电制冷机组的供冷系统经济效益不能令人满意，为了提高高层建筑供冷系统的经济 效益，多供冷机组互补协同系统应用得越来越普遍。典型的多供冷机组互补协同系统一般 包括：冰蓄冷、常规电制冷机，冷热电三联供及地源热泵等，其中常规电制冷机和冰蓄冷一 般共同承担90%以上的冷负荷。由于能够自主控制储能与用能时段，冰蓄冷系统越来越受 到人们的青睐，有效地对其与常规电制冷机进行协调控制是确保该协同系统节能降耗并提 尚经济效益的关键。
[0003] 目前冰蓄冷系统与常规电制冷机组的联合控制多停留在定性层面，即在电价谷时 制冰蓄冷，在电价峰时融冰放冷，缺乏对整个系统定量的分析与控制，这样导致了其控制方 法并不优越，而且由于电制冷机从启动达到设定的供冷功率有一个动态过程，所以仅仅考 虑稳态的控制方法难以确保整个系统的实际供冷功率迅速跟上预测负荷。因此，需要研宄 一种优化控制方法使得系统在动态过程中保证总负荷的前提下高效经济运行。
[0004] 对于高层建筑供冷系统这样一个分布式系统，通常可采用一般稳态优化方法对其 进行调度，然而由于电制冷机动态响应问题，进一步采用分布式预测控制优化电制冷机供 冷功率动态性能。对于分布式系统的动态控制，目前已经由集中式向分布式控制结构转化， 因为集中式控制会导致方法复杂且计算量偏大。在实际工业应用中，模型预测控制（MPC) 可以通过被控对象的脉冲或阶跃响应序列预测对象未来变化并且可自然引入过程时滞，对 于多变量系统可以获得比常规控制方法更加优秀的控制效果。在大型系统中，传统的模型 预测控制衍生出分布式预测控制（DMPC)这种形式，将大型系统分解成多个子系统，子系统 与子系统之间数据相互交流，每个子系统利用自身以及其余子系统的数据来完成自身的 MPC问题的求解，显著降低了计算量并可以达到整个系统的最优状态。
[0005] 因此，本发明拟借鉴以上对耦合系统的协调方法，研宄目标耦合但子系统间不耦 合的分布式系统的动态优化，使得整个系统在确保提供总预测冷负荷的前提下高效经济运 行。
【发明内容】
[0006] 针对现有技术存在的上述不足，本发明提出了一种联合供冷系统的递阶优化方 法，该方法对电制冷机和蓄冰槽进行经济性建模并且运用混合整数规划的方法得到经济指 标最优条件下各电制冷机与蓄冰槽的稳态运行时的供冷功率；下层运用分布式预测控制方 法保证各电制冷机和蓄冰槽在动态过程中尽可能跟踪上层经济优化得出的稳态设定值并 且总供冷功率跟踪预测负荷，解决目标耦合但子系统不耦合的分布式预测控制问题。
[0007] 为实现上述的目的，本发明提供一种联合供冷系统的递阶优化方法，所述联合供 冷系统是由常规电制冷机和冰蓄冷系统构成的联合供冷系统，所述方法具体包括如下步 骤：
[0008] 步骤S1 :利用电制冷机运行数据拟合其耗电功率关于制冷功率和冷却水进口温 度的函数；
[0009] 步骤S2 :利用上一步骤中得到的电制冷机耗电功率关于其制冷功率和冷却水进 口温度的函数建立电制冷机供冷的经济模型并进而建立冰蓄冷系统供冷的经济模型，最终 形成稳态经济优化问题的目标函数和约束条件并进行求解，得到各电制冷机和蓄冰槽在经 济效益最优条件下的启停状态和供冷功率；
[0010] 步骤S3 :在上一步骤的基础上，设计目标耦合的协调分布式预测控制方法重新优 化整个供冷系统在每一采样时刻各电制冷机的实时供冷功率设定值以提高整个系统响应 负荷的动态性能；
[0011] 步骤S4 :利用迭代方法求解上述分布式预测控制问题，得到各电制冷机每一采样 时刻的供冷功率设定值。
[0012] 较佳地，所述的步骤S1中：离心式电制冷机最重要的性能参数能效比为制冷功率 与耗电功率的比率。工程实践中，需要根据实际的冷却水进口温度、耗电功率和制冷功率等 数据通过回归分析将耗电功率表示成为冷却水进口温度和制冷功率二元二次函数形式以 便利用数学规划工具进行后续步骤优化问题的求解。
[0013] 较佳地，所述步骤S2中：每台电制冷机有开启和关闭两个工作状态，故采用0-1离 散变量Y来描述电制冷机的工作状态。对于双工况电制冷机空调工况和制冰工况的区分则 需使用两个0-1离散变量来分别描述。电制冷机的制冷功率则用连续变量来描述。根据步 骤S1中进行的回归分析，对于常规电制冷机，其耗电功率可以表示成关于制冷功率与冷却 水进口温度的二元二次函数：
[0015] 于是电制冷机供冷产生的电费花销C为：
[0016] C = YXg(P0Ut,T) XpXt (2)
[0017] 其中p为供冷时段电价，t为供冷时段时间长度。Y为电制冷机启停标志， 为步骤S1中拟合出的电制冷机耗电功率关于制冷功率和冷却水进口温度的函 数。
[0018] 在忽略蓄冰槽储冷量随储存时间损耗的前提下，为了考量蓄冰槽供冷的花销，弓丨 入蓄冰槽冷量均价Ps概念，描述单位冷量的价格，这个价格与充冷量、充冷时段电价和双工 况电制冷机制冰工况输入输出关系相关，由下式定义：
[0020] 其中pinit为冷量均价初值，w init为蓄冰槽冷量初值，t in为冲冷时间，t °ut为放冷时 间，Pm为双工况电制冷机总制冰功率，P s为蓄冰槽供冷功率，Y 3为蓄冰槽供冷启停标志。由 蓄冰槽冷量均价的定义可以看出，在不考虑储冷量随时间损耗的前提下，均价与蓄冰槽是 否放冷无关，只在蓄冰槽冲冷时可能改变。
[0021] 蓄冰槽储冷量增量为：
[0022] A ff = tinP-t0UtPsYs (4)
[0023]目标函数为各种电制冷机空调工况供冷与蓄冰槽供冷所需电费之和：
[0025] p为供冷时段电价，t为供冷时段时间长度，Yi为第i台电制冷剂启停标志， 为第i台电制冷机耗电功率关于制冷功率和冷却水进口温度的函数，其中 为第i台电制冷机供冷功率，TMin为第i台电制冷机冷却水进口温度，蓄冰槽供冷启 停标志，ps为蓄冰槽冷量均价，P s为蓄冰槽供冷功率，Y s为蓄冰槽供冷启停标志。
[0026]约束条件包括供冷功率与预测负荷平衡、各电制冷机额定功率限制、蓄冰槽储冷 量上下限和冰蓄冷系统工作模式唯一：
[0028] 其中I为第i台电制冷剂启停标志，P。^为第i台电制冷机供冷功率，P S为蓄冰 槽供冷功率，Ys为蓄冰槽供冷启停标志，Pp为预测负荷，a iE (0, 1)为第i台电制冷机最 小输出功率比例，Pi Mted为第i台电制冷机额定输出功率，W为蓄冰槽储冷量，Wmax和Wmin为 蓄冰槽储冷量上下限；Y ira为双工况制冷机制冰工况启停标志，默认制冰工况时所有双工况 冷机同时开启，加快制冰速度；YaiM为第i台双工况制冷机空调工况启停标志，共N台双工 况制冷机；
[0029] 最后，在形成稳态经济优化问题的目标函数和约束条件的前提下对上述优化问题 进行求解，得到各电制冷机和蓄冰槽在经济效益最优条件下的启停状态和供冷功率。
[0030] 较佳地，步骤S3中，所述标耦合的协调分布式预测控制方法，具体如下：
[0031] 分布式控制系统按照空调工况电制冷机数量分为N个子系统，每个子系统由模型 预测控制（MPC)和相应控制的电制冷机组成，不同子系统MPC的数据使用通信器进行实时 交换完成计算任务，从而使得各台空调工况电制冷机供冷功率跟踪步骤S2中得到的经济 效益最优条件下的供冷功率并且总供冷功率跟踪各台空调工况电制冷机经济效益最优条 件下的供冷功率之和；
[0032] 电制冷机制冷功率动态过程模型传递函数如下：
[0034]t为一阶惯性环节时间常数，td为延时时间，则第s个子系统离散状态空间模型 形式为：
[0036] 其中us(k)为k时刻电制冷机输出功率设定值，即为步骤S3中优化问题的优化变 量，Ps(k)为k+1时刻电制冷机实际输出功率，xs(k) = [ps(k)ps(k+l)L ps(k+nd)]T为状态 变量，nd为电制冷机延时时间相对于离散系统采样时间的倍数;A s、Bs、Cs为状态空间方程中 的系数矩阵，其中
由离散系统采样时间和电制冷机动态过 程时间常数决定，R表示实数集，设采样时间为At，则有：
[0038] Bs= [0 L 0 l-e_At/T]T (9)
[0039] Cs= [1 0 L 0]
[0040] 对于第s个子系统，s = 1，2,…，N ;MPC协调策略包括电制冷机在未来P个采样时 刻组成的预测时域内预测输出功率尽可能接近步骤S2得到的单台空调工况电制冷机经 济 效益最优条件下的供冷功率r s，并且区别于一般的MPC策略，为了控制所有电制冷机预测输 出功率尽可能跟上各台空调工况电制冷机经济效益最优条件下的供冷功率之和r，在目标 函数中引入预测时域内各台空调工况电制冷机总预测功率与经济效益最优条件下的供冷 功率之和差距的平方和，因此不同子系统MPC目标函数出现耦合；综上，第s个子系统MPC 目标函数如下所示：
[0042] 其中，Us= [1130011；30^+1)~11；30^-1)]1为预测时域内电制冷机输出功率设定值 序列，u s(k)为k时刻电制冷机输出功率设定值，即为优化问题的优化变量；而Q，R G [0, 1] 为权重系数，分别用来衡量单台电制冷机预测负荷跟踪程度与所有电制冷机总负荷跟踪程 度的重要性；rs(k+i)为（k+i)时刻第s台空调工况电制冷机经济效益最优条件下的供冷功 率，r(k+i)为（k+i)时刻各台空调工况电制冷机在经济效益最优条件下的供冷功率之和； /\0 + /|幻为第s台空调工况电制冷机在k时刻通过预测模型预测的其在（k+i)时刻的输 出功率，
为除第s台空调工况电制冷机之外的其余空调工况电制冷机在k时 刻通过预测模型预测的其在（k+i)时刻的输出功率之和；
[0043] MPC预测模型为：
[0045] is(々+ M幻为在k时刻通过预测模型计算出的（k+i)时刻第S台空调工况电制冷 机状态空间方程中状态变量，+ ZI幻为第s台空调工况电制冷机在k时刻通过预测模 型预测的其在（k+i)时刻的输出功率;A/4和A /为状态空间方程中系数矩阵的相应次幂。
[0046] 对于第s个子系统MPC控制器，整个优化问题形式为：
[0049]其中，{us，min，us，max}，{p s，min，ps，J为控制变量的边界约束。
[0050] 较佳地，所述步骤S4具体包括：
[0051] 步骤S41 :在k时刻，根据步骤S3中得到的每台空调工况电制冷机经济效益最优 条件下的供冷功率rs，初始化每个子系统在预测时域内的电制冷机输出功率设定值序列 U s，令迭代次数1 = 0 :
[0053] 上述<〇 + 1)为第1次迭代在（k+1)时刻的电制冷机输出功率设定值，P为预测时 域长度；
[0054] 步骤S42:根据MPC预测模型，计算出每个子系统MPC控制时域内输出量预测值
，并通知给其他N-1个子系统MPC ;
[0055] 步骤S43:对于每个子系统，求解步骤S3中的优化问题，得到最优解，即最优制冷 机输出功率设定值序列Uf 1;
[0056] 步骤S44 :检查所有子系统MPC的收敛条件是否满足，即对给定的精度e sG R(s =1，2，L，N)，是否均有
如果所有的子系统收敛条件满足，令 每个子系统最优控制量a与最优制冷机输出功率设定值序列Uf"相等，即K = uf，转步 骤S45;否则，令
a为介于〇和1之间的常数，引入该 参数为了尽量避免方法不收敛的情况，按需求选取数值大小，1 = 1+1，转步骤S42 ;
[0057] 步骤S45:在k时刻选取控制量<=[丨0 L0]U:(.s'=l二…，A〇作用于相应的 子系统。
[0058] 步骤S46:滚动移位到下一个采样时刻，即k+1 - k，返回步骤S41，重复上述过程。
[0059] 与现有技术相比，本发明具有以下有益效果：
[0060] 本发明针对由常规电制冷机和冰蓄冷系统构成的联合供冷系统的递阶优化方法 中通过建立了电制冷机和蓄冰槽的经济模型，并分别使用混合整数规划和分布式预测控制 解决稳态优化问题和动态优化问题，定量地确定联合供冷系统的调度方法，使得整个系统 在确保提供总预测冷负荷的前提下高效经济运行，为供冷系统能源管理与优化调度提供合 理的建议。
【附图说明】
[0061] 通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、 目的和优点将会变得更明显：
[0062]图1为由常规电制冷机和冰蓄冷系统构成的联合供冷系统结构图；
[0063] 图2为本发明的一实施例的流程图；
[0064] 图3为本发明的一实施例的方法示意图；
[0065] 图4为本发明的一实施例的分布式预测控制迭代求解方法流程图；
[0066] 图5为本发明的一实施例中3900kW冷机回归分析图；
[0067] 图6为本发明的一实施例中2150kW冷机回归分析图；
[0068] 图7为本发明的一实施例中6329kW冷机回归分析图；
[0069] 图8为本发明的一实施例中负荷预测图；
[0070] 图9为本发明的一实施例中两种策略每半小时的电费对比图；
[0071] 图10为本发明的一实施例的动态性能优化仿真结果，其中：(a)为总体，（b)为第 一台3900kW电制冷机，（c)为第二台3900kW电制冷机，（d)为第一台2150kW电制冷机。
【具体实施方式】
[0072] 下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术 人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术 人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进。这些都属于本发明 的保护范围。下面结合附图来对本发明的由常规电制冷机和冰蓄冷系统构成的联合供冷系 统递阶优化方法作进一步详细的说明。
[0073] 如图1所示，为由常规电制冷机和冰蓄冷系统构成的联合供冷系统结构。
[0074] 图2为本发明的一个【具体实施方式】的由常规电制冷机和冰蓄冷系统构成的联合 供冷系统递阶优化方法的流程图。首先根据各冷源经济模型和分时电价政策，采用混合整 数规划方法，优化各冷源的开启状态及最优功率设定值，同时考虑到各冷源的动态调节特 性以及物理约束，为了更好地提高系统动态性能，设计了一种目标耦合的协调分布式预测 控制方法，在分布式框架下重新优化各冷源的设定值，使得各冷源在动态过程中保证总负 荷的同时尽可能跟踪最优制冷功率设定值。
[0075] 具体的，如图2所示，本发明的由常规电制冷机和冰蓄冷系统构成的联合供冷系 统递阶优化方法，包括以下的步骤：
[0076] 步骤S1 :利用电制冷机运行数据拟合其耗电功率关于制冷功率和冷却水进口温 度的函数；
[0077] 离心式电制冷机最重要的性能参数能效比（C0P)为制冷功率与耗电功率的比率。 工程实践中，需要根据实际的冷却水进口温度、耗电功率和制冷功率等数据通过回归分析 将耗电功率表示成为冷却水进口温度和制冷功率二元二次函数形式以便利用数学规划工 具进行后续步骤优化问题的求解。
[0078] 步骤S2 :建立电制冷机和冰蓄冷系统供冷的经济模型，形成稳态经济优化问题的 目标函数和约束条件并进行求解；
[0079] 每台电制冷机有开启和关闭两个工作状态，故采用0-1离散变量Y来描述电制冷 机的工作状态。对于双工况电制冷机空调工况和制冰工况的区分则需使用两个0-1离散变 量来分别描述。电制冷机的制冷功率则用连续变量来描述。根据步骤S1中进行的回归分 析，对于常规电制冷机，其耗电功率可以表示成关于制冷功率与冷却水进口温度的二元二 次函数：
[0081] 于是电制冷机供冷产生的电费花销C为：
[0082] C = YXg(P0Ut,T) XpXt (2)
[0083] 其中p为供冷时段电价，t为供冷时段时间长度。
[0084] 在忽略蓄冰槽储冷量随储存时间损耗的前提下，为了考量蓄冰槽供冷的花销，弓丨 入蓄冰槽冷量均价Ps概念，描述单位冷量的价格，这个价格与充冷量、充冷时段电价和双工 况电制冷机制冰工况输入输出关系相关，由下式定义：
[0086] 其中pinit为冷量均价初值，W init为蓄冰槽冷量初值，t in为冲冷时间，t °ut为放冷时 间，Pm为双工况电制冷机总制冰功率，P s为蓄冰槽供冷功率，Y 3为蓄冰槽供冷启停标志。由 蓄冰槽冷量均价的定义可以看出，在不考虑储冷量随时间损耗的前提下，均价与蓄冰槽是 否放冷无关，只在蓄冰槽冲冷时可能改变。
[0087] 蓄冰槽储冷量增量为：
[0088] A ff = tinP-t0UtPsYs (4)
[0089] 目标函数为各种电制冷机空调工况供冷与蓄冰槽供冷所需电费之和：
[0091] P。^为第i台电制冷机供冷功率，T Min为第i台电制冷机冷却水进口温度。
[0092] 主要的约束条件包括供冷功率与预测负荷平衡、各电制冷机额定功率限制、蓄冰 槽储冷量上下限和冰蓄冷系统工作模式唯一：
[0094] 其中Pp为预测负荷，a iG (〇, 1)为第i台电制冷机最小输出功率比例，P i Mtral为 第i台电制冷机额定输出功率，Wmax和Wmin为蓄冰槽储冷量上下限。Y 为双工况制冷机制 冰工况启停标志，默认制冰工况时所有双工况冷机同时 开启，加快制冰速度。YaiM为第i台 双工况制冷机空调工况启停标志，共N台双工况制冷机。
[0095] 步骤S3:在上一步稳态经济优化得到的各电制冷机最优制冷功率的基础上设计 一种目标耦合的协调分布式预测控制方法优化整个供冷系统响应负荷的动态性能；
[0096] 分布式控制系统按照空调工况电制冷机数量分为N个子系统，每个子系统由MPC 和相应控制的电制冷机组成，不同子系统MPC的数据使用通信器进行实时交换完成计算任 务从而使得各台空调工况电制冷机供冷功率跟踪步骤S2中计算出的期望值并且总供冷功 率跟踪总期望值，具体的分布式控制系统结构如图3所示。
[0097] 电制冷机制冷功率动态过程模型传递函数如下：
[0099] t为一阶惯性环节时间常数，td为延时时间，则第s个子系统离散状态空间模型 形式为：
[0101] 其中U(k)为k时刻电制冷机输出功率设定值，即为优化问题的优化变量，ps(k)为 k+1时刻电制冷机实际输出功率，nd为电制冷机延时时间相对于离散系统采样时间的倍数。 xs(k) = [ps(k)ps(k+l)L ps(k+nd)]T为状态变量。
由采样时间 和电制冷机时间常数决定，设采样时间为At，则有：
[0103] Bs= [0 L 0 l-e_At/T]T (9)
[0104] Cs= [1 0 L 0]
[0105] 对于第s个子系统（s = 1，2,…，N)，MPC协调策略包括电制冷机在未来P个采样 时刻组成的预测时域内预测输出功率尽可能接近步骤S2计算出的期望功率r s，并且区别于 一般的MPC策略，为了控制所有电制冷机预测输出功率尽可能跟上总期望功率r，本策略在 目标函数中引入预测时域内总预测功率与总期望功率差距的平方和，因此不同子系统MPC 目标函数出现耦合。综上，第s个子系统MPC目标函数如下所示：
[0107] 其中，Us= [u(k)u(k+l)-u(k+p-l)] T为预测时域内电制冷机输出功率设定值序 列，而Q，R e [0, 1]为权重系数，分别用来衡量单台电制冷机预测负荷跟踪程度与所有电制 冷机总负荷跟踪程度的重要性。
[0108] MPC预测模型为：
[0110]对于第s个子系统MPC控制器，整个优化问题形式为：
[0113]其中，{us，min，us，max}，{p s，min，ps，J为控制变量的边界约束。
[0114] 步骤S4 :利用迭代方法求解上述分布式预测控制问题，得到各电制冷机每一采样 时刻的供冷功率设定值；
[0115] 迭代方法流程图如图4所示。
[0116] 步骤S41 :在k时刻，根据步骤S3中求解出的每台电制冷机理想制冷功率rs初始 化每个子系统在预测时域内的U s，令迭代次数1 = 0，
[0118]步骤S42 :根据MPC预测模型，计算出每个子系统MPC控制时域内输出量预测值
，并通知给其他N-1个子系统MPC。
[0119] 步骤S43:每个子系统步骤S3中的优化问题，得到最优解lUi+1。
[0120] 步骤S44:检查所有子系统MPC的收敛条件是否满足，即对给定的精 度e sG R(s = 1，2，L，N)，是否均有
，如果所有的 子系统收敛条件满足，令每个子系统最优控制量M=U(/+I，转步骤S45;否则，令
，a为介于〇和1之间的常数，引入该参数为了尽量避 免方法不收敛的情况，按需求选取数值大小，1 = 1+1，转步骤S42。
[0121] 步骤S45:在k时刻选取控制量》:=[1 0 L 0]U〗（5 = 1,2,作用于相应的 子系统。
[0122] 步骤S46 :滚动移位到下一个采样时刻，即k+1 - k，返回步骤S41，重复上述过程。
[0123] 基于上述的实施过程，本发明以上海某超高层建筑低区供冷系统为例进行了仿真 验证。在该供冷系统中，常规电制冷机共三台，两台额定制冷功率3900kW，一台额定制冷功 率2150kW，双工况电制冷机共三台，空调工况额定制冷功率6329kW，制冰工况额定制冷功 率 3868kW。
[0124] 首先，按步骤S1根据电制冷机运行实测数据进行输入输出关系函数的拟合。
[0125] 图5,图6,图7分别为3900kW、2150kW和6329kW离心式电制冷机回归分析结果图。 各图中XI为冷机供冷功率标么值，X2为冷机冷却水进口温度，Y1为冷机耗电功率标幺值。 图中红色曲线包围的区域为置信水平99 %所对应的置信区间
[0126] 从图5-7图中可以看出，二元二次回归分析可以相当精确地建立电制冷机输入电 功率与输出冷功率和冷却水进口温度的函数关系，具体的函数表达式如下列各式，P in为电 制冷机耗电功率，电制冷机制冷功率，2：?为电制冷机冷却水进口温度。
[0127] 对于3900kW离心式电制冷机有：
[0129] 对于2150kW离心式电制冷机有：
[0131] 对于双工况电制冷机6329kW空调工况有：
[0133] 对于双工况电制冷机的制冰工况，为了加快制冰速度，在分配制冰的情况下电制 冷机全功率制冰，所以不存在部分功率制冰的工况，此时制冰功率cold和耗电功率elec分 别仅为关于冷却水进口温度凡^的一元二次函数：
[0135] 然后，按步骤S2进行稳态经济优化仿真，预测负荷和电价政策如下。如图8所示。
[0136] 表1分时电价政策
[0139] 为了经济优化策略节费效果，选取一种常规的稳态调度策略作为对比，该策略在 电价谷值时优先蓄冰槽冲冷，在电价平值与峰值时优先电制冷机供冷，各电制冷机按额定 功率由小到大的顺序启用。两种策略每半小时电费对比如图9。
[0140] 经过计算，稳态经济优化策略在给定的预测负荷条件下，6小时内总电费花销为 3970元，而对比策略总电费花销为4053元，6小时内本发明策略相对节费2%。
[0141] 随后，按照步骤S3和S4进行分布式预测控制优化策略设计和问题求解。动态优 化采样间隔为2. 5分钟，即步骤S2的稳态调度策略每30分钟计算一次各电制冷机与蓄冰 槽最优供冷功率后，而动态性能优化每2. 5分钟更新一次各电制冷机供冷功率设定值以确 保各电制冷机和蓄冰槽在动态过程中尽可能跟踪步骤S2经济优化得出的稳态设定值并且 总供冷功率跟踪预测负荷。冰槽供冷一阶模型时间常数设为2分钟，而各电制冷机空调工 况动态模型时间常数和延时时间设定如下表所示：
[0142] 表2空调工况电制冷机动态参数
[0144] 总体的和各台启用过的空调工况电制冷机的期望承担负荷与分别应用与不应用 DMPC策略导致的实时供冷功率曲线如组图10所示，由于在当前预测负荷情况下6329kW电 制冷机并未开启，所以组图中只包括3900kW和2150kW电制冷机。
[0145] 为了定量描述采用或不采用DMPC策略时整个系统以及每台启用过的空调工况电 制冷机实时供冷功率跟步骤S2中计算出的经济最优条件下的最优供冷功率接近程度的区 另IJ，计算了 6个小时内每个采样时刻偏差的平方和，数值如表3所示，可以看出采用动态性 能优化策略后偏差平方和明显减小，实现了总体系统和每个子系统的迅速响应，体现出良 好的控制效果。
[0146] 表3动态优化效果
[0148] 综合仿真结果，本发明提出的由常规电制冷机和冰蓄冷系统构成的联合供冷系统 的递阶优化方法兼具较好的经济价值和动态控制效果，对于实际供冷系统经济优化运行有 一定的指导意义。
[0149] 综上所述，本发明的针对由常规电制冷机和冰蓄冷系统构成的联合供冷系统的递 阶优化方法中通过建立了电制冷机和蓄冰槽的经济模型，并分别使用混合整数规划和分布 式预测控制解决稳态优化问题和动态优化问题，定量地确定联合供冷系统的调度方法。
[0150] 以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述 特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变形或修改，这并不影 响本发明的实质内容。
【主权项】
1. 一种联合供冷系统的递阶优化方法，所述联合供冷系统是由常规电制冷机和冰蓄冷 系统构成的联合供冷系统，其特征在于，所述方法包括如下步骤： 步骤Sl :利用电制冷机运行数据拟合其耗电功率关于制冷功率和冷却水进口温度的 函数； 步骤S2 :利用上一步骤中得到的电制冷机耗电功率关于其制冷功率和冷却水进口温 度的函数建立电制冷机供冷的经济模型并进而建立冰蓄冷系统供冷的经济模型，最终形成 稳态经济优化问题的目标函数和约束条件并进行求解，得到各电制冷机和蓄冰槽在经济效 益最优条件下的启停状态和供冷功率； 步骤S3 :在上一步骤的基础上，设计目标耦合的协调分布式预测控制方法重新优化整 个供冷系统在每一采样时刻各电制冷机的实时供冷功率设定值以提高整个系统响应负荷 的动态性能； 步骤S4 :利用迭代方法求解上述分布式预测控制问题，得到各 电制冷机每一采样时刻 的供冷功率设定值。2. 如权利要求1所述的联合供冷系统的递阶优化方法，其特征在于，所述的步骤Sl中： 离心式电制冷机最重要的性能参数能效比为制冷功率与耗电功率的比率，根据实际的冷却 水进口温度、耗电功率和制冷功率数据通过回归分析将耗电功率表示成为冷却水进口温度 和制冷功率二元二次函数形式，以便利用数学规划工具进行后续步骤优化问题的求解。3. 如权利要求2所述的联合供冷系统的递阶优化方法，其特征在于，所述步骤S2中： 每台电制冷机有开启和关闭两个工作状态，采用0-1离散变量Y来描述电制冷机的工作状 态；对于双工况电制冷机空调工况和制冰工况的区分则使用两个0-1离散变量来分别描 述；电制冷机的制冷功率则用连续变量来描述；根据步骤Sl中进行的回归分析，对于常规 电制冷机，其耗电功率Pin表示成关于制冷功率Pwt与冷却水进口温度Tf的二元二次函数： Pin=S(PmtXn) ⑴ 于是电制冷机供冷产生的电费花销C为： C^YxgiPm,,J：")xpxt (2) 其中P为供冷时段电价，t为供冷时段时间长度，Y为电制冷机启停标志，If)为 步骤si中拟合出的电制冷机耗电功率关于制冷功率和冷却水进口温度的函数； 在忽略蓄冰槽储冷量随储存时间损耗的前提下，为了考量蓄冰槽供冷的花销，引入蓄 冰槽冷量均价Ps概念，描述单位冷量的价格，这个价格与充冷量、充冷时段电价和双工况电 制冷机制冰工况输入输出关系相关，由下式定义：(3) 其中Pinit为冷量均价初值，W init为蓄冰槽冷量初值，tin为冲冷时间，t °ut为放冷时间，P m为双工况电制冷机总制冰功率，Ps为蓄冰槽供冷功率，Y s为蓄冰槽供冷启停标志；由蓄冰槽 冷量均价的定义看出，在不考虑储冷量随时间损耗的前提下，均价与蓄冰槽是否放冷无关， 只在蓄冰槽冲冷时可能改变； 蓄冰槽储冷量增量ΛW为： Δ W = tinP-t°utPsYs (4) 目标函数为各种电制冷机空调工况供冷与蓄冰槽供冷所需电费之和：(5) P为供冷时段电价，t为供冷时段时间长度，Yi为第i台电制冷剂启停标志， 为第i台电制冷机耗电功率关于制冷功率和冷却水进口温度的函数，其中i为第i台电 制冷机供冷功率，T。，，为第i台电制冷机冷却水进口温度，Y s为蓄冰槽供冷启停标志，p 3为 蓄冰槽冷量均价，Ps为蓄冰槽供冷功率，Y s为蓄冰槽供冷启停标志； 约束条件包括供冷功率与预测负荷平衡、各电制冷机额定功率限制、蓄冰槽储冷量上 下限和冰蓄冷系统工作模式唯一：(6) 其中Yi为第i台电制冷剂启停标志，P 为第i台电制冷机供冷功率，P s为蓄冰槽供 冷功率，Ys为蓄冰槽供冷启停标志，Pp为预测负荷，a (〇, 1)为第i台电制冷机最小输 出功率比例，Pi rated为第i台电制冷机额定输出功率，W为蓄冰槽储冷量，Wmax和Wmin为蓄冰 槽储冷量上下限；Yic^为双工况制冷机制冰工况启停标志，默认制冰工况时所有双工况冷机 同时开启，加快制冰速度；YaiM为第i台双工况制冷机空调工况启停标志，共N台双工况制 冷机； 最后，在形成稳态经济优化问题的目标函数和约束条件的前提下对上述优化问题进行 求解，得到各电制冷机和蓄冰槽在经济效益最优条件下的启停状态和供冷功率。4.如权利要求1所述的联合供冷系统的递阶优化方法，其特征在于，步骤S3中，所述标 耦合的协调分布式预测控制方法，具体如下： 分布式控制系统按照空调工况电制冷机数量分为N个子系统，每个子系统由模型预测 控制（MPC)和相应控制的电制冷机组成，不同子系统MPC的数据使用通信器进行实时交换 完成计算任务，从而使得各台空调工况电制冷机供冷功率跟踪步骤S2中得到的经济效益 最优条件下的供冷功率并且总供冷功率跟踪各台空调工况电制冷机经济效益最优条件下 的供冷功率之和； 电制冷机制冷功率动态过程模型传递函数如下：C7) τ为一阶惯性环节时间常数，Td为延时时间，则第s个子系统离散状态空间模型形式 为：(8) S = 1，2，L，N 其中us(k)为k时刻电制冷机输出功率设定值，即为步骤S3中优化问题的优化变量， PsGO为k+Ι时刻电制冷机实际输出功率，xs(k) = [ps(k)ps(k+l)L ps(k+nd)]T为状态变量， nd为电制冷机延时时间相对于离散系统采样时间的倍数;A s、Bs、Cs为状态空间方程中的系 数矩阵，其中As e R(M"+1)X(M"+1)，Bi e R(""+1|xl由离散系统采样时间和电制冷机动态过程时间 常数决定，R表示实数集，设采样时间为Λ t，则有：Bs= [0 L O 1-θ_Δ?/τ]τ (9) Cs= [I 0 L 0] 对于第s个子系统，s = 1，2, ···，N ;MPC协调策略包括电制冷机在未来P个采样时 刻组成的预测时域内预测输出功率尽可能接近步骤S2得到的单台空调工况电制冷机经济 效益最优条件下的供冷功率rs，并且为了控制所有电制冷机预测输出功率跟上各台空调工 况电制冷机经济效益最优条件下的供冷功率之和r，在目标函数中引入预测时域内各台空 调工况电制冷机总预测功率与经济效益最优条件下的供冷功率之和差距的平方和，因此不 同子系统MPC目标函数出现耦合；综上，第s个子系统MPC目标函数如下所示：其中，Us= [u s (k) Us (k+Ι)…us (k+P-1) ]τ为预测时域内电制冷机输出功率设定值序列， us(k)为k时刻电制冷机输出功率设定值，即为优化问题的优化变量；而Q, R e [〇, 1]为 权重系数，分别用来衡量单台电制冷机预测负荷跟踪程度与所有电制冷机总负荷跟踪程度 的重要性；rs(k+i)为（k+i)时刻第s台空调工况电制冷机经济效益最优条件下的供冷功 率，r(k+i)为（k+i)时刻各台空调工况电制冷机在经济效益最优条件下的供冷功率之和； + 〗㈨为第s台空调工况电制冷机在k时刻通过预测模型预测的其在（k+i)时刻的输 出功率，为除第s台空调工况电制冷机之外的其余空调工况电制冷机在k时 刻通过预测模型预测的其在（k+i)时刻的输出功率之和； MPC预测模型为：i = 1，2，L，P 象认+?μ)为在k时刻通过预测模型计算出的（k+i)时刻第8台空调工况电制冷机状 态空间方程中状态变量，/\认+ /_1幻为第S台空调工况电制冷机在k时刻通过预测模型预 测的其在（k+i)时刻的输出功率;AsH和A/为状态空间方程中系数矩阵的相应次幂； 对于第s个子系统MPC控制器，整个优化问题形式为：i = 1,2, L, P 其中，{us，min，us，max}，{p s，min，ps，max}为控制变量的边界约束。5.如权利要求1-4任一项所述的联合供冷系统的递阶优化方法，其特征在于，所述步 骤S4中，采取迭代方法求解步骤S3中的分布式预测控制问题，具体步骤如下： 步骤S41 :在k时刻，根据步骤S3中得到的每台空调工况电制冷机经济效益最优条件 下的供冷功率rs，初始化每个子系统在预测时域内的电制冷机输出功率设定值序列Us，令 迭代次数1 = 〇 :上述〃 !(A + 1)为第1次迭代在（k+Ι)时刻的电制冷机输出功率设定值，P为预测时域长 度； 步骤S42:根据MPC预测模型，计算出每个子系统MPC控制时域内输出量预测值 β、= [/(U +11幻八认+ 2 I幻…/)你+ PI幻]τ，并通知给其他N-I个子系统MPC ; 步骤S43 :对于每个子系统，求解步骤S3中的优化问题，得到最优解，即最优制冷机输 出功率设定值序列Uf; 步骤S44 :检查所有子系统MPC的收敛条件是否满足，即对给定的精度ε se R(s = 1，2, L，N)，是否均有|υΓ -K|| < & Q = 1 2,L , A〇，如果所有的子系统收敛条件满足，令每 个子系统最优控制量< 与最优制冷机输出功率设定值序列Uf相等，即U! 转步骤 S45;否则，4α为介于〇和1之间的常数，I = 1+1， 转步骤S42 ; 步骤S45:在k时刻选取控制量<=[1 O L 0]U〔（λ' = 1，2,···，Λ〇作用于相应的子系 统； 步骤S46 :滚动移位到下一个采样时刻，即k+Ι - k，返回步骤S41，重复上述过程。
【专利摘要】本发明提供了一种针对由常规电制冷机和冰蓄冷系统构成的联合供冷系统的递阶优化方法，首先根据各冷源经济模型和分时电价政策，采用混合整数规划方法，优化各冷源的开启状态及最优功率设定值，同时考虑到各冷源的动态调节特性以及物理约束，为了更好地提高系统动态性能，设计目标耦合的协调分布式预测控制方法，在分布式框架下重新优化各冷源的设定值，使得各冷源在动态过程中保证总负荷的同时尽可能跟踪最优制冷功率设定值。本发明分别采用混合整数规划和分布式预测控制方法定量解决联合供冷系统的稳态和动态调度问题，为供冷系统能源管理与优化调度提供合理的建议。
【IPC分类】G05B13/04
【公开号】CN104898422
【申请号】CN201510225593
【发明人】蔡旭, 刘楚晖, 郑毅
【申请人】上海交通大学
【公开日】2015年9月9日
【申请日】2015年5月5日
【技术领域】
[0001] 本发明涉及能源网调度领域，特别是一种由常规电制冷机和冰蓄冷系统构成的联 合供冷系统的递阶优化方法。
【背景技术】
[0002] 由于大型城市普遍存在用地紧张的问题，为了利用有限的土地面积尽可能增大城 市容量，高层建筑成为城市设计者们理想的选择。高层建筑往往具有独立的供冷网络，提供 整栋楼各种各样的用冷需求。随着各种新兴技术的出现以及分时电价政策的推广，仅仅使 用常规电制冷机组的供冷系统经济效益不能令人满意，为了提高高层建筑供冷系统的经济 效益，多供冷机组互补协同系统应用得越来越普遍。典型的多供冷机组互补协同系统一般 包括：冰蓄冷、常规电制冷机，冷热电三联供及地源热泵等，其中常规电制冷机和冰蓄冷一 般共同承担90%以上的冷负荷。由于能够自主控制储能与用能时段，冰蓄冷系统越来越受 到人们的青睐，有效地对其与常规电制冷机进行协调控制是确保该协同系统节能降耗并提 尚经济效益的关键。
[0003] 目前冰蓄冷系统与常规电制冷机组的联合控制多停留在定性层面，即在电价谷时 制冰蓄冷，在电价峰时融冰放冷，缺乏对整个系统定量的分析与控制，这样导致了其控制方 法并不优越，而且由于电制冷机从启动达到设定的供冷功率有一个动态过程，所以仅仅考 虑稳态的控制方法难以确保整个系统的实际供冷功率迅速跟上预测负荷。因此，需要研宄 一种优化控制方法使得系统在动态过程中保证总负荷的前提下高效经济运行。
[0004] 对于高层建筑供冷系统这样一个分布式系统，通常可采用一般稳态优化方法对其 进行调度，然而由于电制冷机动态响应问题，进一步采用分布式预测控制优化电制冷机供 冷功率动态性能。对于分布式系统的动态控制，目前已经由集中式向分布式控制结构转化， 因为集中式控制会导致方法复杂且计算量偏大。在实际工业应用中，模型预测控制（MPC) 可以通过被控对象的脉冲或阶跃响应序列预测对象未来变化并且可自然引入过程时滞，对 于多变量系统可以获得比常规控制方法更加优秀的控制效果。在大型系统中，传统的模型 预测控制衍生出分布式预测控制（DMPC)这种形式，将大型系统分解成多个子系统，子系统 与子系统之间数据相互交流，每个子系统利用自身以及其余子系统的数据来完成自身的 MPC问题的求解，显著降低了计算量并可以达到整个系统的最优状态。
[0005] 因此，本发明拟借鉴以上对耦合系统的协调方法，研宄目标耦合但子系统间不耦 合的分布式系统的动态优化，使得整个系统在确保提供总预测冷负荷的前提下高效经济运 行。
【发明内容】
[0006] 针对现有技术存在的上述不足，本发明提出了一种联合供冷系统的递阶优化方 法，该方法对电制冷机和蓄冰槽进行经济性建模并且运用混合整数规划的方法得到经济指 标最优条件下各电制冷机与蓄冰槽的稳态运行时的供冷功率；下层运用分布式预测控制方 法保证各电制冷机和蓄冰槽在动态过程中尽可能跟踪上层经济优化得出的稳态设定值并 且总供冷功率跟踪预测负荷，解决目标耦合但子系统不耦合的分布式预测控制问题。
[0007] 为实现上述的目的，本发明提供一种联合供冷系统的递阶优化方法，所述联合供 冷系统是由常规电制冷机和冰蓄冷系统构成的联合供冷系统，所述方法具体包括如下步 骤：
[0008] 步骤S1 :利用电制冷机运行数据拟合其耗电功率关于制冷功率和冷却水进口温 度的函数；
[0009] 步骤S2 :利用上一步骤中得到的电制冷机耗电功率关于其制冷功率和冷却水进 口温度的函数建立电制冷机供冷的经济模型并进而建立冰蓄冷系统供冷的经济模型，最终 形成稳态经济优化问题的目标函数和约束条件并进行求解，得到各电制冷机和蓄冰槽在经 济效益最优条件下的启停状态和供冷功率；
[0010] 步骤S3 :在上一步骤的基础上，设计目标耦合的协调分布式预测控制方法重新优 化整个供冷系统在每一采样时刻各电制冷机的实时供冷功率设定值以提高整个系统响应 负荷的动态性能；
[0011] 步骤S4 :利用迭代方法求解上述分布式预测控制问题，得到各电制冷机每一采样 时刻的供冷功率设定值。
[0012] 较佳地，所述的步骤S1中：离心式电制冷机最重要的性能参数能效比为制冷功率 与耗电功率的比率。工程实践中，需要根据实际的冷却水进口温度、耗电功率和制冷功率等 数据通过回归分析将耗电功率表示成为冷却水进口温度和制冷功率二元二次函数形式以 便利用数学规划工具进行后续步骤优化问题的求解。
[0013] 较佳地，所述步骤S2中：每台电制冷机有开启和关闭两个工作状态，故采用0-1离 散变量Y来描述电制冷机的工作状态。对于双工况电制冷机空调工况和制冰工况的区分则 需使用两个0-1离散变量来分别描述。电制冷机的制冷功率则用连续变量来描述。根据步 骤S1中进行的回归分析，对于常规电制冷机，其耗电功率可以表示成关于制冷功率与冷却 水进口温度的二元二次函数：
[0015] 于是电制冷机供冷产生的电费花销C为：
[0016] C = YXg(P0Ut,T) XpXt (2)
[0017] 其中p为供冷时段电价，t为供冷时段时间长度。Y为电制冷机启停标志， 为步骤S1中拟合出的电制冷机耗电功率关于制冷功率和冷却水进口温度的函 数。
[0018] 在忽略蓄冰槽储冷量随储存时间损耗的前提下，为了考量蓄冰槽供冷的花销，弓丨 入蓄冰槽冷量均价Ps概念，描述单位冷量的价格，这个价格与充冷量、充冷时段电价和双工 况电制冷机制冰工况输入输出关系相关，由下式定义：
[0020] 其中pinit为冷量均价初值，w init为蓄冰槽冷量初值，t in为冲冷时间，t °ut为放冷时 间，Pm为双工况电制冷机总制冰功率，P s为蓄冰槽供冷功率，Y 3为蓄冰槽供冷启停标志。由 蓄冰槽冷量均价的定义可以看出，在不考虑储冷量随时间损耗的前提下，均价与蓄冰槽是 否放冷无关，只在蓄冰槽冲冷时可能改变。
[0021] 蓄冰槽储冷量增量为：
[0022] A ff = tinP-t0UtPsYs (4)
[0023]目标函数为各种电制冷机空调工况供冷与蓄冰槽供冷所需电费之和：
[0025] p为供冷时段电价，t为供冷时段时间长度，Yi为第i台电制冷剂启停标志， 为第i台电制冷机耗电功率关于制冷功率和冷却水进口温度的函数，其中 为第i台电制冷机供冷功率，TMin为第i台电制冷机冷却水进口温度，蓄冰槽供冷启 停标志，ps为蓄冰槽冷量均价，P s为蓄冰槽供冷功率，Y s为蓄冰槽供冷启停标志。
[0026]约束条件包括供冷功率与预测负荷平衡、各电制冷机额定功率限制、蓄冰槽储冷 量上下限和冰蓄冷系统工作模式唯一：
[0028] 其中I为第i台电制冷剂启停标志，P。^为第i台电制冷机供冷功率，P S为蓄冰 槽供冷功率，Ys为蓄冰槽供冷启停标志，Pp为预测负荷，a iE (0, 1)为第i台电制冷机最 小输出功率比例，Pi Mted为第i台电制冷机额定输出功率，W为蓄冰槽储冷量，Wmax和Wmin为 蓄冰槽储冷量上下限；Y ira为双工况制冷机制冰工况启停标志，默认制冰工况时所有双工况 冷机同时开启，加快制冰速度；YaiM为第i台双工况制冷机空调工况启停标志，共N台双工 况制冷机；
[0029] 最后，在形成稳态经济优化问题的目标函数和约束条件的前提下对上述优化问题 进行求解，得到各电制冷机和蓄冰槽在经济效益最优条件下的启停状态和供冷功率。
[0030] 较佳地，步骤S3中，所述标耦合的协调分布式预测控制方法，具体如下：
[0031] 分布式控制系统按照空调工况电制冷机数量分为N个子系统，每个子系统由模型 预测控制（MPC)和相应控制的电制冷机组成，不同子系统MPC的数据使用通信器进行实时 交换完成计算任务，从而使得各台空调工况电制冷机供冷功率跟踪步骤S2中得到的经济 效益最优条件下的供冷功率并且总供冷功率跟踪各台空调工况电制冷机经济效益最优条 件下的供冷功率之和；
[0032] 电制冷机制冷功率动态过程模型传递函数如下：
[0034]t为一阶惯性环节时间常数，td为延时时间，则第s个子系统离散状态空间模型 形式为：
[0036] 其中us(k)为k时刻电制冷机输出功率设定值，即为步骤S3中优化问题的优化变 量，Ps(k)为k+1时刻电制冷机实际输出功率，xs(k) = [ps(k)ps(k+l)L ps(k+nd)]T为状态 变量，nd为电制冷机延时时间相对于离散系统采样时间的倍数;A s、Bs、Cs为状态空间方程中 的系数矩阵，其中
由离散系统采样时间和电制冷机动态过 程时间常数决定，R表示实数集，设采样时间为At，则有：
[0038] Bs= [0 L 0 l-e_At/T]T (9)
[0039] Cs= [1 0 L 0]
[0040] 对于第s个子系统，s = 1，2,…，N ;MPC协调策略包括电制冷机在未来P个采样时 刻组成的预测时域内预测输出功率尽可能接近步骤S2得到的单台空调工况电制冷机经 济 效益最优条件下的供冷功率r s，并且区别于一般的MPC策略，为了控制所有电制冷机预测输 出功率尽可能跟上各台空调工况电制冷机经济效益最优条件下的供冷功率之和r，在目标 函数中引入预测时域内各台空调工况电制冷机总预测功率与经济效益最优条件下的供冷 功率之和差距的平方和，因此不同子系统MPC目标函数出现耦合；综上，第s个子系统MPC 目标函数如下所示：
[0042] 其中，Us= [1130011；30^+1)~11；30^-1)]1为预测时域内电制冷机输出功率设定值 序列，u s(k)为k时刻电制冷机输出功率设定值，即为优化问题的优化变量；而Q，R G [0, 1] 为权重系数，分别用来衡量单台电制冷机预测负荷跟踪程度与所有电制冷机总负荷跟踪程 度的重要性；rs(k+i)为（k+i)时刻第s台空调工况电制冷机经济效益最优条件下的供冷功 率，r(k+i)为（k+i)时刻各台空调工况电制冷机在经济效益最优条件下的供冷功率之和； /\0 + /|幻为第s台空调工况电制冷机在k时刻通过预测模型预测的其在（k+i)时刻的输 出功率，
为除第s台空调工况电制冷机之外的其余空调工况电制冷机在k时 刻通过预测模型预测的其在（k+i)时刻的输出功率之和；
[0043] MPC预测模型为：
[0045] is(々+ M幻为在k时刻通过预测模型计算出的（k+i)时刻第S台空调工况电制冷 机状态空间方程中状态变量，+ ZI幻为第s台空调工况电制冷机在k时刻通过预测模 型预测的其在（k+i)时刻的输出功率;A/4和A /为状态空间方程中系数矩阵的相应次幂。
[0046] 对于第s个子系统MPC控制器，整个优化问题形式为：
[0049]其中，{us，min，us，max}，{p s，min，ps，J为控制变量的边界约束。
[0050] 较佳地，所述步骤S4具体包括：
[0051] 步骤S41 :在k时刻，根据步骤S3中得到的每台空调工况电制冷机经济效益最优 条件下的供冷功率rs，初始化每个子系统在预测时域内的电制冷机输出功率设定值序列 U s，令迭代次数1 = 0 :
[0053] 上述<〇 + 1)为第1次迭代在（k+1)时刻的电制冷机输出功率设定值，P为预测时 域长度；
[0054] 步骤S42:根据MPC预测模型，计算出每个子系统MPC控制时域内输出量预测值
，并通知给其他N-1个子系统MPC ;
[0055] 步骤S43:对于每个子系统，求解步骤S3中的优化问题，得到最优解，即最优制冷 机输出功率设定值序列Uf 1;
[0056] 步骤S44 :检查所有子系统MPC的收敛条件是否满足，即对给定的精度e sG R(s =1，2，L，N)，是否均有
如果所有的子系统收敛条件满足，令 每个子系统最优控制量a与最优制冷机输出功率设定值序列Uf"相等，即K = uf，转步 骤S45;否则，令
a为介于〇和1之间的常数，引入该 参数为了尽量避免方法不收敛的情况，按需求选取数值大小，1 = 1+1，转步骤S42 ;
[0057] 步骤S45:在k时刻选取控制量<=[丨0 L0]U:(.s'=l二…，A〇作用于相应的 子系统。
[0058] 步骤S46:滚动移位到下一个采样时刻，即k+1 - k，返回步骤S41，重复上述过程。
[0059] 与现有技术相比，本发明具有以下有益效果：
[0060] 本发明针对由常规电制冷机和冰蓄冷系统构成的联合供冷系统的递阶优化方法 中通过建立了电制冷机和蓄冰槽的经济模型，并分别使用混合整数规划和分布式预测控制 解决稳态优化问题和动态优化问题，定量地确定联合供冷系统的调度方法，使得整个系统 在确保提供总预测冷负荷的前提下高效经济运行，为供冷系统能源管理与优化调度提供合 理的建议。
【附图说明】
[0061] 通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、 目的和优点将会变得更明显：
[0062]图1为由常规电制冷机和冰蓄冷系统构成的联合供冷系统结构图；
[0063] 图2为本发明的一实施例的流程图；
[0064] 图3为本发明的一实施例的方法示意图；
[0065] 图4为本发明的一实施例的分布式预测控制迭代求解方法流程图；
[0066] 图5为本发明的一实施例中3900kW冷机回归分析图；
[0067] 图6为本发明的一实施例中2150kW冷机回归分析图；
[0068] 图7为本发明的一实施例中6329kW冷机回归分析图；
[0069] 图8为本发明的一实施例中负荷预测图；
[0070] 图9为本发明的一实施例中两种策略每半小时的电费对比图；
[0071] 图10为本发明的一实施例的动态性能优化仿真结果，其中：(a)为总体，（b)为第 一台3900kW电制冷机，（c)为第二台3900kW电制冷机，（d)为第一台2150kW电制冷机。
【具体实施方式】
[0072] 下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术 人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术 人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进。这些都属于本发明 的保护范围。下面结合附图来对本发明的由常规电制冷机和冰蓄冷系统构成的联合供冷系 统递阶优化方法作进一步详细的说明。
[0073] 如图1所示，为由常规电制冷机和冰蓄冷系统构成的联合供冷系统结构。
[0074] 图2为本发明的一个【具体实施方式】的由常规电制冷机和冰蓄冷系统构成的联合 供冷系统递阶优化方法的流程图。首先根据各冷源经济模型和分时电价政策，采用混合整 数规划方法，优化各冷源的开启状态及最优功率设定值，同时考虑到各冷源的动态调节特 性以及物理约束，为了更好地提高系统动态性能，设计了一种目标耦合的协调分布式预测 控制方法，在分布式框架下重新优化各冷源的设定值，使得各冷源在动态过程中保证总负 荷的同时尽可能跟踪最优制冷功率设定值。
[0075] 具体的，如图2所示，本发明的由常规电制冷机和冰蓄冷系统构成的联合供冷系 统递阶优化方法，包括以下的步骤：
[0076] 步骤S1 :利用电制冷机运行数据拟合其耗电功率关于制冷功率和冷却水进口温 度的函数；
[0077] 离心式电制冷机最重要的性能参数能效比（C0P)为制冷功率与耗电功率的比率。 工程实践中，需要根据实际的冷却水进口温度、耗电功率和制冷功率等数据通过回归分析 将耗电功率表示成为冷却水进口温度和制冷功率二元二次函数形式以便利用数学规划工 具进行后续步骤优化问题的求解。
[0078] 步骤S2 :建立电制冷机和冰蓄冷系统供冷的经济模型，形成稳态经济优化问题的 目标函数和约束条件并进行求解；
[0079] 每台电制冷机有开启和关闭两个工作状态，故采用0-1离散变量Y来描述电制冷 机的工作状态。对于双工况电制冷机空调工况和制冰工况的区分则需使用两个0-1离散变 量来分别描述。电制冷机的制冷功率则用连续变量来描述。根据步骤S1中进行的回归分 析，对于常规电制冷机，其耗电功率可以表示成关于制冷功率与冷却水进口温度的二元二 次函数：
[0081] 于是电制冷机供冷产生的电费花销C为：
[0082] C = YXg(P0Ut,T) XpXt (2)
[0083] 其中p为供冷时段电价，t为供冷时段时间长度。
[0084] 在忽略蓄冰槽储冷量随储存时间损耗的前提下，为了考量蓄冰槽供冷的花销，弓丨 入蓄冰槽冷量均价Ps概念，描述单位冷量的价格，这个价格与充冷量、充冷时段电价和双工 况电制冷机制冰工况输入输出关系相关，由下式定义：
[0086] 其中pinit为冷量均价初值，W init为蓄冰槽冷量初值，t in为冲冷时间，t °ut为放冷时 间，Pm为双工况电制冷机总制冰功率，P s为蓄冰槽供冷功率，Y 3为蓄冰槽供冷启停标志。由 蓄冰槽冷量均价的定义可以看出，在不考虑储冷量随时间损耗的前提下，均价与蓄冰槽是 否放冷无关，只在蓄冰槽冲冷时可能改变。
[0087] 蓄冰槽储冷量增量为：
[0088] A ff = tinP-t0UtPsYs (4)
[0089] 目标函数为各种电制冷机空调工况供冷与蓄冰槽供冷所需电费之和：
[0091] P。^为第i台电制冷机供冷功率，T Min为第i台电制冷机冷却水进口温度。
[0092] 主要的约束条件包括供冷功率与预测负荷平衡、各电制冷机额定功率限制、蓄冰 槽储冷量上下限和冰蓄冷系统工作模式唯一：
[0094] 其中Pp为预测负荷，a iG (〇, 1)为第i台电制冷机最小输出功率比例，P i Mtral为 第i台电制冷机额定输出功率，Wmax和Wmin为蓄冰槽储冷量上下限。Y 为双工况制冷机制 冰工况启停标志，默认制冰工况时所有双工况冷机同时 开启，加快制冰速度。YaiM为第i台 双工况制冷机空调工况启停标志，共N台双工况制冷机。
[0095] 步骤S3:在上一步稳态经济优化得到的各电制冷机最优制冷功率的基础上设计 一种目标耦合的协调分布式预测控制方法优化整个供冷系统响应负荷的动态性能；
[0096] 分布式控制系统按照空调工况电制冷机数量分为N个子系统，每个子系统由MPC 和相应控制的电制冷机组成，不同子系统MPC的数据使用通信器进行实时交换完成计算任 务从而使得各台空调工况电制冷机供冷功率跟踪步骤S2中计算出的期望值并且总供冷功 率跟踪总期望值，具体的分布式控制系统结构如图3所示。
[0097] 电制冷机制冷功率动态过程模型传递函数如下：
[0099] t为一阶惯性环节时间常数，td为延时时间，则第s个子系统离散状态空间模型 形式为：
[0101] 其中U(k)为k时刻电制冷机输出功率设定值，即为优化问题的优化变量，ps(k)为 k+1时刻电制冷机实际输出功率，nd为电制冷机延时时间相对于离散系统采样时间的倍数。 xs(k) = [ps(k)ps(k+l)L ps(k+nd)]T为状态变量。
由采样时间 和电制冷机时间常数决定，设采样时间为At，则有：
[0103] Bs= [0 L 0 l-e_At/T]T (9)
[0104] Cs= [1 0 L 0]
[0105] 对于第s个子系统（s = 1，2,…，N)，MPC协调策略包括电制冷机在未来P个采样 时刻组成的预测时域内预测输出功率尽可能接近步骤S2计算出的期望功率r s，并且区别于 一般的MPC策略，为了控制所有电制冷机预测输出功率尽可能跟上总期望功率r，本策略在 目标函数中引入预测时域内总预测功率与总期望功率差距的平方和，因此不同子系统MPC 目标函数出现耦合。综上，第s个子系统MPC目标函数如下所示：
[0107] 其中，Us= [u(k)u(k+l)-u(k+p-l)] T为预测时域内电制冷机输出功率设定值序 列，而Q，R e [0, 1]为权重系数，分别用来衡量单台电制冷机预测负荷跟踪程度与所有电制 冷机总负荷跟踪程度的重要性。
[0108] MPC预测模型为：
[0110]对于第s个子系统MPC控制器，整个优化问题形式为：
[0113]其中，{us，min，us，max}，{p s，min，ps，J为控制变量的边界约束。
[0114] 步骤S4 :利用迭代方法求解上述分布式预测控制问题，得到各电制冷机每一采样 时刻的供冷功率设定值；
[0115] 迭代方法流程图如图4所示。
[0116] 步骤S41 :在k时刻，根据步骤S3中求解出的每台电制冷机理想制冷功率rs初始 化每个子系统在预测时域内的U s，令迭代次数1 = 0，
[0118]步骤S42 :根据MPC预测模型，计算出每个子系统MPC控制时域内输出量预测值
，并通知给其他N-1个子系统MPC。
[0119] 步骤S43:每个子系统步骤S3中的优化问题，得到最优解lUi+1。
[0120] 步骤S44:检查所有子系统MPC的收敛条件是否满足，即对给定的精 度e sG R(s = 1，2，L，N)，是否均有
，如果所有的 子系统收敛条件满足，令每个子系统最优控制量M=U(/+I，转步骤S45;否则，令
，a为介于〇和1之间的常数，引入该参数为了尽量避 免方法不收敛的情况，按需求选取数值大小，1 = 1+1，转步骤S42。
[0121] 步骤S45:在k时刻选取控制量》:=[1 0 L 0]U〗（5 = 1,2,作用于相应的 子系统。
[0122] 步骤S46 :滚动移位到下一个采样时刻，即k+1 - k，返回步骤S41，重复上述过程。
[0123] 基于上述的实施过程，本发明以上海某超高层建筑低区供冷系统为例进行了仿真 验证。在该供冷系统中，常规电制冷机共三台，两台额定制冷功率3900kW，一台额定制冷功 率2150kW，双工况电制冷机共三台，空调工况额定制冷功率6329kW，制冰工况额定制冷功 率 3868kW。
[0124] 首先，按步骤S1根据电制冷机运行实测数据进行输入输出关系函数的拟合。
[0125] 图5,图6,图7分别为3900kW、2150kW和6329kW离心式电制冷机回归分析结果图。 各图中XI为冷机供冷功率标么值，X2为冷机冷却水进口温度，Y1为冷机耗电功率标幺值。 图中红色曲线包围的区域为置信水平99 %所对应的置信区间
[0126] 从图5-7图中可以看出，二元二次回归分析可以相当精确地建立电制冷机输入电 功率与输出冷功率和冷却水进口温度的函数关系，具体的函数表达式如下列各式，P in为电 制冷机耗电功率，电制冷机制冷功率，2：?为电制冷机冷却水进口温度。
[0127] 对于3900kW离心式电制冷机有：
[0129] 对于2150kW离心式电制冷机有：
[0131] 对于双工况电制冷机6329kW空调工况有：
[0133] 对于双工况电制冷机的制冰工况，为了加快制冰速度，在分配制冰的情况下电制 冷机全功率制冰，所以不存在部分功率制冰的工况，此时制冰功率cold和耗电功率elec分 别仅为关于冷却水进口温度凡^的一元二次函数：
[0135] 然后，按步骤S2进行稳态经济优化仿真，预测负荷和电价政策如下。如图8所示。
[0136] 表1分时电价政策
[0139] 为了经济优化策略节费效果，选取一种常规的稳态调度策略作为对比，该策略在 电价谷值时优先蓄冰槽冲冷，在电价平值与峰值时优先电制冷机供冷，各电制冷机按额定 功率由小到大的顺序启用。两种策略每半小时电费对比如图9。
[0140] 经过计算，稳态经济优化策略在给定的预测负荷条件下，6小时内总电费花销为 3970元，而对比策略总电费花销为4053元，6小时内本发明策略相对节费2%。
[0141] 随后，按照步骤S3和S4进行分布式预测控制优化策略设计和问题求解。动态优 化采样间隔为2. 5分钟，即步骤S2的稳态调度策略每30分钟计算一次各电制冷机与蓄冰 槽最优供冷功率后，而动态性能优化每2. 5分钟更新一次各电制冷机供冷功率设定值以确 保各电制冷机和蓄冰槽在动态过程中尽可能跟踪步骤S2经济优化得出的稳态设定值并且 总供冷功率跟踪预测负荷。冰槽供冷一阶模型时间常数设为2分钟，而各电制冷机空调工 况动态模型时间常数和延时时间设定如下表所示：
[0142] 表2空调工况电制冷机动态参数
[0144] 总体的和各台启用过的空调工况电制冷机的期望承担负荷与分别应用与不应用 DMPC策略导致的实时供冷功率曲线如组图10所示，由于在当前预测负荷情况下6329kW电 制冷机并未开启，所以组图中只包括3900kW和2150kW电制冷机。
[0145] 为了定量描述采用或不采用DMPC策略时整个系统以及每台启用过的空调工况电 制冷机实时供冷功率跟步骤S2中计算出的经济最优条件下的最优供冷功率接近程度的区 另IJ，计算了 6个小时内每个采样时刻偏差的平方和，数值如表3所示，可以看出采用动态性 能优化策略后偏差平方和明显减小，实现了总体系统和每个子系统的迅速响应，体现出良 好的控制效果。
[0146] 表3动态优化效果
[0148] 综合仿真结果，本发明提出的由常规电制冷机和冰蓄冷系统构成的联合供冷系统 的递阶优化方法兼具较好的经济价值和动态控制效果，对于实际供冷系统经济优化运行有 一定的指导意义。
[0149] 综上所述，本发明的针对由常规电制冷机和冰蓄冷系统构成的联合供冷系统的递 阶优化方法中通过建立了电制冷机和蓄冰槽的经济模型，并分别使用混合整数规划和分布 式预测控制解决稳态优化问题和动态优化问题，定量地确定联合供冷系统的调度方法。
[0150] 以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述 特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变形或修改，这并不影 响本发明的实质内容。
【主权项】
1. 一种联合供冷系统的递阶优化方法，所述联合供冷系统是由常规电制冷机和冰蓄冷 系统构成的联合供冷系统，其特征在于，所述方法包括如下步骤： 步骤Sl :利用电制冷机运行数据拟合其耗电功率关于制冷功率和冷却水进口温度的 函数； 步骤S2 :利用上一步骤中得到的电制冷机耗电功率关于其制冷功率和冷却水进口温 度的函数建立电制冷机供冷的经济模型并进而建立冰蓄冷系统供冷的经济模型，最终形成 稳态经济优化问题的目标函数和约束条件并进行求解，得到各电制冷机和蓄冰槽在经济效 益最优条件下的启停状态和供冷功率； 步骤S3 :在上一步骤的基础上，设计目标耦合的协调分布式预测控制方法重新优化整 个供冷系统在每一采样时刻各电制冷机的实时供冷功率设定值以提高整个系统响应负荷 的动态性能； 步骤S4 :利用迭代方法求解上述分布式预测控制问题，得到各 电制冷机每一采样时刻 的供冷功率设定值。2. 如权利要求1所述的联合供冷系统的递阶优化方法，其特征在于，所述的步骤Sl中： 离心式电制冷机最重要的性能参数能效比为制冷功率与耗电功率的比率，根据实际的冷却 水进口温度、耗电功率和制冷功率数据通过回归分析将耗电功率表示成为冷却水进口温度 和制冷功率二元二次函数形式，以便利用数学规划工具进行后续步骤优化问题的求解。3. 如权利要求2所述的联合供冷系统的递阶优化方法，其特征在于，所述步骤S2中： 每台电制冷机有开启和关闭两个工作状态，采用0-1离散变量Y来描述电制冷机的工作状 态；对于双工况电制冷机空调工况和制冰工况的区分则使用两个0-1离散变量来分别描 述；电制冷机的制冷功率则用连续变量来描述；根据步骤Sl中进行的回归分析，对于常规 电制冷机，其耗电功率Pin表示成关于制冷功率Pwt与冷却水进口温度Tf的二元二次函数： Pin=S(PmtXn) ⑴ 于是电制冷机供冷产生的电费花销C为： C^YxgiPm,,J：")xpxt (2) 其中P为供冷时段电价，t为供冷时段时间长度，Y为电制冷机启停标志，If)为 步骤si中拟合出的电制冷机耗电功率关于制冷功率和冷却水进口温度的函数； 在忽略蓄冰槽储冷量随储存时间损耗的前提下，为了考量蓄冰槽供冷的花销，引入蓄 冰槽冷量均价Ps概念，描述单位冷量的价格，这个价格与充冷量、充冷时段电价和双工况电 制冷机制冰工况输入输出关系相关，由下式定义：(3) 其中Pinit为冷量均价初值，W init为蓄冰槽冷量初值，tin为冲冷时间，t °ut为放冷时间，P m为双工况电制冷机总制冰功率，Ps为蓄冰槽供冷功率，Y s为蓄冰槽供冷启停标志；由蓄冰槽 冷量均价的定义看出，在不考虑储冷量随时间损耗的前提下，均价与蓄冰槽是否放冷无关， 只在蓄冰槽冲冷时可能改变； 蓄冰槽储冷量增量ΛW为： Δ W = tinP-t°utPsYs (4) 目标函数为各种电制冷机空调工况供冷与蓄冰槽供冷所需电费之和：(5) P为供冷时段电价，t为供冷时段时间长度，Yi为第i台电制冷剂启停标志， 为第i台电制冷机耗电功率关于制冷功率和冷却水进口温度的函数，其中i为第i台电 制冷机供冷功率，T。，，为第i台电制冷机冷却水进口温度，Y s为蓄冰槽供冷启停标志，p 3为 蓄冰槽冷量均价，Ps为蓄冰槽供冷功率，Y s为蓄冰槽供冷启停标志； 约束条件包括供冷功率与预测负荷平衡、各电制冷机额定功率限制、蓄冰槽储冷量上 下限和冰蓄冷系统工作模式唯一：(6) 其中Yi为第i台电制冷剂启停标志，P 为第i台电制冷机供冷功率，P s为蓄冰槽供 冷功率，Ys为蓄冰槽供冷启停标志，Pp为预测负荷，a (〇, 1)为第i台电制冷机最小输 出功率比例，Pi rated为第i台电制冷机额定输出功率，W为蓄冰槽储冷量，Wmax和Wmin为蓄冰 槽储冷量上下限；Yic^为双工况制冷机制冰工况启停标志，默认制冰工况时所有双工况冷机 同时开启，加快制冰速度；YaiM为第i台双工况制冷机空调工况启停标志，共N台双工况制 冷机； 最后，在形成稳态经济优化问题的目标函数和约束条件的前提下对上述优化问题进行 求解，得到各电制冷机和蓄冰槽在经济效益最优条件下的启停状态和供冷功率。4.如权利要求1所述的联合供冷系统的递阶优化方法，其特征在于，步骤S3中，所述标 耦合的协调分布式预测控制方法，具体如下： 分布式控制系统按照空调工况电制冷机数量分为N个子系统，每个子系统由模型预测 控制（MPC)和相应控制的电制冷机组成，不同子系统MPC的数据使用通信器进行实时交换 完成计算任务，从而使得各台空调工况电制冷机供冷功率跟踪步骤S2中得到的经济效益 最优条件下的供冷功率并且总供冷功率跟踪各台空调工况电制冷机经济效益最优条件下 的供冷功率之和； 电制冷机制冷功率动态过程模型传递函数如下：C7) τ为一阶惯性环节时间常数，Td为延时时间，则第s个子系统离散状态空间模型形式 为：(8) S = 1，2，L，N 其中us(k)为k时刻电制冷机输出功率设定值，即为步骤S3中优化问题的优化变量， PsGO为k+Ι时刻电制冷机实际输出功率，xs(k) = [ps(k)ps(k+l)L ps(k+nd)]T为状态变量， nd为电制冷机延时时间相对于离散系统采样时间的倍数;A s、Bs、Cs为状态空间方程中的系 数矩阵，其中As e R(M"+1)X(M"+1)，Bi e R(""+1|xl由离散系统采样时间和电制冷机动态过程时间 常数决定，R表示实数集，设采样时间为Λ t，则有：Bs= [0 L O 1-θ_Δ?/τ]τ (9) Cs= [I 0 L 0] 对于第s个子系统，s = 1，2, ···，N ;MPC协调策略包括电制冷机在未来P个采样时 刻组成的预测时域内预测输出功率尽可能接近步骤S2得到的单台空调工况电制冷机经济 效益最优条件下的供冷功率rs，并且为了控制所有电制冷机预测输出功率跟上各台空调工 况电制冷机经济效益最优条件下的供冷功率之和r，在目标函数中引入预测时域内各台空 调工况电制冷机总预测功率与经济效益最优条件下的供冷功率之和差距的平方和，因此不 同子系统MPC目标函数出现耦合；综上，第s个子系统MPC目标函数如下所示：其中，Us= [u s (k) Us (k+Ι)…us (k+P-1) ]τ为预测时域内电制冷机输出功率设定值序列， us(k)为k时刻电制冷机输出功率设定值，即为优化问题的优化变量；而Q, R e [〇, 1]为 权重系数，分别用来衡量单台电制冷机预测负荷跟踪程度与所有电制冷机总负荷跟踪程度 的重要性；rs(k+i)为（k+i)时刻第s台空调工况电制冷机经济效益最优条件下的供冷功 率，r(k+i)为（k+i)时刻各台空调工况电制冷机在经济效益最优条件下的供冷功率之和； + 〗㈨为第s台空调工况电制冷机在k时刻通过预测模型预测的其在（k+i)时刻的输 出功率，为除第s台空调工况电制冷机之外的其余空调工况电制冷机在k时 刻通过预测模型预测的其在（k+i)时刻的输出功率之和； MPC预测模型为：i = 1，2，L，P 象认+?μ)为在k时刻通过预测模型计算出的（k+i)时刻第8台空调工况电制冷机状 态空间方程中状态变量，/\认+ /_1幻为第S台空调工况电制冷机在k时刻通过预测模型预 测的其在（k+i)时刻的输出功率;AsH和A/为状态空间方程中系数矩阵的相应次幂； 对于第s个子系统MPC控制器，整个优化问题形式为：i = 1,2, L, P 其中，{us，min，us，max}，{p s，min，ps，max}为控制变量的边界约束。5.如权利要求1-4任一项所述的联合供冷系统的递阶优化方法，其特征在于，所述步 骤S4中，采取迭代方法求解步骤S3中的分布式预测控制问题，具体步骤如下： 步骤S41 :在k时刻，根据步骤S3中得到的每台空调工况电制冷机经济效益最优条件 下的供冷功率rs，初始化每个子系统在预测时域内的电制冷机输出功率设定值序列Us，令 迭代次数1 = 〇 :上述〃 !(A + 1)为第1次迭代在（k+Ι)时刻的电制冷机输出功率设定值，P为预测时域长 度； 步骤S42:根据MPC预测模型，计算出每个子系统MPC控制时域内输出量预测值 β、= [/(U +11幻八认+ 2 I幻…/)你+ PI幻]τ，并通知给其他N-I个子系统MPC ; 步骤S43 :对于每个子系统，求解步骤S3中的优化问题，得到最优解，即最优制冷机输 出功率设定值序列Uf; 步骤S44 :检查所有子系统MPC的收敛条件是否满足，即对给定的精度ε se R(s = 1，2, L，N)，是否均有|υΓ -K|| < & Q = 1 2,L , A〇，如果所有的子系统收敛条件满足，令每 个子系统最优控制量< 与最优制冷机输出功率设定值序列Uf相等，即U! 转步骤 S45;否则，4α为介于〇和1之间的常数，I = 1+1， 转步骤S42 ; 步骤S45:在k时刻选取控制量<=[1 O L 0]U〔（λ' = 1，2,···，Λ〇作用于相应的子系 统； 步骤S46 :滚动移位到下一个采样时刻，即k+Ι - k，返回步骤S41，重复上述过程。
【专利摘要】本发明提供了一种针对由常规电制冷机和冰蓄冷系统构成的联合供冷系统的递阶优化方法，首先根据各冷源经济模型和分时电价政策，采用混合整数规划方法，优化各冷源的开启状态及最优功率设定值，同时考虑到各冷源的动态调节特性以及物理约束，为了更好地提高系统动态性能，设计目标耦合的协调分布式预测控制方法，在分布式框架下重新优化各冷源的设定值，使得各冷源在动态过程中保证总负荷的同时尽可能跟踪最优制冷功率设定值。本发明分别采用混合整数规划和分布式预测控制方法定量解决联合供冷系统的稳态和动态调度问题，为供冷系统能源管理与优化调度提供合理的建议。
【IPC分类】G05B13/04
【公开号】CN104898422
【申请号】CN201510225593
【发明人】蔡旭, 刘楚晖, 郑毅
【申请人】上海交通大学
【公开日】2015年9月9日
【申请日】2015年5月5日
最新回复(0)