Exploring the sparsity of seismic wave propagation in complex media

Owing to continual improvements in computation power, we are getting closer to modelling seismic waves in ways that reflect the realistic complexity of the Earth. Nevertheless, modelling the true multi-scale structure of the Earth at large scale and/or high frequencies is currently still beyond our...

Full description

Bibliographic Details
Main Author: Haindl, CM
Other Authors: Nissen-Meyer, T
Format: Thesis
Language:English
Published: 2021
Subjects:
_version_ 1797106995425181696
author Haindl, CM
author2 Nissen-Meyer, T
author_facet Nissen-Meyer, T
Haindl, CM
author_sort Haindl, CM
collection OXFORD
description Owing to continual improvements in computation power, we are getting closer to modelling seismic waves in ways that reflect the realistic complexity of the Earth. Nevertheless, modelling the true multi-scale structure of the Earth at large scale and/or high frequencies is currently still beyond our reach from a computational cost perspective. Hence, we rely on efficient modelling software to increase the scope of seismic simulations. Azimuthal Complexity Adaptation (ACA) is an innovative way of reducing the computational cost of seismic simulations by exploiting inherent sparsity of the wavefield. In this thesis, AxiSEM3D, an ACA-based modelling solver, is expanded to model wave propagation in local-scale settings. I find that the sparsity of wavefields persists in the presence of complex local-scale structures such as thrust faults and salt bodies and in models with full anisotropy. Further, I develop 3D wavefield scanning to locate complex areas of the wavefield and to investigate how they relate to structures in the subsurface model. The results suggest that high wavefield complexity is typically localised near abrupt wavespeed contrasts, on the low-wavespeed side thereof. I attempt to develop a new modelling solver which overcomes some of the limitations of AxiSEM3D. The resulting method is unstable, but it demonstrates the algorithmic complexity required to perform ACA in settings like coast lines. I conclude that a hybrid method which combines classic 3D modelling and ACA may be a simpler solution for increasing the scope of AxiSEM3D. Finally, I implement fault ruptures to simulate an earthquake scenario in the San Francisco Bay Area. Using ACA, I can include low-wavespeed bay muds which are usually neglected to save computational cost. Neglecting these muds leads to an underestimate of the shaking in the San Francisco Bay, even at low frequencies.
first_indexed 2024-03-07T07:10:22Z
format Thesis
id oxford-uuid:283e030f-e320-46f5-bd83-9e01762b1442
institution University of Oxford
language English
last_indexed 2024-03-07T07:10:22Z
publishDate 2021
record_format dspace
spelling oxford-uuid:283e030f-e320-46f5-bd83-9e01762b14422022-06-20T09:16:16ZExploring the sparsity of seismic wave propagation in complex mediaThesishttp://purl.org/coar/resource_type/c_db06uuid:283e030f-e320-46f5-bd83-9e01762b1442Earthquake hazard analysisSeismic prospectingEquations--Numerical solutionsCalifornia--Hayward FaultSeismic wave propagationEnglishHyrax Deposit2021Haindl, CMNissen-Meyer, THawthorne, JPeter, DOwing to continual improvements in computation power, we are getting closer to modelling seismic waves in ways that reflect the realistic complexity of the Earth. Nevertheless, modelling the true multi-scale structure of the Earth at large scale and/or high frequencies is currently still beyond our reach from a computational cost perspective. Hence, we rely on efficient modelling software to increase the scope of seismic simulations. Azimuthal Complexity Adaptation (ACA) is an innovative way of reducing the computational cost of seismic simulations by exploiting inherent sparsity of the wavefield. In this thesis, AxiSEM3D, an ACA-based modelling solver, is expanded to model wave propagation in local-scale settings. I find that the sparsity of wavefields persists in the presence of complex local-scale structures such as thrust faults and salt bodies and in models with full anisotropy. Further, I develop 3D wavefield scanning to locate complex areas of the wavefield and to investigate how they relate to structures in the subsurface model. The results suggest that high wavefield complexity is typically localised near abrupt wavespeed contrasts, on the low-wavespeed side thereof. I attempt to develop a new modelling solver which overcomes some of the limitations of AxiSEM3D. The resulting method is unstable, but it demonstrates the algorithmic complexity required to perform ACA in settings like coast lines. I conclude that a hybrid method which combines classic 3D modelling and ACA may be a simpler solution for increasing the scope of AxiSEM3D. Finally, I implement fault ruptures to simulate an earthquake scenario in the San Francisco Bay Area. Using ACA, I can include low-wavespeed bay muds which are usually neglected to save computational cost. Neglecting these muds leads to an underestimate of the shaking in the San Francisco Bay, even at low frequencies.
spellingShingle Earthquake hazard analysis
Seismic prospecting
Equations--Numerical solutions
California--Hayward Fault
Seismic wave propagation
Haindl, CM
Exploring the sparsity of seismic wave propagation in complex media
title Exploring the sparsity of seismic wave propagation in complex media
title_full Exploring the sparsity of seismic wave propagation in complex media
title_fullStr Exploring the sparsity of seismic wave propagation in complex media
title_full_unstemmed Exploring the sparsity of seismic wave propagation in complex media
title_short Exploring the sparsity of seismic wave propagation in complex media
title_sort exploring the sparsity of seismic wave propagation in complex media
topic Earthquake hazard analysis
Seismic prospecting
Equations--Numerical solutions
California--Hayward Fault
Seismic wave propagation
work_keys_str_mv AT haindlcm exploringthesparsityofseismicwavepropagationincomplexmedia