Voici les éléments 1 - 8 sur 8
  • Publication
    Accès libre
    Missing data simulation inside flow rate time-series using multiple-point statistics
    The direct sampling (DS) multiple-point statistical technique is proposed as a non-parametric missing data simulator for hydrological flow rate time-series. The algorithm makes use of the patterns contained inside a training data set to reproduce the complexity of the missing data. The proposed setup is tested in the reconstruction of a flow rate time-series while considering several missing data scenarios, as well as a comparative test against a time-series model of type ARMAX. The results show that DS generates more realistic simulations than ARMAX, better recovering the statistical content of the missing data. The predictive power of both techniques is much increased when a correlated flow rate time-series is used, but DS can also use incomplete auxiliary time-series, with a comparable prediction power. This makes the technique a handy simulation tool for practitioners dealing with incomplete data sets.
  • Publication
    Accès libre
    Can one identify karst conduit networks geometry and properties from hydraulic and tracer test data?
    Karst aquifers are characterized by extreme heterogeneity due to the presence of karst conduits embedded in a fractured matrix having a much lower hydraulic conductivity. The resulting contrast in the physical properties of the system implies that the system reacts very rapidly to some changes in the boundary conditions and that numerical models are extremely sensitive to small modifications in properties or positions of the conduits. Furthermore, one major issue in all those models is that the location and size of the conduits is generally unknown. For all those reasons, estimating karst network geometry and their properties by solving an inverse problem is a particularly difficult problem.
    In this paper, two numerical experiments are described. In the first one, 18,000 flow and transport simulations have been computed and used in a systematic manner to assess statistically if one can retrieve the parameters of a model (geometry and radius of the conduits, hydraulic conductivity of the conduits) from head and tracer data. When two tracer test data sets are available, the solution of the inverse problems indicate with high certainty that there are indeed two conduits and not more. The radius of the conduits are usually well identified but not the properties of the matrix. If more conduits are present in the system, but only two tracer test data sets are available, the inverse problem is still able to identify the true solution as the most probable but it also indicates that the data are insufficient to conclude with high certainty.
    In the second experiment, a more complex model (including non linear flow equations in conduits) is considered. In this example, gradient-based optimization techniques are proved to be efficient for estimating the radius of the conduits and the hydraulic conductivity of the matrix in a promising and efficient manner.
    These results suggest that, despite the numerical difficulties, inverse methods should be used to constrain numerical models of karstic systems using flow and transport data. They also suggest that a pragmatic approach for these complex systems could be to generate a large set of karst conduit network realizations using a pseudo-genetic approach such as SKS, and for each karst realization, flow and transport parameters could be optimized using a gradient-based search such as the one implemented in PEST.
  • Publication
    Accès libre
    Stochastic fracture generation accounting for the stratification orientation in a folded environment based on an implicit geological model
    This paper presents a new approach in generating stochastic discrete fracture networks. The particularity of the approach is that it allows us to simulate the theoretical families of fractures that are expected in a folded environment. The approach produces fractures that are consistent with the local stratigraphic orientation. The fractures are modeled as simple rectangular planar objects. When they are modeled, they are rotated according to the local stratigraphic orientation. As the stratigraphy is modeled using an implicit approach, we use the gradient of this geological potential field to retrieve the information about the geological orientation. The fracture number and size are following user-defined probability density functions.
  • Publication
    Accès libre
    Generation of 3D Spatially Variable Anisotropy for Groundwater Flow Simulations
    (2015) ; ;
    Courrioux, Gabriel
    Sedimentary units generally present anisotropy in their hydraulic properties, with higher hydraulic conductivity along bedding planes, rather than perpendicular to them. This common property leads to a modeling challenge if the sedimentary structure is folded. In this paper, we show that the gradient of the geological potential used by implicit geological modeling techniques can be used to compute full hydraulic conductivity tensors varying in space according to the geological orientation. For that purpose, the gradient of the potential, a vector normal to the bedding, is used to construct a rotation matrix that allows the estimation of the 3D hydraulic conductivity tensor in a single matrix operation. A synthetic 2D cross section example is used to illustrate the method and show that flow simulations performed in such a folded environment are highly influenced by this rotating anisotropy. When using the proposed method, the streamlines follow very closely the folded formation. This is not the case with an isotropic model.
  • Publication
    Accès libre
    3D stochastic modeling of karst aquifers using a pseudo-genetic methodology
    Le but de cette thèse est le développement d'une méthodologie de modélisation des aquifères karstiques. Premièrement, la géométrie des conduits karstiques est simulée. La géométrie de ces conduits est contrôlée à large échelle par la géologie (modèle géologique), et à plus petite échelle par la fracturation (modèle stochastique de fracturation).
    Deuxièmement, ces modèles géométriques (formations 3D et conduits ensemble) sont utilisés comme base pour la simulation d'écoulement et de transport.
    En dernière partie, le simulateur de conduits SKS ("Stochastic Karst Simulator") est couplé avec la simulation d'écoulement et transport pour investiguer une approche inverse qui permette d'utiliser cette méthodologie dans une étude d'incertitude.
    La méthodologie de simulation des conduits karstiques développée ici se base sur une approche dite "pseudo-génétique", càd qui mime les résultats des processus de spéléogénèses, sans pour autant simuler toute la dynamique complexe de ces processus, comme la dissolution et le transport réactif de calcite, etc. Dans cette approche pseudo-génétique les conduits karstiques sont simulés par une physique approchée, qui se base sur le principe de minimisation de l'énérgie. L'eau se déplace dans un milieu en cherchant le chemin de moindre résistance. Ce principe est utilisé ici, par l'utilisation d'un algorithme de Fast Marching, qui permet de calculer le chemin de moindre effort., The focus of this thesis is the development of a methodology to model karst aquifers. First the geometry of the karst conduits is simulated. Their geometry is controlled by the geology at large scale (geological model) and by the fracturation at smaller scale (stochastic model of fractures).
    Secondly, these geometrical models (3D geological formation together with conduits) are used as base for flow and transport simulation.
    Finally, the karst conduit generator called SKS ("Stochastic Karst Simulator") is coupled with the physical simulation (flow and transport) to investigate an inverse approach which allows to use this methodology in an uncertainty analysis.
    The karst conduits simulation methodology is called pseudo genetic because it mimic the results of the speleogenetic processes, without simulating all the complex kinetic of karst systems genesis, like reactive transport, calcite dissolution and precipitation. In this approach, the karst conduits are simulated by approaching the physical principle of minimization of energy using a Fast Marching Algorithm. This algorithm allow to compute the minimum effort path, which is assumed to be the one used by water, and consequently the preferential dissolution., Das Ziel dieser Dissertation besteht in der Entwicklung einer Methode zur Modellierung von Karstaquiferen. In einem ersten Schritt wird die Geometrie des Karstsystems simuliert. Im grossen Massstab wird die Geometrie des Karstsystems durch die Geologie bestimmt (geologisches Modell), in kleinerem Massstab durch Klüfte (stochastisches Modell zur Kluftgenese). In einem zweiten Schritt wird dieses Modell als Grundlage für die Modellierung von Grundwasserströmung und Stoffransport im Karstsystem verwendet.
    Schliesslich wird die Simulation des Karstsystems mittels des Karst-Simulator (SKS, "Stochastic Karst Simulator") mit Simulationen von Grundwasserströmung und Stoffransport gekoppelt, um einen inversen Ansatz zu untersuchen, der es erlauben würde diese Methode im Rahmen einer Unsicherheitsanalyse zu verwenden.
    Die hier entwickelte Methodik zur Simulierung von Karstaquiferen basiert auf einem sogenannt "pseudogenetischen" Ansatz, weil sie speleogenetische Prozesse nachbildet, ohne die komplexe Dynamik dieser Prozesse, wie Auflösung und reaktiver Transport von Kalzit etc., im Detail zu simulieren. Die Karstgenese wird in dieser Methodik vielmehr durch einen Ansatz angenähert, der auf dem Prinzip der Energieminimierung beruht. Dabei wird ein Fast Marching Algorithmus verwendet, um zwischen den Stellen mit Wasserzuflüssen ins Karstsystem und den Karstquellen den Weg des geringsten Widerstandes zu berechnen., Lo scopo di questa tesi, è lo svipuppo di una metodologia di modellizzazione realistica degli acquiferi carsici. Nella prima parte viene modellizzata la geometria dei condotti carsici. La geometria di questi ultimi, è controllata a larga scala dalla geologia (modello geologico) e a piu piccola scala dalla fratturazione (modello di fratturazione stocastico).
    In seguito questi modelli geometrici vengono usati come base per la simulazione di fluidi e trasporto di contaminanti. E infine, il simulatore di condotti carsici SKS ("Stochastic Karst Simulator") è usato insieme alla simulazione di fluidi e trasporto per valutare l'applicabilità di un approcio inverso che permetta di utilizzare questa metodologia di simulazione.
    La metodologia di simulazione dei condotti carsici è detta "pseudo genetica", perchè tenta di approssimare la fisica complessa della speleogenesi (come il trasporto reattive, la dissoluzione/precipitazione della calcite) senza dover risolverla numericamente. In effetti si basa sul principio fisico della minimizzazione dell'energia, usando un algoritmo di Fast Marching per calcolare il cammino di minor resistenza (cioè quello che dovrebbe seguire l'acqua). In pratica questa metodologia simula dei sistemi carsici maturi, direttamente nel loro stato finale.
  • Publication
    Accès libre
    A pseudo-genetic stochastic model to generate karstic networks
    (2012) ; ;
    Jenni, Sandra
    In this paper, we present a methodology for the stochastic simulation of 3D karstic conduits accounting for conceptual knowledge about the speleogenesis processes and accounting for a wide variety of field measurements.
    The methodology consists of four main steps. First, a 3D geological model of the region is built. The second step consists in the stochastic modeling of the internal heterogeneity of the karst formations (e.g. initial fracturation, bedding planes, inception horizons, etc.). Then a study of the regional hydrology/hydrogeology is conducted to identify the potential inlets and outlets of the system, the base levels and the possibility of having different phases of karstification. The last step consists in generating the conduits in an iterative manner using a fast marching algorithm. In most of these steps, a probabilistic model can be used to represent the degree of knowledge available and the remaining uncertainty depending on the data at hand.
    The conduits are assumed to follow minimum effort paths in a heterogeneous medium from sinkholes or dolines toward springs. The search of the shortest path is performed using a fast marching algorithm. This process can be iterative, allowing to account for the presence of already simulated conduits and to produce a hierarchical network.
    The final result is a stochastic ensemble of 3D karst reservoir models that are all constrained by the regional geology, the local heterogeneities and the regional flow conditions. These networks can then be used to simulate flow and transport. Several levels of uncertainty can be considered (large scale geological structures, local heterogeneity, position of possible inlets and outlets, phases of karstification).
    Compared to other techniques, this method is fast, to account for the main factors controlling the 3D geometry of the network, and to allow conditioning from available field observations.
  • Publication
    Accès libre
    A method for the stochastic modeling of karstic systems accounting for geophysical data:: an example of application in the region of Tulum, Yucatan Peninsula (Mexico)
    Vuilleumier, Cécile
    ;
    ; ;
    Ottowitz, David
    ;
    Schiller, A.
    ;
    Supper, Robert
    ;
    The eastern coast of the Yucatan Peninsula, Mexico, contains one of the most developed karst systems in the world. This natural wonder is undergoing increasing pollution threat due to rapid economic development in the region of Tulum, together with a lack of wastewater treatment facilities. A preliminary numerical model has been developed to assess the vulnerability of the resource. Maps of explored caves have been completed using data from two airborne geophysical campaigns. These electromagnetic measurements allow for the mapping of unexplored karstic conduits. The completion of the network map is achieved through a stochastic pseudo-genetic karst simulator, previously developed but adapted as part of this study to account for the geophysical data. Together with the cave mapping by speleologists, the simulated networks are integrated into the finite-element flow-model mesh as pipe networks where turbulent flow is modeled. The calibration of the karstic network parameters (density, radius of the conduits) is conducted through a comparison with measured piezometric levels. Although the proposed model shows great uncertainty, it reproduces realistically the heterogeneous flow of the aquifer. Simulated velocities in conduits are greater than 1 cm s−1, suggesting that the reinjection of Tulum wastewater constitutes a pollution risk for the nearby ecosystems.