[发明专利]基于复杂地形的重力场数值模拟方法、装置和计算机设备

说明书

本申请涉及一种基于复杂地形的重力场数值模拟方法、装置、计算机设备和存储介质。所述方法包括：构建包含复杂地形和勘探目标的三维长方体模型，进行网格剖分，得到多个小长方体，设定小长方体的空间域密度值；根据网格剖分的网格参数和预设的高斯参数，得到偏移波数，根据地形信息，将三维长方体模型分为下部分模型和上部分模型，分块叠加计算多个观测面的频率域重力异常场，进行二维高斯傅里叶反变换，得到多个观测面的空间域重力异常场，通过样条插值算法得到复杂地形起伏表面的空间域重力异常场。本发明垂直方向保留在空间域，网格剖分灵活，适合模拟复杂地形及复杂地质体，多个偏移波数便于并行计算，提高了重力异常计算效率。

技术领域

本申请涉及计算机数值模拟技术领域，特别是涉及一种基于复杂地形的重力场数值模拟方法、装置、计算机设备和存储介质。

背景技术

现阶段，由于重力观测数据的增加，重力场正演所需的内存和计算的急剧增长，对于普通计算机难以实现，已经成为三维密度快速成像的难点问题。

目前，空间域方法以其准确性仍然占据主导地位，但是随着标准FFT扩边法和NUFFT在重力异常正演数值模拟中的应用，频率域方法以其简洁性、准确性和高校性成为处理大规模正演的首选方法。文献(Chai Y , Hinze W J . Gravity inversion of aninterface above which the density contrast varies exponentially with depth.Geophysics, 1988, 53(6):837-845.)采用偏移抽样方法有效提高了傅里叶变换的精度，使频率域方法逐渐应用到复杂模型正演模拟中。文献(Tontini, F.C., L. Cocchi, C.Carmisciano. Rapid 3-D forward model of potential fields with application tothe Palinuro Seamount magnetic anomaly (southern Tyrrhenian Sea, Italy).Journal of Geophysical Research, 2009. 114.)给出了任意密度分布情况下的重力异常频率域表达式。

现有技术中，大尺度大规模的航空重力测量已经兴起，但实际航空测量中地形并非是一个平面，而是复杂的起伏地形，传统空间域方法和频率域方法难以满足海量高分辨率数据精细快速反演要求，存在适应性不佳的问题，因此研究一种基于复杂地形的三维重力异常高效高精度数值模拟方法具有现实意义。

发明内容

基于此，有必要针对上述技术问题，提供一种能够适用于复杂地形重力场数值模拟的基于复杂地形的重力场数值模拟方法、装置、计算机设备和存储介质。

一种基于复杂地形的重力场数值模拟方法，所述方法包括：

根据复杂地形的地形信息和勘探目标的目标信息，构建包含所述复杂地形和所述勘探目标的三维长方体模型，将所述三维长方体模型在x轴方向、y轴方向和z轴方向进行网格剖分，得到多个小长方体，每个所述小长方体在x轴方向和y轴方向的尺寸相同，并根据所述地形信息和所述目标信息设定所述小长方体的空间域密度值；

根据所述三维长方体模型网格剖分的网格参数和预设的高斯参数，得到偏移波数；

根据所述空间域密度值和所述网格参数得到所述三维长方体模型的空间域重力位积分表达式，对所述空间域重力位积分表达式沿水平方向进行二维高斯傅里叶变换得到所述三维长方体模型的频率域重力位积分表达式；

根据所述地形信息，将所述三维长方体模型分为下部分模型和上部分模型，根据所述偏移波数分别计算所述下部分模型的第一垂向一维积分和所述上部分模型的第二垂向一维积分；所述下部分模型是复杂地形最低点以下部分的长方体模型，所述上部分是所述最低点以上部分的长方体模型；