发明名称 一种微电网下考虑分布式电源消纳的电动汽车互动响应控制方法
摘要 一种微电网下考虑分布式电源消纳的电动汽车互动响应控制方法,包括如下步骤:S1:将一天连续24h的时间进行离散化处理;S2:由充放电设施记录新入网电动汽车的电池信息和客户充电需求信息,并令初始时段为接入时段;S3:读入当前时刻负荷信息;S4:光伏出力预测;S5:基于光伏出力与负荷供需情况、结合实时电价与倾斜阻塞率发展虚拟电价机制;S6:在虚拟电价的引导下转换最大化光电消纳目标,求解目标函数,制定时长T内电动汽车充放电计划;S7:当前时段下各电动汽车根据控制策略进行用电、闲置或放电操作,更新预测模型信息并将控制信息上传;S8:重复S3~S7直至车辆离开充电设施。本发明光电消纳水平较高、控制效果较好。
申请公布号 CN106230020A 申请公布日期 2016.12.14
申请号 CN201610658639.8 申请日期 2016.08.11
申请人 浙江工业大学 发明人 张有兵;杨晓东;任帅杰;蒋杨昌;谢路耀;翁国庆
分类号 H02J3/38(2006.01)I;B60L11/18(2006.01)I 主分类号 H02J3/38(2006.01)I
代理机构 杭州斯可睿专利事务所有限公司 33241 代理人 王利强
主权项 一种微电网下考虑分布式电源消纳的电动汽车互动响应控制方法,其特征在于:所述控制方法包括以下步骤:步骤1:将一天连续24h的时间进行离散化处理,均分为J个时段,对于任意第k时段,有k∈{1,2,...,J},且第k时段的时长为Δt;步骤2:当电动汽车接入第l(l=1,2,...,n,以下简称电动汽车l)号充放电设施时,充电设施读取电动汽车接入时间,电池的初始状态S<sub>0,l</sub>,且0≤S<sub>0,l</sub>≤1;步骤3:车主输入车辆l的预期离开时间T<sub>out,l</sub>以及离开时期望的荷电状态S<sub>E,l</sub>,且有0≤S<sub>E,l</sub>≤1;步骤4:若电动汽车l持续入网的时长大于将电动汽车l的电池充电至期望电量水平所需的最短时长,则执行步骤5,否则让用户自主选择是否进行修改信息,若用户同意执行修改则跳至步骤3,若用户拒绝执行修改则放弃该用户;步骤5:令初始时段k为车辆接入充电设施的时段;步骤6:读入当前时刻负荷信息,并选择动态优化区间T,根据已有的光伏输出功率的研究结论,以当前时段的光伏出力为起始值,预测未来T时段内的光伏出力;步骤7:基于当前时段微电网内分布式光伏出力与负荷之间的供需情况、结合实时电价与倾斜阻塞率IBR发展了虚拟电价机制,过程如下,步骤7‑1:实时电价机制RTP与系统净负荷关系如下:<maths num="0001"><math><![CDATA[<mrow><msub><mi>RTP</mi><mi>k</mi></msub><mo>=</mo><msub><mi>a</mi><mi>k</mi></msub><msubsup><mi>L</mi><mi>l</mi><mi>k</mi></msubsup><mo>+</mo><msub><mi>b</mi><mi>k</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001075586260000011.GIF" wi="942" he="63" /></maths><maths num="0002"><math><![CDATA[<mrow><msubsup><mi>L</mi><mi>l</mi><mi>k</mi></msubsup><mo>=</mo><mrow><mo>(</mo><msubsup><mi>P</mi><mi>B</mi><mi>k</mi></msubsup><mo>+</mo><msubsup><mi>L</mi><mrow><mi>E</mi><mi>V</mi></mrow><mrow><mi>l</mi><mo>,</mo><mi>k</mi></mrow></msubsup><mo>+</mo><msubsup><mi>P</mi><mrow><mi>E</mi><mi>S</mi><mi>S</mi></mrow><mi>k</mi></msubsup><mo>-</mo><msubsup><mi>P</mi><mrow><mi>P</mi><mi>V</mi></mrow><mi>k</mi></msubsup><mo>)</mo></mrow><mi>&Delta;</mi><mi>t</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001075586260000012.GIF" wi="1046" he="67" /></maths><maths num="0003"><math><![CDATA[<mrow><msubsup><mi>L</mi><mrow><mi>E</mi><mi>V</mi></mrow><mrow><mi>l</mi><mo>,</mo><mi>k</mi></mrow></msubsup><mo>=</mo><munder><mo>&Sigma;</mo><mrow><mi>i</mi><mo>&Element;</mo><msub><mi>M</mi><mrow><mi>l</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></munder><msubsup><mi>P</mi><mi>i</mi><mi>k</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001075586260000013.GIF" wi="916" he="95" /></maths>式中:RTP<sub>k</sub>为k时段的实时电价;<img file="FDA0001075586260000021.GIF" wi="51" he="70" />为电动汽车l接入时,k时段的微电网系统净负荷;a<sub>k</sub>、b<sub>k</sub>为实时电价系数,在不同的时间段取不同的值,取决于用户的需求动态;<img file="FDA0001075586260000022.GIF" wi="58" he="71" />为系统基本负荷,即该微电网中除电动汽车集群负荷之外的所有电力负荷;<img file="FDA0001075586260000023.GIF" wi="86" he="70" />为储能系统在k时段的充放电功率;<img file="FDA0001075586260000024.GIF" wi="70" he="70" />为k时段的光伏出力功率;<img file="FDA0001075586260000025.GIF" wi="78" he="70" />表示车辆l接入微电网时,充放电计划制定已完成的车辆集群负荷;M<sub>l‑1</sub>表示车辆l接入微电网时,充放电计划已完成的车辆结合;步骤7‑2:本发明在IBR中设置三种电价等级:<maths num="0004"><math><![CDATA[<mrow><msub><mi>IBR</mi><mi>k</mi></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>x</mi><mi>k</mi></msub><mo>,</mo></mrow></mtd><mtd><mrow><mn>0</mn><mo>&le;</mo><msubsup><mi>L</mi><mi>l</mi><mi>k</mi></msubsup><mo>&le;</mo><msubsup><mi>&delta;</mi><mi>k</mi><mn>1</mn></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>y</mi><mi>k</mi></msub><mo>,</mo></mrow></mtd><mtd><mrow><msubsup><mi>&delta;</mi><mi>k</mi><mn>1</mn></msubsup><mo>&le;</mo><msubsup><mi>L</mi><mi>l</mi><mi>k</mi></msubsup><mo>&le;</mo><msubsup><mi>&delta;</mi><mi>k</mi><mn>2</mn></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>z</mi><mi>k</mi></msub><mo>,</mo></mrow></mtd><mtd><mrow><msubsup><mi>L</mi><mi>l</mi><mi>k</mi></msubsup><mo>&GreaterEqual;</mo><msubsup><mi>&delta;</mi><mi>k</mi><mn>2</mn></msubsup></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001075586260000026.GIF" wi="1069" he="211" /></maths>式中,<img file="FDA0001075586260000027.GIF" wi="52" he="70" />与<img file="FDA0001075586260000028.GIF" wi="54" he="70" />为不同电价等级之间的界限;x<sub>k</sub>、y<sub>k</sub>与z<sub>k</sub>为三个等级下的电价,具体计算方法如下:<maths num="0005"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>x</mi><mi>k</mi></msub><mo>=</mo><msub><mi>RTP</mi><mi>k</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>y</mi><mi>k</mi></msub><mo>=</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>*</mo><msub><mi>x</mi><mi>k</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>z</mi><mi>k</mi></msub><mo>=</mo><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>*</mo><msub><mi>x</mi><mi>k</mi></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001075586260000029.GIF" wi="933" he="183" /></maths>式中,λ<sub>1</sub>与λ<sub>2</sub>为不同等级下的价格倍率,并且λ<sub>2</sub>>λ<sub>1</sub>>1;步骤7‑3:综上所述,虚拟电价<img file="FDA00010755862600000210.GIF" wi="153" he="71" />的计算方式为:<maths num="0006"><math><![CDATA[<mrow><mi>p</mi><mi>r</mi><mrow><mo>(</mo><msubsup><mi>L</mi><mi>l</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>RTP</mi><mrow><mi>r</mi><mi>e</mi></mrow></msub><mo>,</mo></mrow></mtd><mtd><mrow><msubsup><mi>L</mi><mi>l</mi><mi>k</mi></msubsup><mo>&le;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>RTP</mi><mi>k</mi></msub><mo>,</mo></mrow></mtd><mtd><mrow><mn>0</mn><mo>&lt;</mo><msubsup><mi>L</mi><mi>l</mi><mi>k</mi></msubsup><mo>&le;</mo><msubsup><mi>&delta;</mi><mi>k</mi><mn>1</mn></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>*</mo><msub><mi>RTP</mi><mi>k</mi></msub><mo>,</mo></mrow></mtd><mtd><mrow><msubsup><mi>&delta;</mi><mi>k</mi><mn>1</mn></msubsup><mo>&lt;</mo><msubsup><mi>L</mi><mi>l</mi><mi>k</mi></msubsup><mo>&le;</mo><msubsup><mi>&delta;</mi><mi>k</mi><mn>2</mn></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>*</mo><msub><mi>RTP</mi><mi>k</mi></msub><mo>,</mo></mrow></mtd><mtd><mrow><msubsup><mi>L</mi><mi>l</mi><mi>k</mi></msubsup><mo>&gt;</mo><msubsup><mi>&delta;</mi><mi>k</mi><mn>2</mn></msubsup></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00010755862600000211.GIF" wi="1086" he="263" /></maths>式中,<img file="FDA00010755862600000212.GIF" wi="134" he="70" />时,意味着可再生能源出力过剩,此时,多余的光伏发电量向上级电网倒送,RTP<sub>re</sub>为单位电量的倒送价格;步骤8:在虚拟电价的引导下转换最大化光电消纳目标,制定时长T内电动汽车充放电计划,其目标函数的制定过程如下,基于步骤7所述的虚拟电价模型、以充放电虚拟总成本最小为目标对其进行动态规划:<maths num="0007"><math><![CDATA[<mrow><mi>min</mi><mi> </mi><msub><mi>V</mi><mi>l</mi></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>J</mi></munderover><msubsup><mi>X</mi><mi>l</mi><mi>k</mi></msubsup><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>J</mi></munderover><mi>p</mi><mi>r</mi><mrow><mo>(</mo><msubsup><mi>L</mi><mi>l</mi><mi>k</mi></msubsup><mo>)</mo></mrow><msubsup><mi>P</mi><mi>l</mi><mi>k</mi></msubsup><mi>&Delta;</mi><mi>t</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001075586260000031.GIF" wi="1095" he="112" /></maths>式中,V<sub>l</sub>为电动汽车l的虚拟总成本;P<sub>l</sub><sup>k</sup>为k时段电动汽车l与微电网的交换功率,P<sub>l</sub><sup>k</sup>>0表示车辆l的充电功率;P<sub>l</sub><sup>k</sup><0表示放电功率;P<sub>l</sub><sup>k</sup>=0表示处于闲置状态,EV动力电池模型和约束条件为:<maths num="0008"><math><![CDATA[<mrow><msubsup><mi>S</mi><mi>l</mi><mi>k</mi></msubsup><mo>=</mo><msubsup><mi>S</mi><mi>l</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>+</mo><msubsup><mi>P</mi><mi>l</mi><mi>k</mi></msubsup><msub><mi>&eta;</mi><mrow><mi>E</mi><mi>V</mi></mrow></msub><mi>&Delta;</mi><mi>t</mi><mo>/</mo><msub><mi>Q</mi><mrow><mi>E</mi><mi>V</mi><mo>,</mo><mi>l</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001075586260000032.GIF" wi="1014" he="71" /></maths><maths num="0009"><math><![CDATA[<mrow><msub><mi>S</mi><mrow><mi>E</mi><mi>V</mi><mo>,</mo><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub><mo>&le;</mo><msubsup><mi>S</mi><mi>l</mi><mi>k</mi></msubsup><mo>&le;</mo><msub><mi>S</mi><mrow><mi>E</mi><mi>V</mi><mo>,</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001075586260000033.GIF" wi="1005" he="67" /></maths>‑P<sub>EV,d</sub>≤P<sub>l</sub><sup>k</sup>≤P<sub>EV,c</sub>    (10)<maths num="0010"><math><![CDATA[<mrow><msub><mi>S</mi><mrow><mn>0</mn><mo>,</mo><mi>l</mi></mrow></msub><mo>+</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>J</mi></munderover><msubsup><mi>P</mi><mi>l</mi><mi>k</mi></msubsup><msub><mi>&eta;</mi><mrow><mi>E</mi><mi>V</mi></mrow></msub><mi>&Delta;</mi><mi>t</mi></mrow><msub><mi>Q</mi><mrow><mi>E</mi><mi>V</mi><mo>,</mo><mi>l</mi></mrow></msub></mfrac><mo>&GreaterEqual;</mo><msub><mi>S</mi><mrow><mi>E</mi><mo>,</mo><mi>l</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001075586260000034.GIF" wi="1022" he="171" /></maths><maths num="0011"><math><![CDATA[<mrow><msub><mi>&eta;</mi><mrow><mi>E</mi><mi>V</mi></mrow></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>&eta;</mi><mi>c</mi></msub><mo>,</mo></mrow></mtd><mtd><mrow><msubsup><mi>P</mi><mi>l</mi><mi>k</mi></msubsup><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo></mrow></mtd></mtr><mtr><mtd><mrow><mn>1</mn><mo>/</mo><msub><mi>&eta;</mi><mi>d</mi></msub><mo>,</mo></mrow></mtd><mtd><mrow><msubsup><mi>P</mi><mi>l</mi><mi>k</mi></msubsup><mo>&lt;</mo><mn>0.</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001075586260000035.GIF" wi="1014" he="135" /></maths>式中,<img file="FDA0001075586260000036.GIF" wi="210" he="80" />分别为车辆l在k时段和k‑1时段的电池SOC;Q<sub>EV,l</sub>为车辆动力电池容量;S<sub>EV,max</sub>、S<sub>EV,min</sub>分别为动力电池SOC的上、下限;η<sub>EV</sub>表示电池功率交换效率,与功率交换方向有关,如式(12)所示;η<sub>c</sub>、η<sub>d</sub>分别表示充、放电效率;还需考虑微电网系统功率平衡约束和倒送功率约束其式为:<maths num="0012"><math><![CDATA[<mrow><msubsup><mi>P</mi><mrow><mi>p</mi><mi>v</mi></mrow><mi>k</mi></msubsup><mo>+</mo><msubsup><mi>P</mi><mrow><mi>g</mi><mi>r</mi><mi>i</mi><mi>d</mi></mrow><mi>k</mi></msubsup><mo>=</mo><msubsup><mi>P</mi><mi>B</mi><mi>k</mi></msubsup><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><msubsup><mi>P</mi><mi>l</mi><mi>k</mi></msubsup><mo>+</mo><msubsup><mi>P</mi><mrow><mi>E</mi><mi>S</mi><mi>S</mi></mrow><mi>k</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001075586260000037.GIF" wi="1014" he="111" /></maths><maths num="0013"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>g</mi><mi>r</mi><mi>i</mi><mi>d</mi><mi>o</mi><mi>u</mi><mi>t</mi></mrow></msub><mo>&le;</mo><msubsup><mi>P</mi><mrow><mi>g</mi><mi>r</mi><mi>i</mi><mi>d</mi><mi>o</mi><mi>u</mi><mi>t</mi></mrow><mi>max</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001075586260000038.GIF" wi="916" he="71" /></maths>式中,<img file="FDA0001075586260000039.GIF" wi="86" he="76" />表示k时段微电网与大电网的交互功率,<img file="FDA00010755862600000310.GIF" wi="118" he="77" />为微电网的倒送功率;<img file="FDA00010755862600000311.GIF" wi="120" he="71" />为倒送功率允许的最大值;步骤9:在当前时段,各电动汽车根据控制策略进行具体的用电、闲置或放电操作,同时,更新预测模型信息并将控制信息上传至电能公共服务平台;步骤10:当进入新时段时重复步骤6~9直至车辆离开充电设施。
地址 310014 浙江省杭州市下城区朝晖六区潮王路18号
您可能感兴趣的专利