• <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)帆
    91狼人影院| 少妇裸体淫交视频免费看高清| a在线观看视频网站| 精品午夜福利视频在线观看一区| 国产伦一二天堂av在线观看| 精品人妻视频免费看| 九色成人免费人妻av| 精品久久久久久久末码| 丰满人妻熟妇乱又伦精品不卡| 黄色配什么色好看| 精品久久久久久,| 无遮挡黄片免费观看| 亚州av有码| 色哟哟·www| 久久久久久久精品吃奶| 人妻久久中文字幕网| 久久久久免费精品人妻一区二区| 免费看a级黄色片| 99在线人妻在线中文字幕| 俺也久久电影网| 国产精品永久免费网站| 麻豆成人av在线观看| 搡老岳熟女国产| 国产精品自产拍在线观看55亚洲| xxxwww97欧美| 中文资源天堂在线| 国产男靠女视频免费网站| 久久久久久久午夜电影| 高清毛片免费观看视频网站| 九色国产91popny在线| 又紧又爽又黄一区二区| 黄片小视频在线播放| 最近最新免费中文字幕在线| 亚洲精品乱码久久久v下载方式| 国产伦在线观看视频一区| 身体一侧抽搐| 久久久久久久久久黄片| 别揉我奶头~嗯~啊~动态视频| 国产色爽女视频免费观看| av福利片在线观看| 我的女老师完整版在线观看| 成年人黄色毛片网站| 九色国产91popny在线| 在线天堂最新版资源| 日本 av在线| 99国产极品粉嫩在线观看| 别揉我奶头 嗯啊视频| 婷婷丁香在线五月| 免费一级毛片在线播放高清视频| 亚洲国产日韩欧美精品在线观看| 午夜福利在线观看免费完整高清在 | 国产又黄又爽又无遮挡在线| 国产伦在线观看视频一区| 国产高潮美女av| 久久精品国产清高在天天线| 乱人视频在线观看| 精品熟女少妇八av免费久了| av黄色大香蕉| 真人做人爱边吃奶动态| 午夜免费男女啪啪视频观看 | 亚洲五月天丁香| 国产精品美女特级片免费视频播放器| 欧美中文日本在线观看视频| 日韩大尺度精品在线看网址| 欧美区成人在线视频| 亚洲美女视频黄频| 我要搜黄色片| 非洲黑人性xxxx精品又粗又长| aaaaa片日本免费| 欧美中文日本在线观看视频| 国内久久婷婷六月综合欲色啪| 一本精品99久久精品77| 在线国产一区二区在线| 国产黄片美女视频| 嫩草影视91久久| 久久久成人免费电影| 琪琪午夜伦伦电影理论片6080| 一个人观看的视频www高清免费观看| 熟女人妻精品中文字幕| 国产v大片淫在线免费观看| 久久久精品大字幕| 欧美激情久久久久久爽电影| 成年女人看的毛片在线观看| 亚洲天堂国产精品一区在线| 毛片女人毛片| 1024手机看黄色片| 老司机福利观看| 午夜精品一区二区三区免费看| 久久精品国产亚洲av涩爱 | 久久热精品热| eeuss影院久久| 亚洲内射少妇av| 天天躁日日操中文字幕| 精品无人区乱码1区二区| 人妻久久中文字幕网| 757午夜福利合集在线观看| 国内揄拍国产精品人妻在线| 俺也久久电影网| 天堂网av新在线| 一区二区三区高清视频在线| 超碰av人人做人人爽久久| 精品午夜福利视频在线观看一区| 网址你懂的国产日韩在线| 国产主播在线观看一区二区| 久9热在线精品视频| 免费人成在线观看视频色| 91字幕亚洲| 欧美性猛交黑人性爽| 久久99热6这里只有精品| 特级一级黄色大片| 亚洲av.av天堂| 久久久色成人| 一个人看视频在线观看www免费| 欧美又色又爽又黄视频| 午夜福利欧美成人| 丰满的人妻完整版| 长腿黑丝高跟| 一区二区三区免费毛片| 欧美在线一区亚洲| 成年免费大片在线观看| 国产 一区 欧美 日韩| 国产大屁股一区二区在线视频| 久久久久久久久中文| 亚洲精品一卡2卡三卡4卡5卡| av天堂在线播放| 十八禁网站免费在线| 色精品久久人妻99蜜桃| 亚洲一区高清亚洲精品| 免费在线观看日本一区| .国产精品久久| 久久国产乱子免费精品| 精品久久国产蜜桃| 午夜免费激情av| 色av中文字幕| 性插视频无遮挡在线免费观看| 色精品久久人妻99蜜桃| 午夜精品在线福利| 国产精品美女特级片免费视频播放器| 欧美中文日本在线观看视频| 精品国内亚洲2022精品成人| 又紧又爽又黄一区二区| 精品久久久久久久久亚洲 | 淫秽高清视频在线观看| 久久九九热精品免费| 日本三级黄在线观看| 特级一级黄色大片| 国产男靠女视频免费网站| 深爱激情五月婷婷| 美女 人体艺术 gogo| 十八禁国产超污无遮挡网站| 大型黄色视频在线免费观看| 国产高清有码在线观看视频| 国产大屁股一区二区在线视频| 国产高清视频在线播放一区| 欧美zozozo另类| 午夜久久久久精精品| 成人av一区二区三区在线看| а√天堂www在线а√下载| av黄色大香蕉| 人人妻人人澡欧美一区二区| 亚洲国产精品成人综合色| 真实男女啪啪啪动态图| 亚洲乱码一区二区免费版| 免费人成视频x8x8入口观看| a级一级毛片免费在线观看| 国产午夜精品久久久久久一区二区三区 | 少妇高潮的动态图| 亚洲内射少妇av| 久久精品国产99精品国产亚洲性色| 婷婷精品国产亚洲av| 国产欧美日韩精品亚洲av| 麻豆av噜噜一区二区三区| 激情在线观看视频在线高清| 久久久久久久午夜电影| 一本精品99久久精品77| 亚洲精品日韩av片在线观看| 久9热在线精品视频| 国产精品99久久久久久久久| 国产欧美日韩一区二区三| 欧美日韩福利视频一区二区| 亚洲专区国产一区二区| 偷拍熟女少妇极品色| 欧洲精品卡2卡3卡4卡5卡区| 婷婷精品国产亚洲av在线| 亚洲五月天丁香| 可以在线观看毛片的网站| 成人三级黄色视频| 久99久视频精品免费| 在线十欧美十亚洲十日本专区| 免费电影在线观看免费观看| 51国产日韩欧美| 伦理电影大哥的女人| 男女床上黄色一级片免费看| 校园春色视频在线观看| 中文字幕av成人在线电影| 变态另类丝袜制服| 天天躁日日操中文字幕| 99热这里只有是精品在线观看 | av国产免费在线观看| 亚洲av电影在线进入| 一进一出抽搐动态| 少妇的逼水好多| 国产爱豆传媒在线观看| 国产欧美日韩精品亚洲av| 麻豆国产av国片精品| 日韩大尺度精品在线看网址| 12—13女人毛片做爰片一| 精品午夜福利在线看| 青草久久国产| 国产高清三级在线| 动漫黄色视频在线观看| 变态另类丝袜制服| 久久久国产成人精品二区| 亚洲专区国产一区二区| 久久久久久久久久黄片| 精品熟女少妇八av免费久了| av天堂中文字幕网| 制服丝袜大香蕉在线| av在线老鸭窝| 男女下面进入的视频免费午夜| 欧美色视频一区免费| 国产精品乱码一区二三区的特点| 琪琪午夜伦伦电影理论片6080| 国产成人福利小说| 亚洲国产精品999在线| 十八禁人妻一区二区| netflix在线观看网站| 男女做爰动态图高潮gif福利片| 亚洲av免费高清在线观看| 久久精品国产亚洲av涩爱 | 给我免费播放毛片高清在线观看| 欧美色视频一区免费| 久久久国产成人精品二区| 亚洲av不卡在线观看| 精华霜和精华液先用哪个| 国产一区二区激情短视频| 桃色一区二区三区在线观看| 欧美+日韩+精品| 91av网一区二区| 欧美zozozo另类| 亚洲第一电影网av| 国产视频一区二区在线看| 色吧在线观看| 禁无遮挡网站| 免费在线观看日本一区| 久久久久久九九精品二区国产| 精品久久久久久久久av| 亚洲av免费高清在线观看| 免费一级毛片在线播放高清视频| 日韩免费av在线播放| 欧美成人性av电影在线观看| 国产精品不卡视频一区二区 | 国产激情偷乱视频一区二区| 97超视频在线观看视频| 国产美女午夜福利| 校园春色视频在线观看| 十八禁人妻一区二区| 欧美成人免费av一区二区三区| 熟妇人妻久久中文字幕3abv| 久久久久亚洲av毛片大全| 久久久久久久午夜电影| 欧美乱妇无乱码| 99热精品在线国产| 此物有八面人人有两片| 午夜福利免费观看在线| 国产日本99.免费观看| 在线播放无遮挡| 五月伊人婷婷丁香| 色精品久久人妻99蜜桃| 夜夜躁狠狠躁天天躁| 91午夜精品亚洲一区二区三区 | 免费高清视频大片| 亚洲欧美激情综合另类| 亚洲国产精品合色在线| 宅男免费午夜| 国产精品嫩草影院av在线观看 | 成人国产综合亚洲| 免费观看人在逋| 欧美潮喷喷水| 亚洲国产欧洲综合997久久,| 一本精品99久久精品77| 婷婷六月久久综合丁香| 蜜桃亚洲精品一区二区三区| 国产午夜精品久久久久久一区二区三区 | 成人毛片a级毛片在线播放| 欧美成人一区二区免费高清观看| 天天一区二区日本电影三级| 中文字幕熟女人妻在线| 看片在线看免费视频| 亚洲av美国av| 国产色爽女视频免费观看| 久久精品久久久久久噜噜老黄 | 午夜福利成人在线免费观看| 日日夜夜操网爽| 夜夜夜夜夜久久久久| 精品久久久久久久人妻蜜臀av| 在线观看av片永久免费下载| 欧美性感艳星| 国产成人啪精品午夜网站| 国产精品一区二区性色av| 亚洲无线在线观看| 中出人妻视频一区二区| h日本视频在线播放| 长腿黑丝高跟| 国产亚洲欧美98| 亚洲av五月六月丁香网| 美女黄网站色视频| 欧美xxxx性猛交bbbb| 两性午夜刺激爽爽歪歪视频在线观看| 国产免费av片在线观看野外av| 久99久视频精品免费| 日本熟妇午夜| 日本精品一区二区三区蜜桃| 国产69精品久久久久777片| 又黄又爽又免费观看的视频| 国内毛片毛片毛片毛片毛片| 在线播放无遮挡| 国产亚洲欧美在线一区二区| 成人欧美大片| 久久精品综合一区二区三区| 国产成人a区在线观看| 久久99热这里只有精品18| 特大巨黑吊av在线直播| 成人三级黄色视频| 婷婷丁香在线五月| 97人妻精品一区二区三区麻豆| x7x7x7水蜜桃| 欧美又色又爽又黄视频| 亚洲人成电影免费在线| 国产爱豆传媒在线观看| bbb黄色大片| 两个人的视频大全免费| 国产人妻一区二区三区在| 亚洲电影在线观看av| 一边摸一边抽搐一进一小说| 在线观看66精品国产| 亚洲av熟女| 内射极品少妇av片p| 舔av片在线| 免费一级毛片在线播放高清视频| 亚洲av日韩精品久久久久久密| 免费人成在线观看视频色| 国产精品一区二区三区四区久久| 久久久精品大字幕| 哪里可以看免费的av片| 啦啦啦韩国在线观看视频| 亚洲18禁久久av| 又黄又爽又刺激的免费视频.| 免费av毛片视频| 亚洲精品乱码久久久v下载方式| 3wmmmm亚洲av在线观看| 内地一区二区视频在线| 黄色女人牲交| 观看免费一级毛片| 精品午夜福利在线看| 男女那种视频在线观看| 免费高清视频大片| a级一级毛片免费在线观看| 有码 亚洲区| 99热精品在线国产| 亚洲av日韩精品久久久久久密| 丁香六月欧美| 90打野战视频偷拍视频| 91在线观看av| 搡老妇女老女人老熟妇| 欧美精品啪啪一区二区三区| 69人妻影院| 99国产精品一区二区三区| 一级毛片久久久久久久久女| 人妻制服诱惑在线中文字幕| 男女视频在线观看网站免费| 黄色配什么色好看| 亚洲av熟女| 99久久九九国产精品国产免费| 亚洲aⅴ乱码一区二区在线播放| 日本熟妇午夜| 久久6这里有精品| 99久久九九国产精品国产免费| 日韩av在线大香蕉| 精品不卡国产一区二区三区| 成年人黄色毛片网站| av在线蜜桃| 亚洲国产色片| 91午夜精品亚洲一区二区三区 | 久久伊人香网站| avwww免费| 亚洲久久久久久中文字幕| avwww免费| 色尼玛亚洲综合影院| 热99在线观看视频| 真实男女啪啪啪动态图| 在线看三级毛片| 欧美日韩黄片免| 麻豆成人午夜福利视频| 亚洲av成人不卡在线观看播放网| 99久久精品国产亚洲精品| 人人妻人人澡欧美一区二区| 特大巨黑吊av在线直播| www.色视频.com| 精品久久久久久久久久久久久| 一进一出好大好爽视频| 亚洲三级黄色毛片| 露出奶头的视频| 欧美日本视频| 一夜夜www| 99久久九九国产精品国产免费| 天天躁日日操中文字幕| 亚洲中文字幕日韩| 成人性生交大片免费视频hd| 亚洲精品一区av在线观看| 性欧美人与动物交配| 五月伊人婷婷丁香| 级片在线观看| 有码 亚洲区| 国产亚洲精品综合一区在线观看| 搡老妇女老女人老熟妇| 在线观看66精品国产| 午夜福利视频1000在线观看| 精品人妻偷拍中文字幕| 老司机午夜十八禁免费视频| 深夜a级毛片| 亚洲avbb在线观看| 久久久久久久久久黄片| 波多野结衣巨乳人妻| 18禁在线播放成人免费| 亚洲国产精品合色在线| 午夜免费激情av| 每晚都被弄得嗷嗷叫到高潮| 国内精品久久久久久久电影| 少妇被粗大猛烈的视频| 久久人妻av系列| 国产国拍精品亚洲av在线观看| 免费一级毛片在线播放高清视频| 我要搜黄色片| 国产精品1区2区在线观看.| 一本久久中文字幕| 亚洲 国产 在线| 女人被狂操c到高潮| 色噜噜av男人的天堂激情| 又爽又黄a免费视频| 欧美精品国产亚洲| 99久久99久久久精品蜜桃| 天美传媒精品一区二区| 欧美区成人在线视频| 青草久久国产| 波野结衣二区三区在线| 制服丝袜大香蕉在线| 99久久精品一区二区三区| 久久热精品热| 欧美成狂野欧美在线观看| 亚洲经典国产精华液单 | av在线老鸭窝| av欧美777| 男人狂女人下面高潮的视频| 国产精品久久久久久亚洲av鲁大| 好男人在线观看高清免费视频| 国产精品久久久久久久电影| 天天躁日日操中文字幕| 又黄又爽又免费观看的视频| 深爱激情五月婷婷| 欧美高清成人免费视频www| 久久精品综合一区二区三区| 国产精品乱码一区二三区的特点| 高清毛片免费观看视频网站| 日韩精品青青久久久久久| 亚洲一区高清亚洲精品| 久久精品国产亚洲av香蕉五月| 亚洲在线自拍视频| 嫩草影院精品99| 少妇高潮的动态图| 一区二区三区免费毛片| 国产精品一区二区三区四区免费观看 | 在线观看av片永久免费下载| av视频在线观看入口| 亚洲成av人片在线播放无| 久久精品影院6| 成人美女网站在线观看视频| 色吧在线观看| 欧美国产日韩亚洲一区| 午夜a级毛片| 欧美日韩亚洲国产一区二区在线观看| av天堂在线播放| 在线免费观看不下载黄p国产 | 国产三级中文精品| 亚洲aⅴ乱码一区二区在线播放| 天堂av国产一区二区熟女人妻| 欧美精品国产亚洲| 国产v大片淫在线免费观看| 国产一级毛片七仙女欲春2| 桃红色精品国产亚洲av| 国内精品一区二区在线观看| 精品免费久久久久久久清纯| 国产精品99久久久久久久久| 日韩高清综合在线| 国产三级中文精品| 国产成人啪精品午夜网站| 久久人人爽人人爽人人片va | 国产一区二区激情短视频| 伊人久久精品亚洲午夜| 禁无遮挡网站| 免费看a级黄色片| www.熟女人妻精品国产| 欧美日韩中文字幕国产精品一区二区三区| 国产精品一区二区三区四区免费观看 | 久久精品国产清高在天天线| 色综合婷婷激情| 色播亚洲综合网| 身体一侧抽搐| 老司机午夜福利在线观看视频| 久久中文看片网| 久久久久久久久中文| 久久亚洲真实| 好男人电影高清在线观看| 亚洲av电影不卡..在线观看| 少妇的逼好多水| 国产探花极品一区二区| 午夜免费激情av| 国产男靠女视频免费网站| 尤物成人国产欧美一区二区三区| 九九热线精品视视频播放| 日韩精品中文字幕看吧| 99久久无色码亚洲精品果冻| 网址你懂的国产日韩在线| 色播亚洲综合网| 老熟妇乱子伦视频在线观看| 99国产精品一区二区蜜桃av| 成人高潮视频无遮挡免费网站| 一区二区三区四区激情视频 | 一区二区三区激情视频| 欧美在线一区亚洲| 白带黄色成豆腐渣| 亚洲aⅴ乱码一区二区在线播放| 少妇高潮的动态图| 男女视频在线观看网站免费| av天堂中文字幕网| 嫁个100分男人电影在线观看| 国产精品久久久久久久电影| 亚洲av免费高清在线观看| 两个人的视频大全免费| 国产成人aa在线观看| 国产精品影院久久| 成年免费大片在线观看| 久久国产乱子免费精品| 老司机福利观看| 久久这里只有精品中国| 午夜久久久久精精品| 亚洲自偷自拍三级| 成人一区二区视频在线观看| 国产一区二区激情短视频| 18+在线观看网站| 欧美日韩乱码在线| 亚洲va日本ⅴa欧美va伊人久久| а√天堂www在线а√下载| 国产精品国产高清国产av| 一区二区三区四区激情视频 | 亚洲专区中文字幕在线| 免费在线观看成人毛片| 国产真实伦视频高清在线观看 | av黄色大香蕉| 五月玫瑰六月丁香| 欧美潮喷喷水| 亚洲国产欧美人成| 国产熟女xx| 99国产精品一区二区蜜桃av| 在线a可以看的网站| 久久九九热精品免费| 国产乱人视频| 国产亚洲欧美98| 十八禁人妻一区二区| 免费在线观看影片大全网站| 亚洲最大成人手机在线| 全区人妻精品视频| 一进一出抽搐动态| 亚洲av.av天堂| 国产午夜精品论理片| 色播亚洲综合网| 久久精品91蜜桃| 51国产日韩欧美| 高清日韩中文字幕在线| 国产探花极品一区二区| 能在线免费观看的黄片| 色综合亚洲欧美另类图片| 一本综合久久免费| 99在线视频只有这里精品首页| 午夜a级毛片| 9191精品国产免费久久| 欧美乱妇无乱码| 最新中文字幕久久久久| 亚洲午夜理论影院| 青草久久国产| 男女下面进入的视频免费午夜| 色综合婷婷激情| 国产亚洲精品久久久com| 搡老岳熟女国产| 麻豆一二三区av精品| 高清毛片免费观看视频网站| 五月伊人婷婷丁香| 午夜免费男女啪啪视频观看 | 夜夜爽天天搞| 99国产极品粉嫩在线观看| 淫妇啪啪啪对白视频| 91久久精品国产一区二区成人| 欧美又色又爽又黄视频| 亚洲在线观看片| 亚洲黑人精品在线| 久久久久亚洲av毛片大全| 天堂网av新在线| 最近最新免费中文字幕在线| 亚洲成人精品中文字幕电影| 日日摸夜夜添夜夜添小说| 18禁裸乳无遮挡免费网站照片| 搡老岳熟女国产| 在线国产一区二区在线| 国产精品嫩草影院av在线观看 | 亚洲国产欧美人成|