Jorge Luis Rodriguez

Total Page:16

File Type:pdf, Size:1020Kb

Jorge Luis Rodriguez NUMERICAL SIMULATION OF WATER FLOW THROUGH UNSATURATED SOIL IN VERTICAL AND INCLINED LAYERED SYSTEMS by Jorge Luis Rodriguez A thesis submitted in partial fulfillment of the requirements for the degree of Master in Science In Geotechnical Engineering Department of Civil and Environmental Engineering University of Alberta © Jorge Luis Rodriguez, 2016 Abstract In past decades, there were discussions explaining the mechanisms and factors that influence water flow through soil in unsaturated conditions. This phenomenon is of major importance to the mining industry as water is one of the trigger mechanisms for Acid Rock Generation (ARD) generation. In 1999, Newman evaluated water flow in unsaturated conditions for two column experiments using a vertical layer system of sandy materials and waste rock from the Golden Sunlight mine. Later in 2009, Andrina analyzed water flow using waste rock from the Grasberg Mine in three Meso-scale experiments to understand flow mechanisms for incline layers of waste rock. The current study focuses on modeling the two column experiments by Newman and the three Meso-scale panels by Andrina to analyze the mechanisms controlling water flow in unsaturated soils. Additionally, one of the models is evaluated under three additional materials to compare the effect of different hydraulic properties in an incline layering system. With finite element methods, the experiments are modeled and calculated under equal boundary conditions. The use of climate boundaries recreates the precipitation flux and head pressure boundary to generate the suction from the discharge points. The models are run using SvFlux software, which runs an automatic mesh refinement algorithm and solves the Partial Differential Equations (PDEs) using FlexPDE. The models are validated based on the correlation of the discharge volume from each experiment. The results from the models describe the profile changes of head pressure (hp), flux, flow paths, and matric suction to describe the mechanism of flow in the unsaturated conditions. The validation of the models was achieved through back analyzing the different tests in the experiments, due to the low correlation of the model using the laboratory properties and the experimental results. The back analysis of the material is focused on the air entry value, and the saturated hydraulic conductivity as a means to change the unsaturated flow. The result from the research displays similarities and disagreements between the models and experiments. The models showed that the measured discharge from a particular material does not represent the preferential flow path of water, as a small gradient in the pressure distribution can generate breakthrough at the base of the system. Page | ii Acknowledgements I start to thank my parents for the enormous effort to help me achieve my dreams, giving me their unconditional support in every step of way. They always believe in me and push me to higher goals. I would also like to thank my wife Diana, the most special woman I've ever met, helping me in every day of my Masters, giving me the company, affection, and support. Despite the distance during all this time our love grew unimaginably. Thank you for always giving company and helping me in the darkest days. I also want to give a special acknowledgement to my supervisor Dr. Ward Wilson without whom this project would not have been possible, a great mentor who taught me great things. Thank you very much for giving me the unique opportunity to work under your wing, Thank you very much for believing and trusting me. A very special appreciation to Dr. Macciotta for helping me, even though he did not have any obligation, but gave me guidance nonetheless to proceed with my research in the different stages of the investigation, and also helped me evaluate the early edition of the thesis. Thank you Dr. Andrina for taking the time to clarify the procedures and giving more detailed information from the experiment that she conducted in her PhD. Thank you to Vivian Giang for all the grammatical tips and quick revisions you provided me. Sally and Lorene, thank you very much for everything you did for me, helping me with all department procedures and letters in the department. I want to thank all the people who were part of my academic process, all the friends I made that in one way or another gave me their advice, particularly Ahlam Abdulnabi for all the talk about unsaturated flow mechanics. Finally, I would like to thank the University of Alberta and the Department of Civil & Environmental Engineering. I feel an honor to become a former student at this prestigious institution. Page | iii TABLE OF CONTENTS ABSTRACT II ACKNOWLEDGEMENTS III CHAPTER 1 INTRODUCTION 1 1.1. BACKGROUND 1 1.2. RESEARCH OBJECTIVES 3 CHAPTER 2 LITERATURE REVIEW 4 2.1. WATER IN MINE WASTE ROCK 4 2.2. WASTE ROCK EMBANKMENT IN GOLDEN SUNLIGHT MINE 6 2.3. FLOW PATH OF WATER INFILTRATION FROM THE SURFACE 10 2.4. LATERAL FLOW IN A VERTICAL COLUMN 12 2.4.1. LABORATORY TEST PROGRAM 12 2.4.1. EXPERIMENT CONDITIONS 15 2.4.2. MATERIAL PROPERTIES 16 2.4.3. EXPERIMENT RESULTS 20 2.5. WATER FLOW MODEL IN WASTE ROCK EMBANKMENTS 23 2.6. FLOW IN INCLINE LAYERING SYSTEM 28 2.6.1. SINGLE LAYER SYSTEM 28 2.6.2. DOUBLE LAYER SYSTEM 30 2.6.3. MULTIPLE LAYER SYSTEM 32 2.7. SUMMARY 45 CHAPTER 3 THEORY AND UNSATURATED FLOW 47 3.1. UNSATURATED SOILS 47 Page | iv 3.1.1. SOIL WATER CHARACTERISTIC CURVE 49 3.1.1. ESTIMATION OF THE HYDRAULIC CONDUCTIVITY CURVE 53 3.2. WATER FLOW IN UNSATURATED CONDITION 55 3.3. NUMERICAL MODELING FOR WATER FLOW IN UNSATURATED SOILS 57 3.3.1. FLOW LAWS 59 3.3.2. MASS BALANCE 60 3.3.3. PDE FOR UNSATURATED STEADY STATE WATER FLOW 61 3.4. SUMMARY 62 CHAPTER 4 NUMERICAL MODEL OF WATER FLOW IN UNSATURATED COLUMN LAYERED SYSTEM 64 4.1. INTRODUCTION 64 4.2. OBJECTIVES 64 4.3. LABORATORY COLUMN EXPERIMENT 65 4.3.1. MATERIAL CHARACTERISTICS 68 4.3.2. EXPERIMENT RESULTS 74 4.3.1. SEEP/W SIMULATION RESULTS 76 4.4. NUMERICAL MODEL 79 4.4.1. MODELLING METHODOLOGY 80 4.4.2. MODELING ASSUMPTIONS 84 4.4.3. THEORY ADOPTED BY THE NUMERICAL INTEGRATION FOR A COLUMN LAYERED SYSTEM 84 4.4.4. MODEL BOUNDARY CONDITIONS 87 4.5. WATER FLOW SIMULATION RESULTS 88 4.5.1. WATER FLOW IN SANDY MATERIAL 88 4.5.2. WATER FLOW IN WASTE ROCK MATERIAL 103 4.5.3. INTERNAL RESPONDS FOR SENSITIVITY MODELS 107 4.5.4. HIGHEST CORRELATION CONDITIONS FOR THE COLUMN EXPERIMENT 113 4.6. ANALYSIS AND DISCUSSION 115 4.7. CONCLUSIONS 118 CHAPTER 5 NUMERICAL MODEL OF WATER FLOW MODEL IN INCLINE LAYERED SYSTEM 120 Page | v 5.1. INTRODUCTION 120 5.2. OBJECTIVES 120 5.3. MESO-SCALE PANEL EXPERIMENT 121 5.3.1. MATERIAL PROPERTIES 126 5.3.2. EXPERIMENT RESULTS 128 5.4. NUMERICAL MODEL 131 5.4.1. MODEL METHODOLOGY 132 5.4.2. MODEL ASSUMPTION 135 5.4.3. THEORY ADOPTED FOR THE NUMERICAL INTEGRATION FOR A COLUMN LAYERED SYSTEM 135 5.4.4. MODEL BOUNDARY CONDITIONS 138 5.5. SIMULATION OF PANEL-1 140 5.5.1. SIMULATION APPLYING EXPERIMENTAL CONDITIONS AND MEASURE PROPERTIES 140 5.5.2. SIMULATION OF PANEL-1 WITH CALIBRATED PROPERTIES 144 5.6. SIMULATION OF PANEL-2 152 5.7. SIMULATION OF PANEL-3 157 5.8. ANALYSIS AND DISCUSSION 167 5.8.1. MODELLING OF MESO-SCALE PANELS 168 5.8.2. NUMERICAL MODEL COMPARISON 175 5.9. CONCLUSIONS 177 CHAPTER 6 CONCLUSIONS 180 6.1. FURTHER RESEARCH 181 APPENDIX 183 A. Back Analysis of ksat for BCS in Column-1 183 B. Back Analysis of Ksat for SS in Column-1 186 C. Back Analysis of AEV for BCS in Column-1 188 D. Back Analysis of Ksat For WR in Column-2 190 E. Water Flow Transfer in Sandy Soil 192 F. Water Flow Transfer in Waste rock 198 Page | vi G. Pressure Profiles In Column-1 Using Sandy Soils 200 H. Pressure Profiles In Column-2 Using Waste Rock 203 I. Profile Results from Simulation in Column-1 204 J. Profile Results from Simulation in Column-2 208 K. Calibration of Acid Rock in Panel-1 209 L. Simulation Profiles in Panel-1 218 M. Simulation Profiles in Panel-2 221 N. Simulation Profiles in Panel-3 223 O. Simulation of Panel-1 Under Different Material Properties 226 REFERENCES 246 Page | vii LIST OF FIGURES Figure 2.1 ARD Trigger Mechanisms (INAP, 2009, p. Chapter 4) ____________________________________ 4 Figure 2.2 Types of Hydrostratigraphy in Rock Pile Embankments __________________________________ 5 Figure 2.3 Waste Rock Benches (Wilson, 2003) _________________________________________________ 7 Figure 2.4 Waste Rock Embankment Structure (Wilson, 2001) _____________________________________ 8 Figure 2.5 Water Cycle in a Waste Rock Embankment (Herasymuik, 1996) ___________________________ 9 Figure 2.6 Column Design Experiment for Horton and Hawkins (1965) ______________________________ 11 Figure 2.7 Scheme of Laboratory Column 1 after Newman (1999) _________________________________ 14 Figure 2.8 Scheme of Laboratory Column 2 after Newman (1999) _________________________________ 14 Figure 2.9 Soil Water Characteristic Curve for Materials in Column-1 Wilson (1990), Swanson (1991), Newman (1999) _________________________________________________________________________ 17 Figure 2.10 Hydraulic Conductivity Curve for Material in Column-1 (Newman, 1999) __________________ 18 Figure 2.11 Soil Water Characteristic Curve for Materials in Column-2 (Herasymuik, 1996) _____________ 19 Figure 2.12 Hydraulic Conductivity Curve for Material in Column-1 (Herasymuik, 1996) ________________ 20 Figure 2.13 Total Discharge of the Fine Material for Column-1, after Newman (1999) _________________ 21 Figure 2.14 Total Discharge of the Coarse Material for Column-1, after Newman (1999) _______________ 21 Figure 2.15 Total Discharge for Column-2, after Newman (1999) __________________________________ 23 Figure 2.16 Model Results for EPM and TPROGS at T=65 Days A) Material Zone; B) Pressure Head; C) Saturation (Broda, et al., 2013) _____________________________________________________________ 26 Figure 2.17 Saturation Distribution in 5 Days (a) and Change in Time of Total Flow Entering the Base Compacted Layer (b) In X=35 M (Broda, et al., 2013) ____________________________________________ 27 Figure 2.18 General Characteristics of Soil Tank (Lv, et al., 2013) __________________________________ 29 Figure 2.19 Profile of Vector Flow at 100 Hr for a Tank Inclination of 9° and 19° (Lv, et al., 2013, p.
Recommended publications
  • Numerical Study on the Hydrologic Characteristic of Permeable Friction Course Pavement
    water Article Numerical Study on the Hydrologic Characteristic of Permeable Friction Course Pavement Tan Hung Nguyen 1 and Jaehun Ahn 2,* 1 Faculty of Architectural, Civil and Environmental Engineering, Nam Can Tho University, Can Tho 900000, Vietnam; [email protected] 2 Department of Civil and Environmental Engineering, Pusan National University, Busan 46241, Korea * Correspondence: [email protected]; Tel.: +82-51-510-7627 Abstract: The hydrologic characteristic of a permeable friction course (PFC) pavement is dependent on the rainfall intensity, pavement geometric design, and porous asphalt properties. Herein, the hydrologic characteristic of PFC pavements of various lengths and slopes was determined via numerical analysis. A series of analyses was conducted using length values of 10, 15, 20, and 30 m and slope values of 0.5%, 2%, 4%, 6%, and 8% for the equivalent water flow path. The PFC pavements were simulated for various values of rainfall intensity, which ranged from 10 to 120 mm/h, to determine the time taken for water to flow over the PFC pavement surface. The results show that the time for water overflow decreased when the pavement length or rainfall intensity increased, and it increased when the slope increased. Finally, a series of design charts was developed to determine the time taken for water to flow over the PFC pavement surface for given rainfall intensities. Since this study was conducted based on numerical analysis, further studies are recommended to verify experimentally the results presented. Citation: Nguyen, T.H.; Ahn, J. Keywords: hydrologic characteristic; permeable friction course pavement; geometric design Numerical Study on the Hydrologic Characteristic of Permeable Friction Course Pavement.
    [Show full text]
  • Sensitivity Study of Different Parameters Affecting Design of the Clay Blanket in Small Earthen Dams
    Journal of Himalayan Earth Sciences Volume 48, No. 2, 2015 pp.139-147 Sensitivity study of different parameters affecting design of the clay blanket in small earthen dams Ishtiaq Alam and Irshad Ahmad Department of Civil Engineering, University of Engineering & Technology, Peshawar, Pakistan. Abstract Dams are structures that retain water for human services. Dams may be earthen, concrete, timber, steel or masonry made. On the basis of size, they may be small, medium and large. The main purpose of a dam is to divert the flow of water for the intended use. Flow of water cannot be stopped permanently even by the best dam ever made. Water may seep from dam body, abutments or the foundation bed below the body of the dam. To control seepage from the foundation bed, certain available methods like cutoff trench, cutoff walls, diaphragms, grout curtains, sheet pile walls and upstream impervious blankets are used. Upstream impervious blankets are considered more economical compared with the other methods mentioned above. The key parameters playing role in blanket efficiency are length of blanket, thickness of blanket, clay core width of the dam, foundation bed depth up to impervious zone, reservoir head, permeability of blanket material and permeability of bed material. This study is focused on the effect of these parameters in seepage control. Seep/W, a finite element method based software is used to model all the mentioned parameters within the individually selected ranges. The results based on the software analysis show that when the length of blanket is gradually increased, the seepage quantity reduces gradually until a specific length where the effect of further increase in length become meaningless.
    [Show full text]
  • FHWA/TX-07/0-5202-1 Accession No
    Technical Report Documentation Page 1. Report No. 2. Government 3. Recipient’s Catalog No. FHWA/TX-07/0-5202-1 Accession No. 4. Title and Subtitle 5. Report Date Determination of Field Suction Values, Hydraulic Properties, August 2005; Revised March 2007 and Shear Strength in High PI Clays 6. Performing Organization Code 7. Author(s) 8. Performing Organization Report No. Jorge G. Zornberg, Jeffrey Kuhn, and Stephen Wright 0-5202-1 9. Performing Organization Name and Address 10. Work Unit No. (TRAIS) Center for Transportation Research 11. Contract or Grant No. The University of Texas at Austin 0-5202 3208 Red River, Suite 200 Austin, TX 78705-2650 12. Sponsoring Agency Name and Address 13. Type of Report and Period Covered Texas Department of Transportation Technical Report Research and Technology Implementation Office September 2004–August 2006 P.O. Box 5080 14. Sponsoring Agency Code Austin, TX 78763-5080 15. Supplementary Notes Project performed in cooperation with the Texas Department of Transportation and the Federal Highway Administration. Project Title: Determination of Field Suction Values in High PI Clays for Various Surface Conditions and Drain Installations 16. Abstract Moisture infiltration into highway embankments constructed by the Texas Department of Transportation (TxDOT) using high Plasticity Index (PI) clays results in changes in shear strength and in flow pattern that leads to recurrent slope failures. In addition, soil cracking over time increases the rate of moisture infiltration. The overall objective of this research is to determine the suction, hydraulic properties, and shear strength of high PI Texas clays. Specifically, two comprehensive experimental programs involving the characterization of unsaturated properties and the shear strength of a high PI clay (Eagle Ford clay) were conducted.
    [Show full text]
  • Water-Sciences Software Guide
    Table of Contents In The Name Of God the Compassionate the Merciful ! "#$% Application of Computers in Water-Sciences :&' ' ) #* Seyyed Javad Hoseiny : +, #./ +- Dr. Mhohamad Aflatuni 84 ) ' August 2005 219 3 -01 ! "#$% ,2/ " +, 34 5 6/ Table of Contents 1 Special Thanks 4 1 2 Abstract 5 2 3 Suggestions 6 3 4 Reason of Importance 7 4 5 Similar Researches 8 5 Detailed Guide for The First 6 9 #$% &' 12 !" # 6 Top 12 Software 7 Where to Find the Software 32 +,) # #$% &' () * 7 8 Software-Course List 33 /# .#$% &' -% 8 9 Rating Method 56 0# 1# 9 10 Initials & Expressions 58 23456 #!78 ,94 10 11 Full Software List 59 ,2/ " 7 6/ 11 12 CD Introduction 201 - %: 12 13 Website Introduction 202 - %: 13 To Those Who Are Whishing to # ;<4 0= >< 14 Continue Researching in This Field 204 14 * ? 15 Contact Us 205 / 15 16 Refrences and Resources 206 @8A B 16 17 What Do You Say? 218 + =C D 17 219 4 -01 ! "#$% ,2/ " +, 8 ' 9' Special Thanks ,' =8 =' #= #= , # # E ;= 'F=G H ? # ) &)D # 8 :;= - =: 3 # * ) (# 7-) '=G4% 7) . (* D N' *#) ' 7-) 2HL ,M' ":0 7) . 7D# # 7) =(' . 5< 5< O4 P=QH #@N' ) #= ,-< 1='% / . (;LS.;-# 7:6 N' R=L#=Q7# =7 7H' R=L# (; * # S 7 )) 'H3 * / . ( L7- Griffith N' 7-) Graham Jenkins 7) . (; /# N' 7-) ,: =0 7) . (VS - E.#- N' 7-) 3 #U T L8 7) . (S Utah State N' # S 7-) Wynn R. Walker 7) . ., # '#< ' Back to Contents 219 5 -01 ! "#$% ,2/ " +, :9-' 8; Abstract '# [7E #) VS &= ;XU36 ;-D#) Y;=(' ;) DS ? * Z .- E ? \ )S *=3 # =T= UG
    [Show full text]
  • Numerical Modeling in Geoenvironmental Practice
    Numerical Modeling in Geoenvironmental Practice By Craig H. Benson, Ph.D., P.E., M.ASCE o deling o f non-linear systems is now a regular Software Resources part of geoenvironmental engineering practice M due to Lh e ava il ability of inexpensive, high-speed A variety of powerful codes are available. The most PCs that utilize user-friendly software and graphical user common codes used for analysis o f vadose zone hydrol­ interfaces. Robust and effici ent solvers for no n-linear par­ ogy and unsaturated flow are 1-JYDRUS (www.pc-progress. lial d ifferential equations have also had a major irnpact com), UNSAT-H. (www.hydrology.pnl.gov or www.uwgeo on the practicality of geoenvironmental software packages. soft.org), VADOSE/W (www.geoslope.com), SEEP/W (www. Consequentl y, detailed analyses are being done in practice geoslope.corn), and SVFLUX (www.soilvision.com). that were not possible a decade ago. These include realistic HYDRUS is the most cornmonly used code in vadose zone assessrn ents of potential perfo rrn ance as well as "what if" hydrology worldwide, and can also be used to simulate scenarios. The rnost common analyses are associated with heat flow and contarninant transport in unsaturated rn edia. vadose-zone hydrology to predict the movement of water Other software packages cornrnonly used for contaminant in the vadose zone and i nteractions with the atmosphere. transpo 11 sirnulation in va riably saturated media include However, contarninant transport analyses in variably satu­ CTRAN/W (www.geoslope.com), CHEMFLUX (www.so il rated systerns are also becorning cornmonplace.
    [Show full text]
  • SVSLOPE® Support for Multi-Plane Analysis (MPA™)
    PRODUCT DATA SHEET SVSLOPE® Support for Multi-Plane Analysis (MPA™) SVSLOPE is the most advanced 3D slope stability analysis and slip surface search methods. The entire plane configuration process is designed so software available, with advanced searching methods that it is quick to perform on one or many planes at once. For example, the slope limits that are implemented to correctly determine the location may be defined for all planes at once by simply drawing a polygon that encloses the of the critical slip surface. The software provides you with area of interest on the 3D model. The slip surface search method is then automated, powerful 2D or 3D analysis for increased accuracy when with some options available to the user. calculating the factor of safety. Advanced probabilistic analysis or accommodation of spatial variation is possible Although there are multiple ways to create planes, the most common one, which with the software. SVSLOPE can be combined with was used in this case, is to simply select a point on each of the two banks. Planes SVFLUX™ to import pore-water pressures or SVSOLID™GT are then created along the slope automatically, at configurable distance intervals. to import soil stress conditions. The direction for each plane is automatically set based on the surface geometry. Each plane can be set to use multiple similar directions so that the critical direction Design and Analyze Different Locations Simultaneously is more likely to be found. Slope stability analysis is often targeted at topographically complex sites with features that vary greatly in three dimensions, or seemingly simple Results Collected and Aggregated into the Original 3D Model surface topology with strong and weak internal layers that vary across the for Visualization site.
    [Show full text]
  • Numerical Analysis of Leakage Through Geomembrane Lining Systems for Dams
    The First Pan American Geosynthetics Conference & Exhibition 2-5 March 2008, Cancun, Mexico Numerical Analysis of Leakage through Geomembrane Lining Systems for Dams C.T. Weber, University of Texas at Austin, Austin, TX, USA J.G. Zornberg, University of Texas at Austin, Austin, TX, USA ABSTRACT A numerical simulation was performed to characterize the effects of leakage through defects on the performance of dams with an upstream face lined with a geomembrane. The objective of this study was to determine how leakage through a lining system would affect the design of blanket toe drains in an earth dam. Toe drains decrease the pore pressure at the downstream face and keep the seepage line (or phreatic surface) below the downstream boundary. Simulations were conducted to determine the location of the phreatic surface in a homogeneous dam due to the presence of defects in the liner. Numerical simulations were also conducted to determine the length of the toe drain needed to prevent discharge from occurring on the downstream face of the dam. In addition, the effect of the elevation of the phreatic surface within the dam on the stability of the downstream face of the dam was analyzed. This study provides evidence on the benefits of using a geomembrane liner regarding the stability and toe drain in earth dams. 1. INTRODUCTION Geomembranes have been used as a solution to dam seepage problems since 1959, beginning in Europe (Sembenelli and Rodriguez 1996) and Canada (Lacroix 1984). These thin sheets of polymer have been used to make the upstream face watertight in roller-compacted concrete dams, to retrofit masonry and concrete dams, and as the main impervious layer in fill dams.
    [Show full text]
  • Dra-15. Hydrological Modeling
    DRA-15 DRA-15. HYDROLOGICAL MODELING Murray Fredlund, SoilVision Systems Ltd., May 19, 2008, revised December 23, 2008 Establish the flow system of the final conceptual model based on the GHN rock pile. 1. STATEMENT OF PROBLEM: What are the long-term pore-water pressures within the rock pile in both the upper and lower regions? What are worst-case scenarios for pore-water pressures given extreme climatic events in the next 100 years? Will there be ponding in any particular portion of the waste rock pile? What is the length of residence time of various types of pore-water in the system? What are the percent water saturation levels in the system? The influence of annual climate on the overall flow regime must be established. The distribution of pore-water pressures is needed for the slope stability modeling program. 2. PREVIOUS WORK: Previous related work in this area of the project includes the Phase I hydrological modeling activities. A summary of this work is shown in the following paragraphs. Comprehensive numerical modeling has been completed to provide insight into the possible/probable flow scenarios for the rock piles of the Questa mine. The 1D, 2D, and 3D numerical models all provided contributing pieces to the allow reasonable conclusions to be drawn. The most important conclusions from the hydrologic numerical modeling program are the following: 1D MODELING A large focus of the numerical modeling was to i) determine the sensitivity of various input parameters to the conceptual results as well as ii) to quantify the expected infiltration, and iii) the calibration of results to soil suctions and water contents measured near the surface by New Mexico Tech (NMT) and Golder Associates.
    [Show full text]
  • EVALUATING the RATE of LEAKAGE THROUGH DEFECTS in a GEOMEMBRANE a Thesis Submitted to the College of Graduate and Postdoctoral S
    EVALUATING THE RATE OF LEAKAGE THROUGH DEFECTS IN A GEOMEMBRANE A Thesis Submitted to the College of Graduate and Postdoctoral Studies In Partial Fulfillment of the Requirements For the Degree of Master of Science In the Department of Civil, Geological, and Environmental Engineering University of Saskatchewan Saskatoon By HALEY LOUISE CUNNINGHAM ® Copyright Haley Louise Cunningham, June 2018. All rights reserved. Permission to Use In presenting this thesis/dissertation in partial fulfillment of the requirements for a Postgraduate degree from the University of Saskatchewan, I agree that the Libraries of this University may make it freely available for inspection. I further agree that permission for copying of this thesis/dissertation in any manner, in whole or in part, for scholarly purposes may be granted by the professor or professors who supervised my thesis/dissertation work or, in their absence, by the Head of the Department or the Dean of the College in which my thesis work was done. It is understood that any copying or publication or use of this thesis/dissertation or parts thereof for financial gain shall not be allowed without my written permission. It is also understood that due recognition shall be given to me and to the University of Saskatchewan in any scholarly use which may be made of any material in my thesis/dissertation. Requests for permission to copy or to make other uses of materials in this thesis/dissertation in whole or part should be addressed to: Department of Civil, Geological, and Environmental Engineering 57 Campus Drive University of Saskatchewan Saskatoon, Saskatchewan S7N 5A9 Canada OR Dean College of Graduate and Postdoctoral Studies University of Saskatchewan 116 Thorvaldson Building, 110 Science Place Saskatoon, Saskatchewan S7N 5C9 Canada i Abstract Geomembranes are impermeable barriers when they remain intact.
    [Show full text]
  • Groundwater Modeling in Support of Water Resources Management and Planning Under Complex Climate, Regulatory, and Economic Stresses
    water Article Groundwater Modeling in Support of Water Resources Management and Planning under Complex Climate, Regulatory, and Economic Stresses Emin C. Dogrul *, Charles F. Brush and Tariq N. Kadir California Department of Water Resources, Bay-Delta Office, Room 252A, 1416 9th Street, Sacramento, CA 95814, USA; [email protected] (C.F.B.); [email protected] (T.N.K.) * Correspondence: [email protected]; Tel.: +1-916-654-7018 Academic Editor: M. Levent Kavvas Received: 31 October 2016; Accepted: 2 December 2016; Published: 13 December 2016 Abstract: Groundwater is an important resource that meets part or all of the water demand in many developed basins. Since it is an integral part of the hydrologic cycle, management of groundwater resources must consider not only the management of surface flows but also the variability in climate. In addition, agricultural and urban activities both affect the availability of water resources and are affected by it. Arguably, the Central Valley of the State of California, USA, can be considered a basin where all stresses that can possibly affect the management of groundwater resources seem to have come together: a vibrant economy that depends on water, a relatively dry climate, a disparity between water demand and availability both in time and space, heavily managed stream flows that are susceptible to water quality issues and sea level rise, degradation of aquifer conditions due to over-pumping, and degradation of the environment with multiple species becoming endangered. Over the past fifteen years, the California Department of Water Resources has developed and maintained the Integrated Water Flow Model (IWFM) to aid in groundwater management and planning under complex, and often competing, requirements.
    [Show full text]
  • Hydraulique Souterraine
    GROUNDWATER FLOW and SEEPAGE Engineering Department of Hydraulic Strutures & Environment Master of Hydraulics Academic Year 2012 – 2013 Course n°1: Introduction Dr. Robert WOUMENI FOREWORD Groundwater flow and Seepage problems are encountered in Civil and Environmental Engineering (i.e. flow through or under dams, toward wells or drains, around a sheetpile wall or a clay blanket,…); in Hydrology (i.e. infiltration of rain water in soils, water exchanges between embankment and rivers,…); and also in Hydrogeology (i.e. long term subsurface flows, ground water resource and quality, polluted soils rehabilitation, pumping,…). Robert WOUMENI 2 FOREWORD In all these problems, we have to deal with a fundamental principle which is the Darcy’s Law. The objective over an investigated area can be resumed as the following: If the soil hydrodynamic properties are known, we can calculate the groundwater discharge (Q), the head (H) and pore pressure (P), and then anticipate the occurrence of a breaking (Civil Engineering) or make the assessment of the water resource (Hydrogeology). If the spatial profile of groundwater is known for a given time (i.e. by means of field measurements) then the soil hydrodynamic properties can be estimated. Robert WOUMENI 3 FOREWORD So using Darcy’s Law, we would like through this course, to develop the following skills: Calculate the groundwater discharge, head and pressure, under various configurations (i.e. spatial geometry, anisotropy, saturated conditions, transient flow,…) with different methods graphical, analytical and numerical. Make the critical thinking of a numerical model (i.e. Hydronap, Modflow, Seep-2D); Estimate the hydrodynamic properties of soils (i.e.
    [Show full text]
  • PLAXIS® LE Essential LEM Slope Stability Analysis
    PRODUCT DATA SHEET PLAXIS® LE Essential LEM Slope Stability Analysis Bentley’s PLAXIS LE enables common to complex geotechnical analysis of soil or rock slopes. Large datasets from varied sources can be rapidly interpreted, prototyped, analyzed, and visualized. Perform comprehensive analysis of extensive sites, combined with spatial stability analysis in multiple concurrent locations. Choose the product that best fits your needs: • PLAXIS 2D LE (formerly 2D SVSLOPE, 2D SVFLUX, 2D SVSOLID, Complex open pit mine rapidly analyzed at multiple locations. SV SOILS, and SVSOILS Advanced) covers your 2D workflows. at each location. Consideration of faults, weak planes, and pore-water pressures, • PLAXIS 3D LE (formerly 2D/3D SVSLOPE, 2D/3D SVSLOPE Advanced, 2D/3D SVFLUX, 2D/3D SVSOLID, and 2D/3D along with the most available search methods on the market, provide confidence SVSOLID Advanced) supports your 3D workflows. in the design process. Search methods include Greco, Cuckoo, Wedges, and Fully-Specified for back analysis. Use probabilistic analysis—such as Monte Carlo, Latin Hypercube, and the Alternative Point Estimation Method (APEM)—to build Limit Equilibrium Slope Stability Methods with robust digital twins. Sensitivity analysis and spatial variability features offer Groundwater Analysis further model refinement. 3D analysis gives more rigorous consideration to site geology and increases accuracy when calculating safety factors. This process ensures that you are keeping Visualizing Digital Twins infrastructure safe and reliable. With PLAXIS LE, you can perform limit equilibrium Access state-of-the-art, report-ready graphical presentation of results without method (LEM) analysis choosing from over 15 classic method of slices and newer additional manipulation. The 3D immersive graphics engine provides performance stress-based methods.
    [Show full text]