重力和重力梯度数据联合反演方法研究

来源 :吉林大学 | 被引量 : 0次 | 上传用户:zhoudddmmm
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
重力数据的处理和解释一直是地球物理数据处理和解释中的重要组成部分。近年来,随着全张量重力梯度测量系统(FTG)的发展及应用,重力梯度数据的地位也越来越重要。从理论上说,重力数据是重力位在垂直方向上的偏导数,重力梯度数据则是重力位在各方向的二次求导。在求导计算过程中,低频信号被压制,高频信号被增强。所以,在频率域比较重力数据和重力梯度数据,结果表明,重力数据包含较多的低频信息,重力梯度数据则包含较多的高频信息。因此,将重力和重力梯度数据联合,可以实现信息的互补,有效的降低反演中的多解性。基于这一点,本文针对重力和重力梯度数据联合反演的问题展开讨论。首先,对重力和重力梯度数据进行量化比较分析。通过重力梯度数据在频率域的表达式,推导重力梯度各分量之间的关系,发现在重力梯度各分量中,zzg分量包含的信息更多;随后,通过对比不同深度模型的zg分量和zzg分量的振幅谱,分析zg分量和zzg分量随深度的变化;分别对zg和zzg分量进行反演计算,结果表明,zg反演结果在深部的分辨率较高,zzg反演在浅部具有较高的分辨率。我们通过再加权反演的方法将深度加权函数结合到稳定函数部分,在仅考虑深度加权,不加任何其他约束条件的情况下,分别对重力场和重力梯度数据各分量进行反演计算,依据各分量的反演结果,对各分量在反演计算中的表现有了一个直观的印象,对各分量所包含信息的特点进行对比。随后,我们针对重力和重力梯度联合反演问题深入讨论。与单一分量反演相比,多分量联合反演需要消耗更多的计算时间和内存。首先,在反演计算开始前,我们需要计算并存储各分量的灵敏度矩阵。因为分量种类的增加,存储和计算灵敏度矩阵所消耗的内存空间和计算时间也增加。因此,我们提出了一种快速计算灵敏度矩阵的方法,能够有效的减少计算时间,并基于此提出了一种存储灵敏度矩阵的策略。之后,我们对反演过程中深度加权函数的选择展开了讨论。我们列举了几种深度加权函数,针对其在重力和重力梯度反演中的表现进行了对比,最后选择了基于异常体深度信息的加权函数。这种深度加权函数依据异常体的埋深设置权值,与数据类型无关,因而适用于重力和重力梯度数据联合反演。我们给出了再加权光滑反演的计算流程及相应的公式推导。基于再加权的光滑反演,我们对深度加权函数对异常体埋深信息准确性的要求进行了讨论,发现在埋深信息与实际信息存在一定差异时,采用这种深度加权函数,仍能获得合理的反演结果。之后,我们将其应用到反演计算中,并将这种加权函数在反演中的表现与基于灵敏度矩阵的深度加权函数进行对比。结果表明,基于异常体深度信息的加权函数,因为将深度信息结合到了反演计算中,有效的提高了反演结果的分辨率。接下来,针对稳定函数的选择开展了研究,并对本文所采用的优化算法展开了讨论。针对不同的地质体,需要选择不同的稳定函数。本文中,我们主要针对具有陡峭边界的地质体,因此选择了最小梯度支撑函数作为构成联合反演正则化方程的稳定函数。本文中,我们主要采用非线性共轭梯度算法,这是一种适合解决大规模反演问题的算法。这种算法在三维电磁反演中有着广泛的应用,在重力和重力梯度反演中的应用则较少。通过模型试验,我们将非线性共轭梯度算法与BFGS拟牛顿法进行对比,发现与BFGS拟牛顿法相比,虽然非线性共轭梯度算法收敛速度较慢,但是其消耗的时间更短,对内存的要求更低。为了将密度约束结合到反演计算中,同时,确保不会对反演过程的稳定性造成影响,我们采用了不等式约束条件。我们建立了由多个不同异常体组成的模型,对重力梯度分量组合反演进行对比和讨论。结果表明,多分量组合能够有效的提高反演结果的分辨率,但是,当数据量达到一定程度后,虽然数据能够更好的拟合,反演结果收敛程度更好,但是,反演结果与理论模型的一致性降低,这就需要额外的信息来约束反演过程。综合考虑反演结果和反演效率,我们给出了反演结果和反演效率最佳的分量组合。随后,我们基于同一模型对重力和重力梯度联合反演进行了对比和讨论。结果表明:加入重力数据后,得到的联合反演的结果与重力梯度分量组合反演结果大体一致;对于其中的某些异常体,加入重力数据后,其反演结果的分辨率得到了提升。针对具有不同埋深的异常体,我们将空间梯度加权函数应用到了反演计算中。这一函数的应用,使得我们能够将前次的反演结果中的有效信息作为先验信息结合到下次的反演计算中。通过模型试验,我们发现,采用这种方法,反演结果的分辨率得到了提升。最后,我们将本文中的方法应用到了在美国路易斯安娜州的文顿盐丘测得的实际数据。利用重力和重力梯度数据进行三维反演计算,进而依据反演结果推断盐盖的分布。通过对不同分量组合反演结果对比,我们发现zg|xyg|xzg|yyg|yzg|zzg反演结果最好。我们给出了反演结果的三维分布图,与该地区已有的研究成果对比,可以得出,我们的反演结果是合理的。
其他文献
生源文化基础较差,"专业技能"与"文化基础"之间权重失衡,定位模糊以及教学评价侧重于专业课程是目前中职文化基础课教学存在的主要问题,中职学校应加强对文化基础课学科价值
<正> 随着压铸工艺的广泛应用,各种压铸合金随之而发展。在当前压铸工艺水平不断提高、压铸件质量要求日益严格的情况下,制定我国压铸合金技术条件的工作已迫在眉睫。 本文介
<正>【裁判要旨】货币为种类物,虽然权利人对货币的占有可以认定为所有,但在特定条件下,不能简单根据占有即认定为所有。对于一般账户中的货币,应以账户名称为权属判断的基本
期刊
近年来,全国性微课大赛与各级地方微课大赛相继举办,吸引了广大一线教育工作者的参与。但是,各类大赛中的评价体系差别较大、标准不同,同时在比赛中对微课创作的思路和制作的
就业能力是影响高校毕业生是否成功就业和职业发展的核心因素。从教育生态学角度来看,在毕业生就业能力培养方面,高职教育存在生态定位模糊、教育生态结构失衡、教育生态链不
光石沟铀矿床位于北秦岭东段的商丹地区,该矿床的主要赋矿岩石为黑云母花岗伟晶岩。通过对该区资料的整理综合研究发现:光石沟伟晶岩属于壳源深源岩浆分异伟晶岩;按岩浆源区
2016年3月全国人大四次会议审查通过了《中华人民共和国国民经济和社会发展第十三个五年规划纲要》,对我国的2016-2020年的发展方向进行了规划。该规划如何影响我国精炼铜的
本文结合某公司氧化铝的供货质量情况及目前氧化铝的质量标准,通过对氧化铝的化学成分和物理性能两个方面进行比较分析,阐述了氧化铝的理化性能指标对铝电解生产过程和产品质