重力资料多尺度分析最优小波基的选择

小波变换多尺度分析是重力资料处理中一种常用的方法,其分析的结果与最优小波基的选取有关。笔者结合重力数据特征,首先从理论层面上分析最优小波基的选取方法,再根据理论模型进行不同小波基的对比实验,最后对华北地区实测重力数据进行小波基对比试验,认为db11小波基为重力数据多尺度分析的最优小波基。     

应用奇异值分解与二维经验模型分解提取金矿化致矿重力异常

奇异值分解(singular value decomposition,SVD)和二维经验模型分解(bi-dimensional empirical mode decomposition,BEMD)2种方法用于鲁西铜石金矿田的金矿化致矿重力异常提取.通过对比研究得出如下结论:①基于重力数据,由BEMD和SVD 2种方法获...     

基于经验模态分解的重力异常分离

根据重力异常水平梯度特性,建立了基于经验模态分解的新型重力异常分离方法。通过不同模型实验对比,发现该方法分离所得异常和模型理论异常能很好地吻合,可进行异常定量划分。实测数据处理中,将该方法与平均对数功率谱分析相结合,所分离出的区域异常可以有效描述研究区的区域地质特征。模型试验和实测数据处理都充分表明该方法所分离出的重力异常基本没有畸变和虚假信息产生,具有良好的实际应用价值。     

基于二维经验模态分解的重磁异常分离

由于地质过程的复杂性及成矿过程的多期次叠加性,原始重磁异常往往是多种地质因素的混合信息,既包含区域背景异常信息,也包含与矿床(体)、矿化蚀变带以及隐伏岩体等与找矿密切相关的地质要素所引起的局部重磁异常.如何从复杂的叠加重磁异常中分离出具有找矿意义的局部异常,是当前矿产勘查和资源潜力评价工作中面临的难题之一.采用经验模态分解(EMD)方法来分解重磁异常,为提高分解的稳健性提出了用双调和样条插值(BSI)进行包罗面插值的新方法,并以云南个旧地区重磁数据为例,对其进行非线性多尺度分解,实现对区域异常与局部重磁异常的分离,揭示了深层次找矿信息并拓宽了经验模态分解方法的应用领域.      

基于二维小波分解的MT二维反演方法

对于大地电磁反演,期望得到快速稳定的反演结果并且能够反映清晰的、未发生畸变的地质界面.为了解决这一问题,提出了基于二维小波分解的MT二维反演技术,对小波分解得到的背景场采用光滑反演,细节部分采用聚焦反演.通过模型试验,证明了该方法的有效性,同时对某区的实际资料进行了处理,表明了该方法的适用性和效果.     

重力位场小波多尺度分解性质的分析与应用

小波多尺度分解是重力位场分离的常用方法之一,其最大的优点是突破了传统的二分重力异常的理念,实现了重力异常的多重分解,但对其性质及其存在的局限性还未进行系统研究。为了更好地指导实践,本文从基于剖面和格网的重力位场小波多尺度分解的定义出发,阐述了低阶小波细节不变准则等3个重要性质,通过设计简单和复杂两类典型理论模型,分析了小波多尺度定义及性质,并在应用实例中与插值切割法进行了对比。结果表明,小波多尺度分解可以实现重力位场多层分离并推估源体埋深,同时针对本文指出的异常尺度混叠和比例系数难以确定等局限,给出了改进思路。以上基础研究工作,可为重力位场资料的处理和解释提供一定参考。     

基于二维经验模态分解的地磁基准图制备边界补偿方法

针对采用矩谐分析方法构建地磁基准图的边界震荡问题, 提出一种基于二维经验模态分解的区域地磁异常数据边界补偿方法.采用二维经验模态分解方法对区域地磁异常数据进行多尺度分解, 对分解所得小尺度本征模态函数分量, 利用总体Hilbert变换法进行瞬时频率和瞬时幅值特征提取, 通过自采样和特征匹配进行边界补偿; 将大尺度分量之和作为趋势项, 利用三角函数方法建立模型并计算边界之外的大尺度磁异常值.实验证明, 相比当前已有方法, 可以更加有效地抑制对区域地磁异常数据进行矩谐分析的边界震荡问题, 稳定提高构建地磁基准图的准确性.      

基于经验模态分解的航空重力测量动态误差分离

航空重力测量数据在测量的过程中很容易受到飞机动态运动的影响,飞机动态运动强时,航空重力测量数据的质量也随之降低,虽然经过数字滤波后可以将航空重力异常数据中的大部分噪声去除,但一部分与信号重叠的噪声还是无法被消除.利用经验模态分解对受飞机动态影响较大的一组实测数据进行深化处理.结果表明,在不影响航空重力异常数据分辨率的情况下,该方法很好地消除了由飞机动态运动引起的动态误差,使得航空重力测量的内符合精度由1.43mGal提升到了1.27mGal.     

基于随机介质模型的二维重力正演

岩(矿)石密度的空间分布具有随机性和相关性,采用随机过程的谱分解理论和指数型椭圆自相关函数理论,结合江西德兴某重力勘探剖面资料,建立了随机介质模型,并进行重力异常计算。算例表明,采用随机建模的方法可以很好地模拟实际勘探断面的密度结构,据此计算得到的重力异常特征也与实测的基本一致。因此,在重力资料的解释中,根据实际的地质和密度测试资料建立随机介质模型进行正演计算,有助于对实测异常的识别、分析和解释,提高解释的可靠性。     

基于小波多尺度分析的重力异常分离

在利用小波多尺度分解重磁场中,实际上并非所有的小波都适合于位场分离.这里采用对称性、紧支撑性较好的Daubechies小波对重力场多尺度分解,依照重力背景场的地质意义提取重力异常背景场,并在此基础上计算重力异常异常场.最后通过模型试算和实例分析进一步表明了该方法的有效性和实用性.     

松辽盆地布格重力异常小波多尺度分解与解释

松辽盆地是中国东北盆地群的代表性盆地之一。松辽盆地及其邻区的布格重力异常分布特征以北东—北北东向为 主,其次为东西向,重力场分区分块特征明显,显示了明显的拼贴构造背景,依据布格重力异常特征对研究区构造单元进 行分区,共分成松辽盆地异常区、大兴安岭重力梯级带等 5 个特征区。利用小波多尺度分解方法对松辽盆地重力场进行分 离,选择最佳的分离阶数(4 阶),随后利用对数功率谱确定每阶异常的视深度,提取区域异常,分层研究地壳顶部(D1G)、基 底(D2G,D3G)及深部密度界面(A4G)的起伏特征,共圈定基底凹陷 14 个,并绘制了研究区的 Moho 面埋深图。地球物理特征 揭示了松辽盆地的深部构造,为研究盆地油气勘探及构造演化提供了重要依据。     

多尺度地震资料联合成像

为改善储层目标区域的成像质量,结合地面地震、VSP和井间地震的各自优势,采用归一化后的自适应叠加方法,实现3种资料偏移结果的联合成像,以达到各种资料优势互补、提高成像精度的目的。为了给联合成像提供精度较高的单种资料成像剖面,采用傅里叶有限差分叠前深度偏移算法实现地面地震资料成像;在VSP和井间地震资料的叠前成像中采用了基于起伏地表直接下延法成像思想,基于傅里叶有限差分延拓算子得到VSP和井间地震偏移结果。模型试算和实际资料的处理验证了多尺度地震资料联合成像法的正确性和有效性。     

基于DCT有限元二维重力张量的正演

介绍了国内、外张量测量技术及数据处理技术的发展情况。把离散余弦变换与有限单元方法结合起来,进行重力张量的正演计算。首先,利用有限单元法求解重力场的一阶导数,然后利用离散余弦变换求解重力张量分量。这样既保留了有限单元法对复杂模型的正演优势,又避免了利用有限元法直接求解重力张量时,因多次插值引起的数据量过大的问题。模型试验证明,用基于DCT的有限元法进行重力张量的正演,不但具有很高的计算精度,而且相对于利用有限单元法直接正演,其计算速度得到了显著的提高。     

区域重力异常多尺度密度反演研究

重力测量观测得到的重力异常,反映的是不同深度和不同尺度异常源的叠加效应。在理论上,只有进行适当的异常分解,才能较为准确地反演异常体的物性与几何参数。这里首先对重力异常进行小波多尺度分解,再通过频谱分析求出各层对应的中心深度,然后确定各层的平均密度,用广义线性反演方法反演各层密度的扰动,进而得到各层的密度分布。通过塔里木盆地的重力异常的多尺度分析和密度反演应用实例,显示了多尺度反演方法对区域重力异常的有效性和应用价值。     

基于小波多尺度分析的重力异常分离CSCD

在利用小波多尺度分解重磁场中,实际上并非所有的小波都适合于位场分离。这里采用对称性、紧支撑性较好的Daubechies小波对重力场多尺度分解,依照重力背景场的地质意义提取重力异常背景场,并在此基础上计算重力异常异常场。最后通过模型试算和实例分析进一步表明了该方法的有效性和实用性。     

二维经验模分解（BEMD）技术在提取薄竹山矿集区致矿重力异常中的应用

正滇东南地区位于华南地块右江海西-印支地槽褶皱带,区内蕴含着我国重要的钨、锡、银、铅、锌等多金属矿产,是华南岩浆成矿带的重要组成部分。这些矿产均围绕个旧、薄竹山和老君山3个复式花岗岩体分布,其中已有大量资料证明个旧花岗岩体(秦德先等,2008;毛景文等,2008;程彦博等,2008;Cheng Y B et al., 2013)和老君山花岗岩体(Yan D P et al., 2006;张斌辉等,2012;蓝江波     

南海北部缘重力异常的多尺度分析及其构造讨论

重力异常是由地下各种密度不均匀体所产生的叠加效果。采用二维小波多尺度方法对南海北部缘重力异常进行分析,解译了重力异常所反映的地下源体,对南海北部缘盆地构造进行认识。通过对南海北部缘重力异常进行分解,结合研究区域地形地貌,对盆地构造和地下岩体分布等地质特征进行分析。结果显示,四阶变换逼近主要由莫霍面起伏引起,小波变换细节主要由沉积基底面的起伏引起。利用小波多尺度对重力异常进行分析具有较好的效果,值得进一步推广研究。     

冈底斯带重磁异常二维经验模态分解及地壳结构

地壳结构是青藏高原深部地球物理探测的主要目标之一.研究了结合功率谱分析的二维经验模态分解（Bi-dimensional empirical mode decomposition,BEMD）位场分离方法,并将其应用于冈底斯带重磁异常分离中.模型试验结果表明该方法能够分离叠加异常,并可获得场源平均深度的可靠估计值;实际资料处理表明研究区布格重力异常和化磁极异常均可分为由浅到深的3个等效层场源产生的异常组合.其中深部重、磁场以南北分区、东西分块为特征.证明了青藏高原是南北地体的拼接,并且在形成之后持续受到印度板块不均匀地向北推进,使喜马拉雅地体、冈底斯地体的深部物质属性在东西方向上以约88°E和93°~95°E为界存在显著的差异;推断青藏高原南部区域中下地壳的低密度体东西向是不连续的,且在研究区中、东部,南北向低密度体被雅鲁藏布江缝合带切开.      

二维倾斜断层重力异常的正反演

本文将侯重初同志在文献中提出的上延叠加场思想,有效地用于反演不同密度岩层的倾斜接触面,提高了确定倾斜台阶各参数的精度。并将上延叠加场思想的逆过程,使用于空间域和频率域中倾斜台阶重力异常二阶导数的正演计算。本文另一结果是提供用富里叶交换确定倾斜台阶下底面的深度。     

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

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