A robust numerical method for modeling multiple wells in city-scale geothermal field based on simplified one-dimensional well model. (August 2019)