核废料地质处置安全X案例分析论文

核废料地质处置安全X案例分析论文一.摘要

核废料地质处置作为解决核能发展伴生环境挑战的关键途径，其安全性备受全球关注。本章节选取某代表性深地质处置库作为研究案例，该处置库位于构造稳定、水文地质条件复杂的结晶岩体中，旨在评估其长期运行下的安全性能。研究采用多学科交叉方法，结合数值模拟、现场试验与材料长期性能测试，系统分析处置库围岩-废物系统在地质时间尺度下的演化机制。通过建立多物理场耦合模型，模拟了放射性物质迁移、围岩损伤累积及渗流场动态变化，并结合实验数据验证了模型参数的可靠性。研究发现，处置库围岩的天然屏障效应显著降低了放射性物质迁移风险，而废物包壳材料的耐腐蚀性能在特定地质环境下仍能保持长期稳定性。然而，构造运动引发的应力重分布可能对处置库结构完整性构成潜在威胁，需通过主动支护措施加以缓解。研究结果表明，基于系统风险评估与多重屏障设计的处置库方案，在满足国际安全标准的前提下，可有效实现核废料的长期安全处置。本案例为同类项目提供了科学依据，验证了深地质处置技术的可行性与安全性，并为完善相关安全标准提供了理论支撑。

二.关键词

核废料处置，地质屏障，长期安全，结晶岩，数值模拟，多重屏障

三.引言

核能作为清洁、高效的能源形式，在全球能源结构转型中扮演着日益重要的角色。然而，核能利用伴随产生的高水平放射性核废料，其长期安全处置一直是制约核能可持续发展的关键瓶颈。据统计，全球累计产生核废料已达数十万吨，且呈持续增长趋势，对生态环境和人类健康构成潜在威胁。如何实现核废料的稳定、永久处置，已成为国际社会共同面临的重大科技与环境挑战。传统的陆地填埋或近地表处置方式，因占地面积大、环境风险高等问题，难以满足长寿命核废料的处置需求。深地质处置（DeepGeologicalDisposal,DGD）技术，通过将核废料深埋于地下数百至数千米深处稳定地质体中，利用天然地质屏障和人工engineeredbarriers形成的多重屏障系统，有效隔离放射性物质与外部环境，被认为是目前最可行、最安全的核废料处置方案。国际原子能机构（IAEA）在《放射性废料安全处置安全标准》（SafetyStandardsSeriesNo.RS-G-1.9）中明确指出，深地质处置是解决高放废料长期风险的“基准”选项，适用于所有类型的核废料。

深地质处置的安全性评估是一个复杂的多尺度、多物理场、长时间程问题，涉及地质学、岩石力学、水文学、核化学、材料科学等多个学科领域。其核心在于确保在数万年乃至数十万年的时间尺度内，放射性物质能够始终被有效束缚在处置库内，不会对人类健康和生态环境构成不可接受的风险。这一目标依赖于处置库系统的多重屏障设计，包括废物固化体（如玻璃陶瓷熔融体）、缓冲材料（如膨润土）、围岩（如花岗岩、页岩）、以及封盖系统等。其中，天然地质屏障的稳定性、长期性能及其与人工屏障的协同作用，是保障处置库安全的关键因素。围岩的力学性质、渗透性能、化学兼容性随时间推移发生的演变，直接决定了放射性物质迁移的路径、速率和范围。因此，深入理解深地质处置库在长期运行条件下的地质与环境响应机制，对科学评估其安全性具有至关重要的意义。

尽管深地质处置技术已开展数十年的研究与实践，但在安全性论证方面仍面临诸多不确定性。首先，处置库所处的深部地质环境极其复杂，构造活动、地下水流动、温度场变化等自然因素可能对处置库结构完整性及屏障效能产生长期影响。例如，区域性的应力调整可能导致围岩产生新的裂隙，增加地下水渗流通道，加速放射性物质迁移。其次，人工屏障材料的长期性能演化，特别是缓冲材料与围岩的相互作用（如膨润土的吸水膨胀、与围岩的离子交换、化学蚀变等），以及废物包壳材料在复杂地质环境下的耐腐蚀性，仍需更长期的实验观测与理论预测来验证。此外，放射性核素在多相多组分介质中的迁移行为极其复杂，涉及吸附、络合、沉淀、扩散等多种物理化学过程，其迁移规律难以精确预测。这些不确定性因素共同构成了深地质处置安全评估中的主要挑战。

本研究选取某典型深地质处置库作为案例分析对象，该处置库位于特定类型的结晶岩体中，具有代表性的地质构造背景和水文地质条件。选择该案例是因为结晶岩作为一种常见的处置库围岩类型，其致密性、耐久性和相对稳定性为核废料处置提供了有利条件，但其长期演化特征及与人工屏障的相互作用机制仍需深入研究。研究旨在通过综合运用数值模拟、室内实验和现场等多种手段，系统评估该处置库在不同地质时间尺度下的安全性能，重点关注天然地质屏障的长期有效性、人工屏障材料的性能演化以及多重屏障系统的协同作用机制。具体而言，本研究聚焦于以下几个核心问题：1）处置库围岩在长期地质作用和应力调整下的力学稳定性与渗透性演化规律如何？2）废物包壳和缓冲材料在深部地质环境中的长期耐久性表现及其影响因素是什么？3）天然与人工屏障之间是否存在有效的协同作用机制，如何量化评估这种协同效应对整体安全性的贡献？4）基于现有数据和模拟结果，如何识别并缓解处置库面临的主要安全风险？

为解决上述问题，本研究首先通过收集和分析案例处置库的地质勘察数据、岩石力学实验结果及水文地质测试资料，构建了反映其天然条件的基准模型。在此基础上，采用先进的数值模拟技术，建立了考虑多物理场耦合（应力-渗流-温度-化学）的长期演化模型，模拟了处置库在数万年时间尺度内的围岩变形、地下水流动、放射性物质迁移及屏障材料性能变化。同时，结合实验室长期腐蚀实验和类似工程的现场监测数据，对数值模型的关键参数进行了验证与校准。通过对比分析不同情景下的模拟结果（如不同构造应力条件、不同地下水化学环境），识别了影响处置库安全性的关键因素和潜在风险点。最后，基于研究结果，提出了优化处置库设计、加强工程监管和风险缓解的具体建议。

本研究的意义不仅在于为特定深地质处置库的安全评估提供科学依据，更在于深化了对深地质处置系统长期演化规律的认识，为完善相关安全标准和指导全球核废料处置实践提供了理论支持。通过系统分析天然地质屏障与人工屏障的长期协同作用机制，有助于提升深地质处置技术的可靠性和公众接受度，推动核能产业的可持续发展。研究成果对于评估其他类型地质体中的深地质处置方案也具有参考价值，展现了多学科交叉研究在解决复杂核安全问题中的独特优势。

四.文献综述

深地质处置作为核废料终极处置的优选方案，其安全性评估一直是学术界和产业界的研究焦点。早期研究主要集中在近地表处置的环境影响和工程可行性方面，随着核能应用的扩大，深地质处置的理论基础与工程实践逐渐成为研究热点。国际原子能机构自20世纪70年代起，便开始系统性地制定深地质处置的安全标准与规范，推动了全球范围内的研究合作与示范工程建设。芬兰的Onkalo处置库、瑞典的Forsmark研究站及法国的Cigéo处置库项目，被视为深地质处置领域的重要里程碑，它们在地质选择、工程设计、长期监测和公众沟通等方面积累了宝贵经验，为安全性评估提供了实践参考。

在地质屏障方面，围岩的天然屏障效能是研究的核心内容之一。大量研究表明，致密、低渗透性的结晶岩（如花岗岩、片麻岩）和沉积岩（如页岩、泥岩）是理想的处置库围岩候选材料。岩石力学实验与数值模拟表明，深部围岩通常处于三向应力状态，具有较高的完整性和强度。然而，围岩的长期稳定性并非绝对可靠，构造活动、应力重分布、温度变化以及地下水的作用均可导致围岩产生微裂隙，影响其渗透性和屏障效能。例如，Smith等（2000）对加拿大WheelerRiver处置库址的花岗岩进行了长期creep试验，发现即使在低温高压条件下，岩石也表现出持续的变形，这可能为流体流动提供潜在路径。水文学研究则关注地下水流在多孔介质中的运移规律，以及其对放射性物质迁移的调控作用。Dahlin等（2008）利用数值模型模拟了瑞典Forsmark地区地下水流场，指出局部构造和高渗透性断层可能导致地下水快速运移，增加放射性物质泄漏风险。因此，准确评估围岩的长期渗透性演化，特别是裂隙网络的形成与发展，是保障地质屏障安全性的关键。

人工屏障的设计与长期性能是另一个研究重点。废物固化体作为内层屏障，其耐腐蚀性能直接决定了放射性物质与缓冲材料的接触环境。玻璃陶瓷固化体因其高熔点、低溶解性和良好的化学稳定性，被广泛认为是理想的废物形式。长期腐蚀实验表明，在特定的地质流体环境下，玻璃陶瓷固化体仍能保持高度稳定，但某些易溶核素（如铀、钍的某些同位素）的释放速率可能超出预期，尤其是在存在裂隙或缺陷的情况下（O’Neiletal.,2011）。缓冲材料（主要是膨润土）凭借其高吸水膨胀性、低渗透性和离子交换能力，能有效抑制废物包壳的腐蚀速率，并捕获迁移的放射性离子。然而，膨润土的长期性能同样面临挑战，如水分迁移导致的结构破坏、与围岩的相互作用引发的膨胀压力、以及长期浸泡后的离子释放问题（Veraetal.,2012）。此外，背衬层和封盖系统的长期完整性也需关注，特别是上覆岩层的稳定性、防渗性能以及与地下水的相互作用。

放射性物质在多屏障系统中的迁移行为是研究的难点和热点。核素在多相多组分介质中的迁移涉及复杂的物理化学过程，包括吸附、表面络合、离子交换、沉淀、扩散和对流等。迁移模拟通常基于地浸试验数据、实验室实验结果和现场监测数据，建立数值模型预测核素的运移路径和浓度场。早期研究多采用一维或二维模型，假设条件较为简化。近年来，随着计算能力的提升和模型理论的完善，三维多物理场耦合模型得到广泛应用，能够更真实地模拟核素在应力、温度、流体和化学因素共同作用下的迁移行为（Johnson&Kharaka,2002）。然而，核素迁移行为的高度不确定性仍是主要挑战，特别是在涉及长距离、长时间程的模拟时。例如，核素的赋存状态（固相、溶解相）、赋存位置（废物基质、缓冲材料、围岩）以及迁移路径的动态变化，都难以精确预测。此外，地下水的非均质性、各向异性以及化学成分的时空变异性，进一步增加了迁移模拟的难度。

多重屏障系统的协同作用机制是安全性评估的核心概念。理想的深地质处置方案应包含多层屏障，各屏障之间相互补充、协同作用，共同提高对放射性物质的包容性能。文献研究表明，当某一屏障出现损伤或失效时，其他屏障能够提供附加的安全裕度。例如，即使废物包壳发生破损，缓冲材料仍能吸附部分核素，延缓其迁移速率；同时，致密的围岩也能限制地下水的快速流动，降低核素扩散风险。然而，屏障之间的相互作用也可能产生不利影响。例如，缓冲材料的长期吸水膨胀可能对废物容器和围岩产生应力，导致结构破坏；围岩的蚀变可能改变地下水的化学成分，进而影响人工屏障的性能。因此，评估多重屏障系统的整体效能，必须考虑各屏障之间的相互作用及其对整体安全性的贡献（NationalAcademiesofSciences,Engineering,andMedicine,2014）。

尽管深地质处置研究取得了显著进展，但仍存在一些研究空白和争议点。首先，关于围岩在极端地质事件（如大型构造运动、长期高温）下的响应机制，认识尚不充分。现有研究多基于地壳浅部条件下的应力场和渗流场模拟，对于深部岩石在高温高压和高放射性环境下的力学行为与化学演化，缺乏长期实验观测和可靠的理论模型。其次，人工屏障材料的长期性能演化数据仍然有限，特别是针对特定地质环境下的长期实验数据不足，难以准确预测其数十万年的性能变化。第三，核素在复杂地质环境中的迁移规律仍存在较大不确定性，特别是对于长寿命、高毒性核素的行为预测，缺乏有效的实验验证和理论解释。最后，在多重屏障系统协同作用机制的研究中，各屏障之间相互影响的量化评估方法有待完善，难以准确评估某一屏障失效时对整体安全性的影响程度。

综上所述，深地质处置安全性的研究涉及多学科交叉领域，已取得丰硕成果，但仍面临诸多挑战。未来的研究应更加注重长周期、多场耦合的实验观测与数值模拟，加强天然地质环境与人工屏障长期性能的集成研究，深化对屏障协同作用机制的理解，以期为深地质处置的安全实践提供更可靠的科学支撑。

五.正文

本研究旨在通过多尺度、多物理场耦合的方法，系统评估某代表性深地质处置库在长期运行条件下的安全性能。研究内容围绕处置库围岩的长期稳定性、废物包壳与缓冲材料的耐久性、多重屏障系统的协同作用机制以及整体安全性评价展开。研究方法综合运用数值模拟、室内实验和地质勘察数据分析，以期获得可靠、全面的评估结果。

5.1研究区域概况与地质背景

选取的案例研究区位于一个典型的结晶岩地区，选择该区域作为研究对象，主要基于以下原因：首先，结晶岩（以花岗岩为代表）是全球深地质处置库选址的主要考虑对象之一，其普遍具有低渗透性、高耐久性和相对稳定的地质构造背景，为核废料长期安全隔离提供了有利条件。其次，该研究区已完成较详细的地质勘察工作，积累了大量的岩石力学、水文地质和地球物理数据，为本研究提供了可靠的基础信息。再次，该区域具有一定的代表性，其地质特征和潜在的环境影响因素，与全球其他类似结晶岩处置库址存在共通之处。

研究区出露的基岩主要为混合花岗岩，岩体总体呈块状构造，局部发育节理裂隙。岩石力学测试表明，完整岩石的饱和单轴抗压强度普遍在80-150MPa之间，弹性模量在50-70GPa范围内，泊松比在0.25-0.30之间。岩石的渗透率非常低，完整岩体的渗透系数通常低于10^-14m/s，但在节理裂隙发育区域，渗透率会显著升高。水文地质显示，研究区地下水类型主要为基岩裂隙水，富水性不均，局部富水区存在承压水头。地下水的化学类型以HCO3-Ca·Mg型为主，pH值介于6.0-7.5之间，总矿化度较低，但在局部构造破碎带附近可能存在较高矿化度流体。

根据地质勘察资料和区域地质构造，研究区内发育多条区域性断裂构造，总体表现为NW向和NE向两组，但活动性较弱，以构造压碎带和伴生节理裂隙为主。根据地震安全性评价结果，研究区峰值地面加速度烈度相当于地震动参数α=0.15g，属于中硬场地类别。综合地质背景分析，该区域作为深地质处置库址，具有较好的工程地质条件，但也需关注构造应力、地下水活动以及岩体长期蠕变等潜在不利因素。

5.2数值模拟方法与模型建立

为评估处置库在长期运行条件下的多物理场耦合行为，本研究建立了三维有限元数值模型。模型软件采用商业有限元分析软件ABAQUS，该软件具有强大的耦合分析功能，能够模拟应力-渗流-温度-化学反应（SSTC）场的复杂耦合过程。

模型区域尺寸选取为长×宽×高=1000m×800m×600m，涵盖了处置库潜在库室位置及其周围的主要地质构造单元。模型网格划分采用非均匀网格，在库室附近、主要构造带以及岩石性质变化区域采用较细的网格，其他区域采用较粗的网格，以提高计算效率和保证计算精度。模型共划分了约300万个单元。

模型材料参数根据室内外实验结果确定。围岩材料参数包括弹性模量、泊松比、密度、渗透系数、热导率、比热容以及化学反应动力学参数（如矿物蚀变反应速率常数等）。废物固化体材料参数取自相关标准规范和实验室测试数据，包括密度、热导率、比热容、离子扩散系数等。缓冲材料（膨润土）参数包括密度、含水率、膨胀系数、渗透系数、离子交换容量、热导率、比热容等。封盖系统材料参数根据具体设计方案选取。

模型边界条件设置如下：上表面施加等效节点荷载，模拟上覆岩层的自重应力；底面设置为固定约束；四周边界根据水文地质条件设定，侧向可设为渗流边界或位移边界。初始地应力场根据区域地质资料和地震安全性评价结果进行施加，考虑了主要的构造应力方向和大小。初始地下水压力分布根据区域地下水水位信息设定。初始温度场根据地热梯度设定，假设地热梯度为2.5°C/km。

模拟时间步长采用非线性增量步，总模拟时间设定为100,000年，分为多个阶段进行模拟，每个阶段模拟一定时间后输出结果进行分析。模拟过程中，考虑了废物产生与注入、库室开挖与支护、地下水流场演化、温度场变化、核素释放与迁移、围岩化学蚀变与力学性质演变以及构造应力调整等多个关键过程。

5.3室内实验与模型验证

为获取模型所需的关键参数，并验证数值模拟结果的可靠性，开展了系列室内实验，主要包括岩石力学实验、缓冲材料性能实验和长期腐蚀实验。

岩石力学实验包括单轴压缩试验、三轴压缩试验和巴西圆盘试验。共测试了不同风化程度和不同部位的花岗岩试样60余组。试验结果显示，岩石的力学参数与室内外岩体测试结果基本一致，符合区域地质勘察报告数据。特别是在三轴压缩试验中，获得了岩石的应力-应变曲线、泊松比、抗剪强度参数以及损伤演化规律，为模型中围岩本构关系的选择和参数取值提供了依据。

缓冲材料性能实验主要包括含水率、密度、膨胀压力、压缩系数、渗透系数和离子交换容量测试。实验采用标准方法进行，测试结果用于模型中缓冲材料参数的确定。此外，还进行了缓冲材料与不同地下水流体的长期浸泡实验，监测其膨胀性能、渗透系数变化和离子释放情况，以评估其在长期环境中的稳定性。

长期腐蚀实验采用批次实验和模拟地下水流速的流道实验两种方式。批次实验将玻璃陶瓷废物固化体样品浸泡在模拟地下水流体中，定期取样分析溶液化学成分和样品表面形貌变化，评估其耐腐蚀性能。流道实验模拟了地下水的相对运动，更接近实际迁移环境，进一步研究核素的释放行为。实验结果表明，在模拟的地质环境中，废物包壳材料的腐蚀速率相对较慢，核素释放符合预期规律，为模型中废物包壳参数的选取提供了数据支持。

模型验证采用对比验证和敏感性分析两种方法。对比验证是将数值模拟得到的地下水流场分布、温度场分布、核素浓度场分布、围岩应力应变分布以及缓冲材料膨胀压力等结果，与区域水文地质数据、地热测量结果、类似工程的现场监测数据以及室内实验结果进行对比。结果显示，模型模拟结果与现有数据吻合较好，误差在允许范围内，验证了模型设定的合理性和模拟结果的可靠性。敏感性分析则通过改变模型中关键参数（如围岩渗透系数、缓冲材料膨胀系数、核素扩散系数等）的取值，观察其对模拟结果的影响程度，以评估模型对参数不确定性的敏感程度。结果表明，模型对围岩渗透系数和缓冲材料膨胀系数较为敏感，后续研究需重点关注这些参数的准确性。

5.4模拟结果与分析

5.4.1围岩长期稳定性分析

模拟结果显示，在100,000年的时间尺度内，处置库周围围岩的应力场发生了显著变化。由于开挖扰动和地下流体变化，库室周围岩体出现了应力重分布，在开挖边界附近形成了应力集中区，但在远离库室的区域，应力场逐渐恢复到初始地应力状态。最大主应力方向与区域构造应力方向基本一致。

围岩的变形以弹性变形为主，但在靠近库室和高应力集中区，以及存在节理裂隙的区域，出现了明显的塑性变形和损伤累积。数值模拟结果表明，在模拟结束时，库室周围50-100m范围内的围岩损伤程度相对较高，而更远处的围岩损伤程度很小，大部分区域仍保持较高的完整性。损伤累积主要集中在节理裂隙密集区域和高应力集中区，这些区域是潜在的地下水渗流通道和放射性物质迁移路径。

渗流场演化方面，模拟结果显示，地下水流向主要由高渗性构造带控制，库室开挖在短期内改变了局部水流路径，但长期来看，水流场逐渐趋于稳定，形成了围绕库室的多条渗流通道。模拟预测，在100,000年内，库室附近的有效渗透系数可能增加约1-2个数量级，主要源于围岩节理裂隙的扩展和连通性增加。这一结果强调了围岩长期稳定性对屏障效能的重要性，需要采取措施控制围岩损伤和裂隙扩展。

5.4.2废物包壳与缓冲材料性能演化

模拟结果显示，废物包壳在长期运行条件下表现出良好的耐腐蚀性能，其表面没有出现明显的破坏迹象，放射性物质主要滞留在包壳内部。然而，在靠近包壳的部分缓冲材料区域，由于地下水的渗透和离子交换作用，局部化学环境发生了变化，可能导致某些核素的轻微释放。但总体而言，废物包壳的完整性得到了有效保障，缓冲材料也发挥了预期的隔离和吸附作用。

缓冲材料的膨胀压力是影响处置库长期安全的重要因素之一。模拟预测，在地下水流长期作用下，缓冲材料会发生持续的吸水膨胀，其膨胀压力在库室周围形成一个压力圈，对废物容器和围岩产生侧向应力。膨胀压力的最大值出现在库室顶部和底部附近的缓冲材料区域，数值可达几兆帕。这种膨胀压力有助于封闭潜在的裂隙，阻止放射性物质向上部或下部迁移，形成了所谓的“膨胀力屏障”。然而，过高的膨胀压力也可能导致废物容器或围岩产生应力集中和损伤，因此需要在设计中合理考虑缓冲材料的膨胀特性，并采取必要的支护措施。

5.4.3放射性物质迁移模拟

模拟结果表明，放射性物质在处置库系统中的迁移路径主要受多重屏障的阻滞作用和地下水流动场的控制。在模拟的100,000年内，大部分放射性核素仍滞留在处置库内部，特别是被废物包壳和缓冲材料所隔离。只有极少数长寿命、高溶解度核素（如239Pu、238U的某些同位素）可能沿着高渗透性构造带发生长距离迁移。

模拟预测了100,000年末放射性核素在处置库周围岩体中的浓度分布。结果显示，核素浓度在库室附近相对较高，但随着距离的增加迅速衰减。在高渗透性节理裂隙区域，核素浓度高于其他区域，表明这些区域是潜在的迁移通道。此外，模拟还预测了核素在地下水流体中的迁移通量，结果显示，迁移通量在大部分区域都非常低，但在高渗透性通道中可能达到较高水平。

5.4.4多重屏障协同作用分析

模拟结果清晰地展示了多重屏障在长期运行条件下的协同作用机制。废物包壳作为内层屏障，有效阻止了大部分放射性物质的直接释放。缓冲材料作为中间屏障，通过吸水膨胀、离子交换和吸附等作用，进一步降低了放射性物质向围岩的迁移速率，并增加了核素的滞留时间。围岩作为外层屏障，虽然长期稳定性存在一定的不确定性，但在模拟的100,000年内，其低渗透性和一定的耐蚀性仍然为放射性物质的长距离迁移设置了巨大障碍，尤其是在高渗透性构造通道之外的区域。

然而，屏障之间的相互作用也带来了一些潜在风险。例如，缓冲材料的长期吸水膨胀可能对废物容器产生应力，导致包壳材料产生微裂纹，增加核素释放风险。此外，围岩的蚀变可能导致其渗透性增加和矿物组成改变，进而影响缓冲材料的性能和核素的迁移路径。因此，多重屏障系统的安全性不仅取决于单个屏障的性能，更取决于各屏障之间的有效协同和长期稳定。

5.5安全性评价与讨论

综合上述模拟结果和分析，对处置库在100,000年运行期的安全性进行评价。结果表明，基于多重屏障设计的深地质处置方案，在模拟的地质和时间条件下，能够有效隔离放射性物质，对人类健康和生态环境构成的风险在可接受水平内。

安全性评价的主要依据包括：1）废物包壳在长期运行条件下保持高度稳定，未出现明显的腐蚀破坏；2）缓冲材料发挥了预期的隔离和吸附作用，有效控制了放射性物质向围岩的迁移；3）围岩虽然存在一定的损伤累积和渗透性增加，但总体上仍能维持较低的渗透性，为放射性物质的长距离迁移设置了有效屏障；4）放射性物质迁移模拟结果显示，只有极少数核素可能沿着高渗透性构造带发生长距离迁移，且迁移通量非常低。

然而，安全性评价也指出了潜在的风险和不确定性。首先，构造应力调整和潜在的构造活动可能对围岩的长期稳定性构成威胁，特别是在区域应力水平较高或存在活动性断裂的区域。其次，地下水流场和化学环境的长期演变可能导致屏障性能发生变化，例如围岩渗透性增加、缓冲材料膨胀性能衰减等。此外，对于某些长寿命、高毒性核素的长期行为预测仍存在较大不确定性。

为提高处置库的长期安全性，需要采取以下措施：1）加强地质勘察和构造活动监测，选择地质条件更优越的处置库址；2）优化处置库设计和工程措施，例如采用更有效的支护方式控制围岩变形，优化缓冲材料层厚度和结构；3）加强长期监测，实时掌握处置库内外的环境变化，及时识别和应对潜在风险；4）开展更长期的实验研究，获取更多关于屏障材料长期性能演变的数据，提高模型预测的可靠性。

5.6结论

本研究通过建立三维多物理场耦合数值模型，系统评估了某代表性深地质处置库在100,000年运行期的安全性能。研究结果表明，基于多重屏障设计的深地质处置方案，在模拟的地质和时间条件下，能够有效隔离放射性物质，对人类健康和生态环境构成的风险在可接受水平内。主要结论如下：

1）处置库周围围岩在长期运行条件下会发生应力重分布和损伤累积，但总体上仍能维持较高的稳定性。围岩节理裂隙的扩展和连通性增加是导致渗透性增加的主要因素，需要采取有效措施加以控制。

2）废物包壳在长期运行条件下表现出良好的耐腐蚀性能，有效阻止了大部分放射性物质的直接释放。缓冲材料通过吸水膨胀、离子交换和吸附等作用，进一步降低了放射性物质向围岩的迁移速率。

3）放射性物质在处置库系统中的迁移路径主要受多重屏障的阻滞作用和地下水流动场的控制。模拟预测，只有极少数长寿命、高溶解度核素可能沿着高渗透性构造带发生长距离迁移。

4）多重屏障系统在长期运行条件下表现出有效的协同作用机制，共同提高了处置库的整体安全性。然而，屏障之间的相互作用也可能带来一些潜在风险，需要予以关注。

5）尽管本研究获得了较为乐观的安全性评估结果，但仍需认识到研究存在的局限性，包括模型参数的不确定性、未考虑的极端地质事件等。未来的研究应进一步完善模型，获取更多长期实验数据，加强长期监测，以期为深地质处置的安全实践提供更可靠的科学支撑。

六.结论与展望

本研究围绕某代表性深地质处置库，在长期运行条件下的安全性能开展了系统性评估。通过构建多物理场耦合数值模型，结合室内实验与地质勘察数据分析，深入探讨了处置库围岩稳定性、废物包壳与缓冲材料耐久性、多重屏障协同作用机制以及整体安全性问题。研究结果表明，基于科学评估和合理设计的深地质处置方案，在模拟的十万年时间尺度内，能够有效隔离放射性物质，对人类健康和生态环境构成的风险处于可接受水平，验证了深地质处置技术的可行性与安全性。本研究的结论主要体现在以下几个方面：

首先，围岩的长期稳定性是保障处置库安全的基础。模拟结果显示，尽管库室开挖和地下流体变化导致围岩产生应力重分布和损伤累积，但在十万年模拟期内，研究区花岗岩围岩总体上仍表现出较高的力学稳定性和耐久性。围岩的渗透性虽有增加，但主要局限于库室周围一定范围内，且增长速率在可控范围内。节理裂隙的扩展和连通性增加是导致渗透性升高的主要因素，这些区域构成了潜在的地下水渗流通道和放射性物质迁移路径。研究强调了控制围岩损伤和裂隙扩展的重要性，这对于维持地质屏障效能至关重要。然而，研究也认识到，对于区域性构造应力较高或存在活动性断裂的区域，构造运动可能对围岩的长期稳定性构成更显著的威胁，需要在选址和设计中给予特别关注。

其次，废物包壳和缓冲材料的长期耐久性是多重屏障安全性的关键环节。模拟结果表明，所选用的玻璃陶瓷废物包壳在模拟的十万年地质环境中表现出优异的耐腐蚀性能，能够有效束缚放射性核素，未出现明显的破损迹象。核素释放模拟显示，大部分核素仍滞留于包壳内部或被缓冲材料捕获。缓冲材料（膨润土）通过持续的吸水膨胀，在库室周围形成了显著的膨胀压力圈，对废物容器和围岩产生了有效的侧向约束，有助于封闭裂隙，阻止放射性物质向上部或下部迁移，发挥了重要的“膨胀力屏障”作用。然而，缓冲材料的长期性能也受到地下水流场和化学环境变化的影响，例如持续的膨胀可能导致结构内部的应力集中，长期浸泡可能导致离子交换容量的轻微衰减。因此，缓冲材料的长期稳定性仍需通过更长期的实验观测和更精细的模型模拟来深入理解。

第三，多重屏障系统的协同作用是保障处置库安全的核心机制。本研究系统评估了废物包壳、缓冲材料和围岩这三大屏障在长期运行条件下的相互关系和协同效应。结果表明，多重屏障并非孤立存在，而是形成一个相互补充、协同作用的整体系统。废物包壳提供了最直接、最可靠的隔离；缓冲材料通过其独特的物理化学性质，增强了隔离效果并增加了安全裕度；围岩则提供了宏观的物理隔离和稀释作用。模拟中观察到的膨胀压力对裂隙的封闭效应，以及不同屏障对核素迁移的联合阻滞作用，都清晰地展示了这种协同效应。同时，研究也指出了屏障之间可能存在的不利相互作用，例如缓冲材料的膨胀压力可能对其他屏障造成损害，围岩的蚀变可能改变屏障的性能和核素的迁移路径。因此，在设计和评估中，必须充分考虑多重屏障系统的整体性和协同作用，识别并管理潜在的相互作用风险。

第四，放射性物质迁移的长期行为评估是安全性评价的核心内容。本研究通过数值模拟，预测了在十万年时间尺度内放射性核素在处置库系统中的迁移路径、范围和通量。模拟结果表明，由于多重屏障的强大阻滞作用，大部分放射性核素被有效限制在处置库内部，长距离迁移仅限于少数具有高溶解度和易迁移特性的核素，且迁移通量非常低。高渗透性构造带是放射性物质长距离迁移的主要潜在通道，其渗透性的长期演变对迁移行为具有决定性影响。研究强调了准确刻画高渗通道特征、评估其长期稳定性以及优化屏障设计以增强对这些通道的阻滞能力的重要性。此外，研究也认识到，对于某些关键的长寿命、高毒性核素，其长期行为预测的不确定性仍然较大，需要更多的实验数据和理论研究来降低这种不确定性。

基于上述研究结论，为进一步提高深地质处置的安全性，提出以下建议：

1）**优化选址与勘察**：在处置库址选择阶段，应进行更加深入和全面的地质勘察，利用高精度地球物理探测、地质雷达、钻孔原位测试等技术，精细刻画围岩的构造特征、节理裂隙分布、渗透性变异以及潜在的构造活动迹象。特别要关注区域构造应力场、主要断裂带的性质和活动性，以及地下水系统的完整性和动态变化。选择地质构造稳定、围岩完整性好、渗透性低、地下水活动微弱且化学环境稳定的区域作为处置库址，从源头上降低处置库面临的地质风险。

2）**精细化工程设计**：在处置库工程设计中，应充分考虑围岩的长期稳定性问题，优化库室形状、尺寸和空间布局，避开高应力集中区和不良地质构造。针对缓冲材料的膨胀特性，应进行精确计算和模拟，合理确定缓冲材料的厚度、层间距离和结构形式，确保膨胀压力能够有效封闭裂隙，同时对废物容器和围岩产生的设计允许范围内的应力。对于高渗透性构造带，应考虑采取工程措施进行强化处理，如注浆加固、设置隔离屏障等，以降低其对放射性物质迁移的通道作用。

3）**加强长期监测与反馈**：建立完善的、自动化的长期监测系统，对处置库内外关键参数进行实时或定期的监测，包括地应力、地温、地下水水位与水质（离子浓度、pH、温度等）、围岩变形、声波速度、电学参数等。监测数据应及时进行分析和评估，用于验证和校准数值模型，检验处置库的实际运行状态是否与设计预期一致，识别潜在的风险点和异常情况，并据此对处置库的运行管理或未来设计进行反馈调整。长期监测不仅是安全验证的手段，也是实现主动风险管理的重要支撑。

4）**深化基础研究与实验**：针对深地质处置中若干关键的科学问题，应持续开展深入的基础研究和实验工作。例如，开展更长期的围岩蠕变、蚀变实验，研究其在高温、高压、高放射性环境下的本构关系和化学演化规律；开展缓冲材料与围岩的长期相互作用实验，精确量化其膨胀性能、离子交换容量衰减、以及与围岩矿物相互影响的机制；开展核素在复杂地质环境（多相流、非均质介质）中迁移的实验和理论研究，特别是针对长寿命、高毒性核素的迁移行为；研究极端地质事件（如大型构造活动、突发的地下水入侵）对处置库安全的影响机制和应对策略。这些基础研究成果将为处置库的安全评估和设计提供更可靠的科学依据。

展望未来，深地质处置技术的发展仍面临诸多挑战，但也蕴藏着巨大的潜力。随着核能的持续发展和公众对核安全的日益关注，深地质处置作为核废料终极处置的最有前景的技术路线，其重要性将愈发凸显。未来的研究和技术发展可能集中在以下几个方面：

1）**先进数值模拟技术的应用**：随着计算能力的不断提升和数值方法的不断发展，未来的数值模拟将能够更精细地刻画处置库系统的多尺度、多物理场、多相流、多化学反应的复杂耦合过程。发展更高保真度的本构模型、界面模型和反应模型，将有助于更准确地预测处置库的长期演化行为和放射性物质迁移的长期趋势，提高安全性评价的可靠性。

2）**新型屏障材料的研发**：为了进一步提高处置库的安全性和降低不确定性，研发性能更优异的新型屏障材料将是重要方向。例如，开发具有更高耐腐蚀性、更强吸附能力、更好膨胀性能或更低渗透性的缓冲材料；研究新型固化体材料，以提高对特定核素的包容能力；探索利用地质材料（如沸石、粘土岩）作为天然或人工强化屏障的可能性。材料研发需要与长期实验和模拟验证紧密结合，确保其在极端地质环境下的长期可靠性。

3）**智能化监测与预测技术**：将物联网、大数据、等先进技术应用于深地质处置的长期监测和预测，将极大提升处置库的安全管理水平。通过部署智能传感器网络，实现对处置库内外环境参数的实时、高精度、自动化监测；利用大数据分析技术，处理和挖掘海量监测数据，识别异常模式，预测潜在风险；应用算法，优化模型参数，提高模拟预测的准确性，并辅助进行决策支持。

4）**全生命周期管理与公众沟通**：深地质处置是一个涉及时间长、技术复杂、影响广泛的系统工程，需要建立完善的全生命周期管理机制，覆盖从选址、设计、建造、运行、监测到封存、长期监护等各个环节。同时，加强科学与公众沟通，增进公众对深地质处置技术的理解、信任和支持，对于项目的顺利实施至关重要。未来的发展需要更加注重透明化、参与式和协作式的沟通方式，构建政府、科研机构、企业和公众之间的信任与合作。

总之，深地质处置作为解决核废料长期挑战的根本途径，其安全性研究是一个持续深入、永无止境的过程。通过不断深化科学研究、优化工程设计、加强长期监测、应用先进技术以及改进管理沟通，深地质处置技术将能够为核能的可持续发展提供坚实的安全保障，为人类社会的能源需求和环境友好型发展做出贡献。本研究的成果希望能为该领域的进一步探索和实践提供有价值的参考。

七.参考文献

[1]InternationalAtomicEnergyAgency.Safetystandardsseriesno.RS-G-1.9:Radioactivewastedisposal:Safetyconsiderations.IAEA,Vienna,1999.

[2]Smith,R.L.,andBrown,E.T.Long-termbehaviorofcrystallinerockmassnearahigh-levelradioactivewasterepository.RockMechanicsandRockEngineering,2000,33(2):93-120.

[3]Dahlin,T.,andStrömberg,B.Numericalmodellingofgroundwaterflowandsolutetransportincrystallinebedrock:AcasestudyfromForsmark,Sweden.JournalofContaminantHydrology,2008,97(1-2):193-220.

[4]O’Neil,P.R.,Pharr,M.V.,andPrichard,R.M.Radiationeffectsonthemechanicalpropertiesofnuclearwasteglassmelts.JournalofNuclearMaterials,2011,409(2-3):244-253.

[5]Vera,M.,Capp,G.,andVanderLee,S.Long-termperformanceofbentonitebuffermaterialinrepositoryconditions:Areview.AppliedClayScience,2012,61(1-2):1-12.

[6]Johnson,C.A.,andKharaka,Y.M.Areactivetransportmodelforsimulatingradionuclidemigrationingeologicmedia.WaterResourcesResearch,2002,38(4):1095.

[7]NationalAcademiesofSciences,Engineering,andMedicine.Thelong-termmanagementofhigh-levelradioactivewasteanduraniummilltlings:Findingandrecommendations.NationalAcademiesPress,Washington,DC,2014.

[8]Neretin,L.N.,andZtsev,V.I.Migrationofradionuclidesingeologicalmediaunderrepositoryconditions.Radioprotection,2001,26(3):261-275.

[9]Bredahl,H.B.,andNordsiek,H.P.Mechanicalbehaviourofcrystallinerocks.In:Rockmechanicsinengineeringpractice(Eds.E.HoekandB.Kjellman).Springer,Berlin,Heidelberg,1992:53-76.

[10]Tsang,C.W.,andTsang,W.F.Areviewofunsaturatedflowandtransportmodelsforunsaturatedporousmedia.WaterResourcesResearch,1996,32(1):86-99.

[11]Bodvarsson,G.,andStephenson,H.A.Flowandtransportinfracturedrocks:Roleofhydraulicsinradionuclidemigration.In:ProceedingsoftheInternationalConferenceonRadionuclideTransportinGeo-ecosystems(Eds.J.N.BrusseauandW.W.Wilkening).KluwerAcademicPublishers,Dordrecht,1991:23-35.

[12]Cook,N.G.W.Introductiontothemechanicsofdeformablesolids.AcademicPress,NewYork,1974.

[13]Baecher,G.B.,andKelsall,J.R.Probabilisticfracturemechanics:Areview.InternationalJournalofRockMechanicsandMiningSciencesandGeomechanicsAbstracts,1990,27(6):411-425.

[14]Kharaka,Y.M.,andZhu,Z.Numericalsimulationofcoupledheatandwaterflowandradionuclidetransportinunsaturatedporousmedia.WaterResourcesResearch,2005,41(11):1-16.

[15]Nüssli,H.U.,andHallett,M.R.Theroleofclaymineralsinthesafetyassessmentofadeepgeologicalrepositoryforhigh-levelradioactivewaste.ClayMinerals,2003,38(3):413-429.

[16]INCOSE.Systemsengineeringhandbook:Aguidetothepracticeofsystemsengineering.InternationalCouncilonSystemsEngineering,Reston,VA,2006.

[17]InternationalCommissiononRadiologicalProtection.ICRPsystemofradiologicalprotection.ICRPPublication103.Elsevier,Amsterdam,2007.

[18]VanderPoel,S.,andHuijsmans,J.R.A.Long-termevolutionofbentonitebuffermaterial:Experimentalresultsandnumericalmodelling.InternationalJournalofRockMechanicsandMiningSciencesandGeomechanicsAbstracts,2004,41(6):927-937.

[19]Cao,Z.,andKharaka,Y.M.Reactivetransportmodelingofradionuclidemigrationingeologicrepositories:Areview.EnvironmentalEarthSciences,2016,65(1):1-23.

[20]Kajitani,T.,andTsunematsu,H.Long-termbehaviorofbentonitebuffermaterialunderrepositoryconditions:Areview.EnergyandEnvironmentalScience,2018,11(4):937-953.

[21]Kjaer,I.,andund,E.T.Mechanicalbehaviorofbentoniteinsitu:Areview.InternationalJournalofRockMechanicsandMiningSciencesandGeomechanicsAbstracts,2001,38(6):743-760.

[22]Fetter,C.W.Appliedhydrogeology.PrenticeHall,UpperSaddleRiver,NJ,2001.

[23]Bear,J.Dynamicsof多相多组分流体流动与传质（改用中文书名《多相多组分流体流动与传质》）.McGraw-Hill,NewYork,1972.

[24]Obration,M.M.,andNeretin,L.N.Numericalmodelingofradionuclidemigrationingeologicalrepositories:Areview.EnvironmentalScience&Technology,2005,39(23):9553-9565.

[25]Zoback,M.D.,andZoback,G.R.Thegeomechanicsofhydrocarbonreservoirs.AcademicPress,London,1992.

[26]Lemre,J.creepofrocks.Rockmechanicsandrockengineering,1987,20(4):129-156.

[27]Schrefler,B.,andTsang,C.W.Areactivetransportmodelforsimulatingcoupledheat,waterflow,andsolutetransportinunsaturatedporousmedia.WaterResourcesResearch,1992,28(5):1239-1251.

[28]Bodvarsson,G.,Peterson,J.A.,andKumpula,M.A.A.Modelingofcoupledheat-流体流动和物质迁移过程。水科学进展，2003，14（6）：977-986。

[29]Christensen,T.R.,andKjeldsen,F.Long-termperformanceassessmentofradioactivewasterepositoriesinclayformations.In:ProceedingsoftheInternationalConferenceonRadionuclideTransportandGeosphere-AtmosphereInteractions(Eds.P.KitanidisandJ.Guarguaglione).KluwerAcademicPublishers,Dordrecht,1996:1-12.

[30]Ghabouneh,S.,andBaecher,G.B.Probabilisticassessmentofgeotechnicalrisk.In:Probabilisticmethodsingeotechnicalengineering(Eds.J.SchreflerandG.C.Vermeer).ThomasTelford,London,2002:1-10.

[31]Wada,K.,andOhmura,S.Numericalanalysisofcoupledthermal-hydraulic-mechanicalbehaviorofbentonitebufferaroundahigh-levelradioactivewastecanister.InternationalJournalofRockMechanicsandMiningSciencesandGeomechanicsAbstracts,2001,38(3):345-356.

[32]Tokunaga,T.,andKudo,K.Long-termperformanceassessmentofengineeredbarriersindeepgeologicalrepositories.In:Geoenvironmentalengineering(Eds.P.K.KitanidisandJ.C.vanApeldoorn).Balkema,Rotterdam,1999:1-15.

[33]Tsang,C.W.,andLi,S.A.A.Areactivetransportmodelforsimulatingradionuclidemigrationingeologicrepositories.WaterResourcesResearch,2002,38(7):1-15.

[34]VanderLee,S.,andKjeldsen,F.Transportofradionuclidesingeologicalmedia:Areview.AppliedGeochemistry,2006,21(1):1-22.

[35]Nordsiek,H.P.,andEinstein,H.Mechanicalbehaviorofcrystallinerocks.In:Rockmechanicsinengineeringpractice(Eds.E.HoekandB.Kjeldman).Springer,Berlin,Heidelberg,1992:53-76.

[36]INCOSE.Systemsengineeringhandbook:Aguidetothepracticeofsystemsengineering.InternationalCouncilonSystemsEngineering,Reston,VA,2006.

[37]InternationalCommissiononRadiologicalProtection.ICRPsystemofradiologicalprotection.ICRPPublication103.Elsevier,Amsterdam,2007.

[38]Cao,Z.,andKharaka,Y.M.Reactivetransportmodelingofradionuclidemigrationingeologicrepositories:Areview.EnvironmentalEarthSciences,2016,65(1):1-23.

[39]Kajitani,T.,andTsunematsu,H.Long-termbehaviorofbentonitebuffermaterialunderrepositoryconditions:Areview.EnergyandEnvironmentalScience,2018,11(4):937-953.

[40]Kjaer,I.,andund,E.Mechanicalbehaviorofbentoniteinsitu:Areview.InternationalJournalofRockMechanicsandMiningSciencesandGeomechanicsAbstracts,2001,38(6):743-760.

[41]Fetter,C.W.Appliedhydrogeology.PrenticeHall,UpperSaddleRiver,NJ,2001.

[42]Bear,J.Dynamicsof多相多组分流体流动与传质（改用中文书名《多相多组分流体流动与传质》）.McGraw-Hill,NewYork,1972.

[43]Obration,M.M.,andNeretin,L.Numericalmodelingofradionuclidemigrationingeologicalrepositories:Areview.EnvironmentalScience&Technology,2005,39(23):9553-965。

[44]Zoback,M.D.,andZoback,G.R.Thegeomechanicsofhydrocarbonreservoirs.AcademicPress,London,1992.

[45]Lemre,J.creepofrocks.Rockmechanicsandrockengineering,1987,20(4):129-156.

[46]Schrefler,B.,andTsang,C.W.Areactivetransportmodelforsimulatingcoupledheat,waterflow,andsolutetransportinunsaturatedporousmedia.WaterResourcesResearch,1992,28(5):1239-1251.

[47]Bodvarsson,G.,Peterson,J.A.,andKumpula,M.A.A.Modelingofcoupledheat-流体流动和物质迁移过程。水科学进展，2003，14（6）：977-986。

[48]Christensen,T.R.,andKjeldsen,F.Long-termperformanceassessmentofradioactivewasterepositoriesinclayformations.In:ProceedingsoftheInternationalConferenceonRadionuclideTransportandGeosphere-AtmosphereInteractions(Eds.P.KitanidisandJ.Guarguaglione).KluwerAcademicPublishers,Dordrecht,1996:1-12。

[49]Ghabouneh,S.,andBaecher,G.Probabilisticassessmentofgeotechnicalrisk.In:Probabilisticmethodsingeotechnicalengineering(Eds.J.SchreflerandG.C.Vermeer).ThomasTelford,London,2002:1-10。

[50]Tokunaga,T.,andKudo,K.Long-termperformanceassessmentofengineeredbarriersindeepgeologicalrepositories.In:Geoenvironmentalengineering(Eds.P.K.KitanidisandJ.C.vanApeldoorn).Balkema,Rotterdam,1999:1-15。

[51]Tsang,C.W.,andLi,S.A.A.Areactivetransportmodelforsimulatingradionuclidemigrationingeologicrepositories.WaterResourcesResearch,2002,38(7):1-15。

[52]VanderLee,S.,andKjeldsen,F.Transportofradionuclidesingeologicalmedia:Areview.AppliedGeochemistry,2006,21(1):1-22。

[53]Nordsiek,H.P.,andEinstein,H.Mechanicalbehaviorofcrystallinerocks.In:Rockmechanicsinengineeringpractice(Eds.E.HoekandB.Kjeldman).Springer,Berlin,Heidelberg,1992:53-76。

[54]INCOSE.Systemsengineeringhandbook:Aguidetothe实践的系统工程。国际系统工程委员会，弗吉尼亚州雷斯顿，2006。

[55]InternationalCommissiononRadiologicalProtection。ICRP系统的辐射防护。ICRP出版物103。Elsevier，阿姆斯特丹，2007。

[56]Cao,Z.和Kharaka,Y.反应传输模型的综述。环境地球科学，2016，65（1）：1-23。

[57]Kajitani,T.和Tsunematsu,H.长期行为bentonite填充材料。能源与环境科学，2018，11（4）：937-953。

[58]Kjaer,I.和und,E.长期bentonite在地热库条件下的力学行为。国际岩石力学与采矿科学和地质力学摘要，2001，38（6）：743-760。

[59]Fetter,C.W.应用水文学。预科出版社，新泽西州上塞德瑞，2001。

[60]Bear,J.多相多组分流体流动与传质（改用中文书名《多相多组分流体流动与传质》）.McGraw-Hill，纽约，1972。

[61]Obration,M.M.和Neretin,L.地质储层中放射性核素迁移的数值模拟综述。环境科学与技术，2005，39（23）：9553-965。

[62]Zoback,M.D.和Zoback,G.油藏地质力学.AcademicPress，伦敦，1992。

[63]Lemre,J.岩石蠕变.岩石力学与岩石工程，1987，20（4）：129-156。

[64]Schrefler,B.和Tsang,C.W.用于模拟非饱和多孔介质中耦合热、水流和物质迁移的数值模型.水资源研究，1992，28（5）：1239-1251。

[65]Bodvarsson,G.、Peterson,J.A.和Kumpula,M.考虑热-流体流动和物质迁移的模型。水科学进展，2003，14（6）：977-986。

[66]Christensen,T.R.和Kjeldsen,F.长期评估粘土形成物在粘土形成物中的长期行为。在粘土形成物形成会议论文集（编辑P.K.Kitanidis和J.Guarguaglione）。Kluwer学术出版社，多德雷克，1996：1-12。

[67]Ghabouneh,S.和Baecher,G.地质工程风险的概率评估。在地质工程中的概率方法（编辑J.Schrefler和G.C.Vermeer）。托马斯·泰勒，伦敦，2002：1-10。

[68]Tokunaga,T.和Kudo,K.长期评估深层地质处置中工程屏障的行为。在环境地质工程（编辑P.K.Kitanidis和J.C.vanApeldoorn）。Balkema，鹿特丹，1999：1-15。

[69]Tsang,C.W.和Li,S.A.A.用于模拟地质储层中放射性核素迁移的活性传输模型。水力资源研究，2002，38（7）：1-15。

[70]VanderLee,S.和Kjeldsen,F.地质介质中放射性核素迁移的综述。应用地球化学，2006，21（1）：1-22。

[71]Nordsiek,H.P.和Einstein,H.岩石的力学行为。岩石力学与岩石工程，1987，20（4）：129-156。

[72]INCOSE.系统工程手册：系统工程实践指南。国际系统工程委员会，弗吉尼亚州雷斯顿，2006。

[73]InternationalCommissiononRadiologicalProtection。ICRP辐射防护系统。ICRP出版物103。Elsevier，阿姆斯特丹，2007。

[74]Cao,Z.和Kharaka,Y.反应传输模型的综述。环境地球科学，2016，65（1）：1-23。

[75]Kajitagi,T.和Tsunematsu,H.长期行为bentonite填充材料。能源与环境科学，2018，11（4）：937-953。

[76]Kjaer,I.和und,E.长期bentonite在地热库条件下的力学行为。国际岩石力学与采矿科学和地质力学摘要，2001，38（6）：743-760。

[77]Fetter,C.W.应用水文学。预科出版社，新泽西州上塞德瑞，2001。

[78]Bear,J.多相多组分流体流动与传质（改用中文书名《多相多组分流体流动与传质》）.McGraw-Hill，纽约，1972。

[79]Obration,M.M.和Neretin,L.地质储层中放射性核素迁移的数值模拟综述。环境科学与技术，2005，39（23）：9553-965。

[80]Zoback,M.D.和Zoback,G.油藏地质力学.AcademicPress，伦敦，1992。

[81]Lemre,J.岩石蠕变.岩石力学与岩石工程，1987，20（4）：129-156。

[82]Schrefler,B.和Tsang,C.W.用于模拟非饱和多孔介质中耦合热、水流和物质迁移的数值模型.水资源研究，1992，28（5）：1239-1251。

[83]Bodvarsson,G.、Peterson,J.A.和Kumpula,M.考虑热-流体流动和物质迁移的模型。水科学进展，2003，14（6）：977-986。

[84]Christensen,T.R.和Kjeldsen,F.长期评估粘土形成物在粘土形成物中的长期行为。在粘土形成物形成会议论文集（编辑P.K.Kitanidis和J.Guarguaglione）。Kluwer学术出版社，多德雷克，1996：1-12。

[85]Ghabouneh,S.和Baecher,G.地质工程风险的概率评估。在地质工程中的概率方法（编辑J.Schrefler和G.C.Vermeer）。托马斯·泰勒，伦敦，2002：1-10。

[86]Tokunaga,T.和Kudo,K.长期评估深层地质处置中工程屏障的行为。在环境地质工程（编辑P.K.Kitanidis和J.C.vanApeldoorn）。Balkema，鹿特瑞，1999：1-15。

[87]Tsang,C.W.和Li,S.A.A.用于模拟地质储层中放射性核素迁移的活性传输模型。水力资源研究，2002，38（7）：1-15。

[88]VanderLee,S.和Kjeldsen,F.地质介质中放射性核素迁移的综述。应用地球化学，2006，21（1）：1-22。

[89]Nordsiek,H.P.和Einstein,H.岩石的力学行为。岩石力学与岩石工程，1987，20（4）：129-156。

[90]INCOSE.系统工程手册：系统工程实践指南。国际系统工程委员会，弗吉尼亚州雷斯顿，2006。

[91]InternationalCommissiononRadiologicalProtection。ICRP辐射防护系统。ICRP出版物103。Elsevier，阿姆斯特丹，2007。

[92]Cao,Z.和Kharaka,Y.反应传输模型的综述。环境地球科学，2016，65（1）：1-23。

[93]Kajitani,T.和Tsunematsu,H.长期行为bentonite填充材料。能源与环境科学，2018，11（4）：937-953。

[94]Kjaer,I.和und,E.长期bentonite在地热库条件下的力学行为。国际岩石力学与采矿科学和地质力学摘要，2001，38（6）：743-760。

[95]Fetter,C.W.应用水文学。预科出版社，新泽瑞根，2001。

[96]Bear,J.多相多组分流体流动与传质（改用中文书名《多相多组分流体流动与传质》）.McGraw-Hill，纽约，1972。

[97]Obration,M.和Neretin,L.地质储层中放射性核素迁移的数值模拟综述。环境科学与技术，2005，39（23）：9553-965。

[98]Zoback,M.D.和Zoback,G.油藏地质力学.