Abstract
Abstract
Two approaches are traditionally used to build numerical models for facies distributions within a reservoir. Pixel-based techniques aim at generating simulated realizations that honor the well data values, and reproduce a given variogram which models two-point spatial correlation. However, because the variogram cannot look at spatial continuity between more than two locations at a time, pixel-based algorithms give poor representations of the actual facies geometries. In contrast, object-based techniques allow reproducing crisp geometries, but the conditioning on well data requires iterative "trial-and-error" corrections, which can be time-consuming, particularly when the data are dense with regard to the average object size. This paper presents a new approach that combines the easy conditioning of pixel-based algorithms with the ability to reproduce "shapes" of object-based techniques, without being too time and memory demanding.
In this new approach, the complex geological structures expected to be present in the reservoir are characterized by multiple-point statistics, which express joint variability at many more than two locations at a time. Such multiple-point statistics cannot be inferred from typically sparse well data but could be read from training images depicting the expected subsurface heterogeneities. A training image need not carry any locally accurate information on the reservoir; it need only reflect a prior stationary geological/structural concept. Thus training images can be generated by object-based algorithms freed of the constraint of data conditioning. The multiple-point statistics inferred from the training image(s) are then exported to the reservoir model, where they are anchored to the well data using a pixel-based sequential simulation algorithm.
This algorithm is tested for the simulation of a turbidite system where flow is controlled by meandering channels with cross-bedding. The training image reflecting the channel patterns is an unconditional realization generated by an object-based algorithm. The final simulated numerical models reproduce these channel patterns, and honor exactly all well data values at their locations. The methodology proposed appears to be practical, general, and fast.
Introduction
Most current stochastic reservoir simulation algorithms aim at reproducing statistics inferred from the well data solely. Since these data are typically sparse they can give, at best, only an estimate of the two-point correlation within the reservoir. Hence most traditional simulation techniques, such as sequential Gaussian simulation1, are limited to the reproduction of a variogram model or covariance function inferred from a few sample two-point correlation moments. Identification of two-point statistics is, however, not sufficient to reproduce connected patterns over long ranges, nor clusters of similar values with a characteristic shape (e.g. meandering channels), features which are common occurrences in oil reservoirs. The modeling of such specific patterns calls for characterizing the spatial continuity at three or more locations at a time. Reproduction of multiple-point continuity in the reservoir model is critical, not so much to produce geologically realistic looking maps, but to provide accurate flow performance predictions2.
The quest for more realistic models is motivated by the fact that structural information much beyond the sample two-point correlation calculated from the well data is actually available. Stratigraphic interpretation from well logs, seismic data analysis, and comparison with reservoir analogues do provide information about the type of geological bodies present in the reservoir. Variogram-based algorithms mostly ignore such valuable information.
One straightforward way to reproduce geological bodies with crisp shapes consists in generating objects with similar shapes then distribute these objects over the reservoir model3,4,5,6,7. There are, however, three major limitations to such object-based approach:Each different geometric shape requires its own parameterization, which means that each object-based algorithm can handle only one type of reservoir.