On the estimation of scale of fluctuation in geostatistics

Describing how soil properties vary spatially is of particular importance in stochastic analyses of geotechnical problems, because spatial variability has a significant influence on local material and global geotechnical response. In particular, the scale of fluctuation θ is a key parameter in the correlation model used to represent the spatial variability of a site through a random field. It is, therefore, of fundamental importance to accurately estimate θ in order to best model the actual soil heterogeneity. In this paper, two methodologies are investigated to assess their abilities to estimate the vertical and horizontal scales of fluctuation of a particular site using in situ cone penetration test (CPT) data. The first method belongs to the family of more traditional approaches, which are based on best fitting a theoretical correlation model to available CPT data. The second method involves a new strategy which combines information from conditional random fields with the traditional approach. Both methods are applied to a case study involving the estimation of θ at three two-dimensional sections across a site and the results obtained show general agreement between the two methods, suggesting a similar level of accuracy between the new and traditional approaches. However, in order to further assess the relative accuracy of estimates provided by each method, a second numerical analysis is proposed. The results confirm the general consistency observed in the case study calculations, particularly in the vertical direction where a large amount of data are available. Interestingly, for the horizontal direction, where data are typically scarce, some additional improvement in terms of relative error is obtained with the new approach.     

日照茶树种植区表层土壤营养元素地球化学特征分析

对土壤中的营养元素进行分析,对于茶树种植及宜种性评价意义深远。该研究主要是利用1 km×1 km表层土壤采样数据,借助地统计学和GIS技术对日照茶树种植区表层土壤中主要营养元素的地球化学特征进行分析。结果表明,研究区N,P,K2O,S含量相对较丰富,而其他营养元素相对较贫乏;营养元素富集区大致以西湖、秦楼、岚山头3个镇为顶点呈三角分布。     

条件模拟原理和技术进展及在水资源系统中的应用

在简要介绍条件模拟基本原理基础上,重点评述了国外20世纪90年代新的数学方法和计算技术进展及在国内外水资源系统应用的主要成果,并对今后学科和技术发展作了展望.     

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.     

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.     

Efficiency of kriging estimation for square, triangular, and hexagonal grids

Although several researchers have pointed out some advantages and disadvantages of various soil sampling designs in the presence of spatial autocorrelation, a more detailed study is presented herein which examines the geometrical relationship of three sampling designs, namely the square, the equilateral triangle, and the regular hexagon. Both advantages and disadvantages exist in the use of these designs with respect to estimation of the semivariogram and their effect on the mean square error or variance of error. This research could be used to design optimal sampling strategies; it is based on the theory of regionalized variables, in which the intrinsic hypothesis is satisfied. Among alternative designs, an equilateral triangle design gives the most reliable estimate of the semivariogram. It also gives the minimum maximum mean square error of point estimation of the concentration over the other two designs for the same number of measurements when the nugget effect is small relative to the variance. If the nugget effect is large (.90 ² or more), and the linear sampling density is >0.85r where r is the range, the hexagonal design is best. This study computes and compares the maximum mean square error for each of these designs.     

Stepwise Conditional Transformation for Simulation of Multiple Variables

Most geostatistical studies consider multiple-related variables. These relationships often show complex features such as nonlinearity, heteroscedasticity, and mineralogical or other constraints. These features are not handled by the well-established Gaussian simulation techniques. Earth science variables are rarely Gaussian. Transformation or anamorphosis techniques make each variable univariate Gaussian, but do not enforce bivariate or higher order Gaussianity. The stepwise conditional transformation technique is proposed to transform multiple variables to be univariate Gaussian and multivariate Gaussian with no cross correlation. This makes it remarkably easy to simulate multiple variables with arbitrarily complex relationships: (1) transform the multiple variables, (2) perform independent Gaussian simulation on the transformed variables, and (3) back transform to the original variables. The back transformation enforces reproduction of the original complex features. The methodology and underlying assumptions are explained. Several petroleum and mining examples are used to show features of the transformation and implementation details.     

Efficient generation of conditional simulations by chebyshev matrix polynomial approximations to the symmetric square root of the covariance matrix

Consider the problem of generating a realization y₁ of a Gaussian random field on a dense grid of points ₁ conditioned on field observations y₂ collected on a sparse grid of points ₂. An approach to this is to generate first an unconditional realization y over the grid =_{1 2}, and then to produce y₁ by conditioning y on the data y₂. As standard methods for generating y, such as the turning bands, spectral or Cholesky approaches can have various limitations, it has been proposed by M. W. Davis to generate realizations from a matrix polynomial approximations to the square root of the covariance matrix. In this paper we describe how to generate a direct approximation to the conditional realization y₁, on ₁ using a variant of Davis' approach based on approximation by Chebyshev polynomials. The resulting algorithm is simple to implement, numerically stable, and bounds on the approximation error are readily available. Furthermore we show that the conditional realization y₁ can be generated directly with a lower order polynomial than the unconditional realization y, and that further reductions can be achieved by exploiting a nugget effect if one is present. A pseudocode version of the algorithm is provided that can be implemented using the fast Fourier transform if the field is stationary and the grid ₁ is rectangular. Finally, numerical illustrations are given of the algorithm's performance in generating various 2-D realizations of conditional processes on large sampling grids.     

Evaluation and comparison of spatial interpolators II

The performance of several variations on ordinary kriging and inverse distance estimators is evaluated. Mean squared errors (MSE) were calculated for estimates made on multiple resamplings from five exhaustive data bases representing two distinctly different types of estimation problem. Ordinary kriging, when performed with variograms estimated from the sample data, was more robust than inverse-distance methods to the type of estimation problem, and to the choice of estimation parameters such as number of neighbors.Notice: Although the research described in this article has been funded in part by the United States Environmental Protection Agency through Cooperative Agreement CR818526 to the Harry Reid Center for Environmental Studies, University of Nevada-Las Vegas, it has not been subjected to Agency review. Therefore it does not necessarily reflect the views of the Agency. Mention of trade names or commercial products does not constitute endorsement or recommendation for use.     

含水层渗透系数K的空间变异性研究

本文结合桥梁实例对公路双曲拱桥出现的一些主要病害进行了初步分析，并进而对桥梁状况提出评定意见和处理建议，仅供参考。