The impact of heterogeneity on geothermal production: simulation benchmarks and applications