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

Abstract 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 1-D experiments, it is found that the presence of heterogeneity of porous medium can produce significant P-wave attenuation in the seismic frequency range.

[1]  Q. H. Liu,et al.  A staggered-grid finite-difference method with perfectly matched layers for poroelastic wave equations. , 2001, The Journal of the Acoustical Society of America.

[2]  M. Biot.  MECHANICS OF DEFORMATION AND ACOUSTIC PROPAGATION IN POROUS MEDIA , 1962 .

[3]  M. Biot,et al.  THE ELASTIC COEFFICIENTS OF THE THEORY OF CONSOLIDATION , 1957 .

[4]  P. Ortoleva,et al.  Wave propagation in poroelastic media: A velocity‐stress staggered‐grid finite‐difference method with perfectly matched layers , 2003 .

[5]  J. Geertsma,et al.  SOME ASPECTS OF ELASTIC WAVE PROPAGATION IN FLUID-SATURATED POROUS SOLIDS , 1961 .

[6]  M. Biot.  Theory of Propagation of Elastic Waves in a Fluid‐Saturated Porous Solid. I. Low‐Frequency Range , 1956 .

[7]  S. Shapiro,et al.  Seismic signatures of permeability in heterogeneous porous media , 1999 .

[8]  G. McMechan,et al.  Numerical simulation of seismic responses of poroelastic reservoirs using Biot theory , 1991 .

[9]  T. Plona,et al.  Observation of a second bulk compressional wave in a porous medium at ultrasonic frequencies , 1980 .

[10]  J. Berryman,et al.  Biot slow‐wave effects in stratified rock , 2002 .

[11]  Amos Nur,et al.  High-resolution shallow-seismic experiments in sand, Part I: Water table, fluid flow, and saturation , 1998 .

[12]  Guido Kneib,et al.  Accurate and efficient seismic modeling in random media , 1993 .

[13]  Qing Huo Liu,et al.  The application of the perfectly matched layer in numerical modeling of wave propagation in poroelastic media , 2001 .

[14]  J. Sochacki.  Absorbing boundary conditions for the elastic wave equations , 1988 .

[15]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[16]  B. Engquist,et al.  Absorbing boundary conditions for acoustic and elastic wave equations , 1977, Bulletin of the Seismological Society of America.

[17]  Qing Huo Liu,et al.  PERFECTLY MATCHED LAYERS FOR ELASTODYNAMICS: A NEW ABSORBING BOUNDARY CONDITION , 1996 .

[18]  Robert W. Graves,et al.  Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences , 1996, Bulletin of the Seismological Society of America.

[19]  Nanxun Dai,et al.  Wave propagation in heterogeneous, porous media: A velocity‐stress, finite‐difference method , 1995 .

[20]  J. Kristek,et al.  3D Heterogeneous Staggered-grid Finite-difference Modeling of Seismic Motion with Volume Harmonic and Arithmetic Averaging of Elastic Moduli and Densities , 2002 .

[21]  Wei-Ping Huang,et al.  Application and optimization of PML ABC for the 3-D wave equation in the time domain , 2003 .

[22]  Weng Cho Chew,et al.  A 3D perfectly matched medium from modified maxwell's equations with stretched coordinates , 1994 .

[23]  Full frequency-range transient solution for compressional waves in a fluid-saturated viscoacoustic porous medium , 1996 .

[24]  F. Morgan,et al.  Deriving the equations of motion for porous isotropic media , 1992 .

[25]  J. Carcione,et al.  Velocity and attenuation in partially saturated rocks: poroelastic numerical experiments , 2003 .

[26]  A. Levander.  Fourth-order finite-difference P-SV seismograms , 1988 .

[27]  N. Dutta,et al.  Seismic reflections from a gas‐water contact , 1983 .

[28]  Børge Arntsen,et al.  Numerical simulation of the Biot slow wave in water‐saturated Nivelsteiner Sandstone , 2001 .

[29]  B. Gurevich,et al.  Seismic attenuation in finely layered porous rocks : Effects of fluid flow and scattering , 1997 .

[30]  Jian-Ming Jin,et al.  Perfectly Matched Layers in the Discretized Space: An Analysis and Optimization , 1996 .