@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 {alfaraj2018EAGEasa,
	title = {Automatic statics and residual statics correction with low-rank approximation},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2018},
	note = {(EAGE, Copenhagen)},
	month = {06},
	abstract = {Seismic data, mainly on land, suffers from long-wavelength statics due to the laterally varying and heterogeneous nature of the near-surface weathering layers. We propose an automatic, data-driven and computationally efficient statics correction method based on low-rank approximation to correct for such statics. The method does not require a model to estimate static time shifts, which is the case for other static correction methods; rather it applies the appropriate static corrections on the data such that it becomes low rank in a certain domain. As of now, the method is applicable to data that has been corrected for elevation statics. Due to the near-surface irregularities and due to approximations used by static correction methods that lead to not fully correcting for statics, an iterative residual statics correction becomes necessary. Our proposed method corrects for residual statics without the necessity of the surface consistency assumption and a multi-iterate process. Additional benefits of the method include artifacts and noise suppression. We demonstrate the successful application of our method on several synthetic data examples.},
	keywords = {EAGE, land, Rank, residual statics, statics},
	doi = {10.3997/2214-4609.201801107},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2018/alfaraj2018EAGEasa/alfaraj2018EAGEasa.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2018/alfaraj2018EAGEasa/alfaraj2018EAGEasa_pres.pdf},
	author = {Ali M. Alfaraj and Rajiv Kumar and Felix J. Herrmann}
}
@conference {alfaraj2018SEGswi,
	title = {Seismic waveform inversion using decomposed one-way wavefields},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2018},
	note = {(SEG, Anaheim)},
	month = {10},
	pages = {1379-1383},
	abstract = {Conventional seismic full waveform inversion (FWI) estimates a velocity model by matching the two-way predicted data with the observed data. The non-linearity of the problem and the complexity of the Earth{\textquoteright}s subsurface, which results in complex wave propagation may lead to unsatisfactory inversion results. Moreover, having inaccurate overburden models, i.e water velocity model in the ocean bottom acquisition setting, may lead to erroneous velocity models. In order to obtain better velocity models, we simplify the FWI problem by decomposing the total two-way observed and predicted wavefields into one-way wavefields at the receiver locations using acoustic wavefield decomposition. We then propose to fit the less-complex one-way wavefield to obtain a velocity model that contains valuable information about the overburden. We use this inverted velocity model as a better initial model for the inversion using the other more-complex one-way wavefield. We demonstrate our proposed algorithm on acoustic non-inverse crime synthetic data produced from the Marmousi2 model. The proposed method provides improved inversion results compared with conventional FWI due to properly accounting for separate model updates from the up- and down-going wavefields.},
	keywords = {decomposition, FWI, one-way, SEG},
	doi = {10.1190/segam2018-2998162.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2018/alfaraj2018SEGswi/alfaraj2018SEGswi.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2018/alfaraj2018SEGswi/alfaraj2018SEGswi_pres.pdf},
	author = {Ali M. Alfaraj and Rajiv Kumar and Felix J. Herrmann}
}
@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 {alfaraj2017EAGEswr,
	title = {Shear wave reconstruction from low cost randomized acquisition},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2017},
	note = {(EAGE, Paris)},
	month = {06},
	abstract = {Shear waves travel in the subsurface at a lower speed compared with compressional waves. Therefore, much finer spatial sampling is required to properly record the shear waves. This leads to higher acquisition costs which are typically avoided by designing surveys geared towards only compressional waves imaging. We propose using randomly under-sampled ocean bottom acquisition for recording both compressional and shear waves. The recorded multicomponent data is then interpolated using an SVD-free low rank interpolation scheme that is feasible for large scale seismic data volumes to obtain finely sampled data. Following that, we perform elastic wavefield decomposition at the ocean bottom to recover accurate up- and dow-going S-waves. Synthetic data results indicate that using randomized under-sampled acquisition, we can recover accurate S-waves with an economical cost compared with conventional acquisition designs.},
	keywords = {EAGE, interpolation, randomized acquisition, rank minimization, shear waves},
	doi = {10.3997/2214-4609.201700594},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2017/alfaraj2017EAGEswr/alfaraj2017EAGEswr.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2017/alfaraj2017EAGEswr/alfaraj2017EAGEswr_pres.pdf},
	author = {Ali M. Alfaraj and Rajiv Kumar and Felix J. Herrmann}
}
@conference {esser2015SEGasd,
	title = {Automatic salt delineation {\textendash}- Wavefield Reconstruction Inversion with convex constraints},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2015},
	note = {(SEG, New Orleans)},
	month = {10},
	pages = {1337-1343},
	abstract = {We extend full-waveform inversion by Wavefield Reconstruction Inversion by including convex constraints on the model. Contrary to the conventional adjoint-state formulations, Wavefield Reconstruction Inversion has the advantage that the Gauss-Newton Hessian is well approximated by a diagonal scaling, which allows us to add convex constraints, such as the box- and the edge-preserving total-variation constraint, on the square slowness without incurring significant increases in computational costs. As the examples demonstrate, including these constraints yields far superior results in complex geological areas that contain high-velocity high-contrast bodies (e.g. salt or basalt). Without these convex constraints, adjoint-state and Wavefield Reconstruction Inversion get trapped in local minima for poor starting models.},
	keywords = {convex constraints, Full-waveform inversion, salt, SEG, total-variation norm},
	doi = {10.1190/segam2015-5877995.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2015/esser2015SEGasd/esser2015SEGasd.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2015/esser2015SEGasd/esser2015SEGasd_pres.pdf},
	author = {Ernie Esser and Llu{\'\i}s Guasch and Tristan van Leeuwen and Aleksandr Y. Aravkin and Felix J. Herrmann}
}
@conference {kumar2015IIPFWIrmb,
	title = {Rank minimization based seismic data processing and inversion},
	booktitle = {Inaugural Full-Waveform Inversion Workshop},
	year = {2015},
	note = {(Natal, Brazil)},
	month = {08-09},
	abstract = {We present the low-rank extensions of full-waveform inversion to mitigate the local-minima issues. The purpose formulation is computationally tractable and does not involve any extra computational cost compared to the adjoint-state method.},
	keywords = {FWI, inversion, Processing, rank minimization},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/IIPFWI/kumar2015IIPFWIrmb/kumar2015IIPFWIrmb_pres.pdf},
	author = {Rajiv Kumar and Curt Da Silva and Oscar Lopez and Aleksandr Y. Aravkin and Hassan Mansour and Haneet Wason and Ernie Esser and Felix J. Herrmann}
}
@conference {wason2014SEGsss,
	title = {Source separation via SVD-free rank minimization in the hierarchical semi-separable representation},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2014},
	note = {(SEG)},
	month = {10},
	pages = {120-126},
	abstract = {Recent developments in matrix rank optimization have allowed for new computational approaches in the field of source separation or deblending. In this paper, we propose a source separation algorithm for blended marine acquisition, where two sources are deployed at different depths (over/under acquisition). The separation method incorporates the Hierarchical Semi-Separable structure (HSS) inside rank-regularized least-squares formulations. The proposed approach is suitable for large scale problems, since it avoids SVD computations and uses a low-rank factorized formulation instead. We illustrate the performance of the new HSS-based deblending approach by simulating an over/under blended acquisition, wherein uniformly random time delays (of < 1 second) are applied to one of the sources.},
	keywords = {Acquisition, deblending, HSS, marine, Rank, SEG, source separation},
	doi = {http://dx.doi.org/10.1190/segam2014-1583.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2014/wason2014SEGsss/wason2014SEGsss.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2014/wason2014SEGsss/wason2014SEGsss_pres.pdf},
	author = {Haneet Wason and Rajiv Kumar and Aleksandr Y. Aravkin and Felix J. Herrmann}
}
@conference {kumar2014EAGErank,
	title = {SVD-free low-rank matrix factorization : wavefield reconstruction via jittered subsampling and reciprocity},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2014},
	month = {06},
	abstract = {Recently computationally efficient rank optimization techniques have been studied extensively to develop a new mathematical tool for the seismic data interpolation. So far, matrix completion problems have been discussed where sources are subsample according to a discrete uniform distribution. In this paper, we studied the effect of two different subsampling techniques on seismic data interpolation using rank-regularized formulations, namely jittered subsampling over uniform random subsampling. The other objective of this paper is to combine the fact of source-receiver reciprocity with the rank-minimization techniques to enhance the accuracy of missing-trace interpolation. We illustrate the advantages of jittered subsampling and reciprocity using a seismic line from Gulf of Suez to obtain high quality results for interpolation, a key application in exploration geophysics.},
	keywords = {EAGE, interpolation, low-rank, reciprocity},
	doi = {10.3997/2214-4609.20141394},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2014/kumar2014EAGErank/kumar2014EAGErank.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2014/kumar2014EAGErank/kumar2014EAGErank_pres.pdf},
	author = {Rajiv Kumar and Aleksandr Y. Aravkin and Ernie Esser and Hassan Mansour 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 {vanleeuwen2013EAGErobustFWI,
	title = {In which domain should we measure the misfit for robust full waveform inversion?},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2013},
	month = {06},
	abstract = {Full-waveform inversion relies on minimizing the difference between observed and modeled data, as measured by some penalty function. A popular choice, of course, is the least-squares penalty. However, when outliers are present in the data, the use of robust penalties such as the Huber or Student{\textquoteright}s t may significantly improve the results since they put relatively less weight on large residuals. In order for robust penalties to be effective, the outliers must be somehow localized and distinguishable from the good data. We propose to first transform the residual into a domain where the outliers are localized before measuring the misfit with a robust penalty. This is exactly how one would normally devise filters to remove the noise before applying conventional FWI. We propose to merge the two steps and let the inversion process implicitly filter out the noise. Results on a synthetic dataset show the effectiveness of the approach.},
	keywords = {EAGE, full waveform inversion},
	doi = {10.3997/2214-4609.20130839},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2013/vanleeuwen2013EAGErobustFWI/vanleeuwen2013EAGErobustFWI.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2013/vanleeuwen2013EAGErobustFWI/vanleeuwen2013EAGErobustFWI_pres.pdf},
	author = {Tristan van Leeuwen and Aleksandr Y. Aravkin and Henri Calandra and Felix J. Herrmann}
}
@conference {kumar2013SEGHSS,
	title = {Reconstruction of seismic wavefields via low-rank matrix factorization in the hierarchical-separable matrix representation},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {32},
	year = {2013},
	month = {09},
	pages = {3628-3633},
	abstract = {Recent developments in matrix rank optimization have allowed for new computational approaches in the field of seismic data interpolation. In this paper, we propose an approach for seismic data interpolation which incorporates the Hierarchical Semi-Separable Structure (HSS) inside rank-regularized least-squares formulations for the missing-trace interpolation problem. The proposed approach is suitable for large scale problems, since it avoids SVD computations and uses a low-rank factorized formulation instead. We illustrate the advantages of the new HSS approach by interpolating a seismic line from the Gulf of Suez and compare the reconstruction with conventional rank minimization.},
	keywords = {HSS, interpolation, SEG},
	doi = {10.1190/segam2013-1165.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2013/kumar2013SEGHSS/kumar2013SEGHSS.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2013/kumar2013SEGHSS/kumar2013SEGHSS_pres.pdf},
	author = {Rajiv Kumar and Hassan Mansour and Aleksandr Y. Aravkin and Felix J. Herrmann}
}
@conference {kumar2013EAGEsind,
	title = {Seismic data interpolation and denoising using SVD-free low-rank matrix factorization},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2013},
	month = {06},
	abstract = {Recent developments in rank optimization have allowed new approaches for seismic data interpolation and denoising. In this paper, we propose an approach for simultaneous seismic data interpolation and denoising using robust rank-regularized formulations. The proposed approach is suitable for large scale problems, since it avoids SVD computations by using factorized formulations. We illustrate the advantages of the new approach using a seismic line from Gulf of Suez and 5D synthetic seismic data to obtain high quality results for interpolation and denoising, a key application in exploration geophysics.},
	keywords = {denoising, EAGE, interpolation},
	doi = {10.3997/2214-4609.20130388},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2013/kumar2013EAGEsind/kumar2013EAGEsind.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2013/kumar2013EAGEsind/kumar2013EAGEsind_pres.pdf},
	author = {Rajiv Kumar and Aleksandr Y. Aravkin and Hassan Mansour and Ben Recht and Felix J. Herrmann}
}
@conference {aravkin2013ICASSPssi,
	title = {Sparse seismic imaging using variable projection},
	booktitle = {ICASSP},
	year = {2013},
	month = {05},
	abstract = {We consider an important class of signal processing problems where the signal of interest is known to be sparse, and can be recovered from data given auxiliary information about how this data was generated. For example, a sparse green{\textquoteright}s function may be recovered from seismic experimental data using sparsity optimization when the source signature is known. Unfortunately, in practice this information is often missing, and must be recovered from data along with the signal using deconvolution techniques. In this paper, we present a novel methodology to simulta- neously solve for the sparse signal and auxiliary parameters using a recently proposed variable projection technique. Our main contribution is to combine variable projection with spar- sity promoting optimization, obtaining an efficient algorithm for large-scale sparse deconvolution problems. We demon- strate the algorithm on a seismic imaging example.},
	keywords = {Imaging, Optimization, sparsity, variable projection},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ICASSP/2013/aravkin2013ICASSPssi/aravkin2013ICASSPssi.pdf},
	author = {Aleksandr Y. Aravkin and Tristan van Leeuwen and Ning Tu}
}
@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}
}
@conference {aravkin2012SEGST,
	title = {On non-uniqueness of the Student{\textquoteright}s t-formulation for linear inverse problems},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {31},
	year = {2012},
	month = {11},
	pages = {1-5},
	publisher = {SEG},
	organization = {SEG},
	abstract = {We review the statistical interpretation of inverse problem formulations, and the motivations for selecting non-convex penalties for robust behaviour with respect to measurement outliers or artifacts in the data. An important downside of using non-convex formulations such as the Student{\textquoteright}s t is the potential for non-uniqueness, and we present a simple example where the Student{\textquoteright}s t penalty can be made to have many local minima by appropriately selecting the degrees of freedom parameter. On the other hand, the non-convexity of the Student{\textquoteright}s t is precisely what gives it the ability to ignore artifacts in the data. We explain this idea, and present a stylized imaging experiment, where the Student{\textquoteright}s t is able to recover a velocity perturbation from data contaminated by a very peculiar artifact {\textendash}- data from a different velocity perturbation. The performance of Student{\textquoteright}s t inversion is investigated empirically for different values of the degrees of freedom parameter, and different initial conditions.},
	keywords = {non-convex, robust, SEG, Student{\textquoteright}s t, uniqueness},
	doi = {10.1190/segam2012-1558.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2012/aravkin2012SEGST/aravkin2012SEGST.pdf},
	author = {Aleksandr Y. Aravkin and Tristan van Leeuwen and Kenneth Bube and Felix J. Herrmann}
}
@conference {aravkin2012ICASSProbustb,
	title = {Robust inversion via semistochastic dimensionality reduction},
	booktitle = {ICASSP},
	year = {2012},
	pages = {5245-5248},
	publisher = {ICASSP},
	organization = {ICASSP},
	abstract = {We consider a class of inverse problems where it is possible to aggregate the results of multiple experiments. This class includes problems where the forward model is the solution operator to linear ODEs or PDEs. The tremendous size of such problems motivates the use dimensionality reduction (DR) techniques based on randomly mixing experiments. These techniques break down, however, when robust data-fitting formulations are used, which are essential in cases of missing data, unusually large errors, and systematic features in the data unexplained by the forward model. We survey robust methods within a statistical framework, and propose a sampling optimization approach that allows DR. The efficacy of the methods are demonstrated for a large-scale seismic inverse problem using the robust Student{\textquoteright}s t-distribution, where a useful synthetic velocity model is recovered in the extreme scenario of 60\% corrupted data. The sampling approach achieves this recovery using 20\% of the effort required by a direct robust approach.},
	keywords = {ICASSP},
	doi = {10.1109/ICASSP.2012.6289103},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ICASSP/2012/AravkinFriedlanderLeeuwen/AravkinFriedlanderLeeuwen.pdf},
	author = {Aleksandr Y. Aravkin and Michael P. Friedlander and Tristan van Leeuwen}
}
@conference {aravkin2012EAGErobust,
	title = {Source estimation for frequency-domain FWI with robust penalties},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2012},
	month = {06},
	abstract = {Source estimation is an essential component of full waveform inversion. In the standard frequency domain formulation, there is closed form solution for the the optimal source weights, which can thus be cheaply estimated on the fly. A growing body of work underscores the importance of robust modeling for data with large outliers or artifacts that are not captured by the forward model. Effectively, the least-squares penalty on the residual is replaced by a robust penalty, such as Huber, Hybrid {\textquoteleft}1-{\textquoteleft}2 or Student{\textquoteright}s t. As we will demonstrate, it is essential to use the same robust penalty for source estimation. In this abstract, we present a general approach to robust waveform inversion with robust source estimation. In this general formulation, there is no closed form solution for the optimal source weights so we need to solve a scalar optimization problem to obtain these weights. We can efficiently solve this optimization problem with a Newton-like method in a few iterations. The computational cost involved is of the same order as the usual least-squares source estimation procedure. We show numerical examples illustrating robust source estimation and robust waveform inversion on synthetic data with outliers.},
	keywords = {EAGE},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2012/aravkin2012EAGErobust/aravkin2012EAGErobust.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2012/aravkin2012EAGErobust/aravkin2012EAGErobust_pres.pdf},
	url2 = {http://earthdoc.eage.org/publication/publicationdetails/?publication=59196},
	author = {Aleksandr Y. Aravkin and Tristan van Leeuwen and Henri Calandra 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}
}
@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 {li2011SBGFmgnsu,
	title = {Modified Gauss-Newton with sparse updates},
	booktitle = {SBGF},
	year = {2011},
	publisher = {SBGF},
	organization = {SBGF},
	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 models.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 {\textquoteleft}curse of dimensionality{\textquoteright} exemplified by Nyquist{\textquoteright}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 {\textquoteleft}curse of dimensionality{\textquoteright} by using randomized dimensionality reduction of the FWI problem, coupled with a modified Gauss-Newton (GN) method designed to promote curvelet-domain sparsity of model updates. We solve for these updates using the spectral projected gradient method, implemented in the SPG1 software package. Our approach is successful because it reduces the size of seismic data volumes without loss of information. With this reduction, we can compute Gauss-Newton updates with the reduced data volume at the cost of roughly one gradient update for the fully sampled wavefield},
	keywords = {Full-waveform inversion, SBGF},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SBGF/2011/li11SBGFmgnsu/li11SBGFmgnsu.pdf},
	author = {Xiang Li and Felix J. Herrmann and Tristan van Leeuwen and Aleksandr Y. Aravkin}
}
@conference {herrmann2011SPIEmsp,
	title = {A modified, sparsity promoting, Gauss-Newton algorithm for seismic waveform inversion},
	booktitle = {Proc. SPIE},
	number = {81380V},
	year = {2011},
	month = {08},
	abstract = {Images obtained from seismic data are used by the oil and gas industry for geophysical exploration. Cutting-edge methods for transforming the data into interpretable images are moving away from linear approximations and high-frequency asymptotics towards Full Waveform Inversion (FWI), a nonlinear data-fitting procedure based on full data modeling using the wave-equation. The size of the problem, the nonlinearity of the forward model, and ill-posedness of the formulation all contribute to a pressing need for fast algorithms and novel regularization techniques to speed up and improve inversion results. In this paper, we design a modified Gauss-Newton algorithm to solve the PDE- constrained optimization problem using ideas from stochastic optimization and compressive sensing. More specifically, we replace the Gauss-Newton subproblems by randomly subsampled, -$\ell_1$ regularized subproblems. This allows us us significantly reduce the computational cost of calculating the updates and exploit the compressibility of wavefields in Curvelets. We explain the relationships and connections between the new method and stochastic optimization and compressive sensing (CS), and demonstrate the efficacy of the new method on a large-scale synthetic seismic example.},
	keywords = {Compressive Sensing, Full-waveform inversion, Optimization, SLIM},
	issn = {1},
	doi = {10.1117/12.893861},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SPIE/2011/herrmann2011SPIEmsp/herrmann2011SPIEmsp.pdf},
	author = {Felix J. Herrmann and Xiang Li and Aleksandr Y. Aravkin and Tristan van Leeuwen}
}
@conference {aravkin2011EAGEnspf,
	title = {A nonlinear sparsity promoting formulation and algorithm for full waveform inversion},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2011},
	month = {05},
	abstract = {Full Waveform Inversion (FWI) is a computational procedure to extract medium parameters from seismic data. FWI is typically formulated as a nonlinear least squares optimization problem, and various regularization techniques are used to guide the optimization because the problem is illposed. In this paper, we propose a novel sparse regularization which exploits the ability of curvelets to efficiently represent geophysical images. We then formulate a corresponding sparsity promoting constrained optimization problem, which we call Nonlinear Basis Pursuit Denoise (NBPDN) and present an algorithm to solve this problem to recover medium parameters. The utility of the NBPDN formulation and efficacy of the algorithm are demonstrated on a stylized cross-well experiment, where a sparse velocity perturbation is recovered with higher quality than the standard FWI formulation (solved with LBFGS). The NBPDN formulation and algorithm can recover the sparse perturbation even when the data volume is compressed to 5 \% of the original size using random superposition.},
	keywords = {EAGE, Full-waveform inversion, Optimization},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2011/aravkin11EAGEnspf/aravkin11EAGEnspf.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2011/aravkin11EAGEnspf/aravkin11EAGEnspf_pres.pdf},
	url2 = {http://earthdoc.eage.org/publication/publicationdetails/?publication=50199},
	author = {Aleksandr Y. Aravkin and James V. Burke and Felix J. Herrmann and Tristan van Leeuwen}
}
@conference {aravkin2011SEGrobust,
	title = {Robust full-waveform inversion using the Student{\textquoteright}s t-distribution},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {30},
	year = {2011},
	month = {09},
	pages = {2669-2673},
	publisher = {SEG},
	organization = {SEG},
	abstract = {Full-waveform inversion (FWI) is a computational procedure to extract medium parameters from seismic data. Robust methods for FWI are needed to overcome sensitivity to noise and in cases where modeling is particularly poor or far from the real data generating process. We survey previous robust methods from a statistical perspective, and use this perspective to derive a new robust method by assuming the random errors in our model arise from the Student{\textquoteright}s t-distribution. We show that in contrast to previous robust methods, the new method progres- sively down-weighs large outliers, effectively ignoring them once they are large enough. This suggests that the new method is more robust and suitable for situations with very poor data quality or modeling. Experiments show that the new method recovers as well or better than previous robust methods, and can recover models with quality comparable to standard methods on noise-free data when some of the data is completely corrupted, and even when a marine acquisition mask is entirely ignored in the modeling. The ability to ignore a marine acquisition mask via robust FWI methods offers an opportunity for stochastic optimization methods in marine acquisition.},
	keywords = {Full-waveform inversion, Optimization, SEG},
	doi = {10.1190/1.3627747},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2011/aravkin11SEGrobust/aravkin11SEGrobust.pdf},
	author = {Aleksandr Y. Aravkin and Tristan van Leeuwen and Felix J. Herrmann}
}
@conference {aravkin2011ICIAMrfwiu,
	title = {Robust FWI using Student{\textquoteright}s t-distribution},
	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, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ICIAM/2011/aravkin2011ICIAMrfwiu/aravkin2011ICIAMrfwiu.pdf},
	author = {Aleksandr Y. Aravkin and Tristan van Leeuwen and Felix J. Herrmann}
}
@conference {aravkin2011ICIAMspf,
	title = {Sparsity promoting formulations and algorithms for FWI},
	booktitle = {ICIAM},
	year = {2011},
	note = {Presented at AMP Medical and Seismic Imaging, 2011, Vancouver BC.},
	month = {07},
	publisher = {ICIAM 2011},
	organization = {ICIAM 2011},
	abstract = {Full Waveform Inversion (FWI) is a computational procedure to extract medium parameters from seismic data. FWI is typically formulated as a nonlinear least squares optimization problem, and various regularization techniques are used to guide the optimization because the problem is ill-posed. We propose a novel sparse regularization which exploits the ability of curvelets to efficiently represent geophysical images. We then formulate a corresponding sparsity promoting constrained optimization problem, which we solve using an open source algorithm. The techniques are applicable to any inverse problem where sparsity modeling is appropriate. We demonstrate the efficacy of the formulation on a toy example (stylized cross-well experiment) and on a realistic Seismic example (partial Marmoussi model). We also discuss the tradeoff between model fit and sparsity promotion, with a view to extend existing techniques for linear inverse problems to the case where the forward model is nonlinear.},
	keywords = {Full-waveform inversion, Optimization, Presentation, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ICIAM/2011/aravkin2011ICIAMspf/aravkin2011ICIAMspf_pres.pdf},
	author = {Aleksandr Y. Aravkin and Tristan van Leeuwen and James V. Burke and Felix J. Herrmann}
}
@conference {aravkin2011EAGEspfwi,
	title = {Sparsity promoting formulations and algorithms for FWI},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2011},
	keywords = {EAGE, Full-waveform inversion},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2011/Aravkin2011EAGEspfwi/Aravkin2011EAGEspfwi.pdf},
	author = {Aleksandr Y. Aravkin and Tristan van Leeuwen and James V. Burke and Felix J. Herrmann}
}
@conference {herrmann2003SPIEmsa,
	title = {Multifractional splines: application to seismic imaging},
	booktitle = {Proceedings of SPIE Technical Conference on Wavelets: Applications in Signal and Image Processing X},
	volume = {5207},
	year = {2003},
	pages = {240-258},
	publisher = {SPIE},
	organization = {SPIE},
	abstract = {Seismic imaging commits itself to locating singularities in the elastic properties of the Earth{\textquoteright}s subsurface. Using the high-frequency ray-Born approximation for scattering from non-intersecting smooth interfaces, seismic data can be represented by a generalized Radon transform mapping the singularities in the medium to seismic data. Even though seismic data are bandwidth limited, signatures of the singularities in the medium carry through this transform and its inverse and this mapping property presents us with the possibility to develop new imaging techniques that preserve and characterize the singularities from incomplete, bandwidth-limited and noisy data. In this paper we propose a non-adaptive Curvelet/Contourlet technique to image and preserve the singularities and a data-adaptive Matching Pursuit method to characterize these imaged singularities by Multi-fractional Splines. This first technique borrows from the ideas within the Wavelet-Vaguelette/Quasi-SVD approach. We use the almost diagonalization of the scattering operator to approximately compensate for (i) the coloring of the noise and hence facilitate estimation; (ii) the normal operator itself. Results of applying these techniques to seismic imaging are encouraging although many open fundamental questions remain.},
	keywords = {Presentation, SLIM},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SPIE/2003/herrmann2003SPIEmsa/herrmann2003SPIEmsa.pdf},
	author = {Felix J. Herrmann},
	editor = {Michael A. Unser and Akram Aldroubi and Andrew F. Laine}
}
