Separation of magnetic anomalies into induced and remanent magnetization contributions

Inversion of magnetic data is complicated by the presence of remanent magnetization, and it provides limited information about the magnetic source because of the insufficiency of data and constraint information. We propose a Fourier domain transformation allowing the separation of magnetic anomalies into the components caused by induced and remanent magnetizations. The approach is based on the hypothesis that each isolated source is homogeneous with a uniform and specific Koenigsberger ratio. The distributions of susceptibility and remanent magnetization are subsequently recovered from the separated anomalies. Anomaly components, susceptibility distribution and distribution of the remanent and total magnetization vectors (direction and intensity) can be achieved through the processing of the anomaly components. The proposed method therefore provides a procedure to test the hypotheses about target source and magnetic field, by verifying these models based on available information or a priori information from geology. We test our methods using synthetic and real data acquired over the Zhangfushan iron-ore deposit and the Yeshan polymetallic deposit in eastern China. All the tests yield favourable results and the obtained models are helpful for the geological interpretation.     

New digital linear filters for Hankel J0 and J1 transforms

The numerical evaluation of certain integral transforms is required for the interpretation of some geophysical exploration data. Digital linear filter operators are widely used for carrying out such numerical integration. It is known that the method of Wiener–Hopf minimization of the error can be used to design very efficient, short digital linear filter operators for this purpose. We have found that, with appropriate modifications, this method can also be used to design longer filters. Two filters for the Hankel J₀ transform (61-point and 120-point operators), and two for the Hankel J₁ transform (47-point and 140-point operators) have been designed. For these transforms, the new filters give much lower errors compared to all other known filters of comparable, or somewhat longer, size. The new filter operators and some results of comparative performance tests with known integral transforms are presented. These filters would find widespread application in many numerical evaluation problems in geophysics.     

Detection and identification of north–south trending magnetic structures near the magnetic equator

Long, structurally undeformed north–south trending structures show no magnetic anomaly at the magnetic equator, except at the north and south truncations of the structure. However, folding, faulting, differential erosion or other structural deformation can produce detectable magnetic anomalies in a generally north–south trending equatorial structure. Spatial variation in magnetic susceptibility or remanent magnetization can also produce anomalies in equatorial north–south structures. These anomaly patterns are often more complicated than patterns produced by similar structures at high latitudes, but interpretational insight can be gained through numerical modelling of common structures. Reduction-to-pole and analytic signal filters can aid in interpretation of equatorial anomalies, but these must be applied carefully because of instabilities deriving from filter design and noise amplification.     

A CONTROL OF TWO-DIMENSIONAL MAGNETIC INTERPRETATION BY THREE-DIMENSIONAL MODEL BODY ANOMALIES*

In magnetic routine interpretation the comparison of two-dimensional model curves with measured magnetic anomalies is widely used for an approximate evaluation of the position and depth of magnetic models. Before starting an interpretation of a survey by means of two-dimensional models, it is very useful to have an idea of the shape of anomalies caused by extended but finite bodies, taking into account various strike directions: Three sets of anomalies of thin plates (horizontal length 19, downward length 9, width 1) dipping 30°, 60°, and 90° resp. for various strike directions and an inclination of 20° were computed. Some of these anomalies, e.g. those with nearly N-S strike direction look rather complicated, and at the first glance one would not expect that they are caused by such simple bodies. Several profiles crossing the computed anomalies perpendicularly were interpreted two-dimensionally. For less extended anomalies the depths determined for the top of the plates are 10-20% too small, the magnetization amounts to 50–75 % of the value of the finite bodies. The interpretation of the profiles covering more extended anomalies gave very accurately the same values for the position, depth and magnetization for the two-dimensional body as for the original three-dimensional model. Anomalies of vertical prisms with varying extensions in the y-direction were computed. Their differences in amplitude and in the distance maximum-minimum show that interpretation of short anomalies by two-dimensional methods yields depth errors of up to 20 percent. To see the possibilities of the separation of superimposed anomalies dike anomalies were added to the anomaly of a broad body in great depth and several attempts were made to interpret parts of the composite anomalies. The interpreted bodies lie too deep. In complicated cases the depth values can have large errors, but experienced interpreters should be able to keep the errors in the range of one third of the depth values.     

南海北部盆地基底岩性地震-重磁响应特征与识别

针对性选取东南沿海露头剖面18条,采集245件南海盆地基底可能出现的岩性样品,测定其密度和磁化率,建立各种岩性的密度-磁化率交会图版,以此约束过井地震剖面和重磁异常的地质解释,总结出南海北部盆地基底火山岩、侵入岩、变质岩和沉积岩4大类11亚类岩性的地震-重磁响应特征.应用重磁震-岩性解释模型逐一对南海盆地北部主干剖面进行地质-地球物理综合解释,从而实现了盆地基底岩性的平面填图.这种从盆缘剖面到盆地内部、从岩石物性测量到地质-地球物理综合解释的方法,在资料获取难度大、地质条件复杂的南海盆地基底地质研究中,业已证明是行之有效的,相信在其他盆地研究中也会有借鉴意义.     

Depth determination from a non-stationary magnetic profile for scaling geology

The conventional spectral analysis method for interpretation of magnetic data assumes stationary spatial series and a white‐noise source distribution. However, long magnetic profiles may not be stationary in nature and source distributions are not white. Long non‐stationary magnetic profiles can be divided into stationary subprofiles following Wiener filter theory. A least‐squares inverse method is used to calculate the scaling exponents and depth values of magnetic interfaces from the power spectrum. The applicability of this approach is demonstrated on non‐stationary synthetic and field magnetic data collected along the Nagaur–Jhalawar transect, western India. The stationarity of the whole profile and the subprofiles of the synthetic and field data is tested. The variation of the mean and standard deviations of the subprofiles is significantly reduced compared with the whole profile. The depth values found from the synthetic model are in close agreement with the assumed depth values, whereas for the field data these are in close agreement with estimates from seismic, magnetotelluric and gravity data.     

同时考虑退磁和剩磁的有限体积法正演模拟

为分析同时考虑退磁和剩磁对磁测数据解释的影响,探讨了利用有限体积法求静磁场数值解的方法.从静磁场中的麦克斯韦方程出发,导出了有限体积法控制方程的离散表达式,对边界条件近似处理后求解方程组得到磁异常.通过与退磁改正计算结果对比,验证了方法的正确性,并分析得到忽略剩磁的相对误差与科尼斯布格比（Q）相关;利用有限体积法计算长方体模型在无地磁场情况下的磁异常和内部磁化强度,从数值模拟上说明剩磁也需要进行退磁改正,并表明退磁作用对剩磁的影响不仅与磁化率相关,而且与剩磁的方向和磁性体的形态相关;组合模型的计算结果对比表明,退磁作用对剩磁的影响还会因为临近强磁性体的作用发生改变.在青海灶火河西工区的应用说明,开展同时考虑退磁和剩磁的解释方法对准确识别强磁性岩体具有实用价值.     

A magnetic model for deep plutonic bodies beneath the central Taupo Volcanic Zone, North Island, New Zealand

The central Taupo Volcanic Zone (TVZ) is a region of intense Quaternary rhyolitic volcanism and geothermal activity in the North Island of New Zealand from which about 14,000 km³ of pyroclastics and lavas have been erupted during the last 1.6 Ma. Analysis of aeromagnetic surveys over the TVZ showed the presence of long-wavelength (10 to 25 km) magnetic anomalies which roughly follow the trend of the currently active eastern TVZ, from the north of Lake Taupo to the east of Lake Rotorua. An interpretation of the long-wavelength magnetic anomalies using 3-D magnetic modelling suggests that these anomalies are caused by the magnetic effects of < 3 km thick sequence of volcanic rocks and deeper magnetised bodies within the non-magnetic upper crust (4–7 km depth) beneath the young (age < 0.7 Ma), currently active eastern TVZ. The deep magnetised bodies are interpreted as solidified rhyolitic sub-volcanic plutons that have cooled down to below their Curie temperature.Although the existence of plutonic bodies beneath the TVZ has been postulated prior to this study, this magnetic interpretation result appears to be the first geophysical model of such bodies.     

Magnetic and gravimetric modeling of the central Adriatic region

The high-amplitude wide magnetic anomaly that covers a large area of the north-eastern – central Adriatic Sea in the Croatian offshore is ∼100 km wide and extends NW–SE for ∼350 km. The anomaly is located between the Dinarides and Apennines chains, in an interesting geodynamic scenario. The presence of intruded gabbroid rocks in the Croatian archipelago also contributes to making this intriguing and still not extensively investigated anomaly potentially significant to the geologic and geophysical context in which it is located. In this work, we model the Bouguer and magnetic anomalies across the Adriatic Sea. The 2D geophysical modeling was produced across four cross sections considering surface heat flow data to calculate the Curie depth. The magnetic susceptibilities and densities used for the synthetic bodies are in agreement with the literature and with those derived by previous models. The results suggest the presence of an uplifted magnetized basement with high magnetic susceptibility (0.075 SI units) to be the main contributor to the observed magnetic anomaly. This magnetic susceptibility is interpreted as representative of a gabbroid-intruded basement. The high-susceptibility basement is in lateral continuity with a relatively low susceptibility basement (0.025–0.038 SI units). The results of the geophysical modeling are compared with a conceptual geological model realized from the integration of surface, well and geophysical data, the latter concerning seismic, tomographic, magnetic and gravimetric anomalies and heat flow data. These data have been merged in an integrated data-base using MOVE software, and the geophysical modeling was performed using GM-SYS. The comparison allowed to confirm the hypothesis that the magnetic anomaly is related to the basement and to its position in the complex geodynamic evolution of the Apennines–Adriatic–Dinarides system.     

Use of non-linear filtering for the regional–residual separation of potential field data

Regional–residual separation is essential in gravity and magnetic data interpretation and a variety of techniques have been proposed. Graphical determination of the regional allows geological information to be taken into account. Upward continuation can be used to obtain the regional field either empirically or using some hypothesis about the geology. In some cases, a matched filter can be designed and used to separate deep and shallow sources. Simple low pass filtering has also been used but without much success. Here we propose to use a non-linear filter approach to remove gravity and magnetic anomalies smaller than a given width. This technique attempts to mimic the graphical separation method. The results from synthetic models are presented as well as the results from a case study in eastern Canada and compared to regional gravity and magnetic anomalies obtained by other techniques. Contrary to the regional fields obtained by upward continuation, non-linear filtering does not have any physical meaning. However, its main advantage is that it gives a regional component of the gravity or magnetic field similar to the one obtained from a graphical separation.     

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

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

Interpretation of total intensity magnetic anomalies when anomalous field is arbitrarily large

Summary Not infrequently, in mining geophysics, the anomalous field of the magnetized body is appreciably large and it varies from the direction of the earth's normal field within the vicinity of the magnetized body. Total magnetic intensity data collected on the ground over shallow magnetized bodies cannot be interpreted quantitatively, since all the available methods of interpretation assume that the resultant field lies in the direction of earth's normal field as the anomalous field is small. For geological bodies that can be magnetically represented by infinite line pole, point pole, infinite line dipole, and point dipole, this paper gives quantitative interpretation procedures for total magnetic intensity anomalies without the above assumption using the characteristic points from the theoretical curves (not reproduced). These characteristic curves can be used for the direct determination of depths and offsets by using information from suitably chosen field profiles.Contribution No. N.G.R.I.-71-233.     

Analysis and interpretation of anomalous conductivity and magnetic permeability effects in time domain electromagnetic data: Part I: Numerical modeling

Time domain electromagnetic (TDEM) response is usually associated with eddy currents in conductive bodies, since this is the dominant effect. However, other effects, such as displacement currents from dielectric processes and magnetic fields associated with rock magnetization, can contribute to TDEM response. In this paper we analyze the effect of magnetization on TDEM data. We use a 3-D code based on finite-difference method, developed by Wang and Hohmann [Geophysics 58 (1993) 797], to study transient electromagnetic field propagation through a medium containing bodies with both anomalous conductivity and anomalous magnetic permeability. The remarkable result is that the combination of anomalous conductivity and permeability within the same body could increase significantly the anomalous TDEM response in comparison with purely conductive or purely magnetic anomalies. This effect has to be taken into account in interpretation of TDEM data over electrical inhomogeneous structures with potentially anomalous magnetic permeability.     

Mapping basement structures in the northwestern offshore of Abu Dhabi from high‐resolution aeromagnetic data

This paper presents a case study of mapping basement structures in the northwestern offshore of Abu Dhabi using high‐resolution aeromagnetic data. Lineament analysis was carried out on the derivatives of the reduced‐to‐the‐pole magnetic data, along with supporting information from published geologic data. The lineament analysis suggests three well‐defined basement trends in the north–south, northeast–southwest, and northwest–southeast directions. The reduced‐to‐the‐pole magnetic data reveal high positive magnetic anomalies hypothesized to be related to intra‐basement bodies in the deep seated Arabian Shield. Depth to basement was estimated using spectral analysis and Source Parameter Imaging techniques. The spectral analysis suggests that the intruded basement blocks are at the same average depth level (around 8.5 km). The estimated Source Parameter Imaging depths from gridded reduced‐to‐the‐pole data are ranged between 4 km and 12 km with a large depth variation within small distances. These estimated depths prevent a reliable interpretation of the nature of the basement relief. However, low‐pass filtering of the horizontal local wavenumber data across two profiles shows that the basement terrain is characterized by a basin‐like structure trending in the northeast–southwest direction with a maximum depth of 10 km. Two‐dimensional forward magnetic modelling across the two profiles suggests that the high positive magnetic anomalies over the basin could be produced by intrusion of mafic igneous rocks with high susceptibility values (0.008 to 0.016 SI.     

Use of the Two Dimensional Gabor Filter to Interpret Magnetic Data Over the Marmara Sea, Turkey

Image processing techniques have been frequently used in processing of potential field data. The most important aim in processing potential field data is to figure out the edges of the geological bodies. The Gabor filter was used to determine the discontinuity boundaries of the source bodies causing the magnetic anomalies in this study. The synthetic total magnetic anomaly of three prism bodies with different depths and orientations was used to test the effect of the Gabor filter on the magnetic data. In addition, the Gabor filter was applied on the reduction to the pole aeromagnetic field data of the Marmara Sea and the results were compared with the fault distribution map of the area prepared by various authors.     

SIMULTANEOUS INTERACTIVE MAGNETIC AND GRAVITY INVERSION*

Gravity and magnetic anomalies may have the same source but it is always difficult to achieve correlated solutions if interpretations are carried out separately. Therefore it is useful to invert both anomalous gravity and magnetic profiles at the same time, so that the solution mav be more constrained. Existing inversion techniques do not lend themselves equally to this kind of extension, since the parameters—such as density and susceptibility contrasts—should not be related. The algorithms more easily adaptable to simultaneous inversion are those which look for the shape and the position of anomalous bodies with constant unknown density or susceptibility contrasts. In this study, we use the generalized inverse method in the 2 1/2 D case. Parameters are the coordinates of the vertices of the polygonal cross-sections of the anomalous bodies and both density and susceptibility contrasts of each body. The two types of profile to be processed must be superposable but may have different lengths, spacings, and coordinate origins. Both synthetic examples and field data from geologically known areas have been processed, and it appears that the simultaneous inversion technique may provide an important additional degree of control in the interactive interpretation process.     

Estimating magnetic dike parameters using a non‐linear constrained inversion technique: an example from the Särna area,west central Sweden

In this paper, we describe a non‐linear constrained inversion technique for 2D interpretation of high resolution magnetic field data along flight lines using a simple dike model. We first estimate the strike direction of a quasi 2D structure based on the eigenvector corresponding to the minimum eigenvalue of the pseudogravity gradient tensor derived from gridded, low‐pass filtered magnetic field anomalies, assuming that the magnetization direction is known. Then the measured magnetic field can be transformed into the strike coordinate system and all magnetic dike parameters – horizontal position, depth to the top, dip angle, width and susceptibility contrast – can be estimated by non‐linear least squares inversion of the high resolution magnetic field data along the flight lines. We use the Levenberg‐Marquardt algorithm together with the trust‐region‐reflective method enabling users to define inequality constraints on model parameters such that the estimated parameters are always in a trust region. Assuming that the maximum of the calculated g_zz (vertical gradient of the pseudogravity field) is approximately located above the causative body, data points enclosed by a window, along the profile, centred at the maximum of g_zz are used in the inversion scheme for estimating the dike parameters. The size of the window is increased until it exceeds a predefined limit. Then the solution corresponding to the minimum data fit error is chosen as the most reliable one. Using synthetic data we study the effect of random noise and interfering sources on the estimated models and we apply our method to a new aeromagnetic data set from the Särna area, west central Sweden including constraints from laboratory measurements on rock samples from the area.     

叠前三维地震近地表吸收补偿技术在大庆油田的应用研究(英文)

与深层致密岩石相比,相对疏松的近地表地层严重吸收了地震波的高频成分,降低了地震数据的分辨率和高频成分的信噪比。本文利用井中雷管激发、地面检波器接收的微测井直达波资料,通过分析近地表地层不同传播距离的地震直达波频谱信息差异,采用维纳滤波方法,进行近地表吸收补偿反滤波器的求取,将不同微测井测量点对应的近地表吸收补偿反滤波器应用于相应的叠前共检波点道集地震数据,完成了叠前地震数据的空变吸收衰减补偿,克服了叠后地震数据无法实现空变补偿的难题。叠前三维地震近地表吸收补偿后的数据,较补偿前地震数据的优势信噪比频带宽度明显拓宽,低频成分基本保持,反射信息量有较大程度增加,而且与合成记录吻合更好,能更好的满足地质解释的需要,提高解释精度。     

三维井地磁测联合约束反演

井中磁测在铁矿等磁性金属矿床勘查中具有举足轻重的作用.针对磁性金属矿床深部资源勘查中存在的有效信号弱、矿体形态复杂等技术难题,本文利用井中三分量磁测资料和地面磁测资料直接反演地下空间磁化率的分布情况,充分发挥纵、横向分辨率高的优势,通过磁化率的变化来确定矿（化）体的范围.分析了磁化率和磁场H_ax、H_ay、Z_a、ΔT各参量之间的关系,导出了磁化率反演方程;将钻孔编录、磁化率测井等结果作为约束条件构建了目标函数.以组合倾斜板状体模型为例,验证了方法的有效性,并在青海野马泉铁多金属矿区开展应用试验,试验结果与已有地质资料基本一致.该方法有效抑制了因反演数据源单一、缺乏约束条件而造成反演结果的多解性,实现了磁测资料的精细反演解释,为查明地下矿体空间位置、形态及规模,提高找矿效果,提供了一种新的解决方案.     

A new technique of interpreting gravity and vertical magnetic anomalies due to infinite horizontal cylindrical ore bodies

Summary A novel method of interpreting gravity and magnetic anomalies is presented here. Two diagrams of master curves, one for gravity and the other for vertical magnetic interpretation are presented. They are useful in calculating the depth of burial and the radius of infinite horizontal cylindrical bodies irrespective of density contrast or strength and direction of magnetisation. This method also enables us to infer the direction of magnetisation, resultant intensity of magnetisation and from them the susceptibility contrast; and density contrast. ThoughHenderson [8]²) reported the applicability of continuations and derivatives in the interpretation of magnetic anomalies for the first time, the authors have treated them in a more exhaustive manner in the present paper to obtain valuable relationships.