第29卷 第12期
2007 年 12 月
武 汉 理 工 大 学 学 报
JOURNAL OF WUHAN UNIVERSITY OF TECHNOLOGY
Vo l
.2 9
No
.1 2
Dec
.2 0 0 7
M T Occam
并行反演方案及性能分析
刘 羽
(桂林工学院电子与计算机系 ,桂林 541004)
摘 要 : 针对
MT Occam
反演速度慢的问题 ,分析了实现其并行计算的可行性和途径 ,提出了基于频点和拉格朗日乘
子(μ
-1
)扫描点划分的大粒度并行反演方案和具体实现步骤。对方案的加速性能进行了理论分析 ,得出了理论加速比
曲线 ,明确了方案的适用性。采用基于消息传递的
PVM
并行平台 ,以主-从编程模式 ,实现了偏导数矩阵和 μ值求取的
并行计算。算例表明 ,综合并行效率较高 ,反演速度提高显著。
关键词 :
MT Occam
反演 ; 并行计算 ; 机群
中图分类号 :
TP
393 文献标识码 :
A
文章编号 :1671-4431(2007)12-0136-05
Parallel MT Occam Inversion Scheme and Its Performance Analysis
LIU Yu
(
Departme nt of Electronican d Computer Scie nce
,
G uilin University o f Techno logy
,
Guilin
541004 ,
Ch i na
)
Abstract
:
This paper analy ze d the feasibility and ways for parall el imple mentation of the MT
(
magneto telluric
)
Occam in
-
ver s io n in order to inc rea se the inv ersio n s pe e d
.
Acoarse
-
grained parallel inversion scheme based on frequency point and La
-
grange multiplier
(μ
-1
)
scanning point division was proposed and a concrete parallel procedure was put forward
.
Theoretical
performan ce o f th e sche me was analy zed and the speed
-
up c u rveobtai n e d
,
and th e applicab ility of the sche me was clarified
.
Par
-
allel comp utation s for the partia l derivative m atrix and the
μ
seeking were realized by using the master
-
slav e pr ogra m mi ng m ode l
under the message
-
passing PVM parallel plat form
.
The experiment result indicated that the all
-
aroun d parall e l eff i cie ncy w a s
ideal and the inversion speed increased obviously
.
Key words
:
MT Occa m inversion
;
parallel c om putation
;
CO W
收稿日期 :2007-08-22 .
基金项目 :广西自然科学基金(桂科基 0339070)和广西教育厅科研项目(桂教科研 2003-22) .
作者简介:刘 羽(1961-),男,博士,副教授 .
E
-
mail
:
lewis
-
5709 @ 163 .
com
MT Occam
反演采用了特有的自适应算法求取拉格朗日乘子(μ
-1
) ,方法稳定 ,且无需更多的经验 ,但
速度很慢。文献[1]对方法的 μ值求取策略提出了一些改进措施 ,提高了计算速度 ,但在一定程度上改变了
Occam
反演机理 ,且方案中的初始 μ值和变化步长λ需凭经验确定 ,合理性还值得探讨
[2]
。提早结束极为
耗时的模型光滑过程(反演的后段)也可以缩短计算时间 ,但此时反演接近收敛而实际并未收敛。
Occam
反
演文档指出 ,如果要一味获得快的反演速度 ,不宜采用
Occam
反演方法 ,一些试图通过改变方法本身和收敛
条件来提高速度的做法并不可取 ,因为那样得不到最光滑的模型 ,与
Occam
反演思想相悖。
通过并行计算可以在大幅提高反演速度的同时保持
MT Occam
反演的特色。偏导数的计算和 μ值的
求取具有良好的并行特性 ,是实现并行反演的有利条件。采用大粒度的并行方式 ,可以很容易将计算任务分
解并分发到计算节点完成 ,编程相对简便 ,也容易在小型机群上实现。给出
MT Occam
反演的并行实现方
案,首先讨论并行方案的设计过程 ,然后进行性能分析 ,最后给出一个计算实例。