Ce travail de thèse porte sur la modélisation d'écoulements hétérogènes de fluides à seuil. Dans une première partie nous proposons un modèle phénoménologique qui couple les propriétés mécanique du fluide à une loi de Maxwell généralisée. Le système simple d'équations obtenu décrit une phénoménolgie complexe. Nous établissons ensuite un deuxiéme modèle, de type élasto plastique afin d'appéhender l'écoulement à un niveau plus microscopique. Pour construire ce modèle, nous procédons d'abord à une étude d'élasticité fine, afin d'obtenir l'expression analytique exacte du champ de perturbation du à un événement plastique, dans un milieu élastique, homogène, isotrope, dans une géométrie infinie et dans une géométrie finie. Nous avons ainsi caractérisé les effets de taille finie et les effets de parois. Ensuite, nous choisissons une dynamique de la plasticité. Enfin, le modèle élasto-plastique est résolu numériquement dans une géométrie bipériodique et dans une géométrie confinée par des parois. A fort taux de cisaillement, l'écoulement est homogène, à faibles taux de cisaillement, une organisation spatio temporelle complexe est mise en évidence. En particulier, dans ce régime l'écoulement se manifeste par bouffées de plasticité dont la dimension (dans le temps et dans l'espace) diminue lorque le taux de cisaillement augmente, et une localisation intermittente de l'écoulement près des parois apparaît.