This study presents a characterization of fractured basement reservoir by using statistical and fractal methods with outcrop data, seismic data, as well as FMI log data. In the statistical method, fracture intensity and length have been calculated from various outcrop data. The optimum statistical distribution functions of fracture length for outcrops have been identified with the use of discriminant equation derived from Crofton's theory. The Fisher distribution constant, representing the fracture orientation, has been computed from FMI log data. With the statistical values and distribution functions, a 3D fracture network system has been generated. The result shows that there is no distinction in orientation of the fracture network system, and it excellently matches with the outcrop data. In the fractal method, fractal dimensions of fracture length and strike for the seismic fracture network in areal distribution were calculated; a greater value in fractal dimension means that the fracture network system has intensive fractal characteristics. Meanwhile, vertical distribution and dip angle of the fracture system have been evaluated from FMI log data. The resulting 3D fracture system presents that the overall strike and distribution of the fracture system are excellently matched with those of seismic data.