CN109617526A - 一种基于小波多分辨分析和svm的光伏发电阵列故障诊断和分类的方法 - Google Patents
一种基于小波多分辨分析和svm的光伏发电阵列故障诊断和分类的方法 Download PDFInfo
- Publication number
- CN109617526A CN109617526A CN201811561729.0A CN201811561729A CN109617526A CN 109617526 A CN109617526 A CN 109617526A CN 201811561729 A CN201811561729 A CN 201811561729A CN 109617526 A CN109617526 A CN 109617526A
- Authority
- CN
- China
- Prior art keywords
- signal
- array
- classification
- fault diagnosis
- sample
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000003745 diagnosis Methods 0.000 title claims abstract description 34
- 238000004458 analytical method Methods 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000010248 power generation Methods 0.000 title claims abstract description 23
- 239000013598 vector Substances 0.000 claims abstract description 39
- 238000012549 training Methods 0.000 claims abstract description 32
- 230000008859 change Effects 0.000 claims abstract description 28
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 13
- 238000012360 testing method Methods 0.000 claims abstract description 13
- 238000012706 support-vector machine Methods 0.000 claims abstract description 12
- 230000009466 transformation Effects 0.000 claims abstract description 10
- 238000005286 illumination Methods 0.000 claims abstract description 7
- 230000007257 malfunction Effects 0.000 claims description 7
- 238000010606 normalization Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 4
- 230000005611 electricity Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 230000006870 function Effects 0.000 description 28
- 239000000306 component Substances 0.000 description 14
- 238000005457 optimization Methods 0.000 description 10
- 238000001514 detection method Methods 0.000 description 6
- 241001269238 Data Species 0.000 description 5
- 230000009977 dual effect Effects 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000002310 reflectometry Methods 0.000 description 3
- 238000001931 thermography Methods 0.000 description 3
- 230000004888 barrier function Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000011664 signaling Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- WSGCGMGMFSSTNK-UHFFFAOYSA-M 1-methyl-4-phenyl-1-propan-2-ylpiperidin-1-ium;iodide Chemical compound [I-].C1C[N+](C(C)C)(C)CCC1C1=CC=CC=C1 WSGCGMGMFSSTNK-UHFFFAOYSA-M 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 238000003066 decision tree Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02S—GENERATION OF ELECTRIC POWER BY CONVERSION OF INFRARED RADIATION, VISIBLE LIGHT OR ULTRAVIOLET LIGHT, e.g. USING PHOTOVOLTAIC [PV] MODULES
- H02S50/00—Monitoring or testing of PV systems, e.g. load balancing or fault identification
- H02S50/10—Testing of PV devices, e.g. of PV modules or single PV cells
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
Landscapes
- Testing Of Individual Semiconductor Devices (AREA)
Abstract
本发明涉及一种基于小波多分辨分析和SVM的光伏发电阵列故障诊断和分类的方法，首先采集光伏阵列不同工作状态下不同温照度的电流和电压样本信号；接着进行归一化处理得到电流变化率、电压变化率、功率变化率和电导变化率；然后对得到的四个样本信号进行滑动取窗口信号，并进行小波变换的多分辨率分析；计算四个窗口信号第N层分解的高频信号的2范数，进而构成维度为4的特征向量；然后由不同工作状态的多组样本信号得到多个特征值向量，并分为训练数据与测试数据；最后设置支持向量机SVM的参数，从训练数据中训练故障诊断模型，用测试数据检测该模型的分类的准确性。本发明环境适用性强，有较强的创新性，训练的模型能实现较高精度的故障诊断和分类。
Description
技术领域
本发明涉及光伏发电故障诊断和分类技术领域，特别是一种基于小波多分辨分析和SVM的光伏发电阵列故障诊断和分类的方法。
背景技术
发展新型能源是当今解决化石能源枯竭和环境污染等问题的重要战略，而太阳能因其清洁，能源丰富等优势成为发展新能源的首选。光伏发电是目前利用太阳能的主要方式，而采集太阳能的光伏发电阵列作为光伏发电***的核心部件通常工作在复杂多变的户外环境中，易受风霜雨雪等恶劣因素的影响而产生诸如短路，开路，阴影等故障。故障的产生会降低光伏发电效率，造成意外的能源浪费和经济损失，更有甚者会留下安全隐患，引起火灾，危害人身安全。所以对光伏***的工作状态进行监控，对出现的故障及时地检测和分类并发出警告，就能够有效地进行维护，降低因光伏阵列故障带来的能量损失，防止出现安全事故，具有良好的经济和社会效益。
国内外学者们已经先后提出了一系列的故障诊断方法，实现对不同的故障类型的检测和定位。典型的故障检测方法有对地电容检测法，时域反射分析法，红外热成像等。此外，随着人工智能的快速发展，学者们提出了运用支持向量机，神经网络，决策树等基于机器学习算法的故障诊断方案。对地电容检测法是根据对光伏组串的对地电容的测量来判断其是否发生断路并定位故障。时域反射分析法是向光伏组串注入一个脉冲，对返回信号的形状和延迟时间进行分析以此判断光伏组串中是否存在故障。对地电容测量法和时域反射分析法都需要进行离线检测，缺乏实时性，这样会耗费大量的人力和财力。而工作于正常状态和故障状态下的太阳电池会有存在明显的温差，故还可以采用红外热成像分析法进行故障诊断。红外热成像分析法虽然可以高效的进行故障诊断，但必须配备大量的红外摄像仪，经济效益差，难以得到推广。目前应用最多的是基于机器学习算法的故障诊断和分类方法。该方法具有较强的自学习能力，鲁棒性强，准确率高，已经成为研究的热点。目前大多数的机器学习算法是基于光伏阵列最大功率点电流IMPP,最大功率点电压VMPP,短路电流ISC，开路电压VOC和温照度G，环境温度T等变量作为特征进行训练学习，探索新的特征向量是研究光伏阵列故障诊断的重要研究方向。精确的故障诊断模型一般要求多维的数据，模型训练时间长，环境的不断变换也给模型的准确性带来了挑战。
发明内容
有鉴于此，本发明的目的是提出一种基于小波多分辨分析和SVM的光伏发电阵列故障诊断和分类的方法，环境适用性强，有较强的创新性，训练的模型能实现较高精度的故障诊断和分类。
本发明采用以下方案实现：一种基于小波多分辨分析和SVM的光伏发电阵列故障诊断和分类的方法，包括以下步骤：
步骤S1：采集光伏阵列不同工作状态下不同温照度的电流和电压样本信号；
步骤S2：对采集的电流和电压信号进行归一化处理，并得到电流变化率、电压变化率、功率变化率和电导变化率四个样本信号；
步骤S3：设置窗口大小为L的滑动窗对步骤S2得到的四个样本信号进行滑动取窗口信号，对取得的四种窗口信号进行小波变换的多分辨率分析；
步骤S4：计算四个窗口信号第N层分解的高频信号的2范数；
步骤S5：取四个样本信号在故障发生时的窗口信号的第N层分解的高频信号的2范数作为特征值，构成维度为4的特征向量；
步骤S6：由不同工作状态的多组样本信号得到多个特征值向量，其中将m个特征向量作为训练数据，将n个特征向量作为测试数据；
步骤S7：设置支持向量机SVM的参数，从训练数据中训练故障诊断模型，用测试数据检测该模型的分类的准确性。
进一步地，步骤S1中，光伏阵列的不同工作状态包括：正常、单组串1个组件短路、单组串一个组件开路以及单组串2个组件短路；其中，单组串1个组件短路、单组串一个组件开路以及单组串2个组件短路为故障状态，故障状态的样本信号包含了光伏阵列由正常到故障再到故障稳定的过程，捕获了故障发生时阵列电流和电压的变化特性。其中，电流和电压样本信号的采集时间窗口t＝10s，采样频率为200Hz。
进一步地，步骤S2包括以下步骤：
步骤S21：对采集的电流和电压信号进行归一化处理：将阵列电流除以短路电流，将阵列电压除以开路电压；经过归一化处理能够克服不同光照度和不同温度的影响，归一化的公式如下：
iPV(t)＝IPV(t)/ISC(t)；
vPV(t)＝VPV(t)/VOC(t)
式中，IPV(t)表示采集的阵列电流样本信号，ISC(t)表示阵列短路电流信号，iPV(t)表示归一化后的阵列电流样本信号,VPV(t)表示采集的阵列电压样本信号，VOC(t)表示阵列开路信号,vPV(t)表示归一化后的阵列电压样本信号；经过归一化处理之后的电流和电压信号反映了不同工作状态下阵列信号的变化趋势，突出了变化特性；
步骤S22：采用下式得到电流变化率、电压变化率、功率变化率和电导变化率四个样本信号：
di(t)＝[iPV(n+1)-iPV(n)]/t；
dv(t)＝[vPV(n+1)-vPV(n)]/t；
dp(t)＝[iPV(n+1)v(n+1)-iPV(n)v(n)]/t；
dg(t)＝[iPV(n+1)/vPV(n+1)-iPV(n)/iPV(n)]/t；
式中，di(t)、dv(t)、dp(t)、dg(t)分别描述了光伏阵列不同工作状态的电流变化率，电压变化率，功率变化率和电导变化率的信号。
进一步地，步骤S3中，所述窗口大小L为20，采用DB小波对取得的四种窗口信号进行6层多分辨率分解的小波变换；所述多分辨率分析为用高通和低通滤波器将信号分解成多个级别的高频细节和近似值。
较佳的，小波变换是一种应用广泛的数字信号处理方法，适用于分析非周期信号，可以同时对信号进行时间和频率分析。信号x(t)离散化为x(n),该信号的DWT可以由下式求解：
式中，ψ(t)是母小波，n是信号x(n)的点数，a＝a0 j是尺度因子，b＝kb0a0 j平移因子。a0＞0且为常数，b0≠0，j,k∈Z。
多分辨率分析是用高通和低通滤波器将信号分解成多个级别的高频细节和近似值。时域信号f(t)可以由尺度函数Φ(t)和小波函数Ψ(t)表示,如下式所示：
式中，dj(k)是各级分解的细节系数，表达高频信息，而aN(k)最后一级N级的近似系数，表达低频信息。
进一步地，步骤S4中，计算四个窗口信号第N层分解的高频信号的2范数采用下式：
式中，djn表示，N表示，dj表示第j层的高频系数,fi表示第i个信号提取的特征。
进一步地，步骤S5中，所述维度为4的特征向量为f＝[f1,f2,f3,f4]。
进一步地，步骤S7中，所述设置支持向量机SVM的参数具体为：将SVM的惩罚因子C设置为1000，将两个不同类别的支持向量到超平面的距离之和γ设置为5。具体为：
支持向量机的基本思想是通过线性可分的训练样本集找到一个最优分类超平面以实现不同类别的样本数据的划分。对于非线性的分类问题，可以将样本从原始空间映射到一个线性可分的高维特征空间，从而找到一个合适的最优分类超平面。
给定训练样本数据集，D＝{xi,yi}，i＝1,2,3,...,m,yi∈{-1,1}，其中xi为样本数据，m为训练样本总数，d为样本空间的维度，yi为样本相应的标签。它们能够被一个最优超平面所分离，该超平面可以表示为wTx+b＝0，其中w∈Rd，为法向量，决定了超平面的方向；b∈R，为位移项，决定了超平面与原点间的距离，分类的阈值。
假设超平面能将训练样本正确分类，则对于{xi,yi}∈D，如果yi＝+1，则有wTx+b＞0；如果yi＝-1，则有wTx+b＜0。令：
则距离超平面最近的若干训练样本点使得上式的等号成立，它们被称之为支持向量(Support vectors，SVs)，则两个不同类别的支持向量到超平面的距离之和为该距离被称为分类间隔。如果要使γ最大，则要求||w||2最小，同时要求分类面对所有的样本正确分类，需满足
yi(wTxi+b)≥1,i＝1,2,3,...,l
因此，求最优分类超平面问题，可以转化为一个二次规划问题，优化目标可以写为
s.t.yi(wTxi+b)≥1,i＝1,2,3,...,l
对于在样本空间线性可分的训练样本，其可以被最优分类超平面进行划分。然后，在现实任务中，常常存在线性不可分的情况，此时在训练样本中存在部分样本数据不满足公式(1)，存在一定分类的误差。因此，通过引入松弛变量ξi(ξi≥0)来解决此问题。因此，公式(1)可写为
yi(wTxi+b)≥1-ξi,i＝1,2,3,...,l
在最大化间隔的同时，要求不满足约束的样本要尽可能少。因此，优化目标函数可以改写为
其中，称为惩罚项，C为惩罚因子。因此，可以得到线性不可分时的最优分类面，称为广义分类超平面，可以表示为下面的优化问题
s.t.yi(wTxi+b)≥1-ξi
ξi≥0,i＝1,2,3,...,l
求解该优化问题可以借助拉格朗日乘子法，得到它的对偶问题如下式：
0≤αi≤C,i,j＝1,2,3,...,l
根据Karush–Kuhn–Tucker(KKT)条件，可以得到αi(yi(wTxi+b)-1)＝0。如果αi＞0，则所对应的样本点位于最大间隔边界上，就是支持向量。则可以通过求解出w，并根据yi(wTxi+b)-1＝0求解出b，其中xi为支持向量，n为支持向量的数量。在确定w和b之后，得到分类决策函数如下所示
其次，针对非线性分类问题，SVM采用核函数将样本从原始空间映射到更高维的特征空间，使得样本在这个特征空间内线性可分。令φ(x)表示将x映射后的特征向量，则在高维空间中超平面所对应的优化目标可表示为
s.t.yi(wTφ(xi)+b)≥1-ξi
ξi≥0,i＝1,2,3,...,l
对应的对偶问题为
0≤αi≤C,i,j＝1,2,3,...,l
为避开计算样本xi和xj在高维空间中的内积运算，构造核函数K(·,·)，xi和xj在特征空间的内积转换为在原始样本空间中通过该函数计算的结果。K(·,·)表示如下
K(xi,xj)＝φ(xi)Tφ(xj)
则将前一个公式变为
0≤αi≤C,i,j＝1,2,3,...,l
此时分类决策函数变为
常用的核函数主要有多项式核函数、径向基(Radial based kernel，RBF)核函数、双曲正切(Sigmoid)核函数等。本发明中采用RBF核函数对SVM进行应用。
RBF核函数表示为
K(xi,x)＝exp(-γ||xi-x||2)
支持向量机的分类效果将很大程度取决于C和γ的选择，本发明中将C设置为1000，γ设置为10。运用设置好参数SVM，从构造的训练数据中训练出故障分类器，实现光伏发电阵列的故障诊断和分类。
与现有技术相比，本发明有以下有益效果：
1、本发明只采集了不同温照度下的正常和故障的电流和电压信号，由此得到电流变化率，电压变化率，功率变化率和电导变化率四个样本信号，突出了故障发生时的特性。不需要其他的数据特征，且能在不影响光伏发电***工作的情况下进行故障检测与分类。
2、本发明基于小波多分辨分析的方法提取了4维的特征向量，用SVM训练故障分类模型，实现了光伏发电阵列的故障诊断和分类。这种检测方案环境适用性强，有较强的创新性，训练的模型能实现较高精度的故障诊断和分类。本发明利用小波多分辨率分析技术可以提取信号的细节系数，突出故障发生的特征。用SVM基于提取的特征向量训练的故障诊断模型，环境适用性强，实现了光伏发电阵列的准确地故障检测和分类。本发明的分类准确率达95％以上。
附图说明
图1为本发明实施例的方法流程示意图。
图2为本发明实施例的光伏发电***拓扑图。
图3为本发明实施例的光伏发电***实验平台图。
图4为本发明实施例的故障故障分类结果。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
应该指出，以下详细说明都是示例性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
如图1所示，本实施例提供了一种基于小波多分辨分析和SVM的光伏发电阵列故障诊断和分类的方法，图2为本实施例的光伏发电***拓扑图，***由S×P个太阳能组件组成，通过逆变器与电网进行连接实现并网发电，通过模拟光伏发电阵列出现的故障状况，包括开路1，短路1和短路2故障三种故障。在不同温照下，针对每种故障情况进行实时故障诊断，具体包括以下步骤：
步骤S1：采集光伏阵列不同工作状态下不同温照度的电流和电压样本信号；
步骤S2：对采集的电流和电压信号进行归一化处理，并得到电流变化率、电压变化率、功率变化率和电导变化率四个样本信号；
步骤S3：设置窗口大小为L的滑动窗对步骤S2得到的四个样本信号进行滑动取窗口信号，对取得的四种窗口信号进行小波变换的多分辨率分析；
步骤S4：计算四个窗口信号第N层分解的高频信号的2范数；
步骤S5：取四个样本信号在故障发生时的窗口信号的第N层分解的高频信号的2范数作为特征值，构成维度为4的特征向量；
步骤S6：由不同工作状态的多组样本信号得到多个特征值向量，其中将m个特征向量作为训练数据，将n个特征向量作为测试数据；
步骤S7：设置支持向量机SVM的参数，从训练数据中训练故障诊断模型，用测试数据检测该模型的分类的准确性。
较佳的，本实施例的***参数如下表所示。
表1***详细参数
本发明运用小波多分辨分析提取不同工作状态的光伏阵列电流变化率，电压变化率，功率变化率和电导变化率四个信号的细节系数作为特征，设置SVM参数从构造的4维特征向量中训练故障分类器，实现光伏阵列的故障诊断和分类。
在本实施例中，步骤S1中，光伏阵列的不同工作状态包括：正常、单组串1个组件短路(以下称短路1)、单组串一个组件开路(以下称开路1)以及单组串2个组件短路(以下称短路2)；其中，单组串1个组件短路、单组串一个组件开路以及单组串2个组件短路为故障状态，故障状态的样本信号包含了光伏阵列由正常到故障再到故障稳定的过程，捕获了故障发生时阵列电流和电压的变化特性。其中，电流和电压样本信号的采集时间窗口t＝10s，采样频率为200Hz。
本实施例的方法能够对不同光照度下的短路1，开路1和短路2故障进行检测。不同环境下光伏阵列的电信号在同种故障时有相似的变化特征，所提出的方法在组串式光伏发电***较强的环境适用性。特别地，本实施例模拟光伏发电系正常、短路1，开路1和短路2四种工作状态进行数据采集。在2018年8月份分多个时间段，在不同的温照度下进行样本信号随机采集，每种工作状态采集150个样本信号，运用小波多分辨率分析得到150个特征向量，随机选取100个特征向量作为训练数据，剩余的50个特征向量作为测试数据。图3为本实施例中光伏发电***实验平台图。采集的电流样本信号具体信息如表2所示。
表2样本信号采集信息
在本实施例中，步骤S2包括以下步骤：
步骤S21：对采集的电流和电压信号进行归一化处理：将阵列电流除以短路电流，将阵列电压除以开路电压；经过归一化处理能够克服不同光照度和不同温度的影响，归一化的公式如下：
iPV(t)＝IPV(t)/ISC(t)；
vPV(t)＝VPV(t)/VOC(t)
式中，IPV(t)表示采集的阵列电流样本信号，ISC(t)表示阵列短路电流信号，iPV(t)表示归一化后的阵列电流样本信号,VPV(t)表示采集的阵列电压样本信号，VOC(t)表示阵列开路信号,vPV(t)表示归一化后的阵列电压样本信号；经过归一化处理之后的电流和电压信号反映了不同工作状态下阵列信号的变化趋势，突出了变化特性；
步骤S22：采用下式得到电流变化率、电压变化率、功率变化率和电导变化率四个样本信号：
di(t)＝[iPV(n+1)-iPV(n)]/t；
dv(t)＝[vPV(n+1)-vPV(n)]/t；
dp(t)＝[iPV(n+1)v(n+1)-iPV(n)v(n)]/t；
dg(t)＝[iPV(n+1)/vPV(n+1)-iPV(n)/iPV(n)]/t；
式中，di(t)、dv(t)、dp(t)、dg(t)分别描述了光伏阵列不同工作状态的电流变化率，电压变化率，功率变化率和电导变化率的信号。
在本实施例中，步骤S3中，所述窗口大小L为20，采用DB小波对取得的四种窗口信号进行6层多分辨率分解的小波变换；所述多分辨率分析为用高通和低通滤波器将信号分解成多个级别的高频细节和近似值。
较佳的，在本实施例中，小波变换是一种应用广泛的数字信号处理方法，适用于分析非周期信号，可以同时对信号进行时间和频率分析。信号x(t)离散化为x(n),该信号的DWT可以由下式求解：
式中，ψ(t)是母小波，n是信号x(n)的点数，a＝a0 j是尺度因子，b＝kb0a0 j平移因子。a0＞0且为常数，b0≠0，j,k∈Z。
多分辨率分析是用高通和低通滤波器将信号分解成多个级别的高频细节和近似值。时域信号f(t)可以由尺度函数Φ(t)和小波函数Ψ(t)表示,如下式所示：
式中，dj(k)是各级分解的细节系数，表达高频信息，而aN(k)最后一级N级的近似系数，表达低频信息。
在本实施例中，步骤S4中，计算四个窗口信号第N层分解的高频信号的2范数采用下式：
式中，djn表示，N表示，dj表示第j层的高频系数,fi表示第i个信号提取的特征。本发明提取4个样本信号的第6层分解的高频系数作为特征值，故i＝1,2,3,4，j＝6。
在本实施例中，步骤S5中，所述维度为4的特征向量为f＝[f1,f2,f3,f4]。
在本实施例中，步骤S7中，所述设置支持向量机SVM的参数具体为：将SVM的惩罚因子C设置为1000，将两个不同类别的支持向量到超平面的距离之和γ设置为5。具体为：
支持向量机的基本思想是通过线性可分的训练样本集找到一个最优分类超平面以实现不同类别的样本数据的划分。对于非线性的分类问题，可以将样本从原始空间映射到一个线性可分的高维特征空间，从而找到一个合适的最优分类超平面。
给定训练样本数据集，D＝{xi,yi}，i＝1,2,3,...,m,yi∈{-1,1}，其中xi为样本数据，m为训练样本总数，d为样本空间的维度，yi为样本相应的标签。它们能够被一个最优超平面所分离，该超平面可以表示为wTx+b＝0，其中w∈Rd，为法向量，决定了超平面的方向；b∈R，为位移项，决定了超平面与原点间的距离，分类的阈值。
假设超平面能将训练样本正确分类，则对于{xi,yi}∈D，如果yi＝+1，则有wTx+b＞0；如果yi＝-1，则有wTx+b＜0。令：
则距离超平面最近的若干训练样本点使得上式的等号成立，它们被称之为支持向量(Support vectors，SVs)，则两个不同类别的支持向量到超平面的距离之和为该距离被称为分类间隔。如果要使γ最大，则要求||w||2最小，同时要求分类面对所有的样本正确分类，需满足
yi(wTxi+b)≥1,i＝1,2,3,...,l
因此，求最优分类超平面问题，可以转化为一个二次规划问题，优化目标可以写为
s.t.yi(wTxi+b)≥1,i＝1,2,3,...,l
对于在样本空间线性可分的训练样本，其可以被最优分类超平面进行划分。然后，在现实任务中，常常存在线性不可分的情况，此时在训练样本中存在部分样本数据不满足公式(1)，存在一定分类的误差。因此，通过引入松弛变量ξi(ξi≥0)来解决此问题。因此，公式(1)可写为
yi(wTxi+b)≥1-ξi,i＝1,2,3,...,l
在最大化间隔的同时，要求不满足约束的样本要尽可能少。因此，优化目标函数可以改写为
其中，称为惩罚项，C为惩罚因子。因此，可以得到线性不可分时的最优分类面，称为广义分类超平面，可以表示为下面的优化问题
s.t.yi(wTxi+b)≥1-ξi
ξi≥0,i＝1,2,3,...,l
求解该优化问题可以借助拉格朗日乘子法，得到它的对偶问题如下式：
0≤αi≤C,i,j＝1,2,3,...,l
根据Karush–Kuhn–Tucker(KKT)条件，可以得到αi(yi(wTxi+b)-1)＝0。如果αi＞0，则所对应的样本点位于最大间隔边界上，就是支持向量。则可以通过求解出w，并根据yi(wTxi+b)-1＝0求解出b，其中xi为支持向量，n为支持向量的数量。在确定w和b之后，得到分类决策函数如下所示
其次，针对非线性分类问题，SVM采用核函数将样本从原始空间映射到更高维的特征空间，使得样本在这个特征空间内线性可分。令φ(x)表示将x映射后的特征向量，则在高维空间中超平面所对应的优化目标可表示为
s.t.yi(wTφ(xi)+b)≥1-ξi
ξi≥0,i＝1,2,3,...,l
对应的对偶问题为
0≤αi≤C,i,j＝1,2,3,...,l
为避开计算样本xi和xj在高维空间中的内积运算，构造核函数K(·,·)，xi和xj在特征空间的内积转换为在原始样本空间中通过该函数计算的结果。K(·,·)表示如下
K(xi,xj)＝φ(xi)Tφ(xj)
则将前一个公式变为
0≤αi≤C,i,j＝1,2,3,...,l
此时分类决策函数变为
常用的核函数主要有多项式核函数、径向基(Radial based kernel，RBF)核函数、双曲正切(Sigmoid)核函数等。本发明中采用RBF核函数对SVM进行应用。
RBF核函数表示为
K(xi,x)＝exp(-γ||xi-x||2)
支持向量机的分类效果将很大程度取决于C和γ的选择，本发明中将C设置为1000，γ设置为10。运用设置好参数SVM，从构造的训练数据中训练出故障分类器，实现光伏发电阵列的故障诊断和分类。用50个测试数据检测该分类的精度，图4为所提方案的故障分类结果。
相对应，短路1的标签为1，开路1的标签为2，短路2的标签为3，短路2的标签为4。在分类结果图中，如果预测标签和实际标签重合，则该数据的预测结果是准确的。如图4所示，50个短路1测试数据中有1个数据的预测标签和实际标签不一致，诊断精度为98％，开路1测试数据中有4个数据的预测出错，预测精度为92％，开路2测试数据中有2个数据的预测出错，预测精度为96％，而正常测试数据的预测准确率为100％。该方案实现了整体96.5％精度的故障诊断和分类。综上，本实施例中的故障诊断与分类效果如表3所示。
表3故障检测与分类结果
以上所述仅为本发明的较佳实施例，凡依本发明申请专利范围所做的均等变化与修饰，皆应属本发明的涵盖范围。
Claims (7)
1.一种基于小波多分辨分析和SVM的光伏发电阵列故障诊断和分类的方法，其特征在于：包括以下步骤：
步骤S1：采集光伏阵列不同工作状态下不同温照度的电流和电压样本信号；
步骤S2：对采集的电流和电压信号进行归一化处理，并得到电流变化率、电压变化率、功率变化率和电导变化率四个样本信号；
步骤S3：设置窗口大小为L的滑动窗对步骤S2得到的四个样本信号进行滑动取窗口信号，对取得的四种窗口信号进行小波变换的多分辨率分析；
步骤S4：计算四个窗口信号第N层分解的高频信号的2范数；
步骤S5：取四个样本信号在故障发生时的窗口信号的第N层分解的高频信号的2范数作为特征值，构成维度为4的特征向量；
步骤S6：由不同工作状态的多组样本信号得到多个特征值向量，其中将m个特征向量作为训练数据，将n个特征向量作为测试数据；
步骤S7：设置支持向量机SVM的参数，从训练数据中训练故障诊断模型，用测试数据检测该模型的分类的准确性。
2.根据权利要求1所述的一种基于小波多分辨分析和SVM的光伏发电阵列故障诊断和分类的方法，其特征在于：步骤S1中，光伏阵列的不同工作状态包括：正常、单组串1个组件短路、单组串一个组件开路以及单组串2个组件短路；其中，单组串1个组件短路、单组串一个组件开路以及单组串2个组件短路为故障状态，故障状态的样本信号包含了光伏阵列由正常到故障再到故障稳定的过程，捕获了故障发生时阵列电流和电压的变化特性。
3.根据权利要求1所述的一种基于小波多分辨分析和SVM的光伏发电阵列故障诊断和分类的方法，其特征在于：步骤S2包括以下步骤：
步骤S21：对采集的电流和电压信号进行归一化处理：将阵列电流除以短路电流，将阵列电压除以开路电压；经过归一化处理能够克服不同光照度和不同温度的影响，归一化的公式如下：
iPV(t)＝IPV(t)/ISC(t)；
vPV(t)＝VPV(t)/VOC(t)
式中，IPV(t)表示采集的阵列电流样本信号，ISC(t)表示阵列短路电流信号，iPV(t)表示归一化后的阵列电流样本信号,VPV(t)表示采集的阵列电压样本信号，VOC(t)表示阵列开路信号,vPV(t)表示归一化后的阵列电压样本信号；
步骤S22：采用下式得到电流变化率、电压变化率、功率变化率和电导变化率四个样本信号：
di(t)＝[iPV(n+1)-iPV(n)]/t；
dv(t)＝[vPV(n+1)-vPV(n)]/t；
dp(t)＝[iPV(n+1)v(n+1)-iPV(n)v(n)]/t；
dg(t)＝[iPV(n+1)/vPV(n+1)-iPV(n)/iPV(n)]/t；
式中，di(t)、dv(t)、dp(t)、dg(t)分别描述了光伏阵列不同工作状态的电流变化率，电压变化率，功率变化率和电导变化率的信号。
4.根据权利要求1所述的一种基于小波多分辨分析和SVM的光伏发电阵列故障诊断和分类的方法，其特征在于：步骤S3中，所述窗口大小L为20，采用DB小波对取得的四种窗口信号进行6层多分辨率分解的小波变换；所述多分辨率分析为用高通和低通滤波器将信号分解成多个级别的高频细节和近似值。
5.根据权利要求1所述的一种基于小波多分辨分析和SVM的光伏发电阵列故障诊断和分类的方法，其特征在于：步骤S4中，计算四个窗口信号第N层分解的高频信号的2范数采用下式：
式中，djn表示，N表示，dj表示第j层的高频系数,fi表示第i个信号提取的特征。
6.根据权利要求5所述的一种基于小波多分辨分析和SVM的光伏发电阵列故障诊断和分类的方法，其特征在于：步骤S5中，所述维度为4的特征向量为f＝[f1,f2,f3,f4]。
7.根据权利要求1所述的一种基于小波多分辨分析和SVM的光伏发电阵列故障诊断和分类的方法，其特征在于：步骤S7中，所述设置支持向量机SVM的参数具体为：将SVM的惩罚因子C设置为1000，将两个不同类别的支持向量到超平面的距离之和γ设置为5。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811561729.0A CN109617526A (zh) | 2018-12-20 | 2018-12-20 | 一种基于小波多分辨分析和svm的光伏发电阵列故障诊断和分类的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811561729.0A CN109617526A (zh) | 2018-12-20 | 2018-12-20 | 一种基于小波多分辨分析和svm的光伏发电阵列故障诊断和分类的方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109617526A true CN109617526A (zh) | 2019-04-12 |
Family
ID=66009694
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811561729.0A Pending CN109617526A (zh) | 2018-12-20 | 2018-12-20 | 一种基于小波多分辨分析和svm的光伏发电阵列故障诊断和分类的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109617526A (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111046790A (zh) * | 2019-12-11 | 2020-04-21 | 江苏大学镇江流体工程装备技术研究院 | 一种泵轴承故障诊断方法 |
| CN113746132A (zh) * | 2021-08-11 | 2021-12-03 | 国网江苏省电力有限公司 | 一种基于云边协同的光伏电站及其控制方法 |
| CN113916536A (zh) * | 2021-09-03 | 2022-01-11 | 山东科技大学 | 一种正弦滚轴筛联轴器故障诊断方法及*** |
| CN114993438A (zh) * | 2022-05-30 | 2022-09-02 | 济南建和机电设备有限公司 | 一种利用分布式压力传感器监测粮仓库存的方法 |
| CN115758284A (zh) * | 2022-11-14 | 2023-03-07 | 南京工业大学 | 基于融合核熵与信息增益的光伏热斑故障检测方法及*** |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104638668A (zh) * | 2014-06-12 | 2015-05-20 | 华中科技大学 | 一种光伏发电并网控制方法及*** |
| CN104753461A (zh) * | 2015-04-10 | 2015-07-01 | 福州大学 | 基于粒子群优化支持向量机的光伏发电阵列故障诊断与分类方法 |
| CN105375461A (zh) * | 2015-07-10 | 2016-03-02 | 华北电力大学（保定） | 基于预测技术的主动配电网供电能力实时评估方法 |
| CN108875796A (zh) * | 2018-05-28 | 2018-11-23 | 福州大学 | 基于线性判别分析和支持向量机的光伏阵列故障诊断方法 |
-
2018
- 2018-12-20 CN CN201811561729.0A patent/CN109617526A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104638668A (zh) * | 2014-06-12 | 2015-05-20 | 华中科技大学 | 一种光伏发电并网控制方法及*** |
| CN104753461A (zh) * | 2015-04-10 | 2015-07-01 | 福州大学 | 基于粒子群优化支持向量机的光伏发电阵列故障诊断与分类方法 |
| CN105375461A (zh) * | 2015-07-10 | 2016-03-02 | 华北电力大学（保定） | 基于预测技术的主动配电网供电能力实时评估方法 |
| CN108875796A (zh) * | 2018-05-28 | 2018-11-23 | 福州大学 | 基于线性判别分析和支持向量机的光伏阵列故障诊断方法 |
Non-Patent Citations (2)
| Title |
|---|
| KUN XIA等: ""Wavelet Packet and Support Vector Machine Analysis of Series DC Arc Fault Detection in Photovoltaic System"", 《IEEJ TRANSACTIONS ON ELECTRICAL AND ELECTRONIC ENGINEERING》 * |
| ZHEHAN YI等: ""Line-to-Line Fault Detection for Photovoltaic Arrays Based on Multiresolution Signal Decomposition and Two-Stage Support Vector Machine"", 《IEEE TRANSACTIONS ON INDUSTRIAL ELECTRONICS》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111046790A (zh) * | 2019-12-11 | 2020-04-21 | 江苏大学镇江流体工程装备技术研究院 | 一种泵轴承故障诊断方法 |
| CN113746132A (zh) * | 2021-08-11 | 2021-12-03 | 国网江苏省电力有限公司 | 一种基于云边协同的光伏电站及其控制方法 |
| CN113746132B (zh) * | 2021-08-11 | 2024-01-30 | 国网江苏省电力有限公司 | 一种基于云边协同的光伏电站及其控制方法 |
| CN113916536A (zh) * | 2021-09-03 | 2022-01-11 | 山东科技大学 | 一种正弦滚轴筛联轴器故障诊断方法及*** |
| CN114993438A (zh) * | 2022-05-30 | 2022-09-02 | 济南建和机电设备有限公司 | 一种利用分布式压力传感器监测粮仓库存的方法 |
| CN115758284A (zh) * | 2022-11-14 | 2023-03-07 | 南京工业大学 | 基于融合核熵与信息增益的光伏热斑故障检测方法及*** |
| CN115758284B (zh) * | 2022-11-14 | 2023-05-16 | 南京工业大学 | 基于融合核熵与信息增益的光伏热斑故障检测方法及*** |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108647716B (zh) | 一种基于复合信息的光伏阵列故障诊断方法 | |
| CN106777984B (zh) | 一种基于密度聚类算法实现光伏阵列工作状态分析与故障诊断的方法 | |
| CN109617526A (zh) | 一种基于小波多分辨分析和svm的光伏发电阵列故障诊断和分类的方法 | |
| CN104753461B (zh) | 基于粒子群优化支持向量机的光伏发电阵列故障诊断与分类方法 | |
| CN109738776A (zh) | 基于lstm的风机变流器开路故障识别方法 | |
| CN109460618A (zh) | 一种滚动轴承剩余寿命在线预测方法及*** | |
| CN109660206B (zh) | 一种基于Wasserstein GAN的光伏阵列故障诊断方法 | |
| CN109672406B (zh) | 一种基于稀疏表示和svm的光伏发电阵列故障诊断与分类的方法 | |
| CN110336534A (zh) | 一种基于光伏阵列电气参数时间序列特征提取的故障诊断方法 | |
| Zaki et al. | Deep‐learning–based method for faults classification of PV system | |
| CN109871976A (zh) | 一种基于聚类及神经网络的含分布式电源配电网电能质量预测方法 | |
| CN108983749A (zh) | 基于k-svd训练稀疏字典的光伏阵列故障诊断方法 | |
| CN109766952A (zh) | 基于偏最小二乘法和极限学习机的光伏阵列故障检测方法 | |
| CN107807860B (zh) | 一种基于矩阵分解的电力故障分析方法及*** | |
| Yurtseven et al. | Sensorless fault detection method for photovoltaic systems through mapping the inherent characteristics of PV plant site: Simple and practical | |
| Ying et al. | Automatic detection of photovoltaic module cells using multi-channel convolutional neural network | |
| Iturrino-García et al. | An innovative single shot power quality disturbance detector algorithm | |
| Zu et al. | A simple gated recurrent network for detection of power quality disturbances | |
| He et al. | Compound fault diagnosis for photovoltaic arrays based on multi-label learning considering multiple faults coupling | |
| Serrano-Luján et al. | Case of study: Photovoltaic faults recognition method based on data mining techniques | |
| Hussain et al. | Deployment of AI-based RBF network for photovoltaics fault detection procedure | |
| CN113222067B (zh) | 基于SVM-Adaboost算法的孤岛智能检测方法 | |
| Zgraggen et al. | Physics informed deep learning for tracker fault detection in photovoltaic power plants | |
| Lu et al. | Efficient fault diagnosis approach for solar photovoltaic array using a convolutional neural network in combination of generative adversarial network under small dataset | |
| Gnetchejo et al. | Faults diagnosis in a photovoltaic system based on multivariate statistical analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190412 |
|
| RJ01 | Rejection of invention patent application after publication |