As underground caverns have many advantages such as safety and operation, they can also be used for gas storage purpose. When liquefied gas is stored underground, the cryogenic temperature of the gas affects the stability of the storage cavern. In order to store the liquefied gas successfully, it is essential to estimate the exact temperature distribution of the rock mass around the caverns. The main purpose of this study is the development of theoretical solution to be able to estimate the temperature distribution around storage caverns and the assessment of the solution. In this study, a theoretical solution and a conceptual model for estimating two and three dimensional temperature distribution around the storage caverns are suggested. Based on the multi-dimensional transient heat transfer theory, the theoretical solution is successfully derived by assuming the caverns shape as simplified geometry. In order to assess the theoretical solution, by performing numerical experiments with this multi-dimensional model, the temperature distribution of the theoretical solution is compared with that of numerical analysis. Furthermore, the effects of the caverns size are investigated.