Research output: Chapter in Book/Report/Conference proceeding › Conference contribution › Research › peer-review
Bayesian inference of covariance parameters in spectral approach to geostatistical simulation. / Ismagilov, N.; Azangulov, I.; Borovitskiy, V.; Lifshits, M.; Mostowsky, P.
ECMOR 2020 - 17th European Conference on the Mathematics of Oil Recovery. European Association of Geoscientists and Engineers, 2020. p. 1-10 (Conference Proceedings; Vol. 2020).Research output: Chapter in Book/Report/Conference proceeding › Conference contribution › Research › peer-review
}
TY - GEN
T1 - Bayesian inference of covariance parameters in spectral approach to geostatistical simulation
AU - Ismagilov, N.
AU - Azangulov, I.
AU - Borovitskiy, V.
AU - Lifshits, M.
AU - Mostowsky, P.
N1 - Publisher Copyright: Copyright © ECMOR 2020. All rights reserved. Copyright: Copyright 2021 Elsevier B.V., All rights reserved.
PY - 2020
Y1 - 2020
N2 - The spectral simulation approach (described in Ismagilov and Lifshits (ECMOR XVI)) is a relatively new geostatistical method of stochastic reservoir property simulation. It is based on Fourier analysis of well log data and simulation of Fourier expansion coefficients in the interwell space. The key advantage of this method is its ability to automatically recognize and reproduce vertical non-stationarities observed in well data (Ismagilov et al. (ATCE 2019)). This comes at a price of having many parameters: while usual geostatistical approaches like kriging or sequential Gaussian simulation require estimating one covariance function or variogram (in practice, estimated parameters are variogram model type and ranges in three directions), the spectral approach requires estimating a lot of them (typically, 100-200 covariance functions). Obviously, automatic covariance estimation becomes crucial in this setting. While assuming parametric models for the aforementioned covariance functions and estimating their parameters by maximizing the likelihood works reasonably well in practice, this strategy has some drawbacks. First, in cases when likelihood surface turns out to be multi-modal or flat, the point estimation of parameters may lead to many problems such as incorrect uncertainty estimation or even choosing a wrong model. Second, maximum likelihood estimation usually does not provide a way to incorporate prior knowledge about parameters - a typical example is constraining the resulting variogram range to be in geologically reasonable limits. We argue that Bayesian inference of parameters is a way to overcome both these challenges. Treating covariance parameters as random variables avoids limitations of deterministic point estimations while introducing prior distributions for parameters is the most natural way of incorporating the prior knowledge. We develop and implement in software a version of spectral approach where covariance parameters are treated in a Bayesian way. We show via computations on practical examples that Bayesian inference enables to build better models in cases with complex likelihood surfaces and to account for prior knowledge about covariance parameters.
AB - The spectral simulation approach (described in Ismagilov and Lifshits (ECMOR XVI)) is a relatively new geostatistical method of stochastic reservoir property simulation. It is based on Fourier analysis of well log data and simulation of Fourier expansion coefficients in the interwell space. The key advantage of this method is its ability to automatically recognize and reproduce vertical non-stationarities observed in well data (Ismagilov et al. (ATCE 2019)). This comes at a price of having many parameters: while usual geostatistical approaches like kriging or sequential Gaussian simulation require estimating one covariance function or variogram (in practice, estimated parameters are variogram model type and ranges in three directions), the spectral approach requires estimating a lot of them (typically, 100-200 covariance functions). Obviously, automatic covariance estimation becomes crucial in this setting. While assuming parametric models for the aforementioned covariance functions and estimating their parameters by maximizing the likelihood works reasonably well in practice, this strategy has some drawbacks. First, in cases when likelihood surface turns out to be multi-modal or flat, the point estimation of parameters may lead to many problems such as incorrect uncertainty estimation or even choosing a wrong model. Second, maximum likelihood estimation usually does not provide a way to incorporate prior knowledge about parameters - a typical example is constraining the resulting variogram range to be in geologically reasonable limits. We argue that Bayesian inference of parameters is a way to overcome both these challenges. Treating covariance parameters as random variables avoids limitations of deterministic point estimations while introducing prior distributions for parameters is the most natural way of incorporating the prior knowledge. We develop and implement in software a version of spectral approach where covariance parameters are treated in a Bayesian way. We show via computations on practical examples that Bayesian inference enables to build better models in cases with complex likelihood surfaces and to account for prior knowledge about covariance parameters.
UR - http://www.scopus.com/inward/record.url?scp=85099543837&partnerID=8YFLogxK
U2 - 10.3997/2214-4609.202035092
DO - 10.3997/2214-4609.202035092
M3 - Conference contribution
AN - SCOPUS:85099543837
T3 - Conference Proceedings
SP - 1
EP - 10
BT - ECMOR 2020 - 17th European Conference on the Mathematics of Oil Recovery
PB - European Association of Geoscientists and Engineers
T2 - 17th European Conference on the Mathematics of Oil Recovery, ECMOR 2020
Y2 - 14 September 2020 through 17 September 2020
ER -
ID: 75053297