• <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.

    【責任編輯?崔瀟菡】

    久久久色成人| 日韩欧美国产一区二区入口| 免费在线观看亚洲国产| 变态另类丝袜制服| 欧美bdsm另类| 搡老岳熟女国产| av中文乱码字幕在线| 国产爱豆传媒在线观看| 波多野结衣高清无吗| 国产亚洲精品久久久久久毛片| 一区二区三区国产精品乱码| 欧美成人免费av一区二区三区| 少妇裸体淫交视频免费看高清| 丁香欧美五月| 亚洲va日本ⅴa欧美va伊人久久| 久久国产精品人妻蜜桃| 国产成人aa在线观看| 制服丝袜大香蕉在线| 丰满的人妻完整版| 丰满的人妻完整版| 精华霜和精华液先用哪个| 在线观看免费视频日本深夜| 小说图片视频综合网站| 窝窝影院91人妻| 日本黄色片子视频| 精品熟女少妇八av免费久了| 久久久色成人| 午夜精品久久久久久毛片777| 美女cb高潮喷水在线观看| 亚洲精品美女久久久久99蜜臀| 亚洲精品成人久久久久久| 免费看a级黄色片| 亚洲成人精品中文字幕电影| 欧美日韩综合久久久久久 | 久久午夜亚洲精品久久| 五月玫瑰六月丁香| 亚洲精品一区av在线观看| 一本综合久久免费| 少妇的逼水好多| 狠狠狠狠99中文字幕| 欧美zozozo另类| 有码 亚洲区| 国产三级中文精品| 国产精品1区2区在线观看.| 欧美日韩国产亚洲二区| 高清在线国产一区| 久久久国产成人免费| 美女cb高潮喷水在线观看| 俺也久久电影网| 久久久久国内视频| 偷拍熟女少妇极品色| 久久久久久久午夜电影| av在线天堂中文字幕| 我要搜黄色片| 精品免费久久久久久久清纯| 精品国产三级普通话版| 日本成人三级电影网站| 波多野结衣巨乳人妻| 成人亚洲精品av一区二区| 国产精品美女特级片免费视频播放器| 日本免费a在线| 国产成人福利小说| 88av欧美| 嫩草影视91久久| 亚洲成人免费电影在线观看| 99久久99久久久精品蜜桃| 亚洲av不卡在线观看| 丝袜美腿在线中文| 亚洲一区高清亚洲精品| 久久精品91无色码中文字幕| 嫁个100分男人电影在线观看| 天天躁日日操中文字幕| 天美传媒精品一区二区| 香蕉丝袜av| 国产一区二区激情短视频| 男女下面进入的视频免费午夜| 动漫黄色视频在线观看| 麻豆成人av在线观看| 欧美精品啪啪一区二区三区| 一级毛片高清免费大全| 麻豆成人午夜福利视频| 夜夜爽天天搞| 99热只有精品国产| 国产精品乱码一区二三区的特点| 国产在视频线在精品| 亚洲欧美日韩高清专用| 欧美黑人巨大hd| 淫秽高清视频在线观看| 美女高潮喷水抽搐中文字幕| 男女下面进入的视频免费午夜| 午夜免费成人在线视频| 日韩人妻高清精品专区| 麻豆成人av在线观看| 身体一侧抽搐| 国产精品98久久久久久宅男小说| 757午夜福利合集在线观看| 久久久久国产精品人妻aⅴ院| 99久久99久久久精品蜜桃| 一区二区三区国产精品乱码| 欧美日本视频| 欧美日韩国产亚洲二区| 日韩国内少妇激情av| 精品不卡国产一区二区三区| 人妻夜夜爽99麻豆av| 婷婷丁香在线五月| 亚洲av电影不卡..在线观看| 真实男女啪啪啪动态图| 午夜日韩欧美国产| 久久国产精品人妻蜜桃| 嫩草影院精品99| 舔av片在线| 乱人视频在线观看| 一区福利在线观看| avwww免费| 亚洲国产精品999在线| 国内少妇人妻偷人精品xxx网站| 成人国产综合亚洲| 国产亚洲精品一区二区www| 色综合婷婷激情| 在线国产一区二区在线| 成熟少妇高潮喷水视频| 特大巨黑吊av在线直播| 非洲黑人性xxxx精品又粗又长| 午夜激情欧美在线| 性欧美人与动物交配| 97超级碰碰碰精品色视频在线观看| 亚洲人成伊人成综合网2020| 成年女人毛片免费观看观看9| 国产成人啪精品午夜网站| 午夜精品一区二区三区免费看| 一区二区三区国产精品乱码| 亚洲av五月六月丁香网| 身体一侧抽搐| 悠悠久久av| av国产免费在线观看| avwww免费| 色吧在线观看| 精品国产三级普通话版| 男女下面进入的视频免费午夜| 美女黄网站色视频| 日本精品一区二区三区蜜桃| 国产一级毛片七仙女欲春2| 精品一区二区三区av网在线观看| 国产av一区在线观看免费| 国产97色在线日韩免费| 国产亚洲欧美98| 嫁个100分男人电影在线观看| 久久九九热精品免费| 成人无遮挡网站| 午夜精品在线福利| 一本久久中文字幕| 国产高清有码在线观看视频| 成年女人永久免费观看视频| 有码 亚洲区| 免费电影在线观看免费观看| 午夜福利18| 亚洲乱码一区二区免费版| 久久久久免费精品人妻一区二区| 欧美国产日韩亚洲一区| 国产久久久一区二区三区| 男女视频在线观看网站免费| 麻豆国产97在线/欧美| 免费观看人在逋| 三级毛片av免费| 一级作爱视频免费观看| 日韩欧美精品免费久久 | 国产高清视频在线播放一区| av在线蜜桃| 婷婷亚洲欧美| 成人高潮视频无遮挡免费网站| 少妇高潮的动态图| 老熟妇乱子伦视频在线观看| 精品久久久久久久久久免费视频| 久久久久久久亚洲中文字幕 | 久久久久性生活片| 一级作爱视频免费观看| 级片在线观看| 亚洲欧美日韩高清在线视频| 国产激情偷乱视频一区二区| 啦啦啦观看免费观看视频高清| 日韩大尺度精品在线看网址| 啦啦啦韩国在线观看视频| 日本黄色视频三级网站网址| 少妇高潮的动态图| 99久久久亚洲精品蜜臀av| 精品人妻1区二区| 男女下面进入的视频免费午夜| 欧美黄色片欧美黄色片| АⅤ资源中文在线天堂| 99国产综合亚洲精品| 亚洲国产精品成人综合色| 日韩高清综合在线| 1024手机看黄色片| 久久草成人影院| 法律面前人人平等表现在哪些方面| 成人一区二区视频在线观看| 国产 一区 欧美 日韩| 欧美日本视频| 国产毛片a区久久久久| 999久久久精品免费观看国产| 白带黄色成豆腐渣| 又紧又爽又黄一区二区| 身体一侧抽搐| 亚洲真实伦在线观看| 精品电影一区二区在线| ponron亚洲| 宅男免费午夜| 99国产综合亚洲精品| 欧美日韩一级在线毛片| 别揉我奶头~嗯~啊~动态视频| 亚洲久久久久久中文字幕| 久久久久免费精品人妻一区二区| 看黄色毛片网站| 欧美日韩综合久久久久久 | 国产av不卡久久| 久久久久久久精品吃奶| 日韩亚洲欧美综合| 久99久视频精品免费| 一本久久中文字幕| 在线观看免费午夜福利视频| 色av中文字幕| 中文字幕久久专区| 成人精品一区二区免费| 国产探花极品一区二区| 全区人妻精品视频| 亚洲精品一区av在线观看| 男女做爰动态图高潮gif福利片| 91麻豆精品激情在线观看国产| 有码 亚洲区| 亚洲欧美日韩卡通动漫| 国产中年淑女户外野战色| 精品久久久久久成人av| 首页视频小说图片口味搜索| 免费看日本二区| 亚洲国产精品久久男人天堂| 淫妇啪啪啪对白视频| 久久久久免费精品人妻一区二区| 麻豆成人av在线观看| 日韩精品青青久久久久久| 一a级毛片在线观看| 久久久久久国产a免费观看| 欧美成人性av电影在线观看| 亚洲av免费高清在线观看| 精品人妻一区二区三区麻豆 | av视频在线观看入口| 免费看十八禁软件| 成年免费大片在线观看| 一区二区三区国产精品乱码| 日本黄大片高清| 亚洲国产精品sss在线观看| www.999成人在线观看| 亚洲欧美一区二区三区黑人| 99在线视频只有这里精品首页| 精品熟女少妇八av免费久了| 成年女人毛片免费观看观看9| 午夜视频国产福利| 亚洲av二区三区四区| 欧美一区二区精品小视频在线| 无遮挡黄片免费观看| 国产一区二区亚洲精品在线观看| 桃红色精品国产亚洲av| 夜夜夜夜夜久久久久| 久久久久久久亚洲中文字幕 | 午夜a级毛片| 亚洲美女视频黄频| 19禁男女啪啪无遮挡网站| 舔av片在线| 欧美乱妇无乱码| 久久久久精品国产欧美久久久| 日韩欧美在线乱码| 欧美日韩综合久久久久久 | 搡老妇女老女人老熟妇| 欧美成狂野欧美在线观看| 欧美另类亚洲清纯唯美| 精品国内亚洲2022精品成人| 国产伦一二天堂av在线观看| 精品一区二区三区人妻视频| 欧美av亚洲av综合av国产av| 色综合站精品国产| 亚洲五月婷婷丁香| 国产视频内射| 黑人欧美特级aaaaaa片| 波野结衣二区三区在线 | 美女黄网站色视频| 一个人免费在线观看电影| 国产日本99.免费观看| 国产精品野战在线观看| 在线a可以看的网站| 美女免费视频网站| 最近视频中文字幕2019在线8| 国产亚洲av嫩草精品影院| 法律面前人人平等表现在哪些方面| 一本综合久久免费| 在线视频色国产色| 欧美一区二区亚洲| 一级黄片播放器| 午夜视频国产福利| 一区二区三区高清视频在线| 少妇的逼好多水| 国产精品嫩草影院av在线观看 | 少妇高潮的动态图| 亚洲专区中文字幕在线| 五月玫瑰六月丁香| 一本久久中文字幕| 国产成人啪精品午夜网站| 在线播放无遮挡| 三级国产精品欧美在线观看| 欧美午夜高清在线| 一个人观看的视频www高清免费观看| 欧美丝袜亚洲另类 | 亚洲狠狠婷婷综合久久图片| 亚洲成人久久爱视频| 欧美大码av| 老熟妇乱子伦视频在线观看| 国产三级在线视频| 丰满乱子伦码专区| 久久久色成人| 天堂影院成人在线观看| 韩国av一区二区三区四区| 免费在线观看影片大全网站| 18禁黄网站禁片免费观看直播| 色综合站精品国产| 一级作爱视频免费观看| 亚洲人成网站高清观看| 国产精品国产高清国产av| 亚洲真实伦在线观看| 国产三级黄色录像| 国产精华一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 久久国产精品人妻蜜桃| 精品一区二区三区人妻视频| 丰满人妻一区二区三区视频av | 美女高潮的动态| www日本在线高清视频| 国产爱豆传媒在线观看| 色吧在线观看| 欧美成狂野欧美在线观看| 国产成人av激情在线播放| 欧美日韩精品网址| 少妇丰满av| 18禁裸乳无遮挡免费网站照片| 99久久成人亚洲精品观看| 亚洲 欧美 日韩 在线 免费| 有码 亚洲区| 桃色一区二区三区在线观看| 中文在线观看免费www的网站| 成年版毛片免费区| 深夜精品福利| 午夜激情福利司机影院| 久久久久免费精品人妻一区二区| 美女大奶头视频| 亚洲精华国产精华精| 九色国产91popny在线| 国内精品美女久久久久久| 女生性感内裤真人,穿戴方法视频| 亚洲成av人片免费观看| 日韩欧美在线乱码| 3wmmmm亚洲av在线观看| 18美女黄网站色大片免费观看| 国产在视频线在精品| 精品国产超薄肉色丝袜足j| 有码 亚洲区| 午夜激情欧美在线| 国产亚洲精品一区二区www| 日韩人妻高清精品专区| 噜噜噜噜噜久久久久久91| 欧美色欧美亚洲另类二区| 网址你懂的国产日韩在线| 草草在线视频免费看| 禁无遮挡网站| 九色国产91popny在线| 日韩欧美 国产精品| 最近视频中文字幕2019在线8| 亚洲欧美激情综合另类| 特级一级黄色大片| 无遮挡黄片免费观看| 免费人成视频x8x8入口观看| 欧美国产日韩亚洲一区| 日本黄大片高清| 国产一区二区三区在线臀色熟女| 在线免费观看不下载黄p国产 | 日本 欧美在线| 国产aⅴ精品一区二区三区波| 日本免费一区二区三区高清不卡| 国产黄a三级三级三级人| 少妇人妻精品综合一区二区 | 午夜福利18| 最近最新免费中文字幕在线| 国产又黄又爽又无遮挡在线| 1000部很黄的大片| 久久久久久久精品吃奶| 小蜜桃在线观看免费完整版高清| 亚洲精品日韩av片在线观看 | 国产成人欧美在线观看| 757午夜福利合集在线观看| 精品乱码久久久久久99久播| 欧美色视频一区免费| 91九色精品人成在线观看| 3wmmmm亚洲av在线观看| 欧美在线一区亚洲| 日本一本二区三区精品| 午夜免费成人在线视频| 夜夜爽天天搞| av专区在线播放| 村上凉子中文字幕在线| 欧美成人一区二区免费高清观看| 窝窝影院91人妻| 国产一区在线观看成人免费| 亚洲成av人片免费观看| 中文资源天堂在线| 搡老岳熟女国产| 最近最新免费中文字幕在线| 亚洲av成人精品一区久久| 中文字幕av成人在线电影| 日本a在线网址| xxx96com| 97碰自拍视频| 男女之事视频高清在线观看| 亚洲欧美精品综合久久99| 变态另类成人亚洲欧美熟女| 97超视频在线观看视频| 国产精品嫩草影院av在线观看 | 久久国产精品影院| 男女床上黄色一级片免费看| 中文字幕精品亚洲无线码一区| 欧美乱色亚洲激情| 免费看美女性在线毛片视频| 欧美另类亚洲清纯唯美| 成人特级黄色片久久久久久久| 19禁男女啪啪无遮挡网站| 欧美黄色淫秽网站| 欧美日韩黄片免| 日本a在线网址| 国产又黄又爽又无遮挡在线| 国产精品影院久久| 亚洲五月婷婷丁香| 757午夜福利合集在线观看| 天堂动漫精品| 国产麻豆成人av免费视频| 国产乱人伦免费视频| 亚洲av日韩精品久久久久久密| 美女高潮的动态| 国产精品久久久久久久电影 | 欧美中文综合在线视频| av在线蜜桃| а√天堂www在线а√下载| 怎么达到女性高潮| 色在线成人网| 男插女下体视频免费在线播放| 亚洲av中文字字幕乱码综合| 亚洲自拍偷在线| 色精品久久人妻99蜜桃| 精品久久久久久久久久免费视频| 无遮挡黄片免费观看| 久99久视频精品免费| 亚洲电影在线观看av| 日日干狠狠操夜夜爽| 看片在线看免费视频| 99热只有精品国产| 97碰自拍视频| 女生性感内裤真人,穿戴方法视频| 亚洲久久久久久中文字幕| 在线观看av片永久免费下载| 久久亚洲真实| 久久精品国产综合久久久| 免费av毛片视频| 欧美中文日本在线观看视频| 亚洲av成人不卡在线观看播放网| 一个人观看的视频www高清免费观看| 看片在线看免费视频| 久久婷婷人人爽人人干人人爱| 欧美3d第一页| 欧美日韩黄片免| 男女视频在线观看网站免费| 午夜免费男女啪啪视频观看 | 欧美av亚洲av综合av国产av| 真人一进一出gif抽搐免费| 在线十欧美十亚洲十日本专区| 麻豆国产97在线/欧美| 久久伊人香网站| 欧美+亚洲+日韩+国产| 俺也久久电影网| 久久精品国产亚洲av涩爱 | 亚洲国产精品久久男人天堂| 一进一出抽搐gif免费好疼| 99国产精品一区二区蜜桃av| 国产探花在线观看一区二区| 国产不卡一卡二| 国产中年淑女户外野战色| 国产亚洲精品久久久久久毛片| 18禁美女被吸乳视频| 中文字幕av在线有码专区| 草草在线视频免费看| 亚洲人与动物交配视频| 99视频精品全部免费 在线| 国产私拍福利视频在线观看| 久久国产乱子伦精品免费另类| 18禁裸乳无遮挡免费网站照片| 一本综合久久免费| 五月玫瑰六月丁香| 国产成人啪精品午夜网站| 一a级毛片在线观看| 97超级碰碰碰精品色视频在线观看| 日本黄大片高清| 99riav亚洲国产免费| 成人国产综合亚洲| 少妇裸体淫交视频免费看高清| 国产欧美日韩精品一区二区| 国产精品久久久久久精品电影| 特大巨黑吊av在线直播| 欧美日本亚洲视频在线播放| 欧美日韩一级在线毛片| 亚洲人成网站在线播放欧美日韩| 69av精品久久久久久| 韩国av一区二区三区四区| 国产黄a三级三级三级人| 国产综合懂色| 3wmmmm亚洲av在线观看| 桃色一区二区三区在线观看| 日本黄大片高清| 69人妻影院| netflix在线观看网站| 欧美性猛交黑人性爽| 日日摸夜夜添夜夜添小说| 91麻豆精品激情在线观看国产| 欧美日韩亚洲国产一区二区在线观看| 成人国产一区最新在线观看| 给我免费播放毛片高清在线观看| 免费无遮挡裸体视频| 女人十人毛片免费观看3o分钟| 国产探花极品一区二区| 亚洲精品美女久久久久99蜜臀| 国产精品永久免费网站| 欧美bdsm另类| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品日韩av片在线观看 | svipshipincom国产片| 日本与韩国留学比较| 亚洲最大成人手机在线| 久久国产精品影院| 国产精品 国内视频| 国内精品一区二区在线观看| 午夜福利免费观看在线| 国产在线精品亚洲第一网站| 婷婷六月久久综合丁香| 欧美一级毛片孕妇| 亚洲真实伦在线观看| 国产真实乱freesex| 在线a可以看的网站| 日本 欧美在线| 老熟妇乱子伦视频在线观看| 一个人免费在线观看电影| 亚洲无线观看免费| 99国产综合亚洲精品| 亚洲午夜理论影院| 午夜福利免费观看在线| netflix在线观看网站| 久久6这里有精品| av天堂中文字幕网| 欧美成人a在线观看| 国产精品一区二区三区四区免费观看 | 国产视频内射| 欧美又色又爽又黄视频| 国产日本99.免费观看| 色在线成人网| 婷婷丁香在线五月| 一区福利在线观看| 日本三级黄在线观看| 久久久久久久久大av| 亚洲天堂国产精品一区在线| 国产免费av片在线观看野外av| 欧美极品一区二区三区四区| 日韩欧美精品v在线| 欧美3d第一页| 国产午夜精品久久久久久一区二区三区 | 午夜免费成人在线视频| x7x7x7水蜜桃| 欧美又色又爽又黄视频| 国产亚洲欧美98| 午夜福利18| 久久久国产成人精品二区| 久久国产乱子伦精品免费另类| 精品午夜福利视频在线观看一区| 久久国产乱子伦精品免费另类| 欧美成人免费av一区二区三区| avwww免费| 亚洲天堂国产精品一区在线| 亚洲18禁久久av| 黄色丝袜av网址大全| 无遮挡黄片免费观看| 欧美中文综合在线视频| 亚洲激情在线av| 色av中文字幕| 久久精品国产亚洲av香蕉五月| 天堂影院成人在线观看| 国产精品免费一区二区三区在线| 精品一区二区三区av网在线观看| 国产av不卡久久| 国产av麻豆久久久久久久| 淫秽高清视频在线观看| 成人永久免费在线观看视频| 香蕉丝袜av| 日本熟妇午夜| 日韩av在线大香蕉| 欧美最新免费一区二区三区 | 亚洲国产精品999在线| 久久精品国产99精品国产亚洲性色| 亚洲精品日韩av片在线观看 | 日本成人三级电影网站| 欧美日韩一级在线毛片| 亚洲精品日韩av片在线观看 | 国产伦在线观看视频一区| 一夜夜www| 白带黄色成豆腐渣| www.www免费av|