A multiscale method for Discrete Fracture Modeling (DFM) using unstructured grids is developed. The fine-scale discrete system is obtained by imposing tetrahedron (triangular for 2D domains) shaped grid cells, while lower-dimensional fractures are imposed at the grid interfaces. The DFM approach is then used to describe the transmissibility coefficients for all the interfaces, including those with the lower-dimensional fractures. On this fine-scale discrete system, a new algebraic multiscale formulation is developed, which first imposes two sets of coarseand dual-coarse grids. The former grid is essential for conservative multiscale formulations, while the second one is used for the calculation of local multiscale basis functions. The coarse-scale partitioning of the fracture and matrix domains is flexible and totally independent. Moreover, the multiscale basis functions are constructed for both the matrix and fracture domains. By construction, the basis functions are a partition of unity. For 2D and 3D test cases, the performance of the multiscale method is systematically assessed. It is shown that the method (with no multiscale iterations) provides accurate results, even for complex fractured systems. The presented multiscale method is a promising framework for real-field application of DFM models.