Search Documents

By M Armony, J F. C L Costa, J C. Koppe

Linear kriging methods fail to estimate local grades and local reserves for highly skewed variables due to three reasons: it uses only part of the available data information, it has no physical principles associated with the mathematical proceedings to minimise the problems, and it works with continuous variables while data are discrete because of the limited precision of any equipment. The techniques commonly used up to now to solve the highly skewed variables estimation problems were unsuccessful, unless mathematical rigor is abandoned. Even the method of multiple indicator kriging is not satisfactory. This paper presents a new framework, the field parametric geostatistics (FPG) that transforms noisy variograms into well-behaved variograms and justifies mathematically empirical procedures commonly used, as trimming or capping arbitrarily very high values. The method, when applied to non-skewed variables, yields similar results to classic kriging. The methodology is illustrated in a case study at a gold deposit in Brazil, and compared to results obtained by usual techniques.

Sep 26, 2011

The cloud transform simulation is a technique used to generate realisations of correlated attributes. The principle of this technique is to partition the scatter diagram between the input attributes into several interval classes and derive for each class a cumulative distribution function (CDF) using only the data pairs of the scatter diagram in that class. The simulation then proceeds like a traditional Monte Carlo simulation but using the derived CDF in each class.The cloud transform simulation has the advantage of being a simple and fast technique that is able to reproduce the uncertainty in the relationship between the two input attributes regardless of whether the relationship is linear or not. However, the technique has three major drawbacks. First, the number of classes has an impact on final simulated results and there is not a robust approach to define such number. Second, simulated results exhibit some saw tooth variations that are unrealistic and third the reproduction of the spatial continuity of the simulated attribute is not guaranteed.This paper presents the use of the quantile regression approach in order to overcome the first two disadvantages of the cloud transform simulation. Details of the technique and how it can be coupled with the cloud transform simulation are presented along with a real case study.

Sep 26, 2011

By Q Qian, M Karakus, J Wang, N Melkoumian, D Goodchild, Q Sui, C Wang

Probabilistic analysis is widely used in slope stability assessment because of its capability to take into account uncertainties associated with geotechnical parameters. There are, however, still many unresolved practical issues to the application of the method. One of these issues is the integration of assessments when multiple sets of fractures are considered for the entire slope. In the present study, this issue is addressed within the framework of Monte Carlo simulation. Probability of failure for pairs of fracture sets are assessed independently and the final assessment of the probability of failure for the entire slope is done using conditional probabilities. Prominent Hill open pit, in particular, the problematic fresh Bulldog shale formation, was used as a case study to test the method proposed. The SWEDGE program of the RocScience software package was used for the analysis.

Sep 26, 2011

By R Razanatsimba, A Galli

A mining company has entered into contractual commitments to supply certain quantities of its product to clients in each time period. Failure to fulfil its obligations would damage its reputation as a reliable supplier. The planned production would normally be sufficient unless major problems such as highwall failures, roof falls, natural disaster, or strikes that interrupt production or delays in obtaining environmental permits, occur. To overcome these difficulties, the company may be able to obtain more of the commodity: from its strategic stockpile (if there is sufficient there); by buying it on the spot market; or for some commodities such as gold and uranium, by leasing it. Two sources of uncertainty, the spot price of the commodity and the occurrence of serious incidents that interrupt production, are considered. For simplicity it is assumed that the commodity price follows a geometric Brownian motion, that productions incidents have a binomial distribution, that they are independent from one time period to another and that they do not affect prices. A branching tree structure is used to model the joint evolution of prices and production incidents. The aim of this research project was to evaluate the companyÆs options for fulfilling its contractual commitments. When evaluating these options, multi-stage stochastic programming with recourse was used to solve this problem because in addition to providing the dollar value of the project, it gives decision-makers the æroadmapÆ to reach the optimal value. That is, how much material should be bought/sold on the spot market, how much to add to/subtract from the stockpile, how much to lease and in the worst case what the shortfall would be, for each time period as a function of the current spot price and breakdown status. A case-study over a five-year period is used to illustrate the proposed procedure. In this study the production would be sufficient to cover the companyÆs contractual commitments unless production incidents occur. Focus is on evaluating the impact of the quantity initially in the stockpile on the value of the project when the commodity market is tight, that is, when only small quantities can be bought on the spot market or by leasing.

Sep 26, 2011

The paper provides an overview of the development and content of APCOM symposia since their inception in 1961. This development is set within the context of the academic environment in which much of the research published at APCOM symposia is conducted. This environment includes rigorous external assessment of research quality at national levels and against international benchmarks, which carries significant consequences for funding and for the continued existence of the academic minerals industry disciplines. Adverse consequences that undermine the academic provision of mining and resource engineering will impact negatively on the minerals industry by reducing the level of research and education resources available to it.The traditional reductionist approach to research is contrasted with the inter-disciplinary approach that is required to address the complex systems problems that the minerals industry will have to address to ensure its technical, economical and environmental sustainability. Such an approach offers a means of raising research profiles, of establishing mutually beneficial academic-industry engagement and provides an opportunity for APCOM to take a lead role in identifying the research and development agenda for the future

Sep 26, 2011

By B F. J Kelly

Multiple-point simulation is a family of methods used to generate conditional facies simulations that present the same spatial dependence as a user-defined training image. The main advantage of using a training image is the easy integration of geological concepts in numerical subsurface models. The strength of multiple-point methods is that it allows making the best use of geological expertise when data are too sparse for inferring a model of spatial continuity. However, despite being appealing to geologists, the multiple-point methods are often not used in practice due to the difficulty of building complex 3D training images that contain realistic features.A new framework for multiple-point simulation is presented that involves small and simple training images. The main idea is to use random transformations to expand the range of structures available in the simple training image. The training image is not regarded any more as a global conceptual geological model, but rather a basic structural element of the subsurface. Complex geological structures are obtained, whose spatial structure can be parameterised by adjusting the statistics of the random transformations, based on field data or geological context. In most cases, such parameterisation is possible by adjusting two numbers only. One advantage is that the training images are so simple that they can be easily built even in 3D. The method allows building 3D models that reproduce shapes corresponding to a desired prior geological concept and are in phase with different types of field observations, such as orientation, facies or geophysical measurements. The flexibility of the approach allows for naturally integrating different types of data. This is demonstrated on an example involving seismic data assimilation. Examples where complex 3D folded structures are built from a simple training image and a few orientation data are shown.

Sep 26, 2011

By Q Wang

Net present value (NPV) has been the main criterion in mining project evaluation and production plan optimisation. It is a common practice in NPV calculations to take into account only the tangible cost items in the corporate accounting system. The costs incurred due to damage to the local ecosystem and environment are left out and treated as externalities. Such practice on the one hand underestimates the total cost of a mining operation and, on the other hand, ignores the difference in ecological/environmental impacts between different mine plan alternatives. This paper attempts to measure the ecological costs of mining and to incorporate such costs into the evaluation scheme in mining project evaluation and mine planning. The ecological costs of mining are estimated based on the areas and types of land (eg forest land, grazing land, crop land) damaged by mining activities and the carbon emission from fossil fuel consumption. For the land areas (and their associated ecosystems) directly damaged by mining activities, the ecological costs are measured by the prevention and restoration costs, the lost values of eco-services of the damaged land, including bio-product production, soil erosion control, water run-off control, oxygen release, air pollutant absorption and formation of soil nutrients. The ecological cost of carbon emission from the energy consumption of a mining operation is measured by the cost of cultivating the equivalent forest land required to uptake the emitted carbon. A case study on a large-scale open pit mining project is presented to demonstrate the influence of ecological costs on the mine planning outcome.

Sep 26, 2011

By P Hatherly, A Melkumyan, D Wyman

Chemical species distributions in rock sequences can be modelled using Gaussian processes (GPs) to predict the weight percentage between known data points. This method can be improved by dividing the modelled area into regions of different mineralogy. In the context of many iron ore deposits, marker shales give the initial guidance of where these regions may lie. These shales give rise to distinctive peaks in the natural gamma downhole logs, which are conventionally identified by hand. A GP method has been developed to automate this task.Once the appropriate boundary has been located on the basis of shale occurrences, chemical assays from exploration drill holes are used to find the exact boundaries of interest. These boundaries divide the drill hole stratigraphy into regions of different mineralogy, each of which display their own distinct correlations between the main elements and oxides (Fe, SiO2 and Al2O3). Iron ore shows a negative correlation between Fe and Al2O3, but in banded iron formation (BIF) there is a positive correlation between these species. Similarly, SiO2 and Al2O3 have a positive correlation in the shales and ore but a negative correlation in the BIF. Correlations obtained within ore-, BIF- and shale-dominated regions are therefore better than those obtained using the entire log and can be used to improve the results obtained when modelling.

Sep 26, 2011

By D Cai, E Baafi

Various enabling technologies and systems have recently been introduced in the mining industry. Such systems include miner-bolters with concurrent cutting and bolting, employment of two shuttle cars, adoption of self drilling bolts and utilisation of continuous haulage or monorail systems. The paper discusses how Flexsim discrete event simulation models were used to assess how these proposed mining systems can best be used to achieve and sustain high roadway development rates of 10 m per operating hour for 20 h/day.

Sep 26, 2011

By E Cunha, D Roldao, TáTamantini, J F. C L Costa, T Webber, C Monteiro

The main problem of geostatistical simulations is to reproduce the real characteristics of the mineral deposit when there are few conditioning data. In this case, the geological interpretation, based on structures, geophysics and geological maps, plays the major role in terms of resource evaluation.The paper investigates the use of principal component analysis (PCA) to help classify multiple geostatistical categorical simulations considering the similarity among them and the interpreted model, herein referred as training image (TI).An iron deposit with four lithological facies, in a specific horizontal plane, was used as case study for two-dimensional (2D) simulations. Sequential indicator simulations (SIS) was used to obtain ten equally probable 2D images, which were generated and compared against the interpreted model (TI).In this study, 11 random variables are comprised by the vector of ten SIS realisations plus the training image (TI). The similarity between the TI and each SIS realisation is measured using the angle between the TI and SIS vectors projected in the new coordinate systems (PCA factors). It is possible to order the SIS realisations from the highest to the lowest PCA factor taking into account the TI similarity.The visual analysis of the images shows coherence in the PCA results as a method to measure similarity. This methodology can be used as a post-processing tool to rank simulation results and to select the ones which better represent the prior known geology or to compare the quality of different methods of simulation.

Sep 26, 2011

By G van Niekerk, M Wortley

Screening is a critical part of many mineral processing circuits and keeping the aperture of the screening media within specification is vital to ensure efficient performance of the processing plant. Typical ways of estimating and/or measuring the wear of the screening media vary from very subjective methods (for instance, using the human eye) to manually checking the width of a few of the more severely worn gaps with a measurement device. Subjective methods can be flawed by human error while manual checking can distort the polyurethane or rubber materials that are used in several types of screening media. All of these errors can cause either premature replacement of panels or, even worse, poor screening of the material.Metso Process Technology and Innovation (PTI) and Metso Screening Media Solutions (SMS) have developed an innovative method that measures screen aperture using advanced image processing techniques. Field trials of the device are presented in this paper, which shows the comparison between the different methods of measuring apertures. Case studies demonstrate how the system, when combined with ScreenTrack, MetsoÆs industry leading maintenance system, and screening media, is providing users with the best possible screening solutions.

Sep 26, 2011

By E Nettleton, A Melkumyan

A probabilistic system for autonomous tracking of excavated ore in mining would be highly valuable. Currently, mine sites can have a vast array of sensors; examples include excavator load cells, bucket volume estimators and haul truck strut pressures which can be used to gather information about the excavated ore at different locations. These sensors can be utilised efficiently to increase the accuracy of mine estimates by introducing a real-time ore tracking system. This requires the introduction of a probabilistic ore tracking framework to consistently fuse information from any sensor input. Current modelling practices are manually intensive, deterministic (no recognition of the quality of data) and generally not a complete end-to-end system. This paper will provide a novel method for tracking extensive ore properties (mass, volume) and intensive ore properties such as chemical composition (eg Fe per cent, SiO2 per cent, Al2O3 per cent). This method is proved to ensure that correct correlations between states are conserved as material is transferred, combined and/or separated at different locations with any configuration of sensor inputs. This ensures that consistency of estimates is maintained and other benefits such as real-time probabilistic reconciliation can be implemented. A method for combining material of different mass and chemical compositions is also developed to enable accurate and efficient operation of the system. Representative experiments on a small-scale and a larger field trial which mimic the mining processes of excavation, haul and unload to ROM Stockpile will be used to verify the approach.

Sep 26, 2011

By J Berg, M Hahn, R Winkel

Imaging sensor systems designed for the mining and minerals industry incorporate different physical measurement technologies that can be used to detect the machine environment, such as laser, ultrasonic/sonar, video and RF-sensors. The underlying physical measurement principles cannot be employed in most applications near the face without regular maintenance due to the harsh environmental conditions like dust, dirt or moisture. However, that is the area, where imaging sensors are most valuable to the machine operator. In the domain of 2D radar sensors these have been either restricted by the angular coverage area or by limited adaptability of signal processing algorithms. The unique full digital solid state design will establish a generation step in mining applications.In order to fill this gap of terrain imaging sensors, the Institute for Mining and Metallurgy Machinery (IMR) of the RWTH Aachen University, in cooperation with the radar company indurad, developed an imaging sensor based on radar technology employing the digital beam forming (DBF) concept. Through multiple transmitting and receiving antennas the sensor delivers high accuracy range and angular information. Using adaptive digital signal processing with multiple imaging algorithms and target detection routines, data of the surrounding environment can be used to enable automation functions in mining.A comprehensive range of radar scans in different environments have been evaluated and benchmarked by using laser scanners. Accuracy measurements in range and angular direction have shown minimal differences. Experimental measurements have been performed at slopes, on conveyor belts and hot cement clinker silos. The experimental system uses a frequency modulated continuous wave (FMCW) front-end combined with a processing unit for signal analysis and target detection.

Sep 26, 2011

By S D. Dessureault

In large coal strip mines, there are many different processes at work: stripping, overburden removal using trucks/shovels, cast blasting, slot dozer pushing, more fragmentation and further truck/shovel benches. A key component to each of these steps is the use of machines. In modern operations, most machine costs are tracked within Enterprise systems that transact costs and repair work order management. Similarly, in contemporary mines most machines are tracked using technology that provides precise performance measures for each key process. Such technologies include productivity monitors on draglines, dozers and drills and fleet management systems for trucks and shovels. Linking these two sources of seemingly disparate data û costs and performance by machine û would allow for the precise calculation of unit costs by specific machine. This is similar to an activity-based costing system, except in this machine-based system, calculations are made at an even finer granular level. There are many uses for such an integrated information system in mine planning and performance management, both in real-time and through historical analysis.A real-time fleet costs (FLC) measurement process was developed as part of a larger blending-related, multi-million dollar project centred on a centralised remote control room developed at the University of ArizonaÆs Mining Information Systems and Operations Management Mine Control Center. The FLC component of the overall blending system would allow a mine controller to manage the equipment fleet by viewing its cost in real time and, when mixed with other revenue-related components, the exact profit. For example, a mine controller may show high utilisation while still maintaining a production surplus. This could occur, for example, if the equipment is driven long distances. This might appear in traditional performance metrics as æhigh production performanceÆ, but in reality may prove less profitable. This paper discusses the technical development of the FLC, its deployment and its impact on controller performance, and compares mine profitability before and after its application.

Sep 26, 2011

By Ivo Dreiseitl

Bottom sediment sampling by means of a box corer, as a direct method of deep seabed exploration for polymetallic nodules, is the most comprehensive way to obtain various types of geological information on points forming a network of stations. However, direct information is lacking from areas separating those stations and has to be approximated by interpolation or extrapolation, which leads to uncertainties. Sediment sampling with a gravity or piston corer can provide information on the superficial, up to 4 m thick, part of the sedimentary cover. As opposed to box coring, gravity or piston coring provides no information on nodule abundance. Information for unsampled areas can be gained by application of remote techniques of exploration. The techniques of choice in visualising nodule-rich fields, nodule-free areas, hard rock outcrops and other obstacles for seafloor mining involve using sonar and continuous deep sea (CDS) camera. Seafloor mapping can be rendered with two major sonar types: the side scan sonar used to obtain seafloor imagery and multibeam sonar for bathymetric data; the two types are usually combined and work with an intermediate frequency profiler. Photo-profiling by means of a deep sea towed camera provides a series of seafloor photographs, each supplemented with information on the area covered by each frame, including the nodule abundance. Remote techniques have been quite frequently applied by the Interoceanmetal Joint Organization (IOM) in nodule field exploration. Recently, geoacoustic techniques (side scan sonar + profiler) and CDS were applied for exploration of IOM?s first prime area H11 (5372 km2) on 6 transects. During the IOM-2009 cruise, a total of 295 km of geoacoustic transects and 344 km of photo-transects were surveyed. The data obtained were processed with the aid of earlier information produced by multibeam sonar surveys. The aim of this paper is to demonstrate the advantages of remote techniques, particularly with regard to application of the towfish device, for geological and geotechnical applications in exploration for polymetallic nodules.

Sep 14, 2011

By Sven Petersen

The southern Mid-Atlantic Ridge (SMAR) is one of the least-explored spreading centers in the world with respect to hydrothermal activity. In the past 8 years research focused mainly on the ridge segments close to Ascension Island, discovering hydrothermal fields as far south as 14°20?S (Haase et al., 2007; German et al., 2008; Tao et al., 2011). However, the vast majority of the SMAR ridge axis has never been thoroughly investigated for hydrothermal activity. In the spring of 2013 we used segment-scale AUV deployments as well as conventional CTD surveys to identify areas of hydrothermal activity along 16 mid-ocean ridge segments between 13° and 33°S, thereby covering close to 5 % of the global ridge axis in a single cruise. During cruise MSM25 of the German research vessel Maria S. Merian, we tested longrange (>100 km lines per dive) AUV deployments of GEOMAR?s AUV ABYSS (a Remus 6000 vehicle) along the spreading axis for their capability to locate active venting in a fast and efficient way on regional surveys. Since the missions covered long distances,

Sep 14, 2011

By Carsten Rühlemann

Between 2008 and 2013, the German Federal Institute for Geosciences and Natural Resources (BGR) carried out five exploration cruises to the German license area in the eastern Pacific Nodule Belt. The first two expeditions were mainly dedicated to multibeam mapping to obtain an overview of the seafloor topography and acoustic backscatter strength. These data were used to identify ten to fourteen prospective nodule fields for potential future mining, which together cover about 16% of the total license area of 75,000 km2. During the last three cruises, three of these potential mining areas were explored in detail. For this purpose, the BGR developed a suite of exploration methods to map nodule size distributions and nodule abundances. These methods include acoustic surveys with vessel-based multibeam systems as well as near-bottom video mapping complemented by in situ sampling. The economically most valuable field has a size of ca. 2,000 km2, of which 34% is covered by medium to large nodules (> 4 cm long axis of nodules) with an average abundance of 22.4 kg/m2 and 44% is covered by small nodules (< 4 cm) with a mean abundance of 17.5 kg/m2. The total mass of nodules in this field comprises more than 30 million tons wet weight and could sustain deep-sea mining for at least 10 years.

Sep 14, 2011

By Ye Jun

Until now, more than 200 sites of sea-floor hydrothermal activities are known on the midocean ridge but few of them are documented on South Mid-Atlantic Ridge (SMAR). This on some extent restricts our directly understanding of the hydrothermal mechanism on SMAR. In 2011, During the Cruise 22 of R/V Dayangyihao, a masive sulfide field was discovered in the center part of the first-order ridge segment that bounded between Cardno Fracture Zone and Saint Helena Fracture Zone, near 15.16°S/13.35°W, which give us an opportunity to recognize the mechanism of hydrothermal systems on SMAR directly. This newly discovered hydrothermal field, which informally called SMAR 15°S, is situated in a SE-NW offset depression of an axial neovolcanic ridge, with water depth of about 2793m. In the water column around this field, strong turbidity and temperature anomalies were detected. Ten meters high venting chimneys, hydrothermal mounds as well as a large amount of hydrothermal vent shrimps were also observed and recorded by ROV and Deep tow camera on the ocean floor. The samples we collected from this field include sulfide chimney fragments, massive sulfides and basaltic rocks. In this study, we are focusing on the mineralogy and sulfur isotope geochemistry of sulfide samples for preliminarily uncovering the geological characteristics and hydrothermal processes of this hydrothermal field.

Sep 14, 2011

By Robert G. Ditchburn

Msjflkcisokng During the past sixteen years, GNS Science has developed a number of radiometric methods for dating massive sulfide samples from submarine arc and back-arc volcanoes. This has improved our understanding of the formation of these hydrothermal systems, the frequency and duration of their activity, and assisted marine mineral exploration companies in assessing the economic potential of these deposits that can be rich in copper, zinc and gold (up to 90 ppm Au in some chimneys). The radiometric dating methods that we employ utilize radium isotopes that are extracted from the host rocks without their parent thorium isotopes, by ascending, hot, acidic hydrothermal fluids. When the hydrothermal fluid discharges into cold seawater, chemically similar barium and radium co-precipitate as barite, a common mineral contained within the black smoker chimneys that grow upwards from the seafloor. From the onset of mineralization, the isotopes 228Ra and 226Ra (with half-lives of 5.75 yrs and 1600 yrs, respectively) begin to decay, so that the 228Ra/226Ra and 226Ra/Ba values decrease with time. Thus, samples can be dated using these ratios once the initial values, i.e., the values at the onset of mineralization, have been established, ideally from active chimneys. 228Ra decays to 228Th (half-life 1.91 yrs) and from the 228Th/228Ra values, the growth-rate of individual chimneys can be estimated. The isotope 210Pb (half-life 22.3 yrs) is occasionally used for dating the mineralization, although several samples from a chimney are needed to determine an age from the 210Pb grown from 226Ra, using the isochron technique.

Sep 14, 2011

By Brian Parsons

In June of 2013, Odyssey Marine Exploration conducted sea trials of the newly installed SeaBat 7150 full-ocean multibeam sonar system on the M/V Dorado Discovery over the Fierberling Guyot in the northeast Pacific Ocean. Fieberling Guyot is located approximately 990 kilometers west of San Diego, California and it rises abruptly from the deep ocean floor (~4300 meters) to a minimum depth of approximately 438 meters. The SeaBat 7150 operates at two distinct frequencies, 24 kHz (shallow to 3000 meters) and 12 kHz (full-ocean depth greater than 2000 meters). Onboard operational data are managed by a suite of state-of-the-art software. A SQL Server 2012 database provides the highly specific management and visualization tools as ArcGIS 10, Micromine 2013, and Geobank 2013 database. The 7150 is mounted on a gondola one meter lower than on previous system and greatly reduces bubble interference. The 7150 generates 880 beams

Sep 14, 2011