一种基于SBR和平均类间距离的场景数选择方法技术

本发明专利技术公开了一种基于SBR和平均类间距离的场景数选择方法，包括：步骤1、计算初始场景集中每两个场景组合之间的距离；步骤2、采用SBR算法进行多次场景缩减，进行场景缩减操作时满足以下条件：场景缩减前后场景集之间的概率距离最小；经场景缩减后得到的场景概率等于自身概率加上其他删除的场景对应概率之和；步骤3、得到目标场景数及其对应的概率；步骤4、根据场景概率矩阵进行场景初筛；步骤5、计算初筛后场景间的平均类间距离；步骤6、绘制提取场景数与平均类间距离之间的关系曲线，通过关系曲线来确定合适的场景数。本发明专利技术首次在SBR算法中提出平均类间距离的概念和计算方法，并基于场景出现概率辅助确定最终提取场景数。

A Scene Number Selection Method Based on SBR and Average Class Distance

【技术实现步骤摘要】
一种基于SBR和平均类间距离的场景数选择方法
本专利技术涉及聚类分析领域，尤其涉及一种基于SBR和平均类间距离的场景数选择方法。
技术介绍
聚类分析法主要可分为传统聚类分析法和人工智能聚类分析法两大类。传统聚类分析法主要包括划分聚类和层次聚类，层次聚类弥补了划分聚类的部分缺点，能够识别不同形状的时间序列。人工智能聚类算法主要利用了机器模拟人脑学习功能，代表方法包括：智能搜索聚类、并行聚类、分布式聚类算法等方法。人工智能算法很大程度上提高了聚类速度以及灵活性，但容易产生聚类结果不稳定，过度学习等问题。随着大数据时代的到来，聚类分析法被广泛应用于医疗、海洋以及电力系统等方面。总体看来，大多数聚类分析法主要采用特征变换处理原始数据，但基于特征选择的聚类分析法并不多。其中，同步回代缩减(SimultaneousBackwardReduction,SBR)是一种可以从大量场景中提取典型场景的分析方法。但是，SBR法在应用时，其提取场景数的选择缺乏一定客观性。因此，提出一种基于SBR和平均类间距离的场景数选择方法，对于完善和发展SBR法，具有比较重要的意义，并且可为电力部门各机组出力以及用户负荷过程典型场景提取提供技术支撑。
技术实现思路
本专利技术要解决的技术问题在于针对现有技术中的缺陷，提供一种基于SBR和平均类间距离的场景数选择方法。本专利技术解决其技术问题所采用的技术方案是：本专利技术提供一种基于SBR和平均类间距离的场景数选择方法，该方法包括以下步骤：步骤1、获取初始场景集数据，计算初始场景集中每两个场景组合之间的距离；步骤2、采用SBR算法进行多次场景缩减，进行场景缩减操作时满足以下条件：场景缩减前后场景集之间的概率距离最小；经场景缩减后得到的场景概率等于自身概率加上其他删除的场景对应概率之和；步骤3、进行多次场景缩减后，得到目标场景数及其对应的概率；步骤4、根据场景概率矩阵进行场景初筛；步骤5、计算初筛后场景间的平均类间距离；步骤6、绘制提取场景数与平均类间距离之间的关系曲线，通过关系曲线来确定合适的场景数。进一步地，本专利技术的步骤1中的具体方法为：记N天出力过程线初始场景集为X＝[X1,X2,...,XN]'，其中第i个场景为Xi＝[Xi1,Xi2,...,XiT]，i＝1,2,3…N；采样时间节点为t，t＝1,2,3…T；初始场景对应概率为P＝[P1,P2,...,PN]'，删除的场景为Y，删除场景数为N-M，最终得到的场景为X-Y；计算初始场景集中每两个场景组合之间的距离公式为：进一步地，本专利技术的步骤2的具体方法为：采用SBR算法进行第k次场景缩减，k＝1,2,3…N-M；为了确定合适的场景集X-Y使场景集X删除场景集Y后缩减效果最佳，要求场景缩减前后场景集之间的概率距离Dl最小，其计算公式为：其中，Pi表示初始场景Xi对应的概率；为保证场景缩减后，所有场景的概率之和仍为1，要求经场景缩减后得到的场景概率等于自身概率加上其他因距离接近而删除的场景对应概率之和；其公式为：其中，表示第k次场景缩减后场景Xj更新后的概率，其中：第k次场景缩减得到：将最小的对应的场景lk从场景集合X-Yk-1中删除，并入已删除的场景Yk-1中，得到Yk＝Yk-1∪{lk}；根据场景缩减原则将场景lk对应的概率加到场景集合X-Yk-1中与场景lk距离最近的场景对应概率上；再次计算更新后场景集X-Yk中每两个场景组合之间的距离。进一步地，本专利技术的步骤3的具体方法为：重复步骤2直至经过N-M次场景缩减后，最终得到M个场景及其对应的概率；场景Z及其对应的概率P’为：P'＝[P1',P2',...,PM']'。进一步地，本专利技术的步骤4的具体方法为：根据步骤3得到的场景Z及其对应的概率P’，若M个场景概率中某一场景对应概率低于0.05，则忽略该场景，场景提取数变为M-1，依次类推直至剩下的场景概率均大于0.05，最终得到满足条件的场景Z’以及对应场景概率P”；P″＝[P1″,P2″,...,PM′″]'。进一步地，本专利技术的步骤5中的具体方法为：计算M’个场景间的平均类间距离W(M’)，其公式为：其中，aef表示Z’中任意两个场景Z’e和Z’f之间的欧式距离，共有种组合方式。进一步地，本专利技术的步骤6的具体方法为：以M’为横坐标，W(M’)为纵坐标作出M’与W(M’)之间的关系图，选取最大的W(M’)对应的M’作为合适的场景数。本专利技术产生的有益效果是：本专利技术的基于SBR和平均类间距离的场景数选择方法，首次将平均类间距离这一指标应用到SBR场景数选择上，其意义在于：平均类间距离W(M’)越大，说明经过场景缩减后得到的场景集合之间的平均距离越大，差异性就越大，因而得到的结果更具有代表性；同时，本方法还采用场景出现概率来辅助确定最终缩减的场景数，避免最终场景集中因出现较多小概率场景造成对应平均类间距离过大的情况。本方法从相对客观角度改进了SBR算法，使得在实际场景缩减中场景缩减个数的确定更有依据，有助于SBR法的推广使用。附图说明下面将结合附图及实施例对本专利技术作进一步说明，附图中：图1是本专利技术实施例的方法流程图；图2是本专利技术实施例的以某电网的1月风电出力数据经SBR法缩减后得到效果图；图3是本专利技术实施例的1月份提取场景数与平均类间距离间的关系曲线。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本专利技术，并不用于限定本专利技术。如图1所示，本专利技术实施例包括以下步骤：步骤1，计算初始场景集中每两个场景组合之间的距离本专利技术以某电网2018年风电出力逐日24点资料为实例，事先对数据进行归一化处理，使得所有数据均在[0,1]区间内，取1月份31天(即31个场景)数据进行计算说明。记N天出力过程线初始场景集为X＝[X1,X2,...,XN]'，其中第i(i＝1,2,3…N)个场景为Xi＝[Xi1,Xi2,...,XiT]，采样时间节点为t(t＝1,2,3…T)，初始场景对应概率为P＝[P1,P2,...,PN]'，删除的场景为Y，删除场景数为N-M，最终得到的场景为X-Y。计算初始场景集中每两个场景组合之间的距离公式如下所示：步骤2，采用SBR算法进行第k(k＝1,2,3…N-M)次场景缩减，实现如下，为了确定合适的场景集X-Y使场景集X删除场景集Y后缩减效果最佳，要求场景缩减前后场景集之间的概率距离Dl最小，用公式表示如下：式中，Pi表示初始场景Xi对应的概率；为保证场景缩减前后，所有场景的概率之和仍为1，要求经场景缩减后得到的场景概率等于自身概率加上其他因距离接近而删除的场景对应概率之和。用公式表示如下：式中，表示第k次场景缩减后场景Xj更新后的概率，其中第k(k＝1,2,3…N-M)次缩减得到：将最小的对应的场景lk从场景集合X-Yk-1中删除，并入已删除的场景Yk-1中，得到Yk＝Yk-1∪{lk}。根据场景缩减原则将场景lk对应的概率加到场景集合X-Yk-1中与场景lk距离最近的场景对应概率上。再次计算更新后场景集X-Yk中每两个场景组合之间的距离。步骤3，重复步骤2直至经过N-M次场景缩减后，最终得到M个场景及其对应的概率，实现如下，场景数可以人为本文档来自技高网...

【技术保护点】
1.一种基于SBR和平均类间距离的场景数选择方法，其特征在于，该方法包括以下步骤：步骤1、获取初始场景集数据，计算初始场景集中每两个场景组合之间的距离；步骤2、采用SBR算法进行多次场景缩减，进行场景缩减操作时满足以下条件：场景缩减前后场景集之间的概率距离最小；经场景缩减后得到的场景概率等于自身概率加上其他删除的场景对应概率之和；步骤3、进行多次场景缩减后，得到目标场景数及其对应的概率；步骤4、根据场景概率矩阵进行场景初筛；步骤5、计算初筛后场景间的平均类间距离；步骤6、绘制提取场景数与平均类间距离之间的关系曲线，通过关系曲线来确定合适的场景数。

【技术特征摘要】
1.一种基于SBR和平均类间距离的场景数选择方法，其特征在于，该方法包括以下步骤：步骤1、获取初始场景集数据，计算初始场景集中每两个场景组合之间的距离；步骤2、采用SBR算法进行多次场景缩减，进行场景缩减操作时满足以下条件：场景缩减前后场景集之间的概率距离最小；经场景缩减后得到的场景概率等于自身概率加上其他删除的场景对应概率之和；步骤3、进行多次场景缩减后，得到目标场景数及其对应的概率；步骤4、根据场景概率矩阵进行场景初筛；步骤5、计算初筛后场景间的平均类间距离；步骤6、绘制提取场景数与平均类间距离之间的关系曲线，通过关系曲线来确定合适的场景数。2.根据权利要求1所述的基于SBR和平均类间距离的场景数选择方法，其特征在于，步骤1中的具体方法为：记N天出力过程线初始场景集为X＝[X1,X2,...,XN]'，其中第i个场景为Xi＝[Xi1,Xi2,...,XiT]，i＝1,2,3…N；采样时间节点为t，t＝1,2,3…T；初始场景对应概率为P＝[P1,P2,...,PN]'，删除的场景为Y，删除场景数为N-M，最终得到的场景为X-Y；计算初始场景集中每两个场景组合之间的距离公式为：3.根据权利要求2所述的基于SBR和平均类间距离的场景数选择方法，其特征在于，步骤2的具体方法为：采用SBR算法进行第k次场景缩减，k＝1,2,3…N-M；为了确定合适的场景集X-Y使场景集X删除场景集Y后缩减效果最佳，要求场景缩减前后场景集之间的概率距离Dl最小，其计算公式为：其中，Pi表示初始场景Xi对应的概率；为保证场景缩减后，所有场景的概率之和仍为1，要求经场景缩减后得到的场景概率等于自身概率加上其他因距离接近...

【专利技术属性】
技术研发人员：余姚果，梅亚东，王现勋，朱迪，吴贞晖，
申请(专利权)人：武汉大学，
类型：发明
国别省市：湖北,42