Re´sume´-Une approximation efficace par XFEM pour e´coulements de Darcy dans les milieux poreux arbitrairement fracture´s -Les e´coulements souterrains sont fortement influence´s par la pre´sence de failles et de grandes fractures qui modifient la perme´abilite´du milieu agissant comme des barrie`res ou des conduits pour l'e´coulement. Une description pre´cise des proprie´te´s hydrauliques des fractures est donc essentielle pour la mode´lisation de la migration du pe´trole ou de l'exploitation de ressources non conventionnelles. Toutefois, la largeur de fractures est souvent petite par rapport a`la taille typique des e´le´ments du maillage. Pour re´soudre le proble`me par une approximation nume´rique obtenue sans raffinement du maillage, on remplace les fractures par des surfaces immerge´es dans la matrice poreuse. Par ailleurs, on veut permettre aux surfaces d'eˆtre inde´pendantes par rapport a`la structure du maillage, en manipulant les discontinuite´s à l'inte´rieur des e´le´ments par la me´thode des e´le´ments finis e´tendus (XFEM Extended Finite Element Method). On utilise une technique similaire pour obtenir un outil plus flexible, apte a`la re´solution du proble`me dans les re´gions d'intersection potentielles entre les fractures.Abstract -An Efficient XFEM Approximation of Darcy Flows in Arbitrarily Fractured Porous Media -Subsurface flows are influenced by the presence of faults and large fractures which act as preferential paths or barriers for the flow. In literature models were proposed to handle fractures in a porous medium as objects of codimension 1. In this work we consider the case of a network of intersecting fractures, with the aim of deriving physically consistent and effective interface conditions to impose at the intersection between fractures. This new model accounts for the angle between fractures at the intersections and allows for jumps of pressure across the intersection. This latter property permits to describe more accurately the flow when fractures are characterised by different properties, than other models that impose pressure continuity. The main mathematical properties of the model, derived in the two-dimensional setting, are analysed. As concerns the numerical discretization we allow the grids of the fractures to be independent, thus in general non-matching at the intersection, by means of the Extended Finite Element Method (XFEM), to increase the flexibility of the method in the case of complex geometries characterized by a high number of fractures.