Abstract: In this paper, the computational formulas of fractured networks model and continuously homogeneous field models of fractures intersected either orthogonally or non-orthogonlly are derived systematically. The program for numerical analysis methods is compiled and has been applied to practical examples and checked by physical models.After numerous computations for comparison with seepage flow in fractures and in porous media, the phreatic surface is much higher in fractures than that in porous medium. Consequently, in the study of seepage flow in fractured rock masses and the stability of rock slope, the formulas of fracture seepage should be adopted and it is believed more reliable.