Module JHVIT.JHVI_Tetra

JHVI_Tetra.jntHypoVel_T(data, caldata, Vinit, cells, nodes, rcv, Hypo0, par, threads=1, vPoints=array([], dtype=float64), basename='Vel')

Joint hypocenter-velicoty inversion from P wave arrival time data parametrized using the velocity model.

Parameters

  • data (np.ndarray, shape(arrival time number, 3)) – Arrival times and corresponding receivers for each event..

  • caldata (np.ndarray, shape(number of calibration shots, 3)) – Calibration shot data.

  • Vinit (np.ndarray, shape(nnodes,1) or (1,1)) – Initial velocity model.

  • cells (np.ndarray of int, shape (cell number, 4)) – Indices of nodes forming the cells.

  • nodes (np.ndarray, shape (nnodes, 3)) – Node coordinates.

  • rcv (np.ndarray, shape (receiver number,3)) – Coordinates of receivers.

  • Hypo0 (np.ndarray, shape(event number, 5)) – First guesses of the hypocenter coordinates (must be all diffirent).

  • par (instance of the class Parameters) – The inversion parameters.

  • threads (int, optional) – Thread number. The default is 1.

  • vPoints (np.ndarray, shape(point number,4), optional) – Known velocity points. The default is np.array([]).

  • basename (string, optional) – The filename used to save the output file. The default is ‘Vel’.

  • Returns –

  • ------- –

  • output (python dictionary) – It contains the estimated hypocenter coordinates and their origin times, static correction values, velocity model, convergence states, parameter uncertainty and residual norm in each iteration.

JHVI_Tetra.jntHyposlow_T(data, caldata, Vinit, cells, nodes, rcv, Hypo0, par, threads=1, vPoints=array([], dtype=float64), basename='Slowness')

Joint hypocenter-velicoty inversion from P wave arrival time data parametrized using the slowness model.

Parameters

  • data (np.ndarray, shape(arrival time number, 3)) – Arrival times and corresponding receivers for each event.

  • caldata (np.ndarray, shape(number of calibration shots, 6)) – Calibration shot data.

  • Vinit (np.ndarray, shape(nnodes,1) or (1,1)) – Initial velocity model.

  • Cells (np.ndarray of int, shape (cell number, 4)) – Indices of nodes forming the cells.

  • nodes (np.ndarray, shape (nnodes, 3)) – Node coordinates.

  • rcv (np.ndarray, shape (receiver number,3)) – Coordinates of receivers.

  • Hypo0 (np.ndarray, shape(event number, 5)) – First guesses of the hypocenter coordinates (must be all diffirent).

  • par (instance of the class Parameters) – The inversion parameters.

  • threads (int, optional) – Thread number. The default is 1.

  • vPoints (np.ndarray, shape(point number,4), optional) – Known velocity points. The default is np.array([]).

  • basename (string, optional) – The filename used to save the output files. The default is ‘Slowness’.

Returns

output – It contains the estimated hypocenter coordinates and their origin times, static correction values, velocity model, convergence states, parameter uncertainty and residual norm in each iteration.

Return type

python dictionary

JHVI_Tetra.jntHypoVelPS_T(obsData, calibdata, Vinit, cells, nodes, rcv, Hypo0, par, threads=1, vPnts=(array([], dtype=float64), array([], dtype=float64)), basename='Vel')

Joint hypocenter-velocity inversion from P- and S-wave arrival time data parametrized using the velocity models.

Parameters

  • obsData (tuple of two np.ndarrays (shape(observed data number, 3))) – Observed arrival time data of P- and S-waves.

  • calibdata (tuple of two np.ndarrays (shape (number of calibration shots, 5))) – Calibration data of P- and S-waves.

  • Vinit (tuple of np.ndarrays (shape (nnodes, 1) or (1,1))) – Initial velocity models of P- and S-waves.

  • cells (np.ndarray of int, shape (cell number, 4)) – Indices of nodes forming the cells.

  • nodes (np.ndarray, shape (nnodes, 3)) – Node coordinates.

  • rcv (np.ndarray, shape (receiver number,3)) – Coordinates of receivers.

  • Hypo0 (np.ndarray, shape(event number, 5)) – First guesses of the hypocenter coordinates (must be all diffirent).

  • par (instance of the class Parameters) – The inversion parameters.

  • threads (int, optional) – Thread number. The default is 1.

  • vPnts (tuple of two np.ndarrays, optional) – Known velocity points of P- and S-waves. The default is (np.array([]), np.array([])).

  • basename (string, optional) – The filename used to save the output files. The default is ‘Vel’.

Raises

ValueError – If the Vs/Vp ratio is inverted instead of Vs model and some known velocity points are given for the S wave and not for the P wave.

Returns

output – It contains the estimated hypocenter coordinates and their origin times, static correction values, velocity models of P- and S-waves, hypocenter convergence states, parameter uncertainty and residual norm in each iteration.

Return type

python dictionary

JHVI_Tetra.jntHyposlowPS_T(obsData, calibdata, Vinit, cells, nodes, rcv, Hypo0, par, threads=1, vPnts=(array([], dtype=float64), array([], dtype=float64)), basename='Slowness')

Joint hypocenter-velocity inversion from P- and S-wave arrival time data parametrized using the slowness models.

Parameters

  • obsData (tuple of two np.ndarrays (shape(observed data number, 3))) – Observed arrival time data of P- and S-waves.

  • calibdata (tuple of two np.ndarrays (shape (calibration shot number, 6))) – Calibration data of P- and S-waves.

  • Vinit (tuple of np.ndarrays (shape (nnodes, 1) or (1,1))) – Initial velocity models of P- and S-waves.

  • cells (np.ndarray of int, shape (cell number, 4)) – Indices of nodes forming the cells.

  • nodes (np.ndarray, shape (nnodes, 3)) – Node coordinates.

  • rcv (np.ndarray, shape (receiver number,3)) – Coordinates of receivers.

  • Hypo0 (np.ndarray, shape(event number, 5)) – First guesses of the hypocenter coordinates (must be all diffirent).

  • par (instance of the class Parameters) – The inversion parameters.

  • threads (int, optional) – Thread number. The default is 1.

  • vPnts (tuple of two np.ndarrays, optional) – Known velocity points of P- and S-waves. The default is (np.array([]),np.array([])).

  • basename (string, optional) – The filename used to save the output files. The default is ‘Slowness’.

Raises

ValueError – If the Ss/Sp ratio is inverted instead of Vs model and some known velocity points are given for the S wave and not for the P wave.

Returns

output – It contains the estimated hypocenter coordinates and their origin times, static correction values, velocity models of P- and S-waves, hypocenter convergence states, parameter uncertainty and residual norm in each iteration.

Return type

python dictionary

JHVI_Tetra.jointHypoVel_T(inputFileParam, model='slow')

Joint hypocenter-velocity inversion using P wave data.

Parameters

  • inputFileParam (string) – Text file containing inversion parameters and data filenames.

  • model (string) – Sought model : ‘vel’ for an inversion problem parametrized using the velocity model,’slow’ for an inversion problem parametrized using the slowness model. The default is ‘slow’.

Returns

It contains the estimated hypocenter coordinates and their origin times, static correction values, velocity model, convergence states, parameter uncertainty and residual norm in each iteration.

Return type

python dictionary

JHVI_Tetra.jointHypoVelPS_T(inputFileParam, model='slow')

Joint hypocenter-velocity inversion using P- and S-wave arrival time data.

Parameters

  • inputFileParam (string) – Text file containing inversion parameters and data filenames.

  • model (string) – Sought models: ‘vel’ for an inversion problem parametrized using the velocity model, ‘slow’ for an inversion problem parametrized using the slowness model. The default is ‘slow’.

Returns

It contains the estimated hypocenter coordinates and their origin times, static correction values, velocity models of P and S waves, hypocenter convergence states, parameter uncertainty and residual norm in each iteration.

Return type

python dictionary

JHVI_Tetra.readEventsFiles(time_file, waveType=False)

Read a list of seismic events and corresponding data from a text file.

Parameters

  • time_file (string) – Event data filename.

  • waveType (bool) – True if the seismic phase of each event is identified. The default is False.

Returns

data – Event arrival time data

Return type

np.ndarray or a list of two np.ndarrays