A Combined Three-Dimensional Geological-Geostatistical-Numerical Model of Underground Excavations in Rock   总被引：1，自引：1，他引：1  

M. Stavropoulou  G. Exadaktylos  G. Saratsis 《Rock Mechanics and Rock Engineering》2007,40(3):213-243

Summary. This paper exploits geological and borehole geotechnical data obtained in the exploratory phase of a tunneling project to investigate in a first place if the kriging interpolation scheme may effectively reproduce the spatial variability of rock mass quality (Rock Mass Rating, RMR) in the vicinity of tunnels. For this purpose a quick solver in Fortran has been developed that performs variography analysis of 3D spatial data, fast kriging estimations of RMR between borehole sampling locations at the centroids of the elements of the numerical model, and model validation. For the purpose of an integrated underground excavation design, a step further is made by incorporating into the 3D mechanical numerical model of the rock mass, the three-dimensional (3D) solid geological model, thus coupling the geology with the ground (geotechnical) model (i.e. each element of the numerical model is assigned a geological material). The mechanical properties of each finite difference cell (or Representative Elementary Volume) of the ground model were then prescribed according to its geological type, the spatial heterogeneity of the rock mass expressed quantitatively with the kriging model, and the upscaling calculations of the mechanical properties of the intact rocks determined in the laboratory, based on the size-effect (strength dependence on size) and Damage Theory. Furthermore, a preliminary numerical simulation of the advance of unsupported tunnels in the model of the heterogeneous rock mass was performed for illustration purposes.  相似文献   

Hybrid Estimation of Semivariogram Parameters   总被引：1，自引：0，他引：1  

Hao Zhang  Dale L. Zimmerman 《Mathematical Geology》2007,39(2):247-260

Two widely used methods of semivariogram estimation are weighted least squares estimation and maximum likelihood estimation. The former have certain computational advantages, whereas the latter are more statistically efficient. We introduce and study a “hybrid” semivariogram estimation procedure that combines weighted least squares estimation of the range parameter with maximum likelihood estimation of the sill (and nugget) assuming known range, in such a way that the sill-to-range ratio in an exponential semivariogram is estimated consistently under an infill asymptotic regime. We show empirically that such a procedure is nearly as efficient computationally, and more efficient statistically for some parameters, than weighted least squares estimation of all of the semivariogram’s parameters. Furthermore, we demonstrate that standard plug-in (or empirical) spatial predictors and prediction error variances, obtained by replacing the unknown semivariogram parameters with estimates in expressions for the ordinary kriging predictor and kriging variance, respectively, perform better when hybrid estimates are plugged in than when weighted least squares estimates are plugged in. In view of these results and the simplicity of computing the hybrid estimates from weighted least squares estimates, we suggest that software that currently estimates the semivariogram by weighted least squares methods be amended to include hybrid estimation as an option.  相似文献   

Geostatistical incorporation of spatial coordinates into supervised classification of hyperspectral data   总被引：1，自引：0，他引：1  

P. Goovaerts 《Journal of Geographical Systems》2002,4(1):99-111

 This paper presents a methodology to incorporate both hyperspectral properties and spatial coordinates of pixels in maximum likelihood classification. Indicator kriging of ground data is used to estimate, for each pixel, the prior probabilities of occurrence of classes which are then combined with spectral-based probabilities within a Bayesian framework. In the case study (mapping of in-stream habitats), accounting for spatial coordinates increases the overall producer's accuracy from 85.8% to 93.8%, while the Kappa statistic rises from 0.74 to 0.88. Best results are obtained using only indicator kriging-based probabilities, with a stunning overall accuracy of 97.2%. Significant improvements are observed for environmentally important units, such as pools (Kappa: 0.17 to 0.74) and eddy drop zones (Kappa: 0.65 to 0.87). The lack of benefit of using hyperspectral information in the present study can be explained by the dense network of ground observations and the high spatial continuity of field classification which might be spurious. Received: 12 April 2001 / Accepted: 7 September 2001  相似文献   

松潘地区岩金统计成矿予测中泛克立格法的应用     

刘智光  汪西鸣 《矿物岩石》1992,(1)

川西北地区是国内重要的岩金普查靶区。本文对松潘地区的地质特征,泛克立格法原理,及在本区的应用情况作了介绍,对主要计算步骤,给出基本结果,对计算获得的有意义异常区,进行地质解释,并将此方法与趋势面分析进行对比。  相似文献   

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

Estimation of regionalized compositions: A comparison of three methods   总被引：1，自引：0，他引：1  

V. Pawlowsky  R. A. Olea  J. C. Davis 《Mathematical Geology》1995,27(1):105-127

A regionalized composition is a random vector function whose components are positive and sum to a constant at every point of the sampling region. Consequently, the components of a regionalized composition are necessarily spatially correlated. This spatial dependence—induced by the constant sum constraint—is a spurious spatial correlation and may lead to misinterpretations of statistical analyses. Furthermore, the cross-covariance matrices of the regionalized composition are singular, as is the coefficient matrix of the cokriging system of equations. Three methods of performing estimation or prediction of a regionalized composition at unsampled points are discussed: (1) the direct approach of estimating each variable separately; (2) the basis method, which is applicable only when a random function is available that can he regarded as the size of the regionalized composition under study; (3) the logratio approach, using the additive-log-ratio transformation proposed by J. Aitchison, which allows statistical analysis of compositional data. We present a brief theoretical review of these three methods and compare them using compositional data from the Lyons West Oil Field in Kansas (USA). It is shown that, although there are no important numerical differences, the direct approach leads to invalid results, whereas the basis method and the additive-log-ratio approach are comparable.  相似文献   

地质统计学方法在井间储层参数预测中的应用     

信荃麟  杜庆龙 《地质论评》1993,39(4):283-291

本文探讨了利用地质统计学中的三种克立格方法进行井间储层参数预测的原理及方法,并在某油田中应用,通过最优性检验,建立起了适合研究区地质特点的各种储层参数的最佳估计模型,揭示了储层参数在三维空间的变化规律,为建立三维储层预测模型开辟了新的方法。  相似文献   

扇形沉积体生长的分形几何特征分析     

周江羽  刘常青 《沉积学报》2000,18(1):95-99

扇形沉积体是我国内陆及沿海地区中、新生代含油气盆地中的重要油气勘探目标,对其平面分布形态和内部结构参数特征的定量表征是油气勘探部门极为关注的问题。随着分形几何理论在沉积学领域应用的日趋广泛,使人们可从新的角度来考虑这一问题的解决方法。扇形沉积体的生长是一个复杂的非线性动力学过程,且具有分形特征。分维数的大小反映了沉积体外部形态和内部结构的复杂程度,分维数的变化幅度预示着构造和沉积背景的某些信息。将分形和混沌理论相结合,也许是实施扇形沉积体非线性动力学建模和模拟的重要途经之一。  相似文献   

利用克里金技术标定新疆及附近地区P/S震相幅值比及其在地震事件识别中的应用          下载免费PDF全文

潘常周  靳平  肖卫国 《地震学报》2007,29(6):625-634

针对中国的乌鲁木齐台（WMQ）和哈萨克斯坦的马坎奇台（MAK），采用贝叶斯克里金技术建立了新疆及附近地区地震事件的震相幅值比的校正曲面，并分析了校正曲面与地震波传播路径的关系，以及校正曲面对于改善震相幅值比判据识别效果的有效性.分析结果表明，建立的校正曲面基本合理，结合传播路径差异和地震波传播规律，能够较好地解释校正曲面的起伏特征. 而且，在传播路径差异较大，使得震相幅值比经震中距校正后离散度仍然较高的情况下，利用校正曲面修正震相幅值比，可以进一步降低其离散度，从而改善判据的识别效果.经校正曲面修正后，对天然地震的误识率由16.3%下降为5.2％.   相似文献   

Kriging Regionalized Positive Variables Revisited: Sample Space and Scale Considerations   总被引：1，自引：0，他引：1  

Raimon Tolosana-Delgado  Vera Pawlowsky-Glahn 《Mathematical Geology》2007,39(6):529-558

Frequently, regionalized positive variables are treated by preliminarily applying a logarithm, and kriging estimates are back-transformed using classical formulae for the expectation of a lognormal random variable. This practice has several problems (lack of robustness, non-optimal confidence intervals, etc.), particularly when estimating block averages. Therefore, many practitioners take exponentials of the kriging estimates, although the final estimations are deemed as non-optimal. Another approach arises when the nature of the sample space and the scale of the data are considered. Since these concepts can be suitably captured by an Euclidean space structure, we may define an optimal kriging estimator for positive variables, with all properties analogous to those of linear geostatistical techniques, even for the estimation of block averages. In this particular case, no assumption on preservation of lognormality is needed. From a practical point of view, the proposed method coincides with the median estimator and offers theoretical ground to this extended practice. Thus, existing software and routines remain fully applicable.  相似文献