We present a new 3‐D seismic structural model of the eastern Indonesian region and its surroundings from full‐waveform inversion (FWI) that exploits seismic data filtered at periods between 15–150 s.
Non-invasive subsurface imaging using full waveform inversion (FWI) has the potential to fundamentally change near-surface (<30 m) site characterization by enabling the recovery of high-resolution (metre-scale) 2-D/3-D maps of subsurface elastic material properties. Yet, FWI results are quite sensitive to their starting model due to their dependence on local-search optimization techniques and inversion non-uniqueness. Starting model dependence is particularly problematic for near-surface FWI due to the complexity of the recorded seismic wavefield (e.g. dominant surface waves intermixed with body waves) and the potential for significant spatial variability over short distances. In response, convolutional neural networks (CNNs) are investigated as a potential tool for developing starting models for near-surface 2-D elastic FWI. Specifically, 100 000 subsurface models were generated to be representative of a classic near-surface geophysics problem; namely, imaging a two-layer, undulating, soil-over-bedrock interface. A CNN has been developed from these synthetic models that is capable of transforming an experimental wavefield acquired using a seismic source located at the centre of a linear array of 24 closely spaced surface sensors directly into a robust starting model for FWI. The CNN approach was able to produce 2-D starting models with seismic image misfits that were significantly less than the misfits from other common starting model approaches, and in many cases even less than the misfits obtained by FWI with inferior starting models. The ability of the CNN to generalize outside its two-layered training set was assessed using a more complex, three-layered, soil-over-bedrock formation. While the predictive ability of the CNN was slightly reduced for this more complex case, it was still able to achieve seismic image and waveform misfits that were comparable to other commonly used starting models, despite not being trained on any three-layered models. As such, CNNs show great potential as tools for rapidly developing robust, site-specific starting models for near-surface elastic FWI.
more » « less- Award ID(s):
- 2120155
- NSF-PAR ID:
- 10412842
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 231
- Issue:
- 1
- ISSN:
- 0956-540X
- Page Range / eLocation ID:
- p. 72-90
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract SASSY21 —a recent 3‐D FWI tomographic model of Southeast Asia—is used as a starting model, and our study region is characterized by particularly good data coverage, which facilitates a more refined image. We use the spectral‐element solverSalvus to determine the full 3‐D wavefield, accounting for the fluid ocean explicitly by solving a coupled system of acoustic and elastic wave equations. This is computationally more expensive but allows seismic waves within the water layer to be simulated, which becomes important for periods ≤20 s. We investigate path‐dependent effects of surface elevation (topography and bathymetry) and the fluid ocean on synthetic waveforms, and compare our final model to the tomographic result obtained with the frequently used ocean loading approximation. Furthermore, we highlight some of the key features of our final model—SASSIER22 —after 34 L‐BFGS iterations, which reveals detailed anomalies down to the mantle transition zone, including a convergent double‐subduction zone along the southern segment of the Philippine Trench, which was not evident in the starting model. A more detailed illumination of the slab beneath the North Sulawesi Trench reveals a pronounced positive wavespeed anomaly down to 200 km depth, consistent with the maximum depth of seismicity, and a more diffuse but aseismic positive wavespeed anomaly that continues to the 410 km discontinuity. -
Full waveform inversion (FWI) and distributed acoustic sensing (DAS) are powerful tools with potential to improve how seismic site characterization is performed. FWI is able to provide true 2D or 3D images of the subsurface by inverting stress wave recordings collected over a wide variety of scales. DAS can be used to efficiently collect high-resolution stress wave recordings from long and complex fiber optic arrays and is well-suited for large-scale site characterization projects. Due to the relative novelty of combining FWI and DAS, there is presently little published literature regarding the application of FWI to DAS data for near-surface (depths < 30 m) site characterization. We perform 2D FWI on DAS data collected at a well-characterized site using four different, site-specific 1D and 2D starting models. We discuss the unique benefits and challenges associated with inverting DAS data compared to traditional geophone data. We examine the impacts of using the various starting models on the final 2D subsurface images. We demonstrate that while the inversions performed using all four starting models are able to fit the major features of the DAS waveforms with similar misfit values, the final subsurface images can be quite different from one another at depths greater than about 10 m. As such, the best representation(s) of the subsurface are evaluated based on: (1) their agreement with borehole lithology logs that were not used in the development of the starting models, and (2) consistency at shallow depths between the final inverted images derived from multiple starting models. Our results demonstrate that FWI applied to DAS data has significant potential as a tool for near-surface site characterization while also emphasizing the significant impact that starting model selection can have on FWI results.more » « less
-
Sinkhole collapse may result in significant property damage and even loss of life. Early detection of sinkhole attributes (buried voids, raveling zones) is critical to limit the cost of remediation. One of the most promising ways to obtain subsurface imaging is 3D seismic full-waveform inversion. For demonstration, a recently developed 3D Gauss-Newton full-waveform inversion (3D GN-FWI) method is used to detect buried voids, raveling soils, and characterize variable subsurface soil/rock layering. It is based on a finite-difference solution of 3D elastic wave equations and Gauss-Newton optimization. The method is tested first on a data set constructed from the numerical simulation of a challenging synthetic model and subsequently on field data collected from two separate test sites in Florida. For the field tests, receivers and sources are placed in uniform 2D surface grids to acquire the seismic wavefields, which then are inverted to extract the 3D subsurface velocity structures. The inverted synthetic results suggest that the approach is viable for detecting voids and characterizing layering. The field seismic results reveal that the 3D waveform analysis identified a known manmade void (plastic culvert), unknown natural voids, raveling, as well as laterally variable soil/rock layering including rock pinnacles. The results are confirmed later by standard penetration tests, including depth to bedrock, two buried voids, and a raveling soil zone. Our study provides insight into the application of the 3D seismic FWI technique as a powerful tool in detecting shallow voids and other localized subsurface features.more » « less
-
more » « less
ABSTRACT The use of surface wave testing for near‐surface engineering applications has increased in recent years. Typical surface wave analysis is based on the dispersion of surface waves in one‐dimensional layered models. One‐dimensional models are inappropriate for measurements at sites with appreciable lateral variability, a likely scenario in many engineering applications. Use of such models can subsequently undermine the reliability and accuracy of the surface wave results. Full waveform inversion (FWI) is a high‐resolution imaging technique that is proven to outperform the conventional dispersion‐based analysis of surface waves. Much of the near‐surface literature has focused on full waveform inversion of Rayleigh waves developed by the interaction of primary‐ and vertically polarized shear waves (P‐SV), and the capabilities of surface waves generated by horizontally polarized shear waves (Love waves) in mapping near‐surface spatial variabilities have not been fully explored. In this numerical study, full waveform inversion of Rayleigh and Love waves was performed on two different spatially correlated Gaussian random fields (mean
VS of 200 and 500 m/s) to mimic the natural spatial variability of geologic materials. Each soil structure was produced at a low and high level of stiffness variability. Two sources with different frequency contents, 25 and 50 Hz, were used to evaluate the effects of source characteristics on the resolution of Rayleigh and Love waveform inversions. The inverted results from the high‐velocity domain demonstrated that Love waveform inversion using high‐frequency seismic sources outperforms Rayleigh full waveform inversion in detecting the shape and the velocities of horizontally deposited geologic materials. Results from the low‐velocity domain also confirmed that Love full waveform inversion was comparable or superior to Rayleigh full waveform inversion, though the performance difference was less significant. However, the 25‐Hz frequency inversions yielded superior results than the 50‐Hz frequency inversions for the low‐velocity domain because the dominant wavelength of the high‐frequency signals becomes so small that it offers an impractically small investigation depth. -
more » « less
Abstract We present the first continental‐scale seismic model of the lithosphere and underlying mantle beneath Southeast Asia obtained from adjoint waveform tomography (often referred to as full‐waveform inversion or FWI), using seismic data filtered at periods from 20 to 150 s. Based on >3,000 hr of analyzed waveform data gathered from ∼13,000 unique source‐receiver pairs, we image isotropic
P ‐wave velocity, radially anisotropicS ‐wave velocity and density via an iterative non‐linear inversion that begins from a 1‐D reference model. At each iteration, the full 3‐D wavefield is determined through an anelastic Earth, accommodating effects of topography, bathymetry and ocean load. Our data selection aims to maximize sensitivity to deep structure by accounting for body wave arrivals separately.SASSY21 , our final model after 87 iterations across seven period bands, is able to explain true‐amplitude data from events and receivers not included in the inversion. The trade‐off between inversion parameters is estimated through an analysis of the Hessian‐vector product.SASSY21 reveals detailed anomalies down to the mantle transition zone, including multiple subduction zones. The most prominent feature is the (Indo‐)Australian plate descending beneath Indonesia, which is imaged as one continuous slab along the 180° curvature of the Banda Arc. The tomography confirms the existence of a hole in the slab beneath Mount Tambora and locates a highS ‐wave velocity zone beneath northern Borneo that may be associated with subduction termination in the mid‐late Miocene. A previously undiscovered feature beneath the east coast of Borneo is also revealed, which may be a signature of post‐subduction processes, delamination or underthrusting from the formation of Sulawesi.
