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

    避風(fēng)型漁港規(guī)劃問(wèn)題的啟發(fā)式算法研究

    2012-08-10 12:21:15于紅馮艷紅李放孫庚欒曙光
    關(guān)鍵詞:進(jìn)港漁港漁船

    于紅馮艷紅李放孫庚欒曙光

    (1.大連海洋大學(xué)信息工程學(xué)院,遼寧大連116023;2.大連海洋大學(xué)海洋與土木工程學(xué)院,遼寧大連116023)

    避風(fēng)型漁港規(guī)劃問(wèn)題的啟發(fā)式算法研究

    于紅1,馮艷紅1,李放2,孫庚1,欒曙光2

    (1.大連海洋大學(xué)信息工程學(xué)院,遼寧大連116023;2.大連海洋大學(xué)海洋與土木工程學(xué)院,遼寧大連116023)

    避風(fēng)型漁港的合理規(guī)劃對(duì)減少臺(tái)風(fēng)造成漁業(yè)的損失具有重要意義。建立了中國(guó)東南沿海避風(fēng)型漁港規(guī)劃問(wèn)題的數(shù)學(xué)模型,并對(duì)該模型的復(fù)雜度進(jìn)行了分析,利用漁船就近避風(fēng)的原則,提出了一種優(yōu)先考慮最短回港時(shí)間的啟發(fā)式算法,并用中國(guó)東南沿海實(shí)際的漁港及漁船數(shù)據(jù)進(jìn)行了計(jì)算,用調(diào)研數(shù)據(jù)對(duì)算法的運(yùn)行效果進(jìn)行了評(píng)估。試驗(yàn)結(jié)果表明,本算法具有較好的準(zhǔn)確度。

    避風(fēng)型漁港;漁港規(guī)劃;啟發(fā)式算法

    臺(tái)風(fēng)和強(qiáng)風(fēng)是影響漁業(yè)生產(chǎn)安全的重要因素之一,中國(guó)東南沿海是受臺(tái)風(fēng)和強(qiáng)風(fēng)影響最嚴(yán)重的地區(qū),1990—2010年,平均每年約4個(gè)臺(tái)風(fēng)登陸中國(guó)東南沿海,其中1個(gè)屬于破壞力較強(qiáng)的強(qiáng)臺(tái)風(fēng)及超強(qiáng)臺(tái)風(fēng),給國(guó)家的漁業(yè)生產(chǎn)造成了巨大的損失。為了降低臺(tái)風(fēng)和強(qiáng)風(fēng)帶來(lái)的損失,漁船需要在臺(tái)風(fēng)和強(qiáng)風(fēng)到來(lái)前找到合理的避風(fēng)場(chǎng)所,目前中國(guó)東南五省的漁港或避風(fēng)錨地可以為45%的漁船提供安全避風(fēng)場(chǎng)所[1],覆蓋面偏低。農(nóng)業(yè)部漁業(yè)局在“十二五”漁港建設(shè)與發(fā)展規(guī)劃中,將漁港防災(zāi)減災(zāi)作為漁業(yè)可持續(xù)發(fā)展的重中之重,擬新建、擴(kuò)建一批避風(fēng)型漁港和避風(fēng)錨地,以保證為70%的漁船提供安全避風(fēng)場(chǎng)所??茖W(xué)地規(guī)劃漁港布局,使盡可能多的漁船能夠在最短的時(shí)間內(nèi)到達(dá)港內(nèi)避風(fēng),最大限度地減少臺(tái)風(fēng)和強(qiáng)風(fēng)給漁業(yè)生產(chǎn)帶來(lái)的損失,是目前亟待解決的問(wèn)題。

    1 問(wèn)題描述與避風(fēng)型漁港的規(guī)劃模型

    1.1 避風(fēng)型漁港規(guī)劃的主要因素

    臺(tái)風(fēng)和強(qiáng)風(fēng)到來(lái)前,漁船一般分布在漁場(chǎng)中捕魚作業(yè),出于安全和經(jīng)濟(jì)的原因,在接到臺(tái)風(fēng)和強(qiáng)風(fēng)預(yù)報(bào)時(shí),漁船會(huì)選擇離自己最近的漁港避風(fēng)。避風(fēng)型漁港的規(guī)劃問(wèn)題描述為:有n個(gè)位置可以建設(shè)漁港 (包括已建漁港和備選港址),總共有m條漁船分布在各個(gè)漁場(chǎng),要從n個(gè)可建漁港中選擇n′個(gè)漁港,使得建設(shè)這n′個(gè)漁港之后,所有m條漁船中b%(b為常數(shù))的漁船到達(dá)建設(shè)的n′個(gè)漁港中避風(fēng),總避風(fēng)成本最小。

    避風(fēng)型漁港規(guī)劃要考慮的因素包括:漁港的位置、漁港的容量、漁船的位置、漁船的速度、漁船的自持力等,避風(fēng)型漁港規(guī)劃問(wèn)題的總目標(biāo)就是能夠讓盡可能多的漁船以最小的成本在最短的時(shí)間內(nèi)回港避風(fēng),即避風(fēng)漁船的總航行距離或航行時(shí)間最短。

    1.2 避風(fēng)型漁港布局的數(shù)學(xué)模型

    設(shè)xj(j=1,2,…,n)表示漁港pj(j=1, 2,…,n)是否被選中,yij(i=1,2,…,m;j= 1,2,…,n)表示漁船si(i=1,2,…,m)是否進(jìn)漁港pj(j=1,2,…,n)避風(fēng),tij(i=1,2,…,m;j=1,2,…,n)表示漁船si從其作業(yè)位置到達(dá)漁港pj所用的時(shí)間,漁港pj的容量為cpj,避風(fēng)型漁港規(guī)劃問(wèn)題的目標(biāo)是避風(fēng)漁船總航行時(shí)間最短,因?yàn)槊織l漁船最終只能進(jìn)入一個(gè)漁港避風(fēng),漁港pj只能為cpj條漁船避風(fēng)。因此,目標(biāo)函數(shù)及約束條件如下:

    2 漁港規(guī)劃的啟發(fā)式算法

    2.1 問(wèn)題及背景分析

    按照問(wèn)題的描述,本算法是在已知漁船數(shù)量、每條漁船位置、每條漁船航速、每條漁船類型、漁港數(shù)量 (包括已建和未建)、每個(gè)漁港位置、每個(gè)漁港類型、每個(gè)漁港容量的基礎(chǔ)上,根據(jù)漁船避風(fēng)的實(shí)際需求,計(jì)算在保證70%的漁船避風(fēng)的情況下,建設(shè)哪些漁港成本最小。從問(wèn)題的實(shí)際需求和所建立的數(shù)學(xué)模型可以看出,避風(fēng)型漁港合理規(guī)劃問(wèn)題屬于二級(jí)0-1規(guī)劃問(wèn)題,漁港選擇是一個(gè)0-1規(guī)劃問(wèn)題,規(guī)劃的目標(biāo)函數(shù)為所給出的漁港選擇方案的漁船進(jìn)港避風(fēng)成本最小,而對(duì)于任何一種漁港選擇方案,其漁船進(jìn)港避風(fēng)的成本計(jì)算均為一個(gè)0 -1背包問(wèn)題,因此,本研究中建立的模型是一個(gè)多級(jí)組合優(yōu)化問(wèn)題,屬于復(fù)雜的組合優(yōu)化問(wèn)題。其中涉及到兩組決策變量:漁港的選擇向量X和漁船進(jìn)港避風(fēng)矩陣Y,模型的最終目標(biāo)是求出一組使得整個(gè)目標(biāo)函數(shù)值最小的決策變量X;而計(jì)算決策變量X的代價(jià)需要求解使得該方案的目標(biāo)函數(shù)值最小的決策變量Y。組合優(yōu)化問(wèn)題理論上已經(jīng)被證明屬于NP-Hard問(wèn)題,到目前為止,找不到求解這類問(wèn)題最優(yōu)解的多項(xiàng)式時(shí)間算法,而解決這類問(wèn)題最有效的算法是啟發(fā)式算法 (Heuristic Algorithm)[2-4],比較典型的算法策略包括蟻群算法[2-3]、遺傳算法[4]、禁忌搜索算法[5]和模擬退火算法[6]等。為此,本研究中作者利用問(wèn)題的背景知識(shí),提出了解決該問(wèn)題的啟發(fā)式算法。

    確定求解該問(wèn)題的算法,首先要考慮漁船避風(fēng)的實(shí)際運(yùn)行情況。避風(fēng)型漁港規(guī)劃建設(shè)的基本原則是在臺(tái)風(fēng)來(lái)臨時(shí)盡可能減少臺(tái)風(fēng)所帶來(lái)的損失,讓盡可能多的漁船能夠進(jìn)入離自己作業(yè)區(qū)最近的漁港避風(fēng)。由于實(shí)際運(yùn)行過(guò)程中采用就近避風(fēng)的原則,因此,在建立避風(fēng)型漁港規(guī)劃模型時(shí)也模擬臺(tái)風(fēng)到來(lái)時(shí)的實(shí)際運(yùn)行模式,即漁船總是選擇離自己最近的漁港避風(fēng)。

    2.2 算法的基本思想

    本研究中采用啟發(fā)式搜索策略進(jìn)行解空間的搜索。算法的基本思想是:根據(jù)漁船和漁港的位置,計(jì)算每條漁船到達(dá)其自持力范圍內(nèi)每個(gè)漁港的時(shí)間,按照航行時(shí)間 (每個(gè)航行時(shí)間關(guān)聯(lián)一個(gè)漁船和一個(gè)漁港)由小到大排序,構(gòu)成航行序列 (航行時(shí)間、漁船、漁港);依次從航行序列中取出漁船和漁港,如果漁船未找到合適的漁港且被選擇的漁港有空位,則該船進(jìn)入該港;否則從航行序列中找下一組漁船和漁港,直到70%的漁船找到避風(fēng)港;此時(shí)查看所有漁港的進(jìn)船情況,如果備選漁港的進(jìn)船數(shù)超過(guò)容量的50%,就認(rèn)為在該地區(qū)周圍需要建設(shè)避風(fēng)漁港。

    在實(shí)際運(yùn)行過(guò)程中,中型船和大型船需要到一級(jí)漁港和中心漁港避風(fēng),而小型船可以到一級(jí)以下的漁港避風(fēng)。如果完全按照就近避風(fēng)原則,小型船離港較近,航行時(shí)間較短,中型船和大型船離港較遠(yuǎn),航行時(shí)間較長(zhǎng),小型船將會(huì)先到達(dá)漁港,進(jìn)而占滿中心漁港和一級(jí)漁港,導(dǎo)致最需要進(jìn)港避風(fēng)的中型船和大型船無(wú)處避風(fēng)。因此,本算法考慮漁船進(jìn)港避風(fēng)的實(shí)際情況,在中心漁港和一級(jí)漁港的規(guī)劃過(guò)程中優(yōu)先考慮中型船和大型船,在中型船和大型船均已進(jìn)港避風(fēng)的情況下,再允許小型船進(jìn)入中心漁港和一級(jí)漁港避風(fēng)。

    2.3 算法描述

    簡(jiǎn)單的漁港規(guī)劃啟發(fā)式算法Algorithm NHPPA描述如下。

    輸入:漁船集合S={s1,s2,…,sm};漁港集合P={p1,p2,…,pn};航行時(shí)間集合T= {t11,t12,…,tij,…tmn},tij(i=1,2,…,m;j=1,2,…,n)表示漁船si(i=1,2,…,m)從其作業(yè)位置到達(dá)漁港pj(j=1,2,…,n)所用的時(shí)間。

    輸出:新增的漁港集合NPS。

    變量說(shuō)明:SA和SB表示漁船集合,漁船si的主要參數(shù)為類型type和是否進(jìn)港Flag,進(jìn)港標(biāo)志Flag的初值均為false;PA和PB表示漁港集合,漁港pj的主要參數(shù)為容量c和進(jìn)船數(shù)num,進(jìn)船數(shù)num的初值均為0;NS1和NS2分別表示大中型漁船的航行序列集合和小型漁船的航行序列集合,每個(gè)航行ns包括航行時(shí)間t、漁船s、漁港p3個(gè)參數(shù)。

    算法描述如下:

    (1)Set Navigation Set NS1 and NS2 to empty

    (2)index1=index2=1

    (3)For each siin S{//按照船型對(duì)航行序列進(jìn)行分類

    (4) If si.Type="大型"or si.Type="中型"then{

    (5) For each pjin P{

    (6) If pj.Type="中心"or pj.Type="一級(jí)" then{

    (7) nsindex1.t=tij:nsindex1.p=pj:nsindex1.s=si

    (8) NS1=NS1∪{nsindex1}

    (9) index1++

    (10) }

    (11) }

    (12) }

    (13) Else{

    (14) For each pjin P{

    (15) nsindex2.t=tij:nsindex2.p=pj:nsindex2.s=si

    (16) NS2=NS2∪{nsindex2}

    (17) index2++

    (18) }

    (19) }

    (20)}

    (21)Sort NS1 and NS2 //按照航行時(shí)間由小到大對(duì)每個(gè)航行序列排序

    (22)shipCount=0

    (23)For each nsiin NS1{//讓中型船和大型船優(yōu)先進(jìn)港

    (24) If nsi.s.flag=false then{

    (25) If nsi.p.num<nsi.p.c then{

    (26) nsi.s.flag=true:nsi.p.num++

    (27) shipCount++

    (28) If shipCount>=α*m then break//α是漁船進(jìn)港比例系數(shù)

    (29) }

    (30) }

    (31)}

    (32)For each nsiin NS2{//讓小型船進(jìn)港

    (33) If nsi.s.flag=false then{

    (34) If nsi.p.num<nsi.p.c then{

    (35) nsi.s.flag=true:nsi.p.num++

    (36) shipCount++

    (37) If shipCount>=α*m then break;// α是漁船進(jìn)港比例系數(shù)

    (38) }

    (39) }

    (40)}

    (41)Set New Port Set NPS to empty

    (42)For each pjin P{

    (43) If pj.num>=pj.c*β{//β是符合建港的容量系數(shù)

    (44) NPS=NPS∪{pj}

    (45) }

    (46)}

    (47)輸出NPS

    2.4 算法的改進(jìn)

    在進(jìn)行初步的算法設(shè)計(jì)之后,通過(guò)觀察和試驗(yàn)發(fā)現(xiàn),由于中國(guó)東南沿海海岸線較長(zhǎng),考慮航行成本和航行時(shí)間等因素,一般漁船會(huì)到離作業(yè)位置較近的漁港避風(fēng),不會(huì)到遠(yuǎn)離自己的漁港避風(fēng),因此不需要計(jì)算每一條漁船到達(dá)每個(gè)漁港的時(shí)間,為了減少計(jì)算量,縮小搜索空間,將漁船分布區(qū)域和漁港分布區(qū)域劃分為5個(gè)對(duì)應(yīng)的區(qū)域,區(qū)域示意圖如圖1所示。

    圖1 漁港及漁船分區(qū)示意圖Fig.1 M ap of fishing port and fishing-boat com part mentation

    圖1中每條漁船只能到與自己的區(qū)域?qū)?yīng)的漁港區(qū)及其左右相鄰區(qū)域避風(fēng),如:對(duì)于位于SA3區(qū)的漁船,只需要計(jì)算其到達(dá)PA2、PA3和PA4區(qū)域漁港的航行時(shí)間,而對(duì)于位于SA5區(qū)域的漁船,則只需要計(jì)算其到達(dá)PA4和PA5區(qū)域漁港的航行時(shí)間。

    進(jìn)行區(qū)域劃分以及根據(jù)區(qū)域確定漁船的避風(fēng)漁港,可以大大減少漁港規(guī)劃問(wèn)題的時(shí)間和空間開銷,從而提高計(jì)算速度。

    3 實(shí)驗(yàn)

    3.1 實(shí)驗(yàn)設(shè)計(jì)

    為了檢驗(yàn)算法的效果設(shè)計(jì)一組實(shí)驗(yàn),實(shí)驗(yàn)環(huán)境采用 ThinkPad(Intel I3,CPU 2.13 GHz,RAM 8 GB,64位Windows 7),用Java實(shí)現(xiàn)該算法。實(shí)驗(yàn)中漁船數(shù)據(jù)參照中國(guó)東南沿海實(shí)際擁有的漁船數(shù);漁船位置主要根據(jù)中國(guó)東南沿海漁業(yè)資源分布密度,使?jié)O船在漁場(chǎng)中按密度均勻分布;一級(jí)漁港和中心漁港的位置和容量用中國(guó)東南沿海一級(jí)和中心漁港的實(shí)際數(shù)據(jù);一級(jí)以下漁港沿海岸線均勻分布,另外在中國(guó)東南沿海岸邊選擇了一百多個(gè)天然避風(fēng)條件較好的岙口,作為擬建避風(fēng)型漁港的港址。漁港和漁船區(qū)域參照漁場(chǎng)分布以及各省的海岸線進(jìn)行劃分。

    3.2 實(shí)驗(yàn)結(jié)果的評(píng)價(jià)

    為了避免偶然因素對(duì)算法運(yùn)行結(jié)果的影響,將算法運(yùn)行5次,用5次運(yùn)行算出的擬建港址交集作為最后的擬建港址集合,用5次的平均運(yùn)算時(shí)間作為算法的運(yùn)行時(shí)間,算法的運(yùn)行時(shí)間為5 min。

    算法運(yùn)行前,根據(jù)地理位置的特點(diǎn)選擇一些位置作為候選擬建港址;算法運(yùn)行時(shí),根據(jù)作業(yè)漁船的分布輸出一個(gè)擬建漁港集合。為了評(píng)價(jià)算法的效果,作者利用在浙江省進(jìn)行漁港調(diào)研的機(jī)會(huì),對(duì)計(jì)算結(jié)果中12個(gè)擬建港址進(jìn)行調(diào)研 (受調(diào)研經(jīng)費(fèi)的限制,不能對(duì)所有港址進(jìn)行調(diào)研),了解這12個(gè)漁港與實(shí)際建港需求的符合程度。調(diào)研時(shí),通過(guò)現(xiàn)場(chǎng)考察和與當(dāng)?shù)貪O業(yè)部門領(lǐng)導(dǎo)及漁民座談等形式,了解當(dāng)?shù)乇茱L(fēng)漁港建設(shè)的實(shí)際需求。結(jié)果了解到在多處急需建港的地方,浙江省正在進(jìn)行漁港的改擴(kuò)建,因此,本研究中決定以這12個(gè)擬建港址與實(shí)際改擴(kuò)建漁港的符合程度評(píng)價(jià)模型和算法的效果??紤]到確定候選擬建漁港時(shí)并不了解實(shí)際的漁港改擴(kuò)建情況,候選擬建漁港港址和實(shí)際改擴(kuò)建漁港港址之間存在一定誤差,而某個(gè)擬建港址周圍一定距離內(nèi)有在建漁港,說(shuō)明該擬建漁港周圍有建港需求,因此,用擬建港址與周邊正在改擴(kuò)建的漁港之間的距離來(lái)描述擬建漁港與在建漁港的符合程度,調(diào)研結(jié)果表明,擬建港址與周邊正在改擴(kuò)建漁港的最近距離分別為 6.5、 7.0、 7.0、 2.8、 0.4、11.0、2.0、5.0、0、20.0、6.5、13.0 km。征求當(dāng)?shù)貪O業(yè)管理部門的意見后,認(rèn)為二者之間距離小于10 km可以認(rèn)為二者相符。據(jù)此,本算法的準(zhǔn)確率達(dá)到75%。

    4 結(jié)語(yǔ)

    為了減少臺(tái)風(fēng)及強(qiáng)臺(tái)風(fēng)給中國(guó)漁業(yè)生產(chǎn)造成的影響,首次對(duì)避風(fēng)型漁港和避風(fēng)錨地合理布局問(wèn)題進(jìn)行了研究,建立了中國(guó)東南沿海避風(fēng)型漁港規(guī)劃問(wèn)題的數(shù)學(xué)模型,在此基礎(chǔ)上提出了一種優(yōu)先考慮最短回港時(shí)間的啟發(fā)式算法,用中國(guó)東南沿海實(shí)際的漁港及漁船數(shù)據(jù)驗(yàn)證了該模型和算法的有效性,用現(xiàn)場(chǎng)調(diào)研數(shù)據(jù)對(duì)算法的運(yùn)行結(jié)果進(jìn)行了評(píng)估。結(jié)果表明,本研究中建立的模型和算法取得了較好的效果。但是由于計(jì)算結(jié)果存在一定的誤差,因此不能直接應(yīng)用于實(shí)際漁港布局規(guī)劃,下一步的工作需要根據(jù)實(shí)際應(yīng)用背景,綜合考慮休漁期和漁業(yè)資源分布不均等因素,對(duì)模型和算法進(jìn)行修正,進(jìn)一步改進(jìn)模型和算法的效果。

    [1] 農(nóng)業(yè)部漁業(yè)局.2010中國(guó)漁業(yè)統(tǒng)計(jì)年鑒[M].北京:中國(guó)農(nóng)業(yè)出版社,2010.

    [2] 熊偉清,魏平.二進(jìn)制蟻群進(jìn)化算法[J].自動(dòng)化學(xué)報(bào),2007,33 (3):259-264.

    [3] 袁軍良,熊偉清,江寶釧.混合二元蟻群算法求解集裝箱裝載問(wèn)題[J].計(jì)算機(jī)工程與應(yīng)用,2010,46(36):222-225.

    [4] 鐘金宏,黃玲.帶外包受限批量模型的啟發(fā)式遺傳算法[J].系統(tǒng)仿真學(xué)報(bào),2011,23(12):2623-2628.

    [5] Fred Glover.Tabu Search—Part I[J].ORSA Journal on Computing,1989,1(3):190-206.

    [6] Kirkpatrick S,Gelatt Jr C D,VecchiM P.Optimization by simulated annealing[J].Science,1983,220:671-680.

    Heuristic algorithm for planning problem of fishing port sheltered from typhoon

    YU Hong1,FENG Yan-hong1,LIFang2,SUN Geng1,LUAN Shu-guang2
    (1.College of Information Engineering,Dalian Ocean University,Dalian 116023,China;2.College of Marine and Civil Engineering, Dalian Ocean University,Dalian 116023,China)

    The rational planning of fishing port sheltered from typhoon plays an important role in reducing fishery loss from typhoon.Themathematic model of fishing port planning in the southeastern coastal region is built.The complexity of themodel is analyzed.Following the principle of“Fishing vessles should shelter from typhoon in the nearest fishing port”,the paper puts forward a heuristic algorithm,which firstly considers the shortest time to return fishing port.The actual data of fishing ports and vessels in the southeast region are utilized in calculation and the survey data are used to envaluate computing result of heuristic algorithm.The result shows that the heuristic algorithm carries accepable accuracy.

    fishing port sheltered from typhoon;planning of fishing port;heuristic algorithm

    TP311

    A

    2012-05-28

    農(nóng)業(yè)部項(xiàng)目 (農(nóng)財(cái)發(fā) [2012]26號(hào))

    于紅 (1968-),女,教授。E-mail:yuhong@dlou.edu.cn

    2095-1388(2012)04-0373-04

    猜你喜歡
    進(jìn)港漁港漁船
    漁船
    千舟競(jìng)發(fā)
    廉政瞭望(2021年15期)2021-08-23 04:55:22
    大型滿載油輪使用魚山作業(yè)區(qū)南部進(jìn)港航道航行方法探討
    航海(2020年6期)2020-12-23 06:53:37
    漁港
    船舶進(jìn)靠浙能臺(tái)二電煤炭碼頭風(fēng)險(xiǎn)的研究
    開漁后的博賀漁港總是忙碌而又充滿生機(jī)
    國(guó)內(nèi)新型遠(yuǎn)洋金槍魚圍網(wǎng)漁船首航
    商周刊(2018年19期)2018-12-06 09:49:41
    相聚在王浩兒漁港
    黃河之聲(2017年1期)2017-04-27 13:21:56
    漁船驚魂
    靜靜的漁港 遠(yuǎn)航的風(fēng)帆
    久久精品影院6| 最后的刺客免费高清国语| 俺也久久电影网| 色综合欧美亚洲国产小说| 99热这里只有是精品50| 亚洲av五月六月丁香网| 精品人妻偷拍中文字幕| 国产真实伦视频高清在线观看 | 午夜老司机福利剧场| 一级av片app| 9191精品国产免费久久| av在线观看视频网站免费| 久久精品影院6| 男插女下体视频免费在线播放| 亚洲美女搞黄在线观看 | 99精品在免费线老司机午夜| 精品人妻偷拍中文字幕| 长腿黑丝高跟| 麻豆成人午夜福利视频| 精品国内亚洲2022精品成人| 日韩免费av在线播放| 男女床上黄色一级片免费看| 好男人在线观看高清免费视频| 国产精品永久免费网站| 成人特级黄色片久久久久久久| 亚洲成人中文字幕在线播放| 色播亚洲综合网| 国产精品久久久久久久久免 | 一二三四社区在线视频社区8| 在线观看午夜福利视频| 男女那种视频在线观看| 两个人的视频大全免费| 亚洲自偷自拍三级| 97热精品久久久久久| 99热这里只有精品一区| 国产精品亚洲一级av第二区| 国产精品一区二区三区四区久久| 可以在线观看的亚洲视频| 亚洲三级黄色毛片| 999久久久精品免费观看国产| 色噜噜av男人的天堂激情| 日韩成人在线观看一区二区三区| 高清日韩中文字幕在线| 亚洲av第一区精品v没综合| bbb黄色大片| 99热这里只有精品一区| 夜夜躁狠狠躁天天躁| 白带黄色成豆腐渣| 国产私拍福利视频在线观看| 亚洲欧美日韩高清专用| 久久久精品大字幕| 男人舔奶头视频| 国产主播在线观看一区二区| 国产三级中文精品| 特大巨黑吊av在线直播| 婷婷丁香在线五月| 99精品在免费线老司机午夜| 欧美不卡视频在线免费观看| 亚洲精品亚洲一区二区| 亚州av有码| 日本a在线网址| 少妇人妻精品综合一区二区 | 女生性感内裤真人,穿戴方法视频| 国产一区二区在线av高清观看| 午夜免费激情av| 国产aⅴ精品一区二区三区波| 在线观看66精品国产| 男女做爰动态图高潮gif福利片| 亚洲国产精品成人综合色| 老司机福利观看| 欧美精品啪啪一区二区三区| 他把我摸到了高潮在线观看| 桃红色精品国产亚洲av| 免费在线观看日本一区| 国产男靠女视频免费网站| 欧美性感艳星| 内地一区二区视频在线| 一区二区三区激情视频| 亚洲熟妇熟女久久| aaaaa片日本免费| 亚洲国产日韩欧美精品在线观看| 国产一区二区在线观看日韩| 成人午夜高清在线视频| 国产aⅴ精品一区二区三区波| 亚洲 欧美 日韩 在线 免费| 久久久久久久午夜电影| 国产精品亚洲av一区麻豆| 熟妇人妻久久中文字幕3abv| 国产蜜桃级精品一区二区三区| 一区二区三区免费毛片| 欧美潮喷喷水| 亚洲欧美日韩无卡精品| 欧美xxxx黑人xx丫x性爽| 亚洲无线观看免费| 99国产极品粉嫩在线观看| 国产三级黄色录像| 两人在一起打扑克的视频| 国产精品亚洲美女久久久| 婷婷色综合大香蕉| 十八禁人妻一区二区| av在线老鸭窝| 亚洲真实伦在线观看| 一级毛片久久久久久久久女| 直男gayav资源| 少妇高潮的动态图| 日韩欧美三级三区| 很黄的视频免费| 国产成+人综合+亚洲专区| 亚洲国产精品999在线| 久久久久久久午夜电影| 日本一二三区视频观看| 亚洲第一电影网av| 日韩免费av在线播放| av在线天堂中文字幕| 国产精品伦人一区二区| 久久久久国产精品人妻aⅴ院| 在线观看一区二区三区| 在线观看午夜福利视频| 看黄色毛片网站| 国产三级在线视频| 成人无遮挡网站| 九色国产91popny在线| 免费在线观看影片大全网站| 观看美女的网站| 欧美高清性xxxxhd video| 精品不卡国产一区二区三区| 亚洲无线在线观看| 亚洲精品乱码久久久v下载方式| 在线天堂最新版资源| 国产av不卡久久| 能在线免费观看的黄片| 91狼人影院| 在线观看一区二区三区| 男人的好看免费观看在线视频| 国内精品久久久久精免费| 大型黄色视频在线免费观看| 欧美在线一区亚洲| 久久久久久久亚洲中文字幕 | 久久草成人影院| 婷婷精品国产亚洲av| 亚州av有码| 非洲黑人性xxxx精品又粗又长| 久久精品国产亚洲av涩爱 | 国产色婷婷99| 国产男靠女视频免费网站| 女人十人毛片免费观看3o分钟| 天堂影院成人在线观看| 欧美又色又爽又黄视频| 国产成人影院久久av| 国产精品电影一区二区三区| 亚洲国产精品sss在线观看| 欧美黑人巨大hd| 精品久久久久久久久久免费视频| 日韩成人在线观看一区二区三区| 午夜精品在线福利| 少妇熟女aⅴ在线视频| 香蕉av资源在线| 国产真实乱freesex| 中文资源天堂在线| 国产精品免费一区二区三区在线| 热99re8久久精品国产| 天天一区二区日本电影三级| 少妇的逼好多水| 97超级碰碰碰精品色视频在线观看| 黄色配什么色好看| 国产国拍精品亚洲av在线观看| 很黄的视频免费| 成年版毛片免费区| 一二三四社区在线视频社区8| 97超级碰碰碰精品色视频在线观看| 激情在线观看视频在线高清| 国产一区二区亚洲精品在线观看| 内地一区二区视频在线| 日本免费一区二区三区高清不卡| 一区二区三区免费毛片| 亚洲中文日韩欧美视频| 我的女老师完整版在线观看| 欧美丝袜亚洲另类 | 国产主播在线观看一区二区| 亚洲成人久久爱视频| 亚洲国产精品合色在线| 综合色av麻豆| 观看免费一级毛片| 国产野战对白在线观看| 一本一本综合久久| 国产精品98久久久久久宅男小说| 女生性感内裤真人,穿戴方法视频| 国产v大片淫在线免费观看| 全区人妻精品视频| 亚洲七黄色美女视频| 又黄又爽又刺激的免费视频.| 久久欧美精品欧美久久欧美| 一本久久中文字幕| 伊人久久精品亚洲午夜| 精品乱码久久久久久99久播| 一级黄片播放器| 网址你懂的国产日韩在线| 久久精品影院6| 亚洲第一欧美日韩一区二区三区| 三级国产精品欧美在线观看| 欧美黑人巨大hd| 麻豆一二三区av精品| 午夜福利欧美成人| 一级黄片播放器| 免费av不卡在线播放| 久久久久久久午夜电影| 男女那种视频在线观看| 亚洲美女黄片视频| 国产成人aa在线观看| 看十八女毛片水多多多| 一级毛片久久久久久久久女| 亚洲人成网站高清观看| 久久国产乱子伦精品免费另类| 亚洲欧美清纯卡通| 国产精品国产高清国产av| 变态另类丝袜制服| 精品久久久久久久久久久久久| av天堂在线播放| 久久久久国内视频| 国产精品国产高清国产av| 精品久久久久久久久av| 一本精品99久久精品77| 少妇被粗大猛烈的视频| 久久久久性生活片| 九九在线视频观看精品| 成人午夜高清在线视频| 亚洲无线观看免费| 欧美一区二区国产精品久久精品| 在线播放国产精品三级| 亚洲成人中文字幕在线播放| 亚洲欧美激情综合另类| 亚洲成人久久性| 丝袜美腿在线中文| 好看av亚洲va欧美ⅴa在| 免费电影在线观看免费观看| 午夜福利成人在线免费观看| 免费看a级黄色片| 青草久久国产| 国产精品98久久久久久宅男小说| 最近中文字幕高清免费大全6 | 亚洲一区二区三区不卡视频| 琪琪午夜伦伦电影理论片6080| 别揉我奶头 嗯啊视频| 又粗又爽又猛毛片免费看| 亚洲av二区三区四区| 欧美一区二区国产精品久久精品| 淫妇啪啪啪对白视频| 搡老岳熟女国产| 别揉我奶头~嗯~啊~动态视频| 精品福利观看| 精品免费久久久久久久清纯| 脱女人内裤的视频| 亚洲av免费在线观看| 特大巨黑吊av在线直播| 午夜福利18| 日韩亚洲欧美综合| 日韩中文字幕欧美一区二区| 性色av乱码一区二区三区2| 日韩大尺度精品在线看网址| 亚洲激情在线av| 成人鲁丝片一二三区免费| 日本黄大片高清| 亚洲国产精品成人综合色| 18+在线观看网站| 日韩欧美三级三区| 国产精品久久电影中文字幕| 91av网一区二区| 狂野欧美白嫩少妇大欣赏| 午夜两性在线视频| 欧美黑人巨大hd| 国产熟女xx| 少妇的逼好多水| 在线播放国产精品三级| 国产高清激情床上av| xxxwww97欧美| 欧美+亚洲+日韩+国产| 国产亚洲av嫩草精品影院| 国内精品一区二区在线观看| 国产精品精品国产色婷婷| 首页视频小说图片口味搜索| 日本成人三级电影网站| 丝袜美腿在线中文| 脱女人内裤的视频| 国产一区二区在线观看日韩| 蜜桃亚洲精品一区二区三区| 成人国产一区最新在线观看| 亚洲人成电影免费在线| 亚洲第一区二区三区不卡| 毛片一级片免费看久久久久 | 黄色日韩在线| 亚洲欧美日韩卡通动漫| 久久精品国产亚洲av涩爱 | 伦理电影大哥的女人| 噜噜噜噜噜久久久久久91| 国内精品久久久久精免费| 日韩欧美一区二区三区在线观看| 欧洲精品卡2卡3卡4卡5卡区| avwww免费| 亚洲一区高清亚洲精品| 国产精品影院久久| 亚洲成人久久爱视频| 午夜福利在线观看免费完整高清在 | 97超级碰碰碰精品色视频在线观看| 人妻夜夜爽99麻豆av| 国模一区二区三区四区视频| 国产精品国产高清国产av| 极品教师在线免费播放| 午夜日韩欧美国产| 一个人免费在线观看的高清视频| 中文字幕av在线有码专区| 国产高清视频在线观看网站| 天堂网av新在线| 亚洲激情在线av| 蜜桃久久精品国产亚洲av| 久久国产乱子伦精品免费另类| 无人区码免费观看不卡| 国产三级中文精品| 亚洲人成伊人成综合网2020| 亚洲av成人不卡在线观看播放网| 久久午夜亚洲精品久久| 身体一侧抽搐| 亚洲最大成人手机在线| 亚洲 欧美 日韩 在线 免费| 欧美xxxx黑人xx丫x性爽| 成人国产综合亚洲| 啪啪无遮挡十八禁网站| 国产乱人伦免费视频| 51午夜福利影视在线观看| 亚洲午夜理论影院| 国产黄a三级三级三级人| 亚洲av成人精品一区久久| 欧美午夜高清在线| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 一级黄片播放器| 午夜福利高清视频| 久99久视频精品免费| 三级国产精品欧美在线观看| 国产淫片久久久久久久久 | 毛片一级片免费看久久久久 | 午夜福利在线观看免费完整高清在 | 国产午夜精品论理片| 欧美日韩中文字幕国产精品一区二区三区| 窝窝影院91人妻| 免费高清视频大片| 男女那种视频在线观看| 欧美精品国产亚洲| 国产成人a区在线观看| 我的女老师完整版在线观看| 国产色爽女视频免费观看| 国产91精品成人一区二区三区| 蜜桃久久精品国产亚洲av| 3wmmmm亚洲av在线观看| 蜜桃亚洲精品一区二区三区| 欧美bdsm另类| 色尼玛亚洲综合影院| 色噜噜av男人的天堂激情| 特大巨黑吊av在线直播| 中文字幕人成人乱码亚洲影| 国产成人啪精品午夜网站| 男女下面进入的视频免费午夜| 精品日产1卡2卡| 亚洲国产精品999在线| 国产精品亚洲av一区麻豆| 亚洲美女黄片视频| 色综合欧美亚洲国产小说| 五月伊人婷婷丁香| 久久精品国产亚洲av天美| 久久这里只有精品中国| 久久久久精品国产欧美久久久| 男女做爰动态图高潮gif福利片| 亚洲男人的天堂狠狠| 午夜a级毛片| 亚洲av不卡在线观看| 91av网一区二区| 色综合欧美亚洲国产小说| 久久精品国产自在天天线| 精品午夜福利视频在线观看一区| 亚洲成人中文字幕在线播放| 欧洲精品卡2卡3卡4卡5卡区| 我的女老师完整版在线观看| 啦啦啦观看免费观看视频高清| 亚洲熟妇中文字幕五十中出| 亚洲av中文字字幕乱码综合| av在线天堂中文字幕| 亚洲一区二区三区不卡视频| 两个人视频免费观看高清| 国产精品乱码一区二三区的特点| 嫩草影院入口| 成人特级黄色片久久久久久久| 国产亚洲精品久久久com| 动漫黄色视频在线观看| 麻豆国产av国片精品| 成人av一区二区三区在线看| x7x7x7水蜜桃| 久久国产乱子伦精品免费另类| 老司机午夜福利在线观看视频| 亚洲精品成人久久久久久| 如何舔出高潮| 黄色女人牲交| 天堂网av新在线| 变态另类丝袜制服| 欧美日本视频| 69人妻影院| 激情在线观看视频在线高清| a在线观看视频网站| 999久久久精品免费观看国产| 最近在线观看免费完整版| 两个人的视频大全免费| 免费在线观看亚洲国产| 欧美在线黄色| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 人妻久久中文字幕网| 久久亚洲精品不卡| 真实男女啪啪啪动态图| 有码 亚洲区| 午夜免费男女啪啪视频观看 | 久久精品综合一区二区三区| 91狼人影院| 日本a在线网址| 免费无遮挡裸体视频| 日韩欧美精品免费久久 | 国产美女午夜福利| 村上凉子中文字幕在线| 国产精品久久电影中文字幕| 欧美黄色片欧美黄色片| 天堂av国产一区二区熟女人妻| a级毛片免费高清观看在线播放| 亚洲无线观看免费| 夜夜爽天天搞| 亚洲 欧美 日韩 在线 免费| 日韩 亚洲 欧美在线| 国产男靠女视频免费网站| 99久国产av精品| 亚洲美女视频黄频| 成人美女网站在线观看视频| 综合色av麻豆| 黄色配什么色好看| 免费观看的影片在线观看| 12—13女人毛片做爰片一| 国产黄色小视频在线观看| 亚洲成av人片在线播放无| 制服丝袜大香蕉在线| 欧美日韩瑟瑟在线播放| 欧美成人一区二区免费高清观看| 永久网站在线| 午夜激情欧美在线| 美女大奶头视频| 久久伊人香网站| 国产精品,欧美在线| 免费人成在线观看视频色| 久久久久久国产a免费观看| 国产亚洲精品综合一区在线观看| 51午夜福利影视在线观看| 免费av不卡在线播放| 久久久色成人| av在线观看视频网站免费| 午夜福利成人在线免费观看| 精品乱码久久久久久99久播| 国产在线男女| 一区二区三区四区激情视频 | 少妇人妻精品综合一区二区 | 很黄的视频免费| 3wmmmm亚洲av在线观看| 久久人妻av系列| 国产一区二区在线av高清观看| 好看av亚洲va欧美ⅴa在| 夜夜爽天天搞| 我的老师免费观看完整版| 99久久精品国产亚洲精品| 五月伊人婷婷丁香| 十八禁国产超污无遮挡网站| 成年女人永久免费观看视频| 亚洲一区二区三区不卡视频| 欧美黑人巨大hd| 人妻夜夜爽99麻豆av| 久久99热6这里只有精品| 久久亚洲精品不卡| 免费在线观看日本一区| 免费观看人在逋| 又紧又爽又黄一区二区| 麻豆成人午夜福利视频| 国产真实伦视频高清在线观看 | av天堂中文字幕网| 国产淫片久久久久久久久 | 可以在线观看的亚洲视频| 国产午夜精品久久久久久一区二区三区 | 欧美黄色片欧美黄色片| 99久久99久久久精品蜜桃| 久久久久久久精品吃奶| 国内精品美女久久久久久| 国内精品久久久久久久电影| 久久这里只有精品中国| 可以在线观看毛片的网站| 中文字幕熟女人妻在线| 波野结衣二区三区在线| 精品熟女少妇八av免费久了| 欧美在线黄色| 人妻丰满熟妇av一区二区三区| 嫩草影院精品99| 国产乱人伦免费视频| 亚洲国产色片| 亚洲国产精品sss在线观看| 老女人水多毛片| 噜噜噜噜噜久久久久久91| 国产黄片美女视频| 午夜精品久久久久久毛片777| 久久国产乱子免费精品| 亚洲在线自拍视频| 俄罗斯特黄特色一大片| 亚洲人成网站在线播放欧美日韩| 亚洲欧美激情综合另类| 亚洲综合色惰| 日韩欧美在线乱码| 中国美女看黄片| 嫩草影视91久久| 蜜桃久久精品国产亚洲av| 亚洲熟妇熟女久久| 黄色一级大片看看| 亚洲人成电影免费在线| avwww免费| av专区在线播放| а√天堂www在线а√下载| 日本黄色视频三级网站网址| 夜夜夜夜夜久久久久| 他把我摸到了高潮在线观看| 亚洲成av人片在线播放无| 一区二区三区免费毛片| 久久久久亚洲av毛片大全| 亚洲精品乱码久久久v下载方式| 国产一区二区在线av高清观看| 黄色配什么色好看| 成人性生交大片免费视频hd| av欧美777| 美女被艹到高潮喷水动态| 亚洲欧美日韩卡通动漫| 国产精品综合久久久久久久免费| 校园春色视频在线观看| 亚洲国产欧美人成| 欧美最新免费一区二区三区 | 亚洲第一区二区三区不卡| 欧美丝袜亚洲另类 | 欧美成狂野欧美在线观看| 日韩成人在线观看一区二区三区| 天堂网av新在线| 日日干狠狠操夜夜爽| 国产高潮美女av| 精品久久久久久成人av| 琪琪午夜伦伦电影理论片6080| 看十八女毛片水多多多| 午夜激情福利司机影院| 日韩欧美一区二区三区在线观看| 在线播放无遮挡| 男插女下体视频免费在线播放| 国产午夜精品久久久久久一区二区三区 | 有码 亚洲区| 免费无遮挡裸体视频| 51国产日韩欧美| 久久精品国产亚洲av香蕉五月| 久久久久久久久中文| 精品久久久久久久末码| 亚洲精品一区av在线观看| 真实男女啪啪啪动态图| 精品免费久久久久久久清纯| 一区二区三区高清视频在线| 成人亚洲精品av一区二区| 久久久久久国产a免费观看| 免费av不卡在线播放| 97超级碰碰碰精品色视频在线观看| 宅男免费午夜| 久久精品国产亚洲av涩爱 | 91av网一区二区| 免费人成在线观看视频色| 毛片一级片免费看久久久久 | 国产精品av视频在线免费观看| 日韩欧美在线乱码| 国产精华一区二区三区| 欧美在线黄色| 赤兔流量卡办理| www.www免费av| 村上凉子中文字幕在线| 看免费av毛片| 欧美成人a在线观看| av在线观看视频网站免费| 亚洲在线自拍视频| 日韩av在线大香蕉| 成人国产一区最新在线观看| 老女人水多毛片| 最近中文字幕高清免费大全6 | 日韩中字成人| av专区在线播放| 日韩高清综合在线| 久久国产乱子免费精品| 国产激情偷乱视频一区二区| www.www免费av| 亚洲成人中文字幕在线播放| 看片在线看免费视频| 俺也久久电影网| 婷婷色综合大香蕉| 波多野结衣巨乳人妻| 桃红色精品国产亚洲av| 欧美在线黄色| 18+在线观看网站| 一本一本综合久久| 草草在线视频免费看| 韩国av一区二区三区四区| 毛片女人毛片| 国产主播在线观看一区二区| 国产熟女xx| 成人av在线播放网站| 国产精品一及| 成人高潮视频无遮挡免费网站| 亚洲成人久久爱视频| 国产精品免费一区二区三区在线| 久久久久国产精品人妻aⅴ院| 一级黄片播放器|