In this article, new software has been developed for the simulation of naturally fractured reservoirs using discrete fracture modeling. Discrete fracture modeling in contrary to the double porosity modeling is able to simulate the real fracture distribution of a reservoir. In this modeling, non-ordered gridding has been used. Fractures are placed on the grids’ lines. Therefore, the geometry and the orientation of fractures are preserved. Also, in more heterogeneous region of a reservoir, smaller grids can be defined in order to increase the degree of accuracy in the calculation. In discrete fracture modeling, all connections of matrix-matrix, matrix-fracture, and fracture-fracture are used and the transfer function of double porosity model is not necessary and excluded. Also, a new approach has been used for the construction of triangular grids. This approach is the combination of Delaunay triangulation and the frontal advance method. For the numerical solution of equations, the finite volume method has been used. Finite volume method can be used in non-ordered grids with different size and geometry. The software in this work is developed for both 2-D and 3-D models and is tested by material balance for each grid. Also, this software is compared with commercial software in a case study. The accuracy of the developed software is confirmed.
)