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

    1522號(hào)臺(tái)風(fēng)“彩虹”在雷州半島東部海域的風(fēng)暴潮研究

    2017-12-27 02:16:35熊夢(mèng)婕章衛(wèi)勝張金善殷成團(tuán)
    海洋預(yù)報(bào) 2017年6期
    關(guān)鍵詞:模型

    熊夢(mèng)婕,章衛(wèi)勝,張金善,殷成團(tuán)

    (1.南京水利科學(xué)研究院,江蘇南京210024;2.水文水資源與水利工程國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇南京210024)

    1522號(hào)臺(tái)風(fēng)“彩虹”在雷州半島東部海域的風(fēng)暴潮研究

    熊夢(mèng)婕1,2,章衛(wèi)勝1,2,張金善1,2,殷成團(tuán)1,2

    (1.南京水利科學(xué)研究院,江蘇南京210024;2.水文水資源與水利工程國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇南京210024)

    建立了三重嵌套的高分辨率二維風(fēng)暴潮與天文潮耦合數(shù)值模型,研究了雷州半島東部海域風(fēng)暴增水的時(shí)空分布規(guī)律及其產(chǎn)生原因。模擬結(jié)果表明:1522號(hào)強(qiáng)臺(tái)風(fēng)“彩虹”在雷州半島東部海域引起的風(fēng)暴增水呈現(xiàn)由南向北逐漸增加的趨勢(shì),分布大致以臺(tái)風(fēng)路徑為界:臺(tái)風(fēng)路徑右側(cè)的南三水道、鑒江等地發(fā)生2.8 m以上的強(qiáng)烈增水,向岸東南風(fēng)的持續(xù)影響與向西凹進(jìn)的海灣地形是引起當(dāng)?shù)仫L(fēng)暴增水的兩個(gè)主要因素;臺(tái)風(fēng)路徑左側(cè)的雷州灣西岸則發(fā)生了顯著的減水現(xiàn)象,南渡站的實(shí)測(cè)最大減水值達(dá)到2.11 m,與臺(tái)風(fēng)登陸前后當(dāng)?shù)厥茈x岸風(fēng)的作用有關(guān)。

    臺(tái)風(fēng);雷州半島;數(shù)值模擬;風(fēng)暴增水;風(fēng)暴減水

    1 引言

    雷州半島地處中國(guó)大陸最南端,西臨北部灣,南與海南島隔瓊州海峽相望,夏秋季熱帶氣旋活動(dòng)頻繁,是我國(guó)遭受臺(tái)風(fēng)暴潮災(zāi)害最嚴(yán)重的岸段之一。其特殊的地理位置以及向西凹進(jìn)的彎曲型海岸形態(tài)十分有利于風(fēng)暴增水的發(fā)展。

    據(jù)統(tǒng)計(jì),1949—2005年間,粵西地區(qū)共發(fā)生風(fēng)暴潮災(zāi)害27次,其中特大潮災(zāi)9次,而廣東省其他岸段的特大潮災(zāi)發(fā)生次數(shù)之和也僅有6次[1]。從歷史最大增水上看,南渡站在8007號(hào)強(qiáng)臺(tái)風(fēng)中測(cè)得585 cm的風(fēng)暴增水,是我國(guó)有驗(yàn)潮記錄以來的最大值,居世界第三;湛江港的歷史最大增水也達(dá)到497 cm,居廣東省第二[2]。從災(zāi)害程度上看,20世紀(jì)70年代以來廣東省死亡人數(shù)最多(8007號(hào)強(qiáng)臺(tái)風(fēng),共造成455人死亡)和經(jīng)濟(jì)損失最嚴(yán)重(9615號(hào)強(qiáng)臺(tái)風(fēng),直接經(jīng)濟(jì)損失201.25億元)的兩次風(fēng)暴潮災(zāi)害均發(fā)生在湛江地區(qū)。

    影響雷州半島的臺(tái)風(fēng)大體可分為3種類型。一是在湛江以南登陸的西北方向行進(jìn)的臺(tái)風(fēng),雷州半島處于這類臺(tái)風(fēng)行進(jìn)路徑的右側(cè)風(fēng)場(chǎng),常遭受強(qiáng)烈增水[3-4];二是在湛江或湛江以東登陸并北上的臺(tái)風(fēng),這類臺(tái)風(fēng)引起的增水通常要小于第一類路徑[5-6],但由于雷州半島部分地區(qū)處于臺(tái)風(fēng)行進(jìn)路徑的左側(cè),可能發(fā)生一定程度的減水;三是在海南島以南或繞過海南島在廣西沿海北上登陸的臺(tái)風(fēng),這類臺(tái)風(fēng)的發(fā)生頻次較少,對(duì)湛江的影響程度弱于第一類和第二類臺(tái)風(fēng)路徑。1522號(hào)臺(tái)風(fēng)“彩虹”屬于第二類臺(tái)風(fēng)路徑?!安屎纭痹谡拷衅骂^區(qū)沿海登陸,致使廣東、廣西和海南三地因?yàn)?zāi)直接經(jīng)濟(jì)損失27.02億元、死亡7人、受災(zāi)369.11萬人,是造成2015年中國(guó)經(jīng)濟(jì)損失最嚴(yán)重的臺(tái)風(fēng)。

    本文以雷州半島東部的實(shí)測(cè)氣象和潮位數(shù)據(jù)為驗(yàn)證資料,建立了大-中-小三層嵌套、局部高分辨率的二維風(fēng)暴潮與天文潮耦合數(shù)值模型。大模型覆蓋了整個(gè)南中國(guó)海和東中國(guó)海,與傳統(tǒng)的中尺度模擬區(qū)域相比,這種利用大模型提供邊界條件、局部模型精細(xì)計(jì)算的方法有利于提高局部模型邊界上的風(fēng)暴潮位精度,并加快計(jì)算效率。本文應(yīng)用該模型計(jì)算并分析了在1522號(hào)強(qiáng)臺(tái)風(fēng)“彩虹”作用下,雷州半島東部海域風(fēng)暴潮的時(shí)空分布規(guī)律,并研究了當(dāng)?shù)氐娘L(fēng)暴潮特征及產(chǎn)生原因。

    2 臺(tái)風(fēng)特征分析

    1522號(hào)臺(tái)風(fēng)“彩虹”生成于菲律賓以東洋面,于2015年10月2日02時(shí)(北京時(shí),下同)由熱帶低壓加強(qiáng)為熱帶風(fēng)暴,10時(shí)左右進(jìn)入南海并不斷向西北方向行進(jìn)。3日下午14時(shí)增強(qiáng)為臺(tái)風(fēng),晚上23時(shí)增強(qiáng)為強(qiáng)臺(tái)風(fēng),并于4日14時(shí)10分以強(qiáng)臺(tái)風(fēng)級(jí)別在廣東湛江市坡頭區(qū)登陸,登陸時(shí)近中心最大風(fēng)力達(dá)到15級(jí)(50 m/s),中心最低氣壓為935 hPa,是廣東省有氣象記錄以來10月份登陸的最強(qiáng)臺(tái)風(fēng)。“彩虹”登陸后繼續(xù)向西北方向移動(dòng),受下墊面摩擦的影響,強(qiáng)度不斷減弱,于5日11時(shí)在廣西省境內(nèi)減弱為熱帶低壓。圖1為臺(tái)風(fēng)“彩虹”的路徑示意圖。此次臺(tái)風(fēng)具有移動(dòng)速度快、登陸強(qiáng)度強(qiáng)等特點(diǎn),臺(tái)風(fēng)中心在大洋中的最大移動(dòng)速度約為25 m/s,并在登陸后造成兩廣多地發(fā)生暴雨[7]和龍卷風(fēng)[8]災(zāi)害。

    3 風(fēng)暴潮模型的建立與驗(yàn)證

    3.1 臺(tái)風(fēng)模型

    臺(tái)風(fēng)氣壓場(chǎng)由Holland模型給出:

    式中:p(r)是距離臺(tái)風(fēng)中心r處的氣壓;pc是臺(tái)風(fēng)中心最低氣壓;?p是外圍氣壓與中心最低氣壓的差值。該模型包含兩個(gè)可調(diào)參數(shù)Rmax和B:Rmax為最大風(fēng)速半徑,指臺(tái)風(fēng)水平結(jié)構(gòu)內(nèi)部風(fēng)速最大處到臺(tái)風(fēng)中心的距離;B為Holland引入的用來描述氣壓剖面形狀的參數(shù),根據(jù)臺(tái)風(fēng)實(shí)況選擇不同的參數(shù)B,可以增強(qiáng)模型的適用性。本文應(yīng)用林偉等[9]根據(jù)2001—2009年美國(guó)聯(lián)合臺(tái)風(fēng)警報(bào)中心整編的6 426條熱帶氣旋參數(shù)資料建立的最大風(fēng)速半徑與中心氣壓差的回歸方程作為經(jīng)驗(yàn)公式計(jì)算Rmax。

    圖1 1522號(hào)臺(tái)風(fēng)“彩虹”路徑圖(世界時(shí))

    B系數(shù)選用Vickery等根據(jù)大西洋海域的數(shù)據(jù)提出的計(jì)算公式給出:

    式中:φ為臺(tái)風(fēng)中心所在緯度。

    風(fēng)暴潮模型中的風(fēng)場(chǎng)采用背景風(fēng)場(chǎng)與臺(tái)風(fēng)經(jīng)驗(yàn)?zāi)P偷暮铣伞T谂_(tái)風(fēng)中心附近約幾百公里的范圍以內(nèi),經(jīng)驗(yàn)?zāi)P涂梢暂^好地反映臺(tái)風(fēng)大風(fēng)區(qū)的風(fēng)場(chǎng)特征。模型風(fēng)場(chǎng)為梯度風(fēng)風(fēng)場(chǎng)Vg與移行風(fēng)場(chǎng)Vt的疊加。

    式中:c1和c2是訂正系數(shù),c1=1.0,c2=0.8;θ是計(jì)算點(diǎn)與臺(tái)風(fēng)中心的連線與正東方向的夾角;β為流入角修正。梯度風(fēng)風(fēng)場(chǎng)Vg可表示為:

    式中:ρa(bǔ)為空氣密度,取1.15 kg/m3;f=2ωsinφ為科氏力系數(shù),ω=7.272×10-5rad/s為地球自轉(zhuǎn)角速度。移行風(fēng)場(chǎng)Vt采用宮崎正衛(wèi)公式,其分布形式為:

    圖2 研究區(qū)域地形及測(cè)站位置

    式中:Vx和Vy分別是臺(tái)風(fēng)中心移動(dòng)速度的正東和正北分量。而在臺(tái)風(fēng)外圍,由于同時(shí)受到其他天氣系統(tǒng)的影響,風(fēng)場(chǎng)與經(jīng)驗(yàn)?zāi)P筒顒e較大,需考慮疊加背景風(fēng)場(chǎng)。本文采用日本氣象廳提供的第二代全球氣候再分析數(shù)據(jù)庫(kù)JRA-55給出的風(fēng)速場(chǎng)作為背景風(fēng)場(chǎng)。背景風(fēng)場(chǎng)與臺(tái)風(fēng)經(jīng)驗(yàn)?zāi)P偷暮铣煞椒椋?/p>

    式中:VQ為背景風(fēng)場(chǎng);經(jīng)驗(yàn)風(fēng)場(chǎng)VM由式給出;e為權(quán)重系數(shù),e=c4/(1+c4),c=r/(10×Rmax)。

    3.2 風(fēng)暴潮模型

    本文的風(fēng)暴潮模型采用大、中、小三層嵌套的方法。大模型范圍包括東中國(guó)海和南中國(guó)海,模型采用球面坐標(biāo)系下方程,網(wǎng)格尺度為2'×2',離散采用DSI法,邊界采用8個(gè)主要調(diào)和常數(shù)(M2、S2、K1、O1、N2、K2、P1、Q1)預(yù)報(bào)潮位[10]。中模型范圍包括北部灣和粵西海岸,北邊界取在黃茅島(E113°02'03″,N21°56'52″)附近,南邊界取在越南會(huì)安(E108°24'08″,N15°52'07″)附近。模型采用無結(jié)構(gòu)三角形網(wǎng)格,網(wǎng)格尺寸在0.8~3 km之間,利用大模型提供邊界條件。小模型計(jì)算范圍包括湛江灣和雷州灣,采用無結(jié)構(gòu)三角形網(wǎng)格,開邊界處的網(wǎng)格大小約為1 km,雷州灣的岸邊界、島嶼附近的網(wǎng)格分辨率約為400 m,湛江灣及以北區(qū)域的分辨率約為150 m,海域地形如圖2所示,網(wǎng)格剖分見圖3。小范圍模型的潮位邊界條件由中模型給出。大、中、小模型的計(jì)算均同時(shí)考慮風(fēng)場(chǎng)、氣壓場(chǎng)作用。直角坐標(biāo)系下的二維潮流運(yùn)動(dòng)控制方程:

    式中:ζ水面高度;H=ζ+h為總水深;U、V為x、y方向的垂線平均流速;f為科氏系數(shù);pa為大氣壓強(qiáng);ρ為海水密度;g為重力加速度;εx、εy為x、y方向紊動(dòng)粘性系數(shù),按Smagorinsky公式計(jì)算;τsx、τsy為海面風(fēng)應(yīng)力在x、y方向的分量,表示為:

    式中:ρa(bǔ)為空氣密度為海上10 m處的風(fēng)速大小;wx、wy為x、y方向的風(fēng)速分量;cd為風(fēng)拖曳力系數(shù)。

    圖3 中尺度模型和局地模型網(wǎng)格剖分示意圖

    式中:τbx、τby為底摩擦應(yīng)力分量,表達(dá)式為:

    3.3 模型驗(yàn)證

    圖4 風(fēng)場(chǎng)模型驗(yàn)證結(jié)果

    水值與實(shí)測(cè)相比誤差較小,模擬得到的增水發(fā)生時(shí)刻也對(duì)應(yīng)較好。表明數(shù)值模型對(duì)真實(shí)風(fēng)暴潮過程的模擬具有較高精度。

    4 風(fēng)暴潮特征分析

    4.1 特征站點(diǎn)的風(fēng)暴潮位和增水分析

    在此次強(qiáng)臺(tái)風(fēng)過程中,湛江站出現(xiàn)較大的增水,而南渡則出現(xiàn)明顯的減水現(xiàn)象。湛江站的最高潮位為2.72 m(珠江基面),出現(xiàn)在4日03時(shí),與潮汐表顯示的該站當(dāng)日最高天文潮位的發(fā)生時(shí)刻一致。湛江站的最大增水為2.12 m,發(fā)生在登陸時(shí)刻4日14時(shí)。盡管增水值較大,但由于臺(tái)風(fēng)登陸時(shí)刻恰為天文潮平潮期,較大的增水并未引發(fā)超高水位。從增水過程上看,湛江站為標(biāo)準(zhǔn)型增水,曲線的先兆波、主振和余振3個(gè)階段十分明顯,主振有明顯的大突起峰值,發(fā)生在臺(tái)風(fēng)登陸前后1 h內(nèi)。

    南渡站的最高潮位為2.66 m(珠江基面),也出現(xiàn)在4日03時(shí)。南渡站的增水過程與湛江站不同,為波動(dòng)型增水,最大增水值較小,僅有1.48 m,發(fā)生在臺(tái)風(fēng)登陸前的4日12時(shí);而在4日15時(shí),即臺(tái)風(fēng)登陸后一個(gè)小時(shí),南渡站出現(xiàn)了2.11 m左右的減水。

    4.2 風(fēng)暴增水空間分布特征

    圖6和圖7分別為1522臺(tái)風(fēng)在雷州半島東部海域引起的最大風(fēng)暴增水與風(fēng)暴減水的分布圖。對(duì)比兩幅圖可知,受臺(tái)風(fēng)移動(dòng)方向和海灣地形的影響,1522號(hào)臺(tái)風(fēng)期間,雷州半島東部海域的風(fēng)暴潮增水分布情況大致以臺(tái)風(fēng)路徑為界,臺(tái)風(fēng)前進(jìn)路徑的右側(cè)以增水為主,左側(cè)則以減水現(xiàn)象較為顯著。這與臺(tái)風(fēng)登陸前后的風(fēng)向關(guān)系較為密切。當(dāng)臺(tái)風(fēng)中心登陸湛江灣時(shí),位于臺(tái)風(fēng)路徑右側(cè)的南三島、南三水道、鑒江等地受向岸的東南風(fēng)控制,海水易堆積于海灣和河道內(nèi);而位于臺(tái)風(fēng)路徑左側(cè)的雷州灣、南渡等地則受強(qiáng)烈的西風(fēng)作用,離岸的風(fēng)拖曳力使岸邊水位急劇下降;而在臺(tái)風(fēng)中心經(jīng)過的湛江灣內(nèi),受到登陸前向岸風(fēng)的吹刮以及登陸時(shí)氣壓降低的雙重作用,也發(fā)生了較大增水。

    圖5 計(jì)算風(fēng)暴潮過程與實(shí)測(cè)對(duì)比圖

    雷州半島東部海域的最大風(fēng)暴增水呈現(xiàn)由南向北逐漸增加的趨勢(shì)。由圖6可知,南部雷州灣內(nèi)的最大增水僅為1~1.6 m;湛江灣內(nèi)的增水在2~2.2 m之間,灣頂略有增加;全海域的最大增水發(fā)生在北部的南三水道和鑒江,鑒江內(nèi)有2.8 m以上的增水,而南三水道內(nèi)則有局部超過3.2 m的增水。

    圖6 1522號(hào)臺(tái)風(fēng)增水分布圖

    臺(tái)風(fēng)風(fēng)場(chǎng)與當(dāng)?shù)氐匦问怯绊戯L(fēng)暴增水分布的兩個(gè)重要因素。臺(tái)風(fēng)風(fēng)場(chǎng)很大程度上依賴于最大風(fēng)速半徑Rmax的值,該參數(shù)能夠在一定程度上反映了風(fēng)速高值區(qū)的影響范圍。由于臺(tái)風(fēng)氣旋結(jié)構(gòu)具有顯著的差異性,Rmax的值與中心氣壓差、緯度等許多因素有關(guān),需要根據(jù)臺(tái)風(fēng)氣象資料進(jìn)行擬合才可給定。在臺(tái)風(fēng)“彩虹”中,由Holland模型擬合湛江、陽(yáng)江、電白等多個(gè)氣象站資料得到的臺(tái)風(fēng)登陸時(shí)的最大風(fēng)速半徑約為20 km,從圖8中可以看出:北部的南三水道和鑒江恰位于臺(tái)風(fēng)路徑右側(cè)最大風(fēng)速半徑處,向岸風(fēng)作用最為強(qiáng)烈。且南三水道是向內(nèi)束窄的地形,向岸風(fēng)的作用使能量在水道內(nèi)集中、水體發(fā)生堆積并雍高。可見,南三水道和鑒江內(nèi)的劇烈增水反映了風(fēng)場(chǎng)和地形的雙重影響。

    圖7 1522號(hào)臺(tái)風(fēng)減水分布圖

    圖8 臺(tái)風(fēng)登陸前(4日11時(shí))和登陸時(shí)(4日14時(shí))的風(fēng)場(chǎng)示意圖

    “彩虹”引起的風(fēng)暴減水現(xiàn)象則集中發(fā)生在雷州灣西岸和北岸。雷州灣西岸南渡附近的減水值較大,達(dá)到2.1 m以上,北岸的減水稍小,普遍在1 m以內(nèi)。造成這些岸段減水的原因與強(qiáng)臺(tái)風(fēng)登陸前和登陸時(shí)的離岸風(fēng)有關(guān)。圖8顯示,臺(tái)風(fēng)登陸前,偏北風(fēng)的作用造成雷州灣北側(cè)岸段減水;隨著臺(tái)風(fēng)中心靠近,當(dāng)?shù)仫L(fēng)向轉(zhuǎn)而向西,造成西岸的南渡附近發(fā)生減水,且由于靠近臺(tái)風(fēng)中心,該海域的減水現(xiàn)象更為嚴(yán)重。陳奕德等[11]曾統(tǒng)計(jì)了湛江港在1953—1982年期間117場(chǎng)風(fēng)暴潮增水資料,其中減水超過0.7 m的共12場(chǎng),占資料總數(shù)的10.3%;尤其在6706、7220等在湛江以東登陸并北上的臺(tái)風(fēng)中,湛江港出現(xiàn)了大于1 m的風(fēng)暴減水[12]。風(fēng)暴減水會(huì)使大片海灘露出,嚴(yán)重影響艦船和大型油輪、貨輪的航行和錨泊[13]。需要指出的是,在統(tǒng)計(jì)熱帶氣旋在雷州半島引起的風(fēng)暴潮時(shí),尤其是對(duì)于臺(tái)風(fēng)中心在湛江附近及其以東登陸的情況,應(yīng)當(dāng)注意對(duì)風(fēng)暴減水值的統(tǒng)計(jì)。

    5 結(jié)論

    本文建立了三重嵌套的高分辨率二維風(fēng)暴潮與天文潮耦合數(shù)值模型,對(duì)1522號(hào)臺(tái)風(fēng)“彩虹”在雷州半島東部海域引起的風(fēng)暴增水特征進(jìn)行分析,得到以下結(jié)論:

    (1)臺(tái)風(fēng)中心正面登陸于湛江灣,湛江港出現(xiàn)較大增水,最大增水值為2.12 m,發(fā)生在臺(tái)風(fēng)登陸時(shí)刻;南渡則出現(xiàn)了較大減水,最大減水值為-2.11 m,發(fā)生在臺(tái)風(fēng)登陸后一個(gè)小時(shí)。由于臺(tái)風(fēng)登陸時(shí)刻恰為天文潮平潮期,較大的增水未引發(fā)超高水位;

    (2)臺(tái)風(fēng)期間,雷州半島東部海域的風(fēng)暴潮增水大致以臺(tái)風(fēng)路徑為界:位于臺(tái)風(fēng)前進(jìn)路徑右側(cè)的南三水道、鑒江等地,受向岸的東南風(fēng)控制,以增水為主;而位于臺(tái)風(fēng)前進(jìn)路徑左側(cè)的雷州灣內(nèi),受離岸西風(fēng)的控制,出現(xiàn)了顯著的減水現(xiàn)象;

    (3)研究海域的最大風(fēng)暴增水呈現(xiàn)由南向北逐漸增加的趨勢(shì),全海域的最大增水發(fā)生在南三水道內(nèi)。臺(tái)風(fēng)風(fēng)場(chǎng)與海灣地形是影響風(fēng)暴潮增水的兩個(gè)主要因素。對(duì)于向內(nèi)束窄的河道地形,向岸風(fēng)的吹刮使水體易于流入水道并雍高,造成嚴(yán)重的增水;

    (4)研究海域的風(fēng)暴減水主要發(fā)生在雷州灣的西岸和北岸,西岸的南渡河附近尤為嚴(yán)重。今后在統(tǒng)計(jì)熱帶氣旋經(jīng)過雷州半島引起的風(fēng)暴潮時(shí),尤其是對(duì)于臺(tái)風(fēng)中心恰好在湛江灣或其以東登陸的情況,應(yīng)當(dāng)注意對(duì)風(fēng)暴減水值的統(tǒng)計(jì);

    雷州半島是我國(guó)受臺(tái)風(fēng)影響最頻繁的地區(qū)之一,受北半球右旋風(fēng)場(chǎng)作用和南海開闊陸架地形的影響,臺(tái)風(fēng)路徑東側(cè)的海岸地區(qū)在遭受風(fēng)暴潮災(zāi)害的同時(shí),也易遭受嚴(yán)重的臺(tái)風(fēng)浪災(zāi)害,危及沿岸堤壩安全。今后需要重視對(duì)臺(tái)風(fēng)浪的模擬預(yù)報(bào),在本模型的基礎(chǔ)上進(jìn)一步發(fā)展風(fēng)-浪-流全耦合模型,提高湛江灣海洋災(zāi)害安全預(yù)警的研究及應(yīng)用水平。

    [1]張俊香,黃崇福,劉旭攏.廣東沿海臺(tái)風(fēng)暴潮災(zāi)害的地理分布特征和風(fēng)險(xiǎn)評(píng)估(1949—2005)[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào),2008,16(3):393-402.

    [2]于福江,董劍希,葉琳,等.中國(guó)風(fēng)暴潮災(zāi)害史料集(1949—2009)[M].北京:海洋出版社,2015.

    [3]李文歡,石?,?黃厚衡.0312號(hào)臺(tái)風(fēng)(科羅旺)風(fēng)暴潮預(yù)報(bào)過程分析[J].海洋預(yù)報(bào),2004,21(2):64-69.

    [4]張敏,陳鈺祥,趙雪,等.臺(tái)風(fēng)移動(dòng)方向和速度對(duì)湛江市沿海風(fēng)暴潮影響的數(shù)值分析[J].海洋預(yù)報(bào),2015,32(5):45-52.

    [5]王欣睿,孫波濤,陳強(qiáng),等.0606號(hào)臺(tái)風(fēng)“派比安”風(fēng)暴潮特征分析與總結(jié)[J].海洋預(yù)報(bào),2008,25(2):99-105.

    [6]張娟,馮偉忠,李廣敏,等.1003號(hào)臺(tái)風(fēng)“燦都”風(fēng)暴潮特征分析與模擬[J].海洋通報(bào),2011,30(4):367-370.

    [7]周渭,蔣平凡,黎馨.2015臺(tái)風(fēng)“彩虹”強(qiáng)降雨過程分析[J].氣象研究與應(yīng)用,2016,37(2):16-19.

    [8]李麗云,張錦華,諶志剛,等.1522號(hào)強(qiáng)臺(tái)風(fēng)“彩虹”環(huán)流中的龍卷風(fēng)過程分析[J].廣東氣象,2015,37(6):1-5.

    [9]林偉,方偉華.西北太平洋臺(tái)風(fēng)風(fēng)場(chǎng)模型中Holland B系數(shù)區(qū)域特征研究[J].熱帶地理,2013,33(2):124-132.

    [10]章衛(wèi)勝,陳晨,李鑫,等.海岸工程對(duì)渤海灣風(fēng)暴潮高潮位影響分析[J].海洋工程,2012,30(2):72-78.

    [11]陳奕德,董兆俊,蔣國(guó)榮,等.湛江港風(fēng)暴增水特征分析[J].海洋預(yù)報(bào),2002,19(3):44-52.

    [12]石靜慧,楊支中,沙文鈺.影響湛江的熱帶氣旋路徑、登陸地點(diǎn)和增水關(guān)系探討[J].海洋預(yù)報(bào),2004,21(1):43-51.

    [13]董勝,于亞群,余海靜.海岸帶風(fēng)暴潮減水的統(tǒng)計(jì)分析[J].自然災(zāi)害學(xué)報(bào),2004,13(4):70-74.

    Study on strom surge of the Typhoon“Mujigae”(1522)in the East Leizhou Peninsula

    XIONG Meng-jie1,2,ZHANG Wei-sheng1,2,ZHANG Jin-shan1,2,YIN Cheng-tuan1,2
    (1.Nanjing Hydraulic Research Institute,Nanjing 210024 China;2.State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering,Nanjing Hydraulic Research Institute,Nanjing 210024 China;3.Department of the PLA navy military facilities construction,Southern Theater,Zhanjiang China)

    A two-dimensional numerical model coupled with storm surge and astronomical tide,with a high resolution and three domain nested grid,is established to study the temporal and spatial distribution characteristic of the strom surge in typhoon Mujigae(1522).The simulation results show that the strom surge pattern in the east of Leizhou Peninsula is divided by the typhoon track.The right-side regions of the typhoon track,such as the Nansan Channel and the Jian River,are characterized by intense surge.The surge depends on the special topography of westward-curved bay shape and the sustained onshore wind before landing.While,the left-side regions of the typhoon track,such as the west coast of the Leizhou Bay,where suffered offshore west wind around the landfall of typhoon,are characterized by significant set-down.The maximum set-down reaches-2.11 m at Nandu.

    typhoon;Leizhou Peninsula;numerical simulation;strom surge;set-down

    P731.23

    A

    1003-0239(2017)06-0057-08

    10.11737/j.issn.1003-0239.2017.06.007

    2017-02-04;

    2017-03-29。

    國(guó)家重點(diǎn)研發(fā)計(jì)劃(2016FYC1402000)。

    熊夢(mèng)婕(1990-),女,助理工程師,碩士,從事河口海岸動(dòng)力學(xué)及海岸災(zāi)害數(shù)值模擬方面的研究。E-mail:mjxiong@nhri.cn

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點(diǎn)
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    国产av一区二区精品久久| 91大片在线观看| av电影中文网址| 国产无遮挡羞羞视频在线观看| 国产真人三级小视频在线观看| 视频区欧美日本亚洲| 波多野结衣一区麻豆| 精品久久久久久成人av| 不卡av一区二区三区| 99久久人妻综合| 一二三四社区在线视频社区8| 中国美女看黄片| 一级片'在线观看视频| 国产亚洲精品一区二区www| 自线自在国产av| 欧美午夜高清在线| 午夜福利在线观看吧| 自拍欧美九色日韩亚洲蝌蚪91| 精品国产国语对白av| 999久久久国产精品视频| 亚洲精品一二三| 亚洲五月天丁香| 精品日产1卡2卡| 亚洲自偷自拍图片 自拍| 少妇粗大呻吟视频| 久久精品人人爽人人爽视色| 国产成人av教育| 窝窝影院91人妻| 侵犯人妻中文字幕一二三四区| 国产又色又爽无遮挡免费看| а√天堂www在线а√下载| 十分钟在线观看高清视频www| 色播在线永久视频| 亚洲国产精品sss在线观看 | 亚洲免费av在线视频| 亚洲精品国产精品久久久不卡| 在线观看午夜福利视频| 免费久久久久久久精品成人欧美视频| 欧美精品亚洲一区二区| av在线播放免费不卡| 亚洲欧美一区二区三区黑人| 法律面前人人平等表现在哪些方面| 又黄又爽又免费观看的视频| 亚洲精品在线观看二区| 大陆偷拍与自拍| 国产国语露脸激情在线看| xxxhd国产人妻xxx| 91字幕亚洲| 亚洲精品国产精品久久久不卡| 欧美日韩一级在线毛片| 啦啦啦在线免费观看视频4| 亚洲精品中文字幕在线视频| videosex国产| 十分钟在线观看高清视频www| 国产成+人综合+亚洲专区| 成人国产一区最新在线观看| 淫秽高清视频在线观看| 99久久精品国产亚洲精品| 看黄色毛片网站| 亚洲av成人不卡在线观看播放网| 大陆偷拍与自拍| 国产区一区二久久| 欧美日韩av久久| 午夜a级毛片| 久久亚洲真实| 9色porny在线观看| 天天躁夜夜躁狠狠躁躁| 丝袜美腿诱惑在线| 多毛熟女@视频| 国产精品免费视频内射| 日本精品一区二区三区蜜桃| 午夜精品久久久久久毛片777| av国产精品久久久久影院| 女人被躁到高潮嗷嗷叫费观| 国产一区二区激情短视频| 精品熟女少妇八av免费久了| 亚洲国产欧美网| 亚洲伊人色综图| 天天添夜夜摸| 久久久国产一区二区| 国产精品国产高清国产av| 美女午夜性视频免费| 狠狠狠狠99中文字幕| 国产亚洲欧美精品永久| 中文亚洲av片在线观看爽| 免费av毛片视频| 国产有黄有色有爽视频| 国产精品久久视频播放| 首页视频小说图片口味搜索| 亚洲精品国产一区二区精华液| xxxhd国产人妻xxx| 国产精品乱码一区二三区的特点 | 久久国产精品影院| 亚洲成人国产一区在线观看| 一进一出抽搐动态| 又黄又粗又硬又大视频| 日韩精品免费视频一区二区三区| 国产精品偷伦视频观看了| 我的亚洲天堂| 好男人电影高清在线观看| 国产免费现黄频在线看| 中出人妻视频一区二区| 在线观看www视频免费| 长腿黑丝高跟| 亚洲欧美激情综合另类| 亚洲国产精品合色在线| www.999成人在线观看| 一区在线观看完整版| 美女福利国产在线| 美女午夜性视频免费| 1024视频免费在线观看| 成人精品一区二区免费| 国产xxxxx性猛交| 久久久久久久久久久久大奶| 日日爽夜夜爽网站| 天天躁狠狠躁夜夜躁狠狠躁| 中亚洲国语对白在线视频| 国产精品一区二区在线不卡| 久久欧美精品欧美久久欧美| 国产精品久久久久久人妻精品电影| 免费av中文字幕在线| 亚洲成人免费电影在线观看| 亚洲性夜色夜夜综合| 国产成人欧美在线观看| 国产精品亚洲一级av第二区| 免费高清视频大片| 99久久综合精品五月天人人| 超碰97精品在线观看| av有码第一页| 成年女人毛片免费观看观看9| 岛国在线观看网站| 日本精品一区二区三区蜜桃| 无人区码免费观看不卡| 高清欧美精品videossex| 精品国产一区二区三区四区第35| 一区在线观看完整版| 久久久久久久久中文| 色在线成人网| a在线观看视频网站| 日韩大尺度精品在线看网址 | 十八禁网站免费在线| 一进一出抽搐动态| 国产精品一区二区精品视频观看| 久久久久久久午夜电影 | 国产精品久久久人人做人人爽| 久久性视频一级片| 亚洲午夜理论影院| 欧美午夜高清在线| 国产精品美女特级片免费视频播放器 | www.999成人在线观看| 国产成人av教育| 免费日韩欧美在线观看| 亚洲,欧美精品.| av福利片在线| 咕卡用的链子| 欧美黄色片欧美黄色片| 久99久视频精品免费| 久久精品亚洲精品国产色婷小说| 精品乱码久久久久久99久播| e午夜精品久久久久久久| 日韩精品中文字幕看吧| 亚洲第一青青草原| 性欧美人与动物交配| 日本成人三级电影网站| 欧美性感艳星| 日韩大尺度精品在线看网址| 国产爱豆传媒在线观看| 丰满人妻一区二区三区视频av| 久久午夜福利片| 黄色丝袜av网址大全| 在线免费观看的www视频| 亚洲国产欧美人成| 久久久久国内视频| 欧美黄色淫秽网站| 国产大屁股一区二区在线视频| 熟妇人妻久久中文字幕3abv| 身体一侧抽搐| 男女下面进入的视频免费午夜| 99国产综合亚洲精品| 欧美日韩福利视频一区二区| 国产又黄又爽又无遮挡在线| 免费无遮挡裸体视频| 午夜福利18| 露出奶头的视频| 国产亚洲av嫩草精品影院| 国产成人啪精品午夜网站| 国产探花极品一区二区| 午夜影院日韩av| 最近视频中文字幕2019在线8| 蜜桃久久精品国产亚洲av| 国产一区二区三区视频了| 偷拍熟女少妇极品色| 国产 一区 欧美 日韩| av天堂在线播放| 一区二区三区免费毛片| 啦啦啦观看免费观看视频高清| 国产精品久久久久久亚洲av鲁大| 欧美成狂野欧美在线观看| 国产成+人综合+亚洲专区| 久久中文看片网| 久久99热6这里只有精品| 欧美区成人在线视频| 国产精品自产拍在线观看55亚洲| 精品一区二区三区av网在线观看| 日本与韩国留学比较| www.www免费av| 免费av观看视频| 高清日韩中文字幕在线| 亚洲av电影不卡..在线观看| 免费看a级黄色片| 久久天躁狠狠躁夜夜2o2o| 可以在线观看毛片的网站| 一级a爱片免费观看的视频| 男人和女人高潮做爰伦理| 99国产精品一区二区三区| 国产蜜桃级精品一区二区三区| 在线免费观看不下载黄p国产 | avwww免费| 日韩有码中文字幕| 一区二区三区激情视频| 欧美zozozo另类| 可以在线观看毛片的网站| 午夜a级毛片| 无人区码免费观看不卡| av在线蜜桃| 99国产精品一区二区蜜桃av| 国产伦在线观看视频一区| 亚洲av美国av| 国内久久婷婷六月综合欲色啪| 亚洲av成人不卡在线观看播放网| 99在线视频只有这里精品首页| 757午夜福利合集在线观看| 亚洲国产欧洲综合997久久,| 高清在线国产一区| 波多野结衣高清无吗| 日韩欧美免费精品| 亚洲欧美清纯卡通| 亚洲无线在线观看| 国产成人欧美在线观看| 激情在线观看视频在线高清| 草草在线视频免费看| 性插视频无遮挡在线免费观看| 91麻豆av在线| 久久香蕉精品热| 色哟哟·www| 男人狂女人下面高潮的视频| 波野结衣二区三区在线| 亚洲美女视频黄频| 国产欧美日韩一区二区精品| 国产欧美日韩精品一区二区| av在线蜜桃| 免费观看的影片在线观看| 一个人看的www免费观看视频| 国产午夜福利久久久久久| 亚洲综合色惰| 欧美xxxx黑人xx丫x性爽| bbb黄色大片| 久久精品国产99精品国产亚洲性色| 啦啦啦韩国在线观看视频| 中文亚洲av片在线观看爽| 免费av毛片视频| 欧美激情在线99| 99久久精品国产亚洲精品| 国产免费男女视频| 久久精品久久久久久噜噜老黄 | 中文字幕高清在线视频| 亚洲真实伦在线观看| 亚洲av第一区精品v没综合| 性欧美人与动物交配| 中文在线观看免费www的网站| 波多野结衣巨乳人妻| 三级国产精品欧美在线观看| 国产亚洲欧美在线一区二区| 好男人在线观看高清免费视频| 亚洲第一电影网av| 少妇高潮的动态图| 大型黄色视频在线免费观看| 亚洲综合色惰| 男女之事视频高清在线观看| av黄色大香蕉| 一a级毛片在线观看| 三级毛片av免费| 国产视频一区二区在线看| 国产伦人伦偷精品视频| 免费av观看视频| 亚洲18禁久久av| 亚洲男人的天堂狠狠| 精品人妻1区二区| 国产精华一区二区三区| 女同久久另类99精品国产91| 女生性感内裤真人,穿戴方法视频| 又爽又黄无遮挡网站| 一区二区三区高清视频在线| 亚洲最大成人手机在线| 女同久久另类99精品国产91| 亚洲五月天丁香| 国产成+人综合+亚洲专区| 色哟哟·www| 最近最新中文字幕大全电影3| 少妇被粗大猛烈的视频| 非洲黑人性xxxx精品又粗又长| 成年女人看的毛片在线观看| 高清毛片免费观看视频网站| 一级黄色大片毛片| 一区二区三区免费毛片| 少妇高潮的动态图| 99国产精品一区二区蜜桃av| 人妻制服诱惑在线中文字幕| 日本免费a在线| 美女被艹到高潮喷水动态| 久久久成人免费电影| 国产v大片淫在线免费观看| 男人舔奶头视频| 99在线视频只有这里精品首页| 在线观看免费视频日本深夜| 嫩草影院入口| 日韩精品中文字幕看吧| 女人十人毛片免费观看3o分钟| 国产熟女xx| 色哟哟·www| 国产成人a区在线观看| 午夜a级毛片| 亚洲成人精品中文字幕电影| av欧美777| 亚洲无线观看免费| 国产高清视频在线播放一区| 国产一区二区在线观看日韩| 亚洲av一区综合| 在线播放国产精品三级| www.色视频.com| 午夜福利成人在线免费观看| eeuss影院久久| 久久久国产成人精品二区| 国产三级中文精品| 国内精品久久久久久久电影| 人人妻,人人澡人人爽秒播| 麻豆成人午夜福利视频| 国产免费一级a男人的天堂| 日韩欧美精品免费久久 | 舔av片在线| 亚洲av中文字字幕乱码综合| 久久久久久国产a免费观看| 夜夜躁狠狠躁天天躁| 色吧在线观看| 久久久久久久亚洲中文字幕 | 99热精品在线国产| 亚洲熟妇熟女久久| 99热精品在线国产| 亚洲国产欧美人成| 人人妻人人看人人澡| 精品久久国产蜜桃| 中文字幕高清在线视频| 亚洲狠狠婷婷综合久久图片| 午夜日韩欧美国产| 亚洲最大成人av| av专区在线播放| 国产高清视频在线观看网站| 脱女人内裤的视频| 色综合站精品国产| 亚洲成av人片在线播放无| 成人特级av手机在线观看| 99久久精品一区二区三区| 免费高清视频大片| 免费一级毛片在线播放高清视频| 成年女人永久免费观看视频| 午夜免费男女啪啪视频观看 | 欧美黄色片欧美黄色片| 午夜精品一区二区三区免费看| 国产精品免费一区二区三区在线| 又黄又爽又刺激的免费视频.| 亚洲激情在线av| 久久精品国产亚洲av天美| 日本 av在线| 色av中文字幕| av欧美777| 亚洲美女搞黄在线观看 | 色在线成人网| 婷婷丁香在线五月| 国产主播在线观看一区二区| 级片在线观看| 亚洲,欧美精品.| 真实男女啪啪啪动态图| 亚洲自偷自拍三级| 真实男女啪啪啪动态图| АⅤ资源中文在线天堂| 午夜视频国产福利| 中文字幕熟女人妻在线| 欧美乱色亚洲激情| 免费人成视频x8x8入口观看| 久久国产精品影院| 亚洲av免费高清在线观看| 亚洲 欧美 日韩 在线 免费| 国产亚洲av嫩草精品影院| 成年女人看的毛片在线观看| 国产人妻一区二区三区在| 亚洲专区中文字幕在线| 国产亚洲欧美98| 深夜a级毛片| 成年女人毛片免费观看观看9| 18禁黄网站禁片免费观看直播| 色综合欧美亚洲国产小说| 精华霜和精华液先用哪个| 老司机午夜福利在线观看视频| www日本黄色视频网| 亚洲男人的天堂狠狠| 婷婷六月久久综合丁香| 又黄又爽又刺激的免费视频.| 亚洲18禁久久av| 免费人成在线观看视频色| 琪琪午夜伦伦电影理论片6080| 精品久久久久久久人妻蜜臀av| 2021天堂中文幕一二区在线观| 深夜a级毛片| 亚洲成人精品中文字幕电影| 久久久久久大精品| 亚洲人成网站高清观看| 五月伊人婷婷丁香| 亚洲人与动物交配视频| 看十八女毛片水多多多| 黄色视频,在线免费观看| 国产精品一区二区免费欧美| 少妇丰满av| 成年免费大片在线观看| 日韩免费av在线播放| 国产精品久久久久久人妻精品电影| 丁香六月欧美| 国产精品久久久久久精品电影| 国产精华一区二区三区| 搡老熟女国产l中国老女人| 五月玫瑰六月丁香| 99热这里只有精品一区| 国产野战对白在线观看| 青草久久国产| 亚洲精品一卡2卡三卡4卡5卡| 国产黄色小视频在线观看| 中文字幕av在线有码专区| 久久久久久久久大av| 极品教师在线视频| 日本五十路高清| 精品人妻熟女av久视频| 国内精品美女久久久久久| 嫩草影院入口| 欧美日韩瑟瑟在线播放| 淫秽高清视频在线观看| 亚洲欧美日韩高清在线视频| 精品一区二区免费观看| 亚洲avbb在线观看| 91av网一区二区| 老司机午夜十八禁免费视频| 成人毛片a级毛片在线播放| 欧美乱妇无乱码| 日本成人三级电影网站| 69av精品久久久久久| 特级一级黄色大片| 乱人视频在线观看| 国产在线男女| 亚洲av成人av| 久久久久久久亚洲中文字幕 | 少妇人妻精品综合一区二区 | 中文字幕熟女人妻在线| av女优亚洲男人天堂| 又紧又爽又黄一区二区| 岛国在线免费视频观看| 韩国av一区二区三区四区| 天堂√8在线中文| 在线播放国产精品三级| 国产不卡一卡二| 国产单亲对白刺激| 国产男靠女视频免费网站| 久99久视频精品免费| 亚洲国产高清在线一区二区三| 女人十人毛片免费观看3o分钟| 国产午夜精品论理片| a级毛片免费高清观看在线播放| 成人美女网站在线观看视频| 国产高清激情床上av| 中出人妻视频一区二区| 日本与韩国留学比较| 一级黄片播放器| 成人高潮视频无遮挡免费网站| 欧美bdsm另类| 一区二区三区免费毛片| 少妇的逼水好多| 亚洲自偷自拍三级| 在线看三级毛片| h日本视频在线播放| 亚洲,欧美精品.| 精品无人区乱码1区二区| 欧美最新免费一区二区三区 | 国产精品久久视频播放| 可以在线观看毛片的网站| 日韩中字成人| 午夜精品久久久久久毛片777| 亚洲va日本ⅴa欧美va伊人久久| 深夜精品福利| 尤物成人国产欧美一区二区三区| 黄色女人牲交| 午夜视频国产福利| 天天躁日日操中文字幕| 91久久精品电影网| 人妻久久中文字幕网| 亚洲国产欧洲综合997久久,| 国产精品99久久久久久久久| 十八禁人妻一区二区| 国产精品久久久久久亚洲av鲁大| 别揉我奶头 嗯啊视频| 国产三级黄色录像| 99久久九九国产精品国产免费| 欧美色视频一区免费| 国产午夜精品久久久久久一区二区三区 | 白带黄色成豆腐渣| 国产亚洲精品久久久com| 久久久久亚洲av毛片大全| 免费黄网站久久成人精品 | 久久久久久久久中文| 夜夜爽天天搞| 91字幕亚洲| 内地一区二区视频在线| 亚洲精品粉嫩美女一区| 在线十欧美十亚洲十日本专区| 日本 欧美在线| 亚洲专区中文字幕在线| 国产成人啪精品午夜网站| 麻豆国产97在线/欧美| 久久久久国产精品人妻aⅴ院| 1024手机看黄色片| 3wmmmm亚洲av在线观看| 可以在线观看毛片的网站| 久久6这里有精品| 日韩欧美 国产精品| 亚洲av电影不卡..在线观看| 色哟哟·www| 成年版毛片免费区| 51国产日韩欧美| 国产主播在线观看一区二区| 小蜜桃在线观看免费完整版高清| 中亚洲国语对白在线视频| 成人国产一区最新在线观看| 国内久久婷婷六月综合欲色啪| 琪琪午夜伦伦电影理论片6080| 精品福利观看| 三级毛片av免费| 久久99热这里只有精品18| 日韩欧美在线乱码| 18禁裸乳无遮挡免费网站照片| 最好的美女福利视频网| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 中出人妻视频一区二区| 久久人人爽人人爽人人片va | 亚洲最大成人av| 麻豆成人av在线观看| 亚洲人成伊人成综合网2020| 色吧在线观看| 久久精品国产亚洲av涩爱 | 91av网一区二区| 免费观看精品视频网站| 国产精品久久久久久人妻精品电影| 精品人妻偷拍中文字幕| 高清在线国产一区| 偷拍熟女少妇极品色| 国产极品精品免费视频能看的| 极品教师在线视频| 国产三级黄色录像| 在线a可以看的网站| 久久99热6这里只有精品| 美女高潮喷水抽搐中文字幕| 天天一区二区日本电影三级| 亚洲人成电影免费在线| 男女做爰动态图高潮gif福利片| 免费av观看视频| 嫁个100分男人电影在线观看| 午夜影院日韩av| 亚洲欧美日韩高清在线视频| 91久久精品国产一区二区成人| 国产三级中文精品| 真实男女啪啪啪动态图| 两个人视频免费观看高清| 日韩免费av在线播放| 精品久久久久久久久久久久久| 美女xxoo啪啪120秒动态图 | 成人鲁丝片一二三区免费| 亚洲av成人av| h日本视频在线播放| 99久久精品国产亚洲精品| 男女床上黄色一级片免费看| 欧美激情久久久久久爽电影| 在线免费观看不下载黄p国产 | 久久午夜亚洲精品久久| 最新中文字幕久久久久| 免费av观看视频| 亚洲av免费高清在线观看| 午夜a级毛片| 在线观看舔阴道视频| 好男人在线观看高清免费视频| 搡女人真爽免费视频火全软件 | 最新在线观看一区二区三区| 亚洲国产欧洲综合997久久,| 一个人看的www免费观看视频| 亚洲成a人片在线一区二区| 亚洲av免费在线观看| 老熟妇乱子伦视频在线观看| 少妇人妻精品综合一区二区 | 丰满人妻熟妇乱又伦精品不卡| 一级毛片久久久久久久久女| 成人毛片a级毛片在线播放| 又爽又黄a免费视频| 51午夜福利影视在线观看| 亚洲精品亚洲一区二区| 精品久久久久久久久久免费视频| 免费搜索国产男女视频| 国产三级在线视频| 美女被艹到高潮喷水动态| 高清毛片免费观看视频网站| 国产不卡一卡二| 91九色精品人成在线观看|