@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 {alfaraj2017SEGrsw,
	title = {Reconstruction of S-waves from low-cost randomized and simultaneous acquisition by joint sparse inversion},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2017},
	note = {(SEG, Houston)},
	month = {09},
	pages = {2533-2538},
	abstract = {Spatial Nyquist sampling rate, which is proportional to the apparent subsurface velocity, is the key deciding cost factor for conventional seismic acquisition. Due to the lower shear wave velocity compared with compressional waves, finer spatial sampling is required to properly record the earlier, which increases the acquisition costs. This is one of the reasons why shear waves are usually not considered in practice. To avoid having the Nyquist criterion as the deciding cost factor and to utilize multicomponent data to its available full extent, we propose acquiring randomly undersampled ocean bottom seismic data. Each component can then be interpolated separately, followed by elastic decomposition to recover up- and down-going S-waves. Instead, we jointly interpolate and decompose the recorded multicomponent data by solving one sparsity promoting optimization problem. This way we ensure that the relative amplitudes of the multicomponent data is preserved. We compare two sparsifying transforms: the curvelet transform and the frequency-wavenumber transform. Another key cost deciding factor for seismic acquisition is the efficiency of acquiring data. This calls for simultaneous acquisition, which requires a source separation step. Similarly, instead of taking a two-step approach, we perform a sparsity-promoting joint source separation decomposition. Results on economically and efficiently acquired synthetic data of both joint methods show their ability of reconstructing accurate up- and down-going S-waves.},
	keywords = {decomposition, elastic, interpolation, randomized acquisition, SEG, shear wave, source separation, sparse},
	doi = {10.1190/segam2017-17681376.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2017/alfaraj2017SEGrsw/alfaraj2017SEGrsw.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2017/alfaraj2017SEGrsw/alfaraj2017SEGrsw_pres.pdf},
	author = {Ali M. Alfaraj and Rajiv Kumar and Felix J. Herrmann}
}
@conference {li2013SEGodmvdaiedwawe,
	title = {Optimization driven model-space versus data-space approaches to invert elastic data with the acoustic wave equation},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {32},
	year = {2013},
	month = {09},
	pages = {986-990},
	abstract = {Inverting data with elastic phases using an acoustic wave equation can lead to erroneous results, especially when the number of iterations is too high, which may lead to over fitting the data. Several approaches have been proposed to address this issue. Most commonly, people apply "data-independent" filtering operations that are aimed to deemphasize the elastic phases in the data in favor of the acoustic phases. Examples of this approach are nested loops over offset range and Laplace parameters. In this paper, we discuss two complementary optimization-driven methods where the minimization process decides adaptively which of the data or model components are consistent with the objective. Specifically, we compare the Student{\textquoteright}s t misfit function as the data-space alternative and curvelet-domain sparsity promotion as the model-space alternative. Application of these two methods to a realistic synthetic lead to comparable results that we believe can be improved by combining these two methods.},
	keywords = {elastic, Full-waveform inversion, least-squares, SEG},
	doi = {10.1190/segam2013-1375.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2013/li2013SEGodmvdaiedwawe/li2013SEGodmvdaiedwawe.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2013/li2013SEGodmvdaiedwawe/li2013SEGodmvdaiedwawe_pres.pdf},
	author = {Xiang Li and Anais Tamalet and Tristan van Leeuwen and Felix J. Herrmann}
}
@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 {min2012CSEGrgfe,
	title = {Re-establishment of gradient in frequency-domain elastic waveform inversion},
	booktitle = {CSEG Annual Conference Proceedings},
	year = {2012},
	month = {02},
	abstract = {To obtain solutions close to global minimum in waveform inversion, the gradients computed at each frequency need to be weighted to appropriately describe the residuals between modeled and field data. While the low-frequency components of the gradients should be weighted to recover the long-wavelength structures, the high-frequency components of the gradients need to be weighted when the short-wavelength structures are restored. However, the conventional elastic waveform inversion algorithms cannot properly weight the gradients computed at each frequency. When gradients are scaled using the pseudo-Hessian matrix inside the frequency loop, gradients obtained at high frequencies are over-emphasized. When the gradients are scaled outside the frequency loop, gradients are weighted by the source spectra. In this study, we propose applying weighting factors to the gradients obtained at each frequency so that gradients can properly reflect the differences between the true and assumed models satisfying the general inverse theory. The weighting factors are composed by the backpropagated residuals. Numerical examples for the simple rectangular-shaped model and the modified version of the Marmousi-2 model show that the weighting method enhances gradient images and inversion results compared to the conventional inversion algorithms.},
	keywords = {elastic, frequency-domain, waveform inversion, weighting factors},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/CSEG/2012/min2012CSEGrgfe/min2012CSEGrgfe.pdf},
	author = {Ju-Won Oh and Dong-Joo Min and Felix J. Herrmann}
}
