姚進國，李旭東，楊曉強，李載鵬，楊高升

(中國核電江蘇核電有限公司，江蘇 連云港 222042)

中國核電有限公司旗下的田灣核電站一期工程兩臺機組采用俄羅斯的WWER-1000/428型反應(yīng)堆裝置，并于2007年5月和8月兩臺機組相繼投入商業(yè)運行，設(shè)計循環(huán)長度7 000 h，換料周期為年度換料，負(fù)荷因子80%，與國內(nèi)外核電機組相比，其運行經(jīng)濟性偏低。為了提高田灣核電站機組運行的負(fù)荷因子及核電站經(jīng)濟性，田灣核電站通過引進高性能TVS-2M燃料組件，采用18個月燃料循環(huán)換料策略，并確定在一、二號機組第8燃料循環(huán)逐步開始向18個月燃料循環(huán)過渡。

2010年12月，TVS-2M先導(dǎo)組件的方案設(shè)計及堆芯論證工作完成，并獲得國家核安全局的批準(zhǔn)。一號機組第5燃料循環(huán)反應(yīng)堆堆芯由157個AFA燃料組件和6個TVS-2M先導(dǎo)燃料組件組成，目前6個先導(dǎo)組件正運行在第7燃料循環(huán)，運行考驗情況良好，為田灣核電站一、二號機組實施長周期換料策略奠定了基礎(chǔ)。鑒于TVS-2M燃料組件與目前田灣核電站使用的AFA組件在結(jié)構(gòu)上存在較大差異，為了確保過渡堆芯的安全和可靠運行，通過在機組上的試驗測量來驗證TVS-2M燃料組件的性能及與AFA燃料組件的相容性，同時驗證熱工水力設(shè)計程序計算模型的合理性以及計算結(jié)果與試驗結(jié)果的符合性。

1 TVS-2M先導(dǎo)組件特點

與AFA組件相比，TVS-2M組件在設(shè)計上的改進主要包括：1) 定位格架與導(dǎo)向管采用焊接方式；2) 定位格架數(shù)量減少，格架外條帶的厚度和寬度增加；3) 燃料棒活性段兩端長度增加；4) 燃料芯塊外徑增加，芯塊中心孔減小等。通過上述改進，提高了TVS-2M組件的結(jié)構(gòu)穩(wěn)定性，并增加了燃料裝載量。燃料最大運行時間延長到40 000 h，設(shè)計的最大燃耗為60 MW·d/kgU[1-2]。TVS-2M燃料組件與AFA燃料組件參數(shù)對比列于表1。

表1 TVS-2M與AFA燃料組件主要參數(shù)對比

2 穩(wěn)態(tài)熱工水力分析

熱工水力分析的目的是驗證反應(yīng)堆穩(wěn)態(tài)工況下的參數(shù)是否滿足熱工水力設(shè)計準(zhǔn)則，即在滿足“雙95%”條件下，燃料棒外表面不發(fā)生DNB；燃料包殼外表面最大溫度不超過352 ℃；燃料芯塊最高溫度應(yīng)低于熔化溫度。穩(wěn)態(tài)熱工分析的主要參數(shù)包括通過反應(yīng)堆和燃料組件的流量、燃料組件的冷卻劑溫升、堆芯壓降、燃料棒包殼外表面最大溫度、堆芯出口冷卻劑過冷度以及最小DNBR等。

2.1 計算程序簡介

穩(wěn)態(tài)熱工水力分析采用STAR-1和PUCHOK-1000穩(wěn)態(tài)計算程序。STAR-1用于計算一回路冷卻劑溫度、流量、壓降、燃料組件冷卻劑流量及比焓分布等總體參數(shù)。STAR-1將反應(yīng)堆和一回路視作1個水力網(wǎng)絡(luò)，網(wǎng)絡(luò)每段范圍內(nèi)冷卻劑都是一維流動，在堆芯段和蒸汽發(fā)生器管束段內(nèi)都是一維熱交換[3]。

PUCHOK-1000計算程序是計算燃料棒束的穩(wěn)態(tài)子通道熱工水力分析程序，根據(jù)沿棒束截面和高度的冷卻劑局部參數(shù)(比焓、質(zhì)量流密度)，計算臨界熱流密度和表面?zhèn)鳠嵛C裕度及棒束內(nèi)表面?zhèn)鳠嵛C裕度因子的最小值，并確定其所在的計算坐標(biāo)(網(wǎng)格序號和沿高度的段序號)[4]。

PUCHOK-1000程序考慮堆芯燃料組件的徑向、軸向功率分布，非均勻幾何，子通道的湍流攪混以及燃料棒與冷卻劑之間的傳熱。在求解計算棒束內(nèi)局部熱工水力參數(shù)的方程組時，采用了平分間隔法以及利用松弛因子的簡單迭代法，以保證求解能量方程所需的穩(wěn)定性。

采用的臨界熱流密度關(guān)系式是根據(jù)帶格架的棒束實驗數(shù)據(jù)歸納得出的Gidropress關(guān)系表達式：

qDNB=0.95(1-x)-0.95+0.105p·

式中：qDNB為均勻功率分布條件下用Gidropress關(guān)系式計算的均勻臨界熱流密度，MW/m2；x為發(fā)生DNB點的相對含汽率，范圍為-0.07～0.4；p為壓力，范圍為7.45～16.7 MPa；ρw為質(zhì)量流速，范圍為700～3 800 kg/(m2·s)。

2.2 穩(wěn)態(tài)計算輸入

計算中考慮了最保守的參數(shù)偏差組合，4臺冷卻劑泵運行工況下，考慮參數(shù)偏差的反應(yīng)堆熱工水力主要輸入?yún)?shù)列于表2。

表2 反應(yīng)堆熱工水力設(shè)計主要輸入?yún)?shù)

3 計算結(jié)果

3.1 AFA燃料組件與TVS-2M先導(dǎo)燃料組件堆芯熱工參數(shù)對比

AFA燃料組件堆芯和裝入6組先導(dǎo)燃料組件堆芯的計算結(jié)果表明：裝入6組TVS-2M先導(dǎo)燃料組件堆芯較全AFA燃料組件堆芯的反應(yīng)堆流量減少0.06%；通過TVS-2M燃料組件的冷卻劑流量較AFA燃料組件減少7%左右[5]。考慮保守的參數(shù)偏差組合，第5～8燃料循環(huán)裝有6組TVS-2M先導(dǎo)組件堆芯的熱工水力參數(shù)計算結(jié)果列于表3。

3.2 試驗測量結(jié)果

為驗證從第5燃料循環(huán)裝入6組TVS-2M先導(dǎo)燃料組件堆芯設(shè)計和運行安全準(zhǔn)則的滿足性，跟蹤測量第5～7燃料循環(huán)的堆芯熱工水力參數(shù)，在換料大修后反應(yīng)堆達到滿功率運行狀態(tài)下完成了熱工水力參數(shù)的試驗測量。測量的參數(shù)包括反應(yīng)堆功率，冷卻劑溫度，冷卻劑流量，以及反應(yīng)堆、一回路、蒸汽發(fā)生器壓力損失等熱工水力參數(shù)，熱工水力測量結(jié)果[6]表明，各項參數(shù)均滿足驗收準(zhǔn)則，詳細(xì)試驗結(jié)果列于表4。

表3 計算結(jié)果

表4 熱工水力測量試驗結(jié)果

續(xù)表4

另外，在機組大修卸料期間完成對燃料組件的變形檢查測量，測量結(jié)果表明，燃料組件變形的形式主要有S形和C形兩種，不同燃料循環(huán)下組件的最大變形量均小于參考限值，TVS-2M組件在同等條件下較AFA組件的變形量小[7]，該測量結(jié)果與俄羅斯核電站TVS-2M組件應(yīng)用實踐基本相符[8-9]。測量結(jié)果列于表5。

表5 燃料組件最大變形量測量結(jié)果

4 結(jié)論

田灣核電站一號機組從第5燃料循環(huán)裝入6組TVS-2M先導(dǎo)燃料組件，反應(yīng)堆穩(wěn)態(tài)熱工水力分析結(jié)果表明，堆芯熱工水力設(shè)計參數(shù)滿足穩(wěn)態(tài)設(shè)計準(zhǔn)則要求；機組運行中堆芯熱工水力參數(shù)實測結(jié)果以及燃料組件變形檢查表明，堆芯各熱工水力參數(shù)滿足試驗驗收準(zhǔn)則要求，TVS-2M燃料組件在同等條件下較AFA燃料組件的變形量小，燃料組件變形量遠(yuǎn)低于參考限值。通過對裝入先導(dǎo)燃料組件堆芯的熱工水力分析及試驗測量，驗證了TVS-2M燃料組件的運行性能及與AFA燃料組件的相容性，保證堆芯在過渡循環(huán)下的運行可靠性和安全性，同時驗證了熱工水力設(shè)計程序計算模型的合理性以及計算結(jié)果與試驗結(jié)果的符合性。

參考文獻：

[1]MEDVEDEV V S. Justification of operation of six pilot TVS-2M starting from the fifth fuel cycles at TNPP, Unit 1[R]. Moscow: OKB Gridropress, 2009.

[2]VOLKOV E, PAVLOVICHEV A, PAVLOV V, et al. Calculation of neutron physics characteristics of transient fuel cycles, beginning from the 8th fuel cycle, with reaching the equilibrium fuel cycle of TNPS, Units 1&2[R]. Moscow: NRC Kurchatov Institute, 2011.

[3]Tianwan NPP. Stationary thermo-hydraulic calculation of the reactor and primary circuit “STAR-1”: Description of computer code, 428-Pr-281[R]. Moscow: OKB Gidropress, 1998.

[4]Tianwan NPP. Calculation of DNBR in fuel rod bundles “PUCHOK-1000”: Description of computer code, 428-Pr-289[R]. Moscow: OKB Gidropress, 1998.

[5]Tianwan NPP. Topical report: Steady state conditions with operation of four, three and two RCP sets 428-Pr-743[R]. Moscow: OKB Gidropress, 2011.

[6]姚進國，冷金珍. 田灣核電廠反應(yīng)堆熱工水力設(shè)計[J]. 核動力運行研究，2005，18(4)：10-13.

YAO Jinguo, LENG Jinzhen. The reactor thermohydraulic design of Tianwan Nuclear Power Station[J]. Nuclear Power Operation Research, 2005, 18(4): 10-13(in Chinese).

[7]姚進國. WWER-1000燃料組件特點及棒彎曲分析[J]. 核動力工程，2006，27(1)：43-46.

YAO Jinguo. Characteristic of fuel assemblies and analysis of fuel rod bowing of WWER-1000[J]. Nuclear Power Engineering, 2006, 27(1): 43-46(in Chinese).

[8]Tianwan NPP. Topical report: Experimental study of TVS-2M, 428-Pr-736[R]. Moscow: OKB Gidropress, 2009.

[9]Tianwan NPP. Topical report: Operational experience of TVS-2M at Russian NPPs, 428-Pr-747[R]. Moscow: OKB Gidropress, 2011.