光伏发电系统智能数据清洗方法与装置与流程

1.本发明属于光伏发电技术领域，具体涉及一种光伏发电系统智能数据清洗方法与装置。


背景技术：

2.在世界范围内能源环境危机下，作为化石燃料的一大消耗产业，电力行业的能源转型迫在眉睫。需探究各可再生能源的特性和可用性，广泛发展基于可再生新能源的发电技术实现发电领域的化石能源替代。基于此，风能、太阳能、地热能、潮汐能等可再生新能源成为世界各国关注的焦点。其中，太阳能能量巨大，几乎堪称“取之不尽，用之不竭”，且基本不受地域的限制，在地球绝大部分地方都可实现太阳能的开发利用。除此之外，光伏发电也是产生最早，发展最为成熟的新能源发电技术之一。虽然作为一种易于获取的可再生新能源，太阳能具有极大的环境友好性，然而其易受环境影响也成为其进一步发展过程中的一大困难和挑战。因此，需综合考虑光伏系统的建模和控制优化问题以促进电网的安全稳定。
3.当今世界，大数据技术飞速发展，而光伏发电系统的控制优化与电网电力调度都离不开对数据间特征和联系的挖掘与分析。然而，实际的光伏发电过程往往受到外界环境变化和光伏阵列局部遮挡等的影响，其运行数据中往往受到噪声等的干扰出现数据无效或缺失的情况。与此同时，不同运行工况下的光伏系统的数据一般也具有不同的特性。因此，通过设计合理、性能优越的数据清洗技术实现光伏系统的数据处理和优化是实现其稳定并网的关键一步。
4.数据清洗技术涉及到数据聚类、数据插值以及数据模型辨识等内容。以往的数据清洗技术中需通过不同的算法解决上述问题。随着群智能优化算法的发展，其在建模和控制等多个领域均取得了成功的应用。因此，其具有同时实现上述功能的开发潜力，并且可以从极大程度上提高光伏发电系统数据清洗过程中的智能化和灵活性，对其余复杂工业系统的数据清洗技术也具有极大适用性。
5.因此，基于上述问题，本发明提出一种光伏发电系统智能数据清洗方法与装置。


技术实现要素：

6.本发明旨在至少解决现有技术中存在的技术问题之一，提供一种光伏发电系统智能数据清洗方法与装置。
7.本发明的一方面，提供一种光伏发电系统智能数据清洗方法，包括下述具体步骤：
8.采集光伏发电系统的运行数据并对其进行变量选择；
9.基于所述运行数据，对igwo算法进行原理分析；
10.根据所述igwo算法对所述光伏发电系统的运行数据进行聚类；
11.对所得的每个聚类中的运行数据进行缺失值插补和无效值修正，以完成对所述运行数据清洗过程。
12.可选的，所述采集光伏发电系统的运行数据并对其进行变量选择，包括：
13.以采样步长t采集光伏系统在不同季节、月份、天数和时间点下的n组运行数据{x(1),x(2),
…
,x(n)}，每组运行数据的维度为m，代表光伏系统的m个数据变量，x(t)＝{x1(t),x2(t),
…
,xm(t)}，t表示当前采样点；
14.从m个数据变量中选出m个最终选用的数据变量。
15.可选的，所述基于所述运行数据，对igwo算法进行原理分析，包括：
16.对灰狼种群的混沌镜像进行初始化；
17.引入参数的余弦更新；
18.局部最优的判定及相应的变异操作。
19.可选的，所述对灰狼种群的混沌镜像进行初始化，包括：
20.在igwo的混沌镜像初始化过程中引入混沌序列，通过正弦公式迭代生成混沌矢量：
21.ch1＝rand(0,1)
22.ch
k+1
＝sin(πchk)(i＝1,2,
…
,n
p
；j＝1,2,
…
,d)
[0023][0024]
其中：d是搜索维度，chk代表第k时刻的混沌矢量且ch1为混沌矢量的初值，d
ij
指第i个灰狼在第j个搜索维度的位置向量， ch
kj
代表第k时刻的混沌矢量在第j个搜索维度的分量，和分别表示第j个搜索维度的上下界；
[0025]
基于所得灰狼种群执行镜像操作，并获得镜像解位置向量
[0026]
计算初始种群向量di＝{d
i,1
,d
i,2
,
…
,d
i,d
}和镜像种群的个体适应度并选择n
p
个具有最优适应度值的解最终形成初始种群向量其中d
′i＝{d
′
i,1
,d
′
i,2
,
…
,d
′
i,d
}(i＝1,
…
,n)。
[0027]
可选的，所述引入参数的余弦更新，包括：
[0028]
引入下式所示的余弦公式更新a：
[0029][0030]
可选的，所述局部最优的判定及相应的变异操作：
[0031]
引入阈值θ并通过下式判断优化过程是否陷入局部最优：
[0032][0033]
其中：fit(d
*
(k))为当前最优解的适应度值；
[0034]
当确认优化过程陷入早熟收敛后，通过如下变异操作扩大种群多样性，跳出当前的局部最优：
[0035]
[0036][0037]
j＝(1,2,
…
,d)
[0038]
其中，是第j维的更新最优解，是第j维的原最优解，w是选自[0,1]的随机数，η是分布指数，参数β1和β2通过下式获得：
[0039][0040][0041]
j＝(1,2,
…
,d)。
[0042]
可选的，所述igwo算法是基于gwo算法进行改进得到的。
[0043]
可选的，所述根据所述igwo算法对所述光伏发电系统的运行数据进行聚类，包括：
[0044]
定义聚类个数为nc，在数据对集合中随机选取nc个数据向量作为初始聚类中心；
[0045]
设置每个聚类中元素个数的阈值为th(th≥n/nc)，然后依次选取与聚类中心ci(i＝1,2,
…
,nc)距离最近的ni(i＝1,2,
…
,nc)个数据向量；
[0046]
若ni≥th，则自动转入下一个聚类，直到确定每个聚类的元素；
[0047]
以各聚类中心和本聚类中各元素间的平均距离最小为适应度函数f1，将聚类中心位置视为灰狼个体的位置，采用igwo算法进行最优聚类中心的获取，并将最优聚类中心到各元素的最远距离作为聚类半径
[0048]
可选的，所述对所得的每个聚类中的运行数据进行缺失值插补和无效值修正，以完成对所述运行数据清洗过程，包括：
[0049]
采用传递函数或状态空间模型的形式表示光伏系统每个聚类中数据变量间的关系模型，各聚类对应模型中的未知参数向量集合为其中代表第i(i＝1,2,
…
,nc)个模型所对应的未知参数向量且nb为未知参数个数；
[0050]
模型输出与实际输出的偏差作为适应度函数f2，将参数向量视为灰狼种群中个体的位置向量，通过改进型灰狼算法(igwo)算法进行模型参数的优化；
[0051]
重复模型参数的优化过程直到获取每个数据聚类所对应模型的最优参数向量；
[0052]
采用预测模型的输出填补数据序列中的缺失值，查找数据空间中与模型输出偏差极大的无效值，并用模型输出对该无效值进行替换或修正以完成数据清洗过程。
[0053]
本发明的另一方面，提供一种光伏发电系统智能数据清洗装置，包括采集模块、分析模块、聚类模块以及修正模块；其中，
[0054]
所述采集模块，用于采集光伏发电系统的运行数据并对其进行变量选择；
[0055]
所述分析模块，用于基于所述运行数据，对igwo算法进行原理分析；
[0056]
所述聚类模块，用于根据所述igwo算法对所述光伏发电系统的运行数据进行聚类；
[0057]
所述修正模块，用于对所得的每个聚类中的运行数据进行缺失值插补和无效值修正，以完成对所述运行数据清洗过程。
[0058]
本发明提供一种光伏发电系统智能数据清洗方法，包括下述具体步骤：采集光伏发电系统的运行数据并对其进行变量选择；基于所述运行数据，对igwo算法进行原理分析；根据所述igwo 算法对所述光伏发电系统的运行数据进行聚类；对所得的每个聚类中的运行数据进行缺失值插补和无效值修正，以完成对所述运行数据清洗过程。本发明结合igwo智能优化算法提供了一种光伏发电系统的数据清洗方法，通过数据采样、变量选择、数据聚类、数据插值和修正有效提高了光伏系统运行数据的质量，为后续系统的建模和控制优化设计奠定了基础。
附图说明
[0059]
图1为本发明一实施例的光伏发电系统智能数据清洗方法原理图；
[0060]
图2为本发明另一实施例的光伏发电系统智能数据清洗方法流程框图；
[0061]
图3为本发明另一实施例的改进型灰狼算法的流程图；
[0062]
图4为本发明另一实施例的光伏发电系统智能数据清洗装置示意图。
具体实施方式
[0063]
为使本领域技术人员更好地理解本发明的技术方案，下面结合附图和具体实施方式对本发明作进一步详细描述。显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于所描述的本发明的实施例，本领域普通技术人员在无需创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护范围。
[0064]
如图1和图2所示，本发明的一方面，提供一种光伏发电系统智能数据清洗方法s100，包括下述具体步骤s110～s140。
[0065]
s110、采集光伏发电系统的运行数据并对其进行变量选择。
[0066]
需要说明的是，光伏系统是用于收集、利用和转换太阳能的主要系统，光伏阵列作为光伏系统中的太阳能捕获装置，由许多光伏组件串并联连接而成。除此之外，光伏系统还包括dc-dc电压转换器、功率跟踪控制器和负载等组成部分，基于此，本实施例的步骤s110可具体化为：
[0067]
s1101、以采样步长t＝15min采集光伏系统在不同季节、月份、天数和时间点下的n＝10000组运行数据{x(1),x(2),
…
,x(n)}，每组运行数据的维度为m＝8，代表光伏系统的m个数据变量， x(t)＝{x1(t),x2(t),
…
,xm(t)}，t表示当前采样点。
[0068]
s1102、光伏系统的数据采样中一般涉及到时间因素、气象条件、外界环境变化、系统输出功率和内部其余运行变量，m的取值一般较大。而在光伏系统控制优化设计中为简便起见，往往仅会用到m个变量中与控制需求相关性较强的有限项。因此，通过主元分析等方法从m中选出m个最终选用的数据变量。
[0069]
基于s1采集所得光伏系统运行数据，在步骤s120中对所采用的改进型灰狼算法igwo原理做出具体说明。
[0070]
s120、基于运行数据，对igwo算法进行原理分析。
[0071]
具体的，s1201、为进一步突出igwo的优越性，首先对经典灰狼优化(gwo)算法进行介绍。
[0072]
灰狼算法(gwo算法)源于灰狼种群社会等级和群体狩猎行为的启发。在gwo中，将
整个种群按降序分为四个等级，并分别称之为α狼、β狼、δ狼和ω狼。其中α狼是种群所有重要事务的最高决策者，如狩猎行为、休息地点、活动时间等。β狼负责决策制定辅助、命令强化和信息反馈。处于种群第三级的δ狼则负责提交信息给α狼和β狼。δ狼群体由哨兵、年长者、猎人以及看护者组成。除α狼、β狼和δ狼之外，ω狼处于灰狼种群的最低级别，受到处于高级别的α狼、β狼和δ狼的驱使。因此，在狼群中低等级狼必须严格服从高等级狼。狼的捕猎过程包括跟踪阶段、包围阶段和攻击阶段，类似于优化问题求最优解的过程。
[0073]
具体的，s1201-1、设置灰狼种群规模大小为n
p
＝30，优化过程最大迭代次数为nm＝50，随机初始化种群为d(k)，其中k是当前迭代时间。
[0074]
s1201-2、定义最优解为α，次优解和再次优解分别为β和δ，根据实际需求定义适应度函数f，计算种群中每一个个体的适应度值，并根据适应度值选取α、β和δ。
[0075]
s1201-3、设灰狼个体与猎物间的距离d(k)为：
[0076][0077]
其中，d
p
(k)是猎物的位置，d(k)代表灰狼的位置，c是系数向量，r1是从[0,1]中随机选择的向量。
[0078]
s1201-4、通过下式更新灰狼个体的位置向量：
[0079][0080]
其中：a是系数向量，a是随着迭代次数的增加从2线性减少到0的收敛因子，r2也是从[0,1]中随机选择的向量。
[0081]
s1201-5、考虑到α狼、β狼和δ狼带有更多与猎物位置有关的信息，则当前搜索代理的位置根据α狼、β狼和δ狼的位置进行更新：
[0082][0083][0084]
其中：d
α
(k+1)、d
β
(k+1)和d
δ
(k+1)分别代表更新后的α狼、β狼和δ狼位置，a
α
、a
β
和a
δ
分别为α狼、β狼和δ狼对应的系数向量。α狼、β狼和δ狼与猎物间的距离d
α
、d
β
和d
δ
由下式计算得到：
[0085][0086]
s1201-6、计算灰狼个体位置更新后的适应度值，并根据贪婪选择选定新的α狼、β狼和δ狼。
[0087]
s1201-7、令k＝k+1，返回步骤s1201-3进行下一次迭代，直到满足终止条件，即k达到最大迭代次数或满足搜索精度的要求。
[0088]
s1202、然后基于gwo的主要原理，结合以下几大改进点对 igwo进行介绍，具体流
程参见图3。
[0089]
s1202-1、狼种群的混沌镜像初始化。考虑到混沌映射的确定性、遍历性和随机性，与单纯基于随机种群初始化的gwo相比，在igwo的混沌镜像初始化过程中引入混沌序列可增加种群多样性。通过正弦公式迭代生成混沌矢量：
[0090]
ch1＝rand(0,1)
[0091]
ch
k+1
＝sin(πchk) (i＝1,2,
…
,n
p
；j＝1,2,
…
,d)
[0092][0093]
其中，d是搜索维度，chk代表第k时刻的混沌矢量且ch1为混沌矢量的初值，d
ij
指第i个灰狼在第j个搜索维度的位置向量， ch
kj
代表第k时刻的混沌矢量在第j个搜索维度的分量，和分别表示第j个搜索维度的上下界。
[0094]
然后，基于所得灰狼种群执行镜像操作，并获得镜像解位置向量如果镜像位置向量的任何维度超出上下界，则引入正则化和圆弧正则化操作更新向量。
[0095]
计算初始种群向量di＝{d
i,1
,d
i,2
,
…
,d
i,d
}和镜像种群的个体适应度并选择n
p
个具有最优适应度值的解最终形成初始种群向量其中d
′i＝{d
′
i,1
,d
′
i,2
,
…
,d
′
i,d
}(i＝1,
…
,n)。
[0096]
s1202-2、参数的余弦更新。参数更新策略是优化算法搜索过程中的一个重要因素。在gwo中的所有参数中，a的值决定了其他参数值的变化。a的值越大，算法的全局搜索能力越强，反之，算法的局部搜索能力越强。对于高维复杂函数的优化，传统gwo 的全局搜索能力会随着a从2到0的线性下降而下降，难以避免早熟问题的出现。因此，引入下式所示的余弦公式更新a。
[0097][0098]
由上式可知，a的值在优化过程的早期阶段缓慢减小，以便进行充分的全局探索。在优化过程的后期，a的值会因局部开发而迅速降低。
[0099]
s1202-3、局部最优的判定及相应的变异操作。陷入局部最优是大多数优化算法中不可避免的问题。此时，如何正确判断搜索过程是否陷入局部最优并采取相应的修正措施是需要考虑的重要问题。
[0100]
先引入阈值θ(本实施例中取值为0.0001)并通过下式判断优化过程是否陷入局部最优：
[0101][0102]
其中：fit(d
*
(k))为当前最优解的适应度值。
[0103]
根据式(8)，如果连续n代的历史最优适应度的变化率小于某个阈值θ，则表明优化过程陷入局部最优。在这个判断过程中，阈值θ的选择会产生较大影响，如果θ的值非常大，则过早地判断搜索陷入局部最优。反之，若θ的值太小，则种群多样性将明显降低。
[0104]
当确认优化过程陷入早熟收敛后，通过如下变异操作扩大种群多样性，跳出当前
的局部最优。
[0105][0106]
其中，是第j维的更新最优解，是第j维的原最优解，w是选自[0,1]的随机数。η是分布指数，β1和β2通过下式获得：
[0107][0108]
根据s120构建的igwo算法，在s130中实现光伏系统实际采样数据的聚类。
[0109]
s130、根据所述igwo算法对所述光伏发电系统的运行数据进行聚类。
[0110]
具体的，s1301、定义聚类个数为nc，在数据对集合中随机选取nc个数据向量作为初始聚类中心。
[0111]
s1302、设置每个聚类中元素个数的阈值为th(th≥n/nc)，然后依次选取与聚类中心ci(i＝1,2,
…
,nc)距离最近的ni(i＝1,2,
…
,nc) 个数据向量。
[0112]
s1303、若ni≥th，则自动转入下一个聚类，直到确定每个聚类的元素。
[0113]
s1304、以各聚类中心和本聚类中各元素间的平均距离最小为适应度函数f1，将聚类中心位置视为灰狼个体的位置，采用igwo 算法进行最优聚类中心的获取，进一步地，将最优聚类中心到各元素的最远距离作为聚类半径
[0114]
接下来在步骤s140中步骤s130中所得的每个聚类中的数据进行缺失值插补和无效值修正。主要通过igwo算法建立数据变量间的关系模型实现。
[0115]
s140、对所得的每个聚类中的运行数据进行缺失值插补和无效值修正，以完成对所述运行数据清洗过程。
[0116]
具体的，s1401、采用传递函数或状态空间模型的形式表示光伏系统每个聚类中数据变量间的关系模型，各聚类对应模型中的未知参数向量集合为其中代表第i(i＝1,2,
…
,nc)个模型所对应的未知参数向量且nb为未知参数个数。
[0117]
s1402、以模型输出与实际输出的偏差作为适应度函数f2，将参数向量视为灰狼种群中个体的位置向量，通过igwo算法进行模型参数的优化。
[0118]
s1403、重复s1402直到获取每个数据聚类所对应模型的最优参数向量；
[0119]
s1404、采用预测模型的输出填补数据序列中的缺失值，查找数据空间中与模型输出偏差极大的无效值，并用模型输出对该无效值进行替换或修正以完成数据清洗过程。
[0120]
为挖掘光伏发电系统不同运行情况下的数据特征并对其运行数据中的数据缺失及数据异常等问题进行处理，本发明旨在提供一种基于改进灰狼算法的光伏发电系统智能
数据清洗方法以去除光伏系统运行数据中的噪声及无效信息。该方法首先采集光伏系统季节、月份和天数等不同时间尺度和外界环境下的运行数据并对数据变量进行筛选，然后引入了一种改进的灰狼优化(igwo)算法，通过基于混沌镜像的种群初始化、基于余弦变化的参数更新和局部最优判断有效平衡了算法的搜索和开发能力，使参数收敛速度和精度都有明显的提升。通过igwo算法进行数据聚类和数据变量间的关系辨识，可同时实现数据空间的划分、数据插补和数据修正。有效提高了数据质量，为后续光伏发电系统的控制优化设计奠定了基础。
[0121]
如图4所示，本发明的另一方面，提供一种光伏发电系统智能数据清洗装置200，包括采集模块210、分析模块220、聚类模块230以及修正模块240；其中，所述采集模块210，用于采集光伏发电系统的运行数据并对其进行变量选择；所述分析模块220，用于基于所述运行数据，对igwo算法进行原理分析；所述聚类模块230，用于根据所述igwo算法对所述光伏发电系统的运行数据进行聚类；所述修正模块240，用于对所得的每个聚类中的运行数据进行缺失值插补和无效值修正，以完成对所述运行数据清洗过程。
[0122]
需要说明的是，本实施例的数据清洗装置所采用的具体方法参考前文记载，在此不再赘述。
[0123]
下面将以具体实施例进行说明光伏发电系统智能数据清洗方法，具体包括：
[0124]
请参阅图1，图1为本发明所提供的基于改进灰狼算法的光伏发电系统智能数据清洗方法原理图。本实施例中算法的实现及仿真测试均依托matlab软件平台，具体由以下4个步骤组成：
[0125]
s1：光伏发电系统的运行数据采集和变量选择；
[0126]
s2：igwo的原理分析及描述；
[0127]
s3：基于igwo的光伏发电系统数据聚类；
[0128]
s4：基于igwo的光伏发电系统数据插值和修正。
[0129]
光伏系统是用于收集、利用和转换太阳能的主要系统。光伏阵列作为光伏系统中的太阳能捕获装置，由许多光伏组件串并联连接而成。除此之外，光伏系统还包括dc-dc电压转换器、功率跟踪控制器和负载等组成部分。基于此，步骤s1可具体化为：
[0130]
s1.1：以采样步长t＝15min采集光伏系统在不同季节、月份、天数和时间点下的n＝10000组运行数据{x(1),x(2),
…
,x(n)}，每组运行数据的维度为m＝8，代表光伏系统的m个数据变量，即 x(t)＝{x1(t),x2(t),
…
,xm(t)}，t为当前采样点。
[0131]
s1.2：光伏系统的数据采样中一般涉及到时间因素、气象条件、外界环境变化、系统输出功率和内部其余运行变量，m的取值一般较大。而在光伏系统控制优化设计中为简便起见，往往仅会用到m个变量中与控制需求相关性较强的有限项。因此，通过主元分析等方法从m中选出m＝6个最终选用的数据变量。
[0132]
基于s1采集所得光伏系统运行数据，在步骤s2中对所采用的igwo原理做出具体说明。
[0133]
s2.1：为进一步突出igwo的优越性，首先对经典灰狼优化 (gwo)算法进行介绍。
[0134]
gwo算法源于灰狼种群社会等级和群体狩猎行为的启发。在 gwo中，将整个种群按降序分为四个等级，并分别称之为α狼、β狼、δ狼和ω狼。其中α狼是种群所有重要事务的最高决策者，如狩猎行为、休息地点、活动时间等。β狼负责决策制定辅助、命令强化和信息反
馈。处于种群第三级的δ狼则负责提交信息给α狼和β狼。δ狼群体由哨兵、年长者、猎人以及看护者组成。除α狼、β狼和δ狼之外，ω狼处于灰狼种群的最低级别，受到处于高级别的α狼、β狼和δ狼的驱使。因此，在狼群中低等级狼必须严格服从高等级狼。狼的捕猎过程包括跟踪阶段、包围阶段和攻击阶段，类似于优化问题求最优解的过程。
[0135]
s2.1.1：设置灰狼种群规模大小为n
p
＝30，优化过程最大迭代次数为nm＝50，随机初始化种群为d(k)，其中k是当前迭代时间。
[0136]
s2.1.2：记最优解为α，次优解和再次优解分别为β和δ，根据实际需求定义适应度函数f，计算种群中每一个个体的适应度值，并根据适应度值选取α、β和δ。
[0137]
s2.1.3：设灰狼个体与猎物间的距离d(k)为：
[0138][0139]
其中：d
p
(k)是猎物的位置，d(k)代表灰狼的位置，c是系数向量，r1是从[0,1]中随机选择的向量，即r1＝rand[0,1]。
[0140]
s2.1.4：通过下式更新灰狼个体的位置向量：
[0141][0142]
其中：a是系数向量，a是随着迭代次数的增加从2线性减少到0的收敛因子，r2也是从[0,1]中随机选择的向量，即r2＝rand[0,1]。
[0143]
s2.1.5：考虑到α狼、β狼和δ狼带有更多与猎物位置有关的信息，则当前搜索代理的位置根据α狼、β狼和δ狼的位置进行更新：
[0144][0145][0146]
其中：d
α
(k+1)、d
β
(k+1)和d
δ
(k+1)分别代表更新后的α狼、β狼和δ狼位置，a
α
、a
β
和a
δ
分别为α狼、β狼和δ狼对应的系数向量。α狼、β狼和δ狼与猎物间的距离d
α
、d
β
和d
δ
由下式计算得到：
[0147][0148]
s2.1.6：计算灰狼个体位置更新后的适应度值，并根据贪婪选择选定新的α狼、β狼和δ狼。
[0149]
s2.1.7：令k＝k+1，返回步骤s2.1.3进行下一次迭代，直到满足终止条件，即k达到最大迭代次数或满足搜索精度的要求。
[0150]
s2.2：然后基于gwo的主要原理，结合以下几大改进点对 igwo进行介绍，具体流程参见图3。
[0151]
s2.2.1：灰狼种群的混沌镜像初始化。考虑到混沌映射的确定性、遍历性和随机
性，与单纯基于随机种群初始化的gwo相比，在igwo的混沌镜像初始化过程中引入混沌序列可增加种群多样性。通过正弦公式迭代生成混沌矢量：
[0152]
ch1＝rand(0,1)
[0153]
ch
k+1
＝sin(πchk) (i＝1,2,
…
,n
p
；j＝1,2,
…
,d)
[0154][0155]
其中：d是搜索维度，本实施例中d＝m＝6，和分别表示第j个搜索维度的上下界。
[0156]
然后，基于所得灰狼种群执行镜像操作，并获得镜像解位置向量如果镜像位置向量的任何维度超出上下界，则引入正则化和圆弧正则化操作更新向量。
[0157]
计算初始种群向量di＝{d
i,1
,d
i,2
,
…
,d
i,d
}和镜像种群的个体适应度并选择n
p
个具有最优适应度值的解最终形成初始种群向量其中d
′i＝{d
′
i,1
,d
′
i,2
,
…
,d
′
i,d
}(i＝1,
…
,n)。
[0158]
s2.2.2：参数的余弦更新。参数更新策略是优化算法搜索过程中的一个重要因素。在gwo中的所有参数中，a的值决定了其他参数值的变化。a的值越大，算法的全局搜索能力越强，反之，算法的局部搜索能力越强。对于高维复杂函数的优化，传统gwo 的全局搜索能力会随着a从2到0的线性下降而下降，难以避免早熟问题的出现。因此，引入下式所示的余弦公式更新a。
[0159][0160]
由上式可知，a的值在优化过程的早期阶段缓慢减小，以便进行充分的全局探索。在优化过程的后期，a的值会因局部开发而迅速降低。
[0161]
s2.2.3：局部最优的判定及相应的变异操作。陷入局部最优是大多数优化算法中不可避免的问题。此时，如何正确判断搜索过程是否陷入局部最优并采取相应的修正措施是需要考虑的重要问题。
[0162]
先引入阈值θ(本实施例中取值为0.0001)并通过下式判断优化过程是否陷入局部最优：
[0163][0164]
其中：fit(d
*
(k))为当前最优解的适应度值。
[0165]
根据式(8)，如果连续n＝10代的历史最优适应度的变化率小于某个阈值θ，则表明优化过程陷入局部最优。
[0166]
当确认优化过程陷入早熟收敛后，通过如下变异操作扩大种群多样性，跳出当前的局部最优。
[0167][0168]
其中，是第j维的更新最优解，是第j维的原最优解，w是选自[0,1]的随机数，即w＝rand[0,1]。η是分布指数，β1和β2通过下式获得：
[0169][0170]
根据s2构建的igwo算法，在s3中实现光伏系统实际采样数据的聚类。
[0171]
s3.1：定义聚类个数为nc＝5，在数据对集合中随机选取nc个数据向量作为初始聚类中心。
[0172]
s3.2：设置每个聚类中元素个数的阈值为th＝4000(th≥n/nc)，然后依次选取与聚类中心ci(i＝1,2,
…
,nc)距离最近的ni(i＝1,2,
…
, nc)个数据向量。
[0173]
s3.3：若ni≥th，则自动转入下一个聚类，直到确定每个聚类的元素。
[0174]
s3.4：以各聚类中心和本聚类中各元素间的平均距离最小为适应度函数f1，将聚类中心位置视为灰狼个体的位置，采用igwo 算法进行最优聚类中心的获取，进一步地，将最优聚类中心到各元素的最远距离作为聚类半径
[0175]
接下来在步骤s4中步骤s3中所得的每个聚类中的数据进行缺失值插补和无效值修正。主要通过igwo算法建立数据变量间的关系模型实现。
[0176]
s4.1：用传递函数或状态空间模型的形式表示光伏系统每个聚类中数据变量间的关系模型，各聚类对应模型中的未知参数向量集合为其中代表第i(i＝1,2,
…
,nc)个模型所对应的未知参数向量且nb为未知参数个数且本实施例中 nb＝27。
[0177]
s4.2：以模型输出与实际输出的偏差为此时的适应度函数f2，将参数向量视为灰狼种群中个体的位置向量，然后通过igwo进行模型参数的优化。
[0178]
s4.3：重复s4.2直到获取每个数据聚类所对应模型的最优参数向量即完成不同数据变量间的关系表征过程。
[0179]
s4.4：首先，用预测模型的输出填补数据序列中的缺失值。然后，查找数据空间中与模型输出偏差极大的无效值，并用模型输出对该无效值进行替换或修正以完成数据清洗过程。
[0180]
本发明提供一种光伏发电系统智能数据清洗方法与装置，与现有技术相比具有以下有益效果：
[0181]
第一、本发明结合igwo智能优化算法提供了一种光伏发电系统的数据清洗方法，通过数据采样、变量选择、数据聚类、数据插值和修正有效提高了光伏系统运行数据的质
量，为后续系统的建模和控制优化设计奠定了基础。
[0182]
第二、本发明所用的igwo算法通过基于混沌镜像的种群初始化、基于余弦变化的参数更新和局部最优判断有效平衡了算法的搜索和开发能力，使优化速度和精度都有明显的提升。
[0183]
第三、本发明同时将igwo算法用于数据聚类和数据变量间的关系辨识，可同时实现数据空间的划分、数据插补和数据修正。极大地体现了igwo算法的可扩展性和应用灵活性。
[0184]
可以理解的是，以上实施方式仅仅是为了说明本发明的原理而采用的示例性实施方式，然而本发明并不局限于此。对于本领域内的普通技术人员而言，在不脱离本发明的情况下，可以做出各种变型和改进，这些变型和改进也视为本发明的保护范围。

技术特征：
1.一种光伏发电系统智能数据清洗方法，其特征在于，包括下述具体步骤：采集光伏发电系统的运行数据并对其进行变量选择；基于所述运行数据，对igwo算法进行原理分析；根据所述igwo算法对所述光伏发电系统的运行数据进行聚类；对所得的每个聚类中的运行数据进行缺失值插补和无效值修正，以完成对所述运行数据的清洗过程。2.根据权利要求1所述的方法，其特征在于，所述采集光伏发电系统的运行数据并对其进行变量选择，包括：以采样步长t采集光伏系统在不同季节、月份、天数和时间点下的n组运行数据{x(1),x(2),
…
,x(n)}，每组运行数据的维度为m，代表光伏系统的m个数据变量，x(t)＝{x1(t),x2(t),
…
,x
m
(t)}，t表示当前采样点；从m个数据变量中选出m个最终选用的数据变量。3.根据权利要求1所述的方法，其特征在于，所述基于所述运行数据，对igwo算法进行原理分析，包括：对灰狼种群的混沌镜像进行初始化；引入参数的余弦更新；局部最优的判定及相应的变异操作。4.根据权利要求3所述的方法，其特征在于，所述对灰狼种群的混沌镜像进行初始化，包括：在igwo的混沌镜像初始化过程中引入混沌序列，通过正弦公式迭代生成混沌矢量：ch1＝rand(0,1)ch
k+1
＝sin(πch
k
)(i＝1,2,
…
,n
p
；j＝1,2,
…
,d)其中：d是搜索维度，ch
k
代表第k时刻的混沌矢量且ch1为混沌矢量的初值，d
ij
指第i个灰狼在第j个搜索维度的位置向量，ch
kj
代表第k时刻的混沌矢量在第j个搜索维度的分量，和分别表示第j个搜索维度的上下界；基于所得灰狼种群执行镜像操作，并获得镜像解位置向量计算初始种群向量d
i
＝{d
i,1
,d
i,2
,
…
,d
i,d
}和镜像种群的个体适应度并选择n
p
个具有最优适应度值的解最终形成初始种群向量其中d
′
i
＝{d
′
i,1
,d
′
i,2
,
…
,d
′
i,d
}(i＝1,
…
,n)。5.根据权利要求3所述的方法，其特征在于，所述引入参数的余弦更新，包括：引入下式所示的余弦公式更新a：6.根据权利要求3所述的方法，其特征在于，所述局部最优的判定及相应的变异操作：引入阈值θ并通过下式判断优化过程是否陷入局部最优：
其中：fit(d
*
(k))为当前最优解的适应度值；当确认优化过程陷入早熟收敛后，通过如下变异操作扩大种群多样性，跳出当前的局部最优：部最优：j＝(1,2,
…
,d)其中，是第j维的更新最优解，是第j维的原最优解，w是选自[0,1]的随机数，η是分布指数，参数β1和β2通过下式获得：通过下式获得：(j＝1,2,
…
,d)。7.根据权利要求1至6任一项所述的方法，其特征在于，所述igwo算法是基于gwo算法进行改进得到的。8.根据权利要求1至6任一项所述的方法，其特征在于，所述根据所述igwo算法对所述光伏发电系统的运行数据进行聚类，包括：定义聚类个数为n
c
，在数据对集合中随机选取n
c
个数据向量作为初始聚类中心；设置每个聚类中元素个数的阈值为th(th≥n/n
c
)，然后依次选取与聚类中心c
i
(i＝1,2,
…
,n
c
)距离最近的n
i
(i＝1,2,
…
,n
c
)个数据向量；若n
i
≥th，则自动转入下一个聚类，直到确定每个聚类的元素；以各聚类中心和本聚类中各元素间的平均距离最小为适应度函数f1，将聚类中心位置视为灰狼个体的位置，采用igwo算法进行最优聚类中心的获取，并将最优聚类中心到各元素的最远距离作为聚类半径9.根据权利要求1至6任一项所述的方法，其特征在于，所述对所得的每个聚类中的运行数据进行缺失值插补和无效值修正，以完成对所述运行数据清洗过程，包括：采用传递函数或状态空间模型的形式表示光伏系统每个聚类中数据变量间的关系模型，各聚类对应模型中的未知参数向量集合为其中代表第i(i＝1,2,
…
,n
c
)个模型所对应的未知参数向量且n
b
为未知参数个数；模型输出与实际输出的偏差作为适应度函数f2，将参数向量视为灰狼种群中个体的位置向量，通过igwo算法进行模型参数的优化；重复模型参数的优化过程直到获取每个数据聚类所对应模型的最优参数向量；采用预测模型的输出填补数据序列中的缺失值，查找数据空间中与模型输出偏差极大的无效值，并用模型输出对该无效值进行替换或修正以完成数据清洗过程。
10.一种光伏发电系统智能数据清洗装置，其特征在于，包括采集模块、分析模块、聚类模块以及修正模块；其中，所述采集模块，用于采集光伏发电系统的运行数据并对其进行变量选择；所述分析模块，用于基于所述运行数据，对igwo算法进行原理分析；所述聚类模块，用于根据所述igwo算法对所述光伏发电系统的运行数据进行聚类；所述修正模块，用于对所得的每个聚类中的运行数据进行缺失值插补和无效值修正，以完成对所述运行数据的清洗过程。

技术总结
本发明提供一种光伏发电系统智能数据清洗方法与装置，属于光伏发电技术领域。其中，本发明的方法包括下述具体步骤：采集光伏发电系统的运行数据并对其进行变量选择；基于所述运行数据，对改进型灰狼算法(IGWO)算法进行原理分析；根据所述IGWO算法对所述光伏发电系统的运行数据进行聚类；对所得的每个聚类中的运行数据进行缺失值插补和无效值修正，以完成对所述运行数据清洗过程。本发明结合IGWO智能优化算法提供了一种光伏发电系统的数据清洗方法，通过数据采样、变量选择、数据聚类、数据插值和修正有效提高了光伏系统运行数据的质量，为后续系统的建模和控制优化设计奠定了基础。续系统的建模和控制优化设计奠定了基础。续系统的建模和控制优化设计奠定了基础。


技术研发人员：曾凡春 崔源 李丹阳 王军 陈婷婷 李涛 陈岩磊 张澈 王传鑫 田长凤
受保护的技术使用者：北京华能新锐控制技术有限公司
技术研发日：2021.11.24
技术公布日：2022/3/7