共查询到20条相似文献，搜索用时 31 毫秒
1.
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.
Antoine Bertoncello Tao Sun Hongmei Li Gregoire Mariethoz Jef Caers 《Mathematical Geosciences》2013,45(7):873-893
Geostatistical simulation methods aim to represent spatial uncertainty through realizations that reflect a certain geological concept by means of a spatial continuity model. Most common spatial continuity models are either variogram, training image, or Boolean based. In this paper, a more recent spatial model of geological continuity is developed, termed the event, or surface-based model, which is specifically applicable to modeling cases with complex stratigraphy, such as in sedimentary systems. These methods rely on a rule-based stacking of events, which are mathematically represented by two-dimensional thickness variations over the domain, where positive thickness is associated with deposition and negative thickness with erosion. Although it has been demonstrated that the surface-based models accurately represent the geological variation present in complex layered systems, they are more difficult to constrain to hard and soft data as is typically required of practical geostatistical techniques. In this paper, we develop a practical methodology for constraining such models to hard data from wells and thickness data interpreted from geophysics, such as seismic data. Our iterative methodology relies on a decomposition of the parameter optimization problem into smaller, manageable problems that are solved sequentially. We demonstrate this method on a real case study of a turbidite sedimentary basin. 相似文献
3.
Geostatistical Seismic Inversion with Direct Sequential Simulation and Co-simulation with Multi-local Distribution Functions 总被引：1，自引：0，他引：1
Rúben Nunes Amílcar Soares Leonardo Azevedo Pedro Pereira 《Mathematical Geosciences》2017,49(5):583-601
Stochastic sequential simulation is a common modelling technique used in Earth sciences and an integral part of iterative geostatistical seismic inversion methodologies. Traditional stochastic sequential simulation techniques based on bi-point statistics assume, for the entire study area, stationarity of the spatial continuity pattern and a single probability distribution function, as revealed by a single variogram model and inferred from the available experimental data, respectively. In this paper, the traditional direct sequential simulation algorithm is extended to handle non-stationary natural phenomena. The proposed stochastic sequential simulation algorithm can take into consideration multiple regionalized spatial continuity patterns and probability distribution functions, depending on the spatial location of the grid node to be simulated. This work shows the application and discusses the benefits of the proposed stochastic sequential simulation as part of an iterative geostatistical seismic inversion methodology in two distinct geological environments in which non-stationarity behaviour can be assessed by the simultaneous interpretation of the available well-log and seismic reflection data. The results show that the elastic models generated by the proposed stochastic sequential simulation are able to reproduce simultaneously the regional and global variogram models and target distribution functions relative to the average volume of each sub-region. When used as part of a geostatistical seismic inversion procedure, the retrieved inverse models are more geologically realistic, since they incorporate the knowledge of the subsurface geology as provided, for example, by seismic and well-log data interpretation. 相似文献
4.
As additional multiple-point statistical (MPS) algorithms are developed, there is an increased need for scientific ways for comparison beyond the usual visual comparison or simple metrics, such as connectivity measures. In this paper, we start from the general observation that any (not just MPS) geostatistical simulation algorithm represents two types of variability: (1) the within-realization variability, namely, that realizations reproduce a spatial continuity model (variogram, Boolean, or training-image based), (2) the between-realization variability representing a model of spatial uncertainty. In this paper, it is argued that any comparison of algorithms needs, at a minimum, to be based on these two randomizations. In fact, for certain MPS algorithms, it is illustrated with different examples that there is often a trade-off: Increased pattern reproduction entails reduced spatial uncertainty. In this paper, the subjective choice that the best algorithm maximizes pattern reproduction is made while at the same time maximizes spatial uncertainty. The discussion is also limited to fairly standard multiple-point algorithms and that our method does not necessarily apply to more recent or possibly future developments. In order to render these fundamental principles quantitative, this paper relies on a distance-based measure for both within-realization variability (pattern reproduction) and between-realization variability (spatial uncertainty). It is illustrated in this paper that this method is efficient and effective for two-dimensional, three-dimensional, continuous, and discrete training images. 相似文献
5.
In previous studies, the groundwater flow models formulated for the Hat Yai Basin were conventional and deterministic because the geologic heterogeneity of the alluvial aquifer system in the basin had not yet been assessed. This paper describes an effort to develop hydrofacies models, such that the spatial variability of the aquifer system can be represented in a systematic way. Variogram parameters that characterize the alluvial aquifer heterogeneity were determined. Based on these variogram parameters, an indicator-based geostatistical approach was used to develop hydrofacies models using sequential indicator simulation. The hydrofacies models indicate three distinct aquifer units, namely Hat Yai, Khu Tao, and Kho Hong aquifers, which is in good agreement with a conceptual model, and incorporates spatial variability as observed in field data from borehole logs. The hydrofacies models can be used in groundwater modeling and simulations. 相似文献
6.
Wenlong Xu 《Mathematical Geology》1996,28(7):937-949
In geology, structures displaying differing local directions of continuity are widespread, a typical example being a flusial depositional system. Conventional pixel-based geostatistical algorithms, may fail to reproduce such curvilinear structures. Conversely, object-based algorithms can reproduce curvilinear shapes but are difficult to condition to dense local data. Local depositional directions as obtained from dipmeter data. 3D seismic data, and geological interpretation represent critical information. An improved pixel-based geostatistical algorithm is proposed to account for such directional information. Case studies demonstrate the potential and limitations of the algorithm. 相似文献
7.
Robust Resampling Confidence Intervals for Empirical Variograms 总被引：1，自引：0，他引：1
The variogram function is an important measure of the spatial dependencies of a geostatistical or other spatial dataset. It
plays a central role in kriging, designing spatial studies, and in understanding the spatial properties of geological and
environmental phenomena. It is therefore important to understand the variability attached to estimates of the variogram. Existing
methods for constructing confidence intervals around the empirical variogram either rely on strong assumptions, such as normality
or known variogram function, or are based on resampling blocks and subject to edge effect biases. This paper proposes two
new procedures for addressing these concerns: a quasi-block-bootstrap and a quasi-block-jackknife. The new methods are based
on transforming the data to decorrelate it based on a fitted variogram model, resampling blocks from the decorrelated data,
and then recorrelating. The coverage properties of the new confidence intervals are compared by simulation to a number of
existing resampling-based intervals. The proposed quasi-block-jackknife confidence interval is found to have the best properties
of all of the methods considered across a range of scenarios, including normally and lognormally distributed data and misspecification
of the variogram function used to decorrelate the data. 相似文献
8.
Using Sequential Indicator Simulation as a Tool in Reservoir Description: Issues and Uncertainties 总被引：3，自引：0，他引：3
Flow simulation studies require an accurate model of the reservoir in terms of its sedimentological architecture. Pixel-based reservoir modeling techniques are often used to model this architecture. There are, however, two problem areas with such techniques. First, several statistical parameters have to be provided whose influence on the resulting model is not readily inferable. Second, conditioning the models to relevant geological data that carry great uncertainty on their own adds to the difficulty of obtaining reliable models and assessing model reliability. The Sequential Indicator Simulation (SIS) method has been used to examine the impact of such uncertainties on the final reservoir model. The effects of varying variogram types, frequencies of lithology occurrence, and the gridblock model orientation with respect to the sedimentological trends are illustrated using different reservoir modeling studies. Results indicate, for example, that the choice of variogram type can have a significant impact on the facies model. Also, reproduction of sedimentological trends and large geometries requires careful parameter selection. By choosing the appropriate modeling strategy, sedimentological principles can be translated into the numerical model. Solutions for dealing with such issues and the geological uncertainties are presented. In conclusion, each reservoir modeling study should begin by developing a thorough quantitative sedimentological understanding of the reservoir under study, followed by detailed sensitivity analyses of relevant statistical and geological parameters. 相似文献
9.
Is the ocean floor a fractal? 总被引：1，自引：0，他引：1
The topographic structure of the ocean bottom is investigated at different scales of resolution to answer the question: Can the seafloor be described as a fractal process? Methods from geostatistics, the theory of regionalized variables, are used to analyze the spatial structure of the ocean floor at different scales of resolution. The key to the analysis is the variogram criterion: Self-similarity of a stochastic process implies self-similarity of its variogram. The criterion is derived and proved here: it also is valid for special cases of self-affinity (in a sense adequate for topography). It has been proposed that seafloor topography can be simulated as a fractal (an object of Hausdorff dimension strictly larger than its topological dimension), having scaling properties (self-similarity or self-affinity). The objective of this study is to compare the implications of these concepts with observations of the seafloor. The analyses are based on SEABEAM bathymetric data from the East Pacific Rise at 13°N/104°W and at 9°N/104°W and use tracks that run both across the ridge crest and along the ridge flank. In the geostatistical evaluation, the data are considered as a stochastic process. The spatial continuity of this process is described by variograms that are calculated for different scales and directions. Applications of the variogram criterion to scale-dependent variogram models yields the following results: Although the seafloor may be a fractal in the sense of the definition involving the Hausdorff dimension, it is not self-similar, nor self-affine (in the given sense). Mathematical models of scale-dependent spatial structures are presented, and their relationship to geologic processes such as ridge evolution, crust formation, and sedimentation is discussed. 相似文献
10.
Statistical characterization and stochastic modeling of pore networks in relation to fluid flow 总被引：18，自引：0，他引：18
R. D. Hazlett 《Mathematical Geology》1997,29(6):801-822
11.
The topographic structure of the ocean bottom is investigated at different scales of resolution to answer the question: Can the seafloor be described as a fractal process? Methods from geostatistics, the theory of regionalized variables, are used to analyze the spatial structure of the ocean floor at different scales of resolution. The key to the analysis is the variogram criterion: Self-similarity of a stochastic process implies self-similarity of its variogram. The criterion is derived and proved here: it also is valid for special cases of self-affinity (in a sense adequate for topography). It has been proposed that seafloor topography can be simulated as a fractal (an object of Hausdorff dimension strictly larger than its topological dimension), having scaling properties (self-similarity or self-affinity). The objective of this study is to compare the implications of these concepts with observations of the seafloor. The analyses are based on SEABEAM bathymetric data from the East Pacific Rise at 13°N/104°W and at 9°N/104°W and use tracks that run both across the ridge crest and along the ridge flank. In the geostatistical evaluation, the data are considered as a stochastic process. The spatial continuity of this process is described by variograms that are calculated for different scales and directions. Applications of the variogram criterion to scale-dependent variogram models yields the following results: Although the seafloor may be a fractal in the sense of the definition involving the Hausdorff dimension, it is not self-similar, nor self-affine (in the given sense). Mathematical models of scale-dependent spatial structures are presented, and their relationship to geologic processes such as ridge evolution, crust formation, and sedimentation is discussed. 相似文献
12.
为解决以往油藏描述工作单孔资料在深度和广度上存在明显局限性的问题,提出了一套井震联合地质建模的方法：首先以地震精细解释的层位数据和断层数据为约束,结合测井资料建立构造模型;然后采用地震反演体的概率面在平面上作为约束,纵向上在测井曲线上分类统计各微相的概率分布曲线,在变差函数分析的基础上,采用序贯指示模拟算法模拟出储层微相的空间展布,建立沉积相模型;最后采用相控技术和地质统计学理论,分析参数区域化变量的特征,建立储层属性模型。在大庆萨尔图油田的B1DD区块进行验证的结果表明,井震资料的联合应用揭示了大量的小断层,使断点组合率提高了15%,大于2 m的砂岩预测符合率达到90%以上,提高了井间断层和砂岩的预测精度,进而提高了地质模型、尤其是井间模型的精度。 相似文献
13.
Conditional Simulation of Complex Geological Structures Using Multiple-Point Statistics 总被引：68，自引：0，他引：68
Sebastien Strebelle 《Mathematical Geology》2002,34(1):1-21
In many earth sciences applications, the geological objects or structures to be reproduced are curvilinear, e.g., sand channels in a clastic reservoir. Their modeling requires multiple-point statistics involving jointly three or more points at a time, much beyond the traditional two-point variogram statistics. Actual data from the field being modeled, particularly if it is subsurface, are rarely enough to allow inference of such multiple-point statistics. The approach proposed in this paper consists of borrowing the required multiple-point statistics from training images depicting the expected patterns of geological heterogeneities. Several training images can be used, reflecting different scales of variability and styles of heterogeneities. The multiple-point statistics inferred from these training image(s) are exported to the geostatistical numerical model where they are anchored to the actual data, both hard and soft, in a sequential simulation mode. The algorithm and code developed are tested for the simulation of a fluvial hydrocarbon reservoir with meandering channels. The methodology proposed appears to be simple (multiple-point statistics are scanned directly from training images), general (any type of random geometry can be considered), and fast enough to handle large 3D simulation grids. 相似文献
14.
E. S. Zakirov I. M. Indrupskiy O. V. Liubimova I. M. Shiriaev D. P. Anikeev 《Doklady Earth Sciences》2017,476(2):1120-1124
An approach for geostatistically consistent matching of 3D flow simulation models and 3D geological models is proposed. This approach uses an optimization algorithm based on identification of the parameters of the geostatistical model (for example, the variogram parameters, such as range, sill, and nugget effect). Here, the inverse problem is considered in the greatest generality taking into account facies heterogeneity and the variogram anisotropy. The correlation dependence parameters (porosity-to-log permeability) are clarified for each single facies. 相似文献
15.
G. Bourgault 《Mathematical Geology》1997,29(3):315-334
Parametric geostatistical simulations such as LU decomposition and sequential algorithms do not need Gaussian distributions.
It is shown that variogram model reproduction is obtained when Uniform or Dipole distributions are used instead of Gaussian
distributions for drawing i. i.d. random values in LU simulation, or for modeling the local conditional probability distributions
in sequential simulation. Both algorithms yield simulated values with a marginal normal distribution no matter if Gaussian,
Uniform, or Dipole distributions are used. The range of simulated values decreases as the entropy of the probability distribution
decreases. Using Gaussian distributions provides a larger range of simulated normal score values than using Uniform or Dipole
distributions. This feature has a negligible effect for reproduction of the normal scores variogram model but have a larger
impact on the reproduction of the original values variogram. The Uniform or Dipole distributions also produce lesser fluctuations
among the variograms of the simulated realizations. 相似文献
16.
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. 相似文献
17.
《International Geology Review》2012,54(8):676-682
Knowledge of the spatial variability of grade is important for commercial evaluation of mineral deposits. The low-grade, Kizilyuksek-Yataardic chromite deposit in Karsanti, Adana is one of the greatest reserves in western Anatolia, Turkey and offers great potential. However this deposit lacks geologic information needed to provide a better understanding of the spatial variability of Cr2O3%. In this paper, geostatistical tools are used to improve the understanding of the vertical variability of Cr2O3%. The vertical variability is characterized by an average down-hole variogram from core measurements. The variogram exhibits nested structures at 10, 50, and 100 m. Because geologic information is not adequate to explain the source of these features, geostatistical simulation of Cr2O3% with depth is used. The results indicate that the chromite mineralization is characterized by zonal structure, with average size of the zones being 50 m, and average distance between them being 100 m. A study of the most variable holes also shows that these holes have the greatest effect on the average down-hole variogram. 相似文献
18.
地质统计学从解决固体矿产资源评价起源，应用领域迅速扩展。文章简要介绍了地质统计学的理论研究现状和应用现状。为了减轻人为因素在创建模型中的作用，非参数的地质统计学得到发展，文章重点介绍了非参数地质统计学的4种变异函数模型。此外，文章阐述了多点地质统计学的概念，这是一种用训练映像代替变异函数的地质统计学方法，使得地质现象的解释更直观。 相似文献
19.
Challenges in reservoir forecasting 总被引：3，自引：0，他引：3
The combination of geostatistics-based numerical geological models and finite difference flow simulation has improved our ability to predict reservoir performance. The main contribution of geostatistical modeling has been more realistic representations of reservoir heterogeneity. Our understanding of the physics of fluid flow in porous media is reasonably captured by flow simulators in common usage. Notwithstanding the increasing application and success of geostatistics and flow simulation there remain many important challenges in reservoir forecasting. This application has alerted geoscientists and physicists that geostatistical/flow models in many respects, are, engineering approximations to thereal spatial distribution andreal flow processes. This paper reviews current research directions and presents some new ideas of where reserach could be focused to improve our ability to model geological features, model flow processes, and, ultimately, improve reservoir performance predictions. 相似文献
20.
Non-stationary Geostatistical Modeling Based on Distance Weighted Statistics and Distributions 总被引：3，自引：0，他引：3
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. 相似文献