fracspy.detection.stacking.estimate_origin_times#

fracspy.detection.stacking.estimate_origin_times(idp, ds_full, tt, dt, x, y, z)[source]#

Estimate origin times of the events detected with diffraction stacking.

Estimate origin times of the events using traveltimes and the time indices determined from the maximum stack function (MSF) using STA/LTA analysis.

Parameters:

idpnumpy.ndarray: Integer array of time indices of determined peaks
ds_fullnumpy.ndarray: Diffraction stack full 4D array of shape (nt, nx, ny, nz)
ttnumpy.ndarray: Traveltime table of size $n_{r} \times n_{x} \times n_{y} \times n_{z}$
dtfloat: Time step
xnumpy.ndarray: Imaging area grid vector in X-axis
ynumpy.ndarray: Imaging area grid vector in Y-axis
znumpy.ndarray: Imaging area grid vector in Z-axis

Returns:

eotnumpy.ndarray: Array of estimated origin times for each of the determined peaks

Notes

For each determined MSF peak, the origin time $t_{e s t}$ of the corresponding event can be determined from peak time $t_{p e a k}$ as

t_{e s t} = t_{p e a k} - min_{R} T_{R} (r_{p e a k}) .

where $T_{R} (r_{p e a k})$ is the traveltime to receiver $R$ from $r_{p e a k}$ - the location determined from the maximum of the 4D imaging function $F (r, t)$ at the time moment $t_{p e a k}$ :

r_{p e a k} = \arg max_{r} F (r, t_{p e a k}) .