Leveraging generative adversarial networks with spatially adaptive denormalization for multivariate stochastic seismic data inversion

TL;DR

This paper introduces SPADE-GANInv, an innovative iterative inversion method combining a spatially adaptive GAN with geostatistical simulation to accurately predict multiple coupled subsurface properties from seismic data.

Contribution

The paper presents SPADE-GANInv, a novel inversion framework that integrates a pre-trained SPADE-GAN with geostatistical simulation for multivariate seismic property prediction.

Findings

Accurately predicts facies, porosity, and impedance from seismic data.

Mitigates bias in prior data and incorporates well logs.

Effective on synthetic and field data scenarios.

Abstract

Probabilistic seismic inverse modeling often requires the prediction of both spatially correlated geological heterogeneities (e.g., facies) and continuous parameters (e.g., rock and elastic properties). Generative adversarial networks (GANs) provide an efficient training-image-based simulation framework capable of reproducing complex geological models with high accuracy and comparably low generative cost. However, their application in stochastic geophysical inversion for multivariate property prediction is limited, as representing multiple coupled properties requires large and unstable networks with high memory and training demands. A more recent variant of GANs with spatially adaptive denormalization (SPADE-GAN) enables the direct conditioning of facies spatial distributions on local probability maps. Leveraging on such features, an iterative geostatistical inversion algorithm is proposed, SPADE-GANInv, integrating a pre-trained SPADE-GAN with geostatistical simulation, for the prediction of facies and multiple correlated continuous properties from seismic data. The SPADE-GAN is trained to reproduce realistic facies geometries, while sequential stochastic co-simulation predicts the spatial variability of the facies-dependent continuous properties. At each iteration, a set of subsurface realizations is generated and used to compute synthetic seismic data. The realizations providing the highest similarity coefficient to the observed data are used to update the subsurface probability models in the next iteration. The method is demonstrated on both 2-D synthetic scenarios and field data, targeting the prediction of facies, porosity, and acoustic impedance from full-stack seismic data. Results show that the algorithm enables accurate multivariate prediction, mitigates the impact of biased prior data, and accommodates additional local conditioning such as well logs.