Seismic Ray Impedance Inversion

This thesis investigates a prestack seismic inversion scheme implemented in the ray parameter domain. Conventionally, most prestack seismic inversion methods are performed in the incidence angle domain. However, inversion using the concept of ray impedance, as it honours ray path variation following the elastic parameter variation according to Snell’s law, shows the capacity to discriminate different lithologies if compared to conventional elastic impedance inversion. The procedure starts with data transformation into the ray-parameter domain and then implements the ray impedance inversion along constant ray-parameter profiles. With different constant-ray-parameter profiles, mixed-phase wavelets are initially estimated based on the high-order statistics of the data and further refined after a proper well-to-seismic tie. With the estimated wavelets ready, a Cauchy inversion method is used to invert for seismic reflectivity sequences, aiming at recovering seismic reflectivity sequences for blocky impedance inversion. The impedance inversion from reflectivity sequences adopts a standard generalised linear inversion scheme, whose results are utilised to identify rock properties and facilitate quantitative interpretation. It has also been demonstrated that we can further invert elastic parameters from ray impedance values, without eliminating an extra density term or introducing a Gardner’s relation to absorb this term. Ray impedance inversion is extended to P-S converted waves by introducing the definition of converted-wave ray impedance. This quantity shows some advantages in connecting prestack converted wave data with well logs, if compared with the shearwave elastic impedance derived from the Aki and Richards approximation to the Zoeppritz equations. An analysis of P-P and P-S wave data under the framework of ray impedance is conducted through a real multicomponent dataset, which can reduce the uncertainty in lithology identification.Inversion is the key method in generating those examples throughout the entire thesis as we believe it can render robust solutions to geophysical problems. Apart from the reflectivity sequence, ray impedance and elastic parameter inversion mentioned above, inversion methods are also adopted in transforming the prestack data from the offset domain to the ray-parameter domain, mixed-phase wavelet estimation, as well as the registration of P-P and P-S waves for the joint analysis. The ray impedance inversion methods are successfully applied to different types of datasets. In each individual step to achieving the ray impedance inversion, advantages, disadvantages as well as limitations of the algorithms adopted are detailed. As a conclusion, the ray impedance related analyses demonstrated in this thesis are highly competent compared with the classical elastic impedance methods and the author would like to recommend it for a wider application

Advancements in seismic tomography with application to tunnel detection and volcano imaging

Thesis (Ph.D.) University of Alaska Fairbanks, 1998Practical geotomography is an inverse problem with no unique solution. A priori information must be imposed for a stable solution to exist. Commonly used types of a priori information smooth and attenuate anomalies, resulting in 'blurred' tomographic images. Small or discrete anomalies, such as tunnels, magma conduits, or buried channels are extremely difficult imaging objectives. Composite distribution inversion (CDI) is introduced as a theory seeking physically simple, rather than distributionally simple, solutions of non-unique problems. Parameters are assumed to be members of a composite population, including both well-known and anomalous components. Discrete and large amplitude anomalies are allowed, while a well-conditioned inverse is maintained. Tunnel detection is demonstrated using CDI tomography and data collected near the northern border of South Korea. Accurate source and receiver location information is necessary. Borehole deviation corrections are estimated by minimizing the difference between empirical distributions of apparent parameter values as a function of location correction. Improved images result. Traveltime computation and raytracing are the most computationally intensive components of seismic tomography when imaging structurally complex media. Efficient, accurate, and robust raytracing is possible by first recovering approximate raypaths from traveltime fields, and then refining the raypaths to a desired accuracy level. Dynamically binned queuing is introduced. The approach optimizes graph-theoretic traveltime computation costs. Pseudo-bending is modified to efficiently refine raypaths in general media. Hypocentral location density functions and relative phase arrival population analysis are used to investigate the Spring, 1996, earthquake swarm at Akutan Volcano, Alaska. The main swarm is postulated to have been associated with a 0.2 km\sp3 intrusion at a depth of less than four kilometers. Decay sequence seismicity is postulated to be a passive response to the stress transient caused by the intrusion. Tomograms are computed for Mt. Spurr, Augustine, and Redoubt Volcanoes, Alaska. Relatively large amplitude, shallow anomalies explain most of the traveltime residual. No large amplitude anomalies are found at depth, and no magma storage areas are imaged. A large amplitude low-velocity anomaly is coincident with a previously proposed geothermal region on the southeast flank of Mt. Spurr. Mt. St. Augustine is found to have a high velocity core

Seismic and magnetic constraints on the structure of upper oceanic crust at fast and slow spreading ridges

Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 1998The upper ocean crust contains a comprehensive record of the shallow geological processes active along the world's mid-ocean ridge system. This thesis examines the magnetic and seismic structure of the upper crust at two contrasting ridges-the fast spreading East Pacific Rise (EPR) and the slow spreading Mid-Atlantic Ridge (MAR)-to build a more complete understanding about the roles of volcanic emplacement, tectonic disruption and hydrothermal alteration in the near-ridge environment. A technique that inverts potential field measurements. directly from an uneven observation track is developed and applied to near-bottom magnetic data from the spreading segments north of the Kane transform on the MAR. It is concluded that the central anomaly magnetization high marks the locus of focused volcanic emplacement. A cyclic faulting model is proposed to explain the oscillatory magnetization pattern associated with discrete blocks of crust being transported out of the rift valley between intensely altered fault zones. Seismic waveform and amplitude analyses of the magma sill along the EPR reveal it to be a thin (<100 m) body of partial melt. These characteristics have important implications for melt availability and transport within the cycle of eruption and replenishment. A genetic algorithm-based seismic waveform inversion technique is developed and applied to on- and near-axis multichannel data from 17°20'S on the EPR and the spreading segment south of the Oceanographer transform (MAR) to map and compare for the first time the detailed velocity structure of the upper crust at two different spreading rates. Combined with conventionally processed seismic profiles, our results show that, while final extrusive thickness is comparable at all spreading ridges (300-500 m), the style of thickening may vary. While a thin (≤100 m) extrusive carapace quadruples in thickness within 1-4 km of the EPR crest, the extrusive section at the MAR achieves its final thickness within the inner valley. Both show evidence for a narrow zone of volcanic emplacement. Vigorous hydrothermalism at the EPR may produce a more rapid increase in basement velocities relative to the MAR. Rapid modification of the extrusive/dike transition at both ridges indicates that hydrothermalism is enhanced in this interval. Along-axis transport of lavas may thicken the extrusive pile at slow spreading segment ends, strengthening the magnetic highs generated by lava chemistry.ONR graduate fellow. The magnetics portion of this thesis was also supported by NSF grants OCE-9204141, OCE-9200905 (M. A. Tivey & H. Schouten) and NERC GR3/7702 (R. C. Searle). Research for Chapter 4 was partially funded by NSF grant OCE-9402933 (R. S. Detrick). The final two science chapters were supported by NSF grants OCE-9012707, OCE-9300450 (R. S. Detrick), OCE-9401717 (G. M. Kent & R. S. Detrick), OCE-9400623 (M. A. Tivey) and the Education office

Advanced geophysical studies of accretion of oceanic lithosphere in Mid-Ocean Ridges characterized by contrasting tectono-magmatic settings

Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution February 2012The structure of the oceanic lithosphere results from magmatic and extensional processes taking place at mid-ocean ridges (MORs). The temporal and spatial scales of the variability of these two processes control the degree of heterogeneity of the oceanic lithosphere, represented by two end-member models: the classical Penrose Model exemplified by layered magmatic crust formed along fast-spreading MORs, e.g., East Pacific Rise (EPR); and the recently defined Chapman Model describing heterogeneous mafic and ultramafic lithosphere formed in settings of oceanic detachment faulting common along slow-spreading MORs, e.g., Mid-Atlantic Ridge (MAR). This thesis is using advanced marine geophysical methods (including finite-difference wave propagation modeling, 3D multi-channel seismic reflection imaging, waveform inversion, streamer tomography, and near-bottom magnetics) to study lithospheric accretion processes in MORs characterized by contrasting tectono-magmatic settings: the magmatically dominated EPR axis between 9°30'-10°00’N, and the Kane Oceanic Core Complex (KOCC), a section of MAR lithosphere (23°20’-23°38’N) formed by detachment faulting. At the EPR study area, I found that the axial magma chamber (AMC) melt sill is segmented into four prominent 2-4-km-long sections spaced every ~5- 10 km along the ridge axis characterized by high melt content (>95%). In contrast, within the intervening sections, the AMC sill has a lower melt content (41-46%). The total magma volume extracted from the AMC sill was estimated of ~46 × 106 m3, with ~24 × 106 m3 left unerupted in the upper crust as dikes after 2005-06 eruption. At the KOCC, I used streamer tomography to constrain the shallow seismic velocity structure. Lithological interpretation of the seismic tomographic models provides insights into the temporal and spatial evolution of the melt supply at the spreading axis as the KOCC formed and evolved. Investigation of a magnetic polarity reversal boundary in crosssection at the northern boundary of KOCC suggests that the boundary (representing both a frozen isotherm and an isochron) dips away from the ridge axis along the Kane transform fault scarp, with a west-dipping angle of ~45° in the shallow (<1 km) crust and <20° in the deeper crust.This thesis was funded by National Science Foundation grants OCE-9987004, OCE- 0621660 and OCE-0327885, WHOI Academic Program Office and WHOI Deep Ocean Exploration Institute

Advanced geophysical studies of accretion of oceanic lithosphere in Mid-Ocean Ridges characterized by contrasting tectono-magmatic settings

Thesis (Ph. D.)--Joint Program in Oceanography (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2012.Cataloged from PDF version of thesis.Includes bibliographical references.The structure of the oceanic lithosphere results from magmatic and extensional processes taking place at mid-ocean ridges (MORs). The temporal and spatial scales of the variability of these two processes control the degree of heterogeneity of the oceanic lithosphere, represented by two end-member models: the classical Penrose Model exemplified by layered magmatic crust formed along fast-spreading MORs, e.g., East Pacific Rise (EPR); and the recently defined Chapman Model describing heterogeneous mafic and ultramafic lithosphere formed in settings of oceanic detachment faulting common along slow-spreading MORs, e.g., Mid-Atlantic Ridge (MAR). This thesis is using advanced marine geophysical methods (including finite-difference wave propagation modeling, 3D multi-channel seismic reflection imaging, waveform inversion, streamer tomography, and near-bottom magnetics) to study lithospheric accretion processes in MORs characterized by contrasting tectono-magmatic settings: the magmatically dominated EPR axis between 9°30'-10°00°N, and the Kane Oceanic Core Complex (KOCC), a section of MAR lithosphere (23°20°-23°38°N) formed by detachment faulting. At the EPR study area, I found that the axial magma chamber (AMC) melt sill is segmented into four prominent 2-4-km-long sections spaced every -5- 10 km along the ridge axis characterized by high melt content (>95%). In contrast, within the intervening sections, the AMC sill has a lower melt content (41-46%). The total magma volume extracted from the AMC sill was estimated of ~46 x 106 M3, with ~24 x 106 M3 left unerupted in the upper crust as dikes after 2005-06 eruption. At the KOCC, I used streamer tomography to constrain the shallow seismic velocity structure. Lithological interpretation of the seismic tomographic models provides insights into the temporal and spatial evolution of the melt supply at the spreading axis as the KOCC formed and evolved. Investigation of a magnetic polarity reversal boundary in crosssection at the northern boundary of KOCC suggests that the boundary (representing both a frozen isotherm and an isochron) dips away from the ridge axis along the Kane transform fault scarp, with a west-dipping angle of ~45° in the shallow (<1 km) crust and <20° in the deeper crust.by Min Xu.Ph.D

Seismic and magnetic constraints on the strucutre of upper oceanic crust and fast and slow spreading ridges

Thesis (Ph. D.)--Joint Program in Oceanography (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), September 1998.Includes bibliographical references (leaves 226-246).The upper ocean crust contains a comprehensive record of the shallow geological processes active along the world's mid-ocean ridge system. This thesis examines the magnetic and seismic structure of the upper crust at two contrasting ridges-the fast spreading East Pacific Rise (EPR) and the slow spreading Mid-Atlantic Ridge (MAR)-to build a more complete understanding about the roles of volcanic emplacement, tectonic disruption and hydrothermal alteration in the near-ridge environment. A technique that inverts potential field measurements directly from an uneven observation track is developed and applied to near-bottom magnetic data from the spreading segments north of the Kane transform on the MAR. It is concluded that the central anomaly magnetization high marks the locus of focused volcanic emplacement. A cyclic faulting model is proposed to explain the oscillatory magnetization pattern associated with discrete blocks of crust being transported out of the rift valley between intensely altered fault zones. Seismic waveform and amplitude analyses of the magma sill along the EPR reveal it to be a thin (<100 m) body of partial melt. These characteristics have important implications for melt availability and transport within the cycle of eruption and replenishment. A genetic algorithm-based seismic waveform inversion technique is developed and applied to on- and near-axis multichannel data from 17'20'S on the EPR and the spreading segment south of the Oceanographer transform (MAR) to map and compare for the first time the detailed velocity structure of the upper crust at two different spreading rates. Combined with conventionally processed seismic profiles, our results show that, while final extrusive thickness is comparable at all spreading ridges (300-500 m), the style of thickening may vary. While a thin (<100 m) extrusive carapace quadruples in thickness within 1-4 km of the EPR crest, the extrusive section at the MAR achieves its final thickness within the inner valley. Both show evidence for a narrow zone of volcanic emplacement. Vigorous hydrothermalism at the EPR may produce a more rapid increase in basement velocities relative to the MAR. Rapid modification of the extrusive/dike transition at both ridges indicates that hydrothermalism is enhanced in this interval. Along-axis transport of lavas may thicken the extrusive pile at slow spreading segment ends, strengthening the magnetic highs generated by lava chemistry.by Stefan Anthony Hussenoeder.Ph.D

Seismic Waves

The importance of seismic wave research lies not only in our ability to understand and predict earthquakes and tsunamis, it also reveals information on the Earth's composition and features in much the same way as it led to the discovery of Mohorovicic's discontinuity. As our theoretical understanding of the physics behind seismic waves has grown, physical and numerical modeling have greatly advanced and now augment applied seismology for better prediction and engineering practices. This has led to some novel applications such as using artificially-induced shocks for exploration of the Earth's subsurface and seismic stimulation for increasing the productivity of oil wells. This book demonstrates the latest techniques and advances in seismic wave analysis from theoretical approach, data acquisition and interpretation, to analyses and numerical simulations, as well as research applications. A review process was conducted in cooperation with sincere support by Drs. Hiroshi Takenaka, Yoshio Murai, Jun Matsushima, and Genti Toyokuni