深埋地下洞室施工期硬岩应力型灾害风险动态评估方法：理论、模型与实践

深埋地下洞室施工期硬岩应力型灾害风险动态评估方法：理论、模型与实践一、引言1.1研究背景与意义随着全球基础设施建设的快速推进，深埋地下洞室工程在交通、水利、能源等领域得到了广泛应用。例如，在铁路和公路建设中，为了穿越山脉、河流等复杂地形，常常需要修建深埋隧道；在水利水电工程中，深埋输水隧洞和地下厂房是实现水资源合理利用和水能开发的关键设施；在能源领域，深埋地下储气库和石油储备库对于保障能源安全具有重要意义。然而，深埋地下洞室在施工过程中面临着诸多挑战，其中硬岩应力型灾害是最为突出的问题之一。硬岩应力型灾害主要包括岩爆、板裂破坏等，这些灾害的发生往往具有突发性和强烈的破坏性。当岩体中的应力超过其自身的强度极限时，就会引发岩爆，导致岩石碎片从岩体中突然弹射出来，对施工人员和设备造成严重威胁。板裂破坏则表现为岩体在高应力作用下，沿平行于洞壁的方向发生分层破裂，形成板状剥落，这不仅会影响洞室的稳定性，还可能引发坍塌事故。这些灾害不仅会对施工人员的生命安全构成严重威胁，导致人员伤亡，还会对施工设备造成巨大的破坏，增加工程成本。灾害发生后，修复受损的洞室和设备需要耗费大量的人力、物力和时间，从而导致工程进度延误，给工程带来巨大的经济损失。此外，灾害还可能引发次生地质灾害，进一步加剧工程的安全风险。在锦屏二级水电站引水隧洞施工中，由于地应力极高，岩体完整性好，岩爆灾害频繁发生。强烈的岩爆导致岩石碎片高速弹射，对施工人员的生命安全造成了极大威胁，同时也严重损坏了施工设备，如TBM掘进机的刀盘、护盾等部件，导致施工进度严重滞后，工程成本大幅增加。又如，在挪威的一些深埋公路隧道施工中，也曾发生过严重的岩爆灾害，造成了人员伤亡和工程延误，给当地的交通建设带来了巨大的困难。因此，开展深埋地下洞室施工期硬岩应力型灾害风险动态评估方法的研究具有重要的现实意义。准确评估硬岩应力型灾害风险，能够为工程设计和施工提供科学依据，帮助工程师制定合理的施工方案和防护措施，有效降低灾害发生的概率和危害程度，保障施工人员的生命安全和工程的顺利进行，具有显著的经济和社会效益。1.2国内外研究现状硬岩应力型灾害风险评估一直是岩土工程领域的研究热点，国内外学者在这方面开展了大量研究工作。在岩爆风险评估方面，早期的研究主要集中在单指标经验判据法，如以洞周最大切向应力为主要因素的hoek判据、russenes判据、turchaninov判据和二郎山隧洞判据，以及仅以原岩应力场的最大主应力为特征参量的barton判据、陶振宇判据和国标gb50218-94判据。这些判据简单直观，易于应用，但普遍存在局限性，它们大多采用单个应力分量或两个应力分量与岩石强度建立联系，未充分考虑深部地下工程围岩处于高围压状态下，三个方向主应力对围岩能量积聚的综合影响。随着研究的深入，多指标集结法逐渐被提出，通过综合考虑多个影响因素，如岩石强度、地应力、岩体结构等，来提高评估的准确性。例如，有学者将岩石的弹性模量、泊松比、单轴抗压强度等地力学参数与地应力相结合，建立了岩爆风险评估指标体系。不确定性分析法也得到了应用，该方法考虑了评估过程中各种因素的不确定性，如地应力测量误差、岩石力学参数的离散性等，采用概率统计、模糊数学等方法对岩爆风险进行评估。文献[X]利用模糊综合评判法，对深埋隧道岩爆风险进行评估，通过建立模糊关系矩阵和权重向量，综合考虑多个因素的影响，得出岩爆风险等级。综合排序法通过对多个评估指标进行排序和加权，确定岩爆风险的相对大小。机器学习法近年来在岩爆风险评估中得到了广泛应用，如神经网络、支持向量机等。神经网络具有强大的非线性映射能力，能够自动学习输入数据与岩爆风险之间的复杂关系。支持向量机则在小样本、非线性问题上表现出良好的性能。有学者利用基于粗糙集理论的遗传-RBF神经网络对岩爆进行预测，通过粗糙集理论对数据进行约简，去除冗余信息，提高了神经网络的训练效率和预测精度。数值模拟法通过建立岩体的力学模型，模拟隧道开挖过程中围岩的应力、应变分布，预测岩爆发生的可能性。常用的数值模拟方法包括有限元法、离散元法等。在板裂破坏风险评估方面，研究主要围绕板裂破坏的机理、判据和影响因素展开。通过对现场工程实例的观察和分析，发现硬岩的板裂、片帮、V形槽、岩芯饼化等张性破坏特征。理论研究方面，基于弹性力学、断裂力学等理论，建立了板裂破坏的力学模型。有学者根据弹性薄板理论，分析了洞室围岩在高应力作用下的板裂破坏机理，得出了板裂破坏的判据公式。试验研究则通过室内模型试验，模拟板裂破坏过程，研究其破坏特征和影响因素。尽管国内外在硬岩应力型灾害风险评估方面取得了一定的成果，但仍存在一些不足之处。现有评估方法大多基于特定的工程背景和数据，缺乏普适性，难以直接应用于不同地质条件和工程类型的深埋地下洞室。对灾害发生的微观机理研究还不够深入，导致评估模型的物理基础不够坚实。在实际工程中，地应力、岩石力学参数等数据的获取存在一定困难，且准确性难以保证，这也影响了风险评估的精度。目前的研究主要侧重于灾害发生的可能性评估，对灾害发生后的危害程度评估相对较少，不利于制定全面的风险应对策略。1.3研究内容与方法1.3.1研究内容硬岩应力型灾害机理深入剖析：系统分析岩爆和板裂破坏的发生机理。对于岩爆，从能量角度深入研究岩体在开挖卸荷过程中，弹性应变能的积聚、释放机制，以及能量释放对岩体破坏的影响。结合岩石力学理论，探讨地应力、岩石强度、岩体结构等因素在岩爆发生过程中的相互作用关系。针对板裂破坏，基于弹性力学和断裂力学理论，研究洞室围岩在高应力作用下，板裂破坏的力学机制，分析板裂破坏的发生条件和发展过程，以及其与岩石微观结构的关联。风险评估指标体系构建：全面考虑影响硬岩应力型灾害风险的各种因素，构建科学合理的风险评估指标体系。这些因素包括地应力，详细分析其大小、方向和分布特征对灾害风险的影响；岩石力学参数，如弹性模量、泊松比、单轴抗压强度等，研究它们在评估中的作用；岩体结构，包括节理、裂隙的发育程度和分布规律对灾害风险的影响；地下水，分析其对岩体强度和稳定性的影响；工程施工因素，如开挖方法、支护措施等对灾害风险的影响。通过对这些因素的综合考虑，确保评估指标体系能够全面、准确地反映灾害风险。动态评估模型建立：引入人工智能算法，如神经网络、支持向量机等，建立深埋地下洞室施工期硬岩应力型灾害风险动态评估模型。利用工程现场监测数据和室内试验数据对模型进行训练和验证，不断优化模型的参数和结构，提高模型的准确性和可靠性。模型应能够根据施工过程中实时获取的监测数据，动态更新风险评估结果，及时反映灾害风险的变化情况。同时，考虑模型的可解释性，以便工程人员能够理解模型的评估过程和结果，为工程决策提供有力支持。工程应用与验证：将建立的风险动态评估模型应用于实际深埋地下洞室工程，如锦屏二级水电站引水隧洞、挪威某深埋公路隧道等。结合工程现场的地质条件、施工情况和监测数据，对模型的评估结果进行验证和分析。通过实际应用，检验模型的有效性和实用性，及时发现模型存在的问题和不足，并进行改进和完善。同时，根据评估结果，为工程施工提供科学合理的建议，如优化施工方案、加强支护措施等，降低灾害风险，保障工程的安全顺利进行。1.3.2研究方法理论分析：运用岩石力学、弹性力学、断裂力学等相关理论，深入分析硬岩应力型灾害的发生机理，为风险评估提供坚实的理论基础。通过理论推导，建立灾害发生的力学模型，分析各种因素对灾害发生的影响机制，为后续的研究提供理论指导。数值模拟：利用有限元软件ANSYS、ABAQUS等，建立深埋地下洞室的数值模型，模拟施工过程中围岩的应力、应变分布和变形破坏过程，预测硬岩应力型灾害的发生可能性。通过数值模拟，可以直观地观察到岩体在不同工况下的力学响应，分析灾害的发生过程和影响范围，为风险评估和工程设计提供参考依据。室内试验：开展岩石力学试验，如单轴抗压试验、三轴压缩试验、巴西劈裂试验等，获取岩石的力学参数，研究岩石的力学特性和破坏机制。进行岩爆和板裂破坏的模拟试验，通过控制试验条件，观察灾害的发生过程和特征，为理论分析和数值模拟提供验证数据。案例研究：收集国内外多个深埋地下洞室工程的案例，分析其施工过程中硬岩应力型灾害的发生情况、风险评估方法和应对措施，总结经验教训，为本文的研究提供实践参考。通过对实际案例的研究，可以深入了解灾害的发生规律和工程应对方法，为建立更有效的风险评估模型提供依据。1.4研究技术路线本研究技术路线如图1所示，主要包括以下几个关键环节：数据采集与分析：收集工程现场的地质勘察资料，包括地层岩性、地质构造、地下水等信息，为后续研究提供基础地质背景。通过现场原位测试，如钻孔水压致裂法、应力解除法等，获取地应力数据，包括地应力的大小、方向和分布特征。同时，开展岩石力学试验，如单轴抗压试验、三轴压缩试验、巴西劈裂试验等，测定岩石的弹性模量、泊松比、单轴抗压强度、抗拉强度等力学参数。对采集到的数据进行整理和分析，初步了解工程区域的地质条件和岩石力学特性。机理研究与指标体系构建：基于岩石力学、弹性力学、断裂力学等理论，深入分析岩爆和板裂破坏的发生机理。从能量角度研究岩爆过程中弹性应变能的积聚与释放机制，以及地应力、岩石强度、岩体结构等因素的相互作用关系；运用弹性力学和断裂力学理论，探讨板裂破坏的力学机制、发生条件和发展过程。综合考虑地应力、岩石力学参数、岩体结构、地下水、工程施工等因素，构建科学合理的硬岩应力型灾害风险评估指标体系。通过理论分析和实际案例研究，确定各指标的权重，确保指标体系能够全面、准确地反映灾害风险。模型建立与训练：选用合适的人工智能算法，如神经网络、支持向量机等，建立深埋地下洞室施工期硬岩应力型灾害风险动态评估模型。利用前期采集的数据对模型进行训练，通过不断调整模型的参数和结构，使模型能够准确地学习到输入数据与灾害风险之间的复杂关系。采用交叉验证等方法对训练好的模型进行验证，评估模型的准确性和可靠性。根据验证结果，对模型进行优化和改进，提高模型的性能。工程应用与验证：将建立的风险动态评估模型应用于实际深埋地下洞室工程，如锦屏二级水电站引水隧洞、挪威某深埋公路隧道等。结合工程现场的地质条件、施工情况和实时监测数据，利用模型对硬岩应力型灾害风险进行动态评估。将评估结果与实际发生的灾害情况进行对比分析，验证模型的有效性和实用性。根据验证结果，对模型进行进一步的完善和优化，为工程施工提供更加准确、可靠的风险评估结果，指导工程决策和施工安全管理。通过以上技术路线，本研究旨在建立一套科学、有效的深埋地下洞室施工期硬岩应力型灾害风险动态评估方法，为深埋地下洞室工程的安全施工提供有力的技术支持。\\二、深埋地下洞室施工期硬岩应力型灾害概述2.1灾害类型及表现形式2.1.1岩爆岩爆，也称冲击地压，是一种岩体中聚积的弹性变形势能在一定条件下突然猛烈释放，导致岩石爆裂并弹射出来的现象。它是深埋地下工程在施工过程中常见的动力破坏现象，当岩体中聚积的高弹性应变能大于岩石破坏所消耗的能量时，破坏了岩体结构的平衡，多余的能量导致岩石爆裂，使岩石碎片从岩体中剥离、崩出。岩爆通常具有以下特征：突发性：在未发生前，并无明显的预兆，甚至可能听不到空响声，一般认为不会掉落石块的地方，也会突然发生岩石爆裂声响，石块有时应声而下，有时暂不坠下。例如，在锦屏二级水电站引水隧洞施工中，部分地段在开挖过程中突然发生岩爆，施工人员毫无防备，给施工安全带来了极大威胁。部位集中性：大部分岩爆均发生在新开挖的工作面附近，常见的岩爆部位以拱部或拱腰部位为多。这是因为这些部位在开挖后，围岩的应力集中现象较为明显，容易引发岩爆。时间集中性与延续性：岩爆在开挖后陆续出现，多在爆破后24h内发生，延续时间一般为1-2个月，有的延长1年以上，事前一般无明显预兆。在一些深埋隧道施工中，岩爆在开挖后的一段时间内频繁发生，持续影响施工进度和安全。弹射性：岩爆时，岩块自洞壁围岩母体弹射出来，一般呈中厚边薄的不规则片状。弹射出来的岩块速度较快，具有较大的冲击力，容易对施工人员和设备造成伤害。岩爆的危害极大，往往造成开挖工作面的严重破坏、设备损坏和人员伤亡。轻微的岩爆仅有剥落岩片，无弹射现象；严重的可测到4.6级的震级，烈度达7-8度，使地面建筑遭受破坏，并伴有很大的声响。在挪威的某深埋公路隧道施工中，发生了强烈的岩爆，导致隧道局部坍塌，施工设备被严重损坏，多名施工人员受伤，工程进度受到了严重影响。2.1.2塌方塌方是指路面、山体、建筑物或矿井在非人为的自然力情况下，出现塌陷下坠的自然现象，具体表现为道路、堤坝等旁边的陡坡因风化、水浸、震动等影响或坑道、隧道、矿井的顶部因土质岩层松软突然坍塌。在深埋地下洞室施工中，塌方的形成原因较为复杂，主要包括以下几个方面：地质条件：岩体的结构、强度和完整性对塌方的发生起着关键作用。例如，当岩体中存在断层、节理等地质构造时，这些部位的岩体完整性受到破坏，强度降低，在施工扰动和地应力作用下，容易发生塌方。若岩体为破碎的碎裂结构或软弱的页岩等，其承载能力较低，也容易引发塌方。地应力：高地应力状态下，岩体承受的压力过大，当施工开挖改变了岩体的应力分布，使得局部应力超过岩体的强度时，岩体就会发生破坏，从而导致塌方。在一些深埋地下洞室中，由于地应力较高，施工过程中经常出现洞顶坍塌、边墙失稳等塌方现象。地下水：地下水的存在会对岩体产生多种不利影响。它会软化岩体，降低岩体的强度；增加岩体的重量，使岩体所受的重力荷载增大；产生动水压力和孔隙水压力，进一步破坏岩体的稳定性。在富水地层中进行地下洞室施工时，若地下水处理不当，很容易引发塌方事故。施工方法：不合理的施工方法也是导致塌方的重要原因之一。如开挖顺序不当，可能会导致岩体的应力分布不均匀，引发局部失稳；爆破参数选择不合理，爆破震动过大，会对岩体造成过度扰动，破坏岩体的完整性；支护不及时或支护强度不足，无法有效支撑岩体，也会使岩体在自重和地应力作用下发生塌方。塌方的表现形式主要有洞顶坍塌、边墙失稳等。洞顶坍塌时，顶部岩体突然垮落，形成塌腔，严重时可能导致整个洞室被掩埋；边墙失稳则表现为边墙岩体向外鼓出、滑落，影响洞室的正常使用和后续施工。在梅大高速茶阳路段“5・1”塌方灾害中，由于长时间持续性降水与多种因素叠加耦合作用，导致地下水持续累积、水位升高，路堤底部及基底软化，抗剪强度降低，最终引发路堤上部填土塌方，造成了重大人员伤亡和财产损失。2.1.3大变形大变形是指围岩在地下洞室开挖后发生的显著位移和应力重分布，导致围岩出现较大的变形和破坏。这种变形通常表现为围岩的塌落、片帮、挤压、拱起等现象，同时伴随着应力状态的显著改变，对地下洞室的结构安全和稳定性产生严重影响。大变形对洞室结构和施工的影响主要体现在以下几个方面：影响洞室结构安全：大变形可能导致洞室的形状发生改变，使洞室的净空尺寸减小，影响洞室的正常使用功能。严重的大变形甚至会导致洞室失稳，出现坍塌等事故，威胁施工人员的生命安全和工程的顺利进行。在一些软岩隧道施工中，由于岩体的强度较低，在高地应力作用下，容易发生大变形，导致隧道衬砌开裂、变形，甚至出现坍塌现象。增加施工难度：大变形会使施工过程中的支护难度加大，需要采用更加强劲的支护措施来控制围岩变形。同时，大变形还可能导致施工设备无法正常作业，影响施工进度。在处理大变形问题时，往往需要进行多次支护和加固，增加了工程成本和施工周期。影响工程质量：大变形可能会对洞室的衬砌结构造成破坏，降低衬砌的承载能力和防水性能，从而影响工程的质量和使用寿命。大变形还可能导致洞室周围的岩体产生松动，增加了后期运营的安全隐患。2.2灾害发生机理2.2.1应力集中与释放在深埋地下洞室施工过程中，应力集中与释放是硬岩应力型灾害发生的关键因素。当洞室开挖时，原有的岩体平衡状态被打破，应力重新分布。在洞室周边，由于岩体的约束条件改变，应力会出现集中现象。根据弹性力学理论，对于圆形洞室，在均匀地应力场作用下，洞周的切向应力会显著增大，其表达式为：\sigma_{\theta}=\frac{\sigma_{1}+\sigma_{3}}{2}(1+\frac{r_{0}^{2}}{r^{2}})+\frac{\sigma_{1}-\sigma_{3}}{2}(1+3\frac{r_{0}^{4}}{r^{4}})\cos2\theta其中，\sigma_{\theta}为洞周切向应力，\sigma_{1}和\sigma_{3}分别为最大和最小主应力，r_{0}为洞室半径，r为计算点到洞室中心的距离，\theta为计算点与最大主应力方向的夹角。从公式可以看出，在洞壁处（r=r_{0}），切向应力会达到最大值，是原岩应力的数倍。当洞周应力集中程度超过岩石的强度极限时，岩石就会发生破坏。在岩爆发生过程中，岩体中积聚的弹性应变能会突然释放。根据能量守恒定律，岩体变形过程中储存的弹性应变能U可以表示为：U=\frac{1}{2}\int_{V}\sigma_{ij}\epsilon_{ij}dV其中，\sigma_{ij}和\epsilon_{ij}分别为应力张量和应变张量，V为岩体体积。当岩体破坏时，这些积聚的弹性应变能会迅速转化为动能，导致岩石碎片高速弹射，从而引发岩爆灾害。在一些深埋隧道施工中，由于地应力较高，洞室开挖后，洞周应力集中，岩体中的弹性应变能大量积聚。当能量达到一定程度时，就会突然发生岩爆，岩石碎片像炮弹一样从洞壁射出，对施工人员和设备造成严重威胁。2.2.2岩石力学性质岩石的力学性质对硬岩应力型灾害的发生有着重要影响。岩石的强度是抵抗破坏的能力，包括抗压强度、抗拉强度和抗剪强度等。一般来说，抗压强度大于抗拉强度和抗剪强度。在深埋地下洞室中，岩石受到高地应力作用，当洞室开挖引起的应力变化超过岩石的强度时，就容易发生灾害。对于抗压强度较低的岩石，在高围压下，容易发生塑性变形，导致洞室围岩的大变形；而对于抗拉强度低的岩石，在洞室周边受拉区域，容易出现拉伸裂缝，进而引发岩石的破坏。岩石的脆性也是影响灾害发生的重要因素。脆性岩石在受力时，变形较小，当应力达到一定程度时，会突然发生破裂，释放出大量能量。岩爆多发生在脆性岩石中，这是因为脆性岩石能够储存较多的弹性应变能，一旦破坏，能量瞬间释放，造成强烈的冲击。相比之下，塑性岩石在受力时会发生较大的塑性变形，能够消耗部分能量，不易发生突然的脆性破坏，因此发生岩爆的可能性相对较小。岩石的弹性模量和泊松比也与灾害发生密切相关。弹性模量反映了岩石抵抗弹性变形的能力，弹性模量越大，岩石越坚硬，在相同的应力作用下，变形越小。泊松比则表示岩石在横向应变与纵向应变的比值，泊松比越大，岩石在受力时的横向变形越大。在深埋地下洞室中，高弹性模量和低泊松比的岩石更容易积聚弹性应变能，增加了岩爆等灾害发生的风险。2.2.3地质构造因素地质构造如断层、节理等对硬岩应力型灾害的发生起着重要作用。断层是岩石中的破裂面，两侧岩石发生了相对位移。在断层附近，地应力分布往往不均匀，存在应力集中现象。当洞室开挖遇到断层时，由于断层带岩体破碎，强度降低，容易引发塌方、大变形等灾害。断层还可能导致地下水的富集和运移，进一步降低岩体的稳定性。在某深埋地下洞室施工中，遇到一条正断层，断层带内岩体破碎，呈糜棱岩化，开挖过程中，洞室顶部和边墙出现了严重的坍塌，支护难度极大。节理是岩石中的裂隙，它们将岩体切割成不同的块体。节理的发育程度、间距、方向等都会影响岩体的力学性质和稳定性。节理密集的区域，岩体的完整性被破坏，强度降低，在洞室开挖时，容易沿着节理面发生滑动、脱落等破坏现象。节理的方向与洞室轴线的夹角也会影响灾害的发生。当节理方向与洞室轴线平行时，在洞室周边的切向应力作用下，岩体容易沿着节理面发生张裂破坏；而当节理方向与洞室轴线垂直时，岩体在轴向应力作用下，容易发生剪切破坏。2.3灾害影响因素2.3.1地应力地应力是影响深埋地下洞室施工期硬岩应力型灾害风险的关键因素之一。地应力的大小和方向对灾害风险有着显著的影响。在深埋地下洞室中，高地应力状态较为常见，当岩体所受的地应力超过其自身强度时，就容易引发岩爆、板裂破坏等灾害。地应力的大小直接关系到岩体中储存的弹性应变能的多少。一般来说，地应力越大，岩体中储存的弹性应变能就越多，发生灾害的可能性和强度也就越大。在锦屏二级水电站引水隧洞施工中，最大水平主应力高达30MPa以上，高应力使得岩体中的弹性应变能大量积聚，从而导致岩爆灾害频繁发生。通过对该工程岩爆发生区域的地应力监测和分析发现，岩爆发生时，地应力集中区域的应力值明显高于周边区域，且岩爆的强度与地应力的大小呈正相关关系。地应力的方向也会对灾害风险产生重要影响。当洞室轴线与最大主应力方向夹角较小时，洞室周边的应力集中程度相对较低，灾害发生的可能性也较小；而当夹角较大时，应力集中程度会显著增加，灾害风险也随之增大。对于圆形洞室，当洞室轴线与最大主应力方向平行时，洞周切向应力相对较小；当洞室轴线与最大主应力方向垂直时，洞周切向应力会达到最大值，此时岩爆等灾害发生的风险显著增加。在某深埋地下洞室工程中，通过数值模拟分析了不同地应力方向下洞室周边的应力分布情况，结果表明，当地应力方向与洞室轴线夹角为90°时，洞室拱顶和拱腰处的应力集中系数明显高于其他夹角情况，该区域发生岩爆的可能性也更大。2.3.2岩石特性岩石的特性与硬岩应力型灾害密切相关，其中抗压强度、弹性模量等特性对灾害的发生有着重要影响。岩石的抗压强度是指岩石在单向压力作用下抵抗破坏的能力，它是衡量岩石强度的重要指标。抗压强度较低的岩石在高地应力作用下，更容易发生塑性变形和破坏，从而增加了灾害发生的风险。在一些软弱岩石中，如页岩、泥岩等，由于其抗压强度相对较低，在深埋地下洞室施工中，容易出现洞室围岩的大变形、坍塌等灾害。弹性模量反映了岩石抵抗弹性变形的能力，弹性模量越大，岩石越坚硬，在相同的应力作用下，变形越小。高弹性模量的岩石在受力时，能够储存更多的弹性应变能，一旦发生破坏，释放出的能量也更大，容易引发强烈的岩爆等灾害。在硬岩中，如花岗岩、石英岩等，其弹性模量较高，当岩体中的应力超过其强度时，积聚的弹性应变能会突然释放，导致岩石爆裂，引发岩爆。岩石的脆性也是影响灾害发生的重要因素。脆性岩石在受力时，变形较小，当应力达到一定程度时，会突然发生破裂，释放出大量能量。岩爆多发生在脆性岩石中，这是因为脆性岩石能够储存较多的弹性应变能，一旦破坏，能量瞬间释放，造成强烈的冲击。相比之下，塑性岩石在受力时会发生较大的塑性变形，能够消耗部分能量，不易发生突然的脆性破坏，因此发生岩爆的可能性相对较小。2.3.3洞室形状与尺寸洞室的形状和尺寸对深埋地下洞室施工期硬岩应力型灾害风险有着显著影响，它们主要通过改变应力分布来影响灾害风险。不同形状的洞室在相同的地应力条件下，其周边的应力分布存在明显差异。圆形洞室的应力分布相对较为均匀，在均匀地应力场中，圆形洞室周边的切向应力在洞壁处达到最大值，且沿圆周方向分布较为对称；而矩形、马蹄形等非圆形洞室，由于其形状的不规则性，在洞室的拐角、边缘等部位容易出现应力集中现象。在矩形洞室的拐角处，应力集中系数可达到圆形洞室的数倍，这使得这些部位更容易发生岩石破坏，从而增加了灾害发生的风险。洞室尺寸的大小也会影响应力分布和灾害风险。随着洞室尺寸的增大，洞室周边的应力集中程度会相应增加。对于大跨度的洞室，其顶板和边墙所承受的压力更大，更容易出现变形和破坏。在某深埋地下厂房的施工中，由于洞室跨度较大，顶板在高地应力作用下发生了板裂破坏，导致部分岩体脱落，对施工安全造成了严重威胁。通过数值模拟分析发现，随着洞室跨度的增加，洞室顶板的最大拉应力和最大压应力均显著增大，当应力超过岩石的强度时，就会引发灾害。2.3.4施工方法不同的施工方法如钻爆法、TBM法对深埋地下洞室施工期硬岩应力型灾害的发生有着不同程度的影响。钻爆法是通过钻孔、装药、爆破等工序来开挖洞室，在爆破过程中，会产生强烈的震动和冲击，对围岩造成扰动。这种扰动可能会破坏岩体的原有结构，使岩体的强度降低，从而增加了灾害发生的风险。爆破震动可能会引发岩体中的节理、裂隙扩展，导致岩体的完整性受损，进而引发岩爆、塌方等灾害。在某深埋隧道采用钻爆法施工时，由于爆破参数不合理，爆破震动过大，导致洞室周边的岩体出现了大量裂缝，随后发生了小规模的岩爆和塌方事故。TBM法是利用隧道掘进机进行洞室开挖，与钻爆法相比，TBM法对围岩的扰动相对较小，能够较好地保持岩体的完整性。但TBM法在施工过程中，刀具与岩石的摩擦会产生热量，可能会使岩石的力学性质发生变化。在硬岩中，刀具的高速切削会使岩石局部温度升高，导致岩石的脆性增加，在高地应力条件下，更容易引发岩爆。TBM法在遇到地质条件复杂的地段时，如断层、破碎带等，施工难度较大，也可能会增加灾害发生的风险。三、硬岩应力型灾害风险评估方法基础3.1风险评估基本理论3.1.1风险的定义与度量风险的定义在不同领域和研究中存在多种表述，但核心都围绕不确定性和不利影响。国际标准化组织（ISO）将风险定义为“不确定性对目标的影响”，这种影响可能是正面的，也可能是负面的，但在灾害风险评估中，我们主要关注负面的影响，即可能导致灾害发生并造成人员伤亡、财产损失、工程延误等不利后果的情况。从数学角度看，风险表明的是各种不利结果发生的可能性，它与不确定性密切相关，不确定性越大，风险往往也就越高。风险度量是对风险进行量化描述的过程，常用的指标包括概率和损失程度。概率用于衡量风险事件发生的可能性大小，它可以通过历史数据统计、专家经验判断或基于概率模型计算得出。在岩爆风险评估中，可以通过对类似地质条件和工程环境下的岩爆发生次数进行统计分析，得到岩爆发生的概率。损失程度则是指风险事件发生后所造成的各种损失的大小，包括直接经济损失，如设备损坏、工程修复费用等；间接经济损失，如工程延误导致的额外成本、生产中断造成的经济损失等；以及人员伤亡、环境破坏等非经济损失。在衡量损失程度时，需要将不同类型的损失进行量化处理，以便进行综合评估。对于人员伤亡，可以采用生命价值评估方法，将其转化为经济价值；对于环境破坏，可以通过评估环境修复成本等方式进行量化。除了概率和损失程度，方差和标准差也是常用的风险度量指标。方差是各个数据与其算术平均数的离差平方和的平均数，标准差是方差的平方根。在投资风险评估中，常用方差和标准差来衡量投资收益的波动程度，波动越大，风险越高。在硬岩应力型灾害风险评估中，方差和标准差可以用于衡量风险评估结果的不确定性程度。如果不同评估方法或不同评估人员得到的岩爆风险评估结果的方差或标准差较大，说明评估结果的不确定性较高，需要进一步分析和验证。3.1.2风险评估的流程风险评估是一个系统的过程，一般包括风险识别、风险分析、风险评价等步骤，这些步骤相互关联，共同构成了风险评估的整体框架。风险识别是风险评估的首要步骤，其目的是找出可能导致灾害发生的各种风险因素。在深埋地下洞室施工期硬岩应力型灾害风险评估中，风险识别需要全面考虑地质、工程、环境等多方面因素。地质因素包括地应力、岩石特性、地质构造等；工程因素涵盖洞室形状与尺寸、施工方法等；环境因素有地下水、地震活动等。通过收集工程地质勘察资料、现场监测数据、历史案例分析等方式，对这些因素进行逐一排查和梳理，确定潜在的风险源。在收集某深埋地下洞室的地质勘察资料时，发现存在一条大型断层，断层带岩体破碎，这就识别出了一个可能导致塌方、大变形等灾害的风险因素。还可以采用头脑风暴、专家访谈等方法，充分发挥专家的经验和知识，尽可能全面地识别风险因素。风险分析是在风险识别的基础上，对已识别的风险因素进行深入研究，确定其发生的可能性和可能造成的损失程度。对于每个风险因素，需要分析其在不同条件下发生的概率，以及一旦发生可能引发的灾害类型和损失规模。对于地应力这一风险因素，通过现场地应力测量和数值模拟分析，确定在当前施工条件下，由于地应力导致岩爆发生的概率；同时，结合岩爆灾害的历史案例和数值模拟结果，评估岩爆发生后可能对施工人员、设备和工程进度造成的损失程度。风险分析可以采用定性和定量相结合的方法，定性分析主要依靠专家经验和判断，对风险因素进行主观评价；定量分析则运用数学模型和统计方法，对风险进行量化计算，如概率分析、敏感性分析等。风险评价是将风险分析的结果进行综合评估，确定风险的等级和重要性，以便制定相应的风险应对策略。通常使用风险矩阵或风险优先级矩阵来进行风险评价，将风险划分为高、中、低等级。风险矩阵是一种二维矩阵，横坐标表示风险发生的可能性，纵坐标表示风险造成的损失程度，通过将不同风险因素在矩阵中定位，确定其风险等级。对于风险等级较高的因素，应优先采取措施进行防范和控制；对于风险等级较低的因素，可以进行适当的监测和管理。在某深埋地下洞室风险评估中，通过风险矩阵分析，确定岩爆风险为高等级，塌方风险为中等级，大变形风险为低等级，从而针对不同等级的风险制定相应的应对策略，如对岩爆风险采取加强支护、优化施工工艺等措施，对塌方风险加强地质监测和预警，对大变形风险进行定期观测和评估。三、硬岩应力型灾害风险评估方法基础3.2传统风险评估方法3.2.1经验类比法经验类比法是一种基于过往工程经验和相似案例的风险评估方法，其原理是依据已完成的类似工程在施工过程中所遭遇的灾害情况，以及所采取的应对措施和效果，来推断当前工程可能面临的风险。在评估深埋地下洞室施工期硬岩应力型灾害风险时，若有地质条件、洞室尺寸和施工方法等方面相似的工程案例，就可以参考该案例中岩爆、板裂破坏等灾害的发生情况，对当前工程的风险进行初步判断。在某新建深埋铁路隧道的风险评估中，参考了附近已建成的一条地质条件类似的公路隧道的施工经验。该公路隧道在施工过程中，由于地应力较高，且岩石为脆性花岗岩，发生了多次岩爆灾害。基于此经验，对新建铁路隧道进行风险评估时，判断其在相同地质条件和施工方法下，也有较高的岩爆风险，并提前制定了相应的防治措施。经验类比法的优点在于简单易行，不需要复杂的理论计算和专业知识，能够快速得出初步的风险评估结果。它充分利用了已有的工程实践经验，对于一些常见的地质条件和工程类型，具有较高的参考价值。然而，该方法也存在明显的局限性。由于不同工程之间的地质条件、施工工艺等不可能完全相同，即使是相似的工程案例，也可能存在一些关键差异，这些差异可能导致风险评估结果的不准确。经验类比法缺乏对风险的定量分析，无法准确评估风险发生的概率和可能造成的损失程度，难以满足现代工程对风险评估精确性的要求。3.2.2定性评估法定性评估法是一种基于经验和主观判断的风险评估方法，其特点是不依赖于精确的数据计算，而是通过对风险因素的分析和判断，来确定风险的性质和程度。在深埋地下洞室施工期硬岩应力型灾害风险评估中，定性评估法能够快速地对风险进行初步判断，为后续的风险应对提供方向。专家打分法是定性评估法中常用的一种方法。该方法通过邀请相关领域的专家，对影响硬岩应力型灾害风险的各个因素，如地应力、岩石特性、洞室形状与尺寸、施工方法等，按照一定的标准进行打分。通常采用1-5分或1-10分的评分标准，分数越高表示风险越大。然后根据各个因素的得分，综合评估硬岩应力型灾害的风险等级。假设邀请了5位专家对某深埋地下洞室的岩爆风险进行评估，对于地应力因素，专家们分别给出了4分、3分、4分、5分、3分的评分，经过统计计算，地应力因素的平均得分为3.8分。按照同样的方法，对其他风险因素进行评分，并根据预先设定的风险等级标准，确定该洞室的岩爆风险等级为较高。定性评估法的优点是操作简单、成本较低，能够快速地对风险进行评估，适用于对风险进行初步筛选和判断。它能够充分发挥专家的经验和知识，考虑到一些难以量化的因素对风险的影响。但定性评估法也存在主观性强的缺点，不同专家由于经验、知识背景和判断标准的不同，可能会给出不同的评分结果，导致评估结果的可靠性和一致性较差。该方法缺乏精确的量化分析，无法准确地评估风险的大小和发生概率，对于一些对风险评估精度要求较高的工程，定性评估法的结果可能无法满足需求。3.2.3定量评估法定量评估法是一种通过量化指标和数学模型来评估风险的方法，它能够对风险进行精确的度量和分析。在深埋地下洞室施工期硬岩应力型灾害风险评估中，定量评估法可以更加准确地评估风险发生的概率和可能造成的损失程度，为工程决策提供科学依据。数值模拟法是定量评估法中常用的一种方法，它利用有限元软件ANSYS、ABAQUS等，建立深埋地下洞室的数值模型。通过模拟施工过程中围岩的应力、应变分布和变形破坏过程，预测硬岩应力型灾害的发生可能性。在建立数值模型时，需要输入地应力、岩石力学参数、岩体结构等数据，这些数据可以通过现场测试和室内试验获取。利用ANSYS软件对某深埋地下洞室进行数值模拟，模拟结果显示，在洞室开挖过程中，洞室周边的应力集中区域出现了较大的塑性变形，根据设定的破坏准则，判断该区域有较高的岩爆风险。通过数值模拟，还可以分析不同施工方案和支护措施对风险的影响，为优化施工方案提供参考。可靠性分析法也是定量评估法的一种，它通过对风险因素的不确定性进行分析，计算风险发生的概率和可靠性指标。在深埋地下洞室施工期硬岩应力型灾害风险评估中，地应力、岩石力学参数等风险因素存在一定的不确定性，可靠性分析法可以考虑这些不确定性因素，更加准确地评估风险。采用蒙特卡罗模拟法，对某深埋地下洞室的岩爆风险进行可靠性分析。通过随机生成地应力、岩石力学参数等数据，模拟多次开挖过程，统计岩爆发生的次数，从而计算出岩爆发生的概率。根据计算结果，评估该洞室的岩爆风险水平，并为制定风险应对措施提供依据。定量评估法的优点是评估结果准确、可靠，能够为工程决策提供科学依据。它可以考虑到各种风险因素的不确定性，对风险进行全面、系统的分析。但定量评估法也存在一些缺点，它需要大量的基础数据和专业知识，数据的准确性和可靠性对评估结果影响较大。建立和求解数学模型的过程较为复杂，需要耗费较多的时间和计算资源，对于一些小型工程或时间紧迫的项目，可能不太适用。3.3动态评估方法的优势与需求在深埋地下洞室施工中，动态评估方法相较于传统风险评估方法具有显著优势，能够更好地适应施工过程中的复杂情况。传统风险评估方法多基于施工前的地质勘察和设计资料，对施工过程中风险因素的动态变化考虑不足。而动态评估方法能够实时反映风险变化，这是其最突出的优势之一。在施工过程中，地应力、岩石力学性质等风险因素会随着洞室开挖、支护等作业而发生改变，动态评估方法通过实时监测这些因素的变化，及时调整风险评估结果，为施工决策提供更准确的依据。以某深埋隧道施工为例，在施工初期，根据传统评估方法，基于前期地质勘察数据，评估该地段岩爆风险为中等。但随着施工的推进，通过实时监测发现，由于洞室开挖导致地应力重新分布，洞周局部区域的应力集中程度远超预期，且岩石的脆性特征在开挖扰动下更加明显。此时，动态评估方法及时捕捉到这些变化，重新评估该地段岩爆风险为高，施工方根据新的评估结果，立即调整施工方案，加强了支护措施，从而有效避免了岩爆灾害的发生。若采用传统评估方法，无法及时发现风险变化，可能会导致施工安全事故。在深埋地下洞室施工中采用动态评估方法具有重要的必要性。深埋地下洞室施工环境复杂，地质条件多变，存在诸多不确定性因素。在施工过程中，可能会遇到未探明的断层、破碎带等地质构造，这些情况会显著改变风险状况。施工工艺的调整、施工进度的变化等也会对风险产生影响。因此，需要一种能够及时响应这些变化的评估方法，动态评估方法正好满足了这一需求。随着施工技术的不断发展，对施工安全和工程质量的要求越来越高。动态评估方法能够为施工安全提供有力保障，通过实时评估风险，及时发现潜在的安全隐患，提前采取防范措施，降低灾害发生的可能性和危害程度，从而保障施工人员的生命安全和工程的顺利进行。动态评估方法还能够优化施工方案，根据风险评估结果，施工方可以合理调整施工顺序、选择合适的施工工艺和支护方式，提高施工效率，降低工程成本。四、深埋地下洞室施工期硬岩应力型灾害风险动态评估模型构建4.1数据采集与处理4.1.1地应力测量地应力测量是获取深埋地下洞室施工区域地应力信息的关键环节，对于评估硬岩应力型灾害风险至关重要。目前，常用的地应力测量方法包括应力解除法、水压致裂法等，这些方法各有其特点和适用范围。应力解除法是通过从原岩中取出岩芯，使岩芯解除应力，测量岩芯的应变恢复，从而计算出原岩应力。根据测量应变的方式不同，应力解除法又可细分为孔径变形法、孔底应变法和空心包体应变法等。孔径变形法利用孔壁应变计测量钻孔孔径的变形，进而计算地应力；孔底应变法则通过测量钻孔底部的应变来推算地应力；空心包体应变法是将空心包体应变计安装在钻孔中，测量解除应力过程中的应变变化，该方法能够测量三维地应力，精度相对较高。在某深埋地下洞室工程中，采用空心包体应变法进行地应力测量，在钻孔中安装空心包体应变计，通过套钻解除应力，测量应变恢复值，利用弹性力学理论计算出地应力的大小和方向。这种方法适用于完整岩体，能够较为准确地获取地应力信息，但操作相对复杂，对技术要求较高。水压致裂法是向钻孔内注入高压水，使孔壁岩石产生裂缝，通过测量水压和裂缝的扩展情况来确定地应力。该方法不需要预先知道岩石的力学参数，适用于深部岩体地应力测量。在实施过程中，将封隔器安装在钻孔中，封隔一段钻孔，然后向封隔段内注入高压水，当水压达到一定值时，孔壁岩石产生裂缝，记录此时的水压（破裂压力）和裂缝重新张开时的水压（重张压力），根据相关理论公式计算地应力。水压致裂法的优点是能够测量深部地应力，且对岩体的完整性要求相对较低，但测量结果受到钻孔条件、岩石渗透性等因素的影响。在某深埋隧道工程中，采用水压致裂法测量地应力，通过多次试验，准确获取了深部岩体的地应力数据，为工程设计和施工提供了重要依据。4.1.2岩石力学参数测试岩石力学参数测试是评估深埋地下洞室施工期硬岩应力型灾害风险的重要基础，通过测试获取岩石的力学参数，能够深入了解岩石的力学特性和破坏机制。测试内容主要包括岩石的单轴抗压强度、三轴压缩试验、抗拉强度等。单轴抗压强度测试是将岩石试件置于压力试验机上，在无侧向约束的条件下，以规定的加载速度施加轴向压力，直至试件破坏，记录破坏时的最大荷载，根据试件的横截面积计算单轴抗压强度。单轴抗压强度是衡量岩石强度的基本指标之一，它反映了岩石在单向压力作用下抵抗破坏的能力。在某深埋地下洞室工程的岩石力学参数测试中，采用电子万能压力试验机对岩石试件进行单轴抗压强度测试，严格按照试验标准控制加载速度，得到了准确的单轴抗压强度数据。三轴压缩试验是在有侧向压力的条件下对岩石试件进行加载，模拟岩石在地下实际受力状态。试验时，将岩石试件置于三轴压力室内，先施加一定的围压，然后逐渐增加轴向压力，直至试件破坏。通过三轴压缩试验，可以得到岩石在不同围压下的强度和变形特性，为分析深埋地下洞室围岩的力学行为提供重要依据。在某深埋隧道工程的岩石力学参数测试中，进行了三轴压缩试验，研究了不同围压下岩石的抗压强度、变形模量等参数的变化规律，发现随着围压的增加，岩石的抗压强度显著提高，变形模量也有所增大。抗拉强度测试是测定岩石抵抗拉伸破坏的能力，常用的方法有直接拉伸法和巴西劈裂法。直接拉伸法是将岩石试件两端固定，施加轴向拉力，直至试件拉断，测量破坏时的拉力，计算抗拉强度。但直接拉伸法由于试件制备和试验操作难度较大，实际应用较少。巴西劈裂法是将圆盘形岩石试件置于压力机上，在直径方向施加线性分布的压力，使试件在直径平面内产生拉应力，当拉应力达到岩石的抗拉强度时，试件沿直径方向劈裂破坏，根据破坏时的荷载和试件尺寸计算抗拉强度。在某深埋地下洞室工程的岩石力学参数测试中，采用巴西劈裂法测试岩石的抗拉强度，该方法操作相对简单，得到的抗拉强度数据能够满足工程分析的需要。4.1.3监测数据收集监测数据收集是深埋地下洞室施工期硬岩应力型灾害风险动态评估的重要环节，通过收集各种监测数据，能够实时了解洞室施工过程中的围岩状态和风险变化情况。监测数据来源广泛，类型多样，主要包括位移监测、应力监测等。位移监测是通过在洞室围岩表面或内部布置位移传感器，测量围岩的位移变化。常用的位移监测仪器有全站仪、收敛计、多点位移计等。全站仪可以测量洞室周边的三维位移，通过定期测量，能够获取洞室周边不同位置的位移变化情况，分析位移的发展趋势。收敛计主要用于测量洞室周边两点之间的相对位移，如洞室的拱顶下沉和周边收敛，通过测量收敛值，可以判断洞室围岩的稳定性。多点位移计则可以测量围岩内部不同深度的位移，了解围岩内部的变形分布情况。在某深埋地下洞室施工中，在洞室周边布置了全站仪监测点，定期测量洞室周边的位移，同时在拱顶和边墙布置收敛计，监测拱顶下沉和周边收敛情况，还在围岩内部安装了多点位移计，监测围岩内部的位移变化。通过这些位移监测数据，能够及时发现围岩的变形异常，为评估灾害风险提供重要依据。应力监测是通过在洞室围岩中埋设应力传感器，测量围岩的应力变化。常用的应力监测仪器有压力盒、应变计等。压力盒可以测量围岩的接触压力，通过在洞室衬砌与围岩之间埋设压力盒，能够监测围岩对衬砌的压力大小和分布情况。应变计则可以测量围岩的应变，根据应变与应力的关系，计算出围岩的应力。在某深埋隧道施工中，在围岩中埋设了压力盒和应变计，监测围岩的应力变化。当发现应力突然增大或出现异常分布时，及时进行分析，判断是否存在灾害风险，以便采取相应的措施进行处理。4.1.4数据预处理与分析数据预处理是确保监测数据质量和可靠性的关键步骤，能够提高后续数据分析的准确性和有效性。数据清洗是数据预处理的重要环节，主要是去除数据中的噪声、异常值和缺失值。噪声数据是由于测量仪器的误差、环境干扰等因素导致的数据偏差，通过滤波、平滑等方法可以去除噪声。异常值是与其他数据明显不同的数据点，可能是由于测量错误、设备故障等原因产生的。对于异常值，可以采用统计方法进行识别，如3σ准则，即数据点与均值的偏差超过3倍标准差时，可视为异常值，然后根据具体情况进行处理，如删除异常值或进行修正。缺失值的处理方法有多种，如均值填充、中位数填充、插值法等。均值填充是用该变量的均值来填充缺失值；中位数填充则是用中位数进行填充；插值法是根据相邻数据点的关系，通过插值公式计算缺失值。在某深埋地下洞室监测数据处理中，采用3σ准则识别并删除了部分异常值，对于缺失值，根据数据的时间序列特点，采用线性插值法进行了填充，保证了数据的完整性和可靠性。数据分析技术在深埋地下洞室施工期硬岩应力型灾害风险评估中起着重要作用，能够从大量的监测数据中提取有价值的信息，为风险评估提供支持。趋势分析是数据分析的常用方法之一，通过对监测数据随时间的变化趋势进行分析，能够判断洞室围岩的稳定性和灾害风险的发展趋势。在位移监测数据的趋势分析中，通过绘制位移-时间曲线，可以直观地看出位移的变化情况，如果位移呈现持续增长且增长速率加快的趋势，可能预示着洞室围岩存在失稳的风险。相关性分析用于研究不同监测数据之间的相关性，找出影响灾害风险的关键因素。在某深埋地下洞室的监测数据中，通过相关性分析发现，地应力与位移之间存在显著的正相关关系，即地应力越大，位移也越大，这为进一步分析灾害风险提供了重要线索。主成分分析是一种数据降维方法，能够将多个相关变量转化为少数几个不相关的综合变量，即主成分。通过主成分分析，可以减少数据的维度，提取数据的主要特征，便于后续的分析和处理。在某深埋地下洞室的监测数据处理中，采用主成分分析对位移、应力、岩石力学参数等多个变量进行降维，得到了几个主成分，这些主成分包含了原始数据的大部分信息，且彼此之间相关性较低，为建立风险评估模型提供了更有效的输入变量。四、深埋地下洞室施工期硬岩应力型灾害风险动态评估模型构建4.2风险评估指标体系建立4.2.1指标选取原则指标选取遵循科学性、全面性、可操作性等原则。科学性原则要求所选指标能够准确反映硬岩应力型灾害风险的本质特征，基于岩石力学、地质工程等相关科学理论，确保指标的合理性和可靠性。在选取地应力指标时，应充分考虑地应力的大小、方向和分布对灾害风险的影响，依据弹性力学和岩石力学理论，确定能够准确描述地应力状态的指标，如最大主应力、最小主应力、应力比等。全面性原则强调指标体系要涵盖影响硬岩应力型灾害风险的各个方面，包括地质条件、岩石特性、洞室特征、施工方法等。地质条件方面，要考虑地层岩性、地质构造、地下水等因素；岩石特性包括抗压强度、弹性模量、脆性等；洞室特征涉及洞室形状、尺寸、埋深等；施工方法涵盖钻爆法、TBM法等不同施工工艺对风险的影响。通过全面考虑这些因素，确保风险评估的完整性和准确性。可操作性原则要求所选指标的数据易于获取和测量，在实际工程中具有可行性。对于地应力测量，可以采用应力解除法、水压致裂法等成熟的测量方法；岩石力学参数测试可通过室内试验，如单轴抗压试验、三轴压缩试验等获取。指标的计算和分析方法应简单明了，便于工程技术人员理解和应用，能够为工程决策提供及时有效的支持。4.2.2确定评估指标用于风险评估的具体指标包括地应力指标、岩石强度指标等。地应力指标中，最大主应力\sigma_{1}和最小主应力\sigma_{3}是关键指标，它们直接反映了岩体所受应力的大小和方向，对硬岩应力型灾害的发生起着决定性作用。当最大主应力超过岩石的抗压强度时，容易引发岩爆等灾害。应力比\lambda=\frac{\sigma_{1}}{\sigma_{3}}也具有重要意义，它反映了地应力的不均匀程度，应力比越大，岩体的受力状态越复杂，灾害发生的风险越高。在某深埋地下洞室工程中，通过地应力测量得到最大主应力为25MPa，最小主应力为5MPa，应力比为5，该地段岩爆风险较高，施工过程中需要特别关注。岩石强度指标中，单轴抗压强度\sigma_{c}是衡量岩石抵抗单向压力能力的重要指标，它反映了岩石的基本强度特性。单轴抗压强度较低的岩石，在高地应力作用下更容易发生破坏，增加了灾害发生的风险。在某深埋隧道工程中，岩石的单轴抗压强度为40MPa，相对较低，在施工过程中，该地段出现了较多的塌方和大变形现象。弹性模量E反映了岩石抵抗弹性变形的能力，弹性模量越大，岩石越坚硬，在相同应力作用下，变形越小。但高弹性模量的岩石在受力时，也更容易积聚弹性应变能，一旦发生破坏，释放出的能量更大，容易引发强烈的岩爆等灾害。在某硬岩深埋地下洞室工程中，岩石的弹性模量高达50GPa，施工过程中，该地段发生了多次强烈的岩爆，对施工安全造成了严重威胁。4.2.3指标权重确定方法确定指标权重的方法有层次分析法、熵权法等。层次分析法（AHP）是一种将与决策总是有关的元素分解成目标、准则、方案等层次，在此基础上进行定性和定量分析的决策方法。在硬岩应力型灾害风险评估中，运用层次分析法确定指标权重时，首先需要构建层次结构模型，将风险评估目标作为目标层，地应力、岩石强度、岩体结构等因素作为准则层，具体的评估指标作为指标层。邀请相关领域的专家，采用1-9标度法对各层次指标进行两两比较，构造判断矩阵。通过计算判断矩阵的特征向量和最大特征值，得到各指标的相对权重。对判断矩阵进行一致性检验，确保权重的合理性。假设在某深埋地下洞室硬岩应力型灾害风险评估中，通过层次分析法确定地应力指标的权重为0.4，岩石强度指标的权重为0.3，岩体结构指标的权重为0.2，其他指标的权重为0.1，这表明在该工程中，地应力对灾害风险的影响最大。熵权法是一种根据指标数据的离散程度来确定权重的方法，数据的离散程度越大，熵值越小，该指标提供的信息量越大，权重也就越大。在硬岩应力型灾害风险评估中，运用熵权法确定指标权重时，首先需要对各指标数据进行标准化处理，消除量纲和数量级的影响。然后计算各指标的熵值和熵权。假设在某深埋地下洞室风险评估中，通过熵权法计算得到地应力指标的熵权为0.35，岩石强度指标的熵权为0.25，岩体结构指标的熵权为0.2，其他指标的熵权为0.2，这说明地应力指标在该工程风险评估中提供的信息量较大，对风险评估结果的影响较为重要。4.3动态评估模型原理与算法4.3.1基于机器学习的模型基于机器学习的风险评估模型在深埋地下洞室施工期硬岩应力型灾害风险评估中具有重要作用，其中神经网络和支持向量机是常用的模型。神经网络是一种模拟人类大脑神经元结构和功能的计算模型，它由输入层、隐藏层和输出层组成，各层之间通过权重连接。在硬岩应力型灾害风险评估中，输入层可以接收地应力、岩石力学参数、洞室形状与尺寸、施工方法等风险评估指标数据，隐藏层对这些数据进行非线性变换和特征提取，输出层则输出风险评估结果，如岩爆、板裂破坏等灾害发生的概率或风险等级。以某深埋地下洞室工程为例，采用BP神经网络进行岩爆风险评估，通过收集大量的工程数据，包括地应力、岩石抗压强度、弹性模量等指标，对神经网络进行训练。训练过程中，不断调整权重和阈值，使网络的输出结果与实际的岩爆情况尽可能接近。经过多次训练和优化，该神经网络能够准确地评估岩爆风险，为工程施工提供了有效的决策依据。支持向量机是一种基于统计学习理论的分类算法，它通过寻找一个最优分类超平面，将不同类别的样本分开。在硬岩应力型灾害风险评估中，支持向量机可以将风险分为不同的等级，如高、中、低风险。对于岩爆风险评估，将发生岩爆的样本和未发生岩爆的样本作为训练数据，支持向量机通过学习这些样本的特征，构建分类模型。当有新的样本数据输入时，模型可以判断该样本属于哪个风险等级。在某深埋隧道工程中，利用支持向量机对岩爆风险进行评估，通过对历史数据的学习和训练，建立了岩爆风险评估模型。该模型在实际应用中，能够准确地预测岩爆风险，为工程施工提供了及时的预警。4.3.2基于数值模拟的模型基于数值模拟的模型在深埋地下洞室施工期硬岩应力型灾害风险评估中发挥着重要作用，有限元法和离散元法是其中常用的方法。有限元法是一种将连续体离散化为有限个单元的数值计算方法，它通过对每个单元进行力学分析，然后将这些单元组合起来，得到整个连续体的力学响应。在硬岩应力型灾害风险评估中，利用有限元软件ANSYS、ABAQUS等，建立深埋地下洞室的数值模型。在建立模型时，需要考虑地应力、岩石力学参数、岩体结构等因素，将这些因素作为输入条件，模拟施工过程中围岩的应力、应变分布和变形破坏过程。通过对模拟结果的分析，可以预测硬岩应力型灾害的发生可能性。在某深埋地下洞室工程中，采用ANSYS软件建立数值模型，模拟洞室开挖过程。通过模拟，得到了洞室周边的应力分布云图和位移矢量图，发现洞室拱顶和边墙处存在应力集中现象，且位移较大，根据这些结果，判断该区域有较高的岩爆和大变形风险，为工程施工提供了重要的参考依据。离散元法是一种适用于非连续介质的数值模拟方法，它将岩体看作是由离散的块体组成，块体之间通过接触力相互作用。在硬岩应力型灾害风险评估中，离散元法可以较好地模拟岩体的节理、裂隙等非连续结构，以及岩体在受力过程中的破裂、滑移等现象。利用离散元软件UDEC、PFC等，建立深埋地下洞室的离散元模型，考虑岩体的节理、裂隙分布，模拟施工过程中岩体的力学行为。在某深埋隧道工程中，采用UDEC软件建立离散元模型，模拟洞室开挖过程中岩体的破坏过程。通过模拟，观察到岩体沿着节理面发生滑移和破裂，最终导致洞室塌方，根据模拟结果，评估该洞室的塌方风险，并提出相应的防治措施。4.3.3模型的动态更新机制模型的动态更新机制是确保深埋地下洞室施工期硬岩应力型灾害风险动态评估准确性的关键。在施工过程中，新的数据不断产生，这些数据反映了施工条件的变化和风险因素的动态演变，因此需要根据新数据对模型进行动态更新。当有新的地应力测量数据、岩石力学参数测试数据或监测数据时，首先对这些新数据进行预处理，包括数据清洗、归一化等操作，以确保数据的质量和一致性。然后，将预处理后的数据输入到已建立的风险评估模型中，根据模型的类型和算法，对模型的参数进行调整和优化。对于基于机器学习的模型，如神经网络，采用反向传播算法等优化方法，根据新数据调整网络的权重和阈值，使模型能够更好地适应新的数据分布和风险特征；对于基于数值模拟的模型，根据新数据更新模型中的参数，如地应力、岩石力学参数等，重新进行数值模拟计算，得到新的风险评估结果。通过动态更新机制，模型能够不断学习和适应新的情况，提高评估的准确性。在某深埋地下洞室施工过程中，随着洞室的开挖，地应力发生了变化，通过实时监测得到了新的地应力数据。将这些新数据输入到基于机器学习的风险评估模型中，模型自动调整参数，重新评估岩爆风险，发现风险等级有所提高。施工方根据更新后的风险评估结果，及时调整施工方案，加强了支护措施，有效地降低了岩爆发生的可能性。动态更新机制还可以通过定期对模型进行重新训练和验证，确保模型的性能和可靠性。定期收集一段时间内的新数据，对模型进行全面的重新训练，使模型能够充分利用新的数据信息，提高评估的准确性和稳定性。五、案例分析5.1工程背景介绍5.1.1项目概况本案例选取锦屏二级水电站引水隧洞工程。锦屏二级水电站位于四川省凉山彝族自治州木里、盐源、冕宁三县交界处的雅砻江干流锦屏大河湾上，是雅砻江下游河段的重要梯级电站之一。该工程主要任务是发电，通过长约16.67km的引水隧洞将雅砻江锦屏大河湾的水引至锦屏二级水电站厂房，利用落差发电。引水隧洞采用4条平行布置，每条隧洞的开挖直径约为12.4m，衬砌后内径约为11.8m。隧洞埋深一般为1500-2000m，最大埋深达到2525m，属于典型的深埋地下洞室工程。工程于2005年11月开工，2013年8月首台机组投产发电，2014年11月全部机组投产发电。5.1.2地质条件项目所在地地质构造复杂，处于川滇南北构造带与康滇歹字形构造的复合部位，区内断裂构造发育，主要有F5、F8、F13等断层。这些断层的存在使得岩体完整性受到破坏，增加了施工难度和风险。岩石类型主要为大理岩，大理岩质地坚硬，抗压强度较高，一般在100-150MPa之间。但在高地应力作用下，大理岩具有较强的脆性，容易发生岩爆等应力型灾害。地应力状态以水平应力为主，最大水平主应力可达30MPa以上，且方向与隧洞轴线夹角较大。高应力状态使得岩体中积聚了大量的弹性应变能，一旦条件满足，就会引发岩爆等灾害。地下水主要为基岩裂隙水，受地形、地质构造和岩石透水性等因素影响，地下水分布不均。在断层破碎带和节理密集区，地下水较为丰富，对岩体的稳定性产生不利影响。5.1.3施工方案项目采用钻爆法和TBM法相结合的施工方法。在地质条件较好的地段，采用TBM法进行快速施工，提高施工效率；在地质条件复杂，如断层破碎带、节理密集区等，采用钻爆法，便于灵活调整施工参数，确保施工安全。施工顺序遵循先主洞后支洞，自上而下、分层分段的原则。先进行主洞的开挖和支护，然后进行支洞的施工。在主洞开挖过程中，根据围岩的稳定性和地质条件，合理划分开挖分层和分段，及时进行支护，确保洞室的稳定。在施工过程中，还采用了超前地质预报、监控量测等技术手段，及时掌握围岩的变化情况，为施工决策提供依据。5.2风险评估实施过程5.2.1数据采集与整理在本项目中，数据采集工作全面且细致。地应力测量采用水压致裂法，在引水隧洞沿线多个关键位置布置测量钻孔，共进行了50次有效测量，获取了不同深度、不同位置的地应力数据，包括最大主应力、最小主应力及其方向等信息。岩石力学参数测试通过室内试验完成，对采集的大理岩岩样进行了100组单轴抗压强度测试、80组三轴压缩试验和60组巴西劈裂法抗拉强度测试，得到了岩石的抗压强度、弹性模量、泊松比、抗拉强度等参数。监测数据收集涵盖位移监测和应力监测。位移监测方面，在洞室周边布置了30个全站仪监测点、20个收敛计监测断面和15个多点位移计监测孔，定期测量洞室周边的三维位移、拱顶下沉和周边收敛以及围岩内部不同深度的位移；应力监测通过在围岩中埋设20个压力盒和30个应变计，监测围岩的接触压力和应变，进而计算出应力。数据整理过程中，首先对采集到的原始数据进行分类存储，建立了详细的数据档案。然后进行数据清洗，去除因测量仪器故障、环境干扰等原因产生的噪声数据和异常值，对于少量缺失值，采用插值法进行补充。通过数据整理，确保了数据的准确性和完整性，为后续的风险评估工作提供了可靠的数据支持。5.2.2指标计算与分析根据采集的数据，计算风险评估指标的值。对于地应力指标，计算最大主应力\sigma_{1}、最小主应力\sigma_{3}和应力比\lambda=\frac{\sigma_{1}}{\sigma_{3}}。经计算，引水隧洞部分地段最大主应力达到32MPa，最小主应力为6MPa，应力比约为5.33，表明该地段地应力不均匀程度较高。岩石强度指标方面，单轴抗压强度\sigma_{c}平均值为120MPa，弹性模量E平均值为60GPa。分析这些指标与硬岩应力型灾害风险的关系，地应力指标中，最大主应力和应力比越大，岩爆等灾害发生的风险越高；岩石强度指标中，单轴抗压强度较低时，塌方、大变形等灾害风险增加，而弹性模量较高时，岩爆风险增大。通过对指标的分析，初步判断部分高地应力、岩石脆性较大的地段岩爆风险较高，部分岩石强度相对较低的地段塌方和大变形风险不容忽视。5.2.3风险评估结果采用构建的基于机器学习的神经网络风险评估模型，结合计算得到的风险评估指标值，对锦屏二级水电站引水隧洞施工期硬岩应力型灾害风险进行评估。评估结果显示，在部分深埋地段和地质构造复杂区域，岩爆风险等级为高，发生概率约为0.7；在一些岩石强度相对较低、节理裂隙发育的地段，塌方风险等级为中，发生概率约为0.4；大变形风险在整体施工区域相对较低，大部分地段风险等级为低，发生概率约为0.2。将风险评估结果与实际施工情况进行对比验证，在评估为岩爆高风险的地段，施工过程中确实发生了多次岩爆现象，且发生时间和位置与评估结果基本相符；塌方风险评估为中的地段，也出现了小规模的塌方情况。这表明所构建的风险评估模型具有较高的准确性和可靠性，能够为工程施工提供有效的风险预警和决策支持。五、案例分析5.3评估结果验证与分析5.3.1与实际灾害情况对比将风险评估结果与锦屏二级水电站引水隧洞施工期实际发生的硬岩应力型灾害情况进行对比。在评估为岩爆高风险的地段，施工过程中发生了多次不同程度的岩爆事件。在某深埋地段，评估岩爆风险等级为高，发生概率约为0.7。实际施工中，该地段在开挖后不久就发生了强烈岩爆，岩石碎片高速弹射，对施工设备造成了严重损坏，施工被迫中断进行处理。在其他评估为岩爆高风险的区域，也陆续出现了岩爆现象，且发生频率和强度与评估结果具有一定的相关性，高风险区域的岩爆发生次数明显多于低风险区域。对于塌方风险评估，在评估为塌方风险等级为中的地段，施工过程中出现了小规模的塌方情况。在某岩石强度相对较低、节理裂隙发育的地段，评估塌方风险为中，发生概率约为0.4。实际施工中，该地段在进行支护作业时，由于岩体的局部失稳，出现了部分洞顶坍塌和边墙滑落的塌方现象，虽然规模较小，但对施工进度和安全也产生了一定影响。而在评估为塌方风险等级为低的地段，施工过程中基本未出现明显的塌方情况，仅有少量的局部掉块现象。5.3.2评估结果的准确性分析分析评估结果与实际情况的符合程度，发现总体上评估结果能够较好地反映实际灾害发生的情况，但仍存在一定的偏差。对于岩爆风险评估，部分地段的评估结果与实际岩爆发生情况较为吻合，准确预测了岩爆的高风险区域和发生概率。由于地质条件的复杂性和不确定性，以及监测数据的局限性，仍有个别地段的岩爆发生情况与评估结果存在一定差异。在某些地质构造复杂的区域，虽然评估为高风险，但实际岩爆发生的时间和强度与预测略有不同，这可能是由于在评估过程中，对一些微观地质构造和岩石力学性质的局部变化考虑不够全面。塌方风险评估方面，评估结果与实际塌方情况的符合程度也较高，准确识别了大部分可能发生塌方的地段。然而，在一些特殊情况下，如遇到未探明的小型断层或地下水异常富集区域，评估结果未能准确预测塌方的发生。这可能是因为在数据采集和分析过程中，对这些特殊地质情况的信息获取不足，导致评估模型对这些因素的考虑不够充分。影响准确性的因素主要包括地质条件的复杂性和不确定性，深埋地下洞室所处的地质环境复杂多变，地应力分布、岩石力学性质、地质构造等因素在空间上存在较大的不均匀性，难以通过有限的监测数据全面准确地掌握。监测数据的准确性和完整性也对评估结果有重要影响，如果监测设备的精度不够、监测点的布置不合理或监测数据存在缺失，都会导致评估模型输入数据的误差，从而影响评估结果的准确性。评估模型本身的局限性也是一个重要因素，虽然本文构建的风险评估模型考虑了多个因素的影响，但仍然无法完全涵盖所有可能影响硬岩应力型灾害发生的因素，模型的算法和参数设置也可能存在一定的优化空间。5.3.3模型的改进建议根据验证结果，提出对风险评估模型的改进建议。在数据采集方面，应进一步优化监测方案，增加监测点的数量和分布范围，提高监测数据的准确性和完整性。特别是对于地质条件复杂的区域，要加密监测点，加强对地质构造、地下水等关键因素的监测，确保能够及时准确地获取相关信息。采用多种监测手段相结合的方式，如地应力监测、位移监测、应力监测、地下水监测等，综合分析监测数据，提高对地质条件变化的敏感度。在模型算法和参数优化方面，应进一步研究和改进模型的算法，提高模型对复杂地质条件和风险因素的适应性。引入更先进的机器学习算法或对现有算法进行改进，如采用深度学习算法对风险因素进行更深入的特征提取和分析，提高模型的预测能力。根据实际工程数据和验证结果，对模型的参数进行进一步优化，通过大量的试验和分析，确定更合理的参数取值，使模型能够更准确地反映硬岩应力型灾害风险与各因素之间的关系。还应加强对特殊地质条件和罕见灾害情况的研究，将这些特殊情况纳入评估模型的考虑范围。针对未探明的断层、破碎带、地下水异常等特殊地质条件，建立相应的风险评估指标和评估方法，提高模型对这些特殊情况的预测能力。对罕见的硬岩应力型灾害情况，如大规模的岩爆集群发生、突发的严重塌方等，通过收集和分析相关案例，总结其发生规律和特征，建立相应的预警机制和应对策略，使模型能够在面对这些罕见灾害时，也能提供有效的风险评估和预警信息。六、风险应对措施与建议6.1风险应对策略6.1.1风险规避风险规避是一种通过改变工程设计方案、施工方法等措施，直接消除或避免硬岩应力型灾害风险的策略。在工程设计阶段，应充分考虑地质条件和地应力状态，合理选择洞室的位置和走向。若工程区域存在高地应力区或地质构造复杂带，应尽量避开这些区域，选择地质条件相对稳定、地应力较低的地段进行洞室布置。在某深埋地下洞室工程规划时，通过详细的地质勘察，发现一处高地应力集中区域，且存在多条断层破碎带。经过评估，若在此区域修建洞室，岩爆和塌方等灾害风险极高。于是，设计团队调整了洞室的走向，绕过了该高风险区域，从而有效降低了灾害发生的可能性。在施工方法选择上，对于硬岩应力型灾害风险较高的地段，应避免采用对围岩扰动较大的施工方法。钻爆法在爆破过程中会产生强烈的震动和冲击，容易引发岩爆等灾害。在高地应力硬岩地段，可优先选择TBM法进行施工，TBM法对围岩的扰动相对较小，能够较好地保持岩体的完整性，降低灾害风险。若必须采用钻爆法施工，则应优化爆破参数，采用微差爆破、光面爆破等技术，控制爆破震动和冲击，减少对围岩的损伤。6.1.2风险减轻风险减轻是通过采取一系列措施，降低硬岩应力型灾害发生的可能性和危害程度的策略。支护措施是减轻灾害风险的重要手段之一。在深埋地下洞室施工中，及时有效的支护能够增强围岩的稳定性，防止岩爆、塌方等灾害的发生。对于岩爆风险较高的地段，可采用高强度的锚杆、锚索支护，将锚杆、锚索深入到岩体内部，与岩体形成一个整体，提高岩体的抗剪强度和抗拉强度，限制岩体的变形和破裂。还可采用喷射混凝土支护，喷射混凝土能够及时封闭围岩表面，防止岩石风化和剥落，同时增加围岩的抗变形能力。在某深埋隧道施工中，针对岩爆风险较高的地段，采用了锚杆、锚索和喷射混凝土联合支护的方式，有效地控制了岩爆的发生，保障了施工安全。应力释放也是减轻灾害风险的有效方法。对于高地应力硬岩，可采用超前钻孔、爆破等方式进行应力释放。超前钻孔是在洞室开挖前，在掌子面前方一定距离处钻设若干钻孔，使岩体中的应力在钻孔周围得到释放，从而降低洞室开挖时的应力集中程度。爆破应力释放则是