Introduction to Geostatistics | Course Notes

1 Introduction to Geostatistics Course NotesYe ZhangDept. of Geology & GeophysicsUniversity of WyomingDraft date January 12, Overview..iv1 Why Geostatistics ?.. Geostatistical Prediction.. Geostatistics versus Simple Interpolation.. Limitations.. This Class.. References.. Outline.. Homework, Grades, Tools.. Suggested Reading..122 Probability Theory Nomenclature and Notation.. Univariate Analysis.. Introduction .. Formal Definitions.. Random Variable Arithmetic.. Bivariate Analysis.. Introduction .. Bivariate Random Variables.

2 Bivariate Arithmetics.. Multivariate Analysis.. Linear Combination of Many .. Multivariate Arithmetics.. Gaussian Distribution.. Central Limit Theorem.. Chapter Project.. Suggested Reading..31iiiCONTENTS3 Spatial Conventional Analysis (Nongeostatistical).. Data Posting.. Contour Map.. Symbol Map.. Indicator Map.. Moving Window.. Spatial Continuity Analysis (Geostatistical).. Experimental Variogram.. h-Scatterplot.. Variogram versus Univariate Statistics.. Higher Dimensions & Statistical Anisotropy.

3 Pure Nugget Variogram.. Standard Deviation of (h).. Irregular Data: Variogram Search Envelope.. Exploring Anisotropy.. Spatial Continuity Analysis: Outline.. Chapter Project.. Suggested Reading..744 Variogram Basic Permissible Models.. Positive Definiteness.. Basic Variogram Functions.. 1D Variogram Model Statistical Isotropy.. Higher Dimensional Model Statistical Anisotropy.. Nested Model.. Model-Fitting Rules of Thumb .. Chapter Project.. Suggested Reading..895 The Problem of Estimation.. Nongeostatistical (Deterministic) Estimation.

4 Global Estimation.. Local Estimation.. Estimation Criteria.. Geostatistical (Probabilistic) Estimation.. Random Function Models.. Ordinary Kriging.. Cokriging & Collocated Cokriging.. Kriging with moving neighborhood.. Chapter Project.. Suggested Reading..119 CONTENTSiii6 Advanced Cross Validation.. Other Estimators.. Simple Kriging.. Indicator Kriging.. Block Kriging.. Geostatistical Simulation.. Cholesky Decomposition.. Conditional Simulation.. Sequential Gaussian Simulation (SGS).. Chapter Project.

5 Suggested Reading..1417 Final Overview.. What to hand in.. Additional Notes ..1488 Geostatistics Final Exam149ivCONTENTSThis is the lecture note written & assembled by Ye Zhang for an introductorycourse in 2010 GEOL 54463 CREDITSA-F GRADINGPre-requisite: Calculus I Linear Algebra; Probability & Statistics;Matlab programming languageLocation: ESB1006 Times: TTh (9:35 am 10:50 pm)Office hour: M(4:00 5:30 pm), F(3:00 4:30 pm), GE 220 Email: 307-766-2981 The syllabus: see : The lecture note do not include: (1) solutions to the ex-ercises and projects.

6 (2) proofs to theories and equation will be presented only during lectures. So, please do not relyon the Notes for everything class attendance and participation arekey to doing OverviewGeoscientists often face interpolation and estimation problems when analyzingsparse data from field observations. Geostatistics is an invaluable tool that canbe used to characterize spatial or temporal phenomena1. Geostatistics orig-inated from the mining and petroleum industries, starting with the work byDanie Krige in the 1950 s and was further developed by Georges Matheron inthe 1960 s.

7 In both industries, Geostatistics is successfully applied to solve caseswhere decisions concerning expensive operations are based on interpretationsfrom sparse data located in space. Geostatistics has since been extended tomany other fields in or related to the earth sciences, , hydrogeology, hydrol-ogy, meteorology, oceanography, geochemistry, geography, soil sciences, forestry,landscape ecology. In this class, both fundamental development of geostatisticsand simple, practical applications in the earth sciences will be presented. Ex-ercises and projects are designed to help elucidate the fundamental assignments will be given illustrating the applications of geostatisticsin the particular field of reservoir characterization and this class, we re concerned only with spatial analysis; temporal phenomena might bebetter understood in a separate class on time series 1 OverviewWhat is Geostatistics ?

8 Data analysis and spatial continuity modeling (Journel,1989). Establish quantitative measure of spatial correlation to be used for sub-sequent estimation and simulation (Deutsch, 2002). The following introductionand overview materials are based on compilation of several source materials (seefull references in Sec. ). Why Geostatistics ?Classic statistics is generally devoted to the analysis and interpretation of un-certainties caused by limited sampling of a property under study. Geostatisticshowever deviates from classic statistics in that Geostatistics is not tied to apopulation distribution model that assumes, for example, all samples of a pop-ulation are normally distributed and independent from one another.

9 Most ofthe earth science data ( , rock properties, contaminant concentrations) oftendo not satisfy these assumptions as they can be highly skewed and/or possessspatial correlation ( , data values from locations that are closer together tendto be more similar than data values from locations that are further apart).To most geologists, the fact that closely spaced samples tend to be similar isnot surprising since such samples have been influenced by similar physical andchemical depositional/transport to the classic statistics which examine the statistical distributionof a set of sampled data, Geostatistics incorporates both the statistical distri-bution of the sample dataandthe spatial correlation among the sample of this difference, many earth science problems are more effectively ad-dressed using geostatistical methods.

10 As stated by Marc Cromer (inGeostatis-tics for environmental and geotechnical applications, 1996, ASTM International,edited by Rouhani et al.):Geostatistical methods provide the tools to capture, through rigorous exami-nation, the descriptive information on a phenomenon from sparse, often biased,and often expensive sample data. The continued examination and quantitative12 CHAPTER 1. OVERVIEW rigor of the procedure provide a vehicle for integrating qualitative and quantita-tive understanding by allowing the data to speak for themselves.