Abstract: The long-term evolution of temperature field in high-level radioactive waste disposal repository is one of the important bases for the safety evaluation and design of buffer layer and the spatial optimization of disposal container. A two-dimensional axisymmetric thermal analysis model with two layers of disposal unit is established to reveal the near-field temperature evolution of multi-barrier system. The Laplace and finite Fourier sine transforms are applied to the governing equations for heat transfer to obtain the temperature solutions to the buffer layer and the surrounding rock layer in the Laplace domain. The corresponding solutions in time domain are obtained by numerical inversion of the Laplace domain solutions with the help of the Crump method. The correctness of the model is verified by comparing with the linear heat source solution and the numerical solution. The temporal and spatial distribution characteristics of the temperature field in the near-field multi-barrier system of the disposal unit are analyzed, and the effects of different parameters on the peak temperature of the buffer layer are discussed. The semi-analytical solution of the model is used to determine the minimum disposal container spacing and predict the results of the in-situ tests. The results show that the semi-analytical solution of the model can simulate the temperature change of buffer layer more accurately by comparing with the linear heat source solution. When the tunnel spacing exceeds 40 m, it has small effects on the peak temperature of buffer layer. When the tunnel spacing is 40 m, the container spacing is 7.7 m. The semi-analytical solution of the model can well predict the results of the in-situ test.