Numerical Evaluation and Analysis of the Occurrence of Earth Fissures in Faulted Sedimentary Basins
Total Page:16
File Type:pdf, Size:1020Kb
Numerical Evaluation and Analysis of the Occurrence of Earth Fissures in Faulted Sedimentary Basins Martin Hernandez-Marin Dissertation submitted to the faculty of the Virginia Polytechnic Institute and State University in partial fulfillment of requirements for the degree of Doctor of Philosophy in Geosciences Committee Thomas J. Burbey, Chair Madeline E. Schreiber Mark A. Widdowson Daniel L. Gallagher November 16, 2009 Blacksburg, VA Keywords: Earth fissuring, land subsidence, finite element method, ABAQUS Copyright © Martin Hernandez-Marin Numerical Evaluation and Analysis of the Occurrence of Earth Fissures in Faulted Sedimentary Basins Martin Hernandez-Marin ABSTRACT This dissertation describes the occurrence of pumping-induced earth fissures associated with quaternary faulting using numerical simulations. The Eglington Fault located in Las Vegas valley has been selected as the prototype fault described herein. The finite-element software program ABAQUS is used for the numerical simulations. The Eglington fault area is chosen because it represents one of the best examples displaying the complex relationship between fissuring, faulting and pumping-induced stress. This fault is known to influence both the vertical and horizontal deformation patterns through the accumulation of stress in its vicinity. The result is that fissures are observed on both sides of the fault and in close proximity to the fault plane. In addition to the complex fault-fissure connection, a thick caliche-rich vadose zone with weak mechanical strength allows for the initiation and propagation of fissures. The numerical analysis a) investigates the geometrical and hydromechanical features of the zone of influence surrounding the Eglington Fault; b) identifies the zones of accumulated stress on the surface and at depth that can lead to fissuring; and c) simulates the onset and propagation of tensile-induced fissures. Three-dimensional numerical simulations of this fault indicate that a 100-meter wide fault-zone composed by sand-like material best reproduces the conditions of stress that may lead to fissuring in the vicinity of the fault. Additionally, two-dimensional models reveal that two main mechanisms promote the accumulation of stress in the vicinity of the fault zone: one is the counterclockwise rotation of the unsaturated portion of the fault zone; the other is the differential compaction caused by the difference in the accumulated thickness of compressible layers. Tensile stress is concentrated on the surface in the hanging wall, but maximum shear stress zones are simulated to occur on both sides of the fault at the contact between the saturated aquifer and the vadose zone. A final analysis of the initiation and propagation of tensile-induced fissures demonstrates that fissures commence and propagate only within the vadose zone, and that the propagation path is influenced by the mechanical properties of the medium and the location of the main load, which in this case is pumping. iii DEDICATION This dissertation and all work behind it is dedicated to two exceptional persons: My dear Norma, and loved Mother. Mother, I am sure you are always taking care of me in whatever part of heaven you are now. Norma, thanks for every little thing you do for me and for every little moment you spend with me. iv ACKNOWLEDGEMENTS I want to thank my friend and advisor Thomas Burbey. Tom, thanks for accepting me as your student in the first place and thanks for guiding me toward the completion of this degree. I owe you a lot! Thanks for everything. I also thank Dr. Madeline Schreiber for her encouraging advice. Thanks, Maddy, for all your help! I also appreciate the valuable guidance and advice of my committee members Dr. Mark Widdowson and Dr. Dan Gallagher. I thank my excellent first friends of my office in Derring 3051: Isaac D. Jeng, who gave me an excellent welcome in Blacskburg along with his wonderful wife Audrey; Rachel Lauer who taught me the basics on “How to survive in the USA” along with David Rugh; I also thank my good friend Tingting from China. I will never forget the moments in class and discussion from homeworks that I spent with them, particularly with Rachel and Tingting. I thank my other old friends that left our department earlier than me: JP, Ben, Ankan, Niki and Brad, and those who continue making research: Yinka, Jeanne, Meijing, Youquan, Beth, Jonathan and Daniel. I received enormous help and friendship from them. I am very grateful to my family, specially to my lovely father and Mother. I owe them a big part of this PhD. Dear brothers, sisters and parents, thanks for all your unconditional support in whatever I decide to do in my life. The essential help on the administrative issues at Virginia Tech was provided by Connie Lowe, Mary McMurray, Linda Bland, and Carolyn Williams. I thank and congratulate all these excellent persons for doing their job very efficiently. I am especially grateful to the research group of Dr. Romesh Batra, particularly to Alireza Chadegani for spending some time giving me some crucial advice on ABAQUS. Finally, I recognize the economic support from CONACyT (National Council of Science and Technology) and our department, during most of my PhD. Also I recognize the economic help I received from the Graduate School at the end of my PhD, and I thank Dean Karen DePauw for this. v ATTRIBUTIONS This dissertation consists of three journal papers. Chapter two was published in the Hydrogeology Journal (DOI 10.1007/s10040-009-0501-8), “Hernandez-Marin, M., and Burbey. T.J., The role of faulting on surface deformation patterns from pumping-induced groundwater flow (Las Vegas Valley, USA). The concept of the project, model simulations and interpretations, writing and figure preparations were performed by M. Hernandez-Marin. T.J. Burbey clarified text and reviewed the manuscript. Chapter three was submitted for publication to the Journal of Hydrology: “Hernandez- Marin, M., and Burbey. T.J., Fault-controlled deformation and stress from pumping-induced ground-water flow in a hydrogeologic setting. The concept of the project, model simulations and interpretations, writing and figure preparations were performed by M. Hernandez-Marin. T.J. Burbey clarified text and reviewed the manuscript. Chapter four is being prepared for submission to the Hydrogeology Journal. “Hernandez- Marin, M., and Burbey. T.J., Controls on initiation and propagation of pumping-induced earth fissures: insights from numerical simulations”. The concept of the project, model simulations and interpretations, writing, and figure preparations were performed by M. Hernandez-Marin. T.J. Burbey clarified text, reviewed the manuscript and contributed with the main conclusion. vi Table of Contents Abstract……………………………………………………………………………I Dedication…………………………………………………………………………III Acknowledgments………………………………………………………………... IV Attributions……………………………………………………………………….. V Table of contents………………………………………………………………….. VI List of figures……………………………………………………………………...VIII List of tables……………………………………………………………………….X Chapter 1. Introduction…………………………………………………………… 1 Chapter 2. The role of faulting on surface deformation patterns from pumping- induced groundwater flow (Las Vegas Valley, USA)…………………………….. 3 Abstract……………………………………………………………………………4 2.1. Introduction………………………………………………………………….5 2.2. Las Vegas Valley as a prototype model……………………………………...7 2.3. Numerical modeling…………………………………………………………10 2.3.1. Conceptual models using ABAQUS……………………………………..10 2.3.2. Loads, boundary conditions and time stepping…………………………..11 2.4. Modeling results and discussion……………………………………………. 16 2.4.1. Distribution of vertical displacements at the land surface………………...16 2.4.2. Distribution of magnitudes of horizontal displacements at the land surface……………………………………………………………………20 2.4.3. Direction and evolution of horizontal displacements at the land surface.. 23 2.4.4. Surface zones under compression and tension…………………………...28 2.4.5. Model limitations and future work……………………………………… 32 2.5. Summary and conclusions………………………………………………….. 32 References………………......……………………………………………………34 Chapter 3. Fault-controlled deformation and stress from pumping-induced ground-water flow in a hydrogeologic setting…………………………………….37 Abstract……………………………………………………………………………38 3.1 Introduction…………………………………………………………………..39 3.2. Elasto-Plastic mechanical behavior………………………………………… 41 3.2.1 Poroelasticity theory……………………………………………………... 41 3.2.2 Elasto-plastic mechanical response……………………………………….44 3.3. Model Development…………………………………………………………48 3.4. Model scenarios…………………………………………………………….. 51 3.4.1. Scenarios I and II: Purely-elastic hydrostratigraphic layers………...…… 51 3.4.2. Scenarios III and VI: Elasto-plastic aquitards………………………...…. 51 3.5. Relevance of a cyclical groundwater pumping rate…………………………52 3.6. Simulated Vertical Displacement at the Land Surface………………………55 3.7. Simulated horizontal displacements………………..……………………… 56 3.7.1 Horizontal displacements on the surface and at the saturated/unsaturated interface…………………………………………………………………. 57 3.7.2 Vertical distribution of horizontal displacements in the vicinity of the fault……………………………………………………………………… 59 vii 3.8. Simulated stress…………………………………………………………….. 62 3.8.1 Distribution of tensile and compressive stress……………………………62 3.8.2 Distribution of shear stress………………………………………………..64 3.9. Discussion………………………………………………………………...….65 3.10. Conclusions…………………………………………………………………67 References………………….….……………….......………………………..……69 Chapter 4. Controls on Initiation and Propagation of pumping-induced