Towards Stochastic Time-Varying Geological Modeling

The modeling of subsurface geometry and properties is a key element to understand Earth processes and manage natural hazards and resources. In this paper, we suggest this field should evolve beyond pure data fitting approaches by integrating geological concepts to constrain interpretations or test their consistency. This process necessarily calls for adding the time dimension to 3D modeling, both at the geological and human time scales. Also, instead of striving for one single best model, it is appropriate to generate several possible subsurface models in order to convey a quantitative sense of uncertainty. Depending on the modeling objective (e.g., quantification of natural resources, production forecast), this population of models can be ranked. Inverse theory then provides a framework to validate (or rather invalidate) models which are not compatible with certain types of observations. We review recent methods to better achieve both stochastic and time-varying geomodeling and advocate that the application of inversion should rely not only on random field models, but also on geological concepts and parameters.     

Spatial Association of Mineral Occurrences and Curvilinear Geological Features

Qualitative and quantitative knowledge about the spatial association between mineral occurrences and geological features are important in mineral potential mapping. Two existing methods for quantifying spatial association between mineral occurrences and curvilinear geological features are applied to the Baguio district of the Philippines. An experimental method is described and applied to the study area as well. The results of the three methods are highly similar, which demonstrates the effectiveness of the experimental method presented here for quantifying spatial association between mineral occurrences and curvilinear geological features. It is shown that gold occurrence in the Baguio district are strongly spatially associated with northeasterly trending faults/fractures rather than with northwesterly trending faults/fractures. It is also shown that the spatial association between the gold occurrence and older batholithic intrusives is stronger than the spatial association between gold occurrence and younger porphyry intrusives. These spatial geoinformation characteristics can be used as evidential data layers in GIS-based mineral potential mapping.     

Stochastic Structural Modeling

A consistent stochastic model for faults and horizons is described. The faults are represented as a parametric invertible deformation operator. The faults may truncate each other. The horizons are modeled as correlated Gaussian fields and are represented in a grid. Petrophysical variables may be modeled in a reservoir before faulting in order to describe the juxtaposition effect of the faulting. It is possible to condition the realization on petrophysics, horizons, and fault plane observations in wells in addition to seismic data. The transmissibility in the fault plane may also be included in the model. Four different methods to integrate the fault and horizon models in a common model is described. The method is illustrated on an example from a real petroleum field with 18 interpreted faults that are handled stochastically.     

Non-stationary Geostatistical Modeling Based on Distance Weighted Statistics and Distributions

A common assumption in geostatistics is that the underlying joint distribution of possible values of a geological attribute at different locations is stationary within a homogeneous domain. This joint distribution is commonly modeled as multi-Gaussian, with correlations defined by a stationary covariance function. This results in attribute maps that fail to reproduce local changes in the mean, in the variance and, particularly, in the spatial continuity. The proposed alternative is to build local distributions, variograms, and correlograms. These are inferred by weighting the samples depending on their distance to selected locations. The local distributions are locally transformed into Gaussian distributions embedding information on the local histogram. The distance weighted experimental variograms and correlograms are able to adapt to local changes in the direction and range of spatial continuity. The automatically fitted local variogram models and the local Gaussian transformation parameters are used in spatial estimation algorithms assuming local stationarity. The resulting maps are rich in nonstationary spatial features. The proposed process implies a higher computational effort than traditional stationary techniques, but if data availability allows for a reliable inference of the local distributions and statistics, a higher accuracy of estimates can be achieved.     

基于钻孔数据的滑坡三维地质建模研究

对现有建模方法进行分析,结合千将坪滑坡地质特点,首先进行钻孔数据处理,然后提出引入虚拟钻孔,再利用GOCAD软件提供的Kriging插值技术建立滑坡界面和实体模型.该方法不但为滑坡三维地质建模提供了参考,也为区域性、地质状况复杂、钻孔资料欠缺地区的三维地质建模提供了借鉴.     

一种基于物性数据的深部三维地质建模方法

地质采样信息不足是制约深部三维地质建模的重要因素, 深部物性探测数据由于其易于获取而能够有效形成可视化模型.结合这一特点, 在地质调查项目工作中探索出了一种基于物性探测数据的三维地质建模方法.它首先利用岩石样品物性测量实验数据提取出物性参数及其对应地质属性的映射关系, 然后将不同地球物理方法所获取到的物性数据进行综合建模并解释, 最后将解释后的可视化模型转换为地质三维模型.实践证明, 该方法能够针对性地解决项目中的一些深部地质三维建模问题.      

曲线坐标系下的完全匹配层吸收边界条件

在地震波数值模拟中,需要采用吸收边界条件以吸收人为边界反射。本文针对曲线坐标系下的二阶弹性波方程提出了一种完全匹配层(PML)吸收边界条件。与直角坐标系下的PML吸收边界条件类似,曲线坐标系下的PML吸收边界条件是一种在频率域中给出的人工边界条件,由相应的复坐标变换得到。在变换到时间域后,完全匹配层中将出现复杂的卷积运算。为了避免这些卷积运算,引入了4个中间变量。为了简化自由边界条件,采用正交贴体网格对起伏地表模型进行网格剖分。数值算例表明,该方法可以有效消除人为边界反射。     

基于线框架模型的三维复杂地质体建模方法

提出了一种基于线框架模型的复杂地质体三维建模方法,在准确表达复杂地质体各要素空间几何形态与相互关系的基础上,实现地层交错情况下断层面、地层面模型及复杂地质体模型的准确快速构建.建模实践表明:基于线框架模型的复杂地质体三维建模方法是准确高效的,能适应常见复杂地质条件下三维复杂地质体模型的构建.     

基于三维激光扫描技术的矿山地质建模与应用

激光扫描技术(light detection and ranging technology)作为一门新兴的测绘技术, 目前已应用于诸多领域, 然而对于海量测站数据的拼接、大地坐标的定位等方面还存在一些不足.针对这2个问题, 提出分部式拼接方法以解决大数据量难处理的问题, 并利用罗德里格矩阵算法, 编程实现对地物的绝对坐标定位, 完成矿山的地质建模, 是对矿山地质建模新方法的探讨.通过与已知坐标数据的对比, 表明该方法的误差在矿山建模的允许值范围内, 具有较好的适用性.并以内蒙古白音诺尔铅锌矿地质建模为例, 可为覆盖区矿产综合预测提供3D矿床模型及预测要素空间形态特征等方面的参考.      

基于线框单元体的三维闭合地质块体构建方法

为能充分利用空间数据模型中的拓扑关系来实现三维闭合块体的构建,引入了表示简单块体轮廓的线框单元体概念来组织三维模型间各要素之间的拓扑关系;并以此为基础给出了以方向边和方向三角形作为基本识别单元进行三维简单形体的自动识别方法。应用实例表明：利用该方法不仅可以准确构建复杂地质块体,还可以准确描述块体内各要素的拓扑关系,再现局部特殊地质现象,如小断层（悬面）等。为三维地质模型在地震正演模拟、射线追踪等分析应用方面提供了基础。     

基于地质与重磁数据集成的栾川钼多金属矿区三维地质建模

三维地质建模是当前国内外矿产资源勘查与评价研究的重要途径和热点.以栾川钼多金属矿区为例,针对研究区地质背景复杂、成矿地质条件多样性特征,开展了基于地质(岩性地层及其物性特征、地质体产状数据)与重磁正反演集成技术组合的三维地质模型建模研究,其研究内容为:(1)根据研究区地质背景和成矿条件,确定研究区矿床地质特征以及近地表地质体的空间属性特征;(2)利用重磁数据,在定性解译地质信息的基础上,开展二度半正演方法定量解译研究,并运用三维位场概率成像方法反演和解译深部地质体;(3)根据研究区地质与重磁反演解译资料综合分析,确定晚侏罗世斑岩岩体和控矿构造的时-空-因关系.研究结果表明,基于地质岩性属性特征的重磁正反演联合解译,能够综合利用地质体之间的物性差异、成因关系界定其几何形态,从而在三维环境中建立它们的时空关系,即三维地质模型,为研究区深部找矿、矿体定位预测以及金属矿产资源定量评价提供技术支撑和新途径      

基于钻孔数据的矿山三维地质建模及可视化过程研究

矿山三维地质建模与可视化技术研究,是"数字矿山"的核心组成部分,是现代矿山信息化研究的热点和重点.文章从三维地质建模的概念出发,探讨了基于钻孔数据的矿山三维地质建模及可视化的一般过程,包括地质数据库的建立,实体模型、品位块体模型的创建等,重点探讨了建模过程中地质数据库的设计,地质解译、矿体储量估算的方法与技术,并借助OpenGL三维图形开发库,在VC++6.0编程环境下,实现了三维地质的建模及可视化,为我国矿山三维地质建模及可视化软件的设计提供参考,促进了我国矿山信息化的建设.     

Stochastic Simulation of Patterns Using Distance-Based Pattern Modeling

The advent of multiple-point geostatistics (MPS) gave rise to the integration of complex subsurface geological structures and features into the model by the concept of training images. Initial algorithms generate geologically realistic realizations by using these training images to obtain conditional probabilities needed in a stochastic simulation framework. More recent pattern-based geostatistical algorithms attempt to improve the accuracy of the training image pattern reproduction. In these approaches, the training image is used to construct a pattern database. Consequently, sequential simulation will be carried out by selecting a pattern from the database and pasting it onto the simulation grid. One of the shortcomings of the present algorithms is the lack of a unifying framework for classifying and modeling the patterns from the training image. In this paper, an entirely different approach will be taken toward geostatistical modeling. A novel, principled and unified technique for pattern analysis and generation that ensures computational efficiency and enables a straightforward incorporation of domain knowledge will be presented.     

Mapping Curvilinear Structures with Local Anisotropy Kriging

Interpolating geo-data with curvilinear structures using geostatistics is often disappointing. Channels, for example, become disconnected sets of lakes when interpolated from point data. In order to improve the interpolation of geological structures (e.g., curvilinear structures), we present a new form of kriging, local anisotropy kriging (LAK). Local anisotropy kriging combines a gradient algorithm from image analysis with kriging in an iterative way. After an initial standard kriging interpolation, the gradient algorithm determines the local anisotropy for each cell in the grid using a search area around the cell. Subsequently, kriging is carried out with the spatially varying anisotropy. The anisotropy calculation and subsequent kriging steps will then succeed until the result is satisfactory in the way of reproducing the curvilinear structures. Depending on the size of the search area more or less detail in the geological structures can be reproduced with LAK. Using test examples we show that LAK interpolates data with curvilinear structures more realistically than standard kriging. In a real world case, using bathymetric data of the Oosterschelde estuary, LAK also proves to be quantitatively superior to standard kriging. Absolute interpolation errors are decreased by 23%. Local anisotropy kriging only uses information from point data, which makes the method very objective, it only presents “what the data can tell.”     

Surface-Based 3D Modeling of Geological Structures

Building a 3D geological model from field and subsurface data is a typical task in geological studies involving natural resource evaluation and hazard assessment. However, there is quite often a gap between research papers presenting case studies or specific innovations in 3D modeling and the objectives of a typical class in 3D structural modeling, as more and more is implemented at universities. In this paper, we present general procedures and guidelines to effectively build a structural model made of faults and horizons from typical sparse data. Then we describe a typical 3D structural modeling workflow based on triangulated surfaces. Our goal is not to replace software user guides, but to provide key concepts, principles, and procedures to be applied during geomodeling tasks, with a specific focus on quality control.     

四川甘孜隆措沟富银铅锌矿地质特征及找矿意义

隆措沟富银铅锌矿位于川西义敦岛弧造山带北段火山—沉积盆地内。矿体产出受层间破碎带和断裂控制,容矿岩石为上三叠统拉纳山组浅变质岩。金属矿物主要有方铅矿、闪锌矿、黄铁矿,少量砷黝铜矿、磁黄铁矿等。其矿化特征与呷村大型银多金属矿床相似,属同一成矿系列。该富银铅锌矿的发现与勘查,进一步丰富了该区成矿系列的矿床类型,为银多金属矿的找矿拓宽了思路,扩大了找矿范围和远景。     

绵竹市“天池冰川遗迹”地质特征

天地冰川位于四川绵竹市天池乡 ,属于第四纪中早期的冰川遣迹 ,距今 1～ 2Ｍａ。天池冰川位于龙门山褶断带中段太平推覆构造体的天池向斜中 ,天池向斜呈纺缍形 ,轴的走向北东 ,两翼由二叠系地层组成 ,中部是三叠系下统飞仙关组地层。北西翼地层倒转倾向北西、倾角 5 0°～ 70° ,南东翼倾向北西、倾角 4 5°～ 6 0° ,向斜的中西部褶皱呈斜卧的∑形 ,天池冰川受天地向斜中的高家坪向斜控制 ,冰川长近 5ｋｍ ,宽 10 0～ 10 0 0ｍ ,走向北东 4 0° ,根据它的结构特征可分为三段 ,即冰川源头段 ,中部冰碛丘陵段 ,前缘冰舌、冰槽谷段。冰川源头…     

Stochastic Modeling of Progressive Failure in Heterogeneous Soil Slope

Accurate and efficient simulation of the slope stability, certainly gives the reliable approach for the hazard preparedness. However, inherent difficulties associated with the slope failure make the situation very difficult. First, the computation of fracture path based on the strict mathematical formalism is cumbersome especially for the progressive failure phenomenon. Second, due to the lack of reliable knowledge on material properties, stochastic simulation is necessary. Stochastic modeling together with the progressive failure phenomenon is apparently difficult task. In this paper, heterogeneous slope is simulated considering the progressive failure phenomenon using the modified finite element method introducing the fracture along the edge of the meshes thus incorporating the time evolution of the fracture surface. This is achieved by restructuring the tessellation in every fracture stage. Unlike remeshing this technique only increases the number of total nodes while number of meshes remaining same, which is simple and natural. Further, in this research, effect of uncertainty in the material properties upon the uncertainty in the response of soil slope has been evaluated, hence giving the reliable probabilistic estimation of the factor of safety, failure surface and deformation of the slope. Thus the simulation based on the more reliable understanding of the material properties and the efficient numerical procedure for the progressive failure phenomenon can give the reliable result of the simulation for the heterogeneous soil slope thus enabling the more accurate way for hazard preparedness.     

多点地质统计学在河流相储层建模中的应用

多点地质统计学综合了基于象元方法以及基于目标方法两者的优点,对于河流相等具有复杂地质形态的储层精确建模具有较强的优势.在对传统建模方法综合分析的基础上,介绍了多点地质统计学的基本理论及SNESIM算法,并应用该技术对大牛地气田某开发井区的辫状分流河道相进行了实际建模.研究结果表明,在河流相储层建模中,该方法比传统的建模方法更具优越性.最后,进一步综合讨论了多点地质统计学目前面临的主要问题(包括训练图像、目标体连续性、数据样板选择、综合地震信息等方面)的改进方法.