Runhai Feng, Klaus Mosegaard, Dario Grana, Tapan Mukerji, Thomas Mejer Hansen
{"title":"Stochastic Facies Inversion with Prior Sampling by Conditional Generative Adversarial Networks Based on Training Image","authors":"Runhai Feng, Klaus Mosegaard, Dario Grana, Tapan Mukerji, Thomas Mejer Hansen","doi":"10.1007/s11004-023-10119-0","DOIUrl":null,"url":null,"abstract":"<p>Probabilistic methods for geophysical inverse problems allow the use of arbitrarily complex prior information in principle. Geostatistical techniques, such as multiple-point statistics (MPS), for describing spatial correlation models and higher-order statistics have been proposed to achieve this inversion task, in which stochastic algorithms such as Markov chain Monte Carlo (McMC) are incorporated. However, stochastic sampling and optimization often require a large number of iterations, and thus geostatistical sampling of the prior model can become computationally demanding. To overcome this challenge, a deep learning model, namely conditional generative adversarial networks (CGANs), is proposed, which allows one to perform a random walk to sample the complex prior distribution. CGANs simulate conditional realizations conditioned to the available hard conditioning data, that is, direct measurements, while preserving the geometrical structure of the model parameters of interest and replicating the sequential Gibbs sampling algorithm. Despite the need for a training step, for a large number of simulations, CGANs are more efficient than traditional geostatistical simulation algorithms such as single normal equation simulation (SNESIM). The proposed methodology is used as part of the extended Metropolis algorithm to predict the distributions of categorical facies in two examples, a dune environment in the Gobi Desert and a channel system in an idealized subsurface reservoir, from indirect observational data such as acoustic impedance. The inversion results are compared to the extended Metropolis algorithm using standard MPS sampling.</p>","PeriodicalId":2,"journal":{"name":"ACS Applied Bio Materials","volume":null,"pages":null},"PeriodicalIF":4.6000,"publicationDate":"2023-11-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"ACS Applied Bio Materials","FirstCategoryId":"89","ListUrlMain":"https://doi.org/10.1007/s11004-023-10119-0","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"MATERIALS SCIENCE, BIOMATERIALS","Score":null,"Total":0}
引用次数: 0
Abstract
Probabilistic methods for geophysical inverse problems allow the use of arbitrarily complex prior information in principle. Geostatistical techniques, such as multiple-point statistics (MPS), for describing spatial correlation models and higher-order statistics have been proposed to achieve this inversion task, in which stochastic algorithms such as Markov chain Monte Carlo (McMC) are incorporated. However, stochastic sampling and optimization often require a large number of iterations, and thus geostatistical sampling of the prior model can become computationally demanding. To overcome this challenge, a deep learning model, namely conditional generative adversarial networks (CGANs), is proposed, which allows one to perform a random walk to sample the complex prior distribution. CGANs simulate conditional realizations conditioned to the available hard conditioning data, that is, direct measurements, while preserving the geometrical structure of the model parameters of interest and replicating the sequential Gibbs sampling algorithm. Despite the need for a training step, for a large number of simulations, CGANs are more efficient than traditional geostatistical simulation algorithms such as single normal equation simulation (SNESIM). The proposed methodology is used as part of the extended Metropolis algorithm to predict the distributions of categorical facies in two examples, a dune environment in the Gobi Desert and a channel system in an idealized subsurface reservoir, from indirect observational data such as acoustic impedance. The inversion results are compared to the extended Metropolis algorithm using standard MPS sampling.