Search

Distributed acoustic sensing (DAS) is a relatively new technology for recording the propagation of seismic waves, with promising applications in both engineering and geophysics. DAS's ability to simultaneously collect high spatial resolution waveforms over long arrays suggests that it is well-suited for near-surface imaging applications such as 2-D multichannel analysis of surface waves (MASWs), which require, at a minimum, long, linear arrays of single-component receivers. The 2-D MASW method uses a large number of sensor subarrays deployed along a linear alignment to produce 1-D shear-wave velocity (VS) profiles beneath each subarray. The 1-D VS profiles are then combined to form a pseudo-2-D VS image beneath the entire linear alignment that can be used for the purpose of identifying and characterizing lateral variations in subsurface layering. Traditionally, 2-D MASW is conducted using arrays consisting of either 24 or 48 geophones. While additional receivers could easily be incorporated into the testing configuration, it is rare for researchers and practitioners to have access to greater numbers of seismographs and geophones. When a limited number of geophones are available for deployment, there is a need to pre-determine the geophone spacing and subarray length prior to field data acquisition. Studies examining how the choice of subarray geometry impacts the resulting pseudo-2-D VS cross-sections have been largely limited to synthetic data. In response, this study utilizes DAS data to examine the effects of using various subarray lengths by comparing pseudo-2-D VS cross-sections derived from active-source waveforms collected at a well-characterized field site. DAS is particularly useful for 2-D MASW applications because the subarray geometry does not need to be determined prior to field data acquisition. We organize the DAS waveforms into multiple sets of overlapping MASW subarrays of differing lengths, ranging from 11 to 47 m, along the same alignment, allowing for direct comparison of the derived pseudo-2-D VS results at the site. We show that the length of the individual MASW subarrays has a significant effect on the resulting VS cross-sections, including the resolved location of a strong impedance contrasts at our study site, and evaluate the results relative to ground truth from invasive testing. Our results suggest that the choice of subarray length is important and should be carefully chosen to meet project-specific goals. Furthermore, analysts may consider using multiple subarray geometries during the data processing stage, as is made possible by DAS, to properly evaluate the uncertainty of 2-D MASW results. This study demonstrates the potential of using DAS to collect data for 2-D MASW in a manner that is efficient and flexible, and can be easily scaled up for use with very long arrays. [ABSTRACT FROM AUTHOR]

SWinvert is a workflow developed at The University of Texas at Austin for the inversion of surface wave dispersion data. SWinvert encourages analysts to investigate inversion uncertainty and non-uniqueness in shear wave velocity (Vs) by providing a systematic procedure and specific actionable recommendations for surface wave inversion. In particular, the workflow encourages the use of multiple layering parametrizations to address the inversion's non-uniqueness, multiple global searches for each parametrization to address the inverse problem's non-linearity and quantification of Vs uncertainty in the resulting profiles. While the workflow uses the Dinver module of the popular open-source Geopsy software as its inversion engine, the principles presented are of relevance to analysts using other inversion programs. To illustrate the effectiveness of the SWinvert workflow and to develop a set of benchmarks for use in future surface wave inversion studies, synthetic experimental dispersion data for 12 subsurface models of varying complexity are inverted. While the effects of inversion uncertainty and non-uniqueness are shown to be minimal for simple subsurface models characterized by broad-band dispersion data, these effects cannot be ignored in the Vs profiles derived for more complex models with band-limited dispersion data. To encourage adoption of the SWinvert workflow, an open-source Python package (SWprepost), for pre- and post-processing of surface wave inversion data, and an application on the DesignSafe-Cyberinfrastructure (SWbatch), for performing batch-style surface wave inversions with Dinver using high-performance computing, have been developed and released in conjunction with this work. The SWinvert workflow is shown to provide a methodical procedure and a powerful set of tools for performing rigorous surface wave inversions and quantifying the uncertainty in the resulting Vs profiles. [ABSTRACT FROM AUTHOR]

Surface wave methods provide a cost effective means of developing shear wave velocity (Vs) profiles for applications such as dynamic site characterization and seismic site response analyses. However, the inverse problem involved in obtaining a realistic layered earth model from surface wave dispersion data is inherently ill-posed, non-linear and mix-determined, without a unique solution. When available, a priori information such as geotechnical boreholes or geologic well logs should be used to aid in constraining site-specific inversion parameters. Unfortunately, a priori information is often unavailable, particularly at significant depths, and a 'blind analysis' must be performed. In these situations, the analyst must decide on an appropriate number of layers and ranges for their corresponding inversion parameters (i.e. trial number of layers and ranges in their respective thicknesses, shear wave velocities, compression wave velocities and mass densities). Selection of these parameters has been shown to significantly impact the results of an inversion. This paper presents a method for conducting multiple inversions utilizing systematically varied inversion layering parametrizations in order to identify and encompass the most reasonable layered earth models for a site. Each parametrization is defined by a unique layering ratio, which represents a multiplier that systemically increases the potential thickness of each layer in the inversion parametrization based on the potential thickness of the layer directly above it. The layering ratio method is demonstrated at two sites associated with the InterPacific Project, wherein it is shown to significantly aid in selecting reasonable Vs profiles that are close representations of the subsurface. While the goal of the layering ratio inversion methodology is not necessarily to find the 'optimal' or 'best' Vs profile for a site, it may be successful at doing so for certain sites/datasets. However, the primary reason for using the layering ratio method is to find Vs profiles that realistically represent the uncertainty in Vs resulting from surface wave inversion, and to avoid selection of Vs profiles that are unrealistic and adversely influenced by the choice of inversion parametrization. [ABSTRACT FROM AUTHOR]