A two-dimensional finite difference numerical model was developed in order to simulate two-phase fluid flow in a single fracture. In the model, variation of viscosity with pressure and that of relative permeability with water saturation can be treated. For the numerical solution, IMPES method was used, from which the pressure and the saturation of water and gas were computed one by one. Seven cases of model test using parallel plates for a single fracture were performed in order to obtain the characteristic equation of relative permeability which would be used in the numerical model. it was difficult to match the characteristic curves of relative permeability from the model tests with the existing emperical equations, consequently a logistic equation was proposed. As the equation is composed of the parameters involving aperture size, it can be applied to any fracture.