En el presente trabajo se simula la propagación 3D de ondas elásticas usando el modelo de shale fracturada dado por Schoenberg y Helbig. Éste considera la superposición de la matriz de rigidez y de la compliancia de las fracturas, donde el movimiento de apertura o cerrado de las mismas es independiente de movimientos en las direcciones del plano de fractura. Además, los elementos de esta matriz dan cuenta de la densidad y parámetros geométricos de fracturas; por lo tanto, la estimación de la compliancia de las fracturas es relevante en la caracterización del reservorio y el plan de perforación para producción. El modelado se realiza usando Abaqus-CAE que resuelve la ecuación de onda elástica mediante la técnica de elementos finitos e implementando un procedimiento explícito. La validación de los resultados se lleva a cabo comparando con una solución analítica tridimensional.