地磁正常场的选取与地磁异常场的计算

根据2003年中国地磁观测数据（包括135年地磁测点和35个地磁台）以及我国邻近地区38个IGRF计算点的地磁数据,计算中国地磁异常场的分布。选取两种地磁场模型作为地磁正常场,一是国际参考地磁场的球谐模型,二是中国地磁场泰勒多项式模型。根据各个测点的地磁异常值（观测值减去模型计算值）,用球冠谐分析方法计算地磁异常场的球冠谐模型,并绘制2003年中国地磁异常（△D,△I,△F,△X,△Y,△Z）。分析和讨论了中国地磁异常场。     

Methods for reconstructing harmonic functions from the magnetic field ΔT and V.N. Strakhov’s function ΔS: A review

We give an overview of the methods designed for reconstructing close-to-harmonic functions from the magnetic field ΔT. The formula of Yu.P. Tafeev is refined. It is shown that this refined formula directly leads to the relation derived by V.M. Gordin and his colleagues that allows isolating the harmonic component in the function ΔT. V.N. Strakhov’s linearized representation of the function ΔT is immediately derived from the main approximate Tafeev formula for Q ΔT. The experience of using Strakhov’s ΔS function in the interpretation of the magnetic anomaly ΔT generated by the Krivoi Rog structure is described. It is noted that the problem of reconstructing the corresponding harmonic functions from the data of magnetic and gravity surveys has much in common. The specific features of measuring the magnetic field H and magnetic induction B in the material media are considered, and the physical interpretation of these fields is presented.     

基于等效源的总强度磁异常非线性处理方法

总强度磁异常(ΔT1)常规处理方法通常将其近似当作磁异常矢量在地磁正常场方向上的投影(ΔT2).然而对于高磁环境如磁铁矿处的磁异常场幅值可达104 nT甚至更大的情况,上述近似假设则不再成立,如采用常规处理方法可能会带来明显的误差.本文针对该问题,提出一种基于等效源的总强度磁异常非线性处理方法;该方法根据总强度磁异常获取流程,直接反演实测地磁场总强度幅值与正常场强度幅值之差来求取等效场源.本文首先通过模型试验,分析ΔT1与ΔT2的差异;然后将ΔT1当作ΔT2采用常规处理方法以及ΔT1采用新方法得到的结果作对比分析,结果表明新方法处理结果与理论值的差异为常规方法处理结果的1/5甚至更小,充分说明了新方法的有效性;最后将该方法应用于铁矿实测总强度磁异常处理实例中来转换计算磁异常总模量,其实际应用效果进一步体现了高幅值总强度磁异常数据处理过程中采用新方法的必要性.     

中国地区地磁场的球冠谐和分析

本文根据中国及其邻近地区的地磁三分量绝对测量资料,利用球冠谐和分析方法,计算出中国地区地磁剩余场的冠谐模型,地磁剩余场△X,△Y,△Z的模型均方根偏差分别为106.9,89.7,137.6 nT.提出由地磁场的国际参考地磁场和地磁剩余场的冠谐模型组成的联合模型作为中国参考地磁场的模型,它能较好地表示中国地磁场的分布.以地磁场的联合模型为正常背景场,计算出三分量磁异常的冠谐模型,并分析了磁异常的基本特征,它将为深入研究中国岩石层结构提供新证据.     

Analysis of the lithospheric magnetic anomalies and tectonics in continental China and the adjacent regions using CHAMP satellite data

In this paper, the CHAMP satellite data from the last decade are used to map the lithospheric magnetic field in continental China at altitudes of 300 and 400 km using the spline function model technique. The results show that our model has higher amplitude than the spherical harmonic magnetic field model MF7. The resulting lithospheric anomaly distribution of continental China agrees with the geotectonic structure and crustal thermal structure. The magnetic anomaly distribution clearly indicates extrusion of the Tibetan Plateau towards the northeast and northwest directions. Northward movement is impeded by the relative rigid Tarim Basin featured with positive magnetic anomalies. The positive magnetic anomalies in Tarim basin extend to the northern Tibetan Plateau, suggesting that the southern part of Tarim Basin underthrusted the northern Tibetan Plateau.     

The calculation of perturbation and induction arrows for a three-dimensional conductivity model and dipole source fields

A numerical method is used to calculate the electromagnetic fields associated with a three-dimensional conductivity anomaly. The source field is due to horizontal magnetic dipoles placed at two different positions with respect to the conductivity anomaly. The transfer functions and related perturbation and induction arrows associated with the fields are calculated and compared with the arrows obtained from a uniform source calculation. The results show the source effect on the induction arrows and indicate that the perturbation arrows provide a method of outlining the spatial extent of the anomaly. The transfer function calculations are made for both exact and approximate normal fields. In the transfer function calculation the anomalous fields are correlated with a normal field as suggested by Schmucker (1970) and Cochrane and Hyndman (1970).     

Accuracy estimates for some global analytical models of the Earth’s main magnetic field on the basis of data on gradient magnetic surveys at stratospheric balloons

Two global analytical models of the main magnetic field of the Earth (MFE) have been used to determine their potential in deriving an anomalous MFE from balloon magnetic surveys conducted at altitudes of ～30 km. The daily mean spherical harmonic model (DMSHM) constructed from satellite data on the day of balloon magnetic surveys was analyzed. This model for the day of magnetic surveys was shown to be almost free of errors associated with secular variations and can be recommended for deriving an anomalous MFE. The error of the enhanced magnetic model (EMM) was estimated depending on the number of harmonics used in the model. The model limited by the first 13 harmonics was shown to be able to lead to errors in the main MFE of around 15 nT. The EMM developed to n = m = 720 and constructed on the basis of satellite and ground-based magnetic data fails to adequately simulate the anomalous MFE at altitudes of 30 km. To construct a representative model developed to m = n = 720, ground-based magnetic data should be replaced by data of balloon magnetic surveys for altitudes of ～30 km. The results of investigations were confirmed by a balloon experiment conducted by Pushkov Institute of Terrestrial Magnetism, Ionosphere, and Radio Wave Propagation of the Russian Academy of Sciences and the Moscow Aviation Institute.     

Evidence of electrical anisotropic structures in the lower crust and the upper mantle beneath the Rhenish Shield

We present the results of a long-period electromagnetic investigation of the crust and upper mantle beneath the Rhenish Shield in Western Germany. The magnetotelluric phase data reveal a frequency-dependent regional strike that varies only smoothly across the array. At short periods (t<100 s) the striking of the maximum phase splitting is N45°E, which can be explained with an electrical anisotropic lower crust. At long periods (t>1000 s) there is a consistent striking in WE direction, which provides strong evidence of an anisotropic structure in the upper mantle, too. Geomagnetic data were also used, whereas we reference the magnetic field components of all sites to an arbitrarily chosen field site. That provides a direct view of the anomalous current flow. We show that in case of a non-one-dimensional electrical substructure of the reference site all magnetic transfer functions of the other sites can be affected by lateral conductivity contrasts beneath the reference site. A simple method to remove such effects is introduced. Applying this method the magnetic data show a distinct anomaly in the northwestern part of the area. Finally, we present a 3D model of the conductivity structure beneath the Rhenish Shield which contains a superposition of two anisotropic structures in the lower crust (2000 S) and the upper mantle (20?000 S), respectively, and a local anomaly in the upper crust (4000 S).     

中国及邻区卫星岩石圈磁场与近10年地震活动性

利用NGDC-720全球卫星岩石圈磁场模型给出的球谐系数,计算得到中国及邻区的岩石圈磁场分布;并与该地区活动块体及近10年（2003-2012年）MS 5．0以上地震活动性进行对比,发现岩石圈磁场同活动块体分布及地震活动性存在较为密切的对应关系。该相关性有望在未来地震研究和监测中发挥作用。     

基于印模法的地面磁电阻率数据三维非线性共轭梯度反演

本文提出了能提高异常体分辨能力,同时得到绝对电导率的地面磁电阻率数据三维反演方法.磁电阻率响应用准直流的低频磁场代替;数值模拟由频率域电场满足的Helmholtz方程出发,采用三维交错网格有限差分法;长直导线源作为发射源,其中源的计算包含在背景场中;结合地面磁电阻率数据各分量的特点,选择y分量进行反演研究;反演采用三维非线性共轭梯度反演技术,为了提高异常体的深度分辨能力,进行迭代重构反演;用印模法对初始模型进行重构,采用的是辅模型在浅部,元模型在深部的组合方式.从合成数据和实际数据的反演结果可以得到以下的认识:(1)由频率域麦克斯韦方程组出发,低频磁场数据反演可以直接得到电导率,而不是相对电导率之比;(2)采用印模法组合初始模型,进行迭代重构反演,可以提高地面磁电阻率数据反演对异常体的分辨能力,确定埋深位置,同时不会丧失对于浅部异常体的分辨能力;(3)在结合印模法的地面磁电阻率数据三维反演中,深部异常体的分辨能力受地表不均匀导电体影响较小;(4)确定印模深度可以采用上一次重构反演结束时的模型变化量,通过相邻两次重构反演结束时的模型变化量之差来确定迭代重构是否终止.因为静磁场与重力场在数学上的相似性,本文的反演方法可以被运用到重力场等位场的地面数据的反演中.     

An analysis of geomagnetic response functions prior to the Tohoku,Japan earthquake

Geomagnetic records from 20 Japanese observatories have been used to yield time series of response function (RF) components for 20 years at periods of between 2.5 and 60 min. Six observatories showed anomalous variations lasting 3–5 years in the short period part of the above range of periods prior to the March 11, 2011 Tohoku earthquake. The variations could have been intermediate-term precursors. We made a detailed analysis of how noise affects the results using coherence criteria, visual control, and the remote-reference technique. We clarified the conditions that make response functions dependent on geomagnetic activity. For 19 observatories we constructed the tensor of the anomalous magnetic field with Kakioka as the base site. An anomaly in electrical conductivity striking WNW–ESE has been identified beneath the Boso Peninsula near Tokyo in the conditions of strong noise. We sought to corroborate the reality of the anomaly by visual control and processing of nighttime records with minimum noise. We advanced idea that precursors can be monitored using the DC noise field in the presence of a shallow conductivity anomaly. We provided a tectonic interpretation of the obtained RF anomalies. The Boso conductivity anomaly is interpreted as being due to a graben-shaped structure of the sediments and possibly to a deeper plate-tectonics structure, that is, the Sagami Trough. We examine similarities and differences between the Boso anomaly and the Avacha anomaly in Kamchatka, and provided recommendations for further study of the Boso anomaly and for using the Avacha anomaly to monitor EM precursors in Kamchatka.     

Depth Estimates for Slingram Electromagnetic Anomalies from Dipping Sheet-like Bodies by the Normalized Full Gradient Method

The Normalized Full Gradient (NFG) method was proposed in the mid–1960s and was generally used for the downward continuation of the potential field data. The method eliminates the side oscillations which appeared on the continuation curves when passing through anomalous body depth. In this study, the NFG method was applied to Slingram electromagnetic anomalies to obtain the depth of the anomalous body. Some experiments were performed on the theoretical Slingram model anomalies in a free space environment using a perfectly conductive thin tabular conductor with an infinite depth extent. The theoretical Slingram responses were obtained for different depths, dip angles and coil separations, and it was observed from NFG fields of the theoretical anomalies that the NFG sections yield the depth information of top of the conductor at low harmonic numbers. The NFG sections consisted of two main local maxima located at both sides of the central negative Slingram anomalies. It is concluded that these two maxima also locate the maximum anomaly gradient points, which indicates the depth of the anomaly target directly. For both theoretical and field data, the depth of the maximum value on the NFG sections corresponds to the depth of the upper edge of the anomalous conductor. The NFG method was applied to the in-phase component and correct depth estimates were obtained even for the horizontal tabular conductor. Depth values could be estimated with a relatively small error percentage when the conductive model was near-vertical and/or the conductor depth was larger.     

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

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

Regional spatiotemporal variations of a nondipole magnetic field over the Chinese mainland and neighboring regions in millennial scale

Spatiotemporal variations of the nondipole (ND) magnetic field over the Chinese mainland and neighboring regions from 10000 BC to 1990 AD were analyzed using the latest global geomagnetic models CALS10K.1b, CALS3K.4, and IGRF11. Moreover, for field sources, we investigated 2 ⁿ (n = 2–10) pole ND fields and their energies. The results suggest that the study period can be divided into three. The intensity of the ND field has been mainly positive since 10000 BC and lasted almost 7500 years, then gradually decreased to negative in 2500 BC to 1500 AD, and finally sharply increased to positive. The anomaly areas of the ND field in East Asia took shape for n = 3, when the anomaly areas in East Asia were shaped into closed circles in the mainland. This suggests that the first three harmonic degrees account for most of the ND field. The energy of the ND field rapidly attenuates at the core–mantle boundary and is stable at the surface.     

地磁匹配用于水下载体导航的初步分析

随着地磁场模型的日趋完善和滤波技术的不断成熟,地磁匹配方法在近几年得到迅速发展,它一般用于近地卫星,还未用于海洋环境.针对这一问题,本文结合地磁场特点和海洋地磁测量的现状,分析了水下地磁匹配的可行性.分别从匹配的地磁模型和测量值中存在干扰磁场两方面,指出了制约该技术发展的几个因素.在此基础上,提出了地磁用于水下载体定位的两条技术途径,一是用小波变换滤去地磁异常,与大尺度地磁模型匹配定位;另外,也提出了利用地磁异常,结合同时定位与构图(SLAM)算法实现自主导航.     

河北省数字地电场数据分析

通过分析总结河北省数字地电场的正常场和异常信息，发现正常记录除了显示日变周期外，还记录到了半日和三分之一日周期，比较长的周期还有半月周期等；干扰信息主要与电极极化及磁暴等有关，大磁暴干扰一般会出现高频信号，破坏正常日变形态，其主要表现为多台同步等幅信号；对台站的电场方向性研究发现，在正常时段电场记录方向性显著，而磁暴记录期内方向性较差，显示外空异常源方向比较复杂。最后对一些可能与地震有关的异常现象进行了分析探讨。我们认为，目前的多极距地电场观测对排除干扰及异常识别带来很大方便，但共用电极的观测方式为异常识别增加了障碍。     

Changes of lithospheric magnetic anomalies with altitude (According to the CHAMP satellite)

Maps of the magnitude of the full vector and the vertical component of an anomalous lithospheric magnetic field over the Voronezh anticline (VA) for the three high-altitude observation levels were compiled based on geomagnetic measurements from the CHAMP satellite. The isometric positive anomaly centered at about 50° N and 37° E stands out. Its amplitude decreases with increasing observation altitude without changing the form. Comparison of the parameters of the detected anomaly with data obtained for this site by other methods confirms that it really exists and that its spatial position is accurately determined, which indicates the reliability of the values of the selected field of lithospheric anomalies. The change in the parameters of the magnetic anomaly with respect to the observation level over the Earth’s surface is consistent with the concepts of geological structural features of the lithosphere in the region. The anomaly offset to the south on the satellite altitudes apparently indicates an uplift of crystalline basement rocks and a more southern position of VA deep roots relative to that accepted in the global magnetization model. The use of satellite data obtained at different altitudes opens up additional possibilities for the application of gradient methods in the interpretation of the magnetic fields of lithospheric anomalies.     

位场异常分离的递减半径迭代法

本文提出了位场区域-剩余异常分离的空间域递减半径迭代法.由给定半径的圆周上八点位场值的算术平均导出一个新的八点圆周平均公式,它是一个由重磁数据计算区域异常的滤波器.该滤波器的传递函数有一主瓣和多个旁瓣,且半径越大,旁瓣越多,滤波器特性越差.由从大到小不同半径对应的传递函数的连乘积构造了递减半径迭代传递函数,它以大半径为参变量.递减半径迭代传递函数类似低通滤波器,其截止波数与大半径成反比.由递减半径迭代传递函数,给出空间域分离区域异常的递减半径线性迭代法,由重磁数据减去区域异常求得剩余异常.进一步,通过构造非线性修正系数,把空间域递减半径线性迭代法中线性迭代公式变成非线性迭代公式,得到空间域递减半径非线性迭代法.重磁理论模型及安徽泥河铁矿重磁资料的试验表明,空间域递减半径非线性迭代法可有效地压制假异常和高频干扰以及减少异常畸变,并能有效地分离出区域异常和剩余异常.     

ELECTROMAGNETIC ANOMALOUS PULSES BEFORE AND AFTER EARTHQUAKES(M>6.0) IN CHINA CONTINENT IN 2008

In this paper,we provided the electromagnetic anomaly phenomena prior to five earthquakes with magnitude M_W>6.0 occurring in China continent in 2008.The electromagnetic data in frequency bands 800~0.1Hz are recorded at two stations located in Urumqi,Xinjiang Uygur Autonomous Region (URQ) and Tonghai,Yunnan Province (TH).The time series of four horizontal components of the electromagnetic field at the stations are converted to the spectra using Fast Fourier Transform.The relationship between the electromagnetic anomalous phenomena and the earthquake events is studied through analyzing the temporal variation of electromagnetic spectra and comparing them with earthquake events.This study provides the new examples for electromagnetic anomaly phenomena before the earthquakes.The following features can be found.①The obvious anomalous power spectrum density (PSD) of electromagnetic fields in frequency band of 128~0.5Hz appeared before earthquakes.The anomalous PSD is related to the earthquake magnitude and the epicenter distance.The anomaly size of magnetic PSD is about 1-3 orders of magnitude bigger than the background field before Wenchuan M_S8.0 earthquake.② The changes of parameters are also related with frequency of the signal;The ELF signal strength observed at seismic stations is significantly related with distance between the stations or the measured magnetic and electrical field components and the source of earthquake.We also found that the field strength attenuated faster in area closer to the source than in far area.It is clear that the magnetic signals have higher signal noise ratio than the electric data.③ The amplitudes are also related with the azimuth of the earthquake to station and with the frequency of the data.④ The anomalous pulses showed cluster phenomenon in time.     

3D inversion of gravity data by separation of sources and the method of local corrections: Kolarovo gravity high case study

We present a novel methodology for 3D gravity/magnetic data inversion. It combines two algorithms for preliminary separation of sources and an original approach to 3D inverse problem solution. The first algorithm is designed to separate sources in depth and to remove the shallow ones. It is based on subsequent upward and downward data continuation. For separation in the lateral sense, we approximate the given observed data by the field of several 3D line segments. For potential field data inversion we apply a new method of local corrections. The method is efficient and does not require trial-and-error forward modeling. It allows retrieving unknown 3D geometry of anomalous objects in terms of restricted bodies of arbitrary shape and contact surfaces. For restricted objects, we apply new integral equations of gravity and magnetic inverse problems. All steps of our methodology are demonstrated on the Kolarovo gravity anomaly in the Danube Basin of Slovakia.