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

    基于DE-CSO混合算法的多目標(biāo)無功優(yōu)化

    2020-11-06 09:03:30盧維佳黃海濤魏明磊
    黑龍江電力 2020年3期
    關(guān)鍵詞:父代交叉種群

    盧維佳,黃海濤,魏明磊

    (1.廣東電網(wǎng)汕尾供電局,廣東 汕尾 516600; 2.國網(wǎng)安徽省電力有限公司定遠(yuǎn)供電局,安徽 滁州 233200)

    0 引 言

    多目標(biāo)無功潮流優(yōu)化是一個(gè)多變量、多約束混合的非線性規(guī)劃問題[1],它通過調(diào)整系統(tǒng)中發(fā)電機(jī)端電壓、變壓器分接頭及無功補(bǔ)償設(shè)備的投切來改變無功潮流,從而降低系統(tǒng)的有功損耗、提高電壓穩(wěn)定裕度,同時(shí)使系統(tǒng)的節(jié)點(diǎn)電壓保持在額定范圍內(nèi)。常規(guī)的無功優(yōu)化求解智能方法已經(jīng)有很多,包括粒子群算法[2-4]、差分進(jìn)化算法[5]、遺傳算法[6-8]等。但這些算法自身在優(yōu)化過程容易導(dǎo)致粒子早熟,陷入局部最優(yōu)。另外,單一的智能優(yōu)化算法不能夠提高粒子的多樣性,優(yōu)化后的多目標(biāo)函數(shù)可行解分布不夠均勻、廣泛,限制了不同需求者對多目標(biāo)無功優(yōu)化可行解的選擇。

    本文提出一種融合CSO算法[9-11]和DE算法的DE-CSO算法。CSO算法參數(shù)少,調(diào)整策略相對單一,但該算法通過縱橫兩部分交叉操作產(chǎn)生的新種群不容易陷入早熟收斂。DE算法收斂速度快,種群粒子多樣性好,但是該算法容易導(dǎo)致早熟收斂,收斂精度不夠。而將兩種算法相結(jié)合的DE-CSO算法不僅保留了CSO算法種群粒子不易陷入早熟收斂的特點(diǎn),而且通過DE-CSO算法中的變異交叉操作,能夠有效地提高粒子多樣性,對于無功優(yōu)化這樣的非線性問題求解有很好的適用性。另外,針對多目標(biāo)函數(shù)優(yōu)化問題,本文通過快速非支配排序和計(jì)算粒子擁擠距離的方法來實(shí)現(xiàn)對粒子的優(yōu)勝劣汰,產(chǎn)生均勻、廣泛的Pareto前沿分布。

    1 多目標(biāo)無功優(yōu)化的數(shù)學(xué)模型

    1.1 多目標(biāo)函數(shù)

    多目標(biāo)無功優(yōu)化數(shù)學(xué)模型是在滿足系統(tǒng)和發(fā)電機(jī)組運(yùn)行約束的前提下,通過改變發(fā)電機(jī)端電壓、變壓器分接頭及無功補(bǔ)償設(shè)備的投切,使系統(tǒng)的有功網(wǎng)損及電壓偏移最小,電壓穩(wěn)定裕度最大。

    1)系統(tǒng)網(wǎng)損最小

    式中:Ploss為系統(tǒng)的有功網(wǎng)損;NE為系統(tǒng)總支路數(shù);Gk(i、j)為支路k電導(dǎo);Vi、Vj、δij分別為支路k兩端節(jié)點(diǎn)i、j的電壓幅值及其相角差。

    2)電壓偏差量最小

    式中:NPQ為系統(tǒng)負(fù)荷節(jié)點(diǎn)總數(shù);Vspc為節(jié)點(diǎn)i額定電壓;Vimax、Vimin分別為節(jié)點(diǎn)i的最大、最小允許偏差。

    3)靜態(tài)電壓穩(wěn)定裕度最大

    靜態(tài)電壓穩(wěn)定裕度是指系統(tǒng)當(dāng)前運(yùn)行點(diǎn)到電壓崩潰點(diǎn)間的距離。本文采用最小特征值模[2]作為靜態(tài)電壓穩(wěn)定裕度指標(biāo),這是因?yàn)樽钚√卣髦的D軌蛴沙绷餮趴杀染仃囍苯忧蟮?,方便算法運(yùn)行。當(dāng)電力系統(tǒng)由正常工作點(diǎn)向電壓崩潰點(diǎn)過渡時(shí),系統(tǒng)雅可比矩陣最小特征值模單調(diào)地趨向于0。為了使多目標(biāo)函數(shù)中每個(gè)函數(shù)都具有最小化的形式,對靜態(tài)電壓穩(wěn)定裕度取負(fù)。因此,靜態(tài)電壓穩(wěn)定裕度最大的目標(biāo)函數(shù)為

    minF3=-(min(|eig(J)|))

    式中:J為收斂潮流雅可比矩陣;eig(J)為雅可比矩陣的所有特征值。

    綜上所述,本文建立的目標(biāo)函數(shù)為

    1.2 約束條件

    1)等式約束條件

    式中:PGi、QGi分別為節(jié)點(diǎn)i注入的有功功率和無功功率;Pdi、Qdi分別為節(jié)點(diǎn)i的負(fù)荷有功功率和無功功率;ui、uj分別為節(jié)點(diǎn)i、j的電壓;n為系統(tǒng)節(jié)點(diǎn)總數(shù);Qci為節(jié)點(diǎn)i的電容器無功補(bǔ)償量;Gij、Bij分別為節(jié)點(diǎn)i、j之間的電導(dǎo)、電納;θij為節(jié)點(diǎn)i、j的電壓相角差。

    2)不等式約束條件:

    式中:Tk為變壓器分接頭檔位。

    2 DE-CSO優(yōu)化算法及在多目標(biāo)無功優(yōu)化中的應(yīng)用

    2.1 CSO算法原理

    假設(shè)CSO算法的種群用矩陣B表示,矩陣中的每一行用B(i)表示,代表問題的一個(gè)解,B(i,j)表示矩陣中的每一個(gè)元素,種群規(guī)模M即矩陣的行數(shù),解空間維數(shù)D即矩陣列數(shù)。橫向交叉和縱向交叉得到的解稱為折中解,分別用MShc和MSvc表示。

    1)橫向交叉

    橫向交叉是在種群中兩個(gè)不同個(gè)體所有維之間進(jìn)行的一種算數(shù)交叉。父代個(gè)體X(i)和X(j)的第d維進(jìn)行橫向交叉,產(chǎn)生的第d維后代的折中解MShc(i,d)和MShc(j,d)采用式(1)進(jìn)行計(jì)算:

    (1)

    其中:r1、r2∈[0,1];c1∈[-1,1];X(i,d)、X(j,d)分別是父代種群X(i)和X(j)的第d維。

    隨機(jī)抽取矩陣DSvc中不同行的2個(gè)粒子,抽取共進(jìn)行M/2次,被選中的個(gè)體(如X(no1)和X(no2))通過式(1)產(chǎn)生他們的后代(MShc(no1)和MShc(no2))。

    2)縱向交叉

    縱向交叉是所有個(gè)體都會參與自身個(gè)體不同維之間的一種算數(shù)交叉。假定父代個(gè)體X(i)的第d1維和第d2維進(jìn)行縱向交叉,產(chǎn)生的第d1維后代[12]MSvc(i,d1)采用式(2)進(jìn)行計(jì)算:

    MSvc(i,d1)=r·X(i,d1)+(1-r)·X(i,d2)

    (2)

    式中:d1、d2∈[1,D];i∈[1,M];r∈[0,1]。

    由于縱向交叉法是對于所有種群不同維的運(yùn)算操作,故先對父代種群即經(jīng)橫向交叉后保留下來的矩陣DShc進(jìn)行歸一化處理,然后隨機(jī)抽取不同行的2列數(shù)據(jù),每個(gè)列數(shù)有且只有一次被選中(如d1列和d2列),再通過式(2)產(chǎn)生新的子代(MSvc(i,d1)),最后對新的種群進(jìn)行反歸一化處理??v向交叉抽取共進(jìn)行D/2次,只改變D/2列的數(shù)據(jù),其目的是促使可能陷入停滯不前的維(如d1)跳出局部最優(yōu),同時(shí)不破壞另外可能正常進(jìn)化的維(如d2)。

    2.2 DE算法原理

    DE算法主要通過變異操作和交叉操作來產(chǎn)生新的后代。

    1)變異操作

    DE算法從群體中任意選取3個(gè)互不相同的染色體進(jìn)行加權(quán)變異,產(chǎn)生子代染色體,計(jì)算式如下:

    yi=xno1+F·(xno2-xno3)

    (3)

    式中:xno1、xno2、xno3為3個(gè)互不相同的染色體;F為縮放因子。

    標(biāo)準(zhǔn)的DE算法中F為固定數(shù)值,F(xiàn)值對算法的全局尋優(yōu)能力有很大影響。F選擇過小,算法容易陷入早熟收斂;F選擇過大,算法收斂速度會變慢。為了能夠解決易早熟和迭代收斂速度較慢的問題,提出一種自適應(yīng)縮放因子Fi,其計(jì)算式為

    2)交叉操作

    DE算法將父代粒子維度按一定概率與變異粒子維度執(zhí)行交叉操作,計(jì)算式如下:

    (4)

    式中:CR為交叉概率;yi,j是變異操作后第i條染色體的第j個(gè)基因。

    DE算法通過交叉操作保留父代粒子的特征,使算法在優(yōu)化過程中能夠有效地保持粒子多樣性。

    2.3 多目標(biāo)函數(shù)的處理

    由于多目標(biāo)無功優(yōu)化的目標(biāo)函數(shù)之間存在相互矛盾性,因此很難找到一個(gè)可以使各目標(biāo)函數(shù)同時(shí)達(dá)到最優(yōu)的解,如何評價(jià)這些解的優(yōu)劣性是該問題求解的難點(diǎn)之一。目前,絕大部分算法采用加權(quán)疊加比較的評價(jià)方法,但是該方法不能很好地反應(yīng)各目標(biāo)函數(shù)的優(yōu)劣性。本文通過多目標(biāo)快速非支配排序法[12-13]和比較粒子擁擠距離的方式得到所有解中的Pareto最優(yōu)前沿[14],并采用模糊隸屬度尋找最優(yōu)解。

    1)快速非支配排序法

    對于目標(biāo)函數(shù)fi的兩個(gè)可行解a和b,滿足任意fi(a)≤fi(b)(i=1,2,…,n)且存在fi使得fi(a)

    2)擁擠距離

    (5)

    式中:f[i]p為第i個(gè)種群中第p個(gè)目標(biāo)函數(shù)中的解。

    比較每個(gè)粒子的擁擠距離,擁擠距離越大,粒子具有越好的多樣性和尋優(yōu)性,經(jīng)過多次迭代,能夠使Pareto前沿分布更加均勻。

    2.4 多目標(biāo)DE-CSO優(yōu)化算法及應(yīng)用

    針對CSO算法和DE算法在優(yōu)化過程中的缺點(diǎn),DE-CSO優(yōu)化算法混合了CSO算法中縱向交叉操作、橫向交叉操作及DE算法中變異、交叉操作。在更新粒子種群時(shí),分別對父代粒子進(jìn)行變異、交叉和縱橫交叉操作,然后再將產(chǎn)生的子代粒子與父代粒子進(jìn)行比較,從而實(shí)現(xiàn)粒子的優(yōu)勝劣汰。通過多種更新子代粒子的方式,DE-CSO算法能夠融合兩種算法的優(yōu)點(diǎn),使新種群能夠跳出局部最優(yōu),提高收斂精度和粒子多樣性。為了將DE-CSO算法應(yīng)用到多目標(biāo)函數(shù)中,混合算法不保留基本縱橫交叉算法中競爭算子的運(yùn)算,而是采用多目標(biāo)快速非支配排序法和比較粒子擁擠距離的方式來評價(jià)子代粒子和父代粒子的優(yōu)劣性。通過多次迭代后,粒子的Pareto前沿分布更加廣泛、均勻。將DE-CSO算法應(yīng)用在多目標(biāo)無功優(yōu)化問題上,其算法步驟如下所述,流程圖如圖1所示。

    圖1 多目標(biāo)DE-CSO優(yōu)化算法流程圖

    1)設(shè)置算法參數(shù)、控制變量,初始種群為F1。

    2)將種群F1粒子帶入式(3)、(4)進(jìn)行粒子變異和交叉操作,產(chǎn)生子代種群,并合并父代和子代,形成種群規(guī)模為2n的混合種群F2。

    4)運(yùn)用快速非支配排序法和計(jì)算粒子擁擠距離的擇優(yōu)策略從混合種群F2中選擇種群數(shù)為n的優(yōu)秀粒子,形成新種群F3。

    5)執(zhí)行CSO算法運(yùn)算。將種群F3的粒子帶入式(1)、(3)進(jìn)行橫向交叉操作和縱向交叉操作,產(chǎn)生種群F6。

    6)合并混合種群F3和F6,執(zhí)行步驟4)產(chǎn)生新種群F7,并用種群F7覆蓋種群F1。

    7)更新Pareto前沿分布,檢驗(yàn)是否滿足迭代次數(shù),滿足則退出運(yùn)行,輸出結(jié)果;否則轉(zhuǎn)回步驟2)。

    3 算例分析

    本文使用IEEE14、IEEE57標(biāo)準(zhǔn)節(jié)點(diǎn)系統(tǒng)進(jìn)行仿真實(shí)驗(yàn)。算法參數(shù)設(shè)置如下:種群規(guī)模為50,最大迭代次數(shù)為200,CSO算法中的橫向交叉概率為1,縱向交叉概率為0.6,差分算法的交叉率為 0.5,變異率為 0.5。在測試中除平衡節(jié)點(diǎn)外,電源節(jié)點(diǎn)都作為PV節(jié)點(diǎn)處理,負(fù)荷節(jié)點(diǎn)都作為PQ節(jié)點(diǎn)處理。變壓器變比上限為1.1,下限為0.95,變量步長為0.05;無功補(bǔ)償為0.0~0.5 p.u.,調(diào)節(jié)步長為0.1 p.u.;發(fā)電機(jī)節(jié)點(diǎn)電壓上、下極限分別為1.1、0.9,變量步長為0.05。

    3.1 IEEE14節(jié)點(diǎn)算例

    該系統(tǒng)包含20條支路,有3臺變壓器、5臺發(fā)電機(jī)、1個(gè)無功補(bǔ)償點(diǎn),其中節(jié)點(diǎn)1為平衡節(jié)點(diǎn),節(jié)點(diǎn)2、3、6、8為PV節(jié)點(diǎn),系統(tǒng)數(shù)據(jù)詳見文獻(xiàn)[15]。為了更加清晰地對比不同算法的仿真效果,IEEE14節(jié)點(diǎn)系統(tǒng)只考慮有功網(wǎng)損和電壓穩(wěn)定裕度兩個(gè)目標(biāo)。運(yùn)用CSO算法、DE-CSO算法對該系統(tǒng)進(jìn)行無功優(yōu)化,仿真效果對比如圖2所示。

    圖2 IEEE14節(jié)點(diǎn)系統(tǒng)中兩種算法的Pareto前沿分布

    對比CSO算法、DE-CSO算法仿真的曲線可以看出,經(jīng)過DE-CSO優(yōu)化后的Pareto前沿分布更加均勻、廣泛,其收斂精度更高。這是因?yàn)镈E-CSO算法融合了DE算法和CSO算法各自的優(yōu)點(diǎn),它一方面彌補(bǔ)了DE算法在優(yōu)化過程中容易陷入局部最優(yōu)的缺點(diǎn),通過CSO算法縱橫交叉體制產(chǎn)生的新種群能夠有效跳出局部最優(yōu),提高收斂精度;另一方面融合了DE算法中的交叉操作和變異操作,改變了CSO算法收斂速度慢、調(diào)整策略單一的缺點(diǎn),有效地提高了種群的多樣性。所以,運(yùn)用DE-CSO算法優(yōu)化后的多目標(biāo)無功優(yōu)化Pareto前沿分布明顯優(yōu)于傳統(tǒng)的CSO算法,數(shù)據(jù)對比如表1所示。

    表1 IEEE14節(jié)點(diǎn)系統(tǒng)結(jié)果對比

    3.2 IEEE57節(jié)點(diǎn)算例

    該系統(tǒng)包含41條支路,有15臺變壓器、7臺發(fā)電機(jī)、3個(gè)無功補(bǔ)償點(diǎn),其中節(jié)點(diǎn)1為平衡節(jié)點(diǎn),節(jié)點(diǎn)2、3、6、8、9、12為PV節(jié)點(diǎn)。分別運(yùn)用CSO算法、DE-CSO算法對該系統(tǒng)進(jìn)行無功優(yōu)化,系統(tǒng)數(shù)據(jù)詳見文獻(xiàn)[15]。仿真效果如圖3、4所示。

    圖3 IEEE57節(jié)點(diǎn)CSO算法Pareto前沿分布

    對比CSO算法和DE-CSO算法優(yōu)化后的三維仿真效果圖可以看出,DE-CSO算法產(chǎn)生的Pareto前沿分布均勻、廣泛,同時(shí)收斂精度更高,其優(yōu)化效果明顯優(yōu)于傳統(tǒng)的CSO算法,充分地顯示了混合算法良好的收斂能力及尋優(yōu)效率,達(dá)到了預(yù)期的效果,數(shù)據(jù)對比如表2所示。

    圖4 IEEE57節(jié)點(diǎn) DE-CSO算法Pareto前沿分布

    表2 IEEE57節(jié)點(diǎn)系統(tǒng)結(jié)果對比

    4 結(jié) 語

    為了解決電力系統(tǒng)無功優(yōu)化多約束、多變量的問題,創(chuàng)新性地提出了一種融合DE算法和CSO算法的多目標(biāo)DE-CSO算法。新算法改進(jìn)了標(biāo)準(zhǔn)DE算法和CSO算法,提高了算法的收斂效率,能有效地防止粒子陷入局部最優(yōu)。同時(shí),在采用非支配排序法的基礎(chǔ)上,通過DE-CSO算法的運(yùn)算增加了粒子的多樣性。為提高Pareto分布的均勻性,通過計(jì)算粒子間的擁擠距離來控制目標(biāo)函數(shù)的解集。在IEEE14和IEEE57節(jié)點(diǎn)的系統(tǒng)中進(jìn)行仿真,產(chǎn)生Pareto多目標(biāo)分布前沿。仿真結(jié)果表明,新算法能有效地降低電力系統(tǒng)有功網(wǎng)損和電壓偏移量、提高靜態(tài)電壓穩(wěn)定裕度,為解決電力系統(tǒng)無功優(yōu)化問題提供了新的思路。

    猜你喜歡
    父代交叉種群
    邢氏水蕨成功繁衍并建立種群 等
    農(nóng)村家庭父代在家庭現(xiàn)代性轉(zhuǎn)型中的作用研究
    中國高等教育的代際傳遞及其內(nèi)在機(jī)制:“學(xué)二代”現(xiàn)象存在嗎?
    延遲退休決策對居民家庭代際收入流動性的影響分析
    ——基于人力資本傳遞機(jī)制
    山西省發(fā)現(xiàn)刺五加種群分布
    “六法”巧解分式方程
    男孩偏好激勵(lì)父代掙取更多收入了嗎?
    ——基于子女?dāng)?shù)量基本確定的情形
    連一連
    基于Fast-ICA的Wigner-Ville分布交叉項(xiàng)消除方法
    雙線性時(shí)頻分布交叉項(xiàng)提取及損傷識別應(yīng)用
    中文字幕免费在线视频6| 国产一区二区在线观看日韩| www.色视频.com| 伦理电影免费视频| 午夜福利影视在线免费观看| 91精品国产九色| 国产精品麻豆人妻色哟哟久久| 一个人看视频在线观看www免费| 成人亚洲精品一区在线观看 | 高清日韩中文字幕在线| 嫩草影院新地址| 丝袜脚勾引网站| 亚洲图色成人| 少妇高潮的动态图| 亚洲精品日韩av片在线观看| 亚洲精品乱久久久久久| 国产成人免费无遮挡视频| 国产精品99久久99久久久不卡 | 成年人午夜在线观看视频| 国产欧美日韩一区二区三区在线 | 如何舔出高潮| 97超碰精品成人国产| 九九在线视频观看精品| www.色视频.com| 校园人妻丝袜中文字幕| 干丝袜人妻中文字幕| 成人综合一区亚洲| 大话2 男鬼变身卡| 精品午夜福利在线看| 午夜福利视频精品| 91精品国产九色| 人人妻人人看人人澡| 欧美xxⅹ黑人| 国产成人精品婷婷| 精品久久久精品久久久| 新久久久久国产一级毛片| 在线精品无人区一区二区三 | 晚上一个人看的免费电影| 国产av一区二区精品久久 | 成人特级av手机在线观看| 精品久久久精品久久久| 久久精品国产亚洲av天美| 亚洲精品一二三| 亚洲第一区二区三区不卡| 99热国产这里只有精品6| 亚洲国产精品一区三区| 韩国高清视频一区二区三区| 尾随美女入室| 亚洲精品日韩在线中文字幕| 亚洲精品国产成人久久av| 亚洲四区av| av在线app专区| 亚洲,欧美,日韩| 精品国产一区二区三区久久久樱花 | 成年人午夜在线观看视频| 久久久久久久久久久丰满| 亚洲国产日韩一区二区| 最黄视频免费看| 成人漫画全彩无遮挡| 亚洲精品日本国产第一区| av在线app专区| 国产成人精品福利久久| 精品国产三级普通话版| 欧美xxxx黑人xx丫x性爽| 国产v大片淫在线免费观看| 久久久成人免费电影| 高清av免费在线| 亚洲欧美清纯卡通| 免费少妇av软件| 精品国产一区二区三区久久久樱花 | 伦理电影免费视频| 建设人人有责人人尽责人人享有的 | 欧美日韩综合久久久久久| 日本色播在线视频| 女人十人毛片免费观看3o分钟| 国产视频内射| 看免费成人av毛片| 国产伦精品一区二区三区视频9| 激情 狠狠 欧美| 狂野欧美激情性bbbbbb| 欧美高清成人免费视频www| 纵有疾风起免费观看全集完整版| 欧美少妇被猛烈插入视频| 卡戴珊不雅视频在线播放| 国产亚洲一区二区精品| 麻豆成人av视频| 亚洲精品自拍成人| 青青草视频在线视频观看| 久久精品夜色国产| 亚洲国产精品999| 日产精品乱码卡一卡2卡三| 亚洲欧洲日产国产| 麻豆成人午夜福利视频| 国产在线男女| av在线观看视频网站免费| 男女国产视频网站| 久久精品久久久久久久性| 欧美老熟妇乱子伦牲交| 亚洲熟女精品中文字幕| 在线观看免费高清a一片| 亚洲精品456在线播放app| 亚洲性久久影院| 日韩一本色道免费dvd| 色婷婷av一区二区三区视频| 日韩欧美精品免费久久| 日本欧美视频一区| 天天躁夜夜躁狠狠久久av| 久久久国产一区二区| 美女cb高潮喷水在线观看| 国产亚洲av片在线观看秒播厂| 免费大片黄手机在线观看| 成年美女黄网站色视频大全免费 | 日本免费在线观看一区| 一本色道久久久久久精品综合| 亚洲综合精品二区| 成人综合一区亚洲| freevideosex欧美| 99九九线精品视频在线观看视频| 久久久久久久大尺度免费视频| 99re6热这里在线精品视频| 日产精品乱码卡一卡2卡三| 国产精品爽爽va在线观看网站| 久久精品夜色国产| 色5月婷婷丁香| 中文欧美无线码| 一级毛片久久久久久久久女| 国产91av在线免费观看| 美女脱内裤让男人舔精品视频| 亚洲国产精品专区欧美| 王馨瑶露胸无遮挡在线观看| 国产精品麻豆人妻色哟哟久久| 国产成人精品久久久久久| 五月玫瑰六月丁香| 天堂俺去俺来也www色官网| 亚洲性久久影院| 亚洲精品日韩av片在线观看| 午夜老司机福利剧场| 免费黄频网站在线观看国产| 日韩av在线免费看完整版不卡| 噜噜噜噜噜久久久久久91| 色视频www国产| 精品少妇黑人巨大在线播放| 精品一区二区免费观看| 在线 av 中文字幕| 欧美激情国产日韩精品一区| 亚洲av男天堂| 狂野欧美白嫩少妇大欣赏| 亚洲无线观看免费| 亚洲怡红院男人天堂| 性色avwww在线观看| 免费少妇av软件| 91久久精品电影网| 大片电影免费在线观看免费| 秋霞在线观看毛片| 看十八女毛片水多多多| 2018国产大陆天天弄谢| 亚洲av综合色区一区| 免费看光身美女| 日韩精品有码人妻一区| 国产片特级美女逼逼视频| 亚洲av中文字字幕乱码综合| 中文欧美无线码| 亚洲色图综合在线观看| 成人美女网站在线观看视频| 国产午夜精品久久久久久一区二区三区| 免费大片18禁| 在线观看人妻少妇| 两个人的视频大全免费| 免费观看性生交大片5| 老熟女久久久| 在线观看av片永久免费下载| 久久久久人妻精品一区果冻| 亚洲图色成人| 22中文网久久字幕| 亚洲欧美日韩卡通动漫| 亚洲国产毛片av蜜桃av| 国产精品成人在线| 国模一区二区三区四区视频| 小蜜桃在线观看免费完整版高清| 菩萨蛮人人尽说江南好唐韦庄| 一本久久精品| 欧美高清成人免费视频www| 精华霜和精华液先用哪个| 欧美xxxx黑人xx丫x性爽| 久久国产精品男人的天堂亚洲 | 国产乱来视频区| 日韩国内少妇激情av| 日本wwww免费看| 美女cb高潮喷水在线观看| 人妻 亚洲 视频| 一区二区三区精品91| 亚洲欧美日韩卡通动漫| 秋霞伦理黄片| 亚洲国产av新网站| 狂野欧美激情性bbbbbb| av福利片在线观看| 18禁在线播放成人免费| 亚洲国产成人一精品久久久| 亚洲精品aⅴ在线观看| 久久精品国产自在天天线| 女人十人毛片免费观看3o分钟| 老女人水多毛片| 亚洲精品亚洲一区二区| 国语对白做爰xxxⅹ性视频网站| 精品一品国产午夜福利视频| 亚洲av成人精品一二三区| 99久久精品国产国产毛片| 80岁老熟妇乱子伦牲交| 我的女老师完整版在线观看| 久久精品久久久久久久性| 一个人看的www免费观看视频| 日韩三级伦理在线观看| 免费看日本二区| 久久国产精品大桥未久av | 国产淫片久久久久久久久| 国产淫语在线视频| 26uuu在线亚洲综合色| 春色校园在线视频观看| 亚洲欧美日韩东京热| 高清日韩中文字幕在线| 亚洲成人一二三区av| 全区人妻精品视频| 亚洲精品久久午夜乱码| 亚洲三级黄色毛片| 日韩不卡一区二区三区视频在线| 秋霞在线观看毛片| 亚洲电影在线观看av| 久久ye,这里只有精品| 熟妇人妻不卡中文字幕| 男女啪啪激烈高潮av片| 插阴视频在线观看视频| 在线播放无遮挡| 久久久成人免费电影| 最近2019中文字幕mv第一页| 久久久久久久久久人人人人人人| a级毛片免费高清观看在线播放| 特大巨黑吊av在线直播| 国产爱豆传媒在线观看| 国产免费一区二区三区四区乱码| 91狼人影院| 一本—道久久a久久精品蜜桃钙片| 亚洲精品成人av观看孕妇| 女人久久www免费人成看片| 国产亚洲91精品色在线| 99久久人妻综合| 中文天堂在线官网| 最后的刺客免费高清国语| 亚洲欧美成人精品一区二区| 国产伦精品一区二区三区视频9| 欧美精品亚洲一区二区| 久久人妻熟女aⅴ| 18禁在线播放成人免费| 成年女人在线观看亚洲视频| 99热网站在线观看| 国产v大片淫在线免费观看| 日本-黄色视频高清免费观看| 丰满迷人的少妇在线观看| 国产乱人偷精品视频| 日韩av在线免费看完整版不卡| 欧美成人a在线观看| 国产在线男女| 五月玫瑰六月丁香| 又爽又黄a免费视频| 亚洲一区二区三区欧美精品| 美女cb高潮喷水在线观看| 亚洲,欧美,日韩| 亚洲国产欧美在线一区| 久久精品国产a三级三级三级| 久久久亚洲精品成人影院| 久久韩国三级中文字幕| 91狼人影院| 亚洲国产欧美在线一区| 国产乱人偷精品视频| 日产精品乱码卡一卡2卡三| 夜夜骑夜夜射夜夜干| 日日摸夜夜添夜夜添av毛片| 国产一区亚洲一区在线观看| 亚洲图色成人| 深夜a级毛片| 亚洲av免费高清在线观看| 国产成人一区二区在线| 黑人猛操日本美女一级片| av黄色大香蕉| av国产免费在线观看| 水蜜桃什么品种好| 午夜福利在线观看免费完整高清在| 天堂8中文在线网| 免费少妇av软件| 在线观看一区二区三区| 啦啦啦啦在线视频资源| 国产成人精品久久久久久| 嫩草影院新地址| 日韩欧美一区视频在线观看 | 久久99热这里只频精品6学生| 精品一区二区三卡| 日韩大片免费观看网站| 男人狂女人下面高潮的视频| 久久99精品国语久久久| a 毛片基地| 波野结衣二区三区在线| 日本av手机在线免费观看| 欧美97在线视频| 九九久久精品国产亚洲av麻豆| 天堂俺去俺来也www色官网| 大话2 男鬼变身卡| 亚洲va在线va天堂va国产| 日本黄色日本黄色录像| 丰满迷人的少妇在线观看| 欧美一级a爱片免费观看看| 看非洲黑人一级黄片| 51国产日韩欧美| 国产av码专区亚洲av| 色视频在线一区二区三区| 亚洲丝袜综合中文字幕| 精品午夜福利在线看| 蜜桃久久精品国产亚洲av| 汤姆久久久久久久影院中文字幕| 亚洲国产精品成人久久小说| 欧美日韩一区二区视频在线观看视频在线| av卡一久久| 日日摸夜夜添夜夜添av毛片| 精品一区在线观看国产| 秋霞伦理黄片| 日韩三级伦理在线观看| 精品亚洲成国产av| 国产精品秋霞免费鲁丝片| 超碰97精品在线观看| 天天躁夜夜躁狠狠久久av| 最近手机中文字幕大全| 黑丝袜美女国产一区| 午夜免费鲁丝| 色5月婷婷丁香| 成人特级av手机在线观看| 亚洲人成网站高清观看| 色婷婷av一区二区三区视频| 久久久久久久久久久丰满| 哪个播放器可以免费观看大片| 国产人妻一区二区三区在| 欧美日韩综合久久久久久| 日日摸夜夜添夜夜爱| 亚洲精品一二三| 少妇的逼好多水| 国产精品麻豆人妻色哟哟久久| 少妇人妻 视频| 成人二区视频| 街头女战士在线观看网站| 亚洲国产精品成人久久小说| 我的老师免费观看完整版| 欧美成人一区二区免费高清观看| 国产白丝娇喘喷水9色精品| 免费大片18禁| 麻豆成人av视频| 男女国产视频网站| 免费看不卡的av| 大片电影免费在线观看免费| 成人综合一区亚洲| 又粗又硬又长又爽又黄的视频| 涩涩av久久男人的天堂| 色吧在线观看| 日日啪夜夜爽| 少妇人妻久久综合中文| 欧美成人午夜免费资源| 搡女人真爽免费视频火全软件| 国产v大片淫在线免费观看| 汤姆久久久久久久影院中文字幕| 亚洲精品成人av观看孕妇| 成人毛片a级毛片在线播放| 亚洲综合精品二区| 亚洲国产高清在线一区二区三| 18禁在线无遮挡免费观看视频| 日韩国内少妇激情av| 免费大片18禁| 伊人久久精品亚洲午夜| 在线观看美女被高潮喷水网站| 国产精品不卡视频一区二区| 欧美一区二区亚洲| 精品少妇黑人巨大在线播放| 国产av国产精品国产| 国产精品99久久久久久久久| 日本av手机在线免费观看| 高清视频免费观看一区二区| 青青草视频在线视频观看| 色5月婷婷丁香| 国产精品一区二区在线不卡| 久久av网站| 久久影院123| 欧美高清成人免费视频www| 一级毛片电影观看| 久久久欧美国产精品| 国产亚洲午夜精品一区二区久久| 日韩视频在线欧美| 麻豆国产97在线/欧美| 丝袜喷水一区| 日韩中文字幕视频在线看片 | 欧美另类一区| 亚洲精华国产精华液的使用体验| 搡老乐熟女国产| 午夜免费男女啪啪视频观看| 激情 狠狠 欧美| 亚洲精华国产精华液的使用体验| 亚洲av电影在线观看一区二区三区| 午夜视频国产福利| 狂野欧美激情性xxxx在线观看| av视频免费观看在线观看| 七月丁香在线播放| 国产精品麻豆人妻色哟哟久久| 国产精品一及| 直男gayav资源| 日本vs欧美在线观看视频 | 久久久久久久久久久免费av| 欧美成人一区二区免费高清观看| 夜夜看夜夜爽夜夜摸| 最近手机中文字幕大全| 在线看a的网站| 国内精品宾馆在线| av网站免费在线观看视频| 精品一区二区三区视频在线| 国内少妇人妻偷人精品xxx网站| 国产淫片久久久久久久久| 高清黄色对白视频在线免费看 | 免费少妇av软件| 大码成人一级视频| 久久久精品免费免费高清| 有码 亚洲区| 国产淫语在线视频| av一本久久久久| av视频免费观看在线观看| 97在线视频观看| 亚洲精品一二三| 18禁在线播放成人免费| 黄色配什么色好看| 美女内射精品一级片tv| 亚洲国产高清在线一区二区三| 夜夜骑夜夜射夜夜干| 制服丝袜香蕉在线| 乱系列少妇在线播放| 亚洲精品乱码久久久久久按摩| 免费在线观看成人毛片| 亚洲欧美日韩无卡精品| 国产精品伦人一区二区| 视频区图区小说| 高清在线视频一区二区三区| 最近最新中文字幕免费大全7| 日韩一区二区三区影片| av网站免费在线观看视频| 欧美精品亚洲一区二区| 亚洲精品久久久久久婷婷小说| 青春草国产在线视频| 久久国产乱子免费精品| 直男gayav资源| 一本色道久久久久久精品综合| 亚洲精品国产av成人精品| 久久国产精品大桥未久av | 极品少妇高潮喷水抽搐| 啦啦啦在线观看免费高清www| 日韩欧美一区视频在线观看 | 人妻系列 视频| 精品人妻视频免费看| 如何舔出高潮| 精品一区二区免费观看| 久久国产精品男人的天堂亚洲 | 亚洲,欧美,日韩| 91在线精品国自产拍蜜月| h日本视频在线播放| 插阴视频在线观看视频| 国产高潮美女av| 少妇 在线观看| 欧美日韩在线观看h| 纯流量卡能插随身wifi吗| 久久99热这里只有精品18| 又黄又爽又刺激的免费视频.| 日韩av免费高清视频| 欧美老熟妇乱子伦牲交| 国精品久久久久久国模美| .国产精品久久| 99久国产av精品国产电影| 久久久久久九九精品二区国产| 国产视频内射| 视频区图区小说| 国产色爽女视频免费观看| 一边亲一边摸免费视频| 日日摸夜夜添夜夜添av毛片| 一二三四中文在线观看免费高清| 永久网站在线| 久久久久性生活片| 少妇被粗大猛烈的视频| 乱码一卡2卡4卡精品| 一级片'在线观看视频| 人人妻人人看人人澡| 久久青草综合色| 精品国产一区二区三区久久久樱花 | 亚洲人成网站在线播| 久久久久久伊人网av| 国产淫语在线视频| 视频中文字幕在线观看| 国产精品一区二区在线观看99| 男人和女人高潮做爰伦理| 毛片一级片免费看久久久久| 夫妻性生交免费视频一级片| 国产精品嫩草影院av在线观看| 高清黄色对白视频在线免费看 | 人妻 亚洲 视频| 欧美日韩视频高清一区二区三区二| 五月开心婷婷网| 午夜福利影视在线免费观看| 亚洲精品乱码久久久久久按摩| 国产精品蜜桃在线观看| 成人国产麻豆网| 国产日韩欧美在线精品| 一级二级三级毛片免费看| 久久久久国产网址| 99久国产av精品国产电影| 2022亚洲国产成人精品| 国产深夜福利视频在线观看| 亚洲av福利一区| 亚洲精品一二三| 六月丁香七月| 日韩欧美 国产精品| 青春草国产在线视频| 亚洲av成人精品一二三区| 夜夜爽夜夜爽视频| 青青草视频在线视频观看| 老司机影院毛片| 黄色视频在线播放观看不卡| 欧美日韩视频精品一区| 中文字幕人妻熟人妻熟丝袜美| 一级a做视频免费观看| 日韩国内少妇激情av| 最近手机中文字幕大全| 日日啪夜夜撸| 99久久人妻综合| 免费久久久久久久精品成人欧美视频 | 免费观看在线日韩| 少妇熟女欧美另类| 久久久成人免费电影| 欧美变态另类bdsm刘玥| 亚洲,一卡二卡三卡| 高清黄色对白视频在线免费看 | 人人妻人人爽人人添夜夜欢视频 | 简卡轻食公司| 成人黄色视频免费在线看| 国产高清三级在线| 国产亚洲最大av| 久久久久久久久久久免费av| 乱系列少妇在线播放| 最黄视频免费看| 男人添女人高潮全过程视频| av国产精品久久久久影院| 亚洲无线观看免费| 一本久久精品| 黄片无遮挡物在线观看| 高清在线视频一区二区三区| 国产高清国产精品国产三级 | 简卡轻食公司| 精品久久久久久久久亚洲| 一级爰片在线观看| 少妇裸体淫交视频免费看高清| 一级毛片黄色毛片免费观看视频| 人妻系列 视频| 久久久久久伊人网av| 午夜免费男女啪啪视频观看| 美女内射精品一级片tv| 日韩一本色道免费dvd| 久久久久久人妻| 夫妻午夜视频| 国产白丝娇喘喷水9色精品| 国产精品久久久久久精品电影小说 | 在线 av 中文字幕| 观看av在线不卡| 国产免费又黄又爽又色| 免费观看性生交大片5| 两个人的视频大全免费| 国产欧美日韩一区二区三区在线 | 国产精品秋霞免费鲁丝片| 成人午夜精彩视频在线观看| 天天躁日日操中文字幕| 22中文网久久字幕| 国产久久久一区二区三区| 国国产精品蜜臀av免费| 亚洲经典国产精华液单| 久久久久性生活片| 亚洲一级一片aⅴ在线观看| 国产伦在线观看视频一区| 久久国内精品自在自线图片| 哪个播放器可以免费观看大片| 日韩电影二区| 亚洲欧美日韩卡通动漫| 久久久亚洲精品成人影院| 欧美成人午夜免费资源| 舔av片在线| 久久久久国产网址| 狂野欧美白嫩少妇大欣赏| 九九久久精品国产亚洲av麻豆| 亚洲欧美成人精品一区二区| 精品久久久精品久久久| 亚洲熟女精品中文字幕| 久久精品熟女亚洲av麻豆精品| 一二三四中文在线观看免费高清| 男女啪啪激烈高潮av片| av黄色大香蕉| 综合色丁香网| 国产免费一级a男人的天堂| 中文字幕亚洲精品专区| 黑人猛操日本美女一级片| 大陆偷拍与自拍| 建设人人有责人人尽责人人享有的 | 久久午夜福利片| 2022亚洲国产成人精品| tube8黄色片| 亚洲aⅴ乱码一区二区在线播放| 日韩av不卡免费在线播放| 日日啪夜夜爽| 人人妻人人添人人爽欧美一区卜 | 国模一区二区三区四区视频| 日本与韩国留学比较|