def rho_gardner(dummy, keys, vals):
input = Struct(zip(keys, vals))
vp = input.vp
vs = input.vs
rho = 310. * vp**0.25
rho_water = 1050
thresh = 1.0e-3
rho[vs < thresh] = rho_water
return rho
def thomsen_voigt_2d(dummy, keys, vals):
input = Struct(zip(keys, vals))
output = Struct()
vp = input.vp
vs = input.vs
rho = input.rho
epslion = input.epsilon
delta = input.delta
output.c11 = rho * vp**2. * (1. + 2. * epsilon)
output.c12 = rho * vp**2. * (1. + 2. * epsilon) - 2. * rho * vs**2. * (
1. + 2. * epsilon)
output.c13 = rho * vp**2. * (1. + 2. * delta) - 2. * rho * vs**2.
output.c33 = rho * vp**2.
output.c55 = rho * vs**2.
return output
def lambda_mu_inverse(dummy, keys, vals):
input = Struct(zip(keys, vals))
output = Struct()
lame1 = input.lame1
lame2 = input.lame2
rho = input.rho
output.vp = ((lame1 + 2. * lame2) / rho)**0.5
output.vs = (lame2 / rho)**0.5
output.rho = rho
return output
def lambda_mu_forward(dummy, keys, vals):
input = Struct(zip(keys, vals))
output = Struct()
vp = input.vp
vs = input.vs
rho = input.rho
output.lame1 = rho * (vp**2. - 2. * vs**2.)
output.lame2 = rho * vs**2.
output.rho = rho
return output
def kappa_mu_inverse(dummy, keys, vals):
input = Struct(zip(keys, vals))
output = Struct()
kappa = input.kappa
mu = input.mu
rho = input.rho
output.vp = ((kappa + (4. / 3.) * mu) / rho)**0.5
output.vs = (mu / rho)**0.5
output.rho = rho
return output
def kappa_mu_forward(dummy, keys, vals):
input = Struct(zip(keys, vals))
output = Struct()
vp = input.vp
vs = input.vs
rho = input.rho
output.kappa = rho * (vp**2. - (4. / 3.) * vs**2.)
output.mu = rho * vs**2.
output.rho = rho
return output
def phi_beta_forward(dummy, keys, vals):
input = Struct(zip(keys, vals))
output = Struct()
vp = input.vp
vs = input.vs
rho = input.rho
kappa = rho * (vp**2. - (4. / 3.) * vs**2.)
output.bulk_c = (kappa / rho)**0.5
output.bulk_beta = vs
output.rho = rho
return output
def phi_beta_inverse(dummy, keys, vals):
input = Struct(zip(keys, vals))
output = Struct()
phi = input.bulk_c
vs = input.bulk_beta
rho = input.rho
kappa = rho * phi**2.
mu = rho * vs**2.
output.vp = ((kappa + (4. / 3.) * mu) / rho)**0.5
output.vs = vs
output.rho = rho
return output
def voigt_voigt_2d(dummy, keys, vals):
input = Struct(zip(keys, vals))
output = input
return output
def tti_voight_2d(dummy, keys, vals):
input = Struct(zip(keys, vals))
output = Struct()
vp = input.vp
vs = input.vs
rho = input.rho
epsilon = input.epsilon
delta = input.delta
theta = input.theta
c11 = rho * vp**2. * (1. + 2. * epsilon)
c12 = rho * vp**2. * (1. + 2. * epsilon) - 2. * rho * vs**2. * (
1. + 2. * epsilon)
c13 = rho * vp**2. * (1. + 2. * delta) - 2. * rho * vs**2.
c33 = rho * vp**2.
c55 = rho * vs**2.
sint = sin(PI / 180. * theta)
cost = cos(PI / 180. * theta)
sin2t = sin(2. * PI / 180. * theta)
cos2t = cos(2. * PI / 180. * theta)
output.c11 = c11 * cost**4 + c33 * sint**4 + 2 * c13 * sint**2 * cost**2 + c55 * sin2t**2
output.c12 = 1.e-6
output.c13 = c13 * (cost**4 + sint**4) + (
c11 + c33) * sint**2 * cost**2 - c55 * sin2t**2
output.c15 = ((c11 - c13) * cost**2 +
(c13 - c33) * sint**2) * sint * cost - c55 * sin2t * cos2t
output.c23 = 1.e-6
output.c25 = 0.0
output.c33 = c11 * sint**4 + c33 * cost**4 + 2 * c13 * sint**2 * cost**2 + c55 * sin2t**2
output.c35 = ((c11 - c13) * sint**2 +
(c13 - c33) * cost**2) * sint * cost + c55 * sin2t * cos2t
output.c55 = (c11 - 2. * c13 + c33) * sint**2 * cost**2 + c55 * cos2t**2
return output
def vp_vs_inverse(dummy, keys, vals):
input = Struct(zip(keys, vals))
output = input
return output
def vp_vs_forward(dummy, keys, vals):
input = Struct(zip(keys, vals))
output = input
return output