Array beamforming With a seismic array, the signal-to-noise ratio (SNR) of a seismic signal can be improved by summing the coherent signals from the individual array sites. The most important point during the
beamforming process is to find the best delay times by which the single traces must be shifted before summation in order to get the largest
amplitudes due to coherent interference of the signals. s.
Weighted stack methods Schimmel and Paulssen introduced another non-linear stacking technique in 1997 to enhance signals through the reduction of incoherent noise, which shows a smaller waveform distortion than the N-th root process. Kennett proposed the use of the
semblance of the signal as a weighting function in 2000 and achieved a similar resolution. An easily implementable weighted stack method would be to weight the amplitudes of the single sites of an array with the SNR of the signal at this site before beamforming, but this does not directly exploit the coherency of the signals across the array. All weighted stack methods can increase the slowness resolution of velocity
spectrum analysis.
Double beam technique A cluster of earthquakes can be used as a source array to analyze coherent signals in the seismic coda. This idea was consequently expanded by Krüger et al. in 1993 by analyzing seismic array data from well-known source locations with the so-called "double beam method". The principle of reciprocity is used for source and receiver arrays to further increase the resolution and the SNR for small amplitude signals by combining both arrays in a single analysis.
Array transfer function The array transfer function describes sensitivity and resolution of an array for seismic signals with different frequency contents and slownesses. With an array, we are able to observe the wavenumber k=2\pi/\lambda=2\pi\cdot f\cdot s of this wave defined by its frequency f and its slowness s. While time-domain
analog-to-digital conversion may give aliasing effects in the time domain, the spatial sampling may give aliasing effects in the
wavenumber domain. Thus the wavelength range of seismic signals and the sensitivity at different wavelengths must be estimated. and further developed to include wide-band analysis, maximum-likelihood estimation techniques, and three-component data in the 1980s. The methodology exploits the deterministic, non-periodic character of seismic wave propagation to calculate the frequency-wavenumber spectrum of the signals by applying the
multidimensional Fourier transform. A monochromatic
plane wave w(x,t) will propagate along the x direction according to equation :w(x,t)=Ae^{i2\pi(f_{0}t-k_{0}x)} It can be rewritten in frequency domain as :W(k_{x},f)=A\delta(f-f_{0})\delta(k_{x}-k_{0}) which suggests the possibility to map a monochromatic plane wave in the frequency-wavenumber domain to a point with coordinates (f, kx) = (f0, k0). Practically, f-k analysis is performed in the frequency domain and represents in principle beamforming in the frequency domain for a number of different slowness values. At
NORSAR slowness values between -0.4 and 0.4 s/km are used equally spaced over 51 by 51 points. For every one of these points the beam power is evaluated, giving an equally spaced grid of 2601 points with power information. Because of the amount of required computations, plane wave fitting is most effective for arrays with a smaller number of sites or for subarray configurations. == Applications ==