L’objectif de ce travail de thèse est de prendre en compte à la fois l’hétérogénéité, l’anisotropie et des incertitudes dans la simulation 3D de la propagation d’ondes élastiques. Pour ce faire, dans un premier temps, on modélise le champ de propriétés mécaniques, ici le champ de tenseur d’élasticité, par un modèle de champ stochastique 3D des matrices définie-positives. La construction de ce modèle de champ est essentiellement fondée sur celle de Soize [2008]. Notre modèle conserve ainsi les propriétés principales du modèle de Soize comme le paramétrage minimal contrôlant l’amplitude de la fluctuation et la taille caractéristique de la variabilités patiale, le comportement local a priori arbitrairement anisotrope (anisotropie triclinique) et les propriétés mathématiques fondamentales. De plus, un nouveau paramètre est introduit dans ce modèle pour imposer un niveau d’anisotropie moyen souhaité. Dans un deuxième temps, on effectue des adaptations du code de calcul d’éléments finis spectraux, à savoir le code parallèle SPEC3D, afin d’une part de générer les réalisations du champ stochastique du tenseur d’élasticité et d’autre part de prendre en compte l’anisotropie dans la résolution numérique du problème élastodynamique. Des études paramétriques utilisant SPEC3D sont ensuite réalisées mettant en évidence les influences de l’anisotropie et des paramètres d’hétérogénéité sur la propagation d’ondes sismiques. En particulier, elles démontrent une dépendance directe entre la longueur de corrélation du champ de propriétés et le temps caractéristique d’apparition de la diffusion. Ce régime se manifeste par l’équipartition d’énergie entre les mouvements irrotationnels et rotationnels. / The aim of this thesis is to take into account the heterogeneity, the anisotropy and the uncertainties within 3D numerical simulation of elastic waves propagation. Firstly, the elasticity tensor field is modeled by means of a stochastic tensor-valued field. Its construction is generalized from the model of Soize [2008]. Hence, our model preserves principle properties of the former : a small set of parameters controlling the whole dispersion and the characteristic size of spatial variability, a local behavior being a priori arbitrary anisotropic (triclinic anisotropy) andothers essential mathematical properties. Moreover, a new parameter is added in order to impose a desired anisotropy mean level. Secondly, we carry out adaptations of an existing spectral finite elements-based elastic waves simulation software, namely the SPEC3D parallel computing code. On the one hand a sample generator of the elasticity random field model is implemented and on the other hand anisotropic material behavior is introduced in the elastodynamic solver. Finally, numerical parametric studies are performed using SPEC3D highlighting influences of heterogeneity and anisotropy on elastic waves behavior. In particular, it is observed that the characteristic time beyond which a multiple scattering pattern can be approximated by a diffusion regime directly depends on the correlation length of elasticity tensor field model. This time is detected by an energy equipartition between rotational and irrotational movements.