Voici les éléments 1 - 7 sur 7
  • Publication
    Accès libre
    Exploring substitution random functions composed of stationary multi-Gaussian processes
    Simulation of random felds is widely used in Earth sciences for modeling and uncertainty quantifcation. The spatial features of these felds may have a strong impact on the forecasts made using these felds. For instance, in fow and transport problems the connectivity of the permeability felds is a crucial aspect. Multi-Gaussian random felds are the most common tools to analyze and model continuous felds. Their spatial correlation structure is described by a covariance or variogram model. However, these types of spatial models are unable to represent highly or poorly connected structures even if a broad range of covariance models can be employed. With this type of model, the regions with values close to the mean are always well connected whereas the regions of low or high values are isolated. Substitution random functions (SRFs) belong to another broad class of random functions that are more fexible. SRFs are constructed by composing (Z = Y◦T) two stochastic processes: the directing function T (latent feld) and the coding process Y (modifying the latent feld in a stochastic manner). In this paper, we study the properties of SRFs obtained by combining stationary multi-Gaussian random felds for both T and Y with bounded variograms. The resulting SRFs Z are stationary, but as T has a fnite variance Z is not ergodic for the mean and the covariance. This means that single realizations behave diferently from each other. We propose a simple technique to control which values (low, intermediate, or high) are connected. It consists of adding a control point on the process Y to guide every single realization. The conditioning to local values is obtained using a Gibbs sampler.
  • Publication
    Accès libre
    A parsimonious parametrization of the Direct Sampling algorithm for multiple-point statistical simulations
    Multiple-point statistics algorithms allow modeling spatial variability from training images. Among these techniques, the Direct Sampling (DS) algorithm has advanced capabilities, such as multivariate simulations, treatment of non-stationarity, multi-resolution capabilities, conditioning by inequality or connectivity data. However, finding the right trade-off between computing time and simulation quality requires tuning three main parameters, which can be complicated since simulation time and quality are affected by these parameters in a complex manner. To facilitate the parameter selection, we propose the Direct Sampling Best Candidate (DSBC) parametrization approach. It consists in setting the distance threshold to 0. The two other parameters are kept (the number of neighbors and the scan fraction) as well as all the advantages of DS. We present three test cases that prove that the DSBC approach allows to identify efficiently parameters leading to comparable or better quality and computational time than the standard DS parametrization. We conclude that the DSBC approach could be used as a default mode when using DS, and that the standard parametrization should only be used when the DSBC approach is not sufficient.
  • Publication
    Accès libre
    Analysis and stochastic simulation of geometrical properties of conduits in karstic networks
    (2020-11)
    Frantz, Yves
    ;
    Collon, Pauline
    ;
    ;
    Viseur, Sophie
    Despite intensive explorations by speleologists, karstic systems remain only partially described as many conduits are not accessible to humans. The classical exploration techniques produce sparse data, leading to various uncertainties about the conduit dimensions, essential parameters for flow simulations. Stochastic simulations offer a possibility to better assess these uncertainties. In this paper, we propose different methods to stochastically simulate the properties (size and shape anisotropy) of karstic conduits on already existing skeletons. These approaches, based on Sequential Gaussian Simulations (SGS), allow taking different conditioning data into account, while respecting the intricacy of the networks. To infer the input parameters, we perform a statistical study of the conduit dimensions of 49 explored karstic networks, focusing on their equivalent radius and width-height ratio. Thanks to the definition of 1D curvilinear variograms, we demonstrate the existence of a spatial correlation along the networks, which is even higher when considering independently each conduit. Finally, using ad hoc algorithms implemented for computing both a conduit hierarchy inside karstic networks and a relative position regarding outputs, we find no evidence of an obvious link between these two entities and the studied metrics. The simulation methods are then demonstrated on the karstic network of Arrestelia (Pyrénées-Atlantiques, France), and show the consistency of the proposed approach with the observations made on the explored natural systems.
  • Publication
    Accès libre
    A geostatistical approach to the simulation of stacked channels
    (2017)
    Rongier, G
    ;
    Collon, P
    ;
    Turbiditic channels evolve continuously in relation to erosion-deposition events. They are often gathered into complexes and display various stacking patterns. These patterns have a direct impact on the connectivity of sand-rich deposits. Being able to reproduce them in stochastic simulations is thus of significant importance. We propose a geometrical and descriptive approach to stochastically control the channel stacking patterns. This approach relies on the simulation of an initial channel using a Lindenmayer system. This system migrates proportionally to a migration factor through either a forward or a backward migration process. The migration factor is simulated using a sequential Gaussian simulation or a multiple-point simulation. Avulsions are performed using a Lindenmayer system, similarly to the initial channel simulation. This method makes it possible to control the connectivity between the channels by adjusting the geometry of the migrating areas. It furnishes encouraging results with both forward and backward migration processes, even if some aspects such as data conditioning still need to be explored.
  • Publication
    Accès libre
    Stochastic simulation of channelized sedimentary bodies using a constrained L-system
    (2017)
    Rongier, G
    ;
    Collon, P
    ;
    Simulating realistic sedimentary bodies while conditioning all the available data is a major topic of research. We present a new method to simulate the channel morphologies resulting from the deposition processes. It relies on a formal grammar system, the Lindenmayer system, or L-system. The L-system puts together channel segments based on user-defined rules and parameters. The succession of segments is then interpreted to generate non-rational uniform B-splines representing straight to meandering channels. Constraints attract or repulse the channel from the data during the channel development. They enable to condition various data types, from well data to probability cubes or a confinement. The application to a synthetic case highlights the method's ability to manage various data while preserving at best the channel morphology.
  • Publication
    Accès libre
    A methodology for pseudo-genetic stochastic modeling of discrete fracture networks
    (2013-1-10)
    Bonneau, François
    ;
    Henrion, Vincent
    ;
    Caumon, Guillaume
    ;
    ;
    Sausse, Judith
  • 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.