吳洲洋，艾 欣，胡俊杰

（新能源電力系統(tǒng)國家重點實驗室，華北電力大學(xué)，北京市102206）

0 引言

隨著風(fēng)電、光伏等新能源的不斷接入，其短期出力難以預(yù)測的特性與負荷側(cè)的不確定性使電力系統(tǒng)面臨源-荷雙重不確定性的挑戰(zhàn)［1］，從而造成短期功率失衡，這對系統(tǒng)的充裕性提出了更高的要求。而電動汽車（electric vehicle，EV）、溫控負荷等靈活性需求側(cè)資源可在短期內(nèi)主動調(diào)節(jié)用電，經(jīng)聚合商統(tǒng)一調(diào)配，可成為潛力巨大的輔助服務(wù)提供方［2-3］，緩解供需失衡，具有良好的社會效益和經(jīng)濟效益。但需求側(cè)資源的調(diào)用需要在滿足用戶需求的前提下進行，如EV 的儲能特性使其可在短時間內(nèi)快速調(diào)整充電功率，但需要滿足車主充電需求［4］；而現(xiàn)代建筑物普遍使用的中央暖通空調(diào)（heating，ventilation，and air conditioning，HVAC）具備較大的熱慣性［5］，也能在短期內(nèi)調(diào)整壓縮機功率而不使溫度大幅波動［6］。

當(dāng)參與調(diào)頻輔助服務(wù)時，負責(zé)統(tǒng)一調(diào)配的聚合商將面臨調(diào)頻和用戶用電2 個方向的需求。文獻［7］以滿足EV 用戶充電需求為前提定義了實時充電可行域，研究了不同調(diào)控策略下EV 集群跟蹤系統(tǒng)調(diào)整信號的準(zhǔn)確度表現(xiàn)；文獻［8］提出的調(diào)控策略在HVAC 參與需求響應(yīng)的同時將室溫約束在了舒適范圍內(nèi)，均取得了良好的調(diào)度效果。

通常，輔助服務(wù)市場會要求參與方在日前上報其不同時段的備用容量，作為參與方在該時段響應(yīng)功率調(diào)整的上限［9］。聚合商在日前上報時，首先需要考慮備用容量與用電需求的耦合性［10］問題，將用電功率和備用容量聯(lián)合調(diào)度［11-12］。同時，為了滿足用戶與電網(wǎng)調(diào)頻的雙方需求，備用容量的優(yōu)化不僅要考慮資源的瞬時功率邊界［13-14］，也要考慮一個時段內(nèi)的電量儲備［15］。由于輔助服務(wù)市場中的調(diào)頻費用要基于跟蹤調(diào)整信號的準(zhǔn)確度結(jié)算，若電量儲備不足以持續(xù)提供調(diào)頻，則響應(yīng)準(zhǔn)確度會下降，從而降低收入［16-17］。文獻［18-19］均基于大量歷史數(shù)據(jù)進一步建立了基于備用容量來估計實時運行表現(xiàn)的模型。

在這方面，可以對不能滿足用戶需求和電網(wǎng)調(diào)頻的風(fēng)險做進一步量化分析。文獻［20］引入了條件風(fēng)險價值（conditional value at risk，CVaR）作為量化手段應(yīng)用于隨機優(yōu)化問題中。該指標(biāo)已被廣泛運用于電力市場研究［21-23］，將風(fēng)險成本量化并作為懲罰項引入優(yōu)化目標(biāo)中。

結(jié) 合 日 前、實 時2 個 時 間 尺 度，儲 能 系 統(tǒng)［24］或HVAC［25-26］均可以設(shè)計為日前分時段能量/備用調(diào)度和實時優(yōu)化控制的調(diào)度模型。然而目前的研究針對需求側(cè)資源缺乏統(tǒng)一的靈活性刻畫方法和調(diào)度策略，使得調(diào)頻服務(wù)的開展有2 點不足:①沒有將日前/實時2 個階段的調(diào)度問題結(jié)合起來，使得日前上報時的備用空間難以評估；②EV 和HVAC 是工業(yè)園區(qū)的常見柔性負荷，便于聚合和統(tǒng)一調(diào)配，但由于用電特性、功率級別的不同，總是各自單獨出現(xiàn)于調(diào)度模型的研究中，缺乏日前和實時結(jié)合的聯(lián)合調(diào)度方法。

因此，本文以含EV 和HVAC 的聚合商參與調(diào)頻輔助服務(wù)為例，提出了一種優(yōu)化調(diào)度模型，將EV和HVAC 聚合并統(tǒng)一調(diào)度，從而在日前階段進行電量、備用的聯(lián)合優(yōu)化，在實時階段進行功率的統(tǒng)一分配，在保證用戶需求的前提下使調(diào)頻收益最大化。

1 聚合商參與調(diào)頻備用的運行機制

本文參考美國PJM 電力市場的輔助服務(wù)交易機制，首先在日前對能量和調(diào)頻備用容量聯(lián)合出清；實時階段聚合商需響應(yīng)系統(tǒng)發(fā)送的調(diào)頻信號，備用結(jié)束后考核聚合商的調(diào)頻響應(yīng)準(zhǔn)確度和里程，并結(jié)算服務(wù)費用。因此在日前階段，聚合商通過評估自身用電靈活性進行備用容量上報（備用容量中的上、下容量相等，一同上報），每個備用的維持時段為1 h，在此期間，該聚合商可提供的向上、向下調(diào)整功率不超過該容量。在運行日的每個備用時段，系統(tǒng)運營商將根據(jù)電網(wǎng)的實時供需不平衡量以及交易完成并激活的備用資源計算發(fā)送給各備用資源的實時調(diào)整信號。所述日前/實時2 個階段的運行機制見圖1。

因此，為了完成輔助服務(wù)的優(yōu)化調(diào)度，聚合商需要處理2 個時間尺度的優(yōu)化問題:在日前備用優(yōu)化中，聚合商須確定每個時段的備用容量；在實時提供調(diào)整功率時，聚合商需實時分配調(diào)整功率給各用電單元，從而響應(yīng)調(diào)整信號。記備用時段i 內(nèi)的全部時刻t 構(gòu)成集合Γi。則從聚合商的視角，在某一時刻t，其需要響應(yīng)的調(diào)整功率PIRt為:

式中:Ri為時段i 內(nèi)聚合商的備用容量；δt為時刻t 系統(tǒng)發(fā)送的調(diào)整信號，取值范圍為[-1，1]。

最終，系統(tǒng)支付給聚合商的費用Qi將分為2 個部分（容量費用與里程費用），即

式中:QCi為時段i 內(nèi)的容量費用；QMi為時段i 內(nèi)的里程費用；πRi和πMLi分別為時段i 內(nèi)的容量價格和里程價格；MLi為時段i 內(nèi)的爬坡里程，由調(diào)整功率PIRt定義，如式（4）所示；ηi為時段i 內(nèi)的調(diào)頻準(zhǔn)確度，由式（5）定義。

因此，在日前調(diào)度階段，聚合商須首先確定分時段的用電功率PEi，在此基礎(chǔ)上估算其負荷的調(diào)整空間從而確定上報的備用容量Ri，本文將以這2 個變量作為優(yōu)化變量同時優(yōu)化。在實時運行中備用時段i 的每個時刻t，聚合商與電網(wǎng)交互的總功率Ptotalt可表示為用電功率和調(diào)整功率的總和，即

在任意時段i 內(nèi)，通過日前規(guī)劃購買的用電功率PEi保持不變，而聚合商響應(yīng)的調(diào)整功率PRRt則隨時刻t 而變化。

在本文場景中，聚合商能提供的靈活性資源由EV 集群和溫控負荷構(gòu)成，則PEi和PRRt均由這2 種靈活性資源的用電功率組成，即

式中:PEVt和PHVACt分別為在時刻t的EV 和HVAC 的用電功率，二者將作為實時功率分配中的優(yōu)化變量。

2 日前優(yōu)化模型——風(fēng)險厭惡的備用與用電功率日前調(diào)度

2.1 面向輔助服務(wù)的聚合商建模與電量可行域評估

針對上報容量時的風(fēng)險評估需要，本節(jié)將基于EV 和HVAC 的運行約束定義每個備用時段的總可調(diào)整電量。

1）EV 可消耗電量范圍

考慮EV 用戶的行為不確定性，對集群EV 的充電運行建模做如下假定:①聚合商與EV 用戶之間的合約對EV 離網(wǎng)時間不做限定；②為滿足用戶充電需求，設(shè)置每個時段的保底充電電量，避免EV 離網(wǎng)時充電不足，保證用戶滿意度。

記時刻t 在網(wǎng)充電EV 數(shù)量為Wt，則基于時段i內(nèi)在網(wǎng)EV 數(shù)的最大值ui和最小值vi進行電量空間核算，可表示為:

從備用時段內(nèi)聚合商交互的總電量視角，記保底充電電量為Sbase，根據(jù)電池的運行約束，EV 集群在一個時段內(nèi)可消耗的總電量應(yīng)既不使電池過充，同時也不低于保底功率下充入的電量。則EV 集群在時段i 內(nèi)可消耗電量的最大、最小值分別記為V，up和V，down。為提高魯棒性，此處假設(shè)最小值以使其最高的ui確定，最大值以使其最低的vi確定，即

2）HVAC 可消耗電量范圍

HVAC 可通過變頻系統(tǒng)調(diào)整制冷/制熱量實現(xiàn)用電功率的連續(xù)調(diào)節(jié)［6］。以夏天HVAC 工作于制冷模式為例，樓內(nèi)室溫隨制冷量變化的公式為［8］:

式中:Tin和Tout分別為室內(nèi)和室外溫度；α，β，γ 為常系數(shù)；PHVAC為HVAC 的功率。

假定2 個時刻之間的溫度變化為線性，對式（11）離散化后可得遞推公式，即

式中:下標(biāo)t 表示在時刻t 的對應(yīng)變量；ΔTint為在時刻t 的室內(nèi)溫度的變化量；Δt 為2 個時刻的間隔，根據(jù)系統(tǒng)發(fā)送調(diào)整信號的時間間隔而定。

考慮樓內(nèi)用戶對溫度的舒適性需求，記用戶可接受的最高、最低溫度分別為Tup和Tdown，且Tint∈[Tdown，Tup]。考慮恒溫運行下，使ΔTint=0，可得時段i 內(nèi)HVAC 用電功率的最大、最小值（PHVAC，upi和PHVAC，downi）分別為:

式中:Tout，ai為時段i 內(nèi)的平均室外溫度。

此時可推導(dǎo)得出HVAC 在時段i 內(nèi)消耗電量的最大、最小值（VAC，up和VAC，down）分別為:

式中:tspan為一個備用時段的持續(xù)時間。

3）聚合商可調(diào)整電量空間

綜上，在每個時段i，基于日前調(diào)度給定的用電計劃PEi，可得該聚合商的總可調(diào)整電量空間（上行、下行空間分別記為，up和，down），即

2.2 隨機變量的不確定性表達

式（3）中變量ηi和MLi需實時運行后結(jié)算時才能計算。本節(jié)將基于歷史數(shù)據(jù)及前文的容量核算，在缺少實時運行模型的前提下對ηi和MLi進行概率分布估計，從而進行日前調(diào)度中的風(fēng)險分析。

1）MLi的不確定性描述

聯(lián)立式（1）和式（4），將里程中的調(diào)整信號與備用容量解耦，可得:

2）ηi的不確定性描述

需強調(diào)，式（20）和式（22）的處理方式能夠成立的原因是本文參與的備用容量市場中上報的上下備用相等，均為Ri，因此才能將Ri提取出來進而解耦，當(dāng)備用市場可以將上下備用分別上報時，這種處理方法將不能成立。

3）各隨機變量的概率分布

以上推導(dǎo)中涉及的隨機變量的概率分布可由各自歷史數(shù)據(jù)分析得到。本文中將基于多維變量核密度 估 計（multi-variable kernel density estimation，mvKDE）理論給出概率分布。將每組隨機變量視為一個多維變量x（維數(shù)為d），mvKDE 方程的變量形式可由式（28）給出。

式中:yj，r為第j 個樣本的第r 維的值；xr為x 第r 維的值；m 為樣本總數(shù)；K(?)為圍繞著每一個樣本點的峰值函數(shù)；hr為第r 維的帶寬。

該方法可由MATLAB 內(nèi)置函數(shù)Ksdensity 實現(xiàn)。由該函數(shù)給出x 的概率分布可保存于概率向量Y（n 維）和場景矩陣X（n 行d 列）中，其數(shù)學(xué)意義為:x 取第k 行X 值的概率為Y 中的第k 個元素。由此可得出3 組隨機變量的3 個概率分布描述。

設(shè)3 種變量得出的場景數(shù)，即3 個場景矩陣的行數(shù)分別為n1，n2，n3。遍歷所有場景應(yīng)有n1n2n3個組合，記全部場景組合為n1n2n3行、7 列的矩陣Xi；將3 個概率向量的對應(yīng)元素相乘可得對應(yīng)的概率向量Yi。

2.3 優(yōu)化目標(biāo)及其不確定性建模

本節(jié)將基于以上聚合商不確定性建模，推導(dǎo)出基于CVaR 的風(fēng)險厭惡模型及求解方法。首先，日前優(yōu)化目標(biāo)應(yīng)使用能成本Ci最小，即購電成本減去式（3）定義的輔助服務(wù)收益，以PEi和Ri為優(yōu)化變量，模型如式（29）所示。

式中:πEi為聚合商的購電價格。變量ηi和MLi需實時運行后結(jié)算時計算，在日前優(yōu)化中視作隨機變量。

根據(jù)前文對相關(guān)隨機變量的概率描述，Xi中每個場景都可通過聯(lián)立式（5）、式（14）—式（19）、式（22）—式（25）求出準(zhǔn)確度ηi，通過聯(lián)立式（4）、式（20）和式（21）求出爬坡里程MLi，繼而得出對應(yīng)的Ci，記Ci可能發(fā)生的場景矩陣為Ci（維度為1），且其對應(yīng)概率向量也為Yi。

Ci的計算可由線性方程式（30）表示，其中Ci，k為Ci的第k 個元素；Xi，k為Xi的第k 行；A1，A2，A3，A4均為常數(shù)矩陣。

對于不確定場景，常規(guī)的處理方法是直接優(yōu)化Ci的 期 望 值，即∑Ci，kYi，k，Yi，k為 概 率 向 量Yi的 第k個元素。但在輔助服務(wù)結(jié)算的場景下，聚合商應(yīng)著重關(guān)注任何出現(xiàn)調(diào)頻誤差（即準(zhǔn)確度ηi＜1）的場景，因為系統(tǒng)會對各個輔助服務(wù)提供方的準(zhǔn)確度進行考核，甚至?xí)腥蹟鄼C制［27］，例如準(zhǔn)確度低于某個下限時需退出市場。

為應(yīng)用CVaR 理論，首先定義固定成本CG(PEi，Ri)，表示假定能夠完全響應(yīng)調(diào)整功率（無論上報多少容量，準(zhǔn)確度ηi都設(shè)為1）情況下的總成本，如式（31）所示。

CVaR 在定義上［28］，指的是給定一個成本閾值，對任一優(yōu)化變量PEi和Ri組合，總成本Ci超過該閾值條件下的期望值，而本文中Ci增加的原因來自于聚合商上報的備用容量超出自己的調(diào)整能力時，運行準(zhǔn)確度ηi會降低，從而可減少輔助服務(wù)費用。本文中取固定成本F 為這個成本閾值，則時段i 內(nèi)的CVaR（記為CVi），以及與不確定場景Xi對應(yīng)的場景矩陣CVi的元素CVi可以表示為:

式中:CVi，k為CVi的k 個元素。

至此將CVaR 作為風(fēng)險成本的表達，則優(yōu)化目標(biāo)應(yīng)為最小化固定成本和CVaR 之和。為消去式（33）中的max 運算，將CVi改為受約束變量的形式，最終得到取代式（29）的優(yōu)化模型，如式（34）所示，可調(diào)用線性優(yōu)化算法求解。

3 實時優(yōu)化模型——標(biāo)準(zhǔn)化裕度評估方法與實時功率分配

3.1 實時運行建模與約束

市場接受聚合商的PEi和Ri日前上報，則在次日的實時運行中，系統(tǒng)除提供給聚合商用電功率PEi之外，還會每隔數(shù)秒（美國PJM 的時間間隔為2 s）給聚合商發(fā)送功率調(diào)整信號PIRt。聚合商需要實時給各EV 及HVAC 分配用電功率及調(diào)整功率，以滿足各用戶的用能需求及系統(tǒng)的調(diào)頻需求。前文對聚合商中的HVAC 和EV 的可調(diào)整電量空間做了定義（式（18）和式（19）），在實時運行中同樣可以以該電量空間作為功率分配的基準(zhǔn)。

首先定義實時運行中的各用戶運行變量。根據(jù)前文對運行機制的介紹，在時刻t，需要進行優(yōu)化的變量包括HVAC 的功率PHVACt和編號為o 的EV 的功率P。

而此時聚合商可以掌握的狀態(tài)變量包括樓宇室內(nèi)溫度Tint和室外溫度Toutt、在網(wǎng)EV 總數(shù)Wt、編號為o 的EV 的 電 量S、本 時 段 的 初 始 電 量S，ini、本時段的已充電時間t。

為滿足用戶約束，應(yīng)使每輛EV 的總充電功率不小于保底功率且不造成過充，以及樓宇內(nèi)室溫處于舒適性區(qū)間，即

式中:μ 為EV 充電效率；PEV，max，PHVAC，max，PHVAC，min分別為實時運行中EV 的最大充電功率和HVAC 的最大、最小用電功率，其中EV 充電功率最小值約束為0，保底電量通過控制電量空間來保證。

式中:Pmax和Pbase分別為最大功率和保底功率

在實時調(diào)度中，聚合商將在確保用戶需求的前提下盡力滿足調(diào)整功率信號，若不能滿足則只能提供其最大調(diào)整功率?？傆秒姽β屎涂傉{(diào)整功率分別滿足式（37）和式（38）。

3.2 基于標(biāo)準(zhǔn)化裕度的實時功率分配

基于式（9）、式（10）、式（16）、式（17）對電量空間的定義，在實時視角下，各個用戶都有一個隨時刻t變動的電量空間，為對不同體量、不同形式的需求側(cè)資源的電量空間進行標(biāo)準(zhǔn)化量度，本文定義了實時運行下各用戶的標(biāo)準(zhǔn)化裕度參數(shù):在此時的用電功率下，將可調(diào)整電量消耗完畢所需的時間。本文中記EV 和HVAC 的實時標(biāo)準(zhǔn)化裕度分別為γ和VAC（以下簡稱裕度），基于時刻t 的優(yōu)化變量計算時刻t+1 的裕度的遞推公式如下。

圖2（a）和（b）分 別 為 實 時 運 行 下 某EV 和HVAC 以最大、最小功率運行時的裕度示意圖，灰色區(qū)域為狀態(tài)變量約束所定義的可行域。圖2（a）中，藍色曲線表示EV 電量，EV 部分的可行域以接入時電量和充電時刻對應(yīng)的點為起點A，灰色區(qū)域邊AB 的斜率為最大功率Pmax，邊AC 的斜率為保底功率Pbase，最大可能電量的直線構(gòu)成的三角形如綠色區(qū)域所示。保底功率Pbase=Sbase/tspan，表示恒功率充電時可充入保底電量的功率。圖2（b）中，藍色曲線表示室內(nèi)溫度，HVAC 最大可能溫度的直角構(gòu)成的三角形如綠色區(qū)域所示。

圖2 EV 和HVAC 的實時運行可行域與標(biāo)準(zhǔn)化裕度Fig.2 Real-time feasible operation region and standardized margin of EV and HVAC

基于本文所提裕度，該優(yōu)化調(diào)度方法的目的為:對可調(diào)電量空間較小的單元（需向上調(diào)整時的上行空間，或需向下調(diào)整時的下行空間），應(yīng)避免進一步消耗其該方向上的電量空間，從而避免出現(xiàn)功率分配不均衡，導(dǎo)致用戶需求不能滿足的情況。

4 模型仿真及結(jié)果分析

本章將以某聚合商參與日前能量市場與輔助服務(wù)市場為例，分析模型有效性并研究靈活性資源的特性。該聚合商可調(diào)度的靈活性資源包括所屬園區(qū)內(nèi)EV 充電站（80 個充電樁，對單個EV 的保底充電功率取3 kW），以及辦公樓宇內(nèi)的HVAC（制冷面積為1 500 m2，溫度舒適區(qū)間為24～28 ℃）。充電EV的概率分布以及模型中的輔助服務(wù)市場、能量市場的價格數(shù)據(jù)見附錄A。聚合商首先在日前市場購買用電功率并上報備用容量，然后在實時運行中響應(yīng)功率調(diào)整信號，以提高輔助服務(wù)效益、保證用電需求為目標(biāo)，實現(xiàn)對各EV 和HVAC 的功率優(yōu)化分配。模型仿真時長為24 h（備用時長為1 h，共劃分為24 個時段）。實時運行中的時間粒度為2 s，系統(tǒng)每隔2 s 發(fā)送一次調(diào)整信號，聚合商隨之調(diào)整總用電功率，并實時將功率分配給各EV 及HVAC。本模型中EV 充電數(shù)據(jù)來自中國某城市充電樁歷史合約，市場價格及調(diào)整信號數(shù)據(jù)來自美國PJM。

4.1 日前調(diào)度結(jié)果

本節(jié)除展示EV 集群和HVAC 構(gòu)成的聚合商的運行結(jié)果外，還將展示EV 集群獨立運行時的優(yōu)化調(diào)度結(jié)果以供對比分析。

表1 展示了EV 集群與HVAC 各自單獨運行時以及以聚合商的形式聯(lián)合運行時的電量成本、調(diào)頻收益與總成本。

表1 電量成本與調(diào)頻收益情況Table 1 Situation of electricity cost and frequency regulation profit

可以看出，EV 與HVAC 各自單獨運行時的總成本為4 262.4 元，大于合并運行時的3 372.1 元。顯然合并運行既降低了電量購買成本，也提高了調(diào)頻收益，從而大大降低了總的運行成本，體現(xiàn)了不同靈活性資源聚合后的相互支撐作用。下面將對調(diào)度結(jié)果進行展示并分析原因。

由于在原模型合并運行時，優(yōu)化調(diào)度只考慮總體的純用電功率和備用容量，而對EV 和HVAC 各自的貢獻不做區(qū)分，此處假定HVAC 的純用電功率應(yīng)為穩(wěn)定運行在最適宜溫度（(Tdown+Tup)/2）下的功率，HVAC 的備用容量則取穩(wěn)定運行的功率邊界（PHVAC，upi或PHVAC，downi）與純用電功率之差的絕對值。由于本文中對HVAC 的線性建模，取PHVAC，upi或PHVAC，downi的備用容量相同，非線性建模取最小值即可。其余用電功率和備用容量則視為EV 所屬。通過采用這種方法，以EV 集群為重點比較EV 單獨參與調(diào)頻時與EV 和HVAC 聯(lián)合參與調(diào)頻時的區(qū)別。

首先圖3（a）展示了僅EV 集群獨立運行時的充電功率購入和容量上報情況，可以看出上報容量基本可以與下行電量空間的期望值對應(yīng)，這說明對EV 來說，由于充電需求的存在，下行調(diào)整是制約調(diào)整容量的主要因素。從圖3（b）中可以看出，實際上報容量總是略高于下行電量空間的期望值，這是由于模型對CVaR 的評估認為通過多報一些容量產(chǎn)生的收益可以覆蓋風(fēng)險成本。

圖3 EV 集群獨立運行時的優(yōu)化結(jié)果Fig.3 Optimization results of independent operation for EV aggregator

圖4（a）和（b）將對合并運行時和EV 集群單獨運行時上報的備用容量與購買的用電功率進行對比。從圖4（a）可以看出，在合并運行時，EV 集群能提供的備用要大于獨立運行時的備用，這是因為HVAC的加入提升了聚合商整體的容量保底能力。具體地，在本文日前不確定性的考量中，HVAC 和EV 集群同時出現(xiàn)低容量的概率，會小于各自出現(xiàn)低容量的概率，從而降低了CVaR。同樣地，在圖4（b）中，EV 集群獨立運行時需購買的能量大于合并運行時供給EV 的能量，這是因為在備用能力提高時，EV 集群提供上行備用本身就可以充入電量，從而降低了購入電量的需求。由此驗證了不同靈活性負荷聯(lián)合運行下降低成本、提高收益的能力。

圖4 合并運行與EV 獨立運行時的調(diào)度結(jié)果對比Fig.4 Comparison of scheduling results between joint operation and independent operation of EV

4.2 實時功率分配結(jié)果

圖5（a）和（b）分別展示了第8 個小時內(nèi)EV 與HVAC 的運行結(jié)果。在這個時段內(nèi)準(zhǔn)確度為100%，備用容量為90 kW。圖5（a）中灰色曲線表示不同EV 的電池電量，起點和終點分別為接入和離網(wǎng)時間，藍色曲線為全部EV 的總功率?？梢钥吹剑趯υ６鹊募s束，所有EV 都能獲得充入電量的穩(wěn)定增長，同時還能響應(yīng)上調(diào)、下調(diào)功率。而圖5（b）中灰色區(qū)域為由式（14）和式（15）定義的、恒溫運行下HVAC 的功率可行域，該可行域表示了最魯棒的功率邊界，即在該范圍內(nèi)持續(xù)以最大或最小功率運行時，室內(nèi)溫度也不會超出用戶舒適范圍?？梢钥闯鲭m然HVAC 經(jīng)常在短時間內(nèi)以超出PHVAC，upi和PHVAC，downi的功率來響應(yīng)調(diào)整信號，但基于裕度評估的控制手段可以在隨后反向調(diào)整功率，從而將溫度控制在可接受范圍內(nèi)。同時注意到，以26 ℃為最理想溫度狀態(tài)，則提供調(diào)整備用后的樓宇溫度總是會稍偏高，可以推斷出相對于EV 而言，HVAC 的功率下調(diào)裕度比較大且容易恢復(fù)；而EV 由于保底充電功率的限制，下調(diào)裕度較小，但上調(diào)功率可以同時滿足充電需求。因此，HVAC 更傾向于響應(yīng)下調(diào)功率，而EV 更傾向于響應(yīng)上調(diào)功率。

圖5 實時功率分配結(jié)果Fig.5 Real-time power distribution results

此外，雖然上報備用只有90 kW，但EV 集群與HVAC 各自的功率波動范圍顯然遠高于這個范圍，這說明在聚合商基于各資源裕度的統(tǒng)一調(diào)配下，二者起到了很好的相互支撐作用。

圖6 給出了日前優(yōu)化時估計的準(zhǔn)確度期望值與實際運行時最終達成的準(zhǔn)確度對比?？紤]到調(diào)頻信號的隨機性，實際準(zhǔn)確度為同一天相同參數(shù)下運行100 次的平均值。由于式（9）、式（10）、式（16）、式（17）、式（24）中的假設(shè)都會低估聚合商的實際裕度從而低估期望準(zhǔn)確度，實際對比也體現(xiàn)出了這一點。本模型通過這種對電量空間的悲觀假設(shè)來降低實時運行中出現(xiàn)小概率風(fēng)險時的損失。

圖6 預(yù)期準(zhǔn)確度與實際準(zhǔn)確度對比Fig.6 Comparison of expected accuracy and actual accuracy

綜上所述，作為靈活性資源，EV 集群與HVAC的不同特性可以在聚合統(tǒng)一調(diào)度的運行模式下實現(xiàn)互補，EV 集群下行裕度的緊缺可由HVAC 補足，而EV 集群通過提供上行備用可以覆蓋用電成本。本文所提調(diào)度模型在日前優(yōu)化階段通過不確定性風(fēng)險分析實現(xiàn)了對聚合商備用能力的發(fā)掘，在實時調(diào)度階段則能夠使2 種資源相互支撐，既滿足用戶需求又能響應(yīng)調(diào)整功率。仿真結(jié)果表明，本文基于風(fēng)險分析的日前調(diào)度方法能夠?qū)崟r階段的不確定性進行預(yù)期，而EV 和HVAC 的互補特性能夠提高聚合商參與輔助服務(wù)的魯棒性。

5 結(jié)語

本文建立了需求側(cè)靈活性資源參與輔助服務(wù)的調(diào)度模型。通過基于通用電量空間的評估方法，在日前階段實現(xiàn)了能量、備用聯(lián)合優(yōu)化的風(fēng)險分析，在實時階段實現(xiàn)了聚合商內(nèi)各用戶的功率優(yōu)化分配，可得出以下結(jié)論。

1）基于CVaR 的風(fēng)險評估方法，可以描述不確定性場景下優(yōu)化目標(biāo)的預(yù)期損失，且便于應(yīng)用在線性規(guī)劃模型中。

2）EV 和HVAC 在參與輔助服務(wù)時相對于對方各有劣勢:EV 備用容量受限于其下行空間，HVAC不能通過提供備用覆蓋其用電成本。

3）基于本文所提標(biāo)準(zhǔn)化裕度方法進行聯(lián)合優(yōu)化與實時分配，可以使EV 和HVAC 發(fā)揮各自優(yōu)勢實現(xiàn)相互支撐，提高了運行效率和經(jīng)濟效益。

然而，本文僅考慮了美國PJM 市場模式下上下備用容量對稱的情況，在美國CAISO 等非對稱容量市場中，基于電量空間評估和CVaR 的優(yōu)化模型的線性化處理將不再適用，需要研究新的線性化處理或求解方法。此外，競價規(guī)則下聚合商報價曲線的優(yōu)化與博弈也將是下一步的主要研究工作。

本文受到北京市科技新星計劃資助項目（Z201100006820106）資助，特此感謝！

附錄見本刊網(wǎng)絡(luò)版（http：//www.aeps-info.com/aeps/ch/index.aspx），掃英文摘要后二維碼可以閱讀網(wǎng)絡(luò)全文。