Search Documents

By D. Komljenovic

This paper presents the use of queuing theory in truck fleet sizing in open pit mining operations. In the first part, a methodology is developed by which an optimal ratio between the size of shovels and trucks may be determined. The second part gives the fundamentals of the queuing theory and its application in analyzing shovel and trucks. The stochastic character of involved times (loading, truck cycle, etc.) is taken into consideration. Using the above theory, the optimal truck fleet size in actual conditions may be estimated. Minimal operating cost of the working system (shovels-trucks) represents the basic requirement of the analyzed process. The methodology was validated on Caterpillar front-shovel ? truck combinations.

Sep 1, 2004

By D. Komljenovic, J. Paraszczak

"This paper presents the use of queuing theory in truck fleet sizing in open pit mining operations. In the first part, a methodology is developed by which an optimal ratio between the size of shovels and trucks may be determined. The second part gives the fundamentals of the queuing theory and its application in analyzing shovel and trucks. The stochastic character of involved times (loading, truck cycle, etc.) is taken into consideration. Using the above theory, the optimal truck fleet size in actual conditions may be estimated. Minimal operating cost of the working system (shovelstrucks) represents the basic requirement of the analyzed process. The methodology was validated on Caterpillar front-shovel – truck combinations. IntroductionLes camions représentent probablement les équipements les plus utilisés pour le transport dans les mines. Le coût de transport peut parfois atteindre 50 % du coût d’exploitation et 25 % du coût total (Collins, 1996; Kennedy, 1990). Il est critique pour une opération minière de le minimiser. C’est un sujet qui a été traité relativement souvent dans la littérature spécialisée. Les premières applications informatiques dans l’industrie minière ont été développées pour la détermination du temps de cycle des camions.Au fil des ans, la capacité des camions a augmenté et, aujourd’hui, les plus gros manufacturiers de camions (Caterpillar, Orenstein & Koppel, Euclid, Terex, Unit Rig, Komatsu, Liebherr, Belaz) fabriquent des camions ayant une capacité de 200 t et plus dont quelques-uns d’une capacité supérieure à 320 t. Pour les plus gros camions, le coût d’opération peut même dépasser 240 $US/heure dans des conditions modérées (Caterpillar, 1999). Les mines modernes de grande taille peuvent posséder plusieurs dizaines de camions dans leur flotte."

Jan 1, 2004

By Rae Simpson

"This paper provides insight into the elements of computerized maintenance systems and identifies some of the factors which influence their effectiveness. The paper describes what to expect in terms of benefits and problems, the nature of effective computerized maintenance systems, some of the systems available and a process for selecting a system for your organization.The paper will be of interest to mining company executives and maintenance or computer personnel who are involved with, or contemplating, the introduction of computerized maintenance systems.IntroductionIn the drive to improve profits in an increasingly competitive world, mining companies are now using some of the most sophisticated equipment available. The introduction of large and complex equipment such as drag lines, haul trucks, and high-capacity milling equipment has resulted in order-of magnitude improvements in productivity. However, this productivity is now, more than ever, dependent on the continued operation of the equipment. A prolonged breakdown of a critical component can have a significant effect on companies' bottom-line results. As more technologically sophisticated equipment is introduced, maintenance costs become an increasingly important factor for mining companies.Fortunately, sound modern management principles are being used more often in maintenance departments. Progressive maintenance managers are turning to computers to help manage the myriad of details associated with running a maintenance organization.Some maintenance departments in larger mining operations would find it difficult, if not impossible, to operate at current productivity levels without their computer systems. Unfortunately, some other companies gain very little from their computerized systems and would probably be better off without them.How can you be sure that computerized maintenance will improve the maintenance picture for your company? What benefits can you expect to achieve? How do you go about selecting a maintenance system for your organization? These are some of the questions facing industry managers and executives today. In the following sections, we will provide some guidelines and insights to help address these challenges."

Jan 1, 1988

By Earl M. Frizzell, Joseph A. Dos Santos

"Conventional belt conveyors offer an economical method for transporting materials at inclines approaching 18 degrees. Sandwich belt conveyors may be used to transport materials at high angles, to 90 degrees, while maintaining the most positive features of conventional conveyors, that is, (1) proven conventional conveyor components. (2) high-speed operation, (3) high-single lifts, (4) continuous cleaning of the smooth surf aced belts, and (5) easy, quick repair by hot or cold vulcanizing. The present discussion deals with the evolution and development of sandwich belt conveyors to satisfy the operational requirements of the open-pit mining and materials-handling industries.A look at the limitation s of conventional shallow angle conveyors reveals the conveying angle at which a cover belt is needed. The relationship between conveying angle and hugging pressure by the cover belt is developed mathematically. A study of past and current developments reveals that the Loop Belt and Beltavator are significant advances in the state-of -the- art. The vertical radius of curvature constraints are discussed in detail and the constraint equations are developed. Criteria for future development is established as a basis for evolutionary developments in sandwich belt high-angle conveyors. The use of intermediate drives in sandwich belt high-angle conveyors is investigated.IntroductionIn connection with U.S. Bureau of Mines funded–study relating to continuous haulage from deep open-pit mines, the authors had occasion to study extensively, the use of high-angle conveyors as a means of hauling mine products continuously from the pit while maintaining optimum mine slopes to minimize total excavation. The many possible high-angle conveying methods investigated included skip hoists, bucket ladders, pocket belts, fin belts, cleated belts, sandwich belts, pipe belts, screw conveyors, slurry pipe lines, and others. The sandwich belt high-angle conveyor appeared to be the most operationally-appropriate and economical solution for the mining industry. The present paper presents the theory, historical development, state-of-the-art, and evolutionary developments in sandwich belt high-angle conveyors."

Jan 1, 1983

By Frands Jensen, Erik Eriksson

"Haldor Topsoe’s Wet gas Sulfuric Acid (WSA) technology has been proven in a number of industries converting sulfurous gasses to concentrated sulfuric acid. In the metallurgical sector plants have been delivered handling SO2 gasses from copper, zinc, molybdenum, and lead ore. WSA has high energy efficiency and allows autothermal operation from feed gas strengths as low as 3.0% (volume). Recent variants of the WSA process have been developed that allow very high SO2 conversion efficiencies of over 99.9%. This paper describes the technology, its advantages and industrial achievements.The WSA technology is well suited for production of concentrated sulfuric acid of commercial quality from the sulfur dioxide (SO2) gas that inevitably emerges from the metallurgical processes. Some of the features that make the WSA technology attractive for this application are:- High SO2 conversion rates (up to 99.99%)- High energy efficiency leading to a low autothermal point (~ 3% SO2)- No need for gas cooling and drying to produce concentrated sulfuric acid- No production of waste- Very low cooling water consumptionThis paper describes the WSA process principles and highlights features that make this process attractive for handling SO2 gases."

Jan 1, 2014

By Peter Dunn

London Mining Plc. is developing an iron ore project in Sierra Leone. The Marampa mine is a brownfield development effort, having been operated by others since the 1930?s until low iron ore prices forced closure. Peak production of the prior operations was approximately 2.5 Mt/a in the 1960?s. Following civil war In Sierra Leone, London Mining obtained the mining license in 2006 and commenced activity to redevelop the Marampa mining and processing operations. There are two ore deposits of interest for development as well as a resource comprised of tailings from previous operations. Phase 1 will be the mining of the tailings resource with a 1.8 Mt/a plant producing 66 per cent Fe sinter feed, expandable to 3.6 Mtpa. The first phase is currently under construction with start-up expected in the third quarter of 2011. Phase 2 will include development of the two ore deposits. These deposits are weathered secondary ore and the main primary ore resource. Phase 2 is currently under prefeasibility study and targets a final 16 Mt/a concentrate production in 2016. Outotec partnered with London Mining on testing efforts to develop the magnetic separation flowsheet using SLon VPHGMS machines to process the Marampa material. This paper explains the flowsheet development process for the weathered and primary ores commencing with sample characterization and detailed parametric investigations to evaluate SLon performance. The final flowsheet configuration is explained, including expected circuit recoveries and concentrate quality, clearly demonstrating the performance of the SLon VPHGMS circuit in iron ore tailings treatment applications.

May 1, 2011

By M. Melashvili

"The SART Process was introduced to the gold mining industry in the late 1990s and several commercial plants were built over the following 10 to 15 years. Most of these plants have operated well and enjoyed plant performance in line with expectations based on the underlying chemistry of the process. However, two of the plants have consistently suffered from high reagent consumption and poor product quality, which has motivated a re-examination of the process fundamentals.The primary product in the SART Process is synthetic chalcocite, Cu2S, and the stoichiometric requirement for sulphide addition is therefore 0.5 moles per mole of copper in solution. It is shown here that when sulphide is added in excess of this amount (generally with the express purpose of lowering residual copper in solution) it not only fails in its intended purpose, but can have a negative impact on the settling and filtration characteristics of the copper sulphide product, as well as the final product quality. Moreover, these problems are exacerbated the higher the selected operating pH in the pH 3 to 5 range, and the longer the residence time of slurry in this pH range.Freshly precipitated copper sulphide compounds are known to be very fine and slimy, with poor settling and filtration characteristics. It is therefore important to engineer the SART Process in such a way as to encourage growth in the average particle size of the precipitate crystals. This was the rationale for recycling a portion of the thickener underflow to the primary reactor when the SART Process was developed. However, the theory of crystallization predicts that compounds of extremely low solubility, such as chalcocite (Cu2S: KSP = 10-48), are unlikely to experience growth in crystal particle size through heterogeneous nucleation and are more likely to grow through aggregation. On this basis, the SART practice of recycling a portion of thickener underflow to the primary reactor in the commercial plants is unlikely to lead to crystal growth via nucleation, and could in fact promote breakage of the aggregates that formed in the thickener. Recycling thickener underflow could therefore have the opposite of the desired outcome. The result of testwork that is reported here lends credence to this theory and questions the wisdom of recycling in the SART process."

Jan 1, 2016

By K. Fumimura, T. Nishimura, M. Kohno

When applying numerical techniques to problems of deformation of rocks, the body to be analyzed is divided into small portions and is represented by a system of points linked together with the neighboring points. In the so-called distinct element modeling, the particles in the system are directly connected by springs and a contact law is used to specify forces on the particles to simulate the mechanical behavior of the body. The relation between the spring stiffness and commonly measured material properties, such as the Young’s modulus and the Poisson’s ratio, should be known a priori. For the general case of a random packing of particles, the relation is found by means of a process in which a system of particles and set of parameters are prepared to simulate a set of specimen-size test and then the input parameters are chosen to match the measured material properties. The validity of the numerical result is only demonstrated by comparison model behaviors with the measured responses. The process includes trial-and-error manner, therefore, no model is complete to reproduce the mechanical properties. In this paper, a relation between the spring stiffness and the elastic constants; the Young’s modulus and the Poisson’s ratio, is introduced to a three-dimensional distinct lattice spring method. Two simple numerical examples will be presented to demonstrate the performance of the method with the procedure for the determination of the spring stiffness.

Jan 1, 2015

By Jeffrey S. Loen

"The size, morphology, surface texture, per cent mineral inclusions, and bulk chemistry (Ag and trace elements) of placer gold particles vary systematically according to transport distance from source rocks, as indicated by a case study from the Pioneer district in Montana, and data from the literature.Placer gold samples show a rapid decrease in mean and maximum particle length, and a gradual increase in flatness index (F.I. = length + breadth/2 x thickness) from about F.I. = 2 in bedrock deposits, to as much as F.I. = 45 over transport distances of several tens of kilometers. Primary mineral inclusions (mainly quartz, sulfide minerals, and iron oxides) are mostly lost within 10 km to 20 km of fluvial transport; this loss of quartz and trace metals causes a major change in bulk chemistry of placer gold samples. Following the loss of inclusions and the development of flakes, edges of flakes are folded and the mean number of folded edges in a sample increases with distance. When viewed under the electron microscope, surface textures of placer gold particles also become progressively more pitted with distance.Analysis of placer gold particles can be used to help evaluate placer source-area problems (placers with no known hard-rock gold deposits). For example, coarse placer gold containing abundant inclusions indicates a proximal source. On the other hand, an abundance of thin flakes is evidence of a long transport distance and possible reworking from older gravels. Placer gold characteristics can therefore provide information on distance to mineralized source rocks."

Jan 1, 1995

By Phil Keele

Canadian Natural ??Business Approach Efficiency and Execution ?Drive to be the low cost producer ?Focus on exploitation ?Augment exploitation with strategic acquisitions ?Maintain Discipline and Strength while Delivering on the Defined Plan ?Maintain flexibility and control allocation of capital ?Strive for balance (asset portfolio) ?Maintain financial strength

May 1, 2010

By H. Tokuichi, H. Ohta, K. Sato

A physical separation technique of tantalum capacitors from waste print circuit boards is in a development phase in Japan. A tantalum capacitor consists of an exterior coating of resin (epoxy resin) mixed with silica fillers and a sintered tantalum compact, the latter of which is impregnated with manganese dioxide (MnO2). In this study, we developed a process to recover tantalum oxide from tantalum capacitors by roasting, fluidized-bed separation, and hydrometallurgical treatments. Our aim was to recover tantalum oxide with high purity from waste tantalum capacitors. In this process, the exterior resin was selectively decomposed by roasting without oxidation and pulverization of the sintered tantalum compact at a proper temperature. Furthermore, the decomposed resin with silica filler on the surface of the sintered tantalum compact was mostly removed by fluidized-bed separation. The surface of the sintered tantalum compact was cleaned by ultrasonic irradiation in an acid solution to remove the MnO2 paste as well as residual Ag paste and silica filler. The impregnated MnO2 in the sintered tantalum compact was leached as Mn(II) ions using acid solutions via the reducing effect of tantalum metal to recover purified tantalum oxide.

Jan 1, 2014

In terms of throughput, product size and power draw, grindcurves provide a measure of the response of the mill performance to changes in mill operating parameters – usually mill filling. Since presenting the principle in 2006, grindcurves have been applied to a wide variety of AG and SAG mills. Mill response has been carefully measured for variations in: mill fillings from 15% to 40%; throughput by a factor of two; power by 50%; feed blends; liner profile; and mill speed from 65% to 80%. The outcomes of this extensive program of work are presented as the measured grindcurves. How they have enhanced our understanding of mill operation and how they may be applied to practical daily mill control is discussed. It is proposed that this technique provides a viable and valuable mill control tool that may find extensive application in the industry.

Jan 1, 2012

By D. H. Wilton, C. I. Godwin, A. J. Sinclair

"A computer-based drill core logging system used to examine computer-based drill core logging system used to examine copper mineralization hosted in volcaniclastic rocks at the Sustut deposit, in north-central British Columbia, has many advantages over conventional logging. A relatively simple format for recording geological information allows complex, but significant, geologic parameters to be recorded consistently, accurately and quickly. This allows graphical presentations and statistical interpretation of the data that lead to rapid and detailed analyses of the control and distribution of copper mineralization at the Sustut deposit.Ore mineralization is concentrated in lens-shaped bodies and is zoned from native copper cores through chalcocite, bornite and chalcopyrite to a margin containing pyrite. Chalcocite, the principal mineral, occurs mainly in volcaniclastic agglomerate characterized by a specific grain size, a narrow range of matrix and fragment proportions, and a specific sorting.IntroductionWe report here a case history of the practical application of a detailed computer-based drill core logging system applied to the Sustut copper deposit of Falconbridge Nickel Mines Ltd. The Sustut deposit, containing a few tens of millions of tons of ore grading slightly more than one per cent copper (Church, 1975), is in north-central British Columbia about 875 km north of Vancouver. The deposit is contained in volcaniclastic sedimentary rocks of the Upper Triassic Moosevale Formation, which is within the upper part of the Takla Group (Monger and Church, 1977) of basic volcanic and associated sedimentary rocks. Takla rocks occur along the eastern fringe of the Intermontane Belt in the Canadian Cordillera (Sutherland Brown et al., 1971). The deposit consists of concentrations of copper sulphides, native copper and pyrite within pods and lenses parallel to bedding in the host sedimentary rocks. These sheet-like deposits exhibit a mineral zoning from pyrite at the fringes successively through chalcopyrite, bornite and chalcocite to native copper at the core (Wilton, 1978; Wilton and Sinclair, 1979)."

Jan 1, 1980

By Tiago Webber

Mineral resources estimation comprises basically two stages: (i) data collection, analysis and validation, and (ii) the estimate of the grades and tones of the deposit. The most applied method for sampling survey in mineral exploration is rotary drilling with total core recovery. High drilling costs make borehole geophysical logging attractive to obtain the maximum possible information from each borehole helping to understand ore intersections, grades variability, limits of the deposits and continuity of the mineralization. In this paper, borehole geophysical logging was used to enhance the sampling survey and to derive grades indirectly (without lab analysis) with a certain amount of uncertainty. Indicator kriging (IK) was the estimation method used, as it is adequate to incorporate the uncertainty associated with the geophysical derived grade values. The results were compared to the conventional procedures using grades determined via core analysis and ordinary kriging. Borehole geophysical logging proved to be efficient to define the essential features of ore deposit and IK provided results with acceptable precision and accuracy if compared to the model provided by ordinary kriging.

Oct 1, 2009

By K. A. Foxford, D. A. Polya, R. P. B. Hebblethwaite

"A variety of techniques for interpreting systems of mineralized echelon veins is discussed using examples from Panasqueira and SW England. As aids for exploration and mine planning, and as a tool for ore reserve evaluation, particular emphasis is placed upon establishing vein morphology in the third dimension, and upon understanding vein dilation and propagation histories. Special attention is also placed upon investigating vein thickness and spacing relationships, and upon understanding vein segmentation. The relevant data for these techniques are all readily obtainable from stope and core sources.The vein systems described were produced in host-rock sequences that underwent incremental dilation in response to fluid overpressures, and individual veins display evidence for mineral infilling synchronous with their propagation. Echelon vein geometries resemble those reported from igneous dikes and sills: vein walls are stepped and veins are divided into segments arranged in echelons. Changes in vein geometry, from arrays to composite structures, reflect both the evolution of the vein system with time, and the changes occurring between the margins (arrays) and core (composite structures) of the system. These characteristic morphologies can be used to infer vein propagation and fluid injection trends, proximity of neighboring veins, and the shape of the system as a whole. The techniques are applicable to systems of echelon veins produced under tensile conditions in rocks where faults were absent at the time of vein formation, and where subsequent deformation has been minimal. Further work is required to assess the application to echelon vein systems formed within more complex tectonic regimes."

Jan 1, 1995

By Emmanuel Henry

Few industries are more challenged by technical, political and market uncertainty than the mineral industry. Technical uncertainty in mining projects concerns the volume and quality of the resource, the processing and metallurgical parameters, the production efficiency, etc. Political uncertainty can take diverse forms and includes possibilities of change in the taxation system, risk of expropriation and risk of social instability or war. Market uncertainty depends on the market structure of the commodity; price and demand volatilities are the most common sources of market uncertainty. Intensive feasibility studies including these uncertainties and risks are required before investing up to billions of dollars in the development of a new mine. Current feasibility studies are often based on a single (or very few) mean geological scenario and simplified price profiles along the life of the project. Yet, strategic decisions, for example, the level of investment in exploration and metallurgical studies, adoption of a flexible, but more costly, mine design and risk-sharing strategies are made at an early stage. Over the last twenty years, Real Option theory has emerged as a unifying framework enabling strategic planning under uncertainty (technical, political and market). Geological and market variances (second order statistics) are pivotal for proper uncertainty quantification. The present study explores the possible benefits of combining modern geostatistical tools and financial techniques for project management. More precisely, the following questions are examined: Is geological uncertainty really an issue? Can we combine geological and price uncertainties in a unifying Real Option framework? A simple example aimed at illustrating the concepts involved is presented.

May 1, 2003

By Neda Zangeneh

Production efficiency from low permeability gas shales require techniques to optimize hydraulic fracture completions. One completion strategy involves the simultaneous or near simultaneous hydraulic fracturing of multiple horizontal wells to maximize the fracture network area and stimulated reservoir volume. However, changes to the in-situ stress field caused by an earlier hydrofrac on subsequent hydraulic fractures are not accounted for in traditional hydrofrac design calculations. Presented here are the results from a set of transient, coupled hydro-mechanical simulations of a rock mass containing two wellbores using the discontinuum-based distinct-element method. The results demonstrate the influence of stress-shadows generated by a hydraulic fracture on the development of subsequent hydrofracs from an adjacent well. It is shown here that these interactions have the potential to change the size and effectiveness of the hydraulic fracture stimulation by changing the extent of the hydraulic fracture around the secondary well.

Jun 1, 2012

By Clifford R. Stanley

"Abstract - Reproducible analyses of representative samples for elements affected by the ‘nugget effect’ is vital in mineral exploration and reserve estimation. Unfortunately, a priori knowledge of what size of sample is representative from a mineral deposit or its dispersion train is commonly unavailable. Use of the simple ‘equant grain model’ and the NUGGET computer program facilitates determination of the sampling parameters that will produce appropriate samples for chemical analysis. Based on Poisson statistics, NUGGET can be used to determine: (1) the ‘effective grain size’ of nuggets in a suite of samples, (2) the level of representativity (precision) that samples with specific nugget grain sizes will exhibit, and (3) the sample size required to obtain a specific level of sample representativity. The program determines these sampling parameters for any element in any mineral, including solid solution minerals (e.g., Au, electrum, Pt-Fe alloy, diamond), and allows assumption of a variety of grain shapes (e.g., ellipsoids, cylinders, parallelepipeds, etc.). Furthermore, it outputs the effective number of nuggets expected in a sample of given size and the percentile confidence bounds on any element concentration affected by the nugget effect.Results from the software allow determination of effective sampling and analysis strategies that minimize the nugget effect in geochemical analysis/assaying used in mineral exploration and reserve estimation. This will not only improve confidence in the resulting exploration results or reserve estimate, but will facilitate estimation of the error associated with any subsequent interpretation or reserve figure."

Jan 1, 1998

By G. Chung, W. Petruk, D. Doyle

"A series of tests was conducted on Carol Lake iron ore, which contains about 38 wt % hematite and 19 wt % magnetite, to establish a flowsheet for high weight yield by spirals. A pilot plant was run in parallel with the IOC concentrator to simultaneously test two flowsheets on the same feed material. The IOC concentrator processes the ore by treating it directly in a spiral circuit, whereas the pilot plant used a laconite flowsheet which had a low intensity magnetic separator (LIMS) at the front of the circuit. The LIMS produced a non-magnetic (hematite-rich) product and a magnetic (magnetite-rich) product prior to treating the ore in two separate spiral circuits. The results from about 150 tests were evaluated by chemical assays and materials balancing techniques, and a selected test was evaluated by the image analysis and materials balancing technique.The image analysis and materials balancing evaluation showed that a higher hematite recovery was obtained by the hematite spiral circuit than by the IOC concentrator spiral circuit, and a slightly lower magnetite recovery by the magnetite spiral circuit than by the IOC concentrator spiral circuit. Furthermore the hematite recovery by hematite spirals was higher than the magnetite recovery by magnetite spirals. The better hematite recovery was achieved because a higher liberation is obtained for hematite than for magnetite, and the hematite in the Carol Lake ore has a narrower grain size range than magnetite. Consequently, better spiral control can be maintained when processing a hematite-rich spiral feed, than when processing a magnetite-rich spiral feed. Since the feed to the IOC concentrator spiral is a mixture of hematite and magnetite, an intermediate hematite recovery is obtained by the IOC concentrator. INTRODUCTIONA study of characteristics that affect recoveries of hematite and magnetite from the Carol Lake ores of the Iron Ore Company of Canada (IOC), was carried out in connection with a project on designing a flowsheet for maximum weight yields. A pilot plant using a taconite flowsheet was run in parallel with the existing IOC concentrator which uses a classical iron ore flowsheet. About 150 tests were run to compare iron recoveries and weight yields by the spirals in the two flowsheets. Better iron recoveries were indicated for the taconite flowsheet, but the results were not consistent and the reasons for the observed differences were not apparent. The material from one test (Test 61), which gave better iron recoveries for the taconite flowsheet, was evaluated by the image analysis and materials balancing technique. A suite of samples was taken from 5 points in the IOC concentrator and from 8 points in a pilot plant (taconite flowsheet). The samples were analyzed to compare the behaviours of the Carol Lake hematite and magnetite in the spirals of the two flowsheets. The results provided a sound basis for further tests in designing a high weight yield flow sheet."

Jan 1, 1993

When underground extractions occur beneath or in close proximity to valleys and other forms of irregular surface topographic features, the observed vertical subsidence at the base of the valley is less than that would be expected in flat terrain, while the valley sides are observed to move inwards towards the centre of the valley. The reduction in subsidence at the bottom of the valley is referred to as upsidence and the convergence of two sides of the valley is termed valley closure. Research coupled with mining experience over the past few decades in Australia has identified valley closure subsidence as being very real and quite significant part of non-conventional subsidence effects, which always have adverse effects on natural features, with risks to water flow and aquatic ecosystems. In this paper, we describe a threedimensional numerical modelling approach to investigate the geomechanical mechanisms which cause the valley closure subsidence behaviour, using distinct element computer package 3DEC. The calibrated models can accurately predict subsidence and horizontal displacements, and they are subsequently extended to specifically study the effect of a series of mining geological and topographical factors on valley closure subsidence phenomenon, by developing a range of hypothetical scenarios. Previous analyses carried out using a two-dimensional numerical modelling technique have identified the potential parameters that affect the valley closure subsidence, which mainly include longwall position, horizontal stresses, valley sloping angle, and valley shape. Further three-dimensional analyses regarding the relationships between the major factors and valley closure subsidence are discussed in this paper. The results have significant implications with regard to both the underlying mechanisms that lead to this nonconventional subsidence behaviour, and how these should be incorporated in future valley closure subsidence prediction and mitigation of its impacts on natural features.

Jan 1, 2015