Karst simulation with Lindenmayer-systems and ODSIM

International audienceKarstic systems are geological structures that strongly impact underground flows. Despite intensive explorations by speleologists, they remain partially described as many conduits are not accessible to humans. Paleokarsts are buried karstic systems with a significant reservoir potential. But they are not easily identifiable on seismic images. In those contexts, a huge uncertainty subsists on the network location and the conduit geometry. Stochastic simulations help to better assess that uncertainty. The difficulty is to reproduce the system connectivity at different scales while integrating as much geological knowledge as possible without involving poorly constrained parameters (e.g. paleo-climate, boundary conditions...). In this paper we propose to work on two aspects and scales of karstic systems. At large scale, we stochastically simulate karst network skeletons with a new method based on a formal grammar, the Lindenmayer-system. Based on an alphabet, an axiom and user-defined rules, the method puts together segments to build the network skeleton. The definition of proper rules and the introduction of karst-dedicated parameters generate curves reproducing the complex architectures encountered in those systems, mixing branchwork and anastomotic patterns. At the conduit scale, we propose to build a 3D envelope around these skeletons with an enhanced Object Distance based Simulation Method. It uses a custom distance field from the skeleton which takes into account geological features influencing karstogenesis (horizons, faults or fractures). This controls the first-order shape of the conduits. It is then combined to a custom random threshold controlling finer-scale features of the conduits. This threshold is generated with several parameter values depending on the involved geological structures. This workflow is demonstrated on a synthetic case, showing the potentialities of the approach at both scales. Results are encouraging and various improvements are in focus. Data conditioning, both to karst observations and local shape information has to be enhanced. The network simulation has the advantage to be grid-free, meaning that no background grid is needed to perform the simulation. Thus, it avoids the stair-step effect that can be observed in other techniques. On the opposite, the method used to simulate the conduit shapes relies on a grid, necessary to compute the distance fields and to perform the threshold geostatistical simulation. For detailed conduit geometry, the grid requires a high resolution, which impacts directly the computational efficiency. Finally, it would be interesting to test the approach on a real dataset and to develop a coupling with a flow simulator to evaluate the impact of the shape and of the network connections on the flow response

Modelling the long-term evolution of groundwater's quality in a flooded iron-ore mine using a reactive transport pipe network model

International audienceOver the past 25 years, underground mine flooding in the lorraine iron Basin (France) has resulted in a high concentration of dissolved sulphate and have made the water unsuitable for human consumption. this problematic issue has led to the development of numerical tools to support waterresource management in mining contexts. as water flows mainly in galleries and collapsed zones, we consider the flooded mine as a network of pipes and tanks. the software used for simulating flow andreactive transport in this network is the ePanet 2 code. a simplified sulphate dissolution-precipitation model, based on previous works, is included as source/sink in the tanks. Flow rates are calculated by processing data records with a rainfall-discharge model. the simulator gives good agreement between the calculated and observed sulphate concentrations

Comparing connected structures in ensemble of random fields

International audienceVery different connectivity patterns may arise from using different simulation methods or sets of parameters, and therefore different flow properties. This paper proposes a systematic method to compare ensemble of categorical simulations from a static connectivity point of view. The differences of static connectivity cannot always be distinguished using two point statistics. In addition, multiple-point histograms only provide a statistical comparison of patterns regardless of the connectivity. Thus, we propose to characterize the static connectivity from a set of 12 indicators based on the connected components of the realizations. Some indicators describe the spatial repartition of the connected components, others their global shape or their topology through the component skeletons. We also gather all the indicators into dissimilarity values to easily compare hundreds of realizations. Heat maps and multidimensional scaling then facilitate the dissimilarity analysis. The application to a synthetic case highlights the impact of the grid size on the connectivity and the indicators. Such impact disappears when comparing samples of the realizations with the same sizes. The method is then able to rank realizations from a referring model based on their static connectivity. This application also gives rise to more practical advices. The multidimensional scaling appears as a powerful visualization tool, but it also induces dissimilarity misrepresentations: it should always be interpreted cautiously with a look at the point position confidence. The heat map displays the real dissimilarities and is more appropriate for a detailed analysis. The comparison with a multiple-point histogram method shows the benefit of the connected components: the large-scale connectivity seems better characterized by our indicators, especially the skeleton indicators

OM-MADE:An open-source program to simulate one-dimensional solute transport in multiple exchanging conduits and storage zones

International audienceOM-MADE (One-dimensional Model for Multiple Advection, Dispersion, and storage in Exchanging zones) is an open-source python code for simulating one-dimensional solute transport in multiple exchanging conduits and storage zones in steady-state flow conditions. It aims at helping the interpretation of multi-peaked skewed breakthrough curves (BTCs) that can be observed in tracer tests conducted in karstic systems. OM-MADE is based on the resolution of classical mass conservation equations. In OM-MADE, all parallel and exchanging flow zones are divided along the direction of flow into reaches, in which all model parameters are kept constant. The total flowrate may be modified through lateral in and outflows. The solute may also be affected by decay processes either in mobile or immobile zones. Each reach is subdivided into discrete segments of equal length. The partial differential equations can be solved using two second order schemes, one based on an operator-split approach, the other on Crank-Nicholson pondered scheme. A verification is performed against analytical solutions, OTIS software [Runkel, 1998], and the dual-advection dispersion equation (DADE) proposed by Field and Leij [2012]. An application to a tracer test carried out in the karstic area of Furfooz (Belgium) is then performed to reproduce the double-peaked positively skewed BTC that has been observed. It constitutes a demonstration of the software capacities in the case of two reaches and three exchanging zones, among which two are mobile ones and one represents a storage zone. It thus permits to verify numerically the consistency of the conceptual interpretation of the observed BTC

Ennoyage des mines de fer lorraines : impact sur la qualité de l'eau

International audienceThe flooding of the abandoned ironmines of Lorraine leads to a degradation of the groundwater quality. Laboratory experiments allowed us to build a kinetic chemical model based on simple chemical mechanisms. During mining operations, pyrite oxidation and carbonate dissolution lead to the precipitation of gypsum. The local decrease of pH favours local dissolutions of minerals and releases ions that are fixed on cationic exchange sites. During the flooding of the mine workings, the dissolution of the newly precipitated gypsum, the precipitation of carbonates and cationic exchange reactions are responsible for the increase of the concentrations of sulphate, magnesium, sodium, potassium and strontium in water. Thereafter, these concentrations would decrease with the natural leaching of the mining reservoir.L’ennoyage des mines de fer abandonnées de Lorraine entraîne une dégradation de la qualité de l’eau souterraine. Des expériences en laboratoire ont permis la construction d’un modèle chimique cinétique, basé sur des schémas réactionnels simples. Lors de l’exploitation, l’oxydation de la pyrite et la dissolution des carbonates conduisent à la précipitation de gypse. Les baisses locales de pH favorisent des dissolutions ponctuelles de minéraux et libèrent des ions qui se fixent sur des sites d’échanges cationiques. Lors de l’ennoyage, la dissolution du gypse néoformé, la précipitation des carbonates et des réactions d’échanges cationiques entraînent une hausse des concentrations de l’eau en sulfate, magnésium, sodium, potassium et strontium qui devraient ensuite diminuer au fur et à mesure du lessivage naturel du réservoir

3D geomodelling combining implicit surfaces and Voronoi-based remeshing: A case study in the Lorraine Coal Basin (France)

International audienceIn this paper we demonstrate how recent geomodelling techniques can be combined and used to build a 3D geological model on a real case study: the former coal mine of Merlebach (France), that is targeted to be exploited for low-temperature geothermal energy production. From geological maps, cross-sections, borehole and mine exploitation data, we build a 3D model in which are identified the rocks and infrastructures having significantly different permeabilities. First, a structural model of the main geological interfaces in our area of interest (2 horizons and 13 faults) is built with classical geomodelling techniques. Then, we propose to model by surfaces the 71 irregularly stacked, very close and very thin, subvertical coal beds. To ease their construction, we use an implicit method which represents 3D surfaces as isovalues of a scalar field defined in a 3D tetrahedral grid of the area. The corresponding triangulated surfaces are remeshed with a recently proposed method based on Voronoi diagrams so that the exploited parts of the coal beds, now filled by sand, can be computed. The 3D surface-based geological model, in which infrastructures can be inserted as piecewise lines, can be volumetrically meshed. It is available for download as supplemental material, as well as a volumetric grid

3D modelling from outcrop data in a salt tectonic context: Example from the Inceyol mini-basin, Sivas Basin, Turkey

International audienceWe propose a 3D modelling strategy of the encased mini-basin of Inceyol in Sivas (Turkey). The challenge lies in the combination of sparse outcrop data and the complex interpretive geometry of geological structures that comes from salt tectonics. We succeeded in modelling the convoluted salt surface using an explicit indirect surface patch construction method followed by a manual mesh improvement. Then, we modelled the mini-basin sediments with an implicit approach. The result highlights the remarkable geometry of the convoluted salt horizon and its associated mini-basin by extending in 3D the geologist's interpretive 2D sections. This case study proves that building complex geometries is feasible with the existing tools and a good expertise in the various geomodelling techniques. The work also underlines the need for new methods to ease the modelling of such tectonic features from sparse data. We propose a 3D view of the model thanks to WebGL technology, as well as downloadable data to constitute a reference case study