Direct Sequential Co-simulation with Joint Probability Distributions
|
|
- Nathaniel Jenkins
- 6 years ago
- Views:
Transcription
1 Math Geosci (2010) 42: DOI /s x Direct Sequential Co-simulation with Joint Probability Distributions Ana Horta Amílcar Soares Received: 13 May 2009 / Accepted: 3 January 2010 / Published online: 5 February 2010 International Association for Mathematical Geosciences 2010 Abstract The practice of stochastic simulation for different environmental and earth sciences applications creates new theoretical problems that motivate the improvement of existing algorithms. In this context, we present the implementation of a new version of the direct sequential co-simulation (Co-DSS) algorithm. This new approach, titled Co-DSS with joint probability distributions, intends to solve the problem of mismatch between co-simulation results and experimental data, i.e. when the final biplot of simulated values does not respect the experimental relation known for the original data values. This situation occurs mostly in the beginning of the simulation process. To solve this issue, the new co-simulation algorithm, applied to a pair of covariates Z 1 (x) and Z 2 (x), proposes to resample Z 2 (x) from the joint distribution F(z 1,z 2 ) or, more precisely, from the conditional distribution of Z 2 (x 0 ), at a location x 0, given the previously simulated value z (l) 1 (x 0) (F(Z 2 Z 1 = z (l) 1 (x 0)). The work developed demonstrates that Co-DSS with joint probability distributions reproduces the experimental bivariate cdf and, consequently, the conditional distributions, even when the correlation coefficient between the covariates is low. Keywords Geostatistics Joint sequential simulation Bivariate distribution 1 Introduction The use of geostatistical models based on stochastic simulation algorithms is a reliable option for addressing problems in environmental and earth sciences if the purpose is to assess the spatial distribution of a certain attribute as well as spatial uncertainty. A. Horta ( ) A. Soares Center for Natural Resources and Environment, Instituto Superior Técnico, Technical University of Lisbon, Lisbon, Lisbon, Portugal ahorta@ist.utl.pt
2 270 Math Geosci (2010) 42: To reproduce the spatial relation of more than one interdependent attribute or variable, joint simulation algorithms have been proposed for different research fields: environment (Goovaerts 1997), petroleum (Caers 2005; Deutsch 2002), and mining (Journel and Kyriakidis 2004). To reproduce joint spatial patterns (given by covariograms and bi-histograms), the following simulation algorithms can be used: sequential multi-gaussian co-simulation (Verly 1993), multi-gaussian co-simulation with collocated co-kriging (Almeida and Journel 1994), co-simulation with LU decomposition method (Myers 1988), and simulation of autocorrelation factors (Desbarats and Dimitrakopoulos 2000; Boucher and Dimitrakopoulos 2009). However, for most practical applications, the common stochastic co-simulation algorithm applied to a set of correlated variables is based on a sequential approach where each variable is simulated in turn with the other previously simulated variables, as secondary information (Almeida and Journel 1994). This is the principle for the sequential Gaussian co-simulation (Almeida and Journel 1994) and the direct sequential co-simulation (Soares 2001). The main advantages of using the sequential approach are its simplicity, efficiency, and the possibility of choosing which variable will condition the next covariate. Normally, this variable (called the primary) is chosen due to its evident spatial continuity or given its relevance for the case study (Almeida and Journel 1994). Let us consider the direct co-simulation (Soares 2001) of just two variables, Z 1 (x) and Z 2 (x), where Z 1 (x) is the primary variable. This joint simulation algorithm can be summarized in two steps. First, the simulation of the value z (l) 1 (x u) (where x u is a node in a regular grid) using the direct sequential simulation (DSS) algorithm (Soares 2001). Second, assuming the previously simulated Z 1 (x) as the secondary variable, the DSS algorithm is again applied to simulate Z 2 (x). Collocated simple co-kriging is used to calculate z 2 (x u ) and σsk 2 (x u) conditioned to nearby data z 2 (x α ) and the collocated datum z 1 (x u ) (Goovaerts 1997). By using direct sequential co-simulation (Co-DSS), we perform the joint simulation of a set of N v variables without any prior or subsequent transformation. Also, the method ensures the reproduction of marginal cdfs (F 1 (z 1 ) and F 2 (z 2 )), variograms (γ 1 (h) and γ 2 (h)) and the joint spatial pattern characterized by the covariograms γ 1,2 (h). The use of Co-DSS in different case studies revealed that, for some situations, the method did not reproduce the bivariate cdf, i.e. the experimental bi-histogram was not respected. As an example of this practical fact, see the case study presented by Soares (2005). Figure 1 shows the experimental scatterplot of two variables, copper (Z 1 (x)) and tin (Z 2 (x)), from a sulfide copper deposit. Figure 2 presents the resulting scatterplot of co-simulated values. By comparing real and simulated biplots, we notice that high simulated values of Z 2 (x) are collocated to low values of Z 1 (x), unlike the experimental data plot. Nevertheless, the joint simulation method ensures the reproduction of the variograms of Z 1 (x) and Z 2 (X), the marginal histograms and the correlation coefficient. Hence, for a particular class of values of Z 1 (x) = z, the simulated values z (l) 2 (x) must not exceed the higher and lower limits of the experimental data class to which z (l) 1 (x) belongs, and should reproduce the experimental conditional histogram of Z 2 (x), givenz 1 (x) = z. This can be considered a limitation of the existing sequential co-simulation algorithm (including sequential Gaussian co-simulation Co-SGS and direct sequential co-simulation). For those methods, after estimating the conditional mean and variance by simple kriging and
3 Math Geosci (2010) 42: Fig. 1 Experimental biplot (Sn vs. Cu) Fig. 2 Simulation biplot using usual Co-DSS (co-simulated Sn vs. simulated Cu)
4 272 Math Geosci (2010) 42: simple collocated co-kriging, z (l) 1 (x) and z(l) 2 (x) are drawn, respectively, from its local Gaussian transform (Co-SGS) or from the marginal global distributions (Co-DSS). Although the marginal histograms of Z 1 (x) and Z 2 (x) are reproduced (as required for both methods), high conditional variances, common at the very beginning of the sequential simulation process, can produce simulated values that do not respect the conditional histograms. This can lead to erroneous conclusions, mostly when nonlinear cost functions (used in environmental and mining applications) are applied to the set of co-simulated values. In this paper, we propose a new approach for Co-DSS based on the idea of resampling Z 2 (x) from the joint distribution F(z 1,z 2 ) or, more precisely, from the conditional distribution of Z 2 (x 0 ), at a location x 0, given the previously simulated value z (l) 1 (x 0) (F(Z 2 Z 1 = z (l) 1 (x 0)). This method succeeds in reproducing the experimental bivariate cdf and, consequently, the conditional distributions. 2 Direct Sequential Co-simulation with Joint Probability Distributions 2.1 Methodology The problem posed by the implementation of the Co-DSS algorithm consists in the fact that the final bi-distribution of simulated and co-simulated values does not reproduce the same experimental distribution of the original data values. Hence a conditional distributions F [Z 2 (x) Z 1 (x)] are not necessarily reproduced at the end of the simulation. The direct sequential co-simulation procedure is based on the application of the Bayes rule in a successive sequence of steps. A first variable Z 1 (x), chosen as the most important for the case study or with the most continuous spatial pattern, is simulated for the entire area. The second variable Z 2 (x) is then co-simulated taking into account the experimental data z 2 (x α ) and previously simulated values z (l) 1 (x) as secondary data. At a given location x 0 of the sequential simulation path, we estimate the local mean and variance using collocated co-kriging, conditioned to the experimental data z 2 (x α ) and previously simulated data z (l) 1 (x 0). Based on the local mean and variance, the simulated value (z (l) 2 (x 0)) is resampled from the global cdf F [Z 2 (x)]. If the local estimated variance (σ 2 [Z 2 (x 0 ) Z 2 (x α ), Z (l) 1 (x 0)]) is sufficiently high in the beginning of simulation process, most probably the values z (l) 2 (x 0), drawn from the global cdf F [Z 2 (x)], will be outside of the boundaries determined by the conditional distribution F [Z 2 (x) Z 1 (x)]. Consequently, the bi-histograms of simulated values (z (l) 2 (x), z(l) 1 (x)) will have higher conditional variances than the corresponding estimates with experimental data, with the same correlation coefficient and marginal histograms. The solution proposed in this paper is based on the use of conditional distributions (estimated from experimental data) to simulate Z 2 (x), given the known Z 1 (x): Z 2 (x) is simulated from F [Z 2 (x) Z 1 (x) = Z (l) 1 (x)], and the known realiza- (x) as secondary data, the conditional tion z (l) 1 (x) at each location x. Assuming Z(l) 1
5 Math Geosci (2010) 42: mean is computed for every location x 0, identified with the simple collocated cokriging estimator (Z 2 (x 0 ) ) [ z2 (x 0 ) ] sck = α λ α,x0 [ z2 (x α ) m 2 ] + λx0 [ z1 (x 0 ) m 1 ] + m2 (1) where λ is the co-kriging weights, z 2 (x α ) is the experimental conditioning data, z 1 (x 0 ) is the collocated datum, and m 1 and m 2 are the means of z 1 and z 2 data, respectively. Also, the conditional variance corresponds to the simple collocated cokriging variance. The value z (l) 2 (x 0) is drawn from the conditional distribution F [ Z 2 (x 0 ) Z1 (x 0 ) = z1 l (x 0) ] = prob { Z 2 (x 0 )<z Z1 (x 0 ) = z (l) 1 (x 0) } K (2) rather than from the global cdf of Z 2 (x) as in the usual Co-DSS. The proposed algorithm continues to ensure that simulated realizations of Z 1 (x) and Z 2 (x) reproduce the variograms (γ 1 (h) and γ 2 (h)), the marginal cdfs (F(Z 1 (x) and F(Z 2 (x)) and the conditional cdf (F(Z 2 (x) Z 1 (x)). The proposed methodology of Co-DSS with joint probability distributions can be summarized in the following sequence, beginning with the estimation of the global bi-distribution from experimental data, followed by the simulation of the first covariate Z 1 (x) using direct sequential simulation. Realizations of Z (l) 1 (x) reproduce the variogram γ 1 (h) and marginal cdf F(Z 1 (x)). Then the co-simulation of Z 2 (x) is calculated at each location x 0, and an estimation of the local mean and variance is made, identified with estimated simple collocated co-kriging and corresponding estimation variance. Then, based on previously simulated z (l) 1 (x 0), computation of the conditional cdf F [Z 2 (x) Z 1 (x) = z (l) 1 (x 0)] from the bi-distribution F [Z 1 (x), Z 2 (x)] occurs. Finally, resampling of the simulated value z (l) 2 (x 0) from the conditional cdf F [Z 2 (x) Z 1 (x) = z (l) 1 (x 0)] as in the usual direct sequential simulation procedure (Soares 2001) is made. 2.2 Corrections for Local Bias and Local Means The corrections for local bias and local means applied in DSS and Co-DSS (Soares 2001) are now applied in the framework of resampling from the conditional histogram. When resampling a histogram F [Z 2 (x) Z 1 (x) = z (l) 1 (x 0)], with little data, mostly in low frequency classes, the following bias can occur: E{z (l) 2 (x 0)} z 2 (x 0 ), where z 2 (x 0 ) is obtained by (1) (see Soares 2001 for details). Therefore, in each sequential step, the following estimator of E{z (l) (x 0 )} is calculated: Z (l) 2 (x 0) = 1 N i Ni i=1 z(l,i) 2 (x 0 ), by drawing, for example, N i values of z (l,i) 2 (x 0 ) at each iteration. The simulated value z (l) (x 0 ) is corrected by the deviation z 2 (x 0 ) z (l) 2 (x 0) z (l) 2 (x 0) = z (l) 2 (x 0) + [ z 2 (x 0 ) z (l) 2 (x 0) ] (3) Note the fact that, in each sequential step, N i values (say N i = 50) are drawn to calculate z (l) 2 (x 0), which does not have a significant impact in terms of computing
6 274 Math Geosci (2010) 42: time. Regarding the correction of local means, suppose that m z is the conditional mean of F(Z 2 (x) Z 1 (x)) to be reproduced in the final simulated maps. A deviation e s can be calculated between m z and the mean of simulated values z (l) 2 (x) at a given step s of the sequential procedure where m (l) z is the mean of all simulated values at step s m (l) z e s = m z m (l) z (4) = 1 N s N s i=1 z (l) 2 (x i) (5) N s being the number of simulated nodes up to step s. Hence, the estimated local means (1) can be corrected with the deviation e s The remaining steps of Co-DSS are unchanged. Z (l) (x 0 ) = z (l) (x 0 ) + e s (6) 2.3 Computing Conditional Distributions: Practical Implementation Based on previously simulated z (l) 1 (x 0), the conditional cdf F [Z 2 (x) Z 1 (x) = z (l) 1 (x 0)] is calculated from the bi-distribution F [Z 1 (x), Z 2 (x)] using moving classes (rather than fixed classes) of Z 2 (x). The user must define the minimum number of data N c of each class of the conditional histogram F [Z 2 (x) Z 1 (x) = z (l) 1 (x 0)]. By ranking the pairs of values (Z 1 (x i ), Z 2 (x i )) by increasing order of Z 1 (x i ), for one given conditioning value z (l) 1 (x 0) = Z 1 (x j ) we have the corresponding Z 2 (x j ),inthe position jth of ordered (Z 1 (x i ), Z 2 (x i )). Hence, F [Z 2 (x) Z 1 (x) = z (l) 1 (x 0)] is composed of the closest N c values of Z 2 (x j ) in the rank-ordered list. Figure 3 shows this procedure schematically. 3 Case Study 3.1 Data Analysis To illustrate the proposed methodology, we chose a case study of a sulfide copper deposit (Neves Corvo, Portugal). Copper (Cu) and tin (Sn) contents were analyzed in 523 boreholes (marginal histograms are presented in Figs. 4 and 5). The Cu/Sn biplot in Fig. 1 shows the existing relationship between both elements. Spatial continuity main patterns of Cu and Sn can both present a similar isotropic behavior, modeled by an exponential model. Sn shows a clear nugget effect, representing 20% of the total variance. As noted in Soares (2005), the biplot in Fig. 1 shows two populations with a different behavior regarding the correlation between Cu and Sn. Splitting the total set of samples using a Cu threshold of 10%, the values with a Cu content lower than 10% show a higher correlation coefficient (r = 0.72; Fig. 6) than the values with a Cu content higher than 10%, which do not present a significant correlation with Sn (r = 0.37; Fig. 7). This set of data was used to test both Co-DSS methods.
7 Math Geosci (2010) 42: Fig. 3 Computing moving classes to define the conditional distribution: for a given z 1 (x 0 ), the conditional distribution F(Z 2 (x) Z 1 (x)) is calculated with the N c selected values of Z 2 (x) 3.2 Co-simulation with Joint Probability Distributions of Cu and Sn Cu and Sn grades were simulated in a regular grid of points ( m). First, Cu was simulated using DSS. Then, Sn values were obtained by Co-DSS assuming previously simulated maps of Cu as a secondary variable. The two algorithms were tested, with and without the improvement proposed to account for joint distributions. Figure 8 presents the result (biplot of co-simulated Sn and simulated Cu) using the proposed Co-DSS with joint probability distributions. Figure 9 shows the result using the usual Co-DSS. The comparison with Fig. 1 (experimental biplot) reveals that the worst matching between simulated Cu and Sn is shown in Fig. 9, when the usual Co-DSS is used. In this case, low value classes of Cu are related with high simulated values of Sn, contradicting the experimental bi-histogram. Also, the marginal histograms obtained for Sn (Figs. 10 and 11) were compared with the experimental histogram (Fig. 5). These results show that the histogram is better reproduced when we use Co-DSS with joint probability distributions.
8 276 Math Geosci (2010) 42: Fig. 4 Cu marginal histogram Fig. 5 Sn marginal histogram
9 Math Geosci (2010) 42: Fig. 6 Experimental biplot (Sn vs. Cu with content lower than 10%) Fig. 7 Experimental biplot (Sn vs. Cu with content higher than 10%)
10 278 Math Geosci (2010) 42: Fig. 8 Simulation biplot using Co-DSS with joint probability distributions (co-simulated Sn vs. simulated Cu) Fig. 9 Simulation biplot using usual Co-DSS (co-simulated Sn vs. simulated Cu)
11 Math Geosci (2010) 42: Fig. 10 Marginal histogram of Sn co-simulated values using Co-DSS with joint probability distributions Fig. 11 Marginal histogram of Sn co-simulated values using usual Co-DSS
12 280 Math Geosci (2010) 42: Another test was done using the two subsets of Cu samples, one corresponding to the classes of Cu values lower than 10% (Fig. 6) and the other to the classes of Cu values equal to or greater than 10% (Fig. 7). Co-simulation results are plotted in Fig. 12(a) and (b) (Co-DSS with joint probability distributions) as well as Fig. 13(a) and (b) (Co-DSS). When compared with the experimental biplot (Figs. 6 and 7), we conclude that even for different correlation coefficients, Co-DSS with joint probability distributions succeeds in better reproducing the experimental bivariate relation. Also, conditional histograms were calculated corresponding to F(Sn Cu < 5%), F(Sn Cu < 7.5%), F(Sn Cu < 10%), F(Sn 5% < Cu < 10%) and F(Sn 10% < Cu < 15%) (see Fig. 14(a c), Fig. 15(a c), Fig. 16(a c), Fig. 17(a c), and Fig. 18(a c)). For these examples, Sn conditional histograms obtained with Co- DSS with joint probability distributions and with the usual Co-DSS were compared with the conditional histogram obtained from the experimental data. The evidence is that conditional histograms are better reproduced with the improved Co-DSS algorithm. It is worth noting that, for low classes of Cu values (less than 5%), there is an excessive frequency of zero values of Sn that are not observed in the experimental conditional histogram. Finally, we compare, in Fig. 19(a) and (b), co-simulated images obtained by both methods. Using low threshold of Sn values (for example, 0.2%) (Fig. 19(c) and (d)) it is noticeable that Co-DSS tends to overestimate low Sn values since it is not conditioned to the experimental bi-histogram. 4 Conclusions The work presented in this paper contributes to the improvement of the Co-DSS algorithm to ensure that simulated covariates reproduce the bivariate distributions as measured or estimated by the experimental data. This was not yet implemented in Co-DSS or any sequential co-simulation algorithm. The advantages of the proposed approach are related with the importance of conditional distributions of simulated covariates for several reasons, such as environmental applications, where the definition of critical areas in contaminated sites, based on the joint concentration of multi-elements, is provided by the product of impact costs of different pollutant concentrations above a certain remediation level. Therefore, the correct characterization of conditional distributions is extremely important, mostly when impact costs are non-linear functions of contaminant contents and consequently can magnify the existing bias of the conditional distributions. This approach can also benefit earth science applications, where the classification of mining reserves of multi-element ore bodies must be based on the joint distributions, mostly when treatment and commercial costs are non-linear functions of the grades and risk of magnifying the costs of misclassification errors. Concerning the practical implementation of the new algorithm, there are no significant differences between both algorithms in what concerns CPU efficiency, although for a large data set, Co-DSS with joint probability distributions could have a slower performance. Other practical issues can be pointed out as limitations in the use of Co-DSS with joint probability distributions. It can be used when little experimental data is available but is advisable to use an optimization algorithm of smoothing
13 Math Geosci (2010) 42: (a) (b) Fig. 12 Simulation biplot using Co-DSS with joint probability distributions: (a) co-simulated Sn vs. simulated Cu with content lower than 10%; (b) co-simulated Sn vs. simulated Cu with content higher than 10%
14 282 Math Geosci (2010) 42: (a) (b) Fig. 13 Simulation biplot using usual Co-DSS: (a) co-simulated Sn vs. simulated Cu with content lower than 10%; (b) co-simulated Sn vs. simulated Cu with content higher than 10%
15 Math Geosci (2010) 42: (a) (b) Fig. 14 Conditional histogram of Sn given Cu content less than 5%: (a) experimental data; (b) Co-DSS with joint probability distributions; (c)co-dss
16 284 Math Geosci (2010) 42: (c) Fig. 14 (Continued) (a) Fig. 15 Conditional histogram of Sn given Cu content less than 7.5%: (a) experimental data; (b)co-dss with joint probability distributions; (c)co-dss
17 Math Geosci (2010) 42: (b) Fig. 15 (Continued) (c)
18 286 Math Geosci (2010) 42: (a) (b) Fig. 16 Conditional histogram of Sn given Cu content less than 10%: (a) experimental data; (b) Co-DSS with joint probability distributions; (c)co-dss
19 Math Geosci (2010) 42: (c) Fig. 16 (Continued) (a) Fig. 17 Conditional histogram of Sn given Cu content between 5% and 10%: (a) experimental data; (b) Co-DSS with joint probability distributions; (c)co-dss
20 288 Math Geosci (2010) 42: (b) Fig. 17 (Continued) (c)
21 Math Geosci (2010) 42: (a) (b) Fig. 18 Conditional histogram of Sn given Cu content between 10% and 15%: (a) experimental data; (b) Co-DSS with joint probability distributions; (c)co-dss
22 290 Math Geosci (2010) 42: (c) Fig. 18 (Continued) (a) Fig. 19 Co-simulated Sn: (a) Co-DSS with joint probability distributions; (b) Co-DSS; (c) Snvalues lower than 0.2% using Co-DSS with joint probability distributions; (d) Sn values lower than 0.2% using Co-DSS
23 Math Geosci (2010) 42: (b) (c) Fig. 19 (Continued) bi-histograms, for example as proposed by Deutsch and Journel (1998). The purpose is to obtain a more representative (or less erratic fluctuations) bi-histogram. Also, the proposed method is based on a sequential approach where each variable is simulated in turn having the others, previously simulated, as secondary information: F [Z 1,Z 2,...,Z N ]=F [Z 1 ] F [Z 2 Z 1 ]...Then, the joint simulation of multielements reproduces the conditional distributions for the pair of covariates which are being co-simulated. The case study presented refers to the particular case for the joint simulation of only two variables. The influence of each covariate in each cosimulation was not evaluated in this work.
24 292 Math Geosci (2010) 42: (d) Fig. 19 (Continued) Acknowledgements The company Somincor Sociedade Mineira de Neves Corvo is gratefully acknowledged for giving permission to publish the data. This paper was produced in the context of Project Soil Contamination Risk Assessment (PTDC/CTE-SPA/69127/2006) (financed by FEDER through the National Operational Science and Innovation Program 2010 with the support of the Foundation for Science and Technology). References Almeida A, Journel A (1994) Joint simulation of multiple variables with a Markov-type coregionalization model. Math Geol 26(5): Boucher A, Dimitrakopoulos R (2009) Block simulation of multiple correlated variables. Math Geosci 41: Caers J (2005) Petroleum geostatistics. Society of Petroleum Engineers, Richardson Desbarats AJ, Dimitrakopoulos R (2000) Geostatistical simulation of regionalized pore-size distributions using min/max autocorrelation factors. Math Geol 32(8): Deutsch CV (2002) Geostatistical reservoir modeling. Oxford University Press, New York Deutsch CV, Journel AG (1998) Gslib: geostatistical software library and user s guide. Oxford University Press, New York Goovaerts P (1997) Geostatistics for natural resources evaluation. Oxford University Press, New York Journel AG, Kyriakidis P (2004) Evaluation of mineral reserves a simulation approach. Oxford University Press, New York Myers DE (1988) Vector conditional simulation. In: Armstrong M (ed) Geostatistics. Kluwer, Dordrecht, pp Soares A (2001) Direct sequential simulation and cosimulation. Math Geol 33(8): Soares A (2005) Classification of mining reserves using direct sequential simulation. In: Leuangthong O Deutsch CV (eds) Geostatistics Banff 2004, vol 1. Springer, Dordrecht, pp Verly GW (1993) Sequential Gaussian cosimulation: A simulation method integrating several types of information. In: Soares A (ed) Geostatistics Troia, vol 1. Kluwer, Dordrecht, pp
Geostatistical Reservoir Characterization of McMurray Formation by 2-D Modeling
Geostatistical Reservoir Characterization of McMurray Formation by 2-D Modeling Weishan Ren, Oy Leuangthong and Clayton V. Deutsch Department of Civil & Environmental Engineering, University of Alberta
More informationSelected Implementation Issues with Sequential Gaussian Simulation
Selected Implementation Issues with Sequential Gaussian Simulation Abstract Stefan Zanon (szanon@ualberta.ca) and Oy Leuangthong (oy@ualberta.ca) Department of Civil & Environmental Engineering University
More informationUsing Blast Data to infer Training Images for MPS Simulation of Continuous Variables
Paper 34, CCG Annual Report 14, 212 ( 212) Using Blast Data to infer Training Images for MPS Simulation of Continuous Variables Hai T. Nguyen and Jeff B. Boisvert Multiple-point simulation (MPS) methods
More informationIntegration of Geostatistical Modeling with History Matching: Global and Regional Perturbation
Integration of Geostatistical Modeling with History Matching: Global and Regional Perturbation Oliveira, Gonçalo Soares Soares, Amílcar Oliveira (CERENA/IST) Schiozer, Denis José (UNISIM/UNICAMP) Introduction
More information11-Geostatistical Methods for Seismic Inversion. Amílcar Soares CERENA-IST
11-Geostatistical Methods for Seismic Inversion Amílcar Soares CERENA-IST asoares@ist.utl.pt 01 - Introduction Seismic and Log Scale Seismic Data Recap: basic concepts Acoustic Impedance Velocity X Density
More informationMultiple Point Statistics with Multiple Training Images
Multiple Point Statistics with Multiple Training Images Daniel A. Silva and Clayton V. Deutsch Abstract Characterization of complex geological features and patterns has been one of the main tasks of geostatistics.
More informationShort Note: Some Implementation Aspects of Multiple-Point Simulation
Short Note: Some Implementation Aspects of Multiple-Point Simulation Steven Lyster 1, Clayton V. Deutsch 1, and Julián M. Ortiz 2 1 Department of Civil & Environmental Engineering University of Alberta
More informationLocal Recoverable Reserves Prediction with Block LU Simulation
Local Recoverable Reserves Prediction with Block LU Simulation Julián Ortiz C. 1 and Clayton V. Deutsch 1 Department of Mining Engineering University of Chile Centre for Computational Geostatistics (CCG)
More informationIndicator Simulation Accounting for Multiple-Point Statistics
Indicator Simulation Accounting for Multiple-Point Statistics Julián M. Ortiz 1 and Clayton V. Deutsch 2 Geostatistical simulation aims at reproducing the variability of the real underlying phenomena.
More informationTensor Based Approaches for LVA Field Inference
Tensor Based Approaches for LVA Field Inference Maksuda Lillah and Jeff Boisvert The importance of locally varying anisotropy (LVA) in model construction can be significant; however, it is often ignored
More informationThe Proportional Effect: What it is and how do we model it?
The Proportional Effect: What it is and how do we model it? John Manchuk Centre for Computational Geostatistics Department of Civil & Environmental Engineering University of Alberta One of the outstanding
More informationSimulating Geological Structures Based on Training Images and Pattern Classifications
Simulating Geological Structures Based on Training Images and Pattern Classifications P. Switzer, T. Zhang, A. Journel Department of Geological and Environmental Sciences Stanford University CA, 9435,
More informationA workflow to account for uncertainty in well-log data in 3D geostatistical reservoir modeling
A workflow to account for uncertainty in well-log data in 3D geostatistical reservoir Jose Akamine and Jef Caers May, 2007 Stanford Center for Reservoir Forecasting Abstract Traditionally well log data
More informationA Geostatistical and Flow Simulation Study on a Real Training Image
A Geostatistical and Flow Simulation Study on a Real Training Image Weishan Ren (wren@ualberta.ca) Department of Civil & Environmental Engineering, University of Alberta Abstract A 12 cm by 18 cm slab
More informationSampling informative/complex a priori probability distributions using Gibbs sampling assisted by sequential simulation
Sampling informative/complex a priori probability distributions using Gibbs sampling assisted by sequential simulation Thomas Mejer Hansen, Klaus Mosegaard, and Knud Skou Cordua 1 1 Center for Energy Resources
More informationA Program for Robust Calculation of Drillhole Spacing in Three Dimensions
A Program for Robust Calculation of Drillhole Spacing in Three Dimensions David F. Machuca Mory and Clayton V. Deutsch Centre for Computational Geostatistics (CCG) Department of Civil & Environmental Engineering
More informationVariogram Inversion and Uncertainty Using Dynamic Data. Simultaneouos Inversion with Variogram Updating
Variogram Inversion and Uncertainty Using Dynamic Data Z. A. Reza (zreza@ualberta.ca) and C. V. Deutsch (cdeutsch@civil.ualberta.ca) Department of Civil & Environmental Engineering, University of Alberta
More informationMarkov Bayes Simulation for Structural Uncertainty Estimation
P - 200 Markov Bayes Simulation for Structural Uncertainty Estimation Samik Sil*, Sanjay Srinivasan and Mrinal K Sen. University of Texas at Austin, samiksil@gmail.com Summary Reservoir models are built
More informationGeostatistics on Stratigraphic Grid
Geostatistics on Stratigraphic Grid Antoine Bertoncello 1, Jef Caers 1, Pierre Biver 2 and Guillaume Caumon 3. 1 ERE department / Stanford University, Stanford CA USA; 2 Total CSTJF, Pau France; 3 CRPG-CNRS
More informationPetrel TIPS&TRICKS from SCM
Petrel TIPS&TRICKS from SCM Knowledge Worth Sharing Petrel 2D Grid Algorithms and the Data Types they are Commonly Used With Specific algorithms in Petrel s Make/Edit Surfaces process are more appropriate
More informationGeostatistical modelling of offshore diamond deposits
Geostatistical modelling of offshore diamond deposits JEF CAERS AND LUC ROMBOUTS STANFORD UNIVERSITY, Department of Petroleum Engineering, Stanford, CA 94305-2220, USA jef@pangea.stanford.edu TERRACONSULT,
More informationsurface but these local maxima may not be optimal to the objective function. In this paper, we propose a combination of heuristic methods: first, addi
MetaHeuristics for a Non-Linear Spatial Sampling Problem Eric M. Delmelle Department of Geography and Earth Sciences University of North Carolina at Charlotte eric.delmelle@uncc.edu 1 Introduction In spatial
More informationFlexible Lag Definition for Experimental Variogram Calculation
Flexible Lag Definition for Experimental Variogram Calculation Yupeng Li and Miguel Cuba The inference of the experimental variogram in geostatistics commonly relies on the method of moments approach.
More informationNon Stationary Variograms Based on Continuously Varying Weights
Non Stationary Variograms Based on Continuously Varying Weights David F. Machuca-Mory and Clayton V. Deutsch Centre for Computational Geostatistics Department of Civil & Environmental Engineering University
More informationAccounting for Estimation Optimality Criteria in Simulated Annealing 1
Mathematical Geology, Vol. 30, No. 5, 1998 Accounting for Estimation Optimality Criteria in Simulated Annealing 1 P. Goovaerts 2 This paper presents both estimation and simulation as optimization problems
More informationIndicator Simulation for Categorical Variables
Reservoir Modeling with GSLIB Indicator Simulation for Categorical Variables Sequential Simulation: the Concept Steps in Sequential Simulation SISIM Program Sequential Simulation: the Concept 2 1 3 1.
More informationConditioning a hybrid geostatistical model to wells and seismic data
Conditioning a hybrid geostatistical model to wells and seismic data Antoine Bertoncello, Gregoire Mariethoz, Tao Sun and Jef Caers ABSTRACT Hybrid geostatistical models imitate a sequence of depositional
More informationAlgorithm-driven and Representation-driven Random Function : A New Formalism for Applied Geostatistics
Algorithm-driven and Representation-driven Random Function : A New Formalism for Applied Geostatistics Alexandre Boucher Dept of Geological and Environmental Sciences Stanford University Abstract This
More informationApplication of MPS Simulation with Multiple Training Image (MultiTI-MPS) to the Red Dog Deposit
Application of MPS Simulation with Multiple Training Image (MultiTI-MPS) to the Red Dog Deposit Daniel A. Silva and Clayton V. Deutsch A Multiple Point Statistics simulation based on the mixing of two
More informationPrograms for MDE Modeling and Conditional Distribution Calculation
Programs for MDE Modeling and Conditional Distribution Calculation Sahyun Hong and Clayton V. Deutsch Improved numerical reservoir models are constructed when all available diverse data sources are accounted
More informationModeling Multiple Rock Types with Distance Functions: Methodology and Software
Modeling Multiple Rock Types with Distance Functions: Methodology and Software Daniel A. Silva and Clayton V. Deutsch The sub division of the deposit into estimation domains that are internally consistent
More informationCONDITIONING FACIES SIMULATIONS WITH CONNECTIVITY DATA
CONDITIONING FACIES SIMULATIONS WITH CONNECTIVITY DATA PHILIPPE RENARD (1) and JEF CAERS (2) (1) Centre for Hydrogeology, University of Neuchâtel, Switzerland (2) Stanford Center for Reservoir Forecasting,
More information2D Geostatistical Modeling and Volume Estimation of an Important Part of Western Onland Oil Field, India.
and Volume Estimation of an Important Part of Western Onland Oil Field, India Summary Satyajit Mondal*, Liendon Ziete, and B.S.Bisht ( GEOPIC, ONGC) M.A.Z.Mallik (E&D, Directorate, ONGC) Email: mondal_satyajit@ongc.co.in
More informationThe Components of Geostatistical Simulation
This file was created by scanning the printed publication. Errors identified by the software have been corrected; however, some errors may remain. The Components of Geostatistical Simulation Carol A. ~otwayl
More informationMPS Simulation with a Gibbs Sampler Algorithm
MPS Simulation with a Gibbs Sampler Algorithm Steve Lyster and Clayton V. Deutsch Complex geologic structure cannot be captured and reproduced by variogram-based geostatistical methods. Multiple-point
More informationSpatial Interpolation & Geostatistics
(Z i Z j ) 2 / 2 Spatial Interpolation & Geostatistics Lag Lag Mean Distance between pairs of points 1 Tobler s Law All places are related, but nearby places are related more than distant places Corollary:
More informationCONDITIONAL SIMULATION OF TRUNCATED RANDOM FIELDS USING GRADIENT METHODS
CONDITIONAL SIMULATION OF TRUNCATED RANDOM FIELDS USING GRADIENT METHODS Introduction Ning Liu and Dean S. Oliver University of Oklahoma, Norman, Oklahoma, USA; ning@ou.edu The problem of estimating the
More informationSpatial Interpolation - Geostatistics 4/3/2018
Spatial Interpolation - Geostatistics 4/3/201 (Z i Z j ) 2 / 2 Spatial Interpolation & Geostatistics Lag Distance between pairs of points Lag Mean Tobler s Law All places are related, but nearby places
More informationSurface Smoothing Using Kriging
1 AutoCAD Civil 3D has the ability to add data points to a surface based on mathematical criteria. This gives you the ability to strengthen the surface definition in areas where data may be sparse or where
More informationSpatial Interpolation & Geostatistics
(Z i Z j ) 2 / 2 Spatial Interpolation & Geostatistics Lag Lag Mean Distance between pairs of points 11/3/2016 GEO327G/386G, UT Austin 1 Tobler s Law All places are related, but nearby places are related
More informationUSE OF DIRECT AND INDIRECT DISTRIBUTIONS OF SELECTIVE MINING UNITS FOR ESTIMATION OF RECOVERABLE RESOURCE/RESERVES FOR NEW MINING PROJECTS
USE OF DIRECT AND INDIRECT DISTRIBUTIONS OF SELECTIVE MINING UNITS FOR ESTIMATION OF RECOVERABLE RESOURCE/RESERVES FOR NEW MINING PROJECTS W Assibey-Bonsu 1 and D G Krige 2 1 Gold Fields Limited, South
More informationCrosswell Tomographic Inversion with Block Kriging
Crosswell Tomographic Inversion with Block Kriging Yongshe Liu Stanford Center for Reservoir Forecasting Petroleum Engineering Department Stanford University May, Abstract Crosswell tomographic data can
More informationUniversity of Alberta. Multivariate Analysis of Diverse Data for Improved Geostatistical Reservoir Modeling
University of Alberta Multivariate Analysis of Diverse Data for Improved Geostatistical Reservoir Modeling by Sahyun Hong A thesis submitted to the Faculty of Graduate Studies and Research in partial fulfillment
More informationChuanrong Zhang 1 and Weidong Li Department of Geography, Kent State University, Kent, Ohio 44242
Comparing a Fixed-Path Markov Chain Geostatistical Algorithm with Sequential Indicator Simulation in Categorical Variable Simulation from Regular Samples Chuanrong Zhang 1 and Weidong Li Department of
More informationGeostatistical Time-Lapse Seismic Inversion
Geostatistical Time-Lapse Seismic Inversion Valter Proença Abstract With the decrease of the oil price and of new discoveries, it is not only important to use new exploration and production methods and
More informationMulti-Objective Stochastic Optimization by Co-Direct Sequential Simulation for History Matching of Oil Reservoirs
Multi-Objective Stochastic Optimization by Co-Direct Sequential Simulation for History Matching of Oil Reservoirs João Daniel Trigo Pereira Carneiro under the supervision of Amílcar de Oliveira Soares
More informationSPE Copyright 2002, Society of Petroleum Engineers Inc.
SPE 77958 Reservoir Modelling With Neural Networks And Geostatistics: A Case Study From The Lower Tertiary Of The Shengli Oilfield, East China L. Wang, S. Tyson, Geo Visual Systems Australia Pty Ltd, X.
More informationD025 Geostatistical Stochastic Elastic Iinversion - An Efficient Method for Integrating Seismic and Well Data Constraints
D025 Geostatistical Stochastic Elastic Iinversion - An Efficient Method for Integrating Seismic and Well Data Constraints P.R. Williamson (Total E&P USA Inc.), A.J. Cherrett* (Total SA) & R. Bornard (CGGVeritas)
More informationHigh Resolution Geomodeling, Ranking and Flow Simulation at SAGD Pad Scale
High Resolution Geomodeling, Ranking and Flow Simulation at SAGD Pad Scale Chad T. Neufeld, Clayton V. Deutsch, C. Palmgren and T. B. Boyle Increasing computer power and improved reservoir simulation software
More informationGeostatistics 2D GMS 7.0 TUTORIALS. 1 Introduction. 1.1 Contents
GMS 7.0 TUTORIALS 1 Introduction Two-dimensional geostatistics (interpolation) can be performed in GMS using the 2D Scatter Point module. The module is used to interpolate from sets of 2D scatter points
More informationUncertainty Quantification for History- Matching of Non-Stationary Models Using Geostatistical Algorithms
Uncertainty Quantification for History- Matching of Non-Stationary Models Using Geostatistical Algorithms Maria Helena Caeiro 1, Vasily Demyanov 2, Mike Christie 3 and Amilcar Soares 4 Abstract The quantification
More informationImprovement of Realizations through Ranking for Oil Reservoir Performance Prediction
Improvement of Realizations through Ranking for Oil Reservoir Performance Prediction Stefan Zanon, Fuenglarb Zabel, and Clayton V. Deutsch Centre for Computational Geostatistics (CCG) Department of Civil
More informationGeostatistics on Unstructured Grids: Application to Teapot Dome
Geostatistics on Unstructured Grids: Application to Teapot Dome John G. Manchuk This work demonstrates some of the developments made in using geostatistical techniques to populate unstructured grids with
More informationMultivariate Standard Normal Transformation
Multivariate Standard Normal Transformation Clayton V. Deutsch Transforming K regionalized variables with complex multivariate relationships to K independent multivariate standard normal variables is an
More informationVARFIT: A Program for Semi-Automatic Variogram Modelling
VARFIT: A Program for Semi-Automatic Variogram Modelling Paula F. Larrondo (larrondo@ualberta.ca), Chad T. Neufeld (neufeld@ualberta.ca), and Clayton V. Deutsch (cdeutsch@ualberta.ca) Department of Civil
More informationTrend Modeling Techniques and Guidelines
Trend Modeling Techniques and Guidelines Jason A. M c Lennan, Oy Leuangthong, and Clayton V. Deutsch Centre for Computational Geostatistics (CCG) Department of Civil and Environmental Engineering University
More informationLocal Linear Approximation for Kernel Methods: The Railway Kernel
Local Linear Approximation for Kernel Methods: The Railway Kernel Alberto Muñoz 1,JavierGonzález 1, and Isaac Martín de Diego 1 University Carlos III de Madrid, c/ Madrid 16, 890 Getafe, Spain {alberto.munoz,
More informationCLASS VS. THRESHOLD INDICATOR SIMULATION
CHAPTER 4 CLASS VS. THRESHOLD INDICATOR SIMULATION With traditional discrete multiple indicator conditional simulation, semivariogram models are based on the spatial variance of data above and below selected
More informationOn internal consistency, conditioning and models of uncertainty
On internal consistency, conditioning and models of uncertainty Jef Caers, Stanford University Abstract Recent research has been tending towards building models of uncertainty of the Earth, not just building
More informationMultiple-point geostatistics: a quantitative vehicle for integrating geologic analogs into multiple reservoir models
Multiple-point geostatistics: a quantitative vehicle for integrating geologic analogs into multiple reservoir models JEF CAERS AND TUANFENG ZHANG Stanford University, Stanford Center for Reservoir Forecasting
More informationLinear interpolation and joint model fitting of experimental transiograms for Markov chain simulation of categorical spatial variables
International Journal of Geographical Information Science Vol. 24, No. 6, June 2010, 821 839 Linear interpolation and joint model fitting of experimental transiograms for Markov chain simulation of categorical
More informationGridding and Contouring in Geochemistry for ArcGIS
Gridding and Contouring in Geochemistry for ArcGIS The Geochemsitry for ArcGIS extension includes three gridding options: Minimum Curvature Gridding, Kriging and a new Inverse Distance Weighting algorithm.
More informationB. Todd Hoffman and Jef Caers Stanford University, California, USA
Sequential Simulation under local non-linear constraints: Application to history matching B. Todd Hoffman and Jef Caers Stanford University, California, USA Introduction Sequential simulation has emerged
More informationFast FILTERSIM Simulation with Score-based Distance Function
Fast FILTERSIM Simulation with Score-based Distance Function Jianbing Wu (1), André G. Journel (1) and Tuanfeng Zhang (2) (1) Department of Energy Resources Engineering, Stanford, CA (2) Schlumberger Doll
More informationUncertainty analysis in spatial environmental modelling. Kasia Sawicka Athens, 29 Mar 2017
Uncertainty analysis in spatial environmental modelling Kasia Sawicka Athens, 29 Mar 2017 Uncertainty propagation overview Input ± U Model (± U) Output ± U Input data Model parameters Model structure Solution
More informationB002 DeliveryMassager - Propagating Seismic Inversion Information into Reservoir Flow Models
B2 DeliveryMassager - Propagating Seismic Inversion Information into Reservoir Flow Models J. Gunning* (CSIRO Petroleum) & M.E. Glinsky (BHP Billiton Petroleum) SUMMARY We present a new open-source program
More informationFMA901F: Machine Learning Lecture 3: Linear Models for Regression. Cristian Sminchisescu
FMA901F: Machine Learning Lecture 3: Linear Models for Regression Cristian Sminchisescu Machine Learning: Frequentist vs. Bayesian In the frequentist setting, we seek a fixed parameter (vector), with value(s)
More informationSPE Copyright 1999, Society of Petroleum Engineers Inc.
SPE 56654 Methodology for Variogram Interpretation and Modeling for Improved Reservoir Characterization E. Gringarten, SPE, Landmark Graphics Corp., and C. V. Deutsch, SPE, University of Alberta Copyright
More informationImprovements in Continuous Variable Simulation with Multiple Point Statistics
Improvements in Continuous Variable Simulation with Multiple Point Statistics Jeff B. Boisvert A modified version of Mariethoz et al s (2010) algorithm for simulating continuous variables using multiple
More informationRotation and affinity invariance in multiple-point geostatistics
Rotation and ainity invariance in multiple-point geostatistics Tuanfeng Zhang December, 2001 Abstract Multiple-point stochastic simulation of facies spatial distribution requires a training image depicting
More informationMethods to define confidence intervals for kriged values: Application on Precision Viticulture data.
Methods to define confidence intervals for kriged values: Application on Precision Viticulture data. J-N. Paoli 1, B. Tisseyre 1, O. Strauss 2, J-M. Roger 3 1 Agricultural Engineering University of Montpellier,
More informationExploring Direct Sampling and Iterative Spatial Resampling in History Matching
Exploring Direct Sampling and Iterative Spatial Resampling in History Matching Matz Haugen, Grergoire Mariethoz and Tapan Mukerji Department of Energy Resources Engineering Stanford University Abstract
More informationkeywords: geostatistical simulation; scenario reduction; machine learning; net present value (NPV).
http://dx.doi.org/10.1590/0370-44672018720135 Ryu Okada 1,2 http://orcid.org/0000-0001-9629-9358 João Felipe Coimbra Leite Costa 1,3 Áttila Leães Rodrigues 1,4 http://orcid.org/0000-0003-4524-4087 Bruno
More informationThe Effect of Changing Grid Size in the Creation of Laser Scanner Digital Surface Models
The Effect of Changing Grid Size in the Creation of Laser Scanner Digital Surface Models Smith, S.L 1, Holland, D.A 1, and Longley, P.A 2 1 Research & Innovation, Ordnance Survey, Romsey Road, Southampton,
More informationA Generalized Markov Chain Approach for Conditional Simulation of Categorical Variables from Grid Samples
: 651 669 Research Article A Generalized Markov Chain Approach for Conditional Simulation of Categorical Variables from Grid Samples Weidong Li Department of Geography Kent State University Chuanrong Zhang
More informationResampling Methods for Dependent Data
S.N. Lahiri Resampling Methods for Dependent Data With 25 Illustrations Springer Contents 1 Scope of Resampling Methods for Dependent Data 1 1.1 The Bootstrap Principle 1 1.2 Examples 7 1.3 Concluding
More informationOn Secondary Data Integration
On Secondary Data Integration Sahyun Hong and Clayton V. Deutsch A longstanding problem in geostatistics is the integration of multiple secondary data in the construction of high resolution models. In
More informationMultivariate Standard Normal Transform: Advances and Case Studies
Multivariate Standard Normal Transform: Advances and Case Studies Ryan M. Barnett and Clayton V. Deutsch The Multivariate Standard Normal Transformation (MSNT) was recently proposed to transform arbitrary
More informationFacies Modeling Using a Markov Mesh Model Specification
Math Geosci (2011) 43:611 624 DOI 10.1007/s11004-011-9350-9 Facies Modeling Using a Markov Mesh Model Specification Marita Stien Odd Kolbjørnsen Received: 28 April 2010 / Accepted: 13 February 2011 / Published
More informationSOME PRACTICAL COMPUTATIONAL ASPECTS OF MINE PLANNING
M. Guarascio et al (eds.), Advanced Geostatistics in the Mining Industry, 391-399. All Rights Reserved. Copyright 1976 by D. Reidel Publishing Company, Dordrecht-Holland. SOME PRACTICAL COMPUTATIONAL ASPECTS
More informationHIERARCHICAL SIMULATION OF MULTIPLE-FACIES RESERVOIRS USING MULTIPLE-POINT GEOSTATISTICS
HIERARCHICAL SIMULATION OF MULTIPLE-FACIES RESERVOIRS USING MULTIPLE-POINT GEOSTATISTICS AREPORT SUBMITTED TO THE DEPARTMENT OF PETROLEUM ENGINEERING OF STANFORD UNIVERSITY IN PARTIAL FULFILLMENT OF THE
More informationCross-validation. Cross-validation is a resampling method.
Cross-validation Cross-validation is a resampling method. It refits a model of interest to samples formed from the training set, in order to obtain additional information about the fitted model. For example,
More informationSurface-based model conditioning using an hybrid optimization: methodology and application
Surface-based model conditioning using an hybrid optimization: methodology and application Antoine Bertoncello, Jef Caers, Hongmei Li and Tao Sun Department of Energy Resources Engineering Stanford University
More informationCOPULA MODELS FOR BIG DATA USING DATA SHUFFLING
COPULA MODELS FOR BIG DATA USING DATA SHUFFLING Krish Muralidhar, Rathindra Sarathy Department of Marketing & Supply Chain Management, Price College of Business, University of Oklahoma, Norman OK 73019
More informationDirect Pattern-based Simulation of Non-Stationary Geostatistical Models
Direct Pattern-based Simulation of Non-Stationary Geostatistical Models Mehrdad Honarkhah and Jef Caers Department of Energy Resources Engineering Stanford University Abstract Non-stationary models often
More informationReservoir Characterization with Limited Sample Data using Geostatistics
Reservoir Characterization with Limited Sample Data using Geostatistics By: Sayyed Mojtaba Ghoraishy Submitted to the Department of Chemical and Petroleum Engineering and the Faculty of the Graduate School
More informationFeature Selection. Department Biosysteme Karsten Borgwardt Data Mining Course Basel Fall Semester / 262
Feature Selection Department Biosysteme Karsten Borgwardt Data Mining Course Basel Fall Semester 2016 239 / 262 What is Feature Selection? Department Biosysteme Karsten Borgwardt Data Mining Course Basel
More informationThe Bootstrap and Jackknife
The Bootstrap and Jackknife Summer 2017 Summer Institutes 249 Bootstrap & Jackknife Motivation In scientific research Interest often focuses upon the estimation of some unknown parameter, θ. The parameter
More informationSoftware Tutorial Session Universal Kriging
Software Tutorial Session Universal Kriging The example session with PG2000 which is described in this and Part 1 is intended as an example run to familiarise the user with the package. This documented
More informationNon-linear Univariate and. Optimal Design
Non-linear Univariate and Multivariate Spatial Modelling and Optimal Design Gnai Nishani Musafer B.Sc (Hons) Principal supervisor: Dr. Helen Thompson Associate supervisors: Prof. Erhan Kozan, Prof. Rodney
More information1. Estimation equations for strip transect sampling, using notation consistent with that used to
Web-based Supplementary Materials for Line Transect Methods for Plant Surveys by S.T. Buckland, D.L. Borchers, A. Johnston, P.A. Henrys and T.A. Marques Web Appendix A. Introduction In this on-line appendix,
More informationAppropriate algorithm method for Petrophysical properties to construct 3D modeling for Mishrif formation in Amara oil field Jawad K.
Appropriate algorithm method for Petrophysical properties to construct 3D modeling for Mishrif formation in Amara oil field Jawad K. Radhy AlBahadily Department of geology, college of science, Baghdad
More informationarxiv: v1 [eess.iv] 25 Dec 2018
Vector Field-based Simulation of Tree-Like Non-Stationary Geostatistical Models arxiv:1812.11030v1 [eess.iv] 25 Dec 2018 Viviana Lorena Vargas, Sinesio Pesco Pontificia Universidade Catolica, Rio de Janeiro
More informationODSIM: An Object-Distance Simulation. method for Conditioning Complex Natural. Structures 1
ODSIM: An Object-Distance Simulation method for Conditioning Complex Natural Structures 1 Vincent Henrion 23 4 Guillaume Caumon 35 Nicolas Cherpeau 36 December 22, 2011 Submitted to Mathematical Geosciences
More informationA 3D code for mp simulation of continuous and
A 3D code for mp simulation of continuous and categorical variables: FILTERSIM Jianbing Wu, Alexandre Boucher & André G. Journel May, 2006 Abstract In most petroleum and geoscience studies, the flow is
More informationImage Segmentation for Image Object Extraction
Image Segmentation for Image Object Extraction Rohit Kamble, Keshav Kaul # Computer Department, Vishwakarma Institute of Information Technology, Pune kamble.rohit@hotmail.com, kaul.keshav@gmail.com ABSTRACT
More information2 P. ABRAHAMSEN Reservoir, Z R 6 Reference, (Z = 0) Top Reservoir, Z TR (reflector) Layer 3, Z L3 Top Layer 2, Z TL2 Layer 2, Z L2 Top Layer 1, Z TL1
COMBINING METHODS FOR SUBSURFACE PREDICTION PETTER ABRAHAMSEN Norwegian Computing Center, PO Bo 114 Blindern, NO-0314 Oslo, Norway Abstract. The depth to subsurfaces in a multi-layer model is obtained
More informationSuper-Resolution Land Cover Mapping with Indicator Geostatistics
Super-Resolution Land Cover Mapping with Indicator Geostatistics Alexandre Boucher 1 and Phaedon C. Kyriakidis 2 1 Department of Geological and Environmental Sciences Stanford University 2 Department of
More informationINFERENCE OF THE BOOLEAN MODEL ON A NON STATIONARY CASE
8 th annual conference of the International Association for Mathematical Geology, 15-20 September 2002, Berlin INFERENCE OF THE BOOLEAN MODEL ON A NON STATIONARY CASE M. BENITO GARCÍA-MORALES and H. BEUCHER
More informationCollaborative Rough Clustering
Collaborative Rough Clustering Sushmita Mitra, Haider Banka, and Witold Pedrycz Machine Intelligence Unit, Indian Statistical Institute, Kolkata, India {sushmita, hbanka r}@isical.ac.in Dept. of Electrical
More information