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

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

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

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

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

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

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

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

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

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

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

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

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

當(dāng)臺風(fēng)頻次為泊松分布，波高符合Gumbel分布時，可簡化為：

HP=u+XP/α。

(1)

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

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

(2)

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

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

1.3 2005 、2012美國兩場颶風(fēng)重災(zāi)及2013中國菲特臺風(fēng)災(zāi)害，驗證了CEVD和MCEVD 1982和2006預(yù)測結(jié)果的正確性

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

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

表1 墨西哥灣A區(qū)卡特里娜颶風(fēng)強度概率分析(按文獻[19]圖8計算)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)的強度和特征，驗證了使用復(fù)合極值分布理論1982年預(yù)測結(jié)果的正確性(見表1)。圖1則為作者使用2002年以來新開發(fā)的MCEVD分析颶風(fēng)及其誘發(fā)的災(zāi)害性海況，顯示2005年以前美國不同學(xué)者公開發(fā)表對上述海域預(yù)測成果，明顯低于MCEVD的預(yù)測值(見文獻[9]中圖6)。MCEVD的預(yù)測結(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)速相近 (見表2)。

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

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

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

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

表2 不同方法預(yù)測結(jié)果比較(見文獻[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 美國大西洋沿岸各區(qū)不同重現(xiàn)期風(fēng)暴增水1982年預(yù)測結(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

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

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

Note：①Model；②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)帶來巨大災(zāi)害的教訓(xùn)，對海岸，港口工程的防災(zāi)設(shè)計標(biāo)準(zhǔn)的確定，特別是濱海核電站防護工程的安全非常重要，全盤翻譯引用國外規(guī)范的隱患，必須引以為鑒[19]。

(2)鑒于沿海港區(qū)大風(fēng)導(dǎo)致港內(nèi)集裝箱垜倒塌事故時有發(fā)生，例如：2014年8月鹽田港幾十個集裝箱遭臺風(fēng)襲擊，很多民房破壞，造成500多萬損失：2010年8月蕪湖港臺風(fēng)導(dǎo)致兩座門吊倒塌，造成2死1傷：2005年Octavia輪?？繉幉ū眮龈?，溫帶氣旋襲擊，波高7～8 m，造成集裝箱損失。建議針對所在海域臺風(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] 劉德輔， 龐亮. 中國沿海臺風(fēng)災(zāi)害區(qū)劃、防臺風(fēng)標(biāo)準(zhǔn)應(yīng)急制定、防臺風(fēng)應(yīng)急評估標(biāo)準(zhǔn)制定。國家防汛抗旱總指揮部委托項目(20060120)[S]. 北京： 國家防辦組織專家委員會， 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] 劉德輔， 龐亮， 史宏達. 臺風(fēng)影響海區(qū)設(shè)計波浪概率預(yù)測理論、計算方法及規(guī)范條文[S]. “國家行業(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]. 中國海洋大學(xué)學(xué)報(自然科學(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.

國家自然科學(xué)基金項目“基于多目標(biāo)四層嵌套概率模式的海岸工程風(fēng)險評估研究”(51379195)；山東省自然科學(xué)基金項目“海岸工程多目標(biāo)三層嵌套聯(lián)合概率預(yù)測及風(fē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é)任編輯 龐 旻