核廢物地質處置緩沖/回填材料研究綜述

馬立平，韓永國

(西南科技大學計算機科學與技術學院，四川綿陽621010)

緩沖/回填材料位于金屬包裝器和地質體之間，是高放廢物處置庫多重屏障系統中最后一道人工屏障，其主要功能有：工程屏障作用，維護處置庫結構的穩定性；水力學屏障作用，阻止地下水的滲流；化學屏障作用，阻滯核素遷移；導體作用，對輻射熱具有良好的熱傳導和擴散作用。因此，緩沖材料性能應具有低透水性，能阻止和延緩地下水向廢物包裝容器滲透流動；良好的膨脹性，緩沖材料經吸水后體積膨脹，可堵塞周圍介質中的所有空隙，降低介質的孔隙度；良好的吸附性，能吸附從高放廢物中泄漏出的核素，在一定厚度的緩沖材料的阻隔下，核素在進入生態圈之前，有足夠長的時間把大部分核素衰變到無害水平；良好的熱穩定性、耐輻射性，不易因輻射或因輻射導致溫度升高發生相變、變質等不利現象而影響其功能；導熱性，能傳導和散失核素衰變熱保持處置庫的穩定；良好的化學穩定性；較好的機械緩沖性，能減弱、消散巖石因熱膨脹、變形等造成的對廢物容器的應力，為廢物容器提供機械屏障；經濟性，處置庫的緩沖材料應經濟價廉、資源豐富，以減輕高放廢物處置的費用。從20世紀70年代已經開始，國內外專家經過大量的科學試驗和理論分析，一致認為以蒙脫石為主要成分的高壓實膨潤土及其集成混合物作為緩沖/回填材料是最適合的材料，并對膨潤土及其混合集成材料作為高放廢物處置緩沖回填材料的基本特性、溫度對礦物相變的影響、核素的遷移、氣體的擴散、緩沖回填材料的制備和放置技術、熱傳導性、滲透性以及多場耦合方面進行了廣泛的研究。

1 國外研究現狀分析

1.1 膨潤土的材料篩選研究

礦物純度、水的滯留特性、可塑性、低滲透性、較高的膨脹壓力和較好的導熱性等是膨潤土選材的主要特性標準，而粘粒含量(小于2 mm)、蒙脫石含量、離子交換能力、有機物質含量、鈣質含量和膨脹力等構成主要篩選指標，其中蒙脫石含量常被作為選料的重要指標。至20世紀80年代末，世界上核能利用較為先進的國家在膨潤土作為緩沖/回填材料的選材方面的工作均已基本完成。目前，世界上較有代表性的有美國的MX80(MX80是火山巖分解形成的Wyoming鈉蒙脫石粘土，由美國Collolds公司經銷。瑞典SKB的大多數緩沖材料的研究采用這種膨潤土，并被國際原子能機構(IAEA)推薦作為高放射性廢物深地質處置緩沖材料的參考標準粘土、法國的Fo-Ca(Fo-Ca是法國巴黎盆地的沉積粘土，主要成分是高嶺石一貝得石間層粘土礦物，由50%貝得石和50%高嶺石組成)粘土、西班牙的S-2(S-2是西班牙Almera地區的粘土，由火山巖熱液分解形成，主要成分是鈣蒙脫石)和日本的Kunigel VI(Kunigel VI是日本Kunimine公司生產的鈉蒙脫石膨潤土)等[1]～[3]。

1.2 膨潤土的滲透性研究

在膨潤土的滲透性研究方面。1996年，Delage和Cui等人在獲取試樣土水特征曲線的基礎上，采用瞬時斷面法計算了高壓實膨潤土的滲透系數。隨后, Delage和Cui等人又發現膨潤土在水飽和條件下，其中水的傳導符合達西定律；而在水飽和之前，水的運動則可分液體運動和蒸汽運動兩種形式；當含水量較低時，水的運動以蒸汽擴散為主，而當含水量較高時，水的運動以液體運動為主；鈉基膨潤土比鈣基膨潤土能更有效地阻止純水或離子溶液的滲透[4]。1996年Hideo Komine采用砂-膨潤土為試樣，待試樣體積膨脹穩定后，運用達西定律測定滲透系數，試驗表明回填材料砂-膨潤土的滲透系數隨著膨潤土的含量增多以及干密度的增大而降低[3]。但由于高壓實膨潤土非飽和滲透自身的復雜性，此項仍有待深入開展研究。

1.3 膨潤土的熱傳導特性研究

在膨潤土的熱傳導特征研究方面。膨潤土的熱傳導系數主要與密度、孔隙率、含水量及其所含的非粘土礦物的種類和數量有關。同一種膨潤土的熱傳導系數和比熱隨壓實密度增大而增大；在相同壓實密度時，含水量越大，其熱傳導系數越高。在1983、1984年Moss和Molecke研究了膨潤土和石英砂的混合物的熱傳導率，結果認為，在石英砂含量為30%的范圍內，混合物中石英砂的比例越大，熱傳導系數越高，這是因為膨潤土的熱傳導系數遠小于石英砂。而在混合物中，石英比例達到70%時，由于石英粒徑大，孔隙增大，熱傳導率將有所下降。Moss和Molecke根據理論推導認為膨潤土的熱傳導率遠大于空氣的熱傳導率，從而壓實膨潤土能達到處置庫所要求的熱傳導率。因此，人們普遍認為：提高其熱傳導系數的最簡單方法就是在膨潤土中加入石英。目前，許多研究機構對添加石英顆粒的大小和數量對熱一水一力學性能的影響正在進行研究。

1.4 膨潤土的多場耦合研究

在多場耦合試驗研究方面。1998年、2000年Cho等人發現壓實膨潤土隨著溫度的升高其滲透性會降低；2001年Romero等人通過試驗發現土的持水能力隨溫度上升而下降；2001年Romero 等人通過試驗發現壓實膨脹土在不同溫度條件下的體積變化特征取決于試驗材料的塑限、膨脹性以及壓實含水量；2008年Tong F.G，JingL等人[5]對緩沖/回填和近場的多場耦合進行了實驗研究；2009年Valyashko VM等人對膨潤土在高溫高壓下的熱傳導和滲透性進行實驗研究[6]。

在多場耦合理論研究方面。1981年Weststic等人研究了作為回填材料的膨潤土的滲透性與膨脹力等特征；1994年Olivella等人、2001年Rutqvist等人提出了水氣質量、能量及力學平衡方程；2000年 Cui等人、2002年Cui等人分別建立了兩場耦合的熱-應力模型和應力-吸力模型。針對非飽和土的應力應變本構模型，不少學者也展開了理論研究，其中1990年Alonso等人提出了巴塞羅那基本模型(BBM)；1992年Gens 等人提出的非線性彈性模型能較好地反映非飽和土諸多重要力學現象，而且參數較易確定；2008年Tang A-M等人、2009年TARONJ等人、 2009年Joshua Taron等人分別對溫度場、滲流場、應力場、化學場的多場耦合的建模理論進行研究，并取得一定的成果。在數值模擬及模擬程序設計方面，國外許多學者也作了大量的工作，他們開發了一些數值模擬程序如：ROCMAS、 FRACON、CODE_BRIGHT、 THAMES 以及ABAQUS-CLAY等，這些都是基于有限元而開發的程序，對于模擬熱、滲流、應力、化學等引起的耦合問題都有一定的局限性。但總的來說，能全面反映耦合作用條件下，高壓實膨潤土的行為特征和核素遷移行為的研究成果仍然有限。

2 國內研究現狀分析

國內在高放廢物地質處置庫緩沖/回填材料方面的研究源于1985年，主要由中國核工業北京地質研究院最早開展研究。自1985年以來，核工業北京地質研究院完成了系列研究工作[7]。主要包括：1996年至2000年間，進行了地表鈣膨潤土基本性能研究；從2001年起開展了深部鈉膨潤土的基本性能研究；2003年劉月妙等人[8]按照高放射性廢物深地質處置庫中緩沖材料的作用及要求，參照8個主要核國家采用的9種蒙脫石類粘土的樣品來源、礦物成分、主要物理化學特征和熱—水—力學性能進行分析比較發現中國內蒙古高廟子(GMZ)膨潤土蒙脫石含量較高，陽離子交換容量和比表面積較大，同時具有很好的不滲透性、高膨脹性和良好的導熱性能；從2003年開始，針對內蒙古興和縣高廟子膨潤土開展進一步的研究工作，并且也把它作為我國高放廢物深地質處置緩沖/回填材料的基料。

2.1 膨潤土的特性研究

在膨潤土的特性研究方面。2001年，核工業北京地質研究院的劉月妙等人[8]對未經提純處理的天然高廟子膨潤土，進行了X射線衍射分析、紅外光譜分析、化學全分析、X射線熒光稀土元素分析等多種物相分析法，并對試樣進行了差熱分析、熱重分析及掃描電子顯微鏡觀察，研究了其壓實、膨脹力和膨脹變形等特性；2003年劉月妙等人[8]又對高廟子膨潤土的物理化學特性進行了分析，結果表明其蒙脫石含量可高達80.9%，具有高陽離子交換能力和膨脹能力；同年，采用高壓固結儀進行膨潤土的膨脹力實驗，其結果表明在干密度為1.46 g/cm3和含水量為23.99%的條件下，該土的膨脹能力可達4.75 MPa。2005年葉為民等人[9]利用MIP和SEM方法研究表明，在高壓實膨脹土水化過程中，體積膨脹主要來源于膨脹土集合體結構中集合體之間的大孔的擴張。

2.2 膨潤土的滲透性研究

在膨潤土的滲透性方面。1995年，周抗寒等人[10]研制了用于膨脹土研究的多功能膨脹滲透儀，并用該儀器對壓實膨潤土進行滲透性試驗研究，發現用于高活度放射性廢物深地質處置的壓實膨潤土，其滲透系數要求低于10～11 cm/s；1998年，沈珍瑤等人[11]進行了高壓實膨潤土脫濕一吸濕過程試驗，結果表明吸濕過程中試樣孔隙比變化大于脫濕過程；從2005年開始，葉為民等人對膨潤土的緩沖性能進行了大量實驗[9]。通過試驗發現，高壓實膨潤土初始壓實密度越大，膨脹勢越大，通過理論分析發現處置庫運行狀態下的密實度是動態變化的，它與緩沖/回填材料的初始密實度、吸力以及人工屏障形成過程中可能發生的變形量等相關，且在吸濕過程中，隨著吸力的變化而變化，密實度的變化會反過來影響到非飽和緩沖/回填材料的滲透性，因此，密實度、土水特征是影響高壓實膨潤土的非飽和滲透性的關鍵影響因素。2006年陳寶等人采用滲析法和水汽平衡法來控制吸力技術，通過土水特征曲線測定試驗、環境掃描電鏡和壓汞試驗，研究高壓密高廟子膨潤土在不同吸力下的持水特性及其微觀結構特征，其研究結果表明在不同吸力作用下膨潤土的持水特性與其微觀結構之間有著密切關系；

2.3 膨潤土的熱傳導特性研究

在膨潤土的熱傳導特性性能方面。2003年劉月妙等人[8]對內蒙古興和縣高廟子(GMZ)膨潤土的熱傳導系數進行了測定，結果表明熱傳導系數與土的壓實壓力和含水量成正比；2006年溫志堅研究了GMZ-1的礦物組成、基本特性參數和GMZ-1在不同干密度、不同含水量條件下的熱傳導、水傳導、力學性能以及GMZ-1在不同干密度條件下的膨脹特性，并指出GMZ-1鈉基膨潤土具有蒙脫石含量高(75%左右)，雜質礦物相對較少的特點。

2.4 膨潤土中放射性核素遷移規律研究

在膨潤土作為緩沖回填材料阻止核素遷移的研究方面。2001年，劉月妙等人[7]進行了開放環境下GMZ膨潤土對核素90Sr、137Cs和239Pu的吸附實驗，結果表明，pH值、反應溫度和核素濃度對核素吸附有重要影響，137Cs和239Pu的吸附比隨溫度的升高有下降的趨勢，但90Sr的吸附比隨溫度的升高而增大，在低濃度時，吸附比隨濃度的增加而增加，高濃度時則相反，三種核素平衡吸附濃度約為10-5mmol/L；2004年，易發成等人[12]對四川三臺鈣基膨潤土的基本物化性能和對Pb2+,Co2+離子的吸附性能作了初步的研究，通過實驗得出膨潤土對重金屬離子Pb2+的吸附性能優于Co2+,Pb2+和Co2+的吸附比(Rd)分別為167.8 ml/g及53.3 ml/g；Pb2+和Co2+去除率(r)分別為45.87%和21.20%；從2007年至今，張玉軍等人[13]對高放廢物處置庫近場熱-滲-力耦合條件下核素遷移及多場耦合問題進行了大量研究，特別考慮了緩沖/回填層中溫度梯度水分擴散、水蒸汽擴散對水連續性及能量守恒的影響，建立了飽和-非飽和介質中熱-滲-力耦合現象的核素遷移控制方程，并針對試驗資料，使用所開發的有限元程序對一個核廢料處置概念庫近場的熱-滲-力耦合過程進行了數值模擬；2009年，西南科技大學的李寧波、易發成等人采用間歇法研究在不同環境條件下(濃度、固液比、溫度、pH值、介質)，新疆阿爾泰膨潤土對Sr2+、Cs+的吸附性能影響；2011年，西南科技大學的易發成、周敏娟等人對新疆阿爾泰膨潤土土水特征進行室內試驗研究，結果表明新疆阿爾泰膨潤土具有極強的吸水性和膨脹性，但其持水性能低于內蒙古高廟子膨潤土，其原因可能與塑性指數大小有關，對于體積受吸力變化影響較顯著的膨潤土，宜用重力含水量表達其土—水特征曲線；2012年[14]，陳正漢等人對混合緩沖/回填材料在高溫-高壓-高吸力條件下的力學特性及其應用進行了的研究。目前國內關于膨潤土中核素遷移的研究數據不多，還需要做大量的理論研究與試驗工作。

3 結語與展望

縱觀國內外研究的現狀及進展，對膨潤土及其集成混合材料作為緩沖/回填材料的研究主要表現在以下幾點：①至20世紀80年代末，世界上核能利用較為先進的國家在膨潤土作為緩沖/回填材料的選材方面的工作均已基本完成；②由于高壓實膨潤土非飽和滲透自身的復雜性，導致膨潤土作為緩沖/回填材料的滲透性研究方面仍有待深入；③在膨潤土作為緩沖/回填材料的熱傳導特征研究方面，許多研究機構對添加石英顆粒的大小和數量對熱一水一力學性能的影響正在進行研究；④深地質處置中，化學特征對高壓實膨潤土緩沖性能的影響是極其復雜的，而國內外對此研究也大多處在不成熟階段；⑤能全面反映熱-滲-力耦合作用條件下，高壓實膨潤土的行為特征的多場耦合研究成果仍然有限，還有待深入進行；⑥由于我國高放廢物深地質處置方面的研究工作總體起步較晚，除基本完成膨潤土的選材工作外,總的還處在膨潤土基本物理化學性質、土水特性、熱傳導、滲流的研究階段；⑦在熱、滲、力多場耦合條件下核素在膨潤土中遷移規律研究方面可參考的研究成果還不多見。

[1] The Disposal of Canada’s Nuclear Fuel Waste: Engineered Barriers Alternatives. 1994, AECL-1078COG-93-8

[2] Japan Nuclear Cycle Development Institute. Second Progress Report on Researeh and Development for the Geological Disposal of HLW in Japan Support Report 2 Repository Desing and Engineering Technology.2000，JNC TN 1410 2000-003

[3] Hideo Komine，Ogata N. Prediction for swelling characteristics of compacted bentonite. Canadian Geotechnical Journal, 1996, vol.33(1):11-22

[4] CUI Y.J., DELAGE P., ALZOGHBI P. 2003. Retention and transport of a hydrocarbon in a silt. Géotechnique, Vol. LIII, No 1, 83-91

[5] Tong F.G., Jing L. The Interim Report on Buffer BMT test studies. Work Package 4: Coupled Processes in Buffer and near Rock-Buffer Interfaces. THERESA Project, 2008

[6] Valyashko VM, editor. Hydrothermal Properties of Materials. Experimental data on aqueous phase equilibria and solution properties at elevated temperatures and pressures. Wiley, London, 2009

[7] 王駒，陳偉明，蘇銳，等．高放廢物地質處置及其若干關鍵科學問題[J]．巖石力學與工程學報，2006，25(4)：801-812

[8] 劉月妙，溫志堅．用于高放射性廢物深地質處置的粘土材料研究[J]．礦物巖石，2003，23(4):42-45

[9] 葉為民, 黃雨, 崔玉軍,等. 自由膨脹條件下高壓密膨脹粘土微觀結構隨吸力變化特征[J]. 巖石力學與工程學報,2005，24(24):4570-4575

[10] 周抗寒，李國鼎，俞珂. 多功能膨脹滲透儀的研制及壓實膨潤土的滲透性試驗[J].環境科學，1995, 16(5)：38-42

[11] 沈珍瑤，李國鼎．高壓實膨潤土熱濕藕合效應試驗室模擬[J]．水文地質工程地質， 1999，(1):11-12

[12] 易發成，楊劍，李玉香，等．四川三臺膨潤土的基本物化性能及其對鉛、鈷的吸附性研究[J]．礦物巖石，2004,24(4)：39-42

[13] 張玉軍．核廢料處置概念庫近場熱-水-應力耦合模型及數值分析[J]．巖土力學，2007，28(1)：17-25

[14] 陳正漢，秦冰．非飽和土的應力狀態變量研究[J]．巖土力學，2012，33(1)：1-11