<jupyter_start><jupyter_text>Modelagem 2D de uma bacia sedimentar  **[Referência bibliográfica]*** <NAME>., <NAME>, and <NAME> (1959), Rapid Gravity Computations for Two-Dimensional Bodies with Application to the Mendocino Submarine Fracture Zone, J. Geophys. Res., 64(1), 49-59, doi:10.1029/JZ064i001p00049  Importando as bibliotecas<jupyter_code>import numpy as np
import matplotlib.pyplot as plt
import talwani as tw<jupyter_output><empty_output><jupyter_text>Calculando as profundidades verdadeiras da bacia através de uma função matemática<jupyter_code>xs = np.arange(0, 100000,100)[::-1]
depths = (-1e-15*(xs - 50000)**4 + 7000. - 2500*np.exp(-(xs - 70000)**2/(10000**2)))
depths -= depths.min()<jupyter_output><empty_output><jupyter_text>Visualização da bacia<jupyter_code>plt.figure(figsize=(20,5))

plt.xlabel('x (m)', fontsize = 16)
plt.ylabel('Depth (m)', fontsize = 16)
plt.fill_between(xs,depths,max(depths),color='black')
plt.fill_between(xs,depths,color='orange')
plt.xlim(min(xs),max(xs))
plt.ylim(max(depths),min(depths))
plt.xticks(fontsize = 12)
plt.yticks(fontsize = 12)

file_name = 'images/basement2D_relief'
plt.savefig(file_name+'.png',dpi=300)

plt.show()<jupyter_output><empty_output><jupyter_text>Gerando as coordenadas de observação<jupyter_code>xp = np.arange(0,100000,step=500)
xp.size
zp = -100.*np.ones_like(xp)<jupyter_output><empty_output><jupyter_text>Calculando o efeito gravitacional ($g_z$) da bacia sedimentar considerando um contraste de densidade $\Delta\rho = -300 \, kg/m^3$<jupyter_code>gz = tw.gz(xp,zp,xs,depths,density=-300.)<jupyter_output><empty_output><jupyter_text>Visualização dos dados observados<jupyter_code>fig, (ax1, ax2) = plt.subplots(2, sharex=True,figsize=(20,10))

## Efeito gravitacional da bacia
ax1.set_title('Basin 2D modelling')
ax1.set_ylabel('$g_z$ (mGal)', fontsize = 16)
ax1.plot(xp,gz,'k.')

## Bacia 2D
ax2.set_xlabel('x (m)', fontsize = 16)
ax2.set_ylabel('Depth (m)', fontsize = 16)
ax2.fill_between(xs,depths,max(depths),color='white')
ax2.fill_between(xs,depths,color='gray')
ax2.plot(xs,depths,'b.')
ax2.set_xlim(min(xs),max(xs))
ax2.set_ylim(max(depths),min(depths))

file_name = 'images/basement2D_relief_grav'
plt.savefig(file_name+'.png',dpi=300)


plt.show()<jupyter_output><empty_output>
{'fraction_non_alphanumeric': 0.16424482606781154, 'fraction_numerical': 0.05856450902686042, 'mean_word_length': 8.088, 'pattern_counts': {'":': 0, '<': 28, '<?xml version=': 0, '>': 28, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'chain_length': '6', 'dir': 'jupyter-structured-clean-dedup', 'id': '402731', 'license': 'CC-BY-4.0', 'n_tokens_mistral': 1045, 'n_tokens_neox': 935, 'n_words': 151, 'path': 'Content/code/3. Grav_Mag_modeling/3.1. Basin_2D_modeling/1. 2D_grav_bacia_sedimentar.ipynb', 'repo_name': 'andrelreis/metodos-potenciais', 'size': '2287'}