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

Two-dimensional (2-D) inversion of magnetotelluric (MT) data for two survey lines having south-north direction from Jeiu Island has been carried out. The 2-D models show a thick layer having around 10 ohm-m in the depth of a few hundred meters throughout the survey area, which can be considered as the unconsolidated sedimentary layer. And they also show a conductive anomaly at the central part of each survey lines. But unfortunately by now, we do not have any further information about the anomaly. Comparison of the 2-D inversion model using MT band only and that using both AMT and MT bands said that it is helpful for us to include AMT band as well as MT band in the inversion to interpret not only the shallow part but also the deep structures.

• Geophysics and Geophysical Exploration
• /
• v.22 no.4
• /
• pp.186-194
• /
• 2019

Deterministic optimization, commonly used to find the geophysical inverse solutions, have its limitation that it cannot find the proper solution since it might converge into the local minimum. One of the solutions to this problem is to use global optimization based on a stochastic approach, among which a large number of particle swarm optimization (PSO) applications have been introduced. In this paper, I developed a geophysical inversion algorithm applying PSO method for the layered-earth resistivity inversion of the small-loop electromagnetic (EM) survey data and carried out numerical inversion experiments on synthetic datasets. From the results, it is confirmed that the PSO inversion algorithm could increase the inversion success rate even when attempting the inversion of small-loop EM survey data from which it might be difficult to find a best solution by applying the Gauss-Newton inversion algorithm.

• Journal of the Korean Geotechnical Society
• /
• v.21 no.4
• /
• pp.155-165
• /
• 2005

Love wave and Rayleigh wave are the major elastic waves belonging to the category of the surface wave. Those waves are used to determine the ground stiffness profile using their dispersion characteristics. The fact that Love wave is not contaminated by P-wave makes Love wave superior to Rayleigh wave and other body waves. Therefore, the information that Love wave carries is more distinct and clearer than that of others. Based on theoretical research, the joint inversion analysis that uses the dispersion information of both Love and Rayleigh wave was proposed. Numerical analysis, theoretical model test, and field test were performed to verify the joint inversion analysis. Results from 2D, 3D finite element analysis were compared with those from the transfer matrix method in the numerical analysis. On the other hand, the difference of results from each inversion analysis was investigated in the theoretical model analysis. Finally, practical applicability of the joint inversion analysis was verified by performing field test. As a result, it is confirmed that considering dispersion information of each wave simultaneously prevents excessive divergence and improves accuracy.

• Journal of the Korean Institute of Telematics and Electronics D
• /
• v.36D no.4
• /
• pp.34-41
• /
• 1999

Unlike the traditional radar system, Synthetic Aperture Radar(SAR) system is capable of imaging a target scene to ceertain degree of cross-range resolution. And this resolution is mainly depends on the size of aperture synthesized. Thus, a good system model and inversion scheme should be developed to actually give effect of synthesizing aperture size, which in turn gives better cross range resolution of reconstructed target scene. Among several inversion schemes for SAR imaging, we used an inversion scheme called wavefront reconstruction which has no approximation in wave propagation analysis, and tried to verify whether the collected data with synthesized aperture actually give the same support as that with physical aperture in the same size. To do this, we performed a signal subspace comparison of two imaging models with physical and synthesized arrays, respectively. Theoretical comparisons and numerical analysis using Gram-Schmidt procedures have been performed. The results showed that the synthesized array data fully span the physical array data with the same system geometry. This result strongly supports the previously proposed inversion scheme valuable in high resolution radar imaging.

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.291-298
• /
• 2002

We have developed a two-dimensional (2-D) magnetotelluric (MT) inversion algorithm, which can include topographic effects in inversion. We use the finite element method (FEM) to incorporate topography into forward calculation. Topography is implemented simply by moving nodes of rectangular elements in z-direction according to the elevation of air-earth interface. In the inversion process, we adopt a spatially variable Lagrangian multiplier algorithm in the smoothness-constrained least-squares inversion. The inversion algorithm developed in this study reconstructs subsurface resistivity structure quite well when topography variation exists. Also, it turns out to be effective in both resolution and stability from a model study and field data application.

• Geophysics and Geophysical Exploration
• /
• v.24 no.4
• /
• pp.180-193
• /
• 2021

The complex resistivity method is an exploration technique that can obtain various characteristic information of underground media by measuring resistivity and phase in the frequency domain, and its utilization has recently increased. In this paper, a three-dimensional inversion algorithm for the CR data was developed to increase the utilization of this method. The Poisson equation, which can be applied when the electromagnetic coupling effect is ignored, was applied to the modeling, and the inversion algorithm was developed by modifying the existing algorithm by adopting comlex variables. In order to increase the stability of the inversion, a technique was introduced to automatically adjust the Lagrangian multiplier according to the ratio of the error vector and the model update vector. Furthermore, to compensate for the loss of data due to noisy phase data, a two-step inversion method that conducts inversion iterations using only resistivity data in the beginning and both of resistivity and phase data in the second half was developed. As a result of the experiment for the synthetic data, stable inversion results were obtained, and the validity to real data was also confirmed by applying the developed 3D inversion algorithm to the analysis of field data acquired near a hydrothermal mine.

• Geophysics and Geophysical Exploration
• /
• v.16 no.4
• /
• pp.217-224
• /
• 2013

We developed a time-lapse inversion using new cross-model constraints based on change ratio and resolution of model parameters. The cross-model constraint based on change ratio imposes the same penalty on the model parameters with equal change ratio. This constraint can emphasize the model parameters with significant change regardless of their increase or decrease. The resolution cross-model constraint imposes a small penalty on the model parameters with poor resolution, but a large penalty on the model parameters with good resolution. Thus, the model parameter with poor resolution can be effectively identified in the inversion result if they are significantly changed with time. Through the numerical tests for 3D resistivity monitoring data sets, the performance of these two cross-model constraints was confirmed. Finally, for the safety estimation of a sea dyke, we applied the developed time-lapse inversion to the 3D resistivity monitoring data that were acquired at a sea dike located in western coastal area of Korea. The result of time-lapse inversion suggested that there were no significant changes at the sea dike during the monitoring period.

• Atmosphere
• /
• v.31 no.4
• /
• pp.405-420
• /
• 2021

Forest fire happens every year at Yeongdong, Gangwon-do, due to the strong local wind during the spring time and it causes a huge damage. This wind is named "Yangganjipung" or "Yanggangjipung" that blows along Yeongdong. However, the occurrence conditions of the wind have been still unclear. To identify the occurrence mechanism of local strong wind through three-dimensional observation data, Gangwon Regional Meteorological Administration performed Joint Gangwon-Yeongdong 3D Observation Project in 2020. The special observation was carried out for 6 times from March to April. The observation data was analyzed by focusing on the structure of synoptic pressure distribution and inversion layer. The result showed that the strength of wind is different depending on the latitude of low pressure, intensity of inversion layer, and changes on height in the south-high and north-low pressure distribution. As the interval of the upper and lower parts of the inversion layer was narrow, the strength of the wind became stronger, which is one of the observational characteristics of the springtime wind pattern at Yeongdong, Gangwon-do. In future, the clear mechanism of the local wind in the Yeongdong during the spring time is expected to be verified based on the accumulative observation data and close analysis.

• Geophysics and Geophysical Exploration
• /
• v.24 no.4
• /
• pp.139-148
• /
• 2021

Negative apparent chargeability data in time-domain induced polarization (IP) survey are very often considered as noise and rejected before the inversion stage. Negative IP data can, however, occur naturally as a consequence of the distribution of chargeable zones in the ground. In some cases, the negative apparent chargeability values may account for most of the data measured. Negative IP data are caused by the geometry of chargeable zones and electrode positions. Negative apparent chargeability data appear most frequently when a dipole-dipole array is used. In this study, the effect of negative apparent chargeability data on inversion results is analyzed through the numerical 2D time-domain IP modeling and nonlinear inversion. The results demonstrate clearly that negative apparent chargeability data have to be included in the inversion as they contain important information on the distribution of subsurface chargeability.

• Geophysics and Geophysical Exploration
• /
• v.19 no.1
• /
• pp.1-10
• /
• 2016

The hypocenter distribution of microseismic events generated by hydraulic fracturing for shale gas development provides essential information for understanding characteristics of fracture network. In this study, we evaluate how inaccurate velocity models influence the inversion results of two widely used location programs, hypoellipse and hypoDD, which are developed based on an iterative linear inversion. We assume that 98 stations are densely located inside the circle with a radius of 4 km and 5 artificial hypocenter sets (S0 ~ S4) are located from the center of the network to the south with 1 km interval. Each hypocenter set contains 25 events placed on the plane. To quantify accuracies of the inversion results, we defined 6 parameters: difference between average hypocenters of assumed and inverted locations, $d_1$; ratio of assumed and inverted areas estimated by hypocenters, r; difference between dip of the reference plane and the best fitting plane for determined hypocenters, ${\theta}$; difference between strike of the reference plane and the best fitting plane for determined hypocenters, ${\phi}$; root-mean-square distance between hypocenters and the best fitting plane, $d_2$; root-mean-square error in horizontal direction on the best fitting plane, $d_3$. Synthetic travel times are calculated for the reference model having 1D layered structure and the inaccurate velocity model for the inversion is constructed by using normal distribution with standard deviations of 0.1, 0.2, and 0.3 km/s, respectively, with respect to the reference model. The parameters $d_1$, r, ${\theta}$, and $d_2$ show positive correlation with the level of velocity perturbations, but the others are not sensitive to the perturbations except S4, which is located at the outer boundary of the network. In cases of S0, S1, S2, and S3, hypoellipse and hypoDD provide similar results for $d_1$. However, for other parameters, hypoDD shows much better results and errors of locations can be reduced by about several meters regardless of the level of perturbations. In light of the purpose to understand the characteristics of hydraulic fracturing, $1{\sigma}$ error of velocity structure should be under 0.2 km/s in hypoellipse and 0.3 km/s in hypoDD.