Proper Generalized Decomposition solution of the parameterized Helmholtz problem: application to inverse geophysical problems.

62/2015

Proper Generalized Decomposition solution of the parameterized Helmholtz problem: application to inverse geophysical problems.

Thursday 10th December 2015

Signorini, M.; Zlotnik, S.; Díez, P.

The identification of the geological structure from seismic data is formulated as an inverse problem. The properties and the shape of the rock formations in the subsoil are described by material and geometric parameters, which are taken as input data for a predictive model. Here, the model is based on the Helmholtz equation, describing the acoustic response of the system for a given wave length. Thus, the inverse problem consists in identifying the values of these parameters such that the output of the model agrees the best with observations. This optimization algorithm requires multiple queries to the model with different values of the parameters. Reduced Order Models are especially well suited to significantly reduce the computational overhead of the multiple evaluations of the model. In particular, the Proper Generalized Decomposition (PGD) produces a solution explicitly stating the parametric dependence, where the parameters play the same role as the physical coordinates. A PGD solver is devised to inexpensively explore the parametric space along the iterative process. This exploration of the parametric space is in fact seen as a post-process of the generalized solution. The approach adopted demonstrates its viability when tested in two illustrative examples.

This report, or a modified version of it, has been also submitted to, or published on
International Journal for Numerical Methods in Engineering