Random media emerge in several applications in reactor physics and safety analysis. Most often, models of stochastic media assume spatial homogeneity, whereas real-world complex materials, such as fuel chunks resulting from core degradation, typically display apparent heterogeneities. In a series of previous works, we have shown that stochastic tessellations can be successfully used in order to describe the material properties of several classes of random media. In this paper we extend these results to the case of heterogeneous random media by using Voronoi tessellations with space-dependent seed distributions, allowing for spatial gradients.