Search Documents

By J Finlay, P Betts, D Milton, H Williams

The aim of this paper is to present an innovative method to help delineate potential target volumes for haematite-goethite mineralisation as opposed to more magnetic ores. This is achieved through interpretation and modelling (2.5D and 3D inversions) of ground magnetics and gravity surveys over the Koolanooka South prospect in Western Australia. The targets are assumed to be dense and weakly magnetic and hosted within banded iron formation (BIF) horizons. Geophysical interpretation (2D map) was performed using a variety of processed magnetic and gravity images (Betts et al, 2007) to determine the distribution of lithologies, to understand the structural architecture and to map out any anomalous areas of excess density (possibly indicating goethite alteration). The created geological map is used to build geological cross-sections (geometrically consistent) that were modelled (in 2.5 D) against both potential field datasets. The crosssections are consequently used as the framework to build a 3D geological model. The model defines 3D volumes of rocks including several BIF horizons. The distribution of densities and magnetic susceptibilities of the BIFs is inverted within the GOCAD«-VPmg environment to reduce the misfit between the response of the model and the observed gravity and magnetics. By combining the results of both petrophysical property distribution it is possible to map out (using classical GIS type queries in 3D) volumes of higher density and lower magnetic susceptibility.

Jan 1, 2009

By B. D. Nichols, W. S. Gregory

SOLA-DM is a three-dimensional finite-difference computer code designed to model the dynamics of an incompressible fluid and the transport of discrete particulate material around obstacles impervious to flow. The numerical methods used in this code are described. SOLA-DM was used to predict the particle flux sampled by the 10-mm Dorr-Oliver Cyclone and MINIRAM dust monitors. Various geometric and dynamic variations of monitor and airflow combinations were tested. The code predictions are shown in computer-generated graphic plots.

Jan 1, 1987

By Tarek M. Hatem, Salah M. Bedair, Khaled H. Khafagy

Our understanding for intrinsic stresses and defects evolution in photovoltaic devices has became an essential part of new developments. In particular, Multi-Junction Photovoltaic (MJ-PV) modules depend on multi-layer structures that may suffer high dislocation-densities as a result of high lattice and thermal expansion coefficient mismatch. These defects limit the performance, reliability, and lifetime of PV devices. In the current study, a three-dimensional multiple-slip crystal-plasticity model and specialized finite-element formulations are used to investigate InGaN growth on Si substrates. The formulation is based on accounting for thermal and intrinsic stresses as a result of different processing conditions and microstructures. Furthermore, the formulation was used to investigate a recently developed technique, Embedded Void Approach (EVA), which can be used to address both the high density of defects and the cracking/bowing of InGaN growth on Si. The current work lays the groundwork for more extensive use of silicon in MJ-PV devices.

Mar 1, 2018

By Bacigalupo-Rose

Simulations of complex physical phenomena, such as hurricanes and fluid flow around space vehicles, can now be effectively explored using modem visualisation techniques which have been developed to examine the results of numerical modelling using supercomputers. These techniques are also being applied extensively in medical investigations, including viewing the interior of the human brain. We have now applied these same techniques to the visualisation of gold deposits.Grade control in open pit gold mining operations may involve the collection of at least half a million samples in each mine. Extensive drilling programs also accompany the early reserve evaluation stages and exploration in advance of mining. Maximum utilisation of these very large amounts of information requires an ability to integrate, examine and present large three-dimensional data sets in order to be able to derive and effectively communicate information on grade distribution and geological controls.

Jan 1, 1993

By Yue Gong, ChuanJie Cai, Jing Chen, Shaobo Zheng, Huigai Li

"MnS-AlN duplex precipitates were successfully extracted from electrical steel slabs by means of non-aqueous solution electrolysis. The three-dimensional duplex morphology became visible using emission scanning electron microscope (SEM) JSM-6700F and energy dispersive spectroscopy (EDS). The precipitation sequence of MnS and AlN in the hot process can be predicted via the three-dimensional duplex morphology. Thermodynamic calculations were carried out on account of the chemical compositions of the steel studied in this paper. The results verified that AlN nucleated on MnS, which precipitated earlier during cooling process of electrical steel. Based on the results, the material compositions and rolling temperature can be adjusted to modify their precipitation timing and morphology. This is conducive to the growth of MnS-AlN duplex precipitates and reducing the obstacles to grain growth. Finally the magnetic property of electrical steel would be improved.IntroductionElectrical steel is a kind of ferrosilicon soft magnetic alloys with extremely low carbon content. It is mainly used in the manufacture of motors, generators, transformers and varieties of teleco=unications equipment. Non-oriented electrical steel requires isotropic high magnetic induction and low iron loss [l, 2]. Precipitates in steel have a great influence on its performance [3]. On the one hand, the indirect effect of precipitates on magnetic properties is to restrict grain growth in high temperature annealing process. On the other hand, the direct effect is pinning of domain wall motion, which leads to increased hysteresis loss."

Jan 1, 2013

By M. J. Beus

A shaft is truly the "lifelinen of an underground mine. Damage to the shaft lining and guides as a result of ground movement can result in serious loss of production and extensive and continual repair. In some mining districts, annual shaft repair and maintenance costs resulting from high rock stresses and excessive displacement are in the millions of dollars per year. For many years, the Bureau of Mines has conducted research to characterize the rock mass and to develop structural guidelines to improve the design of accessways in deep mines. Much of this work has centered on determining the in situ conditions that affect the structural stability of shafts in the Coeur d'Alene Mining District of northern Idaho. Recent work by Beus and Board (1984) has concentrated on measuring the rock and support behavior during sinking of the Silver Shaft at the Lucky Friday Mine. These results are currently being used as a comparative data base for further basic research as well as to establish structural design criteria. The factors affecting the deign and structural stability of deep mine shafts are numerous. A design approach that use field data as input and directly compares field measurements with numerical models can validate the procedure and provide realistic design criteria. Naturally occurring "fixed? conditions such as magnitude, direction, and ratio of in situ stresses, geologic environment, and physical properties and constitutive relationships of the rock mass ere obviously basic to numerical modeling.

Jan 1, 1989

By C. O. Aksoy

"Pillars left to ensure the stability of major underground coal mine infrastructure (shafts, drifts) often contain substantial quantities of coal. Recovery of this coal is an efficient use of the resource, but often presents safety hazards; pillar extraction destabilizes the protected structure, and can make its use dangerous. In this study, the stability of an active shaft at an underground coal mine situated in the Soma- Manisa region of Turkey was investigated during the recovery of the coal in the shaft safety pillar.A proprietary three-dimensional finite element analysis model (ANSYS) was used. An elastic-plastic material model utilizing the Drucker-Prager failure criterion was employed. The three-dimensional Drucker-Prager formulation was preferred over the two-dimensional Mohr-Coulomb model because it was thought to more accurately represent the strain hardening behaviour shown by most geologic materials.Stresses and deformations formed around the shaft support were estimated using the model and the results were verified by in situ observations. The results show that stress distributions around the shaft and the subsidence of the shaft support are influenced by the location location of the extraction front and the break angle of the overlying material."

Jan 1, 2005

By Wedam Nyaaba, Samuel Frimpong, Bruno Ayaga Kansake, Irene Akalanyaba Ateng

Haul truck capacities have increased due to their economies of scale in large-scale surface mines. These high-capacity trucks impose dynamic loads on haul roads. A previous study estimated the dynamic forces imposed by ultra-large trucks on haul roads to be more than three times the static load. These dynamic forces exacerbate haul road stresses and deformation. Current haul road response models are 2D and use static truckloads for low-capacity trucks. No models capture the ultra-large truck dynamic effects on haul road structural response. This study uses 3D finite element (FE) modeling to compute the haul road response to ultra-large truck dynamic loading. The dynamic forces were modeled using multibody dynamics (MBD) in MSC.ADAMS. These forces were used in an FE model developed, verified, and validated in ABAQUS, to simulate the response of the haul road to the dynamic truck forces. The road was modeled using an elastoplastic Mohr-Coulomb model. The results showed that road deformation decreases with increasing layer modulus and increases with increasing payload. Increasing the base modulus from 100 MPa to 450 MPa reduced the maximum deformation from 258 mm to 62 mm. When the subbase elastic modulus increased from 100 MPa to 500 MPa, the maximum deformation decreased from 160 mm to 84 mm. Increasing the payload from 0 to 120% of the rated load increased the maximum wearing surface deformation from 131 mm to 216 mm. This study forms a basis for designing structurally competent haul roads for improving haulage efficiency and truck and operator health.

Jan 7, 2023

The finite element method of stress analysis is nowadays widely used in the fields of mining, mechanical and civil engineering. This paper presents a short description of the three-dimensional method, in which the elements used are arbitrary hexahedrons. Two examples of the application of the method are presented, namely, calculation of the elastic stress distributions in a coal pillar and in an opencast mine.

By A. L. Holloway, N. M. Soonawala, L. S. Collett

"Since the late 1970s, a major R&D program has been under way in Canada to assess the concept of disposing of nuclear fuel wastes in deep plutonic rock. Atomic Energy of Canada Ltd . (AECL), a federally-owned corporation, has been given the responsibility of directing this research (Rummery et al., 1983).The Underground Research Laboratory (URL), constructed in a granite pluton 15 km east of Lac du Bonnet, Manitoba (Fig. 1) is a major component of the Canadian Nuclear Fuel Waste Management Program. A diverse geotechnical experimental program is being carried out during both the construction and operating phases of the URL. This program involves detailed geological mapping, rock mechanics and hydrogeology studies, and geophysical surveying, to characterize the rock mass surrounding the shaft, and it has been described by Simmons et al. (1983). As part of the construction-phase program, ground-penetrating radar (GPR) surveys are being carried out during excavation, in an attempt to delineate fractures both ahead of the excavation and behind the walls of the excavation. This paper presents the results of surveys performed during shaft excavation and describes a procedure for plotting the data.Ulricksen (1982) cites examples of radar surveys to evaluate surficial deposits, such as peat, sands and gravels, and cavity detection in subsurface hard rock. Fowler (1981) used radar for fracture and clay vein detection in coal seams.In Canada, pioneering work in GPR has been done by Annan and Davis since the mid-seventies (Annan and Davies, 1976). Their work concentrated on problems relating to permafrost in the Arctic ."

Jan 1, 1986

By S. Heusermann, S. Fahland, J. Behlau, A. Schäfers, J. Hammer

In the course of the deep subsurface use, the Federal Institute for Geosciences and Natural Resources has been carrying out extensive geological and geomechanical research and practical project work on domal salt structures, to prove their suitability e.g. for the disposal of radioactive wastes. Recent work on disposal sites has been taking into account the 3-D geological and geomechanical modelling of the underground structures. The geological models are generated using the special-construction openGEOTM code to improve the visualisation and interpretation of the geological data basis, e.g. borehole, mine, and geophysical data. For the geomechanical analysis of underground structures the JIFE finite-element tool has being developed to analyze large 3-D structures and to consider coupled THMC processes. The results of finite-element calculations, i.e. stresses, strains, displacements, and dilatant zones, are basic values for the evaluation of the stability of the structures and the long-term integrity of the geological barrier. As an example, the results of geological and geomechanical investigations of the Morsleben repository are presented including three-dimensional modelling of the entire geological structure of the central part of the mine as well as three-dimensional geomechanical model calculations to provide the necessary results for a long-term safety prediction.

Jan 1, 2015

By R Knight

The spatial variability of geometallurgical characteristics in a mineral deposit can have significant impact on financial and mine planning, mineral processing, concentrate quality and tailings composition throughout the life of a project. These characteristics that may determine the success of a project are not readily apparent to geologists unless they are integrated with deposit geology and assay data, modelled in 3D and incorporated into resource block models.One of the primary objectives of geometallurgy is to integrate and incorporate various geologic, geotechnical and geometallurgical parameters into resource block models, which typically are limited to in situ tonnes and grade of economic metals. Critical attributes, such as concentration of deleterious elements that affect the project viability need to be identified and integrated into the resource block model early in the exploration stage to predict project profitability and address potential environmental impacts.Using data from a volcanogenic massive sulfide deposit we present an example of the integration of ore geochemistry (cadmium, arsenic and antimony) with the results of preliminary, bench-scale metallurgical test work in a 3D geometallurgical model that predicts subdomains of poor-quality concentrate and tailings containing environmentally-sensitive elements. This critical information was incorporated into the existing resource block model in order to inform both financial modelling and mine planning.This example demonstrates that geometallurgical data can be effectively integrated with geological and assay data using Leapfrog MiningTM software (Leapfrog Mining) early in the exploration stage of a project. Three-dimensional representations of geometallurgical data can be created directly using primary additive data, as well as from proxies of non-additive parameters. Resource block models can be updated continually with new information as the project passes from exploration through production. In addition, the software allows for 3D visualisation of the characteristics of the materials that will eventually constitute the mill concentrate, tailings, and waste. This study emphasises the critical importance of establishing a systematic geometallurgical approach early in the exploration stage that can be continuously updated with new information throughout the life of the project. In the example that we present, sulfide ore that had otherwise been identified as high-value (in the absence of geochemical and metallurgical data) turns out to be low-value ore due to concentrate penalties and increased environmental impact.

Sep 5, 2011

By Dan J. Van Roosendaal

A closely spaced grid of survey points was installed over the tensile zone of a longwall panel in southern Illinois and monitored daily during subsidence to document the complex ground-surface movements associated with dynamic subsidence. Vertical and horizontal displacements were used to calculate maximum tilts and principal strains within the grid. The dynamic nature of subsidence is made evident by changes in both the direction and magnitude of maximum tilts and principal strains. The maximum tilt (smax) of 0.036 and maximum tensile strain (Emax) of 0.027 were measured after subsidence and were oriented transverse to the longwall panel. Observations of surface fracturing were correlated with strain determinations to estimate a surface strain at incipient cracking (E;) of 0.006 to 0.009. Using shallow inclinometers, a linear variation of horizontal displacement with respect to depth was observed, suggesting a neutral¬axis bending mechanism for near-surface deformation. Horizontal displacements transverse to the longwall panel reduce to zero at a depth of 40 feet. A linear relationship between strain and curvature was used to estimate ground-surface curvatures associated with Emax (0.00068 ft-1) and E; (0.00015 to 0.00023 ft-1). The curvature estimates are consistent with published results.

Jan 1, 1986

Owing to improvements in computer capacity, it is now possible to undertake fully coupled or partially coupled, three-dimensional (3D) hydromechanical Finite Element (FE) simulation of mines, including explicit representations of large numbers of discontinuities. The analysis requires fewer assumptions than uncoupled analysis or two-dimensional (2D) analysis, but may require a change in how hydrological data is collected. The way hydrological data is incorporated in stability assessment and the approach to uncertainty in hydrological information also requires discussion. The current procedures for this type of analysis are discussed, along with some practical difficulties and teachings from early applications. A work flow for integration of 3D hydromechanical modelling with pit planning and some considerations for the management of uncertain data are also proposed.

May 1, 2010

By C A. Foss, D A. Clark, D Hillan

While magnetic surveys are the most common geophysical method used in mineral exploration, magnetics and geology remain to be fully integrated quantitatively in iron ore provinces. Although high-resolution aeromagnetic data has been collected over prospective areas of the Hamersley Province, much of the information in the data sets has hitherto not been fully exploited due to a lack of suitable interpretation methods. Recent advances in the processing and modelling of magnetics in highly magnetic terrains, constrained inversion and 3D visualisation and understanding the magnetic properties of banded iron formations and iron ores potentially create the capability to generate realistic 3D mapping of magnetite and hematite concentrations over large areas. Initial processing of conventional total magnetic intensity data produces true potential field maps that are unaffected by local perturbations of the regional geomagnetic field direction by intense local anomalies. The true potential field can be converted into vector and tensor components, each of which emphasises different structural features of the geology. Scalar quantities derived from these components simplify interpretation, using new methods, developed by CSIRO and elsewhere, that are less sensitive to magnetisation direction. Preliminary interpretations allow interpolation between mapped units and structures and evaluation of remanence and anisotropy effects. Accurate modelling requires rigourous treatment of self-demagnetisation and magnetic interaction effects. Semi-automated constrained inversion, combining magnetics with available geological information and other geophysical data, can be utilised to generate models of disposition of major magnetic units and the 3D magnetisation distribution within them. The final stage of interpretation is to exploit the understanding of magnetic properties (susceptibility, anisotropy and remanence) developed over many studies by CSIRO, the University of Western Australia and others to convert the interpreted magnetisation distributions into 3D maps of magnetite and hematite contents, which are the quantities of direct exploration interest.CITATION:Clark, D A, Foss, C A, Austin, J and Hillan, D, 2015. Three-dimensional mapping of magnetite and hematite concentrations from reprocessing of detailed aeromagnetic data, in Proceedings Iron Ore 2015, pp 229–240 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Jul 13, 2015

By R. Gong, R. Bhaskar

The development of a three-dimensional mathematical model for evaluating airflow fields at a continuous miner face is described. The model utilizes three-dimensional incompressible flow, mass continuity, and momentum (Navier-Stokes) equations to predict the airflow velocity at the face. Obstacles and equipment in the face area are modeled using appropriate boundary conditions. Parameters such as flow velocity and face dimensions can be modified and evaluated for various schemes, such as exhausting and blowing ventilation, using either brattice or tubing. A model such as the one described, with additions for local effects, can be used to evaluate novel face ventilation ideas before initiating expensive full-scale tests or field studies. The model described is a step in the direction of advanced predictor models.

Jan 1, 1993

By X. Jia

The ability to perform three dimensional imaging at micron-scale resolution of real mineral floc sediments and structures using commercially available Micro X-ray Tomography (XMT) instrumentation can be used for the analysis of aggregates and granules, liberation and dewatering of sediments. This paper presents are view of recent work by the authors relevant to the mining and minerals industries. It covers three main topics: X-ray imaging of agglomerates and multi-component mineralic particles,, permeability of sediments, and dissolution of granules with potential application to heap leaching.

Jan 1, 2008

By Chan. C. L., Zehr. R., M. M. Chen, J. Mazumder

"Two three-dimensional models of a surface tension gradient driven flow during laser surface heating are developed. The first model is based on a perturbation solution. The basic solution corresponds to the stationary axisynwnetric case and perturbation is based on a small scanning velocity. The advantage of seeking a perturbation solution, as it turns out, is that the three-dimensional flow is modeled by two sets of two-dimensional equations which are presumably much more tractable than the original three-dimensional equations. Numerical solutions are obtained and discussed. The second model is a full three-dimensional numerical solution of the Navier-Stokes equations, using a point-by-point partially vectorized iteration scheme. Surface shape is also determined in a self-consistent manner. The effect of the presence of convection on pool geometry, cooling rate and solute redistribution is presented and discussed.IntroductionThe existence of fluid flow and convection and their effects were first reported in the fifties (1,2). The role of convection in controlling pool geometry such as pool shape, undercut and ripples, and defects such as variable penetration, porosity and lack of fusion and weld properties such as weld homogeneity and bead formation is studied by a number of researchers and is reviewed by the authors in a recent paper (3).Recently, surface tension gradient driven flow has been identified as one of the forces responsible for convection within the molten pool (4,5). It is postulated that surface tension gradient driven flow is responsible for the surface ripple formation in the weld pool (6-8). Oreper, et al. (9) developed a two-dimensional convection model. Buoyancy, electromagnetic and surface tension gradient forces were considered. The formulation was based on a specified pool profile. It was found that the surface tension gradient is the dominant force in most cases. However, the solid-liquid interface is not known a priori. In fact, the solid-liquid interface, i.e. the pool shape, is a piece of information that the model should be able to predict. A two-dimensional transient self-consistent (i.e. the solid-liquid interface is solved as part of the solution] model for laser surface melted pool was first developed by Chan, Mazumder, and Chen (3)). This model predicted a non-uniform cooling rate, hence non-uniform microstructure. Depending on the Prandtl number (i.e. material) and surface tension number, the shape of the molten region as well as the cooling rate changes. It was also predicted that the recirculating velocity is of one or two order of magnitude higher than that of the scanning velocity. This model provides the detail of the flow field and heat transfer on a plane perpendicular to the scanning direction. However, the front-to-back motion and its effect on heat transfer also play important roles in the physical process, especially for the fact that the recirculating velocity is so much higher than the scanning velocity. Such a motion cannot be obtained from two-dimensional models."

Jan 1, 1986

By J. T. Birkholzer, M. Rutenberg, K. -H. Lux, R. Wolters, J. Rutqvist, L. Blanco Martin

In this paper, we present numerical results corresponding to the Thermal Simulation for Drift Emplacement (TSDE) test, conducted in the Asse salt mine (Germany) between 1990 and 1999 to simulate reference repository conditions for spent nuclear fuel in rock salt. This test is of significant relevance for nuclear waste disposal in salt because extensive research and measurement programs were set up for the experiment. In our numerical study, we use two sequential simulators, TOUGH-FLAC and FLAC-TOUGH, and we investigate two-way coupled processes, focusing on thermo-mechanical analyses due to the main effect of heating. The two simulators used have recently been updated and are able to deal with large strains and creep processes. Moreover, specific constitutive relationships and coupling functions for rock salt and crushed salt under repository conditions are available. Within the framework of a benchmark exercise between Lawrence Berkeley National Laboratory and Clausthal University of Technology, the same scenario has been modeled in three dimensions using the two simulators mentioned. The main objective of the current, first stage of the benchmark is to evaluate the capabilities of these two simulators to predict relevant processes under repository conditions for high-level nuclear waste in saliniferous media. The comparison of the numerical predictions between them is very satisfactory, leading to a successful code-to-code verification. The comparison between measurements and numerical results shows that the most important processes are captured by the codes and constitutive models used. However, some differences exist between measurements and numerical results. These differences will be addressed in the next stage of the benchmark exercise, in which a parameter fitting effort will be undertaken. Overall, from the results obtained it can be concluded that TOUGH-FLAC and FLAC-TOUGH have the necessary capabilities to reproduce relevant processes under repository conditions for high-level nuclear waste in saliniferous formations.

Jan 1, 2015

By K. V. Raj, T. Bhowmick

"Control of fugitive dust in high latitude open-pit mines is challenging due to unique atmospheric phenomenon resulting in complicated flow regimes, and atmospheric inversion due to lack of adequate solar insolation during prolonged winter seasons. In this study, two idealized (One conical and one trapezoidal) three-dimensional open-pit mine geometries are simulated for differing seasonal conditions using the CRADLE-CFD package. The airflow is modeled using Standard Kappa-Epsilon Turbulence Model solved by RANS equations. The concept of particle tracking is used to predict the flow pattern of dust particles. For various climatic conditions, and two different pit geometries, fugitive dust particles varying in size (PM0.1 to PM10) and concentration are generated at various locations of a mine and dispersed by the airflow pattern inside the open-pit mine. The amount and location of dust particles inside the pit is reported at various time intervals. The development of atmospheric inversion during the winter season significantly affected the amount of dust retention inside the open-pit domain. While during the summer, airborne dust particles were transported outside the open-pit domain in a short time. INTRODUCTIONComputational fluid dynamics (CFD) modeling of open-pit mines is challenging. In a recent paper, Bhowmick et. al. [1] presented different constraints and challenges in developing a 3-Dimensional CFD model. In order to obtain realistic results from a simulation, however, appropriate values of various inputs are required. Simulation results of a 3D non-hydrostatic model by Shi et al. [2] showed that both the mechanical and thermal forcing are important mechanisms, controlling the development of the atmosphere inside an open-pit. Alvarez et al [3] used a CFD model for predicting dust dispersions in a medium sized limestone quarry, and presented several iso-surfaces of varying dust concentrations using ANSYS CFX. A similar study by Torno et al [4] introduced a more sophisticated parameterization of the dust sources. Flores et al. [5] used an OpenFOAM CFD simulation package to simulate and predict pollutant dispersion in an idealized and an actual open-pit mine in Chile with an intense insolation condition. The model was simulated for three different conditions: (i) air advection (10 m/s) driven mechanical turbulence; (ii) surface heat flux (240 W/m2) driven thermal buoyancy; and (iii) both air advection and surface heat flux driven conditions. The results of the simulations showed that the buoyant currents contribute to the removal of a large percentage of the dust particles. The buoyant currents modify the purely mechanical-turbulence induced flow pattern and reduce the particle residence time, commonly observed in a purely mechanical turbulent flow (non-buoyant case)"

Jan 1, 2015