发明名称 基于V2G技术的电动汽车参与经济调度优化控制方法
摘要 基于V2G技术的电动汽车参与经济调度优化控制方法提出了一个多目标优化模型,将电动汽车车主成本和经济调度成本作为多目标模型的两个目标函数,并让电动汽车通过有序充放电的方式来作为经济调度时的备用容量。在满足各种约束条件下,采用多目标遗传算法(NSGA‑II)对模型进行求解。通过选择合适的pareto解集中的值,电动汽车采用V2G技术参与经济调度后可以节省车主成本和经济调度的成本,并且可以实现对负荷削峰填谷的功能。与不含电动汽车的经济调度有所区别,本发明充分考虑了电动汽车车主的充电成本和电动汽车集群的负荷特性,能提高电力系统运行的经济性和稳定性。
申请公布号 CN104600729B 申请公布日期 2016.12.07
申请号 CN201410407390.4 申请日期 2014.08.19
申请人 浙江工业大学 发明人 张有兵;陆建丽;杨晓东;任帅杰;周文委;黄直;蒋洁毅;翁国庆
分类号 H02J3/32(2006.01)I 主分类号 H02J3/32(2006.01)I
代理机构 杭州天正专利事务所有限公司 33201 代理人 王兵;黄美娟
主权项 一种基于V2G技术的电动汽车参与经济调度优化控制方法,包括以下步骤:步骤1,向电动汽车充放电控制和经济调度控制总系统输入各种初始信息,包括电动汽车车辆集群的车辆信息、电动汽车动力电池信息、参与经济调度的各台发电机以及风机的参数信息、各种价格信息、24小时各个时段的负荷预测值;所述的电动汽车车辆集群的车辆信息包括:车辆数量N,电动汽车充满电时的里程数R<sub>i</sub>,车辆接入电网时间T<sub>i</sub><sup>1</sup>,车辆离开电网时间T<sub>i</sub><sup>2</sup>,车辆日行驶里程D<sub>i</sub>,电动汽车充放电功率PC<sub>i,t</sub>,i=1,2,3...N,t=1,2,3...T,T为研究周期时间段数;所述的电动汽车动力电池信息包括:电池容量S<sub>i</sub>,电池最大充电上限E<sub>max,i</sub>,电池最小放电下限E<sub>min,i</sub>,电池消耗系数ε<sub>1</sub>,ε<sub>2</sub>,ε<sub>3</sub>,ε<sub>4</sub>,δ,充电和放电效率系数η<sub>c,i</sub>、η<sub>d,i</sub>,i=1,2,3...N;所述的参与经济调度的各台发电机以及风机的参数信息包括:风机的形状参数α,风机的尺度参数β,风机额定输出功率P<sub>r</sub>,风机的切入风速v<sub>ci</sub>,风机的额定风速v<sub>r</sub>,风机的切出风速v<sub>co</sub>,参与调频的机组数n,机组i的成本系数a<sub>i</sub>、b<sub>i</sub>、c<sub>i</sub>,机组i的NO<sub>X</sub>排放量的特征系数α<sub>i</sub>、β<sub>i</sub>、θ<sub>i</sub>、δ<sub>i</sub>、λ<sub>i</sub>,机组出力的最大值和最小值P<sub>gi,max</sub>、P<sub>gi,min</sub>,机组i的强迫停运率γ<sub>i</sub>,i=1,2...n;所述的各种价格信息包括:购买备用容量价格p<sub>r,t</sub>,实际调用备用容量价格p<sub>re,t</sub>,系统购买电动汽车备用价格p<sub>v,t</sub>,失负荷价格V<sub>oll,t</sub>,污染排放价格p<sub>pol,t</sub>,弃风价格p<sub>e,t</sub>,车主支付的充电价格p<sub>g2v,t</sub>,车主获利的放电价格p<sub>v2g,t</sub>;步骤2,根据电动汽车集群的车辆信息得到每辆车的初始荷电状态,采用蒙特卡洛算法对各个时刻的风能、负荷波动和机组停运状态进行预测;每辆车的初始荷电状态E<sub>0,i</sub>的计算公式<img file="FDA0000976029500000011.GIF" wi="434" he="139" />假设风速服从Weibull分布,即<img file="FDA0000976029500000012.GIF" wi="709" he="127" />其中,v<sub>t</sub>是t时刻的风速;由风速预测值可以求得实际风能P<sub>w,t</sub>输出,即<img file="FDA0000976029500000013.GIF" wi="789" he="347" />假定t时刻的负荷波动ΔP<sub>l,t</sub>服从正太分布<img file="FDA0000976029500000014.GIF" wi="267" he="79" />σ<sub>l</sub>为其标准差;t时刻机组i的停运状态数为d<sub>i,t</sub>,随机产生一个服从[0,1]内均匀分布的伪随机数ω<sub>i</sub>,若ω<sub>i</sub>≤γ<sub>i</sub>,则机组i的停运状态为0,表示机组故障停运,否则机组的停运状态为1,表示机组正常运行;步骤3,根据电动汽车车辆集群的车辆信息、电动汽车动力电池信息、参与经济调度的各台发电机以及风机的参数信息、每辆车的初始荷电状态列出电动汽车充放电、发电机和风机出力以及备用购买所需要满足的约束条件,随机初始化个体数为X的满足这些约束条件的群体P;初始化种群,随机产生X个群体P;群体P中每个个体包括的变量有t时刻各台发电机组的出力P<sub>gi,t</sub>、计划风能输出P<sub>ew,t</sub>、购买备用容量D<sub>r,t</sub>、电动汽车充放电控制矩阵C<sub>M×T</sub>;充放电控制矩阵C<sub>M×T</sub>中,C<sub>ij</sub>为第i量电动汽车在当日第j个时间段其电池的充放电功率,C<sub>ij</sub>&gt;0表示此刻该辆车处于充电状态,C<sub>ij</sub>&lt;0表示此刻该辆车处于放电状态;判断群体P的每个个体是否满足如下约束条件:1.系统功率平衡约束:<img file="FDA0000976029500000021.GIF" wi="685" he="159" />式中,D<sub>l,t</sub>是不包含电动汽车的t时刻的负荷量,D<sub>vl,t</sub>为t时刻车辆集群的总充电量;2.发电机运行约束、风机出力约束和购买备用约束:P<sub>gi,min</sub>≤P<sub>gi,t</sub>≤P<sub>gi,max</sub>0≤P<sub>w,t</sub>≤P<sub>r</sub>D<sub>r,t</sub>≥max{P<sub>gi,t</sub>}3.车辆充放电约束:E<sub>min,i</sub>≤E<sub>i,t</sub>+SD<sub>g2v,i,t</sub>·η<sub>c,i</sub>/S<sub>i</sub>‑SD<sub>v2g,i,t</sub>/η<sub>d,i</sub>/S<sub>i</sub>≤E<sub>max,i</sub>式中,E<sub>i,t</sub>是车辆i在t时段的SOC;SD<sub>v2g,i,t</sub>和SD<sub>g2v,i,t</sub>分别为车辆i在t时段的充放电量;4.车辆充电需求约束:<maths num="0001"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><mo>&lsqb;</mo><msub><mi>SD</mi><mrow><mi>g</mi><mn>2</mn><mi>v</mi><mo>,</mo><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>&eta;</mi><mrow><mi>c</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>-</mo><msub><mi>SD</mi><mrow><mi>v</mi><mn>2</mn><mi>g</mi><mo>,</mo><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>/</mo><msub><mi>&eta;</mi><mrow><mi>d</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>&rsqb;</mo><mo>=</mo><msub><mi>S</mi><mi>i</mi></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><mi>E</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>-</mo><msub><mi>E</mi><mrow><mn>0</mn><mo>,</mo><mi>i</mi></mrow></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000976029500000022.GIF" wi="1318" he="166" /></maths>不满足上述约束条件的群体P中的个体将被重新初始化,直到满足上述约束条件,于是就可以得到最终生成的X个初始化群体P;步骤4,根据初始化的群体P计算出电动汽车用户总成本;电动汽车用户成本模型中包括充电成本、放电获益和电池损耗成本,计算公式如下:<maths num="0002"><math><![CDATA[<mrow><mi>min</mi><mi> </mi><msub><mi>F</mi><mn>1</mn></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><mo>&lsqb;</mo><msub><mi>p</mi><mrow><mi>g</mi><mn>2</mn><mi>v</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>SD</mi><mrow><mi>g</mi><mn>2</mn><mi>v</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>-</mo><msub><mi>p</mi><mrow><mi>v</mi><mn>2</mn><mi>g</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>SD</mi><mrow><mi>v</mi><mn>2</mn><mi>g</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>D</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mrow><mo>(</mo><msub><mi>I</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>,</mo><msub><mi>I</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000976029500000023.GIF" wi="1316" he="139" /></maths>D<sub>i,t</sub>(I<sub>i,t</sub>,I<sub>i,t‑1</sub>)=ε<sub>1</sub>(I<sub>i,t</sub>)<sup>2</sup>+ε<sub>2</sub>(I<sub>i,t</sub>‑I<sub>i,t‑1</sub>)<sup>2</sup>+ε<sub>3</sub>(min(E<sub>i,t</sub>‑δ·S<sub>i</sub>,0))<sup>2</sup>+ε<sub>4</sub><maths num="0003"><math><![CDATA[<mrow><msub><mi>SD</mi><mrow><mi>g</mi><mn>2</mn><mi>v</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>SD</mi><mrow><mi>g</mi><mn>2</mn><mi>v</mi><mo>,</mo><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub></mrow>]]></math><img file="FDA0000976029500000024.GIF" wi="470" he="150" /></maths><maths num="0004"><math><![CDATA[<mrow><msub><mi>SD</mi><mrow><mi>v</mi><mn>2</mn><mi>g</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>SD</mi><mrow><mi>v</mi><mn>2</mn><mi>g</mi><mo>,</mo><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub></mrow>]]></math><img file="FDA0000976029500000025.GIF" wi="462" he="148" /></maths>式中,F<sub>1</sub>为电动汽车用户成本,T=24h;SD<sub>g2v,t</sub>和SD<sub>v2g,t</sub>分别是t时段的充电量和放电量;I<sub>i,t</sub>为电池充放电状态,放电为‑1,充电为+1,不充电也不放电为0;D<sub>i,t</sub>(I<sub>i,t</sub>,I<sub>i,t‑1</sub>)表示t‑1时刻到t时刻第i辆电动汽车电池状态变化的损耗成本,D<sub>i,t</sub>(I<sub>i,t</sub>,I<sub>i,t‑1</sub>)的第一部分表示快速充放电对电池的损耗,第二部分表示充放电模式变换对电池的损耗,第三部分表示过度放电会损害电池,第四部分是一个常数;步骤5,根据初始化的群体P计算出系统经济调度总成本;经济调度模型中的总成本包含发电成本、备用成本、失负荷成本、污染排放成本和弃风成本,其中备用成本包含购买备用成本、实际调用备用成本和电动汽车充放电调度成本,系统经济调度总成本计算公式如下:<maths num="0005"><math><![CDATA[<mrow><mi>min</mi><mi> </mi><msub><mi>F</mi><mn>2</mn></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><mrow><mo>(</mo><msubsup><mi>C</mi><mi>t</mi><mi>g</mi></msubsup><mo>+</mo><msubsup><mi>C</mi><mi>t</mi><mi>r</mi></msubsup><mo>+</mo><msubsup><mi>C</mi><mi>t</mi><mi>int</mi></msubsup><mo>+</mo><msubsup><mi>C</mi><mi>t</mi><mrow><mi>p</mi><mi>l</mi></mrow></msubsup><mo>+</mo><msubsup><mi>C</mi><mi>t</mi><mi>w</mi></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000976029500000031.GIF" wi="1038" he="181" /></maths>发电成本<img file="FDA0000976029500000032.GIF" wi="99" he="70" /><img file="FDA0000976029500000033.GIF" wi="630" he="143" />备用成本<img file="FDA0000976029500000034.GIF" wi="86" he="63" /><img file="FDA0000976029500000035.GIF" wi="1013" he="92" />失负荷成本<img file="FDA0000976029500000036.GIF" wi="110" he="74" /><img file="FDA0000976029500000037.GIF" wi="398" he="83" />污染排放成本<img file="FDA0000976029500000038.GIF" wi="107" he="79" /><img file="FDA0000976029500000039.GIF" wi="750" he="138" /><maths num="0006"><math><![CDATA[<mrow><mi>e</mi><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>g</mi><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>&alpha;</mi><mi>i</mi></msub><msubsup><mi>P</mi><mrow><mi>g</mi><mi>i</mi><mo>,</mo><mi>t</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msub><mi>&beta;</mi><mi>i</mi></msub><msub><mi>P</mi><mrow><mi>g</mi><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>+</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>+</mo><msub><mi>&delta;</mi><mi>i</mi></msub><msup><mi>e</mi><mrow><msub><mi>&lambda;</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>P</mi><mrow><mi>g</mi><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub></mrow></msup></mrow>]]></math><img file="FDA00009760295000000310.GIF" wi="885" he="96" /></maths>弃风成本<img file="FDA00009760295000000311.GIF" wi="98" he="70" /><img file="FDA00009760295000000312.GIF" wi="917" he="207" />式中,F<sub>2</sub>为系统经济调度总成本,T=24h;D<sub>re,t</sub>为实际调用备用容量;E<sub>ens,t</sub>为失负荷量;其中,当P<sub>w,t</sub>≤P<sub>ew,t</sub>时,风机输出真实风能,没有弃风成本;当P<sub>w,t</sub>&gt;P<sub>ew,t</sub>时,风机输出计划风能,弃风成本不为零;用D<sub>h,t</sub>变量来描述备用调用关系,假设D<sub>h,t</sub>表达式如下:<maths num="0007"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>D</mi><mrow><mi>h</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>=</mo><msub><mi>D</mi><mrow><mi>l</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>+</mo><msub><mi>&Delta;P</mi><mrow><mi>l</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>+</mo><msub><mi>D</mi><mrow><mi>v</mi><mi>l</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>P</mi><mrow><mi>g</mi><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><msub><mi>d</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>-</mo><msub><mi>SD</mi><mrow><mi>v</mi><mn>2</mn><mi>g</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>-</mo><msub><mi>P</mi><mrow><mi>w</mi><mo>,</mo><mi>t</mi><mo>,</mo></mrow></msub><msub><mi>P</mi><mrow><mi>w</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>&le;</mo><msub><mi>P</mi><mrow><mi>e</mi><mi>w</mi><mo>,</mo><mi>t</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>D</mi><mrow><mi>h</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>=</mo><msub><mi>D</mi><mrow><mi>l</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>+</mo><msub><mi>&Delta;P</mi><mrow><mi>l</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>+</mo><msub><mi>D</mi><mrow><mi>v</mi><mi>l</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>P</mi><mrow><mi>g</mi><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><msub><mi>d</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>-</mo><msub><mi>SD</mi><mrow><mi>v</mi><mn>2</mn><mi>g</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>-</mo><msub><mi>P</mi><mrow><mi>e</mi><mi>w</mi><mo>,</mo><mi>t</mi><mo>,</mo></mrow></msub><msub><mi>P</mi><mrow><mi>w</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>&gt;</mo><msub><mi>P</mi><mrow><mi>e</mi><mi>w</mi><mo>,</mo><mi>t</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA00009760295000000313.GIF" wi="1467" he="330" /></maths>当D<sub>h,t</sub>&lt;0时,表示供大于需,不需要调用备用容量,也没有失负荷量,电网可以储存剩余的电量;D<sub>h,t</sub>&gt;0则调用车辆充放电量与备用容量,若全部调用备用容量与车辆可调度量还不满足负荷需求,将产生失负荷量;步骤6,评价群体P,对P进行非支配排序和拥挤距离计算;步骤7,对种群P进行选择、交叉和变异,生成满足系统各个约束条件的种群Q,评价Q;步骤8,生成新种群R=P∪Q,对R作非支配排序和拥挤距离计算;步骤9,从R中选择X个最优个体组成种群P;步骤10,判断是否满足终止条件,若不满足约束条件,则返回步骤7;若满足终止条件,则生成Pareto最优解集,并确定最终解,最终解包括各时刻各台发电机组的出力、计划风能输出、购买备用容量、电动汽车充放电控制矩阵,最后由最终解来调整各个时刻的系统运行;终止条件的判断是指判断多目标遗传算法的迭代次数是否已经达到最大遗传代数Y代,若还未达到,则返回步骤7;若已达到,则生成最优解集;最终解可以按系统运行的需要,从最优解集中选择其中的一组。
地址 310014 浙江省杭州市下城区潮王路18号