one of the important problem in better understand of Slope stability, Foundation dams, Underground excavations, Geothermal reservoirs, Evaluation of pollutants, Oil reservoirs and Nuclear waste disposal areas is simulation of rock mass fracture. Theoretical Studies has been demonstrated the result of three-dimensional modelling discrete fracture network is closer to reality than other models. There are different methodes to choose the parameters of fracture network for modeling. In general there are two methods for choosing the parameters: Statistics and Geostatistics. In statistics methods, there aren’t any information about the position of components and samples taken in space of the whole society. As a result, will not be the measured values of a specific quantity in a particular instance no information about value of same quantity at to a certain distance another case. In other words in statistical methods, fracture network modeling, generated without Spatial correlation structures, randomly and by Monte Carlo algorithm. Code GFracIUT Has been developed in order to three-dimensional modelling discrete fracture network to geostatistical approach. This program then produces microfractures to geostatistical approach checks Intersection or lack of intersection them together. After that eliminate subcritical microfractures. This code is capable of generate and display discrete fracture network for dimensions and realization different, with the distinction between categories of fractures. This program also are capable of fractures modeling with different geometric shapes Including disc, Parallelogram and generally polygon. This code in order to modelling fluid transmitter network in three-dimensional space, on the basis size, area and fractures geometric shape calculated the intersection of any fracture with other fractures or boundaries of the study area. Then, using depth-first search algorithm eliminates the clusters of fractures that have no connection with border or with only one the border are linked. Therefore identified network of fractures that effective on the transmission network and modeled are transmission network based on the model of the pipe network. Finally, using the discrete element method calculates the flow rate and conductivity of any fracture. Now by identification flow rate of any fracture be calculated output flow from any side and permeability in each direction.