Photovoltaic Cells and Systems: 45-56 SQU Journal for Science, 16 (2011) © 2011 Sultan Qaboos University 45 Spectral Analysis of Magnetic Anomalies Due to a 2-D Horizontal Circular Cylinder: A Hartley Transforms Technique Mansour A. Al-Garni Department of Geophysics, Faculty of Earth Sciences, King Abdulaziz University, Jeddah, Saudi Arabia, Email: Maalgarni@kau.edu.sa. التحليل الطيفي للشواذ المغناطيسية الناتجة من اسطوانة أفقية ذو مقطع دائري ثنائية األبعاد: تقنية محوالت هارتلي عبدهللا القرنيمنصور الناتجة من اسطوانة أفقٌة ذو مقطع دائري ثنائٌة المغناطٌسٌةلشواذ العمودي تخدام محول هارتلً للتحلٌللقد تم اس :صخمل لقد تم حساب تبر محول هارتلً كوسٌلة بدٌلة للتحلٌل الطٌفً على غرار التحلٌل الطٌفً باستخدام محول فورٌر.ٌع األبعاد. ع كما انه تم استخدام مثال مصطن العمق الى مركز األسطوانة األفقٌة باستخدام معادلة رٌاضٌة بسٌطة كدالة فً التردد. الطرٌقة المقترحة بنجاح على شاذة هنه قد تم تطبٌق هذإى ذلك ففة الضا. باإلومدى صالحٌتها التقنٌة هلتوضٌح خطوات هذ تأثٌرمن بلدة كارٌمناغار، الهند. لقد تمت دراسة تٌت مأخوذة من مانغامبالً بالقربماغنٌ -حقلٌة على شرٌط من الكوارتز لطرٌقة المقترحة على مثال الشاذة ن نتائج اأكما .لثقةظهرت مستوى عال من اأو المقترحةالشوشرة العشوائٌة على الطرٌقة الحقلٌة أظهرت تطابقا مع نتائج الطرق األخرى المنشورة. ABSTRACT: The spectral analysis of the vertical effect of magnetic anomalies due to a 2-D horizontal circular cylinder is presented using Hartley transform. Hartley transform is an alternative approach to the famous complex Fourier transform. The depth to the center of the horizontal cylinder can be computed by a simple equation as a function of frequency. A synthetic example has been used to illustrate this technique and the validity of this approach has been proved by applying it to real data of a narrow band of quartz-magnetite in Mangampalli near Karimnagar town, India. The noise analyses were tested on the technique and showed a high level of confidence. The results of the field example are in good agreement with the ones published in the literature. KEYWORDS: Hartley transform; Magnetic; 2-D Horizontal cylinder. 1. Introduction he Hartley transform (Hartley, 1942) has gained an importance in the field of geophysics in the last decade (Bracewell, 1983; Villasenor and Bracewell, 1987; Saatcilar et al. 1990; Saatcilar and Ergintav, 1991; T mailto:Maalgarni@kau.edu.sa MANSOUR A. AL-GARNI 68 Sundararajan 1995, 1997; Sundararajan et al. 2007). The importance of this transform has been ignored not because of the complexity of the transform but because scientists have been overwhelmed by the complex algebra concept (Sundararajan et al. 2007). The Hartley transform is purely real and exactly equivalent to the Fourier transform (Bracewell, 1983; Rajan, 1993; Sundararajan, 1995). The significance of this transform is that it requires no assumptions to be made, unlike the Fourier transform (Mohan and Seshagiri Rao, 1982). The Hartley and Fourier transforms provide two numbers, having the same information at each frequency, which represent a physical oscillation in amplitude and phase. Sundararajan and Brahmam (1998) used the Hartley transform to interpret gravity anomalies caused by slab-like structures. Sundararajan et al. (2007) used the Hartley transform to interpret the deformation of a homogeneous electric field over a thin bed. In this paper, the Hartley transform approach is used to estimate the causative target parameters of a 2-D horizontal circular cylinder from its magnetic anomaly. This approach is applied first to a theoretical example to illustrate the method and then applied to the vertical magnetic anomaly over a narrow band of quartz-magnetite in Mangampalli near Karimnagar town, India, to demonstrate the applicability of the method. Figure 1. Geometry of the 2-D horizontal circular cylinder. 2. The Magnetic effect of a horizontal cylinder The vertical magnetic effect of a buried horizontal circular cylinder extending infinitely in the horizontal direction along the Y-direction, with its normal section parallel to the x-z plane, is considered. The coordinate system origin is taken on the ground surface such that the Z-axis coincides with the diameter (Figure 1). Hence, the vertical magnetic effect at a point P on the surface can be expressed as follows (Mohan et al. 1982).       2 2 2 2 2 sin 2 cosh x x h V x K x h              (1) where 2 2 ,K R I R is the radius of the cylinder, I is the magnetization intensity, is the polarization angle, and h is the depth to the center of the cylinder. 3. Hartley transform of the magnetic effect The Hartley transform of the real function V(x) is defined by Hartley (1942) as: http://www.sciencedirect.com/science?_ob=ArticleURL&_udi=B6VFC-3V5M2BF-5&_user=1366817&_coverDate=05%2F31%2F1998&_rdoc=1&_fmt=full&_orig=search&_cdi=6007&_sort=d&_docanchor=&view=c&_acct=C000052384&_version=1&_urlVersion=0&_userid=1366817&md5=6ccad01ca7dfc2377ad689c7207313f3#b6 http://www.sciencedirect.com/science?_ob=ArticleURL&_udi=B6VFC-3V5M2BF-5&_user=1366817&_coverDate=05%2F31%2F1998&_rdoc=1&_fmt=full&_orig=search&_cdi=6007&_sort=d&_docanchor=&view=c&_acct=C000052384&_version=1&_urlVersion=0&_userid=1366817&md5=6ccad01ca7dfc2377ad689c7207313f3#b7 http://www.sciencedirect.com/science?_ob=ArticleURL&_udi=B6VFC-3V5M2BF-5&_user=1366817&_coverDate=05%2F31%2F1998&_rdoc=1&_fmt=full&_orig=search&_cdi=6007&_sort=d&_docanchor=&view=c&_acct=C000052384&_version=1&_urlVersion=0&_userid=1366817&md5=6ccad01ca7dfc2377ad689c7207313f3#b8 SPECTRAL ANALYSIS OF MAGNETIC ANOMALIES 69      H V x cas x dx      (2) where        cos sin 2 sin / 4cas x x x x        (3) is the kernel which represents a phase-shifted sine wave of 45° and hence takes the harmonics of both cosine and sine functions. The frequency ω has the same physical significance as in the case of the Fourier transform. By analogy with the real and imaginary components of the Fourier transform, the Hartley transform may be expressed in terms of its even and odd components as:      H E O    (4) where           cos 2 H H E V x x dx           (5) and           sin 2 H H o V x dx            (6) where      H V cas x        (7) and        cos sin 2 sin / 4cas x x x x          (8) Hence the amplitude and phase spectra can be expressed, respectively, as:       2 2 A E O    (9)       1 tan O E              (10) Numerically, the amplitude A(ω) is equivalent to the Fourier amplitude; however, the phase-shifted    differs by 45o from that of Fourier phase   F   . Therefore,   4 F      (11) Alternately, the amplitude A(ω) and phase-shifted    can be computed as:       2 2 2 H H A       (12)           1 tan H H H H                   (13) MANSOUR A. AL-GARNI 6: Substituting for V(x) in equation (1) into equation (2), the even and odd components of the Hartley transform for the vertical magnetic effect of the horizontal circular cylinder infinitely extending in the horizontal direction can be easily evaluated as:   sin h E K e       (14)   cos h O K e       (15) Therefore, the Hartley transform   H  (sum of the even and odd components), amplitude  A  and phase- shifted    of the horizontal circular cylinder infinitely extending in the horizontal direction can be given, respectively, as:    sin cos h H K e         (16)   h H A K e      (17) 2     (18) 4. Parameters evaluation At two successive frequencies i  and 1i   ,   i h H i i A K e      (19)   1 1 1 i h H i i A K e         (20) Where 2 / N x   is the fundamental frequency expressed in radian per unit length, N is the total number of samples and x is the station interval. At 1i  and dividing equation (19) by equation (20), one can obtain:       2 1 1 2 2 1 h A e A         (21) Taking the natural logarithm of both sides:     1 2 2 1 2 1 1 1n 1n A h A                (22) The term K is evaluated by substituting the value of h in equation (19)  H hA K e     (23) and  can be computed from equations (14 and 15) as:     1 tan E O            (24) SPECTRAL ANALYSIS OF MAGNETIC ANOMALIES 6; Therefore, based on equations (22) – (24), we can easily estimate the depth h of the polarization angle  and the magnetization intensity related parameter K of the buried cylinder. 5. Synthetic example The Hartley transform approach is illustrated by a synthetic model assuming a depth to the center of the horizontal circular cylinder 10h  units, a polarization angle 60  and 1K  unit (Figure 2). The even component, the odd component, the Hartley transform and the amplitude spectrum are computed and shown in Figures 3a, b, c and d, respectively. Using the method that has been developed throughout the text, the parameters ( K , h and  ) of the horizontal circular cylinder are estimated and the results are shown in Table 1. It can be noticed that the interpreted results, using the proposed technique, agree well with the assumed values. Figure 2. Response of vertical magnetic effect of a 2-D horizontal circular cylinder. 6. Noise analysis To investigate the noise effect on our estimation method, we added a synthetic anomaly with 5% and 10% of white Gaussian noise (WGN) as shown in Figures 4 and 5. The even components, the odd components, the Hartley transforms and the amplitude spectra of the contaminated anomalies are shown in Figures 6 and 7, respectively. The interpreted results are shown in Table 2. It is clear that the present technique produces satisfactory results even though the anomaly was contaminated with up to 10% of WGN. Table 1. Synthetic example in arbitrary units Theoretical Model K h  Assumed values 1 10 60 o Interpreted values 1.071 9.94 57.16 o Percentage of error 7.1 6 4.73 MANSOUR A. AL-GARNI 75 Table 2. Synthetic example in arbitrary units, contaminated with 5% and 10% of WGN. Theoretical Model K h  Assumed values 1 10 60 o Interpreted values with 5% of WGN 0.957 9.686 56.433 o Interpreted values with 10% of WGN 0.912 9.082 57.028 o Figure 3. The even component (a), the odd component (b), the Hartley transform (c) and the amplitude spectrum (d) of the horizontal circular cylinder. SPECTRAL ANALYSIS OF MAGNETIC ANOMALIES 75 Figure 4. Response of the vertical magnetic effect of a horizontal circular cylinder with 5% of WGN. Figure 5. Response of the vertical magnetic effect of a horizontal circular cylinder with 10% of WGN. 7. Field example The proposed technique is tested with an example of vertical magnetic anomaly over a narrow band of quartz-magnetite in Mangampalli near Karimnagar town, India (Murthy et al. 1980), as shown in Figure 8. The anomaly is digitized at 15 ft intervals over 540 ft. The even component, odd component, Hartley transform, and the amplitude spectrum are computed and shown in Figure 9. The interpretation parameters, using the procedures mentioned in the text, are tabulated and shown in Table 3. It shows that the results of the proposed technique inversion are in good agreement with the other published ones. MANSOUR A. AL-GARNI 75 Figure 6. The even component (a), the odd component (b), the Hartley transform (c) and the amplitude spectrum (d) of the horizontal circular cylinder anomaly, contaminated with 5% of WGN. SPECTRAL ANALYSIS OF MAGNETIC ANOMALIES 75 Figure 7. The even component (a), the odd component (b), the Hartley transform (c) and the amplitude spectrum (d) of the horizontal circular cylinder anomaly, contaminated with 10% of WGN. MANSOUR A. AL-GARNI 76 Table 3. Field example over a narrow band of quartz-magnetite in Mangampalli near Karimnagar town, India Method Depth in feet Murthy et al. (1980) 78 Sudhakar et al. (2004) 87 Al-Garni (2009) 82.93 Present technique 83.26 Figure 8. Vertical magnetic anomaly over a narrow band of quartz-magnetite in Mangampalli near Karimnagar town, India (Murthy et al. 1980). 8. Conclusion Spectral analysis, using the Hartley transform, of the magnetic anomalies due to a horizontal circular cylinder has been carried out. This approach was applied first to a synthetic data and then to a real data of the vertical magnetic anomaly over a narrow band of quartz-magnetite in Mangampalli near Karimnagar town. The noises effect on the present technique is tested. This technique shows a level of confidence in the quantitative interpretation of the parameters of the vertical magnetic effect of horizontal cylinder anomalies. Due to the fact that the Hartley transform is purely real, it in general has advantages over the conventional spectral analysis (Fourier transform) in terms of its efficient and economical calculations particularly for more sophisticated problems. It is very interesting to notice that the interpreted results of the real data agree well with those obtained by other techniques, published in the literature. 9. Acknowledgment The author thanks Prof. N. Sundararajan, Department of Earth Sciences, Sultan Qaboos University, Sultanate of Oman, for his suggestions to improve the manuscript. SPECTRAL ANALYSIS OF MAGNETIC ANOMALIES 77 Figure 9. The even component (a), the odd component (b), the Hartley transform (c) and the amplitude spectrum (d) of the vertical magnetic anomaly over a narrow band of quartz-magnetite in Mangampalli near Karimnagar town. MANSOUR A. AL-GARNI 78 6. References AL-GARNI, M.A., 2009. Interpretation of Some Magnetic Bodies using Neural Networks Inversion. Arabian Journal of Geosciences 2: 175-184. BRACEWELL, R.N. 1983. The discrete Hartley Transform. Journal of the Optical Society of America 73: 1832- 1835. MOHAN, N.L., SESHAGIRI RAO, S.V., 1982. Spectral Interpretation of Gravity Anomalies Due to Horizontal Slab like Bodies with Lateral Variation of Density. Proc. Ind. Acad. Sci. Earth Planet Sci. 91: 43-54. MOHAN, N.L., SUNDARARAJAN, N. and SESHAGIRI RAO, S.V., 1982. Interpretation of Some Two- Dimensional Magnetic Bodies using Hilbert Transforms. Geophysics 47: 376-387. MURTHY, R.I.V., RAO, V.C., KRISHNA, G.G., 1980. A Gradient Method for Interpreting Magnetic Anomalies Due to Horizontal Circular Cylinders, Infinite Dykes and Vertical Steps: Journal of Earth System Science 89: 31-42. RAJAN, N.S. 1993. Discussion on the use of Hartley transform in geophysical applications. Geophysics 56: 1058-1059. HARTLEY, R.V.L., 1942. A More Symmetrical Fourier Analysis Applied to Transmission Problems. Proc. IRE 30 2: 144-150. SAATCILAR, R., ERGINTAV, S., CANITEZ, N., 1990. The use of the Hartley Transform in Geophysical Applications. Geophysics 55: 1488-1495. SAATCILAR, S., ERGINTAV, S., 1991. Solving Elastic Wave Equations with the Hartley Method. Geophysics 56: 274-278 SUDHAKAR, K.S., RAO, R.P.R., MURTHY, R.I.V., 2004. Modified Werner Deconvolution for Inversion of Magnetic Anomalies of Horizontal Circular Cylinders. J. Ind. Geophys. Union 8: 179-183. SUNDARARAJAN, N., 1995. 2-D Hartley Transforms. Geophysics 60: 262-267. SUNDARARAJAN, N., 1997. Fourier and Hartley Transforms-a Mathematical Twin. Indian Journal of Pure Applied Mathematics 28: 1361-1365. SUNDARARAJAN, N., AL-GARNI, M.A., RAMABRAHMAM, G., SRINIVAS, Y., 2007. A Real Spectral Analysis of the Deformation of a Homogeneous Electric Field Over a Thin Bed - A Hartley Transform Approach. Geophysical Prospecting 55: 901-910. SUNDARARAJAN, N., BRAHMAM, G., 1998. Spectral Analysis of Gravity Anomalies Caused by Slab-like Structures: A Hartley Transform Technique. Journal of Applied Geophysics 39: 53-61. VILLASENOR, J.D., BRACEWELL, R.N., 1987. Optical Phase Obtained by Analog Hartley Transformation. Nature 330: 735-737. Received: 21 September 2010 Accepted: 12 February 2011