Fluid flow behavior in a porous medium is a function of the geometry and topology of its pore space. The construction of a three dimensional pore space model of a porous medium is therefore an important first step in characterizing the medium and predicting its flow properties. A stochastic technique for reconstruction of the 3D pore structure of unstructured random porous media from a 2D thin section training image is presented. The proposed technique relies on successive 2D multiple point statistics simulations coupled to a multi-scale conditioning data extraction procedure. The Single Normal Equation Simulation Algorithm (SNESIM), originally developed as a tool for reproduction of long-range, curvilinear features of geological structures, serves as the simulation engine. Various validating criteria such as marginal distributions of pore and grain, directional variograms, multiple-point connectivity curves, single phase effective permeability and two phase relative permeability calculations are used to analyze the results. The method is tested on a sample of Berea sandstone for which a 3D micro-CT scanning image is available. The results confirm that the equi-probable 3D realizations obtained preserve the typical patterns of the pore space that exist in thin sections, reproduce the long-range connectivities, capture the characteristics of anisotropy in both horizontal and vertical directions and have single and two phase flow characteristics consistent with those of the measured 3D micro-CT image. (C) 2011 Elsevier Ltd. All rights reserved.