CN111324939A - 一种基于flac3d数值模拟矿产开采的地表下沉量预算方法 - Google Patents
一种基于flac3d数值模拟矿产开采的地表下沉量预算方法 Download PDFInfo
- Publication number
- CN111324939A CN111324939A CN201811529691.9A CN201811529691A CN111324939A CN 111324939 A CN111324939 A CN 111324939A CN 201811529691 A CN201811529691 A CN 201811529691A CN 111324939 A CN111324939 A CN 111324939A
- Authority
- CN
- China
- Prior art keywords
- model
- boundary
- mineral
- mining
- stress
- 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
- 229910052500 inorganic mineral Inorganic materials 0.000 title claims abstract description 33
- 239000011707 mineral Substances 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000004088 simulation Methods 0.000 title claims abstract description 18
- 238000005065 mining Methods 0.000 claims abstract description 30
- 239000011435 rock Substances 0.000 claims abstract description 20
- 238000012360 testing method Methods 0.000 claims abstract description 6
- 238000006073 displacement reaction Methods 0.000 claims description 8
- 230000005484 gravity Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 2
- 239000003245 coal Substances 0.000 description 8
- 238000009412 basement excavation Methods 0.000 description 3
- 238000005553 drilling Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种基于FLAC3D数值模拟的地表下沉量计预算方法，对实际揭露矿产的钻孔进行分析，结合实际开采时矿产的地质条件，以采掘工作面走向、倾向长度及选取钻孔矿层顶底板岩层厚度为模型的X、Y、Z长度，建立数值模型，选取摩尔‑库伦物理模型准则，矿产开采后的模型定义为空模型，并将模型底部定义为全约束边界，左右边界、前后边界为约束边界条件，模拟开采过程中，在模型顶部按上覆岩层的重力施加向下的应力，利用节点号输出地表下沉量数值。本发明结合物理力学测试，利用FLAC3D数值模拟软件，建立分析、调试模型，执行矿产开采模拟求解，获取模型响应结果，量化矿产开采对地表沉陷的影响，操作简单且易于实际应用。
Description
技术领域
本发明涉及煤层、金属矿产等领域矿层开采过程中地表下沉量预算方法，具体涉及一种基于FLAC3D数值模拟矿产开采地表下沉量的预算方法。
背景技术
矿层采出后，采区围岩原始应力平衡遭到破坏，使得岩层产生移动和变形，当开采面积达到一定范围后，岩层移动便会扩展到地表，从而引起地表的下沉。针对矿层开采造成的地表下沉量一般是在实际开采过程中进行的实时监测，难以对其进行事前的预测，而室内相似试验模拟不仅费用大，且耗时耗力。
发明内容
本发明的目的就是提供一种省时省力、成本低，能在未进行矿产开采时即可预测矿层开采后对地表下沉的影响的一种基于FLAC3D数值模拟矿产开采的地表下沉量预算方法。
本发明的一种基于FLAC3D数值模拟矿产开采的地表下沉量预算方法，包括，
步骤1：分析揭露矿产的钻孔，确定实际开采时矿产的采厚、埋深、覆岩厚度、采掘工作面走向、倾向长度地质条件，根据所确定的上述地质条件，以采掘工作面走向、倾向长度及选取钻孔矿层顶底板岩层厚度为模型的X、Y、Z长度；
步骤2：根据矿产及其顶底板岩性力学测试获取各岩层的厚度、密度、内摩擦角、抗拉强度、泊松比、弹性模量基础参数，另外，岩层的体积模量和剪切模量根据泊松比和弹性模量进行计算得到，其计算公式如下：
其中：K—体积模量，G—剪切模量，E—弹性模量，υ—泊松比；
步骤3：根据摩尔-库伦物理模型准则，矿产开采后该岩层的物理模型定义为空模型，并定义如下的模型边界条件并进行模型平衡分析，
(1)位移边界
①将模型底部定义为全约束边界，即节点在X、Y和Z方向上的位移均为0；
②左右边界、前后边界为约束边界条件，即X和Y方向均固定，水平位移为0)，Z方向为自由边界，不予约束；
(2)应力边界
①模型顶部为地表，故上边界初始垂直应力条件为0；
②模拟开采过程中，为便于模型研究，将上覆岩层的重力假设为均匀荷载，在模型顶部按上覆岩层的重力，施加向下的应力，具体公式如下：
q＝Σγh，
式中：q——模型上部边界的应力荷载，
γ——开采矿层上覆岩体的平均容重，
h——开采矿层上覆岩体的厚度，
另，模型X和Y方向上的应力条件按侧压系数λ(一般取泊松比的倒数)与垂直应力的乘积施加水平应力，侧压系数的变化范围如下：
式中：H为采深，
根据上述边界条件、应力条件的设定，对模型进行运行处理；
步骤4，通过对矿产开采的模拟，利用节点编号输出地表下沉量数值，其输入语言为：print gp dis range id xx any id xx any，必要时，可通过获取的下沉量数据，结合Excel作地表下沉曲线，进行综合分析。
本发明的一种基于FLAC3D数值模拟矿产开采的地表下沉量预算方法，依据勘探钻孔实测资料，结合该地区地质条件，利用FLAC3D数值模拟软件建立地质模型，通过对模型调试、加载，满足矿产实际开采情况，模拟矿层开采后对地表沉陷的影响，以数表的形式更直观的显示地表的下沉量，整个模拟操作简单，经济合理，易于实际应用。
附图说明
图1为本发明的流程图；
图2为本发明实例某矿矿层开采地表下沉量模拟地质模型；
图3为模拟开采后地表下沉量预算数据图(id为地表节点号，X-Dis、Y-Dis为地表水平运移，Z-Dis为地表下沉量。
具体实施方式
一种基于FLAC3D数值模拟矿产开采的地表下沉量预算方法，具体步骤是：
1、依据实际地质条件，综合分析揭露矿产的钻孔，建立合适的模型，根据以往资料的收集，分析揭露矿层的钻孔，选取地表以下适当的深度，建立适合的地质模型；实例以ZK1002孔实际揭露岩性为模拟对象，设计模拟尺寸为100m*100m*220m，模拟煤层为5-3#、7#、15-1#、15-2#，煤层厚分别为1.5m、2.5m、3.0m及2.0m，地质模型如图2所示。
2、根据矿产及其顶底板岩性力学测试获取相关参数，对钻孔揭露的矿层顶底板岩性进行采样，测试，获取模拟岩层的物理力学参数(包括：密度ρ、粘聚力C、内摩擦角φ、抗拉强度T、泊松比υ、弹性模量E)，体积模量和剪切模量可根据以下公式计算得到，
体积模量
剪切模量
3、结合所建立的模型，选取合适的物理模型准则，定义恰当的模型边界条件并进行模型平衡分析，具体为：
①模型准则：对模型选取摩尔—库伦准则，矿层开挖后设为空模型；
②模型边界条件：底部边界采用全约束边界条件，即X、Y和Z方向位移均为零；左右边界、前后边界为约束边界条件，即X方向、Y方向均固定(水平位移为零)，Z方向为自由边界，不予约束。
4、根据实际地质条件，施加应力条件，运行处理，
模型上边界为地表，其初始垂直应力为0，煤层开挖过程中按模拟岩体总厚对模型上边界施加应力(即：σz＝γH，γ取模型底部上覆岩层的加权平均)，即5.5MPa，测压系数为0.33。
5、模拟矿产开采，利用节点编号输出地表下沉量数值，
模拟煤层为5-3#、7#、15-1#、15-2#，煤层厚分别为1.5m、2.5m、3.0m及2.0m，设计开挖100m，煤层两端均留有20m的保护煤柱，开挖步距设为20m。
Claims (2)
1.一种基于FLAC3D数值模拟矿产开采的地表下沉量预算方法，其特征在于：它包括以下步骤：
步骤1：分析揭露矿产的钻孔，确定实际开采时矿产的采厚、埋深、覆岩厚度、采掘工作面走向、倾向长度地质条件，根据所确定的上述地质条件，以采掘工作面走向、倾向长度及选取钻孔矿层顶底板岩层厚度为模型的X、Y、Z长度；
步骤2：根据矿产及其顶底板岩性力学测试获取各岩层的厚度、密度、内摩擦角、抗拉强度、泊松比、弹性模量基础参数，另外，岩层的体积模量和剪切模量根据泊松比和弹性模量进行计算得到，其计算公式如下：
其中：K—体积模量，G—剪切模量，E—弹性模量，υ—泊松比；
步骤3：根据摩尔-库伦物理模型准则，矿产开采后该岩层的物理模型定义为空模型，并定义如下的模型边界条件并进行模型平衡分析，
(1)位移边界
①将模型底部定义为全约束边界，即节点在X、Y和Z方向上的位移均为0；
②左右边界、前后边界为约束边界条件，即X和Y方向均固定，水平位移为0)，Z方向为自由边界，不予约束；
(2)应力边界
①模型顶部为地表，故上边界初始垂直应力条件为0；
②模拟开采过程中，为便于模型研究，将上覆岩层的重力假设为均匀荷载，在模型顶部按上覆岩层的重力，施加向下的应力，具体公式如下：
q＝Σγh，
式中：q——模型上部边界的应力荷载，
γ——开采矿层上覆岩体的平均容重，
h——开采矿层上覆岩体的厚度，
根据上述边界条件、应力条件的设定，对模型进行运行处理；
步骤4：通过对矿产开采的模拟，利用节点编号输出地表下沉量数值。
2.根据权利要求1所述的一种基于FLAC3D数值模拟矿产开采的地表下沉量预算方法，其特征在于：步骤3中所述模型X和Y方向上的应力条件按侧压系数λ与垂直应力的乘积施加水平应力，侧压系数λ取泊松比的倒数，侧压系数的变化范围如下：
式中：H为采深。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811529691.9A CN111324939A (zh) | 2018-12-14 | 2018-12-14 | 一种基于flac3d数值模拟矿产开采的地表下沉量预算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811529691.9A CN111324939A (zh) | 2018-12-14 | 2018-12-14 | 一种基于flac3d数值模拟矿产开采的地表下沉量预算方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111324939A true CN111324939A (zh) | 2020-06-23 |
Family
ID=71170077
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811529691.9A Pending CN111324939A (zh) | 2018-12-14 | 2018-12-14 | 一种基于flac3d数值模拟矿产开采的地表下沉量预算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111324939A (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113027529A (zh) * | 2021-03-22 | 2021-06-25 | 中国矿业大学 | 一种基于数值反演的采场冲击危险性评估方法 |
| CN113486575A (zh) * | 2021-06-21 | 2021-10-08 | 武汉科技大学 | 一种地表变形破坏预测预警方法、装置及存储介质 |
| CN115983052A (zh) * | 2023-03-21 | 2023-04-18 | 中国矿业大学(北京) | 一种适用于flac3d中多密度岩层条件下的初始应力生成方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103093048A (zh) * | 2013-01-14 | 2013-05-08 | 中国矿业大学 | 一种岩层移动数值模拟自动化建模方法 |
| CN105631102A (zh) * | 2015-12-24 | 2016-06-01 | 河南理工大学 | 一种深部高应力巷道钻孔卸压参数的数值模拟确定方法 |
| CN106919770A (zh) * | 2017-03-21 | 2017-07-04 | 安徽理工大学 | 一种基于数值模拟的损伤变量确定方法 |
| CN108827233A (zh) * | 2018-09-17 | 2018-11-16 | 中国地质大学（北京） | 一种两层采空区地面沉降的预测方法 |
-
2018
- 2018-12-14 CN CN201811529691.9A patent/CN111324939A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103093048A (zh) * | 2013-01-14 | 2013-05-08 | 中国矿业大学 | 一种岩层移动数值模拟自动化建模方法 |
| CN105631102A (zh) * | 2015-12-24 | 2016-06-01 | 河南理工大学 | 一种深部高应力巷道钻孔卸压参数的数值模拟确定方法 |
| CN106919770A (zh) * | 2017-03-21 | 2017-07-04 | 安徽理工大学 | 一种基于数值模拟的损伤变量确定方法 |
| CN108827233A (zh) * | 2018-09-17 | 2018-11-16 | 中国地质大学（北京） | 一种两层采空区地面沉降的预测方法 |
Non-Patent Citations (1)
| Title |
|---|
| 陈四宝 等: "基于FLAC3D数值模拟的硅石灰开采巷道和采场支护技术探讨" * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113027529A (zh) * | 2021-03-22 | 2021-06-25 | 中国矿业大学 | 一种基于数值反演的采场冲击危险性评估方法 |
| CN113027529B (zh) * | 2021-03-22 | 2022-03-01 | 中国矿业大学 | 一种基于数值反演的采场冲击危险性评估方法 |
| CN113486575A (zh) * | 2021-06-21 | 2021-10-08 | 武汉科技大学 | 一种地表变形破坏预测预警方法、装置及存储介质 |
| CN115983052A (zh) * | 2023-03-21 | 2023-04-18 | 中国矿业大学(北京) | 一种适用于flac3d中多密度岩层条件下的初始应力生成方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Adhikary et al. | Modelling of longwall mining-induced strata permeability change | |
| Gao et al. | Evaluation of coal longwall caving characteristics using an innovative UDEC Trigon approach | |
| Akin et al. | Empirical correlations of shear wave velocity (Vs) and penetration resistance (SPT-N) for different soils in an earthquake-prone area (Erbaa-Turkey) | |
| US9164194B2 (en) | Method for modeling deformation in subsurface strata | |
| US8423337B2 (en) | Method for multi-scale geomechanical model analysis by computer simulation | |
| Chen et al. | Ground characterization using breaking-action-based zoning analysis of rotary-percussive instrumented drilling | |
| Albertz et al. | Critical state finite element models of contractional fault-related folding: Part 2. Mechanical analysis | |
| CN111324939A (zh) | 一种基于flac3d数值模拟矿产开采的地表下沉量预算方法 | |
| Huang et al. | 3-D discontinuum numerical modeling of subsidence incorporating ore extraction and backfilling operations in an underground iron mine in China | |
| Sainoki et al. | Methodology for the interpretation of fault-slip seismicity in a weak shear zone | |
| Jalbout et al. | Rock mechanics tools for mining in high stress ground conditions at Nickel Rim South Mine | |
| Jianjun et al. | Deformation behavior of mining beneath flat and sloping terrains in mountainous areas | |
| CN108304959B (zh) | 提高地层流体压力预测精度的方法 | |
| Wu et al. | Experimental validation of core-based pre-stress evaluations in rock: a case study of Changchikeng sandstone in the Tseng-wen reservoir transbasin water tunnel | |
| US20220050224A1 (en) | Evaluating anisotropic effective permeability in rock formations having natural fracture networks | |
| Ramiah et al. | A 2D geomechanical model of an offshore gas field in the Bredasdorp Basin, South Africa | |
| Sapigni et al. | Engineering geological characterization and comparison of predicted and measured deformations of a cavern in the Italian Alps | |
| Fernandez-Rodriguez et al. | Characterization of intensely jointed rock masses by means of in situ penetration tests | |
| Wang et al. | Investigation of mine pressure and deformation due to phosphate ore body excavation based on Hoek-Brown model | |
| Shrestha | Stability of tunnels subject to plastic deformation-a contribution based on the cases from the Nepal Himalaya | |
| Mandal et al. | Viscoelastic approach to capture varying least principal stress magnitude and the effect of observed stress layering on hydraulic fracturing-An example from shale formations of the Perth Basin | |
| Azarfar | Numerical and experimental analysis of the effect of faults in open-pit mining stability | |
| Li | Optimization of multistage hydraulic fracturing treatment for maximization of the tight gas productivity | |
| Pradhan | Stability analysis of open PIT slope using FLAC | |
| Chandra et al. | Effect of coal seam presence in-between comparatively denser layers over high wall slope stability-a case study of Kistaram Project, Singareni Collieries Company Limited (SCCL) |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20200623 |