Ground Source Heat Pumps (GSHPs) have emerged as a prevalent substitute for traditional heating and cooling systems in buildings due to their remarkable efficiency and energy-saving characteristics. GSHP systems are normally composed of a heat pump connected to a ground heat exchanger, through which heat is injected into or extracted from the subsurface. Currently, borehole heat exchangers (BHEs), which provide more stable thermal resources and require relatively limited land area, are the most commonly adopted configuration in GSHP systems. However, the thermal behavior of BHE fields is governed by complex and strongly coupled heat transfer processes, making accurate and computationally efficient simulation essential for system design and operation. Therefore, the objective of this study is to develop a fast and reliable simulation framework for BHE fields that can accurately capture their thermal response under practical operating conditions, and to further integrate advanced control and optimization strategies to enhance the overall operational performance, energy efficiency, and long-term sustainability of GSHP systems.