Transiogram

Abstract: Traditionally, spatial continuity models for indicator variables are developed by empirical curvefitting to the sample indicator (cross-) variogram. However, geologic data may be too sparse to permit a purely empirical approach, particularly in application to the subsurface. Techniques for model synthesis that integrate hard data and conceptual models therefore are needed. Interpretability is crucial. Compared with the indicator (cross-) variogram or indicator (cross-) covariance, the transition probability is more interpretable. Information on proportion, mean length, and juxtapositioning directly relates to the transition probability: asymmetry can be considered. Furthermore, the transition probability elucidates order relation conditions and readily formulates the indicator (co)kriging equations.

Abstract: The continuous-lag Markov chain provides a conceptually simple, mathematically compact, and theoretically powerful model of spatial variability for categorical variables. Markov chains have a long-standing record of applicability to one-dimensional (1-D) geologic data, but 2- and 3-D applications are rare. Theoretically, a multidimensional Markov chain may assume that 1-D Markov chains characterize spatial variability in all directions. Given that a 1-D continuous Markov chain can be described concisely by a transition rate matrix, this paper develops 3-D continuous-lag Markov chain models by interpolating transition rate matrices established for three principal directions, say strike, dip, and vertical. The transition rate matrix for each principal direction can be developed directly from data or indirectly by conceptual approaches. Application of Sylvester's theorem facilitates establishment of the transition rate matrix, as well as calculation of transition probabilities. The resulting 3-D continuous-lag Markov chain models then can be applied to geo-statistical estimation and simulation techniques, such as indicator cokriging, disjunctive kriging, sequential indicator simulation, and simulated annealing.

Abstract: Multi-dimensional Markov chain conditional simulation (or interpolation) models have potential for predicting and simulating categorical variables more accurately from sample data because they can incorporate interclass relationships. This paper introduces a Markov chain random field (MCRF) theory for building one to multi-dimensional Markov chain models for conditional simulation (or interpolation). A MCRF is defined as a single spatial Markov chain that moves (or jumps) in a space, with its conditional probability distribution at each location entirely depending on its nearest known neighbors in different directions. A general solution for conditional probability distribution of a random variable in a MCRF is derived explicitly based on the Bayes’ theorem and conditional independence assumption. One to multi-dimensional Markov chain models for prediction and conditional simulation of categorical variables can be drawn from the general solution and MCRF-based multi-dimensional Markov chain models are nonlinear.

Abstract: Markov chains, applied to stratigraphic sections, are suitable for the description of geological cycles. Cyclicity is a property which ranges from completely damped to perfect periodic oscillations. A measurement of “cyclicity” is provided by the eigenvalues of the Markov matrix. At the same time, it was found that the method of Markov chains and discrete stochastic chains, in general, is a cumber-some tool when applied to stratigraphic analysis. Many sections require high-order chains and consequently larger transition probability matrices, which are difficult to interpret. Great care must be taken not to oversimplify these chains, because this results in a loss of cyclicity. By applying the methods to a composite section of Upper Pennsylvanian strata of northeastern Kansas, a 40- to 50-ft “cycle” was found. The oscillations are damped if directly measured data are used. Removal of a trend in sedimentation rates reduces damping to such an extent that the hypothesis of time cyclicity is at least feasible.