朱永胜 ,王 杰 ,瞿博阳 ,李 健 ,SUGANTHAN P N
(1.郑州大学 电气工程学院,河南 郑州 450001;2.中原工学院 电子信息学院,河南 郑州 450007;3.南洋理工大学 电机与电子工程学院,新加坡 639798)
由于能源、环境的双重压力,电动汽车EV(Electric Vehicle)受到了各国政府、能源企业以及汽车厂商的广泛关注。研究表明,在中等的发展速度下,至2020年、2030年和2050年,美国汽车总量中电动汽车的比例将分别达到35%、51%和62%[1]。我国纯电动汽车以及插电式混合动力汽车的累计产销量也将在 2020 年突破 500 万辆[2]。
经济调度是电力系统优化运行的关键问题。未来如此大规模的电动汽车充电负荷,受用户出行需求、入网车辆数量、电池设备特性等多种因素的影响,在加重电网负担的同时,又具有典型的随机不确定性[3];特别是车载储能装置以“车-网”互动 V2G(Vehicle to Grid)模式[4]接入电网后,其复杂的充放电行为使调度问题成为一个不同时段间联系紧密的动态耦合系统,这将给传统的经济调度特别是动态经济调度DED(Dynamic Economic Dispatch)带来新的挑战[5],近年来已经成为学者们的研究热点[6-10]。
同时,随着国家对电力工业强制性减排要求的日益严苛以及全社会对环保问题的持续关注,环境因素在调度问题中的地位愈发重要,且这也更符合电动汽车“节能”、“环保”的自身定位。但如果再综合环境调度目标,则上述动态经济调度成为动态环境经济调度DEED(Dynamic Economic Emission Dispatch)。相比于传统的在单一调度时段内以经济性最优为唯一目标的经济调度,DEED是典型的多目标、多时段、高维度、强约束的非线性最优化问题,尤其是在时间和能量上综合考虑电动汽车的行驶及充放电需求后,调度问题的建模和求解将更为复杂。因此,目前对含电动汽车DEED问题的研究较少。
文献[11]构建了计及电动汽车规模化接入的多目标动态经济调度模型,并提出自适应多目标差分进化算法求解,但环境因素在模型中仅作为约束存在,且其只讨论了电动汽车充电的场景;文献[12-13]将燃料费用、排放成本等多目标问题转化为社会总成本单目标问题,并采用双种群带精英学习策略的微分进化算法等方法求解,但仅靠充放电时间很难全面描述电动汽车的V2G情况;文献[14-15]建立了综合考虑经济及环境因素的单目标调度模型,并采用PSO等算法获得最优的调度方案,然而其对电动汽车的建模也仅涉及注册V2G电动汽车数量与基本的电池容量约束,并未详细考虑车主的出行需求及电池的充放电特征;文献[16]建立了计及车辆行驶特征、电池充放电情况的多目标动态经济调度模型,并利用NSGA-Ⅱ算法求解,但环境因素仅作为模型众多调度目标之一,并未做深入讨论,且所建模型忽略了涉及电量非线性平衡问题的电网有功损耗,在增大误差的同时也降低了调度问题本身的难度。
基于以上分析,本文构建了计及电动汽车可调度V2G接入模式的电力多目标DEED模型,该模型能够同时兼顾系统的经济和环境效益,且在常规电力约束外,充分考虑了参与调度电动汽车的出行需求以及车载电池的充放电特性等因素;并在前期多目标进化算法的研究基础上,设计一种采用改进MOEA/D,即IMOEA/D(Improved Multi-Objective Evolutionary Algorithm based on Decomposition)的优化调度求解方法。仿真验证了所提调度模型的合理性以及求解算法的有效性。
(1)燃料费用。
系统的燃料费用目标函数采用下式表示[14]:
其中,T为调度时段数;N为常规机组的台数;Pi,t为常规机组 i在时段 t的有功出力;ai、bi、ci为机组 i的费用系数。
(2)污染排放。
系统的污染气体排放目标函数可表示为[15]:
其中,αi、βi、γi、ζi、φi为常规机组 i的排放系数。
(1)功率平衡约束。
系统的功率平衡以等式约束形式给出:
其中,PCh,t为电动汽车在时段 t的充电负荷;PDch,t为电动汽车在时段t的放电负荷;PD,t为时段t的系统负荷;PL,t为时段t的网损,本文采用B系数法求得,其计算式如式(4)所示[15]。
其中,Bij、Bi0和 B00为网损参数。
(2)电池剩余电量约束。
电动汽车储能电池在时段t的剩余电量St为:
其中,ηC、ηD分别为充、放电效率;Δt为调度时间间隔;STrip,t为电动汽车在时段t行驶过程中消耗的电量。
其中,ΔS为单位距离的平均耗电量;L为行驶里程。
为保证电池的寿命和运行安全,其剩余容量St要满足:
其中,Smax、Smin分别为电池电量的上、下限。
(3)电动汽车充、放电功率约束。
一般情况下,电动汽车充放电功率不能超过其额定充放电功率:
其中,PNCh、PNDch分别为电动汽车的额定充、放电功率,具体受电池特点及线路容量影响。
(4)车主出行需求约束。
电动汽车的基本功能是要满足车主的出行需求。假设电动汽车在一个调度周期内完成一次充放电循环,则应有:
(5)常规机组出力约束。
其中,Pi,max、Pi,min分别为机组 i出力的上、下限。
(6)常规机组爬坡约束。
其中,URi、DRi分别为常规机组i的升、降爬坡速率。
(7)旋转备用约束。
其中,SR,t为时段t系统的旋转备用容量需求。
针对 DEED问题,文献[17]提出了一种基于MOEA/D的调度求解新方法。该方法利用分解算法将多目标环境经济调度问题转化为一定数量的单目标优化子问题,并利用进化算法在指定的邻域内同时对各子问题进行优化计算,并且方法中加入了相关约束处理及进化控制策略,最终通过多次迭代实现对最优前沿的有效逼近。该方法计算速度快、收敛特性好,且能得到在目标空间分布均匀的Pareto最优解集。
(1)分解算法。
利用Tchebycheff分解算法获得一定数量单目标优化子问题的计算方法如下[17-18]:
其中,gte代表分解后的优化子问题;M为目标函数个数,对于 EED 问题,M=2;fm(x)为第 m 个调度目标函数;x 为决策变量;S 为可行解的区域为算法的参考点;λ=[λ1,…,λM]T为算法的权重向量。 对于 m=1,2,…,M,存在且
(2)进化算法。
针对种群中的第 i个个体 xi(i=1,2,…,Np,Np为种群规模),根据单纯形格子点设计法,计算相应的权重向量[17-18],并定义与其权重向量欧氏距离相近的 H 个个体为其进化邻域,记为 B(i)= {i1,…,iH}。随机从其邻域 B(i)内选择父代 r1、r2和 r3,且 r1≠r2≠r3≠i,然后进行差分进化计算产生其子代y。
其中,FS和CR为进化控制参数。
含电动汽车的DEED问题是在传统多目标动态负荷调度的基础上,综合考虑电动汽车的行驶及充放电等因素限制,其模型更复杂且变量及约束间的耦合度更高。为此,本文在上述方法的基础上进行适应性改进,以获得所建模型最优的调度方案。
(1)种群设置。
在本文的调度模型中,决策变量设置为各个调度时段常规机组的出力以及电动汽车的充放电功率,算法的种群x表示为:
个体xi为其中的一个调度方案:
每个个体共有(N+1)×T 维。 其中,Pev,t(t=1,2,…,T)为时段t的充放电功率:当电动汽车在该时段处于充电状态时,Pev,t=PCh,t;而当电动汽车在该时段处于放电状态时,Pev,t=PDch,t。
(2)约束处理。
MOEA/D设计之初是用来解决无约束优化问题的,而所建模型具有各种等式及不等式约束。因此,本文复杂约束的有效处理是解决该调度问题的关键。本文采用惩罚函数法,利用所求解问题的目标函数和约束条件,构造无约束的增广目标函数,把非线性的约束规划问题转化成无约束规划问题求解。
模型的等式约束包括系统功率平衡约束及车主出行需求约束,将其直接计入增广目标函数存在一定困难。为此,本文在算法中设计一种决策变量的两步制处理策略,首先对等式约束进行调整:
a.针对模型的出行需求等式约束,根据式(9)调整不同时段的V2G功率以满足电量平衡;
b.在获得电动汽车充放电功率的基础上,根据式(3)依次调整各个时段的常规机组出力,以保证系统的功率平衡。
以上分步处理过程是通过调用决策变量动态调整模块来实现的,其主要步骤如下。
a.针对不同的等式约束,计算其约束违反量θ,如果θ小于或等于事先设定的阈值ε,或调整次数达到最大调整次数K,则转步骤c;否则执行步骤b。
b.根据所调整决策变量的情况,将θ/n(n为时段数或者机组数)叠加到每个决策变量,并根据变量的出力上下限,进行越界处理。
c.如果是针对车主出行需求约束中的充放电功率,则结束调整;如果是系统功率平衡约束中的机组出力,则待所有时段都调整完成后,结束处理过程。
以上等式约束的具体处理流程见图1。
图1 两步制等式约束处理流程图Fig.1 Flowchart of two-step equality constraint processing
对于调度模型中的不等式约束,本文采用如下分类处理的方法。
a.将其中的电动汽车充放电功率约束、机组出力约束、爬坡约束,纳入到以上动态调整模块的越界处理中[17]。
b.在等式约束的动态调整过程中达到最大调整次数后,仍然可能得到不可行解,算法将这些解的等式约束违反量,以及电池剩余容量约束和旋转备用约束违反量,记为系统的总约束违反量V(x),其计算式如下:
模型的增广目标函数F(x)以基本的调度目标函数f(x)与总约束违反量罚函数之和的形式给出:
其中,s为罚系数。
最终,在对多目标函数分解、进化的基础上,利用基于罚函数的约束处理方法,实现迭代计算过程中对不可行解的有效处理。
针对含电动汽车的DEED问题,采用IMOEA/D的求解步骤如下。
(1)参数设置。
设置电动汽车、电力系统及IMOEA/D的具体参数。
(2) 初始化。
a.初始化算法的种群Np、权重向量λ,并计算每个个体对应的进化邻域;
b.针对每个个体xi,利用两步制算法对等式约束中的决策变量进行动态调整,并对不等式约束进行分类处理,计算包含总约束违反量V(xi)在内的目标函数值F(xi),确定IMOEA/D的参考点z=[z1,z2]T,其中
c.迭代次数gen置零。
(3)算法更新。
a.对每个个体xi,利用差分进化生成新的个体y,并对y进行前述的等式及不等式约束处理,计算新的目标函数 F(y)。 如果 zm>Fm(y),则更新参考点值 zm=Fm(y)。
b.对邻域 B(i)内的个体 xr,利用 Tchebycheff法进行分解计算。如果 gte(y|λr,z)≤gte(xr|λr,z),则更新最优解 xr=y,且 F(xr)=F(y)。
(4)终止迭代。
如果gen等于最大迭代次数,则终止迭代计算;否则,gen=gen+1,返回步骤(3)。
(5)输出结果。
输出目标值、最优解及Pareto最优前沿,并应用模糊集理论[11]确定最优折中解。
本文采用10机系统进行调度研究,调度周期24 h,机组参数及负荷数据见文献[19]。电网内共有50000辆电动汽车,其电池容量为24kW·h(以Nissan Leaf为例),每 100 km 耗电 15 kW·h。
假定电动汽车在每天早上离家时的荷电状态SOC(State Of Charge)为 100%,并且在 07:00以及17:00开始的1 h内行驶在上下班路上(共50 km),其余时间均可参与电网调度。调度周期内的最低SOC限制及额定充放电功率限制均设定为其额定值的20%,车载电池的充、放电效率为0.85,系统的旋转备用需求设为各时段负荷值的10%。
首先,为保证IMOEA/D求解结果的真实可信,同时利用NSGA-Ⅱ算法求解该调度模型,以进行对比。IMOEA/D的种群规模设为100,邻域H为20,迭代次数为 5000 次,ε、K 分别取 10-6、10,罚系数s设为100,其余参数与文献[17]相同。NSGA-Ⅱ算法的交叉、变异概率分别取0.9、0.2,交叉、变异算子的分布指数取20,且采用本文的个体及约束处理机制,种群规模与迭代次数也与本文算法保持一致。2种算法得到的Pareto最优前沿和最优解分别见图2和表1。
图2 Pareto最优前沿比较Fig.2 Comparison of Pareto optimal fronts
表1 极端解和最优折中解比较Table 1 Comparison of extreme solutions and best compromise solutions
由优化结果可见,针对此类多目标问题,并不存在使2个调度目标同时最优的绝对最优解,取而代之的是Pareto最优解。若运行人员仅考虑经济效益而选择经济最优方案时,对环保不利;反之若仅考虑污染气体排放量最小而选择环境最优方案时,则以增加燃料费用为代价。因此,在多目标问题特别是DEED的决策方面,要充分挖掘Pareto最优解集所蕴含的信息,并综合考虑各方面的因素,进行科学合理决策。
同时,由图2可见,由于该调度问题的复杂性,传统的NSGA-Ⅱ算法并不能获得完整的最优前沿,相比之下,本文IMOEA/D得到的最优前沿范围更广,分布更为均匀。且通过表1可看出,NSGA-Ⅱ算法找到的只是局部最优解,明显劣于IMOEA/D的极端解和最优折中解。所以,本文算法能为经济效益和环境保护目标间的相互妥协提供更为丰富的调度信息,充分显示了该算法的优越性。
为验证调度方案的正确性,表2给出IMOEA/D最优折中解的具体情况,图3则给出了不同方案下电动汽车的充放电功率和SOC(图中时段1表示01:00—02:00,其他依此类推)。
表2 IMOEA/D的最优折中解Table 2 Best compromise solutions of IMOEA/D
图3 不同解的电动汽车充放电情况对比Fig.3 Comparison of EV charge/discharge power among different solutions
由图3可以看出,极端解和最优折中解对应的电动汽车充放电规律比较相似,两者只是在具体功率上有所不同,这也就改变了常规机组之间的负荷分配,导致最终的燃料费用及污染气体排放量的差异较大。
综合图3(a)、(b),电动汽车在22:00至次日 06:00,基本处于充电状态,以保证日间的出行及调度需求,至07:00出发时SOC达到100%;07:00—08:00车辆行驶在路上,车载电池放电,SOC有所下降;08:00—15:00为负荷的高峰期,最高达到2150 MW,最低也有1776 MW,在此期间,电动汽车基本处于放电状态,以缓解常规火电机组压力,SOC持续下降;由于17:00—18:00车主的行驶需求,车辆在16:00进行充电,SOC有所上升;而20:00及21:00是夜间的负荷高峰,电动汽车继续放电直至SOC达到下限,然后在夜间负荷低谷时补充电量直至次日出行。
(1)不同场景的调度研究。
在以上工作的基础上,为深入验证所提调度模型及调度方法,分别对以下场景进行动态环境经济调度研究。
场景1:无电动汽车接入系统。
场景2:50000辆电动汽车接入系统,采用直接控制充电方式[3],对充电时间加以限定,即电动汽车在18:00到家后开始充电,至第2天07:00离家时结束。假定在此期间其所充电量满足均匀分布,50 000辆电动汽车共需电能375 MW·h。
场景3:50000辆电动汽车接入系统,采用可调度的V2G模式。
在上述3个场景中,场景1只进行基本的DEED计算,其余的2个场景均有电动汽车参与。所有场景都采用本文的IMOEA/D进行求解,且算法种群、迭代次数以及其他参数与3.2节保持一致。限于篇幅,在此只给出了不同场景下的极端解和折中解对应的负荷曲线,分别见表3和图4。
表3 不同场景的极端解Table 3 Extreme solutions of different scenarios
图4 不同场景的负荷曲线Fig.4 Load curves of different scenarios
对比以上解的情况和负荷曲线,若电动汽车采用场景2的常规充电方式接入电网,相比于无电动汽车接入的场景1,系统的最优燃料费用和最优污染排放分别增加了$28800和6 070 lb。可见,电动汽车虽然能减少车辆在行驶过程中的污染排放,但单一的充电模式,特别是对传统化石能源的过度依赖,反而使充电过程中的经济与环境问题更加突出。同时,此场景下负荷增长集中在18:00至次日07:00,虽然有一定的填谷作用,但与原有的夜间负荷高峰期叠加,增加了系统的容量需求且影响系统安全。
而如果电动汽车采用场景3的可调度V2G接入模式,其最优燃料费用和最优污染排放反而比场景1分别减少了$11300和8310 lb。所以,相比于场景2的单一充电模式,采用可调度的智能充放电方式,通过其有效的双向能量交换,在高峰重负荷时缓解了常规机组的压力,能有效减少总的燃料费用和污染气体排放量;且场景3的充电行为集中在22:00至次日06:00,放电行为集中在08:00—15:00,这明显改善了系统的负荷特性,其峰谷差率降为35.89%,负荷率升至84.85%。3个场景的具体负荷特性指标见表4。
(2)不同电动汽车规模的调度研究。
为研究不同电动汽车接入规模对系统调度结果影响,本文借鉴文献[20]对风电的处理思路,定义电动汽车的渗透率为:
表4 负荷特性指标对比Table 4 Comparison of load characteristic indexes
其中,PD.peak为系统的峰值负荷。不同渗透率下的极端解及其变化趋势见表5和图5、图6。
表5 不同渗透率的极端解Table 5 Extreme solutions of different penetrations
图5 最优燃料费用随不同渗透率的变化趋势Fig.5 Curve of best fuel cost vs.penetration
图6 最优污染排放随不同渗透率的变化趋势Fig.6 Curve of best pollution emission vs.penetration
由以上结果可以看出,随着电动汽车接入规模的不断扩大,其渗透率不断提高,相应的最优燃料费用及最优污染排放也随之不断减小。可见,电动汽车以V2G形式接入电网并参与电网调度,在初期确实能够缓解常规火电机组的经济与环境压力。
然而,在本模型接入60000辆电动汽车(渗透率13.40%)后,随着电动汽车规模的继续扩大,最优燃料费用及最优污染排放反而由最低点向上攀升。这说明,随着电动汽车渗透率逐渐达到拐点,其对经济和环境效益的改善效果不断减弱,取而代之的是大规模电动汽车充电需求的不断增加,系统为应对此额外的充电负荷需牺牲一定的燃料费用及污染气体排放量。所以,在考虑电动汽车接入水平及制定相关规划时,并非接入规模越大越好,而是要根据系统实际情况进行综合评价。
基于电动汽车规模化接入电网这一应用背景,针对电力系统多目标DEED问题,本文主要做了以下工作。
a.构建计及电动汽车可调度V2G接入模式的多目标DEED模型,并在模型中充分考虑了电动汽车能量属性与交通属性的耦合,以及系统的网络损耗等因素。该模型在兼顾经济和环境效益的同时,很好地体现了电动汽车接入对系统调度的影响。
b.设计一种基于改进MOEA/D的优化调度求解方法,提出针对系统多等式约束的决策变量两步制处理策略,并通过对违反量的适当惩罚,最终实现对模型复杂约束的有效处理。
c.以10机电力系统进行DEED仿真计算,并对不同的调度方案进行分析研究,结果验证了所提调度模型及方法的合理性及有效性。
参考文献:
[1]DUVALL M,KNIPPING E,ALEXANDER M,etal.Environmental assessment of plug-in hybrid electric vehicles.volume 1:nationwide greenhouse gas emissions[R].Palo Alto,USA:Electric Power Research Institute,2007.
[2]中华人民共和国国务院. 《节能与新能源汽车产业发展规划(2012—2020 年)》[EB/OL]. (2012-06-28)[2015-05-08].http:∥www.gov.cn /zwgk /2012-07 /09 /content_2179032.htm.
[3]王锡凡,邵成成,王秀丽,等.电动汽车充电负荷与调度控制策略综述[J]. 中国电机工程学报,2013,33(1):1-10.WANG Xifan,SHAO Chengcheng,WANG Xiuli,et al.Survey of electric vehicle charging load and dispatch control strategies[J].Proceedings of the CSEE,2013,33(1):1-10.
[4]KEMPTON W,TOMIC'J.Vehicle-to-grid power fundamentals:calculating capacity and net revenue[J].Journal of Power Sources,2005,144(1):268-279.
[5]郭建龙,文福拴.电动汽车充电对电力系统的影响及其对策[J]. 电力自动化设备,2015,35(6):1-9.GUO Jianlong,WEN Fushuan. Impact of electric vehicle chargingon powersystem and relevantcountermeasures[J].Electric Power Automation Equipment,2015,35(6):1-9.
[6]赵俊华,文福拴,薛禹胜,等.计及电动汽车和风电出力不确定性的随机经济调度[J]. 电力系统自动化,2010,34(20):22-29.ZHAO Junhua,WEN Fushuan,XUE Yusheng,et al.Power system stochastic economic dispatch considering uncertain outputs from plug-in electric vehicles and wind generators[J].Automation of Electric Power Systems,2010,34(20):22-29.
[7]HE Y,VENKATESH B,GUAN L.Optimal scheduling for charging and discharging of electric vehicles[J].IEEE Transactions on Smart Grid,2012,3(3):1095-1105.
[8]DERAKHSHANDEHSY,MASOUMAS,DEILAMIS,etal.Coordination of generation scheduling with PEVs charging in industrial microgrids[J].IEEE Transactions on Power Systems,2013,28(3):3451-3461.
[9]张新松,顾菊平,袁越,等.电动汽车最优充电模型及其近似求解方法[J]. 中国电机工程学报,2014,34(19):3148-3155.ZHANG Xinsong,GU Juping,YUAN Yue,et al.Optimum electric vehicle charging model and its approximate solving technique[J].Proceedings of the CSEE,2014,34(19):3148-3155.
[10]魏大钧,张承慧,孙波,等.基于分时电价的电动汽车充放电多目标优化调度[J]. 电网技术,2014,38(11):2972-2977.WEI Dajun,ZHANG Chenghui,SUN Bo,et al.A time-of-use price based multi-objective optimal dispatching forcharging and discharging of electric vehicles[J].Power System Technology,2014,38(11):2972-2977.
[11]邱威,张建华,刘念,等.计及插电式混合动力汽车规模化接入的多目标发电优化调度[J]. 电工技术学报,2013,28(12):160-167.QIU Wei,ZHANG Jianhua,LIU Nian,etal.Multi-objective optimal generation dispatching with large-scale plug-in hybrid electric vehicles penetration[J].Transactions of China Electrotechnical Society,2013,28(12):160-167.
[12]何明杰,彭春华,曹文辉,等.考虑电动汽车规模化入网的动态经济调度[J]. 电力自动化设备,2013,33(9):82-88.HE Mingjie,PENG Chunhua,CAO Wenhui,etal.Dynamic economic dispatch considering large-scale integration of electric vehicles[J].Electric PowerAutomation Equipment,2013,33(9):82-88.
[13]蔡秋娜,文福拴,薛禹胜,等.基于SCUC的可入网混合电动汽车优化调度方法[J]. 电力系统自动化,2012,36(1):38-46.CAIQiuna,WEN Fushuan,XUE Yusheng,etal.AnSCUC-based optimization approach for power system dispatching with plug-in hybrid electric vehicles [J].Automation of Electric Power Systems,2012,36(1):38-46.
[14]SABER A Y,VENAYAGAMOORTHY G K.Plug-in vehicles and renewable energy sources for cost and emission reductions[J].IEEE Transactions on Industrial Electronics,2011,58(4):1229-1238.
[15]GHOLAMIA,ANSARIJ,JAMEIM,etal.Environmental/economic dispatch incorporating renewable energy sources and plug-in vehicles[J].IET Generation,Transmission&Distribution,2014,8(12):2183-2198.
[16]李惠玲,白晓民,谭闻,等.基于智能电网的动态经济调度研究[J]. 电网技术,2013,37(6):1547-1554.LI Huiling,BAI Xiaomin,TAN Wen,et al. Research on dynamiceconomicdispatch based on smartgrid [J].Power System Technology,2013,37(6):1547-1554.
[17]朱永胜,王杰,瞿博阳,等.含风电场的多目标动态环境经济调度[J]. 电网技术,2015,39(5):1315-1322.ZHU Yongsheng,WANG Jie,QU Boyang,et al.Multi-objective dynamic economic emission dispatching of power grid containing wind farms[J].Power System Technology,2015,39(5):1315-1322.
[18]LI H,ZHANG Q.Multiobjective optimization problems with complicatedParetosets,MOEA /D andNSGA-Ⅱ[J].IEEE Transactions on Evolutionary Computation,2009,13(2):284-302.
[19]BASU M.Dynamic economic emission dispatch using nondominated sorting genetic algorithm-Ⅱ[J].International Journal of Electrical Power&Energy Systems,2008,30(2):140-149.
[20]王彩霞,乔颖,鲁宗相,等.低碳经济下风火互济系统日前发电计划模式分析[J].电力系统自动化,2011,35(22):111-117.WANG Caixia,QIAO Ying,LU Zongxiang,etal. Day-ahead dispatch mode for wind thermal power system in low-carbon economy[J].Automation of Electric Power Systems,2011,35(22):111-117.