img

QQ群聊

img

官方微信

高级检索

黄金科学技术 ›› 2019, Vol. 27 ›› Issue (6): 851-861.doi: 10.11872/j.issn.1005-2518.2019.06.851

• 采选技术与矿山管理 • 上一篇    下一篇

大新锰矿复杂空区群三维数值模型构建方法及胶结充填治理研究

高远(),陈庆发(),蒋腾龙   

  1. 广西大学资源环境与材料学院,广西 南宁 530004
  • 收稿日期:2019-04-19 修回日期:2019-08-11 出版日期:2019-12-31 发布日期:2019-12-24
  • 通讯作者: 陈庆发 E-mail:gygxgjy2018@163.com;chqf98121@163.com
  • 作者简介:高远(1983-),男,山东淄博人,硕士研究生,从事三维数字矿山研究工作。gygxgjy2018@163.com
  • 基金资助:
    广西自然科学基金联合资助培育项目“金属矿床协同采矿方法”(2018GXNSFAA138105)

Research on Three-Dimensional Numerical Model Construction Method and Cemented Filling Treatment of Complex Cavities in Daxin Manganese Mine

Yuan GAO(),Qingfa CHEN(),Tenglong JIANG   

  1. College of Resources,Environment and Materials,Guangxi University,Nanning 530004,Guangxi,China
  • Received:2019-04-19 Revised:2019-08-11 Online:2019-12-31 Published:2019-12-24
  • Contact: Qingfa CHEN E-mail:gygxgjy2018@163.com;chqf98121@163.com

摘要:

矿山岩土工程数值分析中研究复杂空区群三维数值模型构建方法对空区治理具有现实意义。为精确分析矿山地下工程围岩移动和应力变化规律,特别是不同采矿方法形成的复杂空区的稳定性和矿柱应力变化,提出了3DMine-Surfer-Rhino-ANSYS-FLAC3D多软件联合建模方法,解决了复杂空区群三维数值模型构建和大数量小尺寸矿房矿柱网格剖分问题。以广西大新锰矿矿体开采为案例,构建了留矿法与房柱法形成的采空区的三维数值模型,模拟分析空区的稳定性和决策胶结充填治理方案。研究结果表明:(1)大新锰矿留矿法开采的空区部分顶柱有垮落风险;(2)房柱法开采的空区部分矿柱变形大,有塑性破坏;(3)对于胶结充填留矿法的空区,隔一充一充填方案治理效果更好,在治理初期地表沉降更小。

关键词: 空区群建模, 数值模拟, 空区胶结充填治理, 3DMine, FLAC3D, 留矿法采空区, 房柱法采空区

Abstract:

In the numerical analysis of mine geotechnical engineering,it is of practical significance to study the construction method of 3D numerical model of complex cavities for the treatment of cavities.The pre-processing ability of geotechnical engineering simulation software is poor,and the current modeling method can’t well build complex geological body numerical model,and cannot fine divide the grid,resulting in the difficulty and work intensity of geotechnical engineering researchers in the early stage of modeling.In order to accurately analyze the law of surrounding rock movement and stress change of underground engineering, especially the stability of complex goaf formed by different mining methods and the stress change of ore pillar,the paper proposed and improved the 3DMine-surfer-rhino-ANSYS-FLAC3D multi-software joint modeling method.The method can solve the problems of 3D numerical model construction of complex goaf group and mesh generation of large number and small size ore pillars.The paper also described the specific modeling steps.The complex cavity group model constructed is closer to the reality,including the undulating form of surface DEM and more refined mesh of cavity and pillar model. At the same time,the total number of model grid elements is reduced,and the grid size of local key research areas can be adjusted according to the needs,so as to improve the accuracy of the calculation of key areas and accurately describe the stress change of ore pillars.This modeling method can be used to flexibly control the sequence excavation of different layers,stages and chambers of the layered ore bodies,and can be used to simulate the production sequence according to the engineering design at the early stage or the actual engineering production sequence at the later stage,so as to find a safe and reasonable mining production plan.A three-dimensional numerical model of goaf formed by shrinkage stoping method and chamber and pillar method was established by taking the mining of manganese ore body in Daxin mine of Guangxi as an example. Research results show that:(1) The Daxin manganese ore of shrinkage stoping method in mining area,empty part of the needle deformation is large,there are risks caving,the surface subsidence deformation is high;(2) The method of room and pillar mining zone,part of the pillar plastic failure,the surface subsidence deformation amount is larger;(3) Method of shrinkage stoping with cemented filling empty area:Fill in the blanks by stage filling area,after filling in the fourth stage,the surface settlement decreased significantly.In the vertical and vertical consolidation filling scheme,after filling one ore chamber with two ore chambers,the surface subsidence is smaller than that before filling,and the area with large subsidence becomes discontinuous. In the scheme of filling one ore chamber with another ore chamber,the surface subsidence deformation is one order of magnitude lower than the former two schemes,and the area of deformation area decreases and the continuity decreases.Among the three filling schemes,the one-filling and one-filling schemes have better governance effect,and the surface is more stable at the initial stage of governance.

Key words: cavities modeling, numerical simulation, cavities consolidation and filling treatment, 3DMine, FLAC3D, goafs of shrinkage stoping method, goafs of room and pillar mining

中图分类号: 

  • TD32

图1

大新锰矿矿体与采空区分布图1-西北地采采区矿体;2-西北采区矿体;3-中部采区矿体;4-西北地采房柱法开采空区;5-西北地采留矿法开采空区;6-西北采区房柱法开采空区;7-中部采区留矿法开采空区"

图2

大新锰矿采空区稳定性分析建模流程图"

图3

Surfer与Rhino软件处理过程"

图4

Workbench软件中由面模型到实体模型"

图5

Workbench软件中完整模型"

图6

Workbench软件中网格剖分"

图7

FLAC3D西北采区模型网格"

图8

3个采区模拟计算结果"

表1

矿岩体物理力学参数"

名称体积模量/GPa剪切模量/GPa黏聚力/MPa内摩擦角/(°)抗拉强度/MPa密度/(g·m-3)
矿体20.8018.304.6036.726.283.33
上盘23.0010.010.1431.351.392.63
下盘14.246.194.3135.111.632.64
充填体0.5970.330.1025.000.331.75

图9

胶结充填方案对比"

表2

胶结充填成本与充填体强度"

胶固粉与尾砂配比动力与人工费/(元·m-3)尾矿砂/(元·m-3)水/(元·m-3)胶固粉/(元·m-3)充填成本/(元·m-3)90 d抗压强度/MPa
1∶44.554.801.83666.077.1865.65
1∶64.555.161.70146.257.6114.22
1∶84.555.341.67436.347.8642.96
1∶104.555.461.67429.741.3841.82

图10

不同方案的充填效果对比"

1 陈庆发,周科平,古德生. 协同开采与采空区协同利用[J].中国矿业,2011,20(12):77-80.
Chen Qingfa,Zhou Keping,Gu Desheng.Synergetic mining and cavity synergetic utilization[J].China Mining Magazine,2011,20(12):77-80.
2 刘科伟,李夕兵,宫凤强,等. 基于CALS及Surpac-FLAC3D耦合技术的复杂空区稳定性分析[J].岩石力学与工程学报,2008,27(9):1924-1931.
Liu Kewei,Li Xibing,Gong Fengqiang,et al.Stability analysis of complicated 3D technology cavity based on CALS and coupled Surpac-FLAC3D technology[J]. Chinese Journal of Rock Mechanics and Engineering,2008,27(9):1924-1931.
3 寇向宇,贾明涛,王李管,等.基于CMS及DIMINE-FLAC3D耦合技术的采空区稳定性分析与评价[J].矿业工程研究,2010,25(1):31-35.
Kou Xiangyu,Jia Mingtao,Wang Liguan,et al.Evaluation and analysis of stability of mined-out area based on the CMS and DIMINE-FLAC3D coupling technique[J]. Mineral Engineering Research,2010,25(1):31-35.
4 宋卫东,付建新,杜建华,等.基于精密探测的金属矿山采空区群稳定性分析[J].岩土力学,2012,33(12):3781-3787.
Song Weidong,Fu Jianxin,Du Jianhua,et al. Analysis of stability of goaf group in metal mines based on precision detection[J]. Rock and Soil Mechanics,2012,33(12):3781-3787.
5 崔芳鹏,胡瑞林,刘照连,等.基于Surfer平台的FLAC3D复杂三维地质建模研究[J].工程地质学报,2008,16(5):699-702.
Cui Fangpeng,Hu Ruilin,Liu Zhaolian,et al. Surfer software platform based on complex three-dimensional geological models for pre-processing of FLAC3D[J]. Jounal of Engineering Geology,2008,16(5):699-702.
6 闫长斌,徐国元,李夕兵.爆破震动对采空区稳定性影响的FLAC3D分析[J].岩石力学与工程学报,2005,24(16):2894-2899.
Yan Changbin,Xu Guoyuan,Li Xibing. Stability analysis of mined-out areas influenced by blasting vibration with FLAC3D[J]. Chinese Journal of Rock Mechanics and Engineering,2005,24(16):2894-2899.
7 刘德峰.东升庙采空区稳定性的监测、预报及控制技术研究[D].包头:内蒙古科技大学,2014.
Liu Defeng. Research on Monitoring Forecast and Controlling Technology of the Goaf in Dongshengmiao[D].Baotou: Inner Mongolia University of Science and Technology,2014.
8 郑文棠,张勇平,李明卫,等. 基于三维可视化模型的高边坡演化过程分析[J].河海大学学报(自然科学版),2009,37(1):66-70.
Zheng Wentang,Zhang Yongping,Li Mingwei,et al. Evolution process of high slopes based on 3D visualization model[J]. Journal of Hohai University( Natural Sciences),2009,37(1):66-70.
9 廖秋林,曾钱帮,刘彤,等.基于ANSYS平台复杂地质体FLAC3D模型的自动生成[J].岩石力学与工程学报,2005,24(6):1010-1013.
Liao Qiulin,Zeng Qianbang,Liu Tong,et al. Automatic model generation of complex geologic body with FLAC3D based on ANSYS platform[J]. Chinese Journal of Rock Mechanics and Engineering,2005,24(6):1010-1013.
10 马新根,孙健东,赵洪泽,等.基于Rhino3D NURBS的FLAC3D快速建模方法在露天边坡问题中的应用[J].煤矿安全,2015,46(8):138-141.
Ma Xingen,Sun Jiandong,Zhao Hongze,et al. Application of FLAC3D rapid modeling in slope problem of open-pit mine based on Rhino3D NURBS[J]. Safety in Coal Mines,2015,46(8):138-141.
11 李翔,王金安,张少杰.复杂地质体三维数值建模方法研究[J].西安科技大学学报,2012,32(6):676-681.
Li Xiang,Wang Jin’an,Zhang Shaojie.3D modeling method of complicated geological body[J]. Journal of Xi’an University of Science and Technology,2012,32(6):676-681.
12 王树仁,慎乃齐,张海清,等.下伏采空区高速公路隧道变形特征数值分析[J].中国矿业,2008,17(3):76-79,81.
Wang Shuren,Shen Naiqi,Zhang Haiqing,et al. Numerical analysis of deformation characteristics for highway tunnels through the mined-up regions[J]. China Mining Magazine,2008,17(3):76-79,81.
13 罗周全,杨彪,刘晓明,等.基于CMS实测及Midas-FLAC3D耦合的复杂空区群稳定性分析[J].矿冶工程,2010,30(6):1-5.
Luo Zhouquan,Yang Biao,Liu Xiaoming,et al.Stability analysis for goaf group based on CMS and coupling of Midas-FLAC3D[J]. Mining and Metallurgical Engineering,2010,30(6):1-5.
14 王敏,陈晓艳,柯波,等.基于软件数据耦合精细建模的采区稳定性分析研究[J].世界科技研究与发展,2016,38(4):768-772.
Wang Min,Chen Xiaoyan,Ke Bo,et al. Mining area stability analysis based on software data-coupled and refined modeling[J]. World SCI-TECH R&D,2016,38(4):768-772.
15 Hanri1937.基于Gocad-surfer-犀牛-ANSYS的FLAC3D建模步骤[EB/OL].,2018.
Hanri1937. FLAC3D modeling steps based on Gocad-Surfer-Rhinoceros-ANSYS[EB/OL].,2018.
16 陈庆发,杨家彩,高远,等.大型复杂地质体三维数值模型构建方法比较研究[J].岩土力学,2016,37(增2):753-760.
Chen Qingfa,Yang Jiacai,Gao Yuan,et al.Comparative study of construction method for 3D numerical model of large complex geologic body[J]. Rock and Soil Mechanics,2016,37(Supp.2):753-760.
17 Bock S. New open-source ANSYS-SolidWorks-FLAC3D geometry conversion programs[J].Journal of Sustainable Mining,2015,14:124-132.
18 杨宁,尹贤刚,钟勇,等.基于3DMine-FLAC3D的复杂采空区稳定性分析[J].化工矿物与加工,2019(4):22-24,28.
Yang Ning,Yin Xiangang,Zhong Yong,et al. Stability analysis of complex goaf based on 3DMine-FLAC3D[J]. Industrial Minerals and Processing,2019(4):22-24,28.
19 李家琪,杨帆,杜森.基于复杂地形的等高线FLAC模拟矿区沉陷预计[J].测绘与空间地理信息,2018,41(2):200-202,206.
Li Jiaqi,Yang Fan,Du Sen.Prediction of mining subsidence based on FLAC simulation of complex terrain[J]. Geomatics and Spatial Information Technology,2018,41(2):200-202,206.
20 王飞飞,邹平,孟中华,等. 基于三维数值分析的大桥磷矿采空区稳定性研究[J]. 有色金属(矿山部分),2019,71(2):66-70,81.
Wang Feifei,Zou Ping,Meng Zhonghua,et al.Stability of goaf of Daqiao phosphorite mine based on 3D numerical analysis[J]. Nonferrous Metals(Mining Section),2019,71(2):66-70,81.
21 卢欣奇,李雪峰,张勤斌,等. 基于PS-InSAR技术的老采空区地表沉陷监测与分析[J]. 中国矿业,2019,28(4):104-110,114.
Lu Xinqi,Li Xuefeng,Zhang Qinbin,et al.Surface subsidence monitoring of old goaf based on the PS-InSAR technology[J]. China Mining Magazine,2019,28(4):104-110,114.
22 马凤山,卢蓉,郭捷,等. 金川二矿区大体积充填体变形的三维数值模拟[J]. 工程地质学报,2019,27(1):14-20.
Ma Fengshan,Lu Rong,Guo Jie,et al. Deformation analysis of large backfill by three-dimensional numerical simulation in No.2 zone of Jinchuan Mine[J]. Journal of Engineering Geology,2019,27(1):14-20.
[1] 谢也真,曹平,陈昊然. 滥泥坪铜矿三维地应力测量及巷道布置优化研究[J]. 黄金科学技术, 2019, 27(6): 862-870.
[2] 田军, 刘建坡, 杨勇, 张长银. 进路充填法爆破扰动诱发充填体破坏规律研究[J]. 黄金科学技术, 2019, 27(5): 687-695.
[3] 宋恩祥, 李强, 张静, 彭康. 蚀变带内矿体开采中人工假底的应用研究[J]. 黄金科学技术, 2019, 27(5): 722-730.
[4] 唐玉成, 方翔, 张中雷, 崔晓荣. 多规格石高强度开采矿山开拓系统的三维规划设计与施工[J]. 黄金科学技术, 2019, 27(5): 755-761.
[5] 杨仕教,王志会. 上覆公路浅埋采空区群稳定性数值模拟[J]. 黄金科学技术, 2019, 27(4): 505-512.
[6] 梁瑞,俞瑞利,周文海,黄小彬,王建勇,熊征宇. 空气间隔装药爆破模型在矿房回采中的应用[J]. 黄金科学技术, 2019, 27(3): 358-367.
[7] 秦世康,陈庆发,尹庭昌. 岩石与岩体冻融损伤内涵区别及研究进展[J]. 黄金科学技术, 2019, 27(3): 385-397.
[8] 刘锋,王昭坤,马凤山,王波,王剑波,董春蕾. 矿山深部卸压技术研究现状及展望[J]. 黄金科学技术, 2019, 27(3): 425-432.
[9] 王春,王成,熊祖强,程露萍,王怀彬. 动力扰动下深部出矿巷道围岩的变形特征[J]. 黄金科学技术, 2019, 27(2): 232-240.
[10] 方传峰,王晋淼,李剡兵,贾明涛. 基于PFC2D-DFN的自然崩落法数值模拟研究[J]. 黄金科学技术, 2019, 27(2): 189-198.
[11] 吕闹,汪海波. 厚硬顶板弱化前后垮落致灾数值模拟研究[J]. 黄金科学技术, 2019, 27(2): 257-264.
[12] 王从陆,李童. 基于Ventsim的深井线性热源对有效通风量的影响分析[J]. 黄金科学技术, 2019, 27(2): 271-277.
[13] 李响,怀震,李夕兵,张倬瑶. 基于裂纹扩展模型的脆性岩石破裂特征及力学性能研究[J]. 黄金科学技术, 2019, 27(1): 41-51.
[14] 陈冲,李夕兵,冯帆. 诱导巷道的围岩松动破坏区数值研究[J]. 黄金科学技术, 2018, 26(6): 771-779.
[15] 刘连生,钟清亮,闫雷,钟文,梁龙华. 频繁生产爆破加载下饱水岩体累积损伤效应声波测试研究[J]. 黄金科学技术, 2018, 26(6): 750-760.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!