Three-dimensional analysis of tunnel face stability in spatially variable soils