一种中低纬电网中地磁感应电流的识别方法技术

本申请涉及电力技术研究领域，尤其涉及一种中低纬电网中地磁感应电流的识别方法。该方法包括：对中低纬电网的样本数据进行聚类经验模态分解，获得第一本征模态函数组,该第一本征模态函数组包括不同时间尺度的本征模态函数。对所述第一本征模态函数组中的本征模态函数进行筛选和组合，获得第二本征模态函数组。对所述第二本征模态函数组进行希尔伯特变换，获得所述第二本征模态函数组的时频谱。根据所述时频谱，识别所述样本数据中是否具有地磁感应电流。该方法能够识别样本数据中较弱的地磁感应电流，并提高对地磁感应电流的识别精度。

【技术实现步骤摘要】
一种中低纬电网中地磁感应电流的识别方法
本申请涉及电力技术研究领域，尤其涉及一种中低纬电网中地磁感应电流的识别方法。
技术介绍
太阳的大尺度爆发性活动会对全球的电网系统造成广泛的影响。其影响机理是：携带高速等离子体的太阳风与地球磁层相互作用，引起地磁场短时间内的大幅度扰动，在跨越地面大范围区域内产生强大的地磁感应电场。该电场在地表的人工导电体中驱动的电流称为地磁感应电流(GeomagneticallyInducedCurrent,GIC)。当GIC通过接地的中性点流入变压器时会导致变压器铁芯半周饱和，在激励电流中产生大量谐波，不仅增加了变压器的无功消耗，还可能引起保护继电器误操作。同时，铁芯的高度饱和还会引起因金属结构件和油箱局部过热及绝缘老化，严重的情况下会导致变压器损毁而产生大面积断电事故。相比于高纬地区，中低纬地区关于GIC的研究尚存在很大差距。原因主要有以下几个方面：首先，超级磁暴出现的频率很低(约11的年周期)，而GIC的统计和预报研究需要大量样本事件。其次，中低纬GIC的平均水平比高纬低一个量级，在中低纬电网的电流信号中很可能无法直接获取到GIC。因此，中低纬电网GIC风险评估的前提是要能有效地从实际电流信号中识别出GIC。目前常用的GIC识别方法有：①自然正交向量分解法；②小波分析法。前者能成功地分解出地磁场中不同来源的特征向量，但主要适用于缓变信号，处理强磁暴扰动数据时结果不可靠。后者由于受Heisenberg测不准原理的制约，不能在时间和频率上同时达到较高的精度，加之小波基函数的选取有很大灵活性，不便于系统化的快速操作。由此可见，以上两种方法对于分析中低纬电网的GIC信号识别精度均不高。
技术实现思路
本申请提供了一种中低纬电网中地磁感应电流的识别方法，以提高地磁感应电流的识别精度。一种中低纬电网中地磁感应电流的识别方法，包括：步骤一、对中低纬电网的样本数据进行聚类经验模态分解，获得第一本征模态函数组，所述第一本征模态函数组包括不同时间尺度的本征模态函数；步骤二、对所述第一本征模态函数组中的本征模态函数进行筛选和组合，获得第二本征模态函数组；步骤三、对所述第二本征模态函数组进行希尔伯特变换，得到所述第二本征模态函数组的时频谱；所述时频谱为时间和频率的能量函数；步骤四、根据所述时频谱，识别所述样本数据中是否具有地磁感应电流。可选的，所述样本数据包括：中低纬电网的中性点电流信号，和与所述中性点电流信号对应的太阳风能量耦合函数。可选的，所述对所述第一本征模态函数组中的本征模态函数进行筛选和组合，获得第二本征模态函数组，包括：将所述第一本征模态函数组中的本征模态函数随机组合，获得至少一个本征模态函数组；将所述至少一个本征模态函数组与所述太阳风能量耦合函数进行相关性分析，获得最优的本征模态函数组；将与所述太阳风能量耦合函数相关性最优的本征模态函数组，作为第二本征模态函数组。可选的，所述根据所述时频谱，识别所述样本数据中是否具有地磁感应电流，包括：将所述时频谱划分为磁暴前相、磁暴初相、磁暴主相和磁暴恢复相四个阶段；将所述四个阶段分别对应的能量值沿相同长度的时间尺度积分，获得磁暴前相阶段能量指标、磁暴初相阶段能量指标、磁暴主相阶段能量指标和磁暴恢复相阶段能量指标；根据所述能量指标，识别所述样本数据中是否具有地磁感应电流。可选的，所述根据所述能量指标，识别所述样本数据中是否具有地磁感应电流，包括：分别计算所述磁暴初相阶段能量指标、所述磁暴主相阶段能量指标和所述磁暴恢复相阶段能量指标与所述磁暴前相阶段能量指标的比值；根据所述比值识别所述样本数据中是否具有地磁感应电流。本申请提供的技术方案包括以下有益技术效果：在本申请提供的方法中，首先，聚类经验模态分解能够依据样本数据自身的时间尺度特征来进行分解，无须预先设定任何基函数，有利于处理非平稳及非线性数据。其次，对本征模态函数进行筛选和组合，能够去除样本数据中，太阳辐射昼夜变换在电离层产生的电流，以减小该电流对GIC识别的影响。最后，通过希尔伯特变换可以获得时频谱，该时频谱能准确地表达电流信号在时间和频率空间的各类信息，能够识别样本数据中较弱的GIC并提高对GIC的识别精度。附图说明为了更清楚地说明本申请的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，对于本领域普通技术人员而言，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1为本申请实施例提供的一种中低纬电网中地磁感应电流的识别方法流程图。图2为本申请实施例提供的聚类经验模态分解流程图。图3为本申请实施例提供的本征模态函数筛选和组合流程图。具体实施方式此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本申请的实施例，并与说明书一起用于解释本申请的原理。为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，对于本领域普通技术人员而言，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。实施例1：图1为本申请实施例提供的一种中低纬电网中地磁感应电流的识别方法流程图，该方法包括如下所示的步骤：步骤一、对中低纬电网的样本数据进行聚类经验模态分解，获得第一本征模态函数组，所述第一本征模态函数组包括不同时间尺度的本征模态函数。可选的，所述样本数据包括：中低纬电网的中性点电流(NeutralPointCurrent，NPC)信号，和与所述NPC信号对应的太阳风能量耦合函数(Ein)。可选的，所述太阳风能量耦合函数根据太阳风基础参数计算得到。该太阳风相关基础参数是预设的，可以从美国航空航天局公布的开源数据库WIND卫星数据网址下载获得。聚类经验模态分解能使复杂信号分解为有限个本征模态函数(IntrinsicModeFunction，IMF)，所分解出来的各IMF分量包含了原样本数据的不同时间尺度的局部特征信号，具有自适应性。聚类经验模态分解能使非平稳数据进行平稳化处理，适合于分析非线性、非平稳信号序列。本实施例中以NPC信号(即NPC(t))为例，对聚类经验模态分解的流程进行说明，具体如图2所示。(1)在该NPC信号中添加白噪声后，提取该NPC信号的极值点，包括极大值点和极小值点。(2)采用样条函数分别拟合该NPC信号的上包络线和和下包络线，其中，上包络线为Emax(t)；下包络线为Emin(t)。(3)求上包络线和下包络线的平均值，即平均值Mij(t)＝[Emax(t)+Emin(t)]/2。(4)用NPC信号减去包络平均值获得差值Hij(t)，即差值Hij(t)＝NPC(t)-Mij(t)。(5)判断Hij(t)是否满足IMF条件。其中，IMF条件为：a.极值点(极大值或极小值)数目和过零点数目相等或最多相差一个。b.上包络线和下包络线的平均值为零。(6)如果Hij(t)满足IMF条件，则重复步骤(1)-步骤(5)X次，取Hij(t)的平均值作为第i个IMF，即Ci(t)＝mean(Hij(t))。随后执行步骤(8)。(7)如果Hij(t)不满足IMF条件，则重复步骤(1)-步骤(5)，直至Hij(t)满足IMF条件。随后执行步骤(6)。(8)该NPC信号减去Ci(t)得残差Ri(本文档来自技高网...

【技术保护点】
一种中低纬电网中地磁感应电流的识别方法，其特征在于，该方法包括：对中低纬电网的样本数据进行聚类经验模态分解，获得第一本征模态函数组，所述第一本征模态函数组包括不同时间尺度的本征模态函数；对所述第一本征模态函数组中的本征模态函数进行筛选和组合，获得第二本征模态函数组；对所述第二本征模态函数组进行希尔伯特变换，得到所述第二本征模态函数组的时频谱；所述时频谱为时间和频率的能量函数；根据所述时频谱，识别所述样本数据中是否具有地磁感应电流。

【技术特征摘要】
1.一种中低纬电网中地磁感应电流的识别方法，其特征在于，该方法包括：对中低纬电网的样本数据进行聚类经验模态分解，获得第一本征模态函数组，所述第一本征模态函数组包括不同时间尺度的本征模态函数；对所述第一本征模态函数组中的本征模态函数进行筛选和组合，获得第二本征模态函数组；对所述第二本征模态函数组进行希尔伯特变换，得到所述第二本征模态函数组的时频谱；所述时频谱为时间和频率的能量函数；根据所述时频谱，识别所述样本数据中是否具有地磁感应电流。2.根据权利要求1所述的方法，其特征在于，所述样本数据包括：中低纬电网的中性点电流信号，和与所述中性点电流信号对应的太阳风能量耦合函数。3.根据权利要求2所述的方法，其特征在于，所述对所述第一本征模态函数组中的本征模态函数进行筛选和组合，获得第二本征模态函数组，包括：将所述第一本征模态函数组中的本征模态函数随机组合，获得至少一个本征模态函数组；将所述至少一个本征模态函数...

【专利技术属性】
技术研发人员：杨家全，刘瑾，
申请(专利权)人：云南电网有限责任公司电力科学研究院，
类型：发明
国别省市：云南,53