@conference {herrmann2025SEGSPEcua,
	title = {Context- and uncertainty-aware Digital Twins for Integrated Reservoir Management},
	booktitle = {SEG/SPE 2nd Integrated Reservoir Characterization, Surveillance, and Management Workshop: Advances and Challenges},
	year = {2025},
	month = {7},
	abstract = {We present a new framework for Digital Twins (DTs) that supports uncertainty-aware monitoring, optimization, and risk management of underground storage operations, with a particular focus on Geological Carbon Storage (GCS) and Enhanced Oil Recovery (EOR). The proposed DT addresses key challenges in reservoir characterization, surveillance, and management by integrating generative AI methods with simulation-based inference, enabling scalable, probabilistic decision support under uncertainty. Our approach combines sensitivity-aware amortized Bayesian inference (SA-ABI) with ensemble-based data assimilation to efficiently incorporate multimodal time-lapse observations, such as seismic and well data, into dynamic reservoir models. By embedding context-awareness into neural networks, we account for variability in rock-physics relationships linking pressure, saturation, and seismic responses, and quantify uncertainties arising from model misspecifications{\textendash}-e.g. usage of patchy versus uniform saturation models. Computational efficiency is maintained through weight sharing across the network trained on different rock-physics scenarios, enabling rapid evaluation of the impact changing the underlying rock-physics model. The DT is built on industry-strength open-source software: {\textquoteleft}JutulDarcy.jl{\textquoteleft} for fully coupled multiphase flow simulations and {\textquoteleft}JUDI.jl{\textquoteleft}, a seismic inversion package leveraging the {\textquoteleft}Devito{\textquoteleft} finite-difference framework. This computational foundation allows scalable integration of physics-based simulations and data assimilation across large reservoir models, while retaining flexibility to adapt to varying physical contexts such as the inclusion of pressure-dependence of the compressional wavespeed. The resulting DT not only delivers posterior distributions for the reservoir{\textquoteright}s state (pressure/saturation) conditioned on available observations but also enables rapid {\textquotedblleft}what-if{\textquotedblright} scenario analysis by drawing posterior samples across a range of physical contexts without retraining. This capability supports proactive storage management, optimization of injection operations, and improved assurance of storage conformance.},
	keywords = {Amortized Variational Inference, Bayesian inference, conditional normalizing flows, data assimilation, deep learning, digital twin, experimental design, GCS, Imaging, Inverse problems, permeability, reservoir simulation, rock physics, RTM, SEG, SPE, Summary Statistics, Uncertainty quantification, WISE},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEGSPE/2025/herrmann2025SEGSPEcua},
	author = {Abhinav Prakash Gahlot and Ipsita Bhar and Haoyun Li and Felix J. Herrmann}
}
@conference {herrmann2025AIDT4EScad,
	title = {Context-aware Digital Twin for Underground Storage},
	booktitle = {Artificial Intelligence and Digital Twins for Earth Systems 2025, Austin, TX},
	year = {2025},
	month = {9},
	keywords = {AIDT4ES, Amortized Variational Inference, Bayesian inference, conditional normalizing flows, data assimilation, deep learning, digital twin, experimental design, GCS, Imaging, Inverse problems, permeability, reservoir simulation, rock physics, RTM, Summary Statistics, Uncertainty quantification},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/AIDT4ES/2025/herrmann2025AIDT4EScad},
	author = {Abhinav Prakash Gahlot and Haoyun Li and Ipsita Bhar and Huseyin Tuna Erdinc and Rafael Orozco and Ziyi Yin and Felix J. Herrmann}
}
@conference {herrmann2025IMSIcad,
	title = {Context-Aware Digital Twin for Underground Storage Operations and Decision Making},
	booktitle = {Data Assimilation and Inverse Problems for Digital Twins Workshop 205, IMSI, Chicago},
	year = {2025},
	month = {10},
	keywords = {Amortized Variational Inference, Bayesian inference, conditional normalizing flows, data assimilation, decision making, digital twin, experimental design, GCS, Imaging, Inverse problems, permeability, reservoir simulation, rock physics, RTM, Summary Statistics, Uncertainty quantification},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/IMSI/2025/herrmann2025IMSIcad},
	author = {Abhinav Prakash Gahlot and Haoyun Li and Ipsita Bhar and Huseyin Tuna Erdinc and Rafael Orozco and Ziyi Yin and Felix J. Herrmann}
}
@conference {park2025ML4SEISMICftn,
	title = {Fisher-Informed Training of Neural Operators for Reliable PDE Inversion},
	booktitle = {ML4SEISMIC Partners Meeting},
	year = {2025},
	month = {11},
	abstract = {Neural operators have shown strong performance for PDE-solution learning; however, their effectiveness in PDE-constrained optimization, such as inversion, remains limited, in part because inaccurately learned derivatives can mislead gradient-based updates. Standard training typically optimizes only the forward least-squares fit, so during iterative inversion the optimization variable can drift outside the training distribution, further degrading quality of derivative information and inversion performance. To address this limitation, we propose a training algorithm designed for inversion: the surrogate is trained not only to predict forward outputs but also to learn observation-relevant gradients, which leads to Fisher-information directed parameter updates. We evaluate forward and gradient accuracy, as well as the inversion trajectory, on problems governed by PDEs{\textemdash}Darcy flow and laminar (incompressible) Navier{\textendash}Stokes{\textemdash}which are relevant to porous-media fluid flow and energy-storage applications. Our results suggest that training for gradient fidelity with forward accuracy unlocks a pathway to reliable, efficient, neural operator{\textendash}based inversion.},
	keywords = {deep learning, derivatve-informed, Fisher-information, FNO, GCS, Inverse problems, ML4SEISMIC, permeability, reduced-order modeling, reservoir simulation, SLIM, Uncertainty quantification},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ML4SEISMIC/2025/park2025ML4SEISMICftn},
	author = {Jeongjin Park and Grant Bruer and Huseyin Tuna Erdinc and Richard Rex Arockiasamy and Nisha Chandramoorthy and Felix J. Herrmann}
}
@conference {gahlot2025SIAMiac,
	title = {It{\textquoteright}s All in the Context: Sensitivity-Aware Digital Twins for Underground Storage},
	booktitle = {SIAM Conference on Mathematical \& Computational Issues in the Geosciences, Baton Rouge},
	year = {2025},
	month = {10},
	keywords = {Amortized Variational Inference, Bayesian inference, conditional normalizing flows, data assimilation, decision making, digital twin, experimental design, GCS, Imaging, Inverse problems, permeability, reservoir simulation, rock physics, RTM, SIAM, Summary Statistics, Uncertainty quantification},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SIAM/2025/gahlot2025SIAMiac},
	author = {Abhinav Prakash Gahlot and Ipsita Bhar and Felix J. Herrmann}
}
@conference {gahlot2025ML4SEISMICmrm,
	title = {Mitigating Rock-Physics Model Misspecification in Digital Shadows via Amortized Bayesian Inference},
	booktitle = {ML4SEISMIC Partners Meeting},
	year = {2025},
	month = {11},
	abstract = {Digital shadows for subsurface monitoring rely on simulation-based inference to map seismic observations to evolving reservoir states. However, when the simulations are generated using simplified or different rock physics models than those governing the observed data, model misspecification leads to biased posterior estimates and degraded uncertainty quantification. We address this challenge through Bayesian frameworks including rock-physics marginalization and a context-aware amortized inference using conditional normalizing flows. The framework learns a mapping from seismic responses to pressure and saturation while accounting for variability across different randomly drawn models for the rock-physics. To enhance sensitivity to fluid and stress effects, we also incorporate pressure-dependent rock physics into the seismic forward model and perform inference in the Radon domain of Common Image Gathers (CIGs) instead of conventional RTM images. Because fluids exhibit distinct angle-dependent signatures in AVA responses, Radon-transformed CIGs are expected to provide improved separation between pressure and saturation effects, leading to more accurate and physically consistent state estimates. Thus, it enables an uncertainty-aware, physics-informed digital shadows under rock-physics model mismatch.},
	keywords = {Amortized Variational Inference, Bayesian inference, CIG, conditional normalizing flows, data assimilation, deep learning, digital twin, experimental design, GCS, Imaging, Inverse problems, ML4SEISMIC, permeability, reservoir simulation, rock physics, RTM, SLIM, Summary Statistics, Uncertainty quantification, WISE},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ML4SEISMIC/2025/gahlot2025ML4SEISMICmrm},
	author = {Abhinav Prakash Gahlot and Ipsita Bhar and Felix J. Herrmann}
}
@conference {park2025IMAGErdn,
	title = {A reduced-order derivative-informed neural operator for subsurface fluid-flow},
	booktitle = {International Meeting for Applied Geoscience and Energy},
	year = {2025},
	month = {8},
	keywords = {deep learning, derivatve-informed, Fisher-information, FNO, GCS, IMAGE, Inverse problems, permeability, reduced-order modeling, reservoir simulation, SEG, Uncertainty quantification},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2025/park2025IMAGErdn/IMAGE_Jayjay.pdf},
	author = {Jeongjin Park and Grant Bruer and Huseyin Tuna Erdinc and Abhinav Prakash Gahlot and Felix J. Herrmann}
}
@conference {herrmann2025SIAMrdt,
	title = {Robust Digital Twin for Geological Carbon Storage},
	booktitle = {SIAM Conference on Computational Science and Engineering (CSE25)},
	year = {2025},
	note = {(SIAM CSE25, Dallas)},
	month = {03},
	abstract = {While recent work has shown that conditional neural networks can be used to effectuate prior-to-posterior mappings during ensemble Kalman filtering, their accuracy relies on the quality of the training ensemble, consisting of the simulated time-advanced nonlinear state and possibly nonlinear observations. Its precision also depends on the generative neural networks{\textquoteright} ability to generalize. In this work, we will report progress on how to make our Digital Twin more robust with respect to distribution shifts and modeling errors, e.g. due to erroneous choices for the physics within the Forecast Step. To this end, we will follow a two-pronged approach consisting of augmenting the Forecast step, e.g. by including different models for the rock physics that map the Digital Twin{\textquoteright}s state to seismic properties, and by adding more physics during the prior-to-posterior mapping. In its current implementation, training of our Digital Twin{\textquoteright}s neural networks relies on relative low-fidelity amortization of the Forecast Ensemble, which is susceptible to distribution shifts. By conducting non-amortized inference informed by the physics and low-fidelity prior-to-posterior mapping, we propose an approach that is more robust than the current implementation of the Digital Twin.},
	keywords = {Bayesian inference, conditional normalizing flows, CSE25, digital twin, rock physics, sequential Bayes, SIAM, Uncertainty quantification},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SIAMCSE/2025/herrmann2025SIAMrdt},
	author = {Felix J. Herrmann and Abhinav Prakash Gahlot and Rafael Orozco and Ziyi Yin and Grant Bruer}
}
@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}
}
@article {bruer2024smpd,
	title = {Seismic monitoring of CO2 plume dynamics using ensemble Kalman filtering},
	journal = {Transactions on Geoscience and Remote Sensing},
	year = {2025},
	note = {(Accepted in Transactions on Geoscience and Remote Sensing)},
	month = {5},
	abstract = {Monitoring carbon dioxide (CO2) injected and stored in subsurface reservoirs is critical for avoiding failure scenarios and enables real-time optimization of CO2 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 CO2 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-CO2 monitoring domain by applying this scalable DA algorithm to a high-dimensional CO2 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 CO2 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 CO2 reservoir monitoring.},
	keywords = {Ensebmle Kalman Filter, GCS, Inverse problems, non-linear dynamical systems, seismic},
	doi = {10.48550/arXiv.2409.05193},
	author = {Grant Bruer and Abhinav Prakash Gahlot and Edmond Chow 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}
}
@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}
}
@conference {herrmann2024SIAMUQdt,
	title = {An Uncertainty-Aware Digital Twin for Geological Carbon Storage},
	booktitle = {SIAM Conference on Uncertainty Quantification},
	year = {2024},
	note = {(SIAM UQ, Trieste)},
	month = {03},
	abstract = {Arguably, Geological Carbon Storage constitutes the only truly scalable net-negative carbon emission technology. To mitigate its risks and optimize its operations, an uncertainty-aware Digital Twin is being developed. To leverage existing fluid-flow and seismic simulation and imaging capabilities, the envisioned twin combines techniques from sequential and simulation-based Bayesian inference to train its deep generative neural networks to draw samples from the posterior of the Digital Twin{\textquoteright}s state. Because these samples are conditioned on observed time-lapse field data, these twins are capable of capturing the dynamics of CO2 plumes and their uncertainty.},
	keywords = {Bayesian inference, CCS, deep learning, digital twin, GCS, Imaging, monitoring, SIAM, Uncertainty quantification},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SIAMUQ/2024/digital-twin},
	author = {Abhinav Prakash Gahlot and Rafael Orozco and Haoyun Li and Grant Bruer and Ziyi Yin and Mathias Louboutin and Felix J. Herrmann}
}
@presentation {bruer2023ML4SEISMICcrm,
	title = {CO2 reservoir monitoring through Bayesian data assimilation},
	year = {2023},
	month = {11},
	abstract = {Carbon capture and storage can be implemented by injecting supercritical carbon dioxide (CO2) into geological carbon reservoirs for long-term containment. Monitoring the pressure and saturation of the CO2 is necessary to optimize the injection amount without causing CO2 leakage or seismic activity. Directly measuring the CO2 at locations within the reservoir requires expensive drilling procedures that may damage the reservoir, so direct measurements are sparse and usually lie along injection and production wells. Indirect measurements such as seismic data are typically noisy, and inverting for the CO2 state is ill-posed. Bayesian data assimilation techniques allow us to integrate known physics for CO2 flow into this inversion process. The most well-established data assimilation algorithms are the family of Kalman filters. The ensemble Kalman filter is designed to efficiently work with large problem sizes and nonlinearity. In this work, we apply the ensemble Kalman filter to seismic measurements of a CO2 reservoir, yielding an estimate of the saturation and pressure fields with quantified uncertainties. This method models the CO2 plume state as a random field with a known distribution and assimilates information from seismic measurements with information from a physics model describing the CO2 flow. We show that the data assimilation strategy is a valuable contribution to advancing reservoir monitoring technology.},
	keywords = {bayes, CCS, data assimilation, kalman, ML4SEISMIC, monitoring, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ML4SEISMIC/2023/bruer2023ML4SEISMICcrm},
	author = {Grant Bruer and Felix J. Herrmann and Edmond Chow}
}
@conference {bisbas2020IPDPStbf,
	title = {Temporal blocking of finite-difference stencil operators with sparse "off-the-grid" sources},
	booktitle = {IEEE International Parallel and Distributed Processing Symposium},
	year = {2021},
	note = {(IPDPS, virtual)},
	month = {5},
	abstract = {Stencil kernels dominate a range of scientific applications including seismic and medical imaging, image processing, and neural networks. Temporal blocking is a performance optimisation that aims to reduce the required memory bandwidth of stencil computations by re-using data from the cache for multiple time steps. It has already been shown to be beneficial for this class of algorithms. However, optimising stencils for practical applications remains challenging. These computations often include sparsely located operators, not aligned with the computational grid ("off-the-grid"). For example, our work is motivated by sources that inject a wavefield and measurements interpolating grid values. The resulting data dependencies make the adoption of temporal blocking much more challenging. We propose a methodology to inspect these data dependencies and reorder the computation, leading to performance gains in stencil codes where temporal blocking has not been applicable. We implement this novel scheme in the Devito domain-specific compiler toolchain. Devito implements a domain-specific language embedded in Python to generate optimised partial differential equation solvers using the finite-difference method from high-level symbolic problem definitions. We evaluate our scheme using isotropic acoustic, anisotropic acoustic and isotropic elastic wave propagators of industrial significance. Performance evaluation, after auto-tuning, shows that this enables substantial performance improvement through temporal blocking, over highly-optimised vectorized spatially-blocked code of up to 1.6x.},
	keywords = {finite-difference, HPC, performance, wave-equation},
	doi = {10.1109/IPDPS49936.2021.00058},
	url = {https://arxiv.org/pdf/2010.10248.pdf},
	author = {George Bisbas and Fabio Luporini and Mathias Louboutin and Rhodri Nelson and Gerard Gorman and Paul H. J. Kelly}
}
@techreport {louboutin2020SCsta,
	title = {Scaling through abstractions {\textendash} high-performance vectorial wave simulations for seismic inversion with Devito},
	number = {TR-CSE-2020-3},
	year = {2020},
	month = {4},
	institution = {Georgia Institute of Technology},
	abstract = {[Devito] is an open-source Python project based on domain-specific language and compiler technology. Driven by the requirements of rapid HPC applications development in exploration seismology, the language and compiler have evolved significantly since inception. Sophisticated boundary conditions, tensor contractions, sparse operations and features such as staggered grids and sub-domains are all supported; operators of essentially arbitrary complexity can be generated. To accommodate this flexibility whilst ensuring performance, data dependency analysis is utilized to schedule loops and detect computational-properties such as parallelism. In this article, the generation and simulation of MPI-parallel propagators (along with their adjoints) for the pseudo-acoustic wave-equation in tilted transverse isotropic media and the elastic wave-equation are presented. Simulations are carried out on industry scale synthetic models in a HPC Cloud system and reach a performance of 28TFLOP/s, hence demonstrating Devito{\textquoteright}s suitability for production-grade seismic inversion problems.},
	keywords = {devito, elastic, finite-difference, HPC, large-scale, RTM, TTI},
	url = {https://slim.gatech.edu/Publications/Public/TechReport/2020/louboutin2020SCsta/louboutin2020SCsta.html},
	author = {Mathias Louboutin and Fabio Luporini and Philipp A. Witte and Rhodri Nelson and George Bisbas and Jan Thorbecke and Felix J. Herrmann and Gerard Gorman}
}
@conference {moldoveanu2018SEGcsbdl,
	title = {Compressive sensing based design for land and OBS surveys: The noise issue},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2018},
	note = {(SEG, Anaheim)},
	month = {10},
	pages = {102-106},
	abstract = {Two important developments in seismic acquisition in the last decade were introduction of simultaneous shooting and reconstruction via inversion of the complete seismic wavefield. As the seismic wavefield is typically under sampled in at least one of the four spatial coordinates, both developments could contribute to improve the sampling and, in addition, to increase acquisition efficiency. Herein, we look at the impact of seismic noise when simultaneous shooting and seismic wavefield reconstruction are implemented for a land or ocean bottom seismic survey.},
	keywords = {SEG},
	doi = {10.1190/segam2018-2997759.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2018/moldoveanu2018SEGcsbdl/moldoveanu2018SEGcsbdl.pdf},
	author = {Nick Moldoveanu and Philip Bilsby and John Quigley and Rajiv Kumar and Felix J. Herrmann}
}
@presentation {bougher2016SINBADFaaa,
	title = {Amplitude vs. angle analysis as an unsupervised learning problem},
	year = {2016},
	publisher = {SINBAD},
	abstract = {Amplitude vs. angle analysis (AVA) of pre-stack seismic data is a commonly used method for inferring petrophysical information from seismic data. Conventionally, a two-term linearized rock physics model (Shuey equation) is used to invert angle-domain common-image gathers. Multivariate analysis of the inverted terms leads to a background of siliciclastic interfaces, where outlying points are associated with hydrocarbon saturated sands. The acquisition and processing of seismic data does not result in highly-calibrated measurements that adhere to the rock physics model, which often inhibits the success of AVA analysis. We offer an alternative approach that uses PCA-based methods to learn projections directly from the data without the need of a physical model. Results on synthetic and field data show that PCA-based projections can improve segmentation of potential reservoirs in seismic data.},
	keywords = {Presentation, SINBAD, SINBADFALL2016, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/bougher2016SINBADFaaa/bougher2016SINBADFaaa.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/bougher2016SINBADFaaa/bougher2016SINBADFaaa.mov},
	author = {Ben B. Bougher and Felix J. Herrmann}
}
@conference {bougher2016SEGava,
	title = {AVA classification as an unsupervised machine-learning problem},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2016},
	note = {(SEG, Dallas)},
	month = {10},
	pages = {553-556},
	abstract = {Much of AVA analysis relies on characterizing background trends and anomalies in pre-stack seismic data. Analysts reduce a seismic section into a small number of these trends and anomalies, suggesting that a low-dimensional structure can be inferred from the data. We describe AVA-attribute characterization as an unsupervised-learning problem, where AVA classes are learned directly from the data without any prior assumptions on physics and geological settings. The method is demonstrated on the Marmousi II elastic model, where a gas reservoir was successfully delineated from a background trend in a depth migrated image.},
	keywords = {AVA, machine learning, SEG},
	doi = {10.1190/segam2016-13874419.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/bougher2016SEGava/bougher2016SEGava.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/bougher2016SEGava/bougher2016SEGava_pres.pdf},
	author = {Ben B. Bougher and Felix J. Herrmann}
}
@mastersthesis {bougher2016THmla,
	title = {Machine learning applications to geophysical data analysis},
	year = {2016},
	note = {(MSc)},
	month = {08},
	school = {The University of British Columbia},
	type = {masters},
	address = {Vancouver},
	abstract = {The sedimentary layers of the Earth are a complex amorphous material formed from chaotic, turbulent, and random natural processes. Exploration geophysicists use a combination of assumptions, approximate physical models, and trained pattern recognition to extract useful information from complex remote sensing data such as seismic and well logs. In this thesis I investigate supervised and unsupervised machine learning models in geophysical data analysis and present two novel applications to exploration geophysics. Firstly, interpreted well logs from the Trenton-Black River study are used to train a classifier that results in a success rate of 67\% at predicting stratigraphic units from gamma ray logs. I use the scattering transform, a multiscale analysis transform, to extract discriminating features to feed a K-nearest neighbour classifier. A second experiment frames a conventional pre-stack seismic data characterization workflow as an unsupervised machine learning problem that is free from physical assumptions. Conventionally, the Shuey model is used to fit the angle dependent reflectivity response of seismic data. I instead use principle component based approaches to learn projections from the data that improve classification. Results on the Marmousi II elastic model and an industry field dataset show that unsupervised learning models can be effective at segmenting hydrocarbon reservoirs from seismic data.},
	keywords = {AVA, machine learning, MSc, PCA, scattering transform, thesis, well logs},
	url = {https://slim.gatech.edu/Publications/Public/Thesis/2016/bougher2016THmla/bougher2016THmla.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Thesis/2016/bougher2016THmla/bougher2016THmla_pres.pdf},
	author = {Ben B. Bougher}
}
@article {bougher2015CSEGust,
	title = {Using the scattering transform to predict stratigraphic units from well logs},
	journal = {CSEG Recorder},
	volume = {41},
	number = {1},
	year = {2016},
	note = {(CSEG Recorder)},
	month = {01},
	pages = {22-25},
	abstract = {Much of geophysical interpretation relies on trained pattern recognition of signals and images, a workflow that can be modeled by supervised machine learning. A challenge of supervised learning is determining a physically meaningful feature set that can successfully classify the data. Defined by a network of cascading wavelets, the scattering transform provides a non-linear multiscale analysis that has deep connections to the fractal statistics of the signal. Interestingly, the scattering transform takes the form of a pre-trained convolutional neural network. This paper uses the scattering transform to extract features from well logs in order to train a classifier that can predict stratigraphic units. The methodology is tested on interpreted well logs from Trenton-Black River project and initial results are presented.},
	keywords = {machine learning, scattering transform, well logs},
	url = {https://slim.gatech.edu/Publications/Public/Journals/CSEGRecorder/2016/bougher2015CSEGust/bougher2015CSEGust.html},
	url2 = {http://csegrecorder.com/articles/view/using-the-scattering-transform-to-predict-strat-units-from-well-logs},
	author = {Ben B. Bougher and Felix J. Herrmann}
}
@conference {bougher2015CSEGpsu,
	title = {Prediction of stratigraphic units from spectral co-occurance coefficients of well logs},
	booktitle = {CSEG Annual Conference Proceedings},
	year = {2015},
	note = {(CSEG, Calgary)},
	month = {05},
	abstract = {Well logging is the process of making physical measurements down bore holes in order to characterize geological and structural properties. Logs are visually interpreted and correlated to classify regions that are similar in structure, a process that can be modelled with machine learning. This project applies supervised learning methods to labelled well logs from the Trenton Black River data set in order to classify major stratigraphic units. Spectral co-occurance coefficients were used for feature extraction, and a k-nearest-neighbours approach was used for classification. This novel approach was applied to real field data in a high-impact domain, yielding promising results for future research.},
	keywords = {CSEG, machine learning, scattering transform, well logs},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/CSEG/2015/bougher2015CSEGpsu/bougher2015CSEGpsu.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/CSEG/2015/bougher2015CSEGpsu/bougher2015CSEGpsu_pres.pdf},
	author = {Ben B. Bougher and Felix J. Herrmann}
}
@conference {fang2015IIPFWIsoa,
	title = {Stochastic optimization and its application to seismic inversion},
	booktitle = {Inaugural Full-Waveform Inversion Workshop},
	year = {2015},
	note = {(Natal, Brazil)},
	month = {08-09},
	keywords = {FWI, inversion, Optimization, stochastic},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/IIPFWI/fang2015IIPFWIsoa/fang2015IIPFWIsoa_pres.pdf},
	author = {Zhilong Fang and Tim Burgess and Felix J. Herrmann}
}
@conference {aravkin2012SEGST,
	title = {On non-uniqueness of the Student{\textquoteright}s t-formulation for linear inverse problems},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {31},
	year = {2012},
	month = {11},
	pages = {1-5},
	publisher = {SEG},
	organization = {SEG},
	abstract = {We review the statistical interpretation of inverse problem formulations, and the motivations for selecting non-convex penalties for robust behaviour with respect to measurement outliers or artifacts in the data. An important downside of using non-convex formulations such as the Student{\textquoteright}s t is the potential for non-uniqueness, and we present a simple example where the Student{\textquoteright}s t penalty can be made to have many local minima by appropriately selecting the degrees of freedom parameter. On the other hand, the non-convexity of the Student{\textquoteright}s t is precisely what gives it the ability to ignore artifacts in the data. We explain this idea, and present a stylized imaging experiment, where the Student{\textquoteright}s t is able to recover a velocity perturbation from data contaminated by a very peculiar artifact {\textendash}- data from a different velocity perturbation. The performance of Student{\textquoteright}s t inversion is investigated empirically for different values of the degrees of freedom parameter, and different initial conditions.},
	keywords = {non-convex, robust, SEG, Student{\textquoteright}s t, uniqueness},
	doi = {10.1190/segam2012-1558.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2012/aravkin2012SEGST/aravkin2012SEGST.pdf},
	author = {Aleksandr Y. Aravkin and Tristan van Leeuwen and Kenneth Bube and Felix J. Herrmann}
}
@conference {aravkin2011SIAMfwi,
	title = {Full waveform inversion with compressive updates},
	booktitle = {SIAM},
	year = {2011},
	publisher = {SIAM CS\&E 2011},
	organization = {SIAM CS\&E 2011},
	abstract = {Full-waveform inversion relies on large multi-experiment data volumes. While improvements in acquisition and inversion have been extremely successful, the current push for higher quality models reveals fundamental shortcomings handling increasing problem sizes numerically. To address this fundamental issue, we propose a randomized dimensionality-reduction strategy motivated by recent developments in stochastic optimization and compressive sensing. In this formulation conventional Gauss-Newton iterations are replaced by dimensionality-reduced sparse recovery problems with source encodings.},
	keywords = {Full-waveform inversion, SLIM},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SIAM/2011/aravkin2011SIAMfwi/aravkin2011SIAMfwi.pdf},
	author = {Aleksandr Y. Aravkin and Felix J. Herrmann and Tristan van Leeuwen and James V. Burke and Xiang Li}
}
@conference {aravkin2011EAGEnspf,
	title = {A nonlinear sparsity promoting formulation and algorithm for full waveform inversion},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2011},
	month = {05},
	abstract = {Full Waveform Inversion (FWI) is a computational procedure to extract medium parameters from seismic data. FWI is typically formulated as a nonlinear least squares optimization problem, and various regularization techniques are used to guide the optimization because the problem is illposed. In this paper, we propose a novel sparse regularization which exploits the ability of curvelets to efficiently represent geophysical images. We then formulate a corresponding sparsity promoting constrained optimization problem, which we call Nonlinear Basis Pursuit Denoise (NBPDN) and present an algorithm to solve this problem to recover medium parameters. The utility of the NBPDN formulation and efficacy of the algorithm are demonstrated on a stylized cross-well experiment, where a sparse velocity perturbation is recovered with higher quality than the standard FWI formulation (solved with LBFGS). The NBPDN formulation and algorithm can recover the sparse perturbation even when the data volume is compressed to 5 \% of the original size using random superposition.},
	keywords = {EAGE, Full-waveform inversion, Optimization},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2011/aravkin11EAGEnspf/aravkin11EAGEnspf.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2011/aravkin11EAGEnspf/aravkin11EAGEnspf_pres.pdf},
	url2 = {http://earthdoc.eage.org/publication/publicationdetails/?publication=50199},
	author = {Aleksandr Y. Aravkin and James V. Burke and Felix J. Herrmann and Tristan van Leeuwen}
}
@conference {aravkin2011EAGEspfwi,
	title = {Sparsity promoting formulations and algorithms for FWI},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2011},
	keywords = {EAGE, Full-waveform inversion},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2011/Aravkin2011EAGEspfwi/Aravkin2011EAGEspfwi.pdf},
	author = {Aleksandr Y. Aravkin and Tristan van Leeuwen and James V. Burke and Felix J. Herrmann}
}
@conference {aravkin2011ICIAMspf,
	title = {Sparsity promoting formulations and algorithms for FWI},
	booktitle = {ICIAM},
	year = {2011},
	note = {Presented at AMP Medical and Seismic Imaging, 2011, Vancouver BC.},
	month = {07},
	publisher = {ICIAM 2011},
	organization = {ICIAM 2011},
	abstract = {Full Waveform Inversion (FWI) is a computational procedure to extract medium parameters from seismic data. FWI is typically formulated as a nonlinear least squares optimization problem, and various regularization techniques are used to guide the optimization because the problem is ill-posed. We propose a novel sparse regularization which exploits the ability of curvelets to efficiently represent geophysical images. We then formulate a corresponding sparsity promoting constrained optimization problem, which we solve using an open source algorithm. The techniques are applicable to any inverse problem where sparsity modeling is appropriate. We demonstrate the efficacy of the formulation on a toy example (stylized cross-well experiment) and on a realistic Seismic example (partial Marmoussi model). We also discuss the tradeoff between model fit and sparsity promotion, with a view to extend existing techniques for linear inverse problems to the case where the forward model is nonlinear.},
	keywords = {Full-waveform inversion, Optimization, Presentation, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ICIAM/2011/aravkin2011ICIAMspf/aravkin2011ICIAMspf_pres.pdf},
	author = {Aleksandr Y. Aravkin and Tristan van Leeuwen and James V. Burke and Felix J. Herrmann}
}
@article {herrmann2009GEOPcbm,
	title = {Curvelet-based migration preconditioning and scaling},
	journal = {Geophysics},
	volume = {74},
	year = {2009},
	month = {07-08},
	pages = {A41-A46},
	abstract = {The extremely large size of typical seismic imaging problems has been one of the major stumbling blocks for iterative techniques to attain accurate migration amplitudes. These iterative methods are important because they complement theoretical approaches that are hampered by difficulties to control problems such as finite-acquisition aperture, source-receiver frequency response, and directivity. To solve these problems, we apply preconditioning, which significantly improves convergence of least-squares migration. We discuss different levels of preconditioning that range from corrections for the order of the migration operator to corrections for spherical spreading, and position and reflector-dip dependent amplitude errors. While the first two corrections correspond to simple scalings in the Fourier and physical domain, the third correction requires phase-space (space spanned by location and dip) scaling, which we carry out with curvelets. We show that our combined preconditioner leads to a significant improvement of the convergence of least-squares {\textquoteleft}wave-equation{\textquoteright} migration on a line from the SEG AA{\textquoteright} salt model.},
	keywords = {Imaging, migration, SLIM},
	doi = {10.1190/1.3124753},
	url = {http://geophysics.geoscienceworld.org/content/74/4/A41.abstract},
	url2 = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2009/herrmann2009GEOPcbm/herrmann2009GEOPcbm.pdf},
	author = {Felix J. Herrmann and Cody R. Brown and Yogi A. Erlangga and Peyman P. Moghaddam}
}
@conference {moghaddam2008SEGcbm,
	title = {Curvelet-based migration preconditioning},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {27},
	year = {2008},
	month = {11},
	pages = {2211-2215},
	publisher = {SEG},
	organization = {SEG},
	abstract = {In this paper, we introduce a preconditioner for seismic imaging{\textendash}-i.e., the inversion of the linearized Born scattering operator. This preconditioner approximately corrects for the {\textquoteleft}{\textquoteleft}square root{\textquoteright}{\textquoteright} of the normal{\textendash}-i.e., the demigration-migration operator. This approach consists of three parts, namely (i) a left preconditoner, defined by a fractional time integration designed to make the migration operator zero order, and two right preconditioners that apply (ii) a scaling in the physical domain accounting for a spherical spreading, and (iii) a curvelet-domain scaling that corrects for spatial and reflector-dip dependent amplitude errors. We show that a combination of these preconditioners lead to a significant improvement of the convergence for iterative least-squares solutions to the seismic imaging problem based on reverse-time migration operators.},
	keywords = {Presentation, SEG, SLIM},
	doi = {10.1190/1.3059325},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2008/moghaddam08SEGcbm/moghaddam08SEGcbm.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2008/moghaddam08SEGcbm/moghaddam08SEGcbm_pres.pdf},
	author = {Peyman P. Moghaddam and Cody R. Brown and Felix J. Herrmann}
}
@conference {herrmann2008SINBADitc,
	title = {Introduction to compressive (wavefield) computation},
	booktitle = {SINBAD 2008},
	year = {2008},
	keywords = {Presentation, SINBAD, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2008/herrmann2008SINBADitc/herrmann2008SINBADitc.pdf},
	author = {Felix J. Herrmann and Yogi A. Erlangga and Tim T.Y. Lin and Cody R. Brown}
}
@article {herrmann2007GJInlp,
	title = {Non-linear primary-multiple separation with directional curvelet frames},
	journal = {Geophysical Journal International},
	volume = {170},
	number = {2},
	year = {2007},
	month = {08},
	pages = {781-799},
	abstract = {Predictive multiple suppression methods consist of two main steps: a prediction step, during which multiples are predicted from seismic data, and a primary-multiple separation step, during which the predicted multiples are {\textquoteright}matched{\textquoteright} with the true multiples in the data and subsequently removed. This second separation step, which we will call the estimation step, is crucial in practice: an incorrect separation will cause residual multiple energy in the result or may lead to a distortion of the primaries, or both. To reduce these adverse effects, a new transformed-domain method is proposed where primaries and multiples are separated rather than matched. This separation is carried out on the basis of differences in the multiscale and multidirectional characteristics of these two signal components. Our method uses the curvelet transform, which maps multidimensional data volumes into almost orthogonal localized multidimensional prototype waveforms that vary in directional and spatio-temporal content. Primaries-only and multiples-only signal components are recovered from the total data volume by a non-linear optimization scheme that is stable under noisy input data. During the optimization, the two signal components are separated by enhancing sparseness (through weighted l1-norms) in the transformed domain subject to fitting the observed data as the sum of the separated components to within a user-defined tolerance level. Whenever, during the optimization, the estimates for the primaries in the transformed domain correlate with the predictions for the multiples, the recovery of the coefficients for the estimated primaries will be suppressed while for regions where the correlation is small the method seeks the sparsest set of coefficients that represent the estimation for the primaries. Our algorithm does not seek a matched filter and as such it differs fundamentally from traditional adaptive subtraction methods. The method derives its stability from the sparseness obtained by a non-parametric (i.e. not depending on a parametrized physical model) multiscale and multidirectional overcomplete signal representation. This sparsity serves as prior information and allows for a Bayesian interpretation of our method during which the log-likelihood function is minimized while the two signal components are assumed to be given by a superposition of prototype waveforms, drawn independently from a probability function that is weighted by the predicted primaries and multiples. In this paper, the predictions are based on the data-driven surface-related multiple elimination method. Synthetic and field data examples show a clean separation leading to a considerable improvement in multiple suppression compared to the conventional method of adaptive matched filtering. This improved separation translates into an improved stack.},
	keywords = {Processing, signal separation, SLIM},
	doi = {10.1111/j.1365-246X.2007.03360.x},
	url = {https://slim.gatech.edu/Publications/Public/Journals/GeophysicalJournalInternational/2007/herrmann2007GJInlp/herrmann2007GJInlp.pdf},
	author = {Felix J. Herrmann and Urs Boeniger and D. J. Verschuur}
}
@manual {rossross07MNsda,
	title = {SLIMpy development and programming interface for seismic processing},
	year = {2007},
	abstract = {Inverse problems in (exploration) seismology are known for their large to very large scale. For instance, certain sparsity-promoting inversion techniques involve vectors that easily exceed unknowns while seismic imaging involves the construction and application of matrix-free discretized operators where single matrix-vector evaluations may require hours, days or even weeks on large compute clusters. For these reasons, software development in this field has remained the domain of highly technical codes programmed in low-level languages with little eye for easy development, code reuse and integration with (nonlinear) programs that solve inverse problems.Following ideas from the Symes{\textquoteright} Rice Vector Library and Bartlett{\textquoteright}s C++ object-oriented interface, Thyra, and Reduction/Transformation operators (both part of the Trilinos software package), we developed a software-development environment based on overloading. This environment provides a pathway from in-core prototype development to out-of-core and MPI {\textquoteright}production{\textquoteright} code with a high level of code reuse. This code reuse is accomplished by integrating the out-of-core and MPI functionality into the dynamic object-oriented programming language Python. This integration is implemented through operator overloading and allows for the development of a coordinate-free solver framework that (i) promotes code reuse; (ii) analyses the statements in an abstract syntax tree and (iii) generates executable statements. In the current implementation, we developed an interface to generate executable statements for the out-of-core unix-pipe based (seismic) processing package RSF-Madagascar (rsf.sf.net). The modular design allows for interfaces to other seismic processing packages and to in-core Python packages such as numpy. So far, the implementation overloads linear operators and element-wise reduction/transformation operators. We are planning extensions towards nonlinear operators and integration with existing (parallel) solver frameworks such as Trilinos.},
	keywords = {SLIM, software},
	url = {https://slim.gatech.edu/Software/SLIM/SLIMpy/},
	author = {Sean Ross-Ross and Henryk Modzelewski and Cody R. Brown and Felix J. Herrmann}
}
@conference {yarham2006SEGcgrr,
	title = {Curvelet-based ground roll removal},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {25},
	year = {2006},
	pages = {2777-2782},
	publisher = {SEG},
	organization = {SEG},
	abstract = {We have effectively identified and removed ground roll through a two-step process. The first step is to identify the major components of the ground roll through various methods including multiscale separation, directional or frequency filtering or by any other method that identifies the ground roll. Given this estimate for ground roll, the recorded signal is separated during the second step through a block-coordinate relaxation method that seeks the sparsest set for weighted curvelet coefficients of the ground roll and the sought-after reflectivity. The combination of these two methods allows us to separate out the ground roll signal while preserving the reflector information. Since our method is iterative, we have control of the separation process. We successfully tested our algorithm on a real data set with a complex ground roll and reflector structure.},
	keywords = {SEG},
	doi = {10.1190/1.2370101},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2006/yarham06SEGcgrr/yarham06SEGcgrr.pdf},
	author = {Carson Yarham and Urs Boeniger and Felix J. Herrmann}
}
@conference {herrmann2006SINBADpms,
	title = {Primary-multiple separation by curvelet frames},
	booktitle = {SINBAD 2006},
	volume = {170},
	year = {2006},
	pages = {781-799},
	publisher = {Geophysical Journal International},
	organization = {Geophysical Journal International},
	abstract = {Predictive multiple suppression methods consist of two main steps: a prediction step, during which multiples are predicted from seismic data, and a primary-multiple separation step, during which the predicted multiples are {\textquoteright}matched{\textquoteright} with the true multiples in the data and subsequently removed. The last step is crucial in practice: an incorrect separation will cause residual multiple energy in the result or may lead to a distortion of the primaries, or both. To reduce these adverse effects, a new transformed-domain method is proposed where primaries and multiples are separated rather than matched. This separation is carried out on the basis of differences in the multiscale and multidirectional characteristics of these two signal components. Our method uses the curvelet transform, which maps multidimensional data volumes into almost orthogonal localized multidimensional prototype waveforms that vary in directional and spatio-temporal content. Primaries-only and multiples-only signal components are recovered from the total data volume by a nonlinear optimization scheme that is stable under noisy input data. During the optimization, the two signal components are separated by enhancing sparseness (through weighted l1-norms) in the transformed domain subject to fitting the observed data as the sum of the separated components to within a user-defined tolerance level. Whenever the prediction for the two signal components in the transformed domain correlate, the recovery is suppressed while for regions where the correlation is small the method seeks the sparsest set of coefficients that represent each signal component. Our algorithm does not seek a matched filter and as such it differs fundamentally from traditional adaptive subtraction methods. The method derives its stability from the sparseness obtained by a non-parametric multiscale and multidirectional overcomplete signal representation. This sparsity serves as prior information and allows for a Bayesian interpretation of our method during which the log-likelihood function is minimized while the two signal components are assumed to be given by a superposition of prototype waveforms, drawn independently from a probability function that is weighted by the predicted primaries and multiples. In this paper, the predictions are based on the data-driven surface-related multiple elimination (SRME) method. Synthetic and field data examples show a clean separation leading to a considerable improvement in multiple suppression compared to the conventional method of adaptive matched filtering. This improved separation translates into an improved stack.},
	keywords = {Presentation, SINBAD, SLIM},
	doi = {10.1111/j.1365-246X.2007.03360.x},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2006/herrmann2006SINBADpms/herrmann2006SINBADpms.pdf},
	author = {Felix J. Herrmann and Urs Boeniger and D. J. Verschuur}
}
@conference {beyreuther2005SEGcot,
	title = {Computation of time-lapse differences with 3-D directional frames},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {24},
	year = {2005},
	pages = {2488-2491},
	publisher = {SEG},
	organization = {SEG},
	abstract = {We present an alternative method of extracting production related differences from time-lapse seismic data sets. Our method is not based on the actual subtraction of the two data sets, risking the enhancement of noise and introduction of artifacts due to local phase rotation and slightly misaligned events. Rather, it mutes events of the monitor survey with respect to the baseline survey based on the magnitudes of coefficients in a sparse and local atomic decomposition. Our technique is demonstrated to be an effective tool for enhancing the time-lapse signal from surveys which have been cross-equalized. \copyright2005 Society of Exploration Geophysicists},
	keywords = {SEG, SLIM},
	doi = {10.1190/1.2148227},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2005/Beyreuther05SEGcot/Beyreuther05SEGcot.pdf},
	author = {Moritz Beyreuther and Jamin Cristall and Felix J. Herrmann}
}
@conference {beyreuther2004EAGEcdo,
	title = {Curvelet denoising of 4-D seismic},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2004},
	month = {06},
	abstract = {With burgeoning world demand and a limited rate of discovery of new reserves, there is increasing impetus upon the industry to optimize recovery from already existing fields. 4D, or time-lapse, seismic imaging is an emerging technology that holds great promise to better monitor and optimise reservoir production. The basic idea behind 4D seismic is that when multiple 3D surveys are acquired at separate calendar times over a producing field, the reservoir geology will not change from survey to survey but the state of the reservoir fluids will change. Thus, taking the difference between two 3D surveys should remove the static geologic contribution to the data and isolate the time- varying fluid flow component. However, a major challenge in 4D seismic is that acquisition and processing differences between 3D surveys often overshadow the changes caused by fluid flow. This problem is compounded when 4D effects are sought to be derived from vintage 3D data sets that were not originally acquired with 4D in mind. The goal of this study is to remove the acquisition and imaging artefacts from a 4D seismic difference cube using Curveket processing techniques.},
	keywords = {EAGE, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2004/Beyreuther04EAGEcdo/beyreuther2004EAGEcdo_paper.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2004/Beyreuther04EAGEcdo/Beyreuther04EAGEcdo_pres.pdf},
	url2 = {https://circle.ubc.ca/bitstream/handle/2429/453/EAGE4D2004.pdf?sequence=1},
	author = {Moritz Beyreuther and Felix J. Herrmann and Jamin Cristall}
}
@conference {cristall2004CSEGcpa,
	title = {Curvelet processing and imaging: 4-D adaptive subtraction},
	booktitle = {CSEG Annual Conference Proceedings},
	year = {2004},
	month = {05},
	publisher = {CSEG},
	organization = {CSEG},
	abstract = {With burgeoning world demand and a limited rate of discovery of new reserves, there is increasing impetus upon the industry to optimize recovery from already existing fields. 4D, or time-lapse, seismic imaging holds great promise to better monitor and optimise reservoir production. The basic idea behind 4D seismic is that when multiple 3D surveys are acquired at separate calendar times over a producing field, the reservoir geology will not change from survey to survey but the state of the reservoir fluids will change. Thus, taking the difference between two 3D surveys should remove the static geologic contribution to the data and isolate the time-varying fluid flow component. However, a major challenge in 4D seismic is that acquisition and processing differences between 3D surveys often overshadow the changes caused by fluid flow. This problem is compounded when 4D effects are sought to be derived from legacy 3D data sets that were not originally acquired with 4D in mind. The goal of this study is to remove the acquisition and imaging artefacts from a 4D seismic difference cube using Curvelet processing techniques.},
	keywords = {SLIM},
	url = {http://www.cseg.ca/assets/files/resources/abstracts/2004/059S0201-Cristall_J_Curvelet_4D.pdf},
	author = {Jamin Cristall and Moritz Beyreuther and Felix J. Herrmann}
}
@article {bernabe2004JGRpas,
	title = {Permeability and storativity of binary mixtures of high-and low-porosity materials},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {109},
	year = {2004},
	month = {12},
	pages = {B12207},
	abstract = {As a first step toward determining the mixing laws for the transport properties of rocks, we prepared binary mixtures of high- and low-permeability materials by isostatically hot-pressing mixtures of fine powders of calcite and quartz. The resulting rocks were marbles containing varying concentrations of dispersed quartz grains. Pores were present throughout the rock, but the largest ones were preferentially associated with the quartz particles, leading us to characterize the material as being composed of two phases, one with high permeability and the second with low permeability. We measured the permeability and storativity of these materials using the oscillating flow technique, while systematically varying the effective pressure and the period and amplitude of the input fluid oscillation. Control measurements performed using the steady state flow and pulse decay techniques agreed well with the oscillating flow tests. The hydraulic properties of the marbles were highly sensitive to the volume fraction of the high-permeability phase (directly related to the quartz content). Below a critical quartz content, slightly less than 20 wt \%, the high-permeability volume elements were disconnected, and the overall permeability was low. Above the critical quartz content the high-permeability volume elements formed throughgoing paths, and permeability increased sharply. We numerically simulated fluid flow through binary materials and found that permeability approximately obeys a percolation-based mixing law, consistent with the measured permeability of the calcite-quartz aggregates.},
	keywords = {Modeling, permeability, porosity, SLIM},
	doi = {10.1029/2004JB003111},
	url = {https://slim.gatech.edu/Publications/Public/Journals/JournalOfGeophysicalResearch/2004/bernabe04JGRpas/bernabe04JGRpas.pdf},
	url2 = {http://onlinelibrary.wiley.com/doi/10.1029/2004JB003111/abstract},
	author = {Y. Bernab{\'e} and U. Mok and B. Evans and Felix J. Herrmann}
}
@article {herrmann2004GJIssa,
	title = {Seismic singularities at upper-mantle phase transitions: a site percolation model},
	journal = {Geophysical Journal International},
	volume = {159},
	number = {3},
	year = {2004},
	month = {12},
	pages = {949-960},
	abstract = {Mineralogical phase transitions are usually invoked to account for the sharpness of globally observed upper-mantle seismic discontinuities. We propose a percolation-based model for the elastic properties of the phase mixture in the coexistence regions associated with these transitions. The major consequence of the model is that the elastic moduli (but not the density) display a singularity at the percolation threshold of the high-pressure phase. This model not only explains the sharp but continuous change in seismic velocities across the phase transition, but also predicts its abruptness and scale invariance, which are characterized by a non-integral scale exponent. Using the receiver-function approach and new, powerful signal-processing techniques, we quantitatively determine the singularity exponent from recordings of converted seismic waves at two Australian stations (CAN and WRAB). Using the estimated values, we construct velocity{\textendash}depth profiles across the singularities and verify that the calculated converted waveforms match the observations under CAN. Finally, we point out a series of additional predictions that may provide new insights into the physics and fine structure of the upper-mantle transition zone.},
	keywords = {Modelling, percolation, SLIM},
	doi = {10.1111/j.1365-246X.2004.02464.x},
	author = {Felix J. Herrmann and Y. Bernab{\'e}}
}
