rockphypy.utils

Module Contents

Classes

utils

Basic calculations for velocities, moduli and stiffness matrix.

class rockphypy.utils.utils

Basic calculations for velocities, moduli and stiffness matrix.

static V(K, G, rho)

Compute velocity given density and elastic moduli.

Parameters:

• K (float or array) – (GPa): bulk modulus

• G (float or array) – (GPa): shear moulus

• rho (float or array) – (g/m3): density of the frame

Returns:

float or array – Vp, Vs (m/s): velocity

static poi(K, G)

Compute poisson’s ratio from K an G

Parameters:

• K (float or array) – (GPa): bulk modulus

• G (float or array) – (GPa): shear moulus

Returns:

float or array – Poisson’s ratio

static lame(K, G)

Compute lame constant lamdba from K an G

Parameters:

• K (float or array) – (GPa): bulk modulus

• G (float or array) – (GPa): shear moulus

Returns:

float or array – Poisson’s ratio

static M_from_V(den, vp, vs)

_summary_

Parameters:

• den (float or array) – (g/cm3): bulk density

• vp (float or array) – (m/s): p wave velocity

• vs (float or array) – (m/s): s wave velocity

Returns:

float or array – K, G (GPa):bulk and shear moduli

static write_HTI_matrix(C11, C33, C13, C44, C55)

formulate HTI stiffness matrix

Parameters:

• C11 (float) – (GPa): stiffness

• C33 (float) – (GPa): stiffness

• C13 (float) – (GPa): stiffness

• C44 (float) – (GPa): stiffness

• C55 (float) – (GPa): stiffness

Returns:

2d array – C: 6x6 stiffness matrix

static write_VTI_compliance(S11, S12, S13, S33, S44)

formulate VTI compliance matrix

Parameters:

• S11 (float) – (GPa): stiffness

• S12 (float) – (GPa): stiffness

• S13 (float) – (GPa): stiffness

• S33 (float) – (GPa): stiffness

• S44 (float) – (GPa): stiffness

Returns:

2d array – S: 6x6 compliance matrix

static write_VTI_matrix(C11, C33, C13, C44, C66)

formulate VTI stiffness matrix

Parameters:

• C11 (float) – (GPa): stiffness

• C33 (float) – (GPa): stiffness

• C13 (float) – (GPa): stiffness

• C44 (float) – (GPa): stiffness

• C66 (float) – (GPa): stiffness

Returns:

2d array – C: 6x6 stiffness matrix

static write_matrix(C11, C22, C33, C12, C13, C23, C44, C55, C66)

formulate general 6x6 stiffness matrix in Voigt notation

Parameters:

• C11 (float) – (GPa): stiffness

• C22 (float) – (GPa): stiffness

• C33 (float) – (GPa): stiffness

• C12 (float) – (GPa): stiffness

• C13 (float) – (GPa): stiffness

• C23 (float) – (GPa): stiffness

• C44 (float) – (GPa): stiffness

• C55 (float) – (GPa): stiffness

• C66 (float) – (GPa): stiffness

Returns:

2d array – C: 6x6 stiffness matrix

static write_iso(K, G)

formulate isotropic 6x6 stiffness matrix in Voigt notation

Parameters:

K –

Gfloat or array

(GPa): shear moulus

Returns:

2d array – C: 6x6 stiffness matrix

static crack_por(crd, alpha)

compute crack porosity from crack aspect ratio and crack density

Parameters:

• crd (float or array) – (unitless): crack density

• alpha (float or array) – crack aspect ratio

Returns:

float or array – cpor (frac): crack porosity

static v_to_c_VTI(Vp0, Vp45, Vp90, Vs0, Vsh90, den)

compute stiffness matrix given velocity measurements along different directions

Parameters:

• Vp0 (float or array) – (km/s): incident angle dependent velocity measurements

• Vp45 (float or array) – (km/s): incident angle dependent velocity measurements

• Vp90 (float or array) – (km/s): incident angle dependent velocity measurements

• Vs0 (float or array) – (km/s): incident angle dependent velocity measurements

• Vsh90 (float or array) – (km/s): incident angle dependent velocity measurements

• den (float or array) – (g/cm3):density of the sample

Returns:

2d array – C: VTI stiffness matrix