Hoek - Brown强度准则在岩体工程中的多维度应用与深度剖析

Hoek-Brown强度准则在岩体工程中的多维度应用与深度剖析一、引言1.1研究背景与意义岩体工程作为现代工程建设的重要组成部分，广泛涉及到水利水电、交通、矿山开采、地下空间开发等诸多领域。在这些工程中，岩体的强度特性直接关系到工程的稳定性、安全性以及经济效益。准确评估岩体强度，对于合理设计工程结构、选择施工工艺以及确保工程长期稳定运行起着决定性作用。例如，在大型水利水电工程中，坝基岩体的强度决定了大坝能否承受巨大的水压力；在隧道工程里，围岩强度影响着隧道的开挖方式和支护措施；在矿山开采中，岩体强度与开采效率、边坡稳定性密切相关。Hoek-Brown强度准则自1980年由E.Hoek和E.T.Brown提出以来，凭借其独特优势在岩体工程领域得到广泛应用。该准则能够有效反映岩石和岩体固有的非线性破坏特点，充分考虑结构面、应力状态对强度的影响，并且可以解释低应力区、拉应力区和最小主应力对强度的作用，还适用于各向异性岩体描述。这使得它在解决复杂岩体工程问题时具有重要价值，为岩体强度分析提供了一种行之有效的工具，极大地推动了岩体工程理论与实践的发展。然而，随着工程建设规模的不断扩大和地质条件的日益复杂，传统的Hoek-Brown强度准则在应用中逐渐暴露出一些局限性，如不能考虑中间主应力的影响、难以准确确定准则中的参数、对各向异性明显的节理岩石适用性差等。这些问题在一定程度上限制了该准则的应用效果，影响了工程设计的准确性和可靠性。因此，深入研究Hoek-Brown强度准则，探讨其在不同岩体工程条件下的应用方法，分析其存在的问题并寻求改进措施，具有重要的理论意义和实际工程价值。通过对该准则的研究，可以进一步完善岩体强度理论体系，为岩体工程设计和施工提供更科学、更准确的理论依据；同时，也有助于提高工程建设的安全性和经济性，减少工程事故的发生，推动岩体工程领域的技术进步。1.2国内外研究现状自1980年Hoek-Brown强度准则提出以来，国内外学者围绕该准则展开了广泛而深入的研究。在准则的发展历程方面，E.Hoek和E.T.Brown于1980年首次提出该准则，用于描述岩石破坏时极限主应力间的非线性经验关系。此后，在1992年对其进行改进，引入反映岩体特征的经验参数，使其可同时应用于岩石和岩体，称为广义Hoek-Brown岩体强度准则，实现了研究对象从岩石到工程岩体的转变。1997年和2002年又对岩体参数的取值方法等进行完善，引入地质强度指标（GSI）和扰动参数D等，不断提高准则对不同岩体条件的适应性。在三维Hoek-Brown强度准则研究方面，由于传统准则未考虑中间主应力对强度的影响，而大量研究和工程实践证实这种影响确实存在。为此，X.Pan等学者先后提出多种三维Hoek-Brown强度准则。如X.Pan和J.A.Hudson基于原准则提出一种三维表达式；B.Singh等基于广义准则提出可考虑中间主应力影响的经验公式；昝月稳等结合该准则和统一强度理论，提出适用于岩石的非线性统一强度准则；张永兴、S.D.Priest等分别通过外接广义准则三轴拉伸破坏点的方法，确定Drucker-Prager强度准则参数，进而提出三维准则。对于Hoek-Brown强度准则的参数研究，准确选取合适的岩石和岩体参数是准则准确应用的关键。国内外学者在这方面进行诸多探索，提出多种参数确定方法，包括基于现场试验、经验公式以及数值模拟等手段。在工程应用方面，该准则在地下采矿、深部空间开挖与支护、露天采矿及边坡稳定性、隧道工程等领域广泛应用。如在某一级水电站工程中，利用该准则计算边坡和洞室岩体强度参数，为工程稳定性分析提供依据；在露天矿山边坡稳定性分析中，结合广义准则估算岩体力学参数，并通过极限平衡法分析整体稳定性。尽管取得上述成果，但目前研究仍存在一些问题。部分三维准则在特定应力条件下无法简化为二维初始准则，不能直接使用原准则参数，其适用性和准确性有待进一步验证。在参数确定方面，不同方法得到的参数可能存在差异，如何更精准、便捷地确定参数仍是研究难点。此外，对于复杂地质条件下岩体，如含有多种结构面、强烈各向异性岩体等，准则的应用效果和改进方向还需深入研究。1.3研究内容与方法本文旨在深入研究Hoek-Brown强度准则在岩体工程中的应用，主要研究内容如下：详细阐述Hoek-Brown强度准则原理：系统介绍该准则的提出背景、基本表达式以及各参数的物理意义，分析其反映岩体非线性破坏特性的机制，深入探讨准则从最初形式到广义形式的发展过程，以及不同阶段参数取值方法的演变，为后续研究奠定理论基础。全面分析准则在岩体工程中的应用案例：广泛收集地下采矿、隧道工程、边坡工程等领域应用Hoek-Brown强度准则的实际案例，详细剖析在不同工程条件下，如不同地质构造、岩体类型和应力状态，该准则如何用于岩体强度计算、稳定性分析以及工程设计，总结应用过程中的经验和存在的问题。深入探讨准则存在的问题及改进措施：针对传统Hoek-Brown强度准则不能考虑中间主应力影响、参数确定困难等问题，深入分析其产生原因，梳理国内外学者提出的改进方法，如改进模型结构、加入裂隙影响因素、优化参数选择方法等，并对这些改进措施的有效性和局限性进行评价，提出进一步改进的思路和方向。开展数值模拟研究：运用数值模拟软件，建立基于Hoek-Brown强度准则的岩体工程模型，模拟不同工况下岩体的力学响应，对比模拟结果与实际工程数据或理论分析结果，验证准则在数值模拟中的适用性和准确性，通过数值模拟分析不同参数对岩体力学行为的影响规律，为工程参数优化提供参考。本文采用以下研究方法：文献研究法：广泛查阅国内外关于Hoek-Brown强度准则的相关文献，包括学术论文、研究报告、工程案例等，全面了解该准则的研究现状、发展历程、应用成果以及存在问题，为本文研究提供理论支持和研究思路。案例分析法：选取具有代表性的岩体工程案例，详细分析Hoek-Brown强度准则在实际工程中的应用过程和效果，从实践角度深入研究该准则的优点和不足，总结经验教训，为准则的改进和更广泛应用提供实际依据。数值模拟法：利用专业的岩土工程数值模拟软件，如FLAC、ANSYS等，建立符合实际工程地质条件的岩体模型，运用Hoek-Brown强度准则进行数值模拟分析。通过改变模型参数，模拟不同工况下岩体的变形、破坏过程，分析岩体的力学特性，与实际工程情况或理论结果进行对比验证，深入研究准则在数值模拟中的应用效果和参数敏感性。二、Hoek-Brown强度准则概述2.1准则的发展历程1980年，EvertHoek和E.T.Brown在研究地下开挖工程时，为解决当时缺乏合适的岩体强度估算理论的问题，依据1936年的混凝土强度理论，推导出Hoek-Brown强度准则的最初公式。该公式主要用于节理岩体，并需根据莫尔圆包络线确定参数。其表达式为：\sigma_{1}=\sigma_{3}+\sqrt{m_{c}\sigma_{c}\sigma_{3}+s\sigma_{c}^{2}}其中，\sigma_{1}、\sigma_{3}为岩石破裂时的最大、最小主应力（压应力为正）；\sigma_{c}为岩样的单轴抗压强度；m_{c}反映岩石的软硬程度，取值范围在0.0000001-25之间，理想完整坚硬岩体取25；s表示岩体破碎程度，取值范围在0-1之间，0代表极破碎岩体，1代表理想完整岩石。该准则初步建立了岩块和岩体破坏时极限主应力之间的非线性关系，为岩体强度分析提供了新的思路，但在实际应用中存在一定局限性，如对岩体质量描述不够全面等。1983年，Hoek对常数m和s进行改进，引入地下水分级（固定取值10）和节理方位分类的概念，使公式能更好地应用于各向异性节理岩体。在1988年，Hoek和Brown将“未扰动”和“扰动”岩体概念引入准则，将近地表岩体的质量等级降低一个级别，使准则不仅适用于地下工程，也适用于边坡工程。并在1989年的RMR公式中，把地下水分级的值修改为15。1992年，Hoek、Wood和Shah等引进材料参数a，考虑节理岩体抗拉强度为0的情况，并应用Balmer的理论计算应力对数目，对Hoek-Brown强度准则进行改进，使其可同时应用于岩石和岩体，称为广义Hoek-Brown岩体强度准则，实现了研究对象从岩石到工程岩体的转变。其表达式为：\sigma_{1}=\sigma_{3}+\left(m_{b}\sigma_{c}\sigma_{3}+s\sigma_{c}^{2}\right)^{a}其中，m_{b}、s、a为反映岩体特征的经验参数，m_{b}针对不同岩体，反映岩体破碎程度，取值范围0.0-1.0，对于完整的岩体（即岩石），s=1.0。1994年和1995年，Hoek和Brown改进并提出Hoek-Brown强度准则的一般形式。考虑岩体是否扰动的情况下，引进地质强度指标（GSI）参数，把岩体分为一般-好（GSI较高）和很差-差（GSI较低）。当GSI\geq25时，m_{b}=m_{i}exp\left(\frac{GSI-100}{28-14D}\right)；当GSI\lt25时，m_{b}=m_{i}exp\left(\frac{GSI-100}{9-3D}\right)。GSI的引入，使得对岩体质量的评价更加全面和准确，进一步提高了准则的适用性。2002年，Hoek、Carranza-Torres及Corkum再次改进Hoek-Brown强度准则，引入GSI、m_{b}、s和a之间的关系，还引入扰动系数D（取值范围为0.0-1.0，现场无扰动岩体为0.0，而非常扰动岩体为1.0），用于表示应力释放和爆炸破坏等对岩体的扰动作用。并且在研究隧道和边坡时，引入摩尔库伦准则中的参数内聚力C和内摩擦角\varphi。当针对隧道等地下硐室时，\sigma_{3max}=\gammaH\left(0.44-0.044\sqrt{\frac{\sigma_{c}}{\gammaH}}\right)；当针对边坡等露天采矿时，\sigma_{3max}=\gammaH\left(0.65-0.065\sqrt{\frac{\sigma_{c}}{\gammaH}}\right)，式中H为硐室埋藏深度或边坡高度，\gamma为容重。2006年，Hoek和Diederichs等针对中国地区的岩土工程现状，引入新的岩体变形模量E_{m}和GSI值，同时引入Sigmoid函数。在无实验室岩石力学数据时，给出评估岩体变形模量的公式E_{m}=MR\times10^{\frac{GSI-100}{40}}，并给出换算系数MR。这些改进使得准则能够更好地适应不同地区的地质条件和工程需求，进一步完善了Hoek-Brown强度准则体系。2.2基本原理与公式推导Hoek-Brown强度准则基于岩石的非线性破坏特性，认为岩石的破坏不仅与应力大小有关，还与岩石内部的结构特征、节理裂隙分布等因素密切相关。其基本原理源于对大量岩石三轴试验资料和岩体现场试验成果的统计分析，结合岩石破坏机理的理论研究，通过试错法导出岩块和岩体破坏时极限主应力之间的关系式。从岩石破坏机理来看，岩石在受力过程中，内部会产生微裂纹，随着应力的增加，微裂纹逐渐扩展、贯通，最终导致岩石破坏。Hoek-Brown强度准则通过引入反映岩石软硬程度的参数m和反映岩体破碎程度的参数s，来考虑这些因素对岩石强度的影响。对于最初的Hoek-Brown强度准则公式\sigma_{1}=\sigma_{3}+\sqrt{m_{c}\sigma_{c}\sigma_{3}+s\sigma_{c}^{2}}，其推导过程基于对岩石三轴试验数据的拟合和分析。假设岩石破坏时的强度可以表示为最小主应力\sigma_{3}和单轴抗压强度\sigma_{c}的函数，通过大量试验数据确定m_{c}和s的值，使得公式能够较好地描述岩石的破坏条件。广义Hoek-Brown岩体强度准则公式\sigma_{1}=\sigma_{3}+\left(m_{b}\sigma_{c}\sigma_{3}+s\sigma_{c}^{2}\right)^{a}的推导，在最初公式的基础上，考虑了岩体的特征。引入参数a是为了更好地适应不同质量岩体的破坏特性，通过对不同岩体的试验研究和理论分析，确定a以及m_{b}、s与岩体性质的关系。例如，对于质量较好的岩体，m_{b}和s的取值相对较大，反映出岩体的完整性较好，强度较高；而对于破碎岩体，m_{b}和s的取值较小，体现出岩体的破碎程度对强度的削弱作用。在公式中，各参数具有明确的物理意义。\sigma_{1}和\sigma_{3}分别为最大和最小主应力，它们决定了岩体的应力状态。\sigma_{c}是岩块的单轴抗压强度，是岩石本身的固有属性，反映了岩石抵抗轴向压力的能力。m_{b}反映了岩体中岩石的软硬程度以及结构面对强度的影响，m_{b}值越大，说明岩石越坚硬，岩体的强度相对较高；s体现了岩体的破碎程度，s值越接近1，岩体越完整，强度越高，当s=0时，表示岩体极为破碎。参数a则调整了公式对不同岩体破坏特性的适应性，使得准则能够更准确地描述不同质量岩体的强度。2.3参数确定方法准确确定Hoek-Brown强度准则中的参数m、s等，是保证准则在岩体工程中准确应用的关键。以下介绍几种常见的参数确定方法。三轴试验法：对于完整岩块，一般先假定s=1，通过进行三轴压缩试验，得到不同围压下岩块破坏时的最大主应力\sigma_{1}和最小主应力\sigma_{3}数据。令x=\sigma_{3}和y=\sigma_{1}-\sigma_{3}，利用数据回归计算得到m值。对于破碎岩石，通常先利用已求得的完整岩石m、\sigma_{c}值，再根据试验数据计算s值。若计算结果s\lt0，则令s=0。例如，在某岩石三轴试验研究中，对花岗岩岩块进行不同围压下的三轴压缩试验，通过对试验数据的处理和回归分析，得到该花岗岩的m值，为后续岩体强度分析提供了重要参数。岩体RMR分类法：Hoek和Brown将m、s和Bieniawski的岩石分类指标RMR相联系，建议利用RMR直接确定m、s。对于未扰动岩体，m_{b}=m_{i}exp\left(\frac{RMR-100}{14}\right)；对于扰动岩体，m_{b}=m_{i}exp\left(\frac{RMR-100}{10}\right)，其中m_{i}为完整岩块的m值。通过对岩体进行现场调查，评估岩体的岩石强度、RQD值、节理间距、节理条件和地下水状况等因素，确定岩体的RMR值，进而计算m、s参数。在某边坡工程中，通过对边坡岩体的现场勘查和RMR分类，确定了岩体的RMR值，再根据上述公式计算出m、s参数，用于后续的边坡稳定性分析。地质强度指标（GSI）法：GSI是一种综合考虑岩体结构、岩块的嵌锁状态和岩体中不连续面质量等因素的岩体质量评价指标。当GSI\geq25时，m_{b}=m_{i}exp\left(\frac{GSI-100}{28-14D}\right)，s=exp\left(\frac{GSI-100}{9-3D}\right)；当GSI\lt25时，m_{b}=m_{i}exp\left(\frac{GSI-100}{9-3D}\right)，s=0，其中D为扰动系数。通过对岩体进行详细的地质调查，根据岩体结构特征、节理面的粗糙度、充填物等因素，确定岩体的GSI值，再结合扰动系数D计算m、s参数。在某地下洞室工程中，通过对洞室围岩的地质勘查和GSI评价，考虑到洞室开挖过程中的爆破扰动，确定了扰动系数D，进而准确计算出m、s参数，为洞室的支护设计提供了依据。断裂力学法：从断裂力学角度出发，通过研究岩石内部微裂纹的扩展和贯通机制，确定m、s参数与岩石断裂特性的关系。例如，利用断裂韧性试验，测量岩石的断裂韧性K_{IC}，通过理论推导和经验公式，建立K_{IC}与m、s之间的联系，从而确定m、s参数。在某岩石力学研究中，通过对砂岩进行断裂韧性试验，结合断裂力学理论，建立了m、s与断裂韧性的关系模型，为砂岩岩体强度分析提供了参数确定的新方法。三、在地下采矿工程中的应用3.1工程案例选取与概况湖北兴山后坪矿段磷矿层作为重要的地下磷矿资源开采区域，具有典型的地质特征和开采特点，为研究Hoek-Brown强度准则在地下采矿工程中的应用提供了理想案例。该矿段位于宜昌磷矿聚集区西部，属树崆坪矿区树崆坪矿段往北（深部）延伸地段，矿段面积达41.59km^{2}，大地构造处于扬子准地台北缘龙门—大巴台缘褶皱带的东端，北与东秦岭褶皱带的南端相连，处于次一级构造“黄陵背斜”北缘。从地层分布来看，矿段内出露地层较为复杂，涵盖第四系（Q）、寒武系下统石龙洞组（\in_{1}sh）、天河板组（\in_{1}t）、石牌组（\in_{1}sp）、牛蹄塘组（\in_{1}n）、震旦系上统灯影组（Z_{2}dn）第一至第三段，震旦系下统陡山沱组（Z_{1}d）第一至第四段、南华系上统南沱组（Nh_{2}n）及中元古界崆岭群（PtK）等地层，部分地层通过钻探和坑探工程可见。矿段内地层整体呈单斜层状产出，略有波状起伏，产状较为平缓，倾向285°～65°，倾角2°～13°之间，南部倾角2°～5°，北部倾角5°～13°。在构造方面，后坪矿段发现16条断层，均为正断层，依据走向可分为四组：北西西（NWW）向正断层、北北西（NNW）向正断层、近东西（EW）向正断层及北北东（NNE）向正断层。其中，北西西（NWW）向断层数量较多，对矿层破坏最为严重；北北西（NNW）向断层规模虽大，但数量相对较少，整体破坏作用稍逊于北西西向断层；近东西和北北东向断层规模较小，最大垂直断距一般小于20m，数量较少，对矿层破坏作用较小。这些断层均为张性正断层，断层面倾角较陡，一般大于70°，仅F_{5}断层在矿段中南部往深部变缓，局部仅47°。该矿段磷矿层赋存于震旦系陡山沱组第一段（Z_{1}d_{1}），是区域上Ph_{1}矿层的赋存层位。含磷岩系从上至下依次分布有上白云岩、Ph_{1}^{3}磷矿层、含磷泥岩、Ph_{1}^{2}矿层、含磷泥岩、Ph_{1}^{1}矿层层位、下白云岩。含磷岩系厚度在8.67～58.99m之间，平均厚27.52m。其中，Ph_{1}^{3}矿层最为发育，是矿段内主要工业磷矿层。在开采方式上，由于该矿段磷矿体呈缓倾斜产出，开采技术要求高，为保障开采安全和资源回收率，采用浅孔房柱法开采。在开采过程中，需预留矿柱来支撑顶板，以防止采空区围岩失稳。然而，这种开采方式下，采空区围岩的稳定性成为关键问题，受到岩体的地质条件、开采顺序、矿柱尺寸等多种因素影响。例如，若矿柱尺寸过小，无法有效支撑顶板，可能导致顶板垮落；而若矿柱尺寸过大，则会造成资源浪费。同时，断层等地质构造的存在，也会改变岩体的应力分布，增加采空区围岩失稳的风险。因此，准确评估采空区围岩的稳定性，对于保障矿山的安全开采至关重要。3.2基于准则的数值模拟分析为深入研究湖北兴山后坪矿段磷矿层开采过程中采空区的稳定性，采用有限差分软件FLAC3D，基于Hoek-Brown强度准则对开采过程进行数值模拟。FLAC3D软件能够较好地模拟岩土体的大变形和非线性力学行为，适用于复杂地质条件下的地下工程模拟。在数值模拟前，首先进行单轴压缩模拟试验，以求得各岩体的单轴抗压强度、弹性模量和泊松比等力学参数。通过模拟试验得到的结果与室内试验结果基本一致，验证了模拟方法的可靠性。在模拟过程中，依据Hoek-Brown强度准则确定岩体的强度参数。根据矿段的地质勘查资料，确定岩体的地质强度指标（GSI）值，考虑到开采过程中爆破等因素对岩体的扰动，确定扰动系数D。然后，利用公式m_{b}=m_{i}exp\left(\frac{GSI-100}{28-14D}\right)（当GSI\geq25时）或m_{b}=m_{i}exp\left(\frac{GSI-100}{9-3D}\right)（当GSI\lt25时）计算m_{b}值，再根据s=exp\left(\frac{GSI-100}{9-3D}\right)（当GSI\geq25时）或s=0（当GSI\lt25时）确定s值。在建立数值模型时，根据矿段的实际地质条件和开采设计，对模型的边界条件、初始应力场等进行合理设定。模型的边界条件采用位移约束，底部边界固定，四周边界限制水平位移。初始应力场根据上覆岩层厚度和岩体容重进行计算，考虑到自重应力和构造应力的影响。在模拟开采过程时，按照实际的开采顺序，逐步开挖磷矿层，并实时监测采空区围岩的应力、位移及矿柱屈服区分布情况。从模拟结果来看，在磷矿层开采过程中，采空区围岩的应力分布发生显著变化。随着开采的进行，采空区周边岩体的应力集中现象逐渐明显，尤其是在矿柱与围岩的交界处，应力值大幅增加。例如，在某一开采阶段，采空区顶板围岩的最大主应力达到[X]MPa，而在未开采区域，岩体的最大主应力仅为[Y]MPa。上覆岩体的位移也呈现出一定规律，在采空区上方，岩体产生向下的沉降位移，且位移量随着离采空区距离的增大而逐渐减小。在矿柱屈服区分布方面，模拟结果显示，部分矿柱的边角部位出现屈服现象，随着开采的持续，屈服区有逐渐扩大的趋势。当开采到一定程度时，矿柱屈服区的面积占矿柱总面积的[Z]%，这表明矿柱的承载能力受到一定程度的削弱。这些模拟结果为后续分析采空区的稳定性提供了重要依据。3.3应用效果与问题分析通过基于Hoek-Brown强度准则的数值模拟分析，得到的结果对于评价采空区稳定性和开采可行性具有重要作用。从采空区稳定性评价角度来看，根据模拟得到的采空区围岩应力、位移及矿柱屈服区分布情况，可以直观地了解采空区在开采过程中的稳定性状态。例如，当采空区围岩的应力超过岩体的强度时，可能会发生岩体破坏，导致顶板垮落等事故；而当矿柱屈服区过大时，矿柱可能无法有效支撑顶板，增加采空区失稳的风险。通过对这些参数的分析，可以及时发现潜在的安全隐患，为采取相应的支护措施提供依据。在开采可行性分析方面，模拟结果可以帮助判断当前的开采方案是否合理。如果在模拟过程中发现开采过程中采空区稳定性难以保证，或者资源回收率较低等问题，可以对开采方案进行优化，如调整矿柱尺寸、改变开采顺序等，以提高开采的可行性和安全性。然而，在应用Hoek-Brown强度准则进行数值模拟过程中，也遇到一些问题。首先，准则中参数的确定存在一定困难。虽然有多种参数确定方法，但不同方法得到的参数值可能存在差异。例如，通过三轴试验法和地质强度指标（GSI）法确定的m、s参数，在某些情况下会有较大偏差。这是因为三轴试验法主要基于实验室岩块试验，难以完全反映现场岩体的复杂地质条件；而GSI法虽然考虑了岩体的结构和质量等因素，但在确定GSI值时，存在一定的主观性。其次，对于复杂地质条件下的岩体，如含有多种结构面、强烈各向异性的岩体，Hoek-Brown强度准则的适用性有待进一步验证。在该矿段中，由于存在多条断层和节理裂隙，岩体的各向异性较为明显，准则在描述这种岩体的力学行为时，可能存在一定误差。针对这些问题，可以采取以下解决方法。在参数确定方面，可以综合多种方法来确定参数值。例如，先通过现场试验获取初步的参数值，再结合地质勘查资料，利用GSI法等进行修正，以提高参数的准确性。同时，随着机器学习和人工智能技术的发展，可以利用这些技术对大量的岩体试验数据和地质资料进行分析，建立更准确的参数预测模型。对于复杂地质条件下岩体的适用性问题，可以对准则进行改进，考虑更多的影响因素。例如，在准则中加入反映结构面方向、间距和粗糙度等因素的参数，以更好地描述各向异性岩体的力学行为。此外，还可以结合其他强度准则或数值模拟方法，如离散元法等，对复杂地质条件下的岩体进行分析，相互验证结果，提高分析的可靠性。四、在隧道工程中的应用4.1隧道工程实例介绍龙南隧道1号斜井作为赣深高铁龙南隧道的重要施工通道，其建设对于整个隧道工程的顺利推进至关重要。龙南隧道位于江西省赣州市龙南和全南县境内，是赣深高铁重难点控制性工程，线路全长10.24km，最大埋深580m，最小埋深12m，双线单洞。1号斜井承担龙南隧道正洞3624m（占隧道全长近三分之一）的施工任务，施工任务异常繁重。从地质条件来看，龙南隧道1号斜井区域地质情况极为复杂。隧道共穿越11条断层，其中有长度120m的大断层破碎带，1号斜井在施工过程中不可避免地受到这些断层的影响。断层的存在使得岩体破碎，完整性遭到严重破坏，节理裂隙发育，导致岩体的力学性质变得极为复杂且不稳定。例如，在断层破碎带附近，岩体的强度大幅降低，其单轴抗压强度可能仅为正常岩体的[X]%。同时，该区域洞身存在危岩落石、突泥涌水、岩爆、高地温、岩溶、浅埋、围岩失稳塌方等诸多风险。在高地温地段，洞内温度常年高达36-40度，给施工人员和设备带来极大挑战，高温环境不仅影响施工人员的身体健康和工作效率，还可能导致设备故障频发。而在岩溶发育区域，溶洞和溶蚀裂隙的存在增加了突泥涌水的风险，一旦发生突泥涌水事故，可能会淹没隧道，造成严重的人员伤亡和财产损失。在施工难点方面，首先，由于斜井井身长度大、正洞施工长，导致1号斜井施工区域作业环境恶劣。各类灰尘、柴油废气等形成环形带，一直不能有效排出，严重影响施工人员的呼吸系统健康。洞内高温高湿，在遇到高温时，水变成水汽，与各类灰尘等结合，使得视线差、胸闷及窒息感强烈，极大地降低了施工人员的工作效率和安全性。其次，复杂的地质条件给施工带来了技术难题。在穿越断层破碎带时，如何保证围岩的稳定性成为关键问题。传统的施工方法难以满足要求，需要采用特殊的支护措施和施工工艺。例如，在某断层破碎带施工时，采用常规的锚杆支护无法有效控制围岩变形，不得不增加锚索支护，并采用超前小导管注浆加固等措施，才确保了施工的安全进行。此外，隧道施工过程中的爆破作业对围岩稳定性也产生了较大影响。爆破产生的振动和应力波可能导致围岩进一步破碎，增加了施工风险。4.2考虑爆破损伤的准则修正在龙南隧道1号斜井的施工过程中，爆破作为主要的开挖方式，不可避免地对岩体造成损伤。爆破产生的应力波在岩体中传播，导致岩体内部产生微裂纹，随着爆破次数的增加，这些微裂纹逐渐扩展、贯通，形成宏观裂缝，从而降低岩体的强度和稳定性。例如，通过现场声波检测发现，在爆破区域附近，岩体的纵波波速明显降低，降低幅度可达[X]%，这表明岩体的完整性受到了破坏。为了准确评估爆破损伤对岩体强度的影响，提出基于爆破累积损伤效应修正Hoek-Brown强度准则经验常数的方法。该方法考虑了爆破次数、爆破药量、岩体特性等因素对损伤的影响。引入爆破累积损伤因子D_{b}，其表达式为：D_{b}=1-exp\left(-\alphaNQ^{\beta}\right)其中，N为爆破次数，Q为单段最大爆破药量，\alpha和\beta为与岩体性质相关的参数，可通过现场试验或经验公式确定。根据爆破累积损伤因子D_{b}，对Hoek-Brown强度准则中的经验常数m_{b}和s进行修正。修正后的m_{b}和s表达式为：m_{b}^{*}=m_{b}\left(1-D_{b}\right)^{\gamma}s^{*}=s\left(1-D_{b}\right)^{\delta}其中，m_{b}^{*}和s^{*}为修正后的经验常数，\gamma和\delta为与岩体性质和损伤程度相关的参数，可通过试验研究确定。一般来说，\gamma和\delta的值越大，说明爆破损伤对经验常数的影响越显著。通过这种基于爆破累积损伤效应修正Hoek-Brown强度准则经验常数的方法，可以更准确地考虑爆破对岩体强度的影响，为隧道施工过程中的围岩稳定性分析提供更可靠的依据。在龙南隧道1号斜井的实际应用中，利用该修正方法对不同施工阶段的岩体强度进行评估，与未修正前的结果相比，更符合现场实际情况。例如，在某施工段，未修正前预测的围岩变形量与实际监测值存在较大偏差，而采用修正后的准则进行计算，预测的围岩变形量与实际监测值基本吻合，误差在可接受范围内，从而为施工决策提供了更准确的参考。4.3围岩参数反分析与工程指导为了准确获取龙南隧道1号斜井围岩的力学参数，利用遗传—粒子群反分析程序进行围岩力学参数反演。遗传算法具有全局搜索能力强、鲁棒性好等优点，能够在复杂的参数空间中寻找最优解。粒子群算法则具有收敛速度快、计算简单等特点，能够快速逼近最优解。将两者结合，充分发挥各自的优势，提高了反演的效率和准确性。在反演过程中，以现场监测的围岩位移数据为依据，建立目标函数。目标函数定义为：f(X)=\sum_{i=1}^{n}\left(u_{i}^{cal}(X)-u_{i}^{obs}\right)^{2}其中，X为待反演的围岩力学参数向量，包括弹性模量E、泊松比\mu、岩石单轴抗压强度\sigma_{c}、地质强度指标GSI等；u_{i}^{cal}(X)为根据参数X计算得到的第i个监测点的围岩位移值；u_{i}^{obs}为第i个监测点的实际监测位移值；n为监测点的数量。通过遗传—粒子群反分析程序不断迭代计算，调整围岩力学参数，使得目标函数f(X)达到最小值，此时的参数X即为反演得到的围岩力学参数。根据反演得到的围岩力学参数，利用数值模拟软件对隧道开挖过程进行模拟，预测围岩变形及应力状态。在模拟过程中，考虑了隧道的施工顺序、支护措施等因素。例如，模拟结果显示，在某一施工阶段，隧道拱顶的最大沉降量为[X]mm，边墙的最大水平位移为[Y]mm，通过与现场监测数据对比，验证了模拟结果的准确性。基于预测结果，为隧道施工提供科学指导。当预测到围岩变形量超过允许范围时，及时调整施工方案，加强支护措施。如增加锚杆长度和密度、加大喷射混凝土厚度等。在某施工段，预测到围岩变形较大，通过增加锚杆长度和密度后，现场监测数据表明围岩变形得到了有效控制，确保了施工安全。同时，根据预测的应力状态，合理安排施工进度，避免在高应力区进行大规模开挖，减少了施工风险。五、在边坡工程中的应用5.1边坡工程案例分析清江隔河岩水电站厂房开挖高边坡是一个极具代表性的工程案例，其在水电工程建设中具有重要地位。该边坡前沿总长约350m，正面坡高110m，侧面坡高170m，是一个大型的人工高边坡。其上部为灰岩直立边坡，下部为页岩斜坡，按多级直立台阶式布置。这种特殊的地质结构使得边坡的稳定性面临诸多挑战。从地质条件来看，该边坡为逆向坡，上部为灰岩（上硬）、下部为页岩（下软）的不利地层结构，岩层走向70°～80°，倾南东，倾角25°～30°。虽然这种逆向坡结构在一定程度上对边坡稳定性有利，但10余条断层、夹层，4组裂隙和多组裂隙断层切割，以及两个危岩体及岸溶塌陷体等地质缺陷的存在，使得局部地区岩体较破碎，极大地影响了边坡的整体稳定性。例如，断层的存在破坏了岩体的完整性，降低了岩体的强度，使得边坡在受到外部荷载作用时，更容易发生变形和破坏。在开挖情况方面，为满足边坡工程稳定性的需要，施工时采取了预裂光面爆破方式。将高程240、195、180、165m减载开挖为直立边坡，高程165m平台以下150、128、120、110、100m为斜坡段。上部边坡岩体为灰岩，下部为页岩，斜坡段高差为10m，平台宽2～3m。这种开挖方式虽然有利于整体边坡稳定，但爆破施工不可避免地对岩体造成一定损伤，进一步影响了岩体的力学性质。通过现场声波检测发现，爆破区域附近岩体的纵波波速明显降低，降低幅度可达[X]%，这表明岩体的完整性受到了破坏，强度有所降低。在边坡加固措施上，采取了开挖减载和锚索加固等方法。通过开挖减载，将自然边坡开挖达到多级台阶式人工边坡，减少了边坡的荷载。预应力锚索加固则对灰岩中岸剪裂隙按抗倾要求作了相应加固处理，主动限制了拉裂作用，且群体补偿作用明显。在某段边坡加固中，通过布置预应力锚索，有效地控制了边坡岩体的变形，使得该段边坡在后续的运行中保持了稳定。5.2与其他准则对比分析为了更全面地评估Hoek-Brown准则在清江隔河岩水电站厂房开挖高边坡工程中的适用性，将其与Mohr-Coulomb准则进行对比分析。在模拟卸荷岩体力学性质方面，两者存在显著差异。Mohr-Coulomb准则是一种基于线性假设的强度准则，它假设岩石的抗剪强度与正应力呈线性关系，表达式为\tau=c+\sigma\tan\varphi，其中\tau为抗剪强度，c为黏聚力，\sigma为正应力，\varphi为内摩擦角。该准则在处理简单地质条件下的岩体力学问题时具有一定的优势，计算相对简便。然而，在面对复杂地质条件下的岩体，如清江隔河岩水电站厂房开挖高边坡这种存在大量断层、裂隙和破碎带的岩体时，其局限性就凸显出来。由于该准则没有考虑岩体的非线性破坏特性以及中间主应力的影响，在模拟卸荷岩体力学性质时，往往会高估岩体的强度，导致计算结果与实际情况存在较大偏差。相比之下，Hoek-Brown准则能够更好地模拟卸荷岩体力学性质。该准则充分考虑了岩石和岩体固有的非线性破坏特点，通过引入反映岩体特征的经验参数，如m_{b}、s和a等，能够更准确地描述岩体在不同应力状态下的强度特性。在清江隔河岩水电站厂房开挖高边坡工程中，Hoek-Brown准则能够考虑到断层、裂隙等地质缺陷对岩体强度的影响，更真实地反映岩体的力学行为。通过数值模拟分析发现，采用Hoek-Brown准则计算得到的边坡位移和应力分布与现场监测数据更为吻合。例如，在某一工况下，采用Mohr-Coulomb准则计算得到的边坡最大位移为[X]mm，而实际监测值为[Y]mm；采用Hoek-Brown准则计算得到的最大位移为[Z]mm，与实际监测值更为接近，误差在可接受范围内。这表明Hoek-Brown准则在模拟复杂地质条件下卸荷岩体力学性质方面具有明显优势，能够为边坡工程的设计和施工提供更可靠的依据。5.3边坡稳定性评价与预测基于Hoek-Brown强度准则对清江隔河岩水电站厂房开挖高边坡的稳定性进行评价，具有重要的工程意义。通过该准则，可以准确分析边坡在不同工况下的稳定性变化。在正常工况下，根据边坡的地质条件、岩体力学参数以及荷载情况，利用Hoek-Brown强度准则进行计算分析。结果表明，边坡整体处于稳定状态，但在局部区域，如断层附近和边坡下部页岩层，存在一定的潜在不稳定因素。这些区域的岩体由于受到地质构造和岩石性质的影响，强度相对较低，在长期的荷载作用下，可能会发生变形和破坏。通过数值模拟可以直观地看到，在这些区域，岩体的应力集中现象较为明显，位移也相对较大。当考虑地震、暴雨等特殊工况时，边坡的稳定性会受到更大的挑战。地震作用会使边坡岩体受到强烈的震动，增加岩体的动应力，导致岩体的强度降低。暴雨则会使边坡岩体的含水量增加，岩体的重度增大，同时孔隙水压力也会升高，进一步降低岩体的抗剪强度。利用Hoek-Brown强度准则，结合地震和暴雨的相关参数，对边坡在这些特殊工况下的稳定性进行预测。模拟结果显示，在地震和暴雨的共同作用下，边坡的潜在不稳定区域范围扩大，部分区域可能会发生滑坡等地质灾害。例如，在某一地震烈度和暴雨强度下，边坡的某一段可能会出现沿软弱结构面滑动的情况，滑动体的体积约为[X]立方米。针对不同工况下边坡的稳定性分析结果，提出相应的加固措施。对于正常工况下存在潜在不稳定因素的区域，可以采用锚杆、锚索等支护措施，增强岩体的整体性和强度。在断层附近，可以加密锚杆和锚索的布置，提高岩体的锚固力。对于页岩层等软弱部位，可以采用注浆加固的方法，填充岩体的孔隙和裂隙，提高岩体的抗剪强度。在特殊工况下，除了加强上述支护措施外，还需要采取一些特殊的防护措施。例如，在边坡顶部设置截水沟，防止雨水大量流入边坡；在坡脚设置挡土墙，增加边坡的抗滑力。通过这些加固措施的实施，可以有效提高边坡在不同工况下的稳定性，确保水电站厂房的安全运行。六、应用优势与局限性分析6.1应用优势探讨Hoek-Brown强度准则在岩体工程应用中展现出多方面显著优势。首先，该准则能全面考虑多种岩石破坏机制。岩石的破坏是一个复杂过程，涉及微裂纹的产生、扩展和贯通等多种机制。Hoek-Brown强度准则通过引入反映岩石软硬程度的参数m和反映岩体破碎程度的参数s等，充分考虑了岩石内部结构特征、节理裂隙分布以及岩石本身的力学性质对破坏机制的影响。例如，在含有大量节理裂隙的岩体中，参数s的取值会较小，反映出岩体的破碎程度高，强度较低，从而更准确地描述这种岩体在受力时的破坏过程。其次，准则适用范围广泛。无论是完整的岩石还是节理发育、破碎的岩体，Hoek-Brown强度准则都能适用。在地下采矿工程中，面对不同地质条件下的矿体和围岩，该准则可用于分析采场的稳定性；在隧道工程里，能对不同围岩条件下的隧道进行稳定性评估；在边坡工程中，对于各种复杂地质结构的边坡，也能利用该准则进行稳定性分析。如在某复杂地质条件的隧道工程中，通过Hoek-Brown强度准则准确评估了围岩的稳定性，为隧道的支护设计提供了科学依据。再者，Hoek-Brown强度准则使用相对方便。其基本公式形式简洁，虽然参数确定可能需要一定的试验和分析，但一旦确定参数，在实际工程计算中应用较为便捷。与一些复杂的强度准则相比，不需要进行大量复杂的数学推导和计算，降低了工程应用的难度。在某小型矿山的开采设计中，工程人员利用Hoek-Brown强度准则，通过简单的参数确定和计算，快速评估了采场的稳定性，提高了设计效率。最后，该准则在一定程度上具有较高的精度。通过与大量的岩石三轴试验资料和岩体现场试验成果对比验证，Hoek-Brown强度准则能够较好地反映岩体在不同应力状态下的强度特性，计算结果与实际情况较为吻合。在某大型水电工程的边坡稳定性分析中，利用该准则计算得到的边坡稳定性系数与现场监测结果基本一致，误差在可接受范围内，为工程的安全运行提供了可靠保障。6.2局限性分析尽管Hoek-Brown强度准则具有诸多优势，但在实际应用中也存在一些局限性。首先，准则中常数取值依赖实测数据。准确确定参数m、s等对于准则的应用至关重要，但这些参数的取值往往需要通过现场试验、三轴试验等实测手段来确定。在实际工程中，进行大量的现场试验和三轴试验不仅成本高昂，而且耗时费力。例如，在一些偏远地区的岩体工程中，进行现场试验可能面临交通不便、试验设备难以运输等困难。同时，不同的试验方法和试验条件可能导致参数取值存在差异，增加了参数确定的不确定性。其次，Hoek-Brown强度准则对岩石破坏上限假设存在局限性。该准则假设岩石破坏时的强度上限是固定的，然而在实际情况中，岩石的破坏上限可能受到多种因素的影响，如温度、加载速率等。在高温环境下，岩石的力学性质会发生变化，其破坏上限也会相应改变，而Hoek-Brown强度准则无法准确反映这种变化。在一些深部岩体工程中，由于地温较高，岩石的强度特性与常温下有很大不同，此时该准则的应用就存在一定误差。再者，准则无法考虑任意曲面强度。Hoek-Brown强度准则主要基于主应力空间来描述岩石的强度特性，对于一些复杂的受力情况，如岩体中存在任意方向的结构面时，该准则难以准确考虑结构面的强度和方向对岩体整体强度的影响。在某岩体工程中，存在一组倾斜角度较大的节理面，传统的Hoek-Brown强度准则在分析该岩体的强度时，无法充分考虑节理面的影响，导致计算结果与实际情况存在偏差。6.3改进方向与展望针对Hoek-Brown强度准则存在的局限性，可从以下几个方面进行改进