地震层析成像LSQR算法的并行化_刘劲松.pdfVIP

49 2 　　　　　 Vol.49, No.2 2006 3 CHINESE　JOURNAL　OF　GEOPHYSICS Mar., 2006 , , 　. LSQR., 2006, 49(2):540～ 545 Liu J S, Liu F T, Liu J, et al.Parallel LSQR algorithms used in seismic tomography ChineseJ .Geophys.(in Chinese), 2006, 49(2): 540 ～ 545 LSQR 2 刘劲松, 刘福田, 刘　俊, 郝天珧 , 　 00029 2 , 　 00080 　　LSQR (QR )., , , , LSQR ., , , GB GB , LSQR . , , , , LSQR , 6800 .LSQR . . 　, , LSQR , 　000 -5733(2006)02-0540-06 　P63 　2005-04-06, 2005- 2-26 Parallel LSQR algorithms used in seismic tomography 2 LIU Jin-Song , LIU Fu-Tian , LIU Jun , HAO Tian-Yao Instituteof Geology and Geophysics, Chinese Academy of Sciences, Beijing 00029, China 2 Supercomputing Center of ComputerNetwork Information Center, Chinese Academy of Sciences, Beijing 00080, 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. To solve this problem, we have studied the distribution property of non-zero elements in the partial deriv