刘洋等：基于主要目标法的堑沟底部结构位置优化 ·317· 水平时，除堑沟巷道与矿体间距为18.35m方案，其他 采出金属量和采出品位先上升后下降，但两者不同时 方案均不满足采出品位约束条件：当堑沟巷道所在水 达到极大值，当堑沟巷道与矿体间距过小时不满足采 平一定时，随着在水平方向从矿体边界向矿体侧移动， 出品位约束条件，且采出金属量下降至较低水平. 表1以1m步距运算结果表 由运算结果可以看出：在105、104和103m水平 Table 1 Computation results by Im per step 满足采出品位约束条件的对应L值分别为21.35~ H/m 17.35、20.35~17.35和18.35m.由于104m水平金属 回采量普遍大于105m水平且根据品位确定的可行 L/m 105 104 103 方案范围差别不大，所以在垂直方向上确定运算范 0/t G/% QIt G1% QIt G/% 围为104~103m;在水平方向上，远离矿体侧采出金 22.352741.4517.512809.4916.152850.0415.04 属量较大，为避免遗漏最优方案，选择21.35m边界： 21.352782.3318.422870.4517.192924.3915.86 靠近矿体侧采出金属量已较小，选择17.35m为边界 20.352768.1319.112905.1018.112986.7016.73 即可.根据上述条件，修改约束条件的最优化问题描 19.352643.6419.392870.5818.703017.05 17.60 述如下： 18.352399.2219.072716.2318.672961.0918.06 maxQ=f(L,Hl）, 17.352101.5018.402453.6518.232777.5317.94 s.t.17.35≤L≤21.35， (6) 16.351833.8717.752161.3817.572504.0917.43 103≤H≤104， 15.351594.0217.131910.6216.992227.4616.81 G≥18%. 14.351387.4416.611682.0316.421989.8116.26 SURPAC软件通过运算4个顶点的小四面体的体 13.351250.0216.521498.0916.041775.0415.76 积进行累加获得实体体积，且任意两顶点间距不小于 12.351161.5716.771381.9316.081618.7215.53 0.05m,可进行精度不高于0.05m的运算.对新的约 11.351115.0917.181311.3416.391526.2615.68 束条件，将移动步距修改为0.1m进行精细运算，计算 10.351106.4617.671284.4716.811474.4516.01 结果见图5 根据图5(b)可知满足采出品位约束条件的可行 a 3000 2900 2800 2700 2600Y 2500 2400 21.1 2 20.75 0. 104.0 102R /103.6 m 18.75 11251030 =2400-2500=2500-2600■2600-27002700-2800■2800-2900=2900-3000 0.19 0.18 毫017 0.16 0.15 20.9 0. 2055 . 7104.0 20. 19.95 .75 7103.8 /103.6 1/m 185 8.3 8.15 17.7 3510a0 55 ■0.18-0.19■0.17-0.18■0.16-0.17■0.15-0.16 图5采出金属量()及采出品位(b)与堑沟底部结构空间位置关系分布图 Fig.5 Scattergrams of trench position vs.mineout metal amount (a)and trench position vs.mine-out grade (b)刘 洋等: 基于主要目标法的堑沟底部结构位置优化 水平时，除堑沟巷道与矿体间距为 18. 35 m 方案，其他 方案均不满足采出品位约束条件; 当堑沟巷道所在水 图 5 采出金属量( a) 及采出品位( b) 与堑沟底部结构空间位置关系分布图 Fig． 5 Scattergrams of trench position vs． mine-out metal amount ( a) and trench position vs． mine-out grade ( b) 平一定时，随着在水平方向从矿体边界向矿体侧移动， 表 1 以 1 m 步距运算结果表 Table 1 Computation results by 1 m per step L /m H/m 105 104 103 Q/t G/% Q/t G /% Q/t G /% 22. 35 2741. 45 17. 51 2809. 49 16. 15 2850. 04 15. 04 21. 35 2782. 33 18. 42 2870. 45 17. 19 2924. 39 15. 86 20. 35 2768. 13 19. 11 2905. 10 18. 11 2986. 70 16. 73 19. 35 2643. 64 19. 39 2870. 58 18. 70 3017. 05 17. 60 18. 35 2399. 22 19. 07 2716. 23 18. 67 2961. 09 18. 06 17. 35 2101. 50 18. 40 2453. 65 18. 23 2777. 53 17. 94 16. 35 1833. 87 17. 75 2161. 38 17. 57 2504. 09 17. 43 15. 35 1594. 02 17. 13 1910. 62 16. 99 2227. 46 16. 81 14. 35 1387. 44 16. 61 1682. 03 16. 42 1989. 81 16. 26 13. 35 1250. 02 16. 52 1498. 09 16. 04 1775. 04 15. 76 12. 35 1161. 57 16. 77 1381. 93 16. 08 1618. 72 15. 53 11. 35 1115. 09 17. 18 1311. 34 16. 39 1526. 26 15. 68 10. 35 1106. 46 17. 67 1284. 47 16. 81 1474. 45 16. 01 采出金属量和采出品位先上升后下降，但两者不同时 达到极大值，当堑沟巷道与矿体间距过小时不满足采 出品位约束条件，且采出金属量下降至较低水平． 由运算结果可以看出: 在 105、104 和 103 m 水平 满足采出品位约束条件的对应 L 值分别为 21. 35 ～ 17. 35、20. 35 ～ 17. 35 和 18. 35 m． 由于 104 m 水平金属 回采量普遍大于 105 m 水平且根据品位确定的可行 方案范围差别不大，所以在垂直方向上确定运算范 围为 104 ～ 103 m; 在水平方向上，远离矿体侧采出金 属量较大，为避免遗漏最优方案，选择 21. 35 m 边界; 靠近矿体侧采出金属量已较小，选择 17. 35 m 为边界 即可． 根据上述条件，修改约束条件的最优化问题描 述如下: max Q = f1 ( L，H) ， s． t． 17. 35≤L≤21. 35， 103≤H≤104， G≥18 { % ． ( 6) SUＲPAC 软件通过运算 4 个顶点的小四面体的体 积进行累加获得实体体积，且任意两顶点间距不小于 0. 05 m，可进行精度不高于 0. 05 m 的运算． 对新的约 束条件，将移动步距修改为 0. 1 m 进行精细运算，计算 结果见图 5． 根据图 5( b) 可知满足采出品位约束条件的可行 ·317·