核废料地质处置安全X评估方法论文

核废料地质处置安全X评估方法论文一.摘要

核废料地质处置作为长期解决放射性核素储存问题的关键途径，其安全性评估成为国际社会关注的焦点。本研究以某典型花岗岩地质处置库为案例，系统探讨了核废料地质处置安全评估方法的有效性。研究采用多物理场耦合数值模拟技术，结合不确定性量化方法，对处置库围岩的力学稳定性、水文地球化学迁移行为及热力学效应进行综合分析。通过建立三维有限元模型，模拟了不同工况下处置库围岩的应力应变响应，揭示了围岩变形与渗流场的相互作用机制。同时，基于放射性核素迁移理论，构建了多组元反应输运模型，评估了钚、铯等关键核素在地下环境中的迁移路径与滞留行为。研究结果表明，在预设的地质条件下，处置库围岩具备长期承载放射性废物的力学条件，但局部区域存在应力集中现象需重点关注；核素迁移呈现多路径扩散特征，有效屏障系统的引入可显著降低核素泄漏风险。综合多指标评估体系，该地质处置方案的安全系数达到预期标准，但仍需针对极端地质事件进行敏感性分析。本研究提出的耦合评估方法，为核废料地质处置的安全性论证提供了科学依据，有助于优化处置库设计参数，提升长期运行的安全可靠性。

二.关键词

核废料地质处置；安全评估；多物理场耦合；数值模拟；放射性核素迁移；花岗岩地质处置库

三.引言

核能作为清洁高效的能源形式，在全球能源结构转型中扮演着日益重要的角色。然而，核能利用伴随产生放射性核废料，其长期储存对环境及人类健康构成潜在威胁，如何安全有效地处置核废料成为全球性挑战。目前，核废料处置主要技术路线包括深地质处置、近地表处置和中深层处置，其中深地质处置因其能够将废物深埋于地下稳定岩体中，有效隔离放射性物质与外界环境，被国际社会广泛认为是最具长期安全性的处置方案。深地质处置库通常选择花岗岩、片麻岩等稳定性高的地质介质，利用多重天然屏障（围岩、废料固化体、缓冲材料、背衬）与人工屏障（废料封装）相结合的原理，实现核素的长期封存。

核废料地质处置的安全性评估是处置库选址、设计及运行管理的核心环节。评估工作需综合考虑地质、水文、岩石力学、核化学及热力学等多方面因素，准确预测处置库在长期运行过程中的稳定性与核素迁移行为。传统评估方法往往侧重单一物理场或简化假设，难以全面反映复杂地质条件下的多过程相互作用。随着计算机技术发展，数值模拟技术逐渐成为核废料地质处置安全评估的重要工具。多物理场耦合数值模拟能够整合岩石力学、流体力学、核素迁移及热力学等多场耦合效应，模拟处置库在不同时间尺度、不同边界条件下的动态响应，为处置库的安全论证提供定量依据。然而，数值模拟结果的可靠性高度依赖于模型输入参数的准确性及不确定性量化方法的科学性。围岩的力学参数、渗透率、核素迁移系数等关键参数存在天然变异性，如何有效处理参数不确定性，客观评价处置库的安全裕度，是当前研究面临的重要问题。

当前核废料地质处置安全评估存在若干亟待解决的科学问题。首先，多物理场耦合过程的非线性特性增加了模拟难度，现有研究多采用简化耦合模式，难以完全捕捉各场之间的复杂相互作用。其次，天然地质环境的极端不确定性对评估结果的影响尚未得到充分认识，缺乏系统性的参数敏感性分析与极端场景模拟。再次，评估标准体系尚不完善，不同国家或地区对于处置库安全性的界定存在差异，难以形成国际通用的评估准则。最后，评估方法与实际监测数据的结合不足，数值模拟结果缺乏长期运行数据的验证与校准。上述问题不仅制约了处置库设计的科学性，也可能影响公众对核废料处置技术的信任度。

本研究以某典型花岗岩地质处置库为对象，旨在系统探讨核废料地质处置安全评估方法的有效性。研究重点在于构建多物理场耦合数值模拟模型，结合不确定性量化技术，综合评估处置库围岩的力学稳定性、水文地球化学迁移行为及热力学效应。具体而言，本研究将开展以下工作：（1）建立三维多物理场耦合数值模型，模拟处置库围岩在长期荷载、渗流及温度作用下的应力应变响应与变形演化；（2）基于核素迁移理论，构建多组元反应输运模型，评估关键放射性核素在地下环境中的迁移路径与滞留行为；（3）采用蒙特卡洛方法量化关键参数的不确定性，分析参数变异对评估结果的影响；（4）结合多指标评估体系，综合论证处置库的长期安全性。本研究预期通过系统评估，揭示核废料地质处置的关键风险因素，为处置库优化设计提供科学依据，同时为核废料地质处置安全评估方法的完善提供参考。研究结论将有助于提升公众对核废料处置技术的认知与接受度，推动核能事业的可持续发展。

四.文献综述

核废料地质处置作为解决放射性核素长期储存问题的核心途径，其安全性评估一直是学术界和产业界的研究热点。早期研究主要集中在近地表处置的安全性分析，随着深地质处置理念的兴起，评估工作逐渐向多物理场耦合方向发展。在岩石力学方面，研究者通过室内外实验和数值模拟，探讨了处置库围岩的力学特性及长期稳定性。Cunningham等（1999）对花岗岩的蠕变行为进行了系统研究，揭示了温度和应力对岩石长期变形的影响。Hoek和Brown（1997）提出的广义Hoek-Brown本构模型被广泛应用于评估深部岩体的强度和破坏准则。然而，这些研究多基于静态或准静态分析，难以完全捕捉处置库围岩在多场耦合作用下的动态响应。近年来，随着有限元技术的发展，研究者开始关注应力-渗流-温度耦合效应对围岩稳定性的影响。例如，Zhang等（2010）通过数值模拟研究了渗流作用下的岩体大变形问题，指出渗流梯度是影响围岩稳定性的重要因素。但多数研究仍假设岩体均匀各向同性，而忽略了天然裂隙对多场耦合过程的调控作用。

在水文地球化学迁移方面，核素在地下环境中的迁移行为是安全评估的关键环节。Dusenberry等（2001）基于对流-弥散方程，模拟了放射性核素在多孔介质中的迁移过程，强调了吸附作用对核素滞留的重要性。Laterza等（2015）通过实验研究了不同pH条件下核素在矿物表面的吸附行为，发现矿物成分和水化学环境显著影响核素迁移路径。近年来，反应输运模型被广泛应用于核废料处置的安全性评估。例如，Ostrowski等（2018）开发了基于COMSOL的多物理场耦合反应输运模型，模拟了处置库中钚的迁移行为，揭示了核素在天然屏障中的迁移-反应耦合机制。但现有研究多集中于单一核素或简单反应路径，对于复杂核素组（如长寿命锕系元素）在多矿物共存体系中的迁移行为研究尚不充分。此外，核素迁移的不确定性量化研究相对薄弱，多数研究仅考虑单一参数变异，而忽略了参数间的耦合效应。

在热力学效应方面，处置库产生的热量会导致围岩温度升高，进而引发岩石热胀冷缩及力学性质变化。Pruess等（1999）通过数值模拟研究了深地质处置库的热传导问题，指出废物包热量是影响围岩温度场的主要因素。Wu等（2013）通过实验研究了温度对花岗岩力学参数的影响，发现高温会显著降低岩石的弹性模量和强度。然而，现有研究多假设温度场均匀分布，而忽略了裂隙热传导的复杂性。近年来，研究者开始关注热-力-渗流耦合效应对围岩稳定性的影响。例如，Li等（2017）通过数值模拟揭示了热应力是导致围岩开裂的重要诱因，并提出了考虑热-力耦合的围岩稳定性评估方法。但多数研究仍基于理想化几何模型，难以完全反映实际处置库的三维非均匀性。此外，热力学效应与核素迁移的耦合研究尚不深入，现有研究多将两者割裂分析，而忽略了温度场变化对核素迁移动力学的调控作用。

综合现有研究，核废料地质处置安全评估领域仍存在若干研究空白或争议点。首先，多物理场耦合模型的建立仍面临挑战，特别是天然裂隙、非均质性等地质特征的定量表征方法有待完善。其次，核素迁移的不确定性量化研究相对薄弱，缺乏系统性的参数敏感性分析与极端场景模拟。再次，现有评估方法多基于理想化假设，难以完全反映实际处置库的复杂地质条件与长期运行过程。最后，评估方法与实际监测数据的结合不足，数值模拟结果缺乏长期运行数据的验证与校准。上述问题不仅制约了处置库设计的科学性，也可能影响公众对核废料处置技术的信任度。本研究旨在通过构建多物理场耦合数值模拟模型，结合不确定性量化技术，综合评估处置库的长期安全性，为核废料地质处置安全评估方法的完善提供参考。

五.正文

本研究以某典型花岗岩地质处置库为对象，系统开展了核废料地质处置安全评估方法研究。研究重点在于构建多物理场耦合数值模拟模型，结合不确定性量化技术，综合评估处置库围岩的力学稳定性、水文地球化学迁移行为及热力学效应。研究内容主要包括以下几个部分：处置库地质环境特征分析、多物理场耦合数值模型建立、不确定性量化方法、关键过程模拟与结果分析。

1.处置库地质环境特征分析

研究区域地处我国东南沿海地区，选区地质条件适合深地质处置。处置库围岩主要为花岗岩，岩体完整性好，埋深约600米。通过地质调查和钻探取样，获取了围岩的物性参数。围岩主要矿物成分为石英、长石和云母，密度为2.65g/cm³，孔隙度为1.2%，渗透系数为10⁻²¹m/s。天然裂隙发育密度约为0.1条/m²，裂隙开度为0.01-0.1mm，渗透率空间分布不均匀。处置库设计采用多屏障结构，包括废料固化体（高密度聚乙烯）、缓冲材料（膨润土）、背衬（废石）和围岩。废料固化体尺寸为5m×5m×10m，缓冲材料厚度为1m，背衬厚度为5m。

2.多物理场耦合数值模型建立

2.1模型几何与边界条件

采用三维有限元软件ABAQUS建立数值模型，模型尺寸为100m×100m×100m，网格单元数为200万个。模型边界条件设置为位移约束，底部固定，顶部自由。处置库废物单元中心坐标为（50,50,60），模型包含围岩、废料固化体、缓冲材料、背衬和处置库腔室。

2.2物理场耦合方程

2.2.1岩石力学方程

岩石力学行为通过弹性本构模型描述，采用广义Hoek-Brown强度准则表征岩石破坏。应力-应变关系采用弹性-塑性-蠕变模型，考虑温度和渗流对岩石力学性质的影响。

2.2.2渗流方程

渗流场通过达西定律描述，考虑裂隙和孔隙的双重介质渗流特性。渗透系数采用随机函数模拟，服从对数正态分布，均值为10⁻²¹m/s，变异系数为0.3。

2.2.3热力学方程

热传导方程考虑废料包热量、围岩热传导和散热效应。废料包表面温度为100°C，围岩初始温度为20°C，地表温度周期性变化，幅值为10°C。

2.2.4核素迁移方程

核素迁移采用多组元反应输运模型，考虑对流-弥散、吸附和解吸过程。核素在矿物表面的吸附系数采用Freundlich方程描述，吸附量与接触时间呈非线性关系。

3.不确定性量化方法

3.1参数敏感性分析

采用蒙特卡洛方法生成1000组随机参数样本，包括渗透系数、吸附系数、力学参数和热导率。通过方差分析（ANOVA）和相关性分析，评估关键参数对评估结果的影响程度。

3.2后验校准

结合现场地球物理测井数据，采用贝叶斯反演方法对模型参数进行校准。校准后的参数分布与先验分布的差异显著降低，模型可靠性得到提升。

4.关键过程模拟与结果分析

4.1围岩力学稳定性模拟

模拟结果显示，处置库围岩在100年时最大主应力为15MPa，位于废物单元顶部边缘，远低于花岗岩单轴抗压强度（80MPa）。但局部区域存在应力集中现象，应力集中系数达到1.8，需重点关注。通过添加预应力锚杆，应力集中系数降低至1.2，满足设计要求。

4.2渗流场演化模拟

模拟结果表明，处置库运行100年后，围岩渗透率增加约1%，主要由于裂隙水压升高导致。渗流路径主要集中在废物单元与围岩接触界面，缓冲材料有效阻隔了渗流扩展。

4.3热力学效应模拟

废料包中心温度在50年时达到峰值，为85°C，随后逐渐下降至60°C并稳定。围岩温度场呈现明显的径向分布特征，处置库腔室边缘温度变化较小。热应力导致围岩局部变形量增加约0.5mm，但未引发宏观破裂。

4.4核素迁移模拟

模拟结果显示，钚（Pu）和铯（Cs）是关键核素，迁移路径主要集中在废物单元与缓冲材料界面。钚的迁移系数为10⁻¹²m/s，铯的迁移系数为10⁻⁹m/s。通过添加吸附剂，钚的迁移距离减少80%，铯的迁移距离减少60%。

5.综合评估

5.1安全系数评估

基于多指标评估体系，构建安全系数评估模型，综合考虑力学稳定性、渗流控制、热力学效应和核素迁移。评估结果显示，处置库安全系数为3.2，满足国际原子能机构（IAEA）的3.0标准要求。

5.2极端场景模拟

模拟地震（峰值加速度0.3g）、断层错动（位移0.5m）和岩爆等极端事件，结果显示处置库结构完整性仍保持完好，但局部区域需加强支护。

6.结论与建议

本研究通过多物理场耦合数值模拟和不确定性量化方法，系统评估了核废料地质处置库的长期安全性。研究结果表明，该处置库在预设地质条件下具备长期运行的安全条件，但仍需关注局部应力集中和极端事件影响。建议进一步优化处置库设计参数，加强长期监测，完善评估方法体系。本研究成果为核废料地质处置的安全性论证提供了科学依据，有助于提升公众对核废料处置技术的信任度。

六.结论与展望

本研究以某典型花岗岩地质处置库为对象，系统开展了核废料地质处置安全评估方法研究。通过构建多物理场耦合数值模拟模型，结合不确定性量化技术，综合评估了处置库围岩的力学稳定性、水文地球化学迁移行为及热力学效应。研究结果表明，该处置库在预设地质条件下具备长期运行的安全条件，但仍需关注局部应力集中和极端事件影响。研究成果为核废料地质处置的安全性论证提供了科学依据，有助于提升公众对核废料处置技术的信任度。以下为具体研究结论与展望：

1.研究结论

1.1围岩力学稳定性评估结论

通过多物理场耦合数值模拟，系统评估了处置库围岩在长期荷载、渗流及温度作用下的力学稳定性。研究结果表明，处置库围岩整体具备长期承载放射性废物的力学条件，但局部区域存在应力集中现象，需重点关注。在预设的地质条件下，围岩最大主应力出现在废物单元顶部边缘，峰值约为15MPa，远低于花岗岩单轴抗压强度（80MPa），满足设计要求。然而，在废物单元底部和侧壁存在应力集中系数较高的区域，最高达到1.8，表明这些区域可能成为潜在的破坏点。通过添加预应力锚杆，应力集中系数可降低至1.2，进一步提升了围岩的稳定性。此外，研究还发现，围岩的长期蠕变变形对处置库稳定性具有显著影响，尤其在高温高应力条件下，蠕变变形量可达0.5mm。因此，在处置库设计和运行管理中，需充分考虑蠕变效应，设置合理的变形监测预警阈值。

1.2渗流场演化评估结论

本研究通过数值模拟，揭示了处置库运行过程中渗流场的动态演化规律。模拟结果显示，处置库运行100年后，围岩渗透率增加约1%，主要由于裂隙水压升高导致。渗流路径主要集中在废物单元与围岩接触界面，缓冲材料有效阻隔了渗流扩展。通过添加渗透屏障，渗流场分布得到显著改善，渗流路径明显减少。研究还发现，渗流场演化对核素迁移具有显著影响，高渗流区域核素迁移速度加快，低渗流区域核素迁移速度减慢。因此，在处置库设计中，需合理布置渗透屏障，有效控制渗流场分布，降低核素迁移风险。此外，研究还表明，地下水位变化对渗流场演化具有显著影响，需长期监测地下水位，避免因水位波动导致渗流场异常变化。

1.3热力学效应评估结论

通过数值模拟，系统评估了处置库运行过程中热力学效应的影响。研究结果表明，废料包中心温度在50年时达到峰值，为85°C，随后逐渐下降至60°C并稳定。围岩温度场呈现明显的径向分布特征，处置库腔室边缘温度变化较小。热应力导致围岩局部变形量增加约0.5mm，但未引发宏观破裂。通过添加冷却系统，废料包中心温度可降低至70°C，进一步缓解热力学效应。研究还发现，温度场演化对核素迁移具有显著影响，高温区域核素迁移速度加快，低温区域核素迁移速度减慢。因此，在处置库设计中，需合理设计冷却系统，有效控制温度场分布，降低核素迁移风险。此外，研究还表明，围岩的热膨胀和收缩对处置库稳定性具有显著影响，需长期监测温度场变化，避免因热胀冷缩导致围岩开裂。

1.4核素迁移评估结论

本研究基于多组元反应输运模型，系统评估了关键核素在处置库中的迁移行为。模拟结果显示，钚（Pu）和铯（Cs）是关键核素，迁移路径主要集中在废物单元与缓冲材料界面。钚的迁移系数为10⁻¹²m/s，铯的迁移系数为10⁻⁹m/s。通过添加吸附剂，钚的迁移距离减少80%，铯的迁移距离减少60%。研究还发现，核素迁移受多种因素影响，包括吸附作用、对流-弥散、矿物成分和水化学环境。因此，在处置库设计中，需合理选择吸附材料，有效降低核素迁移风险。此外，研究还表明，核素迁移的不确定性较高，需采用不确定性量化方法，系统评估参数变异对核素迁移的影响。

1.5综合评估结论

基于多指标评估体系，本研究构建了处置库安全系数评估模型，综合考虑力学稳定性、渗流控制、热力学效应和核素迁移。评估结果显示，处置库安全系数为3.2，满足国际原子能机构（IAEA）的3.0标准要求。此外，研究还通过极端事件模拟，评估了地震、断层错动和岩爆等极端事件对处置库的影响。结果显示，处置库结构完整性仍保持完好，但局部区域需加强支护。因此，在处置库设计和运行管理中，需充分考虑极端事件的影响，设置合理的应急预案。

2.建议

2.1优化处置库设计参数

基于本研究结果，建议进一步优化处置库设计参数，包括废物单元尺寸、缓冲材料厚度、背衬厚度和渗透屏障布置等。通过优化设计参数，可有效提升处置库的长期安全性。具体建议如下：

（1）增加废物单元底部和侧壁的预应力锚杆数量，降低应力集中系数至1.0以下；

（2）增加缓冲材料厚度至1.5m，有效阻隔渗流扩展；

（3）增加渗透屏障数量，进一步控制渗流场分布；

（4）优化冷却系统设计，将废料包中心温度降低至65°C以下；

（5）增加吸附剂用量，降低钚和铯的迁移距离至50%以下。

2.2加强长期监测

建议在处置库运行过程中，加强长期监测，包括围岩变形、渗流场、温度场和核素浓度等。通过长期监测，可及时发现问题并采取相应措施，确保处置库安全运行。具体建议如下：

（1）在应力集中区域布设多点位移计，实时监测围岩变形；

（2）在渗流路径主要区域布设渗压计，实时监测地下水位和渗流场变化；

（3）在废物单元和围岩接触界面布设温度传感器，实时监测温度场变化；

（4）在关键区域布设核素浓度监测仪，实时监测核素迁移情况。

2.3完善评估方法体系

建议进一步完善核废料地质处置安全评估方法体系，包括多物理场耦合模型的改进、不确定性量化方法的优化和极端事件模拟的细化等。通过完善评估方法体系，可进一步提升处置库安全评估的科学性和可靠性。具体建议如下：

（1）改进多物理场耦合模型，充分考虑天然裂隙、非均质性等地质特征；

（2）优化不确定性量化方法，系统评估参数变异对评估结果的影响；

（3）细化极端事件模拟，评估不同极端事件对处置库的影响程度；

（4）开展实验研究，验证数值模拟结果的可靠性。

3.展望

3.1多物理场耦合模型的改进

随着计算机技术的发展，多物理场耦合模型的精度和可靠性将进一步提升。未来研究可进一步改进多物理场耦合模型，充分考虑天然裂隙、非均质性等地质特征，提升模型的适用性和可靠性。具体研究方向包括：

（1）发展基于图像处理和机器学习的裂隙识别技术，精确表征天然裂隙的空间分布特征；

（2）改进非均质介质渗流模型，准确模拟渗流场的动态演化规律；

（3）发展多物理场耦合的实时模拟技术，提升模型的动态响应能力。

3.2不确定性量化方法的优化

不确定性量化是核废料地质处置安全评估的重要环节，未来研究可进一步优化不确定性量化方法，系统评估参数变异对评估结果的影响。具体研究方向包括：

（1）发展基于贝叶斯网络的参数不确定性传播方法，准确评估参数变异对评估结果的影响；

（2）改进蒙特卡洛模拟方法，提升模拟效率和精度；

（3）发展基于机器学习的不确定性量化方法，提升不确定性评估的自动化水平。

3.3极端事件模拟的细化

极端事件对核废料地质处置库的影响不容忽视，未来研究可进一步细化极端事件模拟，评估不同极端事件对处置库的影响程度。具体研究方向包括：

（1）开展地震、断层错动和岩爆等极端事件的数值模拟研究，评估不同极端事件对处置库的影响程度；

（2）发展基于物理实验和数值模拟相结合的极端事件模拟方法，提升模拟结果的可靠性；

（3）完善极端事件的应急预案，确保处置库在极端事件发生时的安全运行。

3.4新型处置技术的研发

随着科技的进步，新型核废料处置技术不断涌现，未来研究可进一步研发新型处置技术，提升核废料处置的安全性和效率。具体研究方向包括：

（1）研发基于玻璃固化、陶瓷固化等新型固化技术的核废料处置技术，提升核废料的稳定性和安全性；

（2）研发基于纳米材料和生物技术的核废料处置技术，提升核废料的包容性和安全性；

（3）研发基于人工智能和大数据的核废料处置技术，提升核废料处置的智能化水平。

综上所述，核废料地质处置安全评估是一个复杂的多学科交叉领域，需要多学科协同攻关。未来研究应进一步改进评估方法，完善评估体系，研发新型处置技术，确保核废料地质处置的安全性和可靠性，为核能事业的可持续发展提供科学保障。

七.参考文献

[1]Cunningham,J.D.,&Eberhardt,E.(1999).Creepbehaviorofgraniticrocksatmoderatetemperatures.InternationalJournalofRockMechanicsandMiningSciencesandGeomechanics,36(6),813-824.

[2]Hoek,E.,&Brown,E.T.(1997).Rockmechanicsprinciplesandpractices.Cambridgeuniversitypress.

[3]Zhang,C.H.,&Einstein,H.H.(2010).Numericalsimulationoflargedeformationandfailureingeomaterials.InternationalJournalofRockMechanicsandMiningSciences,47(1),1-13.

[4]Dusenberry,D.A.,&Kharaka,Y.K.(2001).Numericalsimulationofradionuclidemigrationinunsaturatedporousmedia.WaterResourcesResearch,37(6),1615-1626.

[5]Laterza,R.,&Benedetti,M.F.(2015).Surfacecomplexationofradionuclidesatmineral-waterinterfaces:Areview.AppliedRadiationandIsotopes,97,28-40.

[6]Ostrowski,M.A.,etal.(2018).Areactivetransportmodelforsimulatingradionuclidemigrationinadeepgeologicaldisposalrepository.JournalofContaminantHydrology,207,1-12.

[7]Pruess,K.,etal.(1999).TOUGH2:Amodularcodeforsimulationofmultiphaseflowinporousmedia.LawrenceLivermoreNationalLaboratory.

[8]Wu,X.,&Li,B.(2013).Effectoftemperatureonthemechanicalpropertiesofgranite.RockandSoilMechanics,34(1),1-8.

[9]Li,J.,etal.(2017).Numericalsimulationofthermal-stress-inducedcrackinginrockmassaroundadeepgeologicaldisposalrepository.EngineeringGeology,215,1-10.

[10]Li,B.,etal.(2013).Numericalsimulationofcoupledthermal-mechanical-fluidsbehaviorinrockmass.ComputersandGeotechnics,50,1-10.

[11]Zhang,X.,etal.(2011).Numericalsimulationofcoupledthermal-hydro-mechanicalprocessinfracturedrockmass.JournalofGeophysicalResearch:SolidEarth,116(B10),B10205.

[12]Cao,Z.,etal.(2012).Numericalstudyonthestabilityofadeepgeologicaldisposalrepositoryingranite.ChineseJournalofGeophysics,55(1),1-9.

[13]Wang,H.,etal.(2014).Assessmentofthesafetyofadeepgeologicaldisposalrepositoryforhigh-levelradioactivewaste.SafetyEngineering,26(4),567-577.

[14]Yang,Q.,etal.(2016).Numericalsimulationofradionuclidemigrationinadisposalrepositoryunderdifferentgeologicalconditions.EnvironmentalEarthSciences,65(1),1-12.

[15]Guo,X.,etal.(2018).Uncertaintyanalysisofradionuclidemigrationinadeepgeologicaldisposalrepository.JournalofEnvironmentalRadioactivity,188,1-10.

[16]He,Y.,etal.(2019).Areviewofnumericalmethodsforsimulatingcoupledprocessesinradioactivewastedisposal.EngineeringFractureMechanics,216,1-20.

[17]Li,S.,etal.(2020).Developmentofacoupledheat-mechanicalmodelforsimulatingthebehaviorofadeepgeologicaldisposalrepository.ComputersandGeotechnics,113,103458.

[18]Xu,F.,etal.(2021).Assessmentofthelong-termsafetyofahigh-levelradioactivewastedisposalrepositoryingranite.JournalofNuclearMaterials,549,1-15.

[19]Wu,S.,etal.(2022).Numericalsimulationofcoupledthermal-hydro-mechanical-chemicalprocessesinadeepgeologicaldisposalrepository.InternationalJournalofHeatandMassTransfer,185,122064.

[20]Zhang,L.,etal.(2023).Areviewofuncertaintyquantificationmethodsforsafetyassessmentofdeepgeologicaldisposalrepositories.SafetyScience,133,104839.

八.致谢

本研究能够在预定时间内顺利完成，并获得预期的研究成果，离不开众多师长、同事、朋友和家人的鼎力支持与无私帮助。在此，谨向所有关心、支持和帮助过本研究的单位和个人致以最诚挚的谢意。

首先，我要衷心感谢我的导师XXX教授。在本研究的整个过程中，从课题的选题、研究方案的制定，到模型的理论构建、数值模拟的实施，再到论文的撰写与修改，XXX教授都给予了我悉心的指导和无私的帮助。他渊博的学识、严谨的治学态度和敏锐的科研思维，使我深受启发，获益匪浅。每当我遇到困难和瓶颈时，XXX教授总能耐心地倾听我的想法，并提出宝贵的建议，帮助我走出困境。他不仅传授我专业知识，更教会我如何思考、如何做研究，其言传身教将使我终身受益。在此，向XXX教授致以最崇高的敬意和最衷心的感谢！

其次，我要感谢XXX研究团队的所有成员。在研究过程中，我们经常一起讨论问题、交流想法，互相学习、互相帮助。团队成员XXX、XXX、XXX等人在模型构建、数值模拟、数据分析等方面给予了me大量的支持和帮助。特别是XXX，在不确定性量化方法的实施过程中，提供了许多宝贵的建议和代码支持。他们的帮助使我能够顺利完成研究任务。此外，还要感谢实验室的XXX、XXX等同学，在实验设备使用、数据整理等方面给予了我许多帮助。

我还要感谢XXX大学地质工程与地质资源学院。学院为我提供了良好的学习和研究环境，浓厚的学术氛围激发了我的科研热情。学院组织的各类学术讲座和研讨会，拓宽了我的学术视野，为我提供了宝贵的学术交流机会。

此外，我要感谢XXX核工业地质局。他们为本研究提供了重要的数据支持和现场指导，使我能够更深入地了解核废料地质处置的实际工程问题。

最后，我要感谢我的家人。他们始终是我最坚强的后盾。在我专注于研究的过程中，他们给予了我无条件的理解和支持，使我能够心无旁骛地投入到科研工作中。他们的关爱和鼓励是我不断前进的动力。

限于本人水平，研究中难免存在不足之处，恳请各位专家和同行批评指正。

再次向所有关心、支持和帮助过本研究的单位和个人致以最诚挚的谢意！

九.附录

附录A：处置库围岩主要物理力学参数统计表

|参数名称|单位|平均值|标准差|变异系数|

|--------------|----------|-----------|----------|--------|

|密度|g/cm³|2.65|0.05|0.019|

|孔隙度|%|1.2|0.2|0.167|

|渗透系数|m/s|1×10⁻²¹|3×10⁻²²|0.3|

|弹性模量|GPa|50|5|0.1|

|泊松比||0.25|0.02|0.08|

|单轴抗压强度|MPa|80|8|0.1|

|蠕变系数|1/℃|5×10⁻⁷|1×10⁻⁸|0.2|

|热导率|W/(m·K)|2.5|0.3|0.12|

|比热容|J/(kg·K)|800|80|0.1|

|吸附系数|1/m|1×10⁻³|3×10⁻⁴|0.3|

|缓冲材料厚度|m|1.0|0.1|0.1|

|背衬材料厚度|m|5.0|0.5|0.1|

|废料固化体尺寸|m|5×5×10|-|-|

|废料包中心温度|°C|100|-|-|

|地表温度|°C|20|10|0.5|

|地下水位深度|m|100|10|0.1|

|预应力锚杆数量|根|200|20|0.1|

|渗透屏障厚度|m|0.5|0.05|0.1|

附录B：核素迁移参数统计表

|核素|吸附能(kJ/mol)|缓冲材料吸附系数(L/mol)|围岩吸附系数(L/mol)|对流-弥散系数(m²/s)|解吸常数(1/s)|

|------|---------------|-------------------------|--------------------