3D frequency-domain finite-difference viscoelastic-wave modeling using weighted average 27-point operators with optimal coefficients

被引:38
作者
Gosselin-Cliche, Benjamin [1 ]
Giroux, Bernard [2 ]
机构
[1] Coll Francois Xavier Garneau, Quebec City, PQ, Canada
[2] Ctr Eau Terre Environm, Inst Natl Rech Sci, Quebec City, PQ, Canada
关键词
PERFECTLY MATCHED LAYER; ABSORBING BOUNDARY-CONDITIONS; PARALLEL DIRECT SOLVER; GRAZING-INCIDENCE; FORM INVERSION; LINEAR-SYSTEMS; ELASTIC-WAVES; PROPAGATION; EQUATIONS; MEDIA;
D O I
10.1190/GEO2013-0368.1
中图分类号
P3 [地球物理学]; P59 [地球化学];
学科分类号
0708 ; 070902 ;
摘要
Experimental data suggest that the viscoelastic behavior of rocks is more easily and accurately described in the frequency domain than in the time domain, supporting the idea of simulating seismic wave propagation in the frequency domain. We evaluated weighted-averaged 27-point finite-difference operators for 3D viscoelastic wave modeling in the frequency domain. Within the proposed framework, we developed general equations for normalized phase velocities that can be used with arbitrary finite-difference operators. Three sets of weighting coefficients for second-order central finite-difference operators that minimize the numerical dispersion for up to five grid points per wavelength were found using a damped least-squares (LS) criterion as well as a global optimization scheme based on l(1)- and l(2)-norm criteria. The three sets produced very similar dispersion curves, and improvement provided by global optimization appeared marginal in this respect. We also evaluated a discrete form for the heterogeneous formulation of the 3D viscoelastic equations with a perfectly match layer (PML). Heuristic performance assessment of frequency-dependent PML absorption coefficients provided a simple rule giving good results for eight PMLs at all frequencies. The proposed formalism was implemented with a massively parallel direct solver. Modeling results were compared with an analytic solution and a time-domain finite-difference code, and they gave good agreement when using LS and l(2)-norm optimal coefficients. On the other hand, l(1)-norm coefficients produced noisy results, indicating that minimizing the difference between analytic and numerical phase velocities, although necessary, is not a sufficient condition to guarantee low-numerical noise. Finally, analysis of the computational resources required to factorize the impedance matrix revealed that the memory complexity of the factorization is O(292N)(4) for an N-3 grid, compared to O(30N(4)) for the viscoacoustic case.
引用
收藏
页码:T169 / T188
页数:20
相关论文
共 55 条
[1]   Seismic wave attenuation in carbonates [J].
Adam, L. ;
Batzle, M. ;
Lewallen, K. T. ;
van Wijk, K. .
JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH, 2009, 114
[2]  
Ajo-Franklin J. B., 2005, TECHNICAL REPORT
[3]   Hybrid scheduling for the parallel solution of linear systems [J].
Amestoy, PR ;
Guermouche, A ;
L'Excellent, JY ;
Pralet, S .
PARALLEL COMPUTING, 2006, 32 (02) :136-156
[4]   A fully asynchronous multifrontal solver using distributed dynamic scheduling [J].
Amestoy, PR ;
Duff, IS ;
L'Excellent, JY ;
Koster, J .
SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS, 2001, 23 (01) :15-41
[5]  
[Anonymous], 1998, A Practical Guide to Pseudospectral Methods
[6]  
[Anonymous], 2005, QUANTITATIVE SEISMIC, DOI [10.1017/CBO9780511600074, DOI 10.1017/CBO9780511600074, DOI 10.1017/CBO9780511600074.001]
[7]  
[Anonymous], TECHNICAL REPORT
[8]   Fluid mobility and frequency-dependent seismic velocity - Direct measurements [J].
Batzle, ML ;
Han, DH ;
Hofmann, R .
GEOPHYSICS, 2006, 71 (01) :N1-N9
[9]   A PERFECTLY MATCHED LAYER FOR THE ABSORPTION OF ELECTROMAGNETIC-WAVES [J].
BERENGER, JP .
JOURNAL OF COMPUTATIONAL PHYSICS, 1994, 114 (02) :185-200
[10]   Parallel 3-D viscoelastic finite difference seismic modelling [J].
Bohlen, T .
COMPUTERS & GEOSCIENCES, 2002, 28 (08) :887-899