This paper describes new techniques and algorithms for building consistent and accurate structural models of the subsurface in structurally and stratigraphically complex areas, and proposes an innovative workflow for iteratively building large-scale models in presence of sparse or noisy data.Contrary to traditional 3D modeling approaches, which consist in modeling first geological interfaces independently from each other before gluing them together to form a watertight representation of geological layers, the proposed technique creates first a 3D unstructured mesh representing the volume of interest, which is then sliced and diced by geological horizons and faults. It is based on a global interpolation method in which a scalar property representing the relative stratigraphic age of the formations is interpolated everywhere in the volume of interest. In the process, the geometry of each horizon is effectively constrained by interpretation data of all other horizons belonging to the same conformable sequence.The interpolation is practically insensitive to the complexity of the fault network (X faults, Y faults, thrusts) and yields a high resolution, 3D representation of geological layers, allowing the extraction of main horizons based on seismic input and of closely spaced intermediate horizon surfaces defined by well data.The proposed technique ensures the preservation of four key quality features:1. Geometrical accuracy: surfaces are as smooth as possible while fitting all input data geometry. 2. Topological correctness: all model contacts are perfect and watertight. 3. Stratigraphic consistency: horizon surfaces cannot cross each other; layering and truncation patterns dictated by the input stratigraphic column are honored everywhere in the model. 4. Structural coherency: variations of fault displacement should be consistent with their size; layer thickness changes should be as regular as possible.These benefits are illustrated through the step-by-step construction and refinement of a large scale geological model located offshore Western Australia. It is also demonstrated that this technique and methodology can successfully be applied to a large variety of structural and stratigraphic contexts.