Seismic activity and related rockbursts are currently one of the most dangerous threats negatively affecting work safety and continuity of operation in Polish underground copper mines. Taking into account the experience to date, it can be clearly stated that current technical and organizational tools do not make it possible to completely eliminate mining tremors, however, by taking appropriate actions, it is possible to minimize the threat and partially control it, e.g. by active de-stressing the rock mass. However, developing an appropriate schedule of preventive activities is still a burning issue. This is due to the random nature of mining-induced seismic phenomena which makes it impossible to accurately predict the place, energy and time of the seismic event (Fuławka et al., 2018). A breakthrough could be the use of stage-based large-scale numerical modelling based on which it would be possible to track stress changes with ongoing exploitation and locate areas of increased seismic risk. This study presents the results of large-scale three-dimensional FEM-based numerical modelling enabling tracking of changes in stress distribution with the progress of exploitation. Then stress distribution and areas identified as prone to instability occurrence were correlated with the areas of high-energy seismic tremors manifestation. Models were prepared at monthly intervals, and validated with the use of measurements obtained with underground geomechanical monitoring systems. As preliminary analyses show, well-validated numerical models can be the basis for estimating seismic risk and may be useful at the stage of designing methods for preventing active rockburst and seismicity prevention.