Interferometry by deconvolution: Part 1 — Theory for acoustic waves and numerical examples

Geophysics ◽  
2008 ◽  
Vol 73 (3) ◽  
pp. S115-S128 ◽  
Author(s):  
Ivan Vasconcelos ◽  
Roel Snieder
Keyword(s):  
Boundary Condition ◽  
Acoustic Waves ◽  
Time Function ◽  
Source Time Function ◽  
Field Quantity ◽  
Representation Theorems ◽  
Direct Wave ◽  
Wavefield Separation ◽  
Acoustic Media ◽  
Deconvolution Interferometry

Interferometry allows for synthesis of data recorded at any two receivers into waves that propagate between these receivers as if one of them behaves as a source. This is accomplished typically by crosscorrelations. Based on perturbation theory and representation theorems, we show that interferometry also can be done by deconvolutions for arbitrary media and multidimensional experiments. This is important for interferometry applications in which (1) excitation is a complicated source-time function and/or (2) when wavefield separation methods are used along with interferometry to retrieve specific arrivals. Unlike using crosscorrelations, this method yields only causal scattered waves that propagate between the receivers. We offer a physical interpretation of deconvolution interferometry based on scattering theory. Here we show that deconvolution interferometry in acoustic media imposes an extra boundary condition, which we refer to as the free-point or clamped-point boundary condition, depending on the measured field quantity. This boundary condition generates so-called free-point scattering interactions, which are described in detail. The extra boundary condition and its associated artifacts can be circumvented by separating the reference waves from scattered wavefields prior to interferometry. Three wavefield-separation methods that can be used in interferometry are direct-wave interferometry, dual-field interferometry, and shot-domain separation. Each has different objectives and requirements.

scholarly journals Fully probabilistic seismic source inversion – Part 1: Efficient parameterisation

2013 ◽  
Vol 5 (2) ◽  
pp. 1125-1162 ◽  
Author(s):  
S. C. Stähler ◽  
K. Sigloch
Keyword(s):  
Moment Tensor ◽  
Source Parameters ◽  
Practical Importance ◽  
Principal Component ◽  
Seismic Source ◽  
Empirical Orthogonal Functions ◽  
Time Function ◽  
Source Inversion ◽  
Source Time Function ◽  
Earthquake Parameters

Abstract. Seismic source inversion is a non-linear problem in seismology where not just the earthquake parameters themselves, but also estimates of their uncertainties are of great practical importance. Probabilistic source inversion (Bayesian inference) is very adapted to this challenge, provided that the parameter space can be chosen small enough to make Bayesian sampling computationally feasible. We propose a framework for PRobabilistic Inference of Source Mechanisms (PRISM) that parameterises and samples earthquake depth, moment tensor, and source time function efficiently by using information from previous non-Bayesian inversions. The source time function is expressed as a weighted sum of a small number of empirical orthogonal functions, which were derived from a catalogue of >1000 STFs by a principal component analysis. We use a likelihood model based on the cross-correlation misfit between observed and predicted waveforms. The resulting ensemble of solutions provides full uncertainty and covariance information for the source parameters, and permits to propagate these source uncertainties into travel time estimates used for seismic tomography. The computational effort is such that routine, global estimation of earthquake mechanisms and source time functions from teleseismic broadband waveforms is feasible.

The 14 September 1995 (M = 7.3) Copala, Mexico, earthquake: A source study using teleseismic, regional, and local data

1997 ◽  
Vol 87 (4) ◽  
pp. 999-1010
Author(s):  
F. Courboulex ◽  
M. A. Santoyo ◽  
J. F. Pacheco ◽  
S. K. Singh
Keyword(s):  
Time Window ◽  
Time Function ◽  
Source Process ◽  
Source Time Function ◽  
Nonlinear Method ◽  
Local Data ◽  
Plate Interface ◽  
Linear Inversion ◽  
Field Source ◽  
Mexico Earthquake

Abstract We analyze source characteristics of the 14 September 1995, Copala, Mexico, earthquake (M = 7.3) using teleseismic, regional, and local seismograms. In the analysis of the teleseismic and the regional seismic waves, we apply the empirical Green's function (EGF) technique. The recording of an appropriate aftershock is taken as the EGF and is used to deconvolve the mainshock seismogram, thus obtaining an apparent far-field source-time function at each station. The deconvolution has been done using surface waves. For teleseismic data, we apply a spectral deconvolution method that enables us to obtain 37 apparent source-time functions (ASTFs) at 29 stations. In the analysis of the regional broadband seismograms, we use two different aftershocks as EGF, and the deconvolution is performed in the time domain with a nonlinear method, imposing a positivity constraint, and the best azimuth for the directivity vector is obtained through a grid-search approach. We also analyze two near-source accelerograms. The traces are inverted for the slip distribution over the fault plane by applying a linear inversion technique. With the aid of a time-window analysis, we obtain an independent estimation of the source-time function and a more detailed description of the source process. The analysis of the three datasets permits us to deduce the main characteristics of the source process. The rupture initiated at a depth of 16 km and propagated in two directions: updip along the plate interface toward 165° N and toward 70° N. The source duration was between 12 and 14 sec, with the maximum of energy release occurring 8 sec after the initiation of the rupture. The estimated rupture dimension of 35 × 45 km is about one-fourth of the aftershock area. The average dislocation over the fault was 1.4 m (with a maximum dislocation of 4.1 m located 10 km south of the hypocenter), which gives roughly 1 MPa as the average static stress drop.

On sound energy scattered by a rigid body near a compliant surface

1995 ◽  
Vol 451 (1943) ◽  
pp. 543-552 ◽  
Keyword(s):  
Acoustic Waves ◽  
Specular Reflection ◽  
Simple Relation ◽  
Optical Theorem ◽  
Energy Calculation ◽  
Far Field ◽  
Sound Energy ◽  
Compliant Surface ◽  
Scattered Energy ◽  
Acoustic Media

A theorem is derived which generalizes the classical extinction theorem (also known as the optical theorem) to cases where a rigid scatterer of arbitrary shape is lo­cated near a large compliant surface which has quite general mechanical properties, including dissipation and wave-bearing features, and where the acoustic media on both sides of the compliant surface may have different densities and sound speeds. The theorem relates the sound energy scattered from incident planar acoustic waves to the far field pressure in the specular reflection direction and that in the transmis­sion direction, determined by the Snell’s law. From this simple relation, the scattered energy can be found almost trivially from the far field pressures in these two partic­ular directions; the energy calculation then completely avoids integration of energy flux over control surfaces.

scholarly journals Flow in the exit of open pipes during acoustic resonance

1980 ◽  
Vol 99 (2) ◽  
pp. 293-319 ◽  
Author(s):  
J. H. M. Disselhorst ◽  
L. Van Wijngaarden
Keyword(s):  
Boundary Layer ◽  
Mathematical Model ◽  
Boundary Condition ◽  
Acoustic Wave ◽  
Acoustic Waves ◽  
Acoustic Radiation ◽  
Acoustic Resonance ◽  
The Other ◽  
Open Tube ◽  
Theoretical Predictions

The flow near the mouth of an open tube is examined, experimentally and theoretically, under conditions in which resonant acoustic waves are excited in the tube at the other end. If the edge of the tube is round, separation does not occur at high Strouhal numbers, which enables us to verify theoretical predictions for dissipation in the boundary layer and for acoustic radiation. Observation with the aid of schlieren pictures shows that in the case of a sharp edge vortices are formed during inflow. The vortices are shed from the pipe during outflow. Based on these observations a mathematical model is developed for the generation and shedding of vorticity. The main result of the analysis is a boundary condition for the pressure in the wave, to be applied near the mouth. The pressure amplitudes in the acoustic wave measured under resonance are compared with theoretical predictions made with the aid of the boundary condition obtained in the paper.

Sign in / Sign up

Export Citation Format

Share Document