三度体梯度磁异常全方位正反演方法

首创了磁法勘探学科中三度体任意方向梯度磁异常全方位正、反演方法。首先 ,在球坐标系中 ,导出了三度体任意方向梯度磁异常的两个球谐级数通用正演表达式 ,为实现全方位反演奠定了理论基础。然后 ,由第 1个梯度磁异常球谐级数表达式出发 ,提出了三度体总磁化方向、磁矩模值反演方法 ;由第 2个梯度磁异常球谐级数表达式出发 ,提出了三度体磁矩重心、近似形态与分布范围 ,以及有广泛实用价值的“均磁凸形体”边界点位置等全方位反演方法。最后 ,给出了 1个三度体梯度磁异常全方位反演的理论模型验证实例。     

三度体梯度磁异常全方位正反淀方法

首创了磁法勘探学科中三度体任意方向梯度磁异常全方位正、反演方法。首先,在球坐标系中,导出了三度体任意方向梯度磁异常的两个球谐级数通用正演表达式,为实现全方位反演奠定了理论基础。然后,由第１个梯度磁异常球谐级数表达式出发,提出了三度体总磁化方向、磁矩模值反演方法;由第２个梯度磁异常球谐级数表达式出发,提出了三度体磁矩重心、近似形态与分布范围,以及有广泛实用价值的“均磁凸形体”边界点位置等全方位反演方     

局部重磁场源全方位成像理论概要(续)

作者近年来首创了"复杂条件下局部重磁场源全方位成像"理论体系。在本文中,概要阐述该理论体系中下延有限二度体复重磁场级数正演通式、复场模值全方位延拓和全方位反演等。     

起伏地形上规则二度体复重磁场正演和直接反演

"起伏地形上规则二度体复重磁场正演和直接反演"由4篇连续文章组成.第一篇和第二篇阐述复重磁场正演问题,主要采用前苏联学者的方法,在复数域内,严格而简明地推导出了规则二度体及其组合体复重磁场表达式.这些复场表达式与直角坐标系内对应形体重磁场表达式相比,具有非常简单的形式.第三篇和第四篇阐述复重磁场直接反演问题,是本系列文章的核心,从规则二度体及其组合体复重磁场表达式出发,完全在复坐标系内,采用独创的线性化方法,建立起求解中间变量的线性方程组,再结合复高阶方程解法,成功地解决了起伏地形上(或者观测线处于二维场源任意方位甚至封闭时)几乎所有规则二度体重磁场的直接反演问题,推动了重磁场解释理论的发展.     

起伏地形上规则二度体复重磁场正演和直接反演

"起伏地形上规则二度体复重磁场正演和直接反演"由4篇连续文章组成。第一篇和第二篇阐述复重磁场正演问题,主要采用前苏联学者的方法,在复数域内,严格而简明地推导出了规则二度体及其组合体复重磁场表达式。这些复场表达式与直角坐标系内对应形体重磁场表达式相比,具有非常简单的形式。第三篇和第四篇阐述复重磁场直接反演问题,是本系列文章的核心,从规则二度体及其组合体复重磁场表达式出发,完全在复坐标系内,采用独创的线性化方法,建立起求解中间变量的线性方程组,再结合复高阶方程解法,成功地解决了起伏地形上(或者观测线处于二维场源任意方位甚至封闭时)几乎所有规则二度体重磁场的直接反演问题,推动了重磁场解释理论的发展。     

三度体外磁场三分量面元量板的计算方法

三度体外磁场三分量面元量板(简称三度体量板)与桂林冶金地质研究所似二度体量板的数理推导类似,都以小球的磁场解析式为基础;其差别是先制成小扇形块的Ⅰ、Ⅱ、Ⅲ三个面元量板,算出三度体各截面的面元磁场,然后再沿走向进行近似积分.三度体量板可以制成与似二度体量板完全相同的形式.为了解决计算点到不同距离的面元都能通用一套量板的问题,对量板相邻半径的环带间隔给予了不同环带系数.用面元量板计算出三度体的各截面(一般只需要计算包括极大、极小和转折面等3～5个截面即可)面元磁场值后,可用任何近似积分方法计算出磁场的三个分量.所介绍的这种方法虽仍较繁,但目前在物探野外工作中还没有简单实用的计算三度体磁场的手算方法,因此还是可用的.经过试算,对较简单的三度体,两个熟练的计算人员每日可以计算5～7个点的磁场三分量,相当于对一条野外实测剖面概略地算出     

三度体球冠模型重力异常计算方法及其重力归一化总梯度场特征

提出利用铅垂柱体薄片作为面元的面元积分法,计算三度体球冠重力异常的计算公式,用三度体球冠模型模拟三度背斜体,计算了均匀密度球冠模型、非均匀密度储油球冠模型的重力异常。最后通过误差分析,验证了该方法的可行性。采用波数域重力归一化总梯度计算方法,计算了均匀密度球冠模型、非均匀密度储油气球冠模型G^H场,发现似三度背斜体与二度背斜体有相似的G^H场等值线特征,似三度贮油气藏背斜构造的G^H场亦表现出明显的“两高夹一低”的储油典型标志。所以,可以用非均匀密度三度体储油球冠模型模拟三度贮油气藏背斜构造。     

计算非均匀磁化条件下磁体内外磁场的简化方法

非均匀磁化、任意形状的磁性体的磁场计算是经典场论还一直未解决的问题.当磁性体磁性较强、形态复杂、埋藏较浅时,特别是穿过强磁性体的井中磁测结果,不均匀磁化对磁异常影响很大,不可以忽略. 最近在我国先后发表过计算任意形状强磁性三度体磁场的论文其中黄诚的泊松公式面积元素法,是根据离散化概念提出的,可用来计算三度体全空间的磁场.用近似程度不同的两套公式     

重磁异常的人机联作校正－迭代反演

本文提出了三度体重磁异常的人机联作校正－迭代反演方法。该方法用二度半组合多边形棱柱体来逼近三度体，从而把三度体重磁异常反演问题转化为二度半体的反演问题；为了消除组合体迭加场的影响，该方法采用了一种校正－迭代技术。理论模型反演计算表明，该方法实际可行。     

重磁异常的人机联作校正－迭代反演

本文提出了三度体重磁异常的人机联作校正－迭代反演方法。该方法用二度半组合多边形棱柱体来逼近三度体，从而把三度体重磁异常反演问题转化为二度半体的反演问题；为了消除组合体迭加场的影响，该方法采用了一种校正－迭代技术。理论模型反演计算表明，该方法实际可行。     

重磁方法在金川铜镍矿东延M-15异常勘查中的应用

M-15磁异常(东湾异常)是国内外学者关注的寻找铜镍矿的重点预测区之一。该异常处于全覆盖区,前人推断它是由超基性岩体引起的。本次采用重、磁方法和已知白家咀子超基性岩体磁异常特征对比,认为M-15系非同源重磁异常,以磁力高、重力低为组合特征,异常体埋深较大,磁异常可能为1200 m以下的超基性岩引起,重力低异常系埋深1 000 m左右的龙首山群元古界基地凹陷引起。这一结论对于扩大铜镍矿源起到了积极的作用。     

基于GPU的任意三维复杂形体重磁异常快速计算

提出了基于图形处理单元的任意三维复杂形体的重磁异常快速正演计算方法。将地下半空间剖分为大小相等规则排列的一组长方体单元,任意三维复杂形体可以表示成很多不同体积和密度(磁性)的长方体的近似组合。用解析方法计算出所有这些长方体在计算点的重力(磁力)异常,并累加求和,就可以得到整个模型体在计算点引起的重(磁)异常值。为了提高近似程度,需将地下半空间剖分得很细,用传统的CPU串行程序计算相当耗时。GPU在处理能力和存储器带宽上相对CPU有明显优势,采用GPU并行算法,可大大提高计算速度。相关试验结果表明,用GPU实现的正演快速算法计算结果正确,效率明显提高,为重磁异常三维物性反演提供了基础。     

均质多面体重磁异常正演计算表达式的一致性

根据重力位、重力异常、磁位和磁异常公式,应用泊松关系,通过严格的公式推导,给出了重力位、重力异常、磁位和磁异常之间的数学关系,从而使重磁异常正演计算表达式包含于简单的规则之中,揭示出均质多面体重磁异常正演理论的一致性。     

利用平台磁异常的综合地球物理特征推断矿体的成因类型

简要介绍了接触交代—热液型铁矿床和沉积变质铁矿床形成的地质环境、控制条件,矿体形态特点和经济利用价值,通过平台磁异常Za-Ha双分量参量图研究,拟合反演,提出平台磁异常可能对应陡立接触带形态的磁性体。在重磁相关分析中,由相关系数的大小判断出平台磁异常重磁不同源,说明平台磁异常主体部分对应的磁性体不是沉积变质铁矿床。在CSAMT测深法探测结果分析中,提出平台磁异常深部可能存在岩浆岩侵入活动的物探证据。从而说明平台磁异常区具备形成接触交代—热液型铁矿床的可能性和几个必要条件,综合地质、物探资料推测平台磁异常成矿类型可能是接触交代—热液型铁矿床。     

大兴安岭中段地球物理特征及其地质解释

以区域地质资料和物性资料为依据 ,对大兴安岭中段地区重、磁异常特征及其对应的地质情况进行综合解释 ;探讨区域矿产分布与重磁场的关系 ,以及不同属性矿产的磁异常特点 ,指出该区重磁异常寻找多金属矿床有指导作用     

The Glueckstadt Graben of the North-German Basin: new insights into the structure from 3D and 2D gravity analyses

The structure of the Glueckstadt Graben has been investigated by use of 3D gravity backstripping technique and by 2D gravity and magnetic modelling. Subtracting the gravity effects of the Meso-Cenozoic sediments together with Permian salt reveals a positive residual anomaly within the Glueckstadt Graben. This anomaly includes two local maxima over the Westholstein and Eastholstein Troughs. The 2D gravity models point to the presence of a high-density body within the lower crust of the Glueckstadt Graben. In addition, the results of 2D magnetic modelling indicate that the central part of the high-density body is overlain by an area with high susceptibility. Most probable, the formation of this high-density body is a result of complex poly-phase tectonic history of the study area. Finally, the results of gravity modelling indicate that Permian salt is not homogeneous. 3D gravity analysis and, especially, 2D gravity modelling have distinguished the differences in degree of salt saturation in salt-rich bodies, and elucidate the proportion of Rotliegend salt.     

重磁2.5D/3D互相约束反演技术在三维地质填图中的应用研究

三维地质填图为我国启动的新一轮地质调查项目,为提高重磁资料在三维地质填图中的应用效果,笔者提出重磁资料2.5D/3D相互约束重磁反演技术方案:利用重磁资料2.5D剖面反演结果、3D物性反演结果作为彼此反演约束条件,并通过了理论模型试验。试验结果表明,该技术方案使反演结果中物性参数、空间位置更接近理论模型体。通过对本溪—临江地区思山岭铁矿磁异常及酸性岩体重力异常进行反演实践——估算磁异常铁矿资源量、研究酸性侵入岩深部展布形态的效果良好,可为大面积三维填图提供有效途径。     

海南石碌地区重磁异常特征及深部找矿

通过对石碌铁矿区重磁异常特征的分析,总结了该矿区重磁异常特征与矿体之关系.探索重磁方法在该矿区找矿的有效性,建立其找矿模型,指导该矿区深部及外围的地质找矿工作.从而扩大矿区远景规模,增加矿区资源储量,为目前正在开展的全国矿产资源潜力评价项目提供资料,具有一定的实际意义.     

在波数域计算一维重磁异常导数的Matlab语言算法

利用Matlab内建的快速傅氏变换函数可以方便地在波数域计算重磁异常导数。介绍了基于Matlab语言的波数域求导的算法,给出了程序源代码,讨论了一些有助于提高计算精度的编程技巧。通过模型试验和数据分析,发现在计算垂向导数时波数域求导算法的精度比傅氏级数的精度有明显改善;而水平导数的计算,2种方法的精度相当。在某区钾盐勘探中,用该方法处理高精度重力剖面数据,取得了较好效果。