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

    基于Benders分解的獨(dú)立型微電網(wǎng)魯棒優(yōu)化容量配置模型

    2017-12-22 10:02:37劉思夷汪湘晉邱海峰林湘寧
    電力系統(tǒng)自動(dòng)化 2017年21期
    關(guān)鍵詞:臺(tái)數(shù)魯棒儲(chǔ)能

    劉思夷, 趙 波, 汪湘晉, 邱海峰, 林湘寧

    (1. 強(qiáng)電磁工程與新技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室(華中科技大學(xué)), 湖北省武漢市 430074; 2. 國(guó)網(wǎng)浙江省電力公司電力科學(xué)研究院, 浙江省杭州市 310014)

    基于Benders分解的獨(dú)立型微電網(wǎng)魯棒優(yōu)化容量配置模型

    劉思夷1, 趙 波2, 汪湘晉2, 邱海峰2, 林湘寧1

    (1. 強(qiáng)電磁工程與新技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室(華中科技大學(xué)), 湖北省武漢市 430074; 2. 國(guó)網(wǎng)浙江省電力公司電力科學(xué)研究院, 浙江省杭州市 310014)

    為了提高獨(dú)立型微電網(wǎng)電源優(yōu)化配置的魯棒性,綜合考慮了風(fēng)、光和負(fù)荷的間歇性和隨機(jī)性,建立了獨(dú)立型微電網(wǎng)的雙層魯棒優(yōu)化容量配置模型。外層為初始投資費(fèi)用主問(wèn)題,內(nèi)層為基于風(fēng)、光和負(fù)荷不確定集的運(yùn)行優(yōu)化子問(wèn)題。采用對(duì)偶的方法對(duì)子問(wèn)題進(jìn)行解耦,利用大M法將其線性化,并通過(guò)Benders分解的方法將主子問(wèn)題進(jìn)行交互迭代,直至經(jīng)濟(jì)成本達(dá)到最優(yōu)。在算例中將該方法與改進(jìn)的非支配排序遺傳算法(NSGA-Ⅱ)進(jìn)行對(duì)比,結(jié)果表明所提模型的魯棒性更好,求解效率更高,驗(yàn)證了模型的合理性和有效性。

    微電網(wǎng)(微網(wǎng)); 優(yōu)化配置; 魯棒優(yōu)化; Benders分解; 不確定集

    0 引言

    微電網(wǎng)由于能夠靈活接入分布式電源,其應(yīng)用得到了廣泛關(guān)注。目前在一些偏遠(yuǎn)山區(qū)和遠(yuǎn)離城市的海島,由于很難建立與大電網(wǎng)的互聯(lián),在這些地區(qū)建立供電可靠的獨(dú)立型微電網(wǎng)是有效利用當(dāng)?shù)乜稍偕Y源、解決電力不足、減少環(huán)境污染的有效方法[1]。

    分布式電源的優(yōu)化配置是微電網(wǎng)規(guī)劃的主要前期工作,其涉及電源容量配置和運(yùn)行優(yōu)化的耦合性,解決該類問(wèn)題的一種途徑是將電源容量和運(yùn)行變量聯(lián)合建模優(yōu)化。文獻(xiàn)[2]采用枚舉法聯(lián)立求解所有可能機(jī)組臺(tái)數(shù)和運(yùn)行變量的組合方案,但存在解空間大、求解時(shí)間較長(zhǎng)的缺點(diǎn)。另一種途徑是將問(wèn)題分為兩層交互求解,外層為容量配置問(wèn)題,內(nèi)層為運(yùn)行優(yōu)化問(wèn)題。文獻(xiàn)[3]基于粒子群算法,以電網(wǎng)等年值費(fèi)用最低為目標(biāo),對(duì)風(fēng)/光/儲(chǔ)獨(dú)立型供電系統(tǒng)優(yōu)化配置問(wèn)題進(jìn)行了求解;文獻(xiàn)[4]采用改進(jìn)細(xì)菌覓食算法以系統(tǒng)的經(jīng)濟(jì)性為目標(biāo)函數(shù),求解了風(fēng)/光/儲(chǔ)獨(dú)立型微電網(wǎng)的優(yōu)化配置問(wèn)題。

    以上方法都是基于智能算法的啟發(fā)式隨機(jī)搜索算法,該類算法通常需要大量初始種群,不能保證找到全局最優(yōu)解。同時(shí),以上方法均沒(méi)有充分考慮容量配置和運(yùn)行優(yōu)化的耦合性,不能較好體現(xiàn)運(yùn)行優(yōu)化對(duì)容量配置的影響。

    此外,由于風(fēng)、光資源和負(fù)荷具有較強(qiáng)的間歇性和波動(dòng)性,會(huì)影響微電網(wǎng)優(yōu)化配置的經(jīng)濟(jì)性和可靠性,因此應(yīng)充分考慮其不確定性因素的影響。減小不確定性因素影響的方法有引入儲(chǔ)能系統(tǒng)和負(fù)荷側(cè)響應(yīng)。文獻(xiàn)[5]研究了引入電動(dòng)汽車和水泵儲(chǔ)水系統(tǒng)以應(yīng)對(duì)獨(dú)立儲(chǔ)能系統(tǒng)的可再生能源的波動(dòng)性;文獻(xiàn)[6]考慮引入負(fù)荷側(cè)響應(yīng)來(lái)減小大規(guī)模風(fēng)電并網(wǎng)不確定性的影響。針對(duì)微電網(wǎng)中不確定因素影響的研究方法主要有如下三種。

    1)基于歷史數(shù)據(jù)的序貫仿真方法。文獻(xiàn)[7]利用確定的數(shù)據(jù)序貫仿真模擬風(fēng)、光出力,建立了考慮電網(wǎng)經(jīng)濟(jì)性和可靠性的微電網(wǎng)優(yōu)化配置模型。但基于歷史數(shù)據(jù)的序貫仿真方法不能很好地體現(xiàn)不確定性對(duì)系統(tǒng)魯棒性的影響。

    2)基于概率論的隨機(jī)規(guī)劃方法。文獻(xiàn)[8]考慮系統(tǒng)的不確定因素,采用隨機(jī)機(jī)會(huì)約束規(guī)劃方法。文獻(xiàn)[9]采用拉丁超立方抽樣生成風(fēng)電場(chǎng)景并進(jìn)行縮減?;诟怕收摰碾S機(jī)規(guī)劃方法存在分布函數(shù)不能準(zhǔn)確擬合的缺點(diǎn),并且需要大量抽樣以獲得隨機(jī)場(chǎng)景并進(jìn)行場(chǎng)景縮減,增加了問(wèn)題的復(fù)雜性。

    3)魯棒優(yōu)化方法。文獻(xiàn)[10]采用盒式集合表示風(fēng)速的不確定性對(duì)并網(wǎng)型微電網(wǎng)的風(fēng)機(jī)容量進(jìn)行規(guī)劃;文獻(xiàn)[11]將魯棒優(yōu)化理論引入到電力系統(tǒng)優(yōu)化調(diào)度中,構(gòu)建了大規(guī)模間歇式電源接入的多目標(biāo)魯棒優(yōu)化調(diào)度模型。由于魯棒優(yōu)化不需要獲取不確定數(shù)據(jù)的概率分布,而是采用集合的方式表示不確定因素的取值范圍[12]。因此采用魯棒優(yōu)化方法構(gòu)建目標(biāo)函數(shù),有利于增強(qiáng)系統(tǒng)應(yīng)對(duì)最惡劣自然條件時(shí)的性能,提高系統(tǒng)穩(wěn)定性。

    本文以包含風(fēng)/光/柴/儲(chǔ)的獨(dú)立型微電網(wǎng)為研究對(duì)象,構(gòu)建了考慮風(fēng)、光出力和負(fù)荷不確定性的二層魯棒優(yōu)化混合整數(shù)規(guī)劃模型。目標(biāo)函數(shù)的外層為微電網(wǎng)的初始投資費(fèi)用,內(nèi)層為包含不確定集的微電網(wǎng)的運(yùn)行費(fèi)用、停電懲罰費(fèi)用和棄電費(fèi)用;采用Lagrange對(duì)偶方法將內(nèi)層最小—最大魯棒優(yōu)化模型進(jìn)行解耦,引入大M法將對(duì)偶問(wèn)題線性化;用Benders分解的方法將外層問(wèn)題和內(nèi)層問(wèn)題交互迭代求解直到目標(biāo)函數(shù)達(dá)到最優(yōu)。在算例分析中,將該模型與改進(jìn)的非支配排序遺傳算法(NSGA-Ⅱ)進(jìn)行對(duì)比,驗(yàn)證了所提模型的合理性。

    1 獨(dú)立型交流微電網(wǎng)模型

    本文采用的獨(dú)立型微電網(wǎng)的拓?fù)浣Y(jié)構(gòu)如圖1所示,其中蓄電池和光伏通過(guò)換流設(shè)備連接在交流母線上。微電網(wǎng)采用孤島運(yùn)行方式,依靠?jī)?nèi)部電源滿足負(fù)荷供電需求。當(dāng)風(fēng)機(jī)和光伏發(fā)電過(guò)多或負(fù)荷較小時(shí),微電網(wǎng)向儲(chǔ)能充電,若儲(chǔ)能容量達(dá)到上限則丟棄部分多余發(fā)電功率。當(dāng)風(fēng)電和光伏發(fā)電不足時(shí),儲(chǔ)能向負(fù)荷放電或者開啟柴油發(fā)電機(jī)向負(fù)荷供電。若仍不能滿足負(fù)荷需求,則切除部分負(fù)荷。

    圖1 獨(dú)立型交流微電網(wǎng)結(jié)構(gòu)Fig.1 Structure of stand-alone alternating current microgrid

    2 微電源模型

    2.1 風(fēng)力發(fā)電模型

    風(fēng)機(jī)的發(fā)電功率由風(fēng)機(jī)轂輪高度處的風(fēng)速?zèng)Q定,用指數(shù)定律表達(dá)式對(duì)風(fēng)速進(jìn)行變換[13]。

    (1)

    式中:α為地表粗糙程度的描述因素,本文取值為0.22;vref(t)為測(cè)風(fēng)點(diǎn)風(fēng)速;Href為測(cè)風(fēng)點(diǎn)高度;H為風(fēng)機(jī)轂輪高度。

    風(fēng)機(jī)的輸出特性曲線采用變漿距風(fēng)機(jī)輸出特性曲線[14],當(dāng)風(fēng)速小于切入風(fēng)速vci時(shí),風(fēng)機(jī)發(fā)電功率為0;當(dāng)風(fēng)速大于切入風(fēng)速vci小于額定風(fēng)速vcr時(shí),其輸出功率Pr(v)由變漿距風(fēng)機(jī)曲線擬合得到;當(dāng)風(fēng)速達(dá)到額定風(fēng)速vcr小于切出風(fēng)速vco時(shí),風(fēng)機(jī)以額定功率Pr運(yùn)行;當(dāng)風(fēng)速大于切出風(fēng)速vco時(shí),風(fēng)機(jī)停止運(yùn)行,輸出功率為0。

    (2)

    2.2 光伏發(fā)電模型

    本文采用的光伏發(fā)電的輸出功率模型為[15]:

    (3)

    式中:fpv為光伏陣列能量轉(zhuǎn)換效率,本文取0.9;Prate為光伏陣列的額定功率;A為光伏的實(shí)際輻照度;As為光伏的額定輻照度;αp為功率溫度系數(shù),一般取-0.004 7 ℃-1;TSTC為標(biāo)準(zhǔn)測(cè)試條件下光伏組件的溫度,一般取25 ℃;Tpv為光伏組件的實(shí)際溫度。

    2.3 柴油發(fā)電機(jī)模型

    (4)

    柴油發(fā)電機(jī)的耗油量V(t)與運(yùn)行功率之間的關(guān)系為[16]:

    (5)

    式中:F0為柴油發(fā)電機(jī)燃料曲線的截距系數(shù),即柴油機(jī)的空載油耗量;F1為柴油發(fā)電機(jī)燃料曲線的斜率。

    2.4 儲(chǔ)能模型

    儲(chǔ)能具有最大充放電功率限制,其最大充放電功率與電池額定容量Erate成正比,儲(chǔ)能實(shí)際充放電功率Pbat(t)為正表示儲(chǔ)能放電,若為負(fù)則表示儲(chǔ)能充電:

    -εErate≤Pbat(t)≤εErate

    (6)

    式中:ε為儲(chǔ)能最大充放電功率與儲(chǔ)能額定容量的比值。

    儲(chǔ)能的荷電狀態(tài)與充放電功率的關(guān)系為:

    (7)

    式中:SSOC為儲(chǔ)能的荷電狀態(tài);Pdch(t)為儲(chǔ)能的放電功率;Δt為時(shí)間步長(zhǎng),本文取1 h;ηd為電能轉(zhuǎn)化為放電功率的效率;Pch(t)為儲(chǔ)能充電功率;ηc為儲(chǔ)能的充電效率。

    儲(chǔ)能電池的荷電狀態(tài)具有上下限,若荷電狀態(tài)太低或太高都會(huì)縮短儲(chǔ)能電池的使用壽命:

    (8)

    2.5 風(fēng)、光出力及負(fù)荷的不確定集模型

    由于在獨(dú)立型微電網(wǎng)中風(fēng)機(jī)和光伏發(fā)電的比例較大,而風(fēng)速、光照強(qiáng)度和負(fù)荷通常具有間歇性、波動(dòng)性和不可預(yù)測(cè)性,風(fēng)機(jī)、光伏發(fā)電和負(fù)荷功率的大小會(huì)影響系統(tǒng)的運(yùn)行狀態(tài),因此需考慮不確定因素對(duì)系統(tǒng)的影響。魯棒優(yōu)化采用集合的方式表示變量的不確定性,其關(guān)鍵在于不確定集的設(shè)計(jì)。根據(jù)長(zhǎng)期統(tǒng)計(jì)信息,風(fēng)機(jī)、光伏的輸出功率和負(fù)荷的大小可以表示成在預(yù)測(cè)值附近上下浮動(dòng),浮動(dòng)量由系統(tǒng)的不確定性決定。整個(gè)周期內(nèi)擾動(dòng)量的取值具有上下限值。不確定集的表示形式如式(9)所示。

    (9)

    3 優(yōu)化配置模型的建立和求解算法

    3.1 優(yōu)化配置的目標(biāo)函數(shù)

    由于微電網(wǎng)電源容量配置在很大程度上受到運(yùn)行優(yōu)化的影響,在本質(zhì)上是雙層優(yōu)化,并且要考慮風(fēng)、光、負(fù)荷等不確定因素對(duì)系統(tǒng)穩(wěn)定性的作用,因此可建立雙層魯棒優(yōu)化容量配置模型,以提高系統(tǒng)的魯棒性。雙層魯棒優(yōu)化原問(wèn)題的目標(biāo)函數(shù)分為兩層:外層為分布式電源的最小凈現(xiàn)值投資成本,內(nèi)層為包含風(fēng)、光和負(fù)荷不確定集的運(yùn)行費(fèi)用、停電懲罰費(fèi)用和棄電費(fèi)用。內(nèi)層運(yùn)行成本的目標(biāo)函數(shù)為在所有不確定集作用下使得目標(biāo)函數(shù)盡可能大的情況下,在運(yùn)行變量作用下使目標(biāo)函數(shù)盡可能小,在數(shù)學(xué)上可表示帶有約束的最小—最大魯棒優(yōu)化模型[17]。

    Ctotal=minCint+min(maxCopt)

    (10)

    (11)

    Copt=Cde+Clsh+Cdump=

    (12)

    3.2 廣義Benders分解的基本原理

    以上優(yōu)化問(wèn)題是一個(gè)典型的混合整數(shù)規(guī)劃問(wèn)題,適用于Benders分解的方法求解。Benders分解最開始用于求解混合整數(shù)規(guī)劃問(wèn)題,即實(shí)型變量和整型變量同時(shí)出現(xiàn)時(shí)的極值問(wèn)題[18],其基本原理見附錄A。

    3.3 Benders分解流程

    利用Benders分解原理求解微電網(wǎng)優(yōu)化配置問(wèn)題的基本原理見附錄B圖B1,求解步驟見附錄B圖B2。

    求解過(guò)程中將原問(wèn)題分解為主問(wèn)題和子問(wèn)題。主問(wèn)題是以機(jī)組臺(tái)數(shù)為決策變量的電源容量規(guī)劃問(wèn)題,變量類型為整型變量。子問(wèn)題是以各機(jī)組出力、棄電量和切負(fù)荷量為決策變量的最優(yōu)運(yùn)行子問(wèn)題,包含風(fēng)、光出力和負(fù)荷的不確定集,變量類型為整型變量和實(shí)型變量。主問(wèn)題初始化機(jī)組臺(tái)數(shù)代入子問(wèn)題,子問(wèn)題首先進(jìn)行可行性檢驗(yàn),判斷主問(wèn)題的決策變量是否能滿足所有約束條件。若不能滿足則返回不可行割,對(duì)主問(wèn)題優(yōu)化臺(tái)數(shù)進(jìn)行調(diào)整,直到滿足所有約束條件為止。若主問(wèn)題的決策變量能夠滿足所有約束條件,子問(wèn)題進(jìn)一步進(jìn)行經(jīng)濟(jì)優(yōu)化計(jì)算,向主問(wèn)題返回該臺(tái)數(shù)配置附近的邊際增量信息(即對(duì)偶變量的大小),對(duì)主問(wèn)題中的臺(tái)數(shù)配置進(jìn)行調(diào)整,直到找到最優(yōu)解。采用Benders分解方法求解電源容量?jī)?yōu)化配置問(wèn)題的具體步驟如下。

    步驟1:初始化原問(wèn)題的上界和下界。

    步驟2:求解機(jī)組臺(tái)數(shù)配置主問(wèn)題,并更新原問(wèn)題的下界。

    步驟3:將主問(wèn)題中求解得到的機(jī)組臺(tái)數(shù)代入可行性檢驗(yàn)子問(wèn)題,判斷該機(jī)組臺(tái)數(shù)配置是否能夠滿足所有約束。若可行性檢驗(yàn)問(wèn)題中引入的松弛變量之和為零,說(shuō)明問(wèn)題可行,進(jìn)行下一步。若松弛變量之和不為零,說(shuō)明該機(jī)組臺(tái)數(shù)配置不滿足所有約束條件。則給主問(wèn)題添加不可行割,返回步驟2重新求解主問(wèn)題,對(duì)主問(wèn)題的機(jī)組臺(tái)數(shù)進(jìn)行調(diào)整。

    步驟4:若主問(wèn)題初始化臺(tái)數(shù)滿足所有約束,則計(jì)算運(yùn)行優(yōu)化子問(wèn)題的對(duì)偶問(wèn)題,更新原問(wèn)題的上界。

    步驟5:若上界與下界的差小于閾值則說(shuō)明原問(wèn)題已經(jīng)達(dá)到最優(yōu)解。否則給主問(wèn)題增加可行割,重新求解主問(wèn)題,對(duì)機(jī)組臺(tái)數(shù)進(jìn)行調(diào)整。

    子問(wèn)題和主問(wèn)題的數(shù)學(xué)形式如下。

    1)可行性檢驗(yàn)子問(wèn)題的目標(biāo)函數(shù)是為了滿足所有約束條件引入的松弛變量之和最小,其在第w次迭代時(shí)的求解形式:

    min(maxS)=S1+S2+S3

    (13)

    (14)

    (15)

    (16)

    (17)

    (18)

    (19)

    (20)

    (21)

    Plsh(t)≥0

    (22)

    Pdump(t)≥0

    (23)

    (24)

    (25)

    (26)

    (27)

    S1≥0

    (28)

    S2≥0

    (29)

    S3≥0

    (30)

    式(14)表示功率平衡約束;式(15)表示柴油發(fā)電機(jī)組的出力范圍約束;式(16)為可再生能源機(jī)組的發(fā)電功率;式(17)至式(19)為儲(chǔ)能約束,由于初始時(shí)刻儲(chǔ)能的容量E0為常數(shù)而不是變量,所以單獨(dú)列出;式(22)表示切負(fù)荷量大于等于零;式(23)表示棄電量大于等于零;式(24)表示全年切負(fù)荷不超過(guò)總負(fù)荷的比例,取值為μL;式(25)表示棄電量不超過(guò)總負(fù)荷的λ倍,否則會(huì)造成可再生能源的浪費(fèi);式(26)表示可再生能源功率滲透率要達(dá)到總負(fù)荷ζ倍;式(27)表明在一個(gè)運(yùn)行周期后儲(chǔ)能電池的容量保持在初值;式(28)至式(30)表示引入的對(duì)偶變量大于等于零。

    2)經(jīng)濟(jì)最優(yōu)運(yùn)行子問(wèn)題

    經(jīng)濟(jì)最優(yōu)運(yùn)行子問(wèn)題的目標(biāo)函數(shù)為最大不確定性下的運(yùn)行費(fèi)用最小:

    (31)

    經(jīng)濟(jì)最優(yōu)子問(wèn)題的約束條件與可行性檢驗(yàn)子問(wèn)題相同,只是約束條件(式(24)至式(26))中不含有松弛變量。采用Lagrange對(duì)偶方法對(duì)最小—最大模型進(jìn)行解耦,其對(duì)偶問(wèn)題的目標(biāo)函數(shù)為:

    (32)

    式中:yi為序號(hào)為i的約束條件引入的對(duì)偶變量。

    原問(wèn)題的上界為:

    (33)

    式中:Λ為運(yùn)行優(yōu)化子問(wèn)題的目標(biāo)函數(shù)值。

    3)主問(wèn)題第w次迭代時(shí)的求解形式為:

    (34)

    可行性檢驗(yàn)子問(wèn)題增加的不可行割為:

    (35)

    最優(yōu)運(yùn)行子問(wèn)題增加的可行割為:

    (36)

    返回的可行割通過(guò)對(duì)偶變量來(lái)表示原問(wèn)題的邊際增量,調(diào)整后的機(jī)組臺(tái)數(shù)決策變量能夠得到更經(jīng)濟(jì)的運(yùn)行子問(wèn)題,給定原問(wèn)題的下界為Z。

    在求解子問(wèn)題對(duì)偶問(wèn)題的過(guò)程中,存在對(duì)偶變量與不確定集中0-1變量的乘積,采用大M法[19]將其線性化進(jìn)行求解。

    4 算例分析

    本文采用實(shí)際風(fēng)、光資源和負(fù)荷數(shù)據(jù)進(jìn)行計(jì)算,其中風(fēng)速、光照幅度和負(fù)荷曲線見附錄B圖B3,全年平均風(fēng)速為5.70 m/s,負(fù)荷均值為73.16 kW,日平均輻照為3.07(kW·h)/m2。

    系統(tǒng)的具體參數(shù)如下:全年切負(fù)荷率小于0.3%,棄電率小于分布式電源總發(fā)電量的10%,可再生能源功率滲透率大于等于70%。儲(chǔ)能容量荷電狀態(tài)初值為0.7,荷電狀態(tài)的設(shè)定范圍為[0.4,0.95],儲(chǔ)能充放電效率為0.95。整個(gè)工程的壽命為20年,通貨膨脹率為0.03,棄電費(fèi)用為0.6元/kW,切負(fù)荷費(fèi)用為1.3元/kW。分布式電源和儲(chǔ)能的參數(shù)見附錄B表B1和表B2。

    4.1 不同不確定水平下優(yōu)化配置結(jié)果的比較

    對(duì)全年每個(gè)月構(gòu)造典型日,采用每個(gè)月各時(shí)刻的平均值作為該月典型日預(yù)測(cè)值,即不考慮不確定性。此時(shí),風(fēng)、光發(fā)電及負(fù)荷數(shù)據(jù)采用每月預(yù)測(cè)值時(shí)的機(jī)組臺(tái)數(shù)優(yōu)化配置結(jié)果詳見附錄B表B3場(chǎng)景B1,經(jīng)濟(jì)成本見表1場(chǎng)景B1。設(shè)置不確定集的波動(dòng)范圍為預(yù)測(cè)值的±20%,改變不確定集取端點(diǎn)值個(gè)數(shù)的多少,得到機(jī)組配置臺(tái)數(shù)詳見附錄B表B3場(chǎng)景B2和B3,得到系統(tǒng)成本見表1場(chǎng)景B2和B3。β為在12個(gè)典型日的計(jì)算周期內(nèi),不確定集取上界(下界)的個(gè)數(shù),當(dāng)其為0,表示不考慮不確定性。

    表1 不同不確定性下經(jīng)濟(jì)成本比較Table 1 Comparison of optimal cost under various uncertainty

    從表1可以看出,將只采用柴油發(fā)電機(jī)供電模式下各場(chǎng)景的經(jīng)濟(jì)成本與采用微電網(wǎng)供電模式下各場(chǎng)景的經(jīng)濟(jì)成本比較可知,只采用柴油組發(fā)電機(jī)發(fā)電,需配置5臺(tái)柴油發(fā)電機(jī),全年的總成本高于含風(fēng)機(jī)和光伏發(fā)電的微電網(wǎng)總成本。說(shuō)明雖然可再生能源初始投資成本較高,但由于運(yùn)行費(fèi)用較小,采用風(fēng)機(jī)、光伏等可再生能源機(jī)組不僅具有環(huán)保性還具有經(jīng)濟(jì)性。

    從機(jī)組臺(tái)數(shù)的優(yōu)化配置結(jié)果可以看出,隨著系統(tǒng)的不確定性增加,微電網(wǎng)系統(tǒng)會(huì)減小可再生能源機(jī)組配置,增加柴油發(fā)電機(jī)組的配置。原因是隨著系統(tǒng)的不確定性增加,棄電率增加,為保證棄電率不超過(guò)10%,減小可再生能源出力而增大柴油機(jī)的發(fā)電量,以提高系統(tǒng)的魯棒性。從表1中的經(jīng)濟(jì)成本可以看出,當(dāng)風(fēng)、光和負(fù)荷的不確定性增大,無(wú)論是柴油發(fā)電機(jī)供電系統(tǒng)還是微電網(wǎng)系統(tǒng),系統(tǒng)的運(yùn)行費(fèi)用都隨之增大。

    在不考慮不確定性的情況下,確定優(yōu)化配置的機(jī)組臺(tái)數(shù)后,以運(yùn)行成本最小為目標(biāo)函數(shù)優(yōu)化各運(yùn)行變量,給出在2月、6月、9月和12月4個(gè)月份典型日內(nèi)的各機(jī)組出力曲線,如圖2所示。圖中,Pload為負(fù)荷功率,Pwt為風(fēng)機(jī)的發(fā)電功率,Ppv為光伏發(fā)電功率。

    從圖2可以看出,在2月和12月,風(fēng)力發(fā)電充足,會(huì)產(chǎn)生較多的棄電量,柴油發(fā)電機(jī)處于最小出力。而在6月和9月,全天負(fù)荷增大,風(fēng)力發(fā)電的減小,光伏發(fā)電增大,在白天,光伏發(fā)電較好地彌補(bǔ)了風(fēng)力發(fā)電的缺額,減小了柴油發(fā)電機(jī)的開啟和儲(chǔ)能的充放,說(shuō)明風(fēng)、光具有良好的互補(bǔ)性;在下午和傍晚,開啟柴油發(fā)電彌補(bǔ)由于風(fēng)機(jī)和光伏發(fā)電量減小引起的負(fù)荷缺額。

    圖2 不考慮不確定性時(shí)四季典型日各機(jī)組出力Fig.2 Output power of each unit in classic days of four seasons without considering uncertainty

    在求解過(guò)程中發(fā)現(xiàn),單位臺(tái)數(shù)的風(fēng)機(jī)經(jīng)濟(jì)性高于光伏發(fā)電的經(jīng)濟(jì)性。雖然光伏發(fā)電和風(fēng)電的單位功率運(yùn)行成本相同,求解算出的單位功率的邊際增量(對(duì)偶變量)也相同,但由于風(fēng)機(jī)全年功率的最大值為13.61 kW,最小值為2.01 kW,而光伏發(fā)電在早晨和晚間功率出力為0 kW,全年最大發(fā)電功率為11.18 kW,所以單位臺(tái)數(shù)風(fēng)機(jī)發(fā)電的總功率更大,邊際增量更大,導(dǎo)致采用Benders分解方法配置的風(fēng)機(jī)臺(tái)數(shù)多于光伏的臺(tái)數(shù),經(jīng)濟(jì)性更好。

    當(dāng)系統(tǒng)的不確定性如場(chǎng)景B3所示時(shí),確定機(jī)組臺(tái)數(shù)后,以運(yùn)行成本最小為目標(biāo)函數(shù),對(duì)各運(yùn)行變量進(jìn)行優(yōu)化,給出在2月、6月、9月和12月4個(gè)月份的典型日內(nèi)的各機(jī)組出力曲線,如圖3所示。

    圖3 考慮不確定性時(shí)四季典型日各機(jī)組出力Fig.3 Output power of each unit in classic days of four seasons considering uncertainty

    將圖3各時(shí)刻的功率與圖2各時(shí)刻功率進(jìn)行比較,從圖3可以看出,雖然減少了風(fēng)機(jī)的配置臺(tái)數(shù),但在2月和12月,由于風(fēng)、光及負(fù)荷不確定性的增加,負(fù)荷和風(fēng)、光出力的波動(dòng)更大,在09:00至17:00時(shí)仍有較多的棄電量。在6月和9月,對(duì)比柴油發(fā)電機(jī)09:00至14:00的發(fā)電量,發(fā)現(xiàn)柴油發(fā)電機(jī)發(fā)電量有所增加,發(fā)電功率的波動(dòng)性也更大。

    4.2 采用NSGA-Ⅱ的優(yōu)化配置結(jié)果

    采用NSGA-Ⅱ,取棄電量和系統(tǒng)總成本為目標(biāo)函數(shù),不確定集波動(dòng)范圍為20%,不確定集取邊界值的個(gè)數(shù)同場(chǎng)景B3,取各機(jī)組臺(tái)數(shù)為優(yōu)化變量,種群個(gè)數(shù)為50,進(jìn)化代數(shù)為20,迭代收斂后得到8組Pareto最優(yōu)解,優(yōu)化配置結(jié)果詳見附錄B表B4,經(jīng)濟(jì)成本見表2。

    表2 NSGA-Ⅱ優(yōu)化配置成本Table 2 Optimal configuration result of NSGA-Ⅱ

    從機(jī)組臺(tái)數(shù)的優(yōu)化配置結(jié)果可以看出,柴油機(jī)、光伏、風(fēng)機(jī)的臺(tái)數(shù)相對(duì)固定,隨著儲(chǔ)能臺(tái)數(shù)的增加,系統(tǒng)投資成本逐漸增加,運(yùn)行成本逐漸降低,系統(tǒng)的總成本逐漸增加,棄電量逐漸減小。

    比較分析配置N4和配置N5,當(dāng)配置的儲(chǔ)能和柴油發(fā)電機(jī)臺(tái)數(shù)相同,可再生能源機(jī)組的總臺(tái)數(shù)也相同,投資成本相同,但是配置N4的棄電量為35 943 kW·h,高于配置N5 的棄電量12 217 kW·h;配置N4的運(yùn)行成本為15.16萬(wàn)元,低于配置N5的運(yùn)行成本21.55萬(wàn)元??芍獑挝慌_(tái)數(shù)風(fēng)機(jī)發(fā)電比光伏發(fā)電更加經(jīng)濟(jì),但棄電量有所增加。原因是單位臺(tái)數(shù)的風(fēng)機(jī)發(fā)電量高于光伏的發(fā)電量。當(dāng)配置的風(fēng)機(jī)臺(tái)數(shù)比例較高時(shí),可再生能源的發(fā)電量更高,柴油發(fā)電機(jī)的運(yùn)行時(shí)間更少,因此經(jīng)濟(jì)性更好,但是棄電率較高。從配置N1和配置N5對(duì)比也可發(fā)現(xiàn),要達(dá)到相同的可再生能源功率滲透率,配置風(fēng)機(jī)臺(tái)數(shù)比例較高時(shí),總可再生能源機(jī)組臺(tái)數(shù)較少,投資成本也較小,更具經(jīng)濟(jì)性,與Benders分解方法的配置結(jié)果一致。

    當(dāng)系統(tǒng)具有相同的不確定性時(shí),將NSGA-Ⅱ求解結(jié)果中最具經(jīng)濟(jì)性的方案配置N1與本文提出的二層魯棒優(yōu)化模型場(chǎng)景B3中的配置結(jié)果對(duì)比,可知兩種配置儲(chǔ)能臺(tái)數(shù)相同,可再生能源和柴油發(fā)電機(jī)的總?cè)萘肯嗤?都為390 kW,說(shuō)明基于本文模型算法的配置結(jié)果具有合理性。通過(guò)對(duì)比發(fā)現(xiàn)在系統(tǒng)的不確定性較大時(shí),采用Benders分解算法配置結(jié)果會(huì)減小可再生能源機(jī)組的配置,增大柴油發(fā)電機(jī)機(jī)組的配置以提高系統(tǒng)的魯棒性。

    從經(jīng)濟(jì)成本上來(lái)看,采用NSGA-Ⅱ配置的運(yùn)行成本為16.23萬(wàn)元,小于Benders分解算法的運(yùn)行成本30.91萬(wàn)元。從棄電量上來(lái)看,采用Benders分解方法的棄電量為51 739 kW·h,小于NSGA-Ⅱ的N1配置的棄電量53 452 kW·h。采用Benders分解方法的配置由于減少了可再生能源機(jī)組的臺(tái)數(shù),系統(tǒng)配置的棄電量有所下降。將配置N1機(jī)組臺(tái)數(shù)帶入Benders分解方法中的可行性檢驗(yàn)子問(wèn)題中,得到松弛變量的值為226.56,說(shuō)明該方法不滿足在二層魯棒優(yōu)化模型中不確定集作用下的所有約束條件。原因是本文構(gòu)建了最小—最大模型的運(yùn)行子問(wèn)題,并采用Lagrange對(duì)偶的方法對(duì)問(wèn)題進(jìn)行了解耦,求解的結(jié)果為最大不確定性情況下的運(yùn)行成本,而采用NSGA-Ⅱ求解的運(yùn)行成本是在不確定集作用下的最小運(yùn)行成本,所以采用本文提出的二層魯棒優(yōu)化模型的運(yùn)行費(fèi)用更大,柴油發(fā)電機(jī)組臺(tái)數(shù)更多。因此,本文提出的模型能較好地體現(xiàn)系統(tǒng)的魯棒性,能提高系統(tǒng)的抗風(fēng)險(xiǎn)能力,對(duì)運(yùn)行子問(wèn)題的建模也更符合實(shí)際。

    仿真在CPU為2.30 GHz的IBM-PC,MATLAB 2012a的環(huán)境下運(yùn)行,優(yōu)化通過(guò)調(diào)用CPLEX和YALMIP優(yōu)化工具箱完成。從求解時(shí)間上看,采用Benders分解求解模型經(jīng)過(guò)11次迭代即可獲得結(jié)果,所需時(shí)間為106 s,而NSGA-Ⅱ迭代需要求解1 000次優(yōu)化才能獲得結(jié)果,所需時(shí)間為444 min。采用Benders分解求解方法采用對(duì)偶變量返回外層決策變量的邊際增量信息,充分體現(xiàn)了內(nèi)層運(yùn)行問(wèn)題對(duì)外層投資問(wèn)題的影響,提高了模型的求解效率。

    5 結(jié)語(yǔ)

    本文以獨(dú)立型交流微電網(wǎng)的投資費(fèi)用、運(yùn)行費(fèi)用、棄電費(fèi)用和停電費(fèi)用之和最小為目標(biāo)函數(shù),構(gòu)建了包含風(fēng)、光和負(fù)荷不確定集的二層魯棒優(yōu)化模型;采用最小—最大模型表征表示不確定性對(duì)內(nèi)層運(yùn)行子問(wèn)題的影響,用Benders分解方法對(duì)模型內(nèi)外層進(jìn)行迭代求解,提高了模型的求解效率。在算例分析中,將基于Benders分解方法的配置結(jié)果與采用NSGA-Ⅱ的配置結(jié)果進(jìn)行對(duì)比,采用本文提出的模型能根據(jù)風(fēng)、光及負(fù)荷不確定性的大小提升微電網(wǎng)的魯棒性,增強(qiáng)系統(tǒng)抵抗風(fēng)險(xiǎn)的能力,證明了模型的有效性和合理性。

    由于風(fēng)、光、負(fù)荷偏離預(yù)測(cè)值的不確定性隨著時(shí)間的推移有可能發(fā)生變化,在本文研究的基礎(chǔ)上,后續(xù)還需要計(jì)及該因素影響,研究更準(zhǔn)確的機(jī)組容量配置方法。

    本文的研究工作受到國(guó)網(wǎng)浙江省電力公司科技項(xiàng)目(5211DS15002A)資助,謹(jǐn)此致謝。

    附錄見本刊網(wǎng)絡(luò)版(http://www.aeps-info.com/aeps/ch/index.aspx)。

    [1] 丁明,王波,趙波,等.獨(dú)立風(fēng)光柴儲(chǔ)微網(wǎng)系統(tǒng)容量?jī)?yōu)化配置[J].電網(wǎng)技術(shù),2013,37(3):575-581.

    DING Ming, WANG Bo, ZHAO Bo, et al. Configuration optimization of capacity of standalone PV-wind-diesel-battery hybrid microgrid[J]. Power System Technology, 2013, 37(3): 575-581.

    [2] XU L, RUAN X B, MAO C X, et al. An improved optimal sizing method for wind-solar-battery hybrid power system[J]. IEEE Trans on Sustainable Energy, 2013, 4(3): 774-785.

    [3] 江全元,石慶均,李興鵬,等.風(fēng)光儲(chǔ)獨(dú)立供電系統(tǒng)電源優(yōu)化配置[J].電力自動(dòng)化設(shè)備,2013,33(7):19-26.

    JIANG Quanyuan, SHI Qinjun, LI Xinpeng, et al. Optimal configuration of standalone wind-solar-storage power supply system[J]. Electric Power Automation Equipment, 2013, 33(7): 19-26.

    [4] 馬溪原,吳耀文,方華亮,等.采用改進(jìn)細(xì)菌覓食算法的風(fēng)/光/儲(chǔ)混合微電網(wǎng)電源優(yōu)化配置[J].中國(guó)電機(jī)工程學(xué)報(bào),2011,31(25):17-25.

    MA Xiyuan, WU Yaowen, FANG Hualiang, et al. Optimal sizing of hybrid solar-wind distributed generation in an islanded microgrid using improved bacterial foraging algorithm[J]. Proceedings of the CSEE, 2011, 31(25): 17-25.

    [5] RAMIREZ-DIAZ A, RAMOS-REAL F J, MARRERO G A. Complementarity of electric vehicles and pumped-hydro as energy storage in small isolated energy systems: case of La Palma, Canary Islands[J]. Journal of Modern Power Systems & Clean Energy, 2016, 4(4): 604-614.

    [6] BIE Z H, XIE H P, HU G W, et al. Optimal scheduling of power systems considering demand response[J]. Journal of Modern Power Systems and Clean Energy, 2016, 4(2): 180-187.

    [7] CHEN H C. Optimum capacity determination of stand-alone hybrid generation system considering cost and reliability[J]. Applied Energy, 2013, 103: 155-164.

    [8] 余舟子,郭力,劉文建,等.獨(dú)立微網(wǎng)隨機(jī)優(yōu)化規(guī)劃軟件及其實(shí)現(xiàn)[J].電力系統(tǒng)自動(dòng)化,2015,39(6):18-23.DOI:10.7500/AEPS20140207002.

    YU Zhouzi, GUO Li, LIU Wenjian, et al. A stochastic optimal planning software for stand-alone microgrid and its implementation[J]. Automation of Electric Power Systems, 2015, 39(6): 18-23. DOI: 10.7500/AEPS20140207002.

    [9] 李豐,張粒子,舒雋.考慮間歇式電源與儲(chǔ)能的隨機(jī)柔性優(yōu)化調(diào)度方法[J].電力系統(tǒng)自動(dòng)化,2014,38(5):1-7.DOI:10.7500/AEPS20130717011.

    LI Feng, ZHANG Lizi, SHU Jun. A flexible stochastic optimal scheduling approach considering intermittent power and energy storage[J]. Automation of Electric Power Systems, 2014, 38(5): 1-7. DOI: 10.7500/AEPS20130717011.

    [10] 李斯,周任軍,童小嬌,等.基于盒式集合魯棒優(yōu)化的風(fēng)電并網(wǎng)最大裝機(jī)容量[J].電網(wǎng)技術(shù),2011,35(12):208-213.

    LI Si, ZHOU Renjun, TONG Xiaojiao, et al. Robust optimization with box set for maximum installed capacity of wind farm connected to grid[J]. Power System Technology, 2011, 35(12): 208-213.

    [11] 謝鵬,彭春華,于蓉.大規(guī)模間歇式電源接入電網(wǎng)多目標(biāo)魯棒優(yōu)化調(diào)度[J].電網(wǎng)技術(shù),2014,38(6):1479-1484.

    XIE Peng, PENG Chunhua, YU Rong. Multi-objective robust optimized scheduling of power grid connected with large-scale intermittent power sources[J]. Power System Technology, 2014, 38(6): 1479-1484.

    [12] 向月,劉俊勇,魏震波,等.考慮可再生能源出力不確定性的微電網(wǎng)能量?jī)?yōu)化魯棒模型[J].中國(guó)電機(jī)工程學(xué)報(bào),2014,34(19):3063-3072.

    XIANG Yue, LIU Junyong, WEI Zhenbo, et al. Robust model of microgrid energy optimization with uncertain renewable energy sources[J]. Proceedings of the CSEE, 2014, 34(19): 3063-3072.

    [13] 茆美琴,孫樹娟,蘇建徽.包含電動(dòng)汽車的風(fēng)/光/儲(chǔ)微電網(wǎng)經(jīng)濟(jì)性分析[J].電力系統(tǒng)自動(dòng)化,2011,35(14):30-35.

    MAO meiqing, SUN shujuan, SU jianhui. Economic analysis of a microgrid with wind/photovoltaic/storages and electric vehicles[J]. Automation of Electric Power Systems, 2011, 35(14): 30-35.

    [14] 趙波.微電網(wǎng)優(yōu)化配置關(guān)鍵技術(shù)及應(yīng)用[M].北京:科學(xué)出版社,2015.

    [15] 王成山,洪博文,郭力.不同場(chǎng)景下的光蓄微電網(wǎng)調(diào)度策略[J].電網(wǎng)技術(shù),2013,37(7):1775-1782.

    WANG Chengshan, HONG Bowen, GUO Li. Dispatch strategies of PV-battery microgrid in different scenarios[J]. Power System Technology, 2013, 37(7): 1775-1782.

    [16] 郭力,富曉鵬,李霞林,等.獨(dú)立交流微網(wǎng)中電池儲(chǔ)能與柴油發(fā)電機(jī)的協(xié)調(diào)控制[J].中國(guó)電機(jī)工程學(xué)報(bào),2012,32(25):70-78.

    GUO Li, FU Xiaopeng, LI Xialing, et al. Coordinated control of battery storage and diesel generators in isolated AC microgrid systems[J]. Proceedings of the CSEE, 2012, 32(25): 70-78.

    [17] 梅生偉,郭文濤,王瑩瑩,等.一類電力系統(tǒng)魯棒優(yōu)化問(wèn)題的博弈模型及應(yīng)用實(shí)例[J].中國(guó)電機(jī)工程學(xué)報(bào),2013,33(19):47-56.

    MEI Shengwei, GUO Wentao, WANG Yingying, et al. A game model for robust optimization of power systems and its application[J]. Proceedings of the CSEE, 2013, 33(19): 47-56.

    [18] BRNDERS J F. Partitioning procedures for solving mixed-variables programing problems[J]. Numerische Mathmatik, 1962, 4(1): 238-252.

    [19] GUPT A, AHMED S, CHEON M S, et al. Solving mixed integer bilinear problems using MILP formulations[J]. Siam Journal on Optimization, 2013, 23(2): 721-744.

    CapacityConfigurationModelforRobustOptimizationofStand-aloneMicrogridBasedonBendersDecomposition

    LIUSiyi1,ZHAOBo2,WANGXiangjin2,QIUHaifeng2,LINXiangning1

    (1. State Key Laboratory of Advanced Electromagnetic Engineering and Technology (Huazhong University of Science and Technology), Wuhan 430074, China; 2. Electric Power Research Institute of State Grid Zhejiang Electric Power Company, Hangzhou 310014, China)

    To improve the robustness of power source configuration for stand-alone microgrid, a double layer configuration model for robust power source capacity is built considering the intermittency and volatility of wind, solar radiation and load. The outer layer is the master problem regarding power source capacity optimization, the inner layer is the subproblem regarding operational optimization including uncertainty sets of wind, solar radiation and load. The dual problem is established to decouple the subproblem and the bilinear terms in the subproblem are converted to linear terms using the big M method, iterate the master problem and subproblem alternatively until the optimal solution is obtained. The proposed method is compared with improved non-dominant sorting genetic algorithm (NSGA-Ⅱ) in the test case. The testing result shows the robust performance and solving efficiency of the proposed model is higher than usual. The rationality and effectiveness of the model is verified.

    This work is supported by National High Technology Research and Development Program of China (863 Program) (No. 2015AA050104).

    microgrid; optimal configuration; robust optimization; Benders decomposition; uncertainty sets

    2016-12-28;

    2017-04-11。

    上網(wǎng)日期: 2017-06-27。

    國(guó)家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)資助項(xiàng)目(2015AA050104)。

    劉思夷(1993—),女,碩士研究生,主要研究方向:分布式發(fā)電及微電網(wǎng)優(yōu)化配置。E-mail: liusiyi@163.com

    趙 波(1977—),男,通信作者,博士,教授級(jí)高級(jí)工程師,主要研究方向:分布式能源及微電網(wǎng)關(guān)鍵技術(shù)。E-mail: zhaobozju@163.com

    汪湘晉(1992—),男,碩士研究生,主要研究方向:微電網(wǎng)優(yōu)化配置及儲(chǔ)能模型。E-mail: hfwangxj@126.com

    (編輯萬(wàn)志超)

    猜你喜歡
    臺(tái)數(shù)魯棒儲(chǔ)能
    相變儲(chǔ)能材料的應(yīng)用
    煤氣與熱力(2021年6期)2021-07-28 07:21:24
    基于學(xué)習(xí)的魯棒自適應(yīng)評(píng)判控制研究進(jìn)展
    最優(yōu)開機(jī)臺(tái)數(shù)對(duì)水利泵站的節(jié)能研究
    儲(chǔ)能技術(shù)在電力系統(tǒng)中的應(yīng)用
    2016年XX市工業(yè)鍋爐內(nèi)部檢驗(yàn)缺陷統(tǒng)計(jì)及原因分析
    儲(chǔ)能真要起飛了?
    能源(2017年12期)2018-01-31 01:42:59
    目標(biāo)魯棒識(shí)別的抗旋轉(zhuǎn)HDO 局部特征描述
    直流儲(chǔ)能型準(zhǔn)Z源光伏并網(wǎng)逆變器
    基于Cauchy魯棒函數(shù)的UKF改進(jìn)算法
    目標(biāo)軌跡更新的點(diǎn)到點(diǎn)魯棒迭代學(xué)習(xí)控制
    高清日韩中文字幕在线| 黑人高潮一二区| 亚洲精品国产成人久久av| 国内精品宾馆在线| 涩涩av久久男人的天堂| 伦精品一区二区三区| 亚洲电影在线观看av| 亚洲欧美精品专区久久| 大片电影免费在线观看免费| 精品一区二区三区视频在线| 成人综合一区亚洲| 亚洲欧美精品自产自拍| 99久久中文字幕三级久久日本| 身体一侧抽搐| av又黄又爽大尺度在线免费看| 成人漫画全彩无遮挡| 激情五月婷婷亚洲| 夫妻性生交免费视频一级片| 久久国产精品大桥未久av | 日本欧美视频一区| 亚洲精品一区蜜桃| 中文字幕久久专区| 亚洲久久久国产精品| 国产白丝娇喘喷水9色精品| 永久网站在线| 午夜福利在线在线| 日本vs欧美在线观看视频 | 精品久久久久久电影网| 国产精品一及| 亚洲av.av天堂| av在线蜜桃| 精品国产一区二区三区久久久樱花 | 超碰av人人做人人爽久久| 深爱激情五月婷婷| 欧美日韩视频高清一区二区三区二| 女人十人毛片免费观看3o分钟| 久久午夜福利片| 亚洲婷婷狠狠爱综合网| 少妇猛男粗大的猛烈进出视频| 成人无遮挡网站| 久久国内精品自在自线图片| 少妇 在线观看| 欧美 日韩 精品 国产| 伊人久久国产一区二区| 夫妻午夜视频| 亚洲综合精品二区| 免费黄频网站在线观看国产| 国产免费福利视频在线观看| 2018国产大陆天天弄谢| 少妇 在线观看| 午夜免费男女啪啪视频观看| 久久热精品热| 丰满迷人的少妇在线观看| 韩国av在线不卡| 免费看不卡的av| 日韩av不卡免费在线播放| 国产亚洲午夜精品一区二区久久| 韩国av在线不卡| 蜜桃在线观看..| 亚洲一区二区三区欧美精品| 在线观看美女被高潮喷水网站| 国产中年淑女户外野战色| 午夜免费观看性视频| 人妻夜夜爽99麻豆av| 国产亚洲91精品色在线| 51国产日韩欧美| 赤兔流量卡办理| 中文天堂在线官网| 成人18禁高潮啪啪吃奶动态图 | 又粗又硬又长又爽又黄的视频| 亚洲四区av| 国产一区有黄有色的免费视频| 黑丝袜美女国产一区| 国产精品一及| 国产精品99久久99久久久不卡 | 亚洲av成人精品一区久久| 全区人妻精品视频| 大码成人一级视频| 建设人人有责人人尽责人人享有的 | 午夜日本视频在线| 国产精品免费大片| 国产淫片久久久久久久久| 性色av一级| 特大巨黑吊av在线直播| 国产精品女同一区二区软件| 国产av国产精品国产| 少妇裸体淫交视频免费看高清| 99久久精品热视频| 色哟哟·www| 成人无遮挡网站| 亚洲精品中文字幕在线视频 | 97热精品久久久久久| 少妇熟女欧美另类| 一个人免费看片子| 国产一区有黄有色的免费视频| 夜夜爽夜夜爽视频| 国产午夜精品一二区理论片| 伦理电影大哥的女人| 久久人妻熟女aⅴ| 久久99蜜桃精品久久| 国产精品国产三级国产av玫瑰| 国产精品精品国产色婷婷| 亚洲综合精品二区| 91精品国产九色| 免费少妇av软件| 日韩国内少妇激情av| av网站免费在线观看视频| 男女下面进入的视频免费午夜| 多毛熟女@视频| av国产久精品久网站免费入址| 国产视频内射| 欧美97在线视频| 免费少妇av软件| 激情 狠狠 欧美| 日韩亚洲欧美综合| 丝袜喷水一区| 久久av网站| 日本与韩国留学比较| 激情五月婷婷亚洲| 国产中年淑女户外野战色| 国产亚洲av片在线观看秒播厂| 国产av一区二区精品久久 | 网址你懂的国产日韩在线| 国产爱豆传媒在线观看| 久久久久精品久久久久真实原创| 免费人成在线观看视频色| 成人美女网站在线观看视频| 成年女人在线观看亚洲视频| 五月开心婷婷网| 边亲边吃奶的免费视频| 黄片wwwwww| 97在线视频观看| 久久人人爽av亚洲精品天堂 | 成人美女网站在线观看视频| 日本与韩国留学比较| 国产成人免费无遮挡视频| 国产 一区 欧美 日韩| 午夜福利在线观看免费完整高清在| 久久人人爽人人爽人人片va| 久热久热在线精品观看| 人体艺术视频欧美日本| 狂野欧美激情性xxxx在线观看| 亚洲成人一二三区av| 女性生殖器流出的白浆| 国产亚洲91精品色在线| 特大巨黑吊av在线直播| 亚洲国产高清在线一区二区三| 自拍欧美九色日韩亚洲蝌蚪91 | 中国美白少妇内射xxxbb| 亚洲天堂av无毛| 国产精品人妻久久久影院| 2021少妇久久久久久久久久久| 美女中出高潮动态图| 亚洲真实伦在线观看| 国精品久久久久久国模美| 高清毛片免费看| 伦精品一区二区三区| 成人亚洲欧美一区二区av| 建设人人有责人人尽责人人享有的 | 久久精品久久精品一区二区三区| av国产久精品久网站免费入址| 日韩 亚洲 欧美在线| 一个人免费看片子| 亚洲无线观看免费| 亚洲av日韩在线播放| 一区二区av电影网| 日韩强制内射视频| 乱系列少妇在线播放| 中文在线观看免费www的网站| 色综合色国产| 中文字幕av成人在线电影| 十分钟在线观看高清视频www | 欧美日韩国产mv在线观看视频 | 日本av免费视频播放| 夫妻午夜视频| 国产精品熟女久久久久浪| 久久精品国产鲁丝片午夜精品| 在线精品无人区一区二区三 | 少妇熟女欧美另类| 久久久久久久久久久免费av| 精品一区在线观看国产| 老女人水多毛片| 亚洲精品视频女| 久久精品久久精品一区二区三区| 亚洲第一av免费看| 国产午夜精品久久久久久一区二区三区| 毛片一级片免费看久久久久| 小蜜桃在线观看免费完整版高清| 少妇人妻久久综合中文| 韩国高清视频一区二区三区| 午夜福利在线观看免费完整高清在| 18禁在线播放成人免费| 丝袜喷水一区| 99热国产这里只有精品6| 久久久色成人| 国产精品一区二区三区四区免费观看| 欧美另类一区| 最近最新中文字幕大全电影3| 亚洲av国产av综合av卡| 美女内射精品一级片tv| 久久这里有精品视频免费| 日本午夜av视频| 亚洲欧美清纯卡通| 色5月婷婷丁香| 国产成人91sexporn| 水蜜桃什么品种好| 肉色欧美久久久久久久蜜桃| 91精品国产国语对白视频| 青春草视频在线免费观看| 日本午夜av视频| 国产淫片久久久久久久久| 岛国毛片在线播放| av一本久久久久| 精品亚洲乱码少妇综合久久| 亚洲一区二区三区欧美精品| 国产精品熟女久久久久浪| 久久精品国产亚洲av涩爱| 爱豆传媒免费全集在线观看| 亚洲图色成人| 国产精品无大码| 在现免费观看毛片| 黄色欧美视频在线观看| 99视频精品全部免费 在线| 自拍偷自拍亚洲精品老妇| 免费不卡的大黄色大毛片视频在线观看| 视频中文字幕在线观看| 国产免费一区二区三区四区乱码| 亚洲av在线观看美女高潮| 尾随美女入室| 日本午夜av视频| 黄色怎么调成土黄色| 嘟嘟电影网在线观看| 欧美高清成人免费视频www| 激情 狠狠 欧美| 久久影院123| 黄色配什么色好看| 亚洲精品久久久久久婷婷小说| 在线观看三级黄色| 一级毛片我不卡| 久久久久久久久久久丰满| 免费大片黄手机在线观看| 日本色播在线视频| 青青草视频在线视频观看| 亚洲欧美一区二区三区国产| 另类亚洲欧美激情| av网站免费在线观看视频| 男女国产视频网站| 少妇熟女欧美另类| 熟女av电影| 国产一区二区三区av在线| 国产av国产精品国产| 亚洲人成网站在线播| 人体艺术视频欧美日本| 卡戴珊不雅视频在线播放| 国产在线男女| 五月玫瑰六月丁香| 免费黄频网站在线观看国产| 久久久成人免费电影| 亚洲成人一二三区av| 成人毛片a级毛片在线播放| 国产精品久久久久久av不卡| 久久av网站| 九色成人免费人妻av| 精品99又大又爽又粗少妇毛片| 久久久亚洲精品成人影院| 99久久综合免费| 欧美高清成人免费视频www| 97超视频在线观看视频| 亚洲精品一区蜜桃| 国产精品av视频在线免费观看| 国产亚洲一区二区精品| 人妻夜夜爽99麻豆av| 国产精品久久久久久av不卡| 国产欧美亚洲国产| 激情 狠狠 欧美| 大又大粗又爽又黄少妇毛片口| 99久久精品热视频| 成年av动漫网址| 亚洲精品久久久久久婷婷小说| 在现免费观看毛片| av又黄又爽大尺度在线免费看| 国产精品伦人一区二区| 国产精品成人在线| 久久久久精品性色| 两个人的视频大全免费| 99久久精品国产国产毛片| 亚洲色图av天堂| 99九九线精品视频在线观看视频| 亚洲精品一二三| 精品午夜福利在线看| 久久精品国产a三级三级三级| 亚洲精品日韩av片在线观看| 亚洲av综合色区一区| 91午夜精品亚洲一区二区三区| 一级毛片我不卡| 五月开心婷婷网| 天堂8中文在线网| 精品一品国产午夜福利视频| 久久久国产一区二区| videos熟女内射| 最黄视频免费看| 色视频www国产| 成人毛片a级毛片在线播放| 春色校园在线视频观看| 一级毛片黄色毛片免费观看视频| 一本久久精品| 一个人看的www免费观看视频| 少妇裸体淫交视频免费看高清| 国产午夜精品一二区理论片| 国产亚洲最大av| 男女无遮挡免费网站观看| 中国三级夫妇交换| 中文在线观看免费www的网站| 少妇人妻 视频| 一级av片app| 只有这里有精品99| 久久精品久久久久久久性| 久久热精品热| 狂野欧美激情性bbbbbb| 国产永久视频网站| 亚洲av免费高清在线观看| 欧美日韩综合久久久久久| 美女福利国产在线 | 国产一区二区在线观看日韩| 中文字幕制服av| 国产永久视频网站| 日韩av免费高清视频| 国产精品.久久久| 亚洲精品自拍成人| 少妇猛男粗大的猛烈进出视频| 国产伦精品一区二区三区四那| 日本色播在线视频| 黑人猛操日本美女一级片| 国产69精品久久久久777片| 国产欧美亚洲国产| 亚州av有码| 免费在线观看成人毛片| 女人久久www免费人成看片| 丰满人妻一区二区三区视频av| 亚洲国产色片| 色5月婷婷丁香| 一级毛片 在线播放| 99热6这里只有精品| 亚洲国产日韩一区二区| 各种免费的搞黄视频| 亚洲精品第二区| 国产成人精品福利久久| 各种免费的搞黄视频| 亚洲精品成人av观看孕妇| 国产精品一区二区在线观看99| 亚洲内射少妇av| 内地一区二区视频在线| 亚洲欧美精品自产自拍| 日本一二三区视频观看| 国产精品一区二区性色av| 日本欧美国产在线视频| 三级国产精品欧美在线观看| 久久久久精品性色| 男人舔奶头视频| 精品人妻熟女av久视频| 国产极品天堂在线| 欧美3d第一页| 看非洲黑人一级黄片| 国产精品伦人一区二区| 亚洲欧美精品自产自拍| 黄片无遮挡物在线观看| 亚洲伊人久久精品综合| 在线免费十八禁| 一个人看的www免费观看视频| 久久影院123| 97在线人人人人妻| 婷婷色麻豆天堂久久| 国产成人a∨麻豆精品| 久久精品国产鲁丝片午夜精品| 我的老师免费观看完整版| 欧美zozozo另类| 成年女人在线观看亚洲视频| 又大又黄又爽视频免费| 国产精品99久久99久久久不卡 | 一个人看的www免费观看视频| 日本vs欧美在线观看视频 | 中文字幕精品免费在线观看视频 | 国产亚洲5aaaaa淫片| 好男人视频免费观看在线| 欧美高清成人免费视频www| 亚洲综合色惰| 久久人人爽人人片av| 久久久久久久久久人人人人人人| 亚洲国产欧美人成| 狠狠精品人妻久久久久久综合| 欧美最新免费一区二区三区| 国产精品女同一区二区软件| 国产成人a区在线观看| 五月伊人婷婷丁香| 亚洲欧美中文字幕日韩二区| 亚洲经典国产精华液单| 色吧在线观看| 亚洲av在线观看美女高潮| 国产精品av视频在线免费观看| 尤物成人国产欧美一区二区三区| av又黄又爽大尺度在线免费看| 男人添女人高潮全过程视频| 国产精品爽爽va在线观看网站| 国产伦理片在线播放av一区| 国产av国产精品国产| 成人亚洲精品一区在线观看 | 免费看光身美女| 久久国产亚洲av麻豆专区| 人人妻人人澡人人爽人人夜夜| 午夜老司机福利剧场| 国产熟女欧美一区二区| 在线 av 中文字幕| 汤姆久久久久久久影院中文字幕| 极品教师在线视频| 少妇人妻久久综合中文| 亚洲三级黄色毛片| 亚洲精品国产色婷婷电影| 国产成人一区二区在线| av卡一久久| 熟女电影av网| 国产精品99久久99久久久不卡 | 交换朋友夫妻互换小说| 韩国高清视频一区二区三区| 男女下面进入的视频免费午夜| 国产乱人偷精品视频| 18禁裸乳无遮挡动漫免费视频| 亚洲三级黄色毛片| 成人漫画全彩无遮挡| 亚洲国产日韩一区二区| 午夜福利影视在线免费观看| 精品人妻偷拍中文字幕| 一级av片app| 草草在线视频免费看| 亚洲国产精品国产精品| 国产男女内射视频| 中国美白少妇内射xxxbb| 啦啦啦中文免费视频观看日本| 三级经典国产精品| 日韩免费高清中文字幕av| 两个人的视频大全免费| 麻豆乱淫一区二区| 日韩中字成人| 精品亚洲成a人片在线观看 | 亚洲国产av新网站| 亚洲不卡免费看| 成人国产av品久久久| 久久国内精品自在自线图片| av网站免费在线观看视频| 日韩一区二区三区影片| 视频区图区小说| 免费黄频网站在线观看国产| 亚洲精品国产成人久久av| 伊人久久精品亚洲午夜| 精品久久久精品久久久| 人妻夜夜爽99麻豆av| 三级国产精品片| av国产免费在线观看| 中文字幕免费在线视频6| 国产成人免费观看mmmm| 国产欧美日韩精品一区二区| 久久久久人妻精品一区果冻| 国产精品.久久久| 免费观看的影片在线观看| 在线观看av片永久免费下载| 男人添女人高潮全过程视频| 欧美国产精品一级二级三级 | 国产 一区精品| 亚洲av在线观看美女高潮| 自拍欧美九色日韩亚洲蝌蚪91 | 黄片无遮挡物在线观看| 又爽又黄a免费视频| 国产精品免费大片| 亚洲最大成人中文| 亚洲精品国产av成人精品| 久久热精品热| 纯流量卡能插随身wifi吗| 爱豆传媒免费全集在线观看| 久久亚洲国产成人精品v| 国产高清有码在线观看视频| 亚洲第一区二区三区不卡| 狂野欧美白嫩少妇大欣赏| 五月伊人婷婷丁香| 欧美一级a爱片免费观看看| 亚洲精品国产色婷婷电影| 黄色一级大片看看| 国内精品宾馆在线| 91aial.com中文字幕在线观看| 伊人久久国产一区二区| 激情五月婷婷亚洲| 一个人看的www免费观看视频| 中文乱码字字幕精品一区二区三区| 日产精品乱码卡一卡2卡三| 中文字幕免费在线视频6| av又黄又爽大尺度在线免费看| 日韩在线高清观看一区二区三区| 美女内射精品一级片tv| 日韩一本色道免费dvd| 欧美日韩亚洲高清精品| 欧美精品国产亚洲| 99热这里只有是精品50| 日韩三级伦理在线观看| 18禁动态无遮挡网站| 国产久久久一区二区三区| 免费高清在线观看视频在线观看| 久久久久国产精品人妻一区二区| 男女免费视频国产| 久久国内精品自在自线图片| 免费看光身美女| 久久人人爽人人爽人人片va| 免费人成在线观看视频色| 久久久久久久久大av| 女的被弄到高潮叫床怎么办| 国产精品成人在线| 中文字幕免费在线视频6| 日本vs欧美在线观看视频 | 在线观看人妻少妇| 久久久成人免费电影| 干丝袜人妻中文字幕| 欧美日韩精品成人综合77777| 中文字幕亚洲精品专区| 深夜a级毛片| 国产欧美日韩一区二区三区在线 | 国产在线免费精品| av国产久精品久网站免费入址| 高清午夜精品一区二区三区| 日日摸夜夜添夜夜爱| 全区人妻精品视频| 亚洲四区av| 国产爱豆传媒在线观看| 亚洲精品第二区| 日日啪夜夜撸| 国产69精品久久久久777片| 婷婷色综合大香蕉| 丝瓜视频免费看黄片| 国产精品国产三级国产av玫瑰| 在线亚洲精品国产二区图片欧美 | 国产成人精品一,二区| av又黄又爽大尺度在线免费看| 免费少妇av软件| 色哟哟·www| 免费看不卡的av| 欧美性感艳星| 不卡视频在线观看欧美| 国精品久久久久久国模美| .国产精品久久| 99久国产av精品国产电影| 日韩欧美 国产精品| 国产精品一区www在线观看| 在线观看av片永久免费下载| 亚洲精品第二区| 亚洲欧美日韩无卡精品| 91精品国产九色| 久久99热这里只有精品18| 久久久a久久爽久久v久久| 又粗又硬又长又爽又黄的视频| 秋霞伦理黄片| 国产高清国产精品国产三级 | 在线观看美女被高潮喷水网站| 亚洲最大成人中文| 国产免费一区二区三区四区乱码| 欧美三级亚洲精品| 99久久精品国产国产毛片| 亚洲精品亚洲一区二区| videos熟女内射| 日韩成人av中文字幕在线观看| 欧美xxⅹ黑人| 国产91av在线免费观看| 伦理电影免费视频| 熟女电影av网| 久热久热在线精品观看| 亚洲国产欧美人成| 欧美区成人在线视频| 热re99久久精品国产66热6| 日韩三级伦理在线观看| 中文资源天堂在线| 国产免费视频播放在线视频| 99精国产麻豆久久婷婷| 日韩人妻高清精品专区| 热99国产精品久久久久久7| 亚洲电影在线观看av| 乱系列少妇在线播放| 汤姆久久久久久久影院中文字幕| 国产无遮挡羞羞视频在线观看| 国内揄拍国产精品人妻在线| 国产大屁股一区二区在线视频| 国产国拍精品亚洲av在线观看| 欧美变态另类bdsm刘玥| 超碰av人人做人人爽久久| 三级经典国产精品| 国产精品嫩草影院av在线观看| 日韩一本色道免费dvd| 国产精品久久久久久久久免| 高清日韩中文字幕在线| 美女中出高潮动态图| 欧美一区二区亚洲| 女性生殖器流出的白浆| 99久国产av精品国产电影| 汤姆久久久久久久影院中文字幕| 欧美日韩亚洲高清精品| 丰满人妻一区二区三区视频av| 国产av一区二区精品久久 | 国产精品一区www在线观看| 自拍偷自拍亚洲精品老妇| 国产91av在线免费观看| 舔av片在线| 三级国产精品欧美在线观看| 丝瓜视频免费看黄片| 在线观看美女被高潮喷水网站| 午夜福利高清视频| 久久久a久久爽久久v久久| 亚洲aⅴ乱码一区二区在线播放| 2022亚洲国产成人精品| 日韩视频在线欧美| 亚洲国产最新在线播放| 五月伊人婷婷丁香|