Inverting hydraulic heads in an alluvial aquifer constrained with ERT data through MPS and PPM: a case study
|
|
- Basil Daniels
- 6 years ago
- Views:
Transcription
1 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 2 Stanford University, Energy Resources Department GEOENV 2014
2 The construction of hydrogeological models is generally done with sparse data There are generally two types of data : borehole log description (lithology) and hydraulic head levels 2
3 What happens between boreholes???? Geophysical methods provide spatially distributed information Clay Sand Bed-rock Gravel 3
4 The integration of hydrogeological data in geostatistics is not straightforward Geostatistical models are useful tools to integrate various types of data and to impose a priori spatial constraints However, integrating hydrogeological data remains challenging Hydraulic heads and tracer tests data result from a dynamical process which requires to solve flow (and transport) equations Possible solutions : the probability perturbation method, pilot-point method, gradual deformation, MCMC, etc. 4
5 Context and choice of the training image Probability Perturbation Method Synthetic Benchmark Field case : Hermalle-sous-Argenteau (Belgium) 5
6 Context and choice of the training image Probability Perturbation Method Synthetic Benchmark Field case : Hermalle-sous-Argenteau (Belgium) 6
7 We investigate the site of Hermalle-sous-Argenteau (Belgium) in the alluvial plain of the Meuse River 7
8 From our knowledge of the deposits, we propose four different scenarios for training images 3 facies based on high, intermediate and low hydraulic conductivity 2 scenarios for the gravel facies : channels or bars with two sizes proportions : gravel = 42%, sand = 40%, clay = 18% (from borehole) 8
9 We would like to ensure that training images and geophysical models are somehow compatible This would ensure that conditioning with soft data will not create conflict with features of the TI. 9
10 For each training images, 2D sections are randomly extracted to be compared with field models For each model, facies are assigned a value of electrical resistivity ERT data are simulated and inverted to obtain geophysical inverted models that can be compared to true field data 10
11 All the models (field and synthetic) are mapped using multidimensional scaling Hermans et al., IAMG 2013 Park et al., 2013 If the TI-based scenarios are consistent with ERT, field models should fall in the distribution of synthetic models 11
12 The probability of each scenario according to field data is computed using krenel density estimation P(TI ti D d ) k All scenarios have the same prior probability P = 0.25 obs f (d obs ti k )P(TI ti k) f (d ti )P(TI ti ) k obs k k (Park et al., 2013) Training Images conditional probabilities SC MC SB BB Mean SC = small channels SB = small bars MC = medium channels BB = big bars Scenarios with channels have lower probabilities of occurrence than scenarios with bars 12
13 The electrical resistivity distribution can be transformed into probability maps to constrain MPS simulations f ( A)P(A) P(A ) f ( ) = Distribution of resistivity ρ given the facies = Probability to observe a facies given the resistivity ρ 13
14 Context and choice of the training image Probability Perturbation Method Synthetic Benchmark Field case : Hermalle-sous-Argenteau (Belgium) 14
15 The probability perturbation method (PPM) is a stochastic, iterative inversion method optimizing the geometry of the deposits Generate initial model with MPS Run flow simulation with HGS 1. Initialization with prior model constraints 1 calibrated realization Yes N i 1 (h h ) obs calc 2 i i N 3. Inner loop to optimize the perturbation No Change the random seed to generate new random numbers Choose a value of r D the perturbation parameter Generate a model with MPS Run flow simulation with HGS Converge to best r D? 2. Outer loop to ensure the convergence The structure of the algorithm is divided in three steps 15
16 The perturbation is computed through a single parameter in the inner loop We introduce : D = the observed hydraulic heads P(A D) is the probability to have a facies distribution according to these data. It is estimated by P(A D) = (1-r D ) I 0 (u) + r D P(A) I 0 (u) = best matching model from MPS (constrained by TI, hard data and geophysical data) = realization of the joint probability P(A B,C) r D =0 : no perturbation r D =1 : new independent realization P(A B,C,D) is calculated through the tau model 16
17 The process converges towards the target level of the objective function and thus to a model explaining all available data At each iteration r D is optimized The model is coherent with direct geological observations (imposed), near-surface geophysical data (soft data) and hydrogeological data (dynamic data) The process is repeated to obtain several equiprobable models 17
18 Context and choice of the training image Probability Perturbation Method Synthetic Benchmark Field case : Hermalle-sous-Argenteau (Belgium) 18
19 The synthetic benchmark was designed to simulate an alluvial aquifer TI and prior probability Hard Data Data to reproduce Soft Data = ERT data Clay lenses 22% : K = 10-6 m/s Gravel channels 20% : K = 10-2 m/s Sand matrix 58% : K = 10-4 m/s The aim is to study the role of geophysical data and the choice of the training image on final results 19
20 We calculate the posterior distributions based on the mean of 50 realizations With geophysical soft data Mean OF after 10 iterations = m Without geophysical soft data Mean OF after 10 iterations = 0.06 m Geophysical data (ERT) help to build plausible models by giving information on facies between boreholes Through this improvement of models, geophysical data increase the convergence of the method 20
21 The training image influences the geometry of the deposits and consequently modifies groundwater flows and posterior distributions Gravel posterior probability Small channels Gravel posterior probability Medium channels Bigger channels yield more continuous high probability gravel zones in the posterior distribution 21
22 Geophysical data reduce the «error» introduced by the wrong training image (here bigger channels) Gravel posterior probability Medium channels Without geophysics Gravel posterior probability Medium channels With geophysics Geophysical data correct the posterior distribution where the uncertainty related to the presence of gravel should be higher 22
23 Context and choice of the training image Probability Perturbation Method Synthetic Benchmark Field case : Hermalle-sous-Argenteau (Belgium) 23
24 The methodology is applied to the site of Hermalle-sous- Argenteau with a pumping followed in 8 piezometers and the pumping well 24
25 Electrical resistivity models are used to constrain MPS simulations 12 parallel ERT profiles Probability map for each facies (here sand) 25
26 PPM is applied 50 times with a different initial model, posterior distributions are used to classify each cell of the model Posterior probability of gravel Classification 26
27 PPM is applied 50 times with a different initial model, posterior distributions are used to classify each cell of the model Posterior probability of gravel Classification 27
28 The results for the 4 TI are linearly combined using their respective probability to produce final maps Posterior probability of gravel 28
29 The results for the 4 TI are linearly combined using their respective probability to produce final maps Classification 29
30 Conclusion and perspectives 30
31 Development of a global methodology to integrate nearsurface geophysical, geological and hydrogeological data with multiple-point geostatistics TI Uncertainty Soft data Models Hard data + = Hydrogeological data Posterior distributions Classification Application to 3 facies models with a potential geophysical method (ERT) 31
32 The method can now be extended to other geological contexts, other scales and other geophysical methods Fractured limestones Water catchment scale SIP (Robert, 2012) 32
33 THANK YOU VERY MUCH FOR YOUR ATTENTION 33
34 APPENDICES 34
35 Traditional geostatistics (variogram-based, kriging) fail to reproduce multimodal distributions and complex structures 3 different kinds of heterogeneity having similar variograms (Caers and Zhang, 2002) Multiple-point geostatistics (MPS) were developed to overcome this problem 35
36 MPS uses a training image (TI) to describe the expected geological heterogeneity, generally with discrete variables TI representing sinuous fluvial channels partially interconnected 500 m The training image is scanned with a given template of neighbors The reproductions of each possible combination are counted and stored in a search tree 36
37 Sequential simulations are used to produce equiprobable realizations (SNESIM algorithm) The neighborhood of the simulated point is extracted, and compared to the search tree to derive the conditional probability 2 examples of equiprobable realizations Training image 100 m 100 m 500 m 37
38 There is a lack of sedimentological data concerning the alluvial aquifer of the Meuse, to build accurate training image A lot of boreholes are available but poorly described (only basic description of the types of sediments) and far away from each other More than 30 different descriptions of facies e.g.: gravelly loam, clayey gravel, loamy sands, sandy gravel, gravelly sand, etc. Proportions vary from site to site From the sedimentological point of view, we can expect elements from meandering and braided systems 38
39 We will find sand/gravel bars or channels and point bars, loam/clay crevasse splays and floodplain deposits Meandering systems 1 channel Braided systems 1 main channel + secondary channels Channel lag deposits Point bar deposits Levees Crevasse splays Floodplain deposits Channel bars Longitudinal bars Cut and fill structures Floodplain deposits Old channel fill 39
40 The final aim is to provide hydrogeological models of the deposits The facies must be related to the hydrogeological properties 3 facies description 1) Gravel = high hydraulic conductivity = channel or bars 2) Clay and loam = low hydraulic conductivity = floodplain, filling and crevasse splays 3) Sand and sandy gravel = intermediate hydraulic conductivity = bars, levees, low energy channels 40
41 Consistency 41
42 MPS simulations use a TI to depict the expected geological heterogeneity Size of channels? Orientation? Size of lenses? Considerable uncertainty often remains concerning the geometry and architecture of facies elements. 42
43 Geophysical data may provide spatial information on subsurface properties clay gravel sand Classification based on electrical resistivity Electrical resistivity tomography (ERT) yields a model of subsurface resistivity, high values being linked with sand and low values with clay. 43
44 The use of geophysical image to extract TI is limited due to the resolution of geophysics and the philosophy of MPS 1. Geophysical methods are not able to differentiate unequivocally different facies: small scale features may be hidden by larger elements 2. Geophysical methods may be used to constraint spatially MPS simulations, TI and soft data should be independent Training image 44
45 They are compared to inverted field data Example of an ERT profile acquired on the site of Hermalle-sous-Argenteau 45
46 Synthetic inverted sections and true field sections are compared using multi-dimensional scaling (MDS) 1. The Euclidean distance between any two inverted models (electrical resistivity of each cell) is calculated d ( ) ( ) T ij i j i j 2. This defines a metric space in which the L models are defined in terms of distance to each other 3. A change of variable is done such that T B HAH Size L x L with aij dij and H I 11 2 L 4. The eigenvalue decomposition of B offers a way to project the models into a d- dimension map, using the d first eigenvalues and eigenvectors. B V V T B B B MDS X V 1/2 d B,d B,d 5. X d are the coordinates of the model in the MDS map of d-dimensions 46
47 The probability analysis of different scenarios can be seen as a conditional probability problem P(TI ti D d ) k obs f (d obs ti k )P(TI ti k) f (d ti )P(TI ti ) k obs k k d obs is the projection of the response from each synthetic model indicating a location on the MDS map We can estimate f (d obs ti k ) using the adaptive kernel density estimation technique with a bivariate normal function applied on the MDS map (x ) (x ) f (x, x ) exp( [ ]) 2 X x 2 y x y 2 x y (Park et al., 2013) x 1 and x 2 are the coordinates of the models, μ x and μ y correspond to the coordinates of the reference value (the field data for example), σ x and σ y are bandwith parameters depending on the cluster 47
48 The proposed methodology enables to verify the «consistency» of trainingimage based scenarios with geophysical data using synthetic benchmarks This method is a way to justify and validate the choice of a training image and to propose alternatives for uncertainty analysis and stochastic simulations 48
49 Important remarks and limitations 1. The choice of resistivity distribution plays an important role, the geophysicist must ensure that reliable values are chosen 2. The choice of the Euclidean distance is locally dependent which may be a drawback in MPS since training image are not locally dependent 3. It would be possible to use other types of distance measurements or other attributes than resistivity 4. The choice of the dimensions of the map is often 2D or 3D for visual inspection, here it represents about 60% of the total variance, higher dimensions are possible 5. 3D geophysical surveys are more costly and time consuming, but the methodology can be extended to 3D models as well, it would be useful to distinguish between channels and elongated bars 6. The methodology shows that there is some consistency between TI and geophysical data, but not that the TI is perfect or that other scenarios are impossible 49
50 The initial and the history matched models both respect the same prior constraints, i.e. the same geostatistical parameters Some features are quite similar (e.g. red zone), other features are completely different (e.g. blue zone) 50
51 ERT data are simulated and a 3D model of resistivity is created and transformed in probability maps Gravel probability map Clay probability map High probability of gravel are related to the presence of gravel channels as imaged by ERT 51
Multiple-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 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 informationApplication of multiple-point geostatistics on modelling groundwater flow and transport in the Brussels Sands
Application of multiple-point geostatistics on modelling groundwater flow and transport in the Brussels Sands Marijke Huysmans (1) and Alain Dassargues (1,2) (1) Earth and Environmental Sciences, K.U.Leuven,
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 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 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 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 informationUSING GEOSTATISTICS TO DESCRIBE COMPLEX A PRIORI INFORMATION FOR INVERSE PROBLEMS THOMAS M. HANSEN 1,2, KLAUS MOSEGAARD 2 and KNUD S.
USING GEOSTATISTICS TO DESCRIBE COMPLEX A PRIORI INFORMATION FOR INVERSE PROBLEMS THOMAS M. HANSEN 1,2, KLAUS MOSEGAARD 2 and KNUD S. CORDUA 1 1 Institute of Geography & Geology, University of Copenhagen,
More informationReducing Uncertainty in Modelling Fluvial Reservoirs by using Intelligent Geological Priors
Reducing Uncertainty in Modelling Fluvial Reservoirs by using Intelligent Geological Priors Temístocles Rojas 1, Vasily Demyanov 2, Mike Christie 3 & Dan Arnold 4 Abstract Automatic history matching reservoir
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 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 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 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 information2. REGIS II: PARAMETERIZATION OF A LAYER-BASED HYDROGEOLOGICAL MODEL
HYDRAULIC PARAMETERIZATION OF 3D SUBSURFACE MODELS: FROM MEASUREMENT-SCALE TO MODEL-SCALE Jan L. Gunnink, Jan Stafleu, Denise Maljers and Jan Hummelman TNO Geological Survey of the Netherlands, P.O. Box
More informationGY 111: Physical Geology
UNIVERSITY OF SOUTH ALABAMA GY 111: Physical Geology Lecture 14: Sedimentary Rocks Part 2: Alluvial Fans and Rivers Instructor: Dr. Douglas W. Haywick Last Time 1) Types of Sediment 2) Sedimentary Rock
More informationGY 402: Sedimentary Petrology
UNIVERSITY OF SOUTH ALABAMA GY 402: Sedimentary Petrology Lecture 17: Sandy Fluvial Depositional Environments Instructor: Dr. Douglas W. Haywick Last Time Volcaniclastic Sedimentary Rocks 1. Origin of
More informationModeling uncertainty in metric space. Jef Caers Stanford University Stanford, California, USA
Modeling uncertainty in metric space Jef Caers Stanford University Stanford, California, USA Contributors Celine Scheidt (Research Associate) Kwangwon Park (PhD student) Motivation Modeling uncertainty
More informationMulti-point Modeling of Clay Lenses and its Impact on Aquifer Vulnerability
Multi-point Modeling of Clay Lenses and its Impact on Aquifer Vulnerability Whitney Trainor and Jef Caers Earth, Energy, and Environmental Sciences Stanford Center for Reservoir Forecasting April 30 &
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 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 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 informationA008 THE PROBABILITY PERTURBATION METHOD AN ALTERNATIVE TO A TRADITIONAL BAYESIAN APPROACH FOR SOLVING INVERSE PROBLEMS
A008 THE PROAILITY PERTURATION METHOD AN ALTERNATIVE TO A TRADITIONAL AYESIAN APPROAH FOR SOLVING INVERSE PROLEMS Jef AERS Stanford University, Petroleum Engineering, Stanford A 94305-2220 USA Abstract
More informationStreams. Water. Hydrologic Cycle. Geol 104: Streams
Streams Why study streams? Running water is the most important geologic agent in erosion, transportation and deposition of sediments. Water The unique physical and chemical properties of water make it
More informationNATURAL RIVER. Karima Attia Nile Research Institute
NATURAL RIVER CHARACTERISTICS Karima Attia Nile Research Institute NATURAL RIVER DEFINITION NATURAL RIVER DEFINITION Is natural stream of water that flows in channels with ih more or less defined banks.
More informationComparing the gradual deformation with the probability perturbation method
Comparing the gradual deformation with the probability perturbation method Jef Caers Petroleum Engineering Department Stanford University, California USA 1. Introduction This paper is intended as a comparison
More informationWe LHR3 04 Realistic Uncertainty Quantification in Geostatistical Seismic Reservoir Characterization
We LHR3 04 Realistic Uncertainty Quantification in Geostatistical Seismic Reservoir Characterization A. Moradi Tehrani* (CGG), A. Stallone (Roma Tre University), R. Bornard (CGG) & S. Boudon (CGG) SUMMARY
More informationAppendix D. Sediment Texture and Other Soil Data
5 6 7 8 Appendix D. Sediment Texture and Other Soil Data This appendix describes the sediment texture of the aquifer system in the Restoration Area. The contents of this appendix describe the: Importance
More informationGeostatistics 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 informationAPPLICATION OF ELECTRICAL RESISTIVITY TOMOGRAPHY FOR SAND UNDERWATER EXTRACTION
International Scientific Conference GEOBALCANICA 2018 APPLICATION OF ELECTRICAL RESISTIVITY TOMOGRAPHY FOR SAND UNDERWATER EXTRACTION Maya Grigorova Ivaylo Koprev University of Mining and Geology St. Ivan
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 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 informationApplication of Transition Probability Geostatistics in a Detailed Stratigraphic Framework. Gary Weissmann University of New Mexico
Application of Transition Probability Geostatistics in a Detailed Stratigraphic Framework Gary Weissmann University of New Mexico ACKNOWLEDGEMENTS National Science Foundation American Chemical Society
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 informationINTEGRATION OF BOREHOLE INFORMATION AND RESISTIVITY DATA FOR AQUIFER VULNERABILITY
INTEGRATION OF BOREHOLE INFORMATION AND RESISTIVITY DATA FOR AQUIFER VULNERABILITY ANDERS V. CHRISTIANSEN 1, ESBEN AUKEN 1 AND KURT SØRENSEN 1 1 HydroGeophysics Group, Aarhus University, Finlandsgade 8,
More informationBasin Analysis Applied to Modelling Buried Valleys in the Great Lakes Basin
EARTH SCIENCES SECTOR GENERAL INFORMATION PRODUCT 35 Basin Analysis Applied to Modelling Buried Valleys in the Great Lakes Basin Sharpe, D R; Russell, H A J 2004 Originally released as: Basin Analysis
More informationAssessing uncertainty on Net-to-gross at the Appraisal Stage: Application to a West Africa Deep-Water Reservoir
Assessing uncertainty on Net-to-gross at the Appraisal Stage: Application to a West Africa Deep-Water Reservoir Amisha Maharaja April 25, 2006 Abstract A large data set is available from a deep-water reservoir
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 informationStochastic Simulation of Inclined Heterolithic Stratification with Streamline-based Stochastic Models
Stochastic Simulation of Inclined Heterolithic Stratification with Streamline-based Stochastic Models Michael J. Pyrcz and Clayton V. Deutsch Department of Civil and Environmental Engineering, University
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 information27. Running Water I (p ; )
27. Running Water I (p. 424-436; 440-444) Hydrosphere How much of the Earth s surface is covered by water? Earth's water is collectively called the and is stored in a number of so-called as follows: 1.
More informationBootstrap confidence intervals for reservoir model selection techniques
Bootstrap confidence intervals for reservoir model selection techniques Céline Scheidt and Jef Caers Department of Energy Resources Engineering Stanford University Abstract Stochastic spatial simulation
More informationDirect forecasting without full model inversion Jef Caers
Direct forecasting without full model inversion Jef Caers Professor of Geological Sciences Stanford University, USA Why moving to Geological Sciences? To attract a greater diversity of students To be able
More informationLecture 19: Fluvial Facies
GEOL 440 Sedimentology and stratigraphy: processes, environments and deposits Lecture 19: Fluvial Facies Aims Examine controls on alluvial channel pattern List the types of channel pattern Examine meandering
More informationChallenges in conditioning a stochastic geological model of a heterogeneous glacial aquifer to a comprehensive soft data set
doi:10.5194/hess-18-2907-2014 Author(s) 2014. CC Attribution 3.0 License. Challenges in conditioning a stochastic geological model of a heterogeneous glacial aquifer to a comprehensive soft data set J.
More informationHierarchical Geostatistical Analysis of an Experimental Stratigraphy
DOI 10.1007/s11004-008-9180-6 Hierarchical Geostatistical Analysis of an Experimental Stratigraphy Ye Zhang Received: 31 August 2006 / Accepted: 28 January 2008 International Association for Mathematical
More informationRADIONUCLIDE DIFFUSION IN GEOLOGICAL MEDIA
GEOPHYSICS RADIONUCLIDE DIFFUSION IN GEOLOGICAL MEDIA C. BUCUR 1, M. OLTEANU 1, M. PAVELESCU 2 1 Institute for Nuclear Research, Pitesti, Romania, crina.bucur@scn.ro 2 Academy of Scientists Bucharest,
More informationGY 402: Sedimentary Petrology
UNIVERSITY OF SOUTH ALABAMA GY 402: Sedimentary Petrology Lecture 13: Immature Siliciclastic Sedimentary Environments Alluvial Fans, Braided Streams Instructor: Dr. Douglas W. Haywick Last Time Immature
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 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 informationHeterogeneity: Why important? What to do?
Heterogeneity: Why important? What to do? Jesus Carrera Institut Jaume Almera (IJA) for Earth Sciences Higher Council for Scientific Research (CSIC) Barcelona, Spain 1 Transport Effect of space and time
More informationAnalyzing the effects of geological and parameter uncertainty on prediction of groundwater head and travel time
Model Development ess doi:10.5194/hess-17-3245-2013 Author(s) 2013. CC Attribution 3.0 License. Hydrology and Earth System Sciences Open Access Ocean Science Analyzing the effects of geological and parameter
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 informationSASKATCHEWAN STRATIGRAPHY GLACIAL EXAMPLE BOULDERS IN GLACIAL DEPOSITS
SASKATCHEWAN STRATIGRAPHY GLACIAL EXAMPLE BOULDERS IN GLACIAL DEPOSITS 51 SASKATCHEWAN STRATIGRAPHY GLACIAL SURFICIAL STRATIFIED DEPOSITS 52 SASKATCHEWAN STRATIGRAPHY GLACIAL EXAMPLE OF SEDIMENT DEPOSITION
More informationA scale/resolution problem
Near Surface 2005 September 4-7, 2005, Palermo Italy WORKSHOP on Hydrogeophysics a tool for sustainable use of groundwater resources Integrated Seismic and GPR characterization of fractured rocks Pipan,
More informationRESERVOIR SEISMIC CHARACTERISATION OF THIN SANDS IN WEST SYBERIA
www.senergyltd.com RESERVOIR SEISMIC CHARACTERISATION OF THIN SANDS IN WEST SYBERIA Erick Alvarez, Jaume Hernandez, Bolkhotivin E.A., Belov A.V., Hakima Ben Meradi,Jonathan Hall, Olivier Siccardi, Phil
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 informationImpact of the Danube River on the groundwater dynamics in the Kozloduy Lowland
GEOLOGICA BALCANICA, 46 (2), Sofia, Nov. 2017, pp. 33 39. Impact of the Danube River on the groundwater dynamics in the Kozloduy Lowland Peter Gerginov Geological Institute, Bulgarian Academy of Sciences,
More informationSeismic tomography with co-located soft data
Seismic tomography with co-located soft data Mohammad Maysami and Robert G. Clapp ABSTRACT There is a wide range of uncertainties present in seismic data. Limited subsurface illumination is also common,
More informationEntropy of Gaussian Random Functions and Consequences in Geostatistics
Entropy of Gaussian Random Functions and Consequences in Geostatistics Paula Larrondo (larrondo@ualberta.ca) Department of Civil & Environmental Engineering University of Alberta Abstract Sequential Gaussian
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 information' International Institute for Land Reclamation and Improvement. 2 Groundwater Investigations. N.A. de Ridder'? 2.1 Introduction. 2.
2 Groundwater Investigations N.A. de Ridder'? 2.1 Introduction Successful drainage depends largely on a proper diagnosis of the causes of the excess water. For this diagnosis, one must consider: climate,
More informationQuantifying uncertainty of geological 3D layer models, constructed with a-priori
Quantifying uncertainty of geological 3D layer models, constructed with a-priori geological expertise Jan Gunnink, Denise Maljers 2 and Jan Hummelman 2, TNO Built Environment and Geosciences Geological
More informationGroundwater Sustainability at Wadi Al Bih Dam, Ras El Khaimah, United Arab Emirates (UAE) using Geophysical methods
Groundwater Sustainability at Wadi Al Bih Dam, Ras El Khaimah, United Arab Emirates (UAE) using Geophysical methods Ahmed Murad, Amir Gabr, Saber Mahmoud, Hasan Arman & Abdulla Al Dhuhoori Geology Department
More informationReview of permeability in buried-valley aquifers: Centimeter to kilometer scales
Review of permeability in buried-valley aquifers: Centimeter to kilometer scales R. W. RITZI, Z. DAI & D. F. DOMINIC Department of Geological Sciences, Wright State University, Dayton, Ohio 45435, USA
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 informationQuantitative Seismic Interpretation An Earth Modeling Perspective
Quantitative Seismic Interpretation An Earth Modeling Perspective Damien Thenin*, RPS, Calgary, AB, Canada TheninD@rpsgroup.com Ron Larson, RPS, Calgary, AB, Canada LarsonR@rpsgroup.com Summary Earth models
More information10. Facies Modeling Sequential Indicator Simulation (SIS)
10. Facies Modeling 10.1 Sequential Indicator Simulation (SIS) In this section the process of building a basic facies model conditioned to well observations using SIS. The variogram type, ranges, and azimuth
More informationINTRODUCTION Water resources are vitally important for the future of humankind. Groundwater from karst aquifers is among the most important drinking w
Vulnerability Mapping for the Protection of Carbonate (Karst) Aquifers (Ramallah-Al Bireh District) Eng. Amjad da Assi House of fw Water and de Environment May, 2007 INTRODUCTION Water resources are vitally
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 informationResistivity & IP methods
International PhD Course in HYDROGEOPHYSICS Resistivity & IP methods Andrew Binley Lancaster University Overview We have demonstrated links between hydrological and geophysical properties and show the
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 informationQuantitative Relation of the Point BarWidth and Meander Belt Width of Subsurface Reservoir
Advances in Petroleum Exploration and Development Vol. 13, No. 1, 2017, pp. 17-24 DOI:10.3968/9465 ISSN 1925-542X [Print] ISSN 1925-5438 [Online] www.cscanada.net www.cscanada.org Quantitative Relation
More informationReservoir Modeling with GSLIB. Overview
Reservoir Modeling with GSLIB Overview Objectives of the Course What is Geostatistics? Why Geostatistics / 3-D Modeling? Uncertainty Quantification and Decision Making Heterogeneous Reservoir Modeling
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 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 informationAnalysis of the temporal and spatial variations of the chloride concentration in the Río Mojotoro basin
Analysis of the temporal and spatial variations of the chloride concentration in the Río Mojotoro basin Jacqueline Köhn, Eduardo Kruse & Jerónimo Ainchil Facultad de Ciencias Astronómicas y Geofísicas
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 informationFUTURE MEANDER BEND MIGRATION AND FLOODPLAIN DEVELOPMENT PATTERNS NEAR RIVER MILES 241 TO 235, SACRAMENTO RIVER
FUTURE MEANDER BEND MIGRATION AND FLOODPLAIN DEVELOPMENT PATTERNS NEAR RIVER MILES 241 TO 235, SACRAMENTO RIVER Eric W. Larsen University of California, Davis With the assistance of Evan Girvetz REPORT
More informationEvaluation of Subsurface Formation of Pabna District, Bangladesh
IOSR Journal of Applied Geology and Geophysics (IOSR-JAGG) e-issn: 2321 0990, p-issn: 2321 0982.Volume 1, Issue 4 (Sep. Oct. 2013), PP 30-36 Evaluation of Subsurface Formation of Pabna District, Bangladesh
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 informationBoreholes. Implementation. Boring. Boreholes may be excavated by one of these methods: 1. Auger Boring 2. Wash Boring 3.
Implementation Boreholes 1. Auger Boring 2. Wash Boring 3. Rotary Drilling Boring Boreholes may be excavated by one of these methods: 4. Percussion Drilling The right choice of method depends on: Ground
More informationFUTURE MEANDER BEND MIGRATION AND FLOODPLAIN DEVELOPMENT PATTERNS NEAR RIVER MILES 200 TO 191 OF THE SACRAMENTO RIVER PHASE III REPORT
FUTURE MEANDER BEND MIGRATION AND FLOODPLAIN DEVELOPMENT PATTERNS NEAR RIVER MILES 200 TO 191 OF THE SACRAMENTO RIVER PHASE III REPORT Eric W. Larsen REPORT FOR DUCKS UNLIMITED March 31, 2006-1 - Contents
More informationModeling of Atmospheric Effects on InSAR Measurements With the Method of Stochastic Simulation
Modeling of Atmospheric Effects on InSAR Measurements With the Method of Stochastic Simulation Z. W. LI, X. L. DING Department of Land Surveying and Geo-Informatics, Hong Kong Polytechnic University, Hung
More informationScholars Research Library
Available online at www.scholarsresearchlibrary.com Scholars Research Library Archives of Physics Research, 2010, 1 (2):37-45 (http://scholarsresearchlibrary.com/archive.html) ISSN 0976-0970 2-D Resistivity
More informationIPMO2-1. Groundwater Modelling of Chiang Rai Basin, Northern Thailand. Sattaya Intanum* Dr.Schradh Saenton**
IPMO2-1 Groundwater Modelling of Chiang Rai Basin, Northern Thailand Sattaya Intanum* Dr.Schradh Saenton** ABSTRACT Chiang Rai basin, situated in Chiang Rai and Phayao provinces covering an area of 11,000
More informationInternational Workshop on Groundwater Systems in Europe, August, 2013 Berlin
Hydrogeological Mapping in Albania From the IHME contribution to larger scale national maps R. Eftimi, I. Tafilaj, G. Bisha & Xh. Sheganaku Mali me Gropa karst plateau Selita spring The hydrogeological
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 informationA NOVEL APPROACH TO GROUNDWATER MODEL DEVELOPMENT. Thomas D. Krom 1 and Richard Lane 2
A NOVEL APPROACH TO GROUNDWATER MODEL DEVELOPMENT Thomas D. Krom 1 and Richard Lane 2 1 Touch Water Ltd., P.O. Box 143, Lincoln, Christchurch, New Zealand; email: touchwater@gmail.com 2 Applied Research
More informationSummary. Introduction. Theory and/or Method
Case studies from Fort Mc Murray, Horn River basin and Manitoba displaying advances in technology, new approach and updated interpretation of AEM data for unconventional hydrocarbons and groundwater mapping
More informationFLUVIAL LANDFORMS. Floodplains
FLUVIAL LANDFORMS Floodplains fairly flat & continuous surface occupying much of a valley bottom normally underlain by unconsolidated sediments subject to periodic flooding (usually once every year or
More informationAppendix B: Geophysical Data (Thesis Appendix, 2013)
Utah State University From the SelectedWorks of David J Richey 2013 Appendix B: Geophysical Data (Thesis Appendix, 2013) David J Richey, Utah State University Available at: https://works.bepress.com/david_richey/2/
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 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 informationGEOL 652. Poudre River Fieldtrip
GEOL 652. Poudre River Fieldtrip One of the more difficult variables to measure and/or estimate when studying flow in natural channels is that of roughness. Roughness, usually approximated with Manning
More informationLocation of protection zones along production galleries: an example of methodology
Tracers and Modelling in Hydrogeology (Proceedings of the TraM'2000 Conference held at Liège, Belgium, May 2000). IAHS Publ. no. 262, 2000. 141 Location of protection zones along production galleries:
More informationRESISTIVITY IMAGING AND BOREHOLE INVESTIGATION OF THE BANTING AREA AQUIFER, SELANGOR, MALAYSIA. A.N. Ibrahim Z.Z.T. Harith M.N.M.
JOURNAL OF ENVIRONMENTAL HYDROLOGY The Electronic Journal of the International Association for Environmental Hydrology On the World Wide Web at http://www.hydroweb.com VOLUME 11 2003 RESISTIVITY IMAGING
More informationRIVERS, GROUNDWATER, AND GLACIERS
RIVERS, GROUNDWATER, AND GLACIERS Delta A fan-shaped deposit that forms when a river flows into a quiet or large body of water, such as a lake, an ocean, or an inland sea. Alluvial Fan A sloping triangle
More informationEstimating the hydrocarbon volume from elastic and resistivity data: A concept
INTERPRETER S CORNER Coordinated by Rebecca B. Latimer Estimating the hydrocarbon volume from elastic and resistivity data: A concept CARMEN T. GOMEZ, JACK DVORKIN, and GARY MAVKO, Stanford University,
More informationKick-off Meeting ENIGMA ITN
Kick-off Meeting ENIGMA ITN 30/01/2017, Paris General Presentation This project has received funding from European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie
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 information