本发明属于光伏发电技术领域，具体讲涉及一种用于故障诊断的光伏组件参数提取方法。
背景技术：
近年来，光伏发电作为国内热门的研发投资领域，引发了一股研究太阳光伏发电的热潮。对于太阳能光伏发电的研究从光生伏特效应发现以来，一直不断进步与发展，国内对光伏领域的研究也日趋深入。
一般来说，一个光伏阵列由光伏组串并联而成，光伏组串又由光伏组件串联而成，一个光伏组件又由固定个数的太阳能电池片串联而成；这类光伏组件之间的串并联连接结合后续电能转换器件，构成了一个初步的光伏发电系统。
然而，光伏组件的发电会受到外界环境的干扰，具体而言主要为光强与环境温度的影响，光照强度降低、环境温度上升，都会导致组件输出电能降低。此外，组件内部由于生产工艺水平的不同而导致的串联失配问题，亦会导致整体输出电能降低。由于光伏组串由若干个光伏组件串联而成，因此组件与组件之间的失配会带来组串整体输出功率的损失，甚至导致部分组件发生热斑等故障；为减少这类故障的发生，一般在光伏组件封装时即会在若干个串联的太阳能电池片上反向并联一个旁路二极管，其目的在于当组串的整体输出电流大于旁路二极管所并联的太阳能电池片输出的电流时，为电池片提供短路通路，以防电流反向灌入太阳能电池片引发发热甚至烧毁等问题。
旁路二极管一定程度上解决了组件之间因失配而产生的热斑问题，但造成了组串整体输出呈现了非连续变化的性质，表现在旁路二极管导通与关断状态前后，光伏组串的i-v输出特性曲线在该状态点发生突变，光伏组串的p-v输出特性在该状态点出现波谷；综上，一旦光伏组串中有旁路二极管导通，其p-v输出特性即出现多个功率峰值，且单个峰值区间上p-v曲线不再为该工作状态下的组件从短路电流至开路电压的完整曲线。
目前，一般通过解析法求解光伏组件五参数(光生电流iph、二极管理想因子n、串联电阻rs、并联电阻rsh和二极管反向饱和电流io)均需要实测以下数据：短路电流isc、开路电压uoc、最大功率点(umpp,impp)以及开路电压处与短路电流处的du/di；获得这些数据要求光伏组件提供出完整的从短路电路至开路电压的i-v曲线，然而如前文所述，当组串中因部分组件失配而导致旁路二极管工作时，组件无法提供完整的i-v曲线，因此以目前解析法的手段无法求出光伏组串中各组件的五参数。
光伏组件的五参数较为全面的描述了光伏组件的工作状态与使用性能，因此对于光伏组件的故障诊断来说，提取其五参数是至关重要的，本发明提供了一种组件五参数提取方法，该方法在组件无法提供完整的i-v曲线时，仍可以较好地提取出光伏组件在该工作状态下的五参数，为光伏组件的故障诊断提供了依据。
技术实现要素：
光伏组串的五参数是表征光伏组串工作性能的重要指标，本发明提出一种利用解析法提取光伏组件五参数的方法，其特征在于该光伏组件i-v特性曲线不完整时，仍可以通过采样该光伏组件的部分电压电流数据点求解出表征光伏组件物理性质的五参数(光生电流iph、二极管反向饱和电流io、二极管理想因子n、串联电阻rs、并联电阻rsh)，具体步骤如下：
(1)根据光伏组件单二极管物理模型，推导出新的光伏组件五参数解析式；由于在该解析式中不再需要最大功率点、开路电压值以及开路电压处的du/di，因此减少了对光伏组件完整i-v数据的依赖性；
(2)将单块光伏组件连接至电压调制模块，采样计算五参数所需的7个数据点；
(3)根据光伏组件五参数解析式采样相应电压电流数据点，计算出该光伏组件五参数。
步骤(1)中的光伏组件五参数解析式推导如下：
(一)根据光伏组件单二极管等效电路模型，并基于该模型利用基尔霍夫电流定律(kcl)得到光伏组件电流特性方程；
其中iph为光生电流，id为通过二极管的电流，ir为通过并联电阻的电流，io为二极管反向饱和电流，rs为串联电阻，rsh为并联电阻，n为二极管理想因子，i为光伏组件输出电流，u为光伏组件输出电压，q为电子电荷常数，t为环境温度，k为玻尔兹曼常数，后文中，令
(二)对(1)式求导，可得：
为替换(2)式中的幂指数式，将(1)式带入(2)式，且iph≈isc得：
(三)为获得更多组件参数信息，分别将短路电流点(i＝isc，u＝0)与开路电压点(i＝0，u＝uoc)带入(1)和(2)式，其中isc为光伏组件短路电流，uoc为光伏组件开路电压，可得：
(四)由(4)式与(6)式可以消去iph，可得式(7)：
(五)将(5)式变形可得：
(六)对(7)式进行化简：
化简原则：rsh≈rsho＞＞rs(9)
得到：
(七)由(3)式可得rs与任一点rso*的关系，结合(9)式化简如下：
其中i*与u*为光伏组件的电压电流点，由于(12)式涉及化简，u*的取值接近uoc时计算值较为准确，其中uoc为组件开路电压值；
(八)将(9)式带入(12)式化简得：
(九)设光伏组件电压点(ua，ia)带入(1)式，可得：
将(14)式带入(4)式，消去iph，可得：
并将(8)式带入(15)式，消去
(十)对(16)式进行化简：
化简原则：
得到：
由于已经对(16)式做出了适当化简，因此实际采样点时，ua的取值一般接近uoc时计算结果较为准确；
(十一)将(11)带入(17)消去io，可得：
(十二)将(18)式两边取对数可得：
引入另一个电压调制点(ub，ib)，该点不能与(ua，ib)重合，且靠近uoc为宜，最终可消去(19)式中等号右边的量，可得：
(十三)将(12)式带入(20)式消去rs，
将n·vth看做整体变量，得到一个关于n·vth的一元二次方程：
形如：a(n·vth)2+b(n·vth)+c＝0(21)
其中：
故：
(十四)求出二极管理想因子n后，可以利用式(13)求出串联电阻rs：
(十五)为求得io，引入另一电压调制点(uc，ic)，该点不能与(ua，ia)、(ub，ib)重合，根据(1)式可推导出(27)式：
式(27)为二元一次方程，结合rsho≈rsh故可得到io与iph的表达式：
(十六)结合(5)式求出rsh表达式：
结合上文公式(25)(26)(28)(29)(30)的光伏组件五参数解析式，需采样的电压电流数据点如下：
①采样求解n·vth，
结合上式，求出nvth所需的采样点为：
(1)短路电流点(0，isc)；
(2)短路电流附近点(δv，isc+δi)，其中，δv为电压调制模块调制电压值得最小步长，从而
(3)(ua，ia)点，该点电压ua应靠近uoc，实际ua可取0.8uoc；
(4)(ub，ib)点，该点电压ub应靠近uoc，且不能等于ua；
(5)(u*，i*)点，该点电压u*不能等于ua和ub，实际u*可取0.7uoc；
(6)(u*，i*)附近点(δv+u*，i*+δi)，其中，δv为电压调制模块调制电压值得最小步长，从而
③求解rs，
③求解io，
④求解iph，
⑤求解rsh，
a.短路电流点(0，isc)；
b.短路电流附近点(δv，isc+δi)，其中，δv为电压调制模块调制电压值得最小步长，从而
c.(ua，ia)点，该点电压ua应靠近uoc，实际ua可取0.8uoc；
d.(ub，ib)点，该点电压ub应靠近uoc，且不能等于ua；
e.(u*，i*)点，该点电压u*不能等于ua和ub，实际u*可取0.7uoc；
f.(u*，i*)附近点(δv+u*，i*+δi)，其中，δv为电压调制模块调制电压值得最小步长，从而
g.(uc，ic)点，该点的电压值uc不能等于ua。
光伏组串由若干个光伏组件串联，实际过程中，由于部分阴影遮挡或组件之间参数失配等问题，光伏组串的输出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)得到光伏组件电流特性方程；
其中iph为光生电流，id为通过二极管的电流，ir为通过并联电阻的电流，io为二极管反向饱和电流，rs为串联电阻，rsh为并联电阻，n为二极管理想因子，i为光伏组件输出电流，u为光伏组件输出电压，q为电子电荷常数，t为环境温度，k为玻尔兹曼常数，后文中，将
(二)、对(1)式求导，可得：
为解决(2)式中的幂指数式，由公式(1)可得：
又iph≈isc，将其带入式(2)(3)可得出式(4)：
(三)、为获得更多组件参数信息，将短路电流点(i＝isc，u＝0)和开路电压点(i＝0，u＝uoc)带入(1)和(2)式，可得：
(四)、由(5)式与(7)式可以消去iph，可得下式：
(五)、将(6)式变形可得：
(六)、对(8)式进行化简：
化简原则：rsh≈rsho＞＞rs
得到：
(七)、由(4)式可得rs与任一点的
由于io数量级相对于isc等可以忽略，且
其中i*与u*为光伏组件的电压电流点，
(八)、继续化简式(13)
由于rsh≈rsho，故式(13)可化简为：
(九)、设光伏组件电压点(ua，ia)带入(1)式，可得：
将(15)式带入(4)式，消去iph，可得：
并将(8)式带入(16)式，消去
(十)、对(17)式进行化简：
化简原则：
得到：
由于已经对(17)式做出了适当化简，因此实际采样点时，ua的取值一般接近uoc时计算结果较为准确；
(十一)、将(15)带入(18)消去io，可得：
(十二)、将(19)式两边取对数可得：
引入另一个电压调制点(ub，ib)，该点不能与(ua，ia)重合，且靠近uoc为宜，得到：
(十三)、联立式(20)与式(21)可得：
(十四)、将(13)式带入(22)式消去rs，
将n·vth看做整体变量，得到一个关于n·vth的一元二次方程：
形如：a(n·vth)2+b(n·vth)+c＝0(23)
其中：
故：
(十五)、求出二极管理想因子n后，可以利用式(13)求出串联电阻rs：
(十六)、为求得io，引入另一电压调制点(uc，ic)，该点不能与(ua，ia)、(ub，ib)重合，根据(1)式可推导出(29)式：
式(29)为二元一次方程，结合rsho≈rsh故可得到io与iph的表达式：
(十七)、结合(5)式求出rsh表达式：
步骤二中，根据推导公式，我们可以明确求出光伏组件五参数所需要的点如下：
①短路电流点(0，isc)；
②短路电流附近点(δv，isc+δi)，其中，δv为电压调制模块调制电压值得最小步长，从而
③(ua，ia)点，该点电压ua应靠近uoc，实际ua可取0.8uoc；
④(ub，ib)点，该点电压ub应靠近uoc，且不能等于ua；
⑤(u*，i*)点，该点电压u*不能等于ua和ub，实际u*可取0.7uoc；
⑥(u*，i*)附近点(δv+u*，δi+i*)，其中，δv为电压调制模块调制电压值得最小步长，从而
⑦(uc，ic)点，该点的电压值uc不能等于ua。
根据采样所需的七个点结合式(27)(28)(30)(31)(32)可以求出光伏组件的五参数。
综上所述，在光伏组件i-v曲线缺失开路电压附近i-v输出特性曲线时，依旧可以采用上述七个采样点依据五参数解析公式计算得出光伏组件在该状态下的五参数。七个采样点位置示意图如图2所示。
为验证该方法的准确性与可靠性，根据经典实验数据求出按该方法求出曲线数据与测试数据做对比，结果验证了该方法的可行性，实验数据与仿真数据对比如图3所示。由图3可知该方法的准确性较高，可以较为准确地求出光伏组件五参数，因此本发明所提出方法具有一定实用性。