• <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ì)有漲潮和落潮
    神奇的潮汐
    超碰97精品在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产在线精品亚洲第一网站| e午夜精品久久久久久久| 热99re8久久精品国产| 18在线观看网站| 亚洲一码二码三码区别大吗| 大陆偷拍与自拍| 欧美日韩中文字幕国产精品一区二区三区 | 欧美性长视频在线观看| 色综合婷婷激情| 久久国产精品人妻蜜桃| 国产亚洲欧美精品永久| 日韩一卡2卡3卡4卡2021年| 国产欧美日韩综合在线一区二区| 18在线观看网站| 久久九九热精品免费| cao死你这个sao货| 黑丝袜美女国产一区| 国产成人精品在线电影| 国产午夜精品久久久久久| 91精品三级在线观看| 亚洲精品美女久久久久99蜜臀| 男女床上黄色一级片免费看| 91精品三级在线观看| 国产一区二区激情短视频| 在线天堂中文资源库| 国产精品九九99| 亚洲熟妇熟女久久| 国产免费现黄频在线看| 国产男女超爽视频在线观看| 成人免费观看视频高清| 国产一区二区激情短视频| 天天添夜夜摸| 日韩欧美免费精品| 国产精品.久久久| 91成年电影在线观看| 国产一区二区三区在线臀色熟女 | 黑人巨大精品欧美一区二区蜜桃| av线在线观看网站| 在线观看一区二区三区激情| 精品亚洲成a人片在线观看| 亚洲少妇的诱惑av| 涩涩av久久男人的天堂| 男女午夜视频在线观看| av天堂在线播放| 天天影视国产精品| 久久久精品免费免费高清| 热99久久久久精品小说推荐| 丰满的人妻完整版| 一本一本久久a久久精品综合妖精| 19禁男女啪啪无遮挡网站| 成人亚洲精品一区在线观看| 天堂中文最新版在线下载| 日韩欧美免费精品| 免费女性裸体啪啪无遮挡网站| 国产精品久久久久久人妻精品电影| 大陆偷拍与自拍| 成人亚洲精品一区在线观看| 国产av精品麻豆| 久久中文字幕一级| 亚洲欧美一区二区三区黑人| 精品国产一区二区三区久久久樱花| 90打野战视频偷拍视频| 男女之事视频高清在线观看| 亚洲七黄色美女视频| 精品国产乱码久久久久久男人| 中文字幕人妻丝袜制服| 一夜夜www| 视频在线观看一区二区三区| 国产黄色免费在线视频| 香蕉国产在线看| 一进一出抽搐动态| 午夜精品久久久久久毛片777| 欧美日韩一级在线毛片| 制服诱惑二区| 国产成人av教育| 大码成人一级视频| 无人区码免费观看不卡| 久久久精品国产亚洲av高清涩受| 很黄的视频免费| 精品一区二区三区四区五区乱码| 18禁裸乳无遮挡动漫免费视频| 精品欧美一区二区三区在线| 一级片'在线观看视频| 亚洲精品国产一区二区精华液| 夫妻午夜视频| 欧美久久黑人一区二区| 丁香六月欧美| 一边摸一边抽搐一进一出视频| 一边摸一边做爽爽视频免费| 中文字幕精品免费在线观看视频| 精品高清国产在线一区| 久久精品国产综合久久久| 9191精品国产免费久久| 超碰成人久久| 777久久人妻少妇嫩草av网站| 丝袜在线中文字幕| 婷婷精品国产亚洲av在线 | 这个男人来自地球电影免费观看| 午夜精品在线福利| 精品熟女少妇八av免费久了| 免费在线观看亚洲国产| 女人精品久久久久毛片| 成人三级做爰电影| 中文字幕最新亚洲高清| 亚洲情色 制服丝袜| 在线观看免费高清a一片| 在线国产一区二区在线| 窝窝影院91人妻| netflix在线观看网站| 欧美激情久久久久久爽电影 | 母亲3免费完整高清在线观看| 国产色视频综合| 无人区码免费观看不卡| 欧美中文综合在线视频| 亚洲精品av麻豆狂野| 精品高清国产在线一区| 亚洲午夜精品一区,二区,三区| 国产视频一区二区在线看| 12—13女人毛片做爰片一| 国产精品久久久av美女十八| 91麻豆精品激情在线观看国产 | 另类亚洲欧美激情| 大香蕉久久网| 午夜福利在线观看吧| 在线观看免费午夜福利视频| 国产成人欧美在线观看 | av在线播放免费不卡| 人妻 亚洲 视频| 黄色女人牲交| 成在线人永久免费视频| 亚洲国产欧美一区二区综合| 别揉我奶头~嗯~啊~动态视频| 首页视频小说图片口味搜索| 亚洲精品国产色婷婷电影| 真人做人爱边吃奶动态| 国产免费男女视频| 999久久久精品免费观看国产| 久久精品熟女亚洲av麻豆精品| 久久久国产精品麻豆| 一级片免费观看大全| 很黄的视频免费| 国产精品一区二区在线不卡| 成人永久免费在线观看视频| 午夜两性在线视频| 亚洲国产精品sss在线观看 | 亚洲欧美日韩另类电影网站| 成人特级黄色片久久久久久久| 精品电影一区二区在线| 日本欧美视频一区| 亚洲熟妇中文字幕五十中出 | 久久精品人人爽人人爽视色| 久久人妻熟女aⅴ| 日韩大码丰满熟妇| 免费在线观看亚洲国产| 伦理电影免费视频| 免费日韩欧美在线观看| 一级a爱视频在线免费观看| 纯流量卡能插随身wifi吗| 久久国产精品男人的天堂亚洲| 亚洲精品美女久久av网站| 国产精品久久久久久人妻精品电影| 别揉我奶头~嗯~啊~动态视频| 国产精品 欧美亚洲| 午夜91福利影院| 亚洲精品中文字幕在线视频| 亚洲成人免费av在线播放| 午夜福利,免费看| 男女之事视频高清在线观看| 日本一区二区免费在线视频| 高清av免费在线| 色播在线永久视频| 999久久久精品免费观看国产| 国产又色又爽无遮挡免费看| 亚洲av电影在线进入| 青草久久国产| 丰满人妻熟妇乱又伦精品不卡| 韩国av一区二区三区四区| 亚洲专区中文字幕在线| 亚洲第一青青草原| 亚洲黑人精品在线| 国产97色在线日韩免费| 色综合欧美亚洲国产小说| 99久久99久久久精品蜜桃| 超碰97精品在线观看| 久久精品aⅴ一区二区三区四区| 99国产精品一区二区蜜桃av | 电影成人av| 色综合欧美亚洲国产小说| 久久影院123| 又紧又爽又黄一区二区| 国产日韩欧美亚洲二区| 19禁男女啪啪无遮挡网站| 国产极品粉嫩免费观看在线| 亚洲专区国产一区二区| 国产亚洲精品第一综合不卡| 90打野战视频偷拍视频| 午夜精品久久久久久毛片777| 人人妻人人添人人爽欧美一区卜| 亚洲精品成人av观看孕妇| 黄色女人牲交| 亚洲一码二码三码区别大吗| 精品午夜福利视频在线观看一区| 热99国产精品久久久久久7| 国产成人系列免费观看| 国产成人精品在线电影| 桃红色精品国产亚洲av| 两个人看的免费小视频| 不卡一级毛片| 捣出白浆h1v1| 精品国产亚洲在线| 老汉色∧v一级毛片| 午夜免费观看网址| 久久精品成人免费网站| 伦理电影免费视频| 欧美日韩精品网址| 国产在线观看jvid| 日韩欧美免费精品| 日韩一卡2卡3卡4卡2021年| 啦啦啦视频在线资源免费观看| 欧美成人午夜精品| 80岁老熟妇乱子伦牲交| 欧美日韩亚洲综合一区二区三区_| 欧美成狂野欧美在线观看| 国产aⅴ精品一区二区三区波| 精品一区二区三区视频在线观看免费 | 成在线人永久免费视频| 亚洲一区二区三区欧美精品| 中文字幕色久视频| 老司机深夜福利视频在线观看| 91成人精品电影| 日韩视频一区二区在线观看| 中文字幕人妻熟女乱码| 丰满人妻熟妇乱又伦精品不卡| 一本综合久久免费| 极品少妇高潮喷水抽搐| 亚洲av日韩精品久久久久久密| 操出白浆在线播放| 日日摸夜夜添夜夜添小说| 国产一区二区三区视频了| 午夜免费观看网址| 一夜夜www| 自拍欧美九色日韩亚洲蝌蚪91| 欧美激情久久久久久爽电影 | 校园春色视频在线观看| 高清欧美精品videossex| 亚洲av成人不卡在线观看播放网| 天堂√8在线中文| 国产精品偷伦视频观看了| 色尼玛亚洲综合影院| 在线观看66精品国产| 国产淫语在线视频| 大型av网站在线播放| 91麻豆av在线| 国产免费男女视频| 欧美 日韩 精品 国产| av在线播放免费不卡| 国产有黄有色有爽视频| 日本一区二区免费在线视频| 啦啦啦视频在线资源免费观看| 亚洲全国av大片| 免费高清在线观看日韩| 久久久精品国产亚洲av高清涩受| 在线观看舔阴道视频| 新久久久久国产一级毛片| 日韩视频一区二区在线观看| 免费在线观看亚洲国产| 欧美日韩国产mv在线观看视频| 天天躁夜夜躁狠狠躁躁| 露出奶头的视频| 免费不卡黄色视频| 又黄又粗又硬又大视频| 日韩精品免费视频一区二区三区| 亚洲五月天丁香| 免费在线观看日本一区| 国产高清视频在线播放一区| 中文字幕人妻丝袜一区二区| 久久精品亚洲av国产电影网| 脱女人内裤的视频| 身体一侧抽搐| 女人久久www免费人成看片| x7x7x7水蜜桃| 91九色精品人成在线观看| 亚洲少妇的诱惑av| 少妇裸体淫交视频免费看高清 | 又黄又粗又硬又大视频| 一区二区三区激情视频| a级片在线免费高清观看视频| 亚洲精华国产精华精| 久久午夜综合久久蜜桃| 三上悠亚av全集在线观看| 露出奶头的视频| x7x7x7水蜜桃| 久久精品亚洲精品国产色婷小说| 午夜福利欧美成人| 欧美在线一区亚洲| 国内久久婷婷六月综合欲色啪| 久久国产精品影院| 18禁国产床啪视频网站| a在线观看视频网站| 精品国内亚洲2022精品成人 | 欧美av亚洲av综合av国产av| 欧美激情高清一区二区三区| 成熟少妇高潮喷水视频| 国产精品98久久久久久宅男小说| 久久久久精品人妻al黑| 欧美丝袜亚洲另类 | 日本撒尿小便嘘嘘汇集6| 亚洲国产欧美网| 十八禁高潮呻吟视频| 中文字幕人妻丝袜制服| 大码成人一级视频| 女性被躁到高潮视频| 午夜福利一区二区在线看| 欧美亚洲 丝袜 人妻 在线| 怎么达到女性高潮| 久久人人97超碰香蕉20202| 国产男女超爽视频在线观看| 久久精品国产99精品国产亚洲性色 | 91老司机精品| 每晚都被弄得嗷嗷叫到高潮| 老司机影院毛片| 亚洲精品国产一区二区精华液| 精品国产亚洲在线| 俄罗斯特黄特色一大片| 欧美成人免费av一区二区三区 | 五月开心婷婷网| 国产97色在线日韩免费| 午夜激情av网站| a级片在线免费高清观看视频| 女警被强在线播放| 欧美 亚洲 国产 日韩一| 国产高清激情床上av| 69精品国产乱码久久久| 成人永久免费在线观看视频| 国产高清激情床上av| 日本一区二区免费在线视频| 免费日韩欧美在线观看| 欧美老熟妇乱子伦牲交| 久久午夜亚洲精品久久| 国精品久久久久久国模美| 村上凉子中文字幕在线| 国产精品久久久久久人妻精品电影| 精品第一国产精品| 乱人伦中国视频| 啪啪无遮挡十八禁网站| 久久久久精品国产欧美久久久| 天天躁夜夜躁狠狠躁躁| 国产成人精品无人区| 亚洲欧美日韩高清在线视频| 黑人巨大精品欧美一区二区蜜桃| 18禁美女被吸乳视频| 久久久久久久久免费视频了| 丝瓜视频免费看黄片| 狠狠狠狠99中文字幕| 久久人妻av系列| 女人高潮潮喷娇喘18禁视频| 久久久久久久久久久久大奶| 精品少妇久久久久久888优播| 啦啦啦免费观看视频1| tocl精华| 亚洲国产精品合色在线| 热99久久久久精品小说推荐| 无人区码免费观看不卡| 黄色成人免费大全| 女人久久www免费人成看片| 老司机午夜福利在线观看视频| 欧美日韩黄片免| 高清黄色对白视频在线免费看| 怎么达到女性高潮| 日韩 欧美 亚洲 中文字幕| 免费观看a级毛片全部| 好看av亚洲va欧美ⅴa在| 黄色女人牲交| 久久久国产成人免费| a级毛片黄视频| 少妇被粗大的猛进出69影院| 精品久久久久久久毛片微露脸| 国产精品二区激情视频| 亚洲精品在线美女| 精品久久久精品久久久| 午夜免费鲁丝| 久久99一区二区三区| 嫁个100分男人电影在线观看| 91老司机精品| 一边摸一边做爽爽视频免费| 国产一区二区三区在线臀色熟女 | 色综合欧美亚洲国产小说| 在线观看www视频免费| 黄色视频,在线免费观看| 国产精品一区二区免费欧美| 亚洲国产毛片av蜜桃av| 桃红色精品国产亚洲av| 国产精品久久久久成人av| 日本黄色日本黄色录像| 亚洲avbb在线观看| 国产一区二区三区视频了| 亚洲精品中文字幕一二三四区| 免费在线观看日本一区| 精品一区二区三卡| 99热国产这里只有精品6| 在线十欧美十亚洲十日本专区| 91精品三级在线观看| 大码成人一级视频| 美女视频免费永久观看网站| 欧美日韩亚洲综合一区二区三区_| 欧美大码av| 国产精品久久久av美女十八| 在线观看免费日韩欧美大片| 一本综合久久免费| 欧美日本中文国产一区发布| 色尼玛亚洲综合影院| 91字幕亚洲| 丝袜人妻中文字幕| 中文字幕精品免费在线观看视频| 国产精品98久久久久久宅男小说| 日韩熟女老妇一区二区性免费视频| 熟女少妇亚洲综合色aaa.| 老司机靠b影院| 最近最新中文字幕大全免费视频| 脱女人内裤的视频| 精品一区二区三卡| 无人区码免费观看不卡| 变态另类成人亚洲欧美熟女 | xxxhd国产人妻xxx| 精品亚洲成a人片在线观看| 18禁裸乳无遮挡免费网站照片 | 亚洲中文av在线| 操出白浆在线播放| 老熟妇乱子伦视频在线观看| 亚洲av成人av| 欧洲精品卡2卡3卡4卡5卡区| 亚洲 国产 在线| 欧美亚洲 丝袜 人妻 在线| 国产亚洲精品一区二区www | 日韩欧美在线二视频 | 国产欧美日韩综合在线一区二区| 国产精品1区2区在线观看. | 操出白浆在线播放| 亚洲精品一卡2卡三卡4卡5卡| 大陆偷拍与自拍| 精品一区二区三区av网在线观看| tube8黄色片| 99热国产这里只有精品6| 亚洲成人免费电影在线观看| 久久精品国产亚洲av高清一级| 日日摸夜夜添夜夜添小说| 国产精品国产av在线观看| 免费不卡黄色视频| 悠悠久久av| 丝袜人妻中文字幕| 在线十欧美十亚洲十日本专区| 大陆偷拍与自拍| 窝窝影院91人妻| 国产欧美日韩一区二区三| 久久草成人影院| av天堂在线播放| 久久久国产成人免费| 国产视频一区二区在线看| 天天操日日干夜夜撸| 亚洲专区字幕在线| 色94色欧美一区二区| 国产成人av激情在线播放| 国产精品香港三级国产av潘金莲| 亚洲免费av在线视频| 自线自在国产av| 欧美日韩一级在线毛片| 9热在线视频观看99| 一夜夜www| 午夜91福利影院| 捣出白浆h1v1| 国产精品一区二区免费欧美| 欧美黑人精品巨大| 可以免费在线观看a视频的电影网站| 国产深夜福利视频在线观看| 久99久视频精品免费| 亚洲精品国产区一区二| 亚洲五月天丁香| 757午夜福利合集在线观看| 国产精品自产拍在线观看55亚洲 | 亚洲成人国产一区在线观看| 欧美另类亚洲清纯唯美| 亚洲美女黄片视频| 黄色片一级片一级黄色片| 亚洲精品中文字幕在线视频| 日韩大码丰满熟妇| 久热爱精品视频在线9| 亚洲av第一区精品v没综合| 女人爽到高潮嗷嗷叫在线视频| 日韩熟女老妇一区二区性免费视频| 黄色a级毛片大全视频| 欧美老熟妇乱子伦牲交| 国产精品欧美亚洲77777| 国产xxxxx性猛交| 亚洲精品乱久久久久久| 女人高潮潮喷娇喘18禁视频| 国产成人av教育| 18禁观看日本| 免费观看a级毛片全部| 国产亚洲精品久久久久久毛片 | 超碰成人久久| 黄频高清免费视频| 一本综合久久免费| 久久久久久久午夜电影 | 香蕉久久夜色| 成在线人永久免费视频| 久久精品国产综合久久久| 亚洲av电影在线进入| 色综合婷婷激情| 五月开心婷婷网| 美女福利国产在线| av有码第一页| 真人做人爱边吃奶动态| 一二三四社区在线视频社区8| 夜夜躁狠狠躁天天躁| 亚洲久久久国产精品| 亚洲欧美日韩高清在线视频| 国产精品二区激情视频| 久久热在线av| 人妻一区二区av| 国产男女内射视频| 999久久久国产精品视频| 嫁个100分男人电影在线观看| 亚洲成人免费电影在线观看| 国产精品久久久久久人妻精品电影| 激情在线观看视频在线高清 | 国产精品 欧美亚洲| 国产在线观看jvid| 女人精品久久久久毛片| 精品国产超薄肉色丝袜足j| 久久人妻熟女aⅴ| 中文字幕精品免费在线观看视频| 国产av精品麻豆| 精品高清国产在线一区| 亚洲精品粉嫩美女一区| 黄色毛片三级朝国网站| 亚洲av日韩精品久久久久久密| 中文字幕av电影在线播放| 国产午夜精品久久久久久| 精品国内亚洲2022精品成人 | 午夜精品在线福利| 国产一区二区三区视频了| 一区在线观看完整版| 免费在线观看影片大全网站| 后天国语完整版免费观看| 欧美成人午夜精品| 最近最新免费中文字幕在线| 侵犯人妻中文字幕一二三四区| 国产成人精品无人区| 欧美日韩精品网址| 亚洲五月天丁香| 变态另类成人亚洲欧美熟女 | 成人国产一区最新在线观看| 精品一区二区三区四区五区乱码| av天堂在线播放| 亚洲成a人片在线一区二区| 亚洲美女黄片视频| av片东京热男人的天堂| 极品少妇高潮喷水抽搐| 亚洲成国产人片在线观看| 日本欧美视频一区| 日韩欧美在线二视频 | 精品一区二区三区四区五区乱码| 18禁观看日本| 操出白浆在线播放| 久久狼人影院| 欧美黄色淫秽网站| 村上凉子中文字幕在线| 99热网站在线观看| 亚洲全国av大片| 视频区图区小说| 天天躁夜夜躁狠狠躁躁| 成人黄色视频免费在线看| 国产一区二区三区视频了| 嫁个100分男人电影在线观看| 侵犯人妻中文字幕一二三四区| 最近最新中文字幕大全免费视频| 欧美 亚洲 国产 日韩一| 人人妻人人爽人人添夜夜欢视频| 中文字幕最新亚洲高清| 国产精品香港三级国产av潘金莲| 十分钟在线观看高清视频www| 国产男靠女视频免费网站| 日本vs欧美在线观看视频| 黑人巨大精品欧美一区二区蜜桃| 在线观看免费午夜福利视频| 免费在线观看影片大全网站| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美 日韩 精品 国产| 操出白浆在线播放| 男人舔女人的私密视频| 无人区码免费观看不卡| 久久人人97超碰香蕉20202| 精品无人区乱码1区二区| 纯流量卡能插随身wifi吗| 中文字幕人妻熟女乱码| 天堂俺去俺来也www色官网| 两性夫妻黄色片| 久久婷婷成人综合色麻豆| 一区二区三区精品91| 一边摸一边做爽爽视频免费| 国产黄色免费在线视频| 亚洲精品国产区一区二| 一级作爱视频免费观看| 女人被躁到高潮嗷嗷叫费观| 久久香蕉激情| 成人国产一区最新在线观看| x7x7x7水蜜桃| av天堂在线播放| 性少妇av在线| 亚洲美女黄片视频| 婷婷丁香在线五月| 80岁老熟妇乱子伦牲交| 最近最新中文字幕大全电影3 |