coupled poroelastic stressing and pore-pressure accumulation along pre-existing faults in deep basement contribute to recent occurrence of seismic events at subsurface energy exploration sites. Our coupled fluid-flow and geomechanical model describes the physical processes inducing seismicity corresponding to the sequential stimulation operations in pohang, South Korea. Simulation results show that prolonged accumulation of poroelastic energy and pore pressure along a fault can nucleate seismic events larger than M w 3 even after terminating well operations. In particular the possibility of large seismic events can be increased by multiple-well operations with alternate injection and extraction that can enhance the degree of pore-pressure diffusion and subsequent stress transfer through a rigid and low-permeability rock to the fault. This study demonstrates that the proper mechanistic model and optimal well operations need to be accounted for to mitigate unexpected seismic hazards in the presence of the site-specific uncertainty such as hidden/undetected faults and stress regime. Over the past decade elevated levels of seismic activities have been observed at the sites related to subsurface energy exploration activities, including wastewater injection for conventional and unconventional oil/gas development 1-4 and geothermal stimulation 5-7. Despite the recent progress in statistical and physics-based investigation of induced seismicity 8-10 , recent unexpected moderate to large magnitude earthquakes (M 3 w ≥) after shut-in (e.g., 2006 M w 3.2 Basel, Switzerland 11 , 2017 M w 5.5 Pohang, South Korea 7) show the need of the mechanistic study to understand underlying physical mechanisms. Since fluid injection-extraction associated with these elevated earthquakes has often been operated with multiple wells, interactions of well operations and other hydrogeological features need to be further investigated. Large earthquakes require large seismogenic faults, and pressure perturbation and shear stressing are two primary factors to trigger fault slip by reducing fault strength 12-14. Pore-pressure accumulation along conductive faults has been considered as the principal mechanism for inducing seismicity 3,15 in which diffusive propagation of pressure plumes is essential, but controlled by hydraulic connectivity from faults to the fluid-injection reservoir. Another primary mechanism is the poroelastic stressing in which the volumetric changes of the pressurized zone perturb the stress field of the surrounding rock by transmitting elastic forces to longer distances even beyond the hydraulically affected region such as distant and disconnected basement faults 16-19. Temporal changes in stress states at frictional faults will determine the onset of fault slip corresponding to rate-and-state friction mechanisms 20,21 that can generate delayed surge of seismic events along faults even after shut-in. Site-specific features of geological formation and/or operational controls govern spatio-temporal patterns, rates, and magnitude...