重力异常及其梯度张量数据快速解释技术研究

被引量 : 0次 | 上传用户:tjh2088
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
重力勘探数据由重力异常数据及其梯度数据集构成,可以通过相关重力仪或重力梯度仪在卫星、航空、地面和海洋等重力勘探过程中获得。重力勘探数据从不同侧面反映了产生重力异常的场源的空间分布形态和密度变化规律,是由地球物理场现象解释地质对象的基本数据和研究出发点。为了更好地利用数据进行快速解释,通常需要针对不同的地质研究对象,选取不同的梯度张量分量或结合重力异常形成数据组合,建立与之相适应的解释方法。不同的梯度张量分量反映了同一地质体在不同方向的变化规律,其分辨能力和探测能力也随方向变化而存在差异。因此,本文在理论和实践环节开展解释方法技术研究工作:研究不同数据类型和组合的地质目标发现能力;研究不同数据组合针对场源的快速解释方法。首先从理论上,通过模型正演,对重力异常数据和梯度数据进行了评估分析,为解释工作中对重力异常以及梯度张量分量的选取提供依据。在重力垂向一阶导数和重力异常之间关系的对比研究基础上,对其他张量分量与重力异常分辨能力以及各个张量分量之间的分辨能力展开了系统对比分析。研究结果表明重力梯度张量分量在相对浅部范围内具有更高的分辨能力,而相对深部探测中并不具有优势,高精度梯度张量分量的优势主要在于可以提供更多高频的细节信息,从而更加详细地对地质体进行解释。另外,从应用角度,文中分析了测线布设位置和测线间距对重力异常及其梯度的影响,发现随着测线间距的增大和测线位置偏离地质体时,梯度张量分量仍然可以较完整地给出场源的信息,而重力异常及计算梯度异常均不能提供较完整的信息。但相较于重力异常,计算所得的张量梯度分量可以从不同侧面对地质体进行更加精细的描述。可以看出,在重力异常及其梯度数据的解释工作中,梯度张量数据以及张量多分量联合解释方法具有更大的优势。在传统解释方法基础上,本文开展快速解释方法的研究,强调同时利用更多的张量分量,优化梯度张量分量组合,从而推进多类数据联合解释技术的完善,加强解释过程中的信息融合。重力异常及其梯度数据的解释工作,通常需要从数据中获取有关场源的三个方面信息:水平位置、深度范围和物性参数。因此本论文主要从这三种参数的确定展开新方法的研究。水平位置的确定方法通常利用边界识别方法;深度范围参数的确定需要深度计算方法来完成;物性参数的获取除了传统物性反演之外,可由快速成像方法给出与物性相关的量值。近年迅速发展的边界识别方法由于可以快速且准确地确定地质体边界位置、构造水平位置等而被广泛关注。分析表明,最常用的传统方法,如Tilt梯度(斜导数)、Theta图等方法,在深部地质体边界位置确定中误差较大,且对细节信息的提取能力不足。文中首先对二者的相互关系进行分析,论述了这两种方法的相关关系,将其统一归类为“角度边界识别方法”,在此基础上,建立其他相类似识别方法,如反正弦等角度边界识别方法。由模型计算和对比分析发现,这些方法在地质体边界识别中效果基本相同,对解释对象的细节信息提取能力不足。为了改善细节信息提取能力,提出了基于标准差计算的改进型的“角度边界识别方法”,其计算过程是:先求取重力异常三个方向梯度在窗口内的标准差,然后再按照常规角度边界识别方法求取结果。模型检验和分析表明,改进后的方法可以探测识别更多和实际情况相吻合的细节特征。将该方法应用于实际数据处理中,取得了良好的效果。为了改善传统方法在深度地质体边界识别中的偏移问题,提出了“平面全张量梯度角度边界识别方法”:将水平面视作整个与垂向方向相对的分量,根据梯度分量矢量特征,类似于角度边界识别方法进行比值计算。模型计算结果表明,由于利用了更多的张量分量,使得边界位置的确定更加准确。文中还对反双曲斜导数边界识别方法在复杂的正、负剩余密度同时存在时产生虚假边界的缺陷进行了改进,给出了重力异常和梯度分量的数据联合解释的一种途径。模型和实测数据处理表明,改进后的方法在应用中更具优势。重力异常场源深度计算的方法经历了漫长的发展历史,形成了针对不同的场源类型和地质研究对象的不同深度计算方法。但是,在实际应用中,目前还无法找到一种包罗万象的万能方法来解决复杂的地质问题,每一种方法均有完善和改进的研究空间。文中首先基于Tilt-depth(斜导数深度计算方法)方法,推导建立了可应用于球体(点源)等三度体模型的重力异常深度计算方法,理论模型表明了该方法具有较准确的计算结果。为了更进一步体现张量分量解释手段的优势,将张量曲率矩阵引入到深度计算中,建立了基于张量曲率矩阵的深度计算方法,该方法类似于Tilt-depth方法的计算,只是将Tilt-depth方法中垂向导数替换为张量曲率矩阵的特征值,提高了水平分辨能力,从而使得深度计算结果也更加准确。另外,Tilt-depth方法以及张量曲率矩阵深度计算方法均是单独针对正或负剩余密度体进行计算,为了更快捷地对复杂地质情况进行计算,同样利用了重力异常和梯度相结合的办法,改进了张量曲率矩阵的特征值,从而可以对正、负剩余密度体同时存在的情况直接计算深度。其后建立了基于全张量不变量的深度计算方法,利用重力异常和全张量不变量的比值求取深度范围,发挥了张量不变量的优势特征,可以简单快捷地计算出地质体的埋深。建立了不同的模型试验并选取了多个实际重力数据对这些方法分别进行检验,均可以得到很好的应用结果,表明这些新方法的优越性和实用性。快速成像方法在物性参数的计算中具有快速、耗内存小等优势,此类方法虽然不能给出真正的物性参数,但这种与物性相关的参数在地质解释中具有重要的价值。因此,本文针对传统的极值点深度计算方法(DEXP)进行了相应的评价,对其优点、误差产生的原因等进行了分析,发现DEXP方法在针对重力一阶导数计算成像时误差更小。而广义线性反演成像方法是一种基于垂向导数的成像方法,因此将DEXP方法深度加权因子引入到广义线性反演成像计算中,使得广义线性反演成像可以得到更加合理、准确的成像结果。为了更进一步发挥梯度张量分量解释的优势,在广义线性反演成像基础上推导建立了梯度张量的广义线性反演成像方法。模型试验和实测数据处理表明其分辨能力更高,结果也更加准确。这些新的解释方法旨在实现重力异常及其梯度数据的快速解释工作,文中的模型数据和实际数据分析表明,利用这些方法可以从水平位置、深度范围和物性参数的计算三个方面对场源体进行快速的解释,且均可以得到高精度的结果,说明文中所建立的方法具有很好的应用价值。
其他文献
本文是基于日本MMG学派分离模型思想建立的冰区航行的船舶运动数学模型。《STCW公约》马尼拉修正案中明确提出了新时代对航海模拟器的要求——冰区航行将是航海模拟器中不可
镉(Cd)是动植物非必需的有毒重金属,容易通过土壤-植物体系迁移并进入食物链,对人体健康构成威胁。柑橘作为我国栽植面积最大的果树,部分果园土壤已受到Cd威胁,柑橘的多年生
马克思基于时代顺应时代并超越时代创立了马克思主义的意识形态领导权思想,马克思主义意识形态领导权思想是基于时代的领导权,其根本宗旨是为了人的自由解放,马克思主义意识形态
提出一种适用于同类别的图像子集的类别判断方法。同类别的图像子集经过PCA主特征提取后,选择较大的p个特征值对应的线性无关的特征向量,组成特征矩阵,则同类别的图像子集可
河南省自2008年8月开始实施新课程改革,作为第一批尝试新课程教学的高一教师来说,大家在新奇与兴奋的同时也明白即将面临诸多挑战,当年我也是其中的一员。新课程强调学校要致
目的 考察地黄在加工炮制过程中糖类成分的变化。方法 在加工炮制的 4个不同阶段 (6 5℃烘焙 0 ,1和6 d,以及烘干样品切片后常压蒸制 4 h)取样 ,热水提取。色谱条件为 :Sug
压力性溃疡(pressure ulcer,PU)又称压疮,是指身体局部软组织持续受压、血液循环障碍,导致组织细胞缺血、低氧、营养代谢障碍而发生变性、坏死[1]。常见于老年患者、长期卧床患
从2004年《中华人民共和国农业机械化促进法》公布实施以来,到2013年这十年我国农业机械化发生了历史性巨变,这是我国农机化发展史上很不平凡的十年,凸显出四大特点。十年巨
近代中国体育建筑经历了从被动植入到内化吸收、从租界启蒙到华界壮大的发展过程,目前国内针对此阶段体育建筑的研究较为匮乏且价值认知不足。通过探索近代中国体育建筑在四
局部二值模式(LBP)是一种对灰度图像的纹理进行描述的方式。然而经典的LBP算法存在数据使用不充分的问题,因此提出一种分层的LBP算子(HLBP)。HLBP选取多层阈值对图像进行LBP序列