Dans une démarche de modélisation hydrogéologique, la méthode des éléments analytiques présente une alternative intéressante par rapport aux méthodes numériques usuellement employées pour l'hydrodynamique souterraine. C'est une méthode aux frontières qui a été développée spécifiquement pour la modélisation des écoulements souterrains. Après avoir identifié sa place parmi les autres méthodes existantes et discuté de ses avantages et inconvénients, il s'agit d'approfondir cette comparaison sur un cas d'étude réel. C'est pourquoi, cette méthode est comparée à la méthode des différences finies à travers une méthodologie envisageant toutes les étapes de la modélisation ; de la création du modèle conceptuel jusqu'aux résultats de simulation. Cette comparaison fait apparaître l'importance de l'étape de la discrétisation et l'importance cruciale qu'elle représente pour l'utilisateur sur la qualité générale de la modélisation. Après avoir identifié les approximations inhérentes à la méthode des éléments analytiques, une analyse mathématique de problèmes simplifiés et localisés sur quelques éléments analytiques permet de quantifier ces erreurs et d'y extraire une analyse de leur comportement numérique. Cela permet d'introduire par la même occasion des éléments de théorie de la méthode des éléments analytiques. La représentation objet et vectorielle des éléments analytiques permet de bien valoriser les interactions possibles avec les SIG. En particulier, il s'agit d'élaborer des algorithmes de discrétisation adaptés aux spécificités et aux besoins de la méthode, mais aussi de mettre en valeur des critères d'évaluation de la qualité de cette discrétisation. La synthèse de ces outils et méthodes permet à l'utilisateur de disposer d'un outil d'aide à la discrétisation pour réduire au mieux ses effets sur la qualité de modélisation. Il s'agit ensuite d'étudier le comportement numérique des éléments analytiques adaptés à la simulation en régime transitoire. Pour cela, des éléments analytiques 1D sont d'abord développés pour faciliter cette étude. Leur comportement numérique reste similaire tout en réduisant les erreurs numériques qui sont induites par le calcul des éléments analytiques euxmême. Cette analyse permet de mettre en évidence une méthode de résolution adaptée au régime transitoire qui est ensuite transposée au problème 2D. Elle y expose en particulier une méthode pour estimer le champ d'influence de l'élément transitoire, utilisant un critère de discrétisation spatial qui doit être corrélé avec le pas de discrétisation temporel.