Direct Multiple-Point Geostatistical Simulation of Edge Properties for Modeling Thin Irregularly Shaped Surfaces

Thin, irregularly shaped surfaces such as clay drapes often have a major control on flow and transport in heterogeneous porous media. Clay drapes are often complex, curvilinear three-dimensional surfaces and display a very complex spatial distribution. Variogram-based stochastic approaches are also often not able to describe the spatial distribution of clay drapes since complex, curvilinear, continuous, and interconnected structures cannot be characterized using only two-point statistics. Multiple-point geostatistics aims to overcome the limitations of the variogram. The premise of multiple-point geostatistics is to move beyond two-point correlations between variables and to obtain (cross) correlation moments at three or more locations at a time using training images to characterize the patterns of geological heterogeneity. Multiple-point geostatistics can reproduce thin irregularly shaped surfaces such as clay drapes, but this is often computationally very intensive. This paper describes and applies a methodology to simulate thin, irregularly shaped surfaces with a smaller CPU and RAM demand than the conventional multiple-point statistical methods. The proposed method uses edge properties for indicating the presence of thin irregularly shaped surfaces. Instead of pixel values, edge properties indicating the presence of irregularly shaped surfaces are simulated using snesim. This method allows direct simulation of edge properties instead of pixel properties to make it possible to perform multiple-point geostatistical simulations with a larger cell size and thus a smaller computation time and memory demand. This method is particularly valuable for three-dimensional applications of multiple-point geostatistics.     

Application of multiple-point geostatistics on modelling groundwater flow and transport in a cross-bedded aquifer (Belgium)

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.     

多点地质统计学在苏里格气田地质建模的应用

多点地质统计学是一种建立地质模型的统计学方法,该方法比传统的两点地质统计学更适合河流相沉积体系。本文综合地震、测井、录井及生产资料,绘制了长庆油田苏里格气田苏48区块盒8下段辫状河训练图像,并在此基础上利用多点地质统计学方法,加入三维地震资料作为约束,以水平井整体开发为研究对象,建立了该区的地质模型,优化了水平井整体部署,指导水平井导向。利用多口水平井实钻效果验证了地质模型的精度,以苏19-62井为例得到:录井显示气层钻遇率为78.3%,测井显示气层钻遇率为69%。这与地震反演剖面预测相似,且与实钻效果有很好的对应性。     

在储层建模中利用多点地质统计学整合地质概念模型及其解释(英文)

基于三维空间中稀疏的观测数据,地质学家和储层建模人员尝试预测井间的地质沉积相的空间非均质性时,地质概念模型和先验认识在其中扮演着重要的角色。这种整合先验模型或解释的过程有时是隐蔽或不易察觉的,正如在手工绘等值线图中的情形;它也能够被显式地运用到某种算法当中,比如数字绘图中的算法。新近兴起的多点地质统计学为地质学家和储层建模人员提供了一种有力工具,它强调使用训练图像把先验模型明确而定量地引入到储层建模当中。先验地质模型包含了被研究的真实储层中确信存在的样式,而训练图像则是该模型的定量化表达。通过再现高阶统计量,多点算法能够从训练图像中捕捉复杂的(非线性)特征样式并把它们锚定到观测的井位数据。文中描述了多点地质统计学原理,以突出训练图像概念重要性为主线,描述了多点地质统计学在建立三维储层模型中的应用。     

Projection Pursuit Multivariate Transform

Transforming complex multivariate geological data to a Gaussian distribution is an important and challenging problem in geostatistics. A variety of transforms are available for this goal, but struggle with high dimensional data sets. Projection pursuit density estimation (PPDE) is a well-established nonparametric method for estimating the joint density of multivariate data. A central component of the PPDE algorithm transforms the original data toward a multivariate Gaussian distribution. The PPDE approach is modified to map complex data to a multivariate Gaussian distribution within a geostatistical modeling context. Traditional modeling may then take place on the transformed Gaussian data, with a back-transform used to return simulated variables to their original units. This approach is referred to as the projection pursuit multivariate transform (PPMT). The PPMT shows the potential to be an effective means for modeling high dimensional and complex geologic data. The PPMT algorithm is developed before discussing considerations and limitations. A case study compares modeling results against more common techniques to demonstrate the value and place of the PPMT within geostatistics.     

A multiple-point geostatistical method based on geological vector information

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.     

Stochastic analysis of the effect of spatial variability of diffusion parameters on radionuclide transport in a low permeability clay layer

Most studies that incorporate subsurface heterogeneity in groundwater flow and transport models only analyze and simulate the spatial variability of hydraulic conductivity. Heterogeneity of the other flow and transport parameters are usually neglected. This approach is often justified, but there are, however, cases in which disregarding the heterogeneity of the other flow and transport parameters can be questionable. In low permeability media, for instance, diffusion is often the dominant transport mechanism. It therefore seems logical to incorporate the spatial variability of the diffusion parameters in the transport model. This study therefore analyses and simulates the spatial variability of the effective diffusion coefficient and the diffusion accessible porosity with geostatistical techniques and incorporates their heterogeneity in the transport model of a low permeability formation. The formation studied was Boom clay (Belgium), a candidate host rock for the deep geological disposal of high-level radioactive waste. The calculated output radionuclide fluxes of this model are compared with the fluxes calculated with a homogeneous model and a model with a heterogeneous hydraulic conductivity distribution. This analysis shows that the heterogeneity of the diffusion parameters has a much larger effect on the calculated output radionuclide fluxes than the heterogeneity of hydraulic conductivity in the low permeability medium under study.     

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.     

Challenges in reservoir forecasting

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.     

Hierarchical simulation of aquifer heterogeneity: implications of different simulation settings on solute-transport modeling

The fine-scale heterogeneity of porous media affects the large-scale transport of solutes and contaminants in groundwater and it can be reproduced by means of several geostatistical simulation tools. However, including the available geological information in these tools is often cumbersome. A hierarchical simulation procedure based on a binary tree is proposed and tested on two real-world blocks of alluvial sediments, of a few cubic meters volume, that represent small-scale aquifer analogs. The procedure is implemented using the sequential indicator simulation, but it is so general that it can be adapted to various geostatistical simulation tools, improving their capability to incorporate geological information, i.e., the sedimentological and architectural characterization of heterogeneity. When compared with a standard sequential indicator approach on bi-dimensional simulations, in terms of proportions and connectivity indicators, the proposed procedure yields reliable results, closer to the reference observations. Different ensembles of three-dimensional simulations based on different hierarchical sequences are used to perform numerical experiments of conservative solute transport and to obtain ensembles of equivalent pore velocity and dispersion coefficient at the scale length of the blocks (meter). Their statistics are used to estimate the impact of the variability of the transport properties of the simulated blocks on contaminant transport modeled on bigger domains (hectometer). This is investigated with a one-dimensional transport modeling based on the Kolmogorov-Dmitriev theory of branching stochastic processes. Applying the proposed approach with diverse binary trees and different simulation settings provides a great flexibility, which is revealed by the differences in the breakthrough curves.     

非参数地质统计学进展

地质统计学从解决固体矿产资源评价起源,应用领域迅速扩展。文章简要介绍了地质统计学的理论研究现状和应用现状。为了减轻人为因素在创建模型中的作用,非参数的地质统计学得到发展,文章重点介绍了非参数地质统计学的4种变异函数模型。此外,文章阐述了多点地质统计学的概念,这是一种用训练映像代替变异函数的地质统计学方法,使得地质现象的解释更直观。     

Nodal and Mixed Finite Elements for the Numerical Homogenization of 3D Permeability

The aim of upscaling is to determine equivalent homogeneous parameters at a coarse-scale from a spatially oscillating fine-scale parameter distribution. To be able to use a limited number of relatively large grid-blocks in numerical oil reservoir simulators or groundwater models, upscaling of the permeability is frequently applied. The spatial fine-scale permeability distribution is generally obtained from geological and geostatistical models. After upscaling, the coarse-scale permeabilities are incorporated in the relatively large grid-blocks of the numerical model. If the porous rock may be approximated as a periodic medium, upscaling can be performed by the method of homogenization. In this paper the homogenization is performed numerically, which gives rise to an approximation error. The complementarity between two different numerical methods – the conformal-nodal finite element method and the mixed-hybrid finite element method – has been used to quantify this error. These two methods yield respectively upper and lower bounds for the eigenvalues of the coarse-scale permeability tensor. Results of 3D numerical experiments are shown, both for the far field and around wells.     

多点地质统计学建模研究进展

从精细油藏描述中地质建模的意义和现状入手,介绍了多点地质统计学建模研究现状及其与传统地质建模方法的差 异。以辽河盆地西部凹陷某蒸汽驱试验区为例,分析了多点地质统计学建模中建模的基础、训练图像的建立、多点地质统计学 建模与传统地质建模相比所具有的优势等内容。指出多点地质统计学在井间预测方面具有明显优于其他传统建模方法的特 点。在文献调研基础上,结合自身工作实践,探讨了多点地质统计学建模目前存在的问题和未来的发展方向。指出未来多点 地质统计学建模的发展方向主要包括多信息综合地质成因分析基础上的训练图像获取、多点地质统计学算法进行改进和完善 和多点地质统计学建模方法应用领域的扩大等。     

Conditioning of Multiple-Point Statistics Facies Simulations to Tomographic Images

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.     

Upscaling Hydraulic Conductivities in Cross-Bedded Formations

Modern geostatistical techniques allow the generation of high-resolution heterogeneous models of hydraulic conductivity containing millions to billions of cells. Selective upscaling is a numerical approach for the change of scale of fine-scale hydraulic conductivity models into coarser scale models that are suitable for numerical simulations of groundwater flow and mass transport. Selective upscaling uses an elastic gridding technique to selectively determine the geometry of the coarse grid by an iterative procedure. The geometry of the coarse grid is built so that the variances of flow velocities within the coarse blocks are minimum. Selective upscaling is able to handle complex geological formations and flow patterns, and provides full hydraulic conductivity tensor for each block. Selective upscaling is applied to a cross-bedded formation in which the fine-scale hydraulic conductivities are full tensors with principal directions not parallel to the statistical anisotropy of their spatial distribution. Mass transport results from three coarse-scale models constructed by different upscaling techniques are compared to the fine-scale results for different flow conditions. Selective upscaling provides coarse grids in which mass transport simulation is in good agreement with the fine-scale simulations, and consistently superior to simulations on traditional regular (equal-sized) grids or elastic grids built without accounting for flow velocities.     

Conditional Simulation of Complex Geological Structures Using Multiple-Point Statistics

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.     

Integration of Tracer Test Data to Refine Geostatistical Hydraulic Conductivity Fields Using Sequential Self-Calibration Method

On the basis of local measurements of hydraulic conductivity,geostatistical methods have been found to be useful in heterogeneity characterization of a hydraulic conductivity field on a regional scale. However,the methods are not suited to directly integrate dynamic production data,such as,hydraulic head and solute concentration,into the study of conductivity distribution. These data,which record the flow and transport processes in the medium,are closely related to the spatial distribution of hydraulic conductivity. In this study,a three-dimensional gradient-based inverse method-the sequential self-calibration (SSC) method-is developed to calibrate a hydraulic conductivity field,initially generated by a geostatistical simulation method,conditioned on tracer test results. The SSC method can honor both local hydraulic conductivity measurements and tracer test data. The mismatch between the simulated hydraulic conductivity field and the reference true one,measured by its mean square error (MSE),is reduced through the SSC conditional study. In comparison with the unconditional results,the SSC conditional study creates the mean breakthrough curve much closer to the reference true curve,and significantly reduces the prediction uncertainty of the solute transport in the observed locations. Further,the reduction of uncertainty is spatially dependent,which indicates that good locations,geological structure,and boundary conditions will affect the efficiency of the SSC study results.     

A Comparison of Methods for the Stochastic Simulation of Rock Fractures

Methods reported in the literature for rock fracture simulations include approaches based on stochastic geometry, multiple-point statistics and a combination of geostatistics for fracture density and object-based modelling for fracture geometries. The advantages and disadvantages of each of these approaches are discussed with examples. By way of review, the authors begin with the geostatistical indicator simulation method, based on the truncated–Gaussian algorithm; this is followed by multiple-point statistical simulation and then the stochastic geometry approach, which is based on marked point process simulation. A new approach, based on pluriGaussian structural simulation, is then introduced. The new approach incorporates in the simulation the spatial correlation between different sets of fractures, which in general, is very difficult, if not impossible, to accomplish in the three methods reviewed. Each simulation method is summarised together with detailed simulation procedures for each. A published two-dimensional fracture dataset is used as a means of assessing the performance of each simulation method and of demonstrating the concepts discussed in the text.     

Improving Sequential Simulation with a Structured Path Guided by Information Content

Multiple-point simulation is a newly developed geostatistical method that aims at combining the strengths of two mainstream geostatistical methods: object-based and pixel-based methods. It maintains the flexibility of pixel-based algorithms in data conditioning, while enhancing its capability of reproducing realistic geological shapes, which is traditionally reserved to object-based algorithms. However, the current snesim program for multiple-point simulation has difficulty in reproducing large-scale structures, which have a significant impact on the flow response. To address this problem, we propose to simulate along a structured path based on an information content measure. This structured path accounts for not only the information from the data, but also some prior structural information provided by geological knowledge. Various case studies show a better reproduction of large-scale structures. This concept of simulating along a structured path guided by information content can be applied to any sequential simulation algorithms, including traditional variogram-based two-point geostatistical algorithms.     

High-order Statistics of Spatial Random Fields: Exploring Spatial Cumulants for Modeling Complex Non-Gaussian and Non-linear Phenomena

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.