La simulation numérique de la propagation des ondes électromagnétiques sur de longues distances et au-dessus de terrain a ces dernières années reçu une attention particulière du fait de son fort impact sur les systèmes radar et de télécommunications. Habituellement, il est considéré une modélisation bidimensionnelle pour traiter ces problématiques, cependant par une telle approche il est impossible de considérer les effets transverses au plan vertical passant par l'émetteur et le récepteur ainsi que la dépolarisation de l'onde. Pour pallier ces problèmes, une approche tridimensionnelle doit obligatoirement être considérée.<br /><br />La méthode de modélisation proposée est basée sur l'Equation Parabolique 3D (EP3D). Deux résolutions de celle-ci ont été considérées : nommées Split-Step Fourier (SSF) et Différences Finies (DF). La résolution SSF est basée sur une décomposition en un spectre angulaire d'ondes planes par l'intermédiaire d'une transformée de Fourier. La résolution de l'EP3D par DF développée utilise quant à elle un algorithme dit de Crank-Nicholson. Afin d'optimiser le temps de calcul et l'espace mémoire nécessaire, la méthode des directions alternées a été appliquée pour résoudre cette équation de propagation. Toutes deux ont été couplées avec la condition aux limites de Léontovich pour pouvoir prendre en compte le relief 3D.<br /><br />Ces deux méthodes ont été implémentées et validées sur différents cas tests canoniques. On a ainsi pu constater la capacité de ces méthodes à modéliser les phénomènes de réflexion, diffraction et réfraction. Celles-ci ont ensuite été appliquées au-dessus de scènes tridimensionnelles réalistes. Ces applications ont permis de comparer les deux méthodes développées ainsi que de mettre en relief les effets 3D dus au terrain et souligné les avantages d'une résolution tridimensionnelle.