Constraining the recharge area of a hydrothermal system in fractured carbonates by numerical modelling