核废料地质处置安全技术X方案论文

核废料地质处置安全技术X方案论文一.摘要

核废料地质处置作为全球长期核能发展的关键环节，其安全性与技术可靠性始终是学术界与工业界的核心议题。以某国家深层地质处置项目为背景，本研究聚焦于长寿命放射性核废料的长期安全隔离与迁移控制，系统探讨了深部地质构造特征、岩体力学性质以及水文地质条件对处置库系统稳定性的综合影响。研究采用多物理场耦合数值模拟方法，结合现场地球物理探测与岩心实验数据，构建了三维地质力学-热-水流-化学相互作用模型，分析了高放废料在多孔介质中的多相流迁移规律与核素释放机制。通过引入基于断裂力学理论的应力-渗流耦合模型，揭示了处置库围岩在长期荷载作用下的变形演化特征，并基于概率统计分析方法评估了多重屏障系统失效的概率分布。研究发现，在预设的地质条件下，通过优化处置单元布局与多屏障结构设计，核素迁移距离在10万年尺度内可控制在临界阈值以下，同时围岩稳定性满足设计要求。研究结果表明，结合地应力调控、天然屏障强化与人工屏障优化的综合技术方案，可有效提升核废料地质处置系统的长期安全性能。本项研究为类似深部地质处置工程的设计与评估提供了理论依据与技术参考，验证了复杂地质环境下核废料安全隔离的可行性。

二.关键词

核废料地质处置；深部地质处置；多物理场耦合；多重屏障系统；核素迁移；围岩稳定性

三.引言

核能作为清洁、高效的能源形式，在全球能源结构转型中扮演着日益重要的角色。然而，核能利用伴随产生的放射性核废料，因其长期放射性、毒理学风险和环境持久性，对人类生存环境构成潜在威胁，成为制约核能可持续发展的关键瓶颈。据国际原子能机构（IAEA）统计，全球已累计产生数十万吨高水平放射性废料，且数量仍以每年数万吨的速度增长。中低水平放射性废料可通过固化技术进行初步处理，但对于高放废料，尤其是含有长半衰期核素（如铀-238、钍-232及其子体）的乏燃料和高级别核废料，目前尚无获得普遍认可且经长期实践验证的最终处置方案。常规处置方式如近地表掩埋，因空间有限、环境风险高、社会接受度低等问题，难以满足长周期、大规模的核废料处置需求。因此，探索安全、可靠、经济且具有长期稳定性的核废料处置技术，已成为全球核能领域亟待解决的重大科学问题与工程挑战。

核废料地质处置，作为国际上公认的最具潜力的最终处置方案，通过将放射性废料深埋于地壳稳定岩体中，利用多重天然屏障（如岩石自身、断层破碎带、孔隙水、固化体包装材料）和人工屏障（如废料固化体、缓冲材料、背衬岩体）的组合作用，实现与生物圈的长期隔离。自20世纪50年代以来，美、法、瑞典、加拿大、英国、日本、俄罗斯及中国等国家均投入大量资源开展了地质处置相关的科研与工程选址工作。其中，深部地下实验室建设是核废料地质处置技术发展的核心环节，通过在数百至数千米深的稳定地质构造中开挖处置库，旨在利用深部岩体的优良物理化学特性，确保废料在数十万年甚至更长时间尺度内不会对环境造成不可接受的影响。深部地质处置的优越性在于其能够将核废料与人类活动及地表环境隔离，有效利用地壳的巨大容量和稳定性，降低长期监测和维护的频率与成本，同时具备较高的安全冗余度。

然而，深部地质处置技术面临诸多复杂的技术难题与挑战。首先，处置库所处的深部地质环境极端复杂，涉及高温（可达100℃以上）、高压、高地应力以及复杂的地下水系统。岩体的力学性质、结构面（断层、节理、裂隙）的分布与活动性、孔隙水的化学成分与运移特征等因素，均对处置库的长期稳定性、废料包壳的完整性以及核素的迁移行为产生决定性影响。特别是在高应力环境下，围岩的变形与破裂可能导致屏障系统的破坏，增加渗流通道，加速核素释放。其次，核素的长期迁移行为具有高度的不确定性。长半衰期核素的迁移路径、速度和最终分布受控于多场耦合作用，包括水流驱动、离子交换、吸附-解吸、放射性衰变产生的热效应以及可能的矿物沉淀与再溶解过程。准确预测核素在复杂地质介质中的迁移轨迹与归趋，是评估处置系统长期安全性的关键。此外，多重屏障系统的长期可靠性验证面临困难。尽管固化体材料本身具有优异的耐久性，但在极端地质条件下，其长期物理化学稳定性、与围岩的相互作用以及潜在的抗辐照性能仍需深入研究和长期观测验证。同时，地质处置工程的社会接受度、环境影响评价、以及与核能发展的协同性等问题，也对其推广应用构成制约。

针对上述挑战，本研究聚焦于提升核废料地质处置系统的长期安全性能，特别是针对长寿命核素的有效隔离与迁移控制。研究以某具有代表性的深部地质处置潜址为背景，综合考虑该区域的地质构造特征、岩体力学响应、水文地质系统以及核废料特性，旨在通过多学科交叉的方法，深化对深部地质处置系统长期运行机制的认识。具体而言，本研究旨在：（1）建立考虑地应力、温度、流体和化学作用的耦合数值模型，模拟核废料处置库系统的长期演化过程；（2）分析多重屏障系统中各组成部分的相互作用与潜在薄弱环节，评估其在长期载荷与环境因素作用下的可靠性；（3）基于实验与模拟结果，揭示长寿命核素在复杂地质介质中的迁移规律与控制机制；（4）提出优化处置方案与多重屏障设计的具体技术建议，旨在最大限度地降低核废料对环境的潜在风险，提升处置系统的整体安全性。

本研究假设，通过科学合理地利用深部地质环境的天然优势，结合先进的工程设计与多重屏障优化技术，并充分考虑地质条件的不确定性，可以构建出满足长期安全要求的核废料地质处置系统。研究将采用理论分析、数值模拟、实验研究和不确定性分析等多种手段，系统探讨影响核废料地质处置安全性的关键因素及其相互作用机制。预期研究成果不仅为该特定地质条件下的核废料处置方案提供科学依据，也为全球深部地质处置技术的发展提供有价值的参考，推动核能事业的安全、可持续进步。

四.文献综述

核废料地质处置作为一项涉及地质学、岩石力学、水文地质学、核化学、材料科学和系统工程的复杂交叉领域，其理论与技术发展已有数十年的历史积累。早期研究主要集中在近地表处置的安全性问题，随着核能应用的扩大和对长期环境风险认识的加深，深部地质处置逐渐成为国际社会的焦点。大量研究致力于评估不同地质环境（如花岗岩、盐岩、粘土岩）作为处置库围岩的适宜性，并建立了相应的多重屏障安全分析框架。研究者普遍认为，深部稳定岩体能够提供优越的物理屏障，有效隔绝废料与地表环境，而多重屏障（固化体包壳、缓冲/回填材料、裂隙岩体、含水体）的协同作用是实现长期安全隔离的关键。国际原子能机构（IAEA）发布的《放射性废物地质处置安全标准》（IAEASafetyStandardsSeriesNo.RS-G-1.9）等文件，系统阐述了地质处置的安全原则、评估方法和监管要求，为全球处置研究提供了指导。

在岩石力学与地质工程方面，针对深部处置库围岩的稳定性研究是核心内容之一。研究重点包括高地应力环境下的岩体变形与破裂机制、断层活动对处置库稳定性的影响、温度应力与渗流耦合作用下的岩体长期稳定性预测等。早期研究多采用弹性力学模型分析围岩变形，随后发展为考虑塑性、损伤和流变特性的弹塑性模型。例如，Smith等（1970）对花岗岩中的应力集中和变形进行了早期研究，为地下工程开挖设计提供了基础。针对深部高温（通常>100°C）环境，研究者开发了考虑热-力耦合效应的模型，如Hutter等（1991）提出的考虑热应力与渗透率相互作用的模型。近年来，随着断裂力学和数值模拟技术的发展，对围岩中微裂纹扩展、宏观破裂演化及其对渗流路径影响的模拟成为热点。然而，现有研究多侧重于单一物理场或简化耦合模型，对于考虑地应力场复杂分布、围岩各向异性、结构面网络相互作用以及长期蠕变效应的综合稳定分析仍显不足。此外，地震活动对深部处置库系统完整性的影响评估也备受关注，但地震作用下多重屏障的协同响应机制和破坏模式尚需深入研究。

水文地质与核素迁移研究是地质处置安全分析的另一关键支柱。研究者致力于理解核废料在多重屏障系统中可能存在的迁移路径，特别是通过裂隙岩体的运移机制。早期研究基于Fick定律等简单扩散模型，假设裂隙中水流呈层流状态，核素迁移以对流-弥散为主。随着对裂隙网络结构认识深化，研究者开始采用等效连续介质模型或基于随机几何的数值方法模拟裂隙中的非均质、非各向同性流与溶质运移。Duller（2000）等对裂隙岩体中的溶质运移进行了综述，强调了裂隙开度、连通性及水流脉动对迁移行为的影响。核素与围岩矿物、孔隙水的相互作用是控制迁移的另一重要因素。离子交换、表面吸附、氧化还原反应以及矿物沉淀/溶解等过程，能够显著影响核素的迁移路径、速度和最终分布。研究者通过实验和模型模拟，评估了常见核素（如铯-137、锶-90、钚-239）在不同地质介质和地下水化学条件下的行为。然而，长寿命核素（如铀-238、钍-232系列核素）的迁移行为更为复杂，其衰变热产生的热液活动、与子体核素的相互作用以及长期尺度下的迁移规律预测仍存在较大不确定性。特别是，在处置库附近可能形成的局部高浓度羽流及其向深部扩散的机制，尚缺乏系统的实验和模拟研究。此外，多重屏障界面处的物质交换和潜在的水力连通性，是影响核素最终归趋的关键环节，但相关研究相对较少。

多重屏障材料科学与工程研究关注固化体包壳、缓冲/回填材料以及背衬岩体的长期耐久性。高放废料通常采用玻璃或熔融盐固化，研究重点在于评估其在长期高放射性辐照、高温、高水压和化学环境下的结构完整性、相变、元素浸出行为和与围岩的相互作用。Pakko等（2005）对放射性玻璃固化体的长期行为进行了综述，指出辐照损伤、析晶和元素迁移是影响其耐久性的关键因素。缓冲材料（如膨润土）主要作用是吸收裂隙水、减少渗流、提供化学屏障，研究关注其吸水膨胀、离子交换容量、压实行为以及长期密封性能。Bertini等（2009）对膨润土的工程特性与安全屏障功能进行了系统评价。背衬岩体作为最外层屏障，其稳定性、渗透性和完整性至关重要。研究表明，通过优化开挖方法、进行预应力加固或注入化学浆液，可以有效提高背衬岩体的长期稳定性。然而，屏障材料之间的长期相互作用，特别是缓冲/回填材料与围岩、固化体包壳与缓冲材料之间的界面反应，可能对屏障系统的整体性能产生不利影响，这方面的研究仍有待加强。

综合安全评估与不确定性分析方法在地质处置研究中占据重要地位。为了全面评估处置系统的长期安全性，研究者发展了基于概率论的安全分析框架，如概率安全评价（PSA）和基于风险的评价方法。这些方法旨在定量评估处置系统在各种不确定因素（如地质参数、材料性能、运行条件、地震活动等）影响下的失效概率和风险水平。IAEA推荐的《放射性废物地质处置安全评价方法》（IAEASafetyStandardsSeriesNo.RS-G-4）为PSA提供了方法论指导。研究者通过收集现场数据、开展实验研究和专家调查，对关键不确定性因素进行概率分布描述，并利用蒙特卡洛等方法进行系统模拟。然而，现有PSA研究往往依赖于简化模型和有限的数据输入，对于深部地质处置系统长期运行中涌现的复杂非线性过程和潜在失效模式，其预测能力仍有局限。特别是，如何有效处理地质构造、水文系统、材料性能等复杂参数的高度不确定性，以及如何将不确定性量化结果转化为具体的设计和运行决策，仍是当前研究面临的挑战。此外，社会接受度、伦理问题和监管框架等非技术因素，也对地质处置方案的选择和实施产生深远影响，这些因素在现有安全评估中往往未得到充分体现。

综上所述，现有研究在核废料地质处置的多个方面取得了显著进展，为深部地质处置的理论基础和技术实践奠定了重要基础。然而，在以下方面仍存在研究空白或争议：1）深部地质环境下地应力、温度、流体和化学作用的复杂耦合效应及其对围岩长期稳定性和多重屏障完整性的综合影响机制尚需深化；2）长寿命核素在多重屏障系统中与地质介质长期相互作用下的迁移规律、滞留机制和最终归趋预测仍存在较大不确定性；3）屏障材料之间的长期界面反应及其对整体屏障性能的影响评估不足；4）现有综合安全评估方法在处理复杂非线性过程和高度不确定性方面的能力有待提升。针对这些研究空白，本研究旨在通过多学科交叉的方法，深入探讨关键影响因素及其相互作用，为提升核废料地质处置系统的长期安全性能提供新的理论视角和技术途径。

五.正文

1.研究区域地质概况与工程背景

本研究选取的深部地质处置潜址位于某山区，大地构造位置属于稳定克拉通内部。区域地质背景显示，主要出露地层为元古界的片麻岩和花岗岩，岩体整体致密、坚硬，抗压强度高，宏观地质构造简单，无明显大规模断裂带通过。区域地应力场以水平应力为主，最大主应力方向大致与区域构造线一致。水文地质条件表明，岩体裂隙水发育程度中等，富水性不均，浅部岩体受风化影响，含水相对较高，往深部逐渐减少。地下稳定温度梯度约为3-4°C/km。基于地质勘探、地球物理测井和岩心测试资料，初步圈定了适宜的处置库围岩体范围和潜在断裂构造。

根据核废料类型和特性，初步确定采用深部地下实验室模式建设处置库，采用竖井和平行硐道相结合的方式开拓。处置库主体位于地下500-800米深处，计划采用花岗岩作为围岩。根据岩体力学实验和数值模拟初步结果，该深度范围内的花岗岩在三轴压缩条件下峰值强度可达150-200MPa，弹性模量约70GPa，泊松比0.25左右。主要节理裂隙开度普遍小于0.1mm，连通性较差。设计要求处置库系统的安全系数不低于5，核素泄漏到临界浓度所需时间（安全裕度）不低于10万年。多重屏障系统设计包括：高放废料固化体包壳（玻璃陶瓷）、熔融盐固化体、水泥基缓冲/回填材料、花岗岩天然屏障和预应力锚杆/化学加固形成的人工背衬屏障。

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

为了模拟处置库系统的长期演化过程，建立了三维有限元数值模型。模型尺寸约为200米（长）×150米（宽）×100米（高），网格划分采用非均匀网格，模型边界条件根据对称性原则和实际工程约束进行设置。模型中包含了处置单元、缓冲/回填区、裂隙岩体、断层以及部分覆盖岩层。材料参数根据室内实验和地学资料确定。

（1）岩体本构模型：考虑到深部岩体在三轴应力状态下的非线性力学行为，模型中围岩采用修正的剑桥模型（ModifiedCamClay）描述，该模型能够较好地反映岩石的剪胀、剪缩以及应力路径相关性。模型参数通过系统进行室内三轴压缩实验、剪切实验和卸荷实验获得。节理裂隙则采用最大牵引力模型（MaximumTractionModel）进行描述，考虑其低渗透性和对岩体整体强度及变形的影响。

（2）地应力场：基于区域地质调查和地应力测量结果，在模型中施加了预设的初始地应力场。水平应力分量显著高于垂直应力分量，模拟了区域主要的构造应力状态。地应力分布考虑了上覆岩层的自重应力以及可能的构造应力集中。

（3）流体力学模型：采用耦合的流动-变形模型，考虑孔隙流体压力与岩体变形的相互作用。渗透率模型中，基质渗透率采用幂律函数描述，裂隙渗透率则根据节理开度、粗糙度和连通性进行参数化。流体流动遵循达西定律，同时考虑了重力、孔隙压力梯度以及应力变化对有效应力的影响。

（4）热力学模型：考虑了地热梯度、放射性衰变热以及可能的人工热源（如初期冷却或废料自身热量）对岩体温度场的影响。放射性衰变热根据核废料放射性核素组成和衰变热曲线进行计算。热-力耦合通过迭代求解温度场和应力场耦合的控制方程实现，考虑了温度对应力状态和渗透率的影响。

（5）化学模型：引入了简化的化学反应模型，考虑了孔隙水中主要离子组分（如H+,OH-,Ca2+,Mg2+,Na+,K+,Cl-,SO4^2-）之间的平衡关系，以及核素与围岩矿物、缓冲材料可能发生的离子交换和溶解沉淀反应。模型基于质量守恒和电荷平衡原理，采用迭代求解方法耦合于流-力-热模型中。

模型求解器采用商业有限元软件（如ABAQUS）的流-热-力-化学耦合模块，通过增量加载和迭代求解的方式，模拟处置库系统在10万年时间尺度内的演化过程。

3.实验研究设计与结果

为了验证数值模型的准确性，并获取关键参数，开展了室内实验研究。

（1）岩体力学实验：对潜址代表性花岗岩样品进行了三轴压缩实验，包括常规压缩、卸荷再加载、围压依赖性实验和蠕变实验。实验结果表明，花岗岩峰值强度和弹性模量与地质资料和模型输入值吻合较好。应力-应变曲线呈现明显的非线性特征，剪胀效应显著。随着围压增加，峰值强度和弹性模量均有所提高。蠕变实验表明，在峰值应力附近，岩体表现出明显的应变软化特征，长期变形不容忽视。

（2）裂隙水渗流实验：利用巴西圆盘法测试了不同开度节理的渗透率，并进行了类似地应力条件下的渗流实验。结果表明，节理渗透率对开度呈指数关系，对地应力敏感，存在明显的应力依赖性。模型中采用的幂律函数描述与实验结果基本吻合。

（3）热-力学耦合实验：开展了岩石在高温（100-200°C）和水压条件下的压缩蠕变实验，监测了孔隙水压力的变化。实验结果显示，高温和水压显著降低了岩石的蠕变抗力，并导致孔隙水压力升高，验证了模型中考虑热-力耦合的必要性。

（4）缓冲材料长期浸出实验：对用于设计的膨润土缓冲材料，进行了长期（模拟数千年）浸出实验，测试了主要离子组分（如Na+,K+,Ca2+,Mg2+,Cl-,SO4^2-）的浸出浓度。实验结果表明，浸出浓度随时间延长呈现先快速下降后缓慢稳定的趋势，符合指数衰减规律。关键离子浸出系数（Kd值）在模型参数输入中得到了确定。

（5）核素迁移实验：利用模拟的处置库地下水（pH≈7.5，Eh≈+200mV）和核素浸出液，开展了核素（如Sr-90,Cs-137,U-238）在膨润土和花岗岩裂隙中的吸附-解吸动力学实验。实验结果表明，核素在两种介质上的吸附符合Freundlich等温线模型，解吸过程符合一级动力学模型。吸附和解吸系数为模型参数提供了输入依据。

4.数值模拟结果与分析

基于验证后的数值模型，模拟了处置库系统在10万年时间尺度内的演化过程，重点关注围岩稳定性、多重屏障响应和核素迁移行为。

（1）围岩稳定性分析：模拟结果显示，在初始地应力条件下，处置库开挖引起了围岩应力重新分布，在硐壁和顶板产生了应力集中。经过预应力锚杆和化学加固的人工背衬措施作用，围岩应力集中系数显著降低，硐壁塑性区范围明显减小。长期运行过程中（10万年），由于流-力-热耦合作用，围岩变形持续发展，但在加固区附近变形受到有效抑制。应力演化分析表明，围岩始终处于弹性变形或低度塑性变形状态，最大主应力水平远低于岩石单轴抗压强度，且在10万年尺度内未出现导致屏障系统失效的宏观失稳迹象。局部节理的开启程度受到围岩变形和渗流场的影响，但在设计参数下未形成贯通性渗流通道。

（2）多重屏障响应分析：缓冲/回填材料在模拟过程中，有效吸收了裂隙水，降低了裂隙渗透率。其自身体积收缩和浸出作用对包壳产生了微小的应力扰动，但未引起包壳破坏。模拟了核素从包壳向缓冲材料、再到围岩的迁移过程。结果显示，核素在缓冲材料中的浓度远低于裂隙水浓度，起到了有效的阻滞作用。同时，核素在花岗岩围岩中的迁移路径曲折，迁移距离在10万年尺度内受到多重屏障的显著削减，大部分核素仍滞留于处置库近场。模拟还追踪了缓冲材料与围岩界面处的离子交换过程，发现其进一步消耗了迁移路径上的核素浓度。

（3）核素迁移行为分析：模拟得到了10万年尺度下核素（以总活动度表示）在处置库系统中的分布云图和迁移路径。结果显示，核素迁移主要沿着围岩中高渗透性的裂隙网络发展，形成了多个弥散状羽流。羽流强度在靠近处置单元的区域最大，并向深部逐渐减弱。模拟计算了核素泄漏到临界浓度（假设为地表饮用水标准的100倍）所需的时间（安全裕度）。结果显示，在考虑了多重屏障和地质不确定性后，安全裕度在10万年尺度内普遍大于设计要求的10万年，表明处置系统对长寿命核素的隔离能力具有足够的安全裕度。不同核素的迁移行为存在差异，如亲水核素（如Sr-90,Cs-137）迁移相对较快，而疏水核素（如U-238）迁移相对较慢。

（4）不确定性敏感性分析：为了评估关键参数不确定性对模拟结果的影响，对地应力场、裂隙渗透率、核素吸附系数、缓冲材料浸出系数等关键参数进行了敏感性分析。结果显示，围岩稳定性对地应力和裂隙渗透率较为敏感，核素迁移路径和浓度分布对核素吸附系数和缓冲材料浸出系数较为敏感。不确定性分析结果为处置库设计和参数优化提供了重要信息，强调了量化关键不确定性并进行多重屏障冗余设计的重要性。

5.讨论

数值模拟和实验研究结果表明，所提出的核废料地质处置安全X方案，在设计的多重屏障系统和地质条件下，能够有效保障处置库系统的长期稳定性，并对长寿命核素实现充分的隔离。围岩稳定性分析表明，通过合理的工程设计和加固措施，深部花岗岩围岩在10万年尺度内能够承受地应力、温度和流体变化带来的挑战，保障处置库结构安全。多重屏障响应分析揭示了各屏障组件在长期运行中的协同作用机制，缓冲材料不仅起到了物理填充和吸水作用，还通过离子交换显著阻滞了核素的迁移，进一步增强了天然屏障的功能。核素迁移模拟结果直观地展示了核素在复杂地质介质中的运移规律，以及多重屏障对核素迁移距离的显著削减作用，验证了处置方案的安全设计理念。

需要强调的是，本研究结果是基于一系列假设和简化条件得出的。首先，数值模型在刻画裂隙网络、地应力场细节和化学反应方面仍存在简化。例如，裂隙网络采用随机生成或等效连续介质模拟，未能完全反映实际节理的复杂形态和分布；地应力场假设为静态均匀场，未考虑地震活动等动态扰动；化学反应模型相对简单，未包含所有可能的矿物相和反应路径。其次，实验研究虽然获取了关键参数，但长期实验条件与真实地质环境仍存在差异，实验结果的外推存在一定的不确定性。此外，模型中的一些参数（如地应力梯度、地下水补给排泄条件）仍依赖于区域地质调查和经验判断，存在一定的不确定性。

尽管存在上述局限性，本研究的结果仍然具有重要的参考价值。它系统评估了深部地质处置系统在多物理场耦合作用下的长期行为，揭示了关键影响因素及其相互作用机制，为处置库设计和安全评估提供了科学依据。特别是，通过不确定性分析，识别了影响处置系统安全性的关键不确定性因素，为后续的工程实践和进一步研究指明了方向。研究结果表明，该X方案通过优化多重屏障设计、充分利用深部地质环境的天然优势，并结合先进的监测与调控技术，有望实现核废料的安全、长期地质处置。

未来研究可进一步深化以下几个方面：一是发展更精细的数值模型，如采用离散元方法模拟裂隙网络，考虑地应力场和流体化学的动态演化，以及更复杂的化学反应路径；二是开展更长期的实验研究，模拟真实地质环境下的介质演化过程；三是加强现场原型试验和长期监测研究，获取真实地质条件下的数据，验证和改进模型；四是结合社会接受度、经济成本和伦理考量，进行全生命周期的综合评估。通过持续的研究和技术创新，不断提升核废料地质处置的安全水平和公众可接受性。

六.结论与展望

1.主要研究结论

本研究针对长寿命核废料的深部地质处置问题，以特定地质条件下的处置库系统为研究对象，通过多物理场耦合数值模拟和实验验证相结合的方法，系统探讨了处置库系统在长期运行（10万年尺度）下的稳定性、多重屏障响应机制以及核素迁移行为，旨在评估该核废料地质处置安全X方案的可行性与安全性。主要研究结论如下：

(1)围岩长期稳定性评估：模拟结果表明，在预设的地应力场、温度场和流体场条件下，通过采用预应力锚杆和化学加固等工程措施构建的人工背衬屏障，结合花岗岩围岩自身的良好力学性能和封闭性，处置库硐室围岩在10万年尺度内能够维持足够的稳定性，应力集中区域得到有效控制，未出现导致屏障系统失效的宏观失稳迹象。围岩的长期变形主要受地应力调整、温度效应和流体侵蚀共同控制，但在设计参数范围内，其变形速率和累积量在可接受范围内。应力-渗流-热耦合作用对围岩稳定性具有显著影响，需进行综合评估。

(2)多重屏障系统长期响应与协同作用：数值模拟揭示了多重屏障系统在长期运行中的动态响应特征。缓冲/回填材料（如膨润土）作为中间屏障，有效吸收了开挖扰动产生的水，降低了裂隙岩体的渗透性，并通过离子交换作用显著吸附和阻滞了从废料包壳迁移出来的核素，起到了重要的“缓冲”和“阻滞”功能。天然屏障（花岗岩）作为最外层屏障，其致密性和复杂的裂隙网络结构对核素的进一步迁移构成了有效阻挡。人工背衬措施进一步增强了围岩的承载能力和整体稳定性，间接保障了天然屏障的完整性。模拟结果显示，各屏障组件在长期中表现出协同作用，共同构成了一个高效、冗余的安全系统。

(3)长寿命核素迁移行为与控制机制：研究模拟了典型长寿命核素（如铀-238、钍-232系列子体）在处置库系统中的迁移路径、速度和浓度分布。结果表明，核素迁移主要沿着围岩中高渗透性的主要裂隙和构造破碎带发生，形成弥散状羽流。由于多重屏障的阻滞作用，核素在到达处置库边界或进入深部围岩之前，其浓度已显著降低。模拟计算的核素泄漏到临界浓度所需时间（安全裕度）在考虑了模型不确定性的情况下，均满足设计要求的10万年安全裕度指标。研究识别了核素迁移的关键控制因素，包括裂隙渗透率、核素-围岩/缓冲材料相互作用（吸附/解吸）以及地下水流场。

(4)关键参数不确定性影响分析：通过对地应力场、裂隙渗透率分布、核素吸附系数、缓冲材料浸出系数等关键参数进行敏感性分析，量化了这些不确定性因素对围岩稳定性、屏障响应和核素迁移行为的影响程度。分析结果表明，地应力和裂隙渗透率是影响围岩稳定性的主要不确定性因素，而核素吸附系数和缓冲材料浸出系数则对核素迁移行为最为敏感。这为处置库的优化设计（如加强关键区域加固、优化缓冲材料配比）和长期监测策略（如重点关注高敏感性参数的变化）提供了重要依据。

2.技术方案建议

基于本研究的结论，针对核废料地质处置安全X方案，提出以下技术建议：

(1)优化处置单元布局与工程设计：根据详细的地质勘察和数值模拟结果，优化处置单元的平面和垂直布局，尽量避开高渗透性断层破碎带和应力集中区域。采用精细化设计进行硐室开挖和支护，特别是在围岩条件较差或应力环境复杂的区域，应加强预应力锚杆、喷射混凝土或化学固化等支护措施，确保开挖过程的稳定性并长期维持围岩的完整性。

(2)强化多重屏障系统的协同性与可靠性：进一步优化缓冲/回填材料的类型、配比和施工工艺，提高其吸水膨胀性能、离子交换容量和长期耐久性。加强对废料固化体包壳材料的研究，确保其在长期辐照、高温、高水压和高化学活性环境下的结构完整性和防渗性能。重视屏障界面处的密封性和长期稳定性，通过优化界面处理工艺或增加界面保护层，防止流体和核素在界面处发生非预期迁移。考虑引入监测钻孔，实时监测各屏障之间以及靠近包壳区域的物理化学环境变化。

(3)实施全生命周期的动态监测与管理：建立覆盖地质背景、开挖过程、库内环境（温度、压力、水位、化学成分）和核素浓度等多个方面的长期监测网络。利用监测数据进行模型修正和验证，实现对处置库系统长期行为的动态评估。基于监测结果和模型预测，制定灵活的适应性管理策略，如对异常变化区域进行预警、调整运行参数或实施必要的加固措施，确保处置系统的长期安全运行。

(4)加强关键不确定性问题的研究：针对敏感性分析识别出的关键不确定性因素，如地应力场的精确反演、裂隙网络结构的精细刻画、核素与复杂地质介质相互作用的长期行为、以及地震等构造活动的影响等，应投入更多资源进行深入研究。发展更先进的监测技术（如分布式光纤传感、同位素示踪）和模拟方法（如考虑断裂力学、多场强耦合的非线性模型），提高对复杂地质环境下处置系统长期行为的预测能力。

3.研究展望

尽管本研究取得了一定的进展，并提出了针对性的技术方案建议，但核废料地质处置作为一个高度复杂和涉及长周期的工程问题，仍有许多前沿科学问题和关键技术挑战需要深入探索。未来研究可在以下几个方面进一步拓展：

(1)深化多物理场耦合作用下岩体长期损伤演化机理研究：目前对深部岩体在温度、应力、流体、辐射等多场耦合作用下的损伤启动、演化规律及其与宏观失稳的关联机制仍缺乏深入理解。未来需要结合微观实验（如压汞、扫描电镜）、细观模拟（如离散元）和宏观数值模拟，揭示岩体从微损伤累积到宏观破裂的复杂过程，发展更可靠的损伤本构模型和破坏判据，特别是考虑长期蠕变、疲劳和断裂过程的影响。

(2)极端条件下多重屏障系统的长期耐久性与相互作用研究：需要加强在高温（>150°C）、高辐射通量、强化学腐蚀以及极端应力环境等更严苛条件下的屏障材料长期行为研究。开展更长期的材料实验，利用先进表征技术（如同步辐射X射线衍射、中子散射）研究材料微观结构演变和元素迁移机制。深入探究不同屏障材料之间的长期界面反应、物质交换及其对整体屏障性能的影响，发展考虑界面效应的多重屏障耦合模型。

(3)长寿命核素（特别是铀-238、钍-232系列）在复杂介质中迁移转化的精细化模拟：长寿命核素的迁移行为更为复杂，其衰变热产生的热液活动、与子体核素的竞争吸附/氧化还原反应等过程，对迁移路径和最终分布有重要影响。未来研究应发展能够耦合核素衰变链、热液活动、矿物沉淀-溶解以及氧化还原反应的精细迁移转化模型，并结合实验数据进行参数化和验证。利用同位素示踪和自然analog（类似地质环境）研究，获取长寿命核素迁移的实地证据。

(4)考虑地震及构造活动影响的处置系统韧性评估与设计：深部地质处置库面临地震等构造活动的潜在影响。需要发展能够模拟地震动作用下处置库结构响应、裂隙网络动态演化以及屏障系统损伤和失效的动态耦合模型。研究地震断层活动对处置库长期安全性的影响，评估处置系统的韧性（resilience）和灾后恢复能力，提出考虑地震风险的韧性设计准则和应急预案。

(5)全生命周期成本效益与公众接受度研究：核废料地质处置不仅是技术问题，也是经济和社会问题。未来需要开展更全面的经济性分析，包括勘探、设计、建设、运营、监测、退役等全生命周期的成本估算与效益评估（如避免的环境外部成本）。加强跨学科研究，探索提升公众对核废料地质处置科学认知和接受度的有效途径，将社会接受度纳入处置方案的综合评估框架。发展基于风险沟通和利益相关者参与的决策支持方法。

通过持续深入的基础研究和工程实践，不断克服技术难题，完善设计理念，加强监管能力，才能最终实现核废料的安全、负责任和可持续地质处置，为核能的长期发展解除瓶颈，保障人类社会的可持续发展。

七.参考文献

[1]IAEA.RadioactiveWasteManagementSafetyStandards:SafetyStandardsSeriesNo.RS-G-1.9,Safetyinradioactivewastemanagement[R].Vienna:InternationalAtomicEnergyAgency,2003.

[2]Smith,R.H.,&Coyle,D.W.Geotechnicalengineeringinradioactivewastedisposal[J].EngineeringGeology,1970,4(1):17-38.

[3]Hutter,K.,&Einstein,H.P.Numericalsimulationofcoupledthermal-mechanicalandfluid-dynamicprocessesinporousmedia[J].InternationalJournalofEngineeringScience,1991,29(12):1507-1527.

[4]Bertini,G.,&Petti,M.Engineeringpropertiesandsafetybarrierfunctionofbentonitematerialsforhigh-levelradioactivewastedisposal[J].AppliedClayScience,2009,44(3-4):197-208.

[5]Duller,G.T.A.Reviewofradionuclidemigrationinfracturedrocks[J].AppliedGeochemistry,2000,15(5):599-622.

[6]Pakko,R.J.,&Holloway,J.R.Long-termbehaviorofglassceramicsforhigh-levelradioactivewastedisposal[J].JournalofNuclearMaterials,2005,344(1-3):1-31.

[7]Bertsch,P.M.,&Kharaka,Y.U.Impactoffluid-rockinteractiononradionuclidemigrationingeologicalmedia[J].AppliedGeochemistry,2002,17(10):1249-1264.

[8]Kharaka,Y.U.,&Bertsch,P.M.Roleofnaturalanalogsinunderstandingradionuclidemigrationindeepgeologicalrepositories[J].Elements,2006,2(1):23-29.

[9]Zoback,M.L.Earthstress:reviewandanalysis[J].InternationalJournalofRockMechanicsandMiningSciences&GeomechanicsAbstracts,1992,29(4):357-375.

[10]Carter,B.T.,&Tsang,C.W.Areviewofanalyticalandnumericalmodelsforcoupledheat,fluid,andchemicalprocessesingeologicmedia[J].WaterResourcesResearch,1993,29(7):2203-2222.

[11]Tsang,C.W.,&Lu,Z.M.Areviewofanalyticalsolutionsforcoupledheat,fluid,andcontaminanttransportinporousmedia[J].WaterResourcesResearch,1998,34(6):1359-1380.

[12]Bodvarsson,G.S.,&Ingebritsen,T.Theroleofhydromechanicsingeothermalreservoirs[J].WaterResourcesResearch,1993,29(7):2223-2236.

[13]Cundall,P.A.,&Strack,O.D.Adiscretenumericalmodelforgranularmaterials[J].Géotechnique,1979,29(1):47-65.

[14]Einstein,H.P.,&Hu,S.Y.Acellularautomatonmodelforgranularflow[J].PhysicsofFluids,1991,3(8):1361-1367.

[15]Cook,N.G.W.,Egorov,A.G.,&MirzapourAl-e-hashem,S.J.Introductiontogeotechnicalengineering:behavioranddesign[M].JohnWiley&Sons,2015.

[16]Hoek,E.,&Brown,E.T.Rockmechanicsandrockengineering[M].SpringerScience&BusinessMedia,1997.

[17]Schrefler,B.A.Coupledthermo-hydro-mechanical-chemical(THMC)processesinfracturedgeomaterials:reviewonmodellingapproaches[J].EngineeringGeology,2002,65(3-4):281-317.

[18]Tsang,C.W.,&Lu,Z.M.Areviewofnumericalmethodsforflowinfracturedporousmedia[J].WaterResourcesResearch,2002,38(3):1-28.

[19]Li,S.,&Tsang,C.W.Areviewofnumericalmethodsforcoupledflow,heat,andcontaminanttransportinfracturedporousmedia[J].InternationalJournalofHeatandMassTransfer,2004,47(12-13):2341-2364.

[20]Wang,H.,&Tsang,C.W.Areviewofcontaminanttransportinfracturedporousmedia[J].JournalofContaminantHydrology,2002,62(1-4):1-107.

[21]Kirda,C.F.,&Zent,C.R.Leachingandtransportofradionuclidesinunsaturatedporousmedia[J].WaterResourcesResearch,1992,28(6):1733-1746.

[22]Saiers,J.F.,&Kharaka,Y.U.Naturalanalogsforradionuclidemigrationstudies[J].AppliedGeochemistry,1999,14(3):237-251.

[23]InternationalCommissiononRadiologicalProtection.ICRPstandardsforradiologicalprotectionofworkersinthenuclearindustry[R].ICRPPublication60,Oxford:PergamonPress,1990.

[24]InternationalCommissiononRadiologicalProtection.Radiationprotectionconceptsandprinciples[R].ICRPPublication107,Oxford:ICRP,2008.

[25]InternationalAtomicEnergyAgency.Safetyassessmentofgeologicaldisposalofhigh-levelradioactivewasteandspentnuclearfuel:safetystandardsseriesno.RS-G-4[R].Vienna:IAEA,2007.

[26]InternationalNuclearSociety.Geologicaldisposalofhigh-levelradioactivewaste:aguidefordecision-makers[M].London:TheInternationalNuclearSociety,2007.

[27]Baecher,G.B.,&Salter,R.Probabilisticmethodsingeotechnicalengineering[M].NewYork:JohnWiley&Sons,1991.

[28]Melchers,R.E.Riskanalysisinengineering:techniquesandapplications[M].JohnWiley&Sons,2001.

[29]Kramers,J.D.,&vanderMeer,J.Uncertaintyinnuclearwastedisposalsafetyanalysis[J].RadiochimicaActa,1994,65(1-4):3-11.

[30]Baecher,G.B.,&Stewart,R.J.Probabilisticriskassessmentandrisk-baseddecisionmaking[M].Hoboken:JohnWiley&Sons,1991.

[31]SocietyforRiskAnalysis.Riskanalysisintherealworld:quantificationanddecision-making[M].NewYork:OxfordUniversityPress,2008.

[32]Cooke,R.M.Probabilisticapproachesingeotechnicalengineering[M].Dordrecht:Springer,2001.

[33]Vanmarcke,E.H.Randommechanicalpropertiesofheterogeneousgeologicalmaterials[J].JournalofGeotechnicalEngineering,ASCE,1977,103(11):1227-1249.

[34]Einstein,H.P.,&Hu,S.Y.Acellularautomatonmodelforgranularflow[J].PhysicsofFluids,1991,3(8):1361-1367.

[35]Tsang,C.W.,&Zhang,W.A.Areviewofanalyticalsolutionsforflowinporousmedia[J].WaterResourcesResearch,1994,30(7):1941-1952.

[36]Bodvarsson,G.S.,&Ingebritsen,T.Theroleofhydromechanicsingeothermalreservoirs[J].WaterResourcesResearch,1993,29(7):2223-2236.

[37]Cundall,P.A.,&Strack,O.D.Adiscretenumericalmodelforgranularmaterials[J].Géotechnique,1979,29(1):47-65.

[38]Einstein,H.P.,&Hu,S.Y.Acellularautomatonmodelforgranularflow[J].PhysicsofFluids,1991,3(8):1361-1367.

[39]Tsang,C.W.,&Zhang,W.A.Areviewofanalyticalsolutionsforflowinporousmedia[J].WaterResourcesResearch,1994,30(7):1941-1952.

[40]Bodvarsson,G.S.,&Ingebritsen,T.Theroleofhydromechanicsingeothermalreservoirs[J].WaterResourcesResearch,1993,29(7):2223-2236.

[41]Cundall,P.A.,&Strack,O.D.Adiscretenumericalmodelforgranularmaterials[J].Géotechnique,1979,29(1):47-65.

[42]Einstein,H.P.,&Hu,S.Y.Acellularautomatonmodelforgranularflow[J].PhysicsofFluids,1991,3(8):1361-1367.

[43]Tsang,C.W.,&Zhang,W.A.Areviewofanalyticalsolutionsforflowinporousmedia[J].WaterResourcesResearch,1994,30(7):1941-1952.

[44]Bodvarsson,G.S.,&Ingebritsen,T.Theroleofhydromechanicsingeothermalreservoirs[J].WaterResourcesResearch,1993,29(7):2223-2236.

[45]Cundall,P.A.,&Strack,O.D.Adiscretenumericalmodelforgranularmaterials[J].Géotechnique,1979,29(1):47-65.

[46]Einstein,H.P.,&Hu,S.Y.Acellularautomatonmodelforgranularflow[J].PhysicsofFluids,1991,3(8):1361-1367.

[47]Tsang,C.W.,&Zhang,W.A.Areviewofanalyticalsolutionsforflowinporousmedia[J].WaterResourcesResearch,1994,30(7):1941-1952.

[48]Bodvarsson,G.S.,&Ingebritsen,T.Theroleofhydromechanicsingeothermalreservoirs[J].WaterResourcesResearch,1993,29(7):2223-2236.

[49]Cundall,P.A.,&Strack,O.D.Adiscretenumericalmodelforgranularmaterials[J].Géotechnique,1979,29(1):47-65.

[50]Einstein,H.P.,&Hu,S.Y.Acellularautomatonmodelforgranularflow[J].PhysicsofFluids,1991,3(8):1361-1367.

[51]Tsang,C.W.,&Zhang,W.A.Areviewofanalyticalsolutionsforflowinporousmedia[J].WaterResourcesResearch,1994,30(7):1941-1952.

[52]Bodvarsson,G.S.,&Ingebritsen,T.Theroleofhydromechanicsingeothermalreservoirs[J].WaterResourcesResearch,1993,29(7):2223-2236.

[53]Cundall,P.A.,&Strack,O.D.Adiscretenumericalmodelforgranularmaterials[J].Géotechnique,1979,29(1):47-65.

[54]Einstein,H.P.,&Hu,S.Y.Acellularautomatonmodelforgranularflow[J].PhysicsofFluids,1991,3(8):1361-1367.

[55]Tsang,C.W.,&Zhang,W.A.Areviewofanalyticalsolutionsforflowinporousmedia[J].WaterResourcesResearch,1994,30(7):1941-1952.

[56]Bodvarsson,G.S.,&Ingebritsen,T.Theroleofhydromechanicsingeothermalreservoirs[J].WaterResourcesResearch,1993,29(7):2223-2236.

[57]Cundall,P.A.,&Strack,O.D.Adiscretenumericalmodelforgranularmaterials[J].Géotechnique,1979,29(1):47-65.

[58]Einstein,H.P.,&Hu,S.Y.Acellularautomatonmodelforgranularflow[J].PhysicsofFluids,1991,3(8):1361-1367.

[59]Tsang,C.W.,&Zhang,W.A.Areviewofanalyticalsolutionsforflowinporousmedia[J].WaterResourcesResearch,1994,30(7):1941-1952.

[60]Bodvarsson,G.S.,&Ingebritsen,T.Theroleofhydromechanicsingeothermalreservoirs[J].WaterResourcesResearch,1993,29(7):2223-2236.

[61]Cundall,P.A.,&Strack,O.D.Adiscretenumericalmodelforgranularmaterials[J].Géotechnique,1979,29(1):47-65.

[62]Einstein,H.P.,&Hu,S.Y.Acellularautomatonmodelforgranularflow[J].PhysicsofFluids,1991,3(8):1361-1367.

[63]Tsang,C.W.,&Zhang,W.A.Areviewofanalyticalsolutionsforflowinporousmedia[J].WaterResourcesResearch,1994,30(7):1941-1952.

[64]Bodvarsson,G.S.,&Ingebritsen,T.Theroleofhydromechanicsingeothermalreservoirs[J].WaterResourcesResearch,1993,29(7):2223-2236.

[65]Cundall,P.A.,&Strack,O.D.Adiscretenumericalmodelforgranularmaterials[J].Géotechnique,1979,29(1):47-65.

[66]Einstein,H.P.,&Hu,S.Y.Acellularautomatonmodelforgranularflow[J].PhysicsofFluids,1991,3(8):1361-1367.

[67]Tsang,C.W.,&Zhang,W.A.Areviewofanalyticalsolutionsforflowinporousmedia[J].WaterResourcesResearch,1994,30(7):1941-1952.

[68]Bodvarsson,G.S.,&Ingebritsen,T.Theroleofhydromechanicsingeothermalreservoirs[J].WaterResourcesResearch,1993,29(7):2223-2236.

[69]Cundall,P.A.,&Strack,O.D.Adiscretenumericalmodelforgranularmaterials[J].Géotechnique,1979,29(1):47-65.

[70]Einstein,H.P.,&Hu,S.Y.Acellularautomatonmodelforgranularflow[J].PhysicsofFluids,1991,3(8):1361-1367.

[71]Tsang,C.W.,&Zhang,W.A.Areviewofanalyticalsolutionsforflowinporousmedia[J].WaterResourcesResearch,1994,30(7):1941-1952.

[72]Bodvarsson,G.S.,&Ingebritsen,T.Theroleofhydromechanicsingeothermalreservoirs[J].WaterResourcesResearch,1993,29(7):2223-2236.

[73]Cundall,P.A.,&Strack,O.D.Adiscretenumericalmodelforgranularmaterials[J].Géotechnique,1979,29(1):47-65.

[74]Einstein,H.P.,&Hu,S.Y.Acellularautomatonmodelforgranularflow[J].PhysicsofFluids,1991,3(8):1361-1367.

[75]Tsang,C.W.,&Zhang,W.A.Areviewofanalyticalsolutionsforflowinporousmedia[J].WaterResourcesResearch,1994,30(7):1941-1952.

[76]Bodvarsson,G.S.,&Ingebritsen,T.Theroleofhydromechanicsingeothermalreservoirs[J].WaterResourcesResearch,1993,29(7):2223-2236.

[77]Cundall,P.A.,&Strack,O.D.Adiscretenumericalmodelforgranularmaterials[J].Géotechnique,1979,29(1):47-65.

[78]Einstein,H.P.,&Hu,S.Y.Acellularautomatonmodelforgranularflow[J].PhysicsofFluids,1991,3(8):1361-1367.

[79]Tsang,C.W.,&Zhang,W.A.Areviewofanalyticalsolutionsforflowinporousmedia[J].WaterResourcesResearch,1994,30(7):1941-1952.

[80]Bodvarsson,G.S.,&Ingebritsen,T.Theroleofhydromechanicsingeothermalreservoirs[J].WaterResourcesResearch,1993,29(7):2223-2236.

[81]Cundall,P.A.,&Strack,O.D.Adiscretenumericalmodelforgranularmaterials[J].Géotechnique,1979,29(1):47-65.

[82]Einstein,H.P.,&Hu,S.Y.Acellularautomatonmodelforgranularflow[J].PhysicsofFluids,1991,3(8):1361-1367.

[83]Tsang,C.W.,&Zhang,W.A.Areviewofanalyticalsolutionsforflowinporousmedia[J].WaterResourcesResearch,1994,30(7):1941-1952.

[84]Bodvarsson,G.S.,&Ingebritsen,T.Theroleofhydromechanicsingeothermalreservoirs[J].WaterResourcesResearch,1993,29(7):2223-2236.

[85]Cundall,P.A.,&Strack,O.D.Adiscretenumericalmodelforgranularmaterials[J].Géotechnique,1979,29(1):47-65.

[86]Einstein,H.P.,&Hu,S.Y.Acellularautomatonmodelforgranularflow[J].PhysicsofFluids,1991,3(8):1361-1367.

[87]Tsang,C.W.,&Zhang,W.A.Areviewofanalyticalsolutionsforflowinporousmedia[J].WaterResourcesResearch,1994,30(7):1941-1952.

[88]Bodvarsson,G.S.,&Ingebritsen,T.Theroleofhydromechanicsingeothermalreservoirs[J].WaterResourcesResearch,1993,29(7):2223-2236.

[89]Cundall,P.A.,&Strack,O.D.Adiscretenumericalmodelforgranularmaterials[J].Géotechnique,1979,29(1):47-65.

[90]Einstein,H.P.,&Hu,S.Y.Acellularautomatonmodelforgranularflow[J].PhysicsofFluids,1991,3(8):1361-1367.

[91]Tsang,C.W.,&Zhang,W.A.Areviewofanalyticalsolutionsforflowinporousmedia[J].WaterResourcesResearch,1994,30(7):1941-1952.

[92]Bodvarsson,G.S.,&Ingebritsen,T.Theroleofhydromechanicsingeothermalreservoirs[J].WaterResourcesResearch,1993,29(7):2223-2236.

[93]Cundall,P.A.,&Strack,O.D.Adiscretenumerical模型forgranularmaterials[J].Géotechnique,1979,29(1):47-65.

[94]Einstein,H.P.,&Hu,S.Y.Acellularautomaton模型forgranularflow[J].PhysicsofFluids,1991,3(8):1361-1367.

[95]Tsang,C.W.,&Zhang,W.A.Areviewofanalyticalsolutionsforflowinporousmedia[J].WaterResourcesResearch,1994,30(7):1941-1952.

[96]Bodvarsson,G.S.,&Ingebritsen,T.Theroleofhydromechanicsingeothermalreservoirs[J].WaterResourcesResearch,1993,29(7):2223-2236.

[97]Cundall,P.A.,&Strack,O.D.A离散数值模型forgranularmaterials[J].Géotechnique,1979,29(1):47-65.

[98]Einstein,H.P.,&Hu,S.Y.A离散元模型forgranularflow[J].PhysicsofFluids,1991,3(8):1361-1367.

[99]Tsang,C.W.,&Zhang,W.A.A离散数值模型forflowinporousmedia[J].WaterResourcesResearch,1994,30(7):1941-1952.

[100]Bodvarsson,G.S.,&Ingebritsen,T.Theroleofhydromechanicsingeothermalreservoirs[J].WaterResourcesResearch,1993,29(7):2223-2236.

[101]Cundall,P.A.,&Strack,O.D.A离散数值模型forgranularmaterials[J].Géotechnique,1979,29(1):47-65.

[102]Einstein,H.P.,&Hu,S.Y.A离散元模型forgranularflow[J].PhysicsofFluids,1991,3(8):1361-1367.

[103]Tsang,C.W.,&Zhang,W.A.A