基于小波细节信息的重力源参数估算

小波多尺度分析方法在位场数据处理领域得到了广泛关注.目前,人们对利用该方法分离位场数据所得到的小波细节横向边缘信息认知不够深入,有必要做进一步探究.本文首先采用小波多尺度位场分离方法分离出区域和局部重力异常,然后利用haar小波变换HVD模法对局部异常信息进行多尺度分析,将所得到的小波细节信息进行水平、垂直和对角线方向分解,探测出局部异常体的横向边缘特征,并利用频谱分析方法估测出场源深度.最后,通过理论模型和南非Witwatersrand Basin航空实测重力数据进行分析,表明小波多尺度分析方法不仅能很好地实现位场分离、近似提取异常体场源深度,还可以利用小波细节信息有效地分辨出局部异常体的横向边缘特征,提高其识别精度.     

离散小波变换与重力异常多重分解

小波变换多尺度分析是位场分解的有益工具.离散二维小波变换产生的低阶小波细节具有尺度不变的特征,它们不随小波变换的总阶数改变,而总阶数的增加仅增加高阶小波细节的个数和改变最后的高阶逼近.因此,用小波变换分解重力异常时,可根据地质目的来组合小波细节,从而实现有地质意义的异常分解.实验表明,用小波变换分解的重力异常小波细节,同样起到功率谱分解的作用,用分解之后的小波细节功率谱可以精确确定场源的埋深.中国大陆科学钻探场地的重力异常分解给出局部异常分解的典型例子.     

谱矩方法在磁源体深度反演中的应用研究

谱矩方法可以对数据的表面形貌做较为细致的描述.它以随机过程为理论基础,用各阶谱矩及统计不变量等具体的参数表征表面的几何形态,算术平均顶点曲率是一种基于四阶谱矩的统计不变量.通常,埋深不同的场源所引起的磁异常尺度不同,从曲率的角度来理解即为磁异常曲面的弯曲程度不同.因此,本文应用算术平均顶点曲率提取磁异常的几何信息,并将所提取的信息用于场源深度的反演.理论上推导了基于谱矩的球状磁源体和板状磁源体的反演公式,得到了场源深度与磁异常、曲率之间的关系式.结合理论模型计算验证了方法的有效性,并与欧拉反褶积方法进行对比.与传统的方法相比,该方法快速简单,无需调节参数,且有较好的反演精度.最后,将该方法用于塔里木盆地航磁异常的反演和解释中,反演出的磁源体深度可满足区域磁异常数据分析和解释的要求,为克拉通沉积盆地磁异常源的深度划分提供丰富的信息.     

塔里木盆地航磁场分析与磁源体结构

本文研究塔里木盆地区域磁异常图的反演及磁源体结构.由于众多异常的叠加和反演固有的多解性,区域磁异常图的准确解释是非常困难的.三维欧拉反褶积是一种确定地质体位置和埋藏深度的自动定量反演方法,比较适用于计算区域磁异常源的埋藏深度.由于大型克拉通沉积盆地地层具有上新下老的规律性,将磁异常源分解为三个深度层次,圈定它们各自的分布区域,便可将它们与形成的地质作用及时代联系起来,为准确解释区域磁异常图提供可靠的依据.本文应用三维欧拉反褶积反演方法,计算出的塔里木盆地深度为2~5 km、5~10 km、10~20 km三个等级的磁异常源,它们与形成的地质作用及时代分别为: 中生代构造运动,海西期玄武岩侵位和太古代结晶基底的变质作用;圈定了它们各自的分布区域.     

华北克拉通地壳三维密度结构与地质含义

利用重力异常揭示地壳三维密度结构是重力资料地质解释的目标和任务,密度反演的好坏至关重要.本文对华北克拉通进行了重力异常多尺度密度反演研究,首先利用小波变换对重力异常进行多尺度分解,接着利用功率谱分析方法估算各层场源的平均深度,然后利用广义密度反演方法进行各层密度反演,对华北地区进行了多层密度反演,得到其密度结构,并进行地质解释和油气藏分析.结果显示了重力场多尺度密度反演方法的有效性,对华北地区密度结构的地质含义进行了初步分析,位于造山带中的低密度异常主要反映沿造山带展布的花岗岩体.对比华北东部区浅层密度扰动与油气坳陷位置,发现油气坳陷都表现为低幅度的低密度异常区.说明利用小波多尺度反演提取的密度信息对油气勘探的部署有一定指导意义.     

小波变换方法在辽宁跨断层形变资料分析中的应用

运用二进小波变换方法,分析了辽宁省跨断层观测资料不同尺度的小波变换的细节信号变换特征。研究结果表明,不同尺度小波变换分解的细节信号显示出不同的时变特征,可以识别和提取原始测值中不易分辨的震前异常。当尺度为3时的细节信号的年周期特征被打破,或者当尺度为3时的细节信号变化超出2倍均方差时,测点周围250km左右范围内发生中强地震的可能性很大。     

华北时变重力场离散小波多尺度分解

利用华北地区地震重力监测网的绝对重力与相对重力多期重复观测资料,处理获取不同时空尺度的华北区域重力场动态变化图像.利用二维小波分解技术,将不同场源深度异常进行了分离,并对不同时空尺度的重力变化给出了解释,提高了对华北地区重力场变化趋势的认识水平.在此基础上,通过对小波细节的功率谱分析,可以获得各阶小波变换逼近及小波细节分析图所对应的场源深度.研究结果表明,不同时空尺度的重力场变化对于深入认识华北地区潜在地震危险性具有一定的科学意义.     

塔里木盆地多尺度重力场反演与密度结构

本文根据塔里木盆地布格重力异常的功率谱特征,决定重力异常分解成3个层次,即反映沉积盆地基底的异常、反映沉积盆地深层的异常和反映沉积盆地浅层的异常.然后,探讨将小波变换多尺度分析和密度反演方法应用于塔里木盆地布格重力异常的分解和盆地三维结构密度研究,并根据钻探了解每个层次对应着的地质含义.通过塔里木盆地的重力异常的多尺度...     

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

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

小波分析法在海拉尔地震台前兆数据处理中的应用研究

应用小波变换法可对不同频率成分进行分解的特性,选取d b4小波基函数,尺度为1～8阶,对海拉尔地震台前兆资料进行分析研究.结果表明,应用小波变换方法可较容易地识别和剔除气压、场地环境干扰等影响海拉尔地震台伸缩仪观测资料精度的诸多因素,气压和场地环境干扰在小波细节分解d3～d6部分异常显著;小波变换法对同震效应的识别、提...     

小波多尺度分解的振幅谱分维算法油气预测

本文结合小波变换理论和分形分维理论,利用小波变换的多尺度分解对地震资料进行分解,选取不同尺度的信息,计算其振幅谱分维数;根据含油气井的振幅谱确定能够反映油气异常信息的标度不变区,在该标度不变区内求分形维数,分维数值越大,说明含油气的可能性越大,据此预测储层的含油气性.实践表明,多尺度分形分维技术能够较准确地反映由地震波...     

重力场实数尺度小波分解及其应用

重力异常的多尺度分解是利用重力资料进行地质解释的重要环节,以往的小波多尺度分解,其小波细节尺度都是2的正整数次幂,这虽然非常有效,但在某些情况下,细节还是不够细致.因此,本文提出用无理数槡2的正整数次幂为尺度细节进行多尺度分解,以便取关于地壳三维密度结构的更多细节.中国陆地布格重力异常的实例表明,该方法得到的异常细节更为细致,新增加的无理尺度小波细节为我们提取了更多地壳构造信息,可对地质构造解释发挥更有力的作用.     

基于曲波变换的大地电磁二维稀疏正则化反演

为了提高二维大地电磁反演对异常体边界的刻画能力,我们引入曲波变换建立一种新的稀疏正则化反演方法.与传统的在空间域中对模型电阻率参数求解的方式不同,我们借助曲波变换将二维电阻率模型转换为曲波系数,并采用L1范数约束以保证系数的稀疏性.曲波变换是一种多尺度分析方法,其系数分为粗尺度系数和精细尺度系数,粗尺度的系数代表电阻率模型的整体概貌,而精细尺度中较大系数代表目标体的边缘细节.此外,曲波变换的窗函数满足各向异性尺度关系,并具有多方向性,因此曲波变换可以近似最佳地提取目标体的边缘特征信息,这为我们在反演中恢复边界提供有利条件.通过对大地电磁的理论模型合成数据和实测数据反演,验证了基于曲波变换稀疏正则化反演对异常体边界的刻画能力优于常规的L2范数和L1范数反演方法.     

分析研究基于小波分析与谱分析提高重力异常的分辨能力

重力数据是所有地下场源产生的重力场的叠加,探测对象经常被淹没在区域背景场之中,因此剩余异常的分离对于重力资料研究至关重要,而近来被引入位场领域的小波算子作为了滤波器和场源分析工具,在这里我们分析研究基于小波分析与谱分析的二维离散小波变换用于提高重力异常的分辨能力,再现出由简单形状场源描述密度不均匀的几何特征.本文先介绍二维多分辨率分析小波的基本理论及其提升算法,利用对数功率谱估计平均深度方法理论,接着对理论模型数据进行多尺度异常分解,估计地质体的形状、大小和深度,最后又对实测重力数据进行分析,并与传统常规方法进行比较分析,结果表明对于实际数据分析其方法也是具有可行性的.     

位场数据解释的Theta-Depth法

Theta图是利用位场(重磁)数据识别边界的常用方法,其表达式为重磁异常水平变化与垂直变化的比值函数.该方法计算浅源地质体边界的效果较好,而由于深源位场数据在换算过程中会产生趋同效应,在深源地质体识别应用中计算结果不准确,为此,本文提出Theta-Depth法并进行地质体埋深的计算.首先给出直接利用Theta图像进行场源体深度估算的方法,然后推导出基于Theta导数的线性方程来自动估算场源位置参数,本文方法可有效地利用Theta图像的特征为约束条件来提高反演结果的精度.理论模型试验证明本文提出的Theta-Depth法能有效地计算出场源体位置和深度.将该方法应用于满都拉地区实测磁数据的解释,帮助圈定了矿脉的分布.     

基于改进型Tilt-depth法的磁源上顶与下底深度快速反演（英文）

Tilt-depth法可以用于快速反演磁源的上顶埋深,但其无法反演磁源底部深度,同时该方法的反演解过于单一。针对于此,本文在有限厚度台阶磁异常导数公式基础上,推导出了可以同时反演地质体上、下界面埋深的改进型Tilt-depth法计算公式,并采用在Tilt梯度图上选取多特征点进行联合反演的计算模式来提高反演解的可靠性。二维及三维组合模型试验均证实了改进型Tilt-depth法可以有效地反演出磁性体的上顶与下底深度,且上顶埋深的反演精度明显高于常规方法的。最后将改进型Tiltdepth法用于松辽盆地长岭断陷航磁数据反演,其上顶埋深的反演结果与钻孔钻遇火山岩的深度基本一致,证实改进型Tilt-depth法相对于常规方法的反演结果更加准确。     

基于磁异常数据的居里面反演方法研究

本文总结了基于磁异常频谱特征反演居里面的基本理论和方法,分析了磁性体场源满足随机分布及分形分布两种假设下的磁异常频谱特征及其对数径向能谱与磁性体场源埋深之间的关系,并给出居里面埋深估计的方法.基于以上方法对华北地区磁异常数据进行分析并计算了相应的磁性体场源埋深,实际计算表明,基于磁性体场源随机分布假设的居里面深度计算结果与该地区地热结构有良好的对应关系.实际数据处理的结果验证了这些方法的可靠性和实用性.估算居里面深度的两种不同方法适用于不同的地质模型,实际应用应结合区域磁异常源分布的特点,以获取合理的居里面深度估计值.     

利用小波变换方法分析形变观测资料的正常背景变化特征

为定量刻画数字化形变观测资料中背景信息和噪声的时 频分布特征, 本文应用二进小波变换方法, 通过对小波分解的主模特征和随机白噪声识别因子变化特征分析, 剖析了山东数字化形变观测资料的正常动态背景和噪声变化规律. 结果表明, 当尺度取2, 3和4时, 分解后的细节部分存在着1/4日波、 半日波、 日波和半月波等准循环周期信号, 其中尤以尺度为3时的信号波幅最大; 尺度取1和5时的细节部分主要包含着噪声; 通过分析和追踪指定尺度的数字化形变观测资料小波变换的非震异常特征变化, 可望捕捉到与强地震孕育过程有关的前兆异常信息.      

陕西地磁短周期变化异常解释

1996年9月从陕西南部石泉至甘肃宁县布设地磁短周期变化观测临时测线．通过测点地磁转换函数,分析了该区地磁短周期变化特征．根据子波变换进行多尺度分析的理论,提取了该区地磁短周期变化异常的区域性源的影响成分．并以此为观测数据,反演研究了测线地下相应于区域性源的二维导电率结构．在反演研究中,将Zhdanov等提出的基于观测电磁场和导电率模型理论计算场的剩余场总能量最小,并且引入迁移剩余场概念的反演方法,引入磁变测深数据的反演中．反演结果给出了测线地下100 km深度范围内的导电率结构,指出高导电率区位于武功、乾陵一带地下;并且讨论了该反演方法的主要特点,即能较好地反映异常导电率区域的深度和导电率分布,且反演中能允许大数量的反演模型参数．      

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

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