| 摘 要： | 密度界面反演方法在油气勘探、推断区域构造、获取地壳结晶基底面、莫霍面起伏形态等方面具有重要意义。现有的密度界面反演方法大多基于笛卡尔坐标系统,当涉及到大区域乃至全球尺度的密度界面反演时,地球曲率的影响将不可忽视,需考虑基于球坐标系Tesseroid模型的密度界面反演方法。然而,受计算精度和效率制约,已有的基于Tesseroid模型的密度界面反演方法并不能很好地适用于地表重力观测数据的反演计算。笔者基于前人研究,给出了一种适用于地表观测数据的球坐标系密度界面反演方法。该方法首先将常规的球坐标系高斯—勒让德重力积分公式进行简化,提高了重力正演计算效率。随后,引入并改进了前人的自适应剖分方案,提高了重力正演计算精度。在此基础上,采用Cordell迭代优化算法,得到了适用于地表观测数据的球坐标系密度界面反演方法。通过模型数据试验,对本文球坐标系密度界面反演方法进行了验证。结果表明,笔者对高斯—勒让德积分公式加以改进和引入改进的自适应剖分方案后,很好地克服了计算精度和效率对地表观测重力计算的掣肘,并且,基于球坐标系的密度界面反演结果优于基于笛卡尔坐标系的密度界面反演结果。在华南大陆实际数据试验中,利用文中方法得到的华南大陆的莫霍面深度与前人得到的莫霍面深度具有较高的吻合度,莫霍面自西向东逐渐抬升,西部抬升剧烈,东部抬升平缓,沿武陵山—贵桂交界一带呈现明显的NNE向莫霍面梯级带,验证了文中方法的科学有效性。
|
