CN114820949B - 深切河谷区地下洞室群地应力场的预测方法 - Google Patents
深切河谷区地下洞室群地应力场的预测方法 Download PDFInfo
- Publication number
- CN114820949B CN114820949B CN202210613271.9A CN202210613271A CN114820949B CN 114820949 B CN114820949 B CN 114820949B CN 202210613271 A CN202210613271 A CN 202210613271A CN 114820949 B CN114820949 B CN 114820949B
- Authority
- CN
- China
- Prior art keywords
- ground stress
- stress field
- cavern group
- group
- calculation model
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- 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
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/04—Architectural design, interior design
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/20—Hydro energy
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)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明深切河谷区地下洞室群地应力场的预测方法,涉及地应力的预测,通过实测地应力场进行三维地应力场反演,获得三维数值计算模型的边界位移和重力加速度,利用所述三维数值计算模型获得初始洞室群地应力场,然后利用三维数值计算模型对洞室群第一层进行开挖模拟,预测围岩的变形破坏特征,然后,观测开挖完成后的围岩的变形破坏特征,以实际围岩的变形破坏特征反算三维数值计算模型的边界条件,以此来获得精确的三维数值计算模型的边界条件,从而获得准确的洞室群地应力场,解决了现有技术对深切河谷区地下洞室群地应力场预测困难的问题,本发明适用于洞室群地应力场的预测。
Description
技术领域
本发明涉及地应力的预测,特别涉及深切河谷区地下洞室群地应力场的预测方法。
背景技术
深切河谷是指一种长期受流水侵蚀作用影响形成的河谷地带。
大型水电工程地下洞室群的稳定性受到所处区域地应力场的影响,准确的地下洞室群初始地应力场是进行洞室群开挖支护方案设计的基础。水电工程地下洞室群初始地应力场受到河谷、重力和构造运动,以及各级结构面的影响,准确地预测水电工程地下洞室群初始地应力场是一个重要而困难的任务。
发明内容
本发明所解决的技术问题:提供一种深切河谷区地下洞室群地应力场的预测方法解决现有技术对深切河谷区地下洞室群地应力场预测困难的问题。
本发明解决上述技术问题采用的技术方案:深切河谷区地下洞室群地应力场的预测方法,包括以下步骤:
S01、建立地下洞室群的三维数值计算模型,并利用地下洞室群所在地的实测地应力数据进行三维地应力场反演,所述反演结果为三维数值计算模型的边界位移组合,对三维数值计算模型施加所述边界位移组合,利用所述三维数值计算模型获得初始洞室群地应力场;
S02、依据初始洞室群地应力场,设计地下洞室群第一层的开挖方案,并利用所述三维数值计算模型对洞室群第一层进行开挖模拟,预测洞室群第一层开挖完成后围岩的变形破坏特征;
S03、根据地下洞室群第一层的开挖方案,进行开挖,在开挖结束后,观测围岩的变形破坏特征,如果观测到的围岩的变形破坏特征与预测的洞室群第一层开挖完成后洞室围岩的变形破坏特征在预设的误差内,则地下洞室群地应力场为所述初始洞室群地应力场,否则,进行洞室群初始地应力场二次反演;
所述二次反演过程为:基于正交试验设计,构建不同的边界位移组合样本,对三维数值计算模型施加不同的边界位移样本,模拟洞室群第一层开挖过程,得到每个位移组合样本对应的围岩的变形破坏特征,构建边界位移组合与围岩的变形破坏特征的神经网络,以观测的围岩的变形破坏特征作为反演目标,获得新的边界位移组合,再以所述新的边界位移组合施加到三维数值计算模型,得到的洞室群区域地应力场。
进一步的,所述三维数值计算模型包括地下洞室群、河谷地形、剥蚀分层以及二、三级结构面。
进一步的,所述边界位移组合包括重力加速度。
进一步的,所述围岩的变形破坏特征包括围岩片帮的位置和深度。
进一步的,所述三维数值计算模型的尺寸边界距离洞室群中最近的洞室的距离不小于该洞室尺寸的三倍。
进一步的,在开挖模拟过程中,主应力差的等值线与洞壁轮廓所围区域为片帮,等值线与洞壁轮廓的距离为深度,所述主应力差是第一主应力与第三主应力之差,所述主应力差的等值线是所述主应力差等于起裂应力时对应的等值线。
进一步的,所述起裂应力的获取方式为:通过对岩样进行单轴压缩试验,当岩样起裂时,对应的轴向应力为起裂应力。
本发明的有益效果:本发明深切河谷区地下洞室群地应力场的预测方法,通过实测地应力进行三维地应力场反演,获得三维数值计算模型的边界位移组合,对三维数值计算模型施加所述边界位移组合,利用所述三维数值计算模型获得初始洞室群地应力场,然后利用三维数值计算模型对洞室群第一层进行开挖模拟,预测围岩的变形破坏特征,在实际开挖完成后,观测围岩的变形破坏特征是否在预设的误差内,如果在,则初始洞室群地应力场为地下洞室群地应力场,否则,进行二次反演,基于正交试验设计,构建不同的边界位移组合样本,对三维数值计算模型施加不同的边界位移样本,模拟洞室群第一层开挖过程,得到每个位移组合样本对应的围岩的变形破坏特征,构建边界位移组合与围岩的变形破坏特征的神经网络,以观测的围岩的变形破坏特征作为反演目标,获得新的边界位移组合,再以所述新的边界位移组合施加到三维数值计算模型,得到的洞室群区域地应力场,解决了深切河谷区地下洞室群地应力场预测困难的问题。
附图说明
附图1是本发明深切河谷区地下洞室群地应力场的预测方法的流程示意图。
具体实施方式
本发明深切河谷区地下洞室群地应力场的预测方法,如附图1所示,包括以下步骤:
S01、建立地下洞室群的三维数值计算模型,并利用地下洞室群所在地的实测地应力数据进行三维地应力场反演,所述反演结果为三维数值计算模型的边界位移组合,对三维数值计算模型施加所述边界位移组合,利用所述三维数值计算模型获得初始洞室群地应力场;
具体的,所述三维数值计算模型包括地下洞室群、河谷地形、剥蚀分层以及二、三级结构面;边界位移组合包括重力加速度;围岩的变形破坏特征包括围岩片帮的位置和深度;三维数值计算模型的尺寸边界距离洞室群中最近的洞室的距离不小于该洞室尺寸的三倍,在利用所述三维数值计算模型获得初始洞室群地应力场的过程中,采用的本构模型即可以为弹性模型也可以为弹塑性模型。
S02、依据初始洞室群地应力场,设计地下洞室群第一层的开挖方案,并利用所述三维数值计算模型对洞室群第一层进行开挖模拟,预测洞室群第一层开挖完成后围岩的变形破坏特征;
具体的,在开挖模拟过程中,采用的本构模型为弹性模型,主应力差的等值线与洞壁轮廓所围区域为片帮,等值线与洞壁轮廓的距离为深度,所述主应力差是第一主应力与第三主应力之差,所述主应力差的等值线是所述主应力差等于起裂应力时对应的等值线;所述起裂应力的获取方式为:通过对岩样进行单轴压缩试验,当岩样起裂时,对应的轴向应力为起裂应力,所述岩样是以地下洞室群区域的岩芯制作成的标准岩样,在判断岩样起裂时,可以采用声发射检测技术。
S03、根据地下洞室群第一层的开挖方案,进行开挖,在开挖结束后,观测围岩的变形破坏特征,如果观测到的围岩的变形破坏特征与预测的洞室群第一层开挖完成后洞室围岩的变形破坏特征在预设的误差内,则地下洞室群地应力场为所述初始洞室群地应力场,否则,进行洞室群初始地应力场二次反演;
所述二次反演过程为:基于正交试验设计,构建不同的边界位移组合样本,对三维数值计算模型施加不同的边界位移样本,模拟洞室群第一层开挖过程,得到每个位移组合样本对应的围岩的变形破坏特征,采用遗传算法-智能反演的方法,构建边界位移组合与围岩的变形破坏特征的神经网络,以观测的围岩的变形破坏特征作为反演目标,获得新的边界位移组合,再以所述新的边界位移组合和所述重力加速度施加到三维数值计算模型,得到的洞室群区域地应力场。
具体的,所述预设的误差是指观测到的围岩片帮的位置与S02中预测的位置偏差在设定的范围内,且观测到的相近围岩片帮的深度差不超过设定差值,如果观测到的围岩的变形破坏特征在预设的误差内,则说明初始洞室群地应力场符合实际要求,即可以将初始洞室群地应力场作地下洞室群地应力场,进行后续的开挖工作,否则,初始洞室群地应力场不符合实际要求,需要重新计算洞室群地应力场,所述新的边界位移组合是以实际围岩的变形破坏特征反算三维数值计算模型的边界条件,以此来获得精确的三维数值计算模型的边界条件,避免了S01中三维地应力反演结果误差引起的初始洞室群地应力场不准确的问题,从而可以避免后续的开挖过程中出现需要调整支护措施等问题,在二次反演过程中,采用的本构模型为弹性模型。
Claims (7)
1.深切河谷区地下洞室群地应力场的预测方法,其特征在于,包括以下步骤:
S01、建立地下洞室群的三维数值计算模型,并利用地下洞室群所在地的实测地应力数据进行三维地应力场反演,所述反演结果为三维数值计算模型的边界位移组合,对三维数值计算模型施加所述边界位移组合,利用所述三维数值计算模型获得初始洞室群地应力场;
S02、依据初始洞室群地应力场,设计地下洞室群第一层的开挖方案,并利用所述三维数值计算模型对洞室群第一层进行开挖模拟,预测洞室群第一层开挖完成后围岩的变形破坏特征;
S03、根据地下洞室群第一层的开挖方案,进行开挖,在开挖结束后,观测围岩的变形破坏特征,如果观测到的围岩的变形破坏特征与预测的洞室群第一层开挖完成后洞室围岩的变形破坏特征在预设的误差内,则地下洞室群地应力场为所述初始洞室群地应力场,否则,进行洞室群初始地应力场二次反演;
所述二次反演过程为:基于正交试验设计,构建不同的边界位移组合样本,对三维数值计算模型施加不同的边界位移样本,模拟洞室群第一层开挖过程,得到每个位移组合样本对应的围岩的变形破坏特征,构建边界位移组合与围岩的变形破坏特征的神经网络,以观测的围岩的变形破坏特征作为反演目标,获得新的边界位移组合,再以所述新的边界位移组合施加到三维数值计算模型,得到的洞室群区域地应力场。
2.根据权利要求1所述的深切河谷区地下洞室群地应力场的预测方法,其特征在于,所述三维数值计算模型包括地下洞室群、河谷地形、剥蚀分层以及二、三级结构面。
3.根据权利要求1所述的深切河谷区地下洞室群地应力场的预测方法,其特征在于,所述边界位移组合包括重力加速度。
4.根据权利要求1所述的深切河谷区地下洞室群地应力场的预测方法,其特征在于,所述围岩的变形破坏特征包括围岩片帮的位置和深度。
5.根据权利要求1所述的深切河谷区地下洞室群地应力场的预测方法,其特征在于,所述三维数值计算模型的尺寸边界距离洞室群中最近的洞室的距离不小于该洞室尺寸的三倍。
6.根据权利要求1-5任一项所述的深切河谷区地下洞室群地应力场的预测方法,其特征在于,在开挖模拟过程中,主应力差的等值线与洞壁轮廓所围区域为片帮,等值线与洞壁轮廓的距离为深度,所述主应力差是第一主应力与第三主应力之差,所述主应力差的等值线是所述主应力差等于起裂应力时对应的等值线。
7.根据权利要求6所述的深切河谷区地下洞室群地应力场的预测方法,其特征在于,所述起裂应力的获取方式为:通过对岩样进行单轴压缩试验,当岩样起裂时,对应的轴向应力为起裂应力。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210613271.9A CN114820949B (zh) | 2022-05-31 | 2022-05-31 | 深切河谷区地下洞室群地应力场的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210613271.9A CN114820949B (zh) | 2022-05-31 | 2022-05-31 | 深切河谷区地下洞室群地应力场的预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114820949A CN114820949A (zh) | 2022-07-29 |
CN114820949B true CN114820949B (zh) | 2023-04-28 |
Family
ID=82519619
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210613271.9A Active CN114820949B (zh) | 2022-05-31 | 2022-05-31 | 深切河谷区地下洞室群地应力场的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114820949B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103233741A (zh) * | 2013-04-18 | 2013-08-07 | 中国水电顾问集团成都勘测设计研究院 | 定量化评价超大型地下洞室群施工期围岩稳定性的方法 |
CN105675384A (zh) * | 2016-01-11 | 2016-06-15 | 中国地质大学(北京) | 储油环境下深部硬石膏洞室群稳定性评价及布局优化方法 |
CN113221228A (zh) * | 2021-06-04 | 2021-08-06 | 中国电建集团成都勘测设计研究院有限公司 | 水电站地下洞群围岩力学参数反演方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130191044A1 (en) * | 2011-03-21 | 2013-07-25 | Schlumberger Technology Corporation | Method and system for presenting stress domain |
FR3023641A1 (zh) * | 2014-07-11 | 2016-01-15 | Schlumberger Services Petrol |
-
2022
- 2022-05-31 CN CN202210613271.9A patent/CN114820949B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103233741A (zh) * | 2013-04-18 | 2013-08-07 | 中国水电顾问集团成都勘测设计研究院 | 定量化评价超大型地下洞室群施工期围岩稳定性的方法 |
CN105675384A (zh) * | 2016-01-11 | 2016-06-15 | 中国地质大学(北京) | 储油环境下深部硬石膏洞室群稳定性评价及布局优化方法 |
CN113221228A (zh) * | 2021-06-04 | 2021-08-06 | 中国电建集团成都勘测设计研究院有限公司 | 水电站地下洞群围岩力学参数反演方法 |
Non-Patent Citations (2)
Title |
---|
Annan Jiang等.Back-Analysis of Stress Field for Underground Powerhouse Using Difference Evolutionary Self-Adaptive BP.《2008 Fourth International Conference on Natural Computation》.2008,第314-318页. * |
裴书锋等.考虑洞室岩体应力型破坏特征的局部地应力反演方法及应用.《岩土力学》.2020,第4093-4104页. * |
Also Published As
Publication number | Publication date |
---|---|
CN114820949A (zh) | 2022-07-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109271738B (zh) | 一种用于获取巷道围岩Weibull分布参数的数值反演方法 | |
CN106383172B (zh) | 一种基于能量释放系数的围岩损伤预测方法 | |
CN110005407B (zh) | 盐穴储气库稳定性评价方法 | |
CN111553101B (zh) | 一种隧道开挖上覆岩层开裂预报方法及围岩支护方法 | |
CN112035937B (zh) | 一种隧道施工围岩变形分析及信息化施工方法 | |
CN111504252B (zh) | 一种长距离隧洞膨胀性围岩变形超前预测预报方法 | |
CN106529092A (zh) | 钻井套管的剪切变形力有限元计算方法 | |
CN111119902B (zh) | 一种基于bp神经网络的隧道动态施工方法 | |
CN115267905B (zh) | 一种复杂构造区钻井工程中裂缝性井漏预测方法 | |
CN110096824B (zh) | 盐穴储油库稳定性评价方法 | |
Rimoy | Ageing and axial cyclic loading studies of displacement piles in sands | |
CN115238553B (zh) | 一种地埋管线渗漏浸蚀的危险区域划分方法和系统 | |
CN113221462A (zh) | 一种pso-de算法和神经网络联合的基坑土体参数反演方法 | |
Le et al. | Subsurface ground movements due to circular shaft construction | |
CN110263448B (zh) | 一种双隧道开挖对邻近地埋管线影响的预测方法 | |
CN114820949B (zh) | 深切河谷区地下洞室群地应力场的预测方法 | |
CN109711063A (zh) | 一种适用于上软下硬地层的车站暗挖进尺调整方法与装置 | |
Liang et al. | Displacement back analysis for a high slope of the Dagangshan hydroelectric power station based on BP neural network and particle swarm optimization | |
Yin et al. | Spudcan penetration simulation using the coupled Eulerian-Lagrangian method with thermo-mechanical coupled analysis | |
Trinh et al. | Continued rock stress and displacement measurements combined with numerical modeling as an active, realistic rock engineering tool | |
CN116680934A (zh) | 基于卸荷率的地层开挖引起下部隧道竖向变形计算方法 | |
CN113495013A (zh) | 一种用于多场耦合下沉管隧道结构变形模型的实验方法 | |
Wu et al. | In Situ Monitoring and Numerical Experiments on Vertical Deformation Profiles of Large‐Scale Underground Caverns in Giant Hydropower Stations | |
CN115081294B (zh) | 高水头钢筋混凝土衬砌压力隧洞充水过程分析方法 | |
Vrouwenvelder et al. | General report 4: Reliability of predictions |
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 |