Hydrological characterization of cave drip waters in a porous limestone: Golgotha Cave, Western Australia

Cave drip water response to surface meteorological conditions is complex due to the heterogeneity of water movement in the karst unsaturated zone. Previous studies have focused on the monitoring of fractured rock limestones that have little or no primary porosity. In this study, we aim to further understand infiltration water hydrology in the Tamala Limestone of SW Australia, which is Quaternary aeolianite with primary porosity. We build on our previous studies of the Golgotha Cave system and utilize the existing spatial survey of 29 automated cave drip loggers and a lidar-based flow classification scheme, conducted in the two main chambers of this cave. We find that a daily sampling frequency at our cave site optimizes the capture of drip variability with the least possible sampling artifacts. With the optimum sampling frequency, most of the drip sites show persistent autocorrelation for at least a month, typically much longer, indicating ample storage of water feeding all stalactites investigated. Drip discharge histograms are highly variable, showing sometimes multimodal distributions. Histogram skewness is shown to relate to the wetter-than-average 2013 hydrological year and modality is affected by seasonality. The hydrological classification scheme with respect to mean discharge and the flow variation can distinguish between groundwater flow types in limestones with primary porosity, and the technique could be used to characterize different karst flow paths when high-frequency automated drip logger data are available. We observe little difference in the coefficient of variation (COV) between flow classification types, probably reflecting the ample storage due to the dominance of primary porosity at this cave site. Moreover, we do not find any relationship between drip variability and discharge within similar flow type. Finally, a combination of multidimensional scaling (MDS) and clustering by k means is used to classify similar drip types based on time series analysis. This clustering reveals four unique drip regimes which agree with previous flow type classification for this site. It highlights a spatial homogeneity in drip types in one cave chamber, and spatial heterogeneity in the other, which is in agreement with our understanding of cave chamber morphology and lithology. © Author(s) 201

Analog-based meandering channel simulation

Characterizing the complex geometries and the heterogeneity of the deposits in meandering river systems is a long-standing issue for the 3-D modeling of alluvial formations. Such deposits are important sources of accessible groundwater in alluvial aquifers throughout the world and also play a major role as hydrocarbons reservoirs. In this paper, we present a method to generate meandering river centerlines that are stochastic, geologically realistic, connected, and conditioned to local observations or global geomorphological characteristics. The method is based on fast 1-D multiple-point statistics in a transformed curvilinear domain: the succession in directions observed in a real-world meandering river (the analog) is considered as statistical model for multiple-point statistics simulation. The integration of local data is accomplished by an inverse procedure ensuring that the channels pass through a given set of locations while conserving the high-order spatial characteristics of an analog. The methodology is applied on seven real-world case studies. This work demonstrates the flexibility and the applicability of multiplepoint statistics outside the standard paradigm that considers the simulation of a 2-D or 3-D variable with spatial coordinates

High resolution multi-facies realizations of sedimentary reservoir and aquifer analogs

Geological structures are by nature inaccessible to direct observation. This can cause difficulties in applications where a spatially explicit representation of such structures is required, in particular when modelling fluid migration in geological formations. An increasing trend in recent years has been to use analogs to palliate this lack of knowledge, i.e., exploiting the spatial information from sites where the geology is accessible (outcrops, quarry sites) and transferring the observed properties to a study site deemed geologically similar. While this approach is appealing, it is difficult to put in place because of the lack of access to well-documented analog data. In this paper we present comprehensive analog data sets which characterize sedimentary structures from important groundwater hosting formations in Germany and Brazil. Multiple 2-D outcrop faces are described in terms of hydraulic, thermal and chemical properties and interpolated in 3-D using stochastic techniques. These unique data sets can be used by the wider community to implement analog approaches for characterizing reservoir and aquifer formations

Image synthesis with graph cuts: a fast model proposal mechanism in probabilistic inversion

Geophysical inversion should ideally produce geologically realistic subsurface models that explain the available data. Multiple-point statistics is a geostatistical approach to construct subsurface models that are consistent with site-specific data, but also display the same type of patterns as those found in a training image. The training image can be seen as a conceptual model of the subsurface and is used as a non-parametric model of spatial variability. Inversion based on multiple-point statistics is challenging due to high nonlinearity and time-consuming geostatistical resimulation steps that are needed to create new model proposals. We propose an entirely new model proposal mechanism for geophysical inversion that is inspired by texture synthesis in computer vision. Instead of resimulating pixels based on higher-order patterns in the training image, we identify a suitable patch of the training image that replace a corresponding patch in the current model without breaking the patterns found in the training image, that is, remaining consistent with the given prior. We consider three cross-hole ground-penetrating radar examples in which the new model proposal mechanism is employed within an extended Metropolis Markov chain Monte Carlo (MCMC) inversion. The model proposal step is about 40 times faster than state-of-the-art multiple-point statistics resimulation techniques, the number of necessary MCMC steps is lower and the quality of the final model realizations is of similar quality. The model proposal mechanism is presently limited to 2-D fields, but the method is general and can be applied to a wide range of subsurface settings and geophysical data types