Search

We present a technique to automatically classify the wave type of seismic phases that are recorded on a single six-component recording station (measuring both three components of translational and rotational ground motion) at the earth's surface. We make use of the fact that each wave type leaves a unique 'fingerprint' in the six-component motion of the sensor. This fingerprint can be extracted by performing an eigenanalysis of the data covariance matrix, similar to conventional three-component polarization analysis. To assign a wave type to the fingerprint extracted from the data, we compare it to analytically derived six-component polarization models that are valid for pure-state plane wave arrivals. For efficient classification, we make use of the supervised machine learning method of support vector machines that is trained using data-independent, analytically-derived six-component polarization models. This enables the rapid classification of seismic phases in a fully automated fashion, even for large data volumes, such as encountered in land-seismic exploration or ambient noise seismology. Once the wave-type is known, additional wave parameters (velocity, directionality, and ellipticity) can be directly extracted from the six-component polarization states without the need to resort to expensive optimization algorithms. We illustrate the benefits of our approach on various real and synthetic data examples for applications such as automated phase picking, aliased ground-roll suppression in land-seismic exploration, and the rapid close-to real time extraction of surface wave dispersion curves from single-station recordings of ambient noise. Additionally, we argue that an initial step of wave type classification is necessary in order to successfully apply the common technique of extracting phase velocities from combined measurements of rotational and translational motion., Comment: Preprint submitted to Geophysical Journal International

Mitigating the energy requirements of artificial intelligence requires novel physical substrates for computation. Phononic metamaterials have a vanishingly low power dissipation and hence are a prime candidate for green, always-on computers. However, their use in machine learning applications has not been explored due to the complexity of their design process: Current phononic metamaterials are restricted to simple geometries (e.g. periodic, tapered), and hence do not possess sufficient expressivity to encode machine learning tasks. We design and fabricate a non-periodic phononic metamaterial, directly from data samples, that can distinguish between pairs of spoken words in the presence of a simple readout nonlinearity; hence demonstrating that phononic metamaterials are a viable avenue towards zero-power smart devices., Comment: 13 pages, 11 figures

Rendering objects invisible to impinging acoustic waves (cloaking) and creating acoustic illusions (holography) has been attempted using active and passive approaches. While passive methods are applicable only to narrow frequency bands, active approaches attempt to respond dynamically, interfering with broadband incident or scattered wavefields by emitting secondary waves. Without prior knowledge of the primary wavefield, the signals for the secondary sources need to be estimated and adapted in real-time. This has thus far impeded active cloaking and holography for broadband wavefields. We present experimental results of active acoustic cloaking and holography without prior knowledge of the wavefield so that objects remain invisible and illusions intact even for broadband moving sources. This opens novel research directions and facilitates practical applications including architectural acoustics, and stealth.

Finite-difference (FD) modeling of seismic waves in the vicinity of dipping interfaces gives rise to artifacts. Examples are phase and amplitude errors, as well as staircase diffractions. Such errors can be reduced in two general ways. In the first approach, the interface can be anti-aliased (i.e., with an anti-aliased step-function, or a lowpass filter). Alternatively, the interface may be replaced with an equivalent medium (i.e., using Schoenberg \& Muir (SM) calculus or orthorhombic averaging). We test these strategies in acoustic, elastic isotropic, and elastic anisotropic settings. Computed FD solutions are compared to analytical solutions. We find that in acoustic media, anti-aliasing methods lead to the smallest errors. Conversely, in elastic media, the SM calculus provides the best accuracy. The downside of the SM calculus is that it requires an anisotropic FD solver even to model an interface between two isotropic materials. As a result, the computational cost increases compared to when using isotropic FD solvers. However, since coarser grid spacings can be used to represent the dipping interfaces, the two effects (an expensive FD solver on a coarser FD grid) equal out. Hence, the SM calculus can provide an efficient means to reduce errors, also in elastic isotropic media., Comment: 34 pages, 18 figures

We describe a method for removing the numerical errors in the modeling of linear evolution equations that are caused by approximating the time derivative by a finite difference operator. The method is based on integral transforms realized as certain Fourier integral operators, called time dispersion transforms, and we prove that, under an assumption about the frequency content, it yields a solution with correct evolution throughout the entire lifespan. We demonstrate the method on a model equation as well as on the simulation of elastic and viscoelastic wave propagation., Comment: 41 pages, 8 figures

A major focus of research in the seismic industry of the past two decades has been the acquisition and subsequent separation of seismic data using multiple sources fired simultaneously. The recently introduced method of {\it signal apparition} provides a new take on the problem by replacing the random time-shifts usually employed to encode the different sources by fully deterministic periodic time-shifts. In this paper we give a mathematical proof showing that the signal apparition method results in optimally large regions in the frequency-wavenumber space where exact separation of sources is achieved. These regions are diamond-shaped and we prove that using any other method of source encoding results in strictly smaller regions of exact separation. The results are valid for arbitrary number of sources. Numerical examples for different number of sources (three resp.~four sources) demonstrate the exact recovery of these diamond-shaped regions. The theoretical proofs' implementation in the field is illustrated by the results of a conducted field test., Comment: 16 pages, 5 figures. Expanded Section 3 with an additional numerical experiment and the results of a field test. Added references

We analyse a linear lattice Boltzmann (LB) formulation for simulation of linear acoustic wave propagation in heterogeneous media. We employ the single-relaxation-time Bhatnagar-Gross-Krook (BGK) as well as the general multi-relaxation-time (MRT) collision operators. By calculating the dispersion relation for various 2D lattices, we show that the D2Q5 lattice is the most suitable model for the linear acoustic problem. We also implement a grid-refinement algorithm for the LB scheme to simulate waves propagating in a heterogeneous medium with velocity contrasts. Our results show that the LB scheme performance is comparable to the classical second-order finite-difference schemes. Given its efficiency for parallel computation, the LB method can be a cost effective tool for the simulation of linear acoustic waves in complex geometries and multiphase media.

UNKNOWN