Search Documents

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

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 Lindsay A. Robertson

The long term management of acid generating tailings at mine closure is a common issue, however, management of ARD is significantly impacted by local conditions, including climate. Acid generating tailings, over time, will produce acid rock drainage (ARD) if left exposed to water and air, which both contribute to the sulphide oxidation process. The long term closure of tailings storage areas is a concern from both environmental and engineering perspectives, with long term acid generation and dam stability the main risks. Several approaches can be taken to minimize ARD from tailings facilities using cover systems. Three tailings facility cover case histories, in areas of net positive water balances, were investigated using unsaturated flow modelling software to determine the long term potential for limiting ARD by increasing saturation of the tailings and decreasing oxygen ingress. Each case study demonstrates the importance of integrating the tailings dam engineering into tailings cover closure design. Model predictions are compared to measured saturation levels to assess the accuracy of the calculations.

May 1, 2011

By J. B. Boisvert, W. Board, C. R. Mooney

"Traditional geostatistical modelling methods are difficult to apply to heavy-tailed deposits because the extreme positive skew of the grade distribution leads to bias. The focus here is on an innovative method of modelling the high-grade population by conditionally simulating discrete volumes (nuggets) of gold with marked Poisson and Cox processes. Capping and other ad-hoc manipulation of the assay data is not required; however, an intensity map for the point processes is required. The method is demonstrated on a dataset from an epithermal deposit and demonstrates an improvement in reducing smearing and controlling the skew of the grade distribution.RÉSUMÉ Les méthodes traditionnelles de modélisation géostatistique sont difficiles à appliquer dans des gisements dont la distribution de la teneur est à queue lourde car l’extrême asymétrie positive de la distribution de la teneur génère des biais. Dans le présent article, l’accent est mis sur une méthode novatrice de modélisation de la population à haute teneur par la simulation conditionnelle des volumes discrets d’or (pépites) au moyen de procédés marqués de Poisson et de Cox. La limite d’utilisation et d’autres manipulations ad hoc des données ne sont pas requises; toutefois, il faut une carte des intensités pour le processus ponctuel. La méthode est démontrée sur un ensemble de données d’essais provenant d’un gisement épithermal et elle montre une amélioration dans le lissage du flou et dans le contrôle de l’asymétrie de la distribution de la teneur."

Jan 1, 2015

By D. Hodouin

Jan 1, 1978

By Frédéric Flament, Mary Goldman, Gille Barbery

"Pilot plant testwork on autogenous and semi-autogenous grinding was conducted at the Centre des Recherches Minérales du Québec (CRM) to investigate the effect of ball addition, solids concentration in the feed, specific energy consumption and feed size distribution.A mathematical model that incorporates breakage and surface mechanisms of loss of mass was developed to predict load size distribution of autogenous and semi-autogeneous mills. An empirical model was used to predict product size distribution of autogenous and semi-autogeneous mills.IntroductionAutogenous (AG) and semi-autogenous (SAG) grinding circuits offer potential capital and operating cost savings when compared with conventional grinding. Empirical relationships have been developed to calculate specific power consumption for SAG mills using pilot plant test results (Mular et al., 1989; Barrett and Allan, 1986). Empirical models can be obtained from plant practice data but are applicable only to the particular ore and set of operating conditions used in the model. Grinding mechanisms have to be taken into account in order to develop a model that can be applied to different ores and mills.In AG and SAG, the ore itself acts as grinding media, with the finer sizes being crushed by the coarser ones. Surface mechanisms of loss of mass, in which the inner core of the particle remains intact, also take place. In SAG milling, the effect of balls also has to be taken into account. The role of each one of these mechanisms changes with particle size and the model has to account for such changes."

Jan 1, 1991

By C. Bazin, D. Hodouin, M. Ounpuu

Sensors are now available in mineral processing circuits to monitor pulp flowrate and density, as well as ore characteristics and other process variables. The amount of information provided by these instruments is considerable but in many cases poorly exploited. This paper gives the results of the application of time series analysis techniques to signals generated by the sensors installed on an industrial Pb dezincing flotation circuit in New-Brunswick. Autocorrelation and cross correlation functions appear to be interesting tools for a preliminary circuit analysis, however the calibration of single input single output discrete transfer functions using collected sensors signals leads to conflicting results. For this reason a phenomenological approach, based on the kinetic and transport properties of the flotation process, is also considered. It has several advantages over the empirical transfer function modelling but also limitations in its applicability to industrial processes. A combination of the transfer function and kinetic modelling approches appears to be a promising avenue for better process modelling and control.

Jan 1, 1985

By S. Esmaeili

The yield strength model developed by Esmaeili, Lloyd and Poole [1] has been adapted to model isothermal aging in Cu-containing AA6xxx alloys, with (a) a moderate and (b) a very low Cu content in the solution-treated, as-quenched, condition, Isothermal calorimetry has been used to analyze the kinetics of precipitation in these alloys, While the kinetic parameters for the low Cu alloy has been derived from the present experiments, the previously determined kinetic parameters for the high Cu alloy AA6111 has been used to model the precipitation kinetics in the moderate Cu alloy. The age hardening responses of both alloys have been well predicted by the model.

Jan 1, 2005

By Ngonidzashe Chimwani, David Glasser, Diane Hildebrandt

"Milling circuits play an important role in achieving the desired product; tonnage, size distribution and properties, since they are directly linked to the comminution energy consumption for the whole process. A grinding model with a mathematical simulation was developed for an open milling circuit with a scalped feed from a hydro-cyclone. The circuit was optimized by determining the optimal conditions in terms of energy consumption which would produce a floatable product having a desired size range (i.e. -75 +9 µm) for a platinum-bearing ore. The simulation results showed an operating region comprising the best combination of the ball filling and mill speed to achieve the highest production of a floatable size class, with the minimum energy consumption.INTRODUCTIONBall milling is usually conducted in closed circuits in order to achieve a desired product size without material overgrinding. The technological importance of optimising milling circuits cannot be overemphasised (Toneva and Peukert, 2007). Previous work has brought significant insights into the benefits of modelling milling circuits (Austin et al., 1984; King, 2012; Toneva and Peukert, 2007). However, the ability to model these circuits is dependent on the performance of the classifiers and the fidelity of the grinding mill model used compared to the actual circuit (Toneva and Peukert, 2007). Austin et al. (1984) identified five types of milling circuits as shown in Figure 1 namely; the grinding circuits employing mill in an open circuit, normal closed circuit, reverse closed circuit, open circuit with scalped feed, and, the combined closed circuit. Depending on the presence of the first classifier, the second, or both as well as their corresponding loops (streams), the combined circuit (e) can reduce to; an open (d), a normal closed (b), an open circuit with scalped feed (c) or a reverse closed circuit (a). In the case of an open circuit, the mill is operated without being directly influenced by a classifier. The normal closed circuit does not have the first classifier and therefore the make-up feed is the actual feed stream. Conversely, the reverse closed circuit keeps only the first classifier. The open circuited with scalped feed has only the first classifier but does not have the second one. The overflow stream from the classifier is mixed with the mill product to form the final product for the downstream process."

Jan 1, 2016

By Ali Mortazavi

In underground mining applications, successful blasting results in less damage, which may lower support costs, lower dilution and produce a uniform fragmentation. The explosive detonation process creates a pressure front that propagates a stress wave into the surrounding rock mass at a high velocity. The stress wave can induce varying degrees of damage, depending on rock strength and physical properties, to the excavated area. This paper presents a summary of a research project that was conducted to investigate the propagation of explosive-induced stress waves in hard rock media and the mechanism of rock damage and fracturing due to shock wave propagation. The primary objective of this work was to numerically simulate the experimental study of blast-induced rock fracturing conducted within Inco Limited, Mines Research department. ITASCA?s discontinuum code UDEC (Universal Distinct Element Code), was used to simulate the results of the controlled blasting experiments. In brief, the focal point of the paper can be summarized as follows: ? Investigate the size of the stress-wave-induced fractured zone around the blasthole, ? Examine the extent and direction of the radial fractures that form outside the blasthole, and ? Explore the mechanisms that cause rock breakage.

May 1, 2002

By J. L. Feng, F. Zhang, S. K. Hu

"Cracking is one of the most serious damages in tunnel linings, which is often observed by in-situ inspection to tunnel engineering. It is difficult to better characterise for the deformations from continua to discontinua by the traditional continuum mechanics due to emerging of crack initiation, crack propagation, and coalescence, at several scales. In this paper, the process of failure that was from micro damage to macro crack in concrete tunnel linings under experimental load conditions, was quantitatively simulated by the combined finite-discrete element method (FDEM). Comparison of the simulating results with that obtained by the full-scale test of concrete lining under different loading conditions in the literature was respectively good agreed and shows that such FDEM technique is appropriate to model the continuousdiscontinuous problems. The present study indicates that the FDEM technique could be benefit to quantitatively evaluate the stability and damage level of tunnel lining system. INTRODUCTIONMost damage to tunnel linings is caused during the construction phase, which can emerge in a form of single macro cracks or even large-scale spalling. The cause of such damage is mostly the occurrence of, either impermissibly large contact stresses resulting from high thrust loads combined with production and installation tolerances in the segments resulting in a geometrical lack of fit (Maidl et al., 2012; Transportation Research Board, Basic Research and Emerging Technologies Related to Concrete Committee, 2006), or unallowable tensile/shear stresses appearing in the tunnel linings under complicated loading conditions (Ansell, 2010). Such a macro failure phenomenon may also be regarded to as a result of growth, coalescence, initiation of microcracks in tunnel linings due to the complicated loading action, and even potentially lead to an instability or collapse of the local or entire tunnel system. Much work has been devoted to investigate for the characters of micro and macro cracks, the influence factors of damage patterns, the settlement of a tunnel, the width, length and density of lining cracks, especially from classifiedly and qualitatively (Wang, 2010; Huang et al., 2013). However, the safety evaluation of tunnel is definitely required to not only inspect in situ, but further also quantify computing and analyzing for tunnel lining systems, so as to provide a global viewpoint of process of the tunnel deformation, which aims to help for subsequent decisions, e.g. controlling traffic flow or immediately repairing task. Additionally, the design of lining thickness and its structural analysis that is mostly determined by previous experience and similar construction examples is also needed to be improved."

Jan 1, 2015

By Viraj Parekh, Kinnor Chattopadhyay, Sina Mostaghel, Xiao Zhang

Electric arc furnace steelmaking is a dominant approach in North America and controlling EAF emissions is a key issue today because of stringent permit limits imposed by various regulatory bodies. In the present study, an attempt has been made to mathematically model EAF emissions, and use the models as a predictive tool for fine tuning the process. CFD modelling has been used to model various fume release scenarios, and optimize the exhaust conditions for maximum capture. Thermodynamic calculations have been used to analyze various chemical reactions which lead to the emissions of hazardous air pollutants with special emphasis on Zinc, Lead and Manganese.

Jan 1, 2015

By M. Kozlowski, I. Candy, S. Kashani-Nejad

The production of ferronickel and matte from nickel laterite ores is dominated by the Rotary Kiln, Electric Furnace (RKEF) process. Over time, major improvements have been made to kiln operation. In particular, kilns are now required to perform substantial amounts of reduction of the metallic oxides, thereby reducing the load on the electric furnace and substituting less expensive fossil fuel for electrical energy. To the basic refractory lined cylinder have been added: internal features such as dams and lifters, aimed at changing residence times in different sections of the kiln and enhancing gas/solids interaction, on-board fans to introduce tertiary air and on-board scoop mechanisms to allow reductant coal additions to be made at advantageous locations along the kiln. Many of these advances in design have been made on a trial and error basis and by 'borrowing' technology from related kiln processes, supplemented by proprietary modeling to various levels of sophistication. This paper presents a conceptual review of modeling components for the development of a computational model for granular material flow, combustion and heat and mass transfer in a nickel laterite kiln. These components include: granular material flow and segregation models for the simulation of bed motion, heat transfer models that simulate heat transfer by turbulent diffusion and radiation and the distribution of this thermal energy within the bed, and fluid flow models for the prediction of combustion aerodynamics.

Jan 1, 2005

By M. W. Bates, A. L. Mular

This paper, based on work carried out at Queen's University, describes the development and use of mathematical models of the parallel cyclone classifiers being used in the Strathcona mill of Falconbridge Nickel Min.es Ltd.

Jan 1, 1971

By G. P. Demopoulos

The solubilities of calcium sulfate in HCl-CaC1z-H20 solutions up to 100°C were successfully modelled with the aid of OLI software after bringing proper modification to it. Modelling, more specifically, was carried out with the help of the Environmental Simulation Program, V-6.6.0.4 (www.olisystems.com). New model parameters were obtained by regressing the solubility data. A private databank was established by adding the newly defined model parameters and the calcium sulfate hemihydrate phase that is not included in the original OLI database. The McGill-modified OLI software was successfully applied not only to modelling solubilities but also the transformation of the various CaS04 phases, namely dihydrate, hemihydrate and anhydrite. In this paper, the modelling methodology will be discussed along with a comparison of the capabilities of the newly developed model vs the original OLI model.

Jan 1, 2005

By X. Z. Yuan

Contamination introduced by the impurities from the feed stream can impact the PEM fuel cell performance dramatically. The presence of these unwanted trace species such as CO, H2S, NH3, can adversely affect the fuel cell functioning. It has been reported that the major impact of CO and H2S contamination on the fuel cell performance is kinetic while the effect of NH3 contamination is speculated as membrane conductivity reduction only. In this paper, a model regarding the NH3 contamination is described and the mechanism of the NH3 contamination is illuminated by the model calculations.

Jan 1, 2005

By Carl Bertrand, Patrice Nadeau, Claude Bazin

"Rio Tinto Fer et Titane processes the ore extracted from the Lac Tio orebody to concentrate the hemo-ilmenite (FeO.TiO2-Fe2O3) mineral that is subsequently smelted to produce pig iron and titanium oxide slag. The ore is crushed, heat treated (roasted) and concentrated using dry magnetic separators. This paper describes a model to predict the magnetic susceptibility of the roasted ore and its behaviour in the magnetic separation plant. The model uses a classification of the ore as a function of the size and magnetic susceptibility of the particles. A subsequent division into mineral species allows the simulation of the magnetic separation process. The model is used to assess various operating options when process disturbances occur.INTRODUCTIONOre processing operations consist of ore size reduction followed by the recovery of valuable species in a concentrate that is subsequently sold for its minerals or metal content. The concentration processes are based on the physical and/or chemical properties of the minerals contained in the ore. Flotation is typically used to recover minerals of base metals such as copper, zinc, lead or nickel. Iron ore and other dense minerals are typically concentrated by gravity separation processes. Magnetic separation uses the magnetic properties of the minerals to concentrate magnetite and other minerals, to clean coal concentrate and to recycle ferrosilicon for dense separation processes (Brophy, 1985).Magnetic separators can be divided into low and high intensity separators and high gradient magnetic separators (Gillet, 2008). Most of the low gradient magnetic separators are built as rotary drums that are operated in wet or dry conditions. Dry magnetic drum separators are used by Rio Tinto Fer et Titane in Sorel-Tracy, Canada to concentrate hemo-ilmenite from the ore extracted at the Lac Tio Lake mine near Havre St-Pierre, Canada. The analysis of production data suggests that process improvement could be achieved by identifying key variables for the control of the magnetic separation circuit and/or by modification to the circuit flowsheet. The use of a process model was considered to be an effective approach to conduct a preliminary analysis of the circuit with a view toward process optimization. However, the literature contains few detailed studies dealing with the modelling of magnetic separation processes."

Jan 1, 2011

By M. R. Khalesi, S. Bellec, C. Bazin, D. Hodouin

Batch grinding tests conducted on a gold ore indicated that gold is not randomly distributed within the size fractions of the ground ore. Since gold leaching performances were found to be dependent of the size of the gold carrying particles an integrated gold ore grinding-leaching simulator should incorporate a method to predict the gold size distribution of a ground ore. The paper describes the results obtained from batch grinding tests showing the non random distribution of gold within the size intervals of a ground ore and the initial developments of a grinding model to predict the gold size distribution in a batch grinding product.

Jan 1, 2011

CFD (Computational Fluid Dynamic) models have been developed and validated for a pulsed column with disk-doughnut internals, as used for solvent extraction. First a single phase model was developed to simulate the turbulent oscillating flow around the internals of the column. This model was successfully validated using detailed PIV measurements of velocity taken at eight different phases through the pulse cycle. Second, the CFD model was extended to allow for a continuous aqueous phase and a dispersed organic phase (kerosene) flowing counter-currently. The model was a two-fluid model that accounts for droplet breakup and coalescence, and solves for droplet size distribution using a population balance model incorporated into the CFX4 package. PDIA (Particle Droplet Image Analysis) measurements were taken on a pilot scale solvent extraction pulsed column to determine droplet size distribution for calibration and validation of the CFD model. The measurements were made under pulsing intensities ranging from 10 mm/s to 32.5 mm/s.

Jan 1, 2007

By Nigel Saunders, Zhanli Guo

"Quench sensitivity of heat-treatable aluminium alloys is closely related to the precipitation process taking place during quenching. Faster cooling results in less precipitation of coarse phases during cooling, leaving more solutes in the solution before ageing. An alloy in such state would be of greater hardening potential as larger amounts of hardening phases may precipitate out during ageing. Hence, to understand the quench sensitivity of an alloy, it is essential to understand its precipitation process. Precipitation is a diffusion-controlled process and the diffusion is made complicated by the so-called quenched-in vacancies. These vacancies form during solution treatment, become “excess” when temperature goes down, and annihilate during the following cooling and ageing treatments, making diffusion now a function of both temperature and time. This paper first investigates the formation and annihilation of quenched-in vacancies and their effect on diffusion. The diffusion affected by quenched-in vacancies is considered in the kinetic models to realise the calculation of TTT/CCT diagrams for aluminium alloys. The calculated CCT diagrams have been used to explain the observed quench sensitivity and age hardening behaviour of various commercial alloys. The transformation from GP zones to other hardening phases during ageing is also discussed.INTRODUCTION The heat treatment of aluminium alloys mainly consists of three stages: a solution treatment at an elevated temperature, quenching from this solution treatment temperature to room temperature, followed by an ageing treatment to allow the precipitation of strengthening phases(Davis, 1993).The following statements generally describes how quench sensitivity is observed in heat-treatable Al-alloys: faster cooling results in,(i) less precipitation of coarse phases formed at high temperature during cooling, (ii) more solutes kept in solution before ageing, and (iii) larger amounts of hardening phases being formed, i.e. a greater hardening potential. Thus, the quench sensitivity of an alloy is closely related to its precipitation kinetics."

Jan 1, 2018