This numerical study investigates hydraulic fracturing and induced seismicity in intact and fractured reservoirs under anisotropic in situ stress using hydro-mechanical coupled discrete particles joints model. A 2 km x 2 km reservoir model with granitic rock and joints properties is constructed. Various injection scenarios are tested which involve continuous and cyclic styles of pressure controlled and flow rate controlled injections. Results are compared which include: spatial and temporal evolution of induced seismic events in relation with fluid pressure distribution, moment magnitudes of the induced events, occurrence of post-shut-in large magnitude events, etc. Several field observations on induced seismicity phenomena are simulated which include creation of new fractures, re-activation of the pre-existing joints, post-shut-in seismicity and large magnitude event with non-double-couple source, Kaiser phenomenon, moment magnitude vs. frequency distribution of the induced events following the Gutenberg-Richter law, etc. Cyclic injection results in larger volume of injected fluid but less number of total events and larger magnitude events; hence less seismic energy radiated by the induced events, slower relaxation of the fluid pressure after shut-in, longer and thinner propagated fractures with larger fluid saturated area. The major conclusions of this study are that the presented modeling is capable of simulating induced seismicity phenomena in Enhanced Geothermal System and fluid injection in fractured reservoirs in cyclic way has potential in mitigating the effects of larger magnitude induced events. (C) 2014 Elsevier Ltd. All rights reserved.