发明名称 用混合储能平抑风电功率波动的模型预测控制方法
摘要 本发明涉及用混合储能平抑风电功率波动的模型预测控制方法,属于新能源发电中的控制技术领域。该方法包括:首先确定控制对象的拓扑结构,即包含风电、电池储能和超级电容储能的系统结构。其次,利用平均值法预测风电在控制预测区间内的出力值。然后,确定优化控制的目标函数和三个约束条件。在每一时刻求解上述优化问题,仅取得到的当前时刻的优化控制量分别输入电池和超级电容进行控制。最后,将电池和超级电容输出的功率与风电的实际输出功率加和,即可得到功率波动平抑后整个系统的并网功率。本发明提出的控制方法,可以更加经济、可靠的实现风电的功率波动平抑,而且计算简单,易于工程实现。
申请公布号 CN104104107B 申请公布日期 2016.04.06
申请号 CN201410267717.2 申请日期 2014.06.16
申请人 清华大学 发明人 陈来军;陈天一;梅生伟;郑天文
分类号 H02J3/38(2006.01)I;H02J3/24(2006.01)I;H02J3/28(2006.01)I;H02J3/32(2006.01)I 主分类号 H02J3/38(2006.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 廖元秋
主权项 用混合储能平抑风电功率波动的模型预测控制方法,其特征在于,包括如下步骤:1)确定待预测控制的分布式发电系统的拓扑结构,该分布式发电系统包括风电、电池储能、超级电容储能三个子系统,三个子系统均通过开关和变压器与电网相连;2)对风电子系统的风电功率进行超短期预测:设获取风电出力数据的时间间隔为Δt,设预测控制时间长度为M,利用当前预测时刻之前获取预测控制时间长度内的M/Δt个风电出力历史数据的平均值对当前预测时刻之后的M时间长度内的风电功率进行预测,如式(1)所示:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>P</mi><mi>W</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>10</mn></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>10</mn></munderover><msub><mi>P</mi><mi>W</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>i</mi><mo>-</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000840905190000011.GIF" wi="724" he="127" /></maths>式(1)中i=0,1,…,M‑1,k时刻是当前时刻,式子左边的P<sub>w</sub>(k+i)为风电出力的预测值,式子右边的P<sub>w</sub>(k+i‑j)为获取的风电出力历史数据;3)构建优化控制的目标函数:为了综合考虑储能平抑功率波动的经济性和输出功率的平抑效果,定义优化控制目标为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>min</mi><mi>J</mi><mo>=</mo><msub><mi>f</mi><mi>B</mi></msub><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>M</mi><mo>-</mo><mn>1</mn></mrow></munderover><msubsup><mi>u</mi><mi>B</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>f</mi><mrow><mi>u</mi><mi>c</mi></mrow></msub><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>M</mi><mo>-</mo><mn>1</mn></mrow></munderover><msubsup><mi>u</mi><mrow><mi>u</mi><mi>c</mi></mrow><mn>2</mn></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>f</mi><mi>R</mi></msub><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>M</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mrow><mo>&lsqb;</mo><msub><mi>P</mi><mi>o</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><msub><mi>P</mi><mi>o</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000840905190000012.GIF" wi="1407" he="142" /></maths>式(2)中J为优化目标;f<sub>B</sub>和f<sub>uc</sub>分别是电池和超级电容的经济性惩罚值;u<sub>B</sub>(k)和u<sub>uc</sub>(k)分别为k时刻电池和超级电容的出力大小,将这两个值作为优化变量;f<sub>R</sub>是对输出功率波动的惩罚值;P<sub>o</sub>(k)是分布式发电系统在k时刻的风电、电池储能和超级电容储能三个子系统总的输出,如式(3)所示:P<sub>o</sub>(k)=P<sub>W</sub>(k)+u<sub>B</sub>(k)+u<sub>uc</sub>(k)    (3)表达式<img file="FDA0000840905190000013.GIF" wi="407" he="125" />为该系统在预测控制区间内输出功率波动的大小;4)确定所述目标函数的约束条件:该分布式发电系统运行共有三个约束条件:4‑1)功率波动约束条件:即风电并网标准中对预测控制时间长度内输出功率波动的限制,该功率波动约束条件如式(4)所示:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mrow><mi>j</mi><mo>=</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>,</mo><mn>...</mn><mo>,</mo><mn>9</mn></mrow></munder><msub><mi>P</mi><mi>o</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>i</mi><mo>-</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mrow><mi>j</mi><mo>=</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>,</mo><mn>...</mn><mo>,</mo><mn>9</mn></mrow></munder><msub><mi>P</mi><mi>o</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>i</mi><mo>-</mo><mi>j</mi><mo>)</mo></mrow><mo>&le;</mo><mi>&gamma;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000840905190000014.GIF" wi="957" he="86" /></maths>式(4)中i=0,1,…,M‑1,γ为设定的预测控制时间长度内功率波动的限制值;4‑2)功率约束条件:即电池和超级电容的出力不能超过它们充放电的额定功率,如式(5)所示:‑P<sub>Bc</sub>≤u<sub>B</sub>(k)≤P<sub>Bd</sub>‑P<sub>ucc</sub>≤u<sub>uc</sub>(k)≤P<sub>ucd</sub>    (5)式(5)中P<sub>Bc</sub>、P<sub>Bd</sub>分别为储能的额定充放电功率;P<sub>ucc</sub>、P<sub>ucd</sub>分别为超级电容的额定充放电功率;4‑3)容量约束条件:定义k时刻电池的荷电状态(SOC)为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>S</mi><mi>O</mi><mi>C</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>Q</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mi>Q</mi></mfrac><mo>&times;</mo><mn>100</mn><mi>%</mi></mrow>]]></math><img file="FDA0000840905190000021.GIF" wi="422" he="125" /></maths>式中Q(t)是储能在t时刻的剩余电量,Q为储能可储存的总电量;则储能在预测控制区间内需不超过各自的容量限制,即容量约束条件如式(6)所示:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>S</mi><mi>B</mi></msub><mo>&lsqb;</mo><msub><mi>SOC</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>SOC</mi><mrow><mi>B</mi><mo>,</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>&rsqb;</mo><mo>&le;</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mi>i</mi></munderover><msub><mi>u</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>j</mi><mo>)</mo></mrow><mi>&Delta;</mi><mi>t</mi><mo>&le;</mo><msub><mi>S</mi><mi>B</mi></msub><mo>&lsqb;</mo><msub><mi>SOC</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>SOC</mi><mrow><mi>B</mi><mo>,</mo><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000840905190000022.GIF" wi="1229" he="125" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>S</mi><mrow><mi>u</mi><mi>c</mi></mrow></msub><mo>&lsqb;</mo><msub><mi>SOC</mi><mrow><mi>u</mi><mi>c</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>SOC</mi><mrow><mi>u</mi><mi>c</mi><mo>,</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>&rsqb;</mo><mo>&le;</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mi>i</mi></munderover><msub><mi>u</mi><mrow><mi>u</mi><mi>c</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>j</mi><mo>)</mo></mrow><mi>&Delta;</mi><mi>t</mi><mo>&le;</mo><msub><mi>S</mi><mrow><mi>u</mi><mi>c</mi></mrow></msub><mo>&lsqb;</mo><msub><mi>SOC</mi><mrow><mi>u</mi><mi>c</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>SOC</mi><mrow><mi>u</mi><mi>c</mi><mo>,</mo><mi>min</mi></mrow></msub><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000840905190000023.GIF" wi="1573" he="127" /></maths>其中i=0,1,…,M‑1;S<sub>B</sub>和S<sub>uc</sub>分别是电池和超级电容的容量;SOC<sub>B,max</sub>和SOC<sub>B,min</sub>分别为电池SOC的最大和最小值;SOC<sub>uc,max</sub>和SOC<sub>uc,min</sub>分别为超级电容SOC的最大和最小值;SOC<sub>B</sub>(k)为储能当前时刻的SOC;SOC<sub>uc</sub>(k)为超级电容当前时刻的SOC;Δt为时间间隔;5)将所述目标函数及约束条件化为二次规划标准形进行求解,即:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mi>&theta;</mi></munder><mi>J</mi><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>&theta;</mi><mi>T</mi></msup><mi>&phi;</mi><mi>&theta;</mi></mrow>]]></math><img file="FDA0000840905190000024.GIF" wi="292" he="111" /></maths>s.t.f(θ)≤ω    (7)其中<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>&theta;</mi><mo>=</mo><msup><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msub><mi>u</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msub><mi>u</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd><mtd><mo>...</mo></mtd><mtd><mrow><msub><mi>u</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>M</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>u</mi><mrow><mi>u</mi><mi>c</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msub><mi>u</mi><mrow><mi>u</mi><mi>c</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd><mtd><mo>...</mo></mtd><mtd><mrow><msub><mi>u</mi><mrow><mi>u</mi><mi>c</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>M</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mi>T</mi></msup></mrow>]]></math><img file="FDA0000840905190000025.GIF" wi="814" he="150" /></maths>为优化变量,φ是系数矩阵,ω是一维列向量;φ、ω和函数f由步骤3)中的目标函数和步骤4)的三个约束条件得到;根据约束条件优化每一时刻k的目标函数,得到预测控制时间长度内的M/Δt个电池储能出力的优化变量u<sub>B</sub>(k),u<sub>B</sub>(k+1),…,u<sub>B</sub>(k+M‑1)和M/Δt个超级电容储能出力的优化变量u<sub>uc</sub>(k),u<sub>uc</sub>(k+1),…,u<sub>uc</sub>(k+M‑1),从中仅取得到的当前时刻的优化控制量u<sub>B</sub>(k)和u<sub>uc</sub>(k)分别对电池储能子系统和超级电容储能子系统进行预测控制;6)将每一时刻电池和超级电容输出的功率与风电的实际输出功率相加,即达到整个系统的并网功率波动平抑后的效果。
地址 100084 北京市海淀区清华园1号
您可能感兴趣的专利