@article {orozco2023invnet,
	title = {InvertibleNetworks.jl: A Julia package for scalable normalizing flows},
	journal = {Journal of Open Source Software},
	volume = {9},
	year = {2024},
	month = {7},
	abstract = {InvertibleNetworks.jl is a Julia package designed for the scalable implementation of normalizing flows, a method for density estimation and sampling in high-dimensional distributions. This package excels in memory efficiency by leveraging the inherent invertibility of normalizing flows, which significantly reduces memory requirements during backpropagation compared to existing normalizing flow packages that rely on automatic differentiation frameworks. InvertibleNetworks.jl has been adapted for diverse applications, including seismic imaging, medical imaging, and CO2 monitoring, demonstrating its effectiveness in learning high-dimensional distributions.},
	keywords = {Bayesian inference, computing, conditional normalizing flows, deep learning, HPC, Inverse problems, memory, Normalizing flows, software, Uncertainty quantification},
	doi = {10.21105/joss.06554},
	url = {https://doi.org/10.21105/joss.06554},
	author = {Rafael Orozco and Philipp A. Witte and Mathias Louboutin and Ali Siahkoohi and Gabrio Rizzuti and Bas Peters and Felix J. Herrmann}
}
@article {louboutin2023lmi,
	title = {Learned multiphysics inversion with differentiable programming and machine learning},
	journal = {The Leading Edge},
	volume = {42},
	year = {2023},
	note = {(The Leading Edge)},
	month = {07},
	pages = {452-516},
	abstract = {We present the Seismic Laboratory for Imaging and Modeling/Monitoring (SLIM) open-source software framework for computational geophysics and, more generally, inverse problems involving the wave-equation (e.g., seismic and medical ultrasound), regularization with learned priors, and learned neural surrogates for multiphase flow simulations. By integrating multiple layers of abstraction, our software is designed to be both readable and scalable. This allows researchers to easily formulate their problems in an abstract fashion while exploiting the latest developments in high-performance computing. We illustrate and demonstrate our design principles and their benefits by means of building a scalable prototype for permeability inversion from time-lapse crosswell seismic data, which aside from coupling of wave physics and multiphase flow, involves machine learning.},
	keywords = {CCS, deep learning, Fourier neural operators, GCS, Imaging, inversion, invertible networks, machine learning, monitoring, multiphysics, Normalizing flows, SLIM, software, time-lapse},
	url = {https://library.seg.org/doi/10.1190/tle42070474.1},
	presentation = {https://slim.gatech.edu/Publications/Public/Journals/TheLeadingEdge/2023/louboutin2023lmi/le_software.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Journals/TheLeadingEdge/2023/louboutin2023lmi/le_software.html},
	author = {Mathias Louboutin and Ziyi Yin and Rafael Orozco and Thomas J. Grady II and Ali Siahkoohi and Gabrio Rizzuti and Philipp A. Witte and Olav M{\o}yner and Gerard J. Gorman and Felix J. Herrmann}
}
@article {grady2022SCtll,
	title = {Model-Parallel Fourier Neural Operators as Learned Surrogates for Large-Scale Parametric PDEs},
	journal = {Computers \& Geosciences},
	volume = {178},
	year = {2023},
	note = {(Computers and Geosciences)},
	month = {07},
	pages = {105402},
	abstract = {Fourier neural operators (FNOs) are a recently introduced neural network architecture for learning solution operators of partial differential equations (PDEs), which have been shown to perform significantly better than comparable deep learning approaches. Once trained, FNOs can achieve speed-ups of multiple orders of magnitude over conventional numerical PDE solvers. However, due to the high dimensionality of their input data and network weights, FNOs have so far only been applied to two-dimensional or small three-dimensional problems. To remove this limited problem-size barrier, we propose a model-parallel version of FNOs based on domain-decomposition of both the input data and network weights. We demonstrate that our model-parallel FNO is able to predict time-varying PDE solutions of over 2.6 billion variables on Perlmutter using up to 512 A100 GPUs and show an example of training a distributed FNO on the Azure cloud for simulating multiphase CO2 dynamics in the Earth{\textquoteright}s subsurface.},
	keywords = {CCS, Computers and Geosciences, deep learning, Fourier neural operators, HPC, large-scale, Model Parallelism, Multiphase Flow, Operator Learning},
	doi = {10.1016/j.cageo.2023.105402},
	url = {https://slim.gatech.edu/Publications/Public/Journals/ComputersAndGeosciences/2023/grady2022SCtll/grady2022SCtll.pdf},
	software = {https://github.com/slimgroup/dfno},
	author = {Thomas J. Grady II and Rishi Khan and Mathias Louboutin and Ziyi Yin and Philipp A. Witte and Ranveer Chandra and Russell J. Hewett and Felix J. Herrmann}
}
@conference {louboutin2022RHPCafa,
	title = {Abstractions for at-scale seismic inversion},
	booktitle = {Rice Oil and Gas High Performance Computing Conference 2022},
	year = {2022},
	note = {Rice Oil and Gas High Performance Computing Conference 2022},
	month = {03},
	pages = {Thursday Workshop: Devito Training and Hackathon},
	abstract = {We present the SLIM open-source software framework for computational geophysics, and more generally, inverse problems based on the wave-equation (e.g., medical ultrasound). We developed a software environment aimed at scalable research and development by designing multiple layers of abstractions. This environment allows the researchers to easily formulate their problem in an abstract fashion, while still being able to exploit the latest developments in high-performance computing. We illustrate and demonstrate the benefits of our software design on many geophysical applications, including seismic inversion and physics-informed machine learning for geophysics (e.g., loop unrolled imaging, uncertainty quantification), all while facilitating the integration of external software.},
	keywords = {CCS, devito, FWI, HPC, inversion, JUDI, machine learning, RHPC, software, Uncertainty quantification},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/RHPC/2022/louboutin2022RHPCafa/RiceHPC22.pdf},
	url2 = {https://youtu.be/scRTbP8w6Wk?t=4542},
	author = {Mathias Louboutin and Ali Siahkoohi and Ziyi Yin and Rafael Orozco and Thomas J. Grady II and Yijun Zhang and Philipp A. Witte and Gabrio Rizzuti and Felix J. Herrmann}
}
@conference {louboutin2022SEGais,
	title = {Accelerating innovation with software abstractions for scalable computational geophysics},
	booktitle = {International Meeting for Applied Geoscience and Energy Expanded Abstracts},
	year = {2022},
	note = {(IMAGE, Houston)},
	month = {05},
	abstract = {We present the SLIM open-source software framework for computational geophysics, and more generally, inverse problems based on the wave-equation (e.g., medical ultrasound). We developed a software environment aimed at scalable research and development by designing multiple layers of abstractions. This environment allows the researchers to easily formulate their problem in an abstract fashion, while still being able to exploit the latest developments in high-performance computing. We illustrate and demonstrate the benefits of our software design on many geophysical applications, including seismic inversion and physics-informed machine learning for geophysics(e.g., loop unrolled imaging, uncertainty quantification), all while facilitating the integration of external software.},
	keywords = {CCS, devito, FWI, HPC, inversion, JUDI, machine learning, SEG, software, Uncertainty quantification},
	doi = {10.1190/image2022-3750561.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2022/louboutin2022SEGais/louboutin_seg22.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2022/louboutin2022SEGais/index.html},
	author = {Mathias Louboutin and Philipp A. Witte and Ali Siahkoohi and Gabrio Rizzuti and Ziyi Yin and Rafael Orozco and Felix J. Herrmann}
}
@presentation {grady2022ML4SEISMICesn,
	title = {Effective scaling of numerical surrogates via domain-decomposed Fourier neural operators},
	year = {2022},
	month = {11},
	abstract = {Numerical surrogates are models which learn to mimic a complex physical process (such as the solution to a PDE produced by a solver) from a set of input/output pairs. Fourier neural operators (FNOs) are a specific type of numerical surrogate which use a learned matched filter to quickly approximate solutions to relatively smooth complex physical processes. In the case of carbon capture sequestration (CCS) technology, FNOs have been shown to well-approximate solutions to the two-phase flow equations, with speedups of 1,000 to 10,000 times at inference time versus a tradtitional solver. This speed combined with the fact that FNOs are differentiable with respect to their input parameters allows for inverse and uncertainty quantification problems to theoretically be solved on real 3D data, a previously intractible task. However, due to the size of the input data, network weights, and optimizer state, FNOs have thus far been limited to small to medium 2D and 3D problems, well below the size of an industry standard such as the Sleipner benchmark. Here we alleviate this problem by proposing a model-parallel FNO which makes use of domain decomposition of the input data and network weights, and exploits architectural features of FNOs to also include a natural form of asynchronous pipeline parallelism. Our network can scale to arbitrary problem sizes on CPU and GPU systems.},
	keywords = {CCS, deep learning, Fourier neural operators, HPC, ML4SEISMIC, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ML4SEISMIC/2022/grady2022ML4SEISMICesn/index.html},
	author = {Thomas J. Grady II and Rishi Khan and Mathias Louboutin and Ziyi Yin and Philipp A. Witte and Ranveer Chandra and Russell J. Hewett and Felix J. Herrmann}
}
@presentation {yin2022TRANSFORMjulia,
	title = {Julia for Geoscience},
	year = {2022},
	month = {04},
	abstract = {In this tutorial, we will introduce the Julia programming language to the geoscience community, covering topics such as I/O, data processing, inversion, and machine learning. We will begin by installing Julia and relevant packages. Through a series of tutorials, we will demonstrate Julia{\textquoteright}s abstraction power and show how to load and plot data, write your own functions/operators, form and solve a geophysical inverse problem, and demonstrate how to integrate wave-equation solvers in Julia with machine learning frameworks. The intent of this presentation is to provide an introductory level tutorial that will be useful to members of the geoscience community.},
	keywords = {JOLI, JUDI, Julia, machine learning, SLIM, software, TRANSFORM},
	url = {https://transform.softwareunderground.org/2022-julia-for-geoscience},
	url2 = {https://www.youtube.com/watch?v=HyWfp3NzIbg},
	software = {https://github.com/slimgroup/SLIMTutorials},
	author = {Ziyi Yin and Mathias Louboutin and Philipp A. Witte and Ali Siahkoohi and Gabrio Rizzuti and Rafael Orozco and Henryk Modzelewski and Felix J. Herrmann}
}
@article {kukreja2020lcc,
	title = {Lossy Checkpoint Compression in Full Waveform Inversion},
	journal = {Geoscientific Model Development},
	volume = {15},
	number = {9},
	year = {2022},
	note = {(Geoscientific Model Development)},
	month = {5},
	pages = {3815-3829},
	abstract = {This paper proposes a new method that combines check-pointing methods with error-controlled lossy compression for large-scale high-performance Full-Waveform Inversion (FWI), an inverse problem commonly used in geophysical exploration. This combination can significantly reduce data movement, allowing a reduction in run time as well as peak memory. In the Exascale computing era, frequent data transfer (e.g., memory bandwidth, PCIe bandwidth for GPUs, or network) is the performance bottleneck rather than the peak FLOPS of the processing unit. Like many other adjoint-based optimization problems, FWI is costly in terms of the number of floating-point operations, large memory footprint during backpropagation, and data transfer overheads. Past work for adjoint methods has developed checkpointing methods that reduce the peak memory requirements during backpropagation at the cost of additional floating-point computations. Combining this traditional checkpointing with error-controlled lossy compression, we explore the three-way tradeoff between memory, precision, and time to solution. We investigate how approximation errors introduced by lossy compression of the forward solution impact the objective function gradient and final inverted solution. Empirical results from these numerical experiments indicate that high lossy-compression rates (compression factors ranging up to 100) have a relatively minor impact on convergence rates and the quality of the final solution.},
	keywords = {checkpointing, full waveform inversion, GMD, Lossy compression, memory},
	doi = {10.5194/gmd-15-3815-2022},
	url = {https://arxiv.org/pdf/2009.12623.pdf},
	author = {Navjot Kukreja and Jan Hueckelheim and Mathias Louboutin and John Washbourne and Paul H. J. Kelly and Gerard J. Gorman}
}
@conference {herrmann2021KAUSTdbi,
	title = {Deep Bayesian Inference for Task-based Seismic Imaging},
	booktitle = {KAUST},
	year = {2021},
	note = {(KAUST, virtual)},
	month = {03},
	abstract = {High dimensionality, complex physics, and lack of access to the ground truth make seismic imaging one of the most challenging inversion problems in the computational imaging sciences. If these challenges were not bad enough, modern applications of computational imaging increasingly call for the assessment of uncertainty on the image itself and on subsequent tasks conducted on these images. By making use of a Markov-chain Monte Carlo (McMC) sampling technique, we demonstrate how uncertainty in the data can be propagated to the task of seismic horizon tracking. While this example shows that uncertainty can in principle be handled in a systematic manner by using neural nets as deep priors, it unfortunately also reveals that the McMC method fundamentally lacks ability to scale to relevant problem sizes. To overcome this shortcoming, we introduce a variational formulation based on normalizing flows. In this approach, invertible neural networks are trained to generate samples from the posterior. There are strong indications that this approach has the capability to scale better to problems where the forward operator is expensive to evaluate.},
	keywords = {Bayesian inference, Inverse problems, seismic imaging, Uncertainty quantification},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/KAUST/2021/herrmann2021KAUSTdbi/Kaust.pdf},
	author = {Ali Siahkoohi and Gabrio Rizzuti and Mathias Louboutin and Philipp A. Witte and Felix J. Herrmann}
}
@article {rizzuti2020dfw,
	title = {A dual formulation of wavefield reconstruction inversion for large-scale seismic inversion},
	journal = {Geophysics},
	volume = {86},
	number = {6},
	year = {2021},
	note = {(Geophysics)},
	month = {10},
	pages = {1ND-Z3},
	abstract = {Many of the seismic inversion techniques currently proposed that focus on robustness with respect to the background model choice are not apt to large-scale 3D applications, and the methods that are computationally feasible for industrial problems, such as full waveform inversion, are notoriously limited by convergence stagnation and require adequate starting models. We propose a novel solution that is both scalable and less sensitive to starting models or inaccurate parameters (such as anisotropy) that are typically kept fixed during inversion. It is based on a dual reformulation of the classical wavefield reconstruction inversion, whose empirical robustness with respect to these issues is well documented in the literature. While the classical version is not suited to 3D, as it leverages expensive frequency-domain solvers for the wave equation, our proposal allows the deployment of state-of-the-art time-domain finite-difference methods, and is potentially mature for industrial-scale problems.},
	keywords = {3D, finite-difference, Full-waveform inversion, wave equation},
	doi = {10.1190/geo2020-0743.1},
	url = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2021/rizzuti2020dfw/rizzuti2020dfw.html},
	author = {Gabrio Rizzuti and Mathias Louboutin and Rongrong Wang and Felix J. Herrmann}
}
@presentation {rizzuti2021ML4SEISMICdfw,
	title = {A dual formulation of wavefield reconstruction inversion for large-scale seismic inversion},
	year = {2021},
	month = {11},
	abstract = {Many of the seismic inversion techniques currently proposed that focus on robustness with respect to the background model choice are not apt to large-scale 3D applications, and the methods that are computationally feasible for industrial problems, such as full waveform inversion, are notoriously limited by convergence stagnation and require adequate starting models. We propose a novel solution that is both scalable and less sensitive to starting models or inaccurate parameters (such as anisotropy) that are typically kept fixed during inversion. It is based on a dual reformulation of the classical wavefield reconstruction inversion, whose empirical robustness with respect to these issues is well documented in the literature. While the classical version is not suited to 3D, as it leverages expensive frequency-domain solvers for the wave equation, our proposal allows the deployment of state-of-the-art time-domain finite-difference methods, and is potentially mature for industrial-scale problems.},
	keywords = {FWI, ML4SEISMIC, SLIM, wave-equation, WRI},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ML4SEISMIC/2021/rizzuti2021ML4SEISMICdfw/Tue-10-55-Rizzuti.pdf},
	author = {Gabrio Rizzuti and Mathias Louboutin and Rongrong Wang and Felix J. Herrmann}
}
@conference {siahkoohi2021EarthMLfar,
	title = {Fast and reliability-aware seismic imaging with conditional normalizing flows},
	booktitle = {Intelligent illumination of the Earth},
	year = {2021},
	note = {(KAUST, virtual)},
	month = {09},
	abstract = {The posterior probability distribution provides a comprehensive description of the solution in ill-posed inverse problems. Sampling from the posterior distribution in the context of seismic imaging is challenged by the high-dimensionality of the unknown and the expensive-to-evaluate forward operator. These challenges limit the applicability of Markov Chain sampling methods due to the costs associated with the forward operator. Moreover, explicitly choosing a prior distribution that captures the true heterogeneity exhibited by the Earth{\textquoteright}s subsurface further complicates casting seismic imaging into a Bayesian framework. To handle this situation and to assess uncertainty, we propose a data-driven variational inference approach based on conditional normalizing flows (NFs). The proposed scheme leverages existing data, which are in the form of low- and high-fidelity migrated image pairs, to train a conditional NF capable of characterizing the posterior distribution. After training, the NF can be used to sample from the posterior distribution associated with a previously unseen seismic survey, which is in some sense close, e.g., data from a neighboring survey area. In our numerical example, we obtain high-fidelity images from the Parihaka dataset and low-fidelity images are derived from these images through the process of demigration, followed by adding band-limited noise and migration. During inference, given shot records from a new neighboring seismic survey, we first compute the reverse-time migration image. Next, by feeding this low-fidelity migrated image to the NF we gain access to samples from the posterior distribution virtually for free. We use these samples to compute a high-fidelity image including a first assessment of the image{\textquoteright}s reliability.},
	keywords = {Normalizing flows, seismic imaging, Uncertainty quantification},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/KAUST/2021/siahkoohi2021EarthMLfar/siahkoohi2021EarthMLfar.pdf},
	author = {Ali Siahkoohi and Rafael Orozco and Gabrio Rizzuti and Philipp A. Witte and Mathias Louboutin and Felix J. Herrmann}
}
@conference {witte2021JULIACONmedlj,
	title = {InvertibleNetworks.jl - Memory efficient deep learning in Julia},
	booktitle = {JuliaCon},
	year = {2021},
	note = {(JuliaCon, virtual)},
	month = {07},
	abstract = {We present InvertibleNetworks.jl, an open-source package for invertible neural networks and normalizing flows using memory-efficient backpropagation. InvertibleNetworks.jl uses manually implement gradients to take advantage of the invertibility of building blocks, which allows for scaling to large-scale problem sizes. We present the architecture and features of the library and demonstrate its application to a variety of problems ranging from loop unrolling to uncertainty quantification.},
	keywords = {deep learning, Invertible network, Julia, Normalizing flows, segmentation},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/JuliaCon/2021/witte2021JULIACONmedlj/witte2021JULIACONmedlj.pdf},
	url2 = {https://www.youtube.com/watch?v=9M-zEGHY4i4},
	author = {Philipp A. Witte and Mathias Louboutin and Ali Siahkoohi and Gabrio Rizzuti and Bas Peters and Felix J. Herrmann}
}
@conference {herrmann2021Delftlwi,
	title = {Learned wave-based imaging - variational inference at scale},
	booktitle = {Delft},
	year = {2021},
	note = {(Delft, virtual)},
	month = {06},
	abstract = {High dimensionality, complex physics, and lack of access to the ground truth rank medical ultrasound and seismic imaging amongst the most challenging problems in the computational imaging sciences. If these challenges were not bad enough, modern applications of computational imaging increasingly call for the assessment of uncertainty on the image itself and on subsequent tasks. During this talk, I will show how recent developments in Normalizing Flows, a new type of invertible neural networks, can be used to cast wave-based imaging into a scalable Bayesian framework. Contrary to conventional methods, where sample images are drawn from the posterior distribution during inversion, our approach trains Normalizing Flows capable of generating samples from the posterior. Aside from greatly reducing the computational cost, this approach gives us access to the image itself (via Maximum a posteriori or mean estimation) and its multidimensional statistical distribution including its pointwise variance.},
	keywords = {medical imaging, Normalizing flows, seismic imaging, Uncertainty quantification},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/Delft/2021/herrmann2021Delftlwi/herrmann2021Delftlwi.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/Delft/2021/herrmann2021Delftlwi/herrmann2021Delftlwi.mp4},
	author = {Felix J. Herrmann and Ali Siahkoohi and Rafael Orozco and Gabrio Rizzuti and Philipp A. Witte and Mathias Louboutin}
}
@conference {herrmann2021SEGWSlts,
	title = {Low-cost time-lapse seismic imaging of CCS with the joint recovery model},
	booktitle = {SEG Workshop on Geophysical Challenges in Presalt Carbonates; virtual},
	year = {2021},
	note = {(IMAGE Workshop, virtual)},
	month = {10},
	keywords = {born, CCS, Imaging, JRM, multiple, SEG, workshop},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2021/herrmann2021SEGWSlts/herrmann2021SEGWSlts_pres.pdf},
	author = {Felix J. Herrmann and Mathias Louboutin and Ziyi Yin and Philipp A. Witte}
}
@techreport {louboutin2021NIPSmcte,
	title = {Low-memory stochastic backpropagation with multi-channel randomized trace estimation},
	number = {TR-CSE-2021-1},
	year = {2021},
	month = {06},
	abstract = {Thanks to the combination of state-of-the-art accelerators and highly optimized open software frameworks, there has been tremendous progress in the performance of deep neural networks. While these developments have been responsible for many breakthroughs, progress towards solving large-scale problems, such as video encoding and semantic segmentation in 3D, is hampered because access to on-premise memory is often limited. Instead of relying on (optimal) checkpointing or invertibility of the network layers{\textendash}-to recover the activations during backpropagation{\textendash}-we propose to approximate the gradient of convolutional layers in neural networks with a multi-channel randomized trace estimation technique. Compared to other methods, this approach is simple, amenable to analyses, and leads to a greatly reduced memory footprint. Even though the randomized trace estimation introduces stochasticity during training, we argue that this is of little consequence as long as the induced errors are of the same order as errors in the gradient due to the use of stochastic gradient descent. We discuss the performance of networks trained with stochastic backpropagation and how the error can be controlled while maximizing memory usage and minimizing computational overhead.},
	keywords = {Convolutions, HPC, Low Memory, machine learning, randomized linear algebra},
	url = {https://slim.gatech.edu/Publications/Public/TechReport/2021/louboutin2021NIPSmcte/louboutin2021NIPSmcte.pdf},
	software = {https://github.com/slimgroup/XConv},
	author = {Mathias Louboutin and Ali Siahkoohi and Rongrong Wang and Felix J. Herrmann}
}
@presentation {siahkoohi2021ML4SEISMICmcn,
	title = {Multifidelity conditional normalizing flows for physics-guided Bayesian inference},
	year = {2021},
	month = {11},
	abstract = {We introduce a scalable Bayesian inference approach that combines techniques from deep learning with a physic-based variational inference formulation. Bayesian inference for ill-posed inverse problems is challenged by the high-dimensionality of the unknown, computationally expensive forward operator, and choosing a prior distribution that accurately encodes prior knowledge on the unknown. To handle this situation and to assess uncertainty, we propose to approximate the posterior distribution using a pretrained conditional normalizing flow, which is trained on existing low- and high-fidelity estimations of the unknown. To further improve the accuracy of this approximation, we use transfer learning and finetune this normalizing flow by minimizing the Kullback-Leibler divergence between the predicted and the desired high-fidelity posterior density. This amounts to minimizing a physic-based variational inference objective with respect to the network weights, which we believe might scale better than Bayesian inference with Markov Chain sampling methods. We apply the proposed Bayesian inference approach to seismic imaging where we use quasi-real data obtained from the Parihaka dataset.},
	keywords = {deep learning, ML4SEISMIC, Normalizing flows, seismic imaging, SLIM, Variational Inference},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ML4SEISMIC/2021/siahkoohi2021ML4SEISMICmcn/Mon-11-45-Siahkoohi.pdf},
	author = {Ali Siahkoohi and Rafael Orozco and Gabrio Rizzuti and Philipp A. Witte and Mathias Louboutin and Felix J. Herrmann}
}
@conference {siahkoohi2021AABIpto,
	title = {Preconditioned training of normalizing flows for variational inference in inverse problems},
	booktitle = {3rd Symposium on Advances in Approximate Bayesian Inference},
	year = {2021},
	note = {AABI},
	month = {01},
	abstract = {Obtaining samples from the posterior distribution of inverse problems with expensive forward operators is challenging especially when the unknowns involve the strongly heterogeneous Earth. To meet these challenges, we propose a preconditioning scheme involving a conditional normalizing flow (NF) capable of sampling from a low-fidelity posterior distribution directly. This conditional NF is used to speed up the training of the high-fidelity objective involving minimization of the Kullback-Leibler divergence between the predicted and the desired high-fidelity posterior density for indirect measurements at hand. To minimize costs associated with the forward operator, we initialize the high-fidelity NF with the weights of the pretrained low-fidelity NF, which is trained beforehand on available model and data pairs. Our numerical experiments, including a 2D toy and a seismic compressed sensing example, demonstrate that thanks to the preconditioning considerable speed-ups are achievable compared to training NFs from scratch.},
	keywords = {Inverse problems, Normalizing flows, Uncertainty quantification},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/AABI/2021/siahkoohi2021AABIpto/siahkoohi2020ABIfab.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/AABI/2021/siahkoohi2021AABIpto/siahkoohi2021AABIpto_pres.pdf},
	url2 = {https://openreview.net/forum?id=P9m1sMaNQ8T},
	software = {https://github.com/slimgroup/Software.siahkoohi2021AABIpto},
	author = {Ali Siahkoohi and Gabrio Rizzuti and Mathias Louboutin and Philipp A. Witte and Felix J. Herrmann}
}
@presentation {witte2021ML4SEISMICrtc,
	title = {Redwood {\textendash} towards clusterless supercomputing in the cloud},
	year = {2021},
	month = {11},
	abstract = {We present Redwood, a Julia framework for clusterless supercomputing in the cloud. Redwood provides a set of distributed programming macros that enable users to remotely execute Julia functions in parallel through cloud services for batch and serverless computing. We present the architecture and design of Redwood, as well as its application to existing Julia packages for machine learning and inverse problems.},
	keywords = {cloud computing, clusterless, HPC, Julia, ML4SEISMIC, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ML4SEISMIC/2021/witte2021ML4SEISMICrtc/Mon-10-30-Witte.pdf},
	author = {Philipp A. Witte}
}
@conference {ren2021AGUsvi,
	title = {Seismic Velocity Inversion and Uncertainty Quantification Using Conditional Normalizing Flows},
	booktitle = {AGU Annual Meeting},
	year = {2021},
	note = {(AGU 2021, New Orleans)},
	month = {12},
	pages = {U12A-03},
	abstract = {In this work, we use a conditional normalizing flow (CNF) to address the seismic velocity inversion problem. Considering the large dimension difference between seismic data and velocity model, we reduce the data dimension by calculating its reverse time. After that, we train the CNF on pairs of migrated data and velocity. During inference, given a new seismic data, feeding the corresponding migrated image into the trained CNF will lead to posterior samples of the velocity inversion distribution. In addition, uncertainty quantification of the inverted results can be achieved by statistical metrics like mean and standard deviation. In our numerical example, the implementation is based on open-sourced software InvertibleNetworks.jl (Witte et al., 2021), JUDI.jl (Witte et al., 2019) and Devito (Louboutin et al., 2019). The training dataset are built based on the SEG/EAGE Overthrust model. For an unseen seismic data, the posterior samples of inversion results given by the trained CNF can be considered as good estimates of the true velocity. Especially, judging from the metrics like MAE, MSE, PSNR, SSIM, et al., the posterior mean is usually closer to the true velocity and the standard deviation indicates that the velocity value is more reliable within the subsurface layers than that on layer edges. Moreover, the inverted results, either the posterior samples or posterior mean, can be used as an initial model in the subsequent FWI for a more accurate result.},
	keywords = {AGU, inversion, machine learning, Normalizing flows, Uncertainty quantification},
	url = {https://agu.confex.com/agu/fm21/meetingapp.cgi/Paper/815883},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/AGU/2021/ren2021AGUsvi/ren2021AGUsvi_pres.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/AGU/2021/ren2021AGUsvi/ren2021AGUsvi.mp4},
	software = {https://github.com/slimgroup/INN_Velocity-Migration},
	author = {Yuxiao Ren and Philipp A. Witte and Ali Siahkoohi and Mathias Louboutin and Ziyi Yin and Felix J. Herrmann}
}
@conference {rizzuti2021SIAMGSwrid,
	title = {Time-domain Wavefield Reconstruction Inversion for large-scale seismic inversion},
	booktitle = {SIAM Conference on Mathematical and Computational Issues in the Geosciences},
	year = {2021},
	note = {(SIAM GS, virtual)},
	month = {06},
	abstract = {Data fitting based on a wave equation model has driven the modern development in seismic imaging, and it is now relevant for large scale inversion. A wider adoption is currently being limited by available computational resources and the classical limitations of non-linear optimization, which often produces suboptimal local minima. A class of methods based on the inexact solution of the wave equation, known as wavefield reconstruction inversion, have shown remarkable robustness towards local minima, however the solvers required for the inexact wave equation scale badly to 3D problems. We propose a reformulation of wavefield reconstruction inversion which involves the solution of the classical time-domain wave equation, and it is practical for large-scale problems. Applications to tilted-transverse isotropy synthetic datasets clearly demonstrate the advantage over traditional full-waveform inversion.},
	keywords = {FWI, SIAM, time domain, WRI},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SIAMGS/2021/rizzuti2021SIAMGSwrid/rizzuti2021SIAMGSwrid.pdf},
	author = {Gabrio Rizzuti and Mathias Louboutin and Rongrong Wang and Felix J. Herrmann}
}
@article {luporini2018aap,
	title = {Architecture and performance of Devito, a system for automated stencil computation},
	journal = {ACM Trans. Math. Softw.},
	volume = {46},
	number = {1},
	year = {2020},
	note = {(ACM Trans. Math. Softw.)},
	month = {04},
	abstract = {Stencil computations are a key part of many high-performance computing applications, such as image processing, convolutional neural networks, and finite-difference solvers for partial differential equations. Devito is a framework capable of generating highly-optimized code given symbolic equations expressed in Python, specialized in, but not limited to, affine (stencil) codes. The lowering process {\textendash} from mathematical equations down to C++ code {\textendash} is performed by the Devito compiler through a series of intermediate representations. Several performance optimizations are introduced, including advanced common sub-expressions elimination, tiling and parallelization. Some of these are obtained through well-established stencil optimizers, integrated in the back-end of the Devito compiler. The architecture of the Devito compiler, as well as the performance optimizations that are applied when generating code, are presented. The effectiveness of such performance optimizations is demonstrated using operators drawn from seismic imaging applications.},
	keywords = {compiler, finite difference method, performance optimization, stencil, structured grid, symbolic processing},
	doi = {10.1145/3374916},
	url = {https://slim.gatech.edu/Publications/Public/Journals/ACMTOMS/2020/luporini2018aap/luporini2018aap.pdf},
	author = {Fabio Luporini and Mathias Louboutin and Michael Lange and Navjot Kukreja and Philipp A. Witte and Jan H{\"u}ckelheim and Charles Yount and Paul H. J. Kelly and Felix J. Herrmann and Gerard J. Gorman}
}
@article {witte2019TPDedas,
	title = {An Event-Driven Approach to Serverless Seismic Imaging in the Cloud},
	journal = {IEEE Transactions on Parallel and Distributed Systems},
	volume = {31},
	number = {9},
	year = {2020},
	note = {(IEEE Transactions on Parallel and Distributed Systems)},
	month = {03},
	pages = {2032-2049},
	abstract = {Adapting the cloud for high-performance computing (HPC) is a challenging task, as software for HPC applications hinges on fast network connections and is sensitive to hardware failures. Using cloud infrastructure to recreate conventional HPC clusters is therefore in many cases an infeasible solution for migrating HPC applications to the cloud. As an alternative to the generic lift and shift approach, we consider the specific application of seismic imaging and demonstrate a serverless and event-driven pproach for running large-scale instances of this problem in the cloud. Instead of permanently running compute instances, our workflow is based on a serverless architecture with high throughput batch computing and event-driven computations, in which computational resources are only running as long as they are utilized. We demonstrate that this approach is very flexible and allows for resilient and nested levels of parallelization, including domain decomposition for solving the underlying partial differential equations. While the event-driven approach introduces some overhead as computational resources are repeatedly restarted, it inherently provides resilience to instance shut-downs and allows a significant reduction of cost by avoiding idle instances, thus making the cloud a viable alternative to on-premise clusters for large-scale seismic imaging.},
	keywords = {cloud, event-driven, Imaging, LSRTM, serverless},
	doi = {10.1109/TPDS.2020.2982626},
	url = {https://slim.gatech.edu/Publications/Public/Journals/IEEETPDS/2020/witte2019TPDedas/witte2019TPDedas.html},
	author = {Philipp A. Witte and Mathias Louboutin and Henryk Modzelewski and Charles Jones and James Selvage and Felix J. Herrmann}
}
@conference {yin2020SEGesi,
	title = {Extended source imaging {\textendash}- a unifying framework for seismic and medical imaging},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2020},
	note = {(SEG, virtual)},
	month = {09},
	pages = {3502-3506},
	abstract = {We present three imaging modalities that live on the crossroads of seismic and medical imaging. Through the lens of extended source imaging, we can draw deep connections among the fields of wave-equation based seismic and medical imaging, despite first appearances. From the seismic perspective, we underline the importance to work with the correct physics and spatially varying velocity fields. Medical imaging, on the other hand, opens the possibility for new imaging modalities where outside stimuli, such as laser or radar pulses, can not only be used to identify endogenous optical or thermal contrasts but that these sources can also be used to insonify the medium so that images of the whole specimen can in principle be created.},
	keywords = {medical imaging, SEG, seismic imaging, variable projection},
	doi = {10.1190/segam2020-3426999.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2020/yin2020SEGesi/yin2020SEGesi.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2020/yin2020SEGesi/yin2020SEGesi_pres.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2020/yin2020SEGesi/yin2020SEGesi_pres.mp4},
	software = {https://github.com/slimgroup/Software.SEG2020},
	author = {Ziyi Yin and Rafael Orozco and Philipp A. Witte and Mathias Louboutin and Gabrio Rizzuti and Felix J. Herrmann}
}
@techreport {siahkoohi2020TRfuqf,
	title = {Faster Uncertainty Quantification for Inverse Problems with Conditional Normalizing Flows},
	number = {TR-CSE-2020-2},
	year = {2020},
	month = {07},
	institution = {Georgia Institute of Technology},
	abstract = {In inverse problems, we often have access to data consisting of paired samples $(x,y)\sim p_X,Y(x,y)$ where $y$ are partial observations of a physical system, and $x$ represents the unknowns of the problem. Under these circumstances, we can employ supervised training to learn a solution $x$ and its uncertainty from the observations $y$. We refer to this problem as the "supervised" case. However, the data $y\sim p_Y(y)$ collected at one point could be distributed differently than observations $y{\textquoteright}\sim p_Y{\textquoteright}(y{\textquoteright})$, relevant for a current set of problems. In the context of Bayesian inference, we propose a two-step scheme, which makes use of normalizing flows and joint data to train a conditional generator $q_θ(x|y)$ to approximate the target posterior density $p_X|Y(x|y)$. Additionally, this preliminary phase provides a density function $q_θ(x|y)$, which can be recast as a prior for the "unsupervised" problem, e.g. when only the observations $y{\textquoteright}\sim p_Y{\textquoteright}(y{\textquoteright})$, a likelihood model $y{\textquoteright}|x$, and a prior on $x{\textquoteright}$ are known. We then train another invertible generator with output density $q{\textquoteright}_φ(x|y{\textquoteright})$ specifically for $y{\textquoteright}$, allowing us to sample from the posterior $p_X|Y{\textquoteright}(x|y{\textquoteright})$. We present some synthetic results that demonstrate considerable training speedup when reusing the pretrained network $q_θ(x|y{\textquoteright})$ as a warm start or preconditioning for approximating $p_X|Y{\textquoteright}(x|y{\textquoteright})$, instead of learning from scratch. This training modality can be interpreted as an instance of transfer learning. This result is particularly relevant for large-scale inverse problems that employ expensive numerical simulations.},
	keywords = {deep learning, invertible networks, Uncertainty quantification},
	url = {https://slim.gatech.edu/Publications/Public/TechReport/2020/siahkoohi2020TRfuqf/siahkoohi2020TRfuqf.html},
	author = {Ali Siahkoohi and Gabrio Rizzuti and Philipp A. Witte and Felix J. Herrmann}
}
@conference {rizzuti2020SEGuqavp,
	title = {Parameterizing uncertainty by deep invertible networks, an application to reservoir characterization},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2020},
	note = {(SEG, virtual)},
	month = {09},
	pages = {1541-1545},
	abstract = {Uncertainty quantification for full-waveform inversion provides a probabilistic characterization of the ill-conditioning of the problem, comprising the sensitivity of the solution with respect to the starting model and data noise. This analysis allows to assess the confidence in the candidate solution and how it is reflected in the tasks that are typically performed after imaging (e.g., stratigraphic segmentation following reservoir characterization). Classically, uncertainty comes in the form of a probability distribution formulated from Bayesian principles, from which we seek to obtain samples. A popular solution involves Monte Carlo sampling. Here, we propose instead an approach characterized by training a deep network that "pushes forward" Gaussian random inputs into the model space (representing, for example, density or velocity) as if they were sampled from the actual posterior distribution. Such network is designed to solve a variational optimization problem based on the Kullback-Leibler divergence between the posterior and the network output distributions. This work is fundamentally rooted in recent developments for invertible networks. Special invertible architectures, besides being computational advantageous with respect to traditional networks, do also enable analytic computation of the output density function. Therefore, after training, these networks can be readily used as a new prior for a related inversion problem. This stands in stark contrast with Monte-Carlo methods, which only produce samples. We validate these ideas with an application to angle-versus-ray parameter analysis for reservoir characterization.},
	keywords = {Full-waveform inversion, SEG, Uncertainty quantification},
	doi = {10.1190/segam2020-3428150.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2020/rizzuti2020SEGuqavp/rizzuti2020SEGuqavp.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2020/rizzuti2020SEGuqavp/rizzuti2020SEGuqavp_pres.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2020/rizzuti2020SEGuqavp/rizzuti2020SEGuqavp_pres.mp4},
	software = {https://github.com/slimgroup/Software.SEG2020},
	author = {Gabrio Rizzuti and Ali Siahkoohi and Philipp A. Witte and Felix J. Herrmann}
}
@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 {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}
}
@conference {witte2019RHPCssi,
	title = {Serverless seismic imaging in the cloud},
	booktitle = {Rice Oil and Gas High Performance Computing Conference 2020},
	year = {2020},
	note = {Rice Oil and Gas High Performance Computing Conference 2020},
	month = {03},
	abstract = {This talk presents a serverless approach to seismic imaging in the cloud based on high-throughput containerized batch processing, event-driven computations and a domain-specific language compiler for solving the underlying wave equations. A 3D case study on Azure demonstrates that this approach allows reducing the operating cost of up to a factor of 6, making the cloud a viable alternative to on-premise HPC clusters for seismic imaging.},
	keywords = {cloud, HPC, reverse-time migration, serverless},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/RHPC/2019/witte2019RHPCssi/witte2019RHPCssi.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/RHPC/2019/witte2019RHPCssi/witte2019RHPCssi_pres.pdf},
	url2 = {https://www.youtube.com/watch?v=epilCpa078A\&list=PLcsG4X8Zn_UA3m7tgIMbOdXD8eOEPiUYD\&index=17},
	author = {Philipp A. Witte and Mathias Louboutin and Charles Jones and Felix J. Herrmann}
}
@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}
}
@article {yang2020tdsp,
	title = {Time-domain sparsity promoting least-squares reverse time migration with source estimation},
	journal = {Geophysical Prospecting},
	volume = {68},
	number = {9},
	year = {2020},
	note = {(Geophysical Prospecting)},
	month = {08},
	pages = {2697-2711},
	abstract = {Least-squares reverse time migration is well-known for its capability to generate artifact-free true-amplitude subsurface images through fitting observed data in the least-squares sense. However, when applied to realistic imaging problems, this approach is faced with issues related to overfitting and excessive computational costs induced by many wave-equation solves. The fact that the source function is unknown complicates this situation even further. Motivated by recent results in stochastic optimization and transform-domain sparsity-promotion, we demonstrate that the computational costs of inversion can be reduced significantly while avoiding imaging artifacts and restoring amplitudes. While powerfull, these new approaches do require accurate information on the source-time function, which is often lacking. Without this information, the imaging quality deteriorates rapidly. We address this issue by presenting an approach where the source-time function is estimated on the fly through a technique known as variable projection. Aside from introducing negligible computational overhead, the proposed method is shown to perform well on imaging problems with noisy data and problems that involve complex settings such as salt. In either case, the presented method produces high resolution high-amplitude fidelity images including an estimates for the source-time function. In addition, due to its use of stochastic optimization, we arrive at these images at roughly one to two times the cost of conventional reverse time migration involving all data.},
	keywords = {penalty, source estimation, sparsity inversion},
	doi = {10.1111/1365-2478.13021},
	url = {https://slim.gatech.edu/Publications/Public/Journals/GeophysicalProspecting/2020/yang2020tdsp/yang2020tdsp.html},
	author = {Mengmeng Yang and Zhilong Fang and Philipp A. Witte and Felix J. Herrmann}
}
@conference {rizzuti2020EAGEtwri,
	title = {Time-domain wavefield reconstruction inversion for large-scale seismics},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2020},
	note = {Accepted in EAGE, withdrawn due to logistical challenge},
	month = {1},
	abstract = {Wavefield reconstruction inversion is an imaging technique akin to full-waveform inversion, albeit based on a relaxed version of the wave equation. This relaxation aims to beat the multimodality typical of full-waveform inversion. However it prevents the use of time-marching solvers for the augmented equation and, as a consequence, cannot be straightforwardly employed to large 3D problems. In this work, we formulate a dual version of wavefield reconstruction inversion amenable to explicit time-domain solvers, yielding a robust and scalable inversion technique.},
	keywords = {3D, EAGE, Full-waveform inversion, Time-domain},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2020/rizzuti2020EAGEtwri/rizzuti2020EAGEtwri.html},
	author = {Gabrio Rizzuti and Mathias Louboutin and Rongrong Wang and Felix J. Herrmann}
}
@conference {herrmann2019EAGEHPCaii,
	title = {Accelerating ideation and innovation cheaply in the Cloud the power of abstraction, collaboration and reproducibility},
	booktitle = {4th EAGE Workshop on High-performance Computing},
	year = {2019},
	note = {(EAGE HPC Workshop, Dubai)},
	month = {10},
	keywords = {cloud, devito, EAGE, HPC, JUDI},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGEHPC/2019/herrmann2019EAGEHPCaii/herrmann2019EAGEHPCaii_pres.pdf},
	author = {Felix J. Herrmann and Charles Jones and Gerard Gorman and Jan H{\"u}ckelheim and Keegan Lensink and Paul H. J. Kelly and Navjot Kukreja and Henryk Modzelewski and Michael Lange and Mathias Louboutin and Fabio Luporini and James Selvages and Philipp A. Witte}
}
@conference {witte2019SIAMCSEclsmft,
	title = {Compressive least squares migration with on-the-fly Fourier transforms},
	booktitle = {SIAM Conference on Computational Science and Engineering},
	year = {2019},
	note = {(SIAM CSE)},
	month = {03},
	abstract = {Least-squares seismic imaging is an inversion-based approach for accurately imaging the earth{\textquoteright}s subsurface. However, in the time-domain, the computational cost and memory requirements of this approach scale with the size and recording length of the seismic experiment, thus making this approach often prohibitively expensive in practice. To overcome these issues, we borrow ideas from compressive sensing and signal processing and introduce an algorithm for sparsity-promoting seismic imaging using on-the-fly Fourier transforms. By computing gradients and functions values for random subsets of source locations and frequencies, we considerably limit the number of wave equation solves, while on-the-fly Fourier transforms allow computing an arbitrary number of monochromatic frequency-domain wavefields with a time-domain modeling code and without having to solve large-scale Helmholtz equations. The memory requirements of this approach are independent of the number of time steps and solely depend on the number of frequencies, which determine the amount of crosstalk and subsampling artifacts in the image. We show the application of our approach to several large-scale open source data sets and compare the results to a conventional time-domain approach with optimal checkpointing.},
	keywords = {Fourier, Imaging, least-squares, migration, on-the-fly},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SIAMCSE/2019/witte2019SIAMCSEclsmft/witte2019SIAMCSEclsmft.pdf},
	author = {Philipp A. Witte and Mathias Louboutin and Fabio Luporini and Gerard J. Gorman and Felix J. Herrmann}
}
@article {witte2018cls,
	title = {Compressive least-squares migration with on-the-fly Fourier transforms},
	journal = {Geophysics},
	volume = {84},
	number = {5},
	year = {2019},
	note = {(Geophysics)},
	month = {08},
	pages = {R655-R672},
	abstract = {Least-squares reverse-time migration is a powerful approach for true amplitude seismic imaging of complex geological structures, but the successful application of this method is currently hindered by its enormous computational cost, as well as high memory requirements for computing the gradient of the objective function. We tackle these problems by introducing an algorithm for low-cost sparsity-promoting least-squares migration using on-the-fly Fourier transforms. We formulate the least-squares migration objective function in the frequency domain and compute gradients for randomized subsets of shot records and frequencies, thus significantly reducing data movement and the number of overall wave equations solves. By using on-the-fly Fourier transforms, we can compute an arbitrary number of monochromatic frequency-domain wavefields with a time-domain modeling code, instead of having to solve individual Helmholtz equations for each frequency, which quickly becomes computationally infeasible when moving to high frequencies. Our numerical examples demonstrate that compressive imaging with on-the-fly Fourier transforms provides a fast and memory-efficient alternative to time-domain imaging with optimal checkpointing, whose memory requirements for a fixed background model and source wavelet is independent of the number of time steps. Instead, memory and additional computational cost grow with the number of frequencies and determine the amount of subsampling artifacts and crosstalk. In contrast to optimal checkpointing, this offers the possibility to trade both memory and computational cost for image quality or a larger number of iterations and is advantageous in new computing environments such as the cloud, where compute is often cheaper than memory and data movement.},
	keywords = {Fourier, least squares migration, sparsity-promotion},
	doi = {10.1190/geo2018-0490.1},
	url = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2019/witte2018cls/witte2018cls.pdf},
	author = {Philipp A. Witte and Mathias Louboutin and Fabio Luporini and Gerard J. Gorman and Felix J. Herrmann}
}
@article {louboutin2018dae,
	title = {Devito (v3.1.0): an embedded domain-specific language for finite differences and geophysical exploration},
	journal = {Geoscientific Model Development},
	year = {2019},
	note = {(Geoscientific Model Development)},
	abstract = {We introduce Devito, a new domain-specific language for implementing high-performance finite difference partial differential equation solvers. The motivating application is exploration seismology where methods such as Full-Waveform Inversion and Reverse-Time Migration are used to invert terabytes of seismic data to create images of the earth{\textquoteright}s subsurface. Even using modern supercomputers, it can take weeks to process a single seismic survey and create a useful subsurface image. The computational cost is dominated by the numerical solution of wave equations and their corresponding adjoints. Therefore, a great deal of effort is invested in aggressively optimizing the performance of these wave-equation propagators for different computer architectures. Additionally, the actual set of partial differential equations being solved and their numerical discretization is under constant innovation as increasingly realistic representations of the physics are developed, further ratcheting up the cost of practical solvers. By embedding a domain-specific language within Python and making heavy use of SymPy, a symbolic mathematics library, we make it possible to develop finite difference simulators quickly using a syntax that strongly resembles the mathematics. The Devito compiler reads this code and applies a wide range of analysis to generate highly optimized and parallel code. This approach can reduce the development time of a verified and optimized solver from months to days.},
	keywords = {finite-differences, HPC, Modeling, wave-equation},
	doi = {10.5194/gmd-12-1165-2019},
	url = {https://slim.gatech.edu/Publications/Public/Journals/GMD/2019/louboutin2018dae/louboutin2018dae.pdf},
	author = {Mathias Louboutin and Michael Lange and Fabio Luporini and Navjot Kukreja and Philipp A. Witte and Felix J. Herrmann and Paulius Velesko and Gerard J. Gorman}
}
@presentation {witte2019HOTCSEdsagip,
	title = {Domain-specific abstractions for large-scale geophysical inverse problems},
	year = {2019},
	note = {(HotCSE)},
	abstract = {During these times of sustained low oil prices, it is essential to look for new innovative ways to collect (time-lapse) seismic data at reduced costs and preferably also at reduced environmental impact. By now, there is an increasing body of corroborating evidence {\textemdash} whether these are simulated case studies or actual acquisitions on land and marine {\textemdash} that seismic acquisition based on the principles of compressive sensing delivers on this premise by removing the need to acquire replicated dense surveys. Up to ten-fold increases in acquisition efficiency have been reported by industry while there are indications that this breakthrough is only the beginning of a paradigm shift where full-azimuth time-lapse processing will become a reality. To familiarize the audience with this new technology, I will first describe the basics of compressive sensing, how it relates to missing-trace interpolation and simultaneous source acquisition, followed by how this technology is driving innovations in full-azimuth (time-lapse) acquisition, yielding high-fidelity data with a high degree of repeatability and at a fraction of the costs.},
	keywords = {full-waveform-inversion, Imaging, inversion, Julia, large-scale, software},
	url = {https://slim.gatech.edu/Publications/Public/Lectures/HotCSE/2019/witte2019HOTCSEdsagip/witte2019HOTCSEdsagip.pdf},
	author = {Philipp A. Witte and Mathias Louboutin and Felix J. Herrmann}
}
@conference {rizzuti2019SEGadf,
	title = {A dual formulation for time-domain wavefield reconstruction inversion},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2019},
	note = {(SEG, San Antonio)},
	month = {09},
	pages = {1480-1485},
	abstract = {We illustrate a dual formulation for full-waveform inversion potentially apt to large 3-D problems. It is based on the optimization of the wave equation compliance, under the constraint of data misfit not exceeding a prescribed noise level. In the Lagrangian formulation, model and wavefield state variables are complemented with multipliers having the same dimension of data ("dual data" variables). Analogously to classical wavefield reconstruction inversion, the wavefield unknowns can be projected out in closed form, by solving a version of the augmented wave equation. This leads to a saddle-point problem whose variables are only model and dual data. As such, this formulation represents a model extension, and is potentially robust against local minima. The classical downsides of model extension methods and wavefield reconstruction inversion are here effectively mitigated: storage of the dual variables is affordable, the augmented wave equation is amenable to time-marching finite-difference schemes, and no continuation strategy for penalty parameters is needed, with the prospect of 3-D applications.},
	keywords = {3D, Full-waveform inversion, SEG, Time-domain},
	doi = {10.1190/segam2019-3216760.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2019/rizzuti2019SEGadf/rizzuti2019SEGadf.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2019/rizzuti2019SEGadf/rizzuti2019SEGadf_pres.pdf},
	software = {https://github.com/slimgroup/Software.rizzuti2019SEGadf},
	author = {Gabrio Rizzuti and Mathias Louboutin and Rongrong Wang and Emmanouil Daskalakis and Felix J. Herrmann}
}
@conference {rizzuti2019SIAMGStwri,
	title = {A dual formulation for time-domain wavefield reconstruction inversion},
	booktitle = {SIAM Conference on Mathematical and Computational Issues in the Geosciences},
	year = {2019},
	note = {(SIAM GS)},
	month = {03},
	abstract = {Wavefield reconstruction inversion (WRI) is based on the combination of data misfit and PDE-residual penalty terms, and it aims at the reconstruction of the optimal property model and field pair. A full-space optimization would require the storage of field variables for every frequency or time sample (and every source). Remedies via variable projection are well studied, but they entail the solution of a modified wave equation, which is costly both in frequency and time domain. Here, we explore an efficient alternative to WRI, where properties are complemented, instead, by elements of the dual data space. The starting point of our proposal is the denoising reformulation of WRI: the PDE compliance is optimized under the constraint of a given data fit. Since the field variable can be eliminated from the associated saddle-point problem, we obtain a new objective functional of only properties and Lagrange multipliers. Its evaluation, and gradient computation, involves the solution of the original wave equation, which is amenable to time-domain discretization. The optimization scheme of choice can now leverage on affordable storage of the iterates.},
	keywords = {WRI},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SIAMGS/2019/rizzuti2019SIAMGStwri/rizzuti2019SIAMGStwri.pdf},
	author = {Gabrio Rizzuti and Mathias Louboutin and Rongrong Wang and Emmanouil Daskalakis and Felix J. Herrmann}
}
@conference {witte2019SEGedw,
	title = {Event-driven workflows for large-scale seismic imaging in the cloud},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2019},
	note = {(SEG, San Antonio)},
	month = {09},
	pages = {3984-3988},
	abstract = {Cloud computing has seen a large rise in popularity in recent years and is becoming a cost effective alternative to on-premise computing, with theoretically unlimited scalability. However, so far little progress has been made in adapting the cloud for high performance computing (HPC) tasks, such as seismic imaging and inversion. As the cloud does not provide the same type of fast and reliable connections as conventional HPC clusters, porting legacy codes developed for HPC environments to the cloud is ineffective and misses out on an opportunity to take advantage of new technologies presented by the cloud. We present a novel approach of bringing seismic imaging and inversion workflows to the cloud, which does not rely on a traditional HPC environment, but is based on serverless and event-driven computations. Computational resources are assigned dynamically in response to events, thus minimizing idle time and providing resilience to hardware failures. We test our workflow on two large-scale imaging examples and demonstrate that cost-effective HPC in the cloud is possible, but requires careful reconsiderations of how to bring software to the cloud.},
	keywords = {cloud, Imaging, large-scale, LS-RTM, RTM, SEG, workflow},
	doi = {10.1190/segam2019-3215069.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2019/witte2019SEGedw/witte2019SEGedw.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2019/witte2019SEGedw/witte2019SEGedw_pres.pdf},
	author = {Philipp A. Witte and Mathias Louboutin and Henryk Modzelewski and Charles Jones and James Selvage and Felix J. Herrmann}
}
@article {sharan2018fsp,
	title = {Fast sparsity-promoting microseismic source estimation},
	journal = {Geophysical Journal International},
	volume = {216},
	number = {1},
	year = {2019},
	note = {(Geophysical Journal International)},
	month = {01},
	pages = {164-181},
	abstract = {Microseismic events are generated during hydraulic fracturing of unconventional reservoirs and carry information on fracture locations and the origin times associated with these microseismic events. For drilling purposes and to prevent hazardous situations, we need to have accurate knowledge on the fracture locations as well as on their size, and their density. Because microseismic waves can travel far distances, microseismic data collected at the surface and or in boreholes can help us to monitor hydraulic fracturing. While so-called back propagation or time-reversal methods are able to focus recorded energy back onto the sources when a reasonable velocity model is available, these methods suffer from blurring especially in situations where the data acquisition suffers from lack of aperture, sparse sampling, and noise. As a result, these methods typically cannot resolve sources in close proximity, a desired feature since we need this information if we want to follow the fracture evolution in space and time. In that situation, we need to estimate the locations and the associated source-time functions for closely spaced microseismic sources along the active fractures. To overcome the limitations of time-reversal methods, we propose a wave-equation based inversion approach where we invert for the complete source wavefield in both space and time. By promoting sparsity on the source wavefield in space, we negate the effects of non-radiating sources during the inversion and obtain high-resolution intensity plots and high-fidelity estimates for the source-time functions. We obtain these results relatively fast by accelerating the linearized Bregman method with a dual formulation. Through experiments, we demonstrate that our method is computationally feasible, robust to noise, and works for closely spaced sources with overlapping source-time functions in complex geological settings.},
	keywords = {induced seismicity, Joint inversion, waveform inversion},
	doi = {10.1093/gji/ggy415},
	url = {https://slim.gatech.edu/Publications/Public/Journals/GeophysicalJournalInternational/2018/sharan2018fsp/sharan2018fsp.pdf},
	author = {Shashin Sharan and Rongrong Wang and Felix J. Herrmann}
}
@article {witte2019ecl,
	title = {Geophysics Bright Spots: Efficient coding of large-scale seismic inversion algorithms},
	journal = {The Leading Edge},
	volume = {38},
	number = {6},
	year = {2019},
	note = {(The Leading Edge)},
	pages = {482-484},
	abstract = {In {\textquotedblleft}A large-scale framework for symbolic implementations of seismic inversion algorithms in Julia,{\textquotedblright} Witte et al. describe new developments in how to code complex geophysical algorithms in a concise way. Subsurface seismic imaging and parameter estimation are among the most computationally challenging problems in the scientific community. Codes for solving seismic inverse problems, such as FWI or least-squares reverse time migration (LS-RTM), need to be highly optimized, but at the same time, facilitate the implementation of complex optimization algorithms. Traditionally, production-level codes in the oil and gas industry were exclusively written in low-level languages, such as C or Fortran, with extensive amounts of manual performance optimizations, thus making code maintenance, debugging, and adoption of new algorithms prohibitively challenging. Witte et al. present a paradigm of software engineering for seismic inverse problems based on symbolic user interfaces and code generation with automated performance optimization. Inspired by recent deep learning frameworks, the Julia Devito inversion framework (JUDI; an open-source software package) combines high-level abstractions for expressing seismic inversion algorithms with a domain-specific language compiler called Devito for solving the underlying wave equations. Devito{\textquoteright}s generated code is compiled just in time and outperforms codes with manual performance optimizations. JUDI utilizes Julia{\textquoteright}s high-level parallelization, making the software easily adaptable to a variety of computing environments such as densely connected HPC clusters or the cloud. The numerical examples (Figure 3) demonstrate the ability to implement a variety of complex algorithms for FWI and LS-RTM in a few lines of Julia code and run it on large-scale 3D models. The paper concludes that abstractions and performance are not mutually exclusive, and use of symbolic user interfaces can facilitate the implementation of new and innovative seismic inversion algorithms.},
	keywords = {bright spots, inversion, Julia, large-scale},
	doi = {10.1190/tle38060482.1},
	url = {https://library.seg.org/doi/10.1190/tle38060482.1},
	author = {Philipp A. Witte and Mathias Louboutin and Navjot Kukreja and Fabio Luporini and Michael Lange and Gerard J. Gorman and Felix J. Herrmann}
}
@article {witte2018alf,
	title = {A large-scale framework for symbolic implementations of seismic inversion algorithms in Julia},
	journal = {Geophysics},
	volume = {84},
	number = {3},
	year = {2019},
	note = {(Geophysics)},
	month = {03},
	pages = {F57-F71},
	abstract = {Writing software packages for seismic inversion is a very challenging task, since problems such as full-waveform inversion or least-squares imaging are both algorithmically and computationally demanding due to the large number of unknown parameters and the fact that we are propagating waves over many wavelengths. Software frameworks therefore need to combine both versatility and performance to provide geophysicists with the means and flexibility to implement complex algorithms that scale to exceedingly large 3D problems. Following these principles, we introduce the Julia Devito Inversion framework, an open-source software package in Julia for large-scale seismic modeling and inversion based on Devito, a domain-specific language compiler for automatic code generation. The framework consists of matrix-free linear operators for implementing seismic inversion algorithms that closely resembles the mathematical notation, a flexible resilient parallelization and an interface to Devito for generating optimized stencil code to solve the underlying wave equations. In comparison to many manually optimized industry codes written in low-level languages, our software is built on the idea of independent layers of abstractions and user interfaces with symbolic operators, making it possible to manage both the complexity of algorithms and performance optimizations, while preserving modularity, which allows for a level of expressiveness needed to formulate a broad range of wave-equation-based inversion problems. Through a series of numerical examples, we demonstrate that this allows users to implement algorithms for waveform inversion and imaging as simple Julia scripts that scale to large-scale 3D problems; thus providing a truly performant research and production framework.},
	keywords = {FWI, inversion, LSRTM, Modeling, software},
	doi = {10.1190/geo2018-0174.1},
	url = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2019/witte2018alf/witte2018alf.pdf},
	author = {Philipp A. Witte and Mathias Louboutin and Navjot Kukreja and Fabio Luporini and Michael Lange and Gerard J. Gorman and Felix J. Herrmann}
}
@conference {sharan2018SEGada,
	title = {A debiasing approach to microseismic},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2018},
	note = {(SEG, Anaheim)},
	month = {10},
	pages = {2942-2946},
	abstract = {Microseismic data is often used to locate fracture locations and their origin in time created by fracking. Although surface microseismic data can have large apertures and is easier to acquire than the borehole data, it often suffers from poor signal to noise ratio (S/R). Poor S/R poses a challenge in terms of estimating the correct location and source-time function of a microseismic source. In this work, we propose a denoising step in combination with a computationally cheap debiasing based approach to locate microseismic sources and to estimate their source-time functions with correct amplitude from extremely noisy data. Through numerical experiments, we demonstrate that our method can work with closely spaced microseismic sources with source-time functions of different peak amplitudes and frequencies. We have also shown the ability of our method with the smooth velocity model.},
	keywords = {Full-waveform inversion, Imaging, microseismic, noise, SEG, wave equation},
	doi = {10.1190/segam2018-2997252.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2018/sharan2018SEGada/sharan2018SEGada.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2018/sharan2018SEGada/sharan2018SEGada_pres.pdf},
	author = {Shashin Sharan and Rajiv Kumar and Rongrong Wang and Felix J. Herrmann}
}
@conference {louboutin2018SEGeow,
	title = {Effects of wrong adjoints for RTM in TTI media},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2018},
	note = {(SEG, Anaheim)},
	month = {10},
	pages = {331-335},
	abstract = {In order to obtain accurate images of the subsurface, anisotropic modeling and imaging is necessary. However, the twenty-one parameter complete wave-equation is too computationally expensive to be of use in this case. The transverse tilted isotropic wave-equation is then the next best feasible representation of the physics to use for imaging. The main complexity arising from transverse tilted isotropic imaging is to model the receiver wavefield (back propagation of the data or data residual) for the imaging condition. Unlike the isotropic or the full physics wave-equations, the transverse tilted isotropic wave-equation is not not self-adjoint. This difference means that time-reversal will not model the correct receiver wavefield and this can lead to incorrect subsurface images. In this work, we derive and implement the adjoint wave-equation to demonstrate the necessity of exact adjoint modeling for anisotropic modeling and compare our result with adjoint-free time-reversed imaging.},
	keywords = {Adjoint, anisotropy, finite-differences, Imaging, RTM, SEG, TTI},
	doi = {10.1190/segam2018-2996274.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2018/louboutin2018SEGeow/louboutin2018SEGeow.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2018/louboutin2018SEGeow/louboutin2018SEGeow_pres.pdf},
	author = {Mathias Louboutin and Philipp A. Witte and Felix J. Herrmann}
}
@article {louboutin2017fwip2,
	title = {Full-Waveform Inversion - Part 2: adjoint modeling},
	journal = {The Leading Edge},
	volume = {37},
	number = {1},
	year = {2018},
	note = {(The Leading Edge)},
	month = {1},
	pages = {69-72},
	abstract = {This tutorial is the second part of a three part tutorial series on full-waveform inversion (FWI), in which we provide a step by step walk through of setting up forward and adjoint wave equation solvers and an optimization framework for inversion. In part 1 (Louboutin et al., 2017), we demonstrated how to discretize the acoustic wave equation and how to set up a basic forward modeling scheme using Devito, a domain-specific language (DSL) in Python for automated finite-difference (FD) computations (Lange et al., 2016). Devito allows us to define wave equations as symbolic Python expressions (Meurer et al., 2017), from which optimized FD stencil code is automatically generated at run time. In part 1, we show how we can use Devito to set up and solve acoustic wave equations with (impulsive) seismic sources and sample wavefields at the receiver locations to model shot records.},
	keywords = {acoustic, devito, finite-difference, FWI, tutorial},
	doi = {10.1190/tle37010069.1},
	url = {https://slim.gatech.edu/Publications/Public/Journals/TheLeadingEdge/2018/louboutin2017fwip2/louboutin2017fwip2.html},
	author = {Mathias Louboutin and Philipp A. Witte and Michael Lange and Navjot Kukreja and Fabio Luporini and Gerard Gorman and Felix J. Herrmann}
}
@article {witte2018fwip3,
	title = {Full-Waveform Inversion - Part 3: optimization},
	journal = {The Leading Edge},
	volume = {37},
	number = {2},
	year = {2018},
	note = {(The Leading Edge)},
	month = {1},
	pages = {142-145},
	abstract = {This tutorial is the third part of a full-waveform inversion (FWI) tutorial series with a step-by-step walkthrough of setting up forward and adjoint wave equations and building a basic FWI inversion framework. For discretizing and solving wave equations, we use Devito, a Python-based domain-specific language for automated generation of finite-difference code (Lange et al., 2016). The first two parts of this tutorial (Louboutin et al., 2017, 2018) demonstrated how to solve the acoustic wave equation for modeling seismic shot records and how to compute the gradient of the FWI objective function using the adjoint-state method. With these two key ingredients, we will now build an inversion framework that can be used to minimize the FWI least-squares objective function.},
	keywords = {devito, finite-differences, FWI, inversion, Modeling, tutorial},
	doi = {10.1190/tle37020142.1},
	url = {https://slim.gatech.edu/Publications/Public/Journals/TheLeadingEdge/2018/witte2018fwip3/witte2018fwip3.html},
	author = {Philipp A. Witte and Mathias Louboutin and Keegan Lensink and Michael Lange and Navjot Kukreja and Fabio Luporini and Gerard Gorman and Felix J. Herrmann}
}
@conference {herrmann2018SMCtpa,
	title = {The power of abstraction in Computational Exploration Seismology},
	booktitle = {Smoky Mountains Computational Sciences and Engineering Conference},
	year = {2018},
	month = {08},
	abstract = {The field of Computational Exploration Seismology has over the years benefited tremendously from developments in HPC. Back in the 80{\textquoteright}s and early 90{\textquoteright}s, the oil \& gas industry was among the main drivers of HPC technology with a resurgence about ten years ago with the advent of full-waveform inversion-i.e., an adjoint-state method for inverse problems that involve the wave equation. While the combination of hardware developments and tedious manual coding efforts have resulted in major improvements, the rate of innovations is slow in part because of the inherent complexities of the mostly monolithic code bases. As a result, innovations make it too slowly into these optimized codes. This stifles innovation and could prevent adaptation of important developments such as machine learning with deep convolutional neural networks into existing workflows. In comparison, the technological advances in the academic as well as commercial machine learning communities have been much faster in part because of robust and well abstracted code bases such as PyTorch and TensorFlow. Since both machine learning and adjoint-state methods rely on back propagation strong similarities exist ready to be exploited. During my talk, I will discuss lessons we learned in finding the appropriate abstractions, read Domain Specific Language, for time-stepping stencil-based finite differences in Devito, and how this is relevant for machine learning with deep networks. I will also share a vision on how ideas from machine learning can be merged into the development of the next-generation of wave-equation based imaging codes. In particular, I will demonstrate the remarkable ability of deep convolutional networks to map low-fidelity dispersed numerical wave simulations to high-fidelity ones. I feel that this apparently generalizable capability of neural nets opens a complete new approach to tackle fundamental problems in computational exploration seismology.},
	keywords = {SMC},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SMC/2018/herrmann2018SMCtpa/herrmann2018SMCtpa_pres.pdf},
	author = {Felix J. Herrmann and Gerard J. Gorman and Jan H{\"u}ckelheim and Keegan Lensink and Paul H. J. Kelly and Navjot Kukreja and Henryk Modzelewski and Michael Lange and Mathias Louboutin and Fabio Luporini and Ali Siahkoohi and Philipp A. Witte}
}
@conference {yang2018SEGrde,
	title = {Removing density effects in LS-RTM with a low-rank matched filter},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2018},
	note = {(SEG, Anaheim)},
	month = {10},
	pages = {4251-4255},
	abstract = {Least-squares reverse-time migration faces difficulties when it inverts the data containing strong components related to density variation with velocity-only Born modeling operator. The strong density perturbation will be inverted as strong dummy velocity perturbations, which influence the amplitudes and phase of the velocity perturbations in the inverted model. The traditional method is to invert the additional density variations by developing Born operator with respect to both density and velocity or modify the image condition. In this work, we develop a matched-filter based LS-RTM for velocity-only Born modeling operator, which removes the artifacts in the imaging created by the strong density variation. This method doesn{\textquoteright}t call for extra work of finite difference stencil and is more general. In the experiment part, we use a complex discontinuous layered medium with strong density variations at the ocean bottom, and show the efficacy of the propose formulation.},
	keywords = {density effect, Least square, low rank, SEG},
	doi = {10.1190/segam2018-2997814.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2018/yang2018SEGrde/yang2018SEGrde.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2018/yang2018SEGrde/yang2018SEGrde_pres.pdf},
	author = {Mengmeng Yang and Rajiv Kumar and Rongrong Wang and Felix J. Herrmann}
}
@article {fang2017sewri,
	title = {Source estimation for wavefield-reconstruction inversion},
	journal = {Geophysics},
	volume = {83},
	number = {4},
	year = {2018},
	note = {(Geophysics)},
	pages = {R345-R359},
	abstract = {Source estimation is essential to all the wave-equation-based seismic inversions, including full-waveform inversion and the recently proposed wavefield-reconstruction inversion. When the source estimation is inaccurate, errors will propagate into the predicted data and introduce additional data misfit. As a consequence, inversion results that minimize this data misfit may become erroneous. To mitigate the errors introduced by the incorrect and pre-estimated sources, an embedded procedure that updates sources along with medium parameters is necessary for the inversion. So far, such a procedure is still missing in the context of wavefield-reconstruction inversion, a method that is, in many situations, less prone to local minima related to the so-called cycle skipping, compared to full-waveform inversion through exact data-fitting. While wavefield-reconstruction inversion indeed helps to mitigate issues related to cycle skipping by extending the search space with wavefields as auxiliary variables, it relies on having access to the correct source functions. In this paper, we remove the requirement of having the accurate source functions by proposing a source estimation technique specifically designed for wavefield-reconstruction inversion. To achieve this task, we consider the source functions as unknown variables and arrive at an objective function that depends on the medium parameters, wavefields, and source functions. During each iteration, we apply the so-called variable projection method to simultaneously project out the source functions and wavefields. After the projection, we obtain a reduced objective function that only depends on the medium parameters and invert for the unknown medium parameters by minimizing this reduced objective. Numerical experiments illustrate that this approach can produce accurate estimates of the unknown medium parameters without any prior information of the source functions.},
	keywords = {estimation, FWI, source, variable projection, WRI},
	doi = {10.1190/geo2017-0700.1},
	url = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2018/fang2017sewri/fang2017sewri.html},
	author = {Zhilong Fang and Rongrong Wang and Felix J. Herrmann}
}
@conference {sharan2018IEEEIUSspoa,
	title = {Sparsity-promoting photoacoustic imaging with source estimation},
	booktitle = {2018 IEEE International Ultrasonics Symposium (IUS)},
	year = {2018},
	note = {(IEEE IUS, Kobe)},
	month = {10},
	pages = {206-212},
	abstract = {Photoacoustics has emerged as a high-contrast imaging modality that provides optical absorption maps inside of tissues, therefore complementing morphological information of conventional ultrasound. The laser-generated photoacoustic waves are usually envelope-detected, thus disregarding the specific waveforms generated by each photoabsorber. Here we propose a sparsity-promoting image reconstruction method that allows the estimation of each photoabsorber{\textquoteright}s source-time function. Preliminary studies showed the ability to reconstruct the optical absorption map of an in silico vessel phantom. By using a sparsity-promoting imaging method, absorption maps and source-time functions can still be recovered even in situations where the number of transducers is decreased by a factor of six. Moreover, the recovery is able to attain higher resolution than conventional beamforming methods. Because our method recovers the source-time function of the absorbers, it could potentially also be used to distinguish different types of photoabsorbers, or the degree of aggregation of exogenous agents, under the assumption that these would generate different source-time functions at the moment of laser irradiation.},
	keywords = {Acceleration, estimation, IEEE, Image reconstruction, Phantoms, photoacoustic imaging, source-time function, sparsity promoting inversion, Transducers, wave equation},
	doi = {10.1109/ULTSYM.2018.8580037},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/IEEEIUS/2018/sharan2018IEEEIUSspoa/sharan2018IEEEIUSspoa.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/IEEEIUS/2018/sharan2018IEEEIUSspoa/sharan2018IEEEIUSspoa_pres.pdf},
	author = {Shashin Sharan and Rajiv Kumar and Diego S. Dumani and Mathias Louboutin and Rongrong Wang and Stanislav Emelianov and Felix J. Herrmann}
}
@techreport {witte2017TRcls,
	title = {Compressive least-squares migration with source estimation},
	number = {TR-EOAS-2017-3},
	year = {2017},
	institution = {UBC},
	abstract = {Least-squares reverse-time migration is a powerful approach for true-amplitude seismic imaging of complex geological structures. The successful application of this method is hindered by its exceedingly large computational cost and required prior knowledge of the generally unknown source wavelet. We address these problems by introducing an algorithm for low-cost sparsity-promoting least-squares migration with source estimation. We adapt a recent algorithm from sparse optimization, which allows to work with randomized subsets of shots during each iteration of least-squares migration, while still converging to an artifact-free solution. We modify the algorithm to incorporate on-the-fly source estimation through variable projection, which lets us estimate the wavelet without additional PDE solves. The algorithm is easy to implement and allows imaging at a fraction of the cost of conventional least squares reverse-time migration, requiring only around two passes trough the data, making the method feasible for large-scale imaging problems with unknown source wavelets.},
	keywords = {LSRTM, migration, source estimation, sparsity, time domain},
	url = {https://slim.gatech.edu/Publications/Public/TechReport/2017/witte2017TRcls/witte2017TRcls.html},
	author = {Philipp A. Witte and Mengmeng Yang and Felix J. Herrmann}
}
@conference {wang2017SEGdff,
	title = {A denoising formulation of full-waveform inversion},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2017},
	note = {(SEG, Houston)},
	month = {09},
	pages = {1594-1598},
	abstract = {We propose a wave-equation-based subsurface inversion method that in many cases is more robust than the conventional Full-Waveform Inversion. The new formulation is written in a denoising form that allows the synthetic data to match the observed ones up to a small error. Compared to the Full-Waveform Inversion, our method treats the noise arising from the data measuring/recording process and that from the synthetic modelling process separately. Comparing to the Wavefields Reconstruction Inversion, the new formulation mitigates the difficulty of choosing the penalty parameter λ. To solve the proposed optimization problem, we develop an efficient frequency domain algorithm that alternatively updates the model and the data. Numerical experiments confirm strong stability of the proposed method by comparisons between the results of our algorithm with that from both plain FWI and a weighted formulation of the FWI.},
	keywords = {denoising, FWI, SEG},
	doi = {10.1190/segam2017-17794690.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2017/wang2017SEGdff/wang2017SEGdff.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2017/wang2017SEGdff/wang2017SEGdff_pres.pdf},
	author = {Rongrong Wang and Felix J. Herrmann}
}
@article {louboutin2017fwi,
	title = {Full-Waveform Inversion - Part 1: forward modeling},
	journal = {The Leading Edge},
	volume = {36},
	number = {12},
	year = {2017},
	note = {(The Leading Edge)},
	month = {12},
	pages = {1033-1036},
	abstract = {Since its re-introduction by Pratt (1999), full-waveform inversion (FWI) has gained a lot of attention in geophysical exploration because of its ability to build high resolution velocity models more or less automatically in areas of complex geology. While there is an extensive and growing literature on the topic, publications focus mostly on technical aspects, making this topic inaccessible for a broader audience due to the lack of simple introductory resources for newcomers to geophysics. We will accomplish this by providing a hands-on walkthrough of FWI using Devito (Lange et al. 2016), a system based on domain-specific languages that automatically generates code for time-domainfinite-differences.},
	keywords = {devito, finite-differences, FWI, Modeling, tutorial},
	doi = {10.1190/tle36121033.1},
	url = {https://slim.gatech.edu/Publications/Public/Journals/TheLeadingEdge/2017/louboutin2017fwi/louboutin2017fwi.html},
	author = {Mathias Louboutin and Philipp A. Witte and Michael Lange and Navjot Kukreja and Fabio Luporini and Gerard Gorman and Felix J. Herrmann}
}
@conference {sharan2017SEGhrf,
	title = {High resolution fast microseismic source collocation and source time function estimation},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2017},
	note = {(SEG, Houston)},
	month = {09},
	pages = {2778-2783},
	abstract = {Sparsity promotion based joint microseismic source collocation and source time function estimation, using Linearized Bregman algorithm, although simple to implement, suffers from slow convergence. This is due to the fact that Linearized Bregman algorithm has only first order of convergence. In this work, we propose to accelerate the existing Linearized Bregman algorithm using the L-BFGS algorithm. Without any initial guess for the source location or source time function, our method is able to estimate the source location and source time function for kinematically correct smooth velocity model. We demonstrate the effectiveness of our method for multiple sources spaced within half a wavelength. We also compare our results with Linearized Bregman based method in {\textquoteleft}{\textquoteleft}2.5{\textquoteleft}{\textquoteleft}D instead of {\textquoteleft}{\textquoteleft}2{\textquoteleft}{\textquoteleft}D.},
	keywords = {algorithm, Full-waveform inversion, microseismic, passive imaging, SEG, wave equation},
	doi = {10.1190/segam2017-17787749.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2017/sharan2017SEGhrf/sharan2017SEGhrf.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2017/sharan2017SEGhrf/sharan2017SEGhrf_pres.pdf},
	author = {Shashin Sharan and Rongrong Wang and Felix J. Herrmann}
}
@article {kumar2017hrc,
	title = {Highly repeatable 3D compressive full-azimuth towed-streamer time-lapse acquisition {\textendash}- a numerical feasibility study at scale},
	journal = {The Leading Edge},
	volume = {36},
	number = {8},
	year = {2017},
	note = {(The Leading Edge)},
	month = {08},
	pages = {677-687},
	abstract = {Most conventional 3D time-lapse (or 4D) acquisitions are ocean-bottom cable (OBC) or ocean-bottom node (OBN) surveys since these surveys are relatively easy to replicate compared to towed-streamer surveys. To attain high degrees of repeatability, survey replicability and dense periodic sampling has become the norm for 4D surveys that renders this technology expensive. Conventional towed-streamer acquisitions suffer from limited illumination of subsurface due to narrow azimuth. Although, acquisition techniques such as multi-azimuth, wide-azimuth, rich-azimuth acquisition, etc., have been developed to illuminate the subsurface from all possible angles, these techniques can be prohibitively expensive for densely sampled surveys. This leads to uneven sampling, i.e., dense receiver and coarse source sampling or vice-versa, in order to make these acquisitions more affordable. Motivated by the design principles of Compressive Sensing (CS), we acquire economic, randomly subsampled (or compressive) and simultaneous towed-streamer time-lapse data without the need of replicating the surveys. We recover densely sampled time-lapse data on one and the same periodic grid by using a joint-recovery model (JRM) that exploits shared information among different time-lapse recordings, coupled with a computationally cheap and scalable rank-minimization technique. The acquisition is low cost since we have subsampled measurements (about 70\% subsampled), simulated with a simultaneous long-offset acquisition configuration of two source vessels travelling across a survey area at random azimuths. We analyze the performance of our proposed compressive acquisition and subsequent recovery strategy by conducting a synthetic, at scale, seismic experiment on a 3D time-lapse model containing geological features such as channel systems, dipping and faulted beds, unconformities and a gas cloud. Our findings indicate that the insistence on replicability between surveys and the need for OBC/OBN 4D surveys can, perhaps, be relaxed. Moreover, this is a natural next step beyond the successful CS acquisition examples discussed in this special issue.},
	keywords = {3D, CS, marine, rank minimization, simultaneous long offset, time-lapse seismic},
	doi = {10.1190/tle36080677.1},
	url = {https://slim.gatech.edu/Publications/Public/Journals/TheLeadingEdge/2017/kumar2017hrc/kumar2017hrc.html},
	author = {Rajiv Kumar and Haneet Wason and Shashin Sharan and Felix J. Herrmann}
}
@presentation {herrmann2017SINBADFhrc,
	title = {Highly repeatable 3D compressive full-azimuth towed-streamer time-lapse acquisition {\textendash}- a numerical feasibility study at scale},
	year = {2017},
	publisher = {SINBAD},
	abstract = {Most conventional 3D time-lapse (or 4D) acquisitions are ocean-bottom cable (OBC) or ocean-bottom node (OBN) surveys since these surveys are relatively easy to replicate compared to towed-streamer surveys. To attain high degrees of repeatability, survey replicability and dense periodic sampling has become the norm for 4D surveys that renders this technology expensive. Conventional towed-streamer acquisitions suffer from limited illumination of subsurface due to narrow azimuth. Although, acquisition techniques such as multi-azimuth, wide-azimuth, rich-azimuth acquisition, etc., have been developed to illuminate the subsurface from all possible angles, these techniques can be prohibitively expensive for densely sampled surveys. This leads to uneven sampling, i.e., dense receiver and coarse source sampling or vice-versa, in order to make these acquisitions more affordable. Motivated by the design principles of Compressive Sensing (CS), we acquire economic, randomly subsampled (or compressive) and simultaneous towed-streamer time-lapse data without the need of replicating the surveys. We recover densely sampled time-lapse data on one and the same periodic grid by using a joint-recovery model (JRM) that exploits shared information among different time-lapse recordings, coupled with a computationally cheap and scalable rank-minimization technique. The acquisition is low cost since we have subsampled measurements (about 70\% subsampled), simulated with a simultaneous long-offset acquisition configuration of two source vessels travelling across a survey area at random azimuths. We analyze the performance of our proposed compressive acquisition and subsequent recovery strategy by conducting a synthetic, at scale, seismic experiment on a 3D time-lapse model containing geological features such as channel systems, dipping and faulted beds, unconformities and a gas cloud. Our findings indicate that the insistence on replicability between surveys and the need for OBC/OBN 4D surveys can, perhaps, be relaxed. Moreover, this is a natural next step beyond the successful CS acquisition examples discussed during this session.},
	keywords = {Presentation, SINBAD, SINBADFALL2017, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/herrmann2017SINBADFhrc/herrmann2017SINBADFhrc.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/herrmann2017SINBADFhrc/herrmann2017SINBADFhrc.mov},
	author = {Felix J. Herrmann and Rajiv Kumar and Haneet Wason and Shashin Sharan and Felix Oghenekohwo}
}
@presentation {witte2017SINBADFals,
	title = {A large-scale framework in Julia for fast prototyping of seismic inversion algorithms},
	year = {2017},
	publisher = {SINBAD},
	abstract = {We present our progress on a large-scale seismic modeling workflow in Julia for wave-equation based inversion. The software offers a range of high-level abstractions to easily express PDE constrained optimization problems in terms of linear algebra expressions, while utilizing the DSL Devito to symbolically express the underlying PDEs and to generate fast and parallel code for solving them. Data containers and linear operators can be set up without much effort from input SEG-Y data and scale to large-scale 3D applications. This talk provides an overview of the basic functionalities of our software and applications to least squares imaging and 3D FWI.},
	keywords = {Presentation, SINBAD, SINBADFALL2017, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/witte2017SINBADFals/witte2017SINBADFals.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/witte2017SINBADFals/witte2017SINBADFals.mov},
	author = {Philipp A. Witte and Mathias Louboutin and Felix J. Herrmann}
}
@techreport {witte2016OGHPClst,
	title = {A large-scale time-domain seismic modeling and inversion workflow in Julia},
	number = {TR-EOAS-2017-1},
	year = {2017},
	institution = {UBC},
	abstract = {We present our initial steps towards the development of a large-scale seismic modeling workflow in Julia that provides a framework for wave equation based inversion methods like full waveform inversion or least squares migration. Our framework is based on the Devito, a finite difference domain specific language compiler that generates highly optimized and parallel code. We develop a flexible workflow that is based on abstract matrixfree linear operators and enables developers to write code that closely resembles the underlying math, while at the same time leveraging highly optimized wave equation solvers, allowing us to solve large-scale three-dimensional inverse problems.},
	keywords = {HPC, inversion, Modelling},
	url = {https://slim.gatech.edu/Publications/Public/TechReport/2017/witte2016OGHPClst/witte2016OGHPClst.pdf},
	author = {Philipp A. Witte and Mathias Louboutin and Gerard Gorman and Felix J. Herrmann}
}
@conference {witte2017SIAMlsw,
	title = {Large-scale workflows for wave-equation based inversion in Julia},
	booktitle = {SIAM Conference on Computational Science and Engineering},
	year = {2017},
	note = {(SIAM Conference on Computational Science and Engineering, Atlanta)},
	month = {02-03},
	abstract = {Full-waveform inversion is a PDE-constrained optimisation problem involving massive amounts of data (petabytes) and large numbers of unknowns (O(10^9)). This well known compute-intensive and data-intensive is extremely challenging for several reasons. First, there is the complexity of having to handle extremely large data volumes with metadata related to experimental details in the field, and the discretization of the unknown earth parameters and approximate physics. Second, reduced or adjoint-state methods call for computationally intensive PDE solves for each source experiment (of which there are thousands) for each iteration of a gradient-based optimization scheme. The talks will give an overview how carefully chosen layers of abstraction can help manage both the complexity and scale of inversion while still achieving the high degree of computational performance required to make full-waveform a practical tool. Specifically, the presentations will focus on domain specific stencil language for time-stepping methods to solve various types of wave equations and on abstracts for large-scale parallel optimization frameworks.},
	keywords = {Full-waveform inversion, Julia, large scale, least-squares migration, SIAM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SIAM/2017/witte2017SIAMlsw/witte2017SIAMlsw_pres.mov},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SIAM/2017/witte2017SIAMlsw/witte2017SIAMlsw_pres.pdf},
	author = {Philipp A. Witte and Mathias Louboutin and Felix J. Herrmann}
}
@presentation {louboutin2017SINBADFldi,
	title = {Latest developments in Devito},
	year = {2017},
	publisher = {SINBAD},
	abstract = {We present an overview of the latest developments in Devito. We introduced Devito in the previous meeting as a prototype finite-difference DSL for seismic modelling and inversion. We are presenting here the latest improvements and functionalities of Devito. We will also discuss the current future plans as well as non-supported features that the audience may be interested in. This presentation will be followed by/mixed with a hands-in tutorial if the time and resources allows it.},
	keywords = {Presentation, SINBAD, SINBADFALL2017, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/louboutin2017SINBADFldi/louboutin2017SINBADFldi.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/louboutin2017SINBADFldi/louboutin2017SINBADFldi.mov},
	author = {Mathias Louboutin and Michael Lange and Fabio Luporini and Navjot Kurjeka and Jan Hueckelheim and Gerard Gorman and Philipp A. Witte and Felix J. Herrmann}
}
@article {oghenekohwo2016GEOPctl,
	title = {Low-cost time-lapse seismic with distributed compressive sensing{\textendash}-Part 1: exploiting common information among the vintages},
	journal = {Geophysics},
	volume = {82},
	number = {3},
	year = {2017},
	note = {(Geophysics)},
	month = {05},
	pages = {P1-P13},
	abstract = {Time-lapse seismic is a powerful technology for monitoring a variety of subsurface changes due to reservoir fluid flow. However, the practice can be technically challenging when one seeks to acquire colocated time-lapse surveys with high degrees of replicability among the shot locations. We have determined that under "ideal" circumstances, in which we ignore errors related to taking measurements off the grid, high-quality prestack data can be obtained from randomized subsampled measurements that are observed from surveys in which we choose not to revisit the same randomly subsampled on-the-grid shot locations. Our acquisition is low cost because our measurements are subsampled. We have found that the recovered finely sampled prestack baseline and monitor data actually improve significantly when the same on-the-grid shot locations are not revisited. We achieve this result by using the fact that different time-lapse data share information and that nonreplicated (on-the-grid) acquisitions can add information when prestack data are recovered jointly. Whenever the time-lapse data exhibit joint structure{\textendash}-i.e., they are compressible in some transform domain and share information{\textendash}-sparsity-promoting recovery of the "common component" and "innovations," with respect to this common component, outperforms independent recovery of the prestack baseline and monitor data. The recovered time-lapse data are of high enough quality to serve as the input to extract poststack attributes used to compute time-lapse differences. Without joint recovery, artifacts{\textendash}-due to the randomized subsampling{\textendash}-lead to deterioration of the degree of repeatability of the time-lapse data. We tested this method by carrying out experiments with reliable statistics from thousands of repeated experiments. We also confirmed that high degrees of repeatability are achievable for an ocean-bottom cable survey acquired with time-jittered continuous recording.},
	keywords = {Acquisition, joint recovery method, marine, random sampling, time-lapse seismic},
	doi = {10.1190/geo2016-0076.1},
	url = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2017/oghenekohwo2016GEOPctl/oghenekohwo2016GEOPctl.html},
	author = {Felix Oghenekohwo and Haneet Wason and Ernie Esser and Felix J. Herrmann}
}
@article {wason2016GEOPctl,
	title = {Low-cost time-lapse seismic with distributed compressive sensing{\textendash}-Part 2: impact on repeatability},
	journal = {Geophysics},
	volume = {82},
	number = {3},
	year = {2017},
	note = {(Geophysics)},
	month = {05},
	pages = {P15-P30},
	abstract = {Irregular or off-the-grid spatial sampling of sources and receivers is inevitable in field seismic acquisitions. Consequently, time-lapse surveys become particularly expensive since current practices aim to replicate densely sampled surveys for monitoring changes occurring in the reservoir due to hydrocarbon production. We demonstrate that under certain circumstances, high-quality prestack data can be obtained from cheap randomized subsampled measurements that are observed from nonreplicated surveys. We extend our time-jittered marine acquisition to time-lapse surveys by designing acquisition on irregular spatial grids that render simultaneous, subsampled and irregular measurements. Using the fact that different time-lapse data share information and that nonreplicated surveys add information when prestack data are recovered jointly, we recover periodic densely sampled and colocated prestack data by adapting the recovery method to incorporate a regularization operator that maps traces from an irregular spatial grid to a regular periodic grid. The recovery method is, therefore, a combined operation of regularization, interpolation (estimating missing fine-grid traces from subsampled coarse-grid data), and source separation (unraveling overlapping shot records). By relaxing the insistence on replicability between surveys, we find that recovery of the time-lapse difference shows little variability for realistic field scenarios of slightly nonreplicated surveys that suffer from unavoidable natural deviations in spatial sampling of shots (or receivers) and pragmatic compressed-sensing based nonreplicated surveys when compared to the "ideal" scenario of exact replicability between surveys. Moreover, the recovered densely sampled prestack baseline and monitor data improve significantly when the acquisitions are not replicated, and hence can serve as input to extract poststack attributes used to compute time-lapse differences. Our observations are based on experiments conducted for an ocean-bottom cable survey acquired with time-jittered continuous recording assuming source equalization (or same source signature) for the time-lapse surveys and no changes in wave heights, water column velocities or temperature and salinity profiles, etc.},
	keywords = {joint recovery method, Marine acquisition, off-the-grid recovery, Optimization, random sampling, time-lapse seismic},
	doi = {10.1190/geo2016-0252.1},
	url = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2017/wason2016GEOPctl/wason2016GEOPctl.html},
	author = {Haneet Wason and Felix Oghenekohwo and Felix J. Herrmann}
}
@presentation {wang2017SINBADFnra,
	title = {Noise robust and time-domain formulations of Wavefield Reconstruction Inversion},
	year = {2017},
	publisher = {SINBAD},
	abstract = {We propose a wave-equation-based subsurface inversion method that in many cases is more robust than conventional Full-Waveform Inversion. The new formulation is written in a denoising form that allows the synthetic data to match the observed ones up to a small error. Compared to regular Full-Waveform Inversion, our method treats the noise arising from the data meassuring/recording process and that from the synthetic modelling process separately. Compared to Wavefields Reconstruction Inversion, the new formulation mitigates the difficulty of choosing the penalty parameter λ. To solve the proposed optimization problem, we develop an efficient frequency domain algorithm that alternatively updates the model and the data. Numerical experiments confirm strong stability of the proposed method by comparisons between the results of our algorithm with that from both plain FWI and a weighted formulation of the FWI. We also discuss a new memory efficient time-domain formulation for Wavefield Reconstruction Inversion based on duality.},
	keywords = {Presentation, SINBAD, SINBADFALL2017, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/wang2017SINBADFnra/wang2017SINBADFnra.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/wang2017SINBADFnra/wang2017SINBADFnra.mov},
	author = {Rongrong Wang and Mathias Louboutin and Bas Peters and Emmanouil Daskalakis and Felix J. Herrmann}
}
@mastersthesis {wason2017THsss,
	title = {Simultaneous-source seismic data acquisition and processing with compressive sensing},
	year = {2017},
	note = {(PhD)},
	month = {08},
	school = {The University of British Columbia},
	type = {phd},
	address = {Vancouver},
	abstract = {The work in this thesis adapts ideas from the field of compressive sensing (CS) that lead to new insights into acquiring and processing seismic data, where we can fundamentally rethink how we design seismic acquisition surveys and process acquired data to minimize acquisition- and processing-related costs. Current efforts towards dense source/receiver sampling and full azimuthal coverage to produce high-resolution images of the subsurface have led to the deployment of multiple sources across survey areas. A step ahead from multisource acquisition is simultaneous-source acquisition, where multiple sources fire shots at near-simultaneous/random times resulting in overlapping shot records, in comparison to no overlaps during conventional sequential-source acquisition. Adoption of simultaneous-source techniques has helped to improve survey efficiency and data density. The engine that drives simultaneous-source technology is simultaneous-source separation {\textendash}- a methodology that aims to recover conventional sequential-source data from simultaneous-source data. This is essential because many seismic processing techniques rely on dense and periodic (or regular) source/receiver sampling. We address the challenge of source separation through a combination of tailored simultaneous-source acquisition design and sparsity-promoting recovery via convex optimization using l1 objectives. We use CS metrics to investigate the relationship between marine simultaneous-source acquisition design and data reconstruction fidelity, and consequently assert the importance of randomness in the acquisition system in combination with an appropriate choice for a sparsifying transform (i.e., curvelet transform) in the reconstruction algorithm. We also address the challenge of minimizing the cost of expensive, dense, periodically-sampled and replicated time-lapse surveying and data processing by adapting ideas from distributed compressive sensing. We show that compressive randomized time-lapse surveys need not be replicated to attain acceptable levels of data repeatability, as long as we know the shot positions (post acquisition) to a sufficient degree of accuracy. We conclude by comparing sparsity-promoting and rank-minimization recovery techniques for marine simultaneous-source separation, and demonstrate that recoveries are comparable; however, the latter approach readily scales to large-scale seismic data and is computationally faster.},
	keywords = {Acquisition, Compressive Sensing, marine, Optimization, PhD, simultaneous source, source separation},
	url = {https://slim.gatech.edu/Publications/Public/Thesis/2017/wason2017THsss/wason2017THsss.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Thesis/2017/wason2017THsss/wason2017THsss_pres.pdf},
	author = {Haneet Wason}
}
@conference {witte2017EAGEspl,
	title = {Sparsity-promoting least-squares migration with the linearized inverse scattering imaging condition},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2017},
	note = {(EAGE, Paris)},
	month = {06},
	abstract = {Reverse-time migration (RTM) with the conventional cross-correlation imaging condition suffers from low-frequency artifacts that result from backscattered energy in the background velocity models. This problem translates to least-squares reverse-time migration (LS-RTM), where these artifacts slow down the convergence, as many of the initial iterations are spent on removing them. In RTM, this problem has been successfully addressed by the introduction of the so-called inverse scattering imaging condition, which naturally removes these artifacts. In this work, we derive the corresponding linearized forward operator of the inverse scattering imaging operator and incorporate this forward/adjoint operator pair into a sparsity-promoting (SPLS-RTM) workflow. We demonstrate on a challenging salt model, that LS-RTM with the inverse scattering imaging condition is far less prone to low-frequency artifacts than the conventional cross-correlation imaging condition, improves the convergence and does not require any type of additional image filters within the inversion. Through source subsampling and sparsity promotion, we reduce the computational cost in terms of PDE solves to a number comparable to conventional RTM, making our workflow applicable to large-scale problems.},
	keywords = {EAGE, imaging condition, least-squares migration, linearized Bregman},
	doi = {10.3997/2214-4609.201701125},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2017/witte2017EAGEspl/witte2017EAGEspl.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2017/witte2017EAGEspl/witte2017EAGEspl_poster.pdf},
	author = {Philipp A. Witte and Mengmeng Yang and Felix J. Herrmann}
}
@presentation {sharan2017SINBADFtts,
	title = {Tracking the spatial-temporal evolution of fractures by microseismic source collocation},
	year = {2017},
	publisher = {SINBAD},
	abstract = {Unlike conventional reservoirs, unconventional plays are not naturally viable for economical production of oil and gas. They require stimulation by injecting high-pressure fluid causing fractures in the rocks. These fractures make the medium more permeable, hence, the extraction of oil and gas becomes feasible. For drilling purposes and to prevent potentially hazardous situations, we need to have good knowledge of the location of these fractures. Also, we need to have good knowledge about how these fractures originated in time. Hydraulic fracturing changes stress in rocks, which results in the emission of microseismic waves. The opening of cracks due to high pressure fluid injection during hydraulic fracturing mainly causes this change in stress in the rocks. Therefore, microseismic events are mostly localized along these fractures and have finite energy along time. To accurately track the evolution of fractures in both space and time, we need to locate closely spaced microseismic events along these fractures activating at very small time intervals. A naive approach can be the back propagation of the observed data to find out a point in space and time where maximum focusing of back propagating energy occurs. This point corresponds to the location and origin time of a microseismic source. This approach, although simpler, suffers from low resolution and requires scanning of complete 4D volume (3D in space and 1D in time). Hence, this method can be challenging when there are multiple closely spaced microseismic sources originating at different times. We in this work propose a sparsity promotion based method that can locate closely spaced microseismic events, with spatial separation as low as within half a wavelength, activating at small time intervals. We simultaneously estimate the origin time of microseismic events by estimating their source time functions. Our method exploits the fact that microseismic events are localized in space and have finite energy. We use accelerated Linearized Bregman algorithm with a preconditioning operator to arrive at a computationally feasible scheme.},
	keywords = {Presentation, SINBAD, SINBADFALL2017, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/sharan2017SINBADFtts/sharan2017SINBADFtts.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/sharan2017SINBADFtts/sharan2017SINBADFtts.mov},
	author = {Shashin Sharan and Rongrong Wang and Felix J. Herrmann}
}
@presentation {oghenekohwo2016SINBADFast,
	title = {Asymmetric sampling for time-lapse surveys},
	year = {2016},
	publisher = {SINBAD},
	abstract = {In this talk, we use the joint recovery model (JRM), which derives from distributed compressed sensing, to reconstruct time-lapse data acquired via time-jittered marine sources that are significantly subsampled in the monitor survey(s) in comparison to the baseline survey, hence, saving acquisition cost. We analyze the effects of two or more asymmetrically subsampled monitor surveys on wavefield reconstruction via JRM. We leverage the common information shared amongst the vintages during our sparsity-promoting recovery.},
	keywords = {Presentation, SINBAD, SINBADFALL2016, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/oghenekohwo2016SINBADFast/oghenekohwo2016SINBADFast.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/oghenekohwo2016SINBADFast/oghenekohwo2016SINBADFast.mov},
	author = {Felix Oghenekohwo and Haneet Wason and Felix J. Herrmann}
}
@conference {herrmann2016SEGWSctl,
	title = {Compressive time-lapse marine acquisition},
	booktitle = {SEG Workshop on Low cost geophysics: How to be creative in a cost-challenged environment; Dallas},
	year = {2016},
	note = {(SEG Workshop, Dallas)},
	month = {10},
	keywords = {Acquisition, SEG, workshop},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/herrmann2016SEGWSctl/herrmann2016SEGWSctl_pres.pdf},
	author = {Felix J. Herrmann and Rajiv Kumar and Felix Oghenekohwo and Shashin Sharan and Haneet Wason}
}
@article {esser2016TLEcwi,
	title = {Constrained waveform inversion for automatic salt flooding},
	journal = {The Leading Edge},
	volume = {35},
	number = {3},
	year = {2016},
	note = {(The Leading Edge)},
	month = {03},
	pages = {235-239},
	abstract = {Given appropriate data acquisition, processing to remove nonprimary arrivals, and use of an accurate migration algorithm, it is the quality of the subsurface velocity model that typically controls the quality of imaging that can be obtained from salt-affected seismic data. Full-waveform inversion has the potential to improve the accuracy, resolution, repeatability, and speed with which such velocity models can be generated, but, in the absence of an accurate starting model, that potential is difficult to realize in practice. Presented are successful inversion results, obtained from synthetic subsalt models, using a robust full-waveform inversion code that includes constraints upon the set of allowable earth models. These constraints include limitations on the total variation of the velocity of the model and, most significantly, on the asymmetric variation of velocity with depth such that negative velocity excursions are limited. During the iteration, these constraints are relaxed progressively so that the final model is driven principally by the seismic data, but the constraints act to steer the inversion path away from local minima in its early stages. This methodology is applied to portions of the 2004 BP benchmark and Phase I SEAM salt models, recovering an accurate model of the salt body, including its base and flanks, and an accurate model of the subsalt velocity structure, starting from one-dimensional velocity models that are severely cycle skipped. This approach removes entirely the requirement to pick salt boundaries from migrated seismic data, and acts as a form of automatic salt and sediment flooding during full-waveform inversion.},
	keywords = {FWI, high-contrast and high velocity, salt, total-variation minimization, velocity-model building},
	doi = {10.1190/tle35030235.1},
	url = {http://dx.doi.org/10.1190/tle35030235.1},
	author = {Ernie Esser and Llu{\'\i}s Guasch and Felix J. Herrmann and Mike Warner}
}
@presentation {kumar2016SINBADFels,
	title = {Efficient large-scale 5D seismic data acquisition and processing using rank-minimization},
	year = {2016},
	publisher = {SINBAD},
	abstract = {Seismic data collection is becoming challenging because of increased demands for high-quality, long-offset and wide-azimuth data. Leveraging ideas from CS, in this work we establish a cost effective acquisition and processing techniques, which are no longer dominated by survey area size but by the sparsity of seismic data volumes. In the first part of abstract, we establish connections between random time dithering and jittered sampling in space. Specifically, we recover high-quality 5D seismic data volumes from time-jittered marine acquisition where the average inter-shot time is significantly reduced, leading to cheaper surveys due to fewer overlapping shots. The time-jittered acquisition, in conjunction with the shot separation by Singular-Value Decomposition (SVD)-free factorization based rank-minimization approach, allows us to recover high quality 5D seismic data volumes. Results are illustrated for simulations of simultaneous time-jittered continuous recording for a 3D ocean-bottom cable survey, where we outperforms existing techniques, by an order of magnitude computational speedup and using 1/20th of the memory, that use sparsity in transforms domains. The second part of abstract focussed on leveraging low-rank structure in seismic data to solve extremely large data recovery (interpolation) problems. We introduced a large{\textendash}scale SVD-free optimization framework that is robust with respect to outliers and that uses information on the support. We test the efficacy of the proposed interpolation framework on a large-scale 5D seismic data, generated from the geologically complex synthetic 3D Compass velocity model, where 80\% of the data has been removed. Our findings show that major computational and memory gains are possible compared to curvelet-based reconstruction.},
	keywords = {Presentation, SINBAD, SINBADFALL2016, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/kumar2016SINBADFels/kumar2016SINBADFels.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/kumar2016SINBADFels/kumar2016SINBADFels.mov},
	author = {Rajiv Kumar and Shashin Sharan and Haneet Wason and Felix J. Herrmann}
}
@conference {wang2016SEGfde,
	title = {Frequency down-extrapolation with TV norm minimization},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2016},
	note = {(SEG, Dallas)},
	month = {10},
	pages = {1380-1384},
	abstract = {In this work, we present a total-variation (TV)-norm minimization method to perform model-free low-frequency data extrapolation for the purpose of assisting full-waveform inversion. Low-frequency extrapolation is the process of extending reliable frequency bands of the raw data towards the lower end of the spectrum. To this end, we propose to solve a TV-norm based convex optimization problem that has a global minimum and is equipped with a fast solver. The approach takes into account both the sparsity of the reflectivity series associated with a single trace, as well as the inter-trace correlations. A favorable byproduct of this approach is that it allows one to work with coarsely sampled trace along time (as coarse as 0.02s per sample), hence substantially reduces the size of the proposed optimization problem. Last, we show the effectiveness of the method for frequency-domain FWI on the Marmousi model.},
	keywords = {extrapolation, FWI, SEG, sparse},
	doi = {10.1190/segam2016-13879674.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/wang2016SEGfde/wang2016SEGfde.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/wang2016SEGfde/wang2016SEGfde_pres.pdf},
	author = {Rongrong Wang and Felix J. Herrmann}
}
@presentation {sharan2016SINBADFhrm,
	title = {High resolution microseismic source collocation},
	year = {2016},
	publisher = {SINBAD},
	abstract = {Microseismic waves generated by hydraulic fracturing are used to locate the fractures created in the reservoir. Most of the existing source collocation methods have their limitations in terms of resolution of very closely spaced microseismic events. In this work we propose a method to collocate microseismic events with high resolution along with simultaneous estimation of source time function. In the context of waves, the finite-difference modelling operator acts as an analysis operator for wavefields. This implies that wavefields under the action of a finite-difference modelling kernel can be focused to sources generating these wavefields. We exploit this property of wavefields to locate microseismic events and estimate the source time function jointly. To arrive at a computational feasible scheme, we use the accelerated version of Linearized Bregman (LBR) algorithm to solve the above mentioned problem. This is joint work with Rongrong Wang.},
	keywords = {Presentation, SINBAD, SINBADFALL2016, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/sharan2016SINBADFhrm/sharan2016SINBADFhrm.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/sharan2016SINBADFhrm/sharan2016SINBADFhrm.mov},
	author = {Shashin Sharan and Rongrong Wang and Felix J. Herrmann}
}
@presentation {witte2016SINBADFlst,
	title = {A large-scale time-domain modeling and inversion workflow in Julia},
	year = {2016},
	publisher = {SINBAD},
	abstract = {We present our initial steps towards the development of a large-scale seismic modeling workflow in Julia which provides a framework for wave-equation based inversion methods like full waveform inversion or least squares migration. Our framework is based on the Devito tool for optimized finite difference computations that generates highly optimized code from symbolic PDEs. We aim at developing a flexible workflow which is based on abstract operators and allows coding that is {\textquotedblleft}close to the math{\textquotedblright}, while at the same time relying on extremely fast and efficient, state of the art wave equation solvers.},
	keywords = {Presentation, SINBAD, SINBADFALL2016, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/witte2016SINBADFlst/witte2016SINBADFlst.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/witte2016SINBADFlst/witte2016SINBADFlst.mov},
	author = {Philipp A. Witte and Felix J. Herrmann}
}
@conference {chai2016SEGlbm,
	title = {A linearized Bregman method for compressive waveform inversion},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2016},
	note = {(SEG, Dallas)},
	month = {10},
	pages = {1449-1454},
	abstract = {We present an implementation of a recently developed and relatively simple linearized Bregman method to solve the large-scale $\ell_1$-norm sparsity-promoting Gauss-Newton (GN) full-waveform inversion (FWI) problem. Numerical experiments demonstrate that: the much simpler linearized Bregman method does a comparable and even superior job compared to a state-of-the-art $\ell_1$-norm solver SPGL1, which is previously used in the modified Gauss-Newton FWI; the linearized Bregman method is also more efficient (faster) than SPGL1; the FWI result with the linearized Bregman method solving $\ell_1$-norm sparsity-promoting problems to get the model updates is better than that obtained by solving $\ell_2$-norm least-squares problems to get the model updates.},
	keywords = {Bregman method, SEG, waveform inversion},
	doi = {10.1190/segam2016-13848105.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/chai2016SEGlbm/chai2016SEGlbm.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/chai2016SEGlbm/chai2016SEGlbm_pres.pdf},
	author = {Xintao Chai and Mengmeng Yang and Philipp A. Witte and Rongrong Wang and Zhilong Fang and Felix J. Herrmann}
}
@presentation {wason2016SINBADFlcr,
	title = {Low-cost, randomized 3D towed-marine time-lapse seismic acquisition},
	year = {2016},
	publisher = {SINBAD},
	abstract = {Most often conventional 3D towed-marine seismic data acquisition suffers from limited illumination of subsurface due to narrow azimuth. Although, acquisition techniques such as Ocean Bottom Cable (OBC), Ocean Bottom Node (OBN), Rich-Azimuth (RAZ) acquisition, Multi-Azimuth (MAZ) acquisition have been developed to illuminate the subsurface from all possible angles, these acquisition techniques are expensive. Motivated by the design principles of compressed sensing, we acquire randomly subsampled and simultaneous (e.g., simultaneous long offset) 3D towed-marine time-lapse seismic data for few (random) azimuths, and recover densely sampled multiple azimuth baseline and monitor data via sparsity promotion and rank minimization. Our acquisition is low cost since we have subsampled measurements.},
	keywords = {Presentation, SINBAD, SINBADFALL2016, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/wason2016SINBADFlcr/wason2016SINBADFlcr.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/wason2016SINBADFlcr/marine_acq_mask.pdf},
	author = {Rajiv Kumar and Felix Oghenekohwo and Shashin Sharan and Haneet Wason and Felix J. Herrmann}
}
@conference {witte2016SEGpve,
	title = {Phase velocity error minimizing scheme for the anisotropic pure P-wave equation},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2016},
	note = {(SEG, Dallas)},
	month = {10},
	pages = {452-457},
	abstract = {Pure P-wave equations for acoustic modeling in transverse isotropic media are derived by approximating the exact pure Pwave dispersion relation. In this work, we present an alternative approach to the approximate dispersion relation of Etgen and Brandsberg-Dahl, in which we approximate the exact dispersion relation through a polynomial expansion and determine its coefficients by solving a linear least squares problem that minimizes the phase velocity error over the entire range of phase angles. The coefficients are also optimized over a pre-defined range of Thomsen parameters, so that the phase error is small for models with spatially varying anisotropy. Phase velocity error analysis shows that the optimized pure P-wave equation is up to one order of magnitude more accurate than other popular pure P-wave equations, even for highly non-elliptic anisotropy. The optimized equation can be easily turned into a time-domain forward modeling scheme and comparisons of the modeled waveforms with analytical travel times once more illustrate its high accuracy. We also provide an efficient implementation of our approach for 3D tilted TI media that limits the count of fast Fourier transforms per time step to a number that is comparable to other pure P-wave equations.},
	keywords = {anisotropy, least squares, Modelling, SEG, TTI},
	doi = {10.1190/segam2016-13844850.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/witte2016SEGpve/witte2016SEGpve.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/witte2016SEGpve/witte2016SEGpve_poster.pdf},
	author = {Philipp A. Witte and Christiaan C. Stolk and Felix J. Herrmann}
}
@presentation {witte2016SINBADFpve,
	title = {Phase velocity error minimizing scheme for the anisotropic pure P-wave equation},
	year = {2016},
	publisher = {SINBAD},
	abstract = {Pure P-wave equations for acoustic modeling in transverse isotropic media are derived by approximating the exact pure P-wave dispersion relation. In this work, we present an alternative to the common Taylor expansion based equations, in which we approximate the exact dispersion relation through a polynomial expansion and determine its coefficients by solving a least squares problem that minimizes the phase velocity error over the complete range of phase angles. We show that this approach is up to an order of magnitude more accurate than other pure P-wave equations at a comparable computational cost.},
	keywords = {Presentation, SINBAD, SINBADFALL2016, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/witte2016SINBADFpve/witte2016SINBADFpve.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/witte2016SINBADFpve/witte2016SINBADFpve.mov},
	author = {Philipp A. Witte and Christiaan C. Stolk and Felix J. Herrmann}
}
@conference {sharan2016SEGspj,
	title = {Sparsity-promoting joint microseismic source collocation and source-time function estimation},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2016},
	note = {(SEG, Dallas)},
	month = {10},
	pages = {2574-2579},
	abstract = {In this work, we propose a new method to simultaneously locate microseismic events (e.g., induced by hydraulic fracturing) and estimate the source signature of these events. We use the method of linearized Bregman. This algorithm focuses unknown sources at their true locations by promoting sparsity along space and at the same time keeping the energy along time in check. We are particularly interested in situations where the microseismic data is noisy, sources have different signatures and we only have access to the smooth background-velocity model. We perform numerical experiments to demonstrate the usability of the proposed method. We also compare our results with full-waveform inversion based microseismic event collocation methods. Our method gives flexibility to simultaneously get a more accurate source image along with an estimate of the source-time function, which carries important information on the rupturing process and source mechanism.},
	keywords = {microseismic, noise, SEG, sources, sparse, wave equation},
	doi = {10.1190/segam2016-13871022.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/sharan2016SEGspj/sharan2016SEGspj.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/sharan2016SEGspj/sharan2016SEGspj_pres.pdf},
	author = {Shashin Sharan and Rongrong Wang and Tristan van Leeuwen and Felix J. Herrmann}
}
@presentation {witte2016SINBADFtdl,
	title = {Time-domain least-squares RTM with sparsity promotion on field data},
	year = {2016},
	publisher = {SINBAD},
	abstract = {Least squares reverse time migration (LSRTM) is an inversion based imaging algorithm which can surpress acquisition related artifacts and provides images with increased frequency content, decreased noise level and balanced amplitudes. In this work demonstrate our time-domain LSRTM workflow using various pre-conditioners, source-subsampling and sparsity promotion and show its application to the 2D Machar field data set.},
	keywords = {Presentation, SINBAD, SINBADFALL2016, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/witte2016SINBADFtdl/witte2016SINBADFtdl.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/witte2016SINBADFtdl/witte2016SINBADFtdl.mov},
	author = {Philipp A. Witte and Felix J. Herrmann}
}
@conference {yang2016SEGtds,
	title = {Time-domain sparsity-promoting least-squares migration with source estimation},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2016},
	note = {(SEG, Dallas)},
	month = {10},
	pages = {4225-4229},
	abstract = {Traditional reverse-time migration (RTM) gives images with wrong amplitudes and low resolution. Least-squares RTM (LS-RTM) on the other hand, is capable of obtaining true-amplitude images as solutions of $\ell_2$-norm norm minimization problems by fitting the synthetic and observed reflection data. The shortcoming of this approach is that solutions of these $\ell_2$ problems are typically smoothed, tend to be overfitted, and computationally too expensive because it requires compared to standard RTM too many iterations. By working with randomized subsets of data only, the computational costs of LS-RTM can be brought down to an acceptable level while producing artifact-free high-resolution images without overfitting the data. While initial results of these "compressive imaging" methods were encouraging various open issues remain including guaranteed convergence, algorithmic complexity of the solver, and lack of on-the-fly source estimation for LS-RTMs with wave-equation solvers based on time-stepping. By including on-the-fly source-time function estimation into the method of Linearized Bregman (LB), on which we reported before, we tackle all these issues resulting in a easy-to-implement algorithm that offers flexibility in the trade-off between the number of iterations and the number of wave-equation solves per iteration for a fixed total number of wave-equation solves. Application of our algorithm on a 2D synthetic shows that we are able to obtain high-resolution images, including accurate estimates of the wavelet, for a single pass through the data. The produced image, which is by virtue of the inversion deconvolved with respect to the wavelet, is roughly of the same quality as the image obtained given the correct source function.},
	keywords = {least squares, RTM, SEG, source estimation, time domain},
	doi = {10.1190/segam2016-13850609.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/yang2016SEGtds/yang2016SEGtds.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/yang2016SEGtds/yang2016SEGtds_pres.pdf},
	author = {Mengmeng Yang and Philipp A. Witte and Zhilong Fang and Felix J. Herrmann}
}
@presentation {yang2016SINBADFtds,
	title = {Time-domain sparsity-promoting least-squares migration with source estimation},
	year = {2016},
	publisher = {SINBAD},
	abstract = {Compared to traditional reverse-time migration (RTM), least-squares RTM (LS_RTM) is able to obtain true amplitude images as solutions of $\ell_2$-norm minimization problems by fitting the synthetic and observed reflection data. The shortcoming is that solutions of these $\ell_2$ problems tend to be overfitted and computationally too expensive. By working with randomized subsets of data only, the computational costs of LS-RTM can be brought down to an acceptable level, producing artifact-free high-resolution images. By including on-the-fly source-time function estimation into the method of Linearized Bregman (LB), we tackle the open issues of these "compressive imaging" in the aspects of algorithmic complexity of solver, guaranteed convergence and source estimation. We also inbvestigate whether the algorithm can be accelerated when we combine LB with the Nesterov method. Application of our algorithm on a 2D synthetic shows that we are able to get high-resolution images, with accurate estimates of the wavelet, for one single data pass.},
	keywords = {Presentation, SINBAD, SINBADFALL2016, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/yang2016SINBADFtds/yang2016SINBADFtds.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/yang2016SINBADFtds/yang2016SINBADFtds.m4v},
	author = {Mengmeng Yang and Philipp A. Witte and Zhilong Fang and Felix J. Herrmann}
}
@conference {kumar2016SEGtjm,
	title = {Time-jittered marine acquisition{\textendash}-a rank-minimization approach for 5D source separation},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2016},
	note = {(SEG, Dallas)},
	month = {10},
	pages = {119-123},
	abstract = {Simultaneous source marine acquisition has been recognized as an economic way of improving spatial sampling and speedup acquisition time, where a single- (or multiple-) source vessel fires at jittered source locations and time instances. Consequently, the acquired simultaneous data volume is processed to separate the overlapping shot records resulting in densely sampled data volume. It has been shown in the past that the simultaneous source acquisition design and source separation process can be setup as a compressed sensing problem, where conventional seismic data is reconstructed from simultaneous data via a sparsity-promoting optimization formulation. While the recovery quality of separated data is reasonably well, the recovery process can be computationally expensive due to transform-domain redundancy. In this paper, we present a computationally tractable rank-minimization algorithm to separate simultaneous data volumes. The proposed algorithm is suitable for large-scale seismic data, since it avoids singular-value decompositions and uses a low-rank based factorized formulation instead. Results are illustrated for simulations of simultaneous time-jittered continuous recording for a 3D ocean-bottom cable survey.},
	keywords = {5D, marine, SEG, source separation, time-jittered acquisition},
	doi = {10.1190/segam2016-13878249.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/kumar2016SEGtjm/kumar2016SEGtjm.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/kumar2016SEGtjm/kumar2016SEGtjm_pres.pdf},
	author = {Rajiv Kumar and Shashin Sharan and Haneet Wason and Felix J. Herrmann}
}
@presentation {wang2016SINBADFtmf,
	title = {Two methods for frequency down-extrapolation},
	year = {2016},
	publisher = {SINBAD},
	abstract = {It is well known that low frequency data are essential to the success of the Full-Waveform Inversion. Field data from a usual marine acquisition, however, typically lack reliable low frequencies due to the existing physical limitations. One research direction focuses on studying low-frequency extrapolation, the process of extending reliable frequency bands of the raw data towards the lower end of the spectrum. We propose two optimization problems, by solving which low frequencies can be made out of the high SNR high frequency component. We present both theoretical explanation and numerical experiments to demonstrate the efficacy as well as limitation of these methods.},
	keywords = {Presentation, SINBAD, SINBADFALL2016, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/wang2016SINBADFtmf/wang2016SINBADFtmf.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/wang2016SINBADFtmf/wang2016SINBADFtmf.mov},
	author = {Rongrong Wang and Felix J. Herrmann}
}
@presentation {witte2015SINBADFara,
	title = {Anisotropic RTM applied to field data},
	year = {2015},
	publisher = {SINBAD},
	abstract = {In this talk we present our latest reverse time-migration results of the BP Machar data set, a 2D seismic line from the North sea featuring a large chalk dome. In contrast to our prior imaging attempts, we now account for anisotropic effects in the data using our new TTI pure quasi-P wave modeling code in the time domain.},
	keywords = {Presentation, SINBAD, SINBADFALL2015, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/witte2015SINBADFara/witte2015SINBADFara.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/witte2015SINBADFara/witte2015SINBADFara.mov},
	author = {Philipp A. Witte and Felix J. Herrmann}
}
@conference {wason2015EAGEcsm,
	title = {Compressed sensing in 4-D marine{\textendash}-recovery of dense time-lapse data from subsampled data without repetition},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2015},
	note = {(EAGE, Madrid)},
	month = {06},
	abstract = {We present an extension of our time-jittered marine acquisition for time-lapse surveys by working on more realistic field acquisition scenarios by incorporating irregular spatial grids without insisting on repeatability between the surveys. Since we are always subsampled in both the baseline and monitor surveys, we are interested in recovering the densely sampled baseline and monitor, and then the (complete) 4-D difference from subsampled/incomplete baseline and monitor data.},
	keywords = {EAGE, NFFT, off-the-grid, simultaneous acquisition, time-lapse},
	doi = {10.3997/2214-4609.201413088},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2015/wason2015EAGEcsm/wason2015EAGEcsm.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2015/wason2015EAGEcsm/wason2015EAGEcsm_poster.pdf},
	author = {Haneet Wason and Felix Oghenekohwo and Felix J. Herrmann}
}
@presentation {wang2015SINBADFifw,
	title = {Improving full waveform inversion with spectral extrapolation},
	year = {2015},
	publisher = {SINBAD},
	abstract = {We develop a sparsity based frequency extrapolation method to obtain low frequency data from the high frequency ones. As other frequency extrapolation methods, the extrapolation introduces noise due to dispersion and possibly densely distributed reflectors. To mitigate such noise, we incorporate Curvlet coefficient regularization in the extrapolation algorithm. Numerical results demonstrate the effectiveness of the extrapolation in frequency domain FWI in models of moderate complexity.},
	keywords = {Presentation, SINBAD, SINBADFALL2015, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/wang2015SINBADFifw/wang2015SINBADFifw.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/wang2015SINBADFifw/wang2015SINBADFifw.mov},
	author = {Rongrong Wang and Felix J. Herrmann}
}
@presentation {wason2015SINBADldr,
	title = {Latest developments in randomized marine acquisition and source separation},
	year = {2015},
	publisher = {SINBAD},
	abstract = {We present the latest developments in randomized (simultaneous) marine acquisition for time-lapse surveys. We also show a comparison between two source separation algorithms {\textemdash} sparsity-promotion and rank-minimization {\textemdash} for simultaneous data.},
	keywords = {Presentation, SINBAD, SINBADSPRING2015, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Spring/wason2015SINBADldr/wason2015SINBADldr.pdf},
	author = {Haneet Wason and Rajiv Kumar and Felix Oghenekohwo and Felix J. Herrmann}
}
@conference {esser2015EAGElcs,
	title = {A lifted $\ell_1/\ell_2$ constraint for sparse blind deconvolution},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2015},
	note = {(EAGE, Madrid)},
	month = {06},
	abstract = {We propose a modification to a sparsity constraint based on the ratio of $\ell_1$ and $\ell_2$ norms for solving blind seismic deconvolution problems in which the data consist of linear convolutions of different sparse reflectivities with the same source wavelet. We also extend the approach to the Estimation of Primaries by Sparse Inversion (EPSI) model, which includes surface related multiples. Minimizing the ratio of $\ell_1$ and $\ell_2$ norms has been previously shown to promote sparsity in a variety of applications including blind deconvolution. Most existing implementations are heuristic or require smoothing the $\ell_1/\ell_2$ penalty. Lifted versions of $\ell_1/\ell_2$ constraints have also been proposed but are challenging to implement. Inspired by the lifting approach, we propose to split the sparse signals into positive and negative components and apply an $\ell_1/\ell_2$ constraint to the difference, thereby obtaining a constraint that is easy to implement without smoothing the $\ell_1$ or $\ell_2$ norms. We show that a method of multipliers implementation of the resulting model can recover source wavelets that are not necessarily minimum phase and approximately reconstruct the sparse reflectivities. Numerical experiments demonstrate robustness to the initialization as well as to noise in the data.},
	keywords = {blind deconvolution, EAGE, EPSI},
	doi = {10.3997/2214-4609.201413420},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2015/esser2015EAGElcs/esser2015EAGElcs.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2015/esser2015EAGElcs/esser2015EAGElcs_poster.pdf},
	author = {Ernie Esser and Tim T.Y. Lin and Rongrong Wang and Felix J. Herrmann}
}
@techreport {witte2015TRoam,
	title = {Overview on anisotropic modeling and inversion},
	number = {TR-EOAS-2015-6},
	year = {2015},
	month = {08},
	institution = {UBC},
	abstract = {This note provides an overview on strategies for modeling and inversion with the anisotropic wave equation. Since linear and non-linear inversion methods like least squares RTM and Full Waveform Inversion depend on matching observed field data with synthetically modelled data, accounting for anisotropy effects is necessary in order to accurately match waveforms at long offsets and propagation times. In this note, the two main strategies for anisotropic modelling by solving either a pseudo acoustic wave equation or a pure quasi-P-wave equation are discussed and an inversion workflow using the pure quasi-P-wave equation is provided. In particular, we derive the exact adjoint of the anisotropic forward modelling and jacobian operator and give a detailled describtion of their implementation. The anistropic FWI workflow is tested on a sythetic data example.},
	keywords = {anisotropy, full waveform inversion, Modeling},
	url = {https://slim.gatech.edu/Publications/Public/TechReport/2015/witte2015TRoam/witte2015TRoam.html},
	author = {Philipp A. Witte and Mathias Louboutin and Felix J. Herrmann}
}
@presentation {wason2015SINBADFrtl,
	title = {Randomization of time-lapse marine surveys},
	year = {2015},
	publisher = {SINBAD},
	abstract = {We present an extension of our time-jittered (simultaneous) marine acquisition for time-lapse surveys, where we work on more realistic field acquisition scenarios by incorporating irregular spatial grids without insisting on repeatability between the surveys. Since we are always subsampled in both the baseline and monitor surveys, we recover the densely sampled baseline and monitor and then the (complete) 4-D difference from subsampled baseline and monitor data.},
	keywords = {Presentation, SINBAD, SINBADFALL2015, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/wason2015SINBADFrtl/wason2015SINBADFrtl.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/wason2015SINBADFrtl/wason2015SINBADFrtl.mov},
	author = {Haneet Wason and Felix Oghenekohwo and Felix J. Herrmann}
}
@conference {kumar2015IIPFWIrmb,
	title = {Rank minimization based seismic data processing and inversion},
	booktitle = {Inaugural Full-Waveform Inversion Workshop},
	year = {2015},
	note = {(Natal, Brazil)},
	month = {08-09},
	abstract = {We present the low-rank extensions of full-waveform inversion to mitigate the local-minima issues. The purpose formulation is computationally tractable and does not involve any extra computational cost compared to the adjoint-state method.},
	keywords = {FWI, inversion, Processing, rank minimization},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/IIPFWI/kumar2015IIPFWIrmb/kumar2015IIPFWIrmb_pres.pdf},
	author = {Rajiv Kumar and Curt Da Silva and Oscar Lopez and Aleksandr Y. Aravkin and Hassan Mansour and Haneet Wason and Ernie Esser and Felix J. Herrmann}
}
@presentation {oghenekohwo2015SINBADrdc,
	title = {Recent developments in compressive sensing for time-lapse studies},
	year = {2015},
	publisher = {SINBAD},
	abstract = {We present an overview of the advances we have made with respect to extending compressive sensing ideas to time-lapse studies. We show applications in compressive time-lapse seismic data acquisition, imaging and full-waveform inversion.},
	keywords = {Presentation, SINBAD, SINBADSPRING2015, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Spring/oghenekohwo2015SINBADrdc/oghenekohwo2015SINBADrdc.pdf},
	author = {Felix Oghenekohwo and Rajiv Kumar and Haneet Wason and Ernie Esser and Ning Tu and Felix J. Herrmann}
}
@conference {esser2015CAMSAPrsa,
	title = {Resolving scaling ambiguities with the $\ell_1/\ell_2$ norm in a blind deconvolution problem with feedback},
	booktitle = {IEEE International Workshop on Computational Advances in Multi-Sensor Adaptive Processing},
	year = {2015},
	note = {(IEEE CAMSAP Workshop, Canc{\'u}n, Mexico)},
	month = {12},
	pages = {365-368},
	abstract = {Compared to more mundane blind deconvolution problems, blind deconvolution in seismic applications involves a feedback mechanism related to the free surface. The presence of this feedback mechanism gives us an unique opportunity to remove ambiguities that have plagued blind deconvolution for a long time. While beneficial, this feedback by itself is insufficient to remove the ambiguities even with $\ell_1$ constraints. However, when paired with an $\ell_1/\ell_2$ constraint the feedback allows us to resolve the scaling ambiguity under relatively mild assumptions. Inspired by lifting approaches, we propose to split the sparse signal into positive and negative components and apply an $\ell_1/\ell_2$ constraint to the difference, thereby obtaining a constraint that is easy to implement. Numerical experiments demonstrate robustness to the initialization as well as to noise in the data.},
	keywords = {$\ell_1/\ell_2$, blind deconvolution, CAMSAP, lifting, method of multipliers},
	doi = {10.1109/CAMSAP.2015.7383812},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/CAMSAP/2015/esser2015CAMSAPrsa/esser2015CAMSAPrsa.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/CAMSAP/2015/esser2015CAMSAPrsa/esser2015CAMSAPrsa_poster.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/CAMSAP/2015/esser2015CAMSAPrsa/esser2015CAMSAPrsa_ext.pdf},
	author = {Ernie Esser and Rongrong Wang and Tim T.Y. Lin and Felix J. Herrmann}
}
@presentation {esser2015SINBADFrsa,
	title = {Resolving scaling ambiguities with the L1/L2 norm in a blind deconvolution problem with feedback},
	year = {2015},
	publisher = {SINBAD},
	abstract = {Compared to more mundane blind deconvolution problems, blind deconvolution in seismic applications involves a feedback mechanism related to the free surface. The presence of this feedback mechanism gives us an unique opportunity to remove ambiguities that have plagued blind deconvolution for a long time. While beneficial, this feedback by itself is insufficient to remove the ambiguities even with L1 constraints. However, when paired with an L1/L2 constraint the feedback allows us to resolve the scaling ambiguity under relatively mild assumptions. Inspired by lifting approaches, we propose to split the sparse signal into positive and negative components and apply an \ell_1/\ell_2 constraint to the difference, thereby obtaining a constraint that is easy to implement. Numerical experiments demonstrate robustness to the initialization as well as to noise in the data.},
	keywords = {Presentation, SINBAD, SINBADFALL2015, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/esser2015SINBADFrsa/esser2015SINBADFrsa.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/esser2015SINBADFrsa/esser2015SINBADFrsa.mov},
	author = {Ernie Esser and Rongrong Wang and Tim T.Y. Lin and Felix J. Herrmann}
}
@presentation {dasilva2015SINBADFsss,
	title = {Scaling SINBAD software to 3-D on Yemoja},
	year = {2015},
	publisher = {SINBAD},
	abstract = {We present early results on the scalability of SINBAD{\textquoteright}s wavefield reconstruction and wave-equation based inversion technologies on Yemoja, a 17k core cluster made available to us by BG Group at SENAI CIMATEC Supercomputing Centre in Brazil.},
	keywords = {Presentation, SINBAD, SINBADFALL2015, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/dasilva2015SINBADFsss/dasilva2015SINBADFsss.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/dasilva2015SINBADFsss/dasilva2015SINBADFsss.mov},
	author = {Curt Da Silva and Haneet Wason and Mathias Louboutin and Bas Peters and Shashin Sharan and Zhilong Fang and Felix J. Herrmann}
}
@presentation {sharan2015SINBADFscu,
	title = {Source collocation using the method of Linearized Bregman},
	year = {2015},
	publisher = {SINBAD},
	abstract = {In this work, we present a method to collocate sources from seismic reflections recorded at receiver locations. We use the method of linearized Bregman. This algorithm focuses unknown sources by promoting sparsity in the {\textquotedblleft}Helmholtz domain{\textquotedblright}{\textemdash}i.e., the wavefield under the action of the the Helmholtz system{\textemdash}under the constraint of a data misfit within a particular tolerance. We extend this method to noisy measurements. We are in particular interested in situations where the noise is coherent because it is in the range of the Helmholtz equation.},
	keywords = {Presentation, SINBAD, SINBADFALL2015, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/sharan2015SINBADFscu/sharan2015SINBADFscu.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/sharan2015SINBADFscu/sharan2015SINBADFscu.mov},
	author = {Shashin Sharan and Rongrong Wang and Felix J. Herrmann}
}
@presentation {wason2015SINBADFsss,
	title = {Source separation for simultaneous towed-streamer acquisition via compressed sensing},
	year = {2015},
	publisher = {SINBAD},
	abstract = {We show a comparison between two compressed-sensing based source-separation algorithms, i.e., sparsity promotion and rank minimization, for simultaneous data.},
	keywords = {Presentation, SINBAD, SINBADFALL2015, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/wason2015SINBADFsss/wason2015SINBADFsss.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/wason2015SINBADFsss/wason2015SINBADFsss.mov},
	author = {Haneet Wason and Rajiv Kumar and Felix J. Herrmann}
}
@article {kumar2015sss,
	title = {Source separation for simultaneous towed-streamer marine acquisition {\textendash}- a compressed sensing approach},
	journal = {Geophysics},
	volume = {80},
	number = {6},
	year = {2015},
	note = {(Geophysics)},
	month = {11},
	pages = {WD73-WD88},
	abstract = {Simultaneous marine acquisition is an economic way to sample seismic data and speed up acquisition, wherein single or multiple source vessels fire sources at near-simultaneous or slightly random times, resulting in overlapping shot records. The current paradigm for simultaneous towed-streamer marine acquisition incorporates {\textquotedblleft}low variability{\textquotedblright} in source firing times, i.e., 0 <= 1 or 2 s because the sources and receivers are moving. This results in a low degree of randomness in simultaneous data, which is challenging to separate (into its constituent sources) using compressed-sensing-based separation techniques because randomization is key to successful recovery via compressed sensing. We have addressed the challenge of source separation for simultaneous towed-streamer acquisitions via two compressed-sensing-based approaches, i.e., sparsity promotion and rank minimization. We have evaluated the performance of the sparsity-promotion- and rank-minimization-based techniques by simulating two simultaneous towed-streamer acquisition scenarios, i.e., over/under and simultaneous long offset. A field data example from the Gulf of Suez for the over/under acquisition scenario was also developed. We observed that the proposed approaches gave good and comparable recovery qualities of the separated sources, but the rank-minimization technique outperformed the sparsity-promoting technique in terms of the computational time and memory. We also compared these two techniques with the normal-moveout-based median-filtering-type approach, which had comparable results.},
	keywords = {2D, Acquisition, inversion, marine, Optimization, Rank, source separation, sparsity},
	doi = {10.1190/geo2015-0108.1},
	url = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2015/kumar2015sss/kumar2015sss_revised.pdf},
	url2 = {http://library.seg.org/doi/abs/10.1190/geo2015-0108.1},
	author = {Rajiv Kumar and Haneet Wason and Felix J. Herrmann}
}
@conference {witte2015IIPFWIspl,
	title = {Sparsity-promoting least-square migration with linearized Bregman and compressive sensing},
	booktitle = {Inaugural Full-Waveform Inversion Workshop},
	year = {2015},
	note = {(Natal, Brazil)},
	month = {08-09},
	abstract = {We present a novel adaptation of a recently developed relatively simple iterative algorithm to solve large-scale sparsity-promoting optimization problems. Our algorithm is particularly suitable to large-scale geophysical inversion problems, such as sparse least-squares reverse-time migration or Kirchoff migration since it allows for a tradeoff between parallel computations, memory allocation, and turnaround times, by working on subsets of the data with different sizes. Comparison of the proposed method for sparse least-squares imaging shows a performance that rivals and even exceeds the performance of state-of-the art one-norm solvers that are able to carry out least-squares migration at the cost of a single migration with all data.},
	keywords = {Bregman, compressed sensing, FWI, migration, sparsity promotion, time domain},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/IIPFWI/witte2015IIPFWIspl/witte2015IIPFWIspl_pres.pdf},
	author = {Philipp A. Witte and Ning Tu and Ernie Esser and Mengmeng Yang and Mathias Louboutin and Felix J. Herrmann}
}
@conference {lopez2015SEGWSsfm,
	title = {SVD-free matrix completion in randomized marine acquisition},
	booktitle = {SEG Workshop on Rank-Reduction and Other Sparse Transform Methods with Application to Data Reconstruction, De-Noising, De-Blending and Imaging; New Orleans},
	year = {2015},
	note = {(SEG, New Orleans)},
	month = {10},
	abstract = {During this presentation, we will illustrate how insights from matrix-free matrix completion can be used to solve problems in (on- and off-the-grid) missing trace interpolation and time-jittered (simultaneous) marine acquisition. In particular, we will demonstrate the importance of data organization to favor recovery by rank minimization and the use of deterministic spectral properties to predict the recovery. We also make make comparisons between wavefield reconstructions based on transform-domain sparsity promotion and matrix completion.},
	keywords = {Acquisition, interpolation, matrix completion, rank reduction, SEG, workshop},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2015/lopez2015SEGWSsfm/lopez2015SEGWSsfm_pres.pdf},
	author = {Oscar Lopez and Haneet Wason and Curt Da Silva and Rajiv Kumar and Felix J. Herrmann}
}
@presentation {witte2015SINBADFtdf,
	title = {Time-domain FWI in TTI media},
	year = {2015},
	publisher = {SINBAD},
	abstract = {We develop an inversion workflow for tilted transverse isotropic (TTI) media using a purely acoustic formulation of the wave equation. The anisotropic modeling kernel is used for the forward modeling operator, as well as for the adjoint Jacobian to back propagate the data residual, thus providing the true gradient of the FWI objective function. We apply this workflow on a synthetic FWI example and perform RTM on a field data set. Furthermore, we discuss alternatives for the pseudo-spectral Laplacian operator of the current implementation.},
	keywords = {Presentation, SINBAD, SINBADFALL2015, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/witte2015SINBADFtdf/witte2015SINBADFtdf.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/witte2015SINBADFtdf/witte2015SINBADFtdf.mov},
	author = {Philipp A. Witte and Felix J. Herrmann}
}
@conference {witte2015IIPFWItdf,
	title = {Time-domain FWI in TTI media},
	booktitle = {Inaugural Full-Waveform Inversion Workshop},
	year = {2015},
	note = {(Natal, Brazil)},
	month = {08-09},
	abstract = {We develop an inversion workflow for tilted transverse isotropic (TTI) media using a purely acoustic formulation of the wave equation. The anisotropic modeling kernel is used for the forward modeling operator, as well as for the adjoint Jacobian to back propagate the data residual, thus providing the true gradient of the FWI objective function.},
	keywords = {FWI, time domain, TTI},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/IIPFWI/witte2015IIPFWItdf/witte2015IIPFWItdf_pres.pdf},
	author = {Philipp A. Witte and Mathias Louboutin and Felix J. Herrmann}
}
@techreport {kumar2015EAGEtjm,
	title = {Time-jittered marine acquisition: low-rank v/s sparsity},
	number = {TR-EOAS-2015-2},
	year = {2015},
	month = {01},
	institution = {UBC},
	abstract = {Time-jittered simultaneous marine acquisition has been recognized as an economic way of improving the spatial sampling, and speedup acquisition, where a single (or multiple) source vessel fires at {\textendash} jittered source locations and instances in time. It has been shown in the past that this problem can be setup as a {\textendash} compressed sensing problem, where conventional seismic data is reconstructed from blended data via a sparsity-promoting optimization formulation. While the recovery quality of deblended data is very good, the recovery process is computationally very expensive. In this paper, we present a computationally efficient rank-minimization algorithm to deblend the seismic data. The proposed algorithm is suitable for large-scale seismic data, since it avoids SVD computations and uses a low-rank factorized formulation instead. Results are illustrated with simulations of time-jittered marine acquisition, which translates to jittered source locations for a given speed of the source vessel, for a single source vessel with two airgun arrays.},
	keywords = {deblending, low-rank},
	url = {https://slim.gatech.edu/Publications/Public/TechReport/2015/kumar2015EAGEtjm/kumar2015EAGEtjm.html},
	author = {Rajiv Kumar and Haneet Wason and Felix J. Herrmann}
}
@presentation {lopez2015SINBADFumc,
	title = {Universal matrix completion: applications to seismic data acquisition},
	year = {2015},
	publisher = {SINBAD},
	abstract = {This talk discusses the potential applications of universal matrix completion for infill management and acquisition design. We consider recent developments in the theory of matrix completion that produce practical tools to quantify how successful a given sub sampling mask will be for seismic trace interpolation via nuclear norm minimization. This approach provides an on the fly infill management system as well as instruments to design acquisition schemes that favor interpolation via rank penalization techniques.},
	keywords = {Presentation, SINBAD, SINBADFALL2015, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/lopez2015SINBADFumc/lopez2015SINBADFumc.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/lopez2015SINBADFumc/lopez2015SINBADFumc.mov},
	author = {Oscar Lopez and Haneet Wason and Rajiv Kumar and Felix J. Herrmann}
}
@conference {wang2015SIAMwds,
	title = {Wavefield-denoising and source encoding},
	booktitle = {SIAM Conference on Mathematical and Computational Issues in the Geosciences},
	year = {2015},
	note = {(SIAM Conference on Mathematical and Computational Issues in the Geosciences, Stanford University, California)},
	month = {06-07},
	keywords = {SIAM, source encoding, WRI},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SIAM/2015/wang2015SIAMwds/wang2015SIAMwds_pres.pdf},
	author = {Rongrong Wang and Ozgur Yilmaz and Felix J. Herrmann}
}
@presentation {warner2014SINBADawi,
	title = {Adaptive waveform inversion - FWI without cycle skipping},
	year = {2014},
	publisher = {SINBAD},
	abstract = {Conventional FWI minimises the direct differences between observed and predicted seismic datasets. Because seismic data are oscillatory, this approach will suffer from the detrimental effects of cycle skipping if the starting model is inaccurate. We reformulate FWI so that it instead adapts the predicted data to the observed data using Wiener filters, and then iterates to improve the model by forcing the Wiener filters towards zero-lag delta functions. This adaptive FWI scheme is demonstrated on synthetic data where it is shown to be immune to cycle skipping, and is able to invert successfully data for which conventional FWI fails entirely. The new method does not require low frequencies or a highly accurate starting model to be successful. Adaptive FWI has some features in common with wave-equation migration velocity analysis, but it works for all types of arrivals including multiples and refractions, and it does not have the high computational costs of WEMVA in 3D.},
	keywords = {Presentation, SINBAD, SINBADSPRING2014, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2014/Spring/warner2014SINBADawi/warner2014SINBADawi.pdf},
	author = {Mike Warner}
}
@presentation {wang2014SINBADdwi,
	title = {Denoising the wavefield inversion problem through source blending and penalty method},
	year = {2014},
	publisher = {SINBAD},
	abstract = {Previous FWI denoising techniques, such as source blending and stacking, are designed applied to the seismic traces. Hence no physical information is used to improve the fidelity of the reconstruction, even when a good initial model is available. Moreover, traditional source blending requires the existence of common source gathers, which might be missing from the data. In this talk, we incorporate the source blending idea into the penalty method, and extract the synthetic sources corresponding to the principle directions of the gradient update. By blending the sources, we can not only achieve an acceleration and storage efficiency of the penalty method, but also make the reconstruction more robust to the random noise. This method is applicable to any source receiver position but is more efficient if common source gathers are available.},
	keywords = {Presentation, SINBAD, SINBADFALL2014, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2014/Fall/wang2014SINBADdwi/wang2014SINBADdwi.pdf},
	author = {Rongrong Wang and Felix J. Herrmann}
}
@techreport {wang2014SEGfwi,
	title = {Full waveform inversion with interferometric measurements},
	number = {TR-EOAS-2014-5},
	year = {2014},
	month = {04},
	institution = {UBC},
	abstract = {In this note, we design new misfit functions for full-waveform inversion by using interferometric measurements to reduce sensitivity to phase errors. Though established within a completely different setting from the linear case, we obtain a similar observation: the interferometry can improve robustness under certain modeling errors. Moreover, in order to deal with errors on both source and receiver sides, we propose a higher order interferometry, which, as a generalization of the usual definition, involves the cross correlation of four traces. A proof of principle simulations is included on a stylized example.},
	keywords = {FWI},
	url = {https://slim.gatech.edu/Publications/Public/TechReport/2014/wang2014SEGfwi/wang2014SEGfwi.html},
	author = {Rongrong Wang and Ozgur Yilmaz and Felix J. Herrmann}
}
@presentation {esser2014SINBADlcs,
	title = {A lifted $\ell_1$/$\ell_2$ constraint for sparse blind deconvolution},
	year = {2014},
	publisher = {SINBAD},
	abstract = {We propose a modification to a sparsity constraint based on the ratio of l1 and l2 norms for solving blind seismic deconvolution problems in which the data consist of linear convolutions of different sparse reflectivities with the same source wavelet. No assumptions are made about the location of the support of either the wavelet or the sparse signals. Minimizing the ratio of l1 and l2 norms has been previously shown to promote sparsity in a variety of applications including blind deconvolution. Most existing implementations are heuristic or require smoothing the l1/l2 penalty. Lifted versions of l1/l2 constraints have also been proposed but are still challenging to implement. Inspired by the lifting approach, we propose to split the sparse signals into positive and negative components and apply an l1/l2 constraint to the difference, thereby obtaining a constraint that is easy to implement without smoothing the l1 or l2 norms. We show that a method of multipliers implementation of the resulting model can recover source wavelets that are not necessarily minimum phase and approximately reconstruct the sparse reflectivities. It appears to be robust to the initialization and to small amounts of noise in the data. We also discuss extensions to the Estimation of Primaries by Sparse Inversion (EPSI) convolution model.},
	keywords = {Presentation, SINBAD, SINBADFALL2014, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2014/Fall/esser2014SINBADlcs/esser2014SINBADlcs.pdf},
	author = {Ernie Esser and Tim T.Y. Lin and Rongrong Wang and Felix J. Herrmann}
}
@conference {wason2014SEGrrt,
	title = {Randomization and repeatability in time-lapse marine acquisition},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2014},
	note = {(SEG)},
	month = {10},
	pages = {46-51},
	abstract = {We present an extension of our time-jittered simultaneous marine acquisition to time-lapse surveys where the requirement for repeatability in acquisition can be waived provided we know the acquisition geometry afterwards. Our method, which does not require repetition, gives 4-D signals comparable to conventional methods where repeatability is key to their success.},
	keywords = {Acquisition, deblending, marine, SEG, time-laspe},
	doi = {http://dx.doi.org/10.1190/segam2014-1677.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2014/wason2014SEGrrt/wason2014SEGrrt.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2014/wason2014SEGrrt/wason2014SEGrrt_pres.pdf},
	author = {Haneet Wason and Felix Oghenekohwo and Felix J. Herrmann}
}
@conference {herrmann2014EAGEWSrrt,
	title = {Randomization and repeatability in time-lapse marine acquisition},
	booktitle = {EAGE Workshop on Land and Ocean Bottom; Broadband Full Azimuth Seismic Surveys; Spain},
	year = {2014},
	note = {(EAGE Workshop, Spain)},
	month = {04},
	keywords = {4D seismic, EAGE, Marine acquisition, workshop},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2014/herrmann2014EAGEWSrrt.pdf},
	author = {Haneet Wason and Felix Oghenekohwo and Felix J. Herrmann}
}
@presentation {wason2014SINBADrrt,
	title = {Randomization and repeatability in time-lapse marine acquisition},
	year = {2014},
	publisher = {SINBAD},
	abstract = {During this talk, we will show how randomization techniques from (distributed) compressed sensing affect the way we think about time-lapse surveys.},
	keywords = {Presentation, SINBAD, SINBADFALL2014, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2014/Fall/wason2014SINBADrrt/wason2014SINBADrrt.pdf},
	author = {Haneet Wason and Felix Oghenekohwo and Felix J. Herrmann}
}
@presentation {oghenekohwo2014SINBADrsw,
	title = {Randomized sampling without repetition in time-lapse seismic surveys},
	year = {2014},
	publisher = {SINBAD},
	abstract = {In this talk, we will show a method for acquiring time-lapse data, where we do not have to repeat the survey geometry. Our method works provided our acquisition is randomized, where sources and receivers are at random locations on a computational grid, and provided we know the spatial locations of the shots and receivers. In addition, we show the implications of either (sub)sampling more for the baseline and less for the monitor or vice-versa, and how this sampling scheme affects the 4-D signal.},
	keywords = {Presentation, SINBAD, SINBADSPRING2014, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2014/Spring/oghenekohwo2014SINBADrsw/oghenekohwo2014SINBADrsw.pdf},
	author = {Felix Oghenekohwo and Haneet Wason and Ernie Esser and Felix J. Herrmann}
}
@presentation {oghenekohwo2014SINBADrst,
	title = {Randomized subsampling in time-lapse surveys and recovery techniques},
	year = {2014},
	publisher = {SINBAD},
	abstract = {This talk centres on new insights on how we think of acquisition of time-lapse data using ideas from compressed sensing. Specifically, we will show how simultaneous processing of time-lapse data from randomized subsampling gives more reliable estimates of the vintages and time-lapse signal, compared to parallel processing.},
	keywords = {Presentation, SINBAD, SINBADFALL2014, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2014/Fall/oghenekohwo2014SINBADrst/oghenekohwo2014SINBADrst.pdf},
	author = {Felix Oghenekohwo and Rajiv Kumar and Haneet Wason and Ernie Esser and Felix J. Herrmann}
}
@conference {wason2014SEGsss,
	title = {Source separation via SVD-free rank minimization in the hierarchical semi-separable representation},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2014},
	note = {(SEG)},
	month = {10},
	pages = {120-126},
	abstract = {Recent developments in matrix rank optimization have allowed for new computational approaches in the field of source separation or deblending. In this paper, we propose a source separation algorithm for blended marine acquisition, where two sources are deployed at different depths (over/under acquisition). The separation method incorporates the Hierarchical Semi-Separable structure (HSS) inside rank-regularized least-squares formulations. The proposed approach is suitable for large scale problems, since it avoids SVD computations and uses a low-rank factorized formulation instead. We illustrate the performance of the new HSS-based deblending approach by simulating an over/under blended acquisition, wherein uniformly random time delays (of < 1 second) are applied to one of the sources.},
	keywords = {Acquisition, deblending, HSS, marine, Rank, SEG, source separation},
	doi = {http://dx.doi.org/10.1190/segam2014-1583.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2014/wason2014SEGsss/wason2014SEGsss.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2014/wason2014SEGsss/wason2014SEGsss_pres.pdf},
	author = {Haneet Wason and Rajiv Kumar and Aleksandr Y. Aravkin and Felix J. Herrmann}
}
@presentation {wason2014SINBADsss,
	title = {Source separation via SVD-free rank-minimization in the hierarchical semi-separable representation},
	year = {2014},
	publisher = {SINBAD},
	abstract = {During this talk, we will show a source separation algorithm for blended marine acquisition, where two sources are deployed at different depths (over/under acquisition). The separation method incorporates the Hierarchical Semi-Separable (HSS) structure inside rank-regularized least-squares formulations. We also compare this deblending scheme with the sparsity-promoting one-norm minimization scheme.},
	keywords = {Presentation, SINBAD, SINBADFALL2014, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2014/Fall/wason2014SINBADsss/wason2014SINBADsss.pdf},
	author = {Haneet Wason and Rajiv Kumar and Felix J. Herrmann}
}
@presentation {kumar2014SINBADtjm,
	title = {Time-jittered marine acquisition: low-rank v/s sparsity},
	year = {2014},
	publisher = {SINBAD},
	abstract = {In this work we will show how simultaneous or blended acquisition can be setup as a rank-minimization (deblending) problem. The simultaneous acquisition scenario is a pragmatic time-jittered marine acquisition scheme, where a source vessel (with two airgun arrays) sails across an ocean-bottom array firing at jittered source locations and instances in time, resulting in better spatial sampling and acquisition speedup. We make comparisons with sparsity-promoting based techniques and demonstrate that rank-minimization based debelending techniques are computationally faster and memory efficient. This is joint work with Haneet Wason.},
	keywords = {Presentation, SINBAD, SINBADFALL2014, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2014/Fall/kumar2014SINBADtjm/kumar2014SINBADtjm.pdf},
	author = {Rajiv Kumar and Haneet Wason and Felix J. Herrmann}
}
@conference {lin2013SEGdss,
	title = {Dense shot-sampling via time-jittered marine sources},
	booktitle = {SEG Workshop on Simultaneous Sources; Houston},
	year = {2013},
	month = {09},
	keywords = {Acquisition, SEG, simultaneous sources, workshop},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2013/lin2013SEGdss/lin2013SEGdss_pres.pdf},
	author = {Tim T.Y. Lin and Haneet Wason and Felix J. Herrmann}
}
@article {herrmann2013TLEffwi,
	title = {Frugal full-waveform inversion: from theory to a practical algorithm},
	journal = {The Leading Edge},
	volume = {32},
	number = {9},
	year = {2013},
	month = {09},
	pages = {1082-1092},
	abstract = {As conventional oil and gas fields are maturing, our profession is challenged to come up with the next-generation of more and more sophisticated exploration tools. In exploration seismology this trend has led to the emergence of wave-equation based inversion technologies such as reverse-time migration and full-waveform inversion. While significant progress has been made in wave-equation based inversion, major challenges remain in the development of robust and computationally feasible workflows that give reliable results in geophysically challenging areas that may include ultra-low shear velocity zones or high-velocity salt. Moreover, sub-salt production carries risks that needs mitigation, which raises the bar from creating sub-salt images to inverting for sub-salt overpressure.},
	keywords = {Optimization, waveform inversion},
	doi = {10.1190/tle32091082.1},
	url = {https://slim.gatech.edu/Publications/Public/Journals/TheLeadingEdge/2013/herrmann2013ffwi/herrmann2013ffwi.html},
	author = {Felix J. Herrmann and Andrew J. Calvert and Ian Hanlon and Mostafa Javanmehri and Rajiv Kumar and Tristan van Leeuwen and Xiang Li and Brendan R. Smithyman and Eric Takam Takougang and Haneet Wason}
}
@presentation {li2013SINBADgom,
	title = {Lessons learned from Chevron Gulf of Mexico data set},
	year = {2013},
	publisher = {SINBAD},
	abstract = {Chevron Gulf of Mexico data set is very challenging for FWI because of elastic phases, limited offset, lack of low frequencies and salt structure. To overcome these issue, we first use ray-based tomography on the hand-picked first breaks to generate initial model for FWI, and then we apply curvelet-denosing techniques to improve the poor signal-to-noise ratio of the observed data at low frequencies. Finally, Curvelet domain sparsity promoting Gauss-Newton FWI helps to suppress model space artifacts caused by elastic phases. This is joint work with Andrew J. Calvert, Ian Hanlon, Mostafa Javanmehri, Rajiv Kumar, Tristan van Leeuwen, Brendan R. Smithyman, Haneet Wason and Felix J. Herrmann},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/li2013SINBADgom/li2013SINBADgom.pdf},
	author = {Xiang Li and Andrew J. Calvert and Ian Hanlon and Mostafa Javanmehri and Rajiv Kumar and Tristan van Leeuwen and Brendan R. Smithyman and Haneet Wason and Felix J. Herrmann}
}
@presentation {wang2013SINBADnru,
	title = {Noise reduction by using interferometric measurements},
	year = {2013},
	publisher = {SINBAD},
	abstract = {The interferometric formulation of linear wave-based inversion problem was proposed by Demanet and Jugnon recently. Instead of directly fitting the data, they proposed to fit a subset of the data{\textquoteright}s cross-correlation. It can be verified that if the full cross-correlation is used, then the problem is equivalent to the usual least square problem. The subsampling, which is usually considered to cause instability to the solution, is surprisingly useful in this setting. Numerical experiments for the inverse source problem and the inverse scatting problem have both suggested that a {\textquoteright}good{\textquoteright} sampling strategy can actually increase the stability under modeling error caused by the uncertainty of a kinetic nature. We will study the mathematical mechanism behind this phenomenon, and try to see whether or not there exists a universally {\textquoteright}good{\textquoteright} sampling strategy independent of the types of forward operators we use.},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/wang2013SINBADnru/wang2013SINBADnru.pdf},
	author = {Rongrong Wang and Bas Peters and Felix J. Herrmann}
}
@conference {wason2013EAGEobs,
	title = {Ocean bottom seismic acquisition via jittered sampling},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2013},
	month = {06},
	abstract = {We present a pragmatic marine acquisition scheme where multiple source vessels sail across an ocean-bottom array firing at airgunsjittered source locations and instances in time. Following the principles of compressive sensing, we can significantly impact the reconstruction quality of conventional seismic data (from jittered data) and demonstrate successful recovery by sparsity promotion. In contrast to random (under)sampling, acquisition via jittered (under)sampling helps in controlling the maximum gap size, which is a practical requirement of wavefield reconstruction with localized sparsifying transforms. Results are illustrated with simulations of time-jittered marine acquisition, which translates to jittered source locations for a given speed of the source vessel, for two source vessels.},
	keywords = {Acquisition, blended, deblending, EAGE, interpolation, marine},
	doi = {10.3997/2214-4609.20130379},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2013/wason2013EAGEobs/wason2013EAGEobs.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2013/wason2013EAGEobs/wason2013EAGEobs_pres.pdf},
	author = {Haneet Wason and Felix J. Herrmann}
}
@presentation {warner2013SINBADrfwi,
	title = {Reflection FWI with a poor starting model},
	year = {2013},
	publisher = {SINBAD},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/warner2013SINBADrfwi/warner2013SINBADrfwi.pdf},
	author = {Mike Warner}
}
@presentation {wason2013SINBADtjm,
	title = {Time-jittered marine sources},
	year = {2013},
	publisher = {SINBAD},
	abstract = {Current efforts towards dense shot (and/or receiver) sampling and full azimuthal coverage to produce higher-resolution images have led to the deployment of multiple source vessels across the survey area. A step ahead from multi-source seismic acquisition is simultaneous or blended acquisition where different source arrays/vessels fire shots at near-simultaneous or slightly random times. Seismic data acquisition with simultaneous (or blended) sources has helped improve acquisition efficiency and mitigate acquisition related costs. Deblending then aims to recover unblended data, as acquired during conventional acquisition, from blended data since many processing techniques rely on full, regular sampling. We present a simultaneous/blended marine acquisition setup where shots fire at significantly jittered instances in time resulting in jittered shot locations for a given speed of the source vessel. The conventional, unblended data is recovered from the blended, jittered/irregular data by sparsity-promoting inversion using the non-equispaced fast discrete curvelet transform. The optimization scheme aims to deblend the blended data along with regularization and interpolation to a (finer) regular grid.},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/wason2013SINBADtjm/wason2013SINBADtjm.pdf},
	author = {Haneet Wason and Felix J. Herrmann}
}
@conference {wason2013SEGtjo,
	title = {Time-jittered ocean bottom seismic acquisition},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {32},
	year = {2013},
	month = {09},
	pages = {1-6},
	abstract = {Leveraging ideas from the field of compressed sensing, we show how simultaneous or blended acquisition can be setup as a {\textendash} compressed sensing problem. This helps us to design a pragmatic time-jittered marine acquisition scheme where multiple source vessels sail across an ocean-bottom array firing airguns at {\textendash} jittered source locations and instances in time, resulting in better spatial sampling, and speedup acquisition. Furthermore, we can significantly impact the reconstruction quality of conventional seismic data (from jittered data) and demonstrate successful recovery by sparsity promotion. In contrast to random (under)sampling, acquisition via jittered (under)sampling helps in controlling the maximum gap size, which is a practical requirement of wavefield reconstruction with localized sparsifying transforms. Results are illustrated with simulations of time-jittered marine acquisition, which translates to jittered source locations for a given speed of the source vessel, for two source vessels.},
	keywords = {Acquisition, blending, deblending, interpolation, jittered sampling, marine, OBC, SEG},
	doi = {10.1190/segam2013-1391.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2013/wason2013SEGtjo/wason2013SEGtjo.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2013/wason2013SEGtjo/wason2013SEGtjo_pres.pdf},
	author = {Haneet Wason and Felix J. Herrmann}
}
@presentation {wason2013SINBADtjo,
	title = {Time-jittered ocean bottom seismic acquisition},
	year = {2013},
	publisher = {SINBAD},
	abstract = {Leveraging ideas from the field of compressed sensing, we show how simultaneous or blended acquisition can be setup as a {\textendash} compressed sensing problem. This helps us to design a pragmatic time-jittered marine acquisition scheme where multiple source vessels sail across an ocean-bottom array firing airguns at {\textendash} jittered source locations and instances in time, resulting in better spatial sampling, and speedup acquisition. Furthermore, we can significantly impact the reconstruction quality of conventional seismic data (from jittered data) and demonstrate successful recovery by sparsity promotion. In contrast to random (under)sampling, acquisition via jittered (under)sampling helps in controlling the maximum gap size, which is a practical requirement of wavefield reconstruction with localized sparsifying transforms. Results are illustrated with simulations of time-jittered marine acquisition, which translates to jittered source locations for a given speed of the source vessel.},
	keywords = {Presentation, SINBAD, SINBADSPRING2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Spring/wason2013SINBADtjo/wason2013SINBADtjo.pdf},
	author = {Haneet Wason and Felix J. Herrmann}
}
@presentation {warner2012SINBADafw,
	title = {Anisotropic 3D full-waveform inversion of the Tommeliten Alpha field},
	year = {2012},
	publisher = {SINBAD},
	abstract = {We have implemented a robust and practical scheme for anisotropic 3D acoustic full-waveform inversion. We demonstrate this scheme on a field data set, applying it to a four-component ocean-bottom survey over the Tommeliten Alpha field in the North Sea. This shallow-water data set provides good azimuthal coverage to offsets of 7 km, with reduced coverage to a maximum offset of about 11 km. The reservoir lies at the crest of a high-velocity antiformal chalk section, overlain by about 3000 m of clastics within which a low-velocity gas cloud produces a seismic obscured area. We inverted only the hydrophone data, and we retained free-surface multiples and ghosts within the field data. We invert in six narrow frequency bands, in the range 3 to 6.5 Hz. At each iteration, we selected only a subset of sources, using a different subset at each iteration; this strategy is more efficient than inverting all the data every iteration. Our starting velocity model was obtained using standard PSDM model building including anisotropic reflection tomography, and contained epsilon values as high as 20\%. We have also attempted full-elastic inversion of these data to recover a shallow isotropic model of both p and s-wave velocities. The final FWI velocity model shows a network of shallow high-velocity channels that match similar features in the reflection data. Deeper in the section, the FWI velocity model reveals a sharper and more-intense low-velocity region associated with the gas cloud in which low-velocity fingers match the location of gas-filled faults visible in the reflection data. The resulting velocity model provides a better match to well logs, and better flattens common image gathers, than does the starting model. Reverse-time migration, using the FWI velocity model, provides significant uplift to the migrated image, simplifying the planform of the reservoir section at depth. The workflows, inversion strategy, and algorithms that we have used have broad application to invert a wide-range of analogous field data sets.},
	keywords = {Presentation, SINBAD, SINBADFALL2012, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2012/Fall/warner2012SINBADafw/warner2012SINBADafw_pres.pdf},
	author = {Mike Warner}
}
@conference {herrmann2012EAGEcsm,
	title = {Compressive sensing in marine acquisition and beyond},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2012},
	month = {06},
	abstract = {Simultaneous-source marine acquisition is an example of compressive sensing where acquisition with a single vessel is replaced by simultaneous acquisition by multiple vessels with sources that fire at randomly dithered times. By identifying simultaneous acquisition as compressive sensing, we are able to design acquisitions that favour recovery by sparsity promotion. Compared to conventional processing that yields estimates for sequential data, sparse recovery leads to significantly improved results for simultaneous data volumes that are collected in shorter times. These improvements are the result of proper design of the acquisition, selection of the appropriate transform domain, and solution of the recovery problem by sparsity promotion. During this talk, we will show how these design principles can be applied to marine acquisition and to other problems in exploration seismology that can benefit from compressive sensing.},
	keywords = {Acquisition, EAGE, marine, workshop},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2012/herrmann2012EAGEcsm/herrmann2012EAGEcsm.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2012/herrmann2012EAGEcsm/herrmann2012EAGEcsm_pres.pdf},
	url2 = {http://earthdoc.eage.org/publication/publicationdetails/?publication=59854},
	author = {Felix J. Herrmann and Haneet Wason}
}
@conference {herrmann2012SEGfwi,
	title = {FWI from the West Coasts: lessons learned from "Gulf of Mexico Imaging Challenges: What Can Full Waveform Inversion Achieve?"},
	booktitle = {SEG Workshop on FWI; Las Vegas},
	year = {2012},
	publisher = {SEG},
	organization = {SEG},
	keywords = {FWI, SEG, workshop},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2012/herrmann2012SEGfwi/herrmann2012SEGfwi_pres.pdf},
	author = {Andrew J. Calvert and Ian Hanlon and Mostafa Javanmehri and Rajiv Kumar and Tristan van Leeuwen and Xiang Li and Brendan R. Smithyman and Eric Takam Takougang and Haneet Wason and Felix J. Herrmann}
}
@conference {vanderneut2012SEGirs,
	title = {Interferometric redatuming with simultaneous and missing sources using sparsity promotion in the curvelet domain},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {31},
	year = {2012},
	month = {11},
	pages = {1-7},
	publisher = {SEG},
	organization = {SEG},
	abstract = {Interferometric redatuming is a velocity-independent method to turn downhole receivers into virtual sources. Accurate redatuming involves solving an inverse problem, which can be highly ill-posed, especially in the presence of noise, incomplete data and limited aperture. We address these issues by combining interferometric redatuming with transform-domain sparsity promotion, leading to a formulation that deals with data imperfections. We show that sparsity promotion improves the retrieval of virtual shot records under a salt flank. To reduce acquisition costs, it can be beneficial to reduce the number of sources or shoot them simultaneously. It is shown that sparse inversion can still provide a stable solution in such cases.},
	keywords = {Imaging, Interferometry, Optimization, Processing, SEG},
	doi = {10.1190/segam2012-0566.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2012/vanderneut2012SEGirs/vanderneut2012SEGirs.pdf},
	author = {Joost van der Neut and Felix J. Herrmann and Kees Wapenaar}
}
@presentation {wason2012SINBADobs,
	title = {Ocean bottom seismic acquisition via jittered sampling},
	year = {2012},
	publisher = {SINBAD},
	abstract = {We present a pragmatic marine acquisition scheme where a single (or multiple) vessel sails across an ocean-bottom array firing airguns at ? optimally jittered source locations and instances in time. Following the principles of compressive sensing, we can significantly impact the reconstruction quality of conventional seismic data (from jittered data) and demonstrate successful recovery by sparsity promotion. In contrast to random (under)sampling, acquisition via jittered (under)sampling helps in controlling the maximum gap size, which is a practical requirement of wavefield reconstruction with localized sparsifying transforms. Results are illustrated with simulations of optimally jittered marine acquisition, and periodic time-dithering marine acquisition.},
	keywords = {Presentation, SINBAD, SINBADFALL2012, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2012/Fall/wason2012SINBADobs/wason2012SINBADobs_pres.pdf},
	author = {Haneet Wason}
}
@conference {wason2012CSEGode,
	title = {Only dither: efficient simultaneous marine acquisition},
	booktitle = {CSEG Annual Conference Proceedings},
	year = {2012},
	month = {02},
	abstract = {Simultaneous-source acquisition is an emerging technology that is stimulating both geophysical research and commercial efforts. Simultaneous marine acquisition calls for the development of a new set of design principles and post-processing tools. The focus here is on simultaneous-source marine acquisition design and sparsity-promoting sequential-source data recovery. We propose a pragmatic simultaneous-source, randomized marine acquisition scheme where multiple vessels sail across an ocean-bottom array firing airguns at {\textendash}- sequential locations and randomly time-dithered instances. By leveraging established findings from the field of compressive sensing, where the choice of the sparsifying transform needs to be incoherent with the compressive sampling matrix, we can significantly impact the reconstruction quality, and demonstrate that the compressive sampling matrix resulting from the proposed sampling scheme is sufficiently incoherent with the curvelet transform to yield successful recovery by sparsity promotion. Results are illustrated with simulations of {\textquotedblleft}purely{\textquotedblright} random marine acquisition, which requires an airgun to be located at each source location, and random time-dithering marine acquisition with one and two source vessels. Size of the collected data volumes in all cases is the same. Compared to the recovery from the former acquisition scheme (SNR = 10.5dB), we get good results by dithering with only one source vessel (SNR = 8.06dB) in the latter scheme, which improve at the cost of having an additional source vessel (SNR = 9.85dB).},
	keywords = {Acquisition, CSEG, marine, simultaneous},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/CSEG/2012/wason2012CSEGode/wason2012CSEGode.pdf},
	author = {Haneet Wason and Felix J. Herrmann}
}
@conference {wason2012EAGEode,
	title = {Only dither: efficient simultaneous marine acquisition},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2012},
	month = {06},
	abstract = {Simultaneous-source acquisition is an emerging technology that is stimulating both geophysical research and commercial efforts. The focus here is on simultaneous-source marine acquisition design and sparsity-promoting sequential-source data recovery. We propose a pragmatic simultaneous-source, randomized marine acquisition scheme where multiple vessels sail across an ocean-bottom array firing airguns at {\textendash}- sequential locations and randomly time-dithered instances. Within the context of compressive sensing, where the choice of the sparsifying transform needs to be incoherent with the compressive sampling matrix, we can significantly impact the reconstruction quality, and demonstrate that the compressive sampling matrix resulting from the proposed sampling scheme is sufficiently incoherent with the curvelet transform to yield successful recovery by sparsity promotion. Results are illustrated with simulations of {\textquoteleft}{\textquoteleft}purely" random marine acquisition, which requires an airgun to be located at each source location, and random time-dithering marine acquisition with one and two source vessels. Size of the collected data volumes in all cases is the same. Compared to the recovery from the former acquisition scheme (SNR = 10.5dB), we get good results by dithering with only one source vessel (SNR = 8.06dB) in the latter scheme, which improve at the cost of having an additional source vessel (SNR = 9.44dB).},
	keywords = {Acquisition, EAGE, marine, simultaneous},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2012/wason2012EAGEode/wason2012EAGEode.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2012/wason2012EAGEode/wason2012EAGEode_pres.pdf},
	url2 = {http://earthdoc.eage.org/publication/publicationdetails/?publication=58915},
	author = {Haneet Wason and Felix J. Herrmann}
}
@presentation {herrmann2012SINBADfwi,
	title = {Our findings on the Chevron benchmark dataset},
	year = {2012},
	publisher = {SINBAD},
	abstract = {During this presentation, we will review our findings working with the synthetic GOM data released as part of the post SEG workshop:"Gulf of Mexico Imaging Challenges: What Can Full Waveform Inversion Achieve?". This is joint work with Andrew J. Calvert, Ian Hanlon, Mostafa Javanmehri, Rajiv Kumar, Tristan van Leeuwen, Xiang Li, Brendan R. Smithyman, Eric Takam Takougang, and Haneet Wason.},
	keywords = {Presentation, SINBAD, SINBADFALL2012, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2012/Fall/herrmann2012SINBADfwi/herrmann2012SINBADfwi_pres.pdf},
	author = {Andrew J. Calvert and Ian Hanlon and Mostafa Javanmehri and Rajiv Kumar and Tristan van Leeuwen and Xiang Li and Brendan R. Smithyman and Eric Takam Takougang and Haneet Wason and Felix J. Herrmann}
}
@presentation {wason2012SINBADrma,
	title = {Randomized marine acquisition for ocean-bottom surveys},
	year = {2012},
	publisher = {SINBAD},
	keywords = {Presentation, SINBAD, SINBADSPRING2012, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2012/Spring/wason2012SINBADrma/wason2012SINBADrma_pres.pdf},
	author = {Haneet Wason}
}
@article {Mansour11TRssma,
	title = {Randomized marine acquisition with compressive sampling matrices},
	journal = {Geophysical Prospecting},
	volume = {60},
	number = {4},
	year = {2012},
	month = {07},
	pages = {648-662},
	abstract = {Seismic data acquisition in marine environments is a costly process that calls for the adoption of simultaneous-source or randomized acquisition - an emerging technology that is stimulating both geophysical research and commercial efforts. Simultaneous marine acquisition calls for the development of a new set of design principles and post-processing tools. In this paper, we discuss the properties of a specific class of randomized simultaneous acquisition matrices and demonstrate that sparsity-promoting recovery improves the quality of reconstructed seismic data volumes. We propose a practical randomized marine acquisition scheme where the sequential sources fire airguns at only randomly time-dithered instances. We demonstrate that the recovery using sparse approximation from random time-dithering with a single source approaches the recovery from simultaneous-source acquisition with multiple sources. Established findings from the field of compressive sensing indicate that the choice of the sparsifying transform that is incoherent with the compressive sampling matrix can significantly impact the reconstruction quality. Leveraging these findings, we then demonstrate that the compressive sampling matrix resulting from our proposed sampling scheme is incoherent with the curvelet transform. The combined measurement matrix exhibits better isometry properties than other transform bases such as a non-localized multidimensional Fourier transform. We illustrate our results with simulations of "ideal" simultaneous-source marine acquisition, which dithers both in time and space, compared with periodic and randomized time-dithering.},
	keywords = {curvelet transform, Fourier, Marine acquisition},
	doi = {10.1111/j.1365-2478.2012.01075.x},
	url = {http://onlinelibrary.wiley.com/doi/10.1111/j.1365-2478.2012.01075.x/abstract},
	url2 = {https://slim.gatech.edu/Publications/Public/Journals/GeophysicalProspecting/2012/Mansour11TRssma/Mansour11TRssma.pdf},
	author = {Hassan Mansour and Haneet Wason and Tim T.Y. Lin and Felix J. Herrmann}
}
@article {herrmann2011RECORDERcsse2,
	title = {Compressive sensing in seismic exploration: an outlook on a new paradigm},
	journal = {CSEG Recorder},
	volume = {36},
	number = {6},
	year = {2011},
	note = {Part 1 was published in April and Part 2 was published in June},
	month = {06},
	pages = {34-39},
	abstract = {Many seismic exploration techniques rely on the collection of massive data volumes that are subsequently mined for information during processing. While this approach has been extremely successful in the past, current efforts toward higher resolution images in increasingly complicated regions of the Earth continue to reveal fundamental shortcomings in our workflows. Chiefly amongst these is the so-called "curse of dimensionality" exemplified by Nyquist{\textquoteright}s sampling criterion, which disproportionately strains current acquisition and processing systems as the size and desired resolution of our survey areas continues to increase. We offer an alternative sampling method leveraging recent insights from compressive sensing towards seismic acquisition and processing for data that, from a traditional point of view, are considered to be undersampled. The main outcome of this approach is a new technology where acquisition and processing related costs are decoupled the stringent Nyquist sampling criterion. At the heart of our approach lies randomized incoherent sampling that breaks subsampling-related interferences by turning them into harmless noise, which we subsequently remove by promoting sparsity in a transform-domain. Acquisition schemes designed to fit into this regime no longer grow significantly in cost with increasing resolution and dimensionality of the survey area, but instead its cost ideally only depends on transform-domain sparsity of the expected data. Our contribution is twofold. First, we demonstrate by means of carefully designed numerical experiments that ideas from compressive sensing can be adapted to seismic acquisition. Second, we leverage the property that seismic data volumes are well approximated by a small percentage of curvelet coefficients. Thus curvelet-domain sparsity allows us to recover conventionally-sampled seismic data volumes from compressively-sampled data volumes whose size exceeds this percentage by only a small factor. Because compressive sensing combines transformation and encoding by a single linear encoding step, this technology is directly applicable to seismic acquisition and therefore constitutes a new paradigm where acquisitions costs scale with transform-domain sparsity instead of the gridsize. We illustrate this principle by showcasing recovery of a real seismic line from simulated compressively sampled acquisitions.},
	url = {http://csegrecorder.com/articles/view/compressive-sensing-in-seismic-exploration-an-outlook-on-a-new-paradigm},
	url2 = {https://slim.gatech.edu/Publications/Public/Journals/CSEGRecorder/2011/herrmann2011RECORDERcsse/herrmann2011RECORDERcsse.pdf},
	author = {Felix J. Herrmann and Haneet Wason and Tim T.Y. Lin}
}
@article {herrmann2011RECORDERcsse1,
	title = {Compressive sensing in seismic exploration: an outlook on a new paradigm},
	journal = {CSEG Recorder},
	volume = {36},
	number = {4},
	year = {2011},
	note = {Part 1 was published in April and Part 2 was published in June},
	month = {04},
	pages = {19-33},
	abstract = {Many seismic exploration techniques rely on the collection of massive data volumes that are subsequently mined for information during processing. While this approach has been extremely successful in the past, current efforts toward higher resolution images in increasingly complicated regions of the Earth continue to reveal fundamental shortcomings in our workflows. Chiefly amongst these is the so-called "curse of dimensionality" exemplified by Nyquist{\textquoteright}s sampling criterion, which disproportionately strains current acquisition and processing systems as the size and desired resolution of our survey areas continues to increase. We offer an alternative sampling method leveraging recent insights from compressive sensing towards seismic acquisition and processing for data that, from a traditional point of view, are considered to be undersampled. The main outcome of this approach is a new technology where acquisition and processing related costs are decoupled the stringent Nyquist sampling criterion. At the heart of our approach lies randomized incoherent sampling that breaks subsampling-related interferences by turning them into harmless noise, which we subsequently remove by promoting sparsity in a transform-domain. Acquisition schemes designed to fit into this regime no longer grow significantly in cost with increasing resolution and dimensionality of the survey area, but instead its cost ideally only depends on transform-domain sparsity of the expected data. Our contribution is split into two part.},
	url = {http://csegrecorder.com/articles/view/compressive-sensing-in-seismic-exploration-outlook-on-a-new-paradigm},
	url2 = {https://slim.gatech.edu/Publications/Public/Journals/CSEGRecorder/2011/herrmann2011RECORDERcsse/herrmann2011RECORDERcsse.pdf},
	author = {Felix J. Herrmann and Haneet Wason and Tim T.Y. Lin}
}
@conference {mansour2011SBGFcspsma,
	title = {A compressive sensing perspective on simultaneous marine acquisition},
	booktitle = {SBGF},
	year = {2011},
	publisher = {SBGF},
	organization = {SBGF},
	abstract = {The high cost of acquiring seismic data in marine environments compels the adoption of simultaneous- source acquisition - an emerging technology that is stimulating both geophysical research and commercial efforts. In this paper, we discuss the properties of randomized simultaneous acquisition matrices and demonstrate that sparsity-promoting recovery improves the quality of the reconstructed seismic data volumes. Simultaneous marine acquisition calls for the development of a new set of design principles and post-processing tools. Leveraging established findings from the field of compressed sensing, the recovery from simultaneous sources depends on a sparsifying transform that compresses seismic data, is fast, and reasonably incoherent with the compressive sampling matrix. To achieve this incoherence, we use random time dithering where sequential acquisition with a single airgun is replaced by continuous acquisition with multiple airguns firing at random times and at random locations. We demonstrate our results with simulations of simultaneous Marine acquisition using periodic and randomized time dithering.},
	keywords = {Acquisition, Compressive Sensing, Presentation, SBGF},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SBGF/2011/Mansour11SBGFcspsma/Mansour11SBGFcspsma.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SBGF/2011/Mansour11SBGFcspsma/Mansour11SBGFcspsma.pdf},
	author = {Hassan Mansour and Haneet Wason and Tim T.Y. Lin and Felix J. Herrmann}
}
@presentation {wason2011SINBADode,
	title = {Only dither: efficient marine acquisition {\textquoteleft}{\textquoteleft}without" simultaneous sourcing},
	year = {2011},
	publisher = {SINBAD},
	keywords = {Presentation, SINBAD, SINBADFALL2011, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2011/Fall/wason2011SINBADode/wason2011SINBADode_pres.pdf},
	author = {Haneet Wason}
}
@conference {wason2011SEGsprsd,
	title = {Sparsity-promoting recovery from simultaneous data: a compressive sensing approach},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {30},
	year = {2011},
	month = {09},
	pages = {6-10},
	publisher = {SEG},
	organization = {SEG},
	abstract = {Seismic data acquisition forms one of the main bottlenecks in seismic imaging and inversion. The high cost of acquisition work and collection of massive data volumes compel the adoption of simultaneous-source seismic data acquisition - an emerging technology that is developing rapidly, stimulating both geophysical research and commercial efforts. Aimed at improving the performance of marine- and land-acquisition crews, simultaneous acquisition calls for development of a new set of design principles and post-processing tools. Leveraging developments from the field of compressive sensing the focus here is on simultaneous-acquisition design and sequential-source data recovery. Apart from proper compressive sensing sampling schemes, the recovery from simultaneous simulations depends on a sparsifying transform that compresses seismic data, is fast, and reasonably incoherent with the compressive-sampling matrix. Using the curvelet transform, in which seismic data can be represented parsimoniously, the recovery of the sequential-source data volumes is achieved using the sparsity-promoting program {\textemdash} SPGL1, a solver based on projected spectral gradients. The main outcome of this approach is a new technology where acquisition related costs are no longer determined by the stringent Nyquist sampling criteria.},
	keywords = {Acquisition, Compressive Sensing, Presentation, SEG},
	doi = {10.1190/1.3628174},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2011/wason11SEGsprsd/wason11SEGsprsd.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2011/wason11SEGsprsd/wason11SEGsprsd_pres.pdf},
	author = {Haneet Wason and Felix J. Herrmann and Tim T.Y. Lin}
}
@presentation {wason2010SINBADssd,
	title = {Sequential source data recovery from simultaneous acquisition through transform-domain sparsity promotion - curvelet versus shearlet transform},
	year = {2010},
	publisher = {SINBAD},
	keywords = {Presentation, SINBAD, SINBADFALL2010, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2010/Fall/wason2010SINBADssd/wason2010SINBADssd_pres.pdf},
	author = {Haneet Wason}
}
@article {herrmann2008GEOPacd,
	title = {Adaptive curvelet-domain primary-multiple separation},
	journal = {Geophysics},
	volume = {73},
	number = {3},
	year = {2008},
	month = {08},
	pages = {A17-A21},
	publisher = {SEG},
	abstract = {In many exploration areas, successful separation of primaries and multiples greatly determines the quality of seismic imaging. Despite major advances made by surface-related multiple elimination (SRME), amplitude errors in the predicted multiples remain a problem. When these errors vary for each type of multiple in different ways (as a function of offset, time, and dip), they pose a serious challenge for conventional least-squares matching and for the recently introduced separation by curvelet-domain thresholding. We propose a data-adaptive method that corrects amplitude errors, which vary smoothly as a function of location, scale (frequency band), and angle. With this method, the amplitudes can be corrected by an elementwise curvelet-domain scaling of the predicted multiples. We show that this scaling leads to successful estimation of primaries, despite amplitude, sign, timing, and phase errors in the predicted multiples. Our results on synthetic and real data show distinct improvements over conventional least-squares matching in terms of better suppression of multiple energy and high-frequency clutter and better recovery of estimated primaries. \copyright2008 Society of Exploration Geophysicists},
	keywords = {Geophysics, Processing, SLIM},
	doi = {10.1190/1.2904986},
	url = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2008/herrmann08GEOacd/herrmann08GEOacd.pdf},
	author = {Felix J. Herrmann and Deli Wang and D. J. Verschuur}
}
@article {wang2008GEOPbws,
	title = {Bayesian wavefield separation by transform-domain sparsity promotion},
	journal = {Geophysics},
	volume = {73},
	number = {5},
	year = {2008},
	month = {07},
	pages = {1-6},
	abstract = {Successful removal of coherent noise sources greatly determines the quality of seismic imaging. Major advances were made in this direction, e.g., Surface-Related Multiple Elimination (SRME) and interferometric ground-roll removal. Still, moderate phase, timing, amplitude errors and clutter in the predicted signal components can be detrimental. Adopting a Bayesian approach along with the assumption of approximate curvelet-domain independence of the to-be-separated signal components, we construct an iterative algorithm that takes the predictions produced by for example SRME as input and separates these components in a robust fashion. In addition, the proposed algorithm controls the energy mismatch between the separated and predicted components. Such a control, which was lacking in earlier curvelet-domain formulations, produces improved results for primary-multiple separation on both synthetic and real data.},
	keywords = {curvelet transform, Geophysics, Optimization, Processing, SLIM},
	doi = {10.1190/1.2952571},
	url = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2008/wang08GEObws/wang08GEObws.pdf},
	html_version = {https://slim.gatech.edu/Publications/Public/Journals/2008/wang08TRbss/paper_html/paper.html},
	author = {Deli Wang and Rayan Saab and Ozgur Yilmaz and Felix J. Herrmann}
}
@article {herrmann2008GEOPcbs,
	title = {Curvelet-based seismic data processing: a multiscale and nonlinear approach},
	journal = {Geophysics},
	volume = {73},
	number = {1},
	year = {2008},
	month = {03},
	pages = {A1-A5},
	publisher = {SEG},
	abstract = {Mitigating missing data, multiples, and erroneous migration amplitudes are key factors that determine image quality. Curvelets, little {\textquoteleft}{\textquoteleft}plane waves,{\textquoteright}{\textquoteright} complete with oscillations in one direction and smoothness in the other directions, sparsify a property we leverage explicitly with sparsity promotion. With this principle, we recover seismic data with high fidelity from a small subset (20\%) of randomly selected traces. Similarly, sparsity leads to a natural decorrelation and hence to a robust curvelet-domain primary-multiple separation for North Sea data. Finally, sparsity helps to recover migration amplitudes from noisy data. With these examples, we show that exploiting the curvelet{\textquoteright}s ability to sparsify wavefrontlike features is powerful, and our results are a clear indication of the broad applicability of this transform to exploration seismology. \copyright2008 Society of Exploration Geophysicists},
	keywords = {Acquisition, curvelet transform, Processing, SLIM},
	doi = {10.1190/1.2799517},
	url = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2008/herrmann08GEOcbs/herrmann08GEOcbs.pdf},
	author = {Felix J. Herrmann and Deli Wang and Gilles Hennenfent and Peyman P. Moghaddam}
}
@conference {herrmann2008SEGcdm,
	title = {Curvelet-domain matched filtering},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {27},
	number = {3643-3649},
	year = {2008},
	month = {08},
	publisher = {SEG},
	organization = {SEG},
	abstract = {Matching seismic wavefields and images lies at the heart of many pre-/post-processing steps part of seismic imaging{\textendash}- whether one is matching predicted wavefield components, such as multiples, to the actual to-be-separated wavefield components present in the data or whether one is aiming to restore migration amplitudes by scaling, using an image-to-remigrated-image matching procedure to calculate the scaling coefficients. The success of these wavefield matching procedures depends on our ability to (i) control possible overfitting, which may lead to accidental removal of energy or to inaccurate image-amplitude corrections, (ii) handle data or images with nonunique dips, and (iii) apply subsequent wavefield separations or migraton amplitude corrections stably. In this paper, we show that the curvelet transform allows us to address all these issues by imposing smoothness in phase space, by using their capability to handle conflicting dips, and by leveraging their ability to represent seismic data and images sparsely. This latter property renders curvelet-domain sparsity promotion an effective prior.},
	keywords = {Presentation, SEG, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2008/herrmann08SEGcdm/herrmann08SEGcdm.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2008/herrmann08SEGcdm/herrmann08SEGcdm_pres.pdf},
	author = {Felix J. Herrmann and Peyman P. Moghaddam and Deli Wang}
}
@conference {wang2008SINBADrri,
	title = {Recent results in curvelet-based primary-multiple separation},
	booktitle = {SINBAD 2008},
	year = {2008},
	abstract = {We present a nonlinear curvelet-based sparsity-promoting formulation for the primary-multiple separation problem. We show that these coherent signal components can be separated robustly by explicitly exploiting the locality of curvelets in phase space (space-spatial frequency plane) and their ability to compress data volumes that contain wavefronts. This work is an extension of earlier results and the presented algorithms are shown to be stable under noise and moderately erroneous multiple predictions.},
	keywords = {Presentation, SINBAD, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2008/wang2008SINBADrri/wang2008SINBADrri.pdf},
	author = {Deli Wang and Rayan Saab and Ozgur Yilmaz and Felix J. Herrmann}
}
@conference {herrmann2008SEGswi,
	title = {Seismic wavefield inversion with curvelet-domain sparsity promotion},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {27},
	year = {2008},
	month = {11},
	pages = {2497-2501},
	publisher = {SEG},
	organization = {SEG},
	abstract = {Inverting seismic wavefields lies at the heart of seismic data processing and imaging{\textendash}- whether one is applying {\textquoteleft}{\textquoteleft}a poor man{\textquoteright}s inverse{\textquoteright}{\textquoteright} by correlating wavefields during imaging or whether one inverts wavefields as part of a focal transform interferrometric deconvolution or as part of computing the {\textquoteright}data verse{\textquoteright}. The success of these wavefield inversions depends on the stability of the inverse with respect to data imperfections such as finite aperture, bandwidth limitation, and missing data. In this paper, we show how curvelet domain sparsity promotion can be used as a suitable prior to invert seismic wavefields. Examples include, seismic data regularization with the focused curvelet-based recovery by sparsity-promoting inversion (fCRSI), which involves the inversion of the primary-wavefield operator, the prediction of multiples by inverting the adjoint of the primary operator, and finally the inversion of the data itself {\textendash}- the so-called {\textquoteright}data inverse{\textquoteright}. In all cases, curvelet-domain sparsity leads to a stable inversion.},
	keywords = {Presentation, SEG, SLIM},
	doi = {10.1190/1.3063862},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2008/herrmann08SEGswi/herrmann08SEGswi.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2008/herrmann08SEGswi/herrmann08SEGswi_pres.pdf},
	author = {Felix J. Herrmann and Deli Wang}
}
@conference {saab2007SEGcbp,
	title = {Curvelet-based primary-multiple separation from a Bayesian perspective},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {26},
	year = {2007},
	pages = {2510-2514},
	publisher = {SEG},
	organization = {SEG},
	abstract = {In this abstract, we present a novel primary-multiple separation scheme which makes use of the sparsity of both primaries and multiples in a transform domain, such as the curvelet transform, to provide estimates of each. The proposed algorithm utilizes seismic data as well as the output of a preliminary step that provides (possibly) erroneous predictions of the multiples. The algorithm separates the signal components, i.e., the primaries and multiples, by solving an optimization problem that assumes noisy input data and can be derived from a Bayesian perspective. More precisely, the optimization problem can be arrived at via an assumption of a weighted Laplacian distribution for the primary and multiple coefficients in the transform domain and of white Gaussian noise contaminating both the seismic data and the preliminary prediction of the multiples, which both serve as input to the algorithm. \copyright2007 Society of Exploration Geophysicists},
	keywords = {Presentation, SEG, SLIM},
	doi = {10.1190/1.2792988},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2007/saab07SEGcbp/saab07SEGcbp.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2007/saab07SEGcbp/saab07SEGcbp_pres.pdf},
	author = {Rayan Saab and Deli Wang and Ozgur Yilmaz and Felix J. Herrmann}
}
@conference {verschuur2007SEGmmp,
	title = {Multiterm multiple prediction using separated reflections and diffractions combined with curvelet-based subtraction},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {26},
	year = {2007},
	pages = {2535-2539},
	publisher = {SEG},
	organization = {SEG},
	abstract = {The surface-related multiple elimination (SRME) method has proven to be successful on a large number of data cases. Most of the applications are still 2D, as the full 3D implementation is still expensive and under development. However, the earth is a 3D medium, such that 3D effects are difficult to avoid. Most of the 3D effects come from diffractive structures, whereas the specular reflections normally have less of a 3D behavior. By separating the seismic data in a specular reflecting and a diffractive part, multiple prediction can be carried out with these different subsets of the input data, resulting in several categories of predicted multiples. Because each category of predicted multiples can be subtracted from the input data with different adaptation filters, a more flexible SRME procedure is obtained. Based on some initial results from a Gulf of Mexico dataset, the potential of this approach is investigated. \copyright2007 Society of Exploration Geophysicists},
	keywords = {Presentation, SEG, SLIM},
	doi = {10.1190/1.2792993},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2007/verschuur07SEGmmp/verschuur07SEGmmp.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2007/verschuur07SEGmmp/verschuur07SEGmmp_pres.pdf},
	author = {D. J. Verschuur and Deli Wang and Felix J. Herrmann}
}
@conference {wang2007SEGrri,
	title = {Recent results in curvelet-based primary-multiple separation: application to real data},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {26},
	year = {2007},
	pages = {2500-2504},
	publisher = {SEG},
	organization = {SEG},
	abstract = {In this abstract, we present a nonlinear curvelet-based sparsity-promoting formulation for the primary-multiple separation problem. We show that these coherent signal components can be separated robustly by explicitly exploting the locality of curvelets in phase space (space-spatial frequency plane) and their ability to compress data volumes that contain wavefronts. This work is an extension of earlier results and the presented algorithms are shown to be stable under noise and moderately erroneous multiple predictions. \copyright2007 Society of Exploration Geophysicists},
	keywords = {Presentation, SEG, SLIM},
	doi = {10.1190/1.2792986},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2007/wang07SEGrri/wang07SEGrri.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2007/wang07SEGrri/wang07SEGrri_pres.pdf},
	author = {Deli Wang and Rayan Saab and Ozgur Yilmaz and Felix J. Herrmann}
}
@conference {herrmann2007SEGsdp,
	title = {Seismic data processing with curvelets: a multiscale and nonlinear approach},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {26},
	year = {2007},
	pages = {2220-2224},
	publisher = {SEG},
	organization = {SEG},
	abstract = {In this abstract, we present a nonlinear curvelet-based sparsity-promoting formulation of a seismic processing flow, consisting of the following steps: seismic data regularization and the restoration of migration amplitudes. We show that the curvelet{\textquoteright}s wavefront detection capability and invariance under the migration-demigration operator lead to a formulation that is stable under noise and missing data. \copyright2007 Society of Exploration Geophysicists},
	keywords = {Presentation, SEG, SLIM},
	doi = {10.1190/1.2792927},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2007/herrmann07SEGsdp/herrmann07SEGsdp.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2007/herrmann07SEGsdp/herrmann07SEGsdp_pres.pdf},
	author = {Felix J. Herrmann and Deli Wang and Gilles Hennenfent and Peyman P. Moghaddam}
}
@conference {herrmann2007SEGsnt,
	title = {Seismic noise: the good, the bad and the ugly},
	booktitle = {SEG Summer Research Workshop: Seismic Noise: Origins Preventions, Mitigation, Utilization},
	year = {2007},
	note = {Presented at SEG Summer Research Workshop: Seismic Noise: Origins, Prevention, Mitigation, Utilization},
	abstract = {In this paper, we present a nonlinear curvelet-based sparsity-promoting formulation for three problems related to seismic noise, namely the {\textquoteright}good{\textquoteright}, corresponding to noise generated by random sampling; the {\textquoteright}bad{\textquoteright}, corresponding to coherent noise for which (inaccurate) predictions exist and the {\textquoteright}ugly{\textquoteright} for which no predictions exist. We will show that the compressive capabilities of curvelets on seismic data and images can be used to tackle these three categories of noise-related problems.},
	keywords = {SEG, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2007/herrmann07SEGsnt/herrmann07SEGsnt.pdf},
	author = {Felix J. Herrmann and D. Wilkinson}
}
