A Method of Fundamental Solutions in Poroelasticity to Model the Stress Field in Geothermal Reservoirs