Room and pillar sizes are key factors for safe mining and ore recovery in open-stope mining. To investigate the influence of room and pillar configurations on stope stability in highly fractured and weakened areas, an orthogonal design with two factors, three levels and nine runs was proposed, followed by three-dimensional numerical simulation using ANSYS and FLAC3~. Results show that surface settlement after excavation is concentrically ringed, and increases with the decrease of pillar width and distances to stope gobs. In the meantime, the ore-control fault at the ore-rock boundary and the fractured argillaceous dolomite with intercalated slate at the hanging wall deteriorate the roof settlement. Additionally, stope stability is challenged due to pillar rheological yield and stress concentration, and both are induced by redistribution of stress and plastic zones after mining. Following an objective function and a constraint function, room and pillar configuration with widths of 14 m and 16 m, respectively, is presented as the optimization for improving the ore recovery rate while maintaining a safe working environment.
Methods and procedures of three-dimensional （3D） characterization of the pore structure features in the packed ore particle bed are focused. X-ray computed tomography was applied to deriving the cross-sectional images of specimens with single particle size of 1-2, 2-3, 3-4, 4-5, 5-6, 6-7, 7-8, 8-9, 9-10 ram. Based on the in-house developed 3D image analysis programs using Matlab, the volume porosity, pore size distribution and degree of connectivity were calculated and analyzed in detail. The results indicate that the volume porosity, the mean diameter of pores and the effective pore size （d50） increase with the increasing of particle size. Lognormal distribution or Gauss distribution is mostly suitable to model the pore size distribution. The degree of connectivity investigated on the basis of cluster-labeling algorithm also increases with increasing the particle size approximately.
The images of granular ore media were captured by X-ray CT scanner. Combined with digital image processing and finite element techniques, the three-dimensional geometrical model, which represents the realistic pore structure of the media, was constructed. With this model, three dimensional pore scale fluid flow among particles was simulated. Then the distributions of fluid flow velocity and pressure were analyzed and the hydraulic conductivity was calculated. The simulation results indicate the fluid flow behaviors are mainly dominated by the volume and topological structure of pore space. There exist obvious preferential flow and leaching blind zones simultaneously in the medium. The highest velocities generally occur in those narrow pores with high pressure drops. The hydraulic conductivity obtained by simulation is the same order of magnitude as the laboratory test result, which denotes the validity of the model. The pore-scale and macro-scale are combined and the established geometrical model can be used for the simulations of other phenomena during heap leaching process.
This paper dealt with the development of a two-dimensional （2D） mathematical model for column leaching and confirmed the important simulation parameters through experiment. The unsaturated state of the variably saturated flow column and the solute transport of copper ions were studied during leaching. The fluid flow problem was handled using the Richards equation on the premise of an ambient pressure column air, where the van Genuchten formulas were applied to define the nonlinear relationships of pressure head with the retention and permeability properties. The ore column permeability test gave a varied hydraulic conductivity, which was analyzed in the model. In the solute transport problem, the copper ion concentration was solved using the advection-diffusion-reaction equation whose reaction term was determined by the joint analysis of experimental copper leaching rate and the shrinking core model. Particle-and column-scale leaching tests were carried out to illustrate the difference and connection of copper extraction in both processes. This fluid flow and solute transport cou-pled model was determined through the finite element method using the numerical simulation software, COMSOL Multiphysics.
Xiu-xiu MiaoAi-xiang WuBao-hua YangJin-zhi LiuSheng-hua YinHong-jiang Wang
