Search Documents

By B King

A large number of mineral deposits contain regions in which it is economically beneficial to extract ore both via open pit and underground mining techniques. Current industry practice lacks an effective means to establish the extent of the open pit and underground mining areas for an individual deposit. We present a mixed integer programming model for scheduling both the open pit and underground mine that can be used to create a holistic extraction sequence for a specific deposit. The model considers both spatial and temporal elements of the open pit and underground mines in order to maximise the net present value of the entire deposit, while adhering to its production, processing and block sequencing constraints. By doing so, we can schedule the open pit and underground in concert.CITATION:King, B, 2014. Modelling open pit and underground production scheduling in concert, in Proceedings Orebody Modelling and Strategic Mine Planning Symposium 2014 , pp 369–374 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Nov 24, 2014

By A. Krause

Shovel-truck systems for loading and hauling material in open pit mines are now routinely analysed using simulation models or off-the-shelf simulation software packages, which can be very expensive for once-off or occasional use. The simulation models invariably produce different estimations of fleet sizes due to their differing estimations of cycle time. No single model or package can accurately estimate the required fleet size because the fleet operating parameters are characteristically random and dynamic. In order to improve confidence in sizing the fleet for a mining project, at least two estimation models should be used. This paper demonstrates that the Machine Repair Model can be modified and used as a model for estimating truck fleet size in an open pit shovel-truck system. The modified Machine Repair Model is first applied to a virtual open pit mine case study. The results compare favourably to output from other estimation models using the same input parameters for the virtual mine. The modified Machine Repair Model is further applied to an existing open pit coal operation, the Kwagga Section of Optimum Colliery as a case study. Again the results confirm those obtained from the virtual mine case study. It is concluded that the Machine Repair Model can be an affordable model compared to off-the-shelf generic software because it is easily modelled in Microsoft Excel, a software platform that most mines already use. This paper reports part of the work of a MSc research study submitted to the University of Witwatersrand, Johannesburg, South Africa. Keywords: simulation, bunching, probability distributions, cycle time, queuing, matching, shovel-truck, OEM.

Jan 1, 2007

By R Shaw, M Hidding

The design of primary sample cutters for high-capacity conveyors is critical to achieving a correct sample, both in collecting the sample and ensuring smooth sample discharge flow. Modelling usually occurs at the design stage and requires extensive experience to achieve the correct result. However, it is difficult to confirm the effectiveness of the cutter design until after installation. Even with high-speed video, confirmation of complete sample capture is difficult, especially if the material characteristics change over time. Recent investigation using discrete element method (DEM) simulation to model the flow of the ore into the cutter and observe how it interacts with the cutter has revealed some interesting design requirements not observed in the past. By including DEM modelling, the correct design and optimisation of the cutter can be more accurately determined prior to commissioning.CITATION:Hidding, M and Shaw, R, 2017. Modelling ore flow in the design of high-capacity sample cutters, in Proceedings Eighth World Conference on Sampling and Blending , pp 145–150 (The Australasian Institute of Mining and Metallurgy: Melbourne).

May 9, 2017

By D. Xie, S. Sun

An overview of recent development of CSIRO?s Multi-Phase reaction models will be presented with examples comparing the calculated physico-chemical properties of metallurgical slags with the existing experimental data. The paper will also outline the development of a kinetically based model for molten bath reactors and its application to submerged combustion of coal in a slag bath as well as the reaction between the slag and injected coal particles. The plume model combines established models for the combustion of pulverized fuel with models for the growth and rise of bubbles from a lance or tuyere. This model predicts the extent of de-volatilization and combustion of the coal in the plume, the utilization of oxygen and the exchange rate of oxygen between the gas and the bath as well as the heat transfer rate between the gas bubbles and the bath. The slag reduction model brings together the recently measured rate data and fundamental understanding of the kinetics of CO-CO2reaction with iron oxide containing slags and predicts the changes resulting from injection of coal particles into molten slags. The top-space model covers reactions between gases and molten droplets in the post-combustion zone, where heat and mass transfer between phases results in higher utilization of the fuel values in the partly oxidized gases and transfer of heat to the molten bath. Keywords: Phase equilibria, viscosity, electrical conductivity, slags, kinetics, reaction rate.

Jan 1, 2004

By E. Tabosa, S. Vianna, K. Duffy, W. Valery, B. Andrade, L. Pyle, P. Holtham

Pneumatic flotation cells are high-intensity flotation machines with rapid kinetics, resulting in short residence times. Examples include the Jameson, PneuFlot and ImhoFlot cells. These cells do not use impeller-like mechanical cells. Rather, the air and pulp are mixed in a continuous stream of high fluid velocity through a downcomer or venturi to disperse the air into fine (approximately 0.5 mm) bubbles and maintain particle suspension. Froth wash water is often used to enhance concentrate grade, minimising recovery by entrainment. With the growing application of these types of cells, not only in cleaner stages (their traditional use) but also in pre-flotation and scalping duties, there is a need to predict their performance more accurately for circuit design and optimisation. No widely accepted model for pneumatic flotation cells is available and their metallurgical performance is often scaled up from laboratory and/or pilot plant testing by assuming fixed values for metal recovery, concentrate grade and concentrate percent solids. This paper describes a new approach to simulate pneumatic flotation cell performance and presents case studies demonstrating its applicability. Keywords: Pneumatic flotation, flotation modelling, flotation circuit design, flotation circuit optimisation, Jameson cell

Jan 1, 2020

By Boland J

Prediction of rock cuttability by various excavation machines such as tunnel boring machines, road headers, rippers, picks, roller disc cutters, excavators, etc is of major concern to both civil and mining engineers and rock mechanics. In principle, the design of a rock cutting tool is based on achieving maximum excavation with minimum input energy, minimum applied force, minimum equipment weight and volume and, more importantly, minimum tool wear. When a rock material is indented by the force of a cutting tool and then unloaded several material zones form underneath and around the tool. As the tool moves into the rock the volume of its associated rock cavity (ie volume of displaced rock) expands and depending on the indenter's shape various phenomena may occur. Immediately beneath the tool a small zone in which the stress state is more or less hydrostatic, is formed which is underlain by an inelastic, or damaged, material zone. Both size and shape of this damaged zone control the magnitude of the indentation force and the chip formation mechanism during unloading. The rock is unable to store the whole strain energy induced by the tool and consequently the extra energy is expended either in forming new fracture surfaces and/or in plastic deformation in the damaged zone. Therefore, it is not surprising to see that both perfect-elastic and rigid-plastic constitutive models are not generally suitable for modelling indentation of rocks. In contrast elastoplastic constitutive models are more appropriate as they can include both elastic and plastic deformations (Kral et al, 1993). Elastic solutions for several half-space indentation problems can be found in most text books on elasticity and in particular in Timoshenko and Goodier (1951), Sneddon, (1951), Poulos and Davis (1974) and Johnson (1985a). If a failure criterion, eg the Mohr-Coulomb failure criterion, is applied to delimit elastic solutions for indentation problems, a lower bound limit of the applied indentation force (for a given indentation depth) can be estimated. It may be noted that elastic analyses, without any delimiting failure criteria, normally overestimate the indentation force. Plastic solutions for a few indentation problems may be found in Hill (1950), Johnson, (1985a) and Alehossein et al (1992). When analysing the behaviour of indented rocks the choice of an appropriate constitutive model describing effectively the behaviour of the rock at various stages and phases during indentation is a challenging issue. Several questions need to be answered before starting an effective analysis. The first question is of course: Which one of these is more appropriate for modelling zones underneath and around a cutter, elasticity or plasticity theory, continuum damage mechanics or fracture mechanics, and what is the bridging link between these various theories? Which mode is more dominant in fracturing rock, tensile or shear mode? What governs the mechanism of chip formation? What is the function, shape and the behaviour of the pulverised zone that forms just underneath a cutter? Is the hydrostatic stress state assumption reasonable for this zone? How does the rock change during the absorption and dissipation of the input energy and how does the indentation energy vary with rock type and depth of penetration? How can we model the damaged zone under the indenter? Furthermore does a plastic material form beneath the pulverised zone? When and how do radial and ring cracks form around and under an indenter?

Jan 1, 1995

By Pierre Montmitonnet, Yukio Shigaki, Jonatas Mezêncio Silva

A precise roll stack deformation model of a flat rolling mill is crucial for strip profile assessment as it depends strongly on the roll’s deformation. The need for more and more advanced products with stringent both geometrical and mechanical tolerances (e.g., Advanced High Strength Steels – AHSS), with thinner and harder steel strips, demands more advanced simulation models. The present paper addresses assumptions made in Finite Element Method (FEM)-based roll stack models on some important, yet barely discussed features: What is the impact on precision of the simplifications such as symmetry and anti-symmetry? Should roll neck boundary conditions in bearings be considered as fixed or simply supported? In this paper a hybrid model developed by the authors for a general flat rolling mill is used to assess those questions. This model runs much faster (100x in some cases) than a whole-FEM-modeled 6-high rolling mill with the elastoplastic strip deformation, thanks to its simplifications adopted. This model assumes a steady state rolling condition, and the roll stack is modeled with tridimensional finite elements, and the loads are calculated via a slab method adopting a plane strain state for the strip deformation (especially true for thin strips). An industrial 6-high cold rolling mill is modeled and results on the application of symmetries and different boundary conditions on the roll’s necks are presented. It is concluded that the application of symmetries and anti-symmetries for the 6-high rolling mill is possible for thin strips. And the way the necks are modeled does not impact on the final strip shape, though it influences the magnitude of the roll’s displacement. The coupling of the strip with the work roll applied assuming a symmetry plane of the arc of contact has little influence on the strip shape.

Oct 1, 2019

By J. A. Meech, R. J. Tucker

"A dynamic computer model has been developed to simulate the operation of ore crushing plants. The model was created in a general form to allow imitation of a wide variety of secondary crushing circuits with only minor program modifications. The system can be run on either a microcomputer or a large time-sharing main frame computer as required.The size reduction and separation characteristics of crushers and screens are based on standard manufacturers' design and performance data. These calculations are isolated within various subroutines to facilitate tailoring of the model to emulate any particular plant of interest. Input and output data are presented in a variety of formats. Crusher, screen and conveyor amperages and surge bin levels can be presented in a tabular form along with appropriate alarms arising from improper circuit operation or process interruptions. Graphical representation of the circuit together with strip chart data output are features of the system when operated in a microcomputer environment.The program is currently being used at Queen's University as a teaching aid to demonstrate plant operating characteristics, control strategies and design concepts. The acceptance of, and reaction to the program at this level indicates a good potential for it to provide a valuable tool for operator training."

Jan 1, 1985

By Anna Franzkowiak, Dean Gregurek, Alfred Spanring, Angelika Ressler

In order to determine the most suitable refractory products and improve the lining lifetime for the diverse furnaces used in the nonferrous metal industry, corrosion tests are performed at RHI’s Technology Center. The practical facilities include the cup test, induction furnace, rotary kiln, and drip slag test described in this paper, which enable a comprehensive understanding of the chemo-thermal brick wear on a pilot scale. The corrosion trials are performed with actual slags generated during operations at a customer’s plant. To determine the highest influencing wear parameter, every single test is combined with a detailed mineralogical investigation and thermochemical calculations performed using FactSage. Based on the results, tailored refractory solutions for the nonferrous metal industry can be provided in combination with trials conducted at the customer’s site.

Jan 1, 2014

By Q. H. Liao, B. G. Monaghan

The options for designing slopes in complex geological conditions are reviewed. It is shown that modelling, especially numerical modelling, allows a degree of analysis that cannot otherwise be achieved. Examples are presented of how numerical modelling can be used to investigate fundamental aspects of rock mass behaviour. It is demonstrated that an important consequence of change in scale is that the mechanisms of failure also change. The complex failure of the footwall slope at Aznalcollar, a large open-pit mine in southern Spain, is discussed. Numerical modelling with UDEC (universal distinct-element code) has allowed the realistic simulation of displacements measured in the field

Jun 18, 1905

By D Craw, P Upton

"An EXTENDED ABSTRACT is available for download. A full-length paper was not prepared for this presentation. The Southern Alps of New Zealand and the mountains of Taiwan are both young transpressive orogens characterised by rapid uplift and exhumation, high heat flow and vigorous surface processes. Comparing and contrasting the two provides an ideal framework through which to explore the structural and lithological controls on fluid and gold mobility during collisional deformation. Both have two fluid flow systems centred about their drainage divide with shallow topographically-driven meteoric water above a deeper flow system dominated by mineralising rock-exchanged fluids. Despite similar tectonic settings, the distribution of heat flow, hot springs and veins in the two orogens are different. Taiwan has higher heat flow, distributed hot springs and localised veining whereas the Southern Alps have a narrow heat flow anomaly, localised warm springs and widespread veining. We use three-dimensional numerical models to evaluate key parameters that affect fluid flow and gold mineralisation. Models show that differences in crustal strength and thickness affect the locus, scale and rate of crustal fluid flow and hence influence gold deposition at a variety of scales. We explore how surface topography can lead to localised changes in the stress state, facilitating the switching of principal stresses and hence vein orientations.CITATION:Upton, P and Craw, D, 2015. Modelling structural and lithological controls on the mobility of fluids and gold in orogenic belts – examples from New Zealand and Taiwan, in Proceedings PACRIM 2015 Congress, pp 189–192 (The Australasian Institute of Mining and Metallurgy: Melbourne)."

Mar 18, 2015

By Mike Dry

While McCabe-Thiele plots are useful for designing single-component solvent extraction systems, not least because of their simplicity, they are less useful for multi-component solvent extraction. This paper presents a relatively simple method, compared to other approaches, for modelling solvent extraction as chemical equilibrium reactions between the various aqueous and organic species. The essence of this approach is how values for the equilibrium constants can be extracted from standard laboratory data. Synergistic solvent extraction of nickel is used to illustrate the method.

By F. Bello, R. Yanez, M. Pinto, T. Errazuriz

"Throughout its operation, Minera Esperanza (ME) will deposit one of the largest tonnages rates of high density tailings in the world with a solids concentration of close to 67% by weight. The tailings presently being deposited have a solids concentration of approximately 62% - 64%, and so cannot be used to directly characterize the flow and deposition characteristics of the higher concentration tailings. To increase the knowledge of the behaviour of their tailings, ME, with the support of Golder Associates, developed a large scale tailings deposition test with a capacity of 30,000 tonnes per day (1,200 m3/h). The test was carried out by depositing tailings into two large cells or paddocks in five stages. Multiple spigots were used to allow different flows to be deposited. At each stage, topographical surveys of the natural terrain and the geometry of the deposited tailings, as well as continuous flow measurement and periodical rheological characterization of the discharged tailings were conducted to create a data set against which beach slope models could be calibrated and then used to predict performance of the higher concentration tailings under different deposition scenarios. This study adopted the beach slope prediction model proposed by A.L. Li (Li, 2011). The model was successfully adjusted to the observations made during this large scale field test by applying the values measured in-situ. This paper presents the application of the model, the results obtained, the findings, and the analysis of possible future scenarios. The model was applied, using the discharge flow rate, the geometrical characteristics of the discharge area, and the rheology with its variation along the beach to estimate the development of the beach formation and the steady state of the beach slopes. Through the use of dimensionless variables, the model was applied to higher solids concentration scenarios, obtaining results that would allow beach slope estimation at different deposition stages in multiple scenarios of varying parameters such as flow, rheology and geometry of discharge area."

Jan 1, 2014

By A. V. Rose

A redevelopment in the City of London has been constructed with a piled foundation system that is believed to be the first of its kind. The purpose of the research is to explore this novel foundation method, which involves the installation of closely spaced piles in a perimeter arrangement with a body of soil enclosed in the centre. Both centrifuge modelling and numerical modelling techniques are used.

Jan 1, 2011

A model has been developed as part of the AMIRA P420 Gold Processing Technology Project to simulate the recovery of gold by centrifugal gravity devices. The model (which is available online) uses the simulation methodology developed by Laplante et al (1995) to predict gravity recoverable gold (GRG) recovery under varying operating conditions and circuit configurations. In addition to predicting GRG recovery, the model also predicts the GRG content and size distribution of the cyclone overflow and underflow. The predicted overflow GRG content and size distribution can be used as inputs for a leaching model in which size dependent leaching rates are applied for each size class of GRG and non-gravity recoverable gold (non-GRG) in the cyclone overflow, allowing size-by-size leach recoveries to be determined for GRG and non-GRG. The application of the gravity model and leaching model are demonstrated using a case study. The simulated gravity recovery and leaching performance showed excellent agreement with the measured values. The developed models can be used to assess the sensitivity of gravity gold and leach recoveries to varying feed GRG content, GRG size distribution, gravity circuit feed rate, classifier performance and leach residence time. The described methods can be used to optimise new and existing gravity and leaching circuits.

Aug 1, 2010

By V. G. Papangelakis, M. Cheung, M. Dry

The high pressure acid leaching of laterite was modelled and simulated using a combination of two commercially available software packages, namely Aspen Plus 11.1 and the OLI Stream Analyzer. New kinetic parameters for the acid-controlled metal dissolution and precipitation reactions at temperature were first determined from previous experimental data on the leaching of limonite. A grain model with cylindrical pellets was used to model the solid-liquid dissolution reactions. Precipitation reactions were modeled via aqueous equilibrium chemistry. A dedicated thermodynamic database previously developed at the University of Toronto was used in conjunction with the OLI software to model the solution chemistry and calculate the concentration of hydrogen ions in solution at the leaching temperature. Mass balances were handled by Aspen Plus. The simulation gave good predictions of nickel dissolution under industrial and laboratory conditions, for the dissolution of limonite and limonite/saprolite blends. Computational

Jan 1, 2005

By P Holtham, R Dunne, D Stoll

Flash flotation is a process designed to remove coarse, liberated, fast floating minerals from hydrocyclone underflow. Irrespective of liberation, coarse particles would normally remain in a grinding circuit until they are reground fine enough to report to the hydrocyclone overflow. Additionally, minerals of high density also report to hydrocyclone underflow despite being fine in size and/or liberated because hydrocyclones separate on the characteristics of both particle size and density. Dense particles continually reporting to the hydrocyclone underflow can lead to over-grinding and the production of ultra fine particles. These particles exhibit slow flotation kinetics and thus are difficult to float in conventional flotation circuits. Therefore, removing coarse, liberated particles from the hydrocylcone underflow by flash flotation has the potential to improve recovery and can improve dewatering unit operations and reduce grinding costs.The impact on the overall performance of the flotation circuit by removing coarse, fast floating minerals prior to conventional flotation is not well understood. Ultimately, predicting the overall recovery improvement provides the justification for the installation of a flash flotation cell within the milling circuit. This paper will describe attempts to estimate flash flotation recovery improvements to the overall circuit using kinetic data from simplified batch cell flotation tests combined with limn-based modelling techniques.CITATION:Stoll, D, Dunne, R and Holtham, P, 2012. Modelling the impact of flash flotation from batch cell flotation tests using limn, in Proceedings 11th AusIMM Mill OperatorsÆ Conference, pp 179-184 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Oct 29, 2012

By S S. Kanchibotla, I Brunton

Most blast fragmentation models assume the rock mass properties, explosive properties and blast design variables to be constants and uniformly distributed within a blast. However, in reality all these input variables vary within a blast resulting in variation in the resulting fragmentation size distribution. A stochastic modelling approach is introduced in this paper to quantify this variation. This technique takes the input variables as statistical distributions rather than constants and through several thousand iterations, generates a statistical representation of the expected fragmentation resulting from a production blast. A case study of three production blasts from a large open pit mine are presented and the modelled fragmentation æenvelopeÆ shows good agreement with the fragmentation æenvelopeÆ estimated from Split image analysis. The various blast-related parameters influence different parts of the fragmentation distribution, eg rock strength and explosive velocity of detonation have most impact on the fines. The technique is used to identify the parameters that have the greatest influence on various size fractions. Such an analysis will be useful to direct resources to efficiently minimise the variation.

Jan 1, 2001

By Markus A. Reuter, Yongxiang Yang, James R. Post

The packed coke bed in the lower part of an iron making blast furnace, the dead man, dissolves slowly in a pool of hot metal, which is tapped regularly. The coke particles of the dead man not only exhibit a size distribution, but also their size distribution is non-uniform throughout the dead man and as a result, the porosity of the dead man is non-uniform. This causes the flow of hot metal through the dead man to be irregular and non-uniform, bringing about channelling in the coke bed. Channelling is detrimental to the refractory lining life of the hearth of a blast furnace. Using population balance modelling and a porosity model it is possible to keep track of the local changes in size distribution and bed porosity during dissolution of coke particles. Numerical simulations show that carbon dissolution from a packed coke bed (with particles which exhibit a size distribution) increases the average porosity, decreases the mean particle size and changes the size distribution to a distribution with a smaller size range. Influence of the type of initial size distribution (mean, variance) is investigated. The calculated changes in porosity are used to dynamically adapt the flow and the mass transfer (kinetics) during calculations. This approach provides the possibility to implement dynamic porosity calculations for CFD applications, and thus simulate the results of various coke types and particle size distributions on dissolution and flow behaviour.

Jan 1, 2003

By R Wolff

One critical factor in open pit mining projects is the estimation of the recoverable reserves. The reason for this is that at the valuation stage there is a lack of information about future metal prices and production costs. Consequently, the estimation of the recoverable reserves needs to be done based on a fixed break-even cut-off grade (BECG), in which fixed expected metal prices and production costs are assumed throughout the operating life of the mine (OLM). In this paper, an alternative technique for estimating the recoverable reserves of an open pit mining project is presented and explained in detail. The main characteristic of this technique is that it estimates the recoverable reserves of the project using mining costs as reference: this is done due to the fact that mining costs are easier to model than metal prices. Another important characteristic of this technique is that it uses the uncertainty of the geology of the orebody, given by simulations of the deposit, to make final strategic decisions, such as the selection of the operating life of the mine and the optimal production rate, which minimise the risk of not achieving future production targets while maximising the upside potential of future rewards. As it will be shown, the final outcome of the proposed technique is the generation of a robust open pit design, that include cutbacks and ultimate pit limits, called the marginal open pit mine design (MOPMD). One characteristic of the MOPMD is that, at each production period, it will be characterised by their respective geological uncertainty, given as probability distributions of ore, waste and metal quantities, and by marginal economical project indicators, such as cut-off grades and metal prices, among others. The details follow in the paper.

Jan 1, 2006