{"nbformat":4,"nbformat_minor":0,"metadata":{"colab":{"name":"template_instruções_Ex1_modelagem_tx_2022.ipynb","provenance":[{"file_id":"1h2j1Vraei7l2WBlzYHhV4vvO_ARDWvE2","timestamp":1661825626216}],"collapsed_sections":[]},"kernelspec":{"name":"python3","display_name":"Python 3"}},"cells":[{"cell_type":"markdown","metadata":{"id":"0_fXg6U1h547"},"source":["**DESCREVA AQUI: Informações sobre este programa; o autor (seu nome); a data; e \n","o que mais achar importante**"]},{"cell_type":"markdown","source":["Sejam as equações de tempo-distância (t(x)), respectivamente, para as ondas direta e refratadas na primeira e segunda interfaces planas:\n","$$ t_1(x) = \\frac{x}{V_1} $$\n","\n","$$ t_2(x) = \\frac{2h_1cos(i_{12})}{V_1} + \\frac{x}{V_2} $$\n","\n","$$ t_3(x) = \\frac{2h_1cos(i_{13})}{V_1} + \\frac{2h_2cos(i_{23})}{V_2} + \\frac{x}{V_3} $$\n","\n"],"metadata":{"id":"n4qvQ04XjcSn"}},{"cell_type":"markdown","source":["**1.1 COMPLETE a descrição sobre qual é o significado dos termos:**\n","\n","  X   \n","\n","  V1, V2 e V3  \n","\n","  h1 e h2\n","\n","  i12  \n","\n","  i13  \n","\n","  i23  \n","\n","$ t_{i_1} = \\frac{2h_1cos(i_{12})}{V_1} $\n","\n","$ t_{i_2} = \\frac{2h_1cos(i_{13})}{V_1} + \\frac{2h_2cos(i_{23})}{V_2} $\n"],"metadata":{"id":"S-l4_BeMlXI2"}},{"cell_type":"markdown","source":["1.2) Do ponto de vista analítico (físico-matemático), quais são as variáveis destas equações (ou seja, que não dependem do modelo geológico)?\n","\n","R."],"metadata":{"id":"Tj6zgJITmHXk"}},{"cell_type":"markdown","source":["1.3) Quais dos termos acima são os parâmetros que descrevem o modelo geológico?\n","\n","R."],"metadata":{"id":"oyBcXShomMrH"}},{"cell_type":"markdown","source":["1.4) Como cada um dos ângulos iab são obtidos? ( i12, i23, i13 )\n","\n","R."],"metadata":{"id":"WUVeqwe3mQ8n"}},{"cell_type":"markdown","metadata":{"id":"S_I0C4lRA3HM"},"source":["**1.5) Na célula a seguir, defina os parâmetros do modelo geológico para duas interfaces**\n","(são variáveis de entrada para o programa). \n","\n","Procure escolher valores que tenham significado geológico\n","(se utilizou a informação de algum trabalho ou tabela para definir o modelo, cite a referência)."]},{"cell_type":"code","metadata":{"id":"pX_Da6tTQOvf"},"source":["v1=\n","v2=\n","v3=\n","h1=\n","h2="],"execution_count":null,"outputs":[]},{"cell_type":"code","metadata":{"id":"MAgwRtPHmnfw","executionInfo":{"status":"ok","timestamp":1661828419044,"user_tz":180,"elapsed":26,"user":{"displayName":"Liliana Alcazar Diogo","userId":"17867606949541679083"}}},"source":["import numpy as np\n","import matplotlib.pyplot as plt\n","from google.colab import files "],"execution_count":1,"outputs":[]},{"cell_type":"markdown","metadata":{"id":"hg9Uv5FjBTmJ"},"source":["Na próxima célula do código:\n","\n","- Defina a geometria de aquisiçao:\n","  - afastamento mínimo fonte receptor (xmin)\n","  - intervalo entre geofones (dx)\n","  - Numero de geofones (NG)\n"," \n","- Crie o vetor (x) das distâncias fonte-receptor com o comando:\n","  **np.arange**    ou    **np.linspace**\n","\n","- Confira se o tamanho do vetor é igual ao NG que você definiu, para tal:\n","  print (len(x))"]},{"cell_type":"code","metadata":{"id":"OQIY-7zrPrgg"},"source":["xmin=\n","dx=\n","NG= \n","\n","x=np.()\n","print (len(x))"],"execution_count":null,"outputs":[]},{"cell_type":"markdown","metadata":{"id":"RmGD8YARRFJg"},"source":["Na  próxima célula, programe as fórmulas da distância crítica para cada uma das interfaces:\n","\n","$$Xc_1 = 2 h_1 tan(i_{12}) $$\n","\n","$$Xc_2 = 2 h_1 tan(i_{13}) + 2 h_2 tan(i_{23}) $$\n"]},{"cell_type":"code","metadata":{"id":"fZlbwLiuBOR6"},"source":["xc1=\n","xc2="],"execution_count":null,"outputs":[]},{"cell_type":"markdown","metadata":{"id":"YPlCzJIqRXAK"},"source":["A seguir, programe as fórmulas das equações de tempo-distância das ondas direta e refratadas:\n","\n","$$ t_1(x) = \\frac{x}{V_1} $$\n","\n","$$ t_2(x) = \\frac{2h_1cos(i_{12})}{V_1} + \\frac{x}{V_2} $$\n","\n","$$ t_3(x) = \\frac{2h_1cos(i_{13})}{V_1} + \\frac{2h_2cos(i_{23})}{V_2} + \\frac{x}{V_3} $$\n"]},{"cell_type":"code","source":["#Onda Direta\n","t1="],"metadata":{"id":"-fUXa4HVqIDA"},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":["**Para as Ondas Refratadas, o DESAFIO é calcular os tempos das ondas refratadas só para x maior do que a distância crítica**"],"metadata":{"id":"DHJM_lYJqPXt"}},{"cell_type":"code","metadata":{"id":"D-REOKEjR8RB"},"source":["#Ondas Refratadas\n","#Uma opção é selecionar os valores de x maiores do que a distância crítica de cada onda refratada\n","# e calcular os tempos com esses vetores x modificados (xr1) e (xr2)\n","xr1=\n","xr2=\n","\n","t2=\n","t3="],"execution_count":null,"outputs":[]},{"cell_type":"markdown","metadata":{"id":"ZJoJYpwQS1aG"},"source":["**Para gerar o gráfico das curvas de tempo-distância, complete as linhas na célula a seguir**"]},{"cell_type":"code","metadata":{"id":"6f0O73ggSWWP"},"source":["plt.plot(          ,label='')\n","plt.plot(          ,label='')\n","plt.plot(           ,label='')\n","plt.xlabel('')\n","plt.ylabel('')\n","plt.legend()\n","plt.grid()\n","plt.show()\n"],"execution_count":null,"outputs":[]},{"cell_type":"markdown","metadata":{"id":"DeWd-DGfSHLI"},"source":["Para os próximos Exercícios, precisaremos importar esses valores de tempo calculados.\n","Então, salve os tempos calculados em arquivos de texto: \n","\n","Complete o nome do arquivo (\" \") nas linhas da célula a seguir.\n","\n","Para fazer o download do que foi salvo com o np.savetxt, use o comando:\n","\n","`files.download('nome do arquivo')`"]},{"cell_type":"code","metadata":{"id":"_L83eqgKS9Fp"},"source":["x,t=np.savetxt(\"nome_do_arquivo\", np.c_[x,t1])\n","np.savetxt(\" \", np.c_[xr1,t2])\n","np.savetxt(\" \", np.c_[xr2,t3])"],"execution_count":null,"outputs":[]}]}