Search Documents

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

Sep 13, 2010

By K Williams, T Donohue, W Chen, T Roessler, A Katterfeld

Wet sticky iron ore poses particular difficulties in every aspect of the material handling process. Material hang-ups in bins, pipes and chutes; carry-back and spillage due to insufficient belt cleaning; clogging of screens; malfunctioning of crushers as well as rat holing of gravity reclaimed stockpiles are all predominantly due to the inherent cohesion and adhesion characteristics of the wet, sticky ore (Roberts, 1993). Limited understanding and quantitative measures on the stickiness of such a material significantly lower the mining operation efficiency.This research aims to investigate the cohesion and adhesion characteristics of the wet, sticky ore on physical modelling, experimental techniques and numerical implementation. Initially, the suitable bridging mechanism that contributes to the cohesive and adhesive forces within the iron ore material was identified (Hering, Martin and Stohrer, 1989; Rumpf, 1958). A cohesion/adhesion energy density factor based contact model was then introduced to describe the stickiness of the material. Secondly, an experimental system was proposed to directly test the stickiness of the material. Such a system is able to obtain both the cohesive/adhesive shear force and tensile force of the material at a specific moisture content. Lastly, the implementation of the proposed cohesion/adhesion model into the discrete element modelling (DEM) software was also discussed. A cohesion/adhesion calibration for modelling wet sticky iron ore in DEM was then conducted to ensure the material modelled in the simulation is reflective of the reality.From the study, understanding of the cohesion/adhesion characteristics of the wet iron ore material is significantly enhanced. Experimental and numerical outcomes can be directly applied for better designing or optimising the material handling processes that involve wet sticky iron ore. Consequently, the mining operating costs can be greatly reduced.CITATION:Chen, W, Donohue, T, Williams, K, Katterfeld, A and Roessler, T, 2015. Modelling cohesion and adhesion of wet, sticky iron ores in discrete element modelling for material handling processes, in Proceedings Iron Ore 2015, pp 387–392 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Jul 13, 2015

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

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

By R. P. Jensen

This is a sampling of the papers presented at the conference, the proceedings can be bought from the SAIMM Arbitrary Lagrangian Eulerian (ALE) computational techniques allow treatment of gases, liquids and solids in the same simulation. ALE methods include the ability to treat shock waves in gases, liquids and solids and the interaction of shock waves with each other and with media from one of the other categories. ALE codes can also treat explosive detonation and the expansion of the explosive gases and their interaction with air and solids. ALEGRA is a 3-D ALE code that has been developed at Sandia National Laboratories over the past few years. ALEGRA has been applied to a 2-D simulation of presplitting using decoupled explosives in rock blasting with very interesting results. The detonation of the explosive at the bottom of the hole sends a shock wave up the borehole driven by the explosive gas expanding into air. The explosive gas compresses the air against the stemming column where it rebounds and recompresses at the bottom of the borehole. This type of ringing takes several cycles to damp out. The explosively induced expansion of the borehole is also treated by ALEGRA as well as the shock wave imparted to the rock. The presentation of this paper will include several computer animations to aid in understanding this complex phenomenon.

Jan 1, 2000

By W Valery, S Morrell

In the mining industry blasts are usually designed to fracture the in situ rock mass and prepare it for excavation and subsequent transport. The run of mine (ROM) fragmentation is considered good when it is fine enough and loose enough to ensure efficient digging and loading operations. Mining optimisation strategy is hence usually focussed on minimising total mining costs and maintaining these ROM fragmentation characteristics. Although this approach ensures an efficient mining operation it ignores the potential impact on crushing and grinding. Investigations by several researchers have shown that designing blasts to produce ROM fragmentation to optimise crushing and grinding performance enhances the overall efficiency and profitability (Eloranta, 1995; Kojovic et al, 1995; Bulow et al, 1998; Kanchibotla et al, 1998; Scott et al, 1998; Simkus and Dance, 1998). Field experimentation in this area is often found to be difficult because of high implementation costs and insufficient understanding about the effect of fragmentation at different stages in a comminution circuit. Modelling and simulation of blasting and comminution processes provide a more economic alternative to explore the impact of blast design changes on the down stream operations (Kanchibotla et al, 1998; Kojovic et al, 1998; Adam and Siddall, 1998; Bulow et al, 1998). However, the reliability of such simulations depends on:the ability of the models to reproduce the physical processes that take place; characterisation of rock mass in terms of its blastability and comminution characteristics; and in the case of blast modelling the ability to predict fines (-10 mm) has been found to be particularly important as this fraction has a significant effect on mill throughput.

Jan 1, 1999

By Emöke Imre

A point-symmetrical linear coupled consolidation model with displacement boundary condition is suggested for the description of the consolidation around the pile tip after pile penetration. Analytical solution is given. Main features of the solutions are predicted by the analysis of the system of differential equations and, the boundary conditions.

Jan 1, 2002

Ore deposit modelling is becoming an increasingly accepted tool in mineral exploration. Successful models have already been proposed for several ore deposits of base metals and gold, but tin has so far received comparatively little attention, perhaps because of an excessive preoccupatioi with alluvial operations. This situation is rapidly changing and some of the world's most profitable and newly commissioned tin mines are large, primary tin deposits. Some basic principles of ore deposit modelling are briefly reviewed and the fallacy of an excessive reliance on grade tonnage diagrams based on the tin industry worldwide is discussed. The four proposed models and submodels of large, hard-rock, tin deposits focus on basic concepts such as the spatial relationship with "tin granites", depth of emplacement and the level of erosion. General concepts are drawn. The modelling of large primary tin deposits is reviewed not only in its geological and exploration aspects but also in its mining and economic implications.

Jan 1, 1986

By Harald Ulrik Sverdrup, Anna Hulda Olafsdottir

The long-term supply of nickel to society was assessed with theWORLD7 model for the global nickel cycle, using new estimates of nickel reserves and resources, indicating that the best estimate of the ultimately recoverable resources for nickel is in the range of 650–720 million ton. This is significantly larger than earlier estimates. The extractable amounts were stratified by extraction cost and ore grade in the model, making them extractable only after price increases and cost reductions. The model simulated extraction, supply, ore grades, and market prices. The assessment predicts future scarcity and supply problems after 2100 for nickel. The model reconstructs observed extraction, supply and market prices for the period 1850–2020, and is used to simulate development for the period 2020–2200. The quality of nickel ore has decreased significantly from 1850 to 2020 and will continue to do so in the future according to the simulated predictions from the WORLD7 model. For nickel, extraction rates are suggested to reach their maximum value in 2050, and that most primary nickel resources will have been exhausted by 2130. After 2100, the supply per capita for nickel will decline towards exhaustion if business-as-usual is continuing. This will be manifested as reduced supply and increased prices. The peak year can be delayed by a maximum of 100 years if recycling rates are improved significantly and long before scarcity is visible.

Jan 26, 2021

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

One of the challenges in the development of mineral processing simulators is the modelling of liberation characteristics of particles produced by comminution circuits. Understanding the mineral characteristics of these particles is one of the keys to predicting their behaviour in separation circuits. While more fundamental research seeks to model the non-random breakage of ores this work examines whether we can move forward now by using observed characteristics of ores to model the composition of particles generated in comminution processes. The current work examines a range of ores and identifies which of these exhibit the pattern of similarity of composition for particles in a given size class that has been produced as the result of a variety of breakage conditions. The impact of mineralogy and texture on the evolution of liberation characteristics for both valuable and gangue minerals is discussed.

Jan 1, 2009

By A Lyon

The minerals industry, as many other industries, has a wealth of indicators and a thriving æindustryÆ in their development, measurement and reporting. It is natural for this to occur for a new endeavour. However, perhaps it is time for the community concerned with minerals industry indicators to ask whether we have achieved a certain maturation whereby proliferation of indicators may be replaced by a more disciplined and constrained approach. In many natural and engineered systems, thermodynamics provides well-understood and tested constraints. Indicators that are uninterpretable in this framework are not considered acceptable. This paper provides initial evidence that such a formalism may be possible for mine sustainability. The formalism developed does not claim to have the rigour of thermodynamics. However, it does provide a first step towards formalising selection of indicators on the basis of their true comprehensibility. A set of simple equations is developed that represent the conventional five capitals (actually we use six because we separate renewable and non-renewable natural capital) and how they are linked. This is coded into a simulation model. The base equations are then used to develop a system level indicator of sustainability performance, which is asserted to be a reflection of the mineÆs reputation. We define a least sustainable state, termed æmaximum revenue mining (MRM)Æ in the system and define the degree of system sustainability as how far the system is from the MRM. We then introduce a simple mechanism whereby different actors can express differing views on the emphasis for the distance from MRM. Sustainability thereby becomes a non-unique thing that is defined by a distance from MRM and a set of preferences across the capitals.

Jan 1, 2009

By A Lines

This paper addresses a gap in longwall multi-seam interactions and their predictions, providing a case study for future multi-seam operations for both anticipated conditions and a method for understanding the changes in stress regime. Multiple panels were modelled with the results analysed and compared to panel experiences to validate the parameters and properties. The software chosen (LaModel) is a numerical modelling package designed to model stresses and displacements on thin deposits. It utilises a displacement-discontinuity variation of the boundary-element method and has been used extensively throughout both the United States and Australia in multi-seam operations. Analysis of the results demonstrates a clear correlation with the above workings and regions of significant strata instability. This provides the opportunity going forward to analyse new models and incorporate into mine planning and geotechnical designs. Continuous site-based assessment should be undertaken at multi-seam operations to better understand and quantify likely vertical stress levels within the workings. With a detailed analysis of models prior to extraction and a review of models post extraction. Further research is required as LaModel is limited to vertical stress and it is likely that interactions impact both the horizontal and shear stress regime. 3-dimensional modelling could further assist in this understanding, combined with a thorough stress monitoring regime. CITATION:Lines, A, 2018. Modelling multi-seam interactions for longwall mining, in Proceedings The Fourth Australasian Ground Control in Mining Conference (AusRock), pp 185–195 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Nov 30, 2018

By Lars Arne Granlund

Due to the coarse structure of mining explosives, one has to model a finite reaction zone,contrary to C-J calculations, describing their detonations, i.e they become dependent of diameter and confinement, An effort within the DYNO group has been to investigate if the DSD theory ( originally developed by Bdzil, Stewart / Los Alamos ) can be used to describe major groups of civil explosives.The groups of explosives tested are dynamite, anfo, aluminised emulsion, heavy anfo and straight emulsions.

Jan 1, 1995

By E. Lefeuvre

Precambien in age, the phospate deposit under study which outcrops to the south forms a synclinal consisting of only one mineralized horizon. The synclinal strikes SSW/NNE, dips from 30 to 60 NNE and is about 5 km wide. The deposit has been explored by : - trenches and pits about 50 m apart in the outcropping area ; - core drill holes forming lines perpendicular to the outcrop. These lines are about 100 m apart. On each line, the drill holes are about 50 m apart The phosphate layer is clearly visible and has then been defined from the geologicia1 description of drill hole logs. The layer has been divided into two parts : a very rich, altered part at the hanging wall and a poor part displaying opposite physical characteristics, i.e. a fresh and hard rock. With open pit mining in view, the deposit has been divided into a 3-dimension network of blocks, each block measuring 50 X 50 X 100 m. The tonnages of rich and poor ore, PO25, MgO, Si02, CaO and R2O3 grades and the tonnages of barren waste at the hanging wall and the footwall have been calculated for each block. Given the dip of the phosphate layer ranging from 30 to 60° and its relative thickness (between 10 and 15 m), it was impossible or inappropriate to carry out direct calculation by 3-dimension kriging : too many neighboring blocks to obtain sufficient date and impossibility to draw up the morphology of the layer.

Jan 1, 1986

By S. Witzke, J. Lehmann, H. Gaye

The modelling work presented here concerns Ni-bearing steels produced by Imphy Ugine Précision. In such grades, the contents of elements such as silicon or titanium must be precisely controlled by slag/metal treatments. The slag phase is described by IRSID slag model that has been slightly revised in order to better fit recent SiO2activity measurements in the low SiO2content domain of the CaO-Al2O3-CaO system. The metallic alloys are described by the unified interaction parameter model (UIPM). In this formalism, the interaction coefficients, giving the effect of nickel on Al, Si, Mn and Ti activities, have been determined by electrochemical measurements and/or extracted from industrial data on slag/metal treatments concerning a wide range of alloy compositions: [ ] Keywords: slag, Ni alloyed steels, activities, interaction coefficient

Jan 1, 2004

By J. M. Font, R. G. Reddy

Jan 1, 2004