On the elastic-wave imaging and characterization of fractures with specific stiffness

TL;DR

This paper extends topological sensitivity methods to 3D elastic-wave imaging, enabling simultaneous fracture reconstruction and interface characterization, including specific stiffness, with no iterations needed for geometry and qualitative interfacial condition assessment.

Contribution

It introduces a novel TS-based imaging framework that reconstructs fracture geometry and characterizes interface stiffness without iterative procedures, using elastic waves and asymptotic analysis.

Findings

TS accurately reconstructs fracture geometry regardless of boundary conditions.

The method qualitatively identifies shear and normal stiffness ratios.

Numerical experiments validate the approach's effectiveness.

Abstract

The concept of topological sensitivity (TS) is extended to enable simultaneous 3D reconstruction of fractures with unknown boundary condition and characterization of their interface by way of elastic waves. Interactions between the two surfaces of a fracture, due to e.g. presence of asperities, fluid, or proppant, are described via the Schoenberg's linear slip model. The proposed TS sensing platform is formulated in the frequency domain, and entails point-wise interrogation of the subsurface volume by infinitesimal fissures endowed with interfacial stiffness. For completeness, the featured elastic polarization tensor - central to the TS formula - is mathematically described in terms of the shear and normal specific stiffness (ks,kn) of a vanishing fracture. Simulations demonstrate that, irrespective of the contact condition between the faces of a hidden fracture, the TS (used as a waveform imaging tool) is capable of reconstructing its geometry and identifying the normal vector to the fracture surface without iterations. On the basis of such geometrical information, it is further shown via asymptotic analysis -- assuming "low frequency" elastic-wave illumination, that by certain choices of (ks,kn) characterizing the trial (infinitesimal) fracture, the ratio between the shear and normal specific stiffness along the surface of a nearly-panar (finite) fracture can be qualitatively identified. This, in turn, provides a valuable insight into the interfacial condition of a fracture at virtually no surcharge -- beyond the computational effort required for its imaging. The proposed developments are integrated into a computational platform based on a regularized boundary integral equation (BIE) method for 3D elastodynamics, and illustrated via a set of numerical experiments.