24 research outputs found

    On the stress-wave imaging of cavities in a semi-infinite solid

    No full text
    The problem of mapping underground cavities from surface seismic measurements is investigated within the framework of a regularized boundary integral equation (BIE) method. With the ground modeled as a uniform elastic half-space, the inverse analysis of elastic waves scattered by a three-dimensional void is formulated as a task of minimizing the misfit between experimental observations and theoretical predictions for an assumed void geometry. For an accurate treatment of the gradient search technique employed to solve the inverse problem, sensitivities of the predictive BIE model with respect to cavity parameters are evaluated semi-analytically using an adjoint problem approach and a continuum kinematics description. Several key features of the formulation, including the rigorous treatment of the radiation condition for semi-infinite solids, modeling of an illuminating seismic wave field, and treatment of the prior information, are highlighted. A set of numerical examples with spherical and ellipsoidal cavity geometries is included to illustrate the performance of the method. It is shown that the featured adjoint problem approach reduces the computational requirements by an order of magnitude relative to conventional finite-difference estimates, thus rendering the three-dimensional elastic-wave imaging of solids tractable for engineering applications

    A computational basis for elastodynamic cavity identification in a semi-infinite solid

    No full text
    The focus of this paper is a computational platform for the non-intrusive, active seismic imaging of subterranean openings by means of an elastodynamic boundary integral equation (BIE) method. On simulating the ground response to steady-state seismic excitation as that of a uniform, semi-infinite elastic solid, solution to the 3D inverse scattering problem is contrived as a task of minimizing the misfit between experimental observations and BIE predictions of the surface ground motion. The forward elastodynamic solution revolves around the use of the half-space Greenrsquos functions, which analytically incorporate the traction-free boundary condition at the ground surface and thus allow the discretization and imaging effort to be focused on the surface of a hidden cavity. For a rigorous approach to the gradient-based minimization employed to resolve the cavity, sensitivities of the trial boundary element model with respect to (geometric) void parameters are evaluated using an adjoint field approach. Details of the computational treatment, including the regularized (i.e. Cauchy principal value-free) boundary integral equations for the primary and adjoint problem, the necessary evaluation of surface displacement gradients and their implementation into a parallel code, are highlighted. Through a suite of numerical examples involving the identification of an ellipsoidal cavity, a parametric study is presented which illustrates the importance of several key parameters on the imaging procedure including the prior information, ldquomeasurementrdquo noise, and the amount of experimental input
    corecore