熱堆中釷鈾轉(zhuǎn)化規(guī)律

張海青 林 俊 曹長(zhǎng)青 朱天寶 朱智勇

（中國(guó)科學(xué)院上海應(yīng)用物理研究所 嘉定園區(qū) 上海 201800）

熱堆中釷鈾轉(zhuǎn)化規(guī)律

張海青 林 俊 曹長(zhǎng)青 朱天寶 朱智勇

（中國(guó)科學(xué)院上海應(yīng)用物理研究所 嘉定園區(qū) 上海 201800）

釷鈾燃料循環(huán)以其優(yōu)異的物理和化學(xué)特性，受到核能界的廣泛關(guān)注。本文利用單群的點(diǎn)燃耗計(jì)算程序ORIGEN，分別研究了釷燃料在沸水堆(Boiling Water Reactor, BWR)、壓水堆(Pressurized Water Reactor, PWR)和加拿大重水鈾反應(yīng)堆(Canada Deuterium Oxide Uranium, CANDU，又稱坎杜堆)能譜中輻照時(shí)，232Th、233Th、233Pa、233U等核素生成量隨中子注量率和中子能譜的變化規(guī)律，并探索了多次“輻照-冷卻”循環(huán)對(duì)釷鈾轉(zhuǎn)化率的影響。計(jì)算結(jié)果表明，能譜相同時(shí)，233Th和233Pa存量的最大值與注量率有關(guān)；233U存量的最大值與注量率無(wú)關(guān)，大概在注量（注量率×?xí)r間）為4×1016n·cm?2左右；注量率相同時(shí)，能譜越硬，233U存量的最大值越大。采取循環(huán)“輻照-冷卻”可以提高233Th-233U的轉(zhuǎn)化率，對(duì)于相同的總輻照時(shí)間，每次循環(huán)周期內(nèi)的輻照時(shí)間越短，相對(duì)于總輻照時(shí)間相同的單次輻照，轉(zhuǎn)化率增量提高越明顯；當(dāng)總輻照時(shí)間超過(guò)兩個(gè)月時(shí)，循環(huán)輻照對(duì)轉(zhuǎn)化率增量的作用較小，與單次輻照不冷卻相比，轉(zhuǎn)化率相對(duì)增量不超過(guò)1倍。

中子注量率，中子能譜，釷鈾轉(zhuǎn)化率

232Th的天然同位素豐度接近100%，它與238U一樣，都是可增殖核材料。近年來(lái)，隨著鈾資源消耗量的急劇增加，釷鈾燃料循環(huán)研究越來(lái)越為世界核能界所關(guān)注。我國(guó)釷資源相對(duì)豐富，據(jù)統(tǒng)計(jì)，國(guó)內(nèi)釷資源約為鈾資源的6倍[1]。所以，中國(guó)科學(xué)院上海應(yīng)用物理研究所啟動(dòng)了“未來(lái)先進(jìn)核裂變能-釷基熔鹽核能系統(tǒng)”(Thorium Molten Salt Reactor, TMSR)固態(tài)燃料釷基熔鹽堆的建設(shè)。該反應(yīng)堆將采用高溫氣冷堆的球形燃料元件，堆芯由釷球和鈾球混合組成。

近年來(lái)世界各國(guó)很重視釷的利用[2]。美國(guó)、法國(guó)、德國(guó)、加拿大等都開(kāi)展了不少研究工作，取得一定進(jìn)展，其中包括釷在動(dòng)力堆、非傳統(tǒng)動(dòng)力堆[3?4]中的應(yīng)用，釷的應(yīng)用方法[5?6]和燃料循環(huán)技術(shù)等[7?9]。但都是在特定能譜下進(jìn)行的研究。由于中子能譜、中子注量率不同，釷鈾轉(zhuǎn)化情況就不同。本文對(duì)二氧化釷燃料球在燃燒過(guò)程中相關(guān)核素的釷鈾轉(zhuǎn)化規(guī)律進(jìn)行了研究，包括不同熱堆模型下，中子能譜、中子注量率、輻照方式等對(duì)釷鈾轉(zhuǎn)化過(guò)程中相關(guān)核素生成規(guī)律的影響，為將來(lái)釷球燃料元件在堆中的放置方法、增殖周期、添料換料和后處理周期的確定提供了參考數(shù)據(jù)。

1 計(jì)算工具和模型

計(jì)算工具為ORIGEN[10]，它是由美國(guó)橡樹(shù)嶺國(guó)家實(shí)驗(yàn)室(Oak Ridge National Laboratory, ORNL)研發(fā)的一個(gè)多功能點(diǎn)燃耗及放射性衰變計(jì)算機(jī)程序。該程序能模擬核燃料循環(huán)過(guò)程中燃料與中子反應(yīng)后，各放射性物質(zhì)的積累、衰變，核素或元素的添加和元素的連續(xù)移除等各種處理過(guò)程，其中包含1119種獨(dú)立裂變產(chǎn)物，129種錒系核素，698種結(jié)構(gòu)材料或活化材料數(shù)據(jù)，包含衰變庫(kù)和光子數(shù)據(jù)庫(kù)。本文所選的數(shù)據(jù)庫(kù)分別是ORIGEN中的能量為0.0253 eV的熱中子數(shù)據(jù)庫(kù)(THERMAL.LIB)、壓水堆數(shù)據(jù)庫(kù)(PWRU.LIB)、沸水堆(BWRU.LIB)和重水堆(CANDUSEU.LIB)核素中子反應(yīng)截面數(shù)據(jù)庫(kù)。其中PWRU和BWRU的燃料為低富集度鈾，卸載燃耗為33.0 MWd?kgHM?1和27.5 MWd?kgHM?1，CANDUSEU的燃料為1.2 wt%低富集度濃縮鈾(Slightly Enriched Uranium, SEU)。

計(jì)算模型是直徑6 cm的釷燃料球[11?12]（含釷量為6.5 g，232Th富集度為100%），該燃料球由燃料區(qū)和非燃料區(qū)組成，燃料區(qū)半徑2.5 cm，是由釷燃料的三結(jié)構(gòu)同向性型(Tri-structural Iso-tropic, TRISO)包覆顆粒均勻分布在石墨中組成的；非燃料區(qū)的內(nèi)、外半徑分別為2.5 cm和3 cm，材料為石墨。TRISO包覆顆粒由核心和4層包覆層組成，其幾何參數(shù)和成分如表1所示。

表1 TRISO包覆燃料顆粒設(shè)計(jì)參數(shù)Table 1 Design parameters of TRISO coated fuel particles.

2 結(jié)果與討論

2.1 中子注量率對(duì)釷鈾轉(zhuǎn)化相關(guān)核素存量的影響

為研究注量率對(duì)釷鈾轉(zhuǎn)化的影響，本文選定4種注量率：1×1012n·cm?2·s?1、1×1013n·cm?2·s?1、5×1013n·cm?2·s?1、1×1014n·cm?2·s?1；運(yùn)用同一截面庫(kù)THERMAL.LIB分別計(jì)算了在這4種不同注量率下釷球燃料元件中232Th、233Th、233Pa、233U存量的變化情況。

圖1給出了232Th的核反應(yīng)鏈，其反應(yīng)主要為(n,γ)、(n,2n)、α衰變。其貝特曼方程可寫(xiě)為：

在熱中子輻照下，(n,2n)反應(yīng)可忽略不計(jì)，由于232Th的半衰期非常長(zhǎng)，所以α衰變也可忽略不計(jì)。最后解方程(1)，可得232Th隨時(shí)間的變化關(guān)系式(2)。即232Th隨時(shí)間和注量率按指數(shù)關(guān)系下降。圖2為ORIGEN計(jì)算的232Th質(zhì)量隨時(shí)間的變化情況。

圖1 232Th的核反應(yīng)鏈Fig.1 Nuclear chain reaction of 232Th.

圖2 不同注量率下232Th質(zhì)量隨時(shí)間的變化Fig.2 Inventory of 232Th vs. time at different neutron fluence rates.

圖3為233Th的存量占最初232Th (232Th0)摩爾數(shù)的份額隨時(shí)間的變化。由圖3可知，233Th在120 min左右達(dá)到最大值，其最大值主要決定于中子注量率的大小。而達(dá)到最大值的時(shí)間與注量率大小無(wú)關(guān)。這主要是因?yàn)椋?33Th的生成主要由232Th的(n,γ)反應(yīng)而來(lái)，所以注量率大時(shí)生成量就大。其消耗主要為自身的β衰變，與注量率大小無(wú)關(guān)，所以注量率大時(shí)，峰值就大。但半衰期非常短，僅為22.3 min，累積質(zhì)量受時(shí)間的影響很小。所以不同注量率下達(dá)到最大值的時(shí)間基本相同。達(dá)到最大值后，隨著232Th的消耗，峰值也緩慢地減小。

圖4是233Pa的存量占最初232Th(232Th0)摩爾數(shù)的份額隨時(shí)間的變化。由圖4可知，233Pa與233Th的變化趨勢(shì)相同，這也是因?yàn)槠浒胨テ诒容^短的原因(T1/2=26.95 d)。當(dāng)中子注量率大時(shí)，233Pa峰值也高。由圖1可見(jiàn)，233Pa主要來(lái)自233Th的β衰變，233Th含量大時(shí)，233Pa生成量也較多。對(duì)于消耗部分，由于233Pa的自身β衰變比其(n,γ)反應(yīng)截面要大的多，所以233Pa的消耗主要決定于自身衰變，其與注量率大小也基本無(wú)關(guān)。

圖3 不同注量率下233ThInv/232Th0的變化Fig.3 233ThInv/232Th0 vs. time at different neutron fluence rates.

圖4 不同注量率下233PaInv/232Th0的變化Fig.4 233PaInv/232Th0 vs. time at different neutron fluence rates.

由圖5可見(jiàn)，注量率越高時(shí)，233U達(dá)到飽和的時(shí)間越短。但是同一能譜下，233U峰值大小受注量率大小影響很小。為證明這個(gè)結(jié)論，本文專(zhuān)門(mén)做了在與圖5相同的4種注量率下，中子注量相同時(shí)233U的生成情況，如圖6所示，可見(jiàn)不同中子注量率時(shí)233U生成量幾乎在同一注量下達(dá)到峰值。

圖5 不同注量率下233UInv/232Th0的變化Fig.5 233UInv/232Th0 vs. time at different neutron fluence rates.

圖6 相同輻照注量下233UInv/232Th0的變化Fig.6 233UInv/232Th0 vs. time at the same fluencies.

綜上所述，在反應(yīng)堆中利用釷，233Th和233Pa生成量的最大值與注量率有關(guān)，而233U的最大值與注量率無(wú)關(guān)，大概在4×1016n·cm?2左右。這也是想要提取最大質(zhì)量233U時(shí)的熱中子注量。

2.2 中子能譜對(duì)釷鈾轉(zhuǎn)化相關(guān)核素存量的影響

在中子注量率4×1014n·cm?2·s?1下，采用不同中子截面數(shù)據(jù)庫(kù)，計(jì)算233Pa和233U的存量占初始232Th(232Th0)摩爾數(shù)的份額隨時(shí)間的變化，結(jié)果如圖7所示。

在不同能譜下，233Pa生成量與233U生成規(guī)律不同。233Pa在全熱中子譜下，所達(dá)到的飽和值最大；而能譜較硬的壓水堆和沸水堆，飽和值較小。不同數(shù)據(jù)庫(kù)中，232Th的中子俘獲反應(yīng)截面不同，而233Th半衰期只有22.3 min，很快能達(dá)到飽和值，所以232Th的中子俘獲截面(n,γ)大小是影響233Pa最大質(zhì)量的主要因素。由表2可知，能譜越軟，232Th的(n,γ)反應(yīng)截面越大。當(dāng)233Pa質(zhì)量達(dá)到最大后，使233Pa生成速率減小，233Pa飽和值也慢慢減小。

在如圖7(b)所示的注量率下，熱中子數(shù)據(jù)庫(kù)和重水堆數(shù)據(jù)庫(kù)計(jì)算結(jié)果表明，233U生成量在300 d左右達(dá)到最大值，其中重水堆和熱中子庫(kù)的最大值分別約為1.5%和1.2%。而壓水堆和沸水堆的233U質(zhì)量在1 200 d左右達(dá)到最大值，取值在3.5%?4%。這里，采用熱中子截面庫(kù)計(jì)算的233U飽和值最小，因233U與熱中子發(fā)生(n,f)反應(yīng)截面(表2)最大，比CANDUSEU.LIB大20 b左右，而生成233U的232Th(n,γ)反應(yīng)，THERMAL.LIB庫(kù)截面只比CANDUSEU大2 b左右，致使CANDUSEU.LIB比THERMAL.LIB數(shù)據(jù)庫(kù)下233U的飽和值稍大。

圖7 不同能譜時(shí)233PaInv/232Th0 (a)和233UInv/232Th0 (b)的變化Fig.7 233PaInv/232Th0 (a) and 233UInv/232Th0 (b) vs. time at different neutron spectra.

表2 不同截面數(shù)據(jù)庫(kù)相關(guān)核素與中子反應(yīng)的有效截面Table 2 Effective cross-section of neutron-induced reactions in different database.

由此可見(jiàn)，中子能譜越軟，233U/232Th0在越短時(shí)間內(nèi)達(dá)到最大，且最大值較小。要想提取最大質(zhì)量的233U，需要在能譜較硬的地方輻照到飽和，而要想使更多的233U為有效增值因子(keff)做貢獻(xiàn)，即消耗較多的233U，則需要在較軟的地方輻照。

2.3 循環(huán)輻照對(duì)釷鈾轉(zhuǎn)化率的影響

釷鈾轉(zhuǎn)化率主要是針對(duì)232Th轉(zhuǎn)化生成233U的總量來(lái)說(shuō)的，包括其存量和消耗量，定義如式(3)：

運(yùn)用同一中子能譜(THERMAL.LIB)，采用不同中子注量率，連續(xù)輻照下232Th-233U轉(zhuǎn)化率的計(jì)算結(jié)果見(jiàn)圖8。從圖8中可以看出，連續(xù)輻照相同時(shí)間，在相同能譜下，232Th-233U的轉(zhuǎn)化率與注量率幾乎成正比。注量在1.90×1022n·cm?2內(nèi)，其轉(zhuǎn)化率不超過(guò)12.7%，這表明一次通過(guò)式燃料輻照方案釷的利用率很低，所以需要尋求一種增加釷-鈾轉(zhuǎn)化率的輻照方案。

圖8 232Th-233U的轉(zhuǎn)化率Fig.8 Conversion rate of 232Th to 233U.

由圖1可知，233U由233Pa衰變而來(lái)，如果在輻照適當(dāng)時(shí)間后，將燃料從堆中拿出來(lái)，就減少了233Pa與中子發(fā)生的其它反應(yīng)路徑的消耗，233Pa將可以完全衰變?yōu)?33U，從而增加釷-鈾轉(zhuǎn)化率。所以本節(jié)采用燃料輻照模式為“輻照-冷卻-輻照-冷卻……”的開(kāi)循環(huán)模式，對(duì)循環(huán)過(guò)程中轉(zhuǎn)化率的影響因素進(jìn)行分析，包括注量率、冷卻時(shí)間、輻照時(shí)間、循環(huán)次數(shù)等，以求得到最佳的輻照循環(huán)方案。

2.3.1 注量率大小的影響

以釷燃料球每次輻照兩個(gè)月，冷卻233Pa的5個(gè)半衰期的循環(huán)方式循環(huán)6次為例，將不同注量率下的轉(zhuǎn)化率與單次輻照相同時(shí)間時(shí)的轉(zhuǎn)化率相比，結(jié)果值見(jiàn)表3。在第6次循環(huán)結(jié)束時(shí)，若注量率相差4倍(5×1013n·cm?2·s?1和2×1014n·cm?2·s?1)，轉(zhuǎn)化率和轉(zhuǎn)化率的增量也相差4倍。所以，注量率對(duì)轉(zhuǎn)化率及轉(zhuǎn)化率增量（循環(huán)相對(duì)于不循環(huán)時(shí)的增量）的影響與注量率的變化成正比。

表3 循環(huán)輻照轉(zhuǎn)化率及其相對(duì)于單次輻照的轉(zhuǎn)化率增量Table 3 Conversion rate of cyclic irradiation, in comparison with that of continuous irradiation (n·cm?2·s?1).

2.3.2 冷卻時(shí)間的影響

冷卻時(shí)間越長(zhǎng)，由233Pa轉(zhuǎn)化到233U的量就越大。但從實(shí)用性考慮，不能無(wú)限長(zhǎng)，要選擇一個(gè)合適的冷卻時(shí)間。冷卻4個(gè)半衰期的時(shí)候，核素的衰變率已達(dá)到96.9%，于是我們選擇4個(gè)、5個(gè)、7個(gè)半衰期做比較，以注量率為5×1013n·cm?2·s?1為例，輻照兩個(gè)月后，分別求出衰變4個(gè)、5個(gè)半衰期時(shí)，232Th-233U的轉(zhuǎn)化率與衰變7個(gè)半衰期時(shí)的差別，將該數(shù)值列于表4。冷卻時(shí)間選4或5個(gè)半衰期即可達(dá)到與最大值非常相近的值。所以，可根據(jù)實(shí)際情況選擇冷卻4個(gè)或5個(gè)半衰期。

表4 每次循環(huán)的衰變時(shí)間對(duì)轉(zhuǎn)化率的影響Table 4 Effect of decay time in each cycle on conversion rate.

2.3.3 輻照時(shí)間、循環(huán)次數(shù)的影響

圖4中，233Pa在連續(xù)輻照200 d時(shí)存量達(dá)到峰值，在前30 d的時(shí)候，濃度上升非?？?；30?60 d上升較快；60?100 d上升速度減慢；100?200 d曲線變得比較平緩，直到達(dá)到最大值。所以，本文分別選輻照一個(gè)月、兩個(gè)月、6個(gè)月和200 d來(lái)研究循環(huán)輻照中輻照時(shí)間對(duì)轉(zhuǎn)化率的影響。

圖9給出注量率為5×1013n·cm?2·s?1、總輻照時(shí)間相同時(shí)(180 d)，采用以下前三種不同循環(huán)方式下的轉(zhuǎn)化率：(1) 輻照一個(gè)月，冷卻5個(gè)半衰期，循環(huán)6次，共用時(shí)33個(gè)月；(2) 輻照兩個(gè)月，冷卻5個(gè)半衰期，循環(huán)三次，共用時(shí)26個(gè)月；(3) 輻照6個(gè)月，冷卻5個(gè)半衰期，循環(huán)1次，共用10.5個(gè)月。

由圖9可知，即使循環(huán)方式不同，在圖9的增量很小，經(jīng)計(jì)算在0.01%?0.02%。這主要是因?yàn)樵趫D示注量范圍內(nèi)，轉(zhuǎn)化率本來(lái)就很小，不超過(guò)0.6%。

表5給出注量率為5×1013n·cm?2·s?1時(shí)，5種循環(huán)方式下，輻照一段時(shí)間再冷卻后的值與只輻照不冷卻相比，各個(gè)輻照階段轉(zhuǎn)化率的相對(duì)增量。相對(duì)增量的計(jì)算式如下：

圖9 總輻照時(shí)間相同時(shí)不同循環(huán)方式的轉(zhuǎn)化率—：輻照過(guò)程，---：冷卻過(guò)程Fig.9 Conversion rate of different cyclic manner at the same total irradiation time.—: Irradiation process, ---: Decay process

表5 注量率為5×1013 n·cm?2·s?1時(shí)，循環(huán)輻照相對(duì)于單次輻照232Th-233U轉(zhuǎn)化率的相對(duì)增量Table 5 Conversion rate increment of 232Th to 233U compare different cyclic manner to single irradiation with the same neutron fluence rate 5×1013 n·cm?2·s?1.

由表5可看出，循環(huán)輻照中，每次循環(huán)的輻照時(shí)間越短，相對(duì)于相同總輻照時(shí)間的連續(xù)輻照，轉(zhuǎn)化率的增量越大，每次循環(huán)的輻照時(shí)間越長(zhǎng)，轉(zhuǎn)化率相對(duì)增量越不明顯。這是因?yàn)椋?33Pa半衰期比233U短得多，輻照時(shí)間越長(zhǎng)，前邊生成的233Pa已經(jīng)都衰變成233U，達(dá)到了與冷卻同樣的效果，所以轉(zhuǎn)化率增量就比短期輻照小。

可見(jiàn)，采用循環(huán)輻照冷卻的方法可以提高釷鈾轉(zhuǎn)化率，對(duì)于相同的總輻照時(shí)間，每次循環(huán)的輻照時(shí)間越短，相對(duì)于連續(xù)輻照，轉(zhuǎn)化率增量提高的倍數(shù)越多，作用越明顯。但由于圖9和表5條件下絕對(duì)轉(zhuǎn)化率本來(lái)就很小，所以采用循環(huán)輻照，與連續(xù)輻照相比絕對(duì)增量提高很小。因此，是否采取循環(huán)輻照，要根據(jù)實(shí)際情況來(lái)考慮。如果輻照時(shí)間比較長(zhǎng)（大于兩個(gè)月），就沒(méi)必要。如果輻照時(shí)間較短，循環(huán)輻照不僅可以提高轉(zhuǎn)化率，還可以節(jié)省中子注量、降低成本，所以是不錯(cuò)的選擇。

3 結(jié)語(yǔ)

通過(guò)用ORIGEN點(diǎn)燃耗計(jì)算程序，分析不同條件下釷鈾轉(zhuǎn)化規(guī)律，可得以下結(jié)論：

(1) 同一能譜，注量率越大，233U生成速率越快；但233U的最大值與注量率無(wú)關(guān)。注量相同時(shí)，最大值幾乎相同(圖6)，大概在4×1016n·cm?2。

(2) 與233U不同，能譜相同時(shí)，233Th在120 min左右達(dá)到最大值，其最大值主要決定于中子注量率的大小。而達(dá)到最大值的時(shí)間與注量率大小無(wú)關(guān)。233Pa與233Th的變化趨勢(shì)相同，即最大值的大小決定于注量率，達(dá)到最大值的時(shí)間與注量率大小無(wú)關(guān)。

(3) 在不同能譜下，能譜越硬，233U飽和值越大(圖7)，如果要提取最大質(zhì)量的233U，應(yīng)該在能譜較硬的地方輻照；能譜越軟，消耗掉的233U越多，在反應(yīng)堆中利用233U做裂變?nèi)剂?，?yīng)該選用較軟的能譜。

(4) 在不同能譜下，233Pa生成量與233U生成規(guī)律不同。233Pa在全熱中子譜下，所達(dá)到的飽和值最大；能譜越硬，飽和值越小。

(5) 采取循環(huán)輻照冷卻可以提高233Th-233U的轉(zhuǎn)化率，每次循環(huán)周期內(nèi)輻照時(shí)間越短，相對(duì)增量提高的越明顯，當(dāng)輻照時(shí)間超過(guò)兩個(gè)月時(shí)，循環(huán)輻照對(duì)轉(zhuǎn)化率增量的作用較小，與單次輻照不冷卻相比，轉(zhuǎn)化率相對(duì)增量不超過(guò)1倍。

1 張家驊, 包伯榮. 我國(guó)釷資源調(diào)查及釷-鈾燃料循環(huán)研究[A]. 張家驊. 釷-鈾核燃料循環(huán)研究[C]. 上海: 中國(guó)科學(xué)院上海應(yīng)用物理研究所, 2009: 22?26

ZHANG Jiahua, BAO Borong. Investigation of thorium resources in China and study of thorium-uranium fuel cycle[A]. ZHANG Jiahua. Study of thorium-uranium nuclear fuel cycle[C]. Shanghai: Shanghai Institute of Applied Physics, Chinese Academy of Sciences, 2009: 22?26

2 Troyanov M F. Thorium fuel utilization: options and trends[R]. Vienna: IAEA-TECDOC-1319, 2002

3 Furukawa K, Lecocq A, Kato Y, et al. Thorium molten-salt nuclear energy synergetics[J]. Journal of Nuclear Science and Technology, 1990, 27(12): 1157?1178

4 Furukawa K, Mitachi K, Kato Y, et al. Small molten-salt reactors with a rational thorium fuel-cycle[J]. Nuclear Engineering and Design, 1992, 136(1): 157?165

5 Murogov V M, Zinin A J, Ilyunin V G, et al. Fast reactors with different fuels in uranium-plutonium and combined fuel cycle[G]. IPPE-1920, 1998, 23(1): 203?206

6 Rutten H J, Teuchert E. Advanced safety features of pebble bed HTR's with Th-utilization[R]. Julich: KFA, 1993

7 Yu J Y, Wang K, Sollychin R, et al. Thorium fuel cycle of a thorium-based advanced nuclear energy system[J]. Progress in Nuclear Energy, 2004, 45(1): 71?84

8 Chang J J, Chang J P, Won I K. Dynamic analysis of a thorium fuel cycle in CANDU reactors[J]. Annals of Nuclear Energy, 2008, 35(10): 1842–1848

9 Yamamoto T, Suwarno H, Kayano H, et al. Development of new reactor fuel materials[J]. Journal of Nuclear Science and Technology, 1995, 32(3): 260?262

10 Croff A G. A user's manual for the ORIGEN2 computer code[M]. USA: Oak Ridge National Laboratory, 1980

11 熊文綱, 李文新, 王敏. 基于熱堆的釷鈾轉(zhuǎn)換過(guò)程中232U生成的模擬計(jì)算[J]. 核技術(shù), 2012, 35(5): 395?400

XIONG Wengang, LI Wenxin, WANG Min. Simulative calculation of232U productions in thorium-uranium transform process based on thermal reactor[J]. Nuclear Techniques, 2012, 35(5): 395?400

12 尹文靜, 張海青, 曹長(zhǎng)青, 等. TRISO釷鈾包覆燃料顆粒裂變氣體生成規(guī)律研究[J]. 核技術(shù), 2014, 37(1): 010601

YIN Wenjing, ZHANG Haiqing, CAO Changqing, et al. Investigation of fission product gas formation in TRISO Th/U coated fuel particles[J]. Nuclear Techniques, 2014, 37(1): 010601

CLC TL249

Simulation study on232Th-233U conversion in thermal reactors

ZHANG Haiqing LIN Jun CAO Changqing ZHU Tianbao ZHU Zhiyong

(Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Jiading Campus, Shanghai 201800, China)

Background: Thorium-uranium fuel cycle is attracting more and more attention because of its unique physical and chemical characteristics. Purpose: Realizing the utilization of thorium fuel in thermal reactors can save the valuable natural uranium resources and produce more fissile fuel, which is conducive to the nuclear energy sustainable development. Methods: ORIGEN code was used to simulate the build-up characteristics of232Th,233Th,233Pa and233U in neutron fluency rates and neutron spectra of typical Boiling Water Reactor (BWR), Pressurized Water Reactor (PWR) and Canada Deuterium Oxide Uranium (CANDU) reactor, and the effects of multiple“irradiation-cooling” cycles on the thorium uranium conversion rate under various irradiation time were analyzed by numerical comparison. Results: Simulation results showed that the maximum inventory of233Th and233Pa is irrelevant to the neutron fluency rates when neutron spectrum is fixed. When the neutron fluency rates are preset, the harder the neutron spectrum is, the larger the maximum inventory of233U will be obtained. Multiple“irradiation-cooling” cycles can improve the conversion rate of232Th-233U, but the relative increment is becoming smaller and smaller, compared to the continuous irradiation conversion rate. Conclusion: This work provides relevant theoretical basis for thermal reactors thorium-uranium fuel cycle research.

Neutron fluence rate, Neutron spectrum, Conversion rate of232Th-233U

TL249

10.11889/j.0253-3219.2015.hjs.38.050601

中國(guó)科學(xué)院戰(zhàn)略性先導(dǎo)科技專(zhuān)項(xiàng)(No.XDA02030200)資助

張海青，女，1981年出生，2011年于中國(guó)原子能科學(xué)研究院獲博士學(xué)位，從事釷鈾燃料技術(shù)研究，E-mail: zhanghaiqing@sinap.ac.cn

2014-12-09，

2015-03-05