برآورد مرز بي‌هنجاري‌هاي ميدان پتانسيل با مشتق قائم سيگنال تحليلي

Edge detection of potential field anomalies using vertical derivative of analytic signal

The analytic signal for magnetic anomalies was initially defined as a "complex field deriving from a complex potential" (Nabighian, 1972). This function can be computed easily in the frequency domain, its real part is the horizontal derivative of the field and its imaginary part is the vertical derivative. Analytic signal processing and interpretation requires few initial assumptions regarding the source body geometry and magnetization and is particularly efficient at an early stage of the interpretation even if constraints are not available. For 2-D structures (Nabighian, 1972), the method assumes that the causative bodies have a polygonal cross-section with uniform magnetization. Such structures can also be considered as the superimposition of a finite number of magnetic steps. Narrow dikes and thin sheets can also be taken into account using a lower order of derivation; for example, the field itself instead of the horizontal derivatives. Nabighian (1972) demonstrates that the analytic signal has simple poles at each corner of the structures. The amplitude of the analytic signal is a bell-shaped symmetric function maximizing exactly over the top of each contact, with the width of the amplitude curve being related directly to the depth of the contact. This is also true for any of the derivatives of the signal (Nabighian, 1974); these properties can be used to locate the magnetic contacts and to estimate their depths. E-mail: kamal.alamdar@gmail.com .via •r^)^-Ar^.\\t> cr ojU-i CTr iLai j j iJL>«j» Extension of the 2-D analytic signal to three dimensions will allow more general interpretation procedures to be developed, the two-dimensionality assumption being no longer required. The relationship between the horizontal and vertical derivatives for the 3- D case was first derived by Nabighian (1984). As with the 2-D case, the vertical and horizontal derivatives are Hilbert transforms of each other and so the concept of the analytic signal in 2-D can be extended to three dimensions. The analytic signal amplitude can then be defined as "the square root of the squared sum of the vertical and the two horizontal first derivatives of the magnetic field" (Roest et al., 1992). (1) Because of interference effects, the use of the simple analytic signal in the 3-D case seems insufficient to detect geologic boundaries. Since the existing interference is usually inevitable, improving resolution becomes a requirement. In the 2-D case, Nabighian (1974) suggested using the following bell-shaped function to enhance the analytic signal from shallow sources: S"Gh) ( 8"G. "I — + — J [dh» J _(l2 x22 x3 2 x...xn2)a2 (d2+h2fX where Gh and Gz are the horizontal and vertical gradients of the potential-field anomaly, respectively; h is the distance along the horizontal axis which is perpendicular to the strike of the 2-D structure; n is a positive integer; d is the depth to the top surface of the source, while the lower surface is at infinity; a is the ambient parameter and is equal to a = 2KFc sin 0 when the step model of magnetic anomaly is applied; k is the susceptibility contrast of the step model; F is the earthʹs magnetic field magnitude; 0 is the dip angle of the step model; c = 1 - cos2 i sin2 \|/ for total magnetic field anomalies; i is the inclination of the earthʹs magnetic field; and V|/ is the angle between magnetic north and positive h - axis. All the above parameters are presented in Fig. 1. In the 3-D case, the simple analytic signal is defined in Nabighian (1984) as: