地球物理勘探中几种二维插值方法的误差分析

为了研究采样和网格化方法对地球物理数据成图精度的影响,为野外数据采集布设提供一定的依据,采用数值模拟确定重力异常场场值,通过不同采样间距和不同插值方法计算重力异常绝对误差均方根值和节点处的绝对误差值,对比不同插值方法的误差,得到了如下认识：1）对于同一插值方法而言,存在小间距绝对误差均方根值小于大间距绝对误差均方根值的关系。2）对不同的插值方法而言:当采样间距小于最小异常地质体尺度时,绝对误差均方根值由小到大的顺序是径向基函数法、改进的谢别德法、克里金插值法、自然邻点法、反距离加权插值法、最近邻点法、最小曲率法,并且线性插值三角网法与自然邻点法具有几乎相同的数值;当采样间距大于最小异常地质体尺度时,绝对误差均方根值由小到大的顺序是径向基函数法、改进的谢别德法、克里金插值法、自然邻点法、最小曲率法、最近邻点法、反距离加权插值法,并且线性插值三角网法和自然邻点法具有几乎相同的数值。3）从绝对误差均方值看,径向基函数方法、改进的谢别德方法和克里金方法数值较小,其中径向基函数值绝对误差均方根值最小。4）从节点处绝对误差值来看,径向基函数方法、克里金方法、改进的谢别德方法相对其他插值方法具有更小的误差,不存在局部误差较小或较大的情况,是相对较好的插值方法,并且径向基函数方法是最好的。     

基于地形分区IDW的复垦排土场土壤养分空间分布

准确掌握土壤养分的空间分布，对于评价复垦排土场土壤质量，并制定相应的复垦措施具有重要意义，但是排土场的复杂地形导致一般插值方法的预测效果较差。以内蒙古锡林浩特盟胜利一号露天矿北排土场为例，采用基于地形分区的反距离权重法对土壤养分（速效氮、速效钾、有效磷和有机质）进行空间插值并与常规反距离权重法和克里金插值法进行比较。结果表明，使用分区反距离权重法插值得到的交叉验证精度均高于克里金插值法与反距离权重法，除此之外，有效磷的R2值为0.22，速效氮的R2值为0.37，速效钾以及有机质的R2值分别为0.34与0.14，同样高于克里金插值法与反距离权重法。插值结果显示，排土场速效氮含量整体由北至南逐渐减少；有效磷含量仅在中心地区稍高；速效钾在研究区东北部含量偏高；土壤有机质在中心地区含量相对偏低。研究成果为露天矿排土场土壤养分评价提供借鉴。      

Geostatistical Mapping with Continuous Moving Neighborhood

An issue that often arises in such GIS applications as digital elevation modeling (DEM) is how to create a continuous surface using a limited number of point observations. In hydrological applications, such as estimating drainage areas, direction of water flow is easier to detect from a smooth DEM than from a grid created using standard interpolation programs. Another reason for continuous mapping is esthetic; like a picture, a map should be visually appealing, and for some GIS users this is more important than map accuracy. There are many methods for local smoothing. Spline algorithms are usually used to create a continuous map, because they minimize curvature of the surface. Geostatistical models are commonly used approaches to spatial prediction and mapping in many scientific disciplines, but classical kriging models produce noncontinuous surfaces when local neighborhood is used. This motivated us to develop a continuous version of kriging. We propose a modification of kriging that produces continuous prediction and prediction standard error surfaces. The idea is to modify kriging systems so that data outside a specified distance from the prediction location have zero weights. We discuss simple kriging and conditional geostatistical simulation, models that essentially use information about mean value or trend surface. We also discuss how to modify ordinary and universal kriging models to produce continuous predictions, and limitations using the proposed models.     

An Experimental Comparison of Ordinary and Universal Kriging and Inverse Distance Weighting

A factorial, computational experiment was conducted to compare the spatial interpolation accuracy of ordinary and universal kriging and two types of inverse squared-distance weighting. The experiment considered, in addition to these four interpolation methods, the effects of four data and sampling characteristics: surface type, sampling pattern, noise level, and strength of small-scale spatial correlation. Interpolation accuracy was measured by the natural logarithm of the mean squared interpolation error. Main effects of all five factors, all two-factor interactions, and several three-factor interactions were highly statistically significant. Among numerous findings, the most striking was that the two kriging methods were substantially superior to the inverse distance weighting methods over all levels of surface type, sampling pattern, noise, and correlation.     

Accounting for exhaustive secondary data into the mapping of water table elevation

This paper presents the incorporation of a digital elevation model into the spatial prediction of water table elevation in Mazandaran province (Iran) using a range of interpolation techniques. The multivariate methods used are: linear regression (LR), cokriging (COK), kriging with an external drift (KED) and regression kriging (RK). The analysis is performed on 3 years (1987, 1997 and 2007) of water table elevation data from about 260 monitoring wells. Prediction performances of the different algorithms are compared with two univariate techniques, i.e. inverse distance weighting and ordinary kriging (OK), through cross validation and examination of the consistency of the generated maps with the natural phenomena. Significantly smaller prediction errors are obtained for four multivariate algorithms but, in particular, KED and RK outperform LR and COK for 3 years. The results show the potential for using elevation for a more precise mapping of water table elevation.     

Matheronian geostatistics—Quo vadis?

Components of geostatistical estimation, developed as a method for ore deposit assessment, are discussed in detail. The assumption that spatial observations can be treated as a stochastic process is judged to be an inappropriate model for natural data. Problems of semivariogram formulation are reviewed, and this method is considered to be inadequate for estimating the function being sought. Characteristics of bivariate interpolation are summarized, highlighting kriging limitations as an interpolation method. Limitations are similar to those of inverse distance weighted observations interpolation. Attention is drawn to the local bias of kriging and misplaced claims that it is an “optimal” interpolation method. The so-called “estimation variance,” interpreted as providing confidence limits for estimation of mining blocks, is shown to be meaningless as an index of local variation. The claim that geostatistics constitutes a “new science” is examined in detail. Such novelties as exist in the method are shown to transgress accepted principles of scientific inference. Stochastic modeling in general is discussed, and purposes of the approach emphasized. For the purpose of detailed quantitative assessment it can provide only prediction qualified by hypothesis at best. Such an approach should play no part in ore deposit assessment where the need is for local detailed inventories; these can only be achieved properly through local deterministic methods, where prediction is purely deductive.     

基于两种插值算法的三维地质建模对比

三维地质建模过程中,插值算法对模型准确性具有显著影响。为评价不同插值算法对三维建模准确性的影响,本文选取反距离权重插值法和自然邻域插值法开展对比研究。通过理论分析和案例研究,从统计学原理、插值误差和可视化效果等三方面进行了对比分析。结果表明：反距离权重插值法在建模中精度较高,适应面更广。与自然邻域插值法相比,反距离权重插值法更加适用于地层缺失严重的层位,能够更好地保留地层缺失的特征;同时,反距离权重插值法能够更好地处理断层构造,对于地层的错断起伏情况表现效果更好;反距离权重插值法在沉积地层中误差更小,与实际情况更接近。     

Application of kriging technique to areal precipitation mapping in Arizona

The classical methods for interpolating and spatial averaging of precipitation fields fail to quantify the accuracy of the estimate. On the other hand, kriging is an interpolation method for predicting values of regionalized variables at points (punctual kriging) or average values over an area (block kriging).This paper demonstrates the use of the kriging method for mapping and evaluating precipitation data for the State of Arizona. Using 158 rain gauge stations with 30 years or more of record, the precipitation over the state has been modeled as a realization of a two dimensional random field taking into consideration the spatial variability conditions.Three data sets have been used: (1) the mean annual precipitation over the state; (2) the mean summer rainy season; and (3) the mean winter rainy season. Validation of the empirical semi-variogram for a constant drift case indicated that the exponential model was appropriate for all the data sets. In addition to a global kriging analysis, the data have been examined under an anisotropic assumption which reflects the topographic structure of the state.     

An Alternative Measure of the Reliability of Ordinary Kriging Estimates

This paper presents an interpolation variance as an alternative to the measure of the reliability of ordinary kriging estimates. Contrary to the traditional kriging variance, the interpolation variance is data-values dependent, variogram dependent, and a measure of local accuracy. Natural phenomena are not homogeneous; therefore, local variability as expressed through data values must be recognized for a correct assessment of uncertainty. The interpolation variance is simply the weighted average of the squared differences between data values and the retained estimate. Ordinary kriging or simple kriging variances are the expected values of interpolation variances; therefore, these traditional homoscedastic estimation variances cannot properly measure local data dispersion. More precisely, the interpolation variance is an estimate of the local conditional variance, when the ordinary kriging weights are interpreted as conditional probabilities associated to the n neighboring data. This interpretation is valid if, and only if, all ordinary kriging weights are positive or constrained to be such. Extensive tests illustrate that the interpolation variance is a useful alternative to the traditional kriging variance.     

Advantages of using the kriging interpolator to estimate the gravity surface, comparison and spatial variability of gravity data in the El Kef-Ouargha region (northern Tunisia)

A comparative study on two interpolator methods for Bouguer anomaly mapping was performed in the El Kef-Ouargha region, Tunisia. In particular, we compare the results of minimum curvature gridding and kriging. We undertake a geostatistical study using both anisotropic and isotropic analysis for studding the spatial variability and kriging the random variable: “Bouguer anomaly”. The model parameters of each method are calculated to compare the results. The anisotropic analysis most closely approximates the measured data. Consequentially, the Bouguer anomaly and vertical gradient maps using anisotropic kriging have the best correlation with the geological map, best represent density distributions at depth and most closely approximate the measured anomaly amplitudes.     

Incorporation of principal component analysis,geostatistical interpolation approaches and frequency-space-based models for portraying the Cu-Au geochemical prospects in the Feizabad district,NW Iran

Delineation of mineralization-related geochemical anomalies of stream sediment data is an essential stage in regional geochemical exploration. In this study, principal component analysis (PCA) was applied to 12 selected elements to acquire a multi-element geochemical signature associated with Cu-Au mineralization in Feizabad district, NE Iran. The spatial distribution of enhanced multi-element geochemical signature of the second component (PC2) was modeled by different geostatistical procedures including variogram calculation, ordinary kriging (OK) and inverse distance weighting (IDW) interpolation techniques. Concentration-area (C-A) fractal and U-spatial statistics models were then applied to the continuous-value interpolated models for delineation of geochemical anomalies. Quantitative comparison of results based on the known mineral occurrences in the study area was carried out using normalized density index and success-rate curves. All generated models represent a high positive relation with known Cu (±Au) deposits in the study area, although, comparison of the results revealed that the OK-based U-spatial statistics model was superior to the rest of models. Besides, the low, moderate and high-intensity anomalies are spatially associated with geological-structural features in the study area.     

Large-scale spatial interpolation of soil pH across the Loess Plateau, China

Soil pH plays an important role in biogeochemical processes in soils. The spatial distribution of soil pH provides basic and useful information relevant to soil management and agricultural production. To obtain an accurate distribution map of soil pH on the Loess Plateau of China, 382 sampling sites were investigated throughout the region and four interpolation methods, i.e., inverse distance weighting (IDW), splines, ordinary kriging, and cokriging, were applied to produce a continuous soil pH surface. In the study region, soil pH values ranged from 6.06 to 10.76, with a mean of 8.49 and a median of 8.48. Land use type had a significant effect (p < 0.01) on soil pH; grassland soils had higher pHs than cropland and forestland soils. From a regional perspective, soil pH showed weak variation and strong spatial dependence, indicated by the low values of the coefficient of variation (0.05) and the nugget-to-sill ratios (<0.25). Indices of cross-validation, i.e., average error, mean absolute error, root mean square error, and model efficiency coefficient were used to compare the performance of the four different interpolation methods. Kriging methods interpolated more accurately than IDW and splines. Cokriging performed better than ordinary kriging and the accuracy was improved using soil organic carbon as an auxiliary variable. Regional distribution maps of soil pH were produced. The southeastern part of the region had relatively low soil pH values, probably due to higher precipitation, leaching, and higher soil organic matter contents. Areas of high soil pH were located in the north of the central part of the region, possibly associated with the salinization of sandy soils under inappropriate irrigation practices in an arid climate. Map accuracy could be further improved using new methods and incorporating other auxiliary variables, such as precipitation, elevation, terrain attributes, and vegetation types.     

Generalized sequential kriging (GSK) for heterogeneous cell size: property modeling without matrices

Sequential kriging avoids the use of matrices and resolves the issue of unstable solutions. It allows for stepwise ways to get joint estimations and cosimulations that are equivalent to the simultaneous solution. The approach is proposed as the solution for geocellular modeling with variable cell size from heterogeneous structural properties (HSPs) as required for modeling with structural constraints. Rock properties are controlled by structural domains, regions, and structural geology parameters. In some cases, rock properties are cross-correlated to formation thickness, curvature of structures, and other structural attributes. Cell thickness may be proportional to formation thickness and may enter as a conditioning property in the estimation of rock property parameters for simulation. In addition, cell volume controls the upscaling of covariance structures (i.e., regularized variograms). Structural properties are priorly modeled. Perturbation response functions (PRFs) are computed for each cell vs all possible sample point locations to facilitate sequential kriging. Upscaled PRFs are modified following conditional updating after each new data value is included in the estimation of parameters. Generalized sequential kriging is expected to become the main tool for real-time spatial modeling of 3D cellular models with HSP. In addition, some new developments related to the sequential kriging algorithm are included. Sequential kriging can be used for the estimation of parameters for simulation in the so-called unstructured grids.     

The equivalence of predictions from universal kriging and intrinsic random-function kriging

A proof is provided that the predictions obtained from kriging based on intrinsic random functions of orderk are identical to those obtained from anappropriate universal kriging model. This is a theoretical result based on known variability measures. It does not imply that people performing traditional universal kriging will get the same predictions as those using intrinsic random functions, because traditionally these methods differ in how variability is modeled. For intrinsic random functions, the same proof shows that predictions do not depend on the specific choice of the generalized covariance function. It is argued that the choice between these methods is really one of modeling and estimating the variability in the data.     

Spatial interpolation of severely skewed data with several peak values by the approach integrating kriging and triangular irregular network interpolation

It was not unusual in soil and environmental studies that the distribution of data is severely skewed with several high peak values, which causes the difficulty for Kriging with data transformation to make a satisfied prediction. This paper tested an approach that integrates kriging and triangular irregular network interpolation to make predictions. A data set consisting of total Copper (Cu) concentrations of 147 soil samples, with a skewness of 4.64 and several high peak values, from a copper smelting contaminated site in Zhejiang Province, China. The original data were partitioned into two parts. One represented the holistic spatial variability, followed by lognormal distribution, and then was interpolated by lognormal ordinary kriging. The other assumed to show the local variability of the area that near to high peak values, and triangular irregular network interpolation was applied. These two predictions were integrated into one map. This map was assessed by comparing with rank-order ordinary kriging and normal score ordinary kriging using another data set consisting of 54 soil samples of Cu in the same region. According to the mean error and root mean square error, the approach integrating lognormal ordinary kriging and triangular irregular network interpolation could make improved predictions over rank-order ordinary kriging and normal score ordinary kriging for the severely skewed data with several high peak values.     

Evaluation of spatial interpolation methods for groundwater level in an arid inland oasis,northwest China

The Wuwei oasis, situated in the upper reaches of the Shiyang River basin in the arid inland of northwest China, is intensively cultivated using both groundwater and irrigation water originating from the Qilian Mountains. Groundwater levels are declining due to overuse of irrigation water. To estimate the decline over the entire Wuwei oasis, eight different interpolation methods were used for interpolating groundwater levels over 3 years, i.e. starting in 1983, followed by 1988 and ending with 1992. Cross-validation and orthogonal-validation were applied to evaluate the accuracy of the different methods. Root mean squared error and the correlation coefficient (R ²) were calculated for each of the interpolation methods and years. Three kriging methods (simply, ordinary, and universal) gave the best fit. Modified ordinary kriging was found better than simple and universal kriging methods with a smaller number of points having large differences (>50 m) between estimated and predicted values. Based on the groundwater surfaces determined by the ordinary kriging as modified by Yamamoto, the groundwater decline was found from 1983 to 1992 to be a modest 2.1 m in average.     

Effects of search parameters on kriged reserve estimates

Summary Reliable ore reserve estimates for deposits with highly skewed grade distributions are difficult tasks to perform. Although some recent geostatistical techniques are available to handle problems with these estimations, ordinary kriging or conventional interpolation methods are still widely used to estimate the ore reserves for such deposits. The estimation results can be very sensitive to the search parameters used during the interpolation of grades with these methods.This paper compares the ore reserve estimates from ordinary kriging using several cases in which certain search parameters are varied. The comparisons are extended to different mineralizations to show the changing effects of these parameters.     

半变异函数及取样间距对克里金法在海洋地层分析中的影响研究

以浅剖数据为源数据,钻孔实测数据为验证数据,利用普通克里金法对海底地层厚度进行空间插值得到地层分布特征,采用3种半变异函数模型和不同取样间距对某井场3组地层厚度进行普通克里金插值并验证其插值效果。结果表明:普通克里金是一种有效的海底地层厚度预测方法;结构分析最佳的模型不一定是误差最小的模型,应对不同模型下的插值结果进行综合分析来选择最合适的模型,并提出球状模型在该井场厚度估计中最优,高斯模型次之;对于球状模型,增大取样间距对地层厚度变化剧烈的地层回归效果影响较小,对地层厚度变化不大的地层回归效果影响较大;同时,SE预测值变化率分析表明对于地层厚度变化剧烈的地层,减小取样间距可以大幅度地减少插值误差,而对于地层厚度变化不大的地层,减小取样间距对插值精度提高的意义不大。     

Comparison between different kriging estimators

Six different geostatistical estimators (linear kriging, lognormal kriging, and disjunctive kriging, each with and without a nonbias, i.e., universality condition) were compared using data from a polymetallic deposit in Algeria. The differences between estimators with and without the nonbias condition were far more pronounced than between the different kriging methods. This highlights the importance of choosing an appropriate stationarity model for the data. The criterion concerning kriging weight of the mean in simple kriging, proposed by Remacre (1984, 1987) and Rivoirard (1984) was found to be helpful for determining blocks where the choice of the stationarity hypothesis was critical.     

Comparison of statistical prediction methods for characterizing the spatial variability of apparent electrical conductivity in coastal salt-affected farmland

Soil salinity has been known to be problematic to land productivity and environment in the lower Yellow River Delta due to the presence of a shallow, saline water table and marine sediments. Spatial information on soil salinity has gained increasing importance for the demand of management and sustainable utilization of arable land in this area. Apparent electrical conductivity, as measured by electromagnetic induction instrument in a fairly quick manner, has succeeded in mapping soil salinity and many other soil physical and chemical properties from field to regional scales. This was done based on the correlation that existed between apparent electrical conductivity and many other soil properties. In this paper, four spatial prediction methods, i.e., local polynomial, inverse distance weighed, ordinary kriging and universal kriging, were employed to estimate field-scale apparent electrical conductivity with the aid of an electromagnetic induction instrument (type EM38). The spatial patterns estimated by the four methods using EM38 survey datasets of various sample sizes were compared with those generated by each method using the entire sample size. Spatial similarity was evaluated using difference index (DI) between the maps created using various sample sizes (i.e., target maps) and the maps generated with the entire sample size (i.e., the reference map). The results indicated that universal kriging had the best performance owing to the inclusion of residuals and spatial detrending in the kriging system. DI showed that spatial similarity between the target and reference maps of apparent electrical conductivity decreased with the reduction in sample size for each prediction method. Under the same reduction in sample size, the method retaining the most spatial similarity was universal kriging, followed by ordinary kriging, inverse distance weighed, and local polynomial. Approximately, 70 % of total survey data essentially met the need for retaining 90 % details of the reference map for universal kriging and ordinary kriging methods. This conclusion was that OK and UK were two most appropriate methods for spatial estimation of apparent electrical conductivity as they were robust with the reduction in sample size.