共查询到20条相似文献，搜索用时 31 毫秒
1.
Teacher''s Aide Variogram Interpretation and Modeling 总被引：13，自引：0，他引：13
The variogram is a critical input to geostatistical studies: (1) it is a tool to investigate and quantify the spatial variability of the phenomenon under study, and (2) most geostatistical estimation or simulation algorithms require an analytical variogram model, which they will reproduce with statistical fluctuations. In the construction of numerical models, the variogram reflects some of our understanding of the geometry and continuity of the variable, and can have a very important impact on predictions from such numerical models. The principles of variogram modeling are developed and illustrated with a number of practical examples. A three-dimensional interpretation of the variogram is necessary to fully describe geologic continuity. Directional continuity must be described simultaneously to be consistent with principles of geological deposition and for a legitimate measure of spatial variability for geostatistical modeling algorithms. Interpretation principles are discussed in detail. Variograms are modeled with particular functions for reasons of mathematical consistency. Used correctly, such variogram models account for the experimental data, geological interpretation, and analogue information. The steps in this essential data integration exercise are described in detail through the introduction of a rigorous methodology. 相似文献
2.
The variogram is a critical input to geostatistical studies: (1) it is a tool to investigate and quantify the spatial variability of the phenomenon under study, and (2) most geostatistical estimation or simulation algorithms require an analytical variogram model, which they will reproduce with statistical fluctuations. In the construction of numerical models, the variogram reflects some of our understanding of the geometry and continuity of the variable, and can have a very important impact on predictions from such numerical models. The principles of variogram modeling are developed and illustrated with a number of practical examples. A three-dimensional interpretation of the variogram is necessary to fully describe geologic continuity. Directional continuity must be described simultaneously to be consistent with principles of geological deposition and for a legitimate measure of spatial variability for geostatistical modeling algorithms. Interpretation principles are discussed in detail. Variograms are modeled with particular functions for reasons of mathematical consistency. Used correctly, such variogram models account for the experimental data, geological interpretation, and analogue information. The steps in this essential data integration exercise are described in detail through the introduction of a rigorous methodology. 相似文献
3.
湛江组灰色黏土空间展布规律研究 总被引：1，自引：0，他引：1
分布于湛江地区的湛江组灰色黏土因其特有的强结构性在实际工程中产生一系列岩土工程问题,20世纪80年代已为国内众学者所关注。由于湛江组地层在雷州半岛区域内,具有不同沉积环境而表现出不同空间展布规律,实际工程中难以掌握其分布特性。根据现场地质调查资料,结合雷州半岛的地形地貌、地质构造以及湛江组地层的沉积环境,在雷州半岛布置27个钻孔进行地质钻探取样。不同区域典型地层剖面图说明,灰色黏土在不同的沉积环境有着不同沉积韵律、埋藏深度、地层厚度等,分布具有区域性,根据层顶埋深及地层厚度分析灰色黏土的结构性对实际工程的影响。同时,根据钻孔地层数据,利用基础 GIS软件平台MAPGIS构建湛江组结构性黏土空间展布规律三维地质模型,全面展示出灰色黏土在雷州半岛的整个展布规律,实现灰色黏土空间展布的可视化,最后简单地介绍该模型的基本应用。 相似文献
4.
Susan J. Webb Lewis D. Ashwal R. Grant Cawthorn 《Contributions to Mineralogy and Petrology》2011,162(1):101-107
The three-dimensional shapes of mafic layered intrusions have to be inferred from surface outcrops, in some cases aided by drilling and/or geophysical data. However, geophysical models are often equivocal. For the 2.06?Ga Bushveld Complex of South Africa, early geological models proposed a shape of a single, gently inward dipping lopolith. Subsequent resistivity and gravity data were interpreted to suggest that the eastern and western limbs were discrete, dipping wedge-shaped intrusions separated by ~150?km. A more recent gravity model that takes crustal flexure into account allows continuity and the reversal to the original model. Distinguishing between these possibilities is difficult from surface-based studies because the central regions of the Complex are obscured by large volumes of younger granites and sedimentary/volcanic cover rocks. Here, we describe xenoliths from the Cretaceous Palmietgat kimberlite pipe, located mid-way between the exposed western and eastern lobes of the Complex. They are chromite-bearing feldspathic pyroxenites considered equivalent to those of the typical outcropping Critical Zone of the Bushveld Complex. This result provides strong support for a regionally interconnected Bushveld Complex, implying its emplacement as a single sill-like body. Confirming the continuity of the Bushveld Complex greatly expands exploration opportunities and implies that other layered mafic intrusions could have similar geometry. 相似文献
5.
《地学前缘(英文版)》2022,13(6):101435
The past two decades have seen a rapid adoption of artificial intelligence methods applied to mineral exploration. More recently, the easier acquisition of some types of data has inspired a broad literature that has examined many machine learning and modelling techniques that combine exploration criteria, or ‘features’, to generate predictions for mineral prospectivity. Central to the design of prospectivity models is a ‘mineral system’, a conceptual model describing the key geological elements that control the timing and location of economic mineralisation. The mineral systems model defines what constitutes a training set, which features represent geological evidence of mineralisation, how features are engineered and what modelling methods are used. Mineral systems are knowledge-driven conceptual models, thus all parameter choices are subject to human biases and opinion so alternative models are possible. However, the effect of alternative mineral systems models on prospectivity is rarely compared despite the potential to heavily influence final predictions. In this study, we focus on the effect of conceptual uncertainty on Fe ore prospectivity models in the Hamersley region, Western Australia. Four important considerations are tested. (1) Five different supergene and hypogene conceptual mineral systems models guide the inputs for five forest-based classification prospectivity models model. (2) To represent conceptual uncertainty, the predictions are then combined for prospectivity model comparison. (3) Representation of three-dimensional objects as two-dimensional features are tested to address commonly ignored thickness of geological units. (4) The training dataset is composed of known economic mineralisation sites (deposits) as ‘positive’ examples, and exploration drilling data providing ‘negative’ sampling locations. Each of the spatial predictions are assessed using independent performance metrics common to AI-based classification methods and subjected to geological plausibility testing. We find that different conceptual mineral systems produce significantly different spatial predictions, thus conceptual uncertainty must be recognised. A benefit to recognising and modelling different conceptual models is that robust and geologically plausible predictions can be made that may guide mineral discovery. 相似文献
6.
层状地质体的三维模拟与可视化 总被引：51，自引：1，他引：50
文中探讨了层状地质体的三维模拟与可视化方法和技术。地质数据通常是有限的、空间上分布不均匀的点数据 ,需要大量使用插值技术来形成合理的模型。简单层状实体可在二维插值的基础上 ,按照这些面的空间位置叠加 ,形成三维体模型。复杂的层状体包括沉积相以及一些物性参数的空间变化 ,断层等形成的不连续以及复杂的褶皱等。这些实体的模拟与可视化 ,需要进行三维插值 ,并根据不同现象的特点选择合适的方法。合理的专家知识或地质解释的交互加入 ,是形成合理的三维模型的关键。 相似文献
7.
利用相控模型进行井间参数预测 总被引：18，自引：0，他引：18
储层物性参数在空间的展布，一直是储层非均质性研究的难点，利用分形和蒙特卡洛方法，在精细的沉积微相及小层对比研究基础之上，预测了井间的岩石相分布，并根据取心井资料，建立了各种岩石相的渗透率相控模型，从而预测出井间渗透率的分布。该方法以地质特征为基础，将地质资料与数学方法相结合，充分体现了渗透率在井间既具有连续性，同时变化又非常剧烈的特点，为储层物性空间预测提供了新的思路。 相似文献
8.
Application of multiple-point geostatistics on modelling groundwater flow and transport in a cross-bedded aquifer (Belgium) 总被引：3，自引：1，他引：2
Sedimentological processes often result in complex three-dimensional subsurface heterogeneity of hydrogeological parameter values. Variogram-based stochastic approaches are often not able to describe heterogeneity in such complex geological environments. This work shows how multiple-point geostatistics can be applied in a realistic hydrogeological application to determine the impact of complex geological heterogeneity on groundwater flow and transport. The approach is applied to a real aquifer in Belgium that exhibits a complex sedimentary heterogeneity and anisotropy. A training image is constructed based on geological and hydrogeological field data. Multiple-point statistics are borrowed from this training image to simulate hydrofacies occurrence, while intrafacies permeability variability is simulated using conventional variogram-based geostatistical methods. The simulated hydraulic conductivity realizations are used as input to a groundwater flow and transport model to investigate the effect of small-scale sedimentary heterogeneity on contaminant plume migration. Results show that small-scale sedimentary heterogeneity has a significant effect on contaminant transport in the studied aquifer. The uncertainty on the spatial facies distribution and intrafacies hydraulic conductivity distribution results in a significant uncertainty on the calculated concentration distribution. Comparison with standard variogram-based techniques shows that multiple-point geostatistics allow better reproduction of irregularly shaped low-permeability clay drapes that influence solute transport. 相似文献
9.
J A Vargas-Guzmán 《中国地质大学学报(英文版)》2008,19(1):85-92
An analysis of statistical expected values for transformations is performed in this study to quantify the effect of heterogeneity on spatial geological modeling and evaluations. Algebraic transformations are frequently applied to data from logging to allow for the modeling of geological properties. Transformations may be powers, products, and exponential operations which are commonly used in well-known relations (e.g., porosity-permeability transforms). The results of this study show that correct computations must account for residual transformation terms which arise due to lack of independence among heterogeneous geological properties. In the case of an exponential porosity-permeability transform, the values may be positive. This proves that a simple exponential model back-transformed from linear regression underestimates permeability. In the case of transformations involving two or more properties, residual terms may represent the contribution of heterogeneous components which occur when properties vary together, regardless of a pair-wise linear independence. A consequence of power- and product-transform models is that regression equationswithin those transformations need corrections via residual cumulants. A generalization of this result isthat transformations of multivariate spatial attributes require multiple-point random variable relations. This analysis provides practical solutions leading to a methodology for nonlinear modeling using correct back transformations in geology. 相似文献
10.
11.
Characterization of complex geological features and patterns remains one of the most challenging tasks in geostatistics. Multiple point statistics (MPS) simulation offers an alternative to accomplish this aim by going beyond classical two-point statistics. Reproduction of features in the final realizations is achieved by borrowing high-order spatial statistics from a training image. Most MPS algorithms use one training image at a time chosen by the geomodeler. This paper proposes the use of multiple training images simultaneously for spatial modeling through a scheme of data integration for conditional probabilities known as a linear opinion pool. The training images (TIs) are based on the available information and not on conceptual geological models; one image comes from modeling the categories by a deterministic approach and another comes from the application of conventional sequential indicator simulation. The first is too continuous and the second too random. The mixing of TIs requires weights for each of them. A methodology for calibrating the weights based on the available drillholes is proposed. A measure of multipoint entropy along the drillholes is matched by the combination of the two TIs. The proposed methodology reproduces geologic features from both TIs with the correct amount of continuity and variability. There is no need for a conceptual training image from another modeling technique; the data-driven TIs permit a robust inference of spatial structure from reasonably spaced drillhole data. 相似文献
12.
准噶尔盆地车排子凸起新近系沙湾组为湖盆缓坡多物源、复杂沉积体系沉积,其沉积体系认识争议很大,时空演化难以刻画.准确刻画其沉积体系类型、演化规律和砂体分布特征,是精细刻画单砂体及内部构型等储层的基石,对后期油气开发具有重要指导意义.在岩心、测井、测试分析及地震资料,充分利用相位转换、地层切片等地震沉积学技术方法,对其沉积体系特征进行深入研究.建立研究区精细层序地层格架,将沙湾组划分为2个三级层序、4个四级层序;识别出冲积扇、辫状河、曲流河、曲流河三角洲、扇三角洲及滨浅湖6种沉积体系类型,并确定沉积体系边界及演化规律;建立了多物源冲积扇-辫状河、扇三角洲-曲流河三角洲沉积模式及砂体分布模式.辫状河道、心滩、点坝复合体及冲积扇砂体均钻遇较好油气显示,为重要储集体.此次研究丰富了地震沉积学在湖盆缓坡多物源、多沉积体系类型中的应用. 相似文献
13.
全面系统地研究了沉积地层系统三维实体构模的原理、方法及实施过程。针对沉积地层系统的几何结构特征,提出了一种基于钻孔数据的三维实体模型生成方法“钻孔-层面-实体法”。该方法首先使用二分拓扑结构将钻孔数据离散化为一系列散点,基于这些散点插值拟合出各个地层顶、底界面的初始高程,自动推算出缺失地层的地质成因;然后按照不同的地质成因,自动对缺失地层及其控制界面进行交切处理、高程调整与一致性处理;最后生成用三维块体或三棱柱体元充填的三维实体模型。通过上海城建工程中的一个应用实例验证了“钻孔-层面-实体法”的有效性。 相似文献
14.
厚结壳的形成条件及控制因素分析 总被引：2，自引：0，他引：2
当前对于大洋厚结壳的研究较少,而厚结壳是未来开发利用的重点,同时它记录的地质历史长,控制和影响结壳生长发育的各种因素必然能在厚结壳上充分体现出来。因此,开展厚结壳的形成条件及控制因素的研究,具有重要的理论与现实意义。文章利用GIS空间分析技术与地质统计方法,对厚结壳的定义、结壳厚度与水深、结壳厚度与基岩、结壳生长与构造活动、结壳厚度与沉积速率、结壳厚度与经纬度、结壳厚度与主成矿元素之间的关系进行了分析,结果表明,结壳的厚度变化受控于水深、基岩、构造活动、沉积速率、经纬度等因素,同时结壳厚度与主成矿元素之间存在明显的相关性。 相似文献
15.
FLAC3D平台下前处理阶段复杂地质体三维数值计算模型的建立一直是岩土工程数值模拟工作者面临的难题,而三维激光扫描技术在获取复杂物体的精细空间数据的方面具有独特的优势。因此,本文提出以逆向工程为过渡平台实现复杂三维地质体FLAC3D数值模型的构建。利用三维激光扫描技术精细几何数据的获取功能,对复杂三维地质体空间形状几何信息进行提取;利用Geomagic Studio良好的点云数据处理及曲面建模功能,准确地重构了复杂地质体的CAD曲面模型;利用Hypermesh强大的几何处理能力及网格划分能力,将复杂地质体的CAD曲面模型进行实体化和网格化;最后由FLAC3D软件内嵌的Fish语言编制的文件转换为FLAC3D能够识别的文件格式,从而实现了该软件前处理过程中复杂三维地质模型的高效、精准的建立。结合工程实例,对该方法的可行性与仿真效果进行了检验。结果表明该方法具有方便、快捷、合理且实用性强的优点。 相似文献
16.
Application of EM algorithms for seismic facices classification 总被引：1，自引：0，他引：1
Identification of the geological facies and their distribution from seismic and other available geological information is
important during the early stage of reservoir development (e.g. decision on initial well locations). Traditionally, this is
done by manually inspecting the signatures of the seismic attribute maps, which is very time-consuming. This paper proposes
an application of the Expectation-Maximization (EM) algorithm to automatically identify geological facies from seismic data.
While the properties within a certain geological facies are relatively homogeneous, the properties between geological facies
can be rather different. Assuming that noisy seismic data of a geological facies, which reflect rock properties, can be approximated
with a Gaussian distribution, the seismic data of a reservoir composed of several geological facies are samples from a Gaussian
mixture model. The mean of each Gaussian model represents the average value of the seismic data within each facies while the
variance gives the variation of the seismic data within a facies. The proportions in the Gaussian mixture model represent
the relative volumes of different facies in the reservoir. In this setting, the facies classification problem becomes a problem
of estimating the parameters defining the Gaussian mixture model. The EM algorithm has long been used to estimate Gaussian
mixture model parameters. As the standard EM algorithm does not consider spatial relationship among data, it can generate
spatially scattered seismic facies which is physically unrealistic. We improve the standard EM algorithm by adding a spatial
constraint to enhance spatial continuity of the estimated geological facies. By applying the EM algorithms to acoustic impedance
and Poisson’s ratio data for two synthetic examples, we are able to identify the facies distribution. 相似文献
17.
The reproduction of the non-stationary distribution and detailed characteristics of geological bodies is the main difficulty of reservoir modeling. Recently developed multiple-point geostatistics can represent a stationary geological body more effectively than traditional methods. When restricted to a stationary hypothesis, multiple-point geostatistical methods cannot simulate a non-stationary geological body effectively, especially when using non-stationary training images (TIs). According to geologic principles, the non-stationary distribution of geological bodies is controlled by a sedimentary model. Therefore, in this paper, we propose auxiliary variables based on the sedimentary model, namely geological vector information (GVI). GVI can characterize the non-stationary distribution of TIs and simulation domains before sequential simulation, and the precision of data event statistics will be enhanced by the sequential simulation’s data event search area limitations under the guidance of GVI. Consequently, the reproduction of non-stationary geological bodies will be improved. The key features of this method are as follows: (1) obtain TIs and geological vector information for simulated areas restricted by sedimentary models; (2) truncate TIs into a number of sub-TIs using a set of cut-off values such that each sub-TI is stationary and the adjacent sub-TIs have a certain similarity; (3) truncate the simulation domain into a number of sub-regions with the same cut-off values used in TI truncation, so that each sub-region corresponds to a number of sub-TIs; (4) use an improved method to scan the TI or TIs and construct a single search tree to restore replicates of data events located in different sub-TIs; and (5) use an improved conditional probability distribution function to perform sequential simulation. A FORTRAN program is implemented based on the SNESIM. 相似文献
18.
A Distance-based Prior Model Parameterization for Constraining Solutions of Spatial Inverse Problems 总被引：4，自引：3，他引：1
Spatial inverse problems in the Earth Sciences are often ill-posed, requiring the specification of a prior model to constrain
the nature of the inverse solutions. Otherwise, inverted model realizations lack geological realism. In spatial modeling,
such prior model determines the spatial variability of the inverse solution, for example as constrained by a variogram, a
Boolean model, or a training image-based model. In many cases, particularly in subsurface modeling, one lacks the amount of
data to fully determine the nature of the spatial variability. For example, many different training images could be proposed
for a given study area. Such alternative training images or scenarios relate to the different possible geological concepts
each exhibiting a distinctive geological architecture. Many inverse methods rely on priors that represent a single subjectively
chosen geological concept (a single variogram within a multi-Gaussian model or a single training image). This paper proposes
a novel and practical parameterization of the prior model allowing several discrete choices of geological architectures within
the prior. This method does not attempt to parameterize the possibly complex architectures by a set of model parameters. Instead,
a large set of prior model realizations is provided in advance, by means of Monte Carlo simulation, where the training image
is randomized. The parameterization is achieved by defining a metric space which accommodates this large set of model realizations.
This metric space is equipped with a “similarity distance” function or a distance function that measures the similarity of
geometry between any two model realizations relevant to the problem at hand. Through examples, inverse solutions can be efficiently
found in this metric space using a simple stochastic search method. 相似文献
19.
Simultaneous Estimation of Geologic and Reservoir State Variables Within an Ensemble-Based Multiple-Point Statistic Framework 总被引：1，自引：1，他引：0
Liangping Li Sanjay Srinivasan Haiyan Zhou J. Jaime Gómez-Hernández 《Mathematical Geosciences》2014,46(5):597-623
Assessment of uncertainty due to inadequate data and imperfect geological knowledge is an essential aspect of the subsurface model building process. In this work, a novel methodology for characterizing complex geological structures is presented that integrates dynamic data. The procedure results in the assessment of uncertainty associated with the predictions of flow and transport. The methodology is an extension of a previously developed pattern search-based inverse method that models the spatial variation in flow parameters by searching for patterns in an ensemble of reservoir models. More specifically, the pattern-searching algorithm is extended in two directions: (1) state values (such as piezometric head) and parameters (such as conductivities) are simultaneously and sequentially estimated, which implies that real-time assimilation of dynamic data is possible as in ensemble filtering approaches; and (2) both the estimated parameter and state variables are considered when pattern searching is implemented. The new scheme results in two main advantages—better characterization of parameters, especially for delineating small scale features, and an ensemble of head states that can be used to update the parameter field using the dynamic data at the next instant, without running expensive flow simulations. An efficient algorithm for pattern search is developed, which works with a flexible search radius and can be optimized for the estimation of either large- or small-scale structures. Synthetic examples are employed to demonstrate the effectiveness and robustness of the proposed approach. 相似文献
20.
High-order Statistics of Spatial Random Fields: Exploring Spatial Cumulants for Modeling Complex Non-Gaussian and Non-linear Phenomena 总被引：9，自引：7，他引：2
The spatial distributions of earth science and engineering phenomena under study are currently predicted from finite measurements
and second-order geostatistical models. The latter models can be limiting, as geological systems are highly complex, non-Gaussian,
and exhibit non-linear patterns of spatial connectivity. Non-linear and non-Gaussian high-order geostatistics based on spatial
connectivity measures, namely spatial cumulants, are proposed as a new alternative modeling framework for spatial data. This
framework has two parts. The first part is the definition, properties, and inference of spatial cumulants—including understanding
the interrelation of cumulant characteristics with the in-situ behavior of geological entities or processes, as examined in
this paper. The second part is the research on a random field model for simulation based on its high-order spatial cumulants.
Mathematical definitions of non-Gaussian spatial random functions and their high-order spatial statistics are presented herein,
stressing the notion of spatial cumulants. The calculation of spatial cumulants with spatial templates follows, including
anisotropic experimental cumulants. Several examples of two- and three-dimensional images, including a diamond bearing kimberlite
pipe from the Ekati Mine in Canada, are analyzed to assess the relations between cumulants and the spatial behavior of geological
processes. Spatial cumulants of orders three to five are shown to capture directional multiple-point periodicity, connectivity
including connectivity of extreme values, and spatial architecture. In addition, they provide substantial information on geometric
characteristics and anisotropy of geological patterns. It is further shown that effects of complex spatial patterns are seen
even if only subsets of all cumulant templates are computed. Compared to second-order statistics, cumulant maps are found
to include a wealth of additional information from underlying geological patterns. Further work seeks to integrate this information
in the predictive capabilities of a random field model. 相似文献