The article is devoted to Circle covering problem for a bounded set in a two-dimensional metric space with a given amount of circles. Here we focus on a more complex problem of constructing reserve and multiply coverings. Besides that, we consider the case where covering set is a multiply-connected domain. The numerical algorithms based on fundamental physical principles, established by Fermat and Huygens, are suggested and implemented. This allows us to solve the problems for the cases of non-convex sets and non-Euclidean metrics. Preliminary results of numerical experiments are presented and discussed. Calculations show the applicability of the proposed approach.