提出了基于图形处理单元的任意三维复杂形体的重磁异常快速正演计算方法。将地下半空间剖分为大小相等规则排列的一组长方体单元,任意三维复杂形体可以表示成很多不同体积和密度(磁性)的长方体的近似组合。用解析方法计算出所有这些长方体在计算点的重力(磁力)异常,并累加求和,就可以得到整个模型体在计算点引起的重(磁)异常值。为了提高近似程度,需将地下半空间剖分得很细,用传统的CPU串行程序计算相当耗时。GPU在处理能力和存储器带宽上相对CPU有明显优势,采用GPU并行算法,可大大提高计算速度。相关试验结果表明,用GPU实现的正演快速算法计算结果正确,效率明显提高,为重磁异常三维物性反演提供了基础。
This paper presents a fast forward calculation method of gravity and magnetic anomalies for 3D arbitrary bodies based on GPU.The surface is subdivided into a large number of regular voxels.For any 3D arbitrary body,the total volume can be combined by these voxels with fixed size and unknown physical properties(density or susceptibility).The gravity(magnetic) anomaly of these voxels at the surface points can be computed by the analytical method.The sum of the above gravity(magnetic) anomaly of each voxel can represent the anomaly of the whole arbitrary body.For the purpose of improving the approximation,the subsurface should be divided into a larger number of voxels.The authors have adopted GPU parallel method for forward calculation because the traditional CPU program is time-consuming.The result shows that the forward calculation by GPU has correct results and improved efficiency.The idea can provide the basis for 3D inversion of gravity and magnetic anomaly.