Carbonate reservoirs are a highly heterogeneous type of reservoir characterized by the presence of a large amount of vugs and pores. During two-phase displacement, the two-phase flow regime in the vugs might be gravity segregated. The distribution pattern of two-phase fluid in the vugs would accelerate the water flow in downward and horizontal directions, meanwhile decelerating in an upward direction, resulting in a different oil recovery ratio. This gives rise to the question of whether the relative permeability should be modeled as a directional dependent in a vugular porous medium since it is usually treated as an isotropic quantity. In this study, via both experiment and numerical simulation, we demonstrate that the relative permeability of vugular porous medium is dependent on the angle between the flow direction and the horizontal plane and should be considered for oil recovery estimation for carbonate reservoirs. Using the transmissibility-weighted upscaling method and a single-vug model, the relative permeability curves for different flow directions are obtained by numerical simulation. A directional relative permeability model for a vugular porous medium is also proposed.