2026-09-22T17:54:44Zhttps://dspace.mit.edu/server/oai/requestoai:dspace.mit.edu:1721.1/93672022-01-13T07:54:24Zcom_1721.1_7582com_1721.1_7581col_1721.1_131023
M. Nafi Toksöz.
Hayashi, Koichi, 1967-
Massachusetts Institute of Technology. Dept. of Earth, Atmospheric, and Planetary Sciences.
Massachusetts Institute of Technology. Department of Earth, Atmospheric, and Planetary Sciences
2005-08-22T20:38:43Z
2005-08-22T20:38:43Z
1999
http://hdl.handle.net/1721.1/9367
44623119
Thesis (S.M.)--Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences, September 1999.
"August 6, 1999."
Includes bibliographical references (leaves 188-190).
We have developed a two-dimensional viscoelastic finite-difference modeling method for highly complex surface topography and subsurface structures. Realistic modeling of seismic wave propagation in the near surface region is complicated by many factors, such as strong heterogeneity, topographic relief and large attenuation. In order to account for these complications, we use a velocity-stress staggered grid and employ an 0(2,4) accurate viscoelastic finite-difference scheme. The implementation includes an irregular free surface condition for topographic relief and a variable grid technique in the shallow parts of the model. Several methods of free surface condition are bench marked, and an accurate and simple condition is proposed. In the proposed free surface condition, stresses are calculated so that the shear and normal stresses perpendicular to the boundary are zero. The calculation of particle velocities does not involve any specific calculations, and the particle velocities are set to zero above the free surface. A stable variable grid method is introduced, where we use a three times finer grid in the near surface or low velocity region compared to the rest of the model. In order to reduce instability, we apply averaging or weighting to the replacement of the coarse grid components within the fine grid. The method allows us to avoid any limitation of the shape of the grid size boundary. Numerical tests indicate that approximately ten grid-points per shortest wavelength with the variable grid method results in accurate calculations. The method requires a stair-shaped discretization of a free surface. We investigated the stair-shaped structures, and found that the cause of the dispersion from irregular free surface is mainly a numerical error due to the large grid sizes rather than the Rayleigh waves scattering due to the stair-shaped boundary. The finite-difference modeling is applied to the investigation of near surface wave propagation. Several numerical simulations are performed to show the characters of wave propagation in the near surface region. The simulations show that the low velocity thin layers just below the surface and anelastic attenuation have significant effect on surface seismic record. The 2-D modeling of near surface structure beneath a 2-D refraction survey line is carried out. The comparison of the observed data with theoretical waveforms is performed. The characters in the observed data can be explained by a subsurface model constructed by P-wave traveltime tomography.
by Koichi Hayashi.
S.M.
190 leaves
11865902 bytes
11865660 bytes
application/pdf
application/pdf
eng
Massachusetts Institute of Technology
M.I.T. theses are protected by copyright. They may be viewed from this source for any purpose, but reproduction or distribution in any format is prohibited without written permission. See provided URL for inquiries about permission.
http://dspace.mit.edu/handle/1721.1/7582
Earth, Atmospheric, and Planetary Sciences.
Variable grid finite-difference modeling including surface topography
Thesis
application/pdf
false
Publication
unknown
unknown
<Publication xmlns="https://www.openaire.eu/cerif-profile/1.1/" id="895c9041-81f6-4348-92f9-d7a74b8cc05d">
<Type xmlns="https://www.openaire.eu/cerif-profile/vocab/COAR_Publication_Types">http://purl.org/coar/resource_type/c_1843</Type>
<Language>eng</Language>
<Title>Variable grid finite-difference modeling including surface topography</Title>
<PublishedIn>
<Publication>
</Publication>
</PublishedIn>
<PublicationDate>1999</PublicationDate>
<Authors>
<Author>
<DisplayName>Hayashi, Koichi, 1967-</DisplayName>
<Affiliation>
<OrgUnit>
</OrgUnit>
</Affiliation>
</Author>
</Authors>
<Editors>
</Editors>
<Publishers>
<Publisher>
<DisplayName>Massachusetts Institute of Technology</DisplayName>
<OrgUnit />
</Publisher>
</Publishers>
<License>http://dspace.mit.edu/handle/1721.1/7582</License>
<Keyword>Earth, Atmospheric, and Planetary Sciences.</Keyword>
<Abstract>We have developed a two-dimensional viscoelastic finite-difference modeling method for highly complex surface topography and subsurface structures. Realistic modeling of seismic wave propagation in the near surface region is complicated by many factors, such as strong heterogeneity, topographic relief and large attenuation. In order to account for these complications, we use a velocity-stress staggered grid and employ an 0(2,4) accurate viscoelastic finite-difference scheme. The implementation includes an irregular free surface condition for topographic relief and a variable grid technique in the shallow parts of the model. Several methods of free surface condition are bench marked, and an accurate and simple condition is proposed. In the proposed free surface condition, stresses are calculated so that the shear and normal stresses perpendicular to the boundary are zero. The calculation of particle velocities does not involve any specific calculations, and the particle velocities are set to zero above the free surface. A stable variable grid method is introduced, where we use a three times finer grid in the near surface or low velocity region compared to the rest of the model. In order to reduce instability, we apply averaging or weighting to the replacement of the coarse grid components within the fine grid. The method allows us to avoid any limitation of the shape of the grid size boundary. Numerical tests indicate that approximately ten grid-points per shortest wavelength with the variable grid method results in accurate calculations. The method requires a stair-shaped discretization of a free surface. We investigated the stair-shaped structures, and found that the cause of the dispersion from irregular free surface is mainly a numerical error due to the large grid sizes rather than the Rayleigh waves scattering due to the stair-shaped boundary. The finite-difference modeling is applied to the investigation of near surface wave propagation. Several numerical simulations are performed to show the characters of wave propagation in the near surface region. The simulations show that the low velocity thin layers just below the surface and anelastic attenuation have significant effect on surface seismic record. The 2-D modeling of near surface structure beneath a 2-D refraction survey line is carried out. The comparison of the observed data with theoretical waveforms is performed. The characters in the observed data can be explained by a subsurface model constructed by P-wave traveltime tomography.</Abstract>
<Access xmlns="http://purl.org/coar/access_right"
>
</Access>
</Publication>