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: Reciprocity-gap misfit functional for distributed acoustic sensing, combining data from passive and active sources
Quantitative imaging of subsurface earth properties in elastic media is performed from distributed acoustic sensing data. A new misfit functional based upon the reciprocity gap is designed, taking crosscorrelations of displacement and strain, and these products further associate an observation with a simulation. In comparison with other misfit functionals, this functional has the advantage of only requiring little a priori information on the exciting sources. In particular, the misfit criterion enables the use of data from regional earthquakes (teleseismic events can be included as well), followed by exploration data to perform a multiresolution reconstruction. The data from regional earthquakes contain the low-frequency content that is missing in the exploration data, allowing for the recovery of the long spatial wavelength, even with very few sources. These data are used to build prior models for the subsequent reconstruction from the higher frequency exploration data. This results in the elastic full reciprocity-gap waveform inversion method, and we illustrate its performance with a pilot experiment for elastic isotropic reconstruction.  more » « less
Award ID(s):
1815143
PAR ID:
10499959
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1190
Date Published:
Journal Name:
GEOPHYSICS
Volume:
86
Issue:
2
ISSN:
0016-8033
Format(s):
Medium: X Size: p. R211-R220
Size(s):
p. R211-R220
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    The quantitative reconstruction of subsurface earth properties from the propagation of waves follows an iterative minimization of a misfit functional. In marine seismic exploration, the observed data usually consist of measurements of the pressure field, but dual-sensor devices also provide the normal velocity. Consequently, a reciprocity-based misfit functional is specifically designed, and it defines the full reciprocity-gap waveform inversion (FRgWI) method. This misfit functional provides additional features compared to the more traditional least-squares approaches, in particular, in that the observational and computational acquisitions can be different. Therefore, the positions and wavelets of the sources from which the measurements are acquired are not needed in the reconstruction procedure and, in fact, the numerical acquisition (for the simulations) can be chosen arbitrarily. Based on 3D experiments, FRgWI is shown to behave better than full-waveform inversion in the same context. It allows for arbitrary numerical acquisitions in two ways: when few measurements are given, a dense numerical acquisition (compared to the observational one) can be used to compensate. However, with a dense observational acquisition, a sparse computational one is shown to be sufficient, for instance, with multiple-point sources, hence reducing the numerical cost. FRgWI displays accurate reconstructions in both situations and appears more robust with respect to crosstalk than least-squares shot stacking. 
    more » « less
  2. We introduce a novel approach to waveform inversion based on a data-driven reduced order model (ROM) of the wave operator. The presentation is for the acoustic wave equation, but the approach can be extended to elastic or electromagnetic waves. The data are time resolved measurements of the pressure wave gathered by an acquisition system that probes the unknown medium with pulses and measures the generated waves. We propose to solve the inverse problem of velocity estimation by minimizing the square misfit between the ROM computed from the recorded data and the ROM computed from the modeled data, at the current guess of the velocity. We give a step by step computation of the ROM, which depends nonlinearly on the data and yet can be obtained from them in a noniterative fashion, using efficient methods from linear algebra. We also explain how to make the ROM robust to data inaccuracy. The ROM computation requires the full array response matrix gathered with colocated sources and receivers. However, we find that the computation can deal with an approximation of this matrix, obtained from towed-streamer data using interpolation and reciprocity on-the-fly. Although the full-waveform inversion approach of nonlinear least-squares data fitting is challenging without low-frequency information, due to multiple minima of the data fit objective function, we find that the ROM misfit objective function has better behavior, even for a poor initial guess. We also find by explicit computation of the objective functions in a simple setting that the ROM misfit objective function has convexity properties, whereas the least-squares data fit objective function displays multiple local minima. 
    more » « less
  3. Abstract Acoustic non-reciprocity, referring to the phenomenon of path-dependent propagation, has diverse applications in mechanical devices. This paper presents a numerical study on a periodic triangle-shape structure that breaks reciprocity in a passive manner over a broad range of frequency and energy. The proposed structure contains strong nonlinearity and geometric asymmetry, which contributes to a direction-dependent dispersion relationship. When the signal frequency falls in the band pass in one direction, and band gap in the other, a unidirectional wave propagation results. The system achieves giant non-reciprocity with minimal distortion in the frequency content of the signal. 
    more » « less
  4. Abstract Seismic signals, whether caused by earthquakes, other natural phenomena, or artificial noise sources, have specific characteristics in the time and frequency domains that contain crucial information reflecting their source. The analysis of seismic time series is an essential part of every seismology-oriented study program. Enabling students to work with data collected from their own campus, including signals from both anthropogenic and natural seismic sources, can provide vivid, practical examples to make abstract concepts communicated in classes more concrete and relevant. Data from research-grade broadband seismometers enable us to record time series of vibrations at a broad range of frequencies; however, these sensors are costly and are often deployed in remote places. Participation in the Raspberry Shake citizen science network enables seismology educators to record seismic signals on our own campuses and use these recordings in our classrooms and for public outreach. Yale University installed a Raspberry Shake three-component, low-cost seismometer in the Earth and Planetary Sciences department building in Summer 2022, enabling the detection of local, regional, and teleseismic earthquakes, microseismic noise, and anthropogenic noise sources from building construction, an explosive event in a steam tunnel, and general building use. Here, we discuss and illustrate the use of data from our Raspberry Shake in outreach and education activities at Yale. In particular, we highlight a series of ObsPy-based exercises that will be used in courses taught in our department, including our upper-level Introduction to Seismology course and our undergraduate classes on Natural Disasters and Forensic Geoscience. 
    more » « less
  5. Abstract We investigate non-Hermitian elastic lattices characterized by non-local feedback interactions. In one-dimensional lattices, proportional feedback produces non-reciprocity associated with complex dispersion relations characterized by gain and loss in opposite propagation directions. For non-local controls, such non-reciprocity occurs over multiple frequency bands characterized by opposite non-reciprocal behavior. The dispersion topology is investigated with focus on winding numbers and non-Hermitian skin effect, which manifests itself through bulk modes localized at the boundaries of finite lattices. In two-dimensional lattices, non-reciprocity is associated with directional wave amplification. Moreover, the combination of skin effect in two directions produces modes that are localized at the corners of finite two-dimensional lattices. Our results describe fundamental properties of non-Hermitian elastic lattices, and suggest new possibilities for the design of meta materials with novel functionalities related to selective wave filtering, amplification and localization. The considered non-local lattices also provide a platform for the investigation of topological phases of non-Hermitian systems. 
    more » « less