Search Documents

By K. Becerra, D. G. Almond, T. D. Maru

"Nkomati Nickel Mine in South Africa processes one of the world’s hardest ores in a 375k TPM autogenous grinding circuit. The plant employs an unconventional processing circuit, comprising a primary autogenous mill (AG or FAG) and pebble crushing circuit followed by a secondary hybrid pebble/ball mill and classification cyclone cluster.Advanced Process Control (APC) or Expert Control is now generally accepted as an advantage in most modern high capacity concentrator grinding circuits. This paper reviews the installation of model predictive control at Nkomati and examines the results and benefits of installation in an autogenous mill during plant ramp-up. The paper also demonstrates how predefining process control objectives together with good plant and supplier cooperation can achieve a robust and reliable APC implementation on a very short site program.The results conclusively demonstrate operational and economic benefits during circuit operation with the advanced process control system ON. Production performance and circuit stability were significantly better than during manual control."

Jan 1, 2012

By Donald Leroux, Martin Suichies, Andy Trusiak, Chris Dechert

"In modem concentrators, the operation of flotation machines is often stabilized by simple control loops. Froth level, slurry density and aeration rate are all physical parameters typically controlled through classical PID loops. Few concentrators, however make use of advanced computing technologies for controlling the grade of flotation product streams. Technical challenges such as long process delay and the long and uneven sampling interval of on-stream XRF analyzers often restrain their implementation. Model Predictive Control (MPC) is known to have better performance than PID for controlling processes with long delays (such as flotation circuits). This paper describes the implementation of a multi-input/single-output (MISO) MPC algorithm and an FIR model identification routine using an object-oriented graphical computing environment.The MPC module (controller and identification) was implemented on an industrial control system and applied to many sulfide flotation circuits in the Brunswick Mining concentrator (New Brunswick, Canada). In most cases, the manipulated variable was the collector addition. As long as it is provided with decent process models, the performance of the MPC controller is good on all the flotation circuits on which it was applied. Significant reductions of grade variance are reported and the controller responds well to set-point changes. The user-friendly, graphical nature of the computing environment led to a high degree of acceptance by both the process engineers and the flotation operators.1. INTRODUCTIONOver the past three decades, much instrumentation has been developed for flotation machines. The introduction of instruments such as pulp density gauges and cell level probes was particularly well received by the plants where pneumatic and electronic controllers were available for implementing PID control loops. In a way, these simple loops revolutionized the mineral industry. The operation of flotation circuits could then be stabilized through the use of automatic controllers."

Jan 1, 1999

By F. Vasilikou

Commercial deployment of Carbon Capture Utilization and Storage (CCUS) requires field testing of a scale that can stress the geologic reservoirs. One such reservoir of interest consists of unminable coal seams that exhibit favorable characteristics and depositional environments and lower pressure and temperature than other, deeper, reservoirs. Such conditions can reduce compression costs while utilizing the action of adsorption that offers a more effective carbon dioxide (CO2) bonding than free storage or solution. To ensure the success of such tests, a number of parameters should be accurately determined, such as the seam geometry and stratigraphy, coal porosity and permeability parameters as well as optimum injection conditions. In essence the CO2 injection model should be calibrated for the reservoir characteristics using all available data. History matching calibrations are based on gas and water production from existing wells prior to injection. This paper will present reservoir models that were developed for a pilot--scale 907 tonnes (1000 tons) CO2 injection test that was performed in 2009 through one legacy coalbed methane production well in Russell County, Virginia, USA. The model incorporates a number of individual coal seams about 0.3 m (1ft) in thickness located at depths ranging from 300 to 700 m (1,000 to 2,200 feet). Model calibration was performed through history matching to prior production data, and subsequently the model was utilized to develop different injection scenarios. The developed model was used after injection to calculate CO2 plume distribution patterns that were monitored at the injection test. The paper will present model data and assumptions with a special emphasis on the effect of hydraulic fractures and the skin factor to the coalbed methane (CBM) production model.

Aug 1, 2013

By André Desbiens, Jocelyn Bouchard, René del Villar

"The metallurgical performance of the column flotation process is determined by the concentrate grade and recovery. Although the former can be continuously measured using an on-stream analyzer, the latter must be estimated by a material balance calculation based on seldom verified steady-state assumptions. Consequently, the automatic control and optimization of flotation columns should be performed using secondary variables having a strong influence on the metallurgical performance, such as froth depth, gas hold-up or bubble surface area flux and bias. The aim of this paper is to present a model-based automatic control strategy using froth depth and bias measurements obtained from virtual sensors developed at Laval University. Preliminary pilot plant results on two-phase (water– air) and three-phase system (water – minerals – air) will be presented and eventual industrial implementation of such a strategy as well as the possibility of achieving effective real-time optimization from bias and froth depth measurements will also be discussed.INTRODUCTIONFlotation columns have been widely studied and used in mineral processing for about twenty years. Although a great deal of work has been accomplished to understand the influence of different variables such as froth depth (H), bubble surface area flux (Sb) and bias (Jb) on the metallurgical performance, the measurement of some of these variables is still difficult to be performed – except at the laboratory level – because in most cases, no commercial device yet exists.As a result of this, it is almost impossible to achieve neither efficient control nor optimization of the process at the industrial level, even in open-loop operation. This paper aims at presenting a simple control strategy of froth depth and bias. Description of the sensors used to achieved measurements, preliminary pilot plant results for froth depth control, eventual industrial applications and possibilities for process optimization are also treated."

Jan 1, 2003

By S. Banisi, M. Pourkani, A. R. Laplante, A. Partoazar

"Canadian experience in industrial modeling of grinding, classification, and flotation kinetics was applied at the Sarcheshmeh copper mine, which processes 40,000 t/d of ore at 0.9 % Cu, in an effort to improve capacity and metallurgical performance. Mcivor's functional performance analysis was applied to ball mill cyclone circuits to maximize circuit capacity. Plitt's cyclone model was applied to predict how apex wear would affect classification performance, as part of a predictive maintenance program whereby the optimum wear period is determined for cast iron apexes. The flotation flowsheet was improved using Agar's optimization criteria. The impact of these technologies has been nothing short of spectacular, with an improved capacity of 1 S to 20""/o, and increased recoveries.1-IntroductionThe Sarcheshmeh copper ore body, which may rank as third or fourth largest in the world, contains 1 billion tonnes averaging 0.90 % copper and 0.03 % molybdenum. It is located southeast of Iran and has been processing 40,000 t/d of ore since 1982. The mine produces 100,000 tonnes of copper and 2,200 tonnes of molybdenum concentrate per year. After three stages of crushing, the ore with a Fso of 12. 7 mm is fed to 8 parallel ball mills in closed circuit with cyclones to produce a product 70""/o finer than 7 S microns. The concentrate of the first stage of flotation is reground and the tailings constitute the main component of the plant final tailings. A concentrate with an average grade of 31 % Cu is obtained after cleaning and re-cleaning stages. Total recovery of the plant varies between 83 and 87 % depending on the operating conditions and ore type.To improve the overall performance of the plant, it was first decided to evaluate the performance of grinding, classification and flotation. Measures were then taken to improve performance using the following criteria: (i) The improvements should be of the types that have resulted in significant improvements in other plants (e.g. Shannon et al, 1993); and (ii) The proposed improvements must be achievable with the manual operation presently used at Sarcheshmeh. The outcome of this study could also help in making decisions on the implementation of automatic control systems in the Sarcheshmeh copper mine."

Jan 1, 1998

By Tuula Mäkinen, Lauri Pesonen, Peter Björklund, Jouni Pihlasalo, Mikko Korpi

Outotec Direct Blister Flash Smelting technology has proved to be an efficient and environmentally friendly way to treat copper concentrates from the Zambian Copperbelt at Konkola Copper Mines' (KCM) Nchanga smelter in Zambia. The modeling of slag chemistry and thermodynamics with MTDATA software played a central role in deciding the operation parameters for the smelter. Outotec FSF Online Process Advisor is used to calculate the process parameters and to give valuable information of the process state for the operating personnel. The pre-calculated slag liquidus graphs are also integrated to the Online Advisor.

Jan 1, 2011

By Michael L. Free, Shukla. Abhijeet

The short-circuiting tendency of a copper electrowinning cell (EW) based on the roughness of cathode deposits was investigated with varying parameters such as Cu concentration, guar concentration, current density, and temperature. A bench scale EW cell was set up to produce copper under conditions similar to industrial Cu EW. In addition, a statistical analysis and finite element analysis were performed to analyze the deposited thickness, surface roughness and maximum feature height. Interaction data for the cathode deposits suggest that current density and Cu concentration together affect the roughness of deposit more than the combined effect of guar and temperature. Model predictions of roughness were similar to experimental values. The predicted time to short-circuiting was similar to the measured time. This study provides an industrial relevant baseline against which new operating parameters can be compared.

Jan 1, 2014

By Peter A. Amelunxen, Richard LaDouceur, Courtney A. Young

"Many of the rare earth elements, the lanthanide series on the periodic table, are on the U.S. Department of Defense Strategic and Critical Materials list. In this regard, the U.S. Army Research Laboratory (ARL) has funded several projects at Montana Tech focusing on rare earth element recovery from natural resources. Current froth flotation operations utilize octyl hydroxamate for rare earth mineral extraction, however recovery is not ideal and flotation performance can be improved. Salicylhydroxamic acid (C7H7NO3) was identified as a novel collector due to similarities with n-octyl hydroxamate, notably the amine head group. Salicylhydroxamic acid is primarily used in the medical industry but due to the amine head group is ideally suited as a collector for rare earth mineral flotation.Parametric studies were performed for the flotation of rare earth minerals by modifying a wide range of variables. Mechanical variables of a 6L bench scale flotation cell, such as air inlet flow, rotor speed, and milling time for liberation were modified to elucidate energy input parameters for the desired minerals. Chemical parameters such as collector concentration and pH were varied to understand solution chemistry considerations of the salicylhydroxamic-rare earth mineral system. The study obtained a statistically significant experimental model for the flotation of rare earth minerals. A two stage compartmental model was also utilized for the determination of the first order flotation rate constants for each mineral. Optimized recovery of rare earth minerals was found to occur at near neutral pH ranges using salicylhydroxamic acid with maximum recoveries approaching ideal industrial standards although rare earth mineral grades were still not sufficient.INTRODUCTIONRare earth elements (REEs) are found in minerals that contain several REE which makes costeffective separation of the rare earths difficult. Significant scientific challenges have hindered mineral separation. The molecular similarity of rare earth elements results in minerals containing several REEs, such as bastnaesite [(Ce, La) CO3F]. This prevents efficient separation during concentration steps, for instance flotation, and has a direct impact on successive metallurgical processes. Consequently, concentration during flotation is of rare earth minerals (REM) as opposed to REE."

Jan 1, 2016

Building block cave operations involves significant commitments to capital and time before arriving at completion or full scale production. The focus of the mine planners is habitually to minimize the upfront capital and achieve target development rates and construction dates. However, management of the construction of block cave projects requires understanding of the strategic/tactical decisions that the project planners can take to increase the value of the project and cope with schedule risks. As block cave projects involve development of vertical and lateral access to the underground footprint and ore body, more attention to uncertainties is required as most of those projects need to confront delays leading to a decrease in value. The standard process of development and construction scheduling is often not planned with provisions to realise uncertainties in the construction phase. Uncertainties and risks such as poor ground conditions and equipment breakdowns are the primary reasons for schedule slippage. As a result most of these projects experience delays. In order to mitigate risk to the project schedule, we centre our attention on increasing the value of the project by spending more money upfront and improving the schedule rather than focus on minimizing the upfront capital. Scenario-based simulations along with project appraisal methods hold great promise to address those challenges and enable mine planners to make better decision and enhancing project value. This paper will present a hybrid simulation based approach to model the process of development and construction based on discrete event (DES) and continuous based simulations (CBS). Such simulations will be conducted for scheduling multiple development headings (Apex, undercut, extraction, and haulage levels). Two scenarios are utilized to investigate the impact of risks on schedule and how flexibility in construction can be incorporated to minimize delays: (1) The bench mark scenario is created based on actual data from an existing case study. (2) The flexible scenario is an alternative way to use primary production infrastructure (Ore Handling System) including conveyer drifts, conveyer belt, and a gyratory crusher in the construction process and that requires the system to be in place before starting lateral development in the footprint. As a result, a truck-load system versus a conveying system will be studied to observe the gain in construction speed and to estimate the premium cost for mucking rock from headings and drawbells and transporting it to surface. Such strategic flexibility would be part of the original budget to make the tactical flexibility available during construction. Continuous cash flow models are used to evaluate and compare the scenarios. This study is an example of how new techniques in planning can be used to support decision-making processes in block cave design.

Aug 1, 2013

By H. Chaouki

The aim of this work is to simulate the forming process of a green anode. For this purpose, a nonlinear compressible viscoelastic constitutive law is presented. The concept of natural reference configuration is considered. Within an isothermal thermodynamic framework, a Helmholtz free energy is proposed to take into account the nonlinear compressible deformation process occurring between natural reference configuration and current configuration. A dissipation potential is introduced in order to characterize the irreversible aspect of compaction process. Thus, the constitutive law is formulated through two equations: 1) an expression of Cauchy stress tensor and 2) a differential equation characterizing the evolution of the natural reference configuration. Material parameters are assumed to be a function of the bulk density. This allows capturing the micro-mechanical phenomena within the chosen macroscopic approach. An experimental study based on compaction tests is carried out in order to characterize the mechanical behaviour of the anode paste. A user's material VUMAT subroutine for finite element dynamic explicit analysis has been developed and implemented in the ABAQUS commercial software. To evaluate the model predictive capability, numerical simulations of the compaction forming process of anode paste were performed. Simulation results show that the constitutive law predicts the experimental trends and gives interesting physical responses. This constitutes a first step toward making a benchmark with experimental results on the forming process of anode paste.

Jan 1, 2011

By Luis Bernal, Pertti Lamberg

"SkimAir®, Outotec’s Flash Flotation unit, is placed in the grinding circuit to float the valuable liberated mineral particles from the circulating load of the grinding process while still relatively coarse. Esperanza copper-gold deposit located near Sierra Gorda, Chile, hosts 535 Mtn ore grading 0.57 % copper, 0.27 g/t gold. A series of laboratory and pilot flotation tests were conducted to find out the flotation performance of different ore types both with and without the Flash Flotation. This data was used in modeling and simulating to find the optimal flow sheet for the highest Au and Cu recoveries with required Cu concentrate grade. Simulations with the HSC Chemistry 6.1 software showed that with the SkimAir® and Flash Cleaning it is possible to achieve the required copper grade of 27% Cu with 5% higher recovery for gold and 2% higher recovery for copper compared to the circuit without the Flash Flotation. The pilot tests with simulations were the back up for the installation of four Outotec SkimAir® SK 2400 Flash Cells in the Esperanza circuit for the processing 95.000 tpd of copper-gold ore.INTRODUCTIONEsperanza copper-gold deposit is located near the town of Sierra Gorda, 2nd Administrative Region in Chile. Deposit hosts 535 Mtn ore grading 0.57 % copper, 0.27 g/t gold. The main copper minerals are chalcopyrite and bornite. Gold is carried by native gold and solid solution of copper.Antofagasta Minerals S.A. conducted a series of flotation studies at SGS Minerals Services. The studies included both laboratory and pilot scale work. In the pilot scale study the applicability of the Flash Flotation technology was studied with High-grade Andesite and High-grade Porphyry samples (SGS Mineral Services). The Flash Flotation technology gave significantly higher recovery for gold than the conventional flotation circuit. The recoveries for copper were slightly higher using the conventional flotation."

Jan 1, 2009

By Nader Golchinfar

The V-shaped breakout zones around excavation boundary have been captured by a variety of numerical methods such as FEM and DEM methods and failure criteria such as linear Mohr-Coulomb and non-linear Hoek-Brown failure criteria. Diederichs (2007) used Hoek-Brown parameters in his spalling parameter approach; however, it is convenient to use the Mohr-Coulomb failure criterion in numerical modeling because the material parameters can easily be determined from triaxial tests. It is shown in this study that when peak and residual strength parameters as well as dilation angle are properly defined, the brittle Mohr-Coulomb failure criterion can be used to model depth of failure near underground excavations satisfactorily.

Jun 1, 2012

By B. R. Rakishev

When transferring from the railway transport to the combined automobile-conveyor transport by stripping works at coal mines, it was proposed to form at the working zone a concentration horizon for reloading part of the stripping volumes at the railway transport. The procedure is adduced for substantiating parameters of bloc-panels, depending on the step of transferring of reloading points in interrelation with the intensity of the working zone?s development with observing the balanced movement of parts of the stripping zone relative to the temporary concentration level.

Aug 1, 2013

By R. Panchadhara, P. A. Gordon

"A non-local formulation of the classical continuum mechanics theory called peridynamics is used to study driven dynamic fractures in the near wellbore region. The controlled application of pressure on the inner walls of the borehole can create multiple fractures that initiate from the borehole and propagate radially outwards. The number of fractures and their final length is primarily controlled by the loading conditions and the in-situ stress state. It is the objective of this study to use a novel numerical technique called peridynamics to investigate: 1) the influence of loading conditions on the fracture pattern, and 2) the effect of in-situ stress anisotropy on fracture propagation. In peridynamics, the momentum equation from the classical theory of solid mechanics is replaced with the integro-differential momentum equation of peridynamics. As a result, discontinuities in the body, such as fractures, occur spontaneously. The material is represented by a continuum body with material points that interact non-locally with all other points within a specified distance. Fractures evolve using a damage evolution law, wherein, damage at any point is represented as the ratio of the number of broken bonds in the current configuration to the total number of bonds in the undeformed configuration.Simulations are performed on simplified test cases and the results are compared to relevant experimental and numerical studies found in the literature. It is shown that the loading rate significantly influences the number of fractures initiated from the borehole as well as their extent. Low loading rate produces fewer fractures, whereas, high loading rates produce longer and a greater number of fractures from the borehole. The numerical method successfully predicts the fracture pattern obtained from a wide range of loading conditions. Predictions also indicate that the resultant fractures are shorter when in-situ stress is applied. Numerical results show that the fracture extension is largest in the direction perpendicular to the minimum compressive stress for the range of loading rates used in this study."

Jan 1, 2015

By A. Saigal

A finite element analysis based investigation using DEFORM software was undertaken to model the extrusion of 15% alumina and 15% silicon carbide reinforced aluminum metal matrix composites. The study dealt with predicting the extrusion load as a function of ram displacement and initial billet temperature. The finite element simulation was also used to determine the maximum extrusion load, temperature distribution, maximum temperature, strain rate, velocity and stress distributions in the billet during the hot extrusion process. The constitutive equation used was of the type: Z = E exp(Q RT) = A{sinh(ao-)}" where Z is the Zener-Hollman parameter (temperature compensated strain rate), s is the strain rate, Q is the activation energy, R is the gas constant, T is the absolute temperature, a is the flow stress and A, a and n are constants. Based on the work of Davies et al and McQueen and Konopleva, the values of Q (KJ/mol), a (MPa-1), n and A for aluminum 6061 composites reinforced with 15% alumina and 15% silicon carbide used are 181, 0.024, 4.15 and 1.33x10'2 and 233, 0.052, 2.60 and 1.26xl0'4, respectively. Even though the individual parameters are quite different, it was found that the predicted maximum extrusion load, temperature distribution, maximum temperature, strain rate, velocity and stress distributions are quite similar and differ by no more than 1-2 percent between the two composites indicating that the deformation behavior is primarily controlled by the matrix material and is not a function of the type of reinforcing ceramic particulates.

Jan 1, 2004

By Raj K. Rajamani, Nikhil Dhawan

"Agglomeration has become common pre-treatment step to improve heap leach percolation thereby recovery, the fundamental understanding of the agglomeration process for crushed ores is still lacking. Agglomeration is required for ores which either contain excessive amounts of clay or an excessive quantity of fines generated during mining and crushing. Numerous experiments have been carried out in laboratory scale batch drums using different ores. It has been found out that experimental agglomerate size distributions exhibit self-preserving spectra. The self-preserving/ similarity distribution of the agglomerates is found to be independent of operating conditions and is uniquely determined by the size enlargement mechanism. An empirical model based on the absolute difference between moisture applied and optimum moisture was developed. This model uses self-preserving spectra in the form of a Rosin-Rammler distribution function to predict size distribution at different operating conditions. The empirical model was considerably easier to work with plant scale agglomeration drums. Also, a test known as dip test is developed which improves fundamental understanding of the agglomerate formation and growth mechanisms. This test also provides the progeny size distribution of each agglomerate size. The progeny size distribution was analyzed with a model that partitions the particles into a host and guest category. The ensuing partition coefficient identifies the behaviour of a particle as a host for other particles or a guest layering on other particles. This concept helps in identifying the type of agglomerate that would result for a specific feed size distribution. INTRODUCTIONHeap leaching is one of the hydrometallurgical processes that have served as a green and economically viable process for the treatment of low grade complex ores of copper, uranium, nickel, silver and gold (Dhawan et al., 2013). Heap leaching appears to be a favorable candidate for the least energy consuming process for low grade ores. In heap leach operations nowadays a pretreatment step known as crushed ore agglomeration has been implemented (Dhawan et al., 2013). A bed of agglomerates improves the uniform percolation of solution through the heaps of ore. Therefore, it was possible to treat wastes and mill tailings via agglomeration and heap leaching. Agglomeration pretreatment is required for ores which either contain excessive amounts of clay or an excessive quantity of fines generated during mining and crushing (Bouffard, 2005; Dhawan et al., 2013).While agglomeration has become common practice pretreatment step in heap leach operations, the fundamental understanding of the crushed ore agglomeration process is still lacking. The author reviewed publications up to 2012 on the crushed ore agglomeration process and did not come up with modeling studies on the topic (Dhawan et al., 2013). Some of the possible reasons can be the stochastic nature of the process, experimental difficulties, and lack of instrumentation for fundamental studies."

Jan 1, 2016

By F. Andrieux

"Extruded aluminum profiles are increasingly used for lightweight vehicle construction. Since fracture strains of aluminum profiles are relatively low, damage modeling is crucial for reliable crash simulation. For aluminum profiles, not only the stress state but also the orientation influenced both the deformation and damage behavior. For the material characterization smooth tensile tests were performed in six orientations, both stress-strain curves and r-values were registered. Moreover notched tensile and shear-tension tests were performed in three directions. Finally punch specimens were tested to investigate biaxial loading. Digital image correlation (DIC) analyses were performed to determine local strain values such as critical values at failure. Numerical investigations were conducted with the strain based damage model GISSMO, which is well suited and widely applied for crash simulation. For a strain based failure model the choice of the deformation model is crucial. Deformation models with increasing complexity were investigated to find the best compromise over all experiments. It was found that there is still a great requirement on development of a material model which describes the orientation dependence of the hardening behavior over all stress states in order to be able to use a strain based failure model reliably.INTRODUCTION Weight reduction is an important step for less energy consumption and as a consequence aluminium is increasingly used in automotive parts. An important field of application is the well-established aluminium extrusion which allows the realisation of quite complex shapes of profiles required for an innovative lightweight design with integrated functions. Typical applications are chassis parts, bumpers, crash elements, air bags, etc. Especially the Al–Mg–Si medium strength alloys of the 6000 series can be easily extruded to form complex profiles and prevail for extrusion. Since fracture strains of aluminum profiles are relatively low, damage modeling is crucial for reliable crash simulation. For aluminum profiles not only the influence of the stress state but also the orientation with respect to the extrusion direction have an influence on both the deformation and damage behavior and need to be characterized. The effect of orientation on deformation behavior has been intensively investigated in the last decades and material models haves been proposed (e.g. Barlat & Lian, 1989; Barlat, Lege, & Brem, 1991; Barlat et al., 2003) to describe accurately anisotropic behavior. However it can be noticed that the parameters of anisotropic material models often rely on the yield stress in different directions, which means that the effect of strain hardening is not automatically accounted for. The accuracy of the models at large deformation is not guaranteed and has to be verified. In contrast to deformation models both the characterisation and the modeling of anisotropic damage behavior are not systematically studied. Isotropic failure models based on a critical strain as function of stress triaxiality (ratio of the hydrostatic stress to von Mises equivalent stress) are widely used for crash application (e.g., Bai & Wierzbicki, 2010; Sun, Andrieux & Feucht, 2009; Sun et al., 2013). The mechanisms of failure are manifold and not fully understood; to develop an anisotropic failure model it is mandatory to also have experimental data in different orientations for loading scenarios other than uniaxial tension, e. g. shear, biaxial tension, which implies an extensive experiment plan. This makes the development of failure models more challenging than deformation models for which it is possible to develop a reasonable anisotropic model using the framework of yield function and calibrate the material parameters based on results of tensile tests in different directions. Moreover the development of deformation and failure models are intricate as the variables governing the damage evolution are predicted by the def"

Jan 1, 2018

By K. Baris

Methane (CH4) present in coal seams, generally known as Coalbed Methane (CBM), Coal Mine Methane (CMM) or Coal Seam Gas (CSG), is a major energy source, as it is the primary constituent of natural gas. At the same time, released methane is a potent greenhouse gas (GHG) and a major safety hazard in underground mining operations. Therefore, a considerable effort has been undertaken to develop coal bed methane extraction and utilization technologies that can achieve both environmental and safety benefits. Countries with significant coal reserves and gassy seams, i.e., USA, Australia, Canada, China, India, Poland and Germany have developed mature methodologies and practices to extract, before and during mining, and utilize significant amounts of methane from coal seams. Nevertheless, on a global scale, methane recovery remains untested and many technical and economic challenges require urgent attention. The research presented in this paper, sponsored by the United States Environmental Protection Agency (EPA), is focused on the Zonguldak Basin, the only hard coal region in Turkey. The basin is known to contain very gassy coal seams and consists of three major regions: Armutcuk, Zonguldak and Amasra. Among the three, the Amasra coalfield is considered as the most favorable for CBM recovery, as it comprises of deeper coal seams and where a private company, Hema Enerji A.S. (Hema) is already developing coalbed methane production facilities. This study aims at developing a reservoir model to predict the potential CBM recovery from the CBM-2 well, the only producing CBM well drilled by HEMA Enerji A.S in the Amasra coalfield. History matching was utilized to model the 15-day production data and uncertainty analysis was performed to estimate the optimum values for the reservoir parameters. The produced model yielded an average production rate of 143 m3/day for 10 years. The total amount of CH4 to be produced during this period was 0.53 million m3. Furthermore, the total amount of GHG to be avoided was estimated at 6585 tCO2e.

Aug 1, 2013

By C. J. Li, J. L. Feng

In the present study we focuses on developing the poroelasticity model that deals with adsorptioninduced permeability change mostly occurred in porous media constituted several phases such as coal. The underlying pressure induced by adsorption is obtained by Virial expansion form, which is substantially founded on statistical physics. The corresponding thermodynamic pressure applying on the pore surface is then used to detail the pore pressure and further to evaluate the macroscopic volumetric strain and the porosity. In the process of physical adsorption, the adsorbed phase is in thermodynamic equilibrium with the gas phase regardless of pore size, and the corresponding chemical potentials are all equal. The relationship between surface strain and volumetric strain is derived under infinitesimal deformation conditions to gain an explicit formulation of pore surface area with respect to volumetric strain and porosity. After the formulation of surface stress modified by the Langmuir isothermal adsorption is derived, the differential micromechanical constitutive equations of porous media are obtained. By supposing that the ratio of current permeability to initial is equal to the cubic ratio of current porosity to initial porosity, the differential equations for permeability change is consequently derived. The numerical simulations by the present model for permeability change of coal induced by adsorption of CO2 and CH4, respectively, are presented under confined stress conditions. The predictions from the model are good in agreement with the experimental data of sorption-induced permeability change of coals in the literature.

Jan 1, 2015

By M Reza Afshar

In this work a mathematical model was developed to study the effect of electromagnetic force on the magnesium melt containing nonmetallic inclusion in an electrical isolated channel. The electromagnetic force exerted on the inclusion induced by DC current field was calculated. A three dimensional turbulent flow equations were solved numerically to compute the velocity field within the channel. On the basis of computed velocity and magnetic fields, the trajectories of inclusions were calculated using the equations of motion for inclusions. Parametric studies were carried out to evaluate the effect of various parameters such as current density, channel diameter and particle density on the removal efficiency of inclusions from magnesium melts.

Jan 1, 2004