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

    港口及海岸工程防災(zāi)理論體系的研發(fā)及其應(yīng)用?

    2018-12-22 05:28:33劉德輔劉桂林王鳳清王永發(fā)
    關(guān)鍵詞:颶風(fēng)極值臺(tái)風(fēng)

    劉德輔, 劉桂林, 王鳳清, 王永發(fā)

    (1.中國(guó)海洋大學(xué)工程學(xué)院,山東 青島 266003; 2.青島新前灣集裝箱碼頭公司,山東 青島 266555)

    港口及海岸工程防災(zāi)理論體系的研發(fā)及其應(yīng)用?

    劉德輔1, 劉桂林1, 王鳳清1, 王永發(fā)2

    (1.中國(guó)海洋大學(xué)工程學(xué)院,山東 青島 266003; 2.青島新前灣集裝箱碼頭公司,山東 青島 266555)

    本文介紹了中國(guó)自主研發(fā)的復(fù)合極值分布理論和多維復(fù)合極值理論30多年來(lái)在國(guó)內(nèi)外的引用和應(yīng)用,特別是2005年Katrina和2012年Sandy颶風(fēng)對(duì)新奧爾良市和新澤西州及紐約周邊帶來(lái)的災(zāi)難性破壞,以及2013年菲特臺(tái)風(fēng)誘發(fā)上海超警戒水位,都驗(yàn)證了1982及2006年使用上述理論預(yù)測(cè)結(jié)果的正確性。本文就該理論在海洋,海岸、港口等工程防災(zāi)設(shè)防標(biāo)準(zhǔn)中的應(yīng)用,對(duì)國(guó)內(nèi)外有關(guān)規(guī)范的防災(zāi)標(biāo)準(zhǔn)進(jìn)行了風(fēng)險(xiǎn)評(píng)估和建議。

    臺(tái)風(fēng)災(zāi)害;復(fù)合極值分布;多維復(fù)合極值分布;海岸和港口工程;規(guī)范校核

    隨著全球氣候變暖、海平面上升等因素加劇的趨勢(shì),臺(tái)風(fēng)、巨浪、暴潮、暴雨、洪水巨災(zāi)已成為影響海岸、海洋、水利以及港口工程安全的重大問(wèn)題。世界氣象組織指出:百分之九十的自然災(zāi)害為極端氣象災(zāi)害,例如臺(tái)風(fēng)、颶風(fēng),熱帶風(fēng)暴災(zāi)害。1975年尼娜臺(tái)風(fēng)在福建登陸,經(jīng)歷108 h北上至河南省,臺(tái)風(fēng)帶來(lái)的暖濕氣流遇到了北方的冷空氣,導(dǎo)致3天之內(nèi)降水量達(dá)到1 m,板橋水庫(kù)及下游62個(gè)水庫(kù)的潰決,17萬(wàn)多人死亡,1 200多萬(wàn)人受災(zāi),百萬(wàn)公頃農(nóng)田被淹,京廣鐵路100多公里被毀。2006年桑美和碧利斯兩場(chǎng)臺(tái)風(fēng)共導(dǎo)致1 600多人死亡和失蹤,5 000多億經(jīng)濟(jì)損失。桑美臺(tái)風(fēng)誘發(fā)了7 m巨浪和3.8 m風(fēng)暴增水,導(dǎo)致沙埕港952只船舶沉沒(méi),1 594只船舶毀壞。如果臺(tái)風(fēng)延后2 h登陸,則恰逢天文大潮,這種“三碰頭”的組合可完全淹沒(méi)包括若干個(gè)核電站在內(nèi)的福建、浙江大部分土地,其災(zāi)難性后果難以估量。這種小概率極端事件的發(fā)生是完全可能的。2013年10月8日臺(tái)風(fēng)菲特導(dǎo)致暴雨洪水、風(fēng)暴潮,受天文大潮頂托的共同影響,黃浦江沿線(xiàn)水位上漲迅速,下午2時(shí)05分,吳淞口實(shí)測(cè)水位達(dá)5.15 m,超過(guò)按國(guó)家規(guī)范計(jì)算的警戒水位0.35 m線(xiàn)[1]。

    長(zhǎng)期以來(lái),由于概率統(tǒng)計(jì)理論和方法上的限制,在港口、海岸、近海及河口城市防護(hù)工程的設(shè)計(jì)仍停留在對(duì)實(shí)測(cè)資料年最大值概率分析基礎(chǔ)上,假定年最大值服從某種分布(Gumbel,Weibull或P-Ⅲ型分布等),估計(jì)不同重現(xiàn)期對(duì)應(yīng)的設(shè)計(jì)波高、水位、風(fēng)速作為設(shè)計(jì)標(biāo)準(zhǔn),對(duì)惡劣環(huán)境條件下多種荷載共同作用組合的問(wèn)題一直未得到妥善解決。一些沒(méi)有理論依據(jù)的模式和規(guī)定影響到概率預(yù)測(cè)的正確性和防災(zāi)標(biāo)準(zhǔn)的合理性,導(dǎo)致極端氣候因素誘發(fā)的潰壩,決堤巨災(zāi)時(shí)有發(fā)生[2]。

    1 復(fù)合極值分布理論(CEVD)及多維復(fù)合極值分布理論(MCEVD) 的建立

    自1978年代初交通部“海港水文規(guī)范編寫(xiě)組”成立以來(lái),進(jìn)行了全國(guó)各地港口的調(diào)查研究,特別是1972年3號(hào)臺(tái)風(fēng)重創(chuàng)大連港,有效波高6.8 m,超過(guò)了1896年的臺(tái)風(fēng)巨災(zāi),導(dǎo)致大連海域的空前損失,更進(jìn)一步加深了對(duì)臺(tái)風(fēng)波浪長(zhǎng)期概率預(yù)測(cè)重要性的認(rèn)識(shí),啟發(fā)我們開(kāi)展了針對(duì)臺(tái)風(fēng)特征進(jìn)行臺(tái)風(fēng)影響海區(qū)設(shè)計(jì)波浪(及其他致災(zāi)因素)概率預(yù)測(cè)理論和計(jì)算方法的研究。

    鑒于國(guó)內(nèi)外慣用年極值系列按照某種概率模式(如P-Ⅲ型、Gumbel,Weibull分布)外延推求不同重現(xiàn)期設(shè)計(jì)值的作法存在的弊端——遺漏了每年各次臺(tái)風(fēng)致災(zāi)因素提供的重要信息,影響到概率預(yù)測(cè)結(jié)果的置信度,考慮到臺(tái)風(fēng)影響不同海區(qū)每年出現(xiàn)的頻次各不相同,可用離散型隨機(jī)變量及相應(yīng)的離散型概率分布表示;臺(tái)風(fēng)誘發(fā)的致災(zāi)因素(風(fēng)速、增水、波高等)可用連續(xù)型極值分布表示,二者的組合在順序統(tǒng)計(jì)學(xué)和測(cè)度論的理論基礎(chǔ)上經(jīng)過(guò)嚴(yán)格的數(shù)學(xué)推導(dǎo),可構(gòu)成一種新的極值分布模式——復(fù)合極值分布(簡(jiǎn)稱(chēng)CEVD)。根據(jù)統(tǒng)計(jì)檢驗(yàn)結(jié)果,臺(tái)風(fēng)(颶風(fēng))影響各海區(qū)的頻次符合Poisson分布,臺(tái)風(fēng)(颶風(fēng))波浪的不同特征可分別采用Gumbel和Weibull分布,可構(gòu)成Poisson-Gumbel CEVD分布(用于中國(guó)臺(tái)風(fēng)影響海域)[3]和Poisson-Weibull CEVD分布(用于美國(guó)墨西哥灣和大西洋沿岸)[4]。

    1.1 Poisson-Gumbel復(fù)合極值分布(CEVD)

    當(dāng)臺(tái)風(fēng)頻次為泊松分布,波高符合Gumbel分布時(shí),可簡(jiǎn)化為:

    HP=u+XP/α。

    (1)

    1.2 多維復(fù)合極值分布(MCEVD)的理論求解及應(yīng)用

    在上述復(fù)合極值分布理論的基礎(chǔ)上進(jìn)一步建立了一個(gè)新的理論分布模式:多維復(fù)合極值分布模式,既考慮了不同海區(qū)每年臺(tái)風(fēng)(颶風(fēng),熱帶風(fēng)暴)出現(xiàn)的頻次,又考慮了它們誘發(fā)的多種極端海況(風(fēng)速、波高、周期、風(fēng)暴潮、海流等)出現(xiàn)的概率特性。前者可構(gòu)成一種離散型分布,后者則可用不同的多維聯(lián)合極值分布表達(dá),從而構(gòu)成離散型分布與多維聯(lián)合極值分布組合而成的新分布—多維復(fù)合極值分布。以三維復(fù)合極值分布為例,經(jīng)過(guò)嚴(yán)格的數(shù)學(xué)推導(dǎo),當(dāng)臺(tái)風(fēng)頻次為Poisson分布時(shí),三維泊松復(fù)合極值分布的表達(dá)式如下:

    (2)

    式中:g(x,y,z)可為任何一種三維聯(lián)合概率密度函數(shù);G1(u)為分布函數(shù)。

    新理論自1980年在國(guó)外重要學(xué)術(shù)刊物公開(kāi)發(fā)表后,得到了國(guó)內(nèi)外學(xué)術(shù)界和工程界的重視、引用和應(yīng)用。30多年來(lái),CEVD模式不僅正式列入了中國(guó)教育部港工專(zhuān)業(yè)、海工專(zhuān)業(yè)全國(guó)統(tǒng)編教材的正式內(nèi)容,并在近40項(xiàng)工程項(xiàng)目的設(shè)計(jì)波浪推求中得到應(yīng)用。 在國(guó)家防汛抗旱總指揮部辦公室委托項(xiàng)目“中國(guó)沿海臺(tái)風(fēng)災(zāi)害區(qū)劃、防臺(tái)風(fēng)標(biāo)準(zhǔn)應(yīng)急制定、防臺(tái)風(fēng)應(yīng)急評(píng)估標(biāo)準(zhǔn)制定”中,研發(fā)了基于多維復(fù)合極值分布理論的“雙層嵌套多目標(biāo)聯(lián)合概率模式”(簡(jiǎn)稱(chēng)DLNMOPM),用于沿海各區(qū)防臺(tái)風(fēng)設(shè)防標(biāo)準(zhǔn)的概率預(yù)測(cè)及災(zāi)害區(qū)劃,2008年經(jīng)國(guó)家防辦組織專(zhuān)家驗(yàn)收,下達(dá)執(zhí)行。同年,交通部在國(guó)家行業(yè)標(biāo)準(zhǔn)“海港水文規(guī)范”(JTJ213-98)修訂稿中(2013年出版)正式采用了“復(fù)合極值分布理論”作為設(shè)計(jì)波高概率預(yù)測(cè)的理論和方法,首次取代了該國(guó)家行業(yè)標(biāo)準(zhǔn)中持續(xù)使用近四十年的P-Ⅲ型曲線(xiàn)[5-9];此外,在美國(guó)墨西哥灣及大西洋沿岸、加拿大、尼羅河三角洲的波候、地中海及新加坡的波況和風(fēng)況,挪威和韓國(guó)極端海況研究等多項(xiàng)工程中都得到應(yīng)用[10-18]。

    1.3 2005 、2012美國(guó)兩場(chǎng)颶風(fēng)重災(zāi)及2013中國(guó)菲特臺(tái)風(fēng)災(zāi)害,驗(yàn)證了CEVD和MCEVD 1982和2006預(yù)測(cè)結(jié)果的正確性

    2005年美國(guó)發(fā)生的Katrina颶風(fēng)重災(zāi),對(duì)新奧爾良的毀滅性破壞,導(dǎo)致新奧爾良市毀滅性破壞,1 800人死亡,800億美元經(jīng)濟(jì)損失。颶風(fēng)的強(qiáng)度和特征不僅驗(yàn)證了使用復(fù)合極值分布理論模式1982年預(yù)測(cè)結(jié)果的正確性(見(jiàn)表1),也顯示了2005災(zāi)后使用MCEVD復(fù)核結(jié)果的合理性(見(jiàn)文獻(xiàn)[19]中圖8),同時(shí)也表明美國(guó)國(guó)家海洋和大氣管理局(NOAA)制訂的上述海區(qū)防災(zāi)標(biāo)準(zhǔn)—可能最大颶風(fēng)(PMH)和標(biāo)準(zhǔn)設(shè)計(jì)颶風(fēng)(SPH)過(guò)低,是巨災(zāi)的主要誘因,這一事實(shí)已得到公認(rèn)[1-2]。2008年美國(guó)佛羅里達(dá)國(guó)際颶風(fēng)研究中心在海岸防護(hù)颶風(fēng)災(zāi)害工程設(shè)計(jì)中,引用了作者的4篇論文[1,3-5]正式用作防護(hù)颶風(fēng)災(zāi)害設(shè)計(jì)標(biāo)準(zhǔn)的依據(jù)[20]

    2012年10月30—31日桑迪颶風(fēng)重現(xiàn)堤壩沖毀,特拉華州,新澤西州及紐約周邊大面積淹沒(méi),大量人員傷亡,數(shù)百億美元經(jīng)濟(jì)損失,再次驗(yàn)證了傳統(tǒng)的設(shè)計(jì)標(biāo)準(zhǔn)和工程結(jié)構(gòu)是致災(zāi)的主要誘因。作者對(duì)上述海區(qū)颶風(fēng)強(qiáng)度預(yù)測(cè)成果再次得到證實(shí)(見(jiàn)文獻(xiàn)[19]圖8)。

    表1 墨西哥灣A區(qū)卡特里娜颶風(fēng)強(qiáng)度概率分析(按文獻(xiàn)[19]圖8計(jì)算)Table 1 Comparison between NOAA and CEVD predicted value (by [19],Fig.8 calculated) /hPa

    Note:①Zone number;②Hurrieane example;③Gulf of Mexico Zone A

    颶風(fēng)的強(qiáng)度和特征,驗(yàn)證了使用復(fù)合極值分布理論1982年預(yù)測(cè)結(jié)果的正確性(見(jiàn)表1)。圖1則為作者使用2002年以來(lái)新開(kāi)發(fā)的MCEVD分析颶風(fēng)及其誘發(fā)的災(zāi)害性海況,顯示2005年以前美國(guó)不同學(xué)者公開(kāi)發(fā)表對(duì)上述海域預(yù)測(cè)成果,明顯低于MCEVD的預(yù)測(cè)值(見(jiàn)文獻(xiàn)[9]中圖6)。MCEVD的預(yù)測(cè)結(jié)果顯示,重現(xiàn)期100年一遇的風(fēng)速3區(qū)(新奧爾良區(qū))與2005卡特里娜颶風(fēng)狀況相近,明顯高于設(shè)防標(biāo)準(zhǔn);8-9區(qū)(新澤西,紐約沿岸)則與2012年桑迪颶風(fēng)最大風(fēng)速相近 (見(jiàn)表2)。

    圖1 100年一遇的極端風(fēng)速對(duì)比PNLTCED 及美國(guó)學(xué)者的概率預(yù)測(cè)[19.圖.8]Fig.1 100 years return period for extreme wind velocities by different predicted results

    同樣,如圖2所示,1982年預(yù)測(cè)的風(fēng)暴增水(見(jiàn)文獻(xiàn)[4]中的圖8),費(fèi)城周邊(圖中虛線(xiàn))100年一遇風(fēng)暴增水達(dá)到10Ft(3 m),完全為2012年桑迪颶風(fēng)對(duì)新澤西,紐約,特拉華州大面積淹沒(méi)的事實(shí)所證實(shí)。

    2 多維復(fù)合極值分布理論在河口、海岸城市防災(zāi)設(shè)計(jì)標(biāo)準(zhǔn)中的應(yīng)用

    多維復(fù)合極值分布模型在推算上海市防洪設(shè)計(jì)標(biāo)準(zhǔn)中的應(yīng)用:上海地處長(zhǎng)江入海口處,防洪水位極受重視,有關(guān)部門(mén)多次提高防洪標(biāo)準(zhǔn),現(xiàn)采用千年一遇水位5.86 m。針對(duì)‘中國(guó)臺(tái)風(fēng)年鑒’1949—1987年資料,吳淞受臺(tái)風(fēng)影響期間的實(shí)測(cè)水位經(jīng)過(guò)風(fēng)險(xiǎn)分析得出,在極端環(huán)境條件下影響上海市防洪水位的主要因素為長(zhǎng)江下泄洪峰和風(fēng)暴增水與天文大潮同時(shí)發(fā)生的“三碰頭”組合。計(jì)算結(jié)果考慮臺(tái)風(fēng)暴潮,暴雨洪水與天文大潮的組合,2006年使用MCEVD計(jì)算的百年一遇預(yù)測(cè)結(jié)果如表3、4所示,已超過(guò)防護(hù)工程500年一遇超警戒水位設(shè)防標(biāo)準(zhǔn)。10月8日,受天文大潮頂托、上游洪水下泄,本地大暴雨和臺(tái)風(fēng)暴潮的共同影響,水位上漲迅速,14:05吳淞實(shí)測(cè)水位達(dá)5.15 m,超過(guò)警戒水位0.35 m。

    表2 不同方法預(yù)測(cè)結(jié)果比較(見(jiàn)文獻(xiàn)[19].圖.8)Table 2 Comparison Predicted Results by Different Mrthods([19],Fig.8)

    Note:①Comparison Predicted Results by Different Marthods;②Predicted 100-yrs return period wind speed for New Orleans(3 zone);③Predicted 100-yrs return period wind speed for New Jersey and New York coastal areas(9 zone)

    圖2 美國(guó)大西洋沿岸各區(qū)不同重現(xiàn)期風(fēng)暴增水1982年預(yù)測(cè)結(jié)果([4]中圖8)

    表3 吳淞港防洪水位Table 3 Design water level at Wuson datum plane

    Note:①Joint return period;②Flood;③Storm surge;④Astronomical tide;⑤Design water level

    由表可見(jiàn),按單因素外延的千年一遇上海防洪水位5.86 m僅相當(dāng)于考慮臺(tái)風(fēng)頻次,風(fēng)暴潮、天文大潮、長(zhǎng)江上游洪水的聯(lián)合重現(xiàn)期百年一遇值。正如Kirby和Moss在“美國(guó)洪水概率分析總結(jié)”一文中指出:“將基于事件取樣和聯(lián)合概率結(jié)合起來(lái)將大大推動(dòng)特大洪水概率預(yù)測(cè)的發(fā)展”[10], 而MCEVD理論的立論依據(jù)與Kirby和Moss所期待的特大洪水預(yù)測(cè)模式發(fā)展方向不謀而合。 即每年臺(tái)風(fēng)過(guò)程取樣較傳統(tǒng)理論方法的年最大值取樣,統(tǒng)計(jì)數(shù)據(jù)遠(yuǎn)遠(yuǎn)超過(guò)傳統(tǒng)取樣方法,概率預(yù)測(cè)結(jié)果的置信區(qū)間遠(yuǎn)遠(yuǎn)小于年最大值取樣預(yù)測(cè)值的置信區(qū)間,由此可見(jiàn),應(yīng)用復(fù)合極值理論預(yù)測(cè)結(jié)果在相應(yīng)的置信區(qū)間范圍內(nèi),必然高于傳統(tǒng)理論年最大值預(yù)測(cè)結(jié)果。因此, MCEVD預(yù)測(cè)結(jié)果較傳統(tǒng)方法更為合理的結(jié)論也就不言而喻了。

    表4 臺(tái)風(fēng)暴潮,暴雨洪水與天文大潮MCEVD和傳統(tǒng)外延法預(yù)測(cè)結(jié)果比較Table 4 Cpmparison between MCEVD and traditional method predicted extreme water level for Shanghai

    Note:①M(fèi)odel;②Return period;③Design value;④Tradional method(P-3 model);⑤Shanghai warning water level(P-3 model);⑥Typhoon Fitow induced Wuson port observd sea level

    3 建議

    (1)Katrina颶風(fēng)帶來(lái)巨大災(zāi)害的教訓(xùn),對(duì)海岸,港口工程的防災(zāi)設(shè)計(jì)標(biāo)準(zhǔn)的確定,特別是濱海核電站防護(hù)工程的安全非常重要,全盤(pán)翻譯引用國(guó)外規(guī)范的隱患,必須引以為鑒[19]。

    (2)鑒于沿海港區(qū)大風(fēng)導(dǎo)致港內(nèi)集裝箱垜倒塌事故時(shí)有發(fā)生,例如:2014年8月鹽田港幾十個(gè)集裝箱遭臺(tái)風(fēng)襲擊,很多民房破壞,造成500多萬(wàn)損失:2010年8月蕪湖港臺(tái)風(fēng)導(dǎo)致兩座門(mén)吊倒塌,造成2死1傷:2005年Octavia輪??繉幉ū眮龈?,溫帶氣旋襲擊,波高7~8 m,造成集裝箱損失。建議針對(duì)所在海域臺(tái)風(fēng)或極端海況可能誘發(fā)影響集裝箱垜安全的外荷研究和制訂港區(qū)安全規(guī)范。

    [1] Liu D F, Pang L, Xie B T. Typhoon disaster zoning and prevention criteria—A double layer nested multi-objective probability model and its application[J]. Science in China (E), 2008, 51(7): 1038-1048.

    [2] Liu D F, Li H J, Liu G L, et al. Book Title: Natural Disasters: Prevention, Risk Factors and Management[M]. USA: NOVA Science Publishers, Inc. Chapter 1: Typhoon/ Hurricane/ Tropical Cyclone Disasters: Prediction USA: Prevention and Mitigation, 2013: 1-72.

    [3] Liu D F, Ma F S. Prediction of extreme wave heights and wind velocities[J]. Journal of the Waterway Port Coastal and Ocean Division, ASCE, 1980, (WW4): 469-479.

    [4] Liu D F, Long Term Distributions of Hurricane Characteristics [C]. Huston: Proc Offshore Technology Conference, OTC4 325, 1982: 305-313.

    [5] Liu D F, Pang L. Joint Probability Analysis of Hurricane Katrina 2005[C]. Proc Intern Offshore & Polar Eng. Conf(ISOPE), USA: San Francisco, USA: 2006, 3: 74-80.

    [6] Liu D F, Pang L, Xie B T. Typhoon Disaster in China—Prediction, Prevention and Mitigation, Natural Hazards[C]. Germany: Springer-Verlag, 2009, 49: 421-436.

    [7] Liu D F, Wang F Q, Wang Y F. Prediction of Typhoon Triggered Sea Hazards in China[C]. Bremen, Germany: Proc OCEANS’ 09, IEEE, 090204-001. 2009.

    [8] 劉德輔, 龐亮. 中國(guó)沿海臺(tái)風(fēng)災(zāi)害區(qū)劃、防臺(tái)風(fēng)標(biāo)準(zhǔn)應(yīng)急制定、防臺(tái)風(fēng)應(yīng)急評(píng)估標(biāo)準(zhǔn)制定。國(guó)家防汛抗旱總指揮部委托項(xiàng)目(20060120)[S]. 北京: 國(guó)家防辦組織專(zhuān)家委員會(huì), 2008.

    Liu D F, Pang L. “Typhoon Disaster Prediction, Prevention and Guideline of Emergency Plan for China Sea Coasts”: China State Flood Control and Drought Relief Headquarters 2008 established project (20060120)[S]. Beijing: Experts committee of National Disaster Provention Office, 2008.

    [9] 劉德輔, 龐亮, 史宏達(dá). 臺(tái)風(fēng)影響海區(qū)設(shè)計(jì)波浪概率預(yù)測(cè)理論、計(jì)算方法及規(guī)范條文[S]. “國(guó)家行業(yè)規(guī)范”海港水文“2008修訂版. 北京: 人民交通出版社, 2013.

    Liu D F, Pang L, Shi H D. Design Wave Criteria for Typhoon Influenced Sea Areas: Probability Prediction Theory, Calculation Method and Design Calibrations[S]. Published “ China Code for Sea Port Hydrology”. Beijing: Press of People’s Transportation, 2013.

    [10] Kirby W H, Moss M E. Summary of flood-frequency analysis in the United States[J]. J Hydrology, 1987, 96: 5-14.

    [11] Ochi M K. Stochastic analysis and probabilistic prediction of random seas[J]. Adv Hydro-Science, 1982, 13: 217-315.

    [12] Langley R M, El-Shaarawi A H. On the calculation of extreme wave heights: A review [J]. Ocean Eng, 1986, 13(1): 93-118.

    [13] Naffa M G, Fanos A M, Elganainy M A. Characteristics of waves off the Mediterranean Coast of Egypt [J]. J Coast Res, 1991, 7(3): 665-676.

    [14] Quek S T, Cheong H F. Prediction of extreme 3 s gusts accounting for seasonal effects[J]. Structure Safety, 1992, 11(2): 121-129.

    [15] Cheong H F. Estimation of the minimum pressure coefficient due to gusts [J]. Structure Safety, 1995, 17: 1-16.

    [16] Naffaa M G. Wave climate along the Nile delta coast[J]. Journal of Coastal Research, 1995, 11: 219-229.

    [17] Bea R. Reliability assessment and management:lessons from hurricane Katrina[C].[s.1]: Proceedings of the Offshore Mech. & Arc. Eng. 2007: OMAE2007-29650.

    [18] Chowdhury A G, Huang E, Erwin J. Novel full-scale wind- structure interaction experimentation for mitigating hurricane induced coastal disasters[J]. Far East Journal of Ocean Research, 2009, 2: 1-27.

    [19] Liu D F, Shi H D, Liu G L. Risk Assessment for Nuclear Power Plants Against Natural Disasters: Probability Prediction and Disaster Prevention[M]. USA:NOVA Science Publisher Inc. Free Charge Published, 2015.

    [20] Coles S, Simiu E. Estimating uncertainty in the extreme value analysis of data generated by a hurricane simulation model[J]. J of Engineering Mech. ASCE, 2003, 129(11): 1288-1294.

    [21] Casson E, Coles S. Simulation and extremenal analysis of hurricane events[J]. Appl Statist, 2000, 49(2): 227-245.

    [22] Georgiou P N, Davenport A G, Vickery P J. Design wind speeds in regions dominated by tropical cyclones[J]. Journal Wind Eng. Indu Aero. Dyn, 1983(13): 139-152.

    StudyonHarbourandCoastalEngineeringDisasterPrediction:TheoryandApplication

    LIU De-Fu, LIU Gui-Lin, WANG Feng-Qing, WANG Yong-Fa

    (1.College of Engineering, Ocean University of China,Qingdao 266003,China; 2.Qingdao New Qianwan Container Terminal Co.,Ltd.Qingdao 266555,China)

    This paper introduced our proposed unit-variant and multivariate compound extreme value distribution (CEVD and MCEVD) models. It is showed that all the planned, designed and constructed coastal infrastructures accepted the traditional safety regulations are menaced by possibility of future typhoon/hurricane disasters and cannot satisfy the safety requirements with the increasing tendency of the extreme natural hazards in the past 35 years. 2005 hurricane Katrina, Rita and 2012 hurricane Sandy induced disasters proved 1982 CEVD and 2006 MCEVD predicted extreme hazards in New Orleans, Gulf of Mexico and Philadelphia areas. 2013 typhoon Fitow induced disaster in China also proved MCEVD 2006 predicted results. This paper shows that all the planned, designed and constructed coastal infrastructures accepted the traditional safety regulations are menaced by possibility of future typhoon/hurricane disasters and cannot satisfy the safety requirements with the increasing tendency of the extreme natural hazards.

    Typhoon disaster; CEVD;MCEVD; coastal and port engineering; design code calibration

    U656.2

    A

    1672-5174(2018)02-091-05

    10.16441/j.cnki.hdxb.20170162

    劉德輔, 劉桂林, 王鳳清, 等. 港口及海岸工程防災(zāi)理論體系的研發(fā)及其應(yīng)用[J]. 中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2018, 48(2): 91-95.

    LIU De-Fu, LIU Gui-Lin, WANG Feng-Qing, et al. Study on harbour and coastal engineering disaster prediction:Theory and application[J]. Periodical of Ocean University of China, 2018, 48(2): 91-95.

    國(guó)家自然科學(xué)基金項(xiàng)目“基于多目標(biāo)四層嵌套概率模式的海岸工程風(fēng)險(xiǎn)評(píng)估研究”(51379195);山東省自然科學(xué)基金項(xiàng)目“海岸工程多目標(biāo)三層嵌套聯(lián)合概率預(yù)測(cè)及風(fēng)險(xiǎn)評(píng)估模式研究”(ZR2013EEM034)資助

    Supported by National Natural Science Foundation of China for the Project “Coastal Engineering and Risk Assessment Based on A Four-Layer Nested Multi-Objective Probability Model” (51379195);Natural Science Foundation of Shandong Province for the Project “Three-Layer Nested multi-Objective Probability Prediction and Risk Assessment in Coastal Engineering” (ZR2013EEM034)

    2017-04-10;

    2017-05-11

    劉德輔(1936-),男,博導(dǎo)。E-mail:liu@ ouc.edu.cn

    責(zé)任編輯 龐 旻

    猜你喜歡
    颶風(fēng)極值臺(tái)風(fēng)
    颶風(fēng)來(lái)襲
    臺(tái)風(fēng)過(guò)韓
    極值點(diǎn)帶你去“漂移”
    極值點(diǎn)偏移攔路,三法可取
    臺(tái)風(fēng)來(lái)了
    小讀者(2020年4期)2020-06-16 03:33:46
    一類(lèi)“極值點(diǎn)偏移”問(wèn)題的解法與反思
    臺(tái)風(fēng)愛(ài)搗亂
    臺(tái)風(fēng)來(lái)時(shí)怎樣應(yīng)對(duì)
    揭開(kāi)殺手真面目
    瘋狂的“殺手”
    国产一区二区在线观看日韩 | 99久久精品国产亚洲精品| 久久精品国产清高在天天线| 不卡一级毛片| 欧美一级毛片孕妇| 日韩高清综合在线| 男女之事视频高清在线观看| 国产一区二区在线av高清观看| 免费av毛片视频| 国产一区在线观看成人免费| 最好的美女福利视频网| 不卡一级毛片| 亚洲自偷自拍图片 自拍| www.精华液| 国内精品久久久久久久电影| 老鸭窝网址在线观看| 亚洲精品色激情综合| 免费看a级黄色片| 中文资源天堂在线| 精品国内亚洲2022精品成人| 色视频www国产| 免费人成视频x8x8入口观看| 深夜精品福利| 久久精品夜夜夜夜夜久久蜜豆| 亚洲真实伦在线观看| 国产激情偷乱视频一区二区| 18禁黄网站禁片午夜丰满| 欧美黑人欧美精品刺激| 两个人的视频大全免费| 少妇的逼水好多| 色视频www国产| 欧美在线一区亚洲| 1024香蕉在线观看| 久久精品91蜜桃| 99国产精品99久久久久| 精品无人区乱码1区二区| 国产精品女同一区二区软件 | 大型黄色视频在线免费观看| 精品国产超薄肉色丝袜足j| 亚洲一区二区三区色噜噜| 淫妇啪啪啪对白视频| 黄色片一级片一级黄色片| 天天添夜夜摸| 18禁黄网站禁片免费观看直播| h日本视频在线播放| 国产激情偷乱视频一区二区| 美女高潮的动态| 欧美成人一区二区免费高清观看 | 男女那种视频在线观看| 亚洲第一电影网av| 精品久久久久久成人av| 少妇的丰满在线观看| 在线看三级毛片| 国产精品美女特级片免费视频播放器 | 午夜福利18| 男人舔女人的私密视频| 午夜福利高清视频| 日本黄色视频三级网站网址| 老司机深夜福利视频在线观看| 人妻丰满熟妇av一区二区三区| 免费看日本二区| 国产淫片久久久久久久久 | 色精品久久人妻99蜜桃| 九九久久精品国产亚洲av麻豆 | 亚洲精品一区av在线观看| 白带黄色成豆腐渣| 亚洲av成人精品一区久久| 美女 人体艺术 gogo| 日韩三级视频一区二区三区| 午夜影院日韩av| 久久国产精品影院| 天堂网av新在线| 国产69精品久久久久777片 | 日韩人妻高清精品专区| 午夜影院日韩av| 国产精品1区2区在线观看.| 国产精品日韩av在线免费观看| 成年女人永久免费观看视频| 色播亚洲综合网| 国产精华一区二区三区| 狂野欧美白嫩少妇大欣赏| 天堂影院成人在线观看| 免费在线观看视频国产中文字幕亚洲| 国产精品一区二区精品视频观看| 国产视频一区二区在线看| 蜜桃久久精品国产亚洲av| 亚洲欧美日韩无卡精品| 男人舔女人的私密视频| 日韩精品中文字幕看吧| 小蜜桃在线观看免费完整版高清| 国产伦精品一区二区三区四那| 日本三级黄在线观看| 99久久成人亚洲精品观看| 亚洲第一电影网av| 久久热在线av| 亚洲欧美日韩高清在线视频| 中文字幕av在线有码专区| 99久久99久久久精品蜜桃| 欧美乱色亚洲激情| 亚洲欧美日韩高清在线视频| 国产高清三级在线| 特大巨黑吊av在线直播| 欧美极品一区二区三区四区| 色综合欧美亚洲国产小说| 欧美黄色淫秽网站| 国产美女午夜福利| 美女cb高潮喷水在线观看 | 久久人妻av系列| 天天添夜夜摸| АⅤ资源中文在线天堂| 国产不卡一卡二| 男女那种视频在线观看| 黑人欧美特级aaaaaa片| 99热这里只有精品一区 | 免费观看的影片在线观看| 18禁裸乳无遮挡免费网站照片| 午夜精品在线福利| 精品欧美国产一区二区三| 久久久精品欧美日韩精品| 香蕉av资源在线| 日韩成人在线观看一区二区三区| 国内精品美女久久久久久| 好男人电影高清在线观看| 成人三级黄色视频| 999精品在线视频| 手机成人av网站| 韩国av一区二区三区四区| 久久久精品大字幕| 99热这里只有是精品50| 黄色丝袜av网址大全| 真实男女啪啪啪动态图| 99热这里只有是精品50| 看黄色毛片网站| 老司机深夜福利视频在线观看| 一个人看的www免费观看视频| 在线观看66精品国产| 国产高潮美女av| 色在线成人网| 国产一区二区在线观看日韩 | 18美女黄网站色大片免费观看| 禁无遮挡网站| 国产欧美日韩一区二区精品| xxxwww97欧美| 美女扒开内裤让男人捅视频| 亚洲人成伊人成综合网2020| 老司机午夜十八禁免费视频| 狂野欧美激情性xxxx| 麻豆国产av国片精品| 麻豆成人午夜福利视频| 国产97色在线日韩免费| 国产高清videossex| 一级毛片精品| 99久久综合精品五月天人人| 嫁个100分男人电影在线观看| av视频在线观看入口| 国产精品女同一区二区软件 | 午夜免费成人在线视频| 国产精品综合久久久久久久免费| 三级男女做爰猛烈吃奶摸视频| 亚洲欧美精品综合一区二区三区| 欧美日本视频| 亚洲国产欧美网| 久久精品91蜜桃| 日韩欧美精品v在线| 久久中文看片网| 在线观看美女被高潮喷水网站 | 国内久久婷婷六月综合欲色啪| 可以在线观看毛片的网站| 超碰成人久久| 又紧又爽又黄一区二区| 怎么达到女性高潮| 高潮久久久久久久久久久不卡| 亚洲午夜理论影院| 精品国产超薄肉色丝袜足j| 国产男靠女视频免费网站| 九九久久精品国产亚洲av麻豆 | 国产精品日韩av在线免费观看| 亚洲无线在线观看| АⅤ资源中文在线天堂| 国产成人福利小说| 91在线精品国自产拍蜜月 | 成人午夜高清在线视频| 91av网站免费观看| 精品久久久久久久毛片微露脸| 亚洲熟妇熟女久久| 深夜精品福利| 岛国视频午夜一区免费看| 黄片小视频在线播放| 丰满的人妻完整版| 大型黄色视频在线免费观看| 无限看片的www在线观看| 好看av亚洲va欧美ⅴa在| 亚洲av片天天在线观看| 国产精品99久久久久久久久| 一区二区三区高清视频在线| 亚洲专区中文字幕在线| 91在线观看av| 成人特级黄色片久久久久久久| 欧美激情久久久久久爽电影| 久久精品人妻少妇| 国产免费av片在线观看野外av| 又粗又爽又猛毛片免费看| 亚洲av成人一区二区三| 亚洲中文日韩欧美视频| 天天一区二区日本电影三级| 亚洲成人久久性| 搡老岳熟女国产| 国产高清视频在线观看网站| 偷拍熟女少妇极品色| 90打野战视频偷拍视频| 亚洲精品粉嫩美女一区| 亚洲精品粉嫩美女一区| 久久香蕉国产精品| 国产精品九九99| 久久久久国产精品人妻aⅴ院| 中文字幕av在线有码专区| netflix在线观看网站| 在线观看舔阴道视频| 91九色精品人成在线观看| 国产精品亚洲美女久久久| 两性午夜刺激爽爽歪歪视频在线观看| 特级一级黄色大片| 亚洲午夜精品一区,二区,三区| 亚洲色图av天堂| 免费av不卡在线播放| 成人特级黄色片久久久久久久| 欧美日韩黄片免| 免费在线观看影片大全网站| 欧美日韩黄片免| 91字幕亚洲| 日韩欧美精品v在线| 男女下面进入的视频免费午夜| 久久精品国产亚洲av香蕉五月| 黄色女人牲交| 精品国产乱码久久久久久男人| 老鸭窝网址在线观看| 桃红色精品国产亚洲av| 精品久久久久久成人av| 少妇的逼水好多| 三级男女做爰猛烈吃奶摸视频| 婷婷精品国产亚洲av在线| 丰满人妻一区二区三区视频av | 脱女人内裤的视频| 一个人看的www免费观看视频| 1000部很黄的大片| 国产精品久久久久久亚洲av鲁大| 日韩成人在线观看一区二区三区| 国产一区二区在线av高清观看| 久久久久九九精品影院| 特大巨黑吊av在线直播| 99久久99久久久精品蜜桃| 国内精品一区二区在线观看| 国内少妇人妻偷人精品xxx网站 | 久久精品国产99精品国产亚洲性色| 日韩欧美在线二视频| 中文字幕高清在线视频| 精品国产亚洲在线| 欧美在线黄色| 亚洲中文av在线| 亚洲成av人片免费观看| 欧美一级毛片孕妇| 一卡2卡三卡四卡精品乱码亚洲| 最新美女视频免费是黄的| 99久久精品一区二区三区| 搡老熟女国产l中国老女人| 成人国产综合亚洲| 在线观看免费视频日本深夜| 免费在线观看视频国产中文字幕亚洲| 国产精品 欧美亚洲| 一进一出抽搐动态| 日韩av在线大香蕉| 成年女人永久免费观看视频| 在线观看一区二区三区| 午夜视频精品福利| 国产精品国产高清国产av| 亚洲 欧美一区二区三区| 欧美成人免费av一区二区三区| 精品一区二区三区视频在线 | 精品国产亚洲在线| 亚洲最大成人中文| 在线国产一区二区在线| www国产在线视频色| svipshipincom国产片| 国内精品一区二区在线观看| 亚洲 国产 在线| 男女床上黄色一级片免费看| 国产成人精品久久二区二区91| 在线看三级毛片| 小蜜桃在线观看免费完整版高清| 12—13女人毛片做爰片一| 日韩三级视频一区二区三区| 国产欧美日韩一区二区精品| 精品国产乱子伦一区二区三区| 国产精品99久久99久久久不卡| 99久久久亚洲精品蜜臀av| 亚洲国产看品久久| 亚洲av第一区精品v没综合| 午夜激情欧美在线| 久久久久久久久中文| 欧美日韩黄片免| 久久精品国产清高在天天线| 欧美性猛交黑人性爽| 黄片大片在线免费观看| 高清在线国产一区| 9191精品国产免费久久| 日本免费一区二区三区高清不卡| 久久久国产精品麻豆| 俺也久久电影网| 国产精品永久免费网站| 不卡av一区二区三区| 亚洲真实伦在线观看| 色综合亚洲欧美另类图片| 成在线人永久免费视频| 久久性视频一级片| 99久久综合精品五月天人人| 成人一区二区视频在线观看| 国产黄片美女视频| 亚洲真实伦在线观看| 五月玫瑰六月丁香| 91久久精品国产一区二区成人 | 看片在线看免费视频| 一个人看视频在线观看www免费 | 亚洲色图av天堂| 成人三级黄色视频| 国产伦人伦偷精品视频| 老司机午夜福利在线观看视频| 超碰成人久久| 国产主播在线观看一区二区| 男女视频在线观看网站免费| 高清在线国产一区| 夜夜看夜夜爽夜夜摸| 深夜精品福利| 99热6这里只有精品| 一级毛片女人18水好多| 不卡av一区二区三区| 国产精品久久久久久亚洲av鲁大| 哪里可以看免费的av片| 熟女人妻精品中文字幕| 午夜日韩欧美国产| 女人高潮潮喷娇喘18禁视频| 免费大片18禁| 亚洲欧美激情综合另类| 啦啦啦观看免费观看视频高清| 午夜福利在线观看免费完整高清在 | 国产亚洲欧美98| 亚洲无线观看免费| 国产成人系列免费观看| 午夜福利在线观看吧| 级片在线观看| 欧美精品啪啪一区二区三区| av中文乱码字幕在线| 日韩三级视频一区二区三区| 亚洲五月婷婷丁香| 日韩欧美在线二视频| 又黄又粗又硬又大视频| 18禁观看日本| 国产高清视频在线观看网站| 丝袜人妻中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 久久精品亚洲精品国产色婷小说| 久久精品人妻少妇| 国产成人福利小说| 精品欧美国产一区二区三| 黄色丝袜av网址大全| 白带黄色成豆腐渣| 成人三级黄色视频| 亚洲性夜色夜夜综合| 午夜激情福利司机影院| 日本一本二区三区精品| 天天一区二区日本电影三级| 日本与韩国留学比较| 国产高清有码在线观看视频| 欧美中文综合在线视频| 亚洲 欧美 日韩 在线 免费| 床上黄色一级片| aaaaa片日本免费| 超碰成人久久| 免费av毛片视频| 麻豆国产av国片精品| 欧美乱码精品一区二区三区| 国产精品影院久久| 91av网站免费观看| 成熟少妇高潮喷水视频| 亚洲 欧美 日韩 在线 免费| 两性午夜刺激爽爽歪歪视频在线观看| 黄色 视频免费看| 舔av片在线| 蜜桃久久精品国产亚洲av| 一级毛片女人18水好多| 男女午夜视频在线观看| 欧美日韩乱码在线| 99久久综合精品五月天人人| 可以在线观看的亚洲视频| 欧美一级a爱片免费观看看| 亚洲av成人不卡在线观看播放网| 制服丝袜大香蕉在线| 丰满人妻熟妇乱又伦精品不卡| 我的老师免费观看完整版| 舔av片在线| 蜜桃久久精品国产亚洲av| 日韩中文字幕欧美一区二区| 国产91精品成人一区二区三区| 欧美色视频一区免费| 国产亚洲精品一区二区www| 国产蜜桃级精品一区二区三区| 欧美成人一区二区免费高清观看 | 搡老熟女国产l中国老女人| 一a级毛片在线观看| 久久午夜亚洲精品久久| netflix在线观看网站| 国产精品香港三级国产av潘金莲| 黄片大片在线免费观看| 在线国产一区二区在线| 最新中文字幕久久久久 | 亚洲人成伊人成综合网2020| 国产亚洲精品av在线| 村上凉子中文字幕在线| 91在线精品国自产拍蜜月 | 欧美日韩一级在线毛片| 麻豆一二三区av精品| 亚洲精品在线观看二区| 在线观看美女被高潮喷水网站 | 国产激情欧美一区二区| 91字幕亚洲| 免费看美女性在线毛片视频| 亚洲欧美精品综合久久99| 一个人看的www免费观看视频| 国产免费av片在线观看野外av| 十八禁人妻一区二区| 亚洲人成网站高清观看| 国产精品综合久久久久久久免费| a级毛片a级免费在线| 欧美丝袜亚洲另类 | 脱女人内裤的视频| 成人国产一区最新在线观看| 99国产综合亚洲精品| 少妇的丰满在线观看| 网址你懂的国产日韩在线| 日韩三级视频一区二区三区| 国产精品亚洲av一区麻豆| 久久久久精品国产欧美久久久| av在线蜜桃| 91在线精品国自产拍蜜月 | 国内久久婷婷六月综合欲色啪| 亚洲专区中文字幕在线| 久久精品91无色码中文字幕| 少妇熟女aⅴ在线视频| 日韩精品中文字幕看吧| 免费高清视频大片| 精品久久蜜臀av无| 国产高清有码在线观看视频| 最近最新中文字幕大全免费视频| 97碰自拍视频| 成人特级av手机在线观看| 桃色一区二区三区在线观看| 少妇熟女aⅴ在线视频| 美女免费视频网站| 又大又爽又粗| 国产精品电影一区二区三区| 日本与韩国留学比较| 国产黄色小视频在线观看| 国产av一区在线观看免费| 午夜福利在线在线| 色播亚洲综合网| 久久热在线av| 日本免费一区二区三区高清不卡| 69av精品久久久久久| 高潮久久久久久久久久久不卡| 国产一区二区在线观看日韩 | 国产成人av激情在线播放| 国产男靠女视频免费网站| 两性夫妻黄色片| 老司机在亚洲福利影院| 亚洲av日韩精品久久久久久密| 国产综合懂色| 国产不卡一卡二| 国产视频内射| 国语自产精品视频在线第100页| 九色成人免费人妻av| 精品午夜福利视频在线观看一区| 国产一级毛片七仙女欲春2| 两个人视频免费观看高清| 国产成人精品久久二区二区免费| 国产私拍福利视频在线观看| 特大巨黑吊av在线直播| 国产av在哪里看| 久久午夜亚洲精品久久| 国产免费男女视频| 亚洲18禁久久av| 亚洲av电影不卡..在线观看| 男人舔奶头视频| 男插女下体视频免费在线播放| 亚洲av片天天在线观看| 男女做爰动态图高潮gif福利片| 免费在线观看成人毛片| 全区人妻精品视频| 欧美激情在线99| 99久久99久久久精品蜜桃| 国产成人福利小说| 成人18禁在线播放| 成人国产一区最新在线观看| 午夜福利在线观看吧| 在线国产一区二区在线| 九九久久精品国产亚洲av麻豆 | 亚洲av电影不卡..在线观看| 老熟妇仑乱视频hdxx| 两个人视频免费观看高清| 啦啦啦韩国在线观看视频| 最新中文字幕久久久久 | svipshipincom国产片| 亚洲 国产 在线| 中亚洲国语对白在线视频| 国产伦在线观看视频一区| 97超视频在线观看视频| 性欧美人与动物交配| 悠悠久久av| 亚洲精华国产精华精| 两性午夜刺激爽爽歪歪视频在线观看| 成人特级av手机在线观看| 亚洲在线自拍视频| 久久精品国产亚洲av香蕉五月| 久久草成人影院| 神马国产精品三级电影在线观看| 青草久久国产| 亚洲精品一卡2卡三卡4卡5卡| 男女之事视频高清在线观看| 麻豆一二三区av精品| 久久中文看片网| 久久草成人影院| 国产精品久久久久久亚洲av鲁大| 一本一本综合久久| 变态另类成人亚洲欧美熟女| 最近最新中文字幕大全免费视频| 在线国产一区二区在线| 国产av不卡久久| 淫妇啪啪啪对白视频| 热99在线观看视频| 久久精品国产99精品国产亚洲性色| 三级国产精品欧美在线观看 | 黑人操中国人逼视频| 亚洲精品粉嫩美女一区| 97超级碰碰碰精品色视频在线观看| 欧美不卡视频在线免费观看| 免费大片18禁| 成人精品一区二区免费| 一二三四社区在线视频社区8| 男女那种视频在线观看| 久久久成人免费电影| 日本在线视频免费播放| 国产一区二区在线观看日韩 | 麻豆国产97在线/欧美| 国产 一区 欧美 日韩| 精华霜和精华液先用哪个| 精品一区二区三区四区五区乱码| 午夜福利在线观看免费完整高清在 | 午夜福利成人在线免费观看| 黄色 视频免费看| 美女cb高潮喷水在线观看 | or卡值多少钱| 国产激情欧美一区二区| 深夜精品福利| 亚洲欧洲精品一区二区精品久久久| 精华霜和精华液先用哪个| 欧美色欧美亚洲另类二区| 日韩欧美 国产精品| 国产亚洲精品久久久com| 亚洲最大成人中文| 香蕉丝袜av| 久久精品国产综合久久久| 亚洲av日韩精品久久久久久密| 午夜精品一区二区三区免费看| 午夜免费观看网址| www.精华液| 国产一区在线观看成人免费| 久久久久性生活片| 日韩欧美 国产精品| 久久久色成人| 每晚都被弄得嗷嗷叫到高潮| 1024手机看黄色片| 国产高清视频在线播放一区| 一本久久中文字幕| 在线a可以看的网站| 这个男人来自地球电影免费观看| 九九在线视频观看精品| 亚洲精品久久国产高清桃花| 黑人巨大精品欧美一区二区mp4| 中国美女看黄片| 熟女电影av网| a在线观看视频网站| 精品国产超薄肉色丝袜足j| 久久性视频一级片| 国产一级毛片七仙女欲春2| 亚洲成av人片免费观看| 国产99白浆流出| 1024手机看黄色片| 伊人久久大香线蕉亚洲五| 非洲黑人性xxxx精品又粗又长| 亚洲 国产 在线| 亚洲人成伊人成综合网2020| 日韩欧美国产一区二区入口| 1024手机看黄色片| 亚洲 国产 在线| 桃色一区二区三区在线观看| 色在线成人网| 成人三级做爰电影| 男女做爰动态图高潮gif福利片| 国产高清videossex| 亚洲av免费在线观看| 淫妇啪啪啪对白视频| 五月玫瑰六月丁香| 中文字幕久久专区| 亚洲性夜色夜夜综合| 国产视频一区二区在线看| 亚洲成av人片免费观看| 12—13女人毛片做爰片一|