Browsing by Author "Bencomo, Mario J."

Now showing 1 - 2 of 2
  • Thumbnail Image
    Item
    Discretization of Multipole Sources in a Finite Difference Setting for Wave Propagation Problems
    (2018-06-20) Bencomo, Mario J.; Symes, William
    Seismic sources are commonly idealized as point-sources due to their small spatial extent relative to seismic wavelengths. The acoustic isotropic point-radiator is inadequate as a model of seismic wave generation for seismic sources that are known to exhibit directivity. Therefore, accurate modeling of seismic wavefields must include source representations generating anisotropic radiation patterns. Such seismic sources can be modeled as linear combinations of multipole point-sources. In this paper we present a method for discretizing multipole sources in a finite difference setting, an extension of the moment matching conditions developed for the Dirac delta function in other applications. We also provide the necessary analysis and numerical evidence to demonstrate the accuracy of our singular source approximations. In particular, we develop a weak convergence theory for the discretization of a family of symmetric hyperbolic systems of first-order partial differential equations, with singular source terms, solved via staggered-grid finite difference methods. Numerical experiments demonstrate a stronger result than what is presented in our convergence theory, namely, optimal convergence rates of numerical solutions are achieved point-wise in space away from the source if an appropriate source discretization is used.
  • Thumbnail Image
    Item
    Representation and Estimation of Seismic Sources via Multipoles
    (2017-05) Bencomo, Mario J.
    Accurate representation and estimation of seismic sources are essential to the seismic inversion problem. General sources can be approximated by a truncated series of multipoles depending on the source anisotropy. Most research in the joint inversion of source and medium parameters assumes seismic sources can be modeled as isotropic point-sources resulting in an inability to fit the anisotropy observed in data, ultimately impacting the recovery of medium parameters. In this thesis I lay the groundwork for joint source-medium parameter inversion with potentially anisotropic seismic sources via full waveform inversion through three key contributions: a mathematical and computational framework for the modeling and inversion of sources via multipoles, construction and analysis of discretizations of multipole sources on regular grids, and preconditioners based on fractional time derivative/integral operators for the ill-conditioned source estimation subproblem. As an application of my multipole framework, I also study the efficacy of multipoles in modeling the airgun array source, the most common type of active source in marine seismic surveying. Inversion results recovered a dominating isotropic component of the multipole source model that accounted for 84% of the observed radiation pattern. An extra 10% of the observed output pressure field can be explained when incorporating dipole terms in the source representation, thus motivating the use of multipoles to capture source anisotropy.