Two-dimensional Conditional Simulations Based on the Wavelet Decomposition of Training Images

Scale dependency is a critical topic when modeling spatial phenomena of complex geological patterns that interact at different spatial scales. A two-dimensional conditional simulation based on wavelet decomposition is proposed for simulating geological patterns at different scales. The method utilizes the wavelet transform of a training image to decompose it into wavelet coefficients at different scales, and then quantifies their spatial dependence. Joint simulation of the wavelet coefficients is used together with available hard and or soft conditioning data. The conditionally co-simulated wavelet coefficients are back-transformed generating a realization of the attribute under study. Realizations generated using the proposed method reproduce the conditioning data, the wavelet coefficients and their spatial dependence. Two examples using geological images as training images elucidate the different aspects of the method, including hard and soft conditioning, the ability to reproduce some non-linear features and scale dependencies of the training images.     

Modeling complex reservoir geometries with multiple-point statistics

Large-scale reservoir architecture constitutes first-order reservoir heterogeneity and dietates to a large extent reservoir flow behavior. It also manifests geometric characteristics beyond the capability of traditional geostatistical models conditioned only on single-point and two-point statistics. Multiple-point statistics, as obtained by scanning a training image deemed representative of the actual reservoir, if reproduced properly provides stochastic models that better capture the essence of the heterogeneity. A growth algorithm, coupled with an optimization procedure, is proposed to reproduce target multiple-point histograms. The growth algorithm makes an analogy between geological accretion process and stochastic process and amounts to restricting the random path of sequential simulation at any given stage to a set of eligible nodes (immediately adjacent to a previously simulated node or sand grain). The proposed algorithm, combined with a multiple-grid approach, is shown to reproduce effectively the geometric essence of complex training images exhibiting long-range and curvilinear structures. Also, by avoiding a rigorous search for global minimum and accepting local minima, the proposed algorithm improves CPU time over traditional optimization procedures by several orders of magnitude. Average flow responses run on simulated realizations are shown to bracket correctly the reference responses of the training image.     

Conditional Simulation with Patterns

An entirely new approach to stochastic simulation is proposed through the direct simulation of patterns. Unlike pixel-based (single grid cells) or object-based stochastic simulation, pattern-based simulation simulates by pasting patterns directly onto the simulation grid. A pattern is a multi-pixel configuration identifying a meaningful entity (a puzzle piece) of the underlying spatial continuity. The methodology relies on the use of a training image from which the pattern set (database) is extracted. The use of training images is not new. The concept of a training image is extensively used in simulating Markov random fields or for sequentially simulating structures using multiple-point statistics. Both these approaches rely on extracting statistics from the training image, then reproducing these statistics in multiple stochastic realizations, at the same time conditioning to any available data. The proposed approach does not rely, explicitly, on either a statistical or probabilistic methodology. Instead, a sequential simulation method is proposed that borrows heavily from the pattern recognition literature and simulates by pasting at each visited location along a random path a pattern that is compatible with the available local data and any previously simulated patterns. This paper discusses the various implementation details to accomplish this idea. Several 2D illustrative as well as realistic and complex 3D examples are presented to showcase the versatility of the proposed algorithm.     

Multiple-point geostatistical modeling based on the cross-correlation functions

An important issue in reservoir modeling is accurate generation of complex structures. The problem is difficult because the connectivity of the flow paths must be preserved. Multiple-point geostatistics is one of the most effective methods that can model the spatial patterns of geological structures, which is based on an informative geological training image that contains the variability, connectivity, and structural properties of a reservoir. Several pixel- and pattern-based methods have been developed in the past. In particular, pattern-based algorithms have become popular due to their ability for honoring the connectivity and geological features of a reservoir. But a shortcoming of such methods is that they require a massive data base, which make them highly memory- and CPU-intensive. In this paper, we propose a novel methodology for which there is no need to construct pattern data base and small data event. A new function for the similarity of the generated pattern and the training image, based on a cross-correlation (CC) function, is proposed that can be used with both categorical and continuous training images. We combine the CC function with an overlap strategy and a new approach, adaptive recursive template splitting along a raster path, in order to develop an algorithm, which we call cross-correlation simulation (CCSIM), for generation of the realizations of a reservoir with accurate conditioning and continuity. The performance of CCSIM is tested for a variety of training images. The results, when compared with those of the previous methods, indicate significant improvement in the CPU and memory requirements.     

Multiple-Point Simulation with an Existing Reservoir Model as Training Image

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.     

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.     

Comparison of mineralization pattern of geochemical data in spatial and position-scale domain using new DWT- PCA approach

In the current research to determine the mineralization pattern and discuss the mineralization components, the information of position - scale domain of geochemical data has been analyzed. A new method is proposed based on coupling discrete wavelet transforms (DWT) and principal component analysis (PCA) for mineralization elements forecasting applications. The results of this study indicate the potential of DWT–PCA method for geochemical data processing. Wavelet transform (WT), as a multi-spectral analysis method, can decompose the spatial and temporal signals into different frequencies. The features of mineralization can be identified using the position - scale domain of geochemical data that may not be achievable in spatial domain. The geochemical data from the Dalli region have been processed in the spatial domain using PCA. The surface geochemical data of 30 elements have been transformed to position–scale domain using two-dimensional discrete wavelet transform (2DDWT). Wavelet functions (WFs) of Haar, Coiflet2, Biorthogonal3.3 and Symlet7 have been applied separately to decompose the geochemical data to high and low frequencies in one level. To obtain more accurate and complete information of mineralization, a new index has been presented based on wavelet coefficients. Based on this new index, significant results have been obtained by using PCA of the index. The coefficients distribution map (CDM) as a new exploratory criterion has been generated based on 2DDWT to show the geochemical distribution map (GDM). Finally, the results of WT have been compared with the results of spatial domain and the best method of wavelet for interpretation of geochemical data has been introduced. The results of geochemical data analysis by DWT–PCA approach have been confirmed by the exploratory drillings in the study area.     

High-order Stochastic Simulation of Complex Spatially Distributed Natural Phenomena

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.     

Sub-pixel Mapping of Coarse Satellite Remote Sensing Images with Stochastic Simulations from Training Images

Super-resolution or sub-pixel mapping is the process of providing fine scale land cover maps from coarse-scale satellite sensor information. Such a procedure calls for a prior model depicting the spatial structures of the land cover types. When available, an analog of the underlying scene (a training image) may be used for such a model. The single normal equation simulation algorithm (SNESIM) allows extracting the relevant pattern information from the training image and uses that information to downscale the coarse fraction data into a simulated fine scale land cover scene. Two non-exclusive approaches are considered to use training images for super-resolution mapping. The first one downscales the coarse fractions into fine-scale pre-posterior probabilities which is then merged with a probability lifted from the training image. The second approach pre-classifies the fine scale patterns of the training image into a few partition classes based on their coarse fractions. All patterns within a partition class are recorded by a search tree; there is one tree per partition class. At each fine scale pixel along the simulation path, the coarse fraction data is retrieved first and used to select the appropriate search tree. That search tree contains the patterns relevant to that coarse fraction data. To ensure exact reproduction of the coarse fractions, a servo-system keeps track of the number of simulated classes inside each coarse fraction. Being an under-determined stochastic inverse problem, one can generate several super resolution maps and explore the space of uncertainty for the fine scale land cover. The proposed SNESIM sub-pixel resolution mapping algorithms allow to: (i) exactly reproduce the coarse fraction, (ii) inject the structural model carried by the training image, and (iii) condition to any available fine scale ground observations. Two case studies are provided to illustrate the proposed methodology using Landsat TM data from southeast China.     

Stationarity Scores on Training Images for Multipoint Geostatistics

This research introduces a novel method to assess the validity of training images used as an input for Multipoint Geostatistics, alternatively called Multiple Point Simulation (MPS). MPS are a family of spatial statistical interpolation algorithms that are used to generate conditional simulations of property fields such as geological facies. They are able to honor absolute “hard” constraints (e.g., borehole data) as well as “soft” constraints (e.g., probability fields derived from seismic data, and rotation and scale). These algorithms require 2D or 3D training images or analogs whose textures represent a spatial arrangement of geological properties that is presumed to be similar to that of a target volume to be modeled. To use the current generation of MPS algorithms, statistically valid training image are required as input. In this context, “statistical validity” includes a requirement of stationarity, so that one can derive from the training image an average template pattern. This research focuses on a practical method to assess stationarity requirements for MPS algorithms, i.e., that statistical density or probability distribution of the quantity shown on the image does not change spatially, and that the image shows repetitive shapes whose orientation and scale are spatially constant. This method employs image-processing techniques based on measures of stationarity of the category distribution, the directional (or orientation) property field and the scale property field of those images. It was successfully tested on a set of two-dimensional images representing geological features and its predictions were compared to actual realizations of MPS algorithms. An extension of the algorithms to 3D images is also proposed. As MPS algorithms are being used increasingly in hydrocarbon reservoir modeling, the methods described should facilitate screening and selection of the input training images.     

A Multiple Training Image Approach for Spatial Modeling of Geologic Domains

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.     

A Distance-based Prior Model Parameterization for Constraining Solutions of Spatial Inverse Problems

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.     

基于二维地质剖面的三维地质结构多点统计学模拟方法

利用三维地质模拟技术重构地质现象的三维空间分布,是实现自然资源管理和风险评估的重要基础和前提。多点统计学方法通过探寻多点间的空间结构关系,结合随机模拟方法生成具有差异性的模拟结果,较好地再现了复杂的地质现象。然而,如何构建合适、有效的训练图像一直是基于多点统计学三维地质模拟的核心问题。本文提出了一种改进的多点统计学算法。本方法结合了序贯模拟和迭代的方法,将二维剖面扩展为三维训练图像,再结合EM-Like算法,实现了三维地质结构的优化模拟。建模实例结果表明,本方法能确保训练图像对内部模拟网格的约束,准确模拟研究区的地层层序,并很好地再现二维地质剖面所反映的地层结构关系。     

Stochastic computational modelling of highly heterogeneous poroelastic media with long‐range correlations

The compaction of highly heterogeneous poroelastic reservoirs with the geology characterized by long‐range correlations displaying fractal character is investigated within the framework of the stochastic computational modelling. The influence of reservoir heterogeneity upon the magnitude of the stresses induced in the porous matrix during fluid withdrawal and rock consolidation is analysed by performing ensemble averages over realizations of a log‐normally distributed stationary random hydraulic conductivity field. Considering the statistical distribution of this parameter characterized by a coefficient of variation governing the magnitude of heterogeneity and a correlation function which decays with a power‐law scaling behaviour we show that the combination of these two effects result in an increase in the magnitude of effective stresses of the rock during reservoir depletion. Further, within the framework of a perturbation analysis we show that the randomness in the hydraulic conductivity gives rise to non‐linear corrections in the upscaled poroelastic equations. These corrections are illustrated by a self‐consistent recursive hierarchy of solutions of the stochastic poroelastic equations parametrized by a scale parameter representing the fluctuating log‐conductivity standard deviation. A classical example of land subsidence caused by fluid extraction of a weak reservoir is numerically simulated by performing Monte Carlo simulations in conjunction with finite elements discretizations of the poroelastic equations associated with an ensemble of geologies. Numerical results illustrate the effects of the spatial variability and fractal character of the permeability distribution upon the evolution of the Mohr–Coulomb function of the rock. Copyright © 2004 John Wiley & Sons, Ltd.     

Dimensional Reduction of Pattern-Based Simulation Using Wavelet Analysis

A pattern-based simulation technique using wavelet analysis is proposed for the simulation (wavesim) of categorical and continuous variables. Patterns are extracted by scanning a training image with a template and then storing them in a pattern database. The dimension reduction of patterns in the pattern database is performed by wavelet decomposition at certain scale and the approximate sub-band is used for pattern database classification. The pattern database classification is performed by the k-means clustering algorithm and classes are represented by a class prototype. For the simulation of categorical variables, the conditional cumulative density function (ccdf) for each class is generated based on the frequency of the individual categories at the central node of the template. During the simulation process, the similarity of the conditioning data event with the class prototypes is measured using the L ₂-norm. When simulating categorical variables, the ccdf of the best matched class is used to draw a pattern from a class. When continuous variables are simulated, a random pattern is drawn from the best matched class. Several examples of conditional and unconditional simulation with two- and three- dimensional data sets show that the spatial continuity of geometric features and shapes is well reproduced. A comparative study with the filtersim algorithm shows that the wavesim performs better than filtersim in all examples. A full-field case study at the Olympic Dam base metals deposit, South Australia, simulates the lithological rock-type units as categorical variables. Results show that the proportions of various rock-type units in the hard data are well reproduced when similar to those in the training image; when rock-type proportions between the training image and hard data differ, the results show a compromise between the two.     

基于地质矢量信息的冲积扇储层沉积微相建模： 以克拉玛依油田三叠系克下组为例

冲积扇储层具有典型的非平稳分布特征,传统的基于平稳假设的地质统计建模方法无法准确模拟其沉积微相。在充 分分析其沉积微相分布特征的基础上,以克拉玛依油田三叠系克下组冲积扇储层为例,综合现代沉积、密井网、露头等信 息,应用基于地质矢量信息的多点地质统计学方法建立了精细、合理的储层沉积微相三维模型。首先,综合现代沉积特征 与地下储层沉积微相解剖成果,应用“模式指导、规模约束”的方法,建立符合研究区沉积模式和定量规模的训练图像和 基于冲积扇沉积体系的矢量坐标系统,表征了训练图像的地质矢量信息,提出基于预模拟实现的模拟域地质矢量信息表征 方法,建立了待模拟地质体的三维矢量信息模型;最后,采用基于矢量信息的多点地质统计方法建立了冲积扇储层沉积微 相模型,并采用多种标准对模拟实现进行可靠性分析。模拟结果表明冲积扇相带分布符合地质认识,沉积微相的形态、展 布特征、组合样式、规模与训练图像、地下沉积微相解剖认识相符。     

Integration of Markov mesh models and data assimilation techniques in complex reservoirs

We present a methodology that allows conditioning the spatial distribution of geological and petrophysical properties of reservoir model realizations on available production data. The approach is fully consistent with modern concepts depicting natural reservoirs as composite media where the distribution of both lithological units (or facies) and associated attributes are modeled as stochastic processes of space. We represent the uncertain spatial distribution of the facies through a Markov mesh (MM) model, which allows describing complex and detailed facies geometries in a rigorous Bayesian framework. The latter is then embedded within a history matching workflow based on an iterative form of the ensemble Kalman filter (EnKF). We test the proposed methodology by way of a synthetic study characterized by the presence of two distinct facies. We analyze the accuracy and computational efficiency of our algorithm and its ability with respect to the standard EnKF to properly estimate model parameters and assess future reservoir production. We show the feasibility of integrating MM in a data assimilation scheme. Our methodology is conducive to a set of updated model realizations characterized by a realistic spatial distribution of facies and their log permeabilities. Model realizations updated through our proposed algorithm correctly capture the production dynamics.     

多点地质统计学建模的发展趋势

从算法研究、训练图像处理和实际应用三个方面详细解剖了国内外多点地质统计学的发展历程,在此基础上,分析了多点地质统计学主流的几种算法的核心原理、适用范围及优缺点,以此来对储层建模的发展趋势作出展望。目前,多点地质统计学虽是随机建模的一种前沿研究热点,但由于其尚未成熟,仍需对建模算法进行研究。为此,在前人研究的基础上,重点分析了多点地质统计学的发展趋势：合理处理训练图像;合理利用软信息;选择合适的相似性方法;选择合适的标准化方法;合理利用平稳性;算法间的耦合;选择合适的过滤器;拓展缝洞型碳酸盐岩模拟。最后,提出多点地质统计学在储层建模方面,应从增加储层的模拟区域、提高模拟精度、扩大储层相的模拟范围和提高计算机模拟效率等方面进行改进。     

Representing Spatial Uncertainty Using Distances and Kernels

Assessing uncertainty of a spatial phenomenon requires the analysis of a large number of parameters which must be processed by a transfer function. To capture the possibly of a wide range of uncertainty in the transfer function response, a large set of geostatistical model realizations needs to be processed. Stochastic spatial simulation can rapidly provide multiple, equally probable realizations. However, since the transfer function is often computationally demanding, only a small number of models can be evaluated in practice, and are usually selected through a ranking procedure. Traditional ranking techniques for selection of probabilistic ranges of response (P10, P50 and P90) are highly dependent on the static property used. In this paper, we propose to parameterize the spatial uncertainty represented by a large set of geostatistical realizations through a distance function measuring “dissimilarity” between any two geostatistical realizations. The distance function allows a mapping of the space of uncertainty. The distance can be tailored to the particular problem. The multi-dimensional space of uncertainty can be modeled using kernel techniques, such as kernel principal component analysis (KPCA) or kernel clustering. These tools allow for the selection of a subset of representative realizations containing similar properties to the larger set. Without losing accuracy, decisions and strategies can then be performed applying a transfer function on the subset without the need to exhaustively evaluate each realization. This method is applied to a synthetic oil reservoir, where spatial uncertainty of channel facies is modeled through multiple realizations generated using a multi-point geostatistical algorithm and several training images.     

Facies Modeling Using a Markov Mesh Model Specification

The spatial continuity of facies is one of the key factors controlling flow in reservoir models. Traditional pixel-based methods such as truncated Gaussian random fields and indicator simulation are based on only two-point statistics, which is insufficient to capture complex facies structures. Current methods for multi-point statistics either lack a consistent statistical model specification or are too computer intensive to be applicable. We propose a Markov mesh model based on generalized linear models for geological facies modeling. The approach defines a consistent statistical model that is facilitated by efficient estimation of model parameters and generation of realizations. Our presentation includes a formulation of the general framework, model specifications in two and three dimensions, and details on how the parameters can be estimated from a training image. We illustrate the method using multiple training images, including binary and trinary images and simulations in two and three dimensions. We also do a thorough comparison to the snesim approach. We find that the current model formulation is applicable for multiple training images and compares favorably to the snesim approach in our test examples. The method is highly memory efficient.