楊海峰，高 鑫

（中國核電工程有限公司，北京 100840）

核電廠SCIENCE程序包反射層功能擴展研究

楊海峰，高 鑫

（中國核電工程有限公司，北京 100840）

文章根據(jù)實際工作需要，闡述了基于來自法國的SCIENCE V1程序包進行反射層計算的功能擴展開發(fā)；深入研究部分功能的實現(xiàn)細節(jié)，按需增加部分代碼，擴展開發(fā)出反射層計算功能。應用擴展后的SCIENCE程序包研究秦山二期堆型反射層參數(shù)，并與1、2號機組共19個循環(huán)堆芯的實測值進行比較，結(jié)果符合良好，滿足工程精度要求。通過對SCIENCE V1程序包的功能擴展研究，擴大了SCIENCE V1程序包的適用范圍，同時積累了大型程序的研發(fā)經(jīng)驗。

SCIENCE程序包；功能擴展開發(fā)；反射層計算

引進法國的SCIENCE V1程序包，主要包括圖形界面程序COPILOTE[1]和底層計算程序APOLLO2-F[2]、SMART、SQUALE等，廣泛應用于方形組件和M310堆型的計算研究。由于SCIENCE V1程序包自帶的反射層參數(shù)僅適用于M310堆型，且未提供反射層計算功能，因此，SCIENCE V1程序包用于其他堆型時，特別是圍板反射層結(jié)構(gòu)與M310堆型差別較大時，會引入一定的、無法評估的系統(tǒng)誤差，從而限制了其適用范圍。

為了進行秦山二期600 MW堆型的計算以及新堆型的研究，需要基于SCIENCE V1程序包擴展出反射層計算功能，為相應堆型研究配套的反射層參數(shù)。反射層計算功能擴展開發(fā)工作，滿足了當前的工程需要，同時擴大SCIENCE V1程序包的適用范圍，并且積累大型程序的開發(fā)經(jīng)驗。

1 主體構(gòu)架

在SCIENCE程序包的圖形界面下（見圖1），是底層的計算程序，APOLLO2-F、SMART、SQUALE等。這些程序具有相同的體系構(gòu)架，以組件計算程序APOLLO2-F為例：APOLLO2-F是一個基于模塊化設(shè)計的，功能強大、通用的組件計算程序。每個模塊都是一個選項眾多、功能強大的子程序。APOLLO2-F有一個龐大的模塊集，包括生成幾何模塊，生成材料模塊，通量計算模塊，燃耗計算模塊，信息輸出模塊等。這個模塊集是APOLLO2-F的基礎(chǔ)，其所有功能都是通過這些模塊在不同層面上的組合搭配提供。

模塊集中的每個模塊只完成一個專門的、基本的功能，讓用戶基于這些模塊，安排計算流程，準備輸入文件，是一件繁瑣且易出錯的事情，要求用戶非常熟悉物理計算流程，并且熟練掌握APOLLO2-F龐大的模塊集以及每個模塊的各種用法及眾多選項。因此，基于模塊集，針對常用的組件及各種通用計算功能，SCIENCE程序包為APOLLO2-F配備了一套PROC函數(shù)和BLOC塊。用戶不再需要基于模塊集為各種通用功能重復準備輸入文件，只需要調(diào)用相應的PROC/ BLOC即可，減輕用戶負擔，減少了用戶出錯的機會，降低了APOLLO2-F程序的使用難度。

圖1 燃料組件的標準計算流程Fig.1 Standard fuel assembly calculation flow chart

在這些PROC/BLOC的基礎(chǔ)上，對各種日常計算任務進一步歸納總結(jié)，將每種計算任務封裝成一個程序，并為每個程序配備圖形界面、輸入文件模板等。用戶通過圖形界面輸入?yún)?shù)，執(zhí)行計算，查看計算結(jié)果。在SCIENCE程序包中，組件計算主要用到4個程序：Gentiane生成組件模型、Evolution_apollo2計算組件燃耗、Reprise_ apollo2完成多參數(shù)重啟動計算、Creation_ bibliotheque建立多參數(shù)表格化數(shù)據(jù)庫。

基于這套程序集，SCIENCE程序包給出了組件計算的標準流程，如圖1所示，每個方塊代表一個計算步，最終生成堆芯計算所需要的宏觀截面數(shù)據(jù)庫。

SCIENCE程序包的主體框架如圖2所示。

圖2 SCIENCE程序包主體框架示意圖Fig.2 Main framework of SCIENCE code package

2 反射層計算功能開發(fā)

以鐵水混合物為主要成分的壓水堆反射層，通常采用一維模型進行擴散計算。燃料組件模擬堆芯，鐵水混合物模擬堆芯外圍的圍板、吊籃、含硼水等結(jié)構(gòu)，最后形成一個均勻的反射層組件，使用兩群通量不連續(xù)因子ADF1、ADF2，兩群擴散系數(shù)D1、D2，轉(zhuǎn)移截面r∑，兩群吸收截面來表征中子物理特性，用于堆芯計算程序SMART中求解中子擴散方程。

反射層組件的宏觀截面主要受慢化劑密度和可溶硼濃度的影響，因此，在數(shù)據(jù)庫文件TMPN中存放著隨硼濃度和慢化劑密度變化的宏觀截面。硼濃度的典型值為0、600、1 200、1 800 ppm（1 ppm=10-6），慢化劑密度的典型值為0.907 0、0.798 4、0.713 0、0.644 1 g/cm3，形成的數(shù)據(jù)貯存如圖3所示。在堆芯計算時，根據(jù)實際的慢化劑密度及硼濃度，采用二維三點插值的方法，恢復出相應狀態(tài)下的宏觀截面，用于中子擴散方程的求解。

圖3 TMPN文件中數(shù)據(jù)貯存示意圖Fig.3 Sketch of data storage in TMPN

常規(guī)的燃料組件是1/8或1/4對稱的，因此，在進行計算及結(jié)果輸出時，默認1/8或1/4組件兩種模型。在反射層組件計算模型中，在組件的一側(cè)設(shè)置燃料棒柵元，其余柵元為鐵水混合物。燃料棒柵元提供中子源，鐵水混合物柵元組成反射層。1/8或1/4對稱性不復存在，必須考慮全組件的計算。常規(guī)燃料組件及反射層組件模型如圖4所示。

在建立全組件模型時，1/8或1/4組件模型對稱性不復存在，涉及邊水隙和角水隙的處理、等效幾何的產(chǎn)生等，因此需要對相應的PROC函數(shù)進行擴展，增加對全組件模型的支持。通過功能擴展，使得APOLLO2-F支持全組件的計算，在輸出文件中，能夠正確地顯示整個組件所有柵元的功率、通量、反應率等信息。

圖4 常規(guī)組件與反射層組件示例Fig.4 Models of normal fuel assembly and reflector assembly

在APOLLO2-F完成通量計算后，生成整個組件平均的兩群宏觀截面。對于反射層組件，只考慮鐵水混合物柵元，因此需要研究AOLLO2-F生成宏觀截面的實現(xiàn)細節(jié)。通過研究可知，在APOLLO2-F中，針對每個柵元，統(tǒng)計了對能量積分、對空間平均的各種少群宏觀反應率：總反應率、吸收反應率、裂變反應率、增殖（n2n，n3n，…）反應率、輸運反應率、泄漏反應率、轉(zhuǎn)移反應率等。通過對相應的鐵水混合物柵元的少群宏觀反應率進行體積加權(quán)，除以相應柵元通量的體積加權(quán)，可得到反射層組件的兩群宏觀截面。

基于APOLLO2-F組件計算得到的逐柵元的、非均勻化通量，自編二維擴散程序，得到燃料部分和反射層部分均勻化宏觀截面，以此為基礎(chǔ)求解擴散方程，得到界面處的均勻通量，從而最終得到通量不連續(xù)因子ADF1、ADF2。

借助于組件計算的標準程序，Evolution_ apollo2程序可以完成參考狀態(tài)下組件通量計算，Reprise_apollo2程序可完成不同硼濃度和慢化劑密度組合下的組件通量計算。Gentiane程序因為不能輸入多種鐵水混合物材料而不適合布置反射層組件，同時缺少生成反射層組件宏觀截面庫TMPN文件的程序。因此，深入研究了SCIENCE程序包的圖形界面，創(chuàng)建兩個專門用于反射層組件計算的程序：GenRefl程序生成反射層組件的計算模型，ReflCalc程序生成反射層組件的宏觀截面，計算組件不連續(xù)因子ADF1、ADF2，并調(diào)用數(shù)據(jù)庫文件生成程序MKDB生成TMPN文件。

按照SCIENCE程序包的規(guī)則，為新創(chuàng)建的兩個程序GenRefl和ReflCalc配備了輸入文件模板、執(zhí)行腳本模板、圖形界面，將其并入到SCIENCE程序包中，建立反射層計算的標準流程，如圖5所示。

圖5 反射層計算標準流程Fig.5 Standard reflector calculation flow chart

3 驗證與應用

開發(fā)的反射層計算功能，應用于秦山二期1、2號機組的反射層參數(shù)研究。秦山二期1、2號機組圍板實際厚度為22.22 mm，成形板（輔板）的厚度為35 mm，燃料組件和圍板之間的水隙為0.64 mm。圍板和成形板（輔板）材料為控氮不銹鋼，其中鎘和鎳元素質(zhì)量分數(shù)占30%，其余約為鐵元素和微量元素。通過對反射層組件計算過程中參數(shù)的敏感性研究，包括組件形狀、組件尺寸、燃料棒列數(shù)、燃料富集度、水層厚度、水隙等，考慮秦山二期1、2號機組實際的反射層結(jié)構(gòu)，建立600 MW堆型的反射層組件計算模型[3]：31×31方形組件，左側(cè)11列燃料棒，燃料棒與圍板之間1列水隙、2列柵元描述圍板、17列柵元描述水層。利用計算的600 MW堆型反射層，進行秦山二期1號機組1～10個循環(huán)堆芯、2號機組1～9個循環(huán)堆芯共19個循環(huán)堆芯的驗算，并與業(yè)主提供的堆芯實測數(shù)據(jù)進行比較，結(jié)果如下：

針對組件功率大于0.9的所有組件進行統(tǒng)計：63.1%的組件，其功率誤差在±1%以內(nèi)；99.9%的組件，其功率誤差在±3%以內(nèi)；所有的組件，其功率誤差在±5%以內(nèi)，滿足工程精度要求。其誤差分布如圖6所示。

針對組件功率小于0.9的所有組件進行統(tǒng)計：72.17%的組件，其功率誤差在±3%以內(nèi)；95.3%的組件，其功率誤差在±5%以內(nèi)；所有的組件，其功率誤差在±8%以內(nèi)，滿足工程精度要求。其誤差分布如圖7所示。

圖6 組件功率大于0.9的組件功率誤差分布圖Fig.6 Power deviation of fuel assemblies whose power is greater than 0.9

圖7 組件功率小于0.9的組件功率誤差分布圖Fig.7 Power deviation of fuel assemblies whose power is less than 0.9

堆芯臨界硼濃度方面：5 922個實測點的硼濃度偏差均小于50 ppm，滿足工程精度要求，其誤差分布如圖8所示。

由以上19個循環(huán)堆芯的主要計算結(jié)果與實測值的比較可知，基于擴展的SCIENCE程序包研究的600 MW堆型的反射層參數(shù)，堆芯主要計算結(jié)果與實測值符合良好，滿足工程精度的要求。

圖8 硼濃度誤差分布圖Fig.8 Deviation of boron concentration

4 結(jié)論

基于從法國引進的SCIENCE V1程序包，深入研究部分功能的實現(xiàn)細節(jié)與圖形界面的工作原理，創(chuàng)建兩個專門用于反射層計算的程序，搭建反射層計算的標準流程，為SCIENCE V1程序包添加了反射層計算的功能，生成的反射層截面參數(shù)可無縫地應用到堆芯計算程序SMART中，從而擴大了SCIENCE V1程序包的適用范圍。

利用開發(fā)后的SCIENCE程序包，對秦山二期1、2號機組的反射層參數(shù)進行研究，得到秦山二期1、2號機組的反射層參數(shù)，并對秦山二期1號機組1～10個循環(huán)堆芯、2號機組1～9個循環(huán)堆芯進行跟蹤計算，將主要計算結(jié)果與實測結(jié)果進行對比，符合良好，證明擴展后的SCIENCE程序包可很好地應用于其他堆型如秦山二期600 MW堆型及新堆型等的反射層參數(shù)的計算研究。

[1] T.MARTIN，COPILOTE USER MANUAL，EPDS/DC.0011 Rev.A，1998，F(xiàn)RAMATOME.

[2] A.SARGENI，APOLLO2-F USER MANUAL，EPDN/DC.078 Rev.A，1998，F(xiàn)RAMATOME.

[3] 高鑫，楊海峰. 基于SCIENCE程序的反射層參數(shù)計算. 原子能科學技術(shù)，2012，12（46）.（GAO Xin, YANG Hai-feng. Calculation of reflecting layer parameters based on SCIENCE program. Atomic Energy Science and Technology, 2012, 12(46).）

Study on Ref ector Function Expansion for SCIENCE Code Package

YANG Hai-feng，GAO Xin
(China Nuclear Power Engineering Co.，Ltd.，Beijing 100840，China)

Based on SCIENCE V1 code package imported from France, function expansion for reflector calculation is developed according to the requirements. The function modules beneath the graphical user interface are studied, and then the function of reflector calculation is developed. After the development, SCIENCE code package is used to investigate the reflector parameters for Unit 1 & 2 Qinshan II. Compared with the measurements from the total 19 fuel cycles of Unit 1 & 2, the main theoretical calculation results fit well with the measured ones. Benefit from this study, reflector calculation function is added to SCIENCE V1 code package, the scope usage of SCIENCE code package is expanded, and some R&D experience of large scale code package is obtained.

SCIENCE code package; development of function expansion; reflector calculation

TL37 Article character：A Article ID：1674-1617(2014)03-0201-06

TL37

A

1674-1617(2014)03-0201-06

2014-06-17

楊海峰（1981—），男，河南南陽人，高級工程師，碩士研究生，從事反應堆物理與臨界安全分析工作。