论文部分内容阅读
重力正演是已知地下物性求地上观测面的重力异常值,其反演问题是根据地上观测面的重力异常值求地下场源体的物性以及形态。重力异常反演在地球物理勘探中有着十分重要的作用,是了解地质体场源特征最便捷的途径。采用长方体剖分单元求物性的三维反演技术,因模型参数易于操作、能模拟任意复杂地质体、可使用分离格架技术减少存储量等优点,已成为重力反演的一个主要技术手段。对于线性或者非线性的物性反演方案,即使使用了分离格架等技术,正演的计算量仍然在整个反演流程中占很大比重,当反演达到一定规模,整个反演计算99%的计算量都是由正演计算产生的,因此加快正演计算速度对整个反演流程有重要意义。本文利用NVIDIA CUDA并行计算平台,实现了基于GPU的三维重力正演与反演计算。对于正演计算,本文研究了几何格架的等效存储技术,在此基础上首次提出了几何格架的空间局部性,并应用此性质对GPU正演计算程序进行深入优化,获得了普通版本的GPU正演计算程序的4倍加速。对于反演计算,阐明了基于深度加权的反演方法是一个严格凸二次规划问题,并在CUDA异构计算平台上设计和实现了基于共轭梯度法的整个反演流程,给出了数值模拟结果。