Parallel Simulation of Subsurface Fluid Flow
|
|
- Camron Floyd
- 5 years ago
- Views:
Transcription
1 Parallel Simulation of Subsurface Fluid Flow Scientific Achievement A new mortar domain decomposition method was devised to compute accurate velocities of underground fluids efficiently using massively parallel computers. Significance and Impact The sheer size and heterogeneous nature of a subsurface reservoir makes it difficult to compute accurate fluid flow velocities, which are needed to simulate, e.g., the response of a geosystem to injection of CO 2. Research Details Rock heterogeneity, i.e., fine-scale variation in permeability, must be resolved on a fine mesh to compute accurate flow velocities. The reservoir is decomposed into subdomains. Each subdomain can be simulated efficiently in parallel. They are tied together through subdomain mortar interfaces. Homogenization theory predicts fine-scale variations using only a few quantities per interface, tying the subdomains together by to a small and efficiently computed mortar problem. Fluid flow can be solved accurately over a huge reservoir volume to a fine resolution efficiently by a massively parallel computer. Fine-scale velocity grid Interface = 1040 New method 3 11 subdomains Interface = 156 T. Arbogast & H. Xiao, A multiscale mortar mixed space based on homogenization for heterogeneous elliptic problems, SIAM J. Numer. Anal., T. Arbogast, Z. Tao & H. Xiao, Multiscale mortar mixed methods for heterogeneous elliptic problems, Contemp. Math., J. Li &H. Yang, eds., Amer. Math. Soc., 2013, Chapter 2. T. Arbogast, Mixed multiscale methods for heterogeneous elliptic problems, in Numerical Analysis of Multiscale Problems, I. G. Graham et al., eds., Springer, 2011, pp
2 Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy s National Nuclear Security Administration under contract DE-AC04-94AL85000.
3 Spanning the Scales in CO 2 Sequestration Simulation Basin (100km) 3 Molecular (1nm) Pore (10μm) Geosystems are large and highly heterogeneous. Simulation requires a fine computational mesh and multiscale models. Field (100m) Continuum (10cm) J. A. Bellian, BEG
4 To simulate geosystem response, we must compute accurate fluid flow velocities. 4 Basin (100km) Size & heterogeneity: 1m grid spacing? Basin: 300km 300km 10 3 m = grid-blocks? Field: 10km 10km 10 2 m = or 10 8 at 10m 10m 1m resolution Additional computational load: Continuum(1m) Multiphase and multicomponent Long time simulations (1000 years) Uncertainty quantification Need parallel computation and multiscale techniques!
5 The Multiblock Concept Divide the domain into nonoverlapping subdomains and simulate the relevant processes locally as appropriate. Tie these together through a mortar interface. 5 Compositional model scco 2 oil Fully implicit two-phase IMPES two-phase brine scco 2 Non-matching grids Faults brine Faults
6 6
7 7 Subdomain blocks for parallel processing
8 8 Fine grid for resolution: Many unknowns!
9 9 Coarse mortar grid: Few unknowns.
10 10 p Guess mortar pressure p.
11 11 p Step 1: Compute the pressure and fluid velocity on each block efficiently in parallel. Based on velocity mismatch on interface, update p. Efficient due to small size.
12 12 p Step 2: Compute the pressure and fluid velocity on each block efficiently in parallel. Based on velocity mismatch on interface, update p. Efficient due to small size.
13 13 p Step 3: Compute the pressure and fluid velocity on each block efficiently in parallel. Done: Negligible velocity mismatch.
14 14 Pressure and velocity computed efficiently in parallel.
15 15 Mortars based on Homogenization Theory Difficulty. Rock heterogeneity, i.e., fine-scale variation in permeability, leads to fine-scale variation in pressure and velocity. Problem. Polynomials do not approximate rough functions well. Idea. Homogenization theory tells us that p (1 + ω )p 0 where p 0 is smooth. So approximate λ (1 + ω )q where q(x) is a polynomial. p p 0 x
16 16 Numerical Test Results Fluvial Permeability (channelized) Fine-scale fluid velocity grid Interface = 1040 Homogenized mortar 3 11 subdomains Interface = 156 Error = 11%
Coupling of Multi fidelity Models Applications to PNP cdft and local nonlocal Poisson equations
Coupling of Multi fidelity Models Applications to PNP cdft and local nonlocal Poisson equations P. Bochev, J. Cheung, M. D Elia, A. Frishknecht, K. Kim, M. Parks, M. Perego CM4 summer school, Stanford,
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part II Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationOn some numerical convergence studies of mixed finite element methods for flow in porous media
On some numerical convergence studies of mixed finite element methods for flow in porous media Gergina Pencheva Abstract We consider an expanded mixed finite element method for solving second-order elliptic
More informationUtilizing Adjoint-Based Techniques to Improve the Accuracy and Reliability in Uncertainty Quantification
Utilizing Adjoint-Based Techniques to Improve the Accuracy and Reliability in Uncertainty Quantification Tim Wildey Sandia National Laboratories Center for Computing Research (CCR) Collaborators: E. Cyr,
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part IV Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationT. Arbogast, Zhen Tao, and Hailong Xiao, Multiscale mortar mixed methods for heterogeneous elliptic problems, in Recent Advances in Scientific
T. Arbogast, Zhen Tao, and Hailong Xiao, Multiscale mortar mixed methods for heterogeneous elliptic problems, in Recent Advances in Scientific Computing and Applications, Jichun Li, Hongtao Yang, and Eric
More informationA MULTISCALE METHOD FOR MODELING TRANSPORT IN POROUS MEDIA ON UNSTRUCTURED CORNER-POINT GRIDS
A MULTISCALE METHOD FOR MODELING TRANSPORT IN POROUS MEDIA ON UNSTRUCTURED CORNER-POINT GRIDS JØRG E. AARNES AND YALCHIN EFENDIEV Abstract. methods are currently under active investigation for the simulation
More informationNonlinearity of Joints in Structural Dynamics of Weapons Systems
Nonlinearity of Joints in Structural Dynamics of Weapons Systems Daniel Brown AWE Dan Segalman Sandia National Laboratories Sandia National Laboratories is a multi-program laboratory managed and operated
More informationPhysical and Computational Domain Decompositions for Modeling Subsurface Flows
Contemporary Mathematics Physical and Computational Domain Decompositions for Modeling Subsurface Flows Mary F. Wheeler and Ivan Yotov 1. Introduction Modeling of multiphase flow in permeable media plays
More informationOutline: 1 Motivation: Domain Decomposition Method 2 3 4
Multiscale Basis Functions for Iterative Domain Decomposition Procedures A. Francisco 1, V. Ginting 2, F. Pereira 3 and J. Rigelo 2 1 Department Mechanical Engineering Federal Fluminense University, Volta
More informationAlgebraic Coarse Spaces for Overlapping Schwarz Preconditioners
Algebraic Coarse Spaces for Overlapping Schwarz Preconditioners 17 th International Conference on Domain Decomposition Methods St. Wolfgang/Strobl, Austria July 3-7, 2006 Clark R. Dohrmann Sandia National
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationA Multiscale Mortar Method And Two-Stage Preconditioner For Multiphase Flow Using A Global Jacobian Approach
SPE-172990-MS A Multiscale Mortar Method And Two-Stage Preconditioner For Multiphase Flow Using A Global Jacobian Approach Benjamin Ganis, Kundan Kumar, and Gergina Pencheva, SPE; Mary F. Wheeler, The
More informationAn overview of the. package for non-conformal mesh tying or simple contact problems
An overview of the package for non-conformal mesh tying or simple contact problems Glen Hansen 2012 Trilinos User Group Meeting Oct. 30, 2012 SAND 2012-9202P (UUR) Sandia National Laboratories is a multi-program
More informationMultiscale Finite Element Methods. Theory and
Yalchin Efendiev and Thomas Y. Hou Multiscale Finite Element Methods. Theory and Applications. Multiscale Finite Element Methods. Theory and Applications. October 14, 2008 Springer Dedicated to my parents,
More informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationMiguel Aguilo, Ted Blacker, Andre Claudet, Brett Clark, Ryan Rickerson, Josh Robbins, Louis Vaught, and Tom Voth Sandia National Laboratories
Miguel Aguilo, Ted Blacker, Andre Claudet, Brett Clark, Ryan Rickerson, Josh Robbins, Louis Vaught, and Tom Voth (as told by Corbett Battaile) Sandia National Laboratories Sandia National Laboratories
More informationModeling Interfaces in Structural Dynamics Analysis
Modeling Interfaces in Structural Dynamics Analysis Enriching our joints models and capabilities 2nd Workshop on Joints Modelling Dartington Hall Totnes, Devon, UK 26-29 April 2009 Michael Starr, Org.
More informationHydrogeology of Deep Borehole Disposal for High-Level Radioactive Waste
SAND2014-18615C Hydrogeology of Deep Borehole Disposal for High-Level Radioactive Waste Geological Society of America Annual Meeting October 20, 2014 Bill W. Arnold, W. Payton Gardner, and Patrick V. Brady
More informationOptimal Interface Conditions for an Arbitrary Decomposition into Subdomains
Optimal Interface Conditions for an Arbitrary Decomposition into Subdomains Martin J. Gander and Felix Kwok Section de mathématiques, Université de Genève, Geneva CH-1211, Switzerland, Martin.Gander@unige.ch;
More informationMultiphysics modeling of thermal batteries
Multiphysics modeling of thermal batteries Scott A. Roberts, Ph.D. Thermal/Fluid Component Sciences Department Sandia National Laboratories, Albuquerque, NM The Future of Munitions Batteries Workshop Army
More informationScalability of Partial Differential Equations Preconditioner Resilient to Soft and Hard Faults
Scalability of Partial Differential Equations Preconditioner Resilient to Soft and Hard Faults K.Morris, F.Rizzi, K.Sargsyan, K.Dahlgren, P.Mycek, C.Safta, O.LeMaitre, O.Knio, B.Debusschere Sandia National
More informationMMsFEM and Streamlines
MMsFEM and Streamlines Combining a mixed multiscale FEM with streamline simulation for enhanced reservoir performance prediction on large grid models. Jørg E. Aarnes, Vegard Kippe and Knut-Andreas Lie
More informationSENSITIVITY ANALYSIS OF THE PETROPHYSICAL PROPERTIES VARIATIONS ON THE SEISMIC RESPONSE OF A CO2 STORAGE SITE. Juan E. Santos
SENSITIVITY ANALYSIS OF THE PETROPHYSICAL PROPERTIES VARIATIONS ON THE SEISMIC RESPONSE OF A CO2 STORAGE SITE Juan E. Santos Instituto del Gas y del Petróleo, Facultad de Ingeniería UBA and Department
More informationNonparametric density estimation for elliptic problems with random perturbations
Nonparametric density estimation for elliptic problems with random perturbations, DqF Workshop, Stockholm, Sweden, 28--2 p. /2 Nonparametric density estimation for elliptic problems with random perturbations
More informationUncertainty Quantification for Machine Learning and Statistical Models
Uncertainty Quantification for Machine Learning and Statistical Models David J. Stracuzzi Joint work with: Max Chen, Michael Darling, Stephen Dauphin, Matt Peterson, and Chris Young Sandia National Laboratories
More informationPhotos placed in horizontal p osition with even amount of white sp ace between photos and header
Photos placed in h orizontal position with even amount of white space between photos a nd header Interbed Modeling to Predict Wellbore Damage for ig Hill Strategic Petroleum Reserve Photos placed in horizontal
More informationConvergence Behavior of a Two-Level Optimized Schwarz Preconditioner
Convergence Behavior of a Two-Level Optimized Schwarz Preconditioner Olivier Dubois 1 and Martin J. Gander 2 1 IMA, University of Minnesota, 207 Church St. SE, Minneapolis, MN 55455 dubois@ima.umn.edu
More informationSierra/SD: A Massively Parallel Finite Element Code for Structural Dynamics and Acoustics Analysis
Sierra/SD: A Massively Parallel Finite Element Code for Structural Dynamics and Acoustics Analysis Jerry W. Rouse, Timothy F. Walsh, Greg D. Tipton, Garth M. Reese Sandia National Labs Albuquerque, NM
More informationMixed Multiscale Methods for Heterogeneous Elliptic Problems
Mixed Multiscale Methods for Heterogeneous Elliptic Problems Todd Arbogast Abstract We consider a second order elliptic problem written in mixed form, i.e., as a system of two first order equations. Such
More informationModeling of two-phase flow in fractured porous media on unstructured non-uniform coarse grids
Modeling of two-phase flow in fractured porous media on unstructured non-uniform coarse grids Jørg Espen Aarnes and Vera Louise Hauge SINTEF ICT, Deptartment of Applied Mathematics Applied Mathematics
More informationMultiscale Methods for Multiphase Flow in Porous Media
Multiscale Methods for Multiphase Flow in Porous Media Jan M. Nordbotten Department of Mathematics, University of Bergenjan.nordbotten@math.uib.no Department of Civil and Environmental Engineering, Princeton
More informationTexas at Austin 2 Nazarbayev University. January 15, 2019
Recovery of the Interface Velocity for the Incompressible Flow in Enhanced Velocity Mixed Finite Element Method arxiv:1901.04401v1 [math.na] 14 Jan 2019 Yerlan Amanbek 1,2, Gurpreet Singh 1, and Mary F.
More informationFracture-Matrix Flow Partitioning and Cross Flow: Numerical Modeling of Laboratory Fractured Core Flood
Fracture-Matrix Flow Partitioning and Cross Flow: Numerical Modeling of Laboratory Fractured Core Flood R. Sanaee *, G. F. Oluyemi, M. Hossain, and M. B. Oyeneyin Robert Gordon University *Corresponding
More informationA Θ(n) Approximation Algorithm for 2-Dimensional Vector Packing
A Θ(n) Approximation Algorithm for 2-Dimensional Vector Packing Ekow Otoo a, Ali Pinar b,, Doron Rotem a a Lawrence Berkeley National Laboratory b Sandia National Laboratories Abstract We study the 2-dimensional
More informationLocal Time Step for a Finite Volume Scheme I.Faille F.Nataf*, F.Willien, S.Wolf**
Controlled CO 2 Diversified fuels Fuel-efficient vehicles Clean refining Extended reserves Local Time Step for a Finite Volume Scheme I.Faille F.Nataf*, F.Willien, S.Wolf** *: Laboratoire J.L.Lions **:Université
More informationICES REPORT Adaptive Numerical Homogenization for Non-Linear Multiphase Flow and Transport
ICES REPORT 17-13 June 217 Adaptive Numerical Homogenization for Non-Linear Multiphase Flow and Transport by Gurpreet Singh, Yerlan Amanbek, and Mary F. Wheeler The Institute for Computational Engineering
More informationMultiphysics Modeling of Pre-ignition Damage in Energetic Materials and the Effect on Cookoff Violence
Multiphysics Modeling of Pre-ignition Damage in Energetic Materials and the Effect on Cookoff Violence Dan Turner dzturne@sandia.gov Department of Thermal & Fluid Mechanics Sandia National Laboratories
More informationCharacterization of wind velocity distributions within a full-scale heliostat field
Photos placed in horizontal position with even amount of white space between photos and header Characterization of wind velocity distributions within a full-scale heliostat field Jeremy Sment, Graduate
More informationLocal Mesh Refinement with the PCD Method
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 1, pp. 125 136 (2013) http://campus.mst.edu/adsa Local Mesh Refinement with the PCD Method Ahmed Tahiri Université Med Premier
More informationAn Empirical Comparison of Graph Laplacian Solvers
An Empirical Comparison of Graph Laplacian Solvers Kevin Deweese 1 Erik Boman 2 John Gilbert 1 1 Department of Computer Science University of California, Santa Barbara 2 Scalable Algorithms Department
More informationA FRONT-TRACKING METHOD FOR HYPERBOLIC THREE-PHASE MODELS
A FRONT-TRACKING METHOD FOR HYPERBOLIC THREE-PHASE MODELS Ruben Juanes 1 and Knut-Andreas Lie 2 1 Stanford University, Dept. Petroleum Engineering, USA 2 SINTEF IKT, Dept., Norway ECMOR IX, August 30 September
More informationNumerical Simulation of Flows in Highly Heterogeneous Porous Media
Numerical Simulation of Flows in Highly Heterogeneous Porous Media R. Lazarov, Y. Efendiev, J. Galvis, K. Shi, J. Willems The Second International Conference on Engineering and Computational Mathematics
More informationDomain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions
Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions Ernst P. Stephan 1, Matthias Maischak 2, and Thanh Tran 3 1 Institut für Angewandte Mathematik, Leibniz
More informationA MATRIX ANALYSIS OF OPERATOR-BASED UPSCALING FOR THE WAVE EQUATION
SIAM J. NUMER. ANAL. Vol. 44, No. 2, pp. 586 62 c 2006 Society for Industrial and Applied Mathematics A MATRIX ANALYSIS OF OPERATOR-BASED UPSCALING FOR THE WAVE EQUATION OKSANA KOROSTYSHEVSKAYA AND SUSAN
More informationUnravel Faults on Seismic Migration Images Using Structure-Oriented, Fault-Preserving and Nonlinear Anisotropic Diffusion Filtering
PROCEEDINGS, 44th Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 2019 SGP-TR-214 Unravel Faults on Seismic Migration Images Using Structure-Oriented,
More informationMultiscale method and pseudospectral simulations for linear viscoelastic incompressible flows
Interaction and Multiscale Mechanics, Vol. 5, No. 1 (2012) 27-40 27 Multiscale method and pseudospectral simulations for linear viscoelastic incompressible flows Ling Zhang and Jie Ouyang* Department of
More informationMultiscale Modeling and Simulations of Flows in Naturally Fractured Karst Reservoirs
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 6, No. 1, pp. 162-184 Commun. Comput. Phys. July 2009 Multiscale Modeling and Simulations of Flows in Naturally Fractured Karst Reservoirs Peter Popov 1,, Yalchin
More informationThe Deflation Accelerated Schwarz Method for CFD
The Deflation Accelerated Schwarz Method for CFD J. Verkaik 1, C. Vuik 2,, B.D. Paarhuis 1, and A. Twerda 1 1 TNO Science and Industry, Stieltjesweg 1, P.O. Box 155, 2600 AD Delft, The Netherlands 2 Delft
More informationAvailable online at ScienceDirect
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 114 (17) 3343 3352 www.elsevier.com/locate/procedia 13th International Conference on Greenhouse Gas Control Technologies, GHGT-13,
More informationWESTCARB Regional Partnership
WESTCARB Regional Partnership Subsurface Flow Modeling at King Island Christine Doughty, CADoughty@lbl.gov Curtis Oldenburg, CMOldenburg@lbl.gov Staff Scientists Lawrence Berkeley National Laboratory WESTCARB
More informationDSMC-Based Shear-Stress/Velocity-Slip Boundary Condition for Navier-Stokes Couette-Flow Simulations
DSMC-Based Shear-Stress/Velocity-Slip Boundary Condition for Navier-Stokes Couette-Flow Simulations J. R. Torczynski and M. A. Gallis Engineering Sciences Center, Sandia National Laboratories, P. O. Box
More informationarxiv: v1 [math.na] 11 Jul 2011
Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients arxiv:07.260v [math.na] Jul 20 Blanca Ayuso De Dios, Michael Holst 2, Yunrong Zhu 2, and Ludmil Zikatanov
More informationModeling Reactive Flows in Porous Media
Modeling Reactive Flows in Porous Media Peter Lichtner (lead PI), Los Alamos National Laboratory Glenn Hammond, Pacific Northwest National Laboratory Richard Tran Mills, Oak Ridge National Laboratory NCCS
More informationSarah Jane Riordan. Australian School of Petroleum University of Adelaide March 2009
Managing the Interdisciplinary Requirements of 3D Geological Models Sarah Jane Riordan Australian School of Petroleum University of Adelaide March 2009 Thesis submitted in accordance with the requirements
More information20. A Dual-Primal FETI Method for solving Stokes/Navier-Stokes Equations
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 23 DDM.org 2. A Dual-Primal FEI Method for solving Stokes/Navier-Stokes
More informationParallel Tensor Compression for Large-Scale Scientific Data
Parallel Tensor Compression for Large-Scale Scientific Data Woody Austin, Grey Ballard, Tamara G. Kolda April 14, 2016 SIAM Conference on Parallel Processing for Scientific Computing MS 44/52: Parallel
More informationSUMMARY NUMERICAL SIMULATION OF BRINE-CO 2 FLOW AND WAVE PROPAGATION
Sensitivity analysis of the petrophysical properties variations on the seismic response of a CO 2 storage site Juan E. Santos, Universidad de Buenos Aires, Facultad Ingeniería, Instituto del Gas y del
More informationA Two-grid Method for Coupled Free Flow with Porous Media Flow
A Two-grid Method for Coupled Free Flow with Porous Media Flow Prince Chidyagwai a and Béatrice Rivière a, a Department of Computational and Applied Mathematics, Rice University, 600 Main Street, Houston,
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 informationMortar upscaling for multiphase flow in porous media
Computational Geosciences 6: 73 100, 2002. 2002 Kluwer Academic Publishers. Printed in the Netherlands. Mortar upscaling for multiphase flow in porous media Małgorzata Peszyńska a, Mary F. Wheeler a andivanyotov
More informationNumerical and Laboratory Study of Gas Flow through Unconventional Reservoir Rocks
Numerical and Laboratory Study of Gas Flow through Unconventional Reservoir Rocks RPSEA Piceance Basin Tight Gas Research Review Xiaolong Yin, Assistant Professor Petroleum Engineering, Colorado School
More informationA FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS
Proceedings of ALGORITMY 2005 pp. 222 229 A FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS ELENA BRAVERMAN, MOSHE ISRAELI, AND ALEXANDER SHERMAN Abstract. Based on a fast subtractional
More informationTwo new enriched multiscale coarse spaces for the Additive Average Schwarz method
346 Two new enriched multiscale coarse spaces for the Additive Average Schwarz method Leszek Marcinkowski 1 and Talal Rahman 2 1 Introduction We propose additive Schwarz methods with spectrally enriched
More informationNumerical and Laboratory Study of Gas Flow through Unconventional Reservoir Rocks
Numerical and Laboratory Study of Gas Flow through Unconventional Reservoir Rocks RPSEA Piceance Basin Tight Gas Research Review Xiaolong Yin, Assistant Professor Petroleum Engineering, Colorado School
More informationDownloaded 08/04/16 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 38, No. 4, pp. B593 B618 c 2016 Society for Industrial and Applied Mathematics ACTIVE-SET REDUCED-SPACE METHODS WITH NONLINEAR ELIMINATION FOR TWO-PHASE FLOW PROBLEMS IN POROUS
More informationTWO-SCALE SIMULATION OF MAXWELL S EQUATIONS
ESAIM: PROCEEDINGS, February 27, Vol.6, 2-223 Eric Cancès & Jean-Frédéric Gerbeau, Editors DOI:.5/proc:27 TWO-SCALE SIMULATION OF MAWELL S EQUATIONS Hyam Abboud,Sébastien Jund 2,Stéphanie Salmon 2, Eric
More informationField Scale Modeling of Local Capillary Trapping during CO 2 Injection into the Saline Aquifer. Bo Ren, Larry Lake, Steven Bryant
Field Scale Modeling of Local Capillary Trapping during CO 2 Injection into the Saline Aquifer Bo Ren, Larry Lake, Steven Bryant 2 nd Biennial CO 2 for EOR as CCUS Conference Houston, TX October 4-6, 2015
More informationA Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation
A Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation Tao Zhao 1, Feng-Nan Hwang 2 and Xiao-Chuan Cai 3 Abstract In this paper, we develop an overlapping domain decomposition
More informationPengcheng Fu, Yue Hao, and Charles R. Carrigan
Pengcheng Fu, Yue Hao, and Charles R. Carrigan Math, Science and Computation of Hydraulic Fracturing Workshop Stanford, March 21, 2013 This work was performed under the auspices of the U.S. Department
More informationMéthodes de raffinement espace-temps
Méthodes de raffinement espace-temps Caroline Japhet Université Paris 13 & University of Maryland Laurence Halpern Université Paris 13 Jérémie Szeftel Ecole Normale Supérieure, Paris Pascal Omnes CEA,
More information2 IVAN YOTOV decomposition solvers and preconditioners are described in Section 3. Computational results are given in Section Multibloc formulat
INTERFACE SOLVERS AND PRECONDITIONERS OF DOMAIN DECOMPOSITION TYPE FOR MULTIPHASE FLOW IN MULTIBLOCK POROUS MEDIA IVAN YOTOV Abstract. A multibloc mortar approach to modeling multiphase ow in porous media
More informationANALYSIS OF THE EFFECTS OF HETEROGENEITY ON HEAT EXTRACTION IN AN EGS REPRESENTED WITH THE CONTINUUM FRACTURE MODEL
PROCEEDINGS, Thirty-Seventh Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 30 - February 1, 2012 SGP-TR-194 ANALYSIS OF THE EFFECTS OF HETEROGENEITY ON
More informationNinth International Water Technology Conference, IWTC9 2005, Sharm El-Sheikh, Egypt 673
Ninth International Water Technology Conference, IWTC9 2005, Sharm El-Sheikh, Egypt 673 A NEW NUMERICAL APPROACH FOR THE SOLUTION OF CONTAMINANT TRANSPORT EQUATION Mohamed El-Gamel Department of Mathematical
More informationA010 MULTISCALE RESERVOIR CHARACTERIZATION USING
1 A010 MULTISCALE RESERVOIR CHARACTERIZATION USING RODUCTION AND TIME LASE SEISMIC DATA Mokhles MEZGHANI, Alexandre FORNEL, Valérie LANGLAIS, Nathalie LUCET IF, 1 & 4 av de Bois réau, 92852 RUEIL-MALMAISON
More informationPartial Saturation Fluid Substitution with Partial Saturation
Fluid Substitution with 261 5 4.5 ρ fluid S w ρ w + S o ρ o + S g ρ g Vp (km/s) 4 3.5 K fluid S w K w + S o K o + S g K g Patchy Saturation Drainage 3 2.5 2 Fine-scale mixing 1 = S w + S o + S g K fluid
More informationPreliminary TOUGH2 Model of King Island CO 2 Injection. Modeling Approach
Preliminary TOUGH2 Model of King Island CO 2 Injection Christine Doughty Presented by Curt Oldenburg Earth Sciences Division Lawrence Berkeley National Laboratory October 25, 2011 Modeling Approach Characterization
More informationPeaceman and Thiem well models or how to remove a logarithmic singularity fr. your numerical solution
Peaceman and Thiem well models or how to remove a logarithmic singularity from your numerical solution Department of Mathematics Oregon State University AMC, 3/2/2007 Current work supported by NSF-0511190
More informationMultiscale Computation for Incompressible Flow and Transport Problems
Multiscale Computation for Incompressible Flow and Transport Problems Thomas Y. Hou Applied Mathematics, Caltech Collaborators: Y. Efenidev (TAMU), V. Ginting (Colorado), T. Strinopolous (Caltech), Danping
More informationWP 4.1. Site selection criteria and ranking methodology. Karen Kirk
WP 4.1 Site selection criteria and ranking methodology Karen Kirk 1 Basic site selection criteria Sufficient depth and storage capacity supercritical CO 2 below 700-800 m (rule of thumb) 2 Variation of
More informationM.Sc. Track Petroleum Engineering & Geosciences
M.Sc. Track Petroleum Engineering & Geosciences 3-9-2017 Delft University of Technology Challenge the future Challenging industry (and study) 2 Geothermal Energy Produce energy (heat) from subsurface for
More informationOn Application of the Weak Galerkin Finite Element Method to a Two-Phase Model for Subsurface Flow
J Sci Comput (2016) 66:225 239 DOI 10.1007/s10915-015-0021-8 On Application of the Weak Galerkin Finite lement Method to a Two-Phase Model for Subsurface Flow Victor Ginting 1 Guang Lin 2 Jiangguo Liu
More informationDECOUPLED AND MULTIPHYSICS MODELS FOR NON-ISOTHERMAL COMPOSITIONAL TWO-PHASE FLOW IN POROUS MEDIA
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 9, Number 1, Pages 17 28 c 2012 Institute for Scientific Computing and Information DECOUPLED AND MULTIPHYSICS MODELS FOR NON-ISOTHERMAL COMPOSITIONAL
More informationGenerating Equidistributed Meshes in 2D via Domain Decomposition
Generating Equidistributed Meshes in 2D via Domain Decomposition Ronald D. Haynes and Alexander J. M. Howse 2 Introduction There are many occasions when the use of a uniform spatial grid would be prohibitively
More informationImaging complex structure with crosswell seismic in Jianghan oil field
INTERPRETER S CORNER Coordinated by Rebecca B. Latimer Imaging complex structure with crosswell seismic in Jianghan oil field QICHENG DONG and BRUCE MARION, Z-Seis, Houston, Texas, U.S. JEFF MEYER, Fusion
More informationMaximising the use of publicly available data: porosity and permeability mapping of the Rotliegend Leman Sandstone, Southern North Sea
Maximising the use of publicly available data: porosity and permeability mapping of the Rotliegend Leman Sandstone, Southern North Sea Claire Imrie & Henk Kombrink 09 May 2018 Overview The Oil and Gas
More informationA SIMULATOR WITH NUMERICAL UPSCALING FOR THE ANALYSIS OF COUPLED MULTIPHASE FLOW AND GEOMECHANICS IN HETEROGE-
A SIMULATOR WITH NUMERICAL UPSCALING FOR THE ANALYSIS OF COUPLED MULTIPHASE FLOW AND GEOMECHANICS IN HETEROGE- NEOUS AND DEFORMABLE POROUS AND FRACTURED MEDIA A Dissertation by DAEGIL YANG Submitted to
More informationHomogenization and numerical Upscaling. Unsaturated flow and two-phase flow
Homogenization and numerical Upscaling Unsaturated flow and two-phase flow Insa Neuweiler Institute of Hydromechanics, University of Stuttgart Outline Block 1: Introduction and Repetition Homogenization
More informationApplications of Partial Differential Equations in Reservoir Simulation
P-32 Applications of Partial Differential Equations in Reservoir Simulation Deepak Singh Summary The solution to stochastic partial differential equations may be viewed in several manners. One can view
More informationChapter 2 Overview of Processes Occurring During CO 2 Geological Storage and Their Relevance to Key Questions of Performance
Chapter 2 Overview of Processes Occurring During CO 2 Geological Storage and Their Relevance to Key Questions of Performance Chin-Fu Tsang and Auli Niemi Abstract The objective of this chapter is to provide
More informationDemystifying Tight-gas Reservoirs using Multi-scale Seismic Data
Demystifying Tight-gas Reservoirs using Multi-scale Seismic Data Overview Murray Roth* Transform Software and Services, Inc., Littleton, CO murray@transformsw.com Tom Davis Colorado School of Mines, Golden,
More informationA FAMILY OF MULTISCALE HYBRID-MIXED FINITE ELEMENT METHODS FOR THE DARCY EQUATION WITH ROUGH COEFFICIENTS. 1. Introduction
A FAMILY OF MULTISCALE HYBRID-MIXED FINITE ELEMENT METHODS FOR THE DARCY EQUATION WITH ROUGH COEFFICIENTS CHRISTOPHER HARDER, DIEGO PAREDES 2, AND FRÉDÉRIC VALENTIN Abstract. We aim at proposing novel
More informationModeling pressure response into a fractured zone of Precambrian basement to understand deep induced-earthquake hypocenters from shallow injection
Modeling pressure response into a fractured zone of Precambrian basement to understand deep induced-earthquake hypocenters from shallow injection S. Raziperchikolaee 1 and J. F. Miller 1 Abstract Analysis
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationINTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR GRIDS. 1. Introduction
Trends in Mathematics Information Center for Mathematical Sciences Volume 9 Number 2 December 2006 Pages 0 INTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR
More informationAnomaly Density Estimation from Strip Transect Data: Pueblo of Isleta Example
Anomaly Density Estimation from Strip Transect Data: Pueblo of Isleta Example Sean A. McKenna, Sandia National Laboratories Brent Pulsipher, Pacific Northwest National Laboratory May 5 Distribution Statement
More informationGeometric Multigrid Methods
Geometric Multigrid Methods Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University IMA Tutorial: Fast Solution Techniques November 28, 2010 Ideas
More informationUpscaling Wave Computations
Upscaling Wave Computations Xin Wang 2010 TRIP Annual Meeting 2 Outline 1 Motivation of Numerical Upscaling 2 Overview of Upscaling Methods 3 Future Plan 3 Wave Equations scalar variable density acoustic
More informationAN UNFITTED METHOD FOR TWO-PHASE FLOW IN FRACTURED POROUS MEDIA
XIX International Conference on Water Resources CMWR 2012 University of Illinois at Urbana-Champaign June 17-22,2012 AN UNFIED MEHOD FOR WO-PHASE FLOW IN FRACURED POROUS MEDIA A. Fumagalli and A. Scotti
More information