CN117574024A - 一种破碎岩体的隧道围岩径深介质场计算方法 - Google Patents

一种破碎岩体的隧道围岩径深介质场计算方法 Download PDF

Info

Publication number
CN117574024A
CN117574024A CN202311701180.1A CN202311701180A CN117574024A CN 117574024 A CN117574024 A CN 117574024A CN 202311701180 A CN202311701180 A CN 202311701180A CN 117574024 A CN117574024 A CN 117574024A
Authority
CN
China
Prior art keywords
rock
tunnel
mass unit
rock mass
calculating
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.)
Pending
Application number
CN202311701180.1A
Other languages
English (en)
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.)
China Railway Tunnel Group Co Ltd CRTG
CRTG Survey and Design Institute Co Ltd
Original Assignee
China Railway Tunnel Group Co Ltd CRTG
CRTG Survey and Design Institute Co Ltd
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 China Railway Tunnel Group Co Ltd CRTG, CRTG Survey and Design Institute Co Ltd filed Critical China Railway Tunnel Group Co Ltd CRTG
Priority to CN202311701180.1A priority Critical patent/CN117574024A/zh
Publication of CN117574024A publication Critical patent/CN117574024A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Excavating Of Shafts Or Tunnels (AREA)

Abstract

本发明公开了一种破碎岩体的隧道围岩径深介质场计算方法,包括如下:步骤一、测得隧道周向围岩内某点处岩体单元局部坐标系下,沿三个正交方向的岩石介质尺寸分别为:lx′r、ly′r、lz′r;计算岩体单元局部坐标系下的岩石介质特征长度比:步骤二、所述岩石介质特征长度比是否满足下式,如果满足,则进行步骤三。步骤三、计算岩体单元点处的介质场量值:步骤四、计算岩体单元局部坐标系下的等效岩石特征长度比;步骤五、计算隧道整体坐标系下围岩径深介质场Me(x)。该计算方法为隧道周边围岩一维径深介质场的计算方法,具有简便实用的优点。

Description

一种破碎岩体的隧道围岩径深介质场计算方法
技术领域
本发明属于隧道监测技术领域,具体涉及一种破碎岩体的隧道围岩径深介质场计算方法。
背景技术
隧道工程修建过程中,开挖隧道轮廓内岩体后,周边围岩因改变了原先的赋存状态而产生应力状态调整,从而产生净空内变形,甚至失稳;支护结构用于控制周边围岩变形,确保隧道安全稳定;隧道大变形或者塌方则是由于周边围岩失稳。因此,隧道周边围岩的物理力学参数及所处的环境(围岩强度、变形参数、破碎程度、地质构造、应力环境、地下水赋存等)是隧道工程重点研究和控制的对象。
隧道周边围岩岩体,具有宏观区域性或场属性,其范围一般为3~5倍隧道跨度,高地应力软岩隧道该范围更大。围岩岩体主要由岩石介质组成,充填气相或液相的空隙介质占比较少,因此隧道围岩介质场的意义重大。
关于隧道围岩介质场,目前还没有明确的定义方法,以往研究更多是表征岩体破碎状态,或者完整岩石内部损伤裂纹发展演化等;破碎岩体岩石介质场的理论方法或物理模型研究现状还是空白,更多停留在反映岩体破碎程度的经验方法层面,尚未上升到介质场理论层面,如完整性系数、节理密度、岩体最小特征长度概念等。隧道围岩介质场的计算方法也是空白,亟需对隧道围岩介质场态进行科学合理的定义,并提出简便实用的计算方法,为隧道工程建设提供技术储备。
发明内容
本发明的目的是提供一种破碎岩体的隧道围岩径深介质场计算方法,为隧道周边围岩一维径深介质场的计算方法,具有简便实用的优点。
本发明采用以下技术方案:一种破碎岩体的隧道围岩径深介质场计算方法,该计算方法包括如下:
步骤一、测得隧道周向围岩内某点处岩体单元局部坐标系下,沿三个正交方向的岩石介质尺寸分别为:lx′r、ly′r、lz′r
计算岩体单元局部坐标系下的岩石介质特征长度比:
其中:(x′,y′,z′)为该点处岩体单元的局部坐标系,原点O′在该岩体单元的中心;l为岩体单元的长度;
某点处指的从隧道洞周向围岩内部径深x处;
步骤二、所述岩石介质特征长度比是否满足下式:
如果不满足,则结束;如果满足,则进行步骤三;
步骤三、计算岩体单元点处的介质场量值:
其中:Vr′为岩石介质体积;V为岩体单元体积;
步骤四、计算岩体单元局部坐标系下的等效岩石特征长度比:
其中:mer′为该点处岩体单元内岩石介质体积的等效平均尺寸;
步骤五、计算隧道整体坐标系下围岩径深介质场Me(x):
进一步地,在步骤五中,计算隧道整体坐标系下围岩径深介质场前,建立隧道原点随动整体坐标系,具体为:过岩体单元的原点O′作隧道轮廓的法线,其与隧道轮廓交点为随动整体坐标系的原点O,为坐标轴x轴正向,/>为x的坐标值。
进一步地,该岩体单元为正方体,边长为l,l取值为0.2~2m。
进一步地,该岩体单元为岩-气二相组成。
本发明的有益效果是:1.明确定义了隧道围岩介质场的定义,对隧道围岩介质场的量化计算奠定了理论基础。2.隧道围岩径深介质场的计算方法具有明确的物理意义,采用径深坐标系,等效岩石介质长度比,隧道围岩介质场的分量个数减少,物理意义明确,同时避免了计算介质场时,坐标系旋转导致的复杂验算,工程实用性强。3.提出了科学的计算方法,避免了人工经验的不确定性。
附图说明
图1为岩体单元立方体的平面示意图;
图2为明隧道原点随动坐标系图;
其中:1.岩体单元;2.隧道轮廓;3.周边围岩;4.岩石介质。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明一种破碎岩体的隧道围岩径深介质场计算方法适用于地下水不发育或少量地下水,岩体介质可视为或基本视为岩-气二相组成。
该一种破碎岩体的隧道围岩径深介质场计算方法包括如下步骤:
步骤一、测得隧道周边围岩3内某点处岩体单元1局部坐标系下,沿三个正交方向的岩石介质4尺寸分别为:lx′r、ly′r、Iz′r;岩体单元1可采用波速测试,或者现场取芯波速测试,或电法测试。
岩体单元1为正方体,边长l,l取值为0.2~2m;岩体单元1为:围岩内选定的某点,选择该点微小邻域范围,即该点处的微元体。
其中(x′,y′,z′)为该点处岩体单元1的局部坐标系,原点O′在该岩体单元1的中心,三个坐标轴方向分别平行于岩体单元1的六个面。某点处指的从隧道洞向周边围岩内部径深x处。
岩体单元1的尺寸的确定依据以下原则:根据隧道尺寸跨度,如果隧道跨度不大,如在5m内,则岩体单元1设定边长为0.2米;如果隧道跨度大,如在15米左右,则岩体单元1设定边长为2米。另外,还与周边围岩3相关,如果周边围岩3受到的影响范围交大,则I取较大值;一般l等于岩石芯样的长度,如果超过2米或小于0.2米,则不能真实反应岩体的介质场,因此将l取值限定在0.2~2m。
如图1所示,计算岩体单元1局部坐标系下的岩石介质4特征长度比:
上式表明,在隧道周边围岩的某点(x′,y′,z′)处,岩体单元1三个方向的岩石介质4特征长度比,是三维空间向量。
在图1中,与岩体单元1中的岩石介质4相邻的区域为空隙介质,即岩体单元1视为岩-气二相组成。
步骤二、各向同性条件验算,所述岩石介质(4)特征长度比是否满足下式:
如果三个方向的岩石介质4特征长度比不满足同性条件验算,则各向异性明显,则结束。如果三个方向的岩石介质4特征长度比满足同性条件验算,则进入步骤三。
各向同性条件验算式子中的数值选择,为采用宏观各向同性模型计算得到的结果,与采用各向异性模型计算得到的结果相比,其误差在10%以内,满足工程允许的误差要求。
步骤三、计算岩体单元1该点处的介质场量值:
该点处岩体单元1范围内,岩石介质4体积Vr′与该岩体单元体积V的比值,即单位体积范围内岩石介质4的占比,反映该点处岩体的破碎程度。根据的物理意义,表征该点处岩体单元1在局部坐标系下的介质标量值,不具有场量属性。
破碎岩体满足各向同性条件时,在同一点处,任意旋转岩体单元1的局部坐标系,岩体的破碎程度和岩石介质4占比均不变,即不随坐标系旋转而改变。
步骤四、计算岩体单元1局部坐标系下的等效岩石特征长度比mer′:
该式将岩石介质4特征长度比(mx′r,my′r,mx′r)这一空间三维向量等效为一个与岩体单元1局部坐标系无关的等效特征长度比mer′,从原先的三分量减少为单一量,同时mer′不受岩体单元1局部坐标系旋转的影响,便于计算。mer′的等效原则是:该点处岩体单元的岩石介质体积不变。等效前提条件:破碎岩体具有宏观各向同性。mer′物理意义:该点处岩体单元1内岩石介质体积的等效平均尺寸。
要计算隧道周边围岩4内部的介质场,需建立隧道原点随动整体坐标系,如图2所示,方法是:
过岩体单元1的原点O′作隧道轮廓2的法线,其与隧道轮廓2交点为随动整体坐标系的原点O,为坐标轴x轴正向,/>为x的坐标值,为一维整体坐标系。该坐标系的原点O随岩体单元1在隧道周边围岩4的环向位置不同而变动,即随岩体单元1的位置而动。坐标值x表征从隧道洞周至岩体单元1的径深。
从岩体单元1局部坐标系的(m′xr,m′yr,m′xr)到隧道原点随动整体坐标系的mxr(x),存在两个困难:其一是分量个数不同,物理意义不一致;其二是局部坐标系(O′,x′,y′,z′)到整体坐标系(O,x)的旋转分量计算困难。为此,根据介质场物理意义,引入等效岩石介质4长度比,将mer′赋予整体坐标系的径深坐标x,即可得到隧道原点随动整体坐标系的mer(x)=mer′,且mer(x)=mxr(x)。
步骤五、计算隧道整体坐标系下围岩径深介质场Me(x):

Claims (4)

1.一种破碎岩体的隧道围岩径深介质场计算方法,其特征在于,该计算方法包括如下:
步骤一、测得隧道周边围岩(3)内某点处岩体单元(1)局部坐标系下,沿三个正交方向的岩石介质(4)尺寸分别为:lx′r、ly′r、lz′r
计算岩体单元(1)局部坐标系下的岩石介质(4)特征长度比:
其中:(x′,y′z′)为该点处岩体单元(1)的局部坐标系,原点O′在该岩体单元(1)的中心;l为岩体单元(1)的长度;
某点处指的从隧道洞向周边围岩内部径深x处;
步骤二、所述岩石介质(4)特征长度比是否满足下式:
如果不满足,则结束;如果满足,则进行步骤三。
上式子中,
步骤三、计算岩体单元(1)点处的介质场量值:
其中:Vr′为岩石介质(4)体积;V为岩体单元(1)体积;
步骤四、计算岩体单元(1)局部坐标系下的等效岩石特征长度比:
其中:mer′为该点处岩体单元(1)内岩石介质(4)体积的等效平均尺寸;
步骤五、计算隧道整体坐标系下围岩径深介质场Me(x):
2.如权利要求1所述的一种破碎岩体的隧道围岩径深介质场计算方法,其特征在于,在所述步骤五中,计算隧道整体坐标系下围岩径深介质场前,建立隧道原点随动整体坐标系,具体为:过岩体单元(1)的原点O′作隧道轮廓(2)的法线,其与隧道轮廓(2)交点为随动整体坐标系的原点O,为坐标轴x轴正向,/>为x的坐标值。
3.如权利要求2所述的一种破碎岩体的隧道围岩径深介质场计算方法,其特征在于,所述岩体单元(1)为正方体,边长为l,l取值为0.2~2m。
4.如权利要求3所述的一种破碎岩体的隧道围岩径深介质场计算方法,其特征在于,所述岩体单元(1)为岩-气二相组成。
CN202311701180.1A 2023-12-12 2023-12-12 一种破碎岩体的隧道围岩径深介质场计算方法 Pending CN117574024A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311701180.1A CN117574024A (zh) 2023-12-12 2023-12-12 一种破碎岩体的隧道围岩径深介质场计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311701180.1A CN117574024A (zh) 2023-12-12 2023-12-12 一种破碎岩体的隧道围岩径深介质场计算方法

Publications (1)

Publication Number Publication Date
CN117574024A true CN117574024A (zh) 2024-02-20

Family

ID=89884441

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311701180.1A Pending CN117574024A (zh) 2023-12-12 2023-12-12 一种破碎岩体的隧道围岩径深介质场计算方法

Country Status (1)

Country Link
CN (1) CN117574024A (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040220783A1 (en) * 2001-06-22 2004-11-04 Egmont Lammer Method for determining and evaluating removal methods and device for carrying out said method
CN111695285A (zh) * 2020-06-17 2020-09-22 大连海事大学 一种各向异性岩体应力-损伤-渗流耦合数值模拟方法
CN112648017A (zh) * 2020-12-28 2021-04-13 中国科学院武汉岩土力学研究所 一种用于模拟隧道严重偏压大变形的数值方法
CN113268801A (zh) * 2021-06-04 2021-08-17 中铁隧道局集团有限公司 一种基于隧道场解重构理论的隧道围岩支护方法及系统
CN114996825A (zh) * 2022-07-07 2022-09-02 河北工程技术学院 一种深埋隧道挤压型大变形地质力学模式构建方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040220783A1 (en) * 2001-06-22 2004-11-04 Egmont Lammer Method for determining and evaluating removal methods and device for carrying out said method
CN111695285A (zh) * 2020-06-17 2020-09-22 大连海事大学 一种各向异性岩体应力-损伤-渗流耦合数值模拟方法
CN112648017A (zh) * 2020-12-28 2021-04-13 中国科学院武汉岩土力学研究所 一种用于模拟隧道严重偏压大变形的数值方法
CN113268801A (zh) * 2021-06-04 2021-08-17 中铁隧道局集团有限公司 一种基于隧道场解重构理论的隧道围岩支护方法及系统
CN114996825A (zh) * 2022-07-07 2022-09-02 河北工程技术学院 一种深埋隧道挤压型大变形地质力学模式构建方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
郗鹏程等: "基于统计岩体力学的隧道围岩分级方法", 地下空间与工程学报, vol. 14, no. 01, 15 February 2018 (2018-02-15) *

Similar Documents

Publication Publication Date Title
CN113221393B (zh) 一种基于非结构有限元法的三维大地电磁各向异性反演方法
CN112012712B (zh) 一种嵌入式离散裂缝的注水生长缝数值模拟方法和装置
EP2601642B1 (en) System and method for summarizing data on an unstructured grid
CN113821953B (zh) 基于动态地应力和剩余油特征的暂堵重复压裂设计方法
US20170074770A1 (en) Method for characterizing the fracture network of a fractured reservoir and method for exploiting it
EP2639604A2 (en) Method and system for presenting stress domain
CN105550448A (zh) 基于钻井轨迹设计参数的预钻井三维井眼建模方法及装置
Yang et al. Flow simulation of complex fracture systems with unstructured grids using the fast marching method
CN108875211A (zh) 一种基于fetd与fdtd耦合的二维模型的探地雷达正演方法
CN105156081A (zh) 一种碳酸盐岩稠油油藏酸化模拟评价方法
CN111695228A (zh) 一种缝洞型碳酸盐岩油藏的多尺度裂缝建模方法
CN115324560A (zh) 利用地应力场模拟确定压裂诱发油气套管变形位置的方法
CN113296166A (zh) 裂缝模型的构建方法
CN110969699A (zh) 一种高精度三维地应力建模方法及其应用
CN113779811B (zh) 一种构造变形对现今地应力方位扰动的定量分析方法
US20160162614A1 (en) Advanced Downhole Waveform Interpretation
Abes et al. The Impact of Geometric Attributes of Fractures on Fluid Flow Characteristics of Reservoir: A Case Study in Alrar Field, Algeria
Lu et al. Numerical investigation of unpropped fracture closure process in shale based on 3D simulation of fracture surface
Zhan et al. Comprehensive characterization and clustering of orientation data: A case study from the Songta dam site, China
CN117574024A (zh) 一种破碎岩体的隧道围岩径深介质场计算方法
Marache et al. Characterization and reconstruction of a rock fracture surface by geostatistics
Pandurangan et al. Tiltmeter mapping of measured nonsymmetric hydraulic-fracture growth in a conglomerate/sandstone formation using the implicit level-set algorithm and the extended Kalman filter
CN106869915B (zh) 一种水平井井间隔夹层预测方法及装置
CN107945271B (zh) 基于地质块体追踪的三维压力场建模方法
CN105092107A (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