9:40 AM - 9:55 AM
[SSS11-03] Comparison of the 3D FD simulation and statistical methods for the scalar wavelet propagation through random media
Keywords:Scattering, Finite Difference Simulation, Markov Approximation, Born Approximation
In this study, we conduct 3D FD simulations of the scalar wavelet propagation by using the Earth Simulator, which is a vector supercomputer managed by JAMSTEC. The target frequency is 1.5Hz and the grid spacing is 0.08 km. The average propagation velocity is 4 km/s, so there are 33 grids per a wavelength, which ensure to prevent grid dispersion. Precisions of the differentiation are 4th and 2nd order in space and time, respectively. A model medium is a cube of 307 km each side and the number of grids for each axis is 3840. In order to realize such a large random medium, we smoothly merge small-size random media generated by using different random seeds. We use 6 merged random media in total. We set an isotropic source of a Ricker wavelet at the center of the medium and receivers at propagation distances of 25, 50, 75 and 100km. The velocity fluctuation is characterized by an exponential autocorrelation function where the fractional fluctuation is 0.05 and the correlation distances (a) are 1 or 5 km. First, we compare the stacked mean square (MS) envelopes of FD simulations with those obtained by RTE+Born and those by the improved Markov approximation (Sato, 2016). For the correlation distance of 1km (akc=2.3, where kc is the centered wavenumber), RTE+Born envelopes adequately fit with FD envelopes in entire lapse time. On the other hand, Markov envelopes well model FD envelopes around the peak for the case of the correlation distance of 5 km (akc=12). In this case, late coda of FD envelopes are well reproduced by RTE+Born. Next we investigate the distribution of squared amplitudes of each trace. We find that the distribution changes from the log normal distribution at the onset to the exponential distribution at coda. For the case of the correlation distance of 1 km, squared amplitudes show the exponential distribution soon after the onset. This means that the random scattered waves are dominant. According to the FK analysis of the FD traces, we found that the squared amplitudes obey the exponential distribution even when the energy flux of scattered waves is not isotropic.