共查询到20条相似文献，搜索用时 390 毫秒
1.
Stochastic fractal (fGn and fBm) porosity and permeability fields are conditioned to given variogram, static (or hard), and multiwell pressure data within a Bayesian estimation framework. Because fGn distributions are normal/second-order stationary, it is shown that the Bayesian estimation methods based on the assumption of normal/second-order stationary distributions can be directly used to generate fGn porosity/permeability fields conditional to pressure data. However, because fBm is not second-order stationary, it is shown that such Bayesian estimation methods can be used with implementation of a pseudocovariance approach to generate fBm porosity/permeability fields conditional to multiwell pressure data. In addition, we provide methods to generate unconditional realizations of fBm/fGn fields honoring all variogram parameters. These unconditional realizations can then be conditioned to hard and pressure data observed at wells by using the randomized maximum likelihood method. Synthetic examples generated from one-, two-, and three-dimensional single-phase flow simulators are used to show the applicability of our methodology for generating realizations of fBm/fGn porosity and permeability fields conditioned to well-test pressure data and evaluating the uncertainty in reservoir performance predictions appropriately using these history-matched realizations. 相似文献
2.
Simulation of Non-Gaussian Transmissivity Fields Honoring Piezometric Data and Integrating Soft and Secondary Information 总被引：1，自引：0，他引：1
The conditional probabilities (CP) method implements a new procedure for the generation of transmissivity fields conditional to piezometric head data capable to sample nonmulti-Gaussian random functions and to integrate soft and secondary information. The CP method combines the advantages of the self-calibrated (SC) method with probability fields to circumvent some of the drawbacks of the SC method—namely, its difficulty to integrate soft and secondary information or to generate non-Gaussian fields. The SC method is based on the perturbation of a seed transmissivity field already conditional to transmissivity and secondary data, with the perturbation being function of the transmissivity variogram. The CP method is also based on the perturbation of a seed field; however, the perturbation is made function of the full transmissivity bivariate distribution and of the correlation to the secondary data. The two methods are applied to a sample of an exhaustive non-Gaussian data set of natural origin to demonstrate the interest of using a simulation method that is capable to model the spatial patterns of transmissivity variability beyond the variogram. A comparison of the probabilistic predictions of convective transport derived from a Monte Carlo exercise using both methods demonstrates the superiority of the CP method when the underlying spatial variability is non-Gaussian. 相似文献
3.
Chunsheng Ma 《Mathematical Geosciences》2011,43(2):229-242
This paper studies vector (multivariate, multiple, or multidimensional) random fields in space and/or time with second-order
increments, for which the variogram matrix is an important tool to measure the dependence within each component and between
each pair of distinct components. We introduce an efficient approach to construct Gaussian or non-Gaussian vector random fields
from the univariate random field with higher dimensional index domain, and particularly to generate a class of variogram matrices. 相似文献
4.
This paper describes two new approaches that can be used to compute the two-dimensional experimental wavelet variogram. They
are based on an extension from earlier work in one dimension. The methods are powerful 2D generalizations of the 1D variogram
that use one- and two-dimensional filters to remove different types of trend present in the data and to provide information
on the underlying variation simultaneously. In particular, the two-dimensional filtering method is effective in removing polynomial
trend with filters having a simple structure. These methods are tested with simulated fields and microrelief data, and generate
results similar to those of the ordinary method of moments variogram. Furthermore, from a filtering point of view, the variogram
can be viewed in terms of a convolution of the data with a filter, which is computed fast in O(NLogN) number of operations
in the frequency domain. We can also generate images of the filtered data corresponding to the nugget effect, sill and range
of the variogram. This in turn provides additional tools to analyze the data further. 相似文献
5.
The variogram matrix function is an important measure for the dependence of a vector random field with second-order increments, and is a useful tool for linear predication or cokriging. This paper proposes an efficient approach to construct variogram matrix functions, based on three ingredients: a univariate variogram, a conditionally negative definite matrix, and a Bernstein function, and derives three classes of variogram matrix functions for vector elliptically contoured random fields. Moreover, various dependence structures among components can be derived through appropriate mixture procedures demonstrated in this paper. We also obtain covariance matrix functions for second-order vector random fields through the Schoenberg–Lévy kernels. 相似文献
6.
Tingting Yao 《Mathematical Geology》2004,36(4):487-506
The application of spectral simulation is gaining acceptance because it honors the spatial distribution of petrophysical properties, such as reservoir porosity and shale volume. While it has been widely assumed that spectral simulation will reproduce the mean and variance of the important properties such as the observed net/gross ratio or global average of porosity, this paper shows the traditional way of implementing spectral simulation yields a mean and variance that deviates from the observed mean and variance. Some corrections (shift and rescale) could be applied to generate geologic models yielding the observed mean and variance; however, this correction implicitly rescales the input variogram model, so the variogram resulting from the generated cases has a higher sill than the input variogram model. Therefore, the spectral simulation algorithm cannot build geologic models honoring the desired mean, variance, and variogram model simultaneously, which is contrary to the widely accepted assumption that spectral simulation can reproduce all the target statistics. However, by using Fourier transform just once to generate values at all the cells instead of visiting each cell sequentially, spectral simulation does reproduce the observed variogram better than sequential Gaussian simulation. That is, the variograms calculated from the generated geologic models show smaller fluctuations around the target variogram. The larger the generated model size relative to the variogram range, the smaller the observed fluctuations. 相似文献
7.
8.
This paper is concerned with vector random fields on spheres with second-order increments, which are intrinsically stationary and mean square continuous and have isotropic variogram matrix functions. A characterization of the continuous and isotropic variogram matrix function on a sphere is derived, in terms of an infinite sum of the products of positive definite matrices and ultraspherical polynomials. It is valid for Gaussian or elliptically contoured vector random fields, but may not be valid for other non-Gaussian vector random fields on spheres such as a χ 2, log-Gaussian, or skew-Gaussian vector random field. Some parametric variogram matrix models are derived on spheres via different constructional approaches. A simulation study is conducted to illustrate the implementation of the proposed model in estimation and cokriging, whose performance is compared with that using the linear model of coregionalization. 相似文献
9.
Factorial Kriging (FK) is a data- dependent spatial filtering method that can be used to remove both independent and correlated
noise on geological images as well as to enhance lineaments for subsequent geological interpretation. The spatial variability
of signal, noise, and lineaments, characterized by a variogram model, have been used explicitly in calculating FK filter coefficients
that are equivalent to the kriging weighting coefficients. This is in contrast to the conventional spatial filtering method
by predefined, data-independent filters, such as Gaussian and Sobel filters. The geostatistically optimal FK filter coefficients,
however, do not guarantee an optimal filtering effect, if filter geometry (size and shape) are not properly selected. The
selection of filter geometry has been investigated by examining the sensitivity of the FK filter coefficients to changes in
filter size as well as variogram characteristics, such as nugget effect, type, range of influence, and anisotropy. The efficiency
of data-dependent FK filtering relative to data-independent spatial filters has been evaluated through simulated stochastic
images by two examples. In the first example, both FK and data-independent filters are used to remove white noise in simulated
images. FK filtering results in a less blurring effect than the data-independent fillers, even for a filter size as large
as 9 × 9. In the second example, FK and data-independent filters are compared relative to the extraction of lineaments and
components showing anisotropic variability. It was determined that square windows of the filter mask are effective only for
removing Isotropie components or white noise. A nonsquare windows must be used if anisotropic components are to be filtered
out. FK filtering for lineament enhancement is shown to be resistant to image noise, whereas data-independent filters are
sensitive to the presence of noise. We also have applied the FK filtering to the GLORIA side-scan sonar image from the Gulf
of Mexico, illustrating that FK is superior to the data-independent filters in removing noise and enhancing lineaments. The
case study also demonstrate that variogram analysis and FK filtering can be used for large images if a spectral analysis and
optimal filter design in the frequency domain is prohibitive because of a large memory requirement. 相似文献
10.
Markov chain Monte Carlo methods for conditioning a permeability field to pressure data 总被引：3，自引：0，他引：3
Generating one realization of a random permeability field that is consistent with observed pressure data and a known variogram
model is not a difficult problem. If, however, one wants to investigate the uncertainty of reservior behavior, one must generate
a large number of realizations and ensure that the distribution of realizations properly reflects the uncertainty in reservoir
properties. The most widely used method for conditioning permeability fields to production data has been the method of simulated
annealing, in which practitioners attempt to minimize the difference between the ’ ’true and simulated production data, and
“true” and simulated variograms. Unfortunately, the meaning of the resulting realization is not clear and the method can be
extremely slow. In this paper, we present an alternative approach to generating realizations that are conditional to pressure
data, focusing on the distribution of realizations and on the efficiency of the method. Under certain conditions that can
be verified easily, the Markov chain Monte Carlo method is known to produce states whose frequencies of appearance correspond
to a given probability distribution, so we use this method to generate the realizations. To make the method more efficient,
we perturb the states in such a way that the variogram is satisfied automatically and the pressure data are approximately
matched at every step. These perturbations make use of sensitivity coefficients calculated from the reservoir simulator. 相似文献
11.
用变差函数研究重磁场的区域变化特征.变差函数的变程反映重磁场的相干范围, 块金效应反映随机干扰, 基台值反映变异程度.重磁场的理论模拟说明: 重力场的相干范围大于磁场, 重磁场变程主要取决于场源深度, 浅源重磁场变差函数近似为球状模型或指数模型, 深源重磁场近似为连续性更好的高斯模型.磁场场源深度近似等于变程的一半, 重力场场源深度近似等于变程的四分之一.湖北大冶铁矿垂直分量磁异常具有几何各向异性, 北西-南东走向, 变差函数推测磁铁矿平均深度为250m.磁异常小波多尺度分解细节和逼近部分磁场具有协调几何各向异性, 变差函数的各阶场源深度估计结果与功率谱估计结果吻合. 相似文献
12.
Statistical assignment of upscaled flow functions for an ensemble of geological models 总被引：1，自引：0，他引：1
Upscaled flow functions are often needed to account for the effects of fine-scale permeability heterogeneity in coarse-scale
simulation models. We present procedures in which the required coarse-scale flow functions are statistically assigned to an
ensemble of upscaled geological models. This can be viewed as an extension and further development of a recently developed
ensemble level upscaling (EnLU) approach. The method aims to efficiently generate coarse-scale flow models capable of reproducing
the ensemble statistics (e.g., cumulative distribution function) of fine-scale flow predictions for multiple reservoir models.
The most expensive part of standard coarsening procedures is typically the generation of upscaled two-phase flow functions
(e.g., relative permeabilities). EnLU provides a means for efficiently generating these upscaled functions using stochastic
simulation. This involves the use of coarse-block attributes that are both fast to compute and correlate closely with the
upscaled two-phase functions. In this paper, improved attributes for use in EnLU, namely the coefficient of variation of the
fine-scale single-phase velocity field (computed during computation of upscaled absolute permeability) and the integral range
of the fine-scale permeability variogram, are identified. Geostatistical simulation methods, which account for spatial correlations
of the statistically generated upscaled functions, are also applied. The overall methodology thus enables the efficient generation
of coarse-scale flow models. The procedure is tested on 3D well-driven flow problems with different permeability distributions
and variable fluid mobility ratios. EnLU is shown to capture the ensemble statistics of fine-scale flow results (water and
oil flow rates as a function of time) with similar accuracy to full flow-based upscaling methods but with computational speedups
of more than an order of magnitude. 相似文献
13.
Digital simulation of multivariate two- and three-dimensional stochastic processes with a spectral turning bands method 总被引：3，自引：0，他引：3
Aristotelis Mantoglou 《Mathematical Geology》1987,19(2):129-149
The space domain version of the turning bands method can simulate multidimensional stochastic processes (random fields) having particular forms of covariance functions. To alleviate this limitation a spectral representation of the turning bands method in the two-dimensional case has shown that the spectral approach allows simulation of isotropic two-dimensional processes having any covariance or spectral density function. The present paper extends the spectral turning bands method (STBM) even further for simulation of much more general classes of multidimensional stochastic processes. Particular extensions include: (i) simulation of three-dimensional processes using STBM, (ii) simulation of anisotropic two- or three-dimensional stochastic processes, (iii) simulation of multivariate stochastic processes, and (iv) simulation of spatial averaged (integrated) processes. The turning bands method transforms the multidimensional simulation problem into a sum of a series of one-dimensional simulations. Explicit and simple expressions relating the cross-spectral density functions of the one-dimensional processes to the cross-spectral density function of the multidimensional process are derived. Using such expressions the one-dimensional processes can be simulated using a simple one-dimensional spectral method. Examples illustrating that the spectral turning bands method preserves the theoretical statistics are presented. The spectral turning bands method is inexpensive in terms of computer time compared to other multidimensional simulation methods. In fact, the cost of the turning bands method grows as the square root or the cubic root of the number of points simulated in the discretized random field, in the two- or three-dimensional case, respectively, whereas the cost of other multidimensional methods grows linearly with the number of simulated points. The spectral turning bands method currently is being used in hydrologic applications. This method is also applicable to other fields where multidimensional simulations are needed, e.g., mining, oil reservoir modeling, geophysics, remote sensing, etc. 相似文献
14.
Changik Han Jiyang Wang Mingguo Zheng Ende Wang Jianming Xia GwangSu Li Sunchol Choe 《Earth Science Informatics》2016,9(2):197-213
A critical step for kriging in geostatistics is estimation of the variogram. Traditional variogram modeling comprise of the experimental variogram calculation, appropriate variogram model selection and model parameter determination. Selecting of the variogram model and fitting of model parameters is the most controversial aspect of geostatistics. Shapes of valid variogram models are finite, and sometimes, the optimal shape of the model can not be fitted, leading to reduced estimation accuracy. In this paper, a new method is presented to automatically construct a model shape and fit model parameters to experimental variograms using Support Vector Regression (SVR) and Multi-Gene Genetic Programming (MGGP). The proposed method does not require the selection of a variogram model and can directly provide the model shape and parameters of the optimal variogram. The validity of the proposed method is demonstrated in a number of cases. 相似文献
15.
Estimating Variogram Uncertainty 总被引：10，自引：0，他引：10
The variogram is central to any geostatistical survey, but the precision of a variogram estimated from sample data by the method of moments is unknown. It is important to be able to quantify variogram uncertainty to ensure that the variogram estimate is sufficiently accurate for kriging. In previous studies theoretical expressions have been derived to approximate uncertainty in both estimates of the experimental variogram and fitted variogram models. These expressions rely upon various statistical assumptions about the data and are largely untested. They express variogram uncertainty as functions of the sampling positions and the underlying variogram. Thus the expressions can be used to design efficient sampling schemes for estimating a particular variogram. Extensive simulation tests show that for a Gaussian variable with a known variogram, the expression for the uncertainty of the experimental variogram estimate is accurate. In practice however, the variogram of the variable is unknown and the fitted variogram model must be used instead. For sampling schemes of 100 points or more this has only a small effect on the accuracy of the uncertainty estimate. The theoretical expressions for the uncertainty of fitted variogram models generally overestimate the precision of fitted parameters. The uncertainty of the fitted parameters can be determined more accurately by simulating multiple experimental variograms and fitting variogram models to these. The tests emphasize the importance of distinguishing between the variogram of the field being surveyed and the variogram of the random process which generated the field. These variograms are not necessarily identical. Most studies of variogram uncertainty describe the uncertainty associated with the variogram of the random process. Generally however, it is the variogram of the field being surveyed which is of interest. For intensive sampling schemes, estimates of the field variogram are significantly more precise than estimates of the random process variogram. It is important, when designing efficient sampling schemes or fitting variogram models, that the appropriate expression for variogram uncertainty is applied. 相似文献
16.
The multivariate (co)variogram as a spatial weighting function in classification methods 总被引：1，自引：0，他引：1
The multivariate variogram and the multivariate covariogram are used as spatial weighting functions for forming spatially homogeneous groups automatically. The groups are created after either deflating similarities between distant samples with the multivariate covariogram or by inflating dissimilarities between distant samples with the multivariate variogram. These approaches can be seen as generalization of the Oliver and Webster proposal. Two data sets show the efficiency of the two weighting functions when compared to the classical approach which does not take spatial information into account. In one case study, the weighting of similarities by the multivariate covariogram showed more interpretable results than the weighting of dissimilarities by the multivariate variogram. 相似文献
17.
18.
In the present paper, we propose a new method for the estimation of the variogram, which combines robustness with efficiency under intrinsic stationary geostatistical processes. The method starts by using a robust estimator to obtain discrete estimates of the variogram and control atypical observations that may exist. When the number of points used in the fit of a model is the same as the number of parameters, ordinary least squares and generalized least squares are asymptotically equivalent. Therefore, the next step is to fit the variogram by ordinary least squares, using just a few discrete estimates. The procedure is then repeated several times with different subsets of points and this produces a sequence of variogram estimates. The final estimate is the median of the multiple estimates of the variogram parameters. The suggested estimator will be called multiple variograms estimator. This procedure assures a global robust estimator, which is more efficient than other robust proposals. Under the assumed dependence structure, we prove that the multiple variograms estimator is consistent and asymptotically normally distributed. A simulation study confirms that the new method has several advantages when compared with other current methods. 相似文献
19.
基于饱和渗透系数空间变异结构的斜坡渗流及失稳特征 总被引：1，自引：0，他引：1
以往研究一般采用单随机变量方法（SRV）或基于水平或垂直方向波动范围生成的空间变异随机场来模拟岩土参数的空间变异性,对具有倾斜定向特征的空间变异随机场未有涉及.基于条件模拟相关理论和非侵入式随机有限元的理论框架,提出了利用序贯高斯模拟方法进行斜坡参数条件随机场模拟并运用有限元方法进行斜坡渗流和稳定性分析的方法.针对理想边坡,对各向同性和几何各向异性的共7种空间变异结构的饱和渗透系数（Ks）各进行了200次条件随机场模拟,基于条件随机场模拟结果进行了有限元渗流和稳定性计算,对每种空间变异结构多次计算结果进行了统计分析.结果表明:本文所提出的方法不仅再现了研究区域参数的空间二阶统计特性,通过设定变异函数参数进行不同空间变异类型、变异程度、变异定向性的随机场模拟,同时利用现场观测数据对随机场模拟结果进行条件限制,从而提高了随机场的赋值精度;Ks的空间变异结构对孔隙水压力的分布规律、地下水位线变化范围、稳定性系数和最危险滑动面分布特征均有一定程度的影响.本研究为库岸斜坡稳定性评价提供方法支撑. 相似文献
20.
3D geological modeling should be an effective tool for accurately representing the geometric structure boundary and internal property fields of geological body. However, conventional methods have rarely focused on the expression of geological heterogeneous properties, known as 4D modeling. A volume function is defined as a piecewise mathematic function that describes a parameterized geological property field and created by fitting the property functions of certain 3D voxels. The quadratic generalized tri-prism volume function (QGTPVF) model is proposed for representing geological property fields with a sedimentary strata structure based on the volume function method and generalized tri-prism (GTP) voxel model. A QGTPVF is designed for borehole sample points, and it interpolates geological properties by fitting quadratic volume functions combined with the influence of geological geometric structure constraints, including stratum interfaces and fault planes. This research mainly focuses on the QGTPVF definition and fitting method, single GTP volume function is a continuous quadratic function, and a property smoothing along directionally adjacent bedding GTPs method is also studied, so geological property field is also expressed by the smooth function. A preprocessing method of faults geological structures is given, and the general framework QGTPVF is also discussed for practical applications. The QGTPVF model could be converted to tetrahedron voxel and 3D grid models for calculation and analysis easily. An example on a porosity property field is studied to verify the method’s accuracy and reliability by comparing to Kriging and inverse distance weighted interpolation methods, and the result of accuracy and time complexity are better, so the QGTPVF provides a new solution for modeling geological property fields. 相似文献