发明名称 一种风光互补发电系统最大功率跟踪控制方法
摘要 本发明提供了一种风光互补发电系统最大功率跟踪控制方法,结合扰动观察法和二次插值法计算风光互补发电系统的最大功率,其特点是:(1)扰动观察法测量参数少,结构简单,通过不断扰动趋于系统的最大功率点;(2)二次插值法通过拟合二次曲线得出光伏发电的最大功率点电压以及风力发电的最大功率点转速;(3)本专利将两种方法进行结合可以避免扰动观察法震荡所产生的功率损失,扬长避短,可以快速精确的找到最大功率点。
申请公布号 CN102749955B 申请公布日期 2014.09.10
申请号 CN201210253156.1 申请日期 2012.07.20
申请人 北方民族大学 发明人 张秀霞;张丽霞;刘婷;娄秀丽;刘芳;高岳林;马行
分类号 G05F1/67(2006.01)I 主分类号 G05F1/67(2006.01)I
代理机构 西安通大专利代理有限责任公司 61200 代理人 徐文权
主权项 一种风光互补发电系统最大功率跟踪控制方法,其特征在于:包括以下步骤:(1)给定太阳能发电系统的步长ΔU与初始电压U<sub>0</sub>和初始电流I<sub>0</sub>及风力发电系统的步长Δω与初始转速ω<sub>0</sub>,然后测量太阳能发电系统的电压U<sub>1</sub>和电流I<sub>1</sub>或者风力发电系统的转速ω<sub>1</sub>,得到功率P<sub>1</sub>,采用扰动观察法搜索太阳能发电系统的最大功率点电压U<sub>m</sub>或风力发电系统的最大功率点转速ω<sub>m</sub>,如果该电压U<sub>m</sub>与前一周期的最大功率点电压U<sub>m‑1</sub>的差值小于δ或者该转速ω<sub>m</sub>与前一周期的最大功率点转速ω<sub>m‑1</sub>的差值小于δ,则停止搜索,δ为系统设定的最大功率点电压或者转速的临界值;(2)在步骤(1)的搜索结果中选取最大功率点附近的三个点电压U<sub>1</sub>’、U<sub>2</sub>’和U<sub>3</sub>’或者三个点转速ω<sub>1</sub>’,ω<sub>2</sub>’,ω<sub>3</sub>’,使其满足U<sub>m</sub>‑N<U<sub>1</sub>’、U<sub>2</sub>’、U<sub>3</sub>’<U<sub>m</sub>+N或者ω<sub>m</sub>‑N<ω<sub>1</sub>’、ω<sub>2</sub>’、ω<sub>3</sub>’<ω<sub>m</sub>+N,其中,N是最大功率点电压或者转速附近的一段范围,如果不满足,重新选取U<sub>1</sub>’、U<sub>2</sub>’和U<sub>3</sub>’或ω<sub>1</sub>’、ω<sub>2</sub>’和ω<sub>3</sub>’;(3)判断<img file="FDA0000517041400000011.GIF" wi="1018" he="157" />或者<img file="FDA0000517041400000012.GIF" wi="1009" he="170" />是否成立,其中,P<sub>1</sub>’、P<sub>2</sub>’、P<sub>3</sub>’分别为最大功率点附近的三个点电压U<sub>1</sub>’、U<sub>2</sub>’、U<sub>3</sub>’或者三个点转速ω<sub>1</sub>’,ω<sub>2</sub>’,ω<sub>3</sub>’所对应的功率值,如果成立,则求出二次插值法图像的系数a,b,c,如果不成立,则令U<sub>3</sub>’=U<sub>2</sub>’,U<sub>2</sub>’=U<sub>1</sub>’,P<sub>3</sub>’=P<sub>2</sub>’,P<sub>2</sub>’=P<sub>1</sub>’,或者令ω<sub>3</sub>’=ω<sub>2</sub>’,ω<sub>2</sub>’=ω<sub>1</sub>’,P<sub>3</sub>’=P<sub>2</sub>’,P<sub>2</sub>’=P<sub>1</sub>’,然后则求出二次插值法图像的系数a,b,c,其中,对于太阳能发电系统,二次插值多项式的函数一般式为:L(U)=aU<sup>2</sup>+bU+c,其中,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>a</mi><mo>=</mo><mfrac><msup><msub><mi>P</mi><mn>1</mn></msub><mo>,</mo></msup><mrow><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>1</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>U</mi><mn>2</mn></msub><mo>,</mo></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>1</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>U</mi><mn>3</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow></mfrac><mo>+</mo><mfrac><msup><msub><mi>P</mi><mn>2</mn></msub><mo>,</mo></msup><mrow><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>2</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>U</mi><mn>1</mn></msub><mo>,</mo></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>2</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>U</mi><mn>3</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow></mfrac><mo>+</mo><mfrac><msup><msub><mi>P</mi><mn>3</mn></msub><mo>,</mo></msup><mrow><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>3</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>U</mi><mn>1</mn></msub><mo>,</mo></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>3</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>U</mi><mn>2</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000517041400000021.GIF" wi="1464" he="145" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>b</mi><mo>=</mo><mfrac><mrow><msup><msub><mi>P</mi><mn>1</mn></msub><mo>,</mo></msup><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>2</mn></msub><mo>,</mo></msup><mo>+</mo><msup><msub><mi>U</mi><mn>3</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow><mrow><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>1</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>U</mi><mn>2</mn></msub><mo>,</mo></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>1</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>U</mi><mn>3</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><msub><mi>P</mi><mn>2</mn></msub><mo>,</mo></msup><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>1</mn></msub><mo>,</mo></msup><mo>+</mo><msup><msub><mi>U</mi><mn>3</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow><mrow><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>2</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>U</mi><mn>1</mn></msub><mo>,</mo></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>2</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>U</mi><mn>3</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><msub><mi>P</mi><mn>3</mn></msub><mo>,</mo></msup><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>1</mn></msub><mo>,</mo></msup><mo>+</mo><msup><msub><mi>U</mi><mn>2</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow><mrow><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>3</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>U</mi><mn>1</mn></msub><mo>,</mo></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>3</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>U</mi><mn>2</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000517041400000022.GIF" wi="1462" he="145" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>c</mi><mo>=</mo><mfrac><mrow><msup><msub><mi>P</mi><mn>1</mn></msub><mo>,</mo></msup><msup><msub><mi>U</mi><mn>2</mn></msub><mo>,</mo></msup><msup><msub><mi>U</mi><mn>3</mn></msub><mo>,</mo></msup></mrow><mrow><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>1</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>U</mi><mn>2</mn></msub><mo>,</mo></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>1</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>U</mi><mn>3</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><msub><mi>P</mi><mn>2</mn></msub><mo>,</mo></msup><msup><msub><mi>U</mi><mn>1</mn></msub><mo>,</mo></msup><msup><msub><mi>U</mi><mn>3</mn></msub><mo>,</mo></msup></mrow><mrow><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>2</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>U</mi><mn>1</mn></msub><mo>,</mo></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>2</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>U</mi><mn>3</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><msub><mi>P</mi><mn>3</mn></msub><mo>,</mo></msup><msup><msub><mi>U</mi><mn>1</mn></msub><mo>,</mo></msup><msup><msub><mi>U</mi><mn>2</mn></msub><mo>,</mo></msup></mrow><mrow><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>3</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>U</mi><mn>1</mn></msub><mo>,</mo></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><msub><mi>U</mi><mn>3</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>U</mi><mn>2</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000517041400000023.GIF" wi="1483" he="145" /></maths>对于风力发电系统,二次差值多项式的函数一般式为:L(ω)=aω<sup>2</sup>+bω+c,其中,<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>a</mi><mo>=</mo><mfrac><msup><msub><mi>P</mi><mn>1</mn></msub><mo>,</mo></msup><mrow><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>1</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>2</mn></msub><mo>,</mo></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>1</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>3</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow></mfrac><mo>+</mo><mfrac><msup><msub><mi>P</mi><mn>2</mn></msub><mo>,</mo></msup><mrow><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>2</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>1</mn></msub><mo>,</mo></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>2</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>3</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow></mfrac><mo>+</mo><mfrac><msup><msub><mi>P</mi><mn>3</mn></msub><mo>,</mo></msup><mrow><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>3</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>1</mn></msub><mo>,</mo></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>3</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>2</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000517041400000024.GIF" wi="1422" he="145" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>b</mi><mo>=</mo><mfrac><mrow><msup><msub><mi>P</mi><mn>1</mn></msub><mo>,</mo></msup><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>2</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>3</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow><mrow><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>1</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>2</mn></msub><mo>,</mo></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>1</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>3</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><msub><mi>P</mi><mn>2</mn></msub><mo>,</mo></msup><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>1</mn></msub><mo>,</mo></msup><mo>+</mo><msup><msub><mi>&omega;</mi><mn>3</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow><mrow><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>2</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>1</mn></msub><mo>,</mo></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>2</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>3</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><msub><mi>P</mi><mn>3</mn></msub><mo>,</mo></msup><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>1</mn></msub><mo>,</mo></msup><mo>+</mo><msup><msub><mi>&omega;</mi><mn>2</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow><mrow><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>3</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>1</mn></msub><mo>,</mo></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>3</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>2</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000517041400000025.GIF" wi="1420" he="145" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>c</mi><mo>=</mo><mfrac><mrow><msup><msub><mi>P</mi><mn>1</mn></msub><mo>,</mo></msup><msup><msub><mi>&omega;</mi><mn>2</mn></msub><mo>,</mo></msup><msup><msub><mi>&omega;</mi><mn>3</mn></msub><mo>,</mo></msup></mrow><mrow><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>1</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>2</mn></msub><mo>,</mo></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>1</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>3</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><msub><mi>P</mi><mn>2</mn></msub><mo>,</mo></msup><msup><msub><mi>&omega;</mi><mn>1</mn></msub><mo>,</mo></msup><msup><msub><mi>&omega;</mi><mn>3</mn></msub><mo>,</mo></msup></mrow><mrow><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>2</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>1</mn></msub><mo>,</mo></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>2</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>3</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><msub><mi>P</mi><mn>3</mn></msub><mo>,</mo></msup><msup><msub><mi>&omega;</mi><mn>1</mn></msub><mo>,</mo></msup><msup><msub><mi>&omega;</mi><mn>2</mn></msub><mo>,</mo></msup></mrow><mrow><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>3</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>1</mn></msub><mo>,</mo></msup><mo>)</mo></mrow><mrow><mo>(</mo><msup><msub><mi>&omega;</mi><mn>3</mn></msub><mo>,</mo></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>2</mn></msub><mo>,</mo></msup><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000517041400000026.GIF" wi="1418" he="145" /></maths>(4)根据公式U<sub>ref</sub>=‑2a/b或者ω<sub>ref</sub>=‑2a/b计算二次插值法得到的最大功率点电压或转速;(5)t时间后,判断比较U<sub>ref</sub>‑U<sub>m</sub><δ或ω<sub>ref</sub>‑ω<sub>m</sub><δ是否成立,如果成立,则输出风光互补发电系统的最大功率点U<sub>ref</sub>或者转速ω<sub>ref</sub>,如果不成立,重复步骤(1)至步骤(4)直至U<sub>ref</sub>‑U<sub>m</sub><δ或者ω<sub>ref</sub>‑ω<sub>m</sub><δ成立,其中,U<sub>m</sub>和ω<sub>m</sub>分别为扰动观察法得到的最大功率点电压和转速;采用扰动观察法搜索风光互补发电系统的最大功率点的方法包括以下步骤:(a)对于太阳能发电系统计算P<sub>1</sub>和P<sub>2</sub>,P<sub>1</sub>=U<sub>1</sub>×I<sub>1</sub>,P<sub>2</sub>=U<sub>2</sub>×I<sub>1</sub>,U<sub>2</sub>=(U<sub>1+</sub>ΔU),其中,P<sub>1</sub>为测量的电压U<sub>1</sub>所对应的功率值;对于风力发电系统计算P<sub>1</sub>=kω<sub>1</sub><sup>3</sup>,P<sub>2</sub>=kω<sub>2</sub><sup>3</sup>,ω<sub>2</sub>=(ω<sub>1</sub>+Δω),k=4A×ρ×C<sub>T</sub>×π<sup>3</sup>×R<sup>3</sup>,其中,A为风力发电的扇叶面积,ρ为空气密度,C<sub>T</sub>为风能利用率,R为风轮半径,其中,P<sub>1</sub>为测量的转速ω<sub>1</sub>所对应的功率值;(b)判断P<sub>2</sub>和P<sub>1</sub>的关系,如果P<sub>1</sub>>P<sub>2</sub>,且U<sub>1</sub>>U<sub>2</sub>,或者P<sub>1</sub><P<sub>2</sub>,且U<sub>1</sub><U<sub>2</sub>或者P<sub>1</sub>>P<sub>2</sub>,且ω<sub>1</sub>>ω<sub>2</sub>,或者P<sub>1</sub><P<sub>2</sub>,且ω<sub>1</sub><ω<sub>2</sub>,则最大功率点电压U<sub>m</sub>=U<sub>2</sub>+ΔU,最大功率点转速ω<sub>m</sub>=ω<sub>2</sub>+Δω,如果P<sub>1</sub>>P<sub>2</sub>,且U<sub>1</sub><U<sub>2</sub>,或者P<sub>1</sub><P<sub>2</sub>,且U<sub>1</sub>>U<sub>2</sub>,或者P<sub>1</sub>>P<sub>2</sub>,且ω<sub>1</sub><ω<sub>2</sub>,或者P<sub>1</sub><P<sub>2</sub>,且ω<sub>1</sub>>ω<sub>2</sub>,则最大功率点电压U<sub>m</sub>=U<sub>2</sub>‑ΔU,最大功率点转速ω<sub>m</sub>=ω<sub>2</sub>‑Δω。
地址 750021 宁夏回族自治区银川市西夏区文昌北路204号