@conference {orozco2024IPMSnwi,
	title = {Neural wave-based imaging with amortized uncertainty quantification},
	booktitle = {Inverse Problems: Modelling and Simulation},
	year = {2024},
	month = {5},
	keywords = {Amortized Variational Inference, Bayesian inference, CIG, conditional normalizing flows, deep learning, FWI, Imaging, Inverse problems, IPMS, MVA, RTM, Summary Statistics, Uncertainty quantification, WISE},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/IPMS/2024/orozco2024IPMSnwi},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/IPMS/2024/orozco2024IPMSnwi/herrmann2024.pdf},
	author = {Rafael Orozco and Ziyi Yin and Ali Siahkoohi and Mathias Louboutin and Felix J. Herrmann}
}
@conference {orozco2024EAGEnfb,
	title = {Normalizing Flows for Bayesian Experimental Design in Imaging Applications},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2024},
	month = {6},
	abstract = {Neural density estimators, such as invertible normalizing flows, are capable of estimating the Bayesian posterior distribution in a variety of imaging problems, including medical MRI and seismic imaging/monitoring. So far, few works explore the possibility to make explicit use of probabilistic information contained within the full Bayesian solution of these inverse problems. During our talk, we investigate how a simple modification to the training objective of conditional normalizing flows allows for Bayesian experimental design without modifying the normalizing flow{\textquoteright}s neural architecture itself. By establishing a key relationship between the expected information gain (EIG) and the maximum-likelihood, attained during the training of normalizing flows, we show that optimal experimental design can be achieved. During our talk, we first verify, on a stylized problem, that our method indeed maximizes the expected information gain, followed by demonstrating the advocacy of our methodology on large-scale medical and seismic problems.},
	keywords = {Bayesian inference, deep learning, expected information gain, Normalizing flows, optimal experimental design, Uncertainty quantification},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2024/orozco2024EAGEnfb},
	author = {Rafael Orozco and Abhinav Prakash Gahlot and Peng Chen and Mathias Louboutin and Felix J. Herrmann}
}
@conference {orozco2024SIAMboed,
	title = {Normalizing Flows for Bayesian Experimental Design in Imaging Applications},
	booktitle = {SIAM Conference on Uncertainty Quantification},
	year = {2024},
	note = {(SIAM UQ, Trieste)},
	month = {03},
	abstract = {Neural density estimators such as normalizing flows have shown promise for estimation of the Bayesian posterior in a variety of imaging problems. Few works have explored how to practically exploit the probabilistic information contained in the full Bayesian solution of the inverse problem. Here we explore a simple modification to conditional normalizing flow training that enables Bayesian experimental design without modifying existing architectures. Based on a relationship between the expected information gain and maximum-likelihood training of normalizing flows, we show that experimental design can be achieved with the same training objective. We first verify that our method maximizes the expected information gain using a stylized problem. Then, we demonstrate our method can solve imaging problems in large scale medical and seismic applications.},
	keywords = {expected information gain, Normalizing flows, optimal experimental design, SIAM},
	doi = {10.48550/arXiv.2402.18337},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SIAMUQ/2024/orozco2024SIAMUQboed},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SIAMUQ/2024/orozco2024SIAMUQboed/orozco2024SIAMUQboed.pdf},
	author = {Rafael Orozco and Abhinav Prakash Gahlot and Peng Chen and Mathias Louboutin and Felix J. Herrmann}
}
@conference {herrmann2023SIAMCSEtns,
	title = {The Next Step: Interoperable Domain-Specific Programming},
	booktitle = {SIAM Conference on Computational Science and Engineering},
	year = {2023},
	note = {(SIAM CSE, Amsterdam)},
	month = {02},
	abstract = {Even though domain-specific programming approaches allow for readable, scalable, and maintainable software without sacrificing performance, the new paradigm of learned physics-informed models calls for an interdisciplinary approach typically involving multiple domain-specific languages. Take for example the problem of inverting for the fluid-flow properties from time-lapse seismic data, which entails domain-specific programming on the intersection of wave simulators, matrix-free linear algebra, learned neural surrogates for two-phase flow, and prior and posterior distributions for the fluid-flow properties. While domain-specific solutions exist for each of these sub-disciplines, integrating these approaches {\textendash} which may involve different programming languages {\textendash} into a single coupled scalable inversion framework that supports algorithmic differentiation can be a challenge. However, we show that challenges like this can be met when working with proper abstractions. In our inversion example, this involves math-inspired symbolic abstractions for numerical solutions of the wave equation (Devito), matrix-free implementations for its Jacobians (JUDI.jl), abstractions for Automatic Differentiation (ChainRules.jl), and homegrown implementations for conditional Invertible Neural Networks (InvertibleNetworks.jl) and Fourier Neural Operators (ParametricOperators.jl).},
	keywords = {algorithms, deep learning, devito, end-to-end, Fourier neural operators, GCS, JUDI, Jutul, SIAM, software, workshop},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SIAMCSE/2023/herrmann2023SIAMCSEtns/index.html},
	author = {Felix J. Herrmann and Mathias Louboutin and Thomas J. Grady II and Ziyi Yin and Rishi Khan}
}
@conference {peters2019SEGans,
	title = {A numerical solver for least-squares sub-problems in 3D wavefield reconstruction inversion and related problem formulations},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2019},
	note = {(SEG, San Antonio)},
	month = {09},
	pages = {1536-1540},
	abstract = {Recent years saw a surge of interest in seismic waveform inversion approaches based on quadratic-penalty or augmented-Lagrangian methods, including Wavefield Reconstruction Inversion. These methods typically need to solve a least-squares sub-problem that contains a discretization of the Helmholtz equation. Memory requirements for direct solvers are often prohibitively large in three dimensions, and this limited the examples in the literature to two dimensions. We present an algorithm that uses iterative Helmholtz solvers as a black-box to solve the least-squares problem corresponding to 3D grids. This algorithm enables Wavefield Reconstruction Inversion and related formulations, in three dimensions. Our new algorithm also includes a root-finding method to convert a penalty into a constraint on the data-misfit without additional computational cost, by reusing precomputed quantities. Numerical experiments show that the cost of parallel communication and other computations are small compared to the main cost of solving one Helmholtz problem per source and one per receiver.},
	keywords = {Full-waveform inversion, least-squares, numerical linear algebra, penalty method, quadratic constraints, root finding, SEG},
	doi = {10.1190/segam2019-3216638.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2019/peters2019SEGans/peters2019SEGans.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2019/peters2019SEGans/peters2019SEGans_pres.pdf},
	author = {Bas Peters and Felix J. Herrmann}
}
@conference {oghenekohwo2015PIMSntc,
	title = {A new take on compressive time-lapse seismic acquisition, imaging and inversion},
	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 = {Compressive sensing (CS), a sampling paradigm that is changing how data is acquired and processed in many fields, has attracted tremendous attention in recent years. In exploration seismology, particularly in 3D seismic technology, there have been published work on the application of CS to seismic data acquisition, processing and inversion. However, very little is known about the implications of CS in time-lapse (4D) seismic. Leveraging ideas from distributed CS where more than one signal is acquired and processed, as in 4D seismic, we present a joint recovery algorithm for time-lapse seismic data acquisition, imaging and inversion. The algorithm exploits the common information in the data sets (baseline and monitor) as part of an inversion procedure. Application of our algorithm to realistic synthetic examples highlight its advantages over other conventional approaches for processing time-lapse seismic data.},
	keywords = {PIMS, time-lapse, workshop},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/PIMS/2015/oghenekohwo2015PIMSntc/oghenekohwo2015PIMSntc_pres.pdf},
	author = {Felix Oghenekohwo and Felix J. Herrmann}
}
@conference {leeuwen2014EAGEntf,
	title = {A new take on FWI: wavefield reconstruction inversion},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2014},
	month = {06},
	abstract = {We discuss a recently proposed novel method for waveform inversion: Wavefield Reconstruction Inversion (WRI). As opposed to conventional FWI {\textendash} which attempts to minimize the error between observed and predicted data obtained by solving a wave equation {\textendash} WRI reconstructs a wave-field from the data and extracts a model-update from this wavefield by minimizing the wave-equation residual. The method does not require explicit computation of an adjoint wavefield as all the necessary information is contained in the reconstructed wavefield. We show how the corresponding model updates can be interpreted physically analogously to the conventional imaging-condition-based approach.},
	keywords = {EAGE, Full-waveform inversion, Optimization, penalty method, Wavefield Reconstruction Inversion},
	doi = {10.3997/2214-4609.20140703},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2014/leeuwen2014EAGEntf/leeuwen2014EAGEntf.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2014/leeuwen2014EAGEntf/leeuwen2014EAGEntf_pres.pdf},
	author = {Tristan van Leeuwen and Felix J. Herrmann and Bas Peters}
}
@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 {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 {hennenfent2008SINBADnii2,
	title = {New insights into one-norm solvers from the Pareto curve},
	booktitle = {SINBAD},
	year = {2008},
	abstract = {Several geophysical ill-posed inverse problems are successfully solved by promoting sparsity using one-norm regularization. The practicality of this approach depends on the effectiveness of the one-norm solver used and on its robustness under limited number of iterations. We propose an approach to understand the behavior and evaluate the performance of one-norm solvers. The technique consists of tracking on a graph the data misfit versus the one norm of successive iterates. By comparing the solution paths to the Pareto curve, we are able to assess the performance of the solvers and the quality of the solutions. Such an assessment is particularly relevant given the renewed interest in one-norm regularization.},
	keywords = {SINBAD, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2008/hennenfent2008SINBADnii2/hennenfent2008SINBADnii2.pdf},
	author = {Gilles Hennenfent}
}
@conference {yarham2007EAGEnsw,
	title = {Nonlinear surface wave prediction and separation},
	booktitle = {EAGE 2007},
	year = {2007},
	abstract = {Removal of surface waves is an integral step in seismic processing. There are many standard techniques for removal of this type of coherent noise, such as f-k filtering, but these methods are not always effective. One of the common problems with removal of surface waves is that they tend to be aliased in the frequency domain. This can make removal difficult and affect the frequency content of the reflector signals, as this signals will not be completely separated. As seen in (Hennenfent, G. and F. Herrmann, 2006, Application of stable signal recovery to seismic interpolation) interpolation can be used effectively to resample the seismic record thus dealiasing the surface waves. This separates the signals in the frequency domain allowing for a more precise and complete removal. The use of this technique with curvelet based surface wave predictions and an iterative L1 separation scheme can be used to remove surface waves from shot records more completely that with standard techniques.},
	keywords = {EAGE, Presentation, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2007/yarham2007EAGEnsw/yarham2007EAGEnsw.pdf},
	author = {Carson Yarham}
}
@conference {sastry2007SINBADnor,
	title = {Norm-one recovery from irregular sampled data},
	booktitle = {SINBAD 2007},
	year = {2007},
	abstract = {Seismic traces are sampled irregularly and insufficiently due to practical and economical limitations. The use of such data in seismic imaging results in image artifacts and poor spatial resolution. Therefore, before being used, the measurements are to be interpolated onto a regular grid. One of the methods achieving this objective is based on the Fourier reconstruction, which deals with the under-determined system of equations. The recent pursuit techniques (namely, basis pursuit, matching pursuit etc) admit certain promising features such as faster and simpler implementation even in large scale settings. The presentation discusses the application of the pursuit algorithms to the Fourier-based interpolation problem for the signals that have sparse Fourier spectra. In particular, the objective of the presentation includes: 1). studying the performance of the algorithm if, and how far, the measurement coordinates can be shifted from uniform distribution on the continuous interval. 2). studying what could be the allowable misplacement in the measurement coordinates that does not alter the quality of the reconstruction process},
	keywords = {Presentation, SINBAD, SLIM},
	author = {Challa S. Sastry}
}
@conference {hennenfent2006SINBADtnf,
	title = {The Nonuniform Fast Discrete Curvelet Transform (NFDCT)},
	booktitle = {SINBAD 2006},
	year = {2006},
	abstract = {The authors present an extension of the fast discrete curvelet transform (FDCT) to nonuniformly sampled data. This extension not only restores curvelet compression rates for nonuniformly sampled data but also removes noise and maps the data to a regular grid.},
	keywords = {Presentation, SINBAD, SLIM},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2006/hennenfent2006SINBADtnf/hennenfent2006SINBADtnf.pdf},
	author = {Gilles Hennenfent}
}
@conference {herrmann2005CSEGnld,
	title = {Non-linear data continuation with redundant frames},
	booktitle = {CSEG Annual Conference Proceedings},
	year = {2005},
	month = {05},
	publisher = {CSEG},
	organization = {CSEG},
	abstract = {We propose an efficient iterative data interpolation method using continuity along reflectors in seismic images via curvelet and discrete cosine transforms. The curvelet transform is a new multiscale transform that provides sparse representations for images that comprise smooth objects separated by piece-wise smooth discontinuities (e.g. seismic images). The advantage of using curvelets is that these frames are sparse for high-frequency caustic-free solutions of the wave-equation. Since we are dealing with less than ideal data (e.g. bandwidth-limited), we compliment the curvelet frames with the discrete cosine transform. The latter is motivated by the successful data continuation with the discrete Fourier transform. By choosing generic basis functions we circumvent the necessity to make parametric assumptions (e.g., through linear/parabolic Radon or demigration) regarding the shape of events in seismic data. Synthetic and real data examples demonstrate that our algorithm provides interpolated traces that accurately reproduce the wavelet shape as well as the AVO behavior along events in shot gathers.},
	keywords = {Presentation, SLIM},
	url = {http://www.cseg.ca/assets/files/resources/abstracts/2005/101S0201-Herrmann_F_Non_Linear_Data_Continuation.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/CSEG/2005/Herrmann05CSEGnld/Herrmann05CSEGnld_pres.pdf},
	author = {Felix J. Herrmann and Gilles Hennenfent}
}
@conference {herrmann2005CSEGnlr,
	title = {Non-linear regularization in seismic imaging},
	booktitle = {CSEG Annual Conference Proceedings},
	year = {2005},
	month = {05},
	publisher = {CSEG},
	organization = {CSEG},
	abstract = {Two complementary solution strategies to the least-squares imaging problem with sparseness \& continuity continuity constraints are proposed. The applied formalism explores the sparseness of curvelets coefficients of the reflectivity and their invariance under the demigration-migration operator. We achieve the solution by jointly minimizing a weighted l1-norm on the curvelet coefficients and an anisotropic difussion or total variation norm on the imaged reflectivity model. The l1-norm exploits the sparsenss of the reflectivity in the curvelet domain whereas the anisotropic norm enhances the continuity along the reflections while removing artifacts residing in between reflectors. While the two optimization methods (convex versus non-convex) share the same type of regularization, they differ in flexibility how to handle additional constraints on the coefficients of the imaged reflectivity and in computational expense. A brief sketch of the theory is provided along with a number of synthetic examples.},
	keywords = {Presentation, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/CSEG/2005/Herrmann05CSEGnlr/Herrmann05CSEGnlr.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/CSEG/2005/Herrmann05CSEGnlr/Herrmann05CSEGnlr_pres.pdf},
	author = {Felix J. Herrmann and Peyman P. Moghaddam}
}
