• The Journal of Engineering Geology
• /
• v.14 no.1
• /
• pp.21-28
• /
• 2004

Three-dimensional (3-D) resistivity inversion including a topographic effect can be considered theoretically to be the technique of acquiring the most accurate image in the interpretation of resistivity data, because it includes characteristic image that the actual subsurface structure is 3-D. In this study, a finite-element method was used as the numerical method in modeling, and the efficiency of Jacobian calculation has been maximized with sensitivity analysis for the destination block in inversion process. Also, during the iterative inversion, the resolution of inversion can be improved with the method of selecting the optimal value of Lagrange multiplier yielding minimum RMS(root mean square) error in the parabolic equation. In this paper, we present synthetic examples to compare the difference between the case which has the toprographic effect and the other case which has not the effect in the inversion process.

• Ocean and Polar Research
• /
• v.44 no.3
• /
• pp.191-207
• /
• 2022

Using the results of CTD casts made in Spring from 2017 to 2021, in this study we investigated the water mass distribution and occurrence of temperature inversion in the western seas of Jeju Island in spring. The distribution of water masses was characterized by cold and fresh water in the northwest and warm and saline water in the southeast, forming a strong thermohaline front running in the southwest-to-northeast direction. Strong temperature inversion mainly occurred in the frontal boundary when the cold water intrudes beneath the warm water at depths of 30-50 m. Analysis of the mixing ratio demonstrated that Jeju Warm Water is dominantly distributed in the western seas of Jeju Island, but its ratio can be modified depending on the southward extension of Yellow Sea Cold Water (YSCW). Results of in situ measurement showed that in 2020, the YSCW largely expanded to the western seas of Jeju Island, occupying approximately 40 % of the mixing ratio. Due to the expansion of YSCW, a strong thermohaline front was formed in the study area, thereby causing thick and strong temperature inversion. On the other hand, in 2018 the mixing ratio of YSCW was minimum (~18%) during the study period of 2017-2021, and thus a relatively weak frontal boundary was formed, without the occurrence of temperature inversion. The observational results also suggest that the interannual changes of water mass distribution and the associated temperature inversion in the western seas of Jeju Island are closely related with wind-driven Yellow Sea circulation in spring, which is the summer monsoon transition period.

• Geophysics and Geophysical Exploration
• /
• v.25 no.4
• /
• pp.252-265
• /
• 2022

Subsurface velocity is critical for the accurate resolution geological structures. The estimation of acoustic impedance is also critical, as it provides key information regarding the reservoir properties. Therefore, researchers have developed various inversion approaches for the estimation of reservoir properties. The Markov chain Monte Carlo method, which is a stochastic method, has advantages over the deterministic method, as the stochastic method enables us to attenuate the local minima problem and quantify the uncertainty of inversion results. Therefore, the Markov chain Monte Carlo inversion method has been applied to various kinds of geophysical inversion problems. However, studies on the Markov chain Monte Carlo inversion are still very few compared with deterministic approaches. In this study, we reviewed various types of reversible jump Markov chain Monte Carlo applications and explained the key concept of each application. Furthermore, we discussed future applications of the stochastic method.

• Proceedings of the Korean Ceranic Society Conference
• /
• 2004.10a
• /
• pp.81.1-81.1
• /
• 2004

• Geophysics and Geophysical Exploration
• /
• v.10 no.1
• /
• pp.1-18
• /
• 2007

Non-linear elastic wavefield inversion is a powerful method for estimating elastic parameters for physical constraints that determine subsurface rock and properties. Here, I introduce six elastic-wave velocity models by reconstructing elastic-wave velocity variations from real data and a 2D elastic-wave velocity model. Reflection seismic data information is often decoupled into short and long wavelength components. The local search method has difficulty in estimating the longer wavelength velocity if the starting model is far from the true model, and source frequencies are then changed from lower to higher bands (as in the 'frequency-cascade scheme') to estimate model elastic parameters. Elastic parameters are inverted at each inversion step ('simultaneous mode') with a starting model of linear P- and S-wave velocity trends with depth. Elastic parameters are also derived by inversion in three other modes - using a P- and S-wave velocity basis $('V_P\;V_S\;mode')$; P-impedance and Poisson's ratio basis $('I_P\;Poisson\;mode')$; and P- and S-impedance $('I_P\;I_S\;mode')$. Density values are updated at each elastic inversion step under three assumptions in each mode. By evaluating the accuracy of the inversion for each parameter set for elastic models, it can be concluded that there is no specific difference between the inversion results for the $V_P\;V_S$ mode and the $I_P$ Poisson mode. The same conclusion is expected for the $I_P\;I_S$ mode, too. This gives us a sound basis for full wavelength elastic wavefield inversion.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.68-77
• /
• 2008

We have developed an inversion algorithm for loop-loop electromagnetic (EM) data, based on the localised non-linear or extended Born approximation to the solution of the 2.5D integral equation describing an EM scattering problem. Source and receiver configuration may be horizontal co-planar (HCP) or vertical co-planar (VCP). Both multi-frequency and multi-separation data can be incorporated. Our inversion code runs on a PC platform without heavy computational load. For the sake of stable and high-resolution performance of the inversion, we implemented an algorithm determining an optimum spatially varying Lagrangian multiplier as a function of sensitivity distribution, through parameter resolution matrix and Backus-Gilbert spread function analysis. Considering that the different source-receiver orientation characteristics cause inconsistent sensitivities to the resistivity structure in simultaneous inversion of HCP and VCP data, which affects the stability and resolution of the inversion result, we adapted a weighting scheme based on the variances of misfits between the measured and calculated datasets. The accuracy of the modelling code that we have developed has been proven over the frequency, conductivity, and geometric ranges typically used in a loop-loop EM system through comparison with 2.5D finite-element modelling results. We first applied the inversion to synthetic data, from a model with resistive as well as conductive inhomogeneities embedded in a homogeneous half-space, to validate its performance. Applying the inversion to field data and comparing the result with that of dc resistivity data, we conclude that the newly developed algorithm provides a reasonable image of the subsurface.

• Geophysics and Geophysical Exploration
• /
• v.16 no.1
• /
• pp.27-35
• /
• 2013

In this study, we developed the trial product on 3D magnetic inversion tentatively named 'KMag3D'. Also, we briefly introduced its own function and graphic user interface on which especially focused through the development in the form of user manual. KMag3D is consisted of two fundamental frame for the 3D magnetic inversion. First, algebraic reconstruction technique was selected as a 3D inversion algorithm instead of least square method conventionally used in various magnetic inversion. By comparison, it was turned out that algebraic reconstruction algorithm was more effective and economic than that of least squares in aspect of both computation time and memory. Second, for the effective determination of the 3D initial and a-priori information model required in the execution of our algorithm, we proposed the practical technique based on the assemblage of 2D forward modeling and inversion results for individual user-selected 2D profiles. And in succession, initial and a-priori information model were constructed by appropriate interpolation along the strke direction. From this, we concluded that our technique is both suitable and very practical for the application of 3D magentic inversion problem.

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

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

• Journal of Environmental Science International
• /
• v.3 no.3
• /
• pp.217-228
• /
• 1994

The climatological study of temperature inversion was conducted at the Tsukuba of Japan during the cold half year, the periods of which were both from October 1981 through March 1982 and from October 1982 through March 1983. Meteorological phenomena were observed to appear both from l0m to 200m above the ground and from the surface to 300hPa. The data collected from those phenomena were analyzed in terms of statistics and investigated from a synoptic point of view. The results are as follow. The dispersion of the surface inversion happens right after sunrise through the whole period. The higher the upper layer is, the later that happens. Up to a height of 200m, the wind speed at a height of 25m has the greatest effect on temperature, but on the other hand the wind speed at a height of l0m has the greatest effect on inversion intensity.