{"id":{"repo_id":"mit","oai_identifier":"oai:dspace.mit.edu:1721.1/121757"},"canonical_url":"https://search.dev.ndltd.org/etd/mit/oai:dspace.mit.edu:1721.1/121757","repository":{"repo_id":"mit","name":"MIT","base_url":"https://dspace.mit.edu/oai/request"},"display":{"title":"Applications of a fast helmholtz solver in exploration seismology","abstract":"Seismic imaging techniques rely on a velocity model inverted from noisy data via a non-linear inverse problem. This inferred velocity model may be inaccurate and lead to incorrect interpretations of the subsurface. In this thesis, I combine a fast Helmholtz solver, the field expansion method, with a reduced velocity model parameterization to address the impact of an uncertain or inaccurate velocity model. I modify the field expansion framework to accurately simulate the acoustic field for velocity models that commonly occur in seismic imaging. The field expansion method describes the acoustic field in a periodic medium in which the velocity model and source repeat infinitely in the horizontal direction, much like a diffraction grating. This Helmholtz solver achieves significant computational speed by restricting the velocity model to consists of a number of non-overlapping piecewise layers.","abstract_html":"Seismic imaging techniques rely on a velocity model inverted from noisy data via a non-linear inverse problem. This inferred velocity model may be inaccurate and lead to incorrect interpretations of the subsurface. In this thesis, I combine a fast Helmholtz solver, the field expansion method, with a reduced velocity model parameterization to address the impact of an uncertain or inaccurate velocity model. I modify the field expansion framework to accurately simulate the acoustic field for velocity models that commonly occur in seismic imaging. The field expansion method describes the acoustic field in a periodic medium in which the velocity model and source repeat infinitely in the horizontal direction, much like a diffraction grating. This Helmholtz solver achieves significant computational speed by restricting the velocity model to consists of a number of non-overlapping piecewise layers.","abstract_has_math":false,"creators":["Ely, Gregory Tsiang."],"institution":"Massachusetts Institute of Technology","degree_name":"Doctoral","degree_level":null,"degree_discipline":null,"degree_department":"Massachusetts Institute of Technology. Department of Earth, Atmospheric, and Planetary Sciences","school":null,"contributors":[],"advisors":["Alison Malcolm."],"committee_chairs":[],"committee_members":[],"year":2019,"date_issued":"2019","date_published":"2019","updated_at":"2026-07-22T22:21:40Z","subjects":["Earth, Atmospheric, and Planetary Sciences."],"languages":["eng"],"rights":["MIT theses are protected by copyright. They may be viewed, downloaded, or printed from this source but further reproduction or distribution in any format is prohibited without written permission."],"rights_urls":["http://dspace.mit.edu/handle/1721.1/7582"],"identifier_entries":[]},"links":{"outbound_url":"https://hdl.handle.net/1721.1/121757","outbound_label":"Handle","outbound_source":"dc:identifier.uri"},"metadata_groups":[{"id":"people","label":"People","entries":[{"key":"dc:contributor.advisor","label":"Advisor","values":["Alison Malcolm."]},{"key":"dc:contributor.department","label":"Department","values":["Massachusetts Institute of Technology. Department of Earth, Atmospheric, and Planetary Sciences","EAPS"]},{"key":"dc:contributor.other","label":"Dc Contributor Other","values":["Massachusetts Institute of Technology. Department of Earth, Atmospheric, and Planetary Sciences."]},{"key":"dc:creator","label":"Author","values":["Ely, Gregory Tsiang."]}]},{"id":"academic_context","label":"Academic Context","entries":[{"key":"dc:date.accessioned","label":"Dc Date Accessioned","values":["2019-07-17T21:01:25Z"]},{"key":"dc:date.available","label":"Dc Date Available","values":["2019-07-17T21:01:25Z"]},{"key":"dc:date.issued","label":"Date","values":["2019"]},{"key":"dc:publisher","label":"Institution","values":["Massachusetts Institute of Technology"]},{"key":"dc:type","label":"Dc Type","values":["Thesis"]},{"key":"thesis:degree_name","label":"Degree Name","values":["Doctoral"]}]},{"id":"subjects_keywords","label":"Subjects and Keywords","entries":[{"key":"dc:subject","label":"Dc Subject","values":["Earth, Atmospheric, and Planetary Sciences."]}]},{"id":"language_rights","label":"Language and Rights","entries":[{"key":"dc:language.iso","label":"Language (ISO)","values":["eng"]},{"key":"dc:rights","label":"Dc Rights","values":["MIT theses are protected by copyright. They may be viewed, downloaded, or printed from this source but further reproduction or distribution in any format is prohibited without written permission."]},{"key":"dc:rights.uri","label":"Rights URI","values":["http://dspace.mit.edu/handle/1721.1/7582"]}]},{"id":"identifiers","label":"Identifiers","entries":[{"key":"dc:identifier.uri","label":"Identifier URI","values":["https://hdl.handle.net/1721.1/121757"]}]},{"id":"additional","label":"Additional Metadata","entries":[{"key":"dc:description","label":"Description","values":["Thesis: Ph. D., Massachusetts Institute of Technology, Department of Earth, Atmospheric, and Planetary Sciences, 2019","Cataloged from PDF version of thesis.","Includes bibliographical references (pages 141-150)."]},{"key":"dc:description.abstract","label":"Abstract","values":["Seismic imaging techniques rely on a velocity model inverted from noisy data via a non-linear inverse problem. This inferred velocity model may be inaccurate and lead to incorrect interpretations of the subsurface. In this thesis, I combine a fast Helmholtz solver, the field expansion method, with a reduced velocity model parameterization to address the impact of an uncertain or inaccurate velocity model. I modify the field expansion framework to accurately simulate the acoustic field for velocity models that commonly occur in seismic imaging. The field expansion method describes the acoustic field in a periodic medium in which the velocity model and source repeat infinitely in the horizontal direction, much like a diffraction grating. This Helmholtz solver achieves significant computational speed by restricting the velocity model to consists of a number of non-overlapping piecewise layers.","I modify this restricted framework to allow for the modeling of more complex velocity models with dozens of parameters instead of the thousands or millions of parameters used to characterize pixelized velocity models. This parameterization, combined with the speed of the forward solver allow me to examine two problems in seismic imaging: uncertainty quantification and benchmarking global optimization methods. With the rapid speed of the forward solver, I use Markov Chain Monte Carlo methods to estimate the non-linear probability distribution of a 2D seismic velocity model given noisy data. Although global optimization methods have recently been applied to inversion of seismic velocity model using raw waveform data, it has been impossible to compare various types of algorithms and impacts of parameters on convergence. The reduced forward model presented in this paper allows me to benchmark these algorithms and objectively compare their performance to one another.","I also explore the application of these and other geophysical methods to a medical ultrasound dataset that is well approximated by a layered model."]},{"key":"dc:description.degree","label":"Dc Description Degree","values":["Ph. D."]},{"key":"dc:title","label":"Title","values":["Applications of a fast helmholtz solver in exploration seismology"]}]}],"canonical_facts":{"dc:contributor.advisor":["Alison Malcolm."],"dc:contributor.department":["Massachusetts Institute of Technology. Department of Earth, Atmospheric, and Planetary Sciences","EAPS"],"dc:contributor.other":["Massachusetts Institute of Technology. Department of Earth, Atmospheric, and Planetary Sciences."],"dc:creator":["Ely, Gregory Tsiang."],"dc:date.accessioned":["2019-07-17T21:01:25Z"],"dc:date.available":["2019-07-17T21:01:25Z"],"dc:date.issued":["2019"],"dc:description":["Thesis: Ph. D., Massachusetts Institute of Technology, Department of Earth, Atmospheric, and Planetary Sciences, 2019","Cataloged from PDF version of thesis.","Includes bibliographical references (pages 141-150)."],"dc:description.abstract":["Seismic imaging techniques rely on a velocity model inverted from noisy data via a non-linear inverse problem. This inferred velocity model may be inaccurate and lead to incorrect interpretations of the subsurface. In this thesis, I combine a fast Helmholtz solver, the field expansion method, with a reduced velocity model parameterization to address the impact of an uncertain or inaccurate velocity model. I modify the field expansion framework to accurately simulate the acoustic field for velocity models that commonly occur in seismic imaging. The field expansion method describes the acoustic field in a periodic medium in which the velocity model and source repeat infinitely in the horizontal direction, much like a diffraction grating. This Helmholtz solver achieves significant computational speed by restricting the velocity model to consists of a number of non-overlapping piecewise layers.","I modify this restricted framework to allow for the modeling of more complex velocity models with dozens of parameters instead of the thousands or millions of parameters used to characterize pixelized velocity models. This parameterization, combined with the speed of the forward solver allow me to examine two problems in seismic imaging: uncertainty quantification and benchmarking global optimization methods. With the rapid speed of the forward solver, I use Markov Chain Monte Carlo methods to estimate the non-linear probability distribution of a 2D seismic velocity model given noisy data. Although global optimization methods have recently been applied to inversion of seismic velocity model using raw waveform data, it has been impossible to compare various types of algorithms and impacts of parameters on convergence. The reduced forward model presented in this paper allows me to benchmark these algorithms and objectively compare their performance to one another.","I also explore the application of these and other geophysical methods to a medical ultrasound dataset that is well approximated by a layered model."],"dc:description.degree":["Ph. D."],"dc:identifier.uri":["https://hdl.handle.net/1721.1/121757"],"dc:language.iso":["eng"],"dc:publisher":["Massachusetts Institute of Technology"],"dc:rights":["MIT theses are protected by copyright. They may be viewed, downloaded, or printed from this source but further reproduction or distribution in any format is prohibited without written permission."],"dc:rights.uri":["http://dspace.mit.edu/handle/1721.1/7582"],"dc:subject":["Earth, Atmospheric, and Planetary Sciences."],"dc:title":["Applications of a fast helmholtz solver in exploration seismology"],"dc:type":["Thesis"],"thesis:degree_name":["Doctoral"]},"updated_at":"2026-07-22T22:21:40Z"}