共查询到20条相似文献，搜索用时 31 毫秒
1.
L. Y. Hu Y. Liu C. Scheepens A. W. Shultz R. D. Thompson 《Mathematical Geosciences》2014,46(2):227-240
The multiple-point simulation (MPS) method has been increasingly used to describe the complex geologic features of petroleum reservoirs. The MPS method is based on multiple-point statistics from training images that represent geologic patterns of the reservoir heterogeneity. The traditional MPS algorithm, however, requires the training images to be stationary in space, although the spatial distribution of geologic patterns/features is usually nonstationary. Building geologically realistic but statistically stationary training images is somehow contradictory for reservoir modelers. In recent research on MPS, the concept of a training image has been widely extended. The MPS approach is no longer restricted by the size or the stationarity of training images; a training image can be a small geometrical element or a full-field reservoir model. In this paper, the different types of training images and their corresponding MPS algorithms are first reviewed. Then focus is placed on a case where a reservoir model exists, but needs to be conditioned to well data. The existing model can be built by process-based, object-based, or any other type of reservoir modeling approach. In general, the geologic patterns in a reservoir model are constrained by depositional environment, seismic data, or other trend maps. Thus, they are nonstationary, in the sense that they are location dependent. A new MPS algorithm is proposed that can use any existing model as training image and condition it to well data. In particular, this algorithm is a practical solution for conditioning geologic-process-based reservoir models to well data. 相似文献
2.
利用三维地质模拟技术重构地质现象的三维空间分布,是实现自然资源管理和风险评估的重要基础和前提。多点统计学方法通过探寻多点间的空间结构关系,结合随机模拟方法生成具有差异性的模拟结果,较好地再现了复杂的地质现象。然而,如何构建合适、有效的训练图像一直是基于多点统计学三维地质模拟的核心问题。本文提出了一种改进的多点统计学算法。本方法结合了序贯模拟和迭代的方法,将二维剖面扩展为三维训练图像,再结合EM-Like算法,实现了三维地质结构的优化模拟。建模实例结果表明,本方法能确保训练图像对内部模拟网格的约束,准确模拟研究区的地层层序,并很好地再现二维地质剖面所反映的地层结构关系。 相似文献
3.
Tobias Lochbühler Guillaume Pirot Julien Straubhaar Niklas Linde 《Mathematical Geosciences》2014,46(5):625-645
Geophysical tomography captures the spatial distribution of the underlying geophysical property at a relatively high resolution, but the tomographic images tend to be blurred representations of reality and generally fail to reproduce sharp interfaces. Such models may cause significant bias when taken as a basis for predictive flow and transport modeling and are unsuitable for uncertainty assessment. We present a methodology in which tomograms are used to condition multiple-point statistics (MPS) simulations. A large set of geologically reasonable facies realizations and their corresponding synthetically calculated cross-hole radar tomograms are used as a training image. The training image is scanned with a direct sampling algorithm for patterns in the conditioning tomogram, while accounting for the spatially varying resolution of the tomograms. In a post-processing step, only those conditional simulations that predicted the radar traveltimes within the expected data error levels are accepted. The methodology is demonstrated on a two-facies example featuring channels and an aquifer analog of alluvial sedimentary structures with five facies. For both cases, MPS simulations exhibit the sharp interfaces and the geological patterns found in the training image. Compared to unconditioned MPS simulations, the uncertainty in transport predictions is markedly decreased for simulations conditioned to tomograms. As an improvement to other approaches relying on classical smoothness-constrained geophysical tomography, the proposed method allows for: (1) reproduction of sharp interfaces, (2) incorporation of realistic geological constraints and (3) generation of multiple realizations that enables uncertainty assessment. 相似文献
4.
多点地质统计学建模是近年来储层建模技术的研究热点,其实用性受到训练图像的限制。训练图像的质量决定了多点地质统计学建模的精度和可靠程度,是多点地质统计学建模成功的关键因素。文章阐述了训练图像的特征和意义,从方法的定义、使用情况、实例等方面系统介绍了训练图像建立的方法,包括手工绘制、基于目标模拟、三维地震信息提取或转化、基于原型模型、基于过程模拟和二维图像方法,综合对比了不同训练图像建立方法的数据来源、优势与不足,探讨了多点地质统计学建模依赖训练图像存在的问题。结合文献调研和多点地质统计学建模实践,指出了训练图像及其建立方法的发展方向,为多点地质统计学建模研究者和使用者提供借鉴,为完善多点地质统计学建模方法提供思考。 相似文献
5.
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. 相似文献
6.
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. 相似文献
7.
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. 相似文献
8.
Yu. B. Gladenkov 《Stratigraphy and Geological Correlation》2010,18(6):660-673
The problems of zonal stratigraphy of the early 21st century are discussed. The great advances achieved in recent years in
using zones in geological practice are noted. At the same time, attention is drawn to the controversies existing in the interpretation
of the concepts “biostratigraphic zone” and “chronozone,” in the methods of drawing the boundaries of such zones, and in the
assessment of the spatial scale of zones and in the understanding of them as stratigraphic units. 相似文献
9.
The resolution of measurement devices can be insufficient for certain purposes. We propose to stochastically simulate spatial
features at scales smaller than the measurement resolution. This is accomplished using multiple-point geostatistical simulation
(direct sampling in the present case) to interpolate values at the target scale. These structures are inferred using hypothesis
of scale invariance and stationarity on the spatial patterns found at the coarse scale. The proposed multiple-point super-resolution
mapping method is able to deal with “both continuous and categorical variables,” and can be extended to multivariate problems.
The advantages and limitations of the approach are illustrated with examples from satellite imaging. 相似文献
10.
Training Images from Process-Imitating Methods 总被引：2，自引：2，他引：0
Alessandro Comunian Sanjeev K. Jha Beatrice M. S. Giambastiani Gregoire Mariethoz Bryce F. J. Kelly 《Mathematical Geosciences》2014,46(2):241-260
The lack of a suitable training image is one of the main limitations of the application of multiple-point statistics (MPS) for the characterization of heterogeneity in real case studies. Process-imitating facies modeling techniques can potentially provide training images. However, the parameterization of these process-imitating techniques is not straightforward. Moreover, reproducing the resulting heterogeneous patterns with standard MPS can be challenging. Here the statistical properties of the paleoclimatic data set are used to select the best parameter sets for the process-imitating methods. The data set is composed of 278 lithological logs drilled in the lower Namoi catchment, New South Wales, Australia. A good understanding of the hydrogeological connectivity of this aquifer is needed to tackle groundwater management issues. The spatial variability of the facies within the lithological logs and calculated models is measured using fractal dimension, transition probability, and vertical facies proportion. To accommodate the vertical proportions trend of the data set, four different training images are simulated. The grain size is simulated alongside the lithological codes and used as an auxiliary variable in the direct sampling implementation of MPS. In this way, one can obtain conditional MPS simulations that preserve the quality and the realism of the training images simulated with the process-imitating method. The main outcome of this study is the possibility of obtaining MPS simulations that respect the statistical properties observed in the real data set and honor the observed conditioning data, while preserving the complex heterogeneity generated by the process-imitating method. In addition, it is demonstrated that an equilibrium of good fit among all the statistical properties of the data set should be considered when selecting a suitable set of parameters for the process-imitating simulations. 相似文献
11.
Discrete multifractals 总被引：12，自引：0，他引：12
Qiuming Cheng 《Mathematical Geology》1997,29(2):245-266
The concept of multifractal modeling has been used intensively in various fields of science for characterizing measures with
self- similarity. It has been shown that multifractal modeling provides powerful tools for characterizing patterns in the
spatial distribution of geological quantities and objects. Existing multifractal models were proposed for the purpose of handling
spatially intertwined fractals with continuous fractal spectrum f(α) (or continuous codimension function C(γ)). In this paper, these conventional multifractals are termed “continuous multifractals” whereas multifractals with discrete
fractal dimensions are termed “discrete multifractals.” The properties of discrete multifractals are investigated. It is shown
by various artificial examples and a case study of stratigraphy of Ocean Drilling Program (ODP) site 645 that spatially intertwined
fractals/multifractals indeed can have discrete fractal dimensions. Histogram-and moment-based techniques are proposed for
discrete multifractal modeling and illustrated using the artificial examples. The new concept of discrete multifractals and
associated multifractal modeling yields not only techniques for characterizing multifractals with discrete fractal dimensions
but it also provides insight into the relationships between fractals, bifractals, and multifractals. 相似文献
12.
Avalanche hazard mapping over large undocumented areas 总被引：3，自引：1，他引：2
An innovative methodology to perform avalanche hazard mapping over large undocumented areas is herewith presented and discussed.
The method combines GIS tools, computational routines, and statistical analysis in order to provide a “semi-automatic” definition
of areas potentially affected by avalanche release and motion. The method includes two main modules. The first module is used
to define zones of potential avalanche release, based on the consolidated relations on slope, morphology, and vegetation.
For each of the identified zones of potential release, a second module, named Avalanche Flow and Run-out Algorithm (AFRA),
provides an automatic definition of the areas potentially affected by avalanche motion and run-out. The definition is generated
by a specifically implemented “flow-routing algorithm” which allows for the determination of flow behaviour in the track and
in the run-out zone. In order to estimate the avalanche outline in the run-out zone, AFRA uses a “run-out cone”, which is
a 3D projection of the angle of reach α. The α-value is evaluated by statistical analysis of historical data regarding extreme
avalanches. Pre- and post-processing of the AFRA input/output data is done in an open source GIS environment (GRASS GIS).
The method requires only a digital terrain model and an indication of the areas covered by forest as input parameters. The
procedure, which allows rapid mapping of large areas, does not in principle require any site-specific historical information.
Furthermore, it has proven to be effective in all cases where a preliminary cost-efficient analysis of the territories potentially
affected by snow avalanche was needed. 相似文献
13.
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.” 相似文献
14.
To predict the macroscopic properties (e.g., transport, electromagnetic, and mechanical properties) of porous media, it is necessary to have a three‐dimensional (3D) representation of porous media. We reconstruct the geologically realistic 3D structure of Fontainebleau sandstone based on the two‐dimensional (2D) thin sections by using the multiple‐point statistics method. For this method, the size of template is an important parameter that reflects the perceived scale of spatial structure of porous media. In this paper, we take advantage of entropy method to obtain the appropriate size of the template, which is proven to be correct and feasible. The reconstruction method proposed by us combines successive 2D MPS simulations as well as 3D MPS simulation, which takes account into the pore structure information (e.g., heterogeneity and connectivity) both intralayer and interlayer. This reconstruction method is tested on Fontainebleau sandstone for which 3D images from micro‐CT scanning are available. Applying local percolation theory analysis, this new approach can depict the expected patterns of geological heterogeneities. In addition, it also can well reproduce a high degree of connectivity of the pore space better than other reconstruction methods based on lower‐order statistics. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
15.
Joseph Palis 《GeoJournal》2009,74(3):227-234
My paper aims to ask what space characterizes the various constructs of orientalism and othering in the early short films
of Thomas Edison. Using Lefebvre’s concept that social space “subsumes things produced, and encompasses their interrelationships
in their coexistence and simultaneity” in these early shorts, I will look at three Biograph actualities found at the Library
of Congress-American Memory page to show how space is manifested and negotiated onscreen. I will examine Edison’s “Filipinos
Retreat From Trenches”, “Capture of Trenches at Candaba” and “U.S. Troops and Red Cross in the Trenches Before Caloocan” which
were all released in 1899. These reenacted short films were shot during the tumultuous years of the Spanish-American War.
In the Biograph shorts, the privileged positions of both Spanish and American forces in relation to the annexation of a foreign
land in world history books is indicative of the tendency to de-emphasize the contribution of the native population in the
war. Manthia Diawara has said that “space is related to power and powerlessness, insofar as those who occupy the center of
the screen are usually more powerful than those in the background or completely absent from the screen.” The spatial hierarchies
and spatially situated images in Edison’s short films show how historically configured power relations encrypted oppression
to its external “others” through the cinematic apparatus.
相似文献
| Joseph PalisEmail: |
16.
Joaquín A. Cortés 《Mathematical Geosciences》2009,41(7):817-828
The current theoretical development of the analysis of compositional data in the article by Aitchison and Egozcue neglects
the use of Harker’s variation diagrams and other similar plots as “meaningless” or “useless” on compositional data. In this
work, it is shown that variation diagrams essentially are not a correlation tool but a graphical representation of the mass
actions and mass balances principles in the context of a given geological system, and, when they are used correctly, they
provide vital information for the igneous petrologist. The qualitative validity of the “spurious trends” in these diagrams
is also shown, when they are interpreted in their proper geological framework. The example previously used by Rollinson to
test the usefulness of the log-ratio transformation in the Aitchison and Egozcue article is revisited here in order to fully
illustrate the proper use of this tool. 相似文献
17.
Memory-Efficient Categorical Multi-point Statistics Algorithms Based on Compact Search Trees 总被引：1，自引：1，他引：0
Tuanfeng Zhang Stein Inge Pedersen Christen Knudby David McCormick 《Mathematical Geosciences》2012,44(7):863-879
Multi-point statistics (MPS) has emerged as an advanced geomodeling approach. A practical MPS algorithm named snesim (simple normal equations simulation), which uses categorical-variable training images, was proposed in 2001. The snesim algorithm generates a search tree to store the occurrence statistics of all patterns in the training image within a given set of search templates before the simulation proceeds. The snesim search tree concept makes MPS simulation central processing unit efficient but consumes large amounts of memory, particularly when three-dimensional training images contain complex patterns and when a large search template is required to ensure optimal reproduction of the image patterns. To crack the memory-restriction bottleneck, we have developed a compact search tree that contains the same information but reduces memory cost by one order of magnitude. Furthermore, the compact structure also accelerates MPS simulation significantly. Such remarkable improvement makes MPS a more practical tool to use in building the large and complex three-dimensional facies models required in the oil and gas industry. 相似文献
18.
Stefan M. Luthi 《Mathematical Geology》1994,26(2):181-196
The texture of digital rock images, as recorded, for instance, with borehole imaging devices, is shown to reflect different bedding types. Textural segmentation of borehole images, therefore, subdivides the recorded sequence into bedding units. We show that a textural segmentation algorithm based on the concept of texture energy achieves good results when compared with synthetic as well as real data in which petroleum geologists have performed zonations on cores. Texture energy involves filtering of the original image with a set of texture sensitive masks. The filtering is done as a finite convolution over the size of the masks. On the resulting images the variance is computed over a relatively large sliding window, which, in its practical implementation, covers the full width of the image. The resulting nine one-dimensional curves are then clustered hierarchically into a user-determined number of image texture or lithological bedding classes. Principal component analysis previous to clustering can be used to reduce redundancy in the data. A recurring and relatively ill-defined problem in this field are macro-textures, i.e., the cyclic interbedding of two or more bedding types. We show that sliding Fourier transforms and variable mask scale can successfully address the zonation of macro-textures. In general, the method gives best results with mask sizes equivalent to 2–4 centimeters, reflecting the length scale at which the investigated geological bedding seems to have its highest variation. 相似文献
19.
Spatially distributed and varying natural phenomena encountered in geoscience and engineering problem solving are typically
incompatible with Gaussian models, exhibiting nonlinear spatial patterns and complex, multiple-point connectivity of extreme
values. Stochastic simulation of such phenomena is historically founded on second-order spatial statistical approaches, which
are limited in their capacity to model complex spatial uncertainty. The newer multiple-point (MP) simulation framework addresses
past limits by establishing the concept of a training image, and, arguably, has its own drawbacks. An alternative to current
MP approaches is founded upon new high-order measures of spatial complexity, termed “high-order spatial cumulants.” These
are combinations of moments of statistical parameters that characterize non-Gaussian random fields and can describe complex
spatial information. Stochastic simulation of complex spatial processes is developed based on high-order spatial cumulants
in the high-dimensional space of Legendre polynomials. Starting with discrete Legendre polynomials, a set of discrete orthogonal
cumulants is introduced as a tool to characterize spatial shapes. Weighted orthonormal Legendre polynomials define the so-called
Legendre cumulants that are high-order conditional spatial cumulants inferred from training images and are combined with available
sparse data sets. Advantages of the high-order sequential simulation approach developed herein include the absence of any
distribution-related assumptions and pre- or post-processing steps. The method is shown to generate realizations of complex
spatial patterns, reproduce bimodal data distributions, data variograms, and high-order spatial cumulants of the data. In
addition, it is shown that the available hard data dominate the simulation process and have a definitive effect on the simulated
realizations, whereas the training images are only used to fill in high-order relations that cannot be inferred from data.
Compared to the MP framework, the proposed approach is data-driven and consistently reconstructs the lower-order spatial complexity
in the data used, in addition to high order. 相似文献
20.