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

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

欧拉反褶积法在航空重力勘探中的应用

在论述欧拉反褶积基础理论与方法的基础上,以我国某海域航空重力资料为例,对该海域内地质异常体(侵入岩、断鼻构造)及断裂构造等进行了欧拉反褶积反演计算,获得了场源深度等信息,确定了场源位置及边界,并与航空布格重力异常、垂向导数、上延、水平梯度模、剩余异常以及重磁联合反演结果等资料相互验证,取得了很好的应用效果.基于该地区的反演效果,证明欧拉反褶积法可应用于航空重力勘探中,可为其他航空重力勘探区计算场源位置及深度信息提供可靠并且准确的研究方法,从而达到辅助航空重力勘探测量的效果.     

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

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

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

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

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

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

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

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

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

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

重力全张量数据联合欧拉反褶积法研究及应用

全张量测量技术是在空中或海上用加载了多个加速度计的移动平台技术测量位场的五个独立分量.各张量分量包含不同方向的地下地质体信息,水平张量分量T_(xx)、T_(yy)、T_(xy)、T_(xz)、T_(yz)通常用于识别和映射与地质构造或地层变化有关的测量区域中的目标,垂直张量分量Tzz用于估计深度.然而,这些分量传统上是彼此分开解释,经常遇到错失关键信息的风险.本文所用全张量欧拉反褶积是在单独z方向的欧拉反演基础上发展而来的,它融合了重力异常垂直分量以及其三个方向导数、水平分量以及其三个方向导数.全张量数据信息得以有效应用的同时,欧拉反褶积结果也比常规欧拉反褶积结果更加收敛.最后,结合美国墨西哥湾地区实测航空FTG数据,用重力梯度张量数据进行联合欧拉三维反演研究,有效的识别岩盖的边界信息,划分岩盖范围,为进一步研究盖层底下深部复杂地质情况提供可靠的解释结果.     

欧拉反褶积中构造指数的计算与优化选取

构造指数是欧拉反演中最重要的参数,以往选取构造指数往往根据地质先验信息和人工经验试探选取,这样往往不能快速,准确地确定反演构造指数,从而影响反演效果,甚至反演失败.本文提出利用异常值自动估算构造指数,这样避免了经验选取构造指数产生的不确定因素的干扰,可以快速方便地选出较合理的构造指数,提高反演效率.模型实验显示了该方法良好的反演效果,证明了文中提出的方法的可行性.     

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

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

基于欧拉反褶积方法的航磁三分量应用研究

在中国北部某地区首次开展了航磁三分量矢量飞行测量.本次测量是国内第一次直接测量航磁异常三分量Z_a、H_(ax)、H_(ay)数据,即文中提到的垂直分量V、北分量N、东分量E.由于三分量磁测数据具有矢量空间性,因而凭借其全空间矢量,在反演过程中可以更直接地发挥指向性作用.更为丰富的测量资料参与反演,所得到的最终结果为构造解释提供了更有力的约束,从而有助于减少地质解释过程中的多解性,提升解释结果的可靠性.本文以模型试验为切入点,首先进行三个分量各自的欧拉反褶积处理模型试验,以模型试验结果总结了不同的分量经反演后各自对理论模型边界刻画的特点.而后依据以上特点,将三个分量各自的反演结果进行叠加.与总场的反演结果相比,三个分量反演的叠加结果改善了模型体边界的识别效果,获得了更为准确的模型体位置信息.在总结了模型试验相关规律的基础之上,按照类似步骤分别对实测航磁三分量数据进行了欧拉反褶积处理,再将三个分量实测数据各自的反演结果进行叠加,以叠加后的结果进行构造解释.利用不同分量各自的反演特点,综合发挥多分量的资料优势,最终以三分量反演的叠加结果为依据,有效地识别了区内断裂,更加准确地刻画出断裂的形迹位置,揭示了不同性质断裂的世代关系,显示了良好的应用效果.     

基于重力梯度张量数据的欧拉反演在梅山铁矿的应用

随着在石油勘查、矿产勘查、构造研究等地质工作中提出的应用位场异常研究规模相对较小的地质体及其细节的要求不断提高,传统重力方法在精度和分辨率方面已不能满足要求.重力梯度张量作为重力位二次导数,较之重力异常,能够更直观地反映地下密度体的赋存状态,基于重力梯度张量数据的反演能够提高地质体的定量模拟质量.但是由于受国内仪器的发展水平限制,直接获取实测的重力张量数据难度较大,目前对于重力张量数据的获取主要是通过数值计算的方法,把重力异常测量值变换成张量信息.本文是在样条理论的基础上,提出新的重力梯度张量计算方法,利用三次样条插值方法拟合重力异常,进一步地通过积分导数关系求取样条函数系数,实现了由重力异常求取重力梯度张量的数据计算.在此基础上,利用张量数据进行欧拉反演计算,对比常规欧拉反演结果,进一步分析张量欧拉反演精度.球体及水平板状体模型试验结果显示,张量欧拉反演更能有效地完成目标体的反演,其结果更加准确.最后,将张量欧拉反演应用到梅山铁矿区重力实测数据解释中,准确地反演出梅山主矿体的水平投影范围及矿体中心埋深.同时,相比常规欧拉反演,张量欧拉反演在刘家村附近、梅山村南部及黄林库—吴家洼地段显示出有更加收敛的欧拉解,这为梅山地区进一步找矿工作提供了重要的理论支撑.     

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

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

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

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

全张量磁梯度数据解释的均衡边界识别及深度成像技术

全张量磁梯度数据具有高精度、高分辨率、多参量的优点,能更加清晰地刻画地质体的分布特征,综合利用磁张量梯度数据准确地获得地质体水平位置和深度信息是解释的主要目的.磁张量数据的方向解析信号具有减小倾斜磁化干扰的优点,常被用来圈定磁源体的水平位置,但解析信号强度随着地质体埋深的增加急剧衰减,难以有效识别较深的地质体.张量数据均衡边界识别技术,利用不同方向解析信号的比值函数,能有效地均衡不同深度地质体的响应,同时显示不同深度地质体的边界,提高了对较深地质体的分辨率.磁张量数据深度成像技术根据实测张量数据与假定模型张量数据的相关系数来给定地质体的深度,综合利用多参量数据联合反演提高了反演结果的准确性,且无需进行复杂的反演运算,是大数据量张量数据解释的有效方法.理论模型试验证明：磁张量数据均衡边界识别技术可清晰和准确地识别地质体的水平范围,受倾斜磁化干扰小;磁张量数据深度成像技术可准确地获得地质体的深度信息,具有较强的抗噪性.将上述方法应用于铁矿区实测航磁张量梯度数据解释,获得了铁矿体水平分布与埋深,深度结果与张量欧拉反褶积法计算结果一致.     

磁源体参数反演的快速局部波数法

现有磁异常解释的局部波数法在计算地质体位置时大多是通过复杂的方程求解来完成,且需已知地质体构造指数信息,构造指数是描述地质体类型的参数,在实际数据解释中是难以确定的,因此现有方法计算结果的精度较低。本文提出磁异常解释的快速局部波数法,定义水平与垂直局部波数的平方和为和局部波数,利用和局部波数与其它局部波数的线性组合直接、快速地实现地质体位置参数和属性参数(构造指数)的计算,无需任何先验信息及方程求解运算。通过理论模型试验证明快速局部波数法能准确地完成磁异常的解释,且受背景异常和倾斜磁化干扰较小。将快速局部波数法应用于实际磁数据的解释,获得了地质体的位置和构造指数。     

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

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

郯庐断裂带江苏段及邻区航磁异常与深部构造特征

收集了郯庐断裂带江苏段及周边地区较大比例尺的航磁异常数据,通过水平及垂向求导、欧拉反褶积处理、居里面反演及剖面建模,系统分析了区域地壳磁性结构、断层空间展布及深部构造特征,并揭示了其与地热、地震以及岩浆活动的相互关系.研究表明:(1)区域航磁异常场具有明显的分区性,线性特征明显的航磁异常高值区基本由断裂主控,为岩浆沿断层侵入、灌入导致,而线性特征不明显的航磁异常高值或者低值区主要与磁性基底的埋深有关,其中郯庐断裂带形成NNE向的串珠状高磁性异常条带.(2)航磁场及其方向导数图显示区域发育17条主要断裂,其中郯庐断裂带与周边的10条断裂处岩浆活动剧烈;欧拉反褶积结果显示研究区磁性基底埋深在郯庐断裂带处及其两侧不同构造单元处存在差异;居里面反演结果表明郯庐断裂带为居里面坳陷区,该处为热流值及地温梯度低值区.剖面建模结果揭示郯庐断裂带底部可能存在大范围的基性及超基性岩浆岩,且推测断层剖面延伸形态为张扭性的"负花状构造".(3)以无锡—宿迁断裂为界,郯庐断裂带南北两侧在原始航磁场、航磁梯度场以及欧拉反褶积结果中均体现出了明显的分段性,表现为北侧航磁异常值更高、航磁梯度场波动更加剧烈且磁性基底埋深更浅.     

深度域随机模拟地震反演方法及其应用

常规时间域地质统计学反演方法适用于井控程度高的地区,对合成记录及层位标定要求比较高.对于开发程度较低的工区而言,由于井的密度不高,且可能存在标定效果不理想的井,这样会导致地质统计学反演的缺点被放大,影响储层砂体预测的效果.基于以上常规反演的弱点,本文通过对常规时间域内的地质统计学反演方法进行改进,形成一种新的深度域内的随机模拟反演方法,可有效消除标定误差对反演结果造成的影响.通过反演结果的对比分析,表明这种反演方法对储层的预测能力好于常规反演.     

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

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