李奕杰, 宋恒, 葉晨暉, 申慶祥
(1. 國電南瑞科技股份有限公司, 江蘇 南京 211106;2. 國網江蘇省電力有限公司泰州供電分公司, 江蘇 泰州 225300)
近年來, 電動汽車應用作為能源產業(yè)結構調整和交通電氣化建設的新興解決方案, 得到了發(fā)展和普及[1]。 但隨著電動汽車用戶的激增, 電動汽車用戶的隨機駕駛行為和聚合充電行為不可避免地同時影響電網和交通網絡[2]。 電動汽車與交通系統(tǒng)和電力系統(tǒng)的相互作用成為新的挑戰(zhàn)和機遇[3-4]。
到目前為止, 電動汽車充電負荷預測從方法建模角度而言, 主要有基于模型驅動和基于數(shù)據(jù)驅動兩種類型的研究。 1) 基于模型驅動考慮交通路況, 從電動汽車與電網及交通之間能量和信息交互的角度, 分析掌握電動汽車行駛和充電模式的本質。 例如, 鎖軍[5]等提出一種考慮道路拓撲特性的電動汽車充電需求預測模型, 但未考慮環(huán)境溫度和速度因素對電動汽車電量消耗的影響。 張美霞[6]等提出了一種考慮環(huán)境溫度和車速的電動汽車充電需求預測模型, 采用起止點矩陣法分別獲得私家車和出租車出行的起訖點, 但影響起迄點(origin destination, OD) 分布的因素很多, 非路段交通量因素所能包容, 離實用階段尚有很大距離。2) 基于數(shù)據(jù)驅動意味著采用回歸分析、 神經網絡分析、 深度學習等人工智能算法, 實現(xiàn)基于歷史數(shù)據(jù)的短期預測或基于趨勢的長期預測, 例如, Liu Dunnan[7]提出一種基于反向傳播神經網絡的電動汽車負荷預測方法, 吳丹[8]等提出一種基于極端梯度提升與輕量級梯度提升機融合的電動汽車充電負荷預測模型。
模型驅動的電動汽車充電負荷預測方法需要假定大量模型參數(shù)以表達復雜的充電行為, 數(shù)據(jù)驅動的預測方法需要大量多源異構數(shù)據(jù)支撐訓練學習,二者均具有各自的局限性[9]。 為解決上述問題,本文提出融合模型驅動與數(shù)據(jù)驅動各自優(yōu)勢的充電需求預測模型, 充分考慮道路交通的拓撲特性和環(huán)境溫度等因素的同時, 采用雙向長短期記憶遞歸神經網絡算法對用戶出行時間和空間相對規(guī)律挖掘分析, 精準構建電動汽車充放電負荷預測模型。
電動汽車行駛軌跡數(shù)據(jù)是離散的、 不連續(xù)的和廣泛分布的。 因此, 交通規(guī)劃理論的地理網格建模方法是集中和處理這些軌跡數(shù)據(jù)的有效方法[10]。結合住宅區(qū)、 商業(yè)區(qū)、 工業(yè)區(qū)的分布, 建立如圖1所示的路網-配電網-車網交互模型。
圖1 配電網-交通網-車網交互模型
所建立的交通路網模型共有32 個交通節(jié)點、53 條道路, 進一步使用圖論的方法來進行相應的路網拓撲結構建模, 建立的模型如圖2 所示。
圖2 交通路網拓撲結構
式中,Cu為所選交通網格編號的集合;V為各交通路口的集合;E為各條道路的集合;K為劃分的各個時間段k集合;W為各條道路通行代價, 描述為道路路阻;G為交通路網。
本文研究內容針對城市交通路網, 如圖3所示。
圖3 城市道路路阻示意圖
在城市交通路網中, 電動汽車的出行代價(耗電量) 主要受到道路長度、 道路車流量和交通路口信號燈的影響, 因此引入時間-流量模型[11]對動態(tài)交通路網建模。 城市道路路阻可表示為:
式中,C(t,vij) 為節(jié)點路阻模型;R(t,vij) 為路段路阻模型。
1) 路段路阻模型
式中, 飽和度S是唯一變量, 代表道路的擁堵程度; 路阻影響因子α、β為道路規(guī)劃固定參數(shù)。
2) 節(jié)點路阻模型
式中,c為信號周期;λ為綠信比;q為路段車輛到達率。
電動汽車的行駛速度會直接影響單位里程的耗電量, 并進而影響其充電需求[12]。 因此, 構建一種基于路網實時車流量統(tǒng)計的速度-流量實用模型。 該模型可以描述t時刻車輛行駛于道路R(i,j) 時的車速vij(t)。
式中,vij,m(t) 代表道路零流速度;Cij表示最大通行能力;qij(t) 表示t時刻車流量,qij(t)/Cij為飽和度;β為經驗系數(shù);a、b、n為不同道路等級下的自適應系數(shù)。
電動汽車的單位里程耗電量受到諸多因素的影響, 包括行駛速度和環(huán)境溫度等, 環(huán)境溫度變化會導致車內空調設備附加電量消耗。 此外, 車速下降也會增加電動汽車的實時耗電量[13]。 因此構建了考慮環(huán)境溫度和速度的電動汽車實時電量消耗模型。
式中, 當環(huán)境溫度為Tp時, 車輛以速度vij行駛距離為S時空調所消耗的電量為KT;E為不同車速下的實時單位里程耗電量;FT為單位里程耗電量;Tmin為空調制熱溫度下限;Tmax為空調制冷溫度上限;WL和WR為空調的制冷和制熱功率。
電動汽車一般分為三種類型: 通勤私家車、 出租車、 公交車。 由于公交車有固定的發(fā)車班次, 行駛路線不會受到駕駛員主觀因素的影響而改變, 且充電時間較為固定, 所以本研究主要考慮通勤私家車、 出租車兩種類型。
為了獲得電動汽車的流動模式, 電動汽車充電需求負荷預測在線平臺需要提前預測從一個地區(qū)到另一個地區(qū)的汽車數(shù)量, 將此問題表述為OD 起訖矩陣預測問題。 而一般研究都假設出發(fā)地為住宅區(qū)或通過OD 矩陣反推得到電動汽車OD 起訖矩陣,但由于影響OD 分布的因素很多, 非路段交通量因素所能包容, 故離實用階段尚有很大距離, 并沒有很好地解決這一問題。 為此, 引入能夠考慮更多影響OD 分布因素的神經網絡模型預測電動汽車的起訖矩陣[14-15]。 通過考慮日類型與天氣等多種因素,使用雙向長短期記憶遞歸神經網絡算法對出行數(shù)據(jù)挖掘分析, 預測電動汽車的出行時間和空間相對規(guī)律, 圖4 為交通網格劃分與雙向長短期記憶遞歸神經網絡預測OD 矩陣結果的示意圖。
圖4 交通網格和OD 矩陣示例
以15 min 為一個預測尺度單位, 一天采取96個點進行預測, 始發(fā)地為C1、 目的地為C9的24 h電動汽車預測圖如圖5 所示。
圖5 C1 到C9 的雙向LSTM 預測結果
雙向LSTM 預測次日起訖矩陣結果表明, 居民區(qū)C1到工業(yè)區(qū)C9的預測集中分布在06:30—09:00的時間段。
私家車和出租車駕駛行為會受到路段長度、 行駛速度、 通行時間、 耗油量等權值的影響[16]。 本文基于路段道路出行代價, 使用動態(tài)Dijkstra 算法以最小出行代價為目標對電動汽車的駕駛路徑進行規(guī)劃。 圖6 為采用動態(tài)Dijkstra 算法規(guī)劃網格C1至網格C9的不同時間段的車輛行駛路徑。
圖6 網格C1 至網格C9 的車輛行駛路徑
EV 初始時刻荷電狀態(tài)(SOC) 服從正態(tài)分布[17]。 不同類型電動汽車電池容量服從伽馬分布[18], 二者結合即可得到初始時刻電量Q0(i)。
在該模型中, 當SOC 小于電池容量的0.2 時則視為電動汽車有充電需求, 該用戶執(zhí)行充電操作。 考慮到防止電池過充導致的電池損耗問題, 設置當充至電池容量的0.9 時充電結束, 則t時刻電量Qt(i) 和充電持續(xù)時間Tc為:
式中,η為能耗系數(shù);Ec表示每公里耗電量;Qt-1(i)為t-1 時刻的剩余電量; Δl為t-1 到t時刻第i輛車所行駛的距離;Pc為充電樁充電功率;ηc為充電樁充電效率, 取值為0.8~0.9。
設置快充、 慢充兩個功率, 通過分析充電用戶考慮每日里程、 出發(fā)/返回時間等不確定行為, 為電動汽車分配快速和慢速充電負載。 私家車主要用于工作通勤, 主要在居住地和工作單位往返一次,停車時間較長且停車地點固定, 充電方式主要為目的地慢充。 出租車車輛出行時間和起訖點隨機性較大, 行駛路線不固定且出行次數(shù)較多, 充電方式為前往充電站快充。 慢充功率設置為12 kW, 快充功率設置為45 kW。
若車輛到達各充電站過程服從泊松分布[19],以單位時間前往充電站充電的用戶數(shù)作為參數(shù)l,則用戶排隊等待時間Wq為:
式中,c為充電樁數(shù)量, 本文設置為10;μ為充電樁單位時間可完成充電的車輛數(shù), 設置為2;ρ為充電樁服務強度;Ls為排隊隊長。
在交通路網仿真模型與電動汽車駕駛和充電行為建模的基礎上, 開展電動汽車充電需求時空分布特征預測。 圖7 為根據(jù)一定空間尺度, 將所選的交通平面按住宅區(qū)、 商業(yè)區(qū)、 工業(yè)區(qū)劃分為不同的地理網格作為子區(qū)域。 仿真時間設置為24 h, 各節(jié)點配電網充電需求時空分布如圖8 所示。
圖7 測試區(qū)域路網拓撲示意圖
圖8 配電網充電需求時空分布
由圖8 可知, 該配電網-交通網-車網融合模型的總體充電需求時空分布特征, 用戶充電時間集中在09:00—13:00 和17:00—18:00 時段, 時間分布上呈現(xiàn)與充電需求相一致的“雙高峰” 類型,10、 15、 16、 23、 31 等節(jié)點的充電需求較大, 居民區(qū)節(jié)點的充電需求高于商業(yè)區(qū)和工業(yè)區(qū), 充電需求空間分布不均勻, 表現(xiàn)出較大充電需求的配電網節(jié)點主要集中在居民區(qū)和商業(yè)區(qū)附近。
圖9 為該模型24 h 不同時段下電動汽車充電總需求, 可知該模型下的充電最高峰大致分布在09:00—10:00。
圖9 24 h 該區(qū)域電動汽車充電需求
針對該充電需求時空分布特征預測模型, 通過對比分析不同配電網節(jié)點的充電需求, 可為各配電網節(jié)點制定電網功率調度策略提供重要信息, 并進一步分析充電負荷對電網的影響, 對電力系統(tǒng)負荷調度或充電站預設等研究具有重要意義。
為了更直觀地分析各功能區(qū)域的充電負荷規(guī)律, 圖10—11 給出了不同車輛類型即通勤私家車與出租車的負荷需求分布。
圖10 通勤私家車充電需求
由圖10 (a) 可知, 通勤私家車充電負荷高峰時段跨度較廣, 主要集中在15:00—19:30 時段,最高峰充電需求時段為18:00—19:00; 由圖11(a) 可知, 出租車的充電時間分布呈現(xiàn)與居民基礎用電負荷相仿的“雙高峰” 形式, 充電負荷時空分布特征符合居民出行規(guī)律。 試驗結果分析表明, 充電需求時空分布預測與實際情況相一致, 驗證了所提負荷預測方法的有效性。
圖11 出租車充電需求
最后, 為了應對電動汽車規(guī)?;l(fā)展的需求,圖12 和圖13 給出了24 h 不同交通節(jié)點的電壓變化情況和在19:37 時間斷面下不同EV 滲透率的節(jié)點電壓分布情況。
圖12 交通節(jié)點電壓
圖13 不同滲透率下交通節(jié)點電壓
由圖12 可知, 各節(jié)點電壓情況將受到電動汽車充電負荷大小的影響, 由圖13 可知, 隨著EV滲透率的增加, 節(jié)點電壓發(fā)生跌落。 當電動汽車充電規(guī)模進一步擴大時, 節(jié)點電壓隨著滲透率的增加可能跌落到節(jié)點安全電壓閾值以下。
本文基于神經網絡和動態(tài)交通信息, 充分考慮道路交通的拓撲特性模擬實際道路上的行駛行為,通過對出行數(shù)據(jù)挖掘分析電動汽車充放電的出行時間和空間相對規(guī)律, 精準構建電動汽車充放電負荷預測模型, 提出一種融合模型驅動和數(shù)據(jù)驅動的電動汽車充電負荷時空預測方法, 通過算例仿真, 得到如下結論。
1) 動態(tài)交通模型充分考慮道路交通的拓撲特性和環(huán)境溫度等因素, 引入路網實時車流量統(tǒng)計的速度-流量實用模型, 精準模擬城市路網特性。
2) 通過采用雙向長短期記憶遞歸神經網絡算法對用戶出行時間和空間相對規(guī)律挖掘分析, 使用動態(tài)Dijkstra 算法以最小出行代價為目標對電動汽車的駕駛路徑進行規(guī)劃, 實驗結果分析表明, 充電需求時空分布預測與實際情況相一致。
本文限于篇幅只研究了充電需求時空分布特征, 在接下來的研究中可以融合影響電動汽車充電行為的因素, 考慮影響駕駛用戶充電決策的心理因素, 搭建不同場景對充電負荷時空分布的影響, 并對試驗結果進一步分析, 為評估電動汽車充電負荷對電網影響和電動汽車充電控制策略與充電站擴容、 選址等研究提供參考。