Search Documents

By S. W. Jacobsz, Y. Narainsamy

Several high-impact tailings dam failures around the world in recent years have placed a renewed focus on the stability of tailings dams and pointed to potential shortcomings in traditional drained design and safety evaluation procedures. A need to consider undrained shear strength in the design of tailings facilities has become apparent. However, there are specific requirements that need to be met before undrained shearing occurs. In South Africa, the last major failure was likely the Merriespruit disaster in 1994, which leads to the question of whether the conditions required for undrained shearing are readily applicable to South African tailings dams. This paper describes research at the University of Pretoria that has recently commenced to further investigate the conditions required for undrained failure to occur. The research includes laboratory and field testing to replicate these conditions in the laboratory and relate them to those found in an active tailings dam.

Jun 2, 2022

By Y. Chen, D. He, J. H. Potgieter

The demand for rare earth elements (REE) is increasing rapidly owing to the emergence of new cleanenergy and defense-related technologies. China’s dominates the world production of REE. In this paper, we review the ore characteristics of four typical REE deposits in China and discuss the conventional beneficiation techniques, including gravity and magnetic separation, flotation, and hydrometallurgical processing of these ores. Some of the latest laboratory findings and industrial applications, as well as the development of novel reagents for REE processing in China are presented. The REE beneficiation technologies developed and employed in China could constitute a good case study for many other emerging projects being explored around the world.

Mar 2, 2022

By O. S. Tinkler, K. C. Sole

Approximately 20% of current world copper cathode output is produced using a hydrometallurgical process route, generally referred to as the leach–solvent extraction–electrowinning flowsheet. Since its commercialization in the late 1960s, steady improvements in the performance and efficiency of the solvent-extraction reagents and equipment, combined with significant developments in leaching and electrowinning, have made an ever-widening range of ore types amenable to this technology. Following successful implementation on all major continents, a large proportion of growth in recent years derives from the re-emergence of copper solvent extraction in the Central African Copperbelt. This review provides a brief history of the development and evolution of copper solvent extractants and mixer-settler contactors, and the significance of the Copperbelt region in achieving commercialization and acceptance of the technology. The opportunities and challenges presented by the abundant high-grade oxide ores of the Copperbelt are contrasted with the processing of solutions derived from the low-grade mixed oxide–sulfide ore bodies that are prevalent in other geological regions. The current status of hydrometallurgical copper production in the African Copperbelt, within a global context, and a medium-term outlook for the technology are discussed.

Jul 5, 2023

By A. Yañez, J. Yianatos, P. Vallejos, R. Grau

This paper presents an evaluation of the effect that the mineral feed characteristics have on the metallurgical performance in an industrial flotation circuit. Variables such as mineralogy, particle size, solid content and mineral grade were evaluated using a flotation simulator, which was previously developed and validated using industrial data from different concentrators. The metallurgical performance was evaluated not only from the point of view of recovery and concentrate grade, but also by analysing less common variables such as froth carrying rate, froth and collection recovery, water recovery and others. Then, new data from the industrial flotation circuit was used to complement this analysis and evaluate the simulator predictability by comparing simulated and measured data. The results allowed for the quantification of the significant effect that mineralogy and particle size have on recovery and/or concentrate grade. Therefore, the grinding stage and flotation cells design become important factors to avoid a significant loss in metallurgical performance when feed mineral characteristics change in an industrial operation. In this context, an evaluation of cells design and arrangement should be performed to compensate the loss in recovery and concentrate grade in order to find an optimal operating condition. Keywords: Flotation circuits, mineral feed, modelling, particle size, liberation, mineralogy

Jan 1, 2020

By R. Xie, J. Liu, Y. Zhu, Y. Li

The surface chemistry features of spodumene and feldspar were investigated by XPS, zeta potential measurements and first principle calculation. The results indicated that the zeta potential of feldspar was more negative than that of spodumene. The most possible active sites on spodumene and feldspar were Al sites and O sites. However, the number, type and activity of these sites were different. The flotation behaviour of five new collectors (cationic collector decyloxy-propyl-amine (DPA), anionic collector α-bromododecanoic acid (α-BDDA), N-Dodecyliminodiacetic acid (DIDA), and mixed collector DPA-BDDA and DPA-DIDA) were investigated through micro flotation experiments and analysed based on the surface chemistry features of spodumene and feldspar. The results revealed that cationic collector DPA processed a high collecting capacity for both minerals. The highest recovery of spodumene and feldspar was 91.60% and 95.30% respectively. The anionic collector α-BDDA and DIDA showed poor flotation behaviour. The recovery of spodumene and feldspar when using α-BDDA as collector was less than 5% and when using DIDA as collector was less than 1%. For the mixed collector DPA-BDDA, the recovery of spodumene and feldspar was 82.27% and 32.61% respectively at pH 4.5. For mixed collector DPA-DIDA, the recovery of spodumene and feldspar was 35.69% and 9.08% respectively at pH 5.8. The results indicated that the cationic collector had high collecting capacity, but low selectivity. The collecting capacity of anionic collector for spodumene and feldspar was poor at the same concentration compared with cationic collectors. Mixed collector DPA-BDDA showed good flotation behaviour and selectivity on spodumene. The ―COO- in mixed collector might improve the selectivity. However, too much ―COO- would increase the electrostatic repulsion, which might cause low selectively and poor recovery. Keywords: Spodumene, feldspar, surface chemistry features, mixed collector, flotation

Jan 1, 2020

By Z. Xingrong, L. Jiongtian, Z. Yangge, C. Yijun, L. Liang, H. Yangjia, H. Long

In this present research, a polymer poly (acrylamide-allyl thiourea) (PAM-ATU) with an average molecular weight 1.66 × 104 g/mol was used as a depressant in the flotation separation of molybdenite from chalcopyrite and galena. The micro-flotation experimental results showed that PAM-ATU exhibited a strong influence on the flotation performances of chalcopyrite and galena, but had little effect on that of molybdenite within a pH range of 9.0~11.0 at a dosage of lower than 6.0 mg/L. The flotation results of bulk flotation tests illustrated that PAM-ATU showed a stronger depression power than sodium thioglycollate (TGA) and sodium sulfide. The comparison results demonstrated that PAM-ATU and TGA had an advantage in reducing galena and chalcopyrite impurities in the molybdenum concentrate, respectively. With the combination of PAM-ATU and TGA, a good flotation separation result was achieved through one roughing operation where the recovery of molybdenum in the froth product was over 75.0%, while the recovery of copper and lead were less than 22.0% and 15.0%, individually. It was also found that a molybdenum concentrate with Mo grade of 46.32%, Cu grade of 0.16%, Pb grade of 0.06%, was obtained by a closed circuit flotation flowsheet of one roughing-three scavenging-five cleaning with a Mo recovery of 77.67%. Keywords: PAM-ATU, molybdenite, depressant, chalcopyrite, galena

Jan 1, 2020

By P. J. Le Roux

"Mining companies are under constant pressure to reduce their cost structures to sustain profitability. In mines using an open stoping mining method, dilution often has a significant effect on profitability. At Target mine in South Africa, dilution in some open stopes was found to exceed 10%, falls of ground being a major contributor to the problem. This dilution can reduce the recovered grade from 5.5 to 4.5 g/t, resulting in a potential loss of about R21 million per month based on a gold price of R240 000 per kilogram. In addition to the cost of dilution, the cost of damage to, or loss of, trackless equipment as a direct result of falls of ground in open stopes is significant. Other associated costs include transport, hoisting, secondary blasting, milling, and plant treatment of waste material. Stable open stopes are essential in order to reduce falls of ground. Back-analyses of stope instability at Target mine have indicated that conventional rock mass failure criteria are unsuitable for stope design. An alternative strain-based criterion has been developed, and has proved to be very successful, allowing the stability of open stopes to be calculated reliably. Since its implementation in 2010, dilution and equipment damage have decreased markedly, creating value for the mine. IntroductionMines that make use of an open stoping mining method aim to extract only the ore, leaving the waste behind. Many orebodies have well-defined boundaries between the ore and the waste rock. In such orebodies, the introduction of waste into the ore due to overbreak dilutes the grade. In massive disseminated orebodies, dilution is less problematic. Overblasting, falls of ground, or other causes may result in dilution. A study undertaken in Canada (Pakalnis et al., 1995) found that approximately 51% of all underground metal mines utilized open stoping mining methods and, from surveys conducted at these mines, open stoping operations experienced dilutions of up to 20%, and occasionally more. Dilution of this magnitude has a significant economic impact on any mining venture. Research carried out in Australia by Capes (2009) came to the same conclusion."

Jan 1, 2017

By P. Gonzalez

Minera La Parreña, S.A. de C.V., Milpillas operation is an underground mine designed to process 8 352 tons of ore per day for the recovery of copper through crushing, heap leaching, solvent extraction and electrowinning facilities. In order to optimize copper recovery during leaching, the ore is crushed to 13 mm and subjected to primary sampling every 15 minutes for metallurgical accounting. At the same particle size, primary sample goes to secondary sampling and then, the secondary sample is crushed to 2 mm before tertiary sampling. The tertiary sample is sent to laboratory where a final sample is obtained and pulverized to 106 microns for copper assay. This paper shows a review of sampling strategy focused on minimizing or eliminating sampling errors and an experiment designed to determine copper head grade profile as well as the best sampling frequency.

Jan 1, 2009

By A. McQuillan, N. Bar

Geotechnical models developed during the planning stages of open pit mines are three-dimensional so as to capture the spatial variation in lithological, structural, hydrogeological, and geomechanical conditions. Geological models that describe the lithological and structural (faulting and folding) characteristics of a deposit are always 3D. Likewise, boreholes and piezometers used to develop geomechanical properties and groundwater models are drilled at spatial offsets across the deposit to understand the lateral and vertical characteristics. Yet when geotechnical analysis is completed, often the three-dimensional geological, hydrogeological, and structural models as well as geometrically complex 3D mine designs for optimizing economic mineral recovery and overburden removal are simplified to two-dimensional sections. In this paper we demonstrate that this simplification can lead to the wrong failure mechanism being identified, analysed, and/or a conservative factor of safety being calculated and hence an over-estimation of slope stability. Through case studies we show how three-dimensional analysis methods are more suited to rock slopes, particularly those with anisotropic material strength, when singularities such as geological faults are present, and nonlinear slope geometry. When the same slopes are analysed in two dimensions, the failure mechanism calculated is often fundamentally incorrect. The case studies further reveal that the factor of safety calculated in three dimensions is not always higher than the two-dimensional factor of safety.

Feb 6, 2023

By L. Ntlhakana, K. Khoza-Shangase, G. Nelson, I. Maposa

Occupational noise-induced hearing loss (ONIHL) is associated with various risk factors in South African mines. We investigated the association between standard threshold shifts (STSs) and exposure to noise and platinum mine dust (PMD), using demographic data, five years of annual audiometry screening results, and noise and dust exposure data. Miners’ age, sex, percentage hearing loss (PLH), and dust and noise exposure data were gathered, and a linear mixed effects regression model used to predict STS. Average occupational exposure levels to noise and dust were calculated from recorded measurements. A total of 12 692 records were analysed. Most miners were male (89.6%) and more than 50% were younger than 41 years. More than 70% were exposed to > 85 dBA noise and 58% were exposed to 1.5–2.99 mg/m3 PMD. Changes in hearing levels ranged from 8.3 dBHL at baseline (2014/2015) to 10 dBHL in 2016 (STS), with no changes thereafter. The model showed that the estimated effect on STS for males was 27% and 21% higher than for females, for the left and right ear, respectively The estimated effect of age, PLH, noise exposure and years of exposure on STS was < 10% for each variable. There was no statistically significant association between PMD and STS. The combined effects of age, sex, years of exposure to noise, and noise exposure levels, and strength of associations can be used to predict STS for this group of miners. Our findings may be used to measure the efficacy of the mine’s hearing conservation programme.

Aug 1, 2021

By X. Ndlovu

Underground observations in a coal mine have indicated failure oft he immediate roof above the bords. Careful observations and photographic recording of occurrences of these roof failures, termed ?roof guttering?, have been made. The occurrence of guttering results directly in roof instability, and indirectly, due to the interaction of the guttering with natural planes of weaknesses in the immediate roof rock. In this paper, relevant observations of roof guttering, possible mechanisms of formation of roof gutters, laboratory testing for determination of material strength and deformation properties, and numerical analyses to explain underground observations are described. Mapping of roof failures showed that these took place mainly towards the centre of the roadways. Roof guttering was observed often to occur violently and with little warning. Occurrence of roof guttering has a negative impact on production?some panels are abandoned, production times are increased and safety of workers is compromised. The numerical stress analyses showed that, although none of the constitutive behaviour criteria used predicted the observed underground failures correctly, the extension strain criterion gave the best agreement.

Jan 1, 2007

By K. Isaac, S. O. Bada

This research focuses on the co-firing of discard coal with refuse-derived fuel (RDF) to utilize this abundant resource in South Africa for energy generation and reduce the volume of waste disposed of at landfills. The potential of a coal with a high ash content (> 40%), which is a grade used in some power stations in South Africa, and its combustion compatibility with two different RDFs in terms of emission reductions has been established. Gaseous emissions and ash residues from the combustion and cocombustion of the coal, two different RDFs and coal/RDF blends of different proportion were analysed. One of the RDF samples contained mostly paper (PB) and the other mostly plastic (PL). Co-combustion ash from the discard coal and RDFs showed a decrease in chloride and alkali metal contents as the coal ratio in the blend increased. The slagging propensity of the co-fired blends was found to be very low, while the propensity for fouling decreased from high to medium for all the blends with < 75% RDF. Cocombustion of RDF with coal showed a decrease in SO2 emissions from 387 ppm (discard coal) to 50 ppm for the sample containing 25% coal discard plus 75% PL. A 15% PL to 85% coal blend also reduced NOx emissions from 145 ppm (100% PL) to 88 ppm. The lowest CO2 emission observed was 6000 ppm for the blend of 85% discard coal plus 15% PB. It was established that the most favourable fuel blend that can produce the lowest sulphur emissions if used for power generation is the 25% coal discard plus 75% PL sample.

Aug 2, 2022

Metso Minerals Cisa have developed VisiofrothTM - a system which calculates parameters of the froth surface using vision from remote cameras mounted about individual flotation cells. The parameters calculated from this system can be utilised within process control strategies to achieve stabilisation or optimisation of the flotation process. As part of a JKMRC AMIRA P9 field campaign, froth vision was collected from various cells operated at different air rates and froth depths within Rio Tinto&apos;s Northparkes copper roughing and scavenging circuit. At the same time, samples were collected to determine the cell metallurgical performance at each operating condition. These froth images have been processed using the VisiofrothTM system with froth velocity, bubble size, froth collapse rate and various froth colour parameters calculated. These parameters have been compared to the measured concentrate flow and concentrate grade produced at each condition. This analysis shows that the concentrate tonnage flow is strongly correlated with froth velocity whereas the concentrate grade is best correlated using froth velocity in combination with a bubble stability parameter which affects the bubble size distribution on the top surface of the froth. This paper summarises the results from this analysis as well as an overview of the VisiofrothTM system and the parameters that are measured by the software.

Jan 1, 2007

By I. Zemanek, A. Rozsypal

Mrazovka tunnel was excavated under one building whose subsurface space was not possible to reinforce against additional settlements induced by digging. In the presented contribution the risk analysis, warning levels and results of deformation measurements are described. Lespace souterrain d‘un des bâtiments sous lequel le tunnel Mrazovka a été creusé n‘a pu être renforcée contre l‘affaissement supplémentaire provoqué par l‘emplacement du tunnel. Le matériel présenté décrit l‘analyse des risques, le niveau des limites et les résultats des mesures des déformations. Der Untergrund eines der Gebäude, unter dem der Tunnel "Mrazovka" ausgeschachtet wurde, war es nicht möglich, gegen nachträgliches Setzen zu versteifen, hervorgerufen durch die Untertunnelung. Im vorliegenden Beitrag sind Gefahrenanalyse, Grenzwerte und Ergebnisse von Deformationsmessungen beschrieben.

Jan 1, 2003

By C. Montoya

This paper shows, through a case study, the impact of multivariate grade modelling upon mine design and mine planning. A deposit explored by drill holes is considered, in which the grades of five elements (copper, silver, molybdenum, arsenic, and antimony) are of interest. Forty alternative models of the deposit are constructed by fitting the joint correlation structure of the grade variables and using conditional cosimulation. In addition, a reference model, obtained by averaging the alternative models, is also considered. The study shows that the resulting mine design (final pit characteristics and production schedules) is sensitive to the grade model under consideration, and that the design based on the reference model may not be optimal when compared to the alternative models based on cosimulation. However, when assuming a given long-term plan and extraction sequence, the grades and net present value (NPV) calculated on the reference model are unbiased with respect to those calculated on the alternative models with the same extraction sequence. The latter allow assessing the possible dispersion of the actual grades and NPV around their expected values, and are useful for the planner in order to determine the probability of meeting given production targets and of exceeding or falling short of given threshold grades. Additionally, unlike cosimulation, the separate simulation of each grade variable leads to unrealistic resource models and to biased results in mine design and mine planning. This approach should therefore be avoided, unless the grade variables are spatially uncorrelated.

Jan 1, 2012

By S. Kauchali

"Underground coal gasification (UCG) is recognized as an efficient mining technique capable of chemically converting the coal from deep coal seams into synthesis gas. Depending on the main constituents of the synthesis gas, chemicals, electricity, or heat can be produced at the surface. This paper provides a high-level graphical method to assist practitioners in developing preliminary gasification processes and experimental programmes prior to detailed designs or field trials. The graphical method identifies theoretical limits of operation for sensible gasification within a thermally balanced region, based primarily on the basic coal chemistry. The analyses of the theoretical outputs are compared to actual field trials from sites in the USA and Australia, with very favourable results. A South African coal is studied to determine the possible synthesis gas outputs achievable using various UCG techniques: controlled retractable injection point (CRIP) and linked vertical wells (LVW). For CRIP techniques, an important result suggests that pyrolysis, and subsequent char production, are important intermediate phenomena allowing for increased thermal efficiencies of UCG. The conclusion is that South African coals need to be studied for pyrolysis-char behaviour as part of any future UCG programme. The results also suggest that UCG with CRIP would be a preferred technology choice for Bosjesspruit coal where pyrolysis dynamics are important. Lastly, the use of CO2 as oxidant in the gasification process is shown to produce syngas with significant higher heating value. Introduction and literature surveyCoal is a commonly utilized fossil fuel, providing over 40% of global electricity demand and about 90% of South Africa’s primary energy needs. However, less than 20% of the known world resources are suitable for possible extraction using conventional surface and underground mining techniques (Andrianopolous,, Korre, and Durucan, 2015). Underground coal gasification (UCG) has the potential to recover the energy stored in coal in an environmentally responsible manner by exploiting seams that are deemed unmineable by traditional methods. The UCG process, if successfully developed, can increase coal reserves substantially. For example, in the Limpopo region of South Africa alone, the estimated potential for UCG gas, based on existing geological records, is over 400 trillion cubic feet (TCF) natural gas equivalent – this is about a hundred times more gas than the existing 4-TCF Pande-Temane natural gas field reserve in Mozambique (de Pontes, Mocumbi, and Sangweni, 2014)."

Oct 1, 2018

Spiral concentrators are compact, cost-effective and generally efficient gravity concentration separators for a wide range of applications (for example: coal, beach sands, iron ore, chromite and tantalite). Large mineral processing plants consist of thousands of spiral concentrators resulting in large plant footprints (capital intensive) and the adjustment of splitters is time consuming, impractical and is in many cases neglected?high capacity (HC) spiral concentrators aim to address these shortcomings. As a result Exxaro Namakwa Sands is currently investigating high capacity spiral technology for the spiral circuit upgrade at the primary concentrator plants (PCPs). This article summarizes the rougher spiral performance evaluation that was conducted on different types of spiral concentrators (the traditional MG4 spiral concentrator and the high capacity (HC) spiral concentrator) under different feed conditions. In addition, the effect of slimes on the spiral concentrator performance was also investigated. Slimes rheology was linked to the poor concentrator performance at the higher slimes concentrations. The test campaign shows a sacrifice in recovery under design-feed conditions can be expected when using high capacity spiral concentrators in the rougher stage when compared to traditional spiral technology currently in use. Both spiral concentrators show a detrimental impact of slimes on the performance, but the high capacity spiral concentrator is more sensitive to the higher slimes conditions.

Jan 1, 2010

By J. E. Xu, J. Y. Zhang, M. Q. Sheng, K. Wan

"SynopsisThe densities of selected CaO-Al2O3-MgO slag systems were measured at 1823K by the Archimedean method. Thirteen different slag compositions were chosen based on different levels of the MgO content and the mass ratio of CaO/Al2O3. The results indicated that the density of the slag decreases with increasing MgO content (from 0–3.78 mass%), but increases with further increases in MgO content up to 11.33 mass%. At a fixed MgO content of 5.5%, the trend of density change with CaO/Al2O3 in the slag is similar to that for changes in MgO content. On the basis of the regular solution approximation rules of excess molar quantities, an attempt was made to estimate the molar volume of the slags investigated. The application of the molar volume model confirmed that the present expanded approximation rules are applicable to predict the molar volumes of the melts discussed.IntroductionDensity is one of the most important fundamental properties of slags, which influences metallurgicaly phenomena in many ways. Density is sensitive to the nature of the chemical bonds and the species in the melt, and is of great importance not only for theoretical research on the structural properties of molten slags, but also for industrial applications. Density is also required to estimate other key properties used to assess the behaviour of high-temperature molten oxides, e.g. viscosity, surface tension, and thermal conductivity.Liquid calcium aluminate slags that contain magnesia are the basis of most ladle slags for secondary refining processes. The CaO-Al2O3 binary phase diagram (Hallstedt , 1990) shows the eutectic composition (50 wt% CaO and wt% Al2O3) has the lowest melting temperature. The main emphasis of this study is on the 12CaO·7Al2O3 refining slag with MgO additions. The CaO-Al2O3-MgO ternary system is of considerable importance in industrial metallurgical processes, and is fundamental to the understanding of metallurgical slags, ceramic materials, and geological phenomena."

Jan 1, 2015

By J. F. Lupo

One of the earliest records of metal recovery by solution leaching is described by Agricola as &apos;juice of rock&apos; in the 1550s. In these early years of hydrometallurgy and civil engineering, few controls and systems were employed to enhance recovery and protect the environment. Since that time, leaching (dump or heap) operations have made significant strides in increasing metal and solution recovery while protecting the environment, all goals of a sustainable operation. It is now recognized that the design and operation of modern heap leach facilities requires contributions from many fields of study, including hydrometallurgy, civil engineering, geotechnical engineering, unsaturated-flow hydrology, mine planning, geosynthetics engineering, geochemistry, process engineering, mechanical engineering, and electrical engineering. While advancements in these fields have resulted in more sustainable heap leach operations, challenges in the industry still exist. A number of heap leach operations exhibit poor or lower-than-predicted metal recovery, loss of solution flow and control within the ore heap, loss of ore heap stability under leach, failure of liner and/or solution recovery systems, and overtopping of process water ponds. A number of these issues may be the result of several compounding conditions. For example, poor metal recovery may be due to an inadequate scale-up assumption (scaling laboratory tests to field-size heaps), lack of control of the in-heap geochemical environment, changes in ore mineralogy from the original design, changes in mechanical and hydraulic properties of the ore from the original design, ore handling and pre-treatment, inadequate solution management system, and inadequate solution application. The purpose of this paper is to present and discuss issues that may affect the sustainability of a heap leach operation. Since sustainability encompasses a broad range of topics and issues, the focus of this paper will be on issues affecting metal and solution recovery, solution flow and containment, and stability of the ore heap.

Jan 1, 2012

By M. Cuba

Geostatistical modelling aims at providing unbiased estimates of the grades of elements of economic interest in mining operations, and assessing the associated uncertainty in these resources and reserves. Conventional practice consists of using the data as error-free values and performing the typical steps of data analysis - domaining, semivariogram analysis, and estimation/simulation. However, in many mature deposits, information comes from different drilling campaigns that were sometimes completed decades ago, when little or no quality assurance and quality control (QA/AC) procedures were available. Although this legacy data may have significant sampling errors, it provides valuable information and should be combined with more recent data that has been subject to strict QA/QC procedures. In this paper we show that ignoring the errors associated with sample data considerably underestimates the uncertainty (and consequently the economic risk) associated with a mining project. We also provide a methodology to combine data with different sampling errors, thus preserving the relevant global and local statistics. The method consists of constructing consistent simulated sets of values at the sample locations, in order to reproduce the error of each drilling campaign and the spatial correlation of the grades. It is based on a Gibbs sampler, where at every sample location, the actual sample value (with error) is removed and a conditional distribution is calculated from simulated values at nearby sample locations. A value is drawn from that distribution and kept only if it satisfies some statistical requirements-specifically, the global relative error and local means and variances must be reproduced. All sample locations are visited and simulated sample values are generated iteratively, until the required statistics are satisfactorily attained over all sample locations. This generates one realization of possible sample values, respecting the fact that the actual samples are known to carry an error given by the global relative error. Multiple realizations of simulated sample values can be obtained by repeating the procedure. At the end of this procedure, at every sample location a set of simulated sample values is available that accounts for the imprecision of the information. Furthermore, within each realization, the simulated sample values are consistent with each other, reproducing the spatial continuity and local statistics. These simulated sets of sample values can then be used as input to conventional simulation on a full grid to assess the uncertainty in the final resources over large volumes. The methodology is presented and demonstrated using a synthetic data-set for clarity.

Jan 1, 2012