Ground heat storage systems can play an important role for the reduction of green house gases emissions by increasing the exploitation of renewable energy sources and “waste heat” with a consequent diminution of the use of fossil fuels. A ground heat storage consists in an array of vertical boreholes placed in such a way that promotes the mutual thermal interaction between the ground heat exchangers creating the necessary conditions required to effectively store and retrieve heat. Suitable modeling tools for the estimation of the thermal behavior of these systems are very important to build installations yielding economical performance compatible with what expected during the design phase. This thesis aims at giving a contribution in the development of the thermal modeling of borehole heat storage systems. The main objective is introducing in the modeling process a few features that are not usually considered in state of the art models, with the goal of improving the representation of the physical phenomena. These features are the mathematical description of the topology of the borehole heat exchangers network, and the modeling of borehole fields with arbitrarily oriented boreholes. The detailed modeling of the topology of the borehole heat exchangers is approached with a network model. The overall geothermal system is discretized into smaller systems called components. These are linked between each other in a network fashion to establish the logical relations required to describe a given boreholes connections arrangement. The method showed that the combination of a sufficient level of discretization of the system and of a network representation yields respectively the granularity and the flexibility required to describe any borehole field connections configuration. The modeling of non-vertical borehole fields is approached by developing a method for the calculation of g-functions for these configurations. The method is an extension of a recent work done by Cimmino on the computation of g-functions for vertical borehole fields. This modeling technique is based on describing boreholes as sets of stacked finite line sources and on the superposition principle. This approach requires the computation of response factors relative to couples of finite lines. A procedure for the fast computation of these response factors for the case of arbitrarily oriented lines is given. This yields computational performance that guarantees the practical feasibility of the methodology. The last part of the thesis deals with the modeling of the storage system from a broader perspective. The borehole field is considered as part of a larger system constituted by several interacting components (i.e. heat pump, building, etc.). Interactions play a key role in the resulting overall performance of these systems. The analysis of the mutual relations between building envelope and borehole field design is utilized as an example to highlight advantages and challenges of strategies yielding a more integrated design. / <p>QC 20150507</p>