Geological Modeling of Gold Deposit Based on Grade Domaining Using Plurigaussian Simulation Technique

Mineral resource evaluation requires defining grade domains of an ore deposit. Common practice in mineral resource estimation consists of partitioning the ore body into several grade domains before the geostatistical modeling and estimation at unsampled locations. Many ore deposits are made up of different mineralogical ensembles such as oxide and sulfide zone: being able to model the spatial layout of the different grades is vital to good mine planning and management. This study addresses the application of the plurigaussian simulation to Sivas (Turkey) gold deposits for constructing grade domain models that reproduce the contacts between different grade domains in accordance with geologist’s interpretation. The method is based on the relationship between indicator variables from grade distributions on the Gaussian random functions chosen to represent them. Geological knowledge is incorporated into the model by the definition of the indicator variables, their truncation strategy, and the grade domain proportions. The advantages of the plurigaussian simulation are exhibited through the case study. The results indicated that the processes are seen to respect reproducing complex geometrical grades of an ore deposit by means of simulating several grade domains with different spatial structure and taking into account their global proportions. The proposed proportion model proves as simple to use in resource estimation, to account for spatial variations of the grade characteristics and their distribution across the studied area, and for the uncertainty in the grade domain proportions. The simulated models can also be incorporated into mine planning and scheduling.     

Geochemical Anomaly and Mineral Prospectivity Mapping for Vein-Type Copper Mineralization,Kuhsiah-e-Urmak Area,Iran: Application of Sequential Gaussian Simulation and Multivariate Regression Analysis

Natural Resources Research - In this paper, sequential Gaussian simulation (SGS) and number–size (N–S) fractal modeling were used for copper geochemical anomaly mapping in the western...     

模型驱动数据空间分辨率对模拟生态水文过程的影响

设定两种独立的调参方案,检验将高/低分辨率下的调参值应用到低/高分辨率中驱动模型模拟的可行性;同时分析驱动数据分辨率的差异对模型模拟生态水文过程的影响;此外,设置9种集水面积阈值,分析驱动数据分辨率与集水面积阈值对模拟结果的影响。结果表明两种调参方案下模拟的水文与生态数据相关性强、差异小,说明调参方案受驱动数据分辨率的影响不明显;分辨率对土地覆被和土壤类型面积影响小,而对高程和坡度影响大;驱动数据分辨率降低,模拟径流的均值减小、峰值增大、生态变量值增大;驱动数据分辨与集水面积阈值之间不存在交互作用;集水面积阈值增大,模拟径流的均值、峰值均减小,生态变量值也减小,但阈值变化对峰现时间无影响。     

Cu- and Zn-Soil Anomalies in the NE Border of the South Portuguese Zone (Iberian Variscides,Portugal) Identified by Multifractal and Geostatistical Analyses

Extensive Cu- and Zn-soil geochemical data in the Albernoa/Entradas–S. Domingos region (NE border of the Iberian Pyrite Belt, South Portuguese Zone) were examined to separate anomalies from background using the concentration–area fractal model. Distribution patterns of Cu and Zn concentrations in soil are primarily influenced by bedrock. The regional threshold values of Cu- and Zn-soil contents over metasedimentary sequences are 20–25 and 20–60 ppm, respectively, becoming 30–50 and 20–90 ppm, respectively, when metavolcanic rocks are present. The first-order threshold values for Cu are 80–90 ppm in soils over metasediments and 70–80 ppm in soils over sequences bearing metavolcanics. For Zn, the first-order threshold values are 40–80 and 90–100 ppm in soils over metasediments and metavolcanic rocks, respectively. Metasediments and metavolcanics comprising significant sulphide disseminations are outlined by Cu- and Zn-soil values above 100 and 300 ppm in soil, respectively. On the basis of these results, Alvares and Albernoa/Entradas areas emerge as the first priority targets for exploration. The observed non-coincidence of Cu- and Zn-soil anomalies in soil in the area could reflect difference in element dispersion during weathering, they mostly indicate distinct metal sources related to the original composition of different rock types or to chemical changes developed during Variscan deformation/re-crystallization path. The established regional baseline data can be used as reference for environmental studies.     

Geostatistical Simulation of Hydrofacies Heterogeneity of the West Thessaly Aquifer Systems in Greece

Integrating geological properties, such as relative positions and proportions of different hydrofacies, is of highest importance in order to render realistic geological patterns. Sequential indicator simulation (SIS) and Plurigaussian simulation (PS) are alternative methods for conceptual and deterministic modeling for the characterization of hydrofacies distribution. In this work, we studied the spatial differentiation of hydrofacies in the alluvial aquifer system of West Thessaly basin in Greece. For this, we applied both SIS and PS techniques to an extensive set of borehole data from that basin. Histograms of model versus experimental hydrofacies proportions and indicative cross sections were plotted in order to validate the results. The PS technique was shown to be more effective in reproducing the spatial characteristics of the different hydrofacies and their distribution across the study area. In addition, the permeability differentiations reflected in the PS model are in accordance to known heterogeneities of the aquifer capacity.     

Grade and Tonnage Uncertainty Analysis of an African Copper Deposit Using Multiple-Point Geostatistics and Sequential Gaussian Simulation

Spatial uncertainty analysis is a complex and difficult task for orebody estimation in the mining industry. Conventional models (kriging and its variants) with variogram-based statistics fail to capture the spatial complexity of an orebody. Due to this, the grade and tonnage are incorrectly estimated resulting in inaccurate mine plans, which lead to costly financial decision. Multiple-point geostatistical simulation model can overcome the limitations of the conventional two-point spatial models. In this study, a multiple-point geostatistical method, namely SNESIM, was applied to generate multiple equiprobable orebody models for a copper deposit in Africa, and it helped to analyze the uncertainty of ore tonnage of the deposit. The grade uncertainty was evaluated by sequential Gaussian simulation within each equiprobable orebody models. The results were validated by reproducing the marginal distribution and two- and three-point statistics. The results show that deviations of volume of the simulated orebody models vary from ? 3 to 5% compared to the training image. The grade simulation results demonstrated that the average grades from the different simulation are varied from 3.77 to 4.92% and average grade 4.33%. The results also show that the volume and grade uncertainty model overestimates the orebody volume as compared to the conventional orebody. This study demonstrates that incorporating grade and volume uncertainty leads to significant changes in resource estimates.     

Using a Discrete Fracture Network and Spatial Point Processes to Populate Veins and Model Grade in a Coarse Gold Deposit

Vein-hosted gold deposits are characterized by mineralization, which is spatially restricted to narrow vein structures. Drillholes intersecting a mineralized vein can lead to unreliable and biased assay values compared to selective mining unit scale block grades. In this work, a discrete fracture network is simulated and adapted to model gold mineralization within the veins. Veins are assumed planar and the required inputs are distributions of vein orientation, vein length, and vein intensity (i.e., density). These inputs are collected from drillhole data, geological mapping, and expert knowledge of the deposit. A spatial point process is then applied to model gold grade as discrete events or “nuggets,” which are spatially restricted to the simulated quartz veins for the case of incomplete mineralization of the veins; when the vein is completely mineralized, a vein thickness distribution is required. The methodology is applied to an epithermal gold deposit in northwestern British Columbia, Canada and shows improvement in restricting the influence of the high-grade gold samples without resorting to ad-hoc manipulation of input assays through capping or cutting. The final output of this methodology is a block model of gold grade, which better honors the spatial structure of the veins in the deposit and is suitable for use in mine planning or resource estimation.     

Geometallurgical Approach for Implications of Ore Blending on Cyanide Leaching and Adsorption Behavior of Witwatersrand Gold Ores,South Africa

Gold production in South Africa is projected to continue its decline in future, and prospects for discovery of new high-grade deposits are limited. Many of the mining companies have resorted to mining and processing low-grade and complex gold ores. Such ores are technically challenging to process, which results in low recovery rates, excessive reagent consumption and high operating costs when compared to free-milling gold ores. In the Witwatersrand mines, options of blending low-grade gold ores with high-grade ores exist. Although it is well known that most of the Witwatersrand gold ores are highly amenable to gold cyanidation, not much is known on the leachability of blended ores, especially the effects of mineralogical and metallurgical variability between different gold ores. In this study, we apply a geometallurgical approach to investigate mineralogical and metallurgical factors that influence the leaching of blended ores in a set of bottle shaker and reactor column tests. Three gold-bearing conglomerate units, so-called reefs, i.e., Carbon Leader Reef, Ventersdorp Contact Reef and the Black Reef, all in the Carletonville goldfield, were sampled. The ores were prepared using a terminator jaw crusher followed by vertical spindle pulverizer (20 kg aliquot) and high-pressure grinding rolls (80 kg aliquot). Mineralogical analysis was conducted using a range of complementary tools such as optical microscopy, QEMSCAN and micro–XCT. The results show that Witwatersrand gold ores are amenable to the process of ore blending. Some of the ores, however, contain impervious inert gangue and reactive ore minerals. Leach solution can only access gold locked in impervious gangue minerals through HPGR-induced pores and/or cracks. The optimum ore blending ratio of the bottle shaker experiments (p80?=???75 μm) comprises 60% Carbon Leader Reef, 20% Ventersdorp Contact Reef and 20% Black Reef and yields 92% recovered Au over a leach period of 40 h. Blended ores with high carbonaceous material (>?1 wt% carbonaceous material, (Black Reef?=?36–60%) yield lower recoveries of 60–69% Au). Ore leaching at the mixed-bed reactor column (??75 μm and ??5.6/+?4 mm) yields about 70% over a leach period of two weeks. We therefore suggest that the feasibility of ore blending is strongly controlled by the mineralogy of the constituent ores and that a mixed-bed reactor may be a viable alternative method for leaching of the low-grade Witwatersrand gold ores. Material from certain reefs, such as the Black Reef, has synergistic/antagonistic (nonadditive) blending effects. The overall implication of this study is that ore blending ratios, effects of comminution on mineral liberation, an association of gold with other minerals, and gold adsorption behavior will greatly inform future technology choices in the area of geometallurgy.

     

Selective omission of road features based on mesh density for automatic map generalization

Selection of roads is an intractable generalization operation due to the difficulty in retaining the density difference and connectivity of a road network. This paper proposes a novel approach of selective omission for roads based on mesh density. The density of a road network and its local variations are calculated using meshes as units. Since maps at different scales usually reveal different densities, different density thresholds for road networks are determined on the basis of theoretical analysis and empirical study of mesh densities on maps at different scales. The selection process starts with the identification of the meshes that have a density beyond the threshold. The mesh with the largest density is first treated. Its bounding road segments are ordered according to their relative importance. The least important segment is eliminated. The remaining segments are then merged with the adjacent mesh, thus forming a new mesh. The selection procedure is repeated until none of the meshes has a density beyond the threshold. Such a process of eliminating road segments and merging meshes can ensure the road network connectivity. In this study, the meshes are classified depending on the types of road segment. For the different mesh types, their density thresholds are set to be different, which can be used as an indicator for the preservation of the density difference. This proposed approach considers topological, geometric and semantic properties of the road network. It was applied to two sets of road networks, and the results of selection are convincing. This methodology has now been adopted for the updating of 1:50,000 maps of China.     

Orebody Modelling for Exploration: The Western Mineralisation, Broken Hill, NSW

The Western Mineralisation in the Broken Hill deposit was studied to identify the zonation sequence of lithogeochemical haloes along and across the strike of the orebody. Samples used are from 77 drill holes and the samples were assayed for Pb, Zn, Fe, S, Cu, Ag, Cd, Sb, Bi and As. Variogram analyses were calculated for all the elements and kriging was used to construct the 3D block model. Analysis of cross sections along and across the strike of the orebody shows that Bi and Sb form broader halos around sulphide masses and this suggests that they are pathfinder elements for the Pb and Zn elements of this orebody. The threshold concentrations (minimum anomaly) of the 10 elements were determined using the concentration-area analysis. On east?Cwest vertical cross sections, the values of linear productivity, variability gradient and zonality index were calculated for each element. Based on the maximum zonality index of each element, the sequence of geochemical zonation pattern was determined from top to bottom of the orebody. The result shows that S, Pb, Zn and Cd tend to concentrate in the upper part of the mineralisation whereas Ag, Cu, Bi and As have a tendency to concentrate in the lower part of the mineralised rocks. Also, an empirical product ratio index was developed based on the position of the elements in the zonation sequence. The methods and results of this research are applicable to exploration of similar Zn and Pb sulphide ore deposits.     

Multi-scale Numerical Simulation and 3D Modeling for Deep Mineral Exploration in the Jiaojia Gold District,China

In this study, deposit- and district-scale three-dimensional (3D) fault-and-intrusion structure models were constructed, based on which a numerical simulation was implemented in the Jiaojia gold district, China. The numerical simulation of the models shows the basic metallogenic path and trap of the gold deposits using mineral system theory. The objective of this study was to delineate the uncertainty of the geometry or buffer zones of the ore-forming and ore-controlling fault-and-intrusion domains in 3D environment representing the exploration criteria extraction and the gold potential targeting in the study area. The fast Lagrangian analysis of continua in three dimensions was used as the platform to define the stress deformation fracture ore storage and the hydrothermal seepage channel zone based on the gold deposit features and metallogenic model in the study area. The validity of the numerical simulation was verified by comparing it with robust 3D geological models of the large Xincheng gold deposit. The potential targeting zones are analyzed for uncertainty and then evaluated by Boolean operation in a 3D geological model using the computer-aided design platform. The research results are summarized as follows. (1) In the pre-mineralization period, the Jiaodong fault’s left lateral movement created the Jiaojia network faults and formed a fracture zone with NW- to NNW-trending dips of 20° to 40°. (2) During the mineralization period, hydrothermal flow was associated with the intrusion geometry and features. However, it was constrained by the Jiaojia fault, which blocked the vadose flow into the upper wall rock and made the hydrothermal route close to the fault in the footwall fracture zones. (3) Three gold potential targets were identified by the numerical simulation results in the study area: the NW-trending Sizhuang gold deposit, the NW-trending zone of Jiaojia gold deposit, and the NE-trending zone of the Xincheng gold deposit. (4) The numerical simulation results show the fault-and-intrusion metallogenic domain and the hydrothermal alteration zones, which reflect the main ore-controlling and ore-forming factors of mineralization. The information obtained through the numerical simulation discussed here can be used to define exploration criteria in the study area.

     

Economic-Geological Estimation of Gold and Ore Deposits Using Information-Analytical System for the Selection of Analogues

The computer system is developed for evaluation of gold deposits over the world using data collected by experts in the field of gold deposit geology. With the system in use, gold object of interest can be described on the basis of selected groups of tags. The similarity of the considered object to an object from the database is estimated by designed algorithms on the basis of available information characterizing the considered object. These algorithms allow to apply mathematical methods for analysis of the problem what improves quality and reliability of the results. By the method of analogue determination, user can always evaluate unknown properties of the considered object, as for example, ore formation type, economic geological estimate of the deposit on the basis of known properties of the determined analogue.The software equipped by an intuitive graphical interface is developed with use of GIS technologies and modern methods of object-oriented programming.     

The Fractal Characteristics of Drainage Networks and Erosion Evolution Stages of Ten Kongduis in the Upper Reaches of the Yellow River,China

The fractal characteristics of drainage in the ten kongduis of the upper Yellow River were obtained using the box counting dimension, and the evolution stages of the watershed topography were defined by different ranges of the fractal dimensions of river networks (Dg). The results show that the fractal scaleless range of the Maobula River is 20–370 m based on a combination of artificial judgment, correlation coefficient test and fitting error. Other kongduis show good fractal characteristics in this fractal scaleless range as well. The box counting dimension can be used as a quantitative index of watershed topography fractal characteristics. The fractal dimension of stream networks is independent of the threshold contributing area used for extracting the drainage networks from the DEM. The values of Dg in the upper ten kongduis are in the range of 1.08?1.14. Both the runoff yield and the sediment yield are positively and linearly related with Dg. The positive relation between the sediment yield and Dg reflects the effect of landform features on sediment yield in the young and/or mature stages of landform evolution of the study area. By revising the critical value of Dg, the value of Dg of the basin in the young evolution stage is less than 1.06, while it is more than 1.06 for the basin in mature or old evolution stage. The upper ten kongduis are in the mature stage of landform evolution.     

Integration of Sparse Geologic Information in Gold Targeting Using Logistic Regression Analysis in the Hutti–Maski Schist Belt,Raichur, Karnataka,India—A Case Study

Logistic regression has been used in the study to integrate indicator patterns for estimation of the probability of occurrence of gold deposits in a part of the auriferous Archaean Hutti–Maski schist belt. Data used consist of categorical and continuous variables obtained from a coded lineament map and geochemical anomaly maps of the pathfinder elements of gold in soil and groundwater. Main effects and interactions of the variables studied were used in formulating the logistic regression model. Regression models using lineament-proximity data, combined with soil and groundwater geochemical anomalies were tested on parts of the schist belt with data not used in estimation of model parameters. Predicted probabilities greater than 0.9 identified known deposit locations in the area.     

ArcFractal: An ArcGIS Add-In for Processing Geoscience Data Using Fractal/Multifractal Models

Natural Resources Research - Fractal and multifractal models, including the concentration-area (C–A) fractal model, spectrum-area (S–A) multifractal model, and local singularity...     

Large-scale mineral potential estimation for blind precious metal ore bodies

Conventional evaluation of quantitative mineral potential has focused on target selection at small scales. Mapping at small scales usually results in large-area targets, which may be suitable for grass-roots exploration or regional evaluation of potential. Unfortunately, the estimates in small-scale exploration are commonly associated with large uncertainties. Large-scale estimation is used for optimal in-fill drilling design and step-out drilling target selection. In-fill drilling helps to confirm ore-grade continuities and translate a portion of geological resources into minable reserves, whereas step-out target estimation is useful for finding new orebodies in the vicinity of known ore deposits. Both of these processes are necessary for mine development and production planning. A comprehensive methodology is proposed here, particularly for large-scale mineral exploration. The central information synthesizer is canonical or indicator favorability analysis. A case study is presented to demonstrate the methodology for large-scale target selection. The study involves a gold-mining district where step-out drilling targets are being sought to expand the resource base. Several drilling targets were delineated in the study region. Two of them were tested through surface sampling with positive results.     

Fractal pattern integration for mineral potential estimation

Concepts of fractal/multifractal dimensions and fractal measure were used to derive the prior and posterior probabilities that a small unit cell on a geological map contains one or more mineral deposits. This has led to a new version of the weights of evidence technique which is proposed for integrating spatial datasets that exhibit nonfractal and fractal patterns to predict mineral potential. The method is demonstrated with a case study of gold mineral potential estimation in the Iskut River area, northwestern British Columbia. Several geological, geophysical, and geochemical patterns (Paleozoic-Mesozoic sedimentary and volcanic clastic rocks; buffer zones around the contacts between sedimentary rocks and Mesozoic intrusive rocks; a linear magnetic anomaly; and geochemical anomalies for Au and associated elements in stream sediments) were integrated with the gold mineral occurrences which have fractal and multifractal properties with a box-counting dimension of 1.335±0.077 and cluster dimension of 1.219±0.037.     

引发意大利托斯卡纳区滑坡的区域降水阈值体系的定义——统计与环境分析(英文)

The aim of this work is the determination of regional-scale rainfall thresholds for the triggering of landslides in the Tuscany Region(Italy).The critical rainfall events related to the occurrence of 593 past landslides were characterized in terms of duration(D) and intensity(I).I and D values were plotted in a log-log diagram and a lower boundary was clearly noticeable:it was interpreted as a threshold representing the rainfall conditions associated to landsliding.That was also confirmed by a comparison with many literature thresholds,but at the same time it was clear that a similar threshold would be affected by a too large approximation to be effectively used for a regional warning system.Therefore,further analyses were performed differentiating the events on the basis of seasonality,magnitude,location,land use and lithology.None of these criteria led to discriminate among all the events different groups to be characterized by a specific and more effective threshold.This outcome could be interpreted as the demonstration that at regional scale the best results are obtained by the simplest ap-proach,in our case an empirical black box model which accounts only for two rainfall pa-rameters(I and D).So a set of thresholds could be conveniently defined using a statistical approach:four thresholds corresponding to four severity levels were defined by means of the prediction interval technique and we developed a prototype warning system based on rainfall recordings or weather forecasts.     

Estimating Organic-Rich Shale Fractal Dimensions from Gas Adsorption Isotherms: Combining Different Methods Leads to More Reliable Values and Insight

Fabric and roughness of the pore-size distributions in organic-rich shales determine their fluid flow and storage capabilities. Accurate estimation of the pore-surface fractal dimension (D) provides valuable insight to these qualities in shales. Low-pressure gas adsorption isotherms are widely used for determining D, typically applying the Frenkel–Halsey–Hill (FHH) method. Other D estimation methods, proposed by Neimark (NM) and Wang and Li (WL), are theoretically consistent and mathematically related to the FHH model but yield distinctive D values for many shales. This study evaluates the mathematical relationships between the FHH, NM and WL fractal determination methods, and with the aid of twenty-six published adsorption isotherms from shales around the world, compares their similarities and differences. Uncertainties exist in establishing best-fit lines to curved data trends in the FHH and NM methods, and in fitting power curves to data trends in the NM and WL methods. The FHH and WL D values are found to be more consistent for whole isotherm and isotherm segment analysis than the NM D values, which are systematically higher. The reasons for this are explained in terms of their graphical relationships. This leads to a novel 10-step protocol for a more thorough determination of shale D values that incorporates all three methods and involves graphical analysis that clearly exposes the uncertainties associated with the values determined. Applying this protocol should derive reliable D values to compare with key shale properties such as surface area, surface volume, thermal maturity and organic richness in future research.