Geostatistical Interpolation of Subsurface Properties by Combining Measurements and Models