In this paper, we couple Discrete Fracture Network (DFM) and multi-continuum model with Generalized Multiscale Finite Element Method (GMsFEM) for simulating flow in fractured and vuggy reservoir. Various scales of fractures are treated hierarchically. Fractures that have global effect are modeled by continua while the local ones are embedded as discrete fracture network based on the geologic observation. For independent vugs, a continuum is used to represent their effects with specific configuration that there's no intra-flow of this continua. GMsFEM enables us to systematically develop an approximation space that contains prominent sub-grid scale heterogeneous background information based on the multi-continuum and DFM model. Conforming unstructured mesh is used to surrender the application of random discrete fracture networks. This paper targets on the improvement of the flow simulation performance in complex high-contrast domain by extending the ability of multiscale method to modeling arbitrary discrete fracture network. This advancement by GMsFEM is motivated by the limited capability of Multiscale Finite Element Method (MsFEM) on modeling discrete fractures when multiple fracture networks present in same coarse block. Multiple numerical results are shown to validate the efficiency of our coupled method.

You can access this article if you purchase or spend a download.