MT Occam并行反演方案及性能分析.pdfVIP

维普资讯 第 29卷 12 武 汉 理 工 大 学 学 报 VO1．29 No．12 2007年 12月 JOURNALOFW UHAN UNIVERSITYOFTECHNOI~ ,Y Dec．2007 MTOccam并行反演方案及性能分析 ／刘 羽 (棒林工学院电子 r算机系，桂林 541004) 摘 要： 针对MTOccam反演速度慢的问题，分析 了实现其并行计算的可行性和途径，提 出了基于频点和拉格朗日乘 子( )扫描点划分的大粒度并行反演方案和具体实现步骤。对方案的加速性能进行了理论分析，得出了理论加速 比 曲线，明确了方案的适用性。采用基于消息传递的PVM并行平台，以主一从编程模式，实现了偏导数矩阵和 值求取的 并行计算。算例表明，综合并行效率较高，反演速度提高显著。 关键词： MrI、Occam反演； 并行计算； 机群 中图分类号： II393 文献标识码 ： A 文章编号 ：1671—4431(2007)I2一O136—05 ParallelMTOccam InversionSchemeandItsPerformanceAnalysis LIU Yu (I)epartmentofElectronicandComputerScience，C-uilinUniversityofTechnology，Guilin541004，China) Abstract： IhispaperanalyzedthefeasibilityandwaysforparallelimplementationoftheMT(magnetotelluric)Occamin— versioninordertOincreasetheinversionspeed．A coarse-grainedparallelinversionschemebasedonfrequencypointandLa— grangemultiplier( ),scanningpointdix，isionwasprotxosde andaconcreteparallelprocedureWaSputforward．Theoretical performanceoftheschenlewasanalyzedandthespeed—upcurveobtained，andtheapplicabilityoftheschemewasclarified．Par— aim computationsfo thepartialderivativematrixandthe!tseekingwererealizedbyusingthemaster—slaveprogrammingmodel underIheures,sage-paxsingPVM parallelplatform ．Theexperimentresultindicatedthattheall—aroundparallelefficiencyWS idealandtheinversionspeedincreasedobviously． Keywords： M I、Occam inversion； parallelcomputation； COW MTOccam反演采用了特有的自适应算法求取拉格朗 日乘子 ( )，方法稳定，且无需更多的经验，但 速度很慢。文献[1]对方法的 值求取策略提出了一些改进措施，提高了计算速度，但 一定程度上改变了 Occam反演机理，且方案中的初始 值和变化步长 需凭经验确定 ，合理性还值得探讨l2J。提早结束极为 耗时的模型光滑过程 (反演的后段)也可以缩短计算时问，但此时反演接近收敛而实际并未收敛。Occam反 演文档指出，如果要一味获得快的反演速度，不宜采用Occam反演方法，一些试图通过改变方法本身和收敛 条件来提高速度的做法并不可取，因为那样得不到最光滑的模型，与Occam反演思想相悖。 通过并行计算可以在大幅提高反演速度的同时保持MTOcca