CN102981549A - 一种光伏最大功率点实时跟踪预测控制方法 - Google Patents
一种光伏最大功率点实时跟踪预测控制方法 Download PDFInfo
- Publication number
- CN102981549A CN102981549A CN2012105146405A CN201210514640A CN102981549A CN 102981549 A CN102981549 A CN 102981549A CN 2012105146405 A CN2012105146405 A CN 2012105146405A CN 201210514640 A CN201210514640 A CN 201210514640A CN 102981549 A CN102981549 A CN 102981549A
- Authority
- CN
- China
- Prior art keywords
- photovoltaic
- real
- power point
- time
- prime
- 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.)
- Granted
Links
Images
Classifications
-
- 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
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
Landscapes
- Control Of Electrical Variables (AREA)
Abstract
本发明提供一种光伏最大功率点实时跟踪预测控制方法,S1,建立环境因子与光伏控制器中光伏最大功率点函数模型,并建立对环境因子与光伏控制器中光伏最大功率点函数模型中的实时输出电流进行控制的光伏控制器输出电流和电压的数学模型;S2,建立光伏变换器建立数学模型、状态空间模型,为控制输出电流的光伏变换器建立约束条件;S3,建立实时跟踪预测的目标函数,并建立用于确定目标函数的性能指标;S4,将所述目标函数在所述约束条件下求解,获得用于跟踪光伏最大功率点的最优控制序列。本发明可以实现光伏系统在迅速变化的外部环境条件下实时的模型预测控制,从而提高了对最大功率点的预测和跟踪能力。
Description
技术领域
本发明涉及光伏功率点跟踪的优化控制方法技术领域,特别为涉及一种光伏最大功率点实时跟踪预测控制方法。
背景技术
随着环境问题的突出和传统能源的枯竭,可再生能源如太阳能、风能、潮汐能越来越受到人们的关注,其中,太阳能作为方便、环保、高效的新型可再生能源,有着广阔的应用前景,中国的光伏产业经历了大幅度的升级与扩张。2009年7月21日,财政部、科技部、国家能源局联合发布了《关于实施金太阳示范工程的通知》,决定综合采取财政补助、科技支持和市场拉动方式,加快国内光伏发电的产业化和规模化发展。三部委计划在近年内,采取财政补助方式支持不低于500兆瓦的光伏发电示范项目,据估算,国家将为此投入约100亿元财政资金。光伏发电技术对新型能源技术、智能电网技术的开发与应用有着显著的推动作用。通过引进消化吸收国外先进技术和自主研发,我国已经掌握了世界先进的光伏器件制造技术,但在提高光伏发电效率、友好并网、低电压穿越等方面的研究起步相对较晚,由于我国能源市场的巨大需求和节能减排的要求,迫切需要更加高效、稳定的光伏发电能力。
经对现有技术的公开文献检索发现,N.Femia,D.Granozio,and G.Petrone.Predictive andAdaptive MPPT Perturb and Observe Method.IEEE Transactions on Aerospace and ElectronicSystems,2007,43(3):934–949.(预测与自适应的光伏扰动观察最大功率点跟踪方法,国际期刊:IEEE期刊,航空与电子系统,2007,43(3):934–949),虽然作者将自适应预测控制用于寻找最大功率点,形成了基于电压和功率等自适应参数的模型预测控制方法。但在实际运行中,环境的突变对最大功率点的搜索影响很大,甚至会出现搜索偏离、运行效率降低的状况,而且无法体现外部环境参数变化与最大功率点的内在联系。
发明内容
鉴于以上所述现有技术的缺点,本发明的目的在于提供一种光伏最大功率点实时跟踪预测控制方法,用于在多变环境下搜索光伏最大功率点。
为实现上述目的及其他相关目的,本发明提供一种光伏最大功率点实时跟踪预测控制方法,包括以下步骤:
S1,建立环境因子与光伏控制器中光伏最大功率点函数模型:
p=v×i,pmax=max{p1,p2,…,pn}:
其中:p为输出功率, v为实时输出电压;i为实时输出电流;pmax为输出最大功率;p1,p2,…,pn分别为第一个、第二个、第n个输出功率;C1,C2,C3,C4,C5分别为常数;G为实时测量的光照;T为实时测量的环境温度;
S2,为控制输出电流的光伏变换器建立约束条件:
S3,建立实时跟踪预测的目标函数:J=γTHγ+2fTγ;
其中:J为目标函数,H和f为预测输出电流的系数矩阵,γ为待求解的最优序列,H=ΩTQyΩ+Qm,f=2ζTQyζ-2dTQm,ζ为预测输出模型中状态变量的系数矩阵,Ω为预测输出模型中序列γ的系数矩阵,d为变换器的占空比初始值,Qy和Qm均为正定权矩阵,T为矩阵转置符号;
S4,将所述目标函数在所述约束条件下求解,获得用于跟踪光伏最大功率点的最优控制序列。
可选地,在步骤S1中还包括:对光伏控制器输出电流和电压建立的数学模型以对环境因子与光伏控制器中光伏最大功率点函数模型中的实时输出电流进行控制,所述数学模型为:
其中:I为输出电流;V为输出电压;Ipv为光照产生电流,I0为光伏电池内部等效二极管P-N结反向饱和电流,a为P-N结曲线常数;Rs为串联等效电阻;Rp为并联等效电阻;Vos为开路电压;Is为短路电流;G为实时光照;Gs为1000w/m2的标准环境参数;κ为短路电流温度系数;ΔT为P-N结温度实时变化差;η为开路电压温度系数。
可选地,在步骤S2中还包括根据光伏控制器的等效电路和基尔霍夫定律为控制输出电流的光伏变换器建立数学模型:
1)光伏变换器在断开状态的数学模型为:
2)光伏变换器在闭合状态的数学模型为:
其中:i为输出电流;t为电流输出的当前时间;为输出电流i对电流输出的当前时间t的导数;vpv为光伏阵列两端电压;vo光伏阵列负载电压,为光伏阵列负载电压vo对电流输出的当前时间t的导数;R为负载电阻;C为光伏变换器的并联电容;L为光伏变换器的串联电感。
可选地,还包括:根据所述光伏变换器的数学模型建立状态空间模型:
1)光伏变换器在断开状态的状态空间模型为:
2)光伏变换器在闭合状态的状态空间模型为:
其中: C1=(10);
A1、A2,B1,B2,C1为系数矩阵,x为状态变量,x=(i,v0)T;为状态变量的导数;y为输出电流;i为即时电流;vo光伏阵列负载电压;R为负载电阻;C为光伏变换器的并联电容;L为光伏变换器的串联电感;T为矩阵转置符号。
可选地,将光伏变换器在断开和闭合状态的状态空间模型为转化为线性函数:
x(k+1)=A'×x(k)+B'×d(k)+R×δ(k)+G×z(k);
其中:k为时间状态参数;x(k)为k时刻状态变量;x(k+1)为k+1时刻状态变量;d(k)为k时刻变换器的占空比;δ(k)为k时刻理想开关的开关变量;z(k)为k时刻辅助变量;A'、B'、R、G为系数矩阵。
可选地,根据环境因子与光伏控制器中光伏最大功率点函数模型、所述光伏变换器的空间状态方程、预测时域和控制时域,建立用于确定目标函数的性能指标:
其中:L为性能指标,k,j为时间状态参数;Q=QT≥0,Qy和Qm分别为正定权矩阵;i(k+j/k)为在k时刻对k+j时刻的预测输出电流;i*(j)为k+j时刻的参考电流;Δd(k+j-1/k)为在k时刻对k+j时刻和k+j-1时刻预测占空比之差;P为预测时域,M为控制时域。
可选地,在步骤S3中,建立实时跟踪预测的目标函数的过程为:
1)建立用于获取预测输出的预测输出函数模型:
y(k+j/k)=C×(A',A′2,…,A′P)T×x(k)+Ω×γ(k+j/k);
2)将获取的预测输出带入到所述性能指标,获得目标函数
其中:γ(k+j/k)={d(i)T,δ(i)T,z(i)T},i∈{k,...,k+M-1};
y(k+j/k)为在k时刻对k+j时刻的预测输出;Ω为总序列的系数矩阵;C为输出矩阵;A'、B'为状态变量x(k+1)中的系数矩阵;d(i)T为占空比预测值序列;δ(i)T为开关变量预测值序列;z(i)T为辅助变量预测值序列;γ(k+j/k)为占空比、开关变量和辅助变量的总序列;P为预测时域,M为控制时域;R为负载电阻;G为实时光照;T为矩阵转置符号。
可选地,在步骤S4中,将所述目标函数在所述约束条件下求解具体包括:
如上所述,本发明的一种光伏最大功率点实时跟踪预测控制方法,具有以下有益效果:
本发明可以实现光伏系统在迅速变化的外部环境条件下实时的模型预测控制,从而提高了对最大功率点的预测和跟踪能力,该方法适用于独立式、分布式、孤岛或并网的光伏发电系统最大功率点跟踪控制。
附图说明
图1显示为本发明的光伏最大功率点实时跟踪预测控制方法的流程图。
图2显示为本发明的光伏最大功率点实时跟踪预测控制方法中的具有环境影响因子的最大功率点输出模型效果图。
图3显示为本发明的光伏最大功率点实时跟踪预测控制方法中的输出电压电流曲线和输出电压功率曲线示意图。
图4显示为本发明的光伏最大功率点实时跟踪预测控制方法的控制回路图。
图5显示为本发明的光伏最大功率点实时跟踪预测控制方法中跟踪输出最大功率变化曲线图。
图6显示为本发明的光伏最大功率点实时跟踪预测控制方法在晴天实施时的电流-功率曲线效果图。
图7显示为本发明的光伏最大功率点实时跟踪预测控制方法在晴天实施时的电压-电流曲线效果图。
图8显示为本发明的光伏最大功率点实时跟踪预测控制方法在阴天实施时的电流-功率曲线效果图。
图9显示为本发明的光伏最大功率点实时跟踪预测控制方法在阴天实施时的电压-电流曲线效果图。
图10显示为本发明的光伏最大功率点实时跟踪预测控制方法在雨天实施时的电流-功率曲线效果图。
图11显示为本发明的光伏最大功率点实时跟踪预测控制方法在雨天实施时的电压-电流曲线效果图。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。
需要说明的是,本实施例中所提供的图示仅以示意方式说明本发明的基本构想,遂图式中仅显示与本发明中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制,其实际实施时各组件的型态、数量及比例可为一种随意的改变,且其组件布局型态也可能更为复杂。
本发明提供一种光伏最大功率点实时跟踪预测控制方法,用于在多变环境下搜索光伏最大功率点。
本发明是通过以下技术方案实现的,本发明根据光伏阵列输出电流电压的非线性特性,考虑了引起最大功率点变化的外部环境因素,给出环境参数与最大功率点函数模型,综合被控过程的安全操作约束边界,定义了直流电流控制的优化目标,并求解出在约束条件下性能指标函数的上界和下界,最后通过求解混合整数二次优化问题得出控制律。具体包括具有环境参数因子的最大功率点建模、光伏变换器建模以及安全运行的系统约束和具有约束的动态模型预测控制方法三个步骤,其中具有环境参数因子的最大功率点建模和具有系统约束与优化目标的动态预测控制方法为本发明的创新之处。
以下将详细阐述本发明的一种光伏最大功率点实时跟踪预测控制方法的原理及实施方式,使本领域技术人员不需要创造性劳动即可理解本发明的一种光伏最大功率点实时跟踪预测控制方法。
请参阅图1,显示为本发明的一种光伏最大功率点实时跟踪预测控制方法的流程图。如图1所示,本发明提供一种光伏最大功率点实时跟踪预测控制方法具体包括以下步骤:
S1,建立环境因子与光伏控制器中光伏最大功率点函数模型。
S2,为控制输出电流的光伏变换器建立约束条件。
S3,建立实时跟踪预测的目标函数。
S4,将所述目标函数在所述约束条件下求解,获得用于跟踪光伏最大功率点的最优控制序列。
下面详细对各步骤的建立和实现方法进行说明。
【一】S1,建立环境因子与光伏控制器中光伏最大功率点函数模型。
请参阅图2,显示为本发明的光伏最大功率点实时跟踪预测控制方法中的具有环境影响因子的最大功率点输出模型效果图。
根据外部环境参数如光照、温度对光伏阵列最大功率点的影响,通过曲面拟合建立环境因子与最大功率点的关系并给出函数模型。
环境因子与光伏控制器中光伏最大功率点函数模型为:
p=v×i,pmax=max{p1,p2,…,pn}:
其中:p为输出功率, v为实时输出电压;i为实时输出电流;pmax为输出最大功率;p1,p2,…,pn分别为第一个、第二个、第n个输出功率;C1,C2,C3,C4,C5分别为常数,其中,C1=-0.1754,C2=0.1132,C3=4.0779E-5,C4=0.0032,C5=0.000235;G为实时测量的光照;T为实时测量的环境温度。
在本步骤中,还包括对光伏控制器输出电流和电压建立的数学模型以对环境因子与光伏控制器中光伏最大功率点函数模型中的实时输出电流进行控制。
光伏阵列中半导体二极管暴露于可见光,在P-N节的两侧产生电压,连接到P型材料和N型材料上的电极之间产生电流通过,将光能转化为电能。除了光照强度G,外界温度T的变化对光伏阵列的输出电压V和电流I也有很大影响,一般来说,光照强度的变化率与电压电流的变化率成正比,温度的变化率与电压电流变化率成反比。当光照强度和温度均为变量时,光伏阵列可以看作非线性的直流电源,其输出电压电流V-I关系特性的数学模型可以如下表示:
其中:I为输出电流;V为输出电压;Ipv为光照产生电流,I0为光伏电池内部等效二极管P-N结反向饱和电流,a为P-N结曲线常数(一般取3-5);Rs为串联等效电阻;Rp为并联等效电阻;Vos为开路电压;Is为短路电流;G为实时光照;Gs为1000w/m2的标准环境参数;κ为短路电流温度系数;ΔT为P-N结温度实时变化差;η为开路电压温度系数。
根据图3所示,光伏阵列的电压电流曲线和电压功率曲线呈现出非线性特性,为了使光伏阵列输出能量最大化,应该让光伏电池始终保持工作在唯一的最大功率点,即电压电流乘积最大点。然而在外部环境不断变化的条件下,电压电流输出曲线不断变化,导致最大功率点始终处于变化状态。另外,开路电压和短路电流为反映光伏阵列输出的两个重要变量,决定了输出曲线和最大功率点的位置,这两个变量直接受到外部环境的影响,不同的环境将导致不同的最大功率点。
对环境参数因子和输出电压电流进行数据集合,搜寻对应的最大功率点,利用曲面拟合技术,建立环境因子与最大功率点函数模型。
【二】S2,为控制输出电流的光伏变换器建立约束条件:
光伏变换器为调节光伏阵列输出电流的执行器,用于控制输出电流,通过改变开关环节的占空比,使光伏输出电流达到目标值。在本实施例中,在建立约束条件以前,先根据光伏控制器的等效电路和基尔霍夫定律建立光伏变换器的数学模型:
1)光伏变换器在断开状态的数学模型为:
2)光伏变换器在闭合状态的数学模型为:
其中:i为输出电流,t为电流输出的当前时间,为输出电流i对电流输出的当前时间t的导数;vpv为光伏阵列两端电压;vo光伏阵列负载电压,为光伏阵列负载电压vo对电流输出的当前时间t的导数;R为负载电阻;C为光伏变换器的并联电容;L为光伏变换器的串联电感。
开关频率与脉宽调制信号的频率相等,则在一个PWM周期中能够进行开关器断开和闭合两种状态的切换。根据所述光伏变换器的数学模型建立状态空间模型:
1)光伏变换器在断开状态的状态空间模型为:
2)光伏变换器在闭合状态的状态空间模型为:
其中: C1=(10);
A1、A2,B1,B2,C1为系数矩阵,x为状态变量,x=(i,v0)T;为状态变量的导数;y为模型输出,在本模型中即输出电流;i为即时电流;vo光伏阵列负载电压;R为负载电阻;C为光伏变换器的并联电容;L为光伏变换器的串联电感;T为矩阵转置符号。
为了将此状态切换模型转化为能应用优化控制方法的线性离散数学模型,对一个操作周期进行分段,每一个子周期时间长度为τs=Ts/N,N为分段数(一般N取2),ξ(n)为变换器第n段子周期的状态,则对于一个周期T=[k,k+1],有ξ(0)=x(k),ξ(N)=x(k+1);δn为第n段子周期的开关变量,δn=1表示开关闭合,δn=0表示开关断开,所以变换器状态更新可以表示为:
其中,α=Nd(k)-n,β=1-Nd(k)+n,n∈{1,2,…N},表示第n个子周期; ψ1和ψ2分别表示A1,A2,B1和B2的离散形式;d(k)为k时刻占空比;第三种状态为第一种状态和第二种状态的过渡状态,为两种状态的加权形式。对于上述形式分段离散模型,利用混合整数二次优化方法将分段函数转化为线性函数,通过引进辅助变量z(k),状态变量x(k+1)可以表示为以下形式:
x(k+1)=A'×x(k)+B'×d(k)+R×δ(k)+G×z(k);
式中x(k)为k时刻状态变量;d(k)为k时刻占空比;δ(k)为k时刻开关变量;z(k)为k时刻辅助变量;A',B',R,G为x(k)状态代入ξ(n+1)递推公式得到的系数矩阵。引进的开关变量和辅助变量将分段函数形式转化成线性形式,其中所有逻辑动态的上下界写成不等式形式:
E1×d(k)+E2×δ(k)+E3×z(k)≤E4×x(k)+E5;
另外,变换器输出电流i(k)过低会影响工作效率,过高会影响变换器安全操作值,占空比d(k)作为系统控制输入也有上下界,其约束分别为:
综上获得所需的约束条件。
【三】S3,建立实时跟踪预测的目标函数。
所述目标函数为:J=γTHγ+2fTγ;其中:J为目标函数,H和f为预测输出电流的系数矩阵,具体地,H和f为预测输出电流i(k+j/k)代入性能指标所推得系数矩阵;γ为待求解的最优序列,H=ΩTQyΩ+Qm,f=2ζTQyζ-2dTQm,ζ为预测输出模型中状态变量x(k)的系数矩阵,Ω为预测输出模型中序列γ的系数矩阵,d为光伏变换器的占空比初始值,通常取0.5;Qy和Qm均为正定权矩阵,T为矩阵转置符号;
建立目标函数的过程为:
1)根据环境因子与光伏控制器中光伏最大功率点函数模型、所述光伏变换器的空间状态方程、预测时域和控制时域,建立用于确定目标函数的性能指标:
其中:L为性能指标,k,j为时间状态参数;Q=QT≥0,Qy和Qm分别为正定权矩阵;i(k+j/k)为在k时刻对k+j 时刻的预测输出电流;i*(j)为k+j时刻的参考电流;Δd(k+j-1/k)为在k时刻对k+j时刻和k+j-1时刻预测占空比之差;P为预测时域,M为控制时域。
性能指标表示的物理意义为使输出电流与参考电流的偏差最小,达到跟踪的目的;并且输入变量的变化率最小,避免开关器频繁变动引起能量损失。
结合步骤S2中的状态变量x(k+1)表达式和输出模型,预测输出y(k+j/k)可以表示为:
y(k+j/k)=C·(A',A′2,…,A′P)T·x(k)+Ω·γ(k+j/k);
其中γ(k+j/k)={d(i)T,δ(i)T,z(i)T},i∈{k,...,k+M-1};
y(k+j/k)为在k时刻对k+j时刻的预测输出;Ω为总序列的系数矩阵;C为输出矩阵;A'、B'为状态变量x(k+1)中的系数矩阵;d(i)T为占空比预测值序列;δ(i)T为开关变量预测值序列;z(i)T为辅助变量预测值序列;γ(k+j/k)为占空比、开关变量和辅助变量的总序列;P为预测时域,M为控制时域;R为负载电阻;G为实时光照;T为矩阵转置符号。
2)将预测输出带入到所述性能指标,化成混合整数二次规划形式,获得目标函数:
J=γTHγ+2fTγ。
【四】S4,将所述目标函数在所述约束条件下求解,获得用于跟踪光伏最大功率点的最优控制序列。
具体地,将所述目标函数在所述约束条件下求解具体包括:
本发明的实施步骤:
具有在线约束的最大功率点跟踪模型预测控制算法可以总结如下:
a、离线设计
Step 1.针对光伏阵列输出特性,通过实验数据测得不同光照和温度情况下的最大功率点,结合曲面拟合技术,选用麦夸特法和通用全局优化法选取合适的常数C1,C2,C3,C4,C5,得到环境因子与光伏控制器中光伏最大功率点函数模型。
Step 2.为了得到参考电流,需实时测量光伏阵列输入端电压vpv并由p=v×i计算得i,而在测量电压时会出现储能元件的充放电效应导致电压波动,需对电压的跳变进行限制,所以针对变量Δvpv设置变化上界ξ,取ξ小于1V。
Step 3.选择恰当的MPC参数预测时域P和控制时域M。
b、在线设计
Step 1.实时测量外部环境参数T、G和光伏阵列输入端电压vpv,通过具有环境影响因子的最大功率点模型,计算出此时的最大功率点和最大功率点处的输出电流。
Step 2.在优化时域[t,t+N]内,求解在不同约束边界条件下的目标函数值J,然后通过解混合整数二次规划问题,得到控制序列γ;
Step 3.t=t+1,优化时域移动到下一个区间[t+1,t+N+1],使用新的系统输出测量值重复Step 1得到γt+1更新控制序列,不断进行滚动优化,跟踪最大功率点pmax。
结合具体实例进一步说明本发明的实现效果。
图4为本发明实施的控制回路图。在本实施例中,将本发明中提出的方法在simulink电力电子仿真软件联合低压光伏发电系统进行了实施,成功地实现了光伏阵列在不同天气环境下的最大功率点的自动跟踪控制。运行过程分三种情况,分别为晴天,光照从200w/m2变化到1125w/m2,温度从29.5℃变化到33.7℃;阴天,光照从180w/m2变化到750w/m2,温度从18.3℃变化到24.1℃;雨天,光照从125w/m2变化到225w/m2,温度从14.4℃变化到18.5℃。针对以上不同情况进行仿真和实验,每隔2min进行一次环境数据采样和输入端电压采样,若电压变化量超过限制ξ,则可能是储能元件充放电引起的电压波动导致采样不准,需重新进行电压采样。将计算得参考电流作为参考值,调节变换器占空比改变光伏阵列工作点。
使用本发明进行了光伏系统最大功率点跟踪控制以及抗干扰的调节测试,都取得了较好的控制效果,并通过simulink系统的实时趋势可以将控制效果反映出来。图5为光伏阵列在不同环境下跟踪输出最大功率变化曲线图。图6和图7显示为光伏阵列在阴天的天气情况下电流-功率曲线效果图以及电压-电流曲线效果图;图8和图9显示为光伏阵列在晴天的天气情况下电流-功率曲线效果图以及电压-电流曲线效果图;图10和图11显示为光伏阵列在雨天的天气情况下电流-功率曲线效果图以及电压-电流曲线效果图。从图6至图10可以得出,在外部环境变化较大的情况下,本发明的方法能够较好完成系统的控制,使光伏阵列保持工作在最大功率点处,完全达到了提高发电效率、稳定运行的要求。
综上所述,本发明的一种光伏最大功率点实时跟踪预测控制方法,达到了以下有益效果:
本发明可以实现光伏系统在迅速变化的外部环境条件下实时的模型预测控制,从而提高了对最大功率点的预测和跟踪能力,该方法适用于独立式、分布式、孤岛或并网的光伏发电系统最大功率点跟踪控制。
所以,本发明有效克服了现有技术中的种种缺点而具高度产业利用价值。
上述实施例仅例示性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。
Claims (8)
1.一种光伏最大功率点实时跟踪预测控制方法,其特征在于,包括以下步骤:
S1,建立环境因子与光伏控制器中光伏最大功率点函数模型:
p=v×i,pmax=max{p1,p2,…,pn};
其中:p为输出功率, v为实时输出电压;i为实时输出电流;pmax为输出最大功率;p1,p2,…,pn分别为第一个、第二个、第n个输出功率;C1,C2,C3,C4,C5分别为常数;G为实时测量的光照;T为实时测量的环境温度;
S2,为控制输出电流的光伏变换器建立约束条件:
S3,建立实时跟踪预测的目标函数:J=γTHγ+2fTγ;
其中:J为目标函数,H和f为预测输出电流的系数矩阵,γ为待求解的最优序列,H=ΩTQyΩ+Qm,f=2ζTQyζ-2dTQm,ζ为预测输出模型中状态变量的系数矩阵,Ω为预测输出模型中序列γ的系数矩阵,d为变换器的占空比初始值,Qy和Qm均为正定权矩阵,T为矩阵转置符号;
S4,将所述目标函数在所述约束条件下求解,获得用于跟踪光伏最大功率点的最优控制序列。
5.根据权利要求4所述的光伏最大功率点实时跟踪预测控制方法,其特征在于:将光伏变换器在断开和闭合状态的状态空间模型为转化为线性函数:
x(k+1)=A′×x(k)+B'×d(k)+R×δ(k)+G×z(k);
其中:k为时间状态参数;x(k)为k时刻状态变量;x(k+1)为k+1时刻状态变量;d(k)为k时刻变换器的占空比;δ(k)为k时刻理想开关的开关变量;z(k)为k时刻辅助变量;A'、B'、R、G为系数矩阵。
6.根据权利要求1所述的光伏最大功率点实时跟踪预测控制方法,其特征在于:根据环境因子与光伏控制器中光伏最大功率点函数模型、所述光伏变换器的空间状态方程、预测时域和控制时域,建立用于确定目标函数的性能指标:
其中:L为性能指标,k,j为时间状态参数;Q=QT≥0,Qy和Qm分别为正定权矩阵;i(k+j/k)为在k时刻对k+j时刻的预测输出电流;i*(j)为k+j 时刻的参考电流;Δd(k+j-1/k)为在k时刻对k+j时刻和k+j-1时刻预测占空比之差;P为预测时域,M为控制时域。
7.根据权利要求6所述的光伏最大功率点实时跟踪预测控制方法,其特征在于:在步骤S3中,建立实时跟踪预测的目标函数的过程为:
1)建立用于获取预测输出的预测输出函数模型:
y(k+j/k)=C×(A',A′2,…,A′P)T×x(k)+Ω×γ(k+j/k);
2)将获取的预测输出带入到所述性能指标,获得目标函数;
其中:γ(k+j/k)={d(i)T,δ(i)T,z(i)T},i∈{k,...,k+M-1};
y(k+j/k)为在k时刻对k+j时刻的预测输出;Ω为总序列的系数矩阵;C为输出矩阵;A'、B'为状态变量x(k+1)中的系数矩阵;d(i)T为占空比预测值序列;δ(i)T为开关变量预测值序列;z(i)T为辅助变量预测值序列;γ(k+j/k)为占空比、开关变量和辅助变量的总序列;P为预测时域,M为控制时域;R为负载电阻;G为实时光照;T为矩阵转置符号。
8.根据权利要求1所述的光伏最大功率点实时跟踪预测控制方法,其特征在于:在步骤S4中,将所述目标函数在所述约束条件下求解具体包括:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210514640.5A CN102981549B (zh) | 2012-12-05 | 2012-12-05 | 一种光伏最大功率点实时跟踪预测控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210514640.5A CN102981549B (zh) | 2012-12-05 | 2012-12-05 | 一种光伏最大功率点实时跟踪预测控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102981549A true CN102981549A (zh) | 2013-03-20 |
CN102981549B CN102981549B (zh) | 2014-07-02 |
Family
ID=47855696
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210514640.5A Expired - Fee Related CN102981549B (zh) | 2012-12-05 | 2012-12-05 | 一种光伏最大功率点实时跟踪预测控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102981549B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104020813A (zh) * | 2014-05-13 | 2014-09-03 | 安徽省安泰科技股份有限公司 | 基于fir滤波器预测的mppt滞环控制算法 |
CN104113079A (zh) * | 2014-07-03 | 2014-10-22 | 珠海格力电器股份有限公司 | Mppt控制方法和系统 |
CN105594142A (zh) * | 2013-08-13 | 2016-05-18 | 阿尔斯通技术有限公司 | 用于调节光伏变换器的电力供应的方法及设备 |
CN105974995A (zh) * | 2016-05-13 | 2016-09-28 | 上海电机学院 | 一种最大功率点追踪系统及方法 |
CN106950856A (zh) * | 2017-04-27 | 2017-07-14 | 南通大学 | 基于混合逻辑动态的mppt建模仿真方法 |
CN107147361A (zh) * | 2017-05-27 | 2017-09-08 | 重庆大学 | 局部遮挡情况下多光伏组件效率提升的最优控制方法 |
CN111624883A (zh) * | 2020-06-03 | 2020-09-04 | 广西大学 | 一种高维多步分数阶模型预测的光伏最大功率点跟踪方法 |
CN112274143A (zh) * | 2020-12-29 | 2021-01-29 | 四川写正智能科技有限公司 | 一种人体存在状态的检测方法及检测装置 |
CN112711292A (zh) * | 2021-03-29 | 2021-04-27 | 深圳黑晶光电技术有限公司 | 一种光伏组件最大功率追踪方法、系统及存储介质 |
CN114188989A (zh) * | 2021-12-14 | 2022-03-15 | 福州大学 | 一种考虑相关性的分布式光伏电站出力预测方法 |
CN114661088A (zh) * | 2022-03-08 | 2022-06-24 | 重庆邮电大学 | 不均匀光照下的光伏mppt控制方法及系统 |
CN118214078A (zh) * | 2024-05-20 | 2024-06-18 | 国家电投集团邢台新能源发电有限公司 | 一种光伏逆变器的智能控制方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101651436A (zh) * | 2009-09-16 | 2010-02-17 | 合肥阳光电源有限公司 | 一种高精度最大功率点跟踪方法 |
US20100176771A1 (en) * | 2009-01-15 | 2010-07-15 | Xantrex Technology Inc. | Perturb voltage as a decreasing non-linear function of converter power |
CN102136734A (zh) * | 2010-09-08 | 2011-07-27 | 上海岩芯电子科技有限公司 | 一种光伏微型并网逆变器最大功率点跟踪方法 |
CN102681589A (zh) * | 2012-05-31 | 2012-09-19 | 浙江大学 | 一种光伏发电中追踪最大功率点的装置和方法 |
CN102801363A (zh) * | 2011-05-24 | 2012-11-28 | 中山大学 | 一种基于自适应预测的光伏系统最大功率点跟踪控制方法 |
-
2012
- 2012-12-05 CN CN201210514640.5A patent/CN102981549B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100176771A1 (en) * | 2009-01-15 | 2010-07-15 | Xantrex Technology Inc. | Perturb voltage as a decreasing non-linear function of converter power |
CN101651436A (zh) * | 2009-09-16 | 2010-02-17 | 合肥阳光电源有限公司 | 一种高精度最大功率点跟踪方法 |
CN102136734A (zh) * | 2010-09-08 | 2011-07-27 | 上海岩芯电子科技有限公司 | 一种光伏微型并网逆变器最大功率点跟踪方法 |
CN102801363A (zh) * | 2011-05-24 | 2012-11-28 | 中山大学 | 一种基于自适应预测的光伏系统最大功率点跟踪控制方法 |
CN102681589A (zh) * | 2012-05-31 | 2012-09-19 | 浙江大学 | 一种光伏发电中追踪最大功率点的装置和方法 |
Non-Patent Citations (1)
Title |
---|
黄艺,何毅,何振雨: "基于模型预测控制的MPPT控制策略", 《可再生能源》 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10177706B2 (en) | 2013-08-13 | 2019-01-08 | Alstom Technology Ltd | Method and device for regulating the supply of a photovoltaic converter |
CN105594142A (zh) * | 2013-08-13 | 2016-05-18 | 阿尔斯通技术有限公司 | 用于调节光伏变换器的电力供应的方法及设备 |
CN105594142B (zh) * | 2013-08-13 | 2017-12-08 | 通用电气技术有限公司 | 用于调节光伏变换器的电力供应的方法及设备 |
CN104020813A (zh) * | 2014-05-13 | 2014-09-03 | 安徽省安泰科技股份有限公司 | 基于fir滤波器预测的mppt滞环控制算法 |
CN104113079A (zh) * | 2014-07-03 | 2014-10-22 | 珠海格力电器股份有限公司 | Mppt控制方法和系统 |
CN105974995A (zh) * | 2016-05-13 | 2016-09-28 | 上海电机学院 | 一种最大功率点追踪系统及方法 |
CN106950856A (zh) * | 2017-04-27 | 2017-07-14 | 南通大学 | 基于混合逻辑动态的mppt建模仿真方法 |
CN107147361B (zh) * | 2017-05-27 | 2019-04-30 | 重庆大学 | 局部遮挡情况下多光伏组件效率提升的最优控制方法 |
CN107147361A (zh) * | 2017-05-27 | 2017-09-08 | 重庆大学 | 局部遮挡情况下多光伏组件效率提升的最优控制方法 |
CN111624883A (zh) * | 2020-06-03 | 2020-09-04 | 广西大学 | 一种高维多步分数阶模型预测的光伏最大功率点跟踪方法 |
CN111624883B (zh) * | 2020-06-03 | 2022-12-06 | 广西大学 | 一种高维多步分数阶模型预测的光伏最大功率点跟踪方法 |
CN112274143A (zh) * | 2020-12-29 | 2021-01-29 | 四川写正智能科技有限公司 | 一种人体存在状态的检测方法及检测装置 |
CN112274143B (zh) * | 2020-12-29 | 2021-06-01 | 四川写正智能科技有限公司 | 一种人体存在状态的检测方法及检测装置 |
CN112711292A (zh) * | 2021-03-29 | 2021-04-27 | 深圳黑晶光电技术有限公司 | 一种光伏组件最大功率追踪方法、系统及存储介质 |
CN112711292B (zh) * | 2021-03-29 | 2021-07-09 | 深圳黑晶光电技术有限公司 | 一种光伏组件最大功率追踪方法、系统及存储介质 |
CN114188989A (zh) * | 2021-12-14 | 2022-03-15 | 福州大学 | 一种考虑相关性的分布式光伏电站出力预测方法 |
CN114188989B (zh) * | 2021-12-14 | 2024-04-12 | 福州大学 | 一种考虑相关性的分布式光伏电站出力预测方法 |
CN114661088A (zh) * | 2022-03-08 | 2022-06-24 | 重庆邮电大学 | 不均匀光照下的光伏mppt控制方法及系统 |
CN118214078A (zh) * | 2024-05-20 | 2024-06-18 | 国家电投集团邢台新能源发电有限公司 | 一种光伏逆变器的智能控制方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN102981549B (zh) | 2014-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102981549B (zh) | 一种光伏最大功率点实时跟踪预测控制方法 | |
Aashoor et al. | A variable step size perturb and observe algorithm for photovoltaic maximum power point tracking | |
Reisi et al. | Classification and comparison of maximum power point tracking techniques for photovoltaic system: A review | |
Zainuri et al. | Adaptive P&O-fuzzy control MPPT for PV boost dc-dc converter | |
Chin et al. | Fuzzy logic based MPPT for PV array under partially shaded conditions | |
Abdellatif et al. | A Fuzzy Logic Controller Based MPPT Technique for Photovoltaic Generation System. | |
CN106950856A (zh) | 基于混合逻辑动态的mppt建模仿真方法 | |
Siddique et al. | Maximum power point tracking with modified incremental conductance technique in grid-connected PV array | |
Naick et al. | Fuzzy logic controller based PV system connected in standalone and grid connected mode of operation with variation of load | |
Yang et al. | Photovoltaic array MPPT based on improved variable step size incremental conductance algorithm | |
Kumar et al. | Development of fuzzy logic based MPPT controller for PV system at varying meteorological parameters | |
CN104578091A (zh) | 一种含多源电网的无延时最优无功协调控制系统及方法 | |
Yahiaoui et al. | Experimental validation and intelligent control of a stand-alone solar energy conversion system using dSPACE platform | |
Saber et al. | A comparative study of the FLC, INC and P&O methods of the MPPT algorithm for a PV system | |
Rao et al. | A generalized mppt controlled dc-dc boost converter for pv system connected to utility grid | |
Seba et al. | A comparative analysis of series and parallel topologies of perturb and observe (P&O) and incremental conductance (InC) MPPT algorithms for photovoltaic system | |
CN105824347A (zh) | 改进的基于β参数的变步长最大功率点跟踪控制方法 | |
Kanchev et al. | Modelling and control of a grid-connected PV system for smart grid integration | |
Murtaza et al. | Optimization of the perturb and observe maximum power point tracker for a distributed photovoltaic system | |
Feng et al. | A variable step sizes perturb and observe MPPT method in PV system based on flyback converter | |
Truc et al. | Improved Sliding Mode Control Based on Perturb & Observe for Maximum Power Point Tracking from PV System | |
Zhang | Control simulation and experimental verification of maximum power point tracking based on rt-lab | |
Sagonda et al. | Comparison of three techniques for maximum power point tracking of solar PV | |
Pillai et al. | A PV panel model considering temperature and irradiance effect on panel resistances | |
Khan et al. | Real-time solar power optimization and energy monitoring system with maximum power point tracking |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140702 Termination date: 20161205 |
|
CF01 | Termination of patent right due to non-payment of annual fee |