Les études des phénomènes d'interaction sismique sol-structure pour le dimensionnement d'ouvrages en génie civil sont fondées sur deux hypothèses simplificatrices fortes : l'homogénéité latérale du sol sous la fondation, et la représentation du mouvement sismique par des ondes planes à incidences verticales ou inclinées. Pour de grands ouvrages reposant sur des radiers flexibles étendus, ces hypothèses ne sont plus valides, d'autant que les observations in situ des ondes sismiques font apparaître une variabilité spatiale importante même sur de courtes distances et indépendamment de l'effet du passage d'onde. Dans ce travail sont développés les outils de modélisation et d'analyse numérique permettant de tenir compte de cette variabilité ainsi que de celle des caractéristiques mécaniques du sol. Une approche probabiliste est retenue, et leur influence est directement reliée à la réponse de la structure au séisme par une formulation intégrale incorporant ces deux aspects simultanément. Pour une mise en oeuvre numérique efficace, la dimension aléatoire des fluctuations - grandes ou petites - aléatoires des caractéristiques mécaniques du sol est réduite par l'introduction du développement de Karhunen-Loeve de l'opérateur de raideur dynamique associé. Cette technique est aussi appliquée au champ sismique incident. Les résultats obtenus pour des cas complexes réalistes permettent de mettre en évidence certains phénomènes dont l'appréhension apparaît indispensable dans le cadre d'une étude industrielle. Ils donnent notamment quelques indications utiles sur la sensibilité de la réponse de la structure à la variabilité spatiale du mouvement sismique ou des paramètres du sol.