Change search
ReferencesLink to record
Permanent link

Direct link
A High Order Finite Difference Method for Simulating Earthquake Sequences in a Poroelastic Medium
Uppsala University, Disciplinary Domain of Science and Technology, Mathematics and Computer Science, Department of Information Technology, Division of Scientific Computing.
Uppsala University, Disciplinary Domain of Science and Technology, Mathematics and Computer Science, Department of Information Technology, Division of Scientific Computing.
2016 (English)Independent thesis Advanced level (professional degree), 20 credits / 30 HE creditsStudent thesis
Abstract [en]

Induced seismicity (earthquakes caused by injection or extraction of fluids in Earth's subsurface) is a major, new hazard in the United States, the Netherlands, and other countries, with vast economic consequences if not properly managed. Addressing this problem requires development of predictive simulations of how fluid-saturated solids containing frictional faults respond to fluid injection/extraction. Here we present a numerical method for linear poroelasticity with rate-and-state friction faults. A numerical method for approximating the fully coupled linear poroelastic equations is derived using the summation-by-parts-simultaneous-approximation-term (SBP-SAT) framework. Well-posedness is shown for a set of physical boundary conditions in 1D and in 2D. The SBP-SAT technique is used to discretize the governing equations and show semi-discrete stability and the correctness of the implementation is verified by rigorous convergence tests using the method of manufactured solutions, which shows that the expected convergence rates are obtained for a problem with spatially variable material parameters. Mandel's problem and a line source problem are studied, where simulation results and convergence studies show satisfactory numerical properties. Furthermore, two problem setups involving fault dynamics and slip on faults triggered by fluid injection are studied, where the simulation results show that fluid injection can trigger earthquakes, having implications for induced seismicity. In addition, the results show that the scheme used for solving the fully coupled problem, captures dynamics that would not be seen in an uncoupled model. Future improvements involve imposing Dirichlet boundary conditions using a different technique, extending the scheme to handle curvilinear coordinates and three spatial dimensions, as well as improving the high-performance code and extending the study of the fault dynamics. 

Place, publisher, year, edition, pages
UPTEC F, ISSN 1401-5757 ; 16038
Keyword [en]
numerical analysis, numerical methods, numerical modeling, SBP-SAT, finite differences, high performance computing, geophysics, geomechanics, poroelasticity, fault mechanics, induced seismicity
National Category
Geophysics Physical Sciences
URN: urn:nbn:se:uu:diva-298414OAI: diva2:945945
External cooperation
Stanford University
Educational program
Master Programme in Engineering Physics
Available from: 2016-07-04 Created: 2016-07-04 Last updated: 2016-07-04Bibliographically approved

Open Access in DiVA

fulltext(2008 kB)48 downloads
File information
File name FULLTEXT01.pdfFile size 2008 kBChecksum SHA-512
Type fulltextMimetype application/pdf

By organisation
Division of Scientific Computing
GeophysicsPhysical Sciences

Search outside of DiVA

GoogleGoogle Scholar
Total: 48 downloads
The number of downloads is the sum of all downloads of full texts. It may include eg previous versions that are now no longer available

Total: 250 hits
ReferencesLink to record
Permanent link

Direct link