CN113240184A - 一种基于联邦学习的楼宇空间单元冷负荷预测方法及*** - Google Patents

一种基于联邦学习的楼宇空间单元冷负荷预测方法及*** Download PDF

Info

Publication number
CN113240184A
CN113240184A CN202110562321.0A CN202110562321A CN113240184A CN 113240184 A CN113240184 A CN 113240184A CN 202110562321 A CN202110562321 A CN 202110562321A CN 113240184 A CN113240184 A CN 113240184A
Authority
CN
China
Prior art keywords
data
value
weight
load
model
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.)
Granted
Application number
CN202110562321.0A
Other languages
English (en)
Other versions
CN113240184B (zh
Inventor
刘振杰
黄文君
胡斌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202110562321.0A priority Critical patent/CN113240184B/zh
Publication of CN113240184A publication Critical patent/CN113240184A/zh
Application granted granted Critical
Publication of CN113240184B publication Critical patent/CN113240184B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00Machine learning
    • G06N20/20Ensemble learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/004Artificial life, i.e. computing arrangements simulating life
    • G06N3/006Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/084Backpropagation, e.g. using gradient descent

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Business, Economics & Management (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Computational Linguistics (AREA)
  • Economics (AREA)
  • Molecular Biology (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Marketing (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Development Economics (AREA)
  • Medical Informatics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Air Conditioning Control Device (AREA)

Abstract

本发明公开了一种基于联邦学习的楼宇空间单元冷负荷预测方法及***,该方法通过分布在不同区位的神经网络模型分别训练对应区位楼宇空间单元冷负荷预测模型权重并计算对应区位空调负荷预测值与单位空调负荷预测值;获取有聚合价值的区位的模型权重、空调负荷预测值和单位空调负荷预测值,判断并归类不同区位模型,将归为一类的模型利用横向联邦学习方法聚合权值;利用纵向联邦学习方法将权值存在差异的网络负荷预测值聚合建立回归预测模型,根据回归预测模型对楼宇空间单元冷负荷进行预测;根据聚合的不同区位的模型权值更新对应区位的神经网络模型。本发明方法具有较强的泛化性,加速网络训练效率,降低实现空调负荷预测功能所需的通信成本。

Description

一种基于联邦学习的楼宇空间单元冷负荷预测方法及***
技术领域
本发明涉及分布式空调负荷预测技术,尤其涉及一种基于联邦学习的楼宇空间单元冷负荷预测方法及***。
背景技术
近年来,随着国内环保意识的逐渐加强,环保节能一直是企业乃至社会大环境追求的目标。我国在建筑方面的总能耗占据社会总能耗的30%,而暖通设备产生的能耗占建筑能耗的60%以上,因此通过优化暖通设备的运行负荷曲线来降低不必要的能耗可以在可观的程度上节约资源。对于运用冰蓄冷技术空调设备的企业来说,通过在电价较低时蓄冷取代在电价峰值时从电网获取电能的方法可以为企业节约可观的成本。
但为实现空调机组的优化运行,控制机构及时响应优化策略,需要提前对负荷需求进行准确地预测,进而保障空调设备的安全平稳运行。近年来,常用的空调负荷预测方法有如下几种:时间序列分析法、多元回归分析法和人工神经网络。时间序列分析法对于平稳楼宇空间单元冷负荷的预测效果较好,且计算简单,但针对因非时间因素对负荷大小产生影响的情况时,预测效果下降。多元回归分析法虽然可以将非时间因素考虑进去,但是模型泛化性较低,不同类型建筑物适用效果不同。人工神经网络是一新兴的预测方法,其泛化性强,并且当给予充足的样本时可以达到很高的预测精度,不过该方法运行速度较前两者慢,容易出现因预测速度不足导致不能够满足实时控制的问题,因此还需根据硬件计算单元改进常规的人工神经网络预测方法。
现今,绝大多数的大型楼宇都通过中央暖通设备调节温度,因此在获取楼宇的负荷预测的源数据时需要尽可能分散且均匀地布置传感器,如果传感器采集到的数据都通过通信的方式向预测单元传输时会产生一定的通信压力,并且布置如此多的通信设施成本较大。为了解决通信的开销问题,需要采用云边协同技术,由一边缘计算设备收集周围传感器的数据并进行运算,再将运算结果传输到云端计算设备,通过联邦学习技术优化人工神经网络的预测方法。
现有的一些空调负荷预测的实施方案中,大多数倾向于利用神经网络或者支持向量机预测,但都为单机计算,在规模越来越大的建筑发展趋势下,通信成本越来越高,因此需要一种联邦学习与云边协同技术完成空调负荷的预测。
发明内容
本发明的目的在于针对现有技术的不足,提供一种基于联邦学习的楼宇空间单元冷负荷预测方法及***。
本发明的目的是通过以下技术方案实现的:
第一方面,本发明提供一种基于联邦学习的楼宇空间单元冷负荷预测方法,包括:
步骤S1:获取楼宇不同区位与预测楼宇空间单元冷负荷相关的数据;
步骤S2:根据S1获取的数据,通过分布在不同区位的神经网络模型分别训练出对应区位的楼宇空间单元冷负荷预测模型权重并计算出对应区位空调负荷预测值与对应区位单位空调负荷预测值;
步骤S3:获取有聚合价值的区位的模型权重、空调负荷预测值和单位空调负荷预测值,判断并归类不同区位的模型,将归为一类的模型利用横向联邦学习方法聚合权值;
步骤S4:利用纵向联邦学习方法,将权值存在差异的网络负荷预测值聚合建立回归预测模型,根据回归预测模型对楼宇空间单元冷负荷进行预测;
步骤S5:根据聚合的不同区位的模型权值更新对应区位的神经网络模型。
进一步地,获取楼宇不同区位与预测楼宇空间单元冷负荷相关的数据为当前预测时刻T的室外空气温度、上一个时刻T-1的室外空气温度、当前时刻T室内人流量、上一个时刻T-1室内人流量、当前时刻T太阳辐射强度、上一个时刻T-1太阳辐射强度、当前时刻T相对湿度、当前时刻T的单元空调冷负荷、上一个时刻T-1的单元空调冷负荷、两个时刻前T-2的空调冷负荷和前一日期同时刻T-24的空调冷负荷。
进一步地,上述T时刻为预测时刻,上述预测步长周期为1小时,一个预测周期内获取60次采样值,通过以上数据源预测T时刻的空调冷负荷,进而求得空间单元冷负荷。
进一步地,所述步骤S2中,分别将从不同区位获取的数据进行预处理,处理为可直接输入到神经网络模型的输入数据;数据预处理方法具体包括以下步骤:
步骤S21:剔除采样错误数据。
首先每分钟采集一个数据点,获取1小时采样周期内采样数据的统计量F:
Figure BDA0003078007850000021
其中si为第i次采样点,
Figure BDA0003078007850000022
为采样点的均值,
Figure BDA0003078007850000023
为第i次采样点的回归拟合值,该回归拟合值由一次函数的回归方程获得,通过最小二乘法计算回归方程参数。
随后获取统计量后判断其与Fα(1,58)的大小,Fα(1,58)表示自由度为1、数据点为60的拒绝域的临界值,当F>Fα(1,58)时,则无需剔除数据,当F≤Fα(1,58)时,则剔除采样点中与采样点拟合值差值最大的点,再进行如上判断。
步骤S22:取采样值预测周期内的平均数值作为当前预测周期的代表值。
步骤S23:通过公式
Figure BDA0003078007850000031
归一化方法将采样数据无量纲化。
其中n1为剔除采样错误数据后采样点的数量。
进一步地,所述步骤S2中的神经网络模型的结构为BP神经网络结构,通过粒子群优化的方式更新BP神经网络的权值,粒子适应度V采取如下公式计算:
Figure BDA0003078007850000032
其中,n2为采样点类别数,yi为第i个时序单元冷负荷的真实值,oi为第i个时序单元冷负荷的预测值。
粒子的速度与位置更新公式为:
Figure BDA0003078007850000033
Figure BDA0003078007850000034
其中,
Figure BDA0003078007850000035
分别为第t次迭代第i个粒子的速度和位置,u为惯性权值,c1与c2为加速因子,c3为扰动系数,c3通常取(0,1)之间一常数,
Figure BDA0003078007850000036
为第t次迭代第i个粒子的极值,
Figure BDA0003078007850000037
为第t次迭代中总粒子的极值,
Figure BDA0003078007850000038
为第t次迭代第i个粒子的位置,r1与r2为介于(0,1)之间的随机数。
进一步地,所述步骤S3中的有聚合价值的区域的判断方法为:
当计算出的部署在不同区位的神经网络的权重更新值满足以下公式时则认为有聚合价值:
Figure BDA0003078007850000039
Figure BDA00030780078500000310
Figure BDA00030780078500000311
其中,wt为区位第t次迭代更新的权值,
Figure BDA0003078007850000041
为区位第t次迭代更新的权值梯度,mh为用来拟合的历史权值的个数,mh取值大于5,m为区位总数量,α为一学习率,α通常取(0,1]之间一常数。
进一步地,所述步骤S3中的不同区域的模型判断与归类方法包含以下步骤:
步骤S31:将来自不同区位的权值的历史值按照时序进行排列:
Figure BDA0003078007850000042
其中,
Figure BDA0003078007850000043
为第j个区位的第t次迭代的权值向量,h为用来判别关联性的历史权值的个数,h取值大于10。
步骤S32:将权值归一化:
Figure BDA0003078007850000044
步骤S33:针对被比较区位网络权值产生求差序列:
Figure BDA0003078007850000045
其中
Figure BDA0003078007850000046
为被比较区位o的归一化权值向量。
步骤S34:判断与被比较区位网络权重的关联度:
Figure BDA0003078007850000047
步骤S35:将以上步骤S31-S34遍历mv-1遍形成不同区位权值之间关联度的方阵,mv为有聚合价值的区位数量,将方阵中高于一关联性阈值的区位权值归为一类模型权值。
进一步地,所述步骤S3中的横向联邦学习的方法如下式:
Figure BDA0003078007850000048
其中,
Figure BDA0003078007850000049
为第t轮迭代后聚合的权值,ms为属于同一类模型的区位数量。
进一步地,当进行多轮迭代后,如果
Figure BDA00030780078500000410
出现不收敛的情况时,聚合方法切换为:
Figure BDA00030780078500000411
其中,λ为一速度常数,通常取(0,1)之间一常数。
进一步地,所述步骤S4中的回归预测模型建立方法包括以下步骤:
步骤S41:对于来自不同模型的空调冷负荷预测数据,将其看作是楼宇总空调冷负荷的多元回归自变量,建立两者的多元线性回归模型:
Figure BDA0003078007850000051
其中,βl为第l类区位预测值的回归系数,xl为第l类区位的空调负荷预测值之和,β0为回归模型的常数项系数,ε为噪声项,N为一种高斯分布,σ为白噪声的标准差,mk为不同种类区位数量,y为楼宇的实际总空调负荷。
步骤S42:求取回归系数,公式如下:
Figure BDA0003078007850000052
其中,m1为用来确定回归系数的每个区位的历史空调负荷预测值之和的数据个数,
Figure BDA0003078007850000053
为第l类区位第l1轮空调负荷预测值之和。
进一步地,所述步骤S4中预测方法为:获取不同区位的空调负荷预测值后,经回归预测模型计算出整栋楼宇的空调预测负荷,获取整栋楼宇的室内建筑面积,整栋楼宇的空调预测负荷与室内建筑面积之比即为楼宇空间单元冷负荷预测值。
进一步地,所述步骤S5中的更新不同区域对应的神经网络模型的方法为:将聚合的不同区位的模型权值替换对应区位的神经网络权值。
进一步地,以上模型在预测前需要通过历史采样数据进行训练,随后再通过实时性采样数据进一步修正预测模型的参数并实时预测。
第二方面,本发明提供一种基于联邦学习的楼宇空间单元冷负荷预测***,包括数据采集模块、边缘计算模块和云端计算模块;楼宇每个区位对应一个边缘计算模块,每个边缘计算模块连接有多个数据采集模块,多个边缘计算模块与一个云端计算模块相连。
数据采集模块,采集楼宇不同区位与预测楼宇空间单元冷负荷相关的数据;
边缘计算模块,接收并处理对应区位采集的数据,通过该区位的神经网络模型训练出该区位的楼宇空间单元冷负荷预测模型权重并计算该区位空调负荷预测值与对应区位单位空调负荷预测值;并判断本轮区位的模型权重是否有上传到云端计算模块的聚合价值;
云端计算模块,接收并处理由楼宇不同区位边缘计算模块发送的权值、负荷预测值和单位负荷预测值,判断并归类不同区位的模型,将归为一类的模型利用横向联邦学习方法聚合权值;利用纵向联邦学习方法,将权值存在差异的网络负荷预测值聚合建立回归预测模型,根据回归预测模型对楼宇空间单元冷负荷进行预测;将聚合的不同区位的模型权值回传到不同区位的边缘计算模块。
云端计算模块获取不同区位的边缘计算模块的训练结果并回传模型训练参数,以此减少通信设施及数据传输的成本,提高模型计算速度,保证实时做出相应的控制策略并执行。
进一步地,某一区位的多个数据采集模块采集到数据并转换为数字量后,向边缘计算模块发送模型输入源数据,边缘计算模块在接收到由数据采集模块发送的数据后,对接收的数据进行预处理,并进行训练或预测,边缘计算模块处理完数据后,由云端计算模块聚合不同区位的处理结果,并将聚合后的信息反向发送给边缘计算模块。
进一步地,数据采集模块采集与预测空调冷负荷相关联的数据。数据采集模块每分钟采集一次室外空气温度、室内人流量、太阳辐射强度、相对湿度和空调冷负荷。一个区位的数据采集模块的室外空气温度采集单元、太阳辐射强度采集单元、相对湿度采集单元在该区位内均匀分布,采集同一区域内不同位置的数据,室内人流量采集单元统计该区位内人流量,空调冷负荷采集单元采集本区位消耗的空调冷负荷。数据采集完成后,数据采集模块将其转换为数字量并传输到对应区位的边缘计算模块。
进一步地,边缘计算模块接收并处理对应区位采集的数据并进行预测。边缘计算模块包含数据预处理单元、神经网络模型和权值上传判断单元。当边缘计算模块接收到由数据采集模块发送的源数据后,将数据进行预处理,预处理后的数据作为神经网络模型的输入数据,进而对该神经网络参数寻优。经权值上传判断单元聚合价值判断处理后,边缘计算模块将有聚合价值的网络权值与负荷预测值发送到云端计算模块,并接收由云端计算模块反向发送的权值调整数据。
进一步地,云端计算模块包含模型相似度判别单元、模型聚合单元和多元回归模型。当云端计算模块接收到由边缘计算模块发送的数据后,根据接收的更新权值判别相同类别预测模型的区位,针对同一类预测模型,利用横向联邦学习方法将更新权值聚合并将聚合后的权值反向发送到边缘计算单元。针对不同类预测模型的负荷预测值,通过建立回归预测模型,确定不同区位模型预测值对于总空调负荷的贡献度,云端计算模块计算得到总空调负荷后,获取总建筑面积并计算楼宇空间单元冷负荷预测值。
进一步地,所述基于联邦学习的楼宇空间单元冷负荷预测***首先根据历史数据确定所有预测模型的参数,进而随暖通***运行进行实时参数调整并实时预测楼宇空间单元冷负荷。
进一步地,所述边缘计算模块按照每小时的采样数据计算一次,数据采集单元完成60次数据采集后将采集数据发送给边缘计算模块,当边缘计算模块每计算k次后,将优化后权值与预测值传送到云端计算模块,云端计算模块聚合权值数据后向边缘计算模块返回聚合后的权值。
进一步地,上述单个边缘计算模块负责一个区位的数据处理,需要预测单位冷负荷的楼宇被划分为多个典型区位,当选取的区位与区位之间的差异较大时即被认定为选取了典型区位,所选取的区位可以不覆盖整栋楼宇。
进一步地,所述温度采集单元为温度传感器。
进一步地,所述太阳辐射采集单元为光敏传感器。
进一步地,所述相对湿度采集单元为湿度传感器。
进一步地,所述室内人流量采集单元为一红外线传感器,该传感器布置在人员进出对应区位的通道口,通过计算被遮挡的次数与时间获取人流量。
进一步地,所述空调冷负荷采集单元为温度传感器与风速检测仪。通过获取室内温度与该区位空调出风口的温度差、风口横截面积与风速获取空调冷负荷。
进一步地,所述边缘计算模块为MCU,所述云端计算模块为MCU。
综上所述,本发明具有以下有益效果:
本发明提供的一种基于联邦学习的楼宇空间单元冷负荷预测方法考虑了楼宇不同区位预测模型的差异问题,通过神经网络识别不同区位的模型,通过多元回归模型寻找不同类型的网络模型输出对总负荷的贡献度,与单一网络预测相比提高了预测速度,且该方法泛化性较强。通过联邦学习机制提高了神经网络训练效率,降低了通信成本,并在训练时发现收敛异常状况出现时,调整训练方法加强收敛的稳定性。
本发明提供的一种基于联邦学习的楼宇空间单元冷负荷预测***,减轻了传感设备与计算单元通信的压力。通过云边协同计算的架构增强了预测过程的实时性。只需挑选典型区位布置数据采集模块与边缘计算模块,节省了硬件成本。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例中的***框架图;
图2为本发明实施例中的方法流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚完整地描述,所描述的实施例仅仅是本发明一部分实施例,而不是全部的施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
如图1所示,本发明实施例提供的一种基于联邦学习的楼宇空间单元冷负荷预测***,包括数据采集模块、边缘计算模块和云端计算模块。
数据采集模块,采集楼宇不同区位与预测楼宇空间单元冷负荷相关的数据;
边缘计算模块,接收并处理对应区位采集的数据,通过该区位的神经网络模型训练出该区位的楼宇空间单元冷负荷预测模型权重并计算该区位空调负荷预测值与对应区位单位空调负荷预测值;并判断本轮区位的模型权重是否有上传到云端计算模块的聚合价值;
云端计算模块,接收并处理由楼宇不同区位边缘计算模块发送的权值、负荷预测值和单位负荷预测值,判断并归类不同区位的模型,将归为一类的模型利用横向联邦学习方法聚合权值;利用纵向联邦学习方法,将权值存在差异的网络负荷预测值聚合建立回归预测模型,根据回归预测模型对楼宇空间单元冷负荷进行预测;将聚合的不同区位的模型权值回传到不同区位的边缘计算模块;
其中,楼宇每个区位对应一个边缘计算模块,每个边缘计算模块连接有多个数据采集模块,多个边缘计算模块与一个云端计算模块相连。
单个边缘计算模块负责一个区位的数据处理,需要预测单位冷负荷的楼宇被划分为多个典型区位,当选取的区位与区位之间的差异较大时即被认定为选取了典型区位,所选取的区位可以不覆盖整栋楼宇。
位于同一区位的数据采集模块每分钟采集一次室外空气温度、室内人流量、太阳辐射强度、相对湿度和空调冷负荷。一个区域的数据采集模块的室外空气温度采集单元、太阳辐射强度采集单元、相对湿度采集单元在该区位内均匀分布,采集同一区域内不同位置的数据,室内人流量采集单元统计该区位内人流量,空调冷负荷采集单元采集本区位消耗的空调冷负荷。其中,温度采集单元为温度传感器。太阳辐射采集单元为光敏传感器。相对湿度采集单元为湿度传感器。室内人流量采集单元为一红外线传感器,该传感器布置在人员进出对应区位的通道口,通过计算被遮挡的次数与时间获取人流量。空调冷负荷采集单元为温度传感器与风速检测仪。通过获取室内温度与该区位空调出风口的温度差、风口横截面积与风速获取空调冷负荷。
针对于预测时刻T的空调冷负荷时,每小时获取以下大型建筑不同区位与预测楼宇空间单元冷负荷相关的数据:当前时刻T的室外空气温度、上一个时刻T-1的室外空气温度、当前时刻T室内人流量、上一个时刻T-1室内人流量、当前时刻T太阳辐射强度、上一个时刻T-1太阳辐射强度、当前时刻T相对湿度、当前时刻T的单元空调冷负荷、上一个时刻T-1的单元空调冷负荷、两个时刻前T-2的空调冷负荷和前一日期同时刻T-24的空调冷负荷。一个预测周期内所获取的每类数据有60组,每组数据之间相隔一分钟采样。
数据采集模块采集足够上述数据后发送到对应区位的边缘计算模块,边缘计算模块为一MCU。边缘计算模块接收到上述数据后对数据进行预处理,预处理包括错误采样值的剔除与无量纲化。
剔除错误采样数据的方法为:
首先每分钟采集一个数据点,获取1小时采样周期内采样数据的统计量F:
Figure BDA0003078007850000091
其中si为第i次采样点,
Figure BDA0003078007850000092
为采样点的均值,
Figure BDA0003078007850000093
为第i次采样点的回归拟合值。该回归拟合值由一次函数的回归方程获得,通过最小二乘法计算回归方程参数。
获取统计量后判断其与Fα(1,58)的大小,Fα(1,58)表示自由度为1、数据点为60的拒绝域的临界值,当F>Fα(1,58)时,则无需剔除数据,当F≤Fα(1,58)时,则剔除采样点中与采样点拟合值差值最大的点,再进行如上判断。
数据无量纲化的方法为:
利用公式
Figure BDA0003078007850000094
实现上述采样点的无量纲化。
其中n1为剔除采样错误数据后采样点的数量。
通过以上步骤将数据进行预处理后,作为BP神经网络的输入数据。利用粒子群优化方法调整BP神经网络权值。
首先依据BP神经网络的结构特点构造粒子个体长度:
L=S1*S2+S2*S3+S2+S3 (2)
其中,S1,S2和S3分别为BP神经网络的输入层、隐含层和输出层的节点参数。
其次计算粒子适应度V,适应度采取如下公式计算:
Figure BDA0003078007850000095
其中,n2为采样点类别数,yi为第i个时序单元冷负荷的真实值,oi为第i个时序单元冷负荷的预测值。
随后比较不同粒子的适应度,将具有最佳适应度的粒子记录为该种群本轮迭代的最佳粒子个体。
最后更新每个粒子的速度与位置,更新公式为:
Figure BDA0003078007850000101
Figure BDA0003078007850000102
其中,
Figure BDA0003078007850000103
分别为第t次迭代第i个粒子的速度和位置,u为惯性权值,c1与c2为加速因子,c3为扰动系数,c3通常取(0,1)之间一常数,
Figure BDA0003078007850000104
为第t次迭代第i个粒子的极值,
Figure BDA0003078007850000105
为第t次迭代中总粒子的极值,
Figure BDA0003078007850000106
为第t次迭代第i个粒子的位置,r1与r2为介于(0,1)之间的随机数。
通过迭代重复以上步骤,直至网络权值收敛,每更新一次权值便预测一次空调负荷预测值和单位空调负荷预测值。
不同区位的边缘计算模块在进行k轮上述计算后,首先判断更新的权值是否有聚合价值,当权重更新值满足以下公式时则认为有聚合价值:
Figure BDA0003078007850000107
Figure BDA0003078007850000108
Figure BDA0003078007850000109
其中,wt为区位第T次迭代更新的权值,
Figure BDA00030780078500001010
为区位第t次迭代更新的权值梯度,mh为用来拟合的历史权值的个数,mh取值大于5,m为区位总数量,α为一学习率,α通常取(0,1]之间一常数。
当满足上述判断准则后,边缘计算单元将模型权重、空调负荷预测值和单位空调负荷预测值发送到云端计算模块。
云端计算模块为MCU,云端计算模块接收到不同区位发送的模型权重、空调负荷预测值和单位空调负荷预测值后将模型权重进行聚合。
首先判断来自不同区位的模型重合情况:
判断模型是否重合方法为:将来自不同区位的权值历史值按照时序进行排列:
Figure BDA00030780078500001011
其中,
Figure BDA00030780078500001012
为第j个区位的第t次迭代的权值向量,h为用来判别关联性的历史权值的个数,h取值大于10。
随后将权值按归一化:
Figure BDA0003078007850000111
进而针对被比较区位网络权重产生求差序列:
Figure BDA0003078007850000112
最后判断与被比较区位网络权重的关联度:
Figure BDA0003078007850000113
将以上步骤遍历mv-1遍形成不同区位权值之间关联度的方阵,mv为有聚合价值的区位数量,将方阵中高于一关联性阈值的区位权值归为一类模型权值。关联性阈值的取值规则如表1。
表1
相关系数 相关程度
0.8~1.0 极强相关
0.6~0.8 强相关
0.4~0.6 中等强度相关
0.2~0.4 弱相关
0.0~0.2 极弱相关或无相关
当判断为同一类权值模型时,利用横向联邦学***均方法将同一类模型权值聚合:
Figure BDA0003078007850000114
其中,
Figure BDA0003078007850000115
为第t轮迭代后聚合的权值,ms为属于同一类模型的区位数量。
如果所有区位的模型都相同时,将接收到的单位空调负荷预测值平均后作为大型建筑的楼宇空间单元冷负荷预测值。
如果进行多轮迭代后,出现
Figure BDA0003078007850000116
不收敛的情况时,则切换聚合算法。判断是否收敛的方法为在经过多轮迭代后,
Figure BDA0003078007850000117
是否为一含有噪声的定值。当判断为不收敛时,聚合方法切换为:
Figure BDA0003078007850000121
其中,λ为一速度常数,通常取(0,1)之间一常数。
当云端计算单元检测到存在两种及两种以上的模型权值后,首先寻找属于同一模型的权值并通过公式(13)进行聚合,判别是否属于同一模型的方法还是通过计算公式(12)与关联性阈值进行比较。针对于来自不同模型的空调冷负荷预测数据,将其看作是楼宇总空调冷负荷的多元回归自变量,建立两者的多元线性回归模型:
Figure BDA0003078007850000122
其中,βl为第l类区位预测值的回归系数,xl为第l类区位的空调负荷预测值之和,β0为回归模型的常数项系数,ε为噪声项,N为一种高斯分布,σ为白噪声的标准差,mk为不同种类区位数量,y为楼宇的实际总空调负荷。
步骤S42:求取回归系数,公式如下:
Figure BDA0003078007850000123
其中,ml为用来确定回归系数的每个区位的历史空调负荷预测值之和的数据个数,
Figure BDA0003078007850000124
为第l类区位第l1轮空调负荷预测值之和。
通过上述方法确定回归模型系数后,计算出整栋楼宇的空调预测负荷,获取整栋楼宇的室内建筑面积,整栋楼宇的空调预测负荷与室内建筑面积之比即为楼宇空间单元冷负荷预测值。随后云端计算模块将聚合的不同区位的模型权值替换对应区位的神经网络权值。
该预测***首先根据历史数据确定所有预测模型的参数,进而随暖通***运行进行实时参数调整并预测楼宇空间单元冷负荷预测值。
本具体实施例仅仅是对本发明的解释,其并不是对本发明的限制,本领域技术人员在阅读完本说明书后可以根据需要对本实施例做出没有创造性贡献的修改,但只要在本发明的权利要求范围内都受到专利法的保护。

Claims (10)

1.一种基于联邦学习的楼宇空间单元冷负荷预测方法,其特征在于,包括:
S1:获取楼宇不同区位与预测楼宇空间单元冷负荷相关的数据;
S2:根据S1获取的数据,通过分布在不同区位的神经网络模型分别训练出对应区位的楼宇空间单元冷负荷预测模型权重并计算出对应区位空调负荷预测值与对应区位单位空调负荷预测值;
S3:获取有聚合价值的区位的模型权重、空调负荷预测值和单位空调负荷预测值,判断并归类不同区位的模型,将归为一类的模型利用横向联邦学习方法聚合权值;
S4:利用纵向联邦学习方法,将权值存在差异的网络负荷预测值聚合建立回归预测模型,根据回归预测模型对楼宇空间单元冷负荷进行预测;
S5:根据聚合的不同区位的模型权值更新对应区位的神经网络模型。
2.根据权利要求1所述的一种基于联邦学***均数值作为该预测周期的代表值,之后通过归一化方法将采样数据无量纲化;剔除采样错误数据的方法为:
首先每分钟采集一个数据点,获取1小时采样周期内采样数据的统计量F:
Figure FDA0003078007840000011
其中si为第i次采样点,
Figure FDA0003078007840000012
为采样点的均值,
Figure FDA0003078007840000013
为第i次采样点的回归拟合值,该回归拟合值由一次函数的回归方程获得,通过最小二乘法计算回归方程参数。
获取统计量F后判断其与Fα(1,58)的大小,Fα(1,58)表示自由度为1、数据点为60的拒绝域的临界值,当F>Fα(1,58)时,则无需剔除数据,当F≤Fα(1,58)时,则剔除采样点中与采样点回归拟合值差值最大的点,再进行如上判断。
利用公式
Figure FDA0003078007840000014
实现上述采样点的无量纲化。
其中n1为剔除采样错误数据后采样点的数量。
3.根据权利要求1所述的一种基于联邦学习的楼宇空间单元冷负荷预测方法,其特征在于,所述神经网络模型的结构为BP神经网络结构,通过粒子群优化的方式更新BP神经网络的权值;粒子适应度V采取如下公式计算:
Figure FDA0003078007840000021
其中,n2为采样点类别数,yi为第i个时序单元冷负荷的真实值,oi为第i个时序单元冷负荷的预测值。
粒子的速度与位置更新公式为:
Figure FDA0003078007840000022
Figure FDA0003078007840000023
其中,
Figure FDA0003078007840000024
分别为第t次迭代第i个粒子的速度和位置,u为惯性权值,c1与c2为加速因子,c3为扰动系数,c3取(0,1)之间一常数,Pi t为第t次迭代第i个粒子的极值,
Figure FDA0003078007840000025
为第t次迭代中总粒子的极值,
Figure FDA0003078007840000026
为第t次迭代第i个粒子的位置,r1与r2为介于(0,1)之间的随机数。
4.根据权利要求1所述的一种基于联邦学习的楼宇空间单元冷负荷预测方法,其特征在于,所述有聚合价值的区位的判断方法为:
当计算出的部署在不同区位的神经网络模型的权值更新值满足以下公式时,则认为有聚合价值:
Figure FDA0003078007840000027
Figure FDA0003078007840000028
Figure FDA0003078007840000029
其中,wt为区位第t次迭代更新的权值,
Figure FDA00030780078400000210
为区位第t次迭代更新的权值梯度,mh为用来拟合的历史权值的个数,mh取值大于5,m为区位总数量,α为学习率,α通常取(0,1]之间一常数。
5.根据权利要求1所述的一种基于联邦学习的楼宇空间单元冷负荷预测方法,其特征在于,所述的不同区位的模型判断与归类方法包含以下步骤:
先将来自不同区位的权值的历史值按照时序进行排列:
Figure FDA0003078007840000031
其中,
Figure FDA0003078007840000032
为第j个区位的第t次迭代的权值向量,h为用来判别关联性的历史权值的个数,h取值大于10。
随后将权值归一化:
Figure FDA0003078007840000033
进而针对被比较区位网络权重产生求差序列:
Figure FDA0003078007840000034
其中
Figure FDA0003078007840000035
为被比较区位o的归一化权值向量。
最后判断与被比较区位网络权重的关联度:
Figure FDA0003078007840000036
将以上公式遍历mv-1遍形成不同区位权值之间关联度的方阵,mv为有聚合价值的区位数量,将方阵中高于一关联性阈值的区位权值归为一类模型权值。
判断为同一类模型后,采用下式方法将同一类模型权值聚合:
Figure FDA0003078007840000037
其中,
Figure FDA0003078007840000038
为第t轮迭代后聚合的权值,ms为属于同一类模型的区位数量。
6.根据权利要求5所述的一种基于联邦学习的楼宇空间单元冷负荷预测方法,其特征在于,当进行多轮迭代后,如果
Figure FDA0003078007840000039
出现不收敛的情况时,聚合方法切换为:
Figure FDA00030780078400000310
其中,λ为速度常数,取(0,1)之间一常数。
7.一种基于联邦学习的楼宇空间单元冷负荷预测***,其特征在于,包括:
数据采集模块,采集楼宇不同区位与预测楼宇空间单元冷负荷相关的数据;
边缘计算模块,接收并处理对应区位采集的数据,通过该区位的神经网络模型训练出该区位的楼宇空间单元冷负荷预测模型权重并计算该区位空调负荷预测值与该区位单位空调负荷预测值;并判断本轮区位的模型权重是否有上传到云端计算模块的聚合价值;
云端计算模块,接收并处理由楼宇不同区位边缘计算模块发送的权值、负荷预测值和单位负荷预测值,判断并归类不同区位的模型,将归为一类的模型利用横向联邦学习方法聚合权值;利用纵向联邦学习方法,将权值存在差异的网络负荷预测值聚合建立回归预测模型,根据回归预测模型对楼宇空间单元冷负荷进行预测;将聚合的不同区位的模型权值回传到不同区位的边缘计算模块;
楼宇每个区位对应一个边缘计算模块,每个边缘计算模块连接有多个数据采集模块,多个边缘计算模块与一个云端计算模块相连。
8.根据权利要求7所述的一种基于联邦学习的楼宇空间单元冷负荷预测***,其特征在于:某一区位的多个数据采集模块采集到数据并转换为数字量后,向边缘计算模块发送模型输入源数据,边缘计算模块在接收到由数据采集模块发送的数据后,对接收的数据进行预处理,并进行训练或预测,边缘计算模块处理完数据后,由云端计算模块聚合不同区位的处理结果,并将聚合后的信息反向发送给边缘计算模块。
9.根据权利要求7所述的一种基于联邦学习的楼宇空间单元冷负荷预测***,其特征在于:所述边缘计算模块包含数据预处理单元、神经网络模型和权值上传判断单元。当边缘计算模块接收到由数据采集模块发送的源数据后,将数据进行预处理,预处理后的数据作为神经网络模型的输入数据,进而对该神经网络参数寻优。经权值上传判断单元聚合价值判断处理后,边缘计算模块将有聚合价值的网络权值与负荷预测值发送到云端计算模块,并接收由云端计算模块反向发送的权值调整数据。
10.根据权利要求7所述的一种基于联邦学习的楼宇空间单元冷负荷预测***,其特征在于:所述云端计算模块包含模型相似度判别单元、模型聚合单元和多元回归模型。当云端计算模块接收到由边缘计算模块发送的数据后,根据接收的更新权值判别相同类别预测模型的区位,针对同一类预测模型,利用横向联邦学习方法将更新权值聚合并将聚合后的权值反向发送到边缘计算单元。针对不同类预测模型的负荷预测值,通过建立回归预测模型,确定不同区位模型预测值对于总空调负荷的贡献度,云端计算模块计算得到总空调负荷后,获取总建筑面积并计算楼宇空间单元冷负荷预测值。
CN202110562321.0A 2021-05-21 2021-05-21 一种基于联邦学习的楼宇空间单元冷负荷预测方法及*** Active CN113240184B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110562321.0A CN113240184B (zh) 2021-05-21 2021-05-21 一种基于联邦学习的楼宇空间单元冷负荷预测方法及***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110562321.0A CN113240184B (zh) 2021-05-21 2021-05-21 一种基于联邦学习的楼宇空间单元冷负荷预测方法及***

Publications (2)

Publication Number Publication Date
CN113240184A true CN113240184A (zh) 2021-08-10
CN113240184B CN113240184B (zh) 2022-06-24

Family

ID=77138305

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110562321.0A Active CN113240184B (zh) 2021-05-21 2021-05-21 一种基于联邦学习的楼宇空间单元冷负荷预测方法及***

Country Status (1)

Country Link
CN (1) CN113240184B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114331761A (zh) * 2022-03-15 2022-04-12 浙江万胜智能科技股份有限公司 一种专变采集终端的设备参数分析调整方法及***
CN115037618A (zh) * 2022-06-06 2022-09-09 电子科技大学 一种基于KubeEdge的轻量级边缘智能协同联邦学习平台
CN115773562A (zh) * 2022-11-24 2023-03-10 杭州经纬信息技术股份有限公司 基于联邦学习的统一化暖通空调***故障检测方法
CN116187598A (zh) * 2023-04-28 2023-05-30 深圳市科中云技术有限公司 一种基于楼宇的虚拟电厂负荷预测方法

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106874581A (zh) * 2016-12-30 2017-06-20 浙江大学 一种基于bp神经网络模型的建筑空调能耗预测方法
CN107480815A (zh) * 2017-08-03 2017-12-15 国网河北省电力公司保定供电分公司 一种电力***台区负荷预测方法
CN107704875A (zh) * 2017-09-30 2018-02-16 山东建筑大学 基于改进ihcmac神经网络的建筑负荷预测方法和装置
CN111148066A (zh) * 2018-11-06 2020-05-12 丰田自动车株式会社 高网络负载场景中联网车辆的无线通信确保
CN111401658A (zh) * 2020-04-08 2020-07-10 西安建筑科技大学 一种基于建筑空间单元的并行冷负荷预测方法
US20200279169A1 (en) * 2019-03-01 2020-09-03 Government Of The United States Of America, As Represented By The Secretary Of Commerce Quasi-systolic processor and streaming batch eigenupdate neuromorphic machine
WO2020192896A1 (en) * 2019-03-26 2020-10-01 Huawei Technologies Co., Ltd. Apparatus and method for hyperparameter optimization of a machine learning model in a federated learning system
CN112000988A (zh) * 2020-08-28 2020-11-27 深圳前海微众银行股份有限公司 因子分解机回归模型构建方法、设备及可读存储介质
CN112183730A (zh) * 2020-10-14 2021-01-05 浙江大学 一种基于共享学习的神经网络模型的训练方法
CN112232608A (zh) * 2020-12-16 2021-01-15 浙江大学 一种基于数据驱动的区域建筑空调器短期耗电量预测方法
US20210042628A1 (en) * 2019-08-09 2021-02-11 International Business Machines Corporation Building a federated learning framework
CN112611080A (zh) * 2020-12-10 2021-04-06 浙江大学 一种基于联邦学习的智能空调控制***和方法
CN112634027A (zh) * 2020-12-30 2021-04-09 杭州趣链科技有限公司 一种面向小微企业信用评估的自适应联邦参数聚合方法
CN112686456A (zh) * 2020-12-31 2021-04-20 广东电网有限责任公司 一种结合边缘计算与用能识别的电力负荷预测***和方法
CN112800461A (zh) * 2021-01-28 2021-05-14 深圳供电局有限公司 一种基于联邦学习框架的电力计量***网络入侵检测方法

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106874581A (zh) * 2016-12-30 2017-06-20 浙江大学 一种基于bp神经网络模型的建筑空调能耗预测方法
CN107480815A (zh) * 2017-08-03 2017-12-15 国网河北省电力公司保定供电分公司 一种电力***台区负荷预测方法
CN107704875A (zh) * 2017-09-30 2018-02-16 山东建筑大学 基于改进ihcmac神经网络的建筑负荷预测方法和装置
CN111148066A (zh) * 2018-11-06 2020-05-12 丰田自动车株式会社 高网络负载场景中联网车辆的无线通信确保
US20200279169A1 (en) * 2019-03-01 2020-09-03 Government Of The United States Of America, As Represented By The Secretary Of Commerce Quasi-systolic processor and streaming batch eigenupdate neuromorphic machine
WO2020192896A1 (en) * 2019-03-26 2020-10-01 Huawei Technologies Co., Ltd. Apparatus and method for hyperparameter optimization of a machine learning model in a federated learning system
US20210042628A1 (en) * 2019-08-09 2021-02-11 International Business Machines Corporation Building a federated learning framework
CN111401658A (zh) * 2020-04-08 2020-07-10 西安建筑科技大学 一种基于建筑空间单元的并行冷负荷预测方法
CN112000988A (zh) * 2020-08-28 2020-11-27 深圳前海微众银行股份有限公司 因子分解机回归模型构建方法、设备及可读存储介质
CN112183730A (zh) * 2020-10-14 2021-01-05 浙江大学 一种基于共享学习的神经网络模型的训练方法
CN112611080A (zh) * 2020-12-10 2021-04-06 浙江大学 一种基于联邦学习的智能空调控制***和方法
CN112232608A (zh) * 2020-12-16 2021-01-15 浙江大学 一种基于数据驱动的区域建筑空调器短期耗电量预测方法
CN112634027A (zh) * 2020-12-30 2021-04-09 杭州趣链科技有限公司 一种面向小微企业信用评估的自适应联邦参数聚合方法
CN112686456A (zh) * 2020-12-31 2021-04-20 广东电网有限责任公司 一种结合边缘计算与用能识别的电力负荷预测***和方法
CN112800461A (zh) * 2021-01-28 2021-05-14 深圳供电局有限公司 一种基于联邦学习框架的电力计量***网络入侵检测方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114331761A (zh) * 2022-03-15 2022-04-12 浙江万胜智能科技股份有限公司 一种专变采集终端的设备参数分析调整方法及***
CN114331761B (zh) * 2022-03-15 2022-07-08 浙江万胜智能科技股份有限公司 一种专变采集终端的设备参数分析调整方法及***
CN115037618A (zh) * 2022-06-06 2022-09-09 电子科技大学 一种基于KubeEdge的轻量级边缘智能协同联邦学习平台
CN115037618B (zh) * 2022-06-06 2023-11-07 电子科技大学 一种基于KubeEdge的轻量级边缘智能协同联邦学习平台
CN115773562A (zh) * 2022-11-24 2023-03-10 杭州经纬信息技术股份有限公司 基于联邦学习的统一化暖通空调***故障检测方法
CN116187598A (zh) * 2023-04-28 2023-05-30 深圳市科中云技术有限公司 一种基于楼宇的虚拟电厂负荷预测方法
CN116187598B (zh) * 2023-04-28 2023-07-04 深圳市科中云技术有限公司 一种基于楼宇的虚拟电厂负荷预测方法

Also Published As

Publication number Publication date
CN113240184B (zh) 2022-06-24

Similar Documents

Publication Publication Date Title
CN113240184B (zh) 一种基于联邦学习的楼宇空间单元冷负荷预测方法及***
CN102705957B (zh) 办公建筑中央空调逐时冷负荷在线预测方法及***
CN106920006B (zh) 一种基于isoa-lssvm的地铁站空调***能耗预测方法
CN102779228B (zh) 商场建筑中央空调冷负荷在线预测方法及***
CN110332647B (zh) 地铁地下车站空调***负荷预测方法及空调***
CN107392368B (zh) 一种基于气象预报的办公建筑动态热负荷组合预测方法
CN111486554B (zh) 一种基于在线学习的空调温度无感化控制方法
CN103912966A (zh) 一种地源热泵制冷***优化控制方法
Alamin et al. An Artificial Neural Network (ANN) model to predict the electric load profile for an HVAC system
Song et al. An indoor temperature prediction framework based on hierarchical attention gated recurrent unit model for energy efficient buildings
CN116205508A (zh) 一种分布式光伏发电异常诊断方法和***
CN111815039A (zh) 基于天气分类的周尺度风电功率概率预测方法及***
CN114200839A (zh) 一种耦合环境行为动态监测的办公建筑能耗智能控制模型
CN113887833A (zh) 一种分布式能源用户侧逐时负荷预测方法及***
Yu et al. Predicting indoor temperature from smart thermostat and weather forecast data
Xiang et al. Prediction model of household appliance energy consumption based on machine learning
CN114861555A (zh) 一种基于Copula理论的区域综合能源***短期负荷预测方法
CN115563848A (zh) 一种基于深度学习的分布式光伏总辐射预测方法及***
CN115879190B (zh) 模型构建方法及装置、建筑负荷预测方法及装置
CN115510767B (zh) 基于深度时空网络的区域气温预测方法
CN110070215B (zh) 基于自适应的季节性arima模型的太阳能收集功率预测方法
Lu et al. A baseline model combining physics and data-driven approach for operation evaluation of district heating substation
CN115879190A (zh) 模型构建方法及装置、建筑负荷预测方法及装置
Chen et al. Integrated attention mechanism for GBDT building energy consumption prediction algorithm
CN117114179A (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
GR01 Patent grant
GR01 Patent grant