Data from the World Health Organization indicate that Bulgaria has the second-highest COVID-19 mortality rate in the world and the lowest vaccination rate in the European Union. In this context, to find the crucial epidemiological parameters that characterize the ongoing pandemic in Bulgaria, we introduce an extended SEIRS model with time-dependent coefficients. In addition to this, vaccination and vital dynamics are included in the model. We construct an appropriate Cauchy problem for a system

more »

... of nonlinear ordinary differential equations and prove that its unique solution possesses some biologically reasonable features. Furthermore, we propose a numerical scheme and give an algorithm for the parameters identification in the obtained discrete problem. We show that the found values are close to the parameters values in the original differential problem. Based on the presented analysis, we develop a strategy for short-term prediction of the spread of the pandemic among the host population. The proposed model, as well as the methods and algorithms for parameters identification and forecasting, could be applied to COVID-19 data in every single country in the world.