 | This article may be too technical for most readers to understand.(November 2015) |
A velocity filter removes interfering signals by exploiting the difference between the travelling velocities of desired seismic waveform and undesired interfering signals.
In geophysical applications sensors are used to measure and record the seismic signals. [1] [2] Many filtering techniques are available in which one output waveform is produced with a higher signal-to-noise ratio than the individual sensor recordings. Velocity filters are designed to remove interfering signals by exploiting the difference between the travelling velocities of desired seismic waveform and undesired interfering signals. [3] In contrast to the one dimensional output produced by multi-channel filtering, velocity filters produce a two-dimensional output.
Consider an array of N sensors that receive one desired and M undesired broadband interferences. Let the measurement of nth sensor be modeled by the expression:
where n=1,2,...,N; m=0,1,...,M; sm(t) are signals travelling across the array, and ŋn(t) represents zero-mean white random noise at the nth sensor, uncorrelated from sensor to sensor. The parameters amn and Tmn are amplitude gain and time delays of the signal sm(t) when received at nth sensor.
Without loss of generality, we shall assume that s0(t) is the desired signal and s1(t),s2(t),...,sM(t) are the undesired interferences. Additionally we shall assume that T0n=0, and a0n =1. This essentially means that the data has been time shifted to align the desired seismic signal so that it appears on all sensors at the same time and balanced so that the desired signal appears with equal amplitudes. We assume that the signals are digitized prior to being recorded and that the length K of time sequences of recorded data is large enough for the complete delayed interfering waveforms to be included in the recorded data. In the discrete frequency domain, the nth trace can be expressed as:
where k=0,1,...,K-1; wk=(2π/K) is the sampling angular frequency.
Using matrix notation, (2) can be expressed in the form:
Frequency domain multichannel filters F1(k), F2(k), ..., FN(k) can be applied to the data to produce one single output trace of the form,: [4] [5]
In matrix form, the above expression can be written as:
where F(k) is an N x 1 vector whose elements are the individual channel filters. That is,
By following the procedure discussed in, [6] [7] an optimum filter vector F(k) can be designed to attenuate, in the least square sense, the undesired coherent interferences S1(k),S2(k),...,SM(k) while preserving the desired signal S0(k) in Y(k). This filter can be shown [6],[7] to be of the form:
where h is an arbitrary N x 1 nonzero vector, u = [1,0,...,0], I is the unit matrix, Br(k) is a submatrix of the matrix obtained by dropping all linearly dependent rows, and L(k) is a lower triangular matrix satisfying:
The multichannel processing scheme described by [6]-[11] produces one dimensional output trace. A velocity filter, on the other hand, is a two-dimensional filter which produces a two-dimensional output record.
A two-dimensional record can be generated by a procedure which involves repeatedly applying multichannel optimum filters to a small number of overlapping subarrays of the input data,. [8] [9]
More specifically, consider a subarray of W channels, where W<<N, which slides over the input data as shown in Fig. 1. For every subarray position an optimum multichannel filter based on (9) can be designed so that the undesired interferences are suppressed from its corresponding output trace. In designing this filter we use W instead of N in expression (9). thus traces 1,2,...,W of the input record produce the first trace of the output record, traces K,K+1,...K+W-1 of the input record produce the Kth trace of the output record, and traces N-W+1,N-W+2,...,N of the input record produce the (N-W+1)st trace, which is the last trace, of the output record. For a large N and small W, as is typically the case in geophysical data, the output record can be viewed as comparable in dimensions to the input record. Clearly for such a scheme to work effectively W must be as small as possible; while at the same time it must be large enough to provide the necessary attenuation of the undesired signals. Note that a maximum of W-1 undesired interferences can be totally suppressed by such a scheme. [10] [11]
In computer science, counting sort is an algorithm for sorting a collection of objects according to keys that are small positive integers; that is, it is an integer sorting algorithm. It operates by counting the number of objects that possess distinct key values, and applying prefix sum on those counts to determine the positions of each key value in the output sequence. Its running time is linear in the number of items and the difference between the maximum key value and the minimum key value, so it is only suitable for direct use in situations where the variation in keys is not significantly greater than the number of items. It is often used as a subroutine in radix sort, another sorting algorithm, which can handle larger keys more efficiently.
An adaptive filter is a system with a linear filter that has a transfer function controlled by variable parameters and a means to adjust those parameters according to an optimization algorithm. Because of the complexity of the optimization algorithms, almost all adaptive filters are digital filters. Adaptive filters are required for some applications because some parameters of the desired processing operation are not known in advance or are changing. The closed loop adaptive filter uses feedback in the form of an error signal to refine its transfer function.
A sensor array is a group of sensors, usually deployed in a certain geometry pattern, used for collecting and processing electromagnetic or acoustic signals. The advantage of using a sensor array over using a single sensor lies in the fact that an array adds new dimensions to the observation, helping to estimate more parameters and improve the estimation performance. For example an array of radio antenna elements used for beamforming can increase antenna gain in the direction of the signal while decreasing the gain in other directions, i.e., increasing signal-to-noise ratio (SNR) by amplifying the signal coherently. Another example of sensor array application is to estimate the direction of arrival of impinging electromagnetic waves. The related processing method is called array signal processing. A third examples includes chemical sensor arrays, which utilize multiple chemical sensors for fingerprint detection in complex mixtures or sensing environments. Application examples of array signal processing include radar/sonar, wireless communications, seismology, machine condition monitoring, astronomical observations fault diagnosis, etc.
Direct digital synthesis (DDS) is a method employed by frequency synthesizers used for creating arbitrary waveforms from a single, fixed-frequency reference clock. DDS is used in applications such as signal generation, local oscillators in communication systems, function generators, mixers, modulators, sound synthesizers and as part of a digital phase-locked loop.
Array processing is a wide area of research in the field of signal processing that extends from the simplest form of 1 dimensional line arrays to 2 and 3 dimensional array geometries. Array structure can be defined as a set of sensors that are spatially separated, e.g. radio antenna and seismic arrays. The sensors used for a specific problem may vary widely, for example microphones, accelerometers and telescopes. However, many similarities exist, the most fundamental of which may be an assumption of wave propagation. Wave propagation means there is a systemic relationship between the signal received on spatially separated sensors. By creating a physical model of the wave propagation, or in machine learning applications a training data set, the relationships between the signals received on spatially separated sensors can be leveraged for many applications.
Beamforming or spatial filtering is a signal processing technique used in sensor arrays for directional signal transmission or reception. This is achieved by combining elements in an antenna array in such a way that signals at particular angles experience constructive interference while others experience destructive interference. Beamforming can be used at both the transmitting and receiving ends in order to achieve spatial selectivity. The improvement compared with omnidirectional reception/transmission is known as the directivity of the array.
In signal processing, a filter bank is an array of bandpass filters that separates the input signal into multiple components, each one carrying a single frequency sub-band of the original signal. One application of a filter bank is a graphic equalizer, which can attenuate the components differently and recombine them into a modified version of the original signal. The process of decomposition performed by the filter bank is called analysis ; the output of analysis is referred to as a subband signal with as many subbands as there are filters in the filter bank. The reconstruction process is called synthesis, meaning reconstitution of a complete signal resulting from the filtering process.
Smart antennas are antenna arrays with smart signal processing algorithms used to identify spatial signal signatures such as the direction of arrival (DOA) of the signal, and use them to calculate beamforming vectors which are used to track and locate the antenna beam on the mobile/target. Smart antennas should not be confused with reconfigurable antennas, which have similar capabilities but are single element antennas and not antenna arrays.
Seismic Handler (SH) is an interactive analysis program for preferably continuous waveform data. It was developed at the Seismological Observatory Gräfenberg and is in use there for daily routine analysis of local and global seismic events. In original form Seismic Handler was command line based, but now an interactive version is available.
Space-time adaptive processing (STAP) is a signal processing technique most commonly used in radar systems. It involves adaptive array processing algorithms to aid in target detection. Radar signal processing benefits from STAP in areas where interference is a problem. Through careful application of STAP, it is possible to achieve order-of-magnitude sensitivity improvements in target detection.
Geophysical survey is the systematic collection of geophysical data for spatial studies. Detection and analysis of the geophysical signals forms the core of Geophysical signal processing. The magnetic and gravitational fields emanating from the Earth's interior hold essential information concerning seismic activities and the internal structure. Hence, detection and analysis of the electric and Magnetic fields is very crucial. As the Electromagnetic and gravitational waves are multi-dimensional signals, all the 1-D transformation techniques can be extended for the analysis of these signals as well. Hence this article also discusses multi-dimensional signal processing techniques.
Underwater acoustic communication is a technique of sending and receiving messages below water. There are several ways of employing such communication but the most common is by using hydrophones. Underwater communication is difficult due to factors such as multi-path propagation, time variations of the channel, small available bandwidth and strong signal attenuation, especially over long ranges. Compared to terrestrial communication, underwater communication has low data rates because it uses acoustic waves instead of electromagnetic waves.
In geophysics, seismic inversion is the process of transforming seismic reflection data into a quantitative rock-property description of a reservoir. Seismic inversion may be pre- or post-stack, deterministic, random or geostatistical; it typically includes other reservoir measurements such as well logs and cores.
3D sound localization refers to an acoustic technology that is used to locate the source of a sound in a three-dimensional space. The source location is usually determined by the direction of the incoming sound waves and the distance between the source and sensors. It involves the structure arrangement design of the sensors and signal processing techniques.
Sonar systems are generally used underwater for range finding and detection. Active sonar emits an acoustic signal, or pulse of sound, into the water. The sound bounces off the target object and returns an “echo” to the sonar transducer. Unlike active sonar, passive sonar does not emit its own signal, which is an advantage for military vessels. But passive sonar cannot measure the range of an object unless it is used in conjunction with other passive listening devices. Multiple passive sonar devices must be used for triangulation of a sound source. No matter whether active sonar or passive sonar, the information included in the reflected signal can not be used without technical signal processing. To extract the useful information from the mixed signal, some steps are taken to transfer the raw acoustic data.
This article provides a short survey of the concepts, principles and applications of Multirate Filter Banks and Multidimensional Directional Filter Banks.
Geophysical signal analysis is concerned with the detection and a subsequent processing of signals. Any signal which is varying conveys valuable information. Hence to understand the information embedded in such signals, we need to 'detect' and 'extract data' from such quantities. Geophysical signals are of extreme importance to us as they are information bearing signals which carry data related to petroleum deposits beneath the surface and seismic data. Analysis of geophysical signals also offers us a qualitative insight into the possibility of occurrence of a natural calamity such as earthquakes or volcanic eruptions.
A seismic array is a system of linked seismometers arranged in a regular geometric pattern to increase sensitivity to earthquake and explosion detection. A seismic array differs from a local network of seismic stations mainly by the techniques used for data analysis. The data from a seismic array is obtained using special digital signal processing techniques such as beamforming, which suppress noises and thus enhance the signal-to-noise ratio (SNR).
Multidimensional seismic data processing forms a major component of seismic profiling, a technique used in geophysical exploration. The technique itself has various applications, including mapping ocean floors, determining the structure of sediments, mapping subsurface currents and hydrocarbon exploration. Since geophysical data obtained in such techniques is a function of both space and time, multidimensional signal processing techniques may be better suited for processing such data.
High Resolution Wide Swath (HRWS) imaging is an important branch in Synthetic aperture radar (SAR) imaging, a remote sensing technique capable of providing high resolution images independent of weather conditions and sunlight illumination. This makes SAR very attractive for the systematic observation of dynamic processes on the Earth's surface, which is useful for environmental monitoring, earth resource mapping and military systems.
