基于GPU的任意三维复杂形体重磁异常快速计算

作者：陈召曦; 孟小红; 刘**; 郭良辉

来源：物探与化探, 2012, (01): 117-121.

• 重磁异常正演计算
• 任意三维复杂形体
• GPU并行计算 forward calculation of gravity and magnetic anomalies
• 3D arbitrary complicated bodies
• GPU parallel computation

提出了基于图形处理单元的任意三维复杂形体的重磁异常快速正演计算方法。将地下半空间剖分为大小相等规则排列的一组长方体单元,任意三维复杂形体可以表示成很多不同体积和密度(磁性)的长方体的近似组合。用解析方法计算出所有这些长方体在计算点的重力(磁力)异常,并累加求和,就可以得到整个模型体在计算点引起的重(磁)异常值。为了提高近似程度,需将地下半空间剖分得很细,用传统的CPU串行程序计算相当耗时。GPU在处理能力和存储器带宽上相对CPU有明显优势,采用GPU并行算法,可大大提高计算速度。相关试验结果表明,用GPU实现的正演快速算法计算结果正确,效率明显提高,为重磁异常三维物性反演提供了基础。