Les travaux de cette thèse portent sur la modélisation et la discrétisation des écoulements diphasiques dans les milieux poreux fracturés. On se place dans le cadre des modèles dits dimensionnels hybrides couplant l'écoulement dans la matrice 3D à l'écoulement dans un réseau de fractures modélisées comme des surfaces 2D. La discrétisation s'appuie sur le cadre abstrait des schémas gradients. Dans cette étude nous présentons deux classes de schémas de types Gradient Schemes sur ces modèles en monophasique et en diphasique. Les objectifs sont motivés par l'application cible de la thèse qui concerne les procédés de récupération assistée de gaz par fracturation hydraulique dans les réservoirs de très faibles perméabilités