Search Documents

By W. Robert, J. Wirth, G. Grasselli

ABSTRACT: Strength, deformability and fluid flow properties of rock joints are to a great extent controlled by the surface roughness. The hydraulic conductivity of a rock fracture depends on the aperture distribution, the surface roughness and the contact area, each of which depend on the stresses that act on the fracture plane. However, despite much recent work, a complete understanding of the relationship between void space geometry, contact areas, applied stress, and hydraulic conductivity has yet to be achieved. Moreover, although in situ rock joints are subjected to both normal and shear loading, few studies have been conducted to investigate the effect of shear displacements on fluid flow through a single rock fracture, and which have addressed the issue of flow anisotropy induced by shearing. Experimental flow tests show that with increasing shear displacement, the fracture becomes heterogeneous and anisotropic, and becomes more permeable in the direction parallel to the displacement (Yeo et al., 1998). On the other hand, dry shear tests show that the surface fails in areas facing the shear direction (Grasselli and Egger, 2003). Moreover, the experimental results confirm that the shear strength of rock joints is direction-dependent. In this paper, a three-dimensional morphological description of several a natural fractures that have undergone shear displacement were generated using the Automated Topographic System at the University of Basel. The data analysis method of Grasselli (Grasselli et al., 2002a) is used to determine the surface parameters that are known to control the shear strength. Simulations of fluid flow through this fracture, in various directions relative to the direction of shear, are conducted using an FEM flow simulator. The results constitute the most detailed analysis yet undertaken of the relationship between shear displacement, shear strength, and fracture transmissivity. RESUME: La rugosité des joints en roche joue un rôle primaire sur la réponse hydromécanique du massif fracturé. En effet, la conductivité hydraulique des discontinuités en roche est influencée par la rugosité, et par la distribution des ouvertures et des zones de contact, dont chacun dépend de l’effort normal et du cisaillement qui agissent sur la surface du joint. Les essais d'écoulement en laboratoire montrent que si cisaillée, la fracture devient anisotrope, et plus perméable dans la direction parallèle au déplacement (Yeo et al., 1998). D'autre part, les essais de cisaillement indiquent que seules les zones de la surface qui font face à la direction de cisaillement sont impliquées dans le cisaillement (Grasselli et Egger, 2003). D'ailleurs, les résultats expérimentaux confirment que la résistance au cisaillement des joints de roche est dépendante de la direction. Dans cet article, plusieurs fractures en roche ont été analysées en utilisant un système de mesure topographique à l'Université de Bâle. La méthode d'analyse de données proposée par Grasselli (Grasselli et al., 2002a) a été employée pour calculer les paramètres de la surface qui déterminent la résistance au cisaillement. Plusieurs simulations du coulage, dans diverses directions relativement à la direction du cisaillement, ont été faites à l’aide d’un programme de calcul à éléments finis. Les résultats constituent une des analyses les plus détaillé sur le rapport entre le déplacement pendant le cisaillement, la résistance au cisaillement, et la transmissivité spécifique de la fracture. ZUSAMMENFASSUMG: Die Festigkeit, Verformbarkeit und Hydromechanik von geklüftetem Fels wird hauptsächlich von der Oberflächenrauhigkeit der Bruchflächen bestimmt. Die hydraulische Leitfähigkeit einer Kluft hängt von der Verteilung der Rissbreiten, der Oberflächenrauhigkeit und den Kontaktflächen ab, die ihrerseits von dem auf die Bruchfläche wirkenden Druck beeinflusst werden. Obgleich in letzter Zeit viel in dieser Richtung geforscht wurde, kann der Zusammenhang zwischen Porenraumgeometrie, Kontaktflächen, wirkendem Druck und hydraulischer Leitfähigkeit noch nicht vollständig erklärt werden. Obwohl Klüfte in situ sowohl Normal- als auch Schubspannungen ausgesetzt sind, wurde der Einfluss von Scherverschiebungen auf den Kluftwasserfluss und die daraus resultierende Fliessanisotropie in einer einzelnen Kluft bisher kaum untersucht. Experimentelle Untersuchungen des Durchflusses zeigen, dass die Kluft für größer werdende Scherverschiebungen heterogen und anisotrop wird wobei die Permeabilität parallel zur Verschiebungsrichtung ansteigt (Yeo et al., 1998). Andererseits zeigen trockene Scherversuche, dass die Oberfläche an Stellen versagt, die gegen die Verschiebungsrichtung anstehen (Grasselli und Egger, 2003). Des weiteren bestätigen die experimentellen Ergebnisse eine Richtungsabhängigkeit der Kluftscherfestigkeit. In diesem Artikel wird eine dreidimensionale morphologische Beschreibung mehrerer natürlicher, durch Scherspannungen verschobener Klüfte vorgestellt, die vom Automated Topographic System der Universität von Basel erstellt wurde. Das von Grasselli (Grasselli et al., 2002a) entwickelte Datenanalyseverfahren wird zur Feststellung der die Scherfestigkeit bestimmenden Oberflächenparameter herangezogen. Mittels einer FEM-Simulation wird der Durchfluss durch eine solche Kluft in verschiedenen Richtungen simuliert. Die Ergebnisse stellen die bisher ausführlichste Analyse des Zusammenhangs zwischen Scherverschiebung, Scherfestigkeit und Durchlässigkeit (Transmissivität) einer Kluft dar.

Jan 1, 2003

By K. L. Morton

Groundwater and associated pore pressure represent elements that can have a negative impact on surface mining and slope stability. Groundwater flow and deformation within a slope influence each other. An increase in pore pressure will result in a decrease of effective stress; and conversely, if the pore pressure decreases, the effective stress increases. Slope analysis and slope design calculations require good information on pore pressure distribution around an open pit excavation. Until recently very little pore pressure data were available and assumptions had to be used. Location and design of measuring points depend on geological and structural setting, geotechnical domains and practicality. Good planning and the use of sensitive point piezometers to measure in situ pore pressures means that real-time data can be collected and used for the construction of flow nets, the determination of pore pressure distribution around an open pit mine, the validation of the dewatering/depressurization requirements, the calculation of the efficiency of the dewatering system and finally the use of pore pressure data in slope stability analysis. The recommended layout for a monitoring system, types and design of piezometers, data collection and use of the pore pressure data are presented.

Jan 1, 2008

By H. B. Arntzen, P. O. Grøntvedt, T. Lindstad, R. J. Ishak, M. Syvertsen

Better utilisation of electric power and coke is vital in reducing production costs and CO2-emissions in production of ferromanganese. All FeMn smelting operations consume 10 to 30 % more carbon than the carbon needed to reduce MnO and FeO to metal and the carbon contained in the alloy. Potassium which re-circulates and builds up in the furnace has been shown to be the main culprit causing both excessive carbon consumption, by catalysing the Boudouard reaction, and increased electric power consumption. The reaction between CO2 and ordinary metallurgical coke or coke impregnated with potassium has been studied at temperatures from 800 to 1100°C in a thermobalance apparatus. Ordinary coke did not react to any great extent at temperatures up to 1000°C, whereas the addition of potassium (2-10 %) resulted in a huge increase in reaction rates. The ultimate goal is to fit the results into a reaction rate model containing the equilibrium constant for the solid-gas heterogeneous reaction between CO2 and a free active site on the carbon reactant producing CO and adsorbed oxygen (an occupied site). The model includes the partial pressures of CO and CO2, a reaction rate constant that is a function of both the total number of reaction sites of carbon, and the potassium content. The potassium dependency is taken care of by introducing activity coefficients for the carbon sites. As the first step towards such a model, a regression line for the Boudouard reaction rate has been established.

Jan 1, 2004

By S. E. Olsen, T. Færden, G. Tranell, T. Lindstad, J. Lubetsky

The authors are lead authors for chapter 4.3 ?Emissions from Ferroalloy Production? in ?2006 IPCC Guide-lines for National Greenhouse Gas Inventories, Volume 3?. The paper discusses the methods to calculate and report emissions from the production of ferrosilicon, silicon-metal, ferromanganese, silicomanganese and ferrochromium alloys, as well as the basis for the methods and choice of emission factors. CO2 is the main greenhouse gas emission from ferroalloy production, accounting for about 95 % of the total CO2-equivalents. The most accurate method of determining emissions is to base the calculations on reported amounts and analyses of reducing agents and carbonate ores used in the process. If analyses are not available, typical CO2 emission factors for the different reducing agents are recommended. If neither reliable amounts nor analyses of all reducing agents are known, default generic CO2 emission factors for the type of ferroalloy multiplied with amount produced, can be used. The CO2 produced from bio-carbon (charcoal and wood chips) should be reported but will not be accounted for in national emissions. If any bio-carbon is used, the default generic CO2 emission factors cannot be applied. CH4 accounts for less than 1 % of the total emission of CO2-equivalents. Default emissions factors for CH4 are given for the production in open or semi-closed furnaces of silicon-metal and ferrosilicon. In the production of FeMn, SiMn and FeCr in closed furnaces the gas containing CO,CO2 and volatiles can be combusted in boilers and N2O may be formed. While reporting is not required, IPCC Guidelines recommend that the N2O-content in the off-gas is measured.

Jan 1, 2007

By R. Ribacchi, P. Tommasi, C. Consorti, P. Campedel

The runout of dry blocky rock avalanches produced by planar rockslides affecting a limestone formation with clayey interbeds is analysed by means of distinct element modelling. Numerical analyses are validated by means of the back analysis of a historical rockslide occurred in the investigated area and physical aspects of the phenomenon evidenced by modelling are discussed. On analyse par la méthode aux élements distincts la distance de parcours des avalanches rocheuses sèches produites par des glissements plan interessant une formation de calcaire à intercalaires argilleux. Les analyses numeriques sont validées à l’aide de la back-analysis d’un glissement historique survenue dans le site en examen. On décrit les aspects physiques su phénomène mis en évidence par la modélisation efféctuée. In diesem Beitrag wird die Reichweite von Steinlawinen besthend aus trockenem blockhaftem Fels mittels der Distinct Element Methode untersucht. Die Steinlawinen werden erzeugt von Rutschungen in einer Kalkstein-Formation mit tonhaltige Schichten betreffen. Die numerischen Berechnungen wurden anhand einer Simulation einer historischen Bergsturzen verifiziert, der sich auf den untersuchten Gebiet ereignete. Auf der Basis der erhaltenen Ergebnisse wurden die physikalischen Aspekte der untersuchten Steinlawinen analysiert.

Jan 1, 2003

By W. Biermann

All industrial processes are economically driven, resulting in the ever increasing need for process optimization in order to stay competitive and profitable. Numerous models has been published depicting the optimization of the raw material usage in an arc furnace, as well as the thermal and electrical efficiency; however, to obtain a holistic view of the cost-effectiveness of the process one needs to expand the model to include the cost of all the operational variables. This project therefore aimed to not only include the raw material feeds and power consumption as input variables, but also the labour, maintenance, and environmental costs to generate an all-inclusive financial balance. In the end, the model was able to mathematically optimize the ferrochrome (FeCr) production process through evaluating the effect of a specified combination of input parameters. It could therefore be used as a tool to aid in managerial decision-making; to appraise and correspondingly adjust certain operational variables in order to achieve a maximum return on investment. This paper fully describes the reasoning behind the choice of the essential model components, the assumptions made during the models development, the calculation approach, as well as the advantages, limitations, possible applications, and areas for expansion of the final model. Through running simulations with plant data, it was found that the accuracy of final model remained fairly constant when applied to stable furnace conditions, but was met with a number of limitations when applied to more complex situations. These limitations could, however, be overcome through further development and expansion of the existing model to include a larger number of process variables.

Jan 1, 2012

By A. Bouvard, N. Hoteit, E. Boidy, F. Pellet

In the context of the construction of an Underground Research Laboratory devoted to the feasibility studies of radioactive waste deposits, it appears of the utmost importance to study the long-term behaviour of rock, over a long time period. In this paper, some results of tests performed on the argillite from the Meuse-Haute-Marne site test in Eastern France are presented. Then a numerical simulation of one of this test using the viscoplastic damaging constitutive law proposed by Lemaitre is presented. Finally, this law is used to simulate numerically the behaviour of the rock at the wall of a circular tunnel in order to determine the extension and the change along the time of the Excavated Damaged Zone (EDZ). Advantages and limitations of the model are discussed. Dans le cadre de la réalisation d’un laboratoire de recherche souterrain conçu pour étudier la faisabilité d’installations de stockage de déchets radioactifs, il apparaît fondamental de caractériser non seulement le comportement viscoplastique à très long terme des matériaux mais aussi les phénomènes d’endommagement consécutif à la redistribution de contraintes provoquées par l’excavation des cavités. On présente ici les résultats d’essais réalisés sur les argilites du site de Meuse/Haute-Marne et destinés à caractériser leur comportement visco-endommageable. Il est ensuite procédé à une application de la loi viscoplastique endommageable de Lemaitre pour simuler l'un de ces tests. Les apports et les limites d'un tel modèle sont mis en évidence par des simulations numériques du comportement de la roche au pourtour d'une galerie circulaire permettant de suivre l’évolution au cours du temps de la zone endommagée (EDZ) à la paroi . Im Rahmen des Baus eines unterirdischen Forschungslabors für die Untersuchung der Machbarkeit von Deponien für radioaktive Abfälle erweist es sich als äußerst wichtig, das langfristige Verhalten von Gestein über einen langen Zeitraum zu untersuchen. Im vorliegenden Artikel werden einige Ergebnisse von Untersuchungen vorgestellt, die mit Tonschiefer aus der Meuse-Haute-Marne in Frankreich durchgeführt wurden. Anschließend wird eine numerische Simulation der Untersuchungen unter Verwendung des von Lemaitre vorgeschlagenen viskoplastischen Verformungsgesetzes vorgestellt. Schließlich wird dieses Gesetz auf die numerische Simulation eines Kreistunnels angewandt, um die Ausdehnung und die zeitliche Entwicklung der Auflockerungszone (EDZ) zu bestimmen. Es werden die Vorteile und Grenzen des Modells diskutiert.

Jan 1, 2003

By T. Saito, S. Murata

An original concept of JRC-diagram that can evaluate JRC and its size-effect by using a fractal model was proposed. In order to confirm the applicability of the JRC-diagram, first, JRC of the tensile fracture surfaces of a granite and a sandstone were evaluated by both the JRC-diagram and Barton's JRC-profiles, and then, it was compared with the JRC inversely calculated from the peak shear strength that is obtained from the shear test for each fracture using Barton's empirical equation. As a result, JRC evaluated by the JRC-profiles was much smaller than the inversely calculated JRC. On the other hand, it became clear that JRC evaluated by the JRC-diagram is almost the same with the inversely calculated JRC, if the size of interlocking is evaluated correctly. Therefore, for the JRC evaluation, it is thought that the JRC-diagram is more applicable than the JRC-profiles. Als ein neues Verfahren zur Bewertung des Rauheitsindexes JRC für die Bruchfläche von Felsenschichten anhand eines fraktalen Modells, wobei auch der Maßeffekt berücksichtigt wird, wird ein JRC-Diagramm vorgeschlagen. Um die Angemessenheit der Bewertung vom JRC mit diesem JRC-Diagramm festzustellen, wurden die JRC-Werte an einer Bruchfläche zwischen Granit und Sandstein sowohl mit dem JRC-Diagramm als auch mit dem herkömmlichen JRC-Profilverfahren bewertet. Anschließend wurden die durch die beiden Verfahren ermittelten Werten mit den JRC-Werten verglichen, die von den Schertestergebnissen mit Hilfe von Bartonschen Erfahrungsformel zurückgerechnet werden können. Daraus ergab sich, daß das JRC-Profilverfahren ziemlich kleinere JRC-Werte gibt als bei der Zurückrechnung, während die JRC-Werte durch das JRC-Diagramm unter entsprechender Bewertung des Ineinandergreifungsmaßes an der Bruchfläche beinahe gleich sind wie die zurückgerechneten Werte. Diese Tatsache weist darauf hin, daß das JRC-Diagramm eine bessere JRC-Bewertungsmethode als das JRC-Profilverfahren darstellt. Nous avons proposé le diagramme JRC, nouvelle méthode d'évaluer JRC qui est un indice de rugosité superficielle d'une fracture de roche à l'aide du modèle de fracturation, compte tenu de son effet d'échelle. Pour vérifier la validité de JRC évalué au moyen de ce diagramme, nous avons évalué JRC d'une fracture de roche, granite et grès, par le diagramme JRC et le profil JRC, méthode conventionnelle, et les résultats ainsi obtenus ont été comparés avec la valeur de JRC calculée en sens inverse à l'aide de la formule empirique Barton sur la base des résultats de l'essai de cisaillement des fractures. Il en résulte que la valeur de JRC obtenue par le profil JRC était assez faible en comparaison avec la valeur calculée en sens inverse alors que celle de JRC obtenue par le diagramme JRC correspondait approximativement à la valeur calculée en sens inverse à condition d'évaluer l'échelle d'endentement d'une fracture. Par conséquent, nous considérons que le diagramme JRC est une meilleure méthode d'évaluation de JRC par rapport au profil JRC.

Jan 1, 2003

By J. Werfling, J. R. Kiehl

SUMMARY: The effects of blasting on structures are assessed on the basis of the maximum amplitudes of the occuring vibration velocities. If the standard values in the German Standard DIN 4150 are not exceeded, damages to structures are not to be expected. The estimation of the maximum vibration amplitudes is usually based on experience or on emperical formulae. However, a comparison of expected with measured values shows that these estimations are not conservative, particulary in case of small distances between blasting locations and structures. More reliable predictions are possible on the basis of numerical analyses, which simulate the three-dimensional propagation of waves with consideration of the geomechanical and geometrical conditions. In addition to the expected velocities, these analyses also allow the determination of stresses resulting from dynamic loading. If the dynamic stresses are known, the effects of blasting on structures can be predicted. ZUSAMMENFASSUNG: Die Auswirkungen sprengbedingter Erschütterungen auf bauliche Anlagen werden in der Praxis auf der Grundlage der Maximalamplituden der auftretenden Schwinggeschwindigkeiten beurteilt. In der DIN 4150 werden Anhaltswerte angegeben, bei deren Einhaltung Schäden von Gebäuden nicht zu erwarten sind. Die Abschätzung dieser Schwing-geschwindigkeiten erfolgt meist auf der Grundlage von Erfahrungswerten oder empirischer Formeln. Der Vergleich der prognostizierten Werte mit Meßergebnissen zeigt jedoch, daß diese Ansätze insbesondere dann, wenn die Abstände zwischen den Sprengorten und den Gebäuden gering sind, nicht konservativ sind. Zuverlässigere Prognosen sind auf der Grundlage von numerischen Berechnungen möglich, in denen die räumliche Ausbreitung von Erschütterungen unter Berücksichtigung der geomechanischen und geometrischen Verhältnisse simuliert wird. Mit derartigen Berechnungen können die zu erwartenden Schwinggeschwindigkeiten und die aus der dynamischen Beanspruchung resultierenden Spannungen ermittelt werden. Dies ermöglicht Aussagen über die Auswirkungen sprengbedingter Erschütterungen auf Bauwerke. RÉSUMÉ: Les conséquences de vibrations résultantes d'explosions pour des bâtiments sont estimés dans la pratique à partir des amplitudes maximales des vitesses d'oscillation. Dans la norme allemande DIN 4150 des valeurs d'orientation sont données, avec l'observation desquelles de dommage à des bâtiments n'est pas à présumer. Ces vitesses d'oscillation sont estimés généralement à partir de valeurs d'expérience ou de formules empiriques. La comparaison des valeurs prédits avec des résultats de mesures montre cependant que ces hypothèses ne sont pas conservateurs, en particulier si les distances entre les lieux de sautage et les bâtiments sont petites. Des prédictions plus fiables sont possibles à partir de calculs numériques, dans lesquelles la dissipation spatiale de vibrations est simulée en considération de la situation géoméchanique et géométrique. Avec tels calculs les vitesses de vibration à présumer et les contraintes résultantes de l'effort dynamique peuvent être établies. Cela permet de constater sur les conséquences de vibrations résultantes d'explosions pour des bâtiments.

Jan 1, 2003

By S. Sarac

In this study, a statistical model is suggested in order to determine the methane emission from a working seam. In the proposed method, methane concentrations which are systematically measured in return airway are the sample values, and the calculated statistical parameter is the upper confidence limit of the average of these values. At any given time, the probability of methane concentration in the return airflow to be lower than the calculated upper confidence limit value is 97.5%. The statistical method has been applied to the panel 206 of Armutcuk Mine, Turkish Hardcoal Enterprises. The calculated methane emission value for the considered shift was 15.18 m3/t.

Jan 1, 2003

Changes in reductant and energy inputs which are required to adjust the slag composition in an ilmenite smelter to a higher or lower FeO content are calculated. The calculated changes are based on the assumption that solidification equilibrium with the freeze lining constrains the slag composition to follow a eutectic groove in the TiO2-FeO-TiO1.5phase diagram. Calculated equilibrium conditions along the eutectic groove are presented, and form the basis of the calculations. The inputs of energy and carbon need to deviate strongly from the steady-state values to adjust the slag composition. The inputs required at any given time depend on the desired rate of change of the slag composition, and the current slag composition. The steady-state relationship of a change in energy input of 5 kWh corresponding to a 1 kg change in carbon input also holds for the changed inputs during transient changes in FeO content. Keywords: Ilmenite smelting, freeze lining, solidification, eutectic

Jan 1, 2003

Caving of the second lift of the E26 orebody at Northparkes Mines was monitored with a comprehensive package of instrumentation, including an ISS seismic monitoring system. Seismicity was collected by an array of 19 triaxial accelerometers and 9 uniaxial geophones, giving a high quality seismic data set. The seismic response to block caving at Northparkes has been divided into two papers. The first paper focuses on the seismicity during undercutting. The second paper deals primarily with the seismicity during cave propagation during initial cave production. In this paper, traditional seismic monitoring concepts are applied to seismicity during cave mining. More than 40 000 seismic events were recorded during cave undercutting between February and January 2004. Events up to local magnitude +1.3 were recorded near the mine during undercutting, with event rates of up to 500 per day. Analysis of frequency-magnitude distributions, S-wave to P-wave energy, and temporal variations in event magnitude all gave insight into the seismicity recorded. In addition, a new technique, called apparent stress time history, was shown to be sensitive to stress change from undercut blasting, and useful to identify significant periods during cave development. A significant difference was identified between the cave back location and the halo of seismic events above the cave. Throughout the undercutting, the a seismic zone of loosening was typically 50 to 70 metres in vertical thickness. Based on the seismicity, a stress driven rock mass failure model is proposed. Early in the undercutting, seismicity and movements in the seismogenic zone were directly triggered by stress change from undercutting. Towards the end of undercutting, seismicity was virtually uncorrelated with undercut blasting, and likely due to the onset of cave initiation.

Jan 1, 2008

By D. Limpitlaw

The origins of mine closure practice have influenced the way in which it is implemented by companies and regulators. Mine closure practices essentially started developing in the 1970s in countries with advanced economies and mature mining industries. In these settings, the emphasis was justifiably placed on restoration of the landscape and an attempt was made to return to the ?natural? pre-mining land cover. These practices continued to evolve and incorporated socio-economic and cultural aspects, especially after the Brundtland Report in 1987 and the subsequent Earth Summit in 1992. Today mining is increasingly occurring in remote parts of developing countries where there may be significant need for infrastructure such as roads, clinics, and schools. The costs of returning land to low (economic) value pre-mining use may be far greater than establishing a viable postmining land use that could potentially add value to the community and take pressure off sites for greenfield development elsewhere. Furthermore, natural resource limitations (such as topsoil availability) may limit the degree to which the historical land cover can be re-established. Establishing post-mining land uses may aid in mitigating the loss of employment that is inevitable when mines close. Stakeholder participation in establishing post-mining land cover and land use options is critical for long-term success. Similarly, third parties must be identified to support the development of the post-mining land use This paper draws on the experience of the authors in several developing countries and presents a case for maximizing re-use of mining infrastructure. The paper does not advocate the adoption of poor rehabilitation standards, nor or the wholesale destruction of land capability, but rather leaving key infrastructure in place for post-mining use that may support sustainable development.

Jan 1, 2014

By D. F. Malan

This paper presents an overview of work conducted over several years to investigate the continuous stope closure profiles recorded in tabular excavations and the use of this data as a hazard indicator. In spite of the brittle nature of the rock, measurements in gold mines and intermediate depth platinum mines revealed a significant time-dependent closure component. Various approaches to simulate the time-dependent closure have been investigated. Viscoelastic theory is able to simulate the time-dependent closure behaviour, but there are subtle problems with the use of this theory as it cannot simulate the failure processes in the rock. To overcome these problems, continuum and discontinuum elasto-viscoplastic models were developed by the authors. Discontinuum models have the advantage of being able to simulate the creep behaviour of prominent discontinuities. Development of these models has led to useful insights into rock mass deformation and has made the numerical simulation of parameters, such as rate of mining, possible for the first time. A number of practical applications for the mining industry have been developed from this study. Among these are the delineation of geotechnical areas using closure data, use of the data as a hazard indicator and to give warning of large collapses in the platinum mines, identification of areas to be preconditioned in the gold mines, measurement of the effectiveness of preconditioning and improved da

Jan 1, 2007

By A. Bothma

The Impala No. 1A ventilation shaft is located approximately 15 kms north of Rustenburg and is accessed via the R510road from Rustenburg to Northam/Thabazimbi. This situation is shown on the location plan in Figure 1. The geology at the site consists of an upper layer of black turf of about 1 metre in thickness, grading into weathered norite. The depth of weathering is generally less than 5 metres. Zones of deeper weathering may occur along fracture zones associated with faults and dykes. Groundwater occurs mainly within the weathered/fractured aquifer system. The ventilation shaft is a 5.1 m diameter upcast shaft ,and forms part of a system of ventilation shafts to complement the existing infrastructure as mining progresses down the new twin decline extension of Impala 1 shaft. The principal reason for these shafts is to lessen the risk of an extended RAW and also to counter the effects of the heat build-up as the mining goes deeper, with the resultant effect that the existing cooling centre (fridge plant) becomes more remote from activities. The No. 1A shaft (upcast) is positioned approximately 1,2 km on down dip from the main shaft and the No. 1Bshaft (downcast) approximately 1,2 km further. Geological information was used from a borehole 1,2 km away on strike from the upcast shaft. The first shaft (1A) is a 5,1 m diameter hole to be fitted with two 2,2 MW fans, whilst the second shaft (1B) is 4,1 m diameter to accommodate a 7MW fridge plant.

Jan 1, 2001

By Dolly Chakraborty, S. Ranganathan, S. N. Sinha

Pre-reduction of chromite outside the submerged-arc furnace (SAF) can significantly lower the specific power consumption in the SAF process. Conventional processes for the pre-reduction of the ore operate at relatively high temperatures, at 1400-1500°C. It is possible to enhance the rate of reduction at low temperatures by optimizing the flow rate of inert gas over the reaction site. The influence of the flow rate of inert gas on the rate of reduction of chromite ore was investigated. The rate of reduction initially decreased wit the flow rate of the gas. It increased subsequently at higher flow rates to be followed by a decrease again. There was an optimum flow rate at which maximum reduction could be obtained. It was possible to achieve high levels of reduction at 1200°C by optimising the flow rate of the gas. The degree of reduction obtained was comparable to that obtained at 1400-1500°C, by conventional process routes.

Jan 1, 2007

By J. H. Cloete, S. M. Bradshaw, D. K. Chibwe

"Modern tundishes have evolved as vessels to serve as the final step in refining of molten steel by removing inclusions and promoting thermochemical homogeneity.In this study the flow behaviour in a four-strand tundish was investigated by means of a ½-scale water model as well as numerical modelling. The numerical and physical models were used to characterize residence time distribution and calculate properties pertaining to tundish flow regime. Three different tundish configurations were investigated: a bare tundish with no flow control devices, a tundish with a turbulence inhibitor, and a tundish with both a turbulence inhibitor and a dam.The physical and numerical models showed that a tundish without flow control devices is prone to significant short-circuiting. A tundish with a turbulence inhibitor was shown to be successful in preventing shortcircuiting and provided surface-directed flow that might assist the removal of inclusions from the melt. However, it was also observed that the upward-directed flow caused the maximum turbulence kinetic energy near the surface to increase dramatically. The potential for slag entrainment should therefore be considered during the design and operation of tundishes with turbulence inhibitors."

Jan 1, 2015

By D. J. Fourie

During the production of pure TiO2for the pigment industry, ilmenite, containing 35?60% TiO2, is reduced to high titania slag, containing 70?95% ?TiO2? and pig iron. The high titania slag produced, consists mainly of titanium in different oxidation states and FeO in equilibrium with liquid iron. To predict the phase and chemical reaction equilibria in a furnace, a suitable solution model is required with accurately determined binary interaction parameters, as slag melts tend to behave very non-ideally. This paper investigates two popular slag models, namely, modified quasi-Chemical (MQC) model and the cell model and provides the binary interaction parameters based on recently published data of the ternary system and critically assessed published unary and binary data. Very few modelling results have been published on the ternary Ti-Fe-O system pertaining to high titanium slags. The solution models, with their regressed interaction parameters, can subsequently be implemented in commercial software of the slag-metal phase equilibria to predict how much and which type of solids will precipitate for a given bulk slag chemistry and operating temperature or, alternatively, above which temperature should a furnace be operated to avoid excessive solids precipitation.

Jan 1, 2005

By K. Liddell

Current practice for the recovery of platinum group metals (PGM) from flotation concentrates is by energy-intensive smelting processes followed by separation and refining of the base and precious metals. Energy costs are rapidly increasing in South Africa and this trend is forecast to continue in future years. Moreover, the smelting operations are highly susceptible to deleterious constituents, particularly chromite; and the increasing tonnages of UG2 chromititic ore being mined results in an increasing chromite content of smelter feeds. Hence there is a substantial incentive to assess and implement alternative processes to smelting for the treatment of PGM concentrates. The Kell process is a patented hydrometallurgical alternative to smelting that is totally insensitive to the chromite content of the feed concentrate and comprises three commercially proven unit operations (pressure oxidation leaching of base metals and sulfur, roasting and chlorination of precious metals). The Kell process provides the opportunity for existing and new PGM producers to significantly reduce their energy consumption, and in particular their consumption of electricity, by removing the need to smelt PGM concentrates. Energy costs, greenhouse gas emissions and installed power requirements are also much reduced. Compared with smelting and refining these reductions are calculated as: ? Total energy consumption-50% reduction ? Electrical energy consumption-84% reduction ? Energy consumption costs-76% reduction ? C02 emissions-70% reduction ? Installed power requirement 92% reduction. The data obtained from simple laboratory batch amenability testing can be readily applied to the Simulus mass and energy balance simulation procedure to determine energy and reagent consumptions, equipment selection and process optimization. The applicability of Kell to treating low-grade 'dirty' concentrates with high UG2 content as well as to polymetallic concentrates containing base metals as the major payable metals with subordinate PGM means that these savings can also be expected when the process is used for treating a range of concentrate types and grades. Keywords PGM,UG2, Merensky, smelting, chromite, electricity, hydrometallurgy, Kell Process, Eskom.

Jan 1, 2011

By M. W. Hildyard

A blast was engineered close to a tunnel in a deep-level mine with the purpose of studying the wave interaction with the tunnel. Numerical modelling of seismic wave propagation was used in both the forward planning and the back-analysis of this experiment. This paper presents the forward analysis. A likely model for a propagating blast was investigated and developed. This model predicts the influence of the velocity of detonation of the blast and the rise-time of the blast on the radiation pattern, and hence on the expected positions of maximum velocity at the tunnel surface. An attempt was made to quantify the parameters of the blast model, using waveform data from two other blast experiments. The blast model caused particle velocities which were an order of magnitude lower than these measured velocities. Large values for blast pressure were required to match the amplitudes of the recorded velocities. An alternative ?effective? source with a similar mechanism but increased diameter was proposed to account for the higher velocities. It was established that distant velocities scale with the square of the borehole diameter. A calibration blast at the experimental site was also modelled. Agreement in the waveforms was not sufficient for conclusive statements to be made about the validity of the source model. However, the modelled waveforms show an encouraging correspondence with some measured waveforms. This gave confidence that the model could be used to examine wave propagation around the tunnel, and allowed a scaled source to be developed for forward analysis of the main experiment.

Jan 1, 2001