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

    基于改進蝙蝠算法的梯級水庫發(fā)電優(yōu)化調(diào)度

    2020-10-12 14:29:16張佳丹顧圣平鄭斯水唐鳳珍
    人民黃河 2020年6期

    張佳丹 顧圣平 鄭斯水 唐鳳珍

    摘?要:針對梯級水庫優(yōu)化調(diào)度高維度、多約束以及非線性優(yōu)化的特點,將自適應權重以及連續(xù)禁忌搜索算法引入標準蝙蝠算法,改善標準蝙蝠算法在水庫優(yōu)化調(diào)度應用中出現(xiàn)的早熟收斂且陷入局部最優(yōu)的問題:一方面利用自適應權重避免因更新步長機制導致尋優(yōu)能力不足的問題;另一方面利用連續(xù)禁忌搜索算法避免因種群多樣性差導致陷入局部最優(yōu)的問題。案例分析結果表明,改進蝙蝠算法能有效運用于水庫發(fā)電優(yōu)化調(diào)度中,并且與標準蝙蝠算法相比,具有更強全局尋優(yōu)能力、更高的運行效率,得到的運行調(diào)度結果更優(yōu)。

    關鍵詞:水庫優(yōu)化調(diào)度;蝙蝠算法;自適應權重;連續(xù)禁忌搜索算法

    中圖分類號:TV697.1文獻標志碼:A

    doi:10.3969/j.issn.1000-1379.2020.06.011

    Optimal Operation of Cascade Hydropower Station Reservoirs Based on the Improved Bat Algorithm

    ZHANG Jiadan, GU Shengping, ZHENG Sishui, TANG Fengzhen

    (College of Water Conservancy and Hydropower Engineering, Hohai University, Nanjing 210098, China)

    Abstract:For the optimal operation of cascade reservoirs having the characteristics of high-dimension, multi-constraint and nonlinear optimization, the adaptive weight and continuous tabu search were introduced into the basic bat algorithm to solve the issue of premature convergence and getting the local optimal result easily in the optimal operation: On the one hand, the adaptive weight was used for avoiding the lack of optimizing ability caused by the step update mechanism; on the other hand, continuous tabu search was used for avoiding getting the local optimal result caused by the decline in population diversity. The case analysis shows that the improved bat algorithm can be effectively applied to the optimal operation of reservoir power generation and compared with the basic bat algorithm, it has stronger global optimization ability, higher operation performance and better operation scheduling scheme.

    Key words: optimal operation; bat algorithm; adaptive weight; continuous tabu search

    水庫優(yōu)化調(diào)度作為一個典型的非線性優(yōu)化問題,能根據(jù)徑流資料以及綜合利用要求,運用優(yōu)化方法尋求最優(yōu)的水庫調(diào)度方案,使水庫在調(diào)度周期內(nèi)經(jīng)濟效益最高,提高水能資源利用效率。然而,決策過程的多階段性以及徑流的不確定性使得水庫優(yōu)化調(diào)度求解難度較高且求解精度較低。傳統(tǒng)的水庫優(yōu)化調(diào)度方法主要有動態(tài)規(guī)劃算法[1]、逐步尋優(yōu)算法[2]等,然而隨著水庫調(diào)度計算時段的細分,容易出現(xiàn)“維數(shù)災”、耗時長的問題。近年來,粒子群算法[3]、遺傳算法[4]等智能算法的出現(xiàn)使水庫優(yōu)化調(diào)度進入一個新階段,但這些算法不可避免地出現(xiàn)魯棒性差、尋優(yōu)精度不高等問題。蝙蝠算法[5]是Xin-She Yang基于微型蝙蝠的回聲定位行為于2010年提出的一種啟發(fā)式算法,該算法結合了粒子群算法的學習機制以及模擬退火算法的冷卻原理等,在計算精度、穩(wěn)定性上優(yōu)于其他算法,具有很好的發(fā)展?jié)摿?與此同時,蝙蝠算法仍然會出現(xiàn)早熟收斂、陷入局部最優(yōu)的問題。自適應權重[6]能夠根據(jù)蝙蝠當前的離散程度與種群進化程度,動態(tài)地調(diào)整位置,從而保證最優(yōu)解不易丟失;連續(xù)禁忌搜索算法[7]作為一種亞啟發(fā)式算法,具有很強的局部搜索能力,能通過對一個初始可行解展開鄰域搜索,并利用一種靈活的記憶結構以及相應的規(guī)則指導算法進行,避免出現(xiàn)迂回搜索從而陷入局部最優(yōu)的情況。筆者嘗試在種群更新機制中引入自適應慣性權重,并根據(jù)種群多樣性的優(yōu)劣[8]結合連續(xù)禁忌搜索策略對蝙蝠算法進行改進,以提高算法在水庫優(yōu)化調(diào)度中的全局尋優(yōu)能力。

    1?水庫發(fā)電優(yōu)化調(diào)度模型

    電站運行效益的提高主要依靠發(fā)電量,因此可以在滿足水電站的水位、水量以及出力等約束條件下,根據(jù)水庫初始運行條件和入庫徑流,合理控制各個時段的發(fā)電流量,使得梯級電站在調(diào)度周期內(nèi)發(fā)電量最大。

    1.1?目標函數(shù)

    以各個時段初的水位作為決策變量,建立發(fā)電量最大的目標函數(shù)。

    F=max∑nk=1∑Tt=1KkQk,tHk,tΔt(1)

    式中:n為梯級水庫的個數(shù);T為調(diào)度時段總數(shù);Kk為水庫k電站綜合出力系數(shù);Qk,t為水庫k在t時段發(fā)電流量,m3/s;Hk,t為水庫k在t時段發(fā)電水頭,m;Δt為時段長度,h。

    1.2?約束條件

    (1)水量平衡約束:

    Vk,t+1=Vk,t+(Ik,t-Qk,t-Sk,t)Δt(2)

    (2)水位約束:

    Zk,tmin≤Zk,t≤Zk,tmax(3)

    (3)水庫下泄流量約束:

    Qk,tmin≤Qk,t+Sk,t≤Qk,tmax(4)

    (4)水電站出力約束:

    Nk,tmin≤Nk,t≤Nk,tmax(5)

    式中:Vk,t為水庫k在t時段初水庫庫容,m3;Ik,t為水庫k在t時段入庫流量,m3/s;Sk,t為水庫k在t時段棄水流量,m3/s;Zk,t為水庫k在t時段初的水位,m;Nk,t為水庫k在t段電站出力,kW;Zk,tmin、Zk,tmax分別為水庫k在t時段允許水位上、下限,m;Qk,tmin、Qk,tmax分別為水庫k在t時段允許下泄流量上、下限,m3/s;Nk,tmin、Nk,tmax分別為水庫k在t時段允許電站出力上、下限,kW。

    2?改進蝙蝠算法

    2.1?標準蝙蝠算法

    蝙蝠算法是模擬蝙蝠利用回聲定位原理來捕獲獵物行為的一種啟發(fā)式算法。假定優(yōu)化問題為max F(X),蝙蝠隨機分布在d維搜索空間中,每只蝙蝠的位置代表函數(shù)F(X)的一個解,對應的函數(shù)值即為蝙蝠的適應度值,適應度值越大則表示蝙蝠的位置越優(yōu)。同時,蝙蝠飛行的過程代表解朝著更優(yōu)方向更新的過程,每只蝙蝠可利用自己所發(fā)出的脈沖頻率以及自身與當前最優(yōu)位置的距離來調(diào)整飛行速度;蝙蝠的脈沖頻率通過在[fmin,fmax]中隨機分配獲得,這個范圍可根據(jù)所研究問題的具體情況進行設定,一般取[0,2]。那么,第i只蝙蝠發(fā)出的脈沖頻率為fi,在g時刻的位置為Xgi并以速度vgi飛行,根據(jù)式(6)~式(8)在g+1時刻獲得新的位置Xg+1i和速度vg+1i。

    fi=fmin+(fmax-fmin)β(6)

    vg+1i=vgi+(Xgi-X*)fi(7)

    Xg+1i=Xgi+vg+1i(8)

    式中:fmin、fmax分別為蝙蝠所發(fā)出聲波的最大頻率、最小頻率;β為區(qū)間[0,1]內(nèi)的一個隨機數(shù);X*為當前最優(yōu)位置。

    一旦蝙蝠的最優(yōu)位置被選中,蝙蝠就按照式(9)在原先最優(yōu)位置Xold附近以隨機游走的方式生成一個新位置Xnew。

    Xnew=Xold+εAg(9)

    式中:ε為區(qū)間(-1,1)內(nèi)的一個隨機數(shù);Ag為g時刻所有的蝙蝠發(fā)出脈沖響度Agi的平均值。

    蝙蝠最初會以較小的脈沖發(fā)射速率和較大的脈沖響度來搜尋獵物,以保證有足夠的搜索范圍;在靠近獵物的過程中為能對獵物進行更精細地搜索,會降低發(fā)出的脈沖響度和提高脈沖發(fā)射速率。每只蝙蝠的初始脈沖響度A0i與初始脈沖發(fā)射速率r0i可以根據(jù)具體問題進行設定,一般取A0i=1表示蝙蝠發(fā)出的響度最大,r0i取接近于0的常數(shù)。若第i只蝙蝠在g時刻隨機游走后最優(yōu)蝙蝠位置得到改善并且發(fā)出的響度大于隨機響度,則將蝙蝠當前最優(yōu)位置更新為Xnew,其脈沖響度以及脈沖發(fā)射速率分別按照式(10)、式(11)進行更新;否則保持最優(yōu)位置、脈沖響度以及脈沖發(fā)射速率不變。

    Ag+1i=αAgi(10)

    rg+1i=r0i[1-exp(-γg)](11)

    式中:Ag+1i為第i只蝙蝠在g+1時刻的脈沖響度;rg+1i為第i只蝙蝠在g+1時刻的脈沖發(fā)射速率;α為聲波衰減系數(shù),一般取0.9;γ為脈沖速率增強系數(shù),一般取0.9。

    當?shù)螖?shù)達到最大值G時,蝙蝠位置停止更新,并輸出蝙蝠的最優(yōu)位置作為本次計算的最優(yōu)解。

    2.2?算法的改進

    2.2.1?自適應權重

    在標準蝙蝠算法中,一旦蝙蝠以某種頻率發(fā)出超聲波后,其速度的更新主要依靠當前全局最優(yōu)位置,步長變化過于單一,容易丟失最優(yōu)解,從而影響蝙蝠算法的搜索性能。在速度更新公式中引入自適應慣性權重,可以根據(jù)蝙蝠當前進化程度以及離散程度,動態(tài)改變蝙蝠的更新速度,平衡算法的全局搜索與局部開發(fā)能力。對于整個種群來說:進化初期應分配較大的慣性權重以保證全局搜索性能,在進化后期應分配較小的慣性權重以保證算法的收斂性能。對于蝙蝠個體來說:個體適應度差于整個種群的平均適應度時,應提高尋優(yōu)速度,分配其較大的慣性權重;個體適應度優(yōu)于整個種群的平均適應度時,應提高局部開發(fā)能力,分配其較小的慣性權重。改進的速度更新公式如下:

    vgi=wgivg-1i+(xgi-x*)fi(12)

    wgi=wmax-(wmax-wmin)[g2G+Fgmax-Fgi2(Fgmax-Fgavg)],

    Fgi≥Fgavgwmax,F(xiàn)gi

    式中:wgi為第i只蝙蝠在g時刻的慣性權重;wmax、wmin分別為慣性權重的最大值、最小值;Fgi為第i只蝙蝠在g時刻的適應度值;Fgmax為種群內(nèi)所有蝙蝠在g時刻的最優(yōu)適應度;Fgavg為種群內(nèi)所有蝙蝠在g時刻的平均適應度值。

    2.2.2?連續(xù)禁忌搜索算法

    隨著迭代次數(shù)的增加,蝙蝠個體會不斷靠近當前最優(yōu)位置而導致種群多樣性急速下降,使得尋優(yōu)能力變差。連續(xù)禁忌搜索算法是禁忌搜索算法在連續(xù)優(yōu)化問題上的一種拓展,是通過模仿人類尋找東西行為的亞啟發(fā)式算法:在一段時間內(nèi)不會對已經(jīng)搜尋過的地方進行重復搜索;若沒有找到,則回到原來的地方繼續(xù)尋找。在蝙蝠算法后期已經(jīng)得到相對較優(yōu)解的情況下,可引入連續(xù)禁忌搜索算法來進一步提高尋優(yōu)精度:從蝙蝠搜索到的最優(yōu)解為中心點(即X*)出發(fā)展開鄰域搜索,對于已搜索到的局部解進行標記和篩選,并在接下來的若干次迭代中避開這些解,從而能保證對解的不同區(qū)域的有效搜索。

    連續(xù)禁忌搜索算法有以下幾個重要因素:

    (1)鄰域。假定當前的鄰域中心點為d維變量X,分別以各個分量xj為中心繪制w個同心超矩形,并在每個同心矩形He(e=1,2,…,w)內(nèi)隨機選取1個點以及在中心矩形H0內(nèi)隨機選取s個點,由這w+s個點共同組成X的鄰域;一般情況下,取w=5,s=3。同心超矩形數(shù)學定義如下:

    H0(x,h0)={x′||xj′-xj|

    xj′

    (14)

    He(x,he-1,he)={x′|he-1<|xj′-xj|

    xjmin

    he=2he-1(16)

    h0=0.001(xjmax-xjmin)(17)

    式中:xjmax、xjmin分別為分量xj的最大值、最小值;he為同心矩形He領域半徑;h0為中心矩形H0的領域半徑。

    (2)禁忌表和禁忌規(guī)則。在連續(xù)禁忌搜索算法中,將需要避開的解稱為禁忌對象,換句話說,這些解處于被禁忌的狀態(tài);禁忌表則是用來存儲禁忌對象的結構,禁忌表能存放禁忌對象的最大個數(shù)稱為禁忌長度。算法運行初始,建立一個空的禁忌表,禁忌長度為L,第一個鄰域中心點X*直接放進禁忌表中。

    每次迭代都采用兩重禁忌規(guī)則來判定鄰域中的點X′在本次迭代中是否被禁忌。首先判斷適應度值F(X):若X′與任意禁忌對象Xl(l=1,2,…,L)的關系均滿足式(18),則X′沒有被禁忌;否則表示X′的適應度值接近于某個禁忌對象X′l,判斷這兩個點是否接近,若此時對于任意的分量xj(j=1,2,…,d),均滿足式(19),則點X′被禁忌,否則沒有被禁忌。

    |F(X)-F(X′l)|≥εf1(18)

    |xj-xj′l|≤εf2(19)

    式中:εf1、εf2均為常數(shù),一般取0.001;L為禁忌長度,一般取常數(shù)。

    如果X′沒有被禁忌,則按照“先進先出”原則放置在禁忌表中,并將下一次迭代的鄰域中心點Xnow更新為X′;當禁忌個數(shù)超過禁忌長度時,將最早進入的禁忌對象剔除。如果X′被禁忌,則不把X′放進禁忌表中,并將X′從鄰域中剔除,從剩余鄰域中選取一個最優(yōu)解,再次進行禁忌判斷,直到找到一個解X″沒有被禁忌。

    (3)特赦規(guī)則。為保證算法的優(yōu)化效果,避免因禁忌規(guī)則而導致算法出錯,還需要設置特赦規(guī)則:當解X′優(yōu)于當前最優(yōu)解Xbest時,為保證解X′不會因禁忌規(guī)則而被漏掉,更新當前最優(yōu)解Xbest后,無視禁忌規(guī)則直接判定X′沒有被禁忌;當鄰域內(nèi)所有解都沒有優(yōu)于當前最優(yōu)解Xbest且都被禁忌,為防止算法出現(xiàn)死循環(huán),將禁忌表內(nèi)的所有禁忌對象取線性平均值,直接判定沒有被禁忌。

    (4)終止規(guī)則。迭代步數(shù)達到最大值,則終止迭代,輸出最優(yōu)解。

    除上述幾個因素外,還要考慮進行禁忌搜索的合適時機:若過早進行禁忌搜索則干擾到蝙蝠的尋優(yōu)過程;若太晚進行禁忌搜索,則會影響最終的尋優(yōu)效果。因此,可以設定一個閥值ξ,對蝙蝠的種群多樣性進行判斷。種群在第g時刻的相對多樣性公式如下:

    ξ(g)=div(g)div(0)

    div(g)=?1md∑dj=1∑mi=1[xji(g)-x*j(g)UBji-LBji]2(20)

    式中:xji(g)為第i只蝙蝠在g時刻的第j個分量;x*j(g)為最優(yōu)位置X*蝙蝠在g時刻的第j個分量;UBji、LBji分別為分量xji(g)的最大值和最小值。

    當ξ(g)小于ξ時,表明蝙蝠的種群多樣性比較差,此時,選取最優(yōu)蝙蝠個體作為進行禁忌搜索的初始點展開進一步尋優(yōu)。

    2.3?改進蝙蝠算法計算流程

    在處理水庫優(yōu)化調(diào)度中的約束條件時,水庫各時段的運行水位約束通過限制蝙蝠的位置范圍來實現(xiàn),而其他約束則引用罰函數(shù)概念處理。改進的蝙蝠算法求解水庫優(yōu)化調(diào)度問題的主要步驟如下:

    第一步:確定參數(shù)。蝙蝠的種群規(guī)模m,脈沖頻率范圍[fmin,fmax],初始脈沖響度A0,初始脈沖發(fā)射速率r0,最大迭代次數(shù)G,種群多樣性閾值ξ,禁忌搜索的禁忌長度L。

    第二步:種群初始化。生成m個可行解作為每只蝙蝠的初始位置X0,計算每只蝙蝠適應度F(Xi),選擇當前的最優(yōu)位置X*。

    第三步:種群迭代。按照式(6)、式(8)、式(12)和式(13)對蝙蝠的位置和飛行速度進行自適應更新,并更新當前最優(yōu)位置X*。

    第四步:局部搜索。選擇最優(yōu)的蝙蝠位置并生成一個隨機數(shù)β1。如果ri小于β1,則按照式(9)在附近進行擾動形成一個局部解,并進行下一步判斷;否則,跳到第六步。

    第五步:更新最優(yōu)位置。生成另一個隨機數(shù)β2,如果最優(yōu)位置得到改善并且Ai大于β2,則更新蝙蝠的最優(yōu)位置,并按照式(10)、式(11)調(diào)整蝙蝠的脈沖響度和脈沖發(fā)射速率;否則,保持不變。

    第六步:種群多樣性判斷。當前蝙蝠多樣性若滿足閾值要求,則進行第七步;否則,跳到第八步。

    第七步:如果迭代次數(shù)達到最大,則輸出最優(yōu)結果;否則,回到第三步。

    第八步:禁忌搜索。從當前最優(yōu)解X*出發(fā)進行連續(xù)禁忌搜索,根據(jù)禁忌規(guī)則以及特赦規(guī)則,更新禁忌表、鄰域以及當前最優(yōu)解。

    第九步:如果迭代次數(shù)達到最大,則輸出最優(yōu)結果;否則,回到第八步。

    改進后的蝙蝠算法流程圖見圖1。

    3?實例應用

    我國某梯級系統(tǒng)上游至下游包含2個水電站A、B,兩個電站均以發(fā)電為主要任務,電站間有區(qū)間入流匯入。電站的基本參數(shù)見表1。

    由于2個電站在汛期基本保持在汛限水位運行,因此針對供水期進行優(yōu)化調(diào)度分析。選用梯級電站系統(tǒng)3個典型年供水期(12月至次年5月,共6個月)的徑流資料,以旬為計算時段,分析梯級水庫的發(fā)電優(yōu)化調(diào)度結果。

    設置蝙蝠種群大小m=100,蝙蝠發(fā)出的脈沖頻率范圍[fmin,fmax]為[0,2],初始脈沖響度A0=1,初始脈沖速率r0=0.001,聲波衰減系數(shù)α=0.9,脈沖速率增強系數(shù)γ=0.9;在改進蝙蝠算法中最大慣性權重wmax=1,最小慣性權重wmin=0.4,種群相對多樣性閾值ξ取0.001,禁忌搜索的禁忌長度L=6。標準蝙蝠算法和改進蝙蝠算法計算得到的梯級總發(fā)電量優(yōu)化結果見表2,相應的3個典型年供水期的水庫A、B優(yōu)化調(diào)度水位過程如圖2~圖4所示。

    由以上結果可知,改進蝙蝠算法比標準蝙蝠算法計算得到的梯級發(fā)電量有所增加,枯水年增加1.90×106 kW·h,平水年增加9.57×105 kW·h,豐水年增加1.15×105 kW·h。表明從梯級發(fā)電量考慮,改進蝙蝠算法的優(yōu)化效果更好。

    為了進一步比較這兩種優(yōu)化算法的收斂速度,這里以枯水年供水期計算為例,給出兩種算法的迭代過程,如圖5所示。

    由圖5可知,標準蝙蝠算法達到收斂所需迭代次數(shù)大于改進蝙蝠算法所需迭代次數(shù),表明改進蝙蝠算法在獲得更優(yōu)解的前提下,收斂速度明顯提高。

    4?結?論

    在標準蝙蝠算法中引入自適應權重,改善蝙蝠種群的更新機制,使蝙蝠位置能夠根據(jù)實際情況進行自適應動態(tài)調(diào)整,提高蝙蝠的搜索性能;同時,在算法后期種群多樣性變差時,引入連續(xù)禁忌搜索算法,避免出現(xiàn)迂回搜索而無法跳出局部最優(yōu)的情況。案例分析表明,改進蝙蝠算法與標準蝙蝠算法相比,能以較快的運行效率獲得更優(yōu)的水庫運行方案,為水庫發(fā)電優(yōu)化調(diào)度提供了一個新的途徑。

    參考文獻:

    [1]?楊峰,黃懷禮,張強.用動態(tài)規(guī)劃法對水庫進行優(yōu)化調(diào)度[J].河南科學,2005,23(1):17-19.

    [2]?羅紅兵.POA算法在水庫優(yōu)化調(diào)度中的應用[J].陜西水利,2018(6):127-129.

    [3]?鄧顯羽,彭勇,葉碎高,等.粒子群算法在水庫(群)優(yōu)化調(diào)度研究中的應用綜述[J].水利水電科技進展,2010,30(5):90-94.

    [4]?馮迅,王金文,權先璋,等.遺傳算法在水電站優(yōu)化調(diào)度中的實用研究[J].華中電力,2005(3):12-15.

    [5]?YANG X S. A New Metaheuristic Bat-inspired Algorithm[J].Computer Knowledge &Technology,2010,284:65-74.

    [6]?劉列.自適應粒子群算法在水庫優(yōu)化調(diào)度中的應用研究[J].大眾科技,2017,19(1):11-13.

    [7]?王明興.連續(xù)禁忌搜索算法改進及應用研究[D].杭州:浙江大學,2005:15-27.

    [8]?紀昌明,劉方,喻杉,等.基于鯰魚效應粒子群算法的梯級水庫群優(yōu)化調(diào)度[J].電力系統(tǒng)保護與控制,2011,39(19):63-68.

    【責任編輯?崔瀟菡】

    18禁在线播放成人免费| 99精国产麻豆久久婷婷| 日韩精品免费视频一区二区三区 | 亚洲精品,欧美精品| 精品一区二区免费观看| 王馨瑶露胸无遮挡在线观看| 99热这里只有精品一区| 成年人免费黄色播放视频 | 男人和女人高潮做爰伦理| 亚洲图色成人| 在线精品无人区一区二区三| av在线观看视频网站免费| 有码 亚洲区| 高清欧美精品videossex| 精品人妻一区二区三区麻豆| 99精国产麻豆久久婷婷| 久久久久视频综合| 亚洲色图综合在线观看| 国产亚洲欧美精品永久| 日本91视频免费播放| 99热全是精品| 欧美3d第一页| 成人二区视频| 欧美 日韩 精品 国产| 久久精品国产a三级三级三级| 亚洲人成网站在线播| 久久精品国产亚洲av天美| 欧美日韩亚洲高清精品| 精品99又大又爽又粗少妇毛片| 日本91视频免费播放| av在线app专区| h视频一区二区三区| 日韩成人伦理影院| 亚洲欧美精品专区久久| 噜噜噜噜噜久久久久久91| 下体分泌物呈黄色| 久久毛片免费看一区二区三区| 中文天堂在线官网| 精品一区在线观看国产| 色视频www国产| 蜜桃在线观看..| 美女大奶头黄色视频| 国产成人精品婷婷| 国产美女午夜福利| 在线播放无遮挡| 国产亚洲5aaaaa淫片| 精品一区在线观看国产| 韩国高清视频一区二区三区| 日本黄色日本黄色录像| av免费观看日本| 天美传媒精品一区二区| 嘟嘟电影网在线观看| 欧美人与善性xxx| 日本与韩国留学比较| 国产精品一区二区在线不卡| 国产精品伦人一区二区| 久久久久久久久久久免费av| 人人妻人人澡人人看| 国精品久久久久久国模美| 丝袜在线中文字幕| 人人妻人人澡人人爽人人夜夜| 日本wwww免费看| 久久韩国三级中文字幕| 亚洲美女搞黄在线观看| 久久久久久久久久人人人人人人| 亚洲综合色惰| 国产老妇伦熟女老妇高清| 日韩视频在线欧美| 亚洲精品自拍成人| 天美传媒精品一区二区| 又爽又黄a免费视频| 成人亚洲欧美一区二区av| 菩萨蛮人人尽说江南好唐韦庄| 男男h啪啪无遮挡| 内地一区二区视频在线| 肉色欧美久久久久久久蜜桃| 少妇的逼好多水| 欧美日韩av久久| 久热这里只有精品99| 日韩欧美 国产精品| 九色成人免费人妻av| 人妻少妇偷人精品九色| 国产又色又爽无遮挡免| 五月伊人婷婷丁香| 一区二区三区乱码不卡18| 久久久国产欧美日韩av| 国产永久视频网站| 少妇 在线观看| 国产又色又爽无遮挡免| 亚洲天堂av无毛| 久久精品久久精品一区二区三区| 久久久久久人妻| 黄色日韩在线| xxx大片免费视频| 男人爽女人下面视频在线观看| 美女中出高潮动态图| 欧美精品人与动牲交sv欧美| 成人特级av手机在线观看| 午夜激情福利司机影院| 日日摸夜夜添夜夜添av毛片| 人体艺术视频欧美日本| 一级毛片黄色毛片免费观看视频| 免费黄网站久久成人精品| av黄色大香蕉| 国产黄片美女视频| 亚洲欧美精品自产自拍| 中文欧美无线码| 这个男人来自地球电影免费观看 | 成人18禁高潮啪啪吃奶动态图 | 美女主播在线视频| 七月丁香在线播放| 久久久久人妻精品一区果冻| 99久久综合免费| 国产中年淑女户外野战色| 国产爽快片一区二区三区| 国产高清三级在线| 高清视频免费观看一区二区| 久久久久久久久久人人人人人人| 亚洲精品视频女| 国产日韩一区二区三区精品不卡 | 又爽又黄a免费视频| 国产免费福利视频在线观看| av一本久久久久| 热re99久久精品国产66热6| 日韩,欧美,国产一区二区三区| 久久精品久久久久久噜噜老黄| videossex国产| 精品一区二区三卡| 久久国产乱子免费精品| 我要看黄色一级片免费的| 精品一品国产午夜福利视频| 亚洲真实伦在线观看| 久久6这里有精品| 乱码一卡2卡4卡精品| 国产精品一区二区性色av| 丰满乱子伦码专区| 美女视频免费永久观看网站| 亚洲av成人精品一区久久| 中文字幕精品免费在线观看视频 | 亚洲av成人精品一区久久| 日韩欧美 国产精品| 美女xxoo啪啪120秒动态图| 特大巨黑吊av在线直播| 成人毛片a级毛片在线播放| 国产高清国产精品国产三级| 久久久久久久久久久丰满| 亚洲第一区二区三区不卡| 精品久久国产蜜桃| 99热这里只有精品一区| 亚洲精品视频女| 亚洲,一卡二卡三卡| 极品少妇高潮喷水抽搐| 亚洲精品久久久久久婷婷小说| 最近中文字幕高清免费大全6| 老熟女久久久| 国产黄色免费在线视频| 精品久久久精品久久久| 日韩在线高清观看一区二区三区| 亚洲综合色惰| 另类精品久久| 亚洲国产精品一区二区三区在线| 久久久久久久精品精品| 成人综合一区亚洲| 观看av在线不卡| 高清欧美精品videossex| 中文字幕免费在线视频6| 国产探花极品一区二区| 亚洲国产欧美在线一区| 欧美成人精品欧美一级黄| 波野结衣二区三区在线| 亚洲av成人精品一二三区| 91精品国产国语对白视频| 国产91av在线免费观看| 岛国毛片在线播放| 少妇的逼好多水| 国产精品久久久久久久久免| 纵有疾风起免费观看全集完整版| 最近中文字幕高清免费大全6| 丁香六月天网| 国产真实伦视频高清在线观看| 日韩免费高清中文字幕av| 国产又色又爽无遮挡免| 亚洲精品色激情综合| 水蜜桃什么品种好| 日韩精品免费视频一区二区三区 | 午夜福利视频精品| 在线观看免费视频网站a站| 国产熟女午夜一区二区三区 | 精品亚洲成a人片在线观看| 国产女主播在线喷水免费视频网站| 国产日韩一区二区三区精品不卡 | 国产午夜精品一二区理论片| 九九久久精品国产亚洲av麻豆| 中文乱码字字幕精品一区二区三区| 亚洲无线观看免费| 大码成人一级视频| 成年av动漫网址| 我的老师免费观看完整版| 美女xxoo啪啪120秒动态图| 中文在线观看免费www的网站| 最新的欧美精品一区二区| 精品国产国语对白av| 久热这里只有精品99| 亚洲av成人精品一二三区| 狂野欧美激情性bbbbbb| 大香蕉久久网| 在线观看国产h片| 亚洲精品乱码久久久v下载方式| 国内少妇人妻偷人精品xxx网站| 久久青草综合色| 国产成人精品一,二区| 国产精品成人在线| av女优亚洲男人天堂| 大香蕉久久网| 人人澡人人妻人| 丝袜喷水一区| 26uuu在线亚洲综合色| 2021少妇久久久久久久久久久| 国产伦理片在线播放av一区| 国产亚洲欧美精品永久| 精品亚洲成国产av| 看非洲黑人一级黄片| 中文乱码字字幕精品一区二区三区| 看免费成人av毛片| av天堂久久9| 免费黄频网站在线观看国产| 日韩一本色道免费dvd| 七月丁香在线播放| 国产精品久久久久久精品电影小说| 午夜福利,免费看| 亚洲一级一片aⅴ在线观看| 亚洲成人一二三区av| 午夜福利在线观看免费完整高清在| 亚洲欧洲精品一区二区精品久久久 | 国产精品熟女久久久久浪| 日韩,欧美,国产一区二区三区| 日韩av免费高清视频| 99热这里只有精品一区| 欧美日韩视频精品一区| 人妻系列 视频| 一级毛片黄色毛片免费观看视频| 国产成人a∨麻豆精品| 蜜桃在线观看..| 国产精品国产三级国产av玫瑰| 亚洲精品中文字幕在线视频 | 又黄又爽又刺激的免费视频.| 国产黄色免费在线视频| 亚洲欧美中文字幕日韩二区| 熟女电影av网| videossex国产| 精品少妇黑人巨大在线播放| 国产亚洲午夜精品一区二区久久| 国产亚洲av片在线观看秒播厂| 国产一区二区三区av在线| 久久青草综合色| 亚洲精品中文字幕在线视频 | 精品少妇内射三级| 看免费成人av毛片| 观看美女的网站| 精品一区二区三卡| 我的女老师完整版在线观看| 中文字幕亚洲精品专区| 久久99蜜桃精品久久| 亚洲av免费高清在线观看| 国产探花极品一区二区| 久久人人爽av亚洲精品天堂| 欧美高清成人免费视频www| 亚洲一级一片aⅴ在线观看| 18禁动态无遮挡网站| 亚洲av国产av综合av卡| 国产淫片久久久久久久久| 看十八女毛片水多多多| 校园人妻丝袜中文字幕| 伦精品一区二区三区| 亚洲av福利一区| 内地一区二区视频在线| 日本色播在线视频| 天堂俺去俺来也www色官网| 99视频精品全部免费 在线| 久久综合国产亚洲精品| 欧美xxxx性猛交bbbb| 久久久久久人妻| 午夜视频国产福利| 一级毛片黄色毛片免费观看视频| 国产成人91sexporn| 国产av国产精品国产| a级毛色黄片| 日日摸夜夜添夜夜添av毛片| 如日韩欧美国产精品一区二区三区 | 内射极品少妇av片p| 亚洲天堂av无毛| 热re99久久国产66热| 久久久久网色| 夜夜爽夜夜爽视频| 久久久午夜欧美精品| 天美传媒精品一区二区| 欧美+日韩+精品| 亚洲精品国产成人久久av| 久久久a久久爽久久v久久| 最黄视频免费看| 日韩一区二区三区影片| 亚洲精品视频女| 久久久久网色| 亚洲精品aⅴ在线观看| 午夜老司机福利剧场| 女人久久www免费人成看片| 青青草视频在线视频观看| 另类亚洲欧美激情| 亚洲欧洲精品一区二区精品久久久 | 欧美 日韩 精品 国产| 汤姆久久久久久久影院中文字幕| 最近2019中文字幕mv第一页| 日韩一本色道免费dvd| 女性被躁到高潮视频| 色网站视频免费| 午夜福利视频精品| 亚洲av二区三区四区| 午夜免费鲁丝| 国产成人freesex在线| 噜噜噜噜噜久久久久久91| 午夜免费观看性视频| 97超视频在线观看视频| 国产亚洲av片在线观看秒播厂| 国产精品秋霞免费鲁丝片| 美女视频免费永久观看网站| 69精品国产乱码久久久| 晚上一个人看的免费电影| 久久99精品国语久久久| 韩国av在线不卡| 一级,二级,三级黄色视频| 在现免费观看毛片| 国产女主播在线喷水免费视频网站| 一级毛片我不卡| 国产男女超爽视频在线观看| 国精品久久久久久国模美| 久久精品国产亚洲av涩爱| 偷拍熟女少妇极品色| 日韩电影二区| 国产精品无大码| 亚洲成色77777| 日本wwww免费看| 国产精品嫩草影院av在线观看| 国产日韩欧美在线精品| 国产黄色免费在线视频| 最后的刺客免费高清国语| 另类亚洲欧美激情| 亚洲一级一片aⅴ在线观看| 国产精品一区www在线观看| 一级爰片在线观看| 久久久亚洲精品成人影院| 亚洲高清免费不卡视频| 国产精品.久久久| 欧美性感艳星| 蜜桃在线观看..| 国产视频首页在线观看| 亚洲性久久影院| 亚洲情色 制服丝袜| 免费看不卡的av| 欧美国产精品一级二级三级 | 十分钟在线观看高清视频www | 国产真实伦视频高清在线观看| 建设人人有责人人尽责人人享有的| 日韩av在线免费看完整版不卡| 丝袜脚勾引网站| 久久99一区二区三区| 中文乱码字字幕精品一区二区三区| 男女免费视频国产| 亚洲国产精品成人久久小说| 国产免费一级a男人的天堂| 久久国内精品自在自线图片| 免费看光身美女| 国产精品人妻久久久影院| 色婷婷久久久亚洲欧美| 一二三四中文在线观看免费高清| 青春草视频在线免费观看| 国产成人91sexporn| 男男h啪啪无遮挡| 人妻一区二区av| 97超视频在线观看视频| 国产在线一区二区三区精| 三级经典国产精品| 伦理电影免费视频| 又爽又黄a免费视频| 国产真实伦视频高清在线观看| 简卡轻食公司| 国产免费一级a男人的天堂| 好男人视频免费观看在线| 亚洲精品一二三| 亚洲自偷自拍三级| 国产乱人偷精品视频| 日本av免费视频播放| 美女内射精品一级片tv| 人妻 亚洲 视频| 亚洲美女视频黄频| 日韩 亚洲 欧美在线| 内射极品少妇av片p| 成人国产麻豆网| 熟女人妻精品中文字幕| 色5月婷婷丁香| 极品人妻少妇av视频| 国产美女午夜福利| 欧美+日韩+精品| 国产欧美日韩综合在线一区二区 | 男人爽女人下面视频在线观看| 久久久久久久久久久免费av| 亚洲av福利一区| 国产日韩欧美亚洲二区| av又黄又爽大尺度在线免费看| 夜夜爽夜夜爽视频| av天堂久久9| 亚洲久久久国产精品| 免费观看无遮挡的男女| 色网站视频免费| 亚洲av免费高清在线观看| 纵有疾风起免费观看全集完整版| 国产乱来视频区| 中文欧美无线码| 久久久久久久国产电影| 一本色道久久久久久精品综合| 日韩视频在线欧美| 中国国产av一级| 日本黄色日本黄色录像| 亚洲高清免费不卡视频| 我的老师免费观看完整版| 久久久久国产精品人妻一区二区| 国产精品成人在线| 亚洲精品乱码久久久v下载方式| 一级毛片 在线播放| 黄色欧美视频在线观看| 成人18禁高潮啪啪吃奶动态图 | 制服丝袜香蕉在线| 国产成人aa在线观看| 久久99蜜桃精品久久| 国产欧美日韩一区二区三区在线 | 最后的刺客免费高清国语| 男女免费视频国产| 99热这里只有是精品在线观看| 国产 一区精品| 在线亚洲精品国产二区图片欧美 | 2021少妇久久久久久久久久久| 欧美三级亚洲精品| 一区二区三区免费毛片| 亚洲av成人精品一区久久| 丝袜喷水一区| 亚洲精品日韩在线中文字幕| 91久久精品国产一区二区成人| 久久婷婷青草| 亚洲国产av新网站| 欧美丝袜亚洲另类| 制服丝袜香蕉在线| av国产久精品久网站免费入址| 日韩一区二区视频免费看| 草草在线视频免费看| 国产日韩欧美视频二区| 插逼视频在线观看| 少妇猛男粗大的猛烈进出视频| 日韩欧美一区视频在线观看 | 欧美丝袜亚洲另类| 男女啪啪激烈高潮av片| 99久久精品国产国产毛片| 久久女婷五月综合色啪小说| 欧美日韩国产mv在线观看视频| 国产免费一级a男人的天堂| 国产日韩欧美视频二区| 精品国产一区二区久久| av专区在线播放| 亚洲情色 制服丝袜| 亚洲av二区三区四区| 久久久久网色| 少妇人妻一区二区三区视频| 丁香六月天网| 成人影院久久| 免费av不卡在线播放| 久久久久网色| 免费看av在线观看网站| 日韩伦理黄色片| 老女人水多毛片| 欧美另类一区| 看免费成人av毛片| 女性被躁到高潮视频| 国产精品偷伦视频观看了| 亚洲图色成人| 欧美另类一区| h日本视频在线播放| 日本欧美国产在线视频| 中文欧美无线码| 极品少妇高潮喷水抽搐| 狂野欧美激情性xxxx在线观看| 国产精品一区二区三区四区免费观看| 色婷婷久久久亚洲欧美| 亚洲熟女精品中文字幕| 欧美精品高潮呻吟av久久| 久久影院123| 五月开心婷婷网| 色婷婷久久久亚洲欧美| 哪个播放器可以免费观看大片| 国产色爽女视频免费观看| 国产精品一区二区在线不卡| 国产高清不卡午夜福利| 日韩 亚洲 欧美在线| 少妇被粗大猛烈的视频| 国产一区有黄有色的免费视频| 成人18禁高潮啪啪吃奶动态图 | 久久热精品热| 免费黄网站久久成人精品| 高清av免费在线| 人人妻人人爽人人添夜夜欢视频 | 国产精品.久久久| 精品一区二区三区视频在线| 观看av在线不卡| 国产免费又黄又爽又色| 少妇裸体淫交视频免费看高清| 成人亚洲欧美一区二区av| 最近的中文字幕免费完整| 看非洲黑人一级黄片| 精品一区二区三卡| 久久6这里有精品| 永久免费av网站大全| 婷婷色av中文字幕| 亚洲国产欧美在线一区| 欧美日本中文国产一区发布| 久久午夜福利片| 欧美日韩综合久久久久久| 国语对白做爰xxxⅹ性视频网站| 99精国产麻豆久久婷婷| 国产亚洲精品久久久com| 国产成人精品久久久久久| 亚洲欧美中文字幕日韩二区| a级毛片在线看网站| 久久精品夜色国产| 国产永久视频网站| 人妻系列 视频| 一区在线观看完整版| 亚洲国产成人一精品久久久| 99久久中文字幕三级久久日本| 亚洲欧美成人精品一区二区| 精品久久久精品久久久| 国产精品成人在线| 99热这里只有是精品50| 最近中文字幕高清免费大全6| 大话2 男鬼变身卡| 男人狂女人下面高潮的视频| 国产成人精品婷婷| 伦理电影大哥的女人| 乱系列少妇在线播放| 少妇人妻 视频| 午夜免费男女啪啪视频观看| 王馨瑶露胸无遮挡在线观看| 免费高清在线观看视频在线观看| 97在线人人人人妻| 久久综合国产亚洲精品| 国产免费福利视频在线观看| 国产精品蜜桃在线观看| 多毛熟女@视频| .国产精品久久| 国产精品国产三级国产专区5o| 日韩大片免费观看网站| 国产成人一区二区在线| 老司机影院成人| 最新中文字幕久久久久| 久久毛片免费看一区二区三区| 亚洲国产最新在线播放| 黄片无遮挡物在线观看| 精品一区二区三卡| 毛片一级片免费看久久久久| 日产精品乱码卡一卡2卡三| 一个人免费看片子| av在线播放精品| 天天操日日干夜夜撸| 男人狂女人下面高潮的视频| 亚洲欧美成人综合另类久久久| 老女人水多毛片| 91久久精品国产一区二区三区| 亚洲在久久综合| 国产精品久久久久久av不卡| 99精国产麻豆久久婷婷| 国内精品宾馆在线| 男女边摸边吃奶| 香蕉精品网在线| 少妇人妻久久综合中文| 亚洲丝袜综合中文字幕| 国产综合精华液| 日韩熟女老妇一区二区性免费视频| 六月丁香七月| av.在线天堂| 99热6这里只有精品| 99热国产这里只有精品6| 中文字幕人妻熟人妻熟丝袜美| 免费不卡的大黄色大毛片视频在线观看| 日韩视频在线欧美| 偷拍熟女少妇极品色| 女性生殖器流出的白浆| 在线 av 中文字幕| 亚洲美女黄色视频免费看| 亚洲精品乱码久久久v下载方式| 亚洲一级一片aⅴ在线观看| av线在线观看网站| 中国三级夫妇交换| 人人妻人人看人人澡| 黄色毛片三级朝国网站 | av国产久精品久网站免费入址| 我的老师免费观看完整版| 久久99热6这里只有精品| 亚洲美女黄色视频免费看| 在线亚洲精品国产二区图片欧美 | 一级毛片电影观看| 岛国毛片在线播放| 精品一区二区三区视频在线| 免费观看a级毛片全部| av线在线观看网站| 久久人人爽人人爽人人片va| a级片在线免费高清观看视频| 午夜免费鲁丝| 欧美日韩综合久久久久久| 国产黄色视频一区二区在线观看|