PUREX流程钚還原反萃過程的計算機模擬研究進展

于 婷，葉國安,*，何 輝，陳延鑫，劉金平，張 嵐

1.中國原子能科學研究院 放射化學研究所，北京 102413；2.中國科學院 上海應用物理研究所 放射化學與工程技術部，上海 201800；3.中國科學院 核輻射與核能技術重點實驗室，上海 201800；4.中國科學院 釷基熔鹽堆核能系統卓越創新中心，上海 201800

PUREX流程钚還原反萃過程的計算機模擬研究進展

于 婷1，葉國安1,*，何 輝1，陳延鑫1，劉金平1，張 嵐2,3,4

1.中國原子能科學研究院 放射化學研究所，北京 102413；2.中國科學院 上海應用物理研究所 放射化學與工程技術部，上海 201800；3.中國科學院 核輻射與核能技術重點實驗室，上海 201800；4.中國科學院 釷基熔鹽堆核能系統卓越創新中心，上海 201800

PUREX流程為當前后處理工業的主流流程，其計算機模擬研究為研究熱點。國外一些國家已進行全流程模擬計算，能夠開展工藝條件分析和工藝優化工作，具有重要的應用價值。鈾钚分離工藝單元(1B)和钚反萃單元(2B)是PUREX流程的重要環節，二者計算機模擬的基礎為钚的還原反萃單元模型。本文總結了國外PUREX流程計算模擬程序中的钚還原反萃模型的研究進展，重點對模型的建立和算法做了介紹。

PUREX流程；計算機模擬；反萃；钚還原反應

核燃料后處理是實現核燃料閉式循環的重要環節，是減少需要長期地質處置的放射性廢物的唯一方法和充分利用鈾資源的有效途徑[1]。當前世界各國采用的成熟的溶劑萃取PUREX流程(plutonium uranium reduction extraction)是復雜的多級逆流萃取流程[2]，該流程中的鈾钚分離工藝單元(1B)是整個流程的分水嶺，其運行工藝情況對后續工藝單元有重要的影響，鈾钚分離是通過使用合適還原劑將Pu(Ⅳ)還原成Pu(Ⅲ)，因Pu(Ⅲ)在有機相(30%(體積分數，下同)TBP/煤油)中的分配比非常低，基本不被30%TBP/煤油萃取從而將钚反萃到水中，而鈾可以繼續保留在水相中，從而實現了鈾钚分離[3]；钚還原反萃單元(2B)則是得到合格钚產品的必需步驟，與1B單元一樣對钚進行還原反萃，從而得到了更純凈的钚產品液，以上兩個工藝單元的操作中包含了諸多組分在水相中發生的氧化還原反應，以及在兩相間的分配反應，而Pu(Ⅳ)的還原反萃過程以及過程中系統各組分在兩相間的分配是該工藝單元的關鍵，二者的數學模型——钚還原反萃工藝單元模型是鈾钚分離工藝單元和钚反萃單元進行計算機模擬研究的基礎。

1 PUREX流程

當前世界各國乏燃料后處理采用成熟的溶劑萃取PUREX工藝流程[4]，此流程的主體部分大多是二循環流程——共去污-分離循環、钚凈化循環和鈾凈化循環，而對于鈾钚來講，各自都經歷了兩個循環，其主要的流程[5]示于圖1。

圖1 PUREX流程圖Fig.1 Flow chart of PUREX

1.1 Pu(Ⅳ)的還原反萃工藝過程

整個PUREX流程中涉及钚的還原反應過程單元為鈾钚分離工藝段(1B、1BXX) 和钚線二循環工藝段(2B、2BXX)，其中鈾钚分離是通過使用合適還原劑將Pu(Ⅳ)還原成Pu(Ⅲ)，2B則是將來自2A單元含有鈾和裂變產物的有機相料液進行還原反萃，從而得到更為純凈的钚產品。

1.2 Pu(Ⅳ)還原反萃過程還原劑

PUREX流程中鈾钚分離單元和钚純化過程均為采用合適的還原劑將Pu(Ⅳ)還原成為不易被30%TBP萃取的Pu(Ⅲ)，反萃進入水相，這個過程中還原反萃劑的選擇以及適當的工藝條件的確定是實現鈾钚分離單元的關鍵，目前钚的價態調整主要是依靠化學還原法，電化學還原法尚處于實驗研究階段[3,6]。

在后處理循環中為了獲得滿意的鈾钚分離效果和钚的凈化效果，近些年來，研究者們對可用于钚還原反萃過程的還原劑進行了研究，多集中在應用體系將钚還原至三價、同時不會在反應體系殘留不易分離和不易被破壞的鹽類的化學試劑上，又稱無鹽試劑，主要的研究對象有羥胺及其衍生物、醛類及其衍生物、肼類及其衍生物、脲及衍生物、肟及其衍生物、羥肟酸以及酰肼類等，這些無鹽試劑中只含有C、N、O、H等四種元素，過程中未反應還原劑的氧化產物多為氣體和水等[7-9]。

目前在PUREX流程工廠生產中采用較多的钚還原劑為硝酸羥胺(HAN)-N2H4和U(Ⅳ)-N2H4，如在法國的馬庫爾廠、阿格UP2廠、UP3及UP2-800廠，德國的卡魯斯厄爾廠，俄羅斯的VVER后處理廠以及日本的東海后處理廠和六所村后處理廠，英國的THORP廠以及印度的塔拉普廠均采用U(Ⅳ)-N2H4作為鈾钚分離單元的還原反萃劑。法國的UP2-400、UP3和UP2-800，美國的薩瓦那河后處理廠則采用硝酸羥胺作為钚純化過程的還原劑[3,6]。同時，各國針對PUREX流程钚還原反萃過程的計算機模擬主要以HAN-N2H4和U(Ⅳ)-N2H4還原劑為模擬對象[10-16]。

1.3 萃取設備

PUREX流程鈾钚分離工藝單元和钚反萃單元普遍采用的液液萃取設備為混合澄清槽和脈沖萃取柱，前者是級式萃取設備，后者為微分式萃取設備[2]。

混合澄清槽的每一級由混合室和澄清室兩部分組成，之間采用隔板隔開，級間采用級板間隔，在混合室要外加攪拌設備使得兩相混合，結構較為復雜，兩相在澄清槽中的停留時間較長。因為設備萃取過程中能夠分出明確的萃取級，且兩相在級式混合澄清槽中混合較為充分，此類設備的數學模型在流體狀態上均可以采用全混模型開展模型函數的建立工作，之間的區別為考慮的參數的量不同[5]。

與混合澄清槽不同，脈沖萃取柱的主體為內裝篩板的圓柱筒體，在柱體外部安裝脈沖發生器，使用脈沖實現柱內液體分散均勻，柱體結構相對簡單，兩相流體在脈沖萃取柱內傳質始終處于不平衡狀態，無法劃分出明確的萃取級，組分的傳質與濃度差以及外界能量輸入有很大關系，脈沖萃取柱的數學模型與混合澄清槽也有很大的區別，國外大部分研究者采用擴散模型來描述脈沖萃取柱內的流體流動和傳質行為[2,17]。

英國溫茨凱爾二廠、日本的東海后處理廠、德國的卡魯斯厄爾廠以及俄羅斯的RT-1廠的鈾钚分離工藝單元和钚反萃單元采用混合澄清槽作為萃取設備。美國的薩瓦納河廠，英國的THORP廠，法國的馬庫爾廠、UP3和UP2-800廠，日本的六所村廠等的鈾钚分離單元和钚反萃單元使用脈沖萃取柱作為萃取設備[3]。

2 PUREX流程計算機模擬研究現狀

20世紀20年代在化工生產中開始使用圖表法預測萃取操作的分離效率和產品質量，50年代開始將計算機應用于建立化工萃取過程的數學模型，并逐漸取代了圖表法。在化工萃取中，尤其是在核燃料后處理中使用計算機來實現模擬整個流程，并使用其解出用圖表法無法得出的數值解，不僅可以計算出各液流的出口濃度、萃取器中的各級濃度，有時甚至可以實現透視萃取設備中的動態過程，因此計算機模擬成為了后處理流程設計優化和安全分析的有效工具[18-19]，并且后處理中化工操作的放射性活度高等特點也急切需要使用更安全的方法來預測結果，優化生產條件，進行安全分析，實現放大設計。

國內外對于PUREX流程的模擬已做了大量的研究工作，自20世紀70年代以來，隨著計算機技術的迅猛發展和計算機的普及，大量的后處理模擬程序被開發出來，主要以美國、法國、日本以及印度等國家為代表，這些模擬程序計算出的濃度剖面與實驗得到的數據吻合良好，收斂性也較好。

美國開發的SEPHIS(MOD1—MOD4)[20-23]以及SOLVEX[24]程序主要可以實現PUREX流程混合澄清槽的靜態和動態模擬，SEPHIS-MOD4[21]還可以模擬THOREX流程，橡樹嶺實驗室開發的PUREX流程瞬態模擬軟件(SEPHIS-MOD4)與Sandia實驗室的分離與萃取模擬模型(separations and safeguards performance model, SSPM)[25]相結合，開發形成了一個設計和評估軟件，為未來的安全分析和全流程核材料衡算提供技術幫助[26]。

日本開發的模擬程序主要有PULCO[10]、EXTRA·M[11]、PARC[12]、MIXSET[13]、Revised-MIXSET[14]、MIXSET-X[15]等，PULCO主要模擬PUREX流程脈沖萃取柱中各組分的萃取行為，而EXTRA·M[11]和MIXSET[13]系列軟件可以模擬混合澄清槽中的萃取行為，PARC[12]應用范圍有所擴大，模擬對象為混合澄清槽和脈沖萃取柱。值得說明的是Revised-MIXSET[14]程序可以進行瞬態和穩態模擬，可模擬H+、U(Ⅳ)、U(Ⅵ)、Pu(Ⅲ)、Pu(Ⅳ)、NH2OH、HNO2和N2H4八種組分的化學行為，總共包括Pu(Ⅳ)的還原反應、Pu(Ⅲ)和U(Ⅳ)的氧化反應和HNO2的分解等七種氧化還原反應，并且可以進行萃取過程的工藝條件優化計算，MIXSET-X[15]是Revised-MIXSET[14]的改進版本，可以模擬31種組分的共計45種化學反應。

法國及英國等國家也開發了模擬PUREX的計算程序，分別為PAREX[16]、VISCO[27]，印度甘地原子能研究中心開發了可以模擬233U-232Th、233U-232Th-239Pu分離行為的PUThEX[28]軟件，針對FBR乏燃料后處理流程開發了SIMPSEX[29]程序，但因這些程序以核燃料后處理為應用背景，其相關報道并不詳細。

2.1 钚還原反萃過程單元模型

钚在1B和2B中的萃取分離，主要是利用钚可變價的性質，利用一定的還原劑將其還原成不易被有機相萃取的三價從而使其從有機相反萃進入水相，該過程涉及Pu(Ⅳ)的氧化還原反應，組分Pu(Ⅳ)、Pu(Ⅲ)在兩相中的萃取分配平衡等化學反應，是一個動態變化過程。一般認為，Pu(Ⅳ)的還原反萃過程的模擬計算可以按照還原劑是否溶于有機相來分類：還原劑僅溶于水相和還原劑溶于水相和有機相兩大類，這種分類方式是按照還原反萃過程的實際化學反應的發生位置，即可能發生在水相、兩相界面或水相、有機相和兩相界面兩大類。前一種以硝酸羥胺(HAN)為例，后一種以U(Ⅳ)為例說明，具體示于圖2和圖3。

圖2 Pu(Ⅳ)的還原反萃過程(1)Fig.2 Reductive stripping process (1) of Pu(Ⅳ)

圖3 Pu(Ⅳ)的還原反萃過程(2)Fig.3 Reductive stripping process (2) of Pu(Ⅳ)

HAN在水相中將Pu(Ⅳ)還原：Pu(Ⅳ)a+還原劑a→Pu(Ⅲ)a；在兩相界面上也可能會發生還原反應，即有機相中的Pu(Ⅳ)在兩相界面上與水相中的還原劑反應：Pu(Ⅳ)o+還原劑a→Pu(Ⅲ)o→Pu(Ⅲ)a，生成在兩相中分配比較低的Pu(Ⅲ)，從而反萃進入水相。還原反應的進行使得水相中的Pu(Ⅳ)濃度逐漸降低，打破了其在兩相中的分配平衡，使得以下反應逐漸向右移動，從而實現了Pu的反萃：

U(Ⅳ)在有機相和水相均有一定的分布，還原反應在水相、有機相和兩相界面上發生，過程比還原劑僅溶于水相中的還原反萃過程要復雜， 在水相中Pu(Ⅳ)被U(Ⅳ)還原成為Pu(Ⅲ)，反應如下：

還原劑U(Ⅳ)在有機相中有一定的分布，因此可以將有機相中的Pu(Ⅳ)還原成為Pu(Ⅲ)，而Pu(Ⅲ)基本不被有機相(30%TBP/煤油)萃取，因此所生成大部分的Pu(Ⅲ)被反萃進入水相，此外在兩相界面上可能還存在水相中U(Ⅳ)將有機相中的Pu(Ⅳ)還原成為Pu(Ⅲ)，這三種Pu(Ⅳ)的還原反應使得水相中的Pu(Ⅳ)濃度降低，從而使得下式反應向右進行，有機相中的钚逐漸減少，從而實現了反萃：

但現在公開報道的還原反萃模型大多只考慮了Pu(Ⅳ)在水相中的還原反應，忽略了界面以及有機相中的反應[11-12,21,23]。

2.2 PULCO模擬軟件中的钚還原反萃模型

PULCO是由日本東海研究所開發的一套用于模擬PUREX流程的模擬程序，可模擬萃取設備為脈沖萃取柱的各個工藝單元的運行情況，該程序可以對鈾钚分離工藝單元模擬計算，使用實驗規模的脈沖萃取柱對模擬軟件進行了驗證，計算值與實測值符合良好[10]。

PULCO程序使用擴散模型來模擬脈沖萃取柱中的萃取以及反萃取操作，該模型認為脈沖萃取柱是一連續的萃取設備，柱內兩相中組分濃度的變化由組分間化學還原反應和傳質反應兩部分構成，柱內的傳質過程包括相際間傳質和因濃度梯度導致的擴散傳質，即每一相中存在由高濃度向低濃度的傳質，后者符合Fick擴散定律。對柱內某微分段內兩相分別建立物料守恒平衡方程，具體如下：

(1)

(2)

(3)

式中：φ表示存留分數；t表示時間；y為有機相濃度；x為水相濃度；E為軸向擴散系數；z表示高度；r表示化學反應速率；a表示兩相接觸表面積；D表示分配比；u表示流速；K表示傳質速率；下標w表示組分，x表示水相，y表示有機相；上標*表示平衡濃度。

PULCO[10]程序的鈾钚分離單元操作的模擬，以HAN為還原劑、肼為支持劑為例說明，水相中的Pu(Ⅳ)的濃度變化主要由傳質和钚的氧化還原反應兩部分組成，其中傳質又可以分為單一水相沿脈沖柱軸向方面存在的濃度梯度導致的軸向擴散傳質和兩相間的相際間傳質即有機相向水相傳質，氧化還原反應為水相中的還原劑HAN和肼將Pu(Ⅳ)還原成Pu(Ⅲ)而導致的Pu(Ⅳ)濃度減少，有機相中的Pu(Ⅳ)濃度變化與水相類似，即有機相沿脈沖柱軸向方向存在的濃度梯度導致的擴散傳質和相際間傳質等。而對于Pu(Ⅳ)的還原反萃的實際過程是僅溶于水相的HAN將水相中的Pu(Ⅳ)還原成Pu(Ⅲ)，從而打破了Pu(Ⅳ)在水相和有機相間的分配平衡，使得有機相中的Pu(Ⅳ)源源不斷地被反萃進入水中，繼而又被HAN還原，直至反應達到平衡。其中，Pu(Ⅳ)反萃進入水相的過程實際上為一傳質過程，PULCO對此建立的還原反萃計算模型將相際間傳質、軸向擴散傳質和氧化還原反應視為計算方程中幾個線性無關的變量，從而簡化了數學模型及算法，可以通過差分方法計算得出模型計算結果。具體過程示于圖4。此外模型還假設，在某一微分段中，钚的還原率較低，則氧化還原反應的速率是钚在兩相間的傳質速率的600倍之多，因此在該模型中假定相際間傳質速率為一定值，并忽略了在兩相界面上發生的氧化還原反應，認為還原反萃過程中的氧化還原反應只在水相內發生。其中，模型中涉及的钚在兩相間的傳質系數，是采用上升單液滴法實驗得到的，總傳質系數被表達成為钚的分配比和邊界膜傳質系數構成經驗公式，具體如下：

(4)

其中：Koc,w是組分w的總傳質系數；Diw為組分w在兩相間的分配比；kdw和kcw分別是組分w在分散相和連續相一側的邊界膜傳質系數。

圖4 Pu(Ⅳ)還原反萃過程中的化學反應和傳質Fig.4 Chemical reaction and mass transfer in plutonium stripping unit by reduction

2.3 EXTRA·M模擬軟件中的钚還原反萃模型

日本原子能研究所開發的EXTRA·M程序主要是模擬以混合澄清槽為萃取器的PUREX流程，可以模擬U(Ⅳ)、硝酸羥胺為還原劑的鈾钚分離單元和钚還原工藝單元的運行。該程序對于混合澄清槽的模擬按照實際情況分成混合室和澄清室兩部分進行計算，對每個室都進行體積流量衡算和物質質量衡算。其中體積流量衡算為進入單元操作(如混合室)的水相或有機相體積等于單元中滯留的體積和流出的流體體積之和[11]。兩相中組分的質量衡算則可以用下式表示：

(5)

(6)

每個單元操作中，組分w在水相和有機相間分配平衡，可用下式表示：

其中：x表示水相中組分濃度；y為有機相中組分濃度；V表示總體積流量；F表示每個操作單元(混合室或澄清室)的液體體積；D表示分配比；下標a表示水相，o表示有機相，j表示第j級操作，w表示組分；R表示某萃取級中的化學反應項，即表示因化學反應導致的物質的濃度的變化，可以分為水相和有機相中化學反應，具體如下式表示：

(7)

EXTRA·M模擬程序認為水相中發生的氧化還原反應包括U(Ⅳ)對Pu(Ⅳ)的還原和硝酸羥胺對Pu(Ⅳ)的還原，其中U(Ⅳ)還原钚的氧化還原反應速率公式為：

(8)

硝酸羥胺還原Pu(Ⅳ)的氧化還原反應速率公式為：

(9)

有機相中發生的主要的氧化還原反應為U(Ⅳ)對Pu(Ⅳ)的還原，其化學反應速率公式為：

(10)

式中：c表示各物質在水相或有機相中的濃度；ca,N表示水相中的總硝酸根濃度。 EXTRA·M程序中對于钚的還原反萃過程的計算分成兩部分，以混合澄清槽的一級中的混合室為例說明，首先對進入混合室的兩相分別計算其中發生的氧化還原反應，反應時間為一級的停留時間，然后計算兩相的分配平衡后的濃度，即求解分配比模型，得到了萃取平衡后的平衡濃度，EXTRA·M程序如此計算是基于混合澄清槽中的萃取反應為平衡穩態的假設上的，模型對實際還原反萃過程做了簡化處理，而實際上氧化還原反應和萃取反應是一個動態的耦合化學反應，簡化處理主要目的是方便于計算機求解計算。對于還原反萃過程的模擬計算可以用圖5表示。

圖5 EXTRA·M程序钚還原反萃模型Fig.5 Model of plutonium stripping by reduction unit in EXTRA·M

2.4 PARC模擬軟件中的钚還原反萃模型

PARC程序是日本原子能機構針對PUREX流程中鈾、钚、镎以及裂變產物行為的計算程序(program for advanced extraction process with radiation effect calculation)，模擬設備對象為脈沖萃取柱和混合澄清槽，可進行穩態和瞬態計算[15]。

PARC程序對脈沖萃取柱中钚還原反萃過程建立的模型與PULCO程序類似，即脈沖萃取柱是連續的萃取設備，柱內兩相中組分濃度的變化由組分間化學還原反應和傳質反應兩部分構成，柱內的傳質過程包括相際間傳質和因濃度梯度導致的擴散傳質，氧化還原反應則是在水相和有機相同時發生Pu(Ⅳ)還原反應，模型將同一微分段內發生的傳質反應和钚的還原反應分開考慮，分別計算，然后將這兩方面的影響加和得到脈沖柱某一微分段內組分Pu(Ⅳ)的濃度變化。

對于混合澄清槽內的钚還原反萃過程所建立的模型與EXTRA·M程序一樣，不過與后者不同的是，PARC程序中使用的钚還原反應動力學公式來自于文獻[30]，具體如下所示：

水相中U(Ⅳ)為還原劑對Pu(Ⅳ)的還原反應動力學公式為：

(11)

有機相中U(Ⅳ)為還原劑對Pu(Ⅳ)的還原反應動力學公式為：

(12)

水相中硝酸羥胺為還原劑，Pu(Ⅳ)的還原反應動力學公式為：

(13)

其中：c表示濃度；t表示時間；T表示溫度；下標a表示水相，o表示有機相。

2.5 PAREX模擬軟件中的钚還原反萃模型

法國原子能機構研究開發的PAREX軟件可以對PUREX流程進行全流程模擬，可以對硝酸羥胺、U(Ⅳ)為還原劑的钚反萃單元、鈾钚分離單元進行模擬計算，該模型是在實際分離過程涉及的所有傳質反應和氧化還原反應的實驗結果的基礎上建立的經驗模型，并且使用實驗室規模的臺架實驗數據對其進行了驗證，可用于工廠規模的鈾钚分離單元和钚反萃單元的設計和優化計算[16]。

模型假設鈾钚分離單元中發生了兩類反應，即組分在水相和有機相間的分配和特定組分在兩相匯總發生的氧化還原反應，以下三類數據為基礎：各組分在兩相間的分配比，各種組分在兩相界面上的傳質機理以及各組分間發生的氧化還原反應的相關數據，其中的分配比模型采用了SEPHIS-MOD4[23]程序中的分配比模型，使用各組分萃取反應的表觀平衡常數K和自由TBP濃度來計算分配比D，表觀平衡常數為平衡水相總鹽濃度的函數。傳質模型的傳質機理可以表示為：

其中：F為傳質總量；Vd為分散相體積；K為總傳質系數；a為傳質表面積；y為有機相中組分濃度；y*為水相匯總某組分濃度相平衡的有機相組分濃度；該傳質模型認為K·a值因組分和混合程度而不同。因為PAREX程序所建立的钚還原反萃模型(圖6)由傳質模型和氧化還原反應模型兩部分構成，其氧化還原反應模型與其他模型不同，將發生的化學反應進行了細化，分為以下四類：有效反應、干擾反應、有用的反應(可以阻礙干擾反應進行的反應)以及干擾邊界反應，各個反應的動力學公式均來自于實驗測得，具體細節并無太多文獻披露。

圖6 PAREX程序钚還原反萃模型Fig.6 Model of plutonium stripping by reduction unit in PAREX

2.6 MIXSET-X模擬軟件中的钚還原反萃模型

MIXSET系列軟件是日本東海研究所開發的一系列模擬PUREX萃取流程的模擬軟件，主要可以模擬混合澄清槽的操作，MIXSET-X是20世紀90年代末的最新版本，可以對東海(Tokai)后處理廠的整個溶劑萃取過程進行模擬計算[15]。該軟件可以對硝酸羥胺、U(Ⅳ)為還原劑的钚反萃單元、鈾钚分離單元進行模擬計算，并使用實驗規模的數據和計算結果做了比對驗證。該軟件的還原反萃單元模型獨特之處在于分配比計算中加入了級效率，對分配平衡計算加以校正，最終使用表觀分配比進行萃取計算，模型認為氧化還原反應在兩相中進行，其反應速率分別使用兩相氧化還原反應動力學模型經驗公式對其上一步得到的兩相各組分平衡濃度進行氧化還原反應計算，具體過程示于圖7，采用的計算公式與EXTRA·M軟件[11]基本相同，具體示于圖7。程序中使用的钚還原反應動力學公式具體如下：

水相中U(Ⅳ)為還原劑對Pu(Ⅳ)的還原反應動力學公式為：

(14)

圖7 MIXSET-X程序钚還原反萃模型Fig.7 Model of plutonium stripping by reduction unit in MIXSET-X

水相中硝酸羥胺為還原劑，Pu(Ⅳ)的還原反應動力學公式為：

(15)

有機相中U(Ⅳ)為還原劑對Pu(Ⅳ)的還原反應動力學公式為：

(16)

MIXSET-X軟件的相關報告對硝酸羥胺為還原劑時的钚還原反萃計算結果做了說明，對比結果發現計算結果與實驗結果偏離較大，愈靠近反萃劑進口端(1BW)，計算誤差越大，甚至超過100%，分析原因是，組分的還原反應和其在兩相間的分配是同時進行的動態反應，MIXSET-X系列軟件假設先進行分配平衡，然后認為還原反應在分配平衡的條件下進行，與實際的反應情況不符合，盡管在計算中加入了級效率模型對計算過程進行校正，但從實際計算結果來看，并未起到預期效果。

2.7 SEPHIS-MOD4模擬軟件中的钚還原反萃模型

SEPHIS(solvent extraction processes having interacting solutes)系列軟件[21-23]是由美國橡樹嶺國家實驗室(ORNL)開發的模擬PUREX后處理萃取流程的模擬軟件，第一套開發于1972年，在之后對軟件進行了升級和改進，分別在分配比模型、萃取過程化學反應等幾方面進行了改進，尤其是SEPHIS-MOD4軟件，可以對存在氧化還原反應的反萃過程進行計算。

SEPHIS-MOD4主要可以模擬以混合澄清槽為萃取設備的PUREX流程，主要采用的模型為級模型(stage model)[23]，該軟件可以模擬硝酸、U(Ⅳ)、U(Ⅵ)、Pu(Ⅳ)以及Pu(Ⅲ)等組分在兩相中的分配反應。其中對于涉及钚的氧化還原反萃的過程——1B和2B過程主要考慮了以U(Ⅳ)和HAN為還原劑的情況，還原反萃模型假設化學反應全部在混合室以及澄清室的水相中進行，認為有機相中的有效組分是水相反應的緩沖，通過歸一化計算將有機相中有效組分歸入水相中進行氧化還原反應計算，化學反應計算后再進行分配比模型計算，即進行傳質計算，實際上是對化學反應和傳質反應同時進行的過程作了簡化，模型具體示于圖8。

圖8 SEPHIS-MOD4程序钚還原反萃模型Fig.8 Model of plutonium stripping by reduction unit in SEPHIS-MOD4 program

與其他模型不同的是，SEPHIS-MOD4軟件中钚的還原反萃模型是在混合澄清槽萃取平衡的級效率模型基礎上建立的，該反萃模型認為各個組分在每一級均達到萃取平衡，但是不一定達到氧化還原反應的平衡，氧化還原程度用反應效率來表征，采用反應效率經驗模型來預測反應效率，該經驗模型主要以還原劑濃度為參數。其中該模型中采用的氧化還原反應動力學模型主要如下式所示，其模型與MIXSET-X[15]等軟件相似：

水相中U(Ⅳ)為還原劑對Pu(Ⅳ)的還原反應動力學公式為：

(17)

水相中硝酸羥胺為還原劑，Pu(Ⅳ)的還原反應動力學公式為：

(18)

2.8 其他模型

日本CRIEPI的Tsukada等[31]開發的模擬軟件可以模擬以脈沖萃取柱為萃取設備的PUREX流程，其中對于鈾钚分離單元(1B)的模擬采用擴散模型作為基本數學模型，分別對有機相和水相的濃度進行計算，擴散模型數學表達式如下所示：

(19)

(20)

(21)

式中：φ表示存留分數；t表示時間；y為有機相濃度；x為水相濃度；E為軸向擴散系數；z表示高度；R表示化學反應程度表達式；a表示兩相接觸面積；D表示分配比；u表示流速；上標p表示正向化學反應程度，d表示逆向化學反應程度，f表示脈沖柱外部加料；下標w表示組分，x表示水相，y表示有機相。其中R為化學反應表達式，在1B單元中即為钚的氧化還原反應速率，化學反應項采用的動力學表達式與PULCO和MIXSET相同，具體如下：

水相中U(Ⅳ)為還原劑對Pu(Ⅳ)的還原反應動力學公式為：

(22)

有機相中U(Ⅳ)為還原劑對Pu(Ⅳ)的還原反應動力學公式為：

(23)

采用實驗室規模的小型脈沖萃取柱實驗結果對模型進行了對比驗證，結果發現計算得到的各組分(钚、鈾、硝酸)的濃度剖面與實測值符合良好。

日本崎玉大學的Homma等[32]建立的U(Ⅳ)為還原劑的鈾钚分離過程模型將模型簡化為：假定全部氧化還原反應僅發生在水相中，認為有機相中有效組分并不參與化學反應，在模型計算時首先對水相中的氧化還原反應進行計算，然后再計算兩相間的傳質反應，該模型采用的氧化還原反應動力學模型公式如下式所示：

2k2c(Pu(Ⅳ))c(HNO3)

(24)

式中：r′表示氧化還原反應總量；k1表示钚與U(Ⅳ)的反應速率；k2表示與HNO3的反應速率。Homma所建立的數學模型從本質上分析與 EXTRA·M模型相似，后者只是多進行了有機相中氧化還原反應計算。

2.9 小結

現已公開的模擬PUREX流程的計算機模擬軟件，SEPHIS(MOD1—MOD4)系列程序、PARC軟件、MIXSET系列程序、PULCO程序以及EXTRA·M等均可對U(Ⅳ)-N2H4和HAN-N2H4為還原劑的钚還原反萃工藝單元和鈾钚分離工藝單元進行計算機模擬計算。其中SEPHIS-MOD4是在水相中氧化還原反應經驗數據的基礎上建立的還原動力學模型，并以此作為钚還原反萃模型進行相應的模擬計算。PARC軟件、MIXSET系列程序、PULCO程序以及EXTRA·M等則是參考文獻[30]相關動力學經驗模型建立的钚還原反萃數學模型，這些程序不僅可以計算U(Ⅳ)和HAN對Pu(Ⅳ)還原反應，區別于SEPHIS系列模型的是，后幾種程序建立的還原反萃模型中還包括了還原產物Pu(Ⅲ)的再氧化以及有機相中發生的反應。雖然各個國家建立的單元模型公式形式不一，但基本可以分為兩類。僅拿HAN還原Pu(Ⅳ)過程為例說明，第一種模型是認為氧化還原反應僅在水相中發生，有機相的有效組分如Pu(Ⅳ)等是水相中化學反應的緩沖劑，認為有機相中的有效組分也參與了水相中的氧化還原反應，通過物料守恒歸一化到水相中進行氧化還原反應計算，并假設兩相間的傳質反應瞬間完成，該種模型可以認為是先歸一化，再進行氧化還原反應計算，然后進行傳質計算，主要以美國的SEPHIS(MOD1—MOD4)為代表，這種模型認為有機相中的組分全部參與氧化還原反應，這與實際的傳質與水相氧化還原反應同時進行相互影響的情況偏離過大，從而導致氧化還原反應的計算程度比實際值偏大。另外一種模型雖同樣認為氧化還原反應僅在水相中發生，但與上一種模型不同的是，認為有機相中的組分不參與化學反應，在進行一級計算時，首先進行水相中的氧化還原反應計算，然后進行兩相間的傳質反應計算，該模型忽略了有機相的有效組分在氧化還原反應進行的同時因濃度降低而逐漸向水相傳質的情況，這樣導致氧化還原反應的計算程度比實際情況偏小，與前一種模型屬于兩個極端，并且這種模型在計算有機相中有效組分為Pu(Ⅳ)、水相中有效組分為HAN時，因還原劑和氧化劑存在于不同的相中，導致計算結果僅為組分的傳質反應結果，與實際情況差之甚遠，該種模型主要以Homma等[32]日本科研機構所建立的模型為代表。U(Ⅳ)為還原劑的氧化還原反應雖然有的模型考慮了有機相中發生的化學反應，但也是傳質反應和氧化還原反應先后計算。

這些钚還原反萃單元計算模型的一個共同點是：將整個反應時間認為是一個整體，無論是氧化還原反應僅在水相發生，還是有機相也參與了反應，出于建模方便以及計算方便的考慮，均假設氧化還原反應和傳質反應是兩個過程，這與實際還原反萃單元氧化化學反應與傳質反應相耦合的情況并不符合，并未考慮實際工藝中的反應萃取機理，從這個角度考慮，單元反應結果的計算值與實際值存在一定的誤差也不足為怪。法國的PAREX程序對化學反應進行了細化，對Pu(Ⅳ)的還原反萃過程建模同時考慮了傳質反應和氧化還原反應，其模擬結果與實際運行結果據報道符合良好，但其模型并無太多文獻報道[16]。

除法國PAREX模擬程序外的所有模擬程序均是最小限度考慮還原反萃過程中的化學反應的經驗模型，對實際的萃取反應和氧化還原反應的關聯影響作用并未考察，實際工藝中的化學反應機理并未在模擬軟件中進行考察，據報道[31]一些軟件的鈾钚分離工藝單元的模擬計算結果與實際實驗結果出入較大。另外一些國家的研究人員對鈾钚分離電化學還原過程還進行了計算機模擬研究[33]。

3 展 望

計算機模擬技術已成為后處理流程設計和工藝優化的有效工具和重要手段。PUREX流程中的鈾钚分離工藝單元(1B)和钚還原反萃工藝單元(2B)是流程的重要步驟，二者因過程包含氧化還原反應以及萃取傳質反應等成為流程模擬中的難點單元，二者均以钚的還原反萃過程模擬為基礎。國外很多研究人員針對流程的計算機模擬做了大量的工作，并編寫了相應的計算機模擬程序，這些程序對于工藝研究有重要的意義。但大多數模擬程序對于涉及钚還原反萃的模擬多從經驗數據出發，出于建模以及計算的方便，對模型做了簡化，將氧化還原反應和傳質反應分開考慮，而忽略實際反應中氧化還原反應和傳質反應是互相影響的耦合過程，這使得計算的準確度降低，限制了模擬軟件的應用范圍。國內針對PUREX流程的計算機模擬研究工作尚屬起步階段，目前尚無針對1B、2B的計算機模擬程序，因此有必要開展钚的還原反萃動力學模擬研究，建立還原反萃單元模型，并以此為基礎完成鈾钚分離工藝單元1B和钚還原反萃工藝單元2B的模擬計算。

[1] 章澤甫，王俊峰，張天祥，等.動力堆核燃料后處理工學[M].北京：原子能出版社，2013.

[2] 陳延鑫，何輝，唐洪彬，等.后處理工藝Purex流程計算機模擬研究現狀及展望[J].核化學與放射化學，2012，34(4)：193-200.

[3] 任鳳儀，周鎮興.國外核燃料后處理[M].北京：原子能出版社，2006.

[4] 姜圣階.動力堆核燃料后處理廠設計[M].北京：原子能出版社，1996.

[5] 姜圣階，任鳳儀，馬瑞華，等.核燃料后處理工學[M].北京：原子能出版社，1995.

[6] Marchenko V I, Dvoeglazov K N, Volk V I. Use of redox reagents for stabilization of Pu and Np valence forms in aqueous reprocessing of spent nuclear fuel: chemical and technological aspects[J]. Radiochem, 2009, 51(4): 329-344.

[7] 葉國安.PUREX流程中有機無鹽試劑的應用分析[J].原子能科學技術，2004，38(2)：152-158.

[8] 李颯，高耀斌，歐陽應根.鈾/钚分離過程中無鹽試劑的氧化還原動力學研究進展[J].核技術，2012，35(12)：929-935.

[9] Riddle C, Camphell K, Mausolf E. Overview of reductants utilized in nuclear fuel reprocessing/recycling, INL/CON-12-28 006[R]. USA: Idaho National Laboratory, 2013.

[10]Gonda K, Matsuda T. Calculation code PULCO for PUREX processes in a pulsed column, PNCT-841-82-19[R]. Japan: Power React and Nuclear Fuel Dev Corp, 1982.

[11]Tachimoli S. EXTRA·M: a computing code system for analysis of the PUREX process with mixer settler for reprocessing, JAERI-1331[R]. Japan: Japan Atomic Energy Research Institute, 1993.

[12]Tsubata Y, Asakura T, Morita Y. Development of a computer code PARC for simulation of liquid-liquid extraction process in reprocessing, JAEA-Data/Code 2008-010[R]. Japan: Tokai Works Power Reactor and Nuclear Fuel Development Corporation, 2008.

[13]Kozo G, Shoji F. Calculation code MIXSET for Purex process, PNCT 841-77-60[R]. Japan: Tokai Works Power Reactor and Nuclear Fuel Development Corporation, 1977.

[14]Kozo G, Koichiro O, Shoji F. Calculation code revised MIXSET for PUREX process, PNCT 841-79-26[R]. Japan: Tokai Works Power Reactor and Nuclear Fuel Development Corporation, 1979.

[15]Naito M, Suto T, Asakawa K, et al. A computer code for simulating the PUREX solvent extraction process, JNC T8400-99-005[R]. Japan: Tokai Works Power Reactor and Nuclear Fuel Development Corporation, 1999.

[16]Christian S, Pascal B, Binh D, et al. The simple solution modeling implemented in the PAREX code to simulate solvent extraction operations[C]. Global Nulear Fuel Reprocessing & Recycling Conference, Chiba, Japan, 2011.

[17]陳延鑫.核燃料后處理工業中脈沖萃取柱的計算機模擬研究[D].北京：中國原子能科學研究院，2009.

[18]何輝，李高亮，陳輝,等.Purex流程共去污工藝計算機穩態模擬[J].原子能科學技術，2008，42(9)：784-769.

[19]張春龍.PUREX流程中锝萃取行為的計算機模擬[D].北京：中國原子能科學研究院，2011.

[20]Richardson G L. Effect of high solvent irradiation exposures on TBP processing of spent LMFBR fuels, HEDL-TME 73-51[R]. USA: Oak Ridge National Laboratory, 1973.

[21]Mitchell A D. SEPHIS-MOD4: a user’s manual to a revised model of the PUREX solvent extraction system, ORNL-5471[R]. USA: Oak Ridge National Laboratory, 1979.

[22]Richardson G L, Swanson J L. Plutonium partitioning in the Purex process with hydrazine-stabilized hydroxylamine nitrate, hanford engineering development laboratory, HEDL-TME 75-31[R]. USA: Oak Ridge National Laboratory, 1975.

[23]Watson S B, Rainey R H. Modification of the SEPHIS computer code for calculating the Purex solvent extraction system, ORNL/TM-5123[R]. USA: Oak Ridge National Laboratory, 1975.

[24]Scotton W C. SOLVEX: a computer program for simulation of solvent extraction processes, DP-1391[R]. USA: Savannah River Laboratory, 1975.

[25]Cipiti B B. Separations and safeguards performance model(SSPM), SAND 2009-4896[R]. USA: Sandia National Laboratories, 2009.

[26]Cipiti B B, de Almeida V F, Gauld I C, et al. Coupling a transient solvent extraction module with the separations and safeguards performance model, SAND2009-6307[R]. USA: Sandia National Laboratories, 2009.

[27]Petrich G. The PUREX process computer model “VISCO” and its application[C]. Proc Int Solvent Extraction Conf, Dechema, Munich, West Germany, Vol.1, 1986: 427.

[28]Shekhar K, Koganti S B. Development of a computer code PUTHEX for simulation of233U-232Th and233U-239Pu-232Th separation flowshheets, IGC-255[R]. India: Indira Gandhi Centre for Atomic Research, 2003.

[29]Shekhar K, Koganti S B. Development of a computer code SIMPSEX for FBR fuel reprocessing flowsheets Ⅱ: additional benchmarking results, IGC-251[R]. India: Indira Gandhi Centre for Atomic Research, 2003.

[30]Suzuki A. Handbook on process and chemistry of nuclear fuel reprocessing, JAERI-Review 2001-038[R]. Japan: Japan Atomic Energy Research Institute, 2001.

[31]Tsukada T, Takahashi K, Yoshiki N. Behavior of uranium in the partitioning process of PUREX reprocessing[J]. J Nucl Sci Technol, 2008, 45(2): 179-194.

[32]Homma S, Sakamoto S, Takanashi M, et al. Development of COMPAS, computer aided process flow sheet design and analysis system of nuclear fuel reprocessing[J]. J Nucl Sci Technol, 1995, 32(4): 357-368.

[33]Wickramasinghe S R. Simulation and optimization of the electro-reduction pulsed column for the separation of uranium and plutonium[J]. Solvent Extr Ion Exch, 2002, 20(4-5): 471-490.

Present Situation of Research of Computer Simulation of Stripping of Plutonium by Chemical Reduction in PUREX Process

YU Ting1, YE Guo-an1,*, HE Hui1, CHEN Yan-xin1, LIU Jin-ping1, ZHANG Lan2,3,4

1.China Institute of Atomic Energy, P. O. Box 275(26), Beijing 102413, China;2.Department of Radiochemistry and Engineering, Shanghai Institute of Applied Physics, Chinese Academy of Science, Shanghai 201800, China;3.Key Laboratory of Nuclear Radiation and Nuclear Energy Technology, Chinese Academy of Science, Shanghai 201800, China;4.Center for Excellence TMSR Energy System, Chinese Academy of Science, Shanghai 201800, China

PUREX has been a dominant process in nuclear fuel reprocessing industry. The computer simulation of PUREX has been a research hot spot. In some other countries, computer simulation of complete flow scheme has been realized, and the simulation program can be used in analyzing the influence of input parameter on PUREX and optimizing the input parameter of PUREX. Uranium and plutonium separation unit(1B) and plutonium stripping unit(2B) are critical part of PUREX. The computer simulations of them are based on plutonium stripping chemical reduction model. This paper reviews the recent process of the model of the plutonium stripping by chemical reduction model, focuses on the model and algorithm.

PUREX process; computer simulation; stripping; reduction of plutonium

2016-04-20；

2017-03-22

于 婷(1986—)，女，山東淄博人，博士研究生，核燃料循環與材料專業, E-mail: yuting043703@126.com

*通信聯系人：葉國安(1964—)，男，安徽南陵人，博士，研究員，主要從事核燃料后處理研究，E-mail: yeguoan@ciae.ac.cn

TL249

A

0253-9950(2017)04-0257-11

10.7538/hhx.2017.39.04.0257