Elastic least-squares migration for quantitative reflection imaging of fracture compliances

Abstract

<jats:p> To quantitatively image fractures with high resolution, we develop an elastic least-squares migration (LSM) algorithm coupled with linear-slip theory, which accurately addresses seismic wave interaction with thin structures. We derive a linearized waveform inversion using the Born approximation to the boundary integral equation for scattered waves, including linear-slip interfaces for P-SV and SH wavefields. Numerical modeling tests assuming a laboratory-scale fracture where a 20 cm long fracture is illuminated by waves with 50 KHz center frequency, show that the proposed LSM successfully estimates the fracture compliances. Furthermore, due to the presence of coupling compliances at the fracture, the results using the proposed LSM show better images than those using the conventional LSM estimating Lamé constants. We also numerically illustrate that the proposed LSM can be successfully applied to dipole acoustic borehole logging data with 3 KHz center frequency for single-well reflection imaging of a 10 m long, dipping fracture embedded in a random background. Finally, we apply the LSM to laboratory experimental data, measuring PP reflections from a fluid-filled fracture. We confirm that the estimated fracture compliances correspond well to those estimated by earlier AVO inversion. Furthermore, the LSM resolves the spatially varying fracture compliances due to local filling of water in the fracture. Because the linear-slip theory can be applied to thin structures in a wide range of scales, high-resolution imaging results and estimated fracture compliance distributions will be crucial to further address small-scale properties at fractures, joints and geological faults. </jats:p>

Modelling 3D elastodynamic wave scattering due to density and Lamé parameter contrasts of near‐surface scatterers

“…(2017), Huang (2020) and Minato et al . (2020). To accelerate the calculations further, we make use of parallel programming techniques for the calculation of the wavefields (Harmankaya and Kaslilar, 2019).…”

“…While this method is relatively efficient when compared to more common modelling schemes, such as finite difference or finite element, depending on the size of the model and acquisition geometry, calculation of wavefields with our method can still be computationally expensive, even though the Born approximation already made the calculations significantly faster (Kaslilar, 2007) than the accurate solution of the wave equation (Riyanti and Herman, 2005). The use of wave number integral approach along with the Born approximation for modelling wave scattering due to inhomogeneities is previously examined by several studies, such as Lepage and Schmidt (2003), Malovichko et al (2017), Huang (2020) and Minato et al (2020). To accelerate the calculations further, we make use of parallel programming techniques for the calculation of the wavefields (Harmankaya and Kaslilar, 2019).…”