基于經驗-解析數據的某核電安全殼地震易損性分析

王曉磊，呂大剛，閻衛東

(1.沈陽建筑大學 土木工程學院，遼寧 沈陽 110168；2.河北省地震災害防御與風險評價重點實驗室，河北 三河 065201；3.哈爾濱工業大學 土木工程學院，黑龍江 哈爾濱 150090)

核電廠地震安全評估包括抗震裕量評估(SMA)[1]和地震概率風險評估(SPRA)[2]兩種方法。SPRA是基于全概率公式的核電廠地震安全評估方法，可定量分析核電廠地震安全評估中包含的各種不確定性，分析結果較SMA更為精細。由于地震危險性分析技術的進步，已運行核電廠廠址地震危險性水平可能發生變化，需對該核電廠廠址上的核電廠結構、系統和部件重新進行地震安全評估[3]。

地震易損性分析是SPRA中的重要組成部分。目前安全系數法是核電廠常用的地震易損性分析方法，其數據來源主要包括解析易損性數據、經驗易損性數據、專家判斷易損性數據和混合易損性數據等。文獻[4]指出混合易損性數據由于綜合多種數據源的優點，相對其他數據源更具優勢。我國已運行了一些核電廠，由于地震危險性分析技術進步和實際地震發生等原因，這些已服役核電廠廠址地震危險性(審查水平地震和廠址譜等)可能發生變化?；诨旌蠑祿牡卣鹨讚p性方法研究，對我國已運行核電廠地震安全評估具有重要借鑒意義。

本文基于經驗-解析數據對某核電安全殼地震易損性分析進行研究。

1 核電廠地震易損性分析

1.1 地震易損性安全系數法

地震易損性安全系數法[5]是目前核電廠結構、系統和部件常用方法，基本原理可簡述如下。

核電廠結構、系統和部件抗震能力A[6-9]可表示為：

A=AmeReU

(1)

式中：Am為核電廠結構、系統和部件抗震能力中位值；eR和eU為假設服從對數正態分布的隨機變量(中位值為1，標準差分別為βR和βU)，分別代表本質不確定性和知識不確定性。

安全系數法將抗震能力表示為安全系數與設防地震能力乘積的形式：

A=F·ASSE

(2)

式中：F為抗震能力安全系數；ASSE為安全停堆地震能力。

安全系數F進一步可表示為：

F=FS·Fμ·FRS

(3)

式中：FS為強度系數；Fμ為塑性能吸收系數；FRS為響應系數。

(4)

式中：S為某種失效狀態下單元抗震能力；RN為非地震荷載下的單元響應；RT為所有荷載作用下(包括地震荷載)的單元響應。

FRS=FSA·FSD·FM·FMC·Fδ·FEC·FSSI

(5)

式中：FSA為譜型系數；FSD為隨著深度地震動輸入折減系數；FM為建模系數；FMC為模態組合系數；Fδ為阻尼系數；FEC為地震分量的組合系數；FSSI為土-結相互作用系數。

(6)

(7)

(8)

βFRS=

(9)

1.2 基于混合易損性數據的核電安全殼地震易損性分析方法

本文基于地震易損性安全系數法，采用混合易損性數據對服役核電廠安全殼進行地震易損性分析，具體步驟為：生成核電廠場地相關條件譜，選取匹配場地相關譜的地震動記錄，進行安全殼增量動力分析(IDA)，IDA結果為安全系數FS和FSA的中位值和標準差，安全系數法中的其他系數采用經驗數據，將解析易損性數據和經驗易損性數據整合為混合易損性數據，最終生成核電廠安全殼地震易損性曲線。

2 基于經驗-解析數據的我國某核電安全殼地震易損性分析

2.1 某核電安全殼力學模型

本文分析對象為我國華南地區某核電廠核島廠房，其三維剖面圖和簡化模型[10]如圖1所示，核島中安全殼和筏板基礎簡化模型信息及材料參數參考文獻[10]。本文采用開源程序OpenSees進行有限元建模和分析，模態分析結果列于表1。

a——剖面圖；b——簡化模型

表1 安全殼簡化模型模態分析結果

安全殼是核電廠最后一道安全屏障，其完整性十分重要。本文假設當在地震中安全殼處于塑性狀態時，安全殼即失效。安全殼應力應變關系可表示為三線性骨架曲線模型[11]，如圖2所示，取第一個拐點為安全殼極限狀態，此剪應力極限值τ1[11]為：

(10)

式中：FC為混凝土抗壓強度；σV為豎向壓應力。

2.2 某核電廠廠址條件譜和地震動記錄選取

2.2.1某核電廠廠址條件譜 場地相關譜是地震危險性分析的副產品，其中，一致危險譜是早期被廣泛應用的場地相關譜，有學者發現一致危險譜較為保守，在一致危險譜基礎上，Baker等[12]提出了條件均值譜。條件均值譜相較于一致危險譜進一步考慮了譜型相關性，更符合真實地震動反應譜形式。

圖2 三線性骨架曲線

條件均值譜μln Sa(Ti)|ln Sa(T*)的理論公式[12]為：

μln Sa(Ti)|ln Sa(T*)=μln Sa(M,R,Ti)+

ρTi,T*ε(T*)σln Sa(Ti)

(11)

式中：μln Sa和σln Sa為基于地震動預測方程的ln Sa(Ti)的預測平均值和標準差，Sa(Ti)為周期為Ti的譜加速度；ρTi，T*為相關系數；M為震級；R為距離；T*為條件周期；Ti為預測周期；ε(T*)為條件周期T*處的譜型參數，ε(T*)可以表示為：

(12)

在條件均值譜的基礎上，進一步考慮條件標準差可生成條件譜，條件譜考慮了場址譜的分布(條件標準差)信息。條件標準差σln Sa(Ti)|ln Sa(T*)[13]可表示為：

(13)

將我國華南地區某核電廠廠址作為算例場地，該場地包含1個地震統計區：空間范圍為北緯19°～24°、東經109°～116°，Mmax(最大震級)為8.00，b值(古登堡-里克特公式參數)為0.87，ν4(四級以上地震年平均發生率)為5.60，h(震源深度)為15 km，包含32個潛在震源區，如圖3所示。

圖3 廠址潛在震源區分布

地震動預測方程(GMPE)采用霍俊榮[14]分析給出的華南地區GMPE，可表示為：

lgY=C1+C2M+

C3lg(R+C4exp(C5M))+σlg Yε

(14)

式中：C1、C2、C3、C4和C5為系數；ε為服從標準正態分布的一個隨機變量，代表lgY的觀測變量；σlg Y為GMPE的預測標準差。

圖4 廠址地震危險性曲線

表2 平均值設定地震(萬年一遇)

2.2.2地震動記錄選取 本文采用Baker等[16]提出的基于貪心優化算法的地震動記錄挑選方法，具體理論如下。

選取的地震動與目標譜的誤差可由地震動對數譜加速度與目標譜對數譜加速度差值的平方和SSE表示：

(15)

式中：ln Sa(Tj)為被選地震動調幅后在周期為Tj的對數譜加速度；ln Sa(s)(Tj)為模擬目標譜在周期為Tj的對數譜加速度；p為周期個數。

圖5 算例核電廠條件譜

被選地震動對數譜加速度平均值與目標譜對數譜加速度平均值的最大百分數誤差ERRmean可表示為：

(16)

式中：mln Sa(Tj)為被選地震動在周期為Tj的對數譜加速度的樣本平均值；μln Sa(Tj)為場地目標譜在周期為Tj的對數譜加速度平均值。

被選地震動對數譜加速度標準差與目標譜對數譜加速度標準差的最大百分數誤差ERRstd可表示為：

(17)

式中：sln Sa(Tj)為被選地震動在周期為Tj的對數譜加速度的樣本標準差；σln Sa(Tj)為目標譜在周期為Tj的對數譜加速度的標準差。

上述兩類誤差的權重誤差平方和SSEs可表示為：

w(sln Sa(Tj)-σln Sa(Tj))2]

(18)

式中，w為權重系數，表示平均值和標準差兩類誤差的相對重要性。

基于貪心優化算法的地震動選取方法步驟簡述如下[16]：首先以條件譜(式(11)、式(12)和式(13))為目標譜，基于統計模擬方法生成模擬譜，然后基于最小SSE準則(式(15))初選地震動記錄，判斷選取地震動記錄是否滿足式(16)和式(17)誤差要求(容許誤差為10%)，如不滿足，基于式(18)并運用貪心優化算法進一步選取地震動記錄，反復上述過程，直到選取地震動記錄滿足式(16)和式(17)為止，最終輸出地震動選取結果。

以圖5所示的條件譜分別為目標譜，基于上述地震動記錄選取方法，最終挑選了兩組地震動記錄，挑選的地震動記錄與目標譜擬合情況分別如圖6、7所示，可發現：選取的地震動記錄可很好匹配目標譜的均值和標準差。

2.3 基于IDA的解析易損性數據

本文采用IDA方法得到易損性函數的中位值和標準差，易損性函數中位值和標準差的估計值[17]可表示為：

(19)

(20)

基于本文選取的兩組地震動記錄，采用IDA方法生成擬合易損性曲線，如圖8a、b和圖9a、b所示，得到解析易損性數據，列于表3。將IDA方法計算的地震動解析易損性數據與其他參數經驗易損性數據(表4)結合，可生成易損性曲線，如圖8c和圖9c所示。

a——模擬地震動的反應譜；b——選取地震動的反應譜；c——中位值Sa；d——對數Sa的標準差

a——模擬地震動的反應譜；b——選取地震動的反應譜；c——中位值Sa；d——對數Sa的標準差

a——IDA；b——基于IDA結果擬合易損性曲線；c——地震易損性曲線

a——IDA；b——基于IDA結果擬合易損性曲線；c——地震易損性曲線

表3 IDA擬合的解析易損性模型參數

表4 經驗地震易損性數據[7,18]

2.4 應用建議

我國已經運行了一些核電廠，由于地震危險性分析技術的進步，已運行核電廠廠址的地震危險性水平可能發生變化。采用解析方法(包括IDA等)生成FS和FSA兩個系數(該系數受地震危險性水平影響較大)的解析數據，其他地震危險性水平影響較小的參數采用經驗易損性數據，最終得到核電廠安全殼地震易損性曲線。由于地震危險性分析水平進步或實際地震的發生，我國已運行核電廠廠址地震危險性水平可能發生變化，針對上述情況，建議可采用經驗-解析數據的地震易損性分析方法初步分析我國核電廠安全殼地震易損性。

3 結論

本文生成了我國某核電廠廠址條件譜，選取了符合廠址危險性的地震動記錄，運用IDA方法生成了安全系數FS和FSA的解析易損性數據，連同其他系數的經驗易損性數據，采用經驗-解析易損性分析方法生成了我國已運行某核電廠安全殼地震易損性曲線，并給出了經驗-解析易損性分析方法在我國核電廠安全殼地震易損性中的應用建議，得出以下結論。

1) 條件譜在條件均值譜的基礎上，進一步考慮了廠址譜的變化范圍，能夠表示廠址譜的不確定性。

2) 基于貪心優化算法選取的地震動記錄能夠匹配場地的危險性。

3) IDA方法可生成地震易損性安全系數法中系數FS和FSA的解析易損性數據。

4) 對于場地危險性變化的核電廠廠址，采用基于混合易損性(經驗-解析)數據的地震易損性分析方法可方便生成核電廠安全殼地震易損性曲線。

5) 當我國已運行核電廠廠址地震危險性發生變化時，建議將基于經驗-解析數據的地震易損性分析方法應用在我國核電廠安全殼初步地震易損性分析中。