CN108509736B - 一种用于故障诊断的光伏组件参数提取方法 - Google Patents

Abstract

本发明提出一种用于故障诊断的光伏组件参数提取方法。由于组件的旁路二极管的存在，在组件中部分电池片受到遮挡时，其I‑V输出特性会出现多个膝点，因此其I‑V曲线又可分为若干个子I‑V曲线。采用传统解析法只能针对完整的I‑V曲线提取光伏组件五参数，对于由于多膝点而产生的子I‑V曲线，表征该子I‑V特性下的光伏组件五参数，传统利用解析法求解光伏组件无参数的方法均不再适应。本发明所提方法所需计算点不依赖完整的I‑V特性曲线，可用于解决求解多膝点下子I‑V曲线对应组件五参数的问题，为故障诊断提供依据。

Description

一种用于故障诊断的光伏组件参数提取方法

技术领域

本发明属于光伏发电技术领域，具体讲涉及一种用于故障诊断的光伏组件参数提取方法。

背景技术

近年来，光伏发电作为国内热门的研发投资领域，引发了一股研究太阳光伏发电的热潮。对于太阳能光伏发电的研究从光生伏特效应发现以来，一直不断进步与发展，国内对光伏领域的研究也日趋深入。

一般来说，一个光伏阵列由光伏组串并联而成，光伏组串又由光伏组件串联而成，一个光伏组件又由固定个数的太阳能电池片串联而成；这类光伏组件之间的串并联连接结合后续电能转换器件，构成了一个初步的光伏发电系统。

然而，光伏组件的发电会受到外界环境的干扰，具体而言主要为光强与环境温度的影响，光照强度降低、环境温度上升，都会导致组件输出电能降低。此外，组件内部由于生产工艺水平的不同而导致的串联失配问题，亦会导致整体输出电能降低。由于光伏组串由若干个光伏组件串联而成，因此组件与组件之间的失配会带来组串整体输出功率的损失，甚至导致部分组件发生热斑等故障；为减少这类故障的发生，一般在光伏组件封装时即会在若干个串联的太阳能电池片上反向并联一个旁路二极管，其目的在于当组串的整体输出电流大于旁路二极管所并联的太阳能电池片输出的电流时，为电池片提供短路通路，以防电流反向灌入太阳能电池片引发发热甚至烧毁等问题。

旁路二极管一定程度上解决了组件之间因失配而产生的热斑问题，但造成了组串整体输出呈现了非连续变化的性质，表现在旁路二极管导通与关断状态前后，光伏组串的I-V输出特性曲线在该状态点发生突变，光伏组串的P-V输出特性在该状态点出现波谷；综上，一旦光伏组串中有旁路二极管导通，其P-V输出特性即出现多个功率峰值，且单个峰值区间上P-V曲线不再为该工作状态下的组件从短路电流至开路电压的完整曲线。

目前，一般通过解析法求解光伏组件五参数(光生电流I_ph、二极管理想因子n、串联电阻R_s、并联电阻R_sh和二极管反向饱和电流I_o)均需要实测以下数据：短路电流I_sc、开路电压U_oc、最大功率点(U_mpp,I_mpp)以及开路电压处与短路电流处的dU/dI；获得这些数据要求光伏组件提供出完整的从短路电路至开路电压的I-V曲线，然而如前文所述，当组串中因部分组件失配而导致旁路二极管工作时，组件无法提供完整的I-V曲线，因此以目前解析法的手段无法求出光伏组串中各组件的五参数。

光伏组件的五参数较为全面的描述了光伏组件的工作状态与使用性能，因此对于光伏组件的故障诊断来说，提取其五参数是至关重要的，本发明提供了一种用于故障诊断的光伏组件参数提取方法，该方法在组件无法提供完整的I-V曲线时，仍可以较好地提取出光伏组件在该工作状态下的五参数，为光伏组件的故障诊断提供了依据。

发明内容

光伏组串的五参数是表征光伏组串工作性能的重要指标，本发明提出一种用于故障诊断的光伏组件参数提取方法，其特征在于该光伏组件I-V特性曲线不完整时，通过采样该光伏组件的部分电压电流数据点求解出表征光伏组件物理性质的五参数(光生电流I_ph、二极管反向饱和电流I_o、二极管理想因子n、串联电阻R_s、并联电阻R_sh)，具体步骤如下：

(1)根据光伏组件单二极管物理模型，推导出新的光伏组件五参数解析式；由于在该解析式中不再需要最大功率点、开路电压值以及开路电压处的dU/dI，因此减少了对光伏组件完整I-V数据的依赖性；

(2)将单块光伏组件连接至电压调制模块，采样计算五参数所需的7个数据点；

(3)根据光伏组件五参数解析式采样相应电压电流数据点，计算出该光伏组件五参数。

步骤(1)中的光伏组件五参数解析式推导如下：

(一)根据光伏组件单二极管等效电路模型，并基于该模型利用基尔霍夫电流定律(KCL)得到光伏组件电流特性方程；

其中I_ph为光生电流，I_D为通过二极管的电流，I_R为通过并联电阻的电流，I_o为二极管反向饱和电流，R_s为串联电阻，R_sh为并联电阻，n为二极管理想因子，I为光伏组件输出电流，U为光伏组件输出电压，q为电子电荷常数，T为环境温度，k为玻尔兹曼常数，后文中，令

(二)对(1)式求导，可得：

为替换(2)式中的幂指数式，将(1)式带入(2)式，且I_ph≈I_sc得：

(三)为获得更多组件参数信息，分别将短路电流点(I＝I_sc，U＝0)与开路电压点(I＝0，U＝U_oc)带入(1)和(2)式，其中I_sc为光伏组件短路电流，U_oc为光伏组件开路电压，可得：

(四)由(4)式与(6)式可以消去I_ph，可得式(7)：

(五)将(5)式变形可得：

(六)对(7)式进行化简：

化简原则：R_sh≈R_sho＞＞R_s (9)

得到：

(七)由(3)式可得R_s与任一点R_so ^*的关系，结合(9)式化简如下：

其中R_so*为I-V曲线上任意一点的斜率，I*与U*为光伏组件中I-V曲线上的任意一点处的电压和电流值，由于(12)式涉及化简，U^*的取值接近U_oc时计算值较为准确，其中U_oc为组件开路电压值；

(八)将(9)式带入(12)式化简得：

(九)将光伏组件电压点(U_a，I_a)带入(1)式，可得：

将(14)式带入(4)式，消去I_ph，可得：

并将(8)式带入(15)式，消去

可得：

(十)对(16)式进行化简：

化简原则：

得到：

由于已经对(16)式做出了适当化简，因此实际采样点时，U_a的取值一般接近U_oc时计算结果较为准确；

(十一)将(11)式带入(17)式消去I_o，可得：

(十二)将(18)式两边取对数可得：

引入另一个电压调制点(U_b，I_b)，该点不能与(U_a，I_b)重合，且靠近U_oc为宜，最终可消去(19)式中等号右边的量，可得：

(十三)将(12)式带入(20)式消去R_s，

将n·V_th看做整体变量，得到一个关于n·V_th的一元二次方程：

形如：A(n·V_th)²+B(n·V_th)+C＝0 (21)

其中：

故：

(十四)求出二极管理想因子n后，可以利用式(13)求出串联电阻R_s：

(十五)为求得I_o，引入另一电压调制点(U_c，I_c)，该点不能与(U_a，I_a)、(U_b，I_b)重合，根据(1)式可推导出(27)式：

式(27)为二元一次方程，结合R_sho≈R_sh故可得到I_o与I_ph的表达式：

(十六)结合(5)式求出R_sh表达式：

结合上文公式(25)(26)(28)(29)(30)的光伏组件五参数解析式，需采样的电压电流数据点如下：

①采样求解

其中：

结合上式，求出nV_th所需的采样点为：

(1)短路电流点(0，I_sc)；

(2)短路电流附近点(δV，I_sc+δI)，其中，δV为电压调制模块调制电压值的最小步长，从而

(3)(U_a，I_a)点，该点电压U_a应靠近U_oc，实际U_a可取0.8U_oc；

(4)(U_b，I_b)点，该点电压U_b应靠近U_oc，且不能等于U_a；

(5)(U^*，I^*)点，该点电压U^*不能等于U_a和U_b，实际U^*可取0.7U_oc；

(6)(U^*，I^*)附近点(δV+U^*，I^*+δI)，其中，δV为电压调制模块调制电压值的最小步长，从而

②求解R_s，

由此可见，其中计算I_ph所需的数据点与计算nV_th所需的数据点相同；

③求解I_o，

计算I_o所需的采样点，除了计算nV_th所需的6个点之外，还需采样(U_c，I_c)点，该点的电压值U_c不能等于U_a；

④求解I_ph，

其中计算I_ph所需的数据点与计算nV_th所需的数据点相同；

⑤求解R_sh，

其中计算R_sh所需的数据点与计算nV_th所需的数据点相同；

综上所述，计算光伏组件五参数所需的7个采样点为：

a.短路电流点(0，I_sc)；

b.短路电流附近点(δV，I_sc+δI)，其中，δV为电压调制模块调制电压值的最小步长，从而

c.(U_a，I_a)点，该点电压U_a应靠近U_oc，实际U_a可取0.8U_oc；

d.(U_b，I_b)点，该点电压U_b应靠近U_oc，且不能等于U_a；

e.(U^*，I^*)点，该点电压U^*不能等于U_a和U_b，实际U^*可取0.7U_oc；

f.(U^*，I^*)附近点(δV+U^*，I^*+δI)，其中，δV为电压调制模块调制电压值的最小步长，从而

g.(U_c，I_c)点，该点的电压值U_c不能等于U_a；

光伏组串由若干个光伏组件串联，实际过程中，由于部分阴影遮挡或组件之间参数失配等问题，光伏组串的输出P-V特性曲线会出现多峰，多峰即意味着组件的I-V曲线不再是从短路电流到开路电压的完后只能曲线，因此，采用传统用于处理单峰P-V曲线的组件五参数求解方法不再适用；本发明针对组串输出多峰P-V曲线的情况，提出一种用于故障诊断的光伏组件参数提取方法，该方法可以提取出光伏组串中各种工作状态下组件的五参数，从而作为评价各个光伏组件工作状态与性能的依据。

附图说明

图1为光伏组件等效电路拓扑结构图；

图2为本发明所需采样点位置示意图；

图3为本发明针对提取的光伏组件五参数的仿真I-V曲线与实测I-V曲线对比图；

具体实施方式

本发明提供了一种用于故障诊断的光伏组件参数提取方法，结合光伏组件等效电路，推导出提取组串中组件五参数的解析公式。本发明所提方法适用于在光伏组件无法获知其完整的P-V曲线时提取其五参数，因此拓展了光伏组串P-V曲线输出多峰时的组件五参数提取方法。

太阳电池是利用P-N结的光伏效应将光能转化为电能的，在光照情况下，电流的电流-电压特性表现为二极管的特点，即电流随电压指数变化。如果将太阳能电池视为恒流源，将半导体材料体电阻、电极与半导体材料接触电阻和界面处载流子的复合电阻视为串联电阻，将电池边缘漏电流和电池金属桥漏电流等效为并联电阻，可获得太阳能电池等效电路参见图1。

若干个光伏组件串联称为光伏组串，光伏组串中，旁路二极管之所以导通，是因为该旁路二极管所反向并联的失配光伏组件输出电流小于其余正常组件的输出电流，当组串两端电压逐渐增加，正常组件的输出电流相应减小，当其电流值减小至失配光伏组件的短路电流时，旁路二极管由导通切换为截止状态，此时由正常光伏组件与失配光伏组件共同为负载输出电流；综上所述，光伏组件输出I-V曲线呈现多膝点，每个膝点之间即为不完整的I-V特性曲线，因此传统的解析法求解光伏组件五参数中需要采样的开路电压点、开路电压点处的dU/dI无法获知，从而采用传统的解析法无法求解此时组件的五参数。

为避免上述两个采样点的采集问题，本发明提出了一种用于故障诊断的光伏组件参数提取方法，该方法分为三个步骤：

步骤一：推导光伏组件五参数解析公式；

步骤二：根据五参数解析式采样相应的计算点；

步骤三：处理计算采样点数据，带入步骤一中公式求出五参数；

步骤一中解析推导如下：

(一)、根据光伏组件单二极管等效电路模型，并基于该模型利用基尔霍夫电流定律(KCL)得到光伏组件电流特性方程；

其中I_ph为光生电流，I_D为通过二极管的电流，I_R为通过并联电阻的电流，I_o为二极管反向饱和电流，R_s为串联电阻，R_sh为并联电阻，n为二极管理想因子，I为光伏组件输出电流，U为光伏组件输出电压，q为电子电荷常数，T为环境温度，k为玻尔兹曼常数，后文中，将

(二)、对(1)式求导，可得：

为解决(2)式中的幂指数式，由公式(1)可得：

又I_ph≈I_sc，将其带入式(2)(3)可得出式(4)：

(三)、为获得更多组件参数信息，将短路电流点(I＝I_sc，U＝0)和开路电压点(I＝0，U＝U_oc)带入(1)和(2)式，可得：

(四)、由(5)式与(7)式可以消去I_ph，可得下式：

(五)、将(6)式变形可得：

(六)、对(8)式进行化简：

化简原则：R_sh≈R_sho＞＞R_s

得到：

(七)、由(4)式可得R_s与任一点的

的关系为：

由于I_o数量级相对于I_sc等可以忽略，且

故式(12)可化简为式(13)

其中I^*与U^*为光伏组件的电压和电流值，

为该点处的-dU/dI，由于(13)式涉及化简，U^*的取值接近U_oc时计算值较为准确，其中U_oc为组件开路电压值；

(八)、继续化简式(13)

由于R_sh≈R_sho，故式(13)可化简为：

(九)、设光伏组件电压点(U_a，I_a)带入(1)式，可得：

将(15)式带入(4)式，消去I_ph，可得：

并将(8)式带入(16)式，消去

可得：

(十)、对(17)式进行化简：

化简原则：

得到：

由于已经对(17)式做出了适当化简，因此实际采样点时，U_a的取值一般接近U_oc时计算结果较为准确；

(十一)、将(15)带入(18)消去I_o，可得：

(十二)、将(19)式两边取对数可得：

引入另一个电压调制点(U_b，I_b)，该点不能与(U_a，I_a)重合，且靠近U_oc为宜，得到：

(十三)、联立式(20)与式(21)可得：

(十四)、将(13)式带入(22)式消去R_s，

将n·V_th看做整体变量，得到一个关于n·V_th的一元二次方程：

形如：A(n·V_th)²+B(n·V_th)+C＝0 (23)

其中：

故：

(十五)、求出二极管理想因子n后，可以利用式(13)求出串联电阻R_s：

(十六)、为求得I_o，引入另一电压调制点(U_c，I_c)，该点不能与(U_a，I_a)、(U_b，I_b)重合，根据(1)式可推导出(29)式：

式(29)为二元一次方程，结合R_sho≈R_sh故可得到I_o与I_ph的表达式：

(十七)、结合(5)式求出R_sh表达式：

步骤二中，根据推导公式，我们可以明确求出光伏组件五参数所需要的点如下：

①短路电流点(0，I_sc)；

②短路电流附近点(δV，I_sc+δI)，其中，δV为电压调制模块调制电压值的最小步长，从而

③(U_a，I_a)点，该点电压U_a应靠近U_oc，实际U_a可取0.8U_oc；

④(U_b，I_b)点，该点电压U_b应靠近U_oc，且不能等于U_a；

⑤(U^*，I^*)点，该点电压U^*不能等于U_a和U_b，实际U^*可取0.7U_oc；

⑥(U^*，I^*)附近点(δV+U^*，δI+I^*)，其中，δV为电压调制模块调制电压值的最小步长，从而

⑦(U_c，I_c)点，该点的电压值U_c不能等于U_a。

根据采样所需的七个点结合式(27)(28)(30)(31)(32)可以求出光伏组件的五参数。

综上所述，在光伏组件I-V曲线缺失开路电压附近I-V输出特性曲线时，依旧可以采用上述七个采样点依据五参数解析公式计算得出光伏组件在该状态下的五参数。七个采样点位置示意图如图2所示。

为验证该方法的准确性与可靠性，根据经典实验数据求出按该方法求出曲线数据与测试数据做对比，结果验证了该方法的可行性，实验数据与仿真数据对比如图3所示。由图3可知该方法的准确性较高，可以较为准确地求出光伏组件五参数，因此本发明所提出方法具有一定实用性。

Claims (1)

1.一种用于故障诊断的光伏组件参数提取方法，其特征在于：该光伏组件I-V特性曲线不完整时，通过采样该光伏组件的部分电压电流数据点求解出表征光伏组件物理性质的五参数，所述五参数包括光生电流I_ph、二极管反向饱和电流I_o、二极管理想因子n、串联电阻R_s、并联电阻R_sh，具体步骤如下：

步骤一：根据光伏组件单二极管物理模型，推导出新的光伏组件五参数解析式，具体推导过程如下：

①根据光伏组件单二极管等效电路模型，并基于该模型利用基尔霍夫电流定律(KCL)得到光伏组件电流特性方程：

其中I_ph为光生电流，I_D为通过二极管的电流，I_R为通过并联电阻的电流，I_o为二极管反向饱和电流，R_s为串联电阻，R_sh为并联电阻，n为二极管理想因子，I为光伏组件输出电流，U为光伏组件输出电压，q为电子电荷常数，T为环境温度，k为玻尔兹曼常数，后文中，令

②对(1)式求导，可得：

为替换(2)式中的幂指数式，将(1)式代入(2)式，且I_ph≈I_sc得：

③为获得更多组件参数信息，分别将短路电流点(I＝I_sc，U＝0)与开路电压点(I＝0，U＝U_oc)代入(1)和(2)式，其中I_sc为光伏组件短路电流，U_oc为光伏组件开路电压，可得：

式中R_sho是短路电流处的-dU/dI；

④由(4)式与(6)式可以消去I_ph，可得式(7)：

⑤将(5)式变形可得：

⑥对(7)式进行化简：

化简原则：R_sh≈R_sho＞＞R_s (9)

得到：

⑦由(3)式可得R_s与任一点R_so ^*的关系，结合(9)式化简如下：

其中R_so*为I-V曲线上任意一点的斜率，I^*与U^*为光伏组件中I-V曲线上的任意一点处的电压和电流值，由于(12)式涉及化简，U^*的取值接近U_oc时计算值更为准确，其中U_oc为组件开路电压值；

⑧将(9)式代入(12)式化简得：

⑨将光伏组件电压调制点(U_a，I_a)代入(1)式，可得：

将(14)式代入(4)式，消去I_ph，可得：

并将(8)式代入(15)式，消去

可得：

⑩对(16)式进行化简：

化简原则：

与(9)式；

得到：

将(11)式代入(17)式消去I_o，可得：

将(18)式两边取对数可得：

引入另一个电压调制点(U_b，I_b)，该点不能与(U_a，I_a)重合，最终可消去(19)式中等号右边的量，可得：

将(12)式代入(20)式消去R_s，

将n·V_th看做整体变量，得到一个关于n·V_th的一元二次方程：

形如：A(n·V_th)²+B(n·V_th)+C＝0 (20)

其中：

故：

求出二极管理想因子n后，可以利用式(13)求出串联电阻R_s：

为求得I_o，引入另一电压调制点(U_c，I_c)，该点不能与(U_a，I_a)、(U_b，I_b)重合，根据(1)式可推导出(26)式：

式(26)为二元一次方程，结合R_sho≈R_sh故可得到I_o与I_ph的表达式：

结合(5)式求出R_sh表达式：

其中上式中所引入的电压调制点(U_a，I_a)、(U_b，I_b)、(U_c，I_c)中的U_a、I_a、U_b、I_b、U_c、I_c均为光伏组件中I-V曲线上的任意一点处的电压值与电流值；

由上式(24)(25)(27)(28)(29)可知，所述新的光伏组件五参数解析式中不需要最大功率点、开路电压值以及开路电压处的dU/dI；

步骤二：将单块光伏组件连接至电压调制模块，采样计算五参数所需的7个数据点；所述7个数据点具体为：

①采样求解n·V_th，

其中：

结合上式，求出nV_th所需的数据点为：

(1)短路电流点(0，I_sc)；

(2)短路电流附近点(δV，I_sc+δV)，其中，δV为电压调制模块调制电压值的最小步长，δI为电压调制模块调制电流值的最小步长，从而：

(3)(U_a，I_a)点，该点电压U_a实际取0.8U_oc；

(4)(U_b，I_b)点，该点电压U_b应靠近U_oc，且不等于U_a；

(5)(U^*，I^*)点，该点电压U^*不等于U_a和U_b，实际U^*取0.7U_oc；

(6)(U^*，I^*)附近点(δV+U^*，I^*+δI)，其中，δV为电压调制模块调制电压值的最小步长，δI为电压调制模块调制电流值得的小步长，从而：

②求解R_s，

由此可见，求R_s所需的数据点与求nV_th相同；

③求解I_o，

计算I_o所需的数据点，除了计算nV_th所需的6个点之外，还需采样(U_c，I_c)点，该点的电压值U_c不能等于U_a；

④求解I_ph，

其中计算I_ph所需的数据点与计算nV_th所需的数据点相同；

⑤求解R_sh，

其中计算R_sh所需的数据点与计算nV_th所需的数据点相同；

综上所述，计算所述新的光伏组件五参数所需的7个数据点为：

a.短路电流点(0，I_sc)；

b.短路电流附近点(δV，I_sc+δI)，其中，δV为电压调制模块调制电压值的最小步长，从而：

c.(U_a，I_a)点，该点电压U_a实际取0.8U_oc；

d.(U_b，I_b)点，该点电压U_b应靠近U_oc，且不等于U_a；

e.(U^*，I^*)点，该点电压U^*不等于U_a和U_b，实际U^*取0.7U_oc；

f.(U^*，I^*)附近点(δV+U^*，I^*+δI)，其中，δV为电压调制模块调制电压值的最小步长，δI为电压调制模块调制电流值的最小步长，从而：

g.(U_c，I_c)点，该点的电压值U_c不等于U_a；

步骤三：将采样点数据代入步骤一中所述的新的光伏组件五参数解析式，计算出该光伏组件五参数。

Images