This paper investigates bedrock heterogeneity by applying three different geophysical approaches, in order to study the long-term behaviour and the interaction between closed-loop geothermal systems. The investigated site consists of four boreholes equipped with geothermal pipes on the campus of University of Liege, Belgium. The first approach includes acoustic borehole imaging, gamma-ray logging and cuttings observation and results to a detailed fracture characterisation, rock identification and layer dip angle determination. The second approach consists of measuring the thermal conductivity of cuttings at the laboratory. Study of cuttings thermal conductivity measurements can contribute to bedrock heterogeneity knowledge concerning the transition of one formation to another and the layer dipping. The third approach is based on high-resolution temperature profiles, measured during the hardening of the grouting material and the recovery phase of a Distributed Thermal Response Test. Through this approach a correlation of the temperature profiles to the geological characteristics of the surrounding bedrock is identified. The analysis of this correlation can provide information on fractured zones, alternation of different rock types and layering dipping. This latter approach can be easily applied on closed-loop geothermal systems to characterise the bedrock and investigate its heterogeneity as well as contribute to the their long-term behaviour prediction and to the optimisation of their efficiency.