Esquemas Discretos en Diferencias Finitas para Modelar Propagaciones Sísmicas

Abstract

En este trabajo realizamos una comparación de técnicas en diferencias finitas para la solución numérica de la ecuación elástica 2D. Se obtuvo tres diferentes esquemas, usando el concepto de malla homogénea, intercalada y variable. Estos esquemas difieren en la disposición de los nodos en la malla, el costo, la dispersión numérica y estabilidad. Se obtuvieron registros de la propagación de ondas elásticas donde las propiedades de la roca varía. En estas propiedades se tuvo en cuenta la densidad, el móodulo de rigidez, la compresibilidad del medio y las velocidades de las ondas compreisionales y de corte. Finalmente se muestran resultados de estabilidad, dispersión y el campo de onda teniendo en cuenta las tres técnicas. El objetivo final de este trabajo es determinar o predecir de cada esquema su funcionalidad, y utilidad para el modelado de ondas sísmicas.