seisflows.plugins.preprocess.misfit

Misfit functions used by the ‘default’ preprocess class use to quantify differences between data and synthetics.

All functions defined have four required positional arguments

type syn:

np.array

param syn:

synthetic data array

type obs:

np.array

param obs:

observed data array

type nt:

int

param nt:

number of time steps in the data array

type dt:

float

param dt:

time step in sec

Functions

waveform(syn, obs, nt, dt, *args, **kwargs)

Direct waveform differencing

envelope(syn, obs, nt, dt, *args, **kwargs)

Waveform envelope difference from Yuan et al. 2015 Eq. 9

instantaneous_phase(syn, obs, nt, dt, *args, **kwargs)

Instantaneous phase difference from Bozdag et al. 2011

traveltime(syn, obs, nt, dt, *args, **kwargs)

Cross-correlation traveltime misfit function that is proportional to the

traveltime_inexact(syn, obs, nt, dt, *args, **kwargs)

A faster cc traveltime function but possibly innacurate

amplitude(syn, obs, nt, dt, *args, **kwargs)

Cross-correlation amplitude difference

envelope2(syn, obs, nt, dt, *args, **kwargs)

Envelope amplitude ratio from Yuan et al. 2015 Eq. B-1

envelope3(syn, obs, nt, dt[, eps])

Envelope cross-correlation lag from Yuan et al. 2015, Eq. B-4

instantaneous_phase2(syn, obs, nt, dt[, eps])

Alterative instantaneous phase function

displacement(*args, **kwargs)

velocity(*args, **kwargs)

acceleration(*args, **kwargs)

Module Contents

seisflows.plugins.preprocess.misfit.waveform(syn, obs, nt, dt, *args, **kwargs)

Direct waveform differencing

Parameters:

  • syn (np.array) – synthetic data array

  • obs (np.array) – observed data array

  • nt (int) – number of time steps in the data array

  • dt (float) – time step in sec

seisflows.plugins.preprocess.misfit.envelope(syn, obs, nt, dt, *args, **kwargs)

Waveform envelope difference from Yuan et al. 2015 Eq. 9

Parameters:

  • syn (np.array) – synthetic data array

  • obs (np.array) – observed data array

  • nt (int) – number of time steps in the data array

  • dt (float) – time step in sec

seisflows.plugins.preprocess.misfit.instantaneous_phase(syn, obs, nt, dt, *args, **kwargs)

Instantaneous phase difference from Bozdag et al. 2011

Parameters:

  • syn (np.array) – synthetic data array

  • obs (np.array) – observed data array

  • nt (int) – number of time steps in the data array

  • dt (float) – time step in sec

seisflows.plugins.preprocess.misfit.traveltime(syn, obs, nt, dt, *args, **kwargs)

Cross-correlation traveltime misfit function that is proportional to the squared time shift corresponding to the highest correlation coefficient of the cross-correlation between observed and synthetics

Parameters:

  • syn (np.array) – synthetic data array

  • obs (np.array) – observed data array

  • nt (int) – number of time steps in the data array

  • dt (float) – time step in sec

seisflows.plugins.preprocess.misfit.traveltime_inexact(syn, obs, nt, dt, *args, **kwargs)

A faster cc traveltime function but possibly innacurate

Parameters:

  • syn (np.array) – synthetic data array

  • obs (np.array) – observed data array

  • nt (int) – number of time steps in the data array

  • dt (float) – time step in sec

seisflows.plugins.preprocess.misfit.amplitude(syn, obs, nt, dt, *args, **kwargs)

Cross-correlation amplitude difference

Parameters:

  • syn (np.array) – synthetic data array

  • obs (np.array) – observed data array

  • nt (int) – number of time steps in the data array

  • dt (float) – time step in sec

seisflows.plugins.preprocess.misfit.envelope2(syn, obs, nt, dt, *args, **kwargs)

Envelope amplitude ratio from Yuan et al. 2015 Eq. B-1

Parameters:

  • syn (np.array) – synthetic data array

  • obs (np.array) – observed data array

  • nt (int) – number of time steps in the data array

  • dt (float) – time step in sec

seisflows.plugins.preprocess.misfit.envelope3(syn, obs, nt, dt, eps=0.0, *args, **kwargs)

Envelope cross-correlation lag from Yuan et al. 2015, Eq. B-4

Parameters:

  • syn (np.array) – synthetic data array

  • obs (np.array) – observed data array

  • nt (int) – number of time steps in the data array

  • dt (float) – time step in sec

seisflows.plugins.preprocess.misfit.instantaneous_phase2(syn, obs, nt, dt, eps=0.0, *args, **kwargs)

Alterative instantaneous phase function

Parameters:

  • syn (np.array) – synthetic data array

  • obs (np.array) – observed data array

  • nt (int) – number of time steps in the data array

  • dt (float) – time step in sec

seisflows.plugins.preprocess.misfit.displacement(*args, **kwargs)
seisflows.plugins.preprocess.misfit.velocity(*args, **kwargs)
seisflows.plugins.preprocess.misfit.acceleration(*args, **kwargs)