本发明涉及电网舞动技术领域，尤其涉及一种电网舞动关键影响因素聚类分析方法及系统。
背景技术：
架空线路舞动是影响冬季电网安全稳定运行的重要因素。舞动可能造成线路频繁跳闸，同时舞动还容易大面积发生，造成大面积停电等事故，舞动严重威胁了冬季电网的运行安全。开展准确的舞动预测预警工作，可为电网运检部门提前部署防舞措施，为调度部门提前制定电网防舞应急预案提供准确的科学指导。现有的舞动预警相关专利，如cn201510498199.x《一种输电通道舞动预警方法》和201410448630.5《一种基于adaboost的输电线舞动风险预警》等，均未对舞动的关键影响因素进行筛选，而是笼统的采用所有可能的舞动影响因素进行预警判断，可能出现次要影响因素对舞动主要影响因素的干扰从而影响舞动预测结果的准确程度。为提升舞动预测的准确率，首先必须全面分析舞动的影响因子，辨识出舞动关键影响因子。
通常，架空线路舞动同时受气象、地形等外部条件和线路型号、分裂数等线路结构内部条件的影响，而研究舞动的预测关键因素必须充分考虑外部条件和内部条件的共同影响，由于线路本体结构的影响要素种类多，相互之间影响关系复杂，迫切需要提出线路舞动关键影响因素的辨识方法，提高输电线路覆冰舞动预警准确率，提前部署输电线路防舞措施争取宝贵的时间提供重要技术支撑，提升舞动预防技术水平，保障冬季电网运行安全。
技术实现要素：
本发明目的在于一种电网舞动关键影响因素聚类分析方法及系统，以科学合理分析出电网舞动关键影响因素。
为实现上述目的，本发明公开一种电网舞动关键影响因素聚类分析方法，包括：
步骤s1、获取历史观测数据，所述历史观测数据包括架空线路历史舞动幅值数据以及相对应的舞动影响因素数据；
步骤s2、形成每条历史观测数据的点集s＝{li,xi1,xi2,…,xin},其中，历史观测数据的条数i＝1,2,…,m，舞动影响因子的个数为n；li为第i条历史观测数据中相应线路的历史舞动幅值，xij为对应li的第j个舞动影响因素观测值；
步骤s3、采用非监督聚类方法对所有历史观测数据的点集进行聚类，并设置待分析的舞动幅值区段，计算该舞动幅值区段所对应的聚类关联度；其中，所述聚类关联度为除法运算，分子为聚类所形成各类别中符合所述舞动幅值区段的点集个数，分母为包含有至少一个符合所述舞动幅值区段点集的所有类别的点集个数之和；以及
依次去掉各条历史观测数据所对应点集s中的同一舞动影响因子，并重新进行聚类并计算所述舞动幅值区段所对应的聚类关联度，若聚类关联度变小，则保留该舞动影响因子；若聚类关联度变大，则去掉该舞动影响因子；以此重复迭代，直至不能去掉点集s中所保留的任一舞动影响因素。
与上述方法相对应的，本发明公开一种电网舞动关键影响因素聚类分析系统，包括：
第一处理单元、用于获取历史观测数据，所述历史观测数据包括架空线路历史舞动幅值数据以及相对应的舞动影响因素数据；
第二处理单元、用于形成每条历史观测数据的点集s＝{li,xi1,xi2,…,xin},其中，历史观测数据的条数i＝1,2,…,m，舞动影响因子的个数为n；li为第i条历史观测数据中相应线路的历史舞动幅值，xij为对应li的第j个舞动影响因素观测值；
第三处理单元、用于采用非监督聚类方法对所有历史观测数据的点集进行聚类，并设置待分析的舞动幅值区段，计算该舞动幅值区段所对应的聚类关联度；其中，所述聚类关联度为除法运算，分子为聚类所形成各类别中符合所述舞动幅值区段的点集个数，分母为包含有至少一个符合所述舞动幅值区段点集的所有类别的点集个数之和；以及
依次去掉各条历史观测数据所对应点集s中的同一舞动影响因子，并重新进行聚类并计算所述舞动幅值区段所对应的聚类关联度，若聚类关联度变小，则保留该舞动影响因子；若聚类关联度变大，则去掉该舞动影响因子；以此重复迭代，直至不能去掉点集s中所保留的任一舞动影响因素。
本发明具有以下有益效果：
将聚类分析应用于分析电网舞动关键影响因素，并以聚类关联度筛选对应舞动幅值区段的关键影响因素；其中，相关的历史观测数据可以网罗各地区、各线路、各时期的所有舞动数据，通过大数据分析可确保分析结果的准确度；由于不同舞动幅值区段所对应的关键影响因素可能不同，通过本发明分区段筛选极大提高了分析的精度。
综上，本发明原理清晰，操作方便，具有很高的实用价值；而且，通过为输电线路舞动预测预警模型提供舞动关键影响因素，可有效减少舞动预测模型变量个数，消除次要影响因素对舞动测模型的干扰，提升架空线路舞动预报的准确率。
下面将参照附图，对本发明作进一步详细的说明。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1是本发明实施例公开的电网舞动关键影响因素聚类分析方法流程图。
具体实施方式
以下结合附图对本发明的实施例进行详细说明，但是本发明可以由权利要求限定和覆盖的多种不同方式实施。
实施例一
本实施例公开一种电网舞动关键影响因素聚类分析方法，如图1所示，包括：
步骤s1、获取历史观测数据，该历史观测数据包括架空线路历史舞动幅值数据以及相对应的舞动影响因素数据。
该步骤所获取的舞动影响因素的大类包括舞动发生时的气象因素、地形要素和线路结构因素等，各大类又可以细分成若干小类以确保包容潜在的所有舞动影响因素。例如：舞动发生时的气象因素包括覆冰厚度、气温、风速、风向、湿度、气压、降水等；地形因素包括海拔、地形、起伏度等；线路结构因素包括线路走向、导线线径、导线分裂数、杆塔档距等。
步骤s2、形成每条历史观测数据的点集s。
在该步骤中，s＝{li,xi1,xi2,…,xin}，其中，历史观测数据的条数i＝1,2,…,m，舞动影响因子的个数为n；li为第i条历史观测数据中相应线路的历史舞动幅值，xij为对应li的第j个舞动影响因素观测值。藉此，点集s为高维空间的点据集合，其中，每一个舞动影响因子作为该点集中的一维。
步骤s3、采用非监督聚类方法对所有历史观测数据的点集进行聚类，并设置待分析的舞动幅值区段，计算该舞动幅值区段所对应的聚类关联度；其中，所述聚类关联度为除法运算，分子为聚类所形成各类别中符合所述舞动幅值区段的点集个数，分母为包含有至少一个符合所述舞动幅值区段点集的所有类别的点集个数之和；以及依次去掉各条历史观测数据所对应点集s中的同一舞动影响因子，并重新进行聚类并计算所述舞动幅值区段所对应的聚类关联度，若聚类关联度变小，则保留该舞动影响因子；若聚类关联度变大，则去掉该舞动影响因子；以此重复迭代，直至不能去掉点集s中所保留的任一舞动影响因素。
可选的，在该步骤中，待分析的所述舞动幅值区段可为li＞l0，具体的聚类关联度的计算公式为：
其中，g为聚类的类别数，l0为舞动幅值判别阈值(例如：为0.5m)，subcount(li＞l0)为满足舞动条件点据的个数；allcount(li＞l0)为舞动点据所在的类别中所有点据的个数。
基于同样的原理，进一步地，本实施例还可以细化分析得出【0.5，1.0)、【1.0,1.5)……等不同舞动幅值区段所各自对应的舞动关键影响因素。
综上，本实施例原理清晰，操作方便，具有很高的实用价值；而且，通过为输电线路舞动预测预警模型提供舞动关键影响因素，可有效减少舞动预测模型变量个数，消除次要影响因素对舞动测模型的干扰，提升架空线路舞动预报的准确率。
实施例二
与上述方法实施例相对应的，本实施例公开一种电网舞动关键影响因素聚类分析系统，至少包括以下的第一至第三处理单元。
第一处理单元、用于获取历史观测数据，该历史观测数据包括架空线路历史舞动幅值数据以及相对应的舞动影响因素数据。
第二处理单元、用于形成每条历史观测数据的点集s＝{li,xi1,xi2,…,xin},其中，历史观测数据的条数i＝1,2,…,m，舞动影响因子的个数为n；li为第i条历史观测数据中相应线路的历史舞动幅值，xij为对应li的第j个舞动影响因素观测值。
第三处理单元、用于采用非监督聚类方法对所有历史观测数据的点集进行聚类，并设置待分析的舞动幅值区段，计算该舞动幅值区段所对应的聚类关联度；其中，该聚类关联度为除法运算，分子为聚类所形成各类别中符合该舞动幅值区段的点集个数，分母为包含有至少一个符合该舞动幅值区段点集的所有类别的点集个数之和；以及依次去掉各条历史观测数据所对应点集s中的同一舞动影响因子，并重新进行聚类并计算该舞动幅值区段所对应的聚类关联度，若聚类关联度变小，则保留该舞动影响因子；若聚类关联度变大，则去掉该舞动影响因子；以此重复迭代，直至不能去掉点集s中所保留的任一舞动影响因素。
本实施例中，例如：待分析的该舞动幅值区段为li＞l0，该第三处理单元具体的聚类关联度的计算公式可为：
其中，g为聚类的类别数，l0为舞动幅值判别阈值，subcount(li＞l0)为满足舞动条件点据的个数；allcount(li＞l0)为舞动点据所在的类别中所有点据的个数。
同理，本实施例公开的电网舞动关键影响因素聚类分析系统，原理清晰，操作方便，具有很高的实用价值；而且，通过为输电线路舞动预测预警模型提供舞动关键影响因素，可有效减少舞动预测模型变量个数，消除次要影响因素对舞动测模型的干扰，提升架空线路舞动预报的准确率。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。