Heterogeneity determination of the Delft subsurface for heat flow modelling