• <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)用
    精品国产国语对白av| www.精华液| 一本色道久久久久久精品综合| 中文精品一卡2卡3卡4更新| 国产精品一区二区在线不卡| 精品国产乱子伦一区二区三区 | 亚洲av美国av| 美女脱内裤让男人舔精品视频| 久久免费观看电影| 制服诱惑二区| 国产成人一区二区三区免费视频网站| 夜夜夜夜夜久久久久| 我的亚洲天堂| 色播在线永久视频| 亚洲熟女毛片儿| 午夜福利视频在线观看免费| 一本一本久久a久久精品综合妖精| 亚洲,欧美精品.| 18禁国产床啪视频网站| 免费高清在线观看日韩| 啪啪无遮挡十八禁网站| 在线十欧美十亚洲十日本专区| a 毛片基地| 亚洲欧美激情在线| av线在线观看网站| 少妇被粗大的猛进出69影院| 久久中文字幕一级| 男女下面插进去视频免费观看| 天天操日日干夜夜撸| 两性夫妻黄色片| 欧美人与性动交α欧美精品济南到| 最新在线观看一区二区三区| 欧美大码av| 真人做人爱边吃奶动态| 一边摸一边做爽爽视频免费| 日韩欧美国产一区二区入口| 亚洲国产毛片av蜜桃av| 精品人妻在线不人妻| 老汉色av国产亚洲站长工具| 大香蕉久久网| 天天操日日干夜夜撸| 亚洲美女黄色视频免费看| 中文字幕人妻熟女乱码| 久久久久视频综合| a级毛片黄视频| 美女中出高潮动态图| 肉色欧美久久久久久久蜜桃| 国产精品欧美亚洲77777| 在线十欧美十亚洲十日本专区| av在线老鸭窝| 亚洲精品美女久久久久99蜜臀| 精品国产一区二区三区久久久樱花| 狠狠狠狠99中文字幕| 精品亚洲成国产av| 999久久久精品免费观看国产| 亚洲第一av免费看| 亚洲精品美女久久久久99蜜臀| 亚洲一码二码三码区别大吗| 一区二区三区精品91| 日韩制服骚丝袜av| www.精华液| 日日摸夜夜添夜夜添小说| 老司机影院成人| 日韩欧美国产一区二区入口| 丰满人妻熟妇乱又伦精品不卡| 日日夜夜操网爽| 欧美国产精品va在线观看不卡| 亚洲三区欧美一区| 午夜精品国产一区二区电影| 午夜免费成人在线视频| 在线精品无人区一区二区三| 久久天躁狠狠躁夜夜2o2o| 爱豆传媒免费全集在线观看| 咕卡用的链子| 亚洲av成人不卡在线观看播放网 | 99久久国产精品久久久| 国产精品一区二区免费欧美 | 欧美在线一区亚洲| 极品人妻少妇av视频| 国产免费视频播放在线视频| 热re99久久精品国产66热6| 国产欧美日韩一区二区三 | 中文字幕色久视频| 熟女少妇亚洲综合色aaa.| e午夜精品久久久久久久| 国产精品 欧美亚洲| 久久人人爽av亚洲精品天堂| 色婷婷久久久亚洲欧美| 秋霞在线观看毛片| 91成人精品电影| 亚洲av成人不卡在线观看播放网 | 中文字幕人妻丝袜制服| 日韩视频一区二区在线观看| 亚洲七黄色美女视频| 纵有疾风起免费观看全集完整版| 免费女性裸体啪啪无遮挡网站| 精品免费久久久久久久清纯 | 亚洲成国产人片在线观看| 亚洲第一欧美日韩一区二区三区 | 欧美日韩亚洲综合一区二区三区_| 欧美日韩亚洲综合一区二区三区_| 男女免费视频国产| 最近最新免费中文字幕在线| 十八禁人妻一区二区| 亚洲色图 男人天堂 中文字幕| 精品久久蜜臀av无| xxxhd国产人妻xxx| 在线观看舔阴道视频| 手机成人av网站| 中文字幕制服av| 在线永久观看黄色视频| 国产成人精品在线电影| 国产黄色免费在线视频| 18禁裸乳无遮挡动漫免费视频| 黄色a级毛片大全视频| 岛国在线观看网站| 亚洲全国av大片| 热99久久久久精品小说推荐| 亚洲av电影在线进入| 亚洲人成电影观看| 日韩大片免费观看网站| 日日夜夜操网爽| 国产亚洲一区二区精品| 精品一区二区三卡| 精品国产乱码久久久久久男人| 国产成人啪精品午夜网站| 久久人人爽av亚洲精品天堂| 亚洲欧美激情在线| 欧美在线黄色| 国产亚洲欧美在线一区二区| 十八禁人妻一区二区| 大码成人一级视频| 精品国内亚洲2022精品成人 | 欧美精品高潮呻吟av久久| 俄罗斯特黄特色一大片| 精品久久久久久电影网| 中文字幕最新亚洲高清| 色婷婷久久久亚洲欧美| 久久久国产欧美日韩av| 日韩视频一区二区在线观看| 性少妇av在线| 亚洲中文日韩欧美视频| 91字幕亚洲| 视频区欧美日本亚洲| 麻豆国产av国片精品| 国产一区二区 视频在线| 欧美成人午夜精品| 国产免费av片在线观看野外av| 老鸭窝网址在线观看| 狂野欧美激情性xxxx| 黑人操中国人逼视频| 亚洲国产精品成人久久小说| 亚洲成国产人片在线观看| 一区二区av电影网| 国产亚洲午夜精品一区二区久久| 午夜日韩欧美国产| 中文字幕制服av| 在线永久观看黄色视频| 国产主播在线观看一区二区| 人成视频在线观看免费观看| 一级毛片电影观看| 狂野欧美激情性bbbbbb| 国产一卡二卡三卡精品| 老鸭窝网址在线观看| 成人影院久久| 黄网站色视频无遮挡免费观看| 两个人免费观看高清视频| 狠狠婷婷综合久久久久久88av| 一本—道久久a久久精品蜜桃钙片| 精品国产国语对白av| 国产亚洲av片在线观看秒播厂| 99精国产麻豆久久婷婷| 性少妇av在线| 午夜成年电影在线免费观看| 欧美 亚洲 国产 日韩一| 建设人人有责人人尽责人人享有的| 国产成人欧美在线观看 | www.精华液| 手机成人av网站| 美女高潮喷水抽搐中文字幕| 亚洲色图 男人天堂 中文字幕| xxxhd国产人妻xxx| 亚洲精品一区蜜桃| 久久人妻熟女aⅴ| 欧美中文综合在线视频| 建设人人有责人人尽责人人享有的| 亚洲成av片中文字幕在线观看| 亚洲熟女精品中文字幕| 99久久99久久久精品蜜桃| 久久久久久久国产电影| 黑丝袜美女国产一区| 亚洲欧美成人综合另类久久久| 精品高清国产在线一区| 成人国产一区最新在线观看| 国产欧美日韩一区二区三 | 亚洲久久久国产精品| 亚洲激情五月婷婷啪啪| 欧美另类一区| 午夜激情av网站| 欧美变态另类bdsm刘玥| 免费高清在线观看视频在线观看| 精品一区二区三区av网在线观看 | 国产精品一区二区精品视频观看| 男人添女人高潮全过程视频| 亚洲性夜色夜夜综合| 成人国产av品久久久| 在线av久久热| 淫妇啪啪啪对白视频 | 国产无遮挡羞羞视频在线观看| 女人高潮潮喷娇喘18禁视频| www.自偷自拍.com| 在线观看人妻少妇| 亚洲五月婷婷丁香| 国产免费现黄频在线看| 美女大奶头黄色视频| 午夜福利在线观看吧| 欧美日韩亚洲综合一区二区三区_| 人妻一区二区av| 亚洲成av片中文字幕在线观看| 人妻一区二区av| 亚洲人成77777在线视频| 在线看a的网站| 久久免费观看电影| 日韩欧美免费精品| 色精品久久人妻99蜜桃| 国产男女超爽视频在线观看| 国产不卡av网站在线观看| 悠悠久久av| 黑人巨大精品欧美一区二区mp4| 99九九在线精品视频| 久久影院123| 老司机午夜十八禁免费视频| 欧美日韩成人在线一区二区| 日韩有码中文字幕| 超色免费av| 69av精品久久久久久 | a级毛片在线看网站| 精品一区二区三区四区五区乱码| 欧美一级毛片孕妇| 日本猛色少妇xxxxx猛交久久| 久久久久精品国产欧美久久久 | 国产精品久久久久久精品电影小说| 国产亚洲av高清不卡| 欧美性长视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 欧美另类亚洲清纯唯美| 欧美中文综合在线视频| 国产区一区二久久| 日本vs欧美在线观看视频| 国产男女内射视频| 欧美日韩亚洲高清精品| 久久久久久久大尺度免费视频| 午夜两性在线视频| 日本欧美视频一区| www.999成人在线观看| 51午夜福利影视在线观看| 中文欧美无线码| 国产亚洲午夜精品一区二区久久| 91九色精品人成在线观看| 超碰97精品在线观看| 国产精品久久久人人做人人爽| 丝袜美腿诱惑在线| 一级毛片女人18水好多| 久久毛片免费看一区二区三区| 久久久精品国产亚洲av高清涩受| 日韩一区二区三区影片| 久久99一区二区三区| 日日夜夜操网爽| 久久亚洲精品不卡| 日韩大片免费观看网站| 天天躁夜夜躁狠狠躁躁| 中亚洲国语对白在线视频| 丝袜美足系列| 久久精品国产亚洲av香蕉五月 | 国产91精品成人一区二区三区 | a级片在线免费高清观看视频| 欧美黄色片欧美黄色片| 成年美女黄网站色视频大全免费| 亚洲免费av在线视频| 亚洲精品av麻豆狂野| 亚洲一区中文字幕在线| 99久久人妻综合| 十八禁网站网址无遮挡| 后天国语完整版免费观看| 亚洲欧洲日产国产| 91九色精品人成在线观看| 亚洲欧美色中文字幕在线| 欧美日韩亚洲综合一区二区三区_| 中文欧美无线码| 91大片在线观看| 一本久久精品| 激情视频va一区二区三区| 丝袜在线中文字幕| 精品视频人人做人人爽| 国产成人精品在线电影| 久久久国产成人免费| svipshipincom国产片| 国产精品久久久久成人av| av在线app专区| 黑丝袜美女国产一区| 久久久久国内视频| 久久午夜综合久久蜜桃| 一区二区日韩欧美中文字幕| 成人国语在线视频| 国产深夜福利视频在线观看| 欧美日韩福利视频一区二区| 国产成人欧美| 成人手机av| 9热在线视频观看99| 最近中文字幕2019免费版| xxxhd国产人妻xxx| 人人妻人人爽人人添夜夜欢视频| 一级毛片电影观看| 国产又爽黄色视频| 少妇的丰满在线观看| 成年av动漫网址| 国产欧美日韩一区二区精品| 国产亚洲一区二区精品| 成人手机av| 国产视频一区二区在线看| 最近最新中文字幕大全免费视频| 一本大道久久a久久精品| 久久女婷五月综合色啪小说| 99国产精品99久久久久| 黄色视频不卡| av国产精品久久久久影院| 亚洲av美国av| 9色porny在线观看| 国产伦理片在线播放av一区| 91大片在线观看| 久久毛片免费看一区二区三区| 亚洲国产中文字幕在线视频| 巨乳人妻的诱惑在线观看| 最新在线观看一区二区三区| av在线app专区| 日韩制服骚丝袜av| 老熟女久久久| 国产精品秋霞免费鲁丝片| 国产成人欧美| 99热国产这里只有精品6| 亚洲熟女毛片儿| 老鸭窝网址在线观看| 99国产精品一区二区三区| 欧美激情高清一区二区三区| 一区二区三区乱码不卡18| videosex国产| 丝袜美足系列| 丝袜美腿诱惑在线| 欧美 亚洲 国产 日韩一| 国产精品.久久久| 黄色毛片三级朝国网站| 美女高潮到喷水免费观看| 亚洲欧美日韩高清在线视频 | 满18在线观看网站| www日本在线高清视频| 国产有黄有色有爽视频| 夜夜骑夜夜射夜夜干| 天天添夜夜摸| 国产亚洲精品一区二区www | 19禁男女啪啪无遮挡网站| 淫妇啪啪啪对白视频 | 国产男女内射视频| 在线永久观看黄色视频| av一本久久久久| 久久亚洲精品不卡| 搡老乐熟女国产| 欧美变态另类bdsm刘玥| 久久ye,这里只有精品| 久久精品国产亚洲av香蕉五月 | 国产区一区二久久| 日本欧美视频一区| 多毛熟女@视频| av欧美777| 大片电影免费在线观看免费| 老司机福利观看| 天堂俺去俺来也www色官网| 亚洲av电影在线进入| 精品少妇久久久久久888优播| 午夜福利,免费看| 精品一区二区三区av网在线观看 | 成年人免费黄色播放视频| 色视频在线一区二区三区| 精品少妇黑人巨大在线播放| 脱女人内裤的视频| www.av在线官网国产| 国产免费视频播放在线视频| 成人亚洲精品一区在线观看| 亚洲专区字幕在线| 亚洲国产精品一区二区三区在线| 美女国产高潮福利片在线看| 亚洲欧美精品自产自拍| 国产精品熟女久久久久浪| 黄频高清免费视频| 一边摸一边抽搐一进一出视频| 久久久久久久久久久久大奶| 亚洲国产看品久久| 操出白浆在线播放| 亚洲午夜精品一区,二区,三区| 女人久久www免费人成看片| 日本猛色少妇xxxxx猛交久久| 老熟女久久久| 69精品国产乱码久久久| 热re99久久精品国产66热6| 天天操日日干夜夜撸| 久久热在线av| 国产野战对白在线观看| 亚洲熟女精品中文字幕| 欧美日韩精品网址| 精品国产一区二区三区久久久樱花| 9色porny在线观看| 日韩精品免费视频一区二区三区| 99久久国产精品久久久| 欧美成人午夜精品| 国产精品久久久久久人妻精品电影 | 国产欧美日韩一区二区精品| 亚洲国产成人一精品久久久| 桃花免费在线播放| 国产主播在线观看一区二区| 中文字幕最新亚洲高清| 啦啦啦 在线观看视频| 亚洲av片天天在线观看| 满18在线观看网站| 日韩中文字幕欧美一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 十八禁高潮呻吟视频| 国产成人免费观看mmmm| 欧美黑人精品巨大| 制服人妻中文乱码| 91精品三级在线观看| 国精品久久久久久国模美| 欧美日韩中文字幕国产精品一区二区三区 | 欧美大码av| 精品乱码久久久久久99久播| 国产99久久九九免费精品| 人人妻人人爽人人添夜夜欢视频| 亚洲中文日韩欧美视频| 亚洲av日韩在线播放| 一边摸一边抽搐一进一出视频| 这个男人来自地球电影免费观看| 午夜视频精品福利| 亚洲av男天堂| 成年动漫av网址| 高清在线国产一区| 欧美激情 高清一区二区三区| 一区二区三区四区激情视频| 精品久久久久久久毛片微露脸 | 在线天堂中文资源库| 国产成人免费无遮挡视频| 精品国产乱码久久久久久男人| 国产精品熟女久久久久浪| av线在线观看网站| 日本一区二区免费在线视频| 男女午夜视频在线观看| 久久久国产成人免费| 亚洲少妇的诱惑av| 成人国语在线视频| 国产欧美日韩一区二区三 | 午夜免费观看性视频| 成人影院久久| 国产在线免费精品| 啦啦啦视频在线资源免费观看| 丝袜美腿诱惑在线| 又紧又爽又黄一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜福利视频在线观看免费| 中文字幕色久视频| 亚洲视频免费观看视频| 脱女人内裤的视频| 国产av精品麻豆| 欧美精品啪啪一区二区三区 | 99国产精品一区二区三区| 亚洲一区二区三区欧美精品| 欧美亚洲日本最大视频资源| 成人国语在线视频| 亚洲成人国产一区在线观看| 啦啦啦视频在线资源免费观看| 少妇裸体淫交视频免费看高清 | 国产一区二区三区在线臀色熟女 | 亚洲色图 男人天堂 中文字幕| 人人妻,人人澡人人爽秒播| 大码成人一级视频| 亚洲第一青青草原| 中文字幕高清在线视频| 男女下面插进去视频免费观看| 国产精品99久久99久久久不卡| 国产精品成人在线| 亚洲精品一区蜜桃| 日本猛色少妇xxxxx猛交久久| 99国产精品免费福利视频| 在线观看免费高清a一片| 男女免费视频国产| 在线观看www视频免费| 日本一区二区免费在线视频| 久久精品熟女亚洲av麻豆精品| 亚洲欧美日韩高清在线视频 | 日韩中文字幕视频在线看片| 国产欧美日韩一区二区三 | 人人妻人人爽人人添夜夜欢视频| 一区二区三区乱码不卡18| 一级,二级,三级黄色视频| 久久久久久人人人人人| 午夜激情久久久久久久| 欧美激情 高清一区二区三区| 亚洲国产欧美网| 大型av网站在线播放| 精品少妇一区二区三区视频日本电影| 午夜老司机福利片| avwww免费| 色视频在线一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 丰满饥渴人妻一区二区三| 国产成人一区二区三区免费视频网站| 男女无遮挡免费网站观看| 免费在线观看日本一区| 免费观看人在逋| 久久人人爽人人片av| 亚洲国产av新网站| 在线观看免费视频网站a站| 久久久久久久精品精品| 91av网站免费观看| 不卡一级毛片| 各种免费的搞黄视频| 91麻豆精品激情在线观看国产 | 精品免费久久久久久久清纯 | 下体分泌物呈黄色| 一进一出抽搐动态| 午夜福利一区二区在线看| 美女福利国产在线| 国产成人精品久久二区二区91| 丝袜在线中文字幕| 久久久国产一区二区| 啦啦啦 在线观看视频| 男人爽女人下面视频在线观看| 久久香蕉激情| 精品一品国产午夜福利视频| av国产精品久久久久影院| 国产福利在线免费观看视频| 久久久久久人人人人人| www.精华液| 亚洲久久久国产精品| 咕卡用的链子| 国产色视频综合| 国产成人系列免费观看| 国产精品一区二区免费欧美 | 国产亚洲一区二区精品| 各种免费的搞黄视频| 在线精品无人区一区二区三| 热99国产精品久久久久久7| 日韩大码丰满熟妇| 精品亚洲乱码少妇综合久久| 国产片内射在线| 国产免费现黄频在线看| 黑人巨大精品欧美一区二区mp4| 色婷婷av一区二区三区视频| 午夜福利,免费看| 国产精品久久久久成人av| 亚洲天堂av无毛| 国产在线观看jvid| 在线天堂中文资源库| 亚洲av电影在线观看一区二区三区| 国产精品二区激情视频| 久久99热这里只频精品6学生| 欧美激情极品国产一区二区三区| 国产成人免费无遮挡视频| 99国产综合亚洲精品| 不卡av一区二区三区| 交换朋友夫妻互换小说| 一边摸一边做爽爽视频免费| tocl精华| 国产深夜福利视频在线观看| 一本色道久久久久久精品综合| 日韩视频在线欧美| 91大片在线观看| 亚洲精品美女久久av网站| 一区二区三区乱码不卡18| 久久久久久久久久久久大奶| 啦啦啦啦在线视频资源| 国产成人一区二区三区免费视频网站| 亚洲第一av免费看| 制服人妻中文乱码| 国产精品.久久久| 欧美国产精品va在线观看不卡| 久久精品成人免费网站| 婷婷成人精品国产| 少妇人妻久久综合中文| 国产淫语在线视频| 黑人欧美特级aaaaaa片| 女人久久www免费人成看片| 免费在线观看视频国产中文字幕亚洲 | 精品一品国产午夜福利视频| 亚洲欧洲日产国产| 精品人妻熟女毛片av久久网站| 一区二区三区乱码不卡18| 午夜福利在线观看吧| 国产精品久久久久久精品古装| 国产欧美日韩综合在线一区二区| 搡老乐熟女国产| www.av在线官网国产| 老司机深夜福利视频在线观看 | 免费一级毛片在线播放高清视频 | 亚洲精品美女久久av网站| 久久国产亚洲av麻豆专区| 交换朋友夫妻互换小说| 悠悠久久av| 亚洲综合色网址| av视频免费观看在线观看| 黄色毛片三级朝国网站| 大型av网站在线播放| 亚洲国产欧美日韩在线播放| 国产区一区二久久| 国产精品久久久人人做人人爽| 黄色a级毛片大全视频| 国产免费现黄频在线看| 国产精品熟女久久久久浪| 1024视频免费在线观看|