@conference {park2025ML4SEISMICftn,
	title = {Fisher-Informed Training of Neural Operators for Reliable PDE Inversion},
	booktitle = {ML4SEISMIC Partners Meeting},
	year = {2025},
	month = {11},
	abstract = {Neural operators have shown strong performance for PDE-solution learning; however, their effectiveness in PDE-constrained optimization, such as inversion, remains limited, in part because inaccurately learned derivatives can mislead gradient-based updates. Standard training typically optimizes only the forward least-squares fit, so during iterative inversion the optimization variable can drift outside the training distribution, further degrading quality of derivative information and inversion performance. To address this limitation, we propose a training algorithm designed for inversion: the surrogate is trained not only to predict forward outputs but also to learn observation-relevant gradients, which leads to Fisher-information directed parameter updates. We evaluate forward and gradient accuracy, as well as the inversion trajectory, on problems governed by PDEs{\textemdash}Darcy flow and laminar (incompressible) Navier{\textendash}Stokes{\textemdash}which are relevant to porous-media fluid flow and energy-storage applications. Our results suggest that training for gradient fidelity with forward accuracy unlocks a pathway to reliable, efficient, neural operator{\textendash}based inversion.},
	keywords = {deep learning, derivatve-informed, Fisher-information, FNO, GCS, Inverse problems, ML4SEISMIC, permeability, reduced-order modeling, reservoir simulation, SLIM, Uncertainty quantification},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ML4SEISMIC/2025/park2025ML4SEISMICftn},
	author = {Jeongjin Park and Grant Bruer and Huseyin Tuna Erdinc and Richard Rex Arockiasamy and Nisha Chandramoorthy and Felix J. Herrmann}
}
@conference {zeng2025IMAGEfwi,
	title = {Full-waveform variational inference with common-image gathers},
	booktitle = {International Meeting for Applied Geoscience and Energy},
	year = {2025},
	month = {8},
	keywords = {Bayesian inference, CIG, deep learning, diffusion models, FWI, IMAGE, Imaging, Inverse problems, MVA, RTM, SEG, Summary Statistics, Uncertainty quantification, WISE},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2025/zeng2025IMAGEfwi/IMAGE_Yunlin.pdf},
	author = {Yunlin Zeng and Huseyin Tuna Erdinc and Rafael Orozco and Felix J. Herrmann}
}
@presentation {rex2024ML4SEISMICfca,
	title = {FNO-charged ASPIRE},
	year = {2024},
	month = {11},
	abstract = {During this talk, we will demonstrate how extended re-migrations{\textemdash}i.e, formation of subsurface-offset Common-Image Gathers (CIGs) for a new velocity model, can be avoided altogether by training Fourier Neural Operators during training of ASPIRE {\textemdash} Amortized posteriors with Summaries that are Physics-based and Iteratively REfined. In this approach, FNOs are trained as surrogates capable of mapping CIGs for one migration-velocity model to the other. The approach is computationally feasible because it uses the same training set as used during ASPIRE. As a result, additional training costs are small and the inference costs are reduced by a factor equal to the number of ASPIRE refinements.},
	keywords = {deep learning, FNO, hierarchical tucker tensor, Inverse problems, kronecker product, ML4SEISMIC, SLIM, two-phase flow},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ML4SEISMIC/2024/rex2024ML4SEISMICfca},
	author = {Richard Rex and Yunlin Zeng and Ziyi Yin and Rafael Orozco and Felix J. Herrmann}
}
@conference {herrmann2023IMAGEWSfnf,
	title = {Fast neural FWI with amortized uncertainty quantification},
	booktitle = {International Meeting for Applied Geoscience and Energy},
	year = {2023},
	note = {(IMAGE workshop, Houston)},
	month = {08},
	keywords = {bayesian, experimental design, FWI, Normalizing flows, SEG, Uncertainty quantification},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2023/herrmann2023IMAGEWSfnf},
	author = {Rafael Orozco and Mathias Louboutin 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}
}
@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}
}
@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 {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 {kadu2018SEGMSFWI,
	title = {Full-waveform inversion with Mumford-Shah regularization},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2018},
	note = {(SEG, Anaheim)},
	month = {10},
	pages = {1258-1262},
	abstract = {Full-waveform inversion (FWI) is a non-linear procedure to es- timate subsurface rock parameters from surface measurements of induced seismic waves. This procedure is ill-posed in nature and hence, requires regularization to enhance some structure depending on the prior information. Recently, Total-variation (TV) regularization has gained popularity due to its ability to produce blocky structures. Contrary to this, the earth behave more like a piecewise smooth function. TV regularization fails to enforce this prior information into FWI. We propose a Mumford-Shah functional to incorporate the piecewise smooth spatial structure in the FWI procedure. The resulting opti- mization problem is solved by a splitting method. We show the improvement in results against TV regularization on two synthetic camembert examples.},
	keywords = {FWI, Mumford-Shah, regularization, SEG},
	doi = {10.1190/segam2018-2997224.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2018/kadu2018SEGMSFWI/kadu2018SEGMSFWI.pdf},
	author = {Ajinkya Kadu and Rajiv Kumar and Tristan van Leeuwen}
}
@conference {yang2017SEGfsp,
	title = {Fast sparsity-promoting least-squares migration with multiples in time domain},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2017},
	note = {(SEG, Houston)},
	month = {09},
	pages = {4828-4832},
	abstract = {Based on the latest developments of research in inversion technology with optimization, researchers have made significant progress in the implementation of least-squares reverse-time migration (LS-RTM) of primaries. In Marine data however, these applications rely on the success of a pre-imaging separation of primaries and multiples, which can be modeled as a multi-dimensional convolution between the vertical derivative of the surface-free Green{\textquoteright}s function and the down-going receiver wavefield. Instead of imaging the primaries and multiples separately, we implement the LS-RTM of the total down-going wavefield by combining areal source injection and linearized Born modelling, where strong surface related multiples are generated from a strong density variation at the ocean bottom. The advantage including surface related multiples in LS-RTM is the extra illumination we obtain from these multiples without incurring additional computational costs related to carrying out multi-dimensional convolutions part of conventional multiple prediction procedures. Even though we are able to avert these computational costs, our approach shares the large costs of LSRTM. We reduce these costs by combining randomized source subsampling with our sparsity-promoting imaging technology, which produces artifact-free, high-resolution images, with the surface-related multiples migrated properly.},
	keywords = {least-squares RTM, multiples, SEG, sparsity, time domain},
	doi = {10.1190/segam2017-17741608.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2017/yang2017SEGfsp/yang2017SEGfsp.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2017/yang2017SEGfsp/yang2017SEGfsp_pres.pdf},
	author = {Mengmeng Yang and Emmanouil Daskalakis and Felix J. Herrmann}
}
@presentation {kumar2017SINBADFfas,
	title = {Full-azimuth seismic data processing w/ coil acquisition},
	year = {2017},
	publisher = {SINBAD},
	abstract = {In this work, we will demonstrate the performance of our in-house 5D low-rank based interpolation method on a seismic data acquired using coil shooting full-azimuth acquisition. We will show that we can recover full-azimuthal interpolated data from highly subsampled data, where the subsampling ratio is 4\%. This is the first time, we are testing our interpolation ideas on real 3D marine seismic data acquisition. Our findings show that we can avoid the general practice of windowing the data while performing the interpolation, specially using rank-minimization based framework.},
	keywords = {Presentation, SINBAD, SINBADFALL2017, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/kumar2017SINBADFfas/kumar2017SINBADFfas.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/kumar2017SINBADFfas/kumar2017SINBADFfas.mov},
	author = {Rajiv Kumar and Nick Moldoveanu and Keegan Lensink 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 {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}
}
@mastersthesis {tu2015THfis,
	title = {Fast imaging with surface-related multiples},
	year = {2015},
	note = {(PhD)},
	month = {08},
	school = {The University of British Columbia},
	type = {phd},
	address = {Vancouver},
	abstract = {Surface-related multiples, which are waves that bounce more than once between the water surface and the subsurface reflectors, constitute a significant part of the data acquired in marine seismic surveys. If left untreated, they can lead to misplaced phantom reflectors in the image, and result in erroneous interpretations of the subsurface structure. As a result, these multiples are removed before the imaging procedure in conventional seismic data processing. However, because they interact more with the subsurface medium, they may carry extra information that is not present in the primaries. Therefore instead of removing these multiples, a more desirable alternative is to make active use of them. We derive from the well-established "Surface-Related Multiple Elimination" relation, and arrive at a linearized expression of the wave-equation based modelling that incorporates the surface- related multiples. We then present a computationally efficient approach to iteratively invert this expression to obtain an image of the subsurface from data that contain multiples. We achieve the computational efficiency inside each iteration by (i) using the wave-equation solver to implicitly carry out the expensive multiple prediction; and (ii) reducing the number of wave-equation solves during data simulation by subsampling the monochromatic source experiments. We show that, compared with directly applying the cross-correlation/deconvolutional imaging conditions, the presented approach can suppress the coherent imaging artifacts from multiples more effectively. We also show that, by curvelet-domain sparsity promoting and occasionally drawing new data samples during the inversion, the proposed inversion method gains improved robustness to velocity errors in the background model, as well as modelling errors incurred during linearization of the wave-equation. To combine the information encoded in both the primaries and the multiples, we then propose a highly accurate source estimation method to jointly invert the total upgoing wavefield. We show with field data examples that we can reap benefits from both the relative noise-free primaries and the extra illumination coverage of the multiples. We also demonstrate that the inclusion of multiples help mitigate the amplitude ambiguity during source estimation. We conclude the thesis with an outlook for future research directions, as well as potential extensions of the proposed work.},
	keywords = {inversion, least-squares, multiples, seismic imaging},
	url = {https://slim.gatech.edu/Publications/Public/Thesis/2015/tu2015THfis/tu2015THfis.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Thesis/2015/tu2015THfis/tu2015THfis_pres.pdf},
	author = {Ning Tu}
}
@article {tu2014fis,
	title = {Fast imaging with surface-related multiples by sparse inversion},
	journal = {Geophysical Journal International},
	volume = {201},
	number = {1},
	year = {2015},
	note = {(Geophysical Journal International)},
	month = {04},
	pages = {304-317},
	abstract = {In marine exploration seismology, surface-related multiples are usually treated as noise mainly because subsequent processing steps, such as migration velocity analysis and imaging, require multiple-free data. Failure to remove these wavefield components from the data may lead to erroneous estimates for migration velocity or result in strong coherent artefacts that interfere with the imaged reflectors. However, multiples can carry complementary information compared to primaries, as they interact with the free surface and are therefore exposed more to the subsurface. Recent work has shown that when processed correctly multiples can improve seismic illumination. Given a sufficiently accurate background velocity model and an estimate for the source signature, we propose a new and computationally efficient linearized inversion procedure based on two-way wave equations, which produces accurate images of the subsurface from the total upgoing wavefield including surface-related multiples. Modelling of the surface-related multiples in the proposed method derives from the well-known surface-related multiple elimination method. We incur a minimal overhead from incorporating the multiples by having the wave-equation solver carry out the multiple predictions via the inclusion of an areal source instead of expensive dense matrix-matrix multiplications. By using subsampling techniques, we obtain high-quality true-amplitude least-squares migrated images at computational costs of roughly a single reverse-time migration (RTM) with all the data. These images are virtually free of coherent artefacts from multiples. Proper inversion of the multiples would be computationally infeasible without using these techniques that significantly bring down the cost. By promoting sparsity in the curvelet domain and using rerandomization, out method gains improved robustness to errors in the background velocity model, and errors incurred in the linearization of the wave equation with respect to the model. We demonstrate the superior performance of the proposed method compared to the conventional RTM using realistic synthetic examples.},
	keywords = {approximate message passing, Compressive Sensing, curvelet, inversion, Kaczmarz, multiples},
	url = {https://slim.gatech.edu/Publications/Public/Journals/GeophysicalJournalInternational/2014/tu2014fis/tu2014fis.pdf},
	url2 = {http://gji.oxfordjournals.org/cgi/content/full/ggv020?ijkey=pKzrSXbtxjlOjzd\&keytype=ref},
	author = {Ning Tu and Felix J. Herrmann}
}
@presentation {yang2015SINBADFfis,
	title = {Fast imaging with surface-related multiples: {\textendash} shallow water multiples},
	year = {2015},
	publisher = {SINBAD},
	abstract = {In this talk, we will show surface-related multiples imaging with source estimation. Here the algorithms are based on linearized bregman, and we give the source estimation formulation in this case. With the featured blocked structure of linearized bregman, we can do efficient imaging in straightforward way by randomly choosing sources and frequencies. Another benefit is the simplicity to delivery the algorithm.},
	keywords = {Presentation, SINBAD, SINBADFALL2015, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/yang2015SINBADFfis/yang2015SINBADFfis.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/yang2015SINBADFfis/yang2015SINBADFfis.mov},
	author = {Mengmeng Yang and Ning Tu and Felix J. Herrmann}
}
@article {tu2015TLEfls,
	title = {Fast least-squares imaging with surface-related multiples: application to a North-Sea data set},
	journal = {The Leading Edge},
	volume = {34},
	number = {7},
	year = {2015},
	note = {(The Leading Edge)},
	month = {07},
	pages = {788{\textendash}794},
	abstract = {In marine seismic acquisition, surface-related multiples constitute a significant portion of the acquired data. Typically, multiples are removed during early-stage data-processing as they can lead to phantom reflectors during migration that may result in erroneous geological interpretations. However, if properly dealt with, multiples can provide valuable extra information and complement primaries in illuminating the subsurface. In this article, we demonstrate the limitation of the reverse-time migration in imaging these multiples, and present an alternative inversion procedure that is computationally efficient, that jointly maps both primaries and multiples to the true reflectors, and where the source function is estimated on the fly. As a result, we obtain high-quality, mostly artifact-free, broad-band images where the imprint of the source-function are partly removed at a computationally affordable expense compared to the combined costs of the wave-equation based surface-related multiple elimination and the reverse-time migration. We achieve all this by including the total downgoing wavefields as areal sources in the least-squares migration in combination with curvelet-domain sparsity promotion. We apply the proposed method to a shallow-water marine data set from the North sea, which contains abundant short-period surface-related multiples, and show its efficacy in eliminating coherent imaging artifacts associated with these multiples. We also demonstrate the benefits of joint imaging of primaries and multiples compared to imaging these signal components separately.},
	keywords = {field data, inversion, least squares, migration, multiples},
	doi = {10.1190/tle34070788.1},
	url = {https://slim.gatech.edu/Publications/Public/Journals/TheLeadingEdge/2015/tu2015TLEfls/tu2015TLEfls.html},
	author = {Ning Tu and Felix J. Herrmann}
}
@presentation {herrmann2015SINBADfom,
	title = {Fast {\textquotedblleft}online{\textquotedblright} migration with Compressive Sensing},
	year = {2015},
	publisher = {SINBAD},
	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 = {Presentation, SINBAD, SINBADSPRING2015, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Spring/herrmann2015SINBADfom/herrmann2015SINBADfom.pdf},
	author = {Felix J. Herrmann}
}
@conference {herrmann2015EAGEfom,
	title = {Fast "online" migration with Compressive Sensing},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2015},
	note = {(EAGE, Madrid)},
	month = {06},
	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 = {EAGE, LSRTM},
	doi = {10.3997/2214-4609.201412942},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2015/herrmann2015EAGEfom/herrmann2015EAGEfom.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2015/herrmann2015EAGEfom/herrmann2015EAGEfom_pres.pdf},
	author = {Felix J. Herrmann and Ning Tu and Ernie Esser}
}
@presentation {tu2014SINBADfis2,
	title = {Fast imaging with source estimation},
	year = {2014},
	publisher = {SINBAD},
	abstract = {In seismic imaging, an inaccurate estimate of the source wavelet may result in degraded seismic images. While conventional reverse-time migration requires the knowledge of the source wavelet as prior information, we propose to include source estimation in the fast compressive imaging procedure. Using the proposed method, we can obtain seismic images that are comparable to those imaged with the true source wavelet, with computational costs that are comparable to conventional RTM images with all the data. We also extend the proposed method to image data with surface-related multiples, where these multiples help to mitigate the amplitude ambiguity in source estimation. We verify the proposed method using realistic synthetic examples.},
	keywords = {Presentation, SINBAD, SINBADFALL2014, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2014/Fall/tu2014SINBADfis2/tu2014SINBADfis2.pdf},
	author = {Ning Tu and Aleksandr Y. Aravkin and Tristan van Leeuwen and Tim T.Y. Lin and Felix J. Herrmann}
}
@presentation {tu2014SINBADfis,
	title = {Fast imaging with surface-related multiples by sparse inversion},
	year = {2014},
	publisher = {SINBAD},
	abstract = {When used correctly, surface-related multiples can provide extra illumination coverage compared with primaries. In this talk, I will discuss how to jointly image primaries and surface-related multiples in a computationally efficient fashion. We bring down the computational cost by two means. First we use wave-equation solvers to implicitly carry out the expensive dense matrix-matrix multiplications in the prediction of surface-related multiples. Second we bring down the simulation cost by subsampling the frequencies and monochromatic source experiments together with curvelet-domain sparsity-promoting and rerandomization. As a result, we obtain true-amplitude least-squares migrated seismic images with computational costs that are comparable to a single RTM with all the data. We demonstrate the efficacy of the proposed method using realistic synthetic examples.},
	keywords = {Presentation, SINBAD, SINBADFALL2014, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2014/Fall/tu2014SINBADfis/tu2014SINBADfis.pdf},
	author = {Ning Tu and Felix J. Herrmann}
}
@conference {herrmann2014SEGWSfli,
	title = {Fast linearized inversion with surface-related multiples with source estimation},
	booktitle = {SEG Workshop on Using Multiples as Signal for Imaging; Denver},
	year = {2014},
	note = {(SEG Workshop, Denver)},
	month = {10},
	abstract = {In the well-known SRME relation, multiples are expressed as the multi-dimensional convolution between the subsurface Green{\textquoteright}s function and the downgoing receiver wavefield. Therefore these multiples can be considered as the response of the subsurface to an "areal" source term given by the same downgoing receiver wavefield. This relation can be used to treat these multiples as signals instead of considering them as noise that must removed before seismic imaging. However, when we use conventional reverse-time migration to image these multiple events, it results in acausal imaging artifacts caused by cross-correlations between wrong pairs of up- and downgoing wavefields. We find that these artifacts can be removed by adopting a sparsity-promoting least-squares inversion approach. However, iterative inversions go at the expense of excessive computational costs. By combining the SRME relation and wave-equation based linearized modelling, we are able to significantly reduce the costs by avoiding the dense matrix-matrix multiplications of SRME and the number of wave-equation solves via source subsampling. As a result, we arrive at a method with a cost comparable to that of a single reverse-time migration with all shots.},
	keywords = {inversion, multiples, SEG, workshop},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2014/herrmann2014SEGWSfli/herrmann2014SEGWSfli.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2014/herrmann2014SEGWSfli/herrmann2014SEGWSfli_pres.pdf},
	author = {Ning Tu and Felix J. Herrmann}
}
@article {aravkin2014SISCfmd,
	title = {Fast methods for denoising matrix completion formulations, with applications to robust seismic data interpolation},
	journal = {SIAM Journal on Scientific Computing},
	volume = {36},
	number = {5},
	year = {2014},
	note = {(SISC)},
	month = {10},
	pages = {S237-S266},
	abstract = {Recent SVD-free matrix factorization formulations have enabled rank minimization for systems with millions of rows and columns, paving the way for matrix completion in extremely large-scale applications, such as seismic data interpolation. In this paper, we consider matrix completion formulations designed to hit a target data-fitting error level provided by the user, and we propose an algorithm called LR-BPDN that is able to exploit factorized formulations to solve the corresponding optimization problem. Since practitioners typically have strong prior knowledge about target error level, this innovation makes it easy to apply the algorithm in practice, leaving only the factor rank to be determined. Within the established framework, we propose two extensions that are highly relevant to solving practical challenges of data interpolation. First, we propose a weighted extension that allows known subspace information to improve the results of matrix completion formulations. We show how this weighting can be used in the context of frequency continuation, an essential aspect to seismic data interpolation. Second, we propose matrix completion formulations that are robust to large measurement errors in the available data. We illustrate the advantages of LR-BPDN on collaborative filtering problems using the MovieLens 1M and 10M and Netflix 100M datasets. Then we use the new method, along with its robust and subspace reweighted extensions, to obtain high-quality reconstructions for large-scale seismic interpolation problems with real data, even in the presence of data contamination.},
	keywords = {factorized formulations, matrix completion, robust statistics, seismic trace interpolation, variational optimization},
	doi = {10.1137/130919210},
	url = {https://slim.gatech.edu/Publications/Public/Journals/SIAMJournalOnScientificComputing/2014/aravkin2014SISCfmd/aravkin2014SISCfmd.pdf},
	url2 = {http://epubs.siam.org/doi/abs/10.1137/130919210},
	author = {Aleksandr Y. Aravkin and Rajiv Kumar and Hassan Mansour and Ben Recht and Felix J. Herrmann}
}
@conference {lago2014EAGEfst,
	title = {Fast solution of time-harmonic wave-equation for full-waveform inversion},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2014},
	month = {06},
	abstract = {For many full-waveform inversion techniques, the most computationally intensive step is the computation of a numerical solution for the wave equation on every iteration. In the frequency domain approach, this requires the solution of very large, complex, sparse, ill-conditioned linear systems. In this extended abstract we bring out attention specifically to CGMN method for solving PDEs, known for being flexible (i.e. it is able to treat equally acoustic data as well as visco-elastic or more complex scenarios) efficient with respect both to memory and computation time, and controllable accuracy of the final approximation. We propose an improvement for the known CGMN method by imposing a minimal residual condition, which incurs in one extra model vector storage. The resulting algorithm called CRMN enjoys several interesting properties as monotonically nonincreasing behaviour of the norm of the residual and minimal residual, guaranteeing optimal convergence for the relative residual criterion. We discuss numerical experiments both in an isolated PDE solve and also within the inversion procedure, showing that in a realistic scenario we can expect a speedup around 25\% when using CRMN rather than CGMN.},
	keywords = {CGMN, CRMN, EAGE, FWI, Time-Harmonic Wave Equation},
	doi = {10.3997/2214-4609.20140812},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2014/lago2014EAGEfst/lago2014EAGEfst.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2014/lago2014EAGEfst/lago2014EAGEfst_pres.pdf},
	author = {Rafael Lago and Art Petrenko and Zhilong Fang and Felix J. Herrmann}
}
@conference {fang2014EAGEfuq,
	title = {Fast uncertainty quantification for 2D full-waveform inversion with randomized source subsampling},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2014},
	month = {06},
	abstract = {Uncertainties arise in every area of seismic exploration, especially in full-waveform inversion, which is highly non-linear. In the framework of Bayesian inference, uncertainties can be analyzed by sampling the posterior probability density distribution with a Markov chain Monte-Carlo (McMC) method. We reduce the cost of computing the posterior distribution by working with randomized subsets of sources. These approximations, together with the Gaussian assumption and approximation of the Hessian, leads to a computational tractable uncertainty quantification. Application of this approach to a synthetic leads to standard deviations and confidence intervals that are qualitatively consistent with our expectations.},
	keywords = {EAGE, FWI, Markov chain Monte Carlo, randomized sources subsampling, Uncertainty quantification},
	doi = {10.3997/2214-4609.20140715},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2014/fang2014EAGEfuq/fang2014EAGEfuq.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2014/fang2014EAGEfuq/fang2014EAGEfuq_pres.pdf},
	author = {Zhilong Fang and Curt Da Silva 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 {miao2013SINBADfid,
	title = {Fast imaging via depth stepping with the two-way wave equation},
	year = {2013},
	publisher = {SINBAD},
	abstract = {In this presentation we propose a fast imaging algorithm via depth stepping with the two-way wave equation. Within the framework of survey sinking, a stabilized depth extrapolation operator is computed using a spectral projector which can efficiently split evanescent wave components. The computation of the spectral projector features with an Hierarchically Semi-Seperable (HSS) matrix representation speeded up polynomial recursion, resulting in an accleration from cubic numerical complexity to linear numerical complexity.},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/miao2013SINBADfid/miao2013SINBADfid.pdf},
	author = {Lina Miao and Felix J. Herrmann}
}
@presentation {tu2013SINBADvp,
	title = {Fast imaging with multiples and source estimation},
	year = {2013},
	publisher = {SINBAD},
	abstract = {Multiples are usually treated as unwanted components in seismic data. However, if correctly used, they can provide valuable information about the subsurface. In this presentation, I will talk about how to make use of multiples in $\ell_1$ regularized least-squares imaging, and how to estimate the source wavelet on the fly using multiples. Synthetic examples show that using multiples not only helps to retrieve the true-amplitude seismic image and the source wavelet, it also increases the wavenumber contents in the seismic image. By using dimensionality reduction techniques with rerandomization, we also greatly decrease the simulation cost without compromising the image quality.},
	keywords = {Presentation, SINBAD, SINBADSPRING2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Spring/tu2013SINBADvp/tu2013SINBADvp.pdf},
	author = {Ning Tu and Aleksandr Y. Aravkin and Tristan van Leeuwen and Felix J. Herrmann}
}
@presentation {tu2013SINBADfim,
	title = {Fast imaging with multiples and source estimation},
	year = {2013},
	publisher = {SINBAD},
	abstract = {During this talk, we present a computationally efficient (cost of 1-2 RTMs with all data) iterative sparsity-promoting inversion framework where surface-related multiples are jointly imaged with primaries and where the source signature is estimated on the fly. Our imaging algorithm is computationally efficient because it works during each iteration with small independent randomized subsets of data. The multiples are handled by introducing an areal source term that includes the upgoing wavefield. We update the source signature for each iteration using a variable projection method. The resulting algorithm removes imaging artifacts from surface-related multiples, estimates and removes the imprint of the source, recovers true amplitudes, is fast, and robust to linearization errors by virtue of the statistical independence of the subsets of data we are working with at each iteration.},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/tu2013SINBADfim/tu2013SINBADfim.pdf},
	author = {Ning Tu and Felix J. Herrmann}
}
@conference {tu2013EAGElsm,
	title = {Fast least-squares migration with multiples and source estimation},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2013},
	month = {06},
	abstract = {The advent of modern computing has made it possible to do seismic imaging using least-squares reverse-time migration. We obtain superior images by solving an optimization problem that recovers the true-amplitude images. However, its success hinges on overcoming several issues, including overwhelming problem size, unknown source wavelet, and interfering coherent events like multiples. In this abstract, we reduce the problem size by using ideas from compressive sensing, and estimate source wavelet by generalized variable projection. We also demonstrate how to invert for subsurface information encoded in surface-related multiples by incorporating the free-surface operator as an areal source in reverse-time migration. Our synthetic examples show that multiples help to improve the resolution of the image, as well as remove the amplitude ambiguity in wavelet estimation.},
	keywords = {EAGE, Imaging, multiples, source estimation, sparse},
	doi = {10.3997/2214-4609.20130727},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2013/tu2013EAGElsm/tu2013EAGElsm.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2013/tu2013EAGElsm/tu2013EAGElsm_pres.pdf},
	author = {Ning Tu and Aleksandr Y. Aravkin and Tristan van Leeuwen and Felix J. Herrmann}
}
@conference {herrmann2013EAGEfrtm,
	title = {Fast RTM with multiples and source estimation},
	booktitle = {EAGE/SEG Forum - Turning noise into geological information: The next big step?},
	year = {2013},
	month = {11},
	abstract = {During this talk, we present a computationally efficient (cost of 1-2 RTM{\textquoteright}s with all data) iterative sparsity-promoting inversion framework where surface-related multiples are jointly imaged with primaries and where the source signature is estimated on the fly. Our imaging algorithm is computationally efficient because it works during each iteration with small independent randomized subsets of data. The multiples are handled by introducing an areal source term that includes the upgoing wavefield. We update the source signature for each iteration using a variable projection method. The resulting algorithm removes imaging artifacts from surface-related multiples, estimates and removes the imprint of the source, recovers true amplitudes, is fast, and robust to linearization errors by virtue of the statistical independence of the subsets of data we are working with at each iteration.},
	keywords = {EAGE, multiples, RTM, SEG, source estimation},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2013/herrmann2013EAGEfrtm/herrmann2013EAGEfrtm_pres.pdf},
	author = {Felix J. Herrmann and Ning Tu}
}
@article {VanLeeuwen11TRfwiwse,
	title = {Fast waveform inversion without source encoding},
	journal = {Geophysical Prospecting},
	volume = {61},
	year = {2013},
	note = {Article first published online: 10 JULY 2012},
	month = {06},
	pages = {10-19},
	abstract = {Randomized source encoding has recently been proposed as a way to dramatically reduce the costs of full waveform inversion. The main idea is to replace all sequential sources by a small number of simultaneous sources. This introduces random crosstalk in the model updates and special stochastic optimization strategies are required to deal with this. Two problems arise with this approach: i) source encoding can only be applied to fixed-spread acquisition setups, and ii) stochastic optimization methods tend to converge very slowly, relying on averaging to get rid of the cross-talk. Although the slow convergence is partly offset by the low iteration cost, we show that conventional optimization strategies are bound to outperform stochastic methods in the long run. In this paper we argue that we don{\textlnot}{\o}t need randomized source encoding to reap the benefits of stochastic optimization and we review an optimization strategy that combines the benefits of both conventional and stochastic optimization. The method uses a gradually increasing batch of sources. Thus, iterations are very cheap initially and this allows the method to make fast progress in the beginning. As the batch size grows, the method behaves like conventional optimization, allowing for fast convergence. Numerical examples suggest that the stochastic and hybrid method perform equally well with and without source encoding and that the hybrid method outperforms both conventional and stochastic optimization. The method does not rely on source encoding techniques and can thus be applied to non fixed-spread data.},
	keywords = {FWI, Optimization, SLIM},
	doi = {10.1111/j.1365-2478.2012.01096.x},
	url = {http://onlinelibrary.wiley.com/doi/10.1111/j.1365-2478.2012.01096.x/abstract},
	url2 = {https://slim.gatech.edu/Publications/Public/Journals/GeophysicalProspecting/2013/VanLeeuwen11TRfwiwse/VanLeeuwen11TRfwiwse.pdf},
	author = {Tristan van Leeuwen 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 {herrmann2013SINBADffwi,
	title = {Frugal FWI},
	year = {2013},
	publisher = {SINBAD},
	abstract = {Seismic waveform inversion aims at obtaining detailed estimates of subsurface medium parameters, such as soundspeed, from seismic data. A formulation in the frequency-domain leads to an optimization problem constrained by a Helmholtz equation with many right-hand sides. Application of this technique in 3D precludes the use of factorization techniques to solve the Helmholtz equation due the large number of gridpoints and the bandwidth of the matrix. While many sophisticated pre-conditioned iterative techniques have been developed for the Helmholtz equation, they often include model-specific tuning parameters and are thus not very attractive for inversion since the medium parameters change from one iteration to the next. In this paper, we propose a method for 3D seismic waveform inversion that addresses both the need to efficiently solve the Helmholtz equation as well as the computational cost induced by the many right-hand sides. To solve the Helmholtz equation, we consider a simple generic preconditioned iterative method (CARP-CG) that is well-suited for inversion because of its robustness. We extend this method to a block-iterative method that can efficiently handle multiple right-hand sides. To reduce the computational cost of of the overall optimization procedure, we use recently proposed techniques from stochastic optimization that allow us to work with approximate gradient information. These approximations are obtained by evaluating only a small portion of the right-hand sides and/or by solving the PDE approximately. We propose heuristics to adaptively determine the required accuracy of the PDE solves and the sample-size and illustrate the algorithms on synthetic benchmark models. This is joint work with Tristan van Leeuwen.},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/herrmann2013SINBADffwi/herrmann2013SINBADffwi.pdf},
	author = {Felix J. Herrmann}
}
@presentation {herrmann2013SINBADfwi,
	title = {Frugal FWI},
	year = {2013},
	publisher = {SINBAD},
	abstract = {Seismic waveform inversion aims at obtaining detailed estimates of subsurface medium parameters, such as soundspeed, from seismic data. A formulation in the frequency-domain leads to an optimization problem constrained by a Helmholtz equation with many right-hand- sides. Application of this technique in 3D precludes the use of factorization techniques to solve the Helmholtz equation due the large number of gridpoints and the bandwidth of the matrix. While many sophisticated pre-conditioned iterative techniques have been developed for the Helmholtz equation, they often include model-specific tuning parameters and are thus not very attractive for inversion since the medium parameters change from one iteration to the next. In this paper, we propose a method for 3D seismic waveform inversion that addresses both the need to efficiently solve the Helmholtz equation as well as the computational cost induced by the many right-hand-sides. To solve the Helmholtz equation, we consider a simple generic preconditioned iterative method (CARP-CG) that is well-suited for inversion because of its robustness. We extend this method to a block-iterative method that can efficiently handle multiple right-hand sides. To reduce the computational cost of of the overall optimization procedure, we use recently proposed techniques from stochastic optimization that allow us to work with approximate gradient information. These approximations are obtained by evaluating only a small portion of the right-hand sides and/or by solving the PDE approximately. We propose heuristics to adaptively determine the required accuracy of the PDE solves and the sample-size and illustrate the algorithms on synthetic benchmark models. This is joint work with Tristan van Leeuwen.},
	keywords = {Presentation, SINBAD, SINBADSPRING2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Spring/herrmann2013SINBADfwi/herrmann2013SINBADfwi.pdf},
	author = {Felix J. Herrmann}
}
@presentation {xiang2012SINBADfgn,
	title = {Fast Gauss-Newton full-waveform inversion with sparsity regularization},
	year = {2012},
	publisher = {SINBAD},
	abstract = {Full-waveform inversion (FWI) can be considered as a controlled data fitting process, in which we approximately fit observed data by iteratively updating the initial velocity model, we expect the final model can reveal subsurface structure till the wavefield misfit can converge to designed tolerance. The conventional FWI approach is expensive since it requires the inversion of a linear system, which involves extremely large multi-experiment data volumes. To overcome this issue we percent a curvetlet based sparsity-promoting Gauss-Newton inversion method. In this presentation we invert for the model updates by replacing the normal Gauss-Newton linearized subproblem for subsampled FWI with a sparsity promoting FWI formulation. We speed up the algorithm and avoid over fitting the data by solving the problem approximately. Aside from this, we control wavefield dispersion by gradually increasing grid size as we move to higher frequencies. Our approach is successful because it reduces the size of seismic data volumes without loss of information. With this reduction, we can compute a Newton-like update with the reduced data volume at the cost of roughly one gradient update for the fully sampled wavefield.},
	keywords = {Presentation, SINBAD, SINBADFALL2012, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2012/Fall/xiang2012SINBADfgn/xiang2012SINBADfgn_pres.pdf},
	author = {Xiang Li}
}
@presentation {ning2012SINBADfim,
	title = {Fast imaging with multiples},
	year = {2012},
	publisher = {SINBAD},
	abstract = {If correctly used, multiple energy can be mapped to the correct subsurface locations. However, simply applying the cross-correlation imaging condition will introduce non-causal artifacts into the final image. Here we propose an inversion approach to image primaries and multiples simultaneously that yields an artifact-free image. To address the high computational cost associated with inversion, we propose to: i) have the wave-equation solver carry out the multi-dimensional convolutions implicitly, and ii) reduce the number of PDE solves by randomized subsampling. We then propose to improve the overall performance of this algorithm by a process called rerandomization, which helps to cancel the correlation built between model iterate and the subsampling operator. We show the merits of our approach on a number of examples.},
	keywords = {Presentation, SINBAD, SINBADFALL2012, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2012/Fall/ning2012SINBADfim/ning2012SINBADfim_pres.pdf},
	author = {Ning Tu}
}
@presentation {ning2012SINBADSPRINGfim,
	title = {Fast imaging with multiples by sparse inversion},
	year = {2012},
	publisher = {SINBAD},
	keywords = {Presentation, SINBAD, SINBADSPRING2012, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2012/Spring/ning2012SINBADSPRINGfim/ning2012SINBADSPRINGfim_pres.pdf},
	author = {Ning Tu}
}
@techreport {rajiv2012SEGFRM,
	title = {Fast methods for rank minimization with applications to seismic-data interpolation},
	number = {TR-EOAS-2012-3},
	year = {2012},
	month = {04},
	institution = {Department of Earth and Ocean Sciences},
	address = {University of British Columbia, Vancouver},
	abstract = {Rank penalizing techniques are an important direction in seismic inverse problems, since they allow improved recovery by exploiting low-rank structure. A major downside of current state of the art techniques is their reliance on the SVD of seismic data structures, which can be prohibitively expensive. Fortunately, recent work allows us to circumvent this problem by working with matrix factorizations. We review a novel approach to rank penalization, and successfully apply it to the seismic interpolation problem by exploiting the low-rank structure of seismic data in the midpoint-offset domain. Experiments for the recovery of 2D monochromatic data matrices and seismic lines represented as 3D volumes support the feasibility and potential of the new approach.},
	keywords = {Optimization, Rank, Seismic data Interpolation},
	url = {https://slim.gatech.edu/Publications/Public/TechReport/2012/rajiv2012SEGFRM/rajiv2012SEGFRM.pdf},
	author = {Rajiv Kumar and Aleksandr Y. Aravkin and Felix J. Herrmann}
}
@article {Li11TRfrfwi,
	title = {Fast randomized full-waveform inversion with compressive sensing},
	journal = {Geophysics},
	volume = {77},
	number = {3},
	year = {2012},
	month = {05},
	pages = {A13-A17},
	abstract = {Wave-equation based seismic inversion can be formulated as a nonlinear inverse problem where the medium properties are obtained via minimization of a least- squares misfit functional. The demand for higher resolution models in more geologically complex areas drives the need to develop techniques that explore the special structure of full-waveform inversion to reduce the computational burden and to regularize the inverse problem. We meet these goals by using ideas from compressive sensing and stochastic optimization to design a novel Gauss-Newton method, where the updates are computed from random subsets of the data via curvelet-domain sparsity promotion. Application of this idea to a realistic synthetic shows improved results compared to quasi-Newton methods, which require passes through all data. Two different subset sampling strategies are considered: randomized source encoding, and drawing sequential shots firing at random source locations from marine data with missing near and far offsets. In both cases, we obtain excellent inversion results compared to conventional methods at reduced computational costs.},
	keywords = {Compressive Sensing, FWI, Optimization, SLIM},
	doi = {10.1190/geo2011-0410.1},
	url = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2012/Li11TRfrfwi/Li11TRfrfwi.pdf},
	author = {Xiang Li and Aleksandr Y. Aravkin and Tristan van Leeuwen and Felix J. Herrmann}
}
@conference {aravkin2012ICASSPfastseis,
	title = {Fast seismic imaging for marine data},
	booktitle = {ICASSP},
	year = {2012},
	publisher = {ICASSP},
	organization = {ICASSP},
	abstract = {Seismic imaging can be formulated as a linear inverse problem where a medium perturbation is obtained via minimization of a least-squares misfit functional. The demand for higher resolution images in more geophysically complex areas drives the need to develop techniques that handle problems of tremendous size with limited computational resources. While seismic imaging is amenable to dimensionality reduction techniques that collapse the data volume into a smaller set of "super-shots", these techniques break down for complex acquisition geometries such as marine acquisition, where sources and receivers move during acquisition. To meet these challenges, we propose a novel method that combines sparsity-promoting (SP) solvers with random subset selection of sequential shots, yielding a SP algorithm that only ever sees a small portion of the full data, enabling its application to very large-scale problems. Application of this technique yields excellent results for a complicated synthetic, which underscores the robustness of sparsity promotion and its suitability for seismic imaging.},
	keywords = {ICASSP},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ICASSP/2012/AravkinLiHerrmann/AravkinLiHerrmann.pdf},
	author = {Aleksandr Y. Aravkin and Xiang Li and Felix J. Herrmann}
}
@presentation {herrmann2012SINBADals,
	title = {Fast sparsity-promoting imaging with message passing},
	year = {2012},
	publisher = {SINBAD},
	abstract = {To meet current-day challenges, exploration seismology increasingly relies on more and more sophisticated algorithms that require multiple paths through all data. This requirement leads to problems because the size of seismic data volumes is increasing exponentially, exposing bottlenecks in IO and computational capability. To overcome these bottlenecks, we follow recent trends in machine learning and compressive sensing by proposing a sparsity-promoting inversion technique that works on small randomized subsets of data only. We boost the performance of this algorithm significantly by modifying a state-of-the-art l1-norm solver to benefit from message passing, which breaks the build up of correlations between model iterates and the randomized linear forward model. We demonstrate the performance of this algorithm on a toy sparse-recovery problem and on a realistic reverse-time-migration example with random source encoding. The improvements in speed, memory use, and output quality are truly remarkable.},
	keywords = {Presentation, SINBAD, SINBADFALL2012, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2012/Fall/herrmann2012SINBADals/herrmann2012SINBADals_pres.pdf},
	author = {Felix J. Herrmann}
}
@article {Herrmann11TRfcd,
	title = {Fighting the curse of dimensionality: compressive sensing in exploration seismology},
	journal = {Signal Processing Magazine, IEEE},
	volume = {29},
	number = {3},
	year = {2012},
	month = {05},
	pages = {88-100},
	abstract = {Many seismic exploration techniques rely on the collection of massive data volumes that are mined for information during processing. This approach has been extremely successful, but current efforts toward higher resolution images in increasingly complicated regions of Earth continue to reveal fundamental shortcomings in our typical workflows. The "curse" of dimensionality is the main roadblock and is 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.},
	keywords = {dimensionality curse, Earth, geographic information systems, higher-resolution images, massive data volumes, Nyquist sampling criterion, seismic exploration techniques, seismology, strains current acquisition system, strains current processing system},
	issn = {1053-5888},
	doi = {10.1109/MSP.2012.2185859},
	url = {https://slim.gatech.edu/Publications/Public/Journals/IEEESignalProcessingMagazine/2012/Herrmann11TRfcd/Herrmann11TRfcd.pdf},
	author = {Felix J. Herrmann and Michael P. Friedlander and Ozgur Yilmaz}
}
@techreport {vanleeuwen2012CGMN,
	title = {Fourier analysis of the CGMN method for solving the Helmholtz equation},
	number = {TR-EOAS-2012-1},
	year = {2012},
	institution = {Department of Earth, Ocean and Atmospheric Sciences},
	address = {The University of British Columbia, Vancouver},
	abstract = {The Helmholtz equation arises in many applications, such as seismic and medical imaging. These application are characterized by the need to propagate many wavelengths through an inhomogeneous medium. The typical size of the problems in 3D applications precludes the use of direct factorization to solve the equation and hence iterative methods are used in practice. For higher wavenumbers, the system becomes increasingly indefinite and thus good preconditioners need to be constructed. In this note we consider an accelerated Kazcmarz method (CGMN) and present an expression for the resulting iteration matrix. This iteration matrix can be used to analyze the convergence of the CGMN method. In particular, we present a Fourier analysis for the method applied to the 1D Helmholtz equation. This analysis suggests an optimal choice of the relaxation parameter. Finally, we present some numerical experiments.},
	keywords = {Helmholtz equation, Modelling},
	url = {http://arxiv.org/abs/1210.2644},
	author = {Tristan van Leeuwen}
}
@presentation {peters2012SINBADfde,
	title = {Frequency domain 3D elastic wave propagation in general anisotropic media},
	year = {2012},
	publisher = {SINBAD},
	abstract = {Elastic wave propagation in 3 spatial dimensions is modeled using a wave equation containing the full stiffness tensor consisting of 21 independent components. This allows modeling in general anisotropic media. The wave equation is discretized on several Cartesian and rotated Cartesian staggered finite-difference grids (using a 2nd order approximation). The grids are linearly combined and, in combination with a antilumped mass strategy, minimize numerical dispersion while requiring a low number of grid points per wavelength. In case not all 21 components need to be modeled, an approximation of the stiffness tensor can be used (e.g., orthorhombic anisotropy, TTI, ...). This results in a linear system of equations, which is solved using an iterative method. The modeling of all 21 components of the stiffness tensor (or an approximation) enables the development of new waveform inversion functionalities.},
	keywords = {Presentation, SINBAD, SINBADFALL2012, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2012/Fall/peters2012SINBADfde/peters2012SINBADfde_pres.pdf},
	author = {Bas Peters}
}
@conference {min2012EAGEefwi,
	title = {Frequency-domain elastic waveform inversion using weighting factors related to source-deconvolved residuals},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2012},
	month = {06},
	abstract = {One of the limitations in seismic waveform inversion is that inversion results are very sensitive to initial guesses, which may be because the gradients computed at each frequency are not properly weighted depending on given models.Analyzingthe conventional waveform inversion algorithms using the pseudo-Hessian matrix as a pre-conditioner shows that the gradientsdo not properly describe the feature of given models or high- and low-end frequencies do not contribute the model parameter updates due to banded spectra of source wavelet. For a better waveform inversion algorithm, we propose applying weighting factors to gradients computed at each frequency. The weighting factors are designed using the source-deconvolved back-propagated wavefields. Numerical results for the SEG/EAGE salt model show that the weighting method improves gradient images and its inversion results are compatible with true velocities even with poorly estimated initial guesses.},
	keywords = {EAGE, elastic, frequency-domain, waveform inversion, weighting factors},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2012/min2012EAGEefwi/min2012EAGEefwi.pdf},
	url2 = {http://earthdoc.eage.org/publication/publicationdetails/?publication=59623},
	author = {Ju-Won Oh and Dong-Joo Min and Felix J. Herrmann}
}
@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 {herrmann2011SEGffw,
	title = {Fast full-waveform inversion with compressive sensing},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2011},
	keywords = {Full-waveform inversion, SEG, SLIM},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2011/HerrmannSEG2011fws/HerrmannSEG2011fws.pdf},
	author = {Aleksandr Y. Aravkin and Felix J. Herrmann and Tristan van Leeuwen and Xiang Li}
}
@presentation {vanleeuwen2011SINBADfwi,
	title = {Fast waveform inversion without source encoding},
	year = {2011},
	publisher = {SINBAD},
	keywords = {Presentation, SINBAD, SINBADFALL2011, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2011/Fall/vanleeuwen2011SINBADfwi/vanleeuwen2011SINBADfwi_pres.pdf},
	author = {Tristan van Leeuwen}
}
@conference {herrmann2011SLRAfwi,
	title = {Full waveform inversion with compressive updates},
	booktitle = {SLRA},
	year = {2011},
	publisher = {Sparse and Low Rank Approximation 2011},
	organization = {Sparse and Low Rank Approximation 2011},
	abstract = {Full-waveform inversion relies on large multi-experiment data volumes. While improvements in acquisition and inversion have been extremely successful, the current push for higher quality models reveals fundamental shortcomings handling increasing problem sizes numerically. To address this fundamental issue, we propose a randomized dimensionality-reduction strategy motivated by recent developments in stochastic optimization and compressive sensing. In this formulation conventional Gauss-Newton iterations are replaced by dimensionality-reduced sparse recovery problems with source encodings.},
	keywords = {Full-waveform inversion},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SLRA/2011/herrmann2011SLRAfwi/herrmann2011SLRAfwi.pdf},
	author = {Felix J. Herrmann and Aleksandr Y. Aravkin and Xiang Li and Tristan van Leeuwen}
}
@conference {aravkin2011SIAMfwi,
	title = {Full waveform inversion with compressive updates},
	booktitle = {SIAM},
	year = {2011},
	publisher = {SIAM CS\&E 2011},
	organization = {SIAM CS\&E 2011},
	abstract = {Full-waveform inversion relies on large multi-experiment data volumes. While improvements in acquisition and inversion have been extremely successful, the current push for higher quality models reveals fundamental shortcomings handling increasing problem sizes numerically. To address this fundamental issue, we propose a randomized dimensionality-reduction strategy motivated by recent developments in stochastic optimization and compressive sensing. In this formulation conventional Gauss-Newton iterations are replaced by dimensionality-reduced sparse recovery problems with source encodings.},
	keywords = {Full-waveform inversion, SLIM},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SIAM/2011/aravkin2011SIAMfwi/aravkin2011SIAMfwi.pdf},
	author = {Aleksandr Y. Aravkin and Felix J. Herrmann and Tristan van Leeuwen and James V. Burke and Xiang Li}
}
@conference {li2011EAGEfwirr,
	title = {Full-waveform inversion with randomized L1 recovery for the model updates},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2011},
	month = {05},
	abstract = {Full-waveform inversion (FWI) is a data fitting procedure that relies on the collection of seismic data volumes and sophisticated computing to create high-resolution results. With the advent of FWI, the improvements in acquisition and inversion have been substantial, but these improvements come at a high cost because FWI involves extremely large multi-experiment data volumes. The main obstacle is the "curse of dimensionality" exemplified by Nyquist{\~O}s sampling criterion, which puts a disproportionate strain on current acquisition and processing systems as the size and desired resolution increases. In this paper, we address the "curse of dimensionality" by randomized dimensionality reduction of the FWI problem adapted from the field of CS. We invert for model updates by replacing the Gauss-Newton linearized subproblem for subsampled FWI with a sparsity promoting formulation, and solve this formulation using the SPGl1 algorithm. We speed up the algorithm and avoid overfitting the data by solving for the linearized updates only approximately. Our approach is successful because it reduces the size of seismic data volumes without loss of information. With this reduction, we can compute a Newton-like update with the reduced data volume at the cost of roughly one gradient update for the fully sampled wavefield.},
	keywords = {EAGE, Full-waveform inversion, Presentation},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2011/li11EAGEfwirr/li11EAGEfwirr.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2011/li11EAGEfwirr/li11EAGEfwirr_pres.pdf},
	url2 = {http://earthdoc.eage.org/publication/publicationdetails/?publication=50345},
	author = {Xiang Li and Aleksandr Y. Aravkin and Tristan van Leeuwen and Felix J. Herrmann}
}
@conference {herrmann2011ICIAMconvexcompfwi,
	title = {FWI with sparse recovery: a convex-composite approach},
	booktitle = {ICIAM},
	year = {2011},
	month = {07},
	publisher = {ICIAM 2011},
	organization = {ICIAM 2011},
	abstract = {Iterative inversion algorithms require repeated simulation of 3D time-dependent acoustic, elastic, or electromagnetic wave fields, extending hundreds of wavelengths and hundreds of periods. Also, seismic data is rich in information at every representable scale. Thus simulation-driven optimization approaches to inversion impose great demands on simulator efficiency and accuracy. While computer hardware advances have been of critical importance in bringing inversion closer to practical application, algorithmic advances in simulator methodology have been equally important. Speakers in this two-part session will address a variety of numerical issues arising in the wave simulation, and in its application to inversion.},
	keywords = {Full-waveform inversion, ICIAM, Optimization},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ICIAM/2011/herrmann2011ICIAMconvexcompfwi/herrmann2011ICIAMconvexcompfwi.pdf},
	author = {Felix J. Herrmann and Aleksandr Y. Aravkin and Tristan van Leeuwen and Xiang Li}
}
@conference {li2010SEGfwi,
	title = {Full-waveform inversion from compressively recovered model updates},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {29},
	year = {2010},
	month = {10},
	pages = {1029-1033},
	publisher = {SEG},
	organization = {SEG},
	abstract = {Full-waveform inversion relies on the collection of large multi-experiment data volumes in combination with a sophisticated back-end to create high-fidelity inversion results. While improvements in acquisition and inversion have been extremely successful, the current trend of incessantly pushing for higher quality models in increasingly complicated regions of the Earth reveals fundamental shortcomings in our ability to handle increasing problem size numerically. Two main culprits can be identified. First, there is the so-called {\textquoteleft}{\textquoteleft}curse of dimensionality{\textquoteright}{\textquoteright} exemplified by Nyquist{\textquoteright}s sampling criterion, which puts disproportionate strain on current acquisition and processing systems as the size and desired resolution increases. Secondly, there is the recent {\textquoteleft}{\textquoteleft}departure from Moore{\textquoteright}s law{\textquoteright}{\textquoteright} that forces us to lower our expectations to compute ourselves out of this. In this paper, we address this situation by randomized dimensionality reduction, which we adapt from the field of compressive sensing. In this approach, we combine deliberate randomized subsampling with structure-exploiting transform-domain sparsity promotion. Our approach is successful because it reduces the size of seismic data volumes without loss of information. With this reduction, we compute Newton-like updates at the cost of roughly one gradient update for the fully-sampled wavefield.},
	keywords = {Full-waveform inversion, Presentation, SEG},
	doi = {10.1190/1.3513022},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2010/li10SEGfwi/li10SEGfwi.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2010/li10SEGfwi/li10SEGfwi_pres.pdf},
	author = {Xiang Li and Felix J. Herrmann}
}
@presentation {li2010SINBADfwi,
	title = {Full-waveform inversion with randomized L1 recovery for the model updates},
	year = {2010},
	publisher = {SINBAD},
	keywords = {Presentation, SINBAD, SINBADFALL2010, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2010/Fall/li2010SINBADfwi/li2010SINBADfwi_pres.pdf},
	author = {Xiang Li}
}
@conference {herrmann2007SINBADfrw,
	title = {Focused recovery with the curvelet transform},
	booktitle = {SINBAD 2007},
	year = {2007},
	abstract = {Incomplete data represents a major challenge for a successful prediction and subsequent removal of multiples. In this paper, a new method will be represented that tackles this challenge in a two-step approach. During the first step, the recently developed curvelet-based recovery by sparsity-promoting inversion (CRSI) is applied to the data, followed by a prediction of the primaries. During the second high-resolution step, the estimated primaries are used to improve the frequency content of the recovered data by combining the focal transform, defined in terms of the estimated primaries, with the curvelet transform. This focused curvelet transform leads to an improved recovery, which can subsequently be used as input for a second stage of multiple prediction and primary-multiple separation. This is joint work with Deli Wang and Gilles Hennenfent.},
	keywords = {Presentation, SINBAD, SLIM},
	presentation = {https://slim.gatech.edu/sites/data/Papers/Herrmann2007SINBADfoc.pdf},
	author = {Felix J. Herrmann}
}
@conference {herrmann2007SLIMfsd,
	title = {From seismic data to the composition of rocks: an interdisciplinary and multiscale approach to exploration seismology},
	booktitle = {Berkhout{\textquoteright}s valedictory address: the conceptual approach of understanding},
	year = {2007},
	abstract = {In this essay, a nonlinear and multidisciplinary approach is presented that takes seismic data to the composition of rocks. The presented work has deep roots in the {\textquoteleft}gedachtengoed{\textquoteright} (philosophy) of Delphi spearheaded by Guus Berkhout. Central themes are multiscale, object-orientation and a multidisciplinary approach.},
	keywords = {SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/Misc/herrmann07SLIMfsd/herrmann07SLIMfsd.pdf},
	author = {Felix J. Herrmann}
}
