发明名称 一种基于虚拟电价的电动汽车充放电优化调度方法
摘要 一种基于虚拟电价的电动汽车充放电优化调度方法,包括:电能公共服务平台预测、采样优化时段内目标区域的基本日负荷信息;当有新的EV接入目标区域内的充电桩时,读取其入网信息;用户输入车辆的充电信息;构建EV充放电功率模型;计算虚拟电价信息,间接反映目标区域的负荷水平;构建以充放电功率为优化变量的调度模型;综合小波分析预处理和模糊聚类方法确定用于用户成本计算的动态分时电价;用户的自主响应决策;根据用户决策对EV实施充放电操作,并上传计划。本发明能够在满足用户充电需求及配电变压器容量限制的基础上,实现EV集群负荷的削峰填谷,并降低用户充放电成本。在EV集群规模较大时,本发明仍能满足电网侧期望。
申请公布号 CN105024432A 申请公布日期 2015.11.04
申请号 CN201510458366.8 申请日期 2015.07.30
申请人 浙江工业大学 发明人 张有兵;杨晓东;任帅杰;翁国庆;周文委;谢路耀
分类号 H02J7/00(2006.01)I;H02J7/02(2006.01)I;G06Q10/04(2012.01)I;G06Q50/06(2012.01)I 主分类号 H02J7/00(2006.01)I
代理机构 杭州天正专利事务所有限公司 33201 代理人 王兵;黄美娟
主权项 一种基于虚拟电价的电动汽车充放电优化调度方法,每台充电桩的充电功率额定值为P<sub>c</sub>,放电功率额定值为P<sub>d</sub>,包括如下步骤:1)电能公共服务平台预测、采样优化时段内目标区域的基本日负荷信息,具体包括:S11.预测基本负荷功率曲线:根据已知参数和模型,预测出目标区域的基本日负荷曲线;S12.优化时段内的连续时间离散化:确定最小优化时段Δt,对优化时段进行离散化分析,依次对步骤S11所得的基本日负荷曲线进行采样,得各k∈[1,J]时段内的负荷预测值L<sub>B</sub>(k);其中J为一日内根据设定的确定划分的优化时段数;2)当有新的电动汽车l接入目标区域内的充电桩时,读取车辆l的入网信息,包括入网时间T<sub>in,l</sub>、动力电池容量C<sub>s,l</sub>和初始荷电状态(State of Charge,SOC)S<sub>0,l</sub>等;3)用户输入车辆l的充电信息:预期离网时间T<sub>out,l</sub>,期望荷电状态S<sub>E,l</sub>;结合所述读取信息和所述用户输入信息判断能否满足用户的充电需求,如果不能,则通过用户交互界面告知用户输入不正确,提示用户重新输入预期离网时间或期望SOC信息,若用户接受提示,并输入可实施、正确的充电信息,则进行步骤4)~步骤10)的操作,否则,放弃该用户,执行步骤2);4)构建电动汽车的充放电功率模型:P<sub>l</sub>(k)=p<sub>l</sub>(k)f<sub>m,l</sub>(k);式中,P<sub>l</sub>(k)表示车辆l的充放电功率;p<sub>l</sub>(k)表示车辆l与系统在k时段的功率交换,p<sub>l</sub>(k)&gt;0表示车辆l处于充电状态;p<sub>l</sub>(k)&lt;0表示处于放电状态;p<sub>l</sub>(k)=0表示处于浮充状态;f<sub>m,l</sub>(k)用于表征各时段对车辆电池的可操作性,表达式为:<img file="FDA0000771250550000011.GIF" wi="1574" he="140" />其中,T<sub>m,l</sub>为车辆l接入电网的持续时间T<sub>pe,l</sub>=T<sub>out,l</sub>‑T<sub>in,l</sub>所包含的时段集合;假设参与调度的电动汽车动力电池均为锂电池,根据锂电池的充放电相关特性,在单个时段内,锂电可视为是恒功率充放电,其荷电状态和对应的充放电时间的关系表征为:S<sub>l</sub>(k)=S<sub>l</sub>(k‑1)+P<sub>l</sub>(k)η(P<sub>l</sub>(k))Δt/C<sub>s,l</sub>   (2)式中,S<sub>l</sub>(k‑1)、S<sub>l</sub>(k)分别表示车辆l在第k‑1、k个时段的SOC;η(P<sub>l</sub>(k))表示功率交换效率,具体与功率交换方向有关:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>&eta;</mi><mrow><mo>(</mo><msub><mi>P</mi><mi>l</mi></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>=</mo><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><msub><mi>&eta;</mi><mi>c</mi></msub><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>P</mi><mi>l</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><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><msub><mi>P</mi><mi>l</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&lt;</mo><mn>0.</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000771250550000012.GIF" wi="1595" he="139" /></maths>其中,η<sub>c</sub>、η<sub>d</sub>分别表示充、放电效率;5)从电能公共服务平台读入车辆l接入时刻目标区域的总负荷信息<img file="FDA0000771250550000015.GIF" wi="127" he="67" /><img file="FDA0000771250550000016.GIF" wi="416" he="67" />其中,<img file="FDA0000771250550000017.GIF" wi="109" he="67" />表示车辆l接入时,k时段目标区域的配网总负荷信息;<img file="FDA0000771250550000013.GIF" wi="126" he="75" />表示电动汽车集群负荷:<img file="FDA0000771250550000014.GIF" wi="663" he="110" />其中,M<sub>l‑1</sub>表示车辆l接入电网时,已完成充电计划制定的车辆集合;该计算过程由电能公共服务平台完成;6)计算虚拟电价信息,间接反映目标区域的负荷水平:<img file="FDA0000771250550000021.GIF" wi="1794" he="312" />式中,<img file="FDA00007712505500000211.GIF" wi="214" he="67" />表示车辆l接入时,k时段的虚拟电价;<img file="FDA0000771250550000022.GIF" wi="151" he="57" />和<img file="FDA0000771250550000023.GIF" wi="48" he="48" />表示虚拟电价调整系数;I<sub>R,j</sub>、φ<sub>R,j</sub>分别表示参考电价和参考负荷值;[u]<sup>+</sup>表示max{0,u};<img file="FDA0000771250550000024.GIF" wi="424" he="84" />表示预测总负荷,其中<img file="FDA0000771250550000025.GIF" wi="104" he="80" />表示基本负荷预测值;7)建立电动汽车充放电优化调度模型;结合以上步骤,以最小化虚拟充放电成本为目标,建立电动汽车充放电优化调度模型来优化电动汽车的充放电功率,所建模型如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>min</mi><msub><mi>V</mi><mi>l</mi></msub><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>k</mi><mo>&Element;</mo><msub><mi>T</mi><mrow><mi>m</mi><mo>,</mo><mi>l</mi></mrow></msub></mrow></munder><mo>{</mo><msub><mi>P</mi><mi>l</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>&gamma;</mi><mi>l</mi></msub><mrow><mo>(</mo><msubsup><mi>L</mi><mi>a</mi><mi>l</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mi>&Delta;t</mi><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000771250550000026.GIF" wi="1603" he="110" /></maths>s.t. S<sub>min</sub>≤S<sub>l</sub>(k)≤S<sub>max</sub>   (7)‑P<sub>d</sub>≤P<sub>l</sub>(k)≤P<sub>c</sub>   (8)<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>S</mi><mrow><mn>0</mn><mo>,</mo><mi>l</mi></mrow></msub><mo>+</mo><munder><mo>&Sigma;</mo><mrow><mi>k</mi><mo>&Element;</mo><msub><mi>T</mi><mrow><mi>m</mi><mo>,</mo><mi>l</mi></mrow></msub></mrow></munder><mi>&Delta;</mi><mi>t</mi><mfrac><mrow><msub><mi>P</mi><mi>l</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><msub><mi>C</mi><mrow><mi>s</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>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000771250550000027.GIF" wi="1476" he="125" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>L</mi><mi>a</mi><mi>l</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&le;</mo><msub><mi>&kappa;</mi><mi>T</mi></msub><msub><mi>A</mi><mi>T</mi></msub><mo>,</mo><mover><mrow><msubsup><mi>L</mi><mi>a</mi><mi>l</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mo>^</mo></mover><mo>&le;</mo><msub><mi>&kappa;</mi><mi>T</mi></msub><msub><mi>A</mi><mi>T</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000771250550000028.GIF" wi="1453" he="94" /></maths>k=1,2,…,J.T<sub>pe,l</sub>&gt;T<sub>c,l</sub> l=1,2,…,n   (11)式(6)中,V<sub>l</sub>表示车辆l的虚拟充放电成本;N<sub>l</sub>表示T<sub>m,l</sub>集合的长度;式(7)中,S<sub>max</sub>、S<sub>min</sub>表示允许的SOC的最大值和最小值,防止受控车辆过充电及过放电;式(8)表示充放电功率约束,P<sub>l</sub>(k)具备连续可调的特性,但通常要受到动力电池或充电机的额定充、放电功率限制;式(9)表示充电需求约束,车辆如约离开时,其电池的SOC需满足期望;式(10)表示变压器最大负载约束,κ<sub>T</sub>为变压器效率,A<sub>T</sub>为变压器的额定容量;式(11)表示时间关系约束,n为优化时段内的接入车辆规模,T<sub>c,l</sub>为充电至期望SOC所需的最短时间:T<sub>c,l</sub>=(S<sub>E,l</sub>‑S<sub>0,l</sub>)C<sub>s,l</sub>/P<sub>c</sub>η<sub>c</sub>;求解上述优化调度模型,完成当前接入车辆充放电功率的优化,此时,车辆l的充放电调度计划为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>P</mi><mi>l</mi></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>P</mi><mi>l</mi></msub><mo>(</mo><msubsup><mi>T</mi><mrow><mi>m</mi><mo>,</mo><mi>l</mi></mrow><mn>1</mn></msubsup><mo>)</mo></mrow><mo>,</mo><msub><mi>P</mi><mi>l</mi></msub><mrow><mo>(</mo><msubsup><mi>T</mi><mrow><mi>m</mi><mo>,</mo><mi>l</mi></mrow><mn>2</mn></msubsup><mo>)</mo></mrow><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>P</mi><mi>l</mi></msub><mrow><mo>(</mo><msubsup><mi>T</mi><mrow><mi>m</mi><mo>,</mo><mi>l</mi></mrow><msub><mi>N</mi><mi>l</mi></msub></msubsup><mo>)</mo></mrow><mo>)</mo><mo>,</mo></mrow>]]></math><img file="FDA0000771250550000029.GIF" wi="607" he="76" /></maths>其中,<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mi>T</mi><mrow><mi>m</mi><mo>,</mo><mi>l</mi></mrow><mi>i</mi></msubsup><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>N</mi><mi>l</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA00007712505500000210.GIF" wi="348" he="81" /></maths>表示集合T<sub>m,l</sub>中的第i个元素;8)电动汽车用户成本计算:U<sub>l</sub>=(c<sub>cd,l</sub>|η<sub>c</sub>=η<sub>d</sub>=1)+c<sub>bat,l</sub>+c<sub>loss,l</sub>;式中,U<sub>l</sub>表示用户成本;c<sub>bat,l</sub>表示车辆l的锂电池寿命损耗折算成本;c<sub>loss,l</sub>表示电能损失费用;c<sub>cd,l</sub>|η<sub>c</sub>=η<sub>d</sub>=1,表示不考虑充、放电效率条件下的理想充放电费用:<img file="FDA0000771250550000031.GIF" wi="414" he="106" />式中pri(k)表示电价信息,本发明中,其是一种动态分时电价,即高峰、低谷电价固定,峰谷时段变化的电价,其中,高峰、低谷电价分别表示为pri<sup>h</sup>、pri<sup>l</sup>;利用小波分析和模糊聚类方法实现对分时电价峰谷时段的划分:S81.预处理:将负荷信息进行尺度为3的小波分解,将一、二层的高频分量置零,重构后得到新的负荷信息<img file="FDA0000771250550000032.GIF" wi="141" he="73" />S82.属性表征:采用偏大型和偏小型半梯形模糊分布对重构后的负荷信息进行属性表征,形成属性矩阵A(a)<sub>J×2</sub>,计算方法为:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><msub><mi>a</mi><mrow><mi>k</mi><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mrow><munder><mi>max</mi><mrow><mn>1</mn><mo>&le;</mo><mi>k</mi><mo>&le;</mo><mi>J</mi></mrow></munder><mrow><mo>(</mo><msubsup><mi>L</mi><mi>a</mi><mrow><mi>l</mi><mo>*</mo></mrow></msubsup><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>-</mo><msubsup><mi>L</mi><mi>a</mi><mrow><mi>l</mi><mo>*</mo></mrow></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><munder><mi>max</mi><mrow><mn>1</mn><mo>&le;</mo><mi>k</mi><mo>&le;</mo><mi>J</mi></mrow></munder><mrow><mo>(</mo><msubsup><mi>L</mi><mi>a</mi><mrow><mi>l</mi><mo>*</mo></mrow></msubsup><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>-</mo><munder><mi>max</mi><mrow><mn>1</mn><mo>&le;</mo><mi>k</mi><mo>&le;</mo><mi>J</mi></mrow></munder><mrow><mo>(</mo><msubsup><mi>L</mi><mi>a</mi><mrow><mi>l</mi><mo>*</mo></mrow></msubsup><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>a</mi><mrow><mi>k</mi><mn>2</mn></mrow></msub><mo>=</mo><mfrac><mrow><msubsup><mi>L</mi><mi>a</mi><mrow><mi>l</mi><mo>*</mo></mrow></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><munder><mi>max</mi><mrow><mn>1</mn><mo>&le;</mo><mi>k</mi><mo>&le;</mo><mi>J</mi></mrow></munder><mrow><mo>(</mo><msubsup><mi>L</mi><mi>a</mi><mrow><mi>l</mi><mo>*</mo></mrow></msubsup><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow></mrow><mrow><munder><mi>max</mi><mrow><mn>1</mn><mo>&le;</mo><mi>k</mi><mo>&le;</mo><mi>J</mi></mrow></munder><mrow><mo>(</mo><msubsup><mi>L</mi><mi>a</mi><mrow><mi>l</mi><mo>*</mo></mrow></msubsup><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>-</mo><munder><mi>max</mi><mrow><mn>1</mn><mo>&le;</mo><mi>k</mi><mo>&le;</mo><mi>J</mi></mrow></munder><mrow><mo>(</mo><msubsup><mi>L</mi><mi>a</mi><mrow><mi>l</mi><mo>*</mo></mrow></msubsup><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow></mrow></mfrac></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000771250550000033.GIF" wi="1614" he="326" /></maths>S83.对矩阵A进行平移—标准差变换,用绝对值减数法建立模糊相似矩阵R(r)<sub>J×J</sub>;S84.对相似矩阵求连续二次方,即R→R<sup>2</sup>→R<sup>4</sup>→…→R<sup>2i</sup>→…,直到出现<img file="FDA0000771250550000034.GIF" wi="204" he="68" />(<img file="FDA0000771250550000035.GIF" wi="28" he="33" />表示模糊矩阵合成运算),此时,R<sup>k</sup>即为该相似矩阵的传递闭包t(R);S85.在传递闭包t(R)=(t<sub>ij</sub>)中,令λ为t<sub>ij</sub>中的某一取值,求取t(R)的λ—截矩阵R<sub>ij</sub>,接着,λ由大到小取值就可形成动态聚类,取适当的λ值,确定出峰谷时段划分方案;9)电动汽车用户的自主响应决策;利用步骤7)所述的优化模型得到的电动汽车充放电功率和步骤8)得到的动态分时电价确定用户的充放电成本,并将充放电调度计划与对应的收益情况告知用户,由用户自主响应充电模式,选择响应调度计划或开始无序充电;10)根据用户决策对电动汽车实施充放电操作,并上传计划;若用户选择开始无序充电,充电设施为接入的EV提供持续的恒功率充电服务,直到满足用户的充电需求或车辆离开为止;若用户选择响应调度,则根据调度计划P<sub>l</sub>对电动汽车实施具体充放电操作,据此,便确定了车辆l的充放电计划;将该充放电计划上传至电能公共服务平台,由电能公共服务平台对计划负荷予以整合,完成一次对所述目标区域负荷信息的实时更新,并等待下一辆电动汽车入网;若有新的车辆接入,则跳至步骤2);此过程将持续执行至优化时段结束。
地址 310014 浙江省杭州市下城区潮王路18号