Often, the unknown diffusivity in diffusive processes is structured by piecewise constant patches. This paper is devoted to efficient methods for the determination of such structured diffusion parameters by exploiting shape calculus. A novel shape gradient is derived in parabolic processes. Furthermore, quasi-Newton techniques are used in order to accelerate shape gradient based iterations in shape space. Numerical investigations support the theoretical results.