1887
Volume 44 Number 4
  • E-ISSN: 1365-2478

Abstract

Abstract

Computational seismic modelling (CSM) plays an important role in the geophysical industry as an established aid to seismic interpreters. Numerical solution of the elastic wave equations has proved to be a very important tool for geophysicists in both forward modelling and migration. Among the techniques generally used in CSM, we consider the finite‐element method (FEM) and investigate its computational and visualization requirements. The CSMFEM program, designed for this purpose and developed on an IBM 3090 computer with vector facility, is described in detail. It constitutes a numerical laboratory for performing computer experiments. Two Newmark type algorithms for time integration are compared with other time integration schemes, and both direct and iterative methods for solving the corresponding large sparse system of linear algebraic equations are analysed. Several numerical experiments to simulate seismic energy propagation through heterogeneous media are performed. Synthetics in the form of common shot gathers, vertical seismic profiles and snapshots are suitably displayed, since with the large amounts of data obtained from CSM research, methods for visualization of the computed results must be developed. The FEM is compared with other numerical tools, such as finite‐difference and pseudo‐spectral methods.

Loading

Article metrics loading...

/content/journals/10.1111/j.1365-2478.1996.tb00168.x
2006-04-28
2024-04-28
Loading full text...

Full text loading...

References

  1. BambergerA., ChaventG. and LaillyP.1980. Etude de schemes numeriques pour les equations de I'elastodynamique lineaire. Institute National de Rocherche en Informatique et en Automatique.
  2. BankR.E., ShermanA.H. and WaiserA.1983. Refinement algorithms and data structures for regular local refinement. In: Scientific Computing: Applications of Mathematics and Computing to the Physical Sciences (ed. R.S.Stepleman ), pp. 3–17. North Holland, Publishing Co.
    [Google Scholar]
  3. BatheK.J. and ChaeS.W.1989. On automatic mesh construction and mesh refinement. In: Finite Element Analysis, Computer and Stuctures. 32, (3/4), 911–936.
    [Google Scholar]
  4. BatheK.J. and WilsonE.L.1976. Numerical Methods in Finite Element Analysis. Prentice‐Hall Inc.
    [Google Scholar]
  5. BelytschkoT. and MullenR.1978. On dispersive properties of finite element solutions. In: Modern Problems in Elastic Wave Propagation (eds J.Miklowitz et al.), John Wiley & Sons, Inc.
    [Google Scholar]
  6. CareyG.F. and OdenJ.T.1983. Finite Elements. A Second Course. Prentice‐Hall Inc.
    [Google Scholar]
  7. ChowY.K.1985. Accuracy of consistent and lumped viscous damper in wave propagation problems. International Journal of Numerical Methods in Engineering21, 723–732.
    [Google Scholar]
  8. CiarletP.G.1978. The Finite Element Method for Elliptic Problems. North Holland Publishing Co.
    [Google Scholar]
  9. CullenM.J.P.1982. The use of quadratic finite element methods and irregular grids in the solution of hyperbolic problems. Journal of Computational Physics45, 221–245.
    [Google Scholar]
  10. CuthillE. and MckeeJ.1969. Reducing the bandwidth of sparse symmetric matrices. Technical Note. Naval Ship Research and Development Center, Washington, D.C., Applied Mathematics Laboratory, AML‐40–69.
  11. DuffI.S.ErismanA.M. and ReidJ.K.1986. Direct Methods for Sparse Matices. Oxford University Press.
    [Google Scholar]
  12. FriedY.1980. Irregular finite element meshes in elastodynamics. International Journal of Numerical Methods in Engineering16, 626–628.
    [Google Scholar]
  13. GardnerG., GardnerL. and GregoryA.1974. Formation velocity and density, the diagnostic basis for stratigraphic taps. Geophysics39, 770–780.
    [Google Scholar]
  14. GeradinM.1974. A classification and discussion of integration operators for transient structural response. 12th AIAA Aerospace Sciences meeting, Washington D.C., Expanded Abstracts, 74–105.
  15. GordonW.J. and HallC.A.1973. Construction of curvilinear co‐ordinate systems and applications to mesh generation. International Journal of Numerical Methods in Engineering7, 461–477.
    [Google Scholar]
  16. GoudreauG.L. and TaylorR.L.1972. Evaluation of numerical integration methods in elastodynamics. Computational Methods in Applied Mechanics and Engineering2, 69–97.
    [Google Scholar]
  17. HagemanL. and YoungD.M.1981. Applied lterative Methods. Academic Press, Inc.
    [Google Scholar]
  18. HouboltJ.C.1950. A recurrence matrix solution for the dynamic response of elastic aircraft. Journal of Aeronautical Sciences17, 540–550.
    [Google Scholar]
  19. HughesT.J.R.1987. The Finite Element Method. Prentice‐Hall Inc.
    [Google Scholar]
  20. KamelA., KindelánM. and SguazzeroP.1988. Seismic computations on the IBM 3090 vector multiprocessor. IBM Systems Journal.
  21. KellyK.R., WardR.W., TreitelS. and AlfordR.M.1976. Synthetic seismograms: a finite‐difference approach. Geophysics41, 2–27.
    [Google Scholar]
  22. KincaidD.R., OppeT.C., RespersJ.R. and YoungD.M.1984. ITPACV‐2C user's guide. CNA191, Center Numerical Analysis, University of Texas, Austin, Texas.
  23. KindelánM., SguazzeroP. and KamelA.1987. Elastic modeling with Fourier methods on the IBM 3090 vector multiprocessor. In: Scientific Computing on IBM Vector Multiprocessors (eds R.Benzi and P.Sguazzero ), pp. 635–674. IBM European Center for Scientific and Engineering Computing, Rome .
    [Google Scholar]
  24. MarfurtK.J.1984. Accuracy of finite‐difference and finite‐element modelins of the scalar and elastic wave equation. Geophysics49, 533–549.
    [Google Scholar]
  25. MullenR. and BelytschkoT.1982. Dispersion analysis of finite element semidiscretizations of the two‐dimensional wave equation. International Journal of Numerical Methods in Engineering18, 11–29.
    [Google Scholar]
  26. NewmarkN.M.1959. A method of computation for structural dynamics. American Society of Civil Enginseering, Mechanical Division85, 67–94.
    [Google Scholar]
  27. PittewayM.L.V.1967. Algorithm for drawing ellipses or hyperbolae with a digital plotter. Computer Journal10, 282–289.
    [Google Scholar]
  28. QuarteroniA.1992. Mathematical aspects of domain decomposition methods. Internal Report 79/P, Dipartimento di Matematica, Politecnico di Milano, Italy.
  29. SerónF.J., SanzF.J. and KindelánM.1989. Elastic wave propagation with the finite element method. Internal Report ICE‐0028, IBM European Center for Scientific and Engineering Computing, Rome, Italy.
  30. SochackiJ., KubicheckR., GeorgeJ., FletcherW.R. and SmithsonS.1987. Absorbing boundary conditions and surface waves. Geophysics52, 60–71.
    [Google Scholar]
  31. SzaboB.A. and DunavantD.A.1980. An adaptive procedure based on the p‐version of the flnite element method. In: Finite Element workshop (ed. BabuskaV. ) Lab. Numer. Anal., Inst. Phys. Sci. Technol., University of Maryland, College park, 142–155.
    [Google Scholar]
  32. TaigI.C.1961. Structural analysis by the matrix displacement method. Report S017, English Electric Aviation.
  33. Van AkenJ.R.1984. An efficient ellipse‐drawing algorithm. Computer Graphics and Applications4, 24–35.
    [Google Scholar]
  34. WilsonE.L., FarhoomandY. and BatheK.J.1973. Nonlinear dynamic analysis of complex structures. Earthquake Engineering and Structural Dynamics1, 241–252.
    [Google Scholar]
  35. ZienkiewtczO.C.1977. The Finite Element Method. McGraw‐Hill, Book Co.
    [Google Scholar]
  36. ZlatevZ.1980. On some pivotal strategies in Gaussian elimination by sparse technique. SIAM Journal of Numerical Analysis17, 18–30.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.1111/j.1365-2478.1996.tb00168.x
Loading
  • Article Type: Research Article

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