bruges.rockphysics.elastic module

Elastic impedance.

copyright:2015 Agile Geoscience
license:Apache 2.0
bruges.rockphysics.elastic.elastic_impedance(vp, vs, rho, theta1, k=None, normalize=False, constants=None, use_sin=False, rho_term=False)[source]

Returns the elastic impedance as defined by Connolly, 1999; we are using the formulation reported in Whitcombe et al. (2001). Inputs should be shape m x 1, angles should be n x 1. The result will be m x n.

Parameters:
  • vp (ndarray) – The P-wave velocity scalar or 1D array.
  • vs (ndarray) – The S-wave velocity scalar or 1D array.
  • rho (ndarray) – The bulk density scalar or 1D array.
  • theta1 (ndarray) – Incident angle(s) in degrees, scalar or array.
  • k (float) – A constant, see Connolly (1999). Default is None, which computes it from Vp and Vs.
  • normalize (bool) – if True, returns the normalized form of Whitcombe.
  • constants (tuple) – A sequence of 3 constants to use for normalization. If you don’t provide this, then normalization just uses the means of the inputs. If these are scalars, the result will be the acoustic impedance (see Whitcombe et al., 2001).
  • use_sin (bool) – If True, use sin^2 for the first term, instead of tan^2 (see Connolly).
  • rho_term (bool) – If True, provide alternative form, with Vp factored out; use in place of density in generating synthetics in other software (see Connolly). In other words, the result can be multipled with Vp to get the elastic impedance.
Returns:

The elastic impedance log at the specficied angle or angles.

Return type:

ndarray