CN110276102A - 基于岩体节理特征和爆破振动效应修正围岩压力的方法 - Google Patents
基于岩体节理特征和爆破振动效应修正围岩压力的方法 Download PDFInfo
- Publication number
- CN110276102A CN110276102A CN201910417139.9A CN201910417139A CN110276102A CN 110276102 A CN110276102 A CN 110276102A CN 201910417139 A CN201910417139 A CN 201910417139A CN 110276102 A CN110276102 A CN 110276102A
- Authority
- CN
- China
- Prior art keywords
- rock
- explosion
- pressure
- surrounding rock
- joint
- 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.)
- Granted
Links
Classifications
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种基于岩体节理特征和爆破振动效应修正围岩压力的方法，主要解决现有围岩压力测试未考虑节理特征与爆破效应的影响的技术问题。本发明实现步骤如下：（1）基于现场地质素描获取隧道岩体节理产状；（2）利用统计分析获取优势节理产状特征；（3）基于现场爆破参数建立数值模型，获取爆破振动应力波；（4）基于优势节理特征建立离散元计算模型，施加爆破振动应力波，进行隧道开挖，确定围岩松动圈范围；（5）按照设计围岩压力确定松动圈范围，计算修正系数，修正围岩压力。本发明综合考虑岩体节理特征和爆破振动效应对隧道围岩的影响，对隧道支护结构设计提供支持。
Description
技术领域
本发明涉及隧道围岩压力测试技术领域，特别是涉及一种基于岩体节理特征和爆破振动效应修正围岩压力的方法。
背景技术
围岩压力是隧道支护结构承受的主要荷载，也是隧道支护结构的设计的主要依据。目前隧道围岩压力的主要计算方法包括太沙基极限平衡法、土柱理论、普氏平衡拱理论、深浅埋法等。
专利CN106383945A公开了一种隧道及地下空间深、浅埋垂直围岩压力计算方法。根据围岩地质情况及所处围岩分级，确定物理参数，计算洞室垂直地压。
专利CN106682330A公开了一种深埋硐室围岩压力的计算方法。通过构建破坏机制，建立许可的位移场，获得最优上限解，即围岩压力上限值。
由于岩体内部广泛分布着节理裂隙，且隧道开挖过程中爆破振动使得节理裂隙出现一定的扩展，对围岩压力产生一定的影响，而上述计算方法均未考虑节理特征与爆破效应的影响。因此，综合考虑岩体节理特征和爆破振动效应对隧道围岩的影响，修正围岩压力，对隧道支护结构设计及施工安全具有重要的意义。
发明内容
鉴于现有技术的缺陷，本发明提供一种基于岩体节理特征和爆破振动效应修正围岩压力的方法，包括以下步骤：
（1）基于现场地质素描获取隧道岩体节理产状；
（2）利用统计分析获取优势节理产状特征；
（3）基于现场爆破参数建立数值模型，获取爆破振动应力波；
（4）基于优势节理特征建立离散元计算模型，施加爆破振动应力波，进行隧道开挖，确定围岩松动圈范围；
（5）按照设计围岩压力确定松动圈范围，计算修正系数，修正围岩压力。
进一步，步骤（1）所述隧道岩体节理产状主要包括倾角、倾向、迹长、间距、填充物、填充厚度、张开位移、粗糙度等，且每个掌子面节理统计数量不应低于30组。
进一步，步骤（2）所述统计分析，按照以下步骤分析：
a. 根据倾向、倾角做极值等值线图；
b. 根据图中极点分布密度，按照倾向将节理产状分组；
c. 对每组节理产状各参数做频率直方图，利用概率分布模型拟合，确定各组优势节理产状特征。其中，概率分布模型主要包括均匀分布、负指数分布、正态分布、对数正态分布。
进一步，步骤（3）所述爆破数值模型应根据工程现场实际爆破参数建立，炸药爆炸后压力可按公式（1）计算，空气压力可按公式（2）计算。模型边界条件应设置为无反射边界条件，以减小反射波与入射波的相互叠加带来的误差。
（1）
式中，、、、、为状态方程系数，为初始内力，为初始体积。
（2）
式中，~为气体性质相关常数，，、分别为空气的初始密度、密度。
进一步，步骤（4）所述离散元计算模型建立方法，按照以下步骤：
a. 基于步骤（2）优势节理参数，建立隧道模型；
b. 赋予模型材料及节理物理力学参数；
c. 模型划分网格，模型完成。
进一步，步骤（4）所述爆破应力波施加方法，按照以下步骤：
a. 沿隧道内轮廓线每1m选取一个监测单元；
b. 步骤（3）计算完成后，提取每个监测单元爆破应力波；
c. 隧道开挖的同时，在每个监测单元施加对应的爆破应力波。
进一步，步骤（4）所述确定围岩松动圈方法，可根据极限拉应变、剪应力峰值等方法确定松动圈范围。
进一步，步骤（5）所述按照设计围岩压力计算围岩松动范围按照公式（3）计算。
（3）
式中，为隧道毛洞半跨，为隧道毛洞高度，为围岩内摩擦角，为围岩坚固系数。
进一步，步骤（5）所述修正系数按照公式（4）计算
（4）
式中，为考虑节理特征及爆破振动后的松动圈范围。
进一步，步骤（5）所述修正围岩压力可按公式（5）计算
（5）
式中，为围岩重度。
与传统围岩压力计算方法相比，本发明综合考虑了节理特征与爆破振动效应对围岩的影响，修正后的围岩压力值更加符合工程现场实际，对隧道支护设计及施工安全具有重要的意义。
附图说明
图1为本发明流程图。
图2为本发明极点等值线图。
图3为本发明频率直方图及概率模型。
图4为本发明爆破振动应力波。
图5为本发明爆破振动应力波监测单元示意图。
图6为本发明松动圈示意图。
具体实施方式
以下将参照附图，对本发明的具体实施案例进行详细的描述。
图1说明了基于岩体节理特征和爆破振动效应修正围岩压力的方法，包括：
（1）基于现场地质素描获取隧道岩体节理产状，主要包括倾角、倾向、迹长、间距、填充物、填充厚度、张开位移、粗糙度等，且每个掌子面节理统计数量不应低于30组。
（2）根据倾向、倾角做极值等值线图，如图2所示，其中径向的直线表示倾向，沿着逆时针方向，从E至S至W至N方向分别为0至90°至180°至360°，而周向的线表示倾角，从圆点至最外圈表示0°至90°；根据图中极点分布密度，按照倾向将节理产状分组；如图3所示，对每组节理产状各参数做频率直方图，利用概率分布模型拟合，确定各组优势节理产状特征。其中，概率分布模型主要包括均匀分布、负指数分布、正态分布、对数正态分布。
（3）基于工程现场实际爆破参数建立数值模型，获取爆破振动应力波，如图4所示。其中炸药爆炸后压力可按公式（1）计算，空气压力可按公式（2）计算。模型边界条件应设置为无反射边界条件，以减小反射波与入射波的相互叠加带来的误差；
（1）
式中，、、、、为状态方程系数，为初始内力，为初始体积；
（2）
式中，~为气体性质相关常数，，、分别为空气的初始密度、密度。
（4）如图5所示，沿隧道内轮廓线每1m选取一个监测单元，提取每个监测单元爆破应力波；基于优势节理特征建立离散元计算模型，隧道开挖的同时，在每个监测单元施加对应的爆破应力波，根据极限拉应变或剪应力峰值确定围岩松动圈范围；
所述离散元计算模型建立方法（常规方法），按照以下步骤：
a. 基于步骤（2）优势节理参数，建立隧道模型；
b. 赋予模型材料及节理物理力学参数；
c. 模型划分网格，模型完成。
（5）依照公式（3）计算设计围岩压力下松动圈范围，如图6所示，从而按照公式（4）计算修正系数，最后按照公式（5）修正围岩压力；
（3）
式中，为隧道毛洞半跨，为隧道毛洞高度，为围岩内摩擦角，为围岩坚固系数；
（4）
式中，为考虑节理特征及爆破振动后的松动圈范围；
（5）
式中，为围岩重度。
Claims (5)
1.一种基于岩体节理特征和爆破振动效应修正围岩压力的方法，其特征在于：
（1）基于现场地质素描获取隧道岩体节理产状；
（2）利用统计分析获取优势节理产状特征；
（3）基于现场爆破参数建立数值模型，获取爆破振动应力波；
（4）基于优势节理特征建立离散元计算模型，施加爆破振动应力波，进行隧道开挖，确定围岩松动圈范围；
（5）按照设计围岩压力确定松动圈范围，计算修正系数，修正围岩压力。
2.如权利要求1所述的基于岩体节理特征和爆破振动效应修正围岩压力的方法，其特征在于：步骤（2）所述统计分析，按照以下步骤分析：
a. 根据倾向、倾角做极值等值线图；
b. 根据图中极点分布密度，按照倾向将节理产状分组；
c. 对每组节理产状各参数做频率直方图，利用概率分布模型拟合，确定各组优势节理产状特征。
3.如权利要求2所述的基于岩体节理特征和爆破振动效应修正围岩压力的方法，其特征在于：所述概率分布模型包括均匀分布、负指数分布、正态分布、对数正态分布。
4.如权利要求1所述的基于岩体节理特征和爆破振动效应修正围岩压力的方法，其特征在于：步骤（4）所述爆破应力波施加方法，按照以下步骤：
a. 沿隧道内轮廓线每1m选取一个监测单元；
b. 步骤（3）计算完成后，提取每个监测单元爆破应力波；
c. 隧道开挖的同时，在每个监测单元施加对应的爆破应力波。
所述确定围岩松动圈范围，根据极限拉应变、剪应力峰值确定松动圈范围。
5.如权利要求1所述的基于岩体节理特征和爆破振动效应修正围岩压力的方法，其特征在于：步骤（5）所述按照设计围岩压力确定松动圈范围按照公式（3）获得：
（3）
式中，为隧道毛洞半跨，为隧道毛洞高度，为围岩内摩擦角，为围岩坚固系数；
所述修正系数按照公式（4）获得：
（4）
式中，为考虑节理特征及爆破振动后的松动圈范围；
所述修正围岩压力可按公式（5）获得：
（5）
式中，为围岩重度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910417139.9A CN110276102B (zh) | 2019-05-20 | 2019-05-20 | 基于岩体节理特征和爆破振动效应修正围岩压力的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910417139.9A CN110276102B (zh) | 2019-05-20 | 2019-05-20 | 基于岩体节理特征和爆破振动效应修正围岩压力的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110276102A true CN110276102A (zh) | 2019-09-24 |
| CN110276102B CN110276102B (zh) | 2023-01-17 |
Family
ID=67959443
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910417139.9A Active CN110276102B (zh) | 2019-05-20 | 2019-05-20 | 基于岩体节理特征和爆破振动效应修正围岩压力的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110276102B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110826254A (zh) * | 2019-11-27 | 2020-02-21 | 神华准格尔能源有限责任公司 | 露天矿边坡岩体节理分布参数测定方法、存储介质和系统 |
| CN111581785A (zh) * | 2020-04-18 | 2020-08-25 | 中铁二院工程集团有限责任公司 | 高驼峰应力围岩分级方法 |
| CN112182703A (zh) * | 2020-09-27 | 2021-01-05 | 中电建路桥集团有限公司 | 一种基于离散单元法联合荷载结构法的类岩堆体隧道结构分析方法 |
| CN112649086A (zh) * | 2020-11-10 | 2021-04-13 | 安徽理工大学 | 一种改进的节理岩体爆破模型试验震动监测系统和方法 |
| CN113190891A (zh) * | 2021-02-10 | 2021-07-30 | 国机集团科学技术研究院有限公司 | 一种古建筑木结构残损抗振性能的评估试验方法 |
| CN113642145A (zh) * | 2021-06-23 | 2021-11-12 | 中铁十六局集团第三工程有限公司 | 计算爆破作用下隧道关键块体稳定性的方法及系统 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105787206A (zh) * | 2016-03-23 | 2016-07-20 | 同济大学 | 一种公路隧道锚杆支护设计方法 |
| US20180010316A1 (en) * | 2016-07-06 | 2018-01-11 | PoweChina Huadong Engineering Corporation Limited | Yieldable construction method for early releasing surrounding rock deformation on weak counter-inclined slope |
-
2019
- 2019-05-20 CN CN201910417139.9A patent/CN110276102B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105787206A (zh) * | 2016-03-23 | 2016-07-20 | 同济大学 | 一种公路隧道锚杆支护设计方法 |
| US20180010316A1 (en) * | 2016-07-06 | 2018-01-11 | PoweChina Huadong Engineering Corporation Limited | Yieldable construction method for early releasing surrounding rock deformation on weak counter-inclined slope |
Non-Patent Citations (1)
| Title |
|---|
| J.H. YANG 等: "Dynamic stress adjustment and rock damage during blasting excavation in a deep-buried circular tunnel", 《TUNNELLING AND UNDERGROUND SPACE TECHNOLOGY》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110826254A (zh) * | 2019-11-27 | 2020-02-21 | 神华准格尔能源有限责任公司 | 露天矿边坡岩体节理分布参数测定方法、存储介质和系统 |
| CN111581785A (zh) * | 2020-04-18 | 2020-08-25 | 中铁二院工程集团有限责任公司 | 高驼峰应力围岩分级方法 |
| CN111581785B (zh) * | 2020-04-18 | 2024-08-27 | 中铁二院工程集团有限责任公司 | 高驼峰应力围岩分级方法 |
| CN112182703A (zh) * | 2020-09-27 | 2021-01-05 | 中电建路桥集团有限公司 | 一种基于离散单元法联合荷载结构法的类岩堆体隧道结构分析方法 |
| CN112182703B (zh) * | 2020-09-27 | 2023-08-25 | 中电建路桥集团有限公司 | 一种基于离散单元法联合荷载结构法的类岩堆体隧道结构分析方法 |
| CN112649086A (zh) * | 2020-11-10 | 2021-04-13 | 安徽理工大学 | 一种改进的节理岩体爆破模型试验震动监测系统和方法 |
| CN113190891A (zh) * | 2021-02-10 | 2021-07-30 | 国机集团科学技术研究院有限公司 | 一种古建筑木结构残损抗振性能的评估试验方法 |
| CN113190891B (zh) * | 2021-02-10 | 2023-09-29 | 国机集团科学技术研究院有限公司 | 一种古建筑木结构残损抗振性能的评估试验方法 |
| CN113642145A (zh) * | 2021-06-23 | 2021-11-12 | 中铁十六局集团第三工程有限公司 | 计算爆破作用下隧道关键块体稳定性的方法及系统 |
| CN113642145B (zh) * | 2021-06-23 | 2024-05-10 | 中铁十六局集团第三工程有限公司 | 计算爆破作用下隧道关键块体稳定性的方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110276102B (zh) | 2023-01-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110276102A (zh) | 基于岩体节理特征和爆破振动效应修正围岩压力的方法 | |
| CN105224742B (zh) | 一种分析爆破作用下台阶边坡稳定性的方法 | |
| Baryakh et al. | Sinkhole formation mechanism | |
| CN109583116B (zh) | 基于多源实测信息的地下工程初始应力场动态反演方法 | |
| CN108763725A (zh) | 基于采空区压实效应的工作面开采全过程数值模拟方法 | |
| CN110135030A (zh) | 一种采空区地表沉降的预测方法 | |
| CN108362562B (zh) | 一种复合地层浅埋小净距隧道围岩压力的确定方法 | |
| CN106814391A (zh) | 基于菲涅尔体层析反演的地面微地震事件定位方法 | |
| CN113094778B (zh) | 高地应力互层软岩隧道破坏机理及施工控制研究方法 | |
| CN113255175A (zh) | 一种层状软岩爆破数值模拟方法 | |
| CN111553107A (zh) | 可液化场地桩基随机地震反应分析与安全评价方法 | |
| CN115062368A (zh) | 一种隧道长管棚数字孪生体和精细建模系统及方法 | |
| CN111101999B (zh) | 一种确定松软煤层抽采钻孔最佳塑性区范围的方法 | |
| CN108805365B (zh) | 一种非均质页岩储层裂缝诱导应力预测方法 | |
| Wood et al. | Design of large diameter monopiles in chalk atWestermost Rough offshore wind farm | |
| CN114297824A (zh) | 一种深部高应力硬岩板裂化岩爆释能支护体系设计方法 | |
| CN110543680B (zh) | 一种局部地基沉陷引起地埋管线变形及内力计算方法 | |
| CN106372295B (zh) | 砂土岩溶地层中盾构与溶洞安全水平距离的确定方法 | |
| Gospodarikov et al. | Evaluation of deformation characteristics of brittle rocks beyond the limit of strength in the mode of uniaxial servohydraulic loading | |
| CN115081302A (zh) | 支护构件与硐室围岩接触及相互作用的模拟方法和系统 | |
| CN109612356B (zh) | 基坑开挖爆破时临近混凝土管道动力响应特性的获取方法 | |
| Adachi et al. | Behavior and simulation of soil tunnel with thin cover | |
| CN112883614B (zh) | 一种基于数值模拟的岩溶地层盾构隧道溶洞处理范围判断方法 | |
| CN110552716B (zh) | 一种圆形盾构隧道衬砌结构的拼装方法 | |
| CN109930634B (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 |