A numerical model for transient temperature distribution in an aquifer thermal energy storage system with multiple wells