Numerical Analysis and Scientific Computing Seminar

Untangling the nonlinearity in inverse scattering using data-driven reduced order models

Speaker: Liliana Borcea, University of Michigan

Location: Warren Weaver Hall 1302

Date: Oct. 6, 2017, 10 a.m.


We discuss an inverse problem for the wave equation, where an array of sensors probes an unknown, heterogeneous medium with pulses and measures the scattered waves. The goal in inversion is to determine from these measurements scattering structures in the medium, modeled mathematically by a reflectivity function. Most imaging methods assume a linear mapping between the unknown reflectivity and the array data. The linearization, known as the Born (single scattering) approximation is not accurate in strongly scattering media, so the reconstruction of the reflectivity may be poor. We show that  it is possible to remove the multiple scattering (nonlinear) effects from the data using a reduced order model (ROM). The ROM is defined by an orthogonal projection of the wave propagator operator on the subspace spanned by the time snapshots of the solution of the wave equation. The snapshots are known only at the sensor locations, which is enough information to construct the ROM. The main result discussed in the talk is a novel, linear-algebraic algorithm that uses the ROM to map the data to its Born approximation.