Ce travail est une contribution à l'étude numérique de la propagation des ondes sismiques dans les formations géologiques superficielles. Dans la première partie, un programme de calcul combinant les méthodes d'Aki-Larner et de Thomson-Haskell a été développé. Il permet de calculer la réponse de configurations géologiques géométriquement complexes à des ondes SH planes. Une étude systématique de quelques configurations géologiques-type montre que les effets dûs aux hétérogénéités verticales (contraste d'impédance: amplification) et latérales (irrégularités géométriques: diffraction) sont fortement couplés : il est le plus souvent impossible de décomposer le calcul d'une structure complexe en 2 sous-problèmes, plus simples. En génie parasismique, il est donc nécessaire de traiter chaque site cas par cas pour les études détaillées. Dans la deuxième partie, on étudie l'influence de la saturation en eau dans les formations poreuses très perméables. La représentation du milieu saturé est celle de Biot, avec quelques compléments issus de la théorie de l'homogénéisation. Ce modèle est appliqué au calcul de la réponse sismique de structures simples tricouches (sol poreux sec sur sol poreux saturé sur bedrock élastique imperméable) en 1 D ou 2D, dans la gamme 0 - 25 HZ. On montre des cas théoriques pour lesquels l'onde P2 (générée aux interfaces du milieu poreux, puis très atténuée et non observable en surface) intervient dans le mécanisme d'atténuation des ondes de compression. On insiste sur l'importance de l'interaction fluide-solide aux frontières du milieu poreux: dès que la perméabilité est supérieure à 10-10 m2 environ, il s'avère nécessaire de prendre en compte explicitement la présence de fluide. En particulier, il est alors impossible (sous peine d'erreurs supérieures à 100 %) d'approximer le matériau poreux par un modèle monophasique viscoélastique équivalent.