Abstract The problem of identification of unknown epidemiological parameters (contagiosity, the initial number of infected individuals, probability of being tested) of an agent-based model of COVID-19 spread in Novosibirsk region is solved and analyzed. The first stage of modeling involves data analysis based on the machine learning approach that allows one to determine correlated datasets of performed PCR tests and number of daily diagnoses and detect some features (seasonality, stationarity, data correlation) to be used for COVID-19 spread modeling. At the second stage, the unknown model parameters that depend on the date of introducing of containment measures are calibrated with the usage of additional measurements such as the number of daily diagnosed and tested people using PCR, their daily mortality rate and other statistical information about the disease. The calibration is based on minimization of the misfit function for daily diagnosed data. The OPTUNA optimization framework with tree-structured Parzen estimator and covariance matrix adaptation evolution strategy is used to minimize the misfit function. Due to ill-posedness of identification problem, the identifiability analysis is carried out to construct the regularization algorithm. At the third stage, the identified parameters of COVID-19 for Novosibirsk region and different scenarios of COVID-19 spread are analyzed in relation to introduced quarantine measures. This kind of modeling can be used to select effective anti-pandemic programs.