CN102749955A - 一种风光互补发电系统最大功率跟踪控制方法 - Google Patents
一种风光互补发电系统最大功率跟踪控制方法 Download PDFInfo
- Publication number
- CN102749955A CN102749955A CN2012102531561A CN201210253156A CN102749955A CN 102749955 A CN102749955 A CN 102749955A CN 2012102531561 A CN2012102531561 A CN 2012102531561A CN 201210253156 A CN201210253156 A CN 201210253156A CN 102749955 A CN102749955 A CN 102749955A
- Authority
- CN
- China
- Prior art keywords
- omega
- maximum power
- wind
- power point
- perhaps
- 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
Landscapes
- Control Of Eletrric Generators (AREA)
- Wind Motors (AREA)
Abstract
本发明提供了一种风光互补发电系统最大功率跟踪控制方法,结合扰动观察法和二次插值法计算风光互补发电系统的最大功率,其特点是:(1)扰动观察法测量参数少,结构简单,通过不断扰动趋于系统的最大功率点;(2)二次插值法通过拟合二次曲线得出光伏发电的最大功率点电压以及风力发电的最大功率点转速;(3)本专利将两种方法进行结合可以避免扰动观察法震荡所产生的功率损失,扬长避短,可以快速精确的找到最大功率点。
Description
技术领域
本发明涉及风光互补发电系统,特别涉及一种风光互补发电系统最大功率跟踪控制方法。
背景技术
风光互补发电系统受到光照、温度、风速及方向的影响很大,在某一外界环境条件下对应着相应的最大功率点。因此,需进行最大功率跟踪控制来提高系统的发电效率。目前较为常用的最大功率点跟踪算法之一是扰动观察法,该方法与叶尖速比、增量电导和功率反馈等方法相比,测量参数少,结构简单。其算法原理是通过扰动前后功率的变化量进行扰动方向的判断,先赋予一个输出电压信号,再进行功率测量,若功率变化量大于零,则表示扰动方向正确,可向同一方向继续扰动;若功率变化量小于零,则向反方向扰动。该算法的缺点是跟踪精度、速度与跟踪步长以及初始值的选取密切相关,若跟踪步长太小,当外界环境变化较快时其跟踪速度缓慢;若跟踪步长太大,当到达最大功率点附近时,会加大震荡从而损失较多功率,有时在运行中发生程序控制失序,出现“误判”现象。
鉴于以上缺陷,实有必要提供一种改进的风光互补发电系统最大功率跟踪控制方法。
发明内容
本发明所要解决的技术问题是提供一种风光互补发电系统最大功率跟踪控制方法,避免扰动观察法震荡所产生的功率损失,更精确的跟踪到系统的最大功率点。
为解决以上技术问题,本发明采用以下技术方案:
一种风光互补发电系统最大功率跟踪控制方法,包括以下步骤:
(1)给定太阳能发电系统的步长△U与初始电压U0和初始电流I0及风力发电系统的步长△ω与初始转速ω,然后测量风光互补发电系统的电压U1和电流I1或者风力发电系统的转速ω1,采用扰动观察法搜索风光互补发电系统的最大功率点电压或风力发电系统的转速,如果该电压与前一周期的最大功率点电压的差值小于δ或者该转速与前一周期的最大转速的差值小于δ,则停止搜索;
(2)在步骤(1)的搜索结果中选取最大功率点附近的三个点电压U1’、U2’和U3’或者三个点转速ω1’,ω2’,ω3’,然后判断
是否成立,如果成立,则求出
二次插值法图像的系数a,b,c,如果不成立,则令U3’=U2’,U2’=U1’,P3’=P2’,P2’=P1’,或者令ω3’=ω2’,ω2’=ω1’,P3’=P2’,P2’=P1’,然后则求出二次插值法图像的系数a,b,c,其中,
对于太阳能发电系统,二次插值多项式的函数一般式为:
L(U)=aU2+bU+c,其中,
对于风力发电系统,二次差值多项式的函数一般式为:
L(ω)=aω2+bω+c,其中,
(3)根据公式Uref=-2a/b或者ωref=-2a/b计算二次插值法得到的最大功率点电压或转速;
(4)t时间后,判断比较Uref-Um<δ或ωref-ωm<δ是否成立,如果成立,则输出风光互补发电系统的最大功率点Uref或者转速ωref),如果不成立,重复步骤(1)至步骤(4)直至Uref-Um<δ或者ωref-ωm<δ成立,其中,Um和ωm分别为扰动观察法得到的最大功率点电压和转速。
所述采用扰动观察法搜索风光互补发电系统的最大功率电压的方法包括以下步骤:
(a)对于太阳能发电系统计算P1和P2,P1=U1×I1,P2=U2×I1,U2=(U1+△U),对于风力发电系统计算P1=kω1 3,P2=kω2 3,ω2=(ω1+△ω),k=4A×ρ×CT×π3×R3,其中,A为风力发电的扇叶面积,ρ为空气密度,CT为风能利用率,R为风轮半径;
(b)判断P2和P1的关系,如果P1>P2,且U1>U2,或者P1<P2,且U1<U2或者P1>P2,且ω1>ω2,或者P1<P2,且ω1<ω2,则最大功率点电压Um=U2+△U,最大功率点转速ωm=ω2+△ω,如果P1>P2,且U1<U2,或者P1<P2,且U1>U2,或者P1>P2,且ω1<ω2,或者P1<P2,且ω1>ω2,则最大功率点电压Um=U2-△U,最大功率点转速ωm=ω2-△ω。
本发明风光互补发电系统最大功率跟踪控制方法至少具有以下优点:
(1)二次插值法的多项式图形与系统功率曲线基本吻合,可以精确的进行特定范围内的最大功率点跟踪;
(2)两者结合的方法减小了扰动观察法在最大功率点附近震荡的现象,减小功率损失,能够快速精确的跟踪到最大功率点;
(3)改进的方法可以节省搜索所需时间,避免重复不必要的工作。
附图说明
图1为本发明风光互补发电系统太阳能发电系统最大功率跟踪控制方法的流程图。
图2为本发明风光互补发电系统风力发电系统最大功率跟踪控制方法的流程图;
图3为本发明太阳能发电系统的功率曲线图;
图4为本发明风力发电系统的功率曲线图。
具体实施方式
下面结合附图对本发明风光互补发电系统最大功率跟踪控制方法做详细描述:
本发明风光互补发电系统最大功率跟踪控制方法:首先采用传统的扰动功率法进行光伏发电的最大功率点电压以及风力发电的最大功率点转速搜索;接着,采用二次插值法进行最大功率点精确搜索,得出最大功率点的电压(转速);最后,在工作一个周期后,与前一时刻进行判断,是否进行下一次的搜索任务,以确保系统运行在最大功率点附近。
具体地说,本发明方法包括以下步骤:
(1)确定合适的步长△U(对于风力发电系统为△ω)与初始电压U0,初始电流I0或者初始转速ω0,然后测量风光互补发电系统的电压U1和电流I1(对于风力发电系统,测量转速ω1),得到功率P1,接着,计算U2和P2,U2=(U1+△U),P2=U2×I1,对于风力发电系统,计算ω2和P2,其中,ω2=(ω1+△ω),P2=kω2,k=4A×ρ×CT×π3×R3,其中,A为风力发电的扇叶面积,ρ为空气密度,CT为风能利用率,R为风轮半径;
(2)判断P2和P1的关系,如果P1>P2,且U1>U2,或者P1<P2,且U1<U2,则最大功率点电压Um=U2+△U,如果P1>P2,且U1<U2,或者P1<P2,且U1>U2,则最大功率点电压Um=U2-△U(对于风力发电系统,判断如果P1>P2,且ω1>ω2,或者P1<P2,且ω1<ω2,则最大功率点转速ωm=ω2+△ω,如果P1>P2,且ω1<ω2,或者P1<P2,且ω1>ω2,则最大功率点转速ωm=ω2-△ω);
(3)判断比较本次循环得到的最大功率点电压Um或转速ωm与上次循环得到的最大功率电压Um-1或ωm-1,如果Um-Um-1<δ或者ωm-ωm-1<δ,其中,δ为系统设定的最大功率点电压或者转速的临界值,则在扰动观察法循环测量的电压或转速中选取最大功率点附近的三个电压U1’、U2’和U3’或三个转速ω1’、ω2’和ω3’,如果Um-Um-1>δ或者ωm-ωm-1>δ,则重复步骤(1)至步骤(3),直至Um-Um-1<δ或ωm-ωm-1<δ;
(4)判断U1’、U2’和U3’或者ω1’、ω2’和ω3’,使其满足Um-N<U1’、U2’、U3’<Um-N或者ωm-N<ω1’、ω2’、ω3’<ωm-N,其中,N是最大功率点电压或者转速附近的一段范围,如果不满足,重新选取U1’、U2’和U3’或ω1’、ω2’和ω3’;
(5)判断
插值法图像的系数a,b,c,如果不成立,则令U3’=U2’,U2’=U1’,P3’=P2’,P2’=P1’(对于风力发电系统,令ω3’=ω2’,ω2’=ω1’,P3’=P2’,P2’=P1’),然后则求出二次插值法图像的系数a,b,c;
对于太阳能发电系统:
二次插值多项式的函数一般式为L(U)=aU2+bU+c,其中:
对于风力发电系统:
二次插值多项式的函数一般式为L(ω)=aω2+bω+c,其中:
(6)根据公式Uref=-2a/b或ωref=-2a/b计算二次插值法得到的最大功率点电压或转速;
(7)t时间后,判断比较Uref-Um<δ或ωref-ωm<δ是否成立,如果成立,则输出风光互补发电系统的最大功率点电压Uref或转速ωref,如果不成立,重复步骤(1)至步骤(7)直至Uref-Um<δ或ωref-ωm<δ成立。
以上所述仅为本发明的一种实施方式,不是全部或唯一的实施方式,本领域普通技术人员通过阅读本发明说明书而对本发明技术方案采取的任何等效的变换,均为本发明的权利要求所涵盖。
Claims (2)
1.一种风光互补发电系统最大功率跟踪控制方法,其特征在于:包括以下步骤:
(1)给定太阳能发电系统的步长△U与初始电压U0和初始电流I0及风力发电系统的步长△ω与初始转速ω,然后测量风光互补发电系统的电压U1和电流I1或者风力发电系统的转速ω1,采用扰动观察法搜索风光互补发电系统的最大功率点电压或风力发电系统的最大功率点转速,如果该电压与前一周期的最大功率点电压的差值小于δ或者该转速与前一周期的最大功率点转速的差值小于δ,则停止搜索;
(2)在步骤(1)的搜索结果中选取最大功率点附近的三个点电压U1’、U2’和U3’或者三个点转速ω1’,ω2’,ω3’,然后判断
是否成立,其中,P1’、P2’、P3’分别为最大功率点附近的三个点电压U1’、U2’、U3’或者三个点转速ω1’,ω2’,ω3’所对应的功率值,如果成立,则求出二次插值法图像的系数a,b,c,如果不成立,则令U3’=U2’,U2’=U1’,P3’=P2’,P2’=P1’,或者令ω3’=ω2’,ω2’=ω1’,P3’=P2’,P2’=P1’,然后则求出二次插值法图像的系数a,b,c,其中,
对于太阳能发电系统,二次插值多项式的函数一般式为:
L(U)=aU2+bU+c,其中,
对于风力发电系统,二次差值多项式的函数一般式为:
L(ω)=aω2+bω+c,其中,
(3)根据公式Uref=-2a/b或者ωref=-2a/b计算二次插值法得到的最大功率点电压或转速;
(4)t时间后,判断比较Uref-Um<δ或ωref-ωm<δ是否成立,如果成立,则输出风光互补发电系统的最大功率点Uref或者转速ωref,如果不成立,重复步骤(1)至步骤(4)直至Uref-Um<δ或者ωref-ωm<δ成立,其中,Um和ωm分别为扰动观察法得到的最大功率点电压和转速。
2.如权利要求1所述的风光互补发电系统最大功率跟踪控制方法,其特征在于:所述采用扰动观察法搜索风光互补发电系统的最大功率电压的方法包括以下步骤:
(a)对于太阳能发电系统计算P1和P2,P1=U1×I1,P2=U2×I1,U2=(U1+△U),其中,P1为测量的电压U1所对应的功率值;对于风力发电系统计算P1=kω1 3,P2=kω2 3,ω2=(ω1+△ω),k=4A×ρ×CT×π3×R3,其中,A为风力发电的扇叶面积,ρ为空气密度,CT为风能利用率,R为风轮半径,其中,P1为测量的转速ω1所对应的功率值;
(b)判断P2和P1的关系,如果P1>P2,且U1>U2,或者P1<P2,且U1<U2或者P1>P2,且ω1>ω2,或者P1<P2,且ω1<ω2,则最大功率点电压Um=U2+△U,最大功率点转速ωm=ω2+△ω,如果P1>P2,且U1<U2,或者P1<P2,且U1>U2,或者P1>P2,且ω1<ω2,或者P1<P2,且ω1>ω2,则最大功率点电压Um=U2-△U,最大功率点转速ωm=ωω2-△ω。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210253156.1A CN102749955B (zh) | 2012-07-20 | 2012-07-20 | 一种风光互补发电系统最大功率跟踪控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210253156.1A CN102749955B (zh) | 2012-07-20 | 2012-07-20 | 一种风光互补发电系统最大功率跟踪控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102749955A true CN102749955A (zh) | 2012-10-24 |
CN102749955B CN102749955B (zh) | 2014-09-10 |
Family
ID=47030234
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210253156.1A Expired - Fee Related CN102749955B (zh) | 2012-07-20 | 2012-07-20 | 一种风光互补发电系统最大功率跟踪控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102749955B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105048926A (zh) * | 2015-06-30 | 2015-11-11 | 西安理工大学 | 光伏水泵系统的最大功率跟踪与U-f曲线自寻优结合的方法 |
CN105549676A (zh) * | 2016-02-16 | 2016-05-04 | 上海电机学院 | 基于回测和曲线拟合的光伏最大功率追踪方法 |
CN106712660A (zh) * | 2016-12-28 | 2017-05-24 | 北方民族大学 | 风光互补发电系统的远程监控系统及其监控方法 |
CN106762453A (zh) * | 2016-12-07 | 2017-05-31 | 湖北民族学院 | 具有发电量预测及跟踪控制的风力发电智能网及控制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20100031026A (ko) * | 2008-09-11 | 2010-03-19 | 동아전기부품 주식회사 | 차량용 태양광 발전 시스템 및 차량용 태양광 발전 시스템의 전력제어 방법 |
CN201499003U (zh) * | 2009-09-08 | 2010-06-02 | 国网电力科学研究院武汉南瑞有限责任公司 | 基于z源的风光互补发电装置 |
US20100132757A1 (en) * | 2008-12-01 | 2010-06-03 | Chung Yuan Christian University | Solar energy system |
-
2012
- 2012-07-20 CN CN201210253156.1A patent/CN102749955B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20100031026A (ko) * | 2008-09-11 | 2010-03-19 | 동아전기부품 주식회사 | 차량용 태양광 발전 시스템 및 차량용 태양광 발전 시스템의 전력제어 방법 |
US20100132757A1 (en) * | 2008-12-01 | 2010-06-03 | Chung Yuan Christian University | Solar energy system |
CN201499003U (zh) * | 2009-09-08 | 2010-06-02 | 国网电力科学研究院武汉南瑞有限责任公司 | 基于z源的风光互补发电装置 |
Non-Patent Citations (2)
Title |
---|
QIU PEI-CHUN等: "An Improved MPPT Method for Photovoltaic Power Generation System", 《2010 2ND IEEE INTERNATIONAL SYMPOSIUM ON POWER ELECTRONICS FOR DISTRIBUTED GENERATION SYSTEMS》, 31 December 2010 (2010-12-31), pages 144 - 146 * |
沈利生: "小型风光互补发电系统的建模与仿真研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》, no. 4, 15 April 2012 (2012-04-15) * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105048926A (zh) * | 2015-06-30 | 2015-11-11 | 西安理工大学 | 光伏水泵系统的最大功率跟踪与U-f曲线自寻优结合的方法 |
CN105048926B (zh) * | 2015-06-30 | 2017-11-03 | 西安理工大学 | 光伏水泵系统的最大功率跟踪与U‑f曲线自寻优结合的方法 |
CN105549676A (zh) * | 2016-02-16 | 2016-05-04 | 上海电机学院 | 基于回测和曲线拟合的光伏最大功率追踪方法 |
CN106762453A (zh) * | 2016-12-07 | 2017-05-31 | 湖北民族学院 | 具有发电量预测及跟踪控制的风力发电智能网及控制方法 |
CN106712660A (zh) * | 2016-12-28 | 2017-05-24 | 北方民族大学 | 风光互补发电系统的远程监控系统及其监控方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102749955B (zh) | 2014-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102072083B (zh) | 一种双馈风力发电系统最大风能跟踪方法 | |
CN102809980B (zh) | 基于高效自适应扰动观察法的最大功率跟踪方法 | |
Mala et al. | A twin unidirectional impulse turbine topology for OWC based wave energy plants–Experimental validation and scaling | |
CN101776043B (zh) | 基于误差补偿模型的风电机组最大风能捕获控制方法 | |
CN103019294B (zh) | 一种自适应扰动频率和步长的最大功率点跟踪方法 | |
CN104079226B (zh) | 一种同步坐标系下无锁相环的dfig控制方法 | |
CN103244348A (zh) | 变速变桨风力发电机组功率曲线优化方法 | |
CN103049034A (zh) | 一种基于功率预测的变步长扰动最大功率点跟踪方法 | |
CN102749955B (zh) | 一种风光互补发电系统最大功率跟踪控制方法 | |
CN104378033B (zh) | 一种开关磁阻电机风力发电系统最大功率点跟踪方法 | |
CN108681363A (zh) | 恒压结合变步长电导增量的mppt控制方法 | |
CN105138845B (zh) | 获得风力发电机风速值的方法 | |
CN104899465A (zh) | 风光优化配比的计算方法 | |
Wu et al. | Development of an equivalent wind farm model for frequency regulation | |
CN103557117A (zh) | 风力发电机组功率曲线获取装置 | |
CN109412163A (zh) | 一种分散式风电接入配电网的精确潮流计算方法 | |
CN203892125U (zh) | 一种基于风电机组机舱风速确定风轮前沿风速的装置 | |
CN105116959A (zh) | 基于最优梯度变步长的最大功率点跟踪控制与预测方法 | |
Bibave et al. | A novel maximum power point tracking method for wind energy conversion system: A review | |
CN103375332B (zh) | 变速变桨风力发电机组最优阻力矩动态优化方法 | |
CN103899480A (zh) | 基于Boost类变换器的风力发电系统最大功率点跟踪方法 | |
Tianpei et al. | MPPT method of wind power based on improved particle swarm optimization | |
CN204536976U (zh) | 一种实现新型最大功率跟踪算法的光伏控制装置 | |
CN104466956B (zh) | 一种分布式电源的潮流计算方法及装置 | |
Babu et al. | Optimizing electrical generators of wind energy conversion system for efficient power extraction |
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: 20140910 Termination date: 20210720 |
|
CF01 | Termination of patent right due to non-payment of annual fee |