• 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.11 no.2
• /
• pp.137-147
• /
• 2008

Fault rupture directivity of the Odaesan earthquake, which was inferred to be the main cause of the high PGAvalue (> 0.1 g) unusually observed at the near-source region, was analyzed by using the data from the nearby (R < 100 km) dense seismic stations. The Boatwright's method (2007) was adopted for this purpose in which the azimuth and takeoff angle of the unilateral rupture directivity function could be estimated based on the relative peak ground-motions of seismic stations resulting from the nature of the rupture directivity. In this study, the approximate values of the relative peak ground-motions was derived from the difference between the log residuals of the point-source spectral model (Boore, 2003) for the main and secondary events based on the Random Vibration Theory. In this derivation, the spectral difference for a frequency range between the source corner frequencies of main and secondary events was considered to reflect only the effect of the fault directivity. The inversion result of the model parameters for the fault directivity function showed that the fault-plane of NWW-SEE direction dipping steeply to the North with high rupture velocity near upward in SE direction is responsible for the observed high level of ground-motion at the near-source region.

• 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.

• Geophysics and Geophysical Exploration
• /
• v.18 no.2
• /
• pp.39-53
• /
• 2015

Analysis of multi-channel seismic reflection profiles acquired from the southwestern margin of Ulleung Basin reveals that the cut and fill structures, which show U-shaped or V-shaped morphology, occur on variable size. The cut and fill structure mostly consists of fine-grained sediments on the well data and is characterized by transparent or semitransparent seismic facies on the seismic section. Such cut and fill structures dominantly occur in the syn-compressional megasequence (MSQ3), which was deposited during basin deformation of late Miocene, among the four megasequences of the study area. These cut and fill structures can be divided into three groups based on their size and formation time. The cut and fill structures of Group I were formed when Dolgorae structure was active, and occurred on a small scale. The cut and fill structures of group II were formed when both Dolgorae structure and Gorae V structure were active, and the number and size of those increased compared with group I. The cut and fill structures of group III were formed when Dolgorae structure was weaken gradually but Gorae V structure kept active, and the number and size of those decreased in comparison with group II. Consequently the cut and fill structures in the southwestern margin of Ulleung basin are interpreted as submarine canyon based on spatial distribution, size and fill sediment. They were controlled by the tectonic movement in response to basin closure and tectonic-induced sediment supply variation.

• Geophysics and Geophysical Exploration
• /
• v.15 no.1
• /
• pp.16-22
• /
• 2012

In this study, we setup a ground-based scatterometer using an antenna of which the center frequency is 9.5 GHz (X-band), and measured radar backscatterings from air/ice and ice/water interfaces to extract ice thickness. Both of air/ice and ice/water interfaces make strong radar backscatterings and so we can clearly identify two peaks in measured data by scatterometer. By using the distance of two peaks and refractive index of ice, we confirmed that it is possible to measure ice thickness. Field survey was performed at the downstream of Jiam River flowing into Chuncheon Lake. We measured radar backscattering from river ice along a survey path and extracted ice thickness. The ice thickness map of the downstream of Jiam River was produced by using kriging which is one of well known interpolation methods. The ice thickness was about 50 cm along the mainstream while ice was thin as 30 ~ 40 cm at a fast-flowing meander. Ice thickness was particularly thinner at some locations than that of surrounding areas even in the mainstream region of constant flow. This was because of impurities in ice or artificially formed refrozen holes after fishing. We expect that this study helps to expand utilization field of X-band SAR and airborne scatterometer system.

• Geophysics and Geophysical Exploration
• /
• v.14 no.3
• /
• pp.220-226
• /
• 2011

In the elastic wave equations, both horizontal and vertical displacements are defined. Since we can measure both the horizontal and vertical displacements in field acquisition, these displacements compose a displacement vector. In this study, we propose a frequency-domain elastic waveform inversion technique taking advantage of the magnitudes of displacement vectors to define objective function. When we apply this displacement-vector objective function to the frequency-domain waveform inversion, the inversion process naturally incorporates the back-propagation algorithm. Through the inversion examples with the Marmousi model and the SEG/EAGE salt model, we could note that the RMS error of the solution obtained by our algorithm decreased more stably than that of the conventional method. Particularly, the density of the Marmousi model and the low-velocity sub-salt zone of the SEG/EAGE salt model were successfully recovered. Since the gradient direction obtained from the proposed objective function is numerically unstable, we need additional study to stabilize the gradient direction. In order to perform the waveform inversion using the displacementvector objective function, it is necessary to acquire multi-component data. Hence, more rigorous study should be continued for the multi-component land acquisition or OBC (Ocean Bottom Cable) multi-component survey.

• Geophysics and Geophysical Exploration
• /
• v.17 no.2
• /
• pp.49-57
• /
• 2014

3D magnetic inversion, based on the assemblage of 2D forward modeling and inversion as a practical technique to reflect the a priori information, was conducted to investigate the spatial distribution features of black-shale related and pyrometamorphic uranium deposit and several lithological units of Ogcheon Super Group in an area of Geumsan. By using the 3D visualization technique with suitable susceptibility interval and horizontal slice map, the spatial distribution of magnetic susceptibility corresponded to the black shale related uranium bearing lithological units, Black Slate member was well coincided with a information of uranium deposit. Also, even though it is indirect indicator for the detetction of uranium deposits interbedded in Gray Hornfels member, spatial susceptibility distribution which shows the south-east magnetic linearment corresponding to the Majeon-ri formation and Dark Gray Slate were matched well. From this investigation, we inferred that maximum depth extension which Black Slate member can be separately recognized with respect to adjacent Dark Gray Slate with strong magnetic susceptibility anomaly is about 150m with reference elevation level of 306m. In addition, Majeon-ri formation located south of Black Slate member revels relatively high magnetic susceptibility range but shows high spatial susceptibility fluctation. And, as an intrusive rocks, Jurassic Biotite Granite shows relatively low magnetic susceptibility characteristics. On the contrary, Cretaceous granite distributed in soutthern part of the study area shows the relatively high susceptibility distribution.

• Geophysics and Geophysical Exploration
• /
• v.6 no.4
• /
• pp.165-170
• /
• 2003

Recently, resistivity surveys have been frequently carried out over the irregular terrain such as mountainous area. Such an irregular terrain itself can produce significant anomalies which may lead to misinterpretations. In this study, topographic effects in resistivity survey were studied using the physical scale modeling as well as the numerical one adopting finite element method. The scale modeling was conducted at a pond, so that we could avoid the edge effect, the inherent problem of the scale modeling conducted in a water tank in laboratory. The modeling experiments for two topographic features, a ridge and a valley with various slope angles, confirmed that the results by the two different modeling techniques coincide with each other fairly well for all the terrain models. These experiments adopting dipole-dipole array showed the distinctive terrain effects, such that a ridge produces a high apparent resistivity anomaly at the ridge center flanked by zones of lower apparent resistivity. On the other hand, a valley produces the opposite anomaly pattern, a central low flanked by highs. As the slope of a terrain model becomes steeper, the terrain-induced anomalies become stronger, and moreover, apparent resistivity can become even negative for the model with extremely high slope angle. All the modeling results led us to the conclusion that terrain effects should be included in the numerical modeling and/or the inversion process to interpret data acquired at the rugged terrain area.

• Journal of the Korean Geophysical Society
• /
• v.7 no.2
• /
• pp.99-112
• /
• 2004

The Basin, a marginal basin located between the Antarctic Peninsula and the South Shetland Islands, is consist of three small basins, the Central, Eastern, Western Basins. Seismic data obtained on December 1995 show well-defined spreading ridges, basement highs, faults, morphology of the basin, distribution of sediments, crustal and sedimentary deformation, diapirs, and contourites. The main spreading axis of the Central Bransfield Basin connecting Deception and Bridgeman Islands continues up to the central part of the Eastern Basin, whereas deep basin covered by thick sediments without any spreading structures develops in the northeastern part. This indicates that back-arc spreading along the axis of the Bransfield Basin has been taken place in the southwestern part of the Eastern Basin, not in the northeastern part. Many NW-SE trending faults perpendicular to the axis of the basin would be related with strike-slip movement of the Shackleton Fracture. Zone. Extensinal strutures like deep basin without any spreading structures in the northeastern part, normal faults and diapirs on both continental slopes of the Eastern Basin would be formed by extension as a consequence of the sinistral movement between Antarctic and the Scotia plates.

• Geophysics and Geophysical Exploration
• /
• v.8 no.1
• /
• pp.18-25
• /
• 2005

We derive refraction statics for seismic data recorded in a hard rock terrain, in which there are large and rapid variations in the depth of weathering. The statics corrections range from less than 10 ms to more than 70 ms, often over distances as short as 12 receiver intervals. This study is another demonstration of the importance in obtaining accurate initial refraction models of the weathering in hard rock terrains in which automatic residual statics may fail. We show that the statics values computed with a simple model of the weathering using the Generalized Reciprocal Method (GRM) and the Refraction Convolution Section (RCS) are comparable in accuracy to those computed with a more complex model of the weathering, using least-mean-squares inversion with the conjugate gradient algorithm (Taner et al., 1998). The differences in statics values between the GRM model and that of Taner et al. (1998) systematically vary from an average of 2ms to 4ms over a distance of 8.8 km. The differences between these two refraction models and the final statics model, which includes the automatic residual values, are generally less than 5 ms. The residuals for the GRM model are frequently less than those for the model of Taner et al. (1998). The RCS statics are picked approximately 10 ms later, but their relative accuracy is comparable to that of the GRM statics. The residual statics values show a general correlation with the refraction statics values, and they can be reduced in magnitude by using a lower average seismic velocity in the weathering. These results suggest that inaccurate average seismic velocities in the weathered layer may often be a source of short-wavelength statics, rather than any shortcomings with the inversion algorithms in determining averaged delay times from the traveltimes.