技术摘要：
本发明公开了基于动态时间规整的负荷生产时段检测方法和存储介质，解决了如何利用电力数据开展商业负荷生产时段时长的测算的问题。本发明包括基于处理后的待测负荷数据，形成测试负荷功率曲线，利用生产阶段的参考功率曲线，与所述测试负荷功率曲线进行时间尺度的最优  全部
背景技术：
2019年底新型冠状病毒引发的肺炎疫情发生以来，各地经济受到严重影响。各级 政府一边抓防疫，一边采集措施推动经济有序恢复。在开展商业体恢复情况监测和评价中， 商业负荷有效生产时段时长由于可以反映商业经营者的信心状态，是评价指标之一。通过 商户生产时段时长的测算及其变化情况的刻画，可以为政策制定者掌握城市商业恢复情况 及疫情防控提供有益参考。此外，作为包括商业负荷在内的电力负荷的典型特征，负荷的主 要生产时段也是开展负荷特征标签化、画像等方面研究的重要数据维度。 随着电网智能化进程的加快，智能电表的覆盖率不断提升，采集数据的广度和频 度不断增加，电力数据也呈海量增长趋势。在此背景下，基于负荷特征的分析由于可以更加 精准直观的反映电力负荷、客户的特点，受到越来越多的关注。面向电力负荷的画像技术， 融合了负荷模式，可在多种维度反映负荷特征，为提升客户优质服务水平、增强源网荷互动 水平奠定了基础。因此，利用电力数据开展商业负荷生产时段时长的测算，具有良好的基 础。
技术实现要素：
本发明所要解决的技术问题是：如何利用电力数据开展商业负荷生产时段时长的 测算，本发明提供了解决上述问题的基于动态时间规整的负荷生产时段检测方法和存储介 质。 本发明通过下述技术方案实现： 基于动态时间规整的负荷生产时段检测方法，步骤如下： S1：从用电信息采集系统中抽取待测负荷数据，编号待测负荷数据，对于每个待测 负荷数据进行清洗； S2：基于所述S1中的处理后的待测负荷数据，形成测试负荷功率曲线，利用生产阶 段的参考功率曲线，与所述测试负荷功率曲线进行时间尺度的最优规整路径计算，得到最 优规整路径表，计算过程包括对所述测试负荷功率曲线和参考功率曲线进行拉伸或压缩， 搜索所述最优规整路径表，计算得到测试负荷功率曲线的生产阶段，输出测试负荷功率曲 线的生产阶段的起止时间数据； S3：结合所述S2中测试负荷功率曲线的生产阶段的起止时间数据，起止两时间数 据相减取绝对值得到待测负荷的生产时段时长。 进一步地，所述S1中的对数据进行清洗包括空值处理、异常偏大处理； 其中，数据清洗中的空值和异常偏大数据的清洗处理方法为：将空值点或异常偏 大点数据替换为于时间早于空值点或异常偏大点的最临近功率数据；异常偏大点的识别方 4 CN 111583064 A 说　明　书 2/5 页 法为：若目标点功率数据>待测负荷运行容量，将目标点判定为异常偏大点。 假定测试负荷日功率曲线为Tsig＝{Tsigt} ,t∈[1,2,…,T]，测量点时间间隔Δ T，基于动态时间规整的负荷生产时段提取方法包括以下步骤： 进一步地，构建参考功率曲线Rsig＝{Rsigr},r∈[1,2,…,R]。 其中，rt_start为参考功率曲线设定的生产时段起点时刻t_start对应的数据点编 号，rt_end为终点时刻t_end对应的数据点编号，r为时刻，R为参考功率曲线Rsig的信号点数 量，u为归一化标记。 进一步地，在标幺化的基础上，对测试负荷曲线进行归一化处理，得到归一化的测 试功率曲线Tsigu＝{Tsigu,t},t∈[1,2,…,T]。归一化算法为： 其中，Tsigmax为测试负荷的最大功率，Tsigmin为最小功率。 进一步地，对归一化后的测试负荷功率曲线与参考功率曲线，利用动态时间规整 算法，求取最佳规整路径W＝{wk},wk＝(t,r) ,t∈[1,2,…,T],r∈[1,2,…,R]； 其中，动态时间规整算法获取最佳规整路径的流程如下： 构建原始信号测量点距离矩阵。对两个时间序列信号，计算所有信号测量点对之 间的距离，构建测量点距离矩阵。 CTsigu,Rsig＝{ct,r},t∈[1,2,…,T],r∈[1,2,…,R] 其中，T为归一化后的测试时间序列Tsigu的信号点数量，R为参考时间序列Rsig的 信号点数量。 信号点的距离采用欧氏距离，即： ct,r＝‖Tsigu,t-Rsigr‖ 构建累计距离矩阵DTsigu,Rsig＝{dt,r},t∈[1,2,…,T],r∈[1,2,…,R]。 累计距离矩阵元素dt,r代表的是从两条时间序列信号的起始点对(Tsigu,1,Rsig1) 到当前信号点对(Tsigu,t,Rsigr)之间的信号序列对应的规整后累计距离，也就是两规整信 号序列段的相似性度量。累计距离矩阵元素dt,r计算方法如下： 对于第1行元素： 对于第1列元素： 对于其他元素：dt,r＝ct,r min{dt-1,r,dt-1,r-1,dt,r-1},t∈[2,3,…,T],r∈[2,3,…, R]。 回溯搜索最佳规整路径W＝{wk} ,wk＝(t,r) ,t∈[1,2,…,T] ,r∈[1,2,…,R]。其 中，t代表Tsigu第t个信号点Tsigu,t的序号，r代表Rsig第r个信号点Rsigr的序号。 回溯过程如下： 从累计距离矩阵DTsigu,Rsig右上角元素dT,R开始，即两条时间序列信号的最末端，将 该点位置坐标(T,R)记录为规整映射的最后一个点wend；并将该点记录为当前规整映射点 (i,j)＝(T,R)； 5 CN 111583064 A 说　明　书 3/5 页 找寻当前规整映射点(i,j)的上一级映射点。根据累计距离矩阵的构建方法，当前 规整映射点的上级映射点应为其左侧(i-1,j)、左下方(i-1,j-1)、正下方(i,j-1)三个映射 点中累计距离的最小值，记录该最小值对应的时序信号点的位置坐标为(t,r)；并替换当前 规整映射点(i,j)＝(t,r)； 重复上一步骤，直到抵达两条时间序列信号的起始点，将位置坐标记录为规整映 射的第一个点wstart＝(1,1)。 由于累计距离矩阵的每一个元素，代表的是从两条时序信号的起始点到当前数据 点对的最佳规整路径对应的累计距离，按此方法计算得到终止点累计距离也是两条时序信 号最优规整路径对应的距离。因此，上述回溯得到的路径即是最佳规整路径，在最佳规整路 径W中，搜寻参考功率曲线的生产时段起始时刻rt_start对应的正序第一个数据点编号，记录 为wstart＝(tt_start,rt_start)；搜寻参考功率曲线的生产时段终点时刻rt_end对应的倒序第一 个数据点编号，记录为wend＝(tt_end ,rt_end)。则tt_start、tt_end为测试负荷功率曲线生产时段 起点和终点对应的数据点编号； 进一步地，根据测量点时间间隔ΔT，将测试负荷生产时段起点和终点对应的数据 点编号换算成为时间，得到待测负荷的营业起始时间和终止时间，计算得到待测负荷的生 产时段时长。 进一步的，一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程 序，所述计算机程序被处理器执行时实现本方法中的步骤。本方法的具体使用依赖大量计 算，因此优选的通过计算机程序来实现上述计算过程，所以任何包含本方法中所保护的步 骤的计算机程序及其存储介质也属于本申请的保护范围内。 本发明具有如下的优点和有益效果： 本发明可通过商户用电数据，快速准确地计算营业起始时刻和终止时刻，并据此 计算生产时段时长，为疫情防护和政府掌握商业恢复提供了有价值的参考。 附图说明 此处所说明的附图用来提供对本发明实施例的进一步理解，构成本申请的一部 分，并不构成对本发明实施例的限定。在附图中： 图1为本发明的流程示意图。 图2为本发明的最佳规整路径示意图。 图3为本发明的动态时间规整前的测试曲线与参考曲线图。 图4为本发明的动态时间规整后的负荷曲线图。