@conference {herrmann2009IAPcsisa,
	title = {Compressive seismic imaging with simultaneous acquisition},
	booktitle = {IAP},
	year = {2009},
	note = {Presented at the IAP meeting, Vienna},
	abstract = {The shear size of seismic data volumes forms one of the major impediments for the inversion of seismic data. Turning forward modeling and inversion into a compressive sensing (CS) problem - where simulated data are recovered from a relatively small number of independent sources - can effectively mitigate this high-cost impediment. Our key contribution lies in the design of a sub-sampling operator that commutes with the time-harmonic Helmholtz system. As in compressive sensing, this leads to a reduction of simulation cost. This reduction is commensurate with the transform-domain sparsity of the solution., implying that computational costs are no longer determined by the size of the discretization but by transform-domain sparsity of the solution of the CS problem that recovers the data. The combination of this sub-sampling strategy with our recent work on preconditioned implicit solvers for the time-harmonic Helmholtz equation provides a viable alternative to full-waveform inversion schemes based on explicit time-domain finite-difference methods.},
	keywords = {Presentation},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/AIP/2009/herrmann2009IAPcsisa/Herrmann09AIP1.pdf},
	author = {Felix J. Herrmann and Yogi A. Erlangga and Tim T.Y. Lin}
}
@conference {herrmann2009EAGEcsa,
	title = {Compressive sensing applied to full-waveform inversion},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2009},
	month = {06},
	abstract = {With the recent resurgence of full-waveform inversion, the computational cost of solving forward modeling problems has become{\textendash}-aside from issues with non-uniqueness{\textendash}-one of the major impediments withstanding successful application of this technology to industry-size data volumes. To overcome this impediment, we argue that further improvements in this area will depend on a problem formulation with a computational complexity that is no longer strictly determined by the size of the discretization but by transform-domain sparsity of its solution. In this new paradigm, we bring computational costs in par with our ability to compress seismic data and images. This premise is related to two recent developments. First, there is the new field of compressive sensing (CS in short throughout the paper, Cand{\textquoteleft}es et al., 2006; Donoho, 2006){\textendash}-where the argument is made, and rigorously proven, that compressible signals can be recovered from severely sub-Nyquist sampling by solving a sparsity promoting program. Second, there is in the seismic community the recent resurgence of simultaneous-source acquisition (Beasley, 2008; Krohn and Neelamani, 2008; Herrmann et al., 2009; Berkhout, 2008; Neelamani et al., 2008), and continuing efforts to reduce the cost of seismic modeling, imaging, and inversion through phase encoding of simultaneous sources (Morton and Ober, 1998; Romero et al., 2000; Krohn and Neelamani, 2008; Herrmann et al., 2009), removal of subsets of angular frequencies (Sirgue and Pratt, 2004; Mulder and Plessix, 2004; Lin et al., 2008) or plane waves (Vigh and Starr, 2008). By using CS principles, we remove sub-sampling interferences asocciated with these approaches through a combination of exploiting transform-domain sparsity, properties of certain sub-sampling schemes, and the existence of sparsity promoting solvers.},
	keywords = {EAGE},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2009/Herrmann09EAGEcsa/Herrmann09EAGEcsa.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2009/Herrmann09EAGEcsa/Herrmann09EAGEcsa_pres.pdf},
	url2 = {http://earthdoc.eage.org/publication/publicationdetails/?publication=23961},
	author = {Felix J. Herrmann and Yogi A. Erlangga and Tim T.Y. Lin}
}
@conference {lin2009SEGcsf,
	title = {Compressive simultaneous full-waveform simulation},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {28},
	year = {2009},
	month = {10},
	pages = {2577-2581},
	publisher = {SEG},
	organization = {SEG},
	abstract = {The fact that the computational complexity of wavefield simulation is proportional to the size of the discretized model and acquisition geometry, and not to the complexity of the simulated wavefield, is a major impediment within seismic imaging. By turning simulation into a compressive sensing problem{\textendash}where simulated data is recovered from a relatively small number of independent simultaneous sources{\textendash}we remove this impediment by showing that compressively sampling a simulation is equivalent to compressively sampling the sources, followed by solving a reduced system. As in compressive sensing, this allows for a reduction in sampling rate and hence in simulation costs. We demonstrate this principle for the time-harmonic Helmholtz solver. The solution is computed by inverting the reduced system, followed by a recovery of the full wavefield with a sparsity promoting program. Depending on the wavefield{\textquoteright}s sparsity, this approach can lead to a significant cost reduction, in particular when combined with the implicit preconditioned Helmholtz solver, which is known to converge even for decreasing mesh sizes and increasing angular frequencies. These properties make our scheme a viable alternative to explicit time-domain finite-difference.},
	keywords = {Presentation, SEG},
	doi = {10.1190/1.3255381},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2009/Lin09SEGcsf/Lin09SEGcsf.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2009/Lin09SEGcsf/Lin09SEGcsf_pres.pdf},
	author = {Tim T.Y. Lin and Yogi A. Erlangga and Felix J. Herrmann}
}
@article {herrmann2009GEOPcsf,
	title = {Compressive simultaneous full-waveform simulation},
	journal = {Geophysics},
	volume = {74},
	number = {4},
	year = {2009},
	month = {07-08},
	pages = {A35-A40},
	abstract = {The fact that computational complexity of wavefield simulation is proportional to the size of the discretized model and acquisition geometry, and not to the complexity of the simulated wavefield, is a major impediment within seismic imaging. By turning simulation into a compressive sensing problem{\textendash}-where simulated data is recovered from a relatively small number of independent simultaneous sources{\textendash}-we remove this impediment by showing that compressively sampling a simulation is equivalent to compressively sampling the sources, followed by solving a reduced system. As in compressive sensing, this allows for a reduction in sampling rate and hence in simulation costs. We demonstrate this principle for the time-harmonic Helmholtz solver. The solution is computed by inverting the reduced system, followed by a recovery of the full wavefield with a sparsity promoting program. Depending on the wavefield{\textquoteright}s sparsity, this approach can lead to significant cost reductions, in particular when combined with the implicit preconditioned Helmholtz solver, which is known to converge even for decreasing mesh sizes and increasing angular frequencies. These properties make our scheme a viable alternative to explicit time-domain finite-differences.},
	keywords = {Compressive Sensing, full-waveform, Modelling, SLIM},
	doi = {10.1190/1.3115122},
	url = {http://library.seg.org/doi/abs/10.1190/1.3115122},
	url2 = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2009/herrmann2009GEOPcsf/herrmann2009GEOPcsf.pdf},
	author = {Felix J. Herrmann and Yogi A. Erlangga and Tim T.Y. Lin}
}
@conference {herrmann2009SAMPTAcws,
	title = {Compressive-wavefield simulations},
	booktitle = {SAMPTA},
	year = {2009},
	publisher = {SAMPTA},
	organization = {SAMPTA},
	abstract = {Full-waveform inversion{\textquoteright}s high demand on computational resources forms, along with the non-uniqueness problem, the major impediment withstanding its widespread use on industrial-size datasets. Turning modeling and inversion into a compressive sensing problem{\textendash}-where simulated data are recovered from a relatively small number of independent simultaneous sources{\textendash}-can effectively mitigate this high-cost impediment. The key is in showing that we can design a sub-sampling operator that commutes with the time-harmonic Helmholtz system. As in compressive sensing, this leads to a reduction in simulation cost. Moreover, this reduction is commensurate with the transform-domain sparsity of the solution, implying that computational costs are no longer determined by the size of the discretization but by transform-domain sparsity of the solution of the CS problem which forms our data. The combination of this sub-sampling strategy with our recent work on implicit solvers for the Helmholtz equation provides a viable alternative to full-waveform inversion schemes based on explicit finite-difference methods.},
	keywords = {SAMPTA},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SAMPTA/2009/Herrmann09SAMPTAcws/Herrmann09SAMPTAcws.pdf},
	author = {Felix J. Herrmann and Yogi A. Erlangga and Tim T.Y. Lin}
}
@article {herrmann2009GEOPcbm,
	title = {Curvelet-based migration preconditioning and scaling},
	journal = {Geophysics},
	volume = {74},
	year = {2009},
	month = {07-08},
	pages = {A41-A46},
	abstract = {The extremely large size of typical seismic imaging problems has been one of the major stumbling blocks for iterative techniques to attain accurate migration amplitudes. These iterative methods are important because they complement theoretical approaches that are hampered by difficulties to control problems such as finite-acquisition aperture, source-receiver frequency response, and directivity. To solve these problems, we apply preconditioning, which significantly improves convergence of least-squares migration. We discuss different levels of preconditioning that range from corrections for the order of the migration operator to corrections for spherical spreading, and position and reflector-dip dependent amplitude errors. While the first two corrections correspond to simple scalings in the Fourier and physical domain, the third correction requires phase-space (space spanned by location and dip) scaling, which we carry out with curvelets. We show that our combined preconditioner leads to a significant improvement of the convergence of least-squares {\textquoteleft}wave-equation{\textquoteright} migration on a line from the SEG AA{\textquoteright} salt model.},
	keywords = {Imaging, migration, SLIM},
	doi = {10.1190/1.3124753},
	url = {http://geophysics.geoscienceworld.org/content/74/4/A41.abstract},
	url2 = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2009/herrmann2009GEOPcbm/herrmann2009GEOPcbm.pdf},
	author = {Felix J. Herrmann and Cody R. Brown and Yogi A. Erlangga and Peyman P. Moghaddam}
}
@conference {erlangga2009EAGEmwi,
	title = {Migration with implicit solvers for the time-harmonic Helmholtz equation},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2009},
	month = {06},
	abstract = {From the measured seismic data, the location and the amplitude of reflectors can be determined via a migration algorithm. Classically, following Claerbout{\textquoteright}s imaging principle [2], a reflector is located at the position where the source{\textquoteright}s forward-propagated wavefield correlates with the backward-propagated wavefield of the receiver data. Lailly and Tarantola later showed that this imaging principle is an instance of inverse problems, with the associated migration operator formulated via a least-squares functional; see [6, 12, 13]. Furthermore, they showed that the migrated image is associated with the gradient of this functional with respect to the image. If the solution of the least-squares functional is done iteratively, the correlation-based image coincides up to a constant with the first iteration of a gradient method. In practice, this migration is done either in the time domain or in the frequency domain. In the frequency-domain migration, the main bottleneck thus far, which renders its full implementation to large scale problems, is the lack of efficient solvers for computing wavefields. Robust direct methods easily run into excessive memory requirements as the size of the problem increases. On the other hand, iterative methods, which are less demanding in terms of memory, suffered from lack of convergence. During the past years, however, progress has been made in the development of an efficient iterative method [4, 3] for the frequency-domain wavefield computations. In this paper, we will show the significance of this method (called MKMG) in the context of the frequency-domain migration, where multi-shot-frequency wavefields (of order of 10,000 related wavefields) need to be computed.},
	keywords = {EAGE},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2009/Erlangga09EAGEmwi/Erlangga09EAGEmwi.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2009/Erlangga09EAGEmwi/Erlangga09EAGEmwi_pres.pdf},
	url2 = {http://earthdoc.eage.org/publication/publicationdetails/?publication=23955},
	author = {Yogi A. Erlangga and Felix J. Herrmann}
}
@conference {lin2009DELPHIrwi,
	title = {Randomized wavefield inversion},
	booktitle = {DELPHI},
	year = {2009},
	note = {Presented at the DELPHI meeting. The Hague},
	keywords = {Presentation},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/Delphi/2009/lin2009DELPHIrwi/lin2009DELPHIrwi.pdf},
	author = {Tim T.Y. Lin and Felix J. Herrmann and Yogi A. Erlangga}
}
@conference {erlangga2009SEGswi,
	title = {Seismic waveform inversion with Gauss-Newton-Krylov method},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {28},
	year = {2009},
	month = {10},
	pages = {2357-2361},
	publisher = {SEG},
	organization = {SEG},
	abstract = {This abstract discusses an implicit implementation of the Gauss-Newton method, used for the frequency-domain full-waveform inversion, where the inverse of the Hessian for the update is never formed explicitly. Instead, the inverse of the Hessian is computed approximately by a conjugate gradient (CG) method, which only requires the action of the Hessian on the CG search direction. This procedure avoids an excessive computer storage, usually needed for storing the Hessian, at the expense of extra computational work in CG. An effective preconditioner for the Hessian is important to improve the convergence of CG, and hence to reduce the overall computational work.},
	keywords = {SEG},
	doi = {10.1190/1.3255332},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2009/erlangga09SEGswi/erlangga09SEGswi.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2009/erlangga09SEGswi/erlangga09SEGswi_pres.pdf},
	author = {Yogi A. Erlangga and Felix J. Herrmann}
}
@conference {herrmann2008SIAMcsm,
	title = {Compressive sampling meets seismic imaging},
	booktitle = {SIAM},
	year = {2008},
	abstract = {Compressive sensing has led to fundamental new insights in the recovery of compressible signals from sub-Nyquist samplings. It is shown how jittered subsampling can be used to create favorable recovery conditions. Applications include mitigation of incomplete acquisitions and wavefield computations. While the former is a direct adaptation of compressive sampling, the latter application represents a new way of compressing wavefield extrapolation operators. Operators are not diagonalized but are compressively sampled reducing the computational costs.},
	keywords = {Presentation, SIAM, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SIAM/2008/herrmann2008SIAMcsm/herrmann2008SIAMcsm.pdf},
	author = {Felix J. Herrmann and Yogi A. Erlangga and Tim T.Y. Lin}
}
@conference {lin2008SEGiso,
	title = {Interpolating solutions of the Helmholtz equation with compressed sensing},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {27},
	year = {2008},
	month = {01},
	pages = {2122-2126},
	publisher = {SEG},
	organization = {SEG},
	abstract = {We present an algorithm which allows us to model wavefields with frequency-domain methods using a much smaller number of frequencies than that typically required by the classical sampling theory in order to obtain an alias-free result. The foundation of the algorithm is the recent results on the compressed sensing, which state that data can be successfully recovered from an incomplete measurement if the data is sufficiently sparse. Results from numerical experiment show that only 30\% of the total frequency spectrum is need to capture the full wavefield information when working in the hard 2D synthetic Marmousi model.},
	keywords = {Presentation, SEG, SLIM},
	doi = {10.1190/1.3059307},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2008/lin08SEGiso/lin08SEGiso.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2008/lin08SEGiso/lin08SEGiso.pdf},
	author = {Tim T.Y. Lin and Evgeniy Lebed and Yogi A. Erlangga and Felix J. Herrmann}
}
@conference {herrmann2008SINBADitc,
	title = {Introduction to compressive (wavefield) computation},
	booktitle = {SINBAD 2008},
	year = {2008},
	keywords = {Presentation, SINBAD, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2008/herrmann2008SINBADitc/herrmann2008SINBADitc.pdf},
	author = {Felix J. Herrmann and Yogi A. Erlangga and Tim T.Y. Lin and Cody R. Brown}
}
@conference {erlangga2008SEGaim,
	title = {An iterative multilevel method for computing wavefields in frequency-domain seismic inversion},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {27},
	year = {2008},
	month = {11},
	pages = {1957-1960},
	publisher = {SEG},
	organization = {SEG},
	abstract = {We describe an iterative multilevel method for solving linear systems representing forward modeling and back propagation of wavefields in frequency-domain seismic inversions. The workhorse of the method is the so-called multilevel Krylov method, applied to a multigrid-preconditioned linear system, and is called multigrid-multilevel Krylov (MKMG) method. Numerical experiments are presented for 2D Marmousi synthetic model for a range of frequencies. The convergence of the method is fast, and depends only mildly on frequency. The method can be considered as the first viable alternative to LU factorization, which is practically prohibitive for 3D seismic inversions.},
	keywords = {SEG, SLIM},
	doi = {10.1190/1.3059279},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2008/erlangga08SEGaim/erlangga08SEGaim.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2008/erlangga08SEGaim/erlangga08SEGaim_pres.pdf},
	author = {Yogi A. Erlangga and Felix J. Herrmann}
}
