• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    計(jì)及多種柔性負(fù)荷的虛擬電廠熱電聯(lián)合魯棒優(yōu)化調(diào)度模型

    2021-07-09 09:25:38王志南張宇華黃珂王叢
    電力建設(shè) 2021年7期
    關(guān)鍵詞:魯棒熱電時(shí)段

    王志南,張宇華,黃珂,王叢

    (上海電力大學(xué)電氣工程學(xué)院,上海市 200090)

    0 引 言

    以虛擬電廠(virtual power plant,VPP)形式聚合用戶側(cè)分散的、零碎化的分布式能源和柔性負(fù)荷(flexible load,F(xiàn)L),形成大功率、大容量穩(wěn)定可調(diào)控資源池[1-2],不受地域局限性,打破了多種能源分開(kāi)規(guī)劃、獨(dú)立運(yùn)行的既有模式,是促進(jìn)不同能源之間互補(bǔ)、提高可再生能源消納、加強(qiáng)柔性負(fù)荷管理與優(yōu)化的有效方式[3-5]。

    柔性負(fù)荷參與能源調(diào)度運(yùn)行可以有效提高系統(tǒng)調(diào)峰靈活性、平滑負(fù)荷曲線以及減少新裝機(jī)容量,是VPP在用戶側(cè)的重要可調(diào)度柔性資源。目前對(duì)負(fù)荷側(cè)柔性負(fù)荷的優(yōu)化研究大多只局限于電功率優(yōu)化,對(duì)于熱/電聯(lián)合系統(tǒng)下的多種不同特性柔性負(fù)荷協(xié)調(diào)優(yōu)化運(yùn)行考慮不足[6-9]。文獻(xiàn)[6]在考慮用戶用電經(jīng)濟(jì)性和電網(wǎng)供需關(guān)系下,構(gòu)建了一種以典型冰蓄冷空調(diào)為柔性對(duì)象的虛擬電廠雙層優(yōu)化模型。文獻(xiàn)[7]從柔性負(fù)荷上調(diào)/下調(diào)備用容量角度出發(fā),來(lái)分析柔性負(fù)荷的靈活調(diào)節(jié)作用。文獻(xiàn)[8-9]在計(jì)及負(fù)荷削減特性基礎(chǔ)上,又分別考慮了激勵(lì)負(fù)荷和轉(zhuǎn)移特性負(fù)荷。然而上述文獻(xiàn)中大多只側(cè)重于電力網(wǎng)中柔性負(fù)荷優(yōu)化,缺乏對(duì)熱柔性負(fù)荷的考慮。

    在熱電聯(lián)合運(yùn)行的基礎(chǔ)上,文獻(xiàn)[10]以舒適度模糊性為指標(biāo),將熱負(fù)荷轉(zhuǎn)化為具有彈性可調(diào)的區(qū)間變化,使用自回歸平均模型表示其動(dòng)態(tài)耦合關(guān)系,分析熱量在時(shí)間軸上的可調(diào)性。文獻(xiàn)[11]從供需角度考慮了冷熱電聯(lián)供系統(tǒng)中可平移負(fù)荷的運(yùn)行特性。然而,上述文獻(xiàn)對(duì)柔性負(fù)荷考慮形式較為單一,未涉及熱電系統(tǒng)下多種柔性負(fù)荷之間的協(xié)調(diào)優(yōu)化運(yùn)行,且較少考慮熱電聯(lián)合運(yùn)行中多種不確定因素的影響。

    VPP熱電聯(lián)合運(yùn)行中要面臨風(fēng)電、光伏和柔性負(fù)荷等多種不確定因素,給源、荷雙側(cè)協(xié)調(diào)運(yùn)行造成一定難度。在大多實(shí)際情況中,隨機(jī)變量確切的概率密度分布和隸屬度關(guān)系難以準(zhǔn)確獲得,基于誤差區(qū)間分析的魯棒優(yōu)化能夠較好地解決預(yù)測(cè)誤差對(duì)系統(tǒng)造成的不穩(wěn)定威脅,被更多運(yùn)用到不確定性處理中[12-14]。文獻(xiàn)[6-11]采用確定性模型或基于概率分析的多場(chǎng)景方法來(lái)處理系統(tǒng)中的不確定性,未能提高系統(tǒng)運(yùn)行魯棒性。文獻(xiàn)[12-13]采用魯棒優(yōu)化處理風(fēng)、光出力的不確定性,通過(guò)考慮不確定因素的邊界情況來(lái)提高系統(tǒng)決策的穩(wěn)定性。文獻(xiàn)[14]通過(guò)調(diào)節(jié)魯棒系數(shù)來(lái)平衡VPP運(yùn)行穩(wěn)定性與經(jīng)濟(jì)性的問(wèn)題。然而,現(xiàn)有研究中較少涉及到熱電聯(lián)合運(yùn)行系統(tǒng)下多種柔性負(fù)荷的不確定問(wèn)題分析,且未考慮不同魯棒系數(shù)下多種柔性負(fù)荷協(xié)調(diào)優(yōu)化問(wèn)題。

    基于以上分析,本文首先從用戶側(cè)自主響應(yīng)角度構(gòu)建可削減、可轉(zhuǎn)移、可平移柔性負(fù)荷的數(shù)學(xué)模型,分別以熱/電負(fù)荷曲線標(biāo)準(zhǔn)差值和系統(tǒng)綜合運(yùn)行成本最低為目標(biāo),建立考慮柔性負(fù)荷的VPP熱電聯(lián)合雙層協(xié)調(diào)優(yōu)化運(yùn)行模型。上層通過(guò)優(yōu)化負(fù)荷曲線,來(lái)解決VPP熱電聯(lián)合系統(tǒng)下多種柔性負(fù)荷間的協(xié)調(diào)優(yōu)化問(wèn)題;下層在上層的基礎(chǔ)上,來(lái)實(shí)現(xiàn)VPP熱電聯(lián)合下各個(gè)設(shè)備最優(yōu)出力,從而提高經(jīng)濟(jì)效益。針對(duì)系統(tǒng)中的不確定性,采用魯棒優(yōu)化方式處理,以置信區(qū)間描述風(fēng)、光出力和柔性負(fù)荷的預(yù)測(cè)誤差波動(dòng)。通過(guò)調(diào)節(jié)魯棒系數(shù),分析不同風(fēng)險(xiǎn)水平下多種柔性負(fù)荷協(xié)調(diào)優(yōu)化的效果,提高熱電聯(lián)合系統(tǒng)運(yùn)行的魯棒性。

    1 含柔性負(fù)荷的VPP熱電聯(lián)合運(yùn)行框架

    1.1 系統(tǒng)運(yùn)行框架與優(yōu)化思路

    本文在傳統(tǒng)熱電聯(lián)供的電力流、熱力流和天然氣流三者協(xié)調(diào)運(yùn)行基礎(chǔ)上,建立了VPP熱電聯(lián)合雙層協(xié)調(diào)優(yōu)化模型,上層以多種柔性負(fù)荷之間的協(xié)調(diào)優(yōu)化,來(lái)平滑負(fù)荷曲線;下層通過(guò)合理安排熱電聯(lián)產(chǎn)機(jī)組(combined heat and power,CHP)、燃?xì)忮仩t(gas boiler,GB)、電儲(chǔ)能(battery storage,BS)、熱儲(chǔ)能(heat storage,HS)等機(jī)組的出力來(lái)促進(jìn)不同能源間協(xié)調(diào)互補(bǔ),實(shí)現(xiàn)系統(tǒng)經(jīng)濟(jì)運(yùn)行??紤]到熱功率傳輸具有一定延遲性和用戶對(duì)熱負(fù)荷溫度感知的模糊性,故本文只考慮熱負(fù)荷中具有可削減、可平移特性的柔性負(fù)荷。系統(tǒng)運(yùn)行結(jié)構(gòu)如圖1所示。

    圖1 系統(tǒng)運(yùn)行結(jié)構(gòu)Fig.1 System operation structure

    在熱/電柔性負(fù)荷協(xié)調(diào)優(yōu)化基礎(chǔ)上,充分考慮系統(tǒng)內(nèi)多種不確定因素,采用魯棒優(yōu)化方法來(lái)處理系統(tǒng)內(nèi)風(fēng)、光出力和柔性負(fù)荷的不確定性,以置信區(qū)間來(lái)描述其預(yù)測(cè)誤差波動(dòng)。在熱/電聯(lián)合運(yùn)行基礎(chǔ)上,通過(guò)多種柔性負(fù)荷協(xié)調(diào)運(yùn)行來(lái)優(yōu)化熱/電負(fù)荷曲線,以系統(tǒng)經(jīng)濟(jì)最優(yōu)為目標(biāo),來(lái)求解VPP熱電聯(lián)合系統(tǒng)的日前優(yōu)化方案,具體流程如圖2所示。

    圖2 系統(tǒng)優(yōu)化流程Fig.2 System optimization process

    1.2 用戶側(cè)柔性負(fù)荷特性分析

    為充分挖掘負(fù)荷側(cè)柔性資源的調(diào)度潛力,首先建立多種柔性負(fù)荷數(shù)學(xué)模型,包括可平移、可轉(zhuǎn)移、可削減柔性負(fù)荷。

    1.2.1 可平移負(fù)荷

    可平移負(fù)荷(shiftable loads,SL)以固定工作時(shí)長(zhǎng)按調(diào)度計(jì)劃在時(shí)間軸上進(jìn)行多時(shí)段整體平移,如洗衣機(jī)、熱水器等。

    以24 h為1個(gè)調(diào)度周期T,未響應(yīng)前的可平移負(fù)荷功率分布向量Lshift為:

    Lshift=(0,…,Pts,Pts+1,…,Pts+td-1,…,0)

    (1)

    式中:Pts為起始時(shí)刻的可平移負(fù)荷;ts、td分別為可平移負(fù)荷的起始、持續(xù)時(shí)間。

    以[tsh-,tsh+]為可接受平移區(qū)間,則調(diào)度后的可平移負(fù)荷起始時(shí)段τ為:

    τ∈[tsh-,tsh+-td+1]

    (2)

    為保證其用電的連續(xù)性,應(yīng)滿足:

    (3)

    式中:αt為判斷負(fù)荷是否發(fā)生平移的0-1狀態(tài)變量,αt=1表示負(fù)荷發(fā)生了平移。

    調(diào)度后的可平移負(fù)荷功率分布向量Pshift為:

    Pshift=(0,…,Pτ,Pτ+1,…,Pτ+td-1,…,0)

    (4)

    響應(yīng)調(diào)度后用戶相應(yīng)的補(bǔ)償費(fèi)用Cshift為:

    (5)

    1.2.2 可轉(zhuǎn)移負(fù)荷

    與可平移負(fù)荷相比,可轉(zhuǎn)移負(fù)荷(translational load,TL)的工作時(shí)段、工作時(shí)長(zhǎng)和轉(zhuǎn)移量大小皆可調(diào)節(jié),沒(méi)有時(shí)間連續(xù)性約束,轉(zhuǎn)移后分布時(shí)段較為分散,靈活度更高。調(diào)度前后的可轉(zhuǎn)移負(fù)荷總功率須在可轉(zhuǎn)移區(qū)間[ttr-,ttr+]上滿足一致不變性:

    (6)

    式中:Ltrans,t、Ptrans,t分別為t時(shí)刻調(diào)度前、后的可轉(zhuǎn)移負(fù)荷功率。

    除了對(duì)轉(zhuǎn)移功率進(jìn)行約束外,還需要對(duì)可轉(zhuǎn)移負(fù)荷的最小持續(xù)時(shí)間進(jìn)行約束:

    (7)

    (8)

    響應(yīng)調(diào)度后用戶相應(yīng)的補(bǔ)償費(fèi)用Ctrans為:

    (9)

    1.2.3 可削減負(fù)荷

    可削減負(fù)荷(reducible load,RL)在響應(yīng)調(diào)度指令后對(duì)用戶部分用電量進(jìn)行削減,t時(shí)刻削減后的負(fù)荷功率Pcut,t為:

    Pcut,t=(1-utγt)Lcut,t

    (10)

    式中:Lcut,t為t時(shí)刻調(diào)度前可削減負(fù)荷功率;布爾變量γt表示t時(shí)刻是否發(fā)生削減的狀態(tài)變量,γt∈{0, 1};ut為t時(shí)刻的負(fù)荷削減系數(shù)。

    為了保證用戶用電的舒適性,不能頻繁地對(duì)用戶負(fù)荷進(jìn)行調(diào)度,需要對(duì)削減次數(shù)以及連續(xù)削減時(shí)間進(jìn)行約束:

    (11)

    (12)

    (13)

    響應(yīng)調(diào)度后用戶獲得相應(yīng)的補(bǔ)償費(fèi)用Ccut為:

    (14)

    2 系統(tǒng)不確定性處理

    針對(duì)VPP熱電聯(lián)合運(yùn)行系統(tǒng)中風(fēng)、光出力和柔性負(fù)荷的不確定性,采用魯棒優(yōu)化來(lái)進(jìn)行處理。系統(tǒng)內(nèi)不確定變量為:

    xi,t=[Pw,t,Ppv,t,Pshift,t,Ptrans,t,Pcut,t,Qshift,t,Qcut,t]

    (15)

    式中:Pw,t、Ppv,t分別為風(fēng)電、光伏的輸出功率;Qshift,t、Qcut,t分別為可平移、可削減熱柔性負(fù)荷功率。

    以置信區(qū)間來(lái)描述其預(yù)測(cè)誤差波動(dòng)范圍,考慮形式如下:

    (16)

    (17)

    (18)

    采用魯棒優(yōu)化處理要考慮其“最?lèi)毫印睏l件下發(fā)生的情況,故由式(16)所做出的決策會(huì)變得過(guò)于保守,而喪失經(jīng)濟(jì)性。為避免上述結(jié)果的產(chǎn)生,在此基礎(chǔ)上,引入魯棒控制系數(shù)Γ來(lái)調(diào)節(jié),?!蔥0, |J|],J為魯棒優(yōu)化所處理不確定參數(shù)的集合,此時(shí),不確定變量xi,t取值區(qū)間集合與魯棒系數(shù)的關(guān)系為:

    (19)

    式中:S(Γ)為隨機(jī)變量的取值范圍與魯棒系數(shù)的集合;變量βi由魯棒系數(shù)Γ決定。

    3 VPP熱電聯(lián)合雙層協(xié)調(diào)運(yùn)行模型

    在熱電聯(lián)合運(yùn)行下,充分考慮不同特性柔性負(fù)荷之間的協(xié)調(diào)互補(bǔ)性,且計(jì)及系統(tǒng)內(nèi)多種不確定性因素影響,建立VPP熱電聯(lián)合雙層協(xié)調(diào)運(yùn)行模型。上層以熱/電負(fù)荷曲線標(biāo)準(zhǔn)差最小為目標(biāo)函數(shù),來(lái)優(yōu)化負(fù)荷曲線,解決VPP熱電聯(lián)合系統(tǒng)下多種柔性負(fù)荷間的協(xié)調(diào)優(yōu)化問(wèn)題。下層在上層的基礎(chǔ)上,以系統(tǒng)綜合成本最低為目標(biāo)函數(shù),解決VPP熱電聯(lián)合下各個(gè)設(shè)備最優(yōu)出力問(wèn)題,提高經(jīng)濟(jì)效益。

    3.1 上層協(xié)調(diào)優(yōu)化模型

    3.1.1 目標(biāo)函數(shù)

    上層以系統(tǒng)中熱/電負(fù)荷曲線標(biāo)準(zhǔn)差D最小為目標(biāo),即:

    (20)

    1)電負(fù)荷曲線標(biāo)準(zhǔn)差。

    (21)

    Xt=[Pb,t,Pshift,t,Ptrans,t,Pcut,t]

    (22)

    式中:Xt為t時(shí)刻的電負(fù)荷矩陣;E為全為1的列向量;Pb,t為t時(shí)刻基礎(chǔ)電負(fù)荷功率。

    2)熱負(fù)荷曲線標(biāo)準(zhǔn)差。

    (23)

    Yt=[Qb,t,Qshift,t,Qcut,t]

    (24)

    式中:Yt為t時(shí)刻的熱負(fù)荷矩陣;Qb,t為t時(shí)刻基礎(chǔ)熱負(fù)荷功率。

    3.1.2 約束條件

    熱/電柔性負(fù)荷在滿足式(1)—(14)約束外,還要建立其魯棒優(yōu)化模型下的不確定約束,即:

    (25)

    (26)

    (27)

    (28)

    (29)

    (30)

    (31)

    (32)

    (33)

    (34)

    δshift,t≥1

    (35)

    δtrans,t≥1

    (36)

    δcut,t≥1

    (37)

    δq_shift,t≥1

    (38)

    δq_cut,t≥1

    (39)

    3.2 下層聯(lián)合運(yùn)行模型

    3.2.1 目標(biāo)函數(shù)

    下層以熱/電聯(lián)合系統(tǒng)中能耗、補(bǔ)償、運(yùn)維、環(huán)境以及交互的綜合成本F最低為目標(biāo),即

    (40)

    1)能耗成本。

    (41)

    式中:ηCHP,i、ηGB分別為CHP和GB的效率系數(shù);RN、HN分別為天然氣價(jià)格和熱值。

    2)運(yùn)維成本。

    (42)

    式中:μchp、μes、μqs、μgb、μw、μpv分別為熱電聯(lián)產(chǎn)、蓄電池、儲(chǔ)熱槽、鍋爐、風(fēng)電和光伏機(jī)組的運(yùn)維系數(shù)。

    3)環(huán)境成本。

    (43)

    式中:Vek、Dek分別為CHP和GB第k種污染物的排放量和單位成本系數(shù);Vk為第k種污染物的懲罰成本。

    4)與電網(wǎng)交互成本。

    (44)

    式中:Cprice,t為系統(tǒng)在t時(shí)刻與電網(wǎng)交互電價(jià)。

    3.2.2 約束條件

    1)熱、電能量平衡約束。

    (45)

    (46)

    式中:ηe,i、ηl,i分別為第i個(gè)CHP單元的發(fā)電效率和損失率。

    2)風(fēng)、光不確定性約束。

    (47)

    (48)

    (49)

    (50)

    δw,t≥1

    (51)

    δpv,t≥1

    (52)

    3)儲(chǔ)能設(shè)備約束。

    (53)

    考慮到儲(chǔ)能設(shè)備自身特性以及周期運(yùn)行約束,須滿足如下條件:

    (54)

    (55)

    (56)

    式中:L、I分別為儲(chǔ)能設(shè)備充、放電次數(shù)最大值。

    4)CHP、GB和主電網(wǎng)聯(lián)絡(luò)線功率約束。

    (57)

    5)爬坡率約束。

    -rm,down≤Pm,t-Pm,t-1≤rm,up

    (58)

    式中:rm,down、rm,up分別為設(shè)備m(m=CHP,GB)向下、向上的爬坡率。

    3.3 效果評(píng)價(jià)指標(biāo)

    本文采用平滑性指標(biāo)和經(jīng)濟(jì)性指標(biāo)來(lái)衡量多種柔性負(fù)荷協(xié)調(diào)優(yōu)化后的效果。

    1)平滑性指標(biāo)f1。

    負(fù)荷曲線標(biāo)準(zhǔn)差計(jì)算如式(21)—(24)所示,并以?xún)?yōu)化前后兩者差值作為衡量曲線平滑性指標(biāo)。

    f1=D0-D

    (59)

    式中:D0、D分別為柔性負(fù)荷參與前、后負(fù)荷曲線的標(biāo)準(zhǔn)差。

    2)經(jīng)濟(jì)性指標(biāo)f2。

    以?xún)?yōu)化前后系統(tǒng)綜合成本的差值來(lái)衡量,差值越來(lái)越大,經(jīng)濟(jì)性越優(yōu)。

    f2=F0-F

    (60)

    式中:F0、F分別為柔性負(fù)荷參與前、后的系統(tǒng)綜合成本。

    3.4 問(wèn)題求解

    上述所建立的計(jì)及柔性負(fù)荷的VPP熱電聯(lián)合魯棒優(yōu)化模型為0-1混合整數(shù)線性規(guī)劃(mixed integer linear programming, MILP)模型。本文采用GAMS軟件CPLEX求解器對(duì)其模型進(jìn)行尋優(yōu),求解優(yōu)化流程如圖3所示。

    圖3 系統(tǒng)求解流程Fig.3 System solution flow

    4 算例分析

    4.1 參數(shù)設(shè)置

    本文以北方某園區(qū)為研究對(duì)象,單位調(diào)度時(shí)長(zhǎng)為1 h。其中相關(guān)機(jī)組的污染物排放量、環(huán)境懲罰參數(shù)取自文獻(xiàn)[15];分時(shí)電價(jià)取自文獻(xiàn)[16];GB、CHP機(jī)組參數(shù)見(jiàn)表1;3類(lèi)柔性負(fù)荷的相關(guān)參數(shù)如表2—4所示。可削減負(fù)荷的削減系數(shù)為0.7,天然氣價(jià)格為2.64元/m3,單位熱值為9.7 kW·h/m3。

    表1 可控機(jī)組參數(shù)Table 1 Parameters of controllable units

    表2 可平移負(fù)荷參數(shù)Table 2 Parameters of shiftable loads

    表3 可轉(zhuǎn)移負(fù)荷參數(shù)Table 3 Parameters of transferable loads

    表4 可削減負(fù)荷參數(shù)Table 4 Parameters of reducible loads

    4.2 不同場(chǎng)景下柔性負(fù)荷優(yōu)化效果分析

    為了分析多種柔性負(fù)荷協(xié)調(diào)優(yōu)化對(duì)系統(tǒng)的影響,本文設(shè)置3種不同場(chǎng)景來(lái)進(jìn)行對(duì)比分析,來(lái)驗(yàn)證本文所提方法的優(yōu)越性。場(chǎng)景1不考慮柔性負(fù)荷參與;場(chǎng)景2考慮部分柔性負(fù)荷參與;場(chǎng)景3綜合考慮電/熱多種柔性負(fù)荷參與。具體場(chǎng)景設(shè)置如表5所示。

    表5 場(chǎng)景設(shè)置Table 5 Setting of scenes

    多種柔性負(fù)荷的協(xié)調(diào)優(yōu)化運(yùn)行不僅可以減少凈成本、提高經(jīng)濟(jì)性,還能對(duì)負(fù)荷起到一定的調(diào)節(jié)作用。本文將經(jīng)濟(jì)性和平滑性作為評(píng)估優(yōu)化結(jié)果的2個(gè)指標(biāo),不同場(chǎng)景下電/熱柔性負(fù)荷協(xié)調(diào)優(yōu)化效果對(duì)比如圖4所示。

    圖4 不同場(chǎng)景下優(yōu)化效果對(duì)比Fig.4 Comparison of optimization effect in different scenarios

    可以看出,在綜合考慮多種柔性負(fù)荷協(xié)調(diào)優(yōu)化后,在供電高峰時(shí)段的部分負(fù)荷功率發(fā)生了時(shí)空上的轉(zhuǎn)移、平移和削減,負(fù)荷曲線的平滑性有了明顯的改善,緩解高峰供電壓力的同時(shí)也提高了系統(tǒng)運(yùn)行的穩(wěn)定性。同時(shí),為了保證客觀性與實(shí)用性,多種柔性負(fù)荷優(yōu)化后并不改變?cè)胸?fù)荷的峰谷趨勢(shì)走向,可接受平移區(qū)間和轉(zhuǎn)移區(qū)間設(shè)置與現(xiàn)實(shí)情況相符,避免了柔性負(fù)荷優(yōu)化后出現(xiàn)在響應(yīng)能力較弱的00:00—04:00時(shí)段。其效果評(píng)價(jià)指標(biāo)如表6所示。由于場(chǎng)景1未考慮柔性負(fù)荷參與,故優(yōu)化前后調(diào)度策略未發(fā)生變化,差值為0;而場(chǎng)景2和3,電/熱曲線平滑性都所有上升,經(jīng)濟(jì)性也更優(yōu)。在綜合考慮電/熱柔性負(fù)荷協(xié)調(diào)優(yōu)化后,電負(fù)荷曲線標(biāo)準(zhǔn)差從146.70 kW減少到104.08 kW,平滑性提高了29.05%,熱負(fù)荷曲線標(biāo)準(zhǔn)差從64.766 kW減少到32.427 kW,曲線平滑性提高49.93%,同時(shí),節(jié)約成本為728.17元,系統(tǒng)經(jīng)濟(jì)性提高了7.5%。

    表6 不同場(chǎng)景下優(yōu)化效果指標(biāo)對(duì)比Table 6 Comparison of performance indices of optimization under different scenarios

    以場(chǎng)景3中電負(fù)荷為例,對(duì)多種柔性負(fù)荷協(xié)調(diào)優(yōu)化情況進(jìn)行詳細(xì)分析,優(yōu)化前、后柔性負(fù)荷分布情況如圖5所示。綜合考慮柔性負(fù)荷后,可平移電負(fù)荷1和可平移電負(fù)荷2在保證用電時(shí)間連續(xù)性和可接受平移區(qū)間下,分別從晚高峰時(shí)段20:00—22:00和午高峰時(shí)段12:00—13:00平移到了06:00—08:00和08:00—09:00的用電低谷時(shí)段,以緩解電網(wǎng)高峰的供電壓力。可轉(zhuǎn)移負(fù)荷在滿足最小持續(xù)轉(zhuǎn)移時(shí)間和可接受轉(zhuǎn)移區(qū)間下,從原來(lái)次高峰時(shí)段13:00—16:00轉(zhuǎn)移到電價(jià)平時(shí)和谷時(shí)多個(gè)時(shí)段,從原來(lái)4個(gè)時(shí)刻拆分成了6個(gè)時(shí)刻,其中大部分分布在05:00—07:00的低谷期,小部分分布在15:00—17:00的電價(jià)平時(shí)段,在削峰填谷的同時(shí),也使得負(fù)荷曲線更加平滑??上鳒p負(fù)荷分布在整個(gè)調(diào)度周期,多個(gè)時(shí)刻皆有響應(yīng),主要集中在11:00—14:00和18:00—21:00的2個(gè)高峰時(shí)段,并滿足最大、最小持續(xù)削減時(shí)間??梢钥闯觯喾N柔性負(fù)荷協(xié)調(diào)優(yōu)化后,負(fù)荷曲線的平滑性有了明顯的改善,緩解了高峰供電壓力。

    圖5 柔性負(fù)荷優(yōu)化前后效果對(duì)比Fig.5 Comparison of effect before and after optimization of flexible load

    4.3 優(yōu)化策略分析

    不同魯棒系數(shù)下,柔性負(fù)荷協(xié)調(diào)優(yōu)化效果如圖6所示。魯棒系數(shù)Γ取值不同代表著VPP決策者對(duì)柔性負(fù)荷不確定因素所帶來(lái)的風(fēng)險(xiǎn)意識(shí)不同,可將魯棒系數(shù)水平劃分為4部分,分別代表決策者為激進(jìn)型、較激進(jìn)型、較保守型和保守型,不同魯棒系數(shù)下柔性負(fù)荷優(yōu)化效果如圖6所示。當(dāng)Γ值較小時(shí),決策者更偏向于風(fēng)險(xiǎn)選擇,隨著Γ增大,基于決策的魯棒性,系統(tǒng)中可調(diào)柔性負(fù)荷減少,偏向于不確定因素的風(fēng)險(xiǎn)規(guī)避,系統(tǒng)穩(wěn)定性增強(qiáng)。對(duì)圖6(a)而言,在11:00—14:00和18:00—21:00時(shí)段,表現(xiàn)為可削減電負(fù)荷功率隨魯棒系數(shù)增大而不斷減少,而在05:00—11:00和14:00—17:00時(shí)段表現(xiàn)為優(yōu)化后的可轉(zhuǎn)移負(fù)荷和可平移負(fù)荷的功率變化軌跡,皆為減少趨勢(shì)。與電負(fù)荷相比,熱負(fù)荷曲線較為平緩,可削減熱負(fù)荷隨魯棒控制系數(shù)的變化趨勢(shì)集中體現(xiàn)在11:00—14:00和17:00—20:00的負(fù)荷高峰時(shí)段,而優(yōu)化后可平移熱負(fù)荷變化軌跡集中體現(xiàn)在06:00—10:00時(shí)段,隨著Γ增大有不同程度的減少。

    圖6 不同魯棒系數(shù)下柔性負(fù)荷優(yōu)化效果Fig.6 Optimal configuration effect of flexible load under different robustness coefficients

    不同魯棒系數(shù)下,柔性負(fù)荷協(xié)調(diào)優(yōu)化效果評(píng)價(jià)指標(biāo)如表7所示。可以看出,優(yōu)化后電/熱負(fù)荷曲線的平滑性和系統(tǒng)經(jīng)濟(jì)性都隨著可調(diào)控柔性負(fù)荷含量減少而有所下降。當(dāng)Γ分別從0.1增加至0.9時(shí),與未考慮柔性負(fù)荷場(chǎng)景相比,電/熱負(fù)荷曲線標(biāo)準(zhǔn)差的差值分別從44.52、30.92 kW減小為19.87、20.04 kW,負(fù)荷曲線波動(dòng)程度增大,平滑度指標(biāo)下降;而系統(tǒng)經(jīng)濟(jì)性從6.4%降低至2.31%。這是因?yàn)殡S著魯棒控制系數(shù)增大,系統(tǒng)趨于風(fēng)險(xiǎn)規(guī)避,越多考慮柔性負(fù)荷的不確定因素,運(yùn)行穩(wěn)定性有所提高,而經(jīng)濟(jì)性有所下降。系統(tǒng)魯棒性與經(jīng)濟(jì)性呈負(fù)相關(guān)變化,不同魯棒系數(shù)下的效果評(píng)價(jià)指標(biāo)可以為不同風(fēng)險(xiǎn)偏好的決策者,提供更多靈活的選擇。

    表7 不同魯棒系數(shù)下的優(yōu)化效果對(duì)比Table 7 Comparison of optimization effect under different robustness coefficients

    當(dāng)Γ=0時(shí),綜合考慮電/熱柔性負(fù)荷后的功率優(yōu)化調(diào)度情況如圖7所示。在00:00—07:00電價(jià)較低時(shí)段,CHP1和CHP2發(fā)電量減少,VPP參與電網(wǎng)購(gòu)電來(lái)滿足電負(fù)荷需求,通過(guò)GB和CHP1機(jī)組產(chǎn)熱來(lái)滿足熱負(fù)荷需求。在10:00—18:00電價(jià)較高時(shí)段,可控單元發(fā)電量提高,甚至滿發(fā),在滿足負(fù)荷需求的情況下,為追求利益的最大化,向電網(wǎng)輸送更多功率。在19:00—24:00時(shí)段,VPP與電網(wǎng)交互功率減少,主要通過(guò)優(yōu)化CHP1、CHP2機(jī)組出力來(lái)實(shí)現(xiàn)熱電聯(lián)合運(yùn)行。對(duì)電儲(chǔ)能和熱儲(chǔ)能而言,通過(guò)能量在時(shí)間軸上轉(zhuǎn)移的調(diào)峰方式,來(lái)實(shí)現(xiàn)基于分時(shí)電價(jià)的套利行為,同時(shí)可以在負(fù)荷高峰時(shí)段,減緩電網(wǎng)的供電壓力,起到一定輔助服務(wù)作用。

    圖7 VPP熱電聯(lián)合運(yùn)行功率優(yōu)化Fig.7 VPP combined heat and power operation power optimization diagram

    5 結(jié) 語(yǔ)

    本文從用戶側(cè)角度構(gòu)建了多種柔性負(fù)荷數(shù)學(xué)模型,并充分考慮系統(tǒng)內(nèi)不確定性,提出了一種計(jì)及柔性負(fù)荷的VPP熱電聯(lián)合魯棒優(yōu)化調(diào)度方法,得到如下結(jié)論:

    1)多種柔性負(fù)荷間的協(xié)調(diào)運(yùn)行對(duì)優(yōu)化負(fù)荷曲線、減緩高峰供電壓力,提高系統(tǒng)穩(wěn)定性有著重要作用,與未考慮柔性負(fù)荷場(chǎng)景相比,電/熱負(fù)荷曲線的平滑性分別提高了29.05%和49.93%,系統(tǒng)經(jīng)濟(jì)性提高了7.5%。

    2)魯棒優(yōu)化中魯棒系數(shù)有效地刻畫(huà)了多種柔性負(fù)荷的不確定性,隨著魯棒系數(shù)增加,越多考慮系統(tǒng)中的不確定因素,運(yùn)行魯棒性越高。

    3)不同魯棒系數(shù)下的電/熱柔性負(fù)荷協(xié)調(diào)優(yōu)化效果對(duì)比,為決策者提供一定參考價(jià)值,可以根據(jù)決策者的不同風(fēng)險(xiǎn)偏好,為決策實(shí)施提供更靈活的選擇。

    猜你喜歡
    魯棒熱電時(shí)段
    福州熱電兩臺(tái)660MW熱電聯(lián)產(chǎn)工程核準(zhǔn)獲批
    基于學(xué)習(xí)的魯棒自適應(yīng)評(píng)判控制研究進(jìn)展
    四個(gè)養(yǎng)生黃金時(shí)段,你抓住了嗎
    目標(biāo)魯棒識(shí)別的抗旋轉(zhuǎn)HDO 局部特征描述
    熱電轉(zhuǎn)換材料的開(kāi)發(fā)與應(yīng)用
    電力與能源(2017年6期)2017-05-14 06:19:41
    新型熱電制冷裝置的實(shí)驗(yàn)開(kāi)發(fā)
    基于Cauchy魯棒函數(shù)的UKF改進(jìn)算法
    熱泵在熱電聯(lián)產(chǎn)中的應(yīng)用
    河南科技(2015年15期)2015-03-11 16:25:52
    目標(biāo)軌跡更新的點(diǎn)到點(diǎn)魯棒迭代學(xué)習(xí)控制
    傍晚是交通事故高發(fā)時(shí)段
    亚洲精品日本国产第一区| 18禁裸乳无遮挡免费网站照片| 激情五月婷婷亚洲| 狂野欧美激情性xxxx在线观看| 啦啦啦在线观看免费高清www| 日日撸夜夜添| 国产国拍精品亚洲av在线观看| 91在线精品国自产拍蜜月| 亚洲精品第二区| 国产探花在线观看一区二区| 在线看a的网站| 成人国产麻豆网| 亚洲精品久久午夜乱码| 麻豆乱淫一区二区| 国产精品人妻久久久久久| 国产爱豆传媒在线观看| 久久人人爽人人爽人人片va| 国产黄色视频一区二区在线观看| 免费看不卡的av| 在线观看av片永久免费下载| 啦啦啦在线观看免费高清www| 久久人人爽av亚洲精品天堂 | 久久久久久久久大av| 精品国产三级普通话版| 精华霜和精华液先用哪个| 在线观看av片永久免费下载| 亚洲精品视频女| 国产女主播在线喷水免费视频网站| 中文天堂在线官网| 天堂中文最新版在线下载 | 久久精品久久久久久噜噜老黄| 免费观看a级毛片全部| 黄片wwwwww| www.色视频.com| 久久久久久久久大av| 免费av毛片视频| 免费不卡的大黄色大毛片视频在线观看| 又粗又硬又长又爽又黄的视频| 97在线视频观看| 国产免费视频播放在线视频| 99热这里只有精品一区| 久久国产乱子免费精品| 91久久精品国产一区二区三区| 国产精品99久久99久久久不卡 | 精品一区在线观看国产| 成人鲁丝片一二三区免费| 国产精品不卡视频一区二区| freevideosex欧美| 男女国产视频网站| 一级黄片播放器| 亚洲人成网站在线观看播放| 国产黄色免费在线视频| 国产高清国产精品国产三级 | 九九爱精品视频在线观看| 真实男女啪啪啪动态图| 尾随美女入室| 在线a可以看的网站| 日本黄大片高清| 日韩一本色道免费dvd| 午夜日本视频在线| 中文欧美无线码| 色吧在线观看| 久久久久久久久大av| 亚洲一级一片aⅴ在线观看| 日本猛色少妇xxxxx猛交久久| 美女xxoo啪啪120秒动态图| 69人妻影院| 亚洲欧美中文字幕日韩二区| 亚洲人与动物交配视频| 久久久精品94久久精品| 最后的刺客免费高清国语| 亚洲久久久久久中文字幕| 日韩在线高清观看一区二区三区| 天天躁夜夜躁狠狠久久av| 亚洲av福利一区| 国产成人免费无遮挡视频| 日韩 亚洲 欧美在线| 蜜桃久久精品国产亚洲av| 搞女人的毛片| 欧美亚洲 丝袜 人妻 在线| 国国产精品蜜臀av免费| 亚洲四区av| 一级爰片在线观看| 熟女电影av网| 视频区图区小说| 69人妻影院| 好男人在线观看高清免费视频| 日日啪夜夜爽| 日日摸夜夜添夜夜添av毛片| 亚洲成色77777| 街头女战士在线观看网站| xxx大片免费视频| 久久99热这里只频精品6学生| av国产久精品久网站免费入址| 久久久色成人| 秋霞在线观看毛片| 亚洲不卡免费看| 亚洲av中文字字幕乱码综合| 不卡视频在线观看欧美| 交换朋友夫妻互换小说| 精华霜和精华液先用哪个| 久久99热这里只有精品18| 男女边吃奶边做爰视频| 午夜福利视频精品| 水蜜桃什么品种好| 国产精品蜜桃在线观看| 极品少妇高潮喷水抽搐| 国产午夜精品一二区理论片| 国产一区二区三区综合在线观看 | 男人狂女人下面高潮的视频| 又黄又爽又刺激的免费视频.| 国产欧美亚洲国产| 五月玫瑰六月丁香| 女人被狂操c到高潮| 九草在线视频观看| 国产av不卡久久| 国产精品一区www在线观看| 久久久欧美国产精品| 国产精品无大码| 卡戴珊不雅视频在线播放| 97在线人人人人妻| 永久免费av网站大全| 狂野欧美激情性xxxx在线观看| 亚洲人与动物交配视频| 一级爰片在线观看| 国产成人精品福利久久| 国产永久视频网站| 国产精品国产三级国产av玫瑰| 久久99热这里只频精品6学生| 亚洲综合精品二区| 少妇熟女欧美另类| 真实男女啪啪啪动态图| 色综合色国产| 插逼视频在线观看| 男女边摸边吃奶| 国产成人免费观看mmmm| 99九九线精品视频在线观看视频| 最后的刺客免费高清国语| 国产69精品久久久久777片| 国产免费又黄又爽又色| 看免费成人av毛片| 国产永久视频网站| 高清日韩中文字幕在线| 日日摸夜夜添夜夜爱| 久久精品国产a三级三级三级| 国产午夜福利久久久久久| 欧美精品国产亚洲| 一区二区三区免费毛片| 少妇裸体淫交视频免费看高清| 男女边吃奶边做爰视频| 简卡轻食公司| av在线老鸭窝| 大陆偷拍与自拍| 少妇被粗大猛烈的视频| 黄色视频在线播放观看不卡| 国产综合懂色| 大话2 男鬼变身卡| 国产黄a三级三级三级人| 99热全是精品| 一个人看视频在线观看www免费| 最新中文字幕久久久久| 久久影院123| 久久99热6这里只有精品| 国产精品国产av在线观看| 纵有疾风起免费观看全集完整版| 国产国拍精品亚洲av在线观看| 禁无遮挡网站| 亚洲三级黄色毛片| 我的老师免费观看完整版| 亚洲自偷自拍三级| 国产 精品1| 麻豆成人av视频| av播播在线观看一区| 日日摸夜夜添夜夜爱| 亚洲图色成人| 五月开心婷婷网| 黄片wwwwww| 亚洲婷婷狠狠爱综合网| 97在线人人人人妻| 一级毛片黄色毛片免费观看视频| 国产一区亚洲一区在线观看| 亚洲四区av| 国产日韩欧美在线精品| 亚洲av日韩在线播放| 亚洲综合色惰| 国产精品熟女久久久久浪| 赤兔流量卡办理| 男的添女的下面高潮视频| 亚洲av福利一区| 日韩 亚洲 欧美在线| 亚洲在线观看片| 一级毛片黄色毛片免费观看视频| 国产精品一及| 91精品伊人久久大香线蕉| av线在线观看网站| 国产精品一区二区三区四区免费观看| 国产精品成人在线| 久久久久久久精品精品| 欧美极品一区二区三区四区| 国产成人精品福利久久| 国产一区二区在线观看日韩| 黄色怎么调成土黄色| a级毛色黄片| 夜夜爽夜夜爽视频| 女人久久www免费人成看片| 日本午夜av视频| 国产精品三级大全| 在现免费观看毛片| 国产综合懂色| 日韩成人伦理影院| 国产乱人视频| 最近最新中文字幕免费大全7| 能在线免费看毛片的网站| 特大巨黑吊av在线直播| 特级一级黄色大片| 只有这里有精品99| 国产综合精华液| 女人被狂操c到高潮| 免费观看在线日韩| 毛片一级片免费看久久久久| 欧美人与善性xxx| 亚洲av免费在线观看| 日本黄色片子视频| 久久久亚洲精品成人影院| 亚洲婷婷狠狠爱综合网| 国产一区二区三区综合在线观看 | 国产精品女同一区二区软件| 久热这里只有精品99| 又大又黄又爽视频免费| 国产视频内射| 精品人妻偷拍中文字幕| 欧美xxxx黑人xx丫x性爽| 中文资源天堂在线| 男女边吃奶边做爰视频| 尤物成人国产欧美一区二区三区| 免费观看a级毛片全部| 下体分泌物呈黄色| 欧美三级亚洲精品| 国精品久久久久久国模美| 人妻 亚洲 视频| 亚洲综合色惰| 欧美高清性xxxxhd video| freevideosex欧美| 熟女人妻精品中文字幕| 亚洲第一区二区三区不卡| 久久99热这里只频精品6学生| 成年女人在线观看亚洲视频 | 大片电影免费在线观看免费| 久久精品久久久久久久性| 免费观看av网站的网址| 久久97久久精品| 日韩一区二区视频免费看| 九九久久精品国产亚洲av麻豆| 久久精品国产亚洲网站| 22中文网久久字幕| 在现免费观看毛片| 你懂的网址亚洲精品在线观看| 午夜亚洲福利在线播放| 亚洲国产精品999| 欧美精品国产亚洲| tube8黄色片| 99久久九九国产精品国产免费| 国内精品美女久久久久久| 国产精品爽爽va在线观看网站| 99热全是精品| 午夜视频国产福利| 亚洲精品久久久久久婷婷小说| 免费电影在线观看免费观看| 99视频精品全部免费 在线| 日本一二三区视频观看| 欧美 日韩 精品 国产| 精品亚洲乱码少妇综合久久| 22中文网久久字幕| 欧美激情在线99| 国产在线一区二区三区精| 亚洲精品日韩在线中文字幕| 欧美成人a在线观看| 成年女人在线观看亚洲视频 | 丰满少妇做爰视频| 不卡视频在线观看欧美| 综合色av麻豆| 日韩电影二区| 亚洲av免费在线观看| 色视频www国产| 制服丝袜香蕉在线| 亚洲综合精品二区| 波野结衣二区三区在线| 国产高清国产精品国产三级 | 久久人人爽av亚洲精品天堂 | 日本三级黄在线观看| 成年女人在线观看亚洲视频 | 青春草亚洲视频在线观看| 日韩欧美精品免费久久| 国产视频内射| 欧美变态另类bdsm刘玥| 欧美xxⅹ黑人| 中国三级夫妇交换| 国产av国产精品国产| 欧美xxxx性猛交bbbb| 成人综合一区亚洲| 天堂中文最新版在线下载 | a级一级毛片免费在线观看| 99久久精品热视频| 蜜桃久久精品国产亚洲av| 日韩一本色道免费dvd| 亚洲国产最新在线播放| 黄色一级大片看看| 国产高清有码在线观看视频| 伦精品一区二区三区| 色哟哟·www| 涩涩av久久男人的天堂| 99热网站在线观看| 亚洲国产精品国产精品| 一区二区av电影网| 久久久久久久国产电影| 精华霜和精华液先用哪个| 亚洲欧美日韩另类电影网站 | 精华霜和精华液先用哪个| 久久久久性生活片| 亚洲精品日韩在线中文字幕| 一级黄片播放器| 我的老师免费观看完整版| 欧美日韩国产mv在线观看视频 | 九九久久精品国产亚洲av麻豆| 国产精品伦人一区二区| 国产精品无大码| 在线免费观看不下载黄p国产| 国产午夜福利久久久久久| 王馨瑶露胸无遮挡在线观看| 人人妻人人澡人人爽人人夜夜| 深爱激情五月婷婷| 久久亚洲国产成人精品v| 久久久久久久午夜电影| 国产在线男女| 一边亲一边摸免费视频| 青春草视频在线免费观看| 国产精品一及| 老女人水多毛片| 精品国产三级普通话版| 亚洲美女搞黄在线观看| 日本wwww免费看| 亚洲人成网站在线播| 久久国产乱子免费精品| 嘟嘟电影网在线观看| 只有这里有精品99| 国产69精品久久久久777片| 亚洲精品视频女| 人妻少妇偷人精品九色| 午夜福利高清视频| 精品人妻熟女av久视频| 久久久亚洲精品成人影院| 国产老妇伦熟女老妇高清| 亚洲电影在线观看av| 丰满少妇做爰视频| 免费观看在线日韩| 色婷婷久久久亚洲欧美| 亚洲欧美一区二区三区黑人 | av卡一久久| 久久鲁丝午夜福利片| 欧美一区二区亚洲| 人妻 亚洲 视频| 亚洲av.av天堂| 亚洲av在线观看美女高潮| 精品国产乱码久久久久久小说| 亚洲精品中文字幕在线视频 | 免费黄频网站在线观看国产| 亚洲精品aⅴ在线观看| 男女啪啪激烈高潮av片| 亚洲自拍偷在线| 欧美国产精品一级二级三级 | 国产精品不卡视频一区二区| 欧美xxⅹ黑人| 国产成人免费观看mmmm| 禁无遮挡网站| 久久这里有精品视频免费| 久久热精品热| 激情五月婷婷亚洲| 99热6这里只有精品| 在线观看一区二区三区激情| 午夜福利高清视频| 91久久精品国产一区二区三区| 最近2019中文字幕mv第一页| 国产精品不卡视频一区二区| 3wmmmm亚洲av在线观看| 午夜福利视频1000在线观看| 波多野结衣巨乳人妻| eeuss影院久久| 日本一二三区视频观看| 成人午夜精彩视频在线观看| 国产免费福利视频在线观看| 欧美性感艳星| 亚洲精品中文字幕在线视频 | 亚洲精品亚洲一区二区| 亚洲av免费在线观看| 成年免费大片在线观看| 久久久久国产网址| 26uuu在线亚洲综合色| 国产免费一区二区三区四区乱码| 精品人妻一区二区三区麻豆| 国产老妇女一区| 亚洲综合色惰| 国产淫片久久久久久久久| 久久精品熟女亚洲av麻豆精品| 午夜福利网站1000一区二区三区| 在线观看国产h片| 欧美日韩亚洲高清精品| 插逼视频在线观看| 国产男女内射视频| 如何舔出高潮| 身体一侧抽搐| 日韩成人av中文字幕在线观看| 免费黄频网站在线观看国产| 亚洲精品亚洲一区二区| 丰满少妇做爰视频| 久久亚洲国产成人精品v| 久久久久久伊人网av| 少妇 在线观看| 国产精品人妻久久久久久| 久久久久久伊人网av| av在线老鸭窝| 国产国拍精品亚洲av在线观看| 免费人成在线观看视频色| 国国产精品蜜臀av免费| 国产伦在线观看视频一区| 国产黄片视频在线免费观看| 人人妻人人澡人人爽人人夜夜| 18禁裸乳无遮挡免费网站照片| 亚洲电影在线观看av| av女优亚洲男人天堂| 日韩av不卡免费在线播放| 天堂俺去俺来也www色官网| 国产黄色免费在线视频| 女的被弄到高潮叫床怎么办| 色吧在线观看| 成人一区二区视频在线观看| 听说在线观看完整版免费高清| 伦精品一区二区三区| 色播亚洲综合网| 国产男女超爽视频在线观看| 久久久久网色| av在线观看视频网站免费| 特大巨黑吊av在线直播| 免费av毛片视频| 日韩视频在线欧美| 国产亚洲最大av| 欧美人与善性xxx| 69av精品久久久久久| 老司机影院毛片| 99久国产av精品国产电影| 纵有疾风起免费观看全集完整版| 三级国产精品片| 国产一区二区在线观看日韩| 最近最新中文字幕大全电影3| 国产男人的电影天堂91| 人人妻人人看人人澡| 国产黄片视频在线免费观看| 熟女电影av网| 亚洲成人久久爱视频| 欧美日韩在线观看h| 亚洲久久久久久中文字幕| 亚州av有码| 另类亚洲欧美激情| 国产免费视频播放在线视频| 大话2 男鬼变身卡| 麻豆精品久久久久久蜜桃| 美女主播在线视频| 国产精品麻豆人妻色哟哟久久| 日本猛色少妇xxxxx猛交久久| 日日摸夜夜添夜夜添av毛片| 午夜福利视频精品| 永久免费av网站大全| 青春草国产在线视频| 少妇丰满av| 男女那种视频在线观看| 国产精品三级大全| 国产一区二区三区综合在线观看 | 国产91av在线免费观看| 爱豆传媒免费全集在线观看| 在线播放无遮挡| 女人十人毛片免费观看3o分钟| 人体艺术视频欧美日本| 51国产日韩欧美| av免费观看日本| 免费少妇av软件| 2018国产大陆天天弄谢| 一个人看视频在线观看www免费| 欧美成人一区二区免费高清观看| 国产高潮美女av| 日韩一区二区三区影片| 菩萨蛮人人尽说江南好唐韦庄| 成人综合一区亚洲| 精品少妇久久久久久888优播| 欧美区成人在线视频| 乱系列少妇在线播放| 成人黄色视频免费在线看| av黄色大香蕉| 激情五月婷婷亚洲| 亚洲成色77777| 黄色视频在线播放观看不卡| 亚洲精品日韩av片在线观看| 免费黄频网站在线观看国产| 国产精品精品国产色婷婷| 网址你懂的国产日韩在线| 在现免费观看毛片| av播播在线观看一区| 日本三级黄在线观看| 一级a做视频免费观看| 日韩国内少妇激情av| 最后的刺客免费高清国语| 国产永久视频网站| 免费黄网站久久成人精品| 国产成人精品婷婷| 久久鲁丝午夜福利片| 热99国产精品久久久久久7| 久久99精品国语久久久| av在线app专区| 成人特级av手机在线观看| 少妇高潮的动态图| 亚洲国产色片| 婷婷色综合大香蕉| 黄色欧美视频在线观看| av国产精品久久久久影院| 干丝袜人妻中文字幕| 免费看a级黄色片| 最近2019中文字幕mv第一页| 交换朋友夫妻互换小说| 国产亚洲精品久久久com| 色网站视频免费| 久久韩国三级中文字幕| 日韩,欧美,国产一区二区三区| 大片免费播放器 马上看| 欧美精品人与动牲交sv欧美| 久久久久久久国产电影| 亚洲欧洲日产国产| 亚洲成色77777| 在线 av 中文字幕| 精华霜和精华液先用哪个| 夜夜爽夜夜爽视频| 久久精品久久精品一区二区三区| 亚洲av日韩在线播放| 亚洲精品456在线播放app| 自拍欧美九色日韩亚洲蝌蚪91 | 国产又色又爽无遮挡免| 王馨瑶露胸无遮挡在线观看| 日日啪夜夜撸| 免费黄频网站在线观看国产| 色视频在线一区二区三区| 国内精品美女久久久久久| 特大巨黑吊av在线直播| 亚洲av免费在线观看| 伊人久久精品亚洲午夜| 久久久久精品性色| 免费播放大片免费观看视频在线观看| 午夜视频国产福利| 亚洲国产最新在线播放| 国产午夜精品久久久久久一区二区三区| 亚洲天堂国产精品一区在线| 18禁在线播放成人免费| 大话2 男鬼变身卡| 中文字幕免费在线视频6| 有码 亚洲区| 亚洲经典国产精华液单| 久久人人爽人人片av| 日日撸夜夜添| 国产伦在线观看视频一区| 久久久色成人| 精品久久久久久久久亚洲| 久久久精品免费免费高清| 97热精品久久久久久| 只有这里有精品99| 欧美一级a爱片免费观看看| 久久97久久精品| 男女啪啪激烈高潮av片| 国产亚洲av片在线观看秒播厂| 国产精品精品国产色婷婷| 国产午夜福利久久久久久| 亚洲av不卡在线观看| 精品久久久久久久久av| av天堂中文字幕网| 天堂俺去俺来也www色官网| 亚洲av福利一区| 大香蕉久久网| 韩国av在线不卡| av又黄又爽大尺度在线免费看| 亚洲国产日韩一区二区| 97精品久久久久久久久久精品| 亚洲av.av天堂| 一个人看视频在线观看www免费| 丝袜脚勾引网站| 日韩不卡一区二区三区视频在线| 国产成人午夜福利电影在线观看| 日本一二三区视频观看| 看非洲黑人一级黄片| 国产午夜精品一二区理论片| 婷婷色综合www| 少妇被粗大猛烈的视频| 免费av观看视频| 亚洲四区av| 亚洲欧洲国产日韩| 91狼人影院| 男人爽女人下面视频在线观看| 亚洲美女视频黄频| 中国美白少妇内射xxxbb| 国产片特级美女逼逼视频| 少妇人妻久久综合中文| 内射极品少妇av片p| 乱系列少妇在线播放| 丰满人妻一区二区三区视频av| 日韩,欧美,国产一区二区三区| 在线免费观看不下载黄p国产| 亚洲自偷自拍三级| 伦理电影大哥的女人| 性色av一级| 有码 亚洲区| 日本免费在线观看一区|