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

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 (_κs,_κn) 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 (_κs,_κn) characterizing the trial (infinitesimal) fracture, the ratio between the shear and normal specific stiffness along the surface of a nearly-planar (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 canonical numerical experiments.