CN107315879B - 富水隧道掌子面突水破坏临界水压及临界厚度的确定方法 - Google Patents

富水隧道掌子面突水破坏临界水压及临界厚度的确定方法 Download PDF

Info

Publication number
CN107315879B
CN107315879B CN201710520127.XA CN201710520127A CN107315879B CN 107315879 B CN107315879 B CN 107315879B CN 201710520127 A CN201710520127 A CN 201710520127A CN 107315879 B CN107315879 B CN 107315879B
Authority
CN
China
Prior art keywords
water
formula
critical
tunnel face
rich
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
CN201710520127.XA
Other languages
English (en)
Other versions
CN107315879A (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.)
Hunan University of Science and Technology
Original Assignee
Hunan University of Science and Technology
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 Hunan University of Science and Technology filed Critical Hunan University of Science and Technology
Priority to CN201710520127.XA priority Critical patent/CN107315879B/zh
Publication of CN107315879A publication Critical patent/CN107315879A/zh
Application granted granted Critical
Publication of CN107315879B publication Critical patent/CN107315879B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Civil Engineering (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Architecture (AREA)
  • Excavating Of Shafts Or Tunnels (AREA)
  • Lining And Supports For Tunnels (AREA)

Abstract

本发明公开了一种富水隧道掌子面突水破坏临界水压及临界厚度的确定方法。根据掌子面前方富水的压力、隧道开挖断面尺寸、隔水层厚度、围岩条件等情况，分别计算富水压力所做功率、破裂面能量耗散功率，依据能量守恒原理，求得含有隔水层的函数表达式，基于高等数学中函数求极值方法，求得在既有富水压力下的隔水层临界厚度；改变富水压力的情况，可以得到不同富水程度下掌子面前方临界隔水层厚度，对应的富水压力即临界突水压力；改变掌子面前方围岩条件，可以得到不同围岩条件下掌子面前方突水临界水压与临界隔水层厚度。本发明为判断富水隧道掌子面是否突水提供了计算方法；据此可以评价施工是否安全，隔水层厚度是否满足要求等。

Description

富水隧道掌子面突水破坏临界水压及临界厚度的确定方法
技术领域
本发明属于土木工程中隧道施工技术领域，具体涉及一种富水隧道掌子面三维突水破坏临界水压及临界厚度的确定方法。
背景技术
涌突水灾害是影响隧道施工的最大障碍之一，突水的发生机理和过程较为复杂，不同地区、不同地质条件、不同特点的工程部位突水机理都各不相同。突水的关键是如何确定临界突水压力和隔水层的安全厚度，从而指导设计与施工。开展隧道开挖条件下涌突水研究，对隧道突水预测与安全施工具有理论意义和工程价值。目前，隧道开挖过程中涌突水机理及过程模拟研究取得了一定的进展，但在复杂断层等软弱地质条件下隧道涌突水演化机理研究还有待深入。
发明内容
本发明的目的在于提供一种富水隧道掌子面突水破坏临界水压及临界厚度的确定方法。本发明为判断富水隧道掌子面是否突水提供了计算方法；据此可以评价施工是否安全，隔水层厚度是否满足要求，是否要对前方富水进行泄水降压，或对掌子面进行注浆加固以提升临界突水压力等。
本发明的目的是通过如下的技术方案来实现的：该富水隧道掌子面突水破坏临界水压及临界厚度的确定方法，包括如下步骤：
(1)确定破坏体之间的几个关系如下：
Figure GDA0002413705670000011
Figure GDA0002413705670000021
Figure GDA0002413705670000022
Figure GDA0002413705670000023
式中：dA为破坏面微元面积；r(x)为破坏面到破坏体中心线x轴的垂直距离；a为破坏面上任意一点的切平面与掌子面前方破坏体的速度v的夹角；r与r(x)表达的是同一个意思；r'为r的一阶导数；以下的r'(x)与r'表达的是同一个意思；
(2)隧道掌子面前方富水压力所做功率为：
Figure GDA0002413705670000024
式中：Wp为富水压力所做功率；p为含水层的水压力；d为掌子面前方破裂的宽度；v为掌子面前方破坏体的速度；
(3)掌子面破坏面的能量耗散率为沿整个破坏面的积分：
Figure GDA0002413705670000025
式中：ED为掌子面破坏面的能量耗散率；A为整个破坏面的面积；WA为破坏面微元面积dA上的能量耗散率；
其中：
Figure GDA0002413705670000026
而：
Figure GDA0002413705670000027
所以：
Figure GDA0002413705670000028
式中：σ1为最大主应力；σ3为最小主应力；σt为岩石的抗拉强度；
结合几何关系，可以求得：
Figure GDA0002413705670000029
式中：H为隔水层厚度；
(4)由能量守恒原理，即外力功率与内能耗散功率相等，可以得到：
Wp＝ED；
(5)将步骤(2)公式与步骤(3)公式带入上式，可以得到掌子面前方含水层的水压力p为：
Figure GDA0002413705670000031
(6)为求掌子面临界突水的含水层的水压力p，令泛函数为：
Figure GDA0002413705670000032
式中：ψ[r(x),r'(x),x]称为泛函数；
(7)步骤(6)中的泛函数存在极值时，则由泛函数的变分原理，可得其对应的欧拉方程为：
Figure GDA0002413705670000033
求解可得：
Figure GDA0002413705670000034
式中：r(x)为破坏面到破坏体中心线x轴的垂直距离；ψ为ψ[r(x),r'(x),x]，即泛函数；c1、c2是对欧拉方程求解并进行积分而得的常系数；
(8)步骤(7)中的系数c1、c2由以下条件确定：
Figure GDA0002413705670000035
可得：
Figure GDA0002413705670000036
式中：D为隧道的开挖高度；
(9)将步骤(8)中的系数回代步骤(7)的公式中，即可求得破坏面的方程：
Figure GDA0002413705670000037
对上式转换则得：
Figure GDA0002413705670000041
令
Figure GDA0002413705670000042
则
Figure GDA0002413705670000043
根据上述方程即可绘制掌子面突水的三维破裂面形状；式中，y表示三维破裂面的y坐标值；z表示三维破裂面的z坐标值；即破裂线绕x轴旋转，r不变，但破裂面上的y与z的坐标在变化；
(10)将求得的r(x)代入步骤(5)公式中，即可求得含水层的水压力p：
Figure GDA0002413705670000044
上式中，d为未知量，通过微分求极值的条件可以得到：
Figure GDA0002413705670000045
将上述中微分求得的d，回代入p的表达式中，即为在已知隔水层厚度H下的临界突水的含水层的水压力；
(11)对步骤(10)中p的表达式进行变换可以求得：
Figure GDA0002413705670000046
上式中，d为未知量，通过微分求极值的条件可以得到：
Figure GDA0002413705670000047
将上述中微分求得的d，回代入H的表达式中，即为在已知含水层的水压力p下的临界隔水层厚度H。
本发明与现有研究方法相比的优点在于：突水破裂面形状函数不用人为假定，其根据能量守恒原理求得的关系式，并结合泛函数求极值的方法确定。本发明方法为判断富水隧道掌子面是否突水提供了计算方法；据此可以评价施工是否安全，隔水层厚度是否满足要求，是否要对前方富水进行泄水降压，或对掌子面进行注浆加固以提升临界突水压力等。本发明的方法，亦可以应用于采矿巷道、水工隧洞、地铁等地下建筑结构在富水情况下掌子面或边墙部位是否发生突水破坏的评价，以及临界水压与临界隔水层厚度的确定中。
附图说明
图1是本发明的计算方法示意图。
图2是本发明实施例的含水层的水压力与掌子面前方破裂的宽度的关系曲线图。
图3是本发明实施例的三维突水时的破裂面图。
图中，D为隧道开挖高度；d为掌子面前方破裂的宽度；v为掌子面前方破坏体的速度；a为破坏面上任意一点的切平面与速度v的夹角；τn为围岩破裂面处的剪应力，即抗剪强度；σn为围岩破裂面处的正应力，即法向应力；p为含水层的水压力，在岩溶地段，可以为溶洞/溶腔的水压力；H为隔水层厚度。
具体实施方式
下面结合附图和实施例对本发明作进一步详细的描述。
参见图1，本发明的富水隧道掌子面突水破坏临界水压及临界厚度的确定方法，步骤如下：
(1)确定破坏体之间的几个关系如下：
Figure GDA0002413705670000051
Figure GDA0002413705670000052
Figure GDA0002413705670000053
Figure GDA0002413705670000061
式中：dA为破坏面微元面积；r(x)为破坏面到破坏体中心线x轴的垂直距离；a为破坏面上任意一点的切平面与掌子面前方破坏体的速度v的夹角；r与r(x)表达的是同一个意思；r'为r的一阶导数；以下的r'(x)与r'表达的是同一个意思；
(2)隧道掌子面前方富水压力所做功率为：
Figure GDA0002413705670000062
式中：Wp为富水压力所做功率；p为含水层的水压力；d为掌子面前方破裂的宽度；v为掌子面前方破坏体的速度；
(3)掌子面破坏面的能量耗散率为沿整个破坏面的积分：
Figure GDA0002413705670000063
式中：ED为掌子面破坏面的能量耗散率；A为整个破坏面的面积；WA为破坏面微元面积dA上的能量耗散率；
其中：
Figure GDA0002413705670000064
而：
Figure GDA0002413705670000065
所以：
Figure GDA0002413705670000066
式中：σ1为最大主应力；σ3为最小主应力；σt为岩石的抗拉强度；
结合几何关系，可以求得：
Figure GDA0002413705670000067
式中：H为隔水层厚度；
(4)由能量守恒原理，即外力功率与内能耗散功率相等，可以得到：
Wp＝ED；
(5)将步骤(2)公式与步骤(3)公式带入上式，可以得到掌子面前方含水层的水压力p为：
Figure GDA0002413705670000071
(6)为求掌子面临界突水的含水层的水压力p，令泛函数为：
Figure GDA0002413705670000072
式中：ψ[r(x),r'(x),x]称为泛函数；
(7)步骤(6)中的泛函数存在极值时，则由泛函数的变分原理，可得其对应的欧拉方程为：
Figure GDA0002413705670000073
求解可得：
Figure GDA0002413705670000074
式中：r(x)为破坏面到破坏体中心线x轴的垂直距离；ψ为ψ[r(x),r'(x),x]，即泛函数；c1、c2是对欧拉方程求解并进行积分而得的常系数；
(8)步骤(7)中的系数c1、c2由以下条件确定：
Figure GDA0002413705670000075
可得：
Figure GDA0002413705670000076
式中：D为隧道的开挖高度；
(9)将步骤(8)中的系数回代步骤(7)的公式中，即可求得破坏面的方程：
Figure GDA0002413705670000077
对上式转换则得：
Figure GDA0002413705670000081
令
Figure GDA0002413705670000082
则
Figure GDA0002413705670000083
根据上述方程即可绘制掌子面突水的三维破裂面形状；式中，y表示三维破裂面的y坐标值；z表示三维破裂面的z坐标值；即破裂线绕x轴旋转，r不变，但破裂面上的y与z的坐标在变化；
(10)将求得的r(x)代入步骤(5)公式中，即可求得含水层的水压力p：
Figure GDA0002413705670000084
上式中，d为未知量，通过微分求极值的条件可以得到：
Figure GDA0002413705670000085
将上述中微分求得的d，回代入p的表达式中，即为在已知隔水层厚度H下的临界突水的含水层的水压力；
(11)对步骤(10)中p的表达式进行变换可以求得：
Figure GDA0002413705670000086
上式中，d为未知量，通过微分求极值的条件可以得到：
Figure GDA0002413705670000087
将上述中微分求得的d，回代入H的表达式中，即为在已知含水层的水压力p下的临界隔水层厚度H。
参见图2、图3，本发明的上述计算方法给出了一个算例。对于隧道开挖高度D＝5m，隧道抗拉强度为σt＝50kPa，隔水层厚度为H＝3m，由步骤(10)，可以绘制出含水层的水压力p与掌子面前方破裂的宽度d的关系图，如图2所示；从图2中可以看出在d＝3.2m，水压力p取最小，即临时突水压力，此时p＝8105kPa。把d＝3.2m带入到步骤(9)公式中，即可绘制三维突水时的破裂面图，如图3所示。

Claims (1)

1.一种富水隧道掌子面突水破坏临界水压及临界厚度的确定方法，其特征在于包括如下步骤：
(1)确定破坏体之间的几个关系如下：
Figure FDA0002413705660000011
Figure FDA0002413705660000012
Figure FDA0002413705660000013
Figure FDA0002413705660000014
式中：dA为破坏面微元面积；r(x)为破坏面到破坏体中心线x轴的垂直距离；a为破坏面上任意一点的切平面与掌子面前方破坏体的速度v的夹角；r与r(x)表达的是同一个意思；r'为r的一阶导数；以下的r'(x)与r'表达的是同一个意思；
(2)隧道掌子面前方富水压力所做功率为：
Figure FDA0002413705660000015
式中：Wp为富水压力所做功率；p为含水层的水压力；d为掌子面前方破裂的宽度；v为掌子面前方破坏体的速度；
(3)掌子面破坏面的能量耗散率为沿整个破坏面的积分：
Figure FDA0002413705660000016
式中：ED为掌子面破坏面的能量耗散率；A为整个破坏面的面积；WA为破坏面微元面积dA上的能量耗散率；
其中：
Figure FDA0002413705660000017
而：
Figure FDA0002413705660000021
所以：
Figure FDA0002413705660000022
式中：σ1为最大主应力；σ3为最小主应力；σt为岩石的抗拉强度；
结合几何关系，可以求得：
Figure FDA0002413705660000023
式中：H为隔水层厚度；
(4)由能量守恒原理，即外力功率与内能耗散功率相等，可以得到：
Wp＝ED；
(5)将步骤(2)公式与步骤(3)公式带入上式，可以得到掌子面前方含水层的水压力p为：
Figure FDA0002413705660000024
(6)为求掌子面临界突水的含水层的水压力p，令泛函数为：
Figure FDA0002413705660000025
式中：ψ[r(x),r'(x),x]称为泛函数；
(7)步骤(6)中的泛函数存在极值时，则由泛函数的变分原理，可得其对应的欧拉方程为：
Figure FDA0002413705660000026
求解可得：
Figure FDA0002413705660000027
式中：ψ为ψ[r(x),r'(x),x]，即泛函数；c1、c2是对欧拉方程求解并进行积分而得的常系数；
(8)步骤(7)中的系数c1、c2由以下条件确定：
Figure FDA0002413705660000031
可得：
Figure FDA0002413705660000032
式中：D为隧道的开挖高度；
(9)将步骤(8)中的系数回代步骤(7)的公式中，即可求得破坏面的方程：
Figure FDA0002413705660000033
对上式转换则得：
Figure FDA0002413705660000034
令
Figure FDA0002413705660000035
则
Figure FDA0002413705660000036
根据上述方程即可绘制掌子面突水的三维破裂面形状；式中，y表示三维破裂面的y坐标值；z表示三维破裂面的z坐标值；即破裂线绕x轴旋转，r不变，但破裂面上的y与z的坐标在变化；
(10)将求得的r(x)代入步骤(5)公式中，即可求得含水层的水压力p：
Figure FDA0002413705660000037
上式中，d为未知量，通过微分求极值的条件可以得到：
Figure FDA0002413705660000038
将上式 中微分求得的d，回代入p的表达式中，即为在已知隔水层厚度H下的临界突水的含水层的水压力；
(11)对步骤(10)中p的表达式进行变换可以求得：
Figure FDA0002413705660000039
上式中，d为未知量，通过微分求极值的条件可以得到：
Figure FDA0002413705660000041
将上式 中微分求得的d，回代入H的表达式中，即为在已知含水层的水压力p下的临界隔水层厚度H。
CN201710520127.XA 2017-06-30 2017-06-30 富水隧道掌子面突水破坏临界水压及临界厚度的确定方法 Active CN107315879B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710520127.XA CN107315879B (zh) 2017-06-30 2017-06-30 富水隧道掌子面突水破坏临界水压及临界厚度的确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710520127.XA CN107315879B (zh) 2017-06-30 2017-06-30 富水隧道掌子面突水破坏临界水压及临界厚度的确定方法

Publications (2)

Publication Number Publication Date
CN107315879A CN107315879A (zh) 2017-11-03
CN107315879B true CN107315879B (zh) 2020-05-26

Family

ID=60181256

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710520127.XA Active CN107315879B (zh) 2017-06-30 2017-06-30 富水隧道掌子面突水破坏临界水压及临界厚度的确定方法

Country Status (1)

Country Link
CN (1) CN107315879B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110442948B (zh) * 2019-07-30 2023-04-07 湖南科技大学 上软层厚度与倾角对隧道掌子面稳定性的确定方法
CN111274639B (zh) * 2020-01-20 2022-08-26 西南交通大学 掌子面裂隙水突出破坏的隔水层最小安全厚度的判识方法
CN112131628B (zh) * 2020-07-29 2024-03-26 北京工业大学 一种隧道穿越富水断层防突水最小安全厚度计算方法
CN113030400A (zh) * 2021-03-31 2021-06-25 郑州工程技术学院 一种富水断层隧道掌子面突水演化模拟试验装置及方法
CN114183199A (zh) * 2021-11-11 2022-03-15 中国建筑第七工程局有限公司 一种基于掌子面位移突变理论的安全岩柱厚度确定方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105974056A (zh) * 2016-04-30 2016-09-28 山东大学 隧道突水灾害前兆信息监测模型试验系统及试验方法
CN106570287A (zh) * 2016-11-10 2017-04-19 中国人民解放军理工大学 一种基于三维离散裂隙网络预测隧道突涌水量的方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105974056A (zh) * 2016-04-30 2016-09-28 山东大学 隧道突水灾害前兆信息监测模型试验系统及试验方法
CN106570287A (zh) * 2016-11-10 2017-04-19 中国人民解放军理工大学 一种基于三维离散裂隙网络预测隧道突涌水量的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Study on safe thickness of comparatively intact rock ahead of karst tunnel face;Guo J Q, et al.;《Applied Mechanics & Materials》;20111231;2456-2459 *
岩溶隧道掌子面突水机制及岩墙安全厚度研究;郭佳奇，等;《铁道学报》;20120331;第34卷(第3期);109-115 *

Also Published As

Publication number Publication date
CN107315879A (zh) 2017-11-03

Similar Documents

Publication Publication Date Title
CN107315879B (zh) 富水隧道掌子面突水破坏临界水压及临界厚度的确定方法
CN102306225B (zh) 多线叠交隧道施工过程及对隧道变形影响数值模拟方法
CN104408323A (zh) 一种基于多源信息融合的采场顶板离层水害超前预报方法
CN106593524A (zh) 一种固体充填开采地表沉陷动态预计方法
CN109505628A (zh) 一种富水断层带隧道施工方法
CN103510948A (zh) 一种适用于脆性岩石气体钻井井壁稳定性预测的实验方法
Lei et al. A structural calculation model of shield tunnel segment: heterogeneous equivalent beam model
CN110926288B (zh) 一种考虑空洞效应的浅埋隧道爆破掏槽药量的优化方法
CN107908901A (zh) 一种基坑施工中承压水降水对周边环境影响的计算方法
CN107315880B (zh) 水平地震力作用下隧道直边墙三维破坏形状的定位方法
CN103150443A (zh) 一种微倾斜基坑围护结构的安全系数解算方法
CN111274639B (zh) 掌子面裂隙水突出破坏的隔水层最小安全厚度的判识方法
CN102002962A (zh) 堵水帷幕渗透性反演分析施工确定方法
CN110442948B (zh) 上软层厚度与倾角对隧道掌子面稳定性的确定方法
CN107330272B (zh) 隧道拱顶三维突水破坏临界水压与临界厚度的确定方法
CN103821533B (zh) 端墙式隧道洞门整体抗倾覆设计方法
CN103835281B (zh) 一种盾构穿越软弱地基堤防的加固方案快速确定的方法
CN106192991A (zh) 一种用于低净空条件下的地下连续墙施工方法
CN112664259A (zh) 高突煤层三位一体综合防突的方法
CN104881536B (zh) 一种岩体结构面吻合度系数测量方法
Zhang et al. Excess pore water pressure and ground consolidation settlement caused by grouting of shield tunnelling
Lv et al. Influence of water inrush from excavation surface on the stress and deformation of tunnel-forming structure at the launching-arrival stage of subway shield
CN117313583A (zh) 断层上盘隧道开挖防突厚度与临界水压的确定方法
Zhang et al. Study on the evaluation of surrounding rock stability and limit displacement of high-pressure stratum tunnel
CN117349936A (zh) 水下隧道近断层下盘开挖进尺安全性计算方法

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