“…Recently, with the advance of modern computing power, three-dimensional numerical modeling become widely used such as in Dickinson et al [8], Lee [13], Tenma et al [28], and so on. The main objective of the numerical modeling for the ATES system before 1990s was to predict heat transfer in the model domain by injected water through a single well [16,18,26], and it is now widening into predicting recovery temperature by heat loss [13,16], injection depth in the aquifer [13,28], injection rate [13,28], injection time [28], and regional groundwater flow [17].…”