Empirical Mode Decomposition (EMD) of potential field data: airborne gravity data as an example
Hassan H. Hassan and John W. Peirce
GEDCO, Calgary, Alberta, Canada
In this paper, we introduce a newly developed technique to process potential field data as an alternative to Fourier and wavelet based techniques. This new technique is called the Empirical Mode Decomposition (EMD) and was recently developed by Dr. Norden E. Huang at the NASA Goddard Space Flight Center (Huang et al., 1998). The EMD method is different from the Fourier and wavelet transforms because it handles nonlinear and non-stationary signals.
The Fourier transform (FFT) is designed to work with linear and stationary signals. The wavelet transform, on the other hand, is well-suited to handle non-stationary data, but it is poor at processing nonlinear data. Additionally, the basis functions used in FFT and, to some extent, wavelet trans- forms are fixed and do not necessarily match the varying nature of signals and this will lead to the loss of some useful information in the signal. Since potential field data a re in general nonlinear and non-stationary in nature, we expect limitations in processing the data using FFT or wavelet methods.
This work applies the EMD technique to process potential field data using airborne gravity over the Turner Valley area in the foothills of Alberta, Canada (Figure 1) in order to improve noise removal and thereby enhance the gravity signal.
Airborne gravity surveys are becoming more common in the oil and mining industries. However, the quality of the results is limited by the level of noise introduced into the data during acquisition. We lack processing techniques that effectively remove noise from airborne gravity data. Traditionally, FFT and more recently, wavelet transform have been used to sepa- rate the noise from the signal. However, these data are non- stationary and nonlinear, and therefore neither FFT nor wavelet transform are really appropriate for this application. We have tested the EMD technique because we feel that EMD is more appropriate to process this type of data.
The data selected to test the EMD method is derived from an AIRGrav survey flown over the Turner Valley region of Alberta in the summer of 2001 by Sander Geophysics Ltd. (Peirce et al., 2002). The data set consists of over 12,000 line km of airborne gravity data flown on 250 m spaced east-west traverse lines and 1000 m spaced north-south control lines. The survey was flown with drape elevation that varies from
250 m height in the plains to over 500 m height in the moun- tains. After full processing on a line by line basis the data were leveled and the standard Bouguer reduction correc- tions, including outer terrain corrections using a reduction
density of 2.67 gm/cc, were applied. Figure 1 shows the first vertical derivative of the filtered complete Bouguer anomaly.
The area selected for this study is located within the red box of Figure 2. The study area flanks the eastern edge of the Rocky Mountain where it is dominated by north-south trending faults associated with the foothills region (Fig. 2). The eastern side of the area consists of flat lying sediments. The Turner Valley region in general is a well-established area for oil and gas production and was the site of the first Alberta oil boom in the 1920’s. New discoveries are still being drilled in the structure and in sub-thrust plays where accurate depth mapping from seismic data is a challenge. The principal p roducing zones are from porous Mississippian aged carbonate rocks carried in overthrust structures.
The EMD technique is part of a process known as the Hilbert–Huang Transform (HHT) that consists of two main elements: the EMD and the Hilbert spectral analysis. The EMD generates the intrinsic mode functions (IMFs) from the data, and the Hilbert spectral analysis generates a “time- frequency-energy” representation of the data, based on the
Valley airborne gravity survey draped on NE-shaded topography.
Continued on Page 26
IMFs. In this study we are only concerned with the EMD part of
Hilbert–Huang Transform (HHT).
The EMD is an adaptive decomposition technique with which any complicated signal can be decomposed into a definite number of high-frequency and low frequency components by means of a process called “sifting”.
The sifting process decomposes the original signal, S(x), into a number of intrinsic mode functions (IMFs) according to the following formula:
Figure 2. Index map showing the test area (red block) and the coincident ground and airborne test lines (blue) used in this study.
Figure 3. EMD decomposition of two waves signal, (a) sum of two waves, (b) lower and upper envelopes (red) and their mean (blue), (c) the first IMF and (d) the first residual (after Oonincx and Hermand, 2004).
Figure 4. Synthetic geological model and its magnetic response.
Figure 5. EMD decomposition of the magnetic response (with added random noise)
of the synthetic geological model shown in Figure 4.
Continued on Page 27
Empirical Mode Decomposition…
S(x) = rn(x) + ci(x)
where, ci(x) represents an IMF, and rn(x), is the residual after the
n IMFs have been extracted.
Figure 6. EMD decomposition along a test coincides with the ground gravity data
(see Figure 2).
These IMFs have well-behaved Hilbert transforms and are defined as functions that:
(1) have the same number of zero-crossings and extrema, and
(2) the mean value of the upper and the lower envelopes is equal to zero.
A sifting process extracts IMFs from the signal iteratively in order to obtain a component that satisfies above mentioned conditions. The sifting process separates the IMFs with decreasing order of frequency, i.e., it separates high frequency component first and the low frequency component at the end.
The EMD technique (Huang et al., 1998) is illustrated in Figure 3 for a simple signal consisting of two superposed sine waves. The decomposition of the signal into IMFs is performed as follows:
1. Identify the positive peaks (maxima) and negative peaks
(minima) of the original signal.
2. Construct the lower and the upper envelopes of the signal by the cubic spline method (red).
3. Calculate the mean values (blue) by averaging the upper enve- lope and the lower envelope.
4. Subtract the mean from the original signal to produce the first intrinsic mode function IMF1 component.
5. Calculate the first residual component by subtracting IMF1 from the original signal. This residual component is treated as a new data and subjected to the same process described above to calculate the next IMF.
6. Repeat the steps above until the final residual component becomes a monotonic function and no more IMFs can be extracted.
Figure 3 clearly illustrates that by using EMD we have success- fully reconstructed the two superposed sine waves in the orig- inal signal.
Generally speaking, the sifting process produces a set of IMFs that represent the original data vector broken down into frequency components from highest to lowest frequency. If all of the IMFs for a given signal are added together, the resulting “summation” signal is a near perfect match for the original signal (i.e., with little or no leftover), yielding a high level of confidence in the EMD results.
We have used a Matlab code written by Dr. Patrick Flandrin of Centre National De Recherche Scientifique (CNRS) in Lyon, France (Flandrin et al., 2004) to compute the EMD of the Turner Valley airborne gravity data.
In order to examine the efficiency of the EMD technique to sepa- rate noise from the data and to recover the original signal in the data we have performed a simple test. We have constructed a hypothetical geological model (Fig. 4) and calculated its magnetic response. We have used magnetic data instead of gravity data for demonstration purpose because magnetic
Continued on Page 28
response is sharper than gravity response of the same model. We have added random noise to the generated total magnetic field from the model to a degree that the noise completely masked the original signal (Fig. 5, second panel from the top). Then we ran EMD on the noisy profile and plotted the results. EMD was able to recover (quite well) the original signal from the noisy data (i.e., IMF4 of Fig. 5).
We have selected a segment of a ground gravity line that coin- cides with a portion of Line 902500 of the airborne gravity survey (blue line in Fig. 2) in order to test EMD decomposition. The results are shown in Figure 6. Figure 6 also shows the complete Bouguer gravity of the ground gravity line as well as the unfil- tered and filtered (1500 m and 5300 m lowpass) profiles of the airborne gravity for comparison. The results show that the main components of the ground gravity profile are contained in IMF5, with perhaps some additional contribution from IMF4. These results provide us with some confidence in the technique. Figure
6 also shows that most of the high frequency components, including noise, are contained in IMF1, IMF2 and to some extent IMF3 whereas the low frequency components are confined to IMF4 and IMF5.
Then we used the EMD method to process the line data of the airborne gravity survey that covers the red box area shown on Figure 2. Based on our results of EMD decomposition of the test
lines (Fig. 6) it appears that the useful geological information are probably contained in IMF4 and IMF5. For this reason, we have gridded and contoured the IMF4 and IMF5 components (Figs. 7 and 8, respectively) and compared the results with the 5300 m lowpass (Fig. 9) and the first vertical derivative (Fig. 10) of the original dataset. The results clearly indicate that the EMD decomposition is more resolving and reveal anomalies that coin- cide with the prominent geological structures in the area.
A newly developed technique called EMD was tested on airborne gravity data from the Turner Valley airborne gravity survey in west-central Alberta as an alternative to standard FFT and wavelet analysis techniques. This new technique is powerful at analyzing nonlinear and non-stationary signals such as airborne gravity data. It decomposes the signal to a summation of zero-mean AM-FM components, called Intrinsic Mode Functions (IMF). These IMFs show the main components of the analyzed signal.
The initial results of our test is encouraging and show that there is some application potential in isolating noise from airborne gravity data and to detect meaningful geological information that might have been masked by the amount of noise in the data. This technique can be used as an alternative to lowpass filtering of airborne gravity data because it appears to better preserves
anomaly amplitude and wavelength. R
Figure 7. Map equivalent to IMF4 in Figure 6.
Figure 8. Map equivalent to IMF5 in Figure 6.
Continued on Page 30
Focus Article Cont’d Empirical Mode Decomposition… Continued from Page 28
Flandrin, P., Rilling, G., and Goncalves, P., 2004, ^
Filter Bank, IEEE Signal Processing Letters, 112 – 114.
Huang, N.E., Shen, Z., Long, S.R., Wu, M.C., Shih, H.H., Zheng, Q., Yen, N., Tung, C.C., and Liu, H.H., 1998, The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis, Royal Society London, 903-995.
Oonincx, P.J., and Hermand, J.P., 2004, ^ Proceedings of the Seventh European Conference on Underwater Acoustics, Delft.
Peirce, J.W., Sander, S., Charters, R.A., and Lavoie, V., 2002, Turner Valley, Canada
– A case history in contemporary airborne gravity: Society of Exploration Geophysicists
International Exposition and 72nd Annual Meeting, 783 – 786.
Figure 9. 5300 m lowpass filter of complete Bouguer gravity draped on NE-shaded Figure 10. First vertical derivative of 5300 m lowpass filter of complete Bouguer
topography. gravity draped on NE-shaded topography.
Hassan H. Hassan received his B.Sc. and M.Sc. degrees in geology and geophysics, respectively from the University of Baghdad, Iraq. Following a scholarship from the Iraqi Atomic Energy, Hassan continued his studies in Canada. He obtained a Ph.D. degree (1984) in geology/geophysics from the University of New Brunswick. In 1984 he worked as a research scien- tist in the Department of Surveying Engineering (Geomatics) at the University of New Brunswick on remote sensing data processing and application in geology. In 1985 Hassan joined the Department of Geology at the University of New Brunswick as a postdoctoral research associate to work on a five year project sponsored by the Geological Survey of Canada. Prior to joining GEDCO in July 1996 Hassan worked as a contract geophysicist on several geophysical projects in New Brunswick for the NB Department of Natural Resources & Energy and for the Atlantic Geoscience Center. Currently, Hassan is working as a senior geophysicist with the potential field group at GEDCO. Hassan is a member of CSEG, SEG and APEGGA.
^ is one of GEDCO’s founding partners. He has wide-ranging interpretational experience in integrated basin analysis and the evaluation of new acreage blocks in both Canadian and international settings. He specializes in high resolu- tion aeromagnetics, airborne gravity, and the integration of gravity/magnetics and seismic data into one combined interpre- tation. John’s international experience includes projects in Mexico, South America, Africa, Middle East, Pakistan, Russia, China, Southeast Asia and Indonesia. He is a co-author of the SEG book titled “Planning Land 3D Seismic Surveys”, published in 2000. As of October, 2007, John is “mostly retired” and continues as Chairman of GEDCO, working on strategic issues and research topics of interest.
He is an Associate Editor for Gravity for Geophysics, and a member of APEGGA, APEGBC, CSEG, SEG, CSPG, AAPG and the WHOI Corporation (an advisory body for Woods Hole Oceanographic Institution). John is also a Board member of the Gabriola Land Conservancy in British Columbia and runs the photofinish camera for the Cross Country World Cup Society in Canmore and at other venues in B.C.
30 CSEG RECORDER January 2008
Оригинал: Empirical Mode Decomposition (emd) of potential field data: airborne gravity data as an example. Hassan H. Hassan and John...
1. /John-Lennon/Джон Леннон.doc
1. /John Barth - Frame Tale.pdf
2. /John Barth...
|I. compositions and pieces john Abercrombie||Angel pavement by John B. Priestley|
|John Constable Self-portrait 1806||Dangerous corner by John Boynton Priestley|
|John Sargent was given little regular schooling. As a result of his "Baedeker education," he learned Italian, French, and German||Документы|
1. /John Barth - Chimera.doc