• 한국지구물리탐사학회:학술대회논문집
• /
• 2005.05a
• /
• pp.273-280
• /
• 2005

Robust magneto-telluric (MT) response function estimators are now in standard use in electromagnetic induction research. Properly devised and applied, these methods can reduce the influence of unusual data (outlier) in the response (electric field) variable, but often not sensitive to exceptional predictor (magnetic field) data, which are termed leverage points. A bounded influence estimator is described which simultaneously limits the influence of both outlier and leverage point, and has proven to consistently yield more reliable MT response function estimates than conventional robust approach. The bounded influence estimator combines a standard robust M-estimator with leverage weighting based on the statistics of the hat matrix diagonal, which is a standard statistical measure of unusual predictors. Further extensions to MT data analysis are proposed, including a establishment of data rejection criterion which minimize the influence of both electric and magnetic outlier in frequency domain based on statistical distribution of electromagnetic field. The rejection scheme made in this study seems to have an effective performance on eliminating extreme data, which is even not removed by BI estimator, in frequency domain. The effectiveness and advantage of these developments are illustrated using real MT data.

• Geophysics and Geophysical Exploration
• /
• v.22 no.4
• /
• pp.195-201
• /
• 2019

The logarithmic objective function used in waveform inversion minimizes the logarithmic differences between the observed and modeled data. Laplace-domain waveform inversions usually adopt the logarithmic objective function and the diagonal elements of the pseudo-Hessian for optimization. In this case, we apply the Levenberg-Marquardt algorithm to prevent the diagonal elements of the pseudo-Hessian from being zero or near-zero values. In this study, we analyzed the diagonal elements of the pseudo-Hessian of the logarithmic objective function and showed that there is no zero or near-zero value in the diagonal elements of the pseudo-Hessian for acoustic waveform inversion in the Laplace domain. Accordingly, we do not need to apply the Levenberg-Marquardt algorithm when we regularize the gradient direction using the pseudo-Hessian of the logarithmic objective function. Numerical examples using synthetic and field datasets demonstrate that we can obtain inversion results without applying the Levenberg-Marquardt method.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2005.09a
• /
• pp.93-116
• /
• 2005

SPT-Uphole tomography method was introducedand verified in this paper. In SPT-Uphole method, SPT (Standard Penetration Test) which is common in site investigation, was used as a source and several surface geophones in line were used as receivers. Shear wave velocity (Vs) distribution map which has triangular shape around the boring point can be obtained by tomography inversion. The factors for obtaining reliable result of SPT-Uphole tomography are exact travel time information and accurate inversion method. To establish of the SPT-Uphole tomography procedure, the most reliable method for obtaining exact travel time information and verification of tomography inversion method were studied by using theoretical travel time information and finite element method (FEM) analysis. finally, SPT-Uphole tomography method was performed at the weathered soil site in Kimje. By comparing with several boring data including SPT-N value, feasibility of this method was verified in the field.

• 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.12 no.4
• /
• pp.347-354
• /
• 2009

Physical properties of rocks are strongly dependant on details of pore micro-structures, which can be used for quantifying relations between physical properties of rocks through pore-scale simulation techniques. Recently, high-resolution scan techniques, such as X-ray microtomography and high performance computers make it possible to calculate permeability from pore micro-structures of rocks. We try to extend this simulation methodology to velocity and electrical conductivity. However, the smoothing effect during tomographic inversion creates artifacts in pore micro-structures and causes inaccurate property estimation. To mitigate this artifact, we tried to use sharpening filter and neural network classification techniques. Both methods gave noticeable improvement in pore structure imaging and accurate estimation of permeability and electrical conductivity, which implies that our method effectively removes the smoothing effect in pore structures. However, the calculated velocities showed only incremental improvement. By comparison between thin section images and tomogram, we found that our resolution is not high enough, and it is mainly responsible for the inaccuracy in velocity despite the successful removal of the smoothing effect. In conclusion, our methods can be very useful for pore-scale modeling, since it can create accurate pore structure without the smoothing effect. For accurate velocity estimation, the resolution of pore structure should be at least three times higher than that for permeability simulation.

• Geophysics and Geophysical Exploration
• /
• v.22 no.1
• /
• pp.21-28
• /
• 2019

We perform a 3D tomographic inversion using surface wave dispersion curves to obtain S-velocity model and radial anisotropy beneath Saudi Arabia. The Arabian Peninsula is geologically and topographically divided into a shield and a platform. We used event data with magnitudes larger than 5.5 and epicentral distances shorter than $40^{\circ}$ during 2008 ~ 2014 from the Saudi Geological Survey. We obtained dispersion curves by using the multiple filtering technique after preprocessing the event data. We constructed SH- and SV-velocity models and consequently radial anisotropy model at 10 ~ 60 km depths by inverting Love and Rayleigh group velocity dispersion curves with period ranges of 5 ~ 140 s, respectively. We observe high-velocity anomalies beneath the Arabian shield at 10 ~ 30 km depths and low-velocity anomalies beneath the Arabian platform at 10 km depth in the SV-velocity model. This discrepancy may be caused by the difference between the Arabian shield and the Arabian platform, that is, the Arabian shield was formed in Proterozoic thereby old and cold, while the Arabian platform is covered by predominant Paleozoic, Mesozoic, and Cenozoic sedimentary layers. Also we obtained radial anisotropy by estimating the differences between SH- and SV-velocity models. Positive anisotropy is observed, which may be generated by lateral tension due to the slab pull of subducting slabs along the Zagros belt.

• Journal of the Korean Geophysical Society
• /
• v.3 no.4
• /
• pp.251-260
• /
• 2000

Applicabilities of two numerical methods, the 2-dimensional and the 3-dimensional method, are evaluated to inverse test results obtained from the underwater SASW(Spectral -Analysis-of-Surface-Waves) method. As a result of this study, it has been found that the 2-dimensional method can supposed to be applicable for the cases where stiffness of soil layer increases gradually with depth, and the stiffness is relatively low. For the other cases, however, it has been concluded that the 3-dimensional method needs to be applied to determine realistic theoretical dispersion curves. An example is also shown that in situ soil profile underwater is estimated from experimental dispersion curves using the 3-dimensional method. As a results, it can be concluded that the underwater SASW method can be effectively applied to explore the underwater soil condition.

• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.138-146
• /
• 2007

This paper describes 2.5D induced polarization (IP) modeling and inversion algorithms using complex resistivity. The complex resistivity method has merits for acquiring more valuable information about hydraulic parameters and pore fluid than the conventional IP methods. The IP modeling and inversion algorithms are developed by allowing complex arithmetic in existing DC modeling and inversion algorithms. The IP modeling and inversion algorithms use a 2.5D DC finite-element algorithm and a damped least-squares method with smoothness constraints, respectively. The accuracy of the IP modeling algorithm is verified by comparing its responses of two synthetic models with two different approaches: linear filtering for a three-layer model and an integral equation method for a 3D model. Results from these methods are well matched to each other. The inversion algorithm is validated by a synthetic example which has two anomalous bodies, one is more conductive but non-polarizable than the background, and the other is polarizable but has the same resistivity as the background. From the inverted section, we can cleary identify each anomalous body with different locations. Furthermore, in order to verify its efficiency to the real filed example, we apply the inversion algorithm to another three-layer model which includes phase anomaly in the second layer.

• Journal of the Korean Geophysical Society
• /
• v.6 no.1
• /
• pp.25-29
• /
• 2003

We suggest a method that vectorially composes 3-component earthquake records in the frequency domain, to reduce the uncertainties in the analysis of the Fourier amplitude spectrum. The use of the Fourier amplitude spectrum composed from 3 components has two advantages. First, it provides a more accurate estimate of seismic moment by eliminating the vector-partition term that appears in individual component. Second, it provides more accurate estimates of seismic moment, corner frequency, high-frequency decay constant($χ$) .etc., by enhancing the shape of the Fourier amplitude spectrum. The latter is especially useful in the analysis of small earthquakes with low signal-to-noise ratios.

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

To solve the urban flood problems, it must get the enough channel conveyances and pumping capacities. It needs set up the detention system to control the flow over the channel capacity. Inspite of this detention system, the peak flow may increased by rainfall distribution and the delay of flow. This shows a design model of detention system which can consider the time problems from mutual connections of the detention storages and pumping flow using IDP(Incremental Dynamic Programming) method.