CN110390152B - 一种模拟巷道围岩裂隙演化的离散元方法 - Google Patents
一种模拟巷道围岩裂隙演化的离散元方法 Download PDFInfo
- Publication number
- CN110390152B CN110390152B CN201910635883.6A CN201910635883A CN110390152B CN 110390152 B CN110390152 B CN 110390152B CN 201910635883 A CN201910635883 A CN 201910635883A CN 110390152 B CN110390152 B CN 110390152B
- Authority
- CN
- China
- Prior art keywords
- rock
- roadway
- rock mass
- elastic modulus
- 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
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- 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
- 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/20—Design optimisation, verification or simulation
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0014—Type of force applied
- G01N2203/0016—Tensile or compressive
- G01N2203/0017—Tensile
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0014—Type of force applied
- G01N2203/0016—Tensile or compressive
- G01N2203/0019—Compressive
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/0202—Control of the test
- G01N2203/0212—Theories, calculations
- G01N2203/0218—Calculations based on experimental data
-
- 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)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Road Repair (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种模拟巷道围岩裂隙演化的离散元方法,包括现场煤岩层取芯并记录RQD值,观测巷道变形并统计煤层裂隙分布特征;室内测试煤岩试样的力学参数,根据RQD值计算出岩体的强度;利用UDEC‑Trigon模块建立数值模型调整参数与岩体强度相匹配校正模型参数;建立工程尺度的数值模型调整参数与现场变形特征匹配,最终模拟巷道围岩裂隙演化。本发明为离散元数值模拟巷道变形提供准确的基础力学参数,可确保模拟结果真实可靠。
Description
技术领域
本发明属于煤矿开采技术领域,具体涉及一种模拟巷道围岩裂隙演化的离散元方法。
背景技术
煤柱的稳定性一直以来都是安全高效开采需要研究的重点。随着采深增加地质条件复杂,沿空掘巷出现大变形、采空区泄漏瓦斯等现象,巷道围岩控制效果差。煤柱宽度留设不合理,对煤柱和围岩破坏机理的认识不清,导致不能提出长久有效的控制措施,严重的影响煤炭资源安全高效开采。
数值模拟方法在解决工程问题方面具有明显的优势,主要包括有限元和离散元两种方法。离散元数值模型能够模拟岩体的破坏过程并且显式的模拟裂隙扩展演化规律,广泛应用于解决地下岩土工程问题。UDEC Trigon模型是在泰森多边形的基础上进一步改进的非连续模型,模型由块体和节理面两部分组成,在不同的应力环境下沿着块体之间的节理面发生破坏,能够明显的区别原生裂隙和新产生的裂隙,可以更加真实的模拟巷道围岩的裂隙演化过程。
现有的数值模拟方法和技术在模拟巷道变形时存在:模型中岩体力学参数的选取是直接将实验室测得的力学参数应用到数值模型中;而由于实验室采用的试样通常较小,试样是由原岩经过开采、切割后形成,该过程中,内部应力及结构发生变化,从而与真实的大尺度岩体相比,存在较大差异;现有的数值模拟方法和技术忽略了大尺度的工程岩体与小尺度岩块之间的差异性,最终得到的结果与现场实际有较大的偏差,缺乏可信度,并不能指导现场施工。
发明内容
本发明的目的是提供一种模拟巷道围岩裂隙演化的离散元方法,能够有效真实地还原现场,模拟结果更加真实可靠。
为实现上述目的,本发明采用的一种模拟巷道围岩裂隙演化的离散元方法,包括以下步骤:
S1,在试验巷道内,进行观测巷道变形特征、煤岩样取岩芯;记录巷道围岩变形特征,监测巷道顶板和两帮的裂隙分布特征,在顶板和煤层内取岩芯,记录顶板岩芯及煤层岩芯的RQD值;
S2,将取得的顶板岩芯及煤层岩芯分别制成标准的岩块试样进行力学实验:测试岩块试样的单轴抗压强度、抗拉强度和弹性模量;根据岩体的弹性模量与岩块试样的弹性模量的比值与RQD之间的函数关系,以及岩体的单轴抗压强度和岩块试样单轴抗压强度之间的函数关系,确定岩体的单轴抗压强度和弹性模量;岩体的抗拉强度取岩体的抗压强度的0.1;
S3,应用UDEC中的Trigon模块建立的校正模型,模型中三角形块体的边长的平均值根据步骤S1中统计的裂隙分布特征来确定;以步骤S2中计算得到的岩体的单轴抗压强度、抗拉强度和弹性模量为已知特征值,利用校正模型反演得到校正模型中的三角形块体参数和节理面参数;
S4,利用UDEC建立大尺寸的与现场工作面开采规模一致的数值模型,限定研究区域,在研究区域内利用Trigon模块划分成为三角形块体,三角形块体的边长的平均值以及与三角形块体参数与第三步中校正模型中相同;在研究区域以外的区域划分为长度逐渐增大的矩形块体;以步骤S1中观测得到巷道围岩变形特征和裂隙分布规律为已知状态,反演得到大模型中矩形块体的力学参数。
根据本发明的一些实施例,步骤S3中,反演校正模型中的三角形块体参数和节理面参数,包括三角形块体尺寸、弹性模量和节理面的法向刚度、切向刚度、内聚力、内摩擦角、抗拉强度,具体的反演方法是:
a)分别建立宽高为2m×4m的矩形和直径为2m的圆形校正模型,确保三角形块体的边长的平均值与现场统计的裂隙长度相一致;
b)设置模型中三角形块体的弹性模量与公式(II)计算得到的岩体的弹性模量相同;
c)三角形块体之间节理面的法向刚度Kn根据公式(IV)进行推导,切向刚度Ks与Kn的比值为0.2,泊松比μ根据Ks/Kn的变化确定;公式(IV)的具体表达式为:
式中,K和G分别为三角形块体的体积模量和剪切模量,GPa,根据K=E/3(1-2μ)和G=E/2(1+μ)来计算,其中μ是块体的泊松比;E是块体的弹性模量,GPa;ΔZmin是相邻单元在垂直方向的最小宽度,m;n取10;
d)利用校正模型进行多次单轴压缩和巴西劈裂实验,得到全程应力应变曲线,采用迭代试错的方法使模拟得到的力学参数与步骤S2中计算得到的岩体的单轴抗压强度、抗拉强度和弹性模量相匹配,最终反演出模型的节理面的内聚力、内摩擦角和抗拉强度。
根据本发明的一些实施例,步骤S1中,利用十字交叉布点法安装巷道表面位移测站,采用防爆相机记录巷道围岩变形特征,应用钻孔窥视仪监测巷道顶板和两帮的裂隙分布特征;采用钻孔描图的方法来量化钻孔表面裂隙长度的分布规律。
根据本发明的一些实施例,步骤S1中,RQD是岩石质量指标,具体表达式为:
式中,l为10cm以上(含10cm)的岩芯累计长度,m;L为钻孔总长度。
根据本发明的一些实施例,步骤S2中,岩体的弹性模量与岩块弹性模量比值与RQD之间的函数关系,具体表达式为:
Em/Er=100.0186RQD-1.91 (II)
式中,Em为岩体的弹性模量,GPa;Er为岩块试样的弹性模量,GPa。
根据本发明的一些实施例,步骤S2中,岩体的单轴抗压强度和岩块单轴抗压强度之间的函数关系,具体表达式为:
式中,σcm和σc分别为岩体和岩块试样的单轴抗压强度,MPa;m取值为0.63。
本发明的有益技术成果:考虑工程岩体与岩块试样的差异性,根据岩石质量指标RQD与岩体和岩块弹性模量比值之间的函数关系,以及岩体和岩块单轴抗压强度之间的函数关系,来计算岩体的力学参数,采用UDEC中的Trigon模块建立校正模型,反演得到岩体的数值计算模型参数,能够真实可靠的模拟现场巷道变形特征。因此,本发明提出的离散元模拟巷道变形的参数校正方法填补了目前数值模拟参数确定难的空白,能够真实还原现场情况,使模拟结果更加真实可靠,在一定程度上为安全高效开采提供保障,有更强的现场指导意义。
附图说明
图1煤层的钻孔柱状图;
图2a煤柱帮的裂隙分布特征;
图2b实煤体帮的裂隙分布特征;
图3a煤柱内裂隙长度的概率分布图;
图3b实煤体内裂隙长度的概率分布图;
图4a单轴压缩校正模型;
图4b巴西劈裂校正模型;
图5a校正模型的单轴压缩全程应力应变曲线;
图5b校正模型的巴西劈裂全程应力应变曲线;
图6现场工程尺度的数值模型;
图7数值模拟巷道变形与现场观测对比图;
图8模拟得到的7m煤柱内裂隙分布;
图9本发明一种模拟巷道围岩裂隙演化的离散元方法的流程图。
具体实施方式
下面结合附图及实施例对本发明作进一步详细说明。应当说明,以下实施例仅用于对本发明的说明,并不限定本发明。
实施例以阳煤集团寺家庄矿15106回风巷7m窄煤柱沿空掘巷为工程背景,该煤层平均埋深574m,平均厚度为5.5m,节理裂隙发育,15106回风巷全长1700m,巷道沿煤层顶板掘进,掘进断面为矩形宽×高:4.8m×4.0m,煤层的钻孔柱状图如图1所示。现采用本发明所述的方法对模拟沿巷道变形的数值模型进行参数校正。本实例的模拟巷道围岩裂隙演化的离散元方法,包括以下几个步骤:
1)在15106回风巷内对完整的顶底板进行取样,记录不同岩性的RQD值。采用中国矿业大学研制生产的YTJ20型钻孔窥视仪对窄煤柱侧和实煤体侧进行探测,在距离底板2m位置处施工水平探测孔,钻孔直径29mm,在煤柱侧打4m深的钻孔,在实煤体侧打20m深的钻孔,观测得到裂隙的分布特征分别如图2a和2b所示。由图2a和2b可知,煤柱侧,在观测的4m范围内都出现大量的环形贯通裂隙和竖向裂隙,靠近煤柱帮部的煤体已经发生破坏,成为破碎区,出现塌孔的现象。煤柱的另一侧为采空区,无法采用钻孔窥视仪进行观测,但是根据钻孔窥视结果推断,煤柱靠近采空区侧的破坏程度应该与巷道侧类似或者更加严重。因此,可以确定15106回风巷在巷道掘进后,煤柱内部裂隙发育,整个煤柱已经发生破坏。实煤体侧,在距离煤帮1m范围内煤层破碎较严重,1m-2m范围内存在少量的裂隙,深度超过2m后煤层未受采动影响,几乎保存完整,没有发现较大的裂隙。
采用钻孔描图的方法,对煤柱内和实煤体内的裂隙长度进行统计分析,结果分别如图3a和图3b所示。由图3a和图3b可知,煤柱内和实煤体内的裂隙长度满足对数正态分布规律。在实煤体内,裂隙长度范围是0.01~0.23m,平均值为0.07m,其中长度小于0.2m的裂隙占总数的99.7%,几乎没有大裂隙(长度大于0.2m)。在煤柱内,裂隙长度范围是0.03~0.59m,平均值为0.13m,其中长度小于0.2m的裂隙占总数的80%,大裂隙占到总数的20%。对比分析图3a和3b可知,在未受采动影响的实煤体内不存在大裂隙,而煤柱中20%的大裂隙并不是天然煤体中的原始裂隙,而是煤柱在形成过程中发生破坏而产生的,这些大裂隙并不能代表天然的原生裂隙。因此,煤层内的原生裂隙长度要小于0.2m,为建立校正模型提供基础数据。
2)利用步骤1)取得的岩芯制作成标准试样:直径50mm、高度100mm的圆柱体和直径为50mm、厚度为25mm的小圆柱体进行单轴压缩试验和巴西劈裂实验,主要测试岩样的单轴抗压强度、抗拉强度和弹性模量。根据公式(II)和(III)来计算岩体的强度和弹性模量。岩体的抗拉强度取单轴抗压强度的0.1。测试结果如表1所示。
表1岩块和计算得到的岩体的力学参数
测试结果表明工程岩体与小尺度的岩块试样的力学参数有明显的差异,弹性模量和单轴抗压强度都有不同程度的降低,同时也表明将实验室小石块的力学参数直接带入数值模型中进行计算是不准确的,需要进一步的校正与工程岩体参数相匹配。
其中,岩体的弹性模量与岩块弹性模量比值与RQD之间的函数关系,具体表达式为:
Em/Er=100.0186RQD-1.91 (II)
式中,Em为岩体的弹性模量,GPa;Er为岩块试样的弹性模量,GPa。
岩体的单轴抗压强度和岩块单轴抗压强度之间的函数关系,具体表达式为:
式中,σcm和σc分别为岩体和岩块试样的单轴抗压强度,MPa;m取值为0.63。
3)应用UDEC中的Trigon模块建立小尺寸的校正模型,模型大小为宽2m高4m的矩形模型用于单轴压缩实验,校正模型如图4a所示,直径为2m的圆形模型用于巴西劈裂实验,校正模型如图4b所示。模型中的三角块边长的平均值取值0.2m,能够代表煤层未受扰动前的原生裂隙尺寸,有利于区别新产生的裂隙和原生裂隙。
对4种不同岩性的岩层分别建立模型进行校正。模型内的块体设置为弹性属性,可以发生变形但不会产生破坏。节理采用Mohr-Coulomb残余强度本构模型,岩体只能沿着节理面发生破坏,微观裂隙的产生、扩展、贯通逐渐形成宏观裂隙。
块体的密度与岩体的真实密度相同,通过实验室测得。块体的弹性模量与计算得到的岩体弹性模量相同。节理面的刚度根据模型块体单元的尺寸和Ks/Kn的比例来确定。采用迭代试错的方法来反演出模型节理面的内聚力、内摩擦角和抗拉强度。
利用校正模型进行一系列的单轴压缩和巴西劈裂实验,得到全程应力应变曲线,如图5a和5b所示。
模拟结果与计算得到的岩体力学参数(单轴抗压强度、抗拉强度和弹性模量)相匹配,最终反演出模型中使用的力学参数如表2。节理的残余强度准则认为破坏后的节理内聚力和内摩擦角都减小为0,只保留一定的残余内摩擦角。
表2校正后的UDEC-Trigon模型力学参数
力学参数的校正结果如表3所示。由表3可知,校正后的力学参数与目标值之间的误差都小于10%,也证明该参数校正方法的真实可靠性。
表3力学参数校正结果
4)根据15106工作面的地质条件,采用UDEC离散元软件建立与现场工作面开采规模一致的数值模型,模型长150m,高58m。为提高计算效率,在重点研究的区域利用Trigon模块划分成为三角块体,如图6所示。煤柱区域三角块长度为0.2m,巷道周边煤层区域三角块长度为0.4m,煤层直接顶和直接底区域三角块长度为0.5m。在研究区域以外的区域划分为块体长度不断增大的矩形块体,块体长度1.1m、2m和7m。这种块体尺寸划分方式能够有效地模拟煤柱的力学行为。
模型左右边界在水平方向固定,底部在垂直方向固定,顶部施加上覆岩层重量的等效压力,在模型内部施加原位测量的地应力场。模型参数按照表2中校正的参数进行赋值。开挖计划按照实际的开采顺序进行,先开挖15108工作面,模型计算平衡后再开挖15106回风巷,留设煤柱宽度为7m。巷道的支护方式按照实际的支护参数进行模拟。以步骤1)中观测得到巷道围岩变形特征和裂隙分布规律为已知状态,反演得到大模型中矩形块体的力学参数。
模拟得到的巷道变形与现场观测结果的对比如图7所示。由图可知,巷道在掘进期间就发生严重的变形,顶底板移近量约为592mm,两帮移近量为696mm,与现场距离巷口610m位置处观测的巷道变形(顶底板移近量540mm,两帮移近量645mm,图7b所示)相吻合。图7c~7g对比了数值模拟与现场实测的巷道破坏形式,结果表明现场巷道在采动应力的影响下产生顶板下沉、两帮向巷道空间挤出和底鼓,造成顶板钢带折断和底板铺设的轨道变形,其中煤柱帮变形最大,主要发生在煤柱下半部分,并且挤压运煤皮带。以上现场掘进期间观测到的巷道破坏形式都可以在UDEC Trigon模型中真实的捕捉到,如图7f所示。
UDEC Trigon模型模拟得到15106回风巷掘进后,7m宽的煤柱内裂隙分布情况如图8所示。根据裂隙的分布情况,煤柱被分为大裂隙区和小裂隙区。由图可知,大裂隙区出现在煤柱的两侧,其中距离巷道侧2m范围和距离采空区侧3m范围张拉裂隙发育接近破碎,成为大裂隙区。距离巷道3m~4m范围剪切裂隙发育,成为小裂隙区。模拟结果与现场钻孔窥视结果(图2a)非常吻合。
通过现场尺度的大模型的验证,表明使用校正的参数带入到UDEC Trigon模型中能够真实的模拟现场巷道的变形情况,利用校正后的参数可以进一步预测未来工作面开挖所产生的矿压现象。同时,证明本发明提出方法的可靠性和实用价值。
Claims (4)
1.一种模拟巷道围岩裂隙演化的离散元方法,其特征在于,包括以下步骤:
S1,在试验巷道内,进行观测巷道变形特征、煤岩样取岩芯;记录巷道围岩变形特征,监测巷道顶板和两帮的裂隙分布特征,在顶板和煤层内取岩芯,记录顶板岩芯及煤层岩芯的RQD值;
利用十字交叉布点法安装巷道表面位移测站,采用防爆相机记录巷道围岩变形特征,应用钻孔窥视仪监测巷道顶板和两帮的裂隙分布特征;采用钻孔描图的方法来量化钻孔表面裂隙长度的分布规律;
S2,将取得的顶板岩芯及煤层岩芯分别制成标准的岩块试样进行力学实验:测试岩块试样的单轴抗压强度、抗拉强度和弹性模量;根据岩体的弹性模量与岩块试样的弹性模量的比值与RQD之间的函数关系,以及岩体的单轴抗压强度和岩块试样单轴抗压强度之间的函数关系,确定岩体的单轴抗压强度和弹性模量;岩体的抗拉强度取岩体的抗压强度的0.1;
岩体的弹性模量与岩块弹性模量比值与RQD之间的函数关系,具体表达式为:
Em/Er=100.0186RQD-1.91 (II)
式中,Em为岩体的弹性模量,GPa;Er为岩块试样的弹性模量,GPa;
S3,应用UDEC中的Trigon模块建立的校正模型,模型中三角形块体的边长的平均值根据步骤S1中统计的裂隙分布特征来确定;以步骤S2中计算得到的岩体的单轴抗压强度、抗拉强度和弹性模量为已知特征值,利用校正模型反演得到校正模型中的三角形块体参数和节理面参数;
S4,利用UDEC建立大尺寸的与现场工作面开采规模一致的数值模型,限定研究区域,在研究区域内利用Trigon模块划分成为三角形块体,三角形块体的边长的平均值以及与三角形块体参数与第三步中校正模型中相同;在研究区域以外的区域划分为长度逐渐增大的矩形块体;以步骤S1中观测得到巷道围岩变形特征和裂隙分布规律为已知状态,反演得到大模型中矩形块体的力学参数。
2.根据权利要求1所述的一种模拟巷道围岩裂隙演化的离散元方法,其特征在于,步骤S3中,反演校正模型中的三角形块体参数和节理面参数,包括三角形块体尺寸、弹性模量和节理面的法向刚度、切向刚度、内聚力、内摩擦角、抗拉强度,具体的反演方法是:
a)分别建立宽高为2m×4m的矩形和直径为2m的圆形校正模型,确保三角形块体的边长的平均值与现场统计的裂隙长度相一致;
b)设置模型中三角形块体的弹性模量与公式(II)计算得到的岩体的弹性模量相同;
c)三角形块体之间节理面的法向刚度Kn根据公式(IV)进行推导,切向刚度Ks与Kn的比值为0.2,泊松比μ根据Ks/Kn的变化确定;公式(IV)的具体表达式为:
式中,K和G分别为三角形块体的体积模量和剪切模量,GPa,根据K=E/3(1-2μ)和G=E/2(1+μ)来计算,其中μ是块体的泊松比;E是块体的弹性模量,GPa;ΔZmin是相邻单元在垂直方向的最小宽度,m;n取10;
d)利用校正模型进行多次单轴压缩和巴西劈裂实验,得到全程应力应变曲线,采用迭代试错的方法使模拟得到的力学参数与步骤S2中计算得到的岩体的单轴抗压强度、抗拉强度和弹性模量相匹配,最终反演出模型的节理面的内聚力、内摩擦角和抗拉强度。
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910635883.6A CN110390152B (zh) | 2019-07-15 | 2019-07-15 | 一种模拟巷道围岩裂隙演化的离散元方法 |
PCT/CN2019/109932 WO2021007959A1 (zh) | 2019-07-15 | 2019-10-08 | 一种模拟巷道围岩裂隙演化的离散元方法 |
AU2019443098A AU2019443098B2 (en) | 2019-07-15 | 2019-10-08 | Discrete element method for modelling fracture evolution of roadway surrounding rock |
RU2020133633A RU2746748C1 (ru) | 2019-07-15 | 2019-10-08 | Метод дискретных элементов для моделирования развития разлома в породе, окружающей штрек |
US17/043,707 US11209415B2 (en) | 2019-07-15 | 2019-10-08 | Discrete element method for modelling a fracture evolution of a roadway surrounding rock |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910635883.6A CN110390152B (zh) | 2019-07-15 | 2019-07-15 | 一种模拟巷道围岩裂隙演化的离散元方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110390152A CN110390152A (zh) | 2019-10-29 |
CN110390152B true CN110390152B (zh) | 2021-04-02 |
Family
ID=68286508
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910635883.6A Active CN110390152B (zh) | 2019-07-15 | 2019-07-15 | 一种模拟巷道围岩裂隙演化的离散元方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US11209415B2 (zh) |
CN (1) | CN110390152B (zh) |
AU (1) | AU2019443098B2 (zh) |
RU (1) | RU2746748C1 (zh) |
WO (1) | WO2021007959A1 (zh) |
Families Citing this family (40)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110826254A (zh) * | 2019-11-27 | 2020-02-21 | 神华准格尔能源有限责任公司 | 露天矿边坡岩体节理分布参数测定方法、存储介质和系统 |
CN110887749B (zh) * | 2019-12-04 | 2022-04-15 | 宁波大学 | 结构面抗剪强度尺寸效应试验试样尺寸确定方法 |
CN111008486B (zh) * | 2019-12-27 | 2022-11-01 | 中国矿业大学 | 一种2d模拟综采面滚筒截割煤壁并统计块煤率的方法 |
CN111271129A (zh) * | 2020-01-31 | 2020-06-12 | 天地科技股份有限公司 | 获取回采巷道围岩变形与裂隙扩展规律的方法 |
CN111723421A (zh) * | 2020-06-04 | 2020-09-29 | 中国矿业大学 | 一种沿空留巷巷旁充填体宽度确定方法 |
CN111814242B (zh) * | 2020-07-16 | 2023-10-03 | 新疆工程学院 | 一种利用物联网数据判定沿空掘巷煤柱宽度的方法及系统 |
CN112001086B (zh) * | 2020-08-26 | 2024-01-30 | 湖南科技大学 | 基于再生顶板分层特性确定漏风裂隙胶结材料组分的方法 |
CN112906312A (zh) * | 2021-01-19 | 2021-06-04 | 中国矿业大学(北京) | 变节理刚度的裂隙煤岩体流固耦合模拟参数选取方法 |
CN113240635B (zh) * | 2021-05-08 | 2022-05-03 | 中南大学 | 一种裂缝分辨为基准的结构物检测图像质量测试方法 |
CN113153435B (zh) * | 2021-05-17 | 2022-07-01 | 中国矿业大学 | 一种双巷布置系统复采扰动下煤柱加固参数的确定方法 |
CN113433123B (zh) * | 2021-05-24 | 2024-06-04 | 中国电建集团贵阳勘测设计研究院有限公司 | 一种钻孔内地质参数测试方法及智能化装置 |
CN113536442A (zh) * | 2021-07-30 | 2021-10-22 | 新疆建筑科学研究院(有限责任公司) | 一种模拟互层岩体数值模型的仿真方法、系统及处理终端 |
CN113654477B (zh) * | 2021-08-16 | 2023-02-21 | 中国矿业大学 | 一种煤体变形测试装置、测试系统及测试方法 |
CN114033356B (zh) * | 2021-11-24 | 2023-12-12 | 中海石油(中国)有限公司 | 一种煤系地层地应力计算方法及装置 |
CN114233390B (zh) * | 2021-12-20 | 2023-07-25 | 中铁二院工程集团有限责任公司 | 一种适用于软弱围岩的辅助坑道封堵结构及其计算方法 |
CN114419982B (zh) * | 2021-12-29 | 2024-03-15 | 山东科技大学 | 软岩地层留煤柱采空区巷道变形破坏的模型试验系统及方法 |
CN116106091B (zh) * | 2022-01-20 | 2023-07-07 | 山东科技大学 | 深部含瓦斯煤体失稳灾变时的顶板弹性能贡献率估算方法 |
CN114547805B (zh) * | 2022-02-28 | 2024-06-04 | 华中科技大学 | 一种盾构隧道形状记忆合金柔性环间接头及其设计方法 |
CN114861401B (zh) * | 2022-04-08 | 2024-04-05 | 武汉大学 | 一种层状岩体fdem数值模拟输入参数标定方法 |
CN114893228B (zh) * | 2022-04-25 | 2024-06-07 | 河北工程大学 | 一种基于围岩破裂非对称性的分区分级超前主动锚固方法 |
CN114841003B (zh) * | 2022-05-12 | 2024-04-26 | 福建省水利水电勘测设计研究院有限公司 | 围岩最高温度反演计算的方法 |
CN115203900B (zh) * | 2022-06-17 | 2023-04-18 | 中国科学院武汉岩土力学研究所 | 基于离散元的裂隙岩体建模方法 |
CN115184579B (zh) * | 2022-06-20 | 2023-03-24 | 北京科技大学 | 一种空洞塌陷隐患全过程致灾演变试验装置及方法 |
CN115081302B (zh) * | 2022-07-15 | 2023-07-07 | 中国矿业大学 | 支护构件与硐室围岩接触及相互作用的模拟方法和系统 |
CN115355008B (zh) * | 2022-07-15 | 2023-09-19 | 中国矿业大学 | 一种冲击地压矿井掘进巷道煤岩穿层期间立体式卸压方法 |
CN116502386B (zh) * | 2022-09-09 | 2023-10-31 | 中国矿业大学 | 静动荷载作用下巷道锚固围岩灾变演化的模拟方法及系统 |
CN115600398B (zh) * | 2022-10-10 | 2023-06-16 | 长安大学 | 一种基于蒙特卡洛模拟的大型洞室岩体参数概率估计方法 |
CN115828519B (zh) * | 2022-11-02 | 2023-06-23 | 国网湖北省电力有限公司经济技术研究院 | 一种输电杆塔基础边坡土石混合体单轴抗压强度计算方法 |
CN115859714B (zh) * | 2022-11-23 | 2023-06-23 | 长安大学 | 一种基于fem-dem联合仿真的岩石爆破全过程模拟方法 |
CN115906525B (zh) * | 2022-12-29 | 2023-07-25 | 重庆大学 | 数值模拟岩层运动过程中力学参数映射关系的确定方法 |
CN115964901B (zh) * | 2023-03-16 | 2023-05-16 | 中国科学院地质与地球物理研究所 | 基于离散单元法的水致岩石强度劣化的模拟方法及系统 |
CN116467865B (zh) * | 2023-04-06 | 2023-11-24 | 中国地质大学(北京) | 考虑t应力及裂隙参数的岩体压剪断裂模型及其建立方法 |
CN116385688B (zh) * | 2023-06-01 | 2023-08-15 | 北京畅图科技有限公司 | 三维巷道模型快速构建方法、装置、计算机设备及介质 |
CN116861704B (zh) * | 2023-09-04 | 2023-11-14 | 北京交通大学 | 一种高地应力软岩隧道大变形等级动态快速预测方法 |
CN117172063A (zh) * | 2023-09-06 | 2023-12-05 | 山东能源集团有限公司 | 一种深孔爆破改性邻空巷道悬臂结构效果的评价方法 |
CN117217121B (zh) * | 2023-09-14 | 2024-06-14 | 国网新源集团有限公司 | 基于分布式并行topso-ewoa-de算法的岩体力学参数反演方法及系统 |
CN117251919B (zh) * | 2023-10-13 | 2024-05-17 | 湖南科技大学 | 一种沿空留巷顶板侧向破断位置判定方法 |
CN117763914B (zh) * | 2023-12-27 | 2024-06-21 | 河海大学 | 一种近断层的隧洞应力场和地震动力响应的数值分析方法、终端及介质 |
CN117706068B (zh) * | 2024-02-06 | 2024-04-19 | 湖南省通盛工程有限公司 | 一种桥梁基底岩石rqd测定方法、系统和存储介质 |
CN118095127B (zh) * | 2024-02-28 | 2024-07-26 | 中国矿业大学 | 一种水力压裂硬岩可截割性的综合表征方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2843184A3 (en) * | 2013-08-28 | 2016-01-06 | Services Petroliers Schlumberger | Method for performing a stimulation operation with proppant placement at a wellsite |
CN107290506A (zh) * | 2017-07-28 | 2017-10-24 | 中国石油大学(北京) | 一种定量评价储层成岩演化过程中孔隙度时空演变的方法 |
CN109684785A (zh) * | 2019-03-07 | 2019-04-26 | 湘潭大学 | 一种深部高应力巷道围岩动态损伤破坏演化方法及系统 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9152745B2 (en) * | 2012-03-06 | 2015-10-06 | Ion Geophysical Corporation | Model predicting fracturing of shale |
CN105093349A (zh) * | 2014-05-19 | 2015-11-25 | 中国矿业大学 | 一种实测巷道顶板内部裂隙发展发育规律的方法 |
US20170131192A1 (en) * | 2015-11-06 | 2017-05-11 | Baker Hughes Incorporated | Determining the imminent rock failure state for improving multi-stage triaxial compression tests |
CN106919770A (zh) * | 2017-03-21 | 2017-07-04 | 安徽理工大学 | 一种基于数值模拟的损伤变量确定方法 |
CN108489809B (zh) * | 2018-03-06 | 2020-06-02 | 西南石油大学 | 利用实验手段计算应力作用下粗糙错位裂缝变形量的方法 |
CN109271738B (zh) * | 2018-10-16 | 2023-05-12 | 山东科技大学 | 一种用于获取巷道围岩Weibull分布参数的数值反演方法 |
CN109492262A (zh) * | 2018-10-16 | 2019-03-19 | 山东科技大学 | 一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法 |
US11098582B1 (en) * | 2020-02-17 | 2021-08-24 | Saudi Arabian Oil Company | Determination of calibrated minimum horizontal stress magnitude using fracture closure pressure and multiple mechanical earth model realizations |
-
2019
- 2019-07-15 CN CN201910635883.6A patent/CN110390152B/zh active Active
- 2019-10-08 WO PCT/CN2019/109932 patent/WO2021007959A1/zh active Application Filing
- 2019-10-08 US US17/043,707 patent/US11209415B2/en active Active
- 2019-10-08 AU AU2019443098A patent/AU2019443098B2/en not_active Ceased
- 2019-10-08 RU RU2020133633A patent/RU2746748C1/ru active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2843184A3 (en) * | 2013-08-28 | 2016-01-06 | Services Petroliers Schlumberger | Method for performing a stimulation operation with proppant placement at a wellsite |
CN107290506A (zh) * | 2017-07-28 | 2017-10-24 | 中国石油大学(北京) | 一种定量评价储层成岩演化过程中孔隙度时空演变的方法 |
CN109684785A (zh) * | 2019-03-07 | 2019-04-26 | 湘潭大学 | 一种深部高应力巷道围岩动态损伤破坏演化方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
US20210263003A1 (en) | 2021-08-26 |
US11209415B2 (en) | 2021-12-28 |
AU2019443098A1 (en) | 2021-02-04 |
RU2746748C1 (ru) | 2021-04-20 |
CN110390152A (zh) | 2019-10-29 |
AU2019443098B2 (en) | 2022-09-22 |
WO2021007959A1 (zh) | 2021-01-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110390152B (zh) | 一种模拟巷道围岩裂隙演化的离散元方法 | |
Wu et al. | Numerical study of failure mechanisms and control techniques for a gob-side yield pillar in the Sijiazhuang coal mine, China | |
Li et al. | Numerical investigation of the effect of the location of critical rock block fracture on crack evolution in a gob-side filling wall | |
CN105631102B (zh) | 一种深部高应力巷道钻孔卸压参数的数值模拟确定方法 | |
Xia et al. | Estimating the geological strength index and disturbance factor in the Hoek–Brown criterion using the acoustic wave velocity in the rock mass | |
Zhang et al. | Squeezing failure behavior of roof‐coal masses in a gob‐side entry driven under unstable overlying strata | |
Ding et al. | Stability of large underground caverns excavated in layered rock masses with steep dip angles: a case study | |
CN108827233B (zh) | 一种两层采空区地面沉降的预测方法 | |
CN112035937B (zh) | 一种隧道施工围岩变形分析及信息化施工方法 | |
CN113803066A (zh) | 无煤柱自成巷平衡开采设计方法及系统 | |
WO2023197821A1 (zh) | 一种基于渗透-损伤关系的覆岩损伤分区及高度确定方法 | |
Zhu et al. | Coupled discrete element–finite difference method for analyzing subsidence control in fully mechanized solid backfilling mining | |
Shi et al. | Stability control of gob-side entry retaining in fully mechanized caving face based on a compatible deformation model | |
Liu et al. | Nondestructive testing on cumulative damage of watery fractured rock mass under multiple cycle blasting | |
CN111680896A (zh) | 一种煤矿地下水库安全距离确定方法 | |
Hong et al. | The spatial distribution of excavation damaged zone around underground roadways during blasting excavation | |
Xia et al. | Numerical study of the mechanism of dense drillings for weakening hard roofs and its application in roof cutting | |
Wu et al. | Failure characteristics and cooperative control strategies for gob-side entry driving near an advancing working face: a case study | |
da Cunha | Scale effects in rock masses—A general report on the papers submitted to the Lisbon Workshop | |
Trinh et al. | Continued rock stress and displacement measurements combined with numerical modeling as an active, realistic rock engineering tool | |
CN113090280B (zh) | 一种双巷布置系统巷道底鼓机理的揭示方法 | |
Khosravi et al. | Evaluation of disc cutter performance in rock cutting process using 3D finite element method | |
CN113281176B (zh) | 水压致裂法测量结果的校验方法及加载装置 | |
Mao et al. | Exploration and verification of tunnel stability evolution law under jointed rock mass with various attitudes | |
CN114547746A (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 |