Exemplo n.º 1
0
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
Exemplo n.º 2
0
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