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

    機(jī)組組合問(wèn)題的仿射可調(diào)整魯棒優(yōu)化模型與算法

    2017-07-06 13:01:41李利利涂孟夫謝麗榮
    電力工程技術(shù) 2017年3期
    關(guān)鍵詞:魯棒出力約束

    李利利, 丁 恰, 涂孟夫, 謝麗榮

    (南瑞集團(tuán)公司(國(guó)網(wǎng)電力科學(xué)研究院),江蘇 南京 211106)

    ?

    ·電網(wǎng)技術(shù)·

    機(jī)組組合問(wèn)題的仿射可調(diào)整魯棒優(yōu)化模型與算法

    李利利, 丁 恰, 涂孟夫, 謝麗榮

    (南瑞集團(tuán)公司(國(guó)網(wǎng)電力科學(xué)研究院),江蘇 南京 211106)

    為了有效應(yīng)對(duì)電力系統(tǒng)調(diào)度決策中的不確定因素,尤其是大規(guī)模間歇式能源并網(wǎng)所帶來(lái)的不確定性,提出基于仿射可調(diào)整魯棒優(yōu)化理論的不確定機(jī)組組合求解方法。建立了不確定機(jī)組組合問(wèn)題的仿射可調(diào)整魯棒優(yōu)化模型,利用線性決策規(guī)則建立決策變量與不確定參數(shù)之間的仿射關(guān)系,從而將兩階段問(wèn)題轉(zhuǎn)化為單個(gè)階段優(yōu)化問(wèn)題,在此基礎(chǔ)上,采用對(duì)偶理論將模型轉(zhuǎn)化為可以直接求解的標(biāo)準(zhǔn)混合整數(shù)規(guī)劃模型。通過(guò)標(biāo)準(zhǔn)算例測(cè)試,驗(yàn)證了該方法的有效性。

    發(fā)電計(jì)劃;不確定性;魯棒優(yōu)化;安全約束機(jī)組組合;線性決策規(guī)則

    0 引言

    為應(yīng)對(duì)能源危機(jī)和環(huán)境問(wèn)題,以風(fēng)電為代表的新能源發(fā)電近年來(lái)得到了快速發(fā)展。與常規(guī)電源相比,新能源發(fā)電具有間歇性、波動(dòng)性、不完全可預(yù)測(cè)性等調(diào)度運(yùn)行特點(diǎn)。這些特點(diǎn)顯著增加了電力系統(tǒng)運(yùn)行的不確定性,給安全約束機(jī)組組合(SCUC)帶來(lái)了新的挑戰(zhàn)和要求[1-4]。為此,當(dāng)前迫切需要建立考慮新能源發(fā)電不確定性的調(diào)度計(jì)劃優(yōu)化技術(shù)手段。

    從運(yùn)籌學(xué)的角度,對(duì)不確定問(wèn)題的處理主要有隨機(jī)規(guī)劃、模糊規(guī)劃、魯棒優(yōu)化等方法[5]。其中,魯棒優(yōu)化以其對(duì)問(wèn)題不確定性數(shù)據(jù)的分布信息不作過(guò)多假設(shè)、最優(yōu)解對(duì)不確定集合內(nèi)的任意元素都保證約束可行性的特點(diǎn),具有突出的優(yōu)勢(shì)[6]。目前,魯棒優(yōu)化理論已經(jīng)被用于解決電力系統(tǒng)運(yùn)行中的多種不確定性問(wèn)題,包括:考慮多個(gè)發(fā)電機(jī)組隨機(jī)故障的機(jī)組組合建模[7]、考慮電價(jià)不確定性的電廠競(jìng)價(jià)曲線構(gòu)建[8]、考慮相關(guān)規(guī)劃不確定性的電網(wǎng)支持插入式混合動(dòng)力電動(dòng)汽車(PHEVs)潛力評(píng)估[9]、風(fēng)電并網(wǎng)最大裝機(jī)容量計(jì)算[10]等。

    與此同時(shí),應(yīng)用魯棒優(yōu)化理論解決不確定機(jī)組組合問(wèn)題也逐漸得到國(guó)內(nèi)外學(xué)者的關(guān)注。文獻(xiàn)[11]提出了電力系統(tǒng)魯棒經(jīng)濟(jì)調(diào)度的理論框架,為大規(guī)模新能源并網(wǎng)后電網(wǎng)調(diào)度提供了一種解決方案。發(fā)電計(jì)劃在時(shí)序上可以劃分為機(jī)組的啟停優(yōu)化與實(shí)時(shí)出力分配兩部分內(nèi)容,機(jī)組啟停決策應(yīng)該在不確定因素實(shí)現(xiàn)之前獲得具體的結(jié)果,而出力決策則需要依賴于不確定因素的實(shí)際值,根據(jù)功率平衡實(shí)時(shí)計(jì)算獲得。文獻(xiàn)[12,13]提出了機(jī)組組合的兩階段魯棒優(yōu)化模型,能夠在負(fù)荷或新能源發(fā)電不確定條件下提供魯棒的開(kāi)停機(jī)決策結(jié)果。文獻(xiàn)[14]將兩階段魯棒優(yōu)化方法用于求解故障約束機(jī)組組合問(wèn)題。文獻(xiàn)[15]同時(shí)考慮風(fēng)力發(fā)電與價(jià)格彈性負(fù)荷的不確定性,提出了多階段魯棒機(jī)組組合模型。

    在算法方面,兩階段魯棒優(yōu)化問(wèn)題的模型為混合整數(shù)非線性規(guī)劃且具有較高的復(fù)雜度,難以采用成熟的算法軟件直接求解。為了便于計(jì)算,上述文獻(xiàn)均采用分解算法,如Benders對(duì)偶分解技術(shù)[12-14]或原始分解技術(shù)[15],將問(wèn)題分解為主問(wèn)題與子問(wèn)題后迭代求解。然而,分解迭代算法需要在迭代過(guò)程中巧妙構(gòu)建割平面約束以加速收斂[13],魯棒優(yōu)化問(wèn)題的復(fù)雜性會(huì)使得計(jì)算性能得不到保障。同時(shí),分解后的子問(wèn)題是一個(gè)包含雙線性目標(biāo)函數(shù)的非線性規(guī)劃問(wèn)題,難以求取全局最優(yōu)解。

    本文探索了一條直接求解魯棒機(jī)組組合問(wèn)題的可行途徑。建立了不確定機(jī)組組合問(wèn)題的仿射可調(diào)整魯棒優(yōu)化模型,利用線性決策規(guī)則建立決策變量與不確定參數(shù)之間的仿射關(guān)系,從而將兩階段問(wèn)題轉(zhuǎn)化為單個(gè)階段優(yōu)化問(wèn)題,在此基礎(chǔ)上,采用對(duì)偶理論將模型轉(zhuǎn)化為可以直接求解的標(biāo)準(zhǔn)混合整數(shù)規(guī)劃模型。通過(guò)標(biāo)準(zhǔn)算例測(cè)試,驗(yàn)證了該方法的有效性。

    1 可調(diào)整魯棒機(jī)組組合模型

    可調(diào)整魯棒優(yōu)化的主要思想是把決策變量劃分為可調(diào)整變量和不可調(diào)整變量進(jìn)行分階段決策。由文獻(xiàn)[17],依賴于不確定集并隨著不確定參數(shù)變化而調(diào)整其值的變量稱為可調(diào)整變量,而在不確定參數(shù)實(shí)現(xiàn)前就已確定的變量稱為不可調(diào)整變量。

    對(duì)于不確定機(jī)組組合問(wèn)題而言,機(jī)組狀態(tài)應(yīng)該在不確定因素實(shí)現(xiàn)之前確定具體的結(jié)果,屬于不可調(diào)整變量。機(jī)組出力需要根據(jù)負(fù)荷平衡的等式約束,在不確定因素實(shí)現(xiàn)之后才能獲得其具體的結(jié)果,屬于可調(diào)整變量。

    根據(jù)可調(diào)整魯棒優(yōu)化理論,建立不確定機(jī)組組合的魯棒優(yōu)化模型。為簡(jiǎn)化敘述,模型采用與文獻(xiàn)[13]類似的緊湊表達(dá)形式,具體如下:

    (1)

    約束條件為:

    Fx≤f

    (2)

    Hy(d)≤h(d) ?d∈D

    (3)

    Ax+By(d)≤g?d∈D

    (4)

    Ty(d)+Wd=l?d∈D

    (5)

    式中:x表示與機(jī)組狀態(tài)相關(guān)的二進(jìn)制變量,包括機(jī)組在各時(shí)段的啟停狀態(tài)、開(kāi)機(jī)標(biāo)志與停機(jī)標(biāo)志;y表示與機(jī)組出力相關(guān)的連續(xù)變量;d表示電網(wǎng)運(yùn)行中的不確定參數(shù);D為不確定集合;c,b,F,f,A,B,g,H,h,T,W,l均為緊湊模型表達(dá)中的各類參數(shù)向量。式(1)描述的目標(biāo)函數(shù)為最小化機(jī)組啟停成本與不確定集合中最壞情況下的機(jī)組出力成本;式(2)為與機(jī)組狀態(tài)相關(guān)的約束條件,包括啟停狀態(tài)約束、最小開(kāi)停機(jī)時(shí)間約束;式(3)為與機(jī)組出力相關(guān)的約束條件,包括備用需求約束、機(jī)組爬坡約束、直流潮流電網(wǎng)安全約束;式(4)為與機(jī)組狀態(tài)和出力相耦合的約束條件,如機(jī)組出力上下限約束;式(5)為系統(tǒng)平衡約束。

    根據(jù)上述模型,在機(jī)組組合優(yōu)化時(shí)考慮了不確定集合D所描述的所有可能的情況,機(jī)組啟停決策x均保持可行性,因此結(jié)果是魯棒的。在經(jīng)濟(jì)調(diào)度優(yōu)化中,有功出力y(d)是不確定參數(shù)d的函數(shù),能夠隨著不確定性參數(shù)d的實(shí)際值而變化,因而對(duì)于不確定集合中的任意實(shí)現(xiàn)而言是可調(diào)整的,這與電網(wǎng)運(yùn)行中的實(shí)時(shí)調(diào)度功能是一致的。

    需要說(shuō)明的是,應(yīng)用魯棒優(yōu)化理論解決電力系統(tǒng)實(shí)際問(wèn)題時(shí),其決策的保守性也引起了運(yùn)行人員的擔(dān)憂,如何構(gòu)造合理的不確定集合尤為重要。以新能源發(fā)電為例,現(xiàn)實(shí)中不太可能發(fā)生所有預(yù)測(cè)同時(shí)達(dá)到邊界的情況,不確定集合的構(gòu)建可以綜合考慮新能源發(fā)電的時(shí)空相關(guān)性,加入不確定預(yù)算約束,限制不確定性對(duì)預(yù)測(cè)的總體偏離程度,避免結(jié)果過(guò)于保守。不確定性集合的相關(guān)研究詳見(jiàn)文獻(xiàn)[12,13],本文不再贅述。

    2 仿射可調(diào)整魯棒機(jī)組組合模型

    在可調(diào)整魯棒機(jī)組組合模型中,機(jī)組出力建模為可調(diào)整變量,其與不確定參數(shù)之間的函數(shù)形式y(tǒng)(d)是由優(yōu)化問(wèn)題隱式?jīng)Q定的,該模型難以直接求解。為此,本文采用仿射函數(shù)建立可調(diào)整變量與其所依賴的不確定參數(shù)之間的仿射關(guān)系,從而使模型變得易于求解。

    在包含風(fēng)電的發(fā)電調(diào)度過(guò)程中,基于電力平衡原則,只有當(dāng)不確定性風(fēng)電的實(shí)際值已知之后,通過(guò)實(shí)時(shí)有功調(diào)度,常規(guī)機(jī)組的有功出力才能最終確定。因此,在常規(guī)機(jī)組可調(diào)容量充裕的情況下,可以采用合理的假設(shè),即假設(shè)常規(guī)機(jī)組出力可以根據(jù)對(duì)應(yīng)的風(fēng)電出力進(jìn)行自動(dòng)調(diào)整。為此,引入常規(guī)機(jī)組出力與其所依賴的風(fēng)電出力之間的反饋規(guī)則,設(shè)定常規(guī)機(jī)組出力為相應(yīng)不確定風(fēng)電的仿射函數(shù)。

    線性仿射函數(shù)預(yù)先假設(shè)可調(diào)整變量與不確定參數(shù)之間線性相關(guān),又稱為線性決策規(guī)則,主要用于運(yùn)籌學(xué)領(lǐng)域,其理論與實(shí)例介紹詳見(jiàn)專著[17]。線性仿射函數(shù)下的機(jī)組出力不再是完全可調(diào)整的,而是限于由仿射函數(shù)提供的調(diào)整范圍。雖然相對(duì)自由化仿射函數(shù)而言無(wú)法獲得最優(yōu)的結(jié)果,然而自由化仿射函數(shù)通常難以求解,限制了其實(shí)際應(yīng)用。采用線性仿射函數(shù)的魯棒機(jī)組組合模型是一個(gè)線性優(yōu)化問(wèn)題,相對(duì)于非線性問(wèn)題而言易于求解。在不確定機(jī)組組合問(wèn)題中,機(jī)組出力建模為不確定參數(shù)的線性仿射函數(shù),隨著不確定參數(shù)的變化而進(jìn)行調(diào)整。仿射函數(shù)定義如下:

    (6)

    式中:K為系統(tǒng)中風(fēng)電場(chǎng)的個(gè)數(shù);dk為風(fēng)電場(chǎng)的有功出力;y0與yk為新引入的中間決策變量,其值通過(guò)求解魯棒優(yōu)化模型確定。在優(yōu)化模型中顯式包含該仿射函數(shù),將可調(diào)整魯棒機(jī)組組合模型轉(zhuǎn)化為一個(gè)單個(gè)階段的優(yōu)化問(wèn)題,轉(zhuǎn)化后的模型具體如下,其中Q為新引入的中間變量。

    (7)

    s.t.Fx≤f

    (8)

    (9)

    (10)

    (11)

    (12)

    式中:Q表示機(jī)組的經(jīng)濟(jì)調(diào)度成本;式(7—12)為機(jī)組組合模型式(1—5)中代入仿射函數(shù)轉(zhuǎn)化而來(lái)。式(7)為總成本最低的魯棒優(yōu)化目標(biāo)函數(shù);式(8)包括啟停狀態(tài)約束、最小開(kāi)停機(jī)時(shí)間約束;式(9)為總成本中的經(jīng)濟(jì)調(diào)度成本約束;式(10)包括備用需求約束、機(jī)組爬坡約束、電網(wǎng)安全約束;式(11)為機(jī)組出力上下限約束;式(12)為系統(tǒng)平衡約束。

    3 模型求解

    在上述仿射可調(diào)整魯棒模型中,由于式(9—12)中非線性因素的存在,無(wú)法直接求解。為此,采用對(duì)偶理論進(jìn)一步對(duì)模型進(jìn)行轉(zhuǎn)化。為便于敘述,假設(shè)不確定集合D為“盒式”不等式[11]:

    D={dk|dk0-Δdk≤dk≤dk0+Δdk}

    (13)

    其中,dk,dk0,Δdk分別為風(fēng)電場(chǎng)k的出力實(shí)際值、波動(dòng)區(qū)間中值與波動(dòng)區(qū)間半徑。分別以不等式(9)、等式(12)為例說(shuō)明模型轉(zhuǎn)化過(guò)程,不等式(10)和(11)的轉(zhuǎn)化原理與式(9)相同,不再贅述。

    為了保證不等式(9)的可行性,即最壞情況下的經(jīng)濟(jì)調(diào)度成本不高于Q,需要確保如下優(yōu)化問(wèn)題的目標(biāo)值小于或等于Q:

    (14)

    s.t.dk0-Δdk≤dk≤dk0+Δd

    (15)

    上述優(yōu)化問(wèn)題符合強(qiáng)對(duì)偶理論的要求,因而可以通過(guò)對(duì)偶問(wèn)題模型來(lái)取代極大值模型,其對(duì)偶問(wèn)題如下:

    (16)

    s.t.rk-sk=bTyk

    (17)

    其中,rk,sk為非負(fù)的對(duì)偶變量。根據(jù)強(qiáng)對(duì)偶理論,當(dāng)且僅當(dāng)對(duì)偶問(wèn)題的目標(biāo)值小于或等于Q時(shí),才能保證原始問(wèn)題的目標(biāo)值小于或等于Q。對(duì)偶問(wèn)題是一個(gè)極小值模型,要求其目標(biāo)值小于或等于Q,等價(jià)于存在rk,sk以滿足要求。因此,式(9)可轉(zhuǎn)化為如下線性化的形式:

    (18)

    rk-sk=bTyk

    (19)

    為了保證系統(tǒng)平衡約束式(12)對(duì)于風(fēng)電場(chǎng)的可能出力dk始終保持可行性,分別將式(12)中的常量項(xiàng)與不確定參數(shù)dk相關(guān)項(xiàng)進(jìn)行歸并,當(dāng)且僅當(dāng)如下等價(jià)得到滿足:

    Ty0-l=0

    (20)

    (21)

    至此,將可調(diào)整魯棒機(jī)組組合問(wèn)題轉(zhuǎn)化為標(biāo)準(zhǔn)的混合整數(shù)規(guī)劃模型,該模型由式(7)、式(8)、式(10)、式(11)、式(18—21)組成,可以采用成熟的商用求解器(如CPLEX)快速求解。

    4 算例分析

    4.1 基本算例

    本文采用改進(jìn)的IEEE RTS—1996算例[18],系統(tǒng)原始數(shù)據(jù)中有32臺(tái)發(fā)電機(jī)組,總裝機(jī)容量為3405 MW。系統(tǒng)中有24個(gè)母線節(jié)點(diǎn)、7個(gè)安全約束斷面,峰荷2850 MW出現(xiàn)在第18 h,第19 h。在該系統(tǒng)中加入風(fēng)電機(jī)組,風(fēng)電通過(guò)第16節(jié)點(diǎn)集中接入,其預(yù)測(cè)出力采用國(guó)內(nèi)某電網(wǎng)的區(qū)域匯集風(fēng)電的實(shí)際數(shù)據(jù)。作為系統(tǒng)中的不確定性因素,設(shè)定風(fēng)電的實(shí)際可能出力在其預(yù)測(cè)值上下20%的波動(dòng)區(qū)間內(nèi)變化,如圖1所示。建立安全約束下的24 h日前機(jī)組組合魯棒模型,采用混合整數(shù)規(guī)劃軟件包CPLEX12.3求解。

    圖1 系統(tǒng)負(fù)荷與風(fēng)功率預(yù)測(cè)曲線Fig.1 Load and wind power forecast curve

    將魯棒優(yōu)化(RO)的結(jié)果與當(dāng)前在電力生產(chǎn)中廣泛應(yīng)用的備用調(diào)整方法(REG)的結(jié)果進(jìn)行對(duì)比,備用調(diào)整方法為考慮新能源備用需求的確定性發(fā)電計(jì)劃優(yōu)化。本文將其備用需求量設(shè)定為與風(fēng)電的波動(dòng)量一致。

    4.2 魯棒性分析

    對(duì)大量模擬場(chǎng)景進(jìn)行經(jīng)濟(jì)調(diào)度(ED)計(jì)算,以評(píng)估機(jī)組組合結(jié)果對(duì)波動(dòng)區(qū)間內(nèi)的任意場(chǎng)景是否都能夠保證約束可行性。仿真結(jié)果如表1所示。

    表1 2種算法的ED收斂次數(shù)比較Table 1 Convergence results of ED between two methods

    表1中,RO獲得的機(jī)組組合結(jié)果顯示,計(jì)算10 000個(gè)場(chǎng)景下的經(jīng)濟(jì)調(diào)度問(wèn)題全部收斂,表明該方法的開(kāi)停機(jī)結(jié)果對(duì)于風(fēng)電的隨機(jī)波動(dòng)具有較高的適應(yīng)性。REG得到的機(jī)組組合結(jié)果卻在多個(gè)場(chǎng)景下經(jīng)濟(jì)調(diào)度無(wú)法收斂。

    RO建模考慮了不確定集合內(nèi)的所有可能的情況,機(jī)組組合結(jié)果能夠適應(yīng)風(fēng)電在區(qū)間內(nèi)的任意波動(dòng)變,在該場(chǎng)景下的經(jīng)濟(jì)調(diào)度問(wèn)題是收斂的。

    4.3 經(jīng)濟(jì)性分析

    在經(jīng)濟(jì)調(diào)度模型中引入松弛變量,以保證經(jīng)濟(jì)性分析中優(yōu)化問(wèn)題的收斂性,并設(shè)定松弛變量的懲罰因子為1000 $/(MW·h),算例中機(jī)組平均發(fā)電成本為10 $/(MW·h)。

    RO的特點(diǎn)是考慮了不確定區(qū)間內(nèi)的所有可能情況。為此,對(duì)于機(jī)組組合結(jié)果,分別計(jì)算10 000個(gè)模擬場(chǎng)景下的經(jīng)濟(jì)調(diào)度問(wèn)題,平均成本對(duì)比如表2所示。

    表2 2種算法的ED成本比較Table 2 Cost comparison of ED between two methods $

    由表2可知,基于機(jī)組組合結(jié)果計(jì)算10 000個(gè)模擬場(chǎng)景下的經(jīng)濟(jì)調(diào)度問(wèn)題,RO結(jié)果的經(jīng)濟(jì)調(diào)度平均總成本要低于REG,因?yàn)槠淇紤]了多場(chǎng)景的可能性,不會(huì)產(chǎn)生經(jīng)濟(jì)調(diào)度懲罰成本,從而總體具有較好的經(jīng)濟(jì)效益。

    4.4 計(jì)算性能分析

    測(cè)試算例的可調(diào)整魯棒優(yōu)化模型共有27 000多個(gè)約束條件,14 000多個(gè)決策變量,其中整數(shù)變量為2304個(gè)?;旌险麛?shù)規(guī)劃算法的收斂精度設(shè)置為0.10%。問(wèn)題的求解時(shí)間為26.8 s,優(yōu)化目標(biāo)為459 746.76 $。將本文方法與兩階段分解算法[15]的計(jì)算性能進(jìn)行對(duì)比,如表3所示。

    表3 2種算法的計(jì)算性能比較Table 3 Performance comparison between two methods

    由表3可知,本文方法在計(jì)算效率和計(jì)算結(jié)果上均優(yōu)于兩階段分解算法,顯示出良好的計(jì)算性能。

    5 結(jié)語(yǔ)

    本文針對(duì)間歇式能源接入后電力系統(tǒng)調(diào)度決策環(huán)境不確定性的問(wèn)題,提出了基于可調(diào)整優(yōu)化的魯棒機(jī)組組合求解方法。算例仿真表明魯棒機(jī)組組合結(jié)果可以適應(yīng)風(fēng)電在波動(dòng)區(qū)間內(nèi)的任意變化,降低實(shí)時(shí)調(diào)度的風(fēng)險(xiǎn),提升系統(tǒng)運(yùn)行的整體經(jīng)濟(jì)效益。同時(shí),根據(jù)調(diào)度運(yùn)行需求,通過(guò)設(shè)置合理的不確定集合,可以實(shí)現(xiàn)電網(wǎng)安全與經(jīng)濟(jì)的最優(yōu)協(xié)調(diào)。

    [1] 裴哲義,董 存,辛耀中. 我國(guó)風(fēng)電并網(wǎng)運(yùn)行最新進(jìn)展[J]. 中國(guó)電力,2010,43(6):78-82. PEI Zheyi,DONG Cun,XIN Yaozhong. Review of operation and management of integrating wind power in China[J]. Electric power,2010,43(6):78-82.

    [2] UMMELS B C,GIBESCU M,PELGRUM E,et al. Impacts of wind power on thermal generation unit commitment and dispatch[J]. IEEE Transaction on Energy Conversion,2007,22(1):44-51.

    [3] 李利利,姚建國(guó),耿 建. SCUC/SCED問(wèn)題分析[J]. 江蘇電機(jī)工程,2010,29(3):24-27. LI Lili,YAO Jianguo,GENG Jian.Analysis of the research on SCUC/SCED[J]. Jiangsu Electrical Engineering,2010,29(3):24-27.

    [4] 張文婷,范立新,顧 文. 考慮風(fēng)電預(yù)測(cè)誤差的電力系統(tǒng)有功實(shí)時(shí)調(diào)度[J]. 江蘇電機(jī)工程,2016,35(1):11-15. ZHANG Wenting,F(xiàn)AN Lixin,GU Wen. Real-time active power dispatch of wind integrated system considering wind power forecast reeors[J]. Jiangsu Electrical Engineering,2016,35(1):11-15.

    [5] 趙晉泉,汪 晶. 不確定機(jī)組組合處理方法評(píng)述[J]. 江蘇電機(jī)工程,2010,29(5):1-5. ZHAO Jinquan,WANG Jing. Review on handling methods for uncertain unit commitment[J]. Jiangsu Electrical Engineering,2010,29(5):1-5.

    [6] BEN-TAL A,NEMIROVSKI A.Robust optimization: methodology and applications[J]. Mathematical Programming,2002,92(3):453-480.

    [7] STREET A,OLIVEIRA F,ARROYO J. Contingency constrained unit commitment with n-K security criterion: a robust optimization approach[J]. IEEE Transactions on Power Systems,2011,26(3):1581-1590.

    [8] BARINGO L,CONEJO A.Offering strategy via robust optimization[J]. IEEE Transactions on Power Systems,2011, 26(3):1418-1425.

    [9] HAJIMIRAGHA A,CANIZARES C,F(xiàn)OWLER M. A robust optimization approach for planning the transition to plug-in hybrid electric vehicles[J]. IEEE Transactions on Power Systems,2011,26(4):2264-2274.

    [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] 魏 韡,劉 鋒,梅生偉. 電力系統(tǒng)魯棒經(jīng)濟(jì)調(diào)度:(一)理論基礎(chǔ)[J]. 電力系統(tǒng)自動(dòng)化,2013,37(17):37-43. WEI Wei, LIU Feng,MEI Shengwei.Robust and economical scheduling methodology for power systems Part one theoretical foundation[J]. Automation of Electric Power Systems,2013,37(17):37-43.

    [12] BERTSIMAS D, LITVINOV E, SUN A X. Adaptive robust optimization for the security constrained unit commitment problem[J]. IEEE Transactions on Power Systems,2013,28(1):52-63.

    [13] JIANG R,WANG J,GUAN Y. Robust unit commitment with wind power and pumped storage hydro[J]. IEEE Transactions on Power Systems, 2012,27(2):800-810.

    [14] WANG Q, WATSON J, GUAN Y. Two-stage robust optimization for N-k contingency-constrained unit commitment[J]. IEEE Transactions on Power Systems, 2013,28(3):2366-2375.

    [15] ZHAO C,WANG J,WATSON J,et al. Multi-stage robust unit commitment considering wind and demand response uncertainties[J]. IEEE Transactions on Power Systems, 2013,28(3):2708-2717.

    [16] BEN-TALl A,GORVASHKO A,GUSLITZER E. Adjustable robust solutions of uncertain linear programs[J]. Mathematical Programming,2004,99(2):351-376.

    [17] BEN-TAL A,GHAOUI L,NEMIROVSKI A.Robust optimization[M].New Jersey: Princeton University Press,2009.

    [18] Reliability Test System Task Force. The IEEE reliability test system—1996[J]. IEEE Transactions on Power Systems,1999,27(2):800-810.

    李利利

    李利利(1987 —),男,安徽宿州人,高級(jí)工程師,從事電力系統(tǒng)運(yùn)行和優(yōu)化研究工作;

    丁 恰(1974 —),男,江蘇南京人,研究員級(jí)高級(jí)工程師,從事電力系統(tǒng)自動(dòng)化研究工作;

    涂孟夫(1977 —),男,湖南長(zhǎng)沙人,高級(jí)工程師,從事電力系統(tǒng)運(yùn)行和優(yōu)化研究工作;

    謝麗榮(1977 —),女,河北衡水人,高級(jí)工程師,從事電力系統(tǒng)運(yùn)行和優(yōu)化研究工作。

    (編輯 劉曉燕)

    Affinely Adjustable Robust Optimization Model andAlgorithm for Unit Commitment Problem

    LI Lili, DING Qia, TU Mengfu, XIE Lirong

    (NARI Group Corporation State Grid Electric Power Research Institute, Nanjing 211106, China)

    A unit commitment optimization method based on affinely adjustable robust optimization is proposed for power system scheduling with uncertainty, especially when large-scale intermittent energy integrated. An affinely adjustable robust optimal model of uncertainty unit commitment is developed. The affine relationship between decision variables and uncertain parameters is established by linear decision rules. Therefore, the two-stage model can be reformulated into a single stage optimization problem. On this basis, the problem can be transformed to a standard mixed integer programming model through duality theory. Analysis of test system shows the effectiveness of the proposed algorithm.

    generation scheduling; uncertainty; robust optimization; security constrained unit commitment (SCUC); linear decision rules

    2016-11-30;

    2017-01-15

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51577031);國(guó)家電網(wǎng)公司科技項(xiàng)目(考慮多基地和直流外送的新能源消納調(diào)度決策分析和評(píng)價(jià)方法研究及應(yīng)用)

    TM734

    A

    2096-3203(2017)03-0033-05

    猜你喜歡
    魯棒出力約束
    “碳中和”約束下的路徑選擇
    約束離散KP方程族的完全Virasoro對(duì)稱
    基于學(xué)習(xí)的魯棒自適應(yīng)評(píng)判控制研究進(jìn)展
    目標(biāo)魯棒識(shí)別的抗旋轉(zhuǎn)HDO 局部特征描述
    風(fēng)電場(chǎng)有功出力的EEMD特性分析
    要爭(zhēng)做出力出彩的黨員干部
    河南電力(2016年5期)2016-02-06 02:11:35
    基于Cauchy魯棒函數(shù)的UKF改進(jìn)算法
    風(fēng)電場(chǎng)群出力的匯聚效應(yīng)分析
    目標(biāo)軌跡更新的點(diǎn)到點(diǎn)魯棒迭代學(xué)習(xí)控制
    適當(dāng)放手能讓孩子更好地自我約束
    人生十六七(2015年6期)2015-02-28 13:08:38
    少妇的逼水好多| 3wmmmm亚洲av在线观看| 国产精品,欧美在线| 欧美日韩乱码在线| 少妇人妻精品综合一区二区 | 免费高清视频大片| 欧美黑人巨大hd| 国产v大片淫在线免费观看| 少妇的逼水好多| 日韩欧美 国产精品| 国产精品免费一区二区三区在线| 99riav亚洲国产免费| 成熟少妇高潮喷水视频| 尾随美女入室| 日本五十路高清| 国产免费av片在线观看野外av| 国产高清有码在线观看视频| 精品午夜福利视频在线观看一区| 欧美日韩瑟瑟在线播放| 五月伊人婷婷丁香| 亚洲七黄色美女视频| x7x7x7水蜜桃| netflix在线观看网站| 91久久精品电影网| 久久久久久久午夜电影| 一级黄片播放器| 99精品在免费线老司机午夜| ponron亚洲| 亚洲国产欧洲综合997久久,| 久9热在线精品视频| 久久久久久久亚洲中文字幕| 1000部很黄的大片| 99在线人妻在线中文字幕| 国产亚洲av嫩草精品影院| 精品人妻一区二区三区麻豆 | 日韩欧美在线乱码| 女同久久另类99精品国产91| 亚洲美女黄片视频| 99热网站在线观看| 国产亚洲精品久久久com| 最近最新中文字幕大全电影3| 国语自产精品视频在线第100页| 成人三级黄色视频| 日本爱情动作片www.在线观看 | 不卡视频在线观看欧美| 在线观看一区二区三区| 中国美女看黄片| 亚洲精品乱码久久久v下载方式| 亚洲最大成人中文| 在线观看一区二区三区| 18禁裸乳无遮挡免费网站照片| 色尼玛亚洲综合影院| 老熟妇乱子伦视频在线观看| 午夜福利欧美成人| 在线观看午夜福利视频| 午夜精品一区二区三区免费看| 亚洲一区二区三区色噜噜| 丰满人妻一区二区三区视频av| 久久亚洲真实| 日韩欧美国产在线观看| 国产午夜福利久久久久久| 国产精品,欧美在线| 亚洲va在线va天堂va国产| 国产av一区在线观看免费| 亚洲精品一卡2卡三卡4卡5卡| 精品久久久久久久久亚洲 | 在线观看66精品国产| 欧美bdsm另类| 18禁裸乳无遮挡免费网站照片| 国产 一区 欧美 日韩| 国产老妇女一区| 国产精品,欧美在线| 免费观看的影片在线观看| 乱码一卡2卡4卡精品| 麻豆成人av在线观看| 丰满的人妻完整版| 精品欧美国产一区二区三| 黄色欧美视频在线观看| 国产三级中文精品| 99国产精品一区二区蜜桃av| 欧美日韩亚洲国产一区二区在线观看| 麻豆成人av在线观看| 免费电影在线观看免费观看| bbb黄色大片| 午夜免费男女啪啪视频观看 | 嫩草影院精品99| 啦啦啦啦在线视频资源| 一个人看的www免费观看视频| 色综合亚洲欧美另类图片| 午夜福利在线观看吧| 久久久久精品国产欧美久久久| 中文亚洲av片在线观看爽| 色吧在线观看| 乱系列少妇在线播放| 日韩一区二区视频免费看| 大型黄色视频在线免费观看| 校园人妻丝袜中文字幕| 中文字幕人妻熟人妻熟丝袜美| 国产精品亚洲一级av第二区| 欧美黑人欧美精品刺激| 亚洲成人久久爱视频| 91在线观看av| 日韩欧美精品v在线| 日本免费a在线| 如何舔出高潮| 永久网站在线| 免费大片18禁| 国产又黄又爽又无遮挡在线| 中文字幕av成人在线电影| 综合色av麻豆| 亚洲avbb在线观看| 久久国产精品人妻蜜桃| 中文亚洲av片在线观看爽| 国产精品一区二区三区四区久久| 能在线免费观看的黄片| 欧美日韩瑟瑟在线播放| 亚洲国产欧洲综合997久久,| 三级毛片av免费| 日韩欧美在线乱码| 久久婷婷人人爽人人干人人爱| 18禁裸乳无遮挡免费网站照片| 国内精品久久久久精免费| 久久久久久九九精品二区国产| 噜噜噜噜噜久久久久久91| 精品欧美国产一区二区三| 精品久久久久久久人妻蜜臀av| 久久久午夜欧美精品| а√天堂www在线а√下载| 国内久久婷婷六月综合欲色啪| 国产熟女欧美一区二区| 亚洲狠狠婷婷综合久久图片| 欧洲精品卡2卡3卡4卡5卡区| 一区二区三区高清视频在线| 亚洲精品国产成人久久av| 日本黄色片子视频| 日本黄色视频三级网站网址| 免费看美女性在线毛片视频| 国产精品永久免费网站| 中文字幕av成人在线电影| 观看美女的网站| 内地一区二区视频在线| 国内少妇人妻偷人精品xxx网站| 深爱激情五月婷婷| 精品人妻一区二区三区麻豆 | 国产午夜精品久久久久久一区二区三区 | 久久精品国产自在天天线| 国产精品久久久久久av不卡| 啦啦啦韩国在线观看视频| 亚洲真实伦在线观看| 亚洲七黄色美女视频| 国产又黄又爽又无遮挡在线| 亚洲 国产 在线| 欧美日韩乱码在线| 联通29元200g的流量卡| 午夜影院日韩av| 亚洲aⅴ乱码一区二区在线播放| 中文字幕人妻熟人妻熟丝袜美| 一区二区三区高清视频在线| 夜夜爽天天搞| 亚洲午夜理论影院| 国产成年人精品一区二区| 亚洲欧美激情综合另类| 在线天堂最新版资源| a级毛片免费高清观看在线播放| 国产精品伦人一区二区| 一区二区三区高清视频在线| 国产精品女同一区二区软件 | 亚洲熟妇熟女久久| 校园春色视频在线观看| 日韩人妻高清精品专区| 欧美黑人巨大hd| 国内精品久久久久精免费| 日本与韩国留学比较| 中亚洲国语对白在线视频| 国产色婷婷99| 婷婷六月久久综合丁香| 成年免费大片在线观看| 亚洲成人久久性| 在线观看美女被高潮喷水网站| 国产精品永久免费网站| 免费看a级黄色片| 毛片一级片免费看久久久久 | 国产精品一及| 亚洲欧美精品综合久久99| 亚洲中文字幕日韩| 国产一区二区亚洲精品在线观看| 日韩欧美国产在线观看| 精华霜和精华液先用哪个| 久久九九热精品免费| 色5月婷婷丁香| 久久久久久久久大av| 能在线免费观看的黄片| 亚洲人成网站在线播放欧美日韩| 黄色视频,在线免费观看| 久久精品国产鲁丝片午夜精品 | 成人亚洲精品av一区二区| 国产精品一区二区性色av| 18禁在线播放成人免费| 黄色欧美视频在线观看| 国产高清激情床上av| 熟女人妻精品中文字幕| 白带黄色成豆腐渣| 亚洲性夜色夜夜综合| 亚洲人成网站在线播| 国产爱豆传媒在线观看| 国产精品久久久久久av不卡| 成人永久免费在线观看视频| 女生性感内裤真人,穿戴方法视频| 老司机福利观看| 午夜福利在线观看免费完整高清在 | 国产三级中文精品| 成年人黄色毛片网站| 色5月婷婷丁香| 99热只有精品国产| 欧美极品一区二区三区四区| 此物有八面人人有两片| videossex国产| 麻豆一二三区av精品| 成人国产麻豆网| 日韩亚洲欧美综合| av天堂中文字幕网| 校园人妻丝袜中文字幕| 在线国产一区二区在线| 免费电影在线观看免费观看| 亚洲avbb在线观看| 嫩草影视91久久| 97热精品久久久久久| 欧美日本视频| 国国产精品蜜臀av免费| 久久久久久久久久久丰满 | 日韩亚洲欧美综合| 精品不卡国产一区二区三区| а√天堂www在线а√下载| 欧美最黄视频在线播放免费| 国产精品久久久久久精品电影| 无人区码免费观看不卡| 色哟哟哟哟哟哟| 丝袜美腿在线中文| 搡老熟女国产l中国老女人| 亚州av有码| 婷婷丁香在线五月| 日韩精品青青久久久久久| www.www免费av| 亚洲国产日韩欧美精品在线观看| АⅤ资源中文在线天堂| 日韩欧美国产在线观看| 日韩精品有码人妻一区| 夜夜爽天天搞| 九九爱精品视频在线观看| 日日摸夜夜添夜夜添av毛片 | 亚洲av日韩精品久久久久久密| av女优亚洲男人天堂| 可以在线观看毛片的网站| 91午夜精品亚洲一区二区三区 | 亚洲天堂国产精品一区在线| 国产精品综合久久久久久久免费| 国内揄拍国产精品人妻在线| 免费搜索国产男女视频| 欧美不卡视频在线免费观看| av国产免费在线观看| а√天堂www在线а√下载| 日韩欧美国产在线观看| 亚洲电影在线观看av| 床上黄色一级片| 在线看三级毛片| 欧美绝顶高潮抽搐喷水| 婷婷色综合大香蕉| 色精品久久人妻99蜜桃| 日本三级黄在线观看| 免费av不卡在线播放| 少妇丰满av| 97热精品久久久久久| 国产色婷婷99| 欧美又色又爽又黄视频| 亚洲精品粉嫩美女一区| 国产亚洲精品久久久com| 国产精品野战在线观看| 国产精品久久久久久亚洲av鲁大| 乱人视频在线观看| 国产v大片淫在线免费观看| 天天一区二区日本电影三级| 国产精品一及| 无人区码免费观看不卡| 老司机深夜福利视频在线观看| 97人妻精品一区二区三区麻豆| 成人特级av手机在线观看| 国产真实乱freesex| 国产伦精品一区二区三区四那| 成人av一区二区三区在线看| 在线观看午夜福利视频| 免费av不卡在线播放| 一夜夜www| 麻豆一二三区av精品| 亚州av有码| 国产久久久一区二区三区| 亚洲自拍偷在线| 精品人妻熟女av久视频| 人妻夜夜爽99麻豆av| 亚洲精品久久国产高清桃花| 精品福利观看| 欧美zozozo另类| netflix在线观看网站| 少妇的逼好多水| 国产在视频线在精品| 成人三级黄色视频| 又爽又黄无遮挡网站| 嫩草影视91久久| 最近在线观看免费完整版| h日本视频在线播放| 男女之事视频高清在线观看| 欧美日韩国产亚洲二区| 美女xxoo啪啪120秒动态图| 午夜激情福利司机影院| 小蜜桃在线观看免费完整版高清| 欧美日韩综合久久久久久 | 91久久精品国产一区二区成人| 如何舔出高潮| 国产一区二区三区视频了| 亚洲精品久久国产高清桃花| 在线a可以看的网站| 亚洲av第一区精品v没综合| 日韩欧美在线乱码| 最近视频中文字幕2019在线8| 成人综合一区亚洲| 变态另类丝袜制服| 午夜精品一区二区三区免费看| 成人国产麻豆网| 变态另类成人亚洲欧美熟女| 国产精品久久久久久精品电影| 天堂动漫精品| 高清在线国产一区| av在线观看视频网站免费| 国产精品久久久久久av不卡| 在线播放无遮挡| 精华霜和精华液先用哪个| 久久国内精品自在自线图片| 老司机深夜福利视频在线观看| 麻豆国产97在线/欧美| 乱码一卡2卡4卡精品| 国产一区二区三区av在线 | 网址你懂的国产日韩在线| 男女做爰动态图高潮gif福利片| 日韩欧美免费精品| 能在线免费观看的黄片| 国产aⅴ精品一区二区三区波| 波野结衣二区三区在线| 精品久久久久久久久av| 国内精品久久久久精免费| 日韩欧美一区二区三区在线观看| 久久精品国产鲁丝片午夜精品 | 亚洲精品乱码久久久v下载方式| 美女大奶头视频| 亚洲欧美激情综合另类| 国产午夜精品久久久久久一区二区三区 | 国产中年淑女户外野战色| 高清在线国产一区| 成人国产综合亚洲| 色在线成人网| 久久久久久大精品| 免费观看的影片在线观看| 一卡2卡三卡四卡精品乱码亚洲| 中国美白少妇内射xxxbb| 一级黄片播放器| 免费在线观看日本一区| 国产高清不卡午夜福利| www日本黄色视频网| 黄色视频,在线免费观看| 精品人妻熟女av久视频| av在线老鸭窝| 亚洲成人免费电影在线观看| 色噜噜av男人的天堂激情| 欧美性猛交╳xxx乱大交人| 一本一本综合久久| 中文字幕人妻熟人妻熟丝袜美| 日韩在线高清观看一区二区三区 | 国产精品福利在线免费观看| 成人特级av手机在线观看| 日本一二三区视频观看| 日本在线视频免费播放| 久久久久国产精品人妻aⅴ院| 99国产极品粉嫩在线观看| 国产成人福利小说| 美女cb高潮喷水在线观看| 99热这里只有精品一区| 91狼人影院| 国产精品久久久久久亚洲av鲁大| 亚洲熟妇熟女久久| 1024手机看黄色片| 国产 一区 欧美 日韩| 在线天堂最新版资源| 男女下面进入的视频免费午夜| 一级毛片久久久久久久久女| 啦啦啦啦在线视频资源| 搡老熟女国产l中国老女人| 人人妻,人人澡人人爽秒播| 床上黄色一级片| 日韩亚洲欧美综合| 亚洲在线观看片| 国产老妇女一区| 亚洲精品色激情综合| 在线免费观看不下载黄p国产 | 精品久久国产蜜桃| 亚洲av不卡在线观看| 伦理电影大哥的女人| 亚洲成人中文字幕在线播放| 性色avwww在线观看| 免费在线观看日本一区| 一级a爱片免费观看的视频| 麻豆国产av国片精品| 日韩国内少妇激情av| 变态另类丝袜制服| 99热这里只有精品一区| 日韩一本色道免费dvd| 欧美zozozo另类| 日韩欧美在线乱码| 久久婷婷人人爽人人干人人爱| 免费人成视频x8x8入口观看| 国产高清三级在线| 久久精品人妻少妇| 91久久精品国产一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 此物有八面人人有两片| 在现免费观看毛片| 日韩一本色道免费dvd| 欧美色欧美亚洲另类二区| 深夜a级毛片| 内地一区二区视频在线| 在线免费观看的www视频| 99热精品在线国产| 久久久久久久午夜电影| 搞女人的毛片| 毛片一级片免费看久久久久 | 精品免费久久久久久久清纯| 国产午夜精品论理片| 成人特级av手机在线观看| 国产av一区在线观看免费| 亚洲成人精品中文字幕电影| 成年人黄色毛片网站| 看十八女毛片水多多多| 日本免费a在线| 亚洲精品乱码久久久v下载方式| 亚洲成人精品中文字幕电影| 狂野欧美白嫩少妇大欣赏| 午夜福利在线观看免费完整高清在 | 免费高清视频大片| 日本免费a在线| 国产伦一二天堂av在线观看| 国国产精品蜜臀av免费| 日韩欧美国产在线观看| 18禁黄网站禁片免费观看直播| 麻豆av噜噜一区二区三区| 亚洲专区中文字幕在线| 午夜精品在线福利| 免费av不卡在线播放| 亚洲欧美日韩卡通动漫| 可以在线观看毛片的网站| 国产精品久久久久久精品电影| 校园春色视频在线观看| 午夜影院日韩av| 国产美女午夜福利| 永久网站在线| 最近最新免费中文字幕在线| 黄色一级大片看看| 午夜亚洲福利在线播放| 中文字幕av在线有码专区| 十八禁网站免费在线| 乱码一卡2卡4卡精品| 男女做爰动态图高潮gif福利片| 成人亚洲精品av一区二区| 亚洲成人精品中文字幕电影| 韩国av在线不卡| 国产三级中文精品| 人妻少妇偷人精品九色| 日本黄色视频三级网站网址| 国产精华一区二区三区| 又粗又爽又猛毛片免费看| 国产精品伦人一区二区| 亚洲av.av天堂| 美女被艹到高潮喷水动态| 午夜亚洲福利在线播放| 一区二区三区激情视频| 久久精品国产亚洲网站| 精品国产三级普通话版| 欧美潮喷喷水| 国产伦人伦偷精品视频| 日本成人三级电影网站| 狂野欧美白嫩少妇大欣赏| 婷婷色综合大香蕉| 真实男女啪啪啪动态图| 最近最新中文字幕大全电影3| 久久精品综合一区二区三区| 天天一区二区日本电影三级| 亚洲最大成人手机在线| 亚洲美女黄片视频| 国产 一区 欧美 日韩| 日韩强制内射视频| 亚洲av五月六月丁香网| 中文字幕精品亚洲无线码一区| 国产午夜精品久久久久久一区二区三区 | 少妇猛男粗大的猛烈进出视频 | 在线a可以看的网站| 简卡轻食公司| 美女免费视频网站| 精品一区二区三区av网在线观看| 99九九线精品视频在线观看视频| 国产精品1区2区在线观看.| 好男人在线观看高清免费视频| 国产午夜精品久久久久久一区二区三区 | 国产一区二区激情短视频| 国产成人一区二区在线| 日韩欧美精品免费久久| 嫩草影视91久久| 中文字幕免费在线视频6| 亚洲18禁久久av| 国产精品99久久久久久久久| 久久国内精品自在自线图片| 国产一区二区三区av在线 | 久久精品国产亚洲av天美| 麻豆国产97在线/欧美| 久久久久国产精品人妻aⅴ院| 国产精品嫩草影院av在线观看 | 日韩 亚洲 欧美在线| www日本黄色视频网| 日韩大尺度精品在线看网址| 欧美日韩乱码在线| 国产一区二区三区视频了| 色视频www国产| 色综合亚洲欧美另类图片| 18+在线观看网站| 夜夜爽天天搞| 日本免费一区二区三区高清不卡| 精品久久久久久成人av| 丰满乱子伦码专区| 亚洲经典国产精华液单| 国产美女午夜福利| 日韩欧美一区二区三区在线观看| 亚洲精华国产精华液的使用体验 | 欧美激情在线99| 身体一侧抽搐| 蜜桃亚洲精品一区二区三区| 国产成年人精品一区二区| 久久久久免费精品人妻一区二区| 日韩欧美精品免费久久| 国产精品福利在线免费观看| 99热这里只有是精品50| 国产高清有码在线观看视频| av在线蜜桃| 中文字幕免费在线视频6| 成年女人永久免费观看视频| 波野结衣二区三区在线| 日本精品一区二区三区蜜桃| 亚洲va日本ⅴa欧美va伊人久久| 亚洲最大成人手机在线| 日韩中字成人| 大又大粗又爽又黄少妇毛片口| 精品免费久久久久久久清纯| 成人av一区二区三区在线看| 欧美极品一区二区三区四区| 有码 亚洲区| 精品福利观看| 亚洲人成伊人成综合网2020| 日日撸夜夜添| 国产伦在线观看视频一区| 午夜福利视频1000在线观看| 亚洲成人精品中文字幕电影| 大又大粗又爽又黄少妇毛片口| 亚洲国产精品成人综合色| 国产乱人视频| 欧美国产日韩亚洲一区| 国产精品精品国产色婷婷| 中文亚洲av片在线观看爽| 桃色一区二区三区在线观看| 深夜精品福利| 亚洲va日本ⅴa欧美va伊人久久| 成人性生交大片免费视频hd| 久久久成人免费电影| 免费av观看视频| 桃红色精品国产亚洲av| 精品久久久久久,| 国内精品久久久久精免费| 夜夜爽天天搞| 日本与韩国留学比较| 亚洲av成人精品一区久久| 最近中文字幕高清免费大全6 | 亚洲最大成人手机在线| 少妇丰满av| 国产精品人妻久久久久久| 99久久精品国产国产毛片| 国产不卡一卡二| 免费一级毛片在线播放高清视频| 性插视频无遮挡在线免费观看| 日韩欧美国产一区二区入口| 亚洲天堂国产精品一区在线| 日韩欧美免费精品| 午夜福利高清视频| 又爽又黄a免费视频| 一级毛片久久久久久久久女| 91久久精品电影网| 日韩欧美国产在线观看| 日本a在线网址| 精品人妻偷拍中文字幕| 国产在视频线在精品| 日本a在线网址| 日韩精品青青久久久久久| 亚洲av成人精品一区久久| 色播亚洲综合网| 国产三级在线视频| 成人一区二区视频在线观看| 欧美黑人巨大hd| 国产真实伦视频高清在线观看 | 狠狠狠狠99中文字幕| 亚洲专区中文字幕在线| 国产高清三级在线| 一个人观看的视频www高清免费观看| 久久欧美精品欧美久久欧美|