• 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 Industrial Technology
• /
• v.23 no.A
• /
• pp.63-68
• /
• 2003

Non-destructive method by tomography for safety diagnosis of civil engineering and building structures is tried. There are traveltime tomography that uses traveltime and fullwave tomography that uses the initial shock and seismic amplitude. But these methods have difficulty and weak points in accuracy and selection of initial value. In this study, corrected inversion method which is able to solve the two difficulty and this method is applied to theoretical pier model to calculate tomography.

• Proceedings of the KSEEG Conference
• /
• 1999.04a
• /
• pp.63-65
• /
• 1999

• Journal of the Korean Geophysical Society
• /
• v.8 no.2
• /
• pp.67-74
• /
• 2005

Two-dimensional S-wave velocity sections from SH-wave refraction tomography and surface wave dispersions were obtained by inverting traveltimes of first arrivals and surface wave dispersions, respectively. For the purpose of comparison, a P-wave velocity tomogram was also obtained from a P-wave refraction profiling. P and Rayleigh waves generated by vertical blows on a plate with a sledgehammer were received by 100- and 4.5-Hz geophones, respectively. SH-waves generated by horizontal blows on both sides of a 50 kg timber were received by 8 Hz horizontal geophones. The shear-wave signals were enhanced subtracting data of left-side blows from ones of the right-side blows. Shear-wave velocities from tomography inversion of first-arrival times were compared with ones from inverting dispersion curves of Rayleigh waves. Although the two velocity sections look similar to each other in general, the one from the surface waves tends to have lower velocities. First arrival picking of SH waves is troublesome since P and PS-converted waves arrive earlier than SH waves. Application of the surface wave method, on the other hand, is limited where lateral variation of subsurface tructures is not mild.

• Journal of the Korean earth science society
• /
• v.22 no.3
• /
• pp.179-185
• /
• 2001

Borehole tests are commonly used as a tool to obtain the physical properties of soils and rocks. The results of borehole tests are, however, discontinuous. Interpolation methods are applied to interpret the data gap between the borehole test points. The interpolation is valid only if the horizontal variations of the ground between the test points are small enough to ignore. A surface wave inversion method was used to study the S wave velocity of the very soft soil to provide the continuous 2 dimensional S wave velocity structure. The resolution of the S wave velocity structure was used to interpret the inversion results.

• Proceedings of the Korea Electromagnetic Engineering Society Conference
• /
• 2005.11a
• /
• pp.169-174
• /
• 2005

In this paper, we proposed a new algorithm of the inverse scattering for the reconstruction of unknown dielectric scatterers using the finite-difference time-domain method and the design sensitivity analysis. We introduced the design sensitivity analysis based on the gradient for the fast convergence of the reconstruction. By introducing the adjoint variable method for the efficient calculation, we derived the adjoint variable equation. As an optimal algorithm we used the steepest descent method and reconstructed the dielectric targets using the iterative estimation. To verify our algorithm we will show the numerical examples for the two-dimensional $TM^2$ cases.

• Economic and Environmental Geology
• /
• v.28 no.5
• /
• pp.481-485
• /
• 1995

This paper presents a new approaching method to determine the velocity and geometry of shallow subsurface from seismic refraction events. After picking the first breaks from seismic refraction data, we assume that field refraction seismogram can be replaced by the unit delta function having time shift of first break. Time curves are generated by shooting ray tracing. The partial derivatives seismogram for a damped least squares method is computed analytically at each step of the forward ray tracing. The technique is successfully tested on synthetic and real data. It has the advantage of real full waveform inversion, which is robust at low frequency band even if the initial guess is far from the true model.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.259-267
• /
• 2007

An efficient algorithm for inverting static-shifted magnetotelluric (MT) data has been proposed to produce a three-dimensional (3D) resistivity model. In the Gauss-Newton approach, computational costs associated with construction of a full sensitivity matrix usually make 3D MT inversion impractical. This computational difficulty may be overcome by using approximate sensitivities. We use four kinds of sensitivities in particular orders in the inversion process. These sensitivities are computed 1) analytically for an initial, homogeneous earth, 2) exactly for a current model, 3) approximately by the Broyden method, and 4) approximately using the previous adjoint fields. Inversion experiments with static-shifted synthetic and field MT data indicate that inversion results are highly dependent on characteristics of data and thus applying various combinations of sensitivities is helpful in obtaining a good image of the subsurface structure with reasonable computation time.

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

Two-dimensional velocity tomograms of P- and S-waves were obtained by inverting traveltimes of first arrivals. The two sections of shear-wave velocity show similar features as a whole, with smaller values on the section from surface wave dispersions. Difficulties in picking SH-wave phases due to noise and later arrivals than P waves and PS converted waves are experienced. In addition, a flat layer model based on the surface wave inversion prohibits applications of the method where sgear wave velocities vary strongly in the lateral direction.

• 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.