Search Documents

By H. P. Kritzinger, T. C. Kingsley

Rotary kilns are used for a variety of mineral processing operations. Hatch makes use of a kiln model developed from first principles to evaluate designs for its clients. This tool has been applied to a variety of pyrometallurgical applications. This paper illustrates the application of numerical optimization techniques in combination with the kiln model in the interrogation of a generic iron-ore reduction process. The fundamental modelling concepts are explained, followed by a description of the optimization approach. The results show how the combination of the two methods, with modern computing power, can generate a large number of viable design and operating candidates.

Jan 1, 2014

By H. P. Kritzinger, T. C. Kingsle

"Rotary kilns are used for a variety of mineral processing operations. Hatch makes use of a kiln model developed from first principles to evaluate designs for its clients. This tool has been applied to a variety of pyrometallurgical applications, including ferrovanadium, nickel carbonate, nickel laterite, iron ore reduction, and spodumene (lithium) production.This paper illustrates the application of numerical optimization techniques in combination with the kiln model in the interrogation of a generic iron ore reduction process. The fundamental modelling concepts are explained, followed by a description of the optimization approach. The results show how the combination of the two methods, with modern computing power, can generate a large number of viable design and operating candidates."

Jan 1, 2015

By J. Lehmann, H. Gaye

A brief review of Computational Thermodynamics software concerning the modeling and prediction of reactions involving metals, slags and fluxes is presented. Their application for various steelmaking treatments is discussed. Examples concern phosphorus removal in BOF and EAF, ladle steel desulfurization, and specific treatments for inclusions control in semi-killed or Al-killed carbon and stainless steels.

By Ravi S. Bellur-Ramaswamy, Daniel A. Tortorelli, Robert Haber, Nahil A. Sobh

We optimize continuous quench process parameters to produce functionally graded aluminum alloy extrudates. To perform this task, an optimization problem is defined and solved using a standard nonlinear programming algorithm. Ingredients of this algorithm include 1) the process parameters to be optimized, 2) a cost function: the weighted average of the precipitate number density distribution, 3) constraint functions to limit the temperature gradient (and hence distortion and residual stress) and exit temperature, and 4) their sensitivities with respect to the process parameters. The cost and constraint functions are dependent on the temperature and precipitate size which are obtained by balancing energy to determine the temperature distribution and by using a reaction-rate theory to determine the precipitate particle sizes and their distributions. Both the temperature and the precipitate models are solved via the discontinuous Galerkin finite element method. The energy balance incorporates nonlinear boundary conditions and material properties. The temperature field is then used in the reaction rate model which has as many as 105 degrees of- freedom per finite element node. After computing the temperature and precipitate size distributions we must compute their sensitivities. This seemingly intractable computational task is resolved thanks to the discontinuous Galerkin finite element formulation and the direct differentiation sensitivity method. A three-dimension example is provided to demonstrate the algorithm.

Jan 1, 2001

By S. Zhang, Z. Hua, D. Feng, X. Tong, X. Xie

The separation and beneficiation of low grade collophanite has long been a technical challenge in phosphorous production. It requires the joint efforts of process design, reagents optimisation and operational control. In this study, the statistical response surface method was applied for determining the optimum reagent dosing in the flocculate flotation of a low grade (18.59%) collophanite ore. The dosage of three reagents, including a pH regulator (Na2CO3), depressant (Na2SiO3) and flocculant (modified starch), were input parameters of the modelling design. Two quadratic mathematical models were derived for prediction of grade and recovery based on the statistical design. Statistical analysis of the experimental results validates the reliability of the models. Through further calculation, the optimum dosage of the three reagents were found to be 3,020g/t Na2CO3, 1,990g/t Na2SiO3 and 880g/t modified starch, with which a maximum recovery of 75.81% and P2O5 grade of 30.08% was obtained from flocculate flotation test. Keywords: Process optimisation, modelling, froth flotation, collophanite ore, response surface methodology

Jan 1, 2020

By Ioannis K. Kapageridis

The lignite deposits of the Kozani-Ptolemais-Amyntaio basin in western Macedonia (North-West Greece) present a difficult modelling problem. Each deposit consists of several thin lignite layers ranging from a few centimetres to a few tens of centimeters in thickness, and thin intercalated sterile layers of marly limestones, marls, clays, and sands. These layers are practically impossible to model individually due to fast lateral transition of lignite layers to humus clay and vice versa and the lack of reliable and detailed stratigraphic correlation, leading to compositing methods being applied prior to any interpolation and modelling. The complexity of the deposit is further increased by the presence of a number of faults. This paper presents the modelling and estimation procedures applied during the most recent study of the South Western Field ? a new deposit that is considered for development and lignite production over the next years and will support the operation of a new power station in the area. Special focus is given to drillhole database structuring, fault modelling, stratigraphic correlation, mineable lignite intervals compositing, and resource estimation. The use of specialised mine planning software in all stages of modelling and estimation is thoroughly explained.

Oct 1, 2009

By S. Bellec

A leaching model for gold ore that takes into account the impact of cyanide and oxygen concentrations as well as gold surface exposure on leaching performance is developed. Cyanide consumption is estimated by an empirical model that is a function of ore particles size distribution as well as cyanide, copper and sulphur concentrations. The proposed model has been validated by data from laboratory leaching tests on the Sleeping Giant gold mine ore. Keywords: gold leaching, modelling, exposure surface.

Jan 1, 2009

By G. Ferrara

In process and plant simulation, more reliable prediction of industrial results and evaluation of data for designing purposes may be attained by the combined use of (a) modeling and simulation of ore characteristics and (b) modeling and simulation of integrated plant operations. The systems for modeling and simulation of ore characteristics are discussed as well as the conditions in which their use may be necessary (complex structure of the ore-body and variability of the ore fed to the plant). Two examples, for coal and bauxite run-of-mine (ROM), are reported. The procedure of linking the process unit models and the approaches for the solution of integrated systems are indicated. A general-purpose program for the simulation of integrated plant operations is described, and some examples of its use are reported.

Jan 1, 1984

Modelling and Simulation to Design a Flotation Circuit for Gold Concentration

Sep 13, 2010

By S Hancock, B H. Brady, H Muhlhaus

The aims of this paper are to present the governing science, verify the performance and demonstrate the potential engineering applications of a computational simulator for in situ leach mining of non-ferrous minerals. Assuming that the target mineral is deposited either uniformly within rock matrixes (matrix-hosted) or on fracture walls (fracture-hosted), in situ leaching can be treated as reactive flow and transport in geological media (porous media for matrix-hosted ore deposits and fractured media for fracture-hosted ore deposits). Under this treatment, the processes of in situ leach mining for non-ferrous minerals can be quantified through several sets of equations and these coupled equations have been numerically solved. The computational simulator has been applied to the simulations of one-, two-, and three-dimensional in situ leaching cases. The computational simulator, together with laboratory and in situ leaching tests, can be used to:optimise the leaching chemistry; optimise the in situ scale mining operations; scale up the in situ scale test operation to a full-scale mining operation; and evaluate the post-mining environmental restoration plan.

Jan 1, 2000

By O Handayani, R Barr, E Isaaks

The mining industry has focused a great deal of attention recently on the impact of blast movement on grade control. Although the problem of blast movement and its contribution to dilution and ore loss is well known, solutions to the problem are less than satisfactory. Perhaps the most recent innovation is the application of blast movement monitoring (BMM) devices to quantify blast movement. These devices emit an electronic signal and are planted at surveyed locations within the preblast material. Following the blast, their new locations are found using an electronic scanner. The magnitude and direction of the blast movement is given by the vectors between the pre- and post-blast BMM locations. However, it’s not clear how the BMM vector data can be best used to minimise dilution and ore loss. One common solution is to displace preblast dig line polygons by distances indicated by the BMM vectors. However, studies of BMM vectors show that blast movement not only displaces material laterally but also mixes material internally within the bench (Thornton, 2009a).This paper presents a method for modelling the post-blast muck pile that accounts for displacement and internal dilution using simulated annealing. The preblast ore control block model is re-blocked to 1?×?1?×?1?m sub-blocks. Each sub-block is displaced by simulated BMM vectors conditional to post-blast surface topography and blast initiation sequence data in addition to BMM vector data. The displaced sub-blocks are aggregated into new ore control model blocks whose grades are calculated from the contained sub-blocks. New dig lines are then designed on the post-blast muck pile model using the new ore control block model grades. A case study based on actual data is provided.CITATION:Isaaks, E, Barr, R and Handayani, O, 2014. Modelling blast movement for grade control, in Proceedings Ninth International Mining Geology Conference 2014 , pp 433–440 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Aug 18, 2014

By E J. Sellers, A Cocker

Movement of the fragmented rock mass is an inevitable consequence of open cut blasting. Uncontrolled, unaccounted movement from blasting commonly causes ore loss and dilution from waste mixing. Therefore, the prediction of the post-blast location of both ore and waste is vital in grade control and of direct interest to mining geologists. Direct measurement techniques have matured significantly over the past decade, but their use remains expensive and labour intensive. Access and safety issues are further disincentives to regular movement monitoring. What is needed is a combination of direct measurement to establish baselines and reliable predictive modelling.The W H Bryan Mining and Geology Centre of the University of Queensland has pioneered 3D blast movement modelling using ‘physics engines’ to handle in-flight collisions and predict the final resting place of blasted material. Open cut metalliferous mines are a demanding environment for this style of modelling as they can involve complex geometries and variable buffering, which leads to inconsistent confinement. A version of the model designed to deal with these complexities is under development. A well-controlled field trial was needed to validate this new version.A validation trial was completed at a gold mine in South-East Asia. The aim was to increase throughput by increasing the blast powder factor without losing grade through dilution and mixing with waste. Four blocks were blasted to identify the effect of bench height and face confinement on ore movement, fragmentation and process throughput. The powder factor and blast timing were kept constant. Movement of the broken rock was measured by placing blast movement monitors at known positions in the block and then locating them after the blast. The majority of sensors were recovered, providing a unique insight into blast movement dynamics and generating a substantial data set for model validation.The model (called ‘MetaMove ‘) was calibrated for one blast and applied to investigate the extent of grade loss and possible mitigation strategies. Conventional mining techniques that involved firing across the strike of the orebody, and no post-blast boundary adjustment of the ore/waste boundaries prior to digging, were found to cause up to 44 per cent loss of high-grade ore, depending on the grade distributions. Post-blast adjustment of the ore/waste boundaries can negate this loss. The model indicates that altering timing and ensuring that the block boundaries are well confined can further constrain the movement of ore/waste boundaries. Thus, the model can be successfully applied to minimise grade dilution and limit misclassification of ore as waste whilst using higher powder factors to increase mill throughput. Application of a calibrated model to move ore/waste boundaries prior to digging will enable the production geology team to add significant value to the mining operation.CITATION:Cocker, A and Sellers, E J, 2014. Modelling blast movement for grade control at an open cut gold mine, in Proceedings Ninth International Mining Geology Conference 2014 , pp 377–386 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Aug 18, 2014

Tonnage and grade are the principal reported components of mineral resource estimates but tonnage is a derivative of two more fundamental estimates; volume and bulk density. Typically, most of the effort in a resource estimation study is directed towards grade and volume, yet poor application of bulk density data to resource models can result in significant biases. The impact on tonnage estimates is easily recognised but poor density estimates may also introduce subtle biases into the grade estimates. Methods for measuring bulk density have been discussed previously. This paper focuses on the manipulation of bulk density data and the extension of bulk density from essentially point measurements to block estimates within resource models. Bulk density is a continuous, geologically controlled, spatial variable and should be treated as such. Issues such as combining data from multiple sources, sample support and the interdependence of bulk density and grade estimates must be resolved. The statistical and geostatistical techniques used for grade estimation are, under the right circumstances, equally applicable to bulk density. The importance of bulk density needs to be recognised. With robust raw data and correct estimation procedures, better estimates of bulk density can be obtained and both local and global tonnage and grade estimates can be improved.

Jan 1, 2000

Modelling Channel Scale Flow in Heap Leaching Using Smooth Particle Hydrodynamics (SPH)

Sep 13, 2010

Traditionally within the mining industry, single models for both grade and geology of orebodies are created upon which all mine development decisions are based. These models provide a single interpretation of the extent and continuity of the mineralization envelope based on solids and sections interpreted from relatively wide spaced drilling i.e. the models infer behaviour of grade and geology at smaller/local scales based on data available at wider scales. The variable behaviour of grade and geology cannot be derived from single kriged resource models. In recognition of this, Newmont Mining Corporation typically uses Sequential Gaussian Simulation (SGS) to quantify grade variability and to determine its effect on some of these mine development decisions. In most of these studies the mineralized envelope remains fixed and grade simulations are generated within a fixed volume. To account for uncertainty in the geology/ mineralization envelope Newmont Mining Corporation has been experimenting with Multiple Point Statistics (MPS). Multiple-Point Statistics is an emerging technique, which can be used to generate multiple, geologically realistic, realizations of data that displays complex non-linear features. MPS uses a conceptual model of the geology, called a Training Image (TI), to infer these high-order spatial relationships. The multiple realizations allow for an assessment of the geologic uncertainty, which can be further combined with grade simulations to generate a more complete picture of the true uncertainty of a mineral deposit.

Jan 1, 2011

By S. T. Johansen, C. Colomb, B. Ravary

Combustion and thermal NOx formation have been modelled using the Computational Fluid Dynamics (CFD) commercial code Fluent 6.1. The aim was to predict and give a better understanding of the combustion in electric arc furnaces for the production of ferro-silicon and silicon metal. Another important objective was to assess the impact of a number of operating parameters on the NOx formation. The CFD tools used are de-scribed with a special focus on models for species and turbulence transport. Sensitivity studies have been carried out to study the influence of several parameters: sub-models, process gas composition and distribution at the charge surface and influence of the presence of SiO. The main conclusions are: ?Using the combination of a realistic design, boundary conditions and models, the simulated NOx level is in quantitative agreement with measurements. The presence of water vapor is particularly important because of its influence on the specific heat of the mixture and resulting temperature. ?Release of significant amounts of SiO from the charge surface results in an important increase of NOx, even if we applied a simplified combustion model for SiO gas. However, because of lack of thermal and compositional data for blowing situations (when SiO flows directly from the arc cavity through open channels to the charge surface), no quantitative validation is possible. The way SiO is distributed (uniformly or as a jet) has a limited impact on the result in our computations.

Jan 1, 2007

By Catherine Johnson, Barbara Rutter

The shock tunnel located at the Missouri University of Science and Technology Experimental Mine is 65 feet 9 inches (20.0406 m) in length. The tests performed inside of this tunnel are either arena tests or shock interaction and transfer (SIT) tests. Arena tests simulate a large open air blast, with lower explosive weights, by elongating the impulse of the blast wave. SIT tests are conducted to determine how the produced shock wave interacts and transfers its energy into a target object located at the end of the tunnel. The target used in this research is multilayered, with each layer being a different material. The targets are 2 ft. by 4 ft. (0.6096m by 1.2192m) pieces of Medium Density Fiberboard (MDF), DOW polystyrene foam insulation, and gypsum drywall. The order of these materials varies with each test. This shock tunnel is unique, because the geometry changes along its length. To fully understand how the shock waves are affected by these changes, a working model was needed. The second purpose of this model is to confirm the maximum charge size of 400 grams (14.11 oz.) of C4 that can be detonated inside the tunnel without causing damage to or failure of the tunnel, thus reducing the amount of time needed in downtime and repair. The model predicts pressures in various locations both within the tunnel and the attached target, thus assisting in determining where failures would likely occur. In the process of generating this model, several problems were encountered, including filling the tunnel with air, meshing the complex geometry, and creating the tunnel in a CAD system. This paper will discuss how these problems were solved. A working model was created by using ANSYS Explicit Dynamics and ANSYS Autodyn. The tunnel geometry file was created using SolidWorks and includes the complex geometry changes along with the target and its mount located at the end of the tunnel. The results produced by the model were compared to the field SIT tests. The comparison of the results from both the model and the experimental testing showed a similar trend in the pressure traces.

Jan 1, 2018

By A Keddie

Standard microeconomic theory suggests that the behaviour of mineral extraction and processing firms with respect to their demand for water can be modelled and a market demand for water subsequently derived. Such a model, although an abstraction of the real world, should be able to predict the responsiveness of quantity of water demanded to change in price, and, change in demand to a range of factors including: expectations about reliability of supply, changed commodity prices for mineral production, and change in water technology and storage at mine site. This paper demonstrates that the price of water is particularly responsive to changes in demand and that although price increases would result in some additional supplies, physical constraints on supply suggest that the industry will be paying considerably higher prices for water in the future. It is recommended that the mining industry undertakes quantification of possible price increases to position itself to minimise the impact on their operating costs as well as to make responsible capital investment decisions to facilitate future expansion or commissioning of new mines. Substantial investment is required to commission a new mining venture with expected rates of return subject to risk and uncertainty due to fluctuations in world commodity markets as well as reliability of secure supplies of inputs. If a reliable supply of water was subject to uncertainty and water was a critical input for production, then the initial mine design would typically consider a range of production decisions contingent on levels of uncertainty with respect to water supply. These states or sets of water supply, will affect production decisions and future water demand. Modelling production decisions contingent on different states of nature, termed the state-contingent approach, provides a number of insights about producer behaviour under uncertainty. This paper investigates the assumptions underlying the derivation of the market demand for water by the mineral resource industry and identifies the drivers of demand and the responsiveness of the industry to change in price and supply, particularly in relation to capital investment in water-use technologies. The usefulness of these models to the mining industry is considered by reference to the increasing demand for water by the industry in Queensland, in a climate of water shortages.

Jan 1, 2006

This paper presents a case study where Cenozoic age detritals and an ancient buried landslide have been modelled in the final walls of an iron ore pit. Orebody 25 is located in Western Australia’s Pilbara region and is an operating pit where a redesign of the final wall cutback was required. The existing slope design utilised a simplistic geological model of mixed conglomeratic, calcretic and clay –dominated detritals overlying deformed bedrock sequences. This simplicity resulted in a conservative design that was not reflective of the slopes performance to that point. Further field investigation was carried out with the results assessed in conjunction with historic drilling, mapping and laboratory testing. A detailed geotechnical model was developed including three conglomerate and two clay dominated units and laterstage calcrete overprinting. The clay units were of significant importance to the stability assessment and resulting slope design due to the thickness and poor strength conditions indicated by the laboratory and field testing. However, the drilling investigation shed only limited understanding as to the origin of the clays and thereby how they should be treated as a component in the stability analysis. A comparison of the ratio of various geochemical elements measured by assay was made between the clays and 1) shale bedrock and 2) clays of known detrital origin from other nearby sites. This indicated the clays had a strong affinity with weathered bedrock. This and other textural features resulted in the interpretation that the clay units were relic debris preserved following ancient landslides in bedrock shale. Using laboratory testing data and calibration of results for a landslide origin, stability analyses were run, resulting in a steeper slope design being implemented. Excavation of the final cutback is now in progress with the model confirmed as appropriate.

Jan 1, 2015

By T Ma, F T. Suorineni, Q He, J Oh

Directional hydraulic fracturing was invented by Polish experts and was first applied in coal mining to prevent hard roof rock bursts. Based on its successful application in coal mining, directional hydraulic fracturing has the potential to reorientate the hydraulic fracture trajectory. Recently, this technique has been introduced into cave mining for the purpose of lowering breakdown pressure. The hydraulic fracture is much easier to initiate with the assistance of an initial notch. No systematic theoretical or experimental studies have been carried out to investigate the influencing factors in the hydraulic fracture reorientation process. In this paper, the numerical modelling software Realistic Failure Process Analysis3D (RFPA3D) – Flow Version is used to study the influence of the in situ stress condition and rock mass heterogeneity on directional hydraulic fracturing in brittle heterogeneous rock masses. Firstly, the reliability of RFPA3D in simulating hydraulic fracturing is validated by comparing its results with experimental results and the simulation results from the extended finite element method. The hydraulic fracture trajectory is then observed in a three-dimensional stress environment. The results indicate that for the in situ stress condition, when the borehole lies in the ?2–?3 plane, only the difference between ?2 and ?3 has an important influence on hydraulic fracture reorientation. Furthermore, rock mass heterogeneity is a key factor that must be accounted for when fracture reorientation is required.CITATION:He, Q, Suorineni, F T, Ma, T and Oh, J, 2016. Modelling directional hydraulic fractures in heterogeneous rock masses, in Proceedings Seventh International Conference and Exhibition on Mass Mining (MassMin 2016), pp 369–384 (The Australasian Institute of Mining and Metallurgy: Melbourne).

May 9, 2016