謝全新，耿冰霜，吳建軍

（核工業(yè)理化工程研究院，天津 300180）

二元同位素分離級(jí)聯(lián)的設(shè)計(jì)通常以理想級(jí)聯(lián)［1］為基礎(chǔ)，使實(shí)際級(jí)聯(lián)的結(jié)構(gòu)盡量逼近理想級(jí)聯(lián)。理想級(jí)聯(lián)是在級(jí)的混合管處無豐度混合損耗的級(jí)聯(lián)，即效率最高的級(jí)聯(lián)。在多元即多組分同位素分離級(jí)聯(lián)中，不可能使所有組分在級(jí)的混合管處無豐度混合損耗，但可建立一種級(jí)聯(lián)，使其中某兩個(gè)組分的相對(duì)豐度在級(jí)的混合管處相等，這樣的級(jí)聯(lián)即為豐度比匹配級(jí)聯(lián)（Matched Abundance Ratio Cascade），簡(jiǎn)稱為M級(jí)聯(lián)［2－3］。多元分離的M級(jí)聯(lián)類似于二元分離的理想級(jí)聯(lián)。

某些特殊的分離任務(wù)需借助附加供料級(jí)聯(lián)，如核燃料循環(huán)中回收鈾的再濃縮［4－6］利用了附加供料級(jí)聯(lián)，以天然鈾為主供料，回收鈾為附加供料。本文將以帶附加供料的M級(jí)聯(lián)的數(shù)學(xué)模型為基礎(chǔ)，利用回收鈾濃縮的特殊性，提出一種濃縮回收鈾M級(jí)聯(lián)的近似算法。

1 帶附加供料的M級(jí)聯(lián)的數(shù)學(xué)模型

帶附加供料的M級(jí)聯(lián)（圖1）的外參量有主供料流F、附加供料流E、精料流P及貧料流W，相應(yīng)的組分豐度分別為CiF、CiE、CiP、CiW（i＝1，…，m，m為待分離同位素混合物的組分?jǐn)?shù)），級(jí)聯(lián)總級(jí)數(shù)為N。精料從第N級(jí)取料，貧料從第1級(jí)取料。此外，主供料從第f級(jí)供入，而附加供料從第l級(jí)供入。上述參量之間存在以下關(guān)系［7－8］：

級(jí)聯(lián)總相對(duì)流量可通過下列關(guān)系式來計(jì)算［7－8］：

式（1）～（5）中的gi是第i組分的精料流與貧料流之比［2－8］，可表示為：

式中：q0為基本全分離系數(shù)；Mi為第i組分的摩爾質(zhì)量；M＊為兩關(guān)鍵組分（如第n和第k組分）摩爾質(zhì)量的平均值，即：

多組分同位素分離級(jí)聯(lián)的計(jì)算從功能上可分為兩類：一類是驗(yàn)證計(jì)算，即對(duì)級(jí)聯(lián)結(jié)構(gòu)和運(yùn)行參量已確定的級(jí)聯(lián)進(jìn)行計(jì)算，來驗(yàn)證理論計(jì)算結(jié)果與實(shí)際運(yùn)行數(shù)據(jù)是否相符。對(duì)于帶附加供料的M級(jí)聯(lián)，如級(jí)聯(lián)總級(jí)數(shù)N、主供料級(jí)f、附加供料級(jí)l、主供料豐度CiF、附加供料豐度CiE及gi和E／F已知，由式（1）～（4）可求出相對(duì)流量F／P、F／W、精料組分豐度CiP和貧料組分豐度CiW，之后根據(jù)式（5）可求出M級(jí)聯(lián)的總相對(duì)流量LT／P。另一類是設(shè)計(jì)計(jì)算，即針對(duì)一定分離任務(wù)，確定級(jí)聯(lián)結(jié)構(gòu)和運(yùn)行參量，使級(jí)聯(lián)最經(jīng)濟(jì)。如給定目標(biāo)組分的精料豐度CnP和貧料豐度CnW，確定級(jí)聯(lián)結(jié)構(gòu)，即確定總級(jí)數(shù)N、主供料級(jí)f、附加供料級(jí)l，使級(jí)聯(lián)總流量最小。這類計(jì)算實(shí)際是一以目標(biāo)組分的精料豐度CnP和貧料豐度CnW為約束條件，以級(jí)聯(lián)總相對(duì)流量LT／P為目標(biāo)函數(shù)，以總級(jí)數(shù)N、主供料級(jí)f、附加供料級(jí)l為優(yōu)化變量的帶約束條件的非線性規(guī)劃問題。該類計(jì)算較第1類計(jì)算復(fù)雜得多，必須借助專門的工程優(yōu)化方法來完成。

圖1 帶附加供料的M級(jí)聯(lián)示意圖Fig.1 Scheme of M cascade with additional feed flow

2 用于回收鈾濃縮的M級(jí)聯(lián)的近似算法

回收鈾的濃縮必須同時(shí)考慮232U的凈化和236U 的 補(bǔ) 償 問 題 。 理 論 研 究［4－5］表 明 ，濃 縮回收鈾的最有效的級(jí)聯(lián)形式是帶附加供料的級(jí)聯(lián)。天然鈾作主供料，回收鈾作附加供料。天然鈾和回收鈾的組成及其組分豐度列于表1。表1中回收鈾的組分豐度為俄羅斯輕水堆VVER－440堆后料經(jīng)10年儲(chǔ)存期后的典型豐度。

表1 天然鈾和回收鈾的組成及其組分豐度Table 1 Composition and component abundances in natural and recycled uranium

本文所討論的近似算法主要是針對(duì)第2類級(jí)聯(lián)計(jì)算，即已知目標(biāo)組分的精料豐度和貧料豐度，確定級(jí)聯(lián)結(jié)構(gòu)。毫無疑問，回收鈾濃縮時(shí)的目標(biāo)組分為235U，假設(shè)目標(biāo)組分的精料豐度C3P＝3.5%，貧料豐度C3W＝0.2%。欲確定級(jí)聯(lián)結(jié)構(gòu)，即確定級(jí)聯(lián)總級(jí)數(shù)N、主供料級(jí)f、附加供料級(jí)l，必須利用有效的工程優(yōu)化方法。但回收鈾的濃縮有其自身的特殊性。首先，作為主供料的天然鈾，只含有兩種組分235U和238U，即表1中的C1F＝C2F＝C4F＝0。其次，作為附 加 供 料 的 回 收 鈾，232U、234U 和236U 相 對(duì)235U和238U屬微量組分。這兩個(gè)特殊性允許在回收鈾級(jí)聯(lián)計(jì)算中做某些合理的近似處理，從而大為簡(jiǎn)化計(jì)算。下面將詳細(xì)討論用于回收鈾濃縮的M級(jí)聯(lián)的近似算法。

在帶附加供料的M級(jí)聯(lián)中，第i組分在4股外部流中的組分守恒方程可表示為：

同時(shí)滿足下面的H平衡方程：

其中：Rnk，F(xiàn)、Rnk，E、Rnk，P、Rnk，W分別表示級(jí) 聯(lián) 主供料、附加供料、精料與貧料流中n、k兩關(guān)鍵組分的相對(duì)豐度，即：

式（9）中的λi是與組分有關(guān)的量，可表示為：

H平衡方程由Garza［2］首先提出，起初主要針對(duì)分離系數(shù)小的不帶附加供料的普通擴(kuò)散級(jí)聯(lián)。之所以稱之為H平衡方程，是因?yàn)榉匠讨卸x了一齊次函數(shù)（H函數(shù)）。式（9）中第i組分的 H函數(shù)可表示為Hi＝CiRnk－λi，這樣，式（9）可改寫為：

由于235U 和238U 的豐度在天然鈾和回收鈾中均占有絕對(duì)優(yōu)勢(shì)，所以把他們作為關(guān)鍵組分來進(jìn)行計(jì)算是合理的，即在計(jì)算gi時(shí)，式（7）中的n＝3、k＝5。這樣，根據(jù)式（6）～（7）可計(jì)算出各組分的gi，即：

借助式（13）可把式（9）表達(dá)成如下形式：

其中：rp＝P／F；rw＝W／F；re＝E／F。

組分質(zhì)量守恒方程保留原有形式并可分列如下：

利用式（19）和下列級(jí)聯(lián)流量守恒方程：

可求得rw和rp，即：

式（23）～（24）中的re為自由變量，一般人為給定。

在回收鈾中，235U與238U兩組分的豐度在同位素混合物中占絕對(duì)優(yōu)勢(shì)，因此，可作如下近似：

這樣，通過求解方程組（14）和（17）并考慮到作為主供料的天然鈾中不含232U，即C1F＝0，可得到232U 的精料豐度C1P和貧料豐度C1W：

同理，可求得234U、236U 的精料豐度和貧料豐度：

目標(biāo)組分235U在級(jí)聯(lián)主供料、附加供料、精料和貧料流中的豐度C3F、C3E、C3P、C3W是已知的，附加供料與主供料流量之比re一般預(yù)先給定，因此，根據(jù)式（23）～（24）可求出rw、rp。由式（10）并考慮近似關(guān)系式（式（25））可求得兩關(guān)鍵組分在4股外部流中的相對(duì)豐度Rnk，F(xiàn)、Rnk，E、Rnk，P、Rnk，W。然后根據(jù)式（26）～（31）可求得232U、234U 和236U 3個(gè)組分的精料豐度和貧料豐度。238U的精料豐度和貧料豐度可由豐度歸一化條件＝1直接求出。至此，回收鈾中全部5個(gè)組分的精料豐度和貧料豐度均已確定。

在級(jí)聯(lián)設(shè)計(jì)計(jì)算中，還須確定級(jí)聯(lián)總級(jí)數(shù)N、主供料級(jí)f、附加供料級(jí)l，這3個(gè)參量可由下式直接求出：

而級(jí)聯(lián)中總的相對(duì)流量可由下式確定：

式（32）～（35）在文獻(xiàn)［7－8］中有詳細(xì)推導(dǎo)。

無論是由式（26）～（31）計(jì)算組分豐度，還是由式（32）～（35）確定級(jí)聯(lián)結(jié)構(gòu)和級(jí)聯(lián)總流量，均為簡(jiǎn)單的解析求解，無須復(fù)雜的編程計(jì)算，使用普通的計(jì)算器即可完成計(jì)算。這是近似算法的最大優(yōu)點(diǎn)。

3 近似算法與優(yōu)化算法的比較

為驗(yàn)證近似算法的合理性，將近似算法與優(yōu)化算法的結(jié)果進(jìn)行比較。根據(jù)式（1）～（5）進(jìn)行優(yōu)化計(jì)算時(shí)，先利用罰函數(shù)法將約束優(yōu)化問題轉(zhuǎn)化為無約束優(yōu)化問題，然后用單純形法求解。計(jì)算中天然鈾和回收鈾取表1所列的供料豐度，回收鈾和天然鈾的供料量之比re＝1。兩種方法得到的計(jì)算結(jié)果列于表2。

表2 優(yōu)化算法與近似算法的計(jì)算結(jié)果Table 2 Calculation results obtained by using optimization method and approximate method

由表2可看出，近似算法與優(yōu)化算法的計(jì)算結(jié)果非常接近。通過改變r(jià)e和回收鈾的組分豐度，進(jìn)行了一系列計(jì)算，也得到相同的結(jié)論，即兩種算法的結(jié)果非常接近，最大相對(duì)偏差小于0.1%。因此，對(duì)于濃縮回收鈾的 M級(jí)聯(lián)，采用本文所提供的近似算法是合理的，且近似算法更方便和快捷。

4 結(jié)論

利用濃縮回收鈾的特殊性，提供了一濃縮回收鈾的帶附加供料的M級(jí)聯(lián)的近似算法。得到了求解精料組分豐度、貧料組分豐度、級(jí)聯(lián)總級(jí)數(shù)、主供料級(jí)和附加供料級(jí)以及級(jí)聯(lián)總流量的近似計(jì)算式。所有計(jì)算都無須復(fù)雜的編程，使用最普通的計(jì)算器即可算出結(jié)果。數(shù)值計(jì)算表明，近似算法與工程優(yōu)化算法得到的結(jié)果非常接近，最大相對(duì)偏差不超過0.1%，但近似算法相對(duì)于工程優(yōu)化算法具有簡(jiǎn)單、方便、快捷的優(yōu)點(diǎn)。

［1］COHEN K.The theory of isotope separation as applied to the large scale production of235U［M］.USA：Mc，Graw Hill Book Co.Inc.，1951.

［2］GARZA A D L.A generalization of the matched abundance－ratio cascade for multicomponent isotope separation ［J］.Chemical Engineering Science，1963，18：73－82.

［3］von HALLE E，TENNESSEE O R.Multicomponent isotope separation in matched abundance ratio cascades with losses from each stage［C］∥Proeeding of 1st Workshop on Separation Phenomena in Liquids and Gases.Darmstadt，Germany：［s.n.］，1987：325－356.

［4］SULABERIDZE G A，BORISEVICH V D，XIE Q X.On some separation problems on regenerated uranium application in the fuel cycle ［C］∥Proceedings ofⅨ All－Russia（International）Scientific Conference： Physical and Chemical Processes on Selection of Atoms and Molecules.Zvenigorod，Russia：［s.n.］，2004：70－77.

［5］SULABERIDZE G A，BORISEVICH V D，XIE Q X.Study of regularities of stationary mass transfer in separating cascade for enrichment of regenerated uranium ［J］.Engineering Physics，2005（3）：15－19.

［6］謝全新，李大勇，吳建軍.回收鈾分離級(jí)聯(lián)計(jì)算［J］.核技術(shù)，2009，32（1）：61－65.XIE Quanxin，LI Dayong，WU Jianjun.Cascade calculations for separation of recycled uranium［J］.Nuclear Techniques，2009，32（1）：61－65（in Chinese）.

［7］XIE Quanxin.Stationary molecular selective transfer in ordinary and multiflow cascades with large enrichment factors and losses at stages for the separation of multicomponent isotopic mixture［D］.Moscow：Moscow Engineering Physics Institute，2007.

［8］SULABERIDZE G A，BORISEVICH V D，XIE Q X.Quasi－ideal cascades with an additional flow for separation of multicomponent isotope mixtures［J］.Theoretical Foundations of Chemical Engineering，2006，40（1）：7－16.