A Markov chain model for characterizing medium heterogeneity and sediment layering structure

Abstract

[1] By leveraging use of ''soft'' data (e.g., initial moisture content, q i ), this study applies the transition probability (TP) based Markov chain (MC) model to sediment textural classes for characterizing the medium heterogeneity and sediment layering structure. The TP/MC method is evaluated by simulating the vadose zone moisture movement at a field site, where the stratigraphy consists of imperfectly stratified soil layers. Soil heterogeneity is characterized via spatial variability of the geometry of soil textural classes. When the q i measurements, which carry signature about medium heterogeneity and stratigraphy, are not included in the TP/MC model, it is not possible to identify the horizontal TP. The q i measurements, when transformed into soil classes, are necessary in mapping the soil layering structure prevalent at the site. The soil hydraulic parameters for each soil class are treated deterministically and are estimated on the basis of core samples. To evaluate uncertainty in characterizing geometry of the soil classes, multiple conditional realizations of the soil classes are generated. A Monte Carlo simulation shows that the simulated mean moisture contents agree well with corresponding field observations. The observed splitting of the moisture plume in a coarse sand layer that is sandwiched between two fine-textured layers, the southeastward movement of the plume during the redistribution period, and the near-zero fluid flux below the bottom fine layer are adequately simulated. Spatial variability of the field-measured moisture content is sufficiently captured by the 95% confidence intervals calculated from the Monte Carlo simulations. Investigating the effect of data conditioning on the simulated results shows that a reduction of conditioning data does not necessarily deteriorate simulation results if other conditioning data exist within the mean length of the soil classes. The TP/MC method is flexible so that other types of site characterization data (e.g., geophysical data) can be incorporated as they become available.Citation: Ye, M., and R. Khaleel (2008), A Markov chain model for characterizing medium heterogeneity and sediment layering structure, Water Resour. Res., 44, W09427,

Simulation of industrial-scale CO2 storage: Multi-scale heterogeneity and its impacts on storage capacity, injectivity and leakage

“…Well log data from 14 boreholes outside RSU but in the same sedimentary basin were used to obtain statistical parameters for porosity such as stratum thickness, vertical mean length and volumetric proportion (Table 1). Spatial random fields of porosity in four geological formations above were generated by transition probability based geostatiscal simulations (Carle and Foggs, 1996;Dai et al, 2004;Ye and Khaleel, 2008). Once spatial distributions of porosity in these rocks had been acquired, permeability spatial distributions were obtained on the basis of empirical correlation between porosity and permeability.…”

Doe-Epscor Sponsored Project Final Report

“…By leveraging use of ''soft'' data (e.g., initial moisture content), the work of Ye and Khaleel (2008) applied the transition probability (TP) based Markov chain (MC) model to sediment textural classes for characterizing the medium heterogeneity and sediment layering structure. The TP/MC method was evaluated by simulating the vadose zone moisture movement at the field injection site at Hanford, where the stratigraphy consists of imperfectly stratified soil layers.…”

Challenges in conditioning a stochastic geological model of a heterogeneous glacial aquifer to a comprehensive soft data set

“…, and Ye and Khaleel (2008) use additional spatial information obtained from soil surveys, sequence stratigraphy and soil moisture, respectively, for accessing the complex lateral sedimentary variability and thus improving the quality of the model in terms of spatial variability. It has not been tested whether TProGS is capable of handling abundant soft conditioning data.…”