重磁异常反演的拟BP神经网络方法及其应用

把神经网络与重磁异常反演理论相结合,提出了用于重磁反演的一种拟BP神经网络方法．基于3层神经网络结构,把隐含层神经元设定为三维空间物性（磁化强度或密度）单元．对实测与理论重磁异常经S型函数变换,采用自动修改物性单元物性值的拟BP算法,反演三维空间的物性分布．利用该网络对理论模型数据和内蒙古某花岗岩体上的航磁资料进行了反演计算,取得了满意的反演效果．     

低纬度磁异常的初始模型约束全卷积神经网络化极方法

磁异常化极是磁测数据处理的重要基础工作,低纬度磁异常化极运算不稳定.深度学习是一种有效解决不稳定问题的数据驱动方法.然而,现有基于深度学习的位场数据处理仅以位场响应作为输入,缺乏场景外样本数据的通用性,往往导致深度学习网络预测结果的泛化性能不足.因此,为了克服这一问题以及进一步提高低纬度磁异常化极的精度.受“知识&数据”联合驱动深度学习网络结构的启发,本文提出一种基于初始模型约束的低纬度磁异常化极全卷积神经网络结构.该网络结构的输入包含两部分,一部分为低纬度磁异常,另外一部分为初始模型,网络输出为垂直磁化方向磁异常.训练与测试的磁异常数据采用基于网格点格架函数的空间域快速正演生成,以及对低纬度磁异常采用稳定的频率转换来获取初始模型,另外本文对数据集中10%的随机样本加入了5%的高斯噪声,以此增强全卷积神经网络结构的鲁棒性.理论模型试验验证了本文所提方法的有效性、精确性以及鲁棒性.最后,将本文方法应用于实际磁测数据,取得了良好的效果.     

基于深度学习的重力异常与重力梯度异常联合反演

高效高精度的反演算法在重力大数据时代背景下显得尤为重要,受深度学习卓越的非线性映射能力的启发,本文提出了一种基于深度学习的重力异常及重力梯度异常的联合反演方法.文中首先提出了一种基于网格点几何格架的重力异常及重力梯度异常的空间域快速正演算法,这为本文深度学习反演算法的实现奠定了基础;其次对大量的不同密度模型进行正演计算获得样本数据集;然后设计了一种端到端的深度学习网络结构(GraInvNet),再利用样本数据对该网络结构进行训练;最后进行反演预测.组合模型试验表明,多维度数据联合反演相比单一分量反演其结果更“聚焦”,且与模型边界高度吻合,并且对于复杂模型的姿态与物性预测具有极为显著的优势,以及对于含噪声数据的反演,其质量也不会降低;Vinton岩丘实测重力数据也验证了文中方法的有效性;从而证明了深度学习在重力数据的高效高精度反演方面具有的巨大潜力.     

强剩磁强退磁条件下的二维井中磁测反演

强剩磁、强退磁改变了总磁化强度的大小和方向,给磁测资料解释带来困难.为此,本文利用二维井中磁测数据反演磁化强度矢量的二维分布.首先利用井中磁测的磁异常模量反演磁化强度大小的分布.然后,在已知磁化强度大小分布的前提下,拟合磁场分量,反演磁化强度方向的分布.其中,磁化强度大小和方向均用共轭梯度法求解,并通过预优矩阵改善磁化强度大小的反演效果.理论模拟说明,该方法能准确获得磁化强度矢量分布.磁化强度矢量反演结果包括感磁、剩磁及退磁的影响,这为研究强剩磁、高磁化率矿床提供了一种有效方法.     

重磁数据三维物性反演方法进展

综述了重磁数据三维物性反演方法中的几个关键问题.主要包括正演快速算法、反演框架、约束因子讨论、反演算法实现等方面.正演快速算法主要讨论了等效存储几何格架技术、基于GPU加速的并行计算技术以及小波压缩技术.三维物性反演则是在最小二乘意义下使目标函数达到极小的线性或非线性反演.指出,对于特定地质问题需要谨慎选择不同且合适的约束方法乃至反演算法,才能达到好的效果.最后讨论了重磁数据三维物性反演较好的应用前景及发展方向.     

数据空间磁异常模量三维反演

强剩磁的存在通常导致了总磁化强度方向未知,进而影响了磁异常的反演和解释.磁异常模量是一种受磁化方向影响小的转换量,可以在强剩磁条件下通过反演三维磁化强度大小分布来推测场源分布状态.我们提出了一种数据空间磁异常模量反演算法来减少剩磁的影响.与标准的模型空间L2范数正则化反演方法相比,我们的方法有两个优点:一是无需搜索正则化参数(需要反复求解非线性反演问题),因而可以减少计算时间;二是反演结果更加聚焦,深度分辨率更高,我们对此进行了原因分析.通过模型和实测数据测试证明了该算法的有效性和更好的反演效果.     

磁化强度矢量反演新方法及其在卡拉塔格铜多金属矿床找矿中的应用

磁异常的反演是地球物理勘探的重要手段,三维磁化率反演是磁异常定量解释中的一种重要方法.由于剩磁的存在使得磁化方向与地磁场方向产生偏差,从而影响了磁异常反演与解释的精度.本文基于磁异常模量反演和磁化强度矢量反演方法得到了一种新的磁化强度矢量反演方法.与以往的磁化强度矢量反演方法相比,该方法以磁异常模量反演得到的磁化率模型为约束,采用Lp范数正则化方法求解,提高了磁化强度矢量反演的精度和效率.本文通过模拟试验的反演计算,验证了这种磁化强度矢量反演方法的有效性.最后,将本文方法应用于新疆东天山卡拉塔格地区航磁数据的解释,获得了地下空间不同磁性差异的磁性体的空间分布特征,为进一步分析研究区隐伏矿床提供了重要信息.     

基于瞬变电磁矩变换的快速三维反演方法

瞬变电磁法的严格三维反演计算复杂、占用资源多,在普通计算机上难以实现.本文引入瞬变电磁矩变换的概念,提出一种快速三维反演方法.该方法基于阻性限制(resistive limit)特性,建立包含异常体的三维大地的一阶矩响应正演算法,根据不同约束条件,选择优化的最速下降法实现瞬变电磁快速三维反演.文中通过含异常体的三维大地正演一阶矩与仿真数据一阶矩的对比,验证了快速三维正演算法的有效性,之后在不同约束条件下,利用优化的最速下降法实现了对含噪声的仿真瞬变电磁数据的快速三维反演.结果表明,该方法能够在普通计算机上短时间内较为准确地反演出地下异常体的体积和位置,在瞬变电磁数据的实时解释工作中具有良好的应用前景.     

基于方向tilt-Euler的三维磁数据快速反演

三维磁异常快速反演是磁法勘探进行地质解释的一项重要内容.本文在方向tilt梯度基础上进行导数处理,并进行合理改进,得到了6个改进的方向tilt梯度导数,丰富了不同方向上的磁异常信息.对常规欧拉反褶积方程进行三个方向求导,所构成的方程组与6个改进的方向tilt梯度导数进行结合,从而构建出了两个基于方向tilt梯度导数的快速反演方程组(方向tilt-Euler法),用于反演场源的位置参数.另外,将三个方向导数下的欧拉反褶积进行平方和及均方根处理,得到了场源构造指数的求解公式.单一模型试验表明,6个改进的方向tilt梯度导数异常基本不受磁化方向影响,可从不同方向上反映磁异常信息,方向tilt-Euler能够准确地反演出场源的位置及构造指数.叠加模型试验证实了方向tilt-Euler法相对于tilt-Euler法,具有反演解聚集度更强、连续性更好、反演准确度更高、受噪声影响更小及无虚假反演解等优点.将本文方法应用于内蒙古塔木素航磁资料处理中,获得了丰富的地下磁源分布信息,为了解该地区隐伏断裂构造及岩体分布等提供了有效依据.     

南海北部磁异常特征及对前新生代构造的指示

为了研究南海北部前新生代构造,利用新近的船载磁力测量数据,对磁异常进行变纬度化极,并反演计算视磁化强度和磁源重力异常,以及对三条OBS剖面进行重磁拟合.结果认为东沙隆起高磁异常带是浙闽沿海火山岩带向西的延续,其间被NW向古老的转换边界断裂F10错断;NE向的F2断裂是高磁异常带的南界,并限制了底侵活动的北界;F3断裂在...     

Optimal processing of satellite-derived magnetic anomaly data

The statistical concept of the power spectrum has proven to be of great value in the analysis of time series and linear systems for which the inputs and outputs are functions of time. This paper shows how the concept can be extended to two-dimensional spatial power spectra and illustrates, by example, how the concept can be applied to the determination of optimal data processing methods for satellite-derived magnetic anomaly data and to the planning of missions to obtain such data.The analysis techniques indicated are applied to a data set and data processing procedure described by Mayhew et al. (1980). These authors describe magnetic anomaly data for Australia and surrounding ocean obtained by the polar orbit POGO series satellites. This paper shows that the data processing method used by these authors is approximately equivalent to an invariant two-dimensional linear filter and that it is reasonably close to optimal with respect to accuracy, though some possible improvements are suggested. Nevertheless, as is usual when filtering data, some real “signal” is unavoidably removed along with the “noise”, resulting in errors that can be quite large.A method for reducing these errors by using additional data from a medium inclination orbit satellite (for example, 60° inclination) is suggested.     

关于磁异常的模拟试验研究

通过对不同深度、不同形态、不同倾斜状态的已知模型引起的磁异常垂直分量的测量，进行不同侧面的数据处理，从中探索出不同侧面的数据处理在判断磁异常模型各特征参数中所起的作用，最后总结出测量深度、形形倾斜以及磁化方向对磁异常分布的影响，从中得到了一些有益的结论。     

Spline approximation for TALWANI's method of magnetic anomaly computation

When using TALWANI's method for computing magnetic anomalies of three-dimensional bodies unexpected errors may occur. They are caused by the numerical depth integration. It is shown how those errors can be recognized and reduced by employing spline approximations.     

A study of the Rima Sirsalis lunar magnetic anomaly

The source of the lunar magnetic anomaly associated with the Rima Sirsalis linear rille has been modelled using the vector field intensities due to arbitrary uniform magnetization in a rectangular prism. It is shown that in order to match the Apollo 16 subsatellite data, the lunar surface near the rille must have a vertical magnetization of 6–9 × 10⁻³ G if the anomaly is due to flux leakage from a gap in the crust with the dimensions of the rille. This is more than one order of magnitude larger than the magnetization of any lunar sample, but is comparable with the high magnetization recently deduced for the Reiner γ formation in Oceanus Procellarum. An alternative explanation is that Rima Sirsalis and its surroundings are the site of a vertical magnetization contrast of 10⁻⁵ – 10⁻⁴ G which is at least as wide as the rille and extends to a depth of tens of kilometers in the crust. A wider magnetic source reduces the required magnetization (or depth) proportionately, since to first order the field at high altitude is proportional to the magnetic dipole moment per unit length.     

磁异常△T三维相关成像

本文将重力和重力梯度数据三维相关成像方法推广到磁力勘探领域,推导并建立了磁异常△T的三维相关成像方法,同时提出了基于熵滤波分离异常的三维相关成像算法来提高成像分辨率.组合模型磁异常△T数据和实际磁测资料试验分析表明,本文方法能成像出地下地质体的空间赋存状态和等效磁性分布,具有良好的横向和纵向分辨率.     

非洲磁异常对地磁场结构及其长期变化的控制作用

地球非偶极磁场在主磁场结构及其长期变化中起着重要作用.非偶极磁场主要表现为行星尺度磁异常,它们是南大西洋磁异常、非洲磁异常、欧亚大陆磁异常、澳洲磁异常和北美磁异常.在这5块磁异常中,非洲磁异常对磁赤道的形状和位置以及全球长期变化特征有极大的影响.非洲磁异常的重要性主要表现在3方面：第一,由于异常区位于赤道这一特殊的地理位置,所以它极大地影响磁赤道的形状和位置.相对于偶极场的地磁赤道而言,异常区所在的中北非洲和中大西洋地区的磁赤道向北移动,最大移动量可达约15°.第二,非洲磁异常的快速西漂对全球长期变化的分布起着决定性作用,它在该异常区西边的中美洲形成了全球最主要的长期变化区,在1900~2005年期间,最大年变率Z_max超过200 nT/a.第三,非洲负磁异常区与其南面的南大西洋正磁异常区相结合,〖HJ〗它们的变化使西半球地磁场强度大大减弱,也使全球磁场发生显著畸变.这两块磁异常与深部的反极性斑区有着成因联系.     

墨西哥及邻近地区地壳磁异常场研究

基于墨西哥国内49个地磁测点的3分量数据,利用全球模型CM4,区域模型Taylor多项式及曲面Spline 模型计算并分析了地壳磁异常场,并在境外均匀选取了16个CM4模型的地壳场值作为补充点以控制边界效应.为了进一步比较验证,还选用了最新的EMM2010模型.结果显示将实测值减去CM4地核场的差值所得到的磁异常分布与直接基于测点的CM4磁异常分布在分布形态和强度都有较大的差别,而与基于实测值的Spline网格值减去CM4地核场网格值后非常类似.总强度F在境内以负值居多,强度随纬度的升高其强度由大变小再变大;水平分量H在境内以正值居多,强度随纬度的升高而强度变大;磁偏角D分量除了西北部分区域外,其它区域都为正值,强度随经度的升高而变大.基于现代高精度的CM4及EMM2010模型绘制了墨西哥地区的磁异常网格值进行比较,两者间存在的差异主要由截断阶数与所使用卫星数据的不同所致.通过实测值结合Taylor多项式和曲面Spline两种区域模型研究磁异常,认为若进行地磁场异常的理论研究,可采用现有的高精度全球模型绘制的磁异常图,若进行实际的工程应用,可从实际的测点出发,结合常用的区域模型进行应用,曲面Spline模型较适于模拟区域地区的地壳磁异常.     

烃渗漏引起的磁异常的提取与识别

本文对近年来国内外在烃渗漏与磁异常关系以及磁异常提取与识别技术领域取得的进展进行了论述与总结.利用磁异常寻找油气资源,特别是在圈划异常远景区方面具有经济、快速、有效的优势.但也存在多解性问题.为了有效判别异常源性质，本文提出利用人工神经网络作为以上方法补充的思路，可以作为磁法预测油气远景区技术流程中的一环，有助于减少多解性.     

辽宁省中强以上地震与地磁异常

辽宁省中强以上地震与地磁异常孙涛张桂林（辽宁省地震局沈阳１１００３１）王庆志王波于世昌史素娟（辽中地震站辽中１１０２００）１前言辽宁省自建立地磁台网及流动核旋测网以来，仅发生过几次４．８级以上地震，其中有前兆资料记载的只有３次，即１９７５年海城７．３...