Quasi-static finite element modeling of seismic attenuation and dispersion due to wave-induced fluid flow in poroelastic media

[1] The finite element method is used to solve Biot's equations of consolidation in the displacement-pressure (u − p) formulation. We compute one-dimensional (1-D) and two-dimensional (2-D) numerical quasi-static creep tests with poroelastic media exhibiting mesoscopic-scale heterogeneities to calculate the complex and frequency-dependent P wave moduli from the modeled stress-strain relations. The P wave modulus is used to calculate the frequency-dependent attenuation (i.e., inverse of quality factor) and phase velocity of the medium. Attenuation and velocity dispersion are due to fluid flow induced by pressure differences between regions of different compressibilities, e.g., regions (or patches) saturated with different fluids (i.e., so-called patchy saturation). Comparison of our numerical results with analytical solutions demonstrates the accuracy and stability of the algorithm for a wide range of frequencies (six orders of magnitude). The algorithm employs variable time stepping and an unstructured mesh which make it efficient and accurate for 2-D simulations in media with heterogeneities of arbitrary geometries (e.g., curved shapes). We further numerically calculate the quality factor and phase velocity for 1-D layered patchy saturated porous media exhibiting random distributions of patch sizes. We show that the numerical results for the random distributions can be approximated using a volume average of White's analytical solution and the proposed averaging method is, therefore, suitable for a fast and transparent prediction of both quality factor and phase velocity. Application of our results to frequency-dependent reflection coefficients of hydrocarbon reservoirs indicates that attenuation due to wave-induced flow can increase the reflection coefficient at low frequencies, as is observed at some reservoirs.

[1]  D. L. Johnson,et al.  Theory of frequency dependent acoustics in patchy-saturated porous media , 2001 .

[2]  M. Biot General Theory of Three‐Dimensional Consolidation , 1941 .

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

[4]  Mark Chapman,et al.  The influence of fluid-sensitive dispersion and attenuation on AVO analysis , 2006 .

[5]  B. Gurevich,et al.  Comparative review of theoretical models for elastic wave attenuation and dispersion in partially saturated rocks , 2006 .

[6]  Liu Wing-Kam Introduction to applied mathematics: Gilbert Strang (Wellesley-Cambridge press, Wellesley, MA, 1986) 758 pp , 1987 .

[7]  J. White,et al.  Computed seismic speeds and attenuation in rocks with partial gas saturation , 1975 .

[8]  M. Rapoport,et al.  Direct detection of oil and gas fields based on seismic inelasticity effect , 2004 .

[9]  Jonathan Richard Shewchuk,et al.  Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator , 1996, WACG.

[10]  Marcel Frehner,et al.  Finite-element simulations of Stoneley guided-wave reflection and scattering at the tips of fluid-filled fractures , 2010 .

[11]  T. Mukerji,et al.  The Rock Physics Handbook , 1998 .

[12]  Jonathan Richard Shewchuk,et al.  Delaunay refinement algorithms for triangular mesh generation , 2002, Comput. Geom..

[13]  Douglas R. Schmitt,et al.  Amplitude and AVO responses of a single thin bed , 2003 .

[14]  U. N. Baton “Computational geomechanics-with special reference to earthquake engineering” by O. C. Zienkiewicz, A. H. C. Chan, M. Paston, B. A. Schrefler and T. Shiomi , 2000 .

[15]  José M. Carcione,et al.  P-wave seismic attenuation by slow-wave diffusion: Effects of inhomogeneous rock properties , 2006 .

[16]  C. Zener Elasticity and anelasticity of metals , 1948 .

[17]  Juan E. Santos,et al.  Equivalent viscoelastic solids for heterogeneous fluid-saturated porous rocks , 2009 .

[18]  Y. Masson,et al.  Acoustic attenuation in self-affine porous structures. , 2006, Physical review letters.

[19]  J. Carcione,et al.  Seismic Low-frequency Anomalies In Multiple Reflections From Thinly-layered Poroelastic Reservoirs , 2007 .

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

[21]  Boris Gurevich,et al.  Attenuation and dispersion of P-waves in porous rocks with planar fractures: Comparison of theory and numerical simulations , 2006 .

[22]  Yder J. Masson,et al.  Poroelastic finite difference modeling of seismic attenuation and dispersion due to mesoscopic-scale heterogeneity , 2007 .

[23]  Yder J. Masson,et al.  Finite difference modeling of Biot's poroelastic equations at seismic frequencies , 2006, Journal of Geophysical Research.

[24]  L. C. Wood,et al.  The limits of resolution of zero-phase wavelets , 1982 .

[25]  B. Gurevich,et al.  P‐wave dispersion and attenuation in fractured and porous reservoirs – poroelasticity approach , 2009 .

[26]  Jonathan Richard Shewchuk,et al.  Triangle : Engineering a 2 D Quality Mesh Generator and Delaunay , 2005 .

[27]  J. Altmann,et al.  Anisotropic dispersion and attenuation due to wave‐induced fluid flow: Quasi‐static finite element modeling in poroelastic solids , 2010 .

[28]  D. L. Johnson Erratum: “Theory of frequency dependent acoustics in patchy-saturated porous media” [J. Acoust. Soc. Am.110(2), 682–694 (2001)] , 2007 .

[29]  José M. Carcione,et al.  SOME ASPECTS OF THE PHYSICS AND NUMERICAL MODELING OF BIOT COMPRESSIONAL WAVES , 1995 .

[30]  Thomas M. Daley,et al.  Seismic low-frequency effects in monitoring fluid-saturated reservoirs , 2004 .

[31]  J. Z. Zhu,et al.  The finite element method , 1977 .

[32]  G. Goloshubin,et al.  Reservoir imaging using low frequencies of seismic reflections , 2006 .

[33]  Gilbert Strang,et al.  Introduction to applied mathematics , 1988 .

[34]  R. Knight,et al.  Elastic wave velocities during evaporative drying , 1998 .

[35]  Tobias Malte Müller,et al.  Seismic attenuation due to wave‐induced flow: Why Q in random structures scales differently , 2006 .

[36]  Y. Podladchikov,et al.  “Low-frequency reflections from a thin layer with high attenuation caused by interlayer flow,” GEOPHYSICS, 74, no. 1, N15–N23 , 2009 .

[37]  Andrew N. Norris,et al.  Low‐frequency dispersion and attenuation in partially saturated rocks , 1993 .

[38]  Boris Gurevich,et al.  Wave-induced fluid flow in random porous media: attenuation and dispersion of elastic waves. , 2005, The Journal of the Acoustical Society of America.

[39]  James G. Berryman,et al.  Seismic attenuation due to wave-induced flow , 2004 .

[40]  Roderic S. Lakes,et al.  Viscoelastic measurement techniques , 2004 .

[41]  O. Zienkiewicz,et al.  Dynamic behaviour of saturated porous media; The generalized Biot formulation and its numerical solution , 1984 .

[42]  Yuri Y. Podladchikov,et al.  Low-frequency reflections from a thin layer with high attenuation caused by interlayer flow , 2009 .

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

[44]  J. Harris,et al.  Exact expression for the effective acoustics of patchy-saturated rocks , 2010 .

[45]  J. White,et al.  Low‐frequency seismic waves in fluid‐saturated layered rocks , 1975 .

[46]  N. Dutta,et al.  Attenuation and dispersion of compressional waves in fluid-filled porous rocks with partial gas saturation (White model); Part I, Biot theory , 1979 .

[47]  Roger A. Clark,et al.  Estimation of Q from surface seismic reflection data , 1998 .

[48]  José M. Carcione,et al.  White's model for wave propagation in partially saturated rocks: Comparison with poroelastic numerical experiments , 2003 .

[49]  Holger Steeb,et al.  Comparison of finite difference and finite element methods for simulating two-dimensional scattering of elastic waves , 2008 .

[50]  Tobias M. Müller,et al.  Finite-difference modeling of wave propagation and diffusion in poroelastic media , 2009 .

[51]  M. Lebedev,et al.  Seismic wave attenuation and dispersion resulting from wave-induced flow in porous rocks — A review , 2010 .

[52]  Tadahiko Shiomi,et al.  Practical Programming in Computational Geomechanics: With Special Reference to Earthquake Engineering , 1999 .

[53]  Boris Gurevich,et al.  Velocity and attenuation of elastic waves in finely layered porous rocks , 1995 .