Modal Mineral Composition of a Rock and calculate seismic Anisotropy

Let’s define the modal mineral composition for a rock by assuming 60% Forsterite, 25% Diopside and 15% Enstatite and calculate anisotropy for that modal rock. For that, first import the Material module

from santex.material import Material

material = Material()

Let’s define our model rock with 20% Forsterite, 50% Diopside and 30% Enstatite

rock = ["Forsterite", "Diopside", "Enstatite"]

fraction = [0.6, 0.25, 0.15]

average_tensor, average_density = material.modalRock(rock, fraction, 2, 1000)

Let’s look at the average tensor and the average density

average_tensor

array([[281.828 ,  75.0585,  74.522 ,   0.    ,   2.095 ,   0.    ],
       [ 75.0585, 185.4755,  74.368 ,   0.    ,   2.445 ,   0.    ],
       [ 74.522 ,  74.368 , 226.12  ,   0.    ,  13.1425,   0.    ],
       [  0.    ,   0.    ,   0.    ,  65.984 ,   0.    ,   2.75  ],
       [  2.095 ,   2.445 ,  13.1425,   0.    ,  71.4465,   0.    ],
       [  0.    ,   0.    ,   0.    ,   2.75  ,   0.    ,  74.602 ]])

average_density

3177.1809637264096

Lets import anisotropy class and check anisotropy values and plot it

from santex.anisotropy import Anisotropy

anisotropy = Anisotropy(average_tensor*10**9, average_density)

values = anisotropy.anisotropy_values(return_values = "maxvp")

values

9420.564966931908

values = anisotropy.anisotropy_values()

Max Vp:  9420.564966931908
Min Vp:  7640.509136106333
Max Vs1:  5337.243472897191
Min Vs1:  4664.683873568134
Max Vs2:  4871.668248288787
Min Vs2:  4412.924524052354
Max vs anisotropy percent:  15.906369992412309
Min vs anisotropy percent:  0.10536992746704398
P wave anisotropy percent:  20.866867116046134
S1 Wave anisotropy percent:  13.448599975418523
S2 Wave anisotropy percent:  9.881827571437134
Velocity difference, maxdvs:  779.25752900316
Vp/Vs1 ratio:  24.449690400322446
Mean vp:  8423.657134402454
Mean Vs1:  4933.5859970422935
Mean Vs2:  4588.761738655829

anisotropy.plot()