劉德輔，史宏達，劉桂林，王風(fēng)清

(中國海洋大學(xué)，山東 青島 266003)

我國濱海核電站防護工程設(shè)計標(biāo)準(zhǔn)研究

劉德輔，史宏達，劉桂林，王風(fēng)清

(中國海洋大學(xué)，山東 青島 266003)

面對全球極端氣候頻發(fā)的趨勢，國內(nèi)外核電站海岸防護工程普遍采用“可能最大臺風(fēng)暴潮”、“設(shè)計基準(zhǔn)洪水”等具有含混性的規(guī)定作為確定防災(zāi)設(shè)計標(biāo)準(zhǔn)的依據(jù)，其內(nèi)涵的各種不確定性因素必然會對防災(zāi)決策帶來重大的致災(zāi)隱患。按照國務(wù)院核電站安全規(guī)劃(2011-2020)提出的“必須按照全球最高安全要求”，全面開展濱海核電站防護工程設(shè)計標(biāo)準(zhǔn)研究己勢在必行。使用多維復(fù)合極值分布理論及雙層嵌套多目標(biāo)聯(lián)合概率模式，將核電站水文與氣象的研究標(biāo)準(zhǔn)合二為一，對IAEA及我國濱海核電站防護工程設(shè)計標(biāo)準(zhǔn)進行了風(fēng)險評估，供有關(guān)部門參考。

濱海核電站；防護工程；設(shè)計標(biāo)準(zhǔn)；設(shè)計基準(zhǔn)洪水；多維復(fù)合極值分布；風(fēng)險評估

Abstract: With the tendency of increasing extreme natural hazards, the typhoon, hurricane, and tropical cyclone induced surges, waves and rainstorms menace Nuclear Power Plants (NPP) in coastal areas. National Nuclear Safety Administration of China and International Atomic Energy Agency (IAEA) have recommended safety regulations for NPP site evaluation installation and coastal defence infrastructures. These standards include Probable Maximum Hurricane /Typhoon (PMH/T), Probable Maximum Storm Surge (PMSS), Probable Maximum Flood (PMF) as well as Design Basis Flood (DBF). According to the 2011-2020 safety planning of China State Council for nuclear power plants, it is necessary to do a comprehensive research on design standards for protective engineering and structural technology of the NPP based on the world’s highest safety requirements. This paper discusses the joint probability analysis of meteorological and oceanographic factors based on our proposed Compound Extreme Value Distribution (CEVD), Multivariate Compound Extreme Value Distribution (MCEVD), and Double Layer Nested Multi-Objective Probability Model (DLNMOPM), and ours are compared with IAEA recommended safety regulation design criteria for NPP coastal defence infrastructures in China.

Keywords: nuclear power plants; protection engineering; design standard; design basic floods; multivariate compound extreme value distribution; risk assessment

隨著全球氣候變暖和海平面上升等因素加劇的趨勢，臺風(fēng)、巨浪、風(fēng)暴潮、暴雨和洪災(zāi)等已成為影響濱海核電站防護工程安全的重大問題。日本核電災(zāi)害發(fā)生后，許多學(xué)者對地震-海嘯的研究表明：當(dāng)Manila和Rykyu海溝發(fā)生9級地震時，對我國東南沿海誘發(fā)的海嘯巨浪不超過6 m[1]。世界氣象組織指出：百分之九十的自然災(zāi)害為極端氣象災(zāi)害，例如臺風(fēng)、颶風(fēng)和熱帶風(fēng)暴災(zāi)害。臺風(fēng)是威脅我國濱海核電站的最重要致災(zāi)因素。研究臺風(fēng)、巨浪、風(fēng)暴增水、暴雨、洪水及其與天文大潮組合的小概率極端事件發(fā)生的可能性非常必要。2006年桑美臺風(fēng)誘發(fā)的7 m巨浪和3.8 m風(fēng)暴增水造成沙埕港千艘船只毀壞和多人死亡。如果臺風(fēng)延后兩小時登陸，則恰逢天文大潮，這種“三碰頭”的組合，可完全淹沒包括核電站在內(nèi)的福建和浙江大部分土地，其災(zāi)難性后果難以估量。雖然這屬于小概率極端事件，但是這類事件的發(fā)生是完全可能的。因此，按照國務(wù)院核電站安全規(guī)劃(2011-2020)提出的“必須按照全球最高安全要求”，全面開展濱海核電站防護工程設(shè)計標(biāo)準(zhǔn)研究己勢在必行。必須重新檢驗核電站防護工程設(shè)計標(biāo)準(zhǔn)應(yīng)對各種極端自然災(zāi)害的可靠性。目前我國濱海核電站設(shè)計基準(zhǔn)洪水位的計算分析以確定論方法為主，以“10%超過概率天文高潮位+可能最大風(fēng)暴潮增水及百年一遇波高”的組合作為防護工程設(shè)計標(biāo)準(zhǔn)。但是，在計算防護工程標(biāo)準(zhǔn)時，我國核電站有關(guān)規(guī)范的“可能最大風(fēng)暴潮增水”、“設(shè)計基準(zhǔn)洪水”(DBF)、IAEA2006-2009和U.S.NRC 2011在新的安全規(guī)程中提出的“風(fēng)暴潮、天文潮和巨浪的組合”[2-8]內(nèi)涵各種不確定性因素，這些不確定性因素必然會對防災(zāi)措施的決策帶來重大的致災(zāi)隱患。因此，在港口、海岸、近海及河口城市防護工程的設(shè)計中，由于概率統(tǒng)計理論和方法上的限制，無法考慮惡劣環(huán)境條件下多種荷載共同作用的問題必須得到妥善解決。

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

國內(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)過嚴格的數(shù)學(xué)推導(dǎo)，可構(gòu)成一種新的極值分布模式——復(fù)合極值分布(compound extreme value distribution，簡稱CEVD)。根據(jù)統(tǒng)計檢驗結(jié)果，臺風(fēng)(颶風(fēng))影響各海區(qū)的頻次符合Poisson分布，臺風(fēng)(颶風(fēng))波浪的不同特征，可分別采用Gumbel和Weibull分布，則可構(gòu)成Poisson-Gumbel CEVD分布(用于我國臺風(fēng)影響海域)[9]和Poisson-Weibull CEVD分布(用于美國墨西哥灣和大西洋沿岸)[10]。

近年來，CEVD完成了從一維拓寬到多維復(fù)合極值分布(multivariate compound extreme value distribution，簡稱MCEVD)的理論推導(dǎo)和相應(yīng)工程應(yīng)用。2005年美國Katrina颶風(fēng)對新奧爾良的毀滅性破壞，進一步驗證了CEVD模式應(yīng)用于颶風(fēng)概率預(yù)測結(jié)果的正確性[11-17]。2006年有關(guān)文獻[18-19]在重要國際會議和刊物公開發(fā)表后，得到了國際學(xué)術(shù)界的高度評價。2008年交通部在國家行業(yè)標(biāo)準(zhǔn)“海港水文規(guī)范”(JTJ213-98)修訂稿中，正式采用了“復(fù)合極值分布理論”作為設(shè)計波高概率預(yù)測的理論和方法，首次取代了該國家行業(yè)標(biāo)準(zhǔn)中持續(xù)使用近四十年的P-Ⅲ型曲線。2008年美國佛羅里達國際颶風(fēng)研究中心在海岸防護颶風(fēng)災(zāi)害工程設(shè)計中，引用文獻[9-10,19-20]作為防護颶風(fēng)災(zāi)害設(shè)計標(biāo)準(zhǔn)的依據(jù)[21]。多年來MCEVD用于城市防災(zāi)、三峽防洪、防災(zāi)區(qū)劃、奧帆工程、海洋平臺和濱海核電站防護工程設(shè)計標(biāo)準(zhǔn)的風(fēng)險分析以及暴雨，洪水誘發(fā)滑坡、泥石流等致災(zāi)因素概率預(yù)測研究成果[22-34]；有關(guān)內(nèi)容被美國NOVA科學(xué)出版社收錄在《Natural Disasters: Prevention, Risk Factors and Mitigation》一書的第一章“Typhoon/Hurricane/Tropical Cyclone Disasters: Prediction, Prevention and Mitigation”中[35]；2015年美國NOVA科學(xué)出版社再次將有關(guān)內(nèi)容收錄在學(xué)術(shù)專著《Risk Assessment for Nuclear Power Plants against Natural Disasters: Probability Prediction and Disaster Prevention Infrastructures》中[36]。

2 預(yù)測結(jié)果比較

2005年美國發(fā)生的Katrina颶風(fēng)重災(zāi)，導(dǎo)致新奧爾良市毀滅性破壞，1 800人死亡，800億美元經(jīng)濟損失。這次颶風(fēng)的強度和特征，符合1982年使用復(fù)合極值分布理論模式預(yù)測結(jié)果[10](表1)，顯示2005災(zāi)后使用MCEVD復(fù)核結(jié)果的合理性[20,37]；同時也表明美國國家海洋和大氣管理局(NOAA)制訂的上述海區(qū)防災(zāi)標(biāo)準(zhǔn)的可能最大颶風(fēng)(PMH)和標(biāo)準(zhǔn)設(shè)計颶風(fēng)(SPH)過低，相關(guān)文獻[38-39]也得到相同的結(jié)果。

2012年10月30-31日桑迪颶風(fēng)造成堤壩沖毀，特拉華州、新澤西州及紐約周邊大面積淹沒，大量人員傷亡，數(shù)百億美元經(jīng)濟損失，再次表明傳統(tǒng)的設(shè)計標(biāo)準(zhǔn)和工程結(jié)構(gòu)標(biāo)準(zhǔn)過低，我們對上述海區(qū)颶風(fēng)強度預(yù)測成果再次得到證實[10]。

使用2002年以來新開發(fā)的 MCEVD分析颶風(fēng)及其誘發(fā)的極端風(fēng)速的預(yù)測值明顯高于2005年以前美國不同學(xué)者公開發(fā)表對上述海域預(yù)測成果(圖1)。同時，MCEVD的預(yù)測結(jié)果還顯示，重現(xiàn)期100年一遇的風(fēng)速在3區(qū)(新奧爾良區(qū))的預(yù)測值與2005卡特里娜颶風(fēng)狀況相近，明顯高于設(shè)防標(biāo)準(zhǔn)；8-9區(qū)(新澤西，紐約沿岸)的預(yù)測值則與2012年桑迪颶風(fēng)最大風(fēng)速相近(表2)。

表1 墨西哥灣A區(qū)卡特里娜颶風(fēng)強度概率分析Tab. 1 Probability analysis for hurricane Katrina in Gulf of Mexico area A

表2 不同方法預(yù)測結(jié)果比較Tab. 2 Predicted results by different models

圖1 100年一遇的極端風(fēng)速概率預(yù)測Fig. 1 Predicted 100 a return value of extreme wind velocities

圖2 美國大西洋沿岸各區(qū)不同重現(xiàn)期風(fēng)暴潮增水預(yù)測結(jié)果Fig. 2 1982 predicted storm surges with different return periods in US Atlantic areas

如圖2所示，1982年預(yù)測的費城周邊(圖中虛線)100年一遇風(fēng)暴潮增水為3.048 m，與2012年10月30日8時06分桑迪颶風(fēng)造成的費城周邊3.237 m的實測水位相當(dāng)，而美國國家海洋和大氣管理局(NOAA)對該區(qū)的預(yù)測值僅為2.292 m。

2013年臺風(fēng)菲特帶來重大經(jīng)濟損失，浙江、上海多處出現(xiàn)超警戒水位洪水，多處防護工程倒塌。以臺風(fēng)菲特對上海影響為例(如表3所示)，考慮臺風(fēng)暴潮、暴雨洪水與天文大潮的組合，2006年使用MCEVD計算的100年一遇預(yù)測結(jié)果，已大大超過上海防護工程500年一遇超警戒水位設(shè)防標(biāo)準(zhǔn)。2013年10月8日，受天文大潮頂托、上游洪水下泄、本地大暴雨和臺風(fēng)暴潮的共同影響，黃浦江沿線水位迅速上漲，下午2時05分，吳淞口實測水位達5.15 m，超過警戒水位0.35 m，僅達到MCEVD計算的60年一遇預(yù)測值[36-37，40-42]。

表3 考慮臺風(fēng)暴潮，暴雨洪水與天文大潮使用MCEVD和傳統(tǒng)外延法預(yù)測結(jié)果比較Tab. 3 Comparison between MCEVD predicted combined storm surge, flood and spring tide with traditional method

注：*代表傳統(tǒng)規(guī)范外延法預(yù)測結(jié)果。

3 使用CEVD進行濱海核電站防洪裕度分析

為了彌補用年極值系列按照某種概率模式外延推求不同重現(xiàn)期設(shè)計值作法存在的遺漏每年各次臺風(fēng)波浪提供的重要信息的弊端，CEVD模式將臺風(fēng)每年影響不同海區(qū)的頻次作為一組離散型隨機變量引入模型，增加了可利用的資料信息；同時研究P-III、Gumbel及其與CEVD 對超低概率極值的擬合能力以及置信度水平的對比，以保證濱海核電站防洪安全裕度評估的正確性。

圖3 三種分布的擬合平均差v對比圖Fig. 3 Comparison between three models calculated v

為驗證CEVD理論預(yù)測結(jié)果的正確性和可靠性，我們使用濱海核電站周邊長期實測資料，對Gumbel，P-Ⅲ和CEVD分布的擬合最大偏差Dn，平均差v，標(biāo)準(zhǔn)差d，不同重現(xiàn)期概率預(yù)測值置信區(qū)間ΔH進行了比較。將中國科學(xué)院南海海洋研究所提供的大鵬灣1951-2002實測(缺測部分為后報)近岸風(fēng)、浪資料，分別按各20年為一組進行滑動取樣，對Gumbel，P-Ⅲ和CEVD分布的擬合優(yōu)度，按擬合最大偏差Dn，擬合平均差v和擬合標(biāo)準(zhǔn)差d進行比較，比較結(jié)果如圖3及表4所示。有關(guān)統(tǒng)計特征及擬合優(yōu)度說明，詳見文獻[35,43]。

表4 三種分布擬合優(yōu)度參數(shù)及其變化范圍Tab. 4 Relative differences of predicted returns value between three models

由表4和圖3所示，無論是擬合最大偏差Dn，還是擬合平均差v和擬合標(biāo)準(zhǔn)差d，CEVD的擬合偏差數(shù)均小于P-III分布和Gumbel分布。

表5 大鵬灣不同重現(xiàn)期概率預(yù)測值及置信區(qū)間ΔH對比計算結(jié)果Tab. 5 Comparison between predicted return values with confidence intervals by three models

表5所示P-III、Gumbel以及其與CEVD對超低概率極值的擬合能力以及置信度水平的對比表明：CEVD是迄今為止可保證濱海核電站防洪安全裕度評估正確性的模式。

4 濱海核電站防護工程設(shè)計標(biāo)準(zhǔn)地構(gòu)建

目前國際上多要素聯(lián)合概率理論研究(如Copula理論等)都建立在年極值取樣或過閾取樣(POT)的基礎(chǔ)上，但年極值取樣遺漏了每年臺風(fēng)誘發(fā)的暴潮、波浪等要素的大量信息；POT取樣聯(lián)合概率則無理論依據(jù)[25,44-45]。多維復(fù)合極值分布理論(MCEVD)將極值過程取樣與聯(lián)合概率結(jié)合起來，增加了資料信息，縮小了概率預(yù)測的誤差范圍。既可滿足特大值的概率預(yù)測，又可進行多種致災(zāi)因素聯(lián)合出現(xiàn)的概率預(yù)測。美國學(xué)者在《美國洪水概率分析總結(jié)》中引用了復(fù)合極值分布理論并指出“把事件取樣(過程取樣)和聯(lián)合概率結(jié)合起來，將大大推動特大洪水概率預(yù)測的發(fā)展”[11]。MCEVD正是迄今為止能滿足這種期待的唯一理論模式，現(xiàn)已在多種工程應(yīng)用中得到證實。MCEVD模式簡介如下：

每年的臺風(fēng)頻次λ符合泊松分布：

MCEVD可表示為：

相應(yīng)的概率密度函數(shù)為：

為解決MCEVD的理論解法，引入三元嵌套Logistic模型[46]，其一般表達式為

由此可得應(yīng)用于實際計算的“泊松嵌套邏輯三維復(fù)合極值分布”(poisson nested logistic tri-variant compound extreme value distribution，簡稱PNLTCEVD)：

其中，ξj，μj，σj分別為xj(j=1,2,3)的邊緣分布的形狀參數(shù)、位置參數(shù)和尺度參數(shù)。通過矩估計法，得到相關(guān)參數(shù)α,β的顯式表達式：

式中：ri,j為一般的線性相關(guān)系數(shù)，i

PNLTCEVD成為顯式形式，通過對邊緣分布參數(shù)和相關(guān)參數(shù)的分別估計，可確定該模型的具體表達式。它結(jié)構(gòu)簡單，應(yīng)用方便，考慮了三種環(huán)境荷載下的最不利遭遇組合，同時涉及兩個相關(guān)參數(shù)，含概了變量之間相關(guān)的不對稱性。三個變量的分層結(jié)構(gòu)(α為外層，β為內(nèi)層)表明x1，x2之間比x1，x3和x2，x3之間有著更強的相關(guān)性。

上述理論解法及相應(yīng)的隨機模擬計算軟件，在國內(nèi)外尚無先例。只需輸入臺風(fēng)特征、各種致災(zāi)因素，進行邊緣分布分析，統(tǒng)計檢驗，即可使用相應(yīng)軟件進行臺風(fēng)災(zāi)害概率預(yù)測[20,37]。

4.1臺風(fēng)特征的長期聯(lián)合概率預(yù)測

由于每年臺風(fēng)頻次、路徑各不相同，在其形成、發(fā)展、登陸和消亡過程中，臺風(fēng)各種特征也不斷變化，以致在我國沿海和內(nèi)陸不同區(qū)域的臺風(fēng)特征組合也差異甚大。這些特征與致災(zāi)強度和后果都有直接關(guān)系。因此，研究不同區(qū)域臺風(fēng)特征組合及其相應(yīng)后果分析，是臺風(fēng)災(zāi)害區(qū)劃的重要組成部分。臺風(fēng)特征大多采用中心氣壓差(ΔP)、最大風(fēng)速半徑(Rmax)、臺風(fēng)移動速度(s)、臺風(fēng)中心與影響區(qū)域最小距離(δ)和臺風(fēng)移動方向角(θ)五種變量表示。復(fù)合極值分布理論的特點，在于將臺風(fēng)影響各區(qū)域的頻次(λ)作為一種離散型隨機量考慮。鑒于1975年尼娜臺風(fēng)在福建登陸，經(jīng)歷108個小時北上至河南省，導(dǎo)致暴雨—洪水巨災(zāi)，造成大量人員傷亡的悲劇，以及2006年碧利斯臺風(fēng)登陸后持續(xù)在內(nèi)陸各省造成的傷亡，本模式將臺風(fēng)登陸至消亡的歷時(t)也作為一種臺風(fēng)特征進行研究。對一種離散型分布(λ)和六種連續(xù)型分布的多維聯(lián)合概率研究，用多維復(fù)合極值分布理論的隨機模擬求解法，是可行的[35]。

4.2臺風(fēng)誘發(fā)不同致災(zāi)因素組合的聯(lián)合概率預(yù)測

核電站防護設(shè)施風(fēng)險評估程序如圖4所示。圖4中(b)、(c)分別為以隨機-數(shù)值模型計算建立的資料系列及觀測資料系列，對不同臺風(fēng)特征組合以及不同的核電站地理位置，分別對不同致災(zāi)控制因素(第一位)進行聯(lián)合概率預(yù)測[20]。在臺風(fēng)誘發(fā)風(fēng)暴潮、巨浪及相應(yīng)天文大潮以及暴雨、大風(fēng)等致災(zāi)因素的聯(lián)合概率預(yù)測中，引入整體不確定性，整體敏感性分析(GUA,GUA)[47]，在輸入-輸出過程中，反復(fù)調(diào)整輸入值，計算其對輸出的敏感程度，選擇最不利的組合，進行聯(lián)合概率預(yù)測，作為設(shè)防標(biāo)準(zhǔn)(圖5)。

圖4 核電站防護設(shè)施風(fēng)險評估Fig. 4 Safety analysis of disaster prevention infrastructures for nuclear power plants

4.3海岸帶防護工程結(jié)構(gòu)的可靠度分析

按可接受風(fēng)險分析原則(ALARP)進行防護工程結(jié)構(gòu)可靠度分析(圖4(d))。

按照國際上廣泛采用的ALARP準(zhǔn)則(as low as reasonable practice)[48]，結(jié)構(gòu)失效概率高于10-3為不可接受風(fēng)險，低于10-7為可接受風(fēng)險。使用多維復(fù)合極值分布理論的隨機模擬求解法；針對防護工程結(jié)構(gòu)特征，對防護工程結(jié)構(gòu)抗力因素及風(fēng)、浪、潮、流等外荷，使用多維復(fù)合極值分布理論的隨機模擬求解法，對海洋平臺，核電站進行防護工程結(jié)構(gòu)失效概率分析[30-32]。同樣也可用于越浪程度的分析，以保證其它工程結(jié)構(gòu)的安全[27-29]。

圖5 引入整體不確定性，整體敏感性分析(GUA,GUA)的聯(lián)合概率分析Fig. 5 Joint probability prediction with Global Uncertainty Analysis (GUA) and Global Sensitivity Analysis (GSA)

5 南海嶺澳核電站防護工程設(shè)防標(biāo)準(zhǔn)的風(fēng)險分析

嶺澳核電站位于南海濱。海岸帶臺風(fēng)誘發(fā)風(fēng)暴潮、巨浪與伴生的天文大潮“三碰頭”為主要致災(zāi)因素，在模型計算中對輸入—輸出體系引入整體不確定性(GUA)、整體敏感性(GSA)反復(fù)分析表明，臺風(fēng)誘發(fā)巨浪為致災(zāi)控制因素。按圖4所示程序計算如下。

5.1輸入因素不確定性分析

1) 天文大潮不確定性分析。通過調(diào)和分析和統(tǒng)計分析表明，每年的調(diào)和常數(shù)和平均水位都是不同，25年的調(diào)和常數(shù)和年平均水位進行不確定性和敏感性分析表明，調(diào)和常數(shù)和年平均海平面的變化每年不同，分別計算其平均值μ，均方差σ及變異系數(shù)COV[49]。

調(diào)和常數(shù)不確定性：

μ=166.65；σ=0.88；COV1=σ/μ=0.005 3

平均海平面不確定性：

μ=4.80；σ=3.96；COV2=σ/μ=0.825 6

天文大潮不確定性：

COV=(COV12+COV22)1/2=0.825

2)極值波高概率預(yù)測不確定性分析[44]

概率預(yù)測模式不確定性：COV1=0.08

取樣不確定性：COV2=0.09

統(tǒng)計不確定性：COV3=0.10

極值波高概率預(yù)測總體不確定性：COV=(COV12+COV22+COV32)1/2=0.156

3)置信區(qū)間估計公式[50]

Hp=b+eβ=φ(β)

而

故

ΔHp即為相應(yīng)某種重現(xiàn)值Hp的置信區(qū)間半長公式，置信區(qū)間(置信概率68.27%)為：

HT+ΔHT,HT-ΔHT

5.2南海海區(qū)臺風(fēng)概率特征分析

南海海區(qū)臺風(fēng)各參數(shù)符合的概率分布和概率特征參數(shù)、分布參數(shù)如表6所示。

表6 南海海岸帶臺風(fēng)概率特征Tab. 6 Marginal distribution of typhoon characters for South China Sea

5.3臺風(fēng)特征統(tǒng)計檢驗

臺風(fēng)特征統(tǒng)計檢驗圖如圖6(a)～圖6(e)所示。

圖6 臺風(fēng)特征統(tǒng)計檢驗圖Fig. 6 Statistical test chart of typhoon characteristics

5.4臺風(fēng)誘發(fā)暴潮、巨浪與伴生的天文大潮統(tǒng)計分布診斷檢驗

臺風(fēng)誘發(fā)暴潮、巨浪與伴生的天文大潮統(tǒng)計分布診斷檢驗，如圖7～圖9所示。各圖中圓點為實測資料點，實線為模型曲線；圖(c)中外圍兩條曲線為模型95%置信區(qū)間。

按照表5及圖6、圖7中的臺風(fēng)特征概率統(tǒng)計分析，可使用基于MCEVD的隨機模擬技術(shù)求解不同海區(qū)臺風(fēng)致災(zāi)因素的各種組合。經(jīng)過圖7～圖9的臺風(fēng)誘發(fā)巨浪、暴潮及相應(yīng)天文大潮的統(tǒng)計分布診斷檢驗，使用PNLTCEVD相應(yīng)的公式和軟件，即可獲得表7所示嶺澳核電站海區(qū)以巨浪為致災(zāi)控制因素的極端荷載組合，從而可對核安全規(guī)范制定的防護工程設(shè)計標(biāo)準(zhǔn)進行風(fēng)險評估。

圖7 臺風(fēng)暴潮診斷檢驗圖Fig. 7 Distribution diagnostic testing of storm surge

圖8 天文大潮診斷檢驗圖Fig. 8 Distribution diagnostic testing of spring tide

圖9 波高診斷檢驗圖Fig. 9 Distribution diagnostic testing of wave height

特征值重現(xiàn)期/a1005001000風(fēng)暴增水/m2．79+0．513．49+0．713．85+0．80天文大潮/m2．14+0．302．19+0．352．75+0．46波高/m6．6+0．37．3+0．67．9+0．8

該核電站按HAF0100規(guī)范的設(shè)計標(biāo)準(zhǔn)為[51]：

1)“設(shè)計基準(zhǔn)洪水”(DBF)=可能最大臺風(fēng)暴潮+天文大潮=6.35 m；

2)相應(yīng)最大波高=6.6 m。

雙層嵌套多目標(biāo)聯(lián)合概率模式計算結(jié)果為：

1)500年聯(lián)合重現(xiàn)期最大臺風(fēng)暴潮，天文大潮為4.20+2.54=6.74 m；

2)相應(yīng)最大波高為7.9 m。

圖10 天文大潮和風(fēng)暴增水聯(lián)合概率密度分布 Fig. 10 Probability density distribution of spring tide and storm surge

實例計算表明：使用MCEVD理論，計算取得500年一遇聯(lián)合出現(xiàn)的波高、暴潮和天文大潮的組合，遠遠高于核嶺澳電站按相應(yīng)“可能最大”組合的設(shè)防標(biāo)準(zhǔn)。

6 秦山核電站防護工程設(shè)防標(biāo)準(zhǔn)的風(fēng)險分析

據(jù)1952-2000年統(tǒng)計資料，秦山核電站附近乍浦、澉浦、鹽官等地最大潮差達 7.26～9.00 m。針對秦山核電站海岸防波堤9.76 m高程，進行以天文大潮為控制因素的臺風(fēng)誘發(fā)極端外荷作用下核電站海岸防護工程設(shè)防標(biāo)準(zhǔn)聯(lián)合概率風(fēng)險分析，結(jié)果如圖10～圖12以及表8所示。

圖11 天文大潮和風(fēng)暴潮累積頻率分布Fig. 11 Cumulative probability distribution of spring tide and storm surge

圖12 1 000年一遇天文大潮、風(fēng)暴潮和波高聯(lián)合概率分布Fig. 12 Joint probability distribution of spring tide, storm surge and extreme wave height with 1 000 years joint return period

聯(lián)合概率/a天文潮/m風(fēng)暴潮/m波高/m1004．23．02．55005．03．53．010005．54．03．5

秦山核電站海岸防波堤高度9.76 m，低于100年一遇以天文大潮為控制因素的臺風(fēng)誘發(fā)風(fēng)暴潮，波高極端外荷。

7 結(jié) 語

我國沿海多個核電站的防護工程設(shè)防標(biāo)準(zhǔn)仍然是采用仿照美國(NOAA) 的SPH、PMH的模式，使用“可能最大風(fēng)暴潮”(PMSS)”、“設(shè)計基準(zhǔn)洪水(DBF)”以及國際原子能工程機構(gòu)(IAEA)規(guī)定的“可能最大洪水(PMF)”等方法確定，忽略了這些定義和計算中存在的各種不確定性因素以及各種極端海況聯(lián)合出現(xiàn)的概率。以南海嶺澳核電站和秦山核電站防護工程為例，使用MCEVD理論，計算得出的500年一遇聯(lián)合出現(xiàn)的波高、風(fēng)暴潮和天文大潮的組合，遠遠高于兩個核電站按相應(yīng)“可能最大”組合的設(shè)防標(biāo)準(zhǔn)。按照我國水利工程設(shè)計標(biāo)準(zhǔn)，重要工程都需采用2 000年至3 000年一遇校核值，國際核電組織IAEA新建議采用萬年一遇的風(fēng)暴增水、天文大潮和波高組合作為防災(zāi)標(biāo)準(zhǔn)，而上述核電站的“可能最大”組合，連采用MCEVD計算得出的五百年一遇設(shè)計標(biāo)準(zhǔn)都達不到。面對全球極端氣候頻發(fā)的趨勢，建議對核電防護工程設(shè)防標(biāo)準(zhǔn)進行全面風(fēng)險評估，以達到國務(wù)院核電站安全規(guī)劃(2011-2020)提出的“必須按照全球最高安全要求”的規(guī)定。

致謝：衷心感謝謝世楞院士對本項研究的支持。

[1] LIU Y C, SANTOS A, WANG S M, et al. Tsunami hazards along Chinese coast from potential earthquakes in south China sea[J]. Phys Earth Planet In, 2007, 163:233-244.

[2] IAEA-TECDOC-1341, Extreme external events in the design or assessment of nuclear power plants[S]. Vienna: IAEA, 2003.

[3] IAEA-TECDOC-1487, Advanced nuclear plant design options to cope with extrenal events[S]. Vienna: IAEA, 2006.

[4] IAEA. Probability safety analysis (PSA) of other external events than earthquake[R]. Nuclear Energy Agency, Committee on the Safety of Nuclear Installations, NEA/CSNI, 2009.

[5] Design-basis flood estimation for site characterization at nuclear power plants in the United States of America[S]. U. S. NRC, 2011.

[6] HAD 101/09, 濱海核電廠廠址設(shè)計基準(zhǔn)洪水的確定[S]. 北京: 國家核安全局, 1990. (HAD 101/09, Determination of design basis flood level on coastal nuclear power plant site[S]. Beiing: National Nuclear Safety Administration of China, 1990. (in Chinese))

[7] HAD 101/11, 核電廠設(shè)計基準(zhǔn)熱帶氣旋[S]. 北京: 國家核安全局, 1991. (HAD 101/11. Design tropical cyclone for nuclear power plant[S]. Beiing: National Nuclear Safety Administration of China, 1991. (in Chinese))

[8] HAD 101, 濱海核電廠廠址選擇安全規(guī)定[S]. 北京: 國家核安全局, 1991. (HAD 101, Regulations for selection of nuclear power plant site[S]. Beijing: National Nuclear Safety Administration of China. 1991. (in Chinese))

[9] LIU D F, MA F S. Prediction of extreme wave heights and wind velocities[J]. Journal of the Waterway Port Coastal and Ocean Engineering, ASCE, 1980, 106: 469-479.

[10] LIU D F. Long term distribution of hurricane characteristics[C]//Proceedings of the Offshore Technology Conference. 1982: 305-313.

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

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

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

[14] 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.

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

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

[17] NAFFAA M G. Wave climate along the nile delta coast[J]. Journal of Coastal Research, 1995, 11:219-229.

[18] LIU D F, WANG L P. Theory of multivariate compound extreme value distribution and its application to extreme sea state prediction[J]. Chinese Science Bulletin, 2006, 51(23):2926-2930.

[19] LIU D F, PANG L, FU G, et al. Joint probability analysis of hurricane Katrina 2005[C]//Proceedings of the Intern. Offshore & Polar Eng. Conference (ISOPE2006). 2006: 74-80.

[20] LIU D F, PANG L, Xie B T, et al. Typhoon disaster zoning and prevention criteria -A double layer nested multi-objective probability model and its application[J]. Science in China Series E, Technological Sciences, 2008, 51(7): 1038-1048.

[21] CHOWDHURY A G, HUANG P, ERWIN J. Aerodynamic testing application of a full-scale facility for mitigating hurricane-induced coastal disasters[J]. Far East Journal of Ocean Research, 2009, 2(1):1-27.

[22] LIU D F, SHI H D. Disaster prevention design criteria for the estuarine cities: new orleans and shanghai-the lesson from hurricane katrina [J]. Acta Oceanologica Sinica, 2006, 25(4):124-130.

[23] 劉德輔,龐亮,史宏達. 卡特里娜颶風(fēng)的啟示——有關(guān)海岸和水利工程的風(fēng)險分析[J].中國工程科學(xué), 2007, 9(10):24-29. (LIU Defu, PANG Liang, SHI Hongda. The lesson from hurricane Katrina 2005—Risk analysis for coastal, offshore and hydraulic engineering[J]. Engineering Sciences, 2007, 9(10): 24-30. (in Chinese))

[24] LIU D F, XIE B T, LI H J. Study on the design flood volume of the three gorges project[J]. Journal of Hydrologic Engineering, ASCE, 2011, 16(1): 71-80.

[25] LIU D F, PANG L, LI H J. Extreme value prediction of typhoon induced disasters[C]//Proceedings of the 2nd International Summit on Hurricanes and Climate Changes, Monograph of invited papers in Summit. 2009: Paper No.21.

[26] LIU D F, JIANG Y P, SHI H D, et al. Uncertainty analysis of breakwater wave overtopping volume, wave forces and structure stability[C]//Proceedings of the ISOPE. 2006:387-392.

[27] LIU T F, ZHOU Z G, SHI J G, Optimum design environmental criteria for platforms[C]//Proceedings of the Int. Conference on Behaviour of Offshore Structures (BOSS). 1988, 3:1483-1495.

[28] SHI J G, LIU D F, WANG Z. The joint long-term distribution of design environmental factors for offshore structures multivariable probabilities and simulating method[J]. Acta Oceanologica Sinica, 1993, 12(3):465-474.

[29] LIU D F, YANG Y C, WANG C, et al. Joint probability of environmental loads on marine structures[C]//Proceedings of the OMAE. 1995: 408-415.

[30] LIU D F, DONG S, WANG C. Uncertainty and sensitivity analysis of reliability for marine structure[C]//Proceedings of the OMAE. 1996: 380-386.

[31] LIU D F, LIU G L, WANG F Q. Risk assessment of coastal defense against typhoon attacks for nuclear power plants[C]// Proceedings of the ICAPP. 2011: 2484-2492.

[32] LIU D F, LI H J, LIU G L, et al. Design code calibration of offshore, coastal and hydraulic energy development infrastructures[J]. WSEAS International Journal of Energy and Environment, 2011, 5(6):733-747.

[33] LIU D F, LIU G L, LI H J, et al. Discussion on IAEA and China safety regulation for NPP coastal defense infrastructures against Typhoon/Hurricane Attacks[J]. World Journal of Nuclear Science and Technology, 2012, 2:114-123.

[34] LIU G L,LIU D F, LI T F, et al. Joint probability prediction model of rainfall triggered landslides and debris flows[J]. Open Journal of Geology, 2012, 2(2): 103-110.

[35] LIU D F, LI H J, SHI H D, et al. Natural disasters: prevention, risk factors and management[M]. USA: NOVA Science Publishers, 2012: 1-86.

[36] LIU D F, SHI H D, LIU G L, et al. Risk assessment for nuclear power plants against natural disasters[M]. New York: NOVA Publisher, USA, 2015: 1-131.

[37] LIU D F, PANG L, XIE B T. Typhoon disaster in China-prediction, prevention and mitigation[J]. Natural Hazards, 2009, 49(3): 421-436.

[38] MITTAL A. Army corps of engineers: history of lake pontchartrain and hurricane protection project[R]. GAO-06-244T, 2006.

[39] BEA R. Reliability assessment and management：lessons from hurricane Katrina[C]//Proceedings of the Offshore Mech. & Arc. Eng.. 2007: OMAE2007-29650.

[40] 劉德輔, 史宏達,等. 中國沿海臺風(fēng)災(zāi)害區(qū)劃、防臺風(fēng)標(biāo)準(zhǔn)應(yīng)急制定、防臺風(fēng)應(yīng)急評估標(biāo)準(zhǔn)制定[R]. 青島: 中國海洋大學(xué), 2008. (LIU D F, SHI H D, et al. Typhoon disaster prediction, prevention and guideline of emergency plan for China Sea Coasts [R].Qingdao: Ocean University of China, 2008. (in Chinese))

[41] LIU D F, LIU G L,WANG F Q, et al. Extreme sea hazards statistics and engineering application[C]//Proceedings of the International Conference OCEANS` 14 MTS/IEEE. 2014.

[42] LIU D F, LIU G L,WANG F Q, et al. Risk assessment of design flood for three gorges dam[C]//Proceedings of the 6th International Conference on Larger Asian River. 2014.

[43] 劉德輔,史宏達,王風(fēng)清.對我國兩種海堤規(guī)范設(shè)計波浪標(biāo)準(zhǔn)的風(fēng)險分析[J]. 海洋工程, 2013, 31(3):70-76. (LIU Defu, SHI Hongda, WANG Fengqing, Risk analysis for coastal structure design criteria in China[J]. The Ocean Engineering, 2013,31(3):70-76. (in Chinese))

[44] MADSEN H, ROSBJERG D. The partial duration series method in regional index-flood modeling[J]. Wat Res Res, 1997, 33(4): 737-746.

[45] RANARD B, LANG M. Use of a Gaussian copula for multivariate extreme value analysis: some case studies in hydrology[J]. Adv Wat Res, 2006(30): 897-912.

[46] SHI D J, ZHOU S S. Moment estimation for multivariate extreme value distribution in a nested logistic model[J]. Ann Inst statist math, 1999, 51(2): 253-264.

[47] TARANTOLA S, GIGLIOLI N, JESINGHAUS J, et al. Can global sensitivity analysis steer the implementation of models for environmental assessments and decision-making? [J]. Stochastic Environmental Research and Risk Assessment, 2002, 16: 63-76.

[48] Marine risk assessment by Det Norsk Veritas, for the health and safety executive, offshore technology report[R]. 2001.

[49] LIU D F, KONG L S, ZUO J C, et al. Stochastic-numerical model of tidal current field for Jiaozhou bay of Yellow sea[C]//Proceedings of the ISOPE. 2001:682-685.

[50] 劉德輔,馬逢時.極值分布理論在計算波高多年分布中的應(yīng)用[J].應(yīng)用數(shù)學(xué)學(xué)報, 1976,1:23-37. (LIU D F, MA F S. Theory of extreme value distribution and its application on long term distribution of wave heights[J]. Journal of Applied Mathematics, 1976, 1: 23-37. (in Chinese))

[51] 謝世楞. 核電廠海域工程的設(shè)計標(biāo)準(zhǔn)問題[J].港口工程, 2000(1): 6-9. (XIE S L. Design criteria for nuclear power plants in coastal areas[J]. J. Port Engineering, 2000(1):6-9. (in Chinese))

Discussion on safety regulations for NPP coastal defense

LIU Defu, SHI Hongda, LIU Guilin, WANG Fengqing

(Ocean University of China, Qingdao 266003, China)

P753

A

10.16483/j.issn.1005-9865.2017.04.017

1005-9865(2017)04-0135-12

2016-11-14

國家自然科學(xué)基金資助項目(51379195)

劉德輔(1936-)，男，四川成都人，教授，主要從事海洋工程防災(zāi)減災(zāi)研究。E-mail：liu@ouc.edu.cn