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

    考慮溫控負(fù)荷聚合功率不確定性的樓宇微網(wǎng)魯棒優(yōu)化

    2021-04-09 06:14:22華光輝李晨張勇李丹劉創(chuàng)王程
    電力建設(shè) 2021年4期
    關(guān)鍵詞:魯棒微網(wǎng)溫控

    華光輝,李晨,張勇,李丹,劉創(chuàng),王程

    (1.新能源與儲(chǔ)能運(yùn)行控制國家重點(diǎn)實(shí)驗(yàn)室(中國電力科學(xué)研究院有限公司),南京市 210003;2.國家電力調(diào)度控制中心,北京市 100031;3.國網(wǎng)河北省電力有限公司,石家莊市 050024;4.國家電網(wǎng)有限公司直流建設(shè)分公司,北京市 100052;5.華北電力大學(xué)電氣與電子工程學(xué)院,北京市 102206)

    0 引 言

    隨著中國產(chǎn)業(yè)結(jié)構(gòu)的不斷升級(jí)和城市化水平的顯著提高,樓宇建筑能耗在國家總能耗中的占比持續(xù)攀升[1],其中溫控負(fù)荷(thermostatically controlled loads, TCLs)的份額可達(dá)樓宇建筑總負(fù)荷的50%以上[2]。同時(shí),由于分布式電源及儲(chǔ)能設(shè)備的大量接入,傳統(tǒng)的樓宇能耗系統(tǒng)逐漸演變?yōu)榧啥鄻怯畹亩嗄芪⒕W(wǎng)體系[3],這就為TCLs參與電網(wǎng)調(diào)度的需求響應(yīng)創(chuàng)造了條件[4]。TCLs設(shè)備能在數(shù)分鐘內(nèi)快速響應(yīng),而對用戶的體驗(yàn)影響很小,因而在電網(wǎng)運(yùn)行風(fēng)險(xiǎn)快速防范體系中的地位愈加重要。

    集成樓宇微網(wǎng)系統(tǒng)TCLs具有分布分散、響應(yīng)不確定性強(qiáng)的特點(diǎn)[5],建立實(shí)用的TCLs聚合模型是微網(wǎng)系統(tǒng)對其進(jìn)行調(diào)控的基礎(chǔ)。目前有代表性的聚合方法可以分為2類。第一類是采用近似數(shù)學(xué)方法:文獻(xiàn)[6]通過耦合的Fokker-Planck方程描述TCLs的動(dòng)態(tài)變化,并基于數(shù)值拉普拉斯反變換進(jìn)行求解;文獻(xiàn)[7]將基于一階等效熱參數(shù)(equivalent thermal parameter, ETP)模型的蒙特卡洛法和狀態(tài)隊(duì)列模型相結(jié)合,并對聚合功率的跌落過程進(jìn)行分析;文獻(xiàn)[8]將整個(gè)控制溫度范圍劃分為多個(gè)等間距小段,利用離散化差分法構(gòu)建出聚合TCLs的雙線性狀態(tài)空間模型;文獻(xiàn)[9]將參數(shù)的概率分布引入溫控負(fù)荷的二階ETP動(dòng)態(tài)模型,建立了一種適用于異質(zhì)性TCLs聚合的二維狀態(tài)轉(zhuǎn)移模型??傮w而言,以上文獻(xiàn)均定位在TCLs聚合商的角度,控制模型相對復(fù)雜,不易于集成樓宇微網(wǎng)系統(tǒng)的直接使用。

    第二類是使用物理建模軟件或聚類的方法:文獻(xiàn)[10]通過EnergyPlus軟件模擬建筑物各個(gè)子系統(tǒng)間的復(fù)雜熱交換過程以及對需求響應(yīng)指令的反饋,但當(dāng)應(yīng)用于包含眾多樓宇的集成微網(wǎng)系統(tǒng)時(shí),計(jì)算代價(jià)過大;文獻(xiàn)[11]采用自適應(yīng)的模糊C均值聚類法對區(qū)域TCLs進(jìn)行聚合求取,但沒有在微網(wǎng)調(diào)度層面對不確定性進(jìn)行處理;文獻(xiàn)[12]將EnergyPlus模擬和聚類算法結(jié)合,通過EnergyPlus軟件模擬城市區(qū)域的建筑熱負(fù)荷,再利用k-means聚類算法進(jìn)行聚合分析,并在調(diào)度層面以不確定集模擬其波動(dòng)性,然而該方法不能有效計(jì)及樓宇用戶的舒適度。

    針對以上2類方法所存在的問題,文獻(xiàn)[13]采用核平滑技術(shù)以獲得多個(gè)溫控負(fù)荷聚合后的分布函數(shù);文獻(xiàn)[14]推導(dǎo)出一種可由室外溫度和溫控負(fù)荷一階ETP模型直接求解的TCLs近似聚合模型,無需復(fù)雜的控制原理且易于進(jìn)行需求響應(yīng)潛力評估。雖然這2種方法的控制原理簡單且計(jì)算代價(jià)較低,但都忽略樓宇建筑的具體結(jié)構(gòu),難以適應(yīng)集成樓宇微網(wǎng)中大量建筑的差異化建模需求。

    本文在文獻(xiàn)[14]的基礎(chǔ)上,提出一種基于樓宇熱阻熱容(resistance-capacitance, RC)模型[15]的TCLs聚合模型,并將其集成到區(qū)域樓宇微網(wǎng)系統(tǒng)中進(jìn)行調(diào)度分析。所提出的溫控負(fù)荷模型充分計(jì)及樓宇建筑的物理結(jié)構(gòu),并通過線性假設(shè)將樓宇傳熱的非線性傳熱方程進(jìn)行簡化,得到的TCLs聚合方案較為簡便且能適應(yīng)樓宇建筑的差異化建模需求?;谒酺CLs聚合方案,文章從微網(wǎng)運(yùn)營商的角度,建立一種能有效應(yīng)對可再生分布式電源(distributed generator, DG)出力和TCLs聚合功率不確定性的兩階段魯棒優(yōu)化模型。

    1 樓宇微網(wǎng)的溫控負(fù)荷聚合模型

    1.1 樓宇單區(qū)域RC傳熱模型

    典型的樓宇單區(qū)域RC傳熱模型如圖1所示,將墻體和室內(nèi)、外環(huán)境列為網(wǎng)絡(luò)節(jié)點(diǎn),分別用Cw和Cin表示由墻體厚度和室內(nèi)空間所導(dǎo)致的熱容效應(yīng),狀態(tài)變量為節(jié)點(diǎn)溫度。需要說明的是,由于建筑物內(nèi)部房間的熱交換是相互的,因此對于整棟樓宇的熱/冷需求聚合,可忽略樓宇內(nèi)部的熱交換過程。

    圖1 典型樓宇的單區(qū)域RC網(wǎng)絡(luò)模型Fig.1 Single-user RC network model for building

    單區(qū)域用戶的墻體熱平衡約束條件如下:

    (1)

    單區(qū)域用戶的室內(nèi)熱平衡約束條件如下:

    (2)

    (3)

    式(3)表示夏季制冷條件下,當(dāng)室內(nèi)溫度升至Tset+δ/2時(shí),開啟溫控設(shè)備;降至Tset-δ/2時(shí),關(guān)閉溫控設(shè)備。冬季制熱條件下正好相反。因此,溫控設(shè)備的開關(guān)狀態(tài)會(huì)周期性地發(fā)生跳變,導(dǎo)致室內(nèi)溫度周期性地在設(shè)定界值之間變動(dòng),如圖2所示(以夏季制冷情況為例)。

    圖2 溫控負(fù)荷運(yùn)行狀態(tài)示意圖Fig.2 Schematic diagram of the operation states of thermostatically controlled loads

    1.2 單區(qū)域傳熱模型的線性化

    由式(1)、(2)組成的常微分方程組求解Tin,會(huì)引入二階常系數(shù)非齊次微分方程,導(dǎo)致出現(xiàn)2個(gè)特征值而難以求解溫控負(fù)荷處于開機(jī)或關(guān)機(jī)狀態(tài)的時(shí)間。針對這一情況,對式(1)進(jìn)行有限差分:

    (4)

    式中:Δt表示差分的時(shí)間步長,在本文中即為一個(gè)仿真時(shí)段;Tw(t-Δt)表示在Tw前Δt時(shí)刻的墻體溫度。

    由于式(4)涉及多個(gè)時(shí)段,文章采取文獻(xiàn)[16]中的線性假設(shè),可將Tw改寫為同一時(shí)刻的室內(nèi)、外溫度的函數(shù):

    Tw=a1Tin+a2Tout+a3

    (5)

    (6)

    將式(5)代入到式(2)中,將式(2)改寫成為僅含Tin一個(gè)狀態(tài)變量的一階非齊次線性微分方程:

    (7)

    (8)

    通過求解式(7),可以分別得出溫控負(fù)荷的單次開機(jī)時(shí)間Tc和關(guān)機(jī)時(shí)間Th:

    (9)

    1.3 溫控負(fù)荷聚合模型

    對于微網(wǎng)運(yùn)營商而言,只有大量溫控負(fù)荷聚合后的總功率才有參考價(jià)值。當(dāng)樓宇樣本的數(shù)目足夠多且相互獨(dú)立時(shí),依據(jù)大數(shù)定律,樣本均值逐漸逼近整體的均值。因此,N個(gè)樓宇區(qū)域聚合后的溫控設(shè)備總制熱(冷)功率hagg為式(10)的形式。

    (10)

    (11)

    式中:pon,i表示樓宇區(qū)域i的溫控設(shè)備處于啟動(dòng)狀態(tài)的概率;hint,i、Tc,i、Th,i分別表示樓宇區(qū)域i的溫控設(shè)備制熱(冷)量、單次開機(jī)時(shí)間、單次關(guān)機(jī)時(shí)間。

    將式(9)代入式(11),pon,i的范圍可以表示為[14]:

    (12)

    (13)

    式中:Pint,i、ηi分別為溫控負(fù)荷的耗電功率和轉(zhuǎn)換效率,滿足hint,i=Pint,i·ηi。

    由式(11)—(13)可以近似得到N個(gè)樓宇區(qū)域聚合后的總溫控負(fù)荷的下界Pagg,min和上界Pagg,max:

    (14)

    (15)

    無論對于同質(zhì)性還是異構(gòu)性的樓宇溫控負(fù)荷,式(14)、(15)都能成立,具有較強(qiáng)的可擴(kuò)展性。據(jù)此,智能樓宇溫控負(fù)荷的聚合功率Pagg可由估值區(qū)間[Pagg,min,Pagg,max]表示如下:

    Pagg=αPagg,min+(1-α)Pagg,max,α∈[0,1]

    (16)

    式中:α為加權(quán)系數(shù)。

    2 計(jì)及溫控負(fù)荷聚合功率不確定性的魯棒優(yōu)化

    在樓宇微網(wǎng)的實(shí)際運(yùn)行中,受到通信設(shè)備的通達(dá)度、溫控設(shè)備的可控度以及用戶對參與需求響應(yīng)的可接受度等因素的影響,樓宇集群的實(shí)際TCLs聚合功率存在不確定性。同時(shí),樓宇微網(wǎng)系統(tǒng)中存在的小型風(fēng)機(jī)、光伏等分布式可再生能源設(shè)備的出力與風(fēng)速、光照條件強(qiáng)相關(guān),同樣需要計(jì)及其不確定性的影響。因此,文章建立了一種適用于集成樓宇微網(wǎng)的兩階段魯棒優(yōu)化模型。第一階段問題是根據(jù)日前的可再生DG預(yù)測出力及確定性的TCLs聚合模型,使得微網(wǎng)運(yùn)營商的經(jīng)濟(jì)性最優(yōu)。第二階段問題則利用不確定性集合模擬可能出現(xiàn)的最惡劣運(yùn)行工況,并求解該最惡劣運(yùn)行工況下的聯(lián)絡(luò)線功率的最小調(diào)整量,以檢測第一階段問題所得機(jī)組計(jì)劃的可行性。兩階段之間的交互原理如圖3 所示,每一階段問題的具體構(gòu)成和求解算法詳見2.2至2.5節(jié)。

    圖3 魯棒調(diào)度模型兩階段間的交互原理圖Fig.3 Schematic diagram of the interaction between the two stages of the robust scheduling model

    2.1 集成樓宇微網(wǎng)系統(tǒng)的運(yùn)營模式

    集成樓宇微網(wǎng)系統(tǒng)主要分為微網(wǎng)DG、蓄電池設(shè)備以及樓宇用戶集群3部分,分別對應(yīng)源、儲(chǔ)、荷3種微網(wǎng)資源。微網(wǎng)運(yùn)營商采用“供-儲(chǔ)-需”三者的自平衡集中管理策略對本區(qū)域的資源進(jìn)行全局協(xié)調(diào),建立樓宇微網(wǎng)系統(tǒng)綜合交易平臺(tái)并與上級(jí)配電網(wǎng)進(jìn)行電力交易,其運(yùn)營模式如圖4所示。

    圖4 集成多樓宇微網(wǎng)系統(tǒng)的運(yùn)營模式Fig.4 Integrative operation mode of the microgrid system with multiple smart buildings integrated

    微網(wǎng)運(yùn)營商通過微網(wǎng)能量管理系統(tǒng)對微網(wǎng)內(nèi)的供、儲(chǔ)、需資源進(jìn)行統(tǒng)一管理,并在微網(wǎng)DG側(cè)、儲(chǔ)能設(shè)備側(cè)和樓宇集群側(cè)分設(shè)子能量管理模塊。每一個(gè)子能量管理模塊接收微網(wǎng)能量管理系統(tǒng)下發(fā)的調(diào)度指令,對所屬的微網(wǎng)資源進(jìn)行調(diào)控,同時(shí)監(jiān)測所轄范圍內(nèi)的天氣參數(shù)、樓宇用戶舒適性需求等,并將相關(guān)的信息反饋到上級(jí)微網(wǎng)能量管理系統(tǒng)。微網(wǎng)運(yùn)營商按照“自發(fā)自用,余電上網(wǎng)”的運(yùn)營模式,與上級(jí)配電網(wǎng)進(jìn)行購、售電交易,靈活利用儲(chǔ)能設(shè)備的錯(cuò)峰調(diào)節(jié)和用戶熱(冷)需求彈性最小化運(yùn)營成本。

    2.2 第一階段問題——經(jīng)濟(jì)性最優(yōu)

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

    第一階段問題的目標(biāo)是最小化微網(wǎng)的運(yùn)營成本CFir,包括向配網(wǎng)購電成本、可控DG的燃料及啟動(dòng)成本以及各種設(shè)備的維護(hù)成本等,如式(17)所示。

    (17)

    2.2.2約束條件

    1)微電網(wǎng)功率平衡約束如下:

    (18)

    式中:PD,t表示與熱(冷)需求無關(guān)的非耦合性電負(fù)荷。

    2)微網(wǎng)系統(tǒng)購電約束如下:

    (19)

    3)微網(wǎng)可控DG運(yùn)行約束如下:

    Pi,minui,t≤Pi,t≤Pi,maxui,t

    (20)

    (21)

    (22)

    (23)

    (24)

    -ui,(t-1)+ui,t-zi,t≤0

    (25)

    4)微網(wǎng)蓄電池運(yùn)行約束如下:

    (26)

    (27)

    (28)

    5)樓宇用戶舒適性約束條件如下:

    樓宇溫控負(fù)荷的聚合功率相關(guān)約束條件對應(yīng)于式(14)—(16), 微網(wǎng)系統(tǒng)溫控負(fù)荷聚合功率Pagg,t是關(guān)于溫度設(shè)定值Tset,t的線性函數(shù)。此外,Tset,t的范圍受到如下約束條件的限制:

    (29)

    -ΔTset≤Tset,t+1-Tset,t≤ΔTset

    (30)

    2.3 不確定性建模

    2.3.1可再生DG的不確定性

    建立表征可再生DG不確定性的集合US:

    (31)

    (32)

    (33)

    (34)

    2.3.2溫控負(fù)荷聚合功率不確定性

    第一階段規(guī)劃求解完成后,將所得TCLs聚合功率Pagg,t作為已知參數(shù)傳送給第二階段問題,建立表征TCLs聚合功率不確定性的集合UL,如下所示:

    (35)

    (36)

    (37)

    2.4 第二階段問題——調(diào)整量最少

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

    第二階段問題分為內(nèi)、外兩層,當(dāng)外層max問題的不確定性集合US&UL選取為某一組不確定性參數(shù)時(shí),內(nèi)層min問題的目標(biāo)是求解在該組不確定性參數(shù)下的微網(wǎng)系統(tǒng)不平衡功率的最小值。接著,外層max問題選取US&UL集合內(nèi)的另一組不確定性參數(shù),再次求解內(nèi)層min問題,以此類推。外層max問題的目標(biāo)就是尋找出US&UL集合內(nèi)使得內(nèi)層min問題所解得的不平衡功率最小值取值最大的那一組不確定性參數(shù),即令系統(tǒng)運(yùn)行工況最惡劣的一組不確定性參數(shù),其所對應(yīng)的內(nèi)層min問題的不平衡功率最小值即為第二階段max (min)問題的解。

    (38)

    2.4.2約束條件

    1)微電網(wǎng)功率平衡約束如下:

    (39)

    2)微網(wǎng)系統(tǒng)購電約束如下:

    (40)

    3)其他約束如下:

    第二階段問題的其他約束條件與第一階段問題中相同,包括微網(wǎng)可控DG運(yùn)行約束式(20)—(22)、微網(wǎng)蓄電池運(yùn)行約束式(26)—(28),區(qū)別僅在于用相應(yīng)的考慮不確定性后帶‘Sec’上標(biāo)的修正值代替原有的形式。

    2.5 模型求解

    2.2至2.4節(jié)中所建立的是一個(gè)兩階段的魯棒優(yōu)化模型,第一階段問題是一個(gè)min問題,第二階段問題是一個(gè)max-min問題,兩階段間相互影響,不能直接同時(shí)求解。因此,文章采用廣泛應(yīng)用于求解兩階段魯棒優(yōu)化問題的列和約束生成算法(column and constraint generation, C & CG)進(jìn)行迭代求解,具體的第二階段max-min問題的線性化、兩階段間參數(shù)傳遞可參見文獻(xiàn)[17]。通過C & CG算法,可以將整個(gè)問題轉(zhuǎn)化為一個(gè)混合整數(shù)線性規(guī)劃問題,采用Gurobi等商用求解器即可求解。

    3 算例分析

    3.1 智能樓宇溫控負(fù)荷聚合模型的驗(yàn)證

    3.1.1有效性驗(yàn)證

    在夏季為樓宇用戶供冷的前提條件下,采用5 000臺(tái)空調(diào)設(shè)備,分別為5 000個(gè)智能樓宇區(qū)域供冷。溫控設(shè)備和樓宇區(qū)域的參數(shù)在表1所示的數(shù)值區(qū)間內(nèi)服從均勻分布,而各溫控設(shè)備的初始開關(guān)狀態(tài)隨機(jī),以此進(jìn)行蒙特卡洛仿真的結(jié)果相加后作為實(shí)際功率。然后,通過式(14)、(15)和(16)得到TCLs的聚合功率的估算值,取外界溫度為32 ℃,太陽輻照強(qiáng)度為263 W/m2,加權(quán)系數(shù)系數(shù)α為0.5,智能樓宇溫控負(fù)荷聚合功率估計(jì)值及上下界結(jié)果如圖5所示。

    表1 溫控負(fù)荷和智能樓宇的參數(shù)范圍Table 1 Ranges of parameters of thermostatically controlled loads and Intelligent building

    圖5 智能樓宇溫控負(fù)荷聚合功率估計(jì)值及上下界Fig.5 Estimated aggregated power of TCLs and its upper and lower limits

    由圖5可得,將5 000個(gè)樓宇區(qū)域的溫控負(fù)荷聚合后的實(shí)際穩(wěn)定功率的運(yùn)行區(qū)間約為2.11~2.31 MW,而本文中提出的TCLs聚合方法所求得的功率上限和功率下限能較好地包絡(luò)實(shí)際聚合功率。當(dāng)取加權(quán)系數(shù)α=0.5時(shí),所得的聚合功率估算值為2.18 MW,和實(shí)際聚合功率的相對誤差如圖6所示。

    從圖6(a)中可以看出,所提模型估算的聚合功率和實(shí)際聚合功率的相對誤差在5%以內(nèi);圖6(b) 給出了相對誤差的分布圖,近似滿足正態(tài)分布的規(guī)律。由此,驗(yàn)證了所提聚合模型的有效性,可以滿足微網(wǎng)調(diào)度決策的精度要求。

    圖6 聚合功率的相對誤差分布Fig.6 Distribution of relative error of integrated power

    3.1.2外界天氣參數(shù)不同

    圖7給出了在不同天氣條件下的智能樓宇TCLs的聚合結(jié)果對比(一組只對比一種天氣參數(shù),另一種保持不變)。其中,圖7(a)反映了在32 ℃、30 ℃和28 ℃這3個(gè)室外溫度條件下的估算結(jié)果,圖7(b)則在3種設(shè)定的光照強(qiáng)度(100%、60%、20%)下對所提模型的估算值和實(shí)際聚合功率進(jìn)行對比。

    圖7 不同天氣條件下的樓宇溫控負(fù)荷聚合功率Fig.7 Aggregated power under different weather conditions

    結(jié)果表明:1)本文中所提出的樓宇溫控負(fù)荷聚合方案對于各種不同的天氣狀況都有較強(qiáng)的適應(yīng)性,估計(jì)功率的上下界值總能較好地包絡(luò)實(shí)際功率區(qū)間。2)當(dāng)室外溫度升高或光照強(qiáng)度增加時(shí),樓宇溫控負(fù)荷聚合功率也會(huì)隨之增加。這是因?yàn)楫?dāng)室外溫度或光照強(qiáng)度升高時(shí),樓宇區(qū)域單位面積在單位時(shí)間內(nèi)吸收的外界熱量增加,而溫控設(shè)備的功率恒定,因此室內(nèi)溫度下降到設(shè)定溫度下界的時(shí)間變長,而上升到設(shè)定溫度上界的時(shí)間反而縮短,這就會(huì)使溫控設(shè)備的占空比增加,由式(10)、(11)可知,樓宇用戶的TCLs聚合功率值會(huì)相應(yīng)提高。

    3.2 智能樓宇微網(wǎng)系統(tǒng)魯棒優(yōu)化效果分析

    3.2.1微網(wǎng)系統(tǒng)的算例數(shù)據(jù)

    以中國北方夏季某典型日為例進(jìn)行算例分析,其室外溫度、太陽輻照強(qiáng)度、分時(shí)電價(jià)的日前預(yù)測值以及微網(wǎng)系統(tǒng)的各項(xiàng)經(jīng)濟(jì)參數(shù)見文獻(xiàn)[18]。微網(wǎng)系統(tǒng)中可再生DG的預(yù)測出力如圖8所示,可控DG的運(yùn)行參數(shù)如表2所示,其中風(fēng)機(jī)和光伏各有4臺(tái)。樓宇建筑參數(shù)和溫控設(shè)備的參數(shù)同表1,聚合規(guī)模為5 000。

    圖8 微網(wǎng)系統(tǒng)可再生DG出力預(yù)測值Fig.8 Day-ahead forecasting value of the distributed renewable generators in the microgrid

    表2 微網(wǎng)可控DG的運(yùn)行參數(shù)Table 2 Operation parameters of the dispatchable DGs

    3.2.2魯棒優(yōu)化的結(jié)果分析

    圖9給出了在兩種調(diào)度模型下的微網(wǎng)機(jī)組計(jì)劃對比。由圖9可得,在確定型方案中,蓄電池的充、放電功率與分時(shí)電價(jià)基本同步,在低谷電價(jià)時(shí)(00:00—06:00)充電,在高峰電價(jià)時(shí)(14:00—18:00)放電,且基本都以最大充電或放電功率運(yùn)行;而在魯棒調(diào)度方案中,為了應(yīng)對不確定性帶來的影響,蓄電池的充放電曲線處于非最大充電或放電功率運(yùn)行的過渡時(shí)段的時(shí)長明顯增多。另一方面,11:00—12:00以及14:00—17:00正好是一日內(nèi)溫度較高,用戶TCLs聚合功率較大,且風(fēng)機(jī)和光伏功率也較大的時(shí)段,不確定性的影響最為顯著,因此,可控DG的開機(jī)數(shù)以及運(yùn)行出力相較于確定型的模型明顯增加。

    圖9 兩種調(diào)度運(yùn)行方案下的微網(wǎng)機(jī)組計(jì)劃對比Fig.9 Scheduling strategy of microgrid under different methods

    相較于蓄電池,柴油發(fā)電機(jī)和燃料電池的出力總和變化較小。這是因?yàn)閷τ谙到y(tǒng)中存在的分布式可再生能源出力和溫控負(fù)荷聚合功率的不確定性,蓄電池和可控分布式發(fā)電設(shè)備(柴油發(fā)電機(jī)和燃料電池)都能有效應(yīng)對,兩者的作用有所重疊。相較于柴油發(fā)電機(jī)和燃料電池,蓄電池的運(yùn)行成本更低,無需考慮開機(jī)費(fèi)用且不受爬坡功率約束的限制,經(jīng)濟(jì)性和調(diào)控靈活性上都具有顯著優(yōu)勢。因而在本文所提的魯棒調(diào)度模型下,蓄電池的出力變化比柴油發(fā)電機(jī)和燃料電池更加劇烈。然而,雖然總體的變化幅度較小,單個(gè)分布式發(fā)電設(shè)備在兩種優(yōu)化模型下的啟停狀態(tài)和機(jī)組出力可能相差很大,因此本文所提魯棒調(diào)度模型對可控分布式發(fā)電設(shè)備的影響不可忽略。

    為了驗(yàn)證所提魯棒調(diào)度方案在樓宇微網(wǎng)系統(tǒng)實(shí)際運(yùn)行中的經(jīng)濟(jì)性和魯棒性,利用蒙特卡洛法采樣1 000個(gè)使得TCLs聚合功率、可再生DG預(yù)測出力誤差滿足正態(tài)分布的模擬場景,取其標(biāo)準(zhǔn)差分別為可再生DG預(yù)測出力和TCLs聚合功率日前規(guī)劃值的1/6。兩種調(diào)度方案在1 000個(gè)采樣場景下的仿真結(jié)果如表3所示,由表3可知,采用本文所提的魯棒調(diào)度方法時(shí)不會(huì)有不可行場景出現(xiàn),而確定型的模型則有93個(gè)場景必須調(diào)用緊急備用電源或進(jìn)行強(qiáng)制切去部分負(fù)荷。對于可行的場景而言,魯棒調(diào)度方案的運(yùn)行成本、開機(jī)成本和維護(hù)成本均高于確定型的模型,但其聯(lián)絡(luò)線功率的調(diào)整懲罰量遠(yuǎn)低于確定型模型,而這對于維持微電網(wǎng)的穩(wěn)定運(yùn)行以及參與電力市場運(yùn)營具有顯著意義。

    表3 兩種調(diào)度方案在1 000個(gè)采樣場景下的結(jié)果Table 3 Simulation results with different methods under 1 000 sampled scenarios

    4 結(jié) 論

    通過對樓宇單區(qū)域RC傳熱模型進(jìn)行線性化處理,提出了一種充分計(jì)及樓宇建筑物理結(jié)構(gòu)的溫控負(fù)荷近似聚合方案,并根據(jù)這一方案建立了一種適用于集成樓宇微網(wǎng)系統(tǒng)的兩階段魯棒優(yōu)化模型。通過算例分析得出如下結(jié)論:

    1)所提出的區(qū)域樓宇溫控負(fù)荷聚合方案對于各種不同的室外溫度和光照強(qiáng)度都能適用,其估計(jì)功率的上下界值總能很好地包絡(luò)實(shí)際聚合功率。當(dāng)取加權(quán)系數(shù)為0.5時(shí), 所提聚合模型的功率估算值和實(shí)際聚合功率的相對誤差維持在5%以內(nèi),可以滿足微網(wǎng)調(diào)度決策的精度要求。

    2)建立了一種考慮溫控負(fù)荷聚合功率及分布式可再生能源出力不確定性的兩階段魯棒優(yōu)化模型。相對于確定型的調(diào)度模型,所提模型能夠顯著提高機(jī)組計(jì)劃的可行性并減少再調(diào)度過程中的聯(lián)絡(luò)線功率調(diào)整量,這對于維持微電網(wǎng)的穩(wěn)定以及參與電力市場運(yùn)營具有顯著意義。雖然求解時(shí)間顯著增加,但足夠滿足日前規(guī)劃調(diào)度的時(shí)效性要求。

    猜你喜歡
    魯棒微網(wǎng)溫控
    溫控/光控片上納米裂結(jié)
    基于學(xué)習(xí)的魯棒自適應(yīng)評判控制研究進(jìn)展
    目標(biāo)魯棒識(shí)別的抗旋轉(zhuǎn)HDO 局部特征描述
    基于MSP430和Android溫控距控智能電風(fēng)扇設(shè)計(jì)
    電子制作(2016年15期)2017-01-15 13:39:06
    基于OMAP-L138的微網(wǎng)控制器設(shè)計(jì)
    骨料自動(dòng)溫控技術(shù)
    基于ARM7 和FPGA的風(fēng)扇散熱系統(tǒng)溫控模塊設(shè)計(jì)及實(shí)現(xiàn)
    基于Cauchy魯棒函數(shù)的UKF改進(jìn)算法
    基于改進(jìn)下垂法的微網(wǎng)并網(wǎng)控制策略研究
    目標(biāo)軌跡更新的點(diǎn)到點(diǎn)魯棒迭代學(xué)習(xí)控制
    欧美不卡视频在线免费观看| 日本成人三级电影网站| 亚洲在线自拍视频| 久久99热6这里只有精品| 桃红色精品国产亚洲av| 午夜亚洲福利在线播放| av在线天堂中文字幕| 国产一区二区三区在线臀色熟女| 十八禁网站免费在线| 精品午夜福利在线看| 国产一区二区激情短视频| 免费观看人在逋| 夜夜爽天天搞| 免费观看的影片在线观看| 国产亚洲av嫩草精品影院| 成年人黄色毛片网站| 欧美色视频一区免费| 欧美国产日韩亚洲一区| 中文资源天堂在线| 久久久久久九九精品二区国产| 91麻豆精品激情在线观看国产| 亚洲aⅴ乱码一区二区在线播放| 国产中年淑女户外野战色| 淫妇啪啪啪对白视频| 少妇熟女aⅴ在线视频| 看片在线看免费视频| 日本 欧美在线| av女优亚洲男人天堂| 啦啦啦啦在线视频资源| 久久人妻av系列| 大又大粗又爽又黄少妇毛片口| 老师上课跳d突然被开到最大视频| 久久久久久久久久久丰满 | 狂野欧美白嫩少妇大欣赏| 真人一进一出gif抽搐免费| 成年人黄色毛片网站| 1024手机看黄色片| 非洲黑人性xxxx精品又粗又长| 国产日本99.免费观看| 成人二区视频| 国产av一区在线观看免费| 一本精品99久久精品77| 成人精品一区二区免费| 亚洲中文字幕一区二区三区有码在线看| 久久久午夜欧美精品| 婷婷精品国产亚洲av| 久久精品国产99精品国产亚洲性色| 成人亚洲精品av一区二区| av在线老鸭窝| 国产午夜精品论理片| 日本黄色视频三级网站网址| 2021天堂中文幕一二区在线观| 亚洲av免费在线观看| 免费观看人在逋| 精品欧美国产一区二区三| 蜜桃亚洲精品一区二区三区| 日韩一区二区视频免费看| 久久精品国产亚洲av涩爱 | 日韩av在线大香蕉| www.www免费av| 美女cb高潮喷水在线观看| 男人和女人高潮做爰伦理| 久久久国产成人免费| 97超级碰碰碰精品色视频在线观看| 久久草成人影院| 又爽又黄a免费视频| 久久精品国产鲁丝片午夜精品 | 国产精品一及| 精品久久国产蜜桃| 欧美成人a在线观看| 精品一区二区三区视频在线| 又爽又黄无遮挡网站| 国产 一区 欧美 日韩| 中文资源天堂在线| 婷婷精品国产亚洲av在线| 舔av片在线| 国产69精品久久久久777片| 国产午夜精品论理片| 日韩欧美精品v在线| 欧美日韩乱码在线| 亚洲欧美精品综合久久99| 搞女人的毛片| 91av网一区二区| 狠狠狠狠99中文字幕| 成人亚洲精品av一区二区| 黄色配什么色好看| 色噜噜av男人的天堂激情| 狠狠狠狠99中文字幕| 亚洲人成网站在线播| 国内精品一区二区在线观看| 婷婷色综合大香蕉| 91在线精品国自产拍蜜月| 精品99又大又爽又粗少妇毛片 | 18禁裸乳无遮挡免费网站照片| 99久久成人亚洲精品观看| 欧美区成人在线视频| 国产私拍福利视频在线观看| 给我免费播放毛片高清在线观看| 在线播放国产精品三级| 国产高清不卡午夜福利| 99热只有精品国产| 欧美色视频一区免费| 真人做人爱边吃奶动态| 国产乱人视频| 露出奶头的视频| 美女 人体艺术 gogo| 色5月婷婷丁香| 精品久久久久久成人av| 搞女人的毛片| 色噜噜av男人的天堂激情| 亚洲国产精品成人综合色| 日本在线视频免费播放| 三级男女做爰猛烈吃奶摸视频| 一进一出好大好爽视频| 岛国在线免费视频观看| 成人av在线播放网站| 我要搜黄色片| 国产精品一及| 欧美性感艳星| 好男人在线观看高清免费视频| 久久久久精品国产欧美久久久| 国产精品不卡视频一区二区| 搞女人的毛片| .国产精品久久| 婷婷亚洲欧美| 久久精品91蜜桃| 日本五十路高清| 丰满人妻一区二区三区视频av| 琪琪午夜伦伦电影理论片6080| 九九热线精品视视频播放| 日韩高清综合在线| 日本撒尿小便嘘嘘汇集6| 国产精品一区二区免费欧美| 欧美+日韩+精品| 色哟哟哟哟哟哟| 丰满的人妻完整版| 国产精品精品国产色婷婷| 免费观看的影片在线观看| 精品99又大又爽又粗少妇毛片 | 久久精品国产自在天天线| 在线观看一区二区三区| 久久久久久久久大av| 日本a在线网址| 欧美色欧美亚洲另类二区| 韩国av一区二区三区四区| 少妇人妻一区二区三区视频| 91久久精品国产一区二区成人| 国产在视频线在精品| 悠悠久久av| 亚洲精品国产成人久久av| 干丝袜人妻中文字幕| 国产乱人伦免费视频| www.色视频.com| 禁无遮挡网站| 午夜视频国产福利| 日韩欧美精品v在线| 色精品久久人妻99蜜桃| 精品一区二区三区av网在线观看| 97碰自拍视频| 校园人妻丝袜中文字幕| 精品人妻一区二区三区麻豆 | 精品人妻1区二区| 欧美成人性av电影在线观看| 国产在视频线在精品| 亚洲 国产 在线| 搡老妇女老女人老熟妇| 99久久成人亚洲精品观看| 小说图片视频综合网站| 亚洲无线在线观看| 婷婷丁香在线五月| 国产熟女欧美一区二区| 99在线人妻在线中文字幕| 99热这里只有是精品在线观看| 久久久久国产精品人妻aⅴ院| www.色视频.com| 在线免费十八禁| 不卡视频在线观看欧美| 他把我摸到了高潮在线观看| 国产美女午夜福利| 日本熟妇午夜| 91麻豆精品激情在线观看国产| 国产成人av教育| 久久精品国产亚洲av涩爱 | 黄片wwwwww| 美女cb高潮喷水在线观看| 免费av观看视频| 少妇丰满av| 久9热在线精品视频| 午夜视频国产福利| 亚洲美女视频黄频| 午夜爱爱视频在线播放| 日韩国内少妇激情av| www.色视频.com| 少妇人妻精品综合一区二区 | 欧洲精品卡2卡3卡4卡5卡区| 国产精品一区二区性色av| 欧美性感艳星| 真实男女啪啪啪动态图| 99国产极品粉嫩在线观看| 久久久久国内视频| 欧美高清成人免费视频www| 亚洲一区二区三区色噜噜| 婷婷六月久久综合丁香| 天美传媒精品一区二区| 麻豆精品久久久久久蜜桃| 国产高清三级在线| 一级黄色大片毛片| 搡女人真爽免费视频火全软件 | 97超视频在线观看视频| 99久久无色码亚洲精品果冻| 性欧美人与动物交配| 亚洲精华国产精华液的使用体验 | 性欧美人与动物交配| 久久久久久久久中文| 麻豆久久精品国产亚洲av| 久久99热这里只有精品18| 天美传媒精品一区二区| 最新在线观看一区二区三区| 99久久精品一区二区三区| 午夜福利在线观看免费完整高清在 | 日本a在线网址| 一级黄片播放器| 精品国内亚洲2022精品成人| 日本黄色视频三级网站网址| 国产精品免费一区二区三区在线| 亚洲国产精品成人综合色| 99久国产av精品| 亚洲无线观看免费| 国产女主播在线喷水免费视频网站 | 干丝袜人妻中文字幕| 亚洲人成网站在线播放欧美日韩| 色综合站精品国产| 亚洲精品乱码久久久v下载方式| 91久久精品电影网| 成人午夜高清在线视频| 国产成人a区在线观看| 97人妻精品一区二区三区麻豆| 嫁个100分男人电影在线观看| 韩国av一区二区三区四区| 日日撸夜夜添| 国产单亲对白刺激| 麻豆久久精品国产亚洲av| av视频在线观看入口| 久久久精品大字幕| 国产不卡一卡二| 久久久国产成人免费| 精品福利观看| 淫秽高清视频在线观看| 久久国产精品人妻蜜桃| 国产三级在线视频| 国产亚洲精品久久久com| 老司机深夜福利视频在线观看| 日本精品一区二区三区蜜桃| 人妻少妇偷人精品九色| 看免费成人av毛片| 久久人人爽人人爽人人片va| 欧美高清性xxxxhd video| 亚洲专区中文字幕在线| 国产黄色小视频在线观看| 国产精品久久电影中文字幕| 亚洲电影在线观看av| 一进一出好大好爽视频| 人人妻,人人澡人人爽秒播| eeuss影院久久| 亚洲欧美日韩卡通动漫| 一级黄片播放器| 窝窝影院91人妻| ponron亚洲| 国产主播在线观看一区二区| 夜夜爽天天搞| 中文在线观看免费www的网站| 色综合色国产| 日韩欧美一区二区三区在线观看| 久久久久久久精品吃奶| 校园春色视频在线观看| 伦理电影大哥的女人| 精品久久久久久,| 天天躁日日操中文字幕| 精品久久久久久成人av| 亚洲精品一卡2卡三卡4卡5卡| 日韩中字成人| 性插视频无遮挡在线免费观看| 国产 一区 欧美 日韩| 国内精品美女久久久久久| 亚洲熟妇熟女久久| 日本撒尿小便嘘嘘汇集6| 欧美绝顶高潮抽搐喷水| 国产主播在线观看一区二区| 我要搜黄色片| 国产精品久久久久久久久免| 成人国产综合亚洲| a级毛片免费高清观看在线播放| 精品人妻一区二区三区麻豆 | 国内少妇人妻偷人精品xxx网站| 久久亚洲真实| 老司机深夜福利视频在线观看| 欧美日韩中文字幕国产精品一区二区三区| 91麻豆av在线| 国产免费一级a男人的天堂| 嫩草影院入口| 国产精品乱码一区二三区的特点| 国产一区二区激情短视频| 九色国产91popny在线| 亚洲成人中文字幕在线播放| 日本撒尿小便嘘嘘汇集6| 国内精品美女久久久久久| 国模一区二区三区四区视频| 亚洲最大成人中文| 国产精品女同一区二区软件 | 午夜亚洲福利在线播放| 少妇的逼水好多| 男人舔女人下体高潮全视频| 国模一区二区三区四区视频| 成人一区二区视频在线观看| 亚洲av中文av极速乱 | 国产免费男女视频| 国产主播在线观看一区二区| 春色校园在线视频观看| 99久久精品一区二区三区| 国产男靠女视频免费网站| 91在线精品国自产拍蜜月| av黄色大香蕉| 国产精品综合久久久久久久免费| 极品教师在线免费播放| 亚洲人成网站在线播| 999久久久精品免费观看国产| 国产精品日韩av在线免费观看| 两人在一起打扑克的视频| 婷婷亚洲欧美| 深夜精品福利| 99热这里只有是精品在线观看| 色综合色国产| 国产高清视频在线观看网站| x7x7x7水蜜桃| 老熟妇乱子伦视频在线观看| 国产精品一区二区三区四区久久| 国产在线精品亚洲第一网站| 日韩中字成人| 国产成人aa在线观看| 又黄又爽又免费观看的视频| 色吧在线观看| 欧美中文日本在线观看视频| 亚洲一区二区三区色噜噜| 嫩草影院入口| 国产黄片美女视频| 欧美zozozo另类| 国产精品99久久久久久久久| 综合色av麻豆| 少妇的逼水好多| 国产黄a三级三级三级人| 精品久久久噜噜| 日日撸夜夜添| 乱码一卡2卡4卡精品| 最新中文字幕久久久久| 国产精品久久电影中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 成人二区视频| 91av网一区二区| 最新在线观看一区二区三区| 国产精品无大码| 两人在一起打扑克的视频| 国产三级中文精品| 中文字幕精品亚洲无线码一区| 免费在线观看影片大全网站| 日韩亚洲欧美综合| 日本-黄色视频高清免费观看| 午夜福利视频1000在线观看| 在线a可以看的网站| 亚洲精华国产精华精| 免费无遮挡裸体视频| 欧美zozozo另类| 国内精品一区二区在线观看| 两个人视频免费观看高清| 成人高潮视频无遮挡免费网站| 免费看美女性在线毛片视频| 国产成人一区二区在线| 欧美一区二区亚洲| 在线观看舔阴道视频| 乱人视频在线观看| 在线国产一区二区在线| 最近视频中文字幕2019在线8| 久久久精品大字幕| 五月玫瑰六月丁香| 国产色婷婷99| 日日干狠狠操夜夜爽| 日韩,欧美,国产一区二区三区 | av.在线天堂| 免费观看的影片在线观看| 国产午夜福利久久久久久| 亚洲欧美清纯卡通| 不卡一级毛片| 一本一本综合久久| 三级国产精品欧美在线观看| 麻豆成人午夜福利视频| 一进一出抽搐gif免费好疼| 日本欧美国产在线视频| 国产免费一级a男人的天堂| 国产真实乱freesex| ponron亚洲| 无人区码免费观看不卡| 亚洲真实伦在线观看| 久久久国产成人精品二区| 国产午夜精品久久久久久一区二区三区 | 少妇猛男粗大的猛烈进出视频 | 大又大粗又爽又黄少妇毛片口| АⅤ资源中文在线天堂| 午夜视频国产福利| 最近视频中文字幕2019在线8| 欧美人与善性xxx| 免费看美女性在线毛片视频| 国内精品一区二区在线观看| 日本 欧美在线| 亚洲av二区三区四区| 久久天躁狠狠躁夜夜2o2o| 免费无遮挡裸体视频| 日本 欧美在线| 久久香蕉精品热| 一边摸一边抽搐一进一小说| 亚洲中文字幕日韩| 麻豆成人午夜福利视频| 97碰自拍视频| 最新在线观看一区二区三区| eeuss影院久久| 欧美国产日韩亚洲一区| 我的女老师完整版在线观看| 成人美女网站在线观看视频| 国产成人福利小说| 中文字幕高清在线视频| 深爱激情五月婷婷| 欧美人与善性xxx| 久久久久久久精品吃奶| 女同久久另类99精品国产91| 久久国内精品自在自线图片| 国产精品综合久久久久久久免费| 亚洲人成伊人成综合网2020| 搡老岳熟女国产| 免费黄网站久久成人精品| 国内揄拍国产精品人妻在线| 十八禁国产超污无遮挡网站| 国产亚洲91精品色在线| 小蜜桃在线观看免费完整版高清| 最近视频中文字幕2019在线8| 国产精品99久久久久久久久| 18禁在线播放成人免费| 天天一区二区日本电影三级| 亚洲精品久久国产高清桃花| 麻豆久久精品国产亚洲av| x7x7x7水蜜桃| 色噜噜av男人的天堂激情| 久久精品国产鲁丝片午夜精品 | 亚洲性久久影院| 91在线精品国自产拍蜜月| 国产精品98久久久久久宅男小说| 国产精品av视频在线免费观看| 免费看日本二区| 热99在线观看视频| 淫妇啪啪啪对白视频| 国国产精品蜜臀av免费| 欧美日韩瑟瑟在线播放| 国模一区二区三区四区视频| 久久热精品热| 成人亚洲精品av一区二区| 国产精品久久久久久精品电影| 联通29元200g的流量卡| 舔av片在线| 老熟妇乱子伦视频在线观看| 好男人在线观看高清免费视频| 国产三级中文精品| 蜜桃久久精品国产亚洲av| 精品乱码久久久久久99久播| 成人精品一区二区免费| 丝袜美腿在线中文| 亚洲国产高清在线一区二区三| 国产精品一区www在线观看 | 国产精品一区二区免费欧美| a级毛片a级免费在线| 搞女人的毛片| 国产精华一区二区三区| 国产男靠女视频免费网站| 一卡2卡三卡四卡精品乱码亚洲| 国产精品一及| netflix在线观看网站| 狠狠狠狠99中文字幕| 亚洲中文字幕日韩| 少妇人妻精品综合一区二区 | 麻豆成人av在线观看| 久9热在线精品视频| 国产精品无大码| 天堂影院成人在线观看| 有码 亚洲区| 日韩,欧美,国产一区二区三区 | 色在线成人网| 亚洲男人的天堂狠狠| 亚洲自偷自拍三级| 国产精品99久久久久久久久| 俺也久久电影网| 久久精品人妻少妇| 日韩欧美在线二视频| 国产成人福利小说| 日韩,欧美,国产一区二区三区 | 99久国产av精品| 久久久久久久久大av| 波多野结衣巨乳人妻| 亚洲七黄色美女视频| 国产极品精品免费视频能看的| 简卡轻食公司| 精品久久久久久久久亚洲 | 免费观看人在逋| 国产熟女欧美一区二区| 亚洲成a人片在线一区二区| 婷婷六月久久综合丁香| 午夜视频国产福利| 亚洲五月天丁香| 大型黄色视频在线免费观看| 欧美高清性xxxxhd video| 精品久久久久久久久av| 看免费成人av毛片| 国产真实乱freesex| www日本黄色视频网| 亚洲欧美日韩高清在线视频| 男女之事视频高清在线观看| 色哟哟哟哟哟哟| 91精品国产九色| 国产亚洲91精品色在线| 国内精品一区二区在线观看| 午夜精品久久久久久毛片777| 婷婷色综合大香蕉| 成人一区二区视频在线观看| 丝袜美腿在线中文| 亚洲专区中文字幕在线| 久久精品国产亚洲av天美| 久久精品综合一区二区三区| 亚洲色图av天堂| 国产美女午夜福利| 黄色女人牲交| 在线国产一区二区在线| 久久午夜福利片| 日日啪夜夜撸| 国产又黄又爽又无遮挡在线| 色播亚洲综合网| 成人永久免费在线观看视频| 亚洲七黄色美女视频| 国产一级毛片七仙女欲春2| 亚洲一区高清亚洲精品| 日韩欧美三级三区| 欧美+日韩+精品| 色精品久久人妻99蜜桃| 欧美又色又爽又黄视频| 久久久久久久久久久丰满 | 国产精华一区二区三区| 午夜福利18| 国产大屁股一区二区在线视频| 露出奶头的视频| 麻豆成人av在线观看| 国产日本99.免费观看| 亚洲av.av天堂| 中国美白少妇内射xxxbb| 成人无遮挡网站| 午夜福利在线在线| 国产又黄又爽又无遮挡在线| 超碰av人人做人人爽久久| 亚洲人成网站在线播放欧美日韩| 国产亚洲欧美98| 亚洲精品色激情综合| 一个人看视频在线观看www免费| 日本一二三区视频观看| 国产午夜精品久久久久久一区二区三区 | 成人av在线播放网站| 最新中文字幕久久久久| 中文字幕av在线有码专区| 99国产精品一区二区蜜桃av| 中文字幕av成人在线电影| 久久久久久久久久久丰满 | 99热精品在线国产| 少妇猛男粗大的猛烈进出视频 | 久久精品91蜜桃| 久久草成人影院| 国产一区二区在线av高清观看| 色综合色国产| x7x7x7水蜜桃| 免费不卡的大黄色大毛片视频在线观看 | .国产精品久久| 日本撒尿小便嘘嘘汇集6| 亚洲av中文av极速乱 | 搞女人的毛片| 国产精品自产拍在线观看55亚洲| 国产av麻豆久久久久久久| 99久久中文字幕三级久久日本| 如何舔出高潮| 神马国产精品三级电影在线观看| 国产精品久久视频播放| 色综合色国产| 日韩av在线大香蕉| 亚洲精品一卡2卡三卡4卡5卡| 日本 欧美在线| 国产高清不卡午夜福利| 亚洲va在线va天堂va国产| 尾随美女入室| 亚洲av免费在线观看| 精品午夜福利视频在线观看一区| 国产三级在线视频| 日韩av在线大香蕉| 日韩一区二区视频免费看| 亚州av有码| 又爽又黄a免费视频| 搡女人真爽免费视频火全软件 | 成人特级av手机在线观看| netflix在线观看网站| 一级毛片久久久久久久久女| 伦精品一区二区三区| 日韩强制内射视频| 搡老熟女国产l中国老女人| 欧美性感艳星| 一个人看视频在线观看www免费| 嫩草影院新地址|