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

To estimate the S-wave velocity structure beneath the KS31 broad-band station in Wonju, Korea, we used $H-{\kappa}$ stacking and joint inversion of receiver functions and surface-wave dispersion curves derived from 297 teleseismic events (Mw > 5.5) recorded during the period between 2002 and 2009. We thereby determined that the average depth to a nearly flat Moho is $32.4{\pm}0.5\;km$ within tens of kilometer radius of the seismic station. For the crust at this location, we estimate an average shear-wave velocity of 3.69 km/s and a ratio of P- to S-wave velocities, $V_p/V_s$, of $1.72{\pm}0.04$, as is typical for continental crust. A negative phase in the receiver functions at 1 s indicates the presence of a shearwave low velocity layer in a depth interval of 10 to 18 km in the upper crust beneath the KS31 station.

• 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.15 no.1
• /
• pp.1-7
• /
• 2012

To evaluate methods of determining near-surface shear-wave velocities (${\nu}_s$), we derived dispersion curves of Rayleigh waves generated by both passive and active sources in Chuncheon, Korea. Microtremors were recorded for 5 minutes in each of four triangular arrays with radii of 5 ~ 40 m. Those data were analyzed using the Spatial Autocorrelation method. Rayleigh waves were also generated by a hammer source and recorded in the same area for 2 s using 24 4.5-Hz geophones. Multichannel Analysis of Surface Waves was applied to those data. Velocity spectra were derived with relatively high signal-to-noise ratios in the frequency ranges of 7 ~ 19 and 11 ~ 50 Hz for the microtremors and synthetically generated Rayleigh waves, respectively. The resultant dispersion curves were combined as one and then input to inversion to derive shear wave velocities that were compared with a lithology log from a nearby well. Shearwave velocities in the top soil and soft-rock layers are almost constant with values of 221 and 846 m/s, respectively; while the inverse-modeled ${\nu}_s$ increases linearly in the gravelly sand, cobbles, and weathered-rock layers. If rock type is classified based on shear-wave velocity, the inversion-derived boundary between weathered-rock and soft rock may be about 5 m deeper than in the well log.

• 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.17 no.2
• /
• pp.66-73
• /
• 2014

To estimate the shear-velocity ($v_s$) structure beneath the WIZ station on White Island in New Zealand, we applied receiver function (RF) inversion and H-${\kappa}$ stacking methods to 362 teleseismic events (Mw > 5.5) recorded during April 20, 2007 to September 6, 2013. Using 71 RFs with errors less than 20% after 200 iterative computations, we determined that the depth to Moho of $v_s$ = 4.35 km/s is $24{\pm}1km$ within a 15 km radius of the station. In an 1-d $v_s$ model derived by RF inversions, a 4-km thick low-velocity layer (LVL) at depths of 18 ~ 22 km was identified in the lower crust. This LVL, which is 0.15 km/s slower than the rocks above and below it, may indicate the presence of a deep magma reservoir. The H-${\kappa}$ stacking method yielded an estimate of the depth to the Moho of 24.5 km, which agrees well with the depth determined by RF inversions. The low $v_p/v_s$ ratio of 1.64 may be due to the presence of gas-filled rock or hot crystallizing magma.

• Geophysics and Geophysical Exploration
• /
• v.17 no.2
• /
• pp.74-87
• /
• 2014

Analysis of multi-channel seismic reflection profiles collected from the Ulleung Basin reveals that the Quaternary sequence consists of four stratigraphic units separated by erosional unconformities. Individual stratigraphic unit includes eighteen mass transport deposits which are variable in geometric characteristics and spatial distribution. Each mass transport deposit on the seismic profile is acoustically characterized by chaotic or transparent seismic facies, and shows wedge or lens-shaped external geometry. The mass transport deposits, which comprise a succession of stacked wedges, mainly occur on the southern slope, and their thickness gradually decreases toward the basin plain. The time structure map of erosional unconformities shows that a tectonic-induced structural high and troughs toward the northwest and northeast are developed at the central part of the basin. Based on the isochron map, the mass transport deposits, originated from southern part of the study area, transported to the basin plain and can be divided into two groups by the structural high. Consequently, the mass transport deposits within the Quaternary sequence in the Ulleung Basin are largely controlled by the large amounts of sediment supply, dissociation of gas hydrate during the lowstands, and central structural high.

• Geophysics and Geophysical Exploration
• /
• v.17 no.2
• /
• pp.104-113
• /
• 2014

Electric resistivity monitoring was applied to evaluate the soft ground in reclaimed land in order to figure out the applicability of physical prospecting. For this, electrical resistivity monitoring data were acquired for total three months and analyzed those data with independent inversion, time-lapse inversion, and 4D inversion methods. The result was compared for various inversion methods so as to figure out what showed the soft soil most properly. Moreover, drilling and CPT(Cone Penetration Test) data were also used in order to find out if each of those inversion methods could distinguish either bed rock or the soft soil clearly. And according to the result, time-lapse inversion showed less inversion artifacts than independent inversion, so it could indicate the soft soil better. If data gained for a longer period than three months are used, 4D inversion has been found to be a more efficient analysis method than the time-lapse inversion method. Electrical resistivity monitoring on the soft soil has been found to be a useful method that can analyze the spatio-temporal electric state of the ground serially.

• Geophysics and Geophysical Exploration
• /
• v.15 no.3
• /
• pp.144-149
• /
• 2012

Korea Polar Research Institute (KOPRI) has successfully installed two autonomous very broadband three-component seismic stations at the King George Island (KGI), Antarctica, during the 24th KOPRI Antarctic Summer Expedition (2010 ~ 2011). The seismic observation system is originally designed by the Incorporated Research Institutions for Seismology Program for Array Seismic Studies of the Continental Lithosphere Instrument Center, which is fully compatible with the Polar Earth Observing Network seismic system. The installation is to achieve the following major goals: 1. Monitoring local earthquakes and icequakes in and around the KGI, 2. Validating the robustness of seismic system operation under harsh environment. For further intensive studies, we plan to move and install them adding a couple more stations at ice shelf system, e.g., Larsen Ice Shelf System, Antarctica, in 2013 to figure out ice dynamics and physical interaction between lithosphere and cryosphere. In this article, we evaluate seismic station performance and characteristics by examining ambient noise, and provide operational system information such as frequency response and State-Of-Health information.

• Geophysics and Geophysical Exploration
• /
• v.15 no.3
• /
• pp.136-143
• /
• 2012

In this case study, we present the various and consistent processing techniques for the reasonable interpretation of aeromagnetic data. In the processing stage, we especially focused on the three major respects. First, in the low latitude area, severe artifacts are occurred as a result of reduction to the pole technique. To overcome this problem, variable alternative methods were investigated. From the comparison of each technique, we concluded that energy balancing method gives more fruitful result. Second, because of limited a priori information, it is nearly impossible to employ detailed geological survey due to wide and thick spreading of soils in the survey area. So we especially investigated the new techniques such as extracting slope, curvature and aspect information mainly used in GIS field as well as conventional methods. Finally, by using the Euler deconvolution, we extracted the depth information on the magnetic anomalous body. From the synthetic analysis between depth information and previous discussed results, the detailed future survey area was proposed. We think that a series of processing techniques discussed in this study may perform an important role in the domestic and abroad resource development project as a useful guideline.

• Geophysics and Geophysical Exploration
• /
• v.15 no.3
• /
• pp.129-135
• /
• 2012

The simultaneous-source full waveform inversion improves the applicability of full waveform inversion by reducing the computational cost. Since this technique adopts simultaneous multi-source for forward modeling, unwanted events remain in the residual seismograms when the receiver geometry of field acquisition is different from that of numerical modeling. As a result, these events impede the convergence of the full waveform inversion. In particular, the streamer-type data with limited offsets is the most difficult data to apply the simultaneous-source technique. To overcome this problem, the global-correlation-based objective function was suggested and it was successfully applied to the simultaneous-source full waveform inversion in time domain. However, this method distorts residual wavefields due to the modified objective function and has a negative influence on the inversion result. In addition, this method has not been applied to the frequency-domain simultaneous-source full waveform inversion. In this paper, we apply a timedamping function to the observed and modeled data, which are used to compute global correlation, to minimize the distortion of residual wavefields. Since the damped wavefields optimize the performance of the global correlation, it mitigates the distortion of the residual wavefields and improves the inversion result. Our algorithm incorporates the globalcorrelation-based full waveform inversion into the frequency domain by back-propagating the time-domain residual wavefields in the frequency domain. Through the numerical examples using the streamer-type data, we show that our inversion algorithm better describes the velocity structure than the conventional global correlation approach does.