Elastic waveform inversion for double-couple microseismic source estimation in vertically fractured transversely isotropic media

TL;DR

This paper introduces a waveform inversion method tailored for double-couple microseismic sources in vertically fractured transversely isotropic media, accounting for anisotropy and providing more accurate source characterization.

Contribution

The study develops a novel waveform inversion algorithm that incorporates VFTI media properties, enabling better estimation of source parameters without assuming fault orientation.

Findings

Successfully inverts synthetic data for source location and slip parameters.

Demonstrates improved accuracy over isotropic models in fractured media.

Provides a framework extendable to more complex anisotropic formations.

Abstract

Accurate characterization of microseismic events during fluid injection in sedimentary formations is essential to mitigate environmental risks. The source mechanism for microseismic events related to a slip on a fault plane is given by a double-couple. Waveform inversion has emerged as a promising technique for estimating the moment tensor and the position vector of double-couple sources. In most applications of waveform inversion for the moment tensor of double-couple sources, the formation is typically assumed to be isotropic or, less frequently, transversely isotropic. Modification of the moment-tensor representation to account for anisotropy created by aligned vertical fractures in transversely isotropic formations has not been included while inverting microseismic waveform data. In this study on synthetic microseismic data, we present a waveform inversion algorithm that includes this modification, considering the formation in the focal region to be vertically fractured transversely isotropic (VFTI) and possessing orthorhombic symmetry. Since VFTI media lack rotational symmetry, no assumptions have been made about the orientation of the fault plane where the slip occurred. The moment tensor of double-couple sources is formulated in terms of the elastic parameters of the VFTI medium and geometrical parameters which are slip magnitude, slip angle, fault dip, and azimuth angle of the fault-normal. Source inversion is treated as a local optimization problem, and we invert for the source location and the geometrical parameters. These geometrical parameters are more directly constrained by seismic data than the moment tensor components and offer geologically meaningful insights. This approach enhances microseismic monitoring in fractured formations and can be extended to more complex anisotropic media, such as monoclinic systems.