Heterogeneity: Why important? What to do?
|
|
- Mabel Wilcox
- 6 years ago
- Views:
Transcription
1 Heterogeneity: Why important? What to do? Jesus Carrera Institut Jaume Almera (IJA) for Earth Sciences Higher Council for Scientific Research (CSIC) Barcelona, Spain 1
2 Transport Effect of space and time variability How to characterize. 2
3 Solute transport Advection: v = q/φ (q proportional to K) Dispersion: Proportional to: α q Reactions Mass Conservation c c c f r φ = ( D ) q + t Porosity Dispersion coeff. Water flux Reactions 3
4 Solution of ADE c ( x,t) c 1/ 2 ( 2π ) = + 0 ( ) 18 u e 1 x Vt y + z exp + V d 2 2tDL 2tDT Initial Pulse At early times, little displacement, significant dilution and spreading Later on, dilution and spreading continue, but displacement becomes apparent 4
5 Dispersion and advection scales Longitud (-) Longitud LD LA = = vt 2αvt Tiempo (-) Tiempo c max M M = = φ 2 πdt / φ φ LD Dispersion dominates small distances and times, but later advection becomes dominant. Everything is controlled by D, K and 5 φ
6 Field data Does all the above work? Is transport defined by K, porosity and dispersivity? 6
7 Scale dependence of Hyd. Conductivity (FEBEX) log K (m s -1 ) PULSE HORNER CROSS-HOLE GLOBAL MOD. Fractures Rock Matrix One by one test Inflow SCALE OF INFLUENCE (m) 7
8 Scale dependence of dispersivity Data from tracer tests and pollution plumes worldwide Dispersivity grows with test scale (Lallemand-Barres y Peandecerf, 1978 ) 8
9 Kinematic porosity function of residence time 0 log b (m) Dolomita 5 6 log time (min) Hard rocks Creta pendiente 0.5 Efective porosity in fractured rocks appears to increase with residence time 9
10 Transport predictions are awful!!! Good calibration Good calibration El Cabril (UPC, 1999) 10
11 Reaction rates Observed reaction rates are 2-3 orders of magnitude slower than expected from measurements Physical surface area (log m 2 /kg water) Bruton (1986) White (1986) Liu (1987) Hurd (1973) Herman & Lorah (1987) Delany (1985) Velbel (1985) Claassen & White (1973) Paces (1983) Gislason & Eugster (1987) (White & Peterson,1990) Reactive surface area (log m 2 /kg water) 11
12 The problem MVT deposit outflow Expected Diss/Precip Results Conservative species Resident water Dissolution (Corbera et al, 2003 ) Precip. Brine inflow BUT, Do we know how to simulate mixing? 12
13 1st conclusion: ADE full of problems Scale dependence of K Scale dependence of dispersivity Time (and space, and direction) dependence of porosity Poor prediction of space distribution and time evolution of c and reaction THE CULPRIT: Variability Heterogeneity of K (hydraulic conductivity), and minerals, and... (BUT ESPECIALLY K) Temporal fluctuations of water flux 13
14 Random functions Random variable: concept used for describing variables with an uncertain value e.g.: result from throwing a dice z(ξ) z(ξ) z(ξ 1 )=3, z(ξ 2 )=5, z(ξ 3 )=2, z(ξ 4 )=4, ξ 1, ξ 2, ξ 3, are realizations Random functions: Functions whose values are random variables; or Random entities whose realizations are functions: z(x,ξ) z(x,ξ 1 ) z(x,ξ2 ) z(x,ξ 5 ) 14
15 Random functions need spatial correlation Pdf is not sufficient The uncorrelated field will behave as a homogeneous field (can be treated as constant) Stochastics needed when spatial (and/or temporal) structure (often characterized by correlation function) In passing, this is why recharge is usually viewed as constant (a view that may need to be revised after Howard s presentation) 15
16 Dealing with heterogeneity Deterministic (but uncertain) methods Try to describe reality from geology Stochastic methods Describe reality as one realization of a random function. In practice, work with an ensemble of realizations in the hope that some of them are similar to reality or, rather, that jointly describe the actual uncertainty. 16
17 Stochastic methods Motivation: Deal with variability Procedure: 1) Characterize (stochastically) medium variability of K, minerals, surfaces,... 2) Generate n spatially variable fields 3) Solve them 4) Compute output statistics Solve analytically or numerically 17
18 Methods to describe spatial variability Traditional (stationary) geostatistics Heterogeneity defined by pdf and variogram (or correlation function) Tools: kriging, conditional simulation Sophisticated simulation methods (multipoint geostatistics, Markov chain, etc) Genetic models Very powerful free codes (GSLIB, GEMS) available atstanford web page 18
19 Kriging vs conditional simulation Estimated Y field 1.0 Y Reference X Y field Kriging is an interpolation method. Does not reproduce small scale heterogeneity. 6 1 Y X
20 Traditional (stationary) geostatistics Problem: requires many data or highly uncertain and does not look real 20
21 PluriGaussian Model Problem: requires many trials to be realistic 21
22 Genetic models Try to simulate geological processes. Quite powerful for sedimentary formations 22
23 Genetic models Problem: lots of data, assumptions and CPU, hard to condition to real data 23
24 Conditioning on hydraulic data Problem statement: Find K(x) that Honors K(x i ) point measurements of K Reproduces h m (x i,t), head observations Inverse modelling 24
25 Parametrization Express model properties, which may vary in space and time, as a function of a few unknown scalars (MODEL PARAMETERS) Y Zonation p 1 p 2 p 3 N q 0 i i i= 1 q( x,t) = q ( x,t) + p α ( x,t) x Y Interpolation Y Small scale variability p 1 p 2 p 3 p 1 p 2 p 3 x x 25
26 Optimum parameters defined as minimum of objective function h(p) p 2 h(p 1 ) h(p 0 ) p 0 h(p) t Contour lines of objective function: p 2 p 1 h(p 2 ) t Sum (h-h m ) 2 p 1 26
27 Generic optimization method 1) Initialization: Read input data, i=0, p 0 h(p) h(p 2 ) 2) Solve the simulation problem, h(p i ) 3) Compute an updating vector, d. 4) Update parameters, p i+1 = p i + d. p 2 h(p 1 ) h(p 0 ) t 5) If convergence, stop. Otherwise, i=i+1 and go to 2 p 0 Contour lines of objective function p 2 p 1 p 27 1
28 Step 1: Random Field Example: Stochastic inversion Step 2: Simulate drawdowns from three pump tests to be used as data 28
29 A synthetic example of hydraulic tomography B B Data: Drawdowns at 13 obs. Points, caused by three pumping tests Point measurements of T at 13 points O-1 B-1 B-2 O-2 Estimate assuming: B O-9 B-3 O-4 O-7 O-6 O-3 O-10 O-8 O-5 Known variogram Varying number of pilot points, whose prior estimates come from kriging 29
30 Method Given: - Local T data (e.g. Pulses, sp. capacity..) - Cross hole drawdowns - Conceptual model (i.e., B.C. s, geost,...) Find T field such that: - reproduces cross hole drawdowns - reproduces local T data Approach: - maximize likelihood - adjoint state method 30
31 10 Conditioning on hydraulic data O-1 O-9 B-1 B-3 O-4 O-7 O-3 O-6 B-2 O-10 O-2 O-8 O-5 Joint interpretation 10 of 3 pumping tests Drawdown (m) Observation point B Time (s) Observation point O-5 N p = 41; RMSE d = 0.8 N p =241; RMSE d = 0.6 Measured drawdown O-1 O-9 B-1 B-3 O-4 O-7 O-3 Drawdown (m) O B-2 O Time(s) O-2 O-8 O-5 Observation point O-10 Observation point O-8 N p = 41; RMSE d = 1.29 N p =241; RMSE d = 1.19 Measured drawdown Perfect fit Drawdown (m) 1 Drawdown (m) N p = 41; RMSE d = 2.1 N p =241; RMSE d = 1.1 Measured drawdown N p = 41; RMSE d = 1.9 N p =241; RMSE d = 0.9 Measured drawdown Time(s)
32 Hydraulic data insufficient for unequivocal identification λ=.001 Opt. λ λ = 100 N p =41 N p =65 μ= 10 2 μ= 10 1 Few pilot points do not allow identifying narrow features Small weight to prior information leads to unstable results, but provides food for thought N p =97 N p =241 μ= 10 1 Too much weight to prior information leads to stable but boring results Many parameters (p.p.) OK if properly regularized μ=
33 Blind prediction of tracer tests Each simulation reproduces: Spatial variability Tailing Local concentrations 33
34 Conditioning on geophysical data - Problem statement: Find K(x) that Honors K(x i ) point measurements of K Correlates with geophysical image Methods: External drift Cokriging (cosimulation) Qualitative conditioning Any hope to include basic geological (genetic) info? 34
35 In practice, qualitative use 35
36 Fractured rock site (hydro calibration) Descenso (m) Descenso (m) Descenso (m) SR4-1 (punto de bombeo) Falla Sur - S10 Fractura SR4-2 (Unidad alterada) Tiempo (s) S10 (Fractura S10) Tiempo (s) North Fault 474 Fault South Fault S ` Dike Sr4 F.S10 M I N E Sr1 S5 F.SR1-3 PM 27 Dike 285 Fault Sr3 S5 (Frac.285, 27' Dique) Tiempo (s) S14 Sr3 Dike Sr2 Y Z X Descenso (m) Descenso (m) Descenso (m) PM (Mina) Tiempo (s) SR1-1 (Bandas fracturadas) SR1-2 (Dique 27) SR1-3 (SR1-3 Fractura) SR1-4 (Unidad alterada) Tiempo (s) 36
37 Qualitative integration of geophysics facilitates reasonable long term prediction Descenso (m) Descenso (m) 10 S14: Matriz Externa S10: S10 Fractura SR4-3: Unidad alterada 5 10 SR4-2: Matriz Mina 5 10 SR4-1:Falla Sud - S10 Frac Time (s) PM: Mina (Punto Bombeo) Tiempo (s) Descenso (m) Descenso (m) S14 Filon Sr3 F. Norte Filon 27 Filon 27- Sr2 PM F. 474 M I N A Sr3 Sr1 S5 S10 Sr4 F. Sur 15 SR1-1: 27 Dique SR1-2: Matriz Mina SR1-3: Unidad alterada 5 F. 285 Y X Z Descenso (m) 10 SR2-1: Falla Norte 5 10 SR2-2: Unidad alterada SR3-1: Falla Norte SR3-2: Matriz Mina SR3-3: Unidad alterada S5: 285 Falla - 27' Dique Tiempo (s) Tiempo (s) 37
38 Conditioning on geological data - Geostatistics: Facies simulations Honors point textural information Correlates with known patterns Methods: Traditional geostatistics Indicator simulation Multiple Point Geostatistics (training imag) Transition Probabilities Any hope to include basic geological (genetic) info? (i.e., nonstationarity) 38
39 SUMMARY Spatial variability needed to understand observations (at least, transport observations) Very active research field Work still needed on: Refinement in geophysics (shall we just keep working with images?) Methodological developments in geostatistics (process- Vs. Structure- imitating) Joint inversion Use history of the Earth (i.e., geology) info. 39
Geostatistics for Seismic Data Integration in Earth Models
2003 Distinguished Instructor Short Course Distinguished Instructor Series, No. 6 sponsored by the Society of Exploration Geophysicists European Association of Geoscientists & Engineers SUB Gottingen 7
More informationInverting hydraulic heads in an alluvial aquifer constrained with ERT data through MPS and PPM: a case study
Inverting hydraulic heads in an alluvial aquifer constrained with ERT data through MPS and PPM: a case study Hermans T. 1, Scheidt C. 2, Caers J. 2, Nguyen F. 1 1 University of Liege, Applied Geophysics
More informationCO 2 storage capacity and injectivity analysis through the integrated reservoir modelling
CO 2 storage capacity and injectivity analysis through the integrated reservoir modelling Dr. Liuqi Wang Geoscience Australia CO 2 Geological Storage and Technology Training School of CAGS Beijing, P.
More informationGeostatistical History Matching coupled with Adaptive Stochastic Sampling: A zonation-based approach using Direct Sequential Simulation
Geostatistical History Matching coupled with Adaptive Stochastic Sampling: A zonation-based approach using Direct Sequential Simulation Eduardo Barrela* Instituto Superior Técnico, Av. Rovisco Pais 1,
More informationTraining Venue and Dates Ref # Reservoir Geophysics October, 2019 $ 6,500 London
Training Title RESERVOIR GEOPHYSICS Training Duration 5 days Training Venue and Dates Ref # Reservoir Geophysics DE035 5 07 11 October, 2019 $ 6,500 London In any of the 5 star hotels. The exact venue
More informationContents 1 Introduction 2 Statistical Tools and Concepts
1 Introduction... 1 1.1 Objectives and Approach... 1 1.2 Scope of Resource Modeling... 2 1.3 Critical Aspects... 2 1.3.1 Data Assembly and Data Quality... 2 1.3.2 Geologic Model and Definition of Estimation
More informationSimultaneous use of hydrogeological and geophysical data for groundwater protection zone delineation by co-conditional stochastic simulations
Simultaneous use of hydrogeological and geophysical data for groundwater protection zone delineation by co-conditional stochastic simulations C. Rentier,, A. Dassargues Hydrogeology Group, Departement
More informationCalibration of hydraulic and tracer tests in fractured media represented by a DFN model
Calibration and Reliability in Groundwater Modelling: From Uncertainty to Decision Making (Proceedings of ModelCARE 2005, The Hague, The Netherlands, June 2005). IAHS Publ. 304, 2006. 87 Calibration of
More informationStatistical Rock Physics
Statistical - Introduction Book review 3.1-3.3 Min Sun March. 13, 2009 Outline. What is Statistical. Why we need Statistical. How Statistical works Statistical Rock physics Information theory Statistics
More informationLithium Brines: What can we learn from established brine production regions?
Lithium Brines: What can we learn from established brine production regions? An Overview Of Lithium Brine Exploration for Resource Estimation Camilo de los Hoyos Ph.D, Senior Consultant (Geochemistry &
More informationEnsemble Kalman filter assimilation of transient groundwater flow data via stochastic moment equations
Ensemble Kalman filter assimilation of transient groundwater flow data via stochastic moment equations Alberto Guadagnini (1,), Marco Panzeri (1), Monica Riva (1,), Shlomo P. Neuman () (1) Department of
More informationReservoir characterization
1/15 Reservoir characterization This paper gives an overview of the activities in geostatistics for the Petroleum industry in the domain of reservoir characterization. This description has been simplified
More informationA hybrid Marquardt-Simulated Annealing method for solving the groundwater inverse problem
Calibration and Reliability in Groundwater Modelling (Proceedings of the ModelCARE 99 Conference held at Zurich, Switzerland, September 1999). IAHS Publ. no. 265, 2000. 157 A hybrid Marquardt-Simulated
More informationUnderground nuclear waste storage
Underground nuclear waste storage Groundwater flow and radionuclide transport Jan-Olof Selroos Cargese Summer School, July 5, 2018 Contents: Concept for geological disposal of nuclear waste A few words
More informationIntegration of seismic and fluid-flow data: a two-way road linked by rock physics
Integration of seismic and fluid-flow data: a two-way road linked by rock physics Abstract Yunyue (Elita) Li, Yi Shen, and Peter K. Kang Geologic model building of the subsurface is a complicated and lengthy
More informationGeostatistical History Matching coupled with Adaptive Stochastic Sampling
Geostatistical History Matching coupled with Adaptive Stochastic Sampling A Geologically consistent approach using Stochastic Sequential Simulation Eduardo Barrela Nº 79909 Project Thesis Presentation
More information11280 Electrical Resistivity Tomography Time-lapse Monitoring of Three-dimensional Synthetic Tracer Test Experiments
11280 Electrical Resistivity Tomography Time-lapse Monitoring of Three-dimensional Synthetic Tracer Test Experiments M. Camporese (University of Padova), G. Cassiani* (University of Padova), R. Deiana
More informationReservoir Uncertainty Calculation by Large Scale Modeling
Reservoir Uncertainty Calculation by Large Scale Modeling Naeem Alshehri and Clayton V. Deutsch It is important to have a good estimate of the amount of oil or gas in a reservoir. The uncertainty in reserve
More informationAdvances in Locally Varying Anisotropy With MDS
Paper 102, CCG Annual Report 11, 2009 ( 2009) Advances in Locally Varying Anisotropy With MDS J.B. Boisvert and C. V. Deutsch Often, geology displays non-linear features such as veins, channels or folds/faults
More informationMultiple-Point Geostatistics: from Theory to Practice Sebastien Strebelle 1
Multiple-Point Geostatistics: from Theory to Practice Sebastien Strebelle 1 Abstract The limitations of variogram-based simulation programs to model complex, yet fairly common, geological elements, e.g.
More informationOptimizing Thresholds in Truncated Pluri-Gaussian Simulation
Optimizing Thresholds in Truncated Pluri-Gaussian Simulation Samaneh Sadeghi and Jeff B. Boisvert Truncated pluri-gaussian simulation (TPGS) is an extension of truncated Gaussian simulation. This method
More informationRecent developments in object modelling opens new era for characterization of fluvial reservoirs
Recent developments in object modelling opens new era for characterization of fluvial reservoirs Markus L. Vevle 1*, Arne Skorstad 1 and Julie Vonnet 1 present and discuss different techniques applied
More informationMultiple realizations: Model variance and data uncertainty
Stanford Exploration Project, Report 108, April 29, 2001, pages 1?? Multiple realizations: Model variance and data uncertainty Robert G. Clapp 1 ABSTRACT Geophysicists typically produce a single model,
More informationDevelopment and Application of Groundwater Flow and Solute Transport Models. Randolf Rausch
Development and Application of Groundwater Flow and Solute Transport Models Randolf Rausch Overview Groundwater Flow Modeling Solute Transport Modeling Inverse Problem in Groundwater Modeling Groundwater
More informationAcceptable Ergodic Fluctuations and Simulation of Skewed Distributions
Acceptable Ergodic Fluctuations and Simulation of Skewed Distributions Oy Leuangthong, Jason McLennan and Clayton V. Deutsch Centre for Computational Geostatistics Department of Civil & Environmental Engineering
More information7 Geostatistics. Figure 7.1 Focus of geostatistics
7 Geostatistics 7.1 Introduction Geostatistics is the part of statistics that is concerned with geo-referenced data, i.e. data that are linked to spatial coordinates. To describe the spatial variation
More informationComputational Challenges in Reservoir Modeling. Sanjay Srinivasan The Pennsylvania State University
Computational Challenges in Reservoir Modeling Sanjay Srinivasan The Pennsylvania State University Well Data 3D view of well paths Inspired by an offshore development 4 platforms 2 vertical wells 2 deviated
More informationHydrogeology in clay till. Timo KESSLER, Knud Erik KLINT, Bertel NILSSON and Poul L. BJERG
Hydrogeology in clay till Timo KESSLER, Knud Erik KLINT, Bertel NILSSON and Poul L. BJERG motivation Risk assessment of contaminated sites How do contaminants spread in the subsurface? Is geological heterogeneity
More information23855 Rock Physics Constraints on Seismic Inversion
23855 Rock Physics Constraints on Seismic Inversion M. Sams* (Ikon Science Ltd) & D. Saussus (Ikon Science) SUMMARY Seismic data are bandlimited, offset limited and noisy. Consequently interpretation of
More informationCOLLOCATED CO-SIMULATION USING PROBABILITY AGGREGATION
COLLOCATED CO-SIMULATION USING PROBABILITY AGGREGATION G. MARIETHOZ, PH. RENARD, R. FROIDEVAUX 2. CHYN, University of Neuchâtel, rue Emile Argand, CH - 2009 Neuchâtel, Switzerland 2 FSS Consultants, 9,
More informationProspects for river discharge and depth estimation through assimilation of swath altimetry into a raster-based hydraulics model
Prospects for river discharge and depth estimation through assimilation of swath altimetry into a raster-based hydraulics model Kostas Andreadis 1, Elizabeth Clark 2, Dennis Lettenmaier 1, and Doug Alsdorf
More informationOak Ridge IFRC. Quantification of Plume-Scale Flow Architecture and Recharge Processes
Oak Ridge IFRC Quantification of Plume-Scale Flow Architecture and Recharge Processes S. Hubbard *1, G.S. Baker *2, D. Watson *3, D. Gaines *3, J. Chen *1, M. Kowalsky *1, E. Gasperikova *1, B. Spalding
More informationDownloaded 10/25/16 to Redistribution subject to SEG license or copyright; see Terms of Use at
Facies modeling in unconventional reservoirs using seismic derived facies probabilities Reinaldo J. Michelena*, Omar G. Angola, and Kevin S. Godbey, ireservoir.com, Inc. Summary We present in this paper
More informationStructural and Petrophysical Characterization of Mixed Drain/Barrier Fault Zones in Carbonates: Example from the Castellas Fault (SE France)*
Structural and Petrophysical Characterization of Mixed Drain/Barrier Fault Zones in Carbonates: Example from the Castellas Fault (SE France)* Christophe Matonti 1, Juliette Lamarche 1, Yves Guglielmi 1,
More informationSub-kilometer-scale space-time stochastic rainfall simulation
Picture: Huw Alexander Ogilvie Sub-kilometer-scale space-time stochastic rainfall simulation Lionel Benoit (University of Lausanne) Gregoire Mariethoz (University of Lausanne) Denis Allard (INRA Avignon)
More informationTransiogram: A spatial relationship measure for categorical data
International Journal of Geographical Information Science Vol. 20, No. 6, July 2006, 693 699 Technical Note Transiogram: A spatial relationship measure for categorical data WEIDONG LI* Department of Geography,
More informationConditional Distribution Fitting of High Dimensional Stationary Data
Conditional Distribution Fitting of High Dimensional Stationary Data Miguel Cuba and Oy Leuangthong The second order stationary assumption implies the spatial variability defined by the variogram is constant
More informationCapturing aquifer heterogeneity: Comparison of approaches through controlled sandbox experiments
WATER RESOURCES RESEARCH, VOL. 47, W09514, doi:10.1029/2011wr010429, 2011 Capturing aquifer heterogeneity: Comparison of approaches through controlled sandbox experiments Steven J. Berg 1 and Walter A.
More informationCombining geological surface data and geostatistical model for Enhanced Subsurface geological model
Combining geological surface data and geostatistical model for Enhanced Subsurface geological model M. Kurniawan Alfadli, Nanda Natasia, Iyan Haryanto Faculty of Geological Engineering Jalan Raya Bandung
More informationPetrophysical seismic inversion conditioned to well-log data: Methods and application to a gas reservoir
GEOPHYSICS, VOL. 4, NO. MARCH-APRIL 009 ; P. O1 O1, 11 FIGS. 10.1190/1.30439 Petrophysical seismic inversion conditioned to well-log data: Methods and application to a gas reservoir Miguel Bosch 1, Carla
More informationTraps for the Unwary Subsurface Geoscientist
Traps for the Unwary Subsurface Geoscientist ashley.francis@sorviodvnvm.co.uk http://www.sorviodvnvm.co.uk Presented at SEG Development & Production Forum, 24-29 th June 2001, Taos, New Mexico, USA 24-29
More informationHydraulic tomography: Development of a new aquifer test method
WATER RESOURCES RESEARCH, VOL. 36, NO. 8, PAGES 2095 2105, AUGUST 2000 Hydraulic tomography: Development of a new aquifer test method T.-C. Jim Yeh and Shuyun Liu Department of Hydrology and Water Resources,
More informationTime-lapse filtering and improved repeatability with automatic factorial co-kriging. Thierry Coléou CGG Reservoir Services Massy
Time-lapse filtering and improved repeatability with automatic factorial co-kriging. Thierry Coléou CGG Reservoir Services Massy 1 Outline Introduction Variogram and Autocorrelation Factorial Kriging Factorial
More informationBME STUDIES OF STOCHASTIC DIFFERENTIAL EQUATIONS REPRESENTING PHYSICAL LAW
7 VIII. BME STUDIES OF STOCHASTIC DIFFERENTIAL EQUATIONS REPRESENTING PHYSICAL LAW A wide variety of natural processes are described using physical laws. A physical law may be expressed by means of an
More informationCOLLOCATED CO-SIMULATION USING PROBABILITY AGGREGATION
COLLOCATED CO-SIMULATION USING PROBABILITY AGGREGATION G. MARIETHOZ, PH. RENARD, R. FROIDEVAUX 2. CHYN, University of Neuchâtel, rue Emile Argand, CH - 2009 Neuchâtel, Switzerland 2 FSS Consultants, 9,
More informationTRUNCATED GAUSSIAN AND PLURIGAUSSIAN SIMULATIONS OF LITHOLOGICAL UNITS IN MANSA MINA DEPOSIT
TRUNCATED GAUSSIAN AND PLURIGAUSSIAN SIMULATIONS OF LITHOLOGICAL UNITS IN MANSA MINA DEPOSIT RODRIGO RIQUELME T, GAËLLE LE LOC H 2 and PEDRO CARRASCO C. CODELCO, Santiago, Chile 2 Geostatistics team, Geosciences
More informationIntegration of Rock Physics Models in a Geostatistical Seismic Inversion for Reservoir Rock Properties
Integration of Rock Physics Models in a Geostatistical Seismic Inversion for Reservoir Rock Properties Amaro C. 1 Abstract: The main goal of reservoir modeling and characterization is the inference of
More informationSeismic inversion for reservoir properties combining statistical rock physics and geostatistics: A review
GEOPHYSICS, VOL. 75, NO. 5 SEPTEMBER-OCTOBER 2010 ; P. 75A165 75A176, 8 FIGS. 10.1190/1.3478209 Seismic inversion for reservoir properties combining statistical rock physics and geostatistics: A review
More informationIdentifying sources of a conservative groundwater contaminant using backward probabilities conditioned on measured concentrations
WATER RESOURCES RESEARCH, VOL. 42,, doi:10.1029/2005wr004115, 2006 Identifying sources of a conservative groundwater contaminant using backward probabilities conditioned on measured concentrations Roseanna
More informationReliability of Seismic Data for Hydrocarbon Reservoir Characterization
Reliability of Seismic Data for Hydrocarbon Reservoir Characterization Geetartha Dutta (gdutta@stanford.edu) December 10, 2015 Abstract Seismic data helps in better characterization of hydrocarbon reservoirs.
More informationThe importance of understanding coupled processes in geothermal reservoirs. Thomas Driesner October 19, 2016
The importance of understanding coupled processes in geothermal reservoirs Thomas Driesner October 19, 2016 Findings from natural hydrothermal systems Interaction of permeability and fluid properties The
More informationBest Unbiased Ensemble Linearization and the Quasi-Linear Kalman Ensemble Generator
W. Nowak a,b Best Unbiased Ensemble Linearization and the Quasi-Linear Kalman Ensemble Generator Stuttgart, February 2009 a Institute of Hydraulic Engineering (LH 2 ), University of Stuttgart, Pfaffenwaldring
More informationFred Mayer 1; Graham Cain 1; Carmen Dumitrescu 2; (1) Devon Canada; (2) Terra-IQ Ltd. Summary
2401377 Statistically Improved Resistivity and Density Estimation From Multicomponent Seismic Data: Case Study from the Lower Cretaceous McMurray Formation, Athabasca Oil Sands Fred Mayer 1; Graham Cain
More informationExtreme Value Analysis and Spatial Extremes
Extreme Value Analysis and Department of Statistics Purdue University 11/07/2013 Outline Motivation 1 Motivation 2 Extreme Value Theorem and 3 Bayesian Hierarchical Models Copula Models Max-stable Models
More informationJAPEX s 60 Years Experience Exploring Volcanic Reservoirs in Japan*
Click to view oral presentation (3.00 MB) JAPEX s 60 Years Experience Exploring Volcanic Reservoirs in Japan* Kentaro Takeda 1 and Yasuo Yamada 2 Search and Discovery Article #70297 (2017)** Posted October
More information4th HR-HU and 15th HU geomathematical congress Geomathematics as Geoscience Reliability enhancement of groundwater estimations
Reliability enhancement of groundwater estimations Zoltán Zsolt Fehér 1,2, János Rakonczai 1, 1 Institute of Geoscience, University of Szeged, H-6722 Szeged, Hungary, 2 e-mail: zzfeher@geo.u-szeged.hu
More informationFinding Large Capacity Groundwater Supplies for Irrigation
Finding Large Capacity Groundwater Supplies for Irrigation December 14, 2012 Presented by: Michael L. Chapman, Jr., PG Irrigation Well Site Evaluation Background Investigation Identify Hydrogeologic Conditions
More informationInterpretation and application of hydrogeological concepts to commercial scale brine extraction projects
Interpretation and application of hydrogeological concepts to commercial scale brine extraction projects ABSTRACT Terry Braun, Pablo Cortegoso*, Cristian Pereira* and Vladimir Ugorets* 1. SRK Consulting
More informationWATER RESOURCES RESEARCH, VOL. 40, W01506, doi: /2003wr002253, 2004
WATER RESOURCES RESEARCH, VOL. 40, W01506, doi:10.1029/2003wr002253, 2004 Stochastic inverse mapping of hydraulic conductivity and sorption partitioning coefficient fields conditioning on nonreactive and
More informationThe Generalized Likelihood Uncertainty Estimation methodology
CHAPTER 4 The Generalized Likelihood Uncertainty Estimation methodology Calibration and uncertainty estimation based upon a statistical framework is aimed at finding an optimal set of models, parameters
More informationAnalysis of hydraulic and tracer response tests within moderately fractured rock based on a transition probability geostatistical approach
WATER RESOURCES RESEARCH, VOL. 40, W12404, doi:10.1029/2004wr003188, 2004 Analysis of hydraulic and tracer response tests within moderately fractured rock based on a transition probability geostatistical
More informationNEW GEOLOGIC GRIDS FOR ROBUST GEOSTATISTICAL MODELING OF HYDROCARBON RESERVOIRS
FOR ROBUST GEOSTATISTICAL MODELING OF HYDROCARBON RESERVOIRS EMMANUEL GRINGARTEN, BURC ARPAT, STANISLAS JAYR and JEAN- LAURENT MALLET Paradigm Houston, USA. ABSTRACT Geostatistical modeling of reservoir
More informationOpportunities in Oil and Gas Fields Questions TABLE OF CONTENTS
TABLE OF CONTENTS A. Asset... 3 1. What is the size of the opportunity (size the prize)?... 3 2. Volumetric Evaluation... 3 3. Probabilistic Volume Estimates... 3 4. Material Balance Application... 3 5.
More information11/22/2010. Groundwater in Unconsolidated Deposits. Alluvial (fluvial) deposits. - consist of gravel, sand, silt and clay
Groundwater in Unconsolidated Deposits Alluvial (fluvial) deposits - consist of gravel, sand, silt and clay - laid down by physical processes in rivers and flood plains - major sources for water supplies
More informationSTA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Statistics! rsalakhu@utstat.toronto.edu! http://www.utstat.utoronto.ca/~rsalakhu/ Sidney Smith Hall, Room 6002 Lecture 11 Project
More informationReservoir connectivity uncertainty from stochastic seismic inversion Rémi Moyen* and Philippe M. Doyen (CGGVeritas)
Rémi Moyen* and Philippe M. Doyen (CGGVeritas) Summary Static reservoir connectivity analysis is sometimes based on 3D facies or geobody models defined by combining well data and inverted seismic impedances.
More informationReservoir Characterisation and Modelling for CO 2 Storage
Reservoir Characterisation and Modelling for CO 2 Storage Tess Dance IEA CCS Summer School Perth, December 2015 ENERGY Why build subsurface models? To simulate fluid flow To estimate capacity Predict reservoir
More informationStochastic vs Deterministic Pre-stack Inversion Methods. Brian Russell
Stochastic vs Deterministic Pre-stack Inversion Methods Brian Russell Introduction Seismic reservoir analysis techniques utilize the fact that seismic amplitudes contain information about the geological
More informationMETHODOLOGY WHICH APPLIES GEOSTATISTICS TECHNIQUES TO THE TOPOGRAPHICAL SURVEY
International Journal of Computer Science and Applications, 2008, Vol. 5, No. 3a, pp 67-79 Technomathematics Research Foundation METHODOLOGY WHICH APPLIES GEOSTATISTICS TECHNIQUES TO THE TOPOGRAPHICAL
More information3D geologic modelling of channellized reservoirs: applications in seismic attribute facies classification
first break volume 23, December 2005 technology feature 3D geologic modelling of channellized reservoirs: applications in seismic attribute facies classification Renjun Wen, * president and CEO, Geomodeling
More informationRock Mechanics for Tunneling
Rock Mechanics for Tunneling Discrete Fracture Approach Dr. William Dershowitz Dept of Civil Engineering University of Washington FracMan Technology Group Golder Associates Inc Rock Quality -RQD Rock Quality
More informationChapter 12 - Incorporation of Geostatistics into Flow and Transport Simulations
Chapter 12 - Incorporation of Geostatistics into Flow and Transport Simulations 12.1 Abstract The zeolitic threshold method used in previous chapters is replaced with more rigorous geostatistical modeling.
More informationmonitoring data for the CO2CRC Otway
Calibration of simulation models to monitoring data for the CO2CRC Otway project. Jonathan Ennis-King 1, T. Dance 1, J. Xu 2, C. Boreham 3, B. Freifeld 6, J. Gunning 1, B. Gurevich 4, C.Jenkins 1, L. Paterson
More informationZONAL KRIGING. 5.1: Introduction and Previous Work CHAPTER 5
CHAPTER 5 ZONAL KRIGING Kriging and conditional indicator simulation are valuable tools for evaluating uncertainty in the subsurface, but are limited by the assumption of stationarity (i.e. the mean and
More informationSTA 414/2104: Machine Learning
STA 414/2104: Machine Learning Russ Salakhutdinov Department of Computer Science! Department of Statistics! rsalakhu@cs.toronto.edu! http://www.cs.toronto.edu/~rsalakhu/ Lecture 9 Sequential Data So far
More information2. Governing Equations. 1. Introduction
Multiphysics Between Deep Geothermal Water Cycle, Surface Heat Exchanger Cycle and Geothermal Power Plant Cycle Li Wah Wong *,1, Guido Blöcher 1, Oliver Kastner 1, Günter Zimmermann 1 1 International Centre
More informationBuilding an Integrated Static Reservoir Model 5-day Course
Building an Integrated Static Reservoir Model 5-day Course Prepared by International Reservoir Technologies Lakewood, Colorado http://www.irt-inc.com/ 1 Agenda Day 1 Day 2 Day 3 Day 4 Day 5 Morning Introduction
More informationStochastic Modeling & Petrophysical Analysis of Unconventional Shales: Spraberry-Wolfcamp Example
Stochastic Modeling & Petrophysical Analysis of Unconventional Shales: Spraberry-Wolfcamp Example Fred Jenson and Howard Rael, Fugro-Jason Introduction Recent advances in fracture stimulation techniques
More informationGeologic and Reservoir Characterization and Modeling
Geologic and Reservoir Characterization and Modeling Scott M. Frailey and James Damico Illinois State Geological Survey Midwest Geologic Sequestration Science Conference November 8 th, 2011 Acknowledgements
More informationIJMGE Int. J. Min. & Geo-Eng. Vol.49, No.1, June 2015, pp
IJMGE Int. J. Min. & Geo-Eng. Vol.49, No.1, June 2015, pp.131-142 Joint Bayesian Stochastic Inversion of Well Logs and Seismic Data for Volumetric Uncertainty Analysis Moslem Moradi 1, Omid Asghari 1,
More informationFRACTURED ROCK Characterization and Remediation. Allan Horneman September 30, 2016
FRACTURED ROCK Characterization and Remediation Allan Horneman September 30, 2016 Disclaimers and Notices The materials herein are intended to furnish viewers with a summary and overview of general information
More informationIntroduction to Ensemble Kalman Filters and the Data Assimilation Research Testbed
Introduction to Ensemble Kalman Filters and the Data Assimilation Research Testbed Jeffrey Anderson, Tim Hoar, Nancy Collins NCAR Institute for Math Applied to Geophysics pg 1 What is Data Assimilation?
More informationPETROLEUM GEOSCIENCES GEOLOGY OR GEOPHYSICS MAJOR
PETROLEUM GEOSCIENCES GEOLOGY OR GEOPHYSICS MAJOR APPLIED GRADUATE STUDIES Geology Geophysics GEO1 Introduction to the petroleum geosciences GEO2 Seismic methods GEO3 Multi-scale geological analysis GEO4
More information"Modelling air quality in the city"
"Modelling air quality in the city" Diede Nijmeijer Master thesis University of Twente Applied Mathematics Specialization: Mathematical Physics and Computational Mechanics Chair: Multiscale Modelling and
More informationHyGEM - a new Strategic Research Council project for linking geophysical models to hydrological models
HyGEM - a new Strategic Research Council project for linking geophysical models to hydrological models Anders Vest Christiansen, GEUS Esben Auken, Dept. of Geoscience, Aarhus University HyGEM Integrating
More informationDavid de Courcy-Bower and Samuel Mohr
Applicability and Limitations of LNAPL Transmissivity as a Metric within Bedrock Formations Insert then choose Picture select your picture. Right click your picture and Send to back. David de Courcy-Bower
More informationInverse Modelling for Flow and Transport in Porous Media
Inverse Modelling for Flow and Transport in Porous Media Mauro Giudici 1 Dipartimento di Scienze della Terra, Sezione di Geofisica, Università degli Studi di Milano, Milano, Italy Lecture given at the
More information4.11 Groundwater model
4.11 Groundwater model 4.11 Groundwater model 4.11.1 Introduction and objectives Groundwater models have the potential to make important contributions in the mapping and characterisation of buried valleys.
More informationCHARACTERIZATION OF HETEROGENEITIES AT THE CORE-SCALE USING THE EQUIVALENT STRATIFIED POROUS MEDIUM APPROACH
SCA006-49 /6 CHARACTERIZATION OF HETEROGENEITIES AT THE CORE-SCALE USING THE EQUIVALENT STRATIFIED POROUS MEDIUM APPROACH Mostafa FOURAR LEMTA Ecole des Mines de Nancy, Parc de Saurupt, 54 04 Nancy, France
More informationA comparison of seven geostatistically based inverse approaches to estimate transmissivities for modeling advective transport by groundwater flow
WATER RESOURCES RESEARCH, VOL. 34, NO. 6, PAGES 1373 1413, JUNE 1998 A comparison of seven geostatistically based inverse approaches to estimate transmissivities for modeling advective transport by groundwater
More informationLarge Scale Modeling by Bayesian Updating Techniques
Large Scale Modeling by Bayesian Updating Techniques Weishan Ren Centre for Computational Geostatistics Department of Civil and Environmental Engineering University of Alberta Large scale models are useful
More informationSTATISTICAL MODELS FOR QUANTIFYING THE SPATIAL DISTRIBUTION OF SEASONALLY DERIVED OZONE STANDARDS
STATISTICAL MODELS FOR QUANTIFYING THE SPATIAL DISTRIBUTION OF SEASONALLY DERIVED OZONE STANDARDS Eric Gilleland Douglas Nychka Geophysical Statistics Project National Center for Atmospheric Research Supported
More informationIntroduction to Spatial Data and Models
Introduction to Spatial Data and Models Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. 2 Department of Forestry
More informationKriging in the Presence of LVA Using Dijkstra's Algorithm
Kriging in the Presence of LVA Using Dijkstra's Algorithm Jeff Boisvert and Clayton V. Deutsch One concern with current geostatistical practice is the use of a stationary variogram to describe spatial
More informationA new Hierarchical Bayes approach to ensemble-variational data assimilation
A new Hierarchical Bayes approach to ensemble-variational data assimilation Michael Tsyrulnikov and Alexander Rakitko HydroMetCenter of Russia College Park, 20 Oct 2014 Michael Tsyrulnikov and Alexander
More informationBest Practice Reservoir Characterization for the Alberta Oil Sands
Best Practice Reservoir Characterization for the Alberta Oil Sands Jason A. McLennan and Clayton V. Deutsch Centre for Computational Geostatistics (CCG) Department of Civil and Environmental Engineering
More information3.4 Fuzzy Logic Fuzzy Set Theory Approximate Reasoning Fuzzy Inference Evolutionary Optimization...
Contents 1 Introduction... 1 1.1 The Shale Revolution... 2 1.2 Traditional Modeling... 4 1.3 A Paradigm Shift... 4 2 Modeling Production from Shale... 7 2.1 Reservoir Modeling of Shale... 9 2.2 System
More informationThe Snap lake diamond deposit - mineable resource.
Title Page: Author: Designation: Affiliation: Address: Tel: Email: The Snap lake diamond deposit - mineable resource. Fanie Nel Senior Mineral Resource Analyst De Beers Consolidated Mines Mineral Resource
More informationarxiv: v1 [physics.flu-dyn] 18 Mar 2018
APS Persistent incomplete mixing in reactive flows Alexandre M. Tartakovsky 1 and David Barajas-Solano 1 1 Pacific Northwest National Laboratory, Richland, WA 99352, USA arxiv:1803.06693v1 [physics.flu-dyn]
More informationQUANTITATIVE INTERPRETATION
QUANTITATIVE INTERPRETATION THE AIM OF QUANTITATIVE INTERPRETATION (QI) IS, THROUGH THE USE OF AMPLITUDE ANALYSIS, TO PREDICT LITHOLOGY AND FLUID CONTENT AWAY FROM THE WELL BORE This process should make
More information