CN113240184B - 一种基于联邦学习的楼宇空间单元冷负荷预测方法及*** - Google Patents
一种基于联邦学习的楼宇空间单元冷负荷预测方法及*** Download PDFInfo
- Publication number
- CN113240184B CN113240184B CN202110562321.0A CN202110562321A CN113240184B CN 113240184 B CN113240184 B CN 113240184B CN 202110562321 A CN202110562321 A CN 202110562321A CN 113240184 B CN113240184 B CN 113240184B
- Authority
- CN
- China
- Prior art keywords
- load
- data
- value
- weight
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 77
- 230000002776 aggregation Effects 0.000 claims abstract description 26
- 238000004220 aggregation Methods 0.000 claims abstract description 26
- 238000003062 neural network model Methods 0.000 claims abstract description 20
- 230000004931 aggregating effect Effects 0.000 claims abstract description 14
- 238000012549 training Methods 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims description 61
- 238000005070 sampling Methods 0.000 claims description 33
- 238000004378 air conditioning Methods 0.000 claims description 32
- 239000002245 particle Substances 0.000 claims description 24
- 238000001816 cooling Methods 0.000 claims description 21
- 238000013528 artificial neural network Methods 0.000 claims description 20
- 238000012545 processing Methods 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 238000007781 pre-processing Methods 0.000 claims description 7
- 238000009826 distribution Methods 0.000 claims description 4
- 238000012417 linear regression Methods 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- 238000006116 polymerization reaction Methods 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000013480 data collection Methods 0.000 claims description 3
- 238000012886 linear function Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 2
- 238000004891 communication Methods 0.000 abstract description 9
- 230000006870 function Effects 0.000 abstract 1
- 230000005855 radiation Effects 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 6
- 238000005265 energy consumption Methods 0.000 description 5
- 238000010438 heat treatment Methods 0.000 description 5
- 238000009423 ventilation Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000012731 temporal analysis Methods 0.000 description 2
- 238000000700 time series analysis Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012314 multivariate regression analysis Methods 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
- G06N20/20—Ensemble learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial 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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, 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:
随后获取统计量后判断其与Fα(1,58)的大小,Fα(1,58)表示自由度为1、数据点为60的拒绝域的临界值,当F>Fα(1,58)时,则无需剔除数据,当F≤Fα(1,58)时,则剔除采样点中与采样点拟合值差值最大的点,再进行如上判断。
步骤S22:取采样值预测周期内的平均数值作为当前预测周期的代表值。
其中n1为剔除采样错误数据后采样点的数量。
进一步地,所述步骤S2中的神经网络模型的结构为BP神经网络结构,通过粒子群优化的方式更新BP神经网络的权值,粒子适应度V采取如下公式计算:
其中,n2为采样点类别数,yi为第i个时序单元冷负荷的真实值,oi为第i个时序单元冷负荷的预测值。
粒子的速度与位置更新公式为:
其中,分别为第t次迭代第i个粒子的速度和位置,u为惯性权值,c1与c2为加速因子,c3为扰动系数,c3通常取(0,1)之间一常数,为第t次迭代第i个粒子的极值,为第t次迭代中总粒子的极值,为第t次迭代第i个粒子的位置,r1与r2为介于(0,1)之间的随机数。
进一步地,所述步骤S3中的有聚合价值的区域的判断方法为:
当计算出的部署在不同区位的神经网络的权重更新值满足以下公式时则认为有聚合价值:
进一步地,所述步骤S3中的不同区域的模型判断与归类方法包含以下步骤:
步骤S31:将来自不同区位的权值的历史值按照时序进行排列:
步骤S32:将权值归一化:
步骤S33:针对被比较区位网络权值产生求差序列:
步骤S34:判断与被比较区位网络权重的关联度:
步骤S35:将以上步骤S31-S34遍历mv-1遍形成不同区位权值之间关联度的方阵,mv为有聚合价值的区位数量,将方阵中高于一关联性阈值的区位权值归为一类模型权值。
进一步地,所述步骤S3中的横向联邦学习的方法如下式:
其中,λ为一速度常数,通常取(0,1)之间一常数。
进一步地,所述步骤S4中的回归预测模型建立方法包括以下步骤:
步骤S41:对于来自不同模型的空调冷负荷预测数据,将其看作是楼宇总空调冷负荷的多元回归自变量,建立两者的多元线性回归模型:
其中,βl为第l类区位预测值的回归系数,xl为第l类区位的空调负荷预测值之和,β0为回归模型的常数项系数,ε为噪声项,N为一种高斯分布,σ为白噪声的标准差,mk为不同种类区位数量,y为楼宇的实际总空调负荷。
步骤S42:求取回归系数,公式如下:
进一步地,所述步骤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:
获取统计量后判断其与Fα(1,58)的大小,Fα(1,58)表示自由度为1、数据点为60的拒绝域的临界值,当F>Fα(1,58)时,则无需剔除数据,当F≤Fα(1,58)时,则剔除采样点中与采样点拟合值差值最大的点,再进行如上判断。
数据无量纲化的方法为:
其中n1为剔除采样错误数据后采样点的数量。
通过以上步骤将数据进行预处理后,作为BP神经网络的输入数据。利用粒子群优化方法调整BP神经网络权值。
首先依据BP神经网络的结构特点构造粒子个体长度:
L=S1*S2+S2*S3+S2+S3 (2)
其中,S1,S2和S3分别为BP神经网络的输入层、隐含层和输出层的节点参数。
其次计算粒子适应度V,适应度采取如下公式计算:
其中,n2为采样点类别数,yi为第i个时序单元冷负荷的真实值,oi为第i个时序单元冷负荷的预测值。
随后比较不同粒子的适应度,将具有最佳适应度的粒子记录为该种群本轮迭代的最佳粒子个体。
最后更新每个粒子的速度与位置,更新公式为:
其中,分别为第t次迭代第i个粒子的速度和位置,u为惯性权值,c1与c2为加速因子,c3为扰动系数,c3通常取(0,1)之间一常数,为第t次迭代第i个粒子的极值,为第t次迭代中总粒子的极值,为第t次迭代第i个粒子的位置,r1与r2为介于(0,1)之间的随机数。
通过迭代重复以上步骤,直至网络权值收敛,每更新一次权值便预测一次空调负荷预测值和单位空调负荷预测值。
不同区位的边缘计算模块在进行k轮上述计算后,首先判断更新的权值是否有聚合价值,当权重更新值满足以下公式时则认为有聚合价值:
当满足上述判断准则后,边缘计算单元将模型权重、空调负荷预测值和单位空调负荷预测值发送到云端计算模块。
云端计算模块为MCU,云端计算模块接收到不同区位发送的模型权重、空调负荷预测值和单位空调负荷预测值后将模型权重进行聚合。
首先判断来自不同区位的模型重合情况:
判断模型是否重合方法为:将来自不同区位的权值历史值按照时序进行排列:
随后将权值按归一化:
进而针对被比较区位网络权重产生求差序列:
最后判断与被比较区位网络权重的关联度:
将以上步骤遍历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 | 极弱相关或无相关 |
当判断为同一类权值模型时,利用横向联邦学***均方法将同一类模型权值聚合:
如果所有区位的模型都相同时,将接收到的单位空调负荷预测值平均后作为大型建筑的楼宇空间单元冷负荷预测值。
其中,λ为一速度常数,通常取(0,1)之间一常数。
当云端计算单元检测到存在两种及两种以上的模型权值后,首先寻找属于同一模型的权值并通过公式(13)进行聚合,判别是否属于同一模型的方法还是通过计算公式(12)与关联性阈值进行比较。针对于来自不同模型的空调冷负荷预测数据,将其看作是楼宇总空调冷负荷的多元回归自变量,建立两者的多元线性回归模型:
其中,βl为第l类区位预测值的回归系数,xl为第l类区位的空调负荷预测值之和,β0为回归模型的常数项系数,ε为噪声项,N为一种高斯分布,σ为白噪声的标准差,mk为不同种类区位数量,y为楼宇的实际总空调负荷。
步骤S42:求取回归系数,公式如下:
通过上述方法确定回归模型系数后,计算出整栋楼宇的空调预测负荷,获取整栋楼宇的室内建筑面积,整栋楼宇的空调预测负荷与室内建筑面积之比即为楼宇空间单元冷负荷预测值。随后云端计算模块将聚合的不同区位的模型权值替换对应区位的神经网络权值。
该预测***首先根据历史数据确定所有预测模型的参数,进而随暖通***运行进行实时参数调整并预测楼宇空间单元冷负荷预测值。
本具体实施例仅仅是对本发明的解释,其并不是对本发明的限制,本领域技术人员在阅读完本说明书后可以根据需要对本实施例做出没有创造性贡献的修改,但只要在本发明的权利要求范围内都受到专利法的保护。
Claims (8)
1.一种基于联邦学习的楼宇空间单元冷负荷预测方法,其特征在于,包括:
S1:获取楼宇不同区位与预测楼宇空间单元冷负荷相关的数据;
S2:根据S1获取的数据,通过分布在不同区位的神经网络模型分别训练出对应区位的楼宇空间单元冷负荷预测模型权重并计算出对应区位空调负荷预测值与对应区位单位空调负荷预测值;
S3:获取有聚合价值的区位的模型权重、空调负荷预测值和单位空调负荷预测值,判断并归类不同区位的模型,将归为一类的模型利用横向联邦学习方法聚合权值;所述有聚合价值的区位的判断方法为:
当计算出的部署在不同区位的神经网络模型的权值更新值满足以下公式时,则认为有聚合价值:
所述的不同区位的模型判断与归类方法包含以下步骤:
先将来自不同区位的权值的历史值按照时序进行排列:
随后将权值归一化:
进而针对被比较区位网络权重产生求差序列:
最后判断与被比较区位网络权重的关联度:
将以上公式遍历mv-1遍形成不同区位权值之间关联度的方阵,mv为有聚合价值的区位数量,将方阵中高于一关联性阈值的区位权值归为一类模型权值;
判断为同一类模型后,采用下式方法将同一类模型权值聚合:
S4:利用纵向联邦学习方法,将权值存在差异的网络负荷预测值聚合建立回归预测模型,根据回归预测模型对楼宇空间单元冷负荷进行预测;所述回归预测模型建立方法包括以下步骤:
步骤S41:对于来自不同模型的空调冷负荷预测数据,将其看作是楼宇总空调冷负荷的多元回归自变量,建立两者的多元线性回归模型:
其中,βl为第l类区位预测值的回归系数,xl为第l类区位的空调负荷预测值之和,β0为回归模型的常数项系数,ε为噪声项,N为一种高斯分布,σ为白噪声的标准差,mk为不同种类区位数量,y为楼宇的实际总空调负荷;
步骤S42:求取回归系数,公式如下:
S5:根据聚合的不同区位的模型权值更新对应区位的神经网络模型。
2.根据权利要求1所述的一种基于联邦学***均数值作为预测周期的代表值,之后通过归一化方法将采样数据无量纲化;剔除采样错误数据的方法为:
首先每分钟采集一个数据点,获取1小时采样周期内采样数据的统计量F:
获取统计量F后判断其与Fα(1,58)的大小,Fα(1,58)表示自由度为1、数据点为60的拒绝域的临界值,当F>Fα(1,58)时,则无需剔除数据,当F≤Fα(1,58)时,则剔除采样点中与采样点回归拟合值差值最大的点,再进行如上判断;
其中n1为剔除采样错误数据后采样点的数量。
5.一种基于联邦学习的楼宇空间单元冷负荷预测***,其特征在于,包括:
数据采集模块,采集楼宇不同区位与预测楼宇空间单元冷负荷相关的数据;
边缘计算模块,接收并处理对应区位采集的数据,通过该区位的神经网络模型训练出该区位的楼宇空间单元冷负荷预测模型权重并计算该区位空调负荷预测值与该区位单位空调负荷预测值;并判断本轮区位的模型权重是否有上传到云端计算模块的聚合价值;判断区位聚合价值的过程如下:
当计算出的部署在不同区位的神经网络模型的权值更新值满足以下公式时,则认为有聚合价值:
云端计算模块,接收并处理由楼宇不同区位边缘计算模块发送的权值、负荷预测值和单位负荷预测值,判断并归类不同区位的模型,将归为一类的模型利用横向联邦学习方法聚合权值;利用纵向联邦学习方法,将权值存在差异的网络负荷预测值聚合建立回归预测模型,根据回归预测模型对楼宇空间单元冷负荷进行预测;将聚合的不同区位的模型权值回传到不同区位的边缘计算模块;
不同区位的模型判断与归类过程如下:
先将来自不同区位的权值的历史值按照时序进行排列:
随后将权值归一化:
进而针对被比较区位网络权重产生求差序列:
最后判断与被比较区位网络权重的关联度:
将以上公式遍历mv-1遍形成不同区位权值之间关联度的方阵,mv为有聚合价值的区位数量,将方阵中高于一关联性阈值的区位权值归为一类模型权值;
判断为同一类模型后,采用下式方法将同一类模型权值聚合:
回归预测模型的建立过程如下:
对于来自不同模型的空调冷负荷预测数据,将其看作是楼宇总空调冷负荷的多元回归自变量,建立两者的多元线性回归模型:
其中,βl为第l类区位预测值的回归系数,xl为第l类区位的空调负荷预测值之和,β0为回归模型的常数项系数,ε为噪声项,N为一种高斯分布,σ为白噪声的标准差,mk为不同种类区位数量,y为楼宇的实际总空调负荷;
求取回归系数,公式如下:
楼宇每个区位对应一个边缘计算模块,每个边缘计算模块连接有多个数据采集模块,多个边缘计算模块与一个云端计算模块相连。
6.根据权利要求5所述的一种基于联邦学习的楼宇空间单元冷负荷预测***,其特征在于:某一区位的多个数据采集模块采集到数据并转换为数字量后,向边缘计算模块发送模型输入源数据,边缘计算模块在接收到由数据采集模块发送的数据后,对接收的数据进行预处理,并进行训练或预测,边缘计算模块处理完数据后,由云端计算模块聚合不同区位的处理结果,并将聚合后的信息反向发送给边缘计算模块。
7.根据权利要求5所述的一种基于联邦学习的楼宇空间单元冷负荷预测***,其特征在于:所述边缘计算模块包含数据预处理单元、神经网络模型和权值上传判断单元;当边缘计算模块接收到由数据采集模块发送的源数据后,将数据进行预处理,预处理后的数据作为神经网络模型的输入数据,进而对该神经网络参数寻优;经权值上传判断单元聚合价值判断处理后,边缘计算模块将有聚合价值的网络权值与负荷预测值发送到云端计算模块,并接收由云端计算模块反向发送的权值调整数据。
8.根据权利要求5所述的一种基于联邦学习的楼宇空间单元冷负荷预测***,其特征在于:所述云端计算模块包含模型相似度判别单元、模型聚合单元和多元回归模型;当云端计算模块接收到由边缘计算模块发送的数据后,根据接收的更新权值判别相同类别预测模型的区位,针对同一类预测模型,利用横向联邦学习方法将更新权值聚合并将聚合后的权值反向发送到边缘计算单元;针对不同类预测模型的负荷预测值,通过建立回归预测模型,确定不同区位模型预测值对于总空调负荷的贡献度,云端计算模块计算得到总空调负荷后,获取总建筑面积并计算楼宇空间单元冷负荷预测值。
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 CN113240184A (zh) | 2021-08-10 |
CN113240184B true 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) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114331761B (zh) * | 2022-03-15 | 2022-07-08 | 浙江万胜智能科技股份有限公司 | 一种专变采集终端的设备参数分析调整方法及*** |
CN115037618B (zh) * | 2022-06-06 | 2023-11-07 | 电子科技大学 | 一种基于KubeEdge的轻量级边缘智能协同联邦学习平台 |
CN115773562A (zh) * | 2022-11-24 | 2023-03-10 | 杭州经纬信息技术股份有限公司 | 基于联邦学习的统一化暖通空调***故障检测方法 |
CN116187598B (zh) * | 2023-04-28 | 2023-07-04 | 深圳市科中云技术有限公司 | 一种基于楼宇的虚拟电厂负荷预测方法 |
Citations (13)
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 | 西安建筑科技大学 | 一种基于建筑空间单元的并行冷负荷预测方法 |
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 | 浙江大学 | 一种基于数据驱动的区域建筑空调器短期耗电量预测方法 |
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 | 深圳供电局有限公司 | 一种基于联邦学习框架的电力计量***网络入侵检测方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11651231B2 (en) * | 2019-03-01 | 2023-05-16 | Government Of The United States Of America, As Represented By The Secretary Of Commerce | Quasi-systolic processor and quasi-systolic array |
US20210042628A1 (en) * | 2019-08-09 | 2021-02-11 | International Business Machines Corporation | Building a federated learning framework |
-
2021
- 2021-05-21 CN CN202110562321.0A patent/CN113240184B/zh active Active
Patent Citations (13)
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 | 丰田自动车株式会社 | 高网络负载场景中联网车辆的无线通信确保 |
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 |
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 | 深圳供电局有限公司 | 一种基于联邦学习框架的电力计量***网络入侵检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113240184A (zh) | 2021-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113240184B (zh) | 一种基于联邦学习的楼宇空间单元冷负荷预测方法及*** | |
CN102705957B (zh) | 办公建筑中央空调逐时冷负荷在线预测方法及*** | |
CN106920006B (zh) | 一种基于isoa-lssvm的地铁站空调***能耗预测方法 | |
CN103912966B (zh) | 一种地源热泵制冷***优化控制方法 | |
CN110332647B (zh) | 地铁地下车站空调***负荷预测方法及空调*** | |
CN102779228B (zh) | 商场建筑中央空调冷负荷在线预测方法及*** | |
CN107392368B (zh) | 一种基于气象预报的办公建筑动态热负荷组合预测方法 | |
CN111486554B (zh) | 一种基于在线学习的空调温度无感化控制方法 | |
Alamin et al. | An Artificial Neural Network (ANN) model to predict the electric load profile for an HVAC system | |
CN111461921B (zh) | 一种基于机器学习的负荷建模典型用户数据库更新方法 | |
CN115758912A (zh) | 一种空调能耗优化*** | |
CN112288157A (zh) | 一种基于模糊聚类与深度强化学习的风电场功率预测方法 | |
CN116205508A (zh) | 一种分布式光伏发电异常诊断方法和*** | |
CN111815039A (zh) | 基于天气分类的周尺度风电功率概率预测方法及*** | |
CN113887833A (zh) | 一种分布式能源用户侧逐时负荷预测方法及*** | |
CN114200839A (zh) | 一种耦合环境行为动态监测的办公建筑能耗智能控制模型 | |
CN117553404A (zh) | 一种提高大型水冷中央空调***能效的方法及*** | |
Das et al. | A study on the application of artificial intelligence techniques for predicting the heating and cooling loads of buildings | |
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 | |
CN115563848A (zh) | 一种基于深度学习的分布式光伏总辐射预测方法及*** | |
CN113177675A (zh) | 一种基于天牛群算法优化神经网络的空调冷负荷预测方法 | |
CN115510767B (zh) | 基于深度时空网络的区域气温预测方法 | |
CN115879190B (zh) | 模型构建方法及装置、建筑负荷预测方法及装置 | |
Chen et al. | Integrated attention mechanism for GBDT building energy consumption prediction algorithm |
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 |