CN111551427B - 一种深埋长隧洞软质岩大变形超前量化预报方法 - Google Patents

一种深埋长隧洞软质岩大变形超前量化预报方法 Download PDF

Info

Publication number
CN111551427B
CN111551427B CN202010328971.4A CN202010328971A CN111551427B CN 111551427 B CN111551427 B CN 111551427B CN 202010328971 A CN202010328971 A CN 202010328971A CN 111551427 B CN111551427 B CN 111551427B
Authority
CN
China
Prior art keywords
rock
deformation
tunnel
stress
rock mass
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010328971.4A
Other languages
English (en)
Other versions
CN111551427A (zh
Inventor
黄书岭
丁秀丽
王旭一
张雨霆
何军
付敬
秦洋
赵化蒙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changjiang River Scientific Research Institute Changjiang Water Resources Commission
Original Assignee
Changjiang River Scientific Research Institute Changjiang Water Resources Commission
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Changjiang River Scientific Research Institute Changjiang Water Resources Commission filed Critical Changjiang River Scientific Research Institute Changjiang Water Resources Commission
Priority to CN202010328971.4A priority Critical patent/CN111551427B/zh
Publication of CN111551427A publication Critical patent/CN111551427A/zh
Application granted granted Critical
Publication of CN111551427B publication Critical patent/CN111551427B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0058—Kind of property studied
    • G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
    • G01N2203/0075—Strain-stress relations or elastic constants
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02—Details not specific for a particular testing method
    • G01N2203/0202—Control of the test
    • G01N2203/0212—Theories, calculations
    • G01N2203/0218—Calculations based on experimental data
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02—Details not specific for a particular testing method
    • G01N2203/022—Environment of the test
    • G01N2203/0244—Tests performed "in situ" or after "in situ" use
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02—Details not specific for a particular testing method
    • G01N2203/025—Geometry of the test
    • G01N2203/0252—Monoaxial, i.e. the forces being applied along a single axis of the specimen
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02—Details not specific for a particular testing method
    • G01N2203/025—Geometry of the test
    • G01N2203/0256—Triaxial, i.e. the forces being applied along three normal axes of the specimen
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02—Details not specific for a particular testing method
    • G01N2203/06—Indicating or recording means; Sensing means
    • G01N2203/067—Parameter measured for estimating the property
    • G01N2203/0676—Force, weight, load, energy, speed or acceleration
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02—Details not specific for a particular testing method
    • G01N2203/06—Indicating or recording means; Sensing means
    • G01N2203/067—Parameter measured for estimating the property
    • G01N2203/0682—Spatial dimension, e.g. length, area, angle

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种深埋长隧洞软质岩大变形超前量化预报方法，在施工期，对隧洞掌子面前方的可能出现的软质岩体地层，考虑隧洞掌子面前方远、中、近距离处的量化预报；通过对围岩大变形判别、围岩变形预测，获得直观的判断和预测数据，定量化预报深埋长隧(洞)道软质岩的变形量级及相应的大变形等级，为深埋长隧洞软质岩支护优化设计和施工方案确定提供依据。

Description

一种深埋长隧洞软质岩大变形超前量化预报方法
技术领域
本发明涉及一种隧洞软质围岩大变形的量化预报方法，具体涉及一种深埋长隧洞软质岩大变形超前量化预报方法，在施工期对隧洞掌子面前方的可能出现的软质岩体地层，开展围岩大变形的超前量化预报工作。
背景技术
对于穿山越岭的交通隧道、引调水隧洞等线路工程，线路长，埋深大，通常要穿越岩性和地质条件复杂多变的地层，特别是当隧(道) 洞遭遇泥质岩、蚀变岩、断层挤压带等软弱岩体时，时常面临施工期软岩大变形、支护结构破坏、围岩加固难度大等工程难题。通常，线路工程施工前有限的地勘工作尚不能完全覆盖隧(道)洞穿越的地层信息，因此，在施工期对隧(道)洞掌子面前方的地层，开展软质岩的超前量化预报工作是非常必要的，具有重要的工程意义。
隧(道)洞围岩变形的量化预报一直是地下工程的一项重要课题。开展掌子面前方不良地质体的探测是隧(洞)道工程建设施工过程中必要的环节。一般而言，中远距离(掌子面前方100m)探测多基于地震法或电磁法的超前地质预报原理，需要具有地球物理专业背景的工程师对探测结果进行解读，从而给出掌子面前方是否存在不良地质体的大致判断；近距离(掌子面前方60m内)探测多采用超前钻孔方法，通过钻孔取芯直观获得掌子面前方地质体的直观描述。然而，采用通过中远距离和近距离探测获得的上述前方地质体信息，仅能为工程师提供掌子面前方不良地质体的位置和性状等定性描述，难以用于围岩变形的定量预测预报，也就无法为支护措施的动态调整提供量化依据。因此，建立一种深埋长隧洞软质岩大变形超前量化预报方法，对于隧(道)洞工程保障施工期围岩稳定性具有重要意义。
发明内容
本发明针对上述存在的问题，提供一种深埋长隧洞软质岩大变形超前量化预报方法，在施工期对隧洞掌子面前方的地层，开展软质岩的超前预测预报工作。
为实现上述目的，本发明所设计的一种深埋长隧洞软质岩大变形超前量化预报方法，其特殊之处在于，所述方法包括：
步骤一：中远距离超前地层信息预报技术：采用隧洞超前地质预报***，对掌子面前方设定距离的地层信息进行探测并获得三维空间数据，利用集面向软岩不同尺度变形质构造体的隧洞反射波正演模拟技术、波场分离和偏移成像的数据处理技术、地震数据保幅成像方法、基于深度学习的非线性智能反演方法于一体的信号解译技术，解译并预报隧洞掌子面前方地层的性质、位置及规模大小；
步骤二：近距离超前地质钻探预报技术：在步骤一工作的基础上，当隧洞开挖至软质岩30～60m范围时，在掌子面打超前钻孔；
步骤三：在步骤二钻孔部位进行原位地应力测试，获取该处三维初始应力分布状态、最大水平侧压力系数λ及初始最大主应力值σ0max参数；
步骤四：在步骤二钻孔部位进行岩石取样，开展室内物理力学试验，准确掌握掌子面前方软质岩石物理力学特性；
步骤五：在步骤二钻孔部位进行钻孔电视工作，获得软岩裂隙的埋深、倾向、倾角、宽度、裂隙面的粗糙度、充填物的性质，在步骤四岩石声波测试取样部位进行现场声波测试工作，获得岩体弹性纵波速度Vpm，结合步骤四室内试验获得的岩石弹性纵波速度Vpr，计算岩体完整性系数
Figure BDA0002464250820000021
步骤六：根据室内试验和现场试验，确定岩体基本质量分级，进一步确定掌子面前方软质岩体的物理力学参数：岩体基本质量分级，根据岩体基本质量的定性特征和岩体基本质量指标BQ两者相结合确定，根据岩体的基本质量级别确定岩体黏聚力c，岩体内摩擦角φ和岩体变形模量E；
步骤七：隧洞围岩大变形预测预报技术：计算自重应力场或者构造应力场下软质岩的强度应力比，利用围岩大变形判别方法，对围岩是否会发生挤压大变形进行判断和预测；
步骤八：在步骤七的基础上，若判别为发生围岩大变形，利用考虑支护效应修正的围岩变形和相对变形量化预报公式，通过最速下降法，快速量化预报围岩变形量级和相对变形，给出大变形等级，为软质岩支护优化设计和施工方案提供依据。
优选地，所述步骤三中原位地应力测试，测试仪器采用三向应变计，测试方法采用应力解除法，并向孔内添加水泥砂浆作为耦合介质。
优选地，所述步骤四中岩石室内试验包括：进行岩石饱和单轴压缩试验，获得岩石饱和单轴抗压强度Rc和泊松比μ；进行岩石饱和三轴压缩试验，获得扩容梯度η；进行室内岩石声波测试，获得岩石弹性纵波速度Vpr。
优选地，所述步骤六中岩体基本质量指标BQ的计算公式为 BQ＝100+3Rc+250Kv，使用该式时，遵循以下限制条件：(1)当Rc＞90Kv+30 时，以Rc＝90Kv+30和Kv带入计算BQ值；(2)当Kv＞0.04Rc+0.4时，以 Kv＝0.04Rc+0.4和Rc代入计算BQ值；式中Rc为岩石单轴饱和抗压强度； Kv为岩体完整性系数。
优选地，所述步骤七中自重应力场(λ≤1)或者构造应力场(λ＞1) 下的岩石强度应力比值公式分别为
Figure BDA0002464250820000031
式中，Rc为岩石饱和单轴抗压强度；σ0max为初始最大主应力值。
优选地，所述步骤八中隧洞软质岩洞周变形量化预报，由公式
Figure BDA0002464250820000041
量化计算获得；式中，R′为塑性半径，计算公式为
Figure BDA0002464250820000042
r0为隧洞半径；σc为岩体单轴抗压强度，计算公式为
Figure BDA0002464250820000043
Figure BDA0002464250820000044
c为岩体的黏聚力，
Figure BDA0002464250820000045
为岩体的内摩擦角；其中
Figure BDA0002464250820000046
Figure BDA0002464250820000047
P为铅直向初始地应力；λ为最大水平侧压力系数；η为扩容梯度；反映围岩不同位置的θ角；E为岩体的变形模量；μ为泊松比；Pi为支护抗力。
优选地，所述步骤八中考虑多种支护形式时，支护抗力Pi由公式
Figure BDA0002464250820000048
量化确定；其中Ec是混凝土的弹性模量；R是支护的最外侧半径；S1为钢拱架沿隧洞轴线分布的间距；S2为钢筋沿隧洞轴线分布的间距；Sc为锚杆沿隧洞圆周分布的间距；Sl为锚杆沿隧洞轴线分布的间距；As1为钢拱架的截面面积；As2为隧洞断面中钢筋的截面面积；Es1为钢的弹性模量；Es2为锚杆锚索的弹性模量；vc是混凝土的泊松比；tc是混凝土的厚度；tB为支护钢拱架圈的厚度；B为槽钢或角钢的翼缘宽度；Is为隧洞截面的惯性矩；θ为两个阻块节点夹角的一半弧度；EB为阻块材料的弹性模量；d为超挖区填充混凝土的平均厚度；其中db为锚杆、锚索的直径；l为锚杆、锚索的自由长度；Q为锚杆端部和头部的变形荷载常量。
优选地，所述步骤七中围岩大变形判别方法为对相对变形ε量化计算获得，当相对变形2.5％＜ε＜5％时，为I级围岩大变形，变形程度为中等；当相对变形5％＜ε＜10％时，为II级围岩大变形，变形程度为严重；当相对变形ε＞10％时，为III级围岩大变形，变形程度为极严重。
由于采用了以上技术方案，本发明同时考虑隧洞掌子面前方远、中、近距离处的量化预报工作，是一种深埋长隧洞软质岩大变形超前预测预报方法；通过对围岩大变形判别、围岩变形预测，获得直观的判断和预测数据，定量化预报深埋长隧(洞)道软质岩的变形量级和等级，有效开展软质岩的超前预测预报工作。
附图说明
图1为本发明一种深埋长隧洞软质岩大变形超前量化预报方法的流程图。
图2为本发明的混凝土支护衬砌示意图。
图3为采用本发明获得的不同埋深条件下初始铅直地应力为 4～30MPa时隧洞围岩变形量级和相对变形以及相对应的大变形等级示意图。
具体实施方式
以下结合附图和具体实施例对本发明作进一步的详细描述。
如图1所示，本发明提出的一种深埋长隧洞软质岩大变形超前量化预报方法，包括以下步骤：
步骤一：中远距离超前地层信息预报技术。选用TRT7000隧洞超前地质预报***或者长江科学院研发的TEP隧洞超前地质预报***，对掌子面前方中远距离(100m)的地层信息进行探测并获得三维空间数据，利用集面向软岩不同尺度变形质构造体的隧洞反射波正演模拟技术、波场分离和偏移成像的数据处理技术、地震数据保幅成像方法、基于深度学习的非线性智能反演方法于一体的信号解译技术，解译并预报隧洞掌子面前方地层的性质(如节理裂隙带、软弱带、断层破碎带等)、位置及规模大小。
步骤二：在TEP***反馈分析辅助下进行隧洞开挖，当隧洞开挖至掌子面距软质岩30～60m范围时，采用近距离超前地质钻探预报技术，在掌子面打超前钻孔，一般在掌子面上分散布设3～6个钻孔，沿隧洞洞轴线方向钻进。通过钻进揭露来判定掌子面前方的地质情况，具体地：通过钻孔是否出水和出水量可判断前方地下水赋存情况；钻进速度的高低可以反映前方软质岩的坚硬程度；岩芯采取率可反映前方裂隙发育情况等。
步骤三：在步骤二钻孔部位进行原位地应力测试。测试仪器可采用三向应变计，测试方法可采用应力解除法，并向孔内添加水泥砂浆作为耦合介质。获取该处三维初始应力分布状态、最大水平侧压力系数λ及初始最大主应力值σ0max，并判断当λ≤1时地下洞室的应力环境为自重应力场，λ＞1时地下洞室的应力环境为构造应力场。
步骤四：在步骤二钻孔部位进行岩石取样，开展室内物理力学试验，准确掌握掌子面前方软质岩石物理力学特性。现场取样后将岩芯加工成高径比为1∶2的圆柱形试件，取直径为50mm、高度为100mm。具体的，进行室内岩石声波测试，获得岩石弹性纵波速度Vpr；饱和处理后进行岩石饱和单轴压缩试验，获得岩石饱和单轴抗压强度Rc泊松比μ；饱和处理后进行岩石饱和三轴压缩试验，获得扩容梯度η。
步骤五：在步骤二钻孔部位进行钻孔电视和声波测试工作。先向孔内添加水作为耦合介质，利用钻孔电视自孔口至孔底获取钻孔壁360°图像，数字化处理后获得软岩裂隙的埋深、倾向、倾角、宽度、裂隙面的粗糙度、充填物等性质。声波测试采用单发单收的跨孔声波仪(或单发双收的单孔声波仪)，软岩中宜采用干孔换能器，在步骤四岩石声波测试取样部位进行现场声波测试工作，获得岩体弹性纵波速度Vpm，结合步骤四室内试验获得的岩石弹性纵波速度Vpr，计算岩体完整性系数
Figure BDA0002464250820000071
步骤六：根据室内试验和现场试验，确定岩体基本质量分级，进一步确定掌子面前方软质岩体的物理力学参数。
岩体基本质量分级应根据岩体基本质量的定性特征和岩体基本质量指标(BQ)两者相结合确定。岩体基本质量的定性特征根据现场评价的岩石坚硬程度和岩体完整程度确定；岩体基本质量指标(BQ) 计算公式为BQ＝100+3Rc+250Kv，使用该式时，应遵循以下限制条件：(1) 当Rc＞90Kv+30时，应以Rc＝90Kv+30和Kv带入计算BQ值；(2)当 Kv＞0.04Rc+0.4时，应以Kv＝0.04Rc+0.4和Rc代入计算BQ值。式中Rc为岩石单轴饱和抗压强度即σc；Kv为岩体完整性系数。
根据岩体的基本质量级别确定岩体黏聚力c，岩体内摩擦角
Figure BDA0002464250820000072
和岩体变形模量E。
步骤七：隧洞围岩大变形量化预报技术。在步骤三取得原位最大主应力值σ0max和步骤四取得圆柱体岩样的岩石饱和单轴抗压强度Rc后，可计算自重应力场或者构造应力场下软质岩的强度应力比。自重应力场(λ≤1)或者构造应力场(λ＞1)下的岩石强度应力比值
Figure BDA0002464250820000081
然后采用
Figure BDA0002464250820000082
判断大埋深洞室是否会发生围岩挤压大变形。当SSR大于1时，判别为不会发生围岩挤压大变形，否则判别为会发生围岩挤压大变形。判别为会发生围岩挤压大变形，则需预测预报围岩变形量级；判别为不会发生围岩挤压大变形，则优化支护设计、降低支护强度。
步骤八：在步骤七的基础上，若判别为发生围岩大变形，利用考虑支护效应修正的围岩变形和相对变形量化预报公式，通过最速下降法，快速量化预报围岩变形量级和相对变形，给出大变形等级，为软质岩支护优化设计和施工方案提供依据。
隧洞软质岩洞周变形预测，可由公式
Figure BDA0002464250820000083
计算获得，该式推导过程如下：
假定塑性区中的应变呈轴对称分布，则塑性区几何方程：
Figure BDA0002464250820000084
Figure BDA0002464250820000085
考虑围岩屈服时的扩容效应，并假设扩容梯度η为定值。则根据流动法则ηεθp+εrp＝0，因此有：
Figure BDA0002464250820000086
再利用弹塑***界面处的位移连续条件Urp|r＝R＝Ure|r＝R，即可得解。对于有支护时，隧洞洞周位移求解过程如下：
弹性区应力分量：
Figure BDA0002464250820000091
弹塑***界处径向应力：
Figure BDA0002464250820000092
塑性半径：
Figure BDA0002464250820000093
弹性区的弹性径向位移Ur和环向位移Uθ分别为：
Figure BDA0002464250820000094
Figure BDA0002464250820000095
弹塑***界处弹性位移：
Figure BDA0002464250820000096
塑性区位移：
Figure BDA0002464250820000097
因此在隧洞洞周总的径向位移为：
Figure BDA0002464250820000098
式中，r0为隧洞半径；σc为岩体单轴抗压强度，计算公式为
Figure BDA0002464250820000099
Figure BDA00024642508200000910
c为岩体黏聚力，
Figure BDA00024642508200000911
为岩体的内摩擦角；
Figure BDA00024642508200000912
Figure BDA0002464250820000101
P为铅直向初始地应力；λ为最大水平侧压力系数；η为扩容梯度；反映围岩不同位置的θ角；E为岩体的变形模量；μ为泊松比；Pi为支护抗力。
考虑多种支护形式时，支护抗力Pi由公式Pi＝(Ks1+Ks2+Ks3+Ks4)ur确定，式中Ks1为喷射混凝土的支护刚度系数；Ks2为钢拱架的支护刚度系数；Ks3为钢筋混凝土衬砌的支护刚度系数；Ks4为无粘结锚杆的支护刚度系数。各支护形式的支护刚度系数计算公式如下：
对于喷射混凝土，可将其视为一个封闭的圆环，如图2所示，其支护刚度系数为：
Figure BDA0002464250820000102
其中Ec是混凝土衬砌的弹性模量(MPa)；vc是混凝土的泊松比；R是支护的最外侧半径(m)，通常可采用隧洞半径替代；tc是混凝土的厚度(m)。
对于钢拱架支护，可假设其紧密作用于围岩，支护刚度系数为：
Figure BDA0002464250820000103
其中Es1为钢(拱架)的弹性模量(MPa)；R是支护的最外侧半径(m)，通常可采用隧洞半径替代；tB为支护钢拱架圈的厚度(m)；B为槽钢或角钢的翼缘宽度(m)；As1为钢拱架的截面面积(m2)；Is为隧洞截面的惯性矩(m4)；S1为钢拱架沿隧洞轴线分布的间距(m)；θ为两个阻块节点夹角的一半(弧度)；EB为阻块材料的弹性模量(MPa)。
对于钢筋混凝土衬砌，其类似于喷射混凝土，只不过在混凝土中放置一定的钢筋形成钢筋混凝土圈对隧洞围岩进行支护，其支护刚度系数为：
Figure BDA0002464250820000111
其中Ec为混凝土的弹性模量(MPa)；d为超挖区填充混凝土的平均厚度(m)；Es1为钢(筋)的弹性模量(MPa)；R是支护的最外侧半径 (m)，通常可采用隧洞半径替代；As2为隧洞断面中钢筋的截面面积 (m2)；S2为钢筋沿隧洞轴线分布的间距(m)。
对于无粘结锚杆，可假定锚杆沿圆周是均匀分布的，其支护刚度系数为：
Figure BDA0002464250820000112
其中db为锚杆、锚索的直径(m)；l为锚杆、锚索的自由长度(m)； Q为锚杆端部和头部的变形荷载常量；Es2为锚杆锚索的弹性模量 (MPa)；Sc为锚杆沿圆周分布的间距(m)；Sl为锚杆沿隧洞轴线分布的间距(m)。
在上述基础上，相对变形ε可由
Figure BDA0002464250820000113
计算获得。当相对变形2.5％＜ε＜5％时，为I级围岩大变形，变形程度为中等；当相对变形5％＜ε＜10％时，为II级围岩大变形，变形程度为严重；当相对变形ε＞10％时，为III级围岩大变形，变形程度为极严重。
以半径r0＝2.965m的深埋圆形隧洞为例，围岩力学参数的中变形模量E为5GPa，泊松比为0.3，岩体内摩擦角为30°，黏聚力为0.3MPa，扩容梯度为1，侧压系数为1.25，预报了不同埋深条件下初始铅直地应力为4～30MPa时隧洞围岩变形量级和相对变形，给出了相应初始铅直地应力时对应的大变形等级，见图3所示。
尽管上面结合附图对本发明的优选实施例进行了描述，但是本发明并不局限于上述的具体实方式，上述的具体实施方式仅仅是示意性的，并不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可以作出很多形式的具体变换，这些均属于本发明的保护范围内。

Claims (4)

1.一种深埋长隧洞软质岩大变形超前量化预报方法，其特征在于：所述方法包括：
步骤一：中远距离超前地层信息预报技术：采用隧洞超前地质预报***，对掌子面前方设定距离的地层信息进行探测并获得三维空间数据，利用集面向软岩不同尺度变形质构造体的隧洞反射波正演模拟技术、波场分离和偏移成像的数据处理技术、地震数据保幅成像方法、基于深度学习的非线性智能反演方法于一体的信号解译技术，解译并预报隧洞掌子面前方地层的性质、位置及规模大小；
步骤二：近距离超前地质钻探预报技术：在步骤一工作的基础上，当隧洞开挖至软质岩30～60m范围时，在掌子面打超前钻孔；
步骤三：在步骤二钻孔部位进行原位地应力测试，获取该处三维初始应力分布状态、最大水平侧压力系数λ及初始最大主应力值σ0max参数；
步骤四：在步骤二钻孔部位进行岩石取样，开展室内物理力学试验，准确掌握掌子面前方软质岩石物理力学特性；
步骤五：在步骤二钻孔部位进行钻孔电视工作，获得软岩裂隙的埋深、倾向、倾角、宽度、裂隙面的粗糙度、充填物的性质，在步骤四岩石声波测试取样部位进行现场声波测试工作，获得岩体弹性纵波速度Vpm，结合步骤四室内试验获得的岩石弹性纵波速度Vpr，计算岩体完整性系数
Figure FDA0003898962360000011
步骤六：根据室内试验和现场试验，确定岩体基本质量分级，进一步确定掌子面前方软质岩体的物理力学参数：岩体基本质量分级，根据岩体基本质量的定性特征和岩体基本质量指标BQ两者相结合确定，根据岩体的基本质量级别确定岩体黏聚力c，岩体内摩擦角φ和岩体变形模量E；
步骤七：隧洞围岩大变形预测预报技术：计算自重应力场或者构造应力场下软质岩的强度应力比，利用围岩大变形判别方法，对围岩是否会发生挤压大变形进行判断和预测；
自重应力场(λ≤1)或者构造应力场(λ>1)下的岩石强度应力比值公式分别为
Figure FDA0003898962360000021
式中，Rc为岩石饱和单轴抗压强度；σ0max为初始最大主应力值；
围岩大变形判别方法为对相对变形ε量化计算获得，当相对变形2.5％＜ε＜5％时，为Ⅰ级围岩大变形，变形程度为中等；当相对变形5％＜ε＜10％时，为Ⅱ级围岩大变形，变形程度为严重；当相对变形ε＞10％时，为Ⅲ级围岩大变形，变形程度为极严重；
步骤八：在步骤七的基础上，若判别为发生围岩大变形，利用考虑支护效应修正的围岩变形和相对变形量化预报公式，通过最速下降法，快速量化预报围岩变形量级和相对变形，给出大变形等级，为软质岩支护优化设计和施工方案提供依据；
隧洞软质岩洞周变形量化预报，由公式
Figure FDA0003898962360000022
量化计算获得；式中，R′为塑性半径，计算公式为
Figure FDA0003898962360000023
r0为隧洞半径；σc为岩体单轴抗压强度，计算公式为
Figure FDA0003898962360000024
Figure FDA0003898962360000025
c为岩体的黏聚力，
Figure FDA0003898962360000026
为岩体的内摩擦角；其中
Figure FDA0003898962360000027
Figure FDA0003898962360000028
P为铅直向初始地应力；λ为最大水平侧压力系数；η为扩容梯度；反映围岩不同位置的θ角；E为岩体的变形模量；μ为泊松比；Pi为支护抗力；
考虑多种支护形式时，支护抗力Pi由公式
Figure FDA0003898962360000031
量化确定；其中Ec是混凝土的弹性模量；R是支护的最外侧半径；S1为钢拱架沿隧洞轴线分布的间距；S2为钢筋沿隧洞轴线分布的间距；Sc为锚杆沿隧洞圆周分布的间距；Sl为锚杆沿隧洞轴线分布的间距；As1为钢拱架的截面面积；As2为隧洞断面中钢筋的截面面积；Es1为钢的弹性模量；Es2为锚杆锚索的弹性模量；vc是混凝土的泊松比；tc是混凝土的厚度；tB为支护钢拱架圈的厚度；B为槽钢或角钢的翼缘宽度；Is为隧洞截面的惯性矩；θ为两个阻块节点夹角的一半弧度；EB为阻块材料的弹性模量；d为超挖区填充混凝土的平均厚度；其中db为锚杆、锚索的直径；l为锚杆、锚索的自由长度；Q为锚杆端部和头部的变形荷载常量。
2.根据权利要求1所述的一种深埋长隧洞软质岩大变形超前量化预报方法，其特征在于：所述步骤三中原位地应力测试，测试仪器采用三向应变计，测试方法采用应力解除法，并向孔内添加水泥砂浆作为耦合介质。
3.根据权利要求1所述的一种深埋长隧洞软质岩大变形超前量化预报方法，其特征在于：所述步骤四中岩石室内试验包括：进行岩石饱和单轴压缩试验，获得岩石饱和单轴抗压强度Rc和泊松比μ；进行岩石饱和三轴压缩试验，获得扩容梯度η；进行室内岩石声波测试，获得岩石弹性纵波速度Vpr。
4.根据权利要求1所述的一种深埋长隧洞软质岩大变形超前量化预报方法，其特征在于：所述步骤六中岩体基本质量指标BQ的计算公式为BQ＝100+3Rc+250Kv，使用该式时，遵循以下限制条件：(1)当Rc＞90Kv+30时，以Rc＝90Kv+30和Kv带入计算BQ值；(2)当Kv＞0.04Rc+0.4时，以Kv＝0.04Rc+0.4和Rc代入计算BQ值；式中Rc为岩石单轴饱和抗压强度；Kv为岩体完整性系数。
CN202010328971.4A 2020-04-23 2020-04-23 一种深埋长隧洞软质岩大变形超前量化预报方法 Active CN111551427B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010328971.4A CN111551427B (zh) 2020-04-23 2020-04-23 一种深埋长隧洞软质岩大变形超前量化预报方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010328971.4A CN111551427B (zh) 2020-04-23 2020-04-23 一种深埋长隧洞软质岩大变形超前量化预报方法

Publications (2)

Publication Number Publication Date
CN111551427A CN111551427A (zh) 2020-08-18
CN111551427B true CN111551427B (zh) 2023-01-17

Family

ID=71998167

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010328971.4A Active CN111551427B (zh) 2020-04-23 2020-04-23 一种深埋长隧洞软质岩大变形超前量化预报方法

Country Status (1)

Country Link
CN (1) CN111551427B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113466944A (zh) * 2021-08-13 2021-10-01 中铁二院工程集团有限责任公司 一种隧道软岩变形段落查找能干层的地球物理方法
CN114088524B (zh) * 2021-11-18 2023-04-28 中国电建集团成都勘测设计研究院有限公司 高围压层状岩体地下洞室纵轴线布置方位的设计方法
CN114547731B (zh) * 2022-01-20 2022-11-01 中国电建集团成都勘测设计研究院有限公司 含特定结构面的洞室边墙锚索自由长度确定方法和应用
CN114961776B (zh) * 2022-03-04 2023-04-07 西南交通大学 一种高原铁路高地应力软岩隧道超前支护设计方法
CN116861704B (zh) * 2023-09-04 2023-11-14 北京交通大学 一种高地应力软岩隧道大变形等级动态快速预测方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004279272A (ja) * 2003-03-17 2004-10-07 Tokyo Electric Power Co Inc:The 岩盤物性評価方法およびシステム、プログラムならびに記録媒体
CN102736124A (zh) * 2012-06-14 2012-10-17 北京市市政工程研究院 基于综合参数的隧道开挖围岩动态细化分级方法
CN105350972A (zh) * 2015-10-26 2016-02-24 中铁二十局集团有限公司 一种高地应力软弱围岩隧道开挖施工方法
CN108871262A (zh) * 2018-03-23 2018-11-23 长江水利委员会长江科学院 大埋深洞室挤压型围岩大变形判别方法
CN109506614A (zh) * 2018-10-09 2019-03-22 中铁二院工程集团有限责任公司 一种层状围岩大变形的判定方法
CN110375704A (zh) * 2019-07-26 2019-10-25 华中科技大学 一种深埋圆隧开挖围岩-衬砌受力与变形的预测方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004279272A (ja) * 2003-03-17 2004-10-07 Tokyo Electric Power Co Inc:The 岩盤物性評価方法およびシステム、プログラムならびに記録媒体
CN102736124A (zh) * 2012-06-14 2012-10-17 北京市市政工程研究院 基于综合参数的隧道开挖围岩动态细化分级方法
CN105350972A (zh) * 2015-10-26 2016-02-24 中铁二十局集团有限公司 一种高地应力软弱围岩隧道开挖施工方法
CN108871262A (zh) * 2018-03-23 2018-11-23 长江水利委员会长江科学院 大埋深洞室挤压型围岩大变形判别方法
CN109506614A (zh) * 2018-10-09 2019-03-22 中铁二院工程集团有限责任公司 一种层状围岩大变形的判定方法
CN110375704A (zh) * 2019-07-26 2019-10-25 华中科技大学 一种深埋圆隧开挖围岩-衬砌受力与变形的预测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Analysis of stress and deformation site-monitoring in fault zone of Wushaoling Tunnel under complex geological conditions;Lei Jun等;《Rock and Soil mechanics》;20081231;1367-1368 *
六盘山深埋隧洞工程围岩变形特征及定性预测研究;李国宝;《地下水》;20180525(第03期);135-136 *
复杂地质条件下跨流域调水超长深埋隧洞建设需研究的关键技术问题;钮新强等;《隧道建设(中英文)》;20190505(第04期);6-19 *

Also Published As

Publication number Publication date
CN111551427A (zh) 2020-08-18

Similar Documents

Publication Publication Date Title
CN111551427B (zh) 一种深埋长隧洞软质岩大变形超前量化预报方法
Hsiung et al. Evaluation of the wall deflections of a deep excavation in Central Jakarta using three-dimensional modeling
Fairhurst Stress estimation in rock: a brief history and review
CN107503361B (zh) 一种顺层岩质边坡软弱层应力监测及围岩加固方法
CN109537580A (zh) 一种基于bim技术的工程桩体设计与施工方法
CN204401664U (zh) 一种灌注桩桩底沉渣和桩端地基检测装置
Shen et al. Evaluation of the effect of rolling correction of double-o-tunnel shields via one-side loading
Panet et al. Convergence-confinement method for tunnel design
Chen et al. Dynamic response of segmental lining tunnel
Cheng et al. Analysis, design and construction of foundations
Li et al. Grouting reinforcement of shallow and small clearance tunnel
Poulos Ground movements–A hidden source of loading on deep foundations
Djizanne et al. Some aspects of the hydro-mechanical behaviour of Callovo-Oxfordian (COx) claystone around a gallery parallel to the principal horizontal minor stress
Lin et al. Reinforcement effects and safety monitoring index for high steep slopes: A case study in China
Gao et al. Influence of Benoto bored pile construction on nearby existing tunnel: A case study
Leca et al. Underground works in soils and soft rock tunneling
Deng et al. Rationality determination method and mechanical behavior of underground utility tunnels in a ground fissure environment
Wong et al. Extrusion analysis of a bolt-reinforced tunnel face with finite ground-bolt bond strength
Munir Development of correlation between rock classification system and modulus of deformation
Davies Predicting axially and laterally loaded pile behaviour using in-situ testing methods
Wang et al. Experimental investigation of crack initiation and propagation in the unreinforced masonry specimen subjected to vertical settlement
CN112307606A (zh) 浅埋公路隧道二次衬砌施做时间的确定方法
Yang et al. Large-scale model test for detecting pile defects using the parallel seismic method
Chen et al. Study of longitudinal deformation profiles in high-ground-stress mega-section tunnels based on the Hoek–Brown criterion
Reese et al. Developments in analysis of horizontal capacity of bored piles

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant