A la frontière entre l'atmosphère neutre et l'espace, le plasma ionosphérique est le siège de réactions d'ionisation complexes. Le champ magnétique terrestre et les champs électriques induits causent des fluctuations spatiales et temporelles de la concentration électronique. Ces irrégularités ionosphériques entrainent des variations rapides de l'amplitude et de la phase des signaux radioélectriques les traversant, notamment aux hautes latitudes. Ce phénomène est appelé scintillation ionosphérique et il est particulièrement craint par la communauté utilisatrice d'applications GNSS qui nécessite une disponibilité et une intégrité optimales des signaux. Le travail présenté dans cette thèse propose une modélisation complète, à 3 axes d'anisotropie, de la scintillation ionosphérique. Ce modèle est basé sur une approche numérique 3D et 2D, de type écrans de phase, et sur la résolution analytique des équations de propagation, en 3D et en 2D. Ces dérivations originales des variances et des spectres de log-amplitude et de phase ont mis en relief les limites de validité d'un modèle numérique 2D. L'étude de sensibilité menée sur les variances et les spectres ouvre également des perspectives d'inversion des données GNSS pour remonter aux caractéristiques du milieu ionosphérique. / The ionospheric plasma is located at the border between neutral atmosphere and outer space and many complex ionization reactions occur inside this turbulent medium. The Earth magnetic field and induced electric fields cause rapid fluctuations of electron density, both spatially and temporarily. When crossing this turbulent layer, RF signals show fast variations of amplitude and phase, especially at high latitudes. This phenomenon is called ionospheric scintillation and it is particularly feared by air navigation using GNSS services, since it degrades the availability and the integrity of signals. This PhD dissertation presents a complete modeling of the effects of ionospheric scintillation, with 3 anisotropy axes. It is based on a numerical approach using the multiple phase screens technique, both in 3D and 2D schemes, and on the analytical resolution of electromagnetic propagation equation, also both in 3D and 2D configurations. The limits of use of a 2D numerical scheme have been outlined by these original formulations of phase and log-amplitude variances and spectra. This complete modeling associated with a sensitivity study on these variances and spectra opens up interesting perspectives on data inversion, in order to better estimate the physical characteristics of the ionospheric medium.