573 Hachemi, A. and Remini, R. 2014. A quasi three-dimensional model for the simulation of seawater intrusion into coastal aquifers. Boletín Geológico y Minero, 125 (4): 573-583 ISSN: 0366-0176 A quasi three-dimensional model for the simulation of seawater intrusion into coastal aquifers A. Hachemi (1) and B. Remini (2) (1) National High School for Agronomy, Algiers 16000, Algeria, hachemi_abdelkader2001@yahoo.fr (2) Department of Water Sciences, Blida University, Blida 9000, Algeria, reminib@yahoo.fr ABSTRACT A fnite element software based on a quasi three-dimensional model to simulate seawater intrusion in coastal aquifers has been established. The mathematical model is the same as that developed by Sorek et al. (2001) however, in this study we have used the fnite element method with stabilisation of the convection-diffusion equation by the SUPG method (Streamline Upwind Petrov Galerkin). The code is applied to an unconfned aquifer. The hydraulic head and chlorine concentration contour maps, which also position the interface, are presented. Key words: coastal aquifer, intrusion, modeling, quasi three-dimensional, sea water. Un modelo cuasi-tridimensional para la simulación de intrusión marina en acuíferos costeros RESUMEN Este artículo trata sobre un programa computacional de elementos fnitos basado en un modelo cuasi-tri- dimensional para simular la intrusión marina en acuíferos costeros. El modelo matemático es el mismo al desarrollado por Sorek et al. (2001). Sin embargo, en este trabajo se ha utilizado el método de elementos fnitos con estabilización de la ecuación de convección-difusión por el método SUPG (Streamline Upwind Petrov Galerkin). El código se ha aplicado a un acuífero no confnado. Se muestran los mapas de isolíneas del nivel piezométrico y de la concentración en cloruros, los cuales indican la posición de la interfase. Palabras clave: acuífero costero, agua de mar, cuasi-tridimensional, intrusión, modelado. VERSIÓN ABREVIADA EN CASTELLANO Objetivos El objetivo de este trabajo es presentar un modelo areal 2D basado en la aproximación de Dupuit para simu- lar la intrusión marina en un acuífero costero no confnado. El modelo, que fue propuesto por Sorek et al. (2001), se resuelve utilizando un método de elementos fnitos. La ecuación de transporte de solutos se puede estabilizar utilizando el método SUPG (Streamline Upwind Petrov Galerkin). El código propuesto nos permite calcular las isolíneas de la hidráulica y de la concentración de cloruros así como la profundidad de la interfase entre agua dulce y agua salada. Introducción Hay dos enfoques para el estudio de la intrusión en acuíferos costeros. El primero se denomina modelo de interfase neta y supone que los dos fuidos, agua dulce y agua salada, son inmiscibles y están separados por