康明亮，陳繁榮，楊永強(qiáng)，吳世軍，竇順梅

1.中國(guó)科學(xué)院 廣州地球化學(xué)研究所，廣東 廣州 510640; 2.中國(guó)科學(xué)院 研究生院，北京 100049;3.南京大學(xué) 地球科學(xué)系，江蘇 南京 210093

PHREEQC是由美國(guó)地調(diào)局(USGS)開發(fā)的一款功能強(qiáng)大的水文地球化學(xué)模擬軟件，擁有龐大的數(shù)據(jù)庫(kù)，可進(jìn)行水的混合、溫度變化的影響、液相中的溶解和沉淀反應(yīng)平衡、離子交換平衡、表面絡(luò)合平衡、定壓或定容的氣相平衡、固溶體平衡、動(dòng)力學(xué)控制反應(yīng)以及基于彌散和擴(kuò)散的一維反應(yīng)遷移等的模擬。PHREEQC具有用戶友好性、數(shù)值穩(wěn)定性、兼容性、數(shù)據(jù)格式的直觀性和靈活性等特點(diǎn)，因此在國(guó)際上廣泛使用。

1 反應(yīng)路徑模擬原理

反應(yīng)路徑表示的是一組經(jīng)歷不可逆反應(yīng)的系統(tǒng)，逐漸達(dá)到局域熱力學(xué)平衡狀態(tài)的過程。反應(yīng)路徑模擬是在化學(xué)種態(tài)計(jì)算的基礎(chǔ)上發(fā)展起來(lái)的，主要依據(jù)體系的平衡，目前它已成為研究水-巖作用的強(qiáng)有力工具。它主要基于如下幾種方程：質(zhì)量平衡方程、質(zhì)量作用定律和電荷平衡方程。

PHREEQC對(duì)反應(yīng)路徑的模擬有3種方式：(1) 通過計(jì)算反應(yīng)過程與相圖邊界的單個(gè)交點(diǎn)來(lái)獲取反應(yīng)路徑，該法并不需要知道反應(yīng)的量，但必須通過一系列的模擬去找出相平衡邊界點(diǎn)，找出這些單獨(dú)的點(diǎn)后，在相圖中就能有效地描畫出反應(yīng)路徑；(2) 通過迭代計(jì)算的方式得到反應(yīng)路徑，該法只需要做一次模擬計(jì)算，但事前必須知道合適的反應(yīng)量，這可以通過迭代方式得到，在PHREEQC中通過Reaction命令，逐步增加反應(yīng)量，可以找出反應(yīng)路徑中相平衡邊界點(diǎn)，由此能快速得到反應(yīng)路徑；(3) 通過反應(yīng)動(dòng)力學(xué)過程來(lái)獲得反應(yīng)路徑，該法必須知道反應(yīng)的動(dòng)力學(xué)表達(dá)式，通過調(diào)節(jié)計(jì)算輸出步驟找出相變轉(zhuǎn)折點(diǎn)。

本工作使用PHREEQC(Version 2.15.02)的LLNL.DATA數(shù)據(jù)庫(kù)，采用第二種反應(yīng)路徑模擬法，以甘肅北山五一井地下水化學(xué)成分為實(shí)例，模擬了以黃鐵礦為填料的屏障對(duì)氧化還原敏感的U、Se、Tc的還原阻滯作用。為更好地反映處置庫(kù)的還原環(huán)境，反應(yīng)體系選定為封閉體系，使用的核素濃度為1.0×10-4mol/L，液相中所加入的黃鐵礦量采用mmol/L計(jì)算。地下水化學(xué)組分[1]列于表1，其中pH=7.24,Eh=0.345 V，處于中性環(huán)境，由于地下處置庫(kù)中的Fe多以Fe2+形式存在，因此在無(wú)氧的地下環(huán)境中，黃鐵礦將能穩(wěn)定存在于其中。

表1 北山五一井地下水的化學(xué)組分Table 1 Ground water chemical compositions of Wuyi well in Beishan area

2 結(jié)果和討論

2.1 黃鐵礦與U(Ⅵ)的作用

鈾是乏燃料的主要組成部分，是重要的放射性核素之一。在高放廢物地質(zhì)處置中，鈾的遷移行為受氧化形態(tài)影響，U(Ⅵ)具有較強(qiáng)的遷移性，U(Ⅳ)由于溶解度較低而易沉淀固定[10]。張東等[11]研究了鈾在特定場(chǎng)址下的存在與遷移形態(tài)，得出易遷移的碳酸鈾酰占絕對(duì)優(yōu)勢(shì)。熱力學(xué)上U(Ⅵ)能被黃鐵礦還原，生成穩(wěn)定的瀝青鈾礦UO2，其溶解度很低，易沉淀[2]。初步模擬得出的瀝青鈾礦飽和指數(shù)(a)隨黃鐵礦濃度(cP)的變化列入表2。由表2可知，當(dāng)加入的黃鐵礦濃度為0.26～0.28 mmol/L時(shí)，瀝青鈾礦將達(dá)到飽和。

表2 瀝青鈾礦飽和指數(shù)隨黃鐵礦濃度(cP)的變化Table 2 Saturation index evolution of uraninite as a function of pyrite concentration(cP)

圖1 溶液中總U濃度及pH的變化圖Fig.1 Total U concentration and pH as a function of reaction progress■——總(Total)U，▲——pH

進(jìn)一步細(xì)化分步間隔,并給定附加條件即過飽和的UO2沉淀析出，對(duì)黃鐵礦濃度從0.000～0.284 mmol/L進(jìn)行模擬，結(jié)果示于圖1。從圖1可以看出，反應(yīng)過程中，溶液pH值基本呈上升趨勢(shì)，由初始溶液的7.24上升至8.05，推測(cè)U(Ⅵ)與黃鐵礦的反應(yīng)消耗了H+。當(dāng)cP>0.260 mmol/L時(shí)，瀝青鈾礦開始沉淀，當(dāng)cP=0.280 mmol/L時(shí)，溶液中的總U濃度已由1.0×10-4mol/L降至4.24×10-10mol/L，并且在此過程中所加入的黃鐵礦未達(dá)飽和。

圖2 主要溶解態(tài)鈾的濃度隨黃鐵礦濃度的變化Fig.2 Concentration of uranium species as a function of pyrite concentration■——總U(Ⅵ)，★▲▼

圖3 溶液中和Fe2+以及Fe3+的變化圖Fig.3 Concentration of S, Fe2+ and Fe3+ as a function of reaction progress■★——Fe2+，▲——Fe3+

2.2 黃鐵礦與硒的反應(yīng)

79Se是一種長(zhǎng)壽命裂變產(chǎn)物核素，能表現(xiàn)五個(gè)價(jià)態(tài)：-Ⅱ、-Ⅰ、0、Ⅳ、Ⅵ，其遷移能力隨價(jià)態(tài)升高而增大[16]。從熱力學(xué)上來(lái)看，Se(-Ⅱ)是Se的最穩(wěn)定態(tài)，F(xiàn)e(Ⅱ)能將Se(Ⅳ)和Se(Ⅵ)還原為Se(-Ⅱ)[17]。Bruggeman[4-5]的研究結(jié)果表明，黃鐵礦能還原亞硒酸鹽,降低其在液相中的濃度。Breynaert[18]進(jìn)一步證實(shí)了亞硒酸鹽與黃鐵礦反應(yīng)的還原產(chǎn)物為Se(0)。本工作使用PHREEQC初步模擬得出，當(dāng)黃鐵礦濃度大于0.262 mmol/L時(shí)，反應(yīng)生成的產(chǎn)物Se(0)和FeSe2都將達(dá)到飽和，進(jìn)一步細(xì)化分步區(qū)間，發(fā)現(xiàn)Se(0)先于FeSe2達(dá)到飽和，這從理論上印證了Breynaert的實(shí)驗(yàn)結(jié)果，即亞硒酸與黃鐵礦反應(yīng)的還原產(chǎn)物為Se(0)。若假定只有Se(0)或FeSe2沉淀，結(jié)果顯示溶液中的總Se濃度都能降至低水平，若只有Se(0)沉淀，則在黃鐵礦濃度達(dá)0.295 mmol/L后，F(xiàn)eSe2也達(dá)飽和。因此給定附加條件，即兩種產(chǎn)物只要達(dá)到飽和就可以沉淀，在黃鐵礦濃度為0.0～0.3 mmol/L區(qū)間，對(duì)反應(yīng)路徑進(jìn)行了模擬，整個(gè)模擬過程，所加入的黃鐵礦未達(dá)飽和。

溶液中總Se濃度及pH的變化示于圖4。由圖4知，黃鐵礦濃度為0.260～0.295 mmol/L時(shí)，溶液中的總硒濃度迅速下降，當(dāng)加入的黃鐵礦濃度為0.295 mmol/L時(shí)，溶液中的總硒已由初始溶液的1.0×10-4mol/L降至9.36×10-14mol/L，pH值在硒濃度降低區(qū)間有上升趨勢(shì)，由最初的7.24升至8.22，推測(cè)黃鐵礦還原硒的反應(yīng)過程消耗H+，Bruggeman[4-5]實(shí)驗(yàn)結(jié)果也證實(shí)，反應(yīng)過程中溶液的pH值有上升趨勢(shì)。從溶液中硒濃度的走勢(shì)來(lái)看，黃鐵礦能有效降低液相中Se的濃度。

圖4 溶液中總Se濃度及pH的變化圖Fig.4 Total Se concentration and pH as a function of reaction progress■——總(Total)Se，★——pH

圖5 溶液中各種形態(tài)Se及的濃度變化圖Fig.5 Concentration of various Se and S as a function of reaction progress■——Se(-Ⅱ)，★——Se(Ⅳ)，●——Se(Ⅵ)，▲

2.3 黃鐵礦與锝的反應(yīng)

圖6 溶液中總Tc濃度及pH的變化圖Fig.6 Total Tc concentration and pH as a function of reaction progress■——總(Total)Tc，★——pH

圖7 溶液中Tc(Ⅶ)及濃度的變化圖Fig.7 Concentration of total Tc(Ⅶ) and S as a function of reaction progress■——Tc(Ⅶ)，★

2.4 還原次序的比較

圖8 溶液中Se、Tc和U的濃度變化圖Fig.8 Concentration of total Se, Tc and U as a function of reaction progress■——Se，●——Tc，▲——U

3 結(jié) 論

反應(yīng)路徑模擬能較好地反映化學(xué)反應(yīng)過程中各種物質(zhì)組分的變化情況。本工作使用具有龐大熱力學(xué)數(shù)據(jù)庫(kù)的PHREEQC軟件，以中國(guó)甘肅北山五一井地下水為實(shí)例，對(duì)黃鐵礦還原U、Se、Tc的反應(yīng)路徑進(jìn)行了模擬，并對(duì)核素共存狀態(tài)下的還原次序進(jìn)行了探討。對(duì)單一核素的模擬結(jié)果表明，由于黃鐵礦具有較強(qiáng)的還原性，在其加入量達(dá)到一定程度后(未飽和)，均能將氧化還原敏感的U、Se、Tc的濃度降低至極低水平，并且在此過程中，pH呈現(xiàn)出不同程度的升高，從而能更加有效地阻滯放射性核素的遷移。對(duì)核素還原次序的模擬結(jié)果表明，共存的核素之間存在相互影響，Se優(yōu)先被黃鐵礦還原，U、Tc則因Se的存在還原相對(duì)滯后。模擬結(jié)果不僅在理論上證實(shí)了黃鐵礦還原放射性核素的可行性，而且量化了其反應(yīng)深度，對(duì)今后進(jìn)一步開展這方面的實(shí)驗(yàn)工作以推進(jìn)黃鐵礦在高放廢物地質(zhì)處置庫(kù)中的應(yīng)用具有重要的現(xiàn)實(shí)意義。

[1] 王 駒,范顯華,徐國(guó)慶.中國(guó)高放廢物地質(zhì)處置研究十年進(jìn)展[M].北京：原子能出版社,2004.

[2] Eglizaud N, Miserque F, Simoni E, et al. Uranium(Ⅵ) Interaction With Pyrite (FeS2): Chemical and Spectroscopic Studies[J]. Radiochim Acta, 2006, 94(9-11): 651-656.

[3] Naveau A, Monteil-Rivera F, Guillon E, et al. Interactions of Aqueous Selenium(-Ⅱ) and (Ⅳ) With Metallic Sulfide Surfaces[J]. Environ Sci Technol, 2007, 41(15): 5 376-5 382.

[4] Bruggeman C, Maes A, Vancluysen J, et al. Selenite Reduction in Boom Clay: Effect of FeS2, Clay Minerals and Dissolved Organic Matter[J]. Environ Pollut, 2005, 137(2): 209-221.

[5] Bruggeman C, Vancluysen J, Maes A. New Selenium Solution Speciation Method by Ion Chromatography Plus Gamma Counting and Its Application to FeS2-Controlled Reducing Conditions[J]. Radiochim Acta, 2002, 90(9-11): 629-635.

[6] Maes A, Geraedts K, Bruggeman C, et al. Evidence for the Interaction of Technetium Colloids With Humic Substances by X-Ray Absorption Spectroscopy[J]. Environ Sci Technol, 2004, 38(7): 2 044-2 051.

[7] De Windt L, Schneider H, Ferry C, et al. Modeling Spent Nuclear Fuel Alteration and Radionuclide Migration in Disposal Conditions[J]. Radiochim Acta, 2006, 94(9-11): 787-794.

[8] Aubriet H, Humbert B, Perdicakis M. Interaction of U(Ⅵ) With Pyrite, Galena and Their Mixtures: A Theoretical and Multitechnique Approach[J]. Radiochim Acta, 2006, 94(9-11): 657-663.

[9] Bruggeman C, Maes A, Vancluysen J. The Identification of FeS2as a Sorption Sink for Tc(Ⅳ)[J]. Phys Chem Earth, 2007, 32(8-14): 573-580.

[10] Li W J, Zhao J J, Li C S, et al. Speciation Measurements of Uranium in Alkaline Waters Using Diffusive Gradients in Thin Films Technique[J]. Anal Chim Acta, 2006, 575(2): 274-280.

[11] 張 東,李寬良,康厚軍.鈾在特定場(chǎng)址地下水中存在和遷移形態(tài)研究及其沉積熱力學(xué)分析[J].核化學(xué)與放射化學(xué),2004,26(1):43-47.

[12] Rimstidt J D,Vaughan D J. Pyrite Oxidation: A State-of-the-Art Assessment of the Reaction Mechanism[J]. Geochim Et Cosmochim Acta, 2003, 67(5): 873-880.

[13] Liger E, Charlet L, Van Cappellen P. Surface Catalysis of Uranium(Ⅵ) Reduction by Iron(Ⅱ)[J]. Geochim Cosmochim Acta, 1999, 63(19-20): 2 939-2 955.

[14] O’Loughlin E J, Kelly S D, Cook R E, et al. Reduction of Uranium(Ⅵ) by Mixed Iron(Ⅱ)/Iron(Ⅲ) Hydroxide(Green Rust): Formation of UO2Manoparticies[J]. Environ Sci Technol, 2003, 37(4): 721-727.

[15] Jeon B H, Dempsey B A, Burgos W D, et al. Chemical Reduction of U(Ⅵ) by Fe(Ⅱ) at the Solid-Water Interface Using Natural and Synthetic Fe(Ⅲ) Oxides[J]. Environ Sci Technol, 2005, 39(15): 5 642-5 649.

[16] Chen F R, Burns P C, Ewing R C. Se-79: Geochemical and Crystallo-Chemical Retardation Mechanisms[J]. J Nucl Mater, 1999, 275(1): 81-94.

[17] Myneni S C B, Tokunaga T K, Brown G E. Abiotic Selenium Redox Transformations in the Presence of Fe(Ⅱ, Ⅲ) Oxides[J]. Science, 1997, 278(5 340): 1 106-1 109.

[18] Breynaert E, Bruggeman C, Maes A. XANES-EXAFS Analysis of Se Solid-Phase Reaction Products Formed Upon Contacting Se(Ⅳ) With FeS2and FeS[J]. Environ Sci Technol, 2008, 42(10): 3 595-3 601.

[19] 劉德軍,范顯華,章英杰,等.99Tc在Fe,Fe2O3和Fe3O4上的吸附行為研究[J].核化學(xué)和放射化學(xué),2004,26(1)：23-28.

[20] Cui D Q, Eriksen T E. Reduction of Pertechnetate by Ferrous Iron in Solution: Influence of Sorbed and Precipitated Fe(Ⅱ)[J]. Environ Sci Technol, 1996, 30(7): 2 259-2 262.

[21] 劉德軍,范顯華,章英杰,等.99Tc在模擬地質(zhì)條件下的遷移行為研究概況[J].核科學(xué)與工程,2004,24(1)：6-11.