Spatial structure analysis of regionalized compositions

Like compositions in general, regionalized compositions present the problem of spurious spatial correlation. To avoid this problem, this paper uses the additive-logratio transformation of regionalized compositions, following techniques introduced over the last few years for the statistical analysis of compositional data. It leads to an appropriate definition of a spatial covariance structure to describe spatial dependence between regionalized variables subject to constant-sum constraints in the case of weak stationarity. To illustrate stated problems, simulated data are used.     

Compositional Kriging: A Spatial Interpolation Method for Compositional Data

Compositional data are very common in the earth sciences. Nevertheless, little attention has been paid to the spatial interpolation of these data sets. Most interpolators do not necessarily satisfy the constant sum and nonnegativity constraints of compositional data, nor take spatial structure into account. Therefore, compositional kriging is introduced as a straightforward extension of ordinary kriging that complies with these constraints. In two case studies, the performance of compositional kriging is compared with that of the additive logratio-transform. In the first case study, compositional kriging yielded significantly more accurate predictions than the additive logratio-transform, while in the second case study the performances were comparable.     

不同采样点数量下土壤有机质含量空间预测方法对比

运用普通克里格、泛克里格、协同克里格和回归克里格4种方法,结合由DEM获取的高程因子以及土壤全氮和阳离子交换量（CEC）,预测了黑龙江省海伦市耕地有机质含量的空间分布。不同样点数量下海伦市土壤有机质含量的空间变异结构分析表明,样点数量多并不一定能够识别土壤有机质含量的结构性连续组分,最优化的布置采样点位置可能比单纯增加...     

Estimation of regionalized phenomena by geostatistical methods: lake acidity on the Canadian Shield

 This paper describes a geostatistical technique based on conditional simulations to assess confidence intervals of local estimates of lake pH values on the Canadian Shield. This geostatistical approach has been developed to deal with the estimation of phenomena with a spatial autocorrelation structure among observations. It uses the autocorrelation structure to derive minimum-variance unbiased estimates for points that have not been measured, or to estimate average values for new surfaces. A survey for lake water chemistry has been conducted by the Ministère de l'Environnement du Québec between 1986 and 1990, to assess surface water quality and delineate the areas affected by acid precipitation on the southern Canadian Shield in Québec. The spatial structure of lake pH was modeled using two nested spherical variogram models, with ranges of 20 km and 250 km, accounting respectively for 20% and 55% of the spatial variation, plus a random component accounting for 25%. The pH data have been used to construct a number of geostatistical simulations that produce plausible realizations of a given random function model, while 'honoring' the experimental values (i.e., the real data points are among the simulated data), and that correspond to the same underlying variogram model. Post-processing of a large number of these simulations, that are equally likely to occur, enables the estimation of mean pH values, the proportion of affected lakes (lakes with pH≤5.5), and the potential error of these parameters within small regions (100 km×100 km). The method provides a procedure to establish whether acid rain control programs will succeed in reducing acidity in surface waters, allowing one to consider small areas with particular physiographic features rather than large drainage basins with several sources of heterogeneity. This judgment on the reduction of surface water acidity will be possible only if the amount of uncertainty in the estimation of mean pH is properly quantified. Received: 3 March 1997 · Accepted: 16 November 1998     

A comparison of field and laboratory analytical methods for radon site investigation

Soil-gas radon measurements provide a valuable tool in assessing probable indoor radon levels on a regional basis. However, in Great Britain, seasonal weather changes can cause large changes in soil-gas radon concentration. Although this does not significantly constrain systematic radon potential mapping programmes, it does cause difficulties in responding to ad-hoc requests for site-specific radon investigations. The relationship between soil-gas radon and gamma spectrometry measurements made in the field with radon released from a representative sample of soil in the laboratory has been investigated as part of a program to develop a method of radon potential mapping and site investigation which can be used at any time of the year. Multiple soil and soil-gas samples were collected from sites underlain by bedrocks with widely varying radon potentials. For each geological unit, sites both free of and covered by glacial drift deposits were sampled. Soil and soil-gas samples were taken at the same depth of 60–100 cm. The effectiveness of these radon site investigation procedures has been evaluated by studying the relationship between the soil-gas radon, gamma spectrometry and radon emanation data with an independent estimate of the radon risk. The geologic radon potential (GEORP), which is the proportion of existing dwellings which exceed the UK radon Action Level (200 Bq m⁻³) for a particular combination of solid and drift geology within a defined geographic area, has been used for this study as the independent estimate of radon risk. Soil-gas radon, radon emanation and eU (equivalent uranium by field γ spectrometry) are all good geochemical indicators of radon risk (GEORP) in Derbyshire but only soil-gas radon correlates significantly with GEORP in Northamptonshire. Radon in soil gas discriminates more effectively between sites with different radon potential in Northamptonshire if soil permeability is also taken into account. In general, measurement of soil-gas radon in the field provides the most universally applicable indicator of radon potential. If soil-gas radon concentrations cannot be determined because of climatic factors, for example when the soil profile is waterlogged, measurement of radon emanation in the laboratory or measurement of eU can be used as radon potential indicators in some geological environments. This applies particularly in areas where the soil composition rather than the composition and permeability of the underlying rock or superficial deposits are the dominant controls of radon potential. It appears, therefore, that it may be necessary to use different radon site investigation methods according to the specific factors controlling radon emanation from the ground. In some cases no method will provide a reliable indicator of radon risk under unfavourable climatic conditions.     

Estimation of semivariogram parameters and evaluation of the effects of data sparsity

Semivariogram parameters are estimated by a weighted least-squares method and a jackknife kriging method. The weighted least-squares method is investigated by differing the lag increment and maximum lag used in the fit. The jackknife kriging method minimizes the variance of the jackknifing error as a function of semivariogram parameters. The effects of data sparsity and the presence of a trend are investigated by using 400-, 200-, and 100-point synthetic data sets. When the two methods yield significantly different results, more data may be needed to determine reliably the semivariogram parameters, or a trend may be present in the data.     

A comparison of three gravity inversion methods for crustal thickness modelling in Tibet plateau

Crustal thickness can be determined by gravimetric methods based on different assumptions, e.g. by isostatic hypotheses. Here we compare three gravimetric inversion methods to estimate the Moho depth. Two Moho models based on the Vening Meinesz-Moritz hypothesis and one by using Parker-Oldenburg’s algorithm, which are investigated in Tibet plateau. The results are compared with CRUST2.0, and it will be presented that the estimated Moho depths from the Vening Meinesz-Moritz model will be better than the Parker-Oldenburg’s algorithm.     

Estimation of the generalized covariance function: II. A response surface approach

A procedure is proposed that employs first-moment estimation (kriging), cross-validation, and response surface analysis to estimate parameters of a generalized covariance function. Results from application of this procedure to two data sets are given.     

对数比变换和偏最小二乘法在地球化学组合异常提取中的应用--以湘西北铅锌矿为例

以湘西北水系沉积物测量常量元素含量数据为研究对象,针对湘西北铅锌矿床的层控特征,以研究区含矿层地球化学常量元素含量的统计规律为基础,利用对数比变换( alr、clr、ilr变换)对原始数据进行预处理,并在此基础上进行偏最小二乘( PLS)降维分析,提取了地球化学组合异常,并与常用的主成分分析、因子分析的结果进行对比。分析结果表明,等距对数比变换( ilr变换)相对于其他变换方法有明显的优势,变换后的数据对称特征明显,更趋于正态分布。相对于其他两种降维方法,偏最小二乘法提取的元素组合与研究区铅锌矿的地球化学基本特征更为接近,组合异常与已知矿床分布和构造特征耦合良好,与研究区地质特征吻合,比主成分分析与因子分析得到的结果更易于地质解释。     

三种插值方法在降水等值线绘制中的对比分析

为对比不同插值方法在宁夏回族自治区的适用性,基于宁夏及其周边地区77个气象台站1956—2016年平均降水量数据,运用泰森多边形法、距离平均反比法、克里金插值法对降雨量数据进行空间插值和交叉验证,从插值方法的稳健性和准确度出发进行对比分析,得出适用于宁夏回族自治区的插值方法。研究表明：（1）从准确度来看,克里金插值法与观测值走向、位置较为相似;从稳健性来看,克里金插值法MAE平均值为4.43 mm、RMSE平均值为10.36 mm,显著小于泰森多边形法和距离平均反比法。（2）整体来看,克里金插值法更为理想,可以用于宁夏回族自治区降水量空间插值分析。     

激电法在西北某油田上的应用效果

本文依据前人的理论基础,在总结分析激电法找油的地球物理前提条件下,开展了该方法的试验兼生产工作,取得了令人满意的效果,并进一步证实了激电法找油的异常模式。     

A new approach to null correlations of proportions

Much work on the statistical analysis of compositional data has concentrated on the difficulty of interpreting correlations between proportions with an assortment of tests for nullcorrelations, for independence except for the constraint, F-independence of bounded variables, neutrality in the mean and in the median. This paper questions the appropriateness of characterizing the dependence structure of proportions in terms of such concepts, suggests an alternative method of modeling, develops necessary distribution theory and tests, and illustrates the methodology in applications.     

Estimating fractal dimension of profiles: A comparison of methods

This paper examines the characteristics of four different methods of estimating the fractal dimension of profiles. The semi-variogram, roughness-length, and two spectral methods are compared using synthetic 1024-point profiles generated by three methods, and using two profiles derived from a gridded DEM and two profiles from a laser-scanned soil surface. The analysis concentrates on the Hurst exponent H,which is linearly related to fractal dimension D,and considers both the accuracy and the variability of the estimates of H.The estimation methods are found to be quite consistent for Hnear 0.5, but the semivariogram method appears to be biased for Happroaching 0 and 1, and the roughness-length method for Happroaching 0. The roughness-length or the maximum entropy spectral methods are recommended as the most suitable methods for estimating the fractal dimension of topographic profiles. The fractal model fitted the soil surface data at fine scales but not at broad scales, and did not appear to fit the DEM profiles well at any scale.     

塔河油田叠后地震资料高保真处理技术

塔河油田岩溶储集体(洞穴、礁滩、孔洞、裂缝或裂隙等)的地震波场特征是弱反射背景之上的强弱不同的绕射波,因此,绕射波准确归位成像和振幅保持是研究岩溶储集体的基础.笔者从塔河油田地震资料高保真处理实践中,总结出了一套适应该区碳酸盐岩储层预测要求的以叠后高精度偏移算法和振幅保持处理为主的叠后高保真处理技术,取得显著的效果.     

激电法在西北某油田上的应用效果

本文依据前人的理论基础，在总结分析激电法找油的地球物理前提条件下，开展了该方法的试验兼生产工作，取得了令人满意的效果，并进一步证实了激电法找油的异常模式。     

DDA法和Fisher判别法在潜在滑坡判识中的应用比较

选择地貌特征（包括地形坡度、沟谷切割状况、洼地封闭状况、前缘临空状况）、滑面特征（包括滑面倾角、滑面形态、滑面阻抗比）、滑体结构和近期活动迹象4大类9个因子作为古（老）滑坡类潜在滑坡的判识指标,以三峡库区和岷江上游地区相同数量的典型滑坡作为学习样本,依据距离判别分析（DDA）和Fisher线性判别分析（FLDA）方法,分别建立了两个研究区潜在滑坡的DDA判识模型和FLDA判识模型。实例分析结果表明,DDA方法对两个研究区学习样本和测试样本的误判率均为0,而FLDA方法对两个研究区测试样本的误判率分别为5.56%和11.11%,对学习样本的错判率分别为11.1%和0%。因此,在三峡库区和岷江上游地区潜在滑坡的判识中,DDA法比FLDA法判识准确性更高、适用性更强。     

Evaluation of vertical shaft stability in underground mines: comparison of three weight methods with uncertainty theory

Natural Hazards - Shaft stability evaluation (SSE) is one of the most crucial and important tasks in view of the role of vertical shaft in mining engineering, the accuracy of which determines the...     

A comparison of discrete element and micromechanical methods for determining the effective elastic properties of geomaterials

Data generated from scanning electron microscopy images of oil-shale geomaterials are used to generate configurations and to acquire parameters required for use in homogenization schemes for the determination of the effective elastic properties of the samples. Two alternative homogenization methods are employed: numerical simulation using the Discrete Element Method and the Polycrystalline Self-Consistent Method from micromechanics. The schemes give rise to predictions of the effective elastic properties that are in very good agreement.