CN109271738A - 一种用于获取巷道围岩Weibull分布参数的数值反演方法 - Google Patents
一种用于获取巷道围岩Weibull分布参数的数值反演方法 Download PDFInfo
- Publication number
- CN109271738A CN109271738A CN201811201013.XA CN201811201013A CN109271738A CN 109271738 A CN109271738 A CN 109271738A CN 201811201013 A CN201811201013 A CN 201811201013A CN 109271738 A CN109271738 A CN 109271738A
- Authority
- CN
- China
- Prior art keywords
- rock
- roadway
- surrounding rock
- value
- parameter
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明提供了一种用于获取巷道围岩Weibull分布参数的数值反演方法,涉及采矿岩石力学技术领域,解决了数值模拟过程中巷道非均质分布裂隙岩体参数确定的技术问题,该方法的步骤包括:(1)巷道变形监测;(2)巷道围岩的地质强度指标的区间估计;(3)利用有限元模拟软件建立三维数值模型;(4)非均匀岩体概率分布参数设定;(5)对数值模型的材料参数进行赋值;(6)设置边界条件,模型开挖计算;(7)提取监测数据,监测巷道变形;(8)计算相关系数,校核岩体非均质参数。本方法还具有参数反演精度高,符合现场实际工程条件,可靠性高等优点。
Description
技术领域
本发明涉及采矿岩石力学技术领域,尤其是一种用于获取巷道围岩Weibull分布地质强度指标的平均值和材料的均质性系数参数的数值反演方法。
背景技术
岩土工程的隧道、矿业、水电、核电等领域,尤其是地下工程中,岩体裂隙发育具有随机性、形态多样性和空间组合复杂性等特点,从而使得工程岩体的力学参数如弹性模量等具有较高的非均质性,岩体参数的非均质性对工程稳定产生重要影响,严重时会造成井下巷道等支护工程的意外失效。研究非均匀裂隙围岩稳定性控制是岩土工程领域的重要研究课题,对保证安全、高效的施工和生产具有重要指导意义,但是常规的数值模拟方法往往将岩体视为均匀介质,忽略岩体的非均匀性对工程的影响,因此无法对非均匀裂隙发育巷道围岩控制提供准确可靠的支持。
目前,分析围岩失稳机理、研究围岩控制对策、设计支护方案等领域,岩土工程中广泛采用的数值模拟方法主要有有限元法、有限差分法和离散元法,其中基于连续介质的有限元法和有限差分法凭借其较高的计算效率和准确的计算结果被广泛应用,计算结果比较符合现场的实际情况。岩体中裂隙的产状和发育程度严重影响岩体的强度与变形参数进而控制着巷道变形规律,但是由于岩体裂隙发育的非均匀性造成不同区域巷道围岩力学参数差异性明显,围岩力学参数的非均匀分布对围岩失稳机理、控制对策研究、支护方案设计等产生重要影响,是不可忽略的。因此,基于连续介质数值模拟方法,在保证计算准确性和高效性的前提下,需要考虑围岩裂隙的非均匀分布对岩体参数的影响,在这个过程中反演岩体非均质性参数一直是该领域的难点,但是岩体的非均质性参数对提高数值模拟结果的可靠性具有重要意义,因此需要提出一种获取巷道围岩分布参数的数值反演方法。
发明内容
为解决数值模拟过程中巷道非均质分布裂隙岩体参数确定,尤其是地质强度指标和材料均质性系数参数确定的技术问题,本发明提供了一种用于获取巷道围岩Weibull分布参数的数值反演方法,具体技术方案如下。
一种用于获取巷道围岩Weibull分布参数的数值反演方法,步骤包括:
步骤一.巷道变形监测,包括布置多个巷道变形监测点,统计巷道变形值并将巷道变形值划分为多个区间;
步骤二.巷道围岩的地质强度指标的区间估计,包括钻探取芯获取岩石质量指标和节理条件,地质强度指标的概率分布拟合;
步骤三.利用有限元模拟软件建立三维数值模型,包括数值模型的尺寸设计和网格划分;
步骤四.非均匀岩体概率分布参数设定,包括确定概率分布函数,模拟岩体的非均质性;
步骤五.对数值模型的材料参数进行赋值,所述材料参数包括弹性模量、泊松比、内聚力、内摩擦角、抗拉强度、残余内聚力和塑性应变量;
步骤六.设置边界条件,模型开挖计算,所述开挖计算使用有限差分计算方法;
步骤七.提取数值模型的监测数据,包括监测巷道围岩变形量,统计巷道围岩变形值并将巷道围岩变形值划分为多个区间;
步骤八.计算相关系数,校核岩体非均质参数,包括计算数值模拟监测数据和巷道变形监测数据的相关系数,当相关系数大于0.9时,反演得出地质强度指标的平均值和材料的均质性系数;当相关系数小于0.9时,重复步骤五至步骤八,最后确定地质强度指标的平均值和材料的均质性系数。
优选的是,步骤一中选取的巷道变形监测点根据巷道长度确定测点间隔为2~4mm,将巷道变形值每间隔5~10mm划分为一个统计区间。
优选的是,步骤二中地质强度指标GIS与岩石质量指标RQD和节理条件JCond89的关系为:GSI=1.5JCond89+RQD/2。
进一步优选的是,岩石质量指标RQD取钻探获取的岩芯完整度来表示,具体为长度大于10cm的岩芯占钻孔总长度的比值;所述节理条件JCond89的取值根据巷道围岩裂隙长度、宽度、粗糙度、充填类型和风化程度确定。
进一步优选的是,在巷道围岩完整区域和巷道围岩破碎区域分别钻探获取岩芯,分别计算巷道围岩完整区域和巷道围岩破碎区域的岩石质量指标RQD值以及节理条件JCond89值,从而确定围岩完整区域的地质强度指标GISmax和围岩破碎区域的地质强度指标GISmin,得到地质强度指标GIS的区间[GISmin,GISmax]。
优选的是,步骤三中的三维数值模型为长方体,在三维数值模型中部开挖巷道,巷道中部划分为顶板监测区域。
优选的是,步骤四中概率分布函数取韦布尔概率分布函数:
其中u是单个元素的力学参数值,u0是尺度参数,m为材料的均质性系数。
使用蒙特卡洛法模拟岩体的非均质性,包括:A.在[μmin,μmax]之间随机取一个值μi,计算μi的概率密度函数f(μi);B.在f(μmin)与f(μmax)之间随机产生一个概率密度函数f(μ),当f(μ)<f(μi)时,导出μi,C.继续进行迭代计算,最终得到μi的数据组,μi的数据组服从韦布尔概率分布规律。
优选的是,步骤五中弹性模量E和地质强度指标GIS之间的关系为:
其中σm为岩石单轴抗压强度。
优选的是,步骤六中模型开挖使用有限差分法进行迭代计算,当模型的最大不平衡力当低于1e-5时模拟运算结束。
优选的是,步骤七中选取巷道变形值间隔划分区间和步骤一中的变形值间隔划分相同。
本发明的有益效果包括:
(1)本方法结合现场监测和数值模拟对非均质岩体的巷道裂隙分布和顶板下沉规律进行模拟,利用模拟结果和现场监测相互验证,并根据两组监测数据的相关系数确定地质强度指标的平均值和材料的均质性系数,提高了数值模拟的合理性和可靠性。
(2)根据钻探结果确定地质强度指标从而更好的还原了工程现场的巷道围岩条件,保证了模拟结果的准确性,尤其是对地质强度指标GIS的统计值进行概率分布拟合,确定了顶板的模拟参数,从而能够更准确的确定地质强度指标的取值范围。
另外本方法还具有操作简单,参数反演精度高,符合现场实际工程条件,可靠性高等优点。
附图说明
图1是利用数值模拟分析非均匀分布裂隙巷道稳定性的方法流程结构图;
图2是巷道测点布置情况示意图;
图3是顶板下沉情况统计分布示意图;
图4是三维数值模型结构示意图;
图5是巷道顶底板岩性特征和厚度柱状图;
图6是不同均质性系数的韦布尔分布曲线图;
图7是数值模拟结果巷道顶板下沉示意图;
图8是数值模拟监测数据和巷道变形监测数据对比示意图。
具体实施方式
结合图1至图8所示,本发明提供的一种用于获取巷道围岩Weibull分布参数的数值反演方法,具体实施方式如下。
实施例1
一种用于获取巷道围岩Weibull分布参数的数值反演方法具体的步骤包括:
步骤一.巷道变形监测,包括布置多个巷道变形监测点,统计巷道变形值并将巷道变形值划分为多个区间。
在本步骤中选取的巷道变形监测点根据巷道长度确定测点间隔可以取2~4mm,将巷道变形值每间隔5~10mm划分为一个统计区间,具体根据巷道的长度和巷道的实际情况进行选取,选取尽量多的测点进行测量和统计。
步骤二.巷道围岩的地质强度指标的区间估计,包括钻探取芯获取岩石质量指标和节理条件,地质强度指标的概率分布拟合。
在本步骤中地质强度指标GIS与岩石质量指标RQD和节理条件JCond89的关系为:GSI=1.5JCond89+RQD/2。
其中岩石质量指标RQD取钻探获取的岩芯完整度来表示,具体为长度大于10cm的岩芯占钻孔总长度的比值。节理条件JCond89的取值根据巷道围岩裂隙长度、宽度、粗糙度、充填类型和风化程度确定。在巷道围岩完整区域和巷道围岩破碎区域分别钻探获取岩芯,分别计算巷道围岩完整区域和巷道围岩破碎区域的岩石质量指标RQD值以及节理条件JCond89值,从而确定围岩完整区域的地质强度指标GISmax和围岩破碎区域的地质强度指标GISmin,得到地质强度指标GIS的区间[GISmin,GISmax]。根据钻探结果确定地质强度指标从而更好的还原了工程现场的巷道围岩条件,保证了模拟结果的准确性,尤其是对地质强度指标GIS的统计值进行概率分布拟合,确定了顶板的模拟参数,从而能够更准确的确定地质强度指标的取值范围。
步骤三.利用有限元模拟软件建立三维数值模型,包括数值模型的尺寸设计和网格划分。
在本步骤中的三维数值模型为长方体,在三维数值模型中部开挖巷道,巷道中部划分为顶板监测区域。
步骤四.非均匀岩体概率分布参数设定,包括确定概率分布函数,模拟岩体的非均质性。
在本步骤中概率分布函数取韦布尔(Weibull)概率分布函数:
其中u是单个元素的力学参数值,u0是尺度参数,m为材料的均质性系数。
使用蒙特卡洛法模拟岩体的非均质性,包括:A.在[μmin,μmax]之间随机取一个值μi,计算μi的概率密度函数f(μi);B.在f(μmin)与f(μmax)之间随机产生一个概率密度函数f(μ),当f(μ)<f(μi)时,导出μi,C.继续进行迭代计算,最终得到μi的数据组,μi的数据组服从韦布尔概率分布规律。利用韦布尔概率分布函数模拟岩体的非均质性,提高了模拟的合理性和可靠性,通过FISH语言实现对该方法的数值模拟应用。
步骤五.对数值模型的材料参数进行赋值,所述材料参数包括弹性模量、泊松比、内聚力、内摩擦角、抗拉强度、残余内聚力和塑性应变量。
其中弹性模量E和地质强度指标GIS之间的关系为:
其中σm为岩石单轴抗压强度。另外,其他材料参数根据实际地质情况和表格数据确定。
步骤六.设置边界条件,模型开挖计算,所述开挖计算使用有限差分计算方法。模型开挖使用有限差分法进行迭代计算,当模型的最大不平衡力当低于1e-5时模拟运算结束。
步骤七.提取数值模型的监测数据,包括监测巷道围岩变形量,统计巷道围岩变形值并将巷道围岩变形值划分为多个区间。其中选取巷道变形值间隔划分区间和步骤一中的变形值间隔划分相同。
步骤八.计算相关系数,校核岩体非均质参数,包括计算数值模拟监测数据和巷道变形监测数据的相关系数,当相关系数大于0.9时,反演得出地质强度指标的平均值和材料的均质性系数;当相关系数小于0.9时,重复步骤五至步骤八,最后确定地质强度指标的平均值和材料的均质性系数。
本方法结合现场监测和数值模拟对非均质岩体的巷道裂隙分布和顶板下沉规律进行模拟,利用模拟结果和现场监测相互验证,并根据两组监测数据的相关系数确定地质强度指标的平均值和材料的均质性系数,提高了数值模拟的合理性和可靠性。另外本方法还具有操作简单,参数反演精度高,符合现场实际工程条件,可靠性高等优点。
实施例2
本实施例结合某矿回采巷道的工程实际,对本发明的具体实施做进一步的说明,采用连续介质有限差分软件FLAC3D进行数值模拟,根据数值模拟结果和相关系数确定地质强度指标的平均值和材料的均质性系数的选择是否合理,该参数反演方法的流程如图1所示,具体实施如下。
一种用于获取巷道围岩Weibull分布参数的数值反演方法具体的步骤包括:
步骤一.巷道变形监测,包括布置多个巷道变形监测点,统计巷道变形值并将巷道变形值划分为多个区间。
在回采巷道的顶板设置多个测点,对巷道变形规律进行统计分析,每一个巷道断面设置5个测点,其中测点位置如图2所示,其中同组测点的间隔为0.8m,沿巷道走向方向的测点间隔2m,共测量100m长度的巷道变形量,得到累计测量数据255个。对巷道变形值进行统计处理,将巷道变形值分成若干个统计区间并计算每个区间的比例分数,从60mm~140mm区间没每隔10mm划分为一个统计区间,统计情况如图3所示。
步骤二.巷道围岩的地质强度指标的区间估计,包括钻探取芯获取岩石质量指标和节理条件,地质强度指标的概率分布拟合。
在本步骤中地质强度指标GIS与岩石质量指标RQD和节理条件JCond89的关系为:GSI=1.5JCond89+RQD/2。利用origin软件处理现场巷道变形监测统计数据,进行韦布尔(Weibull)概率分布参数拟合得到地质强度指标的取值。
其中岩石质量指标RQD取钻探获取的岩芯完整度来表示,具体为长度大于10cm的岩芯占钻孔总长度的比值。节理条件JCond89的取值根据巷道围岩裂隙长度、宽度、粗糙度、充填类型和风化程度的实际情况确定。在巷道围岩完整区域和巷道围岩破碎区域分别钻探获取岩芯,分别计算巷道围岩完整区域和巷道围岩破碎区域的岩石质量指标RQD值以及节理条件JCond89值,从而确定围岩完整区域的地质强度指标GISmax和围岩破碎区域的地质强度指标GISmin,得到地质强度指标GIS的区间为[10,80]。根据钻探结果确定地质强度指标从而更好的还原了工程现场的巷道围岩条件,保证了模拟结果的准确性,尤其是对地质强度指标GIS的统计值进行概率分布拟合,确定了顶板的模拟参数,从而能够更准确的确定地质强度指标的取值范围。
步骤三.利用有限元模拟软件建立三维数值模型,包括数值模型的尺寸设计和网格划分。
设置模型长度为200m,模型宽度为80m,高度为62m,在三维数值模型中部开挖巷道。巷道为长方形,巷道长度为200m,巷道两侧到边界的距离取35m巷道中部100m的区域划分为顶板监测区域,巷道顶底板岩层设置根据图4所示的巷道顶底板岩性特征和厚度柱状图设置,三维数值模型结构如图5所示。
步骤四.非均匀岩体概率分布参数设定,包括确定概率分布函数,模拟岩体的非均质性。
在本步骤中概率分布函数取韦布尔概率分布函数:
其中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中运行该程序。利用韦布尔概率分布函数模拟岩体的非均质性,提高了模拟和参数反演的合理性和可靠性。
步骤五.对数值模型的材料参数进行赋值,材料参数包括弹性模量、泊松比、内聚力、内摩擦角、抗拉强度、残余内聚力和塑性应变量,数值模型中各层的材料参数如表1所示。
表1.材料参数取值表
在本步骤中弹性模量E和地质强度指标GIS之间的关系为:
其中σm为岩石单轴抗压强度,根据该公式计算顶板岩体的弹性模量。
利用FLAC3D中的内置Cable结构单元模拟锚杆和锚索支护结构,设置的支护结构力学参数如表2所示。
表2.支护参数取值表
另外,其他材料参数根据实际地质情况和表格数据确定。
步骤六.设置边界条件,模型开挖计算,开挖计算使用有限差分计算方法。
其中,设定边界条件,模型开挖计算,边界条件包括施加重力模拟初始应力,设置应力边界条件和位移边界条件,开挖计算使用有限差分计算方法。在模型顶部施加15MPa的垂直应力,模拟巷道实际埋深650m的实际情况,在水平方向上分别施加水平应力12MPa和18MPa,模型的四周和底部设置位移边界条件。采用显式有限差分计算法进行迭代计算,模型的最大不平衡力低于1e-5时,说明模型达到平衡状态,停止计算,初始应力场模拟完成。
在本步骤中模型开挖使用有限差分法进行迭代计算,具体是:将开挖区域定义为空模型null,巷道的开挖以掘进10m为一个循环,支护紧随巷道开挖,设置支护体并赋予力学参数,运算至模型平衡后,进行下一个10m的开挖循环。当模型的最大不平衡力当低于1e-5时,说明模型达到平衡状态,模拟运算结束。
步骤七.提取数值模型的监测数据,包括监测巷道围岩变形量,统计巷道围岩变形值并将巷道围岩变形值划分为多个区间。其中选取巷道变形值间隔划分区间和步骤一中的变形值间隔划分相同,从60mm~140mm区间没每隔10mm划分为一个统计区间。
步骤八.计算相关系数,校核岩体非均质参数,包括计算数值模拟监测数据和巷道变形监测数据的相关系数,当相关系数大于0.95时,反演得出地质强度指标的平均值和材料的均质性系数;当相关系数小于0.95时,重复步骤五至步骤八,最后确定地质强度指标的平均值和材料的均质性系数。
根据相关系数的大小评价数值模拟结果的可靠程度,相关系数的计算公式为:
其中,r为数值模拟监测数据和巷道变形监测数据的相关系数,n为将巷道变形值的分布范围所划分的区间个数,xdfi为现场巷道变形顶板下沉值各区间的比例分数,现场巷道变形监测顶板下沉值各区间比例分数的平均值,ydni为数值模拟顶板下沉值各区间的比例分数,为数值模拟顶板下沉值各区间比例分数的平均值。
r的取值在-1与+1之间,若r>0,表明两个变量是正相关,即一个变量的值越大,另一个变量的值也会越大;若r<0,表明两个变量是负相关,即一个变量的值越大另一个变量的值反而会越小,r的绝对值越大表明相关性越强。统计数值模拟监测数据和巷道变形监测数据中不同巷道变形值所占百分比,并计算两组数据的相关系数,量化数值模拟计算结果与现场实测结果的相似程度。
根据所需要数值模拟结果的精度要求,人为设定相关系数的取值,当数值模拟结果与现场实测结果计算的相关系数大于设定值时,反演得出GSIave和均质性参数m,当数值模拟结果与现场实测结果计算的相关系数小于设定值时,返回第六步改变GSIave和均质性系数m重新进行数值模拟解算,直到计算的相关性系数大于设定值。
本方法结合现场监测和数值模拟对非均质岩体的巷道裂隙分布和顶板下沉规律进行模拟,利用模拟结果和现场监测相互验证,并根据两组监测数据的相关系数确定地质强度指标的平均值和材料的均质性系数,提高了数值模拟的合理性和可靠性。另外本方法还具有操作简单,参数反演精度高,符合现场实际工程条件,可靠性高等优点。通过不断试错计算,当GSIave取45,均质性系数m取2时,数值模拟监测数据和巷道变形监测数据的相关系数达到0.98,反演得出GSIave为45,均质性系数m为2,同时也证明了利用该方法进行数值模拟与参数反演的步骤精度较高。
当然,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。
Claims (10)
1.一种用于获取巷道围岩Weibull分布参数的数值反演方法,其特征在于,步骤包括:
步骤一.巷道变形监测,包括布置多个巷道变形监测点,统计巷道变形值并将巷道变形值划分为多个区间;
步骤二.巷道围岩的地质强度指标的区间估计,包括钻探取芯获取岩石质量指标和节理条件,地质强度指标的概率分布拟合;
步骤三.利用有限元模拟软件建立三维数值模型,包括数值模型的尺寸设计和网格划分;
步骤四.非均匀岩体概率分布参数设定,包括确定概率分布函数,模拟岩体的非均质性;
步骤五.对数值模型的材料参数进行赋值,所述材料参数包括弹性模量、泊松比、内聚力、内摩擦角、抗拉强度、残余内聚力和塑性应变量;
步骤六.设置边界条件,模型开挖计算,所述开挖计算使用有限差分计算方法;
步骤七.提取数值模型的监测数据,包括监测巷道围岩变形量,统计巷道围岩变形值并将巷道围岩变形值划分为多个区间;
步骤八.计算相关系数,校核岩体非均质参数,包括计算数值模拟监测数据和巷道变形监测数据的相关系数,当相关系数大于0.9时,反演得出地质强度指标的平均值和材料的均质性系数;当相关系数小于0.9时,重复步骤五至步骤八,最后确定地质强度指标的平均值和材料的均质性系数。
2.根据权利要求1所述的一种用于获取巷道围岩Weibull分布参数的数值反演方法,其特征在于,所述步骤一中选取的巷道变形监测点根据巷道长度确定测点间隔为2~4mm,将巷道变形值每间隔5~10mm划分为一个统计区间。
3.根据权利要求1所述的一种用于获取巷道围岩Weibull分布参数的数值反演方法,其特征在于,所述步骤二中地质强度指标GIS与岩石质量指标RQD和节理条件JCond89的关系为:GSI=1.5JCond89+RQD/2。
4.根据权利要求3所述的一种用于获取巷道围岩Weibull分布参数的数值反演方法,其特征在于,所述岩石质量指标RQD取钻探获取的岩芯完整度来表示,具体为长度大于10cm的岩芯占钻孔总长度的比值;所述节理条件JCond89的取值根据巷道围岩裂隙长度、宽度、粗糙度、充填类型和风化程度确定。
5.根据权利要求3或4任一项所述的一种用于获取巷道围岩Weibull分布参数的数值反演方法,其特征在于,在巷道围岩完整区域和巷道围岩破碎区域分别钻探获取岩芯,分别计算巷道围岩完整区域和巷道围岩破碎区域的岩石质量指标RQD值以及节理条件JCond89值,从而确定围岩完整区域的地质强度指标GISmax和围岩破碎区域的地质强度指标GISmin,得到地质强度指标GIS的区间[GISmin,GISmax]。
6.根据权利要求1所述的一种用于获取巷道围岩Weibull分布参数的数值反演方法,其特征在于,所述步骤三中的三维数值模型为长方体,在三维数值模型中部开挖巷道,巷道中部划分为顶板监测区域。
7.根据权利要求1所述的一种用于获取巷道围岩Weibull分布参数的数值反演方法,其特征在于,所述步骤四中概率分布函数取韦布尔概率分布函数:
其中u是单个元素的力学参数值,u0是尺度参数,m为材料的均质性系数。
使用蒙特卡洛法模拟岩体的非均质性,包括:A.在[μmin,μmax]之间随机取一个值μi,计算μi的概率密度函数f(μi);B.在f(μmin)与f(μmax)之间随机产生一个概率密度函数f(μ),当f(μ)<f(μi)时,导出μi,C.继续进行迭代计算,最终得到μi的数据组,μi的数据组服从韦布尔概率分布规律。
8.根据权利要求1所述的一种用于获取巷道围岩Weibull分布参数的数值反演方法,其特征在于,所述步骤五中弹性模量E和地质强度指标GIS之间的关系为:
其中σm为岩石单轴抗压强度。
9.根据权利要求1所述的一种用于获取巷道围岩Weibull分布参数的数值反演方法,其特征在于,所述步骤六中模型开挖使用有限差分法进行迭代计算,当模型的最大不平衡力当低于1e-5时模拟运算结束。
10.根据权利要求1所述的一种用于获取巷道围岩Weibull分布参数的数值反演方法,其特征在于,所述步骤七中选取巷道变形值间隔划分区间和步骤一中的变形值间隔划分相同。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811201013.XA CN109271738B (zh) | 2018-10-16 | 2018-10-16 | 一种用于获取巷道围岩Weibull分布参数的数值反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811201013.XA CN109271738B (zh) | 2018-10-16 | 2018-10-16 | 一种用于获取巷道围岩Weibull分布参数的数值反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109271738A true CN109271738A (zh) | 2019-01-25 |
CN109271738B CN109271738B (zh) | 2023-05-12 |
Family
ID=65196715
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811201013.XA Active CN109271738B (zh) | 2018-10-16 | 2018-10-16 | 一种用于获取巷道围岩Weibull分布参数的数值反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109271738B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109854252A (zh) * | 2019-02-01 | 2019-06-07 | 王�琦 | 地下工程分阶段完整控制方法 |
CN110390152A (zh) * | 2019-07-15 | 2019-10-29 | 中国矿业大学 | 一种模拟巷道围岩裂隙演化的离散元方法 |
CN110987749A (zh) * | 2019-10-30 | 2020-04-10 | 武汉大学 | 多相复合材料的等效渗透系数细观尺度研究方法 |
CN111027236A (zh) * | 2019-10-30 | 2020-04-17 | 武汉大学 | 多相复合材料的等效导热系数的细观尺度研究方法 |
CN111400926A (zh) * | 2020-03-27 | 2020-07-10 | 安徽理工大学 | 一种基于flac3d数值软件的巷道围岩热物理参数反演方法 |
CN111946397A (zh) * | 2020-08-17 | 2020-11-17 | 南京工业大学 | 一种现场评价掘进机掌子面岩土体完整性的快速方法 |
CN112632656A (zh) * | 2020-11-23 | 2021-04-09 | 中交南京交通工程管理有限公司 | 一种基于围岩空间随机特性的隧道位移分区计算方法 |
CN112818442A (zh) * | 2021-01-06 | 2021-05-18 | 河海大学 | 基于围岩空间变异性的隧道轴向差异变形设计值确定方法 |
CN112833807A (zh) * | 2021-01-04 | 2021-05-25 | 清华大学 | 一种适用于双护盾tbm的围岩变形监测方法及预测方法 |
CN113284241A (zh) * | 2021-03-11 | 2021-08-20 | 广东锡源爆破科技股份有限公司 | 一种堆形为Weibull分布的爆堆表面岩块三维模拟方法 |
CN114722478A (zh) * | 2022-04-29 | 2022-07-08 | 山西省交通规划勘察设计院有限公司 | 一种碎屑岩接触模型参数获取方法及边坡稳定性分析方法 |
CN115600398A (zh) * | 2022-10-10 | 2023-01-13 | 长安大学(Cn) | 一种基于蒙特卡洛模拟的大型洞室岩体参数概率估计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100191511A1 (en) * | 2007-08-24 | 2010-07-29 | Sheng-Yuan Hsu | Method For Multi-Scale Geomechanical Model Analysis By Computer Simulation |
CN105550441A (zh) * | 2015-12-12 | 2016-05-04 | 山东科技大学 | 一种基于连续介质的工程岩体破裂劣化数值模拟方法 |
CN105631102A (zh) * | 2015-12-24 | 2016-06-01 | 河南理工大学 | 一种深部高应力巷道钻孔卸压参数的数值模拟确定方法 |
-
2018
- 2018-10-16 CN CN201811201013.XA patent/CN109271738B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100191511A1 (en) * | 2007-08-24 | 2010-07-29 | Sheng-Yuan Hsu | Method For Multi-Scale Geomechanical Model Analysis By Computer Simulation |
CN105550441A (zh) * | 2015-12-12 | 2016-05-04 | 山东科技大学 | 一种基于连续介质的工程岩体破裂劣化数值模拟方法 |
CN105631102A (zh) * | 2015-12-24 | 2016-06-01 | 河南理工大学 | 一种深部高应力巷道钻孔卸压参数的数值模拟确定方法 |
Non-Patent Citations (3)
Title |
---|
田茂霖: "HOEK-BROWN强度准则的参数取值研究及其工程应用", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
罗荣: "基于材料组分的非均质岩土体数值建模方法及试验研究", 《中国优秀博士学位论文全文数据库信息科技辑》 * |
肖术: "基于可靠度理论的露天矿节理岩体边坡稳定性及设计优化研究", 《中国优秀博士学位论文全文数据库信息科技辑》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109854252A (zh) * | 2019-02-01 | 2019-06-07 | 王�琦 | 地下工程分阶段完整控制方法 |
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 | 中国矿业大学 | 一种模拟巷道围岩裂隙演化的离散元方法 |
RU2746748C1 (ru) * | 2019-07-15 | 2021-04-20 | Китайский Университет Горного Дела И Технологии | Метод дискретных элементов для моделирования развития разлома в породе, окружающей штрек |
AU2019443098B2 (en) * | 2019-07-15 | 2022-09-22 | Boxia Mining Technology Ltd. | Discrete element method for modelling fracture evolution of roadway surrounding rock |
CN110987749A (zh) * | 2019-10-30 | 2020-04-10 | 武汉大学 | 多相复合材料的等效渗透系数细观尺度研究方法 |
CN111027236A (zh) * | 2019-10-30 | 2020-04-17 | 武汉大学 | 多相复合材料的等效导热系数的细观尺度研究方法 |
CN111400926A (zh) * | 2020-03-27 | 2020-07-10 | 安徽理工大学 | 一种基于flac3d数值软件的巷道围岩热物理参数反演方法 |
CN111400926B (zh) * | 2020-03-27 | 2023-10-10 | 安徽理工大学 | 一种基于flac3d数值软件的巷道围岩热物理参数反演方法 |
CN111946397A (zh) * | 2020-08-17 | 2020-11-17 | 南京工业大学 | 一种现场评价掘进机掌子面岩土体完整性的快速方法 |
CN111946397B (zh) * | 2020-08-17 | 2021-04-27 | 南京工业大学 | 一种现场评价掘进机掌子面岩土体完整性的快速方法 |
CN112632656A (zh) * | 2020-11-23 | 2021-04-09 | 中交南京交通工程管理有限公司 | 一种基于围岩空间随机特性的隧道位移分区计算方法 |
CN112632656B (zh) * | 2020-11-23 | 2023-12-05 | 中交隧桥(南京)技术有限公司 | 一种基于围岩空间随机特性的隧道位移分区计算方法 |
CN112833807A (zh) * | 2021-01-04 | 2021-05-25 | 清华大学 | 一种适用于双护盾tbm的围岩变形监测方法及预测方法 |
CN112818442A (zh) * | 2021-01-06 | 2021-05-18 | 河海大学 | 基于围岩空间变异性的隧道轴向差异变形设计值确定方法 |
CN112818442B (zh) * | 2021-01-06 | 2022-10-14 | 河海大学 | 基于围岩空间变异性的隧道轴向差异变形设计值确定方法 |
CN113284241A (zh) * | 2021-03-11 | 2021-08-20 | 广东锡源爆破科技股份有限公司 | 一种堆形为Weibull分布的爆堆表面岩块三维模拟方法 |
CN114722478A (zh) * | 2022-04-29 | 2022-07-08 | 山西省交通规划勘察设计院有限公司 | 一种碎屑岩接触模型参数获取方法及边坡稳定性分析方法 |
CN114722478B (zh) * | 2022-04-29 | 2024-01-05 | 山西省交通规划勘察设计院有限公司 | 一种碎屑岩接触模型参数获取方法及边坡稳定性分析方法 |
CN115600398A (zh) * | 2022-10-10 | 2023-01-13 | 长安大学(Cn) | 一种基于蒙特卡洛模拟的大型洞室岩体参数概率估计方法 |
CN115600398B (zh) * | 2022-10-10 | 2023-06-16 | 长安大学 | 一种基于蒙特卡洛模拟的大型洞室岩体参数概率估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109271738B (zh) | 2023-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109271738A (zh) | 一种用于获取巷道围岩Weibull分布参数的数值反演方法 | |
CN109492262A (zh) | 一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法 | |
Jiang et al. | Observe the temporal evolution of deep tunnel's 3D deformation by 3D laser scanning in the Jinchuan No. 2 Mine | |
Liu et al. | Experimental investigation of the effects of heterogeneity and geostress difference on the 3D growth and distribution of hydrofracturing cracks in unconventional reservoir rocks | |
Dai et al. | Deformation forecasting and stability analysis of large-scale underground powerhouse caverns from microseismic monitoring | |
Li et al. | In situ monitoring of rockburst nucleation and evolution in the deeply buried tunnels of Jinping II hydropower station | |
Li et al. | Evolution of the broken rock zone in the mixed ground tunnel based on the DSCM | |
Hu et al. | A DEM-based study of the disturbance in dry sandy ground caused by EPB shield tunneling | |
CN110135030B (zh) | 一种采空区地表沉降的预测方法 | |
CN103293560B (zh) | 一种采动三向应力场的测试方法 | |
Ding et al. | Stability of large underground caverns excavated in layered rock masses with steep dip angles: a case study | |
CN110162882A (zh) | 一种煤矿巷道穿越断层区围岩精准控制方法 | |
Wu et al. | Application of equivalent continuum and discontinuum stress analyses in three-dimensions to investigate stability of a rock tunnel in a dam site in China | |
Cheng et al. | Experimental study on the deformation and failure mechanism of overburden rock during coal mining using a comprehensive intelligent sensing method | |
Wang et al. | The volume of settlement trough change with depth caused by tunneling in sands | |
Zhu et al. | Study on the height of the mining-induced water-conducting fracture zone under the Q2l loess cover of the Jurassic coal seam in northern Shaanxi, China | |
Kanik et al. | Importance of numerical analyses for determining support systems in tunneling: A comparative study from the trabzon-gumushane tunnel, Turkey | |
CN116703244B (zh) | 采矿沉陷区治理效果与综合评价方法 | |
Zhang et al. | Tunnel stability analysis of coral reef limestone stratum in ocean engineering | |
Ma et al. | Experimental study on deformation and failure characteristics and monitoring and early warning of surrounding rock of tunnel crossing sliding surface | |
CN113536414A (zh) | 基于三维建模的岩质边坡稳定性分析方法、系统及介质 | |
Zhang et al. | Numerical simulation study on the influence of the ground stress field on the stability of roadways | |
Lin et al. | Disturbed granite identification by integrating rock mass geophysical properties | |
Sun et al. | Improvement of upper limit of mining under an aquifer of a super thick unconsolidated layer in Huainan based on multi-physics field monitoring | |
Zhang et al. | Stability of a roadway below a coal seam under dynamic pressure: a case study of the 11123 floor gas drainage roadway of a mine in huainan, China |
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 |