Parallel implementation of a velocity-stress staggered-grid finite-difference method for 2-D poroelastic wave propagation


Sheen D., Tuncay K., Baag C., Ortoleva P. J.

COMPUTERS & GEOSCIENCES, cilt.32, sa.8, ss.1182-1191, 2006 (SCI-Expanded) identifier identifier

  • Yayın Türü: Makale / Tam Makale
  • Cilt numarası: 32 Sayı: 8
  • Basım Tarihi: 2006
  • Doi Numarası: 10.1016/j.cageo.2005.11.001
  • Dergi Adı: COMPUTERS & GEOSCIENCES
  • Derginin Tarandığı İndeksler: Science Citation Index Expanded (SCI-EXPANDED), Scopus
  • Sayfa Sayıları: ss.1182-1191
  • Anahtar Kelimeler: Biot's theory, poroelastic media, wave propagation, finite difference, message passing interface (MPI), PERFECTLY MATCHED LAYER, BIOT SLOW-WAVE, NUMERICAL-SIMULATION, FLUID-FLOW, ATTENUATION, EQUATIONS, MOTION, ROCKS, MEDIA
  • Orta Doğu Teknik Üniversitesi Adresli: Hayır

Özet

Numerical simulation of wave propagation in poroelastic media demands significantly more computational capability compared to elastic media simulation. Use of serial codes in a single scientific workstation limits the size of problem. To overcome this difficulty, a parallel velocity-stress staggered-grid finite-difference method is developed for efficient simulation of wave propagation in 2-D poroelastic media. The finite difference formulation of Biot's theory has the properties of fourth order accuracy in space and second order accuracy in time. The model is decomposed into small subdomains for each processor. After each processor updates wavefields within its domain, the processors exchange the wavefields via message passing interface (MPI). The parallel implementation reduces the computational time and also allows one to study larger problems. From our numerical experiment, consistent with other I-D experiments, it is found that the presence of heterogeneity of porous medium can produce significant P-wave attenuation in the seismic frequency range. (c) 2005 Elsevier Ltd. All rights reserved.