Statistics of Analytical solution for multi-borehole heat exchangers field including discontinuous and heterogeneous heat loads