Search Documents

By P B. Hills

This paper discusses the computerisation of the drill hole database, some historical records and the creation of a three-dimensional model of the gold bearing structures for resource estimation and mine design at the Beaconsfield Gold Mine in northern Tasmania. The phase of the work reported here commenced in November 1997 and was completed in March 1998. The historical mining on the Tasmania Reef had reached a depth of 450 m below surface when work ceased in 1914. Production up to this time was 854 000 ounces gold. Recent activity has concentrated on testing depth extensions of the previously mined ore structure, to a depth of 1000 m below surface, and the establishment of an underground mining operation to exploit this resource. The lode structure at Beaconsfield is a typical mesothermal vein style deposit. Gold is contained within narrow quartz, carbonate and minor sulphide reefs or veins. Reef width is variable from a few centimetres to several metres with æsplitsÆ or bifurcations being common. The reef averages 2.7 m true width in the economic zones. The key feature of this style of mineralisation is that reef continuity can be confirmed by drilling to a high level of confidence, but the gold grades are somewhat more erratic in nature. The resource estimates were based entirely on drill data. The drill intercepts on the reef were coded for each split directly into a database and these allocations validated graphically. The distribution of grades and widths show an increasing gold grade with reef width for both the drill hole and historical data. The same histogram patterns, mixed populations and overall arithmetic average grade for the historical mine (at 21.52 g/t > 3 g/t Au) and new drill data (at 21.55 g/t > 3 g/t Au) would indicate that there has been no major change in the overall style or grade of the ore from the historically mined areas to the new resource. The reef, host geology, area of historical workings and resource zones at the Beaconsfield Mine were defined as a series of closed wireframes. The Tasmania Reef was modelled as five mutually exclusive wireframes, one each for the Main Reef, Footwall split, Second footwall split, Third footwall split and Hangingwall split. The data for this exercise was extracted directly from the drill hole database, stored in Microsoft Access format. The formation of the solids was largely automated with manual validation and refinements allowing for multiple runs using different interpolation methods and combinations of drill holes. A rotated 3D-block model in the plane of the Tasmania Reef was created to contain the grade model and allow for tonnage, grade and reef width estimates. An inverse distance squared (ID2) interpolation method was used with ranges of 25 m along strike and 125 m down plunge with blank areas infilled with 50 m by 250 m. Each reef split was filled separately using only drill intercepts from that split. The method was independently validated used MIK methods. Longitudinal sections have been prepared showing gold grade, reef width and width +grade. Estimates were also made for arsenic, sulphur and other minor elements. The long sections show three distinct mining zones. A method is presented for the conversion of the resources to reserves by the addition of mining parameters for minimum mining width, minimum pillar, cut-off grade and internal and external dilution. The selection is done on columns within the block model, and allows for all splits to be considered at once. Multiple runs can be made for varying combinations of parameters to test different mining methods and their impact on reserves, mine life and NPV.

Jan 1, 1999

By S Bhandari

This paper is a step in the direction of developing a computer model of the near source dispersal of suspended particulate matter (SPM) in planetary boundary layer (PBL). Since Surface Mine blasting is done during daytime, study of daytime dispersion in convective boundary layer (CBL) is relevant to the mining engineers. The scope of this paper covers dispersal after the blast in any type of PBL ie stable, convective and neutral but only under weak wind condition. Mathematical analysis of the dispersion of vertical dust column has been carried out during weak wind conditions in a tropical meteorological environment. It is noted that dust concentration reaches its maximum value some distance away from the point of blast, immediately after the blast (t*=1) and thereafter rapidly decreases. After a short interval from the time of blast at t*=11 the dust concentration shows linear decrease in concentration with distance. Dust concentration also rapidly drops during the interval t*=1 to t*=11, thereafter, rate of drop is nearly the same at any distance from the point of blast.

Jan 1, 2001

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 Douglas Morrison

By applying modern 3D numerical modelling techniques to the analysis of a pillar cascade failure in a room and pillar mine, it has been possible to simulate the collapse and develop an understanding of the role of low rock strength areas within the mine. The cascade collapse occurred in a normal strength area but would not have occurred without the influence of a low strength rockmass in one section of the mine. The numerical analysis made it possible to develop a realistic assessment of the residual pillar strength. With the introduction of regular empirical pillar performance assessments, critical pillars or groups of pillars can now be identified proactively. This has led to an aggressive extraction ratio of 90% that manages pillar failure, avoids pillar and roof collapse and prevents surface subsidence. Since the mine also has significant waste products, the use of backfill for ground support has made it possible to significantly improve the performance of strategically located pillars, and so achieve even higher extraction ratios.

Jan 1, 2003

By J. Yianatos

The regrinding stage is necessary to achieve the particle size (liberation) for the final upgrade in the cleaning flotation process. At present, vertical mills (Vertimill) are employed at industrial scale, using small steel balls as grinding media. The regrinding process characterization requires mixing regime knowledge along with the mean residence time of the liquid and solid phases. The residence time distribution (RTD) and the actual mean residence time estimations from theoretical relationships are complex due to the presence of three components (solid mineral, balls and water), mass transport characteristics and segregation affecting solid holdup. Therefore, the experimental RTD determination is suitable for industrial mills. In this work, RTD measurements were performed in two different industrial flotation plants. Both plants operate with Vertimills of 1000-1500 HP in the regrinding process with different circuit arrangements according to the plant treatment capacity. The RTD was measured using the radioactive tracer technique. This technique allows for non-invasive tracer detection where the procedure consists of selecting a liquid (Br82 in aqueous solution) or solid (irradiated solid) tracer for RTD data acquisition in real time. The RTD was modelled by the Large and Small Tanks in Series (LSTS) model, which allowed the characterization of the different mixing regimes for both liquid and solids in the regrinding mills. Effective mean residence times in the range 2-11 min were obtained with RTD shapes similar to those obtained in tumbling ball mills. It was found that residence times of liquid and fine particles were similar, which is consistent with results observed in other mixing systems, e.g. mechanical flotation cells. Also, the mean residence times estimated from the RTD allowed the mass balance adjustment to be validated from metallurgical sampling data. RTD estimations are useful for better understanding the Vertimill behaviour and to identify the actual breakage rates. Keywords: vertical mill, residence time distribution, particle size, modeling

Sep 1, 2012

It is desirable to have a model, or method, for major project development, which will enable flexible modelling of separation circuits. This ability improves the understanding of circuits and how they should be operated, especially in a design phase, where various options are evaluated and decisions are made on the required number of stages and configuration. This paper describes existing gravity separation circuit modelling theory and concepts, as defined by Holland-Batt1 (separation window, separation efficiency and proportionality of efficiency). The crux of the paperis the application of numerical analysis to fit raw separation data by means of a modified Bézier-Spline. The paper proposes an alternative to the theory developed by Holland-Batt that provides more flexibility in the modelling of separation stages and curves, which include S-shaped curves. The paper proposes a solution to the challenge of continuously extrapolating data, and the first difficulty normally encountered, which is the fitting of an appropriate empirical formula to raw data.

Jan 1, 2007

By Mandzic E. H

The Dome Mine of Timmins Mining company is located in South Porcupine, N.W. Ontario. The orebody being mined partially lies in the No. 3 shaft pillar. Stability analysis has been carried out to establish the influence of nearby mining on ground movement in the shaft vicinity. Mathematical modelling of shaft displacement has been done for stope layout at depths from 180-210 m and distances of 10-50 m from the shaft. Analysis indicated that ground displacement around the shaft is below 1 mm/m and one third less for cut and fill than in the case of open stope mining. The result of the investigation was later confirmed by cut and fill mining of this orebody.

Jan 1, 1986

By G. Brooks, J. Naser, M. A. Sattar

CFD models have been developed and numerical simulations have been carried out to predict the formation of foam in oxygen steelmaking. Foam was considered as a separate phase comprising a mixture of gas and liquid. Bubble break up and coalescence models have also been incorporated in a CFD model to predict the number density of individual bubble classes. A population balance equation was used to track the number density of each bubble class. Decarburization with heat generation from chemical reactions was integrated in the process. User subroutines were written in FORTRAN to incorporate the foam formation, the bubble break up and coalescence rate and decarburisation in the main program. The model predicted the foam height, bubble number density, velocity of phases, decarburization, and turbulence. The result from the model has been compared with available data from literature and found to be in reasonable agreement with the experimental and plant data.

Jan 1, 2014

The entrainment of copper matte in the slag during the slag skimming stage of a Peirce-Smith converter was studied with a dam break model. Experiments were carried out with a range of fluid pairs where the lighter fluid was allowed to overflow a weir and the amount of entrained heavier fluid measured. The resulting entrainment is measured as the ratio of the heavy to light fluid withdrawn. The fluids include oil-water, water- glycerol and salt solutions-water pairs to cover a wide range of densities and viscosities. It was found that the Richardson number, Ri, correlates the entrainment well. For the case of the liquid-liquid interface initially at the same height as the weir, the entrainment was found to vary with Ri-1/3. The transient flow was simulated with a volume of fluid (VOF) code written in-house and used a stairstep approach to represent flow obstacles. The code was validated against the experimental results with good agreement for the variation of the heights of the two fluid with time. The simulation showed that the prevailing mechanism for entrainment is due to a pressure differential created across the liquid in the converter during pouring resulting in the heavier liquid being rotated upwards and over the weir.

Jan 1, 1999

By J. L. Liow

Studies have been carried out on the causes and effects of slopping in an annular cylindrical geometry such as the top submerged combustion reactor. Properties of standing waves formed such as the wave profile, period, pressure distribution and velocity have been derived analytically. The wear profile due to wave erosion was calculated based on the assumption of linear motion and that the wear is a dissolution process. The use of the j, factor showed that the wear was proportional to the square root of the wave velocity at any point. Using the dimensions of the lsasmelt reactor given by Errington et al. (1987)[11, the wear profile due to wave erosion was calculated using the analytical solution. The wear profile calculated was shown to have a maximum just below the initial bath surface. The wear was found to be concentrated in the region where the first asymmetric standing wave was most active. The results showed that the refractory near the bath surface will be affected by wave erosion, more so than any other parts of the reactor. The contribution of wave erosion on wear must be accounted for in assessing the total refractory wear of a vessel in which slopping exists.

Jan 1, 1991

By A. Chagnes

This paper presents a modelling tool that can be used to predict the behaviour of solvent extraction processes. In particular, this model allows calculating metal extraction efficiency as a function of pH, solvent flowrate, stripping flowrate and aqueous feed flowrate in the mixers-settlers. This model is applied to the recovery of cobalt from acidic chloride media by using CYANEX® 272 as extractant diluted in kerosene. Keywords: Cobalt, mixer-settler, modelling, solvent extraction

Jan 1, 2020

By J. N. Carras, N. W. Bainbridge

In order to further the understanding of spontaneous heating of longwall goafs a numerical model was developed, based on conservation principles applied to the transport of heat, water vapour and oxygen across the god Advective air flow into the goaf is caused by the differential air pressures across the god The air velocity field within the goaf area was computed using Darcy's law. The boundary and initial conditions concern air pressure, temperature and relative humidity. The oxidation rates and heat generation rates of the goaf materials required for the simulation were obtained from laboratory measurements. The main outputs of the model are temperature, oxygen and water contour maps, the rate of air leakage into the go4 velocity field and pressure distribution across the god. The model while still in its development stage, is able to simulate a heterogenous goafs constructed from different materials with different thermal and mass flow properties. Validation of the model by comparing its output with measurements from within goafs is required before the model can be used with confidence in a predictive sense.

Jan 1, 1995

By P. Myllykoski, J. Nylander

"A commercial program package called Crownoff for the simulation of temper rolling in four high mills has been verified against measured data. The agreement between measured and calculated transversal strain distributions of the strip was found to be good. Effects of the strip crown and work roll bending force on the strain distribution were simulated. The results showed that when a single ground camber of work rolls is used for all strip widths; the strip crown should be parabolic, excluding the edge drop. An increase or decrease of the strip crown by 30 % during a coil can be compensated by a 2 to 5 t change in the work roll bending force. The ability to control the strip crown by work roll bending was insensitive to strip width. IntroductionCold rolled steel strips are usually temper rolled after recrystallization annealing, which follows cold rolling. The purpose of temper rolling is to suppress the yield point, finish the strip surface and correct possible flatness errors. The presence of a yield point might otherwise cause problems during subsequent forming processes. The in homogenous yielding in non-temper rolled steels leads to the formation of Lüders bands in the strip surface.The rolling force is influenced by the rolling speed and changes in frictional conditions in the roll bite. Therefore it is very important to plan carefully the mill settings and control algorithms during rolling. The purpose of mill control algorithms is to ensure constant reduction and the best possible strip flatness during the entire coil.Correctly predicting the rolling force is important for determining the mill presettings. Torque, forward slip and power requirements are not usually critical, if the reductions are small. However, the large elastic deformations occurring in the roll bite make the physical description of the process difficult. Physical rolling models developed for cold rolling contain assumptions that lose their validity as the reductions decrease to the level of one per cent."

Jan 1, 1994

By R. Nilsson

In the present work, an assessment of the sulphide capacities of a number of binary and ternary slag systems have been carried out in the light of a sulphide capacity model. The model is based on the expression: {CS- exp (-A G °/R T) ? ao2 / fs2-} where AGO is the standard Gibbs energy change for the above reaction, ao2 is the activity of oxygen ions and f, is the activity coefficient of the sulphide ion in the slag. In the present model, the ratio (ao2/fs2) is optimized for the experimental data available in literature. The systems assessed are the binaries MnO-SiO2, CaO-SiO2, Al2O3-CaO, as well as the ternaries CaO-Al2O3-SiO2 and CaO-MnO-SiO2. In the case of all these systems, the experimental data available in literature were complemented by measurements carried out at the Division of Theoretical Metallurgy. The sulphide capacities of some FeO- containing binary, ternary and quaternary slags are also examined. The potentiality of the model to estimate the sulphide capacities of complex slag systems from binary and ternary data is demonstrated.

Jan 1, 1996

Modelling of Synergistic Effect of Cyanex 302 and D2EHPA in Zinc Solvent Extraction A comprehensive model is presented to predict the synergistic effect of D2EHPA and Cyanex 302 on co-extraction of zinc and manganese. The experimental results indicated that pH, temperature and the extractants ratio have strong effects on the extraction process. Correlations were proposed to predict the distribution coefficient of zinc and manganese between two phases at different conditions by multiple linear regression. The correlations were simplified by eliminating less important terms by analysis of variance. Optimum conditions (ie pH, temperature, extractants ratio and aqueous to oil phase volumetric ratio) for obtaining maximum separation of zinc from manganese at different zinc concentrations and number of extraction stages were determined. For finding the optimum condition genetic algorithm was used. Comparison the model estimations and the experimental data show that the presented model has good ability to predict the behaviours of zinc extraction system over Cyanex 302 and D2EHPA.

Sep 13, 2010

By Tran T, Young D

During the processing of gold, effluents that are discharged into the tailing dams contain free cyanide, metal cyanide complexes of zinc, copper, iron, etc and their oxidised products such as ammonia, cyanate and thiocyanate. A knowledge of how different cyanides behave in the solution is vital to metallurgists working in the area of cyanide waste treatment. Although the free cyanide ion is not stable in aqueous solutions and is expected to degrade in a natural environment by atmospheric oxidation, the metal cyanide complexes, especially of iron, are very stable. The CHEMIX program from the THERMODATA package, developed by the CSIRO Division of Mineral Products, was used to evaluate various copper-cyanide and iron-cyanide systems. The simulation of these systems show how the free cyanide and cyanide complexes behave at different conditions of solution composition and acidity (pH).

Jan 1, 1991

By T. C. Rao

Oil agglomeration experiments were carried out on coal fines to study the effect of process variables such as pulp density, bridging oil dosage and intensity of agitation on (i) size distribution and (ii) yield of agglomerated product. It has been observed that the growth of agglomerates follows a self-preserving behaviour and using this, a characteristic curve has been derived. This curve has been found to be independent of the process variables. It has also been noted that the growth of agglomerates follows second order kinetics. In an attempt to model the process for predicting the size distribution and yield of agglomerates at different levels of the variables, the cumulative 50 percent of the agglomerates passing size (d50) has been correlated with the variables (i) through the kinetic equation and (ii) directly. These equations, coupled with the characteristic curve can predict the size distribution of the agglomerates at any levels of the variables tested. The yield of agglomerates has been found to follow a definite relationship with d50.

Jan 1, 1984

By Jean-Christophe Gebelin, Mark Jolly, Samantha Jones

"The de-waxing process in investment casting is truly a black art which is carried out in a black hole - the autoclave. The process involves exposing wax, coated by a porous ceramic shell, to steam at a pressure of 10 bar and a temperature of about 180°C. The processes which then occur are truly multi-media and multi-physics. Steam condenses in the porous ceramic, wax melts and flows, the shells expand and sometimes crack. Why the shells crack is one of the most difficult questions asked in the industry.Modelling and validating the de-waxing process is not simple. The IRC at the University of Birmingham is attempting to model this process using Physica+ a code developed in the UK. This paper will present the latest results in this attempt to model one of the casting industries most intractable problems.IntroductionThe de-waxing process using an autoclave involves the following physical phenomena:· Flow of the steam in the autoclave;· Heat transfer from the steam to the wax, possibly through the shell;· Melting and flow of the wax.At the same time, other phenomena may occur, such as the flow of the steam through the shell, the condensation of the steam on the shell surface and on the wax pattern surface, the flow of wax the through the shell, etc ...This shows how complex the de-waxing process is from a physical point of view. The main defect that appears during this process is the cracking of the shell. Previous modelling work carried out by Snow [1] has shown that it is possible to predict the stress in the shell, making numerous assumptions. In this paper, the modelling to predict the state of stress in the shell that is currently in development at the University of Birmingham is presented together with the presentation ofthe physical phenomena to model and the assumptions used. The parameters that are needed in such models will also be discussed, and some experimental work undertaken at the moment presented."

Jan 1, 2001