Spectral simulation has gained application in building geologic models due to the advantage of better honoring the spatial continuity of petrophysical properties, such as reservoir porosity and shale volume. Distinct from sequential simulation methods, spectral simulation is a global algorithm in the sense that a global density spectrum is calculated once and the inverse Fourier transform is performed on the Fourier coefficient also only once to generate a simulation realization. The generated realizations honor the spatial continuity structure globally over the whole field instead of only within a search neighborhood, as with sequential simulation algorithms. However, the disadvantage of global spectral simulation is that it traditionally cannot account for the local information such as the local continuity trends, which are often observed in reservoirs and hence are important to be accounted for in geologic models. This disadvantage has limited wider application of spectral simulation in building geologic models. In this paper, we present ways of conditioning geologic models to the relevant local information. To account for the local continuity trends, we first scale different frequency components of the original model with local-amplitude spectrum ratios that are specific to the local trend. The sum of these scaled frequency components renders a new model that displays the desired local continuity trend. The implementation details of this new method are discussed and examples are provided to illustrate the algorithm.
[1]
Wenlong Xu,et al.
Conditional curvilinear stochastic simulation using pixel-based algorithms
,
1996
.
[2]
Tingting Yao,et al.
Conditional Spectral Simulation with Phase Identification
,
1998
.
[3]
Dulal C. Kar,et al.
On the prime factor decomposition algorithm for the discrete sine transform
,
1994,
IEEE Trans. Signal Process..
[4]
Tingting Yao.
Reproduction of the Mean, Variance, and Variogram Model in Spectral Simulation
,
2004
.
[5]
Andre G. Journel,et al.
Automatic Modeling of (Cross) Covariance Tables Using Fast Fourier Transform
,
1998
.
[6]
Leon E. Borgman,et al.
Three-Dimensional, Frequency-Domain Simulations of Geological Variables
,
1984
.
[7]
Hong Yan,et al.
Computing the inverse DFT with the in-place, in-order prime factor FFT algorithm
,
1995,
IEEE Trans. Signal Process..
[8]
M. Chica-Olmo,et al.
The Fourier Integral Method: An efficient spectral method for simulation of random fields
,
1993
.
[9]
A. Journel,et al.
Geostatistics for natural resources characterization
,
1984
.
[10]
Javier D. Bruguera,et al.
Implementation of the FFT butterfly with redundant arithmetic
,
1996
.
[11]
Tinting Yao,et al.
Specsim: a Fortran-77 program for conditional spectral simulation in 3D
,
1998
.