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

    線性微分遞減自適應(yīng)粒子群算法在水電站優(yōu)化調(diào)度中的應(yīng)用

    2018-01-08 02:26:31萬(wàn)忠海葉生進(jìn)
    水力發(fā)電 2017年9期
    關(guān)鍵詞:微分慣性水電站

    萬(wàn)忠海,葉生進(jìn),鄭 姣

    (四川省都江堰東風(fēng)渠管理處,四川成都610072)

    線性微分遞減自適應(yīng)粒子群算法在水電站優(yōu)化調(diào)度中的應(yīng)用

    萬(wàn)忠海,葉生進(jìn),鄭 姣

    (四川省都江堰東風(fēng)渠管理處,四川成都610072)

    針對(duì)求解水電站優(yōu)化調(diào)度粒子群算法的改進(jìn),分析了粒子群算法在求解水電站優(yōu)化調(diào)度問(wèn)題時(shí)對(duì)關(guān)鍵參數(shù)慣性權(quán)重調(diào)整的需要,提出了線性微分遞減的自適應(yīng)粒子群算法。通過(guò)前期減小緩慢的慣性權(quán)重,增加算法的探索能力跳出局部最優(yōu)解;通過(guò)后期減小較快的慣性權(quán)重,提升算法的開(kāi)發(fā)能力加快算法收斂。以葛洲壩水電站優(yōu)化調(diào)度為例,對(duì)比了改進(jìn)算法和傳統(tǒng)算法。優(yōu)化調(diào)度實(shí)例表明:線性微分遞減自適應(yīng)策略增強(qiáng)了算法的尋優(yōu)能力和穩(wěn)定性。改進(jìn)算法能夠有效改善由于水電站優(yōu)化調(diào)度目標(biāo)函數(shù)非凸性帶來(lái)的粒子群求解易早熟問(wèn)題,為水電站優(yōu)化調(diào)度粒子群算法慣性權(quán)重的改進(jìn)提供了新思路。

    水電站優(yōu)化調(diào)度;粒子群算法;慣性權(quán)重;線性微分遞減

    0 引 言

    水電站優(yōu)化調(diào)度可以在不增加額外投資的情況下顯著增加收益,其數(shù)學(xué)模型求解的傳統(tǒng)方法有:動(dòng)態(tài)規(guī)劃(DP)、增量動(dòng)態(tài)規(guī)劃算法(IDP)和逐步優(yōu)化算法(POA)等。傳統(tǒng)理論較為成熟與完善,卻或多或少的存在“維數(shù)災(zāi)”問(wèn)題[1]。因此傳統(tǒng)算法僅適合單庫(kù)優(yōu)化調(diào)度的求解,對(duì)大規(guī)模水電站水庫(kù)群優(yōu)化調(diào)度的求解存在局限性。

    粒子群算法(PSO)是通過(guò)模擬鳥(niǎo)群覓食行為而發(fā)展起來(lái)的基于群體協(xié)作的隨機(jī)搜索算法。作為一種新興的群體智能進(jìn)化算法,該算法以流程簡(jiǎn)單易實(shí)現(xiàn)和參數(shù)簡(jiǎn)潔著稱[2]。目前已廣泛應(yīng)用于函數(shù)優(yōu)化、數(shù)據(jù)挖掘和模糊系統(tǒng)控制等領(lǐng)域。從可獲文獻(xiàn)知,PSO應(yīng)用于水電站優(yōu)化調(diào)度時(shí)體現(xiàn)了良好的性能[3- 6]。作為群體智能進(jìn)化算法,PSO同遺傳算法一樣在求解水電站優(yōu)化調(diào)度時(shí)也存在易“早熟”現(xiàn)象。基于此,學(xué)者們提出了許多典型的改進(jìn)。文獻(xiàn)[3]引入錦標(biāo)賽選擇機(jī)制和線性遞減自適應(yīng)慣性權(quán)重因子,提出了改進(jìn)微粒群算法,在避免算法早熟方面取得了很好的效果。文獻(xiàn)[4]結(jié)合免疫進(jìn)化算法和PSO,提出了局部和全局搜索性能均表現(xiàn)良好的免疫進(jìn)化粒子群算法。就水電站優(yōu)化調(diào)度PSO算法本身改進(jìn)而言,從可獲文獻(xiàn)知傳統(tǒng)研究多采用線性遞減自適應(yīng)策略作為慣性權(quán)重因子的調(diào)整策略[4,5],缺乏對(duì)求解水電站優(yōu)化調(diào)度時(shí)算法慣性權(quán)重因子調(diào)整需要的分析。同時(shí)線性遞減自適應(yīng)策略提出至今已數(shù)10年,目前迫切需要尋求更最適合求解水電站優(yōu)化調(diào)度問(wèn)題的慣性權(quán)重改進(jìn)策略。基于此,筆者進(jìn)行了深入的探討。

    1 水電站日優(yōu)化調(diào)度模型

    水電站水庫(kù)優(yōu)化調(diào)度有“以水定電”和“以電定水”兩種數(shù)學(xué)模型。研究選用“以水定電”模型,即給定調(diào)度期用水量尋求水庫(kù)調(diào)度期總發(fā)電量最優(yōu),在電價(jià)取為常數(shù)時(shí)總發(fā)電量最大即總出力最大。以單個(gè)水電站為例,以小時(shí)為單位進(jìn)行的日發(fā)電優(yōu)化調(diào)度模型描述如下

    (1)

    主要約束條件

    水量平衡約束Vj+1=Vj+(Ij-Qlj-Qsj)Δtj

    (2)

    水位約束Zt,min≤Zt≤Zt,max

    (3)

    Zo=Z′,Ze=Z″

    (4)

    下泄流量約束Qt,min≤Qt≤Qt,max

    (5)

    出力約束Nmin≤Nt≤Nmax

    (6)

    式中,Qt和Ht為水庫(kù)第t時(shí)段發(fā)電流量和水頭;Vj、Vj+1為水庫(kù)時(shí)段初末庫(kù)容;Ij、Qlj和Qsj分別為水庫(kù)時(shí)段入庫(kù)流量、發(fā)電引用流量和棄水流量;Zmin、Zt和Zmax分別為水庫(kù)時(shí)段允許最低水位、時(shí)段水位、水庫(kù)時(shí)段允許最高水位;Zo和Ze為水庫(kù)調(diào)度期初始和期末水位,Z′和Z″為“以水定電”模型中給定的水庫(kù)調(diào)度期初始和期末水位;Qt,min、Qt和Qt,max分別為水庫(kù)第t時(shí)段的最小下泄流量、時(shí)段下泄流量和時(shí)段最大下泄流量;Nmin、Nt和Nmax分別為水庫(kù)的保證出力、第t時(shí)段出力和裝機(jī)容量。

    2 線性微分自適應(yīng)粒子群算法

    鳥(niǎo)類(lèi)搜尋食物的原則是找尋離食物最近的鳥(niǎo)的周?chē)鷧^(qū)域。PSO算法起源于對(duì)鳥(niǎo)類(lèi)捕食的模擬,將每個(gè)優(yōu)化問(wèn)題的解看作是搜索空間的一只鳥(niǎo)(稱之為“粒子”)。PSO中的粒子們都擁有一個(gè)決定自己飛翔方向和距離的速度,并以追隨當(dāng)前最優(yōu)粒子的方式在解空間中搜索。其求解過(guò)程同遺傳算法類(lèi)似均為多點(diǎn)搜索,但不同之處在于遺傳算法通過(guò)遺傳操作(選擇、交叉和變異)進(jìn)化,而PSO則通過(guò)更新粒子速度和位置來(lái)進(jìn)化。下式為粒子群算法中粒子速度和位置更新的迭代式

    vk+1=w×vk+c1(pbestk-xk)+c2(gbestk-xk)

    (7)

    xk+1=xk+vk+1

    (8)

    式中,vk和vk+1依次代表粒子在第k代和第k+1代的速度;xk和xk+1則依次代表粒子在第k代和第k+1代的位置;參數(shù)w表示慣性權(quán)重;參數(shù)c1和c2為學(xué)習(xí)因子,一般取2;pbest表示粒子本身第k代找到的最優(yōu)解的位置;gbest表示整個(gè)種群到第k代為止目前找到的最優(yōu)解的位置。式(7)和式(8)表明粒子速度和位置的更新是依靠局部值極和全局極值的更新來(lái)實(shí)現(xiàn)的。

    2.1 水電站優(yōu)化調(diào)度粒子群算法慣性權(quán)重調(diào)整需要

    PSO算法與其他智能算法的一個(gè)顯著區(qū)別就是在于PSO待調(diào)整的參數(shù)很少,而這些參數(shù)的設(shè)計(jì)對(duì)算法的精度和效率卻存在顯著影響[2]。慣性權(quán)重w是PSO算法的關(guān)鍵參數(shù),它可以消除對(duì)vmax的需要,起到平衡全局搜索和局部搜索的作用。較大的w利于跳出局部極小點(diǎn);而較小的w則有助于算法收斂?;谶@個(gè)原理,Y. Shi和R. C. Eberhart 提出線性遞減權(quán)重策略[7]

    (9)

    式中,wk為取代傳統(tǒng)固定權(quán)重w的線性遞減自適應(yīng)慣性權(quán)重;wmax和wmin依次為慣性權(quán)重的最大值和最小值,典型取值為0.9和0.4[8];k和kmax則依次代表當(dāng)前進(jìn)化代數(shù)和最大進(jìn)化代數(shù)。

    線性遞減權(quán)重策略是目前水電站優(yōu)化調(diào)度粒子群算法中常常采用的策略[4,5],該策略引導(dǎo)下w線性減小,較固定權(quán)重利于算法的搜索和保持算法收斂性。但該策略的w減小是線性的,當(dāng)目標(biāo)函數(shù)為非凸函數(shù)時(shí),若在進(jìn)化初期PSO未能搜索到最好點(diǎn),隨著w的線性減小易陷入局部極值,導(dǎo)致算法早熟。文獻(xiàn)[9]進(jìn)行了水電站水庫(kù)優(yōu)化調(diào)度數(shù)學(xué)模型凸性研究,指出:當(dāng)水電站遇到出力受阻情況時(shí),水電站優(yōu)化調(diào)度目標(biāo)函數(shù)發(fā)生突變,已不是凸函數(shù)。因此采用的線性遞減權(quán)重策略不適合用于求解水電站優(yōu)化調(diào)度問(wèn)題?;诖?,筆者在水電站優(yōu)化調(diào)度粒子群算法求解中,引入線性微分遞減自適應(yīng)慣性權(quán)重

    (10)

    (11)

    求解水電站優(yōu)化調(diào)度時(shí),由式(11)計(jì)算的w在進(jìn)化初期減小趨勢(shì)緩慢,該策略較線性遞減策略全局搜索能力強(qiáng),利于尋找到優(yōu)秀的優(yōu)化種子;而在進(jìn)化后期,由于前期已經(jīng)找到了合適的優(yōu)化種子,隨著w減小趨勢(shì)加快,算法收斂速度加快。

    2.2 水電站優(yōu)化調(diào)度線性微分自適應(yīng)粒子群算法設(shè)計(jì)

    運(yùn)用粒子群算法求解水電站優(yōu)化調(diào)度時(shí),一個(gè)粒子代表水電站的一種運(yùn)行策略。在本研究中每個(gè)粒子位置均代表一組表征水電站水庫(kù)各時(shí)段末水位的向量(Z1,Z2,…,Zn),同時(shí)每個(gè)粒子都擁有一個(gè)代表水庫(kù)各時(shí)段末水位漲落速度的粒子飛行速度向量(v1,v2,…,vn)。PSO同遺傳算法一樣都是通過(guò)適應(yīng)度函數(shù)衡量粒子解的優(yōu)越程度,但具體的尋優(yōu)卻是通過(guò)粒子“飛行經(jīng)驗(yàn)”在全空間搜索達(dá)到的。粒子尋優(yōu)的時(shí)候在解空間中同時(shí)向局部極值和全局極值兩個(gè)點(diǎn)接近。

    改進(jìn)后的線性微分自適應(yīng)粒子群算法求解步驟如下:

    (2) 計(jì)算粒子適應(yīng)度。構(gòu)造如下適應(yīng)度函數(shù)計(jì)算各粒子適應(yīng)度[10]

    (12)

    (3)更新粒子的位置和速度。在進(jìn)行粒子位置和速度更新之前,需要進(jìn)行局部極值和全局極值的更新。局部極值為Step2中計(jì)算得到的最優(yōu)個(gè)體適應(yīng)度,全局極值則為歷代搜索中達(dá)到的最優(yōu)值。粒子的位置和速度更新按照式(7)和式(8),式(7)中慣性權(quán)重的計(jì)算采用式(11)。進(jìn)行粒子位置和速度更新時(shí),更新后的粒子位置依舊要滿足水位約束式(3)和式(4),每一維粒子的速度也都被限制在最大速度vmax內(nèi)。若不滿足,則按照以下方式處理(4)中止條件判斷。水電站優(yōu)化調(diào)度PSO采用是否達(dá)到最大進(jìn)化代數(shù)作為中止條件,達(dá)到則算法中止,優(yōu)化調(diào)度的解為此時(shí)的全局極值;沒(méi)有達(dá)到最大進(jìn)化代數(shù),則轉(zhuǎn)入步驟(2)繼續(xù)循環(huán)。

    (13)

    (14)

    3 實(shí)例應(yīng)用

    葛洲壩樞紐是長(zhǎng)江干流上建設(shè)的第一座大型水利工程,同時(shí)也是三峽樞紐的組成部分,對(duì)三峽電站的日調(diào)節(jié)進(jìn)行反調(diào)節(jié)。葛洲壩水電站水庫(kù)正常蓄水位為66.50 m,死水位為65.50 m;為了滿足航運(yùn)等要求,水庫(kù)最小下泄流量為3 200 m3/s;電站裝機(jī)容量為271.5萬(wàn)kW,保證出力76.8萬(wàn)kW。

    本研究以葛洲壩水電站為例,進(jìn)行了線性微分遞減自適應(yīng)PSO和傳統(tǒng)PSO對(duì)比。假定“以水定電”模型中,葛洲壩水庫(kù)調(diào)度期初始和期末水位均為66.00 m,初始化粒子時(shí)將水位以0.01 m精度離散。兩種PSO算法學(xué)習(xí)因子參數(shù)c1和c2為,均取為2,粒子最大速度vmax取為0.5。鑒于粒子群算法是基于隨機(jī)理論的進(jìn)化算法,故就葛洲壩水電站日優(yōu)化調(diào)度問(wèn)題,分別運(yùn)用兩種算法進(jìn)行5次模擬調(diào)度,并取最好的一次作為各自最終求解的調(diào)度結(jié)果。兩種算法優(yōu)化調(diào)度出力分配對(duì)比見(jiàn)圖1,5次模擬調(diào)度結(jié)果對(duì)比見(jiàn)表1。兩種算法綜合性能對(duì)比見(jiàn)表2。

    圖1 兩種算法出力分配對(duì)比

    表1 兩種算法5次模擬調(diào)度結(jié)果對(duì)比 萬(wàn)kW

    表2 兩種算法綜合性能對(duì)比 萬(wàn)kW

    從表1可以看出,改進(jìn)算法較傳統(tǒng)算法找到了更優(yōu)的解,減小了無(wú)益棄水,增發(fā)了電。分析表1和圖1可得,兩種算法求得的出力分配結(jié)果主要是調(diào)度期初和調(diào)度期末存在差異,改進(jìn)后的算法在調(diào)度期初找到了更加優(yōu)良的解。從表1和表2可以看出,改進(jìn)算法綜合性能明顯優(yōu)于傳統(tǒng)算法,無(wú)論是其“平均表現(xiàn)”、“最優(yōu)表現(xiàn)”還是“最差表現(xiàn)”均優(yōu)于傳統(tǒng)算法。從兩者的優(yōu)劣差對(duì)比可以看出,改進(jìn)算法穩(wěn)定性優(yōu)于傳統(tǒng)算法。

    以上研究證明改進(jìn)算法較傳統(tǒng)算法全局搜索能力強(qiáng),線性微分遞減自適應(yīng)粒子群算法通過(guò)進(jìn)化初期減小趨勢(shì)緩慢的慣性權(quán)重找到優(yōu)秀的優(yōu)化種子,通過(guò)進(jìn)化后期迅速減小的慣性權(quán)重加快收斂速度。這種改進(jìn)有效改善了由于水電站優(yōu)化調(diào)度函數(shù)非凸性帶來(lái)的粒子群算法求解易早熟問(wèn)題。

    4 結(jié) 語(yǔ)

    (1)以參數(shù)簡(jiǎn)潔著稱的粒子群算法在求解水電站優(yōu)化調(diào)度問(wèn)題時(shí),由于無(wú)法克服水電站優(yōu)化調(diào)度目標(biāo)函數(shù)非凸性,存在易早熟問(wèn)題。線性微分遞減自適應(yīng)策略,通過(guò)構(gòu)造進(jìn)化初期減小趨勢(shì)緩慢和進(jìn)化后期減小趨勢(shì)加快的慣性權(quán)重,改善目標(biāo)函數(shù)非凸性的影響,幫助算法尋求到了優(yōu)秀的進(jìn)化種子,提升了算法的探索能力和開(kāi)發(fā)能力。

    (2)雖然粒子群算法待調(diào)整的參數(shù)很少,但這些參數(shù)的設(shè)計(jì)對(duì)算法的精度和效率卻存在顯著影響。針對(duì)不同的優(yōu)化問(wèn)題,傳統(tǒng)算法難免有其局限性。運(yùn)用粒子群求解優(yōu)化問(wèn)題時(shí),應(yīng)在分析優(yōu)化問(wèn)題對(duì)調(diào)整需要的基礎(chǔ)上,選擇合適的改進(jìn)策略。

    [1] 鄭姣, 楊侃, 倪福全, 等.水庫(kù)群發(fā)電優(yōu)化調(diào)度遺傳算法整體改進(jìn)策略研究[J]. 水利學(xué)報(bào), 2013, 44(2): 205- 211.

    [2] 高尚, 楊靜宇. 群智能算法及其應(yīng)用[M]. 北京: 中國(guó)水利水電出版社, 2006.

    [3] 朱鳳霞, 熊立華, 高仕春, 等. 改進(jìn)微粒群算法在梯級(jí)電站長(zhǎng)期優(yōu)化調(diào)度中的應(yīng)用[J]. 水文, 2007, 27 (5): 42- 45.

    [4] 萬(wàn)芳, 原文林, 黃強(qiáng), 等. 基于免疫進(jìn)化算法的粒子群算法在梯級(jí)水庫(kù)優(yōu)化調(diào)度中的應(yīng)用[J]. 水力發(fā)電學(xué)報(bào), 2010, 29 (1): 202- 206.

    [5] 王少波, 解建倉(cāng), 汪妮. 基于改進(jìn)粒子群算法的水電站水庫(kù)優(yōu)化調(diào)度研究[J]. 水力發(fā)電學(xué)報(bào), 2008, 27(3): 13- 15.

    [6] 張俊, 程春田, 廖勝利, 等. 改進(jìn)粒子群優(yōu)化算法在水電站群優(yōu)化調(diào)度中的應(yīng)用研究[J]. 水利學(xué)報(bào), 2009, 40 (4): 435- 439.

    [7] SHI Y, EBERHART R C. Empirical study of Particle Swarm Optimization[C]∥Proceedings of the Congress on Evolutionary Computation, IEEE, 1999: 1945- 1950.

    [8] 胡建秀, 曾建潮. 微粒群算法中慣性權(quán)重的調(diào)整策略[J]. 計(jì)算機(jī)工程, 2007, 33(11): 193- 195.

    [9] 王濤. 水庫(kù)優(yōu)化調(diào)度理論方法及應(yīng)用[D]. 南京: 河海大學(xué), 2009.

    [10] 鄭姣, 楊侃, 郝永懷, 等. 基于穩(wěn)態(tài)繁殖的遺傳算法在水庫(kù)優(yōu)化調(diào)度中的應(yīng)用[J]. 水電能源科學(xué), 2011, 29 (8): 38- 41.

    ApplicationofParticleSwarmOptimizationwithLinearDifferentialDeclineAdaptiveinOperationOptimizationofHydropowerStation

    WAN Zhonghai, YE Shengjin, ZHENG Jiao

    (Sichuan Dujiangyan Dongfengqu Management Department, Chengdu 610072, Sichuan, China)

    For improving the Particle Swarm Optimization in solving hydropower station optimal scheduling, the Particle Swarm Optimization with Linear Differential Decline Adaptive is proposed by analyzing the adjustment need of key parameters’ inertia weight. The improved method can increase the exploration ability to kip local optimal solution by slowly reducing inertia weight in earlier stage, and can improve the development ability to accelerate algorithm convergence through quickly reducing inertia weight in later stage. Taking the optimal operation of Gezhouba Hydropower Station as an example, the improved method and traditional algorithm are compared. The results show that the Linear Differential Decline Adaptive strategy can enhance the search capability and stability of algorithm and the proposed method can improve the premature problem of Particle Swarm Optimization caused by the non-convex objection function of optimal scheduling. The results provide a new idea for improving the inertia weight of Particle Swarm Optimization in hydropower station optimal scheduling.

    optimal operation of hydropower station; Particle Swarm Optimization; Inertia Weight; Linear Differential Decline

    TV697.1

    A

    0559- 9342(2017)09- 0085- 04

    2017- 03- 25

    萬(wàn)忠海(1963—),男,四川渠縣人,高級(jí)工程師,主要從事水利工程建設(shè)管理和水資源管理工作.

    (責(zé)任編輯高 瑜)

    猜你喜歡
    微分慣性水電站
    你真的了解慣性嗎
    里底水電站工程
    西北水電(2022年1期)2022-04-16 05:20:06
    沖破『慣性』 看慣性
    白鶴灘水電站有多“?!?
    金橋(2021年8期)2021-08-23 01:06:54
    擬微分算子在Hp(ω)上的有界性
    上下解反向的脈沖微分包含解的存在性
    無(wú)處不在的慣性
    小型水電站技術(shù)改造探討
    普遍存在的慣性
    借助微分探求連續(xù)函數(shù)的極值點(diǎn)
    十分钟在线观看高清视频www| 老熟妇乱子伦视频在线观看| 嫩草影视91久久| 91麻豆精品激情在线观看国产 | 91精品三级在线观看| 啪啪无遮挡十八禁网站| 首页视频小说图片口味搜索| 国产成人精品久久二区二区免费| 亚洲精品国产区一区二| 国产伦理片在线播放av一区| 99riav亚洲国产免费| 一级片'在线观看视频| 欧美乱码精品一区二区三区| 日本五十路高清| 久久久久精品人妻al黑| 麻豆国产av国片精品| 搡老乐熟女国产| 手机成人av网站| 在线观看人妻少妇| 露出奶头的视频| 国产av又大| 国产精品欧美亚洲77777| 亚洲va日本ⅴa欧美va伊人久久| av线在线观看网站| 国产男女超爽视频在线观看| aaaaa片日本免费| 一级毛片精品| 大型黄色视频在线免费观看| 女人久久www免费人成看片| 国产av又大| 天堂8中文在线网| 欧美成人午夜精品| 亚洲,欧美精品.| 一级,二级,三级黄色视频| av又黄又爽大尺度在线免费看| 丝袜美足系列| 大型黄色视频在线免费观看| 亚洲性夜色夜夜综合| 一进一出抽搐动态| 亚洲中文日韩欧美视频| 国产在线一区二区三区精| 国产在线免费精品| 无人区码免费观看不卡 | 国内毛片毛片毛片毛片毛片| 在线av久久热| 欧美国产精品va在线观看不卡| 国产男靠女视频免费网站| 国产成人精品在线电影| 人人妻人人澡人人看| 99精品在免费线老司机午夜| 亚洲专区中文字幕在线| 国产亚洲欧美在线一区二区| 女人精品久久久久毛片| 欧美国产精品va在线观看不卡| 欧美精品高潮呻吟av久久| 欧美激情高清一区二区三区| 可以免费在线观看a视频的电影网站| 王馨瑶露胸无遮挡在线观看| 国产精品影院久久| 少妇裸体淫交视频免费看高清 | 天天添夜夜摸| 国产成人欧美在线观看 | 免费少妇av软件| 中文字幕人妻丝袜制服| 亚洲第一欧美日韩一区二区三区 | 一二三四社区在线视频社区8| 国产成人影院久久av| 在线播放国产精品三级| 精品国产亚洲在线| 成年人免费黄色播放视频| 麻豆国产av国片精品| 亚洲 欧美一区二区三区| 丝袜在线中文字幕| 欧美激情极品国产一区二区三区| 午夜老司机福利片| 日韩免费高清中文字幕av| 久久久国产欧美日韩av| 大型av网站在线播放| 亚洲avbb在线观看| 日韩一区二区三区影片| 久久久水蜜桃国产精品网| av福利片在线| 91九色精品人成在线观看| 亚洲少妇的诱惑av| 久久99一区二区三区| 777米奇影视久久| 五月天丁香电影| 国产又色又爽无遮挡免费看| 一级片免费观看大全| 欧美国产精品一级二级三级| 热re99久久精品国产66热6| 一本一本久久a久久精品综合妖精| 黄片大片在线免费观看| 99精品在免费线老司机午夜| 亚洲成人免费av在线播放| 国内毛片毛片毛片毛片毛片| 欧美国产精品一级二级三级| 久久婷婷成人综合色麻豆| 美女高潮到喷水免费观看| 99久久国产精品久久久| 两性午夜刺激爽爽歪歪视频在线观看 | 久久精品国产亚洲av香蕉五月 | 自拍欧美九色日韩亚洲蝌蚪91| 欧美 亚洲 国产 日韩一| 国产真人三级小视频在线观看| 欧美精品人与动牲交sv欧美| 国产精品一区二区在线不卡| 不卡一级毛片| 天天添夜夜摸| 成人18禁在线播放| 国产亚洲精品一区二区www | 亚洲第一青青草原| 国产福利在线免费观看视频| 十八禁网站网址无遮挡| 成人国语在线视频| 黄色成人免费大全| 91麻豆精品激情在线观看国产 | 夜夜爽天天搞| 国产日韩一区二区三区精品不卡| 国产成人欧美| 黄色视频不卡| 一本色道久久久久久精品综合| 99国产精品一区二区三区| 国产淫语在线视频| tocl精华| av电影中文网址| 搡老岳熟女国产| 久久亚洲真实| 国产精品一区二区在线不卡| 成人国产一区最新在线观看| 国产福利在线免费观看视频| 一级毛片女人18水好多| 在线观看人妻少妇| 在线观看免费日韩欧美大片| 亚洲av片天天在线观看| 日韩精品免费视频一区二区三区| 国产成人影院久久av| 国产日韩一区二区三区精品不卡| 99久久人妻综合| 免费女性裸体啪啪无遮挡网站| 欧美成人免费av一区二区三区 | 超碰97精品在线观看| 欧美激情极品国产一区二区三区| 777米奇影视久久| 欧美成人午夜精品| 99精国产麻豆久久婷婷| 国产高清国产精品国产三级| 韩国精品一区二区三区| 欧美在线一区亚洲| 制服人妻中文乱码| 在线观看人妻少妇| 欧美成狂野欧美在线观看| 美女福利国产在线| 国产在线观看jvid| 欧美乱妇无乱码| 精品少妇内射三级| 久久精品国产亚洲av香蕉五月 | av国产精品久久久久影院| av线在线观看网站| 最新的欧美精品一区二区| 一级毛片女人18水好多| 午夜福利欧美成人| 别揉我奶头~嗯~啊~动态视频| 咕卡用的链子| 1024香蕉在线观看| 国产麻豆69| 欧美激情久久久久久爽电影 | 亚洲欧美一区二区三区久久| 久久毛片免费看一区二区三区| 久久午夜亚洲精品久久| 亚洲成人手机| 久久毛片免费看一区二区三区| 69精品国产乱码久久久| www日本在线高清视频| 国产不卡av网站在线观看| 最黄视频免费看| 国产成+人综合+亚洲专区| 国产av一区二区精品久久| 一区在线观看完整版| 午夜91福利影院| 国产淫语在线视频| 日本wwww免费看| 亚洲一区中文字幕在线| 国产在线免费精品| 久久久精品94久久精品| 国产在线一区二区三区精| av超薄肉色丝袜交足视频| 日日夜夜操网爽| 无遮挡黄片免费观看| 国产黄色免费在线视频| 99久久99久久久精品蜜桃| 久久精品国产亚洲av香蕉五月 | 少妇猛男粗大的猛烈进出视频| 久久久国产欧美日韩av| 国产熟女午夜一区二区三区| 九色亚洲精品在线播放| 精品久久久久久久毛片微露脸| 大码成人一级视频| 国产精品久久久久久精品古装| 国产一区二区三区在线臀色熟女 | aaaaa片日本免费| av免费在线观看网站| 性高湖久久久久久久久免费观看| 夫妻午夜视频| 97在线人人人人妻| 国产精品成人在线| 女性被躁到高潮视频| 美女高潮喷水抽搐中文字幕| av又黄又爽大尺度在线免费看| 亚洲精品在线美女| 狠狠狠狠99中文字幕| 久久青草综合色| 国产一区有黄有色的免费视频| 国产黄频视频在线观看| 91九色精品人成在线观看| 18禁国产床啪视频网站| 亚洲第一欧美日韩一区二区三区 | 国产精品 国内视频| 少妇猛男粗大的猛烈进出视频| 一区福利在线观看| 亚洲成a人片在线一区二区| 久久国产精品影院| 日韩视频一区二区在线观看| 黄片播放在线免费| 久久久久久久国产电影| 久久精品亚洲av国产电影网| 欧美激情久久久久久爽电影 | 欧美性长视频在线观看| 性高湖久久久久久久久免费观看| 最近最新中文字幕大全电影3 | 久久亚洲精品不卡| 免费不卡黄色视频| 国产亚洲一区二区精品| 大型av网站在线播放| 男男h啪啪无遮挡| 考比视频在线观看| 可以免费在线观看a视频的电影网站| 另类精品久久| 母亲3免费完整高清在线观看| 亚洲精品av麻豆狂野| 成年动漫av网址| 久久这里只有精品19| 少妇 在线观看| 黑丝袜美女国产一区| 国产精品一区二区在线不卡| 午夜福利乱码中文字幕| 欧美一级毛片孕妇| 日韩欧美国产一区二区入口| 午夜精品久久久久久毛片777| 成年版毛片免费区| 日韩视频一区二区在线观看| 精品一区二区三区av网在线观看 | 狠狠婷婷综合久久久久久88av| av又黄又爽大尺度在线免费看| 亚洲专区国产一区二区| 狠狠婷婷综合久久久久久88av| 亚洲精品久久午夜乱码| 国产精品国产av在线观看| 最新在线观看一区二区三区| 日韩免费av在线播放| a级毛片在线看网站| 王馨瑶露胸无遮挡在线观看| 大香蕉久久网| 午夜免费鲁丝| 91麻豆精品激情在线观看国产 | 成年人黄色毛片网站| 色老头精品视频在线观看| 女人久久www免费人成看片| 国产人伦9x9x在线观看| 国产精品久久久av美女十八| 精品久久久久久电影网| 丝袜在线中文字幕| 香蕉国产在线看| 欧美激情久久久久久爽电影 | 国产精品久久久久久精品电影小说| 久久久国产精品麻豆| 亚洲精品国产色婷婷电影| 在线观看人妻少妇| 桃花免费在线播放| 免费看a级黄色片| 最黄视频免费看| 黄网站色视频无遮挡免费观看| 久久精品aⅴ一区二区三区四区| 国产三级黄色录像| 不卡一级毛片| 国产一区有黄有色的免费视频| 国产精品二区激情视频| 99久久国产精品久久久| 热99国产精品久久久久久7| 人妻 亚洲 视频| 欧美精品一区二区免费开放| 欧美另类亚洲清纯唯美| 97人妻天天添夜夜摸| 无遮挡黄片免费观看| 久久中文字幕一级| 国产一区二区在线观看av| 天天躁日日躁夜夜躁夜夜| 精品一区二区三区视频在线观看免费 | 香蕉久久夜色| 日韩成人在线观看一区二区三区| 极品教师在线免费播放| 国产单亲对白刺激| 手机成人av网站| 精品人妻熟女毛片av久久网站| 熟女少妇亚洲综合色aaa.| 色在线成人网| 两人在一起打扑克的视频| 女人被躁到高潮嗷嗷叫费观| 亚洲色图 男人天堂 中文字幕| 国产精品国产高清国产av | 日本vs欧美在线观看视频| 国产高清视频在线播放一区| 男人操女人黄网站| 另类精品久久| 亚洲天堂av无毛| 国产精品亚洲一级av第二区| 十八禁人妻一区二区| 激情在线观看视频在线高清 | 国产单亲对白刺激| 亚洲国产毛片av蜜桃av| 国产精品亚洲av一区麻豆| 大型黄色视频在线免费观看| 午夜老司机福利片| av线在线观看网站| 性少妇av在线| 亚洲成人免费av在线播放| 国产一区二区在线观看av| 日韩免费av在线播放| 成人黄色视频免费在线看| 黄网站色视频无遮挡免费观看| 日韩人妻精品一区2区三区| 欧美性长视频在线观看| 大型av网站在线播放| 欧美 亚洲 国产 日韩一| 国内毛片毛片毛片毛片毛片| 又黄又粗又硬又大视频| 少妇裸体淫交视频免费看高清 | 欧美黄色淫秽网站| 十八禁高潮呻吟视频| 无人区码免费观看不卡 | 老司机午夜福利在线观看视频 | 在线观看免费午夜福利视频| 久久久久久久大尺度免费视频| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美日韩另类电影网站| 嫁个100分男人电影在线观看| 亚洲色图av天堂| 51午夜福利影视在线观看| 窝窝影院91人妻| 精品人妻在线不人妻| 亚洲成人免费电影在线观看| 日韩 欧美 亚洲 中文字幕| 亚洲精品一二三| 国产黄色免费在线视频| 国产精品久久久久成人av| 午夜福利乱码中文字幕| 黄色片一级片一级黄色片| 国产精品久久久久久精品电影小说| 男女床上黄色一级片免费看| 自拍欧美九色日韩亚洲蝌蚪91| 老汉色∧v一级毛片| 国产一区二区三区在线臀色熟女 | 亚洲欧美色中文字幕在线| tocl精华| a级毛片黄视频| 黑丝袜美女国产一区| 久久久精品94久久精品| 在线亚洲精品国产二区图片欧美| 51午夜福利影视在线观看| 一级毛片女人18水好多| 免费少妇av软件| 欧美日韩成人在线一区二区| 窝窝影院91人妻| 国产精品免费一区二区三区在线 | 黄色怎么调成土黄色| 天堂动漫精品| 一级毛片精品| 极品教师在线免费播放| 国产精品熟女久久久久浪| 亚洲一码二码三码区别大吗| 777久久人妻少妇嫩草av网站| 电影成人av| 久热爱精品视频在线9| 久久久久国内视频| 国产成人欧美| 日本wwww免费看| 国产日韩欧美亚洲二区| 色综合欧美亚洲国产小说| av欧美777| 久久久精品94久久精品| 亚洲久久久国产精品| 制服人妻中文乱码| 在线天堂中文资源库| 99国产精品99久久久久| 天天添夜夜摸| 国产老妇伦熟女老妇高清| 国产福利在线免费观看视频| 麻豆成人av在线观看| 嫁个100分男人电影在线观看| 精品国产亚洲在线| 一进一出抽搐动态| tube8黄色片| 老司机午夜福利在线观看视频 | 亚洲伊人久久精品综合| 天天躁狠狠躁夜夜躁狠狠躁| 97在线人人人人妻| 99热国产这里只有精品6| 18禁观看日本| 99精品久久久久人妻精品| 亚洲精品一卡2卡三卡4卡5卡| 国产亚洲欧美在线一区二区| 国产精品 欧美亚洲| 午夜视频精品福利| 少妇裸体淫交视频免费看高清 | 一级毛片电影观看| 女人被躁到高潮嗷嗷叫费观| 午夜福利在线观看吧| 香蕉久久夜色| 中文字幕人妻丝袜一区二区| 极品人妻少妇av视频| 亚洲欧美激情在线| 丝袜美腿诱惑在线| xxxhd国产人妻xxx| 色在线成人网| 美女主播在线视频| av电影中文网址| 天堂动漫精品| 亚洲欧洲日产国产| 国产免费av片在线观看野外av| 久久久久精品国产欧美久久久| 91麻豆精品激情在线观看国产 | 免费av中文字幕在线| 亚洲视频免费观看视频| 精品少妇久久久久久888优播| 亚洲精品美女久久av网站| 国产免费视频播放在线视频| 18禁美女被吸乳视频| 国产精品久久久久久精品电影小说| 亚洲精品自拍成人| 美女视频免费永久观看网站| 夜夜夜夜夜久久久久| 男女无遮挡免费网站观看| 色视频在线一区二区三区| 黄色视频不卡| 久久精品亚洲精品国产色婷小说| 久久免费观看电影| 欧美日韩亚洲高清精品| 一本色道久久久久久精品综合| 国产不卡一卡二| 久久久久精品国产欧美久久久| 黄片小视频在线播放| 免费一级毛片在线播放高清视频 | 麻豆成人av在线观看| 欧美日韩视频精品一区| 国产在线一区二区三区精| 婷婷丁香在线五月| 91麻豆精品激情在线观看国产 | 91大片在线观看| 十八禁网站网址无遮挡| 亚洲专区字幕在线| 丝袜人妻中文字幕| www.自偷自拍.com| 自拍欧美九色日韩亚洲蝌蚪91| 久久久欧美国产精品| 午夜免费鲁丝| 香蕉丝袜av| 又紧又爽又黄一区二区| bbb黄色大片| 巨乳人妻的诱惑在线观看| www.熟女人妻精品国产| av有码第一页| 国产精品国产高清国产av | 青青草视频在线视频观看| 午夜福利免费观看在线| www日本在线高清视频| 欧美日韩成人在线一区二区| 伊人久久大香线蕉亚洲五| 国产成人欧美在线观看 | 十八禁人妻一区二区| 亚洲国产成人一精品久久久| 一本色道久久久久久精品综合| 久热这里只有精品99| 国产成人av激情在线播放| 大陆偷拍与自拍| 国产精品久久久久久精品电影小说| 亚洲五月婷婷丁香| 极品少妇高潮喷水抽搐| 亚洲欧洲精品一区二区精品久久久| 亚洲精品美女久久久久99蜜臀| 午夜91福利影院| av超薄肉色丝袜交足视频| 每晚都被弄得嗷嗷叫到高潮| 一本—道久久a久久精品蜜桃钙片| 久久精品熟女亚洲av麻豆精品| 成人影院久久| 久久国产精品大桥未久av| 九色亚洲精品在线播放| 一进一出抽搐动态| 久久亚洲精品不卡| 少妇猛男粗大的猛烈进出视频| 天堂中文最新版在线下载| 黑人巨大精品欧美一区二区mp4| 亚洲一卡2卡3卡4卡5卡精品中文| 免费黄频网站在线观看国产| 久久人妻av系列| 国产精品自产拍在线观看55亚洲 | 国产精品久久久人人做人人爽| 窝窝影院91人妻| 最近最新中文字幕大全免费视频| 成人特级黄色片久久久久久久 | 每晚都被弄得嗷嗷叫到高潮| 大码成人一级视频| 捣出白浆h1v1| 国产精品国产高清国产av | 久久 成人 亚洲| 国产成人精品无人区| 丝袜在线中文字幕| 色综合欧美亚洲国产小说| av线在线观看网站| av片东京热男人的天堂| 在线永久观看黄色视频| 亚洲国产欧美网| 啪啪无遮挡十八禁网站| 亚洲 国产 在线| 2018国产大陆天天弄谢| 国产一区二区激情短视频| 欧美精品啪啪一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 啦啦啦免费观看视频1| 精品少妇久久久久久888优播| 中文字幕另类日韩欧美亚洲嫩草| 国产视频一区二区在线看| 99香蕉大伊视频| 精品亚洲乱码少妇综合久久| 亚洲免费av在线视频| 97在线人人人人妻| 国产成人免费观看mmmm| 亚洲精品在线观看二区| av片东京热男人的天堂| 巨乳人妻的诱惑在线观看| 高潮久久久久久久久久久不卡| 久久中文字幕人妻熟女| 亚洲精品在线美女| 精品欧美一区二区三区在线| 日韩三级视频一区二区三区| 成年人午夜在线观看视频| 后天国语完整版免费观看| 欧美在线黄色| 黄色毛片三级朝国网站| 亚洲av欧美aⅴ国产| 亚洲全国av大片| 丰满人妻熟妇乱又伦精品不卡| 国产三级黄色录像| 变态另类成人亚洲欧美熟女 | 国产日韩欧美在线精品| 高潮久久久久久久久久久不卡| 亚洲av第一区精品v没综合| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产中文字幕在线视频| 中文字幕人妻丝袜一区二区| 久久久久网色| 人人妻人人澡人人看| 国产精品国产av在线观看| 欧美乱码精品一区二区三区| 两性夫妻黄色片| 精品国产一区二区三区四区第35| 久久青草综合色| 母亲3免费完整高清在线观看| 自线自在国产av| 男女高潮啪啪啪动态图| 国产一区二区三区在线臀色熟女 | 欧美激情久久久久久爽电影 | avwww免费| 日日夜夜操网爽| 亚洲精华国产精华精| 99久久人妻综合| 国产一区二区 视频在线| 国产免费福利视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 国产精品久久电影中文字幕 | 久久久久久久久久久久大奶| 亚洲伊人色综图| 久久人人爽av亚洲精品天堂| 亚洲av第一区精品v没综合| 久久久久久久国产电影| 成人影院久久| 日本a在线网址| 女人被躁到高潮嗷嗷叫费观| 一级毛片电影观看| 中文字幕另类日韩欧美亚洲嫩草| 久久婷婷成人综合色麻豆| 午夜成年电影在线免费观看| 男女免费视频国产| 中文字幕制服av| svipshipincom国产片| 黑丝袜美女国产一区| 欧美亚洲日本最大视频资源| 老司机午夜福利在线观看视频 | 国产日韩一区二区三区精品不卡| 性色av乱码一区二区三区2| 飞空精品影院首页| 黄色丝袜av网址大全| 国产精品av久久久久免费| 超碰97精品在线观看| 国产精品欧美亚洲77777| 大陆偷拍与自拍| 每晚都被弄得嗷嗷叫到高潮| 丰满人妻熟妇乱又伦精品不卡| 国产免费现黄频在线看| 在线观看www视频免费| 一级片'在线观看视频| 亚洲av日韩精品久久久久久密| 亚洲中文字幕日韩| 亚洲国产欧美在线一区| 一个人免费在线观看的高清视频|