@conference {daskalakis2017SEGdds,
	title = {Developments in the direction of solving extremly large problems in Geophysics},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2017},
	note = {(SEG, Houston)},
	month = {09},
	pages = {4375-4378},
	abstract = {Often in exploration Geophysics we are forced to work with extremely large problems. Acquisitions via dense grids of receivers translate into very large mathematical systems. Usually, depending on the acquisition, the size of the matrix of the system to be solved, can be measured in the millions. The proper way to address this problem is by subsampling. Even though subsampling can reduce the computational efforts required, it can not address stability problems caused by preconditioning and/or instrumental response errors. In this abstract, we introduce a modification of the linearized Bregman solver for these large problems that resolves stability issues.},
	keywords = {cycling, least-squares migration, linearized Bregman, SEG, weighted increment},
	doi = {10.1190/segam2017-17795188.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2017/daskalakis2017SEGdds/daskalakis2017SEGdds.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2017/daskalakis2017SEGdds/daskalakis2017SEGdds_pres.pdf},
	author = {Emmanouil Daskalakis and Rachel Kuske 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}
}
@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}
}
@conference {tu2015PIMSsls,
	title = {Sparse least-squares seismic imaging with source estimation utilizing multiples},
	booktitle = {PIMS Workshop on Advances in Seismic Imaging and Inversion},
	year = {2015},
	note = {(PIMS Workshop on Advances in Seismic Imaging and Inversion, University of Alberta, Edmonton)},
	month = {05},
	abstract = {We present an least-squares seismic imaging method that (1) is efficient by subsampling the monochromatic source experiments; (2) does not require the prior knowledge of the source wavelet, by inverting for the seismic image and the source wavelet simultaneously using variable projection; and (3) makes active use of surface-related multiples to increase the illumination and to resolve the scaling ambiguity in source estimation, by tightly integrating the well-established surface-related-multiple-prediction into the wave-equation based modelling. As a result, we are able to obtain ambiguity-resolved, high-fidelity least-squares inverted seismic images with computational costs that are comparable to conventional reverse-time migrations with all the data. We demonstrate the efficacy of the proposed method using synthetic examples.},
	keywords = {inversion, least-squares migration, multiples, PIMS, seismic imaging, source estimation, workshop},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/PIMS/2015/tu2015PIMSsls/tu2015PIMSsls_pres.pdf},
	author = {Ning Tu and Felix J. Herrmann}
}
