CN116956420A - 一种基于三维地质模型的隧道结构一体化分析方法 - Google Patents
一种基于三维地质模型的隧道结构一体化分析方法 Download PDFInfo
- Publication number
- CN116956420A CN116956420A CN202310923252.0A CN202310923252A CN116956420A CN 116956420 A CN116956420 A CN 116956420A CN 202310923252 A CN202310923252 A CN 202310923252A CN 116956420 A CN116956420 A CN 116956420A
- Authority
- CN
- China
- Prior art keywords
- tunnel
- stratum
- section
- geological
- depth
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000012351 Integrated analysis Methods 0.000 title claims abstract description 13
- 238000004458 analytical method Methods 0.000 claims abstract description 19
- 238000004364 calculation method Methods 0.000 claims abstract description 19
- 230000008859 change Effects 0.000 claims abstract description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 14
- 238000009933 burial Methods 0.000 claims description 9
- 238000013461 design Methods 0.000 claims description 7
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 238000005755 formation reaction Methods 0.000 claims description 6
- 238000005553 drilling Methods 0.000 claims description 5
- 238000011835 investigation Methods 0.000 claims description 5
- 239000002689 soil Substances 0.000 claims description 4
- 230000005484 gravity Effects 0.000 claims description 3
- 239000011435 rock Substances 0.000 claims description 3
- 238000009412 basement excavation Methods 0.000 claims description 2
- 239000000463 material Substances 0.000 claims description 2
- 239000003673 groundwater Substances 0.000 abstract description 6
- 238000013499 data model Methods 0.000 abstract 1
- 238000003696 structure analysis method Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
Abstract
本发明公开了一种基于三维地质模型的隧道结构一体化分析方法,属于隧道工程技术领域,主要解决目前隧道数值计算分析计算三维数据模型大,难以准确获取隧道断面计算分析结果的问题。本发明通过建立三维地质、隧道模型,根据隧道里程提取对应的地质剖面,并获取地质信息。分析判断当前断面地质信息,以隧道埋深、地下水位、地层变化为判断条件,确定隧道危险断面,提取危险断面对应的地质剖面信息,并建立隧道计算分析模型,采用地层结构法进行二维有限元分析得到隧道断面内力、变形等结果。
Description
技术领域
本发明涉及隧道结构的数值分析方法,尤其是公开一种基于三维地质模型的隧道结构一体化分析方法,属于隧道工程技术领域。
背景技术
对交通的需求促使了隧道的修建和发展,使得穿山越岭成为可能。隧道结构设计计算是工程重点问题,隧道结构与地层之间的相互作用方式是计算的关键所在。随着隧道工程快速发展,隧道结构计算理论逐渐形成体系,包括荷载结构法和地层结构法两种基本的计算模型。在荷载结构法将地层对隧道结构的约束视为施加在隧道结构上不同分布形式的弹性抗力,建立隧道结构计算模型,并分析结构变形、内力。地层结构法将地层和隧道隧道结构结构视为受力整体,利用连续介质力学方法进行整体结构计算分析。荷载结构法计算的荷载往往偏大,经济性较差,同时也不能准确反映地层和隧道结构之间真实的相互作用关系及共同作用。
隧道建设、安全运营离不开隧道结构计算分析,地层结构法在反映真实地层和结构作用上具有很大优势,然而建立整体三维地质、三维隧道模型,存在模型复杂、计算效率低等问题。如何通过地质模型数据判断隧道结构危险断面,进行相应断面的隧道结构分析方法的研究具有重要意义。
发明内容
本发明的目的在于,针对隧道结构计算分析提供一种基于三维地质模型的隧道结构一体化分析方法,主要解决现有分析方法存在模型复杂、计算效率低等技术问题;采用地质剖切的方式获取对应里程的地层信息,根据隧道埋深、地下水位和隧道断面地层变化进行危险断面判断,建立危险断面的地层结构法数值分析模型,创建有限元分析模型计算隧道结构内力,实现基于地质数据判断危险断面判断、并进行隧道结构分析。
本发明解决其技术问题所采用的技术方案是:一种基于三维地质模型的隧道结构一体化分析方法,包括以下步骤:
S1:基于地质勘察、工程设计资料建立三维地质、隧道模型;
S2:将隧道设计信息(埋深、里程)与三维地质模型关联,通过设定的判断条件自动计算分析,获取危险断面里程;
S3:根据三维地质数据,自动获取危险断面的边界计算条件,分析计算当前隧道断面的位移、应力等。
所述S1具体步骤如下:
(1)统一地质、隧道模型坐标系,根据设计资料确定建模范围,地质模型宽度一般为8D(D为隧道直径),轴线左右各4D,深度为隧道底部以下2D;
(2)根据地层、钻孔、地下水位等地质勘察数据建立三维地质模型,地层数据包括各地层重度γ、泊松比υ、水平基床系数K、弹性模量E;钻孔数据包括孔口标高、平面位置坐标、平面位置坐标、各地层标高;
(3)根据隧道轴线、隧道断面尺寸信息在三维地质模型中创建隧道三维模型。
所述S2包括根据隧道埋深(隧道拱顶至地面的距离)判断危险断面、根据地下水位判断危险断面、根据地层变化判断危险断面;提取剖面的方法为沿隧道轴线,根据隧道里程每Bm提取地质剖面Geoi,盾构法隧道中B取管片环宽,矿山法隧道B取隧道开挖进尺。
所述根据隧道埋深判断危险断面的步骤如下:
(1)遍历Geoi计算隧道埋深Depthi;
(2)根据Depthi计算求得隧道最大埋深Depthmax,记录对应的隧道里程MDmax、地质剖面GeoDmax;
(3)根据Depthi计算求得隧道最小埋深Depthmin,记录对应的隧道里程MDmin、地质剖面GeoDmin。
所述根据地下水位判断危险断面的步骤如下:
(1)遍历Geoi计算地下水位高度WaterHj;
(2)根据WaterHj计算求得最高地下水位WaterHmax,记录对应的隧道里程MWmax、地质剖面GeoWmax;
(3)根据WaterHj计算求得最低地下水位WaterHmin,记录对应的隧道里程MWmin、地质剖面GeoWmin。
所述根据地层变化判断危险断面的步骤如下:
(1)从提取的Geoi中每x个断面进行一次地层变化判断,一般x取1~10之间的整数;
(2)第k个剖面对应的隧道里程为Milek=Miles+k×x×B,k=x,2x,3x…;
(3)逐个分析S1提取剖面地层参数,采用以下方法进行危险断面判断:
①当第k-1个剖面的地层中有地层a,第k个剖面的地层中无地层a,则第k个剖面对应的隧道断面为危险断面;
②当第k-1个剖面的地层中无地层b,第k个剖面的地层中出现了地层b,则第k个剖面对应的隧道断面为危险断面;
③第k个剖面Depthk~(Depthk+D)的深度范围存在软土、杂填土或软硬岩交替地层时,则k个剖面对应的隧道断面为危险断面;
④第k个剖面对应的隧道顶部埋深Depthk,在Depthk~(Depthk+D)的深度范围内地层数量超过nmax时(一般取nmax=3),则第j个剖面对应的隧道断面为危险断面。
所述S3具体步骤如下:
(1)以危险断面地质剖面GeoDmax、GeoDmin、GeoWmax、GeoWmin、Geok为基础;
(2)根据地质剖面图提取各地层边界线,并将地质剖面数据转换为边界荷载;
(3)将地层面作为平面单元,盾构管片作为梁单元,根据地质模型获取各地层γ、υ、K、E,赋值给相应的平面单元材料参数;
(4)施加荷载:考虑自重,采用均布荷载施加地面超载,地层底部约束刚体位移,进行网格划分得到有限元分析模型;
(5)采用地层结构法做二维有限元分析。
附图说明
图1本发明提供的基于地质数据的隧道结构分析方法流程图。
图2实施例中的地质模型和隧道示意图。
图3实施例中的地质剖面示意图。
具体实施方式
下面结合图1~3和具体实施例对本发明作进一步说明。
如图1所示为本发明中所述的一种基于地质数据的隧道结构分析方法的计算流程。以图2所示的某盾构工法隧道及三维地质模型为例阐述本发明所述的分析方法的计算过程。图中line_geo为地层线,隧道起始里程为Miles=K1+100,管片环宽B=2m,外径D=6.2m,内径d=5.5m。
根据钻孔、地层、隧道信息建立三维地质、隧道模型。地质模型上表面为勘察得到的实际地表面,模型宽度设置为49.6m,深度为隧道底部以下12.4m;根据隧道轴线、隧道尺寸、隧道环宽在三维地质模型中创建隧道模型。本实施例中地质模型各地层面如图2所示,隧道轴线终止里程为Milet=K1+500,地质模型的长度为Milet-Miles+2B,隧道起始、终止位置分别延伸一个环宽。
通过隧道埋深、地下水位、地层变化3种方法自动获取隧道全长中的危险断面,从起始里程开始,每2m提取地质模型剖面Geoi,共计提取100个地质剖面。
根据隧道埋深获取危险断面:遍历100个地质剖面,得到埋深最大的隧道里程为MDmax、埋深最小的隧道里程为MDmin。在本实施例中隧道埋深最深为24m,里程为K1+170,埋深最浅处为16m,里程为K1+420,提取K1+170、K1+420处地质剖面Geo1070、Geo1420。
根据地下水位获取危险断面:遍历100个地质剖面,得到地下水位最高的隧道里程为MWmax、地下水位最低的隧道里程为MWmin。在本实施例中地下水位最高为-14.5,里程为K1+242,地下水位最低为-22m,里程为K1+320,提取K1+242、K1+320处地质剖面Geo1242、Geo1320。
根据地层变化获取危险断面:本实施例采用地层变化判断条件作为危险断面的选取方式,从上述的100个地质剖面中每5个断面进行一次地层变化判断,地质剖面编号为k=5,10,15…100。
(1)将第k个地质剖面中的地层信息和第k-1环进行对比,k环地质剖切得到的地层中没有a地层,则该环断面作为危险断面之一,并提取该环的地质剖面。
(2)判断第k+1环是否比k环是否多出了一个地层,是则将该环作为危险断面之一。
(3)判断第k环管片位置处竖向Depthj~(Depthj+D)的深度范围内地层数量超过3时,将该环作为危险断面之一。
(4)判断第j环管片埋深处是否处于软岩,是则作为危险断面处理。
按照上述方式进行分析得到危险断面列表。在本实施例提取的20个断面中,编号90断面比编号89断面多一个地层a,提取该环的地质剖面Geo1280。
根据地质模型获取管片环和地质剖面计算模型的边界条件:
将上一步中提取得到的地质剖面Geo1170、Geo1242、Geo1280、Geo1320、Geo1420进行地层边界线提取,得到该地质剖面中每一个地层的边界线和地层参数。如图3所示为Geo1280的地质剖面,图中linek表示地层线,对应地层的重度γk、泊松比υk、弹性模量Ek。
将地质剖面处理为平面单元、管片作为梁单元,根据地层、管片参数设置计算参数,建立有限元分析模型。
所选的地质剖面中地面有超载为20kN/m2,将地面超载作为均布荷载施加在分析模型的地表线上,约束分析模型的刚体位移,进行网格划分,得到地层结构法的二维有限元分析模型,进行分析计算可得地层和管片结构的内力和位移。
上述具体实施方式是用来解释说明本发明,仅为本发明的优选实施例,并不能用于限制本发明的范围。应当理解,凡在本发明的精神和权利要求的保护范围内,对本发明做出的任何修改、等同替换、改进等,都应属于本发明的保护范围。
Claims (7)
1.一种基于三维地质模型的隧道结构一体化分析方法,其特征在于:包括以下步骤:
S1:基于地质勘察、工程设计资料建立三维地质、隧道模型;
S2:将隧道设计信息与三维地质模型关联,通过设定的判断条件自动计算分析,获取危险断面里程;
S3:根据三维地质数据,自动获取危险断面的边界计算条件,分析计算当前隧道断面的位移、应力。
2.根据权利要求1所述的一种基于三维地质模型的隧道结构一体化分析方法,其特征在于:所述S1具体步骤如下:
(1)统一地质、隧道模型坐标系,根据设计资料确定建模范围,地质模型宽度为8D,D为隧道直径,轴线左右各4D,深度为隧道底部以下2D;
(2)根据地层、钻孔、地下水位等地质勘察数据建立三维地质模型,地层数据包括各地层重度γ、泊松比v、水平基床系数K、弹性模量E;钻孔数据包括孔口标高、平面位置坐标、平面位置坐标、各地层标高;
(3)根据隧道轴线、隧道断面尺寸信息在三维地质模型中创建隧道模型。
3.根据权利要求1所述的一种基于三维地质模型的隧道结构一体化分析方法,其特征在于:所述S2包括根据隧道埋深判断危险断面、根据地下水位判断危险断面、根据地层变化判断危险断面;提取剖面的方法为沿隧道轴线,根据隧道里程每Bm提取地质剖面Geoi,盾构法隧道中B取管片环宽,矿山法隧道B取隧道开挖进尺。
4.根据权利要求3所述的一种基于三维地质模型的隧道结构一体化分析方法,其特征在于:所述根据隧道埋深判断危险断面的步骤如下:
(1)遍历Geoi计算隧道埋深Depthi;
(2)根据Depthi计算求得隧道最大埋深Depthmax,记录对应的隧道里程MDmax、地质剖面GeoDmax;
(3)根据Depthi计算求得隧道最小埋深Depthmin,记录对应的隧道里程MDmin、地质剖面GeoDmin。
5.根据权利要求3所述的一种基于三维地质模型的隧道结构一体化分析方法,其特征在于:所述根据地下水位判断危险断面的步骤如下:
(1)遍历Geoi计算地下水位高度WaterHj;
(2)根据WaterHj计算求得最高地下水位WaterHmax,记录对应的隧道里程MWmax、地质剖面GeoWmax;
(3)根据WaterHj计算求得最低地下水位WaterHmin;记录对应的隧道里程MWmin、地质剖面GeoWmin。
6.根据权利要求3所述的一种基于三维地质模型的隧道结构一体化分析方法,其特征在于:所述根据地层变化判断危险断面的步骤如下:
(1)从提取的Geoi中每x个断面进行一次地层变化判断,一般x取1~10之间的整数;
(2)第k个剖面对应的隧道里程为Milek=Miles+k×x×B,k=x,2x,3x…;
(3)逐个分析S1提取剖面地层参数,采用以下方法进行危险断面判断:
①当第k-1个剖面的地层中有地层a,第k个剖面的地层中无地层a,则第k个剖面对应的隧道断面为危险断面;
②当第k-1个剖面的地层中无地层b,第k个剖面的地层中出现了地层b,则第k个剖面对应的隧道断面为危险断面;
③第k个剖面Depthk~(Depthk+D)的深度范围存在软土、杂填土或软硬岩交替地层时,则k个剖面对应的隧道断面为危险断面;
④第k个剖面对应的隧道顶部埋深Depthk,在Depthk~(Depthk+D)的深度范围内地层数量超过nmax时,则第j个剖面对应的隧道断面为危险断面。
7.根据权利要求1所述的一种基于三维地质模型的隧道结构一体化分析方法,其特征在于:所述S3具体步骤如下:
(1)以危险断面地质剖面GeoDmax、GeoDmin、GeoWmax、GeoWmin、Geok为基础;
(2)根据地质剖面图提取各地层边界线,并将地质剖面数据转换为边界荷载;
(3)将地层面作为平面单元,盾构管片作为梁单元,根据地质模型获取各地层γ、υ、K、E,赋值给相应的平面单元材料参数;
(4)施加荷载:考虑自重,采用均布荷载施加地面超载,地层底部约束刚体位移,进行网格划分得到有限元分析模型;
(5)采用地层结构法做二维有限元分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310923252.0A CN116956420A (zh) | 2023-07-26 | 2023-07-26 | 一种基于三维地质模型的隧道结构一体化分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310923252.0A CN116956420A (zh) | 2023-07-26 | 2023-07-26 | 一种基于三维地质模型的隧道结构一体化分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116956420A true CN116956420A (zh) | 2023-10-27 |
Family
ID=88454343
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310923252.0A Pending CN116956420A (zh) | 2023-07-26 | 2023-07-26 | 一种基于三维地质模型的隧道结构一体化分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116956420A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117540937A (zh) * | 2024-01-10 | 2024-02-09 | 深圳市云图华祥科技有限公司 | 一种基于bim的中风化膨胀性泥岩盾构隧道施工分析方法 |
-
2023
- 2023-07-26 CN CN202310923252.0A patent/CN116956420A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117540937A (zh) * | 2024-01-10 | 2024-02-09 | 深圳市云图华祥科技有限公司 | 一种基于bim的中风化膨胀性泥岩盾构隧道施工分析方法 |
CN117540937B (zh) * | 2024-01-10 | 2024-03-22 | 深圳市云图华祥科技有限公司 | 一种基于bim的中风化膨胀性泥岩盾构隧道施工分析方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105184864B (zh) | 一种用于天然地基换填工程量计算的场地地层三维地质结构模型生成方法 | |
CN103758160B (zh) | 一种超深地下连续墙变形自动实时监测装置及其工作方法 | |
CN106644206B (zh) | 一种用于浅埋隧道的围岩压力计算方法 | |
CN109577982A (zh) | 壁式连采连充保水采煤及水资源运移监测、水害预警方法 | |
CN109933867A (zh) | 一种岩土工程勘察方案动态优化方法 | |
CN106372297B (zh) | 砂土岩溶地层中盾构与溶洞安全垂直距离的确定方法 | |
CN109492262A (zh) | 一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法 | |
WO2023185734A1 (zh) | 一种基于三维点云数据库的矿山数字模型建立方法 | |
CN105912798A (zh) | 基于超深基坑抽水的地面沉降智能预警方法及监测系统 | |
CN103559372A (zh) | 隧道开挖施工断面大拱脚crd分层进尺动态分析方法及装置 | |
CN108842821B (zh) | 一种钻爆法修建海底隧道合理埋深的计算方法 | |
CN110630274B (zh) | 一种盾构隧道复合地层组段实施系统及实施方法 | |
CN111119902B (zh) | 一种基于bp神经网络的隧道动态施工方法 | |
CN116956420A (zh) | 一种基于三维地质模型的隧道结构一体化分析方法 | |
CN104809266A (zh) | 一种基于样条曲线的工作面煤层赋存情况精确预测方法 | |
CN108241772A (zh) | 考虑多重因素的裂隙承压含水层中隧道涌水量的确定方法 | |
CN105389451B (zh) | 一种用于无人化采掘工作面的适时精准建模方法 | |
CN110197013B (zh) | 基于Morphing的河床基岩面建模方法 | |
CN115690354A (zh) | 一种基于三维实景数值分析的浅埋隧道施工动态控制方法 | |
CN106372295B (zh) | 砂土岩溶地层中盾构与溶洞安全水平距离的确定方法 | |
CN113536414B (zh) | 基于三维建模的岩质边坡稳定性分析方法、系统及介质 | |
CN101845815A (zh) | 一种软土地层含承压水基坑突涌塑性破坏计算处理方法 | |
CN111737895B (zh) | 一种动态评价露天矿地下采空区顶板稳定性的方法 | |
CN105421335B (zh) | 基于场地超孔隙水压比的水泥搅拌桩复合地基抗液化方法 | |
Ma et al. | Subway tunnel construction settlement analysis based on the combination of numerical Simulation and neural network |
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 |