  In micro-seismic event measurements, pinpointing the passive source's exact
spatial and temporal location is paramount. This research advocates for the
combined use of both P- and S-wave data, captured by geophone monitoring
systems, to improve source inversion accuracy. Drawing inspiration from the
secondary source concept in Elastic Reflection Full Waveform Inversion (ERFWI),
we introduce an equivalent source term. This term combines source functions and
source images. Our optimization strategy iteratively refines the spatial
locations of the source, its temporal functions, and associated velocities
using a full waveform inversion framework. Under the premise of an isotropic
medium with consistent density, the source is defined by two spatial and three
temporal components. This offers a nuanced source representation in contrast to
the conventional seismic moment tensor. To address gradient computation, we
employ the adjoint-state method. However, we encountered pronounced
non-linearity in waveform inversion of micro-seismic events, primarily due to
the unknown source origin time, resulting in cycle skipping challenges. To
counteract this, we devised an objective function that is decoupled from the
source origin time. This function is formulated by convolving reference traces
with both observed and predicted data. Through the concurrent inversion of the
source image, source time function, and velocity model, our method offers
precise estimations of these parameters, as validated by a synthetic 2D example
based on a modified Marmousi model. This nested inversion approach promises
enhanced accuracy in determining the source image, time function, and velocity
model.
