Geostatistical modeling of the geological uncertainty in an iron ore deposit

This paper addresses the problem of quantifying the joint uncertainty in the grades of elements of interest (iron, silica, manganese, phosphorus and alumina), loss on ignition, granulometry and rock types in an iron ore deposit. Sampling information is available from a set of exploration drill holes. The methodology considers the construction of multiple rock type outcomes by plurigaussian simulation, then outcomes of the quantitative variables (grades, loss on ignition and granulometry) are constructed by multigaussian joint simulation, accounting for geological domains specific to each quantitative variable as well as for a stoichiometric closure formula linking these variables. The outcomes are validated by checking the reproduction of the data distributions and of the data values at the drill hole locations, and their ability to measure the uncertainty at unsampled locations is assessed by leave-one-out cross validation.Both the plurigaussian and multigaussian models offer much flexibility to the practitioner to face up to the complexity of the variables being modeled, in particular: (1) the contact relationships between rock types, (2) the geological controls exerted by the rock types over the quantitative variables, and (3) the cross-correlations and stoichiometric closure linking the quantitative variables. In addition to this flexibility, the use of efficient simulation algorithms turns out to be essential for a successful application, due to the high number of variables, data and locations targeted for simulation.     

Risk quantification with combined use of lithological and grade simulations: Application to a porphyry copper deposit

The uncertainty in the recoverable tonnages and grades in a mineral deposit is a key factor in the decision-making process of a mining project. Currently, the most prevalent approach to model the uncertainty in the spatial distribution of mineral grades is to divide the deposit into domains based on geological interpretation and to predict the grades within each domain separately. This approach defines just one interpretation of the geological domain layout and does not offer any measure of the uncertainty in the position of the domain boundaries and in the mineral grades. This uncertainty can be evaluated by use of geostatistical simulation methods. The aim of this study is to evaluate how the simulation of rock type domains and grades affects the resources model of Sungun porphyry copper deposit, northwestern Iran. Specifically, three main rock type domains (porphyry, skarn and late-injected dykes) that control the copper grade distribution are simulated over the region of interest using the plurigaussian model. The copper grades are then simulated in cascade, generating one grade realization for each rock type realization. The simulated grades are finally compared to those obtained using traditional approaches against production data.     

A probabilistic parametrization for geological uncertainty estimation using the ensemble Kalman filter (EnKF)

In the past years, many applications of history-matching methods in general and ensemble Kalman filter in particular have been proposed, especially in order to estimate fields that provide uncertainty in the stochastic process defined by the dynamical system of hydrocarbon recovery. Such fields can be permeability fields or porosity fields, but can also fields defined by the rock type (facies fields). The estimation of the boundaries of the geologic facies with ensemble Kalman filter (EnKF) was made, in different papers, with the aid of Gaussian random fields, which were truncated using various schemes and introduced in a history-matching process. In this paper, we estimate, in the frame of the EnKF process, the locations of three facies types that occur into a reservoir domain, with the property that each two could have a contact. The geological simulation model is a form of the general truncated plurigaussian method. The difference with other approaches consists in how the truncation scheme is introduced and in the observation operator of the facies types at the well locations. The projection from the continuous space of the Gaussian fields into the discrete space of the facies fields is realized through in an intermediary space (space with probabilities). This space connects the observation operator of the facies types at the well locations with the geological simulation model. We will test the model using a 2D reservoir which is connected with the EnKF method as a data assimilation technique. We will use different geostatistical properties for the Gaussian fields and different levels of the uncertainty introduced in the model parameters and also in the construction of the Gaussian fields.     

Simulation of mineral grades and classification of mineral resources by using hard and soft conditioning data: application to Sungun porphyry copper deposit

The presence of geological units with different grade characteristics mostly leads to problems during the grade modeling process. In special cases, if the area under study has units with small thickness and low grade with respect to the dominant unit of the area, it is difficult to reproduce different grade contents in these units in the simulated grade models because of the low thickness and lack of data in these units. In this study, the local moment constraints method, based on the definition of soft conditioning data reflecting geological knowledge, is investigated for improving simulated grade models under the mentioned conditions. This method is applied for grade simulation at the 1,750 m level of Sungun porphyry copper mine. The studied area is divided into two rock type domains: Sungun porphyry and Dyke. The Sungun porphyry unit is the dominant rock type in the considered area and has, on average, a higher copper grade, while dykes discontinue Sungun porphyry rock units sporadically and most of them are barren of mineralization. It is demonstrated that the use of soft conditioning data makes the simulated grade model closer to reality and improves the reproduction of grade contents considering the rock type units in the area. In the next step, the results obtained from conditional simulation are used for mineral resources classification. To this end, the conditional coefficient of variation is chosen as a criterion for measuring uncertainty and for defining the resources classes. Then, it is shown that uncertainty can be considerably reduced in the prepared models if soft data are considered; as a result, an increase in measured resource classification is observed.     

Data Assimilation in Truncated Plurigaussian Models: Impact of the Truncation Map

Assimilation of production data into reservoir models for which the distribution of porosity and permeability is largely controlled by facies has become increasingly common. When the locations of the facies bodies must be conditioned to observations, the truncated plurigaussian model has been often shown to be a useful method for modeling as it allows gaussian variables to be updated instead of facies types. Previous experience has also shown that ensemble Kalman filter-like methods are particularly effective for assimilation of data into truncated plurigaussian models. In this paper, some limitations are shown of the ensemble-based or gradient-based methods when applied to truncated plurigaussian models of a certain type that is likely to occur for modeling channel facies. It is also shown that it is possible to improve the data match and increase the ensemble spread by modifying the updating step using an approximate derivative of the truncation map.     

酒西盆地下沟组湖相白烟型喷流岩岩石学特征

酒西盆地下白垩统下沟组首次发现的湖相白烟型喷流岩,以其复杂的矿物组分和丰富的沉积构造业已引起众多地质工作者的重视。本文通过深入系统的岩石学研究,按物质组分的组合特征,从此类喷流岩中划分出6种热水矿物组合,分别为钠长石-铁白云石或铁白云石-钠长石组合、石英-钠长石-铁白云石组合、石英-重晶石-钠长石-铁白云石组合、石英-方沸石-钠长石-铁白云石组合、地开石-铁白云石组合和单一铁白云石组合,根据热水矿物组合和结构构造特征,可将下沟组喷流岩划分成4种成因类型,分别为脉状充填型、水爆角砾型、盆地沉积型和区域扩散型,不同成因类型的喷流岩具有不同的热水矿物组合特征和热水沉积地质意义,在此基础上,建立了下沟组湖相白烟型喷流岩热水沉积模式,     

Incorporating Information About Edge Effects When Simulating Lithofacies

Some years ago, earth scientists came to realize that knowing more about the geology of an orebody or an oil reservoir makes it easier to make the appropriate decisions concerning mine planning or reservoir exploitation. Geostatistical techniques for simulating lithofacies—that is, the geometry of the geology—were developed as a result of this. These methods should be able to produce geological images that respect not only the anisotropies of the different lithofacies but also their spatial layout relative to one another. While indicator variograms ensure that anisotropies are respected, another tool needs to be incorporated in the simulation technique to reflect the relative spatial layout of the different lithofacies. We propose to use the concept of edge effects that define the position of one lithofacies relative to another. Simple tests using direct and cross indicator variograms confirm the presence or absence of edge effects. We investigate if and how edge effect information can be incorporated in the different indicator simulation techniques—sequential indicator simulations, simulated annealing, the truncated Gaussian method and plurigaussian simulations. Results show that the choice of simulation method must be guided by the edge effect characteristics of the experimental lithologic data.     

拉伸剪切条件下岩石的工程地质力学特性

拉伸剪切破坏是岩石边坡工程与岩石隧道工程中常出现的一种岩石失稳类型。结合岩石工程地质条件,分析了拉伸剪切应力状态产生的条件。总结了当前岩石拉伸剪切研究的实验手段,本构关系与破坏准则。指出了岩石拉伸剪切存在的不足。提出如下研究趋势:(1)研制标准的岩石直接拉伸剪切设备,进行专门试验研究; (2)进行更多类型岩石,更多数量的岩石拉伸剪切实验。研究不同岩石类型拉伸剪切本构关系与破坏准则; (3)加强岩石拉伸剪切实验与此类型岩石工程的数值模拟工作,与试验相互印证。     

地层中的硅岩: 研究进展及展望*

硅岩是由自生石英形成的化学沉积岩,通常以层状或非层状(结核状、条带状)独自出现、与其他岩性互层出现及赋存于其他岩性地层中。硅岩不仅记录了大量盆地古地理、古海洋、古气候的地质信息,而且对于岩石大地构造学研究具有重要意义。前人对硅岩的研究取得了丰富的成果,近年来关于硅岩的新观点为硅岩的研究提供了新的思考途径和研究线索。通过调研近年来关于硅岩的岩石学特征、地球化学特征及沉积环境的研究,分析总结了不同地质历史时期、不同沉积环境、不同产状硅岩的形成模式、成因类型以及判别方式,认为硅岩的沉积不是单一因素控制的结果,而是地质构造活动强度、生物类型及活动性、沉积古环境的综合结果。常规地球化学元素结合岩石组构与微区分析、包裹体分析、新同位素(Rb-Sr和Sm-Nd)分析,以及新的地球化学判定标准(Ge/Si)是硅岩研究的必备手段和方法。     

地质断层三维构模技术研究

对包含断层的复杂地质体三维构模技术展开讨论并提出相应的解决方案。分析了实际勘察所获取的断层数据类型,设计了标准的断层数据库格式;系统地评价了两类基本的断层构模方法（即基于地层恢复的断层构模技术和基于分区插值的断层构模技术）,分析了各自的特点和不足,在此基础上提出了断层与地层的统一构模技术;该技术既能应用于同沉积断层,又能适应断层终止于地层内部的情况,还能处理具有多值面的逆断层构造的拟合问题,在理论上能够构造任意复杂的断层系统。通过北京奥运公园场址的断层构模实例验证了该方法的有效性。     

Transitive Geostatistics for Stepwise Modeling Across Boundaries between Rock Regions

A single intrinsic stationary random field may not account for transitional heterogeneity and abrupt dissimilarity of geological properties across boundaries between rock type regions. This paper proposes the stepwise construction of transitive covariance models for modeling continuous properties correlated across boundaries of multiple disjoint physical domains such as rock type bodies. Modeling in geology is usually simplified by splitting the geological space into rock type geo-domains (e.g., strata, sedimentary facies, soil series, diagenetic regions and alteration zones). Due to the limitations of simultaneous solutions, a simplification is to model each domain independently at the cost of losing the conditioning of properties across domains. This paper proposes to organize the modeling process in a triangular array which follows events in the geological time domain; for example, the younger formations are at the top of the pyramid and the older formation at the base. The estimation may go from top to base by assuming that younger events have perturbed older formations. Geology shows the scars of events that cumulate in rock formations before they are finally eroded. In some cases, older formations may be parent material for younger formations. The continuous property within each geo-domain has a conditional covariance in the main diagonal of the array which may belong to a specific event in the geological time. This sequence leads to transitive estimation and simulations in the physical space. If a simultaneous solution is sought (i.e., the future and past are correlated both ways), the complex covariance functions can be constructed stepwise from conditional spectra.     

基于多源数据约束的成矿动力学模拟——以宁芜盆地钟姑矿田典型矿床为例

数据驱动是大数据分析的特点之一。成矿过程的数值模拟即采用数据驱动方式,在基本的物理化学方程及地质模型基础上,通过定量的手段,模拟不同阶段、不同条件下成矿系统的演化过程,从而能够更深入、直观的讨论各种地质要素对成矿系统的影响,验证已有的成矿理论,获得新的认识。姑山铁矿床和杨庄铁矿床是宁芜火山岩盆地钟姑矿田内的典型矿床。前人研究显示,杨庄铁矿床为典型的充填-接触交代型热液矿床,而姑山铁矿床目前在成因上仍存在争议,存在岩浆成因和热液成因等多种观点。为了更好的对比和了解姑山铁矿床和杨庄铁矿床的成因及成矿过程,本文基于多源数据的综合解译剖面,采用力-热耦合数值模拟方法对与热液矿床密切相关的扩容空间的形成机制和特点开展研究。模拟结果显示,姑山铁矿床和杨庄铁矿床可能分属不同的成因类型;随着岩体的逐渐冷却,杨庄铁矿床内岩体和地层接触带会出现明显的扩容空间,可为成矿流体的汇聚及矿质的沉淀提供良好的条件,有利于热液矿床的形成;而姑山铁矿床内岩体和地层接触带只有局部区域出现扩容空间,同时体应变数值较小且形态变化大,并不利于大规模的流体汇聚及矿质的沉淀,姑山铁矿床可能主要为岩浆成因,并受少量热液的叠加改造。     

辽宁省抚顺矿区宏观地质结构及其工程效应分析

抚顺采煤已有百余年历史,由于各种原因造成了一系列地质灾害:滑坡、矿震、地表沉陷和地裂缝等。矿区的中部和西部的形成带状的地表拉裂和建筑物破坏;软岩区形成大面积的地表沉陷,硬岩区形成脆性和线性拉裂,这些现象与矿区地质结构是分不开的。通过对矿区地质结构和矿区地层构造演化研究,认为抚顺矿区地质灾害是发生在花岗岩与火山岩之间的"凹槽"内,是有限的。根据矿区主断层F1、F1A交汇点位置,将矿区地质结构类型划分为三种类型,并通过FLAC二维数值计算对三种类型的典型剖面进行了数值模拟和验证,合理地解释了矿区地质灾害的产生机理和发展趋势,为抚顺矿区地质灾害综合治理和采煤规划提供了科学依据。采用地质结构控制论来研究城市地质灾害提供了新的思路和方法。     

江西银山矿田伴生金银综合回收利用研究

叙述了银山铜铅锌多金属矿田铅锌成矿系列和铜硫成矿系列６个成矿带中伴生金银的赋存状态、富集特征和成矿规律，并进行伴生金银相关性分析。认为伴生金银是以独立矿物存在为主，主要载体矿物有黄铜矿，含砷铜矿物、黄铁矿、方铅矿和闪锌矿。Ａｕ、Ａｇ与Ｐｂ、Ｚｎ、Ｃｕ、Ｓ密切相关。同时分析了金银生产中存在的具体问题，提出综合回收金银的措施和建议。     

中国典型钴矿床地质特征及重点地区矿产资源预测

我国既是钴矿资源消费大国,又是钴矿资源进口大国,受新能源电动汽车工业的影响,近年来钴矿资源得到越来越多的关注。为了科学地评估我国钴矿资源潜力,指导找矿勘查部署工作,立足国内,提高我国钴矿资源的保障能力,笔者2019—2021年开展了9个重点省(自治区)的钴矿资源潜力评价工作。本文从我国钴矿床分布和地质特征角度,先以成矿地质作用为主线,厘定我国找矿预测矿床类型,包括与沉积地质作用有关的风化型、化学沉积型、砂岩型和海底喷流沉积型,与火山地质作用有关的海相火山岩型、陆相火山岩型,与侵入岩浆地质作用有关的岩浆型、矽卡岩型(接触交代型)、热液脉型和斑岩型,与变质地质作用有关的沉积变质型,以及少量与大型变形地质作用有关的变质核杂岩型;进而梳理并总结了各类型钴矿有关的成矿地质体。通过分析我国钴矿的矿床地质特征和时空分布规律,综合地质、物探、化探和遥感信息等预测要素,建立找矿预测模型,并在此基础上圈定找矿远景区,估算资源量。结果表明：我国钴矿包括风化型镍钴矿、海底喷流沉积型铜钴矿、海相火山岩型铁铜钴矿、海相火山岩型块状硫化物铜锌钴矿、岩浆型铜镍钴矿和热液脉型钴矿6种主要找矿预测类型;圈定416个钴矿最...     

三峡库区典型塌岸模式研究

三峡库区塌岸地质条件调查结果表明,不同的岸坡结构,蓄水后库岸的变形破坏形式塌岸模式,往往差别较大。本文根据现场调查成果,结合室内分析,总结归纳出三峡库区几种典型的塌岸模式,即冲蚀磨蚀型、坍塌型、崩塌(落)型、、滑移型和流土型。在此基础上,对不同塌岸模式的发展演进过程、发育分布特征,以及塌岸与岸坡岩土结构、地形地貌等的关系进行了深入系统地研究,提出了三峡库区典型塌岸发育的一般地质条件。     

浙江省新昌县下山村滑坡地质特征与稳定性评价

本文研究了浙江省新昌县下山村滑坡区地质环境条件,介绍了滑坡地质灾害体的地质结构、形成特征。认为该区地质环境复杂,雨量充沛,地形呈阶梯状;岩层以5个旋回的软硬相间为主,覆盖层为第四纪松散层;滑坡体赋存松散岩类孔隙水和基岩孔隙-裂隙水,滑坡区岩土体遇水后工程地质性质变差,其抗滑、抗软化能力较差。滑坡表现为后缘切层、中-前缘顺层型。文章从滑坡成因类型的形成条件、变形破坏机制和诱发滑坡体失稳的主要因素方面,对滑坡的稳定性作出了基本评价,对其发展变化作出了趋势性预测。认为该滑坡体物质结构复杂,分7（大）层,滑坡有多个滑移带,主滑移面（带）2个,滑坡体厚达60m,降雨是诱发滑坡体失稳的主要原因,该滑坡目前仍以蠕滑为主。     

Petrographical,structural, and mechanical analysis of armorican phtanites: A key raw material for the Mesolithic in Western France

Over the last ninety years or so, archaeological field work on outcrops of raw materials in western France has allowed researchers to identify the use of a large number of deformed rock types in lithic industries. These materials include Armorican phtanites, which are cherts that were widely exploited during the Mesolithic, principally during its final phase (5500–5000 yrs cal. B.C.). Because of variations in the mineralogical composition and physical properties of phtanites, different mechanical properties are exhibited by samples from the four main geological formations containing this rock type in the Armorican Massif of western France. These geological and mechanical parameters influenced the quantity of rock used and the knapping methods employed. They also influenced the tool types produced and their transfer, as recorded in archaeological sequences. The types of raw materials used are also controlled by the distance to sources of coastal flint pebbles or other autochthonous siliceous rocks. Therefore, the role of raw material properties in the creation and persistence of techno‐cultural or economic characteristics is different for each type of phtanite and yields fundamental information on the interaction of society with the geological environment. © 2010 Wiley Periodicals, Inc.     

内蒙古乌努格吐山铜钼矿露天采场边坡工程地质特征及破坏类型

通过对乌努格吐山铜钼矿露天采场边坡自然背景条件、岩体结构类型、岩石力学参数、工程地质岩组特征、结构面特征、边坡工程地质分区、边坡稳定性评价、边坡破坏类型及规模等论述,阐明了露天采场边坡工程地质特征.边坡稳定性计算结果表明,A区A-A'剖面、C区E-E'剖面稳定性计算结果不能满足安全系数的要求,其他区边坡总体稳定性较好.露天采场现处于生产建设阶段,边坡正处于剥离形成过程中,边坡几何形态保存较完整,仅局部出现小型崩塌、垮塌、粒状解体、楔形破坏、地裂缝和大气降水冲刷等.尚未出现较大规模破坏,现状边坡总体稳定,基本可代表未来边坡的破坏类型和稳定状态.     

Deep learning based classification of rock structure of tunnel face

The automated interpretation of rock structure can improve the efficiency,accuracy,and consistency of the geological risk assessment of tunnel face.Because of the high uncertainties in the geological images as a result of different regional rock types,as well as in-situ conditions(e.g.,temperature,humidity,and construction procedure),previous automated methods have limited performance in classification of rock structure of tunnel face during construction.This paper presents a framework for classifying multiple rock structures based on the geological images of tunnel face using convolutional neural networks(CNN),namely Inception-ResNet-V2(IRV2).A prototype recognition system is implemented to classify 5 types of rock structures including mosaic,granular,layered,block,and fragmentation structures.The proposed IRV2 network is trained by over 35,000 out of 42,400 images extracted from over 150 sections of tunnel faces and tested by the remaining 7400 images.Furthermore,different hyperparameters of the CNN model are introduced to optimize the most efficient algorithm parameter.Among all the discussed models,i.e.,ResNet-50,ResNet-101,and Inception-v4,Inception-ResNet-V2 exhibits the best performance in terms of various indicators,such as precision,recall,F-score,and testing time per image.Meanwhile,the model trained by a large database can obtain the object features more comprehensively,leading to higher accuracy.Compared with the original image classification method,the sub-image method is closer to the reality considering both the accuracy and the perspective of error divergence.The experimental results reveal that the proposed method is optimal and efficient for automated classification of rock structure using the geological images of the tunnel face.