| 地震层析成像LSQR算法的并行化 |
| |
| 引用本文： | 刘劲松,刘福田,刘俊,郝天珧.地震层析成像LSQR算法的并行化[J].地球物理学报,2006,49(2):540-545. |
| |
| 作者姓名： | 刘劲松 刘福田 刘俊 郝天珧 |
| |
| 作者单位： | 1.中国科学院地质与地球物理研究所, 北京 100029 2 中国科学院计算机网络信息中心超级计算中心，北京 100080 |
| |
| 基金项目： | 中国科学院知识创新工程信息化建设重大专项（INF105-SCE）和国家自然科学基金（40034010）资助. |
| |
| 摘 要： | 讨论了地震层析成像的LSQR算法(最小二乘QR分解). 在建立偏导数矩阵方程组时，对区内地震在方程中保留震源项，引入正交投影算子进行参数分离，对区外远震采用传统的平滑处理方式，用LSQR法求解联立的方程组. 由于区内地震的正交分解处理和区外远震的平滑处理，使得偏导数矩阵中的非零元素成倍增加，对于大型反演问题，这些非零元素常常达到几十GB到几百GB的数量级，巨量的内存占用成为LSQR算法的瓶颈. 针对这一问题，本文研究了偏导数矩阵中非零元素的分布规律，设计出合理的存储结构，采用分布式存储进行矩阵计算，提出了LSQR算法的并行化方案，并在联想深腾6800超级计算机上实现. 导出了LSQR算法的并行效率估算公式. 对两个地区的实际地震层析成像数据进行了效率测试.
|
| 关 键 词： | 地震层析成像 走时反演 LSQR算法 并行计算 |
| 文章编号： | 0001-5733（2006）02-0540-06 |
| 收稿时间： | 2005-04-06 |
| 修稿时间： | 2005-04-062005-12-26 |
Parallel LSQR algorithms used in seismic tomography |
| |
| LIU Jin-Song,LIU Fu-Tian,LIU Jun,HAO Tian-Yao.Parallel LSQR algorithms used in seismic tomography[J].Chinese Journal of Geophysics,2006,49(2):540-545. |
| |
| Authors: | LIU Jin-Song LIU Fu-Tian LIU Jun HAO Tian-Yao |
| |
| Institution: | 1.Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China 2 Supercomputing Center of Computer Network Information Center, Chinese Academy of Sciences, Beijing 100080, China |
| |
| Abstract: | We discuss the LSQR algorithms used in earthquake travel time tomography. We keep the epicenter terms in the equation for regional events, and then use the orthogonal projection method to eliminate the epicenter terms. For tele-events, the classic smoothing process is used. The number of non-zero elements in the partial derivative matrix is increased by several times because of the orthogonal projection and smoothing processes. For a large scale inversion problem, the amount of non-zero elements can be dozens of Gigabytes or hundreds of Gigabytes. The huge amount of memory requirement becomes the bottle neck of LSQR algorithms. matrix, designed an efficient data structure for the sparse matrix, used a distributed memory and computation scheme for matrix computation, and implemented it on a multi-processor super-computer. We have derived an estimation formula of parallel efficiency and tested two real tomography models. |
| |
| Keywords: | Seismic tomography Traveltime inversion LSQR algorithms Parallel computation |
| 本文献已被 CNKI 维普 万方数据 等数据库收录！ |
| 点击此处可从《地球物理学报》浏览原始摘要信息 |
| 点击此处可从《地球物理学报》下载免费的PDF全文 |
