Cette thèse vise à étudier le phénomène dit de " cascades électromagnétiques cosmologiques ". Ces cascades sont typiquement générées dans le milieu intergalactique par l'absorption de rayons gamma sur les photons du fond optique / UV et par la production de paires électron / positron associés. Ces leptons eux-mêmes interagissent avec les photons du fond diffus cosmologique via diffusion inverse Compton pour produire de nouveaux rayons gamma qui eux même peuvent s'annihiler, générant à partir d'un unique photon primaire toute une gerbe de photons et de particules secondaires. D'un point de vue observationnel, le développement de cette cascade introduit trois effets : une déformation du spectre à haute énergie, un retard temporel dans l'arrivée des rayons gamma et une extension de la taille apparente de la source. Les cascades électromagnétiques cosmologiques ont commencé à être étudiées dans les années soixante. Mais ce n'est qu'à partir des années 2010 avec l'arrivée du satellite Fermi (entre autres) et des observations dans la bande au GeV et au TeV que la discipline a explosé. Le phénomène est particulièrement important. D'une part il altère le spectre observé des sources rendant difficile la compréhension de la physique de ces dernières. D'autre part les cascades se développant dans le milieu extragalactique, elles sont très sensibles à la composition de ce dernier (fond diffus de photons, champ magnétique). Or ce milieu étant très ténu, il est difficile à étudier. Les cascades deviennent alors une formidable sonde pour accéder à sa compréhension et pouvoir en comprendre l'origine qui remonte au commencement de l'Univers. Pourtant les cascades cosmologiques sont un phénomène complexe faisant intervenir des interactions difficiles à modéliser (sections efficaces complexes) et le transport de particules dans un Univers en expansion (cosmologie). Face à cette complexité les expressions analytiques sont vite limitées et le passage au numérique devient inévitable. Dans le cadre de cette thèse un code de simulation Monte Carlo a donc été développé visant à reproduire aussi précisément que possible le phénomène des cascades. Ce code a été testé et validé en le confrontant aux expressions analytiques. Grâce à ce code, le rôle des différents paramètres physiques impactant le développement de la cascade a été étudié de manière systématique. Cette étude a permis de mieux comprendre la physique du phénomène. En particulier, l'impact des propriétés du milieu extragalactique (fond diffus extragalactique, champ magnétique extragalactique) sur les observables a été mis en évidence. Finalement, une seconde étude a été menée pour mesurer la contribution des cascades au fond gamma extragalactique. Des travaux récents montrent qu'une grande partie de l'émission diffuse à très haute énergie provient de sources ponctuelles non résolues (blazars en particulier). Ces sources gamma (résolues et non résolues) doivent en principe initier des cascades qui peuvent contribuer au fond diffus. En partant d'une modélisation de l'émission des blazars à différents redshifts, l'absorption et la contribution des cascades ont alors été calculées à l'aide du code Monte Carlo. Les résultats montrent que la contribution des cascades au fond gamma extragalactique pourrait violer les limites Fermi mais l'excès doit encore être confirmé.