基于递推最小二乘法的空调所属建筑物一阶模型实时参数辨识方法
【技术领域】
[0001] 本发明属于电力系统调度部分的参数辨识技术，尤其涉及负荷聚合商对空调负荷 的实时参数辨识。
【背景技术】
[0002] 在硬件支持方面，负荷聚合商可以通过智能电网的高级测量终端和双向通信网络 获取用户空调的开关状态、室内外温度的实时变化情况，并对空调负荷进行实时控制与调 度。在在线辨识方法方面，目前常用的方法有最小二乘法、人工智能法、卡尔曼滤波法等，但 由于人工智能法其本身的复杂性，在实际应用中较少使用，而卡尔曼滤波法在迭代过程中 也要做大量计算工作，因此在多种参数辨识方法中，最小二乘法的原理简单，计算量少，被 广泛应用于各个领域的在线参数辨识。为了减少实时参数辨识过程中的计算量，一种非常 实用的改进的最小二乘法一递推最小二乘法被提了出来。依托于智能电网的硬件支持和相 关的在线辨识技术，负荷聚合商的空调负荷在线辨识对空调负荷参与需求响应起着至关重 要的作用。
[0003] 空调负荷作为一种重要的需求响应资源，其所属建筑物ETP模型已经被广泛应用 于空调负荷控制的各个领域，但其模型中的一些参数与建筑物的墙体的厚度、窗户面积、体 积大小等因素密切相关，无法通过测量获得，因此需要通过一定的参数辨识手段对空调负 荷模型的参数进行识别。但由于我国需求响应工作起步较晚，目前有关空调负荷模型的参 数辨识还未提及，忽视了空调负荷的实时在线建模。
【发明内容】
[0004] 发明目的：为了弥补现有需求响应过程中空调所属建筑物一阶模型参数在线识别 技术的空白，本发明提供一种基于递推最小二乘法的空调模型实时参数辨识方法，将空调 负荷的ETP模型等效为标准差分方程，确定输入输出序列及待辨识序列，通过智能电网的 高级计量终端对输入输出序列进行实时数据采集，并向负荷聚合商进行数据传递，负荷聚 合商利用递推最小二乘法对空调负荷模型的实时参数辨识，实时建立负荷聚合商管辖范围 内的空调负荷模型，用以进行空调负荷的实时控制与调度。
[0005] 技术方案：为实现上述目的，本发明采用的技术方案为：
[0006] -种基于递推最小二乘法的空调所属建筑物一阶模型实时参数辨识方法，实时建 立负荷聚合商管辖范围内的空调负荷模型，用以进行空调负荷的实时控制与调度；将空调 所属建筑物ETP模型等效为标准差分方程，同时确定输入输出序列及待辨识参数序列；通 过智能电网的高级计量终端对输入输出序列进行实时数据采集，并向负荷聚合商进行数据 传递，负荷聚合商利用递推最小二乘法对空调负荷模型的实时参数进行辨识。
[0007] 上述方法具体包括如下步骤：
[0008] (1)将空调所属建筑物ETP模型等效为单输入/输出线性系统，并使用差分方程的 标准形式进行描述，确定k时刻的输入输出序列为[u(k)，y(k)]，k时刻的随机变量序列为e(k)，待辨识参数数量为2n+l，待辨识参数为ap. .ajPb。. . .bn，负荷聚合商管辖范围内的 空调负荷模型为：
[0009]
[0010] 式中：n为差分方程的阶数；
[0011] (2)负荷聚合商通过智能电网的高级计量终端采集输入输出数据，并对采集的输 入输出数据进行预处理形成输入输出序列；
[0012] (3)获取k-2n-3~k-3时刻的N= 2n+l个输入输出序列，根据式⑴列出N个差 分方程为：
[0013] y(n+1) = -a^ (n)----any(1) +b〇u(n+1) +??? +bnu(1)+e(n+1)
[0014] y(n+2)=-aj(n+1)----any(2)+b0u(n+2) + …+bnu(2)+e(n+2) (2)
[0015] ......
[0016] y(n+N) =-aj(n+N_l)----any(N) +b〇u(n+N) +??? +bnu(N)+e(n+N)
[0017] 将式（2)写成矩阵的形式为：
[0018] Y(N) = 〇 (N)n(N)+E(N) (3)
[0019]式中：
[0020] Y(N) = [y(n+l),y(n+2), ???,y(n+N)]T (4)
[0021] n(N) = [a!，…，an,b0,…，bn]T (5)
[0022] E(N) =[e(n+1)，e(n+2)，...，e(n+N)]T (6)
[0023]
[0025] 式中：Y(N)为输出矩阵，〇 (N)为由输入输出序列组成的矩阵，n(N)为待辨识参 数矩阵，E(N)为误差矩阵；
[0026] (4)对采集的数据进行标准最小二乘法的参数估计，获取递推最小二乘法递推过 程的初始值：
[0027]
[0028]
[0029] 式中：P(N)为中间变量，/+V：)为经过标准最小二乘法求得的待辨识参数估计值；
[0030] (5)负荷聚合商提取k-2时刻的输入输出数据，并对输入输出数据进行预处理得 至ljk-l时刻的输入输出序列为[u(n+N+l)，y(n+N+l)];
[0031] (6)负荷聚合商通过k-2时刻的输入输出序列[u(n+N+l)，y(n+N+l)]和P(N)，对 修正系数进行计算：
[0032]
[0033]
[0034] (7)根据#以+ 1)计算更新后的待辨识参数估计值
[0035]
[0036] 根据更新式（1)，实时更新空调负荷模型；
[0037] (8)计算k-2 时刻的P(N+1):
[0038]
[0039] 然后令？(沁=卩(糾1)从=1^1，返回步骤（5)。
[0040] 优选的，所述步骤（2)和步骤（5)中，对采集的输入输出数据进行预处理，具体包 括以下内容：
[0041] (a)对于室内/室外温度，按照下述步骤进行预处理：
[0042] (al)由于室内/室外温度变化比较缓慢，根据统计规律，连续采集的五个点数值 变化应该不大，因此按照如下方法计算室内/室外温度的平均值f和标准差S:
[0043]
[0044]
[0045] 式中：Xl表示连续采集的第i个室内/室外温度；
[0046] (a2)判断|x,.-是否成立：若不成立，则说明Xi是正常数据；否则说明^是 坏数据，需要进行修正；
[0047] (a3)坏数据修正公式为：
[0048]
[0049] 式中：i；表示修正后的第i个室内/室外温度，a，0为自定义的修正系数，且 a +0 = 2；
[0050] (b)对于空调负荷的开关状态量，按照下述步骤进行预处理：
[0051] 在室内温度上升或者下降过程中，空调的开关状态保持不变，只有在到达 所设温度的上下限的时候空调开关才会发生变化；设连续采集的三个开关状态量为{XiXi,xi+1}，对应采集的三个室内温度为{y; y;，yi+1}，进行如下判断：
[0052] ①若X; : ?X; ?xi+1= 1，则表示x;是正常数据；
[0053] ②若Xii?xi+1= 1且x;i?Xi?xi+1= 0,则表示x;是坏数据，需要进行修正，令x; =Xi+1;
[0054] ③若X; : ?xi+1= 0,则：
[0055] 在（yiiiJ(yi+1-yi)>〇 时：若Xi: 1，则表示x;是正常数据；若x; : = 0， 则表示1;是坏数据，需要进行修正，令x;=x; 1;
[0056] 在（yi-y;D(7;+11;)〈0 时：若若Xi?xi+1= 1，则表示x;是正常数据；若x; ?xi+1 = 0,则表示1;是坏数据，需要进行修正，令x;=xi+1;
[0057] 式中：?表示逻辑与计算。
[0058] 有益效果：本发明提供的基于递推最小二乘法的空调所属建筑物一阶模型实时参 数辨识方法，可实时构建空调负荷的等效热参数模型，能够为相关部门提供空调负荷的相 关参数，为空调负荷参与电力系统的削峰、调频、备用、抑制可再生能源波动等提供了依据； 同时本发明弥补了国内空调负荷实时建模方面的研究不足，为其参与需求响应提供了技术 支撑。
【附图说明】
[0059] 图1为本发明方法的总流程图；
[0060] 图2为负荷聚合商系统的硬件组成图。
【具体实施方式】
[0061] 下面结合附图对本发明作更进一步的说明。
[0062] -种基于递推最小二乘法的空调所属建筑物一阶模型实时参数辨识方法，具体实 现流程如图1所示，下面就各个步骤加以具体说明。
[0063] ( -）将空调所属建筑物ETP模型等效为单输入/输出线性系统，并使用差分方程 的标准形式进行描述，确定k时刻的输入输出序列为[u(k)，y(k)]，k时刻的随机变量序列 为e(k)，待辨识参数数量为2n+l，待辨识参数为a:. . .ajPb.bn，负荷聚合商管辖范围内 的空调负何t旲型为：
[0064]
i-l
[0065] 式中：n为差分方程的阶数。
[0066] 式（1)的推导过程为如下：
[0067] (1. 1)空调所属建筑物ETP模型的离散方程为：
[0070] 式中：为t时刻室内温度（单位：r)，f为t时刻室外温度（单位：r)，at
[0068]
[0069] 表示采样时间间隔，e为散热系数，f= 5T。为时间常数；Q为制冷/制热量（单位:W)， R为等效热阻抗（单位：°C/W) ;s为开关状态量，S= 0表示关，s= 1表示开。
[0071](1. 2)根据式（1-1)和（1-2)，将空调负荷的差分方程进一步改写为：
[0072]
[0073] 式中：S(t)表示t时刻的开关状态量，S(t) = 0表示关，S(t) = 1表示开。
[0074] 以n= 1的情况为例进行推导，对式（1-3)进行如下的变量替换：
[0075] r(k) =T；'1 (1-4)
[0076] .v(々一1) = 7, ( 1-5
[0077] =Cl-6.)
[0078] u(k~l) =S(t) (1-7)
[0079] &1= -e(1-8)
[0080] b〇= 1-e(1-9)
[0081] bi=QR(1-e) (1-10)
[0082] 最终得到：
[0083] y(k) =-a^ (k~l)+b〇u(k)+b:u(k~l) (1-11)
[0084] (二）负荷聚合商通过智能电网的高级计量终端采集输入输出数据，并对采集的 输入输出数据进行预处理形成输入输出序列。
[0085](三）获取k-2n-3~k-3时刻的N= 2n+l个输入输出序列，根据式⑴列出N个 差分方程为：
[0086] y(n+1) = -a^ (n)----any(1) +b〇u(n+1) +??? +bnu(1)+e(n+1)
[0087] y(n+2)=-aj(n+1)----any(2)+b〇u(n+2) + …+bnu(2)+e(n+2) (2)
[0088] ......
[0089] y(n+N)=-aj(n+N_l)----any(N)+b〇u(n+N) + …+bnu(N)+e(n+N)
[0090] 将式（2)写成矩阵的形式为：
[0091] Y(N) = 〇 (N)n(N)+E(N) (3)
[0092] 式中：
[0093] Y(N