本发明公开了一种太阳能组件快速MPPT的方法，包括：整点时，改变太阳能组件板外部Buck电路5次，测量5组电压、电流值，确定相应5参数模型；根据5参数模型利用迭代法求解5参数，得到整点的U‑I；对任意时刻，将其中受温度变化影响的3个参数用整点值代，对于另外2个参数，通过测量2组电压电流值确定2个方程，解出2参数的值，得到任意时刻的U‑I；利用功率对电压导数为0确定最大功率点对应的电压电流值，确定Buck电路占空比并调节使得光伏组件板运行在最大功率点处；结合扰动观察法小范围精细调整，跟踪最大功率点。本发明通过直接求解太阳能组件板U‑I方程，求解最大功率点快速确定最大功率点避免电路多次试探跳动。
1.一种太阳能组件快速MPPT方法，其特征在于，包括如下步骤：步骤S1：整点时，改变光伏组件板外部Buck电路5次，测量5组电压、电流值，根据电压情况确定相应的遮蔽情况，从而确定相应的5参数模型；
步骤S2：根据5参数模型和这5组电压、电流值确定5个方程，用迭代法求解5参数，从而得到整点时的U-I关系；
步骤S3：对任意时刻，设定辐照度与上一整点时刻的相同，将5参数中只受温度变化影响的3个参数用整点值代，对于2个受辐照度影响较大的参数，通过测量2组电压、电流值确定2个方程，直接解出2参数的值，从而得到任意时刻的U-I关系；
步骤S4：利用U-I关系式和最大功率点处功率P对电压U的导数为0确定最大功率点对应的电压、电流值，继而确定Buck电路占空比，调节占空比，使得光伏组件板运行在最大功率点处；
步骤S5：结合扰动观察法小范围精细调整，跟踪最大功率点；
所述光伏组件板包括60个电池单体，一个电池单体等效的5参数模型如公式(1)所示：式中：q为电子电量常量，为1.602×10-19C；k为玻尔兹曼常量，为1.381×10-23J/K；A为二极管特性拟合系数，是一个变量，其值在1到2之间；
所述5参数分别是Iph为光生电流、Rsh为等效并联电阻、Is为二极管饱和电流、Rs为等效串联电阻以及 为系数，其中q为电子电量常量，A为二极管特性拟合系数，T为实时的温度，标准工况为辐照度为1000W/m2，电池工作温度为25℃；
其中，光生电流Iph在不同照度和辐照度下的修正量如式(2)所示：式中：Iphref是标准工况下的光生电流；S为实时的辐照度，Sref为标况下的辐照度；CT为温度系数；
二极管饱和电流Is在不同温度下的修正量如公式(3)所示：式中：Isref是标况下的二极管饱和电流，Tref为标况下的温度，Eg为禁带宽度；
等效串联电阻Rs的修正式如公式(4)所示：式中：β＝0.217，Rsref为标准工况下的等效串联电阻，计算可知，辐照度变化，Rs变化不大；
等效并联电阻Rsh的修正式如公式(5)所示：式中：Rshref为标准工况下的等效并联电阻；
所述5参数中，短时间内温度以及辐照度不对参数 Is和Rs产生影响，对参数Iph、Rsh受辐照度影响较大，从而在确定某一特定时刻的五参数值时， Is和Rs参数为整点处的值，只需确定Iph、Rsh的值。
2.根据权利要求1所述的一种太阳能组件快速MPPT方法，其特征在于，所述光伏组件板中，每20个电池单体为1组，每组串联一个旁路二极管，记N＝60，光伏组件板输出电流和电压的关系式如公式(6)所示：考虑遮蔽情况，当一组电池单体中有一个单体或几个单体被遮蔽，电流将直接通过旁路二极管，每组的端电压将是二极管的正向导通电压；
假设一块电池板有N1个硅片没有遮蔽，N2个硅片被遮蔽，从而考虑遮蔽情况下U-I关系式如式(7)所示：式中：UDD为一个二极管的正向导通电压。
3.根据权利要求2所述的一种太阳能组件快速MPPT方法，其特征在于，在整点时求取5个参数的方法是，光伏电池外接buck电路，任意时刻通过改变buck电路占空比，得到同一U-I曲线上不同运行点，在每个整点测量五组U-I值，分别是：(U1，I1)，(U2，I2)，(U3，I3)，(U4，I4)，(U5，I5)通过Newton法迭代确定五个参数，从而确定每个整点的U-I关系式，具体方法如下：首先，由五组值得到五个方程，如式(8)所示：
记关于这五个参数的五个函数分别为式(9)所示：即求多元非线性方程组(9)的零点，其Jacobi矩阵如下：该矩阵非奇异，Newton法局部线性收敛；
记 F＝(F1,F2,F3,F4,F5)′则迭代公式如式(10)所示：
在任意时刻求取5个参数的方法是， Is和Rs参数取最近整点处的值，通过改变外部电路的工作状态2次，测量两组值：(U1，I1)，(U2，I2)
从而可以求出Iph和Rsh，从而确定了任意时刻的五参数值，从而确定任意时刻形如式(7)的U-I关系式。
4.根据权利要求3所述的一种太阳能组件快速MPPT方法，其特征在于，求最大功率点的方法的依据是，一定光照条件和温度条件下I-U曲线和P-U曲线是确定的，功率表达式如式(12)所示：P＝UI (12)
最大功率点处对应输出功率最大值点，P对U的导数为0，即式(13)所示：此时输出功率是最大输出功率；
所述求最大功率点的方法是，将U-I曲线式(7)变形得到式(14)如下：式(12)和式(14)两边对U求导化简可得式(15)如下：将式(15)与式(14)联立，由式(15)求得最大功率点(UMP，IMP)。
5.根据权利要求4所述的太阳能组件快速MPPT方法，其特征还包括，所述步骤S5中，负载模型恒功率、恒电压、恒电阻的组合通过Buck电路与光伏组件板连接，由Buck电路前后功率相等，得系统端电压电流与占空比的关系式如式(16)所示：求出最大功率点对应的U-I值后，求出相应的占空比，调节占空比直接将工作点跳到最大功率点附近。
一种太阳能组件快速MPPT方法
技术领域
[0001] 本发明涉及到光伏发电技术领域，具体地，涉及一种太阳能组件快速MPPT方法。
背景技术
[0002] 太阳能是一种十分具有潜力的新能源，光伏发电是当前利用太阳能的主要方式之一。光伏器件的输出功率受辐照度和温度的影响，在外部环境稳定的条件下，光伏器件存在唯一的最大输出功率点。为了在同样的辐照度和电池结温条件下获得尽可能多的电能，需要对最大功率点实时跟踪。
[0003] 早期对光伏阵列MPPT技术的研究主要是定电压跟踪法(CVT)、光伏阵列组合法以及实际测量法。目前光伏阵列的最大功率跟踪(MPPT)方法主要分为基于数学模型的方法、基于扰动的自寻优法和基于智能技术的方法。
[0004] 基于数学模型的方法是以建立优化的数学模型为出发点，构造求解方法及光伏阵列特性曲线，从而得出光伏阵列的最大功率输出，所以光伏电池的等效电路模型及各种参数的正确性是需要着重考虑的。
[0005] 利用五参数模型时，生产厂家一般只会给出标准工况下，即辐照度为1000W/m2，电池工作温度为25℃(298K)的开路电压UOC，短路电流ISC，最大功率点处的电压VMP和电流IMP值。现有技术中先利用上面的四个值求标况下的五参数，再引入修正量，来确定在一般情况下的五个参数值，但辐照度的实际测量误差可能会很大，可行性不高；且该方法是在假定二极管理想因子不变的情况下确定的，实际二极管理想因子对曲线的影响较大，对所有电池板没有一个普遍适合的值。
[0006] 基于扰动的自寻优法是目前研究最广也是应用较为普遍的控制方法，根据直接测量到的光伏阵列的电压和电流等信息进行最大功率跟踪，括扰动观察法(P&O)、增加电导法和类似P&O法的波动相关控制法(RCC)等。这类方法需要试探性的调整电路的工作状态，逻辑简单，但比较盲目，时间不够快。而智能法包括模糊理论和神经网络等建立在大数据的条件下，目前不易实施。
发明内容
[0007] 针对上述技术问题，本发明的目的在于提供一种太阳能组件快速MPPT方法。
[0008] 为实现上述目的，本发明是根据以下技术方案实现的：
[0009] 一种太阳能组件快速MPPT方法，包括如下步骤：
[0010] 步骤S1：整点时，改变太阳能组件板外部Buck电路5次，测量5组电压、电流值，根据电压情况确定相应的遮蔽情况，从而确定相应的5参数模型；
[0011] 步骤S2：根据5参数模型和这5组电压、电流值确定5个方程，用迭代法求解5参数，从而得到整点时的U-I关系；
[0012] 步骤S3：对任意时刻，设定辐照度与上一整点时刻的相同，将5参数中只受温度变化影响的3个参数用整点值代，对于2个受辐照度影响较大的参数，通过测量2组电压、电流值确定2个方程，直接解出2参数的值，从而得到任意时刻的U-I关系；
[0013] 步骤S4：利用U-I关系式和最大功率点处功率P对电压U的导数为0确定最大功率点对应的电压、电流值，继而确定Buck电路占空比，调节占空比，使得光伏组件板运行在最大功率点处；
[0014] 步骤S5：结合扰动观察法小范围精细调整，跟踪最大功率点。
[0015] 上述技术方案中，所述光伏电池组件板包括60个电池单体，一个电池单体等效的5参数模型如公式(1)所示：
[0016]
[0017] 式中：q为电子电量常量，为1.602×10-19C；k为玻尔兹曼常量，为1.381×10-23J/K；A为二极管特性拟合系数，是一个变量，其值在1到2之间；
[0018] 所述5参数分别是Iph为光生电流、Rsh为等效并联电阻、Is为二极管饱和电流、Rs为等效串联电阻以及 为系数，其中q为电子电量常量，A为二极管特性拟合系数，T为实时的温度，标准工况为辐照度为1000W/m2，电池工作温度为25℃；
[0019] 其中，光生电流Iph在不同照度和辐照度下的修正量如式(2)所示：
[0020]
[0021] 式中：Iphref是标准工况下的光生电流；S为实时的辐照度，Sref为标况下的辐照度；CT为温度系数；
[0022] 二极管饱和电流Is在不同温度下的修正量如公式(3)所示：
[0023]
[0024] 式中：Isref是标况下的二极管饱和电流，Tref为标况下的温度，Eg为禁带宽度；
[0025] 等效串联电阻Rs的修正式如公式(4)所示：
[0026]
[0027] 式中：β＝0.217，Rsref为标准工况下的等效串联电阻，计算可知，辐照度变化，Rs变化不大；
[0028] 等效并联电阻Rsh的修正式如公式(5)所示：
[0029]
[0030] 式中：Rshref为标准工况下的等效并联电阻。
[0031] 所述5参数中，短时间内温度以及辐照度不对参数 Is和Rs产生影响，对参数Iph、Rsh受辐照度影响较大，从而在确定某一特定时刻的五参数值时， Is和Rs参数为整点处的值，只需确定Iph、Rsh的值。
[0032] 上述技术方案中，所述太阳能组件板中，每20个电池单体为1组，每组串联一个旁路二极管，记N＝60，太阳能组件板输出电流和电压的关系式如公式(6)所示：
[0033]
[0034] 所述考虑遮蔽情况，当一组电池单体中有一个单体或几个单体被遮蔽，电流将直接通过旁路二极管，每组的端电压将是二极管的正向导通电压；
[0035] 假设一块电池板有N1个硅片没有遮蔽，N2个硅片被遮蔽，从而考虑遮蔽情况下U-I关系式如式(7)所示：
[0036]
[0037] UDD是一个二极管的正向导通电压。
[0038] 上述技术方案中，所述的在整点时求取5个参数的方法是，光伏电池外接buck电路，任意时刻通过改变buck电路占空比，得到同一U-I曲线上不同运行点，在每个整点测量五组U-I值，分别是：
[0039] (U1，I1)，(U2，I2)，(U3，I3)，(U4，I4)，(U5，I5)
[0040] 通过Newton法迭代确定五个参数，从而确定每个整点的U-I关系式，具体方法如下：
[0041] 首先，由五组值得到五个方程，如式(8)所示：
[0042]
[0043] 记关于这五个参数的五个函数分别为式(9)所示：
[0044]
[0045] 即求多元非线性方程组(9)的零点，其Jacobi矩阵如下：
[0046]
[0047] 该矩阵非奇异，Newton法局部线性收敛；
[0048] 记 F＝(F1,F2,F3,F4,F5)′
[0049] 则迭代公式如式(10)所示：
[0050]
[0051] 取上个整点测得的值作为初始值进行迭代；
[0052] 所述的在任意时刻求取5个参数的方法是， Is和Rs参数取最近整点处的值，通过改变外部电路的工作状态2次，测量两组值：
[0053] (U1，I1)，(U2，I2)
[0054] 得到相应的方程如式(11)所示：
[0055]
[0056] 从而可以求出Iph和Rsh，从而确定了任意时刻的五参数值，从而确定任意时刻形如式(7)的U-I关系式。
[0057] 上述技术方案中，所述求最大功率点的方法的依据是，一定光照条件和温度条件下I-U曲线和P-U曲线是确定的，功率表达式如式(12)所示：P＝UI (12)
[0058] 最大功率点处对应输出功率最大值点，P对U的导数为0，即式(13)所示：
[0059] 此时输出功率是最大输出功率；
[0060] 所述求最大功率点的方法是，将U-I曲线式(7)变形得到式(14)如下：
[0061]
[0062] 式(12)和式(14)两边对U求导化简可得到式(15)如下：
[0063]
[0064] 将式(15)与式(14)联立，由式(15)求得最大功率点(UMP，IMP)；
[0065] 上述技术方案中，所述步骤S5中，负载模型恒功率、恒电压、恒电阻的组合通过Buck电路与光伏组件板连接，由Buck电路前后功率相等，得该系统端电压电流与占空比的关系式如式(16)所示：
[0066]
[0067] 求出最大功率点对应的U-I值后，求出相应的占空比，调节占空比直接将工作点跳到最大功率点附近。
[0068] 本发明与现有技术相比，具有如下有益效果：
[0069] 本发明基于光伏电池组件板的5参数模型，通过调节组件板外部的Buck电路获得电压电流的不同测量值，利用测量值计算任意时刻的5参数的值，求取任意时刻的U-I方程，利用U-I方程直接求解最大功率点。与现有的采用跟踪最大功率点的方法相比，所述方法不需要对电路做多次的调节，一步到位，具有快速、简单易操作的优点。
附图说明
[0070] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还根据这些附图获得其它附图。
[0071] 图1为光伏电池单体等效电路；
[0072] 图2为光伏电池组件板示意图及外部Buck电路等效电路；
[0073] 图3为示意图光伏阵列I-U曲线和P-U曲线；
[0074] 图4为外接Buck电路及负载等效电路图。
具体实施方式
[0075] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。
[0076] 一种太阳能组件快速MPPT方法，包括如下步骤：
[0077] 步骤S1：整点时，改变太阳能组件板外部Buck电路5次，测量5组电压、电流值，根据电压情况确定相应的遮蔽情况，从而确定相应的5参数模型；
[0078] 步骤S2：根据5参数模型和这5组电压、电流值确定5个方程，用迭代法求解5参数，从而得到整点时的U-I关系；
[0079] 步骤S3：对任意时刻，设定辐照度与上一整点时刻的相同，将5参数中只受温度变化影响的3个参数用整点值代，对于2个受辐照度影响较大的参数，通过测量2组电压、电流值确定2个方程，直接解出2参数的值，从而得到任意时刻的U-I关系；
[0080] 步骤S4：利用U-I关系式和最大功率点处功率P对电压U的导数为0确定最大功率点对应的电压、电流值，继而确定Buck电路占空比，调节占空比，使得光伏组件板运行在最大功率点处；
[0081] 步骤S5：结合扰动观察法小范围精细调整，跟踪最大功率点。
[0082] 上述技术方案中，所述光伏电池组件板包括60个电池单体，一个电池单体等效如图1所示，5参数模型如公式(1)所示：
[0083]
[0084] 式中：q为电子电量常量，为1.602×10-19C；k为玻尔兹曼常量，为1.381×10-23J/K；A为二极管特性拟合系数，是一个变量，其值在1到2之间；
[0085] 所述5参数分别是Iph为光生电流、Rsh为等效并联电阻、Is为二极管饱和电流、Rs为等效串联电阻以及 为系数，其中q为电子电量常量，A为二极管特性拟合系数，T为实时的温度，标准工况为辐照度为1000W/m2，电池工作温度为25℃；
[0086] 其中，光生电流Iph在不同照度和辐照度下的修正量如式(2)所示：
[0087]
[0088] 式中：Iphref是标准工况下的光生电流；S为实时的辐照度，Sref为标况下的辐照度；CT为温度系数；
[0089] 二极管饱和电流Is在不同温度下的修正量如公式(3)所示：
[0090]
[0091] 式中：Isref是标况下的二极管饱和电流，Tref为标况下的温度，Eg为禁带宽度；
[0092] 等效串联电阻Rs的修正式如公式(4)所示：
[0093]
[0094] 式中：β＝0.217，Rsref为标准工况下的等效串联电阻，计算可知，辐照度变化，Rs变化不大；
[0095] 等效并联电阻Rsh的修正式如公式(5)所示：
[0096]
[0097] 式中：Rshref为标准工况下的等效并联电阻；
[0098] 所述5参数中，短时间内温度以及辐照度不对参数 Is和Rs产生影响，对参数Iph、Rsh受辐照度影响较大，从而在确定某一特定时刻的五参数值时， Is和Rs参数为整点处的值，只需确定Iph、Rsh的值。
[0099] 所述太阳能组件板中，每20个电池单体为1组，每组串联一个旁路二极管，如图2所示，记N＝60，太阳能组件板输出电流和电压的关系式如公式(6)所示：
[0100]
[0101] 所述考虑遮蔽情况，当一组电池单体中有一个单体或几个单体被遮蔽，电流将直接通过旁路二极管，每组的端电压将是二极管的正向导通电压；
[0102] 假设一块电池板有N1个硅片没有遮蔽，N2个硅片被遮蔽，从而考虑遮蔽情况下U-I关系式如式(7)所示：
[0103]
[0104] 其中，UDD是一个二极管的正向导通电压。
[0105] 所述的在整点时求取5个参数的方法是，光伏电池外接buck电路，任意时刻通过改变buck电路占空比，得到同一U-I曲线上不同运行点，在每个整点测量五组U-I值，分别是：
[0106] (U1，I1)，(U2，I2)，(U3，I3)，(U4，I4)，(U5，I5)
[0107] 通过Newton法迭代确定五个参数，从而确定每个整点的U-I关系式，具体方法如下：
[0108] 首先，由五组值得到五个方程，如式(8)所示：
[0109]
[0110] 记关于这五个参数的五个函数分别为式(9)所示：
[0111]
[0112] 即求多元非线性方程组(9)的零点，其Jacobi矩阵如下：
[0113]
[0114] 该矩阵非奇异，Newton法局部线性收敛；
[0115] 记 F＝(F1,F2,F3,F4,F5)′，则迭代公式如式(10)所示：
[0116] 取上个整点测得的值作为初始值进行迭代；
[0117] 所述的在任意时刻求取5个参数的方法是， Is和Rs参数取最近整点处的值，通过改变外部电路的工作状态2次，测量两组值：
[0118] (U1，I1)，(U2，I2)
[0119] 得到相应的方程如式(11)所示：
[0120]
[0121] 从而可以求出Iph和Rsh，从而确定了任意时刻的五参数值，从而确定任意时刻形如式(7)的U-I关系式。
[0122] 上述技术方案中，所述求最大功率点的方法的依据是，一定光照条件和温度条件下I-U曲线和P-U曲线是确定的，如图3所示，功率表达式如式(12)所示：
[0123] P＝UI (12)
[0124] 最大功率点处对应输出功率最大值点，P对U的导数为0，即式(13)所示：
[0125] 此时输出功率是最大输出功率；
[0126] 所述求最大功率点的方法是，将U-I曲线式(7)变形得到式(14)如下：
[0127]
[0128] 式(12)和式(14)两边对U求导可得化简可得式(15)如下：
[0129]
[0130] 将式(15)与式(14)联立，由式(15)求得最大功率点(UMP，IMP)；
[0131] 上述技术方案中，所述步骤S5中，负载模型恒功率、恒电压、恒电阻的组合通过Buck电路与光伏组件板连接，由Buck电路前后功率相等，如图4所示，得该系统端电压电流与占空比的关系式如式(16)所示：
[0132]
[0133] 求出最大功率点对应的U-I值后，求出相应的占空比，调节占空比直接将工作点跳到最大功率点附近。。
[0134] 本发明的一个具体实施例为:首先，对外接Buck电路的光伏组件板添加实测装置，便于实时反馈组件板的端电流和端电压。利用下式
[0135]
[0136] 设计控制电路。在整点时，通过改变阳能组件板外部Buck电路占空比，运行点跳5次，反馈每次的电压、电流值给控制芯片，芯片根据电压情况判断遮蔽情况，如果U为正常值，则取N1＝60,N2＝0；当有硅片被遮时，如果U的值变化为原来的三分之一左右时，取N1＝40,N2＝20；如果U的值变化为原来的三分之二左右时，取N1＝20,N2＝40；如果U非常小，电池板不工作，无需MPPT。
[0137] 判断出相应的U-I五参数模型后，用迭代法求解这五组值得到的方程组：
[0138] 雅克比矩阵为：
[0139] 迭代方程为：
[0140] 确定了整点时的5参数模型。
[0141] 任意时刻，将 Is和Rs用前一整点值代，改变太阳能组件板外部Buck电路占空比，运行点跳2次，反馈每次的电压、电流值给控制芯片，同整点方法一样先判断考虑遮蔽后的相应模型，用迭代法求解这2组值得到的方程组，得到Iph、Rsh：
[0142]
[0143] 从而确定了任意时刻的5参数模型。
[0144] 得到任意时刻的5参数模型后，利用U-I关系式和最大功率点处功率P对电压U的导数为0，得到关于电压电流的两个方程，
[0145]
[0146]
[0147] 用不动点迭代法用迭代法求解这两个方程，若记
[0148] I＝y
[0149] 则迭代格式为：
[0150] 约束条件是：
[0151]
[0152] 取上一时刻的最大功率点为迭代初值，代入迭代方程，再将x(k+1)带回约束条件，求出y(k+1)。利用U，I大于0这个条件舍取，若y(k+1)有两个值，则比较相应的功率值，取对应功率较大的一组值。迭代至误差在允许范围内，相应的U-I即最大功率处的U-I。
[0153] 最后由最大功率点对应的电压电流值，根据
[0154]
[0155] 确定Buck电路占空比，调节占空比，使得光伏组件板运行在最大功率点处。
[0156] 本发明结合扰动观察法小范围精细调整，跟踪最大功率点。
[0157] 虽然，上文中已经用一般性说明及具体实施例对本发明作了详尽的描述，但在本发明基础上，对之作一些修改或改进，这对本领域技术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的这些修改或改进，均属于本发明要求保护的范围。
