Search Documents

By R. Trammell, O. Eterigho-Ikelegbe, S. O. Bada

In response to the enormous amounts of coal discard generated during coal mining and preparation, the development of an efficient and sustainable strategic use of this resource is essential. Furthermore, the rising urban population over the next decades is confronted with the depletion of quality raw materials for building components. To this end, this study reports new information on the morphology, water absorption, and flexural strength properties of ceramic composites produced from three different coal discards and polysiloxane pre-ceramic polymer (PCP) resin. In addition, test results relating to the continuous operating temperature, chemical resistance, and efflorescence potential of the composites are presented. The results show that the water absorption and flexural strength of the coal composites, up to 1.94% and 36.46 MPa respectively, exceed the requirements for ceramic and clay roof tiles. The continuous operating temperature of composites is found to be more thermally stable than conventional roofing tiles (concrete and ceramic) between ambient temperature and 600°C. In addition, the excellent chemical resistance of the composites (94.43%−99.98%) compared to conventional roofing tiles (67.82%−99.97%) eliminates the need for additional external coatings. The interesting results documented so far suggest that this technique could be used to produce low-temperature application building products such as bricks, panels, roofing tiles, etc. This new recycling technique offers an excellent opportunity to eliminate enormous volumes of coal discard and to advance the circular economy in the coal industry.

Aug 2, 2022

By S. Schutte, J. Pelders, F. Magweregwede, S. Ngobese, W. Botha

A people-centric, systematic approach that involves different stakeholders is considered necessary to successfully address changes associated with modernizing the South African mining industry. In this paper we discuss the inclusion of various South African mining industry stakeholders, particularly the employees, in equipment design and development processes. Data was acquired through a literature review and inputs from South African mining industry stakeholders through focus group discussions, interviews, and an online survey. It was suggested that some of the areas that could be improved by the inclusion of employees in the process are design quality, ergonomics, equipment acceptance by the workforce, safety, efficiency, skills, insight into context of use, and early design iterations and identification of problems. Among the issues considered to be barriers to inclusion were lack of buy-in, lack of resources, difficulty in manufacturers accessing the mines, difficulties in involving too many people, intellectual property concerns, and unavailability of mine employees due to their key responsibilities. Identified critical factors in the equipment design and development process include iterative design, stakeholder participation throughout the process, needs analysis through consultation, and change management. A guideline was developed for worker inclusion in equipment design and development that could be used by the South African mining industry. The study identified several economic and social benefits of including stakeholders in the early stages of design and development. It is recommended that the South African mining industry considers using the developed guideline.

Apr 5, 2023

By N. Nišić, A. Cvjetić, V. Jovanović, D. Knežević, D. Nišić

The risks associated with operating water reservoirs and tailings storage facilities (TSFs) are different because of their different purposes, methods of construction and operation, and characteristics of the materials impounded and their flow behaviour. Regardless of the differences, these two types of structures are often put in the same category when it comes to risk assessment and the application of relevant methodologies, which may result in unrealistic outcomes. In this paper we investigate whether it is possible to apply the risk ranking methodology designed for water reservoirs to TDFs and overview the key differences between these two types of structures. We also provide a comparative analysis of the results obtained by conducting a risk assessment of the operation of a TSF, applying the method recommended by the International Commission on Large Dams (ICOLD), and analyse the results using the failure mode and effects analysis (FMEA) technique. It is concluded that a more realistic evaluation of risk levels can be obtained by adopting a systematic approach to risk assessment in accordance with the traditional risk definition on which the FMEA technique is based.

Nov 30, 2022

By M. M. Bwalya, F. K. Mulenga

"SynopsisThe application of the attainable region (AR) technique to the analysis of ball milling is currently limited to batch data. This paper introduces the use of the technique to continuous milling.To this end, an industrial open milling circuit processing a platinum ore was surveyed. Samples were collected and later characterized by means of laboratory batch testing. On site, several milling parameters were varied systematically so as to collect data for modelling purposes. These paramters included ball filling, slurry concentration, and feed flow rate. After data analysis, a simulation model of the open milling circuit was developed under MODSIM®, a modular simulator for mineral processing operations. The mill was then simulated and the data generated was analysed within the AR framework.Initial findings reveal an opportunity to gain valuable insight by studying milling using the AR technique. From an exploratory perspective and inasmuch as this study is concerned, feed flow rate, ball size, and ball filling were identified as being pivotal for the optimization of open ballmilling circuits. Mill speed, on the other hand, had only a limited effect on the production of particles in the size range -75 +10 µm.IntroductionGlasser and Hildebrandt (1997) propounded the attainable region (AR) as a technique for the analysis of chemical engineering reactor systems. Results have since been produced and tested on both the laboratory and pilot scales. The method results in a graphical description of chemical reactions by considering the fundamental processes taking place in the system, rather than the equipment. From the plotted graphs, the process and the reactors can be synthesized optimally into a flow sheet.The use of the AR method still has a long way to go as far as mineral processing is concerned. For instance, several articles have reported the application of the AR technique to ball milling: Khumalo et al. (2006, 2007, 2008); Khumalo (2007); Metzger et al. (2009, 2012); Metzger (2011); Katubilwa et al. (2011); Hlabangana et al. (2012). Their main shortcoming has been the exclusive use of laboratory batch grinding data.To address this deficiency, Mulenga and Chimwani (2013) proposed a way by which the technique could be extended to continuous milling. In effect, the batch milling characteristics of a platinum-bearing ore (Chimwani et al., 2013) were used and scaled up to an open milling circuit. Then, with simplifying assumptions, an attempt was made to optimize the residence time of particles inside the mill. Later, Chimwani et al. (2014a) presented some optimization examples involving various milling parameters. The sequence of published articles then paved the way for the study of industrial milling systems with the AR methodology."

Jan 1, 2015

By T. Szendrei, S. Tose

The fracturing and movement of rock that occurs in the vicinity of a stemmed borehole charge in open pit mining operations are described by examining the effects of the emitted stress waves – shock and elastic – and the expansion of high-pressure detonation product gases. Three principal modes of momentum transfer to fractured rock are identified, all linked to gas expansion work. This work can be delivered radially (burden), axially (stemming), and in the collar zone (cratering). Flyrock is generated under unusual combinations of blast parameters and rock properties. Although infrequent and seldom predictable, the generation of flyrock can nonetheless be interpreted and modelled in terms of the principal mechanisms of rock projection described in this study. The physical processes underlying these principal mechanisms are identified and will permit the development of predictive models for flyrock velocities.

Dec 13, 2023

By S. A. Villalobos, F. A. Villalobos, and L. E. Aguilera

Rockburst events have been a serious problem for many years in many mines worldwide, and in particular at El Teniente mine in Chile. El Teniente is the largest copper mine in the world, located in the Andes Cordillera where high stress levels are present due to intensing mining activity in addition to complex geology. Consequently, the study and management of the rockburst threat are necessary. In this work, the case of the Diablo Regimiento (DR) mine within El Teniente is studied. The energy capacity of dynamic support systems is determined for different tunnel geometries based on two kinetic methodologies, using data from DR. Initially, rockburst potential is determined by means of a stress analysis around different tunnel geometries through the boundary elements method. In the first methodology a yielding zone (YZ) is estimated for each excavation geometry using the finite element method FEM. The second methodology involves the definition and determination of a critical strain energy (SE) around each excavation geometry using a FEM numerical analysis. In both cases, peak particle velocity PPV is estimated by a scaling law, which is subsequently adjusted due to tunnel amplification effects. According to the results, and knowing the working energy capacity applied in DR mine, it was found that the values of energy capacity for the rock dynamic supports were better estimated by the YZ-PPV approach than by the SE approach.

Sep 1, 2022

By F. Hörtkorn, Y. Narainsamy

Rather quietly, a new generation of European norms dealing with earthworks has been developed over the past few years, in particular a part of EN 16907 pertaining to extractive waste. Whereas the majority of the parts deal with routine geotechnical aspects, one part deals specifically with tailings dams. The standard is currently in draft form but if this becomes a European Standard, members of the European Union are bound to comply with the regulations. As such, the standard has the capacity to influence mining in Europe and potentially in South Africa. In this paper we review the standard and compare it to the South African standard on mine residues and the recently published Global Industry Standard on Tailings Management (GISTM). We conclude found that although the European norm makes great strides towards increasing tailings dam safety, there is room for improvement. The current South African standard has served well since its inception, but can be improved by using concepts from the European and global tailings standards.

Jun 2, 2022

By J. Arzuúa, U. Castro-Filgueira, L. R. Alejano, F. Malan

"Pillar design for underground mining is typically done using empirical formulae or numerical modelling. Practical experience and recent literature, nevertheless, illustrate the shortcomings of these design approaches. Ongoing monitoring of pillars is therefore recommended to minimize the risk associated with these designs.In a mine where a large number of pillars are cut, different pillar strengths can be expected owing to variations in rock mass strength from area to area. This problem is compounded by the fact that in mining environments the pillars are not always cut according to the prescribed dimensions. Although the original design using either empirical methods or modelling may predict stable pillars, unstable pillars will be encountered in reality. Very conservative designs with large factors of safety may circumvent this problem, but this approach is uneconomical.A possible solution to localized stability problems may be to enhance the strength of a few unstable pillars. It may even be hypothesized that reinforcing a few critical pillars may prevent ‘pillar runs’ on a much larger scale. Rockbolting, strapping of pillars, and pillar shotcreting have occasionally been used in the past as possible solutions. It appears that these have not been successful in all cases and large collapses have occurred in spite of the pillar remedial work.The authors investigated the strapping of pillars by conducting laboratory tests on cabled rock specimens. The results were qualitatively compared to actual attempts of pillar reinforcement available in the literature, as well as additional observations in an old haematite roomand- pillar mine in Spain. Based on these results, the value of steel cabling and mesh wrapped around pillars to improve stability is demonstrated. Some cases where this approach will not be successful are also discussed."

Jan 1, 2017

By F. Mwaura

"The aim of the audit was to determine whether the mining project environmental impact assessments (EIAs) in Kenya are undertaken according to international best practice. A sample of 50 EIA reports for the 2007–2016 period was evaluated using 18 criteria in the Environmental Law Alliance Worldwide (ELAW) guidelines. The findings showed that only two criteria were considered ‘excellent’ according to the ELAW guidelines, namely baseline environmental assessment and prediction of physical impacts. Six criteria were found to be ‘satisfactory’, including prediction of social impacts, analysis of alternative options, and impact mitigation. The reports were quite weak in terms of: consideration of all phases in the mining cycle, environmental regulatory framework, prediction of biological impacts, stakeholder consultation and engagement, integration of human right issues, and integration of climate change and cost-benefit analysis. It is therefore recommended that the Kenyan National Environment Management Authority (NEMA) should consider tightening the EIA terms and conditions in the approval of statutory terms of reference (ToR) for full-scale mining EIAs to ensure improved performance of EIA as a tool for environmental protection. IntroductionThe nexus between society and the environment has continued to generate interest because the environment is the source of goods and services that sustain the various needs of human beings as outlined in the UN sustainable development goals (SDGs) (Armstrong and Peart, 2000; Sutton 2004; Chakrabarti 2007; Strange and Bayley 2008). All the five generic categories of societal needs, namely physiological needs (food, water, clothing, medicines, etc.), residential and occupational needs (shelter, transport, energy, etc.), economic needs (tradeable goods and services), leisure and recreational needs (tourism), and cultural and spiritual needs (e.g. sacred resources and ecosystems) (Mwaura et al. 2016) are addressed by the environment. Consequently, a wide range of development sectors have been established around the world in order to address the above needs (Lederer, Galtung, and Antal, 1980; Rammelsberg et al., 2006; Noonan, 2014). All these development sectors rely entirely on the environment for their proper functioning, as shown in Table I."

May 1, 2019

"The management of outlier grades in positively skewed gold distributions is a contentious issue. Incorporating outliers in standard ordinary kriging (OK) estimation procedures in a way that honours the data without smearing extreme grades into surrounding areas has been problematic. Cutting or capping of outliers to mitigate their influence in estimation techniques is common practice, while methods that manipulate the OK system of equations fail to honour the data. We propose a method of post-processing of kriging weights that provides realistic OK estimates and mitigates smearing without manipulating kriging equations or changing the original grades. The method requires that the data is not clustered, is approximately equally spaced, and is of the same support. Positively skewed data is ordered on attribute grade and nonlinearly transformed to a Gaussian histogram of categorical bins whose frequency is based on their likelihood of occurrence and location in the sample distributions. Factors that restrict kriging weights are calculated by dividing the percentage frequency of data in each bin by the percentage frequency of data in the bin with the highest frequency. Restriction factors applied to the kriging weights in the OK estimation restrict the range of influence in proportion to their probability of occurrence in the distribution. Smear reduction post-processing is easy to implement and addresses issues arising from negative kriging weights while considering the spatial location of samples, the sample grades, and their probability of occurrence. The method mitigates both smoothing and conditional bias.IntroductionOutliers are sample observations that deviate considerably from the standard or expected (Hawkins, 1980). These deviations can be high or low in grade depending on the distribution. The presence of outliers in a mine sampling campaign could provoke mixed reactions. Mining engineers might view high-grade outliers as a promise of achieving higher grades, while the geostatisticians view outliers as affecting the accuracy and precision of the overall estimate, among other issues. Having identified outliers, the geostatistician must decide how to account for their presence and treat them so that true underlying estimates are not distorted. Advanced knowledge of the deposit allows the geostatistician to identify and treat outliers acceptably, especially where samples are limited in number and the true grade distribution is uncertain. If ignored, or treated incorrectly, outliers are likely to lead to smearing of extreme grades into the surrounding estimates, thereby steering regions towards potential under- or over-estimation."

Apr 1, 2019

By B. G. Russell

A wide range of minerals is found within the Republic of South Africa, many of which are exploited to a greater extent than in similar deposits found elsewhere. It is often asked whether this is atributed to the high grade of the ores encountered in this country, or whether such factors as fiscal and Infrastructural incentives create an economic climate that has en bled the mining industry to successfully extract and market minerals at a profit. This paper views South Africa's mineral resources in a regional and global context In an attempt to place the country's geological wealth in perspective. It represents the author's own Views, and not necessarily those of the Department of Mine al and Energy Affairs.

Jan 1, 1986

By D. K. Chibwe, P. Taskinen

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

Jan 1, 2014

By E. C. de Beer

"SynopsisA number of important areas in nuclear fuel cycle, at both the front end and back end, offer ideal opportunities for the application of non-destructive evaluation techniques. These techniques do not only provide opportunities for non-invasive testing of e.g. irradiated materials, but also play an important role in the development of new materials in the nuclear sector. The advantage of penetrating radiation used as probe in the investigation and testing of nuclear materials makes X-ray and neutron radiography (2D) and tomography (3D) suitable for various applications in the total nuclear fuel cycle. The unique and different interaction modes of the two radiation probes with materials provide several opportunities. Their complementary nature and non-destructive character makes them most suitable for nuclear material analyses, analytical method development, and the evaluation of the performance of existing nuclear material compositions. This article gives an overview of the X-ray and neutron radiography/tomography applications in the field of nuclear material testing, and highlights a few of the success stories. Several selected areas of application in the nuclear fuel cycle are discussed to illustrate the complementary nature of these techniques as applied to nuclear materials.IntroductionDuring the development of new materials for the nuclear industry, materials testing and characterization are of the outmost importance to maintain safety standards and reliability. No compromise on safety in the workplace in any area within the total nuclear fuel cycle can be tolerated, and therefore most in-situ material characterization and testing is conducted by certified and qualified personnel schooled in destructive and non-destructive testing (DT and NDT) methods. Certification and qualification in NDT can be obtained through many training centres in South Africa in accordance with European, American, and other international standards (SGS, n.d.; SAIW, n.d.; African NTD Centre, n.d.)."

Jan 1, 2015

By E. Lutz, M. Momayez, K. W. Brown Requist

As real-time air quality monitoring becomes more prevalent in US underground mines, it is important to provide the highest data reliability with the fewest possible sensors. Real-time sensors remain costly, and these costs are not exclusively financial; the time required to install, calibrate, and maintain real-time sensors poses a large barrier to widespread implementation. Current atmospheric monitoring systems typically rely on displaying point-specific values. This requires operators to infer real-time airborne contamination distributions. Monitoring and control software utilizing mine ventilation network (MVN) solvers has been implemented in limited cases because of their ability to simulate ventilation systems quickly, but these solvers use a one-dimensional representation of the mine, limiting spatial resolution of estimated distributions. Computational fluid dynamics (CFD) has likewise been considered as a means to improve spatial resolution, but processing times prevent its use as a basis for monitoring and control. For the real-time monitoring of airborne contamination distributions, we propose a spatial interpolation method that can estimate the distribution of airborne contaminants in near-real time. This method provides a middle ground between fast processing times and increased spatial resolution. With the use of a pathfinding algorithm and optimization through absolute percentage error minimization, this method outperforms spatial interpolation with a Euclidean distance. By providing contamination distribution information to operators, this method and its derivatives stand to outperform current atmospheric monitoring systems.

Apr 16, 2024

By M. S. Shishvan, J. R. van Duijvenbode

The results of dig limit delineation in open pit mining are never truly optimized due to gaps in the underlying data, such as insufficient sampling. Aside from the data uncertainty, there is also an influence on the final dig limit by either humans or by the heuristic character of an optimization method like simulated annealing. Several dig limit optimizers have been published, which can replace the manual dig-limits designing process. However, these dig limit designs are generally not adapted to account for this heuristic character. In this paper we present a stochastic analysis tool that can be used with the results of heuristic dig-limit optimization to increase confidence in the obtained results. First, an enhanced simulated annealing algorithm for dig limit optimization is presented. Then, this algorithm is tested on ten different blasts at the Marigold mine, Nevada, USA, as a case study. Finally, the results are analysed with a destination-based ensemble probability map and an analysis conducted of the final solution data distribution. The generated dig-limit designs of the algorithm include high revenue areas that are excluded in comparable manual designs and show improved objective and revenue values. The analysis tool provides block destination probabilities and box plots with the distribution of opportunity value for the dig limit. Furthermore, with the analysis tool, it is possible to make well-informed design decisions in areas of uncertainty.

Dec 13, 2022

By I. Ngo, L. Q. Ma, J. S. Guo

Deformation of gob-side entries has always been a critical concern for ensuring stability in longwall coal mines. This paper addresses the significant deformation and support challenges that arise in thick and hard roof longwall faces (THRLF) due to dynamic pressure. The study aims to elucidate the characteristics and mechanisms of deformation during the retreat of the longwall face. The research findings indicated that the primary cause of deformation was the combination of advanced abutment stress resulting from longwall face mining and the movement of the lateral roof over the chain pillar. To mitigate this issue, we propose a deformation control method known as cutting off the lateral roof (COLR) over the chain pillar. Simulation results demonstrate a significant reduction in roof stress and deformation of the gob-side entry after implementing the lateral roof-cutting technique. These findings provide valuable guidance for effectively managing deformation in gob-side entries, particularly when dealing with thick and hard roof strata.

Aug 3, 2023

By T. Dlokweni, M. Bester, C. Koegelenberg, I. Basson

Certain degrees of safety, economic, and financial risk are implicit in any mining operation. At Thabazimbi Mine, slope stability was one of the major sources of risk, largely due to data uncertainties. Consequently, a risk-based approach in slope design was followed. On 6 June 2015, a slope failure occurred in Thabazimbi Mine’s Kumba pit, involving approximately 65 million tonnes of rock failing into the pit and the valley below. Although the failed mass moved rapidly and resulted in a run-out distance of approximately 900 m, controls had been put in place to successfully evacuate all personnel. Risk assessment during the onset of failure supported the decision to not compromise the safety of personnel in any way, leading to unavoidable equipment loss. This contribution presents the failure as a case study in the context of a risk-based approach to slope design as well as the importance of implementing geotechnical controls to effectively manage slope instability risk. Furthermore, a comprehensive back-analysis was performed, on a fully constrained 3D model, utilising 3DEC software, to gain insights into the failure mechanism. In conclusion, lessons from the Thabazimbi slope failure are valuable as they demonstrate the importance of following a risk-based approach in slope design to effectively manage safety and financial risk in open pit mines.

Sep 1, 2025

By A. Eitzenberger

This paper describes a study of seismic records from the Kiirunavaara mine footwall which were interpreted in relation with numerical models developed outside the study. Seismic data was retrieved from a portion of the mine and filtered with respect to the ratio between energy carried by an event's P (primary) and S (secondary) waves (Es/Ep ratio), local magnitude, and active mining depth. The data was analysed using Es/Ep ratios and Gutenberg-Richter graphs to determine the event origin, mechanisms, and minimum magnitude cut-off. The magnitude completeness was identified by studying the b-value stability and b-value differentiation between origin sets. It was shown that, by separating seismic events into the origin components shear, complex, and tensile based on Es/Ep ratios, a representative value for the magnitude completeness can be identified for a catalogue with a convex cumulative log curve. The majority of the events were shown to be of shear-slip origin based on the recorded Es/Ep ratios, with pure tensile events constituting only about 10% of the recorded data. Spatial and temporal event location patterns were studied and compared with numerical modelling results. The comparison showed a correlation between shear-slip seismic events and volumes experiencing high differential stresses in the lower part of the footwall. In the upper part of the footwall the results did not reveal any clear correlation between observed damage in drifts and seismic event locations. The concentration of seismic events in the lower portion of the footwall is discussed in the context of rock mass displacements. The results indicate a possible connection between mine seismicity at depth and damage observations in the drifts in higher non-seismic areas by seismic softening and subsequent lateral expansion of the rock mass.

Jan 1, 2017

By H. M. Friede, R. H. Steel

Notes on Iron Age copper-smelting technology in the Transvaal by H. M. FRIEDE, Ph.D. and R. H. STEEL (Visitors) Earlier reports on Iron Age copper metallurgy and copper-smelting techniques in the Transvaal are surveyed, and a description and tentative classification of the copper furnaces are given. The analyses of furnace slags and sherds of glazed crucibles from archaeological sites of the Transvaal and the neighbouring Orange Free State pro. vide evidence of copper smelting there in Iron Age times. . An investigation conducted by the Department of Archaeology of the University ofthe Witwatersrand provided information on the details of furnace construction and smelting processes, on the composition of the materials used and produced, and on the making and handling of tools and accessories, such as bellows and tongs, needed by Iron Age smelters. Three models (the KaondejVenda furnace, the Rooiberg crucible furnace, and the Uitkomst furnace) were constructed and operated as a demonstration of Iron Age copper smelting as it was practised in the Transvaal.

Jan 11, 1975

By L. A. Abichequer

ROM grade variability is one of the most important factors influencing plant product recovery. Due to the influence of grade variability on mine profits, it is deemed important to guarantee a proper homogenization of the ore through a well-planned blending and homogenization pile system from which the ore is reclaimed at an acceptable grade fluctuation. This study tests geostatistical simulation efficiency in predicting in situ grade variability. The methodology is illustrated using a data-set from a large phosphate mine in central Brazil. Additionally, real data obtained at the mine site during January and February 2007 were used to derive the actual grades of each pile formed during this period. Each pile was formed by a selection of blocks within a dig line and these lines mined along the year 2007 were also made available. Given the real grades, it was possible to reconciliate and validate the results predicted by the proposed methodology. The main outcome of this study is the opportunity to check the geostatistical simulation as a tool to quantify the variability of the homogenization system.

Jan 1, 2009