def blangy(vp1, vs1, rho1, vp2, vs2, rho2, d1=0, e1=0, d2=0, e2=0, theta1=0):
"""Implements the Blangy equation with the same interface as the other
reflectivity equations. Wraps bruges.anisotropy.blangy(), which you may
prefer to use directly.
Args:
vp1 (ndarray): The upper P-wave velocity; float or 1D array length m.
vs1 (ndarray): The upper S-wave velocity; float or 1D array length m.
rho1 (ndarray): The upper layer's density; float or 1D array length m.
vp2 (ndarray): The lower P-wave velocity; float or 1D array length m.
vs2 (ndarray): The lower S-wave velocity; float or 1D array length m.
rho2 (ndarray): The lower layer's density; float or 1D array length m.
d1 (ndarray): The upper delta; float or 1D array length m.
e1 (ndarray): The upper epsilon; float or 1D array length m.
d2 (ndarray): The lower delta; float or 1D array length m.
e2 (ndarray): The lower epsilon; float or 1D array length m.
theta1 (ndarray): The incidence angle; float or 1D array length n.
Returns:
ndarray. The Blangy approximation for P-P reflectivity at the
interface. Wraps `anisotropy.blangy()`.
"""
_, anisotropic = anisotropy.blangy(
vp1,
vs1,
rho1,
d1,
e1, # UPPER
vp2,
vs2,
rho2,
d2,
e2, # LOWER
theta1)
return anisotropic
def blangy(vp1, vs1, rho1, vp2, vs2, rho2,
d1=0, e1=0, d2=0, e2=0,
theta1=0, terms=False):
"""Implements the Blangy equation with the same interface as the other
reflectivity equations. Wraps bruges.anisotropy.blangy().
Note that the anisotropic parameters come after the other rock properties,
and all default to zero.
:param vp1: The p-wave velocity of the upper medium.
:param vs1: The s-wave velocity of the upper medium.
:param rho1: The density of the upper medium.
:param vp2: The p-wave velocity of the lower medium.
:param vs2: The s-wave velocity of the lower medium.
:param rho2: The density of the lower medium.
:param d1: Thomsen's delta for the upper medium.
:param e1: Thomsen's epsilon for the upper medium.
:param d2: Thomsen's delta for the upper medium.
:param e2: Thomsen's epsilon for the upper medium.
:param theta1: An array of incident angles to use for reflectivity
calculation [degrees].
:returns: a vector of len(theta1) containing the reflectivity
value corresponding to each angle.
:param theta1: An array of incident angles to use for reflectivity
calculation [degrees].
:returns: a vector of len(theta1) containing the reflectivity
value corresponding to each angle.
"""
_, anisotropic = anisotropy.blangy(vp1, vs1, rho1, d1, e1, # UPPER
vp2, vs2, rho2, d2, e2, # LOWER
theta1)
return anisotropic
def blangy(vp1, vs1, rho1, vp2, vs2, rho2,
d1=0, e1=0, d2=0, e2=0,
theta1=0):
"""Implements the Blangy equation with the same interface as the other
reflectivity equations. Wraps bruges.anisotropy.blangy(), which you may
prefer to use directly.
Args:
vp1 (ndarray): The upper P-wave velocity; float or 1D array length m.
vs1 (ndarray): The upper S-wave velocity; float or 1D array length m.
rho1 (ndarray): The upper layer's density; float or 1D array length m.
vp2 (ndarray): The lower P-wave velocity; float or 1D array length m.
vs2 (ndarray): The lower S-wave velocity; float or 1D array length m.
rho2 (ndarray): The lower layer's density; float or 1D array length m.
d1 (ndarray): The upper delta; float or 1D array length m.
e1 (ndarray): The upper epsilon; float or 1D array length m.
d2 (ndarray): The lower delta; float or 1D array length m.
e2 (ndarray): The lower epsilon; float or 1D array length m.
theta1 (ndarray): The incidence angle; float or 1D array length n.
Returns:
ndarray. The Blangy approximation for P-P reflectivity at the
interface. Wraps `anisotropy.blangy()`.
"""
_, anisotropic = anisotropy.blangy(vp1, vs1, rho1, d1, e1, # UPPER
vp2, vs2, rho2, d2, e2, # LOWER
theta1)
return anisotropic
def blangy(vp1,
vs1,
rho1,
vp2,
vs2,
rho2,
d1=0,
e1=0,
d2=0,
e2=0,
theta1=0,
terms=False):
"""Implements the Blangy equation with the same interface as the other
reflectivity equations. Wraps bruges.anisotropy.blangy().
Note that the anisotropic parameters come after the other rock properties,
and all default to zero.
:param vp1: The p-wave velocity of the upper medium.
:param vs1: The s-wave velocity of the upper medium.
:param rho1: The density of the upper medium.
:param vp2: The p-wave velocity of the lower medium.
:param vs2: The s-wave velocity of the lower medium.
:param rho2: The density of the lower medium.
:param d1: Thomsen's delta for the upper medium.
:param e1: Thomsen's epsilon for the upper medium.
:param d2: Thomsen's delta for the upper medium.
:param e2: Thomsen's epsilon for the upper medium.
:param theta1: An array of incident angles to use for reflectivity
calculation [degrees].
:returns: a vector of len(theta1) containing the reflectivity
value corresponding to each angle.
:param theta1: An array of incident angles to use for reflectivity
calculation [degrees].
:returns: a vector of len(theta1) containing the reflectivity
value corresponding to each angle.
"""
_, anisotropic = anisotropy.blangy(
vp1,
vs1,
rho1,
d1,
e1, # UPPER
vp2,
vs2,
rho2,
d2,
e2, # LOWER
theta1)
return anisotropic