@presentation {petrenko2013SINBADaih,
	title = {Accelerating an iterative Helmholtz solver with FPGAs},
	year = {2013},
	publisher = {SINBAD},
	abstract = {Solution of the Helmholtz equation is the main computational burden of full-waveform inversion in the frequency domain. For this task we employ the CARP-CG algorithm (Gordon \& Gordon 2010), an iterative solver that preconditions the original Helmholtz system into an equivalent symmetric positive definite system and then applies the method of conjugate gradients. Forming the matrix for the new system is not necessary as its multiplicative action on a vector is implemented using a series of projections onto the rows of the original system. Our contribution is implementing CARP-CG for a host + accelerator (FPGA) computing environment. The computational paradigm is one of dataflow: vector and matrix elements are streamed from memory through the accelerator which applies the row projections. The advantage of an FPGA to process streams of data is that unless the algorithm is memory bandwidth limited, computation time is directly proportional to the amount of data. The complexity of the algorithm implemented on the FPGA is irrelevant since all the operations programmed onto the FPGA happen in the same clock tick. In contrast, on a CPU, more complex algorithms require more clock ticks as the instructions are executed sequentially, or with only a small amount of parallelism. Ongoing work porting the CARP-CG algorithm to the accelerator is presented.},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/petrenko2013SINBADaih/petrenko2013SINBADaih.pdf},
	author = {Art Petrenko and Tristan van Leeuwen and Felix J. Herrmann and Diego Oriato and Simon Tilbury}
}
@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}
}
@presentation {li2013SINBADgom,
	title = {Lessons learned from Chevron Gulf of Mexico data set},
	year = {2013},
	publisher = {SINBAD},
	abstract = {Chevron Gulf of Mexico data set is very challenging for FWI because of elastic phases, limited offset, lack of low frequencies and salt structure. To overcome these issue, we first use ray-based tomography on the hand-picked first breaks to generate initial model for FWI, and then we apply curvelet-denosing techniques to improve the poor signal-to-noise ratio of the observed data at low frequencies. Finally, Curvelet domain sparsity promoting Gauss-Newton FWI helps to suppress model space artifacts caused by elastic phases. This is joint work with Andrew J. Calvert, Ian Hanlon, Mostafa Javanmehri, Rajiv Kumar, Tristan van Leeuwen, Brendan R. Smithyman, Haneet Wason and Felix J. Herrmann},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/li2013SINBADgom/li2013SINBADgom.pdf},
	author = {Xiang Li and Andrew J. Calvert and Ian Hanlon and Mostafa Javanmehri and Rajiv Kumar and Tristan van Leeuwen and Brendan R. Smithyman and Haneet Wason and Felix J. Herrmann}
}
@presentation {li2013SINBADmsd,
	title = {Model-space versus data-space FWI with the acoustic wave equation},
	year = {2013},
	publisher = {SINBAD},
	abstract = {Inverting data with elastic phases using an acoustic wave equation can lead to erroneous results, especially when the number of iterations is too high, which may lead to over fitting the data. Several approaches have been proposed to address this issue. Most commonly, people data- filtering operations that are aimed to deemphasize the elastic phases in the data in favor of the acoustic phases. Examples of this approach are nested loops over offset range and Laplace parameters. In this presentation, we discuss two complementary optimization-driven methods where the minimization process decides adaptively which of the data or model components are consistent with the objective. Specifically, we compare the s t misfit function as the data-space alternative and curvelet-domain sparsity promotion as the model-space alternative. Application of these two methods to a realistic synthetic lead to comparable results that we believe can be improved by combining these two methods.},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/li2013SINBADmsd/li2013SINBADmsd.pdf},
	author = {Xiang Li and Anais Tamalet and Tristan van Leeuwen and Felix J. Herrmann}
}
@presentation {fang2013SINBADp3dfwi,
	title = {Parallel 3D FWI with simultaneous shots},
	year = {2013},
	publisher = {SINBAD},
	abstract = {In this work, we manage to build the workflow for the parallel 3D full-waveform inversion. In the forward simulation part, we generate the Helmholtz matrix in parallel and use parallel CARPCG to solve the Helmholtz equation. In the inversion process, simultaneous shots are used to reduce the computational costs. Additionally, we propose a method to select the number of simultaneous shots and tolerance of CARPCG dynamically, to reach a compromise between the computational costs and accuracy.},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/fang2013SINBADp3dfwi/fang2013SINBADp3dfwi.pdf},
	author = {Zhilong Fang and Tristan van Leeuwen and Felix J. Herrmann}
}
@presentation {vanleeuwen2013SINBADrpp,
	title = {Relaxing the physics: a penalty method for full-waveform inversion},
	year = {2013},
	publisher = {SINBAD},
	abstract = {Computationally efficient method of solving partial differential equation (PDE)-constrained optimization problems that occur in (geophysical) inversion problems. The method takes measured data from a physical system as input and minimizes an objective function that depends on an unknown model for the physical parameters, fields, and additional nuisance parameters such as the source function. The invention consists of a minimization procedure involving a cost functional comprising of a data-misfit term and a penalty term that measures how accurately the fields satisfy the PDE. The method is composed of two alternating steps, namely the solution of a system of equations forming the discretization of the data-augmented PDE, and the solution of physical model parameters from the PDE itself given the field that solves the data-augmented system and an estimate for the sources. Compared to all-at-once approaches to PDE-constrained optimization, there is no need to update and store the fields for all sources leading to significant memory savings. As in the all-at-once approach, the proposed method explores a larger search space and is therefore less sensitive to initial estimates for the physical model parameters. Contrary to the reduced formulation, the proposed method does not require the solution of an adjoint PDE, effectively halving the number of PDE solves and memory requirement. As in the reduced formulation, fields can be computed independently and aggregated, possibly in parallel.},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/vanleeuwen2013SINBADrpp/vanleeuwen2013SINBADrpp.pdf},
	author = {Tristan van Leeuwen and Felix J. Herrmann}
}
@presentation {vanleeuwen2013SINBADsda,
	title = {Solving the data-augmented wave equation},
	year = {2013},
	publisher = {SINBAD},
	abstract = {The recently proposed penalty method promises to mitigate some of the non-linearity inherent in full-waveform inversion by relaxing the requirement that the wave-equation needs to be solved exactly. The basic workflow of this new method is as follows; i) solve an overdetermined wave-equation (the data-augmented wave-equation), where the data serves as additional constraints for the wavefields, ii) compute the wavefield-residual by substituting this wavefield in the wave-equation, and iii) correlate the wavefield with the wavefield-residual to obtain a model-update. As opposed to the conventional workflow, no explicit adjoint solve is needed to compute the model-update. However, instead of solving a wave-equation, we need to solve a data-augmented wave-equation. In this talk we explore some of the challenges of solving this data-augmented wave-equation and review some possible solution strategies for both time and frequency-domain applications.},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/vanleeuwen2013SINBADsda/vanleeuwen2013SINBADsda.pdf},
	author = {Tristan van Leeuwen and Bas Peters and Felix J. Herrmann}
}
@presentation {fang2013SINBADsfwi,
	title = {Swift FWI},
	year = {2013},
	publisher = {SINBAD},
	abstract = {In 3D case, there is much more data and modeling is much more expensive. As a result, parallel computing is very important for 3D full-waveform inversion. Both domain decomposition and data decomposition need a large number of parallel computing resources. However, programs based on parallel MATLAB suffer from the limitation of licenses. In order to obtain a MATLAB licenses free solution, we use SWIFT, which is a fast and easy parallel scripting language. Once the original Mat file is complied to an executable file, SWIFT can run the code inside the executable file in parallel without using parallel MATLAB. We use the SWIFT to compute the object functions and gradients of different shots in parallel, and test the parallel 3D FWI code with overthrust data. This is joint work with Thomas Lai, Harsh Juneja, Bas Peters.},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/fang2013SINBADsfwi/fang2013SINBADsfwi.pdf},
	author = {Zhilong Fang and Henryk Modzelewski and Tristan van Leeuwen and Thomas Lai and Harsh Juneja and Bas Peters and Felix J. Herrmann}
}
