地震层析成像LSQR算法的并行化.docVIP

第４９卷第２期地球物理学报Ｖ０１．４９，Ｎｏ．２２００６年３月ＣＨＩＮＥＳＥＪＯＵＲＮＡＬＯＦＧＥＯＰＨＹＳＩＣＳＭａｒ．，２００６刘劲松，刘福田，刘俊等．地震层析成像ＬＳＱＲ算法的并行化．地球物理学报，２００６，４９（２）：５４０～５４５ ｕｕＪＳ，ＬｉｕＦＴ，ｕｕＪ，ｅｔａ１．ＰａｒａｌｌｅｌＬＳＱＲａｌｇｏｒｉｔｈｍｓｕｓｅｄｉｎｓｅｉｓｍｉｃｔｏｍｏｇｒａｐｈｙＣｈｉｎｅｓｅＪ．Ｃ．ｅｏｐ蛳．（ｉｎＣｈｉｎｅｓｅ），２００６，４９（２）５４０—５４５ 地震层析成像ＬＳＱＲ算法的并行化 刘劲松１，刘福田１，刘俊２，郝天珧１ １中国科学院地质与地球物理研究所，北京１０００２９ ２中国科学院计算机网络信息中心超级计算中心，北京１０００８０ 摘要讨论了地震层析成像的ＬＳＱＲ算法（最ｄｘ－－乘ＱＲ分解）．在建立偏导数矩阵方程组时，对区内地震在方程 中保留震源项，引入正交投影算子进行参数分离，对区外远震采用传统的平滑处理方式，用ＬＳＱＲ法求解联立的方 程组．由于区内地震的正交分解处理和区外远震的平滑处理，使得偏导数矩阵中的非零元素成倍增加，对于大型反 演问题，这些非零元素常常达到几十ＧＢ到几百ＧＢ的数量级，巨量的内存占用成为ＬＳＱＲ算法的瓶颈．针对这一问 题，本文研究了偏导数矩阵中非零元素的分布规律，设计出合理的存储结构，采用分布式存储进行矩阵计算，提出 了ＬＳＱＲ算法的并行化方案，并