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

    基于風(fēng)險(xiǎn)的多區(qū)互聯(lián)電力系統(tǒng)分布式魯棒動(dòng)態(tài)經(jīng)濟(jì)調(diào)度

    2021-12-12 07:57:38王櫓裕郭創(chuàng)新郭逸豪王勝寒
    電力系統(tǒng)自動(dòng)化 2021年23期
    關(guān)鍵詞:受端聯(lián)絡(luò)線魯棒

    陳 哲,王櫓裕,郭創(chuàng)新,楊 杰,郭逸豪,王勝寒

    (1. 國網(wǎng)浙江省電力有限公司電力科學(xué)研究院,浙江省杭州市 310014;2. 浙江大學(xué)電氣工程學(xué)院,浙江省杭州市 310027)

    0 引言

    隨著電網(wǎng)規(guī)模的不斷擴(kuò)大,現(xiàn)代電力系統(tǒng)通常由多個(gè)子系統(tǒng)互聯(lián)而成。由于中國新能源發(fā)展與負(fù)荷發(fā)展不平衡,大量電力通過跨區(qū)域聯(lián)絡(luò)線從西北地區(qū)傳輸至華東地區(qū)。以中國上海市為例,2018 年接受的外部電力占當(dāng)?shù)刈畲筘?fù)荷的比例已超過50%,7 月和8 月甚至達(dá)到約60%。隨著送端子系統(tǒng)與受端子系統(tǒng)之間的聯(lián)系越來越緊密,為提高系統(tǒng)運(yùn)行整體的經(jīng)濟(jì)性與可靠性,應(yīng)協(xié)調(diào)調(diào)度各區(qū)域電網(wǎng)的發(fā)電資源。文獻(xiàn)[1-3]求解了多區(qū)互聯(lián)電網(wǎng)最優(yōu)潮流問題。文獻(xiàn)[4-5]提出了多區(qū)域電網(wǎng)經(jīng)濟(jì)調(diào)度方法。文獻(xiàn)[6-9]研究了互聯(lián)電網(wǎng)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度問題,優(yōu)化了發(fā)電與備用計(jì)劃。文獻(xiàn)[10-12]協(xié)同優(yōu)化了互聯(lián)電網(wǎng)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度與機(jī)組組合計(jì)劃。以上工作均建立了確定性模型,并按照一定的負(fù)荷比例確定備用容量,以應(yīng)對(duì)電網(wǎng)運(yùn)行所面臨的不確定性因素。然而,隨著以風(fēng)電為代表的可再生能源大規(guī)模并網(wǎng),其高不確定性使得依據(jù)這種簡單規(guī)則獲得的調(diào)度計(jì)劃不經(jīng)濟(jì)或不可靠[13]。為解決這一問題,魯棒優(yōu)化通過不確定集刻畫不確定性因素,并優(yōu)化調(diào)度計(jì)劃以應(yīng)對(duì)各種極端運(yùn)行場(chǎng)景,提升了電力系統(tǒng)運(yùn)行可靠性,故近年來得到了廣泛關(guān)注[14-15]。文獻(xiàn)[16]采用魯棒優(yōu)化模型計(jì)及風(fēng)電不確定性,并優(yōu)化得到了魯棒的多區(qū)互聯(lián)電力系統(tǒng)調(diào)度方案。

    與此同時(shí),由于多區(qū)互聯(lián)電力系統(tǒng)所需調(diào)度的設(shè)備數(shù)目龐大,集中式優(yōu)化調(diào)度模型規(guī)模巨大、難以求解,并且對(duì)通信的極高要求難以滿足。因此,此類問題適合采用分布式方式解決。目前常用的分布式優(yōu)化算法包括拉格朗日松弛法[2]、輔助問題原理法[17]、目 標(biāo) 級(jí) 聯(lián) 分 析 法[4]及 交 替 方 向 乘 子 法(ADMM)[18]等。文獻(xiàn)[19-20]將分布式算法與魯棒優(yōu)化相結(jié)合,對(duì)多區(qū)互聯(lián)電力系統(tǒng)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度問題進(jìn)行了求解。

    盡管國內(nèi)外學(xué)者已經(jīng)針對(duì)多區(qū)互聯(lián)電力系統(tǒng)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度問題開展了廣泛的研究,但目前仍存在3 個(gè)方面問題需要解決。

    1)傳統(tǒng)魯棒優(yōu)化難以確定合適的不確定集。當(dāng)不確定集過大時(shí),調(diào)度策略將會(huì)非常保守;反之,又會(huì)過于激進(jìn)。為解決這一問題,文獻(xiàn)[21]將風(fēng)電概率分布函數(shù)引入兩階段魯棒優(yōu)化,其中不確定集的邊界為優(yōu)化變量而非事先設(shè)定的常數(shù)。文獻(xiàn)[22-24]也在集中式動(dòng)態(tài)經(jīng)濟(jì)調(diào)度問題中應(yīng)用了這一方法。然而,這些研究均要求將預(yù)測(cè)的風(fēng)電出力完全消納,在高比例風(fēng)電并網(wǎng)條件下,這是不經(jīng)濟(jì)甚至不安全的[25]。

    2)由于風(fēng)電出力巨大的不確定性,其實(shí)際出力可能大大偏離預(yù)測(cè)值,而目前的工作僅根據(jù)預(yù)測(cè)風(fēng)電出力安排聯(lián)絡(luò)線功率,不能依據(jù)實(shí)際風(fēng)電出力調(diào)整聯(lián)絡(luò)線傳輸功率。由于需接受大規(guī)模外部電源,對(duì)于本地發(fā)電資源使用不足的受端電網(wǎng)而言,其調(diào)峰困難日益突出。隨著柔性直流技術(shù)的發(fā)展,應(yīng)充分發(fā)揮聯(lián)絡(luò)線功率靈活調(diào)節(jié)的優(yōu)勢(shì),促進(jìn)區(qū)域備用共享。

    3)盡管分布式算法已廣泛用于多區(qū)互聯(lián)電力系統(tǒng)運(yùn)行優(yōu)化問題,但這些方法的收斂性均對(duì)懲罰因子較為敏感。因此,亟待解決參數(shù)選擇的問題。

    基于此,本文提出一種基于風(fēng)險(xiǎn)的兩階段魯棒動(dòng)態(tài)經(jīng)濟(jì)調(diào)度模型。該模型采用可調(diào)不確定集計(jì)及風(fēng)電不確定性,通過優(yōu)化不確定集的邊界實(shí)現(xiàn)運(yùn)行成本與運(yùn)行風(fēng)險(xiǎn)的平衡,同時(shí)可以得到合理的可消納風(fēng)電域。此外,該模型在確定聯(lián)絡(luò)線傳輸功率的同時(shí),也給出了各個(gè)時(shí)間間隔內(nèi)聯(lián)絡(luò)線功率的可調(diào)范圍,以實(shí)現(xiàn)區(qū)域之間備用共享,提高互聯(lián)電網(wǎng)的運(yùn)行經(jīng)濟(jì)性與可靠性。最后,基于ADMM 建立了分布式調(diào)度框架,以減輕區(qū)域間通信負(fù)擔(dān)并保護(hù)區(qū)域內(nèi)信息安全;采用動(dòng)態(tài)乘子更新策略,在分布式求解過程中動(dòng)態(tài)調(diào)整懲罰因子的大小,以解決參數(shù)選擇的難題,同時(shí)提升算法的收斂性與實(shí)用性。

    1 模型建模

    1.1 運(yùn)行風(fēng)險(xiǎn)與可消納風(fēng)電域

    圖1 風(fēng)電出力概率密度曲線Fig.1 Probability density curve of wind power output

    基于CVaR,棄風(fēng)風(fēng)險(xiǎn)與失負(fù)荷風(fēng)險(xiǎn)可分別通過式(1)與式(2)計(jì)算得到。

    1.2 多區(qū)互聯(lián)電力系統(tǒng)分解策略

    圖2 描述了一個(gè)簡單的兩區(qū)域互聯(lián)電力系統(tǒng)。根據(jù)行政調(diào)度區(qū)域?qū)⒃撓到y(tǒng)分為送端子系統(tǒng)和受端子系統(tǒng),送端子系統(tǒng)通過特高壓直流聯(lián)絡(luò)線向受端子系統(tǒng)輸送功率。下文將子系統(tǒng)內(nèi)部稱為區(qū)域內(nèi),將送受端子系統(tǒng)之間稱為區(qū)域間。基于線路撕裂法[27],送端子系統(tǒng)將受端子系統(tǒng)等效為負(fù)荷,受端子系統(tǒng)將送端子系統(tǒng)等效為發(fā)電機(jī)。多區(qū)互聯(lián)電力系統(tǒng)調(diào)度的一個(gè)關(guān)鍵問題是確定聯(lián)絡(luò)線功率,即等效負(fù)荷與發(fā)電機(jī)出力的值。

    圖2 多區(qū)互聯(lián)電力系統(tǒng)Fig.2 Interconnected multi-regional power system

    由于接受大規(guī)模外部電源,受端電網(wǎng)區(qū)域內(nèi)發(fā)電資源使用不足,調(diào)峰困難日益突出。為解決這一問題,本文考慮了區(qū)域間備用互濟(jì),即送端子系統(tǒng)不僅為受端子系統(tǒng)傳輸功率,也同時(shí)通過在風(fēng)電實(shí)際出力場(chǎng)景下確定調(diào)整聯(lián)絡(luò)線傳輸功率的調(diào)整區(qū)間為受端子系統(tǒng)提供備用容量,即聯(lián)絡(luò)線備用。這種情況下,送端子系統(tǒng)將受端子系統(tǒng)等效為不確定性負(fù)荷,受端子系統(tǒng)將送端子系統(tǒng)等效為可調(diào)發(fā)電機(jī)。由于送端子系統(tǒng)發(fā)電與備用成本較受端子系統(tǒng)均更低,因此,協(xié)同兩區(qū)域的發(fā)電與備用資源將取得更好的經(jīng)濟(jì)性。同時(shí)由于備用資源充足,系統(tǒng)運(yùn)行風(fēng)險(xiǎn)也將相應(yīng)降低。

    1.3 基于風(fēng)險(xiǎn)的兩階段魯棒動(dòng)態(tài)經(jīng)濟(jì)調(diào)度模型

    面向多區(qū)互聯(lián)電力系統(tǒng),本文提出一種基于風(fēng)險(xiǎn)的兩階段魯棒動(dòng)態(tài)經(jīng)濟(jì)調(diào)度模型。該模型為日前調(diào)度模型,關(guān)注大規(guī)模風(fēng)電并網(wǎng)后的日內(nèi)調(diào)峰問題,調(diào)度時(shí)間間隔為1 h[28]。同時(shí),由于現(xiàn)有調(diào)度體系中機(jī)組組合方式周前已確定[29],因此本文假定機(jī)組組合方式已知。第1 階段依據(jù)風(fēng)電預(yù)測(cè)出力值安排區(qū)域內(nèi)機(jī)組出力與聯(lián)絡(luò)線傳輸功率計(jì)劃,并確定機(jī)組與聯(lián)絡(luò)線備用容量及可消納風(fēng)電域以應(yīng)對(duì)風(fēng)電不確定性。第2 階段根據(jù)風(fēng)電實(shí)際出力與第1 階段安排的運(yùn)行備用進(jìn)行再調(diào)度,以確保所有在可消納風(fēng)電域內(nèi)的場(chǎng)景均不會(huì)棄風(fēng)與失負(fù)荷。同時(shí),兩個(gè)階段的調(diào)度模型不是孤立的,而是一個(gè)統(tǒng)一的整體。第1 階段制定的調(diào)度計(jì)劃要保證第2 階段最惡劣場(chǎng)景下有足夠的調(diào)整空間來滿足運(yùn)行約束。兩階段問題協(xié)同求解,可以最小化系統(tǒng)運(yùn)行成本、備用成本以及運(yùn)行風(fēng)險(xiǎn),同時(shí)可以得到合理的風(fēng)電消納域。需要說明的是,本節(jié)將給出集中式調(diào)度模型,具體的分布式模型與算法將在2.2 節(jié)中給出。

    1.3.1 第1 階段目標(biāo)函數(shù)

    本文所提模型將多區(qū)互聯(lián)電力系統(tǒng)分為送端子系統(tǒng)與受端子系統(tǒng),其目標(biāo)函數(shù)為最小化各子系統(tǒng)總運(yùn)行成本、備用成本與運(yùn)行風(fēng)險(xiǎn),如式(3)所示。

    本文采用分段線性化函數(shù)描述機(jī)組運(yùn)行成本,如式(4)至式(6)所示。

    式中:B和D分別為聯(lián)絡(luò)線集合和區(qū)域母線集合;Pw,t為時(shí)段t風(fēng)電場(chǎng)w的風(fēng)電消納量;Pb,t為時(shí)段t聯(lián)絡(luò)線b的功率,對(duì)送端子系統(tǒng)相當(dāng)于負(fù)荷,對(duì)受端子系統(tǒng)相當(dāng)于電源;Pd,t為時(shí)段t節(jié)點(diǎn)d的負(fù)荷。

    6)聯(lián)絡(luò)線功率及可調(diào)容量約束

    需要說明的是,由于聯(lián)絡(luò)線調(diào)節(jié)容量的一致性約束,送端子系統(tǒng)可以承受調(diào)整容量內(nèi)受端子系統(tǒng)功率需求的任意波動(dòng),而受端子系統(tǒng)也可以在這個(gè)范圍內(nèi)任意調(diào)節(jié),故實(shí)現(xiàn)了送端與受端子系統(tǒng)間的解耦,同時(shí)也可以獲得所有風(fēng)電出力場(chǎng)景下的最優(yōu)調(diào)整計(jì)劃,保證聯(lián)絡(luò)線的功率平衡。與引入實(shí)際風(fēng)電出力場(chǎng)景下送受端子系統(tǒng)聯(lián)絡(luò)線功率相等的一致性約束相比,本文的建模方法雖然經(jīng)濟(jì)性略有降低,但避免了有關(guān)負(fù)荷和風(fēng)電出力交換的詳細(xì)信息,減輕了通信壓力并提升了求解效率。

    1.3.3 第2 階段目標(biāo)函數(shù)

    第2 階段問題根據(jù)實(shí)際風(fēng)電出力,在第1 階段制定的調(diào)度計(jì)劃基礎(chǔ)上進(jìn)行再調(diào)度。為應(yīng)對(duì)風(fēng)電出力的不確定性,并保證調(diào)度計(jì)劃安全可靠,第2 階段問題要求在可消納風(fēng)電域內(nèi)最惡劣的風(fēng)電出力場(chǎng)景下不會(huì)產(chǎn)生棄風(fēng)與失負(fù)荷。第2 階段目標(biāo)函數(shù)如式(24)所示。

    1.3.4 第2 階段約束條件

    1)功率平衡與線路潮流約束

    第2 階段問題同樣包括節(jié)點(diǎn)功率平衡與線路潮流約束,與第1 階段不同的是引入了松弛變量,具體如式(25)至式(28)所示。

    本模型考慮了兩類不確定性。式(33)、式(34)采用盒式不確定集描述了所有子系統(tǒng)內(nèi)風(fēng)電出力的不確定性。為應(yīng)對(duì)受端子系統(tǒng)對(duì)聯(lián)絡(luò)線功率調(diào)節(jié)的需求,送端子系統(tǒng)將聯(lián)絡(luò)線功率視為具有不確定性的負(fù)荷,并采用式(35)、式(36)進(jìn)行描述。需要說明的是,這兩類不確定集的邊界均需要優(yōu)化,而不是事先設(shè)定的。

    2 模型求解

    2.1 運(yùn)行風(fēng)險(xiǎn)線性化

    根據(jù)式(37)的推導(dǎo)過程,由于任何形式的概率分布函數(shù)均滿足非負(fù)的要求,因此證明了運(yùn)行風(fēng)險(xiǎn)為凸函數(shù),同時(shí)說明了本文所采用的風(fēng)險(xiǎn)概念與線性化方法對(duì)任何風(fēng)電出力概率分布函數(shù)均適用。

    為減小分段線性化運(yùn)行風(fēng)險(xiǎn),首先將式(1)松弛為不等式(38)。

    2.2 分布式調(diào)度體系

    為減小求解模型規(guī)模、降低通信壓力,本文基于ADMM 建立了一種分布式調(diào)度體系,并對(duì)所提模型進(jìn)行了分布式求解。需要說明的是,本文所采用的ADMM 為一致性ADMM。文獻(xiàn)[31]證明了該算法適用于3 個(gè)及以上互聯(lián)調(diào)度主體的分布式優(yōu)化。該體系通過增廣拉格朗日函數(shù)對(duì)一致性約束式(19)、式(20)進(jìn)行松弛。松弛后可將原問題分解為多個(gè)本地子問題,并由本地調(diào)度系統(tǒng)進(jìn)行求解。各本地調(diào)度系統(tǒng)與相鄰系統(tǒng)交換耦合變量數(shù)值,更新本地子問題模型,并不斷迭代求解。當(dāng)相鄰子系統(tǒng)交換的耦合變量數(shù)值足夠接近時(shí),迭代完成并獲得最佳的調(diào)度方案。在第n次迭代過程中,送端子問題與受端子問題分別如式(41)與式(42)所示。

    為使ADMM 算法收斂,在每一次迭代結(jié)束后,各區(qū)域調(diào)度系統(tǒng)需要更新本地問題的拉格朗日乘子,更新方式如式(44)所示。

    ADMM 算法的收斂條件為所有耦合變量相關(guān)的原殘差與對(duì)偶?xì)埐罹∮诮o定限值,如式(45)所示。

    式中:σP為原殘差,表征相連子系統(tǒng)獲得的耦合變量數(shù)值是否足夠接近,εP為原殘差限值;σD為對(duì)偶?xì)埐?表征所得到的結(jié)果是否穩(wěn)定,εD為對(duì)偶?xì)埐钕拗怠?/p>

    對(duì)于傳統(tǒng)ADMM,懲罰因子ρ在求解過程中被固定為常數(shù)。然而,該算法的收斂性對(duì)ρ的數(shù)值較為敏感,不合適的懲罰因子可能導(dǎo)致算法收斂很慢甚至不收斂,但同時(shí)要事先確定一個(gè)合適的ρ的數(shù)值甚至是數(shù)量級(jí)卻很難。例如,對(duì)于不同的問題,文獻(xiàn)[32-33]中合適的懲罰因子分別為0.01 和150 000。因此,該算法的難點(diǎn)在于確定合適的懲罰因子。為解決這個(gè)問題,本文采用了一種動(dòng)態(tài)乘子更新策略[17]。該策略根據(jù)原殘差和對(duì)偶?xì)埐畹南鄬?duì)大小以及兩類殘差的收斂限值,在求解過程中動(dòng)態(tài)調(diào)整懲罰因子ρ的數(shù)值,詳細(xì)策略如式(46)所示。當(dāng)σP相較于σD過大時(shí),懲罰因子ρ將增大μin倍,以利于原殘差滿足收斂要求;而當(dāng)σD相較于σP過大時(shí),懲罰因子ρ將縮小為原來的1/μde。

    步驟2:送端與受端子系統(tǒng)并行求解本地子問題,獲得本地調(diào)度策略與聯(lián)絡(luò)線耦合變量。兩階段魯棒優(yōu)化子問題可采用列約束生成(column-andconstraint generation,C&CG)算法[34]進(jìn)行求解。需要說明的是,C&CG 算法采用主子問題迭代的形式將原問題轉(zhuǎn)化為帶有多場(chǎng)景的線性規(guī)劃問題?;诖?本文附錄A 對(duì)魯棒優(yōu)化問題的可分解性進(jìn)行了說明,即該問題可以采用分布式方法進(jìn)行求解。

    步驟3:送端與受端子系統(tǒng)交換最新的聯(lián)絡(luò)線耦合變量,并各自通過式(43)計(jì)算其平均值。

    步驟4:若收斂判據(jù)式(45)滿足條件,則迭代停止,并獲得最優(yōu)調(diào)度方案。如不滿足條件,則判斷總求解時(shí)間是否超過Tmax,若未超過,則各子系統(tǒng)按照式(46)與式(44)并行依次更新懲罰因子與拉格朗日乘子,設(shè)置n=n+1 并返回步驟2;若超過,則停止迭代,輸出當(dāng)前調(diào)度方案。需要說明的是,本文所提出的模型可以通過縮減風(fēng)電消納域的方式保證輸出調(diào)度方案的可行性。

    此外,在ADMM 算法與C&CG 算法嵌套迭代過程中,可以在每次C&CG 算法求解之后記錄下篩選出的最惡劣場(chǎng)景,并在下次迭代中將這些場(chǎng)景固定在主問題中。ADMM 與C&CG 算法嵌套迭代流程圖如附錄B 圖B1 所示。

    3 算例分析

    3.1 系統(tǒng)描述

    本文設(shè)計(jì)了3 區(qū)354 節(jié)點(diǎn)測(cè)試系統(tǒng)進(jìn)行算例分析。該系統(tǒng)由一個(gè)送端子系統(tǒng)和兩個(gè)受端子系統(tǒng)組成,每個(gè)區(qū)域都是修改的IEEE 118 節(jié)點(diǎn)測(cè)試系統(tǒng)。在送端子系統(tǒng)的節(jié)點(diǎn)5、54、77 和103 分別接入一個(gè)最大容量為800 MW 的風(fēng)電場(chǎng),在兩個(gè)受端子系統(tǒng)的節(jié)點(diǎn)77 處分別接入一個(gè)最大容量為800 MW 的風(fēng)電場(chǎng)。各風(fēng)電場(chǎng)出力假設(shè)服從正態(tài)分布,標(biāo)準(zhǔn)差為預(yù)測(cè)值的10%。

    各子系統(tǒng)24 h 風(fēng)電出力與負(fù)荷預(yù)測(cè)曲線如附錄B 圖B2 至圖B4 所示。該系統(tǒng)包含兩條聯(lián)絡(luò)線,聯(lián)絡(luò)線1 連接送端子系統(tǒng)的節(jié)點(diǎn)25 與受端子系統(tǒng)1的節(jié)點(diǎn)25,聯(lián)絡(luò)線2 連接送端子系統(tǒng)的節(jié)點(diǎn)65 與受端子系統(tǒng)2 的節(jié)點(diǎn)65。本文在MATLAB 2015a 平臺(tái)上采用商業(yè)優(yōu)化求解器Gurobi 9.0.0 對(duì)最優(yōu)問題進(jìn)行求解。

    3.2 不同動(dòng)態(tài)經(jīng)濟(jì)調(diào)度模式的比較

    本文所提的動(dòng)態(tài)經(jīng)濟(jì)調(diào)度模式(模式1)不但決定了各時(shí)段聯(lián)絡(luò)線的傳輸功率,而且確定了聯(lián)絡(luò)線各時(shí)段的功率調(diào)節(jié)區(qū)間。為說明該調(diào)度模式的優(yōu)勢(shì),表1 對(duì)比了其他兩種調(diào)度模式的調(diào)度結(jié)果:僅確定各時(shí)段聯(lián)絡(luò)線功率(模式2)、所有時(shí)段聯(lián)絡(luò)線功率均相同且不可調(diào)節(jié)(模式3)。從表中可以看出,由于本文所提出的調(diào)度模式可以充分發(fā)揮聯(lián)絡(luò)線功率靈活調(diào)節(jié)的特點(diǎn)并促進(jìn)區(qū)域間備用共享,系統(tǒng)發(fā)電與備用資源都能更加充分且合理地被使用,因此,相較于其他兩種調(diào)度模式可以獲得更低的運(yùn)行成本與運(yùn)行風(fēng)險(xiǎn),總成本分別減少了20 995 美元與81 740 美元。

    表1 不同動(dòng)態(tài)經(jīng)濟(jì)調(diào)度模式的結(jié)果比較Table 1 Result comparison of different dynamic economic dispatch modes

    圖3 給出了兩條聯(lián)絡(luò)線的基態(tài)功率與可調(diào)區(qū)間。當(dāng)受端子系統(tǒng)的實(shí)際風(fēng)電出力偏離預(yù)測(cè)值時(shí),送端子系統(tǒng)可以在確定的可調(diào)區(qū)間內(nèi)相應(yīng)調(diào)整聯(lián)絡(luò)線功率,以幫助其減輕風(fēng)電不確定性帶來的影響。采用這種調(diào)度模式,從全局角度考慮可以更多地利用更便宜的備用資源。相較于其他兩種調(diào)度模式,模式1 的總備用成本分別減少了14 795 美元和12 967 美元。

    圖3 聯(lián)絡(luò)線基態(tài)功率與可調(diào)區(qū)間Fig.3 Base-state power and adjustable intervals of tie lines

    同時(shí),本文所提調(diào)度模式得到的棄風(fēng)與失負(fù)荷懲罰均最小,這意味著該模式的可消納風(fēng)電域最大。附錄B 圖B5 給出了3 種調(diào)度模式所得到的送端子系統(tǒng)節(jié)點(diǎn)54 所連風(fēng)電場(chǎng)可消納風(fēng)電域的下界。從圖中可以看出,本文所提模型在時(shí)段9、10、15、16等得到的可消納風(fēng)電域下界明顯低于另外兩種調(diào)度模式,即在風(fēng)電出力因不確定性減小更多時(shí)依然能保證系統(tǒng)不失負(fù)荷,因而更加安全可靠。

    3.3 基于風(fēng)險(xiǎn)的魯棒模型與傳統(tǒng)魯棒模型的比較

    為解決傳統(tǒng)魯棒優(yōu)化模型難以預(yù)先設(shè)置合理不確定集的問題,本文采用了基于風(fēng)險(xiǎn)的魯棒優(yōu)化模型,在確定調(diào)度計(jì)劃的同時(shí),也對(duì)不確定集邊界進(jìn)行優(yōu)化,并得到相應(yīng)的可消納風(fēng)電域。為驗(yàn)證該模型的有效性,本節(jié)采用蒙特卡洛模擬法采樣10 000 個(gè)場(chǎng)景對(duì)4 種調(diào)度模型得到的結(jié)果進(jìn)行驗(yàn)證,其中一種為本文所提模型,另外3 種為傳統(tǒng)魯棒模型。表2比較了4 種模型的運(yùn)行成本與蒙特卡洛模擬產(chǎn)生的棄風(fēng)、失負(fù)荷懲罰。

    表2 含不同不確定集的魯棒模型結(jié)果比較Table 2 Result comparison of robust models with different uncertainty sets

    從表2 中可以看出,對(duì)于傳統(tǒng)魯棒模型,隨著所考慮不確定集的增大,調(diào)度計(jì)劃越來越可靠,即蒙特卡洛仿真得到的棄風(fēng)與失負(fù)荷懲罰越來越小,但同時(shí)發(fā)電與備用成本也更高,調(diào)度計(jì)劃也更保守。相比較而言,本文所采用的基于風(fēng)險(xiǎn)的魯棒優(yōu)化模型通過優(yōu)化不確定集邊界,實(shí)現(xiàn)了調(diào)度計(jì)劃經(jīng)濟(jì)性與魯棒性的平衡,總運(yùn)行成本也最低。同時(shí),采用本文模型也可以得到相應(yīng)的可消納風(fēng)電域。

    圖4 可消納風(fēng)電域Fig.4 Accommodable wind power region

    3.4 動(dòng)態(tài)乘子更新策略的有效性分析

    為了避免ADMM 參數(shù)選擇的難題,提高算法的收斂性,本文采用了動(dòng)態(tài)乘子更新策略。圖5 對(duì)比了該策略與傳統(tǒng)ADMM 在不同初始懲罰因子下的收斂過程。其中,圖5(a)中初始懲罰因子設(shè)置為1,圖5(b)中設(shè)置為0.1;最大殘差是指原殘差與對(duì)偶?xì)埐钪械妮^大值。原殘差與對(duì)偶?xì)埐畹氖諗肯拗稻O(shè)置為1 MW。

    在圖5(a)中,傳統(tǒng)ADMM 無法在100 次迭代中收斂,這意味著設(shè)置的初始懲罰因子的數(shù)量級(jí)不合適。相比之下,使用動(dòng)態(tài)乘子更新策略,該問題在56 次迭代后收斂到全局最優(yōu)解,表明了這一策略的有效性。在圖5(b)中,傳統(tǒng)ADMM 在49 次迭代后收斂,而動(dòng)態(tài)乘子更新策略將迭代次數(shù)減少為30,并行計(jì)算時(shí)間為1 838 s。由于動(dòng)態(tài)乘子更新策略避免了參數(shù)選擇的問題,并且具有更好的收斂性,因此在實(shí)踐中更加適用。

    圖5 不同乘子更新策略下的ADMM 收斂過程Fig.5 Convergence procedure of ADMM with different multiplier update strategies

    4 結(jié)語

    本文提出了一種基于風(fēng)險(xiǎn)的多區(qū)互聯(lián)電力系統(tǒng)分布式魯棒經(jīng)濟(jì)調(diào)度模型,合理配置聯(lián)絡(luò)線功率及區(qū)域內(nèi)發(fā)電與備用資源,以增強(qiáng)系統(tǒng)應(yīng)對(duì)風(fēng)電不確定性的能力,并采用分布式方式進(jìn)行求解。通過算例可得到如下結(jié)論。

    1)同時(shí)優(yōu)化基態(tài)聯(lián)絡(luò)線功率與可調(diào)區(qū)間,可以實(shí)現(xiàn)區(qū)域備用共享,更合理且充分地利用各區(qū)域的發(fā)電及備用資源,降低運(yùn)行成本與運(yùn)行風(fēng)險(xiǎn)。

    2)通過優(yōu)化不確定集邊界,可以實(shí)現(xiàn)調(diào)度策略經(jīng)濟(jì)性與魯棒性的平衡,同時(shí)評(píng)估系統(tǒng)棄風(fēng)與失負(fù)荷風(fēng)險(xiǎn),并獲得合理的可消納風(fēng)電域。在未來的研究工作中,將進(jìn)一步探索風(fēng)電場(chǎng)出力符合不同聯(lián)合分布規(guī)律對(duì)魯棒動(dòng)態(tài)經(jīng)濟(jì)調(diào)度結(jié)果的影響。

    3)采用分布式框架求解該問題可以減輕子系統(tǒng)間的通信負(fù)擔(dān),降低問題求解規(guī)模。采用動(dòng)態(tài)乘子更新策略可以避免ADMM 參數(shù)選擇的難題。同時(shí),算例分析表明,不同數(shù)量級(jí)初始懲罰因子條件下,該策略均大幅減小了迭代次數(shù),提升了算法的收斂性,增強(qiáng)了本文所提模型的有效性與實(shí)用性。

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

    猜你喜歡
    受端聯(lián)絡(luò)線魯棒
    基于短路比約束的受端電網(wǎng)新能源并網(wǎng)出力優(yōu)化方法研究
    特高壓交直流受端電網(wǎng)的穩(wěn)定性分析
    甬臺(tái)溫高速至沿海高速溫嶺聯(lián)絡(luò)線大溪樞紐設(shè)計(jì)
    城市軌道交通聯(lián)絡(luò)線計(jì)軸點(diǎn)的設(shè)置研究
    地鐵聯(lián)絡(luò)線無岔區(qū)段設(shè)計(jì)方案分析
    基于學(xué)習(xí)的魯棒自適應(yīng)評(píng)判控制研究進(jìn)展
    目標(biāo)魯棒識(shí)別的抗旋轉(zhuǎn)HDO 局部特征描述
    考慮短路容量和電壓穩(wěn)定約束的受端電網(wǎng)飽和負(fù)荷規(guī)模研究
    基于Cauchy魯棒函數(shù)的UKF改進(jìn)算法
    目標(biāo)軌跡更新的點(diǎn)到點(diǎn)魯棒迭代學(xué)習(xí)控制
    男人狂女人下面高潮的视频| 精品国产乱码久久久久久小说| 十分钟在线观看高清视频www | 精品少妇黑人巨大在线播放| 久久精品国产鲁丝片午夜精品| 亚洲精品一区蜜桃| 精品久久久噜噜| 97超碰精品成人国产| 波野结衣二区三区在线| 婷婷色麻豆天堂久久| 美女主播在线视频| 有码 亚洲区| 91久久精品国产一区二区三区| 啦啦啦中文免费视频观看日本| 国产黄片视频在线免费观看| 性色avwww在线观看| 熟女人妻精品中文字幕| 伊人久久国产一区二区| 中国美白少妇内射xxxbb| 亚洲欧美精品自产自拍| 91午夜精品亚洲一区二区三区| 国内揄拍国产精品人妻在线| 欧美少妇被猛烈插入视频| 毛片一级片免费看久久久久| 中文字幕精品免费在线观看视频 | 亚洲精品一区蜜桃| 久久av网站| 最近中文字幕2019免费版| 婷婷色麻豆天堂久久| 日日撸夜夜添| 亚洲美女视频黄频| 久久精品久久久久久噜噜老黄| 一边亲一边摸免费视频| xxx大片免费视频| 舔av片在线| 欧美日韩视频高清一区二区三区二| 最近中文字幕高清免费大全6| 久久青草综合色| 国产精品爽爽va在线观看网站| av国产精品久久久久影院| 丝袜脚勾引网站| 欧美成人午夜免费资源| 亚洲国产精品999| 欧美最新免费一区二区三区| 三级国产精品欧美在线观看| 国产精品偷伦视频观看了| 精品久久久久久电影网| 国产精品女同一区二区软件| 日本爱情动作片www.在线观看| 国语对白做爰xxxⅹ性视频网站| 欧美成人一区二区免费高清观看| 亚洲色图综合在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产精品女同一区二区软件| 久久 成人 亚洲| 国产美女午夜福利| av黄色大香蕉| 免费看光身美女| 国产永久视频网站| 一级毛片 在线播放| 91久久精品国产一区二区三区| 在线观看免费高清a一片| 日韩亚洲欧美综合| 亚洲精品中文字幕在线视频 | 亚洲自偷自拍三级| 午夜免费鲁丝| 免费观看的影片在线观看| 秋霞在线观看毛片| 亚洲精品国产av蜜桃| 国产日韩欧美在线精品| 亚洲精品aⅴ在线观看| 免费久久久久久久精品成人欧美视频 | 日日摸夜夜添夜夜爱| 女的被弄到高潮叫床怎么办| 久久鲁丝午夜福利片| 美女主播在线视频| 国产成人a∨麻豆精品| 欧美日韩亚洲高清精品| 最近中文字幕2019免费版| 一级毛片 在线播放| 2021少妇久久久久久久久久久| 亚洲国产毛片av蜜桃av| 国产黄色免费在线视频| 天堂8中文在线网| 蜜桃亚洲精品一区二区三区| 波野结衣二区三区在线| av又黄又爽大尺度在线免费看| 日韩亚洲欧美综合| 另类亚洲欧美激情| 看非洲黑人一级黄片| 国产精品久久久久久av不卡| 精品久久久噜噜| 22中文网久久字幕| av福利片在线观看| av在线观看视频网站免费| 综合色丁香网| av在线app专区| 在线天堂最新版资源| 国产精品人妻久久久影院| 在现免费观看毛片| 亚洲在久久综合| av免费观看日本| 91久久精品国产一区二区成人| 一级av片app| 五月开心婷婷网| 精品一区二区免费观看| 久久国内精品自在自线图片| 国产日韩欧美在线精品| 久久久亚洲精品成人影院| 久久影院123| 国产无遮挡羞羞视频在线观看| 人妻制服诱惑在线中文字幕| 99热网站在线观看| 成年人午夜在线观看视频| 国产av一区二区精品久久 | 男人和女人高潮做爰伦理| 久久这里有精品视频免费| 蜜臀久久99精品久久宅男| 观看美女的网站| 亚洲国产色片| 香蕉精品网在线| 国产亚洲一区二区精品| 亚洲国产成人一精品久久久| 国产黄色视频一区二区在线观看| 麻豆乱淫一区二区| 欧美成人午夜免费资源| av视频免费观看在线观看| 亚洲国产高清在线一区二区三| 五月开心婷婷网| 久久热精品热| 亚洲欧洲国产日韩| 中国三级夫妇交换| av又黄又爽大尺度在线免费看| 色视频在线一区二区三区| 激情五月婷婷亚洲| 丰满少妇做爰视频| 亚洲精品456在线播放app| 久久国产精品男人的天堂亚洲 | 在线观看av片永久免费下载| 久久人人爽av亚洲精品天堂 | 一级毛片 在线播放| 亚洲欧美中文字幕日韩二区| 免费大片黄手机在线观看| 我的女老师完整版在线观看| 人妻系列 视频| 欧美激情极品国产一区二区三区 | 亚洲av中文字字幕乱码综合| 国产精品福利在线免费观看| 国产成人a∨麻豆精品| 又黄又爽又刺激的免费视频.| 男女啪啪激烈高潮av片| 亚洲婷婷狠狠爱综合网| 日韩视频在线欧美| 少妇被粗大猛烈的视频| 一本一本综合久久| 国产精品国产av在线观看| av免费观看日本| 一本久久精品| 99re6热这里在线精品视频| 亚洲欧美成人综合另类久久久| 男男h啪啪无遮挡| 亚洲国产毛片av蜜桃av| 免费观看a级毛片全部| 亚洲精品国产av成人精品| 小蜜桃在线观看免费完整版高清| 欧美日韩精品成人综合77777| 亚洲精品国产av成人精品| 亚洲人成网站在线观看播放| av免费在线看不卡| 欧美丝袜亚洲另类| 国产高清不卡午夜福利| 亚洲精品国产色婷婷电影| 欧美极品一区二区三区四区| 国内少妇人妻偷人精品xxx网站| 校园人妻丝袜中文字幕| 男女无遮挡免费网站观看| 黄片无遮挡物在线观看| 一本久久精品| 十八禁网站网址无遮挡 | 菩萨蛮人人尽说江南好唐韦庄| 嫩草影院新地址| 久久99热6这里只有精品| 一个人看的www免费观看视频| 国产成人精品婷婷| 欧美丝袜亚洲另类| 国产免费一级a男人的天堂| 成人漫画全彩无遮挡| 久久久久久久久久久丰满| 国产爱豆传媒在线观看| 一区二区三区乱码不卡18| 男女无遮挡免费网站观看| 久久国产乱子免费精品| 乱系列少妇在线播放| 久久久久国产网址| av一本久久久久| 久久精品国产亚洲av天美| 我的老师免费观看完整版| 天堂8中文在线网| av在线老鸭窝| 欧美bdsm另类| 精品人妻熟女av久视频| 99视频精品全部免费 在线| 97在线视频观看| 国产精品蜜桃在线观看| 男女啪啪激烈高潮av片| 一区二区三区乱码不卡18| 97超碰精品成人国产| 亚洲伊人久久精品综合| 久久精品国产亚洲网站| 成年美女黄网站色视频大全免费 | 色哟哟·www| 精品99又大又爽又粗少妇毛片| av在线蜜桃| 国精品久久久久久国模美| 亚洲av成人精品一区久久| 国产真实伦视频高清在线观看| 国产精品久久久久成人av| 欧美激情国产日韩精品一区| 欧美少妇被猛烈插入视频| 高清毛片免费看| 在线观看av片永久免费下载| 99久久精品国产国产毛片| 91狼人影院| 亚洲精品色激情综合| 少妇猛男粗大的猛烈进出视频| 男女无遮挡免费网站观看| 在线观看一区二区三区| 韩国高清视频一区二区三区| 亚洲av欧美aⅴ国产| 毛片一级片免费看久久久久| 内射极品少妇av片p| 国产精品国产三级国产av玫瑰| 久久人人爽人人片av| 亚洲成人av在线免费| 2018国产大陆天天弄谢| 精品亚洲成国产av| 久久精品久久久久久噜噜老黄| 狂野欧美白嫩少妇大欣赏| 成人亚洲欧美一区二区av| 国产亚洲最大av| 伊人久久国产一区二区| 一区在线观看完整版| 亚洲欧美中文字幕日韩二区| 你懂的网址亚洲精品在线观看| 少妇裸体淫交视频免费看高清| 国产成人精品一,二区| 婷婷色综合大香蕉| 国产精品国产三级国产专区5o| 国产熟女欧美一区二区| 亚洲,一卡二卡三卡| 欧美一级a爱片免费观看看| 日本午夜av视频| 亚洲人与动物交配视频| 我的老师免费观看完整版| 久久97久久精品| a级一级毛片免费在线观看| 精品一区二区免费观看| 亚洲av综合色区一区| 两个人的视频大全免费| 国产成人免费无遮挡视频| 久久久欧美国产精品| 色吧在线观看| 国产乱人偷精品视频| av播播在线观看一区| 五月天丁香电影| 亚洲av二区三区四区| 深爱激情五月婷婷| 欧美高清成人免费视频www| 97超视频在线观看视频| 尤物成人国产欧美一区二区三区| 成人国产麻豆网| 简卡轻食公司| 人妻夜夜爽99麻豆av| 亚洲高清免费不卡视频| 国产精品一区二区三区四区免费观看| 欧美性感艳星| 18+在线观看网站| 亚洲自偷自拍三级| 最近最新中文字幕大全电影3| 欧美日韩在线观看h| 亚洲欧美中文字幕日韩二区| av在线蜜桃| 久久午夜福利片| 十八禁网站网址无遮挡 | 97精品久久久久久久久久精品| 久久久久久久久久人人人人人人| 久久97久久精品| 久久亚洲国产成人精品v| 亚洲国产欧美在线一区| 伊人久久国产一区二区| 久久99热这里只有精品18| 一级av片app| 久久久午夜欧美精品| 男人和女人高潮做爰伦理| 国产精品一区二区在线观看99| 午夜免费鲁丝| 99re6热这里在线精品视频| a级毛片免费高清观看在线播放| 一级爰片在线观看| 丰满人妻一区二区三区视频av| 欧美 日韩 精品 国产| 亚洲真实伦在线观看| 黄色视频在线播放观看不卡| 成人毛片a级毛片在线播放| 能在线免费看毛片的网站| 亚洲国产日韩一区二区| 免费观看性生交大片5| 97超碰精品成人国产| 亚洲内射少妇av| 在线观看人妻少妇| 成人国产av品久久久| 成人无遮挡网站| 十八禁网站网址无遮挡 | 天堂俺去俺来也www色官网| 久久热精品热| 国产黄色视频一区二区在线观看| 国产乱人视频| 51国产日韩欧美| 国产女主播在线喷水免费视频网站| 欧美成人一区二区免费高清观看| 色网站视频免费| 国产高清有码在线观看视频| 久久久久国产网址| 欧美zozozo另类| 91狼人影院| 久久ye,这里只有精品| 久久这里有精品视频免费| 中文精品一卡2卡3卡4更新| 老女人水多毛片| 日韩不卡一区二区三区视频在线| 亚洲精品国产av蜜桃| 亚洲精品日韩av片在线观看| 久热久热在线精品观看| 国产精品无大码| a级毛色黄片| 日韩av不卡免费在线播放| 麻豆成人午夜福利视频| 久久热精品热| 亚洲精品第二区| 国产乱人视频| 久久人人爽人人片av| 男女啪啪激烈高潮av片| 少妇熟女欧美另类| 亚洲精品,欧美精品| 国产精品久久久久成人av| 国产熟女欧美一区二区| 久久国内精品自在自线图片| av在线蜜桃| 在现免费观看毛片| 91久久精品电影网| a 毛片基地| 中文天堂在线官网| 国产精品一区www在线观看| 身体一侧抽搐| 水蜜桃什么品种好| 国产女主播在线喷水免费视频网站| 2021少妇久久久久久久久久久| 国产有黄有色有爽视频| 亚洲av免费高清在线观看| 你懂的网址亚洲精品在线观看| 亚洲国产最新在线播放| 亚洲电影在线观看av| 中文字幕精品免费在线观看视频 | 久久精品国产亚洲网站| 亚洲欧美精品专区久久| 久热久热在线精品观看| 国内揄拍国产精品人妻在线| 尤物成人国产欧美一区二区三区| 一本—道久久a久久精品蜜桃钙片| 精品人妻一区二区三区麻豆| 爱豆传媒免费全集在线观看| 日韩电影二区| 国产视频首页在线观看| 欧美人与善性xxx| 色视频在线一区二区三区| 99国产精品免费福利视频| 校园人妻丝袜中文字幕| 精品酒店卫生间| 国精品久久久久久国模美| 黄色配什么色好看| 视频中文字幕在线观看| 五月开心婷婷网| 六月丁香七月| 一级毛片黄色毛片免费观看视频| 亚洲精品日本国产第一区| 两个人的视频大全免费| 亚洲欧美日韩另类电影网站 | 国产精品嫩草影院av在线观看| 午夜福利在线观看免费完整高清在| 99久国产av精品国产电影| 成年美女黄网站色视频大全免费 | 欧美xxxx黑人xx丫x性爽| 色网站视频免费| 国产在视频线精品| 人妻夜夜爽99麻豆av| 国产人妻一区二区三区在| av网站免费在线观看视频| 久久精品夜色国产| 精品人妻视频免费看| 性色av一级| av播播在线观看一区| 美女视频免费永久观看网站| 99热国产这里只有精品6| 国产91av在线免费观看| 99久久中文字幕三级久久日本| 中文在线观看免费www的网站| 亚洲精品乱久久久久久| 人妻 亚洲 视频| 成人特级av手机在线观看| 欧美日本视频| 久久青草综合色| 五月开心婷婷网| av国产久精品久网站免费入址| 你懂的网址亚洲精品在线观看| 日韩亚洲欧美综合| 国产黄色免费在线视频| 如何舔出高潮| 高清日韩中文字幕在线| 成人一区二区视频在线观看| 丰满迷人的少妇在线观看| 国产成人freesex在线| 国产成人免费无遮挡视频| 久久国产精品大桥未久av | 在线观看人妻少妇| 麻豆国产97在线/欧美| 久久久欧美国产精品| 青青草视频在线视频观看| 免费大片黄手机在线观看| 亚洲婷婷狠狠爱综合网| 国产精品一区二区性色av| 少妇 在线观看| 国产精品无大码| 国产爽快片一区二区三区| 七月丁香在线播放| a级一级毛片免费在线观看| 国产男女内射视频| 人妻一区二区av| 一区二区三区四区激情视频| 欧美区成人在线视频| 五月天丁香电影| 亚洲一区二区三区欧美精品| 特大巨黑吊av在线直播| 日韩大片免费观看网站| 国产午夜精品一二区理论片| 免费观看性生交大片5| 波野结衣二区三区在线| 精品一区二区三区视频在线| 亚洲真实伦在线观看| 亚洲国产成人一精品久久久| 各种免费的搞黄视频| 亚洲丝袜综合中文字幕| 一级二级三级毛片免费看| 成人高潮视频无遮挡免费网站| 王馨瑶露胸无遮挡在线观看| 一级av片app| 亚洲天堂av无毛| 高清在线视频一区二区三区| 免费久久久久久久精品成人欧美视频 | 在线看a的网站| 一级黄片播放器| av免费观看日本| 精品人妻一区二区三区麻豆| 多毛熟女@视频| 国产伦在线观看视频一区| 91精品伊人久久大香线蕉| 高清午夜精品一区二区三区| 亚洲欧美一区二区三区国产| 婷婷色综合大香蕉| 少妇精品久久久久久久| av一本久久久久| 看免费成人av毛片| 国内精品宾馆在线| 欧美亚洲 丝袜 人妻 在线| 中文字幕免费在线视频6| 我的女老师完整版在线观看| 久久毛片免费看一区二区三区| 亚洲精品久久午夜乱码| 欧美区成人在线视频| 国产美女午夜福利| 欧美日韩精品成人综合77777| 中文在线观看免费www的网站| 色哟哟·www| 国产探花极品一区二区| 亚洲国产最新在线播放| 亚洲一区二区三区欧美精品| 内地一区二区视频在线| 久久久国产一区二区| 国产伦理片在线播放av一区| 好男人视频免费观看在线| 欧美精品一区二区大全| 成人二区视频| 日日撸夜夜添| 97超碰精品成人国产| 美女福利国产在线 | 嫩草影院入口| 中文天堂在线官网| 亚洲国产成人一精品久久久| 777米奇影视久久| 亚洲国产精品999| 欧美成人一区二区免费高清观看| 日韩av不卡免费在线播放| 纯流量卡能插随身wifi吗| 大码成人一级视频| 亚洲性久久影院| 最后的刺客免费高清国语| 91午夜精品亚洲一区二区三区| 中文字幕亚洲精品专区| 欧美精品国产亚洲| 不卡视频在线观看欧美| 99久久精品国产国产毛片| 欧美亚洲 丝袜 人妻 在线| 久久精品国产a三级三级三级| 久久久国产一区二区| 国产探花极品一区二区| 国产精品熟女久久久久浪| 日韩av在线免费看完整版不卡| 亚洲人成网站在线播| 一级片'在线观看视频| 伊人久久精品亚洲午夜| 日韩亚洲欧美综合| av国产免费在线观看| 日本欧美视频一区| 久久青草综合色| 三级国产精品欧美在线观看| 免费看日本二区| 黑人高潮一二区| 高清视频免费观看一区二区| 99久久精品国产国产毛片| av播播在线观看一区| 最近最新中文字幕大全电影3| 成人18禁高潮啪啪吃奶动态图 | 老女人水多毛片| 亚洲国产精品一区三区| 夜夜看夜夜爽夜夜摸| 国产视频内射| 最近最新中文字幕免费大全7| 人人妻人人看人人澡| 王馨瑶露胸无遮挡在线观看| 午夜日本视频在线| 国产久久久一区二区三区| a级毛片免费高清观看在线播放| 国产黄片美女视频| 亚洲精品成人av观看孕妇| 日本vs欧美在线观看视频 | 国产精品偷伦视频观看了| 国产亚洲av片在线观看秒播厂| 美女xxoo啪啪120秒动态图| 夜夜看夜夜爽夜夜摸| 干丝袜人妻中文字幕| 五月玫瑰六月丁香| 免费大片黄手机在线观看| 丝袜脚勾引网站| 精品国产三级普通话版| 小蜜桃在线观看免费完整版高清| videossex国产| 午夜日本视频在线| 精品亚洲成a人片在线观看 | 人体艺术视频欧美日本| 久久精品久久久久久久性| 国产一区亚洲一区在线观看| 伦精品一区二区三区| 男人和女人高潮做爰伦理| 观看av在线不卡| av专区在线播放| 亚洲高清免费不卡视频| 免费高清在线观看视频在线观看| 国产色爽女视频免费观看| 久久99热这里只频精品6学生| 边亲边吃奶的免费视频| 欧美极品一区二区三区四区| 免费av中文字幕在线| 久久国产精品大桥未久av | 在线精品无人区一区二区三 | 亚洲美女黄色视频免费看| 亚洲第一av免费看| 亚洲aⅴ乱码一区二区在线播放| 久久精品人妻少妇| 在线观看人妻少妇| 啦啦啦中文免费视频观看日本| 国产爱豆传媒在线观看| 国产乱人偷精品视频| 青青草视频在线视频观看| 国产免费一级a男人的天堂| 欧美成人午夜免费资源| 亚洲av电影在线观看一区二区三区| 欧美三级亚洲精品| 中文字幕亚洲精品专区| 精品一区在线观看国产| 国产精品久久久久久久电影| 免费看日本二区| 亚洲欧洲日产国产| 美女脱内裤让男人舔精品视频| 国产精品伦人一区二区| 欧美国产精品一级二级三级 | 日韩大片免费观看网站| 日韩,欧美,国产一区二区三区| 国产黄色视频一区二区在线观看| 高清日韩中文字幕在线| 亚洲国产日韩一区二区| 在线免费观看不下载黄p国产| 高清日韩中文字幕在线| 日本黄色日本黄色录像| 久久久久国产精品人妻一区二区| 日本vs欧美在线观看视频 | 黄片wwwwww| 蜜桃在线观看..| 色婷婷久久久亚洲欧美| 九九久久精品国产亚洲av麻豆| 亚洲,欧美,日韩| 精品人妻熟女av久视频| 十分钟在线观看高清视频www | 日本免费在线观看一区| 狠狠精品人妻久久久久久综合| 亚洲内射少妇av| 成人亚洲精品一区在线观看 |