As a key part of rock mass deformation and possible permeable channel，fractures seriously threaten the stability of rock mass and the safety of engineering such as in hydro-engineering project，road and bridge engineering and mining.Xishan mining area of the Sanshandao gold mine in Laizhou City，Shandong Province，China，is a coastal mine with complex geological structure and hydrogeological conditions.The fractures range from as small as micro-crack and as large as fault.And the overlying seawater continues to threaten the safety in mining.Therefore，the identification of the spatial distribution of cracks is an extremely important and indispensable basic work for studying the stability of rock masses in the mining area and the inrush water in the roadway.Combining with monitoring data of hydrological and fractures in Xishan mining area of Sanshandao gold mine for many years，the statistical method was used to describe the rules of geometric parameters such as inclination，dip and spacing of cracks in the main permeable roadways.The rose diagrams of crack strike，pole contour maps of fracture inclination and dip and frequency distribution histograms of each parameter were drawn.The distribution laws of fracture geometrical parameters were fitted.As such，three geometric models of the discrete fracture network of surrounding rock in the main permeable roadways were generated by using self-compiled programs in MATLAB according to the fitting results.The distribution of cracks can be displayed visually and dynamically.Statistics and three-dimensional modeling results show that：\mathrm{①}The main permeable levels include -375 m，-510 m，and -600 m in the Xishan mining area.Each level can be divided into two or three fracture groups according to fracture inclination.The tendencies of two major groups are northwest and southeast，respectively.The inclination of another group mainly developed in the -600 m level is between 0° and 10°.The dips of all the groups are large，mostly between 70° and 90°.\mathrm{②}According to the frequency distribution histogram of each parameter，the fractures could be fitted by the common probability distributions.The inclination of each fracture group obeys the logarithmic normal distribution or normal distribution.The dip of each fracture group obeys the normal distribution (values range from 0° and 90°).The spacing of each fracture group obeys the logarithmic normal distribution.The trace length of each fracture group obeys the negative exponential distribution.The statistical results provide the possibility for the quantitative studies of discrete fractures.\mathrm{③}The three-dimensional fracture network models based on the distribution law of fracture geometric parameters can reflect the characteristics of fractures realistically.The three fracture groups have close ties to the three fault zones named F1，F2 and F3，respectively.According to the characteristic of the three faults，it can be inferred that hydraulic conductivity of the fractures belonging to the tensile and torsional fault F3 are bigger than the other two.And the discrete fracture network models can provide essential references for the stability analysis of surrounding rock and the identification of inrush water passages.