王 俊，田文喜，盧佳楠，蘇光輝，秋穗正

(西安交通大學(xué) 動(dòng)力工程多相流國(guó)家重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710049)

燃料棒包殼氧化是在高溫事故條件下影響堆芯行為的一個(gè)關(guān)鍵現(xiàn)象[1]。從20世紀(jì)開(kāi)始，國(guó)際上在這方面做了大量實(shí)驗(yàn)研究，構(gòu)成了鋯氧化動(dòng)力學(xué)的基礎(chǔ)[2-4]。鋯氧化動(dòng)力學(xué)由拋物線(xiàn)型的氧化定律來(lái)表示，并廣泛應(yīng)用于不同的嚴(yán)重事故分析代碼中[5-7]。目前國(guó)際上用于嚴(yán)重事故綜合分析的系統(tǒng)性程序有STCP、MELCORE、MAAP、ESCADRE、THALES等，機(jī)理性程序有SCDAP-RELAP5、CONTAIN、VICTORIA、CATHARE/ICARE等[8-11]。而國(guó)內(nèi)尚未開(kāi)發(fā)出受?chē)?guó)際廣泛認(rèn)可的嚴(yán)重事故方面的綜合分析程序。受?chē)?guó)家科技重大專(zhuān)項(xiàng)子課題《嚴(yán)重事故綜合分析軟件平臺(tái)開(kāi)發(fā)》的資助，本文目前正在進(jìn)行相關(guān)堆內(nèi)嚴(yán)重事故綜合分析程序的開(kāi)發(fā)，該程序主要包含早期行為、堆芯融化、碎片床、堆內(nèi)保持4部分內(nèi)容。燃料棒包殼氧化是嚴(yán)重事故序列中最早發(fā)生的事故[12-16]。本文主要介紹早期行為氧化模塊MIDAC-OX部分的程序結(jié)構(gòu)及模型，以AP1000[17]為例，計(jì)算全廠(chǎng)斷電事故后燃料芯塊融化、包殼破裂、包殼氧化以及氫氣產(chǎn)生等現(xiàn)象，并分析反應(yīng)堆冷卻劑泵惰轉(zhuǎn)時(shí)期的DNBR、自然循環(huán)時(shí)期對(duì)應(yīng)于不同程度冷卻劑自然循環(huán)流量下燃料棒包殼的完整性，以及燃料棒包殼氧化對(duì)其破裂的延遲作用。

1 氧化分析模塊程序的物理模型

1.1 反應(yīng)序列和氧化模型

反應(yīng)堆嚴(yán)重事故氧化過(guò)程主要受溫度影響，燃料棒在不同溫度下可能發(fā)生事故的序列列于表1。

表1 燃料棒在不同溫度下可能發(fā)生事故的序列[18]

材料氧化模型計(jì)算生成的熱、氫氣量和水蒸氣的變化量。該模型利用與材料溫度相關(guān)的氧化速率方程，假設(shè)材料氧化速率滿(mǎn)足拋物線(xiàn)方程：

(1)

式中：δ為增加的重量或?qū)雍穸龋琸g/m2或m；T為溫度，K；t為時(shí)間，s；A、B為取自MATPRO[19]的拋物線(xiàn)速率常數(shù)。

對(duì)于鋯合金，3個(gè)獨(dú)立的拋物線(xiàn)型方程可解出供氧量、α相鋯和ZrO2層的增加量。

氧化釋熱率Qox可根據(jù)重量的增量算出：

(2)

式中：MO2為氧氣的相對(duì)分子質(zhì)量，D；M為材料的相對(duì)分子質(zhì)量，D；hr為材料發(fā)生反應(yīng)釋放的反應(yīng)熱，J/kg；S為初始表面積，m2；w為單位表面積增加的氧化質(zhì)量，kg/m2。

最初的表面積由氧化前的面積決定，包括形變的影響。

氫氣產(chǎn)生速率ΔMH2和水蒸氣消除速率ΔMH2O可根據(jù)供氧量[19]求出：

(3)

(4)

氧化過(guò)程受3方面的限制。首先，當(dāng)物料完全氧化后，氧化過(guò)程終止，且氧化過(guò)程中，氧化速率會(huì)受到氧化層的抑制[20]。例如在鋯合金中，當(dāng)材料全部轉(zhuǎn)化為ZrO2而不銹鋼中當(dāng)鐵轉(zhuǎn)化為FeO2時(shí)，受攝氧率限制的關(guān)系如下：

(5)

式中：ρ為材料密度，kg/m3；V為材料體積，m3。

其次，氧化速率受限于可用蒸汽量，即：

≤

(6)

式中，m為氧化表面可用的蒸汽質(zhì)量流量。

最后，氧化受水蒸氣擴(kuò)散的影響。水蒸氣摩爾質(zhì)量流量受水蒸氣分壓力大小的驅(qū)動(dòng)，有：

(7)

式中：NH2O為水蒸氣摩爾質(zhì)量流量，kg·mol/s；As為表面積，m2；BH2O為換熱系數(shù)，W/(m2·K)；pH2O為水蒸氣分壓力，Pa。

利用質(zhì)量/熱傳遞類(lèi)比定律，假設(shè)關(guān)于傳熱和傳質(zhì)的Colburnj因素相等，則有：

(8)

式中：Nu為努塞爾數(shù)；Re為雷諾數(shù)；Sh為舍伍德數(shù)；Pr為普朗特?cái)?shù)；Sc為施密特?cái)?shù)。

利用質(zhì)量/熱傳遞類(lèi)比定律，氧化表面的攝氧率受限于如下關(guān)系式：

≤

(9)

式中：kH2O為水蒸氣熱導(dǎo)率，W/(m·K)；h為對(duì)流傳熱系數(shù)，W/(m2·K)；Dv為質(zhì)量擴(kuò)散速率，m2/s；R為氣體常數(shù)，J/(K·mol)；cp，H2O為水蒸氣的比定壓熱容，J/(kg·K)；ρH2O為水蒸氣密度，kg/m3。

質(zhì)量擴(kuò)散速率可根據(jù)下式計(jì)算：

(10)

式中：MH2O為水蒸氣相對(duì)分子質(zhì)量，D；Mi為第i種氣體的相對(duì)分子質(zhì)量，D；μH2O為水蒸氣的分子黏度，kg/(m·s)；μi為第i種氣體的分子黏度，kg/(m·s)；p為總壓力，Pa；X為凝結(jié)氣體總比重；Xi為第i種非凝結(jié)氣體比重；n為非凝結(jié)氣體數(shù)目。

1.2 程序結(jié)構(gòu)

嚴(yán)重事故綜合分析程序氧化模塊的流程圖如圖1所示。

圖1 氧化分析模塊程序流程圖

2 事故初始條件設(shè)定

2.1 流量完全喪失事故

失流事故中，冷卻劑流量與堆功率失配，導(dǎo)致堆芯材料包殼溫度迅速上升。其分類(lèi)包括流量部分喪失、流量完全喪失、主泵卡軸、主泵斷軸。其中，流量完全喪失事故是指由于全部主泵斷電或故障而惰轉(zhuǎn)的情況。本文所研究工況列于表2。

表2 某壓水堆核電廠(chǎng)流量完全喪失事故的事件序列[21]

2.2 冷卻劑流量控制模型

喪失全部流量事故瞬變分兩個(gè)階段。

第一階段：在瞬變開(kāi)始時(shí)，冷卻劑泵惰轉(zhuǎn)，其慣性壓頭比重力壓頭大得多，故此階段冷卻劑流量變化由冷卻劑泵的惰轉(zhuǎn)決定。但在此階段后期，重力壓頭份額逐漸增加，可認(rèn)為瞬態(tài)流量有一保守下限。

本文參考AP1000冷卻劑泵的惰轉(zhuǎn)特性，第一階段的冷卻劑流量變化規(guī)律如圖2所示。

圖2 AP1000冷卻劑泵惰轉(zhuǎn)時(shí)期歸一化流量

第二階段：泵的慣性壓頭消失，冷卻劑完全靠重力壓頭驅(qū)動(dòng)，即自然循環(huán)。自然循環(huán)可由兩種方式實(shí)現(xiàn)，一種是以一回路作為自然循環(huán)回路，其通過(guò)蒸汽發(fā)生器換熱，熱阱為二回路水；另一種是非能動(dòng)余熱排出系統(tǒng)的自然循環(huán)，其通過(guò)內(nèi)置換料水箱中的非能動(dòng)余熱排出熱交換器換熱，熱阱為內(nèi)置換料水箱。

本文假設(shè)第二階段自然循環(huán)能力(即冷卻劑質(zhì)量流量)為某一常數(shù)，不隨時(shí)間變化；且事故嚴(yán)重，出于某些原因，自然循環(huán)能力有不同程度的減弱。

2.3 功率控制模型

2.3.1軸向功率分布 燃料棒軸向歸一化功率分布列于表3。此軸向功率分布數(shù)據(jù)取自《非能動(dòng)安全先進(jìn)核電廠(chǎng)AP1000》[20]中圖2.24，將各控制體軸向高度內(nèi)的相對(duì)功率積分，并歸一化而得。

表3 燃料棒軸向歸一化功率分布

2.3.2徑向功率分布 參考AP1000幾種條件下1/8堆芯平面上的歸一化徑向功率分布，可知其最高功率倍數(shù)為1.321。燃料棒功率取平均功率的1.3倍。文獻(xiàn)[20]的圖2.18～2.23中，各工況下組件的最高功率/平均功率的最高倍數(shù)是1.321，故粗略地取危險(xiǎn)元件功率為平均功率的1.3倍。未考慮工程焓升因子等。

2.3.3緊急停堆控制 當(dāng)冷卻劑流量降至90%時(shí)觸發(fā)停堆保護(hù)，停堆信號(hào)延遲0.767 s。

2.4 相關(guān)設(shè)備和安全設(shè)施在事故中動(dòng)作

2.4.1冷卻劑泵 在核電廠(chǎng)失去外部電力負(fù)荷的Ⅱ類(lèi)事故中，仍有交流電源可運(yùn)行反應(yīng)堆冷卻劑泵，柴油機(jī)亦可在事故中為冷卻劑泵提供電力支持。但在惡劣的自然環(huán)境下，這兩種電力支持系統(tǒng)都有可能失效(如福島事故的地震加海嘯)，導(dǎo)致冷卻劑泵失去動(dòng)力，繼而開(kāi)始惰轉(zhuǎn)。

2.4.2蒸汽發(fā)生器 在核電廠(chǎng)失去外部電力負(fù)荷事故中，汽輪機(jī)旁路系統(tǒng)自動(dòng)打開(kāi)以排出核蒸汽供應(yīng)系統(tǒng)產(chǎn)生的蒸汽。如果凝汽器不可用，過(guò)量蒸汽會(huì)釋放到大氣中，然后啟動(dòng)給水系統(tǒng)來(lái)保持主給水流量。以上步驟的成功執(zhí)行，使得蒸汽發(fā)生器可持續(xù)冷卻一回路冷卻劑，從而建立一定程度的自然循環(huán)[16-17]。

2.4.3穩(wěn)壓器 在核電廠(chǎng)失去外部電力負(fù)荷事故中，堆芯冷卻能力下降，冷卻劑溫度升高，進(jìn)而一回路壓力升高。穩(wěn)壓器噴淋系統(tǒng)和安全閥可維持一回路壓力穩(wěn)定。

2.4.4停堆控制系統(tǒng) 事故中，一系列參數(shù)的超標(biāo)將會(huì)觸發(fā)反應(yīng)堆緊急停堆，如穩(wěn)壓器高壓、穩(wěn)壓器高水位、超溫ΔT、超功率ΔP和反應(yīng)堆冷卻劑泵低轉(zhuǎn)速。停堆后，堆芯功率迅速降至6%，從而減輕事故后果。

2.4.5非能動(dòng)余熱排出系統(tǒng) 在非LOCA事件時(shí)，非能動(dòng)余熱排出熱交換器(PRHR HX)將應(yīng)急排出堆芯余熱。該熱交換器由1組連接在管板上的C型管束和布置在上部(入口)和底部(出口)的封頭組成。PRHR HX的入口管線(xiàn)與反應(yīng)堆冷卻劑系統(tǒng)(RCS)熱管段相連接，出口管線(xiàn)與蒸汽發(fā)生器的下封頭冷腔室連接，他們與RCS熱管段和冷管段組成一非能動(dòng)余熱排出的自然循環(huán)回路。熱交換器位于高于RCS環(huán)路的內(nèi)置換料水箱內(nèi)，從而在反應(yīng)堆冷卻劑泵不可用時(shí)使冷卻劑依靠自然循環(huán)流過(guò)熱交換器。內(nèi)置換料水箱為熱交換器提供熱阱。

3 計(jì)算結(jié)果和討論

3.1 反應(yīng)堆冷卻劑泵惰轉(zhuǎn)時(shí)期燃料棒響應(yīng)

圖3為反應(yīng)堆冷卻劑泵惰轉(zhuǎn)時(shí)期燃料棒的響應(yīng)參數(shù)，包括歸一化冷卻劑流量和功率，不同高度燃料芯塊中心溫度、包殼表面溫度、冷卻劑溫度及最小DNBR等。由圖3可知，惰轉(zhuǎn)時(shí)期的前20 s內(nèi)，燃料芯塊中心最高溫升不到1 ℃，其中最高芯塊中心溫度為2 200 ℃，遠(yuǎn)未達(dá)到芯塊融化溫度；包殼溫度低于400 ℃，處在安全范圍內(nèi)；冷卻劑在3 s左右達(dá)到飽和溫度，出現(xiàn)兩相流，最小DNBR在4.5 s出現(xiàn)，為1.516，未發(fā)生偏離泡核沸騰。圖中各線(xiàn)注釋為燃料棒的軸向高度。如0.213為0.213 m高度處的計(jì)算結(jié)果。

圖3 惰轉(zhuǎn)時(shí)期的燃料棒參數(shù)

3.2 自然循環(huán)階段中的短期燃料棒響應(yīng)

在泵惰轉(zhuǎn)時(shí)期，緊急停堆后，功率迅速下降至滿(mǎn)功率的6%，包殼表面熱流密度急劇減小，偏離泡核沸騰不易發(fā)生。故在自然循環(huán)階段，需著重關(guān)注燃料棒的完整性，即包殼是否破裂。

緊急停堆后，堆芯余熱長(zhǎng)期保持在滿(mǎn)功率的0.3%。由于堆芯余熱長(zhǎng)期保持在0.3%以上，故自然循環(huán)流量須大于等于0.3%方能保證堆芯的長(zhǎng)期安全。本文取最低值來(lái)計(jì)算發(fā)生事故后的短期響應(yīng)，以此說(shuō)明燃料棒的安全性。圖4示出自然循環(huán)能力為正常流量的0.3%時(shí)的燃料棒響應(yīng)計(jì)算結(jié)果。

對(duì)于其他工況(自然循環(huán)流量為正常流量的1%、0.1%、0.01%、0.001%)，燃料棒響應(yīng)計(jì)算結(jié)果如圖5所示。

圖4 0.3%正常流量下燃料棒的參數(shù)變化

圖5 1%、0.1%、0.01%和0.001%正常流量下燃料棒的參數(shù)變化

圖6示出自然循環(huán)流量為正常流量的0.3%時(shí)，MIDAC-OX計(jì)算結(jié)果與其他程序計(jì)算結(jié)果的對(duì)比。結(jié)果顯示吻合程序較高，說(shuō)明MIDAC-OX程序的計(jì)算結(jié)果是合理可信的。

圖6 MIDAC-OX計(jì)算結(jié)果與文獻(xiàn)結(jié)果的對(duì)比

由圖5、6可知：1) 對(duì)于自然循環(huán)流量為正常流量0.01%的工況，雖然最高包殼表面溫度達(dá)2 342 ℃，但由于包殼氧化程度大于60%，使得包殼不易破裂，從而對(duì)包殼起到保護(hù)作用；2) 在自然循環(huán)流量為正常流量0.001%的工況，包殼表面溫度達(dá)1 852 ℃時(shí)未破裂，也是由于包殼氧化的保護(hù)作用；3) 在所計(jì)算的時(shí)間內(nèi)，在包殼破裂之前，芯塊溫度都未達(dá)到其融化溫度。本文未考慮應(yīng)力的影響。

3.3 自然循環(huán)階段中氧化對(duì)包殼的保護(hù)作用

包殼的氧化可使其熔點(diǎn)升高，增加其安全性，但同時(shí)鋯水反應(yīng)也會(huì)產(chǎn)生氫氣，危害安全殼的完整性。分別對(duì)不同工況(自然循環(huán)能力為正常流量0.008%降至0.001%)進(jìn)行分析計(jì)算，研究氧化對(duì)燃料棒的影響，結(jié)果示于圖7。

圖7 0.008%、0.006%、0.004%、0.002%和0.001%正常流量下燃料棒的參數(shù)變化

電站有自然循環(huán)流量的設(shè)計(jì)值，但自然循環(huán)流量的實(shí)際值會(huì)因某些原因而減小，故也需考慮這些情況。

由圖7可知：1) 對(duì)于自然循環(huán)流量為正常流量0.008%的工況，在386 s時(shí)包殼表面最高溫度達(dá)1 852 ℃，對(duì)應(yīng)的氧化程度為36.7%；在524 s時(shí)包殼表面最高溫度達(dá)2 227 ℃，對(duì)應(yīng)的氧化程度為76.5%；在935 s時(shí)包殼表面溫度達(dá)最大值2 573 ℃，此時(shí)包殼已完全氧化，無(wú)破裂風(fēng)險(xiǎn)。2) 對(duì)于自然循環(huán)流量為正常流量0.006%的工況，在347 s時(shí)包殼表面最高溫度達(dá)1 852 ℃，對(duì)應(yīng)氧化程度為33.6%；在439 s時(shí)包殼表面最高溫度達(dá)2 227 ℃，對(duì)應(yīng)氧化程度為64.2%；在647 s時(shí)，包殼表面最高溫度達(dá)2 687 ℃，對(duì)應(yīng)氧化程度為100%，達(dá)氧化層熔點(diǎn)，包殼破裂。3) 對(duì)于自然循環(huán)流量為正常流量0.004%、0.002%、0.001%的工況，包殼表面最高溫度達(dá)1 852 ℃的時(shí)間分別為316、288、276 s，對(duì)應(yīng)包殼氧化程度為31.2%、28.7%、28.2%，大于10%，故包殼不破裂；包殼表面最高溫度達(dá)2 227 ℃的時(shí)間為391、353、335 s，對(duì)應(yīng)包殼氧化程度為58.5%、54.2%、52.1%，未達(dá)60%，包殼破裂。4) 1 200 s內(nèi)單根棒包殼氧化的產(chǎn)氫量不超過(guò)400 μg。

4 結(jié)論

反應(yīng)堆完全喪失流量事故發(fā)生后，在冷卻劑泵惰轉(zhuǎn)時(shí)期，芯塊和包殼的溫度遠(yuǎn)未達(dá)到融化溫度，且最小DNBR在事故后4.5 s出現(xiàn)，為1.516，未發(fā)生偏離泡核沸騰。

在自然循環(huán)時(shí)期，對(duì)應(yīng)于不同程度的自然循環(huán)，包殼溫度先于芯塊溫度達(dá)到其限值。為保證燃料棒完整性，防止放射性物質(zhì)外泄，自然循環(huán)能力應(yīng)高于正常冷卻劑流量的0.01%。

雖然鋯水反應(yīng)會(huì)產(chǎn)熱產(chǎn)氫，但在事故發(fā)生后的較短時(shí)間內(nèi)，產(chǎn)氫量很小，燃料棒的氧化對(duì)包殼完整性起到保護(hù)作用——延遲包殼破裂時(shí)間。自然循環(huán)能力越弱，該延遲時(shí)間越短。

參考文獻(xiàn)：

[1] FICHOT F, ADROGUER B, VOLCHEK A, et al. Advanced treatment of zircaloy cladding high-temperature oxidation in severe accident code calculations, Part Ⅲ: Verification against representative transient tests[J]. Nuclear Engineering and Design, 2004, 232: 97-109.

[2] BOSLAND L, FUNKE F, GIRAUNT N, et al. PARIS project: Radiolytic oxidation of molecular iodine in containment during a nuclear reactor severe accident, Part Ⅰ: Formation and destruction of air radiolysis products—Experimental results and modeling[J]. Nuclear Engineering and Design, 2008, 238: 3 542-3 550

[3] BOSLAND L, FUNKE F, LANGROCK G, et al. PARIS project: Radiolytic oxidation of molecular iodine in containment during a nuclear reactor severe accident, Part Ⅱ: Formation and destruction of iodine oxides compounds under irradiation—Experimental results modeling[J]. Nuclear Engineering and Design, 2011, 241: 4 026-4 044.

[4] STEINBRüCK M, BIRCHLEY J, BOLDYREV A V, et al. High-temperature oxidation and quench behaviour of Zircaloy-4 and E110 cladding alloys[J]. Progress in Nuclear Energy, 2010, 52: 19-36.

[5] SCHANZ G, ADROGUER B, VOLCHEK A. Advanced treatment of zircaloy cladding high-temperature oxidation in severe accident code calculations, Part Ⅰ: Experimental database and basic modeling[J]. Nuclear Engineering and Design, 2004, 232: 75-84.

[6] DZHUSOV Y P, EFANOV A D. The effect of drop let sprinkling on the oxidation kinetics of zirconium cladding in steam[J]. Nuclear Engineering and Design, 2000, 195: 339-351.

[7] STEINBRUCK M, BIRCHLEY J, BOLDYREV A V. High-temperature oxidation and quench behaviour of Zircaloy-4 and E110 cladding alloys[J]. Progress in Nuclear Energy, 2010, 52: 19-36.

[8] KIM K T. The effect of fuel rod oxidation on PCMI-induced fuel failure[J]. Journal of Nuclear Materials, 2011, 418: 249-260.

[9] MITA S, AKEBI N, SAWAMURA O, et al. Development of the advanced thermal reactor in Japan[J]. Nuclear Engineering and Design, 1993, 144: 283-292.

[10] VOLCHEK A, ZVONAREV Y, SCHANZ G. Advanced treatment of zircaloy cladding high-temperature oxidation in severe accident code calculations, Part Ⅱ: Best-fitted parabolic correlations[J]. Nuclear Engineering and Design, 2004, 232: 85-96.

[11] DUREIZ C, STEINBRUCK M, OHAI D, et al. Separate-effect tests on zirconium cladding degradation in air ingress situations[J]. Nuclear Energy and Design, 2009, 239: 244-253.

[12] GROSSE M. Comparison of the high temperature oxidation behavior of different cladding materials[C]∥Jahrestagung Kerntechnik 2007: Tagungsbericht. Karlsruhe: INFORUM GmbH, 2007: 285-288.

[13] STEINBRUCK M. Oxidation of zirconium alloys in oxygen at high temperatures up to 1 600 ℃[J]. Oxidation of Metals, 2008, 70: 317-329.

[14] BEUZET E, LAMY J S, SIMONI E. Modelling of Zr-4 cladding oxidation by air, under severe accident conditions using the MAAP4 code[J]. Nuclear Engineering and Design, 2011, 241: 1 217-1 224.

[15] GAIKWAD A J, VIJAYAN P K, BHARTIYA S, et al. An innovative approach for steam generator pressure control of a nuclear power plant[J]. Nuclear Engineering and Design, 2011, 241: 5 113-5 126.

[16] GROSSE M, STEINBRUECK M, LEHMANN E, et al. Kinetics of hydrogen absorption and release in zirconium alloys during steam oxidation[J]. Oxidation of Metals, 2008, 70: 149-162.

[17] SCHULZ T L. Westinghouse AP1000 advanced passive plant[J]. Nuclear Engineering and Design, 2006, 236: 1 547-1 557.

[18] BUCK M. Modelling of the late phase of core degradation in light water reactors[M]. Stuttgart: University of Stuttgart Press, 2007.

[19] ALLISON C M, BERNA G A, CHAMBERS R, et al. SCDAP/RELAP/MOD3.1 code manual volume Ⅳ, MATPRO: A library of materials properties for light-water-reactor accident analysis[R]. Unite State: Idaho National Engineering Laboratory EG&G Idaho, Inc., 1993.

[20] 林誠(chéng)格. 非能動(dòng)安全先進(jìn)核電廠(chǎng)AP1000[M]. 北京：原子能出版社，2008.

[21] 朱繼洲. 核反應(yīng)堆安全分析[M]. 5版. 西安：西安交通大學(xué)出版社，2010.