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

    渤海海域應(yīng)急救援基地選址優(yōu)化方法

    2017-11-14 03:10:50李歡歡吳克風(fēng)
    關(guān)鍵詞:渤海適應(yīng)度海域

    李歡歡,劉 奕,劉 文,吳克風(fēng)

    (1.武漢理工大學(xué) 航運(yùn)學(xué)院,湖北 武漢 430063; 2.華中科技大學(xué) 自動(dòng)化學(xué)院,湖北 武漢 430074)

    ?

    渤海海域應(yīng)急救援基地選址優(yōu)化方法

    李歡歡1,劉奕1,劉文1,吳克風(fēng)2

    (1.武漢理工大學(xué) 航運(yùn)學(xué)院,湖北 武漢 430063; 2.華中科技大學(xué) 自動(dòng)化學(xué)院,湖北 武漢 430074)

    針對(duì)海上應(yīng)急救援基地規(guī)劃比較單一的問題,結(jié)合渤海海域交通流特征和海上交通事故的統(tǒng)計(jì)結(jié)果,提出了應(yīng)急救援基地雙覆蓋優(yōu)化模型。綜合考慮了渤海海域環(huán)境條件和事故發(fā)生率的影響,基于地理網(wǎng)格劃分方法確定了渤海海域應(yīng)急事故點(diǎn),根據(jù)港口條件選定了應(yīng)急救援基地候選點(diǎn)。應(yīng)用雙覆蓋模型對(duì)渤海海域應(yīng)急救援基地選址進(jìn)行了優(yōu)化分析,利用遺傳算法優(yōu)化求解,得到渤海海域應(yīng)急救援基地選址的最優(yōu)布局方案。研究結(jié)果表明:提出的應(yīng)急基地雙覆蓋模型比現(xiàn)有的選址模型有更好的覆蓋效果,且適應(yīng)度函數(shù)值高達(dá)0.92,充分證明了模型的有效性。

    遺傳算法;海上應(yīng)急救援基地;海上事故風(fēng)險(xiǎn);雙覆蓋模型;布局優(yōu)化

    0 引言

    船舶大型化、高密度化的發(fā)展使得海上事故發(fā)生的概率逐漸增大,導(dǎo)致海上交通事故的后果也越來越嚴(yán)重,因此,海上交通安全和突發(fā)事故的應(yīng)急管理成為研究的重點(diǎn)[1]。船舶交通流量大、交通繁忙、海上石油鉆井平臺(tái)多、自然資源豐富而敏感等特征構(gòu)成了渤海海域復(fù)雜的通航環(huán)境,也使得渤海海域的交通安全形勢(shì)更加嚴(yán)峻。為了進(jìn)一步提高應(yīng)急資源利用效率和應(yīng)急服務(wù)水平,有必要對(duì)海上應(yīng)急救援基地的布局規(guī)劃和覆蓋分配進(jìn)行優(yōu)化研究[2]。

    文獻(xiàn)[3]基于救護(hù)車排隊(duì)系統(tǒng)提出了新的救護(hù)車動(dòng)態(tài)分配模型。文獻(xiàn)[4]基于經(jīng)典覆蓋模型提出了總需求覆蓋最小的設(shè)施選址模型,且在理論上證明了解的唯一性。文獻(xiàn)[5]將緊急醫(yī)療服務(wù)的應(yīng)急設(shè)施選址構(gòu)建成一個(gè)兩階段的隨機(jī)規(guī)劃問題,提出了不確定環(huán)境下的選址模型。文獻(xiàn)[6]根據(jù)每個(gè)需求點(diǎn)的需求量是隨機(jī)的這一特點(diǎn)提出了不同的選址模型,并使得需求點(diǎn)和應(yīng)急中心直線距離最大時(shí)的期望最小化。文獻(xiàn)[7]提出了處理危險(xiǎn)品事故的雙層次應(yīng)急中心選址模型,加入層次管理參數(shù)λ并進(jìn)行了靈敏度分析,從而降低了應(yīng)急成本,提高了應(yīng)急效率。文獻(xiàn)[8]對(duì)美國(guó)芝加哥的大型城市道路交通進(jìn)行了研究,提出了一種醫(yī)療應(yīng)急資源雙層覆蓋模型,配置不同等級(jí)服務(wù)設(shè)施的救護(hù)車,增加覆蓋面積以降低傷亡率。文獻(xiàn)[9]研究了渤海海域溢油應(yīng)急基地的多目標(biāo)優(yōu)化選址問題,并基于遺傳算法求解多目標(biāo)規(guī)劃。文獻(xiàn)[10]構(gòu)建了多目標(biāo)最大覆蓋模型并利用遺傳算法優(yōu)化求解。

    目前,很多選址研究都是針對(duì)公路交通,海上應(yīng)急救援基地選址的研究還不深入,其選址模型大多數(shù)基于較為簡(jiǎn)單的規(guī)劃模型。本文基于文獻(xiàn)[8-9],結(jié)合公路設(shè)施選址和海上應(yīng)急救援經(jīng)典問題,綜合考慮救援時(shí)間、救援距離、雙覆蓋和基地設(shè)置數(shù)量等約束條件,提出了一種基于最短距離和最少應(yīng)急救援基地?cái)?shù)的渤海海域應(yīng)急救援基地雙覆蓋優(yōu)化模型,然后利用遺傳算法進(jìn)行優(yōu)化求解。

    1 現(xiàn)有模型分析

    應(yīng)急基地選址的優(yōu)化目標(biāo)一般有:應(yīng)急時(shí)間最短;在滿足一定時(shí)間緊迫性的條件下,到達(dá)各個(gè)應(yīng)急救援點(diǎn)的路程之和最短;在預(yù)算限制條件下,滿足一定的服務(wù)需求目標(biāo),使應(yīng)急基地設(shè)置數(shù)目最?。辉跐M足距離最短或時(shí)間最短限制條件下,提供應(yīng)急救援的滿意度最大等?,F(xiàn)有模型主要有:

    (Ⅰ)P-中值模型

    主要是考慮在成本最低的基礎(chǔ)上,如何使服務(wù)設(shè)施與服務(wù)點(diǎn)之間的總運(yùn)輸距離或平均運(yùn)輸距離最小,以使需求量和運(yùn)輸距離乘積之和最小,保證經(jīng)濟(jì)效益最大化的目標(biāo),但是對(duì)服務(wù)響應(yīng)的及時(shí)性要求不高。P-中值模型如下:

    (1)

    其中:N為需求點(diǎn)的集合;d為需求點(diǎn)的需求量;M為擬建的設(shè)施候選點(diǎn);c為從服務(wù)設(shè)施點(diǎn)到需求點(diǎn)的運(yùn)輸費(fèi)用;p為可建的設(shè)施總數(shù)。該模型對(duì)服務(wù)反應(yīng)時(shí)間要求高的應(yīng)急設(shè)施選址問題不適用。

    (Ⅱ)P-中心模型

    文獻(xiàn)[11]考慮服務(wù)設(shè)施對(duì)服務(wù)點(diǎn)的平均服務(wù)能力最優(yōu),提出了網(wǎng)絡(luò)上的P-中心問題,即選P個(gè)服務(wù)設(shè)施的位置,使服務(wù)條件最壞的服務(wù)點(diǎn)的服務(wù)狀況盡量好,如縮短最大反應(yīng)時(shí)間、減小最大距離或縮小最大損失等,但是并沒有考慮需求點(diǎn)的具體信息。P-中心模型如下:

    (2)

    (Ⅲ)位置集合覆蓋模型

    文獻(xiàn)[12]提出了設(shè)定服務(wù)設(shè)施的服務(wù)半徑限制標(biāo)準(zhǔn),如為城市救護(hù)車車隊(duì)規(guī)定5 min或10 min的接警響應(yīng)時(shí)間,以此來設(shè)定責(zé)任區(qū)最遠(yuǎn)救援范圍,而此類約束條件在P-中值模型和P-中心模型中并沒有考慮,故產(chǎn)生了服務(wù)設(shè)施選址的覆蓋問題(covering problem,CP)。文獻(xiàn)[13]最早提出了位置集合覆蓋問題(location set covering problem,LSCP),要求使所有應(yīng)急點(diǎn)滿足被覆蓋的條件下,確定應(yīng)急服務(wù)設(shè)施的最小設(shè)置數(shù)量或建設(shè)費(fèi)用。位置集合覆蓋模型如下:

    (3)

    2 本文模型構(gòu)建

    上述模型都存在相應(yīng)的弊端,本文將全部海域考慮在應(yīng)急力量覆蓋范圍之內(nèi),綜合考慮覆蓋問題,建立多個(gè)應(yīng)急救援基地,以達(dá)到對(duì)更大范圍海域事故點(diǎn)的有效救援,則海上應(yīng)急救援基地的選址問題,轉(zhuǎn)化為多個(gè)應(yīng)急事故點(diǎn)對(duì)應(yīng)的應(yīng)急救援基地選址問題。本文的優(yōu)化目標(biāo)是保證應(yīng)急事故點(diǎn)能被最近的應(yīng)急救援基地覆蓋,同時(shí)每個(gè)應(yīng)急救援基地最多覆蓋兩個(gè)應(yīng)急事故點(diǎn),從而達(dá)到海上應(yīng)急救援配置資源的有效利用。應(yīng)急事故點(diǎn)的集合能夠被應(yīng)急救援基地候選點(diǎn)集合中的最少數(shù)量的應(yīng)急救援基地所覆蓋,使其在盡量能夠全面覆蓋海上應(yīng)急事故點(diǎn)的前提下,海上應(yīng)急救援系統(tǒng)的建設(shè)和運(yùn)營(yíng)費(fèi)用盡量低。

    本文構(gòu)建的應(yīng)急救援基地雙覆蓋模型如下:

    (4)

    基于遺傳算法求解的過程中,綜合考慮距離代價(jià)和雙覆蓋代價(jià),從而得到最優(yōu)組合以及應(yīng)急救援基地選址的最優(yōu)布局。

    3 渤海海域應(yīng)急救援基地選址優(yōu)化

    3.1應(yīng)急救援基地候選點(diǎn)位置確定

    海上應(yīng)急救援基地選址布局需要先確定應(yīng)急救援基地候選點(diǎn)的位置和數(shù)量,由于港口附近水域是船舶匯入?yún)R出和聚集的繁忙水域,事故發(fā)生概率比較大。另外,港口水域具有優(yōu)良的航道環(huán)境、較完善的規(guī)劃布局和疏通能力,有利于應(yīng)急資源的儲(chǔ)存、管理和調(diào)度。良好的港航條件和導(dǎo)助航設(shè)施能夠保證快速派出應(yīng)急基地的救助船舶,故將沿海港口作為應(yīng)急救援基地的候選點(diǎn)較為合適。根據(jù)海上應(yīng)急救援基地的實(shí)際建設(shè)情況,應(yīng)急救援基地大多設(shè)置在沿海港口,為港口附近水域的水上交通提供重要保障。

    渤海海域的專業(yè)海上應(yīng)急救援力量主要是交通運(yùn)輸部北海救助局和北海第一救助飛行隊(duì)。其中,北海救助局包括大連、秦皇島、天津、煙臺(tái)、榮成等5個(gè)救助基地和南隍城救助站。同時(shí)在渤海海峽、天津、北海3號(hào)位、石島、大連、煙臺(tái)、長(zhǎng)島、青島等8個(gè)北部海區(qū)設(shè)置了24個(gè)動(dòng)態(tài)值班待命點(diǎn),配備了8~9艘救助船舶進(jìn)行24 h海上應(yīng)急救助和動(dòng)態(tài)待命值班。本文選取的應(yīng)急救援基地候選點(diǎn)分別為大連、營(yíng)口、錦州、葫蘆島、秦皇島、唐山、天津、黃驊、東營(yíng)、蓬萊、煙臺(tái)、威海、萊州和濰坊等14個(gè)港口。

    3.2應(yīng)急點(diǎn)位置確定

    本文對(duì)渤海水域進(jìn)行地理網(wǎng)格化處理,并通過海域事故風(fēng)險(xiǎn)評(píng)價(jià)模型對(duì)網(wǎng)格水域的船舶事故風(fēng)險(xiǎn)進(jìn)行評(píng)價(jià)。通航水域網(wǎng)格化是根據(jù)水域管理要求,按照一定尺度將研究水域劃分為若干尺度相同或不同的網(wǎng)格單元。網(wǎng)格單位尺寸大小取決于研究對(duì)象水域的地理位置和大小。海上事故風(fēng)險(xiǎn)評(píng)價(jià)總體表示發(fā)生事故的風(fēng)險(xiǎn),是對(duì)海域安全形勢(shì)的一個(gè)綜合性評(píng)價(jià)。事故風(fēng)險(xiǎn)表達(dá)式[14]為:

    (5)

    其中:R為風(fēng)險(xiǎn);P為事故概率;I為事故可能產(chǎn)生的影響程度;n為研究區(qū)域類的船舶事故數(shù)量;N為海區(qū)內(nèi)船舶事故的總數(shù);S為船舶類型系數(shù);T為船舶大小系數(shù);A為事故海域系數(shù);D為事故等級(jí)系數(shù)。

    海上事故發(fā)生后,發(fā)生事故的船舶種類、船舶大小、發(fā)生事故的海域位置、發(fā)生事故的等級(jí)均在一定程度上決定了該事故可能造成的影響。為了評(píng)價(jià)海域內(nèi)船舶安全狀況,可以用安全指數(shù)法將船舶事故、船舶大小、事故海域系數(shù)用綜合換算系數(shù)[14]替代,可以將研究船型大小、事故種類換算成對(duì)應(yīng)系數(shù)。

    基于網(wǎng)格水域事故風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果,按照事故風(fēng)險(xiǎn)等級(jí)進(jìn)行排序,篩選出事故風(fēng)險(xiǎn)等級(jí)較高的水域作為應(yīng)急需求水域,并將網(wǎng)格水域中心作為應(yīng)急救援基地選址布局模型中的應(yīng)急事故點(diǎn)。最終根據(jù)渤海海域事故風(fēng)險(xiǎn)分布,篩選出事故風(fēng)險(xiǎn)等級(jí)最高的20個(gè)風(fēng)險(xiǎn)點(diǎn)作為應(yīng)急救援的應(yīng)急事故點(diǎn)。

    海上應(yīng)急救援基地候選點(diǎn)與應(yīng)急事故點(diǎn)之間的距離假定為直線距離,首先將經(jīng)緯度坐標(biāo)轉(zhuǎn)化成平面坐標(biāo)距離,即:

    (6)

    其中:Dis為應(yīng)急救援基地候選點(diǎn)到應(yīng)急事故點(diǎn)之間的距離;R1為赤道平均半徑;LatA為應(yīng)急事故點(diǎn)的北緯度數(shù);LonA為應(yīng)急事故點(diǎn)的東經(jīng)度數(shù);LatB為應(yīng)急救援基地的北緯度數(shù);LonB為應(yīng)急救援基地的東經(jīng)度數(shù)。

    應(yīng)急救援基地候選點(diǎn)和應(yīng)急事故點(diǎn)的經(jīng)緯度坐標(biāo)見表1和表2。

    表1 應(yīng)急救援基地候選點(diǎn)坐標(biāo)

    表2 應(yīng)急事故點(diǎn)坐標(biāo)

    4 基于遺傳算法的優(yōu)化分析

    預(yù)計(jì)到2020年,離岸100 n mile的海上應(yīng)急救助快速響應(yīng)時(shí)間不超過150 min,故本文設(shè)定討論覆蓋半徑為150~170 km,保證在救援范圍之內(nèi),同時(shí)保證應(yīng)急救援基地雙覆蓋最優(yōu)。

    4.1仿真實(shí)驗(yàn)步驟

    為了驗(yàn)證本文提出的選址優(yōu)化方法的可行性,在渤海海域進(jìn)行實(shí)例驗(yàn)證,仿真實(shí)驗(yàn)環(huán)境為Intel i7 2.60 GHz CPU,8 GB內(nèi)存,Windows 10操作系統(tǒng),MATLAB R2012b軟件。

    基于遺傳算法將渤海海域應(yīng)急救援基地布局優(yōu)化問題抽象為應(yīng)急事故點(diǎn)的目標(biāo)選擇問題,在滿足20個(gè)應(yīng)急事故點(diǎn)全部被覆蓋的前提下,尋找應(yīng)急救援基地的最優(yōu)布局方案,主要步驟如下:

    步驟1:確定14個(gè)應(yīng)急救援基地和20個(gè)事故點(diǎn)。

    步驟2:初始化種群,設(shè)置遺傳算法的參數(shù),確定適應(yīng)度函數(shù)。

    步驟3:利用交叉、變異算子進(jìn)行迭代求解。

    步驟4:對(duì)比分析不同覆蓋半徑R(R分別為150 km、155 km、160 km、165 km和170 km)情況下,應(yīng)急救援基地覆蓋情況以及雙覆蓋情況。

    步驟5:根據(jù)適應(yīng)度函數(shù)的最優(yōu)值來確定最優(yōu)覆蓋半徑,在最優(yōu)覆蓋半徑下研究不同種群大小的影響,確定最優(yōu)種群大小。

    步驟6:解碼得到應(yīng)急基地在雙覆蓋條件下的最優(yōu)布局方案。

    4.2仿真實(shí)驗(yàn)遺傳算法參數(shù)設(shè)置

    實(shí)驗(yàn)設(shè)置變化的最大閾值為0.1,以防止閾值太大陷入局部最優(yōu)解。連續(xù)迭代次數(shù)閾值為10 000;最大迭代次數(shù)為10 000,以保證能夠充分搜索到所有可能的解。設(shè)置種群大小(染色體個(gè)數(shù))為50,交叉概率為0.8,變異概率0.8,基因長(zhǎng)度為14×20=280。

    適應(yīng)度函數(shù)為f(x)=w1·log2(x1+1)+w2·(x2+1)-4,w1+w2=1,其中:x1為所有應(yīng)急事故點(diǎn)到各個(gè)應(yīng)急救援基地的最小距離之和與當(dāng)前染色體中應(yīng)急事故點(diǎn)到各個(gè)應(yīng)急救援基地距離和的比值,x1∈[0,1];x2為每個(gè)應(yīng)急救援基地覆蓋多于2個(gè)應(yīng)急事故點(diǎn)的余量和與最壞情況下每個(gè)應(yīng)急救援基地覆蓋多于2個(gè)應(yīng)急事故點(diǎn)的余量和的比值,x2∈[0,1];權(quán)重w1和w2分別為目標(biāo)函數(shù)中的距離要求和約束條件中覆蓋要求的權(quán)重,根據(jù)不同情況來調(diào)節(jié)權(quán)重大小,從而得到全局最優(yōu)解。適應(yīng)度函數(shù)越接近于1,說明實(shí)驗(yàn)結(jié)果越好。實(shí)驗(yàn)設(shè)置w1=0.3,w2=0.7,更加側(cè)重于雙覆蓋的要求,覆蓋半徑為160 km,在該范圍內(nèi)尋找應(yīng)急救援基地能夠?qū)?yīng)急事故點(diǎn)達(dá)到最短距離和雙覆蓋的最優(yōu)解。遺傳算法流程圖如圖1所示。

    圖1 遺傳算法流程圖

    4.3結(jié)果對(duì)比分析

    對(duì)比分析覆蓋半徑R分別為150 km、155 km、160 km、165 km和170 km的情況下對(duì)應(yīng)的應(yīng)急救援基地覆蓋情況,不同覆蓋半徑下適應(yīng)度函數(shù)對(duì)比結(jié)果如圖2所示。由圖2可知:當(dāng)R為150 km和155 km時(shí),適應(yīng)度函數(shù)的最優(yōu)值為0.81;當(dāng)R=160 km時(shí),適應(yīng)度函數(shù)值達(dá)到0.92;當(dāng)R=165 km時(shí),適應(yīng)度函數(shù)值沒有明顯增加,故最優(yōu)的覆蓋半徑為R=160 km。分析當(dāng)R=160 km時(shí),種群大小對(duì)應(yīng)急救援基地優(yōu)化布局的影響,分別對(duì)比分析種群大小T為30個(gè)、40個(gè)、50個(gè)、60個(gè)和70個(gè)時(shí),其適應(yīng)度函數(shù)值的變化。不同種群大小下適應(yīng)度函數(shù)結(jié)果對(duì)比如圖3所示。由圖3可知:當(dāng)種群大小為30個(gè)時(shí),迭代次數(shù)很少時(shí),其適應(yīng)度函數(shù)值已達(dá)到0.90,故其為最優(yōu)種群個(gè)數(shù)。

    圖2不同覆蓋半徑下適應(yīng)度函數(shù)結(jié)果對(duì)比圖

    圖3不同種群大小下適應(yīng)度函數(shù)結(jié)果對(duì)比圖

    不同覆蓋半徑下應(yīng)急事故點(diǎn)對(duì)應(yīng)的應(yīng)急救援基地如表3所示,同一行中同一標(biāo)號(hào)最多出現(xiàn)兩次(即雙覆蓋),同一列中標(biāo)號(hào)重復(fù)率越高說明覆蓋效果越好。不同覆蓋半徑下應(yīng)急救援基地覆蓋的應(yīng)急事故點(diǎn)個(gè)數(shù)如表4所示,表4同一行中數(shù)字2的個(gè)數(shù)越多說明效果越好。經(jīng)過進(jìn)一步分析對(duì)比可知:當(dāng)R=160 km時(shí),應(yīng)急救援基地所達(dá)到的雙覆蓋為最優(yōu),其雙覆蓋率最大。

    表3 不同覆蓋半徑下應(yīng)急事故點(diǎn)對(duì)應(yīng)的應(yīng)急救援基地

    注:表中符號(hào)表示應(yīng)急救援基地標(biāo)號(hào)。

    表4 不同覆蓋半徑下應(yīng)急救援基地覆蓋的應(yīng)急事故點(diǎn)個(gè)數(shù)

    注:0為未覆蓋應(yīng)急事故點(diǎn);1為覆蓋1個(gè)應(yīng)急事故點(diǎn);2為覆蓋2個(gè)應(yīng)急事故點(diǎn)。

    基于最優(yōu)覆蓋半徑R=160 km,不同種群大小下應(yīng)急事故點(diǎn)對(duì)應(yīng)的應(yīng)急救援基地如表5所示,表5中的符號(hào)表示應(yīng)急救援基地標(biāo)號(hào),同一列中標(biāo)號(hào)重復(fù)率越高說明覆蓋效果越好。不同種群大小下應(yīng)急救援基地覆蓋的應(yīng)急事故點(diǎn)個(gè)數(shù)如表6所示,經(jīng)過對(duì)比分析可知,在R=160 km的條件下,不同的種群大小都不影響最終的應(yīng)急救援基地覆蓋結(jié)果,充分驗(yàn)證了本文提出模型的有效性和算法的收斂性,算法收斂得到最優(yōu)值,故選址種群大小為30個(gè),此時(shí)的應(yīng)急救援基地的雙覆蓋為最優(yōu),進(jìn)一步證明了本文方法的可行性。

    表5 不同種群大小下應(yīng)急事故點(diǎn)對(duì)應(yīng)的應(yīng)急救援基地

    注:表中符號(hào)表示應(yīng)急救援基地標(biāo)號(hào);應(yīng)急事故點(diǎn)R=160 km。

    表6 不同種群大小下應(yīng)急救援基地覆蓋的應(yīng)急事故點(diǎn)個(gè)數(shù)

    注:應(yīng)急救援基地R=160 km;0為未覆蓋應(yīng)急事故點(diǎn);1為覆蓋1個(gè)應(yīng)急事故點(diǎn);2為覆蓋2個(gè)應(yīng)急事故點(diǎn)。

    應(yīng)急救援的最優(yōu)雙覆蓋最終結(jié)果為:大連、葫蘆島、秦皇島、唐山、天津、黃驊、東營(yíng)、蓬萊和濰坊共9個(gè)應(yīng)急救援基地都分別覆蓋了2個(gè)應(yīng)急事故點(diǎn),覆蓋情況分別為大連—j18、j19,葫蘆島—j4、j17,秦皇島—j15、j16,唐山—j2、j14,天津—j5、j6,黃驊—j3、j8,東營(yíng)—j1、j13,蓬萊—j11、j20,濰坊—j7、j9;煙臺(tái)和萊州分別覆蓋了1個(gè)應(yīng)急事故點(diǎn),故應(yīng)在大連、葫蘆島、秦皇島、唐山、天津、黃驊、東營(yíng)、蓬萊、濰坊、煙臺(tái)和萊州設(shè)置應(yīng)急救援基地。最優(yōu)雙覆蓋結(jié)果進(jìn)一步證明了應(yīng)急救援基地雙覆蓋的有效性。

    5 結(jié)束語

    本文在保證應(yīng)急事故點(diǎn)被全部覆蓋的基礎(chǔ)上,綜合考慮了應(yīng)急事故點(diǎn)和應(yīng)急救援基地之間的最短距離和雙覆蓋要求,提高了應(yīng)急救援基地布局與海域事故風(fēng)險(xiǎn)分布的適應(yīng)性,構(gòu)建了基于應(yīng)急救援基地的雙覆蓋應(yīng)急事故點(diǎn)模型。使用該模型對(duì)渤海海域應(yīng)急救援基地選址進(jìn)行實(shí)例分析,基于渤海海域環(huán)境條件和事故發(fā)生率,利用地理網(wǎng)格化方法得到20個(gè)事故點(diǎn),選定14個(gè)應(yīng)急救援基地候選點(diǎn),然后基于遺傳算法優(yōu)化求解。本文提出的應(yīng)急基地雙覆蓋模型比現(xiàn)有的選址模型有更大的覆蓋率,且適應(yīng)度函數(shù)值高達(dá)0.92,充分驗(yàn)證了其有效性和可行性。

    [1]李洪成,吳曉平,付鈺,等.海上應(yīng)急保障資源調(diào)度的多目標(biāo)模糊規(guī)劃模型[J].安全與環(huán)境學(xué)報(bào),2015,15(4):172-176.

    [2]詹斌,馮樂,宋文娟.水上突發(fā)事件應(yīng)急資源儲(chǔ)備點(diǎn)選址模型研究[J].武漢理工大學(xué)學(xué)報(bào),2015,37(8):31-36.

    [3]BROTCORNE L,LAPORTE G,SEMET F.Ambulance location and relocation models[J].European journal of operational research,2003,147:451-463.

    [4]BERMAN O,HUANG R.The minimum weighted covering location problem with distance constraints[J].Computers & operations research,2008,35(3):356-372.

    [5]BERALD P,BRUNI M E.A probabilistic model applied to emergency service vehicle location[J].European journal of operational research,2009,196(1):323-331.

    [6]CANBOLAT M S,MASSOW M V.Locating emergency facilities with random demand for risk minimization[J].Expert systems with applications,2011,38(1):10099-10106.

    [7]袁文燕,彭云,楊豐梅.基于危險(xiǎn)化學(xué)品事故的雙層次及中心選址模型[J].系統(tǒng)工程理論與實(shí)踐,2015,35(3):729-735.

    [8]LIU Y,ROSHANDEH A M,LI Z,et al.Heuristic approach for optimizing emergency medical services in road safety within large urban networks[J].Journal of transportation engineering,2014,140(9):04014043.

    [9]李云斌,劉敬賢,魏蕾,等.渤海海域溢油應(yīng)急基地多目標(biāo)優(yōu)化選址方法[J].哈爾濱工程大學(xué)學(xué)報(bào),2016,37(4):1-6.

    [10]艾云飛,呂靖,張麗麗.水上應(yīng)急物資儲(chǔ)備庫選址-分配優(yōu)化模型[J].大連海事大學(xué)學(xué)報(bào),2015,41(2):62-66.

    [11]HAKIMI S L.Optimal locations of switching centers and the absolute centers and medians of a graph[J].Operations research,1970,12(3):450-459.

    [12]REVELLE C S,EISELT H A.Location analysis:a synthesis and survey [J].European journal of operational research,2005,165(1):1-19.

    [13]TOREGAS C,REVELLE C.Optimal location under time or distance constraints[J].Papers in the regional science,1972,28(1):133-144.

    [14]吳兆麟,朱軍.海上交通工程[M].大連:大連海事大學(xué)出版社,2004:291-292.

    國(guó)家自然科學(xué)基金項(xiàng)目(51479156,51179147);湖北省科技支撐基金項(xiàng)目(對(duì)外合作類2015BHE004);武漢理工大學(xué)自主創(chuàng)新基金項(xiàng)目(2016-HY-B1-09)

    李歡歡(1989- ),女,河南鄭州人,博士生,主要研究方向?yàn)榻煌ōh(huán)境與安全保障.

    2016-07-22

    1672-6871(2017)01-0098-07

    10.15926/j.cnki.issn1672-6871.2017.01.020

    X951

    A

    猜你喜歡
    渤海適應(yīng)度海域
    渤海大學(xué)作品精選
    聲屏世界(2022年17期)2022-12-18 01:42:44
    改進(jìn)的自適應(yīng)復(fù)制、交叉和突變遺傳算法
    遺落海域
    埕島海域海上獨(dú)立樁拆除方案探討
    “烈焰”,在渤海邊燃燒
    廣東省海域使用統(tǒng)計(jì)分析
    遼代“斡魯朵”內(nèi)的渤海人移民
    基于空調(diào)導(dǎo)風(fēng)板成型工藝的Kriging模型適應(yīng)度研究
    少數(shù)民族大學(xué)生文化適應(yīng)度調(diào)查
    海域評(píng)估技術(shù)指引
    欧美一级a爱片免费观看看| 嫩草影院精品99| 久久6这里有精品| 午夜激情欧美在线| av专区在线播放| 永久网站在线| 99久久精品热视频| 精品久久久久久久末码| 内地一区二区视频在线| 欧美最新免费一区二区三区| 久久精品国产自在天天线| 国产亚洲最大av| 免费观看的影片在线观看| 中文字幕免费在线视频6| 国产午夜精品一二区理论片| 国产女主播在线喷水免费视频网站 | 国产国拍精品亚洲av在线观看| 成人午夜高清在线视频| 精品久久久久久成人av| 亚洲,欧美,日韩| 日产精品乱码卡一卡2卡三| 亚洲国产欧美在线一区| 国内精品宾馆在线| 夜夜爽夜夜爽视频| av在线蜜桃| 国产成年人精品一区二区| 欧美成人一区二区免费高清观看| 免费黄网站久久成人精品| 熟妇人妻久久中文字幕3abv| eeuss影院久久| 一级二级三级毛片免费看| 亚洲欧美中文字幕日韩二区| 麻豆精品久久久久久蜜桃| 国产欧美日韩精品一区二区| 身体一侧抽搐| 九色成人免费人妻av| 大香蕉97超碰在线| 国产免费又黄又爽又色| 男的添女的下面高潮视频| 亚洲av免费在线观看| 日韩欧美精品免费久久| 免费少妇av软件| 亚洲电影在线观看av| 亚洲av.av天堂| 97热精品久久久久久| 国产久久久一区二区三区| 人妻制服诱惑在线中文字幕| 日韩精品有码人妻一区| 日本爱情动作片www.在线观看| 超碰97精品在线观看| 色5月婷婷丁香| 国产精品国产三级专区第一集| 男女边吃奶边做爰视频| 国产精品一二三区在线看| 欧美高清成人免费视频www| 婷婷色麻豆天堂久久| 激情五月婷婷亚洲| 免费看av在线观看网站| 免费无遮挡裸体视频| 国产一区二区在线观看日韩| 18禁在线无遮挡免费观看视频| 九九爱精品视频在线观看| 看免费成人av毛片| 亚洲国产av新网站| 亚洲欧美清纯卡通| 波多野结衣巨乳人妻| 男女边摸边吃奶| 一级毛片电影观看| 国产单亲对白刺激| 看免费成人av毛片| 精品国产三级普通话版| 毛片女人毛片| 免费黄频网站在线观看国产| 国产高潮美女av| 国产黄a三级三级三级人| 免费观看的影片在线观看| 国产日韩欧美在线精品| 国产精品综合久久久久久久免费| 久久久精品欧美日韩精品| 人妻系列 视频| 亚洲精品一二三| 人人妻人人澡欧美一区二区| 欧美激情在线99| 九九爱精品视频在线观看| 成人亚洲欧美一区二区av| 日日撸夜夜添| 亚洲欧美中文字幕日韩二区| kizo精华| 又大又黄又爽视频免费| 午夜福利视频1000在线观看| 国产成人午夜福利电影在线观看| 日本午夜av视频| 日韩av在线大香蕉| 久久久a久久爽久久v久久| 白带黄色成豆腐渣| 久久久午夜欧美精品| 亚洲成人精品中文字幕电影| 丰满人妻一区二区三区视频av| 久久久精品94久久精品| 69人妻影院| 国产人妻一区二区三区在| 久久精品久久久久久久性| 国产激情偷乱视频一区二区| 一个人看视频在线观看www免费| 免费观看在线日韩| 精品午夜福利在线看| 成人性生交大片免费视频hd| 内射极品少妇av片p| 国产一区二区在线观看日韩| 欧美成人一区二区免费高清观看| 午夜精品一区二区三区免费看| 大话2 男鬼变身卡| 少妇裸体淫交视频免费看高清| 久久精品久久久久久久性| 国产一区二区三区综合在线观看 | 国产高潮美女av| 国产一区二区亚洲精品在线观看| 国产精品久久久久久久久免| 波野结衣二区三区在线| 亚洲国产欧美人成| 99久国产av精品| 人人妻人人澡人人爽人人夜夜 | 午夜精品国产一区二区电影 | 久久热精品热| 嫩草影院新地址| 亚洲美女搞黄在线观看| 国产精品久久久久久久电影| 色综合站精品国产| 久久99热6这里只有精品| 国产探花极品一区二区| 日韩精品有码人妻一区| 一本一本综合久久| 免费av观看视频| 日本色播在线视频| 国产亚洲av嫩草精品影院| 黄色欧美视频在线观看| 2021天堂中文幕一二区在线观| 日本免费a在线| 狠狠精品人妻久久久久久综合| 看免费成人av毛片| 一边亲一边摸免费视频| 狂野欧美白嫩少妇大欣赏| 少妇的逼水好多| 国产精品福利在线免费观看| 免费观看的影片在线观看| 久久这里有精品视频免费| 国产三级在线视频| 国语对白做爰xxxⅹ性视频网站| 欧美一级a爱片免费观看看| 精品人妻偷拍中文字幕| 国产黄色小视频在线观看| 午夜福利成人在线免费观看| 99视频精品全部免费 在线| 日韩三级伦理在线观看| 久久久久性生活片| 精品不卡国产一区二区三区| 久久久精品94久久精品| av线在线观看网站| freevideosex欧美| 午夜免费激情av| 日韩欧美一区视频在线观看 | 国产亚洲av片在线观看秒播厂 | 亚洲精品国产av成人精品| 国产激情偷乱视频一区二区| eeuss影院久久| 不卡视频在线观看欧美| 日日撸夜夜添| 婷婷色av中文字幕| 国产精品一区二区性色av| 综合色av麻豆| 亚洲av不卡在线观看| 三级男女做爰猛烈吃奶摸视频| 亚洲伊人久久精品综合| av在线天堂中文字幕| 久久6这里有精品| 能在线免费观看的黄片| 内地一区二区视频在线| 国产成人一区二区在线| 欧美xxⅹ黑人| 国内揄拍国产精品人妻在线| 日本免费a在线| 色视频www国产| 欧美高清成人免费视频www| 国产永久视频网站| 国产亚洲午夜精品一区二区久久 | 日韩 亚洲 欧美在线| 精品久久久噜噜| 欧美xxⅹ黑人| 久久久久久久久久黄片| 国产在视频线精品| 亚洲av在线观看美女高潮| 欧美精品一区二区大全| 性插视频无遮挡在线免费观看| 黄色欧美视频在线观看| 搡女人真爽免费视频火全软件| 国产人妻一区二区三区在| 国产一区二区三区综合在线观看 | 国产成人91sexporn| 能在线免费观看的黄片| 在线天堂最新版资源| 男女国产视频网站| 精品久久久久久成人av| 国产三级在线视频| 简卡轻食公司| 少妇高潮的动态图| 偷拍熟女少妇极品色| 国产 一区 欧美 日韩| freevideosex欧美| 成年人午夜在线观看视频 | 亚洲av电影在线观看一区二区三区 | 熟妇人妻不卡中文字幕| 亚洲激情五月婷婷啪啪| 亚洲av免费高清在线观看| 伦精品一区二区三区| 亚洲精品乱久久久久久| 久久久久久久国产电影| 免费看a级黄色片| 最近手机中文字幕大全| 不卡视频在线观看欧美| 国产视频首页在线观看| 熟妇人妻不卡中文字幕| av福利片在线观看| 国产一区亚洲一区在线观看| 久久精品国产亚洲网站| 亚洲最大成人中文| 欧美精品国产亚洲| 直男gayav资源| 2022亚洲国产成人精品| 老司机影院成人| 日韩国内少妇激情av| 亚洲av福利一区| 中文字幕制服av| 国产爱豆传媒在线观看| 大片免费播放器 马上看| 波多野结衣巨乳人妻| 国产精品久久久久久精品电影| 麻豆乱淫一区二区| 欧美xxⅹ黑人| 一本久久精品| 亚洲性久久影院| 欧美精品国产亚洲| 老女人水多毛片| 国产午夜精品久久久久久一区二区三区| 久久午夜福利片| 日韩成人伦理影院| 日本猛色少妇xxxxx猛交久久| 丝瓜视频免费看黄片| 国产黄色免费在线视频| 在线 av 中文字幕| 国产精品1区2区在线观看.| 搡老乐熟女国产| 真实男女啪啪啪动态图| 久久这里只有精品中国| 女的被弄到高潮叫床怎么办| 日韩欧美精品v在线| 亚洲av免费高清在线观看| 欧美不卡视频在线免费观看| 一级毛片电影观看| 亚洲丝袜综合中文字幕| 亚洲精品色激情综合| 成人亚洲精品av一区二区| 边亲边吃奶的免费视频| 在现免费观看毛片| 国内精品一区二区在线观看| 日本免费a在线| 可以在线观看毛片的网站| 国产精品久久久久久精品电影| 干丝袜人妻中文字幕| 边亲边吃奶的免费视频| 在现免费观看毛片| 成人亚洲欧美一区二区av| 日本免费a在线| 国产激情偷乱视频一区二区| 天堂√8在线中文| 精品亚洲乱码少妇综合久久| 热99在线观看视频| 精品99又大又爽又粗少妇毛片| 亚洲av国产av综合av卡| 国产精品一区www在线观看| 美女脱内裤让男人舔精品视频| 最近中文字幕高清免费大全6| 国产在线男女| 有码 亚洲区| 丝袜喷水一区| 免费电影在线观看免费观看| 亚洲图色成人| 国内揄拍国产精品人妻在线| 成人亚洲精品av一区二区| 亚洲精品久久久久久婷婷小说| 日本欧美国产在线视频| 亚洲第一区二区三区不卡| 日日摸夜夜添夜夜爱| 欧美高清成人免费视频www| xxx大片免费视频| 国产在线一区二区三区精| 五月伊人婷婷丁香| 欧美成人精品欧美一级黄| 免费高清在线观看视频在线观看| 亚洲欧美一区二区三区国产| 精品人妻偷拍中文字幕| 久99久视频精品免费| 夜夜看夜夜爽夜夜摸| 少妇猛男粗大的猛烈进出视频 | 欧美日韩在线观看h| 最近的中文字幕免费完整| 午夜福利在线观看免费完整高清在| 天天躁日日操中文字幕| 一个人看的www免费观看视频| 国产综合精华液| 午夜免费男女啪啪视频观看| 久久久久久国产a免费观看| 99热这里只有是精品50| 日日撸夜夜添| 国产综合懂色| 国产成人免费观看mmmm| 精品久久久久久久久久久久久| 26uuu在线亚洲综合色| 免费看av在线观看网站| 婷婷色麻豆天堂久久| 能在线免费观看的黄片| 中文字幕av成人在线电影| 三级毛片av免费| 欧美不卡视频在线免费观看| 嫩草影院精品99| 男女那种视频在线观看| 日韩大片免费观看网站| 丝袜喷水一区| 天堂影院成人在线观看| 波野结衣二区三区在线| 亚洲伊人久久精品综合| videossex国产| 亚洲国产成人一精品久久久| 色网站视频免费| 舔av片在线| av线在线观看网站| 国产一级毛片七仙女欲春2| 一级毛片电影观看| 日本三级黄在线观看| 91久久精品电影网| 69人妻影院| 成人综合一区亚洲| av国产免费在线观看| 在现免费观看毛片| 免费看光身美女| 亚洲精品色激情综合| 成人鲁丝片一二三区免费| 夜夜看夜夜爽夜夜摸| 欧美性感艳星| 男女那种视频在线观看| 久久99精品国语久久久| 欧美日韩国产mv在线观看视频 | 国产精品熟女久久久久浪| 精品不卡国产一区二区三区| 高清欧美精品videossex| 熟女电影av网| 久久久色成人| 精品人妻偷拍中文字幕| 免费在线观看成人毛片| 国产av在哪里看| 人妻系列 视频| 日产精品乱码卡一卡2卡三| 99热6这里只有精品| 亚洲精品乱码久久久久久按摩| 午夜爱爱视频在线播放| 神马国产精品三级电影在线观看| 色视频www国产| 日本与韩国留学比较| 人妻系列 视频| 亚洲真实伦在线观看| 久久精品国产亚洲网站| 欧美日韩在线观看h| 超碰97精品在线观看| 深夜a级毛片| 蜜桃亚洲精品一区二区三区| 少妇丰满av| 亚洲av电影不卡..在线观看| 亚洲国产精品专区欧美| 国产综合懂色| 一边亲一边摸免费视频| 欧美高清成人免费视频www| 我要看日韩黄色一级片| 高清视频免费观看一区二区 | 欧美日韩一区二区视频在线观看视频在线 | 国内少妇人妻偷人精品xxx网站| 亚洲精品日本国产第一区| 亚洲欧美一区二区三区黑人 | 在线观看一区二区三区| 别揉我奶头 嗯啊视频| 久久久成人免费电影| av线在线观看网站| 国产成年人精品一区二区| 亚洲三级黄色毛片| 日韩强制内射视频| 不卡视频在线观看欧美| 好男人视频免费观看在线| 99视频精品全部免费 在线| 禁无遮挡网站| 男的添女的下面高潮视频| 久久精品夜色国产| 亚洲自拍偷在线| 晚上一个人看的免费电影| 男人舔女人下体高潮全视频| 久久精品国产鲁丝片午夜精品| av免费在线看不卡| 午夜精品一区二区三区免费看| 日韩av免费高清视频| 亚洲三级黄色毛片| 搡老乐熟女国产| 99re6热这里在线精品视频| 女人十人毛片免费观看3o分钟| 三级经典国产精品| 午夜福利视频1000在线观看| 日韩人妻高清精品专区| a级毛色黄片| av播播在线观看一区| 欧美成人精品欧美一级黄| 搡女人真爽免费视频火全软件| 极品少妇高潮喷水抽搐| 亚洲内射少妇av| 18+在线观看网站| 欧美性感艳星| 网址你懂的国产日韩在线| av国产久精品久网站免费入址| 九九在线视频观看精品| 久99久视频精品免费| 蜜臀久久99精品久久宅男| 成人鲁丝片一二三区免费| eeuss影院久久| 久久精品国产鲁丝片午夜精品| 亚州av有码| 秋霞伦理黄片| 日本av手机在线免费观看| 亚洲色图av天堂| 精品久久久精品久久久| 纵有疾风起免费观看全集完整版 | 两个人的视频大全免费| 免费在线观看成人毛片| 男女下面进入的视频免费午夜| 久久久久久久久久久丰满| 亚洲自偷自拍三级| 97超碰精品成人国产| 一夜夜www| 波野结衣二区三区在线| 精品一区二区三区视频在线| 人人妻人人澡欧美一区二区| 九九爱精品视频在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久这里有精品视频免费| 成人漫画全彩无遮挡| 黄色一级大片看看| 禁无遮挡网站| 亚洲av中文av极速乱| 国产视频内射| 菩萨蛮人人尽说江南好唐韦庄| 精品久久久噜噜| 日韩强制内射视频| 成年版毛片免费区| 精品欧美国产一区二区三| 天天躁夜夜躁狠狠久久av| 婷婷色av中文字幕| 我的女老师完整版在线观看| 亚洲人成网站在线观看播放| 免费观看无遮挡的男女| 日本黄大片高清| 国产av在哪里看| 欧美3d第一页| 七月丁香在线播放| 午夜免费男女啪啪视频观看| 91久久精品国产一区二区成人| 国产老妇女一区| 亚洲无线观看免费| 91精品伊人久久大香线蕉| 亚洲国产高清在线一区二区三| 精品人妻熟女av久视频| 欧美日韩精品成人综合77777| 老女人水多毛片| 九色成人免费人妻av| 久久亚洲国产成人精品v| 久久韩国三级中文字幕| 国产高潮美女av| 婷婷色av中文字幕| 97超视频在线观看视频| 国产精品人妻久久久久久| 91久久精品国产一区二区成人| 免费观看av网站的网址| 国产精品一二三区在线看| 亚洲国产av新网站| 80岁老熟妇乱子伦牲交| 日韩av在线大香蕉| 国产高潮美女av| 久久人人爽人人片av| 日韩精品有码人妻一区| 亚洲丝袜综合中文字幕| 女人十人毛片免费观看3o分钟| 日本-黄色视频高清免费观看| 亚洲精品第二区| 黄片wwwwww| 黄色欧美视频在线观看| 国产黄色小视频在线观看| 99视频精品全部免费 在线| 最近的中文字幕免费完整| 少妇丰满av| 欧美97在线视频| 日本av手机在线免费观看| 亚洲国产高清在线一区二区三| 久久人人爽人人片av| 婷婷六月久久综合丁香| 波多野结衣巨乳人妻| 国产免费视频播放在线视频 | 蜜桃亚洲精品一区二区三区| 亚洲自偷自拍三级| 一个人看视频在线观看www免费| 亚洲欧美清纯卡通| av播播在线观看一区| videos熟女内射| 淫秽高清视频在线观看| 国产高清不卡午夜福利| 少妇熟女欧美另类| 久久久久久伊人网av| 国产精品人妻久久久久久| 美女内射精品一级片tv| 中文字幕av成人在线电影| 18禁动态无遮挡网站| 亚洲av成人av| 久久精品国产亚洲av涩爱| 国产在视频线精品| 亚洲成人精品中文字幕电影| 一级爰片在线观看| 国产精品精品国产色婷婷| 欧美区成人在线视频| 免费在线观看成人毛片| 我的女老师完整版在线观看| 人妻一区二区av| 日韩av在线大香蕉| 国产老妇女一区| 亚洲国产色片| 22中文网久久字幕| 国产片特级美女逼逼视频| 亚洲激情五月婷婷啪啪| 久久精品久久久久久久性| 欧美bdsm另类| 极品少妇高潮喷水抽搐| 久久久久久久久久黄片| 国产亚洲av片在线观看秒播厂 | 国产精品嫩草影院av在线观看| 国产午夜福利久久久久久| 99久国产av精品国产电影| 国内少妇人妻偷人精品xxx网站| 最近的中文字幕免费完整| 日韩欧美精品v在线| 欧美高清性xxxxhd video| 亚洲熟女精品中文字幕| 亚洲精品中文字幕在线视频 | 国产免费福利视频在线观看| 国产精品嫩草影院av在线观看| 亚洲精品自拍成人| av黄色大香蕉| 国产综合精华液| 国内精品美女久久久久久| 免费看日本二区| 日韩av免费高清视频| 成人亚洲精品一区在线观看 | 观看美女的网站| 国产免费福利视频在线观看| 亚洲国产成人一精品久久久| 欧美另类一区| 99视频精品全部免费 在线| 婷婷色av中文字幕| 高清av免费在线| 一个人看的www免费观看视频| 亚洲精品成人av观看孕妇| 最近最新中文字幕免费大全7| 欧美+日韩+精品| 精品一区二区三区视频在线| 日本wwww免费看| av在线亚洲专区| 五月天丁香电影| 少妇被粗大猛烈的视频| 看黄色毛片网站| 婷婷色综合www| 日韩中字成人| 国产又色又爽无遮挡免| 免费观看的影片在线观看| 欧美日韩精品成人综合77777| 天堂av国产一区二区熟女人妻| 大片免费播放器 马上看| 只有这里有精品99| 中文字幕久久专区| 精品久久久久久久末码| 国产午夜精品久久久久久一区二区三区| 久久精品国产亚洲网站| 国产精品国产三级国产av玫瑰| 日本熟妇午夜| 少妇人妻精品综合一区二区| 免费电影在线观看免费观看| 亚洲欧美日韩东京热| 在线免费观看不下载黄p国产| 爱豆传媒免费全集在线观看| 国产精品一区二区三区四区免费观看| 秋霞在线观看毛片| 亚洲av福利一区| 国产精品一区二区在线观看99 | 日韩av不卡免费在线播放| 国产大屁股一区二区在线视频| 嫩草影院新地址| 日韩精品青青久久久久久| 久久久久久久亚洲中文字幕| 午夜福利视频精品| av播播在线观看一区| av专区在线播放| 婷婷色综合www| 久久精品夜夜夜夜夜久久蜜豆| 欧美日本视频| 五月天丁香电影| 久久久久久伊人网av|