膨潤土在高放廢物處置庫緩沖材料中的研究進展

謝朝新 黃祺臨 秦冰

摘 ? ? ?要：在高水平放射性廢物（簡稱高放廢物）深地質處置庫中，由高壓實膨潤土所組成的人工屏障，是深地質處置庫多重屏障系統重要的組成部分。本文從緩沖材料的篩選要求出發，對比了各國膨潤土的理化性質，并對各國膨潤土的持水性能、滲透性能、熱學性能、膨脹性能以及化學作用和熱-水-力-場耦合等方面進行了綜述，最后結合工程實際簡要分析了當前室內試驗及數值模擬存在的局限性，最后對未來的膨潤土研究趨勢進行了展望。

關 ?鍵 ?詞：高放廢物處置;膨潤土;滲透;多場耦合

中圖分類號：TV443+.3 ? ? ? 文獻標識碼： A ? ? ? 文章編號： 1671-0460（2020）10-2281-05

Abstract： In the deep geological repository of high-level radioactive waste， the artificial barrier composed of high-pressure solid bentonite is an important part of the multi barrier system of deep geological repository. In this paper， the physical and chemical properties of bentonites in different countries were compared based on the screening requirements of buffer materials， and the water holding capacity， permeability， thermal performance， expansion performance， chemical action and thermal water force field coupling of bentonites in different countries were summarized. Finally， the limitations of current laboratory tests and numerical simulation were briefly analyzed in combination with engineering practice， and the future research trend of bentonite was prospected.

Key words： High level radioactive waste disposal; Bentonite; Infiltration; Multi field coupling

高水平放射性廢物（簡稱高放廢物）是一種極為危險的特殊廢物，放射性強、發熱量大、毒性大、半衰期長，對其進行最終安全處置難度極大[1]。隨著我國核電事業的飛速發展，產生的高放廢物將會大量累積，根據中國核電發展規劃，我國2020年規劃中核電站（5 800萬kW）全壽期將產生83 000 t乏燃料，根據美國、瑞典、日本等國家核能發展的情況可知，當前高放廢物如何進行安全處置已成為制約各國核能繼續發展的關鍵因素之一，因此高放廢物的安全處置現已成為關系到環境安全、公共安全及核能可持續發展的重大問題，是一個與核安全同等重要的問題，必須進行研究并加以妥善解決[2]。

當前“深地質處置”已被世界公認為高放廢物處置的安全且可行方案，其核心目的在于阻滯有害核素遷移，具體做法為將高放廢物進行減容和穩定處理后封閉在罐體內，而后放入距離生物圈500～100 m的安全處置場所設施中，永久與人類的生存環境隔離，使得未來可能因安全處置漏洞所產生的核素泄露引起的照射劑量不超過國家規定的劑量限制，埋藏高放廢物的地下工程即稱為高放廢物處置庫[3-4]。

當前世界先進國家均采用“多重屏障系統”的設計理念對高放廢物處置庫進行規劃。如圖1所示，從內到外依次是減容固化體、罐體廢物容器、緩沖回填材料以及處置庫天然圍巖。其中減容固化體、罐體廢物容器、緩沖回填材料三者為人工屏障，處置庫天然圍巖則為自然屏障，緩沖/回填材料作為多重屏障系統中最后一層人工屏障，在阻滯有害核素向生物圈遷移中扮演了極其重要的角色。

1 ?緩沖/回填材料的篩選

作為多重屏障系統最后一道人工屏障，緩沖/回填材料處于罐體容器和處置庫圍巖之間，緩沖/回填材料的選取標準對于高放廢物能否進行安全處置意義重大。緩沖/回填材料作用可概括為：①工程屏障作用。減輕圍巖壓力對廢物容器的作用，使得廢物容器保持在處置孔中心處。②水力學屏障作用。密封高壓實的膨潤土塊吸水后土體會膨脹變形，使得廢物容器與圍巖之間的縫隙及圍巖中的孔隙被充分填充，阻滯地下水流向廢物容器。③化學屏障作用。阻滯核素以放射性氣體或水溶化合物的形式從廢物容器向外遷移，避免容器的氧化侵蝕。④導熱作用。將核廢料產生的衰變能及其轉化的熱量及時傳導出去。

膨潤土是一種以蒙脫石為主要成分的黏土礦物，其遇水膨脹劇烈，自封閉性好，比表面積大，陽離子交換能力強，吸附能力強，滲透性極低，熱傳導性良好，且高壓實膨潤土吸水后不僅會膨脹變形，更能造成土中吸力值變化，使得施工縫隙及圍巖孔隙被充實，最終對圍巖中地下水的滲入形成阻障，因此膨潤土被許多國家選作高放廢物地質處置庫的緩沖材料基材[5-8]。

2 ?各國膨潤土的理化性質

當前，膨潤土作為目前國際上首選的緩沖/回填基質材料，代表性的有日本的Kunigel-V1[9]（鈉基）、加拿大的Avon seal（鈉基）、法國的FoCa[10]（鈣基）、美國的MX-80[11]（鈉基）、西班牙的FEBEX（鈣基），中國則通過系統的選取，將內蒙古高廟子膨潤土推選為我國高放廢物處置庫的首選緩沖/回填材料。各國膨潤土主要物理化學性質如表1所示，下面將以理化性質為切入點對美國、西班牙、中國膨潤土詳細介紹。

2.1 ?美國MX-80（鈉基）

美國的高放廢物處置庫采用美國懷俄明州的MX-80膨潤土，MX-80膨潤土相對密度為2.76，液體限度為52%，塑料限度為42%，利用X射線衍射（XRD）對MX-80膨潤土的礦物學組成進行了定量分析，其主要礦物有蒙脫石、石英、長石、方解石、菱鐵礦和黃鐵礦，其中蒙脫石質量分數高達79%，其化學計量公式為 （Si3.95Al0.04） （Al1.52Mg0.26-Fe0.17Ⅲ） Na0.18Ca0.11-O10-（OH）2[12-13]。

2.2 ?西班牙FEBEX（鈣基）

西班牙高放廢物處置庫采用了FEBEX膨潤土，它被ENRESA（西班牙放射性廢物管理機構）選為適合HLW儲存庫回填和密封的材料。FEBEX膨潤土材料的加工過程包括：分解、輕磨、60 °C烘干、 ? ? 5 mm篩分。

在FEBEX項目的大規模試驗中，膨潤土塊是由具有吸濕性含水量的顆粒黏土在干密度接近 ?1.7 Mg·m-3的情況下進行單軸壓實而制成的。FEBEX膨潤土蒙脫石質量分數高于90%。此外，它還含有不同數量的石英、斜長石、鉀長石、方解石和蛋白石。陽離子交換容量（CEC）在 ? ? ? ?96～102 meq·（100 g） -1之間，主要可交換陽離子為：Ca2+ （35～42 meq· （100 g） -1）、 Mg2+ （31～32 meq·?（100 g） -1）、Na+（24～27 meq·（100g） -1）和K+（2～3 meq·（100g） -1）。膨潤土的液限是102 ± 4%，塑性極限是53 ± 3%，相對密度2.70 ± 0.04，吸濕含水量與實驗室環境達到平衡時為13.7 ± 1.3%。BET法測得的外比表面積為32 ± 3 m2·g-1，吸水性法測得的總比表面積約為725 m2·g-1。飽和FEBEX黏土的水電導率與干密度呈指數關系。密度在1.6 g·cm-3左右時，對去離子水的滲透值約為10～13 m·s-1 [14]。

2.3 ?中國高廟子膨潤土（鈉基）

中國高放廢物處置庫采用了內蒙古高廟子膨潤土作為緩沖/回填材料，國內外學者對其各類基本性質進行了實驗測定，當前各實驗室廣泛應用的高廟子膨潤土指標的相對密度為2.66，天然含水率10.53%，比表面積570 m2·g-1，Na+和Ca2+是主要交換陽離子，塑性指數為275[15-16]。劉月妙[17]等通過X射線衍射分析、化學全分析、紅外光譜分析對天然高廟子膨潤土成分進行了測定，高廟子膨潤土主要成分為蒙脫石、方石英、絲光沸石和正長石，其中蒙脫石質量分數高達63.77% ～ 80.92%。

3 ?膨潤土的性能研究

3.1 ?持水性能

在高放廢物處置庫中的緩沖/回填材料在前期都是處于非飽和狀態，膨潤土在非飽和階段的持水性能可通過持水曲線進行表征。沈珍瑤[18]等采用蒸汽平衡法測定了高吸力下膨潤土的持水曲線，得到了土樣的始干密度和初始含水率對最終含水率沒有影響，并推測此與膨潤土中的微裂隙有關。孟德林[19]等還通過濾紙法和壓力板法測量了脫濕過程中砂-膨潤土混合物的持水曲線，研究持水曲線與孔隙比和配合比之間的關系。結果表明，在同一孔隙比下砂-膨潤土混合物的持水曲線隨著膨潤土的比例增加而向右上方移動，即混合物的進氣值隨著膨潤土的比例增加而增大;另外，含砂量以及孔隙比相同時，膨潤土與福建砂的混合物的持水曲線與日本Kunigel-V1膨潤土與豐浦砂的混合物的持水曲線非常接近。DELAGE[20]等針對高壓實非飽和FoCa膨潤土進行了一系列試驗研究，得到了重現性良好的持水曲線。膨潤土體積隨試樣吸力發生變化，變化過程中氣體體積維持恒定，體積變化主要是由土中水和具有活性的膨潤土之間發生復雜理化反應引起的。

由于高放廢物處置庫中核廢料將長期保持衰變狀態，使得產生熱量致使環境溫度升高，因此，溫度對于緩沖/回填材料的持水性能也有較大影響。沈珍瑤[21]等在不同溫度的條件下對膨潤土的持水特性進行了試驗，發現溫度對膨潤土的持水特性有重要影響。葉為民[22]等在控制溫度分別為20、40、80℃的條件下，測定了高廟子膨潤土的持水曲線。結果表明，某一溫度下，在高吸力段，側限狀態和無側限狀態下膨潤土的持水能力基本相同;而在低吸力段，側限狀態下膨潤土的持水能力低于無側限狀態下的持水能力。秦冰[23]等以Le Chatelier原理為基礎，對高吸力段持水曲線的溫度效應進行了研究，建立了土吸附水的吸附熱有關的溫度效應模型，并經過分析及與已有成果的對比，指出該模型可預測全吸力段持水曲線的溫度效應。

目前針對持水特性的研究，主要停留在體積是否恒定、干密度、含水率等方面，并且結合微觀結構對其持水特性變化進行分析的研究僅僅從定性方面出發，沒有數學建模進行定量分析，特別針對非側限狀態施加不同溫度下的GMZ膨潤土持水特性研究不全面，且膨潤土的持水特性與其微觀結構特征密切相關，輔助微觀結構可對其進行更系統的分析。

3.2 ?滲透性能

在高放廢物處置庫中，土體的滲透性能可由其滲透系數進行表征，滲透系數反映了緩沖/回填材料阻滯地下水及氣體泄露擴散的性能。Castellanos[24]等通過實驗得出，干密度是影響滲透性能的重要因素，滲透系數隨干密度增大而減小。劉月妙[25]等通過實驗得出，試驗溫度會影響膨潤土滲透性能，滲透系數會隨試驗溫度的升高而增大，溫度由80 ℃到90 ℃時，則有所降低;同時膨潤土類型不同，滲透性能存在差異，但均能夠滿足處置庫低滲透性的要求。SUZUK [26]等通過實驗得出，入浸液的類型、溶液濃度、pH值等均會影響膨潤土滲透性能，溶液濃度越高，pH值越大，則滲透系數越大。

目前，針對飽和、非飽和土的滲透性能以及相關的影響因素方面的研究已描述了其大致的變化規律，但目前對膨潤土滲透性能的研究儀器均為自主研制，且因膨潤土滲透系數低導致試驗時間較長，宏觀實驗的可重復性較差，其測量結果受到選擇不同的試驗終止標準可能出現不同的測量結果，尤其是在滲濾液濃度較低的情況下更是如此。

3.3 ?熱學性能

GMZ膨潤土具有良好的熱穩定性，一般不會因為在高溫條件下結構發生破壞而影響其工程屏障特性。劉曉東[27]等進行的差熱分析表明，至少要 ?870 ℃才能使其晶格發生破壞。同時作為工程屏障，GMZ膨潤土必須具有良好的導熱性，保證處置庫內外不會因為過大的溫度梯度從而發生破壞。

盡管膨潤土一般不會因為高溫條件對其結構進行破壞，但高溫條件下對于膨潤土的工作性能也會產生一定影響[28]，緩沖回填材料的導熱系數能夠表征高放廢物的衰變熱能是否有效傳遞，因此各國針對導熱系數展開了大量研究。朱國平[29]等通過NK-III100E型雙試件熱導儀測試了GMZ膨潤土的導熱系數，發現膨潤土初始含水率小于塑限時，隨著含水率的增加，導熱系數有增加的趨勢。徐云山[30]等以高廟子膨潤土以及MX-80膨潤土為研究對象，采用熱探針法和KD2Pro型熱性分析儀測定了不同溫度下試樣的熱傳導系數，結果表明二者熱傳導系數均隨溫度增大而增大。OULD-LAHOUCINE等[31]在Kunigel-V1膨潤土中添加不同比例的石英砂配置成混合材料，采用了熱探針法對該混合材料的導熱系數進行了分析，最后運用了數值模擬的預測方法對其適用性進行了研究。TANG[32]等采用了熱線法對MX-80膨潤土的導熱系數進行了研究，分析了干密度、含水率、飽和度、物質成分等因素的影響，得出了MX-80膨潤土中石英質量分數越高，導熱系數越大的結論。劉月妙[33]等研究了將石英砂與石墨材料混合進入的膨潤土的導熱性能進行了研究，結果表明當石英砂與石墨添加適量時可明顯提升膨潤土的導熱、熱擴散性能。謝敬禮[34]等以石英砂、北山花崗巖碎屑組成的混合材料為添加劑，實驗研究了不同質量分數的北山花崗巖碎屑及石英砂對高廟子膨潤土導熱性能的影響，發現兩者均可在不同程度上提升膨潤土的導熱系數，石英砂作用優于北山花崗巖碎屑。

3.4 ?膨脹性能

膨潤土膨脹原理主要有2點：一是晶層膨脹（Crystalline Swelling），二是擴散雙電層理論（Diffuse Double Layer Theory，DDL）。目前，針對膨潤土膨脹力的試驗方法主要有2種：一是加壓平衡法，即通過不斷增大軸向荷載，以平衡試樣膨脹形變;二是恒體積法，即維持試樣體積不變，通過力傳感器量測試樣的膨脹應力變化。深地質處置中緩沖/回填材料總體積變化程度不大，因而恒體積法更符合工程實際。

劉月妙[35]等對GMZ膨潤土分別開展了膨脹力、無荷和有荷膨脹率試驗，發現膨脹力受干密度控制，即試樣的膨脹力隨干密度的增加而增加，同時發現，增加軸向荷載能夠抑制試樣膨脹，相同初始干密度條件下，軸向荷載越大最終膨脹率越小。葉為民[36]等對高壓實高廟子膨潤土膨脹力特性展開了研究，實驗采用了恒體積試驗法，結果表明，膨脹力發展過程曲線與吸水量曲線有明顯的階段性特征，膨脹力隨時間的變化曲線是一條漸近線，并且時間/膨脹力與時間之間存在很好的線性關系，同時發現干密度與膨脹力之間存在指數關系。葉為民[37]等還根據雙電層理論對膨潤土體變特征的預測進行了研究，得到了低吸力范圍內高壓實GMZ膨潤土自由膨脹率的預測方法;并提出采用線性延伸法估算殘余孔隙體積可保證所求膨潤土的孔隙率或孔隙體積更為準確。

3.5 ?熱-水-力-場耦合

在地質處置庫中，膨潤土處于復雜的熱-水力場耦合條件下工作，在多場耦合研究方面，國際上開展了一系列不同尺度的熱-水-力耦合大型室內模型實驗和地下實驗室原位實驗，如日本 BIG-BEN 實驗、西班牙FEBEX Mock-up實驗、比利時 OPHELIE Mock-up實驗、捷克Mock-Up-CZ實驗、FEBEX原位實驗、日本Kamaishi mine原位實驗、瑞典 TBT 原位實驗等。

劉偉[38]等開展了熱-水耦合實驗，結果表明含水率、干密度對溫度場平衡有影響;曹勝飛[39]等對高廟子膨潤土進行了熱-水-力耦合實驗，研究表明隨著濕度及加熱器溫度增大，膨潤土的溫度也在增大。張志紅[40]等首先基于能量、質量守恒分析了單一溶質在飽和土體中運移的熱-水-力-化全耦合的溫度效應，研究表明溫度可顯著影響溶質遷移過程，熱擴散作用加速溶質運移而熱滲透、熱固結對其運移有阻滯作用。吳樹高[41]基于緩沖材料熱-水-力三場耦合理論通過有限元數值分析實現了三場耦合方程組的數值求解。

現有的熱-水-力場研究缺乏對其熱、水、力共同作用下的宏微觀多尺度行為的系統研究，特別是缺乏對其微觀結構變化的系統試驗研究和合理定量描述，因此關于熱-水-力場的耦合作用研究，應更多關注在建立微觀結構變化與膨潤土實際工作性能的聯系，分析其演化規律，并建立相應模型以深化對膨潤土實際工作性能與安全狀態的認識。

4 ?總結與展望

當前，緩沖/回填材料篩選標準已經大致明確，并且各國在室內及地下做了大量試驗，已做出其處置庫的安全評價，但緩沖/回填材料的試驗仍存在一些問題。

1）在對緩沖/回填材料的選擇階段，應結合各國高放廢物處置庫相應的地質條件對膨潤土進行篩選，其篩選條件應不僅局限于高放廢物處置庫的通用標準要求，更應取決于廢棄物類型以及各國的高放廢物處置理念，以此確定所需膨潤土的物化性質再進行實驗。

2）當前膨潤土的實驗器材大多為自主研發，各個實驗室的儀器尺寸大小、實驗控制條件、實驗方案確定均缺少統一的標準進行約束，導致宏觀實驗的可重復性較差。

3）當前針對膨潤土熱學性能的室內試驗大多處于穩定的環境溫度下，但在工程實際中，緩 ? ? 沖/回填材料靠近罐體一側長期處于 100 ℃左右的高溫環境，但與天然巖體接觸一側則為20 ℃左右的自然問題，其環境溫度并不均勻，因此未來對于膨潤土熱學性能研究不僅僅要停留在穩定的環境溫度下進行實驗，更應結合工程實際使其整體處于非穩定的溫度梯度下進行實驗。

4）當前對于高廟子天然鈉基膨潤土的膨脹性、滲透性、導熱性、壓實特性、化學緩沖特性、添加劑、熱-水-力耦合長期性能等各項性質開展了大量室內試驗和數值模擬研究，但缺乏對于宏微觀多尺度行為的系統研究，特別是缺乏對其微觀結構變化的系統試驗研究和數值模擬手段下合理的定量描述。

總之，膨潤土性能研究仍有待進一步探索，隨著研究的不斷深入和對其他相關學科成果的借鑒，我國高放廢物處置的緩沖回填材料技術將會不斷邁上新的臺階。

參考文獻：

[1]WAND J， CHEN L， SU R， et al. The Beishan underground research laboratory for geological disposal of high-level radioactive waste in China： Planning， site selection， site characterization and in situ tests[J].Journal of Rock Mechanics and Geotechnical Engineering. doi;10. 1016/j.jrmge. 2018. 03. 002.

[2]王駒. 中國高放廢物地質處置21世紀進展[J]. 原子能科技技術， 2019（10）：2072-2082.

[3]PUSCH R. Geological storage of highly radioactive waste[M]. Verlag Berlin Heidelberg： Springer， 2008.

[4]陳經明，李尋，羅躍，周澤超. 高放廢物處置庫中膨潤土性質研究進展[J]. 科技與創新，2019（13）：1-2.

[5]ENRESA. Full-scale engineered barriers experiment for a deep geological repository for high level radioactive waste in crystalline host rock[R]. Madrid， 2018.

[6]NAGRA. Project opalinus clay safety report： demonstration of disposal feasibility for spent fuel， vitrified high-level waste and long-lived intermediate-level waste[R]. Wettingen， 2017.

[7]NWMO. Choosing a way forward： The future management of Canadas used nuclear fuel[R]. Toronto， 2015.

[8]王月超，王騰達，楊雙春，等. 膨潤土改性及其在鉆井液方面的應用.[J].當代化工，2016，45 （1）：73-74.

[9]FERN?NDEZ R， RUIZ A I， CUEVAS J. The role of smectite composition on the hyperalkaline alteration [J]. Applied Clay Science， 2014（95）： 83-94.

[10]VILLAR M V， LLORET A， Influence of dry density and water content the swelling of a compacted bentonite [J]. Applied Clay Science， 2018， 39（1）： 38-49.

[11]WANG Q， TANG A M，CUI Y J， et al. Experimental study on the swelling behavior of bentonite/claystone mixture[J]. Engineering Geology， 2018， 124（1）： 38-49.

[12]Montes H G， Duplay J， Martinez L， et al. Influence of interlayer cations on the water sorption and swelling-shrinkage of MX-80 bentonite[J].Appl. Clay Sci.， 2013，23：309-321.

[13]BRADBURY M H， BAEYENS B. Predictive sorption modelling of Ni（II）， Co（II）， Eu（III）， Th（IV） and U（VI） on MX-80 bentonite and opalinus clay： a ‘‘bottom-up approach[J].Appl. Clay Sci. ， 2011，52： 27-33.

[14]LLORET A， VILLAR M V. Advances on the knowledge of the thermo-hydro-mechanical behavior of heavily compacted ‘‘FEBEX bentonite[J]. Phys. Chem. Earth ， 2017，32 （8-14） ：701-715.

[15]He Y， Chen Y G， Ye W M， Equilibrium， kinetic， and thermodynamic studies of adsorption of Sr（II） from aqueous solution onto GMZ bentonite [J]. Environ. Earth Sci.， 2016，75 ：807.

[16]YE W M， CHEN Y G， CHEN B， et al. Advances on the knowledge of the buffer/backfill properties of heavily-compacted GMZ bentonite [J]. Engineering Geology， 2010 （116）： 12-20.

[17]劉月妙，陳璋如.內蒙古高廟子膨潤土作為高放射性廢物處置庫回填材料的可行性[J]. 礦物學報，2011，21（3）：541-543.

[18]沈珍瑤，李國鼎，李書紳.高壓實膨潤土水分特征曲線的測定 [J]. 工程勘察，2006（4）： 27-28.

[19]孟德林，孫德安，劉月妙.高廟子膨潤土與砂混合物的土-水特征曲線[J]. 巖土力學，2012，33（2）：509-514.

[20]HOLTON D ， MYERS S ， CARTA G ， et al. Application of a novel approach to assess the thermal evolution processes associated with the disposal of high-heat-generating waste in a geological disposal facility[J]. Engineering Geology， 2016， 211：102-119.

[21]沈珍瑤，程金茹，馬炳輝. 鹽溶液法測定非飽和膨潤土的土-水特征曲線[J]. 水文地質工程地質，2001（4）： 6-8.

[22]葉為民，申淼，陳寶，等. 溫控高壓實GMZ01膨潤土-砂混合物土水特性[J]. 工程地質學報，2013，21（3）： 385-390.

[23]秦冰，陳正漢，孫發鑫，等. 高吸力下持水曲線的溫度效應及其吸附熱力學模型[J]. 巖土工程學報，2012，34（10）： 1877-1886.

[24]CASTELLANOS， RUIZ A I， CUEVAS J. The role of smectite composition on the hyperalkaline alteration [J]. Applied Clay Science， 2017（95）： 83-94.

[25]劉月妙，王駒，曹勝飛，等.中國高放廢物處置緩沖材料大型試驗臺架和熱-水-力-化學耦合性能研究[J]. 巖土力學，2013，34（10）： 2756–2762.

[26]SUZUK， CHEN B， Thermal conductivity of compacted bentonite as a buffer material for a high-level radioactive waste repository[J]. Annals of Nuclear Energy， 2016，94： 848-855.

[27]劉曉東，羅太安，朱國平，等.緩沖/回填材料內蒙古高廟子膨潤土性能研究[J]. 中國核科技報告，2007（2）： 140-156.

[28]CHEN Z G， TANG C S， SHEN Z， et al. The geotechnical properties high-level of GMZ buffer/backfill material used in radioactive nuclear waste geological repository： a review[J]. Environmental Earth Sciences， 2017， 76（7）： 270-285.

[29]朱國平，劉曉東，羅太安. 內蒙古高廟子膨潤土熱學性能研究[J]. 資源環境與工程，2015，19 （3）：220-222.

[30]徐云山，孫德安，曾召田，等. 膨潤土熱傳導性能的溫度效應[J].巖土力學，2020，41（1）：39-45.

[31]OULD-LAHOUCINE C，SAKASHITA H，KUMADA T. Measure- ment of thermal conductivity of buffer materials and evaluation of existing correlations predicting it[J]. Nuclear Engineering and Design， 2012， 216（1）： 1-11

[32]TANG， RAO D Y， XU T， et al. SPH-FDM boundary for the analysis of thermal process in homogeneous media with a discontinuous interface[J]. International Journal of Heat&Mass Transfer， 2018， 117： 517-526.

[33]劉月妙，蔡美峰，土駒. 高放廢物處置庫緩沖材料導熱性能研究陰[J].巖石力學與工程學報，2017，26 （增刊2）：3891-3896.

[34]謝敬禮，劉月妙，曹勝飛，等. 高廟子膨潤土及其混合材料導熱性能研究[J].地下空間與工程學報，2019，19（6）：1719-1806.

[35]劉月妙，黃立明，葉為民，等.GMZ01 緩沖/回填材料在柔性邊界條件下的膨脹特性[C].第四屆廢物地下處置學術研討會論文集，2018.

[36]葉為民， SCHANZ T，錢麗鑫，等.高壓實高廟子膨潤土 GMZ01 的膨脹力特征[J]. 巖石力學與工程學報，2017，26（2）： 3861-3865.

[37]葉為民，黃偉，陳寶，等. 雙電層理論與高廟子膨潤土的體變特征[J].巖土力學，2019，30（7）：1900-1903.

[38]劉偉，楊仲田，李公平，等. 高壓實緩沖材料熱-濕耦合效應研究[J]. 輻射防護，2012，32（5）：277-282.

[39]曹勝飛，劉月妙，喬蘭，等. 膨潤土熱-水-力耦合性能試驗研究[C]. 第四屆廢物地下處置學術研討會， 2018.

[40]張志紅，韓林，田改壘，等. 飽和土體熱-水-力-化全耦合一維溶質運移模型[J].東南大學學報，2019，49（6）：1001-0505

[41]吳樹高.高放廢物處置庫含接縫緩沖材料砌塊熱-水-力耦合性能研究[D].北京：北京交通大學，2017.