• Journal of the Korean earth science society
• /
• v.26 no.7
• /
• pp.691-697
• /
• 2005

Recently, the development of accurate gravity-meter and GPS make it possible to obtain high resolution gravity data. Though gravity data interpretation like modeling and inversion has significantly improved, gravity data processing itself has improved very little. Conventional gravity data processing removes gravity effects due to mass and height difference between base and measurement level. But, it would be a biased density model when some or whole part of anomalous bodies exist above the base level. We attempted to make a multiquadric surface of the survey area from topography with DEM (Digital Elevation Map) data. Then we constituted rectangular blocks which reflect real topography of the survey area by the multiquadric surface. Thus, we were able to carry out 3-D inversions which include information of topography. We named this technique, 3-D Gravity Terrain Inversion (3DGTI). The model test showed that the inversion model from 3DGTI made better results than conventional methods. Furthermore, the 3-dimensional model from the 3DGTI method could maintain topography and as a result, it showed more realistic geologic model. This method was also applied on real field data in Masan-Changwon area. Granitic intrusion is an important geologic characteristic in this area. This method showed more critical geological boundaries than other conventional methods. Therefore, we concluded that in the case of various rocks and rugged terrain, this new method will make better model than convention ones.

• Geophysics and Geophysical Exploration
• /
• v.20 no.1
• /
• pp.33-42
• /
• 2017

Recently, the grounded electrical-source airborne transient electromagnetic (GREATEM) system with high power source was introduced to achieve deeper investigation depth and to overcome high noise level. Although the GREATEM is a transient electromagnetic system using a long grounded wire as the transmitter, GREATEM data have been interpreted with 1D earth models because 2D or 3D modeling and inversion of vast airborne data are complicated and expensive to calculate. Generally, 1D inversion is subsequently applied to every survey point and combining 1D images together forms the stitched conductivity-depth image. However, the stitched models often result in abrupt variations in neighboring models. To overcome this problem, laterally constrained inversion (LCI) has been developed in inversion of ATEM data, which can yield layered sections with lateral smooth transitions. In this study, we analysed the GREATEM data through 1D numerical modeling for a curved grounded wire source. Furthermore, we developed a laterally constrained inversion scheme for continuous GREATEM data based on a layered earth model. All 1D data sets and models are inverted as one system, producing layered sections with lateral smooth transitions. Applying the developed LCI technique to the GREATEM data, it was confirmed that the laterally constrained inversion can provide laterally smooth model sections that reflect the layering of the survey area effectively.

• Geophysics and Geophysical Exploration
• /
• v.14 no.3
• /
• pp.220-226
• /
• 2011

In the elastic wave equations, both horizontal and vertical displacements are defined. Since we can measure both the horizontal and vertical displacements in field acquisition, these displacements compose a displacement vector. In this study, we propose a frequency-domain elastic waveform inversion technique taking advantage of the magnitudes of displacement vectors to define objective function. When we apply this displacement-vector objective function to the frequency-domain waveform inversion, the inversion process naturally incorporates the back-propagation algorithm. Through the inversion examples with the Marmousi model and the SEG/EAGE salt model, we could note that the RMS error of the solution obtained by our algorithm decreased more stably than that of the conventional method. Particularly, the density of the Marmousi model and the low-velocity sub-salt zone of the SEG/EAGE salt model were successfully recovered. Since the gradient direction obtained from the proposed objective function is numerically unstable, we need additional study to stabilize the gradient direction. In order to perform the waveform inversion using the displacementvector objective function, it is necessary to acquire multi-component data. Hence, more rigorous study should be continued for the multi-component land acquisition or OBC (Ocean Bottom Cable) multi-component survey.

• Geophysics and Geophysical Exploration
• /
• v.2 no.4
• /
• pp.167-173
• /
• 1999

Since the resolution of the 2-D hole-to-hole seismic traveltime tomography is affected by the limited ray transmission angle, various methods were used to improve the resolution. Linear traveltime interpolation(LTI) ray tracing method was chosen for forward-modeling method. Inversion results using the LTI method were compared with those using the other ray tracing methods. As an inversion algorithm, SIRT method was used. In the iterative non-linear inversion method, the cost of ray tracing is quite expensive. To reduce the cost, each raypath was stored and the inversion was performed from this information. Using the proposed method, fast convergence was achieved. Inversion results are likely to be affected by the initial velocity guess, especially when the ray transmission angle was limited. To provide a good initial guess for the inversion, generalized regression neural network(GRNN) method was used. When the transmitted raypath angle is not limited or the geological model is very complex, the inversion results are not affected by initial velocity model very much. Since the raypath angles, however, are limited in most geophysical tomographic problems, the enhancement of resolution in tomography can be achieved by providing a proper initial velocity model by another inversion algorithm such as GRNN.

• 제어로봇시스템학회:학술대회논문집
• /
• 2004.08a
• /
• pp.57-62
• /
• 2004

This paper deals with a waypoint trajectory following problem for the tilt-rotor UAV under development in Korea (TR-KUAV). In this problem, dynamic model inversion based on the linearized model and Sigma-Phi neural network with adaptive weight update are involved to realize the waypoint following algorithm for the vehicle in the helicopter flight mode (nacelle angle=0 deg). This algorithms consists of two main parts: outer-loop system as a command generator and inner-loop system as stabilizing controller. In this waypoint following problem, the position information in the inertial axis is given to the outer-loop system. From this information, Attitude Command/Attitude Hold logic in the longitudinal channel and Rate Command/Attitude Hold logic in the lateral channel are realized in the inner-loop part of the overall structure of the waypoint following algorithm. The nonlinear simulation based on the TR-KUAV is carried out to evaluate the stability and performance of the algorithm. From the numerical simulation results, the algorithm shows very good tracking performance of passing the waypoints given. Especially, it is observed that ACAH/RCAH logic in the inner-loop has the satisfactory performance due to adaptive neural network in spite of the model error coming from the linear model based inversion.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2007.06a
• /
• pp.130-134
• /
• 2007

Classical full waveform inversion for velocity estimation defines the objective function as the $l^2$ -norm of differences between the modeled and the observed wavefields. Although widely used, the results of this method have been less than satisfactory. A moderate improvement of this method is to define the objective function as the $l^2$ -norm of differences between the logarithms of the modeled and observed wavefields. In this paper we propose new objective functions of waveform inversion. They produce better results in sub-salt imaging than those of the classical and the logarithmic objective functions. One objective function defines the residual as the difference between $L^{th}$ power of the modeled wavefields and that of the observed wavefields. Another defines the residual as the difference between the integral of the $L^{th}$ power of the modeled wavefields and that of the observed wavefields. We apply these new objective functions to the synthetic SEG/EAGE salt model, and show that our new waveform inversion algorithms provide more accurate results than those of the classical and logarithmic waveform inversion methods.

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

It is difficult to obtain high-resolution images by 3D gravity inversion, because the problem is extremely underdetermined - there are too many model parameters. In order to reduce the number of model parameters we propose a 3D gravity inversion scheme utilising Euler deconvolution as a priori information. The essential point of this scheme is the reduction of the nonuniqueness of solutions by restricting the inversion space with the help of Euler deconvolution. We carry out a systematic exploration of the growing body process, but only in the restricted space within a certain radius of the Euler solutions. We have tested our method with synthetic gravity data, and also applied it to a real dataset, to delineate underground cavities in a limestone area. We found that we obtained a more reasonable subsurface density image by means of this combination between the Euler solution and the inversion process.

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

• Journal of Aerospace System Engineering
• /
• v.8 no.4
• /
• pp.24-31
• /
• 2014

The modern version of aircrafts is allowed to guarantee the superior handing qualities within the entire flight envelope by imposing the adequate stability and flying qualities on a target aircraft through the various techniques of flight control law design. Generally, the flight control law of the aircraft in service applies the various techniques of the verified control algorithm, such as dynamic inversion and eigenstructure assignment. The supersonic trainer employs the RSS(Relaxed Static Stability) concept in order to improve the aerodynamic performance in longitudinal axis and the longitudinal control laws employ the dynamic inversion with proportional-plus-integral control method. And, lateral-directional control laws employ the blended roll system of both beta-betadot feedback and simple roll rate feedback with proportional control method in order to guarantee aircraft stability. In this paper, the lateral-directional flight control law is designed by applying dynamic inversion control technique as a different method from the current supersonic trainer control technique, where the roll rate command system is designed at the lateral axis for the rapid response characteristics, and the sideslip command system is adopted at the directional axis for stability augmentation. The dynamic inversion of a simple 1st order model is applied. And this designed flight control law is confirmed to satisfy the requirement presented from the military specification. This study is expected to contribute to design the flight control law of KF-X(Korean Fighter eXperimental) which will proceed into the full-scale development in the near future.

• Geophysics and Geophysical Exploration
• /
• v.22 no.2
• /
• pp.49-55
• /
• 2019

In this study, we presented a reference-shot subset method for stable convergence of full waveform inversion using a cyclic-shot subsampling technique. Full waveform inversion needs repetitive modeling of wave propagation and thus its calculation time increases as the number of sources increases. In order to reduce the computation time, we can use a cyclic-shot subsampling method; however, it makes the cost function oscillate in the early stage of the inversion and causes a problem in applying the convergence criteria. We introduced a method in which the cost function is calculated using a fixed reference-shot subset while updating the model parameters using the cyclic-shot subsampling method. Through the examples of full waveform inversion using the Marmousi velocity model, we confirmed that the convergence of cost function becomes stable even under the cyclic-shot subsampling method if using a reference-shot subset.