Cette thèse propose un formalisme général pour modéliser la propagation des ondes acoustiques dans une multicouche composée de toute combinaison de couches liquides, solides élastiques isotropes et poro-élastiques isotropes, la méthode ayant la flexibilité d'être développée pour inclure d'autres natures de couches. Dans un premier lieu, un algorithme stable est développé, basé sur l'approche récursive de la matrice de rigidité, pour modéliser la propagation d'une onde plane incidente sur la multicouche en fonction de son angle d'incidence et de sa fréquence. Cet algorithme fusionne de manière récursive les matrices de rigidité des couches individuelles de la structure en une matrice de rigidité totale et permet ensuite le calcul des coefficients de réflexion et de transmission, ainsi que les composantes de déplacement et de contrainte à l'intérieur de la multicouche pour chaque direction d'incidence des ondes planes. Deuxièmement, pour modéliser la propagation d'un faisceau délimité d'ondes incidentes, la technique du spectre angulaire est utilisée, basée sur la décomposition de ce faisceau en un spectre d'ondes planes se propageant dans des directions différentes. Par la suite, le faisceau d'onde réfléchi dans le milieu d'incidence et le faisceau d'onde transmis dans le milieu de transmission, ainsi que la distribution des champs (composantes de déplacement et de contrainte) à l'intérieur de la multicouche sont obtenus en superposant la contribution de toutes les ondes planes se propageant dans les différentes directions. Comme application numérique, une tri-couche solide-poreuse-solide immergée dans l'eau est simulée. La réflexion et la transmission qui en résultent, ainsi que les composantes de déplacement et de contrainte dans la multicouche, correspondants à l’onde plane incidente et au faisceau limité incident, révèlent la stabilité du procédé et la continuité des déplacements et des contraintes aux interfaces.