CN110159346B - 基于非均匀收敛模式的类矩形隧道地表变形预测方法 - Google Patents
基于非均匀收敛模式的类矩形隧道地表变形预测方法 Download PDFInfo
- Publication number
- CN110159346B CN110159346B CN201910393858.1A CN201910393858A CN110159346B CN 110159346 B CN110159346 B CN 110159346B CN 201910393858 A CN201910393858 A CN 201910393858A CN 110159346 B CN110159346 B CN 110159346B
- Authority
- CN
- China
- Prior art keywords
- tunnel
- section
- quasi
- surface deformation
- rectangular
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000012544 monitoring process Methods 0.000 claims abstract description 7
- 239000002689 soil Substances 0.000 claims abstract description 7
- 238000006073 displacement reaction Methods 0.000 claims description 21
- 238000009412 basement excavation Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000009933 burial Methods 0.000 claims description 2
- 238000010276 construction Methods 0.000 abstract description 9
- 238000010586 diagram Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
Landscapes
- Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Excavating Of Shafts Or Tunnels (AREA)
Abstract
本发明提供一种基于非均匀收敛模式的类矩形隧道地表变形预测方法，涉及隧道施工技术领域。本发明步骤如下：步骤1：将实际的隧道断面类比为类矩形隧道断面，在考虑隧道底部土体隆起的条件下假定隧道断面产生非均匀收敛；步骤2：建立基于非均匀收敛模式的类矩形隧道地表变形的随机介质预测模型；步骤3：依据隧道某断面的地表变形监测值，反分析步骤2公式中的三个未知计算参量；步骤4：将现场实际参数代入基于非均匀收敛模式的类矩形隧道地表变形的随机介质预测模型得到地表沉降值；本方法考虑到在实际工程中的多种因素影响，对地表变形的预测精度较高。
Description
技术领域
本发明涉及隧道施工技术领域，尤其涉及一种基于非均匀收敛模式的类矩形隧道地表变形预测方法。
背景技术
目前，针对于地铁隧道施工引起地表变形预测的经验公式法常常采用随机介质法。基于随机介质理论进行隧道地表变形预测时，通常假定断面收敛模式为均匀收敛。然而在实际工程中，由于受地层的初始应力、土体的软硬不均以及施工工艺等多种因素影响，其收敛形式并非是均匀的。因此原本的随机介质理论方法不能较准确的预测隧道施工引起的地表变形。
发明内容
本发明要解决的技术问题是针对上述现有技术的不足，提供一种基于非均匀收敛模式的类矩形隧道地表变形预测方法，本方法考虑到在实际工程中的多种因素影响，对地表变形的预测精度较高。
为解决上述技术问题，本发明所采取的技术方案是：
本发明提供一种基于非均匀收敛模式的类矩形隧道地表变形预测方法，包括如下步骤：
步骤1：将实际的隧道断面类比为类矩形隧道断面，在考虑隧道底部土体隆起的条件下假定隧道断面产生非均匀收敛，收敛位移大于0mm；
步骤2：建立基于非均匀收敛模式的类矩形隧道地表变形的随机介质预测模型，依据随机介质理论，非均匀收敛模式下的类矩形隧道开挖引起的地表沉降值W(X)为：
式中，二重积分的上下限分别为：
式中：X为隧道断面待求沉降点的横坐标值，β为地层主要影响角；π为圆周率；ξ为随机某单元x轴坐标值；η为随机某单元z轴坐标值；H为隧道中心点埋深；A为隧道断面半圆半径；B为隧道断面中心点至半圆圆心的水平距离；ΔA为隧道顶部收敛位移；ΔB为隧道底部隆起位移；
步骤3：依据隧道某断面的地表变形监测值，反分析步骤2公式中的三个未知计算参量：地层主要影响角β、隧道顶部收敛位移ΔA以及隧道底部隆起位移ΔB；采用方向加速度法取目标函数中最小值的一组参数v＝{ΔA,ΔB,tanβ}为结果；
定义目标函数为:
式中n为地表沉降测点数，Wi 0为第i测点地表沉降实测值，Wi为第i测点地表沉降计算值；v为待反演的未知参数；
步骤4：依据步骤3得到地层影响角β、隧道顶部收敛位移ΔA、隧道底部隆起位移ΔB,同时根据采用数值分析的方法，将现场实际参数代入基于非均匀收敛模式的类矩形隧道地表变形的随机介质预测模型得到地表沉降值W(X)；所述现场实际参数包括隧道中心点埋深H、隧道断面半圆半径A、隧道断面中心点至半圆圆心的水平距离B。
所述的类矩形隧道断面为矩形与两个半圆形的组合图形，所述矩形的宽边为半圆形的直径。
采用上述技术方案所产生的有益效果在于：本发明提供的一种基于非均匀收敛模式的类矩形隧道地表变形预测方法，本发明以随机介质理论为基础，考虑到在实际工程中，受地层的初始应力、土体的软硬不均以及施工工艺等多种因素影响，类矩形隧道断面收敛形式并非是均匀的，提出的类矩形隧道非均匀收敛模式与断面实际收敛规律相符，并通过岩土工程反分析获得计算所需的参数，因此对地表变形的预测精度较高，对实际工程具有一定的指导作用。
附图说明
图1为本发明实施例提供的单元开挖空间坐标图；
图2为本发明实施例提供的类矩形断面非均匀收敛示意图；
图3为本发明具体实施方式的非均匀收敛模式下的类矩形隧道地表沉降曲线及实际监测数据对比图。
具体实施方式
下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
本实施例的方法如下所述。
本实施例采用某市轨道交通3号线试验段采用类矩形盾构施工，此断面的竖向地表沉降监测值，见表1。表1中：X为监测点距隧道中心的水平距离，S为地表沉降实测值。
表1横向地表沉降实测值
本发明提供一种基于非均匀收敛模式的类矩形隧道地表变形预测方法，包括如下步骤：
步骤1：将实际的隧道断面类比为类矩形隧道断面，在考虑隧道底部土体隆起的条件下假定隧道断面产生非均匀收敛，收敛位移不为0mm；
基于随机介质理论进行隧道地表变形预测时，通常假定断面收敛模式为均匀收敛。然而在实际工程中，由于受地层的初始应力、土体的软硬不均以及施工工艺等多种因素影响，其收敛形式并非是均匀的。因此原本的随机介质理论方法不能较准确的预测隧道施工引起的地表变形；
所述的类矩形隧道断面为矩形与两个半圆形的组合图形，所述矩形的宽边为半圆形的直径。
步骤2：建立基于非均匀收敛模式的类矩形隧道地表变形的随机介质预测模型，依据随机介质理论，将隧道整体开挖等效成无穷多个微小单元开挖对地表影响的总和。假设地表以下E′深度任一点(θ′，ν′)处的某无限小单元dθ′dν′dE′开挖并且完全塌落，单元开挖空间坐标如图1所示。则最终的单元下沉值为：
其中：β为地层主要影响角；π为圆周率；θ′为该单元x轴坐标值，ν′为该单元y轴坐标值；E′为该单元z轴坐标值。
依据叠加原理，结合步骤一的假定，类矩形隧道非均匀收敛模式如图2所示。
非均匀收敛模式下的类矩形隧道开挖引起的地表沉降值W(X)为：
式中，二重积分的上下限分别为：
式中：X为隧道断面待求沉降点的横坐标值，β为地层主要影响角；π为圆周率；ξ为随机某单元x轴坐标值；η为随机某单元z轴坐标值；H为隧道中心点埋深；A为隧道断面半圆半径；B为隧道断面中心点至半圆圆心的水平距离；ΔA为隧道顶部收敛位移；ΔB为隧道底部隆起位移；
本实施例中H＝9.7m，A＝3.37m，B＝2.55m；
步骤3：依据岩土工程反分析，即指利用工程实际的监测数据，推断出正分析所需参数的过程。基于随机介质理论进行地表沉降计算时，由于某些计算参数无法用简单的方法确定，通过反分析步骤2公式中的三个未知计算参量：地层主要影响角β、隧道顶部收敛位移ΔA以及隧道底部隆起位移ΔB；采用方向加速度法(Powell法)，取目标函数中最小值的一组参数v＝{ΔA,ΔB,tanβ}为结果；
定义目标函数为:
式中n为地表沉降测点数，Wi 0为第i测点地表沉降实测值，Wi为第i测点地表沉降计算值；v为待反演的未知参数；
反分析是目标函数的参数寻优问题，通过对未知参数的不断修正，使工程实测值与理论计算值的差异到达最小，从而得到最优参数。采用Powell法编程计算，搜索令目标函数取得最小值的一组参数为ΔA＝45.14mm,ΔB＝17.23mm,tanβ＝0.467；
步骤4：依据步骤3得到地层影响角β、隧道顶部收敛位移ΔA、隧道底部隆起位移ΔB,同时将现场实际参数代入基于非均匀收敛模式的类矩形隧道地表变形的随机介质预测模型得到地表沉降值W(X)；由于计算公式的被积函数不可积，需要采用数值分析的方法,用MATLAB软件进行计算。所述现场实际参数包括隧道中心点埋深H、隧道断面半圆半径A、隧道断面中心点至半圆圆心的水平距离B。
如图3所示，将参数代入步骤二的公式进行计算，得到非均匀收敛模式下的类矩形隧道地表沉降曲线及实际监测数据对比图。由此可见本发明的方法预测结果与地表实测曲线变化较吻合，尤其表现于隧道中轴线处最大沉降位置，表明本专利方法能够较好的预测类矩形隧道的地表沉降规律。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明权利要求所限定的范围。
Claims (1)
1.一种基于非均匀收敛模式的类矩形隧道地表变形预测方法，其特征在于：包括如下步骤：
步骤1：将实际的隧道断面类比为类矩形隧道断面，在考虑隧道底部土体隆起的条件下假定隧道断面产生非均匀收敛，收敛位移大于0mm；所述类矩形隧道断面为矩形与两个半圆形的组合图形，所述矩形的宽边为半圆形的直径；
步骤2：建立基于非均匀收敛模式的类矩形隧道地表变形的随机介质预测模型，依据随机介质理论，非均匀收敛模式下的类矩形隧道开挖引起的地表沉降值W(X)为：
式中，二重积分的上下限分别为：
a＝H-A；b＝H+A；
e＝H-A+△A；
f＝H+A-△B；
式中：X为隧道断面待求沉降点的横坐标值，β为地层主要影响角；π为圆周率；ξ为随机某单元x轴坐标值；η为随机某单元z轴坐标值；H为隧道中心点埋深；A为隧道断面半圆半径；B为隧道断面中心点至半圆圆心的水平距离；ΔA为隧道顶部收敛位移；ΔB为隧道底部隆起位移；
步骤3：依据隧道某断面的地表变形监测值，反分析步骤2公式中的三个未知计算参量：地层主要影响角β、隧道顶部收敛位移ΔA以及隧道底部隆起位移ΔB；采用方向加速度法取目标函数中最小值的一组参数v＝{ΔA,ΔB,tanβ}为结果；
定义目标函数为:
式中n为地表沉降测点数，Wi 0为第i测点地表沉降实测值，Wi为第i测点地表沉降计算值；v为待反演的未知参数；
步骤4：依据步骤3得到地层影响角β、隧道顶部收敛位移ΔA、隧道底部隆起位移ΔB,同时采用数值分析的方法，将现场实际参数代入基于非均匀收敛模式的类矩形隧道地表变形的随机介质预测模型得到地表沉降值W(X)；所述现场实际参数包括隧道中心点埋深H、隧道断面半圆半径A、隧道断面中心点至半圆圆心的水平距离B。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910393858.1A CN110159346B (zh) | 2019-05-13 | 2019-05-13 | 基于非均匀收敛模式的类矩形隧道地表变形预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910393858.1A CN110159346B (zh) | 2019-05-13 | 2019-05-13 | 基于非均匀收敛模式的类矩形隧道地表变形预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110159346A CN110159346A (zh) | 2019-08-23 |
| CN110159346B true CN110159346B (zh) | 2020-07-31 |
Family
ID=67634253
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910393858.1A Expired - Fee Related CN110159346B (zh) | 2019-05-13 | 2019-05-13 | 基于非均匀收敛模式的类矩形隧道地表变形预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110159346B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110765630B (zh) * | 2019-10-31 | 2022-05-27 | 莆田学院 | 一种利用地表位移预测隧道收敛位移的方法 |
| CN113190902B (zh) * | 2021-04-30 | 2023-10-17 | 中铁十一局集团有限公司 | 一种隧道施工引起地表位移的预测方法及系统 |
| CN114722578B (zh) * | 2022-03-17 | 2024-06-07 | 中铁第一勘察设计院集团有限公司 | 一种隧道地表沉降计算方法 |
| CN114961751B (zh) * | 2022-05-17 | 2023-03-31 | 浙江大学 | 一种土岩复合地层中盾构掘进引起的土体位移预测方法 |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101457533A (zh) * | 2007-12-13 | 2009-06-17 | 上海隧道工程股份有限公司 | 基坑开挖实时检测方法 |
| CN102967481A (zh) * | 2012-12-13 | 2013-03-13 | 青岛理工大学 | 一种地面建筑与地下结构作用关系测定方法 |
| CN103310039A (zh) * | 2013-05-22 | 2013-09-18 | 青岛理工大学 | 一种地下施工对地面建筑结构影响的检测方法 |
| CN104123454A (zh) * | 2014-07-18 | 2014-10-29 | 青岛理工大学 | 单节理硬岩地下开挖地表沉降预测方法 |
| JP2015119095A (ja) * | 2013-12-19 | 2015-06-25 | 三菱電機株式会社 | 静止誘導機器 |
| CN104766129A (zh) * | 2014-12-31 | 2015-07-08 | 华中科技大学 | 一种基于时空信息融合的地铁盾构施工地表变形的预警方法 |
| CN104897527A (zh) * | 2015-06-11 | 2015-09-09 | 同济大学 | 可视化类矩形盾构同步注浆浆液扩散模式研究平台及应用 |
| CN106223967A (zh) * | 2016-09-20 | 2016-12-14 | 中铁工程装备集团有限公司 | 一种大埋深矩形或类矩形隧道盾构施工装置及方法 |
| CN106934251A (zh) * | 2017-04-25 | 2017-07-07 | 浙江大学城市学院 | 一种类矩形盾构隧道施工中地表沉降计算方法 |
| CN207182783U (zh) * | 2017-06-09 | 2018-04-03 | 浙江大学城市学院 | 一种模拟类矩形盾构推进造成地层损失的装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5302763B2 (ja) * | 2009-05-11 | 2013-10-02 | キヤノン株式会社 | 画像処理装置、画像処理方法 |
-
2019
- 2019-05-13 CN CN201910393858.1A patent/CN110159346B/zh not_active Expired - Fee Related
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101457533A (zh) * | 2007-12-13 | 2009-06-17 | 上海隧道工程股份有限公司 | 基坑开挖实时检测方法 |
| CN102967481A (zh) * | 2012-12-13 | 2013-03-13 | 青岛理工大学 | 一种地面建筑与地下结构作用关系测定方法 |
| CN103310039A (zh) * | 2013-05-22 | 2013-09-18 | 青岛理工大学 | 一种地下施工对地面建筑结构影响的检测方法 |
| JP2015119095A (ja) * | 2013-12-19 | 2015-06-25 | 三菱電機株式会社 | 静止誘導機器 |
| CN104123454A (zh) * | 2014-07-18 | 2014-10-29 | 青岛理工大学 | 单节理硬岩地下开挖地表沉降预测方法 |
| CN104766129A (zh) * | 2014-12-31 | 2015-07-08 | 华中科技大学 | 一种基于时空信息融合的地铁盾构施工地表变形的预警方法 |
| CN104897527A (zh) * | 2015-06-11 | 2015-09-09 | 同济大学 | 可视化类矩形盾构同步注浆浆液扩散模式研究平台及应用 |
| CN106223967A (zh) * | 2016-09-20 | 2016-12-14 | 中铁工程装备集团有限公司 | 一种大埋深矩形或类矩形隧道盾构施工装置及方法 |
| CN106934251A (zh) * | 2017-04-25 | 2017-07-07 | 浙江大学城市学院 | 一种类矩形盾构隧道施工中地表沉降计算方法 |
| CN207182783U (zh) * | 2017-06-09 | 2018-04-03 | 浙江大学城市学院 | 一种模拟类矩形盾构推进造成地层损失的装置 |
Non-Patent Citations (2)
| Title |
|---|
| 《基于ADINA的隧道开挖变形分析》;张向东等;《北方交通》;20081128;第72-74页 * |
| 《类矩形盾构施工引起地下管线竖向位移计算方法研究》;魏纲等;《市政技术》;20180310(第2期);第86-90页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110159346A (zh) | 2019-08-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110159346B (zh) | 基于非均匀收敛模式的类矩形隧道地表变形预测方法 | |
| CN106644206B (zh) | 一种用于浅埋隧道的围岩压力计算方法 | |
| CN104019815B (zh) | 基于铲车载荷监测的gnss平地作业路径动态规划与导航方法 | |
| CN108830012A (zh) | 一种基于全要素挖掘的结构变形预测方法 | |
| CN110319990A (zh) | 基于倾角仪优化布置的桥梁动挠度监测方法 | |
| CN106772639A (zh) | 地下铁质管线埋深磁偶极子构造法优化反演 | |
| CN115358311B (zh) | 地表变形监测多源数据融合处理方法 | |
| CN108303512B (zh) | 一种原位测试土水特征曲线的方法 | |
| CN118010033B (zh) | 深水灌浆钻孔精准定位方法 | |
| CN114372314A (zh) | 预测降压降水引起地面沉降的方法 | |
| CN206772282U (zh) | 海上风电基础绝对沉降监测装置 | |
| CN117188425A (zh) | 一种确定强夯有效加固深度和范围的方法 | |
| CN117553738A (zh) | 一种隧道围岩变形监测及预测的方法 | |
| CN108612075A (zh) | 一种监测深基坑水平位移的方法 | |
| CN104462748B (zh) | 地震过程中覆盖粘土层的液化场地桥梁桩基荷载—位移关系估算新方法 | |
| CN109711063A (zh) | 一种适用于上软下硬地层的车站暗挖进尺调整方法与装置 | |
| CN117824618A (zh) | 一种海底矿区可通行性地图获取方法 | |
| CN110531441B (zh) | 一种利用冷泉气体渗漏计算海洋流场的方法及处理终端 | |
| CN105971037B (zh) | 一种既有工程桩桩底深度确定方法 | |
| CN105220694B (zh) | 一种工程抛石量估算方法 | |
| CN114740540B (zh) | 一种基于方向约束的洋中脊区磁异常图构建方法及系统 | |
| CN113761457B (zh) | 基于测量的重力异常数据提取局部重力异常的方法 | |
| JP2019152005A (ja) | 転圧回数予測システム、及び転圧方法 | |
| CN111504244A (zh) | 海底管道在位状态的检测方法及检测系统 | |
| CN115758696A (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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200731 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |