@article {vanleeuwen2015GPWEMVA,
	title = {Enabling affordable omnidirectional subsurface extended image volumes via probing},
	journal = {Geophysical Prospecting},
	volume = {65},
	number = {2},
	year = {2017},
	note = {(Geophysical Prospecting)},
	month = {03},
	pages = {385-406},
	abstract = {Image gathers as a function of subsurface offset are an important tool for the inference of rock properties and velocity analysis in areas of complex geology. Traditionally, these gathers are thought of as multidimensional correlations of the source and receiver wavefields. The bottleneck in computing these gathers lies in the fact that one needs to store, compute, and correlate these wavefields for all shots in order to obtain the desired image gathers. Therefore, the image gathers are typically only computed for a limited number of subsurface points and for a limited range of subsurface offsets, which may cause problems in complex geological areas with large geologic dips. We overcome increasing computational and storage costs of extended image volumes by introducing a formulation that avoids explicit storage and removes the customary and expensive loop over shots found in conventional extended imaging. As a result, we end up with a matrix{\textendash}vector formulation from which different image gathers can be formed and with which amplitude-versus-angle and wave-equation migration velocity analysis can be performed without requiring prior information on the geologic dips. Aside from demonstrating the formation of two-way extended image gathers for different purposes and at greatly reduced costs, we also present a new approach to conduct automatic wave-equation-based migration-velocity analysis. Instead of focusing in particular offset directions and preselected subsets of subsurface points, our method focuses every subsurface point for all subsurface offset directions using a randomized probing technique. As a consequence, we obtain good velocity models at low cost for complex models without the need to provide information on the geologic dips.},
	keywords = {AVA, migration velocity analysis, Stochastic optimization},
	doi = {10.1111/1365-2478.12418},
	url = {https://slim.gatech.edu/Publications/Public/Journals/GeophysicalProspecting/2016/vanleeuwen2015GPWEMVA/vanleeuwen2015GPWEMVA.pdf},
	author = {Tristan van Leeuwen and Rajiv Kumar and Felix J. Herrmann}
}
@presentation {fang2016SINBADFeaq,
	title = {Efficient approach for quantifying uncertainty of wavefield reconstruction inversion},
	year = {2016},
	publisher = {SINBAD},
	abstract = {Due to the noisy measurements, sparse observations, uncertain forward models, and uncertain prior parameter information, significant uncertainty arises in the inverted velocity model of full-waveform inversion. The uncertain velocity model leads to the uncertainty in the correct positions of events in the migrated images, which subsequently influences the important decision-makings in the oil and gas exploration and production (E\&P) business such as the drilling decisions and economic evaluations. Therefore, understanding the uncertainty in the inverted velocity model is essential for the mitigation of the E\&P risks. A common approach to accounting for uncertainty in seismic inverse problems like full-waveform inversion is to describe the unknown parameters in probabilistic by the Bayesian inference. The Bayesian inference aims at incorporating the information from the data, physics, and one{\textquoteright}s prior knowledge to formulate a posterior distribution function to characterize the statistical properties of the unknown parameters. However, characterizing uncertainties for the large-scale full-waveform inversion in seismic exploration is prohibitively expensive since it requires thousands of function evaluations to sample the parameter space (e.g. via Markov chain Monte Carlo sampling). Apart from this, the existing methods of uncertainty quantification for full-waveform inversion is through applying the Bayesian inference to the standard method such as the adjoint-state method with the assumption that no uncertainty arises from the wave-equation. However, the wave-equation may not be able to reflect all the physics, hence, incorporating uncertainty from wave-equation might benefit the quantification of uncertainty. In this work, we propose to perform uncertainty quantification using the recently proposed wavefield reconstruction inversion technique, which gives an appropriate framework to tackle the uncertainty in wave-equation solve by adding the misfit of the wave-equation to the likelihood distribution and relaxing the wave-equation constraint. Moreover, we propose a computationally feasible approach to analyzing the posterior distribution of the wavefield reconstruction inversion method, which incorporates a Gaussian distribution that approximates the posterior distribution with an optimization-driven Gaussian simulator to sample the Gaussian distribution efficiently. Numerical examples show that comparing to the McMC method, our approach is able to produce comparable results with less computational cost. This is joint work with Curt da Silva and Rachel Kuske.},
	keywords = {Presentation, SINBAD, SINBADFALL2016, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/fang2016SINBADFeaq/fang2016SINBADFeaq.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/fang2016SINBADFeaq/fang2016SINBADFeaq.mov},
	author = {Zhilong Fang and Chia Ying Lee and Curt Da Silva and Tristan van Leeuwen and Felix J. Herrmann and Rachel Kuske}
}
@conference {kumar2014EAGEeia,
	title = {Extended images in action: efficient WEMVA via randomized probing},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2014},
	month = {06},
	abstract = {Image gathers as a function of subsurface offset are an important tool for velocity analysis in areas of complex geology. In this paper, we offer a new perspective on image gathers by organizing the extended image as a function of all subsurface offsets and all subsurface points into a matrix whose (i,j)^th entry captures the interaction between gridpoints i and j. For even small problems, it is infeasible to form and store this matrix. Instead, we propose an efficient algorithm to glean information from the image volume via efficient matrix-vector products. We illustrate how this can be used to construct objective functions for automatic MVA.},
	keywords = {EAGE, extended imaging, MVA, probing},
	doi = {10.3997/2214-4609.20141492},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2014/kumar2014EAGEeia/kumar2014EAGEeia.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2014/kumar2014EAGEeia/kumar2014EAGEeia_pres.pdf},
	author = {Rajiv Kumar and Tristan van Leeuwen and Felix J. Herrmann}
}
@presentation {kumar2014SINBADeia,
	title = {Extended images in action {\textendash}- image gather with surface-related multiples},
	year = {2014},
	publisher = {SINBAD},
	abstract = {Common image gathers are used in building velocity models, invert anisotropy parameters and to analyze reservoir attributes. Often primary reflections are used to form the image gathers and multiples are typically attenuated in processing. However, researchers have shown that, if correctly used, they can provide useful information about the subsurface during reverse-time migration. In this work, i will show how we can use multiples along with primaries to form the image-gathers.},
	keywords = {Presentation, SINBAD, SINBADFALL2014, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2014/Fall/kumar2014SINBADeia/kumar2014SINBADeia.pdf},
	author = {Rajiv Kumar and Ning Tu and Tristan van Leeuwen and Felix J. Herrmann}
}
@conference {kumar2013SEGMVA,
	title = {Efficient WEMVA using extended images},
	booktitle = {SEG Workshop on Advances in Model Building, Imaging, and FWI; Houston},
	year = {2013},
	month = {09},
	abstract = {Image gathers as a function of subsurface offset are an important tool for velocity analysis in areas of complex geology. Here, we offer a new perspective on image gathers by organizing the extended image as a function of all subsurface offsets and all subsurface points in to a matrix whose (i,j)th entry captures the interaction between gridpoints i and j. For even small problems, it is infeasible to form and store this matrix. Instead, we propose an efficient algorithm to glean information from the image volume via efficient matrix-vector products. We illustrate how this can be used to construct objective functions for automated MVA.},
	keywords = {inversion, MVA, poster, SEG, workshop},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2013/kumar2013SEGMVA/kumar2013SEGMVA_poster.pdf},
	author = {Rajiv Kumar and Tristan van Leeuwen and Felix J. Herrmann}
}
@presentation {peters2013SINBADepm,
	title = {Examples from the penalty-method},
	year = {2013},
	publisher = {SINBAD},
	abstract = {A novel penalty method for PDE-constrained optimization was recently proposed by van Leeuwen and Herrmann (2013). The conventional PDE-constrained optimization formulation used in seismic waveform inversion is based on calculating the gradient of the data misfit objective functional via the adjoint-state method, at the cost of two PDE solves. The penalty method requires only one solution of an overdetermined linear system, in the least-squares sense. In this talk some numerical properties of this linear system will be exposed. The penalty method for PDE-constrained optimization involves a parameter, balancing the data misfit and PDE-misfit parts of the objective functional. This talk will address how to select this very important parameter. Some examples will be shown in which the penalty method outperforms the conventional method in non-linear waveform inversion, as well as linearized seismic imaging by migration.},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/peters2013SINBADepm/peters2013SINBADepm.pdf},
	author = {Bas Peters and Tristan van Leeuwen and Felix J. Herrmann}
}
@presentation {kumar2013SINBADava,
	title = {Extended images in action: efficient AVA via probing},
	year = {2013},
	publisher = {SINBAD},
	abstract = {Common image gathers (CIG) are an important tool to perform AVA analysis in areas of complex geology. Unfortunately, it is prohibitively very expensive to compute these CIG for all the subsurface points. In this work, we present an efficient way to compute CIG for all subsurface offsets without explicitly calculating the source and receiver wavefields for all the sources. Because the CIG contain all possible subsurface offsets, we compute the angle-domain image gathers by selecting the subsurface offset that is aligned with the local geologic dip. We propose a method to compute the local dip information directly from common-image-point gathers. To assess the quality of the angle-domain common-image-points gathers we compute the angle-dependent reflectivity coefficients and compare them with theoretical reflectivity coefficients yielded by the (linearized) Zoeppritz equations for a few synthetic models.},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/kumar2013SINBADava/kumar2013SINBADava.pdf},
	author = {Rajiv Kumar and Tristan van Leeuwen and Felix J. Herrmann}
}
@presentation {kumar2013SINBADeia,
	title = {Extended images in action: efficient WEMVA via randomized probing},
	year = {2013},
	publisher = {SINBAD},
	abstract = {Extended images as a function of the full subsurface offset are an important tool to perform wave-equation based migration velocity analysis (WEMVA) in areas of complex geology. Unfortunately, computation \& storage of these extended images is prohibitively expensive. In this work, we present an efficient way to compute extended images for all subsurface offsets without explicitly calculating the source and receiver wavefields for all the sources. Instead, we calculate actions of extended image volumes on probing vectors that live in the image space. The probing can either be defined as vectors from the Dirac basis, which allows us to form the extended image at the location of the point diffractor, or they can be defined in terms of Gaussian noise. The latter corresponds to sources with random weights firing simultaneously at every grid point. We demonstrate that this probing leads to a computationally efficient implementation of WEMVA. This is joint work with Tristan van Leeuwen.},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/kumar2013SINBADeia/kumar2013SINBADeia.pdf},
	author = {Rajiv Kumar and Tristan van Leeuwen and Felix J. Herrmann}
}
@article {aravkin2012IPNuisance,
	title = {Estimating nuisance parameters in inverse problems},
	journal = {Inverse Problems},
	volume = {28},
	number = {11},
	year = {2012},
	month = {10},
	abstract = {Many inverse problems include nuisance parameters which, while not of direct interest, are required to recover primary parameters. Structure present in these problems allows efficient optimization strategies - a well known example is variable projection, where nonlinear least squares problems which are linear in some parameters can be very efficiently optimized. In this paper, we extend the idea of projecting out a subset over the variables to a broad class of maximum likelihood (ML) and maximum a posteriori likelihood (MAP) problems with nuisance parameters, such as variance or degrees of freedom. As a result, we are able to incorporate nuisance parameter estimation into large-scale constrained and unconstrained inverse problem formulations. We apply the approach to a variety of problems, including estimation of unknown variance parameters in the Gaussian model, degree of freedom (d.o.f.) parameter estimation in the context of robust inverse problems, automatic calibration, and optimal experimental design. Using numerical examples, we demonstrate improvement in recovery of primary parameters for several large- scale inverse problems. The proposed approach is compatible with a wide variety of algorithms and formulations, and its implementation requires only minor modifications to existing algorithms.},
	keywords = {full waveform inversion, students t, variance},
	doi = {10.1088/0266-5611/28/11/115016},
	url = {http://arxiv.org/abs/1206.6532},
	url2 = {https://slim.gatech.edu/Publications/Public/Journals/InverseProblems/2012/aravkin2012IPNuisance/aravkin2012IPNuisance.pdf},
	author = {Aleksandr Y. Aravkin and Tristan van Leeuwen}
}
