• Geophysics and Geophysical Exploration
• /
• v.8 no.1
• /
• pp.26-33
• /
• 2005

Inversion of multi-mode surface-wave phase velocity for shallow engineering site investigation has received much attention in recent years. A sensitivity analysis and inversion of both synthetic and field data demonstrates the greater effectiveness of this method over employing the fundamental mode alone. Perturbation of thickness and shear-wave velocity parameters in multi-modal Rayleigh wave phase velocities revealed that the sensitivities of higher modes: (a) concentrate in different frequency bands, and (b) are greater than the fundamental mode for deeper parameters. These observations suggest that multi-mode phase velocity inversion can provide better parameter discrimination and imaging of deep structure, especially with a velocity reversal, than can inversion of fundamental mode data alone. An inversion of the theoretical phase velocities in a model with a low velocity layer at 20 m depth can only image the soft layer when the first higher mode is incorporated. This is especially important when the lowest measurable frequency is only 6 Hz. Field tests were conducted at sites surveyed by borehole and PS logging. At the first site, an array microtremor survey, often used for deep geological surveying in Japan, was used to survey the soil down to 35 m depth. At the second site, linear multichannel spreads with a sledgehammer source were recorded, for an investigation down to 12 m depth. The f-k power spectrum method was applied for dispersion analysis, and velocities up to the second higher mode were observed in each test. The multi-mode inversion results agree well with PS logs, but models estimated from the fundamental mode alone show f large underestimation of the depth to shallow soft layers below artificial fill.

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

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.44 no.4
• /
• pp.503-512
• /
• 2024

In the seismic design standard, input waves for different levels of seismic performance are proposed in the form of design response spectra. At the time of establishing these standards, measured records of significant earthquakes that occurred domestically, such as the 2016 Gyeongju earthquake and the 2017 Pohang earthquake, were not included. Additionally, for the ground response analysis, shear wave velocities representing ground amplification characteristics were derived from the results of standard penetration tests (N-values) and applied in empirical formulas. This approach may not adequately capture sufficient information about the characteristics of domestic ground properties. Therefore, in this study, seismic records from the Gyeongju and Pohang earthquakes were modified to adjust the bedrock standard design response spectra. Ground response analyses were conducted using shear wave velocity profiles obtained from borehole tests in the Sejong City area. The shape of the response spectrum and ground amplification coefficient obtained from the ground response analysis were then compared with those from existing studies and seismic design standard.

• The Journal of Engineering Geology
• /
• v.23 no.4
• /
• pp.457-465
• /
• 2013

We integrated and correlated datasets from surface and subsurface geophysics, drilling cores, and engineering geology to identify geological interfaces and characterize the joints and fracture zones within the rock mass. The regional geometry of a geologically weak zone was investigated via a fence projection of electrical resistivity data and a borehole image-processing system. Subsurface discontinuities and intensive fracture zones within the rock mass are delineated by cross-hole seismic tomography and analyses of dip directions in rose diagrams. The dynamic elastic modulus is studied in terms of the P-wave velocity and Poisson's ratio. Subsurface discontinuities, which are conventionally identified using the N value and from core samples, can now be identified from anomalous reflection coefficients (i.e., acoustic impedance contrast) calculated using a pair of well logs, comprising seismic velocity from suspension-PS logging and density from logging. Intensive fracture zones identified in the synthetic seismogram are matched to core loss zones in the drilling core data and to a high concentration of joints in the borehole imaging system. The upper boundaries of fracture zones are correlated to strongly negative amplitude in the synthetic trace, which is constructed by convolution of the optimal Ricker wavelet with a reflection coefficient. The standard deviations of dynamic elastic moduli are higher for fracture zones than for acompact rock mass, due to the wide range of velocities resulting from the large numbers of joints and fractures within the zone.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.26-33
• /
• 2008

The analysis of wide-angle seismic reflection and refraction data plays an important role in lithospheric-scale crustal structure study. However, it is extremely difficult to develop an appropriate velocity structure model directly from the observed data, and we have to improve the structure model step by step, because the crustal structure analysis is an intrinsically non-linear problem. There are several subjective processes in wide-angle crustal structure modelling, such as phase identification and trial-and-error forward modelling. Because these subjective processes in wide-angle data analysis reduce the uniqueness and credibility of the resultant models, it is important to reduce subjectivity in the analysis procedure. From this point of view, we describe two software tools, PASTEUP and MODELING, to be used for developing crustal structure models. PASTEUP is an interactive application that facilitates the plotting of record sections, analysis of wide-angle seismic data, and picking of phases. PASTEUP is equipped with various filters and analysis functions to enhance signal-to-noise ratio and to help phase identification. MODELING is an interactive application for editing velocity models, and ray-tracing. Synthetic traveltimes computed by the MODELING application can be directly compared with the observed waveforms in the PASTEUP application. This reduces subjectivity in crustal structure modelling because traveltime picking, which is one of the most subjective process in the crustal structure analysis, is not required. MODELING can convert an editable layered structure model into two-way traveltimes which can be compared with time-sections of Multi Channel Seismic (MCS) reflection data. Direct comparison between the structure model of wide-angle data with the reflection data will give the model more credibility. In addition, both PASTEUP and MODELING are efficient tools for handling a large dataset. These software tools help us develop more plausible lithospheric-scale structure models using wide-angle seismic data.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2006.06a
• /
• pp.181-186
• /
• 2006

An imaging method of seismic sources using time-reversal wave propagation is presented. The method is based on the time-reversal invariance and the spatial reciprocity of the wave equation. Time-reversal wave propagation has been used to image anomalous features of a midium in medical imaging, non destructive testing and waveform tomography. Seismogram is the record whose energy is propagated from the seismic source. If time-reversed seismogram propagates back into the medium, seismic energy is concentrated at the origin time of the event and at the source location. In this work, a staggered-grid finite-difference method of the elastic wave equation is parallelized for 3-D wave propagation simulation. With numerical experiments, we show that the time-reversal imaging will enable us to explore the spatio-temporal history of complex earthquake.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.15-25
• /
• 2008

Ocean-bottom seismometer (OBS) positions are one of the key parameters in an OBS-airgun seismic survey for crustal structure study. To improve the quality of these parameters, we have developed a new method of determining OBS positions, using airgun shot data and bathymetric data in addition to available distance measurements by acoustic transponders. The traveltimes of direct water waves emitted by airgun shots and recorded by OBSs are used as important information for determining OBS locations, in cases where there are few acoustic transponder data (<3 sites). The new method consists of two steps. A global search is performed as the first step, to find nodes of the bathymetric grid that are the closest to explaining the observed direct water-wave traveltimes from airgun shots, and acoustic ranging using a transponder system. The use of precise 2D bathymetric data is most important if the bottom topography near the OBS is extremely rough. The locations of the nodes obtained by the first step are used as initial values for the second step, to avoid falling into local convergence minima. In the second step, a non-linear inverse method is executed. If the OBS internal clock shows large drift, a secondary correction for the OBS internal clock is obtained, as well as the OBS location, as final results by this method. We discuss the error and the influence of each measurement used in the determination of OBS location.

• Journal of the Korean Geophysical Society
• /
• v.3 no.1
• /
• pp.37-48
• /
• 2000

In order to investigate in-line ground motions caused by earthquakes, we examine the multicomponent complex trace analysis method (MCTAM) for the synthetic data and apply it to real earthquake data. An experimental result for synthetic data gives correct information on the arrival times, duration of individual phases, and approaching angles for body waves. Rayleigh waves are also easily identified with the MCTAM. A deep earthquake with magnitude of 7.3 was chosen to test various polarization attributes of ground motions. For P waves, instantaneous phase difference between the vertical and the in-line horizontal components ${\phi}(t)$, instantaneous reciprocal ellipticity ${\rho}(t)$, and approaching angle ${\tau}(t)$ are computed to be ${\pm}180^{\circ},\;0{\sim}0.25,\;and\;-30^{\circ}{\sim}-45^{\circ}$, respectively. For S waves, ${\phi}(t)$ tends to vary while ${\rho}(t)$ have values of $0{\sim}0.3\;and\;{\tau}(t)$ remains near vertical, respectively. A relatively low frequency signal registered just prior to the S wave event is interpreted as a P-wave phase based on its polarization characteristics. Velocities of P and S waves are computed to be 8.633 km/s and 4.762 km/s, and their raypath parameters 0.074 s/km and 0.197 s/km. Dynamic Poisson's ratio is obtained as 0.281 from the velocities of P and S waves.

• Geophysics and Geophysical Exploration
• /
• v.12 no.3
• /
• pp.225-232
• /
• 2009

As part of seismic experiments investigating crustal velocity structures of the Korean peninsula, permanent (fixed) seismographs of the Korea Meteorological Administration (KMA) network recorded seismic signals from four and eight large explosions in Korean Crustal Research Team (KCRT) profiles shot in 2004 and 2008, respectively. Among the seismograms recorded by 43 velocity sensors and 103 accelerometers at KMA stations distributed throughout the southern Korean Peninsula, 156 records with epicentral distances less than 120 km and high signal-to-noise ratios were analyzed to determine velocity anisotropy of the Pg phase. Relative elevation corrections of -101.6 to 105.3 ms were made using velocity information derived from the 2004 KCRT profile data and differences in elevation between the permanent KMA stations and the temporary stations in the KCRT profiles at the same source-receiver offsets. To remove site effects, receiver-station corrections of -89.6 to 192.2 ms were additionally made to the KMA station data by subtracting the average differences in traveltimes between KMA stations and portable stations at the same offsets for all available shots with different azimuths. With the exception of anomalously fast velocities along trends of the Chugaryeong fault zone and the Okchon fold belt and anomalously slow velocities in the regions of high terrestrial heat near Yeongduk and Ulsan, the analysis of crustal velocity anisotropy using the Pg phase indicates overall isotropy in the southern half of the Korean peninsula.

• The Korean Journal of Quaternary Research
• /
• v.17 no.1
• /
• pp.7-16
• /
• 2003

A 31m-long sediment core (SSDP-102) was taken from the inner shelf (about 40m water-depth) off the northwestern coast of the Korea Strait. Detailed lithofacies and organic-geochemical analyses were performed to establish a high-resolution stratigraphy in the Korea Strait shelf and to reconstruct the paleoenvironmental changes associated with the Holocene marine transgression. The stratigraphic framework of the core was primarily established using 6 AMS $^{14}C$ ages. The sedimentary record of the core SSDP-102 allows for the reconstruction of the paleoenvironmental changes during the last 12.1 ka BP. According to the high-resolution seismic reflection profiles, lithofacies and organic-geochemical data, the core SSDP-102 can be divided into three units (III to I in ascending order) above the acoustic basement. The three units reflect distinct changes of depositional environments resulted from the post-glacial marine transgression. Therefore, it is suggested that three phases of sea-level change have occurred within the inner shelf of the Korea Strait following the Holocene marine transgression. (1) estuarine environments from ca. 12.1 to 6.2 ka BP; (2) near-shore environments with a period of decreased rising of sea level between 6.2 and 5.1 ka BP; (3) near-shore to modem marine environments after 5.1 ka BP. In particular, the present marine conditions influenced by the warm Tsushima Current have been gradually established after ca. 5.1 ka BP.