CN109492262A - 一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法 - Google Patents
一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法 Download PDFInfo
- Publication number
- CN109492262A CN109492262A CN201811201014.4A CN201811201014A CN109492262A CN 109492262 A CN109492262 A CN 109492262A CN 201811201014 A CN201811201014 A CN 201811201014A CN 109492262 A CN109492262 A CN 109492262A
- Authority
- CN
- China
- Prior art keywords
- crack
- numerical simulation
- uniform distribution
- rock
- 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.)
- Pending
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000004458 analytical method Methods 0.000 title claims abstract description 23
- 238000009828 non-uniform distribution Methods 0.000 title claims abstract description 22
- 239000011435 rock Substances 0.000 claims abstract description 78
- 230000008093 supporting effect Effects 0.000 claims abstract description 30
- 238000009826 distribution Methods 0.000 claims abstract description 22
- 238000013461 design Methods 0.000 claims abstract description 15
- 239000000463 material Substances 0.000 claims abstract description 15
- 238000009412 basement excavation Methods 0.000 claims abstract description 14
- 238000012544 monitoring process Methods 0.000 claims abstract description 12
- 238000013097 stability assessment Methods 0.000 claims abstract description 5
- 238000005315 distribution function Methods 0.000 claims description 16
- 238000005553 drilling Methods 0.000 claims description 7
- 238000006073 displacement reaction Methods 0.000 claims description 5
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 4
- 239000000284 extract Substances 0.000 claims description 4
- 238000009472 formulation Methods 0.000 claims description 4
- 230000005484 gravity Effects 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 4
- 230000003014 reinforcing effect Effects 0.000 claims description 4
- 238000005065 mining Methods 0.000 abstract description 3
- 238000005259 measurement Methods 0.000 abstract 1
- 230000006870 function Effects 0.000 description 7
- 238000011439 discrete element method Methods 0.000 description 5
- 238000011161 development Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000004744 fabric Substances 0.000 description 3
- 230000003245 working effect Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000001687 destabilization Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
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
- 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
- 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
- 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/06—Power analysis or power optimisation
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法，涉及采矿技术领域，解决了通过数值模拟方法分析裂隙非均匀分布的巷道围岩稳定性从而为二次支护提供依据的技术问题，该方法的步骤包括：(1)利用有限元模拟软件建立数值模型；(2)非均匀岩体概率分布参数设定；(3)测量并计算巷道围岩地质强度指标；(4)对数值模型的材料参数进行赋值；(5)设定边界条件，模型开挖计算；(6)提取监测数据，进行巷道稳定性评估和二次支护设计。本方法结合了现场实测和数值模拟，模拟了非均质岩体条件下的巷道开挖变形情况，从而更好的分析了巷道围岩的稳定性，另外该方法还具有操作简单，计算准确，可靠性高等优点。
Description
技术领域
本发明涉及采矿技术领域，尤其是一种利用数值模拟和现场监测方法分析非均匀分布裂隙巷道稳定性的方法。
背景技术
岩土工程的隧道、矿业、水电、核电等领域，尤其是地下工程中，岩体裂隙发育具有随机性、形态多样性和空间组合复杂性等特点，从而使得工程岩体的力学参数如弹性模量等具有较高的非均质性，岩体参数的非均质性对工程稳定产生重要影响，严重时会造成井下巷道等支护工程的意外失效。研究非均匀裂隙围岩稳定性控制是岩土工程领域的重要研究课题，对保证安全、高效的施工和生产具有重要指导意义，但是常规的数值模拟方法往往将岩体视为均匀介质，忽略岩体的非均匀性对工程的影响，因此无法对非均匀裂隙发育巷道围岩控制提供准确可靠的支持。
目前，分析围岩失稳机理、研究围岩控制对策、设计支护方案等领域，岩土工程中广泛采用的数值模拟方法主要有：(1)有限元法，有限元法可以将不规则的网格单元拼接成不同形状的研究对象，有限元法应用性强且计算精度较高，但是计算过程复杂、计算量庞大，因此在模拟复杂的工程问题时有限元法存在一定的局限。(2)有限差分法，有限差分法用均匀的网格划分求解域，用有限个离散点的差分代替微分，具有建模方便计算效率高的优点，尤其适用于范围和规模较大的工程模拟。(3)离散元法，离散元法是解决不连续介质问题的数值模拟方法。该方法将裂隙岩体视为由完整岩体单元和节理面组成，单元可以发生平移、转动和变形，节理面允许被压缩、分离或滑动。离散元法适用于模拟节理面空间分布结构与节理面参数明确的裂隙岩体。但是节理空间分布形态以及参数很难进行合理的估算，而且计算维问题时计算量与存储量巨大。有限元法和有限差分法均采用连续介质力学。离散元法中节理的力学参数需要用户自行定义，但这部分力学参数很难准确估算，节理参数的选取不当会极大影响离散元法的计算结果，无法对工程实际进行准确的模拟和仿真。在不专门针对裂隙发育、扩展的围岩稳定性数值模拟中，基于连续介质的有限元法和有限差分法凭借其较高的计算效率和准确的计算结果较好。
由于岩体中裂隙的产状和发育程度影响岩体的强度与变形参数，岩体裂隙发育状态严重影响巷道的变形规律，但是由于岩体裂隙发育的非均匀性造成不同区域巷道围岩力学参数差异性明显，围岩力学参数的非均匀分布对围岩失稳机理、支护方案设计等产生的影响是不可忽略的。因此，在连续介质数值模拟方法的基础上，并在保证计算准确性和高效性的前提下，考虑围岩裂隙的非均匀分布对岩体参数的影响，并提出反演岩体非均质性参数的方法，从而进一步提高数值模拟过程的合理性和结果可靠性。
发明内容
为提高围岩裂隙的非均匀分布条件下的巷道围岩稳定性数值模拟的合理性和可靠性，本发明提供了一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法，具体技术方案如下。
一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法，其特征在于，步骤包括：
步骤一.利用有限元模拟软件建立三维数值模型，包括数值模型的尺寸设计和网格划分；
步骤二.非均匀岩体概率分布参数设定，包括确定概率分布函数，模拟岩体的非均质性；
步骤三.测量并计算巷道围岩的地质强度指标，包括钻探取芯获取岩石质量指标和节理条件，地质强度指标的概率分布拟合；
步骤四.对数值模型的材料参数进行赋值，所述材料参数包括弹性模量、泊松比、内聚力、内摩擦角、抗拉强度、残余内聚力和塑性应变量；
步骤五.设定边界条件，模型开挖计算，所述边界条件包括施加重力模拟初始应力，设置应力边界条件和位移边界条件，所述开挖计算使用有限差分计算方法；
步骤六.提取监测数据，进行巷道稳定性评估和二次支护设计。
优选的是，步骤一中的三维数值模型为长方体，在三维数值模型中部开挖巷道，巷道中部划分为顶板监测区域。
优选的是，步骤二中概率分布函数取韦布尔概率分布函数：
其中u是单个元素的力学参数值，u0是尺度参数，m为材料的均质性系数。
进一步优选的是，步骤二中使用蒙特卡洛法模拟岩体的非均质性，包括：A.在[μmin，μmax]之间随机取一个值μi，计算μi的概率密度函数f(μi)；B.在f(μmin)与f(μmax)之间随机产生一个概率密度函数f(μ)，当f(μ)<f(μi)时，导出μi，C.继续进行迭代计算，最终得到μi的数据组，μi的数据组服从韦布尔概率分布规律。
优选的是，步骤三中的地质强度指标GIS与岩石质量指标RQD和节理条件JCond89的关系为：GSI＝1.5JCond89+RQD/2。
进一步优选的是，岩石质量指标RQD取钻探获取的岩芯完整度来表示，具体为长度大于10cm的岩芯占钻孔总长度的比值；所述节理条件JCond89的取值根据巷道围岩裂隙长度、宽度、粗糙度、充填类型和风化程度确定。
优选的是，步骤四中弹性模量E和地质强度指标GIS之间的关系为：
其中σm为岩石单轴抗压强度。
优选的是，步骤五中模型开挖使用有限差分法进行迭代计算，当模型的最大不平衡力当低于1e-5时模拟运算结束。
优选的是，步骤六中二次支护设计包括，在巷道变形值达到120mm以上的区域，采取注浆加固和增加锚杆预紧力的支护措施。
本发明的有益效果包括：
(1)本方法结合现场监测和数值模拟对非均质岩体的巷道裂隙分布和顶板下沉规律进行模拟，模拟结果和现场监测相互验证，利用韦布尔概率分布函数模拟岩体的非均质性，提高了模拟的合理性和可靠性。
(2)根据钻探结果确定地质强度指标从而更好的还原了工程现场的巷道围岩条件，保证了模拟结果的准确性，尤其是对地质强度指标GIS的统计值进行概率分布拟合，确定了顶板的模拟参数，从而能够更准确的计算巷道顶板下沉变形值。
(3)根据巷道模拟结果，提取巷道变形量监测数据，对巷道围岩的稳定性和可靠性进行评估分析，从而进一步的指导巷道二次支护保证了巷道的安全和支护的有效性。
另外本方法还具有操作简单，计算准确高效，支护可靠性高等优点。
附图说明
图1是利用数值模拟分析非均匀分布裂隙巷道稳定性的方法流程结构图；
图2是三维数值模型结构示意图；
图3是巷道顶底板岩性特征和厚度柱状图；
图4是不同均质性系数的韦布尔分布曲线图；
图5是韦布尔分布参数拟合示意图；
图6是数值模拟结果巷道顶板下沉示意图。
具体实施方式
结合图1至图6所示，本发明提供的一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法，具体实施方式如下。
实施例1
一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法，其特征在于，步骤包括：
步骤一.利用有限元模拟软件建立三维数值模型，包括数值模型的尺寸设计和网格划分；
在本步骤中的三维数值模型为长方体，在三维数值模型中部开挖巷道，巷道中部划分为顶板监测区域。
步骤二.非均匀岩体概率分布参数设定，包括确定概率分布函数，模拟岩体的非均质性；
在本步骤中概率分布函数取韦布尔概率分布函数：
其中u是单个元素的力学参数值，u0是尺度参数，m为材料的均质性系数。
其中使用蒙特卡洛法模拟岩体的非均质性，包括：A.在[μmin，μmax]之间随机取一个值μi，计算μi的概率密度函数f(μi)；B.在f(μmin)与f(μmax)之间随机产生一个概率密度函数f(μ)，当f(μ)<f(μi)时，导出μi，C.继续进行迭代计算，最终得到μi的数据组，μi的数据组服从韦布尔概率分布规律。利用韦布尔概率分布函数模拟岩体的非均质性，提高了模拟的合理性和可靠性。
步骤三.测量并计算巷道围岩的地质强度指标，包括钻探取芯获取岩石质量指标和节理条件，地质强度指标的概率分布拟合；
在本步骤中的地质强度指标GIS与岩石质量指标RQD和节理条件JCond89的关系为：GSI＝1.5JCond89+RQD/2。岩石质量指标RQD取钻探获取的岩芯完整度来表示，具体为长度大于10cm的岩芯占钻孔总长度的比值。其中节理条件JCond89的取值根据巷道围岩裂隙长度、宽度、粗糙度、充填类型和风化程度确定。钻探结果确定地质强度指标从而更好的还原了工程现场的巷道围岩条件，保证了模拟结果的准确性，尤其是对地质强度指标GIS的统计值进行概率分布拟合，确定了顶板的模拟参数，从而能够更准确的计算巷道顶板下沉变形值。
步骤四.对数值模型的材料参数进行赋值，材料参数包括弹性模量、泊松比、内聚力、内摩擦角、抗拉强度、残余内聚力和塑性应变量；
在本步骤中弹性模量E和地质强度指标GIS之间的关系为：
其中σm为岩石单轴抗压强度。
步骤五.设定边界条件，模型开挖计算，边界条件包括施加重力模拟初始应力，设置应力边界条件和位移边界条件，所述开挖计算使用有限差分计算方法；
在本步骤中模型开挖使用有限差分法进行迭代计算，当模型的最大不平衡力当低于1e-5时，说明模型达到平衡状态，模拟运算结束。
步骤六.提取监测数据，进行巷道稳定性评估和二次支护设计。
在本步骤中的二次支护设计，在巷道变形值达到120mm以上的区域，采取注浆加固和增加锚杆预紧力的支护措施。根据巷道模拟结果，提取巷道变形量监测数据，从而对巷道围岩的稳定性和可靠性进行评估分析，从而进一步的指导巷道二次支护保证了巷道的安全和支护的有效性。
实施例2
本实施例结合某矿回采巷道的工程实际，对本发明的具体实施做进一步的说明，采用连续介质有限差分软件FLAC3D进行数值模拟，根据数值模拟结果指导巷道的二次支护，该方法流程如图1所示，具体实施如下。
在某矿的回采巷道中利用数值模拟分析非均匀分布裂隙巷道稳定性的方法指导二次支护的具体步骤包括：
步骤一.利用有限元模拟软件建立三维数值模型，包括数值模型的尺寸设计和网格划分。
在本步骤中的三维数值模型为长方体，如图2所示，设置模型长度为200m，模型宽度为80m，高度为62m，在三维数值模型中部开挖巷道。巷道为长方形，巷道长度为200m，巷道两侧到边界的距离取35m巷道中部100m的区域划分为顶板监测区域，巷道顶底板岩层根据图3所示的巷道顶底板岩性特征和厚度柱状图设置。
步骤二.非均匀岩体概率分布参数设定，包括确定概率分布函数，模拟岩体的非均质性。
在本步骤中概率分布函数取韦布尔概率分布函数，韦布尔概率分布能够模拟一个区间内的随机分布岩体非均质参数，不同均质性系数m的韦布尔分布曲线如图4所示，韦布尔分布函数如下：
其中u是单个元素的力学参数值，u0是尺度参数，随着尺度参数的改变函数曲线放大或缩小，m为材料的均质性系数，随着均质性系数的改变曲线的形状随之改变。
其中使用蒙特卡洛法模拟岩体的非均质性，步骤包括：A.在[μmin，μmax]之间随机取一个值μi，计算μi的概率密度函数f(μi)；B.在f(μmin)与f(μmax)之间随机产生一个概率密度函数f(μ)，当f(μ)<f(μi)时，导出μi，C.继续进行迭代计算，最终得到μi的数据组，μi的数据组服从韦布尔概率分布规律，将该方法通过FISH语言编写，在FLAC3D中运行该程序。利用韦布尔概率分布函数模拟岩体的非均质性，提高了模拟的合理性和可靠性。
步骤三.测量并计算巷道围岩的地质强度指标，包括钻探取芯获取岩石质量指标和节理条件，地质强度指标的概率分布拟合。在回采巷道内选取52个测点，对测点的地质强度指标GIS进行统计处理，统计处理的步骤包括：A.将地质强度指标GIS统计值制作概率分布直方图，如图5所示；B.利用origin软件处理数据，进行韦布尔Weibull概率分布参数拟合，得到GSIave＝40.72，均质性系数m＝1.96；C.根据拟合结果取GSImax＝10与GSImin＝80。
在本步骤中的地质强度指标GIS与岩石质量指标RQD和节理条件JCond89的关系为：GSI＝1.5JCond89+RQD/2。岩石质量指标RQD取钻探获取的岩芯完整度来表示，具体为长度大于10cm的岩芯占钻孔总长度的比值。其中节理条件JCond89的取值根据巷道围岩裂隙长度、宽度、粗糙度、充填类型和风化程度确定。钻探结果确定地质强度指标从而更好的还原了工程现场的巷道围岩条件，保证了模拟结果的准确性，尤其是对地质强度指标GIS的统计值进行概率分布拟合，确定了顶板的模拟参数，从而能够更准确的计算巷道顶板下沉变形值。
步骤四.对数值模型的材料参数进行赋值，材料参数包括弹性模量、泊松比、内聚力、内摩擦角、抗拉强度、残余内聚力和塑性应变量，各层的材料参数如表1所示。
表1.材料参数取值表
在本步骤中弹性模量E和地质强度指标GIS之间的关系为：
其中σm为岩石单轴抗压强度，根据该公式计算顶板岩体的弹性模量。
利用FLAC3D中的内置Cable结构单元模拟锚杆和锚索支护结构，设置的支护结构力学参数如表2所示。
表2.支护参数取值表
步骤五.设定边界条件，模型开挖计算，边界条件包括施加重力模拟初始应力，设置应力边界条件和位移边界条件，开挖计算使用有限差分计算方法。在模型顶部施加15MPa的垂直应力，模拟巷道实际埋深650m的实际情况，在水平方向上分别施加水平应力12MPa和18MPa，模型的四周和底部设置位移边界条件。采用显式有限差分计算法进行迭代计算，模型的最大不平衡力低于1e-5时，说明模型达到平衡状态，停止计算，初始应力场模拟完成。
在本步骤中模型开挖使用有限差分法进行迭代计算，具体是：将开挖区域定义为空模型null，巷道的开挖以掘进10m为一个循环，支护紧随巷道开挖，设置支护体并赋予力学参数，运算至模型平衡后，进行下一个10m的开挖循环。当模型的最大不平衡力当低于1e-5时，说明模型达到平衡状态，模拟运算结束。
步骤六.提取监测数据，进行巷道稳定性评估和二次支护设计。具体是利用hist命令对巷道顶板的变形进行监测，当巷道开挖结束后提取巷道变形量的监测数据，输出巷道顶板下沉结果如图6所示。
在本步骤中的二次支护设计，在巷道变形值达到120mm以上的区域，采取注浆加固和增加锚杆预紧力的支护措施。根据巷道模拟结果，提取巷道变形量监测数据，从而对巷道围岩的稳定性和可靠性进行评估分析，尤其是巷道变形量差异明显的位置，对巷道变形量大的位置进行重点加强支护，从而提升支护效率和支护效果，进一步的指导巷道二次支护保证了巷道的安全和支护的有效性。
当然，上述说明并非是对本发明的限制，本发明也并不仅限于上述举例，本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换，也应属于本发明的保护范围。
Claims (9)
1.一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法，其特征在于，步骤包括：
步骤一.利用有限元模拟软件建立三维数值模型，包括数值模型的尺寸设计和网格划分；
步骤二.非均匀岩体概率分布参数设定，包括确定概率分布函数，模拟岩体的非均质性；
步骤三.测量并计算巷道围岩的地质强度指标，包括钻探取芯获取岩石质量指标和节理条件，地质强度指标的概率分布拟合；
步骤四.对数值模型的材料参数进行赋值，所述材料参数包括弹性模量、泊松比、内聚力、内摩擦角、抗拉强度、残余内聚力和塑性应变量；
步骤五.设定边界条件，模型开挖计算，所述边界条件包括施加重力模拟初始应力，设置应力边界条件和位移边界条件，所述开挖计算使用有限差分计算方法；
步骤六.提取监测数据，进行巷道稳定性评估和二次支护设计。
2.根据权利要求1所述的一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法，其特征在于，所述步骤一中的三维数值模型为长方体，在三维数值模型中部开挖巷道，巷道中部划分为顶板监测区域。
3.根据权利要求1所述的一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法，其特征在于，所述步骤二中概率分布函数取韦布尔概率分布函数：
其中u是单个元素的力学参数值，u0是尺度参数，m为材料的均质性系数。
4.根据权利要求3所述的一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法，其特征在于，所述步骤二中使用蒙特卡洛法模拟岩体的非均质性，包括：A.在[μmin，μmax]之间随机取一个值μi，计算μi的概率密度函数f(μi)；B.在f(μmin)与f(μmax)之间随机产生一个概率密度函数f(μ)，当f(μ)<f(μi)时，导出μi，C.继续进行迭代计算，最终得到μi的数据组，μi的数据组服从韦布尔概率分布规律。
5.根据权利要求1所述的一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法，其特征在于，所述步骤三中的地质强度指标GIS与岩石质量指标RQD和节理条件JCond89的关系为：GSI＝1.5JCond89+RQD/2。
6.根据权利要求5所述的一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法，其特征在于，所述岩石质量指标RQD取钻探获取的岩芯完整度来表示，具体为长度大于10cm的岩芯占钻孔总长度的比值；所述节理条件JCond89的取值根据巷道围岩裂隙长度、宽度、粗糙度、充填类型和风化程度确定。
7.根据权利要求1所述的一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法，其特征在于，所述步骤四中弹性模量E和地质强度指标GIS之间的关系为：
其中σm为岩石单轴抗压强度。
8.根据权利要求1所述的一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法，其特征在于，所述步骤五中模型开挖使用有限差分法进行迭代计算，当模型的最大不平衡力当低于1e-5时模拟运算结束。
9.根据权利要求1所述的一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法，其特征在于，所述步骤六中二次支护设计包括，在巷道变形值达到120mm以上的区域，采取注浆加固和增加锚杆预紧力的支护措施。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811201014.4A CN109492262A (zh) | 2018-10-16 | 2018-10-16 | 一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811201014.4A CN109492262A (zh) | 2018-10-16 | 2018-10-16 | 一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109492262A true CN109492262A (zh) | 2019-03-19 |
Family
ID=65690269
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811201014.4A Pending CN109492262A (zh) | 2018-10-16 | 2018-10-16 | 一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109492262A (zh) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110390152A (zh) * | 2019-07-15 | 2019-10-29 | 中国矿业大学 | 一种模拟巷道围岩裂隙演化的离散元方法 |
| CN110847969A (zh) * | 2019-10-29 | 2020-02-28 | 中国电建集团华东勘测设计研究院有限公司 | 一种适用于岩体条件地下洞室群变形分级预警指标的确定方法 |
| CN111222094A (zh) * | 2019-10-24 | 2020-06-02 | 中国地质大学（武汉） | Fouché裂隙产状概率分布计算方法运用后残留误差的评价技术 |
| CN111666699A (zh) * | 2020-04-30 | 2020-09-15 | 山东大学 | 基于rev全区域覆盖的岩体工程跨尺度模拟计算方法 |
| CN111695790A (zh) * | 2020-05-26 | 2020-09-22 | 中钢集团山东矿业有限公司 | 一种保安矿柱开采方法 |
| CN111982377A (zh) * | 2020-09-01 | 2020-11-24 | 中国地质科学院地质力学研究所 | 一种基于二次套芯的深孔差应力评估方法 |
| CN112257326A (zh) * | 2020-10-27 | 2021-01-22 | 中国矿业大学(北京) | 一种围岩应力状态的数值模拟方法及系统 |
| CN113343423A (zh) * | 2021-05-08 | 2021-09-03 | 武汉大学 | 基于强度空间变异性的随机裂隙网络生成方法 |
| CN113435087A (zh) * | 2021-06-24 | 2021-09-24 | 中铁二院工程集团有限责任公司 | 一种分析洞穴围岩局部稳定性的方法 |
| CN116882012A (zh) * | 2023-07-11 | 2023-10-13 | 安徽建筑大学 | 深部巷道底鼓影响因素及特征分析方法 |
| CN117290928A (zh) * | 2023-09-25 | 2023-12-26 | 西南交通大学 | 基于随钻参数的隧道围岩力学参数的反演方法及装置 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105550441A (zh) * | 2015-12-12 | 2016-05-04 | 山东科技大学 | 一种基于连续介质的工程岩体破裂劣化数值模拟方法 |
| CN105629308A (zh) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | 相控非均质力学参数地应力方法 |
| CN106326524A (zh) * | 2016-07-28 | 2017-01-11 | 张树森 | 一种非均质地层应力场数值模拟方法 |
| CN108170959A (zh) * | 2017-12-28 | 2018-06-15 | 天地科技股份有限公司 | 基于离散元的岩体力学响应数值分析方法及装置 |
-
2018
- 2018-10-16 CN CN201811201014.4A patent/CN109492262A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105629308A (zh) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | 相控非均质力学参数地应力方法 |
| CN105550441A (zh) * | 2015-12-12 | 2016-05-04 | 山东科技大学 | 一种基于连续介质的工程岩体破裂劣化数值模拟方法 |
| CN106326524A (zh) * | 2016-07-28 | 2017-01-11 | 张树森 | 一种非均质地层应力场数值模拟方法 |
| CN108170959A (zh) * | 2017-12-28 | 2018-06-15 | 天地科技股份有限公司 | 基于离散元的岩体力学响应数值分析方法及装置 |
Non-Patent Citations (4)
| Title |
|---|
| 李朋朋: "深部软岩巷道围岩稳定性分析与控制技术研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
| 罗荣: "基于材料组分的非均质岩土体数值建模方法及试验研究", 《中国优秀博士学位论文全文数据库信息科技辑》 * |
| 肖术: "基于可靠度理论的露天矿节理岩体边坡稳定性及设计优化研究", 《中国优秀博士学位论文全文数据库信息科技辑》 * |
| 郑晨: "近接隧道深基坑开挖及节理破碎岩体对隧道稳定性的影响研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Cited By (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11209415B2 (en) | 2019-07-15 | 2021-12-28 | China University Of Mining And Technology | Discrete element method for modelling a fracture evolution of a roadway surrounding rock |
| CN110390152A (zh) * | 2019-07-15 | 2019-10-29 | 中国矿业大学 | 一种模拟巷道围岩裂隙演化的离散元方法 |
| WO2021007959A1 (zh) * | 2019-07-15 | 2021-01-21 | 中国矿业大学 | 一种模拟巷道围岩裂隙演化的离散元方法 |
| AU2019443098B2 (en) * | 2019-07-15 | 2022-09-22 | Boxia Mining Technology Ltd. | Discrete element method for modelling fracture evolution of roadway surrounding rock |
| CN111222094A (zh) * | 2019-10-24 | 2020-06-02 | 中国地质大学（武汉） | Fouché裂隙产状概率分布计算方法运用后残留误差的评价技术 |
| CN111222094B (zh) * | 2019-10-24 | 2021-08-27 | 中国地质大学（武汉） | Fouché裂隙产状概率分布计算方法运用后残留误差的评价方法 |
| CN110847969A (zh) * | 2019-10-29 | 2020-02-28 | 中国电建集团华东勘测设计研究院有限公司 | 一种适用于岩体条件地下洞室群变形分级预警指标的确定方法 |
| CN110847969B (zh) * | 2019-10-29 | 2021-06-08 | 中国电建集团华东勘测设计研究院有限公司 | 一种适用于岩体条件地下洞室群变形分级预警指标的确定方法 |
| CN111666699A (zh) * | 2020-04-30 | 2020-09-15 | 山东大学 | 基于rev全区域覆盖的岩体工程跨尺度模拟计算方法 |
| CN111666699B (zh) * | 2020-04-30 | 2023-06-02 | 山东大学 | 基于rev全区域覆盖的岩体工程跨尺度模拟计算方法 |
| CN111695790A (zh) * | 2020-05-26 | 2020-09-22 | 中钢集团山东矿业有限公司 | 一种保安矿柱开采方法 |
| CN111982377A (zh) * | 2020-09-01 | 2020-11-24 | 中国地质科学院地质力学研究所 | 一种基于二次套芯的深孔差应力评估方法 |
| CN111982377B (zh) * | 2020-09-01 | 2021-07-23 | 中国地质科学院地质力学研究所 | 一种基于二次套芯的深孔差应力评估方法 |
| CN112257326B (zh) * | 2020-10-27 | 2023-05-30 | 中国矿业大学(北京) | 一种围岩应力状态的数值模拟方法及系统 |
| CN112257326A (zh) * | 2020-10-27 | 2021-01-22 | 中国矿业大学(北京) | 一种围岩应力状态的数值模拟方法及系统 |
| CN113343423A (zh) * | 2021-05-08 | 2021-09-03 | 武汉大学 | 基于强度空间变异性的随机裂隙网络生成方法 |
| CN113343423B (zh) * | 2021-05-08 | 2022-04-29 | 武汉大学 | 基于强度空间变异性的随机裂隙网络生成方法 |
| CN113435087A (zh) * | 2021-06-24 | 2021-09-24 | 中铁二院工程集团有限责任公司 | 一种分析洞穴围岩局部稳定性的方法 |
| CN113435087B (zh) * | 2021-06-24 | 2022-04-22 | 中铁二院工程集团有限责任公司 | 一种分析洞穴围岩局部稳定性的方法 |
| CN116882012A (zh) * | 2023-07-11 | 2023-10-13 | 安徽建筑大学 | 深部巷道底鼓影响因素及特征分析方法 |
| CN116882012B (zh) * | 2023-07-11 | 2024-03-29 | 安徽建筑大学 | 深部巷道底鼓影响因素及特征分析方法 |
| CN117290928A (zh) * | 2023-09-25 | 2023-12-26 | 西南交通大学 | 基于随钻参数的隧道围岩力学参数的反演方法及装置 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109492262A (zh) | 一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法 | |
| CN109271738A (zh) | 一种用于获取巷道围岩Weibull分布参数的数值反演方法 | |
| Jiang et al. | Observe the temporal evolution of deep tunnel's 3D deformation by 3D laser scanning in the Jinchuan No. 2 Mine | |
| Zhang et al. | Estimation of in situ stress along deep tunnels buried in complex geological conditions | |
| He et al. | Discontinuous deformation analysis of super section tunnel surrounding rock stability based on joint distribution simulation | |
| CN103866769B (zh) | 基于3d地质模型及实时监控大坝灌浆工程分析控制方法 | |
| CN103293560B (zh) | 一种采动三向应力场的测试方法 | |
| CN110135030B (zh) | 一种采空区地表沉降的预测方法 | |
| CN103953392B (zh) | 深埋隧道断面上岩爆风险位置判别方法 | |
| Ding et al. | Stability of large underground caverns excavated in layered rock masses with steep dip angles: a case study | |
| Xu et al. | Estimation of the three-dimensional in situ stress field around a large deep underground cavern group near a valley | |
| Wu et al. | Advances in statistical mechanics of rock masses and its engineering applications | |
| Li et al. | Microseismic monitoring and numerical simulation on the stability of high-steep rock slopes in hydropower engineering | |
| CN110162882A (zh) | 一种煤矿巷道穿越断层区围岩精准控制方法 | |
| CN116703244B (zh) | 采矿沉陷区治理效果与综合评价方法 | |
| Cheng et al. | Experimental study on the deformation and failure mechanism of overburden rock during coal mining using a comprehensive intelligent sensing method | |
| Zhang et al. | Tunnel stability analysis of coral reef limestone stratum in ocean engineering | |
| CN111695790A (zh) | 一种保安矿柱开采方法 | |
| CN113536414B (zh) | 基于三维建模的岩质边坡稳定性分析方法、系统及介质 | |
| Sun et al. | Calibration of water inrush channel and numerical simulation of water inrush process in coal mine roof | |
| Jia et al. | Stability analysis of shallow goaf based on field monitoring and numerical simulation: A case study at an Open-Pit Iron Mine, China | |
| Pang et al. | Analysis of main controlling factors of overburden failure in coal mining under thick coal seam geological conditions | |
| SHI et al. | Physical experimental study on excavation disturbance of TBM in deep composite strata | |
| Sun et al. | Failure mechanism and control countermeasures of surrounding rock at deep large section chamber intersection in the Wanfu Coal Mine | |
| Hao et al. | Characteristics of ground surface settlement of double-line adjacent metro construction in sandy cobble stratum: A case study of Beijing Airport Line |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190319 |
|
| RJ01 | Rejection of invention patent application after publication |