Simulating Fluid-Fluid Interfacial Area
|
|
- Branden Skinner
- 6 years ago
- Views:
Transcription
1 Simulating Fluid-Fluid Interfacial Area revealed by a pore-network model V. Joekar-Niasar S. M. Hassanizadeh Utrecht University, The Netherlands July 22, 2009
2 Outline 1 What s a Porous medium 2 Intro to PNMs 3 PNMing technique 4 Two-phase PNMs 5 Applications Case1: idealized porous medium Case2: micromodel Case3: glass-bead column
3 A porous medium at pore scale: How does it look like? Glass beads vs. volcanic Tuff courtesy of D. Wildenschild
4 A porous medium at pore scale: How does it look like? Sandstone courtesy of O. Vizika
5 A porous medium at pore scale: How does it look like? Fibers
6 Characterization of porous media structure Microscopic: Topology, geometry Macroscopic: Porosity, descriptive topology functions, geometry statistics, spatial distribution
7 Quasti-static vs. dynamic PNMs Quasi-static Computationally, very cheap. No pressure field is solved. Pore-scale geometry and morphology are only important. Used extensively, for two-phase and three-phase flow; P c S w, k r S w, S w a nw, reactive transport, etc. Can be used as a predictive tool. Dynamic Computationally, expensive. Pressure field is solved. Network and fluids properties are important. Notbeenusedasextensivelyas quasi-static ones; P c S w, k r S w, ganglia movement, dynamic pressure field Alongwaytogo!
8 Investigation tool Pore Network Models? Yes,because: Physical-based models using pore-scale information. Application for static and dynamic processes. Computationally, not so expensive (especially for quasi-static ones). Capability to provide up-scaled information. But: Translation of topology and geometry is inevitable, and not always straight forward! No detail information within a pore (e.g. pressure field). Local laws/rules are inevitable, and devil! be careful!!!!
9 Definitions Geometry Pore body: large pores among the grains (shape and size) Pore throat: long narrow pores connecting the large pores (shape and size) Topology Coordination number: number of pore bodies connected to a given pore body Fixed coordination number: isotropic, variable coordination number: anisotropic
10 Pore-network modeling steps Step 1. Topological and geometrical data acquirement Step 2. Pore-network structure Step 3. Defining governing equations and local laws Step 4. Numerical experiments
11 Step 1. Topological and geometrical data acquirement Step 1.a Topological data acquirement Detailed network acquired by analyzing images of a real porous medium, e.g. 3DMA-Rock software. Hypothetical to idealized pore-networks. They can be isotropic (fixed coordination number) or anisotropic (variable coordination number). Center of pore bodies can be located on lattice points (structured) or random coordinates (unstructured). Step 1.b Geometrical data acquirement Detailed network acquired by analyzing images of a real porous medium, e.g. 3DMA-Rock software. Statistically equivalent or idealized porous medium. Generating based on statistical parameters with/without spatial correlations.
12 Step 2. Pore-network structure Numbering pore bodies, sequentially. Numbering pore throats, sequentially. Defining connectivity of pore bodies to pore throats (Pb2Pth). Defining connectivity of pore throats to pore bodies (Pth2Pb).
13 Step 2. Pore-network structure Void-Space Network (Pore-Network)
14 Step 2. Pore-network structure Pore Bodies Numbering
15 Step 2. Pore-network structure Pore Throats Numbering
16 Step 2. Pore-network structure Connection Matrices
17 Step 3. Governing equations and local laws One-phase or multiphase flow? Drainage or Imbibition process,...? Dynamic process or static process? Probable physical laws: Young-Laplace Equation, Hagen-Poiseuille Equation, Stokes Equation, Conservation, etc. Probable local assumptions: phase trapping assumptions, hydraulic conductivites of pores, etc.
18 Step 4. Numerical experiments Quasi-Static or dynamic Quasi-static simulations: capillary forces are dominant, no information about time scale of phenomena. Dynamic simulations: solving pressure field, change of saturation with time.
19 Drainage vs. imbibition Let s watch some movies: courtesy of Fabiano G. Wolf ( fgwolf/) Drainage First, larger pore will be filled by the nonwetting phase (blue). Imbibition First, smaller pore will be filled by the wetting phase (red). Simulations done by Lattice-Gas Automata Models
20 Example Quasi-static simulation, drainage process, circular cross section Entry P c, Pi c = 2σ r ij cos θ (for circular cross section). Lower entry pressure: invaded earlier during drainage!
21 Quasi-static PNMs Network: Initial Condition
22 Quasi-static PNMs
23 Quasi-static PNMs
24 Quasi-static example
25 Quasi-static example
26 Quasi-static example
27 Quasi-static example
28 Quasi-static example
29 Quasi-static example P c =2.0, with trapping
30 Quasi-static example movie
31 Computational considerations Adaptive pressure update. Gains: a) Calculating the exact P c -S w curve. figure b) Saving computational time especially at the initial and residual saturation of P c -S w curve. Using flags pore those pores which are in contact with the invading phase, but not filled yet. For strongly anisotropic: save Pb2Pth matrix differently to save memory. The best known suggest approach is compact row storage (CRS). Trapping a the most time-consuming. Robust search algorithm to save time such as depth-first search (Cormen, Leiserson, Rivest, 2000).
32 Computational considerations Compact Row Storage (CRS) for an anisotropic medium or dual porosity medium! Example: dual porosity medium; macro pores connected to micro pores. courtesy of M. Prodanovic
33 Computational considerations Compact Row Storage (CRS) Yellow colors are indices!
34 Computational considerations Compact Row Storage (CRS), if pore throats are included in calculations!
35 Computational considerations Compact Row Storage (CRS), if pore throats are included in calculations!
36 Computational considerations Compact Row Storage (CRS), if pore throats are not included in calculations!
37 Computational considerations Compact Row Storage (CRS), if pore throats are not included in calculations!
38 Geometrical parameters/issues Geometry Cross section geometry Aspect ratio Shape factor Phenomena Piston-like movement Snap-off Capillary diffusion Cooperative filling,...
39 Geometrical evolution??? To have corners, or not to have corners: that is the question! To mimic the grain geometry,sometimes this can be also the question!! Wetting phase in the corners will reduce wetting phase trapping during drainage significantly, however may contribute to residual saturation. Wetting phase in the corners can increase snap-off during imbibition. Snap-off happens if wetting film/ corner flow cuts connectivity of the nonwetting fluid. Including corner flow in dynamic models will increase numerical complexity enormously.
40 Geometrical parameters: cross section geometry (Joekar-Niasar et al,2009) Entry capillary pressure can be calculated for different cross sections using MS-P (Mayer, Stowe and Princen, J. Colloid Interface Sci, 1968, 1969; Ma et al., Colloids Surf, 1996). More complex geometry, more complicated system of equations.
41 Geometrical parameters: shape factor κ = A/P 2 κ: shape factor, A:cross section area, P:cross section perimeter Pore-network models for glass beads (Joekar-Niasar et al,2009)
42 Geometrical parameters: aspect ratio Radius of pore body size to pore throat size Significant effect during imbibition: larger aspect ratio, more snap-off! More snap-off, more nonwetting phase trapping (Lenormand and Zarcone, 1983, 1984; Wardlaw and Yu, 1988; Ioannidis et al. 1991) movie.
43 Geometrical parameters: aspect ratio Overlapping of pore bodies and pore throats size distribution.
44 A porous medium at pore scale: How does it look like? Glass beads vs. volcanic Tuff courtesy of D. Wildenschild
45 A porous medium at pore scale: How does it look like? Sandstone courtesy of O. Vizika
46 A porous medium at pore scale: How does it look like? Fibers
47 What we want: We are interested in investigating extended theories of two-phase flow.
48 Objectives 1. Insights into pore-scale dynamics of two-phase flow 2. Studying validity of extended two-phase flow equations
49 Extended theories of two-phase flow n S α + q α = 0 t q α = K α ( p α ρ α g σ αa a wn σ αs S α) a wn t + (a wn w wn )=r wn (a wn, S w ) w wn = K wn (a wn γ wn ) p n p w = p c τ S w t, pc = f (S w, a wn )
50 Extended theories of two-phase flow in porous media Capillary pressure p c = f (S w, a wn )=(p n p w ) equil Difference in fluid s pressure are equal to static capillary pressure but only at equilibrium. Under non-equilibrium conditions, assuming linearity: p n p w = p c τ S w t, pc = f (S w, a wn )
51 Case1: idealized porous medium Joekar-Niasar, V., Hassanizadeh, S.M., Leijnse A., 2008, Insights into the relationships among capillary pressure, saturation, interfacial area and relative permeability using pore-network modeling, Transport in Porous Media, 74:201219, doi /s
52 Idealized pore-network Specifications Topology: structured, isotropic, lattice network with coordination number of 6 Geometry: sphere pore bodies, circular cross-section pore throats Governing equation: P c = 2σ r cos θ TiPM Joekar-Niasar et al. 2008,
53 Filling process Drainage: Starting from the largest pore throats. If a pore throat is invaded, the connected pore body will be fully filled. Imbibition: Starting from the smallest pore throats. If a pore throat is invaded, the connected pore body will be partially filled.
54 P c -S w and a nw -S w curves: imbibition
55 P c -S w and a nw -S w curves: drainage
56 What s a Porous medium Intro to PNMs PNMing technique Two-phase PNMs Applications P c -S w -anw surfaces Drainage vs. imbibition a y b y y
57 Relative permeabilities Flux passing through the pore throat ij based on Hagen-Poiseuille s law q ij = πr ij 4 8μl ij (P i P j )
58 Relative permeabilities Flux passing through the pore throat ij based on Hagen-Poiseuille s law q ij = πr ij 4 8μl ij (P i P j ) Based on volumetric conservation law for all pore bodies: Ni j=1 q ij =0,N i number of pore throats connected to pore body i.
59 Relative permeabilities Flux passing through the pore throat ij based on Hagen-Poiseuille s law q ij = πr ij 4 8μl ij (P i P j ) Based on volumetric conservation law for all pore bodies: Ni j=1 q ij =0,N i number of pore throats connected to pore body i. Linear system of equations will be resulted AP=B. A is a symmetric banded sparse matrix made of conductivities of pore throats.
60 k r -S w -a nw surfaces c Drainage and imbibition
61 Conclusions We can have unique surfaces for P c -S w -a nw as well as k r -S w -a nw for drainage and imbibition using our idealized pore-network model. Increasing aspect ratio decreases uniqueness of the surface. Further developments should be done; such as including angular cross sections to allow different mechanisms during imbibition and implementing for an unstructured anisotropic cross sections. Continuing for a real porous medium!
62 Case2: micromodel Joekar-Niasar, V., Hassanizadeh, S. M., Pyrak-Nolte, L. J., Berentsen, C., Simulating drainage and imbibition experiments in a high-porosity micro-model using an unstructured pore-network model, 2009, Water Resources Research, 45: W02430, doi: /2007wr
63 Specifications Topology: unstructured, anisotropic. Geometry: square cross-section Joekar-Niasar et al.2008,wrr
64 Challenge How can we use the concept of pore-network modelling for a high-porosity medium? Pore-body and pore throat concept does not hold!!!
65 Challenge How can we use the concept of pore-network modelling for a high-porosity medium? Pore-body and pore throat concept does not hold!!! Skeletonization using medial axis: defining pixels and bottle necks
66 Medial axis transform
67 Governing equation P e = σ nw (a + b)cosθ + (a + b) 2 cos 2 θ +4ab( π 4 θ 2cos( π 4 + θ)cosθ) 4( π 4 θ 2cos( π 4 + θ)cosθ)
68 Another challenge! How can we simulate imbibition process? High porosity and small aspect ration causes significant amount of cooperative filling. No snap-off occurs during imbibition.
69 Another challenge! How can we simulate imbibition process? High porosity and small aspect ration causes significant amount of cooperative filling. No snap-off occurs during imbibition. Defining a rule for no-trapping during imbibition
70 Cooperative filling
71 Cooperative filling
72 Simulations vs. experiments
73 Simulations vs. experiments
74 Conclusions Given geometry, pore-network model can be used to simulate not only P c -S w but also a nw -S w curves. P c -S w -a w surfaces resulted from simulation and experiments showed good agreement. This predictive tool can be useful for assessing the micromodel before manufacturing. It is very hard to simulate imbibition process in high porosity media.
75 Case3: glass-bead column Joekar-Niasar, V., Prodanovic, M., Wildenschild, D., Hassanizadeh, M., Translation of a Granular Porous Medium into a Pore-Network Model: Application for Capillary Pressure-Saturation and Interfacial Area-Saturation Relationships
76 Experiment Air-water system in glass beads. Image obtained through CT scanner.
77 Specifications Topology: unstructured, anisotropic Geometry: mixed hyperbolic polygonal cross sections for pore throats and prolate spheroid for pore bodies. Image analysis and data acquisition using 3DMA-Rock: skeleton, shape factor, pore unit volume, mean inscribed circle at pore throats andporebodiesaswellasmeanradiusofporebodies.
78 Pore throats geometry: irregular hyperbolic For a desired value of cross section area (or R) andκ, onecan calculate radii R 1, R 2, R 3,andanglesα, β, andγ.
79 Pore throats geometry: regular hyperbolic For a desired value of cross section area (or R) andκ, onecan calculate radius R and corner angle ϕ.
80 Entry pressure: irregular hyperbolic
81 Entry pressure: regular hyperbolic
82 Pore throats determination Determination of pore throats cross section based on shape factor distribution.
83 Numerical simulation Determination of cross sections types. Determining pore body geometries based on two min. and mean inscribed radii. Solving a nonlinear system to determine geometry of a cross section based on area and shape factor. Solving a nonlinear system to determine entry capillary pressure. Determining fluid occupancy during drainage and imbibition based on imposed capillary pressure.
84 Results:grain size distribution
85 Results: P c -S w curves
86 Results: a nw -S w curves
87 Conclusions Given geometry, pore-network model can be used to simulate not only the P c -S w but also a nw -S w curves. At least for granular porous medium, including shape factor to generate hyperbolic polygonal cross sections works well. There is enough consistency between grain size distribution and cross sections of pore throats.
88
Pore-Scale Analysis of Dynamics of Two-Phase Flow. in Porous Media
Pore-Scale Analysis of Dynamics of Two-Phase Flow Multi-Process Pore-scale Modelling in Porous Media Vahid J. Niasar Collaborators: S. Majid Hassanizadeh (Utrecht University) Nikos Karadimitriou (University
More informationNetwork Model Investigation of Interfacial Area, Capillary Pressure and Saturation Relationships in Granular Porous Media
WATER RESOURCES RESEARCH, VOL.???, XXXX, DOI:10.1029/, Network Model Investigation of Interfacial Area, Capillary Pressure and Saturation Relationships in Granular Porous Media V. Joekar-Niasar, 1 M. Prodanović
More informationAbility of Darcy s Law for Extension in Two- Phase Flow for Sedimentary Medium in Capillary Non-equilibrium Situations
Research Article imedpub Journals http://www.imedpub.com Resources, Recycling and Waste Management Ability of Darcy s Law for Extension in Two- Phase Flow for Sedimentary Medium in Capillary Non-equilibrium
More informationPaper E A pore network model for calculation of interfacial velocities
Paper E A pore network model for calculation of interfacial velocities Submitted to Advances in Water Resources, fall 21. A PORE NETWORK MODEL FOR CALCULATION OF INTERFACIAL VELOCITIES H.F. Nordhaug a,
More informationPore-scale modeling extension of constitutive relationships in the range of residual saturations
WATER RESOURCES RESEARCH, VOL. 37, NO. 1, PAGES 165 170, JANUARY 2001 Pore-scale modeling extension of constitutive relationships in the range of residual saturations Rudolf J. Held and Michael A. Celia
More informationCoupled free-flow and porous media flow: a numerical and experimental investigation
Coupled free-flow and porous media flow: a numerical and experimental investigation Master s Thesis Pavan Cornelissen 3863514 Supervisors: Kilian Weishaupt, MSc prof. dr. ir. Rainer Helmig prof. dr. ir.
More information2. Modeling of shrinkage during first drying period
2. Modeling of shrinkage during first drying period In this chapter we propose and develop a mathematical model of to describe nonuniform shrinkage of porous medium during drying starting with several
More informationNumerical Treatment of Two-Phase Flow in Porous Media Including Specific Interfacial Area
Procedia Computer cience Volume 51, 2015, Pages 1249 1258 ICC 2015 International Conference On Computational cience Numerical Treatment of To-Phase Flo in Porous Media Including pecific Interfacial Area
More informationA multiscale method coupling network and continuum models in porous media II singleand two-phase flows
A multiscale method coupling network and continuum models in porous media II singleand two-phase flows Jay Chu, Björn Engquist, Maša Prodanović and Richard Tsai 1 Introduction Modeling and computing transport
More informationTHE EFFECTS OF DISPLACEMENT RATE AND WETTABILITY ON IMBIBITION RELATIVE PERMEABILITIES
SCA2005-39 1/12 THE EFFECTS OF DISPLACEMENT RATE AND WETTABILITY ON IMBIBITION RELATIVE PERMEABILITIES V.H. Nguyen 1, A.P. Sheppard 2, M.A. Knackstedt 2 and W.V. Pinczewski 1 1 School of Petroleum Engineering,
More informationThis paper was prepared for presentation at the International Symposium of the Society of Core Analysts held in Toronto, Canada, August 2005
SCA2005-18 1/12 INFLUENCE OF WETTABILITY ON RESISTIVITY OF SEDIMENTARY ROCKS Mikhail Gladkikh 1, The University of Texas at Austin, Institute of Computational and Engineering Sciences Steven L. Bryant,
More informationEvaporation rates from square capillaries limited by corner flow
1 2 3 Evaporation rates from square capillaries limited by corner flow viscous losses Frouke Hoogland, May 2012 4 5 A master thesis for the master program Environmental Hydrogeology at the Department of
More informationPore-Level Bénard Marangoni Convection in Microgravity
Pore-Level Bénard Marangoni Convection in Microgravity Peyman Mohammadmoradi, and Apostolos Kantzas * Chemical and Petroleum Engineering Department, University of Calgary *Corresponding author: 2500 University
More informationComplexity of Two-Phase Flow in Porous Media
1 Complexity of Two-Phase Flow in Porous Media Rennes September 16, 2009 Eyvind Aker Morten Grøva Henning Arendt Knudsen Thomas Ramstad Bo-Sture Skagerstam Glenn Tørå Alex Hansen 2 Declining oil production,
More informationExperimental measurement of parameters governing flow rates and partial saturation in paper-based microfluidic devices
Experimental measurement of parameters governing flow rates and partial saturation in paper-based microfluidic devices Dharitri Rath 1, Sathishkumar N 1, Bhushan J. Toley 1* 1 Department of Chemical Engineering
More informationPORE-SCALE PHASE FIELD MODEL OF TWO-PHASE FLOW IN POROUS MEDIUM
Excerpt from the Proceedings of the COMSOL Conference 2010 Paris PORE-SCALE PHASE FIELD MODEL OF TWO-PHASE FLOW IN POROUS MEDIUM Igor Bogdanov 1*, Sylvain Jardel 1, Anis Turki 1, Arjan Kamp 1 1 Open &
More informationCritical capillary number: Desaturation studied with fast X-ray computed microtomography
GEOPHYSICAL RESEARCH LETTERS, VOL. 41, 55 60, doi:10.1002/2013gl058075, 2014 Critical capillary number: Desaturation studied with fast X-ray computed microtomography Ryan T. Armstrong, 1 Apostolos Georgiadis,
More informationPHYSICS OF FLUID SPREADING ON ROUGH SURFACES
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 5, Supp, Pages 85 92 c 2008 Institute for Scientific Computing and Information PHYSICS OF FLUID SPREADING ON ROUGH SURFACES K. M. HAY AND
More informationAnisotropic permeabilities evolution of reservoir rocks under pressure:
Extended reserves Clean refining Fuel-efficient vehicles Diversified fuels Controlled CO 2 Anisotropic permeabilities evolution : New experimental and numerical approaches (1) Dautriat J. 1-2*, Gland N.
More informationNUMERICAL INVESTIGATION OF THE DEPENDENCE OF RESIDUAL OIL SATURATION ON GEOMETRY, WETTABILITY, INITIAL OIL SATURATION AND POROSITY
SCA2014-042 1/6 NUMERICAL INVESTIGATION OF THE DEPENDENCE OF RESIDUAL OIL SATURATION ON GEOMETRY, WETTABILITY, INITIAL OIL SATURATION AND POROSITY Xiaobo Nie, Y. Mu, R. Sungkorn, V. Gundepalli and J. Toelke
More informationX-Ray Microtomography and NMR: complimentary tools for evaluation of pore structure within a core
SCA2014-038 1/6 X-Ray Microtomography and NMR: complimentary tools for evaluation of pore structure within a core Aleksandr Denisenko, Ivan Yakimchuk Schlumberger This paper was prepared for presentation
More informationDynamic pore-scale network model (PNM) of water imbibition in porous media Li, Juan; McDougall, Steven Robert; Sorbie, Kenneth Stuart
Heriot-Watt University Heriot-Watt University Research Gateay Dynamic pore-scale netork model (PNM) of ater imbibition in porous media Li, Juan; McDougall, Steven Robert; Sorbie, Kenneth Stuart Published
More informationFluid Flow Fluid Flow and Permeability
and Permeability 215 Viscosity describes the shear stresses that develop in a flowing fluid. V z Stationary Fluid Velocity Profile x Shear stress in the fluid is proportional to the fluid velocity gradient.
More informationarxiv: v1 [cond-mat.soft] 5 Jan 2016
Pore-scale simulations of drainage in granular materials: finite size effects and the representative elementary volume arxiv:1601.00840v1 [cond-mat.soft] 5 Jan 2016 Abstract Chao Yuan, Bruno Chareyre,
More informationCORRELATION OF THE MULTIPHASE FLOW COEFFICIENTS OF POROUS MEDIA WITH WETTABILITY: A PORE NETWORK APPROACH
CMWRXVI CORRELATION OF THE MULTIPHASE FLOW COEFFICIENTS OF POROUS MEDIA WITH WETTABILITY: A PORE NETWORK APPROACH CHRISTOS D. TSAKIROGLOU Institute of Chemical Engineering and High Temperature Chemical
More informationPore-scale modelling of Ostwald ripening
Downloaded from https://www.cambridge.org/core. Stanford Libraries, on 13 Feb 218 at 16:51:18, subject to the Cambridge Core terms of use, available at https://www.cambridge.org/core/terms. https://doi.org/1.117/jfm.217.72
More informationPORE NETWORK MODELS TO CALCULATE TRANSPORT AND ELECTRICAL PROPERTIES OF SINGLE OR DUAL-POROSITY ROCKS
SCA2005-35 /6 PORE NETWORK MODELS TO CALCULATE TRANSPORT AND ELECTRICAL PROPERTIES OF SINGLE OR DUAL-POROSITY ROCKS S. Békri, C.Laroche, O. Vizika Institut Français du Pétrole, et 4 avenue de Bois-Préau,
More informationCOMPARISON OF WETTABILITY AND CAPILLARY EFFECT EVALUATED BY DIFFERENT CHARACTERIZING METHODS
18 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS COMPARISON OF WETTABILITY AND CAPILLARY EFFECT EVALUATED BY DIFFERENT CHARACTERIZING METHODS S.K. Wang*, M. Li*, Y.Z. Gu, Y.X. Li and Z.G. Zhang Key
More informationin Mixed-Wet Systems Andrey Ryazanov Heriot-Watt University February 2012
Pore-Scale Network Modelling of Residual Oil Saturation in Mixed-Wet Systems Andrey Ryazanov Submitted for the degree of Doctor of Philosophy Heriot-Watt University Institute of Petroleum Engineering February
More informationAnalysis of Pore Structure Effects on Diffusive Reactive Transport in Bentonite via Pore Network Models 14068
Analysis of Pore Structure Effects on Diffusive Reactive Transport in Bentonite via Pore Network Models 4068 Qingrong Xiong, Simon Kwong 2, Andrey P Jivkov Research Centre for Radwaste & Decommissioning
More informationCapillary pressure as a unique function of electric permittivity and water saturation
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L13306, doi:10.1029/2007gl029674, 2007 Capillary pressure as a unique function of electric permittivity and water saturation Willem-Jan
More informationPUBLICATIONS. Journal of Geophysical Research: Solid Earth
PUBLICATIONS RESEARCH ARTICLE Key Points: Simulations of drainage/imbibition in water-wet pipe networks Approximately universal of the sense percolation theory in network simulations Resistivity index
More information2nd International Conference Mechanics of Unsaturated Soils 7 th 9 th March 2007
2nd International Conference Determination of the Soil Water Retention Curve and the Unsaturated Hydraulic Conductivity from the Particle Size Distribution Alexander Scheuermann & Andreas Bieberstein Motivation
More informationModeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures
Modeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures ENRU LIU, SERAFEIM VLASTOS, and XIANG-YANG LI, Edinburgh
More informationPore-Scale Modeling of Ostwald Ripening
Under consideration for publication in J. Fluid Mech. 1 Pore-Scale Modeling of Ostwald Ripening Jacques A. de Chalendar 1, Charlotte Garing 1 and Sally M. Benson 1 1 Department of Energy Resources Engineering,
More information1 Modeling Immiscible Fluid Flow in Porous Media
Excerpts from the Habilitation Thesis of Peter Bastian. For references as well as the full text, see http://cox.iwr.uni-heidelberg.de/people/peter/pdf/bastian_habilitationthesis.pdf. Used with permission.
More informationPreferential Flow and Extensions to Richards Equation
Preferential Flow and Extensions to Richards Equation David A. DiCarlo Dept. of Petroleum and Geosystems Engineering University of Texas at Austin Preferential flow Uniform infiltration Preferential flow
More informationWATER RESOURCES RESEARCH, VOL. 48, W10519, doi: /2012wr011862, 2012
WATER RESOURCES RESEARCH, VOL. 48, W10519, doi:10.1029/2012wr011862, 2012 A pore network model study of the fluid-fluid interfacial areas measured by dynamic-interface tracer depletion and miscible displacement
More informationThe role of capillary pressure curves in reservoir simulation studies.
The role of capillary pressure curves in reservoir simulation studies. M. salarieh, A. Doroudi, G.A. Sobhi and G.R. Bashiri Research Inistitute of petroleum Industry. Key words: Capillary pressure curve,
More informationMeasurement and prediction of the relationship between capillary pressure, saturation, and interfacial area in a NAPL water glass bead system
WATER RESOURCES RESEARCH, VOL. 46,, doi:10.1029/2009wr007786, 2010 Measurement and prediction of the relationship between capillary pressure, saturation, and interfacial area in a NAPL water glass bead
More informationDRAINAGE AND IMBIBITION CAPILLARY PRESSURE CURVES OF CARBONATE RESERVOIR ROCKS BY DIGITAL ROCK PHYSICS
SCA2012-56 1/6 DRAINAGE AND IMBIBITION CAPILLARY PRESSURE CURVES OF CARBONATE RESERVOIR ROCKS BY DIGITAL ROCK PHYSICS Y. Mu 1, Q. Fang 1, C. Baldwin 1, J. Töelke 1, A. Grader 1, M. Dernaika 2, and Z. Kalam
More informationChapter Seven. For ideal gases, the ideal gas law provides a precise relationship between density and pressure:
Chapter Seven Horizontal, steady-state flow of an ideal gas This case is presented for compressible gases, and their properties, especially density, vary appreciably with pressure. The conditions of the
More informationUNIVERSITY OF CALGARY. Two-Phase Flow at the Pore-Scale Using the Volume of Fluid Method. Carla Jordana Sena Santiago A THESIS
UNIVERSITY OF CALGARY Two-Phase Flow at the Pore-Scale Using the Volume of Fluid Method by Carla Jordana Sena Santiago A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES IN PARTIAL FULFILMENT OF THE
More informationPore Scale Analysis of Oil Shale/Sands Pyrolysis
Pore Scale Analysis of Oil Shale/Sands Pyrolysis C.L. Lin, J.D. Miller, and C.H. Hsieh Department of Metallurgical Engineering College of Mines and Earth Sciences University of Utah Outlines Introduction
More informationPUBLICATIONS. Water Resources Research. Challenges in modeling unstable two-phase flow experiments in porous micromodels
PUBLICATIONS RESEARCH ARTICLE Key Points: We validate our numerical approach by comparison with drainage experiments We demonstrate the impact of three main sources of error We compare the results in a
More informationDissolution and precipitation during flow in porous media
1/25 Class project for GEOS 692: Transport processes and physical properties of rocks Dissolution and precipitation during flow in porous media Gry Andrup-Henriksen Fall 2006 1 2/25 Outline Introduction
More informationROLE OF PORE-SCALE HETEROGENEITY ON REACTIVE FLOWS IN POROUS MATERIALS: VALIDITY OF THE CONTINUUM REPRESENTATION OF REACTIVE TRANSPORT
ROLE OF PORE-SCALE HETEROGENEITY ON REACTIVE FLOWS IN POROUS MATERIALS: VALIDITY OF THE CONTINUUM REPRESENTATION OF REACTIVE TRANSPORT PETER C. LICHTNER 1, QINJUN KANG 1 1 Los Alamos National Laboratory,
More informationTANG Xiaoyan [a],* INTRODUCTION
Advances in Petroleum Exploration and Development Vol. 8, No. 2, 2014, pp. 6-12 DOI:10.3968/5958 ISSN 1925-542X [Print] ISSN 1925-5438 [Online] www.cscanada.net www.cscanada.org A Simulation Experimental
More informationHysteresis of Capillary Pressure, Resistivity Index and Relative Permeability in Different Carbonate Rock Types 1
PETROPHYSICS, VOL. 53, NO. 5 (OCTOBER 2012); PAGE 315 331; 28 FIGURES; 2 TABLES Hysteresis of Capillary Pressure, Resistivity Index and Relative Permeability in Different Carbonate Rock Types 1 Moustafa
More informationEstimating Permeability from Acoustic Velocity and Formation Resistivity Factor
5th Conference & Exposition on Petroleum Geophysics, Hyderabad-2004, India PP 582-587 and Formation Resistivity Factor Majid Nabi-Bidhendi Institute of Geophysics, University of Tehran, P.O. Box 14155-6466,
More informationPore-scale investigation of wettability effects on two-phase flow in porous media
Pore-scale investigation of wettability effects on two-phase flow in porous media A thesis submitted to the University of Manchester for the degree of Doctor of philosophy in the Faculty of Science and
More informationI. Borsi. EMS SCHOOL ON INDUSTRIAL MATHEMATICS Bedlewo, October 11 18, 2010
: an : an (Joint work with A. Fasano) Dipartimento di Matematica U. Dini, Università di Firenze (Italy) borsi@math.unifi.it http://web.math.unifi.it/users/borsi porous EMS SCHOOL ON INDUSTRIAL MATHEMATICS
More informationPore Scale Modeling of Capillary Pressure Curves in 2D Rock Images
Pore Scale Modeling of Capillary Pressure Curves in 2D Rock Images by Yingfang Zhou Thesis submitted in partial fulfillment of the requirements for the degree of Philosophiae Doctor (PhD) Faculty of Science
More informationNetwork Modelling of the Influence of Swelling on the Transport Behaviour of Bentonite
geosciences Article Network Modelling of the Influence of Swelling on the Transport Behaviour of Bentonite Ignatios Athanasiadis, Simon Wheeler and Peter Grassl * School of Engineering, University of Glasgow,
More informationCOMPARING DIFFERENT METHODS FOR CAPILLARY PRESSURE MEASUREMENTS
COMPARING DIFFERENT METHODS FOR CAPILLARY PRESSURE MEASUREMENTS M. Sarwaruddin ), OleTorsæter ), and Arne Skauge 2) ) Norwegian University of Science &Technology 2) Norsk Hydro Abstract Capillary pressure
More informationRELATING CAPILLARY PRESSURE TO INTERFACIAL AREAS
Proceedings of the Project Review, Geo-Mathematical Imaging Group (Purdue University, West Lafayette IN), Vol. 1 (2008) pp. 221-238. RELATING CAPILLARY PRESSURE TO INTERFACIAL AREAS LAURA J. PYRAK-NOLTE,
More informationPore scale consideration in unstable gravity driven finger flow
Pore scale consideration in unstable gravity driven finger flow Steenhuis, T. S., C. E. Baver, B. Hasanpour, C. R. Stoof, D. A. DiCarlo, and J. S. Selker (2013), Pore scale consideration in unstable gravity
More informationOn the origin of Darcy s law 1
Chapter 1 On the origin of Darcy s law 1 Cyprien Soulaine csoulain@stanford.edu When one thinks about porous media, the immediate concepts that come to mind are porosity, permeability and Darcy s law.
More informationSimultaneous Measurement of Capillary Pressure and Dielectric Constant in Porous Media
PIERS ONLINE, VOL. 3, NO. 4, 2007 549 Simultaneous Measurement of Capillary Pressure and Dielectric Constant in Porous Media W. J. Plug, L. M. Moreno, J. Bruining, and E. C. Slob Delft University of Technology,
More informationINFERRING RELATIVE PERMEABILITY FROM RESISTIVITY WELL LOGGING
PROCEEDINGS, Thirtieth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 3-February 2, 25 SGP-TR-76 INFERRING RELATIVE PERMEABILITY FROM RESISTIVITY WELL LOGGING
More informationUnsaturated Flow (brief lecture)
Physical Hydrogeology Unsaturated Flow (brief lecture) Why study the unsaturated zone? Evapotranspiration Infiltration Toxic Waste Leak Irrigation UNSATURATAED ZONE Aquifer Important to: Agriculture (most
More informationA two-dimensional network simulator for two-phase flow in porous media
A two-dimensional network simulator for two-phase flow in porous media Eyvind Aker, Knut Jørgen Måløy Department of Physics University of Oslo, N-0316 Oslo, Norway Alex Hansen Department of Physics, Norwegian
More informationPore-Scale Analysis of Interfacial Instabilities and Impact of Heterogeneity on Relative Permeability by Lattice Boltzmann Method
Louisiana State University LSU Digital Commons LSU Master's Theses Graduate School 3-31-2018 Pore-Scale Analysis of Interfacial Instabilities and Impact of Heterogeneity on Relative Permeability by Lattice
More informationThe Pennsylvania State University. The Graduate School. College of Earth and Mineral Sciences NUMERICAL STUDY OF FRACTURE APERTURE CHARACTERISTICS AND
The Pennsylvania State University The Graduate School College of Earth and Mineral Sciences NUMERICAL STUDY OF FRACTURE APERTURE CHARACTERISTICS AND THEIR IMPACT ON SINGLE-PHASE FLOW AND CAPILLARY-DOMINATED
More informationUNDERSTANDING IMBIBITION DATA IN COMPLEX CARBONATE ROCK TYPES
SCA2014-059 1/6 UNDERSTANDING IMBIBITION DATA IN COMPLEX CARBONATE ROCK TYPES Moustafa Dernaika 1, Zubair Kalam 2, Svein Skjaeveland 3 1 Ingrain Inc.-Abu Dhabi, 2 ADCO, 3 University of Stavanger This paper
More informationSCA : A STRUCTURAL MODEL TO PREDICT TRANSPORT PROPERTIES OF GRANULAR POROUS MEDIA Guy Chauveteau, IFP, Yuchun Kuang IFP and Marc Fleury, IFP
SCA2003-53: A STRUCTURAL MODEL TO PREDICT TRANSPORT PROPERTIES OF GRANULAR POROUS MEDIA Guy Chauveteau, IFP, Yuchun Kuang IFP and Marc Fleury, IFP This paper was prepared for presentation at the International
More informationNumerical Simulation of Diagenetic Alteration and Its Effect on Residual Gas in Tight Gas Sandstones
Transp Porous Med (2013) 96:39 62 DOI 10.1007/s11242-012-0072-3 Numerical Simulation of Diagenetic Alteration and Its Effect on Residual Gas in Tight Gas Sandstones Maša Prodanović Steven L. Bryant J.
More informationSIMULATING IN-SITU CONDITIONS FOR DIGITAL CORE ANALYSIS
SCA2015-032 1/11 SIMULATING IN-SITU CONDITIONS FOR DIGITAL CORE ANALYSIS S. Linden 1,2, T. Cvjetkovic 2, E. Glatt 2, J.-O. Schwarz 2 and A. Wiegmann 2 1 Fraunhofer ITWM, Fraunhofer-Platz 1, 67665 Kaiserslautern,
More informationRelative permeabilities and coupling effects in steady-state gas-liquid flow in porous media: A lattice Boltzmann study
PHYSICS OF FLUIDS 21, 092104 2009 Relative permeabilities and coupling effects in steady-state gas-liquid flow in porous media: A lattice Boltzmann study Haibo Huang a and Xi-yun Lu Department of Modern
More informationThe Physical Origin of Interfacial Coupling in Two-Phase Flow through Porous Media
Transport in Porous Media 44: 109 122, 2001. c 2001 Kluwer Academic Publishers. Printed in the Netherls. 109 The Physical Origin of Interfacial Coupling in Two-Phase Flow through Porous Media RAMON G.
More informationHysteresis and Interfacial Energies in Smooth-walled. Microfluidic Channels
WWR orig. June 12, 2009 Hysteresis and Interfacial Energies in Smooth-walled Microfluidic Channels Yihong Liu 1, D. D. Nolte 1 and L. J. Pyrak-Nolte 1,2 1. Dept. of Physics, Purdue University 2. Dept.
More informationRelative Permeability Measurement and Numerical Modeling of Two-Phase Flow Through Variable Aperture Fracture in Granite Under Confining Pressure
GRC Transactions, Vol. 36, 2012 Relative Permeability Measurement and Numerical Modeling of Two-Phase Flow Through Variable Aperture Fracture in Granite Under Confining Pressure Noriaki Watanabe, Keisuke
More informationPermeability and fluid transport
Permeability and fluid transport Thermal transport: Fluid transport: q = " k # $p with specific discharge (filter velocity) q [m s 1 ] pressure gradient p [N m 3 ] dynamic viscosity η [N s m 2 ] (intrinsic)
More informationMEASUREMENT OF CAPILLARY PRESSURE BY DIRECT VISUALIZATION OF A CENTRIFUGE EXPERIMENT
MEASUREMENT OF CAPILLARY PRESSURE BY DIRECT VISUALIZATION OF A CENTRIFUGE EXPERIMENT Osamah A. Al-Omair and Richard L. Christiansen Petroleum Engineering Department, Colorado School of Mines ABSTRACT A
More informationMicro-Structural Rock Modeling: Methodology and Application in Formation Evaluation
Micro-Structural Rock Modeling: Methodology and Application in Formation Evaluation Guodong Jin 1, Chun Lan 1, Wei Shao 1, and Songhua Chen 1, 2 1 Baker Hughes Incorporated, Houston, Texas, USA 2 Now at
More informationDetermination of Capillary pressure & relative permeability curves
Determination of Capillary pressure & relative permeability curves With Refrigerated Centrifuge Multi speed centrifuge experiments Introduction As the porous diaphragm method see Vinci s CAPRI, the centrifuge
More informationCorrelation Between Resistivity Index, Capillary Pressure and Relative Permeability
Proceedings World Geothermal Congress 2010 Bali, Indonesia, 25-29 April 2010 Correlation Between Resistivity Index, Capillary Pressure and Kewen Li Stanford Geothermal Program, Stanford University, Stanford,
More informationTHE EFFECT OF WATER SATURATION ON GAS SLIP FACTOR BY PORE SCALE NETWORK MODELING
SCA00-53 1/6 THE EFFECT OF WATER SATURATION ON GAS SLIP FACTOR BY PORE SCALE NETWORK MODELING Liu Qingjie, Liu Baohua, Li Xianbing, Yan Shouguo Research Institute of Petroleum Exploration and Development,
More informationPORE-SCALE SIMULATION OF MASS TRANSFER ACROSS IMMISCIBLE INTERFACES
PORE-SCALE SIMULATION OF MASS TRANSFER ACROSS IMMISCIBLE INTERFACES A report submitted to the department of Energy Resources Engineering of Stanford University in partial fulfillment of the requirements
More informationMICRO-SCALE EXPERIMENTAL INVESTIGATION OF THE EFFECT OF HYDROSTATIC STRESS ON PORE- SPACE DEFORMATION AND FLUID OCCUPANCY
SCA2016-026 1/12 MICRO-SCALE EXPERIMENTAL INVESTIGATION OF THE EFFECT OF HYDROSTATIC STRESS ON PORE- SPACE DEFORMATION AND FLUID OCCUPANCY M. Asadollahkhan Vali 1, A.H. Alizadeh 1, M. Piri 1, and J. Wallace
More information1. Water in Soils: Infiltration and Redistribution
Contents 1 Water in Soils: Infiltration and Redistribution 1 1a Material Properties of Soil..................... 2 1b Soil Water Flow........................... 4 i Incorporating K - θ and ψ - θ Relations
More informationDynamic coupling of pore-scale and reservoir-scale models for multiphase flow
WATER RESOURCES RESEARCH, VOL. 49, 5973 5988, doi:10.1002/wrcr.20430, 2013 Dynamic coupling of pore-scale and reservoir-scale models for multiphase flow Qiang Sheng 1 and Karsten Thompson 2 Received 30
More informationPore radius distribution and fractal dimension derived from spectral induced polarization
Pore radius distribution and fractal dimension derived from spectral induced polarization Zeyu Zhang 1 & Andreas Weller 2 1) Southwest Petroleum University, China & Bundesanstalt für Materialforschung
More informationAn Immersed Boundary Method for Computing Anisotropic Permeability of Structured Porous Media
An Immersed Boundary Method for Computing Anisotropic Permeability of Structured Porous Media D.J. Lopez Penha a, B.J. Geurts a, S. Stolz a,b, M. Nordlund b a Dept. of Applied Mathematics, University of
More informationExperimental investigation of dynamic effects in capillary pressure: Grain size dependency and upscaling
WATER RESOURCES RESEARCH, VOL. 46,, doi:10.1029/2009wr008881, 2010 Experimental investigation of dynamic effects in capillary pressure: Grain size dependency and upscaling Geremy Camps Roach, 1 Denis M.
More informationTwo-phase Flow Calculations in Pore Unit-cells Implementing Mixed FEM/Lattice-Boltzmann Simulators
Two-phase Flow Calculations in Pore Unit-cells Implementing Mixed FEM/Lattice-Boltzmann Simulators E. D. Skouras,*, A. N. Kalarakis, M. S. Valavanides 3, and V. N. Burganos Foundation for Research and
More informationUnsteady State Relative Permeability and Capillary Pressure Estimation of Porous Media
CMWRXVI Unsteady tate Relative Permeability and Capillary Pressure Estimation of Porous Media M. H. Ghazanfari 1,2, M. Khodabakhsh 1,2, R. Kharrat 2, D. Rashtchian 1,. Vossoughi 3 1 Chemical and Petroleum
More informationdynamics of f luids in porous media
dynamics of f luids in porous media Jacob Bear Department of Civil Engineering Technion Israel Institute of Technology, Haifa DOVER PUBLICATIONS, INC. New York Contents Preface xvii CHAPTER 1 Introduction
More informationStudy of Multi-phase Flow in Porous Media: Comparison of SPH Simulations with Micro-model Experiments
Transp Porous Med DOI 10.1007/s11242-015-0599-1 Study of Multi-phase Flow in Porous Media: Comparison of SPH Simulations with Micro-model Experiments P. Kunz 1 I. M. Zarikos 2 N. K. Karadimitriou 3 M.
More informationPermeability Estimates & Saturation Height Functions: A talk of two halves. Dr Joanne Tudge LPS Petrophysics 101 Seminar 17 th March 2016
Permeability Estimates & Saturation Height Functions: A talk of two halves Dr Joanne Tudge LPS Petrophysics 101 Seminar 17 th March 2016 Permeability: What is it? How do we measure it? Why do we need it?
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 informationX-RAY MICRO-TOMOGRAPHY OF PORE-SCALE FLOW AND TRANSPORT. University of California Davis. Dorthe Wildenschild & Annette Mortensen
X-RAY MICRO-TOMOGRAPHY OF PORE-SCALE FLOW AND TRANSPORT Jan W. Hopmans University of California Davis Volker Clausnitzer Dorthe Wildenschild & Annette Mortensen ISSUES: Measurements and modeling of water
More informationNumerical Simulation of Single-Phase and Multiphase Non-Darcy Flow in Porous and Fractured Reservoirs
Transport in Porous Media 49: 209 240, 2002. 2002 Kluwer Academic Publishers. Printed in the Netherlands. 209 Numerical Simulation of Single-Phase and Multiphase Non-Darcy Flow in Porous and Fractured
More informationSurface permeability, capillary transport and the Laplace-Beltrami problem.
Surface permeability, capillary transport and the Laplace-Beltrami problem. Penpark Sirimark, Alex V. Lukyanov & Tristan Pryer School of Mathematical and Physical Sciences, University of Reading, Reading
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 informationMicrofluidics 1 Basics, Laminar flow, shear and flow profiles
MT-0.6081 Microfluidics and BioMEMS Microfluidics 1 Basics, Laminar flow, shear and flow profiles 11.1.2017 Ville Jokinen Outline of the next 3 weeks: Today: Microfluidics 1: Laminar flow, flow profiles,
More informationMulti-rate mass transfer modeling of two-phase flow in highly heterogeneous fractured and porous media
Multi-rate mass transfer modeling of two-phase flow in highly heterogeneous fractured and porous media Jan Tecklenburg a,, Insa Neuweiler a, Jesus Carrera b, Marco Dentz b a Institute of Fluid Mechanics
More informationExperiment and Simulation on NMR and Electrical Measurements on Liège Chalk
The Open-Access Journal for the Basic Principles of Diffusion Theory, Experiment and Application Experiment and Simulation on NMR and Electrical Measurements on Liège Chalk Liangmou Li, Igor Shikhov, Yong
More informationpressure-saturations relationships for two-phase flow in 3D porous media: implications micro-heterogeneities
Loughborough University Institutional Repository Dynamic effects in capillary pressure-saturations relationships for two-phase flow in 3D porous media: implications of micro-heterogeneities This item was
More informationMacro-Scale Dynamic Effects in Homogeneous and Heterogeneous Porous Media
Transp Porous Med (2005) 58:121 145 Springer 2005 DOI 10.1007/s11242-004-5472-6 Macro-Scale Dynamic Effects in Homogeneous and Heterogeneous Porous Media SABINE MANTHEY 1,, S. MAJID HASSANIZADEH 2 and
More information