Course website: http://lagex.github.io/geofisica2
Note: This notebook is part of the course "Geofísica 2" of Geology program of the Universidade do Estado do Rio de Janeiro. All content can be freely used and adapted under the terms of the Creative Commons Attribution 4.0 International License.
Esse documento que você está usando é um IPython notebook. É um documento interativo que mistura texto (como esse), código (como abaixo), e o resultado de executar o código (que pode ser números, texto, figuras, videos, etc).
Vamos utilizar as simulações de ondas da biblioteca Fatiando a Terra. Essas simulações utilizam o método de diferenças finitas para calcular soluções da equação da onda.
As tarefas abaixo acompanham as questões que vocês devem responder. Usem as simulações abaixo para te ajudar a responder as perguntas.
Rode a célula abaixo para carregar as componentes necessárias para fazer as simulações. Não se preocupe se aparecer um :0: FutureWarning: IPython widgets are experimental and may change in the future.
abaixo. Isso é consequência de utilizar tecnologia de ponta.
%matplotlib inline
from __future__ import division
import numpy as np
from fatiando import utils
from fatiando.seismic.wavefd import Ricker, ElasticPSV, ElasticSH
from fatiando.vis import mpl
import fatiando
:0: FutureWarning: IPython widgets are experimental and may change in the future.
print(fatiando.__version__)
0.3-219-gb6b00ac
shape = (300, 1000)
spacing = 500
extent = [0, shape[1]*spacing, shape[0]*spacing, 0]
densidade = np.ones(shape)*2600
velocidade = np.ones(shape)*4000
l1 = 40e3//spacing
densidade[l1:,:] = 3500
velocidade[l1:,:] = 5000
Agora vamos criar o nosso simulador de ondas.
sim1 = ElasticSH(velocidade, densidade, spacing=spacing, taper=0.005)
e colocar uma fonte explosiva na superfície.
sim1.add_point_source((20e3//spacing, 50e3//spacing), Ricker(5, 1, 1))
Agora que temos nossa simulação pronta, rode a célcula abaixo para avançar a simulação 400 passos no tempo.
sim1.run(2000)
|##################################################|100% Ran 2000 iterations in 100.698 seconds.
sim1.animate(every=20, embed=True, dpi=50, cutoff=0.5)
sim1.explore(every=20)
<function fatiando.seismic.wavefd.plot>