黃 政，毛亞蔚，石雪垚，王 輝，陳巧艷

（1.中國(guó)核電工程有限公司，北京 100840；2.中國(guó)核工業(yè)二三建設(shè)有限公司，北京 101300）

第三代先進(jìn)壓水堆（如AP1000，APR1400和“華龍一號(hào)”）的設(shè)計(jì)中普遍采用了堆內(nèi)熔融物滯留（IVR）措施，運(yùn)行時(shí)通過(guò)淹沒(méi)堆坑來(lái)冷卻和維持壓力容器的完整性。但如果IVR 的排熱功率不足以帶走熔融物熱量時(shí)，大量高溫熔融物將通過(guò)破損的壓力容器噴射到堆坑，與冷卻水發(fā)生劇烈反應(yīng)，有可能觸發(fā)蒸汽爆炸。業(yè)界普遍認(rèn)為蒸汽爆炸發(fā)生的概率很低，但是一旦觸發(fā)，其產(chǎn)生的沖擊載荷將會(huì)對(duì)安全殼結(jié)構(gòu)的完整性造成威脅。蒸汽爆炸導(dǎo)致安全殼失效主要包括兩類(lèi)情形[1]：

（1）側(cè)壁失效，即堆坑環(huán)墻強(qiáng)度失效；

（2）上沖失效，即壓力容器向上移動(dòng)而引起的連接貫穿件斷裂。

其中，堆坑雖然距離安全殼殼體壁面還有一定距離，但是如果堆坑墻體發(fā)生結(jié)構(gòu)失效而坍塌，一方面會(huì)牽扯關(guān)聯(lián)的貫穿件和管道并導(dǎo)致這些部件的密封失效，引起放射性物質(zhì)的釋放；另一方面也會(huì)引起其所支撐的較重設(shè)備部件（例如蒸汽發(fā)生器等）發(fā)生位移并直接砸向安全殼殼體壁面，從而威脅到安全結(jié)構(gòu)的完整性[2,3]。因此這兩種情形均會(huì)破壞安全殼的密封性。因此有必要對(duì)蒸汽爆炸的載荷進(jìn)行計(jì)算分析，并對(duì)由此導(dǎo)致的安全殼失效的概率進(jìn)行定量評(píng)估。

為計(jì)算和評(píng)估蒸汽爆炸的后果，近年來(lái)業(yè)界開(kāi)發(fā)了多個(gè)計(jì)算分析程序，例如美國(guó)威斯康辛大學(xué)的TEXAS、日本JAERI 的JASMINE、法國(guó)IRSN 的MC3D 程序等。基于這些計(jì)算工具，多位學(xué)者采用數(shù)值模擬的方法對(duì)特定保守工況下的蒸汽爆炸特性和載荷進(jìn)行了分析。例如黃熙等[3]針對(duì)嶺澳二期電站研究了爆炸過(guò)程中堆腔不同位置的壓力和沖量的變化趨勢(shì)。鐘明君等[4]對(duì)1 000 MW 級(jí)壓水堆建立了三維幾何模型，研究了壓力和沖量的分布情況，以及破口位置和大小的影響。張蕊等[5]模擬分析了AP1000 堆外蒸汽爆炸壓力波的傳播特性和峰值強(qiáng)度。在得到爆炸載荷的基礎(chǔ)上，一些學(xué)者還進(jìn)一步開(kāi)展了安全殼結(jié)構(gòu)響應(yīng)的分析，評(píng)估了典型事故工況下的安全殼完整性。例如張娟花和陳鵬[6]，以及Zhang 等[7]用MC3D 計(jì)算得到的CPR1000 大破口事故下沖擊載荷分布作為輸入，利用有限元程序ABQUES 分析了安全殼結(jié)構(gòu)和主要關(guān)鍵設(shè)備的損傷程度。蒸汽爆炸的現(xiàn)象十分復(fù)雜、影響的因素很多，輸入?yún)?shù)和物理模型都存在一定的不確定性，會(huì)對(duì)安全殼完整性評(píng)價(jià)結(jié)論產(chǎn)生較大影響。但是之前開(kāi)展的計(jì)算分析大多是基于確定論方法針對(duì)典型工況下的爆炸載荷進(jìn)行計(jì)算，而對(duì)于輸入和模型參數(shù)的敏感性以及載荷的不確定性范圍的討論還比較少；另外，目前的計(jì)算模擬大都側(cè)重關(guān)注的還是爆炸載荷強(qiáng)度大小，而對(duì)于后續(xù)的安全殼響應(yīng)和后果，特別是對(duì)安全殼失效的概率還缺少定量化的評(píng)價(jià)。當(dāng)前，主要是由日本學(xué)者M(jìn)ORIYAMA 等人開(kāi)展了較為系統(tǒng)的研究[1]，提出了一套概率論分析流程框架方法，并基于JASMINE 程序針對(duì)日本的沸水堆和壓水堆開(kāi)展了失效概率評(píng)估和分析。該學(xué)者提出的研究方法可為后續(xù)針對(duì)其他核電廠設(shè)計(jì)開(kāi)展類(lèi)似的分析工作提供參考和借鑒。此外，當(dāng)前國(guó)內(nèi)也正在大力開(kāi)展第三代1 000 MW 級(jí)的先進(jìn)壓水堆電廠工程設(shè)計(jì)，而針對(duì)該電廠的嚴(yán)重事故條件下的蒸汽爆炸風(fēng)險(xiǎn)評(píng)價(jià)以及預(yù)防緩解措施的論證工作尚未充分開(kāi)展，因此也有必要開(kāi)展相關(guān)的風(fēng)險(xiǎn)評(píng)估工作。

本文采用概率論方法，提出了安全殼不同失效模式（包括側(cè)壁和上沖失效）所導(dǎo)致安全殼失效的概率，以及對(duì)總失效概率貢獻(xiàn)大小的計(jì)算方法。基于MC3D 程序，針對(duì)某1 000 MW先進(jìn)核電廠，對(duì)堆外蒸汽爆炸載荷將導(dǎo)致安全殼失效的概率水平進(jìn)行了定量評(píng)估。需要指出的是，目前關(guān)于蒸汽爆炸觸發(fā)的機(jī)理和判定標(biāo)準(zhǔn)的研究尚不明確，觸發(fā)蒸汽爆炸的概率很難準(zhǔn)確給出，仍存在較大的不確定性。而本文重點(diǎn)關(guān)注的是在已經(jīng)發(fā)生蒸汽爆炸的條件下的爆炸載荷對(duì)安全殼結(jié)構(gòu)完整性的影響和失效概率水平，因此在計(jì)算時(shí)人為觸發(fā)了堆外蒸汽爆炸。首先對(duì)蒸汽爆炸的初始條件和計(jì)算模型參數(shù)進(jìn)行了隨機(jī)采樣，分別計(jì)算得到了安全殼側(cè)壁失效和上沖失效兩種情形的載荷概率分布；然后采用理論推導(dǎo)的方法得到載荷極限值，進(jìn)而得到安全殼失效脆性曲線；最后通過(guò)結(jié)合載荷概率分布和安全殼脆性曲線，計(jì)算得到各個(gè)失效模式下以及總的安全殼失效概率值。本研究的方法和結(jié)論能夠?yàn)橄冗M(jìn)核電廠的安全評(píng)價(jià)和堆坑布置設(shè)計(jì)提供參考。

1 MC3D 程序和計(jì)算模型

1.1 計(jì)算程序

本文采用的是法國(guó)IRSN 開(kāi)發(fā)的多組分、多相流三維熱工水力計(jì)算程序MC3D V3.8[8]。該程序包括PREMIXING 和EXPLOSION 兩個(gè)計(jì)算模塊，分別用于計(jì)算粗混合破碎和爆炸過(guò)程。特別地，爆炸計(jì)算需要以粗混合計(jì)算的結(jié)果作為起始狀態(tài)。作為主流的蒸汽爆炸計(jì)算程序之一，該程序計(jì)算穩(wěn)定且速度較快，物理模型參與了OECD SERENA 項(xiàng)目的基準(zhǔn)實(shí)驗(yàn)驗(yàn)證[9]。

1.2 網(wǎng)格劃分和模型選擇

實(shí)際堆坑結(jié)構(gòu)較為復(fù)雜，本文進(jìn)行了必要簡(jiǎn)化處理，只保留壓力容器本體、堆坑環(huán)墻和上部自由空間。

利用幾何軸對(duì)稱特性，網(wǎng)格劃分采用二維柱坐標(biāo)系，在豎直Z方向劃分64 個(gè)節(jié)點(diǎn)，在半徑R 方向劃分了36 個(gè)節(jié)點(diǎn)，并在重要位置（如軸向中心線和壓力容器下方）進(jìn)行了加密（見(jiàn)圖1）。

圖1 MC3D 模型節(jié)點(diǎn)劃分示意圖Fig.1 Nodalization of the MC3D model

對(duì)于MC3D 的關(guān)鍵物理模型，熔融物物性材料選取的是程序內(nèi)置的CORIUM 材料。粗混合階段，計(jì)算熔融物破碎過(guò)程選用的是CONST模型，并且采用CAPAH2 模型來(lái)計(jì)算破碎過(guò)程中的氧化產(chǎn)氫。如前所述，由于蒸汽爆炸觸發(fā)的機(jī)理和概率尚不明確，本文重點(diǎn)關(guān)注的是發(fā)生蒸汽爆炸條件下的安全殼風(fēng)險(xiǎn)。因此為了保證蒸汽爆炸的發(fā)生，在計(jì)算時(shí)人為添加了壓力脈沖來(lái)觸發(fā)爆炸，選取熔融物液柱觸底時(shí)刻作為觸發(fā)爆炸計(jì)算時(shí)刻。

2 安全殼失效概率分析方法

2.1 計(jì)算方法流程

本文采用的失效概率計(jì)算流程（見(jiàn)圖2）主要分為三個(gè)步驟[10]：

圖2 安全殼失效概率計(jì)算流程Fig.2 The framework of evaluating the failure probability of the containment

（1）確定抽樣的參數(shù)及其分布后進(jìn)行抽樣，將抽樣結(jié)果組合為多組工況，輸入到MC3D 程序計(jì)算，從而得到載荷概率分布；

（2）確定失效模式和對(duì)應(yīng)載荷類(lèi)型后，采用理論分析方法計(jì)算載荷極限值并得到安全殼脆性曲線；

（3）結(jié)合前面得到的載荷概率分布和安全殼脆性曲線，通過(guò)概率方法計(jì)算安全殼失效概率。

2.2 變量選取和抽樣方法

為選取待抽樣物理變量，首先對(duì)蒸汽爆炸的主要進(jìn)程和現(xiàn)象進(jìn)行了梳理（見(jiàn)圖3）。特別關(guān)注的是可能對(duì)爆炸載荷強(qiáng)度造成影響的因素，因此主要是從兩個(gè)方面開(kāi)展分析：

圖3 蒸汽爆炸主要影響參數(shù)識(shí)別Fig.3 Identification of major influential factors of the steam explosion

（1）熔融物液柱破碎生成碎片的產(chǎn)生速率；

（2）熔融物與冷卻劑之間能量差。此外選擇變量時(shí)還同時(shí)參考了其他學(xué)者的選取方式[1,4,10,11]以及MC3D 的粗混合和爆炸物理模型說(shuō)明[8]。通過(guò)上述梳理并結(jié)合之前使用MC3D 開(kāi)展分析的工程經(jīng)驗(yàn)判斷，共選取了10 個(gè)物理量（分為初始條件和模型參數(shù)兩大類(lèi)）進(jìn)行分析。相應(yīng)的分布類(lèi)型和變化范圍的假設(shè)如表1 所示。由于蒸汽爆炸現(xiàn)象復(fù)雜、實(shí)驗(yàn)開(kāi)展和準(zhǔn)確測(cè)量極為困難，因此相關(guān)的實(shí)驗(yàn)研究數(shù)據(jù)稀少，也鮮有學(xué)者給出了重要參數(shù)的不確定概率分布形式。因此本文選取參數(shù)抽樣分布類(lèi)型時(shí)一方面是參考了其他學(xué)者的一般選取方式[1,10]，另一方面是結(jié)合工程經(jīng)驗(yàn)中關(guān)于某個(gè)參數(shù)上下限值范圍、及其分布的分散或集中特性的判斷來(lái)假設(shè)得到的。選取的原則主要包括：

表1 抽樣變量概率分布Table 1 Probability distribution of sampling variables

（1）對(duì)于可能趨向于某個(gè)參考值、且圍繞該值上下均勻波動(dòng)的參數(shù)，假設(shè)為正態(tài)分布；

（2）對(duì)于可能具有一定集中趨向性、但又與參考值有一定偏離特性的參數(shù)，假設(shè)為一些非對(duì)稱的分布形式（例如beta 分布）；

（3）對(duì)于有些具有明確分段特性的（例如一回路壓力），進(jìn)行分段假設(shè)，并且每段都采用均勻分布；

（4）對(duì)于其他難以確定分布特性、信息較少的參數(shù)，均假設(shè)為均勻分布，使之在其覆蓋范圍內(nèi)都等概率的出現(xiàn)[12]。

抽樣范圍的選擇結(jié)合了事故分析的結(jié)果，盡量覆蓋可能出現(xiàn)的變化范圍。特別地，一回路壓力分別對(duì)低壓（0.1～1.0 MPa）和高壓（1.0～5.5 MPa）熔堆工況進(jìn)行分段抽樣；兩段的概率與區(qū)間長(zhǎng)度成反比。此外安全殼壓力抽樣時(shí)需保證小于一回路壓力。對(duì)于熔融物質(zhì)量，一方面其分布形式的確定需要開(kāi)展大量的事故序列計(jì)算，另一方面質(zhì)量過(guò)小有可能無(wú)法觸發(fā)爆炸。因此從簡(jiǎn)化的角度考慮，本文保守假設(shè)事故條件下堆芯完全熔化，初始質(zhì)量取為本文所分析的參考電站的堆芯總質(zhì)量（恒定為152 t），不參與后續(xù)隨機(jī)抽樣計(jì)算。

為提高樣本分布的代表性，同時(shí)減少所需樣本的容量，采用了拉丁超立方抽樣方法（LHS）進(jìn)行隨機(jī)抽樣[13]，共形成了500 組計(jì)算算例樣本。

2.3 安全殼脆性曲線

安全殼脆性曲線（即失效概率-載荷曲線）的確定十分復(fù)雜，需要依賴大量實(shí)驗(yàn)或有限元計(jì)算。本文參考了文獻(xiàn)［1］的工程判斷方法，基于如下規(guī)則來(lái)確定脆性曲線：

（1）假設(shè)極限載荷為X，載荷達(dá)到X時(shí)，失效概率為95%。

（2）載荷達(dá)到0.5X時(shí)，失效概率為50%。

（3）以此作正態(tài)分布曲線，作為該載荷的脆性曲線。

需要指出的是，因?yàn)檎羝ǖ妮d荷均值可預(yù)期是一個(gè)較大的正值，因此使得相應(yīng)的安全脆性曲線分布均大部分落在載荷為正數(shù)的區(qū)間內(nèi)，而小于0 的部分對(duì)應(yīng)的概率則會(huì)非常低（可忽略不計(jì)）。另一方面，在采用公式（1）計(jì)算每種失效模式的概率時(shí)，自變量（載荷）的范圍也都是從有物理意義的正值開(kāi)始取的，因此也能保證本文在采用正態(tài)分布的條件下也不會(huì)產(chǎn)生不符合物理意義的負(fù)數(shù)情形。

該方法的核心是首先確定載荷的極限值X。對(duì)于安全殼的兩種失效情形，本文選取的主要失效模式和對(duì)應(yīng)的極限載荷如表2 所示。

表2 安全殼失效模式和對(duì)應(yīng)極限載荷Table 2 Containment failure modes and limited loads

計(jì)算極限載荷時(shí)，混凝土和鋼筋的材料強(qiáng)度值分別參考了國(guó)家標(biāo)準(zhǔn)混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范GB 50010—2010[14]和鋼結(jié)構(gòu)設(shè)計(jì)規(guī)范GB 50017—2003[15]。此外計(jì)算還考慮了應(yīng)變率效應(yīng)[16,17]，即材料在動(dòng)態(tài)載荷瞬時(shí)作用下強(qiáng)度等特性的提高，通過(guò)引入動(dòng)載增大系數(shù)（DIF）來(lái)體現(xiàn)。

2.4 安全殼失效概率計(jì)算

載荷概率分布曲線表征的是某載荷出現(xiàn)的概率大??；而脆性曲線則表示在某個(gè)水平的載荷作用下發(fā)生失效的概率，因此每種失效模式對(duì)應(yīng)的失效概率Pi可由全概率公式計(jì)算得到：

其中：P(F|Lj)——載荷Lj作用下的安全殼失效概率，可通過(guò)脆性曲線得到；

P(Lj)——載荷Lj發(fā)生的概率，由載荷概率分布曲線給出。

得到每種失效模式的概率后，安全殼的總失效概率Ptot就可通過(guò)以下公式計(jì)算得到：

由此可以進(jìn)一步推導(dǎo)出每一項(xiàng)失效模式對(duì)于總的失效概率的貢獻(xiàn)敏感度為：

該指標(biāo)可以定量反映出每種失效模式概率升高所引起的總失效概率增加的貢獻(xiàn)大小。

3 計(jì)算結(jié)果

3.1 爆炸載荷概率分布

圖4 給出了通過(guò)采樣計(jì)算得到的上述三類(lèi)載荷的概率分布直方圖（樣本容量為500），可以觀察到，在本文的計(jì)算模型和輸入?yún)?shù)分布假設(shè)下，得到的這三種類(lèi)型載荷分布都偏左側(cè)，大體上可以用伽馬分布來(lái)近似（見(jiàn)圖中擬合曲線）。表3 匯總出了載荷分布的主要統(tǒng)計(jì)量，其中最后一行數(shù)據(jù)對(duì)應(yīng)的是典型工況（表1 中的典型工況值）的結(jié)果。可以看出對(duì)于確定論方法通常選取的典型事故工況，由于考慮了較多保守假設(shè)，得到的爆炸后果也更為包絡(luò)；而如果考慮初始狀態(tài)和計(jì)算模型輸入的不確定性，大多數(shù)實(shí)際可能發(fā)生的爆炸情形產(chǎn)生的載荷遠(yuǎn)小于上述典型工況。

表3 蒸汽爆炸載荷概率分布統(tǒng)計(jì)量Table 3 Statistics of the probability distribution of loads

圖4 蒸汽爆炸載荷概率分布和擬合結(jié)果Fig.4 The probability distribution of loads and the fitting curve

3.2 安全殼脆性曲線計(jì)算

本文所分析的堆坑墻體混凝土的厚度為2.23 m，混凝土標(biāo)準(zhǔn)截面1 m×1.17 m 內(nèi)含有Φ40 的箍筋27 根。計(jì)算選取的主管道是尺寸較小的冷管段，外徑為0.85 m。首先分別計(jì)算各載荷極限值。對(duì)于側(cè)壁壓壞失效，采用面積加權(quán)法并乘以相應(yīng)動(dòng)載增大系數(shù)，計(jì)算得到的最大壓強(qiáng)極限值為731.5 MPa。對(duì)于拉斷失效，采用類(lèi)似方法并且考慮壓力波衰減效應(yīng)，得到最大沖量極限值為2.9 MPa·s。對(duì)于上沖失效，假定如果管道向上位移達(dá)到管徑的20%則發(fā)生斷裂失效，對(duì)應(yīng)的極限值為即 0.17 m。由此可根據(jù)2.3 節(jié)的規(guī)則得到相應(yīng)的安全殼脆性曲線（見(jiàn)圖5）；同時(shí)，為便于計(jì)算安全殼失效概率，圖5也同時(shí)畫(huà)出了相應(yīng)載荷的累積概率密度曲線。

圖5 安全殼脆性曲線和載荷累積概率Fig.5 The fragile curve and the cumulative probability of loads

圖5 安全殼脆性曲線和載荷累積概率（續(xù)）Fig.5 The fragile curve and the cumulative probability of loads

3.3 安全殼失效概率

結(jié)合前面的分析結(jié)果，可以得到各個(gè)失效模式對(duì)應(yīng)的安全殼失效概率（見(jiàn)表4）；然后由公式（2）可以得到在假設(shè)已經(jīng)發(fā)生了蒸汽爆炸的條件下，后續(xù)爆炸載荷所引發(fā)的安全殼失效的總概率為0.453。每種失效模式對(duì)安全殼總失效概率貢獻(xiàn)的敏感度則由公式（3）計(jì)算。

表4 各失效模式引發(fā)安全殼失效的概率Table 4 The failure probability and sensitivity

可以看到：側(cè)壁失效的兩種模式導(dǎo)致安全殼失效的概率基本相當(dāng)；而上沖失效對(duì)應(yīng)的失效概率則要高出許多。這是由于一方面?zhèn)缺诘匿摻罨炷翉?qiáng)度比較高，另一方面由于堆坑尺寸比較大，蒸汽爆炸壓力波傳播到側(cè)壁時(shí)強(qiáng)度衰減較為明顯，因而失效概率低。而壓力容器底部距離堆坑底板距離較短（～1.5 m），因而受到的向上的沖擊較大，導(dǎo)致與之連接的主管道被牽引上移的趨勢(shì)比較顯著，因此該失效模式是蒸汽爆炸可能導(dǎo)致安全殼失效的主要模式，敏感性貢獻(xiàn)也最大，需要給予特別關(guān)注。

總體而言，在確定發(fā)生了堆外蒸汽爆炸的條件下，爆炸載荷所引發(fā)安全殼失效的影響不可忽略，即一旦發(fā)生蒸汽爆炸，引起的安全殼失效的概率是比較高的。因此在核電廠設(shè)計(jì)中可以采取一些預(yù)防和緩解措施，避免蒸汽爆炸的發(fā)生或者降低爆炸的后果，例如：采用IVR策略保證事故條件下壓力容器的完整性，從而避免熔融物噴射到已被水淹沒(méi)的堆坑；或者進(jìn)一步改進(jìn)堆坑設(shè)計(jì)，增大尺寸，優(yōu)化幾何布局，針對(duì)壓力波影響較大的薄弱位置強(qiáng)化混凝土結(jié)構(gòu)強(qiáng)度等。針對(duì)本文分析的某先進(jìn)壓水堆，由于已經(jīng)采取了IVR 的預(yù)防策略并且具有較高的可靠性，因此由于IVR 失效導(dǎo)致壓力容器破損并進(jìn)而引發(fā)蒸汽爆炸的概率是很低的，所以堆外蒸汽爆炸風(fēng)險(xiǎn)的總體水平仍可認(rèn)為是很低的。

4 結(jié)論

本研究針對(duì)某1 000 MW 先進(jìn)壓水堆，根據(jù)蒸汽爆炸現(xiàn)象較為系統(tǒng)地梳理了重要影響因素和不確定性來(lái)源參數(shù)，采用概率論方法，提出了不同安全殼失效模式導(dǎo)致安全殼失效的概率計(jì)算方法，對(duì)確定已經(jīng)發(fā)生了蒸汽爆炸條件下，爆炸載荷導(dǎo)致安全殼失效的概率水平進(jìn)行了定量評(píng)價(jià)，并討論了不同失效模式對(duì)失效的貢獻(xiàn)，得到如下結(jié)論：

（1）本文計(jì)算得到了在已經(jīng)發(fā)生了蒸汽爆炸條件下，爆炸載荷將導(dǎo)致安全殼失效的總概率為0.453。因此一旦觸發(fā)了蒸汽爆炸，爆炸載荷造成安全殼損壞的潛在風(fēng)險(xiǎn)不可忽略，在設(shè)計(jì)中需要采取一定的預(yù)防和緩解措施來(lái)避免其發(fā)生或降低其后果影響，例如采用IVR 的事故預(yù)防策略、優(yōu)化堆坑隔室布局等。該1 000 MW先進(jìn)壓水堆設(shè)計(jì)中考慮了IVR，而IVR 失效并進(jìn)而引發(fā)蒸汽爆炸的概率是非常小的，所以該堆型堆外蒸汽爆炸風(fēng)險(xiǎn)的總體水平仍可認(rèn)為是足夠低的。

（2）相比于堆坑側(cè)壁失效，壓力容器上沖位移導(dǎo)致的貫穿件斷裂是主要的失效模式，這是由堆坑的材料、尺寸和布置特性決定的，需要予以特別關(guān)注。

本研究在獲取安全殼脆性曲線時(shí)采用了簡(jiǎn)化的方法，假設(shè)的合理性仍需進(jìn)一步通過(guò)更精細(xì)的有限元計(jì)算或?qū)嶒?yàn)方法來(lái)驗(yàn)證。但總體而言，本研究的方法和結(jié)論仍能夠?yàn)橄冗M(jìn)核電廠針對(duì)蒸汽爆炸風(fēng)險(xiǎn)的安全分析評(píng)價(jià)提供參考。