基于MPI大地电磁二维正则化反演并行计算研究.pdfVIP

第35卷第1期 物探化探计算技术 2013年1月 文章编号：1001—1749(2013)01一0099一08 基于MPI大地电磁二维正则化 反演并行计算研究 刘鹏茂1’2，柳建新1，童孝忠1，孙丽影1 (1．中南大学地球科学与信息物理学院，长沙4l0083； 2．广州市市政工程设计研究院，广州 51∞60) 摘要：运用并行技术能提高运算效率，在大地电磁数据处理中，偏导数矩阵的计算占据了大量 的计算时间，使得反演速度较慢。这里在目标函数中采用基于阻尼高斯一牛顿法的正则化反演 方法，对偏导数计算进行并行处理，其余部份保留串行计算模式，针对不同频率间实现完全对等 的并行策略，模拟在多进程下的并行计算速度。通过两个理论模型试算，验证了并行算法的正确 性、高效性和稳定性。 关键词：大地电磁；MPI；正则化反演；高斯一牛顿法；并行计算 631．3+25 中图分类号：P 文献标识码：A 0前言 1 MPI简介 信息技术是地球物理勘探的支撑，因此，计算 MPI(MessagePassingInterface)能为高性能 机信息技术的发展极大地提高了地球物理领域并 并行计算提供一个方便的环境，具有移植性好，功 行计算技术的快速发展。随着电磁法数据处理的 能强大，效率高等多种优点。它是一个库，而不是 数据量越来越大，对于并行计算的条件也要求越来 一门语言，通过与FORTRAN等高级语言结合，实 越高，并行技术在电磁法数据处理中，将会起到越 现并行性。 来越重要的地位和作用。 MPI是由很多库函数构成的，利用MM并行 在国内方面，陈金窗等[1]在微机网络上实现了 工具进行数据处理时，并行设计的程序在采用何种 并行设计模式下，利用这些函数进行调用，进行通 2．5维CSAMT正演的并行计算；谭捍东等[9]实现 信。在计算过程中，每一个进程都是由很多通信领 了三维模型的大地电磁正演并行计算，随后又实现 了三维快速松弛反演的并行计算，证明了此方法的 域组成的，都是通过进行指定的通信域进行的。在 程序的运行过程中，开始是调用专门的函数MPI— 有效性和稳定性；刘羽[9J实现了二维大地电磁Oc— Init函数进行第一步进入MPI的环境中，当进入 cam反演的并行计算；王月英[1]实现了三维波动方 该环境的时候，首先进行初始化处理，并行程序再 程的有限元并行正演模拟。在国外方面，Zyser— man和Santos[1]实现了大地电磁三维有限元的并 的进程编号，然后调用MPI—Comm—size获得通信 行计算，这些都极大地促进了并行计算技术在地球 的进程个数，再调用其它的相应函数进行计算处 物理学界的发展。因此，通过并行技术会大大加快 理，当计算结束的时候，有专门的函数进行退出并 了地球物理数据处理的速度。 调查局科研项目(资[z011]03一01—64) 收稿日期：2012一06～04改回日期：2012—12—20 万方数据 100 物探化探计算技术