地磁总场测量阵列对磁性目标定位的方法

提出了用地磁总场对磁性目标定位的方法.以远场磁偶极子原理构建了地磁总场的测量阵列,设计了阵列对目标的定位算法,求解中采用二重梯度算法消除了地磁总场随时间变化和空间分布不均匀对目标定位的影响.实验验证了这种二重梯度算法的有效性.利用MATHCAD对定位算法中的非线性方程组进行了仿真,采用LINGO软件求解了阵列方程组.在松花江沿岸某地,用光泵磁传感器阵列对航行中的目标船进行了定位实验,所测磁场数据曲线与理论曲线完全吻合,证明了定位算法物理模型的正确性.提出了改进的磁场梯度偏置定位算法,进一步提高了目标定位精度.实验证明了这种地磁总场目标定位方法的可行性.     

磁总场异常垂直梯度三维相关成像

本文提出了磁总场异常垂直梯度三维相关成像方法,用于成像地下等效磁源分布.它首先将地下待成像空间划分为三维规则网格,然后直接计算每个网格节点磁偶板子在观测面理论磁总场异常垂直梯度与实测磁总场异常垂直梯度的互相关,其相关系数值表征等效磁偶极子分布(即磁偶极子发生的概率).理论长方体组合模型数据和实际某矿区磁测资料试验结果表明本文方法计算得到的相关系数值能基本反映地下的磁源分布,且分辨率明显高于磁总场异常三维相关成像的分辨率,也高于基于熵滤波分离异常的磁总场异常三维相关成像的分辨率.     

Magnetometer array studies

A review of experiments in which arrays of recording magnetometers deployed over an area of land to study induction problems is presented. Intensive array activities take place on a continual basis in North America, Australia, Scotland and Africa and array studies have been conducted in India, Scandinavia, and Russia. The main results are summarised. Analysis, presentation and interpretation of array data have also enjoyed significant developments which are discussed with illustrations. The use of multi-techniques in the analysis of array data has improved the reliability of the interpretational inferences. One-dimensional conductivity profiles can now be deduced by the simple inversion of the inductive response functionC(w, t) estimated from the vertical field Z and spatial gradient of the horizontal field components (ΔX/Δx+ΔY/Δy) of magnetometer array data with large gradient. Bounds can be placed on accepted profiles using the Monte Carlo process just as it is done in magnetotelluric data inversion. The results from array studies continue to improve our understanding of the synthesis of realistic tectonic models of the continents. The structure under some geothermal zones are now known through a number of recent studies.     

A geophysical survey on the archaeological site of Mugardos (NW Iberian peninsula)

Archaeological indications near Mugardos (Ferrol, NW Spain) suggest the existence of a Roman settlement. In fact, in the area were found pavements, walls with north-south and east-west orientations and some structures that endured heating. These remains are covered by soil, more than 1 m thick, and lie over schists. In order to determine the archaeological potential of the area and to delimit future excavations a geophysical survey, consisting of a joint resistivity and magnetic survey, was planned and carried out. The square array of electrodes was used and the data are discussed as apparent square array resistivity maps and azimuthal inhomogeneity ratio (AIR) maps. The magnetic survey included total field measurements using sensor heights of 0.30 and 2.30 m above the ground, so that a magnetic gradient could be computed.A combined interpretation of both resistivity and magnetic data is discussed. Later excavations have confirmed the geophysical interpretation.     

Improved signal discrimination in tectonomagnetism: Discovery of a volcanomagnetic effect at Kilauea,Hawaii

Cancellation of extraterrestrial magnetic disturbances by taking simple differences between total field readings at spaced stations is imperfect. It is shown that improvement is possible when three component observatory data are available from a single station in the general, but not necessarily immediate, vicinity of an array of total field stations used in a tectonomagnetic study. The local effects of a magnetic disturbance field depend upon its orientation, so that local field differences are more effectively generated by certain orientations of the disturbance field. The orientation of the disturbance field which correlates best with a local difference field is determined by a least-squares method, so that the correlated vector signal can be routinely subtracted from the difference field record. Application of the technique to daily averages of records from three synchronised proton magnetometers on Kilauea volcano reveals a 1.5-nT change in the local field at the time of a flank eruption in May, 1973. This effect was obscured by noise in the raw difference field data.     

Location and depth estimation of point-dipole and line of dipoles using analytic signals of the magnetic gradient tensor and magnitude of vector components

The magnetic gradient tensor (MGT) provides gradient components of potential fields with mathematical properties which allow processing techniques e.g. analytic signal techniques. With MGT emerging as a new tool for geophysical exploration, the mathematical modelling of gradient tensor fields is necessary for interpretation of magnetic field measurements. The point-dipole and line of dipoles are used to approximate various magnetic objects. I investigate the maxima of the magnitude of magnetic vector components (MMVC) and analytic signals of magnetic gradient tensor (ASMGT) resulting from point-dipole and line of dipoles sources in determining horizontal locations. I also present a method in which depths of these sources are estimated from the ratio of the maximum of MMVC to the maximum of ASMGT. Theoretical examples have been carried out to test the feasibility of the method in obtaining source locations and depths. The method has been applied to the MMVC and ASMGT computed from the total field data over a basic/ultrabasic body at the emerald deposit of Socoto´, Bahia, Brazil and buried water supply pipe near Jadaguda Township, India. In both field examples, the method produces good correlations with previous interpretations.     

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.     

Application of Hilbert‐like transforms for enhanced processing of full tensor magnetic gradient data

Commonly, geomagnetic prospection is performed via scalar magnetometers that measure values of the total magnetic intensity. Recent developments of superconducting quantum interference devices have led to their integration in full tensor magnetic gradiometry systems consisting of planar‐type first‐order gradiometers and magnetometers fabricated in thin‐film technology. With these systems measuring directly the magnetic gradient tensor and field vector, a significantly higher magnetic and spatial resolution of the magnetic maps is yield than those produced via conventional magnetometers. In order to preserve the high data quality in this work, we develop a workflow containing all the necessary steps for generating the gradient tensor and field vector quantities from the raw measurement data up to their integration into high­resolution, low­noise, and artefactless two‐dimensional maps of the magnetic field vector. The gradient tensor components are processed by superposition of the balanced gradiometer signals and rotation into an Earth‐centred Earth‐fixed coordinate frame. As the magnetometers have sensitivity lower than that of gradiometers and the total magnetic intensity is not directly recorded, we employ Hilbert‐like transforms, e.g., integration of the gradient tensor components or the conversion of the total magnetic intensity derived by calibrated magnetometer readings to obtain these values. This can lead to a better interpretation of the measured magnetic anomalies of the Earth's magnetic field that is possible from scalar total magnetic intensity measurements. Our conclusions are drawn from the application of these algorithms on a survey acquired in South Africa containing full tensor magnetic gradiometry data.     

考虑航磁水平梯度变化的ΔT网格化方法研究

实测航磁横向水平梯度反映垂直于测线的磁场梯度,比传统航磁ΔT数据包含测线之间更多的磁场信息.针对航磁数据网格化问题,采用Hardwick提出的方法,利用航磁水平梯度与ΔT数据构建拟测线,并结合Akima插值法,开展了双方向测线型ΔT网格化方法研究,最终实现了考虑航磁水平梯度变化的ΔT网格化;针对网格化结果中的虚假异常采取了有效滤波方法.通过理论模型数据和实际数据网格化处理,表明该方法可以突出航磁测线之间的异常细节、更清晰地反映线性构造或磁性体走向,提高了网格化的精度和分辨率.     

核磁共振测井在火成岩地层应用的适应性分析

火成岩岩石具有高磁化率特征,与沉积岩差别明显,岩石的高磁化率会对核磁共振信号产生显著影响,开展核磁共振测井在火成岩地层应用的适应性分析研究,为广泛使用核磁共振测井用于火成岩储层评价提供参考.采用理论模拟、实验分析及实际测井资料相结合的方法,分析了高磁化率岩石的核磁共振响应特征.分析结果表明,火成岩岩石具有很高的磁化率,从酸性火成岩到基性火城岩,岩石的磁化率一般是逐渐增大的.高磁化率岩石孔隙内部会产生强梯度磁场,孔隙内部的磁场梯度越强,核磁共振的T2谱前移越明显,核磁信号衰减幅度也越大,核磁分析孔隙度与常规分析的孔隙度误差也越大;相同岩石磁化率情况下,岩石孔径越小,岩石孔隙内部磁场梯度越强,核磁信号衰减越大,核磁共振测井计算的孔隙度也越低.因此,核磁共振测井与岩性有关系,核磁共振测井会受到高磁化率岩石的显著影响.核磁共振测井在部分火成岩,尤其在中基性火成岩和小孔径火成岩储层中应用具有很大的局限性.     

不同方位甚低频电磁波衰减的初步研究*

利用2005—2010年DEMETER卫星记录NWC发射站的19.8kHz电场功率谱数据,采用统计和线性拟合等方法,研究了NWC站发射的电磁波在顶部电离层及磁共轭区激发的电场效应及其在不同方位的衰减特性。结果发现:(1)在NWC发射站上空,卫星记录电场呈椭圆状扩散分布,电场最强中心点位置相对地面人工源位置有一定偏移;(2)研究区上空电场北部衰减梯度最小,东部衰减梯度最大;(3)在其磁共轭区,南部电场衰减梯度最小,北部衰减梯度最大。综合分析认为人工源(NWC)发射的VLF电磁波传播到电离层高度后,受地磁场影响,电磁波主要沿磁力线方向传播,空间电场最强中心点位置相对发射站的位置发生偏移,向磁赤道方向倾斜,而且偏于磁赤道方向的电场衰减梯度最小。     

大连九龙地磁台场地梯度的跟踪测量分析

介绍了大连九龙地磁监测中心建设过程中场地磁场梯度的跟踪测量情况。根据中国“地震及前兆数字观测技术规范”电磁部分的要求,在地磁台站建筑施工过程中,需要用质子旋进式磁力仪在施工现场对场地磁场背景进行监测。     

REMARKS ON THE USE OF THE MAGNETIC GRADIOMETER IN OIL EXPLORATION *

The principle of optical pumping allows the design of magnetometers with high sensitivity. When mapping the field of the total magnetic force, it is no longer possible to make full use of the accuracy of the instruments because the accuracy of the reduction of the diurnal variation is limited. By simultaneously recording with two instruments in different altitudes, the vertical gradient can be measured which doe's not depend on the time variation of the magnetic field. Therefore, the gradiometer seems to be a more adequate tool for oil exploration than the magnetometer. It is investigated in this article whether the results of this gradiometer or the measurements of the total magnetic field by the high sensitive magnetometer are more useful in oil exploration. The article comes to the conclusion that for most problems of oil exploration the total magnetic field is a more valuable unit than the vertical gradient measured directly by the gradiometer. The total magnetic field allows a better investigation of the tectonics than the vertical gradient. The apparent advantages of the gradiometer claimed by its supporters are mainly based on inconsistent mathematical concepts.     

航磁垂直梯度调整ΔT水平方法研究

使用中国自主研发、具有完全自主知识产权的航磁全轴梯度勘查系统首次试飞获得的资料,利用航磁垂直梯度数据进行ΔT水平调整研究,讨论了航磁及航磁梯度资料水平调整的原则,提出梯度资料长波水平的调平方法.对经长波水平调整后的垂直梯度数据进行位场转换,获得不受地磁场随时间变化影响的磁异常,并通过恢复磁异常的长波信息,实现了ΔT的水平调整和增强.实际资料处理表明调平过程具有物理意义、效果明显.     

田庄等测点周围的磁场梯度测量

磁场梯度测量是选建地磁测点过程中的一项重要工作。在田庄、灵山、东兴庄测点周围，进行了磁场梯度的一般测量与详细测量，结果显示，各测桩的水平方向与垂直方向的磁场梯度均不大于１．０ｎＴ／ｍ。完全符合磁测规范的要求。     

A simple method for determining the vertical growth rate of vertical motion at the top of earth's outer core

A simple new method is described for extracting, from magnetic observations taken at Earth's surface, the vertical growth rate of vertical motion, ?u/?r, at special isolated points on the top surface of Earth's liquid core. The technique utilizes only the radial component of the frozen-flux induction equation and it requires information only on the radial magnetic field, B_r, its horizontal gradient, and its secular variations, ?B_r/?t, at the core-mantle boundary.     

2021年5月22日玛多MS7.4地震震前岩石圈磁场变化特征分析

本文利用2019年和2020年两期的全国流动地磁矢量数据,针对玛多M_S7.4地震震中附近的岩石圈磁场空间分布,分析研究了岩石圈磁场各个分量的震前变化特征。结果表明,玛多M_S7.4地震发生前,震中附近的岩石圈磁场分量均发生了不同程度的变化:震中位于各分量的弱变化区域和零变线附近;震中位于各分量的高梯度带和低梯度带之间。本文证实了地震发生前会引起岩石圈磁场变化这一现象,并总结了震前地磁分量的变化特征,为今后强地震（尤其是M≥7.0地震）的震磁关系研究提供了一个震例参考。      

Applying the tilt‐depth and contact‐depth methods to the magnetic anomalies of thin dykes

A new method for the calculation of the depth, location, and dip of thin dykes from pole‐reduced magnetic data is introduced. The depth can be obtained by measuring the distance between chosen values of a tilt angle that is based upon the ratio of the magnetic field and its Hilbert transform over the dyke. Alternatively, it can be obtained from the horizontal derivative of the ratio of the Hilbert transform of the field to the field itself, over the dyke. The latter method also allows the dip of the dyke to be estimated from the gradient of the depth estimates.     

Depth and shape estimates from simultaneous inversion of magnetic fields and their gradient components using differential similarity transforms

The differential similarity transform of a magnetic anomaly is a linear combination of its intensity and gradient components. This transform is sensitive to the distance between a chosen central point of similarity and the source and depends on the degree of homogeneity of the field. Taking advantage of this property, a new field inversion method resulting in the evaluation of source position and shape type is proposed and implemented. The field gradient components are measured directly in magnetic gradiometry, or they can be calculated from the measured field data. Regional and local linear backgrounds are accounted for by the method. The method can be applied on either regularly or irregularly-spaced data sets, on even or uneven surfaces of observation. The solving of the systems of equations is not necessary. A semi-automated inversion for both location and shape of the sources is implemented. Model and field tests illustrate the effectiveness of the proposed inversion technique for depth and shape estimates.