Upscaling of two-phase flow processes in heterogeneous porous media: determination of constitutive relationships
|
|
- Marianna Blair
- 6 years ago
- Views:
Transcription
1 28 Calibration and Reliability in Groundwater Modelling: A Few Steps Closer to Reality (Proceedings of ModclCARH'2002 (Prague. Czech Republic. June 2002). IAHS Publ. no Upscaling of two-phase flow processes in heterogeneous porous media: determination of constitutive relationships RAINER HELMIG University of Stuttgart, Institute of Hydraulic Engineering, Pfaffenwaldring 61, D Stuttgart, Germany helmig(5),i ws.uni-stuftgart.de CHRISTOPHERUS BRAUN Technologieberatung Grundwasser und Umwelt GmbH (TGU), Maria Trost 3, D Koblenz, Germany SABINE MAM IIFY University of Stuttgart, Institute of Hydraulic Engineering, Pfaffenwaldring 61, D Stuttgart, Germany Abstract For the numerical simulation of groundwater contamination by e.g. DNAPL, effective parameters representing the macroscale are needed to ensure reasonable computation times. When we use effective parameters on the macroscale, they should reproduce the influence of local-scale heterogeneities on the multiphase flow processes. Our upscaling approach attempts to incorporate the influence of local-scale heterogeneities on the multiphase flow behaviour at the macroscale in the constitutive relationships, taking into account that, in addition to gravity, capillary forces, which are most noticeable at the interface of heterogeneities, dominate the flow processes. The resultant relative permeability-saturation relationship for the macroscale shows a saturation dependent anisotropy. Moreover residual saturations of the nonwetting phase and hysteresis effects are introduced at the macroscale although they do not have a counterpart at the local scale. Keywords effective parameters; hysteresis; relative permeability-saturation relationship; saturation dependent anisotropy; upscaling INTRODUCTION When subsurface contamination by dense non-aqueous phase liquids (DNAPL) is to be assessed, e.g. with regard to the application of remediation techniques, numerical simulations are often employed. Ideally, effective parameters for the macroscale are required which make a discrete reproduction of local-scale heterogeneities futile. Local-scale heterogeneities can have an important impact on multiphase flow processes in the subsurface. Zones of lower permeability and higher entry pressure than the surrounding material cause a pooling of DNAPL until the critical saturation and thus the entry pressure of the lower permeability region is overcome. The pooling causes a lateral spreading of DNAPL. The residual saturation of the DNAPL is also increased. These effects should be captured by the effective parameters for the macroscale. We achieve a representation of these processes in the constitutive
2 Upscaling of two-phase flow processes in heterogeneous porous media 29 relationships by means of an upscaling procedure from the local to the macroscale which takes into account the structures of the local-scale heterogeneities. At the local-scale, hysteresis and residual saturations are not considered. Possible dynamic effects or functional dependencies of the constitutive relationships, e.g. on the interfacial area, are also neglected. Furthermore, the porous medium is assumed to be isotropic on the local scale. Thus, the intrinsic permeability and the relative permeabilitysaturation relationship are scalar functions. The constitutive relationships on the local scale are parameterized after the well-known approach of Brooks & Corey/Burdine. Description of the upscaling approach In our upscaling approach, one crucial assumption is an equilibrium of (capillary) forces, achieved by reducing the number of variables. The application of one capillary pressure at the boundaries gives a unique saturation distribution. To compute the saturation distribution for a heterogeneous system, made up of different REV, a sitepercolation model (e.g. Yortsos et al, 1993; Kueper & Girgrah, 1994) can be applied, for example. The weighted arithmetic mean of the saturation distribution yields one value of the averaged, macroscopic capillary pressure-saturation relationship. The procedure is repeated with different capillary pressures until the macroscopic P C S U. relationship is obtained with sufficient accuracy. With the saturation distribution, it is also possible to compute a distribution of the relative permeabilities and thus the conductivities. On the macroscale, the conductivities are averaged either arithmetically or harmonically, depending on the structure and the flow direction, with a renormalization method (Williams, 1992; King, 1996). With the renormalization, the size of the REV of the considered medium is also identified. The conductivities thus derived are then set into a relation with the saturated conductivity. As a result, we get a macroscopic relative permeability saturation relationship. One advantage of the approach is its easy applicability. Furthermore, its physical motivation facilitates a physical interpretation of the upscaled parameters. In the following, the approach will be applied to a geostatistical example and to a bench-scale laboratory experiment. Geostatistical example With the help of a geostatistical example we will show how hysteresis effects can arise at the macroscale although, on the local scale, hysteretic behaviour is neglected in the constitutive relationships. A permeability field (log(&o,rc/) = , standard deviation s = , anisotropy factor a = 10) is generated with rf2d (random field 2d; van Lent, 1992; Dykaar & Kitanidis, 1992). Then the Brooks-Corey parameter P c/, the entry pressure, is scaled (with the help of the Leverett function): with the reference values for Âro^/as the average permeability and P ( i.ref= 550 Pa (see Fig. 1). The Brooks-Corey parameter X is constant (k = 2.0) over the whole domain. (1)
3 30 Rainer Helmig et al. Pc [Pa] ^ X Fig. 1 Generated heterogeneity field (64 2 cells) (Braun, 2000) water saturation [-] Fig. 2 Macroscopic P C ~S. for the heterogeneity field (Braun, 2000). The macroscopic P C S W relationship for this heterogeneity field is computed as described above. Although on the local scale no hysteresis effects were assumed in the P C S W relationship for the different REV, process-dependent relationship evolve at the macroscale (see Fig. 2). For low capillary pressures, the nonwetting phase is trapped in the domain, leading to a macroscopic residual saturation of the nonwetting phase. The trapping occurs when high permeability regions with low entry pressure still show nonwetting phase saturations at low capillary pressures whereas surrounding regions with lower permeabilities and higher entry pressure are already fully saturated with water and thus do not provide a mobility for the DNAPL to flow out of the high permeability regions. The relative permeability-saturation relationship for the drainage and the imbibition cycle both show a saturation-dependent anisotropy (see Fig. 3 and 4).
4 Upscaling of two-phase flow processes in heterogeneous porous media S water saturation S [-] Fig. 3 k-s relationships for imbibition of the heterogeneity field (Braun, 2000) water saturation S [-] Fig. 4 kr-s relationships for the drainage cycle of the heterogeneity field (Braun, 2000). The relative permeabilities of the nonwetting phase do not differ from drainage to imbibition because hysteresis effects would arise at small capillary pressures when the saturation of the nonwetting phase is small. But the nonwetting phase only becomes mobile at saturations higher than S = This is a threshold value, below which there is no continuous path for the nonwetting phase in the domain. It corresponds to a residual saturation. The threshold value depends on various parameters, e.g. the dimension of the percolation model or the kind of spatial discretization (Stauffer, 1985). The relative permeabilities of the drainage and imbibition cycle of the wetting phase show a similar behaviour for saturations of up to S w ~ 0.9. Then, the macroscopic residual saturation of the nonwetting phase is reached, i.e. the saturation of the wetting phase cannot be increased further. At the residual saturation of the nonwetting phase, the relative permeability of the wetting phase becomes 1. Thus, a hysteresis effect on the macroscale arises not only for the P c -S w relationship but also for the k n -S w relationship. Verification and application example The upscaling approach was verified with a numerical experiment with a periodically layered medium where DNAPL infiltrated from the top (Braun et al, 1998). Here we
5 32 Rainer Helmig et al. Table 1 Properties of the sands for the VEGAS experiment (Allan, 1998). Fine sand Medium sand Coarse sand Permeability k 0 (m 2 ) 6.38 x lfr x lfr x lfr 10 Entry pressure P rf (Pa) Form factor X (-) Residual saturation S wr (-) Residual saturation S m. (-) gravel source zone constant DNAPL flux fine lenses coarse lenses 120 cm pressure tensiometers very fine sand Fig. 5 Set-up of infiltration experiment: distribution of the lenses (Allan, 1998). will demonstrate the application to an infiltration experiment (Fig. 5) conducted in the VEGAS facility, University of Stuttgart (Allan et al, 1998). Table 1 shows the properties of the sands used. During the experiment, DNAPL (trichlorethylene, TCE, density p = 1460 kg nr 3, dynamic viscosity u = 0.57 x 10" J Pa s) infiltrated for 2970 s with a constant flux of 29 ml min" 1 along 0.02 m of the top. Using this experiment, we compare two different upscaling approaches. First, we apply a simple geometric average of the parameters, taking only the volume fractions of the different sands into account. We then derive the averaged constitutive relationships for a geometric arrangement of the sands by means of our upscaling approach. For the simple geometric average, the resulting constitutive relationships resemble those of the medium sand because it has the highest volume fraction (see Fig. 6 and Fig. 7). The relative permeability-saturation relationships are isotropic. Figure 8 shows the results of the application of these constitutive relationships in a numerical simulation of the experiment (isolines) compared to the experimental results (hatched areas) after one hour. It is obvious that the averaging approach of the weighted geometric mean does not suffice to reproduce the experimental results correctly. The infiltration of the DNAPL in the vertical direction is overestimated whereas the dispersion in the horizontal direction is underestimated. Consequently, in addition to the information about the volume fractions, the geometric arrangement of the materials in the tank (lenses) is taken into account by applying
6 Upscaling of two-phase flow processes in heterogeneous porous media water saturation S [-] Fig. 6 Relative permeability-saturation relationship for the case of the geometric average (Braun, 2000) \ I rd Q. ' ' G water saturation S w [-] Fig. 7 Capillary pressure-saturation relationship for the case of the geometric average (Braun, 2000) m x TCE, observed: Fig. 8 Comparison of the saturation distribution of the DNAPL after one hour: numerical simulation with effective parameters derived by a weighted geometric mean (isolines) and experimental results (hatched areas) (Braun, 2000).
7 T ; 34 Rainer Helmig et al. U if coarse sand lortso [height 2 cm) J 1 1 ] fine sand loose (height 2 cm) crj - 0_ ' ' Q_ - fd s~ to (/) CD Q '5_ nj o I - _!,!,!,!,! i!,!,!, j, f water saturation S w [-] Fig. 9 Capillary pressure-saturation relationship for the case of the upscaling approach applied to a medium with lenses (Braun, 2000). water saturation S w [-] Fig. 10 Relative permeability-saturation relationship for the case of the upscaling approach applied to a medium with lenses (Braun, 2000). our upscaling approach (see Fig. 9). The relative permeability-saturation relation now shows an anisotropy for both phases (see Fig. 10). The anisotropy can be related to characteristic points in the capillary pressuresaturation relations of the different sands. As long as the capillary pressure is higher than the entry pressure of the coarse sand and lower than the entry pressure of the medium sand, the DNAPL can only infiltrate into the coarse sand. However, as the lenses are not continuous, this does not effect the mobility of the DNAPL. With decreasing water saturation, the capillary pressure increases until the entry pressure of the medium sand is exceeded. Only then does the DNAPL start to be mobile. Taken that the fluid is no longer mobile at a relative permeability lower than 1 x l(h, the mobilization of the DNAPL in the horizontal direction sets in at saturations of
8 Upscaling of two-phase flow processes in heterogeneous porous media 35 = and in the vertical direction at saturations of S, m, = These saturations can be seen as residual saturations of the nonwetting fluid. These residual saturations have to be exceeded before the nonwetting fluid becomes mobile and spreads in the respective direction. Moreover, the mobility of the DNAPL in the vertical direction is inhibited by the fine sand lens until its entry pressure is overcome as well. Thus, a larger DNAPL saturation has to be reached for the DNAPL to reach high mobilities in the vertical direction compared to the horizontal direction. The comparison of the numerical simulation with the effective parameters for the REV with the lenses to the experimental results (Fig. 11) shows a good agreement between the simulation and the experiment. The simulation reproduces the spreading in the horizontal direction. The isolines encompass the main part of the infiltrated area and do reproduce the characteristics of the mobilization process correctly m z X TCE, observed: W///W/W/M Fig. 11 Comparison of the saturation distribution of the DNAPL after one hour: numerical simulation with effective parameters derived with the upscaling approach for a medium with lenses (isolines) (Braun, 2000). CONCLUSION During the upscaling process for heterogeneous media, new effects may evolve on the macroscale which do not occur on the local scale. First, a saturation-dependent anisotropy of the relative permeability-saturation relationship can be observed. The twophase flow behaviour amplifies the anisotropy of the effective conductivity as compared to single-phase flow. Furthermore, direction-dependent macroscopic residual saturations evolve, at which the phases are immobile. Residual saturations of the nonwetting phase are an important parameter for assessing the success of remediation processes. Moreover, hysteresis effects can be observed on the macroscale, though on the local scale no process-dependent parameters were applied. The application of our upscaling procedure proves that the structures of a porous medium on the local scale, such as layers or lenses, have an important influence on the effective parameters on the macroscale. The incorporation of the geometry of these structures in the upscaling process enhances the quality of the effective parameters.
9 36 Rainer Helmig et al. REFERENCES Allan, J., Ewing, J., Braun, J. & Helmig, R. (1998) Scale effects in multiphase How modelling. In: Remediation of Chlorinated and Récalcitrant Compounds (ed. by G. B. Wickramanayake & R. E. Hinchee) (First Int. Conf., May 1998, Monterey, California), Battelle Press, Columbus, Ohio, USA. Braun, C. (2000) Ein Upscaling-Verfahren fur Mehrphascnstromungen in potosen Medien (An upscaling approach for multiphase flow in porous media). Scientific report 103, Dissertation, Institute of Hydraulic Engineering, University of Stuttgart, Stuttgart, Germany. Braun, C, Hassanizadeh, S. M. & Helmig, R. (1998) Two-pliase flow in stratified porous media effective parameters and constitutive relationships. In: Computational Methods in Water Resources XII (ed. by V. N. Burganos, G. P. Karatzas & A. C. Payatakes), Vol. I, Comp. Mech. Publ., Southampton, UK and Boston, USA. Dykaar, B. B. & Kitanidis, P. K. (1992) Determination of the effective hydraulic conductivity for heterogeneous porous media using a numerical spectral approach. 1 Method. Water Resour. Res. 28(4), King, P. R. (1996) Upscaling permeability: Error analysis for renormalisalion. Transp. Porous Media 23, Kuepcr, B. H. & Girgrah, B. (1994) Effective large-scale parameters for two-phase flow in heterogeneous porous media. In: Transport and Reactive Processes in Aquifers (ed. by T. H. Dracos & F. StaufFer) (Proc. IAHR Symp., April 1994, Ziirich, Switzerland), Balkema, Rotterdam, The Netherlands. Van Lent, T. (1992) Numerical spectral methods applied to flow in highly heterogeneous aquifers. PhD Thesis, Stanford University, California, USA. Stauffer, D. ( 1985) Introduction to Percolation Theory. Taylor & Francis, London, UK. Williams, J. K. (1989) Simple renormalisalion schemes for calculating effective properties of heterogeneous reservoirs. In: Mathematics of Oil Recovery (ed. by P. R. King) (First Int. Conf., July 1989, Cambridge, UK), Clarendon Press, Oxford, UK. Yortsos, Y. C, Satik, C., Bacri, J. C. & Salin, D. (1993) Large-scale percolation theorv of drainage. Transp. Porous Media 10,
Multiphase Flow in the Subsurface - Flow of a Light Nonaqueous Phase Liquid (LNAPL)
Multiphase Flow in the Subsurface - Flow of a Light Nonaqueous Phase Liquid (LNAPL) March 29, 2011 Wonyong Jang, Ph.D., P.E. Multimedia Environmental Simulations Laboratory (MESL) School of Civil and Environmental
More informationUpscaling sub-core scale heterogeneity
Upscaling sub-core scale heterogeneity Dave Cameron On behalf of the Benson Lab Stanford University Energy Resources Engineering 0 X-ray CT: Sub-core heterogeneity affects CO 2 transport/saturation Perrin,
More informationEXECUTIVE SUMMARY. 2. The effect of remediation on the distribution and mobility of both the LNAPL and water within the zone of interest.
EXECUTIVE SUMMARY For many decades, the oil production industry has recognized that significant limitations exist to complete extraction of oil from geologic formations. Attempts to recover fuels and crude
More informationEffects of heterogeneities on capillary pressure saturation relative permeability relationships
Journal of Contaminant Hydrology 56 (2002) 175 192 www.elsevier.com/locate/jconhyd Effects of heterogeneities on capillary pressure saturation relative permeability relationships Behzad Ataie-Ashtiani
More informationObjective. Technical Approach
Computational and Experimental Investigation of Contaminant Plume Response to DNAPL Source Zone Architecture and Depletion in Porous and Fractured Media (ER-1610) Objective Dense non-aqueous phase liquid
More informationResearch Areas: Subsurface Flow & Transport (modified after Springer; YZ s expertise/interest marked by X )
Research Areas: Subsurface Flow & Transport (modified after Springer; YZ s expertise/interest marked by X ) Classification Description A Fundamental of Transport in Porous Media A1 Phenomenological theory
More informationUPSCALING OF MULTIPHASE FLOW PROCESSES IN HETEROGENEOUS POROUS MEDIA
International Journal of Mechanical Engineering and Technology (IJMET) Volume 6, Issue 10, Oct 2015, pp. 161-170, Article ID: IJMET_06_10_018 Available online at http://www.iaeme.com/ijmet/issues.asp?jtype=ijmet&vtype=6&itype=10
More informationDOMINATING PROCESSES DURING DNAPL REMOVAL FROM THE SATURATED ZONE USING THERMAL WELLS
Lecture Session (LeS): E.11 Modelling DOMINATING PROCESSES DURING DNAPL REMOVAL FROM THE SATURATED ZONE USING THERMAL WELLS Uwe Hiester, Martina Müller, Oliver Trötschler, Hans-Peter Koschitzky, VEGAS
More informationREPORT ON GUIDANCE ON ASSESSMENT OF LIGHT NON- AQUEOUS PHASE LIQUID MOBILITY FOR SITE CLASSIFICATION PURPOSES IN BRITISH COLUMBIA.
REPORT ON GUIDANCE ON ASSESSMENT OF LIGHT NON- AQUEOUS PHASE LIQUID MOBILITY FOR SITE CLASSIFICATION PURPOSES IN BRITISH COLUMBIA Submitted to: BC Ministry of Environment 3 rd Floor, 2975 Jutland Road
More informationAssessing contaminant mass flow rates obtained by the integral groundwater investigation method by using the Virtual Aquifer approach
Calibration and Reliability in Groundwater Modelling: From Uncertainty to Decision Making (Proceedings of ModelCARE 25, The Hague, The Netherlands, June 25). IAHS Publ. 34, 26. 23 Assessing contaminant
More informationSurfactant enhanced aquifer remediation: application of mathematical models in the design and evaluation of a pilot-scale test
Groundwater Quality: Matured and Enhanced Restoration of Groundwater Pollution (Proceedings ofthe Groundwater Quality 2001 Conference held al Sheffield, UK. June 2001). IAHS Publ. no. 275. 2002. 303 Surfactant
More informationComment on the paper Field observations of soil moisture variability across scales by Famiglietti et al.
Comment on the paper Field observations of soil moisture variability across scales by Famiglietti et al. 5 Vereecken H. (1), T. Kamai (2), T. Harter (2), R. Kasteel (1), J.W. Hopmans (2), J.A. Huisman
More informationAn Introduction into Applied Soil Hydrology
Klaus Bohne An Introduction into Applied Soil Hydrology Preface Contents 1 Introduction: The soil as a reactor between the earth's atmosphere and the subsurface 1 2 Mechanical composition of mineral soils
More informationTHE PENNSYLVANIA STATE UNIVERSITY DEPARTMENT OF ENERGY AND MINERAL ENGINEERING GEOEE 408 CONTAMINANT HYDROLOGY
THE PENNSYLVANIA STATE UNIVERSITY DEPARTMENT OF ENERGY AND MINERAL ENGINEERING GEOEE 408 CONTAMINANT HYDROLOGY Mid-term Examination Tuesday March 4 th, 2008 75 minutes Answer all three questions. For water
More informationThis section presents the physics of multi-phase fluid flow in porous media.
Outline for LNAPL Training Module 1 The Basics Approximate time: 2 hours This section presents the physics of multi-phase fluid flow in porous media. I. What is NAPL? a. LNAPL (mix of hydrocarbons) b.
More informationCREATING AND EXPLICITLY CHARACTERIZING HETEROGENEITY AND NAPL MACRO-SCALE ENTRAPMENT AT THE INTERMEDIATE SCALE
CREATING AND EXPLICITLY CHARACTERIZING HETEROGENEITY AND NAPL MACRO-SCALE ENTRAPMENT AT THE INTERMEDIATE SCALE Gilbert R. Barth,, Tissa H. Illangasekare, Harihar Rajaram, and Mary C. Hill, Civil and Environmental
More informationUpscaling retardation factors in 2-D porous media
Calibration and eliability in Groundwater Modelling: From Uncertainty to Decision Making (Proceedings of ModelCAE 5, The Hague, The Netherlands, June 5). IAHS Publ.,. Upscaling retardation factors in -D
More informationLNAPL migration through soil cement barrier with and without flow condition
As. J. Energy Env. 2009, 10(03), 185-193 Asian Journal on Energy and Environment ISSN 1513-4121 Available online at www.asian-energy-journal.info Research Article LNAPL migration through soil cement barrier
More informationEffect of Heterogeneity in Capillary Pressure on Buoyancy Driven Flow of CO 2
CONFERENCE PROCEEDINGS Effect of Heterogeneity in Capillary Pressure on Buoyancy Driven Flow of CO 2 Ehsan Saadatpoor, Steven L. Bryant, Kamy Sepehrnoori The University of Texas at Austin SIXTH ANNUAL
More informationCoupled Simulation of DNAPL Infiltration and Dissolution in Three-Dimensional Heterogeneous Domains: Process Model Validation
The University of San Francisco USF Scholarship: a digital repository @ Gleeson Library Geschke Center Environmental Science College of Arts and Sciences 2013 Coupled Simulation of DNAPL Infiltration and
More informationCHAPTER 4: Risk Assessment Risk in Groundwater Contamination
CHAPTER 4: Risk Assessment Risk in Groundwater Contamination Instructor: Dr. Yunes Mogheir -١ Introduction: Water pollution is nowadays one of the most crucial environmental problems world-wide. Pollution
More informationAvailable online at ScienceDirect. Energy Procedia 63 (2014 ) GHGT-12
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 63 (2014 ) 3253 3260 GHGT-12 CO 2 leakage into shallow aquifers: Modeling CO 2 gas evolution and accumulation at interfaces of heterogeneity
More informationCENTRIFUGE MODELING OF LNAPL TRANSPORT
CENTRIFUGE MODELING OF LNAPL TRANSPORT IN PARTIALLY SATURATED SAND By G. Esposito, 1 Member, ASCE, H. G. B. Allersma, 2 and A. P. S. Selvadurai, 3 Member, ASCE Downloaded from ascelibrary.org by McGill
More informationThe Subsurface Flow and Transport Experimental Laboratory: A New Department of Energy User s Facility for Intermediate-Scale Experimentation
Hydrology days 2004 The Subsurface Flow and Transport Experimental Laboratory: A New Department of Energy User s Facility for Intermediate-Scale Experimentation M. Oostrom 1, T.W. Wietsma, and N.S. Foster
More informationThe Basics. Understanding the Behavior of Light Non- Aqueous Phase Liquids (LNAPLs) in the Subsurface February 2005
The Basics Understanding the Behavior of Light Non- Aqueous Phase Liquids (LNAPLs) in the Subsurface February 2005 1 Welcome. The Remediation Technologies Development Forum (RTDF) Non- Aqueous Phase Liquid
More informationFLOW AND DISTRIBUTION OF WATER IN UNSATURATED, FRACTIONALLY-WET POROUS MEDIA SYSTEMS
Proceedings of the 14 th International Conference on Environmental Science and Technology Rhodes, Greece, 3-5 September 2015 FLOW AND DISTRIBUTION OF WATER IN UNSATURATED, FRACTIONALLY-WET POROUS MEDIA
More informationSimulations of one and two phase flow in porous media: from the pore scale description to the macroscale studies
Simulations of one and two phase flow in porous media: from the pore scale description to the macroscale studies Félix Collin-Bastiani PhD - CERFACS (Master INP) Jacques Franc PhD INP Romain Guibert Research
More informationMultiphase flow and enhanced biodégradation of dense non-aqueous phase liquids: experimental design
Groundwater Quality: Natural and Enhanced Restoration of Groundwater Pollution (Proceedings of die Groundwater Quality 2001 Conference held at Sheffield. UK., June 2001). IAHS Publ. no. 275.2002. 265 Multiphase
More informationFluid Flow in Porous Media
Fluid Flow in Porous Media Petroleum Engineering 524 Fall 2010 Written by Thomas W. Engler, Ph.D., P.E. Professor of Petroleum Engineering New Mexico Tech Copyright 2010 Table of Contents Chapter 1 Introduction
More informationChapter 1. A Simulation and Decision Analysis Approach to Locating. DNAPLS in Subsurface Sediments. Erik K. Webb z. Abstract
Chapter 1 A Simulation and Decision Analysis Approach to Locating DNAPLS in Subsurface Sediments Brian Borchers y Stephen H. Conrad z Roger Cox z Robert J. Glass Jr. z Erik K. Webb z Abstract We present
More informationMineral Dissolution and Wormholing from a Pore-Scale Perspective
Mineral Dissolution and Wormholing from a Pore-Scale Perspective SUPRI-B Annual Meeting, May 1-2, 2017 CYPRIEN SOULAINE (csoulain@stanford.edu) HAMDI TCHELEPI Motivations Target: Understand and model complex
More informationTemperature Plume Migration in Aquifers: The necessary first step to geochemical evaluation of thermally-mobilized constituents
Temperature Plume Migration in Aquifers: The necessary first step to geochemical evaluation of thermally-mobilized constituents Nelson Molina Giraldo, PhD. Gordon MacMillan, P.Geol. Matrix Solutions Inc.
More informationNon-Aqueous Phase Liquids Young Environmental Professionals (YEPs) Presentation 23 February 2017
Non-Aqueous Phase Liquids Young Environmental Professionals (YEPs) Presentation 23 February 2017 Gordon Binkhorst Senior Hydrogeologist at ALTA Environmental Corp. 1998 - Present Brown University, B.S.
More informationGroundwater Modeling Guidance
Groundwater Modeling Guidance Richard J. Mandle Groundwater Modeling Program Michigan Department of Environmental Quality Draft 1.0 10/16/02 Executive Summary The use of groundwater models is prevalent
More informationCO 2 sequestration simulations: A comparison study between DuMuX and Eclipse
Universität Stuttgart - Institut für Wasserbau Lehrstuhl für Hydromechanik und Hydrosystemmodellierung Prof. Dr.-Ing. Rainer Helmig Diplomarbeit CO 2 sequestration simulations: A comparison study between
More informationPrediction of fluid occupancy in fractures using network modeling and x-ray microtomography. II: Results
Prediction of fluid occupancy in fractures using network modeling and x-ray microtomography. II: Results Mohammad Piri* Department of Chemical and Petroleum Engineering, University of Wyoming, Laramie,
More informationSimulation of Pumping Induced Groundwater Flow in Unconfined Aquifer Using Arbitrary Lagrangian-Eulerian Method
Simulation of Pumping Induced Groundwater Flow in Unconfined Aquifer Using Arbitrary Lagrangian-Eulerian Method Y. Jin *1, E. Holzbecher 1, and S. Ebneth 2 1 Applied Geology, Geoscience Centre, Georg-August-University
More informationEnhanced Delivery of Potassium Permanganate Using Horizontal Wells
Paper E-036, in: Bruce M. Sass (Conference Chair), Remediation of Chlorinated and Recalcitrant Compounds 2008. Proceedings of the Sixth International Conference on Remediation of Chlorinated and Recalcitrant
More informationPOROSITY, SPECIFIC YIELD & SPECIFIC RETENTION. Physical properties of
POROSITY, SPECIFIC YIELD & SPECIFIC RETENTION Porosity is the the ratio of the voids to the total volume of an unconsolidated or consolidated material. Physical properties of n = porosity as a decimal
More informationResearch Article Multiphase, Multicomponent Simulation for Flow and Transport during Polymer Flood under Various Wettability Conditions
Applied Mathematics Volume 213, Article ID 1167, 8 pages http://dx.doi.org/1.1155/213/1167 Research Article Multiphase, Multicomponent Simulation for Flow and Transport during Polymer Flood under Various
More informationDesiccation of Porous Media: Intermediate-Scale Flow Cell Experiments and Numerical Simulations
Hydrology Days 2008 Desiccation of Porous Media: Intermediate-Scale Flow Cell Experiments and Numerical Simulations M. Oostrom 1 Hydrology Group, Energy and Environment Directorate, Pacific Northwest National
More informationINVESTIGATION ON SURFACE AND SUBSURFACE FLUID MIGRATION: ENVIRONMENTAL IMPACT
Proceedings of the 13 th International Conference on Environmental Science and Technology Athens, Greece, 5-7 September 2013 INVESTIGATION ON SURFACE AND SUBSURFACE FLUID MIGRATION: ENVIRONMENTAL IMPACT
More informationCHAPTER 2. Objectives of Groundwater Modelling
CHAPTER 2 Objectives of Groundwater Modelling In the last two decades mathematical modelling techniques have increasingly proved their value in furthering the understanding of groundwater systems and,
More informationCHC-REMEDIATION OF THE SATURATED ZONE BY STEAM-AIR INJECTION
CHC-REMEDIATION OF THE SATURATED ZONE BY STEAM-AIR INJECTION Arne M. Färber (arne.faerber@iws.uni-stuttgart.de), Oliver Trötschler, Simon Steidinger, Steffen Ochs, Holger Class, Hans-Peter Koschitzky and
More informationEXPERIMENTAL INVESTIGATION OF FACTORS AFFECTING LABORATORY MEASURED RELATIVE PERMEABILITY CURVES AND EOR
SCA2017-052 1/10 EXPERIMENTAL INVESTIGATION OF FACTORS AFFECTING LABORATORY MEASURED RELATIVE PERMEABILITY CURVES AND EOR Sami M. Aboujafar Libyan Petroleum Institute, Tripoli - Libya This paper was prepared
More informationPropagation of water table waves in unconfined aquifers
Griffith Research Online https://research-repository.griffith.edu.au Propagation of water table waves in unconfined aquifers Author Jazayeri Shoushtari, Amir, Cartwright, Nick Published 203 Conference
More informationhal , version 1-18 Jan 2013
PROCEEDINGS, TOUGH Symposium 2012 Lawrence Berkeley National Laboratory, Berkeley, California, September 17-19, 2012 Author manuscript, published in "TOUGH Symposium 2012, United States (2012)" INFLUENCE
More informationThe Effect of Surface Texture on Evaporation, Infiltration and Storage Properties of Paved Surfaces
The Effect of Surface Texture on Evaporation, Infiltration and Storage Properties of Paved Surfaces M. Mansell* and F. Rollet School of Engineering and Science, University of the West of Scotland, Paisley
More informationFull-Scale Permanganate Remediation of a Solvent DNAPL Source Zone in a Sand Aquifer
Full-Scale Permanganate Remediation of a Solvent DNAPL Source Zone in a Sand Aquifer Beth L. Parker, Ph.D. University of Waterloo 1 Presented at the EPA Seminar: In Situ Treatment of Groundwater Contaminated
More informationALCOHOL FLUSHING IN LABORATORY EXPERIMENTS: IN-SITU REMEDIATION OF DNAPL CONTAMINATED GROUNDWATER
Theme D: Remediation Concepts & Technologies (PhD work: partly, by Klaas Heinrich) ALCOHOL FLUSHING IN LABORATORY EXPERIMENTS: IN-SITU REMEDIATION OF DNAPL CONTAMINATED GROUNDWATER Ulf Mohrlok 1, Klaas
More informationNUMERICAL SIMULATION OF GRAVITY INDUCED AGGREGATE MIGRATION IN LARGE IN-SITU SCC WALL CASTINGS
NUMERICAL SIMULATION OF GRAVITY INDUCED AGGREGATE MIGRATION IN LARGE IN-SITU SCC WALL CASTINGS Jon Spangenberg 1*, Jesper H. Hattel 1, Nicolas Roussel 2 and Mette R. Geiker 3 1 Department of Mechanical
More informationResidual trapping, solubility trapping and capillary pinning complement each other to limit CO[subscript 2] migration in deep saline aquifers
Residual trapping, solubility trapping and capillary pinning complement each other to limit CO[subscript 2] migration in deep saline aquifers The MIT Faculty has made this article openly available. Please
More informationLEAST COST SEARCH ALGORITHM FOR THE IDENTIFICATION OF A DNAPL SOURCE
LEAST COST SEARCH ALGORITHM FOR THE IDENTIFICATION OF A DNAPL SOURCE Z. DOKOU and G. F. PINDER Research Center for Groundwater Remediation Design, University of Vermont, Department of Civil and Environmental
More informationHYSTERESIS EFFECTS IN CAPILLARY PRESSURE, RELATIVE PERMEABILITY AND RESISTIVITY INDEX OF NORTH SEA CHALK
HYSTERESIS EFFECTS IN CAPILLARY PRESSURE, RELATIVE PERMEABILITY AND RESISTIVITY INDEX OF NORTH SEA CHALK M.T. Tweheyo, M.S. Talukdar and O. Torsæter Department of Petroleum Engineering and Applied Geophysics
More informationDocumentation of Groundwater Model
Documentation of Groundwater Model Groundwater Model The primary tool used to evaluate potential impacts on groundwater was a groundwater flow and solute transport model developed by San Benito County
More informationTechniques Used to Evaluate Potential LNAPL Mobility
Techniques Used to Evaluate Potential LNAPL Mobility Submitted by: David A. Fursevich, M.Sc., R. Donald Burnett, M.Sc., P.Eng., and Christopher P. Lach, P.Eng. Morrow Environmental Consultants Inc., Member
More informationDr. J. Wolters. FZJ-ZAT-379 January Forschungszentrum Jülich GmbH, FZJ
Forschungszentrum Jülich GmbH, FZJ ZAT-Report FZJ-ZAT-379 January 2003 Benchmark Activity on Natural Convection Heat Transfer Enhancement in Mercury with Gas Injection authors Dr. J. Wolters abstract A
More informationAre we over-complicating our models?
Are we over-complicating our models? Towards a more efficient and robust representation of sub-surface hydrological processes in Earth System Models Rafael Rosolem rafael.rosolem@bristol.ac.uk Soil Aquifer
More informationAssessing Soil Vapor Extraction Remediation Performance and Closure: A Review
Assessing Soil Vapor Extraction Remediation Performance and Closure: A Review - 12188 M.J. Truex, K.C. Carroll, and M. Oostrom Pacific Northwest National Laboratory, Richland, Washington ABSTRACT Soil
More informationDavid M. Tuck Westinghouse Savannah River Co., Savannah River Technology Center, Aiken SC, USA
WSRC-MS-99-00139 Time-Dependent Interfacial Properties and DNAPL Mobility David M. Tuck Westinghouse Savannah River Co., Savannah River Technology Center, Aiken SC, USA This document was prepared in conjunction
More informationMODELING WATER FLOW AND CONTAMINANT TRANSPORT IN SOILS AND GROUNDWATER USING THE HYDRUS COMPUTER SOFTWARE PACKAGES
MODELING WATER FLOW AND CONTAMINANT TRANSPORT IN SOILS AND GROUNDWATER USING THE HYDRUS COMPUTER SOFTWARE PACKAGES Instructor Dr. Jirka Simunek Department of Environmental Sciences University of California
More informationEnhanced Heap Leaching: II. Applications
Enhanced Heap Leaching: II. Applications Shlomo Orr and Velimir Vesselinov Abstract In our 1 st article (Orr, 2000), we investigated different flow and transport phenomena that could significantly reduce
More informationTransAT for Multiphase Flow in Pipes & Risers
TransAT for Multiphase Flow in Pipes & Risers Djamel Lakehal (July, 2013) ASCOMP GmbH, Zurich www.ascomp.ch; lakehal@ascomp.ch Outline Multiphase flow in pipelines Slug flow in horizontal pipes Slug flow
More informationHEAVY METAL RETENTION WITHIN A POROUS PAVEMENT STRUCTURE
HEAVY METAL RETENTION WITHIN A POROUS PAVEMENT STRUCTURE C. Dierkes, A. Holte, W. F. Geiger Department of Civil Engineering, Urban Water Management, University of Essen, Universitaetsstrasse 15, 45141
More informationMethods for predicting the spreading of steam below the water table during subsurface remediation
Click Here for Full Article WATER RESOURCES RESEARCH, VOL. 46,, doi:10.1029/2007wr006401, 2010 Methods for predicting the spreading of steam below the water table during subsurface remediation S. O. Ochs,
More informationApproaches to characterize the degree of water repellency
Journal of Hydrology 231 232 (2000) 61 65 www.elsevier.com/locate/jhydrol Review Approaches to characterize the degree of water repellency J. Letey*, M.L.K. Carrillo, X.P. Pang University of California,
More informationDesign of a passive hydraulic containment system using FEFLOW modelling
Design of a passive hydraulic containment system using FEFLOW modelling Rémi Vigouroux remi.vigouroux@arteliagroup.com Florence Lenhardt florence.lenhardt@arteliagroup.com Noëlle Doucet noelle.doucet@arteliagroup.com
More informationOverview on in-situ thermal remediation technologies
Overview on in-situ thermal remediation technologies Hans-Peter Koschitzky & Oliver Trötschler VEGAS, University of Stuttgart Experience Exchange Day - pilot test Stuttgart July 19 th, 2013 What you can
More informationDevelopment Of An Information Mapping Framework To Couple Surface Flow And Subsurface Hydraulic Models
City University of New York (CUNY) CUNY Academic Works International Conference on Hydroinformatics 8-1-2014 Development Of An Information Mapping Framework To Couple Surface Flow And Subsurface Hydraulic
More informationAnalysis Fraction Flow of Water versus Cumulative Oil Recoveries Using Buckley Leverett Method
Analysis Fraction Flow of Water versus Cumulative Oil Recoveries Using Buckley Leverett Method Reza Cheraghi Kootiani, and Ariffin Bin Samsuri International Science Index, Physical and Mathematical Sciences
More informationLNAPL Recovery Using Vacuum Enhanced Technology. Theresa Ferguson, R.G. June 2014
LNAPL Recovery Using Vacuum Enhanced Technology Theresa Ferguson, R.G. June 2014 Presentation Overview v LNAPL in the Subsurface v What is Vacuum Enhanced Technology? v Technology Description and Application
More informationJournal of Hazardous Materials
Journal of Hazardous Materials 165 (2009) 332 344 Contents lists available at ScienceDirect Journal of Hazardous Materials journal homepage: www.elsevier.com/locate/jhazmat An environmental screening model
More informationUNIVERSITY OF OSLO Department of Physics. Thesis for the degree Master of Science: Steady-State Upscaling of Polymer Flow. Kristin Asklund Larssen
UNIVERSITY OF OSLO Department of Physics Thesis for the degree Master of Science: Steady-State Upscaling of Polymer Flow Kristin Asklund Larssen June, Abstract The process of polymer flooding is widely
More informationLNAPL Volume and Mobility Estimation to Assess When to Stop Active Recovery. by Louis Sabourin, P.Eng. and David Tarnocai, P.Geo.
LNAPL Volume and Mobility Estimation to Assess When to Stop Active Recovery by Louis Sabourin, P.Eng. and David Tarnocai, P.Geo. Presentation Outline Introduction, Issue and Importance LNAPL Migration
More informationStylistic Modeling of Vadose Zone Transport Insight into Vapor Intrusion Processes
March 15, 2011 Workshop, Addressing Regulatory Challenges In Vapor Intrusion, A State of the Science Update Stylistic Modeling of Vadose Zone Transport Insight into Vapor Intrusion Processes Daniel B.
More informationEVALUATE HORIZONTAL WELL PRODUCTION PERFORMANCE IN HEAVY OIL RESERVOIRS
EVALUATE HORIZONTAL WELL PRODUCTION PERFORMANCE IN HEAVY OIL RESERVOIRS Hisham Kh. Ben Mahmud 1, Ian H. Sheng 2 and Mian U. Shafiq 1 1 Department of Petroleum Engineering, Curtin University, Sarawak, Malaysia
More informationSurfactant Enhanced DNAPL Source Zone Remediation: Results of a Field Demonstration and Implications for Bioavailability
Surfactant Enhanced DNAPL Source Zone Remediation: Results of a Field Demonstration and Implications for Bioavailability Linda M. Abriola The University of Michigan Ann Arbor Presented at: In Situ Treatment
More informationThe construction and use of digital rock cores in accurate flow modelling for CO 2 storage
The construction and use of digital rock cores in accurate flow modelling for CO 2 storage Catriona Reynolds, Sam Krevor Department of Earth Science & Engineering Imperial College London Laboratory reservoir
More informationReport on: Approaches and Methods for Evaluation of Light non- Aqueous - Hydrogeological Assessment Tools Project
Report on: Approaches and Methods for Evaluation of Light non- Aqueous - Hydrogeological Assessment Tools Project Submitted to: Ministry of Environment February 2006 Submitted by: Science Advisory Board
More informationA functional relationship between capillary pressure, saturation, and interfacial area as revealed by a pore-scale network model
WATER RESOURCES RESEARCH, VOL. 32, NO. 8, PAGES 2345 2358, AUGUST 1996 A functional relationship between capillary pressure, saturation, and interfacial area as revealed by a pore-scale network model Paul
More informationSEES 503 SUSTAINABLE WATER RESOURCES GROUNDWATER. Instructor. Assist. Prof. Dr. Bertuğ Akıntuğ
SEES 503 SUSTAINABLE WATER RESOURCES GROUNDWATER Instructor Assist. Prof. Dr. Bertuğ Akıntuğ Civil Engineering Program Middle East Technical University Northern Cyprus Campus SEES 503 Sustainable Water
More informationWATER RESOURCES RESEARCH, VOL. 40, W01511, doi: /2003wr001980, 2004
Dense nonaqueous phase liquid (DNAPL) source zone characterization: Influence of hydraulic property correlation on predictions of DNAPL infiltration and entrapment Lawrence D. Lemke 1 and Linda M. Abriola
More informationAPPLICABILITY OF CARMAN-KOZENY EQUATION FOR WELL CONSOLIDATED SAMPLES
SCA2013-080 1/6 APPLICABILITY OF CARMAN-KOZENY EQUATION FOR WELL CONSOLIDATED SAMPLES Babak Salimifard, University of Manitoba, Douglas Ruth, Faculty of Engineering, University of Manitoba, Craig Lindsay,
More informationPore-Scale micro-ct imaging: Non-Wetting Phase Cluster Size Distribution as a Function of Interfacial Tension
TCCS-6 2011, Trondheim, Norway Pore-Scale micro-ct imaging: Non-Wetting Phase Cluster Size Distribution as a Function of Interfacial Tension Apostolos Georgiadis 1,2, Steffen Berg 1, Geoffrey Maitland
More informationQuestions: What is calibration? Why do we have to calibrate a groundwater model? How would you calibrate your groundwater model?
Questions: What is calibration? Why do we have to calibrate a groundwater model? How would you calibrate your groundwater model? 13-11-2009 1 Uncertainties in groundwater models Conceptual model uncertainty
More informationSite Verified Contaminant Transport Model as a Mathematical Vehicle for Prevention of the Natural Aquifer Contamination
PUBLS. INST. GEOPHYS. POL. ACAD. SC., E-10 (406), 2008 Site Verified Contaminant Transport Model as a Mathematical Vehicle for Prevention of the Natural Aquifer Contamination Andrzej ANISZEWSKI Szczecin
More informationSTEAM ENHANCED EXTRACTION (SEE) AS INNOVATIVE APPROACH FOR TCE REMOVAL
STEAM ENHANCED (SEE) AS INNOVATIVE APPROACH FOR TCE REMOVAL Pavel Dusílek 1, Petr Kvapil 2, Kent S. Udell 3, Craig M. Hampson 4 1 AQUATEST a.s, Prague, Czech Republic, Phone: 420 234 607 151, Fax: 420
More informationInfiltration and Injection Sites and Example Experiments
Hanford 300 A IFC Infiltration and Injection Sites and Example Experiments Mark Rockhold PNNL April 16-19, 2007 Pacific Northwest National Laboratory U.S. Department of Energy Objectives Design a characterization
More informationHYDRUS-1D Computer Session
HYDRUS-1D Computer Session Nonequilibrium Water Flow and Solute Transport In this computer session with HYDRUS-1D we demonstrate the capability of HYDRUS-1D to simulate nonequilibrium water flow and solute
More informationMathematical Modeling of Permeation Grouting and Subsurface Barrier Performance
LBNL-39419 uc-2000 ERNESTORLANDO LAWRENCE BERKELEYNATIONAL LABORATORY Mathematical Modeling of Permeation Grouting and Subsurface Barrier Performance S. Finsterle, C.M. Oldenburg, A.L. James, K. Pruess,
More informationScience Olympiad. Mentor Invitational Hydrogeology ANSWER KEY. Name(s): School Name: Point Totals
Science Olympiad Mentor Invitational Hydrogeology ANSWER KEY Team Number: Raw Score: Rank: Name(s): School Name: Point Totals Possible Part 1: Groundwater Concepts and Vocabulary 30 Part 2: The Hydrogeology
More informationNUMERICAL MODELLING OF THE GROUNDWATER FLOW IN THE LEFT FLOODPLAIN AREA OF THE DANUBE RIVER
Proceedings of ALGORITMY 2002 Conference on Scientific Computing, pp. 237 244 NUMERICAL MODELLING OF THE GROUNDWATER FLOW IN THE LEFT FLOODPLAIN AREA OF THE DANUBE RIVER A. BALÁŽOVÁ, D. BAROKOVÁ, K. MIKULA,
More informationSOURCE ZONE AND MODEL PARAMETER INFLUENCE ON NAPL PLUME DEVELOPMENT IN A HETEROGENEOUS AQUIFER
SOURCE ZONE AND MODEL PARAMETER INFLUENCE ON NAPL PLUME DEVELOPMENT IN A HETEROGENEOUS AQUIFER MICHAEL MARINOVICH SUPERVISOR: DR. DAVID REYNOLDS SCHOOL OF ENVIRONMENTAL SYSTEMS ENGINEERING UNIVERSITY OF
More informationInvestigations in Geologic Carbon Sequestration: Multiphase Flow of CO 2 and Water in Reservoir Rocks. Annual Report 2012
Investigations in Geologic Carbon Sequestration: Multiphase Flow of CO 2 and Water in Reservoir Rocks Annual Report 2012 Sally M. Benson, Boxiao Li, Michael Krause, Sam Krevor, Chia-Wei Kuo, Ronny Pini
More informationInterpreting tracer data to forecast remedial performance
Groundwater Quality: Natural and Enhanced Restoration of Groundwater Pollution (Proceedings ofthe Groundwater Quality 2001 Conference held al Sheffield. UK. June 2001 ). 1AI IS Publ. no. 275. 2002. 11
More informationDEVELOPMENT OF A METHODOLOGY FOR THE DESIGN OF GROUNDWATER QUALITY SAMPLING NETWORKS
DEVELOPMENT OF A METHODOLOGY FOR THE DESIGN OF GROUNDWATER QUALITY SAMPLING NETWORKS Graciela S. Herrera 1, Joseph Guarnaccia 2 & George F. Pinder 3 Abstract - Groundwater quality sampling networks are
More informationRemediation of Brine Spills- What Goes Wrong Kerry Sublette
Remediation of Brine Spills- What Goes Wrong Kerry Sublette University of Tulsa Spills of produced water or brine on soil result in two types of damage: Excess salinity Creates an osmotic imbalance that
More informationEnhanced in Situ Leaching (EISLEACH) Mining by Means of Supercritical Carbon Dioxide - TOUGHREACT Modeling
An Interdisciplinary Response to Mine Water Challenges - Sui, Sun & Wang (eds) 2014 China University of Mining and Technology Press, Xuzhou, ISBN 978-7-5646-2437-8 Enhanced in Situ Leaching (EISLEACH)
More informationAccelerating The Dissolution of CO 2 in Aquifers
Accelerating The Dissolution of CO in Aquifers Yuri Leonenko 1, David W. Keith 1, Mehran Pooladi-Darvish 1, and Hassan Hassanzadeh 1 1 Department of Chemical and Petroleum Engineering, University of Calgary,
More informationChapter 1 INTRODUCTION. 1.1 Motivations and Background CO 2 Sequestration
Chapter 1 INTRODUCTION 1.1 Motivations and Background 1.1.1 CO 2 Sequestration The largely anthropogenic origins of the drivers of climate change over the last 100 years are, according to the fifth assessment
More information