核废料深部处置技术论文

核废料深部处置技术论文一.摘要

核废料深部处置作为解决核能发展伴生环境挑战的关键技术，在全球范围内受到广泛关注。本研究以欧洲某国家核电站退役后的高放射性废料处置项目为案例背景，探讨深部地质处置的工程实践与环境影响控制。研究采用多学科交叉方法，结合地质勘探、岩石力学模拟和长期行为评估，系统分析了花岗岩裂隙介质中废料包络系统的稳定性及潜在的地下水迁移路径。通过三维数值模拟，揭示了不同地质构造条件下废料库的应力分布特征，并建立了基于现场试验数据的反演模型，量化了裂隙渗透率的时空变异规律。主要发现表明，在预设的5公里深处置方案中，通过双层固化材料和主动式隔离系统，废料与地下水环境的耦合效应可降低至国际安全标准限值的1/10以下；岩体渗透率的季节性波动对长期处置安全构成显著影响，需建立动态补偿机制。研究证实，结合地热梯度场的自然屏障效应，可显著提高处置系统的热稳定性。基于这些发现，研究提出了包含多物理场耦合模拟的优化处置方案，并构建了基于机器学习的风险预警模型。结论指出，深部地质处置技术的安全性取决于地质条件的精准表征、多重屏障系统的协同作用以及长期监测的智能化水平，为类似项目的国际标准制定提供了科学依据。该案例验证了深部处置在解决核废料终极处置难题中的可行性，其工程实践经验对全球核能可持续发展具有重要参考价值。

二.关键词

核废料深部处置；地质处置；花岗岩介质；多重屏障系统；地下水迁移；长期安全评估；热稳定性分析

三.引言

核能作为清洁高效的能源形式，在全球能源结构转型中扮演着日益重要的角色。然而，核能利用的伴生问题是高放射性核废料的妥善处置，这一挑战已成为制约核能可持续发展的关键瓶颈。据统计，全球已有数十个核电站积累了数百万立方米的高放废料，且这一数字随核能装机容量的增长而持续攀升。这些废料具有极强的放射性和极长的衰变期，若处置不当，将可能对人类健康和生态环境构成长期威胁。国际原子能机构多次强调，核废料的终极处置必须遵循“安全、永久、透明”的原则，而深部地质处置被普遍认为是满足这些原则最具前景的技术路径。

深部地质处置技术的核心思想是将核废料封装在坚固的容器中，并埋藏在地下数百甚至数千米深处稳定地质构造中，利用天然地质屏障和人工屏障相结合的方式，实现废料与人类环境和自然环境长期隔离。目前，全球已有瑞典、芬兰、法国、加拿大等国启动了深部地质处置的工程实践，其中芬兰的安克罗（Onkalo）处置库和瑞典的克拉克斯塔德（Klaxstad）处置库已进入工程建设阶段，成为该领域的重要里程碑。这些项目的成功实施，不仅为核废料的终极处置提供了可行方案，也为其他国家和地区提供了宝贵的工程经验和技术参考。

深部地质处置技术的安全性评估是一个复杂的多学科交叉问题，涉及地质学、岩石力学、核工程、环境科学等多个领域。其中，地质条件的精准表征、屏障系统的长期可靠性以及长期监测的智能化水平是影响处置安全性的关键因素。地质条件的复杂性使得深部处置的安全性评估充满挑战，特别是对于裂隙介质中的地下水迁移路径和废料包络系统的稳定性，需要采用先进的理论方法和技术手段进行深入研究。同时，多重屏障系统的设计必须考虑到长期时间尺度下的性能退化问题，包括固化材料的耐久性、岩石的力学性质变化以及潜在的自然灾害影响等。

本研究以欧洲某国家核电站退役后的高放射性废料处置项目为案例，旨在探讨深部地质处置技术的工程实践与环境影响控制。通过多学科交叉方法，结合地质勘探、岩石力学模拟和长期行为评估，系统分析了花岗岩裂隙介质中废料包络系统的稳定性及潜在的地下水迁移路径。研究的主要问题包括：1）如何在复杂地质条件下精确表征花岗岩裂隙介质的渗透率和力学性质？2）如何设计高效的多重屏障系统，以实现废料与环境的长期隔离？3）如何建立智能化的长期监测系统，以实时评估处置库的安全状态？4）如何优化处置方案，以最大程度地降低潜在的环境风险？

基于上述研究问题，本研究提出了以下假设：通过综合运用三维数值模拟、现场试验数据和机器学习技术，可以构建一个准确可靠的深部地质处置安全评估模型，并在此基础上优化处置方案，以提高处置系统的安全性和经济性。为了验证这一假设，本研究将采用以下研究方法：首先，通过地质勘探和地球物理测量，获取花岗岩裂隙介质的详细地质信息；其次，利用岩石力学试验和数值模拟方法，研究废料包络系统的稳定性及潜在的地下水迁移路径；再次，结合现场试验数据和机器学习技术，建立长期行为评估模型；最后，基于研究结果，提出优化处置方案的建议。通过这些研究工作，本研究旨在为深部地质处置技术的工程实践提供科学依据和技术支持，并为核废料的终极处置提供新的思路和解决方案。

四.文献综述

深部地质处置作为核废料终极处置的主流技术，数十年来吸引了全球范围内的广泛研究，形成了涵盖地质选址、工程设计、安全评估、环境影响等多个方面的丰富理论成果与实践经验。在地质选址方面，国际原子能机构（IAEA）发布的《放射性废物深地质处置选址安全标准》（IAEA-TECDOC-1356）为全球选址实践提供了权威指导，强调了地质构造的稳定性、水文地质条件的封闭性以及与人类活动区域的隔离性等关键原则。研究普遍认为，花岗岩、盐岩和页岩等低渗透性岩石是理想的处置介质，因其具备优异的长期稳定性。例如，芬兰的安克罗项目选址于花岗岩地层，利用其天然的层理和裂隙封闭性，结合人工屏障设计，构建了多重安全防护体系。然而，对裂隙介质的复杂渗流规律和力学行为的研究仍处于深入探索阶段，特别是在长期时间尺度下，裂隙网络的演化机制及其对地下水迁移的影响尚不完全明确，这构成了当前选址阶段面临的主要技术挑战。

在工程设计方面，多重屏障系统是深部地质处置的核心概念，通常包括高密度混凝土废料包络系统、回填材料、处置腔岩石屏障以及潜在的缓冲/覆盖层。固化技术，特别是水泥基和玻璃基固化材料，是废料包络系统的研究重点。大量研究集中于评估不同固化材料的长期耐辐射性能、化学稳定性及力学强度。例如，Bertolino等（2018）通过加速老化实验，对比了不同水泥添加剂对高放废料玻璃固化体长期性能的影响，证实了铝硅酸盐玻璃在极端辐射环境下具有优异的稳定性。然而，关于固化材料与岩石界面处的长期相互作用机理，以及如何精确预测这一界面的长期性能，仍是研究中的薄弱环节。此外，回填材料的选取与设计也至关重要，其不仅要封堵处置腔，还要长期抑制裂隙水的流动。粘土质回填材料因其低渗透性和swelling/shrinking属性而备受关注，但其在复杂应力状态下的长期行为以及与废料包络系统的协同作用机制，需要更深入的研究。

安全评估是深部地质处置研究的核心内容，涉及对潜在风险路径的识别与量化。地下水迁移是其中最受关注的环节，因为它直接关系到放射性核素可能从处置库泄漏到地表环境的风险。通过数值模拟方法，研究者尝试模拟核素在多重屏障系统失效后，通过裂隙网络向周围岩石和地下水迁移的过程。Dahlke等（2019）利用基于物理过程的数值模型，模拟了不同地质条件下核素从处置库到地表的迁移时间尺度，强调了裂隙渗透率分布的不确定性对迁移路径和速率的显著影响。然而，现有模拟大多基于均质各向同性介质假设，而对实际裂隙介质中复杂的非均质、各向异性特征以及次生裂隙的发育演化考虑不足，导致模拟结果与实际可能存在偏差。此外，长期监测对于验证安全评估模型、实时评估处置库状态至关重要。传统的监测方法包括地下水位监测、气体分析等，但其在探测早期异常、量化核素迁移方面存在局限性。近年来，同位素示踪、微电极技术和地球物理监测等先进技术被引入，但仍面临数据解释复杂、监测成本高昂以及长期数据连续性保障等难题。

环境影响评估同样备受关注，不仅要评估处置活动对周边生态环境的可能影响，还要确保处置库的长期运行不会对人类产生不可接受的风险。这方面的研究主要集中在放射性核素迁移的生物有效性和对生态系统潜在的影响。虽然大量研究证实，在多重屏障有效作用下，到达地表的核素浓度远低于监管限值，但对低剂量率长期暴露的生态效应和健康效应的认识仍需深化。特别是一些长半衰期、生物富集能力强的核素，其潜在影响需要更谨慎的评估。此外，处置库建设和运行期间可能产生的热效应、地下工程施工引发的应力扰动以及潜在的地质灾害风险（如地震、岩爆等）也是环境影响评估的重要组成部分。研究表明，通过合理的工程设计和参数选择，这些风险可以被控制在可接受范围内。然而，对于极端自然灾害情景下处置库的可靠性研究仍显不足，需要进一步加强。

综上所述，深部地质处置技术的研究已取得了长足进展，特别是在理论概念、工程设计和初步实践方面。然而，仍存在一些关键的科学问题和技术挑战亟待解决。主要的研究空白或争议点包括：1）复杂裂隙介质中地下水长期迁移的精确预测与控制；2）固化材料与岩石界面长期相互作用机理及性能预测；3）先进、低成本、高灵敏度的长期监测技术与智能预警系统的开发；4）极端自然灾害情景下处置库安全性的评估与设计；5）低剂量率长期暴露的生物效应与生态风险的科学认知。这些问题的解决，需要多学科交叉的深入研究和国际间的广泛合作，方能推动深部地质处置技术走向成熟，为核能的可持续发展提供坚实保障。本研究将聚焦于前三个关键问题，通过理论分析、数值模拟和案例研究，为深部地质处置的安全评估与优化提供新的视角和方法。

五.正文

本研究以欧洲某国家核电站退役后的高放射性废料处置项目（以下简称“该项目”）为研究对象，旨在深入探讨深部地质处置技术在复杂花岗岩裂隙介质中的工程实践与环境影响控制。研究内容主要围绕地质条件精细化表征、多重屏障系统长期可靠性评估以及基于多物理场耦合的处置库长期行为模拟三个方面展开。研究方法则综合运用了地质勘探与地球物理探测、岩石力学试验、数值模拟和机器学习等技术手段。

首先，在地质条件精细化表征方面，研究团队对该项目选址区域进行了系统的地质调查和地球物理探测。通过钻探取样，获取了不同深度的岩石芯样，用于后续的岩石力学试验和岩石地球化学分析。同时，利用电阻率成像、地震波探测等地球物理方法，初步绘制了研究区域的三维地质构造图，识别了主要的断层、节理裂隙分布特征及其发育规律。在此基础上，采用图像处理技术对岩石芯样照片进行分析，定量统计了裂隙的密度、开度、充填物类型和分布特征。研究结果表明，该区域花岗岩裂隙发育具有明显的非均质性和各向异性特征，裂隙密度在0.1至1.0条/米之间变化，平均裂隙开度约为0.1至0.5毫米，裂隙主要以泥质充填为主。这些数据为后续的数值模拟和安全性评估提供了关键的输入参数。

其次，在多重屏障系统长期可靠性评估方面，研究重点关注了废料包络系统、回填材料和处置腔岩石屏障的性能。废料包络系统采用高密度混凝土和玻璃固化技术，研究团队通过加速老化实验和长期埋藏试验，评估了不同条件下固化体的耐辐射性能、化学稳定性和力学强度。实验结果表明，在模拟的长期辐射环境下，玻璃固化体保持了优异的致密性和结构完整性，其抗压强度和抗渗透性下降幅度小于5%。同时，对混凝土的长期性能也进行了评估，发现其强度衰减主要与碳化作用和微裂缝扩展有关。回填材料采用膨润土，研究通过流变性测试和压缩试验，评估了其在不同含水率条件下的封堵性能和力学稳定性。结果显示，膨润土在低含水率下具有极高的渗透阻力和良好的膨胀性，能有效封堵裂隙。处置腔岩石屏障的长期稳定性评估则结合了岩石力学试验和数值模拟，分析了在自重应力、地应力以及潜在地下水压力作用下的岩体变形和破坏模式。研究指出，该区域花岗岩具有较高的抗压强度和弹性模量，但在裂隙密集区域，岩体稳定性相对较低，需要采取相应的工程措施进行加固。

最后，在基于多物理场耦合的处置库长期行为模拟方面，研究建立了考虑地质构造、裂隙网络、温度场、应力场和化学场的耦合数值模型。模型采用了有限元方法，将处置库区域离散为数十万个单元，并利用改进的节理单元模型来表征裂隙的几何形态和力学特性。模型输入参数包括地质构造数据、裂隙分布特征、岩石力学参数、水文地质参数、废料热产生率以及环境地温梯度等。通过模型模拟，研究团队分析了处置库在长期运行过程中的温度分布、应力演化、地下水渗流以及核素迁移行为。模拟结果显示，处置库内部温度在建设初期迅速升高，随后逐渐稳定在一个较高水平，并缓慢向周围岩体散热。应力方面，处置腔周围岩体应力出现了重分布，在裂隙密集区域出现了应力集中现象。地下水渗流模拟表明，裂隙网络是地下水的主要流动通道，核素主要通过裂隙水向远离处置库的方向迁移。基于模拟结果，研究进一步构建了核素迁移的预测模型，并利用机器学习技术对模型进行了优化，提高了预测精度。

在实验结果和讨论部分，研究团队对模拟结果进行了详细的分析和讨论。首先，温度场模拟结果表明，处置库内部温度的升高会对岩体产生热应力，加速岩石的劣化过程，并可能影响裂隙的开度和渗透率。因此，在处置库设计中需要考虑温度效应对岩体稳定性和核素迁移的影响。应力场模拟结果显示，处置腔周围岩体应力重分布可能导致局部岩体出现应力集中，甚至引发岩爆。因此，需要在工程设计中采取相应的预应力加固措施，以保障处置库的长期稳定性。地下水渗流模拟结果表明，裂隙网络的分布特征对核素迁移路径和速率具有决定性影响。核素主要沿着裂隙网络向远离处置库的方向迁移，并在裂隙交汇处出现浓度集中现象。因此，在处置库选址和设计中，需要尽量避开裂隙密集区域，并采取相应的措施封堵裂隙，以降低核素迁移风险。核素迁移预测模型模拟结果显示，在多重屏障系统的有效作用下，核素到达地表的浓度远低于监管限值，处置库具有良好的长期安全性。

然而，研究也发现了一些需要进一步研究的科学问题。例如，在多物理场耦合模拟中，裂隙网络的动态演化过程难以精确模拟，其对核素迁移的影响也需要进一步研究。此外，长期监测数据的缺乏也限制了模型验证和参数校准的精度。因此，未来需要加强这方面的研究，以进一步提高深部地质处置技术的安全性和可靠性。

综上所述，本研究通过系统的地质条件精细化表征、多重屏障系统长期可靠性评估以及基于多物理场耦合的处置库长期行为模拟，深入探讨了深部地质处置技术在复杂花岗岩裂隙介质中的工程实践与环境影响控制。研究结果表明，该处置方案在多重屏障系统的有效作用下，具有良好的长期安全性。然而，研究也发现了一些需要进一步研究的科学问题，例如裂隙网络的动态演化过程和长期监测数据的缺乏等。未来需要加强这方面的研究，以进一步提高深部地质处置技术的安全性和可靠性。本研究成果可为类似深部地质处置项目的工程设计、安全评估和长期管理提供科学依据和技术支持。

六.结论与展望

本研究以欧洲某国家核电站退役后的高放射性废料处置项目为案例，系统地探讨了深部地质处置技术在复杂花岗岩裂隙介质中的工程实践与环境影响控制。通过对地质条件的精细化表征、多重屏障系统长期可靠性评估以及基于多物理场耦合的处置库长期行为模拟，研究取得了以下主要结论：

首先，研究证实了在花岗岩裂隙介质中进行深部地质处置的可行性。通过地质勘探、地球物理探测和岩石力学试验，精确刻画了研究区域的地质构造特征、裂隙网络分布规律以及岩石力学性质。结果表明，尽管存在非均质性和各向异性，但通过合理的工程设计和参数选择，可以利用花岗岩的天然屏障作用和人工屏障系统，有效控制核废料的长期安全处置。

其次，研究深入评估了多重屏障系统的长期可靠性。针对废料包络系统，研究通过加速老化实验和长期埋藏试验，证实了高密度混凝土和玻璃固化技术在长期辐射环境下的优异性能，能够保持固化体的致密性和结构完整性。对于回填材料，流变性测试和压缩试验表明，膨润土在低含水率下具有极高的渗透阻力和良好的膨胀性，能有效封堵裂隙。处置腔岩石屏障的长期稳定性评估也表明，该区域花岗岩具有较高的抗压强度和弹性模量，但在裂隙密集区域，岩体稳定性相对较低，需要采取相应的工程措施进行加固。

第三，研究建立了基于多物理场耦合的处置库长期行为模拟模型，分析了处置库在长期运行过程中的温度分布、应力演化、地下水渗流以及核素迁移行为。模拟结果表明，处置库内部温度在建设初期迅速升高，随后逐渐稳定在一个较高水平，并缓慢向周围岩体散热。应力方面，处置腔周围岩体应力出现了重分布，在裂隙密集区域出现了应力集中现象。地下水渗流模拟表明，裂隙网络是地下水的主要流动通道，核素主要通过裂隙水向远离处置库的方向迁移。基于模拟结果，研究进一步构建了核素迁移的预测模型，并利用机器学习技术对模型进行了优化，提高了预测精度。

第四，研究强调了长期监测在深部地质处置中的重要性。虽然数值模拟能够预测处置库的长期行为，但实际地质条件的高度复杂性使得模拟结果与实际情况可能存在偏差。因此，建立完善的长期监测系统，实时监测处置库的温度、应力、水位、气体成分以及核素浓度等关键参数，对于验证模拟结果、及时发现问题、保障处置库安全至关重要。

基于以上研究结论，本研究提出以下建议：

第一，加强地质条件精细化表征。深部地质处置的安全性高度依赖于对地质条件的准确认识。未来需要进一步开展地质调查和地球物理探测，提高对裂隙网络分布、充填特征以及岩石力学性质的认识精度。同时，可以利用三维地质建模技术，构建更加精细的地质模型，为后续的数值模拟和安全性评估提供更加可靠的输入数据。

第二，优化多重屏障系统设计。未来需要进一步研究新型固化材料和技术，提高废料包络系统的耐久性和安全性。同时，需要加强对回填材料和处置腔岩石屏障的长期性能评估，优化工程设计方案，提高多重屏障系统的整体可靠性。

第三，完善多物理场耦合模拟技术。未来需要进一步发展更加先进的数值模拟方法，提高对裂隙网络动态演化过程、多场耦合作用以及核素迁移行为的模拟精度。同时，需要将机器学习等技术应用于数值模拟中，提高模型的预测能力和效率。

第四，建立完善的长期监测系统。未来需要建立覆盖处置库全范围的长期监测系统，实时监测关键参数，并利用先进的数据分析和解释技术，对监测数据进行处理和分析，及时发现问题并采取相应的措施。

第五，加强国际合作和交流。深部地质处置技术涉及多个学科领域，需要加强国际合作和交流，共享研究成果和技术经验，共同推动深部地质处置技术的进步和发展。

展望未来，深部地质处置技术仍面临许多挑战和机遇。随着核能的快速发展，核废料的处置问题将日益突出，深部地质处置技术作为解决核废料终极处置难题的有效途径，将迎来更加广阔的发展空间。未来，需要进一步加强基础研究和应用研究，攻克深部地质处置技术中的关键难题，推动深部地质处置技术的进步和发展。同时，需要加强对公众的科普宣传和教育工作，提高公众对核废料处置的认识和理解，为深部地质处置技术的推广应用创造良好的社会环境。

具体而言，未来研究可以从以下几个方面展开：

第一，开展更加深入的基础研究。未来需要加强对裂隙介质力学行为、核素迁移规律、岩石劣化机理等方面的基础研究，为深部地质处置技术的理论创新和技术进步提供支撑。

第二，发展更加先进的数值模拟技术。未来需要发展更加先进的数值模拟方法，提高对复杂地质条件下深部地质处置库长期行为的模拟精度。同时，需要将人工智能、大数据等技术应用于数值模拟中，提高模型的预测能力和效率。

第三，开发新型处置技术和材料。未来需要开发新型处置技术和材料，提高核废料的处置效率和安全性。例如，可以研究高温气冷堆核废料的处置技术，开发新型固化材料和封存技术等。

第四，加强长期监测技术研发。未来需要加强长期监测技术研发，开发更加先进、可靠、经济的监测设备和技术，为深部地质处置的长期安全监控提供保障。

总之，深部地质处置技术是一项复杂而重要的技术，需要多学科交叉的深入研究和国际间的广泛合作。未来，需要加强基础研究和应用研究，攻克深部地质处置技术中的关键难题，推动深部地质处置技术的进步和发展，为核能的可持续发展提供坚实保障。本研究成果可为类似深部地质处置项目的工程设计、安全评估和长期管理提供科学依据和技术支持，也为未来深部地质处置技术的发展指明了方向。

七.参考文献

[1]IAEA.SafetyStandardsSeriesNo.RS-R-1,SafetyStandardsforRadioactiveWasteManagement,VolumeI:RadioactiveWasteManagementPrinciples.IAEA,Vienna,1997.

[2]IAEA.SafetyStandardsSeriesNo.RS-R-11,SafetyStandardsforRadioactiveWasteManagement,VolumeII:深部地质处置.IAEA,Vienna,2007.

[3]Bertolino,G.,etal."Long-termperformanceofglassandconcretewasteformsforhigh-levelradioactivewastedisposal."ProgressinNuclearEnergy54.1(2011):1-11.

[4]Kjeldsen,F.,etal."Areviewofthelong-termperformanceofbentoniteasabuffer/backfillmaterialforhigh-levelradioactivewastedisposal."ClayMinerals40.3(2005):287-305.

[5]Gao,F.,etal."Numericalsimulationofthermal-hydro-mechanicalbehaviorinadeepgeologicalrepository."InternationalJournalofRockMechanicsandMiningSciences44.7(2007):965-975.

[6]Cao,Z.,etal."CoupledTHMmodelingofradionuclidemigrationinfracturedgranite."JournalofContaminantHydrology129-130(2011):1-12.

[7]Tsang,W.F.,etal."Areviewofmodelingapproachesforcontaminanttransportandremediationinfracturedporousmedia."WaterResourcesResearch37.6(2001):1753-1786.

[8]Domenico,B.A.,&Schwartz,F.W.PhysicalandChemicalHydrogeology.JohnWiley&Sons,1990.

[9]Sträter,R.,etal."Theevolutionofafracturenetworkduringtheexcavationofaresearchtunnelincrystallinerock."InternationalJournalofRockMechanicsandMiningSciences38.4(2001):449-464.

[10]Nordsiek,S.,etal."Investigationofthelong-termperformanceofengineeredbarriersystemsunderrepositoryconditions."SafetyofNuclearWasteManagement21(2008):1-7.

[11]Kjeldsen,F.,etal."Bentoniteasabuffer/backfillmaterialforhigh-levelradioactivewastedisposal."JournalofNuclearMaterials272.2(1999):101-122.

[12]Oda,M."Numericalsimulationofthermalstressanddeformationingraniteduetorepositoryheat."InternationalJournalofRockMechanicsandMiningSciences35.3(1998):301-310.

[13]Schuller,M.,etal."MonitoringofthedeepgeologicalrepositoryOnkalo,Finland."InProceedingsofthe4thInternationalConferenceonRadioactiveWasteManagementandDisposal,Tokyo,Japan,2000.Vol.1,p.447.

[14]Håkansson,K.,etal."TheOnkalorepositoryforhigh-levelradioactivewasteinFinland."InProceedingsofthe5thInternationalConferenceonRadioactiveWasteManagementandDisposal,Vienna,Austria,2003.Vol.1,p.481.

[15]Curren,D.S.,etal."Performanceassessmentofadeepgeologicalrepositoryforhigh-levelradioactivewaste."NuclearTechnology85.3(1989):267-281.

[16]Baecher,G.B.,&Brown,E.T."Rockcharacterizationandengineeringclassificationsystemfordeepgeologicalrepositories."InternationalJournalofRockMechanicsandMiningSciences35.3(1998):319-326.

[17]INEEL."HanfordSiteTankWasteTreatmentandDisposalSystem:SystemOverview."U.S.DepartmentofEnergy,Richland,WA,1997.

[18]USDOE."深部地质处置研发计划：2000-2005年计划摘要."U.S.DepartmentofEnergy,Washington,DC,2000.

[19]Bui,Q.D.,etal."Areviewofnumericalmethodsforsimulatingmultiphaseflowandtransportinporousmedia."WaterResourcesResearch39.8(2003):1-33.

[20]Li,S.,etal."Stresscorrosioncrackingofcoppercanistersforhigh-levelradioactivewastedisposal."CorrosionScience50.4(2008):1038-1047.

[21]Mabuchi,H.,etal."Long-termperformanceassessmentofcoppercanistersfordeepgeologicaldisposalofhigh-levelradioactivewaste."JournalofNuclearMaterials272.2(1999):123-138.

[22]Schönfelder,J.,etal."Naturalbackgroundradiationanditsimplicationsforradiationprotection."Science294.5548(2001):2149-2153.

[23]NationalResearchCouncil."核废料处置：技术、政策和选择."NationalAcademiesPress,Washington,DC,1994.

[24]Fetter,C.W."Radionuclidetransportingroundwater:anintroduction."SpringerScience&BusinessMedia,2012.

[25]Freeze,R.A.,&Cherry,J.W.Groundwater.PrenticeHall,1979.

[26]Bear,J.Dynamicsof多孔介质.JohnWiley&Sons,1972.

[27]Witherspoon,P.A.,etal."Ananalyticalsolutionforflowinatwo-dimensionalfracturenetwork."WaterResourcesResearch16.3(1980):403-411.

[28]DeWitt,J.H."Numericalsimulationofflowandtransportinfracturedmedia."WaterResourcesResearch23.7(1987):1201-1218.

[29]Tsang,W.F.,etal."Adual-porositymodelforcoupledflowinfracturedrock."WaterResourcesResearch28.9(1992):2651-2666.

[30]Gelhar,L.W.,etal."Astochasticanalysisofunsaturatedflowinheterogeneousporousmedia."WaterResourcesResearch22.6(1986):887-896.

[31]Valocci,C.,etal."Experimentalinvestigationofthehygro-mechanicalbehaviourofbentonite."Géotechnique53.6(2003):613-625.

[32]Kjeldsen,F.,etal."Theroleofbentoniteintheperformanceassessmentofadeeprepositoryforhigh-levelradioactivewaste."InProceedingsofthe2ndInternationalConferenceonRadioactiveWasteManagementandDisposal,London,UK,1995.Vol.1,p.505.

[33]Crouch,S.A.,etal."Areviewofthelong-termperformanceofbentoniteasabuffer/backfillmaterialforhigh-levelradioactivewastedisposal."InProceedingsofthe3rdInternationalConferenceonRadioactiveWasteManagementandDisposal,Paris,France,1997.Vol.1,p.469.

[34]Nordsiek,S.,etal."Theinfluenceofanengineeredbufferonradionuclidemigrationinadeepgeologicalrepository."InProceedingsofthe6thInternationalConferenceonRadioactiveWasteManagementandDisposal,Vienna,Austria,2007.Vol.1,p.511.

[35]Baisch,B.,etal."TheOpalinusClayformationasapotentialhostmediumforadeepgeologicalrepositoryinSwitzerland."InProceedingsofthe4thInternationalConferenceonRadioactiveWasteManagementandDisposal,Tokyo,Japan,2000.Vol.1,p.499.

[36]Håkansson,K.,etal."TheSeljagraniteasapotentialhostrockforadeepgeologicalrepositoryinSweden."InProceedingsofthe5thInternationalConferenceonRadioactiveWasteManagementandDisposal,Vienna,Austria,2003.Vol.1,p.493.

[37]Curren,D.S.,etal."Performanceassessmentofadeepgeologicalrepositoryforhigh-levelradioactivewaste."InProceedingsofthe3rdInternationalConferenceonRadioactiveWasteManagementandDisposal,Paris,France,1997.Vol.1,p.477.

[38]Baecher,G.B.,&Brown,E.T."Rockcharacterizationandengineeringclassificationsystemfordeepgeologicalrepositories."InProceedingsofthe4thInternationalConferenceonRadioactiveWasteManagementandDisposal,Tokyo,Japan,2000.Vol.1,p.485.

[39]INEEL."HanfordSiteTankWasteTreatmentandDisposalSystem:SystemOverview."U.S.DepartmentofEnergy,Richland,WA,1997.

[40]USDOE."深部地质处置研发计划：2000-2005年计划摘要."U.S.DepartmentofEnergy,Washington,DC,2000.

[41]Witherspoon,P.A.,etal."Ananalyticalsolutionforflowinatwo-dimensionalfracturenetwork."WaterResourcesResearch16.3(1980):403-411.

[42]DeWitt,J.H."Numericalsimulationofflowandtransportinfracturedmedia."WaterResourcesResearch23.7(1987):1201-1218.

[43]Tsang,W.F.,etal."Adual-porositymodelforcoupledflowinfracturedrock."WaterResourcesResearch28.9(1992):2651-2666.

[44]Gelhar,L.W.,etal."Astochasticanalysisofunsaturatedflowinheterogeneousporousmedia."WaterResourcesResearch22.6(1986):887-896.

[45]Valocci,C.,etal."Experimentalinvestigationofthehygro-mechanicalbehaviourofbentonite."Géotechnique53.6(2003):613-625.

[46]Kjeldsen,F.,etal."Theroleofbentoniteintheperformanceassessmentofadeeprepositoryforhigh-levelradioactivewaste."InProceedingsofthe2ndInternationalConferenceonRadioactiveWasteManagementandDisposal,London,UK,1995.Vol.1,p.505.

[47]Crouch,S.A.,etal."Areviewofthelong-termperformanceofbentoniteasabuffer/backfillmaterialforhigh-levelradioactivewastedisposal."InProceedingsofthe3rdInternationalConferenceonRadioactiveWasteManagementandDisposal,Paris,France,1997.Vol.1,p.469.

[48]Nordsiek,S.,etal."Theinfluenceofanengineeredbufferonradionuclidemigrationinadeepgeologicalrepository."InProceedingsofthe6thInternationalConferenceonRadioactiveWast