Efficient updating of kriging estimates and variances   总被引：2，自引：0，他引：2  

R. J. Barnes  A. G. Watson 《Mathematical Geology》1992,24(1):129-133

This short note presents a method for efficiently updating ordinary kriging estimates and variances when one or more additional samples are incorporated into the kriging system. First, the foundation linear algebra result is presented. Then the update equations are derived. Finally, an illustrative application of updating is briefly discussed.  相似文献   

A variance-ratio test for supporting a variable mean in kriging     

Peter K. Kitanidis 《Mathematical Geology》1997,29(3):335-348

This article discusses the issue of whether to use a variable mean and describes a test that can be used to evaluate whether it is justified to add terms to the drift (deterministic part) of a geostatistical model. The basic model could be the intrinsic one, where the deterministic part is a constant, and the alternate model could be any model that includes a constant term in the expression for the drift. Also, differences between constant- and variable-mean models are discussed.  相似文献   

Foliation fields and 3D cartography in geology: Principles of a method based on potential interpolation   总被引：4，自引：0，他引：4  

Christian Lajaunie  Gabriel Courrioux  Laurent Manuel 《Mathematical Geology》1997,29(4):571-584

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.  相似文献   

Is the ocean floor a fractal?   总被引：1，自引：0，他引：1  

Ute C. Herzfeld  Isaac I. Kim  John A. Orcutt 《Mathematical Geology》1995,27(3):421-462

The topographic structure of the ocean bottom is investigated at different scales of resolution to answer the question: Can the seafloor be described as a fractal process? Methods from geostatistics, the theory of regionalized variables, are used to analyze the spatial structure of the ocean floor at different scales of resolution. The key to the analysis is the variogram criterion: Self-similarity of a stochastic process implies self-similarity of its variogram. The criterion is derived and proved here: it also is valid for special cases of self-affinity (in a sense adequate for topography). It has been proposed that seafloor topography can be simulated as a fractal (an object of Hausdorff dimension strictly larger than its topological dimension), having scaling properties (self-similarity or self-affinity). The objective of this study is to compare the implications of these concepts with observations of the seafloor. The analyses are based on SEABEAM bathymetric data from the East Pacific Rise at 13°N/104°W and at 9°N/104°W and use tracks that run both across the ridge crest and along the ridge flank. In the geostatistical evaluation, the data are considered as a stochastic process. The spatial continuity of this process is described by variograms that are calculated for different scales and directions. Applications of the variogram criterion to scale-dependent variogram models yields the following results: Although the seafloor may be a fractal in the sense of the definition involving the Hausdorff dimension, it is not self-similar, nor self-affine (in the given sense). Mathematical models of scale-dependent spatial structures are presented, and their relationship to geologic processes such as ridge evolution, crust formation, and sedimentation is discussed.  相似文献   

多点地质统计学建模研究进展   总被引：1，自引：0，他引：1       下载免费PDF全文

陈欢庆  李文青  洪垚 《高校地质学报》2018,24(4):593

从精细油藏描述中地质建模的意义和现状入手,介绍了多点地质统计学建模研究现状及其与传统地质建模方法的差 异。以辽河盆地西部凹陷某蒸汽驱试验区为例,分析了多点地质统计学建模中建模的基础、训练图像的建立、多点地质统计学 建模与传统地质建模相比所具有的优势等内容。指出多点地质统计学在井间预测方面具有明显优于其他传统建模方法的特 点。在文献调研基础上,结合自身工作实践,探讨了多点地质统计学建模目前存在的问题和未来的发展方向。指出未来多点 地质统计学建模的发展方向主要包括多信息综合地质成因分析基础上的训练图像获取、多点地质统计学算法进行改进和完善 和多点地质统计学建模方法应用领域的扩大等。  相似文献   

渭河河漫滩宏观弥散度计算          下载免费PDF全文

张小筱  李修成 《水文》2017,37(5):21-26

宏观弥散度的确定对地下水溶质运移至关重要。利用竖管法野外现场测定900个点的垂向饱和渗透系数Kv值,然后进行统计分布检验,运用地质统计学中负指数模型在有块金效应影响的条件下拟合出ln Kv的空间相关尺度λ值,进而利用λ与宏观弥散度关系式求解渭河漫滩非饱和带的宏观弥散度。结果表明,研究区整体Kv平均值为0.012cm/s,方差为1.33,表明Kv空间差异性较大;垂向渗透系数基本服从对数正态分布规律,满足宏观弥散度计算条件;估算整体纵向宏观弥散度为21.27m,各排计算结果约2.11~4.65m。  相似文献   

福州市农田土壤养分空间变异特征   总被引：4，自引：0，他引：4  

陈桂香  高灯州  曾从盛  王维奇 《地球信息科学学报》2017,19(2):216-224

研究农田土壤养分空间分布特征及变异规律,分析地形因子对土壤养分空间分布的影响,对土壤养分的有效利用与管理具有重要意义。本文采用GIS和地统计学相结合的方法,揭示东南丘陵区福州市农田土壤养分（有机质、碱解氮、有效磷和速效钾）的空间分布特征和变异规律,并进一步探讨土壤养分含量与地形因子（地形起伏度、高程、地形湿度指数、沉积运输指数和坡度）之间的相关性。结果表明：福州市农田有机质、碱解氮、有效磷和速效钾含量范围分别为1.10~89.5 g/kg、 1.00~461 mg/kg、0.300~298 mg/kg、和4.00~399 mg/kg,其变异系数范围为35.3%~99.0%,均属于中等变异水平;区域内耕地土壤养分存在明显的空间丰缺差异,大部分地区有机质和有效磷含量较为丰富,碱解氮含量处于中等偏上水平,速效钾含量相对较低;有机质、碱解氮、有效磷和速效钾的块金系数分别为32.0%、37.3%、50.0%和50.0%,均呈中等强度空间自相关,表明同时受到结构和随机因素控制;有机质和碱解氮空间自相关尺度较大,且在步长小于0.3 km时,各方向（0°、45°、90°和135°）变化平稳,为各向同性,而有效磷和速效钾变程较小,且各方向变化较复杂,为各向异性。这些结果表明,政府需加强施肥指导,保持氮肥施用量的同时,合理增加钾肥,适当降低有机肥和磷肥施用。此外,在后续调查采样时,样点布设要考虑密度和方向性,适当加密有效磷和速效钾的采样,而有机质和碱解氮采样可以在此基础上适当减少样点。  相似文献   

虎头崖铅锌多金属矿床矿体三维空间变化结构及分布规律          下载免费PDF全文

曹芳  毛先成  邹艳红  赵莹  邹品娟  吴泽权 《地质找矿论丛》2013,28(1):41-49

文章利用Vulcan和Datamine三维可视化软件建立了虎头崖铅锌多金属矿床矿体的三维实体模型、块体模型及地层三维实体模型;运用地质统计学方法分析了虎头崖矿区67线以东和117线以西矿体的Pb,Zn品位空间变化结构,得出搜索椭球体,实现了虎头崖铅锌多金属矿床矿体、地层形态和铅锌品位变化形态的空间分布展示;结合三维建模结果和铅锌矿体空间变化结构,总结出矿体空间分布规律,并用三维可视化软件对矿体分布规律进行显示,指导隐伏矿体的预测,寻找有利的找矿部位。  相似文献   

On the estimation of scale of fluctuation in geostatistics     

《Georisk: Assessment and Management of Risk for Engineered Systems and Geohazards》2013,7(2):129-140

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.  相似文献   

Optimization of mobile radioactivity monitoring networks     

G. B.M. Heuvelink  Z. Jiang  S. De Bruin  C. J.W. Twenhöfel 《International journal of geographical information science》2013,27(3):365-382

In case of a nuclear accident, decision makers rely on high-resolution and accurate information about the spatial distribution of radioactive contamination surrounding the accident site. However, the static nuclear monitoring networks of many European countries are generally too coarse to provide the desired level of spatial accuracy. In the Netherlands, authorities are considering a strategy in which measurement density is increased during an emergency using complementary mobile measuring devices. This raises the question, where should these mobile devices be placed? This article proposes a geostatistical methodology to optimize the allocation of mobile measurement devices, such that the expected weighted sum of false-positive and false-negative areas (i.e. false classification into safe and unsafe zones) is minimized. Radioactivity concentration is modelled as the sum of a deterministic trend and a zero-mean spatially correlated stochastic residual. The trend is defined as the outcome of a physical atmospheric dispersion model, NPK-PUFF. The residual is characterized by a semivariogram of differences between the outputs of various NPK-PUFF model runs, designed to reflect the effect of uncertainty in NPK-PUFF meteorological inputs (e.g. wind speed, wind direction). Spatial simulated annealing is used to obtain the optimal monitoring design, in which accessibility of sampling sites (e.g. distance to roads) is also considered. Although the methodology is computationally demanding, results are promising and the computational load may be considerably reduced to compute optimal mobile monitoring designs in nearly real time.  相似文献