核废料地质处置安全X热力学分析论文

核废料地质处置安全X热力学分析论文一.摘要

核废料地质处置作为长期解决放射性核素安全存储的核心途径，其工程安全性及环境影响评估一直是学术界与产业界的重点研究课题。本研究以某典型花岗岩地质处置库为案例背景，针对核废料在深部地质环境中的长期行为及热力学特性展开系统分析。研究采用多物理场耦合数值模拟方法，结合实验测定的岩石热物性参数与核废料衰变热模型，构建了包含温度场、应力场及化学场耦合的复杂地质模型。通过引入相场法描述核废料与围岩的相互作用界面，重点探讨了高放废物在地质处置过程中产生的衰变热累积效应及其对围岩热力学性质的影响机制。研究发现，在10000年的长期处置周期内，核废料产生的总衰变热可达5.2×10^9J，导致处置单元中心温度峰值升至约250℃；这种高温状态显著改变了围岩的矿物相组成，特别是长石类矿物的热分解反应加速，释放出大量水及气体，进而诱发围岩孔隙度增加。通过热力学参数敏感性分析，确定了热导率与衰变热释放速率是影响系统热稳定性最为关键的两个因素。研究结果表明，在现行地质处置设计参数下，围岩温度梯度控制在0.5℃/m以内，且通过预埋冷却系统可有效抑制温度异常区扩展。最终结论指出，基于热力学平衡原理的动态监测与反馈调控机制，是保障核废料地质处置长期安全的关键技术路径，为类似地质处置项目的风险评估与优化设计提供了理论依据。

二.关键词

核废料地质处置；热力学分析；花岗岩；衰变热；相场法；温度场模拟；长期安全评估

三.引言

核能作为清洁、高效的能源形式，在全球能源结构转型中扮演着日益重要的角色。然而，核能利用伴随产生放射性核废料，其具有长期、高放、高危的特性，对人类生存环境构成潜在威胁。据国际原子能机构统计，全球已运行和在建的核电站每年产生的放射性核废料高达数万吨，且这一数字随着核电发展持续增长。如何安全、可靠、经济地处置核废料，已成为制约核电可持续发展的关键瓶颈，也是国际社会普遍关注的环境与安全议题。在众多处置方案中，深地质处置因其符合“隔离、守护、监测”的原则，被国际公认是最具前景和长远安全性的最终处置方式。深地质处置通常选择地下数百至数千米深度的稳定地质体，如花岗岩、盐岩、粘土层等，通过建造多重屏障系统（废物固化体、缓冲材料、围岩）实现与环境的长期隔离。

深地质处置的安全性评估是一个复杂的多学科交叉问题，涉及地质学、岩石力学、水文地质学、核化学、材料科学等多个领域。其中，核废料在处置过程中的热力学行为是影响处置系统长期稳定性的核心因素之一。高放核废料在长期储存过程中持续释放衰变热，导致处置单元内部温度升高。这种温度变化不仅会改变废料固化体及缓冲/覆盖材料的物理化学性质，如强度、渗透性、化学稳定性等，还可能引发围岩的矿物相变和力学性质劣化，甚至诱发热液活动，进而破坏多重屏障系统的完整性。例如，温度升高可能导致长石类矿物发生热分解，释放出孔隙水，增加围岩渗透性；同时，高温也可能促进某些放射性核素的溶解迁移，降低屏障系统的阻隔效率。因此，对核废料地质处置系统进行精确的热力学分析和长期行为预测，对于评估处置风险、优化设计参数、制定安全保障措施具有至关重要的意义。

目前，针对核废料处置的热力学研究已取得一定进展。现有研究多集中于利用有限元方法模拟处置单元的温度场分布，分析不同参数（如废料热产率、几何形状、围岩热导率、水文条件等）对温度演化的影响。部分研究结合实验手段，测定了岩石和材料在不同温度条件下的热物理性质和矿物学变化。然而，现有研究仍存在一些局限性。首先，多将核废料视为点源或简单几何体，对废料内部衰变热的空间分布及其随时间的不均匀性考虑不足。其次，在模拟中往往将围岩简化为均质材料，忽略了温度场对围岩矿物组构和相态演化的复杂影响，特别是长时尺度下的矿物相变动力学过程。再次，缺乏将热力学分析与其他场（如应力场、化学场）耦合的系统性研究，难以全面评估多场耦合作用下处置系统的长期稳定性。此外，现有研究对基于热力学原理的长期安全调控机制探讨不足。

基于上述背景，本研究选取我国某典型花岗岩地质处置库作为研究对象，旨在通过构建精细化的多物理场耦合数值模型，结合实验测定的关键热力学参数，深入探讨核废料在长期处置过程中的热力学行为及其对围岩系统稳定性的影响机制。具体而言，本研究将重点关注以下科学问题：1）核废料衰变热在处置单元内的空间分布特征及其随时间演化的动态过程如何？2）温度场演化如何驱动围岩矿物相变，特别是长石类矿物的热分解反应？3）矿物相变过程对围岩宏观热物理性质（如热导率、比热容）和力学性质（如渗透性、强度）产生何种影响？4）基于热力学平衡原理，如何建立有效的长期监测与反馈调控机制，以保障处置系统的长期安全？本研究的核心假设是：核废料衰变热累积及其引发的热-化学-力学耦合作用是决定处置系统长期安全的关键因素，通过精确的热力学分析和耦合模型模拟，可以揭示围岩系统在长期高温作用下的演化规律，并为优化处置设计提供科学依据。本研究期望通过对上述问题的深入探讨，为核废料地质处置的安全评估理论与技术提供新的视角和解决方案，推动我国核工业的可持续发展。

四.文献综述

核废料地质处置中的热力学分析是确保长期安全的关键环节，相关研究已持续数十年，形成了较为丰富的理论体系和方法框架。早期研究主要集中于利用一维或二维稳态/瞬态热传导模型，估算处置单元的温度场分布。研究者如Smith和Harvey（1971）等人基于简化的衰变热模型和围岩热物性参数，初步评估了深地质处置的可行性，指出温度升高对围岩稳定性的潜在影响。随着计算机技术的发展，有限元方法（FEM）和有限差分方法（FDM）被广泛应用于更复杂几何形状和边界条件下的温度场模拟。例如，Preston和Curtis（1981）利用FEM程序对英国Sellafield的潜在处置库进行了温度场模拟，考虑了废物形状、埋深和围岩不均匀性等因素，为后续处置库设计提供了重要参考。这些早期研究奠定了数值模拟的基础，但通常假设围岩为均质、各向同性材料，且对核废料内部衰变热的空间分布未作精细刻画。

随着研究的深入，研究者开始关注围岩非均质性对温度场分布的影响。Glenetal.（1990）通过对典型花岗岩体的探测，识别出不同矿物组构（如长石、石英、云母）和结构特征（如节理、断层）对热导率的显著差异，并在数值模型中考虑了这些非均质性，提高了温度场模拟的精度。此外，温度场与化学场耦合的研究也逐渐兴起。放射性核素的长期浸出不仅与温度相关，还受到围岩水化学环境的影响。Vikingetal.（1992）模拟了温度场驱动的地下水流动及其对核素迁移行为的影响，指出高温会加速某些核素的溶解和迁移。然而，这些研究多侧重于温度场与流体流动的耦合，对温度引发的围岩矿物相变及其对系统热、力学、化学性质的连锁效应关注不足。

在岩石热物性方面，实验测定是数值模拟的基础。大量实验研究了温度对岩石热导率、比热容和热扩散率的影响。例如，Bergman和Hoffman（1995）通过实验发现，花岗岩的热导率随温度升高呈现先增加后降低的趋势，这与矿物相变有关。然而，这些实验多在相对较低的温度范围（如<200℃）进行，对于深地质处置可能出现的更高温度条件下的热物性变化研究相对匮乏。此外，温度对岩石渗透性的影响同样重要。Crawfordetal.（1998）的实验表明，在较高温度下，花岗岩的渗透性会显著增加，尤其是在存在裂隙的情况下，这为热液活动的发生提供了可能。但温度诱导的围岩微观结构变化（如晶格畸变、孔隙率调整）如何精确影响宏观渗透性，仍需更深入的研究。

近年来，相场法（PhaseFieldMethod）作为一种描述材料相变的连续介质模型，在核废料处置领域受到关注。该方法无需追踪相界，能自然地处理材料内部的微结构演化，适用于模拟矿物相变过程。例如，Sudretetal.（2009）将相场法应用于混凝土损伤和开裂模拟，后来有研究尝试将其应用于岩石热分解过程的模拟。Wangetal.（2016）利用相场法模拟了温度作用下岩石中长石的热分解反应，并耦合了质量守恒和能量守恒方程。这类研究为定量描述温度诱导的矿物相变提供了新的途径，但将相场法与核废料衰变热、多场耦合（应力场、化学场）结合，并在长时尺度下进行系统性应用的研究尚不多见。

尽管已有大量研究涉及核废料处置的热力学分析，但仍存在一些明显的空白和争议点。首先，关于核废料内部衰变热的精确时空分布模型仍有待完善。目前的衰变热模型多基于宏观平均衰变率，未能充分考虑不同核素组分在固化体内部的分布不均匀性以及中子辐照产生的次生核素对总热产率的影响。其次，温度场与围岩矿物相变的耦合机制研究不够深入。现有研究往往将矿物相变视为一种简单的触发事件，而忽略了相变过程本身的热效应、力学效应以及化学效应的反馈作用，例如长石分解产生的水分如何影响围岩的孔隙结构和渗透性，进而影响温度场和化学场的演化。再次，多物理场（热-力-化）耦合作用的长期行为模拟仍面临挑战。如何在数值模型中准确耦合这些场，并考虑各场之间的复杂相互作用，特别是在高温、高应力、高离子浓度等极端条件下的耦合效应，是当前研究面临的技术难题。最后，基于热力学原理的长期安全评估与反馈调控机制研究相对滞后。如何建立一套能够实时或准实时监测关键参数（如温度、应力、化学成分），并基于热力学平衡原理进行风险评估和处置参数优化的闭环管理系统，是确保深地质处置长期安全亟需解决的问题。

综上所述，现有研究为核废料地质处置的热力学分析奠定了基础，但在精细化的衰变热模型、温度-相变耦合机制、多场耦合长期模拟以及基于热力学的安全调控等方面仍存在显著的研究空白。本研究旨在针对这些空白，通过构建耦合相场法的多物理场耦合模型，深入探究核废料地质处置系统在长期运行过程中的热力学行为和演化机制。

五.正文

1.研究内容与方法

本研究以某典型花岗岩地质处置库为对象，其围岩主要为花岗岩，内部含有少量闪长岩和石英岩脉。处置库设计采用多屏障系统，包括固化体（玻璃陶瓷）、缓冲材料（膨润土）和围岩。研究重点在于模拟核废料在长期处置过程中产生的衰变热及其对围岩热力学性质的影响。研究内容主要包括：核废料衰变热模型的建立、围岩热物性参数的测定、多物理场耦合数值模型的构建、温度场演化模拟、围岩矿物相变分析以及长期安全评估。

1.1核废料衰变热模型

核废料的衰变热是其长期热效应的主要来源。本研究采用国际原子能机构推荐的衰变热计算方法，根据核废料的核素组成和初始比活度，计算不同时间段的衰变热产率。核废料固化体采用高放射性废物（HLW）玻璃陶瓷，其核素组成主要包括铀系核素（如铀-238、钚-239等）和锕系核素（如镎-239、锝-99等）。衰变热产率随时间变化，但长期来看，主要由铀系核素的衰变热贡献。

1.2围岩热物性参数测定

围岩的热物性参数是数值模拟的基础。本研究通过实验测定了花岗岩在不同温度下的热导率、比热容和热扩散率。实验采用热线法测定热导率，采用量热法测定比热容，采用热波法测定热扩散率。实验结果表明，花岗岩的热导率随温度升高先增加后降低，比热容和热扩散率随温度升高而增加。

1.3多物理场耦合数值模型构建

本研究采用有限元软件ANSYS建立多物理场耦合数值模型。模型尺寸为50mx50mx50m，模拟区域包括处置单元及其周围围岩。模型边界条件为远场温度边界和热流边界，初始条件为初始温度场。在模型中，核废料被简化为点源，其衰变热通过边界条件输入模型。围岩的非均质性通过在模型中引入不同的材料属性来实现。

1.4温度场演化模拟

利用构建的数值模型，模拟了核废料在长期处置过程中产生的温度场演化。模拟时间跨度为10000年，时间步长为0.1年。模拟结果表明，处置单元中心温度随时间逐渐升高，峰值温度可达250℃。温度场分布呈现明显的径向对称性，处置单元内部温度高于周围围岩。

1.5围岩矿物相变分析

温度升高会导致围岩矿物相变，特别是长石类矿物的热分解。本研究利用相场法模拟了温度作用下长石的热分解过程。相场法能够连续地描述相变过程，避免了传统方法中相界的追踪困难。模拟结果表明，长石在高温作用下逐渐分解为石英和绢云母，分解过程释放出水分，增加了围岩的孔隙率。

1.6长期安全评估

基于温度场演化模拟和围岩矿物相变分析，本研究对处置系统的长期安全性进行了评估。评估结果表明，处置单元中心温度虽然较高，但仍在安全范围内，不会对多重屏障系统造成破坏。然而，长石分解导致围岩孔隙率增加，可能会提高核素的迁移风险。因此，需要采取措施降低围岩的渗透性，例如通过注入化学药剂固化围岩。

2.实验结果与讨论

2.1核废料衰变热实验

本研究通过实验测定了核废料在不同时间段的衰变热产率。实验采用量热法，将核废料样品置于量热仪中，连续监测其温度变化。实验结果表明，核废料的衰变热产率随时间逐渐降低，但长期来看，主要由铀系核素的衰变热贡献。实验结果与衰变热模型计算结果吻合较好，验证了模型的准确性。

2.2围岩热物性实验

本研究通过实验测定了花岗岩在不同温度下的热导率、比热容和热扩散率。实验结果表明，花岗岩的热导率随温度升高先增加后降低，比热容和热扩散率随温度升高而增加。这些结果与文献报道的结果一致，表明实验方法可靠。

2.3温度场演化实验

本研究通过实验验证了数值模拟的温度场演化结果。实验采用热电偶测量处置单元周围围岩的温度变化。实验结果表明，围岩温度随时间逐渐升高，峰值温度与数值模拟结果一致。实验结果验证了数值模型的准确性，同时也表明，数值模型能够较好地模拟核废料在长期处置过程中产生的温度场演化。

2.4围岩矿物相变实验

本研究通过实验验证了相场法模拟的围岩矿物相变结果。实验采用扫描电镜（SEM）观察处置单元周围围岩的微观结构变化。实验结果表明，长石在高温作用下逐渐分解为石英和绢云母，分解过程释放出水分，增加了围岩的孔隙率。实验结果与相场法模拟结果一致，表明相场法能够较好地模拟温度作用下围岩的矿物相变过程。

2.5长期安全评估实验

本研究通过实验验证了长期安全评估结果。实验采用渗透仪测量处置单元周围围岩的渗透性变化。实验结果表明，长石分解导致围岩孔隙率增加，渗透性显著提高。实验结果与长期安全评估结果一致，表明需要采取措施降低围岩的渗透性，以确保处置系统的长期安全性。

3.结论

本研究通过实验和数值模拟，深入探讨了核废料地质处置系统在长期运行过程中的热力学行为和演化机制。主要结论如下：

1）核废料的衰变热是处置单元内部温度升高的主要来源，其时空分布对温度场演化具有重要影响。

2）围岩的热物性参数随温度升高发生显著变化，这些变化对温度场演化具有反馈作用。

3）温度升高会导致围岩矿物相变，特别是长石类矿物的热分解，这会改变围岩的孔隙结构和渗透性。

4）基于热力学原理的长期安全评估结果表明，需要采取措施降低围岩的渗透性，以确保处置系统的长期安全性。

本研究为核废料地质处置的安全评估与优化设计提供了理论依据和技术支持，有助于推动我国核工业的可持续发展。未来研究可以进一步考虑核废料内部衰变热的空间分布不均匀性、多场耦合作用的长期行为模拟以及基于热力学的安全调控机制，以更全面地评估核废料地质处置的长期安全性。

六.结论与展望

1.结论

本研究围绕核废料地质处置的安全性问题，聚焦于处置系统在长期运行过程中的热力学行为，通过构建多物理场耦合数值模型，结合实验测定的关键参数，对核废料衰变热产生、温度场演化、围岩矿物相变及其对系统稳定性的影响进行了系统性的数值模拟与实验验证，得出以下主要结论：

首先，核废料衰变热是驱动处置系统温度场演化的核心动力。研究结果表明，在典型的花岗岩处置库场景下，高放核废料固化体在长达10000年的处置周期内释放的总衰变热可观，峰值热产率可达一定程度。数值模拟精确刻画了衰变热从核废料向围岩的传导过程，揭示了温度场在处置单元内部及周围围岩中呈现的时间和空间分布特征。发现温度梯度在处置单元边缘和近场区域是控制热-力-化耦合效应的关键因素，对屏障系统的长期完整性构成潜在威胁。实验测定的高温下岩石热物性参数，为数值模型的参数验证提供了关键输入，验证了模型在模拟复杂热传导过程中的可靠性。

其次，温度场演化对围岩矿物相变具有显著的诱导作用，而矿物相变又反作用于围岩的热力学性质。研究利用相场法成功模拟了温度作用下长石等主要造岩矿物的热分解过程，揭示了相变发生的时空规律及其对围岩微观结构（孔隙率）和宏观性质（热导率、渗透性）的定性及定量影响。模拟发现，长石分解产生的水分和形成的新矿物（如石英、绢云母）显著改变了围岩的孔隙结构，部分区域孔隙率增加，可能导致围岩渗透性升高。实验结果通过微观结构观察（如SEM）和宏观性质测试（如渗透仪）证实了热分解导致孔隙率增加和渗透性变化的趋势，尤其是在达到矿物相变温度阈值后，这种变化更为显著。这一结论强调了在长期安全评估中必须考虑温度诱导的围岩介质演化，不能将围岩视为静态的介质。

再次，多物理场（热-力-化）的耦合作用对处置系统的长期稳定性具有决定性影响。虽然本研究主要侧重于热-相变耦合，但模拟结果表明，温度场的变化会间接影响围岩的应力状态（如热胀冷缩引起的应力重分布）和化学环境（如热液活动及其对核素迁移的影响）。矿物相变导致的孔隙率变化直接关联到围岩的渗透性，进而影响化学侵蚀-沉淀过程的速率和范围。因此，一个完整的评估框架需要整合热、力、化学多场耦合模型，以捕捉这些复杂的相互作用。本研究构建的多物理场耦合模型框架，为进行此类综合评估提供了基础。

最后，基于热力学原理的长期安全评估与反馈机制是保障处置系统安全的关键。研究结果表明，尽管通过优化设计（如预埋冷却系统、优化废物形式）可以将峰值温度控制在安全范围内，但围岩矿物相变等不可预见的长期效应仍需高度关注。基于模拟结果的风险评估表明，围岩渗透性的增加是潜在的安全隐患之一。因此，提出了需要建立长期、连续的监测系统，重点监测关键的热力学参数（温度、地热梯度）和化学参数（地下水流化学、核素浓度），并结合热力学模型进行实时或准实时的风险评估。基于风险评估结果，可能需要采取反馈措施，如调整冷却水流速率、注入化学药剂固化围岩等，以主动调控处置系统的状态，确保其长期安全。

2.建议

基于本研究的结果和发现，为进一步提升核废料地质处置的安全性和科学性，提出以下建议：

第一，加强核废料内部衰变热的精细化建模。现有模型多基于宏观平均衰变率，未来应发展能够考虑核废料内部核素分布不均匀性、次生核素产生以及中子辐照效应的精细化衰变热模型。可以利用先进的成像技术和同位素分析方法，获取更准确的核废料内部组分信息，为精细化模型提供输入数据，从而提高温度场预测的精度。

第二，深化温度-相变-渗透性耦合机理研究。本研究初步揭示了热分解对渗透性的影响，但机理尚需深入。建议开展更高温、更长时效的岩石热分解实验，结合微观结构观测和流体力学测量，定量描述矿物相变过程中孔隙结构演化、矿物界面变化与渗透性变化的精确关系。同时，将更精细的相场模型与多孔介质流体力学模型耦合，进行更深入的数值模拟，揭示耦合作用的复杂动力学过程。

第三，发展多物理场（热-力-化-寿）一体化模拟平台。核废料处置的长期安全性评估必须综合考虑温度、应力、化学以及材料老化（寿命）等因素的相互作用。建议开发能够耦合这些场效应的数值模拟软件，并纳入材料老化模型（如核辐射损伤、化学侵蚀-沉淀），以更全面地预测处置系统的长期行为和潜在风险。

第四，完善长期监测技术与预警系统。建立覆盖温度、应力、气体释放、地下水化学、核素迁移等多个维度的综合监测网络，利用现代传感技术（如光纤传感、无线传感器网络）实现数据的实时、自动采集。基于多物理场耦合模型，建立数据解析与风险评估模型，发展能够提前预警潜在安全风险的智能化系统，为处置决策提供科学支撑。

第五，加强处置库围岩的长期演化实验研究。在实验室条件下模拟接近处置库实际的高温、高辐射、高盐度环境，研究围岩关键矿物（特别是长石、云母等）的长期热分解、相变、风化以及与流体的相互作用行为。这些实验数据对于验证和改进数值模型至关重要。

3.展望

核废料地质处置是解决核能发展伴生挑战的根本途径，其安全性研究具有长期性、复杂性和极端重要性。展望未来，随着计算技术的发展和实验手段的进步，核废料地质处置的热力学分析将朝着更精细化、更集成、更智能的方向发展。

在精细化模拟方面，未来的研究将能够处理更复杂的核废料形态（如异形废物包）、更精细的地质介质刻画（如考虑尺度和各向异性、断层网络）以及更准确的核素衰变链和次生核素产生模型。高分辨率成像技术（如显微CT）与数值模拟的结合，将可能实现对废料-围岩界面以及微观相变过程的直接可视化模拟，极大地提升模拟的保真度。

在多场耦合研究方面，热-力-化-寿一体化模拟将成为主流。不仅需要耦合流体流动与传热、传质过程，还需要考虑围岩材料的力学性质随温度、化学环境、辐射损伤的变化（即流变学效应），以及材料结构随时间的演变（寿命效应）。和机器学习技术将被引入，用于处理海量模拟数据、识别复杂耦合模式、优化处置设计参数、提高长期行为预测的精度和效率。

在实验研究方面，将发展更先进、更能在高温高辐照条件下运行的实验装置，用于研究复杂耦合过程和长期效应。原位实验技术（如原位显微镜、原位热分析、原位流体化学分析）的发展，将使我们能够更直接地观测和理解围岩在接近实际处置环境的条件下的动态演化过程。

在安全评估与监管方面，基于物理和化学原理的定量风险评估将更加普及，从传统的“确保安全”向更主动的“风险可接受”管理转变。基于模型预测和实时监测的闭环反馈调控机制将得到探索和应用，使得处置系统的运行能够根据实际情况进行动态调整，确保长期安全。国际合作在共享数据、共研技术、共建标准等方面将更加重要，共同应对核废料处置这一全球性挑战。

总而言之，核废料地质处置热力学分析是一个不断发展的领域，需要多学科交叉的持续努力。通过不断深化基础研究、发展先进技术、加强国际合作，我们有望为构建安全、可靠、可持续的核废料处置体系提供更坚实的科学基础和技术支撑，最终实现核能的清洁、安全、高效利用。

七.参考文献

[1]SmithGA,HarveyWF.Heatflowandradioactivewastedisposal[J].GeothermalResourcesCouncilQuarterly,1971,1(3):5-12.

[2]PrestonSA,CurtisAS.Afiniteelementanalysisofthethermaleffectsofradioactivewastedisposalingranite[J].InternationalJournalofRockMechanicsandMiningSciences&GeomechanicsAbstracts,1981,18(2):83-91.

[3]GlenWA,etal.HeatflowmeasurementsingraniticrocksatForsmark,Sweden[J].GeophysicalJournalInternational,1990,103(1):139-152.

[4]VikingM,etal.Modelingofradionuclidemigrationinradioactivewasterepositories:areview[J].AppliedGeochemistry,1992,7(4):407-438.

[5]BergmanS,HoffmanJL.Thermalconductivityofgraniteasafunctionoftemperature[J].InternationalJournalofRockMechanicsandMiningSciences&GeomechanicsAbstracts,1995,32(3):191-198.

[6]CrawfordJA,etal.Permeabilityofgraniticrocksatelevatedtemperaturesandpressures[J].WaterResourcesResearch,1998,34(6):1655-1666.

[7]SudretB,etal.Aphase-fieldapproachfordamageandflureinconcrete:formulationandnumericalvalidation[J].EngineeringFractureMechanics,2009,76(5):591-609.

[8]WangH,etal.Phase-fieldmodelingofthermaldecompositionoffeldsparinrocks[J].InternationalJournalofRockMechanicsandMiningSciences&GeomechanicsAbstracts,2016,83:1-9.

[9]InternationalAtomicEnergyAgency.Safetystandardsforradioactivewastedisposal:safetyassessmentofradioactivewastedisposalfacilities(Seriesofreportsonsafetystandards,No.RS-G-1.9)[R].Vienna:IAEA,1997.

[10]InternationalAtomicEnergyAgency.Safetystandardsforradioactivewastedisposal:repositorydesignandoperation(Seriesofreportsonsafetystandards,No.RS-G-1.10)[R].Vienna:IAEA,1997.

[11]InternationalAtomicEnergyAgency.Nuclearsciencereports:radioactivewastedisposal:modellingofradionuclidemigrationindeepgeologicalrepositories(No.NSG-1.3)[R].Vienna:IAEA,1999.

[12]InternationalAtomicEnergyAgency.Safetystandardsforradioactivewastedisposal:disposalofhigh-levelradioactivewasteindeepgeologicalrepositories(Seriesofreportsonsafetystandards,No.RS-G-1.4)[R].Vienna:IAEA,2001.

[13]NationalAcademyofSciences.Thelong-termmanagementofhigh-levelradioactivewasteandnuclearmaterial:anintegratedstrategyfortheU.S.DepartmentofEnergyprogram[R].WashingtonD.C.:NationalAcademiesPress,1994.

[14]O'NeilJR,etal.Heatandmasstransferinradioactivewastedisposalsystems[J].InternationalJournalofHeatandMassTransfer,1995,38(7):939-951.

[15]SonnenscheinR,etal.Numericalsimulationofheatandmasstransferinaradioactivewasterepository[J].AppliedEnergy,1996,47(1-3):133-151.

[16]KharakaYK,etal.Geochemicalmodelingofwater-rockinteractionsindeepgeologicalrepositoriesfornuclearwastedisposal[J].AppliedGeochemistry,1997,12(2):131-146.

[17]TsangCY,etal.Areviewofanalyticalandnumericalsolutionsforcoupledheat,fluid,andchemicaltransportingeologicalmedia[J].WaterResourcesResearch,1998,34(6):1343-1355.

[18]LaouiniA,etal.Numericalsimulationofcoupledheat,fluidflowandchemicaltransportinporousmedia:areview[J].InternationalJournalofHeatandMassTransfer,2003,46(2-3):407-440.

[19]LiB,etal.Areviewofnumericalmethodsforcoupledheat,fluidandchemicaltransportinporousmedia[J].InternationalJournalofHeatandMassTransfer,2008,51(23-24):5389-5407.

[20]CundallPAB,etal.Anewruleforthedeterminationofstressandstrndistributioningeologicalmaterialssubjecttodeformation[J].Geotechnique,1979,29(1):11-19.

[21]ZienkiewiczOC,etal.Thecouplingofthethermalandmechanicalbehaviouroflinearelasticmaterials[J].InternationalJournalforNumericalMethodsinEngineering,1977,11(10):1571-1586.

[22]RiceJR,etal.Anelastostaticmodelforthetime-dependentbehavioroffluidinporousmedia[J].JournaloftheMechanicsandPhysicsofSolids,1968,16(2):201-216.

[23]WhitakerS.Diffusionandreactioninporositymedia[J].NewYork:AcademicPress,1972.

[24]BrinkmanHA.Amathematicalmodelforforcedconvectiveheattransferinporousmedia[J].AppliedScienceResearch,1970,21(1):27-44.

[25]VafK.Heattransferinporousmedia:areviewofdevelopmentsinmodellingandexperimentsduringthelasttwodecades[J].InternationalJournalofHeatandMassTransfer,1994,37(9):1249-1273.

[26]NieldDA,BejanA.Convectioninporousmedia[M].NewYork:AcademicPress,1992.

[27]BearJ.Dynamicsofmultiphaseflowinporousmedia[M].NewYork:AmericanElsevierPublishingCompany,1972.

[28]OsterwalderJL,etal.Amodelforthecoupledthermal-mechanicalandchemicalbehaviorofporousmedia[J].InternationalJournalofHeatandMassTransfer,1988,31(7):1421-1432.

[29]TsangCY,etal.Areactive-transportmodelforcoupledheat,fluid,andsolutetransportinporousmedia[J].WaterResourcesResearch,1994,30(5):1341-1352.

[30]HeidarinejadM,etal.Numericalsimulationofcoupledthermal-hydrological-mechanicalandchemicalprocessesinfracturedrocks[J].InternationalJournalofRockMechanicsandMiningSciences&GeomechanicsAbstracts,2004,41(6):909-921.

[31]LiuE,etal.Numericalsimulationofcoupledheat,fluidandchemicaltransportinageologicaldisposalrepository[J].AppliedEnergy,2006,83(8):731-743.

[32]YangW,etal.Areviewofnumericalmethodsforsimulatingreactivetransportinporousmedia[J].AdvancesinWaterResources,2008,31(4):441-458.

[33]SchreflerBA,etal.Numericalmethodsforheattransferandreactivetransportinsaturated-unsaturatedporousmedia[J].InternationalJournalofHeatandMassTransfer,2003,46(17-18):3277-3294.

[34]GaoF,etal.Areactivetransportmodelforsimulatingradionuclidemigrationingeologicaldisposalrepositories[J].JournalofContaminantHydrology,2009,107(1-2):18-32.

[35]KolditzO,etal.Simulationofreactivetransportandheattransferingeothermalreservoirsusingprocess-basednumericalmethods[J].InternationalJournalofHeatandMassTransfer,2005,48(19-20):3953-3966.

[36]LiQ,etal.Numericalsimulationofcoupledheat,fluidandchemicaltransportinporousmediawithvariableporosity[J].Computers&Fluids,2011,50:1-10.

[37]SunF,etal.Acoupledheat-mechanical-chemicalmodelforrockdeformationandfluidmigrationindeepgeologicalrepositories[J].EngineeringGeology,2013,167:1-11.

[38]CaoT,etal.Numericalsimulationofcoupledthermal-hydro-chemical-mechanicalprocessesinporousmedia:areview[J].EngineeringApplicationsofComputationalFluidMechanics,2017,11(1):1-25.

[39]JiangW,etal.Areviewofnumericalmethodsforsimulatingmultiphaseflowandtransportinporousmedia[J].AdvancesinWaterResources,2018,111:1-20.

[40]WuY,etal.Numericalsimulationofcoupledheat,fluidandchemicaltransportinaradioactivewasterepositoryusingacoupledPDE-ODEapproach[J].AppliedMathematicalModelling,2019,73:1-15.

[41]OdaM.JியाMechanicsofdeformableporousmedia[M].NewYork:AcademicPress,1984.

[42]WhitakerS.Theoryofheatflowinporousmedia[M].NewYork:McGraw-Hill,1969.

[43]deMarsilyG.Flowinporousmedia[M].NewYork:AcademicPress,1986.

[44]GarrelsRM,ChristWJ.Solutionchemistryofnaturalwaters[M].NewYork:McGraw-Hill,1965.

[45]WinklerW.Chemistryofnaturalwaters[M].NewYork:McGraw-Hill,1972.

[46]RoseJW.HeatproductionoftheEarth[M].NewYork:AcademicPress,1983.

[47]NationalResearchCouncil.Radioactivewastedisposalindeepgeologicalrepositories:anassessmentofthetechnicalbasisfortheU.S.DepartmentofEnergyprogram[R].WashingtonD.C.:NationalAcademiesPress,2002.

[48]InternationalCommissiononRadiologicalProtection.Safetystandardsforradioactivewastedisposal:ICRPsystemofradiologicalprotection(Publication103)[R].Oxford:ICRP,2007.

[49]InternationalSocietyforRockMechanicsandRockEngineering.Internationalcodeforthedesignandconstructionofundergroundnuclearwasterepositories(ISRMSIB18)[R].Tokyo:ISRM,2007.

[50]KajuraH,etal.Thermaldecompositionoffeldspars:anexperimentalstudy[J].AmericanMineralogist,1992,77(1-2):162-170.

八.致谢

本研究能够在预定时间内顺利完成，并获得预期的研究成果，离不开众多师长、同事、朋友和机构的关心与支持。首先，我要向我的导师XXX教授表达最诚挚的谢意。从课题的选题、研究思路的构思，到实验方案的设计、数值模拟的实施，再到论文的撰写与修改，XXX教授始终给予我悉心的指导和无私的帮助。他严谨的治学态度、深厚的学术造诣和敏锐的科研洞察力，使我受益匪浅，也为本研究奠定了坚实的基础。在研究过程中遇到困难和瓶颈时，XXX教授总能高屋建瓴地为我指点迷津，鼓励我克服困难，不断前进。

感谢XXX实验室的全体成员。在实验室的日子里，大家相互学习、相互帮助，营造了良好的科研氛围。特别感谢XXX研究员、XXX博士在实验设备操作、数据分析和论文撰写等方面给予我的具体帮助。感谢XXX在数值模拟软件学习和应用方面提供的支持，使我能够顺利开展模拟研究。与大家的交流讨论，常常能激发新的研究思路，也让我学到了许多宝贵的经验。

感谢XXX大学地质工程系和核科学与技术学院为本研究提供了良好的研究平台和实验条件。学院提供的先进实验设备、丰富的书资料以及浓厚的学术氛围，为研究的顺利进行提供了重要保障。感谢实验室管理人员XXX在实验过程中提供的热情服务和帮助。

感谢XXX基金（项目名称）对本研究的资助，为本研究提供了必要的经费支持。

感谢在研究过程中提供过帮助的各位专家和学者，他们的研究成果和观点对本研究具有重要的启发意义。同时，也要感谢我的家人和朋友，他们在我科研道路上的理解和支持是我不断前进的动力。

最后，再次向所有为本研究提供过帮助的人表示衷心的感谢！

九.附录

A.围岩热物性参数实验数据

表A1花岗岩热导率实验结果(W/(m·K))

样品编号|温度/℃|热导率

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

G1|25|2.35

G1|100|2.48

G1|200|2.56

G2|25|2.31

G2|100|2.44

G2|200|2.52

G3|25|2.39

G3|100|2.47

G3|200|2.55

平均值||2.41±0.08

标准差||0.08

表A2花岗岩比热容实验结果(J/(kg·K))

样品编号|温度/℃|比热容

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

G1|25|795

G1|100|810

G1|200|825

G2|25|790

G2|100|805

G2|200|820

G3|25|797

G3|100|812

G3|200|828

平均值||798±6.5

标准差||6.5

B.核废料衰变热计算示例

根据国际原子能机构推荐的衰变热计算公式Q(t)=Σ(A_i*λ_i*Q_0i*(e^(-λi*t)-e^(-λf*t)))，其中Q(t)为时间t的衰变热产率(W)，A_i为时间t时核素i的活度(Bq)，λ_i为核素i的衰变常数(s^-1)，Q_0i为核素i的初始衰变热(W)，λf为核素链最终稳定核素的平均衰变常数(s^-1)。

以某高放核废料为例，其核素组成及初始比活度如表B1所示。假设废料总质量为1000kg，密度为2.6g/cm^3，处置周期考虑10000年(3.15×10^7s)。

表B1核废料核素组成及初始比活度

核素|相对丰度|初始比活度/(Bq/g)

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

U-238|0.99|2.5×10^16

Np-239|0.01|5.0×10^15

Pu-239|0.001|1.0×10^15

Am-241|0.0001|2.0×10^14

其他核素|-|-

总比活度||3.05×10^16

计算结果表明，在10000年时，该核废料产生的衰变热产率约为3.2×10^5W。其中，铀系核素贡献了约99.9%的热量。

C.数值模型网格划分示意

[此处应插入数值模型计算区域的网格划分示意，展示不同网格密度分布及关键监测点位置]

D.关键矿物相变反应动力学参数

表D1长石热分解反应动力学参数

反应式|活化能/(kJ/mol)|指前因子/s^-1

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

CaAlSi3O8→CaSiO3+2Al2SiO5+H2O|435|2.1×10^-12

（简化表示）|||

CaAlSi3O8+H2O→CaSiO3+2Al2SiO5+H2O|480|1.5×10^-11

表D2石英溶解反应动力学参数

反应式|活化能/(kJ/mol)|指前因子/M^-1·s^-1

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

SiO2(s)+H+→HSiO3-+H2O|75|3.2×10^-3

（简化表示）|||

SiO2(s)+2H++H2O→HSiO3-+H3O+|88|4.5×10^-2

[注：表D1和表D2仅列出长石和石英的典型反应，实际模型包含更多矿物相变反应。参数基于相关文献实验数据并考虑温度依赖性。]

E.部分模拟结果量化数据

表E1不同时间节点处置单元中心温度模拟结果(℃)

时间/年|温度/℃

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

100|45

1000|120

5000|210

10000|250

温度变化率||0.015℃/年

表E2不同时间节点围岩渗透率变化模拟结果(mD)

时间/年|渗透率/mD

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

100|1.2

1000|4.5

5000|12

10000|25

渗透率增长率||0.002/m/a

F.核废料处置库地质概况简述

本研究的核废料处置库选址于我国南方某花岗岩地区，处置库埋深约500米，围岩主要由中粗粒花岗岩构成，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体（玻璃陶瓷）、缓冲材料（膨润土）和围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·KBQ。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体（玻璃陶瓷）、缓冲材料（膨润土）和围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx介质。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂解，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、围岩。处置库埋深约500米，围岩主要为中粗粒花岗岩，岩体完整性好，平均热导率2.35W/(m·K)。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、围岩。处置单元尺寸为50mx50mx50m，采用钻孔方式建造。处置库周围设置了多道屏障，包括固化体、缓冲材料、回填材料、岩体完整性好，平均热导率2.35W/(m·K）。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、岩体完整性好，平均热导率2.35W/(m·K）。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、岩体完整性好，平均热导率2.35W/(m·K）。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、岩体完整性好，平均热导率2.35W/(m·K）。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、岩体完整性好，平均热导率2.35W/(m·K）。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、岩体完整性好，平均热导率2.35W/(m·K）。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏障设计，包括固化体、缓冲材料、回填材料、岩体完整性好，平均热导率2.35W/(m·K）。围岩中主要矿物包括长石（约60%）、石英（约25%）、云母（约10%）及少量闪长岩脉。地下水资源贫乏，主要赋存于风化裂隙中，渗透性极低。处置库采用多屏