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

    一種適用于機(jī)組組合優(yōu)化的改進(jìn)整數(shù)編碼粒子群算法

    2017-06-12 01:55:57吳和海熊高峰袁晉蓉秦躍杰
    現(xiàn)代電子技術(shù) 2017年11期
    關(guān)鍵詞:粒子群算法

    吳和海++熊高峰+袁晉蓉+秦躍杰

    摘 要: 針對機(jī)組組合這一高維、非線性混合整數(shù)規(guī)劃問題,提出一種結(jié)合修補(bǔ)策略的整數(shù)編碼粒子群(ICPSO) 算法。用正負(fù)整數(shù)分別表示機(jī)組開停機(jī)的時(shí)間長度,有效減少待優(yōu)化變量個(gè)數(shù)?;跈C(jī)組組合問題的特點(diǎn),采用修補(bǔ)策略處理不滿足約束條件的個(gè)體,使算法只在可行解區(qū)域內(nèi)搜索,有效提高收斂速度,通過切除冗余機(jī)組,提高解的質(zhì)量。仿真算例表明,相比整型編碼遺傳(r?ICGA)算法、改進(jìn)粒子群(IPSO)算法、社會(huì)演化(SEP)算法,提出的ICPSO算法能夠更有效地處理大規(guī)模機(jī)組組合優(yōu)化問題,執(zhí)行時(shí)間較短、求解精度更高。

    關(guān)鍵詞: 機(jī)組組合; 粒子群算法; 整數(shù)編碼; 修補(bǔ)策略

    中圖分類號: TN919?34; TP18 文獻(xiàn)標(biāo)識碼: A 文章編號: 1004?373X(2017)11?0167?05

    An modified integer?coded particle swarm optimization algorithm

    for unit commitment optimization

    WU Hehai, XIONG Gaofeng, YUAN Jinrong, QIN Yuejie

    (College of Electrical and Information Engineering, Hunan University, Changsha 410082, China)

    Abstract: Aiming at the high dimension and nonlinear mixed integer programming problems of unit commitment, an integer?coded particle swarm optimization (ICPSO) algorithm combined with repairing scheme is proposed. The positive and negative integers are used to represent the time length of unit startup and shutdown to reduce the quantity of variables under optimization effectively. On the basis of the characteristic of unit commitment, the repairing scheme is adopted to handle the individuals which can′t conform to the constraint conditions, so as to make it only search in the feasible solution region, and improve the convergence rate. The redundancy unit is cut out to improve the quality of solution. The simulation example results show that, in comparison with r?ICGA, IPSO algorithm and SEP algorithm, the ICPSO algorithm can handle the combinatorial optimization of large?scale units effectively, has shorter execution time and higher solving accuracy.

    Keywords: unit commitment; particle swarm optimization; integer coding; repairing scheme

    0 引 言

    機(jī)組組合優(yōu)化是電力系統(tǒng)經(jīng)濟(jì)調(diào)度中一個(gè)非常重要的問題,它的高維數(shù)、非凸、離散、非線性,使得理論上很難求取其最優(yōu)解。但由于它能夠帶來顯著的經(jīng)濟(jì)效益,許多國內(nèi)外學(xué)者一直在積極研究,提出各種方法來解決這個(gè)問題,例如,動(dòng)態(tài)規(guī)劃法、拉格朗日松弛法、優(yōu)先順序法等確定方法以及遺傳算法、蟻群算法、模擬退火等智能優(yōu)化算法[1?2]。但到目前為止,以上方法均存在一定的缺陷,如面臨“維數(shù)災(zāi)”、求解時(shí)間較長、精度不夠高等,不適合實(shí)際機(jī)組組合的優(yōu)化。

    粒子群優(yōu)化算法(PSO)是一種隨機(jī)全局優(yōu)化算法,該算法自提出以來,因?yàn)槠渚哂泻唵?、收斂速度快且效果顯著等優(yōu)點(diǎn),在各個(gè)領(lǐng)域得到廣泛的應(yīng)用[3?6]。目前粒子群算法在求解機(jī)組組合問題時(shí)[7?8],普遍采用二進(jìn)制編碼方式,用0?1表示發(fā)電機(jī)組的開停機(jī)狀態(tài);處理約束條件在適應(yīng)度函數(shù)中加入懲罰項(xiàng)。對于大規(guī)模機(jī)組系統(tǒng),這種編碼方式會(huì)導(dǎo)致種群個(gè)體長度過長,影響算法搜索效率,浪費(fèi)大量求解時(shí)間[9];懲罰項(xiàng)的引入會(huì)使得計(jì)算時(shí)間的極大增加,同時(shí)懲罰乘子的調(diào)整也存在困難[10]。

    針對上述問題,本文提出一種適用于大規(guī)模機(jī)組組合問題的整數(shù)編碼粒子群算法。采用整數(shù)編碼方式來表達(dá)機(jī)組組合問題,有效減少待優(yōu)化變量個(gè)數(shù);同時(shí)采用粒子修補(bǔ)策略處理不滿足約束條件的粒子,以提高搜索速度和解的質(zhì)量。算例仿真驗(yàn)證了文中所提出的算法能夠有效求解電力系統(tǒng)大規(guī)模機(jī)組組合問題。

    1 機(jī)組組合的數(shù)學(xué)模型

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

    電力系統(tǒng)機(jī)組組合問題是指在規(guī)定的時(shí)期內(nèi),在滿足功率平衡等約束條件下,合理的安排機(jī)組的開停機(jī)順序與出力,使得系統(tǒng)的總發(fā)電成本最小。機(jī)組組合問題的目標(biāo)函數(shù)可表示為[11]:

    (1)

    式中:為系統(tǒng)總發(fā)電成本;代表機(jī)組臺(tái)數(shù);為時(shí)段總個(gè)數(shù);為機(jī)組在時(shí)段內(nèi)的輸出功率;為機(jī)組在時(shí)段內(nèi)的開停機(jī)狀態(tài),當(dāng)時(shí)表示機(jī)組開機(jī),時(shí)表示機(jī)組停機(jī);表示機(jī)組在時(shí)段內(nèi)的耗量成本,其中,為機(jī)組的運(yùn)行成本參數(shù);表示機(jī)組在時(shí)段內(nèi)的啟動(dòng)成本,它與機(jī)組啟動(dòng)之前的連續(xù)停機(jī)時(shí)間長短有關(guān)。

    1.2 約束條件

    (1) 功率平衡約束:

    (2)

    (2) 最小開/停機(jī)時(shí)間約束:

    (3)

    (4)

    (3) 機(jī)組出力上下限約束:

    (5)

    (4) 旋轉(zhuǎn)備用約束:

    (6)

    (5) 機(jī)組爬坡約束:

    (7)

    (8)

    (6) 最大開關(guān)機(jī)次數(shù)約束:

    (9)

    式中:表示調(diào)度時(shí)段內(nèi)的系統(tǒng)負(fù)荷需求;表示機(jī)組開機(jī)后的最小開機(jī)時(shí)間;表示機(jī)組的最小停運(yùn)時(shí)間;表示機(jī)組的最小出力;表示機(jī)組的最大出力;為時(shí)段內(nèi)的旋轉(zhuǎn)備用值;分別為機(jī)組的功率最大上調(diào)量、下調(diào)量;表示機(jī)組在調(diào)度周期內(nèi)的開停機(jī)次數(shù);為機(jī)組在調(diào)度周期內(nèi)允許的最大開停機(jī)次數(shù)。

    2 基于整數(shù)編碼粒子群算法的機(jī)組組合問題

    求解

    2.1 整數(shù)編碼

    用正負(fù)整數(shù)分別表示機(jī)組開停機(jī)時(shí)間長度,機(jī)組開停機(jī)狀態(tài)設(shè)計(jì)為矩陣: (10)

    式中:為粒子種群中的第個(gè)個(gè)體;表示機(jī)組的狀態(tài)階段所持續(xù)的時(shí)間,正整數(shù)值代表機(jī)組連續(xù)開機(jī)的小時(shí)數(shù),負(fù)整數(shù)值代表機(jī)組持續(xù)停機(jī)的小時(shí)數(shù);表示機(jī)組開停機(jī)階段數(shù)的上限(每臺(tái)機(jī)組一般每天開停機(jī)次數(shù)不超過5次),一些重負(fù)荷機(jī)組的開停機(jī)階段數(shù)不足時(shí)用數(shù)字0填充,從而保持粒子矩陣列數(shù)恒定。每行表示某機(jī)組在調(diào)度周期內(nèi)的開停機(jī)安排,滿足。第臺(tái)機(jī)組所有時(shí)段的啟停安排由一組正負(fù)相間的整數(shù)串來表示,構(gòu)成矩陣的行向量。所有機(jī)組的整數(shù)串(行向量)組成整個(gè)調(diào)度周期全部機(jī)組的啟停安排,即機(jī)組組合問題的解。

    圖1詳細(xì)演示了如何通過一個(gè)10×5的整數(shù)編碼矩陣表達(dá)10機(jī)組在24 h內(nèi)的一個(gè)調(diào)度安排。若使用二進(jìn)制編碼,同樣的調(diào)度所需的矩陣規(guī)模增大至10×24,計(jì)算量將大大增加。例如第六行向量 [-8,6,-5,4,-1]表示此機(jī)組在1~8時(shí)段內(nèi)停機(jī),9~14時(shí)段內(nèi)開機(jī),15~19時(shí)段內(nèi)停機(jī),20~23時(shí)段內(nèi)開機(jī),24時(shí)段內(nèi)停機(jī)。該編碼方式從特性上就限制了機(jī)組在一個(gè)調(diào)度周期內(nèi)的啟停狀態(tài)次數(shù)最多不超過滿足機(jī)組最大開停機(jī)次數(shù)的約束。

    2.2 粒子種群的初始化

    (1) 時(shí),機(jī)組在第一個(gè)狀態(tài)階段持續(xù)的時(shí)間:

    (11)

    式中:表示機(jī)組在上一個(gè)調(diào)度周期最后一個(gè)狀態(tài)階段持續(xù)的時(shí)間,這樣生成的能夠保證機(jī)組延續(xù)前一個(gè)調(diào)度周期的狀態(tài)且滿足最小開停機(jī)時(shí)間約束。

    (2) 時(shí),考慮到機(jī)組需滿足最小開停機(jī)時(shí)間約束,可按式(12)初始化:

    (12)

    式中:表示第個(gè)狀態(tài)階段后剩余的調(diào)度時(shí)段數(shù)。

    (13)

    (3) 即最后一個(gè)狀態(tài)階段,它所持續(xù)的時(shí)間初始化如下:

    (14)

    由于隨機(jī)數(shù)的原因,若前個(gè)狀態(tài)階段的持續(xù)時(shí)間之和等于調(diào)度周期時(shí)段數(shù)則余下的狀態(tài)階段所持續(xù)的時(shí)間用0來補(bǔ)充,保證所有的粒子矩陣有相同的列數(shù)。根據(jù)式(11)~式(14)初始化得到的粒子種群,滿足機(jī)組最小開停機(jī)時(shí)間約束。

    2.3 粒子進(jìn)化

    粒子在尋優(yōu)過程中,依據(jù)迄今為止找到的最優(yōu)解和整個(gè)粒子群搜索到的最優(yōu)解來調(diào)整粒子的進(jìn)化速度和方向。

    粒子的速度按照下式進(jìn)行變化:

    (15)

    粒子位置的更新公式如下:

    (16)

    式中:為種群粒子個(gè)數(shù);為第次迭代后粒子的速度;為第次迭代后粒子的位置;為粒子個(gè)體的歷史最優(yōu)位置;為種群最好位置;rand1,rand2為[0,1]之間的隨機(jī)數(shù);是加速度常數(shù);為慣性權(quán)重;Int為取整符號。

    2.4 修補(bǔ)粒子

    在粒子群迭代優(yōu)化過程中不可避免地會(huì)產(chǎn)生一些無效的粒子,在這些無效粒子中會(huì)有各個(gè)狀態(tài)階段的持續(xù)時(shí)間之和不等于調(diào)度周期時(shí)段數(shù)整數(shù)串不是正負(fù)相間的情況,從而導(dǎo)致該粒子無效。需要對粒子個(gè)體進(jìn)行修補(bǔ),按行向量的先后順序進(jìn)行。以個(gè)體中行為例,具體措施如下:

    (1) 若保持前個(gè)狀態(tài)階段持續(xù)時(shí)間不變,其中滿足且令段的值變?yōu)?,正負(fù)保持不變;

    (2) 若,考慮到此種情況修補(bǔ)繁瑣、耗時(shí)較多,則重新初始化該粒子;

    (3) 若行不是正負(fù)相間的整數(shù)串,則保持的大小不變,依次調(diào)整它們的正負(fù)號以滿足條件;

    (4) 若有的情況,則交換兩者的位置,否則不滿足編碼方式的特性。

    2.5 約束條件處理

    2.5.1 旋轉(zhuǎn)備用約束處理

    每個(gè)時(shí)段都應(yīng)該對旋轉(zhuǎn)備用約束進(jìn)行檢驗(yàn),若違反約束,則應(yīng)調(diào)整機(jī)組開停機(jī)狀態(tài)。本文采用機(jī)組最小耗量與機(jī)組的最大出力之比作為指標(biāo),對機(jī)組進(jìn)行優(yōu)先排序[12]。當(dāng)時(shí)段所有開機(jī)機(jī)組的最大輸出功率之和小于系統(tǒng)負(fù)荷量和旋轉(zhuǎn)備用時(shí),對按照指標(biāo)排列好的機(jī)組逐一檢查其開停機(jī)狀態(tài),將狀態(tài)為0的機(jī)組開啟,直到滿足旋轉(zhuǎn)備用約束,轉(zhuǎn)至?xí)r段繼續(xù)檢驗(yàn)。

    2.5.2 最小開停機(jī)時(shí)間約束處理

    違背最小開機(jī)時(shí)間約束通常出現(xiàn)在高峰負(fù)荷時(shí)段,因?yàn)榉逯禃r(shí)間往往低于最小開機(jī)時(shí)間。類似地,因低負(fù)荷持續(xù)時(shí)間一般低于最短關(guān)機(jī)時(shí)間,最小關(guān)機(jī)約束的違背常出現(xiàn)于低負(fù)荷時(shí)段。在粒子群進(jìn)化過程中,采用最小開停機(jī)時(shí)間修補(bǔ)策略對粒子編碼矩陣進(jìn)行修復(fù)[13]。

    如圖2所示,若檢查是否滿足最小開機(jī)時(shí)間約束,若不滿足,則將變?yōu)?,機(jī)組繼續(xù)開機(jī);如圖3所示,若檢查~時(shí)刻機(jī)組是否都能滿足最小停機(jī)時(shí)間約束,若不滿足,則機(jī)組繼續(xù)開機(jī),變?yōu)?。

    2.6 切除冗余機(jī)組

    在求解過程中,由于缺少對機(jī)組運(yùn)行實(shí)際特點(diǎn)的考慮,易出現(xiàn)小容量機(jī)組運(yùn)行過多、機(jī)組備用容量過剩的情況,使得發(fā)電成本增加。因此有必要關(guān)閉多余備用機(jī)組,即切除冗余機(jī)組。

    在調(diào)度周期內(nèi),依次檢查每個(gè)時(shí)段內(nèi)所有開機(jī)機(jī)組的最大出力之和是否超過系統(tǒng)負(fù)荷和旋轉(zhuǎn)備用。若超過,則按優(yōu)先表逆序關(guān)閉已開機(jī)組。再重新檢測,直至該時(shí)段機(jī)組輸出總功率不過量。

    3 算例分析

    本文采用經(jīng)典算例對算法的性能進(jìn)行檢驗(yàn),為了便于對比分析,不考慮爬坡約束,可以直接采用等微增率準(zhǔn)則進(jìn)行負(fù)荷經(jīng)濟(jì)分配。計(jì)算中粒子種群大小取20;最大迭代進(jìn)化次數(shù)取100;均取2;權(quán)重系數(shù)從0.9~0.4線性遞減,粒子最大速度為10。計(jì)算周期為一天,全天分為24個(gè)時(shí)段。機(jī)組開停機(jī)階段數(shù)的上限考慮到粒子群算法的隨機(jī)性,應(yīng)用本文ICPSO算法進(jìn)行20次獨(dú)立迭代計(jì)算,取其中最好的解作為最終優(yōu)化結(jié)果。

    采用文獻(xiàn)[8?9]中的10臺(tái)機(jī)組參數(shù)和24 h的負(fù)荷數(shù)據(jù)以及旋轉(zhuǎn)備用,并根據(jù)上述文獻(xiàn)中的復(fù)制方法將系統(tǒng)規(guī)模擴(kuò)大至20臺(tái)、40臺(tái)、60臺(tái)、80臺(tái)和100臺(tái),系統(tǒng)負(fù)荷和備用同時(shí)按比例增加。

    表1給出ICPSO算法求得的10機(jī)組系統(tǒng)最優(yōu)運(yùn)行成本,系統(tǒng)總的運(yùn)行成本為$563 937.7,其中啟動(dòng)總成本費(fèi)用為$4 090,機(jī)組發(fā)電總成本為$559 847.7,最優(yōu)調(diào)度表見表2。

    表3為本文的ICPSO算法與其他優(yōu)化算法[8?9,11]對10~100臺(tái)機(jī)組系統(tǒng)求解結(jié)果的比較??梢钥闯觯瑢τ?0機(jī)組系統(tǒng),本文ICPSO算法和r?ICGA都能獲得相同的總費(fèi)用,優(yōu)于SEP,IPSO算法;隨著機(jī)組臺(tái)數(shù)的增加,ICPSO的求解結(jié)果均優(yōu)于其他三種方法,驗(yàn)證了本文算法求解大規(guī)模機(jī)組組合問題的有效性和高精度性。

    圖4給出了本文算法在不同機(jī)組規(guī)模下的執(zhí)行時(shí)間曲線,隨著機(jī)組臺(tái)數(shù)的增加,執(zhí)行時(shí)間近似線性增長,10機(jī)組時(shí)耗時(shí)約3.2 s,100機(jī)組時(shí)也只需要28.2 s左右。這是由于本文采用修復(fù)策略對不可行粒子進(jìn)行修復(fù),使算法只在可行解區(qū)域內(nèi)搜索最優(yōu)解,有效提高收斂速度;同時(shí)相比二進(jìn)制編碼,本文的整數(shù)編碼方式能有效提高粒子群算法的搜索效率。因此本文提出的整數(shù)編碼粒子群算法更適用于求解大規(guī)模機(jī)組組合問題。

    圖4 不同機(jī)組規(guī)模的執(zhí)行時(shí)間曲線

    4 結(jié) 語

    本文采用整數(shù)編碼粒子群算法來求解機(jī)組組合問題。通過整數(shù)編碼方式,用正負(fù)整數(shù)分別表示機(jī)組開停機(jī)時(shí)間長度,相比于二進(jìn)制編碼它能有效減少待優(yōu)化變量個(gè)數(shù)。求解過程中針對機(jī)組組合問題的特點(diǎn),采用修補(bǔ)策略處理不滿足約束條件的粒子,使算法只在可行解區(qū)域內(nèi)搜索,通過切除冗余機(jī)組提高解的精度。仿真算例表明,相比r?ICGA,IPSO,SEP,本文提出的ICPSO算法求解大規(guī)模機(jī)組組合優(yōu)化問題具有更高的精度,求解時(shí)間大幅度減少。

    參考文獻(xiàn)

    [1] KAZARLIS S A, BAKIRTZIS A G, PETRIDIS V. A genetic algorithm solution to the unit commitment problem [J]. IEEE transactions on power systems, 1996, 11(1): 83?92.

    [2] 黎靜華,蘭飛.機(jī)組組合問題的模型及算法綜述[J].現(xiàn)代電力,2011,28(6):1?10.

    [3] 李丹,高立群,王珂,等.電力系統(tǒng)機(jī)組組合問題的動(dòng)態(tài)雙種群粒子群算法[J].計(jì)算機(jī)應(yīng)用,2008,28(1):104?107.

    [4] 張濤,史蘇怡,徐雪琴.基于二進(jìn)制量子粒子群算法的含分布式電源配電網(wǎng)重構(gòu)[J].電力系統(tǒng)保護(hù)與控制,2016,44(4):22?28.

    [5] 方源,章桐,陳霏霏,等.電動(dòng)車動(dòng)力總成噪聲品質(zhì)粒子群?向量機(jī)預(yù)測模型[J].西安交通大學(xué)學(xué)報(bào),2016,50(1):41?46.

    [6] 王春枝,張會(huì)麗,葉志偉,等.一種基于混沌粒子群算法和支持向量機(jī)的 P2P流量識別方法[J].計(jì)算機(jī)應(yīng)用與軟件,2015(8):288?291.

    [7] 耿曉軍,王剛.粒子群算法應(yīng)用于機(jī)組組合問題的優(yōu)化[J].現(xiàn)代電子技術(shù),2007,30(16):111?113.

    [8] 趙波,曹一家.電力系統(tǒng)機(jī)組組合問題的改進(jìn)粒子群優(yōu)化算法[J].電網(wǎng)技術(shù),2004,28(21):6?10.

    [9] 張偉,趙進(jìn)慧,王寧.帶修復(fù)操作整型編碼遺傳算法求解大規(guī)模機(jī)組組合問題[J].化工學(xué)報(bào),2012,63(9):2972?2979.

    [10] AMJADY N, SHIRZADI A. Unit commitment using a new integer coded genetic algorithm [J]. European transactions on electrical power, 2009, 19(8): 1161?1176.

    [11] 王喆,余貽鑫,張弘鵬.社會(huì)演化算法在機(jī)組組合中的應(yīng)用[J].中國電機(jī)工程學(xué)報(bào),2004,24(4):12?17.

    [12] 黎靜華,韋化.求解機(jī)組組合問題的領(lǐng)域搜索法[J].中國電機(jī)工程學(xué)報(bào),2008,28(13):33?40.

    [13] 蔡振華.考慮電價(jià)的機(jī)組組合問題研究[D].長沙:湖南大學(xué),2013.

    猜你喜歡
    粒子群算法
    幾種改進(jìn)的螢火蟲算法性能比較及應(yīng)用
    基于支持向量機(jī)的短期電力負(fù)荷預(yù)測
    基于云計(jì)算平臺(tái)的資源調(diào)度優(yōu)化研究
    一種基于高維粒子群算法的神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)優(yōu)化研究
    基于PSODE混合算法優(yōu)化的自抗擾控制器設(shè)計(jì)
    蟻群算法的運(yùn)用及其優(yōu)化分析
    電力市場交易背景下水電站優(yōu)化調(diào)度研究
    基于粒子群算法的產(chǎn)業(yè)技術(shù)創(chuàng)新生態(tài)系統(tǒng)運(yùn)行穩(wěn)定性組合評價(jià)研究
    預(yù)測(2016年5期)2016-12-26 10:04:59
    無線傳感器網(wǎng)絡(luò)聯(lián)盟初始結(jié)構(gòu)生成研究
    交通堵塞擾動(dòng)下多車場車輛路徑優(yōu)化
    商(2016年5期)2016-03-28 18:10:26
    国产黄a三级三级三级人| 夜夜躁狠狠躁天天躁| 国产伦一二天堂av在线观看| 又粗又爽又猛毛片免费看| 欧美黑人欧美精品刺激| 亚洲五月婷婷丁香| 日本一本二区三区精品| 久久久久久免费高清国产稀缺| 啦啦啦免费观看视频1| 久久久精品欧美日韩精品| 午夜精品在线福利| 最近最新中文字幕大全电影3| 国产精品国产高清国产av| 村上凉子中文字幕在线| 男女之事视频高清在线观看| 久久久久国产一级毛片高清牌| 人人妻人人看人人澡| 我的老师免费观看完整版| 丰满的人妻完整版| 禁无遮挡网站| 久久久久久久久免费视频了| 国产视频一区二区在线看| 在线观看一区二区三区| 中国美女看黄片| 一卡2卡三卡四卡精品乱码亚洲| 亚洲专区字幕在线| 国产一区在线观看成人免费| 亚洲一码二码三码区别大吗| 亚洲人成77777在线视频| 亚洲熟妇熟女久久| 中文字幕精品亚洲无线码一区| 亚洲 欧美一区二区三区| 一边摸一边做爽爽视频免费| 久久国产精品影院| 日本熟妇午夜| 亚洲av第一区精品v没综合| 制服丝袜大香蕉在线| 啪啪无遮挡十八禁网站| 日本一二三区视频观看| 黄片大片在线免费观看| 日韩高清综合在线| 亚洲五月婷婷丁香| 成人av在线播放网站| 美女扒开内裤让男人捅视频| 国产精品 欧美亚洲| 国产v大片淫在线免费观看| 亚洲av电影不卡..在线观看| 99热这里只有精品一区 | 亚洲无线在线观看| 久久久久久国产a免费观看| 99热只有精品国产| 国产精品九九99| 国产伦人伦偷精品视频| 国产精品久久久久久人妻精品电影| www国产在线视频色| 成人永久免费在线观看视频| 久久人人精品亚洲av| 叶爱在线成人免费视频播放| 亚洲国产中文字幕在线视频| 久久香蕉精品热| 91麻豆精品激情在线观看国产| 亚洲精品久久国产高清桃花| 夜夜躁狠狠躁天天躁| 日本黄大片高清| 日本精品一区二区三区蜜桃| 亚洲人与动物交配视频| 日本熟妇午夜| 国产精品久久久人人做人人爽| 午夜福利在线观看吧| www.熟女人妻精品国产| 身体一侧抽搐| 淫秽高清视频在线观看| 国产v大片淫在线免费观看| 757午夜福利合集在线观看| 啦啦啦免费观看视频1| 欧美zozozo另类| 午夜福利成人在线免费观看| 国语自产精品视频在线第100页| 久久精品aⅴ一区二区三区四区| 啪啪无遮挡十八禁网站| 观看免费一级毛片| 伊人久久大香线蕉亚洲五| 九九热线精品视视频播放| 免费在线观看亚洲国产| 午夜激情av网站| 国产精品av视频在线免费观看| 亚洲五月天丁香| 亚洲欧洲精品一区二区精品久久久| 精品不卡国产一区二区三区| 麻豆国产97在线/欧美 | 久久婷婷成人综合色麻豆| 久久中文看片网| 亚洲五月婷婷丁香| 国产亚洲精品av在线| 亚洲乱码一区二区免费版| 国产伦人伦偷精品视频| 成人18禁高潮啪啪吃奶动态图| 动漫黄色视频在线观看| 中文字幕熟女人妻在线| 色av中文字幕| 午夜两性在线视频| 亚洲狠狠婷婷综合久久图片| 99热这里只有精品一区 | 国产爱豆传媒在线观看 | 一二三四在线观看免费中文在| 丰满的人妻完整版| 一本一本综合久久| 国产伦一二天堂av在线观看| 国产精品,欧美在线| 成年人黄色毛片网站| 日日夜夜操网爽| 国产三级中文精品| 男女之事视频高清在线观看| 午夜精品一区二区三区免费看| 老司机靠b影院| 亚洲国产欧美人成| 色av中文字幕| 老汉色∧v一级毛片| 国产精品日韩av在线免费观看| 亚洲成a人片在线一区二区| 亚洲自偷自拍图片 自拍| 99re在线观看精品视频| 好男人在线观看高清免费视频| 看片在线看免费视频| 精品人妻1区二区| 18禁美女被吸乳视频| 99热这里只有精品一区 | 午夜福利在线在线| 一个人免费在线观看电影 | 国语自产精品视频在线第100页| 免费av毛片视频| 亚洲欧美日韩高清专用| 欧美绝顶高潮抽搐喷水| 岛国视频午夜一区免费看| 亚洲人成电影免费在线| 久久久久国产精品人妻aⅴ院| 两个人免费观看高清视频| 国产精华一区二区三区| 一本综合久久免费| 少妇被粗大的猛进出69影院| 热99re8久久精品国产| 亚洲一区二区三区不卡视频| 欧美黄色片欧美黄色片| 国产成人aa在线观看| 亚洲一区二区三区不卡视频| 美女高潮喷水抽搐中文字幕| 久久精品国产99精品国产亚洲性色| 国语自产精品视频在线第100页| 国产成年人精品一区二区| 黄色片一级片一级黄色片| 一本大道久久a久久精品| 黄色女人牲交| 黄色片一级片一级黄色片| 亚洲va日本ⅴa欧美va伊人久久| 婷婷六月久久综合丁香| 热99re8久久精品国产| 免费在线观看视频国产中文字幕亚洲| 18禁美女被吸乳视频| 成人三级做爰电影| 日韩免费av在线播放| 日本a在线网址| 伦理电影免费视频| 国产男靠女视频免费网站| 大型av网站在线播放| 国产又黄又爽又无遮挡在线| 无遮挡黄片免费观看| 中亚洲国语对白在线视频| 99精品在免费线老司机午夜| 精品人妻1区二区| 久久亚洲真实| 亚洲欧美日韩高清在线视频| av在线播放免费不卡| 禁无遮挡网站| 少妇熟女aⅴ在线视频| 亚洲,欧美精品.| 三级毛片av免费| 亚洲欧美日韩东京热| 久久中文看片网| 精品一区二区三区四区五区乱码| 激情在线观看视频在线高清| 亚洲国产欧美人成| 亚洲天堂国产精品一区在线| 美女免费视频网站| 可以在线观看的亚洲视频| 久久久久国产精品人妻aⅴ院| 亚洲国产欧洲综合997久久,| 午夜福利成人在线免费观看| 日本免费一区二区三区高清不卡| 18美女黄网站色大片免费观看| 国产激情欧美一区二区| 在线视频色国产色| 在线观看日韩欧美| 亚洲精品一卡2卡三卡4卡5卡| 国产av不卡久久| 国产精品亚洲av一区麻豆| 全区人妻精品视频| 最近最新中文字幕大全免费视频| 男男h啪啪无遮挡| 一区福利在线观看| 免费在线观看亚洲国产| 亚洲一区二区三区色噜噜| 欧美日韩亚洲综合一区二区三区_| 亚洲第一电影网av| 午夜福利18| 男插女下体视频免费在线播放| 亚洲av片天天在线观看| 五月玫瑰六月丁香| 国产精品永久免费网站| 亚洲自拍偷在线| 十八禁人妻一区二区| 91在线观看av| 日本精品一区二区三区蜜桃| 欧美国产日韩亚洲一区| 男女午夜视频在线观看| 免费在线观看黄色视频的| 手机成人av网站| 床上黄色一级片| 熟女少妇亚洲综合色aaa.| 少妇粗大呻吟视频| 俺也久久电影网| 亚洲一区二区三区不卡视频| 国产精品九九99| 久久伊人香网站| 男人舔奶头视频| 精品久久久久久久末码| 老司机深夜福利视频在线观看| 99精品欧美一区二区三区四区| 少妇的丰满在线观看| 人妻夜夜爽99麻豆av| 亚洲国产精品成人综合色| 国产精品野战在线观看| 国产精品免费一区二区三区在线| 欧美黄色淫秽网站| 日韩精品青青久久久久久| 亚洲五月天丁香| 欧美日本视频| 亚洲中文av在线| 欧洲精品卡2卡3卡4卡5卡区| 国产精品久久久av美女十八| 很黄的视频免费| 女同久久另类99精品国产91| 狠狠狠狠99中文字幕| 日日夜夜操网爽| 国产av一区二区精品久久| 99热6这里只有精品| 黑人巨大精品欧美一区二区mp4| 叶爱在线成人免费视频播放| 脱女人内裤的视频| 国产私拍福利视频在线观看| 熟女少妇亚洲综合色aaa.| 此物有八面人人有两片| 香蕉国产在线看| 少妇熟女aⅴ在线视频| 久久精品综合一区二区三区| 99热这里只有是精品50| 国产精品一及| 老司机午夜福利在线观看视频| 久久久久久九九精品二区国产 | 日韩欧美 国产精品| 夜夜爽天天搞| www日本在线高清视频| 99在线人妻在线中文字幕| 亚洲avbb在线观看| 一级a爱片免费观看的视频| 亚洲人成电影免费在线| 老汉色∧v一级毛片| 母亲3免费完整高清在线观看| 黑人操中国人逼视频| 久久久久久久久久黄片| 国产高清videossex| 亚洲自拍偷在线| 精品熟女少妇八av免费久了| 国产精品久久久久久久电影 | or卡值多少钱| 精品第一国产精品| 啦啦啦免费观看视频1| 非洲黑人性xxxx精品又粗又长| 亚洲熟女毛片儿| 亚洲一码二码三码区别大吗| 欧美性长视频在线观看| 日本免费a在线| 99国产精品99久久久久| 国产av不卡久久| 国产人伦9x9x在线观看| 97人妻精品一区二区三区麻豆| bbb黄色大片| 国产午夜福利久久久久久| 啪啪无遮挡十八禁网站| 婷婷亚洲欧美| av在线天堂中文字幕| 夜夜躁狠狠躁天天躁| 少妇被粗大的猛进出69影院| 日本在线视频免费播放| 久久这里只有精品中国| 成人亚洲精品av一区二区| 国产男靠女视频免费网站| 精品熟女少妇八av免费久了| 国产高清视频在线播放一区| 最近在线观看免费完整版| 在线a可以看的网站| 51午夜福利影视在线观看| 国产激情偷乱视频一区二区| 深夜精品福利| 最新美女视频免费是黄的| 日本 av在线| 精品久久蜜臀av无| 国产精品久久久av美女十八| 少妇的丰满在线观看| 99国产极品粉嫩在线观看| 国产欧美日韩一区二区精品| 免费在线观看影片大全网站| 久久久国产成人免费| 亚洲熟女毛片儿| 中亚洲国语对白在线视频| www日本在线高清视频| 欧美日韩亚洲综合一区二区三区_| 一级毛片高清免费大全| 九九热线精品视视频播放| 精品国内亚洲2022精品成人| 亚洲,欧美精品.| 国产av一区在线观看免费| 美女午夜性视频免费| 国产69精品久久久久777片 | 我要搜黄色片| 99久久精品国产亚洲精品| ponron亚洲| 给我免费播放毛片高清在线观看| 久久精品91蜜桃| 亚洲精华国产精华精| 岛国在线观看网站| 99久久无色码亚洲精品果冻| av免费在线观看网站| 成人午夜高清在线视频| 又黄又粗又硬又大视频| 黄色视频,在线免费观看| 亚洲中文日韩欧美视频| 亚洲av美国av| 最新在线观看一区二区三区| 激情在线观看视频在线高清| 男女下面进入的视频免费午夜| 国产精华一区二区三区| 成人欧美大片| 精品久久久久久,| 夜夜爽天天搞| 中文字幕久久专区| 97人妻精品一区二区三区麻豆| 大型av网站在线播放| 久久草成人影院| 亚洲第一电影网av| 精品电影一区二区在线| 欧美一区二区国产精品久久精品 | 亚洲国产日韩欧美精品在线观看 | 淫妇啪啪啪对白视频| 91字幕亚洲| 亚洲九九香蕉| 久久久久国产一级毛片高清牌| 在线观看免费视频日本深夜| 露出奶头的视频| 天堂av国产一区二区熟女人妻 | 制服人妻中文乱码| 久久久久久久久久黄片| 免费在线观看日本一区| 变态另类丝袜制服| 男男h啪啪无遮挡| 亚洲无线在线观看| 久久久久免费精品人妻一区二区| 一个人免费在线观看电影 | 亚洲成人中文字幕在线播放| 成熟少妇高潮喷水视频| 婷婷六月久久综合丁香| 国产又黄又爽又无遮挡在线| 一二三四社区在线视频社区8| 亚洲av成人精品一区久久| 深夜精品福利| 免费在线观看日本一区| 久久久久久大精品| 搡老熟女国产l中国老女人| or卡值多少钱| 欧美乱码精品一区二区三区| 一夜夜www| 国产精品乱码一区二三区的特点| 亚洲一码二码三码区别大吗| 色老头精品视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 在线观看免费午夜福利视频| 久久精品人妻少妇| 国产97色在线日韩免费| 18禁黄网站禁片免费观看直播| 欧美一级毛片孕妇| 日韩精品免费视频一区二区三区| 在线播放国产精品三级| 在线免费观看的www视频| 国产aⅴ精品一区二区三区波| 免费电影在线观看免费观看| 日韩大尺度精品在线看网址| 九九热线精品视视频播放| 亚洲午夜理论影院| 国产伦人伦偷精品视频| 777久久人妻少妇嫩草av网站| 午夜福利在线观看吧| 久久香蕉精品热| 欧美成狂野欧美在线观看| 欧美日韩黄片免| 国产真实乱freesex| 日本免费一区二区三区高清不卡| 国产亚洲欧美98| 欧美av亚洲av综合av国产av| 国产激情久久老熟女| 在线永久观看黄色视频| 黄色a级毛片大全视频| 亚洲精品粉嫩美女一区| 香蕉久久夜色| 成人午夜高清在线视频| 一边摸一边做爽爽视频免费| 成人av在线播放网站| 国产精品久久久久久亚洲av鲁大| 亚洲第一欧美日韩一区二区三区| 国产真人三级小视频在线观看| 午夜成年电影在线免费观看| 18禁国产床啪视频网站| 可以在线观看毛片的网站| 国产精品爽爽va在线观看网站| 国产成+人综合+亚洲专区| 最近最新中文字幕大全免费视频| 国产亚洲av高清不卡| tocl精华| 国产精品影院久久| 成人三级做爰电影| 精品乱码久久久久久99久播| 可以在线观看毛片的网站| 中文字幕熟女人妻在线| 成人国产综合亚洲| 黄片大片在线免费观看| 日韩免费av在线播放| e午夜精品久久久久久久| 搡老岳熟女国产| √禁漫天堂资源中文www| 亚洲,欧美精品.| 日日爽夜夜爽网站| 久久中文字幕人妻熟女| 曰老女人黄片| 男插女下体视频免费在线播放| 少妇粗大呻吟视频| 久久久久久亚洲精品国产蜜桃av| 国产精品香港三级国产av潘金莲| 久久热在线av| 露出奶头的视频| 欧美性长视频在线观看| a在线观看视频网站| 老司机靠b影院| 久久久久性生活片| 亚洲片人在线观看| 亚洲欧美日韩高清在线视频| 听说在线观看完整版免费高清| 在线十欧美十亚洲十日本专区| 午夜福利在线观看吧| 老汉色∧v一级毛片| 免费在线观看亚洲国产| 日韩精品免费视频一区二区三区| 国产爱豆传媒在线观看 | 18禁黄网站禁片免费观看直播| 国内久久婷婷六月综合欲色啪| 久久久久久久久久黄片| 亚洲午夜理论影院| 99久久无色码亚洲精品果冻| 久久婷婷成人综合色麻豆| 99久久国产精品久久久| 国产一区二区在线观看日韩 | 亚洲精品久久国产高清桃花| 日本免费一区二区三区高清不卡| x7x7x7水蜜桃| 人人妻,人人澡人人爽秒播| 99久久国产精品久久久| 日韩免费av在线播放| 操出白浆在线播放| 国产一区二区三区视频了| 国产高清有码在线观看视频 | 一区二区三区激情视频| 国产高清视频在线播放一区| 18美女黄网站色大片免费观看| 看黄色毛片网站| 嫩草影院精品99| 夜夜躁狠狠躁天天躁| 草草在线视频免费看| 午夜久久久久精精品| 久久精品国产亚洲av香蕉五月| www.精华液| 久久久国产欧美日韩av| 老司机午夜十八禁免费视频| 午夜视频精品福利| 别揉我奶头~嗯~啊~动态视频| 一二三四在线观看免费中文在| 在线观看午夜福利视频| 精品少妇一区二区三区视频日本电影| 高清毛片免费观看视频网站| 三级国产精品欧美在线观看 | 亚洲av熟女| 亚洲中文字幕日韩| 久久久久久久精品吃奶| 亚洲欧美日韩东京热| 男女之事视频高清在线观看| 国产成+人综合+亚洲专区| av在线天堂中文字幕| 人人妻人人看人人澡| 黄色视频不卡| 啦啦啦观看免费观看视频高清| 深夜精品福利| 亚洲在线自拍视频| 日韩欧美三级三区| 99热这里只有是精品50| 88av欧美| 人妻丰满熟妇av一区二区三区| 免费高清视频大片| 国产在线精品亚洲第一网站| 欧美黑人欧美精品刺激| 黄色a级毛片大全视频| 亚洲熟妇中文字幕五十中出| 免费观看精品视频网站| 亚洲精品美女久久久久99蜜臀| 免费观看精品视频网站| 久久久水蜜桃国产精品网| 日日夜夜操网爽| 亚洲中文av在线| 91九色精品人成在线观看| 1024视频免费在线观看| 色哟哟哟哟哟哟| 国产成人一区二区三区免费视频网站| 五月伊人婷婷丁香| 最近最新中文字幕大全电影3| 亚洲国产欧美人成| 免费在线观看视频国产中文字幕亚洲| 又大又爽又粗| 久久精品人妻少妇| 欧美3d第一页| 亚洲激情在线av| 亚洲在线自拍视频| 黄色视频,在线免费观看| 无限看片的www在线观看| 最近最新免费中文字幕在线| 天天躁夜夜躁狠狠躁躁| 特大巨黑吊av在线直播| 久久99热这里只有精品18| 黄色a级毛片大全视频| 亚洲狠狠婷婷综合久久图片| 精品福利观看| 给我免费播放毛片高清在线观看| 午夜影院日韩av| 国产精品永久免费网站| 国产精品九九99| 国产精品乱码一区二三区的特点| 最近最新免费中文字幕在线| 欧美一区二区精品小视频在线| 少妇的丰满在线观看| 日本 av在线| 看免费av毛片| 国产精品精品国产色婷婷| 精品福利观看| 18禁黄网站禁片午夜丰满| av欧美777| 午夜精品一区二区三区免费看| 舔av片在线| 国语自产精品视频在线第100页| 麻豆国产av国片精品| av中文乱码字幕在线| 日韩欧美在线二视频| 久久精品综合一区二区三区| 欧美日韩国产亚洲二区| 欧美黑人巨大hd| 好男人电影高清在线观看| 国产精品亚洲av一区麻豆| 91麻豆精品激情在线观看国产| 一级毛片高清免费大全| 亚洲黑人精品在线| 每晚都被弄得嗷嗷叫到高潮| 欧美3d第一页| 在线免费观看的www视频| 亚洲国产日韩欧美精品在线观看 | 悠悠久久av| 欧美3d第一页| 国产av一区二区精品久久| 国内精品久久久久精免费| 丝袜美腿诱惑在线| 女同久久另类99精品国产91| 久久这里只有精品19| 两个人免费观看高清视频| 国产1区2区3区精品| 搡老妇女老女人老熟妇| 在线观看免费视频日本深夜| 午夜激情福利司机影院| 两个人免费观看高清视频| 99国产精品99久久久久| 夜夜爽天天搞| 亚洲全国av大片| 正在播放国产对白刺激| 一个人免费在线观看电影 | 久久久国产精品麻豆| 午夜两性在线视频| 91麻豆av在线| 久久精品综合一区二区三区| 亚洲熟女毛片儿| 欧美成人免费av一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 精品不卡国产一区二区三区| 成年女人毛片免费观看观看9| 国产亚洲欧美在线一区二区| 亚洲成人免费电影在线观看| 欧美成人免费av一区二区三区| 91成年电影在线观看| 亚洲电影在线观看av| 精品久久久久久成人av| 天天添夜夜摸| 我的老师免费观看完整版| 成人18禁在线播放| 黄频高清免费视频| 久久精品人妻少妇|