Search Documents

By Luis Dr. Sobral

In this study, the removal of zinc, cadmium and copper from an aqueous sulphate solution, derived from the biomass desorption process, was examined by using a reaction system with a three-dimensional high surface area electrode (RNiF). This reaction system was found to be extremely effective in reducing metal concentrations to the ppb range. The effect of cathode potential and flow rate were studied for three different nickel foam porosities. In addition, current- potential curves were recorded, at a rotating vitreous carbon electrode, to deter- mine the possibility of depositing the aforementioned metals separately.

Jan 1, 1998

By Franklin B. Woodbury

As part of its goal of helping maintain an adequate supply of minerals to meet national economic and strategic needs, the U.S. Department of the Interior's Bureau of Mines is evaluating alternative technologies to treat oxidized taconites. Since process control is an essential element in the application of these process technologies, research was performed with a prototype monitoring system. A Californium-252 neutron source and a thermal neutron capture- gamma ray spectra analysis method was utilized to measure the amount of iron and the percent solids in the process slurries. The concentrate and tailing streams in a 410-kilogram-per-hour (900- pound-per-hour) pilot plant employing a selective flocculation and cationic flotation process were monitored. The results obtained showed that the thermal neutron capture-gamma ray reaction can be used to provide process control parameters.

Jan 1, 1979

By Y. P. Chugh, J. Ma

Higher production, better safety standard, and potential for automation are some of the benefits of longwall mining. Today, longwall face advances at a faster rate exposing many diversified rock layers in a short period of time. It is now a serious challenge to cope with ground control problems such as roof falls, face and floor failure, and excessive shield loading as fast as possible to minimize production and monetary losses. In Illinois coal mines, the existence of weak floor strata blow the coal seam may pose additional problems related to floor heaving, shield base punching, and associated roof and face falls. In this study, the effects of weak floor on longwall ground control are analyzed using two-dimensional finite element models. A two-leg 635.6-ton (700-short-ton) yielding capacity shield is included in the models to evaluate the effects of different thickness and material properties of the weak floor. The study indicates that the thickness and material properties of weak floor have significant effects on shield loading, the distribution and intensity of front abutment stress, failure zones in the surrounding strata, roof-to-floor convergence, and floor punching by the shield base.

Jan 1, 2000

By V. Petr, G. G. W. Mustoe, T. G. Rozgonyi

We introduced a numerical method that is applicable for the analysis of static and dynamic deformations of elastic media. In this numerical study, each elastic body is modeled with a system of several hundred rigid particles that are joined together elastically. Maximum shear stress fields due to external loading are obtained and compared with those obtained by photoelastic experiment. Preliminary results are presented for the stresses generated in circular disk shaped bodies idealized with connected systems of particles. Extension of this simulation method to fragmentation of brittle materials is also illustrated in an example.

Jan 1, 2000

By B Muller

The disposal of tailings has become a greater issue due to the changing economics of mineral extraction. The processing of lower grade deposits and tougher environmental regulations has led to the selection process for tailings disposal options becoming more complicated. Current selection methods are too time consuming or too easily manipulated for inexperienced users to use successfully. This paper examines a different method for the initial selection of a tailings disposal system, using a numerical system that was introduced into the mining method selection field. The numerical method system was compiled by reviewing previous work done on the positive and negative aspects of the tailings disposal methods used by the coal and metallurgical mining industries. Assessing the common trends and design aspects of each of the methods was undertaken and a ranking was designed to explain how well suited the tailings disposal methods are to a situation. The tool successfully analysed sites using common data to indicate methods for disposal that would have provided minimal technical problems, low environmental risk and been practically possible. The tool also proved successful in identifying critical areas that would need further investigation. This new system was tested on some disaster case studies from around the world in order to show that it could have successfully selected a different method for tailings disposal that would have prevented the disaster.

Sep 26, 2011

A newly developed damage model was used in the study of rock fragmentation mechanisms involved in crater blasting. Simulations show that at the detonation of the explosive charge, a crushed zone around the explosive charge is formed due to compressive failure and then a strong pressure wave is sent out forming a pressurized ring that propagates at the dilatation wave velocity. As the pressurized ring leaves the crushed zone, two tensile stress components perpendicular to the radial stress are developed in response to the action of the radial pressure, causing tensile failure and forming a damage zone around the charge. The formation of this damage zone is quite independent of the depth of the explosive charge provided there is no interference by the reflection from the free face. When the pressurized ring reaches the free face, it is reflected there and causes rock failure on the ground surface that forms another damage zone. Under the simulation conditions, it is revealed that the depth of charge is a factor controlling the relationship between the two damage zones. At a depth smaller than the optimum, the two damage zones coalesce into one and there is excessive rock breakage; otherwise, the two damage zones tend to isolate from each other with the increase in the depth of charge. Calculations of fragmentation indicate that the average fragment size increases and the uniformity of fragments decreases when the depth of charge increases as observed in practical blasting.

Jan 1, 1998

Impermeable geologic faults in the coal seam can cause intermittent production problems or can cause unexpected amounts of water or gas to issue from degasification boreholes. These faults also can impact methane emissions into the mine workings, especially if they hinder proper and effective degasification of the coal bed. They may also act as barriers for methane flow in the coal seam. Although this might seem beneficial for advancing mine workings, faulting may also cause gas pressure buildups and result in compartmentalization of the gassy regions from which large quantities of water and methane may rush into mine workings. This study uses reservoir simulations to illustrate the effects of impermeable faults, with and without throws, on the production performance of vertical and horizontal degasification boreholes. These boreholes were drilled from the surface to intercept the coal layer and fault and were assumed to produce methane two years prior to mining of the coal seam. Longwall advances and face position related emissions and reservoir properties were characterized using sector definitions along the path of the panel grids. This work numerically evaluates water and gas productions of each borehole and presents the effects and the impacts of impermeable faults on methane emissions due to an advancing coal face. The results can be used to understand the effects of impermeable faults on borehole productions and face emissions due to these coal seam anomalies.

Jan 1, 2008

By A. Mortazavi, H. Nasab

"The storage of petroleum products on surface is associated with difficulties and limitations. Accordingly, a viable alternative is to excavate underground space in rock and provide a safe way for oil storage. Soft rocks such as salt domes provide good conditions from environmental and operational points of view. The potential for high-volume storage and low permeability are encouraging aspects for oil storage in salt domes. The complicated shape of oil storage caverns, the complicated behavior of salt rocks, and boundary conditions associated with large underground openings are important challenges in the design of oil storage caverns. The aim of this paper is to delve into the stability and deformation mechanisms of storage caverns created in salt domes. A parametric study of effect of the cavern size, orientation, and depth is carried out as a function of ground stress regime. The stress field and deformation mechanisms are calculated for each case aiming at providing design guideline.INTRODUCTIONOne of the most popular soft rock that is used for underground oil storage is salt rock (halite). Salt rock has particular characteristics which makes it suitable for oil storage. Salt rock has several advantages for oil storage, these are: high production rate, high degree of availability, short filling period, high level of safety and low permeability [1]. The low permeability is due to the fact that halite deforms in a plastic manner over a wide range of pressures, temperatures and rates of deformations. This property prevents joints or gaps opening up in the formation, thus preventing the escape of liquids and gases. The state of stresses around the cavern in salt domes depends on the depth of the cavern (initial stresses and rock mass temperatures are increasing with depth), the in-situ stress state within the rock mass, and the internal oil pressure in the cavern. Furthermore, the salt creep behavior and the cavern geometry are influencing stress redistributions during loading and unloading phases.Rock mass structure, cavern geometry, material behavior and cavern operation patterns require complex physical models to describe the load-bearing behavior of the geomechanical structure. The determination of stress, strain and deformation in the rock mass due to cavern construction and operation are complex problems which can be evaluated employing numerical tools such as the Finite Element Method (FEM) or Finite Difference Method (FDM). In this study, the finite difference method was used to evaluate the behavior of oil storage cavern in salt rock."

Jan 1, 2015

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 P Lindsay

A numerical model has been developed to simulate the effectiveness and costs of a number of conceptual remedial measures aimed at reducing acidity and metal concentrations in the Ngakawau River and its tributaries near the Stockton Coal Mine, situated on the West Coast of New Zealand. The Stockton Coal Mine has been in operation over the past 100 years. During that time, seepage from waste rock dumps has contributed to elevated acidity and metal concentrations in the Ngakawau River and some of its tributaries. The yearly mean rainfall at the mine site is 6000 mm. The model, developed using GoldSim software, simulates rainfall, run-off and stream discharge for 26 tributaries. It then calculates acidity and metal concentrations in the mine impacted tributaries and the dilution downstream from the mine. The model has been calibrated to measured flow and acidity and metal concentrations in the Ngakawau River and some of its tributaries. Acidity and metal loads are removed from the system in accordance with specified efficiencies and capacities of the respective conceptual remedial measures. The user can compare specific remedial measures, or a combination of conceptual remedial measures. The model predicts the acidity and metal concentrations at a number of monitoring stations for any remedial strategy selected by the user. It also calculates the capital and operational costs for that particular remedial strategy. This paper presents the development of the numerical model and describes how the end-user can use it to compare the efficiency and costs of various remedial strategies by running different æwhat ifÆ scenarios.

Jan 1, 2006

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 S.K. SHARMA, ANKUSH GALAV, G.S.P. SINGH

This paper proposes a hydromechanically coupled numerical modeling approach to quantify the influence of various parameters on the mechanical and hydraulic performance of protective water barrier pillars (PWBP). The results are validated for two case studies having pillar widths of 16 and 42 m at cover depths of 136 and 190 m and water heads of 25 and 141 m. The outcome can be significantly helpful in the rational design of these pillars in underground coal mines.

Mar 1, 2023

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 Janusz Cygankiewicz, Marian Branny, Józef Wactawik

The paper presents mathematical model of concealed fire. We consider a porous medium (goaf, coal fissures) through which air flows. The set of equations is constituted by transient, two dimensional mass, momentum and energy balance equations with sources describing oxygen loss and heat generation caused by coal oxidation. Paper presents several examples illustrating the phenomenon of spontaneous heating in goaf areas behind longwall faces.

Jan 1, 1997

By Ch. Arvanitidis

This paper presents the design solution for an actual case of settlement tanks, for a wastewater treatment plant under construction, actually founded on swelling clay (the nature of which had not been appreciated beforehand) and suffering considerable damage as a result. Two points of interest prompted the use of a numerical model as an alternative approach to the traditional method of analysis: First, the desire to compare the in situ measured deformation and displacement of the tanks with the figures calculated by the numerical model. Second the need to check the validity of the suggested solution as a remedy, i.e. the replacement of a certain thickness of clay by graded material.

Jan 1, 2005

By Patrick Linzer

The majority of existing models describing (methane) gas transport within, and liberation from, coal structures treat these as continua wherein two main gas flow modes (diffusive and Darcy permeable) prevail, and use finite difference techniques to generate solutions to the resulting multivariate parabolic partial differential equation(s). Isotropy and homogeneity in the physical and geometric coal structure properties are often assumed, and the existence of fundamental random fluctuations in these proper- ties is not catered for. The model proposed herein does not suffer these inadequacies, and its efficiency is enhanced by the replacement of the diffusion equation by an accurate approximation. Dynamic mining activity, gas contributions made by the geology surrounding the mining seam, goafing of overlying strata and the effects of induced strain are all taken into account.

Jan 1, 1995

By Dale S. Preece

The computer program. DMC Qistinct Motion &de), which w&s developed for simulating the rock motion associated with blasting, has been used to study the intluence of row delay timing On rock mOtiOn. The munerical simulations correspond with field observations in that vmy short delays (< 5Oms) and vety long delays(> 3oOms) produce a lower percent-cast thao a medium delay (100 to ZOO ms). The DMC predicted relatio&lp between row delay timing and percent-cast is more complex than expected with a dip in the auve where the optimum timing might be expected. More study is required to gain a full undasranding of this phenomenoo.

Jan 1, 1994

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