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

    欽州灣海上溢油擴(kuò)散數(shù)值模擬*

    2016-06-27 01:32:26牙韓爭(zhēng)李誼純董德信
    廣西科學(xué) 2016年2期
    關(guān)鍵詞:落潮漲潮欽州

    牙韓爭(zhēng),李誼純,董德信

    (廣西科學(xué)院 廣西近海海洋環(huán)境科學(xué)重點(diǎn)實(shí)驗(yàn)室,廣西南寧 530007)

    欽州灣海上溢油擴(kuò)散數(shù)值模擬*

    牙韓爭(zhēng),李誼純,董德信

    (廣西科學(xué)院 廣西近海海洋環(huán)境科學(xué)重點(diǎn)實(shí)驗(yàn)室,廣西南寧530007)

    摘要:【目的】研究欽州灣海上溢油擴(kuò)散特征及其影響因素。【方法】在對(duì)欽州灣水動(dòng)力進(jìn)行數(shù)值模擬的基礎(chǔ)上,利用MIKE 21/3 SA溢油分析模塊,對(duì)其海上溢油擴(kuò)散進(jìn)行數(shù)值模擬?!窘Y(jié)果】常風(fēng)況條件下,漲潮時(shí)刻發(fā)生溢油,溢油先向茅尾海方向漂移,待落潮后退出茅尾海;計(jì)算時(shí)段內(nèi),溢油向茅尾海方向及欽州灣外灣方向漂移的最遠(yuǎn)距離分別約為17.51 km和10.73 km,掃海面積約為71.83 km2。落潮時(shí)刻發(fā)生溢油,溢油則先向欽州灣外灣漂移,待漲潮后轉(zhuǎn)向欽州灣灣頸方向;計(jì)算時(shí)段內(nèi),油膜向欽州灣外灣方向漂移最遠(yuǎn)距離約17.96 km,油污掃海面積約為50.72 km2。【結(jié)論】欽州灣溢油漂移擴(kuò)散特征與溢油發(fā)生時(shí)刻及風(fēng)作用密切相關(guān),溢油時(shí)刻及風(fēng)作用對(duì)欽州灣溢油漂移擴(kuò)散的影響不容忽視。

    關(guān)鍵詞:欽州灣海上溢油數(shù)值模擬

    0引言

    【研究意義】石油對(duì)工業(yè)生產(chǎn)和交通運(yùn)輸業(yè)具有不可代替的作用。隨著海上運(yùn)輸業(yè)的發(fā)展,沿海港口發(fā)生溢油事故的可能性不斷增加[1]。石油中含有烴類(lèi)有毒有害物質(zhì),流入海洋,會(huì)對(duì)海洋環(huán)境造成嚴(yán)重影響,破壞海洋生態(tài)環(huán)境,危害海洋生態(tài)系統(tǒng),因此,開(kāi)展海上溢油數(shù)值模擬研究,科學(xué)預(yù)測(cè)海上溢油擴(kuò)散特征,對(duì)海上溢油防范、擴(kuò)散控制、應(yīng)急處理、減少溢油對(duì)海洋生態(tài)環(huán)境的影響具有重大意義。而欽州灣是我國(guó)西南出??诘闹匾ǖ??!爸袊?guó)-東盟”貿(mào)易合作的加強(qiáng)以及廣西經(jīng)濟(jì)的快速發(fā)展,欽州灣一帶船舶運(yùn)輸業(yè)更為繁忙,使得欽州灣海域發(fā)生溢油事故的可能性大為提高,而欽州灣是一個(gè)半封閉式海灣,水交換能力差,周邊又擁有著豐富的海洋資源,有多種典型生態(tài)系統(tǒng),溢油事故一旦發(fā)生,極有可能對(duì)海洋生態(tài)、環(huán)境造成極大的影響?!厩叭搜芯窟M(jìn)展】自20世紀(jì)60年代,有國(guó)外學(xué)者開(kāi)始對(duì)海上溢油進(jìn)行研究,F(xiàn)ay[2]對(duì)油膜擴(kuò)散做出解釋?zhuān)岢鼋?jīng)典的三階理論,認(rèn)為油膜的擴(kuò)展主要受慣性力、重力、粘性力和表面張力控制,并將溢油擴(kuò)展分為慣性力擴(kuò)展階段、粘性力擴(kuò)展階段和表面張力擴(kuò)展3個(gè)階段。Blokker等[3]建立單組份的溢油風(fēng)化模型,但其未考慮油性質(zhì)隨時(shí)間的變化,故計(jì)算結(jié)果不能很好反應(yīng)溢油變化的真實(shí)情況。Mackay等[4]將風(fēng)速和溢油量加入模型中,建立了考慮風(fēng)速和溢油量的蒸發(fā)經(jīng)驗(yàn)公式,該理論已成為溢油風(fēng)化研究的基礎(chǔ)。Harrison等[5]提出計(jì)算多組分油的蒸發(fā)公式,使得溢油風(fēng)化計(jì)算更接近于實(shí)際。Lehr等[6]在Fay的理論基礎(chǔ)上對(duì)溢油擴(kuò)展進(jìn)行改進(jìn),在油膜的擴(kuò)展中考慮風(fēng)作用。Johansen等[7]基于lagrange追蹤法,改進(jìn)了“油粒子”模型,該模型可直觀地反應(yīng)油膜在風(fēng)、波浪等作用下破碎、分離的實(shí)際情況。Zhang等[8]考慮油膜的輸運(yùn)和風(fēng)化過(guò)程,建立一種海上溢油行為的預(yù)報(bào)模式。張存智等[9]建立三維溢油動(dòng)態(tài)預(yù)報(bào)模型,并將其成功應(yīng)用于渤海灣溢油數(shù)值模擬。楊紅等[10]基于拉格朗日粒子追蹤法及自由走動(dòng)理論建立長(zhǎng)江口溢油模型對(duì)南港水道溢油進(jìn)行模擬,模擬結(jié)果較好地反應(yīng)了溢油在長(zhǎng)江口南港水道的擴(kuò)散情況。另外,黃娟等[11]、李懷明等[12]、陳義中等[13]學(xué)者也進(jìn)行過(guò)類(lèi)似的研究。對(duì)于多種數(shù)值模型,目前應(yīng)用較為廣泛的有美國(guó)研究所(Applied science Associates,ASA)研發(fā)的OILMAP模型、MARKREED建立的OSCAR模型[14],以及丹麥水利研究所(DHI)研究開(kāi)發(fā)的MIKE 21 OS/SA溢油系列模型。【本研究切入點(diǎn)】欽州灣海域的溢油擴(kuò)散特征未曾有報(bào)道?!緮M解決的關(guān)鍵問(wèn)題】采用DHI研發(fā)的MIKE 21/3 SA溢油分析模塊,建立欽州灣海上溢油模型,對(duì)欽州灣海上溢油漂移擴(kuò)散進(jìn)行計(jì)算分析,探索欽州灣海上溢油漂移擴(kuò)散特征及其影響因素,為欽州灣海上溢油事故防范、溢油事件應(yīng)急處理、海洋生態(tài)環(huán)境保護(hù)等提供科學(xué)依據(jù)。

    1水動(dòng)力條件數(shù)值模擬

    水動(dòng)力條件數(shù)值模擬是溢油漂移擴(kuò)散預(yù)測(cè)的基礎(chǔ),其計(jì)算結(jié)果可為溢油擴(kuò)散預(yù)測(cè)模型提供潮流項(xiàng)數(shù)據(jù)。本文采用DHI開(kāi)發(fā)的MIKE 21 HD模塊對(duì)海域水動(dòng)力條件進(jìn)行計(jì)算。

    1.1模型控制方程

    MIKE 21 HD水動(dòng)力模型控制方程包括一個(gè)連續(xù)方程及兩個(gè)動(dòng)量方程。

    連續(xù)方程:

    (1)

    X方向動(dòng)量方程:

    (2)

    Y方向動(dòng)量方程:

    (3)

    式中:ζ為水位;h為水深;H=h+ζ;p、g分別為x、y方向的單寬流量;C為謝才系數(shù),與曼寧系數(shù)的關(guān)系為C=n-1/h1/6;g為重力加速度,f為風(fēng)摩擦系數(shù);V、Vx、Vy分別為風(fēng)速及其在x、y方向的分量;Ω為科氏力系數(shù);ρw為水的密度;x、y為距離;t為時(shí)間; u、v分別為流速在x、y方向的分量。τxx、τxy、τyy分別為各方向剪切應(yīng)力。

    1.2定解條件

    1.3計(jì)算區(qū)域的確定及網(wǎng)格劃分

    欽州灣海域岸線曲折,島嶼眾多,地形變化大,潮波運(yùn)動(dòng)受地形控制明顯,要準(zhǔn)確模擬海域內(nèi)潮波傳播過(guò)程,需建立足夠大的計(jì)算區(qū)域,因此,本次模型計(jì)算區(qū)域西邊界取至企沙鎮(zhèn),東邊界取至中三墩,南邊界取至欽州灣外灣,北邊界取至整個(gè)茅尾海,模型東西方向長(zhǎng)約38 km,南北方向長(zhǎng)約47 km,總面積約820 km2,計(jì)算區(qū)域地形圖如圖1所示。采用適應(yīng)性較好的不規(guī)則結(jié)構(gòu)網(wǎng)格對(duì)計(jì)算區(qū)域進(jìn)行劃分,并就復(fù)雜的島嶼和岸線進(jìn)行局部加密,以保證計(jì)算精度。計(jì)算網(wǎng)格共包涵10 280個(gè)節(jié)點(diǎn)、19 389個(gè)單元,最小邊長(zhǎng)為50 m,最大水深約21 m。

    圖1欽州灣水下地形及潮流要素驗(yàn)證點(diǎn)位

    Fig.1Schematic diagram of underwater topography of Qinzhou Bay and positions of model validation point

    1.4計(jì)算結(jié)果驗(yàn)證

    為驗(yàn)證水動(dòng)力模型的可靠性,采用2012年5月(小潮)及2012年6月(大潮)實(shí)測(cè)水文資料對(duì)計(jì)算結(jié)果進(jìn)行驗(yàn)證,包括2個(gè)潮位驗(yàn)證點(diǎn)及3個(gè)潮流驗(yàn)證點(diǎn),各潮流要素驗(yàn)證點(diǎn)點(diǎn)位如圖1所示,其中犀牛腳、龍門(mén)為潮位驗(yàn)證點(diǎn),V1~V3為潮流驗(yàn)證點(diǎn)。各點(diǎn)位潮流要素驗(yàn)證結(jié)果如圖2~5所示,由于篇幅有限,這里只列舉部分點(diǎn)位的結(jié)果驗(yàn)證。

    圖2犀牛腳潮位驗(yàn)證圖(大潮)

    Fig.2Verification of tidal level in Xiniujiao(spring tide)

    圖3龍門(mén)潮位驗(yàn)證圖(大潮)

    Fig.3Verification of tidal level in Longmen(spring tide)

    圖4V1站流速驗(yàn)證(大潮)

    Fig.4Verification of current velocity in V1(spring tide)

    圖5V1站流向驗(yàn)證(大潮)

    Fig.5Verification of current direction in V1(spring tide)

    由圖2~5可知,各潮位驗(yàn)證點(diǎn)潮位計(jì)算值與實(shí)測(cè)值有一定的偏差,但相差不大,最大潮差不超過(guò)0.3 m,龍門(mén)及犀牛腳兩站潮位過(guò)程與實(shí)際監(jiān)測(cè)結(jié)果幾乎一致,沒(méi)有太大的差別。V1站潮流計(jì)算值高潮位變化過(guò)程稍有偏差,個(gè)別時(shí)段流速計(jì)算值比實(shí)測(cè)值略小,但流速和流向的變化過(guò)程與實(shí)測(cè)值基本吻合。總體上,模型計(jì)算值與實(shí)測(cè)值吻合良好,位相、量值及過(guò)程線與實(shí)際情況基本一致,可以用于溢油擴(kuò)散數(shù)值模擬的計(jì)算研究。

    2溢油輸移擴(kuò)散模型的建立

    溢油進(jìn)入海洋水體后,在自身重力及海洋水體物理化學(xué)的作用下,同時(shí)發(fā)生擴(kuò)展、漂移、擴(kuò)散、蒸發(fā)、乳化、溶解等風(fēng)化過(guò)程。本文采用DHI研發(fā)的MIKE 21/3 SA溢油分析模塊建立欽州灣海上溢油模型,對(duì)欽州灣海洋溢油進(jìn)行預(yù)測(cè)。MIKE 21/3 SA模塊基于拉格朗日理論,采用“油粒子”方法,可模擬溢油在水體中漂移、擴(kuò)散、蒸發(fā)、乳化、溶解等多種過(guò)程。

    2.1擴(kuò)展

    油粒子擴(kuò)展是重力、慣性力、粘性力和表面張力等共同作用的結(jié)果。根據(jù)修正的Fay重力-粘力公式,油膜擴(kuò)散計(jì)算式為

    (4)

    2.2漂移

    油粒子漂移是潮流運(yùn)動(dòng)及風(fēng)拽力共同作用的結(jié)果,其漂移速度與潮流流速及風(fēng)速有關(guān)。油粒子漂移速度可用以下權(quán)重公式計(jì)算:

    Uoil=Vs+CwindUwind。

    (5)

    式中,Uwind為水面以上10 m處的風(fēng)速,Cwind為風(fēng)漂移系數(shù),一般取值在0.03~0.04,Vs為表面流速度。

    2.3擴(kuò)散

    油粒子擴(kuò)散與水流紊動(dòng)有關(guān),根據(jù)隨機(jī)擴(kuò)散理論,假設(shè)各方向擴(kuò)散具有同一性,一個(gè)時(shí)間步長(zhǎng)內(nèi)“油粒子”在a方向的擴(kuò)散距離可表示為

    (6)

    式中,Da為時(shí)間步長(zhǎng)內(nèi)溢油在a方向上的擴(kuò)散距離,R為-1~1的隨機(jī)數(shù),Ca為a方向上的擴(kuò)散系數(shù),Δt為時(shí)間步長(zhǎng)。

    2.4蒸發(fā)

    油膜蒸發(fā)速率與氣溫、水溫、溢油面積、風(fēng)速、太陽(yáng)輻射和油膜厚度等因素有關(guān),可用以下公式表示:

    (7)

    2.5乳化

    溢油在經(jīng)歷一定時(shí)間的漂移擴(kuò)散后,油膜厚度變小,水流紊動(dòng)作用對(duì)油膜的影響逐漸增大。在水流紊動(dòng)下,部分油膜被紊動(dòng)力能量撕裂成油滴,形成水包油的粒子,從而形成一種油水相交乳化物。油膜乳化可用油膜含水率來(lái)表示,其計(jì)算式為

    (8)

    其中,KB=1/YWF,YW為乳化物的含水量,KA為經(jīng)驗(yàn)參數(shù),取值約為4.5×10-6,YWF為油膜最終含水率。

    2.6溶解

    溢油的溶解主要發(fā)生在低碳輕質(zhì)油組分,溶解量通常很小,其溶解率用下式表示:

    (9)

    3欽州灣海上溢油擴(kuò)散預(yù)測(cè)

    3.1溢油工況設(shè)計(jì)及參數(shù)選取

    據(jù)《船舶油污染事故等級(jí)標(biāo)準(zhǔn)》,溢油事故溢油量大于10 t時(shí),稱(chēng)之為重大溢油事故。參照近年來(lái)各地溢油事故溢油量,本次研究溢油源強(qiáng)取為10 t,并假設(shè)1 h內(nèi)全部溢出。溢油點(diǎn)選擇欽州港航道附近,這與溢油事故通常發(fā)生于港口航道等航運(yùn)較繁忙區(qū)域的特點(diǎn)相適應(yīng)。大潮期間,水動(dòng)力能力相對(duì)較強(qiáng),溢油擴(kuò)散影響的范圍較廣,因此,本次溢油時(shí)刻選擇大潮時(shí)期,包括漲潮時(shí)和落潮時(shí)兩個(gè)時(shí)刻。根據(jù)欽州龍門(mén)氣象站資料,欽州灣年平均風(fēng)速不超過(guò)5.0 m/ s,最大風(fēng)速達(dá)28.0 m/s,常風(fēng)向?yàn)镹向風(fēng),頻率為26%,其次為南風(fēng),頻率為7%,綜合考慮,本次計(jì)算風(fēng)況選擇N、W、S、E 4個(gè)風(fēng)向,風(fēng)速取為5 m/s,另選一個(gè)靜風(fēng)情況。溢油模型各參數(shù)參照文獻(xiàn)[15],油源密度取0.94×103kg/m3,縱向擴(kuò)散系數(shù)取0.25 m/s2,橫向擴(kuò)散系數(shù)取0.1 m/s2。模擬時(shí)段設(shè)為24 h,并假定油膜最終擴(kuò)展厚度為0.01 mm[16]。各溢油計(jì)算工況及參數(shù)如表1所示,根據(jù)上述模型及參數(shù),對(duì)欽州灣海上溢油漂移擴(kuò)散進(jìn)行模擬。

    表1事故模擬情景條件

    Table 1Simulated scenes of accidents

    事故情形Accidentcase環(huán)境條件Environmentalcondition溢油量10tOilquantity10t風(fēng)向:N、W、S、E、靜風(fēng)Winddirection:N、W、S、E、staticwind油源密度Oildensity:0.94×103kg/m3風(fēng)速Windspeed:5m/s溢油位置:欽州港附近Oilspilledlocation:NearQinzhouPort氣溫Airtemperature:20℃溢油時(shí)刻:大潮(漲潮、落潮)Oilspilledtime:Springtide(Floodtide、Ebbtide)水溫Watertemperature:25℃溢油時(shí)間Oilspilledperiod:1h縱向擴(kuò)散系數(shù)Longitudinaldispersioncoefficient:0.25m/s2最終擴(kuò)展厚度Extendedthick-ness:0.01mm橫向擴(kuò)散系數(shù)Transversedis-persioncoefficient:0.10m/s2

    3.2結(jié)果與分析

    3.2.1溢油漂移擴(kuò)散特征

    為分析欽州灣海上溢油擴(kuò)散特征,在常風(fēng)向情況下,對(duì)欽州灣海上溢油進(jìn)行預(yù)測(cè)。圖6a為漲潮時(shí)刻發(fā)生溢油24 h溢油掃海范圍圖,可知,當(dāng)漲潮時(shí)刻發(fā)生溢油時(shí),溢油先向茅尾海方向移動(dòng),由欽州灣灣頸進(jìn)入茅尾海中部,待落潮時(shí)退出茅尾海直至欽州灣外灣,受溢油影響的區(qū)域主要為欽州灣灣頸至茅尾海中部以及欽州灣外灣中部一帶海域,以及灣頸中部部分島嶼岸線。在潮流及油膜自身重力的作用下,油膜發(fā)生位置上變化的同時(shí)不斷擴(kuò)散,造成油膜面積的擴(kuò)大。表2為常風(fēng)向條件下各計(jì)算時(shí)刻溢油擴(kuò)散計(jì)算結(jié)果統(tǒng)計(jì)表,可知,漲潮時(shí)刻發(fā)生溢油時(shí),24 h溢油擴(kuò)散面積約為9.80 km2,掃海面積則可達(dá)71.83 km2。整個(gè)計(jì)算時(shí)段,溢油向茅尾海方向漂移的距離約為17.51 km,而向欽州灣外灣方向漂移的距離約為10.73 km。

    圖6常風(fēng)向條件下24 h溢油掃海范圍

    Fig.6Sea swept area of spilled oil under N dominant wind conditions within 24 h

    圖6b為落潮時(shí)刻發(fā)生溢油時(shí)24 h溢油掃海范圍,由圖6b可知,當(dāng)落潮時(shí)刻發(fā)生溢油時(shí),溢油先向欽州灣外灣方向移動(dòng),等漲潮后返回至欽州港附近海域,受溢油影響區(qū)域主要為欽州灣外灣中部深水航道處。影響區(qū)域主要為深水航道附近,其原因可能與航道處水動(dòng)力條件好,潮流流速大,有利于物質(zhì)的輸運(yùn)有關(guān)。由表2可知,落潮時(shí)刻發(fā)生溢油時(shí),24 h油污擴(kuò)散面積及掃海面積分別約為10.09 km2和50.72 km2,油膜最遠(yuǎn)漂移距離約17.96 km,方向?yàn)闅J州灣外灣方向。相比較漲潮時(shí)刻,溢油落潮時(shí)大部分時(shí)刻擴(kuò)散面積稍大,但掃海面積卻較小,其原因在于落潮時(shí)溢油漂移距離相對(duì)較短,從而造成掃海面積較小。

    3.2.2溢油時(shí)刻對(duì)溢油漂移擴(kuò)散的影響

    圖7為靜風(fēng)條件下不同時(shí)刻發(fā)生溢油時(shí)油污漂移軌跡,由圖7可知,當(dāng)漲潮時(shí)刻發(fā)生溢油時(shí),溢油將隨著潮流向西北偏北方向移動(dòng),由欽州灣灣頸中部進(jìn)入茅尾海;而當(dāng)溢油發(fā)生在落潮時(shí)刻時(shí),溢油則向欽州灣外灣方向移動(dòng)。溢油漂移軌跡的不同,也造成了溢油擴(kuò)散程度及掃海范圍的差異,由表2可知,當(dāng)漲潮時(shí)刻發(fā)生溢油時(shí),24 h油污擴(kuò)散面積及掃海面積分別為9.80 km2和71.83 km2,而落潮時(shí)刻發(fā)生溢油時(shí),這兩個(gè)值分別為10.09 km2和50.72 km2,有明顯的差異。由于不同溢油時(shí)刻下海域水動(dòng)力環(huán)境有一定的差異,因此,其漂移路徑及擴(kuò)散程度也不同,其擴(kuò)散結(jié)果也必將不同。說(shuō)明溢油時(shí)刻對(duì)油膜運(yùn)動(dòng)過(guò)程有顯著的影響,溢油漂移擴(kuò)散結(jié)果與溢油時(shí)刻有密切的關(guān)系。

    圖7不同溢油時(shí)刻下溢油漂移擴(kuò)散軌跡

    Fig.7Drift track of spilled oil under different time

    3.2.3風(fēng)況對(duì)溢油漂移擴(kuò)散的影響

    對(duì)不同風(fēng)況下漲潮時(shí)刻發(fā)生的溢油擴(kuò)散特征進(jìn)行計(jì)算,結(jié)果如表3所示。由表3可知,風(fēng)作用對(duì)欽州灣溢油漂移擴(kuò)散有顯著的影響,不同風(fēng)況作用下的溢油擴(kuò)散效果有明顯的區(qū)別。以靜風(fēng)時(shí)及N風(fēng)時(shí)為例,靜風(fēng)作用下,24 h油污擴(kuò)散面、掃海面積分別為4.15 km2和27.04 km2,而N向風(fēng)存在時(shí),對(duì)應(yīng)的兩個(gè)值分別為9.80 km2和71.83 km2,差異明顯。相比較靜風(fēng)時(shí),N風(fēng)作用下,溢油擴(kuò)散面積及掃海面積更大,其原因在于:N風(fēng)作用下,溢油在24小時(shí)內(nèi)先后進(jìn)入茅尾海和欽州灣的中部,漂移所經(jīng)過(guò)的海域水動(dòng)力相對(duì)較強(qiáng),潮流對(duì)溢油擴(kuò)散的推動(dòng)作用明顯;同時(shí),風(fēng)力的作用也使得溢油向欽州灣外灣方向漂移距離明顯增加,從而造成N風(fēng)作用下溢油擴(kuò)散面積及掃海面積相對(duì)較大的特點(diǎn)。

    表2常風(fēng)向條件下不同時(shí)刻溢油擴(kuò)散面積及掃海面積統(tǒng)計(jì)

    Table 2Statistical results of diffusion area and sweeping area at different time under dominant wind condition

    溢油時(shí)刻O(píng)ilspilltime擴(kuò)散面積(km2)Diffusionarea(km2)掃海面積(km2)Sweepingarea(km2)最遠(yuǎn)漂移距離(km)Farthestdriftdistance(km)6h12h18h24h6h12h18h24h茅尾海方向MaoweiSea欽州灣外灣方向OuterQinzhouBay漲潮Floodtide2.933.9010.879.808.2721.4269.7871.8317.5110.73落潮Ebbtide3.434.737.6510.0911.1819.2833.5950.720.0017.96

    表3不同風(fēng)況條件不同時(shí)刻溢油擴(kuò)散面積及掃海面積統(tǒng)計(jì)

    Table 3Statistical results of diffusion area and sweeping area at different time under different wind conditions

    風(fēng)況Windregime擴(kuò)散面積(km2)Diffusionarea(km2)掃海面積(km2)Sweepingarea(km2)最遠(yuǎn)漂移距離(km)Farthestdriftdistance(km)6h12h18h24h6h12h18h24h茅尾海方向MaoweiSea欽州灣外灣方向OuterQinzhouBay靜風(fēng)Staticwind2.462.744.274.157.8414.2422.0827.0416.436.61N風(fēng)Nwind2.933.9010.879.808.2721.4269.7871.8317.5110.73E風(fēng)Ewind4.113.971.652.8911.1620.9527.3339.8219.036.74S風(fēng)Swind2.291.912.214.436.6612.6217.0922.1415.676.01W風(fēng)Wwind1.901.152.254.125.9810.8515.4319.8014.987.04

    由圖8可知,不同風(fēng)況下,溢油影響區(qū)域有明顯的區(qū)別。除W風(fēng)外,其它各向風(fēng)作用下的溢油均進(jìn)入了茅尾海,S風(fēng)作用下,溢油傾向于進(jìn)入茅尾海東部,而N風(fēng)和E風(fēng)作用時(shí)分別進(jìn)入茅尾海中部及西部。W風(fēng)作用下,溢油則直接漂向欽州港一帶岸線,基本沒(méi)有進(jìn)入茅尾海。由于風(fēng)對(duì)溢油的影響主要在于風(fēng)應(yīng)力對(duì)溢油漂移速度項(xiàng)的貢獻(xiàn),故風(fēng)應(yīng)力的存在會(huì)增大或減小溢油在水體中的漂移速度,從而影響溢油漂移擴(kuò)散的整個(gè)過(guò)程,風(fēng)作用對(duì)溢油漂移擴(kuò)散的影響不容忽視。

    圖8不同風(fēng)況下漲潮時(shí)24 h內(nèi)溢油掃海范圍

    Fig.8Sea swept area of spilled oil under different wind conditions within 24 h

    4結(jié)論

    本文利用MIKE 21 HD模塊建立欽州灣海域水動(dòng)力模型對(duì)欽州灣海域潮流特性進(jìn)行模擬,并在此基礎(chǔ)上利用MIKE 21/3 SA溢油分析模塊,分析欽州灣海上溢油漂移擴(kuò)散的特點(diǎn)及其影響因素,結(jié)果表明1)漲潮時(shí)發(fā)生溢油,溢油先向茅尾海方向漂移,待落潮后退出茅尾海;而落潮時(shí)刻發(fā)生溢油,溢油則先向欽州灣外灣漂移,待漲潮后轉(zhuǎn)向欽州灣灣頸方向。2)漲潮時(shí)發(fā)生溢油,24 h內(nèi)溢油擴(kuò)散面積及掃海面積分別約為9.80 km2和71.83 km2,溢油向茅尾海方向及欽州灣外灣方向漂移的最遠(yuǎn)距離分別約為17.51 km和10.73 km;而落潮時(shí)刻發(fā)生溢油,24 h內(nèi)油污擴(kuò)散面積及掃海面積分別約為10.09 km2和50.72 km2,油膜向欽州灣外灣方向漂移最遠(yuǎn)距離約17.96 km。3)油污擴(kuò)散面積、掃海范圍及遷移路徑與溢油發(fā)生時(shí)刻、風(fēng)況密切相關(guān),不同溢油時(shí)刻及風(fēng)況下,其油污擴(kuò)散面積、掃海范圍、遷移路徑有明顯的差異,溢油時(shí)刻及風(fēng)作用對(duì)溢油漂移擴(kuò)散的影響不容忽視。

    參考文獻(xiàn):

    [1]陳憲云,陸海生,陳波.廣西海岸帶海洋環(huán)境污染現(xiàn)狀及防治對(duì)策[J].廣西科學(xué),2014,21(5):555-560.

    CHEN X Y,LU H S,CHEN B.Status and control countermeasures on marine enuironment pollution of Guangxi[J].Guangxi Sciences,2014,21(5):555-560.

    [2]FAY J A.The spread of oil slicks on a calm sea[C]//Oil on the Sea.New York:Plenum Press, 1969:53-63.

    [3]BLOKKER P C.Spreading and evaporation of petroleum

    products on water[C]//Proceedings of the Fourth International Harbor Conference.Antwerp:Internat Harbour Conf,1964:911-919.

    [4]MACKAY D,MATSUGU R.Evaporation rate of liquid hydrocarbon spills on land and water[J].Can J Chem Eng,1973,51(8):434-439.

    [5]HARRISON W,WINNIK M A,KWONG P T Y,et al.Crude oil spills-Disappearance of aromatic and aliphatic components from small sea-surface slicks[J].Environ Science Tech, 1975,9(3):231-234.

    [6]LEHR W J,CEKIRGE H M,F(xiàn)RAGA R J.Empirical studies of the spreading of oil spills[J].Oil & Petrochemical Pollution,1984,2(1):7-11.

    [7]ASCE Task Committee on Modeling of Oil Spills of the Water Resources Engineering Division.State-of-the-art review of modeling transport and fate of oil spills[J].Journal of Hydraulic Engineering,1996,122(11):594-609.

    [8]ZHANG B,ZHANG C Z,OZER J.SURF-A simulation model for the behavior of oil slicks at sea[C]//Ozer J(ed.).Oil Pollution:Environmental Risk Assessment(OPERA).Proceedings of the OPERA Workshop,Dalian,China,1991:61-85.

    [9]張存智,竇振興,韓康,等.三維溢油動(dòng)態(tài)預(yù)報(bào)模式[J].海洋環(huán)境科學(xué),1997,16(1):22-29.

    ZHANG C Z,DOU Z X,HAN K,et al.A three dimensional model to predict the behavior of oil spills [J].Marine Environmental Science,1997,16(1):22-29.

    [10]楊紅,劉成秀,李曰嵩,等.長(zhǎng)江口南港水道溢油數(shù)值模擬[J].海洋通報(bào),2013,32(3):345-351.

    YANG H,LIU C X,LI Y S,et al.Numerical simulation of oil-spill in the south channel of the stuary Yangtze River Estuary[J].Marine Science Bulletin,2013,32(3):345-351.

    [11]黃娟,曹雅靜,高松,等.渤海海上溢油漂移擴(kuò)散數(shù)值模擬研究[J].海洋科學(xué),2014,38(11):100-107.

    HUANG J,CAO Y J,GAO S,et al.Numerical simulation of oil spill drift -diffusion in the Bohai Sea[J].Marine Sciences,2014,38(11):100-107.

    [12]李懷明,婁安剛,王璟,等.蓬萊19-3油田事故溢油數(shù)值模擬[J].海洋科學(xué),2014,38(6):70-77.

    LI H M,LOU A G,WANG J,et al.A numerical simulation of the oil spill accident at Penglai 19-3 oil field[J].Marine Sciences,2014,38(6):70-77.

    [13]陳義中,朱建榮,陳睿,等.渤海原油碼頭溢油漂移擴(kuò)散的數(shù)值模擬[J].華東師范大學(xué)學(xué)報(bào):自然科學(xué)版,2006(6):25-33.

    CHEN Y Z,ZHU J R,CHEN R,et al.Numerical simulation of the drift and diffusion of oil film at the crude oil wharf in the Bohai Sea[J].Journal of East China Normal University:Natural Science,2006(6):25-33.

    [14]臧士文.基于FVCOM模型的二維海上溢油數(shù)值模擬研究[D].大連:大連理工大學(xué),2011.

    ZANG S W.Study on the Numerical Simulation for Two-Dimensional Oil Spills in Seas Based on FVCOM Model[D].Dalian: Dalian University of Technology,2011.

    [15]袁春光,王義剛,黃惠明,等.大亞灣芒洲島海域移動(dòng)溢油數(shù)值模擬研究[J].海洋通報(bào),2014,33(6):659-667.

    YUAN C G,WANG Y G,HUANG H M,et al.Oil spill numerical simulation on the wharf engineering of Mangzhou Island in the Daya Bay[J].Marine Science Bulletin,2014,33(6):659-667.

    [16]吳曉丹,宋金明,李學(xué)剛,等.海洋溢油油膜厚度影響因素理論模型的構(gòu)建[J].海洋科學(xué),2010,34(2):68-74.

    WU X D,SONG J M,LI X G,et al.Construction of theoretical model for thickness of marine oil spilled [J].Marine Sciences,2010,34(2):68-74.

    (責(zé)任編輯:尹闖)

    Numerical Simulation of Oil Spill in the Qinzhou Bay

    YA Hanzheng,LI Yichun, DONG Dexin

    (Guangxi Key Laboratory of Marine Environmental Science,Guangxi Academy of Sciences,Nanning,Guangxi,530007,China)

    Abstract:【Objective】The drift-diffusion feature and its affecting factors of oil spill in Qinzhou Bay were studied.【Methods】The diffusion characteristics of oil spill in the Qinzhou Bay areas was simulated by MIKE 21/3 SA oil spill model based on the hydrodynamic simulation results.【Results】The oil contamination was drifted towards Maowei Sea at first,and then move out Maowei Sea in ebbing tide when oil spill occurred in flood tide.The farthest drift distance of oil contamination towards Maowei Sea and outer Qinzhou Bay was about 17.51 km and 10.73 km,respectively,and the sweeping area of oil spill was about 71.83 km2 within calculation interval.When oil spill occurred in ebb tide ,the oil stain was drifted towards outer Qinzhou Bay and reversed towards the neck area of Qinzhou Bay in tide rising.In this case,the farthest drift distance of oil contamination towards outer Qinzhou Bay was about 17.96 km,and the sweeping area of oil spill was about 59.72 km2.【Conclusion】The drift-diffusion feature of oil spill in the Qinzhou Bay is closely related to the occurring time and the wind conditions,which cannot be ignored in the oil spill influence evaluation in Qinzhou Bay.

    Key words:Qinzhou Bay,marine oil spill,numerical simulation

    收稿日期:2016-03-11

    作者簡(jiǎn)介:牙韓爭(zhēng)(1988-),男,碩士,主要從事河口海岸物質(zhì)輸運(yùn)研究。

    中圖分類(lèi)號(hào):X55

    文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):1005-9164(2016)02-0156-07

    *國(guó)家自然科學(xué)基金項(xiàng)目(41576024),廣西科學(xué)院基本科研業(yè)務(wù)費(fèi)項(xiàng)目(12YJ25HY08,10YJ25HY03)和廣西近海海洋環(huán)境科學(xué)重點(diǎn)實(shí)驗(yàn)室開(kāi)放基金項(xiàng)目(GXKLHY14-01)資助。

    廣西科學(xué)Guangxi Sciences 2016,23(2):156~162

    網(wǎng)絡(luò)優(yōu)先數(shù)字出版時(shí)間:2016-05-12

    網(wǎng)絡(luò)優(yōu)先數(shù)字出版地址:http://www.cnki.net/kcms/detail/45.1206.G3.20160512.0944.006.html

    猜你喜歡
    落潮漲潮欽州
    漲潮
    北部灣港欽州自動(dòng)化集裝箱碼頭/獲獎(jiǎng)證書(shū)
    西部陸海新通道背景下加強(qiáng)欽州水上消防建設(shè)的思考
    水上消防(2021年4期)2021-11-05 08:51:34
    潮 汐
    “水城”被淹
    欽州坭興陶藝造型設(shè)計(jì)初探
    “欽州大蠔”豪氣崛起
    中資出海并購(gòu)“落潮”
    為什么會(huì)有漲潮和落潮
    神奇的潮汐
    国模一区二区三区四区视频| 亚洲精品影视一区二区三区av| 国产男靠女视频免费网站| 高清在线国产一区| 最近在线观看免费完整版| 欧美一区二区亚洲| 午夜免费观看网址| 午夜福利在线在线| 免费av不卡在线播放| 国产av麻豆久久久久久久| 欧美色欧美亚洲另类二区| 亚洲精品影视一区二区三区av| 亚洲18禁久久av| 午夜福利18| 中文资源天堂在线| 亚洲国产精品久久男人天堂| 精品熟女少妇八av免费久了| 91在线精品国自产拍蜜月 | 免费搜索国产男女视频| 宅男免费午夜| 人人妻,人人澡人人爽秒播| 欧美最新免费一区二区三区 | 国产精华一区二区三区| 深爱激情五月婷婷| www日本黄色视频网| 国产99白浆流出| 特级一级黄色大片| 亚洲人成网站高清观看| 国内精品久久久久久久电影| 天天躁日日操中文字幕| 久久久久久久午夜电影| 他把我摸到了高潮在线观看| 日韩精品青青久久久久久| 一个人看的www免费观看视频| 亚洲在线自拍视频| 婷婷丁香在线五月| 麻豆成人午夜福利视频| 亚洲aⅴ乱码一区二区在线播放| 99在线人妻在线中文字幕| 国产精品香港三级国产av潘金莲| 12—13女人毛片做爰片一| 国产黄片美女视频| 9191精品国产免费久久| 99久久成人亚洲精品观看| 日韩亚洲欧美综合| 最近最新中文字幕大全免费视频| 波多野结衣高清作品| 成人三级黄色视频| 99精品在免费线老司机午夜| 非洲黑人性xxxx精品又粗又长| 免费看十八禁软件| 听说在线观看完整版免费高清| 国产伦精品一区二区三区视频9 | 久久久久久久久久黄片| 在线视频色国产色| ponron亚洲| 国产乱人视频| 757午夜福利合集在线观看| 亚洲在线自拍视频| 欧美另类亚洲清纯唯美| 久久久久性生活片| 亚洲av免费高清在线观看| 久久亚洲真实| 欧美性感艳星| 亚洲熟妇熟女久久| 久久精品国产清高在天天线| 在线视频色国产色| 国产三级中文精品| 成人特级黄色片久久久久久久| 国产免费男女视频| 可以在线观看毛片的网站| 99久久成人亚洲精品观看| 99久久综合精品五月天人人| 狂野欧美白嫩少妇大欣赏| 午夜福利视频1000在线观看| 91久久精品国产一区二区成人 | 好男人在线观看高清免费视频| 久久6这里有精品| 90打野战视频偷拍视频| 国产精品,欧美在线| 精华霜和精华液先用哪个| 久久欧美精品欧美久久欧美| 亚洲成人免费电影在线观看| 日本熟妇午夜| 搡老熟女国产l中国老女人| 黑人欧美特级aaaaaa片| 国产精品 欧美亚洲| 国产单亲对白刺激| 99久久精品国产亚洲精品| 国产淫片久久久久久久久 | 丰满人妻一区二区三区视频av | 丝袜美腿在线中文| 中文字幕熟女人妻在线| 他把我摸到了高潮在线观看| 亚洲av成人不卡在线观看播放网| 国产99白浆流出| 99视频精品全部免费 在线| 日韩精品中文字幕看吧| 无限看片的www在线观看| 小蜜桃在线观看免费完整版高清| 国产成人av激情在线播放| 中文字幕人妻熟人妻熟丝袜美 | 2021天堂中文幕一二区在线观| 全区人妻精品视频| 99久久九九国产精品国产免费| 成人永久免费在线观看视频| 免费在线观看日本一区| 午夜久久久久精精品| 校园春色视频在线观看| 欧美中文日本在线观看视频| 99久久精品国产亚洲精品| 午夜激情福利司机影院| 国产欧美日韩一区二区精品| 中国美女看黄片| 99久久无色码亚洲精品果冻| 69人妻影院| 叶爱在线成人免费视频播放| 国产三级黄色录像| 亚洲精品色激情综合| 国产蜜桃级精品一区二区三区| eeuss影院久久| 亚洲精品日韩av片在线观看 | 九九热线精品视视频播放| 色综合站精品国产| 国产一区二区亚洲精品在线观看| 亚洲精品粉嫩美女一区| 免费高清视频大片| 男女床上黄色一级片免费看| 男女床上黄色一级片免费看| 深夜精品福利| 国产一区二区亚洲精品在线观看| 天天添夜夜摸| 亚洲精品美女久久久久99蜜臀| 亚洲美女视频黄频| 精品久久久久久,| 在线观看免费视频日本深夜| 精品久久久久久久人妻蜜臀av| 最近在线观看免费完整版| 性色avwww在线观看| 超碰av人人做人人爽久久 | 国产亚洲av嫩草精品影院| 亚洲国产精品成人综合色| 亚洲七黄色美女视频| 在线观看美女被高潮喷水网站 | 亚洲电影在线观看av| 免费av不卡在线播放| 又黄又爽又免费观看的视频| 亚洲va日本ⅴa欧美va伊人久久| 中文资源天堂在线| 色综合婷婷激情| 在线播放无遮挡| 国产激情偷乱视频一区二区| 国产综合懂色| 欧美国产日韩亚洲一区| 内地一区二区视频在线| 久久精品国产99精品国产亚洲性色| 老司机午夜十八禁免费视频| 国产成人av教育| 在线免费观看不下载黄p国产 | 动漫黄色视频在线观看| h日本视频在线播放| 一个人看的www免费观看视频| 少妇高潮的动态图| 99久久九九国产精品国产免费| 久久久久国产精品人妻aⅴ院| 熟妇人妻久久中文字幕3abv| 成人永久免费在线观看视频| 亚洲国产日韩欧美精品在线观看 | 丰满人妻一区二区三区视频av | 69av精品久久久久久| 啪啪无遮挡十八禁网站| 久久6这里有精品| 99久久无色码亚洲精品果冻| 国产精品 国内视频| 日本 欧美在线| 狠狠狠狠99中文字幕| 观看美女的网站| 欧美成人性av电影在线观看| 欧美在线一区亚洲| 国产精品久久久久久人妻精品电影| 日本免费a在线| 亚洲aⅴ乱码一区二区在线播放| 男女视频在线观看网站免费| 91久久精品国产一区二区成人 | 色视频www国产| 欧美黄色淫秽网站| 少妇的丰满在线观看| 91麻豆av在线| 亚洲avbb在线观看| 国内久久婷婷六月综合欲色啪| 成人18禁在线播放| xxx96com| 高潮久久久久久久久久久不卡| av专区在线播放| 精品久久久久久久人妻蜜臀av| 亚洲国产欧美网| 国产探花极品一区二区| 成年版毛片免费区| 午夜两性在线视频| 人妻夜夜爽99麻豆av| 精品人妻偷拍中文字幕| 亚洲真实伦在线观看| 丁香欧美五月| 色综合亚洲欧美另类图片| 宅男免费午夜| 国产三级黄色录像| 国内毛片毛片毛片毛片毛片| 欧美成人性av电影在线观看| 日韩高清综合在线| 99精品久久久久人妻精品| 在线免费观看不下载黄p国产 | 少妇丰满av| 欧美黄色片欧美黄色片| 观看美女的网站| 一级黄片播放器| 免费人成在线观看视频色| 欧美日韩综合久久久久久 | 中国美女看黄片| 国产97色在线日韩免费| 亚洲av不卡在线观看| 国产精品99久久久久久久久| 成年女人看的毛片在线观看| 国内揄拍国产精品人妻在线| 丰满乱子伦码专区| av女优亚洲男人天堂| av黄色大香蕉| 色老头精品视频在线观看| 人妻夜夜爽99麻豆av| avwww免费| 日日夜夜操网爽| 亚洲一区二区三区不卡视频| 精品电影一区二区在线| eeuss影院久久| 99热只有精品国产| 51午夜福利影视在线观看| 99精品久久久久人妻精品| 欧美国产日韩亚洲一区| 亚洲一区二区三区不卡视频| 国产高清视频在线播放一区| 黄色女人牲交| 亚洲七黄色美女视频| 老司机午夜十八禁免费视频| 窝窝影院91人妻| 在线观看66精品国产| 色在线成人网| 亚洲精品国产精品久久久不卡| 老鸭窝网址在线观看| 久久久国产成人精品二区| 欧美日韩亚洲国产一区二区在线观看| 天堂网av新在线| 在线观看av片永久免费下载| 欧美性猛交黑人性爽| 国产精华一区二区三区| 国产高清三级在线| 黄色成人免费大全| 黄色女人牲交| 欧美黑人巨大hd| 男女那种视频在线观看| 日韩免费av在线播放| 欧美性猛交╳xxx乱大交人| 欧美成人a在线观看| 狂野欧美白嫩少妇大欣赏| 色尼玛亚洲综合影院| 国产激情欧美一区二区| 熟妇人妻久久中文字幕3abv| 特级一级黄色大片| 欧美三级亚洲精品| av中文乱码字幕在线| 国产主播在线观看一区二区| 88av欧美| 久久久精品欧美日韩精品| 国产麻豆成人av免费视频| 亚洲第一欧美日韩一区二区三区| 尤物成人国产欧美一区二区三区| 国产色婷婷99| 亚洲无线观看免费| 老熟妇乱子伦视频在线观看| 亚洲电影在线观看av| 国产精品影院久久| 精品国产三级普通话版| 亚洲无线在线观看| 午夜老司机福利剧场| 两人在一起打扑克的视频| 99久久99久久久精品蜜桃| 亚洲国产高清在线一区二区三| 国产午夜精品久久久久久一区二区三区 | 又爽又黄无遮挡网站| 免费在线观看日本一区| 亚洲专区国产一区二区| 亚洲国产精品sss在线观看| 国产一区二区激情短视频| 一本久久中文字幕| 久久国产精品人妻蜜桃| 亚洲av二区三区四区| 男女下面进入的视频免费午夜| 久久久国产精品麻豆| 亚洲人与动物交配视频| 制服人妻中文乱码| 国内精品一区二区在线观看| 久久精品人妻少妇| 久久中文看片网| 久久亚洲真实| 成人鲁丝片一二三区免费| 中文字幕久久专区| 亚洲精品一卡2卡三卡4卡5卡| 波多野结衣高清无吗| 岛国视频午夜一区免费看| 最近最新中文字幕大全免费视频| 国产精品美女特级片免费视频播放器| 高清在线国产一区| 又黄又爽又免费观看的视频| 亚洲乱码一区二区免费版| 亚洲av日韩精品久久久久久密| 男人和女人高潮做爰伦理| 久久久久精品国产欧美久久久| 757午夜福利合集在线观看| 国产av在哪里看| 丰满乱子伦码专区| 亚洲国产色片| 国产成人a区在线观看| 在线a可以看的网站| 一a级毛片在线观看| 国产伦人伦偷精品视频| aaaaa片日本免费| 午夜视频国产福利| 国产高清视频在线播放一区| 18禁国产床啪视频网站| 久久6这里有精品| 中出人妻视频一区二区| 香蕉久久夜色| 精品人妻1区二区| 露出奶头的视频| 亚洲人与动物交配视频| 啦啦啦韩国在线观看视频| 国产一区二区三区视频了| 丰满人妻一区二区三区视频av | 亚洲第一欧美日韩一区二区三区| 在线观看舔阴道视频| 国产精品影院久久| 久久久久久大精品| 狂野欧美白嫩少妇大欣赏| 天堂av国产一区二区熟女人妻| 99久久久亚洲精品蜜臀av| 深夜精品福利| 欧美成狂野欧美在线观看| 黑人欧美特级aaaaaa片| 全区人妻精品视频| 亚洲av美国av| 久久久久久久午夜电影| 在线观看免费午夜福利视频| 免费观看精品视频网站| 午夜福利在线观看免费完整高清在 | 精品欧美国产一区二区三| 国产成人a区在线观看| 国产极品精品免费视频能看的| 日韩欧美三级三区| 国产精华一区二区三区| 国产精品一区二区免费欧美| 欧美三级亚洲精品| 亚洲精品色激情综合| 久久久色成人| 99热精品在线国产| 精品人妻偷拍中文字幕| 久久精品综合一区二区三区| 熟女人妻精品中文字幕| 看片在线看免费视频| 久久欧美精品欧美久久欧美| 最近视频中文字幕2019在线8| 亚洲欧美激情综合另类| 搡老熟女国产l中国老女人| 久久人人精品亚洲av| 香蕉丝袜av| 亚洲欧美精品综合久久99| 在线十欧美十亚洲十日本专区| av国产免费在线观看| 母亲3免费完整高清在线观看| 国内久久婷婷六月综合欲色啪| 日本成人三级电影网站| 黄片小视频在线播放| 好男人在线观看高清免费视频| 窝窝影院91人妻| 婷婷精品国产亚洲av在线| 免费观看人在逋| 人妻久久中文字幕网| 一个人看的www免费观看视频| 精品国产美女av久久久久小说| 桃红色精品国产亚洲av| 高潮久久久久久久久久久不卡| 国产爱豆传媒在线观看| 噜噜噜噜噜久久久久久91| aaaaa片日本免费| 日韩国内少妇激情av| 久久这里只有精品中国| 午夜老司机福利剧场| 人妻丰满熟妇av一区二区三区| 日韩免费av在线播放| 乱人视频在线观看| 91麻豆av在线| 小说图片视频综合网站| 亚洲熟妇中文字幕五十中出| bbb黄色大片| 在线观看舔阴道视频| 啦啦啦免费观看视频1| 最近最新免费中文字幕在线| xxx96com| 伊人久久大香线蕉亚洲五| 超碰av人人做人人爽久久 | 久久精品国产亚洲av涩爱 | 成人欧美大片| 亚洲美女黄片视频| 亚洲国产精品成人综合色| 国产综合懂色| 精品人妻偷拍中文字幕| 我要搜黄色片| 高清在线国产一区| 99久久久亚洲精品蜜臀av| av在线天堂中文字幕| 99热这里只有精品一区| 欧美日韩综合久久久久久 | 国产日本99.免费观看| 欧美激情久久久久久爽电影| 亚洲激情在线av| 成人av在线播放网站| 亚洲无线观看免费| 亚洲片人在线观看| 欧美成狂野欧美在线观看| 欧美高清成人免费视频www| 午夜福利欧美成人| 久久九九热精品免费| 99久久99久久久精品蜜桃| 午夜福利欧美成人| 免费观看的影片在线观看| 国产伦精品一区二区三区四那| 亚洲欧美日韩无卡精品| 日韩高清综合在线| 中文字幕高清在线视频| 中文字幕人成人乱码亚洲影| 日本一二三区视频观看| 18+在线观看网站| 在线观看一区二区三区| tocl精华| 欧美乱色亚洲激情| 国产91精品成人一区二区三区| 国产麻豆成人av免费视频| 婷婷亚洲欧美| 国产欧美日韩精品亚洲av| 国产成人欧美在线观看| 日本黄色视频三级网站网址| 搡老妇女老女人老熟妇| 日本三级黄在线观看| 久久人人精品亚洲av| 亚洲精品久久国产高清桃花| 99久久九九国产精品国产免费| ponron亚洲| 国产高清有码在线观看视频| 中文字幕高清在线视频| 国产私拍福利视频在线观看| 国产精品永久免费网站| 日韩欧美精品v在线| 热99re8久久精品国产| 熟女电影av网| 国产精华一区二区三区| 波野结衣二区三区在线 | 国产97色在线日韩免费| av在线蜜桃| 亚洲av成人不卡在线观看播放网| 狠狠狠狠99中文字幕| 成熟少妇高潮喷水视频| 久久久久久人人人人人| 免费观看人在逋| 在线看三级毛片| 亚洲成人精品中文字幕电影| 99热6这里只有精品| 精品一区二区三区视频在线观看免费| 老鸭窝网址在线观看| 天天一区二区日本电影三级| 一夜夜www| 欧美日韩国产亚洲二区| 国产一区二区亚洲精品在线观看| 最近最新中文字幕大全电影3| 无遮挡黄片免费观看| 别揉我奶头~嗯~啊~动态视频| 成人国产综合亚洲| 久久久国产精品麻豆| 丰满人妻熟妇乱又伦精品不卡| 人人妻,人人澡人人爽秒播| a在线观看视频网站| 内射极品少妇av片p| АⅤ资源中文在线天堂| 国产精品久久电影中文字幕| 日韩免费av在线播放| 一区二区三区免费毛片| 欧美乱码精品一区二区三区| 亚洲精品色激情综合| 色av中文字幕| 亚洲av成人不卡在线观看播放网| 一边摸一边抽搐一进一小说| 国产午夜精品论理片| АⅤ资源中文在线天堂| 欧洲精品卡2卡3卡4卡5卡区| 精品人妻1区二区| 午夜精品在线福利| 成年版毛片免费区| 国产精品精品国产色婷婷| 亚洲在线观看片| 在线观看一区二区三区| 亚洲人成伊人成综合网2020| 成年女人看的毛片在线观看| 久久久久久久精品吃奶| 亚洲午夜理论影院| 亚洲精品美女久久久久99蜜臀| 内地一区二区视频在线| 男女午夜视频在线观看| 亚洲国产色片| 中文亚洲av片在线观看爽| 性色av乱码一区二区三区2| 99久久精品一区二区三区| 久久精品国产清高在天天线| 精品午夜福利视频在线观看一区| 午夜福利高清视频| 欧美zozozo另类| 午夜日韩欧美国产| 美女高潮喷水抽搐中文字幕| 99久久综合精品五月天人人| 国产精品影院久久| 久9热在线精品视频| 欧美区成人在线视频| 亚洲黑人精品在线| 一个人看视频在线观看www免费 | 淫妇啪啪啪对白视频| 国产一区在线观看成人免费| 亚洲成a人片在线一区二区| 免费在线观看影片大全网站| 一边摸一边抽搐一进一小说| 久久久久国产精品人妻aⅴ院| 亚洲性夜色夜夜综合| 精品国内亚洲2022精品成人| 18禁国产床啪视频网站| 精品99又大又爽又粗少妇毛片 | 99久久综合精品五月天人人| 最好的美女福利视频网| 亚洲人与动物交配视频| 免费在线观看影片大全网站| 男人的好看免费观看在线视频| av福利片在线观看| 久久亚洲真实| 又黄又爽又免费观看的视频| 日本免费a在线| 韩国av一区二区三区四区| 成年女人永久免费观看视频| 人妻丰满熟妇av一区二区三区| 久久中文看片网| 中文字幕人妻熟人妻熟丝袜美 | 成熟少妇高潮喷水视频| 在线视频色国产色| 日本黄色片子视频| 国产三级在线视频| 亚洲成a人片在线一区二区| 欧美又色又爽又黄视频| 国产美女午夜福利| 一进一出抽搐gif免费好疼| 国产欧美日韩一区二区精品| 日本黄大片高清| 亚洲欧美日韩高清专用| 国内精品美女久久久久久| 一本一本综合久久| 麻豆久久精品国产亚洲av| 久久午夜亚洲精品久久| 日本成人三级电影网站| www.www免费av| 欧美日韩综合久久久久久 | 亚洲电影在线观看av| 禁无遮挡网站| 麻豆久久精品国产亚洲av| 国产欧美日韩一区二区精品| 成人av一区二区三区在线看| xxxwww97欧美| 美女大奶头视频| 中出人妻视频一区二区| 日本 av在线| 国产在视频线在精品| 丰满人妻一区二区三区视频av | 九九在线视频观看精品| 一本久久中文字幕| 亚洲一区二区三区色噜噜| 麻豆成人午夜福利视频| 男女午夜视频在线观看| 美女 人体艺术 gogo| 国产高潮美女av| 午夜精品在线福利| 小说图片视频综合网站| 日本 欧美在线| 99国产精品一区二区三区| 宅男免费午夜| 国产单亲对白刺激| 国产三级在线视频| 欧美成狂野欧美在线观看| av专区在线播放| 亚洲人成网站在线播放欧美日韩| 午夜精品在线福利| 午夜福利高清视频| 俺也久久电影网| 国产免费一级a男人的天堂| 欧美日韩福利视频一区二区| 国产精品1区2区在线观看.| 一夜夜www| 欧美日韩亚洲国产一区二区在线观看| 欧美又色又爽又黄视频| 18禁裸乳无遮挡免费网站照片| 国产精品99久久久久久久久| 嫩草影院精品99| 久久人妻av系列| 国产欧美日韩一区二区三| 国产麻豆成人av免费视频| 日韩欧美 国产精品|