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

    珠海市沿海碼頭風(fēng)暴增水及最高潮位情景預(yù)測(cè)

    2022-12-27 11:32:42孫弦王靜夏冬聶會(huì)文高斯蘇燁康吳志棚
    廣東氣象 2022年6期
    關(guān)鍵詞:九洲風(fēng)暴潮山竹

    孫弦,王靜,夏冬,聶會(huì)文,高斯,蘇燁康,吳志棚

    (1.珠海市公共氣象服務(wù)中心,廣東 珠海 5 1 9 0 0 0;2.中山大學(xué)大氣科學(xué)學(xué)院,廣東 珠海 519082;3.南方海洋科學(xué)與工程廣東省實(shí)驗(yàn)室(珠海),廣東 珠海 519082)

    珠海市地處珠江口西側(cè),瀕臨廣闊的南海,幾乎每年都受臺(tái)風(fēng)風(fēng)暴潮增水影響[1]。九洲港作為全國(guó)沿海最大的港口客運(yùn)站之一,在2017年強(qiáng)臺(tái)風(fēng)“天鴿”和2018年強(qiáng)臺(tái)風(fēng)“山竹”期間,受風(fēng)暴潮和天文潮共同影響,受損十分嚴(yán)重。因此,對(duì)臺(tái)風(fēng)所帶來(lái)的風(fēng)暴增水及沿海碼頭可能出現(xiàn)的最高潮位進(jìn)行模擬和預(yù)測(cè),對(duì)防災(zāi)減災(zāi)具有重要意義。

    目前國(guó)內(nèi)外潮位計(jì)算方法主要為極值潮位統(tǒng)計(jì)法、風(fēng)暴潮增減水與天文潮組合法、風(fēng)暴潮增減水與天文潮聯(lián)合概率法以及風(fēng)暴增水模式法共4種[2-3],其中前3種方法需要大量長(zhǎng)期潮位觀測(cè)數(shù)據(jù),對(duì)歷史數(shù)據(jù)依賴度高,因此有很大局限。風(fēng)暴增水模式法從成因分析,建立風(fēng)暴增水模式,根據(jù)歷史臺(tái)風(fēng)資料計(jì)算增水值,并疊加當(dāng)時(shí)的天文潮,得到潮位極值。端義宏等[4]分別利用特定統(tǒng)計(jì)模式和高分辨率ECOM-Si的數(shù)值模式,計(jì)算可能最大臺(tái)風(fēng)對(duì)上海地區(qū)的最大可能增水影響,分析認(rèn)為高分辨率ECOM-Si的數(shù)值模式計(jì)算得到的增水結(jié)果更符合實(shí)際;周才揚(yáng)等[5]采用ERA-Interim數(shù)據(jù)集作為背景風(fēng)場(chǎng)資料建立了雙重嵌套的高分辨率風(fēng)暴潮與天文潮耦合數(shù)學(xué)模型,采用1810號(hào)臺(tái)風(fēng)“安比”實(shí)測(cè)風(fēng)暴增水?dāng)?shù)據(jù)對(duì)模型進(jìn)行了驗(yàn)證,結(jié)果表明風(fēng)暴潮模型能夠較好地模擬河口地區(qū)徑流、天文潮和風(fēng)暴潮的相互作用下的風(fēng)暴潮位變化。研究表明,根據(jù)臺(tái)風(fēng)資料進(jìn)行風(fēng)暴增水的數(shù)值模擬,方法符合實(shí)際,目前得到廣泛研究和越來(lái)越多的應(yīng)用[6-7]。但是,風(fēng)暴增水模式受局地影響較大,針對(duì)珠海沿海碼頭的風(fēng)暴潮模式目前還比較缺乏。

    另外,在利用風(fēng)暴潮模式進(jìn)行增水模擬時(shí),初始強(qiáng)迫場(chǎng)為根據(jù)歷史臺(tái)風(fēng)資料構(gòu)建的可能最大臺(tái)風(fēng),此時(shí)需要用極值模型來(lái)進(jìn)行模擬,目前尚無(wú)統(tǒng)一的極值模型。本研究將引入廣義Pareto分布(GPD)極值模型[8],計(jì)算不同重現(xiàn)期臺(tái)風(fēng)中心最低氣壓和最大風(fēng)速等極值,構(gòu)建出可能出現(xiàn)的最大臺(tái)風(fēng),將其作為初試強(qiáng)迫場(chǎng)進(jìn)行風(fēng)暴潮模擬,確定不同重現(xiàn)期最高增水情況,再疊加天文潮進(jìn)而得出預(yù)測(cè)的最高潮位,從而構(gòu)建出更加符合實(shí)際的珠海沿海碼頭百年一遇的臺(tái)風(fēng)風(fēng)暴潮情景。

    1 資料與方法

    1.1 潮位站資料

    根據(jù)《2017年中國(guó)海洋災(zāi)害公報(bào)》提供的資料,“天鴿”期間的珠海站最大風(fēng)暴增水為2.79 m[9],可作為風(fēng)暴潮模型驗(yàn)證資料;《2018年中國(guó)海洋災(zāi)害公報(bào)》中“山竹”期間的珠海站最大風(fēng)暴增水未給出[10],因此用周邊的三灶站、橫門(mén)站的最大風(fēng)暴增水作為風(fēng)暴潮模型驗(yàn)證資料,分別為3.39和2.89 m。

    1.2 登陸或嚴(yán)重影響珠海的臺(tái)風(fēng)資料

    為計(jì)算不同重現(xiàn)期的潮位極值,通過(guò)歷史臺(tái)風(fēng)的路徑、速度、中心氣壓以及風(fēng)速等資料,構(gòu)建可能最大臺(tái)風(fēng)進(jìn)行風(fēng)暴增水模擬。表1列出了近30年登陸或嚴(yán)重影響珠海的臺(tái)風(fēng)(數(shù)據(jù)來(lái)源于中央氣象臺(tái)臺(tái)風(fēng)網(wǎng)),作為風(fēng)暴潮模式的驅(qū)動(dòng)場(chǎng)。表1中臺(tái)風(fēng)篩選的依據(jù)是登陸或嚴(yán)重影響珠海市。嚴(yán)重影響是指在臺(tái)風(fēng)系統(tǒng)影響下,本地平均風(fēng)力≥8級(jí);或平均風(fēng)力6~7級(jí),且24 h雨量≥80 mm;或48 h雨量≥150 mm。例如1716臺(tái)風(fēng)“瑪娃”在汕尾登陸,2017年9月3—4日受其影響,珠海國(guó)家基本氣象站(區(qū)站號(hào)59488)錄得過(guò)程極大風(fēng)速13.6 m/s(6級(jí)),其中4日錄得最大日雨量205.5 mm(大暴雨),對(duì)珠海造成嚴(yán)重影響。

    表1 近30年登陸或嚴(yán)重影響珠海的臺(tái)風(fēng)資料

    1.3 GPD方法

    關(guān)于極值的理論概率分布,目前尚無(wú)公認(rèn)的統(tǒng)一模型,常用的極值模型有極值Ⅰ型(Gumbel分布)、皮爾遜Ⅲ型分布、Weibull分布、廣義極值分布(GEV),樣本抽取通常采用“年極值抽樣(AM)”,但有的年份影響本地的臺(tái)風(fēng)超過(guò)1個(gè),而有的年份臺(tái)風(fēng)影響又很小,因此之前提到的極值模型以及AM樣本抽取受到限制,容易造成數(shù)據(jù)浪費(fèi)。基于“超門(mén)限峰值POT”的GPD能充分利用有效數(shù)據(jù)信息,不受原始序列樣本量的影響,具有全部取值域的高精度穩(wěn)定擬合,具有高精度的實(shí)用性和穩(wěn)定性。GPD的分布函數(shù)為

    其中,α為尺度參數(shù),起放大或縮小分布面積的作用,描述極值分布的變率,尺度參數(shù)越大,表明極值波動(dòng)范圍越大;β為閾值;k為形狀參數(shù),描述極值的分布情況。GPD的參數(shù)估計(jì)方法有多種,如極大似然估計(jì)、概率權(quán)矩估計(jì)、L矩估計(jì)等,在風(fēng)壓計(jì)算方面,極大似然估計(jì)具有相合性和漸近最優(yōu)性,是參數(shù)估計(jì)最常用的方法[11]。

    1.4 風(fēng)暴潮和天文潮模式

    風(fēng)暴潮模式是基于日本氣象廳JMA研發(fā)的風(fēng)暴潮模式框架,進(jìn)行本地化和模式改進(jìn)建立的,該模式是基于淺水波方程的二維線性正壓模式,目的是預(yù)報(bào)由于風(fēng)攪拌和氣壓抽吸作用引起的水位增量,目前已在業(yè)務(wù)系統(tǒng)中廣泛使用,且模擬效果較好[12]。模式水平分辨率為3 km,覆蓋范圍100°E—130°E,0—30°N。

    本研究使用經(jīng)典調(diào)和分析例行程序包T_tide來(lái)進(jìn)行潮汐預(yù)測(cè)[13]。T_tide根據(jù)古典潮汐計(jì)算理論對(duì)潮汐進(jìn)行分析,針對(duì)大約1年或更短的潮汐資料進(jìn)行調(diào)和分析,并利用節(jié)點(diǎn)的修正解釋未定潮汐成分;T_predic.m利用T_tide.m分析的結(jié)果和計(jì)算出的調(diào)和常數(shù)對(duì)未來(lái)某時(shí)刻的潮汐進(jìn)行預(yù)報(bào)。

    2 風(fēng)暴潮增水?dāng)?shù)學(xué)模型驗(yàn)證計(jì)算

    2017年強(qiáng)臺(tái)風(fēng)“天鴿”和2018年強(qiáng)臺(tái)風(fēng)“山竹”給廣東帶來(lái)了嚴(yán)重的風(fēng)暴潮災(zāi)害,珠海更是受災(zāi)嚴(yán)重。圖1為臺(tái)風(fēng)“天鴿”和“山竹”的路徑示意圖,可以看出,這兩個(gè)臺(tái)風(fēng)均為自東南向西偏北方向移動(dòng),活動(dòng)范圍介于10°N—24°N之間。

    圖1 2017年強(qiáng)臺(tái)風(fēng)“天鴿”(a)和2018年強(qiáng)臺(tái)風(fēng)“山竹”(b)的移動(dòng)路徑

    圖2為風(fēng)暴潮模式模擬的“天鴿”、“山竹”影響期間珠海站及九洲港等站的臺(tái)風(fēng)增水過(guò)程曲線。

    圖2 各站點(diǎn)在強(qiáng)臺(tái)風(fēng)“天鴿”(a)和強(qiáng)臺(tái)風(fēng)“山竹”(b)期間的增水過(guò)程曲線

    從圖2a可以看出,珠海站從23日12:00(北京時(shí),下同)起潮位急增,達(dá)到2.7 m,并且一直維持到當(dāng)天22:00,期間出現(xiàn)了最高2.9 m的風(fēng)暴增水,與實(shí)測(cè)的2.79 m相差較小,模擬效果較好;九洲港風(fēng)暴增水最高出現(xiàn)在23日16:00,為3.7 m,配合當(dāng)天最高天文潮位1.1 m(基于1985國(guó)家高程基準(zhǔn)),九洲港在“天鴿”期間最高潮位可達(dá)為4.8 m。

    從圖2b“山竹”影響期間的橫門(mén)站、三灶站、珠海站和九洲港的臺(tái)風(fēng)增水過(guò)程曲線可以看出,三灶站的風(fēng)暴增水最高,為3.6 m,與觀測(cè)值的相對(duì)誤差為5.9%;橫門(mén)站的風(fēng)暴增水最高,為2.6 m,與觀測(cè)值的相對(duì)誤差為9.3%。珠海站最高增水出現(xiàn)在16日17:00,為2.7 m,九洲港在16:00—17:00也有明顯的風(fēng)暴增水高峰期,可達(dá)2.6m,并且在21:00之后還有一個(gè)峰值區(qū),最大增水接近2.9 m,配合當(dāng)天最高天文潮位0.7 m(基于1985國(guó)家高程基準(zhǔn)),九洲港在“山竹”期間最高潮位可達(dá)為3.6 m。

    3 最大可能風(fēng)暴增水

    在沿海重大工程項(xiàng)目建設(shè)和防災(zāi)規(guī)劃制定時(shí),應(yīng)充分考慮風(fēng)暴潮襲擊的頻繁程度和可能最大臺(tái)風(fēng)增水??赡茏畲笈_(tái)風(fēng)的主要參數(shù)包括:臺(tái)風(fēng)中心氣壓、臺(tái)風(fēng)中心最大風(fēng)速、臺(tái)風(fēng)移動(dòng)速度以及臺(tái)風(fēng)登陸時(shí)與海岸的交角?;诘顷懟驀?yán)重影響珠海的歷史臺(tái)風(fēng)資料(表1),構(gòu)建出對(duì)九洲港可能出現(xiàn)的最大臺(tái)風(fēng)、臺(tái)風(fēng)移速和路徑也依據(jù)統(tǒng)計(jì)方法計(jì)算后得到的[14],最終模擬出可能出現(xiàn)的最大臺(tái)風(fēng)增水。

    ΔP為臺(tái)風(fēng)外圍不受臺(tái)風(fēng)影響處的海平面最大氣壓P∞和臺(tái)風(fēng)中心最低氣壓P0兩者之差,其中P∞取我國(guó)沿海夏季海平面氣壓的平均值(1 008 hPa),對(duì)ΔP進(jìn)行基于閾值抽樣的GPD擬合[15],參數(shù)估計(jì)使用較好的極大似然估計(jì),得到50和100年一遇的臺(tái)風(fēng)中心最低氣壓分別為906和895 hPa。

    同樣對(duì)最大風(fēng)速進(jìn)行基于閾值抽樣的GPD擬合,參數(shù)估計(jì)使用較好的極大似然估計(jì),得到50和100年一遇的最大風(fēng)速分別為64和68 m/s。

    臺(tái)風(fēng)風(fēng)暴數(shù)值模式的研究成果表明,對(duì)于開(kāi)闊水域,登陸臺(tái)風(fēng)或緊靠岸邊又平行于海岸移動(dòng)的臺(tái)風(fēng),臺(tái)風(fēng)移動(dòng)速度越快,臺(tái)風(fēng)風(fēng)暴增水越大[12],但對(duì)于曲折海岸就復(fù)雜多了。在進(jìn)行可能影響最大臺(tái)風(fēng)計(jì)算時(shí),統(tǒng)計(jì)登陸或影響珠海的臺(tái)風(fēng)登陸時(shí)的移動(dòng)速度前10名,取其平均值,為28 km/h。

    確定臺(tái)風(fēng)移動(dòng)的方位和登陸的位置,是設(shè)計(jì)臺(tái)風(fēng)移動(dòng)路徑2個(gè)重要的元素。將登陸地點(diǎn)定為九洲港,移動(dòng)方位為270°到360°,間隔為10°。

    利用風(fēng)暴潮模式模擬重現(xiàn)期為100年的不同移動(dòng)方向的風(fēng)暴增水,圖3a為100年一遇、路徑為270°~360°的臺(tái)風(fēng)增水過(guò)程曲線,發(fā)現(xiàn)從270°到360°,移動(dòng)方向越接近正北,增水越大,100年一遇的最大風(fēng)暴增水為5.1 m。圖3b為重現(xiàn)期為100年的可能最大增水分布,從圖3b可以看出,九洲港北側(cè)的增水要比南側(cè)更高,北側(cè)增水可達(dá)5.2 m以上,南側(cè)增水在4.8到5.0 m之間。

    圖3 重現(xiàn)期為100年、路徑為270°~360°的臺(tái)風(fēng)增水過(guò)程曲線(a)和可能最大臺(tái)風(fēng)增水分布(b)

    4 天文潮調(diào)和分析

    采用九洲港潮位站2018年10月1日00:00到2019年9月30日23:00間隔1 h的潮位資料進(jìn)行調(diào)和分析,T_tide中設(shè)置開(kāi)始時(shí)間為2018年10月1日00:00,緯度為22.241°N,信噪比取2,計(jì)算得到該站點(diǎn)處的潮汐調(diào)和常數(shù),分析回歸估計(jì)值與實(shí)際觀測(cè)值之間的差。該潮位站的回報(bào)準(zhǔn)確率為96.9%,并通過(guò)了95%的置信區(qū)間檢驗(yàn),反演結(jié)果好(圖4)。

    圖4 真實(shí)和反演的潮汐對(duì)比

    在成功實(shí)現(xiàn)了對(duì)當(dāng)前時(shí)次的反演后,對(duì)未來(lái)100年的可能最高潮位進(jìn)行預(yù)測(cè)。從理論數(shù)據(jù)中(圖5)可以看到,未來(lái)100年九洲港可能出現(xiàn)的最高天文潮位接近1.4 m(基于1985國(guó)家高程基準(zhǔn))。

    圖5 未來(lái)100年的最高天文潮位

    5 結(jié)論

    本研究以珠海九洲港為研究對(duì)象,通過(guò)模擬了百年一遇的最強(qiáng)臺(tái)風(fēng)影響下的風(fēng)暴增水以及天文潮位,得出了未來(lái)100年內(nèi)可能出現(xiàn)的最高潮位。對(duì)九洲港來(lái)說(shuō),臺(tái)風(fēng)移動(dòng)方向越接近正北,風(fēng)暴增水越大,路徑為360°的臺(tái)風(fēng)增水模擬最大,當(dāng)中心氣壓和最大風(fēng)速設(shè)為100年一遇的值時(shí),最大增水為5.1 m,九洲港北側(cè)增水要比南側(cè)增水更大。疊加未來(lái)100年可能出現(xiàn)的最高天文潮1.4 m,可以得出基于1985國(guó)家高程基準(zhǔn),未來(lái)100年九州港最高潮位可能達(dá)6.5 m。

    該項(xiàng)工作可作為九洲港海堤高度設(shè)計(jì)的參考,其它碼頭的情況還需再做分析;同時(shí),海堤堤頂高程需綜合考慮設(shè)計(jì)高潮位、波浪爬高以及安全加高值等。另外,風(fēng)暴增水到達(dá)九洲港后,又如何影響下游地區(qū),需更多地考慮地形及建筑等綜合因素。

    猜你喜歡
    九洲風(fēng)暴潮山竹
    廣東省廉江市安鋪中學(xué)九洲江文學(xué)社
    2012年“蘇拉”和“達(dá)維”雙臺(tái)風(fēng)影響的近海風(fēng)暴潮過(guò)程
    1822號(hào)臺(tái)風(fēng)“山竹”演變特征分析
    防范未來(lái)風(fēng)暴潮災(zāi)害的綠色海堤藍(lán)圖
    科學(xué)(2020年4期)2020-11-26 08:27:00
    基于多變量LSTM神經(jīng)網(wǎng)絡(luò)模型的風(fēng)暴潮臨近預(yù)報(bào)
    山竹
    啟蒙(3-7歲)(2019年5期)2019-06-27 07:25:00
    中國(guó)山竹價(jià)格增長(zhǎng)30%
    廣東省風(fēng)暴潮時(shí)空分布特征及重點(diǎn)城市風(fēng)暴潮風(fēng)險(xiǎn)研究
    過(guò)元宵節(jié)啦
    啟蒙(0-3歲)(2012年2期)2012-03-26 06:28:04
    綠色的世界
    99在线视频只有这里精品首页| 欧美乱色亚洲激情| 午夜精品国产一区二区电影| 午夜福利,免费看| а√天堂www在线а√下载| 国产国语露脸激情在线看| 淫妇啪啪啪对白视频| 亚洲专区国产一区二区| av欧美777| 精品免费久久久久久久清纯| 桃色一区二区三区在线观看| 国产精品久久久久成人av| 国产精品综合久久久久久久免费 | 亚洲欧美日韩另类电影网站| 少妇被粗大的猛进出69影院| 久久久久国内视频| 伦理电影免费视频| 日韩中文字幕欧美一区二区| 精品国产美女av久久久久小说| 国产成人精品久久二区二区91| 亚洲一区高清亚洲精品| 久久精品亚洲熟妇少妇任你| 亚洲国产精品一区二区三区在线| 国产单亲对白刺激| 激情在线观看视频在线高清| 搡老岳熟女国产| 亚洲九九香蕉| 免费不卡黄色视频| 中文字幕精品免费在线观看视频| 99re在线观看精品视频| 夜夜爽天天搞| 国产精品一区二区在线不卡| 不卡一级毛片| 女性生殖器流出的白浆| 两人在一起打扑克的视频| 成人精品一区二区免费| 自线自在国产av| 19禁男女啪啪无遮挡网站| 国产色视频综合| 最新美女视频免费是黄的| 九色亚洲精品在线播放| 午夜福利,免费看| 国产精品偷伦视频观看了| 欧美大码av| 亚洲五月色婷婷综合| avwww免费| 亚洲aⅴ乱码一区二区在线播放 | 久久香蕉国产精品| 国产精品美女特级片免费视频播放器 | 高清av免费在线| 一本大道久久a久久精品| 涩涩av久久男人的天堂| 高清黄色对白视频在线免费看| 亚洲第一欧美日韩一区二区三区| 亚洲一区中文字幕在线| 国产一卡二卡三卡精品| 91国产中文字幕| 亚洲情色 制服丝袜| 午夜两性在线视频| 在线天堂中文资源库| 真人一进一出gif抽搐免费| 久久午夜综合久久蜜桃| 国产99白浆流出| 欧美黑人欧美精品刺激| 国产精品秋霞免费鲁丝片| 看片在线看免费视频| 久久久久久人人人人人| av在线天堂中文字幕 | 色哟哟哟哟哟哟| 精品乱码久久久久久99久播| 久99久视频精品免费| 中文字幕av电影在线播放| 老熟妇仑乱视频hdxx| 日韩三级视频一区二区三区| 日本撒尿小便嘘嘘汇集6| 亚洲av成人一区二区三| 在线视频色国产色| 夜夜爽天天搞| 欧美最黄视频在线播放免费 | 嫁个100分男人电影在线观看| 国产aⅴ精品一区二区三区波| 女同久久另类99精品国产91| 久久精品国产清高在天天线| 十八禁人妻一区二区| 精品第一国产精品| 久久人人97超碰香蕉20202| 久久久久久大精品| av网站在线播放免费| 日韩精品中文字幕看吧| 人人妻人人爽人人添夜夜欢视频| 亚洲人成伊人成综合网2020| 欧美日韩瑟瑟在线播放| 久热这里只有精品99| 18美女黄网站色大片免费观看| 国内久久婷婷六月综合欲色啪| 麻豆成人av在线观看| 久久天堂一区二区三区四区| 日韩有码中文字幕| 久久精品亚洲熟妇少妇任你| 日韩欧美一区二区三区在线观看| 99香蕉大伊视频| 看片在线看免费视频| 日本五十路高清| 精品免费久久久久久久清纯| 久久午夜综合久久蜜桃| www.自偷自拍.com| 大型av网站在线播放| 日韩欧美在线二视频| 又紧又爽又黄一区二区| 亚洲精品美女久久av网站| 国产精品美女特级片免费视频播放器 | 欧美人与性动交α欧美软件| 国产欧美日韩一区二区三| 在线观看免费日韩欧美大片| 亚洲精品久久午夜乱码| 村上凉子中文字幕在线| 国产精品 国内视频| 久久久久精品国产欧美久久久| 亚洲 欧美 日韩 在线 免费| 午夜福利免费观看在线| 亚洲精品国产一区二区精华液| 国产亚洲av高清不卡| 深夜精品福利| 国产区一区二久久| 水蜜桃什么品种好| 国产一区二区在线av高清观看| 精品一区二区三区四区五区乱码| 曰老女人黄片| 老汉色av国产亚洲站长工具| 女人精品久久久久毛片| 黄网站色视频无遮挡免费观看| 天堂动漫精品| 久久精品国产99精品国产亚洲性色 | 美女福利国产在线| 精品久久久久久,| 一级片免费观看大全| 亚洲精品国产色婷婷电影| 国产无遮挡羞羞视频在线观看| 美女高潮到喷水免费观看| 亚洲免费av在线视频| 亚洲午夜精品一区,二区,三区| a级片在线免费高清观看视频| 国产熟女午夜一区二区三区| 亚洲欧美一区二区三区久久| 一区福利在线观看| 99热国产这里只有精品6| 日韩免费高清中文字幕av| 在线观看免费视频网站a站| 91国产中文字幕| 男男h啪啪无遮挡| 在线视频色国产色| 国产成人精品久久二区二区91| 成人影院久久| 日本五十路高清| 久久亚洲真实| 99国产精品99久久久久| 久久精品国产亚洲av香蕉五月| 最近最新免费中文字幕在线| 欧美最黄视频在线播放免费 | 日本一区二区免费在线视频| 国产精品一区二区精品视频观看| 啦啦啦 在线观看视频| 欧美丝袜亚洲另类 | 他把我摸到了高潮在线观看| 性少妇av在线| 亚洲三区欧美一区| 欧美激情高清一区二区三区| 黄色丝袜av网址大全| 国产欧美日韩一区二区三区在线| 电影成人av| 黄色怎么调成土黄色| 亚洲 国产 在线| 狠狠狠狠99中文字幕| 国产日韩一区二区三区精品不卡| 波多野结衣高清无吗| 人成视频在线观看免费观看| 麻豆一二三区av精品| 国产成人系列免费观看| 国产色视频综合| netflix在线观看网站| 久热这里只有精品99| 欧美日韩福利视频一区二区| 亚洲精品一卡2卡三卡4卡5卡| 久久精品影院6| 性欧美人与动物交配| 亚洲五月色婷婷综合| 午夜福利免费观看在线| xxx96com| 国产野战对白在线观看| 69精品国产乱码久久久| 亚洲专区国产一区二区| 日韩视频一区二区在线观看| 国产成人av激情在线播放| 国产野战对白在线观看| 91国产中文字幕| 香蕉久久夜色| 国产精品免费一区二区三区在线| 欧美一级毛片孕妇| 男人舔女人的私密视频| 中国美女看黄片| 亚洲中文日韩欧美视频| 国产精品免费一区二区三区在线| 人妻久久中文字幕网| 动漫黄色视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 亚洲一码二码三码区别大吗| 美女福利国产在线| 国产免费av片在线观看野外av| 国产av又大| 欧美乱色亚洲激情| 久久精品影院6| 成人亚洲精品av一区二区 | 国产精品亚洲av一区麻豆| 午夜两性在线视频| 香蕉丝袜av| 最近最新中文字幕大全电影3 | 又黄又粗又硬又大视频| 黄色视频不卡| 成人亚洲精品av一区二区 | av网站在线播放免费| 99久久久亚洲精品蜜臀av| а√天堂www在线а√下载| 久久青草综合色| 美女午夜性视频免费| 欧美色视频一区免费| 热99国产精品久久久久久7| 深夜精品福利| www.熟女人妻精品国产| 国产又爽黄色视频| 操美女的视频在线观看| 欧美日韩黄片免| 在线看a的网站| 看免费av毛片| 国产一卡二卡三卡精品| 国产成人啪精品午夜网站| 大型黄色视频在线免费观看| 黄片大片在线免费观看| 不卡av一区二区三区| 国产精品二区激情视频| 国内久久婷婷六月综合欲色啪| 国产精品影院久久| 成人黄色视频免费在线看| 欧美国产精品va在线观看不卡| 国产欧美日韩精品亚洲av| 欧美另类亚洲清纯唯美| 好男人电影高清在线观看| 一级,二级,三级黄色视频| 欧美日韩国产mv在线观看视频| 黄片小视频在线播放| 久久草成人影院| 三上悠亚av全集在线观看| 久久天躁狠狠躁夜夜2o2o| 好看av亚洲va欧美ⅴa在| 美女高潮喷水抽搐中文字幕| 免费av中文字幕在线| 亚洲美女黄片视频| 免费在线观看亚洲国产| 18禁观看日本| 午夜精品久久久久久毛片777| 老熟妇乱子伦视频在线观看| 日本免费a在线| 亚洲色图综合在线观看| 99热只有精品国产| 多毛熟女@视频| 国产三级在线视频| 后天国语完整版免费观看| 日韩精品中文字幕看吧| 久久人妻福利社区极品人妻图片| 黄色怎么调成土黄色| 男女做爰动态图高潮gif福利片 | 精品一区二区三区av网在线观看| 午夜免费观看网址| 国产精品永久免费网站| 男人舔女人下体高潮全视频| 婷婷丁香在线五月| 精品久久久久久,| 国产欧美日韩综合在线一区二区| 亚洲国产欧美日韩在线播放| 国产精品 国内视频| 午夜福利影视在线免费观看| 午夜福利在线免费观看网站| 丁香六月欧美| 99在线人妻在线中文字幕| 日本三级黄在线观看| 久久人人精品亚洲av| 麻豆国产av国片精品| 亚洲专区国产一区二区| 亚洲国产欧美一区二区综合| 岛国在线观看网站| 亚洲国产精品一区二区三区在线| 欧美 亚洲 国产 日韩一| 久久久国产精品麻豆| 夜夜躁狠狠躁天天躁| 日韩欧美一区二区三区在线观看| √禁漫天堂资源中文www| 欧美av亚洲av综合av国产av| 欧美黄色淫秽网站| 男女午夜视频在线观看| aaaaa片日本免费| 午夜视频精品福利| 国产成人影院久久av| 热99国产精品久久久久久7| 国产精品爽爽va在线观看网站 | 视频区图区小说| 妹子高潮喷水视频| 国产精品成人在线| 热re99久久精品国产66热6| 一级片免费观看大全| 90打野战视频偷拍视频| 国产99久久九九免费精品| 欧美老熟妇乱子伦牲交| 亚洲avbb在线观看| 国产成人系列免费观看| 成年女人毛片免费观看观看9| 婷婷精品国产亚洲av在线| av福利片在线| 久久久久久人人人人人| 麻豆国产av国片精品| 中文亚洲av片在线观看爽| 亚洲aⅴ乱码一区二区在线播放 | 99香蕉大伊视频| 激情视频va一区二区三区| 最近最新免费中文字幕在线| 亚洲专区国产一区二区| 久久精品成人免费网站| 久久香蕉精品热| 夜夜爽天天搞| 国产精品99久久99久久久不卡| 欧美成人性av电影在线观看| 欧美日韩瑟瑟在线播放| 人成视频在线观看免费观看| 亚洲色图av天堂| 一夜夜www| 午夜福利影视在线免费观看| 国产伦人伦偷精品视频| 国产精品久久久久成人av| 久久中文字幕人妻熟女| netflix在线观看网站| 琪琪午夜伦伦电影理论片6080| 我的亚洲天堂| 免费高清视频大片| 一级毛片女人18水好多| 亚洲精品粉嫩美女一区| 国产激情欧美一区二区| 91精品三级在线观看| 两个人看的免费小视频| 免费在线观看完整版高清| 久久久国产成人精品二区 | 成人影院久久| 黄片小视频在线播放| av天堂久久9| 少妇被粗大的猛进出69影院| 制服人妻中文乱码| 亚洲,欧美精品.| 搡老乐熟女国产| 国产精品国产av在线观看| 国产激情久久老熟女| 欧美在线黄色| 性少妇av在线| 日韩成人在线观看一区二区三区| 久久久国产成人精品二区 | 久久国产精品人妻蜜桃| 99热只有精品国产| 国产99久久九九免费精品| 日本黄色日本黄色录像| 国产成人av激情在线播放| 亚洲久久久国产精品| 婷婷丁香在线五月| 国产深夜福利视频在线观看| 视频区图区小说| 丝袜在线中文字幕| 日日干狠狠操夜夜爽| 精品久久久久久久久久免费视频 | 高清av免费在线| 男人舔女人下体高潮全视频| 高清在线国产一区| 黄片小视频在线播放| 亚洲情色 制服丝袜| 99精品欧美一区二区三区四区| 9色porny在线观看| 在线免费观看的www视频| 久久人妻av系列| 欧美老熟妇乱子伦牲交| 看免费av毛片| 精品久久久久久久毛片微露脸| 精品国产国语对白av| 黑丝袜美女国产一区| 亚洲精品国产色婷婷电影| 男女午夜视频在线观看| 精品午夜福利视频在线观看一区| 男女下面进入的视频免费午夜 | 日韩高清综合在线| 999久久久精品免费观看国产| 欧美精品一区二区免费开放| 999久久久国产精品视频| 国产成人一区二区三区免费视频网站| 国产野战对白在线观看| 欧美日韩亚洲国产一区二区在线观看| 少妇被粗大的猛进出69影院| 精品一区二区三区av网在线观看| 美女高潮到喷水免费观看| 精品高清国产在线一区| 叶爱在线成人免费视频播放| 在线十欧美十亚洲十日本专区| а√天堂www在线а√下载| 欧美不卡视频在线免费观看 | 伊人久久大香线蕉亚洲五| 久久这里只有精品19| 最好的美女福利视频网| 婷婷丁香在线五月| 午夜免费观看网址| 亚洲 欧美一区二区三区| 视频区欧美日本亚洲| 亚洲av成人一区二区三| 亚洲精品中文字幕一二三四区| 欧美不卡视频在线免费观看 | 神马国产精品三级电影在线观看 | 国产91精品成人一区二区三区| 国产在线精品亚洲第一网站| 琪琪午夜伦伦电影理论片6080| netflix在线观看网站| 免费女性裸体啪啪无遮挡网站| www日本在线高清视频| 色精品久久人妻99蜜桃| 欧美亚洲日本最大视频资源| 国产精品秋霞免费鲁丝片| 色在线成人网| 亚洲男人的天堂狠狠| 午夜福利欧美成人| 级片在线观看| 欧美日本中文国产一区发布| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品中文字幕一二三四区| 97人妻天天添夜夜摸| 日韩大码丰满熟妇| 国产精品综合久久久久久久免费 | 在线播放国产精品三级| 亚洲精品美女久久av网站| 亚洲av成人一区二区三| 亚洲aⅴ乱码一区二区在线播放 | 日本一区二区免费在线视频| 男女午夜视频在线观看| 国产亚洲欧美98| 一级a爱片免费观看的视频| 久久伊人香网站| 精品乱码久久久久久99久播| 亚洲中文日韩欧美视频| 欧美成人性av电影在线观看| 一区二区三区激情视频| 国产精品美女特级片免费视频播放器 | 国产又色又爽无遮挡免费看| 最新美女视频免费是黄的| 久久精品影院6| 精品国产国语对白av| 最近最新中文字幕大全免费视频| 一区二区三区激情视频| x7x7x7水蜜桃| 午夜福利欧美成人| 淫妇啪啪啪对白视频| 91成年电影在线观看| 国产精品久久久久成人av| 色精品久久人妻99蜜桃| ponron亚洲| 性欧美人与动物交配| 久久精品国产亚洲av香蕉五月| 我的亚洲天堂| 国内毛片毛片毛片毛片毛片| 亚洲专区国产一区二区| 操出白浆在线播放| 国产真人三级小视频在线观看| av天堂久久9| 亚洲中文av在线| 亚洲自偷自拍图片 自拍| 中国美女看黄片| 久久久久久人人人人人| 亚洲第一欧美日韩一区二区三区| 日韩欧美三级三区| 91在线观看av| 最近最新免费中文字幕在线| 久久人人97超碰香蕉20202| 国产成人欧美在线观看| 亚洲精品一二三| 中文字幕高清在线视频| 啪啪无遮挡十八禁网站| 欧美日韩亚洲高清精品| 国产单亲对白刺激| 99热只有精品国产| 国产熟女xx| 亚洲国产欧美一区二区综合| 国产欧美日韩一区二区三区在线| 色综合欧美亚洲国产小说| 免费在线观看完整版高清| 国产亚洲精品综合一区在线观看 | 国产高清视频在线播放一区| 可以免费在线观看a视频的电影网站| 免费av毛片视频| 国产激情欧美一区二区| 亚洲国产精品一区二区三区在线| 不卡一级毛片| 可以免费在线观看a视频的电影网站| 久久天堂一区二区三区四区| 在线国产一区二区在线| 色精品久久人妻99蜜桃| 成人国产一区最新在线观看| 亚洲男人的天堂狠狠| 老汉色av国产亚洲站长工具| 久久精品亚洲熟妇少妇任你| 国产av又大| 久久国产乱子伦精品免费另类| 99久久综合精品五月天人人| 亚洲成av片中文字幕在线观看| 女人爽到高潮嗷嗷叫在线视频| 国产精品久久视频播放| 国产深夜福利视频在线观看| 丰满人妻熟妇乱又伦精品不卡| av有码第一页| 老熟妇乱子伦视频在线观看| 在线观看免费午夜福利视频| 久久狼人影院| 国产精品野战在线观看 | 国产成人欧美| 99久久综合精品五月天人人| 老汉色av国产亚洲站长工具| 久久久久国产精品人妻aⅴ院| 久久性视频一级片| 成人三级黄色视频| 夜夜躁狠狠躁天天躁| 国产伦人伦偷精品视频| 亚洲精品中文字幕在线视频| 久久久久国产一级毛片高清牌| 亚洲成av片中文字幕在线观看| 亚洲专区字幕在线| 成人手机av| 亚洲伊人色综图| 黑人猛操日本美女一级片| 琪琪午夜伦伦电影理论片6080| 精品福利观看| 国产精品秋霞免费鲁丝片| 国产蜜桃级精品一区二区三区| 日日爽夜夜爽网站| 日韩av在线大香蕉| 成人av一区二区三区在线看| 校园春色视频在线观看| 在线天堂中文资源库| 激情在线观看视频在线高清| bbb黄色大片| 国产蜜桃级精品一区二区三区| 午夜日韩欧美国产| 1024视频免费在线观看| 女人高潮潮喷娇喘18禁视频| 黑人欧美特级aaaaaa片| 日本免费a在线| 老汉色∧v一级毛片| 色综合欧美亚洲国产小说| 国产一区二区在线av高清观看| 欧美日韩精品网址| 很黄的视频免费| a级片在线免费高清观看视频| av电影中文网址| 久久人妻福利社区极品人妻图片| 久久人人爽av亚洲精品天堂| 亚洲av第一区精品v没综合| 12—13女人毛片做爰片一| 99久久综合精品五月天人人| 亚洲精品粉嫩美女一区| 黄色成人免费大全| 操美女的视频在线观看| 在线观看一区二区三区| 久久久久久久久中文| 99riav亚洲国产免费| 在线视频色国产色| av国产精品久久久久影院| 日韩欧美免费精品| 黄色 视频免费看| 国产激情久久老熟女| 国产99白浆流出| 99久久综合精品五月天人人| 少妇 在线观看| 国产一卡二卡三卡精品| 久久草成人影院| 欧美+亚洲+日韩+国产| 国产精华一区二区三区| 欧美国产精品va在线观看不卡| 国产精品日韩av在线免费观看 | 久久热在线av| 午夜福利欧美成人| 女人精品久久久久毛片| 精品人妻1区二区| 久99久视频精品免费| 男女午夜视频在线观看| 久久人人爽av亚洲精品天堂| 99精国产麻豆久久婷婷| 激情视频va一区二区三区| 热99re8久久精品国产| 欧美午夜高清在线| 精品国产国语对白av| 91国产中文字幕| 欧美在线黄色| 一级毛片高清免费大全| 1024视频免费在线观看| 午夜免费成人在线视频| 久久人人爽av亚洲精品天堂| 欧美在线一区亚洲| 韩国精品一区二区三区| 国产高清国产精品国产三级| 国产精品av久久久久免费| 久久久久久免费高清国产稀缺| 人妻久久中文字幕网| 五月开心婷婷网| 侵犯人妻中文字幕一二三四区| 免费在线观看黄色视频的| 国产免费av片在线观看野外av| 9191精品国产免费久久| 成人黄色视频免费在线看| av在线播放免费不卡|