A two dimensional analysis program for groundwater flow in fractured network was developed to analyze the influence of discontinuity characteristics on groundwater flow. This program involves the generation of discontinuities and also connectivity analysis. The discontinuities were generated by the probabilistic density function(P.D.F.) reflecting the characteristics of discontinuities. And the fracture network model was completed through the connectivity analysis. This program also involves the analysis of groundwater flow through the discontinuity network. The result of numerical experiment shows that the equivalent hydraulic conductivity increased and became closer to isotropic as the density and trace length increased. And hydraulic head decreased along the fracture zone because of much water-flow. The grouting increased the groundwater head around cavern. An analysis of groundwater flow through discontinuity network was performed around underground oil storage cavern which is now under construction. The probabilistic density functions(P.D.F) were obtained from the investigation of the discontinuity trace map. When the anisotropic hydraulic conductivity is used, the flow rate into the cavern was below the acceptable value to maintain the hydraulic containment. But when the isotropic hydraulic conductivity is used, the flow rate was above the acceptable value.