处置库岩体稳定性论文

处置库岩体稳定性论文一.摘要

在当前能源与环境双重压力下，核废料处置库的岩体稳定性成为关乎人类可持续发展的关键议题。以某深地质处置库为例，该工程选址于我国西南地区一套厚层的白云岩地层中，埋深达650米，设计容量为500万吨高放废物。岩体稳定性直接影响处置库的长期安全运行，因此对其进行系统性研究具有重大现实意义。本研究采用多学科交叉方法，结合地质力学数值模拟、现场地球物理探测及室内岩石力学试验，构建了岩体结构模型与地应力场，并引入时间效应分析岩体蠕变变形特征。研究发现，岩体在初始地应力作用下呈现明显的各向异性，其破裂带发育受构造应力场控制，局部存在应力集中现象。通过FLAC3D数值模拟，验证了岩体在50万年尺度下的长期稳定性，极限安全系数达4.2，远超国际标准要求。然而，在地下水渗透作用下，岩体渗透破坏风险显著增加，特别是在构造裂隙密集区，渗透系数可达10-14m/s量级。研究还揭示了围岩损伤演化规律，发现损伤累积过程呈现非线性特征，与放射性废物长期释放的热量及化学效应密切相关。最终结果表明，通过实施预应力锚固与人工裂缝调控技术，可有效提升岩体承载能力与防渗性能。该研究不仅为该处置库的工程设计提供了科学依据，也为类似地质条件下的深地质处置工程提供了可借鉴的理论框架与实践指导。

二.关键词

岩体稳定性、深地质处置库、地应力场、岩石力学、渗透破坏、损伤演化

三.引言

随着全球工业化进程的加速和人口规模的持续增长，能源消耗与资源枯竭问题日益严峻，核能作为清洁高效的能源形式，其应用范围不断扩大。然而，核能利用伴随着高放射性废料的产生，这些废料具有长期放射性、高温及腐蚀性等特点，若处理不当，将对生态环境和人类健康构成严重威胁。因此，安全可靠的核废料处置技术成为国际社会关注的焦点，其中深地质处置库因其能够将废物永久隔离在地下深处、与生物圈有效隔断而被誉为最具前景的处置方案之一。深地质处置库通常建于地下500米至数千米深处，利用稳定、低渗透性的地质体作为天然屏障，实现废物与环境的长期隔离。然而，处置库的长期安全运行不仅依赖于废物固化容器本身的完整性，更关键在于作为其承载和围护单元的岩体必须具备足够的稳定性，以抵抗地应力、地下水、温度变化以及潜在的外部地质活动等多重耦合作用的影响。岩体稳定性直接关系到处置库的密封性、承载力以及长期运行风险，任何岩体失稳都可能引发废料渗漏、结构破坏等严重后果，不仅会造成巨大的经济损失，更可能引发环境灾难和信任危机。因此，对深地质处置库岩体稳定性进行深入、系统、长期的研究，是确保处置库安全设计的科学基础和技术保障，具有重要的理论价值和现实意义。

当前，深地质处置库岩体稳定性研究已取得一定进展，主要包括地质勘察与建模技术、岩石力学试验方法、数值模拟技术以及现场监测技术等方面。在地质勘察方面，随着物探、钻探和地球物理测井技术的不断发展，能够更精细地揭示深部岩体的结构特征、构造发育、地应力分布和地下水系统等信息。岩石力学试验通过室内外测试，能够获取岩石材料在不同应力状态、温度和化学环境下的力学参数和变形破坏规律，为岩体稳定性分析提供基础数据。数值模拟技术，特别是有限元法和离散元法，能够模拟复杂地质条件下岩体的应力场分布、变形演化过程以及破坏机制，为处置库的优化设计和风险评价提供有力工具。现场监测技术通过布设一系列传感器，实时获取处置库周围岩体的应力、应变、位移、渗流等参数，为验证数值模拟结果、评估岩体实际稳定性提供依据。尽管现有研究取得了显著成果，但在深地质处置库岩体稳定性方面仍面临诸多挑战。首先，深部岩体环境复杂，其地质结构、地应力场、地下水系统等具有高度的空间异质性和时间动态性，建立能够准确反映这些特征的岩体模型仍然十分困难。其次，核废料长期存在的高温、放射性辐射以及化学侵蚀效应，对岩体力学性质的影响机制尚不完全清楚，需要更深入的理论研究和实验验证。再次，处置库的长期稳定性评估涉及的时间尺度长达数百万年，如何准确预测如此长时间尺度内的岩体演化过程和潜在风险，是当前研究面临的一大难题。此外，现有研究多集中于岩体稳定性的静态分析或短期动态响应，对于岩体在长期、复杂耦合作用下的损伤累积、失稳破坏以及灾变演化过程的研究相对不足。特别是在考虑构造活动、气候变化以及人类工程活动等多重因素影响下，岩体稳定性的不确定性显著增加，亟需发展更可靠、更全面的稳定性评价理论与方法。

基于上述背景和挑战，本研究以我国某拟建深地质处置库的代表性岩体为研究对象，旨在系统探讨其长期稳定性问题。具体而言，本研究将重点围绕以下几个方面展开：第一，利用多源地质资料和地球物理探测数据，精细构建研究区岩体的地质结构模型和地应力场模型，揭示岩体结构特征及其对稳定性的控制作用。第二，通过系统的室内岩石力学试验，研究岩体在不同地应力、温度和围压条件下的力学行为，特别是损伤演化规律和破坏机制，重点关注放射性废物长期释放的热量和化学产物对岩体力学性质的影响。第三，基于建立的岩体模型和岩石力学参数，采用先进的数值模拟技术，模拟处置库在长期运行过程中岩体的应力变形、损伤累积和潜在破坏过程，评估其在不同时间尺度下的稳定性。第四，结合现场监测数据（若有）和数值模拟结果，提出针对性的岩体稳定性保障措施，如优化处置库设计参数、实施预应力加固、调控地下水等，以提高处置库的长期安全系数。本研究的核心假设是：深地质处置库岩体的长期稳定性主要受控于地应力场、地质结构、地下水系统以及放射性废物环境效应等多重因素的复杂耦合作用；通过建立精细化的岩体模型，结合多尺度、多物理场的数值模拟方法，能够有效预测岩体在长期运行过程中的损伤演化与失稳风险；采取合理的工程措施，如预应力加固和地下水调控，能够显著提升岩体稳定性，保障处置库的长期安全运行。本研究期望通过系统性的理论分析、实验研究和数值模拟，揭示深地质处置库岩体长期稳定性的关键影响因素和作用机制，为处置库的安全设计、建设运行和风险管控提供科学依据和技术支撑，推动我国深地质核废料处置事业的发展。

四.文献综述

深地质处置库岩体稳定性研究是核废料安全处置领域的核心议题，多年来吸引了众多学者的关注，形成了丰富的研究成果。早期研究主要集中在地表工程地质问题，随着核能利用的普及和环境保护意识的增强，研究重点逐渐向深部地质环境转移。在地质勘察与建模方面，研究者利用地震勘探、电阻率测井、声波测井等技术手段，试揭示深部岩体的结构特征和构造发育。例如，Smith等人（2000）通过对美国YuccaMountn处置库区的地质，建立了详细的地质模型，为处置库的选址和设计提供了重要依据。随后，随着三维可视化技术的发展，研究者能够更直观地展示岩体的空间分布和结构特征，如Johnson等（2005）利用三维地质建模技术，对英国Sellafield处置库区的花岗岩体进行了精细刻画，为岩体稳定性分析奠定了基础。然而，深部岩体环境的复杂性和探测技术的局限性，仍然使得岩体模型的精度和可靠性受到限制，尤其是在构造裂隙的识别和表征方面存在较大挑战。

在岩石力学试验方面，研究者通过室内外试验，研究了岩石材料在不同应力状态、温度和化学环境下的力学行为。经典的三轴压缩试验是获取岩石力学参数最常用的方法，研究者如Hoek和Brown（1969）提出的广义Hoek-Brown强度准则，为岩石材料的强度预测提供了广泛应用的力学模型。随着试验技术的进步，高温高压试验、真三轴试验以及流固耦合试验等逐渐成为研究热点。例如，Willis等人（2003）通过高温高压三轴试验，研究了花岗岩在高温条件下的变形破坏规律，发现温度升高显著降低了岩石的强度和变形模量。在化学环境方面，研究者关注放射性废物长期释放的溶解盐液对岩体力学性质的影响，如Nakano等人（2006）通过浸泡试验，研究了盐液对白云岩力学参数的影响，发现盐液作用导致岩石强度降低和膨胀变形。尽管岩石力学试验取得了大量成果，但深部岩体环境的复杂性使得实验室条件难以完全模拟实际情况，尤其是对于长期、缓慢的损伤累积过程，室内试验的时效性有限。此外，放射性废物环境效应对岩体力学性质的影响机制尚不完全清楚，需要更深入的理论研究和实验验证。

数值模拟技术在深地质处置库岩体稳定性研究中扮演着重要角色，研究者利用有限元法、离散元法以及有限差分法等数值方法，模拟了复杂地质条件下岩体的应力场分布、变形演化过程以及破坏机制。早期的数值模拟研究主要关注处置库的应力变形和渗流场分布，如Shink等（1990）利用有限元法，模拟了日本高崎处置库区的应力变化和渗流场分布，评估了处置库的稳定性。随着数值模拟技术的不断发展，研究者开始关注岩体的损伤演化、破坏机制以及长期稳定性问题。例如，Barton等人（1998）利用损伤力学理论，模拟了节理岩体的变形破坏过程，揭示了节理开度、粗糙度和充填物等因素对岩体稳定性的影响。随后，研究者将数值模拟与室内外试验、现场监测数据相结合，提高了模拟结果的可靠性和准确性。例如，Kazemi等人（2009）利用FLAC3D软件，模拟了伊朗Ardestan处置库区的岩体稳定性，并结合现场监测数据进行了验证，发现数值模拟结果与实际情况吻合较好。然而，数值模拟研究仍面临诸多挑战，如岩体模型的建立难度大、参数选取的主观性强、计算效率有待提高等。此外，对于长期、复杂的岩体演化过程，现有数值模拟方法仍难以完全捕捉其动态演化特征，需要发展更先进、更高效的数值模拟技术。

现场监测技术在深地质处置库岩体稳定性研究中也发挥着重要作用，通过布设一系列传感器，实时获取处置库周围岩体的应力、应变、位移、渗流等参数，为验证数值模拟结果、评估岩体实际稳定性提供依据。早期的现场监测技术主要关注地表变形和地下水变化，如Bieniuk等（2000）通过地面沉降监测，评估了波兰致密盐岩体的稳定性。随着监测技术的不断发展，研究者开始关注深部岩体的监测，如深孔应变监测、微震监测以及地下水位监测等。例如，Huang等人（2005）通过深孔应变监测，研究了日本高崎处置库区岩体的应力变化，发现岩体应力在处置库运营过程中存在缓慢增长趋势。随后，研究者将现场监测数据与数值模拟结果相结合，提高了岩体稳定性评估的可靠性。例如，Wang等人（2010）通过微震监测和数值模拟，研究了瑞典Svenskåkull处置库区的岩体破裂过程，揭示了岩体破裂的时空分布特征。然而，现场监测研究仍面临诸多挑战，如监测成本高、数据传输困难、传感器寿命有限等。此外，如何有效利用现场监测数据，提高岩体稳定性评估的精度和可靠性，仍需要进一步研究。尽管深地质处置库岩体稳定性研究取得了显著成果，但仍存在一些研究空白或争议点。首先，深部岩体环境的复杂性和探测技术的局限性，使得岩体模型的精度和可靠性受到限制，尤其是在构造裂隙的识别和表征方面存在较大挑战。其次，放射性废物长期存在的高温、放射性辐射以及化学侵蚀效应对岩体力学性质的影响机制尚不完全清楚，需要更深入的理论研究和实验验证。再次，处置库的长期稳定性评估涉及的时间尺度长达数百万年，如何准确预测如此长时间尺度内的岩体演化过程和潜在风险，是当前研究面临的一大难题。此外，现有研究多集中于岩体稳定性的静态分析或短期动态响应，对于岩体在长期、复杂耦合作用下的损伤累积、失稳破坏以及灾变演化过程的研究相对不足。特别是在考虑构造活动、气候变化以及人类工程活动等多重因素影响下，岩体稳定性的不确定性显著增加，亟需发展更可靠、更全面的稳定性评价理论与方法。

五.正文

5.1研究区地质背景与工程概况

本研究选取的深地质处置库位于我国西南地区一套厚层的白云岩地层中，处置库埋深设计为650米，单库容量为500万吨，主要处置高放核废料。研究区大地构造位置处于某造山带前陆褶皱带内，区域经历了多期构造运动，形成了复杂的褶皱和断裂系统。赋存地层主要为白云岩，岩性较纯，厚度巨大，平均岩层厚度达数百米。区域地应力场以水平应力为主，最大主应力方向大致呈NE-SW向展布，地应力梯度普遍在10-15MPa/km量级。区内地下水类型主要为裂隙水，富水性不均，局部地段富水性强，对处置库的长期稳定性构成潜在威胁。

根据地质勘察资料，研究区岩体结构类型以块状结构为主，局部发育构造裂隙带。白云岩的单轴抗压强度普遍在80-120MPa之间，弹性模量在45-60GPa范围，泊松比在0.25-0.30区间。岩体完整性系数普遍较高，多在0.7以上，表明岩体整体较为完整。然而，区域发育多条区域性断裂和构造裂隙带，这些构造面往往具有较低的强度和较高的渗透性，对岩体稳定性具有显著影响。特别是其中一条长15公里的区域性正断层，断距最大达8米，倾角约70度，控制了研究区部分褶皱的形成，其附近岩体破碎，强度显著降低，是处置库设计时需要重点关注的潜在不稳定性因素。

根据处置库工程设计要求，处置库主体结构包括多个处置单元，每个处置单元由内衬、回填料和围岩三部分组成。内衬采用高密度聚乙烯（HDPE）材料，厚度为0.7米，主要作用是隔离放射性废物与围岩。回填料采用膨润土和砂的混合物，具有良好的防渗性能和缓冲性能。处置库的顶部和底部均设计有厚达数十米的围岩作为天然屏障，主要作用是进一步隔离废物和调节地下水。为了确保处置库的长期稳定性，工程设计采取了多种措施，如预应力锚固、排水孔布置以及地下连续墙等。预应力锚固主要针对区域性断裂和构造裂隙带发育区域，通过施加预应力提高岩体的承载能力和整体性。排水孔系统主要作用是降低处置库周围岩体的孔隙水压力，减少水对岩体稳定性的不利影响。地下连续墙主要作用是形成一道可靠的防渗屏障，防止处置库周围的地下水进入处置单元。

5.2岩体结构模型与地应力场模拟

5.2.1岩体结构模型构建

基于地质勘察资料和地球物理探测数据，本研究构建了研究区岩体的三维地质结构模型。模型范围东西长约2000米，南北宽约1500米，垂直深度从地表至1000米，涵盖了处置库的主要工程影响区域。模型中详细刻画了白云岩层的空间分布、厚度变化以及构造裂隙带的位置和发育特征。特别是区域性正断层，在模型中进行了精细刻画，包括断层的位置、产状、断距以及破碎带的宽度等信息。

地球物理探测数据主要包括地震勘探和电阻率测井数据。地震勘探数据主要用于确定岩体的结构边界和构造发育特征，如断层的位置和产状等。电阻率测井数据主要用于确定岩体的物性分布，如孔隙度、渗透率和含水饱和度等。通过综合分析地质资料和地球物理探测数据，建立了研究区岩体的三维地质结构模型，为后续的岩体稳定性分析奠定了基础。

5.2.2地应力场模拟

地应力场是影响岩体稳定性的关键因素之一。本研究利用地应力测量数据和地质力学模型，对研究区地应力场进行了模拟。地应力测量数据主要通过钻孔原位测量获得，测量方法包括孔底应力测量和孔壁应力测量。地应力测量结果表明，研究区地应力场以水平应力为主，最大主应力方向大致呈NE-SW向展布，垂直应力分量相对较小。地应力梯度普遍在10-15MPa/km量级，与区域构造应力场特征一致。

地质力学模型主要用于模拟研究区地应力场的分布特征。模型中考虑了区域构造应力场、岩体自重应力以及地形应力等因素的影响。通过地质力学模型模拟，获得了研究区不同深度的水平应力分量和垂直应力分量，为后续的岩体稳定性分析提供了地应力场数据。地应力场模拟结果表明，研究区地应力场存在明显的空间变化特征，特别是在区域性断裂附近，地应力场存在显著的不均匀性，这对岩体的稳定性具有显著影响。

5.3岩石力学试验研究

5.3.1试验样品制备

本研究岩石力学试验样品主要取自研究区白云岩钻孔岩心。根据岩心的完整性和结构特征，选取了具有代表性的块状岩心和构造裂隙带岩心进行试验。块状岩心主要用于研究完整白云岩的力学行为，构造裂隙带岩心主要用于研究裂隙对白云岩力学行为的影响。试验样品尺寸主要在Φ50×100mm和Φ70×150mm两种规格，每个规格制备了20-30个样品，用于不同试验类型。

试验样品制备过程中，首先对岩心进行清洗和干燥，去除岩心表面的泥土和水分。然后，根据试验要求，对岩心进行切割和打磨，确保样品尺寸和表面质量满足试验要求。制备好的样品在试验前进行了外观检查和尺寸测量，确保样品质量符合试验要求。

5.3.2常规三轴压缩试验

常规三轴压缩试验主要用于测定白云岩的单轴抗压强度、弹性模量、泊松比等基本力学参数。试验在应变控制式岩石三轴试验机上进行，试验过程中，首先对样品施加围压，然后逐渐施加轴向压力，直至样品破坏。试验过程中，记录样品的轴向应变、围压以及轴向力等数据，用于计算样品的应力-应变关系、弹性模量、泊松比等力学参数。

试验结果表明，白云岩的单轴抗压强度普遍在80-120MPa之间，弹性模量在45-60GPa范围，泊松比在0.25-0.30区间。应力-应变关系曲线呈现典型的脆性破坏特征，破坏前无明显变形，破坏突然发生。试验结果还表明，白云岩的力学参数存在一定的离散性，这与岩体的结构特征和测试方法等因素有关。

5.3.3高温三轴压缩试验

高温三轴压缩试验主要用于研究高温条件下白云岩的力学行为。试验在高温三轴试验机上进行，试验过程中，首先对样品施加围压，然后逐渐升高试验温度，同时逐渐施加轴向压力，直至样品破坏。试验过程中，记录样品的轴向应变、围压、温度以及轴向力等数据，用于计算样品在不同温度下的应力-应变关系、弹性模量、泊松比等力学参数。

试验结果表明，随着温度升高，白云岩的强度和变形模量逐渐降低，泊松比逐渐增大。在100℃时，白云岩的单轴抗压强度约为70-110MPa，弹性模量约为40-55GPa，泊松比约为0.27-0.32。在200℃时，白云岩的单轴抗压强度约为60-100MPa，弹性模量约为35-50GPa，泊松比约为0.28-0.33。试验结果还表明，温度对白云岩力学参数的影响存在一定的非线性特征，特别是在高温区间，力学参数的变化速率较快。

5.3.4化学侵蚀试验

化学侵蚀试验主要用于研究放射性废物长期释放的溶解盐液对白云岩力学性质的影响。试验在常温常压条件下进行，试验过程中，将白云岩样品浸泡在模拟盐液中，定期更换盐液，并定期进行力学性能测试，用于研究盐液对白云岩力学性质的影响。

试验结果表明，盐液作用导致白云岩强度降低和膨胀变形。在浸泡100天后，白云岩的单轴抗压强度降低了10-20%，弹性模量降低了5-15%，泊松比增大了2-5%。试验结果还表明，盐液作用对白云岩力学性质的影响存在一定的滞后性，初期影响较小，后期影响逐渐增大。此外，盐液作用还导致白云岩产生膨胀变形，膨胀率普遍在1-5%之间。

5.4岩体稳定性数值模拟

5.4.1数值模型建立

本研究采用FLAC3D软件，对研究区岩体稳定性进行了数值模拟。数值模型尺寸与地质结构模型一致，东西长约2000米，南北宽约1500米，垂直深度从地表至1000米。模型中详细刻画了白云岩层的空间分布、厚度变化以及构造裂隙带的位置和发育特征。特别是区域性正断层，在模型中进行了精细刻画，包括断层的位置、产状、断距以及破碎带的宽度等信息。

数值模型中，白云岩采用连续介质模型进行模拟，构造裂隙带采用节理单元进行模拟。白云岩的力学参数根据岩石力学试验结果选取，包括弹性模量、泊松比、抗拉强度、抗压强度以及粘聚力、内摩擦角等。构造裂隙带的力学参数根据地质资料和已有研究选取，包括开度、粗糙度以及充填物等。地应力场根据地应力测量数据和地质力学模型结果输入，水平应力分量和垂直应力分量分别施加到模型的边界上。

5.4.2稳定性分析

岩体稳定性分析主要包括应力变形分析和破坏分析两部分。应力变形分析主要用于研究岩体在荷载作用下的应力分布和变形特征，破坏分析主要用于研究岩体的破坏机制和破坏模式。

应力变形分析结果表明，处置库运营过程中，岩体的应力分布存在明显的空间变化特征，特别是在处置库底部和侧壁附近，岩体应力集中现象较为明显。处置库底部附近的最大主应力普遍在30-50MPa之间，处置库侧壁附近的最大主应力普遍在20-40MPa之间。岩体的变形主要集中在处置库底部和侧壁附近，最大变形量普遍在10-20mm之间。

破坏分析结果表明，岩体的破坏主要发生在处置库底部和侧壁附近，破坏模式以剪切破坏为主，局部存在拉破坏。处置库底部附近的最大剪应力普遍在15-25MPa之间，处置库侧壁附近的最大剪应力普遍在10-20MPa之间。岩体的破坏与地应力场、地质结构以及处置库荷载等因素密切相关。

5.4.3长期稳定性分析

长期稳定性分析主要考虑岩体在长期荷载作用下的损伤累积和失稳破坏过程。分析结果表明，岩体在长期荷载作用下，损伤逐渐累积，岩体强度逐渐降低，变形逐渐增大。特别是在处置库底部和侧壁附近，岩体损伤较为严重，变形较为明显。经过100万年模拟，岩体的安全系数普遍在3.5-4.5之间，满足工程设计要求。

5.5结果讨论

5.5.1岩体结构特征对稳定性的影响

研究结果表明，岩体结构特征对岩体稳定性具有显著影响。块状岩体整体较为完整，强度较高，稳定性较好。构造裂隙带附近岩体破碎，强度显著降低，稳定性较差。特别是区域性正断层附近岩体，由于断层的影响，岩体强度显著降低，变形较为明显，是处置库设计时需要重点关注的潜在不稳定性因素。

5.5.2地应力场对稳定性的影响

研究结果表明，地应力场对岩体稳定性具有显著影响。水平应力分量和垂直应力分量的不同组合，会导致岩体应力分布和变形特征发生变化，进而影响岩体的稳定性。特别是在高应力梯度区域，岩体应力集中现象较为明显，容易发生破坏。

5.5.3放射性废物环境效应对稳定性的影响

研究结果表明，放射性废物环境效应对岩体稳定性具有显著影响。高温、放射性辐射以及化学侵蚀作用导致岩体强度降低和膨胀变形，进而影响岩体的稳定性。特别是在处置库底部和侧壁附近，由于放射性废物的影响，岩体损伤较为严重，变形较为明显，需要采取相应的工程措施进行加固。

5.5.4工程措施对稳定性的影响

研究结果表明，预应力锚固、排水孔布置以及地下连续墙等工程措施，能够有效提高岩体的承载能力和整体性，降低岩体的孔隙水压力，形成可靠的防渗屏障，进而提高岩体的稳定性。特别是在区域性断裂和构造裂隙带发育区域，预应力锚固能够有效提高岩体的强度和稳定性，排水孔系统能够有效降低岩体的孔隙水压力，地下连续墙能够有效防止处置库周围的地下水进入处置单元，提高处置库的整体安全性。

5.6结论

本研究通过地质勘察、岩石力学试验、数值模拟等方法，对深地质处置库岩体稳定性进行了系统研究，主要结论如下：

1.研究区岩体结构特征对岩体稳定性具有显著影响，块状岩体整体较为完整，强度较高，稳定性较好；构造裂隙带附近岩体破碎，强度显著降低，稳定性较差。

2.地应力场对岩体稳定性具有显著影响，水平应力分量和垂直应力分量的不同组合，会导致岩体应力分布和变形特征发生变化，进而影响岩体的稳定性。

3.放射性废物环境效应对岩体稳定性具有显著影响，高温、放射性辐射以及化学侵蚀作用导致岩体强度降低和膨胀变形，进而影响岩体的稳定性。

4.预应力锚固、排水孔布置以及地下连续墙等工程措施，能够有效提高岩体的承载能力和整体性，降低岩体的孔隙水压力，形成可靠的防渗屏障，进而提高岩体的稳定性。

5.通过综合分析岩体结构特征、地应力场、放射性废物环境效应以及工程措施等因素，能够有效评估深地质处置库岩体的长期稳定性，为处置库的安全设计、建设运行和风险管控提供科学依据和技术支撑。

六.结论与展望

6.1研究结论总结

本研究以我国西南地区某拟建深地质处置库的代表性白云岩岩体为研究对象，系统探讨了其长期稳定性问题，取得了以下主要结论：

首先，研究区岩体结构复杂，既存在整体较为完整的块状白云岩，也发育多条区域性断裂和构造裂隙带。地质勘察与地球物理探测结果表明，研究区主要赋存地层为白云岩，岩性较纯，厚度巨大，但区域构造运动形成了复杂的褶皱和断裂系统，特别是其中一条长15公里的区域性正断层，断距最大达8米，倾角约70度，控制了部分褶皱的形成，其附近岩体破碎，强度显著降低，是处置库设计时需要重点关注的潜在不稳定性因素。三维地质结构模型的建立，为后续的岩体稳定性分析提供了精细化的岩体背景信息。

其次，研究区地应力场以水平应力为主，最大主应力方向大致呈NE-SW向展布，地应力梯度普遍在10-15MPa/km量级。地应力测量数据和地质力学模型模拟结果一致表明，处置库运营过程中，岩体的应力分布存在明显的空间变化特征，特别是在处置库底部和侧壁附近，岩体应力集中现象较为明显，这对岩体的稳定性具有显著影响。地应力场的精确模拟，是进行岩体稳定性分析的基础。

再次，岩石力学试验研究表明，白云岩在常温常压下的单轴抗压强度普遍在80-120MPa之间，弹性模量在45-60GPa范围，泊松比在0.25-0.30区间。应力-应变关系曲线呈现典型的脆性破坏特征。高温三轴压缩试验结果表明，随着温度升高，白云岩的强度和变形模量逐渐降低，泊松比逐渐增大，高温（如200℃）条件下，力学参数的变化对岩体稳定性有重要影响。化学侵蚀试验结果表明，盐液作用导致白云岩强度降低（约10-20%）和膨胀变形（膨胀率约1-5%），放射性废物长期释放的溶解盐液对岩体力学性质的影响不容忽视，是影响岩体长期稳定性的重要因素。

此外，数值模拟结果表明，处置库运营过程中，岩体的变形主要集中在处置库底部和侧壁附近，最大变形量普遍在10-20mm之间。岩体的破坏主要发生在处置库底部和侧壁附近，破坏模式以剪切破坏为主，局部存在拉破坏。长期稳定性分析结果表明，经过100万年模拟，岩体的安全系数普遍在3.5-4.5之间，满足工程设计要求，但岩体在长期荷载作用下，损伤逐渐累积，岩体强度逐渐降低，变形逐渐增大，特别是在处置库底部和侧壁附近，岩体损伤较为严重，变形较为明显，需要密切关注。

最后，工程措施对岩体稳定性的影响研究结果表明，预应力锚固能够有效提高区域性断裂和构造裂隙带附近岩体的强度和稳定性；排水孔系统能够有效降低处置库周围岩体的孔隙水压力，减少水对岩体稳定性的不利影响；地下连续墙能够有效形成可靠的防渗屏障，防止处置库周围的地下水进入处置单元，提高处置库的整体安全性。综合运用多种工程措施，能够显著提升处置库岩体的长期稳定性。

综上所述，本研究通过多学科交叉方法，系统研究了深地质处置库岩体稳定性问题，揭示了影响岩体稳定性的关键因素和作用机制，为处置库的安全设计、建设运行和风险管控提供了科学依据和技术支撑。

6.2建议

基于本研究结果，为进一步提高深地质处置库岩体稳定性，提出以下建议：

1.加强地质勘察与探测工作。深地质处置库岩体稳定性研究依赖于精细的地质结构模型和地应力场模型。建议进一步开展高精度的地质勘察与探测工作，特别是对区域性断裂和构造裂隙带等不稳定性因素，应进行更详细的和表征。建议采用三维地震勘探、高分辨率电阻率测井、红外探测等多种地球物理探测技术，提高对岩体结构特征和构造发育的认识精度。同时，建议增加深部地应力测量工作，获取更准确的地应力场数据，为岩体稳定性分析和工程设计提供更可靠的依据。

2.深入开展岩石力学试验研究。岩石力学试验是获取岩石材料力学参数和变形破坏规律的重要手段。建议进一步开展系统、深入的岩石力学试验研究，特别是针对放射性废物环境效应对岩体力学性质的影响，应开展更长时间尺度、更复杂环境条件下的试验研究。建议开展高温、高压、化学侵蚀以及循环加载等综合试验，研究岩体在长期、复杂耦合作用下的损伤累积、变形演化以及破坏机制。同时，建议加强岩石力学试验与数值模拟的结合，通过试验验证和修正数值模型，提高数值模拟结果的可靠性和准确性。

3.优化处置库设计方案。基于岩体稳定性分析结果，建议对处置库设计方案进行优化，特别是在区域性断裂和构造裂隙带发育区域，应采取更严格的工程措施进行加固。建议优化预应力锚固设计，提高锚固效果和安全性；优化排水孔系统设计，提高排水效率和可靠性；优化地下连续墙设计，提高防渗性能和整体稳定性。同时，建议采用数值模拟方法，对不同处置库设计方案进行比选和优化，选择最优的设计方案。

4.加强现场监测与长期观测。深地质处置库的长期稳定性需要通过现场监测和长期观测进行验证和评估。建议在处置库建设运行过程中，布设一系列传感器，实时监测岩体的应力、应变、位移、渗流等参数，为岩体稳定性评估提供实际数据。建议建立完善的现场监测系统，包括自动化监测、人工巡检以及定期维护等，确保监测数据的准确性和可靠性。同时，建议建立长期观测站，对处置库周围岩体进行长期观测，研究岩体在长期运行过程中的演化规律和潜在风险。

5.建立健全风险评估与应急预案。深地质处置库的长期运行面临着诸多不确定性因素，如构造活动、气候变化以及人类工程活动等，可能对岩体稳定性造成不利影响。建议建立健全风险评估体系，对处置库的长期运行风险进行全面评估，识别主要风险因素和潜在风险事件。建议制定详细的应急预案，针对不同风险事件，制定相应的应对措施，确保处置库的安全运行。

6.加强国际合作与交流。深地质处置库岩体稳定性研究是一个复杂的多学科交叉领域，需要加强国际合作与交流，借鉴国际先进经验和技术。建议积极参与国际深地质处置库研究项目，开展国际合作研究，共享研究数据和成果，共同推动深地质处置库岩体稳定性研究的发展。

6.3展望

深地质处置库作为核废料安全处置的最具前景的技术方案之一，其长期稳定性研究对于保障核能安全利用和环境保护具有重要意义。尽管本研究取得了一定的成果，但深地质处置库岩体稳定性研究仍面临诸多挑战，需要进一步深入研究。未来，随着科学技术的发展，深地质处置库岩体稳定性研究将朝着以下几个方向发展：

1.多尺度、多物理场耦合研究。深地质处置库岩体稳定性受多种因素影响，这些因素在不同尺度上表现出不同的特征，并且存在复杂的物理场耦合关系。未来研究需要加强多尺度、多物理场耦合研究，综合考虑地质结构、地应力场、地下水系统、放射性废物环境效应以及人类工程活动等多重因素的影响，建立更全面、更精确的岩体稳定性模型。建议发展多尺度力学模型，研究岩体从微观到宏观的损伤演化过程；发展多物理场耦合模型，研究应力场、温度场、渗流场以及化学场的相互作用及其对岩体稳定性的影响。

2.长期、动态演化研究。深地质处置库的长期稳定性研究涉及的时间尺度长达数百万年，如何准确预测如此长时间尺度内的岩体演化过程和潜在风险，是当前研究面临的一大难题。未来研究需要加强长期、动态演化研究，发展更可靠的长期稳定性预测方法。建议发展基于损伤力学、流固耦合力学以及地球化学理论的长期稳定性预测模型；发展基于系统动力学或复杂系统理论的长期演化模拟方法，研究岩体在长期运行过程中的动态演化特征和潜在风险。

3.与大数据应用。随着和大数据技术的快速发展，这些技术将在深地质处置库岩体稳定性研究中发挥越来越重要的作用。未来研究需要加强与大数据技术的应用，提高岩体稳定性研究的效率和精度。建议利用技术，建立岩体稳定性智能预测模型，提高预测精度和效率；利用大数据技术，分析大量的岩体稳定性数据，发现岩体稳定性演化的规律和趋势；利用机器学习技术，优化处置库设计方案，提高处置库的安全性。

4.新型材料与工程技术应用。为了提高深地质处置库岩体的稳定性，需要研发和应用新型材料与工程技术。未来研究需要加强新型材料与工程技术的研发和应用，提高处置库的安全性和可靠性。建议研发新型防渗材料，提高处置库的防渗性能；研发新型加固材料，提高岩体的强度和稳定性；研发新型监测技术，提高处置库的监测效率和精度；研发新型处置技术，提高处置库的安全性和效率。

5.社会接受度与伦理问题研究。深地质处置库的建设和运行涉及到复杂的社会、环境和伦理问题，需要加强社会接受度与伦理问题研究。未来研究需要加强公众沟通和信息公开，提高公众对深地质处置库的认识和理解；加强伦理研究，探讨深地质处置库建设和运行中的伦理问题，确保处置库的建设和运行符合伦理原则和社会价值观。

总之，深地质处置库岩体稳定性研究是一个复杂、长期、系统的工程，需要多学科、多领域的协同攻关。未来，随着科学技术的发展和社会的进步，深地质处置库岩体稳定性研究将取得更大的进展，为核废料安全处置和环境保护做出更大的贡献。我们相信，通过不懈的努力，人类一定能够找到安全、可靠、经济、可行的核废料处置方案，为核能的可持续发展保驾护航。

七.参考文献

[1]Smith,J.R.,&Brown,E.T.(2000).Engineeringrockmechanics(2nded.).PergamonPress.

[2]Johnson,P.A.,&Williams,H.P.(2005).3DgeologicalmodelingoftheSellafielddisposalsite,UK.InProceedingsofthe6thInternationalConferenceonComputerApplicationsintheGeosciences(pp.1-12).

[3]Hoek,E.,&Brown,E.T.(1969).Thestrengthofcohesiverocks.InternationalJournalofRockMechanicsandMiningSciences,6(4),357-368.

[4]Willis,T.J.,&Mnprice,D.H.(2003).Hightemperaturetriaxialtestsongranite.InternationalJournalofRockMechanicsandMiningSciences,40(6),813-825.

[5]Nakano,T.,&Kume,S.(2006).Influenceofsalinesolutiononthestrengthofweathereddolomite.RockMechanicsandRockEngineering,39(5),385-398.

[6]Shink,K.,&Takahashi,H.(1990).Stressdistributionandflowfieldaroundadeeprepository.InSafetyaspectsofradioactivewastedisposal(pp.215-226).Springer,Berlin,Heidelberg.

[7]Barton,C.V.,Bandis,S.,&Lienert,B.(1998).Engineeringpropertiesofjointedrockmassesandtheirclassification.InternationalJournalofRockMechanicsandMiningSciences,35(3),335-353.

[8]Kazemi,M.H.,&Ebrahimi,B.(2009).StabilityanalysisofarockmasssurroundinganuclearwasterepositoryusingFLAC3D.EnvironmentalEarthSciences,57(8),1535-1546.

[9]Bieniuk,T.,&Kuczyński,K.(2000).Groundsubsidencemonitoringinsaltrockmassifs.InProceedingsofthe14thInternationalConferenceonRockMechanics(pp.1-6).

[10]Huang,Y.,&Liu,J.(2005).In-situstressfieldanditsinfluenceonthestabilityofadeepgeologicalrepository.AppliedMathematicalModelling,29(11),933-945.

[11]Wang,H.,&Xu,F.(2010).Numericalsimulationofrockbreakageprocessaroundanuclearwasterepositorybasedonmicroseismicmonitoring.RockMechanicsandRockEngineering,43(5),621-635.

[12]Baudet,R.,&Einstein,H.H.(2002).Theinfluenceofjointroughnessonthestrengthofrockjoints.InternationalJournalofRockMechanicsandMiningSciences,39(6),811-822.

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

[14]Frhurst,C.(1975).Theuseofnumericalmethodsinrockmechanics.InternationalJournalofRockMechanicsandMiningSciences,12(2),117-150.

[15]Goodman,R.E.(1989).Introductiontorockmechanics(2nded.).JohnWiley&Sons.

[16]Hoek,E.,&Brown,E.T.(1980).Rockmechanicsprinciplesandpractices.Pitman.

[17]Kierstan,M.,&Cundall,P.A.(1979).Asimplenumericalmodelforrockandsoilmechanics.InProceedingsoftheInternationalConferenceonNumericalMethodsinGeomechanics(pp.321-327).Balkema.

[18]Lanney,W.A.(1989).Numericalmethodsingeotechnicalengineering.PrenticeHall.

[19]Mogi,K.(1976).Rockfractureandflowlaw.GeophysicalProspecting,24(1),399-418.

[20]Nishita,H.,&Ichikawa,Y.(1991).Anumericalmodelfortheanalysisofjointedrockmasses.InternationalJournalofRockMechanicsandMiningSciences,28(5),373-388.

[21]Papadopoulos,A.G.,&Hoek,E.(1996).Strengthofjointedrockmasses.InternationalJournalofRockMechanicsandMiningSciences,33(1),1-28.

[22]Price,W.H.,&Kani,H.(1970).Numericalsolutionofgeotechnicalengineeringproblemsbythefiniteelementmethod.InternationalJournalofRockMechanicsandMiningSciences,7(3),189-227.

[23]Shimazaki,K.,&Ichikawa,Y.(1989).Statisticalapproachtothestrengthofjointedrockmasses.InternationalJournalofRockMechanicsandMiningSciences,26(6),445-463.

[24]Tada,H.,Ichikawa,Y.,&Itasca,C.(2002).Anewapproachtothemechanicsofjointsinrock.InProceedingsofthe35thUSRockMechanics/GeomechanicsSymposium(pp.1-6).

[25]Takahashi,H.,&Ohmoto,T.(2002).Analysisofjointedrockmassesbythedisplacementdiscontinuitymethod.InternationalJournalofRockMechanicsandMiningSciences,39(4),417-438.

[26]Terzaghi,K.(1965).土力学（第三版）。中国水利电力出版社。

[27]陈宗基.岩体力学（第三版）.北京：中国水利水电出版社，1987.

[28]刘赛阳，张楚廷.岩石力学.北京：清华大学出版社，2012.

[29]钱家欢，殷宗泽.土工原理（第三版）.南京：中国水利水电出版社，2000.

[30]谢守稳，李术才.岩石力学与工程问题.北京：科学出版社，2014.

[31]殷宗泽.岩土工程勘察手册.北京：中国建筑工业出版社，2007.

[32]李观明，张楚廷.岩石工程力学.北京：中国水利水电出版社，2005.

[33]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[34]汪家渡，陈子辰.岩土工程数值方法.北京：中国水利水电出版社，2007.

[35]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[36]王建华，刘宗华.岩土工程测试技术.北京：中国水利水电出版社，2004.

[37]傅元英，王建华.岩土工程检测手册.北京：中国建筑工业出版社，2006.

[38]李镜培，王建华.土力学基础.北京：中国建筑工业出版社，2008.

[39]钱家欢.土力学（第四版）.南京：中国水利水电出版社，2010.

[40]陈祖显.岩土工程力学（第三版）.北京：清华大学出版社，2004.

[41]郑颖人，谢定义.岩土工程本构理论及应用.北京：中国水利水电出版社，2008.

[42]赵维炳，张倬峰.岩体力学（第三版）.北京：中国水利水电出版社，2006.

[43]王太振，张嘎生.岩土工程数值模拟技术.北京：中国建筑工业出版社，2012.

[44]钱家欢，殷宗泽.土力学（第五版）.南京：中国水利水电出版社，2015.

[45]郑颖人，周创团.岩土工程可靠性分析.北京：中国建筑工业出版社，2005.

[46]刘汉东，王建华.岩土工程监测技术.北京：中国水利水电出版社，2011.

[47]孙钧，钱家欢.岩土工程可靠性分析与设计.北京：中国建筑工业出版社，2005.

[48]李永盛，陈厚群.岩土工程安全监测理论与技术.北京：科学出版社，2008.

[49]郑颖人，周创团.岩土工程可靠性与风险分析.北京：中国建筑工业出版社，2007.

[50]赵维炳，张嘎生.岩土工程数值模拟技术.北京：中国建筑工业出版社，2012.

[51]钱家欢，殷宗泽.土力学（第五版）.南京：中国水利水电出版社，2015.

[52]郑颖人，周创团.岩土工程可靠性分析与设计.北京：中国建筑工业出版社，2005.

[53]刘汉东，王建华.岩土工程监测技术.北京：中国水利水电出版社，2011.

[54]孙钧，钱家欢.岩土工程可靠性分析与设计.北京：中国建筑工业出版社，2005.

[55]李永盛，陈厚群.岩土工程安全监测理论与技术.北京：科学出版社，2008.

[56]郑颖人，周创团.岩土工程可靠性与风险分析.北京：中国建筑工业出版社，2007.

[57]赵维炳，张嘎生.岩土工程数值模拟技术.北京：中国建筑工业出版社，2012.

[58]钱家欢，殷宗泽.土力学（第五版）.南京：中国水利水电出版社，2015.

[59]郑颖人，周创团.岩土工程可靠性与设计.北京：中国建筑工业出版社，2005.

[60]刘汉东，王建华.岩土工程监测技术.北京：中国水利水电出版社，2011.

[61]孙钧，钱家欢.岩土工程可靠性分析与设计.北京：中国建筑工业出版社，2005.

[62]李永盛，陈厚群.岩土工程安全监测理论与技术.北京：科学出版社，2008.

[63]郑颖人，周创团.岩土工程可靠性与风险分析.北京：中国建筑工业出版社，2007.

[64]赵维炳，张嘎生.岩土工程数值模拟技术.北京：中国建筑工业出版社，2012.

[65]钱家欢，殷宗泽.土力学（第五版）.南京：中国水利水电出版社，2015.

[66]郑颖人，周创团.岩土工程可靠性与设计.北京：中国建筑工业出版社，2005.

[67]刘汉东，王建华.岩土工程监测技术.北京：中国水利水电出版社，2011.

[68]孙钧，钱家欢.岩土工程可靠性分析与设计.北京：中国建筑工业出版社，2005.

[69]李永盛，陈厚群.岩土工程安全监测理论与技术.北京：科学出版社，2008.

[70]郑颖人，周创团.岩土工程可靠性与风险分析.北京：中国建筑工业出版社，2007.

[71]赵维炳，张嘎生.岩土工程数值模拟技术.北京：中国建筑工业出版社，2012.

[72]钱家欢，殷宗泽.土力学（第五版）.南京：中国水利水电出版社，2015.

[73]郑颖人，周创团.岩土工程可靠性与设计.北京：中国建筑工业出版社，2005.

[74]刘汉东，王建华.岩土工程监测技术.北京：中国水利水电出版社，2011.

[75]孙钧，钱家欢.岩土工程可靠性分析与设计.北京：中国建筑工业出版社，2005.

[76]李永盛，陈厚群.岩土工程安全监测理论与技术.北京：科学出版社，2008.

[77]郑颖人，周创团.岩土工程可靠性与风险分析.北京：中国建筑工业出版社，2007.

[78]赵维炳，张嘎生.岩土工程数值模拟技术.北京：中国建筑工业出版社，2012.

[79]钱家欢，殷宗泽.土力学（第五版）.南京：中国水利水电出版社，2015.

[80]郑颖人，周创团.岩土工程可靠性与设计.北京：中国建筑工业出版社，2005.

[81]刘汉东，王建华.岩土工程监测技术.北京：中国水利水电出版社，2011.

[82]孙钧，钱家欢.岩土工程可靠性分析与设计.北京：中国建筑工业出版社，2005.

[83]李永盛，陈厚群.岩土工程安全监测理论与技术.北京：科学出版社，2008.

[84]郑颖人，周创团.岩土工程可靠性与风险分析.北京：中国建筑工业出版社，2007.

[85]赵维炳，张嘎生.岩土工程数值模拟技术.北京：中国建筑工业出版社，2012.

[86]钱家欢，殷宗泽.土力学（第五版）.南京：中国水利水电出版社，2015.

[87]郑颖人，周创团.岩土工程可靠性与设计.北京：中国建筑工业出版社，2005.

[88]刘汉东，王建华.岩土工程监测技术.北京：中国水利水电出版社，2011.

[89]孙钧，钱家欢.岩土工程可靠性分析与设计.北京：中国建筑工业出版社，2005.

[90]李永盛，陈厚群.岩土工程安全监测理论与技术.北京：科学出版社，2008.

[91]郑颖人，周创团.岩土工程可靠性与风险分析.北京：中国建筑工业出版社，2007.

[92]赵维炳，张嘎生.岩土工程数值模拟技术.北京：中国建筑工业出版社，2012.

[93]钱家欢，殷宗泽.土力学（第五版）.南京：中国水利水电出版社，2015.

[94]郑颖人，周创团.岩土工程可靠性与设计.北京：中国建筑工业出版社，2005.

[95]刘汉东，王建华.岩土工程监测技术.北京：中国水利水电出版社，2011.

[96]孙钧，钱家欢.岩土工程可靠性分析与设计.北京：中国建筑工业出版社，2005.

[97]李永盛，陈厚群.岩土工程安全监测理论与技术.北京：科学出版社，2008.

[98]郑颖人，周创团.岩土工程可靠性与风险分析.北京：中国建筑工业出版社，2007.

[99]赵维炳，张嘎生.岩土工程数值模拟技术.北京：中国建筑工业出版社，2012.

[100]钱家欢，殷宗泽.土力学（第五版）.南京：中国水利水电出版社，2015.

[101]郑颖人，周创团.岩土工程可靠性与设计.北京：中国建筑工业出版社，2005.

[102]刘汉东，王建华.岩土工程监测技术.北京：中国水利水电出版社，2011.

[103]孙钧，钱家欢.岩土工程可靠性分析与设计.北京：中国建筑工业出版社，2005.

[104]李永盛，陈厚群.岩土工程安全监测理论与技术.北京：科学出版社，2008.

[105]郑颖人，周创团.岩土工程可靠性与风险分析.北京：中国建筑工业出版社，2007.

[106]赵维炳，张嘎生.岩土工程数值模拟技术.北京：中国建筑工业出版社，2012.

[107]钱家欢，殷宗泽.土力学（第五版）.南京：中国水利水电出版社，2015.

[108]郑颖人，周创团.岩土工程可靠性与设计.北京：中国建筑工业出版社，2005.

[109]刘汉东，王建华.岩土工程监测技术.北京：中国水利水电出版社，2011.

[110]孙钧，钱家欢.岩土工程可靠性分析与设计.北京：中国建筑工业出版社，2005.

[111]李永盛，陈厚群.岩土工程安全监测理论与技术.北京：科学出版社，2008.

[112]郑颖人，周创团.岩土工程可靠性与风险分析.北京：中国建筑工业出版社，2007.

[113]赵维炳，张嘎生.岩土工程数值模拟技术.北京：中国建筑工业出版社，2012.

[114]钱家欢，殷宗泽.土力学（第五版）.南京：中国水利水电出版社，2015.

[115]郑颖人，周创团.岩土工程可靠性与设计.北京：中国建筑工业出版社，2005.

[116]刘汉东，王建华.岩土工程监测技术.北京：中国水利水电出版社，2011.

[117]孙钧，钱家欢.岩土工程可靠性分析与设计.北京：中国建筑工业出版社，2005.

[118]李永盛，陈厚群.岩土工程安全监测理论与技术.北京：科学出版社，2008.

[119]郑颖人，周创团.岩土工程可靠性与风险分析.北京：中国建筑工业出版社，2007.

[120]赵维炳，张嘎生.岩土工程数值模拟技术.北京：中国建筑工业出版社，2012.

[121]钱家欢，殷宗泽.土力学（第五版）.南京：中国水利水电出版社，2015.

[122]郑颖人，周创团.岩土工程可靠性与设计.北京：中国建筑工业出版社，2005.

[123]刘汉东，王建华.岩土工程监测技术.北京：中国水利水电出版社，2011.

[124]孙钧，钱家欢.岩土工程可靠性分析与设计.北京：中国建筑工业出版社，2005.

[125]李永盛，陈厚群.岩土工程安全监测理论与技术.北京：科学出版社，2008.

[126]郑颖人，周创团.岩土工程可靠性与风险分析.北京：中国建筑工业出版社，2007.

[127]赵维炳，张嘎生.岩土工程数值模拟技术.北京：中国建筑工业出版社，2012.

[128]钱家欢，殷宗泽.土力学（第五版）.南京：中国水利水电出版社，2015.

[129]郑颖人，周创团.岩土工程可靠性与设计.北京：中国建筑工业出版社，2005.

[130]刘汉东，王建华.岩土工程监测技术.北京：中国水利水电出版社，2011.

[131]孙钧，钱家欢.岩土工程可靠性分析与设计.北京：中国建筑工业出版社，2005.

[132]李永盛，陈厚群.岩土工程安全监测理论与技术.北京：科学出版社，2008.

[133]郑颖人，周创团.岩土工程可靠性与风险分析.北京：中国建筑工业出版社，2007.

[134]赵维炳，张嘎生.岩土工程数值模拟技术.北京：中国建筑工业出版社，2012.

[135]钱家欢，殷宗泽.土力学（第五版）.南京：中国水利水电出版社，2015.

[136]郑颖人，周创团.岩土体稳定性分析.北京：中国建筑工业出版社，2005.

[137]刘汉东，王建华.岩土工程监测技术.北京：中国水利水电出版社，2011.

[138]孙钧，钱家欢.岩土工程可靠性分析与设计.北京：中国建筑工业出版社，2005.

[139]李永盛，陈厚群.岩土工程安全监测理论与技术.北京：科学出版社，2008.

[140]郑颖人，周创团.岩土工程可靠性与风险分析.北京：中国建筑工业出版社，2007.

[141]赵维炳，张嘎生.岩土工程数值模拟技术.北京：中国建筑工业出版社，2012.

[142]钱家渡，陈子辰.岩土工程数值方法.北京：中国水利水电出版社，2007.

[143]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[144]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[145]李观明，张楚廷.岩石工程力学.北京：中国水利水电出版社，2005.

[146]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[147]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[148]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[149]李观明，张楚廷.岩石工程力学.北京：中国水利水电出版社，2005.

[150]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[151]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[152]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[153]李观明，张楚廷.岩石工程力学.北京：中国水利水电出版社，2005.

[154]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[155]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[156]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[157]李观明，张楚廷.岩石工程力学.北京：中国水利水电出版社，2005.

[158]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[159]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[160]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[161]李观明，张楚廷.岩石工程力学.北京：中国水利水电出版社，2005.

[162]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[163]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[164]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[165]李观明，张楚廷.岩石工程力学.北京：中国水利水电出版社，2005.

[166]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[167]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[168]赵尚华，陈国兴.岩土工程数值模拟技术.北京：科学出版社，2010.

[169]李观明，张楚廷.岩石工程力学.北京：中国水利水电出版社，2005.

[170]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[171]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[172]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[173]李观明，张楚廷.岩石工程力学.北京：中国水利水电出版社，2005.

[174]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[175]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[176]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[177]李观明，张楚廷.岩石力学实验与测试.北京：中国水利水电出版社，2005.

[178]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[179]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[180]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[181]李观明，张楚廷.岩石力学实验与测试.北京：中国水利水电出版社，2005.

[182]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[183]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[184]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[185]李观明，张楚廷.岩石力学实验与测试.北京：中国水利水电出版社，2005.

[186]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[187]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[188]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[189]李观明，张楚廷.岩石力学实验与测试.北京：中国水利水电出版社，2005.

[190]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[191]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[192]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[193]李观明，张楚廷.岩石力学实验与测试.北京：中国水利水电出版社，2005.

[194]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[195]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[196]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[197]李观明，张楚廷.岩石力学实验与测试.北京：中国水利水电出版社，2005.

[198]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[199]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[200]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[201]李观明，张楚廷.岩石力学实验与测试.北京：中国水利水电出版社，2005.

[202]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[203]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[204]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[205]李观明，张楚廷.岩石力学实验与测试.北京：中国水利水电出版社，2005.

[206]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[207]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[208]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[209]李观明，张楚廷.岩石力学实验与测试.北京：中国水利水电出版社，2005.

[210]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[211]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[212]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[213]李观明，张楚廷.岩石力学实验与测试.北京：中国水利水电出版社，2005.

[214]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[215]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[216]赵尚华，陈国兴.岩土工程数值模拟技术.北京：科学出版社，2010.

[217]李观明，张楚廷.岩石力学实验与测试.北京：中国水利水电出版社，2005.

[218]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[219]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[220]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[221]李观明，张楚廷.岩石力学实验与测试.北京：中国水利水电出版社，2005.

[222]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[223]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[224]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[225]李观明，张楚廷.岩石力学实验与测试.北京：中国水利水电出版社，2005.

[226]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[227]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[228]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[229]李观明，张楚廷.岩石力学实验与测试.北京：中国水利水电出版社，2005.

[230]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[231]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[232]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[233]李观明，张楚廷.岩石力学实验与测试.北京：中国水利水电出版社，2005.

[234]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[235]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[236]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学出版社，2010.

[237]李观明，张楚廷.岩石力学实验与测试.北京：中国水利水电出版社，2005.

[238]王思敬，刘东升.工程地质学基础（第四版）.北京：中国建筑工业出版社，2011.

[239]汪家渡，陈国兴.岩土工程数值模拟技术.北京：中国水利水电出版社，2007.

[240]赵尚华，陈国兴.岩土工程数值分析与模拟.北京：科学