DOCSLIB.ORG

Particle Retention in Porous Media: Applications to Water Injectivity Decline

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Particle Retention in Porous Media: Applications to Water Injectivity Decline NEI-N0--1119 NO9905156 Kjell Erik Wennberg Particle retention in porous media Applications to water injectivity decline r. €)< X DISCLAIMER Portions of this document may be illegible in electronic image products. Images are produced from the best available original document. Particle Retention in Porous Media: Applications to Water Injectivity Decline by Kjell Erik Wennberg A dissertation for the partial fulfillment of requirements for the degree of doktor ingenidr. Department of Petroleum Engineering and Applied Geophysics The Norwegian University of Science and Technology Trondheim, February 1998 i vtuuiu iuvv iv cApivaa my arnvcic giauiuuc iu my principal supervisor, nr. uie ror- saeter, professor at the Department of Petroleum Engineering and Applied Geophysics, NTNU, and my two technical advisors, Dr. Alex Hansen, professor at the Physics De­ partments at NTNU and Dr. Mukul M. Sharma, professor at the Department of Petro ­ leum and Geosystems Engineering at the University of Texas at Austin. Dedication To my dear wife May Helen and my three lovely children Miriam, Daniel and Sondre for their love, support and understanding throughout this study Acknowledgments I would like to express my sincere gratitude to my principal supervisor, Dr. Ole Tor- saeter, professor at the Department of Petroleum Engineering and Applied Geophysics, NTNU, and my two technical advisors, Dr. Alex Hansen, professor at the Physics De­ partments at NTNU and Dr. Mukul M. Sharma, professor at the Department of Petro ­ leum and Geosystems Engineering at the University of Texas at Austin. With their different backgrounds and research interests, they all gave invaluable contributions during my study. I have had many long and encouraging discussions with each of them. The results obtained in developing the WID simulator for water injectivity decline were accomplished when Dr. Sharma hosted me one year in his research group and gave me access to previous results obtained at UT. The significance of the multidi­ mensional nature of particle deposition in porous media was emphasized by Dr. Hansen who willingly shared his encyclopedia! knowledge in modem physics with me. And with his broad knowledge in Petroleum engineering, Dr. Torsaeter gave inputs all along and helped me to tie everything together in this document. Also, I want to thank my colleagues at IKU Petroleum Research, especially my section manager Erling Fjaer and our project secretary, Reidun Knutsen, for continued support throughout the study. I thank IKU for giving me the opportunity to fulfill this study. Finally, I am grateful to the Norwegian Research Council for the scholarships they provided for my study and the additional expenses related to me and my family staying one year in Austin, Tx. The stay greatly benefited both my scientific work and the per­ sonal development for myself and my family. Abstract The problem of migration and deposition of colloidal particles within porous media is studied theoretically and by computer simulations. Special emphasis is put on the pre­ diction of injectivity decline in water injection wells due to inherent particles in the in­ jection water. The study of particle deposition within porous media requires - among other things — a correct prediction of the deposition rate or filtration coefficient. A thorough review of the modeling approaches used in the past are combined with new ideas in order to arrive at an improved model for the prediction of the filtration coefficient. Also, a new way of determining the transition time, i.e. the time where the dominant deposition mechanism changes from internal deposition (deep filtration) to external cake formation is proposed. From this fundamental theory, equations for water injectivity decline pre­ dictions are given. Based on the theoretical model, a computer program named WID for water injec­ tivity decline predictions was developed. The simulator has a very friendly graphical user interface and runs on a PC under Windows-95. As input parameters, the simulator needs water quality and formation properties as well as injection rate/pressure and com ­ pletion information. Given these input parameters, WID is able to predict decline in vertical and horizontal injection wells with openhole, perforated and fractured comple ­ tions. WID has been tested against field data from different parts of the world. In some cases it matches the field data excellent and in most cases the match is fair. However, a poor match in a few cases indicates that more mechanisms can be responsible for injec­ tivity decline than those presently accounted for by the simulator. The second part of the study is devoted to a theoretical investigation of the multi­ dimensional nature of particle deposition in porous media. A velocity dependent form of the filtration coefficient is proposed. Based on this model, a 2-d numerical code is im­ plemented. It is shown that the permeability decline in a 2-d porous sample is a function of the velocity dependence in the filtration coefficient. In particular, it is shown that under certains conditions, low permeability bands orthogonal to the macroscopic direc­ - IV - tion of flow (localization) will form. For other conditions, high-permeability bands (wormholes) parallel to the direction of flow will form. Attempts to quantify wormhole and band formation are made. Table of Contents 1. Introduction 1 1.1 Objectives of Study 2 2. Literature Review 5 2.1 Particle Mobilization and Retention in Porous Media 5 2.1.1 Types of Particles 5 2.1.2 Overview of Retention Mechanisms. 5 2.1.3 Mobilization and Retention Forces 7 2.2 Formation Damage due to Various Operations 7 2.2.1 Formation Damage due to Drilling Operations 7 2.2.2 Formation Damage due to Completion Operations 8 2.2.3 Formation Damage During Production 9 2.2.4 Formation Damage during Water Injection 10 2.3 Continuum Models for Flow of Particle Suspensions in Porous Media 10 2.3.1 The Convection-Diffusion Equation 10 2.3.2 Approximations to the Convection Diffusion Equation 12 2.4 Water Injectivity Decline Models 13 2.4.1 Iwasaki's Model 13 2.4.2 The Barkman and Davidson Model 14 2.4.3 Eylander's Model 17 2.4.4 The Van Velzen Model 17 2.4.5 The Pang and Sharma Model 18 2.4.6 The IFP Models 19 2.5 General Formation Damage Models 20 2.5.1 The Parallel Pathway Model 20 2.5.2 Other Continuum Models 21 2.5.3 Network Models 21 2.6 Permeability Models 23 2.6.1 Darcy's law 23 2.6.2 Permeability Correlations 23 2.6.3 Experimental Relations 24 2.6.4 Kozeny-Carman Equation 24 2.6.5 Rumpf and Gupte's Equation 26 2.7 Permeability Reduction Models 27 2.7.1 Overview 27 2.7.2 Reduction Models based on the Rumpf and Gupte Equation 28 2.7.3 Reduction Models based on the Kozeny-Carman Equation 28 2.7.4 Other Models 30 3. The Filtration Coefficient and The Transition Time 31 3.1 Solutions for Non-changing Filtration Coefficient 31 3.2 Solutions for Variable Filtration Coefficient 33 3.3 Determination of the Initial Filtration Coefficient 35 3.3.1 Experimental Determination of X q 35 3.3.2 Effect of Fluid Velocity, u. 36 3.3.3 Effect of Injected Particle Size, d p. 39 - vi - 3.3.4 Effect of Formation Grain Size, d g . 40 3.3.5 Effect of Porosity, <|>. 40 3.3.6 Effect of Ion Concentration. 40 3.3.7 Determining Xo from Computer Simulations. 40 3.4 Variation of Filtration Coefficient with Specific Deposit 43 3.5 Simplified Approach for Water Injection Projects 45 3.5.1 The Transition Time 45 3.6 Parameter Estimation from Experimental Data 48 3.6.1 Summary 50 4. Effect of Completion Geometry on Water Injection Well Performance 51 4.1 How Resistance and IPR 51 4.2 How Resistance for Various Well Completions 52 4.3 How Resistance in Linear Geometry (e.g. Cores) 52 4.4 How Resistance in an Openhole Well 53 4.5 How Resistance Calculation for Fractured Wells 54 4.5.1 Raymond and Binder's Method 54 4.5.2 Muskat's Solution for an Infinite-Conductivity Fracture 56 4.5.3 Prats' Method 59 4.5.4 Comparison between the Solutions by R&B and Muskat 59 4.6 Plugging of Fractured Injection Wells 60 4.6.1 Internal Filtration 60 4.6.2 External Filtration 62 4.7 How Resistance around Perforations - Ellipsoidal Geometry 63 4.8 Decline in Horizontal Injection Wells 66 4.9 Layered Reservoirs 69 4.9.1 Unconnected Layers 69 4.9.2 Connected Layers 69 5. Development and Testing of a Simulator for Water Injectivity Decline 73 5.1 Introduction 73 5.2 Development of the Simulator - WID 73 5.2.1 Model Features and Limitations 73 5.2.2 Method of Solution 74 5.3 Comparison of WID with Core Hood Experiments by Todd et al. 77 5.3.1 Alundum Cores 78 5.3.2 Lochaline Sandstone Cores 79 5.3.3 North Sea Cores 81 6. Testing of WID against Field Data 83 6.1 Injectivity Decline in Water Injection Wells: An Offshore Gulf of Mexico Case Study 83 6.1.1 Water Injection Project History 83 6.1.2 Water Quality 85 6.1.3 Results 86 6.1.4 Discussion 90 6.1.5 Summary of Gulf of Mexico Case Study 91 6.2 Application of the WID Simulator to an Offshore Injection Well outside West-Africa with relevant Core How Tests. 101 - vii - 6.2.1 Field history 101 6.2.2 Analysis and treatment of the Injection Water 101 6.2.3 Core Flood tests 102 6.2.4 The TK-M1 Injection Well 108 6.3 Analysis of an Injection Well in the North Sea 111 7.
Load more

Recommended publications
  • Hydrogeology
    Hydrogeology Principles of Groundwater Flow Lecture 3 1 Hydrostatic pressure The total force acting at the bottom of the prism with area A is Dividing both sides by the area, A ,of the prism 1 2 Hydrostatic Pressure Top of the atmosphere Thus a positive suction corresponds to a negative gage pressure. The dimensions of pressure are F/L2, that is Newton per square meter or pascal (Pa), kiloNewton per square meter or kilopascal (kPa) in SI units Point A is in the saturated zone and the gage pressure is positive. Point B is in the unsaturated zone and the gage pressure is negative. This negative pressure is referred to as a suction or tension. 3 Hydraulic Head The law of hydrostatics states that pressure p can be expressed in terms of height of liquid h measured from the water table (assuming that groundwater is at rest or moving horizontally). This height is called the pressure head: For point A the quantity h is positive whereas it is negative for point B. h = Pf /ϒf = Pf /ρg If the medium is saturated, pore pressure, p, can be measured by the pressure head, h = pf /γf in a piezometer, a nonflowing well. The difference between the altitude of the well, H, and the depth to the water inside the well is the total head, h , at the well. i 2 4 Energy in Groundwater • Groundwater possess mechanical energy in the form of kinetic energy, gravitational potential energy and energy of fluid pressure. • Because the amount of energy vary from place-­­to-­­place, groundwater is forced to move from one region to another in order to neutralize the energy differences.
    [Show full text]
  • Continuum Mechanics of Two-Phase Porous Media
    Continuum mechanics of two-phase porous media Ragnar Larsson Division of material and computational mechanics Department of applied mechanics Chalmers University of TechnologyS-412 96 Göteborg, Sweden Draft date December 15, 2012 Contents Contents i Preface 1 1 Introduction 3 1.1 Background ................................... 3 1.2 Organizationoflectures . .. 6 1.3 Coursework................................... 7 2 The concept of a two-phase mixture 9 2.1 Volumefractions ................................ 10 2.2 Effectivemass.................................. 11 2.3 Effectivevelocities .............................. 12 2.4 Homogenizedstress ............................... 12 3 A homogenized theory of porous media 15 3.1 Kinematics of two phase continuum . ... 15 3.2 Conservationofmass .............................. 16 3.2.1 Onephasematerial ........................... 17 3.2.2 Twophasematerial. .. .. 18 3.2.3 Mass balance of fluid phase in terms of relative velocity....... 19 3.2.4 Mass balance in terms of internal mass supply . .... 19 3.2.5 Massbalance-finalresult . 20 i ii CONTENTS 3.3 Balanceofmomentum ............................. 22 3.3.1 Totalformat............................... 22 3.3.2 Individual phases and transfer of momentum change between phases 24 3.4 Conservationofenergy . .. .. 26 3.4.1 Totalformulation ............................ 26 3.4.2 Formulation in contributions from individual phases ......... 26 3.4.3 The mechanical work rate and heat supply to the mixture solid . 28 3.4.4 Energy equation in localized format . ... 29 3.4.5 Assumption about ideal viscous fluid and the effective stress of Terzaghi 30 3.5 Entropyinequality ............................... 32 3.5.1 Formulation of entropy inequality . ... 32 3.5.2 Legendre transformation between internal energy, free energy, entropy and temperature 3.5.3 The entropy inequality - Localization . .... 34 3.5.4 Anoteontheeffectivedragforce . 36 3.6 Constitutiverelations.
    [Show full text]
  • Coupled Multiphase Flow and Poromechanics: 10.1002/2013WR015175 a Computational Model of Pore Pressure Effects On
    PUBLICATIONS Water Resources Research RESEARCH ARTICLE Coupled multiphase flow and poromechanics: 10.1002/2013WR015175 A computational model of pore pressure effects on Key Points: fault slip and earthquake triggering New computational approach to coupled multiphase flow and Birendra Jha1 and Ruben Juanes1 geomechanics 1 Faults are represented as surfaces, Department of Civil and Environmental Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, capable of simulating runaway slip USA Unconditionally stable sequential solution of the fully coupled equations Abstract The coupling between subsurface flow and geomechanical deformation is critical in the assess- ment of the environmental impacts of groundwater use, underground liquid waste disposal, geologic stor- Supporting Information: Readme age of carbon dioxide, and exploitation of shale gas reserves. In particular, seismicity induced by fluid Videos S1 and S2 injection and withdrawal has emerged as a central element of the scientific discussion around subsurface technologies that tap into water and energy resources. Here we present a new computational approach to Correspondence to: model coupled multiphase flow and geomechanics of faulted reservoirs. We represent faults as surfaces R. Juanes, embedded in a three-dimensional medium by using zero-thickness interface elements to accurately model [email protected] fault slip under dynamically evolving fluid pressure and fault strength. We incorporate the effect of fluid pressures from multiphase flow in the mechanical
    [Show full text]
  • Predicting Porosity, Permeability, and Tortuosity of Porous Media from Images by Deep Learning Krzysztof M
    www.nature.com/scientificreports OPEN Predicting porosity, permeability, and tortuosity of porous media from images by deep learning Krzysztof M. Graczyk* & Maciej Matyka Convolutional neural networks (CNN) are utilized to encode the relation between initial confgurations of obstacles and three fundamental quantities in porous media: porosity ( ϕ ), permeability (k), and tortuosity (T). The two-dimensional systems with obstacles are considered. The fuid fow through a porous medium is simulated with the lattice Boltzmann method. The analysis has been performed for the systems with ϕ ∈ (0.37, 0.99) which covers fve orders of magnitude a span for permeability k ∈ (0.78, 2.1 × 105) and tortuosity T ∈ (1.03, 2.74) . It is shown that the CNNs can be used to predict the porosity, permeability, and tortuosity with good accuracy. With the usage of the CNN models, the relation between T and ϕ has been obtained and compared with the empirical estimate. Transport in porous media is ubiquitous: from the neuro-active molecules moving in the brain extracellular space1,2, water percolating through granular soils3 to the mass transport in the porous electrodes of the Lithium- ion batteries4 used in hand-held electronics. Te research in porous media concentrates on understanding the connections between two opposite scales: micro-world, which consists of voids and solids, and the macro-scale of porous objects. Te macroscopic transport properties of these objects are of the key interest for various indus- tries, including healthcare 5 and mining6. Macroscopic properties of the porous medium rely on the microscopic structure of interconnected pore space. Te shape and complexity of pores depend on the type of medium.
    [Show full text]
  • Theoretical and Computational Poromechanics from 14:30 to 16:30 Via Ferrata 3, Pavia Carlo Callari, Università Del Molise
    DIPARTIMENTO DI INGEGNERIA CIVILE E ARCHITETTURA DOTTORATO IN INGEGNERIA CIVILE E ARCHITETTURA UNIVERSITÀ DEGLI STUDI DI PAVIA SHORT COURSE ON JUNE, 22, 23, 24, 25, 26 THEORETICAL AND COMPUTATIONAL POROMECHANICS FROM 14:30 TO 16:30 VIA FERRATA 3, PAVIA CARLO CALLARI, UNIVERSITÀ DEL MOLISE Detailed Outline MOTIVATIONS: The role of poromechanics in civil, environmental and medical engineering: problem analysis and material modeling PART 1: Mechanics of fully saturated porous media with compressible phases Biot's thermodynamics and constitutive equations: • Thermodynamics of barotropic and inviscid fluids • Fluid mass balance in the porous medium • Local forms of first and second principles • Clausius-Duhem inequality • Transport laws • General form of constitutive laws for the porous medium • Linear poroelasticity and extension to poro-elastoplasticity • Fluid-to-solid and solid-to-fluid limit uncoupled influences • Strain-dependent permeability models Finite element formulations for porous media: • Boundary conditions for mechanical and fluid-flow fields • Formulation of mechanical and fluid-flow problems • Extension to non-linear response Research applications to saturated porous media: • 2D response of dam and rock mass to reservoir operations • Tunnel face stability • Poroelastic damage and applications to hydrocarbon wells PART 2: Extension to multiphase fluids Biot's thermodynamics and poroelastic laws: • General hyperelastic laws for three-phase porous media • Generalized Darcy law • Constitutive equations from the theory of mixtures
    [Show full text]
  • POROMECHANICS of POROUS and FRACTURED RESERVOIRS Jishan Liu, Derek Elsworth KIGAM, Daejeon, Korea May 15-18, 2017
    POROMECHANICS OF POROUS AND FRACTURED RESERVOIRS Jishan Liu, Derek Elsworth KIGAM, Daejeon, Korea May 15-18, 2017 1. Poromechanics – Flow Properties (Jishan Liu) 1:1 Reservoir Pressure System – How to calculate overburden stress and reservoir pressure Day 11 1:2 Darcy’s Law – Permeability and its changes, reservoir classification 1:3 Mass Conservation Law – flow equations 1:4 Steady-State Behaviors – Solutions of simple flow problems 1:5 Hydraulic Diffusivity – Definition, physical meaning, and its application in reservoirs 1:6 Rock Properties – Their Dependence on Stress Conditions Day 1 ……………………………………………………………………………………………………………... 2. Poromechanics – Fluid Storage Properties (Jishan Liu) 2:1 Fluid Properties – How they change and affect flow Day 2 2:2 Mechanisms of Liquid Production or Injection 2:3 Estimation of Original Hydrocarbons in Place 2:4 Estimation of Ultimate Recovery or Injection of Hydrocarbons 2:5 Flow – Deformation Coupling in Coal 2:6 Flow – Deformation Coupling in Shale Day 2 ……………………………………………………………………………………………………………... 3. Poromechanics –Modeling Porous Medium Flows (Derek Elsworth) 3:1 Single porosity flows - Finite Element Methods [2:1] Lecture Day 3 3:2 2D Triangular Constant Gradient Elements [2:3] Lecture 3:3 Transient Behavior - Mass Matrices [2:6] Lecture 3:4 Transient Behavior - Integration in Time [2:7] Lecture 3:5 Dual-Porosity-Dual-Permeability Models [6:1] Lecture Day 3 ……………………………………………………………………………………………………………... 4. Poromechanics –Modeling Coupled Porous Medium Flow and Deformation (Derek Elsworth) 4:1 Mechanical properties – http://www.ems.psu.edu/~elsworth/courses/geoee500/GeoEE500_1.PDF Day 4 4:2 Biot consolidation – http://www.ems.psu.edu/~elsworth/courses/geoee500/GeoEE500_1.PDF 4:3 Dual-porosity poroelasticity 4:4 Mechanical deformation - 1D and 2D Elements [5:1][5:2] Lecture 4:5 Coupled Hydro-Mechanical Models [6:2] Lecture Day 4 ……………………………………………………………………………………………………………..
    [Show full text]
  • Porous Medium Typology Influence on the Scaling Laws of Confined
    water Article Porous Medium Typology Influence on the Scaling Laws of Confined Aquifer Characteristic Parameters Carmine Fallico , Agostino Lauria * and Francesco Aristodemo Department of Civil Engineering, University of Calabria, via P. Bucci, Cubo 42B, 87036 Arcavacata di Rende (CS), Italy; [email protected] (C.F.); [email protected] (F.A.) * Correspondence: [email protected]; Tel.: +39-(0)9-8449-6568 Received: 4 March 2020; Accepted: 17 April 2020; Published: 19 April 2020 Abstract: An accurate measurement campaign, carried out on a confined porous aquifer, expressly reproduced in laboratory, allowed the determining of hydraulic conductivity values by performing a series of slug tests. This was done for four porous medium configurations with different granulometric compositions. At the scale considered, intermediate between those of the laboratory and the field, the scalar behaviors of the hydraulic conductivity and the effective porosity was verified, determining the respective scaling laws. Moreover, assuming the effective porosity as scale parameter, the scaling laws of the hydraulic conductivity were determined for the different injection volumes of the slug test, determining a new relationship, valid for coarse-grained porous media. The results obtained allow the influence that the differences among the characteristics of the porous media considered exerted on the scaling laws obtained to be highlighted. Finally, a comparison was made with the results obtained in a previous investigation carried out at the field scale. Keywords: hydraulic conductivity; effective porosity; scaling behavior; grain size analysis. 1. Introduction Among the parameters’ characterizing aquifers, certainly the hydraulic conductivity is one of the most meaningful and useful for the description of the flow and mass transport phenomena in porous media.
    [Show full text]
  • The Porous Medium Equation
    THE POROUS MEDIUM EQUATION Mathematical theory by JUAN LUIS VAZQUEZ¶ Dpto. de Matem¶aticas Univ. Auton¶omade Madrid 28049 Madrid, SPAIN ||||||||- i To my wife Mariluz Preface The Heat Equation is one of the three classical linear partial di®erential equations of second order that form the basis of any elementary introduction to the area of Partial Di®erential Equations. Its success in describing the process of thermal propagation has known a per- manent popularity since Fourier's essay Th¶eorieAnalytique de la Chaleur was published in 1822, [237], and has motivated the continuous growth of mathematics in the form of Fourier analysis, spectral theory, set theory, operator theory, and so on. Later on, it contributed to the development of measure theory and probability, among other topics. The high regard of the Heat Equation has not been isolated. A number of related equa- tions have been proposed both by applied scientists and pure mathematicians as objects of study. In a ¯rst extension of the ¯eld, the theory of linear parabolic equations was devel- oped, with constant and then variable coe±cients. The linear theory enjoyed much progress, but it was soon observed that most of the equations modelling physical phenomena without excessive simpli¯cation are nonlinear. However, the mathematical di±culties of building theories for nonlinear versions of the three classical partial di®erential equations (Laplace's equation, heat equation and wave equation) made it impossible to make signi¯cant progress until the 20th century was well advanced. And this observation applies to other important nonlinear PDEs or systems of PDEs, like the Navier-Stokes equations.
    [Show full text]
  • Manual of Computing and Modeling Techniques and Their Application To
    I . ' ' '. : ... ... ·.· I '· ·.. ' . I : ·_ .. ·· .. -.~.· ~[i)=(Q)ll~~·-··:.· , . : . I . .. -~ . I .. .• . : . ' . ' ~- ' . .. ·· ... ' ~·. ·._· ·.~-- ·.• . •. ' ·. ' ·--- ..... _ .· .. --- -- . ~ ... ' ''. ·i' : . ; . ' .. ' .. ·' , .• . ,' ' ·' . I . •· I : . .. ' . • ' • . • f _: . '" I .. · " . .· . .. ·. I . ; - ; _:;. :. ~ ..· .... : . II' ..·• .·.;. ; - ' . .. ·, I ...~ I ,• . I ' ..' '. ~ . '.. ' .' '·,( ;\ :: •' . ,. !. .. '·: .::- · : ~~~~~~ w~u~!Rl rcorMH%1J~~~~o~ - .· • :· · . ·· .JJ~INIQJJ~IR1'1f·' ll~l5~ . · ' · .. ·, . I . .' ' •• • '·. ••• • • • ' .¥ > ' - .. ,. .. ·' ' . I . LIBRARY --;-­ TEXAS. WATER DEVElOPMENT BVARll ' AUS liN, TEXAS I TEXAS WATER COMMISSION . Joe D. Carter, ~Chairman 0. F. Dent, C~issioner · H. A. Beckwith, Commissioner . ·~:· :.. ; : ; .· _...:._ _:_ ._, , Report LD-0165 __ __ , f ;, I '·'·· ,, ., ' .. ~- •.•' .. ·.::. MANUAL Of ,, COMPUTING AND MODELING TECHNIQUES . I . ' . ·. ' AND THEIR . .. ', :~ ' .. ~·r -(· I ' \• '•' '' ·:. ·: ··:-_ i'.;'' '.,APPLICATION TO HYDROLOGic· STUDIES.:. 1·'" ',' · -:: -: · .. i;. ~ • ' i · .. .-. c , r. <' '_:J' ··i . '~ ;- ., .: '· .. ,-, . ·., . ·-~ .. ,.l ·, •1'' • ., .- , .. By J; Russell Mount :. -· .. ·, ·' --,__ __ -------- ., •.- January· 965 . PREFJICE This manual is written for employees! of the Texas Water Commission .eng_a~ed. in problems 'of water resources, with the hope that by gaining a basic under-·· I . standing of computers and models they rna~ develop new ideas which will lead to increased efficiency and output of worlt:···-Much of· the manual
    [Show full text]
  • A Generalized Lattice Boltzmann Model for Flow Through Tight Porous Media with Klinkenberg’S Effect
    A generalized lattice Boltzmann model for flow through tight porous media with Klinkenberg’s effect Li Chen a,b, Wenzhen Fang a, Qinjun Kang b, Jeffrey De'Haven Hyman b, Hari S Viswanathan b, Wen‐Quan Tao a a: Key Laboratory of Thermo‐Fluid Science and Engineering of MOE, School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an, Shaanxi, 710049, China b: Computational Earth Science, EES‐16, Earth and Environmental Sciences Division, Los Alamos National Laboratory, Los Alamos, New Mexico, 87544, USA Abstract: Gas slippage occurs when the mean free path of the gas molecules is in the order of the characteristic pore size of a porous medium. This phenomenon leads to the Klinkenberg’s effect where the measured permeability of a gas (apparent permeability) is higher than that of the liquid (intrinsic permeability). A generalized lattice Boltzmann model is proposed for flow through porous media that includes Klinkenberg’s effect, which is based on the model of Guo et al. (Z.L. Guo et al., Phys.Rev.E 65, 046308 (2002)). The second‐order Beskok and Karniadakis‐Civan’s correlation (A. Beskok and G. Karniadakis, Microscale Thermophysical Engineering 3, 43‐47 (1999), F. Civan, Transp Porous Med 82, 375‐384 (2010)) is adopted to calculate the apparent permeability based on intrinsic permeability and Knudsen number. Fluid flow between two parallel plates filled with porous media is simulated to validate model. Simulations performed in a heterogeneous porous medium with components of different porosity and permeability indicate that the Klinkenberg’s effect plays significant role on fluid flow in low‐permeability porous media, and it is more pronounced as the Knudsen number increases.
    [Show full text]
  • Accelerated Perturbation Boundary Element Model for Flow Problems in Heterogeneous Reservoirs
    ACCELERATED PERTURBATION BOUNDARY ELEMENT MODEL FOR FLOW PROBLEMS IN HETEROGENEOUS RESERVOIRS a dissertation submitted to the department of petroleum engineering and the committee on graduate studies of stanford university in partial fulfillment of the requirements for the degree of doctor of philosophy By Kozo Sato June 1992 c Copyright 1992 by Kozo Sato All Rights Reserved ii I certify that I have read this thesis and that in my opin- ion it is fully adequate, in scope and in quality, as a dissertation for the degree of Doctor of Philosophy. Roland N. Horne (Principal Adviser) I certify that I have read this thesis and that in my opin- ion it is fully adequate, in scope and in quality, as a dissertation for the degree of Doctor of Philosophy. Khalid Aziz I certify that I have read this thesis and that in my opin- ion it is fully adequate, in scope and in quality, as a dissertation for the degree of Doctor of Philosophy. Henry J. Ramey, Jr. Approved for the University Committee on Graduate Studies: Dean of Graduate Studies iii Abstract The boundary element method (BEM), successfully applied to fluid flow problems in porous media, owes its elegance, the computational efficiency and accuracy, to the existence of the fundamental solution (free-space Green’s function) for the governing equation. If there is no such solution in a closed form, the BEM cannot be applied, and, unfortunately, this is the case for flow problems in heterogeneous media. To overcome this difficulty, the governing equation is decomposed into various order perturbation equations, for which the free-space Green’s function can be found.
    [Show full text]
  • Introduction to the Porous Media Flow Module
    INTRODUCTION TO Porous Media Flow Module Introduction to the Porous Media Flow Module © 1998–2019 COMSOL Protected by patents listed on www.comsol.com/patents, and U.S. Patents 7,519,518; 7,596,474; 7,623,991; 8,457,932; 8,954,302; 9,098,106; 9,146,652; 9,323,503; 9,372,673; 9,454,625; and 10,019,544. Patents pending. This Documentation and the Programs described herein are furnished under the COMSOL Software License Agreement (www.comsol.com/comsol-license-agreement) and may be used or copied only under the terms of the license agreement. COMSOL, the COMSOL logo, COMSOL Multiphysics, COMSOL Desktop, COMSOL Compiler, COMSOL Server, and LiveLink are either registered trademarks or trademarks of COMSOL AB. All other trademarks are the property of their respective owners, and COMSOL AB and its subsidiaries and products are not affiliated with, endorsed by, sponsored by, or supported by those trademark owners. For a list of such trademark owners, see www.comsol.com/ trademarks. Version: COMSOL 5.5 Contact Information Visit the Contact COMSOL page at www.comsol.com/contact to submit general inquiries, contact Technical Support, or search for an address and phone number. You can also visit the Worldwide Sales Offices page at www.comsol.com/contact/offices for address and contact information. If you need to contact Support, an online request form is located at the COMSOL Access page at www.comsol.com/support/case. Other useful links include: • Support Center: www.comsol.com/support • Product Download: www.comsol.com/product-download • Product Updates: www.comsol.com/support/updates • COMSOL Blog: www.comsol.com/blogs • Discussion Forum: www.comsol.com/community • Events: www.comsol.com/events • COMSOL Video Gallery: www.comsol.com/video • Support Knowledge Base: www.comsol.com/support/knowledgebase Part number: CM024802 Contents Introduction .
    [Show full text]