首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 428 毫秒
1.
A number of criteria based on kriging variance calculations may be used for infill sampling design in geologic site characterization. Searching for the best new sample locations from a set of candidate locations can result in excessive computation time if these criteria and the naive rekriging are used. The relative updated kriging estimate and variance for universal kriging estimation are demonstrated as a simple kriging estimate and variance, respectively. The updated kriging variance is demonstrated as the multiplication of two kriging variances. Using these updated kriging variance equations can increase the computational speed for selecting the best new sample locations. The application results for oil rock thickness in an oilfield indicate that minimizing the average relative updated kriging variance is a useful alternative to the other criteria based on kriging variance in optimal infill sampling design for geologic site characterization.  相似文献   

2.
Conclusions The foregoing discussion indicates that geostatistical estimation of ore deposits is not local; it is not objective; it is not sensitive to local data trends; and it is not unrestrained by the range of data values.Kriging, as an interpolation method, is a variant of IDW least squares linear fit. As such, it suffers from the limitations of all IDW linear interpolation methods that employ only data values.The estimation variance, currently used to calculate the confidence limits of values for individual mining blocks, is hypothetical and globally derived. It is more closely related to sampling density than to local variation in the data set.Geostatistical methods, of course, have a real place in ore deposit assessment, e.g. global, comparative evaluation to assist decisions on development and investment. What is questioned here is the validity of employing a global method to assess detail (mining blocks) within an ore deposit.  相似文献   

3.
Sample schemes used in geostatistical surveys must be suitable for both variogram estimation and kriging. Previously schemes have been optimized for one of these steps in isolation. Ordinary kriging generally requires the sampling locations to be evenly dispersed over the region. Variogram estimation requires a more irregular pattern of sampling locations since comparisons must be made between measurements separated by all lags up to and beyond the range of spatial correlation. Previous studies have not considered how to combine these optimized schemes into a single survey and how to decide what proportion of sampling effort should be devoted to variogram estimation and what proportion devoted to kriging An expression for the total error in a geostatistical survey accounting for uncertainty due to both ordinary kriging and variogram uncertainty is derived. In the same manner as the kriging variance, this expression is a function of the variogram but not of the sampled response data. If a particular variogram is assumed the total error in a geostatistical survey may be estimated prior to sampling. We can therefore design an optimal sample scheme for the combined processes of variogram estimation and ordinary kriging by minimizing this expression. The minimization is achieved by spatial simulated annealing. The resulting sample schemes ensure that the region is fairly evenly covered but include some close pairs to analyse the spatial correlation over short distances. The form of these optimal sample schemes is sensitive to the assumed variogram. Therefore a Bayesian approach is adopted where, rather than assuming a single variogram, we minimize the expected total error over a distribution of plausible variograms. This is computationally expensive so a strategy is suggested to reduce the number of computations required  相似文献   

4.
When estimating the mean value of a variable, or the total amount of a resource, within a specified region it is desirable to report an estimated standard error for the resulting estimate. If the sample sites are selected according to a probability sampling design, it usually is possible to construct an appropriate design-based standard error estimate. One exception is systematic sampling for which no such standard error estimator exists. However, a slight modification of systematic sampling, termed 2-step tessellation stratified (2TS) sampling, does permit the estimation of design-based standard errors. This paper develops a design-based standard error estimator for 2TS sampling. It is shown that the Taylor series approximation to the variance of the sample mean under 2TS sampling may be expressed in terms of either a deterministic variogram or a deterministic covariance function. Variance estimation then can be approached through the estimation of a variogram or a covariance function. The resulting standard error estimators are compared to some more traditional variance estimators through a simulation study. The simulation results show that estimators based on the new approach may perform better than traditional variance estimators.  相似文献   

5.
An interpolation method based on a multilayer neural network (MNN), has been examined and tested for the data of irregular sample locations. The main advantage of MNN is in that it can deal with geoscience data with nonlinear behavior and extract characteristics from complex and noisy images. The training of MNN is used to modify connection weights between nodes located in different layers by a simulated annealing algorithm (one of the optimization algorithms of the network). In this process, three types of errors are considered: differences in values, semivariograms, and gradients between sample data and outputs from the trained network. The training is continued until the summation of these errors converges to an acceptably small value. Because the MNN trained by this learning criterion can estimate a value at an arbitrary location, this method is a form of kriging and termed Neural Kriging (NK). In order to evaluate the effectiveness of NK, a problem on restoration ability of a defined reference surface from randomly chosen discrete data was prepared. Two types of surfaces, whose semivariograms are expressed by isotropic spherical and geometric anisotropic gaussian models, were examined in this problem. Though the interpolation accuracy depended on the arrangement pattern of the sample locations for the same number of data, the interpolation errors of NK were shown to be smaller than both those of ordinary MNN and ordinal kriging. NK can also produce a contour map in consideration of gradient constraints. Furthermore, NK was applied to distribution analysis of subsurface temperatures using geothermal investigation loggings of the Hohi area in southwest Japan. In spite of the restricted quantity of sample data, the interpolation results revealed high temperature zones and convection patterns of hydrothermal fluids. NK is regarded as an interpolation method with high accuracy that can be used for regionalized variables with any structure of spatial correlation.  相似文献   

6.
Average kriging variance is a standard tool used in optimization of the location of additional drill holes. However, this tool cannot distinguish between areas with different priorities. This limitation could be eliminated by using weighted average kriging variance. This paper extends the problem of optimal location to three dimensional cases, use grade as a weight and search optimum locations by simulated annealing. Weighted average kriging variance is used as objective function. The method is applied to a copper deposit. Results have shown that weighting of the estimation variance with ??grade?? is effective only when the difference among the grades estimated for different blocks is considerable.  相似文献   

7.
The cumulative semivariogram approach is proposed for modeling regionalized variables in the geological sciences. This semivariogram is defined as the successive summation of half-squared differences which are ranked according to the ascending order of distances extracted from all possible pairs of sample locations within a region. This procedure is useful especially when sampling points are irregularly distributed within the study area. Cumulative semivariograms possess all of the objective properties of classical semivariograms. Classical semivariogram models are evaluated on the basis of the cumulative semivariogram methodology. Model parameter estimation procedures are simplified with the use of arithmetic, semilogarithmic, or double-logarithmic papers. Plots of cumulative semivariogram values vs. corresponding distances may scatter along a straight line on one of these papers, which facilitates model identification as well as parameter estimation. Straight lines are fitted to the cumulative semivariogram scatter diagram by classical linear regression analysis. Finally, applications of the methodology are presented for some groundwater data recorded in the sedimentary basins of the Kingdom of Saudi Arabia.  相似文献   

8.
Sampling design optimization for spatial functions   总被引：4，自引：0，他引：4  
A new procedure is presented for minimizing the sampling requirements necessary to estimate a mappable spatial function at a specified level of accuracy. The technique is based on universal kriging, an estimation method within the theory of regionalized variables. Neither actual implementation of the sampling nor universal kriging estimations are necessary to make an optimal design. The average standard errorand maximum standard error of estimationover the sampling domain are used as global indices of sampling efficiency. The procedure optimally selects those parameters controlling the magnitude of the indices, including the density and spatial pattern of the sample elements and the number of nearest sample elements used in the estimation. As an illustration, the network of observation wells used to monitor the water table in the Equus Beds of Kansas is analyzed and an improved sampling pattern suggested. This example demonstrates the practical utility of the procedure, which can be applied equally well to other spatial sampling problems, as the procedure is not limited by the nature of the spatial function.  相似文献   

9.
The mean daily global solar radiation flux is influenced by astronomical, climatological, geographical, geometrical, meteorological, and physical parameters. This paper deals with the study of the effects of influencing parameters on the mean daily global solar radiation flux, and also with the computation of the solar radiation flux at the surface of the earth in locations without solar radiation measurements. The reference–real data were borrowed from the Iranian Meteorological Organization. The analysis of data showed that the mean daily solar radiation flux on a horizontal surface is related to parameters such as: mean daily extraterrestrial solar radiation, average daily ratio of sunshine duration, mean daily relative humidity, mean daily maximum air temperature, mean daily maximum dew point temperature, mean daily atmospheric pressure, and sine of the solar declination angle. Multiple regression and correlation analysis were applied to predict the mean daily global solar radiation flux on a horizontal surface. The models were validated when compared with the reference–measured data of global solar radiation flux. The results showed that the models estimate the global solar radiation flux within a narrow relative error band. The values of mean bias errors and root mean square errors were within acceptable margins. The predicted values of global solar radiation flux by this approach can be used for the design and performance estimation in solar applications. The model can be used in areas where meteorological stations do not exist and information on solar radiation flux cannot be obtained experimentally.  相似文献   

10.
 The applicability and usefulness of Geostatistics (kriging) as a tool for optimum selection of sites for monitoring groundwater levels has been demonstrated through a case study. The criterion used is the estimation of error variance. Groundwater level data (pre-monsoon 1994) obtained from 32 observation wells of Upper Kongal basin, Nalgonda District, A.P. (India) has been stochastically analyzed. The spatial distribution of water levels and its associated error variance is computed and the locations having maximum error variance are selected as additional sites for augmenting the existing observational well network. Received: 15 june 1998 · Accepted: 14 December 1998  相似文献   

11.
Stability is a key issue in any mining or tunnelling activity. Joint frequency constitutes an important input into stability analyses. Three techniques are used herein to quantify the local and spatial joint frequency uncertainty, or possible joint frequencies given joint frequency data, at unsampled locations. Rock quality designation is estimated from the predicted joint frequencies. The first method is based on kriging with subsequent Poisson sampling. The second method transforms the data to near-Gaussian variables and uses the turning band method to generate a range of possible joint frequencies. The third method assumes that the data are Poisson distributed and models the log-intensity of these data with a spatially smooth Gaussian prior distribution. Intensities are obtained and Poisson variables are generated to examine the expected joint frequency and associated variability. The joint frequency data is from an iron ore in the northern part of Norway. The methods are tested at unsampled locations and validated at sampled locations. All three methods perform quite well when predicting sampled points. The probability that the joint frequency exceeds 5 joints per metre is also estimated to illustrate a more realistic utilisation. The obtained probability map highlights zones in the ore where stability problems have occurred. It is therefore concluded that the methods work and that more emphasis should have been placed on these kinds of analyses when the mine was planned. By using simulation instead of estimation, it is possible to obtain a clear picture of possible joint frequency values or ranges, i.e. the uncertainty.  相似文献   

12.
An approach is proposed to predict the spatial distributions of undiscovered petroleum resources. Each pool is parameterized as a marked-point. The independence chain of the Hastings algorithm is used to generate an appropriate structure for pool combinations in a play. Petroleum-bearing favorability estimated from geological observations is used to represent the sampling probabilities of pool locations. An objective function measuring the distance between characteristics of the realization and constraints is constructed from both the pool size distribution and entropy maximum criterion, in which the entropy criterion places all undiscovered pools in the most favorable positions. The geometrical convergence property of the proposed Hastings algorithm is presented. The method is illustrated by a case study from the Western Canada Sedimentary Basin.  相似文献   

13.
Timeseries of estimated temperature have been combined to create global or hemispheric climate series over periods exceeding 1000 yr. The data used in these studies, however, may be subject to dating errors. It is shown that when timeseries with dating error are combined, the noise in the data smoothes periodic signals but leaves linear trends intact. This means that the effect of dating error of sample data in a timeseries reconstruction is to smooth out any signals (waves, cycles) that may be present. The purpose of this study was to develop signal extraction methods that will work for this type of historical data. The method used was nonlinear estimation of sample series where dating error has been added by Monte Carlo sampling. Several algorithms were tested for handling the dating error problem. Results were that using nonlinear model fitting, the periods of signals can be identified even from the averaged data. In a second stage of the estimation procedure, the cycle magnitudes can be estimated. Very good fits were achieved for two example cases. Temperature estimation error (white noise due to the use of proxies) was also considered and the method was extended to cover this case with quite good results. Using the new estimation methods, the information inherent in multiple series can be used to overcome the problem of dating error.  相似文献   

14.
Monte Carlo simulation was used to examine the error (statistical bias) introduced in estimating a sample semivariogram through application of oriented sampling patterns to variables which are correlated with fracture orientation. Sample semivariograms of the directional components of the water velocity were used to illustrate that oriented sampling schemes can provide biased data sets which result in error in the estimation of the semivariogram, particularly in the estimation of the sill (or variance). Three sampling patterns were used to analyze directional semivariograms of the components of the fluid velocity: sampling along lines parallel to the mean regional hydraulic gradient, sampling among lines perpendicular to the mean regional hydraulic gradient, and sampling along fracture segments. The first two sampling patterns were shown to introduce substantial error in the sills of the velocity variograms. It is argued that this error is due to the combination of unequal sampling of fractures with different orientations (i.e., sampling bias) and the systematic variation in the magnitude of the velocity components with orientation of the fracture. As a consequence, it is suggested that correction factors developed to correct fracture frequency statistics need to be extended to improve estimation of spatial moments of variables which are correlated with fracture orientation.  相似文献   

15.
This study proposes an interactive sampling strategy for locating the hot spot or maximum regions of a concerned attribute in a given area of survey. In the proposed strategy, information analysis is performed based on the ordinary kriging from the existing sample data to suggest a new batch of samples under the criterion of the highest information free energy. The information free energy (F) is a function of information energy (U) and information entropy (S) through F = U - TS, where T is information temperature and is used to coordinate the contribution of U and S to F. Information energy is the value of the concerned attribute, and information entropy is the transformed error variance of kriging and therefore measures the evenness and density of coverage of samples over the area under survey. At early sampling batches, information temperature is high and information entropy dominates the information free energy, and samples are suggested to give an even and dense enough coverage of the whole area under investigation. As samples accumulate, information temperature decreases to enlarge the contribution of information energy, and future samples are taken toward the locations with high attribute values. Two examples demonstrate the efficiency and effectiveness of the proposed sampling strategy in locating the hot spot regions of various fields: (1) a heavy metal contaminated site reproduced by modeling on 55 real field data; (2) a simulated two-dimensional field by the random phase volume (RPV) model. The results show that the proposed strategy, a robust interactive sampling procedure, is able to locate hot spot regions without compromising with the overall profile of an under-survey area.  相似文献   

16.
Drilling of a deep borehole does not normally allow for hydrologic testing during the drilling period. It is only done when drilling experiences a large loss (or high return) of drilling fluid due to penetration of a large-transmissivity zone. The paper proposes the possibility of conducting flowing fluid electrical conductivity (FFEC) logging during the drilling period, with negligible impact on the drilling schedule, yet providing important information on depth locations of both high- and low-transmissivity zones and their hydraulic properties. The information can be used to guide downhole fluid sampling and post-drilling detailed testing of the borehole. The method has been applied to the drilling of a 2,500-m borehole at Åre, central Sweden, firstly when the drilling reached 1,600 m, and then when the drilling reached the target depth of 2,500 m. Results unveil eight hydraulically active zones from 300 m down to borehole bottom, with depths determined to within the order of a meter. Further, the first set of data allows the estimation of hydraulic transmissivity values of the six hydraulically conductive zones found from 300 to 1,600 m, which are very low and range over one order of magnitude.  相似文献   

17.
This paper discusses the combination of kriging variances, which have been considered heretofor unfeasible since linearity of the problem and considerable simplifications which follow were overlooked. A simplified expression for global estimation variance is presented and an algorithm discussed with respect to precision and computer cost. A case study is presented, and, finally, an optimum calculation method is recommended.  相似文献   

18.
Accounting for Estimation Optimality Criteria in Simulated Annealing   总被引：1，自引：0，他引：1  
This paper presents both estimation and simulation as optimization problems that differ in the optimization criteria, minimization of a local expected loss for estimation and reproduction of global statistics (semivariogram, histogram) for simulation. An intermediate approach is proposed whereby an initial random image is gradually modified using simulated annealing so as to better match both local and global constraints. The relative weights of the different constraints in the objective function allow the user to strike a balance between smoothness of the estimated map and reproduction of spatial variability by simulated maps. The procedure is illustrated using a synthetic dataset. The proposed approach is shown to enhance the influence of observations on neighboring simulated values, hence the final realizations appear to be better conditioned to the sample information. It also produces maps that are more accurate (smaller prediction error) than stochastic simulation ignoring local constraints, but not as accurate as E-type estimation. Flow simulation results show that accounting for local constraints yields, on average, smaller errors in production forecast than a smooth estimated map or a simulated map that reproduces only the histogram and semivariogram. The approach thus reduces the risk associated with the use of a single realization for forecasting and planning.  相似文献   

19.
When sampling a random mixture of non-identical particles, the minimum sample size can be determined once the relationship between the sample size and the variance is known. Many different theories have been proposed to describe this relationship. This article provides a novel and systematic comparison of existing theories based on a set of objective criteria. It is found that none of the existing theories meet all criteria, making further research into the development of a suitable theory necessary.  相似文献   

20.
Modern marine sediments can be classified into terrigenous, hemipelagic, and pelagic types according to the distances from the land for the sites recovered and in terms of lithological characteristics of sediments. Ancient cherts are the best sedimentary rocks to reconstruct ancient depositional environments and hemipelagic and pelagic cherts can be identified in terms of lithologic and stratigraphic characteristics. However, geochemical characteristics of modern sediments and ancient cherts can discriminate them into these three types. Particularly, the degree of correlation between Fe and Mg in clay minerals is the most effective criterion for the classification, because the concentration of pelagic clay minerals in which the Fe/Mg ratio is constant may become indices for the classification into the three types. The correlation coefficient and goodness-of-fit to a certain regression line were used as indices of the classification and calculated for several sets of modern sediments and ancient cherts totalling to more than 400 samples. The correlation coefficients and the student-t values for the coefficients become better classification indices than the goodness-of-fit andF-values of the analysis of variance. Goodness-of-fit andF-values can be only used as subsidiary indices for the main classification indices of the correlation coefficients. The classification on the basis of the degree of correlation between Fe and Mg is effective and this approach can be used for other major chemical elements such as Al, K, and Na, if some sequential values of chemical analyses are available.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号