Borehole thermal energy storage (BTES) systems facilitate the subsurface seasonal storage of thermal energy on district heating scales. These systems’ performances are strongly dependent on operational conditions like temperature levels or hydraulic circuitry. Preliminary numerical system simulations improve comprehension of the storage performance and its interdependencies with other system components, but require both accurate and computationally efficient models. This study presents a toolbox for the simulation of borehole thermal energy storage systems in Modelica. The storage model is divided into a borehole heat exchanger (BHE), a local, and a global sub-model. For each sub-model, different modeling approaches can be deployed. To assess the overall performance of the model, two studies are carried out: One compares the model results to those of 3D finite element method (FEM) models to investigate the model’s validity over a large range of parameters. In a second study, the accuracies of the implemented model variants are assessed by comparing their results to monitoring data from an existing BTES system. Both studies prove the validity of the modeling approaches under investigation. Although the differences in accuracy for the compared variants are small, the proper model choice can significantly reduce the computational effort.