Search Documents

By Zoran M. Stevic, Mirjana M. Rajcic-Vujasinovic, Zvonimir D. Stankovic

"Huge electrochemical investigations have been performed on Technical faculty in Bor (Serbia) which main goal was the clarification of the copper sulfide minerals behaviour in hydrometallurgic and flotation processes. Modeling of these processes was a part of the effort for better understanding of the reaction mechanisms.The consideration of electrical analog circuit and mathematical model of processes at interfaces copper sulfide/electrolyte (CuS/(CuSO4+H2SO4), Cu2S/H2SO4 and CuFeS2/H2SO4), based on studies of the overvoltage – time transients, will be presented. The transients were obtained using single pulse galvanostatic method and electrochemical cell with standard three electrode assembly. In that goal it was developed an original system for computer measurements that enables the application of several methods for electrochemical investigations (galvanostatic, potentiostatic, cyclic voltammetry, electrochemical impedance spectrometry etc). The system was tested and adjusted using commercial electrochemical equipment made by AMEL, Milano.In the first stage of consideration the different equivalent electrical circuits have been tested by computer programme OrCAD 9.2 in order to identify the circuit with responses equivalent in shape to experimentally obtained ones. In that way it was established a configuration of an analog electrical circuit for investigated electrochemical system. That phase enabled only very rough determination of values of circuit components and their relative ratio.In the second stage the method for determination of all elements of electrical analog circuit is described. It is based on mathematical consideration of a given circuit on one hand, and valves of the overvoltages at selected times, obtained from charging and discharging parts of galvanostatic transients, on the other."

Jan 1, 2003

By J. L. Feng, F. Zhang, S. K. Hu

"Cracking is one of the most serious damages in tunnel linings, which is often observed by in-situ inspection to tunnel engineering. It is difficult to better characterise for the deformations from continua to discontinua by the traditional continuum mechanics due to emerging of crack initiation, crack propagation, and coalescence, at several scales. In this paper, the process of failure that was from micro damage to macro crack in concrete tunnel linings under experimental load conditions, was quantitatively simulated by the combined finite-discrete element method (FDEM). Comparison of the simulating results with that obtained by the full-scale test of concrete lining under different loading conditions in the literature was respectively good agreed and shows that such FDEM technique is appropriate to model the continuousdiscontinuous problems. The present study indicates that the FDEM technique could be benefit to quantitatively evaluate the stability and damage level of tunnel lining system. INTRODUCTIONMost damage to tunnel linings is caused during the construction phase, which can emerge in a form of single macro cracks or even large-scale spalling. The cause of such damage is mostly the occurrence of, either impermissibly large contact stresses resulting from high thrust loads combined with production and installation tolerances in the segments resulting in a geometrical lack of fit (Maidl et al., 2012; Transportation Research Board, Basic Research and Emerging Technologies Related to Concrete Committee, 2006), or unallowable tensile/shear stresses appearing in the tunnel linings under complicated loading conditions (Ansell, 2010). Such a macro failure phenomenon may also be regarded to as a result of growth, coalescence, initiation of microcracks in tunnel linings due to the complicated loading action, and even potentially lead to an instability or collapse of the local or entire tunnel system. Much work has been devoted to investigate for the characters of micro and macro cracks, the influence factors of damage patterns, the settlement of a tunnel, the width, length and density of lining cracks, especially from classifiedly and qualitatively (Wang, 2010; Huang et al., 2013). However, the safety evaluation of tunnel is definitely required to not only inspect in situ, but further also quantify computing and analyzing for tunnel lining systems, so as to provide a global viewpoint of process of the tunnel deformation, which aims to help for subsequent decisions, e.g. controlling traffic flow or immediately repairing task. Additionally, the design of lining thickness and its structural analysis that is mostly determined by previous experience and similar construction examples is also needed to be improved."

Jan 1, 2015

Sedimentary deposits such as sandstone are formed by the aggregation of cemented sand grains. Phosphates and laterite present similar aggregation of mineral grains. Grains and cement have different mineralogy and chemical characteristics. During size reduction through crushing or coarse grinding, the rocks are crumbled. Mineral grains and cement are then liberated, and the produced size distribution reflects their ?natural size distribution?. Mathematical modelling of these concomitant effects of crumbling and mineral liberation necessitates a specific material description and the adaptation of the commonly used theories of breakage kinetics and population balance. The coarse particles are mainly agglomerates made of mineral grains and cement particles whereas the fine particles are grains and cement with a natural grain size distribution. The population balance between the different particle size and particle types during size reduction takes into account the breakage rate and the progeny size distribution. The breakage rate can be described by a simple probability function or by a first order kinetics whereas the progeny size distribution is described by a breakage matrix. The model proposed here considers that the agglomerates will produce other agglomerated particles, liberated grains and cement particles; the two latter having their natural size distribution. The liberated grains and cement particles are subject to size reduction and follow the most commonly used theory of breakage. The present paper describes the theoretical approach of the concomitant crumbling and liberation effects. The application to a mathematical model for the cone crusher is presented. The advantage of such approach is then demonstrated through a simulation study for the optimisation of a crushing plant. Keywords: mineral liberation, modelling, crushing, sedimentary ore

Sep 1, 2012

By B Yeganeh, A Sam, A Kohsari

The base of separation at Gol-e-Gohar iron ore plant (Sirjan in the south of Iran) is magnetic separation after milling. Magnetic separation in the plant is divided into two parts: dry magnetic separation and wet magnetic separation. Sixty-five per cent of the plant production belongs to the dry magnetic separator section. Therefore, investigation of this part of the plant is necessary. Prediction of plant production is one of the most important subjects in mineral processing plants. That means, before feed enters the plant, specification of production has been calculated with laboratory methods. The purpose of this paper is the prediction concentrate quality has achieved from the dry magnetic separator. The content of iron in the final dry concentrate will be predicted, with laboratory, pilot and modelling methods, before feed enters the plant. Different variables are used in this approach, which consist of: speed of the first dry magnetic separator (rougher), speed of the third dry magnetic separator (cleaner), iron content in plant feed, size of 80 per cent passing of feed, size of 20 per cent passing of feed and iron content in pilot dry magnetic separation concentrate that performs in a standard manner. The artificial neural network is used for modelling and one of its branches is the radial basis network (RBN). Iron content in magnetic separating concentrate in plant will be introduced as the target and will be calculated with the RBN method. In the RBN method, a set of data enters to neural network and they will be trained and simulated. Finally, the target will be predicted. Ultimately, the original target (actual iron content) and the predicted target (predicted iron content) compare and error of them will be calculated. The obtained result is acceptable and satisfactory.

Jan 1, 2007

By Viraj Parekh, Kinnor Chattopadhyay, Sina Mostaghel, Xiao Zhang

Electric arc furnace steelmaking is a dominant approach in North America and controlling EAF emissions is a key issue today because of stringent permit limits imposed by various regulatory bodies. In the present study, an attempt has been made to mathematically model EAF emissions, and use the models as a predictive tool for fine tuning the process. CFD modelling has been used to model various fume release scenarios, and optimize the exhaust conditions for maximum capture. Thermodynamic calculations have been used to analyze various chemical reactions which lead to the emissions of hazardous air pollutants with special emphasis on Zinc, Lead and Manganese.

Jan 1, 2015

By D. K. Chibwe, P. Taskinen

"Typical operation of industrial Peirce-Smith converters (PSCs) used in copper smelting results in several common phenomena such as excessive splashing, slopping, tuyere blockage, and low converter campaign life due to high tuyere line refractory erosion. These phenomena give rise to an inefficient operation from both a process and an energy perspective.This investigation presents numerical and physical models of flow patterns, mixing, solid-liquid mass transfer, and slag-matte phase distribution in a slice model of an industrial PSC. For physical simulations, a 0.2-scale cold model was developed. Water, kerosene, air, and sintered benzoic acid compacts were used to simulate matte, slag, injected gas, and solid additions into the PSC. The two- and three-dimensional numerical simulations of the three-phase system were carried out using volume of fluid (VOF) and realizable k- e (RKE) turbulence models to account for the multiphase and turbulence nature of the flow respectively. These models were implemented using the commercial computational fluid dynamics numerical code FLUENT.Both numerical and physical simulations were able to predict, in agreement, the mixing and dispersion characteristics of the system in relation to various blowing conditions employed in the investigation. Measurement of mass transfer characteristics conclusively indicated that fluid flow in the PSC is stratified. Both blowing configuration and volume of slag in the converters were found to have significant effects on mixing propagation, wave formation, and splashing.The splashing and wave motion in these converters are known to cause losses of metal/matte and potential production time due to the requirement for intermittent cleaning of the converter mouth, resulting in reduced process throughput. To prevent these losses and hence to increase process efficiency, we propose a combined blowing configuration using top lance and lateral nozzles. The numerical simulations were conducted on combined, as well as lateral, blowing conditions. In comparison, the recent results of the combined blowing concept were found to be encouraging."

Jan 1, 2014

By D. K. Chibwe, P. Taskinen

"This investigation consists of a numerical and physical modelling exercise on flow patterns, mixing, solid-liquid mass transfer, and slag-matte phase distribution in a 0.2-scale cold model of an industrial Peirce-Smith converter (PSC). Water, kerosene, air, and sintered benzoic acid compacts were used to simulate matte, slag, injected gas, and solid additions into the PSC. The 2D and 3D numerical simulations were carried out using volume of fluid (VOF) and realizable k-e (RKE) turbulence models to account for the multiphase and turbulence nature of the flow respectively. These models were implemented using the commercial computational fluid dynamics numerical code FLUENT.Numerical and physical simulations were able to predict, in agreement, the mixing and dispersion characteristics of the system in relation to various blowing conditions. Measurement of mass transfer indicated that fluid flow in the PSC is stratified. Blowing configurations and slag volume both had significant effects on mixing propagation, wave formation, and splashing.As a potential process alternative to increase conversion efficiency, we propose a combined blowing configuration using top lance and lateral nozzles. The numerical simulations were conducted on combined as well as lateral blowing conditions, and the results of the combined concept are encouraging."

Jan 1, 2015

By Jose A. Ruiz

INTRODUCTION For heterogeneous non-catalytic solid-gas reactions carried out in fluidized bed reactors, nominal bed temperature is an average value between gas and solid. For highly exothermic reactions, actual particle temperature can be up to hundred degrees above termocouple temperature, since heat dissipate at slower rate than heat generation during some steps of transformation. These micro instabilities at particle level can displace thermodynamic equilibrium from the desired phase to others. This effect is important when relatively narrow range of phase stability exists, like in sulphation roasting of sulphides, where bed discharge contains, in addition to the desired sulphate others metastable phases such as basic sulphate and ferrite. For example, for copper-iron sulphide concentration roasted below 680°C, no cupric ferrite is formed, as shown by Jakob et al., therefore the separate Kellogg's stability phase diagrams Cu-S-O represent the roasting (Fig.l). At temperature above this one, like 760°C for example, phase stability changes from CuSO4 to the basic sulphate Cu0•CuSO4. Solid state reaction between Cu0•CuSO4 and Fe203 also occurs to form cupric ferrite, Cu0•Fe2O3, therefore the Cu-Fe-S-O Kellogg's phase diagram of Fig. 1 (continuous line) represents this condition of roasting. Although both Cu0•CuSO4 and Cu0•Fe2O3 are metastable phases under normal operating conditions of sulphation roasting(1-10% SO2, 1-10% 02 and 680°C, points A and B,Fig. 1) sulphation of these species and particularly ferrite are slow and limited due to diffusional control. This explain why reactors discharge contains several phases: CuSO4, Cu0•CuSO4, CuO•Fe203 plus small amounts of CuS and FeS due to the baclonixed conditions of fluidized bed reactors. The presence of cuprous ferrite due to reducing conditions inside individual particles also occurs, but in very limited extent at the reaction interface, as it is demonstrated by the low soluble iron levels of calcines. On the other hand, the presence of cupric ferrites can be explained only due to the higher temperature achieved by particles, since below 680°C cupric ferrite formation does not take place. [ ] This phenomena also occurs for others non-ferrous metal sulphides concentrates roasted under sulphation conditions, such as 7n, Co and Ni. Therefore the prediction of discharge composition of fluidized beds is not possible based only on kinetics and mechanistics models alone or purely thermodynamical considerations, since all factors should be integrated to describe the overall phenomena. SINGLE PARTICLE TRANSFORMATION Assuming a first order reaction with respect to the gas and zero for solid, such as

Jan 1, 1984

By Amar Prakash, V. M. S. R. Murthy, Ajoy Kumar Singh, Avinash Paul, John Loui Porathur

Weak, layered, and fragile rock mass formation, if not supported properly, is subject to roof failure thereby affecting safety and productivity in underground coal mines. Though Central Mining Research Institute – Indian School of Mines Rock Mass Rating (CMRI-ISM RMR) based well-defined support design guidelines are established, still occurrence of roof failures in underground coal mines is a real matter of concern for mining engineers and researchers. Numerical modeling techniques are successfully used by several researchers by simulating rock mass condition for stability assessment of mine openings. The analysis becomes crucial in case of weak and fragile rock formation. The present research envelops the determination of 31 cases of rock load by CMRI-ISM RMR under different geo-mining conditions followed by the development of modified Rock Mass Rating (RMR), i.e., rock mass rating dynamic ( RMRdyn) by incorporating P-wave velocity as a new parameter. The rock load determined using RMRdyn and numerical models was correlated and found in close agreement. In addition, the deviation in rock load determined by all the three approaches, i.e., CMRI-ISM RMR, numerical modeling, and RMRdyn, was compared with the actual field data. The percentage deviation obtained in RMRdyn and numerical modeling is less compared to CMRI-ISM RMR.

Sep 1, 2023

By J-P Franzidis, E V. Manlap

Froth transportation characteristics have a significant impact on the flotation recovery especially in large flotation cells. The probability of an attached particle being recovered into the concentrate depends strongly on the location where it enters the froth phase. A mathematical model for froth transportation in industrial flotation cells was proposed (Zheng et al, 2004). The model has been tested on the experimental data obtained from testwork conducted at the Western Metals Hellyer lead/zinc concentrator. The results show that the model fits the experimental data well, especially when the froth is well developed (ie for a deep froth). The model fitting can be significantly improved by incorporating an air dispersion factor. However, the models encounter some limitation, ie inconsistent prediction at different froth heights, due to change in froth viscosity which cannot be modelled at the present time.

Jan 1, 2005

By H. Koca, D. O. Aksoy, S. Koca

"In this work, modelling of galena flotation from a disseminated Pb-Zn ore containing 0,83 % Pb was investigated. In the flotation experiments, collector type and amount, frother type and amount, particle size, pH and flotation time were kept constant while solid ratio and depressant dosage were varied. The varied parameters are described with a mathematical model by using statistical design method. Central Composite Design method which is one of the principle response surface methodologies was used in statistical design and analysis of experimental results and process modelling. Two quadratic mathematical models were developed for Pb grade and Pb recovery. Analysis of variance showed that models obtained for Pb grade and Pb recovery were found to be statistically significant in 95 % confidence level. Predicted vs actual values of Pb grade and Pb recovery were also given in graphically. It was found that predicted values represent approximately 96 % and 94 % of actual Pb grade and Pb recovery, respectively. The high values of R-square (R2) indicate that the models developed for two response variables are powerful. INTRODUCTION Lead and zinc metals are widely used in electrical, mechanical, military, metallurgical, chemical, light and medical industries. Sulfide ores are the main resources for lead and zinc minerals (Luo et al., 2016). Froth flotation is the main beneficiation method to achieve lead and zinc sulfide minerals separation. There are two basic approaches to obtain lead and zinc separation by flotation. First, flotation of lead sulphide is carried out and zinc sulphide and gangue minerals are depressed by various depressants and second, the bulk flotation of lead and zinc sulphide is carried out to obtain a concentrate of lead and zinc, followed by the lead and zinc separation. (Yuan et al., 2012). Apart from concentration of lead and zinc by flotation, there are a few literature that reports concentrations of lead and zinc by gravity and hydrometallurgical methods with remarkable success (Goktepe, 2005; Yanbo et al., 2016)."

Jan 1, 2018

By B J. Monaghan, P Zulli, X F. Dong

In the lower zone of blast furnace (BF), slag trickles down through a packed bed of carbonaceous particles in the form of films, rivulets or droplets. During its downward flow, there are significant interactions occurring between slag and other phases, ie gas and packing particles. In particular, the interaction between gas and slag is closely associated with the key flow phenomena such as loading, flooding or channelling in the lower zone of the BF. With gas introduced laterally through tuyeres and slag flowing downward from the cohesive zone to the hearth, this interaction varies throughout the lower zone and in the cohesive zone:  strong cross-flow of gas relative to the slag occurs in front of the raceway  counter-current flow of gas and slag occurs in the upper parts of the lower zone  both strong cross- and counter-current flow of gas and slag around and within the cohesive zone. These interactions can be linked to furnace irregularities, affecting smooth operation and limiting production. Therefore, understanding the gas-slag interaction and its influencing factors remains critical for BF process control. In the current review, the physical behaviour and key influencing factors of gas and slag flow are initially examined. Numerical models describing gas and slag flows in a packed bed at micro-, mesoand macro-scopic levels, are reviewed. Each simulation approach has its advantages and limitations, as well as unique assumptions when used in BF process modelling. Applications of these approaches to gas-slag flow in key identified functional zones will highlight the differences between approaches and modelling outcomes. As no single approach can fully address all aspects of gasslag interactions, an integration of different approaches is considered an effective and beneficial method in providing a more complete understanding of gas-slag flow in the BF process.

Jun 19, 2024

By Du Sichen

The present paper presents the modeling of gas-solid reactions in the case of the decomposition of CaCO3, describing the reaction kinetics on the basis of heat-transfer concepts which could be compared with thermogravimetric results. In the case of solid-solid reactions, the present paper describes the model developed for describing solid-solid reactions from dynamic X-ray diffraction technique, developed by the present team. The model has been successful in explaining the formation of spinel from MgO and Al2O3.

Jan 1, 2006

By O. Konovalov, M. Zhuravkov

In report the results of performance of research works on building up of the connected geoecological and geomechanical monitoring system for modelling and study of a state of mining massif in a region with large-scale underground mining works are discussed. The connected geoecological monitoring requires radical change of methodology of carried out works. To carry out modelling and researches it is impossible without construction of the integrated connected information model of region. Basis of construction of information systems and databases for all sections of Earth Sciences is the set of the approaches, among which, for example, methods of remote sounding, GIS- and GPS-technology etc. Computer technologies are basis on the "supercomputer" technology. The results of performance of some researches for regions of excavation of deposits of potash salts in Russia and Belarus are presented.

Jan 1, 2005

By N. Syna

The selection and preparation of activated carbon suitable for specific application can often involve exhaustive adsorbent preparation and testing. This work examined an emerging technique based on modelling the correlation of observed adsorption capacity to the carbon characteristics as a means of decoupling the roles of carbon textural and surface chemical properties. The insights to these specific roles are then used in tailor-making adsorbents and in establishing the optimal carbon properties for specific application. In this study, a series of activated carbons prepared from bagasse by physical and chemical techniques and commercially available carbons were used in developing the models. These carbons were tested for the adsorption of thiourea-gold (I) chloride complexes. Thiourea is an important alternative gold lixiviant with the potential to address current ecological and safety concerns in the extraction of gold. A four-component model, which correlated the adsorption gold capacity to carbon textural and surface chemical properties, was proposed. The model fit was assessed using multi-linear regression and a partial F-test. This was followed by goodness of fit test achieved through r2, adjusted?r2 and mean square error (MSE) values. The results suggest the adsorption of the gold-thiourea complex is strongly influenced by the internal surface area and pH of the carbon. Tailor-made activated carbon from bagasse showed adsorption gold capacities of up to 198 mgAu/g adsorbent, demonstrating the potential of bagasse as a precursor for gold adsorbents and of this empirical modelling approach in developing adsorbents for metal recovery.

Jan 1, 2014

By Pedro Lobarinhas

In the present work a two dimensional computational model is derived for the description of the time dependent heap leaching process. The flow inside the heap is calculated using the diffusion equations for the flow in porous media. The soil conductivity is also estimated assuming a close packed structure for the grain packing inside the heap. These equations, with appropriate boundary conditions, are solved using a time step splitting scheme. From this calculation the amount of leaching reagent in contact with the ore at any position inside the heap can be calculated and, once coupled with appropriate reaction equations, the amount of metals extracted can be computed and leaching agent concentration updated. The model is tested with available experimental data and the results have shown to be in a very good agreement. The model predicts all the transient features of the data.

Jan 1, 1999

By Roland B. Drugge

The adaptation of pelletizing processes involves the coordination of various factors such as temperatures, gas flows, productions, etc., in order to obtain optimum performance. The consequent problems are considerable due to the large number involved as well as the lack of knowledge of their correlation. It is therefore essential to establish a fundamental theory which takes into account the heat transfer phenomena and the secondary factors which are relevant to the properties of the product. Straight grate and grate kiln pelletizing processes involve heating of a stationary bed on a travelling grate. The heat transfer conditions are therefore of primary importance. This is why this article pays particular attention to conditions which obtain when a stationary bed is heated, in which certain special reactions occur. The conclusions drawn from this penetration of the heat transfer phenomena can be expressed in postulates acting as guide-lines for the establishment of proto-typal processes. Examples of elaborating such processes are also given in the article. Some operational experience at LKAB are also given.

Jan 1, 1975

By D Adhikary, M Khanal, B Poulsen, H Guo

During longwall mining, the ability to reliably estimate the height of longwall mining-induced fractures is critically important in order to effectively manage mining-induced water inflows and aquifer interference. The extent of mining-induced subsurface deformation, fracture, surface subsidence and aquifer interference is largely controlled by local lithological conditions, mining methods and mine layouts. Therefore for the purpose of quantifying mining impacts it is imperative to incorporate the local geology and mine plan parameters in a mine scale model. Two different numerical models based on 3D continuum (Finite Element Method) and 2D discrete (Discrete element method) domains, respectively, are developed to incorporate the mine site parameters of a study site in Australia. Both the models are calibrated and validated with measured mine data. The continuum numerical modelling method is used to investigate the change in permeability of the overburden strata due to underground extraction considering the full 3D nature of geology and mine layout. The discrete element modelling method is used to simulate the initiation and propagation of fractures, allowing the caving of overburden strata to be represented in a realistic manner, and the permeability enhancement of the strata overlying the goaf to be calculated from first principles. Permeability enhancement in the discrete element model is a direct function of the fractures induced in the deformation/caving process; fracture location, extent, aperture and connectedness. From the numerical modelling under the simulated condition for the mine studied here, it is predicted that the permeability of the caving zone from the goaf to five times the mining height increases by almost six orders of magnitude. In the fracture zone to approximately ten times the mining height, the change in permeability is predicted to be around two to three orders of magnitude. At the height of more than ten times the mining height the magnitude of permeability increased is predicted to be one to two orders of magnitude. CITATION:Khanal, M, Poulsen, B, Adhikary, D and Guo, H, 2018. Modelling of height of connective fractures during longwall mining, in Proceedings The Fourth Australasian Ground Control in Mining Conference (AusRock), pp 254–265 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Nov 30, 2018

By B. Csoke

To improve the efficiency of the crushing/screening plant of Bazalt-Középko Kobányák Ltd. (Hungary), the implementation of computer aided quality control is planned. In the current phase, the R & D work aims at the solution of high-priority tasks: 1) Developing the process models, validation of the models by means of industrial measurements. 2) Developing the algorithm for the process control of the crushing/sizing plant. 3) Developing the suitable computer program. To develop the mathematical models of crushing (breakage and probability functions) and to determine the model parameters, the experimental industrial trials were carried out on hydrocone crushers of types Svedala H 4800 EF and H 4800 F, changing the gap size and crushing relatively narrow grain size fractions. Additionally, the parameters of the separation function for the screens operating in the system were determined. It was established that the crushing probability depends on the grain size of the feed (Y) and the gap size (s), so the breakage probability S, for each cone crusher, can be well characterised by the function formed using the ratio of these two parameters: SfYsf==. The distribution of the broken grains can be described reliably using the relative grain size distribution function B(? = x/Y) (usually Schuhmann-Gaudin- or Rosin-Rammler-shaped function). For the screens, the probability function ()2()111/NTx = xd ???---???? is suitable for cha-racterising the separation (where x? grain size, ?-open-area ratio, d ? rectangular screen aperture size, N ? number of appearance of the particle on a screen surface). Using these functions, the comminution circuits were described by means of vector/matrix equations. The calculations and measurements verified the reliability of the developed models and algorithms.

Jan 1, 2014

By H. Lind

The consumption rate of Estonian oil shale increases due to the growing rate of exported oil-shale oil. Estonian oil shale active resources in working mines are predicted until the year of 2025 when the consumption rate stays at the level of 12 tons per year. If the consumption rate stays at that level, new mines should be opened after 20 years, if the rate increases new mines should be opened even earlier. The predicted active resources in Estonian oil shale deposits are calculated for 60 years, based on the parameters of the power station technicoeconomics. As a result, there is a need to open new mines. One of the potential areas for surface mine is the Tammiku-Kose prospect area (Rheinbaun Engineering, 1996).

Jan 1, 2006