CN113742940B - 一种采动卸压边界时空曲线簇确定方法 - Google Patents
一种采动卸压边界时空曲线簇确定方法 Download PDFInfo
- Publication number
- CN113742940B CN113742940B CN202111084544.7A CN202111084544A CN113742940B CN 113742940 B CN113742940 B CN 113742940B CN 202111084544 A CN202111084544 A CN 202111084544A CN 113742940 B CN113742940 B CN 113742940B
- Authority
- CN
- China
- Prior art keywords
- pressure relief
- boundary
- stress
- rock
- overlying strata
- 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
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
-
- 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
-
- 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)
Abstract
本发明公开了一种采动卸压边界时空曲线簇确定方法,具体包括以下步骤:三维数值模型计算、块散组合体本构关系推演、采动卸压位置应力边界时空演化过程确定和卸压边界时空曲线簇绘制。本发明充分考虑了时间和空间因素对卸压边界的影响,可以保证采动卸压边界的真实性和客观性,能够获得更为准确的采动卸压边界,这对于确定被保护层消突范围具有重要意义。
Description
技术领域
本发明涉及安全工程技术领域,特别是涉及一种采动卸压边界时空曲线簇确定方法。
背景技术
在煤炭开采过程中瓦斯灾害依旧是煤矿的主要灾害之一,极易造成群死群伤,社会影响极大,煤与瓦斯突出是其中典型的瓦斯灾害之一,随着煤炭开采向深部推进,煤层地应力、瓦斯压力显著增大,渗透率逐渐降低,岩体的流变特性更加明显,更易发生煤与瓦斯突出。卸压开采是煤与瓦斯突出防治最有效最经济的方法,卸压边界的准确判定对于确定邻近煤层卸压范围有至关重要的影响。目前,卸压边界的确定主要通过特定岩层顶底板相对位移量的测量或者采用弹性波瞬变电磁法来间接测量,上述方法测定的卸压边界不准确,难以有效确定消突范围。
因此,亟需开发一种采动卸压边界时空曲线簇确定方法。
发明内容
本发明的目的是提供一种采动卸压边界时空曲线簇确定方法,以解决现有的卸压边界不准确的问题。
为实现上述目的,本发明提供一种采动卸压边界时空曲线簇确定方法,具体包括以下步骤:
S1、构建工作面的三维数值计算模型,并对所述三维数值计算模型进行开采计算,获取覆岩的初始应力分布;
S2、推演采动下垮落带与裂隙带形成的覆岩块散组合体的考虑时间因素的本构关系,即卸压边界位置不同岩层的本构关系;
S3、计算所述覆岩块散组合体水平推力沿煤层倾向的分布,获得卸压边界位置的应力边界条件;
S4、将所述应力边界条件代入卸压边界位置不同岩层的本构关系中进行解算,得到卸压边界位置的应力随时间演化规律;
S5、取特定时间内垮落带、裂隙带内块散组合体的水平推力分布,进而获得卸压边界位置的应力边界条件,最终获得不同层位岩层的卸压位置,将其连成线,即可获得特定时间的卸压边界分布;
S6、基于所述S5,并取不同时间点,从而获得不同时刻的卸压边界曲线,即卸压边界动态演化时空曲线簇。
优选的,所述S1具体为:
S1.1、根据工作面的地质条件,构建三维数值计算模型;
S1.2、确定所述三维数值计算模型的相关参数,对所述三维数值计算模型进行开挖计算;
S1.3、基于所述S1.1~S1.2,获取煤岩体的初始应力分布。
优选的,所述三维数值计算模型采用摩尔-库伦屈服准则进行相关参数计算。
优选的,所述S2具体为:
S2.1、通过岩石力学实验仪器,得到所述覆岩块散组合体在应力恢复过程中岩层位移随时间的演化关系式;
S2.2、通过3D-DIC装置在像素级别上观察所述覆岩块散组合体在径向变形场的变化,并计算其径向应力,以此确定本构关系中的拟合参数,从而得到所述覆岩块散组合体的本构关系。
优选的,所述S3具体为:
S3.1、基于所述初始应力分布,建立覆岩卸压力学模型;
S3.2、对达到初始平衡后的所述覆岩力学模型进行工作面开采,获得垂直应力云图和位移云图,进而得到上层的垂直应力分布;
S3.3、利用摩擦系数计算上层岩层摩擦力,通过累加岩层厚度计算出下层岩层摩擦力,再通过水平推力,计算出覆岩块散组合体一侧断裂块体的水平推力;然后通过覆岩块散组合体由一侧断裂块体向每层的水平推力进行循环迭代,得到卸压位置应力值空间演化过程,即卸压边界位置的应力边界条件。
优选的,所述S4具体为:
S4.1、将所述应力边界条件代入所述卸压边界位置不同岩层的本构关系中进行解算,得到各岩层边界卸压范围量化表达式;
S4.2、将实际垂直应力作为参量,构建典型组合体扰动损伤流变力学模型,即所述覆岩块散组合体的应力随时间演化的规律。
优选的,所述S6具体为:
S6.1、通过所述典型组合体扰动损伤流变力学模型,逐层计算某一时刻所述裂隙带和垮落带内各层岩层的水平推力;
S6.2、将各层的实际垂直应力进行叠加,得到卸压边界位置处的垂直应力,结合所述S6.1得到的水平推力,计算卸压边界位置处煤岩实际应力,在各层卸压范围内选取实际应力等于原岩应力的点;
S6.3、将所述S6.2得到的点进行连接,得到一条覆岩组合体卸压边界动态演化曲线;
S6.4、取不同时刻得到的覆岩组合体卸压边界动态演化曲线,得到覆岩组合体的卸压边界动态演化曲线簇。
与现有技术相比,本发明具有以下技术效果:
本发明充分考虑了空间和时间因素对卸压边界的影响,以采空区中心点为计算中心分别计算卸压边界位置的应力条件,然后代入本构模型获得不同时期的卸压边界簇,实现对采动卸压边界进行进一步准确估计,基本保障了采动卸压边界的真实性和客观性,能够获得较为真实的采动卸压边界,这对于消除煤与瓦斯突出危险性具有重要意义。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例的方法流程图;
图2为本发明实施例的两带覆岩组合体卸压边界时空演化模型图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例1
参照图1-2所示,本发明提出一种采动卸压边界时空曲线簇确定方法,选取特定矿井作为研究对象,收集相关的矿区地质和开采数据,确定卸压位置应力边界演化特征,代入组合体扰动损伤流变力学模型中,得出特定时间特定采动条件下的卸压边界曲线。
具体包括以下步骤:
S1、构建工作面的三维数值计算模型,并通过所述三维数值计算模型进行开采计算,得到覆岩初始应力分布;
以实际工作面地质条件为依据,建立三维数值计算模型,并通过摩尔-库伦屈服准则来进行相关参数计算,确定该模型的几何尺寸倾向、走向、垂直高度等参数和保护层工作面的采高,并分析获取覆岩的应力变化,即初始应力分布;
S2、将采动下垮落带与裂隙带形成的覆岩散块组合体作为研究对象,推演其考虑时间因素的本构关系,即卸压边界位置不同岩层的本构关系;
在垮落带与裂隙带形成的覆岩组合体中取单个的块散连接组合体试样,通过岩石力学实验仪器,对块散连接组合体试样持续加载轴向压力至固定值,并持续施压一定时间,观察试样位移情况,并将试样位移情况作为因变量,时间作为自变量进行拟合,从而得到应力恢复过程中岩层位移形变随时间演化关系式,即块散连接组合体流变响应特征。同时煤矿井下的开采活动势必引起岩层的移动也有水平推力的作用,可以通过3D-DIC装置在像素级别上精准观察试件在径向上变形场的变化并计算径向应力,再通过煤岩的力学特点,表征煤岩的本构关系,即卸压边界位置不同岩层的本构关系,以此确定本构关系中的拟合参数,从而间接反映出水平推力随时间的演化规律。
S3、计算所述覆岩块散组合体水平推力沿煤层倾向的分布,获得卸压边界位置的应力边界条件;
基于所述S1得到的初始应力分布,建立覆岩卸压力学模型;所述覆岩卸压力学模型选取煤岩采动裂隙带范围内某一岩梁,建立覆岩卸压力学模型。其中,选取的岩梁为理想弹性体、水平无倾角,并且忽略了构造应力、岩梁端部载荷及其反力对移动形态的影响,且只考虑岩梁受载荷、自重应力及其反力的影响。
通过位移函数法求解力学平面问题中函数表达式所需的采空区各数据,并采用数学计算软件对所述覆岩卸压力学模型进行初始平衡解算,通过初始平衡解算使初始应力模型处于原始应力状态。模型达到初始平衡后开始进行工作面的正常开采模拟,模拟开采推进沿模型的长度方向(X轴负方向)进行,推进完成后分别获取相应的垂直应力云图、位移云图,以此算出上层的垂直应力(上下岩层成分相同且不变,摩擦系数恒定);利用摩擦系数计算出上摩擦力后,此时,由于力可以完全从上部传到下部,从而累加岩层厚度计算出下层摩擦力,再结合水平推力即可计算出实际推力,再通过单个块体到边界进行循环累加从而对卸压位置应力边界空间演化过程进行量化,即卸压边界位置的应力边界条件。
S4、将所述应力边界条件代入卸压边界位置不同岩层的本构关系中进行解算,得到卸压边界位置的应力随时间演化规律;
将S3中所得的覆岩组合体的边界条件代入所推演的块散组合体所得的本构关系,推演卸压边界处的水平推力作用,确定卸压边界,代入卸压边界处不同岩石的本构模型,得到各岩层边界卸压范围量化表征公式;将实际垂直应力与卸压范围作为研究参量,推演得到典型组合体扰动损伤流变力学模型,实现对覆岩组合体卸压位置应力边界时间演化进行量化表征,所述覆岩组合体的应力随时间演化的规律。
S5、取特定时间内垮落带、裂隙带内块散组合体的水平推力分布,进而获得卸压边界位置的应力边界条件,最终获得不同层位岩层的卸压位置,将其连成线,即可获得特定时间的卸压边界分布;
S6、基于所述S5,并取不同时间点,从而获得不同时刻不同位置的卸压边界曲线,即卸压边界动态演化时空曲线簇。
通过典型组合体扰动损伤流变力学模型,逐层计算所研究裂隙带、垮落带内各层岩层,依次得到各层卸压边界位置实际垂直应力水平推力;将各层计算所得垂直应力进行叠加,在各层卸压范围内选取垂直应力等于原岩应力的点,并将各层卸压范围内原岩应力点进行连接,从而得到一条组合体卸压边界动态演化曲线,随着时间的变化煤层不同位置的应力会发生变化,卸压效果会不同,与原岩应力大小相等的点的位置也会随之发生变化。将不同时刻各层卸压范围内原岩应力点进行连接会得到多条组合体卸压边界动态演化曲线,从而形成组合体卸压边界动态演化曲线簇。
以上所述的实施例仅是对本发明的优选方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。
Claims (4)
1.一种采动卸压边界时空曲线簇确定方法,其特征在于,具体包括以下步骤:
S1、构建工作面的三维数值计算模型,并对所述三维数值计算模型进行开采计算,获取覆岩的初始应力分布;
S2、推演采动作用下垮落带与裂隙带形成的覆岩块散组合体的考虑时间因素的本构关系,即不同岩层卸压边界位置的本构关系;
S3、计算所述覆岩块散组合体水平推力沿煤层倾向的分布,获得卸压边界位置的应力边界条件;
S4、将所述应力边界条件代入卸压边界位置不同岩层的本构关系中进行解算,得到卸压边界位置的应力随时间演化规律;
S5、取特定时间块散组合体的水平推力分布,进而获得卸压边界位置的应力边界条件,最终获得不同层位岩层的卸压位置,将其连成线,即可获得特定时间的卸压边界分布;
S6、基于所述S5,并取不同时间点,从而获得不同时刻的卸压边界曲线,即卸压边界动态演化时空曲线簇;
所述S3具体为:
S3.1、基于所述初始应力分布,建立覆岩卸压力学模型;
S3.2、对达到初始平衡后的所述覆岩卸压力学模型进行工作面开采,获得垂直应力云图和位移云图,进而得到上层的垂直应力分布;
S3.3、利用摩擦系数计算上层岩层摩擦力,通过累加岩层厚度计算出下层岩层摩擦力,再通过水平推力,计算出覆岩块散组合体一侧断裂块体的水平推力;然后通过覆岩块散组合体由一侧断裂块体向每层的水平推力进行循环迭代,得到卸压位置应力值空间演化过程,即卸压边界位置的应力边界条件;
所述S4具体为:
S4.1、将所述应力边界条件代入所述卸压边界位置不同岩层的本构关系中进行解算,得到各岩层边界卸压范围量化表达式;
S4.2、将实际垂直应力作为参量,构建典型组合体扰动损伤流变力学模型,即所述覆岩块散组合体的应力随时间演化的规律;
所述S6具体为:
S6.1、通过所述典型组合体扰动损伤流变力学模型,逐层计算某一时刻所述裂隙带和垮落带内各层岩层的水平推力;
S6.2、将各层的实际垂直应力进行叠加,得到卸压边界位置处的垂直应力,结合所述S6.1得到的水平推力,计算卸压边界位置处煤岩实际应力,在各层卸压范围内选取实际应力等于原岩应力的点;
S6.3、将所述S6.2得到的点进行连接,得到一条覆岩组合体卸压边界动态演化曲线;
S6.4、取不同时刻得到的覆岩组合体卸压边界动态演化曲线,得到覆岩组合体的卸压边界动态演化曲线簇。
2.根据权利要求1所述的采动卸压边界时空曲线簇确定方法,其特征在于,所述S1具体为:
S1.1、根据工作面条件,构建三维数值计算模型;
S1.2、确定所述三维数值计算模型的相关参数,对所述三维数值计算模型进行开挖计算;
S1.3、基于所述S1.1~S1.2,获取采动后覆岩初始应力分布。
3.根据权利要求2所述的采动卸压边界时空曲线簇确定方法,其特征在于,所述三维数值计算模型采用摩尔-库伦屈服准则进行相关参数计算。
4.根据权利要求1所述的采动卸压边界时空曲线簇确定方法,其特征在于,所述S2具体为:
S2.1、通过岩石力学实验仪器,得到所述覆岩块散组合体在应力恢复过程中岩层位移随时间的演化关系式;
S2.2、通过3D-DIC装置在像素级别上观察覆岩块散组合体在变形场的变化,并确定其对应应力,以此确定本构关系中的拟合参数,从而得到所述覆岩块散组合体的本构关系。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111084544.7A CN113742940B (zh) | 2021-09-16 | 2021-09-16 | 一种采动卸压边界时空曲线簇确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111084544.7A CN113742940B (zh) | 2021-09-16 | 2021-09-16 | 一种采动卸压边界时空曲线簇确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113742940A CN113742940A (zh) | 2021-12-03 |
CN113742940B true CN113742940B (zh) | 2023-09-08 |
Family
ID=78739285
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111084544.7A Active CN113742940B (zh) | 2021-09-16 | 2021-09-16 | 一种采动卸压边界时空曲线簇确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113742940B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105301202A (zh) * | 2015-11-24 | 2016-02-03 | 山东科技大学 | 用于确定上保护层开采卸压范围的测试系统及测试方法 |
CN109374425A (zh) * | 2018-12-17 | 2019-02-22 | 郑州轻工业学院 | 含瓦斯煤体卸压后发生层裂破坏的模拟研究实验装置及实验方法 |
WO2020006818A1 (zh) * | 2018-07-02 | 2020-01-09 | 山东科技大学 | 一种煤岩体分区注水渗流-损伤-应力耦合数值模拟方法 |
CN111812296A (zh) * | 2020-06-16 | 2020-10-23 | 重庆大学 | 一种卸压开采相似模拟方法 |
CN111860952A (zh) * | 2020-06-16 | 2020-10-30 | 重庆大学 | 一种突出煤层关键开采参数快速优选方法 |
-
2021
- 2021-09-16 CN CN202111084544.7A patent/CN113742940B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105301202A (zh) * | 2015-11-24 | 2016-02-03 | 山东科技大学 | 用于确定上保护层开采卸压范围的测试系统及测试方法 |
WO2020006818A1 (zh) * | 2018-07-02 | 2020-01-09 | 山东科技大学 | 一种煤岩体分区注水渗流-损伤-应力耦合数值模拟方法 |
CN109374425A (zh) * | 2018-12-17 | 2019-02-22 | 郑州轻工业学院 | 含瓦斯煤体卸压后发生层裂破坏的模拟研究实验装置及实验方法 |
CN111812296A (zh) * | 2020-06-16 | 2020-10-23 | 重庆大学 | 一种卸压开采相似模拟方法 |
CN111860952A (zh) * | 2020-06-16 | 2020-10-30 | 重庆大学 | 一种突出煤层关键开采参数快速优选方法 |
Non-Patent Citations (1)
Title |
---|
双保护层开采条件下上覆煤岩层裂隙场分布与演化特性研究;刘应科;问小江;蒋名军;;徐州工程学院学报(自然科学版)(第01期);17-23 * |
Also Published As
Publication number | Publication date |
---|---|
CN113742940A (zh) | 2021-12-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gao et al. | Field experiments on fracture evolution and correlations between connectivity and abutment pressure under top coal caving conditions | |
Ju et al. | Numerical analysis of the dynamic evolution of mining-induced stresses and fractures in multilayered rock strata using continuum-based discrete element methods | |
CN103982179B (zh) | 一种油藏储层的古压力定量反演探测方法 | |
Xu et al. | Discrete element modeling of strata and surface movement induced by mining under open-pit final slope | |
CN110135030B (zh) | 一种采空区地表沉降的预测方法 | |
Liu et al. | Sensitivity analysis of the main factors controlling floor failure depth and a risk evaluation of floor water inrush for an inclined coal seam | |
US20130275101A1 (en) | Method For Modeling Deformation In Subsurface Strata | |
Hosseini et al. | Passive seismic velocity tomography on longwall mining panel based on simultaneous iterative reconstructive technique (SIRT) | |
Zeng et al. | Experimental investigation on mining-induced strain and failure characteristics of rock masses of mine floor | |
Albertz et al. | Critical state finite element models of contractional fault-related folding: Part 2. Mechanical analysis | |
Huang et al. | In situ stress distribution and its impact on CBM reservoir properties in the Zhengzhuang area, southern Qinshui Basin, North China | |
Zhang et al. | Cooperative monitoring and numerical investigation on the stability of the south slope of the Fushun west open-pit mine | |
CN108021763A (zh) | 一种老采空区残余沉降预测方法 | |
Fuenkajorn et al. | Prediction of cavern configurations from subsidence data | |
Huang et al. | 3-D discontinuum numerical modeling of subsidence incorporating ore extraction and backfilling operations in an underground iron mine in China | |
Zhao et al. | Dynamic reduction of rock mass mechanical parameters based on numerical simulation and microseismic data–a case study | |
Cheng et al. | Model for the prediction of subsurface strata movement due to underground mining | |
Oyanguren et al. | Stability analysis of Llerin Rockfill Dam: An in situ direct shear test | |
Tan et al. | Influence of complicated faults on the differentiation and accumulation of in-situ stress in deep rock mass | |
CN114692456A (zh) | 一种诱发矿井强动力灾害地应力临界判据的确定方法 | |
CN113742940B (zh) | 一种采动卸压边界时空曲线簇确定方法 | |
CN111767642B (zh) | 薄松散层采煤沉陷区地基稳定性评价方法及装置 | |
CN108732023B (zh) | 一种大采高采场煤壁超前支承压力峰值系数的确定方法 | |
Ma et al. | The influence of backfilling of a caved zone and magma intrusion-type faults on surface deformation in a metal mine caused by sublevel caving mining | |
CN110689612B (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 |