Nonlinear waveform inversion of plane‐wave seismograms in stratified elastic media

Körmendi, F.F.; Dietrich, Michel

doi:10.1190/1.1443083

Cited by 71 publications

(56 citation statements)

References 27 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…For the remainder of this paper we explore an alternative explanation for the origin of the bright spots that is consistent with their distribution and the history of shelf glaciation. data set using the linearized waveform inversion approach of Kormendi and Dietrich [1991]. This type of analysis has been extensively described in the literature, and our particular approach closely follows that of Korenaga et al [1997], so we have placed details of our analysis in Appendix A.…”

Section: Data Description and Observationsmentioning

confidence: 99%

See 1 more Smart Citation

Glacially generated overpressure on the New England continental shelf: Integration of full‐waveform inversion and overpressure modeling

et al. 2014

View full text Add to dashboard Cite

Localized zones of high-amplitude, discontinuous seismic reflections 100 km off the coast of Massachusetts, USA, have P wave velocities up to 190 m/s lower than those of adjacent sediments of equal depth (250 m below the sea floor). To investigate the origin of these low-velocity zones, we compare the detailed velocity structure across high-amplitude regions to adjacent, undisturbed regions through full-waveform inversion. We relate the full-waveform inversion velocities to effective stress and overpressure with a power law model. This model predicts localized overpressures up to 2.2 MPa associated with the high-amplitude reflections. To help understand the overpressure source, we model overpressure due to erosion, glacial loading, and sedimentation in one dimension. The modeling results show that ice loading from a late Pleistocene glaciation, ice loading from the Last Glacial Maximum, and rapid sedimentation contributed to the overpressure. Localized overpressure, however, is likely the result of focused fluid flow through a high-permeability layer below the region characterized by the high-amplitude reflections. These high overpressures may have also caused localized sediment deformation. Our forward models predict maximum overpressure during the Last Glacial Maximum due to loading by glaciers and rapid sedimentation, but these overpressures are dissipating in the modern, low sedimentation rate environment. This has important implications for our understanding continental shelf morphology, fluid flow, and submarine groundwater discharge off Massachusetts, as we show a mechanism related to Pleistocene ice sheets that may have created regions of anomalously high overpressure.

show abstract

Section: Data Description and Observationsmentioning

confidence: 99%

“…We follow the methods of Kormendi and Dietrich [1991], which apply the Fréchet derivatives of Dietrich and Kormendi [1990] and use a conjugate gradient algorithm to iteratively solve the nonlinear inverse problem.…”

Section: Appendix A: Full-waveform Inversion Modelingmentioning

confidence: 99%

Glacially generated overpressure on the New England continental shelf: Integration of full‐waveform inversion and overpressure modeling

et al. 2014

View full text Add to dashboard Cite

show abstract

“…In the Gauss-Newton method, a locally linear misfit functional is assumed (Kormendi & Dietrich, 1991;Menke, 1989;Aster et al, 2005), which allows the second order term of the Hessian to be dropped (Virieux & Operto, 2009). We obtain explicitly the sensitivity matrix J by perturbing each model parameter at each layer depth; the resulting partial derivative wavefield is propagated from the secondary virtual sources to the receivers' position (Rodi, 1976;Sheen et al, 2006;Operto et al, 2013).…”

Section: Gauss-newton Seismic Inversionmentioning

confidence: 99%

“…The most desirable parametrisation guarantees the minimum crosstalk among the unknowns of the inversion (Tarantola, 1986;Kormendi & Dietrich, 1991); ideally, the partial derivative wavefield of one parameter should be uncorrelated with the residual wavefield produced by each other independent parameter (Tarantola, 1986;Operto et al, 2013). In marine reflection seismic data, the presence of only one propagation mode impedes the opportunity to obtain independent estimates of P-wave (Vp) and S-wave velocity (Vs) (Jin et al, 1992;Igel et al, 1996); on the other hand, density is strongly coupled with P-wave velocity at narrow reflection angle; the two parameters can't be effectively resolved and in fact yield a posterior reconstruction of the P-impedance model (Tarantola, 1986;Operto et al, 2013).…”

Section: A Strategy For the Multi-parameter Problemmentioning

confidence: 99%

Pre-stack full waveform inversion of ultra-high-frequency marine seismic reflection data

Provenzano

Vardy

Henstock

2017

Geophysical Journal International

View full text Add to dashboard Cite

SUMMARYThe full waveform inversion (FWI) of seismic reflection data aims to reconstruct a detailed physical properties model of the subsurface, fitting both the amplitude and traveltime of the reflections generated at physical discontinuities in the propagation medium. Unlike reservoirscale seismic exploration, where seismic inversion is a widely adopted remote characterisation tool, ultra high frequency (UHF, 0.2-4.0 kHz) multi-channel marine reflection seismology is still most often limited to a qualitative interpretation of the reflections' architecture. Here we propose an elastic full waveform inversion methodology, custom-tailored for pre-stack UHF marine data in vertically heterogeneous media to obtain a decimetric-scale distribution of Pimpedance, density and Poisson's ratio within the shallow sub-seabed sediments. We address the deterministic multi-parameter inversion in a sequential fashion. The complex trace instantaneous phase is first inverted for the P-wave velocity to make-up for the lack of low-frequency in the data and reduce the non-linearity of the problem. This is followed by a short-offset Pimpedance optimisation and a further step of full offset range Poisson's ratio inversion. Provided that the seismogram contains wide reflection angles (> 40 degrees), we show that it is possible to invert for density and decompose a-posteriori the relative contribution of P-wave velocity and density to the P-impedance. A broad range of synthetic tests is used to prove the potential of the methodology and highlights sensitivity issues specific to UHF seismic. An example application to real data is also presented. In the real case, trace normalisation is applied to minimise the systematic error deriving from an inaccurate source wavelet estimation. G. Provenzano et al.The inverted model for the top 15 meters of the sub-seabed agrees with the local lithological information and core-log data. Thus we can obtain a detailed remote characterisation of the shallow sediments using a multi-channel sub-bottom profiler within a reasonable computing cost and with minimal pre-processing. This has the potential to reduce the need of extensive geotechnical coring campaigns.

show abstract

“…The full waveform inversion method is described in detail by Kormendi and Dietrich [1991], and more details on the inversion procedure can be found in , , and Minshull et al [1994]. We give only a brief outline here.…”

Section: Full Waveform Inversionmentioning

confidence: 99%

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

Xu¹

2012

View full text Add to dashboard Cite

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 × 10 6 m 3 , with ~24 × 10 6 m 3 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. Many thanks to my project advisors and cooperators, Brian Tucholke, RalphStephen, and Maurice Tivey, who have worked with me during the past five and a half years. These research experiences have greatly broadened my horizons on marine geology and geophysics, synthetic seismology, and geomagnetism. Brian introduced me into the geological world of detachment fault system, and his patience and insightful comments on our scientific paper on Kane oceanic core complex are greatly appreciated.I am very thankful for Ralph who introduced me to the world of time-domain finite difference (TDFD) synthetic seismology and Tom Bolmer for setting up the workstation for my TDFD calculations. Ralph's encouragement, enthusiasm, and interactions with me have inspired me to dig out more intere...

show abstract

Nonlinear waveform inversion of plane‐wave seismograms in stratified elastic media

Cited by 71 publications

References 27 publications

Glacially generated overpressure on the New England continental shelf: Integration of full‐waveform inversion and overpressure modeling

Glacially generated overpressure on the New England continental shelf: Integration of full‐waveform inversion and overpressure modeling

Pre-stack full waveform inversion of ultra-high-frequency marine seismic reflection data

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

Contact Info

Product

Resources

About