面向高比例新能源并网场景的风光-电动车协同调度方法
A Wind-Solar-Electric Vehicles Coordination Scheduling Method for High Proportion New Energy Grid-Connected Scenarios
通讯作者: 李莉,女,副教授,电话(Tel.):010-61771855;E-mail:lili@ncepu.edu.cn.
责任编辑: 孙伟
收稿日期: 2022-02-22
基金资助: |
|
Received: 2022-02-22
作者简介 About authors
李林晏(1997-),男,山西省运城市人,硕士生,从事风光出力特性分析、风光互补系统优化调度研究.
风光-电动汽车协同调度能够有效降低风光出力和电动汽车无序充电等多重不确定性对电力系统的不利影响.现有优化调度模型多以等效负荷波动最小为优化目标,仅考虑了等效负荷的整体波动性,无法衡量风光出力与负荷的匹配度且并未考虑不同出力场景下风光出力的差异性.针对上述问题,提出一种面向高比例新能源并网场景的风光-电动车协同调度方法.构建基于蒙特卡罗模拟的电动汽车无序充电模型;基于风光出力预测数据,构建基于Gap statistic和K-means++算法的风光出力典型日划分模型;以等效负荷方差和负荷追踪系数最小为双优化目标,构建风光-电动汽车协同调度模型,并采用NSGA-II算法求解.结果表明:所提模型能够有效提升风光出力与负荷的匹配度,降低等效负荷波动性,从而缓解风光出力和电动汽车无序充电等多重不确定性对电力系统的不利影响.
关键词:
Wind-solar-electric vehicles coordinated optimization scheduling can effectively reduce the adverse effects of multiple uncertainties of wind-solar output and disorderly charging of electric vehicles on the power system. Most of the existing optimization scheduling models take the minimum equivalent load fluctuation as the optimization objective, which, only considering the overall fluctuation of equivalent load, cannot measure the matching degree of output-load, and do not consider the difference of output in different output scenarios. Therefore, a wind-solar-electric vehicles coordination scheduling method for high proportion new energy grid-connected scenarios is proposed. First, the disordered charging model of electric vehicles by Monte Carlo simulation is constructed. Then, a wind-solar output typical day classification model using Gap statistical and K-means++ is constructed based on the forecasting data of wind and solar power. Finally, taking the minimum equivalent load variance and load tracking coefficient as the double optimization objectives, a wind-solar-electric vehicles coordination optimization scheduling model is established, and the NSGA-II algorithm is used to solve it. The results demonstrate that the proposed model can effectively improve the matching degree of wind-solar output and load, and reduce the fluctuation of equivalent load, so as to reduce the adverse effects of multiple uncertainties of wind-solar output and disorderly charging of electric vehicles on the power system.
Keywords:
本文引用格式
李林晏, 韩爽, 乔延辉, 李莉, 刘永前, 阎洁, 刘海东.
LI Linyan, HAN Shuang, QIAO Yanhui, LI Li, LIU Yongqian, YAN Jie, LIU Haidong.
针对风光-电动汽车协同调度,国内外相关学者进行了相应研究.Song等[7]提出的充放电策略考虑了电动汽车与风力发电在平抑负荷波动方面的互补效应;侯建朝等[8]以等效负荷方差和车主支付费用最小为优化目标,构建风电与电动汽车协同调度模型;冉忠等[9]以电动汽车消纳后的风电弃风量最小为目标,建立电动汽车参与受阻风电消纳的源荷优化控制模型;葛少云等[10]通过遗传算法求解以峰谷差率最小为目标的优化模型;郭春林等[11]建立了电动汽车与可再生能源双层协调优化模型,上层模型以等效负荷波动最小和可再生能源渗透率最大为双优化目标,下层模型以实现电动汽车充放电负荷跟踪为目标,并采用遗传算法求解上下层的优化模型.同时,由于风电和光伏出力具有很强的随机不确定性,在不同典型出力场景下,风电和光伏出力特性及其之间的互补关系具有明显差异性,且对电力系统的影响规律不尽相同[12-13],在含大规模新能源的电力系统中,风光出力的随机不确定性对电力系统的不利影响也愈发明显.利用典型日聚类方式从大量风电/光伏历史数据中提取典型出力场景集,以典型出力场景反映风光出力特性,并用于不同出力场景下优化调度的差异化建模分析,是目前电力系统调度分析常用的方法[14-15].
综上所述,现阶段风光-电动汽车协同调度模型多以等效负荷波动最小为优化目标,仅考虑了等效负荷的整体波动性,无法衡量风光出力与负荷的匹配度,且优化调度模型在不同出力场景下的适用性和鲁棒性也并未展开分析.
针对上述问题,提出一种面向高比例新能源并网场景的风光-电动车协同调度方法.首先,通过分析电动汽车行驶规律,构建基于蒙特卡罗模拟的电动汽车无序充电模型;其次,基于风光出力预测数据,构建基于Gap statistic和K-means++算法的风光出力典型日划分模型;再次,以等效负荷方差和负荷追踪系数最小为双优化目标,构建风光-电动汽车协同调度模型,并采用NSGA-II算法进行模型求解;最后,基于中国西北某省风光电站预测出力和电网预测负荷数据,验证所提方法的适用性和鲁棒性.
1 基于蒙特卡罗模拟的电动汽车无序充电模型
1.1 电动汽车行驶规律
电动汽车行驶规律主要包括车辆出发前往工作地点时间、返回住所时间、路程耗时和日行驶里程.2017年美国交通部对全美家庭出行调查结果[16]反映了车辆一天中的行驶规律分布,考虑到中国与美国在工作时间上较为接近,均集中于9:00—17:00,可以假定两国电动汽车具有相同的行驶规律.因此,基于上述调查结果对电动汽车充电曲线进行建模分析,得到电动汽车出发前往工作地点时刻符合混合高斯分布,概率密度表达式为
式中:x为时刻;α1=0.8303,为出发前往工作地点时刻,属于第一个高斯模型概率;μ1=7.1944 和 μ2=13.1477 分别为第一个和第二个高斯模型均值;σ1=1.5080 和σ2=8.7192分别为第一个和第二个高斯模型方差.
工作日电动汽车离开工作地点返回住宅区时刻符合高斯分布,概率密度表达式为
式中:μT,W-H=17.1972,为离开工作地点时刻均值;σT,W-H=2.6197,为离开工作地点时刻方差.
每日出发前往工作地点与从工作地返回住宅耗时均符合对数正态分布,概率密度表达式为
式中:μT=-1.0526,为行程耗时期望均值;σT=0.7491,为行程耗时方差.
每日行驶里程符合对数正态分布,概率密度表达式为
式中:μM=2.8314,为日行驶里程期望均值;σM=1.0575,为日行驶里程方差.
1.2 假设条件
单台电动汽车充电需求主要受电动汽车种类、行驶规律和充电功率等3个方面影响.参考目前电动汽车发展情况,对家用电动汽车作出如下假设:
(1) 百千米耗电为16.7 kW·h.
(2) 电动汽车一天仅需充电一次,电池容量为70 kW·h.
(3) 在工作区,如果时间允许,电动汽车采用慢充方式充电,否则采用快充方式充电;在住宅区,电动汽车采用慢充方式充电;电动汽车充电过程为恒功率,快充、慢充充电功率分别为8 kW和60 kW.
(4) 每次充电均充至满电量.
(5) 电动汽车每日出行时刻、出行耗时、日行驶里程均为独立分布.
1.3 基于蒙特卡罗方法的电动汽车无序充电模拟
假设待研究区域的电动汽车总数为40万辆,在无任何经济措施诱导和限制时,工作日有40%的电动汽车选择白天在办公区域进行充电,60%的车辆选择夜晚回家充电,用户均选择到达目的地后立刻充电,且每辆电动汽车在工作日的通勤规律一致.基于上述电动汽车行驶规律概率分布和假设条件,采用蒙特卡罗模拟法[17]进行电动汽车无序充电模拟,单个用户具体模拟过程如下:① 抽取出发前往工作地点时刻T1;② 抽取前往工作地点耗时,计算到达工作地点时刻T2;③ 抽取从工作地点返家时刻T3;④ 抽取从工作地点返家耗时,计算到家时刻T4;⑤ 抽取总行驶里程,计算所需充电电量Q;⑥ 抽取电动汽车充电地点,判断电动汽车充电方式,生成充电功率曲线.
2 基于Gap statistic和K-means++算法的风光出力典型日划分模型
风电和光伏出力受风能、太阳辐射强度等多种因素影响,不同典型日出力场景下风光出力特性及其对电力系统的影响规律不尽相同;为分析不同典型日下,所提模型和优化调度方法的适用性和鲁棒性,基于风光出力预测数据,提出基于Gap statistic 和K-means++算法的风光出力典型日划分方法,并用于不同典型出力场景下风光-电动汽车协同调度的差异化建模分析.
2.1 Gap statistic算法
Gap statistic算法[18]是通过计算聚类有效指标,确定数据集最佳聚类数目的方法,其结果可以应用到任何一种聚类算法中.算法通过比较平均分布的参考数据集期望值
(1) 对已知观测数据集进行聚类,计算Wk:
式中:Dk为第k簇数据中组内点间的欧氏距离之和;nk为属于第k簇的数据的个数;xi和xj为数据观测点;Ck为第k簇数据集合;k=1, 2, ···, K.
(2) 生成B个参考数据集,并计算每个参考数据集的
式中:
(3) 计算
式中:sk和sdk为中间变量.
2.2 K-means++算法
K-means++聚类算法[19]是基于K-means的改进算法,能够有效克服K-means算法对孤立点较敏感,聚类结果易受初始聚类中心的影响等问题.K-means++算法基本思想为假设已选取k个初始聚类中心(1<k<K),则在选取第k+1个聚类中心时,距离当前k个聚类中心越远的点会有更高的概率被选为第k+1个聚类中心;随机选取第一个聚类中心(k=1). K-means++算法的具体步骤如下:
(1) 给定含有n个样本的数据集X=(x1, x2, ···, xn),从数据集中随机选取一个样本作为初始聚类中心.
(2) 计算每个样本与当前聚类中心之间的最短距离:
(3) 计算每个样本被选为下一个聚类中心的概率:
(4) 在[0, 1]区间生成一个随机数R,用R依次减去P1,P2,···,Pn,取差值首次小于或等于0时对应的对象为下一个聚类中心.
(5) 重复(2)~(4)直至选取K个聚类中心.
(6) 计算每个其他样本数据对象与各类聚类中心的距离,将其赋给最近的类.
(7) 计算每类数据对象的平均值,将其作为每一类新的聚类中心.
(8) 重复(6)~(7)直至聚类中心不再发生变化,从而得到最佳聚类结果.
2.3 风光出力典型日划分流程
基于Gap statistic和K-means++算法的风光出力典型日划分的具体步骤如下:
(1) 数据采集及归一化处理.基于风电场和光伏电站预测出力数据获取新能源预测出力数据,采用最大最小归一化方法将其映射至[0, 1]区间,并将其处理为nd×96维新能源标准出力矩阵(nd为采样天数).
(2) 新能源出力特征值提取.考虑到光伏仅在白天存在出力,且新能源出力对夜间负荷影响较大,将每天6:00—18:00划分为白天,其余时间段划分为夜间,选取白天和夜间新能源出力均值和出力标准差作为特征值,得到两组nd×2维特征值矩阵.
(3) 基于Gap statistic算法确定最佳聚类数.基于白天、夜间出力特征值矩阵,用Gap statistic 算法分别确定最佳聚类数.
(4) 基于K-means++算法进行聚类分析.根据最佳聚类数,以白天、夜间出力特征值矩阵为输入,采用K-means++聚类算法分别对白天、夜间新能源出力进行聚类分析.
(5) 白天和夜间聚类结果融合.采用时间对标法对白天和夜间聚类结果进行融合,获取不同出力场景下新能源出力.
(6) 不同出力场景典型日确定.由于融合后无法确定聚类中心位置,所以本文通过计算不同类别内各样本与其余样本的欧氏距离,以欧氏距离之和最小的样本序列为该类别的聚类中心,即该出力场景下的典型日.
3 风光-电动汽车协同调度模型
3.1 优化目标
等效负荷方差通过计算等效负荷在调度周期内的离散程度来衡量等效负荷整体波动性,负荷追踪系数通过比较新能源出力与负荷的时序波动特性来衡量新能源出力与负荷的匹配程度.考虑等效负荷整体波动和风光出力与负荷的匹配度,构建以等效负荷方差和负荷追踪系数最小为双优化目标的风光-电动汽车协同调度模型.
(1) 等效负荷方差(F1).等效负荷方差主要衡量等效负荷整体波动性:
式中:n1为一个调度周期内数据点个数;L(t)为t时刻负荷;P(t)为新能源在t时刻的实际出力;
(2) 负荷追踪系数.等效负荷方差仅能衡量等效负荷整体波动性,无法衡量风光出力与负荷的匹配度.为了使风光互补发电系统能够更好地跟随负荷变化,保障系统安全稳定运行,从源荷匹配度角度出发,采用能够衡量风光出力与负荷匹配关系的负荷追踪系数[20](F2)为优化目标:
式中:P'(t)和L'(t)分别为最大最小归一化处理后t时刻的新能源出力和负荷;ΔP'(t)和ΔL'(t)分别t时刻为新能源出力和负荷的波动量标幺值;α(t)为t时刻新能源出力与负荷变化量变化的一致性.负荷追踪系数越接近于0,表明新能源出力与负荷在调度周期内的变化特性越一致,新能源联合出力与负荷匹配性越好.
3.2 约束条件
系统运行需满足如下约束.
系统功率平衡约束:
式中:PG(t)为常规机组在t时刻出力.
新能源出力约束:
式中:
电动汽车电量约束:
式中:Q(t)为t时刻电动汽车电量;Qmin和Qmax分别为电动汽车电池储存容量的上下限;pev(t)为t时刻电动汽车充电功率;Δt为充电时长.
3.3 求解方法
所提风光-电动汽车协同调度模型为多目标优化问题,采用NSGA-II多目标优化算法进行模型求解.NSGA-II算法通过引入快速非支配排序技术,提高了算法的速度,具有结构简单、收敛速度快、鲁棒性能好的优点,算法具体实施步骤见文献[21].
在多目标优化问题中,一般不存在一组解使得所有目标函数值同时达到各自的最优值,只能求得非支配解集或 Pareto解集.因此,本文采用模糊隶属度函数[22]分别表示每个Pareto解中各目标函数对应的满意度,通过满意度比较找出最优折衷解.
4 算例分析
基于中国西北某省2019年风光电站预测出力和电网预测负荷数据进行算例分析,数据时间长度为1年,时间分辨率为15 min.风电、光伏以及电网负荷比约为11:9:13,风电场并网装机容量为 11160 MW,光伏电站并网装机容量为 9180 MW,电网最大负荷为 12900 MW,风光装机容量达全省总装机容量的40%,属于高比例新能源电网;假定该省共有40万辆电动汽车可供调度.
4.1 典型日划分
K-means++算法需要事先确定聚类数目,对此采用Gap statistic算法分别确定白天、夜间新能源出力最佳聚类数,聚类范围选定为1~10类.白天、夜间不同聚类数下Gap值如图1所示.根据Gap statistic原理,当聚类数分别为5和4时,白天和夜间的Gap值第一次满足条件,因此确定白天和夜间最佳聚类数分别为5类和4类.
图1
采用K-means++算法对新能源出力特征值进行聚类分析,聚类结果如图2所示.新能源白天出力可划分为5类波动情况,夜间出力可划分为4类波动情况.从均值看,白天对应的低出力小波动、低出力波动、中出力波动、高出力波动、高出力大波动分别分布在(0.05,0.22)、(0.21,0.33)、(0.32,0.43)、(0.44,0.58)、(0.56,0.75)区间;夜间对应的低出力小波动、低出力波动、中出力波动、高出力波动分别分布在(0,0.13)、(0.11,0.24)、(0.19,0.32)、(0.30,0.48)区间.从方差看,白天对应的低出力小波动、低出力波动、中出力波动、高出力波动、高出力大波动的方差逐渐增大,表明新能源出力波动越大;夜间对应的低出力小波动、低出力波动、中出力波动方差逐渐增大,高出力波动时新能源出力波动性有所降低,表明夜间风电高出力时全省风电平滑效应较好.
图2
采用时间对标法对白天和夜间聚类结果进行融合,获取不同典型日对应的新能源出力,如图3所示.图中,时间分辨率为15 min,全年共可划分为16类典型日,不同典型出力场景下新能源出力具有一定差异性,但同时具有一定相似性,且白天出力较大.将其按照出力大小进行划分,主要可划分为全天低出力波动(典型日1~4),主要包括白天、夜间低出力场景,出现率为43.73%;白天低出力波动(典型日5和6),主要包括白天低出力、夜间中出力场景,出现率为5.29%;夜间低出力波动(典型日7~10),主要包括夜间低出力和白天中、高出力场景,出现率为27.30%;中出力波动(典型日11~14),主要包括白天、夜间中出力场景,出现率为15.04%;全天高出力波动(典型日15和16),主要包括白天、夜间高出力波动场景,出现率为8.64%;该地区新能源出力主要为全天低出力类型.
图3
4.2 风光-电动汽车协同调度结果
为验证所提模型的适用性和鲁棒性,根据图3典型日划分结果,随机选取5种典型日出力场景,分别构建以等效负荷方差最小为单优化目标、以等效负荷方差和负荷追踪系数最小为双优化目标的风光-电动汽车协同调度模型,以等效负荷峰谷差(F3)为指标评价不同优化结果.单目标优化模型采用遗传算法进行求解,各项参数与NSGA-II 算法具有相同设置.5种典型日出力场景下优化调度结果如表1所示.不同典型日出力场景下,相比单目标优化调度模型,双目标优化调度模型的等效负荷峰谷差平均降低了19%.表明所提优化模型可以在降低等效负荷整体波动性的基础上,提升新能源出力与电网负荷匹配度,缓解风光出力和电动汽车无序充电等多重不确定性对电力系统的不利影响,且在不同典型出力场景下均具有良好的适用性和鲁棒性.
表1 不同典型日出力场景下优化调度结果
Tab.1
典型日 | 优化 目标 | 优化调度结果 | |||
---|---|---|---|---|---|
F2 | F1/MW2 | F3/MW | α/% | ||
4 | 单目标 | 0.10 | 3395218.02 | 16110.13 | 95.06 |
双目标 | 0.06 | 3099184.40 | 13797.68 | 95.04 | |
5 | 单目标 | 0.11 | 2012967.13 | 14637.51 | 95.19 |
双目标 | 0.07 | 1832434.89 | 14463.81 | 95.03 | |
9 | 单目标 | 0.08 | 6138278.51 | 24960.08 | 95.01 |
双目标 | 0.05 | 5561544.85 | 23093.18 | 95.00 | |
12 | 单目标 | 0.08 | 740009.18 | 13849.89 | 95.19 |
双目标 | 0.05 | 475042.78 | 11241.12 | 95.03 | |
15 | 单目标 | 0.09 | 1481212.29 | 20593.76 | 95.01 |
双目标 | 0.04 | 973148.80 | 9973.25 | 95.00 |
图4
对比不同典型日双目标优化结果可知,新能源出力限电行为均主要发生在新能源出力较大和等效负荷波动较大的时刻;优化后电动汽车充电曲线根据典型日出力特性不同,随新能源出力变化而变化以降低新能源出力对电网负荷的影响.
为进一步分析出力场景对优化调度结果的影响,分别从典型日5和15中任选一天进行风光-电动汽车协同调度,与典型日优化调度结果进行对比分析,优化结果如图5所示.
图5
综上所述,在不同出力场景下,电动汽车充电效果主要与此类典型日固有特征相关;所提优化调度模型在不同典型出力场景下均具有良好的适用性与鲁棒性.
5 结论
提出一种面向高比例新能源并网场景的风光-电动车协同调度方法,所得结论如下:
(1) 基于Gap statistic和K-means++算法的风光出力典型日划分模型能够实现白天和夜间新能源出力典型日的有效划分,新能源白天出力可划分为低出力小波动、低出力波动、中出力波动、高出力波动、高出力大波动共5类波动情况;夜间出力可划分为低出力小波动、低出力波动、中出力波动、高出力波动共4类波动情况,融合后共划分为16类典型出力场景.
(2) 以等效负荷方差和负荷追踪系数最小为双优化目标的源荷协同调度模型能够有效降低等效负荷波动性,提高新能源出力与电网负荷匹配度,从而减轻高比例可再生能源场景下风光出力和电动汽车无序充电不确定性对电力系统的不利影响.
参考文献
含高比例可再生能源电力系统灵活性规划及挑战
[J].
Power system flexibility planning and challenges considering high proportion of renewable energy
[J].
基于典型运行场景聚类的电力系统灵活性评估方法
[J].
Flexibility evaluation method for power system based on clustering of typical operating scenarios
[J].
基于分时电价的电动汽车有序充电控制策略设计
[J].
Coordinated charging strategy for electric vehicles based on time-of-use price
[J].
High efficient valley-filling strategy for centralized coordinated charging of large-scale electric vehicles
[J].DOI:10.1016/j.apenergy.2016.10.117 URL [本文引用: 1]
Reduced grid operating costs and renewable energy curtailment with electric vehicle charge management
[J].DOI:10.1016/j.enpol.2019.111051 URL [本文引用: 1]
Modeling electric vehicle charging station expansion with an integration of renewable energy and Vehicle-to-Grid sources
[J].DOI:10.1016/j.tre.2019.06.006 URL [本文引用: 1]
Integration of plug-in hybrid and electric vehicles: Experience from China
[C]//
计及需求响应的风电-电动汽车协同调度多目标优化模型
[J].
Multi-objective optimization model of collaborative WP-EV dispatch considering demand response
[J].
电动汽车参与受阻风电消纳的源荷优化控制方法
[J].
Source-load optimization control method for electric vehicles participating in blocked wind power consumption
[J].
电动汽车有序充电的峰谷电价时段优化
[J].
Optimization of peak-valley TOU power price time-period in ordered charging mode of electric vehicle
[J].
V2G模式下电动汽车与可再生能源双层协调优化
[J].
Bi-level coordinated optimization of electric vehicles and renewable energy sources in V2G mode
[J].
面向风电接纳能力评价的风资源/负荷典型场景集生成方法与应用
[J].
The generation method and application of wind resources/load typical scenario set for evaluation of wind power grid integration
[J].
生成风电功率时间序列场景的双向优化技术
[J].
A two-dimensional optimal technology for constructing wind power time series scenarios
[J].
充分发挥多能互补作用的风蓄水火协调短期优化调度方法
[J].
Short-term coordinated scheduling of wind-pumped-hydro-thermal power system with multi-energy complementarities
[J].
风光水互补发电系统与需求侧数据中心联动的优化调度方法
[J].
A optimal scheduling method for hybrid wind-solar-hydro power generation system with data center in demand side
[J].
National household travel survey (2017)
[DB/OL]. (
基于出行链随机模拟的电动汽车充电需求分析
[J].
Analysis on charging demand of EV based on stochastic simulation of trip chain
[J].
Estimating the number of clusters in a data set via the gap statistic
[J].DOI:10.1111/1467-9868.00293 URL [本文引用: 1]
K-Means++: The advantages of careful seeding
[C]//
风光水多能互补发电系统日内时间尺度运行特性分析
[J].
Analysis on intraday operation characteristics of hybrid wind-solar-hydro power generation system
[J].
An efficient constraint handling method for genetic algorithms
[J].DOI:10.1016/S0045-7825(99)00389-8 URL [本文引用: 1]
/
〈 |
|
〉 |
