Analysis of the Numerical Simulation to Structure Parameter of Hard-To-Mine Ore Body in Xincheng Gold Mine
-
摘要: 为研究新城金矿极破碎难采矿体用上向水平进路充填法回采时的稳定性，采用ANSYS有限元软件，选择Ⅲ盘区4种结构参数进行了不同开采时期的分析和比较，从而获得了最佳进路规格为5m× 7m，一步盘区超前准一步盘区5个分层，即25m，也与矿山试验采场的参数相吻合。Abstract: In order to study the stability of hard-to-mine ore body when using overhand drift mining method, we made analysis and comparison on four different structure parameter of Ⅲ mining area in different mining period by using the software ANSYS. In the end, we obtain the best standard of the drift 5m× 7m and the distance which first mining area exceed the second one, ie. 25m which are similar with the parameter of the cut-and-try stop.
-
Keywords:
- hard-to-mine ore body /
- numerical simulation /
- ANSYS /
- structure parameter
-
0 引言
新城金矿是我国20世纪70年代发现的蚀变岩型金矿床的典型代表，除具有成矿规模大、工业价值高的优点外，也存在着开采技术条件复杂等问题，尤其是-480m到-530m水平的Ⅴ#矿体属于极破碎难采矿体，经多方论证和比较，采用垂直走向布置的机械化盘区上向进路充填采矿法进行开采。本文采用ANSYS软件对采场结构参数进行数值模拟。
1 数值模拟的必要性
为了研究Ⅴ#极破碎难采矿体回采时的稳定性，并确定合理的回采结构参数，必须采用数值模拟的方法，定量地计算和分析回采过程中围岩中的应力、位移和和破坏的分布状况，确定它们在假定条件下随回采进行的动态变化过程，从而对围岩的稳定性状态做出判断；通过对多种方案的计算和分析，比较不同的回采结构参数下的围岩稳定性状态，做出最优的开采设计，确定最合理的采场结构参数, 从而保证安全生产，提高生产效益。
当前，数值模拟的方法有很多种，如有限元、离散元、边界元、无界元以及流行元等[1]。在已有的数值模拟方法中，有限元方法是比较适合模拟和分析地下开挖稳定性问题的方法。有限元模拟分析的基本思路是：首先将要分析的连续体分割成数目有限的单元，在结点上引进等效力以代替实际作用于单元上的力，对于每个单元根据分块近似的思想，选择一个简单函数来近似地表示位移分量的分布规律，并建立单元结点力和位移之间的关系，再把所有单元的这种特性关系集合起来，就得到一组以结点位移为未知量的代数方程式，由此方程组求出物体上有限个离散结点上的位移分量，进而根据物理方程和几何方程来求得单元的应力和变形。三维有限元数值模拟可以定量地计算和分析回采过程中采场围岩中的应力、位移和塑性区、破坏区的分布状况，确定它们随每步回采的动态变化过程，从而对采场围岩的稳定性状态做出判断[2]。软件ANSYS具有多种有限元分析的能力，包括从简单线性静态分析到复杂的非线性瞬态动力学分析。一个典型的ANSYS分析过程可分为建立模型、加载并求解、查看分析结果3个步骤。
2 计算力学模型
为了达到数值模拟的准确性，在选取参数时采用新城金矿前期工作所测得的矿区岩体的物理力学参数；对矿区工程地质、水文地质特别是对矿区岩体结构面，节理、裂隙的调查、统计和分析结果作为程序的初始条件。岩体的初始应力条件是指由岩体自重而引起的应力状态。模型的垂直与水平应力荷载分别是由原地应力场来决定的。在这个阶段，将要分步开挖的矿岩及不同的岩石性质都将分成不同的块体。在开挖时只要将要挖掉的块体删除或把其材料属性重新设置为零即可。
2.1 岩体与充填体力学参数
岩体物理力学参数是根据矿区内不同的岩体分别选取的，为简化计算，分析时将井下介质分为4类：上盘岩体、矿体、下盘岩体和充填体。由于岩体节理裂隙问题单独考虑，所以用岩石试块平均抗压强度代替岩体的强度参数。其不同类型的岩体物理力学性质如表 1所示。
表 1 不同岩体物理力学参数
2.2 原岩应力场
新城金矿地应力场的规律是：最大主应力为竖直方向，其他两个主应力一个是水平方向垂直矿体走向，另一个是水平方向且与矿体走向一致。模拟中取最大主应力方向为Y方向，水平方向为X方向，另一个为Z方向，原岩应力场经实测并回归如下：
式中：σmax——最大水平主应力；
2.3 计算模型尺寸
为了使数值模拟计算充分反映采场开挖的准确性，拟对两种不同的采场结构模型进行计算，这两个结构模型分别是Ⅲ盘区模型和矿体整体模型，运用ANSYS中的平面应变模型SOLID单元plane42对其进行分析。
Ⅲ盘区模型尺寸为30m×50m，即沿走向长30m，埋藏深度从-480m水平到-530m水平，高50m。如图 1所示。
矿体整体模型尺寸为120m×50m, 即沿走向长120m，，埋藏深度从-480m水平到-530m水平，高50m。如图 2所示。
2.4 计算方案
当采用上向进路充填采矿法进行采场开挖时，破坏了围岩应力固有的平衡，引起应力的重新分布，进而达到新的平衡。但随着开采的不断进行，应力不断改变，围岩应力场不断调整，并处于一个周而复始的动态平衡中，直到采矿活动的停止。从力学的角度看，采矿过程是引起一系列围岩应力动态平衡的原因，因此，围岩中的力学行为状态与采场的结构参数有着密切的关系，同样也可以在对矿区整体和单个采场采用上向进路充填采矿法开采时的应力、应变、位移计算中，研究极破碎矿体的采场结构参数。
2.4.1 Ⅲ盘区模型计算方案
为了确定盘区内合理的采场结构参数，拟对Ⅲ盘区采场4种结构参数（P1~4）进行计算，并认为进路开采时每采一条充填一条。其结构参数如表 2。
表 2 采场结构参数 m
2.4.2 矿体整体模型计算方案
数值模拟的目的之一就是从宏观的角度分析整个矿体在开采过程中因不同的盘区开采顺序而导致的不同应力分布情况。在实际生产中为了增大生产力，往往是Ⅰ、Ⅲ盘区先行开采，Ⅱ、Ⅳ盘区紧随其后。根据经验，Ⅰ、Ⅲ盘区相对Ⅱ、Ⅳ盘区超前的层数对整个矿体在开采过程中的应力分布有很大的影响，故矿体整体模型计算方案制定如下：
方案一：Ⅰ、Ⅲ盘区比Ⅱ、Ⅳ盘区超前3层。即Ⅰ、Ⅲ盘区先用上向进路法进行回采，当采至第4层时，Ⅱ、Ⅳ盘区也开始回采，比Ⅰ、Ⅲ盘区滞后3个分层。
方案二：Ⅰ、Ⅲ盘区比Ⅱ、Ⅳ盘区超前5层。即Ⅰ、Ⅲ盘区先用上向进路法进行回采，当采至第6层时，Ⅱ、Ⅳ盘区也开始回采，比Ⅰ、Ⅲ盘区滞后5个分层。
3 计算结果分析
计算采用ANSYS有限元程序中的平面应变模型，对所选方案进行了模拟研究。
3.1 Ⅲ盘区模型计算分析
由于所要分析的Ⅲ盘区处在整个矿体中，当矿体的其他盘区进行开采时，必然将引起整个应力场的变化，所分析的Ⅲ盘区也将有不同的变形情况出现，所以有必要从整体出发，来分析Ⅲ盘区进行开采时的应力、应变情况。边界条件如下：两侧水平方向位移为零，-530m水平、垂直方向位移为零。上部承受13.2MPa的压力。材料设置为矿体和充填体。
图 3~6为4种模型在同等条件下Y方向上的位移分布图。其最大值分别为0.098、0.089、0.111、0.099m，差别不是很大。图 7~10为4种模型在同等约束条件下的等效应力分布图。
通过数值模拟并对结果进行分析，应力最为集中的地方处于两步进路开采时靠近矿体一侧的顶板上，具体对4种采场结构参数模型位移和应力情况比较见表 3。
表 3 Ⅲ盘区模型位移、应力模拟结果
从表 3可以看出，模型P1和模型P3无论从最大位移还是从最大等效应力方面，都明显比其他两个模型数值上要大些，从这个方面来讲，模型P2和模型P4较模型P1和模型P3要优。通过对模型P2和模型P4进行比较后发现，二者在位移和应力集中方面近似相同，但是结合矿山实际情况，对产量要求较高，故模型P4要较模型一更为切合实际，这也与矿山现阶段的实际情况相符（试验采场采用5m×7m的进路规格）。
3.2 矿体整体模型计算
矿体整体模拟的主要目的是对Ⅰ、Ⅲ盘区超前Ⅱ、Ⅳ盘区的层数进行研究，根据经验，往往是Ⅰ、Ⅲ盘区先进行回采，Ⅱ、Ⅳ盘区紧随其后，便于增大生产能力，但是，超前的层数不同，Ⅱ、Ⅳ盘区进路回采时应力集中状态就不同，根据上文所计算出的最佳进路规格，选取进路为5m×7m规格。边界条件为：矿体两侧水平方向位移为零，矿体底面水平方向和竖直方向位移均为零，矿体承受13.2MPa上覆围岩压力。材料设置为矿体和充填体。
通过数值模拟并对结果分析处理后得知，矿体整体模型在4个盘区同时回采时，其应力主要集中在一步正在开采进路的顶板和侧帮处，随着Ⅰ、Ⅲ盘区比Ⅱ、Ⅳ盘区超前的层数不同，应力集中的程度也有所不同。具体数据如表 4所示。
表 4 矿区整体模型位移、应力模拟结果
由表 4数据可知，方案一除了在Y方向上的位移略小于方案二，其余几个方面均比方案二的数值大，应力较集中，故选择方案二，即Ⅰ、Ⅲ盘区比Ⅱ、Ⅳ盘区超前5层，也就是25m。
4 结论
通过对新城金矿极破碎难采矿体使用有限元程序ANSYS软件进行数值模拟、计算分析并参照现场经验，得出如下结论：
（1）采场的应力分布与其结构参数密切相关，合理的进路参数可以明显地改善应力状态。
（2）通过对单个盘区开采现状的模拟，采用4种不同的进路参数进行计算，得出最佳进路规格，与现场试验采场采用的参数基本相同，验证了数值模拟的正确性。
（3）建立了矿体整体模型，在得出最佳进路规格的基础上，采用了两种超前高度进行模拟计算，最终得出合理的超前高度为25m。
（4）数值模拟作为一种特殊意义下的实验，是科技发展的潮流，不仅减少了试验工作量，更对实际工程具有指导意义。
-
[1] 董卫军, 苏永华, 孙玉堂, 等.乳山金矿深部采场结构参数的数值模拟[J].采矿技术, 2002, (6):15-18. http://www.docin.com/p-7438146.html [2] 乔兰.新城金矿深部采场结构参数和开采顺序优化研究[J].金属矿山, 2001, (6):11-15. doi: 10.3321/j.issn:1001-1250.2001.06.004
下载:
计量
- 文章访问数: 68
- HTML全文浏览量: 45
- PDF下载量: 2
