CN108897368B - 一种适用于部分遮挡条件下的多峰mppt方法 - Google Patents
一种适用于部分遮挡条件下的多峰mppt方法 Download PDFInfo
- Publication number
- CN108897368B CN108897368B CN201810009311.2A CN201810009311A CN108897368B CN 108897368 B CN108897368 B CN 108897368B CN 201810009311 A CN201810009311 A CN 201810009311A CN 108897368 B CN108897368 B CN 108897368B
- Authority
- CN
- China
- Prior art keywords
- point
- power
- value
- dyn
- maximum
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05F—SYSTEMS FOR REGULATING ELECTRIC OR MAGNETIC VARIABLES
- G05F1/00—Automatic systems in which deviations of an electric quantity from one or more predetermined values are detected at the output of the system and fed back to a device within the system to restore the detected quantity to its predetermined value or values, i.e. retroactive systems
- G05F1/66—Regulating electric power
- G05F1/67—Regulating electric power to the maximum power available from a generator, e.g. from solar cell
-
- 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
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Automation & Control Theory (AREA)
- Control Of Electrical Variables (AREA)
Abstract
Description
技术领域
本发明涉及光伏发电领域,具体是一种适用于部分遮挡条件下的多峰MPPT方法。
背景技术
由于能源危机和环境污染,近年来对可再生能源的关注有所增加。其中,光伏发电具有清洁,占地面积小,安装便捷等优势,成为解决环境污染和化石资源枯竭的有效方式之一。然而,光伏阵列有输出效率低,输出特性为非线性的缺点。为了提高效率就必须控制光伏阵列工作在最大功率输出状态。因此,最大功率点跟踪(maximum power pointtracking,MPPT)技术尤为重要。传统MPPT技术包括:固定电压法,电导增量法,扰动观察(P&O)等。
在实际运用中,经常会发生部分光伏模块被云朵、高楼、树荫等遮挡的现象,此时阵列受到的光照不均匀,从而使光伏模块的P-U曲线含多个局域峰值点。而传统的MPPT算法只能追踪单一峰值点,可能无法在部分遮挡情况下追踪到全局最大功率点 (globalmaximum power point,GP)。人工智能方法(如神经网络、粒子群法等)对局部阴影情况时有一定效果,但神经网络法对于不同的光伏阵列系统需要进行长时间有针对性的训练;粒子群法可以提高结果的精度,但实现过程较为复杂,且会导致硬件设备成本的增加。
发明内容
本发明为了解决上述问题,提供一种适用于部分遮挡条件下的多峰MPPT方法。
本发明采取以下技术方案:一种适用于部分遮挡条件下的多峰 MPPT方法,包括以下步骤:
101~在开路电压[0,uoc]范围内,采集光伏模块的电压值u(k),同时采集对应的电流值i(k),并计算出功率:P(u)=u(k)×i(k), k=1~m,式中m为在整个[0,uoc]内总采样点数;
102~调用Dyn测度算法,对光伏阵列的全局最大功率点所对应的电压和功率进行追踪;
103~保持在全局最大功率点处运行,同时监测功率变化;
104~判断部分遮挡情况是否发生变化,当功率变化ΔP>ΔPset时,说明部分遮挡条件发生变化,重新采样计算Dyn测度,返回步骤一。
步骤102中,具体步骤如下,设P(u)在定义区间I=[a,b]上含有 n个极大值,这里区间I=[a,b]为开路电压区间[0,uoc]。
1)确定P(u)曲线上所有极大值点的位置,即电压ui,i=1~n;
2)在极大值点ui的两侧追踪离该点最近的等高点的位置;
3)在该极大值点ui到两侧等高点的区间内找到各自的功率最小值,分别与该点的功率值P(ui)相减,功率差值较小者为该点的Dyn测度值,即Dyn(ui);
4)如果只有一侧有等高点,则该侧的功率差值作为该点的Dyn测度值;如果两侧都没有等高点,则该点的动态测度值定义为该点的功率值;
5)返回2),求取所有极大值点的Dyn测度值;
6)将以上求出的所有的Dyn测度值进行比较,Dyn测度值最大值点就是全局最大功率点uM,其对应的功率值P(uM)为全局最大功率值,即:Dyn(uM)=max(Dyn(ui))。
步骤104中,功率变化的ΔP的计算公式为:
ΔP=Pk+1-Pk
式中,Pk和Pk+1为光伏系统保持在全局最大功率点电压运行后相邻的两点的功率。
与现有技术相比,本发明的有益效果在于:本发明针对部分遮挡条件下光伏阵列的输出特性,提出了基于Dyn测度的光伏阵列多峰 MPPT算法。根据动态测度能够反映信号极值点的结构,并能准确的搜寻到极值点的特点,将其应用到光伏阵列多峰MPPT算法中,避免了陷入局部极值点导致输出功率损失,提高了算法的精度与稳定性。
附图说明
图1为路径Dyn测度;
图2为极大值点的Dyn测度;
图3为本发明的总体流程图;
图4为图3中基于Dyn测度的光伏阵列多峰MPPT算法具体流程图;
图5为光伏组件遮挡情况;
图6为对应图5中(b)的P-U曲线图;
图7为图6中对应的动态测度;
图8为对应图5中(c)的P-U曲线图;
图9为图8中对应的动态测度;
图10为对应图5中(d)的P-U曲线图;
图11为图10中对应的动态测度;
图12为在遮挡条件变化的情况下,基于Dyn测度的多峰MPPT算法和传统MPPT算法的电压曲线比较图;
图13为在遮挡条件变化的情况下,基于Dyn测度的多峰MPPT算法和传统MPPT算法的功率追踪曲线比较图。
具体实施方式
下面结合附图和实施例对本发明作进一步的说明。
如图3、图4所示,一种光伏阵列多峰最大功率群搜索优化跟踪方法,包括以下步骤:
步骤一、在开路电压范围内,采集光伏模块的电压值,同时采集对应的电流值,并计算出功率;
步骤二、调用Dyn测度算法,对光伏阵列的全局最大功率点所对应的电压和功率进行追踪,其中具体包括:
用P-U曲线代替f(x)信号曲线,设P(u)在定义区间I=[a,b]上含有n个极大值,这里区间I=[a,b]为开路电压区间[0,uoc]。
(1)确定P(u)曲线上所有极大值点的位置(即电压ui,i=1~n);
(2)在极大值点ui的两侧追踪离该点最近的等高点的位置;
(3)在该极大值点ui到两侧等高点的区间内找到各自的功率最小值,分别与该点的功率值P(ui)相减,功率差值较小者为该点的Dyn测度值,即Dyn(ui);
(4)如果只有一侧有等高点,则该侧的功率差值作为该点的Dyn测度值;如果两侧都没有等高点,则该点的动态测度值定义为该点的功率值;
(5)返回(2),求取所有极大值点的Dyn测度值;
(6)将以上求出的所有的Dyn测度值进行比较,Dyn测度值最大值点就是全局最大功率点uM,其对应的功率值P(uM)为全局最大功率值,即:
Dyn(uM)=max(Dyn(ui));
步骤三、保持在全局最大功率点处运行,同时监测功率变化;
步骤四、判断部分遮挡情况是否发生变化,当功率变化ΔP>ΔPset时,说明部分遮挡条件发生变化,重新采样计算Dyn测度,返回步骤一。功率变化的ΔP的计算公式为:
ΔP=Pk+1-Pk
式中,Pk和Pk+1为光伏系统保持在全局最大功率点电压运行后相邻的两点的功率,本文中设置ΔPset=100W。
本发明实施例如下:
为验证上述算法,在Matlab环境下建立由4个相同的光伏组件串联而成的光伏阵列和基于Dyn测度算法的多峰MPPT控制模型。光伏组件采用英利公司生产的单晶硅太阳能电池板,其型号为 YL-LW235,在标准测试条件(温度t=25℃,光照强度S=1000W/m2) 下,光伏组件的标称参数如下:组件开路电压Uoc=37.2V,组件短路电流Isc=8.54A,最大功率点工作电压Um=29.6V,最大功率点工作电流Im=8A,最大功率Pm=235W。
仿真的遮挡情况如图5所示,图5(a)所有模块阵列的光照强度都为1000W/m2;(b)M1、M2和M3的光照强度为600W/m2,M4为1000 W/m2;(c)M1的光照强度为200W/m2,M2为600W/m2,M3和M4为1000 W/m2;(d)M1的光照强度为200W/m2,M2为600W/m2,M3为800W/m2, M4为1000W/m2。
如图6~图11所示对应的图5中(b)-(d)所示不同遮挡条件下的 P-U曲线和其对应的动态测度。从图中可以看出,Dyn测度算法能准确地追踪到所有的局部峰值点,同时发现GP的Dyn测度值比其他局部峰值点的Dyn测度值大很多,且GP的Dyn测度值还等于全局最大功率值。通过比较追踪到这个最大的动态测度点,可以直接追踪到 GP电压,同时计算出全局最大功率值。
如图12和图13所示,为本发明在遮挡条件变化的情况下,将基于Dyn测度的多峰MPPT算法和传统MPPT算法的电压和功率追踪曲线作比较,图中遮挡情况按图5中(a)→(b)→(c)→(d)的顺序变化,每种情况持续1s。
0-1s内,光伏阵列如图5(a)所示没有被遮挡,本发明算法和传统MPPT算法分别在0.11s和0.48s追踪到实际的最大功率点;实际 GP电压为123.2V,本发明算法追踪到的GP电压123V,相差0.2V。实际的GP功率值为959.71W,本发明算法追踪到的GP功率值为 958.73W,误差为0.113%。
1-2s内光伏阵列的遮挡方式如图5(b)所示,其P-U曲线如图6 所示。GP在第一个局部峰值点处,所以这段时间内,两种算法依然可以追踪到最大功率点,分别用时0.092s和0.24s;实际GP电压为 91.9V,本发明算法追踪到的GP电压92V,相差0.1V。实际的GP功率值为715.72W,本发明算法追踪到的GP功率值为714.91W,误差为 0.112%。
在3-4s和4-5s时间段内,此时的P-U曲线如图8和图10所示, GP不在第一个局部峰值点的位置,传统的MPPT算法只能追踪到第一个局部极值点,从而导致输出功率大大减小,追踪到的功率分别为 237.91W和237.82W,而在这段时间内,与实际的功率值454.39W和320.73W相差217.02W和82.91W,功率损失高达47.7%和25.8%;本发明算法能够在2.13s和3.091s时准确的找到GP,追踪到的GP功率值为454.39W和320.36W,误差为0.105%和0.115%。光伏模型的平均输出效率可达到99.88%。
综上所述,本发明算法不但能更短时间内追踪到GP,而且在多峰情况下以更小的误差追踪到GP,不会陷入局部峰值点,不会出现传统MPPT算法在GP点功率震荡的现象。
动态(Dynamics,Dyn)测度是基于地质学提出的一种极值点评价测度,该测度能反映信号极值点的结构,并能准确的搜寻到极值点。其具体特性为:
设f(x)为任意随x变化的信号,且f(x)≥0,可以将各点的采样值看作是高度。
定义1两点之间的路径:
设m,n为f(x)上两个不同的点,则f(x)上这两点间的曲线称为路径P(m,n)即
P(m,n)=(p1,p2…,pN) (1)
定义2路径的Dyn测度:
路径P(m,n)的Dyn测度定义为P(m,n)上最大值点和最小值点的高度差,如图1所示,即
Dyn[P(m,n)]={sup(|halt(pi)-halt(pj)|);pi,pj∈P(m,n)} (2)
式中:sup表示上确界;halt表示高度。
定义3极大值点的Dyn测度:
设A为f(x)的一个极大值点,如果存在比A点更高的极大值点时,则极大值点A的Dyn测度等于由点A通向同高度点的所有路径 Dyn测度中的最小值。
Dyn(A)={inf{Dyn[P(A,B)]};halt(A)=halt(B)} (3)
式中:inf表示下确界,halt表示高度。
如图2所示,极大值点A两侧各有一个或多个比点A更高的极大值点时,则点A两侧会存在两个点B1和B2与点A等高。由点A到点B1的路径记为L1,由点A到点B2的路径记为L2。极大值点A 的Dyn测度等于路径L1的Dyn测度和路径L2的Dyn测度中的最小值。图2中路径L2的Dyn测度小于路径L1的Dyn测度,所以极大值点A 的Dyn测度应该等于路径L2的Dyn测度。即
Dyn(A)=min(Dyn(L1),Dyn(L2))=Dyn(L2) (4)
并且极大值点A的Dyn测度与路径L1或路径L2的长度无关。如果极大值点A只有一侧有比点A更高的点时,这时只存在路径L1 或路径L2,则极大值点A的Dyn测度等于路径L1或路径L2的Dyn测度。若点A为信号的最大值点时,设置其Dyn测度等于信号最大值点和最小值点的高度差。因此,最大值点的Dyn测度比其它极大值点的 Dyn测度都大(这一特点有助于我们在多个局部峰值点之中找出全局最大峰值点)。
Claims (2)
1.一种适用于部分遮挡条件下的光伏阵列多峰最大功率点追踪方法,其特征在于:包括以下步骤:
101~在开路电压[0,uoc]范围内,采集光伏模块的电压值u(k),同时采集对应的电流值i(k),并计算出功率:P(u)=u(k)×i(k),k=1~m,式中m为在整个[0,uoc]内总采样点数;
102~调用Dyn测度算法,对光伏阵列的全局最大功率点所对应的电压和功率进行追踪;具体步骤如下,设P(u)在定义区间I=[a,b]上含有n个极大值,这里区间I=[a,b]为开路电压区间[0,uoc];
1)确定P(u)曲线上所有极大值点的位置,即电压ui,i=1~n;
2)在极大值点ui的两侧追踪离该点最近的等高点的位置;
3)在该极大值点ui到两侧等高点的区间内找到各自的功率最小值,分别与该点的功率值P(ui)相减,功率差值较小者为该点的Dyn测度值,即Dyn(ui);
4)如果只有一侧有等高点,则该侧的功率差值作为该点的Dyn测度值;如果两侧都没有等高点,则该点的动态测度值定义为该点的功率值;
5)返回2),求取所有极大值点的Dyn测度值;
6)将以上求出的所有的Dyn测度值进行比较,Dyn测度值最大值点就是全局最大功率点uM,其对应的功率值P(uM)为全局最大功率值,即:Dyn(uM)=max(Dyn(ui));
103~保持在全局最大功率点处运行,同时监测功率变化;
104~判断部分遮挡情况是否发生变化,当功率变化ΔP>ΔPset时,说明部分遮挡条件发生变化,重新采样计算Dyn测度,返回步骤101。
2.根据权利要求1所述的一种适用于部分遮挡条件下的光伏阵列多峰最大功率点追踪方法,其特征在于:步骤104中,功率变化的ΔP的计算公式为:
ΔP=Pk+1-Pk
式中,Pk和Pk+1为光伏系统保持在全局最大功率点电压运行后相邻的两点的功率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810009311.2A CN108897368B (zh) | 2018-01-04 | 2018-01-04 | 一种适用于部分遮挡条件下的多峰mppt方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810009311.2A CN108897368B (zh) | 2018-01-04 | 2018-01-04 | 一种适用于部分遮挡条件下的多峰mppt方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108897368A CN108897368A (zh) | 2018-11-27 |
CN108897368B true CN108897368B (zh) | 2020-01-03 |
Family
ID=64342398
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810009311.2A Active CN108897368B (zh) | 2018-01-04 | 2018-01-04 | 一种适用于部分遮挡条件下的多峰mppt方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108897368B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113359940B (zh) * | 2021-07-02 | 2022-05-17 | 吉林建筑大学 | 一种光伏阵列的最大功率点追踪方法及系统 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101873090B (zh) * | 2010-07-06 | 2012-07-18 | 太原科技大学 | 有部分遮蔽的光伏发电系统最大功率输出优化控制方法 |
CA2902428C (en) * | 2014-10-31 | 2024-01-09 | Majid Pahlevaninezhad | Current sensorless mppt for pv micro-inverters |
CN105930918B (zh) * | 2016-04-11 | 2019-07-02 | 北京交通大学 | 应用于多峰mppt的整体分布-粒子群优化算法 |
CN106527570B (zh) * | 2016-12-20 | 2018-06-15 | 湘潭大学 | 一种光伏阵列多峰最大功率群搜索优化跟踪方法 |
-
2018
- 2018-01-04 CN CN201810009311.2A patent/CN108897368B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108897368A (zh) | 2018-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11190131B2 (en) | Method and system for power prediction of photovoltaic power station based on operating data of grid-connected inverters | |
Mutoh et al. | Prediction-data-based maximum-power-point-tracking method for photovoltaic power generation systems | |
CN105867514B (zh) | 一种光伏系统多峰值最大功率跟踪方法及系统 | |
KR101065862B1 (ko) | 태양전지 어레이의 부분 음영 판단에 따른 태양광 발전 시스템의 최대전력 추정방법 | |
CN104793691B (zh) | 一种基于蚁群算法的局部阴影下光伏阵列全局mppt方法 | |
CN107168447B (zh) | 一种基于改进电导增量法的光伏逆变器多峰值mppt方法 | |
CN112653393B (zh) | 一种光伏系统iv诊断的控制方法及装置 | |
CN110286708B (zh) | 一种光伏阵列的最大功率跟踪控制方法及系统 | |
CN108710402B (zh) | 一种太阳能光伏组串的多峰状态最大功率点跟踪方法 | |
CN108681363A (zh) | 恒压结合变步长电导增量的mppt控制方法 | |
CN110009098A (zh) | 一种光伏电池工作温度与发电功率联合估计方法 | |
KR20170058218A (ko) | 태양광 어레이의 부분 음영에 대비한 하이브리드 최대전력추종 제어방법 | |
CN108897368B (zh) | 一种适用于部分遮挡条件下的多峰mppt方法 | |
CN111984906A (zh) | 考虑光伏与负荷时序相关性的概率潮流快速计算方法 | |
CN114389312B (zh) | 一种配电网分布式状态估计方法 | |
CN110703847A (zh) | 改进型粒子群-扰动观察法的光伏全局最大功率点追踪方法 | |
CN105259971A (zh) | 一种优化的mppt算法 | |
CN109085877B (zh) | 局部阴影条件下基于qapso算法光伏系统mppt控制方法 | |
Su et al. | Parameter extraction of photovoltaic single-diode model using integrated current–voltage error criterion | |
CN105068591B (zh) | 一种光伏阵列局部遮挡下最大功率点跟踪方法 | |
CN105159388B (zh) | 一种用于光伏微网系统中的最大功率点跟踪的方法 | |
CN117572236A (zh) | 一种基于迁移学习的锂电池荷电状态估计方法 | |
Devireddy et al. | A GMPPT algorithm for PV systems using current reference and PI curve under partial shading conditions | |
CN116126086A (zh) | 最大功率点跟踪控制方法、系统及光伏发电输出装置 | |
Rodil et al. | A novel methodology for determination of soiling on PV panels by means of grey box modelling |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |