磁张量梯度数据方向解析信号比值的边界识别和空间位置反演方法

磁张量梯度测量具有高分辨率、多参量的优点,能更准确地描述磁源体的分布特征,在矿产资源勘探中具有广阔的用途.磁异常解析信号具有受倾斜磁化干扰小的特点,且为了增强深部地质体的分辨能力,本文提出磁张量梯度数据的解析信号比值的均衡边界识别及空间位置反演技术.磁张量梯度数据的均衡边界识别方法为不同方向解析信号比值的反正切函数,在降低倾斜磁化干扰的同时能有效地均衡不同深度地质体的响应,提高了对较深地质体的分辨率;空间位置反演技术是建立解析信号比值与地质体位置参数的对应方程,利用解析信号比值与地质体的对应关系作为约束条件来反演获得地质体的水平位置和深度信息,具有无需已知任何先验信息的优势.通过磁性体张量异常试验表明解析信号比值的边界识别方法能清晰和准确地获得不同深度地质体的边界,所建立的反演方程能准确地计算出地质体的范围和深度,具有较高的水平分辨率和精度.将本文方法应用于实测磁张量梯度数据的解释,获得了地下铁矿的分布特征,为区域矿产资源潜力评价提供了翔实的基础资料.     

The effects of anisotropy of marine magnetic anomalies on the Curie point depth estimates from spectral analysis

Azimuthally averaged power spectra are widely used in the Curie point depth (CPD) estimation with the implicit assumption that the magnetization distribution is random and uncorrelated. However, the marine magnetic anomalies are caused by bands of normal and reverse magnetization and show obvious trends. To investigate the effects of the anisotropy of marine magnetic anomalies on the CPD estimates, we develop 3D fractal striped magnetization models to produce lineated marine magnetic anomalies for the first time. We analyze the spectra anisotropy of the lineated magnetic anomalies of the synthetic fractal striped magnetization models and investigate its effects on the CPD estimates. The synthetic models and actual data show that the spectra of the lineated marine magnetic anomalies are directionally anisotropic. The amplitude response is strong and the slope of the logarithmic spectrum is large in a direction perpendicular to the stripes of magnetic anomalies, whereas the amplitude response is weak and the slope of the logarithmic spectrum is small in a direction parallel to the stripes of magnetic anomalies. The depth estimates in the perpendicular direction are close to the actual values, whereas the depths estimates in the parallel direction are significantly lower than the actual values. The actual marine magnetic anomalies of the South China Sea exhibit an anisotropic power spectrum that is consistent with the spectral anisotropy of magnetic anomalies of the synthetic fractal striped magnetization models.     

Interpretation of magnetic and gravity gradient tensor data using normalized source strength – A case study from McFaulds Lake,Northern Ontario,Canada

In this paper, we present a case study on the use of the normalized source strength (NSS) for interpretation of magnetic and gravity gradient tensors data. This application arises in exploration of nickel, copper and platinum group element (Ni‐Cu‐PGE) deposits in the McFaulds Lake area, Northern Ontario, Canada. In this study, we have used the normalized source strength function derived from recent high resolution aeromagnetic and gravity gradiometry data for locating geological bodies. In our algorithm, we use maxima of the normalized source strength for estimating the horizontal location of the causative body. Then we estimate depth to the source and structural index at that point using the ratio between the normalized source strength and its vertical derivative calculated at two levels; the measurement level and a height h above the measurement level. To discriminate more reliable solutions from spurious ones, we reject solutions with unreasonable estimated structural indices. This method uses an upward continuation filter which reduces the effect of high frequency noise. In the magnetic case, the advantage is that, in general, the normalized magnetic source strength is relatively insensitive to magnetization direction, thus it provides more reliable information than standard techniques when geologic bodies carry remanent magnetization. For dipping gravity sources, the calculated normalized source strength yields a reliable estimate of the source location by peaking right above the top surface. Application of the method on aeromagnetic and gravity gradient tensor data sets from McFaulds Lake area indicates that most of the gravity and magnetic sources are located just beneath a 20 m thick (on average) overburden and delineated magnetic and gravity sources which can be probably approximated by geological contacts and thin dikes, come up to the overburden.     

DEPTH ESTIMATION OF MAGNETIC SOURCES BY MEANS OF FOURIER AMPLITUDE SPECTRA*

The Fourier transform of a square-shaped section of a magnetic survey, digitized in a square grid, forms a rectangular matrix of coefficients which can be condensed to a series of average amplitudes dependent only on their frequency and no longer on the direction of the respective partial waves. These average amplitudes together represent a spectrum which–when plotted in a semilogarithmic coordinate system (log amplitude versus frequency)–often shows straight segments which decrease with increasing frequency. By continuing the given field downwards these straight segments become horizontal at a certain depth, the so-called “white depth”. This white depth may be used as a first estimate for the depth of magnetic sources producing the respective part of the field. It is shown that the sources which correspond to such use of the white depth can be expected to be “randomly distributed with some positive autocorrelation”. As an example for such a depth estimation the interpretation of the aeromagnetic survey of NW-Germany by a relief in 8–16 km depth is given. The relief divides the subsurface in an upper nonmagnetic layer and a lower layer with magnetization M= 2 Am^?1.     

Continuous wavelet transform,theoretical aspects and application to aeromagnetic data at the Huanghua Depression,Dagang Oilfield,China

We use the continuous wavelet transform based on complex Morlet wavelets, which has been developed to estimate the source distribution of potential fields. For magnetic anomalies of adjacent sources, they always superimpose upon each other in space and wavenumber, making the identification of magnetic sources problematic. Therefore, a scale normalization factor, a^?n, is introduced on the wavelet coefficients to improve resolution in the scalogram. By theoretical modelling, we set up an approximate linear relationship between the pseudo‐wavenumber and source depth. The influences of background field, random noise and magnetization inclination on the continuous wavelet transform of magnetic anomalies are also discussed and compared with the short‐time Fourier transform results. Synthetic examples indicate that the regional trend has little effect on our method, while the influence of random noise is mainly imposed on shallower sources with higher wavenumbers. The source horizontal position will be affected by the change of magnetization direction, whereas the source depth remains unchanged. After discussing the performance of our method by showing the results of various synthetic tests, we use this method on the aeromagnetic data of the Huanghua depression in central China to define the distribution of volcanic rocks. The spectrum slices in different scales are used to determine horizontal positions of volcanic rocks and their source depths are estimated from the modulus maxima of complex coefficients, which is in good accordance with drilling results.     

全张量磁梯度数据解释的均衡边界识别及深度成像技术

全张量磁梯度数据具有高精度、高分辨率、多参量的优点,能更加清晰地刻画地质体的分布特征,综合利用磁张量梯度数据准确地获得地质体水平位置和深度信息是解释的主要目的.磁张量数据的方向解析信号具有减小倾斜磁化干扰的优点,常被用来圈定磁源体的水平位置,但解析信号强度随着地质体埋深的增加急剧衰减,难以有效识别较深的地质体.张量数据均衡边界识别技术,利用不同方向解析信号的比值函数,能有效地均衡不同深度地质体的响应,同时显示不同深度地质体的边界,提高了对较深地质体的分辨率.磁张量数据深度成像技术根据实测张量数据与假定模型张量数据的相关系数来给定地质体的深度,综合利用多参量数据联合反演提高了反演结果的准确性,且无需进行复杂的反演运算,是大数据量张量数据解释的有效方法.理论模型试验证明：磁张量数据均衡边界识别技术可清晰和准确地识别地质体的水平范围,受倾斜磁化干扰小;磁张量数据深度成像技术可准确地获得地质体的深度信息,具有较强的抗噪性.将上述方法应用于铁矿区实测航磁张量梯度数据解释,获得了铁矿体水平分布与埋深,深度结果与张量欧拉反褶积法计算结果一致.     

On the use of the autocorrelation method for estimating the depth to the magnetic basement

The depth to the top of magnetic dykes can be estimated from total field aeromagnetic data using the relation between the depth to magnetic sources and the autocorrelation function of magnetic data. By using synthetic anomalies we show that in the ideal case, depth can be determined to an accuracy of 10% or better, when the anomaly sources are two-dimensional dykes. However, the estimated depths depend on the width of the dykes. The estimated depth is about 0.6 times the actual depth to the top of thin dykes, and around the true depth for thick dykes having width-to-depth ratio around 3. The depth is considerably overestimated for very thick dykes (e.g., contacts, which is a special case of the thick dyke). Thus, the autocorrelation method requires that the width-to-depth ratio of the dyke is estimated independently to correctly estimate the depths. Alternatively, it must be assumed that the width-to-depth ratio for the two-dimensional source body is between 1.5 and 4.     

Satellite and surface geophysical expression of anomalous crustal structure in Kentucky and Tennessee

An equivalent layer magnetization model obtained from inversion of long-wavelength satellite magnetic anomaly data indicates a very magnetic source region centered in south central Kentucky. The magnetization maximum nearly coincides with a gravity high elongated north-south and extending into Tennessee. Previous refraction profiles suggest that the source of the gravity anomaly is a large mass of rock occupying much of the crustal thickness. The outline of the source delineated by gravity contours is also discernible in aeromagnetic anomaly patterns. Taken together, the geophysical data suggest a large, localized mass of intracrustal rock which is both dense and very magnetic. A simple magnetization/density model is given which accounts for the gravity and long-wavelength aeromagnetic anomalies due to the body. We interpret it as a mafic plutonic complex, and several lines of evidence are consistent with a rift association. The body is, however, clearly related to the inferred position of the Grenville Front. It is bounded on the north by the fault zones of the 38th Parallel Lineament. The inferred mean magnetization (4 A/m) of the body is large, but not inconsistent with values reported by others for deep crustal bodies associated with long-wavelength magnetic anomalies. Such magnetization levels can be achieved with magnetic mineralogies produced by normal oxidation and metamorphic processes and enhanced by viscous build-up, especially in mafic rocks of alkaline character.     

Modelling of magnetic data for scaling geology

Interpretation of magnetic data can be carried out either in the space or frequency domain. The interpretation in the frequency domain is computationally convenient because convolution becomes multiplication. The frequency domain approach assumes that the magnetic sources distribution has a random and uncorrelated distribution. This approach is modified to include random and fractal distribution of sources on the basis of borehole data. The physical properties of the rocks exhibit scaling behaviour which can be defined as P(k) = Ak^?β, where P(k) is the power spectrum as a function of wave number (k), and A and β are the constant and scaling exponent, respectively. A white noise distribution corresponds to β = 0. The high resolution methods of power spectral estimation e.g. maximum entropy method and multi‐taper method produce smooth spectra. Therefore, estimation of scaling exponents is more reliable. The values of β are found to be related to the lithology and heterogeneities in the crust. The modelling of magnetic data for scaling distribution of sources leads to an improved method of interpreting the magnetic data known as the scaling spectral method. The method has found applicability in estimating the basement depth, Curie depth and filtering of magnetic data.     

基于最小反演拟合差的重磁场源深度计算方法

以等效源及位场物性反演为基础,本文提出一种新的求取重磁场源深度的方法.该方法将一层等效源以一定的间隔从浅部向深部移动,并将等效源作为初始模型进行反演,当反演拟合差最小时,停止反演,此时的等效源底深即为所求场源的中心深度.由于仅需要反演一层等效源,比传统的物性反演计算时间大大减少,并且不需要进行深度加权约束.理论模型数据处理结果表明该方法能够获得较准确的场源深度:以长宽比为7.5的薄板模型为例,深度计算误差约为1个点距(25 m);以长宽比为0.5~1.5的厚板模型为例,深度计算误差小于1个点距(25m).将该方法应用于实测航磁梯度数据,计算的磁源中心深度在200~250m之间,钻井资料显示该异常由埋藏深度在200~300m的闪长岩引起,计算结果与钻井资料较吻合.     

JOINT INTERPRETATION OF GRAVITY AND MAGNETIC DATA OVER AXIAL SYMMETRIC BODIES WITH APPLICATION TO THE DARNLEY BAY ANOMALY,NWT CANADA*

A method is presented for determining bounds of the properties of axial symmetric bodies from a finite number of gravity and magnetic observations based on Parker's theory of ideal bodies. Bounds on the density contrast and the intensity of magnetization are calculated as a function of depth to the top of the anomalous source, restricting the range of smallest possible solutions to fit the data. The model studied is approximated by an array of vertical annuli cylinders, each of uniform density and magnetization. Linear programming algorithms based on the ideal body theory were used to calculate the distribution of these parameters within the body. Simultaneous inversion of gravity and magnetic data is performed assuming a constant ratio between the density contrast and the intensity of magnetization and that a common body is responsible for both observed fields. The parameter k(|J|/δp) provides information about the rock type of the structure. Interpretation of gravity and aeromagnetic data from Darnley Bay, NWT, Canada, indicated the presence of a shallow ultrabasic intrusion.     

Interpretation of the Gravity and Magnetic Anomalies of the Cappadocia Region,Central Turkey

The Cappadocia region, located in Central Turkey, is characterized by widespread lava flows and volcanoclastic deposits dating from Miocene to Quaternary. Gravity and aeromagnetic anomalies of the region appear to present similar high and low amplitude regions, although the aeromagnetic anomalies exhibit a rather complex pattern which is thought to be caused by remanent magnetization. The low-pass filtered aeromagnetic map shows a deep-seated magnetic anomaly which may be linked to the widespread volcanic activity at the surface. The pseudogravity transformation of the upward continued anomaly has been constructed. The pseudogravity anomaly demonstrates some form of clockwise rotation. This anomaly was modelled by means of a three-dimensional method. The top and bottom of the body are at 6.3km and 11km (including the flight height) from the ground surface, respectively. This deep body is ellipsoidal and extends along an E-W direction, which is in line with the regional stress direction deduced from GPS measurements. A new mobilistic dynamo-tectonic system appears to explain the body’s E-W elongation. The modelled body may be the source for the inferred geothermal energy of the region. Magnetic measurements were carried out on oriented rock samples collected from outcrops of ignimbrites and basalts, providing directions and intensities of remanent magnetization, susceptibilities and Koeningsberger (Q) ratios. Standard deviations of remanent directions of the Natural Remanent Magnetization (NRM) display a wide scatter implying unreliability of the surface data. Reduction to pole (RTP) transformation of magnetic anomalies was successful with the induced magnetization angle despite the complex pattern of magnetic anomalies.     

全张量磁梯度数据的斜导数特征值边界识别方法研究

本文提出了一种全张量磁梯度数据的斜导数边界识别方法,该方法利用全张量磁梯度数据定义了北向斜导数、东向斜导数和垂向斜导数.模型试验表明,北向和东向斜导数无法有效识别磁性异常体边界,异常形态复杂虚假异常较多,而垂向斜导数在高纬度斜磁化或者垂直磁化条件下可以清晰而准确得识别地质体的边界,具有一定的深度识别能力.与现有磁张量梯度的边界识别方法相比,识别效果好,可以有效均衡不同深度地质体的响应.将上述方法应用于大兴安岭地区实测航磁三分量数据转换得到的全张量数据,获得了研究区浅部磁性体的分布,并与三分量数据进行联合解释,取得良好的地质效果.     

Improved tilt-depth method for fast estimation of top and bottom depths of magnetic bodies

The tilt-depth method can be used to make fast estimation of the top depth of magnetic bodies. However, it is unable to estimate bottom depths and its every inversion point only has a single solution. In order to resolve such weaknesses, this paper presents an improved tilt-depth method based on the magnetic anomaly expression of vertical contact with a finite depth extent, which can simultaneously estimate top and bottom depths of magnetic bodies. In addition, multiple characteristic points are selected on the tilt angle map for joint computation to improve reliability of inversion solutions. Two- and threedimensional model tests show that this improved tilt-depth method is effective in inverting buried depths of top and bottom bodies, and has a higher inversion precision for top depths than the conventional method. The improved method is then used to process aeromagnetic data over the Changling Fault Depression in the Songliao Basin, and inversion results of top depths are found to be more accurate for actual top depths of volcanic rocks in two nearby drilled wells than those using the conventional tilt-depth method.     

龙门山和相邻地域航磁场特征与汶川大地震

本文应用化极、水平及垂向导数、向上延拓、视磁化强度填图及磁性界面反演等方法处理了龙门山及相邻地域最新的航空磁测数据,分析了龙门山及相邻地域的航磁异常展布特征.研究结果表明:1)龙门山造山带与其东、西两侧可划分为三个磁异常区:松潘—甘孜磁异常区、龙门山负磁异常带、四川盆地磁异常区;三个区、带的地壳介质磁性结构存在明显差异.2)根据该区航空磁异常场的分布特征分别研究了,松潘—甘孜地块、龙门山造山带和四川地块的磁场特征.3)除识别前人识别的断层外,还推断鲜水河ES延伸甘洛—雷波北断裂作为四川盆地与滇西的界带.4)航空磁异常,磁性体上、下界面及磁源体深度的空间分布特征与汶川MS8.0大地震及芦山地震发生相关.     

HARMONIC INVERSION OF AN AEROMAGNETIC ANOMALY FROM A REMANENTLY MAGNETIZED PRECAMBRIAN SCHIST BELT (KARNATAKA,INDIA) *

The conventional methods of magnetic interpretation consider in general the induced magnetization only. However, in a computer oriented inverse method, the direction of magnetization can be taken as arbitrary. Following this procedure, this study aims to present a realistic interpretation of an aeromagnetic anomaly from Chitaldrug Schist belt (Mysore, India) assuming a three-dimensional configuration of the causative source possessing remanent magnetization. It further provides the direction of remanent magnetization of the source rock as obtained from the analysis of the magnetic anomaly. The total intensity map obtained delineates the surface and sub-surface structures in the region and provides informations regarding the characteristics of primordial crust.     

Estimation of depth and magnetization of magnetic sources from magnetic data: The linearized continuous inverse problem for 21/2D structures

The estimation of the depth to the top and bottom of a magnetic source from magnetic data defines a nonlinear inverse problem, while the evaluation of the distribution of magnetization determines a linear inverse problem. In this paper, these interpretation problems are resolved in the continuous case of 2_1/2D magnetized bodies with lateral magnetization variations. A formulation of the magnetic problem accounting for different directions of remanent and total magnetization vectors and including a more general definition of apparent susceptibility is presented. Differences between 2D and 2_1/2D formulations are stressed, as regards the anomaly amplitude, shape and zero-level.In order to utilize well-known continuous linear inverse methods, Fréchet derivatives of the data functionals with respect to the depth of the source top and bottom, are analytically described. Thus, using the spectral expansion inverse method (Parker, 1977) and linearizing the problem at several steps of an iterative process, the source depth is obtained within a few iterations, although the starting model is distant from the final solution. The interpretation of an anomaly in the Italian region shows the usefulness of the method.     

Line integrals of potential field data1

Information on the mass and the spatial location of an arbitrary source body can be obtained by performing suitable integrations of 3D gravity and magnetic data along an infinite straight line. No assumptions on the density/magnetization distribution or the shape and location of the source are required. For an oblique borehole, a relationship between the lower limit of the source mass and the distance to the body is obtained. The mass contrast and the magnetic moment of the source can also be estimated. For a vertical borehole, both gravity and vertical magnetic component anomalies have equal areas to the left and right of the depth axis. The particular case of a horizontal gallery not intersecting the body is also studied. If the source is intersected, a lower limit is estimated for the maximum thickness of the body along the gallery. Information on the vertical coordinate of the centre of mass of the source can also be obtained. Numerical tests with synthetic gravity data support the theoretical results.     

基于改进Parker-Oldenburg界面反演方法计算青海省居里面深度

居里面深度的起伏可以大致反映地壳深部温度场的分布特征,从而间接指示干热岩的赋存情况,因此研究居里面的起伏特征十分重要.本文针对经典的Parker-Oldenburg界面反演方法存在的计算不稳定、结果精度低的问题,通过公式推导,将界面反演公式与向下延拓公式进行类比,利用稳定的波数域正则-积分迭代下延法对经典的Parker...     

3D Inversion of magnetic data of grouped anomalies — Study applied to São José intrusions in Mato Grosso,Brazil

The southwest border of the Parecis Basin (central Brazil) presents several occurrences of gold, copper and zinc. Parallel to this border, there is an alignment of magnetic anomalies with varied size and polarities. In particular, five magnetic anomalies are referred to, in this study, as SJ1 to SJ5. The proximity of these anomalies to each other makes it hard to isolate the magnetic component associated with each source. Furthermore, these anomalies have different magnetization directions, which require the use of a technique which is slightly or not affected by the presence of a remanent magnetization, as the amplitude of the anomalous magnetic field. Considering that, in these intrusions no outcrops are observed, the enhanced horizontal derivative technique was used to estimate the edges location and the depth of these sources. The geological context, allied to the results from the magnetic techniques, allowed to establish binds to restrain the interpretation of the results of the 3D inversion. This procedure permitted to compose three hypothesis to explain the magnetic behavior of the region, from which can be conclusively determined with a borehole analysis.