Fractures are widely existed in subsurface. They dominate the fluid flow in many types of media and have a huge impact on many subjects such as environmental engineering, groundwater engineering, geothermal energy and etc (Council, 2001;Rouchier et al., 2012;Watanabe et al., 1998). Hydraulic fracturing techniques are frequently used to improve the productivity of tight reservoirs to make the production of shale gas and oil possible (Medeiros et al., 2007;Roussel & Sharma, 2011). Hence, fractures dominate the geo-mechanical and hydrological behavior of many hydrocarbon reservoirs.A well known and successful method to characterize fractures is the Discrete Fracture Network (DFN) model. This technique was proposed in the 1980s and developed afterward with various applications in civil, environmental and reservoir engineering (Berre et al., 2018;Long & Billaux, 1987;Robinson, 1984). DFN models are a reconstruction of the geometry of real fracture networks which preserve topology and use regular shapes to represent real fractures. The commonly used shapes include rectangular, disc, ellipse etc. The geometry of real fractures are simplified as fracture aperture, length and width (rectangular) or major and minor axes (ellipse). Generally, geometrical properties of fractures can be obtained on the basis of observations and measurements and from techniques such as X-ray micro-Computed Tomography or downhole tools like high resolution resistivity logging and borehole radar techniques (