• Geophysics and Geophysical Exploration
• /
• v.16 no.3
• /
• pp.131-138
• /
• 2013

We developed a method for inverting magnetic data to image 2D susceptibility models. The major difficulty in the inversion of the potential data is the nonuniqueness. Furthermore, generally the number of inversion blocks are greater than the number of the magnetic data available, and thus the magnetic inversion leads to under-determined problem, which aggravates the nonuniqueness. When the magnetic data were inverted by the general least-squares method, the anomalous susceptibility would be concentrated near the surface in the inverted section. To overcome this nonuniqueness problem, we propose a new resolution model constraint that is calculated from the parameter resolution. The model constraint imposes large penalty on the model parameter with good resolution, on the other hand small penalty on the model parameter with poor resolution. Thus, the deep-seated model parameter, generally having poor resolution, can be effectively resolved. The developed inversion algorithm is applied to the inversion of the synthetic data for typical models of magnetic anomalies and is tested on real airborne data obtained at the Okcheon belt of Korea.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2008.10a
• /
• pp.51-56
• /
• 2008

The waveform inversion for isotropic media has ever been studied since the 1980s, but there has been few studies for anisotropic media. We present a seismic waveform inversion algorithm for 2-D heterogeneous transversely isotropic structures. A cell-based finite difference algorithm for anisotropic media in time domain is adopted. The steepest descent during the non-linear iterative inversion approach is obtained by backpropagating residual errors using a reverse time migration technique. For scaling the gradient of a misfit function, we use the pseudo Hessian matrix which is assumed to neglect the zero-lag auto-correlation terms of impulse responses in the approximate Hessian matrix of the Gauss-Newton method. We demonstrate the use of these waveform inversion algorithm by applying them to a two layer model and the anisotropic Marmousi model data. With numerical examples, we show that it's difficult to converge to the true model when we assumed that anisotropic media are isotropic. Therefore, it is expected that our waveform inversion algorithm for anisotropic media is adequate to interpret real seismic exploration data.

• Geophysics and Geophysical Exploration
• /
• v.18 no.1
• /
• pp.9-13
• /
• 2015

All the resistivity data contain various kinds of noise. The major sources of noise in DC resistivity measurement are high contact resistance, measurement errors, and sporadic background noise. Thus, it is required to measure data noise to accurately interpret resistivity data. Reciprocal measurements can provide a measure of data precision and noise. In this study, we proposed a data-weighting method from reciprocity measurement. Furthermore, a data-weighting method using both the reciprocity error and data-misfit in the inversion process was studied. Applying the data-weighting method to the inversion of 3D resistivity data, it was confirmed that local anomalies are slightly suppressed in the final inversion results.

• Economic and Environmental Geology
• /
• v.29 no.3
• /
• pp.345-355
• /
• 1996

The effectiveness of various electrode configurations in horizontal mappings and 1-D inversions of vertical sounding data for delineating 2-D structures was studied. Apparent resistivity values of three point, dipole-dipole, Wenner, and Schlumberger mappings were simulated for such structures as vertical dyke, tabular prism, buried vertical fault, ramp and complex structure by finite difference method (FDM) and they were compared with each other. Also 2-D cross sections for three structures obtained by interpolation of 1-D inverted sounding data in terms of three layers were compared for Schlumberger and Wenner arrays. On these cross sections, horizontal and vertical resistivity interfaces of the 2-D structures are revealed relatively clearly. Apparent resistivity curves of Schlumberger mapping show vertical resistivity discontinuities very well. On the whole, Schlumberger array is superior to the other arrays in electric sounding as well as mapping. This study clearly indicates that interpretations of 2-D structures based on 1-D inversion are possible.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.377-382
• /
• 2007

Several prospecting methods have been used to detect deep seated small tunnel in Korea. Tunnel interpretation of seismic method has been performed mainly by wave traveltime inversion method. But it often gives inacurate solution for the exact tunnel position because of the short distance between two measuring boreholes and picking errors of first arrivals. In this study, "error tomogram" was proposed to detect tunnel position and applied to theoretical and field dat using multi-source amplitude data.

• Proceedings of the Korean Geotechical Society Conference
• /
• 1996.10a
• /
• pp.203-208
• /
• 1996

• Geophysics and Geophysical Exploration
• /
• v.16 no.3
• /
• pp.119-130
• /
• 2013

We developed a method for inverting magnetic data to recover the 3D susceptibility models. The major difficulty in the inversion of the potential data is the non-uniqueness and the vast computing time. The insufficient number of data compared with that of inversion blocks intensifies the non-uniqueness problem. Furthermore, there is poor depth resolution inherent in magnetic data. To overcome this non-uniqueness problem, we propose a resolution model constraint that imposes large penalty on the model parameter with good resolution; on the other hand, small penalty on the model parameter with poor resolution. Using this model constraint, the model parameter with a poor resolution can be effectively resolved. Moreover, the wavelet transform and parallel solving were introduced to save the computing time. Through the wavelet transform, a large system matrix was transformed to a sparse matrix and solved by a parallel linear equation solver. This procedure is able to enormously save the computing time for the 3D inversion of magnetic data. The developed inversion algorithm is applied to the inversion of the synthetic data for typical models of magnetic anomalies and real airborne data obtained at the Geumsan area of Korea.

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

In the geophysical monitoring to understand the change of subsurface material properties with time, the time-invariant static subsurface model is commonly adopted to reconstruct a time-lapse image. This assumption of static model, however, can be invalid particularly when fluid migrates very quickly in highly permeable medium in the brine injection experiment. In such case, the resultant subsurface images may be severely distorted. In order to alleviate this problem, we develop a new least-squares inversion algorithm under the assumption that the subsurface model will change continuously in time. Instead of sampling a time-space model into numerous space models with a regular time interval, a few reference models in space domain at different times pre-selected are used to describe the subsurface structure continuously changing in time; the material property at a certain space coordinate are assumed to change linearly in time. Consequently, finding a space-time model can be simplified into obtaining several reference space models. In order to stabilize iterative inversion and to calculate meaningful subsurface images varying with time, the regularization along time axis is introduced assuming that the subsurface model will not change significantly during the data acquisition. The performance of the proposed algorithm is demonstrated by the numerical experiments using the synthetic data of crosshole dc resistivity tomography.

• Journal of the Korean Geophysical Society
• /
• v.2 no.2
• /
• pp.135-142
• /
• 1999

Interpretations of 3-component magnetic logging data obtained for a reinforced bar as a model of the line of the magnetic dipoles are conducted using a least squared inversion technique. The length of the bar is 1.12 m, sampling interval is 0.05 m, the distance between the bar and the borehole is 0.3 m, and the top of the bar is fixed at 0 m of depth. The bar is set to be approximately vertical. Magnetic anomalies smoothed with FFT are used as input data for the inversion. For the interpretation of magnetic logging data the depth to the top, the length, the magnetic moment per unit length, the direction of the magnetization (declination and inclination), and the bearing and plunge of the line of magnetic dipoles are left as unknown parameters. The comparison of the results obtained from the individual inversion of the horizontal component or the vertical component of the magnetic anomalies, and those from the simultaneous inversion of horizontal and vertical component of the magnetic anomalies shows that there exist some disagreements between each inversion result. The depth to the bottom of the bar, which is actually 1.12 m, is estimated as 1.18 m, and the inclination of the magnetization is estimated as -76°by simultaneous inversion. The negative value of the inclination indicates that the strength of the remnant magnetization is much greater than that of the induced magnetization, so that the direction of the resultant magnetization points to the top of the bar.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.93-101
• /
• 2006

In this study, we propose a joint inversion method, using genetic algorithms, to determine the shear-wave velocity structure of deep sedimentary layers from receiver functions and surface-wave phase velocity. Numerical experiments with synthetic data indicate that the proposed method can avoid the trade-off between shear-wave velocity and thickness that arises when inverting the receiver function only, and the uncertainty in deep structure from surface-wave phase velocity inversion alone. We apply the method to receiver functions obtained from earthquake records with epicentral distances of about 100 km, and Rayleigh-wave phase velocities obtained from a microtremor array survey in the Kanto Plain, Japan. The estimated subsurface structure is in good agreement with the previous results of seismic refraction surveys and deep borehole data.