基于延迟嵌入分析的建筑动态时序数据重构方法及系统与流程

1.本发明属于智能电网技术领域，具体涉及一种基于延迟嵌入分析的建筑动态时序数据重构方法及系统。


背景技术：

2.建筑能源系统以温度、负荷为代表的时序数据分析主要以统计学方法为基础，在分析过程中通过人工经验选取特征数据开展分析与预测，侧重于寻找高维变量之间的非线性联系，对变量的时间演变过程考虑尚不足够，难以满足多场景多类型能源系统实时预测的通用性。
3.因此，基于上述技术问题需要设计一种新的基于延迟嵌入分析的建筑动态时序数据重构方法及系统。


技术实现要素：

4.本发明的目的是提供一种基于延迟嵌入分析的建筑动态时序数据重构方法及系统。
5.为了解决上述技术问题，本发明提供了一种建筑动态时序数据重构方法，包括：
6.获取建筑供暖系统供回水温度的历史运行数据；
7.对历史运行数据进行预处理；
8.确定相空间嵌入维数；
9.确定延迟时间量；
10.构建与原始状态空间同胚的相空间检测二次供水温度变化的工况点；以及
11.基于工况变化点和选取的延迟时间量获取工况变化开始到结束时间段内的供水温度序列和回水温度序列。
12.进一步，所述获取建筑供暖系统供回水温度的历史运行数据的方法包括：
13.搜集建筑供暖系统供回水温度的历史运行数据，包括二级网供水温度ts和回水温度tr。
14.进一步，所述对历史运行数据进行预处理的方法包括：
15.对历史运行数据进行缺失值处理、异常值处理以及通过数据平滑减少噪声；
16.缺失值处理为，以缺失数据的前一个数据和后一个数据的均值作为缺失数据的值；
17.异常值处理为，以历史运行数据作为样本，获取样本的均值μ和标准差σ，样本分布在(μ-3σ，μ+3σ)范围外的点判断为异常点；
18.通过数据平滑减少噪声为，将前后预设个数采集时间点的均值作为当前时间点的值：
[0019][0020]
其中，是t时刻的数据前后总计n个数据点的平均值；是在时刻采集的数据点；n是预设的计算数据点的个数。
[0021]
进一步，所述确定相空间嵌入维数的方法包括：
[0022]
选取任一维度数m，获取m维相空间的序列x1，x2，...，x
l
中任意两向量之差的欧式距离r
ij
(m)：
[0023]rij
(m)＝||x
i-xj||；
[0024]
根据一组数r0，调整r0大小，获取一组lnr0和lnc(r)的值，以获取饱和关联维数关系dm：
[0025][0026][0027]
式中，r0为限定在重构相空间中各点间最小距离与最大距离之间的取值，c(r)为关联积分；l为总相点数，h(x)为heaviside单位函数；
[0028][0029]
构建不同嵌入维数m下的lnr
0-lnc(r)关系图，若lnr0与lnc(r)的曲线中包含一部分直线，则这条直线的斜率为关联维数dm；随着m取值的增大，每一条曲线的斜率随着嵌入维数m的增加到一定m
′
时，不再增加，此时不变的斜率称为饱和关联维数系数，对应的m
′
为重构相空间的嵌入维数。
[0030]
进一步，所述确定延迟时间量的方法包括：
[0031]
x(t)是时序数据{x(t)|-∞＜t＜+∞}，作截尾函数：
[0032][0033]
式中，t为时刻，t为设定的时间周期；
[0034]
则x
t
(t)的傅立叶变换存在：
[0035][0036]
数据在时域向频域变换过程中：
[0037]
[0038]
式中，为功率谱密度；
[0039]
通过开展傅立叶变换的反变换，获得衡量时域能量分布的自相关函数r
x
(τ)：
[0040][0041]
获取时序数据的占据显著能量比例的频率库，以形成时序数据的长短周期特性。
[0042]
进一步，所述构建与原始状态空间同胚的相空间检测二次供水温度变化的工况点的方法包括：
[0043]
对一个变量进行观测，获取时间序列{xi|i＝1，2，...，n}，获取一个时间延迟量τ，以xi，x
i+τ
，x
i+2r
，...，x
i+(m-1
)τ为坐标轴，重构m维相空间；
[0044]
根据时间序列{xi|i＝1，2，...，n}重构相空间：
[0045][0046]
式中，m为相空间嵌入维数；τ为延迟时间；l为总相点数，l＝n-(m-1)τ；x1，x2，...，x
l
为相空间序列；
[0047]
在重构的相空间内，获取围绕各低维稳定可测相点的轨道相对于该点的李雅普诺夫指数λ，通过该指数反映所截取的时序数据是否属于工况突变点：
[0048][0049]
进一步，所述基于工况变化点和选取的延迟时间量获取工况变化开始到结束时间段内的供水温度序列和回水温度序列的方法包括：
[0050]
基于工况变化点和选取的时间延迟量τ获取工况变化开始到结束时间段内的供水温度序列tsi和回水温度序列tri：
[0051]
tsi＝[t
si
，t
s(i+1)
，...，t
s(i+τ)
]
t
；
[0052]
tri＝[t
ri
，t
r(i+1)
，...，t
r(i+τ)
]
t
。
[0053]
另一方面，本发明还提供一种建筑动态时序数据重构系统，包括：
[0054]
数据获取模块，获取建筑供暖系统供回水温度的历史运行数据；
[0055]
预处理模块，对历史运行数据进行预处理；
[0056]
嵌入维数获取模块，确定相空间嵌入维数；
[0057]
延迟时间量获取模块，确定延迟时间量；
[0058]
工况点构建模块，构建与原始状态空间同胚的相空间检测二次供水温度变化的工况点；以及
[0059]
序列获取模块，基于工况变化点和选取的延迟时间量获取工况变化开始到结束时间段内的供水温度序列和回水温度序列。
[0060]
本发明的有益效果是，本发明通过获取建筑供暖系统供回水温度的历史运行数
据；对历史运行数据进行预处理；确定相空间嵌入维数；确定延迟时间量；构建与原始状态空间同胚的相空间检测二次供水温度变化的工况点；以及基于工况变化点和选取的延迟时间量获取工况变化开始到结束时间段内的供水温度序列和回水温度序列，实现了对建筑能源系统供需双侧核心参数的数据重构，实现不同时间尺度下建筑用能特性的精准刻画，有效降低建筑能耗，提高能源利用率，为建筑能源系统的预测调控和能源管理提供技术支撑。
[0061]
本发明的其他特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点在说明书以及附图中所特别指出的结构来实现和获得。
[0062]
为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
[0063]
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0064]
图1是本发明的建筑动态时序数据重构方法的流程图；
[0065]
图2是本发明的建筑动态时序数据重构方法的示意图；
[0066]
图3是本发明的建筑动态时序数据重构系统的原理框图。
具体实施方式
[0067]
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0068]
如图1和图2所示，本实施例提供了一种建筑动态时序数据重构方法，包括：获取建筑供暖系统供回水温度的历史运行数据；对历史运行数据进行预处理；确定相空间嵌入维数；确定延迟时间量；构建与原始状态空间同胚的相空间检测二次供水温度变化的工况点；以及基于工况变化点和选取的延迟时间量获取工况变化开始到结束时间段内的供水温度序列和回水温度序列，实现了对建筑能源系统供需双侧核心参数的数据重构，实现不同时间尺度下建筑用能特性的精准刻画，有效降低建筑能耗，提高能源利用率，为建筑能源系统的预测调控和能源管理提供技术支撑
[0069]
在本实施例中，所述获取建筑供暖系统供回水温度的历史运行数据的方法包括：搜集建筑供暖系统供回水温度的历史运行数据；搜集建筑供暖系统供回水温度的历史运行数据，包括二级网供水温度ts和回水温度tr。
[0070]
在本实施例中，所述对历史运行数据进行预处理的方法包括：对二级网供回水历史数据进行数据预处理；对历史运行数据进行缺失值处理、异常值处理以及通过数据平滑减少噪声；缺失值处理为，以缺失数据的前一个数据和后一个数据的均值作为缺失数据的值；异常值处理为，使用3-sigma原则进行异常点检测，3-sigma方法的原理是在正态分布下
以历史运行数据作为样本，获取样本的均值μ和标准差σ，样本分布在(μ-3σ，μ+3σ)范围外的点判断为异常点，即样本分布在(μ-3σ，μ+3σ)范围的概率为0.9974，反之落入这个区间外的概率只有0.0026，在实际问题中认为这样的小概率事件不会发生，将这类点判定为异常点；通过数据平滑减少噪声为，原始数据中噪声较多，通过数据平滑的方式来消除噪声，将前后预设个数采集时间点的均值作为当前时间点的值：
[0071][0072]
式中，是t时刻的数据前后总计n个数据点的平均值；是在时刻采集的数据点；n是预设的计算数据点的个数。
[0073]
在本实施例中，所述确定相空间嵌入维数的方法包括：采取饱和关联维数(g-p)法确定相空间嵌入维数；选取任一维度数m，获取m维相空间的序列x1，x2，...，x
l
中任意两向量之差的欧式距离r
ij
(m)：
[0074]rij
(m)＝||x
i-xj||；
[0075]
给定一组数r0(取值介于r
ij
(m)的最小值与最大值之间)，调整r0大小，获取一组lnr0和lnc(r)的值，以获取饱和关联维数关系dm：
[0076][0077][0078]
式中，r0为限定在重构相空间中各点间最小距离与最大距离之间的取值，c(r)为关联积分；l为总相点数，h(x)为heaviside单位函数；
[0079][0080]
x为函数的自变量；绘制不同嵌入维数m下的lnr
0-lnc(r)关系图，若lnr0与lnc(r)的曲线中包含一部分直线，则这条直线的斜率为关联维数dm；随着m取值的增大，每一条曲线的斜率随着嵌入维数m的增加到一定m
′
时，便不会增加，此时不变的斜率称为饱和关联维数系数，对应的m
′
为重构相空间的嵌入维数。
[0081]
在本实施例中，所述确定延迟时间量的方法包括：依据功率谱与时域自相关系数等内生特征综合评价确定延迟时间量；
[0082]
x(t)是时序数据{x(t)|-∞＜t＜+∞}，作截尾函数：
[0083][0084]
式中，t为时刻，t为设定的时间周期；
[0085]
则x
t
(t)的傅立叶变换存在：
[0086][0087]
由于数据在时域向频域变换过程中，其表征信息量的广义信号能量保持不变，对于这一数据，其变换前后有：
[0088][0089]
式中，为功率谱密度，表征频域内能量分布；
[0090]
类似的，通过开展傅立叶变换的反变换，获得衡量时域能量分布的自相关函数r
x
(τ)：
[0091][0092]
基于此可得到时序数据的占据显著能量比例的频率库，进一步形成时序数据的长短周期特性。
[0093]
在本实施例中，所述构建与原始状态空间同胚的相空间检测二次供水温度变化的工况点的方法包括：对某一个变量进行观测，获取时间序列{xi|i＝1，2，...，n}，再适当选取一个时间延迟量τ，以xi，x
i+τ
，x
i+2τ
，...，x
i+(m-1)τ
为坐标轴，重构m维相空间；
[0094]
根据时间序列{xi|i＝1，2，...，n}重构相空间：
[0095][0096]
式中，m为相空间嵌入维数；τ为延迟时间；l为总相点数，l＝n-(m-1)τ；x1，x2，...，x
l
为相空间序列；
[0097]
在重构的相空间内，进一步获取围绕各低维稳定可测相点的轨道相对于该点的李雅普诺夫指数λ，通过该指数反映所截取的时序数据是否属于工况突变点：
[0098][0099]
在本实施例中，所述基于工况变化点和选取的延迟时间量获取工况变化开始到结束时间段内的供水温度序列和回水温度序列的方法包括：基于工况变化点和选取的时间延迟量τ获取工况变化开始到结束时间段内的供水温度序列tsi和回水温度序列tri：
[0100]
tsi＝[t
si
，t
s(i+1)
，...，t
s(i+τ)
]
t
；
[0101]
tri＝[t
ri
，t
r(i+1)
，...，t
r(i+τ)
]
t
；
[0102]
由于供水和回水温度的趋势之间存在延迟时间，如果将延迟时间本身添加到回水温度中，就有可能找到一个延迟时间τ，从而使两个序列之间的相关系数最高。为了找到这样的延迟时间，将回水温度的序列tri逐步前移，每次获得一个新的序列。设置一个最大移
动步数k，则通过前移获得的回水温度序列分别为tr
i+1
，tr
i+2
，
…
，tr
i+k
，然后以此构建回水温度时间迟滞矩阵v：
[0103][0104]
其中k是最大迟滞步数；在该方法中，最大迟滞步数k、温度变化预阈值ε和延迟时间量τ都需要根据专家知识或人工经验设定。
[0105]
如图3所示，在本实施例中，本实施例还在建筑动态时序数据重构方法的基础上提供一种建筑动态时序数据重构系统，包括：数据获取模块，获取建筑供暖系统供回水温度的历史运行数据；预处理模块，对历史运行数据进行预处理；嵌入维数获取模块，确定相空间嵌入维数；延迟时间量获取模块，确定延迟时间量；工况点构建模块，构建与原始状态空间同胚的相空间检测二次供水温度变化的工况点；以及序列获取模块，基于工况变化点和选取的延迟时间量获取工况变化开始到结束时间段内的供水温度序列和回水温度序列；各模块中可以编写、存储等记录有建筑动态时序数据重构方法中的各部分对应的方法；并且各模块可以单独存在或集成在控制模块中等。
[0106]
综上所述，本发明通过获取建筑供暖系统供回水温度的历史运行数据；对历史运行数据进行预处理；确定相空间嵌入维数；确定延迟时间量；构建与原始状态空间同胚的相空间检测二次供水温度变化的工况点；以及基于工况变化点和选取的延迟时间量获取工况变化开始到结束时间段内的供水温度序列和回水温度序列，实现了对建筑能源系统供需双侧核心参数的数据重构，实现不同时间尺度下建筑用能特性的精准刻画，有效降低建筑能耗，提高能源利用率，为建筑能源系统的预测调控和能源管理提供技术支撑。
[0107]
在本技术所提供的几个实施例中，应该理解到，所揭露的系统和方法，也可以通过其它的方式实现。以上所描述的系统实施例仅仅是示意性的，例如，附图中的流程图和框图显示了根据本发明的多个实施例的方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现方式中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。
[0108]
另外，在本发明各个实施例中的各功能模块可以集成在一起形成一个独立的部分，也可以是各个模块单独存在，也可以两个或两个以上模块集成形成一个独立的部分。
[0109]
所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。
而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read-only memory)、随机存取存储器(ram，random access memory)、磁碟或者光盘等各种可以存储程序代码的介质。
[0110]
以上述依据本发明的理想实施例为启示，通过上述的说明内容，相关工作人员完全可以在不偏离本项发明技术思想的范围内，进行多样的变更以及修改。本项发明的技术性范围并不局限于说明书上的内容，必须要根据权利要求范围来确定其技术性范围。