Commun. Comput. Phys., 16 (2014), pp. 799-816. A Two-Phase Flow Simulation of Discrete-Fractured Media using Mimetic Finite Difference Method Zhaoqin Huang 1, Xia Yan 1, Jun Yao 1*1 School of Petroleum Engineering, China University of Petroleum (East China), Qingdao 266580, P.R. China. Received 5 April 2013; Accepted (in revised version) 17 March 2014 Available online 8 July 2014 doi:10.4208/cicp.050413.170314a Abstract Various conceptual models exist for numerical simulation of fluid flow in fractured porous media, such as dual-porosity model and equivalent continuum model. As a promising model, the discrete-fracture model has been received more attention in the past decade. It can be used both as a stand-alone tool as well as for the evaluation of effective parameters for the continuum models. Various numerical methods have been applied to the discrete-fracture model, including control volume finite difference, Galerkin and mixed finite element methods. All these methods have inherent limitations in accuracy and applicabilities. In this work, we developed a new numerical scheme for the discrete-fracture model by using mimetic finite difference method. The proposed numerical model is applicable in arbitrary unstructured gridcells with full-tensor permeabilities. The matrix-fracture and fracture-fracture fluxes are calculated based on powerful features of the mimetic finite difference method, while the upstream finite volume scheme is used for the approximation of the saturation equation. Several numerical tests in 2D and 3D are carried out to demonstrate the efficiency and robustness of the proposed numerical model. AMS subject classifications: 76S05, 65N08, 65N55, 35J25 Notice: Undefined variable: pac in /var/www/html/readabs.php on line 165 Key words: Fractured porous media, discrete-fracture model, two-phase flow, mimetic finite difference method. *Corresponding author. Email: huangzhqin@upc.edu.cn (Z. Q. Huang), jsyanxia@163.com (X. Yan), RCOGFR_UPC@126.com (J. Yao)