Numerical generation of the temperature response factors for a Borehole Heat Exchangers field