• 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.9 no.4
• /
• pp.304-315
• /
• 2006

We have examined the applicability of f-k analysis to the GPR direct wave measurement for water content to characterize vadose zone condition. When the vadose zone consists of a dry surface layer over wet substratum, we obtained f-k spectra where most of the energy is bounded by the air and dry soil velocities. In this case, dry soil velocity was successfully estimated by using high frequency data. On the other hands, when wet soil overlies dry substratum, the f-k spectra show a contrasting response where most of the energy travels with the velocity bounded by dry and wet soil velocities. In this case, the radar waves are trapped and guided within wet soil layer, exhibiting velocity dispersion. By adopting modal propagation theory, we could formulae a simple inversion code to find two layer's dielectric constants as well as layer thickness. By inverting the velocity dispersion curve obtained from f-k spectra of synthetic modeling data, we could obtain good estimates of dielectric constants of each layer as well as first layer thickness. Moreover, we could obtain more accurate results by including the higher mode data. We expect this method will be useful to get the quantitative property of real subsurface when the field condition is similar.

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

In exploration seismology, the Kirchhoff hyperbola has been successfully used to migrate reflection seismo-grams. The mathematical basis of Kirchhoff hyperbola has not been clearly defined and understood for the application of prestack or poststack migration. The travel time from the scatterer in the subsurface to the receivers (exploding reflector model) on the surface can be a kinematic approximation of Green's function when the source is excited at position of the scatterer. If we add the travel time from the source to the scatterer in the subsurface to the travel time of exploding reflector model, we can view this travel time as a kinematic approximation of the partial derivative wavefield with respect to the velocity or the density in the subsurface. The summation of reflection seismogram along the Kirchhoff hyperbola can be evaluated as an inner product between the partial derivative wavefield and the field reflection seismogram. In addition to this kinematic interpretation of Kirchhoff hyperbola, when we extend this concept to shallow refraction seismic data, the stacking of refraction data along the straight line can be interpreted as a measurement of an inner product between the first arrival waveform of the partial derivative wavefield and the field refraction data. We evaluated the Kirchhoff hyperbola and the straight line for stacking the refraction data in terms of the first arrival waveform of the partial derivative wavefield with respect to the velocity or the density in the subsurface. This evaluation provides a firm and solid basis for the conventional Kirchhoff migration and the straight line stacking of the refraction data.

• Geophysics and Geophysical Exploration
• /
• v.15 no.2
• /
• pp.85-91
• /
• 2012

Multiple Lapse Time Window (MLTW) analysis for obtaining intrinsic attenuation value require numerous data without directional bias to compensate focal mechanism. The first window of MLTW, therefore, shows large deviation in fitting smoothed theoretical curve. The information on the focal mechanism may reduce burdens of number and distribution. This study combined algorithm of computing focal mechanism to DSMC method by Yoshimoto (2000). However, the MLTW method based on the numerous data was not applicable to this study, because of the limited data to the almost same fault plane solution. This study showed that the available data was too insufficient to construct smoothed theoretical curve, although the deviation of the first window was improved. Instead of conventional solution by more data, the study seems to be needed for new constraints to obtain smoothed curve.

• Korean Journal of Remote Sensing
• /
• v.27 no.1
• /
• pp.17-24
• /
• 2011

Global warming and climatic changes due to human activities impact on marine and terrestrial ecosystems, which feedbacks to climate system. These negative feedbacks amplify or accelerate again global climate change. In particular, it is important to analyze vegetation change. This study attempts to analyze quantitatively vegetation change in Korea peninsula by using harmonic analysis. Harmonic-Analysis based on Fourier Transform is the method to effectively demonstrate for time series data. Especially, Harmonic-Analysis is very suitable method to analyze vegetation change because the vegetation repeats the cycle growth and extinction every year. The result of harmonic-analysis shows vegetation change as time passes. In this study, SPOTNEGETATION S10 MVC NDVI data was used during last 10 years (1999-2008) in Korea Peninsula. Also, land type classification used MODIS Land Cover Map data. The study estimated that phase values moved up approximately 0.5 day per year in cropland and 0.8 day per year in forest.

• Geophysics and Geophysical Exploration
• /
• v.12 no.2
• /
• pp.199-207
• /
• 2009

From 329 earthquake data in the Korean Peninsula, we separated the intrinsic and scattering quality factor ($Q_i^{-1}$ and $Q_s^{-1}$) using the Multiple Lapse Time Window Analysis (MLTW) method. For the homogeneous velocity structure, $Q_i^{-1}$ reduces the amplitudes of both direct and coda waves; $Q_s^{-1}$ diminishes the direct wave amplitude but enlarges the coda wave amplitude. Based on this phenomenon, MLTW method analytically derives theoretical curves and obtains $Q_i^{-1}$ and $Q_s^{-1}$ by least square fit with observation curves. This study is the first approach for the seismic stable region by MLTW method, and show that $Q_i^{-1}$ and $Q_s^{-1}$ in the Korean Peninsula are very low at lower than frequencies of 5 Hz. This low value seems to be related to the inactive tectonism of the Korean Peninsula.

• Geophysics and Geophysical Exploration
• /
• v.3 no.4
• /
• pp.119-124
• /
• 2000

A new fast algorithm for travel time calculation using mono-chromatic one-way wave equation was developed based on the delta function and the logarithms of the single frequency wavefield in the frequency domain. We found an empirical relation between grid spacing and frequency by trial and error method such that we can minimize travel time error. In comparison with other methods, travel time contours obtained by solving eikonal equation and the wave front edge of the snapshot by the finite difference modeling solution agree with our algorithm. Compared to the other two methods, this algorithm computes travel time of directly transmitted wave. We demonstrated our algorithm on migration so that we obtained good section showing good agreement with original model. our results show that this new algorithm is a faster travel time calculation method of the directly transmitted wave for imaging the subsurface and the transmission tomography.

• Geophysics and Geophysical Exploration
• /
• v.27 no.1
• /
• pp.51-56
• /
• 2024

In this paper, the vector gravity and gravity gradient tensor of an elliptical disk are derived. The vector gravity of an elliptical disk is defined by differentiating the gravitational potential due to the elliptical disk expressed by a double integral with respect to each axial direction. The vector gravity defined by the double integral is then transformed into a line integral of a closed curve along the elliptical disk boundary using the complex Green's theorem. Finally, vector gravity due to the elliptical disk is derived by 1D parametric numerical integration along the elliptical disk boundary. The xz, yz, zz components of the gravity gradient tensor due to the elliptical disk are obtained by differentiating the vector gravity with respect to vertical direction. The xx, yy, xy components are derived by differentiating the horizontal components of the vector gravity in the form of a double integral with respect to horizontal directions and then using the complex Green's theorem.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.16 no.5
• /
• pp.121-128
• /
• 2012

In this paper, nonlinear analysis for reinforced concrete structure for power transmission line is performed by considering the characteristics of the failure, which are depend on loading conditions and concrete material models. On the numerical evaluation for the failure behavior, the finite element analysis is applied. For the concrete material model, microplane model based on concrete damage is introduced. However, to describe the crack bridging effect of long and short fiber of steel fiber reinforced concrete (SFRC), tensile softening model is suggested and applied for SFRC. An numerical results by finite element technique are compared with the experiment results for box culvert specimen. Comparing on the experimental and analytical results, validity and reliability of numerical analysis are investigated.

• Geophysics and Geophysical Exploration
• /
• v.17 no.1
• /
• pp.11-20
• /
• 2014

To reveal shear-wave velocities ($v_s$) and site characterization of Wonju, Korea, Rayleigh waves were recorded at 78 sites of lower altitude using 12 to 24 4.5-Hz vertical geophones for 20 days during the period of February to September 2013. Dispersion curves of the Rayleigh waves obtained by the extended spatial autocorrelation method were inverted using the damped least-squares method to derive $v_s$ models. From these 1-D models, the average $v_s$ to a depth of 30 m ($v_s30$), $v_s$ of weathered rocks, depths to these basement rocks, and average $v_s$ of the overburden layer were derived to be $16.3{\pm}0.7m$, $576{\pm}8m/s$, $290{\pm}7m/s$, and $418{\pm}13m/s$, respectively, in the 95% confidence range. To determine adequate proxies for $v_s30$, we computed correlation coefficients of $v_s30$ with topographic slope (r = 0.46) and elevation (r = 0.43). An empirical linear relationship is presented as a combination of individually estimated $v_s30$ with weighting factors of 0.45, 0.45, and 0.1 for topographic slope, elevation, and mapped lithology, respectively. Due to a weak correlation between $v_s30$ obtained from inversion of dispersion curves and the proxy-based estimation (r = 0.50), however, the relatively large error range should be considered for applications of this relationship.