張 琨,史國寶,曹克美,王佳赟,蘆 葦,郭 寧
(上海核工程研究設(shè)計(jì)院有限公司,上海 200233)
熔融物堆內(nèi)滯留(IVR)是非能動核電廠的關(guān)鍵嚴(yán)重事故緩解策略之一。該策略在假想嚴(yán)重事故工況下,通過壓力容器與保溫層之間流道內(nèi)的水對壓力容器內(nèi)的熔融物進(jìn)行充分有效的冷卻,將堆芯熔融物滯留在壓力容器內(nèi),從而保證壓力容器的完整性,防止多數(shù)可能威脅安全殼完整性的堆外現(xiàn)象的發(fā)生。
IVR成功的熱工準(zhǔn)則為熔融物向壓力容器壁面?zhèn)鬟f的熱流密度小于壓力容器下封頭外壁面臨界熱流密度(CHF)。因此,熔融物衰變熱與CHF均為IVR有效性評價(jià)中的關(guān)鍵因素[1],需要分別開展壓力容器內(nèi)的熔池結(jié)構(gòu)傳熱研究以及壓力容器外的CHF研究。
壓力容器下封頭內(nèi)形成穩(wěn)定熔池時(對應(yīng)下封頭蒸干時刻),壓力容器保溫層流道入口流體條件(包括過冷度、壓力等)對CHF有重要影響。根據(jù)已有的試驗(yàn)結(jié)果[2-5],過冷度越大(入口水溫越低),則CHF結(jié)果越大;已有的部分試驗(yàn)結(jié)果表明入口壓力升高對于提高CHF結(jié)果是有利因素。
然而,目前已有的大部分試驗(yàn)主要假設(shè)不同的入口過冷度開展CHF的敏感性試驗(yàn),未對入口過冷度的取值進(jìn)行詳細(xì)論證。同時,壓力容器下封頭內(nèi)形成穩(wěn)定熔池時刻熔融物的衰變熱也是影響IVR有效性的重要因素,衰變熱越大,則熔池向壓力容器壁面?zhèn)鬟f的熱流密度越大,對IVR有效性挑戰(zhàn)也越大。是否存在熔融物衰變熱較高同時入口過冷度又較小的保守事故工況需要開展進(jìn)一步研究。另外,在穩(wěn)定熔池形成后,盡管熔融物衰變熱隨時間降低,但入口過冷度也逐漸減小,是否存在更為危險(xiǎn)的時刻也需要開展研究。
本文以大型非能動核電廠為研究對象,介紹壓力容器保溫層入口水過冷度對CHF影響的敏感性試驗(yàn)結(jié)果,在此基礎(chǔ)上分析典型嚴(yán)重事故工況下,壓力容器下封頭內(nèi)形成穩(wěn)定熔池時及以后堆腔水的過冷度以及熔融物衰變熱,研究兩者之間的對應(yīng)關(guān)系和規(guī)律,為IVR有效性分析提供依據(jù)。
IVR-ERVC全尺寸壓力容器下封頭外壁臨界熱流密度和流道流動試驗(yàn)(以下簡稱ERVC試驗(yàn))的研究目標(biāo)為:在全高度試驗(yàn)裝置上,進(jìn)行反應(yīng)堆下封頭外壁臨界熱流密度的試驗(yàn),掌握反應(yīng)堆壓力容器外部冷卻(ERVC)自然循環(huán)流動、沸騰傳熱特性,驗(yàn)證工程影響因素對CHF的影響機(jī)理,在此基礎(chǔ)上,通過優(yōu)化設(shè)計(jì),獲取相應(yīng)CHF工程數(shù)據(jù),為IVR有效性評價(jià)及相關(guān)安全審評提供支持。
ERVC試驗(yàn)臺架的設(shè)計(jì)特點(diǎn)如下所示。
(1)針對反應(yīng)堆壓力容器、保溫層、堆腔,以及實(shí)施IVR-ERVC時堆腔淹沒途徑的結(jié)構(gòu)與參數(shù),按高度方向1∶1切片形式,橫向重要截面1∶100比例進(jìn)行試驗(yàn)裝置設(shè)計(jì);
(2)關(guān)鍵的壓力容器下封頭—保溫層流道在形狀上完全模擬原型流道形式或可能的變形設(shè)計(jì)形式,而且自然循環(huán)流動回路阻力與原型保持一致;
(3)加熱段采用銅塊加壓力容器下封頭壁面原型材料的爆炸焊接方式加工而成,可模擬原型材料加熱表面;
(4)加熱段沿流道熱流密度分布可調(diào)節(jié),以模擬各種實(shí)際可能的壓力容器下封頭外壁熱流密度分布工況;
(5)壓力容器下封頭—保溫層流道進(jìn)水與排汽組件完全參照原型設(shè)計(jì),同時流通面積設(shè)計(jì)為可調(diào)整;
(6)借助于輔助與支持系統(tǒng),試驗(yàn)裝置還可模擬不同幾何形狀的ERVC加熱流道、不同入口水溫、不同水質(zhì)與水化學(xué)條件、不同加熱壁表面粗糙度與表面老化狀況等。
概括而言,整個試驗(yàn)設(shè)施平臺由以下幾部分組成:
(1)主裝置回路系統(tǒng),包括試驗(yàn)本體、下水箱及其內(nèi)部流道部件、上升管、上水箱及其內(nèi)部的排汽組件、下降水管、入口水室、入口水管等(見圖1);
圖1 ERVC試驗(yàn)臺架Fig.1 ERVC experiment facility
(2)輔助系統(tǒng),包括冷凝與冷卻系統(tǒng)、水化學(xué)調(diào)節(jié)系統(tǒng)、供水與預(yù)熱系統(tǒng)、回路清洗系統(tǒng)等;
(3)支持系統(tǒng),包括循環(huán)冷卻系統(tǒng)、造水系統(tǒng)、廢水處理系統(tǒng)、試驗(yàn)段準(zhǔn)備系統(tǒng)等;
(4)測控系統(tǒng)(測量、數(shù)據(jù)采集與控制系統(tǒng));
(5)電氣系統(tǒng)(大功率電源)等。
入口過冷度敏感性試驗(yàn)工況設(shè)置如下:
(1)試驗(yàn)裝置自由界面處為常壓;
(2)入口水溫分別取75 ℃、80 ℃、85 ℃、90 ℃、100 ℃;
(3)試驗(yàn)段表面為原型材料,長期在去離子水條件下沸騰運(yùn)行并洗滌;
(4)以試驗(yàn)加熱段最低點(diǎn)為0°,CHF測點(diǎn)包括27°、42°、57°、72°、81°、87°;
(5)試驗(yàn)工質(zhì):去離子水;
(6)試驗(yàn)流道幾何尺寸:與原型設(shè)計(jì)一致;
(7)試驗(yàn)段底部至自由液面垂直距離為8 m。
為更好反映入口過冷度敏感性試驗(yàn)結(jié)果的差異,以入口水溫100 ℃工況下最低傾角(27°)處CHF結(jié)果為基準(zhǔn),定義q′CHF如下:
(1)
即基準(zhǔn)點(diǎn)的q′CHF值為1。各試驗(yàn)工況q′CHF結(jié)果如圖2所示。在同一入口水溫條件下,CHF值隨角度增大而基本呈增大的趨勢。
圖2 CHF試驗(yàn)結(jié)果Fig.2 CHF experiment results
從圖2中可以看出,隨著保溫層入口水溫降低(過冷度增大),同一角度對應(yīng)的CHF值呈增大趨勢,表明提高入口過冷度有利于提高CHF值。
另外,在較低角度處(27°~42°),入口過冷度的CHF增強(qiáng)作用較顯著;中間角度處(57°~72°)入口過冷度的CHF增強(qiáng)作用減??;在較高角度處(81°~87°),入口過冷度的CHF增強(qiáng)作用已較小。究其原因,在較低角度處,過冷水直接接觸加熱表面,對于CHF影響較顯著,而在中高角度處,由于流體已經(jīng)過上游試驗(yàn)段加熱,溫度有一定上升,故入口過冷度對于CHF的影響減小。
基于CHF試驗(yàn)結(jié)果,對壓力容器下封頭形成穩(wěn)定熔池時的保溫層入口過冷度進(jìn)行計(jì)算分析,為IVR分析評價(jià)時選取CHF結(jié)果提供依據(jù)。
在IVR措施實(shí)施過程中,堆腔內(nèi)水的來源包括安全殼內(nèi)置換料水箱(IRWST)以及反應(yīng)堆冷卻劑系統(tǒng)(RCS)、堆芯補(bǔ)水箱(CMT)、安注箱內(nèi)的水等,其中最主要的來源為IRWST,由操縱員在堆芯出口溫度達(dá)到650 ℃時手動打開堆腔淹沒管線,將IRWST內(nèi)的水注入堆腔。對于IVR條件下堆腔內(nèi)水的過冷度,主要有以下幾種情況。
(1)IRWST重力注射失效,而且非能動余熱排出系統(tǒng)(PRHR)失效。在這種情況下,IRWST內(nèi)的過冷水通過堆腔淹沒管線注入堆腔,并與堆腔內(nèi)已有的水混合,使堆腔內(nèi)水溫較低,過冷度較高。該種情況的典型工況包括:事故類別3BE(反應(yīng)堆冷卻劑充分降壓但重力注射失效)以及事故類別6ES(蒸汽發(fā)生器傳熱管破裂始發(fā)事件導(dǎo)致的堆芯損傷序列)。
(2)IRWST重力注射失效,同時PRHR有效。在這種情況下,IRWST內(nèi)的水在事故前期受到PRHR的加熱,導(dǎo)致在操縱員實(shí)施堆腔淹沒操作時,注入堆腔內(nèi)的水溫度較高,過冷度較低。該種情況的典型工況包括:事故類別1AP(RCS部分降壓失效從而阻止IRWST重力注射的事故序列,假設(shè)操縱員手動開啟ADS閥門進(jìn)行RCS降壓)。
(3)IRWST重力注射有效,但由于再循環(huán)失效,堆芯在事故后期發(fā)生裸露和熔化。在這種情況下,事故瞬態(tài)進(jìn)程中IRWST中的水已用于對事故進(jìn)行一定的緩解,在堆芯內(nèi)已得到充分加熱,通過RCS破口流入堆腔,或通過自動卸壓系統(tǒng)(ADS)第4級閥門釋放到安全殼空間,使堆腔內(nèi)的水溫較高,過冷度較低。該種情況的典型工況包括:事故類別3BL(RCS充分降壓,IRWST重力注射有效,但長期再循環(huán)失效)。
采用MELCOR程序開展計(jì)算分析,該程序?yàn)橐惑w化程序,主要用于進(jìn)行輕水反應(yīng)堆事故進(jìn)程分析[6]。MELCOR程序的控制體建模方式和物理模型除廣泛用于常規(guī)壓水堆核電廠分析外,也可以靈活模擬非能動核電廠的各種安全系統(tǒng),包括堆芯補(bǔ)水箱、安注箱、安全殼內(nèi)置換料水箱重力注射等。另外,MELCOR程序還包括液膜跟蹤模型,可以有效模擬非能動安全殼冷卻系統(tǒng)(PCS)的運(yùn)行。
MELCOR程序核電廠模型的RCS部分主要包括壓力容器、兩臺蒸汽發(fā)生器、穩(wěn)壓器、穩(wěn)壓器波動管、4臺主冷卻劑泵、4條冷管段、兩條熱管段。二回路系統(tǒng)包括主給水系統(tǒng)、啟動給水系統(tǒng)、主蒸汽隔離閥、蒸汽發(fā)生器安全閥、汽輪機(jī)等。
堆芯及下腔室在軸向上分成15個節(jié)點(diǎn),其中下腔室占3個節(jié)點(diǎn);徑向共分成7環(huán);壓力容器下封頭沿軸向劃分為11個溫度節(jié)點(diǎn)。
專設(shè)安全設(shè)施模型包括兩臺CMT、兩臺安注箱(ACC)、兩條IRWST重力注射管線、兩條再循環(huán)管線、兩條堆腔淹沒管線、ADS第1級至第4級閥門、PRHR、PCS系統(tǒng)等。
安全殼模型對SG隔間、CMT隔間、堆腔、IRWST、上部隔間、非能動堆芯冷卻系統(tǒng)(PXS)隔間、化容系統(tǒng)(CVS)隔間和PCS穹頂進(jìn)行了詳細(xì)模擬。
計(jì)算得到各事故序列中壓力容器下封頭形成穩(wěn)定熔池時刻的堆腔入口過冷度及熔融物衰變熱如表1所示。圖3給出了各事故序列在壓力容器內(nèi)水蒸干時刻及其后續(xù)事故進(jìn)程中的保溫層入口過冷度與堆芯衰變熱的關(guān)系。其中堆腔入口飽和水溫度是根據(jù)安全殼壓力直接計(jì)算得到,未考慮水壓,在安全殼水淹條件下保溫層入口的飽和溫度應(yīng)比上述飽和溫度高10 K左右。
表1 堆芯衰變熱份額及堆腔過冷度Table 1 Decay heat fraction and subcooling degree of reactor cavity inlet fluid
圖3 入口過冷度與堆芯衰變熱的關(guān)系Fig.3 Inlet subcooling degree vs decay heat
從計(jì)算結(jié)果可以看出,壓力容器內(nèi)水蒸干時刻堆腔水過冷度較高的工況是3BE和6ES,但此時堆芯衰變熱相對于其他事故序列也較大。
1AP和3BL的壓力容器內(nèi)水蒸干時刻堆腔水過冷度較低,但同時堆芯衰變熱也較小。
圖3中左下角的空白區(qū)域表明IVR過程中不可能出現(xiàn)堆芯衰變熱很高,同時入口過冷度很低的工況。對于入口過冷度較高的工況,可保守選取入口水溫75 ℃的CHF試驗(yàn)結(jié)果,對于入口過冷度較低的工況,可保守選取入口水溫100 ℃的CHF試驗(yàn)結(jié)果進(jìn)行IVR有效性分析評價(jià)。
對于3BL、1AP這類形成穩(wěn)定熔池時堆腔水已接近飽和的工況,隨后續(xù)事故進(jìn)程發(fā)展,入口水溫保持不變,故CHF值仍取入口水溫100 ℃的試驗(yàn)結(jié)果。同時熔融物衰變熱逐漸減小,因此,在形成穩(wěn)定熔池時q/qCHF最大。
對于3BE、6ES這類形成穩(wěn)定熔池時堆腔水存在較高過冷度的工況,隨著事故進(jìn)程發(fā)展,熔融物衰變熱不斷減小,而過冷度同樣逐漸減小(見圖3),CHF值隨之降低。需要分析q/qCHF隨時間的變化趨勢,以確定形成穩(wěn)定熔池時刻以后是否存在更大的q/qCHF值。與6ES事故類別相比,3BE事故類別形成穩(wěn)定熔池時的熔融物衰變熱大得多,而堆腔入口過冷度相對差異不大(見表1),因此選取3BE事故類別為例,保守考慮如下:
(1)選取形成穩(wěn)定熔池時刻衰變熱最大的ADS第4級閥門誤開事故序列,分析形成穩(wěn)定熔池時刻至堆腔水飽和這段時間內(nèi)q和qCHF的關(guān)系(堆腔水飽和后CHF值不變,而衰變熱繼續(xù)降低,不會出現(xiàn)更危險(xiǎn)的時刻);
(2)假設(shè)為兩層熔池結(jié)構(gòu);
(3)形成穩(wěn)定熔池時刻堆腔水溫保守取75 ℃,假設(shè)過冷度在堆腔水飽和前隨時間線性遞減, CHF值將同步減小(見圖2),因此根據(jù)堆腔入口水溫80 ℃、85 ℃、90 ℃等的CHF試驗(yàn)結(jié)果確定qCHF曲線。
表2給出了壓力容器下封頭壁面熱流最大處(約82°)側(cè)壁的q/qCHF隨時間的變化。結(jié)果表明整個過程中形成穩(wěn)定熔池時刻是最危險(xiǎn)的,可針對該時刻進(jìn)行IVR有效性分析評價(jià)。
表2 形成穩(wěn)定熔池后的傳熱裕度分析結(jié)果Table 2 Analysis results of heat transfer margin after formed stable corium pool
形成穩(wěn)定熔池后q/qCHF的上升與保溫層入口水溫升高及隨之而來的CHF值減小直接相關(guān),而入口水溫的上升速率取決于注入堆腔的過冷水(主要來自IRWST)的量。大型非能動核電廠設(shè)計(jì)的堆腔注水源水量較多,則壓力容器壁面向外傳熱過程中堆腔水溫上升較慢,達(dá)到飽和的時間較晚,CHF值減小較緩慢,q/qCHF不會超過形成穩(wěn)定熔池時刻的數(shù)值。然而,如果堆腔水量較少,堆腔水溫將迅速上升,CHF減小的速率可能顯著超過衰變熱減小速率,則q/qCHF將快速增大,可能超過形成穩(wěn)定熔池時的數(shù)值,需要進(jìn)一步評估形成穩(wěn)定熔池后的IVR有效性。
本文介紹了入口過冷度對于壓力容器外壁面CHF的敏感性試驗(yàn)結(jié)果,并選取典型的嚴(yán)重事故序列,采用MELCOR程序計(jì)算分析了壓力容器下封頭內(nèi)形成穩(wěn)定熔池時刻及以后的保溫層入口過冷度,研究了其與熔融物衰變熱的關(guān)系。相關(guān)研究結(jié)論如下所示。
(1)壓力容器保溫層入口水的過冷度越大,壓力容器外壁面CHF越高,入口過冷度對于提高CHF是有利因素。
(2)入口過冷度對于CHF的作用在較低角度處比較顯著,隨著角度增大,其作用越來越小。
(3)在壓力容器下封頭內(nèi)形成穩(wěn)定熔池時刻,入口過冷度較大工況的熔融物衰變熱較高,而入口過冷度較小工況的熔融物衰變熱較低,不會出現(xiàn)入口過冷度小同時熔融物衰變熱高的情況。
(4)壓力容器下封頭內(nèi)形成穩(wěn)定熔池后,q/qCHF隨時間的變化趨勢與堆腔水量相關(guān),在水量較大的情況下,q/qCHF不會超過形成穩(wěn)定熔池時刻的值,形成穩(wěn)定熔池時刻可作為IVR分析評價(jià)的包絡(luò)狀態(tài)。