Disjunctive Kriging with Hard and Imprecise Data

This paper presents a methodology for assessing local probability distributions by disjunctive kriging when the available data set contains some imprecise measurements, like noisy or soft information or interval constraints. The basic idea consists in replacing the set of imprecise data by a set of pseudohard data simulated from their posterior distribution; an iterative algorithm based on the Gibbs sampler is proposed to achieve such a simulation step. The whole procedure is repeated many times and the final result is the average of the disjunctive kriging estimates computed from each simulated data set. Being data-independent, the kriging weights need to be calculated only once, which enables fast computing. The simulation procedure requires encoding each datum as a pre-posterior distribution and assuming a Markov property to allow the updating of pre-posterior distributions into posterior ones. Although it suffers some imperfections, disjunctive kriging turns out to be a much more flexible approach than conditional expectation, because of the vast class of models that allows its computation, namely isofactorial models.     

GPS载波相位三差观测方程的改化

本文根据相关测量平差原理推导出改化ＧＰＳ载波相位三差观测方程的算法，使改化后的协因数呈分块对角体。这种改化算法不但顾及观测值的相关性，而且节省大量的存储单元。使采用严密平差法平差三差观测值成为一种最简单的数据处理方法之一。     

Dynamic stochastic estimation of physical variables

A fundamental problem facing the physical sciences today is analysis of natural variations and mapping of spatiotemporal processes. Detailed maps describing the space/time distribution of groundwater contaminants, atmospheric pollutant deposition processes, rainfall intensity variables, external intermittency functions, etc. are tools whose importance in practical applications cannot be overestimated. Such maps are valuable inputs for numerous applications including, for example, solute transport, storm modeling, turbulent-nonturbulent flow characterization, weather prediction, and human exposure to hazardous substances. The approach considered here uses the spatiotemporal random field theory to study natural space/time variations and derive dynamic stochastic estimates of physical variables. The random field model is constructed in a space/time continuum that explicitly involves both spatial and temporal aspects and provides a rigorous representation of spatiotemporal variabilities and uncertainties. This has considerable advantages as regards analytical investigations of natural processes. The model is used to study natural space/time variations of springwater calcium ion data from the Dyle River catchment area, Belgium. This dataset is characterized by a spatially nonhomogeneous and temporally nonstationary variability that is quantified by random field parameters, such as orders of space/time continuity and random field increments. A rich class of covariance models is determined from the properties of the random field increments. The analysis leads to maps of continuity orders and covariances reflecting space/time calcium ion correlations and trends. Calcium ion estimates and the associated statistical errors are calculated at unmeasured locations/instants over the Dyle region using a space/time kriging algorithm. In practice, the interpretation of the results of the dynamic stochastic analysis should take into consideration the scale effects.     

Geostatistical regionalization of glacial aquitard thickness in northwestern Germany, based on fuzzy kriging

In many instances hydrogeological parameters obtained by conventional methods for selected localities within an aquifer or an aquitard are not sufficient for adequate regionalization at the scale of the entire layer. Here, we demonstrate an application of the fuzzy kriging method in regionalization of hydrogeological data, in which the set of conventional, crisp values is supplemented by imprecise information subjectively estimated by an expert. It is believed that such an approach eventually may reflect the real-world conditions more closely than a traditional crisp-value approach, because the former does not impose exactness artificially on phenomena which are diffuse by their nature. Spatial interpolation was done for the thickness of one of the major aquitards (till and glaciolacustrine clay) in northwestern Germany. The dataset consists of 329 crisp values from boreholes supplemented by 172 imprecise values defined as fuzzy numbers. It is demonstrated that the reliability of regionalization was higher, compared to regionalization performed with the crisp dataset only. Fuzzy kriging was performed with FUZZEKS (Fuzzy Evaluation and Kriging System) developed at the Ecosystem Research Center at the University of Kiel.     

Foliation fields and 3D cartography in geology: Principles of a method based on potential interpolation

A modeling method that takes into account known points on a geological interface and plane orientation data such as stratification or foliation planes is described and tested. The orientations data do not necessarily belong to one of the interfaces but are assumed to sample the main anisotropy of a geological formation as in current geological situations. The problem is to determine the surfaces which pass through the known points on interfaces and which are compatible with the orientation data. The method is based on the interpolation of a scalar field defined in the space the gradient in which is orthogonal to the orientations, given that some points have the same but unknown scalar value (points of the same interface), and that scalar gradient is known on the other points (foliations). The modeled interfaces are represented as isovalues of the interpolated field. Preliminary two-dimensional tests carried-out with different covariance models demonstrate the validity of the method, which is easily transposable in three dimensions.     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

中国地质统计学（空间信息统计学）发展的回顾与前景

讨论了地质统计学（空间信息统计学）发展的若干问题，其内容包括４部分：１．地质统计学的理论及方法；２．关于软件系统的现状；３．关于实际应用；４．地质统计学发展前景的若干看法。     

普通克里格法在矿产储量计算中的应用

普通克里格法能充分利用样品信息，有效地提高块段平均品位及矿产储量精度。以刁泉银铜矿为实例，介绍了该方法的理论和方法应用以及储量计算流程图。     

Data assimilation (4D-VAR) to forecast flood in shallow-waters with sediment erosion

In this paper, the four-dimensional variational data assimilation technique (4D-VAR) is presented as a tool to forecast floods. Our study is limited to purely hydrological flows and supposes that the weather, here a big rain, has been already forecasted by meteorological services. The technique consists in minimizing, in the sense of Lagrange, the cost function: a measure of the difference between calculated data and available observations, here the water level. This is done under constraints that are the equations of the physical model. In our case, we modified the shallow-water equations to include a simplified sediment transport model. The steepest descent algorithm is then used to find the minimum. This is made possible because we can compute analytically the gradient of the cost function by using the adjoint equations of the model. As an application of the 4D-VAR technique, the overflowing of the Chicoutimi River at the Chute-Garneau dam, during the 1996 flood, is investigated. It is found that the 4D-VAR method reduces the error in the water height forecast even when the erosion model is not activated. In terms of Lyapunov exponents, we estimate the predictability horizon of such an event to be about half-an-hour after a big rain. However, this limit of predictability can be increased by using more observations or by using a finer computational grid.