<p>This study evaluates the equivalent absolute permeabilities and relative permeability curves of a dual-continuum model used to describe two-phase flow in highly heterogeneous two-dimensional porous media that contain fractures and voids embedded in a porous matrix. The Brinkman model, which combines Darcy equation to describe the flow in the porous matrix and Stokes equation to account for the free-flow region, is used to compute the evolution of velocity, pressure, and saturation fields for two-phase flow on the fine scale. The time-evolution of the pressure difference and water flow rate at the outflow boundary predicted by the Brinkman model are used as basis for the solution of optimization procedure to evaluate the equivalent petrophysical properties of the dual-continuum model. Results show that the use of linear relative permeability curves for the fractured domain, usually used in the literature, does not lead to accurate description of the flow. A workflow is presented to determine the relative permeability curve of the fracture domain that is able to accurately describe the flow dynamics.</p>