Skip to main content
Log in

Three-Dimensional Latent Diffusion Models for Parameterizing and History Matching Facies Systems Under Hierarchical Uncertainty

  • Published:
Mathematical Geosciences Aims and scope Submit manuscript

Abstract

Geological parameterization procedures entail the mapping of a high-dimensional geomodel to a low-dimensional latent variable. These parameterizations can be very useful for history matching because the number of variables to be calibrated is greatly reduced, and the mapping can be constructed such that geological realism is automatically preserved. In this work, a parameterization method based on generative latent diffusion models (LDMs) is developed for three-dimensional channel-levee-mud systems. Geomodels with variable geological metaparameters (also referred to as hyperparameters or scenario parameters), specifically mud fraction, channel orientation, and channel width, are considered. A perceptual loss term is included during training to improve geological realism. For any set of scenario parameters, an (essentially) infinite number of realizations can be generated, so our LDM parameterizes over a very wide model space. New realizations constructed using the LDM procedure are shown to closely resemble reference geomodels, both visually and in terms of one- and two-point spatial statistics. Flow response distributions, for a specified set of injection and production wells, are also shown to be in close agreement between the two sets of models. The parameterization method is applied for ensemble-based history matching, with model updates performed in the LDM latent space, for cases involving uncertain scenario parameters. For three synthetic true models, we observe clear uncertainty reduction in both production forecasts and geological scenario parameters. The overall method is additionally shown to provide posterior geomodels consistent with the synthetic true model in each case.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Subscribe and save

Springer+
from €37.37 /Month
  • Starting from 10 chapters or articles per month
  • Access and download chapters and articles from more than 300k books and 2,500 journals
  • Cancel anytime
View plans

Buy Now

Price includes VAT (Norway)

Instant access to the full article PDF.

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8
Fig. 9
Fig. 10
Fig. 11
Fig. 12
Fig. 13
Fig. 14
Fig. 15
Fig. 16
Fig. 17

Similar content being viewed by others

References

Download references

Acknowledgements

We are grateful to the Stanford Doerr School of Sustainability and to the industrial affiliates of the Stanford Smart Fields Consortium for financial support. We also thank the SDSS Center for Computation for providing the computational resources used in this work.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Guido Di Federico.

Ethics declarations

Conflict of interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. The second author is on the editorial board of this journal.

Code availability

The source codes and datasets used in this work are available for download at https://github.com/guidodf09/ldm_3d_geomodel. The codes are based on the implementations provided in the https://github.com/huggingface/diffusers and https://github.com/Project-MONAI/GenerativeModels repositories. The Python libraries diffusers and monai are used in our implementation.

Appendices

Appendix A 3D-LDM Algorithms

Algorithm 1
figure a

VAE training

Algorithm 2
figure b

U-net training

Algorithm 3
figure c

Inference (generation)

Appendix B 3D-LDM Model Architecture

Tables 3 and 4 define the VAE and U-net architectures used in our 3D-LDM implementation.

Table 3 Variational autoencoder architecture used in this work, with input shape \((N_x, N_y, N_z, 1)\). Each ResBlock consists of GroupNorm \(\rightarrow \) SiLU \(\rightarrow \) Conv3D. An AttentionBlock is applied at the bottom level. Downsampling uses Conv3D; upsampling uses ConvTranspose3D
Table 4 Diffusion U-net architecture used in this work, with input shape \((n_x, n_y, n_z, 1)\), that is, the VAE latent space. The discrete step t embedding is applied inside ResBlocks as follows: Linear \(\rightarrow \) SiLU \(\rightarrow \) Linear \(\rightarrow \) 256-dimensional vector. Blocks have the same meaning as for the VAE

Rights and permissions

Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Di Federico, G., Durlofsky, L.J. Three-Dimensional Latent Diffusion Models for Parameterizing and History Matching Facies Systems Under Hierarchical Uncertainty. Math Geosci (2025). https://doi.org/10.1007/s11004-025-10245-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • Version of record:

  • DOI: https://doi.org/10.1007/s11004-025-10245-x

Keywords