论文部分内容阅读
重力数据的处理和解释一直是地球物理数据处理和解释中的重要组成部分。近年来,随着全张量重力梯度测量系统(FTG)的发展及应用,重力梯度数据的地位也越来越重要。从理论上说,重力数据是重力位在垂直方向上的偏导数,重力梯度数据则是重力位在各方向的二次求导。在求导计算过程中,低频信号被压制,高频信号被增强。所以,在频率域比较重力数据和重力梯度数据,结果表明,重力数据包含较多的低频信息,重力梯度数据则包含较多的高频信息。因此,将重力和重力梯度数据联合,可以实现信息的互补,有效的降低反演中的多解性。基于这一点,本文针对重力和重力梯度数据联合反演的问题展开讨论。首先,对重力和重力梯度数据进行量化比较分析。通过重力梯度数据在频率域的表达式,推导重力梯度各分量之间的关系,发现在重力梯度各分量中,zzg分量包含的信息更多;随后,通过对比不同深度模型的zg分量和zzg分量的振幅谱,分析zg分量和zzg分量随深度的变化;分别对zg和zzg分量进行反演计算,结果表明,zg反演结果在深部的分辨率较高,zzg反演在浅部具有较高的分辨率。我们通过再加权反演的方法将深度加权函数结合到稳定函数部分,在仅考虑深度加权,不加任何其他约束条件的情况下,分别对重力场和重力梯度数据各分量进行反演计算,依据各分量的反演结果,对各分量在反演计算中的表现有了一个直观的印象,对各分量所包含信息的特点进行对比。随后,我们针对重力和重力梯度联合反演问题深入讨论。与单一分量反演相比,多分量联合反演需要消耗更多的计算时间和内存。首先,在反演计算开始前,我们需要计算并存储各分量的灵敏度矩阵。因为分量种类的增加,存储和计算灵敏度矩阵所消耗的内存空间和计算时间也增加。因此,我们提出了一种快速计算灵敏度矩阵的方法,能够有效的减少计算时间,并基于此提出了一种存储灵敏度矩阵的策略。之后,我们对反演过程中深度加权函数的选择展开了讨论。我们列举了几种深度加权函数,针对其在重力和重力梯度反演中的表现进行了对比,最后选择了基于异常体深度信息的加权函数。这种深度加权函数依据异常体的埋深设置权值,与数据类型无关,因而适用于重力和重力梯度数据联合反演。我们给出了再加权光滑反演的计算流程及相应的公式推导。基于再加权的光滑反演,我们对深度加权函数对异常体埋深信息准确性的要求进行了讨论,发现在埋深信息与实际信息存在一定差异时,采用这种深度加权函数,仍能获得合理的反演结果。之后,我们将其应用到反演计算中,并将这种加权函数在反演中的表现与基于灵敏度矩阵的深度加权函数进行对比。结果表明,基于异常体深度信息的加权函数,因为将深度信息结合到了反演计算中,有效的提高了反演结果的分辨率。接下来,针对稳定函数的选择开展了研究,并对本文所采用的优化算法展开了讨论。针对不同的地质体,需要选择不同的稳定函数。本文中,我们主要针对具有陡峭边界的地质体,因此选择了最小梯度支撑函数作为构成联合反演正则化方程的稳定函数。本文中,我们主要采用非线性共轭梯度算法,这是一种适合解决大规模反演问题的算法。这种算法在三维电磁反演中有着广泛的应用,在重力和重力梯度反演中的应用则较少。通过模型试验,我们将非线性共轭梯度算法与BFGS拟牛顿法进行对比,发现与BFGS拟牛顿法相比,虽然非线性共轭梯度算法收敛速度较慢,但是其消耗的时间更短,对内存的要求更低。为了将密度约束结合到反演计算中,同时,确保不会对反演过程的稳定性造成影响,我们采用了不等式约束条件。我们建立了由多个不同异常体组成的模型,对重力梯度分量组合反演进行对比和讨论。结果表明,多分量组合能够有效的提高反演结果的分辨率,但是,当数据量达到一定程度后,虽然数据能够更好的拟合,反演结果收敛程度更好,但是,反演结果与理论模型的一致性降低,这就需要额外的信息来约束反演过程。综合考虑反演结果和反演效率,我们给出了反演结果和反演效率最佳的分量组合。随后,我们基于同一模型对重力和重力梯度联合反演进行了对比和讨论。结果表明:加入重力数据后,得到的联合反演的结果与重力梯度分量组合反演结果大体一致;对于其中的某些异常体,加入重力数据后,其反演结果的分辨率得到了提升。针对具有不同埋深的异常体,我们将空间梯度加权函数应用到了反演计算中。这一函数的应用,使得我们能够将前次的反演结果中的有效信息作为先验信息结合到下次的反演计算中。通过模型试验,我们发现,采用这种方法,反演结果的分辨率得到了提升。最后,我们将本文中的方法应用到了在美国路易斯安娜州的文顿盐丘测得的实际数据。利用重力和重力梯度数据进行三维反演计算,进而依据反演结果推断盐盖的分布。通过对不同分量组合反演结果对比,我们发现zg|xyg|xzg|yyg|yzg|zzg反演结果最好。我们给出了反演结果的三维分布图,与该地区已有的研究成果对比,可以得出,我们的反演结果是合理的。