• <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í)控制
    久久精品影院6| 丰满的人妻完整版| 一级毛片女人18水好多| 国产伦人伦偷精品视频| ponron亚洲| 婷婷亚洲欧美| 免费看a级黄色片| 亚洲avbb在线观看| 五月玫瑰六月丁香| 一区二区三区高清视频在线| 国产成人啪精品午夜网站| 国产亚洲精品久久久久5区| 亚洲激情在线av| 久久久久免费精品人妻一区二区| 亚洲国产精品久久男人天堂| 精品福利观看| 夜夜夜夜夜久久久久| 久久精品影院6| 日韩国内少妇激情av| 欧美黑人精品巨大| 欧美大码av| 久99久视频精品免费| 成年人黄色毛片网站| 午夜日韩欧美国产| 蜜桃久久精品国产亚洲av| 99热只有精品国产| 亚洲av片天天在线观看| 一进一出抽搐动态| 婷婷亚洲欧美| 最新在线观看一区二区三区| 欧美在线黄色| 欧美乱码精品一区二区三区| 国内毛片毛片毛片毛片毛片| 久久国产精品影院| 成人18禁高潮啪啪吃奶动态图| 一进一出抽搐gif免费好疼| 精品久久久久久久久久免费视频| 女人高潮潮喷娇喘18禁视频| 亚洲中文av在线| 亚洲一区二区三区色噜噜| 国产99久久九九免费精品| 精品第一国产精品| 欧美日韩亚洲国产一区二区在线观看| 日本免费a在线| 女生性感内裤真人,穿戴方法视频| 嫩草影院精品99| 日本免费a在线| 一个人观看的视频www高清免费观看 | 真人做人爱边吃奶动态| 欧美一级毛片孕妇| 熟女少妇亚洲综合色aaa.| 淫妇啪啪啪对白视频| 久久热在线av| 国产激情偷乱视频一区二区| 国产精品亚洲美女久久久| 亚洲av五月六月丁香网| 国产精品1区2区在线观看.| 妹子高潮喷水视频| 国产av在哪里看| 午夜免费观看网址| 欧美黑人巨大hd| 成熟少妇高潮喷水视频| 亚洲成人免费电影在线观看| 人人妻,人人澡人人爽秒播| 五月玫瑰六月丁香| aaaaa片日本免费| 好看av亚洲va欧美ⅴa在| 国产男靠女视频免费网站| 国产精品爽爽va在线观看网站| 可以在线观看的亚洲视频| 黄色a级毛片大全视频| 亚洲精品在线美女| 久久久久久久午夜电影| 国模一区二区三区四区视频 | 亚洲专区字幕在线| 久久人人精品亚洲av| 视频区欧美日本亚洲| 很黄的视频免费| 亚洲中文日韩欧美视频| 欧美在线黄色| 久久性视频一级片| 天堂av国产一区二区熟女人妻 | 精华霜和精华液先用哪个| 亚洲熟妇中文字幕五十中出| 九色成人免费人妻av| 99国产精品一区二区蜜桃av| 久久午夜亚洲精品久久| 久久久久久久久中文| 日韩免费av在线播放| 国产精品九九99| 99精品久久久久人妻精品| 一级毛片精品| 此物有八面人人有两片| 亚洲一码二码三码区别大吗| 国产av一区在线观看免费| 夜夜看夜夜爽夜夜摸| 亚洲av美国av| 国内精品久久久久久久电影| 老司机在亚洲福利影院| 国产免费av片在线观看野外av| 1024视频免费在线观看| 国产探花在线观看一区二区| 97碰自拍视频| 少妇裸体淫交视频免费看高清 | 嫩草影视91久久| 亚洲性夜色夜夜综合| av片东京热男人的天堂| 国语自产精品视频在线第100页| 91麻豆精品激情在线观看国产| 久久性视频一级片| 色综合婷婷激情| 精品人妻1区二区| 欧美日韩亚洲综合一区二区三区_| 婷婷丁香在线五月| 国产精品1区2区在线观看.| 久久中文看片网| 香蕉国产在线看| 国产1区2区3区精品| 波多野结衣高清无吗| 亚洲性夜色夜夜综合| 老司机午夜福利在线观看视频| 国产精品乱码一区二三区的特点| 国产91精品成人一区二区三区| 在线观看免费日韩欧美大片| 午夜日韩欧美国产| 国产主播在线观看一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 久久精品亚洲精品国产色婷小说| 久久精品国产清高在天天线| 久久久精品国产亚洲av高清涩受| 黄色成人免费大全| 国产蜜桃级精品一区二区三区| 人成视频在线观看免费观看| 欧美一级a爱片免费观看看 | 日本一本二区三区精品| 色在线成人网| 色综合站精品国产| 最近在线观看免费完整版| 蜜桃久久精品国产亚洲av| 99久久99久久久精品蜜桃| 欧美一级毛片孕妇| 女人被狂操c到高潮| 日本 欧美在线| av超薄肉色丝袜交足视频| 日韩有码中文字幕| 午夜福利成人在线免费观看| 两个人看的免费小视频| 丁香六月欧美| 国产乱人伦免费视频| 色哟哟哟哟哟哟| 久久久国产欧美日韩av| 亚洲男人天堂网一区| 亚洲aⅴ乱码一区二区在线播放 | 久久香蕉精品热| 男女午夜视频在线观看| 男人的好看免费观看在线视频 | 国产成人精品无人区| 色尼玛亚洲综合影院| 三级国产精品欧美在线观看 | 亚洲欧美激情综合另类| 变态另类成人亚洲欧美熟女| 久久久久久大精品| 亚洲第一欧美日韩一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | 黄色视频,在线免费观看| 欧美一级毛片孕妇| 亚洲专区中文字幕在线| 国产精品 国内视频| 91老司机精品| 97人妻精品一区二区三区麻豆| 亚洲黑人精品在线| 9191精品国产免费久久| 深夜精品福利| 在线观看免费午夜福利视频| netflix在线观看网站| 国产一区二区在线观看日韩 | 99久久久亚洲精品蜜臀av| 黄色视频不卡| 久久久久国产一级毛片高清牌| 老司机在亚洲福利影院| videosex国产| 一卡2卡三卡四卡精品乱码亚洲| 动漫黄色视频在线观看| 欧美性猛交黑人性爽| 精品欧美一区二区三区在线| 亚洲人成网站在线播放欧美日韩| 国产亚洲欧美在线一区二区| 日韩欧美在线二视频| 日本黄大片高清| 国产精品精品国产色婷婷| 精品午夜福利视频在线观看一区| 亚洲专区中文字幕在线| 一本久久中文字幕| 一区二区三区国产精品乱码| 国产精品久久久久久精品电影| 真人一进一出gif抽搐免费| 日韩欧美在线乱码| 国产又色又爽无遮挡免费看| 99热6这里只有精品| 亚洲午夜精品一区,二区,三区| 1024香蕉在线观看| 成年免费大片在线观看| 午夜视频精品福利| 久久精品91无色码中文字幕| 一本大道久久a久久精品| 久久精品国产清高在天天线| www国产在线视频色| 动漫黄色视频在线观看| 免费人成视频x8x8入口观看| 亚洲av成人一区二区三| 亚洲av美国av| 午夜成年电影在线免费观看| 精品久久久久久,| 毛片女人毛片| 老司机深夜福利视频在线观看| 88av欧美| 熟女少妇亚洲综合色aaa.| 国产爱豆传媒在线观看 | 国产精品免费视频内射| 亚洲九九香蕉| 亚洲在线自拍视频| 少妇被粗大的猛进出69影院| 国产精品亚洲av一区麻豆| 亚洲国产高清在线一区二区三| 久久人妻福利社区极品人妻图片| 亚洲av成人精品一区久久| 村上凉子中文字幕在线| 成人手机av| 亚洲欧美日韩东京热| 中文字幕人妻丝袜一区二区| 特大巨黑吊av在线直播| 欧美+亚洲+日韩+国产| 久久久久久久精品吃奶| 脱女人内裤的视频| 草草在线视频免费看| 黄色片一级片一级黄色片| netflix在线观看网站| 又粗又爽又猛毛片免费看| 午夜a级毛片| 久久中文看片网| 中文字幕熟女人妻在线| 两个人视频免费观看高清| 精品久久久久久久末码| 很黄的视频免费| 国产v大片淫在线免费观看| 成年女人毛片免费观看观看9| 亚洲男人的天堂狠狠| 精品一区二区三区av网在线观看| 亚洲国产欧美网| 欧美三级亚洲精品| 国产精品久久久av美女十八| 深夜精品福利| 国产成人精品无人区| 99在线视频只有这里精品首页| 午夜激情福利司机影院| 女人爽到高潮嗷嗷叫在线视频| 亚洲成人久久性| 亚洲中文av在线| 又爽又黄无遮挡网站| 成人国产综合亚洲| 久久久久久久久久黄片| 久久久久久久精品吃奶| 亚洲天堂国产精品一区在线| 久久久久久九九精品二区国产 | 全区人妻精品视频| 国产蜜桃级精品一区二区三区| 亚洲欧美日韩无卡精品| 高清毛片免费观看视频网站| 亚洲av成人精品一区久久| 不卡一级毛片| 中文字幕熟女人妻在线| 蜜桃久久精品国产亚洲av| 国产精品亚洲一级av第二区| 露出奶头的视频| 在线播放国产精品三级| 免费在线观看影片大全网站| 亚洲国产欧洲综合997久久,| 热99re8久久精品国产| 中文资源天堂在线| 最好的美女福利视频网| 51午夜福利影视在线观看| 美女扒开内裤让男人捅视频| 国产精品1区2区在线观看.| 亚洲精品色激情综合| 老司机午夜十八禁免费视频| 婷婷亚洲欧美| 亚洲精品色激情综合| 俺也久久电影网| 精品第一国产精品| 天天躁夜夜躁狠狠躁躁| 亚洲人成77777在线视频| 精品熟女少妇八av免费久了| 久久伊人香网站| 一本综合久久免费| 此物有八面人人有两片| 亚洲精品美女久久久久99蜜臀| 日韩欧美免费精品| 可以在线观看的亚洲视频| 日日夜夜操网爽| 亚洲一区二区三区色噜噜| 黄色女人牲交| 老鸭窝网址在线观看| 在线观看日韩欧美| 国产激情久久老熟女| 欧美日本视频| 好男人电影高清在线观看| 香蕉av资源在线| 久久香蕉激情| 欧美激情久久久久久爽电影| 久久久久久国产a免费观看| 日本在线视频免费播放| 一进一出抽搐gif免费好疼| 国产精品爽爽va在线观看网站| 美女扒开内裤让男人捅视频| 国内精品一区二区在线观看| 午夜免费成人在线视频| 亚洲18禁久久av| 久久九九热精品免费| 国产成人欧美在线观看| 国产熟女xx| 日韩欧美国产一区二区入口| 国产精品久久久久久亚洲av鲁大| 最好的美女福利视频网| 国产私拍福利视频在线观看| 欧美日韩一级在线毛片| 中文在线观看免费www的网站 | 999久久久国产精品视频| 在线观看www视频免费| 悠悠久久av| 亚洲专区中文字幕在线| av欧美777| 亚洲欧美一区二区三区黑人| 亚洲国产中文字幕在线视频| 国产视频一区二区在线看| av天堂在线播放| 久久中文字幕人妻熟女| 99久久国产精品久久久| 18美女黄网站色大片免费观看| 精品久久久久久久久久久久久| 欧美又色又爽又黄视频| 人成视频在线观看免费观看| 窝窝影院91人妻| 美女免费视频网站| 免费搜索国产男女视频| 在线播放国产精品三级| 人妻夜夜爽99麻豆av| 久久久精品欧美日韩精品| 国产片内射在线| 欧美+亚洲+日韩+国产| 大型av网站在线播放| 麻豆一二三区av精品| 国产野战对白在线观看| 精品熟女少妇八av免费久了| 国产亚洲精品一区二区www| 他把我摸到了高潮在线观看| 狂野欧美激情性xxxx| 妹子高潮喷水视频| 夜夜爽天天搞| 久久精品国产综合久久久| 国产免费av片在线观看野外av| 久久久精品国产亚洲av高清涩受| 一边摸一边做爽爽视频免费| 露出奶头的视频| 观看免费一级毛片| 亚洲国产精品sss在线观看| 久9热在线精品视频| 久久久久国产精品人妻aⅴ院| 黄片小视频在线播放| 国产亚洲欧美98| 999精品在线视频| 亚洲黑人精品在线| 一级毛片高清免费大全| 国产免费av片在线观看野外av| 精品无人区乱码1区二区| 精品国产亚洲在线| 亚洲一码二码三码区别大吗| netflix在线观看网站| 亚洲男人天堂网一区| 免费av毛片视频| 国产精品久久视频播放| 99国产精品一区二区三区| 欧美一区二区国产精品久久精品 | 麻豆一二三区av精品| 国产精品一及| 欧美中文综合在线视频| 久久人人精品亚洲av| 制服诱惑二区| av福利片在线观看| 久久精品影院6| www国产在线视频色| 我的老师免费观看完整版| 日本免费a在线| 丁香六月欧美| 一a级毛片在线观看| 欧美人与性动交α欧美精品济南到| 在线十欧美十亚洲十日本专区| 久久亚洲精品不卡| 夜夜夜夜夜久久久久| 精品一区二区三区视频在线观看免费| 麻豆一二三区av精品| 欧美午夜高清在线| 久久这里只有精品19| 国产精品1区2区在线观看.| 中文字幕熟女人妻在线| 无人区码免费观看不卡| 丁香六月欧美| 欧美在线黄色| 国产熟女午夜一区二区三区| 久久 成人 亚洲| 久久久久久九九精品二区国产 | 欧美一级a爱片免费观看看 | 又大又爽又粗| 757午夜福利合集在线观看| 欧美日韩中文字幕国产精品一区二区三区| 日本撒尿小便嘘嘘汇集6| 一区二区三区高清视频在线| 99久久无色码亚洲精品果冻| 国产免费男女视频| 香蕉丝袜av| 久久 成人 亚洲| 亚洲av成人av| 50天的宝宝边吃奶边哭怎么回事| 午夜免费观看网址| 精品午夜福利视频在线观看一区| 欧美+亚洲+日韩+国产| 在线观看日韩欧美| 悠悠久久av| 高清毛片免费观看视频网站| 一级黄色大片毛片| 亚洲五月天丁香| 亚洲在线自拍视频| 午夜免费观看网址| 黄色丝袜av网址大全| 日韩精品青青久久久久久| 亚洲av五月六月丁香网| 老司机午夜十八禁免费视频| 亚洲最大成人中文| 亚洲九九香蕉| 脱女人内裤的视频| av欧美777| 久久久久国产一级毛片高清牌| 2021天堂中文幕一二区在线观| 欧美日本亚洲视频在线播放| 久久人妻福利社区极品人妻图片| 国产精品电影一区二区三区| 男女下面进入的视频免费午夜| 午夜精品久久久久久毛片777| 一边摸一边做爽爽视频免费| 一个人观看的视频www高清免费观看 | 精品人妻1区二区| 欧美色欧美亚洲另类二区| 伊人久久大香线蕉亚洲五| 美女 人体艺术 gogo| 国产99白浆流出| 亚洲国产欧美一区二区综合| 国产精品乱码一区二三区的特点| 国产成人一区二区三区免费视频网站| 熟女电影av网| 国产精品久久视频播放| 国产乱人伦免费视频| 亚洲国产精品成人综合色| 欧美高清成人免费视频www| 欧美日韩中文字幕国产精品一区二区三区| 午夜精品一区二区三区免费看| 女人高潮潮喷娇喘18禁视频| 给我免费播放毛片高清在线观看| 欧美激情久久久久久爽电影| 亚洲人成77777在线视频| 成熟少妇高潮喷水视频| 在线播放国产精品三级| 午夜福利18| 国产av又大| 日韩国内少妇激情av| 亚洲精品国产精品久久久不卡| 精品久久久久久久久久久久久| 人人妻人人澡欧美一区二区| 黄片大片在线免费观看| 亚洲午夜理论影院| 国产精品久久久人人做人人爽| 精品日产1卡2卡| 少妇人妻一区二区三区视频| 91国产中文字幕| 久久久久免费精品人妻一区二区| 美女高潮喷水抽搐中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 免费看a级黄色片| 无人区码免费观看不卡| 国产三级中文精品| 亚洲成av人片免费观看| 一级黄色大片毛片| 国产精品 国内视频| 亚洲国产精品合色在线| 黄色视频,在线免费观看| 级片在线观看| 国产精品久久久av美女十八| 十八禁人妻一区二区| 在线观看舔阴道视频| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲av电影不卡..在线观看| 老熟妇仑乱视频hdxx| 成人手机av| 999久久久精品免费观看国产| 免费人成视频x8x8入口观看| 国产亚洲精品av在线| 变态另类丝袜制服| 高潮久久久久久久久久久不卡| 久久中文看片网| 熟女少妇亚洲综合色aaa.| 亚洲国产精品999在线| 国产欧美日韩一区二区精品| 亚洲一区中文字幕在线| 欧美日韩亚洲综合一区二区三区_| 欧美一级a爱片免费观看看 | 国产成人一区二区三区免费视频网站| 12—13女人毛片做爰片一| 男插女下体视频免费在线播放| 欧美成人午夜精品| 俺也久久电影网| 一区二区三区激情视频| 日韩大码丰满熟妇| 成年版毛片免费区| 一级片免费观看大全| 在线a可以看的网站| 黄频高清免费视频| 激情在线观看视频在线高清| 青草久久国产| 小说图片视频综合网站| 最近最新免费中文字幕在线| 国产亚洲av高清不卡| 久久性视频一级片| 亚洲成人免费电影在线观看| 久久久久久久精品吃奶| 欧美日韩国产亚洲二区| 国产一区二区三区在线臀色熟女| 日韩欧美在线乱码| av在线天堂中文字幕| 一级毛片精品| 国产一区二区在线av高清观看| 怎么达到女性高潮| 亚洲av成人av| 精品久久久久久久人妻蜜臀av| 19禁男女啪啪无遮挡网站| 亚洲在线自拍视频| 亚洲国产欧洲综合997久久,| 欧美一级a爱片免费观看看 | 欧洲精品卡2卡3卡4卡5卡区| 日韩精品免费视频一区二区三区| 两人在一起打扑克的视频| 女人高潮潮喷娇喘18禁视频| 国产成人影院久久av| 久热爱精品视频在线9| 精品久久久久久久末码| 亚洲色图av天堂| 免费人成视频x8x8入口观看| 色老头精品视频在线观看| 亚洲性夜色夜夜综合| 免费观看精品视频网站| 床上黄色一级片| 日韩大尺度精品在线看网址| 久久久国产精品麻豆| 亚洲国产精品久久男人天堂| 久久99热这里只有精品18| 亚洲国产精品久久男人天堂| 久久婷婷人人爽人人干人人爱| 伦理电影免费视频| 国产亚洲精品第一综合不卡| 女警被强在线播放| 一区二区三区激情视频| 久久性视频一级片| 久久精品91蜜桃| 人成视频在线观看免费观看| 欧美精品亚洲一区二区| 母亲3免费完整高清在线观看| 欧美高清成人免费视频www| 天天添夜夜摸| 两性午夜刺激爽爽歪歪视频在线观看 | 男女视频在线观看网站免费 | 欧美成人性av电影在线观看| 欧美精品亚洲一区二区| 又大又爽又粗| 欧美成人午夜精品| 1024手机看黄色片| 精品电影一区二区在线| 在线观看免费午夜福利视频| 伊人久久大香线蕉亚洲五| 成年免费大片在线观看| 亚洲国产精品999在线| 国产激情久久老熟女| 国产蜜桃级精品一区二区三区| 精品久久久久久久久久免费视频| 国产一区二区在线av高清观看| 一级毛片女人18水好多| 亚洲无线在线观看| 老司机深夜福利视频在线观看| 99国产精品一区二区蜜桃av| 欧美绝顶高潮抽搐喷水| 免费观看人在逋| 白带黄色成豆腐渣| 午夜精品在线福利| 久久久久久九九精品二区国产 | 精品电影一区二区在线| 18禁黄网站禁片午夜丰满| 精品人妻1区二区| 亚洲一区高清亚洲精品| 中文字幕人成人乱码亚洲影| 欧美丝袜亚洲另类 | 精品第一国产精品| 国产精品一区二区三区四区久久| 久久久精品大字幕| 免费高清视频大片| 麻豆一二三区av精品| 欧美日韩乱码在线| 亚洲国产日韩欧美精品在线观看 | av免费在线观看网站| 一个人观看的视频www高清免费观看 |