• 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.10 no.4
• /
• pp.314-321
• /
• 2007

Recently, multi-dimensional joint inversion of geophysical data based on fundamentally different physical properties is being actively studied. Joint inversion can provide a way to obtaining much more accurate image of the subsurface structure. Through the joint inversion, furthermore, it is possible to directly estimate non-geophysical material properties from geophysical measurements. In this study, we developed a new algorithm for jointly inverting dc resistivity and seismic traveltime data based on the multiple constraints: (1) structural similarity based on cross-gradient, (2) correlation between two different material properties, and (3) a priori information on the material property distribution. Through the numerical experiments of surface dc resistivity and seismic refraction surveys, the performance of the proposed algorithm was demonstrated and the effects of different regularizations were analyzed. In particular, we showed that the hidden layer problem in the seismic refraction method due to an inter-bedded low velocity layer can be solved by the joint inversion when appropriate constraints are applied.

• Geophysics and Geophysical Exploration
• /
• v.24 no.4
• /
• pp.149-157
• /
• 2021

Given that induced polarization (IP) and direct current (DC) resistivity surveys are similar in terms of data acquisition, most DC resistivity systems are equipped with a time-domain IP data acquisition function. In addition, the time-domain IP data include the DC resistivity values. As such, IP and DC resistivity data are intimately linked, and the inversion of IP data is a two-step process based on DC resistivity inversions. Nevertheless, IP surveys are rarely applied, in contrast to DC resistivity surveys, as proper inversion software is unavailable. In this study, through numerical modeling and inversion experiments, we analyze the problems with the conventional inverse mapping technique used to invert time-domain IP data. Furthermore, we propose a modified inverse mapping technique that can effectively suppress inversion artifacts. The performance of the technique is confirmed through inversions applied to synthetic IP data.

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

• 제어로봇시스템학회:학술대회논문집
• /
• 1989.10a
• /
• pp.595-599
• /
• 1989

A new simple method for controlling compliant motions of a flexible robot arm is presented. The method aims at controlling translational tip motion, force and moment by directly computing the base motion or torque. A numerical inversion of Laplace transform is used to obtain the results in the time domain. The results show the effectiveness of the method for the hybrid translational position/force control of a flexible robot arm.

• The Journal of the Acoustical Society of Korea
• /
• v.39 no.3
• /
• pp.151-162
• /
• 2020

This paper presents an acoustic inversion method for estimating the bubble size distribution. The estimation error of the attenuation coefficient represented by a Fredholm integral equation of the first kind is defined as an objective function, and an optimal solution is found by applying the Levenberg-Marquardt (LM) method. In order to validate the effectiveness of the inversion method, numerical simulations using two types of bubble distribution are performed. In addition, a series of experiments are carried out in a water tank (1.0 m × 0.54 m × 0.6 m), using bubbles generated by three different generators. Images of the distributed bubbles are obtained by a high-speed camera, and the insertion losses of the bubble layer are measured using a source and a hydrophone. The image is post-processed to glance a distribution characteristics of each bubble generator. Finally, the size distribution of bubbles is estimated by applying the inversion method to the measured insertion loss. From the inversion results, it was observed that the number of bubbles increases exponentially as the bubble size decreases, and then increases again after the local peak at 70 ㎛ - 120 ㎛.

• Proceedings of the KIEE Conference
• /
• 1999.11a
• /
• pp.43-51
• /
• 1999

The power flow calculations(PFc) are the most important and powerful tools in power systems engineering. The conventional power flow problem is solved generally with numerical methods such as Newton-Raphson(NR). The conventional numerical method generally have some convergency problem, which is sensitive to initial value, and numerical stability problem concerned with jacobian matrix inversion. This paper presents a new PFc algorithm based on the improved genetic algorithm (IGA) which can overcome the disadvantages mentioned above. The parameters of GA, with dynamical hierarchy of the coding system, are improved to make GA a practical algorithm in the problem of real system. Some case studies with test bus system also present to show the performance of proposed algorithm. The results of proposed algorithm are compared with the results of PFc obtained using a conventional NR method.

• Structural Engineering and Mechanics
• /
• v.54 no.4
• /
• pp.809-827
• /
• 2015

In this paper, seismic energy response of inelastic steel structures under earthquake excitations is investigated. For this purpose, a numerical procedure based on nonlinear dynamic analysis is developed by considering material, geometric and connection nonlinearities. Material nonlinearity is modeled by the inversion of Ramberg-Osgood equation. Nonlinearity caused by the interaction between the axial force and bending moment is also defined considering stability functions, while the geometric nonlinearity caused by axial forces is described using geometric stiffness matrix. Cyclic behaviour of steel connections is taken into account by employing independent hardening model. Dynamic equation of motion is solved by Newmark's constant acceleration method in the time history domain. Energy response analysis of space frames is performed by using this proposed numerical method. Finally, for the first time, the distribution of the different energy types versus time at the duration of the earthquake ground motion is obtained where in addition error analysis for the numerical solutions is carried out and plotted depending on the relative error calculated as a function of energy balance versus time.

• Steel and Composite Structures
• /
• v.25 no.2
• /
• pp.177-186
• /
• 2017

In this work, the model of magneto-thermoelasticity based on memory-dependent derivative (MDD) is applied to a one-dimensional thermal shock problem for a functionally graded half-space whose surface is assumed to be traction free and subjected to an arbitrary thermal loading. The $Lam{\acute{e}}^{\prime}s$ modulii are taken as functions of the vertical distance from the surface of thermoelastic perfect conducting medium in the presence of a uniform magnetic field. Laplace transform and the perturbation techniques are used to derive the solution in the Laplace transform domain. A numerical method is employed for the inversion of the Laplace transforms. The effects of the time-delay on the temperature, stress and displacement distribution for different linear forms of Kernel functions are discussed. Numerical results are represented graphically and discussed.

• Geophysics and Geophysical Exploration
• /
• v.19 no.3
• /
• pp.136-144
• /
• 2016

Since seismic inversion is based on the wave equation, it is important to calculate the solution of wave equation exactly. In particular, full waveform inversion would produce reliable results only when the forward modeling is accurately performed because it uses full waveform. When we use finite-difference or finite-element method to solve the wave equation, the convergence of numerical scheme should be guaranteed. Although the general proof of convergence is provided theoretically, the consistency and stability of numerical schemes should be verified for practical applications. The implementation of source function is the most crucial factor for the consistency of modeling schemes. While we have to use the sinc function normalized by grid spacing to correctly describe the Dirac delta function in the finite-difference method, we can simply use the value of basis function, regardless of grid spacing, to implement the Dirac delta function in the finite-element method. If we use frequency-domain wave equation, we need to use a conservative criterion to determine both sampling interval and maximum frequency for the source wavelet generation. In addition, the source wavelet should be attenuated before applying it for modeling in order to make it obey damped wave equation in case of using complex angular frequency. With these conditions satisfied, we can develop reliable inversion algorithms.