DOI: 10.18710/qnervt
title: Replication data for: Paleokarst reservoir modelling - A concept-driven approach

description: &lt;p&gt;Paleokarst reservoir model of the Agios Georgios cave system in northern Greece. A cave survey of the cave system is used to simulate collapse and infill processes and to generate paleokarst reservoir analogues. The collapse simulation is based on a continuous Tcrit evaluation throughout the cave system using the equation proposed in White, E. L., 2012, Breakdown, in White, W. B., and Culver, D. C., eds., Encyclopedia of Caves, Elsevier Science & Technology, p. 68-73. The simulation run is based on all input parameters being kept constant, except bedding dip. The input properties of the ceiling bed are:

&lt;p&gt;- Flexural strength: 12 MPa
&lt;p&gt;- Thickness: 5 m
&lt;p&gt;- Bedding dip: 0, 45 and 89 degrees


&lt;p&gt;The file comprises three RMS reservoir models representing different collapse scenarios (varying bedding dip). All models have a global grid cell resolution of 1 m^3, and upscaled versions with 8, 64, 512, 4096, 32768 m^3 grid resolution. Five different paleokarst facies are discretized: sediment-fill + fine chaotic breccia; open/intact cavern; coarse chaotic breccia; highly disturbed strata; disturbed strata.

&lt;p&gt;Readers are referred to Bjarte Lønøys 2021 PhD thesis "Paleokarst reservoir modelling - A concept-driven approach" for a comprehensive description of the workflow and applied parameter values.

&lt;p&gt;The workflow for geocellular rendering of cave surveys in paleokarst reservoir models is described in Lønøy, B., Tveranger, J., Pennos, C., Whitaker, F., & Lauritzen, S. E. (2020). Geocellular rendering of cave surveys in paleokarst reservoir models. Marine and Petroleum Geology, 122, 104652.

&lt;p&gt;The dataset comprises 566 files in 345 folders (Change view: tree structure). The folder structure follow the Emerson/Roxar RMS version 11 standard structure for reservoir models. 
Some of the files can be opened with common text processing softwares (e.g., notepad), whereas others may not. Thus, the files should be used in the Emerson/Roxar RMS software following the presented folder structure.
When the files are used in RMS, the workflow and applied parameters are described within the software. This implies that replication can be carried out following the workflow from stated in the files (as presented by the software).
description: &lt;p&gt;
&lt;p&gt;
ABSTRACT

A significant proportion of the world's hydrocarbon production comes from paleokarst reservoirs. Although these reservoirs boast some of the most productive wells in oil history, the recovery factor is relatively low (RFmean: 32%) compared to other carbonate reservoirs (RFmean: 37 - 51%). The low recovery could relate to current reservoir modelling approaches potentially yielding inaccurate resource estimates or early water-breakthrough. Conventional industry-standard reservoir modelling software suites do not have dedicated workflows or add-ins for handling the complex morphologies commonly associated with paleokarst. Current modelling approaches are often data-driven (conditioned on available seismic and well data) and employ adapted or modified versions of stochastic reservoir modelling workflows used for siliciclastic and carbonate reservoirs. However, many paleokarst features are below seismic resolution, and the representativity of individual well data is often challenging to assess. Consequently, data-driven models often fail to render the connectivity, geometry, and volume of karst features. Karst is the predecessor to paleokarst, and therefore a genetic approach employing existent information from recent karst systems may be a good starting point for generating analogues to paleokarst reservoirs. A concept-driven approach, in combination with current data-driven modelling approaches, may enable model rendering that more closely echoes actual paleokarst reservoir architectures. However, only a few conceptual modelling methods are publicly available and described in the literature. The drawbacks with the available methods are that they under-/overestimate the cave volumes, fail to provide realistic cave morphologies, and forecast clastic sediment infill, and do not differentiate between preserved and collapsed caverns. Consequently, post-collapse reservoir morphologies, volumes and facies distributions may be rendered inaccurately. This thesis aims to address the shortcomings of currently available conceptual methods and present a novel concept-driven workflow for paleokarst reservoir modelling.

&lt;p&gt;A novel methodology for geocellular rendering of karst systems is presented in this thesis. The method utilizes modern cave-survey data to generate dense, equally spaced point-clouds (infilling the cave periphery). These point clouds can be used to discretize the karst systems in a geocellular framework by geometrical modelling. The volumetric and geometric rendering of the method is compared with two pre-established methods and benchmarked against the cave survey. The results show that the new method offers improved volumetric and geometric geocellular rendering compared to the pre-established methods and are comparable to that of the cave survey.

&lt;p&gt;A pilot study using a well-known and pre-established geophysical method, electrical resistivity tomography (ERT), was carried out in the Maaras cave system in northern Greece to evaluate the large-scale volumetric significance and spatial distribution of clastic sediments infilling karst cavities. ERT proved to be a practical and useful method for differentiating mesoscale (&gt;2.5 m2) stratigraphic heterogeneity. Resistivity contrasts allowed the identification of sedimentary thickness variations, interbedded breccias, and cave floor. Results showed that the siliciclastic sediment thickness varied from 25 m to &gt;45 m, occupying a minimum of 69-95 % of the available accommodation space. 

&lt;p&gt;Finally, a novel interactive tool for evaluating cavern stability and forward model collapse and infill processes was developed. The tool employs conventional cave survey data, field measurements and geomechanical data of the host rock to simulate potential post-collapse morphologies and generate spatial output data suitable for geocellular modelling. Collapse propagation, and eventually the volume affected by the collapse, is controlled by user-defined paleokarst facies proportions and associated average porosities following a “mass-balance-principle” (i.e., porosity is final and only redistributed over a larger volume). Three different collapse scenarios were modelled using the Agios Georgios cave system in northern Greece as an analogue. The results show that it is feasible to use cave surveys to simulate collapse and infill processes and estimate the final paleokarst reservoir architecture. The morphology, volume and relative facies-proportions rendered in the reservoir models are comparable to those calculated in the forward collapse modelling tool, indicating that the geocellular model echoes the simulation. The results also show that the vertical continuity and target volume of a reservoir increases significantly with increasing bedding dip. This suggests that improved forecasting of the final reservoir architecture may optimise well positioning, production planning and eventually improve recovery prediction.
description: Emerson/Roxar RMS, 11