This contribution is devoted to modeling damage zones caused by the excavation of tunnels and boreholes (EDZ zones) in connection with the issue of deep storage of spent nuclear fuel in crystalline rocks. In particular, elastic-plastic models with Mohr-Coulomb or Hoek-Brown yield criteria are considered. Selected details of the numerical solution to the corresponding problems are mentioned. Possibilities of elastic and elastic-plastic approaches are illustrated by a numerical example.