f LIDAR-based 3D reconstruction and modelling of a flat-topped non-rimmed carbonate platform: Aptian, Maestrat Basin, Spain

Outcrop-scale reconstruction of depositional geometries and facies distribution of carbonate systems improves our knowledge on their heterogeneity distribution, stacking patterns and stratal architecture. The collected data and derived models can be used as analogues for characterizing and modelling potential subsurface reservoirs. Traditional sedimentological analyses in cropping out carbonate systems have limited accuracy depending on exposure conditions, accessibility or past erosive processes. On the other hand, there is a need to complement classical sedimentological approaches with quantitative characterizations and models of sedimentary bodies. In this respect, processing of three-dimensional (3D) point clouds captured by terrestrial LIght Detection And Ranging (LIDAR) technology combined with real-time kinematic global positioning system offers to field geologists the possibility to construct virtual 3D digital outcrop models (DOMs), which allow for more accurate analyses, reconstructions and quantification of the outcropping facies distribution than conventional digital terrain models. We present a LIDAR 3D DOM of an Aptian flat-topped non-rimmed carbonate platform margin from the western Maestrat Basin (Spain). The DOM served as a departing point to perform a 3D reconstruction that shows the relationship between depositional architecture and facies distribution of the carbonate system. The reconstruction not only highlights the value of digital outcrop models to characterize virtual attributes not observable in the outcrops due to the limitations of the 2D views of the exposures, but also allows to refine outcrop-scale sequence stratigraphic analyses. In addition, the 3D sequence stratigraphic approach obtained together with the 3D facies distribution model generated can be used as an analogue for the characterization of subsurface carbonate reservoirs with similar depositional profiles. The workflow of this study followed these steps: 1) Acquisition of the outcrop 3D point data set using a ground-based terrestrial LIDAR equipped with a differential GPS; 44 overlapping scans were needed to cover the entire outcrops of the flat-topped non-rimmed carbonate system characterized, each scan has associated a high-resolution photograph. 2) Mapping stratigraphic surfaces and pseudowells describing 5 lithofacies onto each individual photograph using a CAD-based tool, the mapping is carried directly onto the photographs because manipulating the images and interpreting the details is easier than directly digitizing onto the point-cloud. 3) The features mapped onto the photographs are projected into the corresponding point-cloud in order to georeference them. 4) Locally georeferenced individual point-clouds and attached interpretations were globally georeferenced by means of the UTM coordinates of each scan. 5) The stratigraphic boundaries mapped are used reconstruct the surfaces bounding stratigraphic units. 6) Population of the internal facies distribution conditioned to the pseudowells. This methodology allows to efficiently extracting information from point clouds, and resulted in the construction of a high-resolution 3D geological model displaying the stratal architecture and facies heterogeneity of sedimentary bodies, confined within a 3D sequence stratigraphical framework.