Search Documents

By André Vervoort, Jef Caers

Compound lognormal distributions serve as new generalized distribution laws for analyzing and estimating ore value distributions or related properties such as the average ore grade. These models are an alternative to the well known lognormal distribution. They are the translation of general geological constraints on the genesis of the deposit cumulated in suitable statistical properties. In this paper a maximum likelihood framework is established for estimating the parameters of such laws or parameter functions thereof. It is shown that under certain conditions maximum likelihood estimation is accurate and efficient in estimating the mean value even when relatively few samples are available. The method is based on a numerical maximization of the likelihood function of these models given the data. Confidence limits on these estimates are set up using resampling techniques such as the bootstrap. By this a complete methodology is presented for inference on ore grade distributions such as is shown in an example of gold data from the Merriespruit mine (South Africa) and the diamond size distribution of the Bougban deposit (Guinea).

Jan 1, 1996

By M. Nejati, A. Paluszny, R. W. Zimmerman

A numerical method is proposed to propagate multiple discrete fractures leading to primary fragmentation in a mine that is being worked using block caving. Deformation is computed using the finite element method, fractures are represented explicitly, and the mine domain is discretized by an unstructured mesh. Bedding planes are represented by systematically varying the elastic modulus of the rock and by defining horizontal weakness planes. Fractures and matrix are represented using parametric surfaces, and tips are defined by their boundary curves. Tip advance is controlled by a failure criterion, and a criterion for propagation direction and magnitude, based on the evaluation of the modal stress intensity factors. A novel domain integral approach is applied to accurately compute stress intensity factors (K) ahead of fracture tips in three dimensions. The method does not require a structured volumetric mesh structure around the crack tip, as integration is performed over a series of virtual surface domains along the crack front. The method is efficient, as it makes direct use of automatically generated, arbitrary tetrahedral meshes, and approximates stress intensity factors (KI, KII, KIII) along each crack front using Interaction-integrals. As opposed to the JIntegral, the method does not decompose K a-posteriori, but instead uses an auxiliary field to directly compute modal K. When using this method, numerical approximations of K do not exhibit dependence on the mesh layout, and require meshes that can generally be ten times coarser than are required by displacement- and stress-based methods. Volumetric meshing requires only, on average, 17% of each computation step. Thus, cracks do not follow any pre-existing mesh structure, and the method is well suited for high-density fracture datasets. The method is demonstrated for primary fragmentation of a mine area covering 110 initial draw points, immediately beneath a 2m high undercut. Displacement is constrained at all boundaries except the surface of the undercut. The growth of the fractures is investigated at the onset of the undercut. Ninety initial disc-shaped fractures are taken into account, each having an initial radius of 10 m. Fractures grow around the undercut, intersect the bedding plane boundaries, and the domain is fragmented as a result.

Jan 1, 2015

By Joshua R. Omer

A numerical model is developed for predicting the load transfer and settlement behaviour of vertically loaded, bored cast in-situ piles. Based on analysis of full-scale pile tests, closed-form equations are formulated to uncouple the shaft and base transfer resistance components of the pile/soil system. Appropriate parameters are introduced to take into account: (i) the shaft resistance distribution in unlined and partially lined piles, (ii) the effects of ?soft toe? (voided toe pile) behaviour on base resistance development. For convenience, an interactive computer program has been developed in Matlab? with a user friendly interface system. The program can be applied to predict the complete load-settlement curve and axial force distribution in bored cast in-situ piles if reliable soil test data are available. To demonstrate its validity, the proposed method is applied to 4 test piles, which are also analysed using two previously published methods: (a) elastic solution method and (b) layer discretisation method. Methods (a) and (b) above are based on load transfer (t-z) and finite element models respectively. It is shown that the proposed method yields axial force distributions and load-settlement curves that are in good agreement with the observed pile test data.

Jan 1, 2010

By S. K. Sharma, ANKUSH GALAV, G. S. P. Singh

Providing a suitable protective water barrier pillar (PWBP) is common to reduce inundation hazards in underground coal mines. The imperative factors influencing its performance include the water head acting on the pillar, cover depth, pillar width, strength properties, and permeability characteristics. The mechanical failure of such pillar is a stress-controlled phenomenon, whereas the hydraulic failure is a strain-based phenomenon. A finite-difference numerical modeling approach was developed to study the hydro-mechanical coupled behavior of protective water barrier pillars. The mechanical stability was evaluated in terms of the percentage of failed (ZoF) and intact zones. The influence of the strain-controlled weakening on the permeability of the flow medium was studied through the coupling of the mechanical and hydraulic effects. The coupled steady-state model was used to estimate the outflow rate and its hydraulic stability. The adequacy of the protective pillar was also investigated by assessing mechanical stability and capability to resist hydraulic pressure against the maximum expected water head. A seepage rate-based classification system has also been proposed to evaluate the seepage potential and assess the hydraulic stability. The model has been validated for two case studies at the cover depth of 136–189.5 m and the existing pillar width of 16–42 m against the water head of 25–141 m.

Aug 18, 2022

By Gang Mei, Nengxiong Xu, Yanlin Zhang

After open-pit mining finishes, an iron mine with steeply inclined ore-bodies will generally transform to underground mining. Non-pillar sublevel caving method, with horizontal slice mining from up to down, is commonly used for mining deep ore-bodies under the open-pit bottom. As horizontal slice mining proceeds, a gradually enlarging goaf will be formed under the open-pit bottom, the rocks surrounding the handing and heading sides of the goaf will collapse continuously, and then these results in rock-mass failure with large deformation and displacement in the open-pit slope. Based on FLAC3D, a modeling approach to predict pit slope deformation and failure induced by transitioning open-pit to underground mining is proposed in this paper. This method mainly consists of three procedures of modeling, i.e., mining modeling, landslide modeling, and filling modeling, of which landslide modeling is the core of the proposed approach. In the procedures of modeling with FLAC3D, iteration non-convergence is taken as a local failure criterion of pit slope. When iteration calculation cannot converge, the sliding modeling will start by deleting some elements with large displacements, and finish until convergence. Before modeling landslide, several displacement contours of rock-masses in the pit slope are searched out, and each contour is an approximate possible sliding surface of the open-pit slope. From these contours, the most proper one is selected as the slope sliding surface by using the iteration calculation convergence and minimum element deletion as the judgment condition. Landslide modeling is then realized by removing the elements outside the slope sliding surface. Finally, to demonstrate the effectiveness, the proposed approach is used to simulate the process of open-pit to underground mining at Yanqianshan Iron Mine, and to predict the deformation and failure characteristics of open-pit slope during underground mining. The case study results show that: the prediction of the slope stability determined by using our method is consistent with that calculated by using the limit equilibrium methods.

Jan 1, 2015

By M. D. Voegele, C. Fairhurst

INTRODUCTION A computer study of excavation support loads in tunnels and other excavations in discontinuous rock masses was undertaken with the distinct element method introduced by Cundall (1971). In the distinct element method the rock mass is modeled as a composition of individual blocks of rock defined by the joint planes. Normal and shear deformation laws are ascribed to the joint surfaces; the blocks of rock themselves are considered to be rigid. The particular version of the program utilized in this study ran a mini computer and used inter, active techniques for the input of geometrical data such as joint orientation and spacing and subsequent display of rock mass deformation and other data. Two major aspects of the behavior of excavations in jointed rock masses are examined in this paper: i) the behavior of unsupported excavations in rock masses with varying joint configuration, joint properties and in situ stress fields; and ii) the modification of the behavior of similar excavations when support is provided. In the latter case the excavation behavior was examined through the use of ground reaction curves. THE DISTINCT ELEMENT METHOD The Distinct Element method introduced by Cundall (1971) is a computer based analysis that simulates the behavior of a system of discrete, semi-rigid rock blocks. Block interactions are governed by realistic friction and stiffness laws; each block may undergo unlimited displacement and rotation while progressive failure is modeled.

Jan 1, 1982

By Saba Gharedash

In this study the results of a series of numerical simulations which are conducted for one of the tunnels of the Mashhad metro line 2 project are presented to evaluate the role of tunnel face support pressure and soil strength parameters on deformation of tunnel lining and face deformation. To achieve this goal, comparative numerical simulations of a mechanized tunnel advance in soft non- homogeneous soil on top of the ground water table and beneath the building are conducted and sensitivity analysis is performed. Development of a step by step tunnel construction process is modeled using a three dimensional finite difference model (FLAC3D) taking into account all the relevant components of shield tunneling. In order to understand the mechanism of the effects of a variety of support pressures on face deformation and lining stresses, numerical simulations are conducted. It is found that the use of different face supports causes a considerable influence on the forces of tunnel lining. In addition, it is revealed that building weight help to change the stress paths and the stress distributions in the tunnel lining.

Aug 1, 2013

By Edson Charikinya, Pierre-Henri Koch

"The extraction of mineral values from ore requires liberation followed by separation steps. Liberation is achieved by size reduction operations which are energy inefficient processes typically dominating the energy consumption in a mineral concentrator. As the grade of ore reserves declines, future viability of mineral operations will be determined by energy costs of comminution. The application of high power microwave energy to secondary crusher products has been suggested as a possible commercially viable thermal treatment method for reducing comminution energy and improving value mineral liberation. Recent studies have shown that microwave pre-treatment of coarse sphalerite ore particles (> 5mm) at specific microwave heating energies (1-3 kWh/t), induces microfractures and creates new crack surfaces. This suggests that subsequent crushing of these microwave treated particles could yield enhanced liberation. However, limited studies have been carried out investigating the mode of breakage and the extent of enhanced liberation in that case. The objective of this study is to develop numerical methods for quantifying the extent of enhanced liberation and mode of breakage in crushed microwave treated and untreated particles. Sphalerite ore particles representing small (-5+4.75) mm, medium (-16+9.5) mm, and large (-25+19) mm HPGR and cone crushed particles were microwave treated at specific energies between 1-3 kWh/t. Cracks in the ore particles before and after microwave treatment were analysed with QEMSCAN and numerical models of the measured particles (before and after microwave treatment) were developed in MATLAB. The propagation of random and non-random cracks was investigated by simulating comminution of the modelled treated and untreated particles. Results of this study demonstrate that the breakage mechanism and the liberation of valuable minerals from gangue in microwave treated and untreated particles is significantly different. This study contributes to the development of numerical tools to quantify crack propagation in heterogeneous ore particles."

Jan 1, 2018

By M. Ghamgosar, Erarslan N.

Rock cutting includes many unknown parameters which result in different consequences of the different cutting conditions. However, based on the current suggested models for rock cutting, there are some critical parameters that control the fracture initiation and propagation in the rock domain. In general, there are two major cutting tools are used in rock cutting: drag bits and indenters. Different mechanical attacks of cutters result in different fracture models. For example drag bits break the rock parallel to the rock surface, which is called ‘undercutting’ while indenters cut rocks by pressing normal to the rock surface. Both tools cause rock failure by inducing tensile cracks in the rock. However, drag bits induce more tensile cracks compared with the indenters. The full benefits of improvement in rock cutting efficiency will arise from a clear understanding of rock fracture mechanics principles. Therefore, fracture mechanics principles will be applied to the rock cutting as the rock fragmentation process is based on the rock fracturing. The Oscillating Disc Cutting (ODC) has more advantages than drag bits and disc cutters. The forces transmitted back through the cutting head are much lower with ODC. Although ODC attacks rock in much the same way as drag bits, a combination of the oscillating action of the cutter and the undercutting mode produces the lower cutting forces and causes rock fatigue. Fatigue phenomenon in rock has advantages for the rock breaking by applying less effort and energy due to the cyclic loading. The aim of this study is using the Extended Finite Elements Analyses (XEFM) numerical modelling for simulating the ODC cutting with different cyclic loading conditions. An accurate prediction of numerical investigation requires detail laboratory experiments. A series of laboratory tests were performed to obtain the fracture toughness values of our Monsonite rock samples under various cyclic amplitude and frequency values. The numerical results showed different types fatigue damage mechanisms due to the various cyclic loading conditions. Thus, modelling and understanding the crack propagation behaviour depending on different amplitude and frequencies are believed very important for the dynamic rock cutting researches and rock cutting machine manufacturers.

Jan 1, 2015

By Juhyung Lee

Recently, large diameter drilled shafts are often used for foundations of large-scale marine bridges with high capacity and required to resist lateral loads. Steel casings are essential in supporting the hole-diving concrete installation of the drill shafts constructed in the sea. The casings, however, are usually not considered as the structural elements and their load bearing capacity is often overlooked in the design of drilled shafts. As the diameter of drilled shafts becomes larger, there are higher demands for more efficient and economical design of large-scaled drill shafts in the sea. Steel-concrete composite columns are popular for the superstructure of bridges, and the outer steel attached to the shaft increases the shaft resistance through confining the concrete. In this study, the lateral resistance of steel-concrete composite drilled shafts was evaluated quantitatively based on the numerical analysis. When steel casings are used as structural elements such as the composite columns. Ultimate lateral bearing capacity of composite drilled shafts with various diameters was numerically calculated through the 3D finite element analysis. For that, elasto-plastic model with the perfect plasticity is involved to capture the ultimate load. A commercial FEA program, MIDAS-GTS, is used in this study. Real field conditions of the west coast, Korea were examined to set up the ground conditions and pile lengths required for this parametric studies. Detailed characteristics of the stress and displacement distributions are evaluated for the better understanding on the mechanisms of the composite shaft behavior.

Jan 1, 2008

By Emile Scheepers, Chris Aldrich, Jacques J. Eksteen

The pneumatic injection of CaC2 powder into molten iron has become a very reliable and effective way to carry out desulphurisation. The key variables that influence the kinetics of this process are the condition and composition of the top slag, the oxygen activity in the hot metal and the injection conditions. Notable injection parameters include carrier gas flowrate, reagent flowrate, lance depth and initial sulphur concentration. A one-dimensional, steady-state model for momentum, heat- and mass transfer in rising gasliquid- powder plumes was developed for conditions relevant to the CaC2 refining process at Saldanha Steel?, South Africa. Combined with a model that predicts the contribution of the topslag to the process, the overall rate of desulphurisation as a function of time was determined, thereby quantitatively demonstrating the influence of the above-mentioned injection parameters on the kinetics of the desulphurisation process.

Jan 1, 2003

By Whiten WJ

A method that allows data to be mass balanced and screen models to be fitted in a single step is described. The inclusion of several variables to describe the feed data in the fitting procedure allows a mass balance to be conducted at the same time as the fitting of the screen model. Analysis of Groote Eylandt Mining Company screening data has demonstrated that this method is ideally suited to applications where complete stream data are not available. Also, certain random fluctuations in screen model parameters arising from the fitting of consecutive mass balances and screening models are significantly reduced. A cubic spline function with fixed knot positions corresponding to particle sizes and variable (fitted) cumulative percent passing values is used to describe the feed size distribution. When the spline function is interpolated at the experimental feed particle sizes a fitted feed size distribution is obtained. This feed size distribution is then used as input to the screen model calculations from which product size distributions are generated. Since both feed and product size distribution data are used to minimise the residual errors between the observed and fitted stream size distributions, a mass balance is Conducted simultaneously with the screen model parameter fitting.

Jan 1, 1984

By Ana E. Bohé, Georgina De Micco, Juan P. Gaviría, FEDERICO J. POMIRO, Gastón G. Fouga

Chlorination and carbochlorination processes have industrial applications in metallurgical activity for extraction and purification of uncommon metals such as rare earth elements (REE). This study reviews the chlorination and carbochlorination reactions of light REE (LREE) oxides with chlorine considering the thermodynamic, stoichiometric, and kinetic aspects of the reactions. Lanthanides between lanthanum to gadolinium were considered as LREE, in agreement with previous literature, and the results were detailed in this review. LREE were analyzed taking into account their oxidation states and were divided into three groups accordingly. The results presented show that, although the LREE are chemically similar, their behaviors during the reactions with Cl2 and Cl2–C, and the different LREE-O–Cl compounds formed, have particularities. It was observed that gadolinium, in both thermodynamic calculations and experimental results, presents some differences, and it can be considered intermediate between light and heavy REE; i.e., reactions between LREE2O3 and Cl2(g) begin at approximately 250 °C, except for Gd2O3 that begins at 327 °C. The present work will contribute to a better understanding of the chemistry of LREE oxides and Cl2 under different conditions and the products obtained in the LREE–O––Cl2(g) and LREE–O–C–Cl2(g) systems.

Oct 15, 2021

By ANA E. BOHE, MICCO GEORGINA DE, GASTON G. FOUGA, JUAN P. GAVIRIA, FEDERICO J. POMIRO

To the best of our knowledge, this paper comprehensively summarizes a thorough study to determine the thermodynamic viability of the reactions between light rare earth oxides and chlorine with and without carbon. The kinetics and mechanisms of light rare earth oxide chlorination and carbochlorination reactions are analyzed and discussed according to results in the literature. Lanthanides between lanthanum (La) and gadolinium (Gd) are considered as light rare earth elements (LREEs) in agreement with previous literature and the results detailed in this review. Both laboratory- and industrial-scale light rare earth chlorinations and carbochlorinations, starting from oxides and ores, are reviewed.

By P. C. Clark

This paper is an overview of the requirements and concerns for minerals and pigments used in papermaking. The first part of this paper deals with mineral or pigment filling in the actual papermaking process. The second part deals with the pigment or mineral properties that are important in today's coating processes.

Jan 1, 1993

By P Viljoen, A Holland, C Parken

In the raisebore industry, there are two known methods of removing large diameter reamer heads once the reaming of an underground shaft is complete: 1. Lower the reamer head to the base of the shaft and use a thermal lance to cut the drill string and drag the reamer away from the hole. This method exposes the person doing the job to falling rocks from the unsupported mineshaft, a hazard which is difficult to control and has caused fatalities. 2. Install breast plates in the backs above the reamer head at the top of the hole, and use air hoists to lift the reamer out of the hole. Steel beams are then placed across the hole which the reamer is lowered down onto, dragged forward away from the hole, disassembled, and removed. This method is time consuming, expensive, and complex in the controls required to ensure safe completion of the job At Carrapateena, the air hoist method had been successfully used, however, with the cost involved and the complexity of this process, development of an improved and fit for purpose reamer removal method was pursued. OZ Minerals collaborated with Pybar to develop an industry first application of Enerpac SL100 lift and shift technology; a track mounted system with hydraulic lifting units that can both lift the reamer out of the hole and tram it back along the tracks away from the hole. A key component of turning the idea into a reality was by establishing a truly collaborative relationship between client and contractor. When the SL100 mobile gantry arrived, it was set up in Adelaide and trialled with the whole raise bore crew attending the two-day workshop and risk assessment; a rare opportunity for drillers. There have now been five successful reamer removals with the SL100 mobile gantry at Carrapateena. It has proven to be an efficient process with the entire process from setup to demobilisation taking less than 48 hours. This new way has completely removed hot works and uncontrolled hazards from working at the bottom of a shaft. It has also saved the reamer stabiliser from being cut and thus wasted. Exposure to the open hole has been minimised significantly as assembly of the SL100 is completed well forward of the open hole. This new process is exponentially safer, quicker and cheaper.

Mar 15, 2021

By Morgan M. Sears, Gamal Rashed, Michael M. Murphy, Brent Slaker

In this study, a parametric study was conducted using FLAC3D numerical models to examine the impact of oblique loading, generated from seam dip, on the strength and the failure propagation pattern of a stone pillar using two simplified geometry types. In type 1, the sidewalls of the pillars were assumed to be perpendicular to the roof and the floor, while in type 2, the sidewalls of pillars were assumed to be vertical. The complex pillar geometry in dipping mines was frequently modeled using these two geometries. To capture a complete picture of the effect of seam dip on pillar stability, the modeled width-to-height (W/H) ratio of the pillars, in situ stress field, and pillars roof/floor interfaces were systematically varied to account for the potential distribution of values for these parameters across the underground stone mines in USA. Results from the numerical modeling indicate that dipping pillars have reduced strength compared with horizontal pillars. Also, an asymmetric failure propagation pattern could be obtained depending on an interaction between the W/H ratio, seam dip, in situ stresses, and pillar geometry.

Feb 5, 2021

By Morgan M. Sears, Gamal Rashed, Michael M. Murphy, Brent Slaker

A parametric study was conducted using FLAC3D numerical models to examine the impact of oblique loading on the strength and yield pattern of a stone pillar using two simplified geometry types. In type-1 the side walls of the pillars were assumed to be perpendicular to the roof and the floor, while in type-2 the side walls were assumed to be vertical. These two geometries are frequently used to simplify/model the creation of pillars in dipping mines. Results from the numerical modeling indicate that dipping pillars have reduced strength compared to horizontal pillars, and they are at an elevated risk of failure. Also, an asymmetric failure propagation pattern could be obtained depending on the interaction between the pillar width/height (W/H) ratio, seam dip, in situ stresses and pillar geometry.

By Neeraj Saxena, Madhu Huggahalli, Jerry Bednarski, Ken Grieshaber, David Stoffel

"The use of oxy-fuel burners in secondary aluminum melting applications offers several advantages over air-fuel burners including reduced fuel consumption, faster charge to tap times and lower NOx emissions. Their successful, safe and economical use in a furnace depends on several factors and considerations such as burner and flue placement, metal circulation, charge practices and the type of refractory used in the furnace.These factors required in the successful conversion of furnaces from air-fuel to oxy-fuel based burners are discussed in this paper. Critical parameters are identified and examined using computational fluid dynamics (CFD) simulations. Parameters are estimated via laboratory testing and validated through trials performed at commercial installations. Guidelines and a simplified approach to estimate a priori the economic impact of converting from air-fuel to oxy-fuel are presentedThe final products of this research are simplified and validated tools for modeling aluminum furnaces. These include improved heat and energy balance models and parameterized CFD solutions which provide rapid customization to individual furnace configurations. These tools provide field engineers with immediate and accurate predictions of performance, allowing for repeated precise scenario analyses."

Jan 1, 2001

A parametric study was conducted to establish the optimum operating conditions for the column flotation of fine coal. The parameters examined in this investigation included collector and frother dosages, feed solids content, and the volumetric flow rates of feed slurry, air and wash water. Despite the wide range of conditions studied the experimental data were found to fall along single grade-recovery curves, just below the maximum separation predicted by the release analysis technique. The deviation of the parametric test data from the release curve was attributed to the nonselective entrainment of mineral matter associated with poor froth washing. The operating point on the characteristic recovery-grade curve was found to be primarily determined by feed rate and aeration rate.

Jan 1, 1992