关于利用△T′异常圈定磁性体边界及断层的一些问题

在电子计算机技术应用于物探资料处理以前,垂向二次导数的计算主要用于消除重力资料中的区域场以突出局部场和分解多体场源的迭加异常。随着频谱分析方法和电子计算机技术在物探资料解释中的逐步推广使用,垂向二次导数异常的计算成了磁异常数据处理和解释主要内容之一。垂向二次导数异常为磁异常的推断解释提供某些与地质体形态更为直接相关的信息,而大量的磁异常资料又为垂向二次导数方法的应用提供广阔的前景。从而,垂向二次导数方法又重新引起人们的兴趣和重视,近年来,国内就发表了不少关于垂向二次导数物     

重力垂直二次导数法的计算、意义和应用

引言众所周知,重力垂直二次导数是布伽异常沿垂直方向的二次导数,也就是重力位沿垂直方向的三次导数,是划分区域场和局部异常的有效方法.重力垂直二次导数有许多计算公式,所得结果不同,原因何在,究竟它们的意义是什么?它们之间有什么异同?能否提出一个准则来选择某一公式?计算精度如何评定?都是值得我们探讨的.本文结合一些实际成果粗略谈几点认识,抛砖引玉,请同志们批评指正.一二次导数法的一般意义大家知道,布伽异常中包含着比局部异常大得多的区域性异常,而和矿产分布有关的局部异常,往往被其掩盖或发生畸变.过去,我们习惯于从实测布伽异常图上来分辨局部异常,就是通常所讲的封闭和凸扭.这种     

辽宁省深部构造特征及其地质意义

本文应用了最新取得的“五统一”区域重力成果，并经过对浅表松散槽的密度亏损进行补偿改正后，求取了全省的区域重力场及其垂向二次导数，计算了莫氏面等深度图。在分析区域重力场及其垂向二次导数异常特征和莫氏面起伏特征的基础上，划分出辽宁深部构造的基本格局。     

位场各阶垂向导数的ISVD算法及其应用

位场垂向导数广泛应用于位场数据处理与解释当中,目前求取位场各阶垂向导数方法大致可以分为两类：①在频率域进行;②在空间域进行。针对这两类方法存在着稳定性较差的问题,这里介绍了计算位场各阶垂向导数的ISVD算法,并通过对单一地质体模型和地质体组合模型产生的布格重力异常进行处理分析,证明了ISVD算法在计算位场各阶垂向导数时,特别是在计算位场高阶垂向导数时,比在频率域及空间域计算结果具有较好的稳定性,ISVD算法能够有效地识别出浅层小尺度地质体边界。为了验证该方法对实际资料的处理效果,对某区块实际布格重力异常数据应用ISVD算法试算各阶垂向导数,结合地质背景等资料划分出断裂体系。结果证明,应用ISVD算法计算的位场各阶垂向导数,对于确定浅层小尺度地质体边界有较好的效果。     

重磁异常的导数

一、什么叫导数异常重(磁)场值沿某个方向两点之差与其距离之比称作场值的变化率,当两点的距离很小时这种变化率称为导数。由这些导数值组成的异常称作导数异常。根据求导方向的不同,主要分作垂直及水平导数两种;此外还有过渡型(如水平导数求垂向导数)。根据求导次数的不同,又可分作一阶、二阶及高阶(三阶以上)导数。     

青藏高原布格重力异常匹配滤波分析及其构造意义

利用匹配滤波的方法处理了青藏高原布格重力异常数据。基于重力异常一阶垂直导数的径向平均对数能谱曲线,青藏高原布格重力异常场可分为三层异常场。其中龙木错—双湖构造带在中部重力异常场和深部重力异常场图中都有显现;革吉—改则—错勤—申扎隐伏断裂在浅部和中部重力异常场中都存在。另一个明显的特征是藏北都显示有巨大且平缓的低重力异常圈闭,这与藏北Pn波速较低和Sn波缺失的特征吻合;在藏南存在一系列近东西向展布的重力异常断块。     

计算重力异常垂向二阶导数的DCT法——模型实验

基于离散余弦变换（DCT）的重力异常垂询二阶导数的计算方法是笔者提出的新方法。以无限长水平圆柱体为例研究了DCT法的计算规律：模型实验证实,该法计算的无限长水平圆柱体重力异常垂向二阶导数的精度与采样闻隔、剩余密度及圆柱体半径大小无关。与圆柱体轴心埋藏深度成正比。通过分析给出了基于DCT的合适的线性滤波方式,滤波后的重力异常垂向二阶导数与理论二阶导数拟合效果好,具有较高的计算精度。     

位场垂向高阶导数的Tikhonov正则化迭代法

在位场数据处理中,垂向导数具有重要的物理意义。其在一定程度上可以划分不同深度和大小异常源产生的叠加异常,且导数的阶次越高,这种分辨能力就越强,但通常认为高阶导数的换算是不稳定的。本文在Tikhonov正则化求位场垂向高阶导数的基础上,结合迭代法进行逐次逼近,提出了位场高阶导数的Tikhonov正则化迭代法,并且得到Tikhonov正则化迭代法的递推公式。通过对该方法的滤波特性分析可以看出,该方法计算的位场垂向高阶导数具有一定的稳定性及保幅性。模型试验和实际数据的处理表明,该方法计算结果较常规FFT求导法有更高的稳定性和实用价值。     

用等效源法进行重磁异常换算的效果

在磁力、重力异常的解释中,常常要进行异常场的换算,例如向上延拓、向下延拓,曲化平(将起伏地形上的观测值换算到某一水平面上的场值),磁异常不同分量之间的换算,求变换磁化方向后的异常(包括计算异常的化到磁极),以及异常的一次、二次垂向导数等等。这些换算已有一些方法可以进行,在上延、不同分量的换算和化到磁极等方面均取得了好的效果,但是在许多时候,特别对于观测面起伏高差较大的情况,常常效果不佳。而且有的方法丢掉边部测点的信息较多。在这里我们介绍重磁异常换算的等效源法,这种方法可以将上述不同方面场的换算问     

重力资料圈定地质构造边界方法的初步探讨

依据理论模型实验的方法,讨论了利用重力资料圈定地质构造边界方法的应用效果,分析了重力场垂向二阶导数、水平总梯度模极值、重力等值线梯度带、重力垂向一阶导数等方法在圈定地质构造边界方面的优缺点,为区域重力资料地质解释提供较为有效的依据.     

多环法计算重力垂向二次导数程序一例

多环法计算重力垂向二次导数程序一例王光杰（长春地质学院）多环平均值加权法是计算重力垂向二次导数常用的有效算法。在多年工作实践中，我们逐步优化了该方法的计算程序。该程序有三个特点：边界不损失，考虑了空区，短小易读。为了边缘不损失，计算平均值时，在边缘采...     

小波二维分解与NVDR-THDR技术在太原盆地重力位场中的应用

重力位场是地下多个不同密度的地质体所产生的场的叠加,太原盆地由于覆盖层较厚,基底断裂构造复杂,单一的位场信息提取技术不能达到划分太原盆地构造断裂、确定覆盖层厚度及基岩埋深的目标要求。笔者通过综合应用小波分解细节变换功率谱研究及利用归一化总水平导数垂向导数(NVDR-THDR)技术对太原盆地1∶5万布格重力异常信息的提取,认为小波分解细节变换技术对确定太原盆地覆盖层厚度及基岩埋深有较大的帮助,得出太原盆地中心覆盖层埋深大于3km的结论;认为NVDR-THDR技术对太原盆地重力位场中断裂构造的识别能够取得较好的效果,圈定断裂构造18条。     

位场垂向导数的稳定算法

垂向导数广泛应用于位场数据处理和解释中,常规的空间域求导公式多以台劳级数为基础,这不可避免会引入截断误差,且不便于求取高阶导数。另一种应用较多的波数域算法。其求导算子在高波数成份产生振荡失真,很不稳定,在实际应用中很难取得理想结果。这里提出了一种计算位场数据垂向导数的稳定算法,在波数域求垂向积分,在空间域计算二个二阶水平方向导数,借助拉普拉斯方程,求取位场异常的各阶垂向导数。算法只在求垂向积分时引入了一次付里叶变换,其它计算全部在空间域中进行,从而有效地避免了吉布斯效应的重复累加传递,有效降低了常规波数域求导算子在高波数成份产生的振荡失真。通过理论模型对比,证明了这里所提出的算法较常规的波数域方法具有良好的稳定性,特别是在计算高阶垂向导数时有明显优势,能有效地提高与垂向导数有关的位场定量解释方法的准确性。     

用有限元插值法计算区域重力异常

有限元插值法利用区域边界上的少量数据可计算区域重力异常。由于只有区域数据参与插值计算,其结果更准确。通过教学模型实验得出结论:在区域场近似二次或三次曲面的条件下,有限元插值法比最小二乘法计算的区域场更接近真实值。在某盆地的重力资料处理中,用有限元插值法分离出该区的局部重力异常.     

直接计算磁性下界面深度的功率谱法

以计算居里等温面深度为目的,本文给出了利用磁异常的垂向导数的功率谱来直接计算磁性下界面深度的方法。同时,借助于异常场的谱与其垂向导数谱之间的关系,还可由用磁异常垂向导数功率谱计算的公式转换出用异常场的功率谱来计算磁性体下底深度的公式。理论模型实验表明,两种计算方法均能得到较好的结果。最后给出了用该方法对泛华北地区航磁异常的试算结果。     

利用重磁水平和垂直二阶导数确定东北地区梯度带

分析了东北地区的重、磁场特征,同时对研究区的布格重力异常和航磁异常数据进行水平和垂直求导。水平和垂直参量图能给出断裂和走向的信息。结合布格重力异常图可以得出东北地区重力梯度带的走向和边界,并对比重力和磁力的垂直和各个角度的水平二阶导数来综合确定东北地区梯度带。利用求重力和磁力的水平和垂直二阶导数来划分东北地区梯度带。     

关于位场垂直二次导数换算公式的频率响应

在重、磁资料数据处理中,经常进行垂直二次导数的换算,以消除区域异常干扰,突出局部异常。位场垂直二次导数换算公式的种类繁多,且使用效果不甚一样。为了评价各公式的计算精度,正确地选用公式及其环半径,已有一些地球物理工作者从各公式的频率响应特性出发讨论了这些问题[1－3],     

利用重、磁资料研究于都—赣县矿集区盘古山地区断裂构造及花岗岩体分布

盘古山地区是于都—赣县矿集区内一个重要的钨锡矿成矿区,研究区断裂构造和花岗岩体与成矿关系密切。利用重、磁资料研究盘古山地区断裂构造及花岗岩体分布,为深部找矿提供地球物理依据。利用重力异常归一化总水平导数垂向导数技术推断断裂构造平面位置;利用欧拉反褶积方法对重力异常垂向一阶导数进行反演计算,推断断裂构造平均深度;利用重磁异常垂向一阶导数技术推断花岗岩体平面位置;利用RGIS软件2.5D重磁剖面人机交互正反演技术推断花岗岩体断面位置,并用钻孔验证解释结果的正确性;利用RGIS软件3D重磁模型编辑模块展示花岗岩体空间位置。研究结果显示盘古山地区既有出露断裂,又有隐伏断裂,既有出露低密度、弱磁性花岗岩体,又有隐伏低密度、强磁性花岗岩体,花岗岩体分布受断裂构造控制,矿体分布受断裂构造和花岗岩体共同控制。     

航空重磁在山东齐河厚覆盖区圈定隐伏岩体的方法及效果

山东齐河地区被大面积巨厚的新生代地层覆盖,但富铁矿资源潜力大,隐伏岩体的圈定对于提高厚覆盖区基础地质研究程度和矿产勘查具有重要意义.通过航空重磁数据处理提高解释的准确性,利用航磁△T化极垂向一阶导数、归一化总水平导数垂向导数、航空重力剩余异常等实现隐伏岩体范围圈定,在研究区共圈定6处中基性岩体,并利用航磁切线法、航空重...     

改进的倾斜角法在位场场源边界识别中的应用

位场场源边界识别是位场数据分析解释中的重点内容之一。相邻异常源由于场的叠加而难以区分,尤其是对于叠加异常中不同深度的场源边界的识别。倾斜角法是位场数据处理中常用的边界增强方法之一,能有效突出低振幅、短波长的异常特征。本文基于倾斜角法定义了新的位场边界探测方法,该方法利用倾斜角方程中位场二阶垂直导数与垂向导数水平导数之比进行计算,能有效区分相邻异常源产生的叠加异常,突出叠加异常中不同深度场源的边界特征。通过综合模型和实测数据,并与现有的基于导数的位场边界探测方法进行比较,评价了该方法的应用效果。结果表明,该方法对场源深度的变化不太敏感,它可以更准确地确定场源边界,区分相邻异常源的叠加异常,尤其对密集型组合地质构造的识别有较好的效果,即使是在场源埋深较大情况下,也能突出叠加异常中不太明显的微弱变化的线性异常特征,以及弱异常源的形态特征,可为精细研究局部地质构造提供技术支持。