Finite-difference seismic wave modeling including surface topography

Simple item page
dc.contributor.advisorTalwani, Maniken_US
dc.creatorHestholm, Stig Ottaren_US
dc.date.accessioned2009-06-04T08:02:43Zen_US
dc.date.available2009-06-04T08:02:43Zen_US
dc.date.issued1999en_US
dc.description.abstractI present synthetics of seismic wave propagation near free surface topography. The velocity-stress formulations of both the full elastic and viscoelastic wave equations are used, and I have derived exact boundary conditions for any arbitrary, smooth topography in terms of the particle velocities. Program codes are developed for 2 and 3 dimensions (2-D and 3-D) using finite-difference (F-D) methods for both spatial and temporal numerical discretizations. An 8th order F-D method is used inside the physical model space, and the spatial F-D order decreases gradually towards the free surface topography. The discretization of the medium equations along the side and bottom boundaries, the free surface topography boundary conditions, and the forward time stepping, are all by 2nd order F-D methods. The leap-frog technique is used for time stepping everywhere except for the memory variable equations in the viscoelastic cases, where an explicit version of the unconditionally stable Crank-Nicholson method is used. I show synthetics applying the schemes to isotropic 2-D and 3-D media covered by topographies that are either described by analytic expressions or by real elevation data. These data are taken from an area in South-Eastern Norway that contains the NORESS seismic receiver array. Domains up to 60 x 60 kilometers are used in 3-D simulations, and the applied sources are plane waves generated by a plane of Ricker type point sources. These sources represent earthquakes or teleseismic explosions. For 2-D simulations I have used both plane waves and point sources, since the larger models permissible in 2-D allow for point sources to represent earthquakes or teleseismic explosions quite well. For 2-D simulations I have also included examples using layered media with randomization by a 2-D von Karman function with and without apparent anisotropy. Synthetic snapshots and seismograms show Rayleigh (Rg)-waves emanating from areas of prominent topography as well as strong surface wave directivity from some topographic features. Full viscoelastic modeling with relatively low Q-values, describing near-surface sedimentary layers, exhibit intrinsic attenuation and physical dispersion of the wavefield. Results coincide with numerous observations. 3-D simulations are performed using domain decomposition parallelization implemented by Message Passing Interface (MPI).en_US
dc.format.extent192 p.en_US
dc.format.mimetypeapplication/pdfen_US
dc.identifier.callnoTHESIS GEOL. 1999 HESTHOLMen_US
dc.identifier.citationHestholm, Stig Ottar. "Finite-difference seismic wave modeling including surface topography." (1999) Diss., Rice University. <a href="https://hdl.handle.net/1911/19392">https://hdl.handle.net/1911/19392</a>.en_US
dc.identifier.urihttps://hdl.handle.net/1911/19392en_US
dc.language.isoengen_US
dc.rightsCopyright is held by the author, unless otherwise indicated. Permission to reuse, publish, or reproduce the work beyond the bounds of fair use or other exemptions to copyright law must be obtained from the copyright holder.en_US
dc.subjectGeophysicsen_US
dc.titleFinite-difference seismic wave modeling including surface topographyen_US
dc.typeThesisen_US
dc.type.materialTexten_US
thesis.degree.departmentEarth Scienceen_US
thesis.degree.disciplineNatural Sciencesen_US
thesis.degree.grantorRice Universityen_US
thesis.degree.levelDoctoralen_US
thesis.degree.nameDoctor of Philosophyen_US
Files
Original bundle
Now showing 1 - 1 of 1
Thumbnail Image
Name:
9928545.PDF
Size:
9.43 MB
Format:
Adobe Portable Document Format
Download
Collections
Rice University Theses and Dissertations