{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# 4 - Sísmica de refração\n", "\n", "Na aula anterior ([3. Refração crítica](3-refracao-critica.ipynb)), vimos o que acontece quando uma onda P incide sobre uma interface no ângulo crítico. A onda que refrata a 90° possui algumas características diferentes das ondas refletida e direta. A sísmica de refração é um método que utiliza essa características para estimar a profundidade da interface e as velocidades dos meios.\n", "\n", "Utilizaremos as simulações de ondas da biblioteca [Fatiando a Terra](http://www.fatiando.org). Essas simulações utilizam o [método de diferenças finitas](http://en.wikipedia.org/wiki/Finite_difference_method) para calcular soluções da equação da onda." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Objetivos\n", "\n", "* Visualizar como seria um dado sísmico de refração (simulado).\n", "* Observar os diversos caminhos que a onda faz para chegar nos receptores e como essas chegadas aparecem nos dados.\n", "* Aprender como usar as medições de primeira chegada para calcular a profundidade da interface e as velocidades." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Questão para entregar\n", "\n", "
\n", "Sua resposta deve conter no mínimo:\n", "
\n", "Failed to display Jupyter Widget of type interactive.
\n", " If you're reading this message in the Jupyter Notebook or JupyterLab Notebook, it may mean\n", " that the widgets JavaScript is still loading. If this message persists, it\n", " likely means that the widgets JavaScript library is either not installed or\n", " not enabled. See the Jupyter\n", " Widgets Documentation for setup instructions.\n", "
\n", "\n", " If you're reading this message in another frontend (for example, a static\n", " rendering on GitHub or NBViewer),\n", " it may mean that your frontend doesn't currently support widgets.\n", "
\n" ], "text/plain": [ "interactive(children=(FloatSlider(value=5.0, description=u'coef. linear (s)', max=54.5945945945946, step=0.5), FloatSlider(value=0.0001, description=u'coef. angular (s/m)', max=0.0005, min=-0.0005, readout_format=u'.6f', step=1e-05), Output()), _dom_classes=('widget-interact',))" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "dados = simul[:, 0, receptores]\n", "times = np.linspace(0, simul.dt*simul.simsize, simul.simsize)\n", "\n", "def reta(linear, angular):\n", " x = receptores*spacing\n", " y_reta = linear + angular*x\n", " plt.figure(figsize=(12, 8))\n", " ax = plt.subplot(111)\n", " ax.pcolormesh(x, times, dados, cmap='Greys', vmin=-50, vmax=50)\n", " ax.plot(x, y_reta, '-b', linewidth=2)\n", " ax.set_xlim(x.min(), x.max())\n", " ax.set_ylim(0, times.max())\n", " ax.set_ylabel('tempo (s)')\n", " ax.set_xlabel(u'Distância (m)')\n", "\n", "ipw.interactive(reta, \n", " linear=ipw.FloatSlider(min=0, max=times.max(), step=0.5, value=5, description='coef. linear (s)'),\n", " angular=ipw.FloatSlider(min=-0.0005, max=0.0005, step=0.00001, value=0.0001,\n", " readout_format='.6f', description='coef. angular (s/m)'))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Para pensar\n", "\n", "* Como aparecem nos dados as ondas direta e refratada no ângulo crítico?\n", "* Tente encontrar o ponto onde a refrata ultrapassa a direta.\n", "* Você pode usar a reta azul para achar os coeficientes angular (em s/m) e linear (em s) de eventos no dado.\n", "* Tente usar os valores obtidos com a reta azul para calcular V1 (de cima), V2 (de baixo) e a profundidade da interface (Moho)." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## License and information\n", "\n", "**Course website**: https://github.com/leouieda/geofisica2\n", "\n", "**Note**: This notebook is part of the course \"Geofísica 2\" of Geology program of the \n", "[Universidade do Estado do Rio de Janeiro](http://www.uerj.br/). \n", "All content can be freely used and adapted under the terms of the \n", "[Creative Commons Attribution 4.0 International License](http://creativecommons.org/licenses/by/4.0/).\n", "\n", "" ] } ], "metadata": { "anaconda-cloud": {}, "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.14" }, "widgets": { "state": { "07bb18f03c5a41dc9ed6a81d8c71cf5d": { "views": [ { "cell_index": 20 } ] } }, "version": "1.2.0" } }, "nbformat": 4, "nbformat_minor": 1 }