skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Geotechnical site characterization with 3D ambient noise tomography
We develop a new 3D ambient noise tomography (3D ANT) method for geotechnical site characterization. It requires recording ambient noise fields using a 2D surface array of geophones, from which experimental crosscorrelation functions (CCFs) are then extracted and directly inverted to obtain an S-wave velocity ([Formula: see text]) structure. The method consists of a forward simulation using 3D P-SV elastic wave equations to compute the synthetic CCF and an adjoint-state inversion to match the synthetic CCFs to the experimental CCFs for extraction of [Formula: see text]. The main advantage of the presented method, as compared with conventional passive-source seismic methods using characteristics of Green’s function (GF), is that it does not require equal energy on both sides of each receiver pair or far-field wavefields to retrieve the true GF. Instead, the source power spectrum density is inverted during the analysis and incorporated into the forward simulation of the synthetic CCFs to account for source energy distribution. After testing on synthetic data, the 3D ANT method is applied to 3 h of ambient noise recordings at the Garner Valley Downhole Array (GVDA) site in California, using a surface array of 196 geophones placed on a 14 × 14 grid with 5 m spacing. The inverted 3D [Formula: see text] model is found to be consistent with previous invasive and noninvasive geotechnical characterization efforts at the GVDA site.  more » « less
Award ID(s):
1930697 2037900
PAR ID:
10457303
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
GEOPHYSICS
Volume:
88
Issue:
4
ISSN:
0016-8033
Page Range / eLocation ID:
KS101 to KS112
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Geotechnical characterization of marine sediments remains an outstanding challenge for offshore energy development, including foundation design and site selection of wind turbines and offshore platforms. We demonstrate that passive distributed acoustic sensing (DAS) surveys offer a new solution for shallow offshore geotechnical investigation where seafloor power or communications cables with fiber-optic links are available. We analyze Scholte waves recorded by DAS on a 42 km power cable in the Belgian offshore area of the southern North Sea. Ambient noise crosscorrelations converge acceptably with just over one hour of data, permitting multimodal Scholte wave dispersion measurement and shear-wave velocity inversion along the cable. We identify anomalous off-axis Scholte wave arrivals in noise crosscorrelations at high frequencies. Using a simple passive source imaging approach, we associate these arrivals with individual wind turbines, which suggests they are generated by structural vibrations. While many technological barriers must be overcome before ocean-bottom DAS can be applied to global seismic monitoring in the deep oceans, high-frequency passive surveys for high-resolution geotechnical characterization and monitoring in coastal regions are easily achievable today. 
    more » « less
  2. Polycrystalline materials consist of grains (crystals) oriented at different angles resulting in a heterogeneous and anisotropic mechanical behavior at that micro-length scale. In this study, a novel method is proposed for the first time to determine the [Formula: see text] crystal orientations of grains in a [Formula: see text] domain, using solely [Formula: see text] deformation fields. The grain boundaries are assumed to be unknown and delineated from the reconstructed changes in the crystallographic orientation. Further, the constitutive equations that describe the mechanical behavior of the domain in [Formula: see text] under plane stress conditions are derived, assuming that the material is transversely isotropic in 3D. Finite element based algorithms are utilized to discretize the inverse problem. The in-house written inverse problem solver is coupled with Matlab-based optimization scripts to solve for the mechanical property distributions. The performance of this method is tested at different noise levels with synthetic displacements that were used as measured data. The reconstructions deteriorate as the noise level is increased. This work presents a first milestone in the verification of this novel technology with synthetic data. 
    more » « less
  3. SUMMARY For seismographic stations with short acquisition duration, the signal-to-noise ratios (SNRs) of ambient noise cross-correlation functions (CCFs) are typically low, preventing us from accurately measuring surface wave dispersion curves or waveform characteristics. In addition, with noisy CCFs, it is difficult to extract relatively weak signals such as body waves. In this study, we propose to use local attributes to improve the SNRs of ambient noise CCFs, which allows us to enhance the quality of CCFs for stations with limited acquisition duration. Two local attributes: local cross-correlation and local similarity, are used in this study. The local cross-correlation allows us to extend the dimensionality of daily CCFs with computational costs similar to global cross-correlation. Taking advantage of this extended dimensionality, the local similarity is then used to measure non-stationary similarity between the extended daily CCFs with a reference trace, which enables us to design better stacking weights to enhance coherent features and attenuate incoherent background noises. Ambient noise recorded by several broad-band stations from the USArray in North Texas and Oklahoma, the Superior Province Rifting EarthScope Experiment in Minnesota and Wisconsin and a high-frequency nodal array deployed in the northern Los Angeles basin are used to demonstrate the performance of the proposed approach for improving the SNR of CCFs. 
    more » « less
  4. ABSTRACT The ability to monitor seismicity and structural integrity of a mine using seismic noise can have great implication for detecting and managing ground-control hazards. The noise wavefield, however, is complicated by induced seismicity and heavy machinery associated with mining operations. In this study, we investigate the nature of time-dependent noise cross-correlations functions (CCFs) across an active underground longwall coal mine. We analyze one month of continuous data recorded by a surface 17 geophone array with an average station spacing of ∼200 m. To extract coherent seismic signals, we calculate CCFs between all stations for each 5-min window. Close inspection of all 5-min CCFs reveals waveforms that can be categorically separated into two groups, one with strong and coherent 1–5 Hz signals and one without. Using a reference station pair, we statistically isolate time windows within each group based on the correlation coefficient between each 5-min CCF and the monthly stacked CCF. The daily stacked CCFs associated with a high correlation coefficient show a clear temporal variation that is consistent with the progression of mining activity. In contrast, the daily stacked CCFs associated with a low correlation coefficient remain stationary throughout the recording period in line with the expected persistent background noise. To further understand the nature of the high correlation coefficient CCFs, we perform 2D and 3D back projection to determine and track the dominant noise source location. Excellent agreement is observed on both short (5-min) and long (daily) time scales between the CCF determined source locations, the overall migration of the active mining operation, and cataloged seismic event locations. The workflow presented in this study demonstrates an effective way to identify and track mining induced signals, in which CCFs associated with background noise can be isolated and used for further temporal structural integrity investigation. 
    more » « less
  5. We have formulated the Fréchet kernel computation using the adjoint-state method based on a fractional viscoacoustic wave equation. We first numerically prove that the 1/2- and the 3/2-order fractional Laplacian operators are self-adjoint. Using this property, we find that the adjoint wave propagator preserves the dispersion and compensates the amplitude, whereas the time-reversed adjoint wave propagator behaves identically to the forward propagator with the same dispersion and dissipation characters. Without introducing rheological mechanisms, this formulation adopts an explicit [Formula: see text] parameterization, which avoids the implicit [Formula: see text] in the conventional viscoacoustic/viscoelastic full-waveform inversion ([Formula: see text]-FWI). In addition, because of the decoupling of operators in the wave equation, the viscoacoustic Fréchet kernel is separated into three distinct contributions with clear physical meanings: lossless propagation, dispersion, and dissipation. We find that the lossless propagation kernel dominates the velocity kernel, whereas the dissipation kernel dominates the attenuation kernel over the dispersion kernel. After validating the Fréchet kernels using the finite-difference method, we conduct a numerical example to demonstrate the capability of the kernels to characterize the velocity and attenuation anomalies. The kernels of different misfit measurements are presented to investigate their different sensitivities. Our results suggest that, rather than the traveltime, the amplitude and the waveform kernels are more suitable to capture attenuation anomalies. These kernels lay the foundation for the multiparameter inversion with the fractional formulation, and the decoupled nature of them promotes our understanding of the significance of different physical processes in [Formula: see text]-FWI. 
    more » « less