@conference {herrmann2025EAGEssf,
	title = {SAGE {\textemdash} Subsurface foundational model with AI-driven Geostatistical Extraction},
	booktitle = {EAGE Workshop},
	year = {2025},
	month = {6},
	abstract = {In this talk, we introduce a novel diffusion-based generative modeling approach for synthesizing high-fidelity subsurface velocity models. Once trained, this foundational model serves as a prior or a generative tool, enabling the sampling of velocity models that closely match the statistical distribution of geologically relevant Earth models. To achieve this, we leverage geostatistical information derived from curated datasets provided by the UK National Data Repository. Our training set comprises check-shot-corrected migrated images intersected by well logs. Numerical experiments on multiple synthetic datasets, including the Compass model and Synthoseis, demonstrate the model{\textquoteright}s ability to generate realistic subsurface velocity models. This capability offers a data-driven prior for full-waveform inversion and facilitates the generation of training samples for amortized neural networks to improve full-waveform inference.},
	keywords = {Bayesian inference, deep learning, EAGE, foundation models, full waveform inference, generative model, Imaging, Inverse problems, RTM, score, Summary Statistics, UK-NDR, Uncertainty quantification},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2025/herrmann2025EAGEssf/EAGE_WS22_Foundation_Herrmann.mp4},
	author = {Huseyin Tuna Erdinc and Thales Souza and Felix J. Herrmann}
}
@conference {erdinc2025ML4SEISMICspm,
	title = {SAGE - Subsurface modeling with AI-driven Geostatistical Extraction and evaluation on North Sea Data},
	booktitle = {ML4SEISMIC Partners Meeting},
	year = {2025},
	month = {11},
	abstract = {Traditional machine learning based approaches often depend on large, high-quality and complete datasets of 2D Earth propety models, which can be difficult to obtain in subsurface applications. In this study, we leverage score-based generative models to synthesize high-fidelity geophysical properties (velocity, density, etc.) trained on a combination of sparse but tied well logs and seismic lines without relying on densely sampled velocity models and adress the sampling-density differences by enabling coherent integration of these complementary sources of information. With experiments on multiple synthetic datasets from subsurface models featuring diverse geological structures (e.g., faults, salt bodies), we demonstrate that our approach can accurately capture long-range geological correlations and align well with unseen ground-truth properties. Finally, we highlight the potential of our method to generalize to field data by conducting fine-tuning on a curated dataset derived from real well-log and seismic data currated from the UK Data Repository for the North Sea. This study represents an important step toward the development of foundational models for inverting physical properties and provides inputs for applications such as full-waveform inference (WISE/ASPIRE), supervised learning algorithms, and enhanced seismic-based subsurface modeling.},
	keywords = {Bayesian inference, deep learning, generative model, Imaging, Inverse problems, ML4SEISMIC, North-sea data, power-scaling, RTM, score, SLIM, Summary Statistics, Uncertainty quantification, WISE},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ML4SEISMIC/2025/erdinc2025ML4SEISMICspm},
	author = {Huseyin Tuna Erdinc and Ipsita Bhar and Thales Souza and Rafael Orozco and Felix J. Herrmann}
}
@conference {bhar2025ML4SEISMICsdc,
	title = {Seismic Dataset Curation from UK National Data Repository to Validate SAGE and WISE},
	booktitle = {ML4SEISMIC Partners Meeting},
	year = {2025},
	month = {11},
	abstract = {This work presents a curated seismic dataset pipeline developed for data residing in the UK National Data Repository (UK NDR) in support of advanced generative geophysical modeling and machine learning. The workflow integrates checkshots and 3D post-stack seismic data volumes to generate accurate depth-domain seismic sections. Custom interpolation algorithms are used to construct velocity models from checkshot data, which are then used for time-to-depth conversion with [OpendTect]https://dgbes.com/software/opendtect. The converted seismic cubes are subsequently processed in Python using [SegySAK]https://segysak.readthedocs.io/en/latest/, a python package, enabling efficient extraction and visualization of 2D lines. These datasets will be used to further validate the workflow of SAGE (Subsurface foundational model for AI-driven Geostatistical Extraction), which focuses on large-scale representation learning for the geosciences. The trained generative models will be used as input to WISE and ASPIRE that form foundation models for seismic inference designed to characterize the subsurface. This curated dataset pipeline streamlines seismic data preparation, enhances reproducibility, and bridges the gap between conventional geophysical workflows and emerging data-driven inference methods.},
	keywords = {Bayesian inference, generative model, Imaging, ML4SEISMIC, North-sea data, RTM, SLIM, WISE},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ML4SEISMIC/2025/bhar2025ML4SEISMICsdc},
	author = {Ipsita Bhar and Huseyin Tuna Erdinc and Thales Souza and Rafael Orozco and Felix J. Herrmann}
}
@presentation {bruer2025GSCSsmc,
	title = {Seismic monitoring of CO2 using ensemble Kalman filtering},
	year = {2025},
	note = {(GSCS, Gatech)},
	month = {02},
	keywords = {CCS, EnKF, fluid flow, GCS, GSCS, Imaging, monitoring, RTM, SLIM, Uncertainty quantification},
	url = {https://slim.gatech.edu/Publications/Public/Posters/GSCS/2025/bruer2025GSCSsmc/2025GSCS_poster.pdf},
	author = {Grant Bruer and Abhinav Prakash Gahlot and Felix Herrmann and Edmond Chow}
}
@conference {gahlot2025IMAGEsrp,
	title = {Sensitivity-aware rock physics enhanced digital shadow for underground-energy storage monitoring},
	booktitle = {International Meeting for Applied Geoscience and Energy},
	year = {2025},
	month = {8},
	keywords = {Amortized Variational Inference, Bayesian inference, conditional normalizing flows, data assimilation, deep learning, digital twin, experimental design, GCS, IMAGE, Imaging, Inverse problems, permeability, reservoir simulation, rock physics, RTM, SEG, Summary Statistics, Uncertainty quantification, WISE},
	doi = {10.48550/arXiv.2504.14405},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2025/gahlot2025IMAGEsrp},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2025/gahlot2025IMAGEsrp/IMAGE_Abhinav.pdf},
	author = {Abhinav Prakash Gahlot and Huseyin Tuna Erdinc and Felix J. Herrmann}
}
@presentation {erdinc2024ML4SEISMICsfm,
	title = {SAGE {\textendash} Subsurface foundational model with AI-driven Geostatical Extraction},
	year = {2024},
	month = {11},
	abstract = {In this study, we present a novel approach for synthesizing diverse subsurface velocity models using diffusion-based generative models. Traditional methods often depend on large, high-quality datasets of 2D velocity models, which can be difficult to obtain in subsurface applications. In contrast, our method leverages incomplete well and seismic data to generate high-fidelity velocity samples without requiring fully sampled training datasets.The results demonstrate that the generative model accurately captures long-range geological structures and aligns well with unseen ground-truth velocity models. Furthermore, it is shown that the diversity of generated velocity models can be increased through prior guidance in the training phase, and model uncertainties can be reduced with well conditioning during inference.Experiments conducted with multiple datasets (BG, Synthoseis, and North Sea data) and velocity models featuring various geological structures (e.g., faults, salt bodies) suggest that our approach facilitates realistic subsurface velocity synthesis, providing valuable inputs for full-waveform inversion and enhancing seismic-based subsurface modeling.},
	keywords = {Amortized Variational Inference, Bayesian inference, deep learning, diffusion models, FWI, Imaging, Inverse problems, ML4SEISMIC, RTM, SLIM, Summary Statistics, Uncertainty quantification, WISE},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ML4SEISMIC/2024/erdinc2024ML4SEISMICsfm},
	author = {Huseyin Tuna Erdinc and Rafael Orozco and Felix J. Herrmann}
}
@presentation {bruer2024ML4SEISMICsmp,
	title = {Seismic monitoring of CO2 plume dynamics using ensemble Kalman filtering},
	year = {2024},
	month = {11},
	abstract = {Monitoring CO 2  injected and stored in subsurface reservoirs is critical for avoiding failure scenarios and enables real-time optimization of CO 2  injection rates. Sequential Bayesian data assimilation (DA) is a statistical method for combining information over time from multiple sources to estimate a hidden state, such as the spread of the subsurface CO 2  plume. An example of scalable and efficient sequential Bayesian DA is the ensemble Kalman filter (EnKF). We improve upon existing DA literature in the seismic-CO 2  monitoring domain by applying this scalable DA algorithm to a high-dimensional CO 2  reservoir using two-phase flow dynamics and time-lapse full waveform seismic data with a realistic surface-seismic survey design. We show more accurate estimates of the CO 2  saturation field using the EnKF compared to using either the seismic data or the fluid physics alone. Furthermore, we test a range of values for the EnKF hyperparameters and give guidance on their selection for seismic CO 2  reservoir monitoring.},
	keywords = {Bayesian inference, data assimilation, ensemble Kalman filter, FWI, Imaging, Inverse problems, ML4SEISMIC, RTM, SLIM, two-phase flow, Uncertainty quantification},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ML4SEISMIC/2024/bruer2024ML4SEISMICsmp},
	author = {Grant Bruer and Abhinav Prakash Gahlot and Edmond Chow and Felix J. Herrmann}
}
@presentation {bhar2024ML4SEISMICshs,
	title = {Sensitivity of Horizontal shear waves on Geological Carbon Storage},
	year = {2024},
	month = {11},
	abstract = {With the growing focus on Carbon Capture and Storage (CCS) activities, the significance of shear wave velocity (Vs) in monitoring CO2 is becoming increasingly important. Shear waves are instrumental in detecting CO2 leakage, evaluating CCS-induced seismicity, and identifying caprock failure. According to Biot-Gassmann rock physics models, changes in seismic velocity are expected during CO2 injection. This work specifically focuses on the horizontal component of shear wave velocity (Vsh). Shear (Sh) waves are more sensitive to density variations than P-wave or Sv wave velocities, leading us to anticipate changes in Sh velocity as CO2 saturation increases. While P-wave velocity tends to decrease with supercritical CO2 injection, the investigation of its effects on Sh waves presents a particularly compelling area of study. Therefore, we simulated the wave equation incorporating both P-wave and Sh wave velocities to better understand their impacts on imaging and CO2 injection processes.},
	keywords = {FWI, Imaging, Inverse problems, ML4SEISMIC, RTM, SH wave, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ML4SEISMIC/2024/bhar2024ML4SEISMICshs},
	author = {Ipsita Bhar and Abhinav Prakash Gahlot and Felix J. Herrmann}
}
@presentation {yin2022ML4SEISMICsfg,
	title = {Simulation-based framework for geological carbon storage monitoring},
	year = {2022},
	month = {11},
	abstract = {While various monitoring modalities exist to track the behavior of CO2 plumes to ensure safe operations and compliance with regulatory requirements, active 3D time-lapse seismic monitoring has proven superior but costly. At SLIM, we aim to reduce the operating costs by optimizing acquisition design, to help drive innovations in seismic monitoring acquisition design and imaging, and to test novel time-lapse acquisition and imaging technologies in silico at scale. In this talk, we will introduce our open-source software platform simulation-based monitoring design framework. We demonstrate how to make use of proxy models for seismic properties derived from real 3D imaged seismic and well data to conduct realistic synthetic geological carbon storage projects. Furthermore, we discuss our proposed sparse non-replicated seismic acquisition and cutting-edge methodology to recover the dense data or to directly image the sparse non-replicated via joint recovery model. This automatic workflow ends with deep neural classifiers to detect potential CO2 leakage over seal through pressure-induced fault openings. We envisage the development of an automatic workflow to handle the large number of continuously monitored CO2 injection sites needed to help combat climate change.},
	keywords = {CAM, CCS, classification, deep learning, explainability, GCS, Imaging, JRM, machine learning, ML4SEISMIC, SLIM, time-lapse},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ML4SEISMIC/2022/yin2022ML4SEISMICsfg/index.html},
	author = {Ziyi Yin and Huseyin Tuna Erdinc and Abhinav Prakash Gahlot and Mathias Louboutin and Felix J. Herrmann}
}
@conference {siahkoohi2020SIAMISsiu,
	title = {Seismic Imaging with Uncertainty Quantification: Sampling from the Posterior with Generative Networks},
	booktitle = {SIAM Conference on Imaging Science},
	year = {2020},
	note = {IS20},
	month = {7},
	abstract = {n/a},
	keywords = {Imaging, invertible networks, IS20, McMC, Uncertainty quantification},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SIAMIS/2020/siahkoohi2020SIAMISsiu/siahkoohi2020SIAMISsiu_pres.pdf},
	author = {Ali Siahkoohi and Philipp A. Witte and Mathias Louboutin and Felix J. Herrmann and Gabrio Rizzuti}
}
@mastersthesis {witte2020THsal,
	title = {Software and algorithms for large-scale seismic inverse problems},
	year = {2020},
	note = {(PhD)},
	month = {02},
	school = {Georgia Institute of Technology},
	type = {phd},
	address = {Atlanta},
	abstract = {Seismic imaging and parameter estimation are an import class of inverse problems with practical relevance in resource exploration, carbon control and monitoring systems for geohazards. The goal of seismic inverse problems is to image subsurface geological structures and estimate physical rock properties such as wave speed or density. Mathematically, this can be achieved by solving an optimization problem in which we minimize the mismatch between numerically modeled data and observed data from a seismic survey. As wave propagation through a medium is described by wave equations, seismic inverse problems involve solving a large number of partial differential equations (PDEs) during numerical optimization using finite difference modeling, making them computationally expensive. Additionally, seismic inverse problems are typically ill-posed, non-convex or ill-conditioned, thus making them challenging from a mathematical standpoint as well. Similar to the field of deep learning, this calls for software that is not only optimized for performance, but also enables geophysical domain specialists to experiment with algorithms in high-level programming languages and using different computing environments, such as high-performance computing (HPC) clusters or the cloud. Furthermore, they call for the adaption of dimensionality reduction techniques and stochastic algorithms to address computational cost from the algorithmic side. This thesis makes three distinct contributions to address computational challenges encountered in seismic inverse problems and to facilitate algorithmic development in this field. Part one introduces a large-scale framework for seismic modeling and inversion based on the paradigm of separation of concerns, which combines a user interface based on domain specific abstractions with a Python package for automatic code generation to solve the underlying PDEs. The modular code structure makes it possible to manage the complexity of a seismic inversion code, while matrix-free linear operators and data containers enable the implementation of algorithms in a fashion that closely resembles the underlying mathematical notation. The second contribution of this thesis is an algorithm for seismic imaging, that addresses its high computational cost and large memory imprint through a combination of on-the-fly Fourier transforms, stochastic sampling techniques and sparsity-promoting optimization. The algorithm combines the best of both time- and frequency-domain inversion, as the memory imprint is independent of the number of modeled time steps, while time-to-frequency conversions avoid the need to solve Helmholtz equations, which involve inverting ill-conditioned matrices. Part three of this thesis introduces a novel approach for adapting the cloud for high-performance computing applications like seismic imaging, which does not rely on a fixed cluster of permanently running virtual machines. Instead, computational resources are automatically started and terminated by the cloud environment during runtime and the workflow takes advantage of cloud-native technologies such as event-driven computations and containerized batch processing. The performance and cost analysis shows that this approach is able to address current shortcomings of the cloud such as inferior resilience, while at the same time reducing operating cost up to an order of magnitude. As such, the workflow provides a strategy for cost effectively running large-scale seismic imaging problems in the cloud and is a viable alternative to conventional HPC clusters.},
	keywords = {algorithms, AWS, FWI, Imaging, PhD, RTM, seismic, software},
	url = {https://slim.gatech.edu/Publications/Public/Thesis/2020/witte2020THsal/witte2020THsal.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Thesis/2020/witte2020THsal/witte2020THsal_pres.pdf},
	author = {Philipp A. Witte}
}
@conference {aravkin2013ICASSPssi,
	title = {Sparse seismic imaging using variable projection},
	booktitle = {ICASSP},
	year = {2013},
	month = {05},
	abstract = {We consider an important class of signal processing problems where the signal of interest is known to be sparse, and can be recovered from data given auxiliary information about how this data was generated. For example, a sparse green{\textquoteright}s function may be recovered from seismic experimental data using sparsity optimization when the source signature is known. Unfortunately, in practice this information is often missing, and must be recovered from data along with the signal using deconvolution techniques. In this paper, we present a novel methodology to simulta- neously solve for the sparse signal and auxiliary parameters using a recently proposed variable projection technique. Our main contribution is to combine variable projection with spar- sity promoting optimization, obtaining an efficient algorithm for large-scale sparse deconvolution problems. We demon- strate the algorithm on a seismic imaging example.},
	keywords = {Imaging, Optimization, sparsity, variable projection},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ICASSP/2013/aravkin2013ICASSPssi/aravkin2013ICASSPssi.pdf},
	author = {Aleksandr Y. Aravkin and Tristan van Leeuwen and Ning Tu}
}
@conference {li2012SEGspmamp,
	title = {Sparsity-promoting migration accelerated by message passing},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {31},
	year = {2012},
	month = {11},
	pages = {1-6},
	publisher = {SEG},
	organization = {SEG},
	abstract = {Seismic imaging via linearized inversion requires multiple iterations to minimize the least-squares misfit as a function of the medium perturbation. Unfortunately, the cost for these iterations are prohibitive because each iteration requires many wave-equation simulations, which without direct solvers require an expensive separate solve for each source. To overcome this problem, we use dimensionality-reduction to decrease the size of seismic imaging problem by turning the large number of sequential shots into a much small number of simultaneous shots. In our approach, we take advantage of sparsifying transforms to remove source crosstalk resulting from randomly weighting and stacking sequential shots into a few super shots. We also take advantage of the fact that the convergence of large-scale sparsity-promoting solvers can be improved significantly by borrowing ideas from message passing, which are designed to break correlation built up between the linear system and the model iterate. In this way, we arrive at a formulation where we run the sparsity-promoting solver for a relatively large number of very iterations. Aside from leading to a significant speed up, our approach had the advantage of greatly reducing the memory imprint and IO requirements. We demonstrate this feature by solving a sparsity-promoting imaging problem with operators of reverse-time migration, which is computationally infeasible without the dimensionality reduction.},
	keywords = {Imaging, inversion, SEG},
	doi = {10.1190/segam2012-1500.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2012/li2012SEGspmamp/li2012SEGspmamp.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2012/li2012SEGspmamp/li2012SEGspmamp_pres.pdf},
	author = {Xiang Li and Felix J. Herrmann}
}
@conference {tu2011EAGEspmsrm,
	title = {Sparsity-promoting migration with surface-related multiples},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2011},
	month = {05},
	abstract = {Multiples, especially the surface-related multiples, form a significant part of the total up-going wave- field. If not properly dealt with, they can lead to false reflectors in the final image. So conventionally practitioners remove them prior to migration. Recently research has revealed that multiples can actually provide extra illumination so different methods are proposed to address the issue that how to use multiples in seismic imaging, but with various kinds of limitations. In this abstract, we combine primary estimation and sparsity-promoting migration into one convex-optimization process to include information from multiples. Synthetic examples show that multiples do make active contributions to seismic migration. Also by this combination, we can benefit from better recoveries of the Greens function by using sparsity-promoting algorithms since reflectivity is sparser than the Greens function.},
	keywords = {EAGE, Imaging, Presentation, Processing},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2011/tu11EAGEspmsrm/tu11EAGEspmsrm.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2011/tu11EAGEspmsrm/tu11EAGEspmsrm_pres.pdf},
	url2 = {http://earthdoc.eage.org/publication/publicationdetails/?publication=50369},
	author = {Ning Tu and Tim T.Y. Lin and Felix J. Herrmann}
}
@article {herrmann2008ACHAsac,
	title = {Sparsity- and continuity-promoting seismic image recovery with curvelet frames},
	journal = {Applied and Computational Harmonic Analysis},
	volume = {24},
	number = {2},
	year = {2008},
	month = {03},
	pages = {150-173},
	abstract = {A nonlinear singularity-preserving solution to seismic image recovery with sparseness and continuity constraints is proposed. We observe that curvelets, as a directional frame expansion, lead to sparsity of seismic images and exhibit invariance under the normal operator of the linearized imaging problem. Based on this observation we derive a method for stable recovery of the migration amplitudes from noisy data. The method corrects the amplitudes during a post-processing step after migration, such that the main additional cost is one ap- plication of the normal operator, i.e. a modeling followed by a migration. Asymptotically this normal operator corresponds to a pseudodifferential operator, for which a convenient diagonal approximation in the curvelet domain is derived, including a bound for its error and a method for the estimation of the diagonal from a compound operator consisting of discrete implementations for the scattering operator and its adjoint the migration operator. The solution is formulated as a nonlinear optimization problem where sparsity in the curvelet domain as well as continuity along the imaged reflectors are jointly promoted. To enhance sparsity, the $ell_1$-norm on the curvelet coefficients is minimized, while continuity is promoted by minimizing an anisotropic diffusion norm on the image. The performance of the recovery scheme is evaluated with a time-reversed {\textquoteright}wave-equation{\textquoteright} migration code on synthetic datasets, including the complex SEG/EAGE AA salt model.},
	keywords = {curvelet transform, Imaging, Processing, SLIM},
	doi = {10.1016/j.acha.2007.06.007},
	url = {https://slim.gatech.edu/Publications/Public/Journals/ACHA/2008/herrmann2008ACHAsac/herrmann2008ACHAsac.pdf},
	author = {Felix J. Herrmann and Peyman P. Moghaddam and Christiaan C. Stolk}
}
