Search Documents

By S. Banisi, M. Maleki-Moghaddam

Direct observation of charge shape and its motion in industrial mills is not possible; therefore, a combination of analytical and physical studies is used to determine charge trajectory. Software packages, which predict charge motion consider only the outermost layer of charge (single ball) and ignore the charge and the interactions of grinding elements. In the new version of the GMT (Grinding Media Trajectory) software, the effect of liners specification and bulk load on charge shape and trajectory are considered. The charge trajectory was measured in a laboratory mill with the diameter of 100 cm and length of 21.5 cm with the transparent end and compared with that of the old version of GMT. Four types of polyurethane rings were accurately machined to scale down the liners arrangements at two industrial mills. The model mill was operated at 55, 70 and 85% of the critical speed for five levels of mill filling (10, 15, 20, 25% and 30% by volume). The experiments indicated that the end walls effect was negligible for the model mill with the length of 10.8 cm and beyond. The proposed relationships to correct the trajectory and charge shape were validated by using the new liner of the Gol-E-Gohar iron ore company AG mill. The positions of shoulder, toe, toe departure and impact points were predicted with the relative errors of 5, 0.7, 2 and 0.2%, respectively. The GMT software was updated with the new relationships. The results of GMT simulation indicated that when the lifter face angle increased from 7 to 30°, the distance between the charge impact point and the toe decreased from 40.1 to 11.2° for 30% filling. The change in the position of the impact point meant increasing the probability of charge impacting the toe, not the liner which favored more efficient comminution practice.

Jan 1, 2016

By R. Mohamad, J. Asfahani

Abstract - Radioactive characteristics of pits drilled in the Al-Sharquieh mine while prospecting for recent phosphatic resources have been investigated. The results indicate that phosphatic layers are characterized generally by high radioactive intensities (more than 800 c.p.s). The lithological section in the study area is very complex sedimentologically, and characterized by high faciological variation in all directions. The Schlumberger geo-electrical configuration with fine separation of AB/2 was successfully applied to determine lithological boundaries, and the resistivities of the layers of the lithological section, especially for thin layers. It was shown that the phosphatic layer resistivities increase or decrease as a function of associated material. As a result, two types of phosphatic layers have been geo-electrically distinguished. The first is characterized by high resistivity and relatively low radioactivity, which could be related to the rock phosphatic layers. The second is characterized by low resistivities and high radioactivity, which could be related to the sandy phosphatic layers.

Jan 1, 2000

By Charles A. Schacht

"Relining, lining maintenance, and ladle repair can represent a significant cost in steelmaking. Therefore, the objective of all steelmaking shops is to minimize this cost and to maximize the lining life. However, the thermomechanical behaviour of the ladle lining is highly complicated and sometimes misunderstood. The ladle design also has a significant influence on lining life. Other factors which have a direct influence on lining life are the shop lining installation practices and manner in which the ladle is used. The objective of this paper is to provide an understanding of the impact of these factors as related to improving the integrity of the lining and increasing lining life.Fundamentals of Steelmaking Ladle Lining Expansion ForcesNumerous studies have been conducted by the author on steelmaking ladles with molten metal capacities ranging from 110 ton s to 350 tons. The steelmaking ladle typically has a sidewall plate with a slight taper and for analytical purposes can be classified as a cylindrical vessel.A typical ladle is described in Figure I. The ladle sidewall lining usually consists of a working lining and a safety lining. Although the type of refractory used in the working lining varies based on the type of steel produced and other metallurgical requirements, the refractories typically used are: 70% alumina brick; direct bond mag-chrome brick; resin bond dolomite brick; and direct bond dolomite brick. Fireclay brick, high alumina brick, or a similar strength brick is used as the safety lining.Normally, the ladle lining is exposed to a through-thickness temperature gradient and restrained by the ladle sidewall plate, is described in Figure 2. The hot face side of the refractory lining develops a compressive thermal force and the ladle sidewall plate develops an equilibrating tensile force. The radial joints open on the cold side of the working lining as the lining joints have little or no tensile strength. A radial compressive load (P) develops between the working lining and the vessel shell due to the thermal restraint imposed on the lining by the shell. The radial compressive load is transferred from the working lining through the safety lining onto the ladle sidewall. Therefore, the safety lining must have the strength to accommodate the radial compressive load (P)."

Jan 1, 1991

By Amir Karami, Derek Kinakin, Catherine Schmid

The Mathews Stability Graph Method (Mathews et al., 1981) has been used by geotechnical engineers and mining engineers for over 40 years to estimate the performance of longhole open stopes in underground mines. The method is widely used and accepted by practitioners because it is relatively simple to use and can be modified to account for observed performance at a given mining operation. While the simplicity of the method increases its accessibility and utilization, it does not intrinsically account for variability in the various inputs, including rockmass quality, structure, stress, or stope orientation. This paper presents a modified approach for the Stability Graph Method, whereby Python scripts are used to process available geotechnical datasets and combine them with actual or potential distributions of Q parameters, induced stress, structure, or stope orientation. The resulting Stability Numbers are then plotted as a contoured distribution on the Stability Graph.  This allows the practitioner to quickly incorporate large datasets into the analyses, while also developing an understanding of how variability in rockmass quality, induced stress, structure, or stope orientation could affect stope stability.  This paper presents the fundamental principles of the Stability Graph method, the Python script and workflow, and two case studies to demonstrate the utility of the tool.

Oct 1, 2024

Thermal barrier coating (TBC) systems provide thermal insulation to high temperature for aero¬propulsion. A fracture mechanical model developed earlier for multilayer systems is applied on TBC to realize superior performance by enhancement of fracture resistance. Functionally gradient materials (FGM) solution concept is the basis of approach followed here. A few realistic FGM TBCs are considered in the analysis in addition to a conventional TBC as the base case and a hypothetical worst case with YSZ-YSZ interface. Parameters like number, composition and thickness of graded layers and elastic modulus are considered for the analysis and quantitative improvement of fracture resistance of TBCs. The analysis is confined to deep crack with crack size to width ratio as 0.5.

Jan 1, 2010

By L. Xu, A. Vazirnegad, R. Bhaskar

"Respirable dust samples are collected by the U.S. Mine Safety and Health Administration (MSHA) and mine operators as part of the requirements of the Federal Coal Mine Health and Safety Act of 1969. About 1.85 million respirable dust samples collected by MSHA were analyzed for the period 197/ to 1990 f or the states of Alabama, Colorado, Illinois, Kentucky, Ohio, Pennsylvania, Utah, Virginia, and West Virginia. Trends in coal dust exposure of continuous miner operators, helpers, and roof bolters were analyzed. The reported production in each of the nine states was analyzed. The relationship between average dust exposure and coal production was evaluated. In addition, the differences between mine operator and MSHA inspector collected data were examined for each state. The average exposure levels for the three occupations - miner operator, miner helper, and roof bolter were also compared.The data indicate that the dust exposure levels have significantly dropped over the last twenty years, while at the same time production has increased, attesting to the ability of the U.S. mining industry to adapt successfully to the stringent dust control regulations. In terms of coal production, it appears that some states, such as Utah, Illinois, and Kentucky, were able to increase their production significantly relative to others, such as Pennsylvania and Virginia. The relationship between dust concentration and coal production is, at best, tenuous. 11 implies that existing dust control technology has assisted several mine operators to meet dust standards, while, at the same time, enabling mines to produce more coal than before. Comparison of data collected by MSHA and mine operators indicate that MSHA inspectors recorded higher dust levels and lower shift production relative to mine operators."

Jan 1, 1993

By A. M. Reyes

Flow field hydrodynamics of serpentine flow channels was investigated numerically using computational fluid dynamics (CFD). The velocity profile simulated by the CFD model was evaluated with respect to experimental results from the literature. Rectangular and tapered flow channels were compared with respect to mixing and water management capabilities, and pressure drop. The CFD model was further developed to account for electrochemical reaction on the cathode side to determine which channel configuration was most effective in reactant conversion. The velocity profile differences between the channels were found to impact reactant conversion, as the rectangular channel had the highest conversion. With respect to water management, although the tapered channel was found to have a higher overall pressure drop, the rectangular channel was found to have higher shear and pressure forces acting on simulated liquid water droplets in the channel.

Jan 1, 2005

By M. Silva, G. Olivares

In a leach pile, the drainage layer or commonly known as the cover layer has been a matter of discussion due to the multiple phenomena or events that occur within it which impacts negatively the leaching process. All these events are directly related to the loss of draining capacity and irregular flow net due to an excessive compaction and entrapment of fines reducing the flow velocity which in turn cause more fines to get trapped. Over time the problem will continue to worsen rapidly and accelerating continuously like a geometric progression. When the cover layer gets to a certain level of deterioration, the only alternative to continue operating safely is to replace the entire cover material, a huge cost for any mining operation. TESRA has developed a solution to counteract the accelerated advance of the deterioration of the drainage layer. The solution consists of using large agricultural tractors coupled to a proprietary design of a heavy- duty subsoiler-ripper for deep tillage for loosening and breaking up the ground at various depth and with different legs spacing from 0.4 m or above. Also, if the ground still requires additional breaking, TESRA proposes the use of a rototiller that breaks any lumps to homogenize the ground to a certain particle size. With the subsoiler-ripper/rototiller combination equipment a layer of 0.5 m deep is properly treated to match closely the original conditions of the ground thus regaining the flow net as close as to the original seepage condition. Once the ripping is completed the subsoiler-ripper tool is replaced with another that is similar to a plow and thus continue with the installation of drainage pipes. The tool digs a Y-type profile trench allowing the pipe to sit precisely at the bottom of the trench, contributing additional lateral and lower support. This tool controls the installation depth of the drainage pipe at different inclination grades and also eliminating any irregularities that may block the flow of the leaching solution. Sometimes the obstructions block completely the flow of the solution. The main result of this methodology is the reduction of groundwater levels and at least 2% less humidity in the leaching gravel due to an increase in the draining capacity of cover. The methodology developed by TESRA applied during a couple of cycles of a dynamic leach pile has proven to be a real alternative to remedy and recuperate the cover close to its original conditions. The methodology requires only one large tractor and two accessories operated in different stages while replacing the more traditional and less efficient machinery such as the Bulldozer, Motor-grader, and backhoes.

Jan 1, 2019

By A. Abadías Llamas, N. Bartie, M. Stelter

Evaluating the economic viability as well as the sustainability of the Circular Economy (CE) system requires a deep understanding of the distribution of all elements, compounds, alloys, materials etc. in flows. In this paper, the circularity of the copper value chain, including primary and secondary processing, is rigorously evaluated. The studied system comprises the metallic copper production from primary sources (from mineral to metal), copper-containing commodity production (copper is mixed with other metals) and copper recycling through secondary smelting to close the loop. This is linked to photovoltaic (PV) panels and battery storage. A simulation model of this system is created using HSC Sim, considering more than 30 elements (and its various compounds), 180 unit operations and 800 flows. From the mass and energy balances obtained through the simulation, an exergy analysis is conducted to evaluate the resource consumption from a second law (entropy) perspective. Additionally, these results are complemented through a Life Cycle Assessment (LCA), the recovery of technology elements and by-products is discussed, while quantifying the losses through the value chain. Through the digitalization of the complete system, a better CAPEX and OPEX understanding of the metal recovery and losses can be obtained, as well as the associated resource consumption and environmental impacts. New flowsheets and technologies can be evaluated. Several scenarios show how the resource consumption and the environmental impacts are affected by the recovery of different materials to produce different products. INTRODUCTION A key challenge of the circular economy (CE) is the linkage and optimization of both resource and energy efficiency of the (renewable) energy infrastructure and the resource cycle depicted by Figure 1 simultaneously. In addition, the CE requires the linkage of all stakeholders as shown by Figure 1. This has been discussed in some detail by Reuter and Bartie (2019) and Reuter et al. (2019). Copper and all its associated technology elements is one of the key enabling metals that is an economical copper colored thread through CE society (Reuter & Kojo, 2014). The Metal Wheel (Figure 2) shows the linkage of and recovery of various elements associated with copper and other carrier metals. This reveals that copper acts as an important carrier metal because of its properties i.e. the wealth of technology elements it carries from geological ores, but also its capacity to act as solvent of valuable elements and thence a key closer of various complexly connected element and material cycles when these technology elements are released during copper refining. The same applies to all other carrier metals of modern CE society as shown for example recently in a lead policy brief (Blanpain, Reuter & Malfliet, 2019).

Jan 1, 2019

By P. Mottahed, J. Ran

Mining practice and research have shown that the jointed roof in stratified rock behaves like a voussoir beam structure with strength greater than that determined by the continuous beam theory. Previous studies have demonstrated the voussoir beam mechanism and developed methods to estimate the crushing and buckling strength of voussoir beams. However, research on voussoir beams with multiple mid joints has not been suffICient and problems such as the voussoir beam support and applications of the voussoir beam mechanism to the mining environment have not been addressed. In this paper, mechanism of voussoir beams containing multiple mid joints and truss support are studied and applications of the voussoir beam mechanism to underground mine openings are discussed based on the numerical and experimental results.

Jan 1, 1995

Mining Impact in Water Quality: Understanding the conflict Mining is an industrial process Metal mines that chemically process ore to concentrate metals such as copper and gold use more water than non-metal mines such as coal, salt and gravel. Tailing Ponds Water contaminated with high concentrations of metals, sulphide minerals, dissolved solids, or salts can negatively affect surface water quality, aquatic ecosystems and groundwater quality.

May 1, 2013

Why innovate? ? How are we doing in Canada? ? How can we do better? ? Why now? ? What?s next? ? Lowering production costs ? Improving productivity ? Increasing value-added ? Discovering new mines, extending life of existing ore bodies ? Reducing environmental footprint ? Improving health and safety

May 1, 2008

By Jan Rosenkranz, Anders Sand, Lisa Malm, Nils-Johan Bolin

Within flotation technology the general trend is towards very large cells. Scale-up is usually done by postulating geometrical similarity, which implies keeping the aspect ratio constant when increasing cell volume. Consequently, the design of flotation circuits typically involves several cells of identical geometry within a bank. Using a few standard sizes in a flotation plant simplifies design, manufacturing and maintenance of the cells, but does not necessarily guarantee optimal performance and selectivity of the flotation process. Geometry parameters for a given cell volume that influence flotation cell performance include cell height, which causes changes in the hydrostatic pressure and suspension hydrodynamics, and also influences the travel distance of particle-bubble agglomerates as well as the homogeneity of the mixing. Also the thickness of the froth layer depends on the aspect ratio of the cell. The objective of this work is to gain a better understanding of material distribution and properties within a flotation cell and, based on this information, elucidate how alterations of the cell geometry can influence separation efficiency. For this purpose systematic measurements and analyses of the spatial distribution of the different phases within a 160 m3 flotation cell have been carried out with a particular focus on an industrial, low-grade copper ore beneficiation process. The concepts for sampling at different vertical and lateral positions of the flotation cell are introduced. The results from the experimental work give an understanding of internal material distribution within flotation cells. Based on the measured phase distribution and concentrations as well as particle properties at various positions within the cell, the implications for an optimized scaling and design of individual flotation cells within a flotation circuit are discussed.

Jan 1, 2016

By Somasundaran P., M. Vasudevan, Z. Dai, D. R. Nagaraj, T. Bhambani, P. Patra

Historically, efficient beneficiation of ultramafic ores has been challenging, partly due to a lack of clear scientific understanding of the effects of Mg silicates, particularly fibrous serpentines, on selective Ni recovery. Recent research at Columbia-Cytec as a part of a Vale JDP project has led to new insights into the mechanistic aspects of the adverse effects of Mg silicates. Our research has provided strong evidence to support the hypothesis that the adverse effects of Mg silicates are mediated through their large influence on pulp rheological properties, far more than any other pathway. The severity of the effect of Mg silicates is predominantly dependent on the aspect ratios of Mg silicates. For example, the pulp viscosity increases dramatically with increase in the fibrous nature of the gangue mineral leading to poor Ni separation. Various ways by which the suspension rheology affects flotation sub-processes, with implications for remedy are discussed. The insights and results have served to develop (or construct) a phenomenological model to adequately explain all of the ore flotation data, particularly the Mg silicate transport to the froth phase and the concentrate.

Jan 1, 2011

By A. B. Irwin

"Discussion at MeetingDR. D. DERRY (Ventures Limited) referring to the diagram (Figure 15) which shows a transverse structure section across the Central Slocan area, asked what the assumed structure at depth was, remarking that dips are steeper and more uniform at depth and that the contortions with gently-dipping axial planes are at higher elevations. He wondered whether the section represented ,part of a very large fold with steeply-dipping axial plane.MR. J. D. LITTLE (Canadian Exploration, Limited) replied that E. B. Mayo (Kelowna Exploration Company) and others had considered that the entire section illustrated faces to the southwest and represents part of the lower limb of a great overturned •syncline which is concave toward the southwest.DR. R. DEWIT (Geological Survey) remarked that the photograph of a ripple mark (Figure 7) showed internal lamination conformable to the ripple mark and wondered if the structure was not due to slumping during consolidation, since the normal mode of formation of ripple mark would seem to preclude the possibility of such internal stratification.MR. LITTLE remarked that in his opinion such structures should not be depended upon alone with-out other confirmatory evidence of the direction of tops of beds."

Jan 1, 1952

By N. Schmitt, E. Blond, F. Hild

Working linings of steel ladles progressively damage during the steel manufacturing process due to both severe thermomechanical and chemical stresses. In the framework of the coupling scheme proposed by Blond et al. (1), this study examines the effects of the thermomechanical loading on the onset of microcracking by considering that the complex microstructure of the slag-impregnated brick is fixed during the steelmaking process. The behaviour of the refractory is identified by using thermoelastic-(asymmetric)viscoplastic constitutive equations. Numerical simulations of a simplified geometry of the lining subjected to cyclic thermal loadings at the inner side are shown. Two cases are analysed, namely, i) the material is a homogeneous continuum with a thermoelastic-(asymmetric)viscoplastic behaviour, ii) the material is porous and saturated by liquid slag, and the skeleton has the previous non-linear behaviour. The stress and strain fields are shown, and for the porous medium the interstitial pressure of slag as well. Comparisons between these assumptions are made to assess the role of slag on the degradation of the working lining.

Jan 1, 2004

By Jin-Hwa Song

Rolling temperatures and cross rolling ratios were varied to investigate the effect of crystallographic texture of alpha-phase on the anisotropy of mechanical properties in modified Ti-6A1-4V alloy. Hot rolling method was unidirectional and cross rolling, during which the rolling temperature was in the range of 830°C and 920°C. Cross rolling ratios were varied from 0.3 to 1.0 at 875°C. For unidirectional rolling, as the rolling temperature increased texture changed from the mixed basal and transverse texture to the simple transverse texture. Contrary to other work, it was observed that the maximum intensity of transverse components increased with the increase in rolling temperature. For cross rolling, on the other hand, the transverse texture component was almost suppressed at 830°C, giving rise to a symmetrical basal texture. With increasing cross rolling temperature, the transverse texture components developed again. A simple transverse texture, which showed two maxima appearing circularly 90° apart, was observed at 920°C. It was found that unidirectional rolling could not eliminate the transverse texture, which resulted in the anisotropy of mechanical properties. The anisotropy of mechanical properties obtained from tensile and V-notch Charpy test was closely related to the crystallographic texture and its intensity at each rolling condition. To obtain the isotropy of mechanical properties by cross rolling, it is suggested that cross-rolling ratios should be controlled above 0.5. The variation of rolling temperatures and directions has great influence on the anisotropy of impact value rather than tensile strength.

Jan 1, 2004

By M. Mahdi Aboutalebi, Roderick I. L. Guthrie, Mihaiela Isac

"A preliminary computational fluid flow model has been developed to simulate the flow field of liquid steel in the mold region of billet mold casting. This study reports the effect of different Submerged Entry Nozzles (SEN) on the fluid flow in round and square billet, curved mold, casters. In this study, ANSYS Fluent 14.5 software was used to model the 3-D turbulent flow of liquid steel. The standard k-0 model was used to simulate turbulence in the mold region. Similarly, the species transport equations were discretized to simulate the dye injection test into the mold cavity. The behavior of flow in the mold region was validated against previous experimental work. Predicted results showed that different nozzle types plus the slight curvatures in the molds can significantly influence the flow patterns generated within the mold region, and thereby affect the finished quality of the billet products.INTRODUCTIONConventional Continuous Casting (C.C) processes are presently used to solidify more than 95% of the 1.5 billion tons of steel now produced in the world annually (World Steel in Figure, 2014). This overwhelming contribution of C.C processes for steel production is the prime motivation for this study, in order to understand and optimize the casting processes and thereby reduce internal defects and minimize energy consumption. The primary objective of all C.C processes is to maintain liquid metal purity during the casting process and to create a uniform flow within the caster so as to avoid internal and external defects. In this regard, studying fluid flows in the mold region of continuous caster has attracted many research groups, in both academia and industry, in recent years. The mold region is one of the most critical parts of the C.C. process, in which the transport phenomena associated with this region can significantly affect the quality of the final products. Thus, proper control of fluid flows within this region can be expected to greatly decrease both internal and surface defects such as Oscillation Marks (OM), inclusion entrapment, and Mold Powder Entrapment (MPE). In the casting process, the fluid flow in the mold region can be affected by a number of variables, such as; the flow control system (stopper rod or slide gate), nozzle geometry, SEN (Submerged Entry Nozzle) depth, casting speed, mold dimensions, gas injection, and the application of electromagnetic fields. One of the most influential and easily-changed fluid flows are the jets of steel issuing from the SEN. Nozzle geometries are variable, consisting of a number of exit ports, different port angles, and/or port sizes, the bottom shape of the nozzle, and the SEN immersion depth (Thomas, 2001). The SEN characteristics can markedly affect meniscus surface velocities, surface level profiles, and turbulence levels that can vary in terms of the frequency and the magnitude of their fluctuations."

Jan 1, 2016

By Z. Pilecki

This work concerns the identification of weak zones in the structure of a rock mass endangered by outburst of gases and rocks in the Rudna copper ore mine at greater depths (approximately 1200 m). A special methodology has been undertaken for the recognition of that threat in the wake of an outburst in the mine in 2009. No similar phenomena have been previously recorded in KGHM Polska Miedz copper ore mines. The research described in this paper was carried out in the new part of the deposit prepared for mining. Geophysical methods were used: seismic tomography and borehole ground penetrating radar (BGPR), along with particular identifications using geological engineering methods. In the researched area, Permian layers of sandstone and shale clay (Weissliegendes formation) and of Zechstein dolomite (top of Rotliegend formation) and anhydrite occur in the copper ore deposit. The rock layers are characterized by variable thickness and intensity of mineralization, varying internal structure, and local fracturing containing different gases and water. It is assumed that accumulation of gases occurs in the tectonically fractured dolomite. The study was focused on the recognition of the structure and properties of the dolomite layer lying directly within the roof of excavations. The first tests were performed with seismic tomography for the recognition of anomalous changes in the elastic properties of the rock mass. The results of these studies were used to design a control borehole for a detailed geological engineering assessment and BGPR surveys with 100MHz antennas. Seismic tomography enabled the identification of changes in Pwave and S-wave velocity, the Vp/Vs ratio, dynamic Poisson's ratio and dynamic Young's modulus. It was clear that seismic anomalies required verification by other methods. The BGPR survey was performed in a specially designed borehole, passing through the selected seismic anomalies. While the borehole was being drilled, the outflow of gases was identified, with varying amounts of CH4, CO2, and H2S detected, as well as water outflow. The results of BGPR confirmed the location of a reservoir of water and gases, about 30 m in length, found in layer of dolomite. The direct cause of the BGPR anomaly was assumed to be intensive fracturing, a greater content of clay material in the dolomite, and the presence of water in the fractured zone. The BGPR image was correlated with information on lithology, fracturing based on the RQD index and additional information acquired during drilling about water and gases surveys and other symptoms. The achieved results also showed a very strong correlation between the location of the zone containing gases and water and the anomalous S-wave velocity and field of dynamic Poisson?s ratio, and to a lesser extent other seismic parameters. It was found that the other measured seismic anomalies may be of the same nature as those identified in the control borehole and accordingly will require the implementation of detailed identification. The performed tests have also shown the usefulness of geophysical methods to identify geodynamic hazards under the geological and mining conditions of the copper ore deposit in the Rudna mine.

Aug 1, 2013

By Miguel A. Cuba

Trends in the mean and/or variance are a longstanding challenge in mineral resource modeling. Geologic modeling and/or grade domaining are often employed to account for heterogeneities in mineral grades; however, this is only practical up to a certain point since the number of data available for reliable model construction can be severely reduced. In such cases, trends remain present in certain geologic/grade domains and directly impact the modeling methodology; the trend must either be handled explicitly or a suitable approach is adopted to implicitly account for it. This is highly dependent on understanding the source of the trend. We propose an approach that detects the source of non-stationarity, which is based on quantifying the impact of non-stationarity in the mean and the variance via the semivariogram. By expanding the methodof- moments approximation of the semivariogram, we decompose the experimental semivariogram is into three terms: 1) trend in mean component, 2) trend in variance component and 3) stationary component. We can then assess the source of trend by calculating the contribution of each component at every lag. Two cases are presented to illustrate this tool. Firstly, a synthetic reference model is simulated and different trend patterns are added; the influence of these patterns in the experimental semivariogram calculation are analyzed and discussed. Secondly, a real dataset is considered where the impact due to the presence of trends in the mean and variance are assessed. The results of this paper reveal the semivariogram to be a tremendously powerful tool in determining the source of the trend. This insight may be useful in the selection of a suitable modeling approach for resource estimation.

Oct 1, 2009