基于欧拉反褶积方法计算辽宁地区重力变化场源特征

以2012—2018年辽宁地区的流动重力观测资料为基础,结合实际地下地质体构造特征构建理论模型,应用欧拉反褶积法对重力变化场源深度及空间分布特征进行反演和解释。结果表明:在辽宁西部辽蒙交界和海城—岫岩断裂带附近场源位置集中,地下介质活动增强;欧拉反褶积方法适用于对流动重力观测资料的重力场反演及以场至源的定量研究应用。     

鄂尔多斯南缘地区重力变化场源特征

以鄂尔多斯南缘地区布格重力异常数据及2014—2017年相对重力观测数据为基础,采用欧拉反褶积方法,对引起重力变化的场源深度进行反演,并对空间分布规律予以解释。通过构建理论模型,反演得到最优模型参数,并对实际数据进行计算和分析。为了减弱和消除欧拉解的发散性,利用水平梯度滤波法优化反演结果。结果表明:①构造指数为1时,适合对相对重力数据进行反演,当滑动窗口选择4—8倍测点间距时,可以获得较为可靠的场源参数;②在主要断裂附近,不同时间尺度的场源呈较好的一致性特征;③引起重力变化的场源深度集中在10—30 km,与鄂尔多斯南缘地区的震源深度基本一致;④不同时间尺度的场源位置相对分散,表明引起重力变化的物质流动具有随机性和波动性。     

基于垂向一阶导数与解析信号比值的欧拉反演方法(英文)

常规欧拉反褶积法中构造指数的选取以及分散解存在较多的问题,提出了基于联立垂向一阶导数与解析信号的欧拉齐次方程的RDAS-Euler反演方法。该方法可以更为精确的估计场源的范围及埋深,且不需考虑构造指数N的影响,避免了因构造指数不当而引起的反演误差。通过对单一地质体及组合地质体模型的实验证明本文方法能有效地完成目标体的反演工作,反演结果与理论值之间的误差小于10%,且相对于常规欧拉反褶积法更加稳定准确,能够更好的得到地质体边界及深度信息。将RDAS-Euler法应用于黑龙江省虎林盆地实测布格重力异常数据,获得了丰富的断裂信息,说明RDASEuler法增强了对断裂平面位置的识别能力。     

基于垂向一阶导数与解析信号比值的欧拉反演方法

常规欧拉反褶积法中构造指数的选取以及分散解存在较多的问题,提出了基于联立垂向一阶导数与解析信号的欧拉齐次方程的RDAS-Euler反演方法。该方法可以更为精确的估计场源的范围及埋深,且不需考虑构造指数N的影响,避免了因构造指数不当而引起的反演误差。通过对单一地质体及组合地质体模型的实验证明本文方法能有效地完成目标体的反演工作,反演结果与理论值之间的误差小于10%,且相对于常规欧拉反褶积法更加稳定准确,能够更好的得到地质体边界及深度信息。将RDAS-Euler法应用于黑龙江省虎林盆地实测布格重力异常数据,获得了丰富的断裂信息,说明RDASEuler法增强了对断裂平面位置的识别能力。     

基于欧拉反褶积方法计算川滇交界重力变化场源特征

本文以川滇地区2012-2014年来的多期流动重力观测网络数据为基础,采用欧拉反褶积方法对引起重力变化的场源深度和空间分布规律进行了反演和解释.本文通过构建理论模型。得到了最优的反演参数,并对实际数据进行了计算和分析.结果表明:选择构造指数(SI)为1时适合对流动重力数据进行反演.实际资料的反演结果显示在2012-2014年间的重力变化场源位置更集中于垂直于昭通断裂的北西向分布,场源深度分布在30±10 km范围,形态与鲁甸地震发震构造位置相似.这可能与2014年8月3日鲁甸M_S 6.5地震前的中上地壳的深部应力场变化相关.依据本文研究结果和方法,可以开展由场至源的定量研究,为流动重力变化信号的反演和解释提供一条新思路.     

基于欧拉反褶积方法估算华北地区重力变化的等效源参数

利用华北地区2009—2014年绝对重力与相对重力多期重复观测资料, 得到不同时空尺度的华北区域重力场动态变化图像。 采用欧拉反褶积方法, 通过对理论模型试算, 获得最优反演参数, 对引起华北地区重力场变化的场源深度和空间分布规律进行了反演计算和分析。 结果表明: 当构造指数为1时, 适合对流动重力场变化数据进行反演。 实际资料的反演结果在2009—2014年间的场源位置集中于河套断裂带。 本文研究方法可以用于重力场反演和以场求源的定量研究。 本工作为流动重力变化信号的反演和解释提供一条新思路, 也为构建地震重力预报指标体系提供了定量依据。     

辽宁地区时变重力场源变化特征分析

地震重力分析通过研究时变重力场变化获取地球内部介质物性变化信息.采用贝叶斯重力网平差方法对辽宁地区2011-2014年共计四年7期的流动重力观测资料进行处理,对研究区重力观测网总体监测能力做出分析,选用欧拉反褶积对研究区的重力变化场源深度及空间分布特征进行反演和解释.通过反演计算发现,在2013年灯塔Ms5.1地震前沈...     

不同高度数据联合欧拉反褶积法研究

基于不同测量高度重力场及其梯度数据可同时对应同一场源并用于反演场源位置的分析原理,拓展不同高度场数据在欧拉反褶积法中的应用范围.首先,立足于对欧拉反褶积方法的理论研究基础,提出不同高度数据融合联合欧拉反演公式.其次,在理论模型上对多种高度数据联合反演做了测试分析计算,验证了不同高度场数据融合联合欧拉反褶积法能够改善位场解释中单一观测面数据计算带来的解的发散问题,收敛过程由此改善.最后,将本文方法应用于龙门山地区实际重力数据的解释,获得了研究区断裂分布特征.     

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

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

基于二阶导数的磁源边界与顶部深度快速反演

为实现磁异常反演,首先提出了磁异常非参数快速反演的概念,即无需提供先验信息.在Nabighian提出的磁场通用梯度公式的基础上,推导实现了基于磁异常垂向二阶导数的非参数快速反演(V2D_depth).它不仅可以获取场源的边界信息,同时可以反演场源的埋深.通过与Tilt_depth方法对比,本文方法计算的场源边界更清晰,反演的深度也更接近真实深度,同时较大程度上克服了Tilt_depth方法受叠加异常的影响.理论模型验证了方法的有效性,并通过准噶尔盆地某区块磁异常数据的处理,提取了受强大的区域背景场掩盖的石炭系火成岩产生的弱异常,突出了构造分区、断裂分布等信息,获取了火成岩的位置及埋深参数,为该区的火成岩油气藏勘探提供了有效的处理途径.     

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

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

磁梯度张量解析信号分析法及其在场源位置识别中的应用

相对于传统磁异常数据,磁梯度张量数据可以提供有关异常体更多的信息,分辨率更高,具有广阔的应用前景.本文给出了基于磁梯度张量解析信号进行边界增强的新方法,该方法是通过对异常垂向一阶导数与z方向解析信号的比值做归一化得到的,其零值线位置可以很准确的反映出地质体边界位置,而且可以有效降低噪声的干扰.文中还证明了异常垂向一阶导数与x、y、z三个方向解析信号的比值同样满足欧拉齐次方程,且在计算过程中不需考虑构造指数N的影响,避免了因构造指数不当而引起的反演误差.通过对单一地质体及组合地质体模型的实验证明:与常规欧拉反褶积法相比,该反演方法能够更好的得到地质体界面及深度信息,所得的解更集中.将其应用到保定实测数据中,获得了更精确的场源信息.     

基于时变重力数据的川滇地区场源模型反演和解释

时变重力信号能反映地球内部介质迁移引起的不均衡变化,本文利用三维欧拉反褶积方法,通过建立合适的场源模型获得最佳反演参数,并采用水平梯度滤波法消除部分虚假发散解,优化反演结果.在此基础上,利用川滇地区2015-2019年流动重力重复观测资料,分析了不同时空尺度的区域重力场动态变化特征,开展对川滇地区重力变化场源特征的定性...     

高精度重力测量在金矿采空区探测中的应用研究

为了确定金矿采空区的位置,对其隐患进行治理,首次应用高精度重力方法.在研究区布设了两条比例尺为1∶100的高精度重力测量剖面,使用加拿大Scintrex公司的CG-5AutoGrav型重力仪, 测取重力异常数据,并采用欧拉反褶积方法实施反演解释,确定了采空塌陷区域及稳定地段.研究结果表明:高精度重力测量应用于金矿采空区探测是行之有效的,用欧拉反褶积方法解释重力异常数据是一种快速、高效、人为干预少的好方法.     

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

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

欧拉反褶积与解析信号相结合的位场反演方法

由于解析信号具有不受(二维)或少受磁化方向影响,能够较好反映磁性体边界的特性,因此受到人们的重视.欧拉反褶积法可以确定场源的位置和深度以及形状因子,具有较强的适应性.因此前人提出将二者相结合的方法.针对前人提出的方法中存在受高频干扰严重的问题,本文提出低阶的欧拉反褶积与解析信号相结合的位场反演方法.本方法在反演中只需计...     

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

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

归一化总水平导数法在位场数据解释中的应用

本文提出归一化总水平导数法,通过对总水平导数进行空间归一化计算实现了异常体水平位置和深度的估计,此外还推导出基于归一化总水平导数的欧拉反褶积法来估算地下地质体的空间位置,两种方法反演结果的相互验证可有效地提高反演结果的可信度.理论模型试验证明空间归一化总水平导数法和归一化总水平导数欧拉反褶积法均能有效地完成异常体的水平位置和深度的估计,所获得的位置参数与理论值相一致.在利用归一化总水平导数法进行磁异常解释时,对数据进行化磁极计算可得到更加准确的结果.将其应用于实际航磁数据的解释,获得了岩脉的大致分布特征.     

南祁连盆地重力异常特征及其地质意义

南祁连盆地具有丰富的油气资源,盆地内部坳陷边界、基底深度及内部构造控制着油气资源的分布.为落实南祁连盆地内部坳陷及构造分布,笔者充分利用多年来在研究区及周缘完成的各类勘探资料,以南祁连盆地的重力数据为主,结合相关资料,开展重力数据处理和综合反演解释,分析提取反映深大断裂的重力异常信息,并基于重力异常完成主要坳陷的密度界面深度反演和剖面2.5维反演,推断南祁连主要沉积盆地和坳陷的基底深度及起伏变化,进而综合讨论研究区主要坳陷的油气远景,为南祁连盆地油气资源战略调查和进一步勘探选区提供了依据.     

空-地-井重力异常正则化协同密度反演方法

围绕综合利用空-地-井多源重力异常联合反演提升精度的目标,提出正则化协同密度反演方法,从而有效利用多维数据的横向和纵向变化特征提高反演精度,且无需先验信息的约束.此外,还提出利用奇异值谱和深度分辨率工具评价航空数据观测高度和层数对反演分辨率的影响.通过理论模型试验本文方法在空-地、地-井等不同情况下的应用效果,结果表明空-地-井重力异常正则化协同密度反演方法能获得更高分辨率、高精度的反演结果,且证明不同的钻孔位置对反演结果有较大的影响,从而可指导实际的空-地-井联合勘探.最后,利用文顿盐丘实际勘探数据进行空地和空地井重力数据协同反演,反演结果垂向分辨率明显优于地面观测数据反演结果,盐盖的顶面及中心埋深与前人研究和地质资料解释相吻合,验证了方法正确性及实用性,为推进我国空-地-井立体重磁勘探提供了重要的技术手段.