• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.241-251
• /
• 2007

The observed spectra from Odaesan earthquake were fitted to a point-source spectral model to evaluate the source spectrum and spatial features of the modelling error. The source spectrum was calculated by removing from the observed spectra the path and site dependent responses (Yun, 2007) that were previously revealed through an inversion process applied to a large accumulated spectral dataset. The stress drop parameter of one-corner Brune's ${\omega}^2$ source model fitted to the estimated source spectrum was well predicted by the scaling relation between magnitude and stress drop developed by Yun et al. (2006). In particular, the estimated spectrum was quite comparable to the two-corner source model that was empirically developed for recent moderate earthquakes occurring around the Korean Peninsula, which indicates that Odaesan earthquake is one of typical moderate earthquakes representative of Korean Peninsula. Other features of the observed spectra from Odaesan earthquake were also evaluated based on the commonly treated random error between the observed data and the estimated point-source spectral model. Radiation pattern of the error according to azimuth angle was found to be similar to the theoretical estimate. It was also observed that the spatial distribution of the errors was correlated with the geological map and the $Q_0$ map which are indicatives of seismic boundaries.

• Geophysics and Geophysical Exploration
• /
• v.18 no.2
• /
• pp.74-84
• /
• 2015

To yield physically meaningful images through elastic reverse-time migration, the wavefield separation which extracts P- and S-waves from reconstructed vector wavefields by using elastic wave equation is prerequisite. For expanding the application of the elastic reverse-time migration to anisotropic media, not only the anisotropic modelling algorithm but also the anisotropic wavefield separation is essential. The anisotropic wavefield separation which uses pseudo-derivative filters determined according to vertical velocities and anisotropic parameters of elastic media differs from the Helmholtz decomposition which is conventionally used for the isotropic wavefield separation. Since applying these pseudo-derivative filter consumes high computational costs, we have developed the efficient anisotropic wavefield separation algorithm which has capability of parallel computing by using GPUs (Graphic Processing Units). In addition, the highly efficient anisotropic elastic reverse-time migration algorithm using MPI (Message-Passing Interface) and incorporating the developed anisotropic wavefield separation algorithm with GPUs has been developed. To verify the efficiency and the validity of the developed anisotropic elastic reverse-time migration algorithm, a VTI elastic model based on Marmousi-II was built. A synthetic multicomponent seismic data set was created using this VTI elastic model. The computational speed of migration was dramatically enhanced by using GPUs and MPI and the accuracy of image was also improved because of the adoption of the anisotropic wavefield separation.

• Tunnel and Underground Space
• /
• v.23 no.3
• /
• pp.180-191
• /
• 2013

The purpose of this research was to find an optimal quarrying for marble by analyzing the applicability and the work efficiency of a chain saw machine newly introduced in the underground Baekwoon mine. From the test results of the physical properties of Baekwoon marble, which affects the efficiency of rock cutting, it was found to have similar physical characteristics as the ones which are now being produced in the other areas in Korea. And especially it shows isotropic property, which can be thought to be advantageous as a dimensional stone. To check the long-term quality of the marble as a stone material, several tests such as corrosion resistance test and abrasion test were carried out. It was found to be vulnerable to acid rain with decrease of weight and seismic wave velocity after applying artificial rain at pH 5.6 for 50 times. The percentage of wear from abrasion test was 22.67%. The working time and cutting speed of the chain saw machine were recorded and analyzed during the test-run at the quarry. The overall work cycle was assorted into 9 unit operations and the operating time per each unit was drawn. The operating times for the two cutting patterns, which could be possibly applicable to the work site, were compared. The results indicated that the pattern B, that the cutting sequence was set to minimize the movement of the machine, showed 6% less working hours than the pattern A, which first cuts the outer boundary. With cutting pattern analysis, the ore body in the Baekwoon mine was 3 dimensionally modeled and a quarrying plan considering the existing conditions of the marble was suggested.

• Geophysics and Geophysical Exploration
• /
• v.8 no.1
• /
• pp.41-49
• /
• 2005

We have developed a random heterogeneous velocity model with bimodal distribution in methane hydrate-bearing Bones. The P-wave well-log data have a von Karman type autocorrelation function and non-Gaussian distribution. The velocity histogram has two peaks separated by several hundred metres per second. A random heterogeneous medium with bimodal distribution is generated by mapping of a medium with a Gaussian probability distribution, yielded by the normal spectral-based generation method. By using an ellipsoidal autocorrelation function, the random medium also incorporates anisotropy of autocorrelation lengths. A simulated P-wave velocity log reproduces well the features of the field data. This model is applied to two simulations of elastic wane propagation. Synthetic reflection sections with source signals in two different frequency bands imply that the velocity fluctuation of the random model with bimodal distribution causes the frequency dependence of the Bottom Simulating Reflector (BSR) by affecting wave field scattering. A synthetic cross-well section suggests that the strong attenuation observed in field data might be caused by the extrinsic attenuation in scattering. We conclude that random heterogeneity with bimodal distribution is a key issue in modelling hydrate-bearing Bones, and that it can explain the frequency dependence and scattering observed in seismic sections in such areas.

• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.98-104
• /
• 2011

Numerical simulation in exploration geophysics provides important insights into subsurface wave propagation phenomena. Although elastic wave simulations take longer to compute than acoustic simulations, an elastic simulator can construct more realistic wavefields including shear components. Therefore, it is suitable for exploration of the responses of elastic bodies. To overcome the long duration of the calculations, we use a Graphic Processing Unit (GPU) to accelerate the elastic wave simulation. Because a GPU has many processors and a wide memory bandwidth, we can use it in a parallelised computing architecture. The GPU board used in this study is an NVIDIA Tesla C1060, which has 240 processors and a 102 GB/s memory bandwidth. Despite the availability of a parallel computing architecture (CUDA), developed by NVIDIA, we must optimise the usage of the different types of memory on the GPU device, and the sequence of calculations, to obtain a significant speedup of the computation. In this study, we simulate two- (2D) and threedimensional (3D) elastic wave propagation using the Finite-Difference Time-Domain (FDTD) method on GPUs. In the wave propagation simulation, we adopt the staggered-grid method, which is one of the conventional FD schemes, since this method can achieve sufficient accuracy for use in numerical modelling in geophysics. Our simulator optimises the usage of memory on the GPU device to reduce data access times, and uses faster memory as much as possible. This is a key factor in GPU computing. By using one GPU device and optimising its memory usage, we improved the computation time by more than 14 times in the 2D simulation, and over six times in the 3D simulation, compared with one CPU. Furthermore, by using three GPUs, we succeeded in accelerating the 3D simulation 10 times.