We present a novel approach for 3D geostatistical simulations. It includes practical strategies for the development of realistic 3D training images, and for incorporating the diverse geological and geophysical inputs, together with their uncertainty levels (due to measurement inaccuracies and scale mismatch). Inputs consist of well-logs, seismics, and an existing 3D geomodel. The simulation domain (45 million voxels) coincides with the Miocene unit over 2810 km2
across the Danish-German border.