AP1000核反應(yīng)堆控制棒價(jià)值特性的MC模擬

謝明亮，于　雷，陳玉清

(海軍工程大學(xué) 核能科學(xué)與工程系，武漢　430033)

?

AP1000核反應(yīng)堆控制棒價(jià)值特性的MC模擬

謝明亮，于雷，陳玉清

(海軍工程大學(xué) 核能科學(xué)與工程系，武漢430033)

摘要：針對(duì)當(dāng)前AP000堆芯采用的兩類控制棒束，基于MCNP5程序建立堆芯仿真計(jì)算模型，分析了含不同硼濃度對(duì)堆芯k(eff)與硼微分價(jià)值的影響，同時(shí)對(duì)AP1000棒組價(jià)值進(jìn)行模擬計(jì)算，對(duì)比分析了黑棒與灰棒插入堆芯對(duì)k(eff)的影響。結(jié)果表明：基于MCNP5程序建立的模型是正確的，硼微分價(jià)值(絕對(duì)值)隨硼濃度增加呈現(xiàn)下降趨勢(shì)，其值在-9.16 ～ -13.60范圍內(nèi)變化，符合反應(yīng)性設(shè)計(jì)要求，有效增殖系數(shù)k(eff)隨控制棒插入呈現(xiàn)非線性變化，得到了控制棒的價(jià)值變化曲線與擬合關(guān)系式，為控制棒在反應(yīng)堆內(nèi)緊急控制與功率調(diào)節(jié)提供參考。

關(guān)鍵詞：AP1000;控制棒價(jià)值;MCNP;灰棒

Citation format:XIE Ming-liang, YU Lei, CHEN Yu-qing.Analysis of Control Rod Worth Features of MC Simulation in AP1000 Reactor[J].Journal of Ordnance Equipment Engineering,2016(3):121-125.

AP1000是美國(guó)西屋公司開發(fā)的大型非能動(dòng)先進(jìn)壓水堆，其堆芯設(shè)計(jì)與在役壓水堆核電廠有較大的差別，控制棒組件由黑棒組件(RCCA)和灰棒組件(GRCA)組成[1]。黑棒控制組件主要吸收體材料為銀-銦-鎘合金，可補(bǔ)償快速的反應(yīng)性變化，控制軸向功率分布[2]，灰棒控制組件可為負(fù)荷跟隨控制提供反應(yīng)性補(bǔ)償，最大程度地減少可溶硼濃度的變化需求[3]，AP1000灰棒組件最初設(shè)計(jì)由12根不銹鋼棒與12根銀-銦-鎘棒組成[4]，國(guó)內(nèi)擬將吸收材料改為鎢和Inconel718 套管，為校核設(shè)計(jì)結(jié)果，系統(tǒng)掌握AP1000控制棒組價(jià)值曲線，本研究基于MCNP5程序建立堆芯仿真計(jì)算模型，分析堆芯含不同硼濃度對(duì)keff與硼微分價(jià)值的影響，對(duì)AP1000棒組價(jià)值進(jìn)行模擬計(jì)算，對(duì)比分析了黑棒與灰棒插入堆芯對(duì)keff的影響，并得到相應(yīng)的價(jià)值曲線與擬合關(guān)系式。

1計(jì)算模型建立

1.1堆芯幾何結(jié)構(gòu)

核反應(yīng)堆堆芯是由一定數(shù)量的燃料組件和相關(guān)組件構(gòu)成，燃料組件主要由燃料棒、上下管座、格架、導(dǎo)向管和儀表管等組成，相關(guān)組件主要包括控制棒組件、可燃毒物組件、中子源與阻力塞組件等[5]。AP1000反應(yīng)堆堆芯活性區(qū)高度(cm)為426.72，堆芯等效直徑(cm)為304.04[6]，冷態(tài)H2O/U為2.40，堆芯第一循環(huán)燃料采用不同的燃料富集度，根據(jù)燃料的富集度不同，堆芯劃分為5個(gè)燃料區(qū)，分別為第一區(qū)到第五區(qū)，同時(shí)在第三、四、五區(qū)軸向方向上，燃料棒的兩端還設(shè)有再生區(qū)，以減少中子的泄露，增加中子利用率[5]，其富集度分別為0.74/--、1.58/--、3.20/1.58、3.776/3.20、4.376/3.20。

堆芯含有157個(gè)按17 cm×17 cm XL Robust燃料(改進(jìn)型)正方形柵格排列的組件[7]，每個(gè)組件包含264根燃料棒，中心位置設(shè)置1根儀表管用于堆內(nèi)測(cè)量，還有24個(gè)位置設(shè)置導(dǎo)向管，這些導(dǎo)向管與格架、上管座和下管座相連接，形成燃料組件骨架，為燃料棒提供支承結(jié)構(gòu)[8]，燃料組件參數(shù)如表1所示。

表1　燃料組件參數(shù)

1.2控制棒的布置

燃料組件的導(dǎo)向管內(nèi)可以放置黑控制棒組件(RCCA)、灰控制棒組件(GRCA)、中子源組件、通水環(huán)狀可燃毒物組件(WABA)或阻力塞組件[6]。AP1000堆芯內(nèi)共布置53束黑棒，其吸收材料為Ag(80%)-In(15%)-Cd(5%)，根據(jù)控制棒用途差異主要分為機(jī)械補(bǔ)償控制黑棒組M1、M2，軸向功率偏移控制黑棒組AO以及停堆棒組S1-S4[9]，M2與S4控制棒組布置在燃料富集度第三區(qū)，其余黑棒布置在第二區(qū)內(nèi)，灰棒的吸收材料主要為鎢(Wu)+套管718，棒組數(shù)量只有16束，其用途為機(jī)械補(bǔ)償控制灰棒組MA-MD，除MB布置在第五區(qū)內(nèi)含124IFBA毒物組件的導(dǎo)向管內(nèi)外，其余灰棒控制組件都布置在第二區(qū)內(nèi)，表2給出含控制棒組件的基本參數(shù)，圖1給出控制棒在堆芯中的布置。

表2　控制棒(黑棒/灰棒)參數(shù)

圖1　RCCA與GRCA在堆芯中的布置

1.3模型建立

根據(jù)AP1000堆芯的幾何結(jié)構(gòu)和控制棒布置情況，基于MCNP5程序建立堆芯仿真計(jì)算模型，由于AP1000堆芯組件個(gè)數(shù)及組件類型比較多，包含兩端含再生區(qū)的組件類型一共有9種，推到柵元的類型就更多了，且柵元、組件具有重復(fù)性，利用MCNP建模中的重復(fù)結(jié)構(gòu)功能，對(duì)幾何中多次出現(xiàn)的結(jié)構(gòu)柵元僅用一次描述，重復(fù)結(jié)構(gòu)主要包含以下4個(gè)柵元參數(shù)卡[3]：Universe卡、TRCL卡、LAT卡、FILL卡。

根據(jù)堆芯幾何結(jié)構(gòu)模塊建立的一般方法，堆芯結(jié)構(gòu)采用Universe卡描述[10]，按照實(shí)際幾何尺寸建立堆芯的幾何模型，圓柱形堆芯半徑為169.86 cm，活性區(qū)冷態(tài)高度為426.72 cm，把這一模型規(guī)定為第1級(jí)，這一級(jí)是真空，沒有材料，只是一個(gè)幾何模塊，等待被填充，采用LAT柵格卡創(chuàng)建被填充無邊界柵格(Universe1)，幾何尺寸根據(jù)實(shí)際燃料組件間距(21.402 cm×21.402 cm)設(shè)定，并等待下一級(jí)燃料組件或相關(guān)組件填充，建立燃料棒柵元空間(Universe2)，柵距為燃料棒間的距離1.26 cm，這一級(jí)模塊依然是沒有材料的幾何模塊，等待被燃料棒等填充，最后將燃料棒、控制棒、可燃毒物棒等柵元描述模塊填入卡中，同時(shí)考慮組件間的水隙，然后再按從里到外的順序構(gòu)造反應(yīng)堆其他部件，堆芯外圍描述還包含吊籃、反射層水和壓力容器，這樣就基本完成堆芯模塊的建立。再將燃料、毒物、慢化劑、包殼等材料描述卡與源的描述卡填入輸入卡中，經(jīng)過上述步驟就初步形成了反應(yīng)堆仿真計(jì)算模型，圖2給出WABA毒物柵元、組件、圍板/反射層布置圖，組件中有燃料柵元，24根控制棒束，同時(shí)可以看到圍板、吊籃、反射層水在堆芯內(nèi)的基本布置。

圖2　WABA柵元、組件、圍板/反射層布置

由于AP1000壓力容器堆內(nèi)構(gòu)件非常復(fù)雜，在對(duì)構(gòu)件建模的過程中，省略了部分對(duì)堆芯計(jì)算結(jié)果影響不大的堆內(nèi)構(gòu)件，堆芯被外層圍板包圍，整個(gè)結(jié)構(gòu)都置于堆芯吊籃中，其余部分被冷卻水填充，圖3給出堆芯剖面示意圖。

圖3　堆芯剖面(堆芯X-Y剖面,堆芯X-Z剖面)

2MC模擬計(jì)算及結(jié)果分析

2.1典型硼濃度驗(yàn)證及硼微分價(jià)值計(jì)算

為了驗(yàn)證模型的準(zhǔn)確性，并開展控制棒價(jià)值特性仿真，基于所建的臨界計(jì)算模型，選取反應(yīng)堆冷態(tài)(20℃)零功率與熱態(tài)(291.7℃)零功率兩種運(yùn)行工況，控制棒全提出情況下，計(jì)算了典型硼溶度下的堆芯有效增殖系數(shù)keff，并與AP1000安全分析報(bào)告參考值進(jìn)行比對(duì)，該參考值由兩群擴(kuò)散理論節(jié)塊模型得到，MCNP計(jì)算結(jié)果記錄如表3。四組典型狀態(tài)的計(jì)算結(jié)果表明，MCNP計(jì)算值與參考值絕對(duì)誤差控制在5‰以內(nèi)，初步驗(yàn)證模型建立的正確性。

表3　典型硼濃度下堆芯keff

圖4給出了不同硼濃度下，堆芯keff與硼微分價(jià)值的變化情況，可以看出硼濃度較低時(shí)，堆芯處于冷態(tài)與熱態(tài)兩工況下，對(duì)堆芯keff有較大的差別，但隨堆芯硼濃度的增加，兩者差別逐漸較少，即隨著堆芯硼濃度的增加，慢化劑溫度系數(shù)絕對(duì)值逐漸減低，在硼濃度為1 350 ppm附近，冷熱工況所對(duì)應(yīng)的堆芯keff值相等，當(dāng)硼濃度大于1 350 ppm時(shí)，將出現(xiàn)正的慢化劑溫度系數(shù)，這在堆芯設(shè)計(jì)中是不允許的，同時(shí)硼微分價(jià)值的絕對(duì)值呈現(xiàn)出下降的趨勢(shì)，熱態(tài)工況下的硼微分價(jià)值(絕對(duì)值)低于冷態(tài)工況對(duì)應(yīng)值，其值在-9.16～-13.60 變化，符合反應(yīng)堆設(shè)計(jì)要求，保證了化學(xué)補(bǔ)償控制的安全裕量。

圖4　堆芯keff與硼微分價(jià)值隨硼濃度的變化情況

2.2AP1000棒組價(jià)值模擬計(jì)算

根據(jù)控制棒控制用途的不同，可將控制棒分為多種類型。不同類型控制棒組，有不同的反應(yīng)性當(dāng)量，棒組價(jià)值直接表征了該控制棒對(duì)反應(yīng)性抵消能力[11]，為此建立插入黑棒(Ag-In-Cd)和灰棒(Wu)堆芯仿真模型，并選取熱態(tài)(291.7℃)工況下，RCCAs提出的硼濃度為1 379 ppm作為臨界基準(zhǔn)，不改變其他控制棒狀態(tài)，將單組控制棒組插入堆芯。表4給出了在熱態(tài)零功率(HZP)，循環(huán)初(BOL)、無氙無重疊下單組控制棒組插入時(shí)，不同棒組類型下有效增殖系數(shù)keff的變化情況與對(duì)應(yīng)的棒組價(jià)值。

表4　棒組價(jià)值隨棒組類型變化情況

可以看出MCNP計(jì)算值與參考值吻合比較好，同時(shí)表明MA、MB、MC、MD控制棒組價(jià)值相當(dāng)，但整體低于M1、AO棒組價(jià)值，AO棒組具有最大的反應(yīng)性控制當(dāng)量。

2.3控制棒價(jià)值特性分析

在反應(yīng)堆堆芯設(shè)計(jì)中不僅要知道控制棒全部插入時(shí)的價(jià)值，還需要知道控制棒在堆芯中插入不同深度對(duì)反應(yīng)性的影響[12]，為此建立對(duì)應(yīng)的堆芯計(jì)算模型，為了單獨(dú)分析黑棒與灰棒對(duì)堆芯keff的影響，計(jì)算模型不考慮可燃毒物布置與化學(xué)補(bǔ)償?shù)挠绊懀捎诙研臼籽h(huán)燃料活性區(qū)高度為426.72 cm，為了更準(zhǔn)確詳細(xì)地描述控制棒在堆芯移動(dòng)情況，將堆芯軸向分成20層(軸向坐標(biāo)分別為213.36、193.05、169.35、145.65、121.95、104.15、86.4、64.8、43.2、21.6，0、-21.6、-43.2、-64.8、-86.4、-108.0、-129.6、-154.95、-172.75、-193.05、-213.36)。

黑棒(Ag-In-Cd)具有很高的吸收棒價(jià)值。圖5給出控制黑棒插入堆芯，對(duì)堆芯有效增殖系數(shù)keff的影響情況。灰棒主要為負(fù)荷跟隨控制提供反應(yīng)性補(bǔ)償?shù)墓δ埽刂品磻?yīng)性當(dāng)量稍小于黑棒，但能抵消部分后備反應(yīng)性以減少對(duì)硼濃度的變化需求。圖6給出兩種灰棒組件類型(12根不銹鋼+12根Ag-In-Cd和24根Wu棒)隨插入堆芯深度對(duì)有效增殖系數(shù)keff的影響情況。

圖5　keff隨黑棒(Ag-In-Cd)插入深度變化情況

可以看出，有效增殖系數(shù)keff隨控制棒插入深度呈現(xiàn)非線性變化，控制棒插入堆芯初始階段對(duì)keff影響較緩慢，但隨著堆芯深度逐漸增加，對(duì)keff影響不斷變大，最后都趨于平緩?？刂坪诎魧?duì)有效增殖系數(shù)keff控制在1.33～1.24之間變化，灰棒的吸收價(jià)值低于黑棒，含不銹鋼+Ag-In-Cd灰棒在控制棒插入初始階段，對(duì)keff影響與含鎢(Wu)灰棒相當(dāng)，但隨著插入堆芯深度逐漸增加，對(duì)keff影響大于含鎢(Wu)灰棒，并控制在1.33～1.31之間變化。

圖6　keff隨灰棒(鎢Wu、不銹鋼+Ag-In-Cd)

由于控制棒價(jià)值隨插入堆芯不同深度而改變，常用控制棒的積分價(jià)值和微分價(jià)值表征[11]。圖7給出黑棒的微分價(jià)值和積分價(jià)值隨控制棒插入深度的變化情況，可以看出在控制黑棒剛插入堆芯時(shí)，控制棒的價(jià)值變化較小，并與控制棒移動(dòng)距離呈現(xiàn)一種非線性關(guān)系。在控制棒插入到300～400 cm 時(shí)，控制棒微分價(jià)值逐漸增大并出現(xiàn)一個(gè)較大的峰值，隨后趨于平緩，歸一化黑棒積分價(jià)值與插入堆芯深度，可以得到黑棒(Ag-In-Cd)的相對(duì)價(jià)值與插入相對(duì)深度的擬合關(guān)系式為

0.338 9×exp(-((h-0.914 9)/0.066)2)

其中:h為控制棒插入相對(duì)深度;ρ(h)為控制棒相對(duì)價(jià)值，擬合R-square=0.999。

圖7　黑棒價(jià)值曲線隨插入深度的變化情況

圖8給出灰棒(12根不銹鋼+12根Ag-In-Cd與24根Wu棒)的微分價(jià)值和積分價(jià)值隨控制棒插入深度的變化情況，可以看出灰棒價(jià)值曲線與黑棒相似，當(dāng)控制棒剛插入堆芯與插入底部附近時(shí)，對(duì)控制棒的價(jià)值變化較小，反應(yīng)性變化不靈敏，但在控制棒插入到250～350 cm時(shí)，出現(xiàn)較大的反應(yīng)性變化率，微分價(jià)值曲線呈現(xiàn)峰值，且含鎢(Wu)灰棒微分價(jià)值峰值低于含不銹鋼+Ag-In-Cd灰棒，因此在基負(fù)荷運(yùn)行和負(fù)荷跟隨控制可以更好補(bǔ)償反應(yīng)性的微小變化，同時(shí)在換料周期末，采用灰棒對(duì)反應(yīng)性進(jìn)行微調(diào)，以提高對(duì)反應(yīng)性控制的精準(zhǔn)度。歸一化灰棒積分價(jià)值與插入堆芯深度，可以得到灰棒(12根Ag-In-Cd與24根Wu棒)與灰棒(24根Wu棒)的相對(duì)價(jià)值與插入相對(duì)深度的擬合關(guān)系式，分別為式(1)與式(2)

(1)

ρ(h)=1.021×exp(-((h-0.980 7)/0.358 5)2)

0≤h≤1, 0≤ρ(h)≤1R-square:0.99

(2)

其中:h為控制棒插入相對(duì)深度;ρ(h)為控制棒相對(duì)價(jià)值。

圖8　灰棒價(jià)值曲線隨插入深度的變化情況

3結(jié)論

針對(duì)當(dāng)前AP000堆芯采用的兩類控制棒，基于MCNP5程序建立堆芯仿真計(jì)算模型，分析堆芯含不同硼濃度對(duì)keff與硼微分價(jià)值的影響，同時(shí)對(duì)AP1000棒組價(jià)值進(jìn)行模擬計(jì)算，對(duì)比分析了控制黑棒與灰棒插入堆芯對(duì)keff的影響，并得到相應(yīng)的價(jià)值曲線。結(jié)果表明：① 典型硼溶度下，基于MCNP5程序的計(jì)算值與參考值吻合很好，初步驗(yàn)證了建立模型的正確性，同時(shí)得到AP1000典型棒組價(jià)值模擬計(jì)算值。② 隨堆芯硼濃度的增加，反應(yīng)堆處于冷熱不同工況下對(duì)keff的影響逐漸減少，當(dāng)硼濃度大于1 350 ppm時(shí)，將出現(xiàn)正的慢化劑溫度系數(shù)，同時(shí)硼微分價(jià)值(絕對(duì)值)呈現(xiàn)出下降的趨勢(shì)，其值在-9.16～-13.60范圍內(nèi)變化，符合反應(yīng)性設(shè)計(jì)要求。③ 有效增殖系數(shù)keff隨控制棒插入呈現(xiàn)非線性變化，控制棒插入堆芯初始階段對(duì)keff影響較緩慢，隨著堆芯深度逐漸增加，對(duì)keff影響不斷變大，最后都趨于平緩，控制黑棒對(duì)有效增殖系數(shù)keff控制在1.33～1.24變化，控制灰棒(鎢)控制在1.33～1.31變化，并得到控制棒價(jià)值變化曲線與擬合關(guān)系式，為控制棒在反應(yīng)堆內(nèi)緊急控制與功率調(diào)節(jié)提供參考。

參考文獻(xiàn)：

[1]魏強(qiáng)林,劉義保,楊波,等.控制黑棒和灰棒對(duì)AP1000反應(yīng)堆keff值影響的MC模擬[J].科學(xué)技術(shù)與工程,2013,13(8):2074-2078．

[2]彭小強(qiáng),鄭養(yǎng)波.AP1000 控制棒棒位的數(shù)字化測(cè)量技術(shù)分析[J].核電子學(xué)與探測(cè)技術(shù),2012,32(5):627-630.

[3]陳德鋒.基于MCNP程序的AP1000反應(yīng)堆屏蔽分析研究[M].哈爾濱:哈爾濱工程大學(xué),2013.

[4]李曉靜,陳義學(xué),張斌,等.AP1000灰棒燃料組件物理特性分析[J].原子能科學(xué)技術(shù),2013,12(47):614-618.

[5]王興華,周四春,張慶賢,等.MCNP程序在反應(yīng)堆可燃毒物控制棒計(jì)算中的應(yīng)用[J].核安全,2013,12(3):53-56.

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

[7]曹博,陳義學(xué),李曉靜,等.AP1000典型燃料組件物理特性初步分析[J].原子能科學(xué)技術(shù),2013,12(47):599-602.

[8]趙夢(mèng)云,楊波,劉義保,等.MCNP在計(jì)算硼酸質(zhì)量分?jǐn)?shù)影響反應(yīng)堆有效增殖系數(shù)keff值中的應(yīng)用[J].能源研究與管理,2011,4(4):23-26.

[9]WANG P F,WAN J S,CHEN Z,et al.Dynamic simulation and study of Mechanical Shim(MSHIM) core control strategy for AP1000 reactor [J].Annals of Nuclear Energy,2014,72(2014):49-62.

[10]鐘兆鵬,施工,胡永明.MCNP程序在反應(yīng)堆臨界計(jì)算中的應(yīng)用[J].核動(dòng)力工程,2003,24(1):8-11.

[11]FAUSTO F,ZHANG B C,LARRY M H,et al.Development of a control rod depletion methodology for the Westinghouse NEXUS system[J].Progress in Nuclear Energy,2013,68(2013):235-242.

[12]張法邦.核反應(yīng)堆運(yùn)行物理[M].北京:原子能出版社,2000.

(責(zé)任編輯楊繼森)

Analysis of Control Rod Worth Features of MC Simulation in AP1000 Reactor

XIE Ming-liang, YU Lei, CHEN Yu-qing

(Department of Nuclear Energy Science and Engineering,Naval University of Engineering, Wuhan 430033, China)

Abstract:In view of current use of two kinds of control rods in AP1000 reactor, the simulation calculation model was set up based on MCNP5 code, k(eff) and differential value of boron with several different boron concentration were analyzed in the reactor, and simulation of group value rods of AP1000 were calculated at the same time, and the contrast analysis of k(eff) effects was accomplished when black and grey rod inserted into the core. Results show that: the model based on MCNP5 code is correct, and the differential value of boron (absolute value) present a downward trend with the increase of boron concentration, and its value changed between 9.16 and 13.60, which is conforms to the requirements of reactive design, and the effective multiplication factor k(eff) present nonlinear variation with rod inserted into the core, and it gets the curves of control rod worth and fitting relation, which provides the reference for emergency control and power regulation with control rod in the reactor.

Key words:AP1000; control rod; MCNP; gray rod

文章編號(hào)：1006-0707(2016)03-0121-05

中圖分類號(hào)：TL315.5

文獻(xiàn)標(biāo)識(shí)碼：A

doi:10.11809/scbgxb2016.03.029

作者簡(jiǎn)介：謝明亮(1990—)，男，碩士研究生，主要從事核反應(yīng)堆安全分析研究。

收稿日期：2015-09-11；修回日期：2015-09-28

本文引用格式：謝明亮，于雷，陳玉清.AP1000核反應(yīng)堆控制棒價(jià)值特性的MC模擬[J].兵器裝備工程學(xué)報(bào)，2016(3):121-125.

【化學(xué)工程與材料科學(xué)】