@conference {park2025IMAGErdn,
	title = {A reduced-order derivative-informed neural operator for subsurface fluid-flow},
	booktitle = {International Meeting for Applied Geoscience and Energy},
	year = {2025},
	month = {8},
	keywords = {deep learning, derivatve-informed, Fisher-information, FNO, GCS, IMAGE, Inverse problems, permeability, reduced-order modeling, reservoir simulation, SEG, Uncertainty quantification},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2025/park2025IMAGErdn/IMAGE_Jayjay.pdf},
	author = {Jeongjin Park and Grant Bruer and Huseyin Tuna Erdinc and Abhinav Prakash Gahlot and Felix J. Herrmann}
}
@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 {herrmann2016SEGWScvp,
	title = {Constraints versus penalties for edge-preserving full-waveform inversion},
	booktitle = {SEG Workshop on Where are we heading with FWI; Dallas},
	year = {2016},
	note = {(SEG Workshop, Dallas)},
	month = {10},
	abstract = {Full-waveform inversion is arguably one of the most challenging inverse problems out there. For this reason, it is remarkable that so much progress has been made over the years. While we can look back with confidence at an increasing number of success stories, the application of full-waveform technologies to deeper targets and complex geologies remains challenging making it an expensive proposition. So far, efforts to overcome some of these challenges have mainly been directed towards extended formulations such as Symes{\textquoteright}s and Biondi{\textquoteright}s subsurface offset/rau-parameters, Mike Warner/Lluis Guasch{\textquoteright}s Wiener filters in their adaptive full-waveform inversion, and van Leeuwen/Herrmann{\textquoteright}s Wavefield Reconstruction Inversion, where the wave-equation appears as a least-squares term in the objective. While this recent work {\textendash}- including other approaches such as different data-misfit objective functions {\textendash}- has removed some of FWI cycle-skip issues it does, with the exception perhaps of box constraints, not incorporate rudimentary information on the geology on the unknown model. Following recent work by Ernie Esser, we present a general framework how to incorporate this type of information in the form of constraints. Contrary to Tikhonov-like regularization or gradient filtering, imposing constraints has several key advantages, namely (i) they can be expressed in simple human understandable terms such as lower and upper limits for the velocity without the necessity to introduce additional parameters and assumptions on the probability distribution of the model; (ii) they can consist of several constraints, e.g. box and smoothness constraints, as long as their intersection in non-empty; and (iii) they can be imposed in a separate inner loop, which does not affect gradients and Hessians of full-waveform inversion itself. By means of examples, we will demonstrate the advocacy of constrained formulations as they apply to full-waveform inversion. This is joint work with Bas Peters.},
	keywords = {Full-waveform inversion, SEG, workshop},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/herrmann2016SEGWScvp/herrmann2016SEGWScvp_pres.pdf},
	author = {Felix J. Herrmann and Bas Peters}
}
@conference {peters2016SEGprs,
	title = {Parallel reformulation of the sequential adjoint-state method},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2016},
	note = {(SEG, Dallas)},
	month = {10},
	pages = {1411-1415},
	abstract = {This abstract is about reducing the total computation time to solve full-waveform inversion problems. A common problem formulation is a nonlinear least-squares problem where the gradient is computed using the adjoint-state algorithm. This formulation offers parallelism for the computation of gradients for different sources and frequencies. The adjoint-state algorithm itself is sequential however and this is a limiting factor when a lot of compute nodes are available and only a few wavefields need to be computed. This situation occurs when stochastic optimization strategies are used to minimize the objective function. We present a parallel reformulation of the sequential adjoint-state algorithm, which allows the forward- and adjoint wavefields to be computed in parallel. Both algorithms are mathematically equivalent but the parallel version is twice as fast in run time. An important characteristic of the proposed algorithm is that one wavefield needs to be computed per source and one per receiver. These fields can be used to apply the (inverse) Gauss-Newton Hessian to a vector without recomputing wavefields. A 2D example shows that good full-waveform inversion results are obtained, even when a small number of sources and receivers is used.},
	keywords = {FWI, numerical linear algebra, parallel computing, PDE-constrained optimization, SEG},
	doi = {10.1190/segam2016-13966771.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/peters2016SEGprs/peters2016SEGprs.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/peters2016SEGprs/peters2016SEGprs_pres.pdf},
	author = {Bas Peters and Tristan van Leeuwen and Felix J. Herrmann}
}
@techreport {peters2015SEGrwi,
	title = {Regularizing waveform inversion by projections onto convex sets {\textendash}- application to the 2D Chevron 2014 synthetic blind-test dataset},
	number = {TR-EOAS-2015-7},
	year = {2015},
	month = {06},
	institution = {UBC},
	abstract = {A framework is proposed for regularizing the waveform inversion problem by projections onto intersections of convex sets. Multiple pieces of prior information about the geology are represented by multiple convex sets, for example limits on the velocity or minimum smoothness conditions on the model. The data-misfit is then minimized, such that the estimated model is always in the intersection of the convex sets. Therefore, it is clear what properties the estimated model will have at each iteration. This approach does not require any quadratic penalties to be used and thus avoids the known problems and limitations of those types of penalties. It is shown that by formulating waveform inversion as a constrained problem, regularization ideas such as Tikhonov regularization and gradient filtering can be incorporated into one framework. The algorithm is generally applicable, in the sense that it works with any (differentiable) objective function and does not require significant additional computation. The method is demonstrated on the inversion of the 2D marine isotropic elastic synthetic seismic benchmark by Chevron using an acoustic modeling code. To highlight the effect of the projections, we apply no data pre-processing.},
	keywords = {blind-test, projection, regularization, SEG, Wavefield Reconstruction Inversion, waveform inversion},
	url = {https://slim.gatech.edu/Publications/Public/TechReport/2015/peters2015SEGrwi/peters2015SEGrwi.html},
	author = {Bas Peters and Zhilong Fang and Brendan R. Smithyman and Felix J. Herrmann}
}
@conference {peters2015SEGWSsew,
	title = {Simultaneous estimation of wavefields and medium parameters - reduced-space versus full-space waveform inversion},
	booktitle = {SEG Workshop on The Limit of FWI in Subsurface Parameter Recovery; New Orleans},
	year = {2015},
	note = {(SEG, New Orleans)},
	month = {10},
	keywords = {full-space, PDE-constrained optimization, quadratic penalty method, SEG, waveform inversion, workshop},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2015/peters2015SEGWSsew/peters2015SEGWSsew_pres.pdf},
	author = {Bas Peters and Felix J. Herrmann}
}
@conference {peters2014SEGsrh,
	title = {A sparse reduced Hessian approximation for multi-parameter wavefield reconstruction inversion},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2014},
	note = {(SEG)},
	month = {10},
	pages = {1206-1210},
	abstract = {Multi-Parameter full-waveform inversion is a challenging problem, because the unknown parameters appear in the same wave equation and the magnitude of the parameters can vary many orders of magnitude. This makes accurate estimation of multiple-parameters very difficult. To mitigate the problems, sequential strategies, regularization methods and scalings of gradients and quasi-Newton Hessians have been proposed. All of these require design, fine-tuning and adaptation to different waveform inversion problems. We propose to use a sparse approximation to the Hessian derived from a penalty-formulation of the objective function. Sparseness allows to have the Hessian in memory and compute update directions at very low cost. This results in decent reconstruction of the multiple parameters at very low additional memory and computational expense.},
	keywords = {Full-waveform inversion, Hessian, Optimization, penalty method, SEG},
	doi = {http://dx.doi.org/10.1190/segam2014-1667.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2014/peters2014SEGsrh/peters2014SEGsrh.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2014/peters2014SEGsrh/peters2014SEGsrh_pres.pdf},
	author = {Bas Peters and Felix J. Herrmann}
}
@conference {herrmann2014SEGWS,
	title = {Wavefield Reconstruction Inversion (WRI) {\textendash} a new take on wave-equation based inversion},
	booktitle = {SEG Workshop on Full Waveform Inversion - Elastic Approaches and Issues with Anisotropy, Nonshallow Inversion, Poor Starting Model; Denver},
	year = {2014},
	note = {(SEG Workshop, Denver)},
	month = {10},
	keywords = {FWI, SEG, workshop, WRI},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2014/herrmann2014SEGWSwri/herrmann2014SEGWSwri_pres.pdf},
	author = {Felix J. Herrmann and Ernie Esser and Tristan van Leeuwen and Bas Peters}
}
