The development and deployment of shale gas formations around the world are relatively recent, starting in the United Stated in the late 1990. From the results obtained with the application of hydraulic fracturing as a method of stimulation, the study and evaluation of other prospects of shale gas in others places in the world was encouraged. However, the analysis, study and characterization of this type of reservoirs are difficult, because it must be taken into account several factors such as geology, mineralogy, petrophysics, geochemistry among others. Hydraulic fracturing is a complicated hydro-mechanical coupled process, with high difficulty degree especially in shale gas reservoir, where natural fractures exist. A numerical study is conducted to investigate the hydromechanical behavior of a natural fracture during fluid injection. UDEC (Universal Distinct Element Code) software based on discrete elements method was employed to numerical modeling development. UDEC has the ability to model the hydro-mechanical behavior of a fracture including phenomena like fracture enlargement, closure, slippage, and dilation under contact or separation condition. In this numerical investigation, numerical reservoir models, with random natural fractures through a distribution voronoi were created aiming to represent the network complexity of natural fractures present in shale gas formations. The numerical study results show that the hydromechanical system behavior is strongly dependent on the parameters variation such as in situ stress, fluid fracturing viscosity and fluid injection rate. Therefore, these results provide a better understanding of fracturing mechanisms and pressure response of a hydraulic fracturing treatment in a non- conventional naturally fractured reservoir.