1887
Volume 37, Issue 2
  • ISSN: 0812-3985
  • E-ISSN: 1834-7533

Abstract

Standard and rotated staggered grid algorithms, to eighth order in space and second order in time, have been developed for simulating seismic wave propagation in heterogeneous viscoelastic/anisotropic media. These algorithms have been validated by comparing their numerical results with analytical solutions for homogeneous models. Numerical simulation results show that high-order staggered finite-difference methods are better than staggered pseudo-spectral methods in modelling strong heterogeneous viscoelastic/anisotropic media. For strong velocity or density contrast interfaces, such as fluid-solid interfaces, staggered pseudo-spectral methods generate apparent artefacts although they give minimum spatial dispersion; while high-order staggered finite-difference methods, whether rotated-staggered or standard-staggered, give reasonable modelling results. This is especially true for fracture modelling, if the fracture is simulated with zero compressional and shear velocities and a very low density such as less than 10-8 kg.m-3.

Based on the proposed algorithms, a uniform random function is used to establish random models with two kinds of components, i.e., shales or clays and pure sandstones. These two media consist of various synthetic shale sand formations to simulate random formations with various shale contents. Elastic wave velocities and attenuations for various shale contents are simulated properly at a frequency range around 10 kHz to 1 MHz. The numerical results show that our algorithms can be used for modelling composite media to study elastic wave attenuation and velocity variation for various rocks, including fracture modelling. The modelling analyses demonstrate that velocity models with shale contents from 30% to 60% generate strong back wave scattering for a wave length around 10 times the minimum size of the synthetic shale lumps. The velocity variations with respect to the shale contents evaluated from synthetic waveforms are in good agreement with the results predicted by Wyllie’s time-average equation at high frequencies.

© 2006 ASEG
Loading

Article metrics loading...

/content/journals/10.1071/EG06160
2006-06-01
2026-01-15

  • From This Site

    /content/journals/10.1071/EG06160
    dcterms_title,dcterms_subject,pub_keyword
    -contentType:Journal -contentType:Contributor -contentType:Concept -contentType:Institution
    10
    5
Loading full text...

Full text loading...

References

  1. Blanch. J.O., Robertsson. J.O.A., and Symes. W.W., 1995. Modelling of a constant Q\ Methodology and algorithm for an efficient and optimally inexpensive viscoelastic technique: Geophysics.60, 176-184.
  2. Carcione. J.M., Kosloff. D., Behle. A., and Seriani. G., 1992. A spectral Scheme for wave propagation simulation in 3-D elastic-anisotropic media: Geophysics.57, 1593-4607.
  3. Carcione. J.M., 1995. Constitutive model and wave equation for linear, viscoelastic, anisotropic media: Geophysics.60, 537-548.
  4. Carcione. J.M., and Wang. J.P., 1993, A Chcbyshcv collocation method for the clastodynamic equation in generalized coordinates: Computational Fluid Dynamics Journal.2, 269-291).
  5. Carcione. J.M., Herman. G., and ten Kroode. F.P.E., 2002. Seismic modelling: Geophysics67, 1304-1325.
  6. Dai. N., Vafidis. A., and Kanasewich. E.R., 1995. Wave propagation in heterogeneous, porous media: A velocity-stress. finite-difference method: Geophysics.60, 327-340.
  7. Dauksher, W., and Emery. A.F., 1997. Accuracy in modelling the acoustic wave equation with Chebyshev spectral finite elements: Finite Elements in Analysis and Design.26, 115-128.
  8. Gazdag, J., 1981. Modelling of the acoustic wave equation with transform methods: Geophysics, 46, 854-859.
  9. Graves. R., 1996. Simulation seismic wave propagation in 3D elastic media using staggered-grid finite differences: Bulletin of the Seismologicul Society of America.86, 1091-1106.
  10. Hestholm. S.O., and Ruud. B.O., 2000, 2D finite-difference viscoelastic wave modelling including surface topography: Geophysical Prospecting48, 341-373.
  11. Holbcrg, O., 1987. Computational aspects of the choice of operator and sampling interval for numerical differentiation in large-scale simulation of wave equation: Geophysical Prospecting, 35, 629-655.
  12. Igel. H., Mora. P., and Riollet, B., 1995. Anisotropic wave propagation through finite- difference grids: Geophysics, 60, 1203-1216.
  13. Kindlean, M., Kamel. A., and Sguazzero, P., 1990. On the construction and efficiency of staggered numerical differentiators for the wave equation: Geophysics, 55, 107-110.
  14. Komatitsch. D., and Vilotte J.P., 1998. The spectral element method: an efficient tool to simulate the seismic response of 2D and 3D geological structures: Bulletin of the Seismologicul Society of America, 88, 368-392.
  15. Komatitsch. D., Vilotte, J.P., Vai. R., Castillo-Covarrubias. J.M., and Sanchez-Sesma. F.J., 1999. The spectral clement method for elastic wave equations -application to 2-D and 3-D seismic problems: International Journal of Numerical Methods in Engineering45, 1139-1164.
  16. Kosloff. D., Kessler. D., Filho. A.Q., Tessmer, E., Bchlc. A., and Strahilevitz,, R., 1990. Solution of the equations of dynamic elasticity by a Chebychev spectral method: Geophysics, 55, 734-748.
  17. Kosloff. R., and Baysal. E., 1982. Forward modeling by a Fourier method: Geophysics, 47, 1402-1412.
  18. Levander, A.R., 1988, Fourth-order finite-difference P-SV seismograms: Geophysics, 53, 1425-1436.
  19. Madariaga, R., 1976, Dynamics of an expanding circular fault: Bulletin of the Seismological Society of America, 65, 163-182.
  20. Marfurt. K.J., 1984, Accuracy of finite-difference and finite-element modelling of the scalar and elastic wave equations: Geophysics. 49. 533-549.
  21. Mu. Y.G., 1984. Elastic wave migration with finite element method: Acta Geophysica Sinica, 11, 268-278.
  22. Orszag. S.A., 1980, Spectral methods for problems in complex geometries: Journal of Computational Physics.37, 70-92.
  23. Ozdenvar.T., and McMechan, G. A., 1996, Causes and reduction of numerical artefacts in pseudospcctral wavefield exploration: Geophysical Prospecting, 126, 819-828.
  24. Patera. A.T., 1984. A spectral element method for fluid dynamics: laminar flow in a channel expansion: Journal of Computational Physics, 54, 468-488.
  25. Priolo. E., Carcione. J., and Seriani. G., 1993, Numerical simulation of interface waves by highorder spectral modelling techniques: Journal of the Acoustic Society of America, 95, 681-693.
  26. Robertsson. J.O.A., 1996. A numerical free-surface condition for elastic/viscoelastic finite-difference modeling in the presence of topography: Geophysics.61, 1921-1934.
  27. Saenger. E.H., Gold. N., and Shapiro. S. A., 2000. Modeling the propagation of clastic waves using a modified finite-difference grid: Wave Motion, 31, 77-92.
  28. Seriani. G., Priolo. E., Carcione. J.M., and Padovani. E., 1992. High-order spectral element method for clastic wave modelling: 62nd Annual International Meeting. Society of Exploration Geophysics. Expanded Abstracts, 1285-1288.
  29. Seriani. G., and Priolo. E., 1994. Spectral element method for acoustic wave simulation in heterogeneous media: Finite Element in Analysis ami Design.16, 337-348.
  30. Seriani. G., 1997. A parallel spectral element method for acoustic wave modelling: Journal of Computational Acoustics.5, 53-69.
  31. Seriani, G., 1998, 3-D large-scale wave propagation modelling by spectral element method on Cray T3E: Computational Methods in Applied Mechanics and Engineering, 164, 235-247.
  32. Tal-Ezcr, H., 1986, Spectral methods in time for hyperbolic problems: SIAM Journal of Numerical Analysis, 23, 12-26.
  33. Teng. Y.C., 1988. Three-dimensional finite element analysis of waves in an acoustic media with inclusion: Journal of the Acoustic Society of America, 86, 414-422.
  34. Tessmer, E., Kosloff. D., and Bohle. A., 1992, Elastic wave propagation simulation in the presence of surface topography: Geophysical Journal International.108, 621-632.
  35. Virieux, J., 1986, P-SV wave propagation in heterogeneous media: Velocity-strcss finite-difference method: Geophysics, 51, 889-901.
  36. Wang, X.M., and Dodds, K., 2003, Viscoelastic wave modelling using a staggered high-order finite-difference in inhomogeneous media: Proceedings of the 6th SEGJ international symposium, Tokyo, Japan.
  37. Wang. X.M., and Zhang. H.L., 2004. Modeling of clastic wave propagation on a curved free surface using an improved finite-difference algorithm: Science in China, Series G-Physics and Astronomy, 47, 633-648.
  38. Wang, X.M., Hornby. B., and Dodds. K., 2002. Dipole sonic response in deviated borehole penetrating an anisotropic formation: 72nd Annual International Meeting, Society of Exploration Geophysicists. Expanded Abstracts, 360-363.
/content/journals/10.1071/EG06160
Loading
  • Article Type: Research Article
Keyword(s): Finite-differences; pseudo-spectral; random velocity media; rotated staggered grids

Most Cited This Month Most Cited RSS feed

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error