Ce travail de thèse a porté sur le développement de méthodes numériques permettant de simuler la propagation des ondes sismiques "large bande" dans des milieux à géométries tridimensionnelles simples: un nombre réduit de discontinuités, surmontés d'une topographie 3D. Le champ incident modélisé doit pouvoir émaner d'une source ponctuelle proche de la surface. Dans l'introduction générale nous précisons les objectifs de cette étude, et passons en revue les méthodes numériques utilisées en sciences de la terre. Une méthode d'éléments de frontière reposant sur une représentation intégrale directe formulée en temps a été développée. Elle utilise un schéma de résolution explicite, qui permet de s'affranchir des limitations rencontrées en 3D par la formulation en fréquence et donne des résultats d'une précision très satisfaisante. Les limitations de la méthode viennent de la possible amplification du bruit numérique, mais surtout d'un temps de calcul important, qui limite en pratique la taille des modélisations réalisables. Les deux premiers chapitres exposent la méthode et ses résultats. Des expressions analytiques des fonctions de Green des milieux élastiques à gradient. de vitesse constant ont été proposées de façon heuristiques par Sanchez-Sesma et al. (1999), Dans le chapitre 3, nous avons testé ces fonctions pour le cas 3D, et cherché à déterminer leur applicabilité aux méthodes intégrales. Pour échapper aux limitations sur la taille des modèles imposées par les méthodes d'éléments de frontière, tout eu utilisant leur capacité à résoudre avec précision les conditions de surface libre, nous avons démarré l'étude d 'une méthode hybride combinant éléments de frontière et différences finies en fréquence. Le dernier chapitre propose essentiellement une étude de l'optimisation d' un schéma de différences finies d'ordre 2, et montre les premiers résu!tats de sa mise en oeuvre.