济南泉域岩溶裂隙介质的多重指示克里格法研究

本文选择研究程度较高,具有代表性的山东济南泉域的含水岩溶裂隙介质,运用多重指示克里格方法研究岩溶裂隙介质的空隙连通性能,并预测岩溶裂隙密集带和强径流带的位置。研究表明,多重指示克里格方法能够通过不同边界域值的设定,有效地克服由于非平稳区域化变量产生的参数"漂移",取得较理想的模拟预测效果。     

A method for the interpolation of nonnegative functions with an application to contaminant load estimation

The objective of this work is to extend kriging, a geostatistical interpolation method, to honor parameter nonnegativity. The new method uses a prior probability distribution based on reflected Brownian motion that enforces this constraint. The work presented in this paper focuses on interpolation problems where the unknown is a function of a single variable (e.g. time), and is developed both for the case with and without measurement error in the available data. The algorithms presented for conditional simulations are computationally efficient, particularly in the case with no measurement error. We present an application to the interpolation of dissolved arsenic concentration data from the North Fork of the Humboldt River, Nevada.     

Geostatistical interpolation of object counts collected from multiple strip transects: Ordinary kriging versus finite domain kriging

Data collected along transects are becoming more common in environmental studies as indirect measurement devices, such as geophysical sensors, that can be attached to mobile platforms become more prevalent. Because exhaustive sampling is not always possible under constraints of time and costs, geostatistical interpolation techniques are used to estimate unknown values at unsampled locations from transect data. It is known that outlying observations can receive significantly greater ordinary kriging weights than centrally located observations when the data are contiguously aligned along a transect within a finite search window. Deutsch (1994) proposed a kriging algorithm, finite domain kriging, that uses a redundancy measure in place of the covariance function in the data-to-data kriging matrix to address the problem of overweighting the outlying observations. This paper compares the performances of two kriging techniques, ordinary kriging (OK) and finite domain kriging (FDK), on examining unexploded ordnance (UXO) densities by comparing prediction errors at unsampled locations. The impact of sampling design on object count prediction is also investigated using data collected from transects and at random locations. The Poisson process is used to model the spatial distribution of UXO for three 5000 × 5000 m fields; one of which does not have any ordnance target (homogeneous field), while the other two sites have an ordnance target in the center of the site (isotropic and anisotropic fields). In general, for a given sampling transects width, the differences between OK and FDK in terms of the mean error and the mean square error are not significant regardless of the sampled area and the choice of the field. When 20% or more of the site is sampled, the estimation of object counts is unbiased on average for all three fields regardless of the choice of the transect width and the choice of the kriging algorithm. However, for non-homogeneous fields (isotropic and anisotropic fields), the mean error fluctuates considerably when a small number of transects are sampled. The difference between the transect sampling and the random sampling in terms of prediction errors becomes almost negligible if more than 20% of the site is sampled. Overall, FDK is no better than OK in terms of the prediction performances when the transect sampling procedure is used.     

Estimation and prediction of a 2D lognormal diffusion random field

This paper describes techniques for estimation and prediction of two-parameter lognormal diffusion random fields. The drift and diffusion coefficients, which characterize a two-parameter lognormal diffusion under certain conditions, are estimated by maximum likelihood. For data on a regular grid, an alternative method is proposed to estimate the diffusion coefficient. Both of these estimates are compared in several situations. The kriging predictors are formulated involving the drift and diffusion coefficients and the predictions are obtained using the estimates of these coefficients.     

A Geostatistical Simulation Approach to a Pollution Case in Northeastern Italy

The study zone is located in the city of Padova (the Veneto region, NE Italy). The industrial activities present in this area since 1950 have produced very high concentrations of Pb, Zn, Cu, Pcb and oils until a depth of 7 m. The unsaturated and the saturated zones are both polluted. Because of the abundance of Pb values it was decided to analyze the lead distribution in the studied zone. In many studies of the polluted sites, the geometry of the investigated volume is highly anisotropic. Generally we have an extension of some hectares in the horizontal plane and of a few meters in depth. It is likely that different horizontal spatial continuity structures in pollution distribution are found at different depths both for the layered nature of the medium and for the transition between unsaturated and saturated conditions. In such conditions the decision to divide a 3D problem into 1D and 2D problems can be useful. Initially the studied volume was divided into seven layers up to 5 m in depth; the study was then approached in two phases. First, the Pb values in the vertical direction were analyzed, considering a derive along z, and estimating the values using the Kriging with Trend (KT) method. Thus it was possible to increase the data in the z direction, especially in the deeper layers. Second, 500 realizations of the Pb distribution for each of the seven layers were simulated using the simulated annealing procedure. Finally, results were presented and discussed for each layer in terms of median and probability maps.     

The Contemporary Tectonic Strain Rate Field of Continental China Predicted from GPS Measurements and its Geodynamic Implications

In Appling studies, the strain rate in China has been computed from using different methods, resulting in quite different estimates of the strain rate the geo-statistics from GPS velocity field of Chinese continent, we obtain the velocity value at each little regularly spaced grid point, by kriging interpolation and the component of strain rate for each volume element, using a method similar to the derivation of shape functions in the finite element algorithm. Therefore the distribution of the strain rate field in whole for the Chinese continent is presented. The result shows that the orientations of principal strain rates are consistent with those of the P and T axes of focal mechanisms. The distribution of maximum shear strain rate clearly delineates some major active fault zones surrounding the Tibetan Plateau. The maximum shear strain rate is comparable with that obtained from analysis of seismic moment release. In part of the Tibetan Plateau containing normal faults and pull-apart grabens, we obtain an extensional state of strain. The absolute value of the strain rate in West China is approximately 5 times larger than that of East China, and the pattern of the strain rate field in most of the Chinese continent is controlled by the India/Eurasia collision.     

基于地形随机场模型的遥感图像几何纠正

从地形平稳随机场模型出发,应用空间统计学中的Kriging插值算法,提出在多项式纠正基础上进行误差改正的图像精纠正方法,并进行相应的实验验证。多项式纠正和基于共线方程的几何成像模型纠正是最常用的两类纠正方法,但多项式纠正无法克服由于地形起伏所引起的投影差。提出的方法和流程,不仅能克服多项式模型中的弊端,而且不额外增加控制点选择和计算的负担,可以作为多项式纠正流程中的求精。     

北京市近郊区土壤镍的空间结构及分布特征

通过对北京市近郊区349处土壤的取样调查，利用地统计学方法研究了北京市土壤中镍的空问结构和分布特征。结果显示，镍的空间结构性较好，具有明显的可迁性特点，块金效应所占的比例较大；实验变异函数可用指数模型拟合，利用所得模型进行克立格插值得到土壤镍含量的分布图，发现北京市土壤中，镍的含量分布具有明显的西北低而东南高的趋势。与北京市土壤镍含量的背景值进行对比，全市土壤并不存在非常严重的镍污染，其空间分布主要受成土母质制约；但是人类活动、风向、河流等人文和自然条件对镍的空间结构和分布格局也有一定程度的影响。     

Précision de l'estimation des précipitations au Sahel selon la densité du réseau d'observation pluviométrique

Many earlier studies have shown the very large spatial variability of rainfall in the Sahel at all time steps, from the event to the season. Often, the meteorological network in these countries is sparse, with one to five rain gauges per 10 000 km². It is thus difficult to calculate accurate estimates of the mean rainfall over such a large area. To improve the knowledge of Sahelian systems and the spatial distribution of rainfall, a dense network was set up in an area of 16 000 km² in southwestern Niger between 1991 and 1996. The aim was to calculate accurate rainfall spatial means over an area of 12 000 km² at different time steps (from the season to the ten-day period). With the spatialisation method used (kriging), it was possible to calculate curves of estimation errors of mean rainfall versus the rain-gauge network density. Operational abacuses of the standard estimation error as a function of the spatial mean of rainfall and the network density are proposed.     

合肥市大兴地区土壤中镉元素的空间分布的多维分形特征研究

土壤中cd元素空间分布的研究是土壤科学的重要内容,也是土壤环境污染评价与防治的重要依据。空间变异研究中最为常用的克里格插值法能够较为精确的进行空间估值,但由于它的低通滤波性,使得对于较高数值的异常点往往不能正确估计。而异常区域的元素分布状况恰是环境污染评价与防治最为关注的。因此,本文采用多维分形地统计学方法,对大兴地区的镉元素进行了空间插值分析。结果表明,该方法能较好的反映cd元素的空间分布奇异性。