Statistics of Thermal impact on the excavation damage zone around a supported drift using the 2nd gradient model