Rotation and affinity invariance in multiple-point geostatistics
|
|
- Marion McCoy
- 5 years ago
- Views:
Transcription
1 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 the geometric patterns deemed present in the study volume. That training image is the equivalent of the variogram model used in traditional 2-point statistics-based simulation. Therefore the training image should be reasonably stationary, that is invariant by any translation. However, even within a homogeneous formation, actual geological patterns may change gradually in space, e.g., in terms of direction and aspect ratio. Such location-dependent patterns can be simulated using location-dependent linear transforms (rotation + ainity) of the stationary patterns read from the training image. The linear transform of training patterns is discussed in this report and its application to simulation of non-stationary fan deposits is shown. 1 Introduction Training images as source of structured data are taking an important role in geostatistics. Under the useful decision of stationarity, multiple-point statistics can be lifted from a training image and used to model a target area. The training image should reflect the patterns of spatial variability deemed to prevail over the target area. These patterns must be repeated often enough over the training image so that they can be captured by multiple-point statistics. Statistics, no matter their order and the number of locations (multiple points) they involve, can only represent stationary features or patterns; they reflect average patterns not location-specific unique features. Geostatistics, as any other spatial statistics, can only capture stationary patterns from a training image and can only model and simulate these stationary patterns. The specific locations and local transforms 1
2 of these stationary patterns can only be obtained through local data: this is the task of conditioning and the essential dierence between non-conditional simulation which can only reproduce stationary patterns and conditional simulation which, in addition, also aims at locating them in space. Any unique (non repeated) pattern, such as a single fracture, is either detected by local data and can be deterministically reproduced, or it is not and no statistics could ever model it. A training image, that is a representation of repeated spatial patterns deemed present over the target area, need not and should not carry any location-specific characteristic. Training patterns are defined relative to each other and their statistics should be independent of the coordinate system used. More precisely, they should be independent of any linear transform of the coordinates. Training images can be adapted to the target area by a series of local rotations and ainity of the coordinates system, much like an isotropic variogram model can be made location-specific by a succession of rotations and ainity transforms of the coordinates, Xu (1996). Consider the training image (TI) of Figure 1a which features EW elongated s with little varying NS width. This TI can be used to model s of any direction, see Figures 1 b-c-d, and of any width, see Figures 2 and 3. The training patterns taken from the TI are made rotation and scale-invariant, more precisely independent of any linear transform of the coordinates system. The impact of stationarity of statistics on stochastic simulation is now investigated. 2 Multiple-point statistics Let fs(u); u 2 D)g be the stationary random function modeling the spatial distribution of an attribute s(u) over an reservoir D. IfS(u) has a finite number of outcomes: s 1 ;s 2 ; :::; s K and a conditional data set fs(u ); =1; 2; :::; ng is available, the probability of the variable S(u) at any unsampled location u can be obtained by considering the correlation between the single-point event S(u) and the known multiple-point (mp) data event B = fs(u ); = 1; 2; :::; ng. The outcome of the multiple-data event B conditions the probability of the single-point event S(u). The mp data is considered as a whole instead of being considered one datum location at a time. Such mp data allows capturing patterns within the data set. This is dierent from the traditional two-point geostatistics approach, in which only the correlation between any two points is considered to construct the probability of S(u). Traditional geostatistics only captures two-point correlation characteristics of the reservoir attributes. In contrast, multiple-point statistics can model curvilinear and large scale structures such as undulating s. Consider the binary random variable I B associated with occurrence of the data 2
3 event B: I B = ( 1 if B occurs 0 if not I B is called an event indicator variable. If the event B involves a single point u in space, I B is the ordinary indicator variable usually denoted by I(u). IfB is a multiple-point event, I B is then an indicator of the specific pattern constituted by all single data constituting B. Consider now the two events A k = fs(u) = s k g and B = fs(u ); = 1; 2; :::; ng. Simple kriging can be used to estimate the probability of the unknown event A k given that data event B has occurred: P (A k jb) =E(I Ak =1ji B )=E(I Ak )+ [1 E(B)] (1) A single normal equation provides the weight : V arfbg = CovfA k ;Bg (2) It is easy to check that equations (1) and (2) exactly identify Bayes relation, Strebelle (2000): P (A k jb) = P (A k;b) (3) P (B) The equivalent relations (1) and (3) show that inference of the conditional probability P (A k jb) requires knowing the probability of the joint occurrence of the two events A k and B as well as that of B alone. In the usual situation of sparse data, it is not possible to infer directly from actual data these two probabilities which are much beyond variograms and any other two-point statistics. An alternative is to identify these probabilities to proportions read from an exhaustively sampled training image, assuming that this training image has the same multi-point statistics as the reservoir under study. An easy and fast way to store and then retrieve the proportions of occurrence of any training pattern is provided by the concept of search tree Strebelle (2000). The training image is scanned only once with a specified template and the corresponding proportions of pattern occurrence are stored. Once the search tree is created, the probabilities attached to all patterns collected from the training image within the specified template can be easily retrieved. 3
4 3 Search Tree The search tree is a powerful tool to store the number of replicates of any mp data event (i.e. pattern) found over the training image within the template retained. Related mp probabilities are then identified to these training proportions and retrieved. The building of a search tree proceeds as follows: l Specify the size and geometry of the template and order the corresponding grid nodes according to their distance to the center of the template. An anisotropic variogram distance can be used. Consider a template fi with n ordered nodes: fi = fh ; =1; 2; :::; ng where h = u u is the vector linking the center u of the template to the -th node location u. l Scan the training image for all occurrences of any data event within the previously specified template and any of its n sub-templates defined as: fi k = fh i ; i =1; 2; :::; Kg; K =1; 2; :::; n The smallest sub-template fi 0 corresponds to a single datum location at the center u of template fi. The next sub-template fi 1 includes that central location u and the next closest node u 1. Finally: fi 0 ρ fi 1 ρ fi 2 ρ ::: ρ fi n = fi l Starting from the smallest fi 0, these sub-templates are centered at each node of the training image: the corresponding training data event is recorded and a proper count is incremented. A data event consists of the sub-template geometry and data values: d k = fs(u + h i )=s i ; i =1; 2; :::; kg Once scanning of the training image is completed, a search tree has been generated. That tree stores the numbers of replicates of any data event found within the specified data template fi The search tree is divided into n +1levels, each level i corresponding to the sub-template fi i ;i = 1; 2; :::; n. The root of the search tree (level i = 0) stores the number of single-point occurrences of each category in the training image, that is the global proportion of each category. This is a single-point statistics. Level 1 of the 4
5 search tree records the number of occurrences of each category at the second node of the sub-template fi 1. From these first two levels of the search tree, one can derive the probability of occurrence of any data event within the sub-template fi 1. This is a twopoint statistics. Similarly, one can extract three-point, four-point,..., up to n-point statistics. Note that not all potential n-point data configurations may exist in any specific training image, in such cases some branches of the search tree will stop before the (n +1)th level. This feature significantly reduces the size of the search tree and RAM demand. In a simulation using multi-point statistics, we need to infer the probability of a specific event conditioned by specific data values arranged along a specific geometric template. Denote this conditioning data event by: d = fs(u + h j )=s j jj =1; 2; :::; mg In most situations, this specific data template fi m does not identify any one of the sub-templates fi i used to generate the search tree. In such case, retain the largest training sub-template fi i included in the conditioning data template fi m, with fi i ρ fi m. The probability associated to the conditioning data event d is then approximated by that corresponding to that largest training sub-template fi i. 4 Rotation and ainity transforms The rotation and ainity (squeezing) transform of a training image is now discussed. These linear coordinates transforms are associated with the rotation and scale invariance of multiple-point statistics. At first, consider a transform with constant rotation angle and constant ainity factors. Next, a transform with location-dependent rotation angles and ainity factors is considered. Application is made to a fan deposit simulated from a training image with straight s. All examples are 2D, but the algorithm is easily extended to 3D. All simulations are non-conditional to focus strictly on reproduction of the mp statistics. 4.1 Rotation Let ft (u); u 2 D L g be the training image defined over the grid D L. Corresponding unconditional simulations should display a single main direction of continuity but possibly dierent from that of the training image. Let be the desired direction of continuity measured in degrees clockwise from the continuity direction (East) of the training image. The simulated image is denoted by fs (u)ju 2 Dg. The necessary rotation can be implemented as follow: 5
6 l Keep the original search tree. Rotate any conditioning data event by - within the data template prior to looking for the number of replicates in the search tree. The coordinates rotation is written: where u new and u old u new = R u old u =(x; y; z) (4) are the old and new coordinates of each single datum component in a data event. R is the rotation matrix. In 2D with u =(x; y), this is written: ( x new y new = x old = x old Because the rotated location (x new cos yold sin sin + yold cos (5) ;ynew ) does not usually fall on a grid node of the training image, it is relocated to the nearest grid node. Note that this relocation entails changing somewhat the original data geometry, see Appendix. A training image (Figure 1a) and a square data template with 7 7 = 49 grid nodes are used. Only two facies are considered: in black with training proportion 0.3 and in white with proportion 0.7. Three rotation angles =0 i ; 45 i, and 90 i are considered. The corresponding unconditional realizations of size are displayed in Figures 1b to 1d. The realizations do reproduce the desired direction of maximum continuity. The lack of full continuity across the simulated area is not related to the rotations done. It is a consequence of using a template of limited size (7 7) 4.2 Ainity We can also squeeze the given training image (an ainity transform) to change the aspect ratio of its geometric patterns. In a fluvial sedimentary environment, this amounts to change the width and sinuosity of the training image s. Two ainity factors a x, a y along the two desired major and minor directions of continuity control the magnitude of the ainity. If the elongation in the training image is horizontal (x-direction), then a y > 1 leads to thinner s and a x > 1 leads to more sinuous s.the corresponding coordinates transform is written: u new = Au old u =(x; y; z) (6) where u new and u old are the old and new coordinates of each single datum location of the mp data event. A is the ainity matrix: 6
7 A = In 2D with u =(x; y), this is written: 0 a x a y a z x new y new = a x x old = a y y old Using the training image of Figure 1a, the following ainity factors were considered: a x =1.0, a y = 2.0, 3.0, 4.0, 5.0. The corresponding four realizations are displayed in Figure 2. The simulated s are increasingly thinner as a y increases. The four realizations corresponding to a x =1.0, a y = 0.9, 0.8, 0.7, 0.6, are displayed in Figure 3. The simulated s are increasingly wider as a y decreases. If a y is fixed, increasing a x leads to increase sinuosity, see Figure 4. 1 C A (7) 4.3 Full linear transform We can simultaneously rotate and squeeze the patterns present in the training image. The corresponding coordinate transform is written: u new = AR u old u =(x; y; z) (8) where u new and u old are the old and new coordinates of each single datum location of the mp data event. R and A are the rotation and ainity matrices. In 2D, this is written: ( x new y new = a x [x old = a y [x old cos yold sin ] (9) sin + yold cos ] Figure 5a shows the sequence of rotation and ainity. The simulated realization, corresponding to a 30 i rotation angle and ainity factors a x =1:0, a y =3:0, is displayed in Figure 5b. That realization is similar to that of Figure 2b except for the 30 i rotation. With a 45 i rotation and a x = a y =3:0, the realization is displayed in Figure 5c; it is seen to be similar to that of Figure 4b. 4.4 Location-dependent transform Stationarity is a statistics requirement rarely exactly met by natural phenomena. Instead, geological patterns typically evolve gradually from one area to another; this 7
8 gradual deformation can be modeled by a series of location-dependent linear transforms (rotations and ainities). If these location-dependent transforms can be evaluated, e.g. from seismic data, the stationary patterns of the training image can be gradually transformed to deliver the locally variable patterns desired. For example, a fan deposit may be seen as consisting of three parts: upper fan, middle fan and lower fan. A main feeder (or a few feed s) split and produce many smaller braided s as they move forward from the upper fan to the lower fan. Hence, the s tend to deviate from the feeder direction and be increasingly smaller. The location-dependent transform is written similarly to expression (8): u new = A(u 0 )R (u0 )u old u =(x; y; z) (10) where u new and u old are the old and new coordinates of each single datum location of the mp data event; u 0 is the center of that data event; R (u0 ) and A(u 0 ) are the rotation and ainity matrices which are location u 0 dependent. In 2D, this is can be developed as: ( x new y new (u) =a x(u 0 )[x old cos( (u 0)) y old (u) =a y (u 0 )[x old sin( (u 0))] sin( (u 0)) + y old cos( (u 0))] Consider the simulation of a fan deposit in which the feeder s are known to move from northwest to southeast. A large stationary training image, twice larger than the simulated grid size of is used, see Figure 6a. At each simulated grid node, a dierent set of rotation angle and ainity factors is used to impose local characteristics of that fan deposit. The maps of local rotation angles and ainity factors used are shown in Figures 6b and 6c. Figure 6d gives one resulting simulated realization. The format of the data file containing rotation angles and ainity factors as well as the corresponding comments are given in Appendix of this report. 8
9 Assume the following notation: l A is the binary (indicator) random variable associated to the occurrence of state s k at location u. l B is the binary random variable associated to the occurrence of the data event constituted by the n hard conditioning data S(u + h )=s k, =1; 2; :::; n, considered jointly. l C is the binary random variable associated to the occurrence of the data event constituted by the n Y soft conditioning data Y (u+h 0 fi )=y fi, fi =1; 2; :::; n Y, considered jointly. If seismic information can be calibrated to provide the partial conditional probability P (AjC), it can be combined with the other partial probability P (AjB) obtained from hard data to give the full conditional probability P (AjB;C), Journel(2000). Two realizations corresponding to dierent combinations of P (AjB) and P (AjC) are shown in Figure 7; the second (Figure 7b) gives more influence to the seismic component P (AjC). A more complex example is shown in Figure 8. In this case, a composite of fan deposits is considered in which the continuity direction and width of the s are dierent in dierent lobes of the composite fan. The spatial distributions of local rotation angles and ainity factors are given in Figure 8b and 8c. The resulting simulated realization is displayed in Figure 8d. 4.5 Non-stationary examples A non-stationary representation of a delta fan cannot be used as a training image because it is too location-specific and its patterns are not repeated over the training area. Figure 9 shows one such non-stationary training image of a delta fan and one resulting simulated realization. The characteristics of the original training delta fan are not reproduced. Indeed, if a non-stationary training image is scanned, the location-dependent patterns of this training image are averaged out by the multiplepoint statistics and results in the average stationary patterns seen on the simulated realization. This fundamental remark is true for all statistics including the variogram models of 2-point geostatistics. 9
10 5 Preliminary conclusions From the above discussion and preliminary test runs, the following conclusions can be made: l Training patterns should be stationary with rotation and ainity invariance. l Location-specific patterns associated to a specific reservoir can be reproduced through coordinates linear transform (rotation and ainity). The multiple-point statistics can be read from a modular stationary training image. l Very dierent facies distributions, although all related to the same geological sedimentary style, can be generated from one single modular training image. Reference [1] S. Strebelle., Sequential simulation drawing structures from training images. SCRF report, 2000, and PhD thesis. [2] W. Xu., Conditional curvilinear stochastic simulation using pixel-based algorithms. Mathematical Geology, 28(7). [3] C.V.Deutsch and A.G.Journel., GSLIB - Geostatistical Software Library and User s Guide. Oxford University Press,
11 2500 a- Training image East b- Realization with angle = c- Realization with angle = 45 East 1500 East d- Realization with angle = 90 East 1500 Figure 1: Training image (a) and simulated realizations with dierent rotation angles (b), (c), and (d) 11
12 1500 a- Realization with ax=1.0, ay= b- Realization with ax=1.0, ay=3.0 East 1500 East c- Realization with ax=1.0, ay= d- Realization with ax=1.0, ay=5.0 East 1500 East 1500 Figure 2: Realizations with decreasing width (a y > 1) 12
13 1500 a- Realization with ax=1.0, ay= b- Realization with ax=1.0, ay=0.8 East 1500 East c- Realization with ax=1.0, ay= d- Realization with ax=1.0, ay=0.6 East 1500 East 1500 Figure 3: Realizations with increasing width (a y < 1) 13
14 1500 a- Realization with ax=2.0, ay= b- Realization with ax=3.0, ay=3.0 East 1500 East c- Realization with ax=4.0, ay= d- Realization with ax=5.0, ay=3.0 East 1500 East 1500 Figure 4: Realizations with increasing sinuosity (a x > 1) 14
15 a - Sequence of rotation and ainity Original data event Ainity u 4 u 3 w 3 u 1 θ u 2 w 4 w 1 w 2 v 3 v 4 v 2 v 1 Rotation b- Realization with angle=30, ax=1.0, ay= c- Realization with angle=45, ax=ay=3.0 East 1500 East 1500 Figure 5: Realizations with dierent rotation angles and ainity factors 15
16 5000 a- Large training image East b- Rotation angles c- Ainity factors East 2500 East d- One realization East 2500 Figure 6: Training image for a delta fan (a), maps of rotation angles (b) and ainity factors (c), one simulated realization (d). 16
17 5000 a- Large training image East b- P(A C) c- One realization (omega = 0.2) East 2500 East d- One realization (omega = 1.0) East 2500 Figure 7: Training image (a), seismic data given as P(AjC) (b), and two realizations with dierent emphasis given to the C-seismic data (c) and (d). 17
18 5000 a- Large training image East b- Rotation angles c- Ainity factors East 2500 East d- one realization East 2500 Figure 8: Training image (a), rotation angles (b), ainity factors (c), and one realization (d). 18
19 2500 Training image East One realization East 2500 Figure 9: Non-stationary delta fan used as a training image and one resulting simulated realization 19
Multiple 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationGeostatistical 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 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 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 informationSTOCHASTIC OBJECT-BASED SIMULATION OF CHANNELS CONSTRAINED ABSTRACT INTRODUCTION 1 DATA AND GEOLOGICAL FRAMEWORK BY HIGH RESOLUTION SEISMIC DATA
STOCHASTIC OBJECT-BASED SIMULATION OF CHANNELS CONSTRAINED BY HIGH RESOLUTION SEISMIC DATA S. Viseur, A. Shtuka, J-L. Mallet ABSTRACT Recent progresses in exploration techniques (sonar images, high resolution
More informationChapter 18. Geometric Operations
Chapter 18 Geometric Operations To this point, the image processing operations have computed the gray value (digital count) of the output image pixel based on the gray values of one or more input pixels;
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 informationHistory matching under training-image based geological model constraints
History matching under training-image based geological model constraints JEF CAERS Stanford University, Department of Petroleum Engineering Stanford, CA 94305-2220 January 2, 2002 Corresponding author
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 informationSIMPAT: Stochastic Simulation with Patterns
SIMPAT: Stochastic Simulation with Patterns G. Burc Arpat Stanford Center for Reservoir Forecasting Stanford University, Stanford, CA 94305-2220 April 26, 2004 Abstract Flow in a reservoir is mostly controlled
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 informationA009 HISTORY MATCHING WITH THE PROBABILITY PERTURBATION METHOD APPLICATION TO A NORTH SEA RESERVOIR
1 A009 HISTORY MATCHING WITH THE PROBABILITY PERTURBATION METHOD APPLICATION TO A NORTH SEA RESERVOIR B. Todd HOFFMAN and Jef CAERS Stanford University, Petroleum Engineering, Stanford CA 94305-2220 USA
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 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 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 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 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 informationGraphics and Interaction Transformation geometry and homogeneous coordinates
433-324 Graphics and Interaction Transformation geometry and homogeneous coordinates Department of Computer Science and Software Engineering The Lecture outline Introduction Vectors and matrices Translation
More informationCOMP30019 Graphics and Interaction Transformation geometry and homogeneous coordinates
COMP30019 Graphics and Interaction Transformation geometry and homogeneous coordinates Department of Computer Science and Software Engineering The Lecture outline Introduction Vectors and matrices Translation
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 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 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 informationComputer Graphics Hands-on
Computer Graphics Hands-on Two-Dimensional Transformations Objectives Visualize the fundamental 2D geometric operations translation, rotation about the origin, and scale about the origin Learn how to compose
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 informationSpatial Analysis and Modeling (GIST 4302/5302) Guofeng Cao Department of Geosciences Texas Tech University
Spatial Analysis and Modeling (GIST 4302/5302) Guofeng Cao Department of Geosciences Texas Tech University 1 Outline of This Week Last topic, we learned: Spatial autocorrelation of areal data Spatial regression
More informationGrid-less Simulation of a Fluvio-Deltaic Environment
Grid-less Simulation of a Fluvio-Deltaic Environment R. Mohan Srivastava, FSS Canada Consultants, Toronto, Canada MoSrivastava@fssconsultants.ca and Marko Maucec, Halliburton Consulting and Project Management,
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 informationRM03 Integrating Petro-elastic Seismic Inversion and Static Model Building
RM03 Integrating Petro-elastic Seismic Inversion and Static Model Building P. Gelderblom* (Shell Global Solutions International BV) SUMMARY This presentation focuses on various aspects of how the results
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 informationHierarchical modeling of multi-scale flow barriers in channelized reservoirs
Hierarchical modeling of multi-scale flow barriers in channelized reservoirs Hongmei Li and Jef Caers Stanford Center for Reservoir Forecasting Stanford University Abstract Channelized reservoirs often
More informationUsing 3D-DEGA. Omer Inanc Tureyen and Jef Caers Department of Petroleum Engineering Stanford University
Using 3D-DEGA Omer Inanc Tureyen and Jef Caers Department of Petroleum Engineering Stanford University 1 1 Introduction With the advance of CPU power, numerical reservoir models have become an essential
More informationTerrain Modeling with ArcView GIS from ArcUser magazine
Lesson 7: Rotating Point Symbols Lesson Goal: Correctly orient geologic symbols by joining adjusted Azimuth and Bearing data to point data. What You Will Need: A Pentium class PC with 32 MB of RAM (minimum)
More informationDetermination of 6D-workspaces of Gough-type parallel. manipulator and comparison between different geometries. J-P. Merlet
Determination of 6D-workspaces of Gough-type parallel manipulator and comparison between different geometries J-P. Merlet INRIA Sophia-Antipolis, France Abstract: We consider in this paper a Gough-type
More informationMODULE - 7. Subject: Computer Science. Module: Other 2D Transformations. Module No: CS/CGV/7
MODULE - 7 e-pg Pathshala Subject: Computer Science Paper: Computer Graphics and Visualization Module: Other 2D Transformations Module No: CS/CGV/7 Quadrant e-text Objectives: To get introduced to the
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 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 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 informationDigital Image Processing ERRATA. Wilhelm Burger Mark J. Burge. An algorithmic introduction using Java. Second Edition. Springer
Wilhelm Burger Mark J. Burge Digital Image Processing An algorithmic introduction using Java Second Edition ERRATA Springer Berlin Heidelberg NewYork Hong Kong London Milano Paris Tokyo 12.1 RGB Color
More informationDynamic data integration and stochastic inversion of a two-dimensional confined aquifer
Hydrology Days 2013 Dynamic data integration and stochastic inversion of a two-dimensional confined aquifer Dongdong Wang 1 Ye Zhang 1 Juraj Irsa 1 Abstract. Much work has been done in developing and applying
More informationGEOG 4110/5100 Advanced Remote Sensing Lecture 4
GEOG 4110/5100 Advanced Remote Sensing Lecture 4 Geometric Distortion Relevant Reading: Richards, Sections 2.11-2.17 Geometric Distortion Geometric Distortion: Errors in image geometry, (location, dimensions,
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 informationScanner Parameter Estimation Using Bilevel Scans of Star Charts
ICDAR, Seattle WA September Scanner Parameter Estimation Using Bilevel Scans of Star Charts Elisa H. Barney Smith Electrical and Computer Engineering Department Boise State University, Boise, Idaho 8375
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 informationAQA GCSE Maths - Higher Self-Assessment Checklist
AQA GCSE Maths - Higher Self-Assessment Checklist Number 1 Use place value when calculating with decimals. 1 Order positive and negative integers and decimals using the symbols =,, , and. 1 Round to
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 informationHierarchical Trend Models Based on Architectural Elements
Hierarchical Trend Models Based on Architectural Elements Michael J. Pyrcz (mpyrcz@ualberta.ca) and Clayton V. Deutsch (cdeutsch@ualberta.ca) Department of Civil & Environmental Engineering University
More informationChapter - 2: Geometry and Line Generations
Chapter - 2: Geometry and Line Generations In Computer graphics, various application ranges in different areas like entertainment to scientific image processing. In defining this all application mathematics
More informationMaximum A Posteriori Selection with Homotopic Constraint
Maximum A Posteriori Selection with Homotopic Constraint Michael J. Pyrcz and Clayton V. Deutsch Department of Civil & Environmental Engineering University of Alberta Abstract The addition of homotopic
More informationFluid flow modelling with seismic cluster analysis
Fluid flow modelling with seismic cluster analysis Fluid flow modelling with seismic cluster analysis Laurence R. Bentley, Xuri Huang 1 and Claude Laflamme 2 ABSTRACT Cluster analysis is used to construct
More informationCamera Model and Calibration
Camera Model and Calibration Lecture-10 Camera Calibration Determine extrinsic and intrinsic parameters of camera Extrinsic 3D location and orientation of camera Intrinsic Focal length The size of the
More informationEffects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines
Effects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY Velocity models used for wavefield-based seismic
More informationChapter 3 Image Registration. Chapter 3 Image Registration
Chapter 3 Image Registration Distributed Algorithms for Introduction (1) Definition: Image Registration Input: 2 images of the same scene but taken from different perspectives Goal: Identify transformation
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 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 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 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 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 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 informationFast stochastic simulation of patterns using a multi-scale search
1 Fast stochastic simulation of patterns using a multi-scale search Pejman Tahmasebi and Jef Caers Department of Energy Resources Engineering Stanford University Abstract Pattern-based spatial modeling
More informationMotion. 1 Introduction. 2 Optical Flow. Sohaib A Khan. 2.1 Brightness Constancy Equation
Motion Sohaib A Khan 1 Introduction So far, we have dealing with single images of a static scene taken by a fixed camera. Here we will deal with sequence of images taken at different time intervals. Motion
More informationModeling Uncertainty in the Earth Sciences Jef Caers Stanford University
Modeling spatial continuity Modeling Uncertainty in the Earth Sciences Jef Caers Stanford University Motivation uncertain uncertain certain or uncertain uncertain Spatial Input parameters Spatial Stochastic
More informationCurriki Geometry Glossary
Curriki Geometry Glossary The following terms are used throughout the Curriki Geometry projects and represent the core vocabulary and concepts that students should know to meet Common Core State Standards.
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 informationComputer Vision. Coordinates. Prof. Flávio Cardeal DECOM / CEFET- MG.
Computer Vision Coordinates Prof. Flávio Cardeal DECOM / CEFET- MG cardeal@decom.cefetmg.br Abstract This lecture discusses world coordinates and homogeneous coordinates, as well as provides an overview
More informationIntroduction to Homogeneous coordinates
Last class we considered smooth translations and rotations of the camera coordinate system and the resulting motions of points in the image projection plane. These two transformations were expressed mathematically
More informationDoyle Spiral Circle Packings Animated
Doyle Spiral Circle Packings Animated Alan Sutcliffe 4 Binfield Road Wokingham RG40 1SL, UK E-mail: nsutcliffe@ntlworld.com Abstract Doyle spiral circle packings are described. Two such packings illustrate
More informationN-Views (1) Homographies and Projection
CS 4495 Computer Vision N-Views (1) Homographies and Projection Aaron Bobick School of Interactive Computing Administrivia PS 2: Get SDD and Normalized Correlation working for a given windows size say
More informationMAT 343 Laboratory 4 Plotting and computer animation in MATLAB
MAT 4 Laboratory 4 Plotting and computer animation in MATLAB In this laboratory session we will learn how to. Plot in MATLAB. The geometric properties of special types of matrices (rotations, dilations,
More informationDirect Sequential Co-simulation with Joint Probability Distributions
Math Geosci (2010) 42: 269 292 DOI 10.1007/s11004-010-9265-x Direct Sequential Co-simulation with Joint Probability Distributions Ana Horta Amílcar Soares Received: 13 May 2009 / Accepted: 3 January 2010
More informationSCRF the art of building on fundamentals
SCRF 2003-2004 the art of building on fundamentals The SCRF crowd at the Getty Museum, Los Angeles, February 2004 It is a recurrent pleasure to report to you every year about research at SCRF. We have
More informationVector Algebra Transformations. Lecture 4
Vector Algebra Transformations Lecture 4 Cornell CS4620 Fall 2008 Lecture 4 2008 Steve Marschner 1 Geometry A part of mathematics concerned with questions of size, shape, and relative positions of figures
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 informationTIPS4Math Grades 4 to 6 Overview Grade 4 Grade 5 Grade 6 Collect, Organize, and Display Primary Data (4+ days)
Collect, Organize, and Display Primary Data (4+ days) Collect, Organize, Display and Interpret Categorical Data (5+ days) 4m88 Collect data by conducting a survey or an experiment to do with the 4m89 Collect
More information1. Carlos A. Felippa, Introduction to Finite Element Methods,
Chapter Finite Element Methods In this chapter we will consider how one can model the deformation of solid objects under the influence of external (and possibly internal) forces. As we shall see, the coupled
More informationSTOCHASTIC SIMULATION OF PATTERNS USING DISTANCE-BASED PATTERN MODELING
STOCHASTIC SIMULATION OF PATTERNS USING DISTANCE-BASED PATTERN MODELING A DISSERTATION SUBMITTED TO THE DEPARTMENT OF ENERGY RESOURCES ENGINEERING AND THE COMMITTEE ON GRADUATE STUDIES OF STANFORD UNIVERSITY
More informationImplicit Boundary Modeling (BOUNDSIM)
Implicit Boundary Modeling (BOUNDSIM) Jason A. McLennan and Clayton V. Deutsch Centre for Computational Geostatistics Department of Civil & Environmental Engineering University of Alberta Boundary modeling
More informationYEAR 7 MATHS SCHEMES OF WORK
YEAR 7 MATHS SCHEMES OF WORK 2018-19 Autumn Term (September-December) SOW: We follow the STP Mathematics course in year 7, published by Oxford University Press. fanez@al-ashraf.gloucs.sch.uk Unit Topic
More informationHow reliable is statistical wavelet estimation? Jonathan Edgar*, BG Group, and Mirko van der Baan, University of Alberta
Jonathan Edgar*, BG Group, and Mirko van der Baan, University of Alberta Summary Well logs are often used for the estimation of seismic wavelets. The phase is obtained by forcing a well-derived synthetic
More informationIf three points A (h, 0), P (a, b) and B (0, k) lie on a line, show that: a b 1.
ASSIGNMENT ON STRAIGHT LINES LEVEL 1 (CBSE/NCERT/STATE BOARDS) 1 Find the angle between the lines joining the points (0, 0), (2, 3) and the points (2, 2), (3, 5). 2 What is the value of y so that the line
More informationDijkstra's Algorithm
Shortest Path Algorithm Dijkstra's Algorithm To find the shortest path from the origin node to the destination node No matrix calculation Floyd s Algorithm To find all the shortest paths from the nodes
More informationTu N Directional Grey Level Co-occurrence Matrixbased Attributes for Fracture Detection
Tu N105 11 Directional Grey Level Co-occurrence Matrixbased Attributes for Fracture Detection C. Eichkitz* (Joanneum Research Forschungesellschaft mbh), J. Amtmann (Joanneum Research Forschungesellschaft
More information